求助——简单的2D对称问题求解不收敛

fruitlover

我建立了一个2D模型，但贴不了图，只好简单描述一下：该模型关于y轴对称，在对称轴的最上方节点施加一个集中力，底边受全部自由度约束。但是为什么求解不收敛呢？材料属性：杨氏模量60E3，泊松比0.45.
代码如下：
/PREP7  
ET,1,PLANE182   
KEYOPT,1,1,0
KEYOPT,1,3,1
KEYOPT,1,6,0
KEYOPT,1,10,0   
KEYOPT,1,1,0
KEYOPT,1,3,1
KEYOPT,1,6,0
KEYOPT,1,10,0   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,60E3   
MPDATA,PRXY,1,,0.45
RECTNG,0,0.0225,0,-0.0045,  
AATT,       1, ,   1,       0,  
FLST,2,1,5,ORDE,1   
FITEM,2,1   
AESIZE,P51X,0.0005,
MSHAPE,0,2D
MSHKEY,1
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
AMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
FINISH  
/SOL
ANTYPE,0
NLGEOM,1
NSUBST,200,1000,10  
OUTRES,ERASE
OUTRES,ALL,ALL  
AUTOTS,1
TIME,100
DL,       4, ,SYMM  
FLST,2,1,4,ORDE,1   
FITEM,2,1   
/GO
DL,P51X, ,ALL,0
FLST,2,1,1,ORDE,1   
FITEM,2,56  
/GO
F,P51X,FY,-0.512
SOLVE   

我又建了该模型的全部模型，但求解可以收敛。代码如下：
/PREP7  
ET,1,PLANE182   
KEYOPT,1,1,0
KEYOPT,1,3,0
KEYOPT,1,6,0
KEYOPT,1,10,0   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,60E3   
MPDATA,PRXY,1,,0.45
RECTNG,-0.0225,0.0225,0,-0.0045,
AATT,       1, ,   1,       0,  
FLST,2,1,5,ORDE,1   
FITEM,2,1   
AESIZE,P51X,0.0005,
MSHAPE,0,2D
MSHKEY,1
CM,_Y,AREA  
ASEL, , , ,       1
CM,_Y1,AREA
CHKMSH,'AREA'   
CMSEL,S,_Y  
AMESH,_Y1   
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  
FINISH  
/SOL
ANTYPE,0
NLGEOM,1
NSUBST,200,1000,10  
OUTRES,ERASE
OUTRES,ALL,ALL  
AUTOTS,1
TIME,100
NSEL,S,LOC,X,0  
FLST,2,1,1,ORDE,1   
FITEM,2,146
/GO
F,P51X,FY,-1.024
ALLSEL,ALL  
FLST,2,1,4,ORDE,1   
FITEM,2,1   
/GO
DL,P51X, ,ALL,0
SOLVE   

请各位帮忙解答好么？这么简单的模型，我自己真的不知如何分析问题了。

fruitlover

已有人帮我解答了，KEYOPT,1,3,1改为KEYOPT,1,3,0就可以求解并收敛。
这里的轴对称是指旋转体的那种轴对称，而非关于坐标轴的平面对称。