求助 [复制链接]

小弟想做一圆柱形试件外加约束，试件和约束之间设置接触面，可是写完之后就是无法运行，麻烦高手帮小弟看看是什么错误，耽误您的宝贵时间，先向您道歉。 XOxm<3gXn  
new [ML4<Eb+x  
gen zone cyl p0  0 0 110 p1 0 25 110 p2 0 0 210 p3 25 0 110 p4 0 25 210 p5 25 0 210   XVY^m}pMe  
gen zone ref nor 0 1 0 origin 0 0 0 A
/'G.H  
gen zone ref nor 1 0 0 origin 0 0 0 Q (`IiV   
group yanshishijian %o5GD  
model ss 9AS,-5;XQ  
prop bulk 1.19e10 shear 1.1e10 tJ 6:$dh  
prop coh 2.72e5 fric 44 ten 2e5 >~6 ;9{@  
prop ctab 1 ftab 2 k(s3~S2h  
table 1 0,2.72e5 1e-4,2e5 2e-4,1.5e5 3e-4,1.03e5 1,1.03e5 HzWZQ6o  
table 2 0,44 1e-4,42 2e-4,40 3e-4,38 1,38 |u03~L9G  
int 1 face range cyl end1 0 0 110 end2 0 0 209.9 r 24 r 26 FC(m)S2  
int 1 prop ks 2e3 kn 2e3 fric 15 [}FP_Su$6  
gen zone cs p0 0 0 10 p1 0 35 10 p2 0 0 110 p3 35 0 10 dim 25 25 25 25 size 4 10 10 10 62.{8Uj  
gen zone ref nor 0 1 0 origin 0 0 0 range z 0 110 ,Qp58u2V  
gen zone ref nor 1 0 0 origin 0 0 0 range z 0 110 g|W~0A@D  
group yueshu range z 0 110 xe 6x!  
model el 2(UT;PSI  
prop bulk 1.67e11 shear 7.7e10 range group yanshishijian not k%hif8y  
fix x y z range z 9.9 10.1 9!o:)99U  
ini z add -100 range group yanshishijian A M# '(k(  
ini zvel -1e-7 range group yanshishijian z 109.9 110.1 8,%y`tUn>u  
hist gp ydisp 0 26 55 %!WQ;(  
hist zone syy 0 26 55 ~=}56yxl[  
hist zone szz 0 0 50 7^S&g.A  
hist gp zdisp 0 0 50 vq x;FAqZ  
set large