在桩顶加动载荷问题？

warlm

我刚学习FLAC3D，现在碰到1个问题，弄了好久解决不了，所以来此请教各位高手。。。

问题描述：模拟桩在循环载荷下下沉。
              建立了桩和土的模型，初始平衡已经OK,在动力分析时出现问题，桩的位移情况根本不对。

动力分析的代码如下，各位看看有什么问题，万分感谢了。。。


;动力计算阶段
rest 1111.sav       加载模型
config dyn
set dyn on
;-------------------------------------------定义正弦力
def setup
freq = 17
  omega = 2.0 * pi * freq
old_time = clock
end
setup
def wave1
  wave1 =20000+5e12*cos(omega * dytime)
end
;-------------------------------------------

free x y z                                     ;去掉模型底部原有的静力条件
apply nquiet squiet dquiet ran x  7.9 8.1  ;静态边界条件
apply nquiet squiet dquiet ran x -8.1 -7.9
apply nquiet squiet dquiet ran y -0.1 0.1
apply nquiet squiet dquiet ran y  7.9  8.1
;apply nquiet squiet dquiet ran z -0.1 0.1
apply nquiet squiet dquiet ran z -8.1 -7.9

sel set damp combined
;set dyn damp local 0.6
;ini state 0
ini xdis 0.0 ydis 0.0 zdis 0.0                     ;位移清零
apply szz 1.0 hist wave1 ran z 0.99 1.01 x -0.3 0.3 y 0 0.3 group pile;加动力荷载  （桩顶的坐标为0 0 -1）

set dyn time = 0                         ;设置动力计算从0s开始
hist id=3 dytime
hist id=4 gp zdis 0 0 0
hist id=5 gp zdis 0 0 -7
hist id=6 zone szz 0 0 0
hist id=7 zone szz 0 0 -2
hist id=8 zone szz 0 0 -8
hist id=9 zone szz 0 0 -8.1
hist id=10 zone szz 0.5 0 -6
hist gp zdis  0 0 0

set dyn multi on  ;动态多步
set hist_rep 1
solve age 0.2
save xwc900.sav

fg888134

我也在苦苦的搞这个问题呢，就是初始条件都有点不平衡下来，苦啊 求交流 515363568