Ansys 关于地震反应谱SRSS分析 APDL源码
我在ANSYS中作地震分解反应谱分析,一次X方向,一次Y方向,他们要求是独立互不干扰的,可是采用直进行一次模态分析的话,他生成的*.mcom文件好像是包含了前面的计算结果,命令流如下:
!进入PREP7并建模
/PREP7
B=15 !基本尺寸
A1=1000 !第一个面积
A2=1000 !第二个面积
A3=1000 !第三个面积
ET,1,beam4 !二维杆单元
R,1,0.25,0.0052,0.0052,0.5,0.5 !以参数形式的实参
MP,EX,1,2.0E11 !杨氏模量
mp,PRXY,1,,0.3
mp,dens,1,7.8e3
N,1,-B,0,0 !定义结点
N,2,0,0,0
N,3,-B,0,b
N,4,0,0,b
N,5,-B,0,2*b
N,6,0,0,2*b
N,7,-B,0,3*b
N,8,0,0,3*b
E,1,3 !定义单元
E,2,4
E,3,5
E,4,6
E,3,4
E,5,6
e,5,7
e,6,8
e,7,8
D,1,ALL,0,,2
FINISH
!
!进入求解器,定义载荷和求解
/SOLU
D,1,ALL,0,,2 !结点UX=UY=0
sfbeam,1,1,PRES,100000,
sfbeam,3,1,PRES,100000,
sfbeam,7,1,PRES,100000,
SOLVE
FINISH
allsel
NMODE=10
/SOL
!*
ANTYPE,2
!*
MSAVE,0
!*
MODOPT,LANB,NMODE
EQSLV,SPAR
MXPAND,NMODE , , ,1
LUMPM,0
PSTRES,0
!*
MODOPT,LANB,NMODE ,0,0, ,OFF
SOLVE
*DIM,FRE,,NMODE
*DO,I,1,NMODE
*GET,FRE(I),MODE,I,FREQ ! OBTAIN MODE FREQENCY FOR MODE I
*ENDDO
FINISH
!地震影响系数
grav=9.81
tg=0.35
amax=0.08
c=0.05
!
*dim,a,,nmode
*dim,t,,nmode
*do,i,1,nmode
t(i)=1.0/fre(i)
*enddo
r=0.9+(0.05-c)/(0.5+5.0*c)
p1=0.02+(0.05-c)/8
p2=1+(0.05-c)/(0.06+1.7*c)
全部命令请查看附件。