关于高速公路交通荷载模拟问题求教

轩辕沥泉

最近在做毕业设计，需要模拟高速公路受交通荷载时的应力情况，根据陈育民老师那本书第七章的例子建立了路堤模型，模拟了路堤施工过程，然后设置了自由边界条件，施加了一个非常简单的正弦波动的荷载，但是出来的位移时间曲线和速度时间曲线毫无规则，请帮忙看看是哪里出了问题，是参数设置错误，还是模型建立过程不对，谢谢了。。。。。。。下面是程序
new
config dyn
set dyn off
gen zone brick p0 0 0 -2.1 p1 32.5 0 -2.1 p2 0 5 -2.1 p3 0 0 0 size 32 1 2 group hard
gen zone brick p0 32.5 0 -2.1 p1 100 0 -2.1 p2 32.5 5 -2.1 p3 32.5 0 0 size 67 1 2 group hard
gen zone brick p0 0 0 -3.9 p1 32.5 0 -3.9 p2 0 5 -3.9 p3 0 0 -2.1 size 32 1 2 group clay
gen zone brick p0 32.5 0 -3.9 p1 100 0 -3.9 p2 32.5 5 -3.9 p3 32.5 0 -2.1 size 67 1 2 group clay
gen zone brick p0 0 0 -10.9 p1 32.5 0 -10.9 p2 0 5 -10.9 p3 0 0 -3.9 size 32 1 7 group clay_1
gen zone brick p0 32.5 0 -10.9 p1 100 0 -10.9 p2 32.5 5 -10.9 p3 32.5 0 -3.9 size 67 1 7 group clay_1
gen zone brick p0 0 0 0 p1 32.5 0 0  p2 0 5 0 p3 0 0 3 p4 32.5 5 0 p5 0 5 3 p6 28 0 3 p7 28 5 3 size 32 1 3 group dam
;边界条件
fix x y z ran z -11.0 -10.8
fix x ran x -.1 .1
fix x ran x 99.9 100.1
fix y
;初始应力
model mohr ran z -10.9 0
model null ran z 0 3
prop bulk 6.89e6 shear 3.55e6 coh 10e10 tension 1e10 ran group hard
ini dens 1863 ran group hard
prop bulk 2.19e6 shear 0.78e6 coh 10e10 tension 1e10 ran group clay
ini dens 1710 ran group clay
prop bulk 5.09e6 shear 2.21e6 coh 10e10 tension 1e10 ran group clay_1
ini dens 1864 ran group clay_1
set grav 0 0 -9.8
solve
prop bulk 6.89e6 shear 3.55e6 coh 21e3 fric 30 ran group hard
prop bulk 2.19e6 shear 0.78e6 coh 6e3 fric 13 ran group clay
prop bulk 5.09e6 shear 2.21e6 coh 16e3 fric 32 ran group clay_1
solve
;施工过程模拟
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
model elastic ran z 0 1
prop bulk 7.8e6 shear 3.0e6 ran z 0 1
ini dens 1500 ran z 0 1
solve
model elastic ran z 1 2
prop bulk 7.8e6 shear 3.0e6 ran z 1 2
ini dens 1500 ran z 1 2
solve  
model elastic ran z 2 3
prop bulk 7.8e6 shear 3.0e6 ran z 2 3
ini dens 1500 ran z 2 3
solve
;动力荷载施加
set dyn on
def wave
   wave=sin(2*pi*dytime)
end
free x y z ran z -11.0 -10.8
apply nquiet squiet dquiet ran z -11.0 -10.8
apply dstress -2e6 hist wave ran z 2.9 3.1 x 0.79 1.03 y 1.2 1.5
apply ff
set dyn time=0
hist reset
hist id=1 unbal
hist id=2 dytime
hist id=3 gp zdisp  0.91 1.35 3
set hist_rep 1
plot add hist  3 vs 2
plot show
solve age 1.0