考慮臺基影響的應縣木塔地震響應有限元分析

薛建陽 吳晨偉 翟磊 王瑞鵬 馬林林

摘 要：為研究臺基對應縣木塔地震響應的影響，采用ANSYS有限元軟件建立了八角形臺基、上部木結構及考慮臺基影響的應縣木塔整體結構有限元模型，采用模態分析法得到了模型的自振頻率及周期。分析得到了上部木結構的弱水平剛度方向，并對上部木結構和考慮臺基影響的整體結構模型輸入地震波激勵，對比分析了兩種模型的地震響應。結果表明：八角形臺基的頻率較高且各階頻率比較接近，當上部木結構的高階頻率與臺基頻率接近時，臺基對上部木結構動力特性有較大影響。木塔東西方向為弱水平剛度方向。相比其他地震波激勵，當長周期Tianjin波激勵時，上部木結構和考慮臺基影響的整體結構模型的加速度和位移響應更大。隨著輸入激勵的峰值加速度增大，模型結構動力系數呈減小趨勢。在相同峰值加速度的同種地震波激勵下，考慮臺基影響的整體結構加速度和位移響應均大于僅考慮上部木結構的情況。

關鍵詞：應縣木塔;臺基;動力特性;地震響應;有限元分析

中圖分類號：TU366.2 文獻標志碼：A 文章編號：2096-6717（2022）02-0022-08

基金項目：國家自然科學基金（51978568）;上海市工程結構安全重點實驗室開放課題（2020-KF07）;陜西省科技創新團隊項目（2019TD-029）;陜西省自然科學基礎研究計劃（2020JZ-50）

作者簡介：薛建陽（1970- ），男，博士，教授，主要從事木結構抗震及鋼混凝土組合結構研究，E-mail：jianyang_xue@163.com。

Abstract： In order to study the effect of stylobate on the seismic responses of Yingxian Wooden Tower ，the models of octagonal stylobate， the upper timber frame and the overall structure are built by ANSYS finite element software. The modal analysis method is used to obtain the natural frequencies and periods of the model. The direction of the weak horizontal stiffness of the upper timber structure is obtained， and seismic wave excitation is input to the upper wooden structure and the overall structure model respectively. Then the seismic responses of the two models are compared and analyzed. Results show that the octagonal stylobate has a higher natural frequency and all frequencies are relatively close. And when the high-order frequency of the upper timber structure is close to that of the stylobate， the stylobate obviously affects the dynamic properties of the upper timber structure. The east-west direction of the wooden tower is the direction of weak horizontal stiffness.Compared with other seismic wave excitations， the acceleration and displacement responses of the upper timber structure and the overall structure model are greater when the long-period Tianjin wave is input. As the peak acceleration of the inputting excitation increases， the dynamic coefficients only considering the structure decrease.Under the excitation of the same type of seismic wave with the same peak acceleration， the acceleration and displacement response of the overall structure are greater than those of the upper timber structure.

Keywords：Yingxian Wooden Tower; stylobate; dynamic properties; seismic response; finite element analysis

應縣木塔全稱佛宮寺釋迦塔，位于山西省應縣縣城西北角，始建于遼清寧二年（公元1056年），是中國現存的最高、最古老、保存最完整的木塔建筑。應縣木塔由上部木結構和臺基組成，如圖1（a）所示。上部木結構平面為八角形，包括5個明層和4個平坐層，共9層，高度達65.878 m。臺基是中國古建筑的獨特組成部分，是指高出地面以承托建筑物的底座，其不僅使建筑物防潮、防腐，還具有重要的文化和建筑價值[1]。應縣木塔的臺基包括八角形臺基以及四方形的外圍臺基（圖1（b））。

一些學者已經對應縣木塔的力學性能進行了研究。陳志勇等[2-3]通過低周往復水平加載和豎向單調加載試驗，研究了其柱腳節點、梁柱節點和斗栱節點的力學性能，并采用有限元軟件建立了精細化有限元模型，得到了其模態參數，同時也對木塔的抗震性能進行了研究。李鐵英等[4]制作了應縣木塔 1∶10結構模型并進行了擬動力試驗研究，計算了木塔在不同地震烈度下的彈塑性地震響應，得到了木塔實體結構非線性恢復力曲線，并提出了木結構古建筑的雙參數地震損壞準則。

