巴振寧 焦鵬 劉冰松 張翔

摘要：以嘉峪關(guān)光化樓為研究對(duì)象，基于榫卯、斗拱節(jié)點(diǎn)、柱墻接觸之間的關(guān)系，利用ANSYS軟件分別建立了純木結(jié)構(gòu)與磚木混合結(jié)構(gòu)的光化樓模型。通過(guò)模態(tài)分析和輸入多種地震波進(jìn)行動(dòng)力時(shí)程分析，結(jié)果表明：①磚墻結(jié)構(gòu)大大提高了光化樓低階自振頻率，在高階振型中自振頻率的差距逐漸減小;②磚墻結(jié)構(gòu)使結(jié)構(gòu)在地震時(shí)程分析下的加速度峰值有了明顯增大，這與結(jié)構(gòu)自振周期降低有關(guān);③在地震作用下，磚墻結(jié)構(gòu)使結(jié)構(gòu)的位移峰值顯著減小，對(duì)防止木構(gòu)架側(cè)移倒塌有很大貢獻(xiàn);④純木結(jié)構(gòu)模型位移反應(yīng)較大，但仍小于木結(jié)構(gòu)建筑彈塑性層間位移角1/50的限值;⑤在Ⅷ度罕遇地震作用下，磚墻結(jié)構(gòu)底部及開(kāi)洞處的塑性應(yīng)變較大，具體可能表現(xiàn)為磚墻底部以及磚墻與木門(mén)連接處出現(xiàn)裂縫。

關(guān)鍵詞：光化樓;純木結(jié)構(gòu);磚木混合結(jié)構(gòu);模態(tài)分析;時(shí)程分析

0 引言

中國(guó)木結(jié)構(gòu)建筑文化源遠(yuǎn)流長(zhǎng)，留有許多寶貴的文化結(jié)晶。嘉峪關(guān)建筑群始建于公元1372年，建成后經(jīng)歷了多次的損毀與修繕。1987年，嘉峪關(guān)建筑群被聯(lián)合國(guó)教科文組織列入《世界文化遺產(chǎn)名錄》。我國(guó)是一個(gè)地震多發(fā)的國(guó)家，嘉峪關(guān)位于世界上最大的地震帶“喜馬拉雅—環(huán)太平洋地震帶”的作用區(qū)。因此對(duì)嘉峪關(guān)木構(gòu)城樓進(jìn)行抗震研究和保護(hù)，具有十分重要的意義。

國(guó)內(nèi)外一些建筑大師，如梁思成（1983）、伊東忠太和陳清泉（2014）對(duì)我國(guó)的古建筑進(jìn)行了系統(tǒng)的研究。隨著力學(xué)學(xué)科體系的發(fā)展，對(duì)結(jié)構(gòu)的整體和不同構(gòu)件的研究也有了迅速發(fā)展。周乾等（2014）建立故宮太和殿有限元模型響應(yīng)譜，分析了在Ⅷ度罕遇地震作用下，太和殿結(jié)構(gòu)的內(nèi)力及變形;謝啟芳等（2018）以西安鐘樓為研究對(duì)象，建立了縮尺比例為1：6的模型進(jìn)行振動(dòng)臺(tái)實(shí)驗(yàn)，得到了在不同幅值、不同地震波作用下模型的加速度、位移和破壞形態(tài);潘毅等（2019）以饒益寺大佛殿和保國(guó)寺的鋪?zhàn)鲗訛閰⒄?，設(shè)計(jì)了4種不同的鋪?zhàn)鲗硬贾媚Ｐ?，進(jìn)行動(dòng)力時(shí)程分析，研究不同鋪?zhàn)鲗硬贾脤?duì)地震激勵(lì)作用下結(jié)構(gòu)產(chǎn)生的動(dòng)力響應(yīng)的影響;Uchida等（2004）、Kawai等（2006）和Fujita等（2008）分別對(duì)日本傳統(tǒng)的多層木塔與廟宇結(jié)構(gòu)進(jìn)行了側(cè)向靜力試驗(yàn)，結(jié)合日本地區(qū)地震記錄數(shù)據(jù)，建立層模型并分析在地震作用下各節(jié)點(diǎn)與結(jié)構(gòu)整體的動(dòng)力響應(yīng)特點(diǎn)。

