初始分析：
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
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
