我刚学习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