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

title ~"-+BG(5  
Oedometric test phi=10 psi=10 VNKtJmt  
round 0.01 yu8xTh$:  
block 0 0 0 1 1 1 1 0 0N02E  
gen quad 2 yhnhORSY;  
prop mat 1 dens 1 (80 Tbi~+  
;change cons 3 r9:Cq  
;prop mat 1 bu 200 sh 200 co 1 fric 10dil 10 Th&-n%r9K  
zone m m .{,
PC  
zone bulk 200 sh 200 co 1 fric 10 Bn>"lDf,  
zone dil 10 [k)xn3[  
;zone dil 0 pU4k/v555;  
step 0 ]ADj9  
def d_sigy |#hj O3  
c_k = bulk_mod ""-#b^DQ  
c_g = shear_mod ktQMkEj#  
e1 = c_k+4.0*c_g/3.0 WDznhMo  
e2 = c_k-2.0*c_g/3.0 ?"oW1a\  
sf = friction * degrad "[W${q+0x  
nf = sin(sf) bvVEV  
nf = (1.0+nf)/(1.0-nf) ,l}mCY  
sp = dilation*degrad ]&l.-0jt  
np = sin(sp) ID+o6/V8  
np = (1.0+np)/(1.0-np) \KmgFyF  
r1 = (e1-e2*nf)/((e1+e2)*nf*np-2.0*e2*(nf+np)+2.0*e1) a%R'x]  
v_dt = -1.0e-5 b5r.N1ms  
velv =v_dt/tdel 3Yu1ZuIR  
dsigy = v_dt*(e1+2.0*r1*(e2*np-e1)) frB~ajXK  
step1 = -2.0*cohesion*sqrt(nf)/((e1-e2*nf)*v_dt) Z/#l~.o[  
end e)
Q{yO  
def esigy .9r+LA{  
while_stepping (sX=#<B%  
if step < step1 then .Z[Bz7  
   a_sy = a_sy+e1*v_dt 6 <&jY  
else <Uz~V;  
   a_sy = a_sy+dsigy h!c6]D4!L  
end_if 0MV^-M   
   n_sy = z_syy(z_pnt) ?|Q[QP  
end #9HQW:
On  
set bulk_mod 200 shear_mod 200 if|j)h&  
set cohesion 1 friction 10 dilation 10 a(|YLN  
d_sigy @|yRo8|  
def z_pnt *R
M?SE6;  
ib = block_head _,S L;*G4|  
z_pnt = b_zone(ib) Y8*k18~  
end 3Zy$NsY3  
z_pnt {Pi+VuLE  
bound xvel 0 yvel 0 ] qT\z<}  
bound yvel velv range 0 1 0.9 1.1 ,k% \f]a  
hist n 50 =P* YwLb  
hist ydisp 0 1 \tL9`RKpg  
hist n_sy cQ:Y@f 9  
hist a_sy HiTn5XNf  
label hist 1 PQ}q5?N  
Vertical displacement xR9<I:^&  
label hist 2 /qf2LO'+  
Vertical stress_udec a[Q\8<  
label hist 3 `R}q&|o7<  
vertical stress_analysis +4[L_  
step 865 >.P/fnvJ  
程序读入后，系统提示“warning：block 202 has no contacts"