上述研究分別對(duì)不同的木結(jié)構(gòu)建筑進(jìn)行了研究分析，但大多集中于受力結(jié)構(gòu)，而忽略了圍護(hù)結(jié)構(gòu)。木結(jié)構(gòu)古建筑的形制復(fù)雜，構(gòu)件的應(yīng)用組合多樣，再加上木材的特性，不同建筑都有其獨(dú)特的結(jié)構(gòu)狀態(tài)，值得研究。嘉峪關(guān)木構(gòu)城樓第一層為磚木混合結(jié)構(gòu)，第二、三層為木結(jié)構(gòu)。一層金柱嵌固在磚砌體墻內(nèi)，因此磚墻對(duì)木結(jié)構(gòu)力學(xué)性能勢(shì)必有所影響。本文采用有限元數(shù)值模擬方法，建立嘉峪關(guān)光化樓木結(jié)構(gòu)有限元分析模型與磚木混合結(jié)構(gòu)有限元分析模型，分別輸入Ⅷ度多遇、設(shè)防與罕遇的地震波激勵(lì)，分析磚墻結(jié)構(gòu)對(duì)抗震性能的影響。

1 數(shù)值模型

1.1 光化樓簡(jiǎn)介

光化樓始建于公元1506年，位于嘉峪關(guān)東城臺(tái)之上，城臺(tái)平面為矩形，臺(tái)面采用條磚十字縫墁地，城臺(tái)下有拱券式門(mén)洞。嘉峪關(guān)建筑群如圖1所示。

光化樓為典型的抬梁木構(gòu)城樓，屋頂樣式采用單檐歇山頂式。由圖2可見(jiàn)，一層柱分內(nèi)外兩周，共計(jì)28根（檐柱18根、金柱10根），為貫穿三層的通柱，層高約4.4 m;二層有外廊，柱分內(nèi)外兩周，廊柱18根，層高約4 m;三層無(wú)外廊，共計(jì)檐柱10根，層高約3.7 m。建筑一層為磚木結(jié)構(gòu)，采用磚墻，二、三層為木結(jié)構(gòu)，采用槅扇，磚墻包裹金柱，厚約1.3 m，周長(zhǎng)約41.5 m，東西兩側(cè)開(kāi)門(mén)，門(mén)高2.55 m、寬2.12 m。

1.2 單元選擇

光化樓模型采用ANSYS有限元分析軟件進(jìn)行模擬。對(duì)于傳統(tǒng)木結(jié)構(gòu)古建筑，結(jié)構(gòu)主體部分主要由梁柱構(gòu)架組成，采用Beam188 單元對(duì)其進(jìn)行數(shù)值模擬。Beam 188單元在每個(gè)節(jié)點(diǎn)上有6個(gè)自由度，包括沿X，Y，Z軸的線位移和繞X，Y，Z軸的角位移，適合分析細(xì)長(zhǎng)到中等長(zhǎng)度的梁結(jié)構(gòu)。

對(duì)于屋蓋、樓板和磚墻結(jié)構(gòu)，采用Shell181單元進(jìn)行建模。該單元內(nèi)各節(jié)點(diǎn)含有6個(gè)自由度，常用于模擬平板和曲殼等形式的結(jié)構(gòu)，能夠傳遞單元彎矩并在非線性計(jì)算中考慮殼厚度的變化。

對(duì)于木結(jié)構(gòu)榫卯、斗拱節(jié)點(diǎn)，采用Combin14彈簧單元進(jìn)行模擬，該彈簧單元可以傳遞軸力、剪力和彎矩，通過(guò)2個(gè)節(jié)點(diǎn)來(lái)定義，可以設(shè)置6個(gè)方向的彈簧剛度。

1.3 材料參數(shù)及邊界條件設(shè)置

對(duì)木結(jié)構(gòu)進(jìn)行結(jié)構(gòu)分析驗(yàn)算時(shí)，材料性能強(qiáng)度不僅要滿(mǎn)足《木結(jié)構(gòu)設(shè)計(jì)規(guī)范》（GB 50005—2017）的規(guī)定，還應(yīng)該考慮荷載長(zhǎng)期作用和木質(zhì)老化的影響，對(duì)木材的各項(xiàng)力學(xué)性能指標(biāo)給予調(diào)整。本文將木材本構(gòu)關(guān)系簡(jiǎn)化為理想彈塑性，遵循Hill屈服準(zhǔn)則。

