模型为四层结构，各层之间粘接连接起来，在顶层SiN薄膜中施加沿半径方向的力，观察各层薄膜滑移与相对弯曲情况。

建立模型为四分之一圆模型，通过接触单元contact174连接，运用TB，，，BILI命令定义接触区域材料属性。由于上面三层薄膜厚度远小于衬底厚度，为了简化模型 定义顶层三层薄膜全部为一体  只有衬底和二氧化硅层之间会滑移。四分之一模型关于截面对称。计算使整个四分之一模型的中轴线保持不动。接触单元的选项和实常数设置为：

!*  

KEYOPT,3,1,0

KEYOPT,3,2,0

KEYOPT,3,4,0

KEYOPT,3,5,2

KEYOPT,3,6,0

KEYOPT,3,7,1

KEYOPT,3,8,0

KEYOPT,3,9,0

KEYOPT,3,10,0   

KEYOPT,3,11,0   

KEYOPT,3,12,2   

KEYOPT,3,14,0   

KEYOPT,3,15,0   

KEYOPT,3,16,0   

!*  

R,3,0,0,1, 0.01,0,0

RMORE,0,0,0, 0,0,0  

RMORE,0,,, ,,   

RMORE,,,0, 0,,  

RMORE,,,,,, 

RMORE,,,,, 

计算的命令流设置如下：

TIME,0.1

AUTOTS,1

NSUBST,1000,1e6,10,1

KBC,0   

!*  

TSRES,ERASE 

!*  

NLGEOM,1

NROPT,FULL, ,OFF

STAOPT,DEFA 

LUMPM,0 

EQSLV, , ,0, ,DELE  

MSAVE,0 

PCGOPT,0, ,AUTO, , ,AUTO

PIVCHECK,1  

PSTRESS,0   

TOFFST,0,   

!*  

NEQIT,25,   

!* 

PRED,OFF 

LNSRCH,1

SAVE

最后爆出来错误如下The L-2 norm of the residual force overflowed. This may be caused by PRED,ON or birth/dead elements. Please use PRED,OFF and try again。

求大神解释