俞 涵，潘加寶，駱乾輝，程 博，陶 韜，施遠(yuǎn)航

1.寧波市房屋建筑設(shè)計研究院有限公司，浙江寧波 315100

2.寧波東衡工程科技發(fā)展有限公司，浙江寧波 315100

0 引言

剪切型軟鋼阻尼器利用低屈服點耗能鋼板平面內(nèi)剪切彈塑性變形來消耗地震能量以達(dá)到消能減震的目的。因其具有易屈服、相對尺寸小、耗能能力強(qiáng)、造價低廉等特性，國內(nèi)外的研究學(xué)者對其展開大量研究，且在工程中得到了廣泛應(yīng)用。LIU Yang 與DENG Kailai［1-2］等對剪切型低屈服點鋼板阻尼器核心板進(jìn)行切削優(yōu)化設(shè)計，對比分析循環(huán)荷載作用下的性能差異，發(fā)現(xiàn)對核心板合理地切削后，阻尼器耗能性能降低較小，但卻能更好地提高其低周疲勞性能。Abebe等［3-4］對剪切型低屈服點鋼阻尼器的失效模式以及耗能性能進(jìn)行了試驗研究和模擬分析，研究結(jié)果表明剪切型阻尼器的損壞失效模式主要有三種類型：一是耗能板屈曲開裂破壞失效；二是核心板與上下端板以及側(cè)翼緣板的焊縫開裂失效破壞；三是端板翼緣板焊縫開裂失效。

阻尼器需兼具優(yōu)越的耗能性能及疲勞性能。不同于常規(guī)經(jīng)歷數(shù)千次循環(huán)荷載作用下不會發(fā)生破壞的低周疲勞破壞，在大震工況下阻尼器局部受到極大應(yīng)變幅往復(fù)循環(huán)作用并屈服，僅經(jīng)歷數(shù)十次甚至十次以內(nèi)的循環(huán)作用即發(fā)生斷裂，其疲勞破壞是一種很典型的延性斷裂破壞，屬于高應(yīng)變超低周疲勞［5］。建筑結(jié)構(gòu)的壽命與安全性評價極為重要，因此對剪切型軟鋼阻尼器的低周疲勞壽命研究是十分必要的。本文應(yīng)用ABAQUS軟件對剪切型軟鋼阻尼器結(jié)構(gòu)進(jìn)行建模分析，根據(jù)阻尼器的應(yīng)力分布云圖和等效塑性應(yīng)變云圖確定危險位置，再經(jīng)FE-SAFE 軟件根據(jù)應(yīng)力及應(yīng)變計算結(jié)果進(jìn)行疲勞分析，得到阻尼器壽命分布云圖及危險區(qū)域分布狀況，對剪切型軟鋼阻尼器疲勞研究具有一定的參考意義，為相關(guān)類型阻尼器進(jìn)一步改進(jìn)提供了參考。

1 材性試驗

阻尼器腹板使用BLY160鋼材，加勁肋及外部約束構(gòu)件均采用Q355b鋼材。為確定鋼材的力學(xué)性能，并為下文提供有限元分析材性數(shù)據(jù)，本文對BLY160進(jìn)行單軸拉伸試驗。依據(jù)《金屬材料 拉伸試驗 第1 部分：室溫試驗方法（GB/T 228.1—2021》制作試驗件并進(jìn)行單軸拉伸試驗（圖1）。

圖1 BLY160材性試驗

對鋼材進(jìn)行單軸拉伸試驗得到的應(yīng)力-應(yīng)變數(shù)據(jù)為名義應(yīng)力-應(yīng)變數(shù)據(jù)，在模擬非線性有限元定義構(gòu)件的彈性行為和塑性行為時，需要將其轉(zhuǎn)換為真實應(yīng)力-應(yīng)變數(shù)據(jù)。材料名義應(yīng)力-應(yīng)變與真實應(yīng)力-應(yīng)變轉(zhuǎn)換采用如下公式：

式（1～2）中：L為試件原始長度；

