方自虎，李向鵬，簡旭陽，洪博愷

(深圳大學(xué) 土木工程學(xué)院，廣東 深圳518060)

鋼筋混凝土梁柱節(jié)點超自由度單元

方自虎，李向鵬，簡旭陽，洪博愷

(深圳大學(xué) 土木工程學(xué)院，廣東 深圳518060)

為分析RC框架結(jié)構(gòu)在循環(huán)荷載作用下的非線性滯回性能，提出一個新的鋼筋混凝土梁柱節(jié)點單元——“超自由度單元”.單元在梁與節(jié)點交界面和柱與節(jié)點交界面被劃分成“節(jié)點截面”和“梁柱截面”，節(jié)點的力學(xué)性能由4節(jié)點單元描述，而梁柱受力鋼筋與節(jié)點核心區(qū)的粘結(jié)滑移由存在于“節(jié)點截面”和“梁柱截面”之間的8根彈簧控制.單元具有4個外節(jié)點和4個內(nèi)節(jié)點，每個內(nèi)節(jié)點具有2個自由度，每個外節(jié)點具有3個自由度，該3個自由度與普通梁單元一致，從而確保本單元適合于同普通一維梁柱單元一起進行鋼筋混凝土結(jié)構(gòu)平面非線性分析.將內(nèi)節(jié)點上的自由度依附到外節(jié)點上，單元在數(shù)值表現(xiàn)上具有4個節(jié)點20個自由度.通過試驗和模擬分析結(jié)果對比，結(jié)果表明：結(jié)構(gòu)極限承載能力和滯回曲線的捏攏特性都能被單元很好模擬，同時計算機的運行時間非常少.本模型適合于進行循環(huán)荷載作用下平面框架結(jié)構(gòu)非線性響應(yīng)分析.

梁柱節(jié)點單元；鋼筋混凝土；框架結(jié)構(gòu)；循環(huán)荷載；非線性分析;滯回性能

鋼筋混凝土結(jié)構(gòu)的梁柱節(jié)點單元[1]，根據(jù)節(jié)點變形和受力機理，將節(jié)點從功能上區(qū)分為節(jié)點核心區(qū)剪切塊的剪切變形、梁柱與核心區(qū)交界面的剪切變形和梁柱受力鋼筋與核心區(qū)混凝土的粘結(jié)滑移三大功能組成部分，并分別用三組共13根一維彈簧來描述該三大組成部分的材料模型.由于該單元的組成方式不同于一般單元，因此，該節(jié)點單元被稱為“超級單元”.其中，核心區(qū)剪切塊假定受純剪切作用，因此用一根彈簧描述剪力與剪切變形之間的關(guān)系；梁柱與核心區(qū)4個交界面剪力與剪切位移關(guān)系用4根彈簧描述；而位于節(jié)點兩邊的梁上下部鋼筋以及位于節(jié)點上下的柱左右兩邊鋼筋共用8根彈簧描述鋼筋與核心區(qū)混凝土之間的粘結(jié)滑移關(guān)系.

對于超級單元核心區(qū)剪切塊材料模型怎樣建立，不同的學(xué)者有不同的主張.超級單元的建立者[1-2]根據(jù)MCFT理論[3]，利用有限元方法計算核心區(qū)的剪力和剪切變形之間的關(guān)系.LaFave等[4]則依據(jù)試驗統(tǒng)計資料，根據(jù)統(tǒng)計分析模型建立的經(jīng)驗公式確定該材料的物理關(guān)系.從而可知，超級單元中核心區(qū)的剪切材料模型怎樣建立，對于超級單元的使用者來說是一個難點.

為了避免核心區(qū)剪切材料物理關(guān)系建立的困難，本文將超級單元的核心區(qū)用平面四節(jié)點單元代替，同時考慮梁柱鋼筋與核心區(qū)混凝土的粘結(jié)滑移作用，建立一種新的基于平面4節(jié)點單元的梁柱節(jié)點單元，并在商業(yè)平臺ABAQUS上開發(fā)了應(yīng)用程序.新單元的核心區(qū)材料模型由鋼筋和混凝土材料根據(jù)平面應(yīng)力狀態(tài)的MCFT理論[5]建立，不需要使用者提供額外的剪切與剪切變形關(guān)系，避免了使用超級單元時不易提供剪切物理關(guān)系的困難.本單元與超級單元一樣，分析結(jié)構(gòu)時，梁柱由一般梁柱單元模擬，因此，本單元的節(jié)點自由度與一般梁柱單元的節(jié)點自由度能夠順利銜接.通過循環(huán)荷載下鋼筋混凝土內(nèi)、邊和角框架3個節(jié)點的試驗結(jié)果數(shù)值分析，驗證了本單元的有效性.