由于對(duì)嘉峪關(guān)光化樓缺乏直接的實(shí)驗(yàn)數(shù)據(jù)，參照祝松濤（2013）對(duì)嘉峪關(guān)關(guān)樓的研究，其木材為紅松，為正交異性材料，木材參數(shù)見(jiàn)表1。將紅松的各個(gè)方向材料性能簡(jiǎn)化為理想彈塑性，此處假定木材的受拉與受壓方向的屈服強(qiáng)度相同，采用Hill準(zhǔn)則與BISO模型結(jié)合的方式進(jìn)行定義，即通過(guò)Hill準(zhǔn)則定義6個(gè)方向不同的屈服強(qiáng)度，采用BISO模型進(jìn)行切線模量的定義，此處切線模量定義為0，各個(gè)方向屈服強(qiáng)度值參考王玉迪（2018）對(duì)紅松木材的研究，見(jiàn)表2，理想彈塑性模型如圖3所示。其中考慮到光化樓修建完成距今已經(jīng)超過(guò)500年，根據(jù)《古建筑木結(jié)構(gòu)維護(hù)與加固技術(shù)規(guī)范》（GB 50165—92）中對(duì)考慮長(zhǎng)期荷載作用和木質(zhì)老化調(diào)整系數(shù)的規(guī)定，對(duì)強(qiáng)度等性能均乘以0.7的折減系數(shù)。

朱才輝等（2019）曾對(duì)同為明代建筑的西安明代城墻以及鐘樓的磚墻砌體進(jìn)行過(guò)研究。結(jié)合其得出的磚墻砌體等效力學(xué)參數(shù)，并參考《砌體結(jié)構(gòu)設(shè)計(jì)規(guī)范》（GB 50003—2011），本文采用的磚砌體墻材料力學(xué)參數(shù)見(jiàn)表3。

式中：σD，εD為考慮損傷因子的應(yīng)力與應(yīng)變;損傷因子D為0～1，數(shù)值越大表示結(jié)構(gòu)的損傷程度越大，光化樓磚墻結(jié)構(gòu)無(wú)明顯破壞且本文主要考慮磚墻結(jié)構(gòu)對(duì)抗震性能的影響，因此D取0;ε0為對(duì)應(yīng)于砌體抗壓平均值的應(yīng)變，根據(jù)朱才輝等（2019）的研究，取0.003。砌體結(jié)構(gòu)的應(yīng)力-應(yīng)變曲線如圖4所示，采用MISO模型進(jìn)行定義。

斗拱是古建筑中主要的結(jié)構(gòu)連接形式，由多層木塊交疊、咬合而成，其動(dòng)力學(xué)特性十分復(fù)雜，在地震等外力荷載作用下，斗拱可以視為一組彈簧，通過(guò)自身擠壓變形，減緩結(jié)構(gòu)的變形，并且吸收掉一部分振動(dòng)能量。本文采用Combin14彈簧來(lái)模擬斗拱半剛性的特性。梁柱榫卯節(jié)點(diǎn)與斗拱，均具有半剛性連接特點(diǎn)，同樣簡(jiǎn)化為Combin14單元。斗拱、榫卯各個(gè)方向自由度的剛度值參考高大峰等（2008）、蘇軍（2008）的試驗(yàn)以及數(shù)值模擬分析結(jié)果，節(jié)點(diǎn)剛度取值見(jiàn)表4。

中國(guó)古建筑結(jié)構(gòu)普遍通過(guò)柱礎(chǔ)與地基相連接，柱根一般擱置在礎(chǔ)石之上，在主體結(jié)構(gòu)的重力作用下，礎(chǔ)石能夠?yàn)橹_提供豎直方向的支撐力和

水平方向上一定的摩擦力，對(duì)柱根的轉(zhuǎn)動(dòng)具有微弱的限制。參考高大峰等（2008）、蘇軍（2008）中對(duì)木結(jié)構(gòu)古建筑的處理方法，本文將柱腳與礎(chǔ)石的連接節(jié)點(diǎn)近似簡(jiǎn)化為鉸接節(jié)點(diǎn);光化樓一層磚墻將金柱包裹在內(nèi)，牢牢嵌固，將墻柱連接設(shè)置為剛接。

