砌體構件受壓承載力可靠性分析

蔡云梅 趙 成 宋 勇 唐 鵬

(1 塔里木大學水利與建筑工程學院, 新疆 阿拉爾 843300)(2 中國人民解放軍69218部隊， 新疆 喀什 844700)

?

砌體構件受壓承載力可靠性分析

蔡云梅1趙 成2宋 勇1唐 鵬1

(1 塔里木大學水利與建筑工程學院, 新疆 阿拉爾 843300)(2 中國人民解放軍69218部隊， 新疆 喀什 844700)

通過試驗研究與理論分析，建立砌體構件有限元模型，并與試驗結果比較，論證ANSYS中建立模型的合理性，在此基礎上，運用實例采用蒙特卡洛法計算構件失效概率，及相關影響因素的敏感程度。結果表明了運用ANSYS概率分析功能實現結構可靠性分析的可行性。

砌體結構；可靠性；有限元；敏感度

目前在我國建筑結構中砌體結構仍占有很大的比重，其中大部分建筑都已進入到了檢測期，需要對結構采用一定的方法進行檢測評定其可靠度，根據可靠指標評定是否滿足安全等級要求。可靠度的計算需要大量的樣本、計算量大，因而人工計算方法的推廣受到一定的限制，如果采用數值模擬方法能夠滿足樣本數量的要求，提高計算效率，同時具有很好的準確性。

1 試驗概況

根據《砌體基本力學性能試驗方法標準》(GBJ 129-90)和《砌墻磚試驗方法》[1]的要求，對外形尺寸為240 mm×115 mm×53 mm的普通磚，高厚比，當墻厚為240 mm時，試件尺寸為240 mm×370 mm×720 mm。為使所加荷載分布均勻需制作厚度為10 mm的鋼墊板若干塊，鋼墊板尺寸在滿足加載儀器的情況下大于240 mm×370 mm，在鋼墊板四角焊接四個吊鉤，方便試件的移動與起吊。同時試件頂部和底部應采用水泥砂漿找平，做漿后養護三天以上進行承載力試驗，試驗裝置見圖1所示，圖2為磚砌體受壓開裂，裂縫主要發生在豎向灰縫及磚塊所在的位置。

按相關文獻[1]規定測定磚塊及砂漿的的材料參數，試驗結果見表1，將此數據用于ANSYS程序中的建模分析。

圖1 磚砌體抗壓試驗 圖2 磚砌體開裂裂縫分布

表1 相關材料參數

2 有限元建模

有限元分析的軟件有很多種，Ansys軟件是其中之一，由于其強大的功能被廣泛應用于各專業中，在結構分析中的應用主要是數值模擬方面，通過數值模擬能夠代替一些試驗過程，為科學研究工作帶來便利。程序中進行分析需要經過前處理、分析計算及后處理三個模塊[2]。前處理模塊中主要是一些參數的定義、模型的建立與網格劃分；分析計算模塊中需要定義分析類型、定義荷載及載荷步然后進行求解；后處理模塊包括通用后處理和時間歷程后處理，經過后處理用戶可以獲得求解過程的計算結果并對其進行顯示輸出。

關于Ansys中一些砌體結構的求解尚未有一般性的結論[3]。前處理砌體結構的建模主要有兩種方式[4]：整體式和分離式。所謂整體式就是將塊體和砂漿作為整體來考慮，這種模型適用于尺寸較大，且不需要考慮塊體與砂漿之間的粘結摩擦的影響；分離式適合于較小尺寸的單個試件，是通過建立彈簧單元或接觸單元考慮砂漿與塊體間的相互作用，這兩種方法可以根據實際要求加以選擇也可相互結合使用。本文選擇分離式建模，考慮塊體和砂漿特性選擇Solid65單元，主要原因在于該單元能夠模擬類似混凝土的開裂及壓碎。

表2 模型材料參數

表3 荷載分布

塊體與砂漿的本構關系參照相關文獻磚材料本構關系近似選用下式[5]：

其中fc取為塊體的單軸抗壓強度ε0=0.0015