1 節(jié)點單元構(gòu)成

本單元由4個外節(jié)點和4個內(nèi)節(jié)點組成，其中外節(jié)點的自由度為3個，與普通梁單元自由度一致，見圖1；而內(nèi)節(jié)點自由度與普通平面應(yīng)力單元的自由度一致，內(nèi)節(jié)點自由度見圖2.圖1中，外部的數(shù)字代表外節(jié)點編號，內(nèi)部的數(shù)字則代表內(nèi)節(jié)點編號.

圖1 單元的節(jié)點和自由度

考慮鋼筋粘結(jié)滑移和梁柱與節(jié)點交界面裂縫影響，本節(jié)點單元將節(jié)點與梁柱交界面處截面，分成“節(jié)點截面”和“梁柱截面”(圖1)，節(jié)點中鋼筋與混凝土之間的粘結(jié)滑移關(guān)系由梁柱截面和節(jié)點截面上的自由度一起反映.

本節(jié)點的構(gòu)成思想是：針對文獻[1]超級單元中反映核心區(qū)剪切，改由平面四節(jié)點單元代替，而反映梁柱受力鋼筋與混凝土粘結(jié)滑移的8根彈簧，以及節(jié)點與核心區(qū)交界面的4個剪切彈簧，本單元保留，即本單元由一個四節(jié)點平面單元組合8根粘結(jié)滑移彈簧和4根剪切彈簧而成.

由于內(nèi)外節(jié)點的自由度不同，不利于在通用商業(yè)平臺上進行二次開發(fā)，因此，本文將內(nèi)節(jié)點的自由度“依附”到外節(jié)點上，形成一種“超自由度”形式，故不妨稱本梁柱節(jié)點單元為“超自由度”單元.經(jīng)過上述自由度依附之后，在數(shù)值表現(xiàn)上，本文的梁柱節(jié)點單元為：“單元具有4個節(jié)點，每個節(jié)點具有5個自由度”的單元幾何構(gòu)造.程序?qū)崿F(xiàn)時，以單元的外節(jié)點作為單元節(jié)點，內(nèi)節(jié)點按照圖1形式由程序內(nèi)定形成.

超自由度單元的節(jié)點位移向量為

(1)

式中Di為i(i=1～4)節(jié)點的超級自由度向量，以節(jié)點1為例：

(2)

式中前3個自由度Δ1、Δ2和Δ3屬于外節(jié)點，后2個自由度u1和v1屬于內(nèi)節(jié)點，即內(nèi)節(jié)點1的水平和豎向位移(圖2)，式(1)中其余3個節(jié)點的自由度向量依次類推.

超級自由度的實現(xiàn)方式，以通用商業(yè)平臺ABAQUS為例，單元節(jié)點前13個自由度具有固定意義，如1至6自由度分別為X、Y和Z坐標方向的位移和轉(zhuǎn)角.式(2)中前3個自由度正好對應(yīng)1、2和6，而式(2)中后2個自由度可以采用14和15，對應(yīng)的輸入命令如下：

*USER ELEMENT, NODES=4, TYPE=U4, ……

1,2,6,14,15

內(nèi)節(jié)點單元為普通的平面應(yīng)力單元，見圖2.hb為節(jié)點的高度，也即梁高；hc為節(jié)點的寬度，即柱寬.下面分別從變形關(guān)系、平衡關(guān)系和物理關(guān)系來介紹本單元，并在該三大關(guān)系的基礎(chǔ)上導(dǎo)出單元剛度矩陣和節(jié)點抗力向量.

圖2 節(jié)點核心區(qū)單元

2 梁柱截面與節(jié)點截面變形關(guān)系

節(jié)點截面與梁柱截面之間的變形關(guān)系見圖3，其中ri(i=1～4)為剪切變形，Si(i=1～8)為鋼筋與混凝土之間的粘結(jié)滑移變形，該變形也就是模擬粘結(jié)滑移彈簧的伸長量.

顯然，剪切變形與節(jié)點自由度之間有如下關(guān)系：

(3)

假定fb和fc分別為梁上下鋼筋間距和柱子兩側(cè)鋼筋間距與截面高度的比值，即fbhb(hb見圖2)為梁上下鋼筋之間的間距[1].

圖3 兩截面之間的位移關(guān)系

令：αc=0.5hcfc、βc=(1+fc)/2、γc=(1-fc)/2，相應(yīng)的αb，βb和γb僅是將式中下標c改為b即可，則彈簧伸長量與節(jié)點自由度有如下關(guān)系：