ε、σ分別為真實應(yīng)變、真實應(yīng)力；

εnom、σnom分別為名義應(yīng)變、名義應(yīng)力。

試驗結(jié)果見表1，名義應(yīng)力-應(yīng)變及真實應(yīng)力-應(yīng)變曲線見圖2。

表1 BLY160鋼材的屈服強(qiáng)度、抗拉強(qiáng)度及伸長率

圖2 BLY160應(yīng)力-應(yīng)變曲線

2 低周加載試驗

2.1 試驗試件尺寸

試驗的兩類阻尼器尺寸見圖3，（MD-250-1B僅在腹板角部引入尺寸為10 mm的倒角）。

圖3 阻尼器平面幾何尺寸（單位：mm）

2.2 試驗設(shè)備

試驗系統(tǒng)由一組平行四邊形的鋼加載框架、1 000 kN伺服作動器、連桿式位移傳感器、YHD-100型位移計（量程±100 mm，測量精度0.01mm）等組成，試驗加載見圖4。

圖4 試驗加載

2.3 試驗方案

試驗加載制度采用位移控制的方式，試驗加載共計5個工況，每個工況循環(huán)三周圈，加載波形采用三角波，各工況詳細(xì)加載位移見圖5。

圖5 加載制度

2.4 試驗結(jié)果及分析

以阻尼器試件耗能板的屈服產(chǎn)生較大面外變形以及各板件焊縫破壞為阻尼器破壞依據(jù)。本次試驗由于試件均焊接了加勁肋板將耗能板分為小的區(qū)格，各個試驗件的面外屈服現(xiàn)象并不明顯。分級加載第一級各試件均處于彈性階段，試件沒有屈服；在第二個加載等級，加載位移達(dá)到了12.5 mm，試件逐漸進(jìn)入屈服階段，在這個階段阻尼器并沒有產(chǎn)生明顯的面外屈服現(xiàn)象。阻尼器也沒有產(chǎn)生破壞，且表現(xiàn)出了明顯的強(qiáng)化現(xiàn)象；進(jìn)入第三個加載等級，阻尼器的加載位移達(dá)到了20 mm，阻尼器的強(qiáng)化現(xiàn)象減弱，剪切荷載峰值的增長幅度變緩；進(jìn)入第四個加載等級，阻尼器的加載位移達(dá)到了25 mm。幾類阻尼器均在第四個加載等級開始產(chǎn)生破壞，焊縫處開始產(chǎn)生裂紋，隨著試驗的進(jìn)行，裂紋擴(kuò)張裂開，破壞主要集中在耗能板邊緣、翼緣板焊縫，以及加勁肋板焊縫處，試件破壞圖見圖6。

圖6 阻尼器試驗破壞

3 有限元模型的建立與驗證

3.1 有限元模型的建立

本文利用非線性有限元分析軟件ABAQUS 對MD-250-1A 和MD-250-1B 進(jìn)行數(shù)值分析并與擬靜力試驗的結(jié)果進(jìn)行對比，驗證該有限元模型的準(zhǔn)確性與可靠性，并為后續(xù)參數(shù)化分析提供依據(jù)。在模型中，鋼材選用Von-Mises屈服準(zhǔn)則及相關(guān)流動法則，BLY160鋼材本構(gòu)選用多線性隨動強(qiáng)化模型（MKIN）和組合強(qiáng)化模型（CHAB &BISO），通常來講，殼單元是指結(jié)構(gòu)中某一個方向的尺寸遠(yuǎn)小于其他方向的尺寸，同時沿該方向的應(yīng)力幾乎可以忽略，阻尼器腹板厚度僅為10 mm，遠(yuǎn)小于其高度和寬度，因此剪切腹板采用精度較高的殼單元S4R，翼緣板及加勁板采用減縮積分實體單元C3D8R，在保障精度滿足要求的同時減少自由度。模型各個板件之間采用Tie連接，網(wǎng)格尺寸取10 mm，厚度方向網(wǎng)格不少于4層。構(gòu)件模型的網(wǎng)格劃分見圖7。材料強(qiáng)化參數(shù)數(shù)據(jù)見表1，連接板件采用理想彈塑性模型，屈服強(qiáng)度為355 MPa。加載方式同試驗。