砂漿的本構關系近似用下式表示[6]:

其中為砂漿的抗壓強度，ε0=0.002εcu=0.0033

對于塊體與砂漿間的相互作用采用接觸單元加以連接，本文采用面—面接觸單元，TARGE170模擬目標單元，CONTA174模擬接觸單元，塊體和砂漿破壞準則均選用Willia—Warnker五參數破壞準則[7](表2)，建立的有限元模型如圖3，邊界條件為底面完全固定的形式，頂面作用均布荷載按分步加載，見表3。

經過求解分析結果如下：X向應變分布如圖4，裂縫分布情況如圖5，由應變等值線圖可以看出較大的應變集中在豎向灰縫對應的位置處，豎向的灰縫首先發生破壞，使得灰縫處的塊體處于復雜應力狀態，裂縫進一步發展導致磚塊開裂，圖5為構件裂縫分布圖，荷載施加至30%左右破壞荷載時出現一些細小裂縫，試驗過程中開裂荷載為破壞荷載的40%左右，通過模擬分析結果與試驗對比可以得出模型的破壞與試驗結果基本一致，差別存在的原因在于程序中的模型施加面部荷載是完全接觸，在試驗中為了保證荷載施加均勻分別在試件的上下做漿后加了一塊鋼板，由于試件在砌筑或制作的過程中一些因素的影響使得其與數值模擬仍存在一定差別。

3 可靠性計算

結構的可靠性是指結構在規定的條件下和規定的時間內完成預定功能的能力。在進行結構設計時有兩種常用的方法，定值設計法與概率設計法。定值設計法假定各設計變量為確定的量，依據一定的安全系數進行結構設計；概率法，早在20世紀初就有學者將概率論與數理統計應用于結構安全度的研究[8]，概率設計方法認為，作用于結構上的真實外荷載及結構的真實承載力都是概率意義上的量，設計時不可能精確的確定他們，他們服從一定的分布。通過對各隨機變量統計參數的運算求得功能函數特征值，并利用與失效概率有一一對應關系的可靠指標來衡量可靠度，這就是近似概率法，也是目前計算結構可靠度的主要方法。

可靠性分析時，常以“極限狀態”為標志，當結構或構件達到不滿足它的某一規定使用功能的某一特定狀態時，這個狀態就稱為極限狀態，文中采用極限狀態方程為：

Z=[σ]-σvon

根據以上步驟在進行結構可靠性分析時首先由分析結果中提取相關參量，如：強度、應變等作為控制結構的失效參量；其次，選擇影響結果參量的各不確定因素作為隨機輸入變量，如材料特性、構件尺寸、荷載大小等，通常假定這一類變量服從某種分布，如：正態分布、對數正態分布、均勻分布等，再選擇近似概率分析法，比如蒙特卡洛法、響應面法等。通過分析得到各控制參量的概率分布函數曲線及各影響因素敏感程度。

圖3 磚砌體抗壓有限元模型 圖4 應變分布圖 圖5 裂縫分布圖

算例：以砌體軸壓試件為例，采用第二部分建立的有限元模型，荷載、磚塊及砂漿回彈模量作為隨機輸入變量，根據相關文獻[5，6]可假定各隨機輸入變量服從正態分布(圖6)，數值特征值見表4。

分析方法選擇蒙特卡洛法，拉丁超立方抽樣法抽樣(圖7)，經過200次循環，概率統計分析結果中可以得到指定控制參量的累積概率分布曲線及各輸入變量對各控制參量的敏感程度如圖8，9。通過概率累積分布曲線可以得到滿足極限方程Z<0的概率(失效概率)，推定可靠指標[10]，并評定結構是否處于安全狀態，本例容許應力[σ]=0.619Mpa，則Pf=1.38922×10-3<5×10-3結構較安全， 查得可靠指標β≈3.0基本滿足《建筑結構可靠度設計統一標準》中關于安全等級為二級的建筑物發生延性破壞時的可靠指標β≈3.2。根據敏感性分析結果可以得到各因素對等效應力影響程度。三個因素影響差別不大，塊體回彈模量影響占37. 5%，砂漿占34. 7%，荷載占27. 8%，塊體回彈模量為主要因素，砂漿回彈模量次之，第三是荷載值影響。

