Re: [問題] 請教Ansys元素種類(SURF154)問題

看板Cad_Cae (電腦輔助設計)作者wapi (青青子衿)時間19年前 (2006/05/09 09:10)推噓0(0推 0噓 0→)

留言0則, 0人參與討論串2/2 (看更多)

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