徐龍河，楊冬玲，李忠獻(xiàn)

(1.北京交通大學(xué) 土木建筑工程學(xué)院，北京 100044;2.天津大學(xué) 建筑工程學(xué)院，天津 300072)

2008年5月12日，我國(guó)四川省汶川縣發(fā)生8.0級(jí)特大地震，截至6月18日12時(shí)，已造成69176人遇難，受災(zāi)人數(shù)已達(dá)4624萬(wàn)，初步統(tǒng)計(jì)損失已達(dá)5000多億人民幣［1]。造成這些巨大經(jīng)濟(jì)損失和大量人員傷亡的根本原因是地震導(dǎo)致大量房屋倒塌，因此房屋安全問(wèn)題必須引起高度的重視。如何正確的評(píng)價(jià)結(jié)構(gòu)的抗震能力及地震過(guò)后房屋的剩余承載力，需要建立一個(gè)合理的損傷模型以定量計(jì)算結(jié)構(gòu)在地震中的損傷指標(biāo)。損傷指標(biāo)是用來(lái)描述結(jié)構(gòu)、構(gòu)件或材料受損程度的變量，一般定義為結(jié)構(gòu)或構(gòu)件反應(yīng)歷程中某一指標(biāo)累計(jì)量與相應(yīng)指標(biāo)極限允許量之比，通常用D表示，并且有以下兩個(gè)特征［2－4]:①損傷指標(biāo) D 的范圍為［0，1] ，當(dāng) D=0時(shí)，對(duì)應(yīng)結(jié)構(gòu)的完好狀態(tài);當(dāng)D=1時(shí)，意味著結(jié)構(gòu)或構(gòu)件完全破壞;②損傷指標(biāo)D為單調(diào)遞增的函數(shù)，即結(jié)構(gòu)的損傷向著增大的方向發(fā)展且不可逆。

結(jié)構(gòu)在地震力的作用下產(chǎn)生不同程度的損傷，并且隨著荷載循環(huán)次數(shù)的增加而不斷地累積，最終導(dǎo)致結(jié)構(gòu)破壞。合理的損傷模型應(yīng)能充分體現(xiàn)出大的地震幅值和重復(fù)加載效應(yīng)聯(lián)合作用導(dǎo)致結(jié)構(gòu)破壞這一事實(shí)。比較有代表性的是Park和Wen［5]提出的最大反應(yīng)變形和累積耗能線性組合的地震損傷模型:

式中，δm和δu分別為結(jié)構(gòu)或構(gòu)件的最大位移和極限位移，d E為滯回耗能的增量，Qy為結(jié)構(gòu)或構(gòu)件的屈服強(qiáng)度，β為能量影響系數(shù)。

變形和能量雙參數(shù)損傷模型較好的體現(xiàn)了地震動(dòng)三要素(幅值，持時(shí)和頻率)對(duì)結(jié)構(gòu)破壞的影響，被學(xué)者廣泛接受和應(yīng)用。很多學(xué)者在Park-Ang模型的基礎(chǔ)上進(jìn)行改進(jìn)，提出了很多不同形式的雙參數(shù)損傷模型。歐進(jìn)萍等［6]同樣采用最大變形與累積耗能的組合，提出了一種適用于鋼結(jié)構(gòu)的雙參數(shù)損傷模型;牛狄濤等［7]指出park模型中變形和能量進(jìn)行線性組合具有不合理性，提出了變形和能量的非線性組合的雙參數(shù)損傷模型，并基于實(shí)際鋼筋混凝土結(jié)構(gòu)的破損分析確定了模型的參數(shù);傅劍平等［8]根據(jù)所收集的近十年來(lái)國(guó)內(nèi)外完成鋼筋混凝土柱的試驗(yàn)結(jié)果，對(duì)Park－Ang損傷模型進(jìn)行了識(shí)別，對(duì)位移貢獻(xiàn)和能量項(xiàng)貢獻(xiàn)進(jìn)行了正確的評(píng)估，提出了一個(gè)更為精確的雙參數(shù)損傷模型，但暫時(shí)只限于柱類(lèi)構(gòu)件。

