如下图，做了一个关于Cohesive单元的算例

图中蓝色与红色单元为线弹性本构模型，单元尺寸为1*1*0.1(mm)。中间有一层Cohesive单元，单元尺寸为1*1*0(mm),采用的本构模型为Mat138，即*MAT_COHESIVE_MIXED_MODE本构模型。在蓝色单元顶部四个节点施加竖直向上位移，红色底部四个节点添加竖直方向约束。

上图为Mat138的拉力位移曲线，其中设置L=0.05mm，u=0.1mm，GC=0.4mJ,QMAX=8N,E=160MPa。

结果如下图所示：

整个过程在0-0.05mm阶段，输出的能量、节点力等参数可以与理论值一一对应，但在0.05mm之后，也就是第二张图中力下降的阶段，模型发生了突变，如下图。

请教发生这种问题的原因。

Mat138材料参数：长度mm-质量t-时间s-模量Mpa-能量mJ-力N

分割线----------------------------------------------------------

以下没有变动单元大小，，看着有点怪...

失效前一个状态

失效时的状态

能量曲线

F-T曲线

顶点有些波动