如圖1和圖2是模型圖，淺綠色柱狀體為鐵鏈，深綠色為水池底部，對鐵鏈和水池底部做接觸分析，施加z軸負方向的重力載荷和y軸負方向的水流載荷。只對如圖鐵鏈最左端施加約束。鐵鏈只有節點，應該是點面接觸。

然而在求解過程中發現Absolute convergence norm會突然增大，如圖3所示。圖4為錯誤文件。我在未添加接觸對時并沒有發現這些問題，所以我相信是在添加了接觸對時出現此類問題。請問如何解決？

接觸對的命令流如下所示。

!******接觸對*******

CM,_NODECM,NODE

CM,_ELEMCM,ELEM

CM,_KPCM,KP

CM,_LINECM,LINE

CM,_AREACM,AREA

CM,_VOLUCM,VOLU

/GSAV,cwz,gsav,,temp

MP,MU,10,0.2

MAT,10  

R,3

REAL,3 

ET,4,170

ET,5,175

R,3,,,0.1,0.1,0,

NROPT,UNSYM

! Generate the target surface  目標面

ASEL,S,,,2 

CM,_TARGET,AREA

TYPE,4 

NSLA,S,1

ESLN,S,0

ESLL,U 

ESEL,U,ENAME,,188,189  

NSLE,A,CT2 

ESURF  

CMSEL,S,_ELEMCM

! Generate the contact surface 接觸面

NSEL,S,,,3559,3648,1  

CM,_CONTACT,NODE

TYPE,5 

ESLN,S,0

ESURF  

*SET,_REALID,3 

ALLSEL 

ESEL,ALL

ESEL,S,TYPE,,4 

ESEL,A,TYPE,,5 

ESEL,R,REAL,,3 

/PSYMB,ESYS,1  

/PNUM,TYPE,1

/NUM,1 

EPLOT  

ESEL,ALL

ESEL,S,TYPE,,4 

ESEL,A,TYPE,,5 

ESEL,R,REAL,,3 

CMSEL,A,_NODECM

CMDEL,_NODECM  

CMSEL,A,_ELEMCM

CMDEL,_ELEMCM  

CMSEL,S,_KPCM  

CMDEL,_KPCM

CMSEL,S,_LINECM

CMDEL,_LINECM  

CMSEL,S,_AREACM

CMDEL,_AREACM  

CMSEL,S,_VOLUCM

CMDEL,_VOLUCM  

/GRES,cwz,gsav 

CMDEL,_TARGET  

CMDEL,_CONTACT