RC Z形柱框架節(jié)點(diǎn)擬靜力試驗(yàn)有限元分析

崔欽淑，郭顏愷

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

RC Z形柱框架節(jié)點(diǎn)擬靜力試驗(yàn)有限元分析

崔欽淑，郭顏愷

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

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

關(guān)鍵詞：鋼筋混凝土Z形柱；框架節(jié)點(diǎn)；ANSYS；本構(gòu)關(guān)系；材料模型

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

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

筆者針對(duì)不同的ANSYS彈塑性材料模型及不同的混凝土本構(gòu)關(guān)系的組合對(duì)混凝土材料特性進(jìn)行一定的模擬，根據(jù)不同模擬結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比，找出與試驗(yàn)結(jié)果最為吻合的模擬模型，并通過(guò)骨架曲線，有限元應(yīng)力云圖等驗(yàn)證模擬的準(zhǔn)確性、可靠性.可為相關(guān)鋼筋混凝土異形柱節(jié)點(diǎn)試件的ANSYS模擬參數(shù)設(shè)置提供參考.

1試驗(yàn)概況

1.1試件設(shè)計(jì)

取文獻(xiàn)[5]中的節(jié)點(diǎn)2(JD2)為研究對(duì)象，JD2的肢高肢厚比為3∶1，翼緣高厚比為2∶1.混凝土強(qiáng)度等級(jí)為C30(fc=17.19 MPa，Ec=2.85×104MPa)，梁、柱箍筋為HPB235，縱筋為HRB400，混凝土的立方體抗壓強(qiáng)度及鋼筋的材料性能按GB/T 50152—2012《混凝土結(jié)構(gòu)試驗(yàn)方法標(biāo)準(zhǔn)》[6]實(shí)測(cè).根據(jù)試驗(yàn)?zāi)康囊螅瑢⒃囼?yàn)試件尺寸按照相似理論，將其進(jìn)行1/2的比例進(jìn)行縮小.縮小后柱截面尺寸肢高300 mm，肢厚100 mm，翼緣高200 mm，翼緣厚100 mm，梁截面尺寸為100 mm×275 mm.鋼筋的材料試驗(yàn)結(jié)果均值，如表1所示.

表1　鋼筋材料試驗(yàn)結(jié)果均值

1.2加載裝置及加載制度

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

圖1　試驗(yàn)實(shí)測(cè)滯回曲線圖Fig.1　Experimental hysteresis curve

2有限元模擬

2.1模型建模方式

對(duì)于要模擬的RC Z形柱框架節(jié)點(diǎn)試件，ANSYS可以用兩種方法來(lái)生成模型：實(shí)體建模和直接建模.

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

對(duì)于直接建模，根據(jù)模型的形狀尺寸，事先編號(hào)所需建立模型的諸多節(jié)點(diǎn)，并對(duì)其進(jìn)行相應(yīng)規(guī)律性編號(hào)，然后將相應(yīng)的節(jié)點(diǎn)進(jìn)行E(單元生成命令)指令，生成所需的鋼筋單元與混凝土單元.

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

對(duì)于要模擬的Z形柱框架節(jié)點(diǎn)，采用直接生成建模方式，因?yàn)閷?duì)于該試件，規(guī)模相對(duì)較小，但需要分析內(nèi)部鋼筋等細(xì)部受力情況，如果采用實(shí)體建模，需要切分出所有縱筋以及箍筋的位置，而且網(wǎng)格劃分也易出錯(cuò)，與直接建模相比占用大量CPU處理時(shí)間，故采用直接生成建模方式，相對(duì)準(zhǔn)確而不易出錯(cuò).

2.2Z形柱節(jié)點(diǎn)ANSYS計(jì)算模型

