ANSYS地震反應(yīng)譜SRSS分析共24頁

1、ANSYS地震反應(yīng)譜SRSS分析我在ANSYS中作地震分解反應(yīng)譜分析，一次X方向，一次Y方向，他們要求是獨(dú)立互不干擾的，可是采用直進(jìn)行一次模態(tài)分析的話，他生成的*.mcom文件好像是包含了前面的計(jì)算結(jié)果，命令流如下：!進(jìn)入PREP7并建模/PREP7B=15 !基本尺寸A1=1000 !第一個(gè)面積A2=1000 !第二個(gè)面積A3=1000 !第三個(gè)面積ET,1,beam4 !二維桿單元R,1,0.25,0.0052,0.0052,0.5,0.5 !以參數(shù)形式的實(shí)參MP,EX,1,2.0E11 !楊氏模量mp,PRXY,1,0.3mp,dens,1,7.8e3N,1,-B,0,0 !定義結(jié)點(diǎn)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!進(jìn)入求解器，定義載荷和求解/SOLUD,1,ALL,0,2 !結(jié)點(diǎn) 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!地震影響系數(shù)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這里在進(jìn)行X方向的反應(yīng)譜分析以后，進(jìn)行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 !注意這里沒有清空數(shù)據(jù)庫(kù)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兩個(gè)方向地震的疊加，可是如果在座Y方向的地震以前把注釋掉的模態(tài)分析在做一下這樣的Y方向的地震的*.mcom就是：/COM, truss.mcomLCOPER,ZERO !注意這里清空數(shù)據(jù)庫(kù)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他不清空數(shù)據(jù)庫(kù)，需要進(jìn)行兩次模態(tài)分析,這很耗時(shí)間對(duì)于大型結(jié)構(gòu)，請(qǐng)大家討論討論如何處理呢？Re:討論：ANSYS地震反應(yīng)譜SRSS分析 本人是學(xué)土木工程的，平時(shí)主要用Patran+Nastran對(duì)結(jié)構(gòu)做線性分析， 偶爾使用Ansys對(duì)結(jié)構(gòu)做地震反應(yīng)譜分析，但對(duì)Ansys的命令流不熟悉。 依據(jù)本人的使用經(jīng)驗(yàn)看

20、，“一次X方向，一次Y方向，他們要求是獨(dú)立互 不干擾的”應(yīng)該這樣完成才對(duì)：1、點(diǎn)擊Solution-) Load Step Opts- Spectrum-Settings，在對(duì)話框中將 SEDX，SEDY，SEDZ 設(shè)置為 1，0，0 后；點(diǎn)擊 Solution-Solve-) Current LS；讀入生成*.mcom 文件完成 X 方向的反應(yīng)譜內(nèi)力計(jì)算。2、點(diǎn)擊 Solution-) Load Step Opts-Spectrum-) Settings，在對(duì)話框中將SEDX，SEDY，SEDZ設(shè)置為0，1，0后；點(diǎn)擊 Solution-Solve-Current LS；讀入生成*.mcom

21、文件完成Y方向的反 應(yīng)譜內(nèi)力計(jì)算。需要注意的是，在對(duì)話框中將SEDX，SEDY，SEDZ為1，1， 0表示地震動(dòng)的方向?yàn)榕cX軸成45度的方向，這顯然不是所希望的Re:討論：ANSYS地震反應(yīng)譜SRSS分析我也知道是這樣，可是他形成的*.mcom好像沒有不是這個(gè)意思啊。按你說的，“一次X方向，一次Y方向，他們要求是獨(dú)立互不干擾的”應(yīng) 該這樣完成才對(duì)：1、點(diǎn)擊 Solution-Load Step Opts-Spectrum-) Settings，在對(duì)話框中將SEDX，SEDY，SEDZ設(shè)置為1，0，0后；點(diǎn)擊 Solution-) Solve-) Current LS；讀入生成*.mcom文件完

22、成X方向的反 應(yīng)譜內(nèi)力計(jì)算。2、點(diǎn)擊 Solution-) Load Step Opts-)Spectrum-) Settings， 在對(duì)話框中將 SEDX，SEDY，SEDZ 設(shè)置為 0，1，0 后；點(diǎn)擊 Solution-)Solve-) Current LS；讀入生成*.mcom文件完成Y方向的反應(yīng)譜內(nèi)力計(jì)算。做后：他把Y方向的同樣寫入到同一個(gè)*.mcom中了，而且也沒有晴空數(shù)據(jù)庫(kù)，形成如下的*.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