圓形土樓夯土結(jié)構(gòu)性能的數(shù)值模擬

葉俊捷， 彭興黔， 施維娟

(華僑大學(xué) 土木工程學(xué)院， 福建 廈門 361021)

?

圓形土樓夯土結(jié)構(gòu)性能的數(shù)值模擬

葉俊捷， 彭興黔， 施維娟

(華僑大學(xué) 土木工程學(xué)院， 福建 廈門 361021)

采用數(shù)值模擬方法，以福建龍巖的振成樓、承啟樓等經(jīng)典圓形土樓為原型，建立圓形土樓夯土結(jié)構(gòu)模型.對(duì)土樓夯土墻體在靜力荷載作用下的結(jié)構(gòu)承載性能進(jìn)行分析，通過(guò)改變土樓各尺寸參數(shù)，分析不同土樓夯土結(jié)構(gòu)模型的受力及變形特征，得出影響受力變形規(guī)律特征的主要因素.結(jié)果表明：不論是在墻厚、直徑還是墻高模擬中，應(yīng)力、位移都不是單調(diào)變化的，有合理的取值區(qū)間.

夯土結(jié)構(gòu)； 圓形土樓； 承載性能； 變形特征； 結(jié)構(gòu)參數(shù)； 數(shù)值模擬

福建土樓作為我國(guó)歷史文化的一張名片，具有極高的歷史、政治、藝術(shù)、旅游、科學(xué)和文化價(jià)值[1].在眾多土樓群中，各個(gè)土樓的尺寸參數(shù)都各不相同，導(dǎo)致結(jié)構(gòu)性能也不同，這也是導(dǎo)致其在長(zhǎng)期相同的環(huán)境作用下所受到的損壞不相同的重要原因之一.合理的尺寸能夠提高土樓結(jié)構(gòu)性能，也使得有些土樓能夠遺留百年，而有些土樓只能拆除或翻建.為了探究影響土樓結(jié)構(gòu)性能的主要尺寸參數(shù)，得出圓形土樓最合理的尺寸，采用有限元軟件進(jìn)行數(shù)值分析[2].采集振成樓、承啟樓等具有代表性的圓形土樓的相關(guān)尺寸參數(shù)，對(duì)這些數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，采用數(shù)理統(tǒng)計(jì)的方法，計(jì)算它們的均值，建立模型.所建立的模型更具有普遍性，更能客觀地反映各項(xiàng)參數(shù)的變化對(duì)土樓帶來(lái)的影響，更能探究土樓的結(jié)構(gòu)性隨著各參數(shù)改變的變化，利于分析土樓的結(jié)構(gòu)可靠性[3-4].本文對(duì)圓形土樓的結(jié)構(gòu)性能進(jìn)行數(shù)值模擬.

1 有限元模型的建立

1.1 幾何模型的描述

為了能更客觀、規(guī)律地將模擬結(jié)果反映到土樓上，模擬尺寸為振成樓、如升樓、承啟樓、僑福樓等圓形土樓的實(shí)際尺寸，包括外部直徑D、墻高h(yuǎn)、每層墻厚d.使用SPSS軟件進(jìn)行數(shù)據(jù)統(tǒng)計(jì)，以偏度-峰度檢驗(yàn)法及W檢驗(yàn)法[5]作為標(biāo)準(zhǔn)，選取平均值進(jìn)行數(shù)值模擬.模擬結(jié)果如下：直徑D為42 m；均勻墻厚d為1 m；層高h(yuǎn)為4.2 m；門的尺寸為1.5 m×2.0 m；窗戶的尺寸為 0.6 m×0.6 m.

1.2 本構(gòu)模型的建立

考慮基于ABAQUS提供的兩種劃分方式[6-7]的缺陷，根據(jù)實(shí)際比較，圓形土樓夯土結(jié)構(gòu)的幾何模型,如圖1所示.由于圓形土樓的特殊結(jié)構(gòu)，為使圓形土樓模型能獲得較高的精度，模型最終選取二次六面體單元，即采用結(jié)構(gòu)化或掃掠的網(wǎng)格劃分技術(shù)，得到的網(wǎng)格較為規(guī)則.網(wǎng)格劃分結(jié)果，如圖2所示.

