歐陽(yáng)清容 王雪亮

(1.武漢理工大學(xué)土木工程與建筑學(xué)院，湖北 武漢 430070； 2.武漢理工大學(xué)道路橋梁與結(jié)構(gòu)工程湖北省重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430070)

1 概述

斗拱是木建筑結(jié)構(gòu)中的關(guān)鍵節(jié)點(diǎn)，在地震災(zāi)害中對(duì)結(jié)構(gòu)的整體穩(wěn)定性會(huì)產(chǎn)生嚴(yán)重影響，所以研究斗拱的耗能機(jī)理及特性就顯得尤為重要，也對(duì)建筑結(jié)構(gòu)的整體抗震性提供了參考。國(guó)內(nèi)學(xué)者謝啟芳等[1，2]通過(guò)試驗(yàn)對(duì)比分析殘損與完好狀態(tài)下斗拱節(jié)點(diǎn)的破壞和耗能情況，得出殘損狀態(tài)下耗能能力增強(qiáng)破壞更重，抵抗承載能力更低的結(jié)論。隋等[3]通過(guò)加載試驗(yàn)對(duì)斗拱的側(cè)向剛度退化的情況進(jìn)行了研究，根據(jù)斗拱在試驗(yàn)加載下的滯回曲線(xiàn)分析建立斗拱的恢復(fù)力模型，得出斗拱在地震作用下的變形特點(diǎn)以及耗能特性。高大峰等[4]對(duì)滯回曲線(xiàn)和骨架曲線(xiàn)的分析得到斗拱的位移延性系數(shù)和耗能指標(biāo)參數(shù)。試驗(yàn)研究表明斗拱在不同受荷方向上力學(xué)性能以及耗能特性差異不大，等效粘滯阻尼系數(shù)及抗側(cè)移剛度與加載位移有很大關(guān)系，通過(guò)數(shù)據(jù)分析得出斗拱層的耗能能力較強(qiáng)的結(jié)論。日本學(xué)者津和佑子等[5，6]對(duì)日本建筑的斗拱為原型的力學(xué)試驗(yàn)得出斗拱的剛度、變形特征以及斗拱下部櫨斗的力學(xué)性能。上述研究對(duì)斗拱的力學(xué)性能和耗能能力進(jìn)行了探討，但是對(duì)于斗拱的低周反復(fù)荷載試驗(yàn)的有限元分析以及斗拱在循環(huán)往復(fù)荷載作用下的耗能特性和剛度退化方面的研究很少。本文研究了仿古木塔典型斗拱的節(jié)點(diǎn)滑移前后的剛度變化和耗能特征系數(shù)并與趙鴻鐵教授的斗拱試驗(yàn)結(jié)果對(duì)比，再根據(jù)單朵斗拱的剛度值和粘滯阻尼系數(shù)推出整個(gè)鋪?zhàn)鲗拥膭偠戎岛秃哪芴卣髦?，為木塔的整體分析奠定了基礎(chǔ)。

2 斗拱耗能有限元分析

2.1 斗拱有限元模型

2.1.1斗拱設(shè)計(jì)

本文選取三層仿古木塔模型中的補(bǔ)間鋪?zhàn)鳛檠芯繉?duì)象，補(bǔ)間五鋪?zhàn)鞫饭俺叽缫?guī)格如表1所示。

表1 補(bǔ)間五鋪?zhàn)鳂?gòu)件尺寸規(guī)格表

2.1.2模型的建立和加載

按照補(bǔ)間五鋪?zhàn)鞫饭暗某叽缫?guī)格和參考《營(yíng)造法式》的相關(guān)做法，使用ABAQUS有限元軟件建立斗拱的實(shí)體模型，如圖1所示。為使荷載能夠更加均勻地施加在斗拱構(gòu)件上，在建立有限元模型時(shí)使斗拱的頂部加裝干枋木，并且使構(gòu)件之間通過(guò)榫卯連結(jié)，達(dá)到增強(qiáng)構(gòu)件之間的咬合力。在LOAD模塊上荷載工況是通過(guò)定義變幅的水平簡(jiǎn)諧荷載達(dá)到，幅值從10 kN加載到40 kN，每個(gè)間隔為5 kN。除此之外在鋪?zhàn)鞯捻敳渴┘?0 kN的豎向荷載以定義頂部枋木之間的耦合約束。為了避免在水平荷載加載時(shí)斗拱發(fā)生傾覆現(xiàn)象，斗拱的邊界條件是對(duì)普拍枋底部固結(jié)，枋木上部施加豎向荷載限制其活動(dòng)。

