第一步地应力平衡
new
config dyn
config fluid
set dyn off
set mech off
gen zone brick size 25,25,25 &
p0(0,0,0) p1(10,0,0) p2(0,10,0) p3(0,0,10)
;边界条件
fix x range x -0.1 0.1
fix x range x 9.9 10.1
fix y range y -0.1 0.1
fix y range y 9.9 10.1
fix z range z -0.1 0.1
; 属性
gro s1 range z=0,5
gro s2 range z=5,10
model fl_iso
prop perm 1e-13 poro 0.7
ini fdens 1000.
ini sat 1.0
ini fmod 2e9
ini pp 100000 grad 0 0 -10000 ran z 0 10
apply pp 100000 grad 0 0 -10000 ran x -0.1 0.1
apply pp 100000 grad 0 0 -10000 ran x 9.9 10.1
set grav 0 0 -10
hist unb
set mech on fluid off
ini fmod 0
model mohr
prop bulk 6.67e6 shear 3.08e6 dens 1800 fric 10 coh 1.0e4 tens 1e10 range gro s1
prop bulk 8.33e6 shear 3.85e6 dens 2000 fric 15 coh 1.2e4 tens 1e10 range gro s2
ini szz 0 grad 0 0 1.8e4 range z 5 10
ini szz 1e4 grad 0 0 2.0e4 range z 0 5
ini syy 0 grad 0 0 9e3 range z 5 10
ini syy 5e3 grad 0 0 1e4 range z 0 5
ini sxx 0 grad 0 0 9e3 range z 5 10
ini sxx 5e3 grad 0 0 1e4 range z 0 5
plot hist 1
pl set plane ori 0 5 0 nor 0 1 0
pl con szz
;pl con pp plane
;pl add flow plane
solve
save geo.sav
第二步动力部分
rest geo.sav
config dyn
ini state 0
ini xvel 0 yvel 0 zvel 0
ini xdisp 0 ydisp 0 zdisp 0
set dyn on mech on fluid on
set fluid time 0
set fluid pcut on
model finn
prop bulk 6.67e6 shear 3.08e6 dens 1800 fric 10 coh 1.0e4 tens 1e10 range gro s1
prop bulk 8.33e6 shear 3.85e6 dens 2000 fric 15 coh 1.2e4 tens 1e10 range gro s2
ini dens 1800 range gro s1
ini dens 2000 range gro s2
ini fmod 2e9
prop ff_latency=50 ff_switch=0 ff_c1=0.8 ff_c2=0.79 ff_c3=0.45 ff_c4=0.73
def dy_load
if dytime<0.05
dy_load=4.55*dytime/0.05
else
if dytime<0.1
dy_load=4.55*(0.1-dytime)/0.05
endif
endif
if dytime>0.1
dy_load=0
endif
end
free x y z ran z -0.1 0.1
apply nquiet squiet dquiet range z -0.1 0.1
apply szz -1e7 hist dy_load range z 9.9 10.1 x 4 6 y 4 6
apply ff
group ff_conner
group ff_side range x 0 10
group ff_side range y 0 10
group main_grid range x 0 10 y 0 10
set dyn time=0
set dyn dt=2e-4
set dyn damp rayleigh 1 22.8
hist gp zdis 5 5 10
hist gp pp 5 5 10
hist dytime
;plot hist 3 vs 4
solve age 1.1
plot set plane ori 0 5 0 norm 0 1 0
plot con pp plane
sav qianghang.sav
计算完成后的地应力阶段正确，动力部分输出的孔压不对
求大神指教