(4)

3 平衡關(guān)系

單元的節(jié)點抗力向量為

F={f1f2f3f4}T.

(5)

式中:f1={F1F2F3X1Y1}T為對應(yīng)節(jié)點1超級自由度向量D1的節(jié)點抗力向量，其余3個節(jié)點的抗力向量依次類推.Fi與圖3節(jié)點位移Δi相對應(yīng)；Xi和Yi與圖2中節(jié)點位移ui和vi相對應(yīng)；8根粘結(jié)滑移彈簧力用Pi表示，Pi與圖3彈簧伸長Si相對應(yīng)；模擬截面剪切變形的4根彈簧力用Qi表示，Qi與圖3剪切位移ri相對應(yīng).

粘結(jié)滑移彈簧力與單元節(jié)點抗力之間的平衡關(guān)系為

(6)

截面剪切彈簧力與單元節(jié)點抗力之間的平衡關(guān)系為

(7)

4 物理關(guān)系

本節(jié)點存在4種材料的物理關(guān)系：(1)平面單元的混凝土材料；(2)平面單元內(nèi)鋼筋材料，包括節(jié)點內(nèi)的箍筋和梁柱穿插于節(jié)點內(nèi)的受力鋼筋；(3)梁柱受力鋼筋與節(jié)點核心區(qū)內(nèi)混凝土之間的粘結(jié)滑移；(4)梁柱截面與節(jié)點截面之間的剪切變形物理關(guān)系.其中，(1)和(2)中混凝土材料采用MCFT理論建立物理關(guān)系[3,5]，鋼筋材料采用根據(jù)Ramberg-Osgood公式建立的滯回鋼筋模型[3，5]，兩種材料在商業(yè)平臺ABAQUS上的UMAT接口程序已經(jīng)過數(shù)值檢驗[5-6].而(3)中的粘結(jié)滑移物理關(guān)系，由于超級單元使用的粘結(jié)滑移關(guān)系[1-2]已經(jīng)過OpenSees大量用戶的檢驗，因此，在本節(jié)點上直接采用該材料模型.對于上述(4)變形關(guān)系[7]，由于本文的重點不在于建立這樣的剪切模型，因此，在本文的數(shù)值實現(xiàn)過程中假定該關(guān)系為線彈性的.

5 剛度矩陣

圖2所示的平面單元剛度矩陣為

(8)

式中：N是高斯積分點數(shù)，W(ξj)和W(ηm)為高斯積分權(quán)重，ξj和ηm為高斯點的自然坐標值，V是節(jié)點的體積，D是根據(jù)MCFT理論得到的混凝土材料應(yīng)力應(yīng)變雅可比矩陣[5]和根據(jù)滯回鋼筋模型得到的節(jié)點箍筋和梁柱受力鋼筋材料應(yīng)力應(yīng)變雅可比矩陣[5]，B是平面單元的應(yīng)變矩陣，如式(9)所示.

(9)

式中

(10)

其余3個類推可得.式(8)建立的剛度矩陣為對應(yīng)單元的內(nèi)節(jié)點自由度部分，而單元的總剛度矩陣還要考慮8根粘結(jié)滑移彈簧和4根剪切彈簧引起的剛度.

粘結(jié)滑移彈簧剛度矩陣，以彈簧1為例，余者依次類推.將式(4)中彈簧1的位移S1寫成矩陣形式，有：

S1=Bs1Δs1,

(11)

式中：Bs1=[1αb-βb-γb]，Δs1是{Δ1Δ3u1u2}T.

此時彈簧的剛度記為Ds1，則其剛度矩陣為

(12)

同理，將式(3)的剪切彈簧變形寫成矩陣形式，以r1為例，則有：

r1=Br1Δr1,

(13)

式中：Br1=[1-0.5-0.5]，Δr1是{Δ2v1v2}T. 若剪切剛度記為Dr1，則其剛度矩陣為

(14)

將式(8)、(12)和(14)組合，即可得到節(jié)點單元的總剛度矩陣.

6 節(jié)點抗力向量

平面單元(圖2)在當前應(yīng)力狀態(tài)下的節(jié)點抗力向量為

(15)

式中：X={N1N2N3N4}T，Ni即為{XiYi}T，σ為平面單元高斯點的應(yīng)力向量.8根彈簧引起的節(jié)點抗力，仍然以彈簧1為例，由式(4)知，若彈簧1當前狀態(tài)的彈簧力為P1，則其引起的節(jié)點抗力向量為

(16)

式中：Fs1是{F1F3X1X2}T，Bs1同式(11).4根剪切變形彈簧引起的節(jié)點抗力向量，以r1為例則有：