相關研究表明，在對有臺基的木結構古建筑進行地震響應分析時，不能忽略臺基對上部木結構的影響[5-6]。目前，對考慮臺基影響的應縣木塔地震響應的研究還較少[7-8]。為了準確評估臺基對應縣木塔地震響應的影響，采用ANSYS有限元軟件建立了臺基、上部木結構及考慮臺基影響的整體結構有限元模型[9]，得到了各模型的自振頻率及周期。通過對上部木結構和考慮臺基影響的整體結構輸入地震波激勵，對比兩模型中參考點的加速度、位移時程曲線和結構動力系數，研究臺基對應縣木塔地震響應的影響。

1 模型建立

1.1 材料參數

對應縣木塔替換梁柱構件的樹種材料進行材性試驗[10]可知，其梁、柱、闌額及支撐等主要受力構件的木材皆為華北落葉松。華北落葉松的密度為548 kg/m3，泊松比為0.34，其他材性參數見表1。木材順紋受壓和順紋受拉的彈性模量取相等值[11]。

八角形臺基是應縣木塔的直接持力層。依據應縣木塔塔區工程地質條件勘察結果[10]可知，八角形臺基為夯實填土，厚度約5 m。夯土土質為密實粉砂質粉土，其材性參數見表2。為避免夯土材料受剪時土顆粒產生體積膨脹現象，膨脹角取為0°。八角形臺基下部的四方形臺基夯實程度較低，主要作用是保護八角形臺基，增強其穩定性。

1.2 上部木結構及構件尺寸

整體木結構建筑空間尺寸、木構件尺寸及細節構造參考陳明達[12]的測繪草圖。木結構主要的結構構件包括木柱、闌額、乳栿和支撐等。雖然構件數量較多，種類豐富，但主要的截面形式相對統一，主要的構件截面尺寸見表3。

1.3 單元及基本參數

采用BEAM188單元模擬木柱、闌額、乳栿和支撐的力學行為;采用COMBIN39單元模擬榫卯節點的半剛性以及斗栱在水平力作用下的非線性特征;采用COMBIN14單元模擬斗栱在豎向荷載作用下的力學行為。八角形臺基夯土部分采用Drucker-Prager模型，單元則采用20節點的SOLID186實體單元。

斗栱節點的水平彈簧和豎向彈簧的實常數采用袁建力等[13]進行的應縣木塔柱頭、補間及轉角鋪作1∶3比例模型的水平擬靜力和豎向靜力試驗中所得數據。采用陳志勇等[2]完成的應縣木塔明層柱闌額1∶3.4縮尺模型的擬靜力試驗結果確定榫卯節點彈簧單元的實常數。

1.4 有限元模型的建立

應縣木塔一層明層外槽與內槽柱間修筑了高大的土墼墻，對木柱有較強的嵌固作用，并增加了木塔1層的抗側剛度。有限元建模時將土墼墻等效為柱間斜撐，并將木塔1層木構件的彈性模量乘以10倍的放大系數[3]，以增大模型1層的抗側剛度。為考慮外圍四方形臺基對八角形臺基的約束作用，對八角形臺基0～1.66 m高度區間內各結點施加UX、UY、UZ、ROTX、ROTY、ROTZ等6個方向的約束。由于應縣木塔1層內外槽木柱皆被土墼墻包裹，木柱柱底與礎石之間不易產生滑移，將模型柱腳約束條件設為鉸接。所建立的上部木結構及考慮臺基影響的整體結構有限元模型如圖2所示。

2 模態分析

采用Block Lanczos法計算得到各模型的振型、自振頻率及周期。圖3分別為八角形臺基1階（東西平動）和3階扭轉振型圖。表4為八角形臺基的前7階自振頻率及周期。由表4和圖3可知，八角形臺基具有較高的自振頻率且各階頻率比較接近，這一方面是因為其內部夯土較密實，另一方面外圍四邊形臺基也提供了較強的側向約束，增加了其抗側剛度。