以上的損傷模型是從構(gòu)件的角度，對(duì)Park模型進(jìn)行改進(jìn)得到的，同時(shí)一些學(xué)者突破park模型的形式，從材料的角度提出的損傷模型也被廣泛應(yīng)用。韋未，李同春［9－10]將多軸應(yīng)力狀態(tài)下的應(yīng)變退化得到等效應(yīng)變，提出了建立在應(yīng)變空間的混凝土四參數(shù)損傷模型。該損傷模型不僅能正確的反映混凝土單軸應(yīng)力狀態(tài)下的強(qiáng)度，而且很好的反映多軸應(yīng)力狀態(tài)下混凝土的強(qiáng)度。沈祖炎［11]采用塑性應(yīng)變作為依據(jù)，并考慮滯回耗能，提出一個(gè)適用于鋼構(gòu)件的損傷模型，并通過(guò)低周疲勞試驗(yàn)確定了其中的參數(shù)。宋振森［12]把一維應(yīng)力的損傷理論應(yīng)用到三維應(yīng)力，引入等效塑性應(yīng)變來(lái)表征損傷模型，得出了三軸應(yīng)力下構(gòu)件的損傷模型。

空間結(jié)構(gòu)受多軸應(yīng)力的作用，將多軸應(yīng)力退化得到的等效應(yīng)力和等效應(yīng)變可以表征空間結(jié)構(gòu)的受力狀態(tài)。本文從材料的角度沿用等效應(yīng)變的思想，并與比能線性組合，提出的損傷模型可以表征三軸應(yīng)力下結(jié)構(gòu)的損傷變化。

1 應(yīng)變和比能雙控?fù)p傷模型

1.1 損傷模型

雙參數(shù)損傷模型以變形和能量組合的方式很好的體現(xiàn)了結(jié)構(gòu)的破壞機(jī)理，本文在此基礎(chǔ)上考慮結(jié)構(gòu)空間三向應(yīng)力的作用，建立了基于等效塑性應(yīng)變和比能雙控的雙參數(shù)損傷模型:

式中:dj為第j個(gè)構(gòu)件的損傷指標(biāo)，pi為第i時(shí)刻第j構(gòu)件等效塑性應(yīng)變;pu為材料的極限等效塑性應(yīng)變;u為循環(huán)過(guò)程中產(chǎn)生的累積比能耗散;σy為材料的屈服應(yīng)力;β為循環(huán)荷載影響系數(shù)。

其中，鋼材的等效塑性應(yīng)變?yōu)?

式中:σi和σi+1分別為第i時(shí)刻和i+1時(shí)刻的應(yīng)力，εi和εi+1分別為第i時(shí)刻和i+1時(shí)刻的應(yīng)變。

根據(jù) Simo和 Ju［12－13]的分析，延性材料的損傷演化過(guò)程可以視為各向同性的，故單向拉伸及循環(huán)拉壓試驗(yàn)得到的結(jié)論在三維應(yīng)力狀態(tài)下仍能成立。把應(yīng)力應(yīng)變用等效應(yīng)力和等效應(yīng)變代替，式(4)可以表征三維應(yīng)力下的累積比能耗散，即:

1.2 參數(shù)確定

文獻(xiàn)［4] 指出從Park模型對(duì)破壞構(gòu)件的計(jì)算結(jié)果中不難看出，能量項(xiàng)在其中占了很小的份額，對(duì)于大多數(shù)構(gòu)件來(lái)言，能量項(xiàng)并沒(méi)有小到可以忽略不計(jì)的地步。因此仍有必要在要求較高的抗震設(shè)計(jì)或評(píng)價(jià)中使用相對(duì)更準(zhǔn)確的雙參數(shù)準(zhǔn)則。本文提出的損傷模型以等效塑性應(yīng)變和比能進(jìn)行線性組合，β為循環(huán)荷載影響系數(shù)，也可稱(chēng)為應(yīng)變項(xiàng)和比能項(xiàng)的組合系數(shù)，需要通過(guò)試驗(yàn)確定。