圖7 有限元非網(wǎng)格劃分以及支撐組成部件

剪切型阻尼器各板件之間焊接連接，因此初始缺陷有焊接殘余應(yīng)力缺陷和幾何外形缺陷且不可避免，通常板件在靠近焊縫區(qū)域的為殘余拉應(yīng)力，遠(yuǎn)離焊縫區(qū)域的為殘余壓應(yīng)力。本文根據(jù)陳之毅等人對剪切鋼板阻尼器的研究成果按照圖8所示通過預(yù)先施加應(yīng)力引入焊接殘余應(yīng)力。本文采用一致模態(tài)法根據(jù)彈性屈曲變形的變形特征來假設(shè)各板件的初始變形大小及分布模式。首先對阻尼器進(jìn)行屈曲模態(tài)分析，選取1～4階模態(tài)以核心板厚度3%的比例作為初始缺陷引入模型中。圖9為初始幾何缺陷模態(tài)引入圖。

圖8 焊接殘余應(yīng)力

圖9 阻尼器缺陷引入屈曲模態(tài)

3.2 試驗與數(shù)值分析結(jié)果對比分析

試驗滯回曲線和數(shù)值模擬滯回曲線對比見圖10，兩者滯回曲線形狀相似程度十分接近，均呈紡錘形，無明顯強(qiáng)度和剛度突變，且差距較小。

圖10 剪切型阻尼器試驗與模擬滯回曲線對比

MD-S-250-1A 試驗最大承載力為454.06 kN，多線性強(qiáng)化模型模擬最大承載力為468.8 kN，組合強(qiáng)化模型模擬最大承載力為446.5 kN，前者與試驗結(jié)果相差3.2%，后者與試驗結(jié)果相差-1.6%；MD-S-250-1B試驗最大承載力為498.06 kN，多線性強(qiáng)化模型模擬結(jié)果為469 kN，與試驗結(jié)果相差0.5%，組合強(qiáng)化模擬結(jié)果為446.06 kN，與試驗相差1.04%。對于阻尼器的極限剪切承載力而言，兩種模型均可準(zhǔn)確地模擬預(yù)測，但多線性強(qiáng)化模型模擬滯回環(huán)呈明顯的對稱特性，雖然反映出了等幅循環(huán)加載歷程中塑性累計剪切力不斷強(qiáng)化最終至平穩(wěn)的強(qiáng)化歷程，但是對于循環(huán)加載的包辛格效應(yīng)卻不能較好的表現(xiàn)，導(dǎo)致滯回曲線過于飽滿，對阻尼器的耗能性能評估過高。而組合強(qiáng)化模型較好地模擬出了試件在彈塑性階段的非線性過渡，在經(jīng)歷正向加載后其反向加載的包辛格效應(yīng)及循環(huán)強(qiáng)化規(guī)律也較準(zhǔn)確，對于阻尼器的耗能性能合理評估相對多線性強(qiáng)化更加準(zhǔn)確。綜合分析，采用組合強(qiáng)化模型可以更準(zhǔn)確地模擬阻尼器循環(huán)加載歷程。

3.3 剪切型軟鋼阻尼器薄膜效應(yīng)分析

試驗加載設(shè)備一般使用多連桿機(jī)構(gòu)，將加載設(shè)備的一端固定，另一端僅限制其轉(zhuǎn)動，在實驗加載過程中，作動頭推動阻尼器剪切變形過程中，阻尼器可以產(chǎn)生豎向變形，為純剪切受力模式，腹板不會產(chǎn)生薄膜內(nèi)力。在實際工程應(yīng)用中，阻尼器的上下兩端分別固定于與上下框架梁相連接的支撐或者懸臂墻上，在阻尼器剪切耗能方向，支撐或懸臂墻的抗側(cè)剛度遠(yuǎn)大于阻尼器剪切剛度。地震作用時，結(jié)構(gòu)層間變形大部分作用于阻尼器，阻尼器和子框架變形不協(xié)調(diào)，在阻尼器循環(huán)剪切變形過程中，結(jié)構(gòu)體系會對阻尼器產(chǎn)生附加軸向薄膜內(nèi)力，即薄膜效應(yīng)。