表5為上部木結構和考慮臺基影響的整體結構前7階平動自振頻率和周期。由表5可知，上部木結構自振頻率變化幅度較大。在低階頻率范圍內，臺基對整體結構的自振頻率的影響可忽略不計，整體結構的自振頻率與上部木結構接近;但當上部木結構的高階頻率與臺基頻率接近時，整體結構的自振頻率開始接近臺基的自振頻率。這說明整體結構的低階頻率主要取決于上部木結構的頻率，而其高階頻率受臺基的影響較大。

上部木結構的一階扭轉頻率為0.852 Hz，略大于1階平動頻率，但小于2階平動頻率。說明木塔上部木結構的抗扭轉剛度較小，在地震作用下，上部木結構發生彎曲變形的同時也容易發生扭轉變形。

為了驗證模擬結果的正確性，將上部木結構平動振型的自振頻率模擬值與文獻[14]實測值進行對比，結果見表6。由表6可知，平動頻率的模擬計算值與實測值的相對誤差較小，說明模擬結果具有一定的準確性。

3 地震響應分析

為提升計算效率，首先對上部木結構輸入三向Taft波以分析其弱水平剛度方向，然后沿弱水平剛度方向和豎直方向分別輸入El Centro波、Taft波及Tianjin波，以研究木塔的地震響應規律。

3.1 上部木結構弱水平剛度方向分析

建立了1∶1的2層明層及3層平座層有限元模型（見圖4），并對其輸入X向（木塔東西方向）、Y向（木塔南北方向）和Z向（木塔豎直方向）Taft波地震激勵，主震方向分別取X向和Y向，加速度峰值取250 cm/s2。X向、Y向及Z向輸入激勵的峰值加速度比分別為1∶0.85∶0.65和0.85∶1∶0.65。

圖5為不同主激振方向參考點（2層明層外槽柱頂測點）的位移時程曲線。當X向為主激振方向時，參考點X向最大位移為207.12 mm，Y向最大位移為190.80 mm。

X向最大位移比Y向大8.60%。當Y向為主激振方向時，參考點X向最大位移為175.14 mm，Y向最大位移為224.40 mm，Y向最大位移較X向大28.13%。參考點在Y向為主激振方向時，Y向的最大位移較X向為主激振方向時X向最大位移大8.15%。說明上部木結構的Y向抗側剛度小于X向，上部木結構的Y向（東西方向）為弱水平剛度方向。

圖6為不同主激振方向參考點加速度時程曲線對比圖。可以看出，X向為主激振方向時，X向最大加速度為196 cm/s2;Y向為主激振方向時，Y向最大加速度為255 cm/s2。Y向最大加速度較X向大30.10%。說明木結構模型的Y向地震響應明顯大于X向，上部木結構的Y向為弱水平剛度方向，這與位移時程分析所得結論一致。

3.2 上部木結構和考慮臺基影響的整體結構地震響應分析

進行結構地震響應分析時，可選擇與木塔所在場地特征相吻合的實際地震記錄。木塔地基土偏軟，基本屬于中硬Ⅱ類場地土[15]，故選用能反映該場地條件且應用較為廣泛的El Centro波、Taft波。由于木塔整體剛度較小，為考察長周期地震激勵的影響，同時選用Tianjin波作為地震激勵。分別對上部木結構和考慮臺基影響的整體結構沿Y向3種地震波和Z向輸入峰值加速度50、100 cm/s2。Y向和Z向輸入峰值加速度比取1∶0.65。各工況均以木塔5層明層外槽木柱柱頂為參考點。

圖7為不同工況上部木結構和考慮臺基影響的整體結構的位移時程曲線。可以看出，在相同峰值加速度作用下，3條地震波所引起的模型結構位移響應絕對值大小排序為Tianjin波>El Centro波≈Taft波。當峰值加速度相同時，Tianjin波激勵下的上部木結構和考慮臺基影響的整體結構參考點的最大位移明顯大于其他地震波激勵時對應的最大位移。這是由于Tianjin波屬于長周期地震波，其主要頻率集中于0～3 Hz，與上部木結構和考慮臺基影響的整體結構的自振頻率分布范圍相重合。由此可見，具有長周期特征的地震波激勵會對應縣木塔產生嚴重影響，在木塔日常振動監測中應給予重點關注。