2.1.3網(wǎng)格劃分

Mesh模塊是對(duì)有限元模型進(jìn)行網(wǎng)格的劃分，網(wǎng)格劃分的精細(xì)程度決定運(yùn)算結(jié)果的準(zhǔn)確度，如果劃分太密集就會(huì)延長(zhǎng)計(jì)算時(shí)間從而影響穩(wěn)定極限，所以網(wǎng)格的劃分也是極其重要的一步。此模型斗拱的網(wǎng)格尺寸大概在50 mm，某些重要構(gòu)件或部位的網(wǎng)格會(huì)細(xì)化，比如受力構(gòu)件和構(gòu)件間的接觸區(qū)域。在ABAQUS/Standard單元庫(kù)中選擇C3D8I實(shí)體單元(8結(jié)點(diǎn)線(xiàn)性六面體單元，非協(xié)調(diào)模式)來(lái)生成六面體單元，另外選擇C3D10M實(shí)體單元(10結(jié)點(diǎn)修正二次四面體單元)來(lái)生成四面體單元，共得到1 396個(gè)實(shí)體單元。斗拱劃分的網(wǎng)格圖如圖2所示。

2.1.4傳力路徑

傳力路徑選擇水平簡(jiǎn)諧荷載為30 kN時(shí)作為斗拱的整體傳力路徑，如圖3所示。

由圖3可以看出，水平向簡(jiǎn)諧荷載主要沿著華拱方向傳遞，而橫拱方向傳遞的荷載較少。由于斗拱頂部加裝若干枋木，使得水平向的作用力主要由各華拱來(lái)承擔(dān)，首先由長(zhǎng)華拱承擔(dān)作用力，后由散斗的暗梢使作用力向短華拱傳遞，再傳至下部的櫨斗，并最后由櫨斗下部的暗梢傳遞至普拍枋。當(dāng)水平向的荷載繼續(xù)增加時(shí)，最易發(fā)生破壞的就是這些主要承力構(gòu)件，如散斗、長(zhǎng)短華拱以及下部的櫨斗等部件。

2.2 簡(jiǎn)諧荷載作用下斗拱的耗能減震機(jī)理分析

2.2.1力—位移曲線(xiàn)

對(duì)斗拱施加水平向的變幅簡(jiǎn)諧荷載，其力—位移滯回曲線(xiàn)如圖4所示。

從圖4可以看出，當(dāng)簡(jiǎn)諧荷載幅值較小時(shí)，力—位移滯回曲線(xiàn)呈梭型，櫨斗與下部普拍枋基本未發(fā)生滑移；當(dāng)簡(jiǎn)諧荷載幅值大于25 kN時(shí)，力—位移滯回曲線(xiàn)呈斜四邊形狀，在力達(dá)到幅值時(shí)，可以發(fā)現(xiàn)有明顯的滑移。從圖4中可知，當(dāng)簡(jiǎn)諧荷載較小時(shí)，力—位移曲線(xiàn)呈線(xiàn)性相關(guān)；當(dāng)水平簡(jiǎn)諧力大于22.5 kN(即層間位移為4.85 mm)時(shí)，斗拱下部櫨斗與普拍枋的靜摩擦力小于水平作用力，下部櫨斗開(kāi)始相對(duì)普拍枋發(fā)生相對(duì)滑移，滯回曲線(xiàn)開(kāi)始發(fā)生非線(xiàn)性變化，斗拱剛度下降，滯回曲線(xiàn)面積開(kāi)始擴(kuò)大，此時(shí)斗拱的耗能增加。另外對(duì)比各幅值簡(jiǎn)諧波荷載的力—位移滯回曲線(xiàn)的面積，簡(jiǎn)諧荷載幅值越大，滯回環(huán)面積越大，說(shuō)明當(dāng)輸入較大外力荷載時(shí)，斗拱耗散的能量也會(huì)隨之增加。

2.2.2豎向剛度

對(duì)斗拱模型頂部的枋木施加豎直向下的均布荷載，根據(jù)有限元分析當(dāng)豎直荷載值大小為715 kN時(shí)斗拱開(kāi)始屈服，可將715 kN視為斗拱的極限承載力值，加載過(guò)程中荷載與位移的關(guān)系圖如圖5所示。

