ANSYS地震反应谱SRSS分析共24页

1、ANSYS地震反应谱SRSS分析我在ANSYS中作地震分解反应谱分析，一次X方向，一次Y方向，他们要求是独立互不干扰的，可是采用直进行一次模态分析的话，他生成的*.mcom文件好像是包含了前面的计算结果，命令流如下：!进入PREP7并建模/PREP7B=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.3mp,dens,1,7.8e3N,1,-B,0,0 !定义结点N,

2、2,0,0,0N,3,-B,0,bN,4,0,0,bN,5,-B,0,2*bN,6,0,0,2*bN,7,-B,0,3*bN,8,0,0,3*bE,1,3 !定义单元E,2,4E,3,5E,4,6E,3,4E,5,6e,5,7e,6,8e,7,8D,1,ALL,0,2FINISH!进入求解器，定义载荷和求解/SOLUD,1,ALL,0,2 !结点 UX=UY=0sfbeam,1,1,PRES,100000,sfbeam,3,1,PRES,100000,sfbeam,7,1,PRES,100000,SOLVEFINISHallselNMODE=10/SOL!*ANTYPE,2!*MSAVE,0!

3、*MODOPT,LANB,NMODEEQSLV,SPARMXPAND,NMODE , , ,1LUMPM,0PSTRES,0!*MODOPT,LANB,NMODE ,0,0, ,OFFSOLVE*DIM,FRE,NMODE*DO,I,1,NMODE*GET,FRE(I),MODE,I,FREQ ! OBTAIN MODE FREQENCY FORMODE I*ENDDOFINISH!地震影响系数grav=9.81tg=0.35amax=0.08c=0.05!*dim,a,nmode*dim,t,nmode*do,i,1,nmodet(i)=1.0/fre(i)*enddor=0.9+(0.05

4、-c)/(0.5+5.0*c)p1=0.02+(0.05-c)/8p2=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,thena(i) = (0.45+(10.0*p2-4.5)*t(i)*amax*grav*elseif,t(i),ge,0.1,and,t(i),le,tga二p2*amax*grav*elseif,t(i),gt,tg,and,t(i),le,5*tga(i) = (tg/t(i)*r*p2*amax*grav*elsea(i) = (p2*0.2*r-p1*(t(i)-5*tg)

5、*amax*grav*endif*enddo! X-方向谱分析 Spectrum analysis along Global X-axisdirection/SOLUANTYPE,SPECTR ! Spectrum analysisSPOPT,SPRS ! Single point spectrumSED,1, ! Global X-axis as spectrum directionSVTYP,2 ! Seismic acceleration response spectrum! Frequency points and Spectrum values for SV vs. freq.tab

6、leFREQ,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/SOLUANTYPE,SPECTRSR

7、SS,0.0,DISP ! Square Root of Sum of Squares Modecombination! with signif=0.0 anddisplacement solution requestedSOLVEFINISH/POST1SET,LIST/INPUT,mcom!*EARTHQUAKE X*ALLSEL,ALLFINISH! Y-方向谱分析 Spectrum analysis along Global X-axisdirection!*!/SOL!*!ANTYPE,2!*!MSAVE,0!*!MODOPT,LANB,NMODE!EQSLV,SPAR!MXPAND

8、,NMODE , , ,1!LUMPM,0!PSTRES,0!*!MODOPT,LANB,NMODE ,0,0, ,OFF!SOLVE!FINISH!.!*/SOLULSCLEAR,LSOPTANTYPE,SPECTR ! Spectrum analysisSPOPT,SPRS ! Single point spectrumSED,1, ! Global Y-axis as spectrum directionSVTYP,2 ! Seismic acceleration response spectrumFREQ! Frequency points and Spectrum values fo

9、r SV vs. freq.tableFREQ,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)SOLVEFINISH!/SOLU!ANTYPE,MODAL ! Mode-frequency analysis!EXPASS,ON!MXPAND,nmode,YES,0.0 ! Expand nmode shapes,calculate element stresses!SOLVE!FINI

10、SH/SOLUANTYPE,SPECTRSRSS,0.0,DISP ! Square Root of Sum of Squares Modecombination! with signif=0.0 anddisplacement solution requestedSOLVEFINISH/POST1SET,LIST/INPUT,mcom!*EARTHQUAKE y*ALLSEL,ALLFINISH这里在进行X方向的反应谱分析以后，进行Y方向的分析，可是他生成的*.mcom文件如下：/COM,ANSYS RELEASE 8.0 UP20030930 09:28:42 07/23/2005/COM

11、, truss.mcomLCOPER,ZEROLCDEFI,1, 1, 1LCFACT,1, 0.263825E-17LCASE,1LCOPER,SQUARELCDEFI,1, 1, 2LCFACT,1, 8.55778LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 3LCFACT,1, -0.188669E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 4LCFACT,1, -0.871099E-15LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 5LCFACT,1, -0.757013LCOPER,ADD,1,MULT,1LCDEFI,

12、1, 1, 6LCFACT,1, 0.967307E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 7LCFACT,1, 0.533141E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 8LCFACT,1, -0.203699LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 9LCFACT,1, 0.445795E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 10LCFACT,1, -0.387808E-13LCOPER,ADD,1,MULT,1LCOPER,SQRT/COM,ANSYS RELEASE 8.0

13、 UP20030930 09:28:42 07/23/2005/COM, truss.mcomLCOPER,SQUARE !注意这里没有清空数据库LCDEFI,1, 1, 1LCFACT,1, 50.7528LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 2LCFACT,1, 0.887017E-14LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 3LCFACT,1, 0.612824E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 4LCFACT,1, -1.96484LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 5LC

14、FACT,1, -0.331613E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 6LCFACT,1, 0.330459E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 7LCFACT,1, 0.366569LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 8LCFACT,1, -0.976991E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 9LCFACT,1, 0.417313E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 10LCFACT,1, 0.401040E-13LCOPER

15、,ADD,1,MULT,1LCOPER,SQRT我感觉这样好像是X和Y两个方向地震的叠加，可是如果在座Y方向的地震以前把注释掉的模态分析在做一下这样的Y方向的地震的*.mcom就是：/COM, truss.mcomLCOPER,ZERO !注意这里清空数据库LCDEFI,1, 1, 1LCFACT,1, 50.7528LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 2LCFACT,1, 0.887017E-14LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 3LCFACT,1, 0.612824E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1

16、, 4LCFACT,1, -1.96484LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 5LCFACT,1, -0.331613E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 6LCFACT,1, 0.330459E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 7LCFACT,1, 0.366569LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 8LCFACT,1, -0.976991E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 9LCFACT,1, 0.417313E-13LCOPER

17、,ADD,1,MULT,1LCDEFI,1, 1, 10LCFACT,1, 0.401040E-13LCOPER,ADD,1,MULT,1LCOPER,SQRT如果在X方向后不作Y方向的地震，他的*.mcom：/COM,ANSYS RELEASE 8.0 UP20030930 08:46:23 07/23/2005/COM, truss.mcomLCOPER,ZEROLCDEFI,1, 1, 1LCFACT,1, 0.263825E-17LCASE,1LCOPER,SQUARELCDEFI,1, 1, 2LCFACT,1, 8.55778LCOPER,ADD,1,MULT,1LCDEFI,1,

18、 1, 3LCFACT,1, -0.188669E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 4LCFACT,1, -0.871099E-15LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 5LCFACT,1, -0.757013LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 6LCFACT,1, 0.967307E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 7LCFACT,1, 0.533141E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 8LCFACT,1, -0.203699LC

19、OPER,ADD,1,MULT,1LCDEFI,1, 1, 9LCFACT,1, 0.445795E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 10LCFACT,1, -0.387808E-13LCOPER,ADD,1,MULT,1LCOPER,SQRT可是在X后作Y他不清空数据库，需要进行两次模态分析,这很耗时间对于大型结构，请大家讨论讨论如何处理呢？Re:讨论：ANSYS地震反应谱SRSS分析 本人是学土木工程的，平时主要用Patran+Nastran对结构做线性分析， 偶尔使用Ansys对结构做地震反应谱分析，但对Ansys的命令流不熟悉。 依据本人的使用经验看

20、，“一次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

21、文件完成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文件完

22、成X方向的反 应谱内力计算。2、点击 Solution-) Load Step Opts-)Spectrum-) Settings， 在对话框中将 SEDX，SEDY，SEDZ 设置为 0，1，0 后；点击 Solution-)Solve-) Current LS；读入生成*.mcom文件完成Y方向的反应谱内力计算。做后：他把Y方向的同样写入到同一个*.mcom中了，而且也没有晴空数据库，形成如下的*.mcom/COM,ANSYSRELEASE 8.0UP2003093009:28:4207/23/2005/COM, truss.mcom！ X 方向LCOPER,ZERO TOC o 1-5

