基于A(yíng)NSYS的鋼筋混凝土力學(xué)分析

2013-08-05 by:廣州Ansys中心來(lái)源:仿真在線(xiàn)

1 引言

由于鋼筋混凝上材料性質(zhì)復(fù)雜,使其表現(xiàn)出明顯的非線(xiàn)性行為。長(zhǎng)期以來(lái)采用線(xiàn)彈性理論的設(shè)計(jì)方法來(lái)研究鋼筋混凝上結(jié)構(gòu)的應(yīng)力或內(nèi)力,顯然不太合理,盡管有此理論是基于人量試驗(yàn)數(shù)據(jù)上的經(jīng)驗(yàn)公式,還是不能準(zhǔn)確反映混凝上的力學(xué)性能,特別是受力復(fù)雜的重要結(jié)構(gòu),必須采用三維鋼筋混凝上非線(xiàn)性有限元方法才能很好地掌握其力學(xué)性能。利用ANSYS對(duì)鋼筋混凝上結(jié)構(gòu)彈塑性的仿真分析,可以對(duì)結(jié)構(gòu)自開(kāi)始受荷載直到破壞的全過(guò)程進(jìn)行分析,獲得不同階段的受力性能。

2 模型的建立

2.1 單元類(lèi)型的選取

2.1.1 混凝土單元

混凝上選用SOLID65單元。SOLID65單元在普通8節(jié)點(diǎn)三維等參元SOLID45單元的基礎(chǔ)上增加了針對(duì)混凝上材料參數(shù)和整體式鋼筋模型,常被用來(lái)模擬鋼筋混凝上和巖石等抗拉能力遠(yuǎn)大于抗拉能力的非均勻材料,可以模擬混凝土材料的開(kāi)裂和壓碎力學(xué)行為。
2.1.2 鋼筋單元

可采用桿元來(lái)模擬縱筋,一般利用空間單元LINK8單元或空間管單元PIPE20建立鋼筋模型,與混凝上SOLID65單元共用節(jié)點(diǎn)。用COMBINE39來(lái)模擬鋼筋和混凝上之間的粘結(jié)。

2.2 材料本構(gòu)關(guān)系模型

2.2.1 混凝土本構(gòu)模型

彈塑性本構(gòu)關(guān)系把服而和破壞而分開(kāi)來(lái)處理。根據(jù)彈塑性理論建立混凝上的本構(gòu)關(guān)系時(shí),必須對(duì)屈服,條件流動(dòng)法則、硬化法則即塑性模型三要素做出基本假定。ANSYS彈塑性本構(gòu)關(guān)系主要使用Mises屈服準(zhǔn)則或Drucker-Prager屈服準(zhǔn)則。
2.2.2 混凝土破壞準(zhǔn)則

一般強(qiáng)度準(zhǔn)則的參數(shù)越多,對(duì)混凝土強(qiáng)度性能的描述就越準(zhǔn)確,多參數(shù)模型大多基于強(qiáng)度試驗(yàn)的統(tǒng)計(jì)而進(jìn)行曲線(xiàn)擬合。ANSYS中的混凝上材料特性用改進(jìn)的William Wamke五參數(shù)破壞破壞準(zhǔn)則和拉應(yīng)力的組合模式,其破壞而子午線(xiàn)和偏平而見(jiàn)相關(guān)文獻(xiàn)。

2.2.3 鋼筋本構(gòu)模型

一般采用雙線(xiàn)形理想彈塑性模型,在A(yíng)NSYS中,鋼筋可以選擇經(jīng)典雙線(xiàn)性隨動(dòng)強(qiáng)化模型(BKIN)和雙線(xiàn)性等向強(qiáng)化模型(BISO)。
2.2.4 本構(gòu)關(guān)系和破壞準(zhǔn)則在A(yíng)NSYS中的實(shí)現(xiàn)

在SOLID65單元中實(shí)現(xiàn)上式的本構(gòu)關(guān)系需通過(guò)定義數(shù)據(jù)表。對(duì)于鋼筋材料,定義一個(gè)應(yīng)力應(yīng)變關(guān)系的數(shù)據(jù)表;而對(duì)于混凝土模型,則要定義兩個(gè)數(shù)據(jù)表,一個(gè)是本構(gòu)關(guān)系的數(shù)據(jù)表,另一個(gè)用于定義混凝上的破壞準(zhǔn)則。通過(guò)混凝上單元材料特性表定義混凝上的裂縫張開(kāi)剪力傳遞系數(shù)閉合剪力傳遞系數(shù)、單軸和多軸抗壓強(qiáng)度等9個(gè)參數(shù),定義混凝上的W illian Wamke強(qiáng)度準(zhǔn)則。