Park等人做了261個(gè)試驗(yàn)認(rèn)為混凝土循環(huán)荷載影響系數(shù)β受剪跨比、軸壓比、配筋率和配箍率的影響，一般在 0.3 ～ 0.12 之間，常取 0.15 左右［2，5]。文獻(xiàn)［14] 中對(duì)歐進(jìn)萍等人提出的鋼結(jié)構(gòu)損傷模型中的參數(shù)進(jìn)行了定義，重要結(jié)構(gòu)取β=2.0，一般結(jié)構(gòu)取為1.0，對(duì)位移項(xiàng)和能量項(xiàng)定義了同樣的組合系數(shù)。這是因?yàn)殇摻Y(jié)構(gòu)的耗能能力比混凝土的強(qiáng)，其損傷模型中的能量項(xiàng)占了損傷指標(biāo)D比較大的份額。

循環(huán)荷載作用下，結(jié)構(gòu)的破壞是由首超破壞和累積破壞相互影響的，在鋼結(jié)構(gòu)中，其耗能能力更強(qiáng)，故循環(huán)荷載對(duì)損傷模型的影響更大，所以本文循環(huán)荷載影響系數(shù)β取為0.3。

1.3 結(jié)構(gòu)層損傷模型

根據(jù)構(gòu)件損傷模型定義如下結(jié)構(gòu)層損傷模型:

式中:Dg為結(jié)構(gòu)層損傷指標(biāo)，Wj為j構(gòu)件的損傷加權(quán)值，體現(xiàn)了第j個(gè)構(gòu)件在整個(gè)結(jié)構(gòu)層中的重要程度。

2 結(jié)構(gòu)損傷程度界定標(biāo)準(zhǔn)

根據(jù)震后建筑或構(gòu)件的破損程度和破壞修復(fù)的難易程度，將震害劃分為五個(gè)等級(jí)，即基本完好，輕微破壞，中等破壞，嚴(yán)重破壞和倒塌。由于綜合損傷模型可以很好的反映變形和能量對(duì)結(jié)構(gòu)的影響，國(guó)內(nèi)外學(xué)者在各自所提出的損傷模型的基礎(chǔ)上進(jìn)一步劃分結(jié)構(gòu)的震害等級(jí)，確定了結(jié)構(gòu)的模糊安全準(zhǔn)則。本文的損傷模型是從等效塑性應(yīng)變角度提出的，當(dāng)結(jié)構(gòu)處在彈性階段時(shí)，塑性應(yīng)變?yōu)?，此時(shí)的損傷指數(shù)為0，在此不考慮結(jié)構(gòu)在彈性階段的損傷，即假定結(jié)構(gòu)處在彈性變形階段時(shí)是完好的。參考《中國(guó)地震烈度表》中給出的震害指數(shù)及文獻(xiàn)［17] ，并結(jié)合本損傷模型的特點(diǎn)，定義了對(duì)應(yīng)結(jié)構(gòu)不同破壞等級(jí)的損傷指數(shù)范圍。表1為國(guó)內(nèi)外學(xué)者損傷模型和本文損傷模型不同的震害等級(jí)對(duì)應(yīng)的損傷指標(biāo)范圍。

表1 不同震害等級(jí)對(duì)應(yīng)的損傷指標(biāo)范圍Tab.1 Ranges of damage index corresponding to different damage grades

3 數(shù)值算例

3.1 有限元模型

為了驗(yàn)證本文所提出的損傷模型的正確性和實(shí)用性，采用 9層 Benchmark結(jié)構(gòu)模型［18]，對(duì)其輸入 El Centro波，進(jìn)行倒塌數(shù)值模擬分析。該結(jié)構(gòu)平面尺寸為45.73 m ×45.73 m，X、Y 方向各五跨，每跨寬 9.15 m;原結(jié)構(gòu)地上九層，地下一層，本文只取地上部分進(jìn)行建模。首層層高為5.49 m，標(biāo)準(zhǔn)層層高為3.96 m。結(jié)構(gòu)的總體質(zhì)量為9×106kg。模型立面圖如圖1所示。梁柱均采用H型鋼，截面尺寸如表2所示。柱的屈服強(qiáng)度為345 MPa，梁的屈服強(qiáng)度為248 MPa，采用與應(yīng)變率有關(guān)的塑性隨動(dòng)強(qiáng)化材料模型，通過(guò)Cowper—Symonds方程考慮應(yīng)變率對(duì)屈服應(yīng)力的影響［19]。采用有限元軟件LS－DYNA建模分析，利用空間梁?jiǎn)卧M梁柱構(gòu)件，板單元模擬樓板，有限元模型如圖2所示。分別輸入2倍，3倍，5倍，6倍，7倍的雙向 El Centro波，持續(xù)時(shí)間為20s。El Centro加速度時(shí)程曲線如圖3所示。

