斗栱梁?jiǎn)卧Ｐ蛿?shù)值分析
——以山西應(yīng)縣木塔為例

施穎，張晨君，陳喆

（江蘇省水利勘測(cè)設(shè)計(jì)研究院有限公司，江蘇 揚(yáng)州225009）

斗栱梁?jiǎn)卧Ｐ蛿?shù)值分析
——以山西應(yīng)縣木塔為例

施穎，張晨君，陳喆

（江蘇省水利勘測(cè)設(shè)計(jì)研究院有限公司，江蘇 揚(yáng)州225009）

將現(xiàn)代數(shù)值模擬技術(shù)運(yùn)用于古建筑研究，構(gòu)建合理的有限元模型并對(duì)斗栱的力學(xué)性能進(jìn)行分析研究，是古建筑保護(hù)的一項(xiàng)基本工作。文章以山西應(yīng)縣木塔斗栱構(gòu)件為實(shí)例，運(yùn)用ANSYS有限元軟件構(gòu)建了能夠真實(shí)反映斗栱構(gòu)造特征及結(jié)構(gòu)性能的梁?jiǎn)卧邢拊Ｐ?，?duì)梁?jiǎn)卧Ｐ驮谪Q向荷載及水平低周反復(fù)荷載作用下的實(shí)驗(yàn)過(guò)程進(jìn)行模擬，對(duì)比分析了其力學(xué)性能、耗能及模擬分析效率。結(jié)果表明:梁?jiǎn)卧Ｐ涂捎行У乇磉_(dá)斗栱的力學(xué)性能，模擬豎向力加載時(shí)相應(yīng)的壓縮變形量為11.60 mm，模擬水平低周反復(fù)荷載加載時(shí)，隨著側(cè)向位移的增加模型由彈性狀態(tài)轉(zhuǎn)變?yōu)閺椝苄誀顟B(tài)，隨著豎向荷載的增加需要施加的水平推力增大，與實(shí)驗(yàn)結(jié)果吻合；梁?jiǎn)卧Ｐ涂捎行У乇磉_(dá)斗栱的實(shí)際耗散特性，其耗能值與實(shí)驗(yàn)結(jié)果基本吻合；梁?jiǎn)卧Ｐ涂商岣叻治鲂?，?jì)算分析耗時(shí)由7 h減少為30 min。

斗栱梁?jiǎn)卧Ｐ停粯?gòu)造特征；力學(xué)性能；耗散特性；數(shù)值分析

0 引言

隨著現(xiàn)代科學(xué)技術(shù)的發(fā)展，計(jì)算機(jī)數(shù)值模擬與分析技術(shù)已經(jīng)成為古建筑保護(hù)研究中的重要工具［1］，而斗栱是中國(guó)傳統(tǒng)木結(jié)構(gòu)建筑最有特色的連接構(gòu)件，被廣泛應(yīng)用于宮殿、寶塔和寺廟等重要的古建筑之中［2］。因此，將現(xiàn)代數(shù)值模擬技術(shù)運(yùn)用于古建筑研究，構(gòu)建合理的有限元模型并對(duì)斗栱的力學(xué)性能進(jìn)行分析研究，是古建筑保護(hù)的一項(xiàng)基本研究工作［3］。目前，運(yùn)用計(jì)算機(jī)模擬分析進(jìn)行古建筑修復(fù)研究已經(jīng)取得了一批較豐碩的成果。袁建力等基于摩擦—剪切耗能特性研制了斗栱實(shí)體單元有限元模型［4］；王玨等根據(jù)現(xiàn)場(chǎng)查勘數(shù)據(jù)，運(yùn)用有限元軟件構(gòu)建了體現(xiàn)木塔“殿堂結(jié)構(gòu)”的“斗栱—柱圈”層疊式有限元模型［5］；徐燊等運(yùn)用應(yīng)縣木塔現(xiàn)有測(cè)繪數(shù)據(jù)，以層為單位，構(gòu)建了整塔復(fù)合現(xiàn)狀的殘損有限元模型［6］；杜雷鳴等在應(yīng)縣木塔抗震性能研究中，將構(gòu)件的節(jié)點(diǎn)域分別簡(jiǎn)化為鉸接和剛接兩種類(lèi)型，在此基礎(chǔ)上建立了整塔有限元模型［7］。