根據(jù)上述試驗(yàn)試件尺寸，建立與上述試件尺寸、配筋等完全一致的ANSYS有限元模型.采用直接生成方式建模，混凝土采用鋼筋混凝土工程中最常用的solid 65單元，該單元為3D加筋實(shí)體單元，用于模擬無(wú)筋或加筋的3D實(shí)體單元，具有受拉開(kāi)裂和受壓破碎的性能，其梁?jiǎn)卧叽鐬?0 mm×10 mm×80 mm，柱單元尺寸為20 mm×20 mm×40 mm，節(jié)點(diǎn)加密區(qū)單元尺寸為20 mm×20 mm×10 mm；鋼筋采用link180單元，該三維桿單元在桿軸的方向?yàn)槔瓑簡(jiǎn)卧?，單元的每個(gè)節(jié)點(diǎn)有三個(gè)自由度：沿著節(jié)點(diǎn)坐標(biāo)系X，Y，Z三個(gè)方向的平動(dòng)位移，具有塑性、蠕變、旋轉(zhuǎn)、大變形和大應(yīng)變等功能.建模所用數(shù)據(jù)嚴(yán)格按照試驗(yàn)數(shù)據(jù)采用，有限元模型如圖2所示.

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

2.3材料屈服準(zhǔn)則與破壞準(zhǔn)則

混凝土常用的本構(gòu)關(guān)系[9]有美國(guó)學(xué)者提出的Hongnestad模型，其表達(dá)式為

(1)

(2)

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

(3)

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

(4)

1964年Saenz提出一個(gè)關(guān)于σ-ε曲線的公式，之后Sargin對(duì)該公式進(jìn)行改進(jìn)，即

(5)

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

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

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

表2　分析模型表

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

圖3　鋼筋的本構(gòu)模型Fig.3　Constitutive model of steel bar

2.4有限元結(jié)果分析

根據(jù)試驗(yàn)與有限元模擬結(jié)果繪制相應(yīng)的骨架曲線，如圖4所示.模擬編號(hào)4與試驗(yàn)骨架曲線對(duì)比，如圖5所示.

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

圖5　模擬4與試驗(yàn)骨架曲線對(duì)比圖Fig.5　Comparison of simulation and test of skeleton curves

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

表3　ANSYS模擬與試驗(yàn)實(shí)測(cè)極限荷載值、屈服荷載值對(duì)比1)

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

針對(duì)誤差相對(duì)較小的模擬編號(hào)4與試驗(yàn)結(jié)果做進(jìn)一步的對(duì)比，繪制其結(jié)構(gòu)變形圖如圖6所示，鋼筋應(yīng)力圖，如圖7，8所示.

圖6　結(jié)構(gòu)變形圖Fig.6　Deformation of the structure

圖7　縱筋應(yīng)力云圖Fig.7　Stress cloud of longitudinal reinforcement

圖8　箍筋應(yīng)力云圖Fig.8　Stress cloud of stirrup

從圖7，8中可以看出：其箍筋與縱筋的應(yīng)力在加載到最后荷載步時(shí)，其梁縱筋在節(jié)點(diǎn)附近的強(qiáng)度達(dá)到了最大值522.79 MPa，根據(jù)實(shí)測(cè)表1中縱筋抗拉強(qiáng)度約為566.25 MPa，屈服強(qiáng)度約為428.23 MPa，從云圖可看出在梁與翼緣相交處鋼筋達(dá)到屈服，而節(jié)點(diǎn)附近箍筋模擬其最大強(qiáng)度為438.5 MPa，從箍筋云圖中也可看出節(jié)點(diǎn)中翼緣處的箍筋達(dá)到屈服強(qiáng)度438.5 MPa，而腹板中的箍筋沒(méi)有屈服，與試驗(yàn)測(cè)試結(jié)果大致相符.