圖1 9層benchmark模型結(jié)構(gòu)立面Fig.1 Elevation of 9-story benchmark model structure

圖2 有限元模型Fig.2 Finite element model

圖3 El Centro波加速度時(shí)程曲線Fig.3 Time histories of El Centro earthquake records

表2 梁柱截面尺寸Tab.2 Sizes of beam and column

3.2 結(jié)果分析

圖4為5倍El Centro波下柱子損傷指標(biāo)時(shí)程曲線。對(duì)比各個(gè)時(shí)刻角柱，邊柱和中柱的損傷演化值，可以看出不同位置的柱子損傷發(fā)展過(guò)程有相同的規(guī)律。在2 s之前柱子處在彈性階段，此時(shí)的損傷為0;在2.6 s和11.4 s時(shí)，損傷突然增大，因?yàn)槟媳毕蚝蜄|西向的El Centro地震波加速度值在這兩個(gè)時(shí)刻分別達(dá)到了最大值;隨著地震波強(qiáng)度的加大和時(shí)間的延長(zhǎng)，損傷逐漸累積，結(jié)構(gòu)的損傷向著增大的方向發(fā)展且損傷不可逆，曲線單調(diào)遞增，這與損傷指標(biāo)特征的第二項(xiàng)是一致的。對(duì)比同時(shí)刻不同位置柱子的損傷可以看出，角柱的損傷大于邊柱，邊柱的損傷大于中柱，與結(jié)構(gòu)在雙向地震中角柱受力最大，其次為邊柱，中柱受力相對(duì)較小相一致。

圖4 5.0倍El Centro地震動(dòng)下柱子損傷發(fā)展過(guò)程Fig.4 Damage evolution process of columns under the excitation of 5.0 times El Centro earthquake

表3是地震波持續(xù)時(shí)間內(nèi)4個(gè)主要時(shí)刻柱子的損傷值和一至三層層間損傷值的變化率。從表上可以清晰的看出，同一時(shí)刻底層柱的損傷值較大，層數(shù)越高損傷值越小?？蚣芙Y(jié)構(gòu)以剪切變形為主，并且首層較高，形成薄弱層，在地震波作用結(jié)束時(shí)，底層柱破壞最為嚴(yán)重，損傷指數(shù)達(dá)到了0.450 7，對(duì)照不同震害等級(jí)對(duì)應(yīng)的損傷指標(biāo)范圍表1，此時(shí)柱子已經(jīng)嚴(yán)重破壞。從柱層間損傷值變化率上來(lái)看，當(dāng)構(gòu)件損傷值達(dá)到最大的時(shí)刻，3層相對(duì)2層的損傷值變化率小于2層相對(duì)1層的，說(shuō)明框架結(jié)構(gòu)層數(shù)越高柱損傷程度變化越小。

表3 主要時(shí)間步柱的損傷值和層間損傷值變化率Tab.3 Column damage value and inter-story damage variation rate at main time steps

圖5為5倍、6倍和7倍El Centro地震輸入時(shí)結(jié)構(gòu)首層損傷指標(biāo)變化曲線。從中看出，結(jié)構(gòu)層的損傷與柱子的損傷發(fā)展有相同的規(guī)律，即結(jié)構(gòu)先是處在彈性階段，隨著地震強(qiáng)度的加大和地震波持續(xù)時(shí)間的增長(zhǎng)而損傷增大。地震強(qiáng)度不同，結(jié)構(gòu)開(kāi)始產(chǎn)生損傷的時(shí)間不同。在輸入的雙向地震波分別達(dá)到其峰值加速度時(shí)，結(jié)構(gòu)層的損傷指標(biāo)發(fā)生了一個(gè)突然增大的過(guò)程，并且在第二個(gè)地震加速度峰值到達(dá)時(shí)，6倍和7倍El Centro波下結(jié)構(gòu)層的損傷指標(biāo)達(dá)到了0.5，結(jié)構(gòu)很快發(fā)生倒塌傾覆。