但大部分斗栱有限元研究都是建立在實(shí)體單元基礎(chǔ)上進(jìn)行的，與一般梁?jiǎn)卧窎砟Ｐ拖啾?，?shí)體單元模型總單元數(shù)量相對(duì)較大，計(jì)算分析也更耗時(shí)。對(duì)于古建筑木結(jié)構(gòu)的力學(xué)分析而言，其承重構(gòu)架主要由柱、闌額、普拍枋等具有梁?jiǎn)卧阅艿臉?gòu)件組成，因此文章針對(duì)斗栱三維實(shí)體有限元模型［8］進(jìn)一步簡(jiǎn)化，以山西應(yīng)縣木塔斗栱構(gòu)件為實(shí)例，構(gòu)建能夠真實(shí)反映斗栱構(gòu)造特征及結(jié)構(gòu)性能的梁?jiǎn)卧邢拊Ｐ停M其在豎向荷載及水平低周反復(fù)荷載作用下的實(shí)驗(yàn)過(guò)程，對(duì)比分析梁?jiǎn)卧Ｐ团c試驗(yàn)?zāi)Ｐ蛯?duì)應(yīng)的力學(xué)性能、耗能及梁?jiǎn)卧Ｐ团c實(shí)體模型的模擬分析效率，驗(yàn)證模型的有效性和可行性。

1 斗栱梁?jiǎn)卧Ｐ偷慕?/h2>

1.1 斗栱構(gòu)件的構(gòu)造特征

山西應(yīng)縣木塔中采用的宋式斗栱的基本構(gòu)造［9］，如圖1所示。采用梁?jiǎn)卧r(shí)，栱的各部分可直接簡(jiǎn)化為同尺寸梁?jiǎn)卧?，?shí)際構(gòu)造中縱橫向栱之間采用槽口卡實(shí)的搭接方式，受力時(shí)基本無(wú)相對(duì)滑動(dòng)，故該部位簡(jiǎn)化模型中做固接處理；斗與栱銷(xiāo)釘連接部位也采用固接處理。為反映真實(shí)受力情況，以更接近實(shí)際構(gòu)造情況，應(yīng)考慮將栱轉(zhuǎn)角部分做倒角處理，防止直角引起應(yīng)力集中。最后，在為梁?jiǎn)卧峙鋯卧獙傩缘倪^(guò)程中，對(duì)部分截面進(jìn)行了軸偏心處理，以使建成的梁?jiǎn)卧Ｐ妥畲蟪潭鹊嘏c實(shí)體模型相似。

圖1 宋式斗栱基本構(gòu)造圖

斗栱的斗部分構(gòu)造較為復(fù)雜，在方形的斗面四個(gè)角點(diǎn)設(shè)置四斗栮，實(shí)際受力時(shí)斗栮主要是起約束過(guò)大位移作用的構(gòu)件［10］，建模時(shí)必須加以考慮。這樣的構(gòu)造特征在實(shí)體模型建模時(shí)，通過(guò)布爾操作挖空塊體底部，很容易實(shí)現(xiàn)斗栮的構(gòu)造，但在梁?jiǎn)卧Ｐ椭袃H用一根桿件很難做到，需要考慮通過(guò)多根桿件組合實(shí)現(xiàn)。在簡(jiǎn)化過(guò)程中將方形的斗面簡(jiǎn)化成由12根梁固接成的“田”字形梁?jiǎn)卧M合，并在“田”字四個(gè)角點(diǎn)分別固接與實(shí)際斗栮尺寸相同的梁?jiǎn)卧讨鶃?lái)體現(xiàn)這種四栮交互斗的構(gòu)造，如圖2所示。另外一種雙栮散斗則是利用相同的原理在兩根梁組成的斗身端部設(shè)置兩梁?jiǎn)卧讨鶚?gòu)成，短柱與上部簡(jiǎn)化斗體均采用剛性連接，以此確保斗栮正常工作。

