RC Z形柱框架節點擬靜力試驗有限元分析

2016-01-19 10:32:00崔欽淑郭顏愷

浙江工業大學學報 2015年6期

崔欽淑，郭顏愷

(浙江工業大學建筑工程學院，浙江杭州 310014)

崔欽淑，郭顏愷

(浙江工業大學建筑工程學院，浙江杭州 310014)

摘要：基于鋼筋混凝土Z形截面柱框架節點低周反復加載試驗結果，運用有限元軟件ANSYS14.5對其進行非線性模擬.針對Z形柱框架節點本身形狀相對規整且內部配筋構造相對復雜的實際情況，根據配筋的實際位置，采用直接生成建模方式建模.利用不同ANSYS的彈塑性材料MKIN模型及MISO模型與不同的混凝土的本構關系Hongnestad模型、Rush模型、Sargin模型進行組合來分別模擬混凝土的材料特性，以研究不同的鋼筋混凝土材料模型對模擬結果的影響.研究表明：MKIN+Sargin的組合模擬結果與試驗結果吻合較好，此外不同的收斂精度控制對模擬結果也有一定的影響.

關鍵詞：鋼筋混凝土Z形柱；框架節點；ANSYS；本構關系；材料模型

Finite element analysis of quasi-static test of Z-shaped column

joints of RC frames

CUI Qinshu, GUO Yankai

(College of Civil Engineering and Architecture, Zhejiang University of Technology, Hangzhou 310014, China)

Abstract:Based on the test results of the Z-shaped column joints of RC frames subjected to low cycle repeated loading, the nonlinear modeling is conducted with ANASYS14.5 software. Since the Z-shaped column joint of frames has a relatively regular shape but comparatively complicated internal reinforcements, the computational model is directly built according to the actual locations of steel bars. The material properties of concrete are modeled using the elasto-plastic MKIN and MISO models with the Hongnestad, Rush, and Sargin constitutive relationships provided by ANASYS to study the effect of different models of RC on the numerical results. It is shown that the numerical results given by the elasto-plastic MKIN model with the Sargin constitutive relationship are in good agreement with the test results. In addition, the control of accuracy and convergence has a certain influence on the numerical results.

Key words:RC Z-shaped column; frame joint; ANSYS; constitutive relationship; material model

隨著計算機技術的迅速發展，利用計算機有限元計算模擬逐漸成為科研的一種有效手段，而且可以有效節約試驗材料用量.ANSYS是目前應用較為廣泛的有限元模擬軟件，在建筑工程領域也有較為廣泛的應用.然而鋼筋混凝土結構，由于混凝土材料的特殊性，使其力學性能較為復雜，給現實有限元模擬研究也帶來了極大的困難性，往往混凝土材料的準確模擬決定了其結果的準確性.國內、外有較多研究者做過相關方面的有限元模擬[1-4]，但并未系統的研究不同的混凝土本構與材料模型組合對試驗結果的影響，而且對于異形柱框架節點結構本身鋼筋配置的復雜性，選擇合理的建模方式與混凝土本構模型，也就決定了其研究結果的可靠性與準確性.

筆者針對不同的ANSYS彈塑性材料模型及不同的混凝土本構關系的組合對混凝土材料特性進行一定的模擬，根據不同模擬結果與試驗結果進行對比，找出與試驗結果最為吻合的模擬模型，并通過骨架曲線，有限元應力云圖等驗證模擬的準確性、可靠性.可為相關鋼筋混凝土異形柱節點試件的ANSYS模擬參數設置提供參考.

1試驗概況

1.1試件設計

取文獻[5]中的節點2(JD2)為研究對象，JD2的肢高肢厚比為3∶1，翼緣高厚比為2∶1.混凝土強度等級為C30(fc=17.19 MPa，Ec=2.85×104MPa)，梁、柱箍筋為HPB235，縱筋為HRB400，混凝土的立方體抗壓強度及鋼筋的材料性能按GB/T 50152—2012《混凝土結構試驗方法標準》[6]實測.根據試驗目的要求，將試驗試件尺寸按照相似理論，將其進行1/2的比例進行縮小.縮小后柱截面尺寸肢高300 mm，肢厚100 mm，翼緣高200 mm，翼緣厚100 mm，梁截面尺寸為100 mm×275 mm.鋼筋的材料試驗結果均值，如表1所示.

