如下圖，做了一個(gè)關(guān)于Cohesive單元的算例

圖中藍(lán)色與紅色單元為線彈性本構(gòu)模型，單元尺寸為1*1*0.1(mm)。中間有一層Cohesive單元，單元尺寸為1*1*0(mm),采用的本構(gòu)模型為Mat138，即*MAT_COHESIVE_MIXED_MODE本構(gòu)模型。在藍(lán)色單元頂部四個(gè)節(jié)點(diǎn)施加豎直向上位移，紅色底部四個(gè)節(jié)點(diǎn)添加豎直方向約束。

上圖為Mat138的拉力位移曲線，其中設(shè)置L=0.05mm，u=0.1mm，GC=0.4mJ,QMAX=8N,E=160MPa。

結(jié)果如下圖所示：

整個(gè)過程在0-0.05mm階段，輸出的能量、節(jié)點(diǎn)力等參數(shù)可以與理論值一一對(duì)應(yīng)，但在0.05mm之后，也就是第二張圖中力下降的階段，模型發(fā)生了突變，如下圖。

請(qǐng)教發(fā)生這種問題的原因。

Mat138材料參數(shù)：長(zhǎng)度mm-質(zhì)量t-時(shí)間s-模量Mpa-能量mJ-力N

分割線----------------------------------------------------------

以下沒有變動(dòng)單元大小，，看著有點(diǎn)怪...

失效前一個(gè)狀態(tài)

失效時(shí)的狀態(tài)

能量曲線

F-T曲線

頂點(diǎn)有些波動(dòng)