圖2 斗的簡(jiǎn)化圖

1.2 斗栱構(gòu)件的力學(xué)性能

根據(jù)對(duì)斗栱模型進(jìn)行豎向荷載試驗(yàn)及水平低周反復(fù)荷載實(shí)驗(yàn)可知:（1）在豎向荷載作用下，壓力沿著豎向軸心傳遞到櫨斗上，彎矩和剪力由泥道栱傳遞至櫨斗，導(dǎo)致櫨斗在彎、剪、壓復(fù)雜受力狀態(tài)下劈裂；（2）在水平低周反復(fù)荷載作用下，泥道栱受到上部傳來(lái)的壓力、斗與栱之間的摩擦力、櫨斗傳來(lái)的水平推力的共同作用，使其處于壓、彎、剪的復(fù)雜應(yīng)力狀態(tài)，導(dǎo)致薄弱面的水平劈裂破壞［11］。因此，斗栱梁?jiǎn)卧?jiǎn)化模型選用BEAM189單元，該單元既可考慮橫截面剪切變形亦可考慮截面的翹曲變形，支持彈性、粘彈、粘塑、塑性和蠕變模型，具備分析橫向、彎曲或扭向的穩(wěn)定問(wèn)題的能力［12］，可以反映斗栱受力分析中需要考慮的各種特性。

斗栱在水平荷載作用下構(gòu)件之間的摩阻力隨水平推力增加而不斷增大，體現(xiàn)出穩(wěn)定的摩擦剪切耗能特性。為了體現(xiàn)這一耗能特性，在斗與栱的接觸范圍定義摩擦接觸面［12］，由CONTA 175單元沿“斗”內(nèi)布置接觸單元，由TARGE 169單元在“栱”滑動(dòng)區(qū)域設(shè)置目標(biāo)單元［12］，CONTA 175接觸單元和TARGE 169目標(biāo)單元的設(shè)置如圖3所示。

1.3 斗栱梁?jiǎn)卧Ｐ偷牟牧咸匦?/p>

建立斗栱梁?jiǎn)卧治瞿Ｐ蜁r(shí)，以力學(xué)實(shí)驗(yàn)采用物理模型的各項(xiàng)參數(shù)為依據(jù)確定模型的材料特性。木材的密度為548 kg／m3，泊松比為0.17。木材為典型的各向異性材料，通過(guò)木材試件軸心受壓試驗(yàn)測(cè)得:木材橫紋的屈服應(yīng)力為1.67 MPa、彈性模量為216 MPa，硬化模量為3.8 MPa，橫紋應(yīng)力—應(yīng)變曲線采用隨動(dòng)雙線性硬化本構(gòu)模型（BISO）輸入；木材順紋的屈服應(yīng)力取為23.08 MPa，彈性模量為907 MPa，硬化模量為44.7 MPa。順紋應(yīng)力—應(yīng)變曲線采用隨動(dòng)多線性硬化本構(gòu)模型（MISO）輸入，輸入的應(yīng)力—應(yīng)變曲線如圖4所示。

圖3 CONTA 175接觸單元和TARGE 169目標(biāo)單元設(shè)置圖

圖4 木材受壓應(yīng)力—應(yīng)變曲線圖

木材的摩擦系數(shù)是基于摩擦耗能有限元模型的關(guān)鍵指標(biāo)?？紤]到斗栱各構(gòu)件在水平荷載作用下之間相對(duì)位移較小，可采用靜摩擦系數(shù)以簡(jiǎn)便結(jié)構(gòu)分析。根據(jù)文獻(xiàn)［4］中木材靜摩擦系數(shù)測(cè)試結(jié)果，確定所用木材橫紋靜摩擦系數(shù)為0.519、順紋靜摩擦系數(shù)為0.455。