1.4 建立數(shù)值模型

光化樓的梁柱框架組成了主要的結(jié)構(gòu)受力體系，主體結(jié)構(gòu)可以劃分為屋架層、鋪?zhàn)鲗右约爸W(wǎng)層。屋架層主要由梁、枋、檁條等組成的梁架體系來(lái)支撐屋頂以及傳遞上部荷載;鋪?zhàn)鲗又饕远饭盀橹?，?duì)各層起承上啟下的作用;柱網(wǎng)層是光化樓結(jié)構(gòu)的豎向支撐部分，將整棟樓所承擔(dān)的荷載傳導(dǎo)到臺(tái)基上。對(duì)主體結(jié)構(gòu)逐層建立無(wú)磚墻結(jié)構(gòu)數(shù)值模型（以下簡(jiǎn)稱(chēng)模型1），如圖5a所示。

磚砌體墻是由砌塊混合砂漿等粘合劑制作成的受力結(jié)構(gòu)，主要受到砌塊材料性能的影響，但由于砂漿等粘合劑的材料特性與磚砌體存在顯著差異，所以砂漿等粘合劑對(duì)結(jié)構(gòu)整體的影響也不可忽視。且結(jié)構(gòu)的整體力學(xué)性能還受到諸如磚砌體尺寸、砌縫的寬度以及磚墻的砌筑方式、工藝水平等因素的影響。

目前，常用的磚砌體墻的有限元建模方法包含分離式和整體式2種。分離式建模方法就是對(duì)磚砌體墻的2種材料分別進(jìn)行建模，并對(duì)2種材料接

觸面做出定義以模擬其工作性能。整體式建模方法就是將磚砌體墻視為均質(zhì)的各向同性材料，忽略結(jié)構(gòu)模型中局部的材質(zhì)差異。采取整體式建模方法時(shí)，可以適當(dāng)簡(jiǎn)化磚墻結(jié)構(gòu)的材料本構(gòu)與力學(xué)參數(shù)，大大簡(jiǎn)化了建模難度和運(yùn)算過(guò)程，適合于將磚砌體墻的看作宏觀整體的有限元結(jié)構(gòu)分析，本文采用該方法將磚墻結(jié)構(gòu)按照?qǐng)D紙進(jìn)行建模，其中東西向磚墻按照門(mén)尺寸進(jìn)行開(kāi)洞，尺寸為2.12 m×2.55 m，最終建立有磚墻結(jié)構(gòu)數(shù)值模型（簡(jiǎn)稱(chēng)模型2），如圖5b所示。

2 模態(tài)分析

模態(tài)分析是研究結(jié)構(gòu)系統(tǒng)特性的一種方法，主要用于求解結(jié)構(gòu)的自振頻率、模態(tài)振型等固有振動(dòng)特性。對(duì)于古建筑來(lái)講，對(duì)其地震荷載作用下振型的分析，是結(jié)構(gòu)動(dòng)力學(xué)分析的基礎(chǔ)。通常情況下高階模態(tài)的權(quán)重遠(yuǎn)低于低階模態(tài)，一般取前幾階模態(tài)進(jìn)行分析。

本文采用子空間法提取數(shù)值模型模態(tài)，子空間法在提取中大型模型且小于40階振型時(shí)具有較高的精度。數(shù)值模型振型頻率值變化如圖6所示。由圖可知，嘉峪關(guān)光化樓數(shù)值模型前20階的頻率均分布至0～10 Hz，自振頻率隨階數(shù)增大逐漸增大，在高階時(shí)增長(zhǎng)逐漸變緩。

磚墻結(jié)構(gòu)的存在使結(jié)構(gòu)低階自振周期顯著減低，在高階振型中自振頻率逐漸趨于一致，對(duì)結(jié)構(gòu)整體剛度有較大的貢獻(xiàn)，對(duì)光化樓整體的動(dòng)力性能有顯著影響。

3 抗震性能分析

3.1 地震波的選取

