小弟想做一圆柱形试件外加约束,试件和约束之间设置接触面,可是写完之后就是无法运行,麻烦高手帮小弟看看是什么错误,耽误您的宝贵时间,先向您道歉。 *"sDsXo- I
new |"k+j_/+
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 Yjix]lUXVf
gen zone ref nor 0 1 0 origin 0 0 0 FWuk@t[<O
gen zone ref nor 1 0 0 origin 0 0 0 nHIW_+<Mf
group yanshishijian (kZ2D
model ss fC!+"g55
prop bulk 1.19e10 shear 1.1e10 Q#sLIZ8=
prop coh 2.72e5 fric 44 ten 2e5 ]L9$JTGF`w
prop ctab 1 ftab 2 ~k/GmH
table 1 0,2.72e5 1e-4,2e5 2e-4,1.5e5 3e-4,1.03e5 1,1.03e5 _rfGn,@BH
table 2 0,44 1e-4,42 2e-4,40 3e-4,38 1,38 );*#s~R
int 1 face range cyl end1 0 0 110 end2 0 0 209.9 r 24 r 26 Rx.dM_S
int 1 prop ks 2e3 kn 2e3 fric 15 Ui05o7xg~p
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 _2nNCu (
gen zone ref nor 0 1 0 origin 0 0 0 range z 0 110 IW- BY =C
gen zone ref nor 1 0 0 origin 0 0 0 range z 0 110 K)TMr"j\
group yueshu range z 0 110 NEcE-7aT
model el fz`\-"f]
prop bulk 1.67e11 shear 7.7e10 range group yanshishijian not @M8vPH
fix x y z range z 9.9 10.1 "[wP1n!G
ini z add -100 range group yanshishijian "oo
j;
ini zvel -1e-7 range group yanshishijian z 109.9 110.1 $D~vuA7
hist gp ydisp 0 26 55 nVv=smVOt
hist zone syy 0 26 55 :DR}lOi`
hist zone szz 0 0 50 QHbjZJ
N
hist gp zdisp 0 0 50 d(;Qe}ok>
set large