1.4 斗栱梁?jiǎn)卧Ｐ偷慕?/p>

根據(jù)上述分析，運(yùn)用ANSYS有限元軟件構(gòu)建完成斗栱柱頭、補(bǔ)間及轉(zhuǎn)角鋪?zhàn)魅N特征性斗栱梁?jiǎn)卧邢拊Ｐ头謩e如圖5所示。

圖5 梁?jiǎn)卧邢拊Ｐ突緲?gòu)造圖

2 斗栱梁?jiǎn)卧Ｐ驮诓煌奢d作用下的模擬結(jié)果

2.1 豎向荷載作用

唐代詩(shī)人宋之問(wèn)，身材高昂、儀表堂堂，進(jìn)士及第，擅寫(xiě)五言詩(shī)，名氣不小。但他的為人卻多為義士譏諷，傳說(shuō)其外甥劉希夷寫(xiě)出新詩(shī)名句“年年歲歲花相似，歲歲年年人不同”。宋之問(wèn)求外甥將這兩句詩(shī)的版權(quán)送給他。劉希夷不肯，再求，還是不肯，這缺德的舅舅一怒之下，竟叫仆人用土沙袋壓死了外甥。此外，宋又是機(jī)會(huì)主義者，為了仕途，攀附武后幸臣張易之、張宗昌，二張被誅，宋之問(wèn)被貶瀧州（今廣東羅定縣）參軍，諸事艱難。慕念昔榮，宋之問(wèn)在次年春秘密逃還洛陽(yáng)，途經(jīng)漢江（襄陽(yáng)附近的漢水），寫(xiě)了一首五言絕句《渡漢江》：“嶺外音書(shū)斷，經(jīng)冬復(fù)歷春。近鄉(xiāng)情更怯，不敢問(wèn)來(lái)人?！?/p>

根據(jù)斗栱實(shí)驗(yàn)過(guò)程分析，去掉了模型試驗(yàn)中試件不斷壓實(shí)時(shí)豎向力N為0～10 kN的過(guò)程，施加豎向力10～40 kN，分15個(gè)荷載步進(jìn)行模擬，加載時(shí)約束模型底部位移，梁?jiǎn)卧Ｐ椭袡径窐?gòu)造類(lèi)似一底部固定，上部自由的短柱，直接施加集中荷載會(huì)引起應(yīng)力集中甚至是失穩(wěn)現(xiàn)象，這里考慮將荷載等分為五份分別施加于櫨斗斗面各交點(diǎn)。以柱頭鋪?zhàn)髁簡(jiǎn)卧Ｐ蜑槔谪Q向荷載作用下的過(guò)程及模擬結(jié)果如圖6所示。

圖6 豎向荷載作用下柱頭鋪?zhàn)髁簡(jiǎn)卧M過(guò)程及結(jié)果圖

2.2 水平低周反復(fù)荷載作用

模擬水平低周反復(fù)荷載作用下的受力過(guò)程時(shí)，加載方式采用水平力單調(diào)加載。分別對(duì)柱頭、補(bǔ)間和轉(zhuǎn)角鋪?zhàn)鞫窎砟Ｐ褪┘迂Q向荷載N為15、22.5 和30 kN，并逐步施加水平荷載，考察當(dāng)水平位移達(dá)到各級(jí)控制位移對(duì)應(yīng)的水平荷載和相應(yīng)的耗能面積。各模型三種豎向荷載（工況）下的模擬結(jié)果分別如圖7～9所示。

3 斗栱梁?jiǎn)卧Ｐ湍M結(jié)果對(duì)比分析

袁建力等以應(yīng)縣木塔中的代表性斗栱為對(duì)象，制作了應(yīng)縣木塔中的斗栱實(shí)驗(yàn)?zāi)Ｐ?，采用擬靜力研究的方法，分別進(jìn)行了模型豎向荷載作用和水平低周反復(fù)荷載作用下的極限承載力試驗(yàn)［2］；袁建力等在此基礎(chǔ)上研制了反應(yīng)斗栱摩擦剪切耗能的實(shí)體單元有限元分析模型，從數(shù)值計(jì)算的角度驗(yàn)證了實(shí)體單元模型可有效代替斗栱實(shí)驗(yàn)?zāi)Ｐ瓦M(jìn)行模擬分析［4］，如圖10所示。