圖1 圓形土樓夯土結(jié)構(gòu)的幾何模型 圖2 土樓夯土結(jié)構(gòu)模型的網(wǎng)格劃分Fig.1 Geometry model of circular Tulou Fig.2 Meshing of Tulou structure rammed earth structure rammed earth structure

通過(guò)對(duì)比巖土中各種彈塑性本構(gòu)模型[8-9]、夯土材料參數(shù)及實(shí)驗(yàn)數(shù)據(jù)，聯(lián)用ABAQUS中的彈性模型[10]與塑性模型中的Mohr-Coulomb模型，結(jié)合前期夯土材料試驗(yàn)數(shù)據(jù)，擬合成適用于土樓夯土材料的彈塑性本構(gòu)模型.在ABAQUS在分析計(jì)算時(shí)，采用連續(xù)光滑的橢圓函數(shù)代替?zhèn)鹘y(tǒng)的Mohr-Coulomb的尖角六邊形函數(shù)，即

(1)

式(1)中:Ψ為剪脹角；c0為初始粘聚力，即未產(chǎn)生塑性變形前的粘聚力;ε為控制子午面上的形狀及函數(shù)漸近線間相似度的偏心率.

假設(shè)Rmw(Θ,e,φ)為控制G在π平面上的形狀，即

(2)

式(2)中:e是π平面上的偏心率，即

(3)

通過(guò)計(jì)算，e使屈服面與塑性勢(shì)面在π平面上的受拉及受壓角點(diǎn)相切，e的大小根據(jù)所擬用的塑性勢(shì)面取值，但必須符合0.5≤e≤1.

2 試驗(yàn)方案選取

為了能夠全面了解土樓參數(shù)對(duì)土樓結(jié)構(gòu)性能的影響，選取如下3個(gè)對(duì)土樓結(jié)構(gòu)性能影響的參數(shù).

1) 變墻厚對(duì)比.土樓墻體作為主要的承重結(jié)構(gòu)，對(duì)土樓整體的結(jié)構(gòu)性能有著不可忽略的影響.此方案探究在不改變其余參數(shù)的情況下，改變墻體厚度，對(duì)土樓夯土墻的整體結(jié)構(gòu)受力變形性能的影響.

2) 變直徑對(duì)比.由于圓形土樓的直徑大小不一，通過(guò)改變土樓的直徑大小，探究土樓直徑的取值除了根據(jù)風(fēng)水及土樓樓主實(shí)際的用地需求外，還要根據(jù)土樓結(jié)構(gòu)的受力影響.

3) 變墻高對(duì)比.過(guò)高的墻體導(dǎo)致土樓在受力過(guò)程中產(chǎn)生失穩(wěn)破壞，并且土樓自身的墻體就有一定的傾斜程度，故不能忽略土樓墻體高度的影響.此方案研究改變土樓墻體高度的大小，對(duì)土樓夯土墻結(jié)構(gòu)受力變形性能的影響.

3 參數(shù)變化對(duì)結(jié)構(gòu)性能的影響

通過(guò)對(duì)土樓原狀土的夯土材料參數(shù)分析，根據(jù)土工試驗(yàn)，夯土墻密度取1 624kg·m-3；彈性模量為130.67MPa；泊松比為0.3；內(nèi)摩擦角為39°；粘聚力為410kPa.在施加荷載時(shí)，按照木梁與夯土墻體接觸面積加載，并針對(duì)不同模擬試驗(yàn)方案，分別施加土樓實(shí)際荷載或極限荷載.

3.1 變墻厚對(duì)比數(shù)值模擬結(jié)果

在保持土樓夯土墻體直徑(42m)、墻體高度(12.6m)為常數(shù)的情況下，變化墻體厚度.按照墻體厚度由小到大的模型為A～G(0.7～1.3m，梯度為0.1).分析不同墻體厚度下，土樓夯土結(jié)構(gòu)模型的極限狀態(tài).不同墻體厚度模型變化曲線，如圖3所示.圖3中：σ為應(yīng)力；s為位移.