根據(jù)《建筑抗震設(shè)計(jì)規(guī)范》（GB 50011—2010），嘉峪關(guān)處于Ⅷ度抗震設(shè)防區(qū)，基本加速度為0.2 g，進(jìn)行時(shí)程分析時(shí)根據(jù)多遇地震、設(shè)防地震、罕遇地震分別調(diào)幅為70 gal，200 gal，400 gal 3種不同加速度峰值的地震波。選取的地震波分別為El-Centro波、Taft波與蘭州波。以200 gal峰值地震波為例，其地震影響系數(shù)如圖7所示。

3.2 阻尼確定

阻尼作為一種能量耗散機(jī)制，使振動(dòng)等運(yùn)動(dòng)隨時(shí)間減弱并最終停止。采用瑞利阻尼進(jìn)行設(shè)置，其阻尼矩陣C由質(zhì)量矩陣M和剛度矩陣K表示為：

根據(jù)俞茂宏等（1991）、Fang等（2001a，b）對(duì)木結(jié)構(gòu)整體動(dòng)力試驗(yàn)測(cè)試與現(xiàn)場(chǎng)環(huán)境脈動(dòng)試驗(yàn)，測(cè)定其阻尼比為2%～9%，一般小于5%，本文阻尼比取3.5%。根據(jù)之前模態(tài)分析結(jié)果，模型1：α=0.241 3，β=0.003 336;模型2：α=0.498 1，β=0.001 605。

3.3 加速度反應(yīng)分析

對(duì)模型輸入水平雙向地震波，由于Y向建筑尺寸小于X向，故輸入地震波時(shí)Y向與X向峰值加速度比為1∶ 0.85，采用前文所述9條地震動(dòng)記錄進(jìn)行輸入，模型加載工況見(jiàn)表5。

為分析不同地震加速度峰值下光化樓地震響應(yīng)，分別選擇一層檐枋、二層檐枋、三層柱頂位置，提取其動(dòng)力響應(yīng)數(shù)據(jù)。所有工況作用下的測(cè)點(diǎn)的振動(dòng)加速度峰值如圖8所示，所有工況作用下的放大系數(shù)平均值見(jiàn)表6，7。

由圖8，表6可知：在水平地震作用下，隨地震波激勵(lì)的增加，各層的加速度峰值也隨之增大。模型1（M1）加速度響應(yīng)放大均小于1，模型2（M2）加速度響應(yīng)值均有所提高，尤其是在Taft波作用下，加速度響應(yīng)峰值增大1倍左右，這與磚墻結(jié)構(gòu)降低了結(jié)構(gòu)的自振周期，使其接近Taft波的卓越周期有關(guān)。

3.4 位移反應(yīng)分析

選擇一層檐枋、二層檐枋、三層柱頂位置，提取其位移響應(yīng)數(shù)據(jù)。由于文章篇幅限制，僅展示模型1與模型2三層柱頂在L1，L2，L3三種工況作用下的位移時(shí)程曲線如圖9所示，所有工況位移峰值見(jiàn)表7。

由圖9，表7可知，在水平地震作用下，考慮磚墻結(jié)構(gòu)的情況下，各層位移的相應(yīng)峰值均減小，這反映了磚墻結(jié)構(gòu)對(duì)木構(gòu)架的位移有很強(qiáng)的約束作用，對(duì)防止木構(gòu)架側(cè)移倒塌有很大貢獻(xiàn)。模型1在峰值為400 gal的El-Centro波作用下短軸向和長(zhǎng)軸向出現(xiàn)最大層架位移角，分別為1/68和1/111，小于木結(jié)構(gòu)建筑彈塑性層間位移角1/50的限值。

3.5 墻體塑性應(yīng)變

光化樓開(kāi)洞處橫墻、未開(kāi)洞處橫墻以及縱墻在400 gal的El-Centro波作用下產(chǎn)生的等效塑性應(yīng)變情況如圖10所示。由圖可知，未開(kāi)洞墻體底部的塑性應(yīng)變明顯高于上部，開(kāi)洞墻體的塑性應(yīng)變?cè)诘撞颗c開(kāi)洞部位更為明顯。這表明墻體在Ⅷ度罕遇地震波作用下，墻體底部以及墻體與木門(mén)連接處容易出現(xiàn)裂縫等震害。

