命令流代码

/PREP7
BLOCK,0,0.15,0,2.85,0,-0.15,
BLOCK,0.15,3,0,2.85,0,-0.15,
BLOCK,3,6,0,2.85,0,-0.15
BLOCK,6,9,0,2.85,0,-0.15
BLOCK,0,0.15,0,2.85,-0.15,-4.5,

BLOCK,0,0.15,2.85,3,0,-0.15,
BLOCK,0.15,3,2.85,3,0,-0.15
BLOCK,3,6,2.85,3,0,-0.15
BLOCK,6,9,2.85,3,0,-0.15
BLOCK,0,0.15,2.85,3,-0.15,-4.5,
BLOCK,0.15,3,2.85,3,-0.15,-4.5
BLOCK,3,6,2.85,3,-0.15,-4.5
BLOCK,6,9,2.85,3,-0.15,-4.5

BLOCK,0,0.15,3,6,0,-0.15,
BLOCK,0.15,3,3,6,0,-0.15,
BLOCK,3,6,3,6,0,-0.15,
BLOCK,6,9,3,6,0,-0.15,
BLOCK,0,0.15,3,6,-0.15,-4.5,

VSEL,S,,,ALL
VGLUE,ALL

/PNUM,AREA,ON
APLOT

VSEL,S,,,ALL
CM,F_SUB,VOLUM !前面部分
CMSEL,S,F_SUB
VSYMM,Z,ALL,,,,0,0

VSEL,S,,,ALL
VSEL,U,,,F_SUB
VGEN,,ALL,,,,,-9,,,1

VSEL,S,,,ALL
VGLUE,ALL

!定义单元
ET,1,SOLID45
MP,EX,1,2.5E8
MP,NUXY,1,0.2
MP,DENS,1,2500

LSEL,S,,,ALL
LSEL,S,LENGTH,,,0.15
LESIZE,ALL,,3
LSEL,S,,,ALL
LSEL,U,LENGTH,,0.15
LESIZE,ALL,0.15
MSHAPE,0,3D
MSHKEY,1
VSEL,S,,,ALL
VMESH,ALL

ASEL,S,LOC,Y,0
CM,BASE_F,AREA !基础面

ASEL,S,,,ALL

ASEL,S,LOC,X,9
CM,L_WALL,AREA !X对称面

ALLSEL
FINISH
/SOLU
ANTYPE,MODAL !求振型和频率
MODOPT,SUBSP,4 !用子空间迭代法求前3阶振型

ASEL,S,,,ALL
CMSEL,S,BASE_F
NSLA,S,1
D,ALL,ALL

ASEL,S,,,ALL
ASEL,S,,,L_WALL
NSLA,S
DSYM,SYMM,X

ALLSEL
MXPAND,5
SOLVE
FINISH

/POST1
SET,LIST,2
SET,FIRST
/VIEW,1,1,1,1,
/ANG,1
/REP,FAST
PLDISP,0
ANMODE,10,0.5E-1,20
/ANFILE,SAVE,"'1'',''AVI'',''''
SET,NEXT
/VIEW,1,1,1,1,
/ANG,1
/REP,FAST
PLDISP,0
ANMODE,10,0.5E-1,20
/ANFILE,SAVE,''2'',''AVI'',''''
SET,NEXT
/VIEW,1,1,1,1,
/ANG,1
/REP,FAST
PLDISP,0
ANMODE,10,0.5E-1,20
/ANFILE,SAVE,''3'',''AVI'',''''
/SOLU
ANTYPE,SPECTR !求振型和频率
SPOPT,SPRS,3,YES !单点谱分析
SVTYPE,2 !加速度谱
SED,1,0,0 !X方向地震作用
!定义反应谱(抗震规范GBJ11-89)：TG=0.4S,AMAX=0.16(8度) !F=A*G=A*M*G A值乘以9.82
FREQ,0.333333333,0.353773168,0.376883475,0.403224181,0.433523506,0.468746338,0.510198876,0.559694183,0.619824466,
SV,,0.256257352,0.270357153,0.286201463,0.304142684,0.324636126,0.348280823,0.375881551,0.408546201,0.447844563,
FREQ,0.694429977,0.789452909,0.914603474,1.086909265,1.339207992,1.744044089,2.5,10,1000
SV,,0.496079267,0.556774157,0.635616298,0.742436687,0.895877629,1.136284746,1.5712,1.5712,0.70704

SOLVE
FINISH
/SOLU
ANTYPE,MODAL !求振型和频率
EXPASS,ON
MXPAND,3,,,YES !扩展前3阶振型，计算单元解和反力
SOLVE
FINISH
/SOLU
ANTYPE,SPECTR
SRSS,1,DISP !应用SRSS方法计算单元应力,该命令将写*.MCOM文件
SOLVE
FINISH
/POST1
SET,LIST
/INP,,MCOM
PLDISP
PLETAB,AVG
PRNSOL,DOF
PRESOL,ELEM
PRRSOL,F