模型為四層結構，各層之間粘接連接起來，在頂層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。

求大神解釋