ANSYS地震反应谱SRSS分析
我在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)
*do,i,1,nmode
*if,t(i),ge,0.0,and,t(i),lt,0.1,then
a(i)=(0.45+(10.0*p2-4.5)*t(i))*amax*grav
*elseif,t(i),ge,0.1,and,t(i),le,tg
a(i)=p2*amax*grav
*elseif,t(i),gt,tg,and,t(i),le,5*tg
a(i)=(tg/t(i))**r*p2*amax*grav
*else
a(i)=(p2*0.2**r-p1*(t(i)-5*tg))*amax*grav
*endif
*enddo
!
! X-方向谱分析 Spectrum analysis along Global X-axis direction
/SOLU
ANTYPE,SPECTR ! Spectrum analysis
SPOPT,SPRS ! Single point spectrum
SED,1,, ! Global X-axis as spectrum direction
SVTYP,2 ! Seismic acceleration response spectrum
! Frequency points and Spectrum values for SV vs. freq. table
FREQ,fre(1),fre(2),fre(3),fre(4),fre(5),fre(6),fre(7),fre(8),fre(9)
FREQ,fre(10)
SV,,a(1),a(2),a(3),a(4),a(5),a(6),a(7),a(8),a(9)
SV,,a(10)
FINISH
!/SOLU
!ANTYPE,MODAL ! Mode-frequency analysis
!EXPASS,ON
!MXPAND,nmode,,,YES,0.0 ! Expand nmode shapes, calculate element stresses
!SOLVE
!FINISH
/SOLU
ANTYPE,SPECTR
SRSS,0.0,DISP ! Square Root of Sum of Squares Mode combination
! with signif=0.0 and displacement solution requested
SOLVE
FINISH
/POST1
SET,LIST
/INPUT,,mcom
!***************EARTHQUAKE X******************
ALLSEL,ALL
FINISH
! Y-方向谱分析 Spectrum analysis along Global X-axis direction
!!**********************************************
!/SOL
!!*
!ANTYPE,2
!!*
!MSAVE,0
!!*
!MODOPT,LANB,NMODE
!EQSLV,SPAR
!MXPAND,NMODE , , ,1
!LUMPM,0
!PSTRES,0
!!*
!MODOPT,LANB,NMODE ,0,0, ,OFF
!SOLVE
!FINISH
!!**********************************************
/SOLU
LSCLEAR,LSOPT
ANTYPE,SPECTR ! Spectrum analysis
SPOPT,SPRS ! Single point spectrum
SED,,1, ! Global Y-axis as spectrum direction
SVTYP,2 ! Seismic acceleration response spectrum
FREQ
! Frequency points and Spectrum values for SV vs. freq. table
FREQ,fre(1),fre(2),fre(3),fre(4),fre(5),fre(6),fre(7),fre(8),fre(9)
FREQ,fre(10)
SV,,a(1),a(2),a(3),a(4),a(5),a(6),a(7),a(8),a(9)
SV,,a(10)
SOLVE
FINISH
!/SOLU
!ANTYPE,MODAL ! Mode-frequency analysis
!EXPASS,ON
!MXPAND,nmode,,,YES,0.0 ! Expand nmode shapes, calculate element stresses
!SOLVE
!FINISH
/SOLU
ANTYPE,SPECTR
SRSS,0.0,DISP ! Square Root of Sum of Squares Mode combination
! with signif=0.0 and displacement solution requested
SOLVE
FINISH
/POST1
SET,LIST
/INPUT,,mcom
!***************EARTHQUAKE Y******************
ALLSEL,ALL
FINISH
这里在进行X方向的反应谱分析以后，进行Y方向的分析，可是他生成的*.mcom文件如下：
/COM,ANSYS RELEASE 8.0 UP20030930 09:28:42 07/23/2005
/COM, truss.mcom
LCOPER,ZERO
LCDEFI,1, 1, 1
LCFACT,1, 0.263825E-17
LCASE,1
LCOPER,SQUARE
LCDEFI,1, 1, 2
LCFACT,1, 8.55778
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 3
LCFACT,1, -0.188669E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 4
LCFACT,1, -0.871099E-15
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 5
LCFACT,1, -0.757013
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 6
LCFACT,1, 0.967307E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 7
LCFACT,1, 0.533141E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 8
LCFACT,1, -0.203699
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 9
LCFACT,1, 0.445795E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 10
LCFACT,1, -0.387808E-13
LCOPER,ADD,1,MULT,1
LCOPER,SQRT
/COM,ANSYS RELEASE 8.0 UP20030930 09:28:42 07/23/2005
/COM, truss.mcom
LCOPER,SQUARE !注意这里没有清空数据库
LCDEFI,1, 1, 1
LCFACT,1, 50.7528
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 2
LCFACT,1, 0.887017E-14
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 3
LCFACT,1, 0.612824E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 4
LCFACT,1, -1.96484
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 5
LCFACT,1, -0.331613E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 6
LCFACT,1, 0.330459E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 7
LCFACT,1, 0.366569
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 8
LCFACT,1, -0.976991E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 9
LCFACT,1, 0.417313E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 10
LCFACT,1, 0.401040E-13
LCOPER,ADD,1,MULT,1
LCOPER,SQRT
我感觉这样好像是X和Y两个方向地震的叠加，
可是如果在座Y方向的地震以前把注释掉的模态分析在做一下这样的Y方向的地震的*.mcom就是：
/COM, truss.mcom
LCOPER,ZERO !注意这里清空数据库
LCDEFI,1, 1, 1
LCFACT,1, 50.7528
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 2
LCFACT,1, 0.887017E-14
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 3
LCFACT,1, 0.612824E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 4
LCFACT,1, -1.96484
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 5
LCFACT,1, -0.331613E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 6
LCFACT,1, 0.330459E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 7
LCFACT,1, 0.366569
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 8
LCFACT,1, -0.976991E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 9
LCFACT,1, 0.417313E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 10
LCFACT,1, 0.401040E-13
LCOPER,ADD,1,MULT,1
LCOPER,SQRT
如果在X方向后不作Y方向的地震，他的*.mcom：
/COM,ANSYS RELEASE 8.0 UP20030930 08:46:23 07/23/2005
/COM, truss.mcom
LCOPER,ZERO
LCDEFI,1, 1, 1
LCFACT,1, 0.263825E-17
LCASE,1
LCOPER,SQUARE
LCDEFI,1, 1, 2
LCFACT,1, 8.55778
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 3
LCFACT,1, -0.188669E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 4
LCFACT,1, -0.871099E-15
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 5
LCFACT,1, -0.757013
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 6
LCFACT,1, 0.967307E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 7
LCFACT,1, 0.533141E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 8
LCFACT,1, -0.203699
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 9
LCFACT,1, 0.445795E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 10
LCFACT,1, -0.387808E-13
LCOPER,ADD,1,MULT,1
LCOPER,SQRT
可是在X后作Y他不清空数据库，需要进行两次模态分析，这很耗时间对于大型结构，请大家讨论讨论如何处理呢？
Re:讨论：ANSYS地震反应谱SRSS分析
本人是学土木工程的，平时主要用Patran+Nastran对结构做线性分析，偶尔使用Ansys对结构做地震反应谱分析，但对Ansys的命令流不熟悉。依据本人的使用经验看，“一次X方向，一次Y方向，他们要求是独立互不干扰的”应该这样完成才对：1、点击Solution-〉Load Step Opts-〉Spectrum-〉Settings，在对话框中将SEDX，SEDY，SEDZ设置为1，0，0后；点击Solution-〉Solve-〉Current LS；读入生成*.mcom文件完成X方向的反应谱内力计算。2、点击Solution-〉Load Step Opts-〉Spectrum-〉Settings，在对话框中将SEDX，SEDY，SEDZ设置为0，1，0后；点击Solution-〉Solve-〉Current LS；读入生成*.mcom文件完成Y方向的反应谱内力计算。需要注意的是，在对话框中将SEDX，SEDY，SEDZ为1，1，0表示地震动的方向为与X轴成45度的方向，这显然不是所希望的
Re:讨论：ANSYS地震反应谱SRSS分析
我也知道是这样，可是他形成的*.mcom好像没有不是这个意思啊。
按你说的\"，“一次X方向，一次Y方向，他们要求是独立互不干扰的”应该这样完成才对：1、点击Solution-〉Load Step Opts-〉Spectrum-〉Settings，在对话框中将SEDX，SEDY，SEDZ设置为1，0，0后；点击Solution-〉Solve-〉Current LS；读入生成*.mcom文件完成X方向的反应谱内力计算。2、点击Solution-〉Load Step Opts-〉Spectrum-〉Settings，在对话框中将SEDX，SEDY，SEDZ设置为0，1，0后；点击Solution-〉Solve-〉Current LS；读入生成*.mcom文件完成Y方向的反应谱内力计算。\"做后：他把Y方向的同样写入到同一个*.mcom中了，而且也没有晴空数据库，形成如下的*.mcom
/COM,ANSYS RELEASE 8.0 UP20030930 09:28:42 07/23/2005
/COM, truss.mcom ！X方向
LCOPER,ZERO
LCDEFI,1, 1, 1
LCFACT,1, 0.263825E-17
LCASE,1
LCOPER,SQUARE
LCDEFI,1, 1, 2
LCFACT,1, 8.55778
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 3
LCFACT,1, -0.188669E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 4
LCFACT,1, -0.871099E-15
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 5
LCFACT,1, -0.757013
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 6
LCFACT,1, 0.967307E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 7
LCFACT,1, 0.533141E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 8
LCFACT,1, -0.203699
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 9
LCFACT,1, 0.445795E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 10
LCFACT,1, -0.387808E-13
LCOPER,ADD,1,MULT,1
LCOPER,SQRT
/COM,ANSYS RELEASE 8.0 UP20030930 09:28:42 /COM, truss.mcom ！Y方向的
LCOPER,SQUARE ！这里没有用zero清空
LCDEFI,1, 1, 1
LCFACT,1, 50.7528
07/23/2005
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 2
LCFACT,1, 0.887017E-14
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 3
LCFACT,1, 0.612824E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 4
LCFACT,1, -1.96484
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 5
LCFACT,1, -0.331613E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 6
LCFACT,1, 0.330459E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 7
LCFACT,1, 0.366569
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 8
LCFACT,1, -0.976991E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 9
LCFACT,1, 0.417313E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 10
LCFACT,1, 0.401040E-13
LCOPER,ADD,1,MULT,1
LCOPER,SQRT
Re:讨论：ANSYS地震反应谱SRSS分析
本版精华贴中好像有关于地震分析的东东，各位可以找找看，希望有用。
Re:讨论：ANSYS地震反应谱SRSS分析
我也发现那个问题，我是利用 /copy 命令去解决的，模态分析完成后
下面命令流：（0_1Ydizh.txt与0_1Zdizh.txt为y、z方向的谱分析命令流。）
/COPY,file,db,,00,db
/COPY,file,full,,00,full
/COPY,file,esav,,00,esav
/COPY,file,mode,,00,mode
/COPY,file,rst,,00,rst
!
/COPY,00,db,,file,db
/COPY,00,full,,file,full
/COPY,00,esav,,file,esav
/COPY,00,mode,,file,mode
/COPY,00,rst,,file,rst
RESUME,file,db,,0,0
!
/INPUT,0_1Ydizh,txt,,, 0
!
LCWRITE,30,fY,L01,
/COPY,file,mcom,,fY,mcom
!
/CLEAR,START
/COPY,00,db,,file,db
/COPY,00,full,,file,full
/COPY,00,esav,,file,esav
/COPY,00,mode,,file,mode
/COPY,00,rst,,file,rst
!
RESUME,file,db,,0,0
/INPUT,0_1Zdizh,txt,,, 0
!
LCWRITE,31,fz,L01,
/COPY,file,mcom,,fz,mcom
!
/CLEAR,START
RESUME,file,db,,0,0
/SOLU
ANTYPE,0
SOLVE
FINISH
.......
Re:讨论：ANSYS地震反应谱SRSS分析
3ks
我后来也是这么实现的
因篇幅问题不能全部显示，请点此查看更多更全内容