(17)

式中：Fr1是{F2Y1Y2}T，而Q1為對應(yīng)的剪力.組合式(15)、(16)和(17)即可得到總節(jié)點抗力向量.

7 力與位移關(guān)系

非線性問題的有限元求解過程中，對于每個荷載增量步每次循環(huán)，都需要求解下式：

K·δ=P-F.

(18)

式中：K即為由式(8)、(12)和(14)求得的當前循環(huán)步的單元剛度矩陣，δ為位移增量向量，P為當前增量步的總荷載向量，F(xiàn)為由式(15)、(16)和(17)求得的節(jié)點抗力向量.

8 驗 證

本文選擇內(nèi)框架節(jié)點UNIT1[8]、邊框架節(jié)點UNIT3[9]和角框架節(jié)點KJ5[10]的試驗結(jié)果來驗證單元的有效性.下列算例中，梁柱單元采用ABAQUS的B21梁單元，梁柱鋼筋用rebar形式，即不考慮梁柱的鋼筋與混凝土間粘結(jié)滑移[5-6].梁柱單元采用分層法，混凝土材料采用OpenSees中的concrete01模型，即忽略混凝土拉應(yīng)力.鋼筋的滯回模型采用Vecchio模型[5].

8.1 算例1-內(nèi)框架節(jié)點UNIT1

本算例[8]是超級單元在OpenSees平臺上驗證的算例，相關(guān)試驗數(shù)據(jù)也可于OpenSees網(wǎng)站獲得.分析中節(jié)點每邊的梁柱都劃分10個單元，其中梁上荷載作用點到節(jié)點之間分4個單元(ABAQUS計算模型見圖4).

圖4 UNIT1分析模型

圖5給出UNIT1柱頂水平荷載和水平位移滯回曲線的試驗[8]與本文分析結(jié)果對比.可以看出，本文提出的單元，在節(jié)點承載能力、滯回曲線的捏攏效果上，基本能夠很好地描述試驗結(jié)果.而其在cpu主頻為2.5 GHz的intel i5內(nèi)存8 G的筆記本電腦上“用戶運算時間”為110.40 s，運算效率應(yīng)該是令人非常滿意的.

8.2 算例2-邊框架節(jié)點UNIT3

分析中梁和上下柱均分10個單元，圖6為梁端荷載與位移滯回曲線的試驗[9]和本文分析結(jié)果.

從圖6可看出，在承載能力、初始剛度和捏攏效果上，本文方法都基本能夠準確描述，說明本單元在邊節(jié)點上有著比較好的表現(xiàn).

8.3 算例3-角框架節(jié)點KJ5

分析中梁柱均分10個單元，原試驗曲線是梁柱交界面上，梁截面彎矩與柱偏轉(zhuǎn)角之間的滯回曲線[10]，本文將該滯回曲線轉(zhuǎn)換為梁加載端剪力與軸向位移之間的滯回曲線.梁端受拉力時，角框架變形為“張開”，對應(yīng)滯回曲線上的正位移與正剪力；梁端受壓力時，角框架變形為“閉合”，對應(yīng)滯回曲線上的負位移與負剪力.滯回曲線的本文方法分析結(jié)果和試驗結(jié)果對比見圖7.

圖5 UNIT1試驗與分析滯回曲線對比

圖6 UNIT3試驗與分析滯回曲線對比

圖7 KJ5試驗與分析滯回曲線對比

從圖7可看出，在承載能力上，本文方法基本能準確描述.“張開”變形的強度描述基本準確，而在“閉合”變形后期加載時，結(jié)構(gòu)強度嚴重退化，本文方法雖有反映，但退化程度與實際差距較大.

9 結(jié) 語

本文將鋼筋混凝土梁柱節(jié)點單元的內(nèi)節(jié)點自由度“依附”到外節(jié)點上形成一種“超自由度”模式的單元，通過編制單元的UEL程序?qū)⑵湟浦驳紸BAQUS通用分析平臺上，使得循環(huán)荷載作用下的鋼筋混凝土框架結(jié)構(gòu)的非線性分析可以真正在通用數(shù)值分析平臺得以進行，為鋼筋混凝土框架結(jié)構(gòu)在通用平臺上進行動力荷載作用下的非線性分析打下了堅實基礎(chǔ).通過對鋼筋混凝土結(jié)構(gòu)內(nèi)、邊和角框架節(jié)點3個梁柱節(jié)點模型試驗的算例分析，也驗證了本文“超自由度”節(jié)點單元的有效性.