由于混凝土選用的多線性隨動(dòng)強(qiáng)化MKIN模型，此模型適用于Mises屈服準(zhǔn)則，故繪制von Mises應(yīng)力.從圖9中可觀察到：節(jié)點(diǎn)核心區(qū)翼緣混凝土達(dá)到了試驗(yàn)測(cè)得的混凝土軸心抗壓強(qiáng)度，而從圖10中可以看出：最終裂縫在節(jié)點(diǎn)核心區(qū)擴(kuò)展，試驗(yàn)中試件的節(jié)點(diǎn)核心區(qū)翼緣混凝土保護(hù)層剝落，節(jié)點(diǎn)翼緣箍筋屈服后，承載力快速降低，也可驗(yàn)證有限元模擬的準(zhǔn)確性.

圖9　Von Mises應(yīng)力Fig.9　Von Mises stress

圖10　節(jié)點(diǎn)裂縫圖Fig.10　Crack graph of the joint

3結(jié)論

對(duì)于鋼筋混凝土Z形柱框架節(jié)點(diǎn)擬靜力試驗(yàn)有限元模擬，利用直接生成建模的方式較為方便、準(zhǔn)確，其建模過(guò)程也相對(duì)簡(jiǎn)單，模型生成較為容易.使用多線性隨動(dòng)強(qiáng)化模型MKIN來(lái)模擬Z形柱框架節(jié)點(diǎn)的混凝土材料，并且利用Sargin的混凝土本構(gòu)關(guān)系來(lái)定義Z形柱框架節(jié)點(diǎn)混凝土的應(yīng)力-應(yīng)變曲線，其對(duì)Z形柱框架節(jié)點(diǎn)的模擬結(jié)果吻合較好.對(duì)于收斂準(zhǔn)則的控制，采用位移與力同時(shí)控制的收斂準(zhǔn)則，可以較為準(zhǔn)確模擬出Z形柱框架節(jié)點(diǎn)試件真實(shí)的試驗(yàn)過(guò)程，研究結(jié)果可為混凝土異形柱框架節(jié)點(diǎn)結(jié)構(gòu)分析提供參考.

參考文獻(xiàn)：

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

[2]趙滇生，張蕾春，楊德桔，等.不等肢L形截面型鋼鋼筋混凝土異形柱極限承載力的試驗(yàn)研究[J].浙江工業(yè)大學(xué)學(xué)報(bào)，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]蘇友文，馬倩雯，李超飛，等.二次受力下預(yù)應(yīng)力碳纖維布加固混凝土梁抗彎性能研究及數(shù)值分析[J].浙江工業(yè)大學(xué)學(xué)報(bào)，2014，42(4):378-387.

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

[6]中華人民共和國(guó)住房和城鄉(xiāng)建設(shè)部.GB/ T 50152—2012混凝土結(jié)構(gòu)試驗(yàn)方法標(biāo)準(zhǔn)[S].北京：中國(guó)建筑工業(yè)出版社，2012.

[7]中國(guó)建筑科學(xué)研究院.JGJ 101—1996建筑抗震試驗(yàn)方法規(guī)程[S].北京：中國(guó)建筑工業(yè)出版社，1997.

[8]江見(jiàn)鯨，陸新征. 混凝土結(jié)構(gòu)有限元分析[M]. 北京：清華大學(xué)出版社，2005.

[9]王新敏. ANSYS工程結(jié)構(gòu)數(shù)值分析[M]. 北京：人民交通出版社，2007.

(責(zé)任編輯：陳石平)

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

中圖分類(lèi)號(hào)：TU375.3

文獻(xiàn)標(biāo)志碼：A

作者簡(jiǎn)介：崔欽淑(1963—)，女，山東萊西人，副教授，國(guó)家一級(jí)注冊(cè)結(jié)構(gòu)工程師，研究方向?yàn)榛炷廉愋沃Y(jié)構(gòu)抗震性能，E-mail:cuiqinshu@zjut.edu.cn．

基金項(xiàng)目：浙江省自然科學(xué)基金資助項(xiàng)目(LY14E080007)

收稿日期：2015-05-07