為分析薄膜效應(yīng)對阻尼器耗能性能及疲勞性能的影響，選用MD-S-250-1A 型阻尼器，低屈服點鋼的強(qiáng)化模型選用混合強(qiáng)化模型，下端板完全固定，上部端板保留加載方向的自由度，其余各向自由度全部約束?。ㄜ浖凶杂啥榷x：U1=U2=0，UR1=UR2=UR3=0）。

阻尼器頂部豎向位移約束釋放與頂部豎向位移進(jìn)行約束工況下，剪切力-剪切位移滯回曲線對比見圖11。

圖11 MD-S-250-1A豎向釋放與豎向約束滯回曲線對比

1）約束阻尼器上端豎向位移，在循環(huán)往復(fù)荷載作用下，因為薄膜效應(yīng)，阻尼器滯回曲線出現(xiàn)捏縮，隨著剪切位移的增大，滯回曲線的捏攏現(xiàn)象愈發(fā)明顯。

2）阻尼器在不同的豎向約束條件下，計算得到的滯回曲線剪切力峰值存在顯著差異。當(dāng)剪切位移達(dá)到30 mm 時，釋放豎向約束的剪切力為446.57 kN，約束住豎向位移條件下的剪切力為514.10 kN，兩者剪切力峰值相差15%。但當(dāng)阻尼器經(jīng)歷過一次大的變形后，再進(jìn)行小位移加載時，達(dá)到相同位移條件下，阻尼力會下降，耗能量也相應(yīng)減少，而且前后兩個水平剪切位移的差值越大，阻尼器的阻尼力下降越明顯，阻尼器的耗能性能越不穩(wěn)定。

由圖12 可知，約束住阻尼器上端豎向自由度后，其附加豎向力也隨剪切變形變大而增大后趨于穩(wěn)定。加載初期薄膜內(nèi)力發(fā)展緩慢，隨著加載圈數(shù)的增長，再次加載到相同的位移值，曲線的斜率越來越大，這也表明在低周循環(huán)加載過程中，塑性累積變形會對阻尼器內(nèi)薄膜力的發(fā)展產(chǎn)生影響。

圖12 薄膜內(nèi)力-水平剪切位移

綜合上述分析說明實際工程應(yīng)用中腹板不僅受剪切力還會承受豎向拉壓作用，耗能板內(nèi)的應(yīng)力-應(yīng)變分布也變得更加復(fù)雜，會影響到其耗能性能。

4 剪切型軟鋼阻尼器超低周疲勞壽命模擬

4.1 超低周疲勞分析過程

本文利用非線性有限元分析軟件ABAQUS 和FE-SAFE 進(jìn)行聯(lián)合仿真模擬。利用ABAQUS 對阻尼器擬靜力分析，得到各個節(jié)點單元的應(yīng)力-應(yīng)變數(shù)據(jù)，將其接駁到FE-SAFE 中，對阻尼器疲勞壽命進(jìn)行模擬，分析流程見圖13。

圖13 疲勞分析流程

4.2 模型前處理

1）定義材料種類及力學(xué)參數(shù)，BLY160鋼材的彈性模量E定義為203 000 MPa，泊松比為0.3，屈服強(qiáng)度163 MPa，極限抗拉強(qiáng)度為262 MPa，F(xiàn)E-SAFE 支持自定義材料的E-N壽命曲線，利用改進(jìn)通用斜率法估算出其疲勞參數(shù)見表2。

2）定義荷載譜：疲勞荷載的加載使用前文分析的位移控制工況，采用正負(fù)交替的等幅循環(huán)交變加載譜。