表1　鋼筋材料試驗結果均值

1.2加載裝置及加載制度

試驗利用電液伺服加載試驗機進行低周往復加載，來模擬地震作用對Z形柱框架節點試件抗震性能的影響，即采用擬靜力試驗方法.根據(JGJ 101—1996)《建筑抗震試驗方法規程》[7]，采用荷載—位移雙控制的加載制度.首先利用千斤頂在柱底端施加恒定軸壓力，使其達到所需的軸壓比，然后再在梁端加載點利用拉壓作動器同步施加反復豎向荷載，以梁端一根縱向鋼筋達到屈服應變時的梁端位移，作為屈服位移.在試件屈服前的彈性階段，采用荷載控制加載，但當試件屈服后，由于梁端塑性鉸的形成，力與變形之間不再是協調的線性關系，故改用位移控制加載，并且一直加載到荷載下降為峰值荷載的85%，即Pu=0.85Pmax，此時認為試件已破壞，對應的位移即為極限位移Δu.試件的試驗滯回曲線，如圖1所示.

圖1　試驗實測滯回曲線圖Fig.1　Experimental hysteresis curve

2有限元模擬

2.1模型建模方式

對于要模擬的RC Z形柱框架節點試件，ANSYS可以用兩種方法來生成模型：實體建模和直接建模.

對于實體建模，則需要首先確定模型的幾何邊界，利用ANSYS的命令將Z形柱部分，框架節點兩端梁部分分別建出，利用布爾運算進行組合，然后利用切分操作(實際切而不分)切分出不同鋼筋的具體位置，設置不同部分鋼筋的相應材料屬性，并對鋼筋單元的大小及形狀進行控制，進行鋼筋的網格劃分，最后設置混凝土的網格尺寸大小劃分網格.

對于直接建模，根據模型的形狀尺寸，事先編號所需建立模型的諸多節點，并對其進行相應規律性編號，然后將相應的節點進行E(單元生成命令)指令，生成所需的鋼筋單元與混凝土單元.

對于龐大或復雜的模型，尤其是對三維實體模型，實體建模更為合適，相對數據處理會少一些，但在模型生成時會需要大量CPU處理時間，而且在處理網格劃分可能出現失敗.對于中小型、簡單的模型，直接生成較為方便，操作者對幾何形狀及每個單元與節點都有完全的把控.

對于要模擬的Z形柱框架節點，采用直接生成建模方式，因為對于該試件，規模相對較小，但需要分析內部鋼筋等細部受力情況，如果采用實體建模，需要切分出所有縱筋以及箍筋的位置，而且網格劃分也易出錯，與直接建模相比占用大量CPU處理時間，故采用直接生成建模方式，相對準確而不易出錯.

2.2Z形柱節點ANSYS計算模型

根據上述試驗試件尺寸，建立與上述試件尺寸、配筋等完全一致的ANSYS有限元模型.采用直接生成方式建模，混凝土采用鋼筋混凝土工程中最常用的solid 65單元，該單元為3D加筋實體單元，用于模擬無筋或加筋的3D實體單元，具有受拉開裂和受壓破碎的性能，其梁單元尺寸為10 mm×10 mm×80 mm，柱單元尺寸為20 mm×20 mm×40 mm，節點加密區單元尺寸為20 mm×20 mm×10 mm；鋼筋采用link180單元，該三維桿單元在桿軸的方向為拉壓單元，單元的每個節點有三個自由度：沿著節點坐標系X，Y，Z三個方向的平動位移，具有塑性、蠕變、旋轉、大變形和大應變等功能.建模所用數據嚴格按照試驗數據采用，有限元模型如圖2所示.