4 結(jié)論

本文以嘉峪關(guān)光化樓為原型，根據(jù)其一樓為磚木混合結(jié)構(gòu)的特點(diǎn)，建立了有限元數(shù)值模型，通過(guò)動(dòng)力時(shí)程分析，討論了磚墻結(jié)構(gòu)對(duì)光化樓抗震性能的影響，得到以下結(jié)論：

（1）金柱與磚墻結(jié)構(gòu)牢牢嵌固在一起，形成環(huán)箍效應(yīng)，大大高了結(jié)構(gòu)的抗側(cè)剛度。在模態(tài)分析時(shí)，磚墻結(jié)構(gòu)的存在對(duì)結(jié)構(gòu)的低階振型影響較大，光化樓結(jié)構(gòu)的低階自振頻率顯著提高，高階自振頻率與無(wú)磚墻結(jié)構(gòu)自振頻率隨振型階數(shù)的增加逐漸相等。

（2）磚墻結(jié)構(gòu)的存在使建筑物加速度響應(yīng)峰值有所提高，這與建筑物自振周期的變小有關(guān)。

（3）磚墻結(jié)構(gòu)使建筑物位移峰值顯著降低，有助于防止木結(jié)構(gòu)位移過(guò)大產(chǎn)生側(cè)移傾覆。

（4）對(duì)層間位移角進(jìn)行了計(jì)算，嘉峪關(guān)光化樓在不同地震波作用下產(chǎn)生的最大層間位移角均符合相關(guān)規(guī)范的限值規(guī)定。

（5）在Ⅷ度罕遇地震作用下，墻體的底部與開(kāi)洞處的塑性變形較大，墻體底部以及與木門(mén)連接處易產(chǎn)生裂縫。

參考文獻(xiàn)：

高大峰，趙鴻鐵，薛建陽(yáng).2008.木結(jié)構(gòu)古建筑中斗栱與榫卯節(jié)點(diǎn)的抗震性能—試驗(yàn)研究[J].自然災(zāi)害學(xué)報(bào)，17（2）：58-64.

梁思成.1983.營(yíng)造法式注釋[M].北京：中國(guó)建筑工業(yè)出版社.

劉桂秋.2005.砌體結(jié)構(gòu)基本受力性能的研究[D].長(zhǎng)沙：湖南大學(xué).

潘毅，袁雙，郭瑞，等.2019.鋪?zhàn)鲗硬贾脤?duì)古建筑木結(jié)構(gòu)抗震性能的影響[J].土木工程學(xué)報(bào)，52（3）：29-40.

蘇軍.2008.中國(guó)木結(jié)構(gòu)古建筑抗震性能的研究[D].西安：西安建筑科技大學(xué).

王玉迪.2018.考慮松動(dòng)殘損直榫節(jié)點(diǎn)的多層木結(jié)構(gòu)古建筑抗震性能分析[D].西安：西安建筑科技大學(xué).

謝啟芳，王龍，張利朋，等.2018.西安鐘樓木結(jié)構(gòu)模型振動(dòng)臺(tái)試驗(yàn)研究[J].建筑結(jié)構(gòu)學(xué)報(bào)，39（12）：128-138.

伊東忠太，陳清泉.2014.中國(guó)建筑史[M].長(zhǎng)沙：湖南大學(xué)出版社.

俞茂宏，劉曉東，方東平，等.1991.西安北門(mén)箭樓靜力與動(dòng)力特性的試驗(yàn)研究[J].西安交通大學(xué)學(xué)報(bào)，25（3）：55-62.

周乾，閆維明，關(guān)宏志，等.2014.罕遇地震作用下故宮太和殿抗震性能研究[J].建筑結(jié)構(gòu)學(xué)報(bào)，35（S1）：25-32.

朱才輝，劉欽佩，周遠(yuǎn)強(qiáng).2019.古建筑磚-土結(jié)構(gòu)力學(xué)性能及裂縫成因分析[J].建筑結(jié)構(gòu)學(xué)報(bào)，40（9）：157-186.

祝松濤.2013.嘉峪關(guān)關(guān)樓抗震性能分析及研究保護(hù)[D].西安：西安建筑科技大學(xué).