減隔振橋梁的三維有限元分析

隨著社會經濟的快速發(fā)展，現代中心城市發(fā)生了毀滅性的地震，造成的經濟損失越來越嚴重。最近二三十年,全球發(fā)生的多次破壞性地震都造成了非常慘重的生命財產損失。一個很重要的原因是,橋梁工程在地震中遭到了嚴重的破壞,切斷了震區(qū)交通生命線,造成救災工作的巨大困難,使次生災害加重,從而導致了非常巨大的經濟損失。因此,對橋梁工程進行正確的抗震設計就顯得尤為重要。以某減隔振橋梁為例,為驗證其抗震設計的可靠性,通過具體建模及計算過程,對該橋進行了地震作用下的強度和變形驗算,并以E2地震作用下的彈塑性地震響應分析作為重點進行深入探討,以期得到一些有益的結論。1橋梁場地類別某減隔振橋梁,采用(25+42+28+20)m預應力混凝土箱梁,截面形式為單箱單室斜腹板箱形截面,箱梁中心處梁高2.2m,頂板寬13.5m,底板寬6.0m,兩側懸臂各長3.0m。采用縱向預應力體系。下部結構中間墩(編號為4、5、6號橋墩)采用雙柱式門形橋墩,墩柱采用矩形截面,距墩頂30cm處的兩柱之間設一道矩形截面橫系梁。交界墩(編號為3、7號橋墩)采用矩形墩上接蓋梁形式,蓋梁高1.8m。橋墩尺寸為2.2m×1.5m,墩中心間距4.5m,承臺尺寸5.5m×6.7m×2.0m。基礎采用鉆孔灌注樁的基礎形式,樁基直徑1.2m。根據《中國地震動參數區(qū)劃圖》,工程沿線的地震動峰值加速度為0.20g,該地震動峰值加速度所對應的地震基本烈度為Ⅷ度。根據《建筑抗震設計規(guī)范》,設計地震分組為第一組。根據波速測試成果和《公路橋梁抗震設計公路橋梁抗震設計細則》的相關規(guī)定計算得出場地現狀地面下20m深度范圍內的土層平均剪切波速值(vse)分別為244m/s和230m/s;場地覆蓋層厚度大于3m且小于50m的分界值。由前述兩項條件判別擬建橋梁場地的場地類別為Ⅱ類。該橋采用了彈塑性減隔振球型鋼支座,因此,依照《公路橋梁抗震設計細則》中6.1.3條,關于規(guī)則橋梁的定義,該橋屬于非規(guī)則B類橋梁。因此E2地震作用下應采用非線性時程計算方法進行驗算。2橋梁彈塑度分析2.1土體作用對結構響應的影響采用MidasCivil中的梁單元建立有限元模型,如圖1、圖2所示。3、7號墩為交界墩,4、5、6號墩為中間墩。上部結構離散為三維彈性空間梁單元;下部結構橋墩墩底采用纖維梁柱單元(塑性鉸),墩身采用彈性梁單元。由于結構考慮土體作用后,會改變上部結構的動力特性,如結構的基本周期延長,結構體系阻尼增大等,本橋采用群樁基礎,樁基與周圍土體相互作用非常顯著,不僅影響樁基的力學響應,也會影響整個橋梁結構的力學行為,且已有研究表明,根據具體相互作用問題的不同,土體對結構響應的影響并不一定是減小的。因此模型中應考慮樁-土相互作用的影響,樁側與土的相互作用采用“m”法模擬,通過MidasCivil中的節(jié)點彈性支承連接實現,樁底采用固定端約束。2.2混凝土結構性能計算模型中混凝土纖維的單軸本構模型采用修正的Kent-Park模型,考慮了混凝土在加、卸載過程中的受壓連續(xù)剛度和強度退化效應、受拉剛化效應以及混凝土開裂后的剛度、強度退化效應;并通過修改混凝土的材料特性參數,有效地考慮箍筋對混凝土的約束作用。鋼筋纖維的單軸本構模型采用了修正的Menegotto-Pinto本構模型,考慮了鋼筋反復加載過程中的Bauschinger效應和等向強化效應。上述2種動力彈塑性本構模型的應力-應變關系如圖3、圖4所示。2.3混凝土施工模型圖5給出了5號橋墩的鋼筋布置,鋼筋坐標為實際位置,外側70mm所對應的區(qū)域為保護層混凝土,中間部分區(qū)域為核心混凝土。保護層混凝土、核心混凝土及鋼筋分別采用不同的材料本構模型模擬,橋墩墩底截面離散化如圖6所示。其他橋墩截面采用相同方法進行纖維離散化處理。2.4水平抗力防線本橋在設計中采用了FLLGZ(Ⅱ)彈塑性減隔振球型鋼支座。該類型支座在正常使用情況下,不但具有普通球型支座承載能力大、轉角大、轉動靈活、壽命長等特點,還具有抗水平剪力、減震耗能的功能。此外,支座設置了3道水平抗力防線,即水平地震作用小于豎向承載力的10%時,剪力銷抵抗地震作用引起的水平力,約束向處于固定狀態(tài);當水平地震作用等于或大于豎向承載力的10%時,剪力銷被剪斷,約束向處于滑動狀態(tài),滑動后會因摩擦而耗能,外側懸臂棒因受彎變形而耗能,直至進入塑性狀態(tài);滑移量達到規(guī)定值時,剛性擋塊限制支座繼續(xù)滑移并承受水平力,保證支座的整體性,從而防止落梁破壞。以上3道水平抗力防線使得該類支座具有可靠的安全性能。在MidasCivil中,該減隔振支座需采用“滯后模型”進行模擬,支座布置如圖7所示。依照單向支座在約束方向提供剛度、雙向和固定支座在順橋向和橫橋向均提供剛度的原則進行設置。3項目結構前模態(tài)分析結果模態(tài)分析是檢驗動力計算模型正確性的重要手段,對本項目橋梁進行模態(tài)分析可以得到結構前5階模態(tài)分析結果,如表1所示。結構基本周期3.71s,頻率0.269Hz,1階振型為順橋向平動,如圖8所示。4地震作用下結構的彈塑化時間的分析4.1力學性能為得到偏于安全的地震位移反應結果,在E2地震作用下,延性構件(墩柱)的有效截面抗彎剛度應采用鋼筋首次屈服時的割線剛度,參照《公路橋梁抗震設計細則》6.1.6條可計算得到各橋墩在橫橋向和順橋向的剛度調整系數,如表2所示,再通過MidasCivil中的截面特性調整系數來實現截面剛度的修正。4.2反應譜設計及結果本橋橋址處的地震動峰值加速度為0.20g,該地震動峰值加速度所對應的地震基本烈度為Ⅷ度。屬Ⅱ類場地,地震動反應譜特征周期為0.35s(1區(qū))。采用傅里葉變化及不斷迭代擬合的方法來生成人工模擬地震動時程,具體生成人工模擬地震波的方法如下:(1)首先按照《公路橋梁抗震設計細則》,根據地震烈度、場地類別等設計參數確定設計反應譜,即期望反應譜;(2)根據期望反應譜近似地計算出人工地震波的功率譜;(3)再由功率譜得到的傅里葉幅值譜加上隨機相位做傅里葉逆變換并加上強度包絡線,得到近似人工地震波;(4)計算近似人工地震波的反應譜,并用期望反應譜與計算反應譜的比值修改傅里葉幅值譜,重新生成人工地震波,不斷進行迭代,直至反應譜在控制頻率點處的誤差處于允許誤差范圍內。采用上述方法得到的人工地震動擬合反應譜曲線如圖9所示。為考慮地震動的隨機性,進一步根據《公路橋梁抗震設計細則》,設計加速度時程不得小于3組,且應保證任意2組間同方向時程的相關系數絕對值小于0.1。經計算得到的每2條地震波的相關系數ρ最大值為0.051,均小于0.1,滿足規(guī)范要求。最終得到的3條E2地震作用下人工擬合地震動時程曲線如圖10所示。4.3驗算結果的分析4.3.1開發(fā)前驗算順橋向及橫橋向的塑性角區(qū)域最大轉角驗算依照《公路橋梁抗震設計細則》7.4.3條進行驗算,結果如表3所示,表中取3條人工波計算結果的最大值與容許值進行比較,來判斷是否滿足規(guī)范要求。經驗算,順橋向及橫橋向塑性角區(qū)域的最大轉角驗算均滿足規(guī)范要求。4.3.2地震波激勵下的柱墩位移墩頂位移的驗算,須按照《公路橋梁抗震設計細則》7.4.6條進行。由于本橋在橫橋向為雙柱墩,因此,橫橋向的墩頂位移容許值Δu應按規(guī)范7.4.8條建立彈塑性靜力分析模型,采用Midas進行pushover分析。Δd的計算參照6.7.6條。全橋模型在地震波激勵下,得到的各橋墩最大墩頂位移列于表4。經驗算,順橋向及橫橋向的墩頂位移均滿足規(guī)范要求。4.3.3橋梁墩抗剪強度E2地震作用下,按照《公路橋梁抗震設計細則》7.3.4條進行墩柱塑性鉸區(qū)域順橋向和橫橋向的斜截面抗剪強度驗算。首先按照公式7.3.4條計算各橋墩順橋向及橫橋向的抗剪強度,并列于表5中。根據纖維截面計算結果,在整個加載過程中,各橋墩墩底、墩頂最外層鋼筋均未進入屈服,因此,參照《公路橋梁抗震設計細則》6.8.1,橋梁墩柱的剪力設計值可用E2地震作用的計算結果。經表5對比后可知,各橋墩塑性鉸區(qū)域順橋向、橫橋向抗剪強度均滿足規(guī)范要求。5e地震作用下動力彈塑性分析本文以某減隔振橋梁為例,采用MidasCivil建立了該橋的三維有限元模型,根據混