模型建立后在模型柱上端約束截面的兩個水平方向自由度，即UX，UY自由度，UZ方向施加柱端壓力，使其達到相應試驗所需的軸壓比，柱底端截面施加三個方向的約束，即約束UX，UY，UZ(自由度，將左、右梁加載點截面的所有節點進行UZ方向的耦合，將耦合點作為試驗時的加載點，然后按實際試驗加載時的梁端位移變化在梁兩端進行低周往復加載.為了使模擬結果能較好收斂，加載過程始終采用位移加載控制，且每一荷載步加載位移大小嚴格按照該試件試驗過程中每一荷載步的位移大小進行相應的加載，而且位移控制加載較力控制加載能更好的幫助收斂且模擬出下降段[8].

2.3材料屈服準則與破壞準則

混凝土常用的本構關系[9]有美國學者提出的Hongnestad模型，其表達式為

(1)

(2)

德國Rush建議的模型，下降段用水平段代替：

(3)

水平段σ=fcε<ε0≤εu

(4)

1964年Saenz提出一個關于σ-ε曲線的公式，之后Sargin對該公式進行改進，即

(5)

式中：A=E0/Es，E0為混凝土初始彈性模量，Es=σ0/ε0為應力峰值時的割線模量；k3=σ0/ε0為側限對強度的影響系數，取k3=1時，適合于無側向約束的素混凝土；D為主要影響下降段的參數，D越小，其下降段越為明顯.

混凝土破壞準則采用五參數Willam-Warnke準則，其中張開裂縫傳遞系數取0.5，閉合裂縫傳遞系數取0.95，為得到較準確的計算結果關閉壓碎功能.為使計算更易收斂，選擇solid65單元的KEYPOT選項，令KEYPOT(1)=1，即不考慮形函數的附加項；令KEYPOT(7)=1，即考慮開裂后的拉應力釋放系數，慢慢釋放拉應力以幫助收斂[8].ANSYS中常用的彈塑性材料模型有雙線性隨動強化模型BKIN，多線性隨動強化模型MKIN與KINH，雙線性等向強化模型BISO、多線性等向強化模型MISO.在ANSYS14.5中，MISO模型如果定義混凝土下降段，則會提示錯誤，故將MISO與Rush模型相結合，組合形式如表2所示.

鋼筋采用理想的彈塑性模型(雙線性隨動強化模型BKIN)，如圖3所示.

表2　分析模型表

模擬結果顯示：在用MISO+Rush模擬時試件并未達到荷載極限位移值便已發生破壞，且極限荷載值也小于試驗極限荷載值，而對于其它三種模擬情況，其極限位移與試驗模擬大致相符，其模擬極限荷載都不同程度大于試驗極限荷載，究其原因是因未充分考慮鋼筋混凝土之間的滑移，以及軟件默認混凝土為理想彈塑性材料所致.

圖3　鋼筋的本構模型Fig.3　Constitutive model of steel bar

2.4有限元結果分析

根據試驗與有限元模擬結果繪制相應的骨架曲線，如圖4所示.模擬編號4與試驗骨架曲線對比，如圖5所示.

圖4　骨架曲線對比分析圖Fig.4　Comparative analysis of skeleton curves

圖5　模擬4與試驗骨架曲線對比圖Fig.5　Comparison of simulation and test of skeleton curves

從圖4可看出：模擬編號1，2，4的骨架曲線與試驗骨架曲線的變化趨勢大體一致.由圖5可知模擬編號4與試驗結果吻合較好.由于模擬時完全按照試驗加載的位移進行加載控制，故對于模擬結果的準確性，利用梁端屈服荷載值與梁端極限荷載值兩項進行比較，其中對于屈服荷載值，通過對骨架曲線采用能量等值法[5]求得其屈服點，從而確定其屈服荷載，對比如表3所示.

表3　ANSYS模擬與試驗實測極限荷載值、屈服荷載值對比1)