表7、表8分別為峰值加速度為50、100 cm/s2的不同地震波作用下上部木結構和考慮臺基影響的整體結構的地震響應。隨輸入地震波峰值加速的增大，上部木結構和考慮臺基影響的整體結構的位移和加速度響應均增大。圖8為在峰值加速度為50 cm/s2和100 cm/s2的不同地震波激勵下上部木結構和考慮臺基影響的整體結構的動力系數，各工況的動力系數基本小于1，木塔結構上測點的地震響應較輸入激勵均減小。這是由于上部木結構中半剛性榫卯節點以及斗栱鋪作層發揮了較明顯的減隔震作用。由圖8、表7及表8可知，除El Centro波激勵外，隨輸入激勵的峰值加速度增大，半剛性榫卯節點和各層鋪作層耗能減震作用不斷提升，上部木結構和考慮臺基影響的整體結構的動力系數均減小。總體來看，應縣木塔具有較良好的耗能減震能力，在較小地震波激勵下，其耗能能力隨著地震動強度的增大而有所增強。

圖9為在峰值加速度為50、100 cm/s2時El Centro波激勵下上部木結構和考慮臺基影響的整體結構的位移和加速度時程曲線。可以看出，在相同峰值加速度的同種地震波激勵下，考慮臺基影響的整體結構的加速度和位移響應均大于上部木結構。因此，對應縣木塔進行地震響應分析時，若不考慮八角形臺基的影響會導致計算結果偏不安全。

4 結論

建立了應縣木塔臺基、上部木結構及考慮臺基影響的整體結構有限元模型，分別對上部木結構和考慮臺基影響的整體結構模型進行了地震響應分析。得到以下主要結論：

1）八角形臺基具有較高的自振頻率且各階頻率比較接近，當上部木結構的高階頻率與臺基頻率接近時，上部木結構自振頻率受臺基影響較大。

2）上部木結構的Y向（木塔東西方向）抗側剛度小于X向（木塔南北方向），因此，Y向為弱水平剛度方向，可在木塔的東西方向布設交叉支撐，提高結構在弱剛度方向的抗側能力。

3）當峰值加速度相同時，不同地震波所引起的地震響應大小排序為：Tianjin波>El Centro波≈Taft波。具有長周期特征的地震波激勵對木塔影響較大，應重點加固破壞嚴重的榫卯節點，適當減小結構自振周期，以減小木塔的位移響應。

4）隨著輸入地震波峰值加速度的增大，上部木結構和考慮臺基影響的整體結構的動力系數均呈減小趨勢，在相同峰值加速度的同種地震波激勵下，考慮臺基影響的整體結構的加速度和位移響應均大于僅考慮上部木結構的情況。

參考文獻：

[1]宋曉濱， 吳亞杰， 羅烈， 等. 傳統樓閣式木結構塔振動臺試驗研究[J]. 建筑結構學報， 2017， 38（2）： 10-19.

SONG X B， WU Y J， LUO L， et al. Shaking table test on a traditional pavilion-style wooden pagoda [J]. Journal of Building Structures， 2017， 38（2）： 10-19. （in Chinese）

[2]陳志勇. 應縣木塔典型節點及結構受力性能研究[D]. 哈爾濱： 哈爾濱工業大學， 2011.

CHEN Z Y. Behaviour of typical joints and the structure of Yingxian wood pagoda [D]. Harbin： Harbin Institute of Technology， 2011. （in Chinese）

[3]陳志勇， 祝恩淳， 潘景龍. 應縣木塔精細化結構建模及水平受力性能分析[J]. 建筑結構學報， 2013， 34（9）： 150-158.

CHEN Z Y， ZHU E C， PAN J L. Lateral structural performance of Yingxian Wood Pagoda based on refined FE models [J]. Journal of Building Structures， 2013， 34（9）： 150-158. （in Chinese）