4 結論及有待解決的問題

本文通過對砌體結構構件數值模擬分析的結果進行概率統計分析，得到了砌體構件承載力可靠度分布函數以及影響砌體結構構件強度的各因素的敏感程度，在設計過程中針對敏感因素加以控制，使其具有較小的波動性，以保證結構具有足夠的可靠性，對結構可靠度影響不顯著的因素則在進行可靠性分析時可把其視為定值來處理，以減少隨機變量的數目，這樣就能大幅度提高結構可靠性分析的效率。另外，由于數值模擬方法具備的優勢使得其在可靠性分析方面的運用有待進一步加強。

表4 隨機輸入變量特征值

圖6 荷載分布 圖7 應力樣本分布 圖8 應力概率分布曲線 圖9 各因素對應力概率分布影響程度

由于本文只考慮了少數幾個因素對結構的可靠度的影響，控制指標也只選擇了等效應力為參量，具有一定的局限性。實際上除了材料自身的屬性參數外，結構的使用情況、測點數據的采集、分布情況也是有影響的，控制指標根據不同的強度理論應采用不同的參量，分析計算的結果都會發生變化。因而，在今后的研究中適當增加影響因素能使分析結果更加全面，為結構的鑒定加固提出更為全面的措施。

[1]GB/T2542-2003.《砌墻磚試驗方法》[S].

[2] 張洪才，何波.ANSYS13.0從入門到實戰[M]. 北京：機械工業出版社，2011. 494-502.

[3] 李英民，韓軍，劉立平.ANSYS在砌體結構非線性有限元分析中的應用研究[J].重慶建筑大學學報，2006，28(5)：90-105.

[4] 杭翠翠.砌體結構抗壓及抗震承載力仿真分析與試驗對比[D].湖南：湖南大學，2011.21-25.

[5] 劉桂秋.砌體結構基本受力性能的研究[D].湖南：湖南大學，2005.13-51.

[6] 過鎮海，時旭東.剛勁混凝土原理和分析[M].北京：清華大學出版社，2003:15-23.

[7] 羅如登.Ansys中砼單元Solid65的裂縫間剪力傳遞系數取值[J].江蘇大學學報(自然科學版)，2008,29(3)：169-172.

[8] 趙國藩，曹居易，張寬權.工程結構可靠度[M]. 北京：科學出版社，2011: 3-22.

[9] 張偉.結構可靠性理論與應用[M].北京：科學出版社，2008:44-48.

[10] 胡乃君.砌體結構[M].北京：高等教育出版社，2008:15-18.

[11]MingLiu;BaoFengSong;DanielDias;JingPan.TheResearchandApplicationofMasonryStructureReliabilityRatingMethodBasedonExtenics[J].AppliedMechanicsandMaterialsvols,2013(291-294)：1029-1038.

Analyses on the Reliability on Bearing Capacity of Masonry structure

Cai Yunmei1Zhao Cheng2Song Yong1Tang Peng1

(1 College of Water Resources and Architectural Engineering， Tarim University, Alar, Xinjiang 843300)

(2 PLA Unit 69218, Kashi, Xinjiang 844700)

Based on the experimental research and theoretical analysis, masonry structure was simulated by using ANSYS. Compared with the experimental results, the simulation results showed the model is rational. Building on this foundation, selecting example and using the monte carlo method to calculate failure probability and the sensitivity of factors. It showed that the technique of ANSYS probabilistic design analyze structure reliability is feasible.

masonry structure; reliability; finite element; sensitivity

2014-06-05

塔里木大學校長基金(TDZKSS201316)

蔡云梅(1984-)，女，講師，碩士，主要從事結構工程研究。E-mail:314899565@qq.com

1009-0568(2015)01-0111-06

TU

ADOI：10.3969/j.issn.1009-0568.2015.01.020