由圖3可知：隨著墻體厚度的增大，土樓夯土結(jié)構(gòu)模型的承載能力呈現(xiàn)增大趨勢(shì)，并且不會(huì)引起較大的豎向位移值.但是較大的厚度會(huì)引起土樓建造成本的提高，在土樓實(shí)際的建造過(guò)程中，土樓樓主所選取的土樓墻體厚度也較為合理，不會(huì)單一地追求墻體厚度大小，而是從土樓的整體考慮，選取合適的土樓墻體厚度，這也是土樓能夠經(jīng)歷悠久歷史屹立不倒的原因之一.

(a) 最大應(yīng)力 (b) 最大豎向位移圖3 不同墻體厚度模型變化曲線Fig.3 Curve of models with the different wall thicknesses

3.2 變直徑對(duì)比數(shù)值模擬結(jié)果

在保持墻體厚度(1 m)、墻體高度(12.6 m)為常數(shù)的情況下，改變土樓夯土結(jié)構(gòu)模型的直徑，分析不同直徑下的極限狀態(tài)，分析結(jié)果如圖4所示.直徑分別由小至大模型編號(hào)為A～F(35.7～46.2 m，梯度為2.1).

(a) 最大應(yīng)力 (b) 最大豎向位移圖4 不同直徑模型變化曲線Fig.4 Curve of models with the different diameters

由圖4可知：直徑較大的F模型，承載力較大，且極限狀態(tài)的位移值反而減小，說(shuō)明增大直徑可以提高土樓的整體性，從而提高土樓的整體承載能力，且極限狀態(tài)下的豎向位移值不會(huì)有明顯的增大.在實(shí)際土樓群中，直徑更大的土樓結(jié)構(gòu)承受能力更大，也直觀地體現(xiàn)出歷史悠久的土樓一般都是直徑較大的原因.

3.3 變墻高數(shù)對(duì)比值模擬結(jié)果

在保持土樓夯土墻體直徑(42 m)、墻體厚度(1 m)為常數(shù)的情況下，變化墻體高度.按照墻體高度由小到大，模型為A～F(8.82～15.12 m，梯度為1.26).分析不同墻體高度下，土樓夯土結(jié)構(gòu)模型的極限狀態(tài).不同墻體高度模型變化曲線，如圖5所示.

(a) 最大應(yīng)力 (b) 最大豎向位移圖5 不同墻體高度模型最大應(yīng)力變化曲線Fig.5 Curve of maximum stress with different wall heights

由圖5可知：針對(duì)給定的其他土樓尺寸參數(shù)，土樓有一個(gè)較為合理的土樓墻體高度，在此高度下，承載能力不會(huì)太低，且極限狀態(tài)下的位移也不是最大.位移是應(yīng)變的累加，當(dāng)高度為13.86 m時(shí)，位移最大，但由于應(yīng)變不是最大，所以應(yīng)力不大.因此，土樓在選取高度時(shí)，不能一味地追求空間開闊，或者為了能有足夠大的承載能力而降低墻體高度，應(yīng)該從多方面考慮，采取科學(xué)的方式選取土樓墻體高度.

4 結(jié)論

1) 采用數(shù)據(jù)統(tǒng)計(jì)處理得到的模型尺寸，得出的模擬結(jié)果更具有普遍性和代表性.

2) 通過(guò)夯土材料的試驗(yàn)數(shù)據(jù)，結(jié)合改進(jìn)Mohr-Coulomb模型的函數(shù)，并將其與ABAQUS中的彈性模型聯(lián)用得到土樓夯土的本構(gòu)模型.