圖5 不同地震強(qiáng)度下首層損傷發(fā)展過(guò)Fig.5 Global damage evolution process of the first floot under different earthquake intensity

圖6 為不同地震強(qiáng)度下結(jié)構(gòu)首層的層間位移時(shí)程曲線圖，表4為不同強(qiáng)度地震動(dòng)下結(jié)構(gòu)層的損傷指標(biāo)和最大的層間位移角。地震動(dòng)強(qiáng)度為2倍El Centro波時(shí)，最大的層間位移角為0.013，結(jié)構(gòu)層的損傷指標(biāo)達(dá)到了0.087，對(duì)照本文模型不同的損傷程度對(duì)應(yīng)的損傷指數(shù)范圍表，此時(shí)結(jié)構(gòu)發(fā)生了輕微的破壞;在地震動(dòng)強(qiáng)度為3倍El Centro波時(shí)，圖6(a)中的層間位移曲線圍繞平衡位置上下擺動(dòng)，此時(shí)的層間位移角達(dá)到了0.02左右，損傷指標(biāo)也達(dá)到了0.18，對(duì)照表1，結(jié)構(gòu)發(fā)生了中等破壞;地震動(dòng)強(qiáng)度為5El Centro波時(shí)，層間位移角達(dá)到了0.093，明顯超過(guò)了我國(guó)抗震設(shè)計(jì)規(guī)范規(guī)定的鋼框架結(jié)構(gòu)彈塑性層間位移角1/50的限值，圖6(b)中層間位移也明顯增大，X、Y向的層間位移偏離平衡位置，達(dá)到最大值后圍繞最大值上下波動(dòng)，說(shuō)明變形未達(dá)到其破壞臨界條件，結(jié)構(gòu)未發(fā)生倒塌傾覆，此時(shí)結(jié)構(gòu)的損傷指數(shù)為0.34，對(duì)照表1本文模型不同的損傷程度對(duì)應(yīng)的損傷指數(shù)范圍，已經(jīng)嚴(yán)重破壞;圖6(c)和圖6(d)中，地震動(dòng)強(qiáng)度分別為6倍和7倍的El Centro波，層間位移曲線持續(xù)增大最后趨于發(fā)散，結(jié)構(gòu)層的損傷指標(biāo)大于0.5，說(shuō)明結(jié)構(gòu)已經(jīng)完全倒塌。

表4 不同地震強(qiáng)度下結(jié)構(gòu)損傷指標(biāo)和層間位移角Tab.4 Global damage and inter-story drift angle at different earthquake intensity

4 結(jié)論

圖6 不同地震強(qiáng)度下層間位移時(shí)程曲線Fig.6 Time history responses of relative displacement at different intensity of El-Centro earthquakes

本文從材料的角度，綜合考慮結(jié)構(gòu)在地震中的破壞機(jī)理，提出了一個(gè)基于等效塑性應(yīng)變和比能雙控的雙參數(shù)損傷模型，并在在前人試驗(yàn)和理論分析的基礎(chǔ)上，參考我國(guó)抗震烈度，結(jié)合本文損傷模型的特點(diǎn)，定出了結(jié)構(gòu)不同破壞程度對(duì)應(yīng)的損傷指標(biāo)范圍。結(jié)合一9層benchmark模型結(jié)構(gòu)應(yīng)用該損傷模型對(duì)其在El Centro波作用下的倒塌過(guò)程進(jìn)行了數(shù)值模擬分析，結(jié)果表明隨著地震強(qiáng)度的加大和地震波持續(xù)時(shí)間的增長(zhǎng)結(jié)構(gòu)層的損傷與柱子的損傷增大，地震強(qiáng)度不同，結(jié)構(gòu)開(kāi)始產(chǎn)生損傷的時(shí)間不同。在輸入的雙向地震波分別達(dá)到其峰值加速度時(shí)，結(jié)構(gòu)層的損傷指標(biāo)發(fā)生了一個(gè)突然增大的過(guò)程，并且在第二個(gè)地震加速度峰值到達(dá)時(shí)結(jié)構(gòu)很快發(fā)生倒塌傾覆。驗(yàn)證了該損傷模型評(píng)估強(qiáng)震下鋼結(jié)構(gòu)豎向構(gòu)件及層的損傷發(fā)展過(guò)程的有效性。

