子模型法實例.docx

子模型法實例 finish/clear/FILNAM,coarse /PREP7 smrt,off!關閉智能化網(wǎng)格劃分 /TITLE, STRESS CONCENTRATION AT A HOLE IN A PLATE /NOPR!下面的操作不被寫入LOG文件 ANTYPE,STATIC ! 靜力分析 ET,1,PLANE2 MP,EX,1,30E6 MP,NUXY,1,0.3 K,1,6 ! 關鍵點 K,2,6,6 K,3,6 K,4,.5 K,5,.5 K,6 L,1,2 L,2,3 L,3,4 LESIZE,3,4,.25 ! 定義線3的分段數(shù) LARC,4,5,6,0.5 LESIZE,4,6 ! 定義線4的分段數(shù) L,5,1 LESIZE,5,4,4 ! 定義線5的分段數(shù) AL,1,2,3,4,5 ESIZE,4 ! 每條線分為4段 AMESH,ALL /AUTO,1!設置自動合適的觀察方式 /PLOPTS,INFO,0!控制圖例的顯示 /PLOPTS,WINS,0! /WINDOW,LTOP!1/4窗口，左上 LSEL,S,LINE,3,5,2 DL,ALL,SYMM!對稱約束 LSEL,S,LINE,1 NSLL,1 SF,ALL,PRES,-1000. ! 施加拉力 LSEL,ALL NSEL,ALL CSYS,1 FINISH /SOLU SOLVE FINISH SAVE ! 存儲文件VM142.DB /POST1 SET,1,1 NSORT,S,X,3!對節(jié)點數(shù)據(jù)排序，按x應力分量排序，由大到小，去前3 PRNSOL,S,COMP!列出節(jié)點結果，所有x,y,z,xy,yz,zx應力分量 /WINDOW,1,OFF!關閉先前定義的窗口，不能再對此窗口操作 /NOERASE!當前顯示的內容被保留，隨后顯示的內容在其上連續(xù)疊加顯示 /DSCALE,2,1!不對窗口2的位移進行縮放 /WINDOW,2,RTOP!1/4窗口，右上，也即定義下面用來顯示內容的新的窗口 PLNSOL,S,X!圖形顯示節(jié)點x應力分量 *GET,CRSESTR,NODE,18,S,X!獲取18號節(jié)點的x向應力分量 *STATUS!列表當前的所有縮寫和參數(shù) *DIM,LABEL,CHAR,1,2!定義CHAR數(shù)組（字符型數(shù)組） *DIM,VALUE,1,3!定義ARRAY數(shù)組（數(shù)值型數(shù)組） LABEL(1,1) = MX STR!賦予char數(shù)組值 LABEL(1,2) = CRS MODEL *VFILL,VALUE(1,1),DATA,3018!賦予array數(shù)組值 *VFILL,VALUE(1,2),DATA,CRSESTR *VFILL,VALUE(1,3),DATA,ABS(CRSESTR/3018) SAVE,TABLE_1!保存數(shù)據(jù) FINISH /CLEAR,NOSTART !CLEAR THE DATABASE /FILNAM,SUBMODEL !定義子模型文件名 /PREP7 smrt,off /NOPR /TITLE, STRESS CONCENTRATION AT A HOLE IN A PLATE ANTYPE,STATIC ET,1,PLANE42 MP,EX,1,30E6 MP,NUXY,1,0.3 CSYS,1 K,10,.5,45 K,11,.5,90 K,12,1.5,45 K,13,1.5,90 A,10,12,13,11 ESIZE,8 MSHK,1 MSHA,0,2D AMESH,1 /WINDOW,2,OFF /NOERASE /PLOPTS,INFO,0 /PLOPTS,WINS,0 /WINDOW,1,LTOP /USER ! 用戶定義的比例 /DIST,1,3.3!設置縮放比例 /FOCUS,1,3,3 !設置縮放焦點 EPLOT LSEL,S,LINE,1,2 NSLL,1 ! 選擇切割邊界的結點 NWRITE ! 將結點寫入文件SUBMODEL.NODE LSEL,ALL NSEL,ALL FINISH SAVE ! 存儲子模型 數(shù)據(jù)庫文件SUBMODEL.DB /POST1 RESUME,coarse,db FILE,coarse,rst CBDOF,0,0 !激活邊界條件插值 FINISH /PREP7 smrt,off RESUME !從文件SUBMODEL.DB中讀入子模型 /NOPR /INPUT,cbdo,:cb1 ! 從文件SUBMODEL.CBDO中讀入插值邊界條件 /GOPR LSEL,S,LINE,3 ! 施加其余的邊界條件 DL,ALL,SYMM FINISH /SOLU SOLVE FINISH /POST1 SET,1,1 NSORT,S,X,3 PRNSOL,S,COMP /WINDOW,1,OFF /AUTO,3 /WINDOW,3,BOT /NOERASE /PLOPTS,MINM,1 /USER,3 /DIST,3,.2 /FOCUS,3,.2,.5 /CONTOUR,3,AUTO PLNSOL,SX *GET,SUBSTR,NODE,18,S,X *DIM,LABEL,CHAR,1,2 *DIM,VALUE,1,3 LABEL(1,1) = MAX STRS LABEL(1,2) = SUBMOD *VFILL,VALUE(1,1),DATA,3018 *VFILL,VALUE(1,2),DATA,SUBSTR *VFILL,VALUE(1,3),DATA,ABS(SUBSTR/3018 ) SAVE,TABLE_2 FINISH /CLEAR,NOSTART /FILNAM,SUBMODEL /PREP7 smrt,off /TITLE, STRESS CONCENTRATION AT A HOLE IN A PLATE /NOPR ANTYPE,STATIC ET,1,PLANE146 MP,EX,1,30E6 MP,NUXY,1,0.3 K,1,6 ! 關鍵點 K,2,6,6 K,3,6 K,4,.5 K,5,.5 K,6 L,1,2 L,2,3 L,3,4 LESIZE,3,4,.25 LARC,4,5,6,0.5 LESIZE,4,6 L,5,1 LESIZE,5,4,4 AL,1,2,3,4,5 ESIZE,4 AMESH,ALL /AUTO,1 /PLOPTS,INFO,0 /PLOPTS,WINS,0 /WINDOW,LTOP LSEL,S,LINE,3,5,2 DL,ALL,SYMM LSEL,S,LINE,1 NSLL,1 SF,ALL,PRES,-1000. LSEL,ALL NSEL,ALL CSYS,1 FINISH /SOLU SOLVE FINISH /POST1 SET,1,1 PRNSOL,S,COMP /WINDOW,1,OFF /NOERASE /DSCALE,2,1 /WINDOW,2,RTOP PLNSOL,S,X *GET,SUBSTR,NODE,18,S,X *DIM,LABEL,CHAR,1,2 *DIM,VALUE,1,3 LABE