地震液化条件下的桩基数值模拟

woaishasha

各位大神能不能帮我看看建好模型之后的静力计算和动力计算有啥问题吗，我输出的加速度，位移曲线这些总感觉不太对劲
;静力分析
model elas
prop shear=1e7 bulk=3e7
ini dens 1850
model fl_iso
prop poro 0.5 perm 1e-11
ini fmod 2e8 fdens 1000
ini pp 0 grad 0 0 -10e3 ran z 0 -20
ini syy -1.25e5 grad 0 0 2.5e4 ran z 0 -20
ini sxx -1e5 grad 0 0 2e4 ran z 0 -20
ini szz -1e5 grad 0 0 2e4 ran z 0 -20
fix z ran z -0.1 0.1
fix x ran x -50.1 -49.9
fix x ran x 49.9 50.1
fix y ran y -50.1 -49.9
fix y ran y 49.9 50.1
set grav 0 0 -10
set fluid off dyn off
ini fmod 0
set mech rat 1e-5
solve
def ini_conf
    _k0=1
    pnt=zone_head
    loop while pnt#null
         val=_k0*z_szz(pnt)-(_k0-1)*z_pp(pnt)
         z_sxx(pnt)=val
         z_syy(pnt)=val
         pnt=z_next(pnt)
    end_loop
end
@ini_conf
solve

;动力计算设置模型参数
set dyn on fluid on
ini fmod 2e8 fden 1000
set fluid pcut on
model finn ran group sand
prop poro 0.5 perm 1e-7 ran group sand
prop shear=1.92e7 bulk=3.2e7 fric=35 cohe=0 ran group sand
ini dens 2080 ran group sand
prop ff_latency=50 ff_switch=0 ff_c1=0.8 ff_c2=0.79 ff_c3=0.45 ff_c4=0.73 ran group sand

ini xvel 0 yvel 0 zvel 0
ini xdis 0 ydis 0 zdis 0
ini state 0

;施加自由场边界条件
apply ff
group ff_corner
group ff_side ran x -50 50
group ff_side ran y -50 50
group main_grid ran x -50 50 y -50 50

;输入0.05g正弦波
def setup
      freq=1
      omega=2*pi*freq
      old_time=clock
end
@setup
def wave
      wave=0.5*sin(omega*dytime)
end
free x
apply xacc=1 hist @wave range z=-20.1 -19.9
apply xacc=1 hist @wave range x=-50.1 -49.9
apply xacc=1 hist @wave range x=49.9 50.1


;局部阻尼
set dyn damp local 0.314


def constants
p_z1=z_near(0.5,0,0)
p_z2=z_near(0.5,0,-5)
p_z3=z_near(0.5,0,-10)
p_z4=z_near(0.5,0,-15)
end
@constants

def eff_stress1
  eff_stress1=(z_sxx(p_z1)+z_syy(p_z1)+z_szz(p_z1))/3+z_pp(p_z1)
end

def eff_stress2
  eff_stress2=(z_sxx(p_z2)+z_syy(p_z2)+z_szz(p_z2))/3+z_pp(p_z2)
end

def eff_stress3
  eff_stress3=(z_sxx(p_z3)+z_syy(p_z3)+z_szz(p_z3))/3+z_pp(p_z3)
end

def eff_stress4
  eff_stress4=(z_sxx(p_z4)+z_syy(p_z4)+z_szz(p_z4))/3+z_pp(p_z4)
end


hist zone pp 0.5 0 0
hist zone pp 0.5 0 -5
hist zone pp 0.5 0 -10
hist zone pp 0.5 0 -15
hist eff_stress1
hist eff_stress2
hist eff_stress3
hist eff_stress4
hist gp xvel 0.5 0 0
hist gp xvel 0.5 0 -5
hist gp xvel 0.5 0 -10
hist gp xvel 0.5 0 -15
hist gp xacc 0.5 0 0
hist gp xacc 0.5 0 -5
hist gp xacc 0.5 0 -10
hist gp xacc 0.5 0 -15
hist gp xdis 0.5 0 0
hist gp xdis 0.5 0 -5
hist gp xdis 0.5 0 -10
hist gp xdis 0.5 0 -15
hist dytime
def tim
      tim=0.01*(clock-old_time)
end
set dyn multi on
solve age 5
print dyn
save mult1.sav

flac3dlaji

留个联系方式讨论一下呗