模拟砂土液化

GZhang2468

动力计算是网格变形过大了，并且动力时间不动了，我在模型四周设置了自由场边界，底部设置了静态边界，并在底部施加应力时程动荷载。请大佬们答疑解惑，是我的动态边界不对还是动荷载施加的不对
rest 1.sav
set dyn on fluid on
ini fmod 2e8
set fluid pcut on;防止负孔隙水压力出现
model finn ran gro sand
prop bulk = 3e7 shear = 1e7 fric =35 ran gro sand
ini dens 2000 ran gro 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 gro sand
;
def setup
    freq =5.0
    ampl =2
    omega =2.0*pi*freq
end
@setup
;
def sine_wave
  vv=9.36e-2*sin(omega*dytime)
  if dytime <2.0
     sine_wave = dytime/2.0*vv
     else
     if dytime < 20.0
        sine_wave = vv
        else
        if dytime <= 30.0
           sine_wave = (30.0-dytime)/10.0*vv
        endif
      endif
  endif
  if dytime > 30.0
      sine_wave =0.0
  endif
end
;
free  x y  z range z -9.9 -10.1
apply dquiet squiet range z -9.9 -10.1
apply dstress =-2.83e5  hist @sine_wave ran z -9.9 -10.1
apply ff
set dyn damp local 0.314  
;
define abc
     pzm1 = z_near(15,0,-5)
end   
@abc
;
def eff_stress_m1
      tmp1 =  z_szz(pzm1)
      eff_stress_m1 = tmp1+ z_pp(pzm1)
end
@eff_stress_m1
;
hist add dytime
hist add fish @eff_stress_m1
hist add zone pp 15 0 -5
hist add fish @sine_wave
;
;set hist_rep 100
set large
set dyn dt 3e-4
set mech rat 1e-50
;
solve age 30
save save_file.sav

QCC

楼主 解决了么 我也是做这一块的 我发现动荷载下网格破坏很明显 我这边是低Hz下 感觉在晃水一样

GZhang2468

我把大变形模式关了就没事了

单一一

GZhang2468 发表于 2020-12-26 22:05
我把大变形模式关了就没事了

您好，您的这个问题解决了吗，我也是这个问题。qq, 2748375059,希望和您交流下

GZhang2468

单一一 发表于 2021-3-10 15:16
您好，您的这个问题解决了吗，我也是这个问题。qq, 2748375059,希望和您交流下 ...

我就是把大变形模式给关了