子模型法实例.docx

子模型法实例 finish/clear/FILNAM,coarse /PREP7 smrt,off!关闭智能化网格划分 /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的分段数 LARC,4,5,6,0.5 LESIZE,4,6 ! 定义线4的分段数 L,5,1 LESIZE,5,4,4 ! 定义线5的分段数 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!对节点数据排序，按x应力分量排序，由大到小，去前3 PRNSOL,S,COMP!列出节点结果，所有x,y,z,xy,yz,zx应力分量 /WINDOW,1,OFF!关闭先前定义的窗口，不能再对此窗口操作 /NOERASE!当前显示的内容被保留，随后显示的内容在其上连续叠加显示 /DSCALE,2,1!不对窗口2的位移进行缩放 /WINDOW,2,RTOP!1/4窗口，右上，也即定义下面用来显示内容的新的窗口 PLNSOL,S,X!图形显示节点x应力分量 *GET,CRSESTR,NODE,18,S,X!获取18号节点的x向应力分量 *STATUS!列表当前的所有缩写和参数 *DIM,LABEL,CHAR,1,2!定义CHAR数组（字符型数组） *DIM,VALUE,1,3!定义ARRAY数组（数值型数组） LABEL(1,1) = MX STR!赋予char数组值 LABEL(1,2) = CRS MODEL *VFILL,VALUE(1,1),DATA,3018!赋予array数组值 *VFILL,VALUE(1,2),DATA,CRSESTR *VFILL,VALUE(1,3),DATA,ABS(CRSESTR/3018) SAVE,TABLE_1!保存数据 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 ! 存储子模型 数据库文件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