摘要：本文通過ANSYS 軟件對混雜纖維混凝土試件模型進行單軸壓縮模擬分析，并將分析結(jié)果與相同情況下的試驗所得數(shù)據(jù)進行比較，得出結(jié)論：正確選擇混雜纖維混凝土的本構(gòu)關(guān)系和破壞準則，就能夠較準確的得到纖維混凝土的計算數(shù)據(jù)。
關(guān)鍵詞：混雜纖維；ANSYS；非線性分析
1 引言
混凝土結(jié)構(gòu)是目前使用最為廣泛的一種結(jié)構(gòu)形式，近年來，各種纖維混凝土在國內(nèi)外都得到了較快的實用研究和工程應用，與普通混凝土相比，纖維混凝土在抑制混凝土裂縫發(fā)展，提高抗拉、抗彎強度，特別是提高混凝土的抗沖擊性能和斷裂韌性方面表現(xiàn)出明顯的優(yōu)越性。但是，對混凝土尤其是纖維混凝土的力學性能還不能說已經(jīng)有了全面的掌握[1]。近年來，隨著有限元數(shù)值方法的發(fā)展和計算機技術(shù)的進步，人們已經(jīng)可以利用有限元分析方法對混凝土結(jié)構(gòu)作比較精確的分析了。由于混凝土材料性質(zhì)的復雜性，對纖維混凝土結(jié)構(gòu)進行有限元分析還存在不少困難，其中選擇正確的本構(gòu)關(guān)系和破壞準則、確定符合實際的纖維混凝土應力－應變?nèi)€就是一個重要的方面。
2 單元的選擇及本構(gòu)關(guān)系
2.1 SOLID65 單元
在ANSYS 中通常用SOLID65 單元來模擬混凝土，SOLID65 單元包括一種實體材料和三種加固材料（一般為鋼筋），可以用MAT 命令定義混凝土材料常數(shù)；而加固材料的常數(shù)可以在實常數(shù)中定義，包括材料號、體積率、方向角。
2.2 本構(gòu)關(guān)系和破壞準則
混凝土單元需要定義破壞準則和本構(gòu)關(guān)系，對于纖維混凝土甚至混雜纖維混凝土而言，纖維充當了骨料的角色，所以，選擇正確的本構(gòu)關(guān)系和破壞準則就是進行計算分析的重要因素。ANSYS中的Concrete材料特性用的是Willam-Wamker五參數(shù)破壞準則和拉應力準則的組合模式，能夠較好的反映混凝土從低到高的靜水壓力作用下的破壞特性。但是ANSYS中默認的混凝土的本構(gòu)關(guān)系是線彈性的，即在開裂前的應力－應變關(guān)系為線性，這不符合實際，即使在較低的應力下混凝土也會表現(xiàn)出明顯的非線性。因此，要在材料性質(zhì)中加入反映其本構(gòu)關(guān)系的特性。ANSYS能夠提供非線性彈性本構(gòu)關(guān)系和彈塑性本構(gòu)關(guān)系，單元應力在70%的極限應力以下時，非線性彈性本構(gòu)關(guān)系與試驗結(jié)果符合較好，在非比例加載甚至是反復加載或者是極限荷載作用的情況下，應該用彈塑性本構(gòu)關(guān)系[2]，建議使用等強硬化模型MISO（Multilinear Isotropic Hardening）或者隨動硬化模型MKIN（Multilinear Kinematic Hardening）, 它能反映材料的包辛閣效應。此外，不少書中還推薦使用Drucker-Prager模型（理想彈塑性模型），筆者不建議使用，因為DP模型需要粘聚力、內(nèi)摩擦角、膨脹角這些比較難測出的參數(shù)。
3 關(guān)于收斂問題
ANSYS 混凝土計算收斂是比較困難的，主要影響因素是網(wǎng)格密度、子步數(shù)、收斂準則等。對于網(wǎng)格密度而言，適當?shù)木W(wǎng)格密度能夠收斂。不是網(wǎng)格越密越好，當然太稀疏也不行，這僅僅是就收斂而言的，但是究竟多少合適，只能靠自己針對情況慢慢調(diào)整試算。子步數(shù)的設置很重要，設置太大或太小都不能達到正常收斂。這點可以從收斂過程圖看出，如果F的范數(shù)曲線在收斂曲線上面走形的很長，可考慮增大子步數(shù)。或者根據(jù)經(jīng)驗慢慢調(diào)正試算。實際上收斂精度的調(diào)正并不能徹底解決收斂的問題，但可以放寬收斂條件以加速。一般不超過5%（缺省是0.5%），且使用力收斂條件即可。 另外混凝土壓碎的設置對最后的計算結(jié)果影響也比較大。不考慮壓碎時，計算相對容易收斂；而考慮壓碎則比較難收斂，即便是沒有達到壓碎應力時。如果是正常使用情況下的計算，建議關(guān)掉壓碎選項（即令單軸抗壓強度UnCompSt=-1）；如果想要更準確的極限計算，建議使用真實試驗得出的混凝土本構(gòu)關(guān)系以及多線性各向同性硬化流動律（即CONCR+MISO）且關(guān)閉壓碎檢查。
4 實例分析
為了驗證上述分析的正確性，筆者采用一組150mm×150mm×150mm的鋼/聚丙烯混雜纖維立方體試件作為分析實例。其中，鋼纖維體積率f ρ =1.5%，聚丙烯纖維參量f =0.9kg/m V 3。試驗采用單軸壓縮，試驗所得的混雜纖維混凝土受壓應力－應變關(guān)系如表1所示。
在進行ANSYS 分析時，筆者結(jié)合表1 的真實試驗數(shù)據(jù)，分別采用CONCR+MISO 和CONCR+MKIN 兩種方式（都關(guān)閉壓碎選項）進行計算。建立模型如圖1 所示：
圖2 和圖3 分別是采用兩種本構(gòu)關(guān)系的Y 方向（受壓方向）的位移等值線圖，對比可以看出，采用CONCR+MISO 本構(gòu)關(guān)系更能形象的模擬試塊的真實受壓內(nèi)部變形情況，而采用CONCR+MKIN 本構(gòu)關(guān)系得出的內(nèi)部變形更加線性。而通過圖4 和圖5 的對比我們可知：雖然兩種本構(gòu)關(guān)系都能夠得到同試驗數(shù)據(jù)比較吻合的曲線，但由于CONCR+MKIN 本構(gòu)關(guān)系最多可以設置6 個參數(shù)，而CONCR+MISO 本構(gòu)關(guān)系可設置的參數(shù)更多，所以圖4 的應力－應變曲線更加接近試驗值，各應變點的應力值更加準確。
5 結(jié)論
用ANSYS 來實現(xiàn)纖維混凝土的非線性分析，只要選擇正確的本構(gòu)關(guān)系和破壞準則，在靜水壓力正常的情況下，關(guān)掉壓碎選項，便能夠得到與試驗數(shù)據(jù)較吻合的計算結(jié)果。
參考文獻：
[1] 王寶庭，高丹盈，王占橋，趙軍.混雜纖維高強混凝土斷裂性能試驗研究[J].第十屆全國纖維混凝土學術(shù)會議論文集，上海：同濟大學出版社，2004
[2] 凌廣，吳同樂，賈永剛.鋼筋混凝土有限元分析[M].工程結(jié)構(gòu)，2003
[3] 宋玉普.多種混凝土材料的本構(gòu)關(guān)系和破壞準則，中國水利水電出版社，2002
