劉陽,鐘沛杰,許一鵬,劉小娟,任彧,蘇龍輝
(1.華僑大學 土木工程學院,福建 廈門,361021;2.華僑大學 福建省結構工程與防災重點實驗室,福建 廈門,361021;3.福建建工集團有限責任公司,福建 福州,350001;4.中建協和建設有限公司,福建 泉州,362700)
目前,結構在遭遇罕遇地震后的“功能可恢復能力”備受關注。工程結構抗震已從抗倒塌能力設計逐漸發(fā)展為結構功能可恢復性能設計。
為實現結構的震后功能可恢復,國內外學者對不同的震后功能可恢復結構及構件開展研究[1]。FORTNEY 等[2]提出了可更換“保險絲”連梁的概念,進行了試驗驗證并提出了相應設計方法。滕軍等[3-7]先后提出了不同的可替換連梁構造并開展了相應的試驗和理論研究,認為所提出的可替換連梁構造均具有良好的抗震性能和震后可恢復性。OH 等[8]將開縫鋼板阻尼器設置于鋼梁端部的下翼緣,發(fā)現在地震作用下,鋼梁保持彈性,損傷集中于阻尼器,震后可方便更換阻尼器。CALADO等[9]在鋼梁端部設置螺栓連接的耗能鋼板,通過耗能鋼板的震后快速替換實現鋼框架的功能快速恢復。郭子雄等[10-11]提出了一種新型裝配式鋼筋混凝土柱-鋼梁框架節(jié)點構造,通過在節(jié)點區(qū)鋼板筒預埋高強螺栓與鋼梁連接,實現震后受損鋼梁的快速替換。
對于框架結構而言,在水平地震作用下,梁端可能產生塑性鉸,底層柱底也會因承受較大的彎矩而出現塑性鉸。因此,為實現框架結構體系的震后性能可恢復,除實現框架梁的可替換外,還需要解決底層柱腳塑性鉸區(qū)的可替換問題,然而,該方面的研究還鮮見報道。目前,已有一些學者針對橋梁結構開展了可更換橋墩的研究。MARRIOTT 等[12]在橋墩底部設置可更換軟鋼阻尼器用于耗散地震能量,同時采用后張拉預應力技術實現震后橋墩的自復位,通過對橋墩進行抗震性能試驗驗證該自復位搖擺橋墩構造的可行性。GUO 等[13]提出了一種適用于預制混凝土橋墩的震后可更換技術,通過在橋墩底部設置鋁合金耗能元件耗散地震能量,同時采用纖維增強復合材料作為預應力筋,張拉預應力以提供恢復力;通過對15 個橋墩試件進行抗震性能試驗驗證該構造的可行性,試驗結果表明,通過外置耗能元件和張拉預應力可以實現橋墩的震后可恢復性,提高預應力水平可以提高試件的自復位能力。呂英婷等[14]提出了一種帶位移型剪切鋼板的新型可原位修復鋼橋墩,通過有限元模擬驗證其震后可恢復性能。已有的關于可恢復性能橋墩的研究可為可恢復性能框架柱的研究提供有益的參考。但是與橋墩相比,框架柱的截面尺寸較小,且底層柱的軸壓比較大,橋墩的可更換技術不適用于框架柱,因此,有必要開展可恢復功能框架柱的研究。
此前,劉陽等[15-16]提出了幾種震損可替換組合柱,通過在柱底設置可更換耗能元件,使得柱身主要承受軸向荷載,而可更換耗能元件承受水平荷載作用下的彎矩和剪力。在地震作用下,塑性變形集中于可更換耗能元件而耗散地震能量,柱身保持彈性,從而實現組合柱的震后性能可恢復。本文在已有研究的基礎上,進一步將開縫鋼板阻尼器設置于柱腳,同時釋放柱身與基礎的直接連接,提出一種新型的震后可恢復功能搖擺柱(earthquake resilience rocking column,簡稱ERR柱),并進行足尺試驗研究[17]。本文在試驗研究的基礎上,對ERR 柱進行精細化有限元建模,對ERR 柱的滯回特性進行進一步的參數分析,在此基礎上提出ERR 柱的恢復力模型,該模型可應用于帶ERR 柱的功能可恢復框架結構的彈塑性地震反應分析。
ERR柱的構造示意圖如圖1所示,其有限元模型如圖2所示。所有構件均采用八節(jié)點線性六面體縮減積分單元(C3D8R)模擬。將開有洞口的鋼構件設置為六面體掃掠網格,并采用中性軸算法進行運算。在螺栓與鋼板連接部位,受力較為復雜,故設置較小網格邊長,最小網格邊長為2 mm 左右。承臺和柱身部分基本處于彈性狀態(tài),為提高計算效率,選擇較大網格邊長,為50~200 mm。
圖1 震后可恢復功能搖擺柱示意圖Fig.1 Schematic diagram of ERR column
圖2 ERR柱有限元模型Fig.2 Finite element model for ERR column
柱身、阻尼器、連接板等鋼材本構均采用理想彈塑性模型,屈服強度和彈性模量按照材性試驗結果取值。高強螺栓屈服強度為900 MPa,抗拉強度為1 000 MPa,采用理想彈塑性應力-應變曲線,泊松比均為0.3。
對于鋼板阻尼器兩側的螺栓,采用實體建模,螺栓與鋼板之間、柱身與承臺之間均采用切向“庫侖摩擦接觸”和法向“硬接觸”進行模擬。螺桿與孔壁側面之間僅考慮法向“硬接觸”,忽略切向摩擦力影響。螺栓預緊力通過“螺栓荷載”施加。當接觸面法向方向為“硬接觸”時,界面僅傳遞壓應力,不傳遞拉應力;切線方向采用罰函數庫侖摩擦模型模擬,根據GB 50017—2017“鋼結構設計規(guī)范”[18],取摩擦因數為0.3。鋼材焊接連接、T 形連接板和承臺之間的螺栓連接簡化為“綁定”約束。
承臺底面固端約束,在柱端設置參考點先施加豎向恒定荷載,按照圖3所示加載制度施加水平往復荷載。
圖3 水平荷載加載制度Fig.3 Horizontal loading protocol
為驗證有限元模型的合理性,選取文獻[17]中的試驗進行模型驗證。試驗中設計軸壓比n、鋼板阻尼器厚度t和剪跨比λ取值見表1,材料屬性見表2。試驗加載裝置及試件幾何尺寸見文獻[17]。表1中,試件編號S20-5.5-0.1 表示其鋼板阻尼器厚度為20 mm,剪跨比為5.5,設計軸壓比為0.1。
表1 試件參數Table 1 Parameters of specimens
表2 鋼材屬性Table 2 Steel properties
圖4 所示為試件S16-5.5-0.1 在位移角θ為1/20 rad 時的破壞形態(tài)與有限元模擬結果對比。在加載過程中,柱和連接板始終保持彈性。在加載后期,試件的損傷主要集中在鋼板阻尼器的水平條帶板上。由圖4可見,有限元模型可以很好地模擬試件的變形形態(tài)和損傷分布。
圖4 試件S16-5.5-0.1試驗與有限元模擬的破壞形態(tài)對比Fig.4 Comparison of failure modes obtained from test and FE models of sample S16-5.5-0.1
試驗得到的荷載-位移角滯回曲線與有限元計算結果對比如圖5 所示。圖5 中,θ為試件加載位移角;P為試件荷載。從圖5可以看出:試驗得到的滯回曲線捏攏,相比之下,有限元計算得到的滯回曲線比較飽滿。有限元計算結果與試驗結果存在差異的原因如下:
1) 試驗時,對高強螺栓施加的預緊力達不到設計要求,高強螺栓的受力機制較早由“摩擦型”向“承壓型”轉變,導致柱腳連接板與開縫鋼板阻尼器之間較早出現轉動滑移,如圖4(b)所示。而有限元模型則能根據設計要求施加預緊力,因此,柱腳連接板與鋼板阻尼器之間的轉動滑移比較小。
2) 實際試件安裝時存在一定的誤差,不同部件之間存在一定的縫隙,在加載過程中虛變形較大,也會導致滯回曲線出現一定程度的捏攏。而有限元建模時各部件的位置是理想的狀態(tài),虛變形較小,滯回曲線相對比較飽滿。
3) 在加載過程中,ERR 柱損傷集中于柱底可更換的鋼板阻尼器,柱身及連接板均保持彈性,為節(jié)約試件制作材料,所有的試件均重復使用同一套鋼柱及連接板,不同試件之間只更換了鋼板阻尼器。經歷多次往復加載后,柱根部連接板表面由于反復摩擦而發(fā)生磨損,摩擦因數逐漸降低,導致后續(xù)加載的試件滑移更加明顯,也造成試驗滯回曲線出現捏攏。對比圖5中4個試件的試驗滯回曲線可以發(fā)現,第一個加載的試件S16-5.5-0.1柱腳連接板與鋼板阻尼器之間的滑移現象還不明顯,其試驗滯回曲線比較飽滿,與有限元模擬結果吻合良好。
圖5 試件滯回曲線對比Fig.5 Comparison of hysteretic curves
總體上,所建立的有限元模型可以很好地模擬ERR柱的受力性能,可用于后續(xù)的參數分析。
軸壓比是影響試件滯回特性的主要參數[17]。由于試驗條件限制,試驗軸壓比最大設計值為0.2,所以采用有限元方法補充參數分析。建立阻尼板厚度為16 mm,剪跨比為5.5 的有限元模型,分析不同軸壓比對試件滯回性能的影響,選取軸壓比為范圍為0.1~0.6。圖6所示為不同軸壓比有限元模型的滯回曲線對比。圖6中,Pmax為模型正向加載時的峰值荷載;Pmin為模型反向加載時的峰值荷載;PT為變形恢復時(即柱底抬起量為0 mm)對應的荷載,即曲線上剛度突變點的荷載;Δ為位移;Δr為卸載至荷載為0 kN時對應的殘余變形。由圖6可見,隨著軸壓比的增加,模型滯回曲線由略顯捏攏的梭型逐漸轉為“旗幟型”。PT和Δr隨軸壓比的增加迅速減小,模型呈現出良好的自復位能力。
圖6 不同軸壓比有限元模型的滯回曲線對比Fig.6 Comparison of hysteretic curves for finite element models with different axial compression ratios
為更加清晰地描述試件滯回環(huán)的特征,以軸壓比不同的模型S16-5.5-0.1 和S16-5.5-0.5 為例進行說明。圖7 所示為模型在位移角為1/20 rad 時的滯回環(huán),其中Ek為滯回環(huán)包圍的面積。滯回環(huán)上特征點處對應的有限元塑性應變云圖如圖8所示。
圖8 不同特征點下ERR柱塑性應變分布Fig.8 Plastic strain distribution of ERR column at different characteristic points
由圖7可見,模型的滯回環(huán)變化可分為如下四個階段。
1) 彈性階段。從開始加載到模型屈服前的階段。由于鋼板阻尼器各水平條帶板不能同時達到屈服,因此,根據等能量法[19]確定試件的名義屈服點,如圖7中的a點和a′點所示。
圖7 滯回環(huán)特征點Fig.7 Characteristic points of hysteretic curve
2) 塑性發(fā)展階段。該階段即為從試件屈服至到達目標位移的階段。鋼板阻尼器屈服后進入強化階段,經歷了較大的塑性變形,使得屈服后模型的承載能力繼續(xù)增大,同時具有較好的耗能能力和變形能力。定義b或b′點為試件的峰值荷載點。
3) 變形復位階段。該階段即為從位移峰值點卸載直至柱頂水平變形恢復的階段。為說明ERR柱的自復位特點,可將變形復位階段分為以下三個子階段:
①階段Ⅰ,即從位移峰值點b(b′)卸載至阻尼板反向屈服點c(c′)的階段。在這個過程中,ERR 柱的轉動點不變,隨著水平力的減小,繞轉動點的傾覆力矩方向發(fā)生了變化,因此,阻尼板的剪力方向發(fā)生改變,甚至出現反向屈服,如圖7 和圖8中的c點和c′點所示。對于軸壓比較小的試件,由于軸力產生的力矩比較小,所以當阻尼板反向屈服時,水平力已經反向。而對于軸壓比較大的試件,由于軸壓力產生的力矩比較大,當阻尼板反向屈服時,水平力方向沒有發(fā)生改變。
②階段Ⅱ,即從阻尼板發(fā)生反向屈服點c(c′)至剛度突變點d(d′)的階段。在這個過程中,由于阻尼板發(fā)生了反向屈服,阻尼板產生較大的塑性變形,所以在這個階段荷載變化較小,變形恢復程度比較明顯,直到柱底抬起量為0 mm。此時,由于柱底剛性轉動點突然轉移到另一側,滯回曲線出現了明顯的拐點,表現為加載剛度的突增。從圖8中d點的受力分析可以看出,當軸壓比較小時,軸壓力不足以提供使試件變形復位的恢復力,這時需要施加反向的水平力以幫助試件復位。當軸壓比較大時,軸壓力足以提供使試件變形復位的恢復力。
③階段Ⅲ,即從柱底抬起量為0 mm到柱頂水平位移恢復到0 mm的階段。當柱底抬起量為0 mm時,柱身還有一些變形沒有恢復,這些變形主要是柱子在軸向壓力和彎矩復合作用下產生的變形。當軸壓力較大時,二階效應比較明顯,柱頂的水平變形比較大(如圖7中的d′點所示)。對于軸壓力小的試件,二階效應不明顯,當柱底抬起量為0 mm時,水平位移也基本恢復到0 mm(如圖7中的d點所示)。
反向加、卸載的過程與正向加、卸載類似。雖然兩個模型的滯回環(huán)形狀有所差異,但是由滯回曲線所包圍的面積(即耗能)差異不大,其耗能能力主要由鋼板阻尼器提供。
為進一步研究螺栓預緊力對試件滯回性能的影響,有限元分析中選取三個螺栓預緊力Tp,Tp最高值按照GB50017—2017“鋼結構設計標準”[18]取125 kN(10.9級M22高強螺栓),其余兩個分別為75 kN和50 kN。
考慮到滯回曲線的對比效果不夠直觀,為清晰展示螺栓預緊力對試件不同位移下承載力的影響,繪制不同螺栓預緊力試件的荷載-位移曲線,如圖9所示。由圖9可見,預緊力的降低對試件初始剛度影響不大,但在相對滑移產生后,模型承載力發(fā)展受到影響。隨著預緊力增大,骨架曲線上名義屈服位移增大,割線剛度減小。當摩擦型螺栓產生一定滑移后,螺桿與孔壁接觸,轉變?yōu)槌袎盒吐菟?,試件剛度有較大提升。最終,不同預緊力試件的荷載-位移曲線趨于一致,即當摩擦型螺栓都轉變成承壓型螺栓后,螺栓預緊力對模型的受力性能影響不大。
圖9 不同螺栓預緊力試件荷載-位移曲線Fig.9 P-Δ curves of specimens with different bolt prestress
由于螺栓滑移影響因素的定量分析需要考慮螺栓布置和孔壁間隙等多個因素,分析過程比較復雜;同時,也考慮到螺栓滑移僅對試件耗能能力有一定影響而對承載力影響不大,因此,后續(xù)建立恢復力模型時均假設螺栓群設計滿足要求,阻尼器與柱連接板之間不會發(fā)生相對滑移。
圖10 所示為阻尼器上部螺栓的螺桿與相鄰孔壁接觸應力τ隨加載位移Δ的變化曲線。由圖10可見,螺栓預緊力越小,螺栓由摩擦型轉變?yōu)槌袎盒蛯奈灰圃叫?,螺桿和孔壁的相互作用應力發(fā)展更快。在實際工程中,應通過加強螺栓群的設計,避免連接板和鋼板阻尼器之間發(fā)生相對滑移,以減小對試件耗能能力的不利影響。
圖10 螺桿與相鄰孔壁的接觸應力τ隨位移Δ的變化Fig.10 Changes of contact stress between screw and the adjacent hole wall with displacement
根據圖5~7所示的滯回曲線,ERR柱的恢復力特性可以用圖11 所示兩個模型進行描述。圖11中,Py和Δy分別為屈服點的荷載和位移;Pm和Δm分別為試件的峰值荷載及對應位移;PT和ΔT分別為反向加載剛度突增點的荷載和位移;括號內數字表示滯回曲線的加載路徑。由圖11 可見,建立恢復力模型的關鍵是上述三個特征點的確定。若卸載到鋼板阻尼器反向屈服時,水平力已改變方向,則恢復力模型如圖11(a)所示,此時ΔT=0 mm。反之,恢復力模型如圖11(b)所示,位移ΔT方向與加載方向相同。
圖11 ERR柱恢復力模型Fig.11 Restoring force model of ERR columns
4.1.1 名義屈服點
根據對試驗結果及有限元獲得的滯回曲線進行分析,發(fā)現ERR 柱名義屈服位移基本在位移角1/100 rad左右,受阻尼板厚度以及軸壓比等參數影響不大。為簡化模型,取Δy=H/100;H為水平加載點到基礎頂面的距離。
取圖12 所示隔離體進行受力分析,發(fā)現柱根除了豎向抬起變形δ,還存在轉動變形,因此兩側的水平條帶板分別存在水平向的拉力ΣTi和壓力ΣCi,如圖12 所示,其中i為水平條帶板編號。試件屈服時位移角較小,水平條帶板的拉、壓力較小,鋼板阻尼器近似處于平剪受力狀態(tài)。由于計算較為復雜,且在水平位移較小時誤差不大,為簡化計算,可認為水平條帶板只受剪力作用。對搖擺柱右下角轉動點取矩,由彎矩平衡可得:
圖12 屈服狀態(tài)下ERR柱受力分析Fig.12 Load resisting mechanism of ERR column under yielding
式中:N為試件所受軸力;a為柱截面邊長;j為阻尼器編號;m為阻尼器個數,本文計算模型僅考慮柱身左右兩側設置阻尼器,故m=2;xj為各阻尼器水平條帶板形心到轉動點的水平投影距離;Fq,s為帶縫鋼板阻尼器的屈服受剪承載力,其計算公式[8]如下:
式中:fy為鋼板阻尼器的屈服強度;L和B分別為鋼板阻尼器水平條帶板長度和寬度;r為倒角半徑;K為水平條帶板數量。
4.1.2 峰值荷載點
由于材料屈服后進入塑性強化階段,試件承載力在屈服后存在一定的提高,試件骨架曲線表現出屈服后強化的特性。試驗中,試件承載力在阻尼器水平條帶板產生撕裂后才開始緩慢降低,峰值荷載對應的位移角為1/20 rad。在有限元模擬中,鋼材本構簡化為理想彈塑性模型,因此,有限元模擬得到的試件滯回曲線未出現下降段,為此,取峰值荷載對應的位移Δm=H/20。
由于達到峰值荷載時,試件的位移角較大,柱腳發(fā)生明顯的轉動變形,因此,兩側鋼板阻尼器分別處于拉彎剪和壓彎剪復合受力狀態(tài),拉彎和壓彎產生的正應力對鋼板阻尼器的強度影響較大,此時式(1)~(4)不再適用。為此,對試件整體采用塑性極限法進行分析,受力簡圖如圖13所示。
根據能量守恒原理可以得到峰值荷載Pm計算公式如下:
式中:Ny和Mp分別為阻尼器條帶板的軸向屈服承載力和塑性極限彎矩;δi和αi分別為第i塊(i=1,2,3,4)水平條帶板的軸向變形和板帶端塑性鉸的轉動變形。對于圖13 中的柱腳抬起側即左側的鋼板阻尼器變形,按照式(6)~(10)進行計算:
式中:α=Δ/Η;xi和yi分別為各條帶板與柱身相連一側在x和y方向上的變形;ai為第i塊條帶板形心距離承臺表面的垂直距離,αiL和αiR分別為條帶板左端和條帶板右端的轉角。
對于柱腳未抬起側即圖13 中右側的阻尼器,其變形按式(11)~(15)進行計算:
圖13 峰值荷載點ERR柱化分析模型Fig.13 Simplified plastic model of ERR column at maximum load point
4.1.3 反向加載剛度突增點
由3.1 節(jié)分析可知,當卸載到位移接近0 mm時,抬起側柱腳底面逐漸接觸基礎頂面,阻尼器近似處于平剪狀態(tài),仍可采用圖12 所示應力狀態(tài)進行受力分析。因此,可以按照式(3)和(4)計算阻尼器的豎向剪力。當反向加載時,轉動點逐漸轉移到另一側,阻尼板的受力方向發(fā)生改變,反向加載剛度出現突變。圖13 所示工況為正向加載,柱腳左側抬起,轉動點位于右側。當反向加載時,轉動點將由右側轉移到左側,條帶板產生的內力將與圖13 所示方向相反。這個狀態(tài)的彎矩計算復雜,為簡化設計公式,引入調整系數η。
PT計算公式如下:
其中,水平荷載以向右為正,彎矩以順時針為正。由式(16)可見,當軸壓比較小時,計算得到的水平力為負值。這說明要使試件變形恢復,需要施加反向的水平荷載(見圖8(a)中c點對應的模型),試件的變形恢復能力弱。反之,當試件軸壓比n>0.3時,試件在軸力的作用下具有良好的變形恢復能力(見圖8(b)中d點對應的模型)。
為確定η的計算公式,以軸壓比(n=0.1~0.6)和阻尼器板厚(t=8~20 mm)為參數,針對24個在低周往復荷載作用下的ERR 柱進行有限元模擬。將計算得到的不同試件變形恢復時的水平承載力PT代入式(16),計算得到η隨不同參數的變化規(guī)律,如圖14 所示。由圖14 可見,η與阻尼器板厚和軸壓比基本呈線性關系。經雙參數線性擬合可得η的計算公式:
圖14 η計算結果Fig.14 Calculation results of η
其中:0.1≤n≤0.6,8≤t≤20。
根據試驗結果及有限元模擬結果可知,構件卸載時的剛度基本與初始彈性剛度相等,且在相同位移幅值下強度的衰減可以忽略不計。因此,在建立滯回規(guī)則時不考慮累積損傷引起的強度衰減及卸載剛度的退化。
ERR柱的滯回規(guī)則描述如下:
1) 在達到屈服點(Py,Δy)之前,沿骨架曲線的彈性段(即圖11 中的Oa段和Oq段)進行加載和卸載;
2) 當超過屈服點后,加載路徑沿骨架曲線前進(圖11 中的al和qs段)。從骨架曲線上任意一點卸載和重新加載,卸載剛度和重新加載剛度均取試件的初始彈性剛度ke;
3) 當卸載到PT時,試件剛度減小為0即荷載-位移曲線斜率為0。對于軸壓比小于0.3 的試件,當繼續(xù)卸載至位移為0 mm 時剛度發(fā)生突變(即圖11(a)中p點)。對于軸壓比大于0.3的試件,則在卸載至與骨架曲線相交時(圖11(b)中p點),試件剛度發(fā)生突變。
4) 反向加載均指向骨架曲線的反向屈服點(圖11中pq段);
5) 反向的加載和卸載規(guī)則同前;
6) 若加載超過最大荷載點時,則骨架曲線負剛度取初始剛度的1/10。
由于試驗中,試件加載到的最大位移角為1/15 rad,雖然此時試件強度衰減仍然較小(<10%),但本模型中的試件極限位移角基于安全考慮取為1/15 rad。
按照本文建議的恢復力模型計算試驗模型與有限元模型的滯回曲線,部分試件的滯回曲線對比如圖15 所示,可見根據恢復力模型得到的計算滯回曲線與試驗結果和有限元分析結果均吻合良好。由此可見,本文建立的恢復力模型可以較好預測ERR 柱的滯回性能,可用于帶ERR 柱結構的彈塑性分析。
圖15 計算滯回曲線與試驗、有限元結果對比Fig.15 Comparison of calculated hysteretic loops with experimental and FE results
1) 建立的三維有限元模型可以較好地反映ERR柱的破壞形態(tài)、應力分布和荷載-位移關系。
2) 軸壓比是影響ERR柱滯回特性的重要參數。隨著軸壓比的增加,ERR 柱滯回曲線由存在一定捏攏的“梭型”逐漸轉變?yōu)椤捌鞄眯汀保嚰目苫謴托阅茈S軸壓比的增加而顯著增加。當軸壓比大于0.3時,試件卸載后變形基本可以完全恢復。
3) 柱身與阻尼器之間的螺栓預緊力對試件剛度和耗能有一定影響。在相同條件下,隨螺栓預緊力的提高,螺栓由摩擦型向承壓型的轉變推遲,試件由此而發(fā)生的剛度降低現象也推遲,試件滯回曲線的“捏攏現象”不明顯。因此,在實際應用中建議采用高強摩擦型螺栓并按照“強連接、弱阻尼器”原則進行試件設計。
4) 本文所提ERR 柱恢復力模型的計算結果與試驗和有限元結果吻合良好,可應用于該功能可恢復結構的彈塑性地震反應分析。