1、結構地震反應譜分析實例在多位朋友的大力幫助下，經(jīng)過半個多月的努力，鄙人終于對結構地震反應譜分析有了一定的了解，現(xiàn)將 其求解步驟整理出來，以便各位參閱，同時，尚有一些問題，歡迎各位討論！為敘述方便，舉一簡單實例：在側水壓與頂部集中力作用下的柱子的地震反應譜分析，譜值為加速度反應譜，考慮X與Y向地震效 應作用。已知地震影響系數(shù)a與周期T的關系：a(T)= 0.4853*(0.4444+2.2222*T) 0T=0.04 秒0.4853*(0.10/T)人(-0.686) 0.04T=0.1 秒0.4853 0.1T = 1.2 秒0.4853*(1.2/T)人1.5 1.2T=4 秒以下是命令流程

2、序/filname,SPEC,1/PREP7!定義單元類型及材料特性ET4,45MPrEX,1/2.8E10MPrDENS,1/2.4E3MPrNUXY,1/0.18!建立模型BLOCK,0,1,0,1,0,5!網(wǎng)格剖分ESIZE,0.5VMESH,all/VIEW,廣0.3廣1,1EPLOTFINISH/SOLU!施加底部約束ASEL,LOC,Z,0DA,ALL,ALLALLSEL!施加自重荷載ACEL,0,0,10!進行模態(tài)求解ANTYPE,MODALMODOPT,LANB,30SOLVEFINISH!進行譜分析/SOLUANTYPE,SPECTRSPOPT,SPRS,30,YESSVTY

3、P，2 !加速度反應譜SED,1,1 !X 與 Y 向FREQ,0.2500,0.2632,0.2778,0.2941,0.3125,0.3333,0.3571,0.3846,0.4167FREQ,0.4545,0.5000,0.5556,0.6250,0.7143,0.8333,1.1111,2.0000,10.0000 FREQ,25.0000,1000.0000SV0.05,0.0797,0.0861,0.0934,0.1018,0.1114,0.1228,0.1362,0.1522,0.1716SV0.05,0.1955,0.2255,0.2642,0.3152,0.3851,0.48

4、53,0.4853,0.4853,0.4853 SV0.05,0.2588,0.2167SOLVEFINISH!進行模態(tài)求解（模態(tài)擴展）/SOLUANTYPE,MODALEXPASS,ONMXPAND,30,YES,0.005SOLVEFINISH!進行譜分析（合并模態(tài)）/SOLUANTYPE,SPECTRSRSS,0.15,dispSOLVEFINISH/POST1SETLIST !結果 1/INPmcomlcwrite,11LCASE,11PRRSOL, !結果 2 SEFIRSTPRRSOL, !結果 3 SETNEXTPRRSOL, !結果 4 SETNEXTPRRSOL, !結果 5

5、 SETNEXTPRRSOL, !結果 6 FINISH!靜力分析/SOLU ANTYPE,STATIC!施加水壓荷載 NSEL,LOC,Y；0NSEL,R,LOC,Z,0,5 SFGRAD,PRES,0Z0,-10000 SF,ALL,PRES,50000!施加集中荷載 NSEL,LOC,Y；0NSEL,R,LOC,z,5 FrALL,FY；10000ALLSEL EPLOT SOLVE FINISH/POST1 set,last lcwrite,12 Lcase,11 Lcoper,add,12 Lcwrite,13LCASE,12PRRSOL, !結果 7LCASE,13PRRSOL,

6、!結果 8FINISH以下是計算的結果 結果 1： （Results Summary）21.64721.647121.51121.51結果2：（單獨譜分析反力LCASE,11）VALUE 2467.9 2290.1 18384.結果3：（單獨譜分析反力SET；FIRST）VALUE 0.13334E+06-0.15785E+07-0.18819E-06結果4：（單獨譜分析反力SETJNEXT）VALUE -0.15785E+07-0.13334E+06-0.48918E-06結果5：（單獨譜分析反力SETJNEXT）VALUE -0.87805E+07 0.27008E+08 0.86846E

7、-07結果6：（單獨譜分析反力SETJNEXT）VALUE 0.27008E+08 0.87805E+07 0.79325E-06結果7：（單獨靜力分析反力LCASE,12）VALUE 0.22901E-08-0.15500E+06 0.12000E+06結果8：（譜分析與靜力分析疊加反力LCASE,13）VALUE 2467.9 -0.15271E+06 0.13838E+06以下是問題的討論1、模態(tài)提取數(shù)為30，即取前30階振型數(shù)，但在譜分析時得到的是4階，這4階是什么意思？2、 在單獨譜分析時，為何結果2、3、4、5、6會相差如此之大？（其應力和位移也是如此。）3、在進行譜分析的合并模態(tài)

8、步驟中，模態(tài)合并方法應該選用SRSS,0.15,DISP還是SRSS,0.15,ACEL （其對 結果影響很大）？4、在譜分析中，似乎應該將水壓力折算成各節(jié)點處的附加質量MASS21，但對于非豎值方向（如本例為Y 方向），MASS21單元如何加上去？5、在選?。x散）頻率和譜值時，似乎沒有一定的標準，而其基本原則是須選取具有代表性的點。6、模態(tài)分析中采用Subspace法與Block Lanczos法所得結果相差2%4%左右，用后一種方法計算更快， 但占用內存大50%左右。7、在網(wǎng)格剖分時采用VMESH和VSWEEP兩種方法，對靜力分析結果沒有影響，但對譜分析結果卻有一定 影響，何也？8、有朋友說可以“做模態(tài)分析之前先做打開預應力的靜力分析，然后再做打開預