以柱頭鋪?zhàn)鳛槔瑢?duì)比分析實(shí)驗(yàn)?zāi)Ｐ团c梁?jiǎn)卧Ｐ土W(xué)性能、耗散性能及實(shí)體單元模型和梁?jiǎn)卧Ｐ偷姆治鲂?，?yàn)證梁?jiǎn)卧Ｐ痛鎸?shí)體模型進(jìn)行簡(jiǎn)化模擬分析的可行性。

圖7 柱頭鋪?zhàn)髁簡(jiǎn)卧Ｐ湍M結(jié)果圖

圖8 補(bǔ)間鋪?zhàn)髁簡(jiǎn)卧Ｐ湍M結(jié)果圖

圖9 轉(zhuǎn)角鋪?zhàn)髁簡(jiǎn)卧Ｐ湍M結(jié)果圖

圖10 實(shí)驗(yàn)?zāi)Ｐ秃蛯?shí)體單元模型圖

3.1 力學(xué)性能分析

3.1.1 豎向荷載作用

梁?jiǎn)卧Ｐ驮谪Q向力加載的過(guò)程中完全處于線彈性狀態(tài)，豎向荷載沿軸心線上的櫨斗、齊心斗向下傳遞，根據(jù)圖6模型相應(yīng)的壓縮變形量為11.60 mm與實(shí)驗(yàn)得到的壓縮變形量11.20 mm［2］基本一致。

3.1.2 水平低周反復(fù)荷載作用

根據(jù)圖7～9，得出3種模型在N為15、22.5、30 kN等3種工況下的水平位移達(dá)到各級(jí)控制位移時(shí)所需的作用力見(jiàn)表1。

表1 三種柱頭鋪?zhàn)髂Ｐ驮诟骺刂莆灰茣r(shí)對(duì)應(yīng)的水平力／kN

從表1中數(shù)據(jù)可以看出:在同一工況下，側(cè)向控制位移較小時(shí)（即3mm增加至12mm時(shí)）梁?jiǎn)卧Ｐ突咎幱趶椥誀顟B(tài)，承載力與側(cè)向位移呈線性比例關(guān)系；隨著側(cè)向控制位移d的增加（由12 mm增加至18 mm），模型呈彈塑性狀態(tài)，其承載力相應(yīng)下降，表明構(gòu)件出現(xiàn)損傷［14］；在同一級(jí)控制位移下，以柱頭鋪?zhàn)髟?5 mm控制位移時(shí)的工況為例，隨著豎向荷載N由15 kN增加至30 kN，每一加載部需要實(shí)加的水平推力由6.79 kN增加至9.01 kN與實(shí)驗(yàn)?zāi)Ｐ?.97 kN增加至8.79 kN的實(shí)驗(yàn)結(jié)果大致相同。數(shù)值對(duì)比結(jié)果表明，梁?jiǎn)卧Ｐ涂捎行У乇磉_(dá)斗栱的力學(xué)性能。

3.2 耗散性能分析

根據(jù)圖7～9，統(tǒng)計(jì)出3種模型在N為15、22.5、30 kN等3種工況下根據(jù)滯回曲線中滯回環(huán)面積計(jì)算的能量耗散值，見(jiàn)表2。

表2 三種柱頭鋪?zhàn)髂Ｐ驮诟骺刂莆灰茣r(shí)耗能值／J

