李海健

龍口市住房保障和房產(chǎn)管理辦公室 山東 煙臺 265701

前言

砌體是由塊體和砂漿砌筑而成的墻或柱[1]，包括磚砌體、砌塊砌體、石砌體和墻板砌體，對于一般的建筑來說，砌體的重量占到整個建筑物自重的一半。砌體材料已經(jīng)采用工業(yè)廢料為原料的綠色建筑材料的砌塊，建筑物的柱子作為建筑物最主要的受力構(gòu)件，它包括磚結(jié)構(gòu)、石結(jié)構(gòu)和其他材料的砌塊結(jié)構(gòu)。

墻體裂縫會嚴(yán)重危害建筑物安全[2]，要對裂縫產(chǎn)生的原因進(jìn)行分析和研究。一般情況下裂縫的產(chǎn)生主要是由于地基沉降產(chǎn)生尤其是不均勻的沉降，不均勻沉降在砌體建筑物引起裂縫事故占有很高的比例[3]，所以研究不均勻沉降和砌體之間的相互影響和相互作用是很有必要和有意義的[4]。

1 有限元模型建立

利用ANSYS有限元對砌體進(jìn)行分析，將混凝土和其內(nèi)部鋼筋結(jié)合成一個整體相互間不產(chǎn)生相對位移，組成同一個單元；在計算時，分別對混凝土和鋼筋按照各自的單元剛度矩陣進(jìn)行，圖1為砌體結(jié)構(gòu)模型。

圖1 砌體結(jié)構(gòu)模型

在ANSYS軟件中，對于墻體材料可采用當(dāng)采用SOLID65三維實(shí)體單元進(jìn)行模擬[5]，這個單元可以模擬變形和應(yīng)變，混凝土的開裂和破碎的模擬可以采用SOLID65和CONCRETE相結(jié)合來進(jìn)行。

2 有限元模型分析

ANSYS軟件是通過各種求解器來計算相互關(guān)聯(lián)的線型方程來分析和判斷建筑結(jié)構(gòu)工程的變化情況，結(jié)構(gòu)是屬于非線性問題，解決非線性問題通常采用線性方程對非線性進(jìn)行擬合和逼近，在這樣的求解過程中會存在誤差的積累，可能最終導(dǎo)致計算無法收斂，為了解決這個問題，ANSYS采用了一種迭代方法，也叫增量迭代算法，可以使每一個增量階段計算結(jié)果先進(jìn)行收斂，這樣對于整個計算過程來說可以起到消除誤差的作用。本文以建筑物結(jié)構(gòu)兩端發(fā)生沉降作用的影響進(jìn)行分析。

建筑物兩端沉降產(chǎn)生的原因主要是由于兩端地基結(jié)構(gòu)土質(zhì)強(qiáng)度不均勻，強(qiáng)度較低，中間地基強(qiáng)度較好，從而產(chǎn)生兩端沉降。對于這種情況，在ANSYS分析中可以改變建筑物模型兩端的彈簧剛度進(jìn)行模擬兩端沉降現(xiàn)象，可以分析計算建筑物的整體變形情況。

圖2 兩端沉降作用下模型整體位移圖

圖3 兩端沉降作用下模型基礎(chǔ)位移圖

圖2,3是ANSYS分析結(jié)果，可以看出在建筑物的整體變形中，兩端下沉，沉降程度從兩側(cè)向中間逐步減小，在壓力作用下，受力最大處位于建筑物底部中間部位，產(chǎn)生了強(qiáng)烈的壓應(yīng)力，建筑物下部產(chǎn)生凸起。整體來說，受力幾中的部分在建筑物中部，圖4為模型的一面墻體。取一面墻體進(jìn)行分析，得到了在兩端沉降作用下墻面的應(yīng)力云圖，如圖5,6所示。

圖4 一面墻體單元

圖5 兩端沉降作用下拉應(yīng)力云圖

圖6 兩端沉降作用下應(yīng)力強(qiáng)度云圖

云圖中不同的顏色代表了不同的應(yīng)力大小，很明顯應(yīng)力集中的地方出現(xiàn)在了建筑物門窗的角部，對整體墻面來說還出現(xiàn)的左右對稱情況。隨著建筑物樓層的升高，應(yīng)力開始減小，應(yīng)力最大位置為最底層的窗戶角部，這是因?yàn)榈撞吭诔两底饔孟伦冃巫畲?。同樣，對于切?yīng)力的變化在圖6有所顯示，是墻體受到剪力的作用，應(yīng)力強(qiáng)度也是集中的地方出現(xiàn)在了建筑物門窗的角部，并隨著建筑物樓層的升高，應(yīng)力開始減小，模擬計算情況符合實(shí)際現(xiàn)象。不管是拉應(yīng)力還是切應(yīng)力，都是上層平緩，底層劇烈。通過模型可以看出，兩端沉降位移大于中間部分沉降位移量，使整個建筑物承受負(fù)彎矩同時承受的剪切的作用。墻體裂縫就是在應(yīng)力集中的地方產(chǎn)生，如圖7,8所示。建筑物門窗開口部分，對整個墻體的強(qiáng)度起到了削弱的作用，同時在開口處的角部存在應(yīng)力集中，當(dāng)應(yīng)力大于墻體材料的強(qiáng)度時，就會產(chǎn)生裂縫，同時裂縫會延展，向著沉降大的方向進(jìn)行。