[1] MITRA N, LOWES L N. Evaluation, calibration, and verification of a reinforced concrete beam-column joint model [J]. Journal of Structural Engineering, ASCE, 2007, 133(1): 105-120. DOI:10.1061/(ASCE)0733-9445(2007)133:1(105).

[2] LOWES L N, ALTOONTASH A. Modeling reinforced concrete beam-column joints subjected to cyclic loading [J]. Journal of Structural Engineering, ASCE, 2003, 129(12): 1686-1697. DOI: 10.1061/(ASCE)0733-9445(2003)129:12(1686).

[3] VECCHIO F J, COLLINS M P. The modified compression field theory for reinforced concrete elements subjected to shear [J]. ACI Structural Journal, 1986, 83(2): 219-231.

[4] LAFAVE J M, KIM J R. Joint shear behavior prediction for RC beam-column connections [J]. International Journal of Concrete Structures and Materials, 2011, 5(1): 57-64.DOI:10.4334/IJCSM.2011.5.1.057.

[5] 方自虎,謝強,周堯,等.循環(huán)荷載下MCFT混凝土結(jié)構(gòu)模型[J].地震工程與工程振動,2015,35(3):73-78. DOI:10.13197/j.eeev.2015.03.073.fangzh.009. FANG Zihu, XIE Qiang, ZHOU Yao, et al. MCFT model of concrete structures subjected to cyclic loading [J]. Earthquake Engineering and Engineering Dynamics,2015,35(3):73-78. DOI:10.13197/j.eeev.2015.03.073.fangzh.009.

[6] 方自虎,簡旭陽,周堯,等.考慮屈曲的鋼筋滯回模型[J].武漢大學(xué)學(xué)報,2016,49(2):254-258. DOI:10-14188/j.1671-8844.2016-02-016. FANGZihu, JIAN Xuyang, ZHOU Yao, et al. Hysteretic model of reinforced bar considering buckling[J]. Engineering Joural of Wuhan University,2016, 49(2):254-258. DOI:10-14188/j.1671-8844.2016-02-016.

[7] SHIOHARA H. New model for shear failure of RC interior beam-column connections [J]. Journal of Structural Engineering, ASCE, 2001, 127(2): 152-160.

[8] PARK R, RUITONG D. A comparison of the behavior of reinforced concrete beam-column joints designed for ductility and limited ductility [J]. Bulletin of the New Zealand National Society of Earthquake Engineering, 1988, 21(4): 255-278.

[9] PARK R, MILBURN J R. Comparison of recent New Zealand and United States seismic design provisions for reinforced concrete beam-column joints and test results from four units designed according to the New Zealand code [J]. Bulletin of the New Zealand national society for earthquake engineering, 1983, 16(1): 3-24.

[10]ANGELAKOS B. The behavior of reinforced concrete knee joints under earthquake loads [D]. Toronto: University of Toronto, 1999.

Super DOF element of RC beam-column joint

FANG Zihu, LI Xiangpeng, JIAN Xuyang, HONG Bokai

(College of Civil Engineering，Shenzhen University，Shenzhen 518060, Guangdong，China)

To analyze the nonlinear hysteretic behaviors of RC frame structures under reversed cyclic loading, a new element of reinforced concrete beam-column joints, super DOF element, was presented. The element edges were divided into “joint plane” and “beam-column plane” at the joint-column interface and joint-beam interface. The inelastic mechanism of joint core was represented by the four-node plane stress element. The anchorage failures of beam and column longitudinal reinforcement embedded in the joint were determined by eight springs between “joint plane” and “beam-column plane”. The proposed element has four exterior nodes and four interior nodes. There are two degrees of freedom on each interior node. There are three degrees of freedom on each exterior node, coinciding with ones of typical beam element, so the element is suitable for use together with typical hysteretic beam-column line elements in two-dimensional nonlinear analysis of reinforced concrete structures. The element is implemented as a four-node twenty-degree-of-freedom element through moving the degrees of freedom on the interior nodes to ones on the exterior nodes. The simulated data were compared with tests, and the results indicate that the ultimate strengths and hysteretic pitching behaviors analyzed by the element are well agreement with the test ones and the computer time is short. It concludes that the proposed element is suitable for use in simulating response of building joints under cyclic loading.

beam-column joint element; reinforced concrete; frame structures; reversed cyclic loading; non-linear analysis; hysteretic behaviors

(編輯 趙麗瑩)

10.11918/j.issn.0367-6234.201510063

2015-10-20

國家自然科學(xué)基金面上項目(51378313，51578336)

方自虎(1962—)，男，副教授，碩士生導(dǎo)師

方自虎，zihu-f@21cn.com

TU375；O242

A

0367-6234(2017)06-0053-05