2.3 鋼筋混凝土結(jié)構(gòu)有限元模型的建立

鋼筋混凝上結(jié)構(gòu)的有限元模型可采用分離式、組合式和整體式,各種模式的優(yōu)缺點(diǎn)見(jiàn)文獻(xiàn)。其中整體式是將鋼筋分布于整個(gè)單元中,并把單元視為連續(xù)均勻材料,這類(lèi)模型假定聯(lián)結(jié)較好,不考慮鋼筋和混凝上之間的滑移,通過(guò)鋼筋和混凝上之間公用節(jié)點(diǎn)來(lái)實(shí)現(xiàn)位移協(xié)調(diào)。

2.4 邊界條件處理

2.4.1 在支座和集中力作用處加剛性墊塊

因?yàn)樗休d荷都是通過(guò)有限人小的壓力面來(lái)添加或傳遞到真實(shí)部件上去的,所以可在支座和集中力處加以處理使其接近真實(shí)的受力狀態(tài),從而消除計(jì)算中的應(yīng)力奇異現(xiàn)象。一般加一個(gè)小尺寸的墊塊,可用SOLID45單元來(lái)模擬定義剛性墊塊,通過(guò)公用節(jié)點(diǎn)使剛性墊節(jié)點(diǎn)與相應(yīng)的混凝上節(jié)點(diǎn)位移協(xié)調(diào)。經(jīng)過(guò)ANSYS計(jì)算,此方法能有效解決應(yīng)力集中現(xiàn)象。

2.4.2 將集中力等效為均布力作用在微元上

選擇一些與作用的中心點(diǎn)相鄰近的小單元來(lái)共同承受邊界條件,將集中力等效成均布力作用在選定的單元上或鄰近的節(jié)點(diǎn)上。這一方法更接近真實(shí)的受力情況。

2.5 網(wǎng)格的劃分和求解

模型的分析是非線(xiàn)性靜力分析過(guò)程,非線(xiàn)性方程組的求解采用New-Raphson迭代法。ANSYS為非線(xiàn)性分析提供了多種增強(qiáng)收斂的辦法,影響計(jì)算結(jié)果收斂卞要因素是荷載水平、網(wǎng)格密度、子步數(shù)、收斂標(biāo)準(zhǔn)等,可通過(guò)一此措施來(lái)改進(jìn)和提高收斂速度與解的精度。

2.5.1 控制網(wǎng)格和單元尺寸

劃分密度適當(dāng)?shù)木W(wǎng)格有助于收斂;在實(shí)際應(yīng)用過(guò)程中,應(yīng)該對(duì)單元?jiǎng)澐诌M(jìn)行有效控制,當(dāng)最小單元尺寸大于5cm時(shí),就可以有效避免應(yīng)力集中帶來(lái)的問(wèn)題。

2.5.2 收斂標(biāo)準(zhǔn)控制

一般位移控制加載最好用位移的co-范數(shù)控制收斂,而用力控制加載時(shí)可以用殘余力的2-范數(shù)控制收斂。在裂紋剛剛出現(xiàn)和接近破壞的階段,可以適當(dāng)放松收斂標(biāo)準(zhǔn),保證計(jì)算的連續(xù)性,誤差控制可以在2%-3%之間,一般不超過(guò)5%。

2.5.3 荷載步和了步的設(shè)置

NSUBST設(shè)置的太大或太小都不能達(dá)到正常收斂,必須合理設(shè)置子步數(shù)并且打開(kāi)線(xiàn)性搜索、預(yù)測(cè)等項(xiàng)以加速收斂。荷載步增量對(duì)收斂有較人的影響,必須合理設(shè)置荷載步以保證結(jié)果收斂穩(wěn)定性和有效逼近問(wèn)題的真實(shí)解。

3 算例

本算例為一鋼筋混凝土簡(jiǎn)支梁的開(kāi)裂問(wèn)題?；炷捎肧OLID65模擬,圖1給出了混凝土簡(jiǎn)支梁的有限元模型。圖2給出了混凝土的參數(shù)的輸入界面。