劉書賢，胡紅珍，魏曉剛，麻鳳海，薛海波

(1.遼寧工程技術(shù)大學(xué)建筑工程學(xué)院，遼寧阜新123000;2.遼寧工程技術(shù)大學(xué)土木與交通學(xué)院，遼寧阜新123000;3.大連大學(xué)建筑工程學(xué)院，遼寧大連116622;4.中建六局建設(shè)發(fā)展有限公司東北分公司，遼寧沈陽110000)

0 引言

據(jù)國有重點煤礦的統(tǒng)計分析，全國“三下”(建筑物、水體、公路鐵路)下壓煤約137.9億噸，其中建筑物下壓煤就高達(dá)87.6億噸 (王金莊，郭增長，2002)。地下煤炭的開采在地下巖土層形成了一個無實物填充承載力的空間，破壞了巖土層的完整性，導(dǎo)致煤礦采動區(qū)上覆巖體發(fā)生變形、冒落、斷裂、彎曲等一系列破壞，嚴(yán)重威脅了礦區(qū)的建筑安全。

目前國內(nèi)專家學(xué)者已經(jīng)對煤礦采動區(qū)建筑物保護(hù)問題開展了研究工作:于廣明等 (2004)利用有限元軟件研究了地震作用下煤礦采動區(qū)建筑物的動力災(zāi)變問題;滕永海和黃樂亭 (1997)根據(jù)現(xiàn)有的實際測量資料庫，指出了抗變形建筑物與其對應(yīng)地表的下沉值在宏觀上是基本一致的;譚志祥和鄧喀中 (2006)為受采動影響的建筑物保護(hù)設(shè)計提供了科學(xué)依據(jù)，總結(jié)出了不同長高比條件下建筑物變形與地表變形之間的綜合關(guān)系;劉書賢等 (2010，2013)提出了煤礦采動區(qū)建筑物抗震抗變形雙重保護(hù)支座，對開采變形和地震耦合作用下的建筑物災(zāi)變機(jī)理進(jìn)行了探討。以上研究沒有系統(tǒng)地從建筑物損傷的角度，深入探討在開采變形和地震聯(lián)合作用下建筑物的損傷演化災(zāi)變過程，缺乏對建筑物在開采變形和地震聯(lián)合作用下?lián)p傷積累的研究 (劉書賢，2014a，b，c)。

現(xiàn)有的研究對于煤礦開采變形與地震二者之間成災(zāi)機(jī)理的研究，更多的是應(yīng)用“耦合”一詞，由于煤礦采動區(qū)工程地質(zhì)環(huán)境的復(fù)雜性，煤礦采動荷載與地震荷載對建筑物所產(chǎn)生的損傷并非簡單的疊加，并且二者之間的成災(zāi)機(jī)理也不相同，所以“耦合”一詞有失偏頗。筆者針對地下煤炭開采對建筑物抗震性能的擾動研究尚不完善的現(xiàn)狀，基于損傷力學(xué)理論，建立煤礦采動損害影響下的建筑物動力損傷判據(jù)，通過建筑結(jié)構(gòu)內(nèi)力變化來探討煤礦采動與地震聯(lián)合作用下建筑物損傷演化災(zāi)變過程。

1 基于損傷因子的煤礦采動區(qū)建筑物的損傷本構(gòu)關(guān)系

煤礦采空區(qū)的地表移動變形破壞了建筑物與地基間的初始平衡狀態(tài)，在建筑物中產(chǎn)生了附加應(yīng)力，導(dǎo)致了建筑物變形破壞。

定義張量 (I-D)為由于損傷而造成的面元d A*的縮減 (李兆霞，2002)，即

式中:I為瞬時損傷;D為損傷張量。

煤礦采動對建筑物的損傷破壞表現(xiàn)在由于建筑自身重力的作用下的徐變損傷過程，而徐變過程中的總損傷包括瞬時損傷和徐變損傷，因此

當(dāng)應(yīng)力水平低于長期強(qiáng)度時，σ/fc'＜0.8徐變最終處于穩(wěn)定階段，對于給定的常應(yīng)力，徐變速率和損傷值都趨于定值。這時，DB=0，Dcr=DB。

損傷變化率關(guān)于損傷是衰減的，并可認(rèn)為是損傷和應(yīng)力的顯函數(shù)，因此=-γDcr+f(σ).

邊界條件為