23、h z LCDEFI,1,1,1LCFACT,1,0.263825E-17LCASE,1LCOPER,SQUARELCDEFI,1,1,2LCFACT,1, 8.55778LCOPER,ADD,1,MULT,1LCDEFI,1,1,3LCFACT,1, -0.188669E-13LCOPER,ADD,1,MULT,1LCDEFI,1,1,4LCFACT,1, -0.871099E-15LCOPER,ADD,1,MULT,1LCDEFI,1,1,5LCFACT,1, -0.757013LCOPER,ADD,1,MULT,1 TOC o 1-5 h z LCDEFI,1,1,6LCFACT,1,0

24、.967307E-13LCOPER,ADD,1,MULT,1LCDEFI,1,1,7LCFACT,1,0.533141E-13LCOPER,ADD,1,MULT,1LCDEFI,1,1,8LCFACT,1, -0.203699LCOPER,ADD,1,MULT,1LCDEFI,1,1,9LCFACT,1,0.445795E-13LCOPER,ADD,1,MULT,1LCDEFI,1,1,10LCFACT,1, -0.387808E-13LCOPER,ADD,1,MULT,1LCOPER,SQRT/COM,ANSYS07/23/20RELEASE 8.0UP2003093009:28:4205/

25、COM, truss.mcom！ Y 方向的LCOPER,SQUARE！这里没有用zero清空 TOC o 1-5 h z LCDEFI,1,1,1LCFACT,1, 50.7528LCOPER,ADD,1,MULT,1LCDEFI,1,1,2LCFACT,1,0.887017E-14LCOPER,ADD,1,MULT,1LCDEFI,1,1,3LCFACT,1,0.612824E-13LCOPER,ADD,1,MULT,1LCDEFI,1,1,4LCFACT,1, -1.96484LCOPER,ADD,1,MULT,1LCDEFI,1,1,5LCFACT,1, -0.331613E-13LCOPER,ADD,1,MULT,1LCDEFI,1,1,6LCFACT,1,0.330459E-13LCOPER,ADD,1,MULT,1LCDEFI,1,1,7LCFACT,1, 0.366569LCOPER,ADD,1,MULT,1LCDEFI,1,1,8LCFACT,1, -0.976991E-13LCOPER,ADD,1,MULT,1LCDEFI,1,1,9LCFACT,1,0.417313E-13LCOPER,ADD,1,MULT,1LCDEFI,1,1,10LCFACT,1,0.401040E-13LCOPER