fruitlover 发表于 2010-4-13 08:31

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

我建立了一个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 发表于 2010-4-13 10:52

已有人帮我解答了,KEYOPT,1,3,1改为KEYOPT,1,3,0就可以求解并收敛。
这里的轴对称是指旋转体的那种轴对称,而非关于坐标轴的平面对称。
页: [1]
查看完整版本: 求助——简单的2D对称问题求解不收敛