金字塔的結構設計,我打算要找第二類地盤,以反應譜法做結構動力分析。根據建築物耐震設計規範及解說(2001),採用其第二類地盤正規化加速度反應譜做計算,依照第36頁所述乘上重力加速度以取得加速度反應譜。1999年09月21日凌晨1:46分,我擁有九二一集集大地震全國600多個測站的地震儀加速度數據,我應該會斟酌是否做歷時結構動力分析。我希望地震工程的教授,得另外提供水平1200GAL的反應譜,以及垂直450GAL的反應譜。
中國大陸的反應譜設置原則,得輸入加速度值求算需求反應譜,此反應譜方法可用於1200GAL。以下,為臺灣各類地盤的結構動力分析的代碼,應該照常採用SRSS對其結果做求算,對於模態週期相近者採用CQC做計算。ANSYS的SRSS和CQC求算,為採用指令名稱SRSS和CQC(無須參數)。此外,反應譜計算要先取得模態數據,以下為反應譜分析前的振態求解代碼。「MODOPT,LANB,20」,選擇常用求解法Block Lanczos,並採用求解出20個振態,以便於在反應譜求解之後,要選用SRSS或CQC。振態不夠用時,應該修改振態數重新求解,我不傾向於求解後擴展振態,也就是振態數不一定是20個。
振態分析
/SOLU
ANTYPE,MODAL
MODOPT,LANB,20
SOLVE
FINISH
鋼結構
第一類地盤(堅實地盤)
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,1
SVTYP,2
FREQ,0,0.03,0.15,0.333,0.4,0.5,0.6,0.7,0.8
FREQ,0.9,1.0,1.1,1.2,1.315,5
SV,, 9.80665, 9.80665,32.68883,32.68883,28.90241,24.90734,22.05666,19.90255,18.20737
SV,,16.83238,15.69064,14.72467,13.89483,13.07553,13.07553
SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,0.288,0.4,0.5,0.6,0.7,0.8
FREQ,0.9,1,1.139,5
SV,, 9.80665, 9.80665,35.95772,35.95772,28.90241,24.90734,22.05666,19.90255,18.20737
SV,,16.83238,15.69064,14.38309,14.38309
SOLVE
FINISH
第二類地盤(普通地盤)
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,1
SVTYP,2
FREQ,0,0.03,0.15,0.465,0.5,0.6,0.7,0.8,0.9
FREQ,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.837
FREQ,5
SV,, 9.80665, 9.80665,32.68883,32.68883,31.13417,27.57083,24.87818,22.75922,21.04048
SV,,19.61330,18.40584,17.36853,16.46601,15.67227,14.96775,14.33741,13.76950,13.07553
SV,,13.07553
SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,0.403,0.5,0.6,0.7,0.8,0.9
FREQ,1,1.1,1.2,1.3,1.4,1.5,1.592,5
SV,, 9.80665, 9.80665,35.95772,35.95772,31.13417,27.57083,24.87818,22.75922,21.04048
SV,,19.61330,18.40584,17.36853,16.46601,15.67227,14.96775,14.38309,14.38309
SOLVE
FINISH
第三類地盤(軟弱地盤)
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,1
SVTYP,2
FREQ,0,0.03,0.2,0.611,0.8,1,1.2,1.4,1.6
FREQ,1.8,2,2.2,2.415,5
SV,, 9.80665, 9.80665,32.68883,32.68883,27.31106,23.53596,20.84224,18.80673,17.20489
SV,,15.90561,14.82673,13.91394,13.07553,13.07553
SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,0.53,0.7,0.9,1.1,1.3,1.5
FREQ,1.7,1.9,2.093,5
SV,, 9.80665, 9.80665,35.95772,35.95772,29.85382,25.24858,22.08700,19.75921,17.96130
SV,,16.52340,15.34250,14.38309,14.38309
SOLVE
FINISH
臺北盆地
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,1
SVTYP,2
FREQ,0,0.03,0.2,1.32,1.5,1.7,1.9,2.1,2.3 FREQ,2.5,2.7,2.9,3.1,3.3,5 SV,, 9.80665, 9.80665,32.68883,32.68883,28.76617,25.38192,22.71014,20.54727,18.76055 SV,,17.25970,15.98121,14.87906,13.91912,13.07553,13.07553 SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,1.32,1.5,1.7,1.9,2.1,2.3
FREQ,2.5,2.7,2.9,3.1,3.3,5
SV,, 9.80665, 9.80665,32.68883,32.68883,28.76617,25.38192,22.71014,20.54727,18.76055
SV,,17.25970,15.98121,14.87906,13.91912,13.07553,13.07553
SOLVE
FINISH
鋼筋混凝土
第一類地盤(堅實地盤)
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,1
SVTYP,2
FREQ,0,0.03,0.15,0.333,0.4,0.5,0.6,0.7,0.8
FREQ,0.9,1.0,1.1,1.2,1.315,5
SV,, 9.80665, 9.80665,24.51663,24.51663,21.67680,18.68050,16.54250,14.92691,13.65553
SV,,12.62429,11.76798,11.04350,10.42112, 9.80665, 9.80665
SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,0.288,0.4,0.5,0.6,0.7,0.8
FREQ,0.9,1,1.139,5
SV,, 9.80665, 9.80665,26.96829,26.96829,21.67680,18.68050,16.54250,14.92691,13.65553
SV,,12.62429,11.76798,10.78732,10.78732
SOLVE
FINISH
第二類地盤(普通地盤)
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,1
SVTYP,2
FREQ,0,0.03,0.15,0.465,0.5,0.6,0.7,0.8,0.9
FREQ,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.837
FREQ,5
SV,, 9.80665, 9.80665,24.51663,24.51663,23.35063,20.67812,18.65864,17.06941,15.78036
SV,,14.70998,13.80438,13.02640,12.34951,11.75420,11.22581,10.75306,10.32712, 9.80665
SV,, 9.80665
SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,0.403,0.5,0.6,0.7,0.8,0.9
FREQ,1,1.1,1.2,1.3,1.4,1.5,1.592,5
SV,, 9.80665, 9.80665,26.96829,26.96829,23.35063,20.67812,18.65864,17.06941,15.78036
SV,,14.70998,13.80438,13.02640,12.34951,11.75420,11.22581,10.78732,10.78732
SOLVE
FINISH
第三類地盤(軟弱地盤)
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,1
SVTYP,2
FREQ,0,0.03,0.2,0.611,0.8,1,1.2,1.4,1.6
FREQ,1.8,2,2.2,2.415,5
SV,, 9.80665, 9.80665,24.51663,24.51663,20.48330,17.65197,15.63168,14.10505,12.90367
SV,,11.92921,11.12004,10.43546, 9.80665, 9.80665
SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,0.53,0.7,0.9,1.1,1.3,1.5
FREQ,1.7,1.9,2.093,5
SV,, 9.80665, 9.80665,26.96829,26.96829,22.39036,18.93643,16.56525,14.81941,13.47097
SV,,12.39255,11.50688,10.78732,10.78732
SOLVE
FINISH
臺北盆地
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,1
SVTYP,2
FREQ,0,0.03,0.2,1.32,1.5,1.7,1.9,2.1,2.3 FREQ,2.5,2.7,2.9,3.1,3.3,5 SV,, 9.80665, 9.80665,24.51663,24.51663,21.57463,19.03644,17.03260,15.41045,14.07041 SV,,12.94478,11.98591,11.15929,10.43934, 9.80665, 9.80665 SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,1.32,1.5,1.7,1.9,2.1,2.3
FREQ,2.5,2.7,2.9,3.1,3.3,5
SV,, 9.80665, 9.80665,24.51663,24.51663,21.57463,19.03644,17.03260,15.41045,14.07041
SV,,12.94478,11.98591,11.15929,10.43934, 9.80665, 9.80665
SOLVE
FINISH
振態疊加
SRSS
/SOLU
ANTYPE,SPECTR
SRSS
SOLVE
FINISH
CQC
/SOLU
ANTYPE,SPECTR
CQC
SOLVE
FINISH
建模方向,以下代碼採取Z向為垂直高度,水平地表加速度欲更換為另一方向,必須變更輸入代碼的參數為「SED,,1」,欲以相反方向做結構動力分析為「SED,,-1」。阻尼比,為舊結構動力分析必須定義的數據,鋼結構是0﹒02,鋼筋混凝土是0﹒05,金字塔是鋼結構建築物,因而阻尼比設定是「DMPRAT,0.05」。「SVTYP,2」,指的是加速度反應譜,此為不變更的指令代碼,因為結構設計都是加速度譜。ANSYS,可以做多點反應譜結構動力分析,但是真實世界幾乎無此類結構設計,此為單點反應譜結構動力分析(SPRS)。
Reference
中國大陸的反應譜設置原則,得輸入加速度值求算需求反應譜,此反應譜方法可用於1200GAL。以下,為臺灣各類地盤的結構動力分析的代碼,應該照常採用SRSS對其結果做求算,對於模態週期相近者採用CQC做計算。ANSYS的SRSS和CQC求算,為採用指令名稱SRSS和CQC(無須參數)。此外,反應譜計算要先取得模態數據,以下為反應譜分析前的振態求解代碼。「MODOPT,LANB,20」,選擇常用求解法Block Lanczos,並採用求解出20個振態,以便於在反應譜求解之後,要選用SRSS或CQC。振態不夠用時,應該修改振態數重新求解,我不傾向於求解後擴展振態,也就是振態數不一定是20個。
振態分析
/SOLU
ANTYPE,MODAL
MODOPT,LANB,20
SOLVE
FINISH
鋼結構
第一類地盤(堅實地盤)
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,1
SVTYP,2
FREQ,0,0.03,0.15,0.333,0.4,0.5,0.6,0.7,0.8
FREQ,0.9,1.0,1.1,1.2,1.315,5
SV,, 9.80665, 9.80665,32.68883,32.68883,28.90241,24.90734,22.05666,19.90255,18.20737
SV,,16.83238,15.69064,14.72467,13.89483,13.07553,13.07553
SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,0.288,0.4,0.5,0.6,0.7,0.8
FREQ,0.9,1,1.139,5
SV,, 9.80665, 9.80665,35.95772,35.95772,28.90241,24.90734,22.05666,19.90255,18.20737
SV,,16.83238,15.69064,14.38309,14.38309
SOLVE
FINISH
第二類地盤(普通地盤)
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,1
SVTYP,2
FREQ,0,0.03,0.15,0.465,0.5,0.6,0.7,0.8,0.9
FREQ,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.837
FREQ,5
SV,, 9.80665, 9.80665,32.68883,32.68883,31.13417,27.57083,24.87818,22.75922,21.04048
SV,,19.61330,18.40584,17.36853,16.46601,15.67227,14.96775,14.33741,13.76950,13.07553
SV,,13.07553
SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,0.403,0.5,0.6,0.7,0.8,0.9
FREQ,1,1.1,1.2,1.3,1.4,1.5,1.592,5
SV,, 9.80665, 9.80665,35.95772,35.95772,31.13417,27.57083,24.87818,22.75922,21.04048
SV,,19.61330,18.40584,17.36853,16.46601,15.67227,14.96775,14.38309,14.38309
SOLVE
FINISH
第三類地盤(軟弱地盤)
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,1
SVTYP,2
FREQ,0,0.03,0.2,0.611,0.8,1,1.2,1.4,1.6
FREQ,1.8,2,2.2,2.415,5
SV,, 9.80665, 9.80665,32.68883,32.68883,27.31106,23.53596,20.84224,18.80673,17.20489
SV,,15.90561,14.82673,13.91394,13.07553,13.07553
SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,0.53,0.7,0.9,1.1,1.3,1.5
FREQ,1.7,1.9,2.093,5
SV,, 9.80665, 9.80665,35.95772,35.95772,29.85382,25.24858,22.08700,19.75921,17.96130
SV,,16.52340,15.34250,14.38309,14.38309
SOLVE
FINISH
臺北盆地
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,1
SVTYP,2
FREQ,0,0.03,0.2,1.32,1.5,1.7,1.9,2.1,2.3 FREQ,2.5,2.7,2.9,3.1,3.3,5 SV,, 9.80665, 9.80665,32.68883,32.68883,28.76617,25.38192,22.71014,20.54727,18.76055 SV,,17.25970,15.98121,14.87906,13.91912,13.07553,13.07553 SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.02
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,1.32,1.5,1.7,1.9,2.1,2.3
FREQ,2.5,2.7,2.9,3.1,3.3,5
SV,, 9.80665, 9.80665,32.68883,32.68883,28.76617,25.38192,22.71014,20.54727,18.76055
SV,,17.25970,15.98121,14.87906,13.91912,13.07553,13.07553
SOLVE
FINISH
鋼筋混凝土
第一類地盤(堅實地盤)
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,1
SVTYP,2
FREQ,0,0.03,0.15,0.333,0.4,0.5,0.6,0.7,0.8
FREQ,0.9,1.0,1.1,1.2,1.315,5
SV,, 9.80665, 9.80665,24.51663,24.51663,21.67680,18.68050,16.54250,14.92691,13.65553
SV,,12.62429,11.76798,11.04350,10.42112, 9.80665, 9.80665
SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,0.288,0.4,0.5,0.6,0.7,0.8
FREQ,0.9,1,1.139,5
SV,, 9.80665, 9.80665,26.96829,26.96829,21.67680,18.68050,16.54250,14.92691,13.65553
SV,,12.62429,11.76798,10.78732,10.78732
SOLVE
FINISH
第二類地盤(普通地盤)
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,1
SVTYP,2
FREQ,0,0.03,0.15,0.465,0.5,0.6,0.7,0.8,0.9
FREQ,1,1.1,1.2,1.3,1.4,1.5,1.6,1.7,1.837
FREQ,5
SV,, 9.80665, 9.80665,24.51663,24.51663,23.35063,20.67812,18.65864,17.06941,15.78036
SV,,14.70998,13.80438,13.02640,12.34951,11.75420,11.22581,10.75306,10.32712, 9.80665
SV,, 9.80665
SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,0.403,0.5,0.6,0.7,0.8,0.9
FREQ,1,1.1,1.2,1.3,1.4,1.5,1.592,5
SV,, 9.80665, 9.80665,26.96829,26.96829,23.35063,20.67812,18.65864,17.06941,15.78036
SV,,14.70998,13.80438,13.02640,12.34951,11.75420,11.22581,10.78732,10.78732
SOLVE
FINISH
第三類地盤(軟弱地盤)
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,1
SVTYP,2
FREQ,0,0.03,0.2,0.611,0.8,1,1.2,1.4,1.6
FREQ,1.8,2,2.2,2.415,5
SV,, 9.80665, 9.80665,24.51663,24.51663,20.48330,17.65197,15.63168,14.10505,12.90367
SV,,11.92921,11.12004,10.43546, 9.80665, 9.80665
SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,0.53,0.7,0.9,1.1,1.3,1.5
FREQ,1.7,1.9,2.093,5
SV,, 9.80665, 9.80665,26.96829,26.96829,22.39036,18.93643,16.56525,14.81941,13.47097
SV,,12.39255,11.50688,10.78732,10.78732
SOLVE
FINISH
臺北盆地
結構動力分析(反應譜法、水平地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,1
SVTYP,2
FREQ,0,0.03,0.2,1.32,1.5,1.7,1.9,2.1,2.3 FREQ,2.5,2.7,2.9,3.1,3.3,5 SV,, 9.80665, 9.80665,24.51663,24.51663,21.57463,19.03644,17.03260,15.41045,14.07041 SV,,12.94478,11.98591,11.15929,10.43934, 9.80665, 9.80665 SOLVE
FINISH
結構動力分析(反應譜法、垂直地表加速度)
/SOLU
ANTYPE,SPECTR
DMPRAT,0.05
SED,,,1
SVTYP,2
FREQ,0,0.03,0.1,1.32,1.5,1.7,1.9,2.1,2.3
FREQ,2.5,2.7,2.9,3.1,3.3,5
SV,, 9.80665, 9.80665,24.51663,24.51663,21.57463,19.03644,17.03260,15.41045,14.07041
SV,,12.94478,11.98591,11.15929,10.43934, 9.80665, 9.80665
SOLVE
FINISH
振態疊加
SRSS
/SOLU
ANTYPE,SPECTR
SRSS
SOLVE
FINISH
CQC
/SOLU
ANTYPE,SPECTR
CQC
SOLVE
FINISH
建模方向,以下代碼採取Z向為垂直高度,水平地表加速度欲更換為另一方向,必須變更輸入代碼的參數為「SED,,1」,欲以相反方向做結構動力分析為「SED,,-1」。阻尼比,為舊結構動力分析必須定義的數據,鋼結構是0﹒02,鋼筋混凝土是0﹒05,金字塔是鋼結構建築物,因而阻尼比設定是「DMPRAT,0.05」。「SVTYP,2」,指的是加速度反應譜,此為不變更的指令代碼,因為結構設計都是加速度譜。ANSYS,可以做多點反應譜結構動力分析,但是真實世界幾乎無此類結構設計,此為單點反應譜結構動力分析(SPRS)。
Reference
- 李立峰、王連華編著,邵旭東、孫全勝主審(2015)ANSYS土木工程實例詳解,人民郵電出版社