假設(shè):混凝土是均質(zhì)、彈性和各向同性的 (李杰，楊衛(wèi)忠，2009)并且徐變損傷是各向同性的，可選取標(biāo)量D作為損傷的度量;對于和時間相關(guān)的應(yīng)力，有效應(yīng)力相應(yīng)地為

將Lemaitre提出的應(yīng)變等效假設(shè)推廣到和時間相關(guān)的徐變問題，推導(dǎo)出彈性徐變耦合損傷的徐變損傷本構(gòu)方程:

式中:ε(t)為鋼筋混凝土的應(yīng)力;σ(t)為鋼筋混凝土的應(yīng)變;E(t)為彈性模量;c(t，τ)為齡期τ下的徐變度;［1-D(t)］為有效應(yīng)力。

由公式 (3)可知:煤礦采動損害影響下建筑物的次生損傷與其應(yīng)力成反比 (克拉夫，彭津，2006)。

耦合性即塊間聯(lián)系，是一種量度，泛指系統(tǒng)結(jié)構(gòu)中各模塊間相互聯(lián)系緊密程度;而聯(lián)合則是一種新的數(shù)據(jù)類型，它是一種特殊形式的變量。由于煤礦采動區(qū)工程地質(zhì)環(huán)境的復(fù)雜性，煤礦采動荷載與地震荷載對建筑物所產(chǎn)生的損傷并非簡單的疊加，根據(jù)混凝土損傷本構(gòu)方程，由于二者之間的成災(zāi)機(jī)理也不相同，所以“耦合”一詞有失偏頗，“聯(lián)合”一詞更加準(zhǔn)確。

2 地震作用下煤礦采動區(qū)建筑物損傷演化的有限元分析模型

2.1 有限元模型的選取

為了探討地震發(fā)生時煤礦采動建筑的次生損傷累積演化過程 (楊偉，歐進(jìn)萍，2009;張春禮，2011)，重點分析煤礦采動與地震聯(lián)合作用下建筑物的損傷演化成災(zāi)機(jī)制，尤其是建筑物從彈性階段逐步發(fā)展到塑性鉸出現(xiàn)的塑性損傷破壞過程，故選取唐山某礦區(qū)一框架結(jié)構(gòu)的辦公樓 (余紅發(fā)等，2006)。

唐山某礦區(qū)地質(zhì)條件良好，抗震設(shè)防烈度為Ⅶ度，該現(xiàn)澆鋼筋混凝土框架結(jié)構(gòu)辦公樓，樓層層高統(tǒng)一為3.3 m，，縱向4跨，橫向2跨，跨度均為6 m，該建筑物總高19.8 m，長24 m，寬12 m。柱子截面尺寸均為600 mm×600 mm，梁截面均為300 mm×650 mm，樓板厚度為120 mm，梁柱均采用C45混凝土，彈性模量E=32.5 GPa，泊松比v=0.2，密度為2 450 kg/m3(圖1)。

該礦區(qū)由于地下煤炭開采致使地表發(fā)生移動和變形 (張永波，2006)，根據(jù)現(xiàn)場監(jiān)測站的實際監(jiān)測數(shù)據(jù)，由于采煤導(dǎo)致辦公樓產(chǎn)生了2 mm/m的傾斜 (顏榮貴，1995)。

2.2 地震波的選取

根據(jù)《建筑抗震設(shè)計規(guī)范》(GB50011—2010)中的相關(guān)內(nèi)容，要求選取地震波要滿足峰值、頻譜特性、地震持時的要求，結(jié)合該辦公樓所處區(qū)域條件及結(jié)構(gòu)的自振周期，為了能夠分析建筑物在采動與地震聯(lián)合作用下的耗散能量的變化過程，需要選取較長時間的強(qiáng)烈地震波。本文選用了15 s的唐山地震波，其地震等級為7級，方向為南北向，最大加速度為0.24 m·s-2(圖2)。

3 有限元數(shù)值計算結(jié)果分析

煤礦采動區(qū)的開采變形是隨著開采過程逐漸發(fā)展的，是長期緩慢發(fā)展的過程，不僅發(fā)生在煤炭開采過程中，而且在開采完成后仍然不斷發(fā)展，對煤礦采動區(qū)的土地及建筑物損害嚴(yán)重;煤礦采動區(qū)的地震則是一個瞬間劇烈破壞的過程，在短時間內(nèi)能量急劇釋放，地表的移動變形強(qiáng)烈，當(dāng)?shù)卣鸾Y(jié)束時地表移動變形幾乎完全消失，對建筑物的損傷破壞也僅存在于地震發(fā)生的過程中。