[4]李鐵英， 魏劍偉， 張善元， 等. 木結構雙參數地震損壞準則及應縣木塔地震反應評價[J]. 建筑結構學報， 2004， 25（2）： 91-98.

LI T Y， WEI J W， ZHANG S Y， et al. Double parameter seismic damage criterion on wooden structure and seismic response appraisement on Yingxian wooden tower [J]. Journal of Building Structures， 2004， 25（2）： 91-98. （in Chinese）

[5]趙鴻鐵， 馬輝， 薛建陽， 等. 高臺基古建筑木結構動力特性及地震反應分析[J]. 地震工程與工程振動， 2011， 31（3）： 115-121.

ZHAO H T， MA H， XUE J Y， et al. Analysis of dynamic characteristics and seismic responses of the ancient timber building on high-platform base [J]. Journal of Earthquake Engineering and Engineering Vibration， 2011， 31（3）： 115-121. （in Chinese）

[6]WU C W， XUE J Y， ZHOU S Q， et al. Seismic performance evaluation for a traditional Chinese timber-frame structure [J/OL]. International Journal of Architectural Heritage. https： //doi.org/10.1080/15583058.2021.1894503.

[7]XIE Q F， WANG L， ZHANG L P， et al. Seismic behaviour of a traditional timber structure： Shaking table tests， energy dissipation mechanism and damage assessment model [J]. Bulletin of Earthquake Engineering， 2019， 17（3）： 1689-1714.

[8]SAKAI H， IZUNO K， SAKO K， et al. Seismic assessment of Japanese traditional wooden structure by dynamic interaction numerical analysis of surrounding ground [J]. Journal of Natural Disaster Science， 2014， 35（1）： 1-20.

[9]WANG J， HE J X， YANG N， et al. Study on aseismic characteristics of Tibetan ancient timber structure [J]. Advances in Materials Science and Engineering， 2017， 2017： 1-15.

[10]侯衛東， 王林安， 永昕群. 應縣木塔保護研究[M]. 北京： 文物出版社， 2016.

HOU W D， WANG L A， YONG X Q. Research and protection of Yingxian Wooden Pagoda [M]. Beijing： Culture Relics Publishing House， 2016. （in Chinese）

[11]王瑞鵬. 考慮臺基影響的應縣木塔地震模擬分析及加固方法的探討[D]. 西安： 西安建筑科技大學， 2018.

WANG R P. Seismic analysis of Yingxian Wooden Tower considered the effect of tower base and the discussion for the new set of enhancement programme [D]. Xian： Xian University of Architecture and Technology， 2018. （in Chinese）

[12]陳明達. 應縣木塔[M]. 北京： 文物出版社， 1966.

CHEN M D. Yingxian Wooden Pagoda [M]. Beijing：Culture Relics Publishing House， 2016. （in Chinese）

[13]袁建力， 陳韋， 王玨， 等. 應縣木塔斗栱模型試驗研究[J]. 建筑結構學報， 2011， 32（7）： 66-72.

YUAN J L， CHEN W， WANG J， et al. Experimental research on bracket set models of Yingxian Timber Pagoda [J]. Journal of Building Structures， 2011， 32（7）： 66-72. （in Chinese）

[14]浩飛虎. 應縣木塔動力特性試驗研究及變形分析[D]. 西安：西安建筑科技大學， 2019.

HAO F H. Study on dynamic characteristics test and deformation analysis of Ying County Wooden Pagoda [D]. Xian： Xian University of Architecture & Technology， 2019. （in Chinese）

[15]魏劍偉， 李鐵英， 張善元， 等. 應縣木塔地基工程地質勘測與分析[J]. 工程地質學報， 2003， 11（1）： 70-78.

WEI J W， LI T Y， ZHANG S Y， et al. Engineering-geological surveying and analysis of Yingxian Wooden Tower foundation [J]. Journal of Engineering Geology， 2003， 11（1）： 70-78. （in Chinese）

（編輯 胡玲）

3743500338281