3）設(shè)定阻尼器分析部件均選擇其全部單元參與分析，而不僅僅選擇表面單元。構(gòu)件表面的粗糙度會對構(gòu)件的疲勞壽命產(chǎn)生影響，在模擬中設(shè)定阻尼器板件的表面粗糙度為1.6 μm

4）分析理論：Miner 損傷累積理論，疲勞壽命分析模型：Morrow總應(yīng)變修正模型。Morrow總應(yīng)變修正模型如下：

b為疲勞強(qiáng)度指數(shù)；

c為疲勞延續(xù)指數(shù)。

改進(jìn)通用斜率法引進(jìn)了新的參數(shù)σb/E，通過與材料的彈性模量以及強(qiáng)度極限相關(guān)聯(lián)估算材料疲勞強(qiáng)度系數(shù)和疲勞延續(xù)系數(shù)，改進(jìn)通用斜率法的表達(dá)式為：

4.3 超低周疲勞分析結(jié)果

利用局部應(yīng)力應(yīng)變壽命分析方法獲得MD-S-250-1A和MD-S-250-1B壽命云圖見圖14，兩類阻尼器在試驗約束條件以及工程約束條件下，塑性累計損傷大的地方主要分布在加勁肋板和翼緣板與腹板連接處，與試驗試件破壞主要集中在翼緣根部彎曲變形及加勁肋與腹板連接處焊縫撕裂相對應(yīng)。阻尼器的模擬與試驗試件破壞模式與Abebe［3-4］的研究結(jié)論也是一致的。工程應(yīng)用中上下兩端完全約束，損傷范圍更大且更嚴(yán)重，表明工程應(yīng)用中阻尼器薄膜效應(yīng)影響不可忽視。

圖14 阻尼器疲勞壽命云圖

MD-S-250-1A和MD-S-250-1B兩類阻尼器超低周疲勞模擬壽命統(tǒng)計結(jié)果見表3。

表3 疲勞分析結(jié)果統(tǒng)計

根據(jù)《建筑抗震設(shè)計規(guī)范（GB 50011—2010）》 第12.3.6條規(guī)定，對位移相關(guān)型阻尼器，在設(shè)計位移幅值下往復(fù)循環(huán)30圈后，阻尼器的主要設(shè)計指標(biāo)誤差和衰減量不應(yīng)超過15%，且不應(yīng)有明顯的低周疲勞現(xiàn)象，即30圈為阻尼器疲勞壽命的限值。通過模擬分析表明，此兩類阻尼器在兩種加載工況下，最低疲勞壽命為39 周次，均大于規(guī)范限值30周次的要求，滿足規(guī)范要求。

5 結(jié)論

本文對剪切型軟鋼阻尼器耗能性能和超低周疲勞性能進(jìn)行了相關(guān)研究，對比分析了在不同約束條件下，剪切型軟鋼阻尼器力學(xué)性能與疲勞性能的差異，主要結(jié)論如下：

1）剪切型軟鋼阻尼器耗能板易發(fā)生疲勞破壞的部位主要集中在對角線板角、加勁肋兩側(cè)和翼緣板根部附近位置，與阻尼器試驗的破壞類型數(shù)值模擬的塑性應(yīng)變集中部位相吻合。

2）約束剪切型軟鋼阻尼器豎向位移產(chǎn)生的薄膜效應(yīng)，會降低其疲勞性能，導(dǎo)致其疲勞壽命下降。

3）數(shù)值模擬表明阻尼器在設(shè)計目標(biāo)位移加載下的循環(huán)壽命最低達(dá)到39周次，超過規(guī)范30周次的要求。

4）對剪切型阻尼器的改進(jìn)優(yōu)化時應(yīng)該注意連接焊縫處殘余應(yīng)力的處理以及新的焊接方式的應(yīng)用，優(yōu)先采用應(yīng)力集中系數(shù)小的對接接頭，盡量少采用角焊縫，且焊縫形狀應(yīng)平穩(wěn)過渡。