初始分析：
impgrid deceng.flac3d
config dyn fluid
model mohr
prop bulk 3.39e7 shear 1.56e7 coh 1.4e4 fri 25 ran gro soil
prop bulk 8e8 shear 5.3e8 coh 1.5e5 fri 40 ran gro rock
ini dens 1430 ran gro soil;计算干密度
ini dens 2307 ran gro rock
model fl_iso;流体模型
prop poro 0.53 perm 7.8e-11 ran gro soil
prop poro 0.128 perm 1.9e-12 ran gro rock
ini fmod 2e8 fdens 1000
ini fmod 0
impgrid shuiwei.flac3d;建立水位面
gen separate water_below
interface 2 wrap water_below water_above
set grav 0 0 -9.81
water den 1000
def water_table
p_i=i_head
p_ie=i_elem_head(p_i)
loop while p_ie # null
p_gp1=ie_vert(p_ie,1)
p_gp2=ie_vert(p_ie,2)
p_gp3=ie_vert(p_ie,3)
x1=in_pos(p_gp1,1)
y1=in_pos(p_gp1,2)
z1=in_pos(p_gp1,3)
x2=in_pos(p_gp2,1)
y2=in_pos(p_gp2,2)
z2=in_pos(p_gp2,3)
x3=in_pos(p_gp3,1)
y3=in_pos(p_gp3,2)
z3=in_pos(p_gp3,3)
command
water table face x1,y1,z1 x2,y2,z2 x3,y3,z3
endcommand
p_ie=ie_next(p_ie)
endloop
end
model null range group water_below
model null range group water_above
water_table
int 2 dele
dele range group water_below
dele range group water_above
gen separate soil;建立滑动面接触面
int 1 wrap soil rock
int 1 prop ks 3.4e8 kn 3.4e8 fric 24 coh 1.2e4
fix x y z ran z -0.1 0.1;边界条件
fix y ran y -0.1 0.1
fix y ran y 689.9 690.1
fix x ran x -0.1 0.1
fix x ran x 0.9 1.1
set grav 0 0 -9.81
set fluid off dyn off
hist id=1 unbal
solve elas
save chushi.sav
动力分析：
restore chushi.sav
set dyn on fluid on
set dyn multi on;动态多步
set large;大变形
ini fmod 2e8
set fluid pcut on
model finn ran gro soil;finn模型
prop bulk 3.39e7 shear 1.56e7 coh=1.4e4 fric=25 ran gro soil;finn参数
ini dens 1430 ran gro soil
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 soil
ini xvel 0 yvel 0 zvel 0;位移速度设为0
ini xdis 0 ydis 0 zdis 0
ini state 0
fix x y z ran z -0.1 0.1
fix y ran y -0.1 0.1
fix y ran y 689.9 690.1
fix x ran x -0.1 0.1
fix x ran x 0.9 1.1
free y ran z -0.1 0.1
table 1 read 1.txt;y方向施加地震波
apply yacc=9.8 hist table 1 ran z -0.1 0.1
set dyn damp local 0.224;局部阻尼
apply ff;自由场边界
group ff_corner
group ff_side ran x 0 1
group ff_side ran y 0 690
group main_grid ran x 0 1 y 0 690
config zextra 20;监测动孔压数据
def get_old_stress
p_z = zone_head
loop while p_z # null
z_extra(p_z,5) = z_pp(p_z) ;初始孔压
p_z = z_next(p_z)
endloop
end
get_old_stress
def get_PPR
whilestepping
p_z = zone_head
loop while p_z # null
z_extra(p_z,7) = z_szz(p_z) + z_pp(p_z) ;竖向有效应力
z_extra(p_z,4) = z_szz(p_z) ;竖向应力
z_extra(p_z,6) = z_pp(p_z) - z_extra(p_z,5) ;超静孔隙水压力
z_extra(p_z,9) = z_szz(p_z)
z_extra(p_z,2) = (-z_extra(p_z,6) ) / (z_extra(p_z,4)+ z_extra(p_z,5)) ;超孔压比
z_extra(p_z,3) = 1.0 - z_extra(p_z,7) / (z_extra(p_z,4)+ z_extra(p_z,5)) ;超孔压比
p_z = z_next(p_z)
endloop
end
get_PPR
def get_ppr_zone
whilestepping
p_z1 = z_near( 1,3.115e2,8.102e1 );监测点
ppr_11 = z_extra(p_z1,2)
ppr_12 = z_extra(p_z1,3)
ppr_13 = z_extra(p_z1,6)
ppr_14 = z_extra(p_z1,7)
end
hist id 2 dytime
hist id 3 ppr_11
hist id 4 ppr_12
hist id 5 ppr_13
hist id 6 ppr_14
set dyn time=0
solve age 10

楼主，这个算出来了吗？

请问楼主，prop ff_latency=50 ff_switch=0
是什么意思？是不是表示 相对密度Dr?

后边的有效应力和超孔隙水压力那里有问题呀，有人算出来了吗

三维砂土液化大变形统一本构模型-张建民院士课题组开发。
包含dll文件，本构模型使用说明以及两个算例（扭剪试验和VELACS试验模拟）。
https://bbs.yantuchina.com/read.php?tid=326726

建立网格的步骤没看见