3) 模型的最大應(yīng)力和最大豎向位移隨著結(jié)構(gòu)直徑、墻體厚度和高度有明顯的變化，所以在土樓實(shí)際尺寸的選取應(yīng)該要有更科學(xué)的方法，合理的尺寸參數(shù)能使土樓的結(jié)構(gòu)性能及材料使用達(dá)到最優(yōu)，從而節(jié)約成本及人力.但由于模型較簡(jiǎn)單，缺少木結(jié)構(gòu)的協(xié)調(diào)機(jī)理，無(wú)法更具體表現(xiàn)出實(shí)際中各參數(shù)對(duì)土樓結(jié)構(gòu)性能帶來(lái)的影響.

[1] 黃漢民.福建土樓[J].福建建筑,2001(71):27-28.

[2] 吳仁偉,彭興黔.圓形土樓屋蓋風(fēng)荷載的數(shù)值模擬分析[J].華僑大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,32(4):427-432.

[3] LUO Yangjun,KANG Zhan,LI A.Structural reliability assessment based on probability and convex set mixed model[J].Computers and Structures,2009,87(21):1408-1415.

[4] 張琰鑫,童麗萍.夯土住宅結(jié)構(gòu)性能分析及加固方法[J].世界地震工程,2012,28(2):72-78.

[5] 俞鐘行.W檢驗(yàn)法程序[J].物探化探計(jì)算技術(shù),1989,11(3):267-270.

[6] 阿肯江·托乎提,亓國(guó)慶.基于Solid 65和Solid 45有限單元的素夯土墻體數(shù)值建模及計(jì)算分析[J].新疆大學(xué)學(xué)報(bào)(自然科學(xué)版),2008,25(4):492-498.

[7] 朱向榮,王金昌.ABAQUS軟件中部分土模型簡(jiǎn)介及其工程應(yīng)用[J].巖石力學(xué),2004,25(增刊1):144-148.

[8] 姚仰平,張丙印,朱俊高.土的基本特性、本構(gòu)關(guān)系及數(shù)值模擬研究綜述[J].土木工程學(xué)報(bào),2012,45(3):127-150.

[9] 趙成,阿肯江·托呼提,陳嘉,等.改性土體材料單軸受壓本構(gòu)關(guān)系研究[J].新疆大學(xué)學(xué)報(bào)(自然科學(xué)版),2010,27(7):123-126.

[10] BUI Q B,HANS S,MOREL J C,et al.First exploratory study on dynamic characteristics of rammed earth buildings[J].Engineering Structures,2011,33(12):3690-3695.

(責(zé)任編輯： 陳志賢 英文審校： 方德平)

Numerical Simulation of Circular Tulou Rammed Earth Structure Performance

YE Junjie, PENG Xingqian, SHI Weijuan

(College of Civil Engineering, Huaqiao University, Xiamen 361021, China)

By the numerical simulation method, taking the classic circular Tulou (rammed earth building) in Longyan, Fujian, such as Zhencheng Tulou, Chengqi Tulou, the circular Tulou structure model was establish. the bearing capacity of Tulou wall was investigated under the static loads. the stress and deformation characteristics of different Tulou structure models were studied for the parameters of various sizes of Tulou, the main factors influencing the stress deformation characteristics were obtained. The results show that for different wall thicknesses, diameters and heights, the stress and displacement do not change monotonously, and it has a reasonable value range.

rammed earth structures; circular Tulou; load-carrying properties; deformation characteristics; structural parameters; numerical simulation

10.11830/ISSN.1000-5013.201606008

2016-01-20

彭興黔(1959-)，男，教授，博士，主要從事防災(zāi)減災(zāi)的研究.E-mail:pxq@hqu.edu.cn.

國(guó)家科技支撐計(jì)劃項(xiàng)目(2014BAK09B01)； 國(guó)家自然科學(xué)基金資助項(xiàng)目(51178196, 51478197)； 福建省科技計(jì)劃重點(diǎn)項(xiàng)目(2013Y0067)； 福建省自然科學(xué)基金資助項(xiàng)目(2013J01193)； 華僑大學(xué)研究生科研創(chuàng)新能力培育計(jì)劃資助項(xiàng)目(2014年度)

TU 361

A

1000-5013(2016)06-0698-04