各位大侠帮忙看下这代码有什么错误 [复制链接]

title 5]kv1nQ  
Oedometric test phi=10 psi=10 W%MS,zkAE  
round 0.01 2PNe~9)*#  
block 0 0 0 1 1 1 1 0 '!MKZKer  
gen quad 2 A^|~>9  
prop mat 1 dens 1 s|[qq7  
;change cons 3 +sq_fd ;'D  
;prop mat 1 bu 200 sh 200 co 1 fric 10dil 10 <&((vrfa  
zone m m =<TJ[,h et  
zone bulk 200 sh 200 co 1 fric 10 eT2Tg5Etc  
zone dil 10 VTX6_&Hc1g  
;zone dil 0 #op0|:/

N  
step 0 bq8h?Q  
def d_sigy .WBp!*4  
c_k = bulk_mod bx-:aC)]2  
c_g = shear_mod v@fy*T\3  
e1 = c_k+4.0*c_g/3.0 _$8:\[J  
e2 = c_k-2.0*c_g/3.0 ExFz@6@  
sf = friction * degrad JPZH%#E(  
nf = sin(sf) 4T~wnTH0Xg  
nf = (1.0+nf)/(1.0-nf) # xX  
sp = dilation*degrad SoFl]^l  
np = sin(sp) puh-\Q/P  
np = (1.0+np)/(1.0-np) [CAFh:o  
r1 = (e1-e2*nf)/((e1+e2)*nf*np-2.0*e2*(nf+np)+2.0*e1) !@arPN$  
v_dt = -1.0e-5 M D&7k,!  
velv =v_dt/tdel OeZ"WO  
dsigy = v_dt*(e1+2.0*r1*(e2*np-e1)) EAC
I>  
step1 = -2.0*cohesion*sqrt(nf)/((e1-e2*nf)*v_dt) HqyAo]{GN  
end 5CfD/}{:#I  
def esigy JZ> (h  
while_stepping U{@2kg-  
if step < step1 then SC3_S.  
   a_sy = a_sy+e1*v_dt s._,IW;   
else UQP>yuSx  
   a_sy = a_sy+dsigy g">^#^hBE  
end_if fL-$wK<p<  
   n_sy = z_syy(z_pnt) yKz%-6cpSl  
end [{>1wJ Pdj  
set bulk_mod 200 shear_mod 200 ]T.+(\I  
set cohesion 1 friction 10 dilation 10 Bq-}BN?pz  
d_sigy rzvKvGd#N  
def z_pnt V8pZr+AJ  
ib = block_head 0q]0+o*%  
z_pnt = b_zone(ib) alsD TQ'  
end L)9Z Op5  
z_pnt
\IqCC h  
bound xvel 0 yvel 0 93,7yZ5#  
bound yvel velv range 0 1 0.9 1.1 \0e`sOS`L  
hist n 50 q(2ZJn13f  
hist ydisp 0 1 {=U*!`D  
hist n_sy d+ [2Sm(7  
hist a_sy ~g@}A  
label hist 1 ZC^NhgX  
Vertical displacement M[u6+`  
label hist 2
7e#|Iq:o  
Vertical stress_udec ]$-<< N{}'  
label hist 3 C/9]TkX}q  
vertical stress_analysis =<K6gC27  
step 865 "kVzN22  
程序读入后，系统提示“warning：block 202 has no contacts"