從圖5可知，在進(jìn)行加載過(guò)程中，斜率在3 mm左右開(kāi)始增大，

并保持斜率不變直到斗拱進(jìn)入塑性狀態(tài)，這是因?yàn)槎饭笆怯稍S多小構(gòu)件組成的，在裝配過(guò)程中難免會(huì)有縫隙存在，在荷載作用下縫隙會(huì)不斷擠緊從而使剛度增加。由于斗拱極限承載力值為715 kN，參考文獻(xiàn)[7]中的方法，在荷載位移曲線(xiàn)中將10%和40%的極限承載力值兩點(diǎn)的割線(xiàn)斜率作為斗拱的豎向剛度值，求解剛度值為52.98 kN/mm。

2.3 抗側(cè)移剛度和阻尼系數(shù)

根據(jù)圖5的荷載—位移滯回曲線(xiàn)，通過(guò)計(jì)算得到斗拱的抗側(cè)移剛度值和阻尼系數(shù)值?？箓?cè)移剛度大小為3 561.2 N/mm，滑移后斗拱的抗側(cè)移剛度為445.15 N/mm；粘滯阻尼系數(shù)取0.391。

3 結(jié)語(yǔ)

1)對(duì)斗拱在簡(jiǎn)諧荷載下的傳力路徑進(jìn)行研究，并得出當(dāng)荷載值繼續(xù)增大時(shí)，斗拱的部分構(gòu)件會(huì)容易發(fā)生破壞。

2)當(dāng)簡(jiǎn)諧荷載較小時(shí)，力—位移曲線(xiàn)呈線(xiàn)性相關(guān)；當(dāng)水平簡(jiǎn)諧力大于22.5 kN(即層間位移為4.85 mm)時(shí)，斗拱下部櫨斗與普拍枋的靜摩擦力小于水平作用力，下部櫨斗開(kāi)始相對(duì)普拍枋發(fā)生相對(duì)滑移，滯回曲線(xiàn)開(kāi)始發(fā)生非線(xiàn)性變化，斗拱剛度下降，滯回曲線(xiàn)面積開(kāi)始擴(kuò)大，此時(shí)斗拱的耗能增加。另外斗拱具有遇強(qiáng)則強(qiáng)的特性，當(dāng)輸入較大外力荷載時(shí)，斗拱耗散的能量也會(huì)隨之增加。通過(guò)分析所得的滯回曲線(xiàn)可知，模擬分析的結(jié)果和趙鴻鐵教授[8]的低周反復(fù)荷載試驗(yàn)結(jié)果類(lèi)似，說(shuō)明通過(guò)簡(jiǎn)諧荷載對(duì)斗拱進(jìn)行有限元分析的計(jì)算結(jié)果符合要求，分析結(jié)果可以作為整體模型分析的參考。

[1] 謝啟芳.殘損古建筑木結(jié)構(gòu)叉柱造式斗栱節(jié)點(diǎn)抗震性能退化規(guī)律研究[J].土木工程學(xué)報(bào)，2014，47(12)：49-51.

[2] 謝啟芳.古建筑木結(jié)構(gòu)叉柱造式斗栱節(jié)點(diǎn)抗震性能試驗(yàn)研究[J].土木工程學(xué)報(bào)，2015，48(8)：20-24.

[4] 高大峰.木結(jié)構(gòu)古建筑斗栱結(jié)構(gòu)層抗震性能試驗(yàn)研究[J].地震工程與工程振動(dòng)，2014，34(1)：132-135.

[5] Hideo KYUKE.Shaking Table Test of“MASUGUMI”Used in Traditional Wooden Architectures.WCTE2008-10th World Conference on Timber Engineering.JAPAN,2008:226.

[6] Iuko Tsuwa.A Study on the Size Effect of Bracket Complexes Used in Traditional Timber Structures on the Vibration Characteristics.WCTE2008-10th World Conference on Timber Engineering,JAPAN,2008:230.

[7] 李合群.中國(guó)傳統(tǒng)建筑構(gòu)造[M].北京：北京大學(xué)出版社，2010：146-148.

[8] 趙鴻鐵.中國(guó)古建筑結(jié)構(gòu)及其抗震：試驗(yàn)、理論及加固方法[M].北京：科學(xué)出版社，2012：50-53.