153﹐154例子 /PREP7 /TITLE, VM38, PLASTIC LOADING OF A THICK-WALLED CYLINDER UNDER PRESSURE ET,1,PLANE42,,1,1 ET,2,SURF153,,,1,1 ! 无中节点的轴对称2D表面效应单元 MP,EX,1,30E6 ﹗材料特性 MP,NUXY,1,.3 TB,BKIN,1,1 ! 双线性随动硬化 TBTEMP,70 TBDATA,1,30000,0 !屈服应力和0正切模数 N,1,4 ! 定义节点 N,6,8 FILL NGEN,2,10,1,6,1,,1 E,11,1,2,12 ! 定义单元 EGEN,5,1,1 CPNGEN,1,UY,11,16 ! 耦合节点 TYPE,2 ! 创建SURF153单元﹐施加表面压力 NSEL,S,LOC,X,4 ESURF NSEL,ALL TREF,70 ! 边界条件及载荷 D,1,UY,,,6 FINISH /SOLU ESEL,S,TYPE,,2 ! 选择SURF153单元施加表面压力 !塑性分析 SFE,ALL,1,PRES,,12990 ESEL,ALL OUTPR,BASIC,1 SOLVE FINISH /POST1 ETABLE,STRS_R,S,X ETABLE,STRS_T,S,Z *GET,SIGR_I,ELEM,1,ETAB,STRS_R *GET,SIGT_I,ELEM,1,ETAB,STRS_T *GET,SIGR_O,ELEM,5,ETAB,STRS_R *GET,SIGT_O,ELEM,5,ETAB,STRS_T *DIM,LABEL,CHAR,4,2 *DIM,VALUE,,4,3 LABEL(1,1) = 'SIGR LFT','SIGT LFT','SIGR RT_','SIGT RT_' LABEL(1,2) = '_END psi','_END psi','END psi ','END psi ' *VFILL,VALUE(1,1),DATA,-9984,18645,-468,9128 *VFILL,VALUE(1,2),DATA,SIGR_I,SIGT_I,SIGR_O,SIGT_O *VFILL,VALUE(1,3),DATA,ABS(SIGR_I/9984),ABS(SIGT_I/18645) *VFILL,VALUE(3,3),DATA,ABS(SIGR_O/468),ABS(SIGT_O/9128) SAVE,TABLE_1 FINISH /SOLU ESEL,S,TYPE,,2 ! 选择SURF153单元施加表面压力 ! 塑性分析 SFE,ALL,1,PRES,,24012 ESEL,ALL SOLVE FINISH /PREP7 EDELE,ALL NDELE,ALL ! 移除模型 ET,1,SOLID45,1 ET,2,SURF154,,,,1 ! 无中节点的轴对称3D表面效应单元 CSYS,1 N,1,4,-2.5 ! 定义节点 N,6,8,-2.5 FILL NGEN,2,6,1,6,1,,5 NGEN,2,12,1,12,1,,,1 NUMCMP,ELEM NUMSTR,ELEM,1 TYPE,1 MAT,1 E,1,2,8,7,13,14,20,19 ! 定义单元 EGEN,5,1,-1 TYPE,2 ! 创建SURF153单元﹐施加表面压力 NSEL,S,NODE,,1,7,6 NSEL,A,NODE,,13,19,6 ESURF NSEL,ALL NROTAT,ALL ! 旋转所有节点到圆柱坐标系下 CPDELE,1,1,1 ! 移除节点耦合 SFDELE,ALL,PRES ! 移除节点压力 D,ALL,UY,0.0 ! 强迫所有节点在切线方向 NSEL,S,LOC,Z,1 ! 选择Z方向为1的节点 CP,1,UZ,ALL ! 在UZ方向耦合选择的节点来模拟普通的3D平面应变形为 NSEL,S,LOC,Z,0 ! 在UZ方向耦合Z坐标等於0的节点 D,ALL,UZ NSEL,ALL FINISH /SOLU ESEL,S,TYPE,,2 ! 选择SURF154单元施加表面压力 ! 塑性分析载荷 SFE,ALL,1,PRES,,12990 ESEL,ALL OUTPR,BASIC,1 SOLVE finish /POST1 ETABLE,STRS_R,S,X ETABLE,STRS_T,S,Y *GET,SIGR_I,ELEM,1,ETAB,STRS_R *GET,SIGT_I,ELEM,1,ETAB,STRS_T *GET,SIGR_O,ELEM,5,ETAB,STRS_R *GET,SIGT_O,ELEM,5,ETAB,STRS_T LABEL(1,1) = 'SIGR LFT','SIGT LFT','SIGR RT_','SIGT RT_' LABEL(1,2) = '_END psi','_END psi','END psi ','END psi ' *VFILL,VALUE(1,1),DATA,-9984,18645,-468,9128 *VFILL,VALUE(1,2),DATA,SIGR_I,SIGT_I,SIGR_O,SIGT_O *VFILL,VALUE(1,3),DATA,ABS(SIGR_I/9984),ABS(SIGT_I/18645) *VFILL,VALUE(3,3),DATA,ABS(SIGR_O/468),ABS(SIGT_O/9128) SAVE,TABLE_3 FINISH /solu ESEL,S,TYPE,,2 ! 选择SURF154单元施加表面压力PRESSURE ! 塑性分析载荷 SFE,ALL,1,PRES,,24012 ESEL,ALL SOLVE FINISH /POST1 RESUME,TABLE_1 /COM /OUT,vm38,vrt /COM,------------------- VM38 RESULTS COMPARISON --------------- /COM, /COM, | TARGET | ANSYS | RATIO /COM, /COM,FULLY ELASTIC, PLANE42 RESULTS: /COM, *VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3) (1X,A8,A8,' ',F10.0,' ',F10.0,' ',1F5.3) /NOPR RESUME,TABLE_3 /GOPR /COM, /COM,FULLY ELASTIC, PLANE45 RESULTS: /COM, *VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3) (1X,A8,A8,' ',F10.0,' ',F10.0,' ',1F5.3) /COM,----------------------------------------------------------- /COM, /COM,----------------------------------------------------------- /COM,NOTE: THERE ARE VERIFIED RESULTS IN VM38 NOT CONTAINED IN /COM,THIS TABLE /COM,----------------------------------------------------------- /OUT FINISH *LIST,vm38,vrt ※ 引述《yunsoon (不断思考)》之铭言： : 小弟有个Ansys元素种类问题想请教板上各位先进， : 就是有没有哪位先进有使用过SURF154这种元素种类， : Surface Effect---->3D Structure 154 : 这种元素在mesh时是否要以面mesh，而不能用体mesh？ : 因为小弟是使用这种元素种类进行整体mesh时， : 会出现下列警告：『Which cannot be used with the VEMS command』， : 请问这句话的含意是什麽，有没有解决的办法， : 小弟恳请各位先进给个提示，谢谢。 -- 酒 楚王帐里醉听歌﹐北戍迎风夜渡河 避乱草堂须典当﹐谪居赤壁借消磨 愁肠一缕红烛泪﹐壮志孤灯宝剑何 古往今来尽如此﹐无情伴我渐蹉跎 --

※ 发信站: 批踢踢实业坊(ptt.cc) ◆ From: 222.171.7.223