滯回曲線中滯回環(huán)的面積可以反映能量吸收和耗散的能力，這些能量是通過(guò)材料的內(nèi)摩阻或局部損傷而散失的，因此，該指標(biāo)是評(píng)定結(jié)構(gòu)耗能的重要指標(biāo)［15］，從表2可以看出，以柱頭鋪?zhàn)髟谪Q向荷載15 kN、側(cè)向控制位移15 mm的工況為例，梁?jiǎn)卧Ｐ秃哪苤禐?125.67 J與實(shí)驗(yàn)?zāi)Ｐ秃哪苤禐?144.03 J基本吻合，說(shuō)明梁?jiǎn)卧Ｐ涂捎行У乇磉_(dá)斗栱的實(shí)際耗散特性

3.3 模型效率分析

對(duì)比梁?jiǎn)卧Ｐ图皩?shí)體單元模型在同一網(wǎng)格精度下各類(lèi)單元數(shù)量，見(jiàn)表3。

表3 梁?jiǎn)卧蛯?shí)體模型構(gòu)件組成 ／個(gè)

從表3可以看出，在同一精度下，以柱頭鋪?zhàn)鳛槔簡(jiǎn)卧Ｐ蛦卧獢?shù)量768個(gè)，遠(yuǎn)小于鋪?zhàn)鲗?shí)體模型142693個(gè)；同一鋪?zhàn)髂Ｐ驮谕缓奢d工況作用下梁?jiǎn)卧Ｐ头治龊臅r(shí)約30 min，而同一實(shí)體模型分析耗時(shí)約需要7 h，模擬速度得到有效的提高，分析結(jié)果占用空間也得以有效減少。這使得在此基礎(chǔ)上建立應(yīng)縣木整塔梁?jiǎn)卧Ｐ?，進(jìn)行動(dòng)力特性分析成為可能。

4 結(jié)論

通過(guò)上述研究表明:

（1）梁?jiǎn)卧Ｐ涂捎行У乇磉_(dá)斗栱的力學(xué)性能，模擬豎向力加載時(shí)相應(yīng)的柱頭鋪?zhàn)髁簡(jiǎn)卧Ｐ偷膲嚎s變形量為11.60 mm；模擬水平低周反復(fù)荷載加載時(shí)，隨著側(cè)向位移由3 mm增加至18 mm，模型由彈性狀態(tài)逐漸轉(zhuǎn)變?yōu)閺椝苄誀顟B(tài)，隨著豎向荷載由15 kN加載至30 kN柱頭鋪?zhàn)髁簡(jiǎn)卧Ｐ托枰┘拥乃酵屏τ?.79 kN增加至9.01 kN，與實(shí)驗(yàn)結(jié)果基本吻合。

（2）梁?jiǎn)卧Ｐ涂捎行У乇磉_(dá)斗栱的實(shí)際耗散特性，經(jīng)對(duì)比分析梁?jiǎn)卧Ｐ偷母骺刂莆灰频暮哪苤蹬c實(shí)驗(yàn)結(jié)果基本吻合。

（3）梁?jiǎn)卧Ｐ涂商岣叻治鲂?，與一般實(shí)體單元斗栱模型相比，柱頭鋪?zhàn)鞯牧簡(jiǎn)卧Ｐ偷目倖卧獢?shù)由142693個(gè)減少為768個(gè)，計(jì)算分析耗時(shí)由7 h減少為30 min，模擬速度得到有效的提高。

［1］ 袁建力，基于現(xiàn)代模擬分析技術(shù)的傳統(tǒng)木結(jié)構(gòu)研究［J］.土木建筑工程信息基礎(chǔ)，2013，5（2）:20-30.

［2］ 袁建力，陳韋，王玨，等.應(yīng)縣木塔斗栱模型試驗(yàn)研究［J］.建筑結(jié)構(gòu)學(xué)報(bào)，2011，32（7）:66-72.

［3］ 李飛.清式大木作古建筑斗栱抗震性能研究［D］.西安:西安建筑科技大學(xué)，2014.

［4］ 袁建力，施穎，陳韋，等.基于摩擦—剪切耗能的斗栱有限元模型研究［J］.建筑結(jié)構(gòu)學(xué)報(bào)，2013，33（6）:151-157.