圖7 兩端沉降作用下最終裂縫圖

圖8 兩端沉降作用下變形圖

為了更好地分析建筑物墻體在兩端沉降的作用下裂紋的產(chǎn)生和相關(guān)應(yīng)力的變化情況，利用后置處理器可以生成沉降位移的變化曲線，整個建筑物是對稱結(jié)構(gòu)，從上面的云圖分析也可以看出受力情況和變形情況也是呈現(xiàn)出對稱結(jié)構(gòu)，因此在繪制位移曲線時，可以取建筑物的一半墻體進(jìn)行分析。先選取中間墻體部分一系列的連續(xù)節(jié)點(diǎn)，定義其為沉降方向的位移，通過增大位移量，其墻體受力也會逐漸增大，可以得到最終產(chǎn)生裂縫時的沉降量，如圖9所示。同樣方法，選取每層的底部變化情況進(jìn)行分析，逐步加大位移量從而產(chǎn)生裂縫，可以得到每層的沉降位移變化曲線，能夠直觀的得到沉降位移變化和裂縫產(chǎn)生的相互關(guān)系。

圖9 兩端沉降作用下半墻的地基沉降量

通過ANSYS分析計算后，可以得到建筑物地基的變化情況，某些部位的傾斜值超過了0.004，這個數(shù)值明顯是不符合國標(biāo)《建筑地基基礎(chǔ)設(shè)計規(guī)范（GB50007-2011）》的要求，國標(biāo)規(guī)定承載砌體建筑物的地基局部的傾斜值不能超過0.003，通過分析，這種兩端沉降已經(jīng)破壞了建筑物的結(jié)構(gòu)。

在ANSYS對建筑物模型的分析中，是逐步增加載荷實(shí)現(xiàn)加載過程，整個分析計算過程一共是進(jìn)行了152步的加載過程，第一條裂縫的出現(xiàn)是在第11步載荷加載，裂縫出現(xiàn)在第一層中間部分的開門洞口的底邊角處。在第58步加載中，建筑物的頂層的門開口處出現(xiàn)了第二條裂縫，同時前面出現(xiàn)的第一條裂縫在不斷地擴(kuò)展。到了第96步，建筑物裂縫開始向二層的洞口發(fā)展，首先是出現(xiàn)在洞口的邊角處，同時可以看出整體墻面變形加大，第110步，建筑物頂層產(chǎn)生的裂縫開始和二層產(chǎn)生的裂縫開始對接，隨著加載步的增多，各樓層的裂紋不斷擴(kuò)展連接，發(fā)生劇烈變化，最嚴(yán)重的是底層裂紋。

圖10～15就是裂縫變化情況圖：

圖10 第11步加載（0.00335）

圖11 第58步加載（0.00346）

圖12 第59步加載（0.00351）

圖13 第60步加載（0.00357）

圖14 第96步加載（0.00378）

圖15 第110步加載（0.00385）

3 結(jié)束語

建立了墻體模型，運(yùn)用ANSYS有限元軟件采用非線性理論分析研究了建筑物沉降變化對墻面裂縫產(chǎn)生的影響和變化規(guī)律，利用非線性有限元，選用對SOLID65和CONCRETE相結(jié)合來進(jìn)行混凝土的開裂和破碎的模擬是可行的。

參考文獻(xiàn)

[1] 徐占發(fā),許大江.砌體結(jié)構(gòu)[M].北京,中國建材工業(yè)出版社,2010:78-80.

[2] 劉斌.地基基礎(chǔ)不均勻沉降的原因分析及應(yīng)對措施[J].河北建筑工程學(xué)院學(xué)報,2013,31(1):20-21.

[3] 潘明楊,韓建.地基不均勻沉降引起裂縫預(yù)防與治理[J].農(nóng)村經(jīng)濟(jì)與技術(shù),2012,21(8):152-154.

[4] 王文平,洪麗仙.砌體結(jié)構(gòu)優(yōu)缺點(diǎn)及設(shè)計要點(diǎn)的分析[J].建筑技術(shù)研究,2012,(10):54-56.

[5] 司炳君,孫治國,艾慶華.Solid65單元在混凝土結(jié)構(gòu)有限元分析中的應(yīng)用[J].工業(yè)建筑,2007,37(1):87-92.