本帖最后由 c_shn 于 2011-7-14 16:49 编辑
没人回答,是我的问题说得不是很清楚吗?现在我把固体单元换成平面单元,发现结果依然一样。现附上我的命令流,希望各位大大看看哪里出错了。还是用大变形分析算法解决小变形问题的确会出错?我发现x方向的正应力是一致的,但剪切力真的很怪异。
FINISH
/CLEAR,START
/TITLE,`ENF`
/PREP7
RECTNG,0,16,-1.7,1.7
RECTNG,16,50,-1.7,1.7
RECTNG,50,100,-1.7,1.7
AGLUE,ALL
NUMCMP,ALL
ET,1,PLANE182
MP,EX,1,128000
MP,EY,1,9450
MP,PRXY,1,0.358
MP,GXY,1,4400
REAL,1
TYPE,1
MAT,1
LSEL,S,LOC,X,0
LSEL,A,LOC,X,16
LSEL,A,LOC,X,50
LESIZE,ALL,,,8
ESIZE,,10
MSHKEY,1
AMESH,1
AMESH,3
ESIZE,0.5
MSHKEY,1
AMESH,2
ALLSEL
FINISH
/CONFIG,NRES,900000
/SOLU
ANTYPE,STATIC
SOLCONTROL,1,1
NLGEOM,1 !选择大变形还是小变形分析
EQSLV,SPARSE
OUTRES,ALL,LAST
ALLSEL,ALL
CSYS,0
NSEL,S,LOC,X,100
NSEL,R,LOC,Y,-1.7
D,ALL,UY,0
NSEL,S,LOC,X,0
NSEL,R,LOC,Y,-1.7
D,ALL,UY,0
!D,ALL,UY,4
D,ALL,UX,0
NSEL,S,LOC,X,50
NSEL,R,LOC,Y,1.7
D,ALL,UY,-4
!D,ALL,UY,0
!D,ALL,UX,0
ALLSEL,ALL
OUTRES,ALL,ALL
NSUBST,50,100,20
SOLVE
FINISH
/POST1 |