请教flac3d中分层赋予本构的不收敛问题！ [复制链接]

lz在做沉降的时候，根据实际情况想对不同的地层想给予不同的本构，或者是同一本构给予不同的参数，这里主要用的是修正剑桥模型。但是在执行运算的过程中，会出现*** Cam-clay:  Yield envelope cannot be reached，而在单层运算的过程中不会出现这个问题。个人猜测是下层的变形快于上层，可能会出现分离的情况因此可能导致无法运行。但是调参往往解决不了问题，而且过度的调参也有失真实。请问大侠们因该如何处理这种问题？是不是这里没有给出地下水的信息，导致的有效应力过大？还是什么其他的原因？下边贴出只有三个层位的命令流，求解答！


new
; --- model geometry ---
gen z b p0 0 0 191 p1 1080 0 191 p2 0 1000 191 p3 0 0 200 size 10 10 9 ratio 1 1 1 group diceng1
gen z b p0 0 0 175.5 p1 1080 0 175.5 p2 0 1000 175.5 p3 0 0 191 size 10 10 3 ratio 1 1 1 group diceng2
gen z b p0 0 0 100 p1 1080 0 100 p2 0 1000 100 p3 0 0 175.5 size 10 10 10 ratio 1 1 1 group diceng3
attach face range z 100 200

;---初始地应力场生成
model mohr range group diceng1
prop bulk 5.271428571e6 shear 2.432967033e6 coh 1e10 fri 16 ten 1e10 range group diceng1
model mohr range group diceng2
prop bulk 6.833333333e6 shear 3.153846154e6 coh 1e10 fri 14 ten 1e10 range group diceng2
model mohr range group diceng3
prop bulk 17.96190476e6 shear 8.29010989e6 coh 1e10 fri 15.5 ten 1e10 range group diceng3

;---边界条件、初始密度
Initial density 1900 range group diceng1
Initial density 2020 range group diceng2
Initial density 2000 range group diceng3
Fix x range x=-0.1 0.1
Fix x range x=1080
fix y range y=-0.1 0.1
fix y range y=1000
fix z range z=100
set gravity 0 0 -9.7966

hist id=1 gp zdis 500 500 200
solve


ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0

;---修正剑桥模型
model cam-clay range group diceng1
Property shear=2.432967033e6 bulk_bound=1000e6 range group diceng1
Property mm=0.60705 lambda=0.0655785 kappa=0.0017372 range group diceng1
Property mpc=0.0924e6 poisson=0.3 range group diceng1
Property mp1=1.0 mv_l=1.99985195 range group diceng1
model cam-clay range group diceng2
Property shear=5.992307692e6 bulk_bound=1000e6 range group diceng2
Property mm=0.526284 lambda=0.0642756 kappa=0.0034744 range group diceng2
Property mpc=0.3e6 poisson=0.3 range group diceng2
Property mp1=1.0 mv_l=2.007570044 range group diceng2
model cam-clay range group diceng3
Property shear=8.29010989e6 bulk_bound=10000e6 range group diceng3
Property mm=0.58674355 lambda=0.061669816 kappa=0.000434294 range group diceng3
Property mpc=0.561e6 poisson=0.3 range group diceng3
Property mp1=1.0 mv_l=1.931061823 range group diceng3

ini sxx -5. syy -5. szz -7000.
def install_prop
pnt=zone_head
loop while pnt # null
  abs_sxx=abs(z_sxx(pnt))
  abs_syy=abs(z_syy(pnt))
  abs_szz=abs(z_szz(pnt))
  p0=(abs_sxx+abs_syy+abs_szz)/3.0
  p0_effective=p0-z_pp(pnt)
  q0=sqrt(((abs_sxx-abs_syy)^2+(abs_syy-abs_szz)^2+(abs_szz-abs_sxx)^2)*0.5)
  z_prop(pnt,'cam_cp')=p0_effective
  pnt=z_next(pnt)
endloop
end
@install_prop
attach face range z 100 200
Fix x range x=-0.1 0.1
Fix x range x=1080
fix y range y=-0.1 0.1
fix y range y=1000
;fix z range z=0
fix z range z=100
set gravity 0 0 -9.7966
hist id=2 gp zdis 500 500 200
solve

另外请问渗流模式中是不是所有的本构关系都可以？
为什么同样的渗流模式，所有的命令流都相同，只是一个用库伦摩尔就可以运行，而另一个修正剑桥就提示
Bad conversion of parameter number 2 (fl_iso).

你的命令流与现在使用的版本不兼容，需要修改

edelmann:你的命令流与现在使用的版本不兼容，需要修改 (2016-08-21 08:50) 

我用的是flac3d 5.0 ，如果是命令流的问题为什么在单层时候运行没有问题？如果是这样请帮忙指出具体是哪里出了问题，谢谢！

