flac2D中fish剑桥模型参数问题

wuyong505

fish目的是定义前期固结压力mpc。如下：

但是结果出现提示：divide by zero
通过print fish 发现：

应该是mm的问题。
请问在flac2d的fish里参数mm如何表示？m_m、cam_mm都试过了。

土体是分层的，砂土采用mc模型，粘土是mcc模型。典型上海土层参数取值。i=1，j=12是第一个MCC模型的grid
。。。。。。。。。。。
通过直接给mm赋值可以解决，但是solve时出现：

为什么啊？？

junfeng2009

好问题  有没有高人指点

wym247

很好 学习下

mazhongyuan9049

fish刚开始的几句就我的理解应该是利用单元内的应力计算偏应力和平均应力，然后将单元的平均应力的值赋给先期固结压力参数mpc，但是如果没有初始应力的话单元内的应力为零，所以先期固结压力就为零值了，所以计算就会出错。