n
conf dyn
gen zone radcyl p0 (0,0,0) p1 (9.5,0,0) p2 (0,0,-14.5) p3 (0,9.5,0) &
p4 (9.5,0,-14.5) p5 (0,9.5,-14.5) p6 (9.5,9.5,0) p7 (9.5,9.5,-14.5) &
p8 (0.5,0,0) p9 (0,0.5,0) p10 (0.5,0,-14.5) p11 (0,0.5,-14.5) &
size 5 30 20 15 ratio 1 1 1 1.2
gen zone radcyl p0 (0,0,-14.5) p1 (9.5,0,-14.5) p2 (0,0,-15) p3 (0,9.5,-14.5) &
p4 (9.5,0,-15) p5 (0,9.5,-15) p6 (9.5,9.5,-14.5) p7 (9.5,9.5,-15) &
p8 (0.5,0,-14.5) p9 (0,0.5,-14.5) p10 (0.5,0,-15) p11 (0,0.5,-15) &
size 5 26 20 15 ratio 1 1 1 1.2
gen zone radcyl p0 (0,0,-15) p1 (9.5,0,-15) p2 (0,0,-30) p3 (0,9.5,-15) &
p4 (9.5,0,-30) p5 (0,9.5,-30) p6 (9.5,9.5,-15) p7 (9.5,9.5,-30) &
p8 (0.5,0,-15) p9 (0,0.5,-15) p10 (0.5,0,-30) p11 (0,0.5,-30) &
size 5 31 20 15 ratio 1 1 1 1.2 fill
group tuceng range z -30 0
;
interface 1 face range cylinder end1 (0,0,0) end2 (0,0,-15) radius 0.51 &
cylinder end1 (0,0,0) end2 (0,0,-15) radius 0.49 not
interface 2 face range cylinder end1 (0,0,-14.9) end2 (0,0,-15) radius 0.51
;
gen zone cyl p0 (0,0,5.5) p1 (0.5,0,5.5) p2 (0,0,5) p3 (0,0.5,5.5) &
p4 (0.5,0,5) p5 (0,0.5,5) &
size 5 26 20
gen zone cyl p0 (0,0,20) p1 (0.5,0,20) p2 (0,0,5.5) p3 (0,0.5,20) &
p4 (0.5,0,5.5) p5 (0,0.5,5.5) &
size 5 30 20
gen zone cyl p0 (0,0,20.1) p1 (0.5,0,20.1) p2 (0,0,20) p3 (0,0.5,20.1) &
p4 (0.5,0,20) p5 (0,0.5,20) &
size 5 1 20
group zhuang range z 5 20.1
;
ini z add -20 range group zhuang
group chenzha range cylinder end1 (0,0,-14.5) end2 (0,0,-15) radius 0.5 group zhuang
;
model elastic range group zhuang
prop bulk 8.3e7 shear 3.8e7 range group zhuang
model mohr range group chenzha
prop bulk 4.17e6 shear 1.92e6 coh 2.5955e5 &
fric 8.73 range group chenzha
model mohr range group tuceng
prop bulk 8.3e7 shear 3.8e7 coh 4.5e5 &
fric 28 range group tuceng
interface 1 prop kn 2.67e8 ks 2.67e8 fric 22.4 coh 3.6e5
interface 2 prop kn 6.2e8 ks 6.2e8 fric 22.4 coh 3.6e5
;
ini dens 2200 range group tuceng
ini dens 2000 range group zhuang
ini dens 740 range group chenzha
model null range z -0.1 0.15
;
fix z range z -29.9 -30.1
fix x range x -0.1 0.1
fix x range x 9.4 9.6
fix y range y -0.1 0.1
fix y range y 9.4 9.6
set gravity 0 0 -10
ini szz 0. grad 0 0 21000 ran z -30 0
ini sxx 0. grad 0 0 12600 ran z -30 0
ini syy 0. grad 0 0 12600 ran z -30 0
hist unbal
solve rat 1.e-5
save zhuang0.sav
def setup
omega=2.0*pi*freq
pulse=1.0/freq
end
set freq=4.0
setup
def wave
if dytime>pulse
wave=0.0
else
wave=0.5*(1.0-cos(omega*dytime))
endif
end
range name bottom z=-0.1 0.1
fix z range z=0.5 35
apply dquiet squiet range bottom
apply sxz -9e5 hist wave syz 0.0 szz 0.0 range bottom
apply nvel 0 plane norm 0,0,1 range bottom
hist gp xvel 0,0,0
hist gp xvel 0,0,15
hist gp xvel 0,0,30
hist dytime
hist wave
plot create hhh
plot add hist 1 2 3 vs 4
plot show
solve age 2
各位老师帮我看一下有什么错误吗 我出来的图形有点不对诶 非常感谢!