bishan
:估计是把边破算屈服吧
(2013-09-23 19:21) 
fix x ran x -0.1 0.1
fix x ran x 99.9 100.1
fix x y z ran y -0.1 0.1
fix z ran z -0.1 0.1
fix z ran z 5.9 6.1
model elas
group 1 ran y -0.1 32.1
prop bulk 7e8 shear 3e8 den 1800
group 2 ran y 31.9 35.1
prop bulk 7e8 shear 3e8 den 1720
group 3 ran y 34.9 50.1
prop bulk 6.5e8 shear 3e8 den 1710
set grav 0 -10 0
solve
ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0
model mohr
group 1 ran y -0.1 32.1
prop bulk 7e8 shear 3e8 den 1800 coh 42100 fric 23.6
group 2 ran y 31.9 35.1
prop bulk 7e8 shear 3e8 den 1720 coh 31600 fric 21.7
group 3 ran y 34.9 50.1
prop bulk 6.5e8 shear 3e8 den 1710 coh 30500 fric 21.5
