求助 [复制链接]

小弟想做一圆柱形试件外加约束，试件和约束之间设置接触面，可是写完之后就是无法运行，麻烦高手帮小弟看看是什么错误，耽误您的宝贵时间，先向您道歉。 ea>\.D-S  
new wR$8drn]Rq  
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   >!.9g  
gen zone ref nor 0 1 0 origin 0 0 0 ah9',((!  
gen zone ref nor 1 0 0 origin 0 0 0 =T}uQ$X  
group yanshishijian 6|cl`}g_j  
model ss 0w. _}Cz  
prop bulk 1.19e10 shear 1.1e10 t+J
6P)=  
prop coh 2.72e5 fric 44 ten 2e5 {~I_rlo n  
prop ctab 1 ftab 2 Wj=ex3K3u.  
table 1 0,2.72e5 1e-4,2e5 2e-4,1.5e5 3e-4,1.03e5 1,1.03e5 Z x%@wH~  
table 2 0,44 1e-4,42 2e-4,40 3e-4,38 1,38 #zs\Z]3#  
int 1 face range cyl end1 0 0 110 end2 0 0 209.9 r 24 r 26 fr2w k}/b  
int 1 prop ks 2e3 kn 2e3 fric 15 l8Qi^<i/  
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 E2kRt'~N  
gen zone ref nor 0 1 0 origin 0 0 0 range z 0 110 J3^ZPW  
gen zone ref nor 1 0 0 origin 0 0 0 range z 0 110 ,x.)L=Cx8  
group yueshu range z 0 110 a\_,_psK  
model el d( g_y m*  
prop bulk 1.67e11 shear 7.7e10 range group yanshishijian not S]=Vr%irX  
fix x y z range z 9.9 10.1 7e[\0:Z  
ini z add -100 range group yanshishijian  F]#fl%  
ini zvel -1e-7 range group yanshishijian z 109.9 110.1 1tz .e\  
hist gp ydisp 0 26 55 gSYX@'Q!  
hist zone syy 0 26 55 1u+(rVQN  
hist zone szz 0 0 50 R_n-&d'PP  
hist gp zdisp 0 0 50 (r\h dLX  
set large