谢谢venture百忙中的计算,我把节点坐标系改成沿着曲面的法向后,也得出合理的结果了,以前球体是向一侧有刚体位移的。
但是现在我把模型改成椭球体用同样的方法计算,又出现向一侧的刚体位移了。模型中椭球长轴2厘米,短轴1厘米,其它条件不变。命令流如下。
/PREP7
ET,1,SOLID95
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2e11
MPDATA,PRXY,1,,0.3
CSWPLA,11,1,0.5,1,
K,1,-0.02,0,0,
K,2,0.02,0,0,
K,3,-0.019,0,0,
K,4,0.019,0,0,
L, 1, 2
L, 3, 4
LSTR, 1, 3
LSTR, 4, 2
FLST,2,4,4
FITEM,2,1
FITEM,2,2
FITEM,2,3
FITEM,2,4
AL,P51X
FLST,2,1,5,ORDE,1
FITEM,2,1
FLST,8,2,3
FITEM,8,1
FITEM,8,2
VROTAT,P51X, , , , , ,P51X, ,90, ,
SMRT,6
SMRT,3
MSHAPE,1,3D
MSHKEY,0
CM,_Y,VOLU
VSEL, , , , 1
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
VMESH,_Y1
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
/UI,MESH,OFF
FLST,2,2,5,ORDE,2
FITEM,2,2
FITEM,2,-3
NORA,P51X,1
CSYS,0,
FINISH
/SOL
FLST,2,2,5,ORDE,2
FITEM,2,1
FITEM,2,4
DA,P51X,SYMM
FLST,2,1,5,ORDE,1
FITEM,2,2
/GO
SFA,P51X,1,PRES,1e6
FLST,2,1,5,ORDE,1
FITEM,2,3
/GO
SFA,P51X,1,PRES,100000
CSYS,0
ALLSEL,ALL
SOLVE |