楼主把渗流模型添加进来，想解决这个问题。结果很差，不仅是多层会出现*** Cam-clay:  Yield envelope cannot be reached。在单层运行的时候也出现了类似的问题。请大神帮忙看看问题到底在哪？下面贴出命令流：
new
config fluid
set fluid biot on
; --- model geometry ---
gen z b p0 0 0 191 p1 1080 0 191 p2 0 1000 191 p3 0 0 200 size 10 10 9 ratio 1 1 1 group diceng1
;---修正剑桥模型
model cam-clay range group diceng1
Property shear=2.432967033e6 bulk_bound=1000e6 range group diceng1
Property mm=0.60705 lambda=0.0655785 kappa=0.0017372 range group diceng1
Property mpc=0.0924e6 poisson=0.3 range group diceng1
Property mp1=1.0 mv_l=1.99985195 range group diceng1
;---重力及土层初始密度
set gravity 0 0 -9.7966
Initial density 1620 range group diceng1
;---渗流模型
set fluid on mech on
model fl_iso
ini ftens 1.0e-15
ini fdens 1000
ini biot_mod 4e9
prop perm 5.92e-11 biot_c 1
;---fish
ini sxx -5. syy -5. szz -7000.
def install_prop
pnt=zone_head
loop while pnt # null
  abs_sxx=abs(z_sxx(pnt))
  abs_syy=abs(z_syy(pnt))
  abs_szz=abs(z_szz(pnt))
  p0=(abs_sxx+abs_syy+abs_szz)/3.0
  p0_effective=p0-z_pp(pnt)
  q0=sqrt(((abs_sxx-abs_syy)^2+(abs_syy-abs_szz)^2+(abs_szz-abs_sxx)^2)*0.5)
  z_prop(pnt,'cam_cp')=p0_effective
pnt=z_next(pnt)
endloop
end
@install_prop
;---边界条件
Fix x range x=-0.1 0.1
Fix x range x=1080
fix y range y=-0.1 0.1
fix y range y=1000
fix z range z=191
app pp 1924e3 grad 0 0 -10e3 range x=-0.1 0.1
app pp 1924e3 grad 0 0 -10e3 range x=1080
app pp 1924e3 grad 0 0 -10e3 range y=-0.1 0.1
app pp 1924e3 grad 0 0 -10e3 range y=1000
ini pp 1924e3 grad 0 0 -10e3
;---观测点记录
hist id=1 gp zdis 500 500 200
solve

看了一下楼主的代码，自己没有模拟，说不了什么，但是楼主算第一应力不变量为啥要取绝对值啊？公式里没说要取绝对值啊。FLAC3D里面压应力为负，拉应力为正，空隙压力反过来的，只能为正，你这样算出来的有效应力p0_effectiv方向可能要反吧，毕竟没算过这个，只是提出我的疑问 ：D

紫色毛笔:看了一下楼主的代码，自己没有模拟，说不了什么，但是楼主算第一应力不变量为啥要取绝对值啊？公式里没说要取绝对值啊。FLAC3D里面压应力为负，拉应力为正，空隙压力反过来的，只能为正，你这样算出来的有效应力p0_effectiv方向可能要反吧，毕竟没算过这个，只是提出我的疑问 ：D (2016-08-26 00:45) 

个人认为，附加应力主要是向下，岩土体中主应力是SZZ，它的方向是向下的为负，syy和sxx很小。三者之和整体的结果为负，有效应力应该是附加应力与孔隙水压力之差。因此有效应力为正，需要对主应力取绝对值。不知道我的理解有没有问题。

最近发现，土体为修正剑桥时，渗流模型如果不给定流体模量或者比奥模量，计算可以进行，计算将只进行至力学平衡。一旦给土层中的流体赋予流体模量则出现*** Cam-clay:  Yield envelope cannot be reached，具体是什么原因还不太理解，希望大侠指明！

代码里面 ini sxx -5. syy -5. szz -7000. 什么意思？计算完初始地应力接着把地应力又全改了？

紫色毛笔:代码里面 ini sxx -5. syy -5. szz -7000. 什么意思？计算完初始地应力接着把地应力又全改了？ (2016-08-26 21:05) 

这个是之前没做初始地应力场时加的，后来忘记删掉了，谢谢！现在最大的问题是渗流模型开启，给定流体模量后，用修正剑桥模型运算不了几步就会提示*** Cam-clay:  Yield envelope cannot be reached。但是用摩尔库伦就不会出现这个问题。

怎么没人解答，求大侠帮助！

科费卡:[表情] 怎么没人解答，求大侠帮助！ (2016-09-03 10:40) 

怎么没人解答，求大侠帮助！