渗流模式下不收敛问题
在渗流模式下，关闭了流体计算，设置ini fmod 0 防止变形改变孔压，那为什么在用摩尔库伦模型进行地应力平衡计算的时候，就不收敛了求大神帮助啊，我已经卡了好久了，建桩之后就不收敛了，是边界设错了，还是少付了什么参数啊
model elas
prop density 1360 bulk 1.25e6 shear 1.67e6 range group top1
prop density 1460 bulk 2.5e6 shear 3.34e6 range group top2
prop density 1500 bulk 2.5e6 shear 3.34e6 range group top3
prop density 1580 bulk 5e6 shear 6.67e6 range group top4
prop density 1660 bulk 5e6 shear 2.34e7 range group top5
prop density 1700 bulk 7.14e6 shear 3.34e7 range group top6
prop density 1700 bulk 1e8 shear 4e7 range group gas
prop density 1700 bulk 1e8 shear 4e7 range group ungas
prop density 4500 bulk 3e13 shear 1e13 range group bottom
prop density 7890 bulk 10.1e10 shear 8.55e10 range group well
prop density 7890 bulk 10.1e10 shear 8.55e10 range group upwell
model fl_iso
ini fdens 1000
ini fmod 2e6
ini ftens -5e5
ini sat 1
prop por 0.4 perm 10e-17 range group gas
prop por 0.4 perm 10e-17 range group ungas
prop por 0.4 perm 10e-18 range group top1
prop por 0.4 perm 10e-18 range group top2
prop por 0.4 perm 10e-18 range group top3
prop por 0.4 perm 10e-18 range group top4
prop por 0.4 perm 10e-18 range group top5
prop por 0.4 perm 10e-18 range group top6
model fl_null range group bottom
model fl_null range group well
model fl_null range group upwell
fix x y z range z -21 21
fix y range y -400.1 -399.9
fix y range y 399.9 400.1
fix x range x -400.1 -399.9
fix x range x 399.9 400.1
set grav 0 0 -10
water dens 1000
water table ori 0 0 1300 norm 0 0 1
fix pp range group top1
fix pp range y -400.1 -399.9
fix pp range y 399.9 400.1
fix pp range x -400.1 -399.9
fix pp range x 399.9 400.1
hist unbal
plo add hist 1
set mech on fluid off
ini fmod 0
set mech ratio 1e-4
solve
save 初始应力场.sav
ini xdis=0 ydis=0 zdis=0
ini xvel=0 yvel=0 zvel=0
model mohr range group top1 any group top2 any group top3 any group top4 any group top5 any group top6 any
prop density 1360 bulk 1.25e6 shear 1.67e6 fric 0 coh 1e4 tens 2e5 range group top1
prop density 1460 bulk 2.5e6 shear 3.34e6 fric 0 coh 2e4 tens 3e5 range group top2
prop density 1500 bulk 2.5e6 shear 3.34e6 fric 0 coh 3e4 tens 5e5 range group top3
prop density 1580 bulk 5e6 shear 6.67e6 fric 0 coh 5e4 tens 9e5 range group top4
prop density 1660 bulk 5e6 shear 2.34e7 fric 0 coh 6e4 tens 1.1e6 range group top5
prop density 1700 bulk 7.14e6 shear 3.34e7 fric 0 coh 8e4 tens 1.5e6 range group top6
model mohr range group gas any group ungas any
prop density 1700 bulk 1e8 shear 4e7 fric 10 coh 6e5 tens 3e5 range group gas any group ungas any
hist unbal
plo add hist 1
set mech on fluid off
ini fmod 0
set mech ratio 1e-4
solve
save 未加载
ini xdis=0 ydis=0 zdis=0
ini xvel=0 yvel=0 zvel=0
sel pile begin 0 0 276 end 0 0 10 nseg 53
sel pile prop density 7890 emod=1.04e11nu=0.3 xcarea=0.1257xcj=0xciy=0 xciz=0 per=1.26&
cs_sk=1.43e10 cs_scoh=0 cs_sfric=0 &
cs_nk=1.43e10 cs_ncoh=0cs_nfric=0 cs_ngap=off
sel shell range cylinder end1 0 0 275 end2 0 0 185 radius 0.51
sel shell prop iso=(1.04e11,0.3) thick=0.01
app nstress -3.2e6 range z 275.9 276.1 group upwell
set mech on fluid off
ini fmod 0
set mech ratio 1e-4
solve
save 建桩后初始应力场.sav 自己顶一下
页:
[1]