通過有限元數(shù)值計算，所得到的不同工況下建筑物的地震動力響應(yīng)如圖3所示。工況一為不考慮煤礦采動作用僅考慮地震荷載，工況二為考慮煤礦采動與地震聯(lián)合作用。

分析圖3a可知:工況1建筑物層間位移變化曲線比較平滑，頂層最大位移為0.045 m;工況2建筑物層間位移變化曲線比較層間位移變化急劇，頂層最大值達(dá)到0.075 m，是工況1的1.67倍。由此可以判斷:煤礦采動作用對建筑物所產(chǎn)生的次生損傷對其地震動力響應(yīng)影響較大，地震發(fā)生時，建筑物的次生損傷進(jìn)一步發(fā)展演化。地震作用下煤礦采動作用加劇了結(jié)構(gòu)的動力反應(yīng)，嚴(yán)重降低了建筑物的抗震安全性能。

分析圖3b可知:工況1建筑物底層的最大加速度為1.225 m·s-2，頂層的最大加速度為7.163 m·s-2;工況2建筑物的底層最大加速度為4.788 m·s-2，頂層的頂層為17.163 m·s-2;分別是工況1底層的加速度3.91倍，頂層的加速度2.5倍。對于煤礦采動與地震聯(lián)合作用下建筑的最大加速度曲線在五層六層發(fā)生突變現(xiàn)象 (加速度變化較大)，說明在煤礦采動損害影響下建筑物產(chǎn)生損傷，其強(qiáng)度和剛度都有一定程度的退化，地震發(fā)生時建筑物鞭梢效應(yīng)明顯，導(dǎo)致其動力響應(yīng)較大。

分析圖3c可知:不考慮煤礦采動對建筑物所產(chǎn)生的次生損傷，地震作用下建筑物的最大層間剪力曲線變化較為光滑均勻;而煤礦采動與地震聯(lián)合作用下建筑物的最大層間剪力曲線變化趨勢非常陡峭，而且其轉(zhuǎn)折點非常明顯 (轉(zhuǎn)折點主要集中在樓層的第二層和第五層)，說明這兩層為建筑物的結(jié)構(gòu)薄弱處，煤礦采動作用對建筑物的抗震能力影響較大，需要根據(jù)煤礦采動對建筑物所產(chǎn)生的次生損傷予以控制，開展建筑物的煤礦采動損傷評估以保證其安全性。

分析圖3d可知:在兩種不同荷載工況影響下的建筑物最大層間軸力變化曲線均較為光滑，最大軸力在底層達(dá)到最大 (與最大剪力相似)，說明地震作用下建筑物的底層位置 (與建筑物是否發(fā)生次生損傷無關(guān))需要予以重視，同時由于地震作用力是按照建筑物的剛度進(jìn)行分配的，建筑物底層一旦發(fā)生損傷破壞，建筑物的質(zhì)量中心和剛度中心不一致會導(dǎo)致建筑物產(chǎn)生扭轉(zhuǎn)效應(yīng)，加劇建筑物的地震動力響應(yīng)，加劇建筑的倒塌破壞。

分析圖3e可知:隨著建筑物高度的增加，建筑物的彎矩峰值整體上呈現(xiàn)出減小的趨勢，但是在建筑物的第二層卻達(dá)到最大。結(jié)合建筑的層間剪力峰值、層間軸力峰值、位移峰值和加速度峰值可知:在建筑物的第二層處為該結(jié)構(gòu)的特殊位置 (薄弱層位置)。在煤礦采動損害影響下，薄弱層位置極易出現(xiàn)大量的塑性鉸，尤其是在建筑物底部梁柱節(jié)點處形成塑性鉸，構(gòu)件強(qiáng)度減小、剛度退化，嚴(yán)重削弱了建筑物的抗震性能，在煤礦采動建筑的損傷安全評估工作中必須對結(jié)構(gòu)的薄弱層位置予以重視，并采取合理的抗震加固措施。

4 結(jié)論

本文從損傷角度出發(fā)，基于損傷結(jié)構(gòu)理論，利用有限元計算軟件ANSYS建立了煤炭開采區(qū)建筑物在地震作用和采動地震聯(lián)合作用下的計算模型，得到各工況下的最大加速度，層間剪力軸力等結(jié)果，分析得到如下研究成果:

(1)由于煤礦采動區(qū)工程地質(zhì)環(huán)境的復(fù)雜性，煤礦采動荷載與地震荷載對建筑物所產(chǎn)生的損傷并非簡單的疊加，并且二者之間的成災(zāi)機(jī)理也不相同，所以“耦合”一詞有失偏頗，“聯(lián)合作用”能夠更好的描述煤礦采動與地震對建筑的破壞。

(2)基于損傷因子所建立的煤礦采動區(qū)建筑物的損傷本構(gòu)關(guān)系為

該本構(gòu)關(guān)系較好的反應(yīng)了煤礦采動損害影響下建筑物的次生損傷與其應(yīng)力成反比，揭示了煤礦采動作用降低建筑物抗震性能的本質(zhì)原因。

(3)地震作用下煤礦采動對建筑物所產(chǎn)生的次生損傷不斷發(fā)展演化，導(dǎo)致建筑物的薄弱層位置形成塑性鉸，嚴(yán)重削弱了建筑物的抗震性能，同時煤礦采動也加劇了建筑的地震動力響應(yīng)，對與煤礦采動建筑必須對煤礦開采變形與地震荷載的成災(zāi)機(jī)理，以及二者的聯(lián)合作用成災(zāi)機(jī)制予以重視，開展煤礦采動建筑物抗開采變形隔震保護(hù)體系相關(guān)的研究工作。

克拉夫R.，彭津J..2006.結(jié)構(gòu)動力學(xué)［M］.北京:高等教育出版社.

李杰，楊衛(wèi)忠.2009.混凝土彈塑性隨機(jī)損傷本構(gòu)關(guān)系研究［J］.土木工程學(xué)報，42(2):31-38.

李兆霞.2002.損傷力學(xué)及其應(yīng)用［M］.北京:科學(xué)出版社.

劉書賢，郭濤，魏曉剛，等.2014b.地震作用下煤礦開采損傷建筑的能量耗散演化致災(zāi)分析［J］.地震研究，37(3):442-449.

劉書賢，王春麗，魏曉剛，等.2014a.煤礦采空區(qū)的地震動力響應(yīng)及其對地表的影響［J］.地震研究，37(4):642-647.

劉書賢，魏曉剛，麻風(fēng)海，等.2013.煤礦采動與地震耦合作用下建筑物災(zāi)變分析［J］.中國礦業(yè)大學(xué)學(xué)報，42(4):526-534.

劉書賢，魏曉剛，魏亞強(qiáng).2010.采動區(qū)建筑物抗震、抗變形雙重保護(hù)作用機(jī)理分析［J］.防災(zāi)減災(zāi)工程學(xué)報，30(增刊):305-308.

劉書賢，魏曉剛，張弛，等.2014c.基于隔震技術(shù)的橋梁結(jié)構(gòu)抗采動變形抗地震保護(hù)［J］.地震研究，37(1):86-93.

譚志祥，鄧喀中.2006.建筑物下采煤研究進(jìn)展［J］.遼寧工程技術(shù)大學(xué)學(xué)報，25(4):485-488.

滕永海，黃樂亭.1997.抗變形建筑物的采動變形特點［J］.東北煤炭技術(shù)，3(3):13-16.

王金莊，郭增長.2002.我國村莊下采煤的回顧與展望［J］.中國煤炭，28(5):28-31.

顏榮貴.1995.地基開采及其地表建筑［M］.北京:中國建筑工業(yè)出版社.

楊偉，歐進(jìn)萍.2009.基于能量原理的Park＆Ang損傷模型簡化計算方法［J］.地震工程與工程振動，29(2):159-165.

于廣明，董春勝，潘永戰(zhàn)，等.2004.混凝土的分形性及其單軸應(yīng)力下裂紋演化的混沌效應(yīng)［J］.青島建筑工程學(xué)院學(xué)報，25(1):1-6.

余紅發(fā)，張偉，金祖權(quán)，等.2006.土木工程結(jié)構(gòu)混凝土壽命預(yù)測的損傷演化方程［J］.東南大學(xué)學(xué)報(自然科學(xué)版)，36(增刊):216-220.

張春禮.2011.采動與地震共同作用下建筑物的破壞過程研究［D］.青島:青島理工大學(xué).

張永波.2006.老采空區(qū)建筑地基穩(wěn)定性評價理論與方法［M］.北京:中國建筑工業(yè)出版社.

GB50011—2010，建筑抗震設(shè)計規(guī)范［S］.