请教flac3d中分层赋予参数的不收敛问题！

科费卡

lz在做沉降的时候，根据实际情况想同一本构给予不同的参数，这里主要用的是修正剑桥模型。考虑到具体的实验数据，没有用fish来按深度赋值，而是直接手动分层赋值。但是在执行运算的过程中楼主只分了三层，就会出现*** 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).

jeaplee1

在吗，我最近也在研究修正剑桥的问题，可以一起讨论吗？

jeaplee1

我QQ11907143，你的呢？

科费卡

lz把渗流模型添加进来，想解决这个问题。结果很差，不仅是多层会出现*** 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

;---边界条件
attach face range z 191

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

科费卡

关掉了fluid biot 单层的地层的修正剑桥模型可以和渗流模型运行，但是具体的原因是什么不太清楚。

科费卡

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

科费卡

现在最大的问题是渗流模型开启，给定流体模量后，用修正剑桥模型运算不了几步就会提示*** Cam-clay:  Yield envelope cannot be reached。但是用摩尔库伦就不会出现这个问题。

科费卡

有没有人啊，求大侠帮助！