本帖最后由 adidas04 于 2013-1-15 15:55 编辑
子模型分析实例(命令行格式) 问题描述 求解矩形平板中心开孔,承受横向拉力时的应力集中情况。材料特性和模型的几何形状见下图。 /FILNAM,coarse /PREP7 smrt,off /TITLE, STRESS CONCENTRATION AT A HOLE IN A PLATE /NOPR 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 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 PRNSOL,S,COMP /WINDOW,1,OFF /NOERASE /DSCALE,2,1 /WINDOW,2,RTOP PLNSOL,S,X *GET,CRSESTR,NODE,18,S,X *STATUS *DIM,LABEL,CHAR,1,2 *DIM,VALUE,,1,3 LABEL(1,1) = 'MX STR ' LABEL(1,2) = 'CRS MODEL' *VFILL,VALUE(1,1),DATA,3018 *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 LABEL(1,1) = 'MAX STRS' LABEL(1,2) = 'CRS MOD ' *VFILL,VALUE(1,1),DATA,3018 *VFILL,VALUE(1,2),DATA,SUBSTR *VFILL,VALUE(1,3),DATA,ABS(SUBSTR/3018 ) SAVE,TABLE_3 FINISH |