案例:

研究水深10m、波高3m的波浪对海床内部孔压的影响。

模型:

使用Abaqus建立了长150m,沉积物厚度为

fssidata.txt
139kB

20m的矩形海床,分为两部分,海床最上方为厚度0.5m的一层,下面为19.5m的第二层。

时间步及材料设置:

step1 计算1s ,两个材料均设置为Elastic,Initial Stress Tensile = Yes;Step2 计算1s,表层设置为Elastic不变,下层设置为General Elastic,Initial Stress Tensile = No; Step3 计算20s,表层为Elastic,下层为PZ3模型,Initial Stress Tensile = No。

水动力、求解器设置:

step1和2 均为1阶斯托克斯波,波高为0,static;step3为3阶斯托克斯波,波高为3m,consolidation。
Elements:30900,Nodes:63210.

计算报错信息如下:

An Error Has Occured In The Running Of The Program: The Message Is:
Stop Due To Negative Hmodu+Cons1 In Pastor_Zienkiewicz_Mark_III
Stop Due To Negative Hmodu+Cons1 In Pastor_Zienkiewicz_Mark_III
Please Enter Termination Code:0: Stop Without Invoking The Debugger
1: Invoke The Debugger

Avec3: -3.900295687173900E-011 -8.449511502541048E-011 1.234980718971495E-010
8.397110637822107E-012
Hmodu: -6610064037.56690
Cons1: 4272872985.42455

$$$$$ Mdivt Exceeded: 1

怀疑是PZ3模型的拉伸问题导致的,但不确定。请问具体是什么原因造成的?需要如何调整?

(附件为FssiData)

我也有这个问题,你设置材料的时候关拉伸了吗

  • Du 回复了此帖

    Du 我刚发现之前的材料更新没关,不知道是不是这个问题

    你好,我检查了文件,发现两个问题:
    (1)边界条件没加上,尽管左侧列表已经显示施加x方向固定边界条件,但是通过右侧边界条件可视化发现,该边界条件没有成功施加。这可能是由于网格节点和边界线没有完全重合。此时,可以选择点施加边界条件,如下图。
    (2)导入的ABAQUS网格有两个Part。目前FssiCAS不支持两种材料之间设置接触、滑移功能(该功能正在开发中),不同材料之间要共节点,采用一套网格即一个Part。

    关于FssiCAS使用和PZIII模型不收敛的问题,有几点注意事项:
    (1)充分利用边界条件可视化功能,检查边界条件是否合理;
    (2)采用弹性模型一般只迭代一次,收敛性很好,如果出现上述弹性模型不收敛的情况,一定是模型设置有问题;
    (3)在多Step计算中,检查初始应力状态是否正确,这里推荐用后处理界面先查看第一步的结果,如果没问题再进行下一步, PZIII模型不收敛一个重要原因是第一个Step获取的初始应力不对。

      KunpengHe
      (1)使用Abaqus只能导入一个part吗?我是先画了两个矩形,合并成为一个整体的。那多层海床应该如何设置part呢,其他复杂的模型如何设置呢?
      (2)使用Case2的设置,只修改了沉积物参数,计算一定时间后仍然出现了“Stop Due To Negative Hmodu+Cons1 In Pastor_Zienkiewicz_Mark_III”这一问题是什么原因呢?基于Case2的设置其他都不变,应该不存在模型设置的问题。

        KunpengHe 老师,直接用case2的模型计算,计算到206s左右时,也是会出现“Stop Due To Negative Hmodu+Cons1 In Pastor_Zienkiewicz_Mark_III”这个问题。这个问题是说的什么呢,或者说是什么原因导致的

        ymbb 可能的原因是200s时候海床已经充分液化了,液化不可能永远算下去!

        • ymbb 回复了此帖

          叶剑红教授 老师咱现在这个软件就是到了充分液化就会停止吗,那我知道是啥问题了。因为在文献里看的用FSSI模型计算的时候,在海床表面充分液化以后还会继续计算一下。

          Du 多个part组装后,在assembly中进行网格共节点操作

          说点什么吧...