注：1)Putes，Puexp，Pytes，Pyexp，分別為有限元模擬極限荷載值、試驗極限荷載值、有限元模擬屈服荷載值及試驗屈服荷載值.

針對誤差相對較小的模擬編號4與試驗結果做進一步的對比，繪制其結構變形圖如圖6所示，鋼筋應力圖，如圖7，8所示.

圖6　結構變形圖Fig.6　Deformation of the structure

圖7　縱筋應力云圖Fig.7　Stress cloud of longitudinal reinforcement

圖8　箍筋應力云圖Fig.8　Stress cloud of stirrup

從圖7，8中可以看出：其箍筋與縱筋的應力在加載到最后荷載步時，其梁縱筋在節點附近的強度達到了最大值522.79 MPa，根據實測表1中縱筋抗拉強度約為566.25 MPa，屈服強度約為428.23 MPa，從云圖可看出在梁與翼緣相交處鋼筋達到屈服，而節點附近箍筋模擬其最大強度為438.5 MPa，從箍筋云圖中也可看出節點中翼緣處的箍筋達到屈服強度438.5 MPa，而腹板中的箍筋沒有屈服，與試驗測試結果大致相符.

由于混凝土選用的多線性隨動強化MKIN模型，此模型適用于Mises屈服準則，故繪制von Mises應力.從圖9中可觀察到：節點核心區翼緣混凝土達到了試驗測得的混凝土軸心抗壓強度，而從圖10中可以看出：最終裂縫在節點核心區擴展，試驗中試件的節點核心區翼緣混凝土保護層剝落，節點翼緣箍筋屈服后，承載力快速降低，也可驗證有限元模擬的準確性.

圖9　Von Mises應力Fig.9　Von Mises stress

圖10　節點裂縫圖Fig.10　Crack graph of the joint

3結論

對于鋼筋混凝土Z形柱框架節點擬靜力試驗有限元模擬，利用直接生成建模的方式較為方便、準確，其建模過程也相對簡單，模型生成較為容易.使用多線性隨動強化模型MKIN來模擬Z形柱框架節點的混凝土材料，并且利用Sargin的混凝土本構關系來定義Z形柱框架節點混凝土的應力-應變曲線，其對Z形柱框架節點的模擬結果吻合較好.對于收斂準則的控制，采用位移與力同時控制的收斂準則，可以較為準確模擬出Z形柱框架節點試件真實的試驗過程，研究結果可為混凝土異形柱框架節點結構分析提供參考.

參考文獻：

[1]洪健.基于ANSYS的混凝土框架非線性分析[D].廣州：中山大學，2009.

[2]趙滇生，張蕾春，楊德桔，等.不等肢L形截面型鋼鋼筋混凝土異形柱極限承載力的試驗研究[J].浙江工業大學學報，2010，38(6)：638-642.

[3]OKAMURA H, MAEKAWA K. Nonlinear analysis and constitutive methods of reinforced concrete [M]. Tokyo：University of Tokyo，1991.

[4]蘇友文，馬倩雯，李超飛，等.二次受力下預應力碳纖維布加固混凝土梁抗彎性能研究及數值分析[J].浙江工業大學學報，2014，42(4):378-387.

[5]崔欽淑，楊俊杰，康谷貽.鋼筋混凝土Z形截面柱框架節點抗震性能試驗研究[J].建筑結構學報，2012，33(6)：86-95 .

[6]中華人民共和國住房和城鄉建設部.GB/ T 50152—2012混凝土結構試驗方法標準[S].北京：中國建筑工業出版社，2012.

[7]中國建筑科學研究院.JGJ 101—1996建筑抗震試驗方法規程[S].北京：中國建筑工業出版社，1997.

[8]江見鯨，陸新征. 混凝土結構有限元分析[M]. 北京：清華大學出版社，2005.

[9]王新敏. ANSYS工程結構數值分析[M]. 北京：人民交通出版社，2007.

(責任編輯：陳石平)

文章編號：1006-4303(2015)06-0690-05