［5］ 王玨.應(yīng)縣木塔扭、傾變形張拉復(fù)位的數(shù)字化模擬和安全性評(píng)價(jià)［D］.揚(yáng)州:揚(yáng)州大學(xué)，2008.

［6］ 徐燊，李子劍，劉魯.高效的建筑動(dòng)態(tài)建模方法研究——在應(yīng)縣木塔三維復(fù)原中的應(yīng)用與分析［J］.華中建筑，2003，21（1）: 54-57.

［7］ 杜雷鳴，李海旺，薛飛，等.應(yīng)縣木塔抗震性能研究［J］.土木工程學(xué)報(bào)，2010（S1）:363-370.

［8］ 石廣斌.實(shí)體單元與梁?jiǎn)卧獢?shù)值分析的差異性剖析及應(yīng)用［J］.西北水電，2012（1）:61-65.

［9］ 包軼楠.古建木結(jié)構(gòu)榫卯節(jié)點(diǎn)的力學(xué)性能研究［D］.南京:東南大學(xué)，2014.

［10］陳明達(dá).應(yīng)縣木塔［M］.北京:文物出版社，2001.

［11］王智華.應(yīng)縣木塔斗栱調(diào)查與力學(xué)性能分析［D］.西安:西安建筑科技大學(xué)，2010.

［12］侯靜，陳金永，萬(wàn)佳，等.基于有限元技術(shù)的傳統(tǒng)木結(jié)構(gòu)關(guān)鍵構(gòu)造節(jié)點(diǎn)的性能研究與分析［J］.科學(xué)技術(shù)與工程，2015，15 （28）:60-66.

［13］宋曉勝.中國(guó)木結(jié)構(gòu)古建筑結(jié)構(gòu)性能與加固保護(hù)研究［D］.北京:北京工業(yè)大學(xué)，2014.

［14］謝啟芳，鄭培君，向偉，等.殘損古建筑木結(jié)構(gòu)單向直榫榫卯節(jié)點(diǎn)抗震性能試驗(yàn)研究［J］.建筑結(jié)構(gòu)學(xué)報(bào)，2014，35（11）:143 -150.

［15］夏海倫.不同松動(dòng)程度下古建筑透榫節(jié)點(diǎn)抗震性能試驗(yàn)研究［D］.西安:西安建筑科技大學(xué)，2015.

（學(xué)科責(zé)編：趙成龍）

Simulation analysis based on Dougong beam elementmodel: A case analysis of Yingxian'sW ood Tower in Shanxi

Shi Ying，Zhang Chenjun，Chen Zhe
（Jiangsu Surveying and Design Institute ofWater Resources Co.，Ltd.，Yangzhou 225009，China）

Applying themodern numerical simulation technology into the study of ancient buildings，and constructing the reasonable FEmodel to analyzemechanical property of Dou-Gong has become the basic work for the protection of ancient buildings.This research proposed the FE model base beam element to truthfully reflect structural characteristics and properties of Dougong in Yingxian's Wood Tower in Shanxi，simulates the experimentation process of the beam elementmodel under the action of vertical loading and horizontal low-cyclic reversed loading，and comparatively analyzes themechanical properties，energy consumption and simulation analysis efficiency.The models have the following features:The beam elementmodel can express the mechanical property of Dougong effectively，and the amount of compression deformation corresponding to simulating vertical loading is11.6mm.When simulating cyclic reverse loading，the model translates elastic state into elastic-plastic state with the increase of the lateral displacement，and the horizontal thrust increases with the increase of vertical load，which is in agreement with the result of experiment.The beam elementmodel can express the actual dispersion features of Dougong effectively，whose energy dissipation is basically identical with the experimental result and it can also improve analysis efficiency，and the time of calculation and analysis is reduced from 7 hours to 30 minutes.

Dougong beam elementmodel；configuration character；mechanical property；dissipativeproperty；numerical analysis

TU317.1

A

1673-7644（2017）03-0295-07

2017-04-25

施穎（1987-），女，工程師，碩士，主要從事建筑結(jié)構(gòu)等方面的研究.E-mail:340272859＠qq.com