［1] 李宏男，肖詩(shī)云，霍林生.汶川地震震害調(diào)查與啟示［J] .建筑結(jié)構(gòu)學(xué)報(bào)，2008，29(4):10 －19.

［2] 張永鑫，凌 煜，張 熠.雙參數(shù)地震損傷模型綜述［J] .山西:山西建筑，2009，35(27):64 －65.

［3] 董 寶，沈祖炎.空間鋼構(gòu)件考慮累積損傷累積效應(yīng)的恢復(fù)力模型及試驗(yàn)驗(yàn)證［J] .上海力學(xué)，1999，20(4):341－347.

［4] 李翌新，趙世春.鋼筋混凝土及勁性鋼筋混凝土構(gòu)件的累積損傷模型［J] .四川:西南交通大學(xué)學(xué)報(bào)，1994，29(4):412－417.

［5] Park Y J，Ang A.Mechanistic seismic damage model for reinforced concrete［J] .Journal of Structural Engineering，ASCE，1985，111(4):740 －756.

［6] 歐進(jìn)萍，牛獲濤，王光遠(yuǎn).多層非線性抗震鋼結(jié)構(gòu)的模糊動(dòng)力可靠性分析與設(shè)計(jì)［J] .地震工程與工程震動(dòng)，1990，10(4):27－37.

［7] 牛狄濤，任利杰.改進(jìn)的鋼筋混凝土結(jié)構(gòu)雙參數(shù)地震破壞模型［J] .地震工程與工程振動(dòng)，1996，16(4):44－54.

［8] 傅劍平，王 敏，白紹良.對(duì)用于鋼筋混凝土結(jié)構(gòu)的Park-Ang雙參數(shù)破壞準(zhǔn)則的識(shí)別和修正［J] .地震工程與工程振動(dòng)，2005，25(5):73 －79.

［9] 韋 未，李同春.一種新的用于各向同性損傷模型的四參數(shù)等效應(yīng)變［J] .工程力學(xué)，2005，22(6):91 －96.

［10] 韋 未，李同春，姚緯明.建立在應(yīng)變空間上的混凝土四參數(shù)破壞準(zhǔn)則［J] .水利水電科技進(jìn)展，2004，24(5):27－29.

［11] Shen Z Y，Dong B.An experiment-based cumulative damage mechanics model of steel under cyclic loading［J] .Advance in Structural Engineering，1997，1(1):39－46.

［12] 宋振森.超高層鋼結(jié)構(gòu)考慮損傷和損傷累積的分析方法［D] .上海:同濟(jì)大學(xué)，2003.

［13] Ju J W.Strain-and stress-based continum damage models-I:formulation［J] . International Journal of Solids and Structures.1987，23(7):821 －840.

［14] Simo JC，Ju，JW.Strain-and stress-based continum damage models-II:computational aspects［J] .International Journal of Solids and Structures.1987，23(7):841 －869.

［15] 陳樹(shù)勛，王光遠(yuǎn).工程結(jié)構(gòu)與系統(tǒng)抗震優(yōu)化設(shè)計(jì)的實(shí)用方法［M] .北京:中國(guó)建筑工業(yè)出版社，1999.

［16] 劉海卿，陳小波，王學(xué)慶.應(yīng)用損傷指數(shù)的框架結(jié)構(gòu)倒塌數(shù)值模擬分析［J] .自然災(zāi)害學(xué)報(bào)，2008，17(5):21 －25.

［17] 姚謙峰，蘇三慶.地震工程［M] .陜西:陜西科學(xué)技術(shù)出版社，2001.

［18] Ohtori Y，Spencer Jr B F，Dyke S J.Benchmark control problems for seismically excited nonlinear buildings［J] .Journal of Engineering Mechanics，ASCE，2004，130(4):366－385.

［19] LS-DYNA.Keyword user's manual［M] .Livermore，California:Livermore Software Technology Corporation，2006.