唐昌輝,倪 佳,程明慧

(湖南大學(xué) 土木工程學(xué)院，湖南 長(zhǎng)沙 410082)

?

無黏結(jié)部分預(yù)應(yīng)力砼梁恢復(fù)力模型的研究

唐昌輝?,倪 佳,程明慧

(湖南大學(xué) 土木工程學(xué)院，湖南 長(zhǎng)沙 410082)

基于OpenSEES軟件的二次開發(fā)平臺(tái)，對(duì)無黏結(jié)部分預(yù)應(yīng)力混凝土梁的截面恢復(fù)力模型進(jìn)行了研究.根據(jù)試驗(yàn)統(tǒng)計(jì)和理論分析，建立了無黏結(jié)部分預(yù)應(yīng)力混凝土梁的截面彎矩-曲率滯回模型，并通過編制Visual C++程序，將其植入到OpenSEES軟件的材料子類中，用于無黏結(jié)部分預(yù)應(yīng)力混凝土結(jié)構(gòu)的非線性計(jì)算，分析過程中無需迭代計(jì)算無黏結(jié)預(yù)應(yīng)力筋的應(yīng)力增量，解決了無黏結(jié)部分預(yù)應(yīng)力混凝土結(jié)構(gòu)受力分析的關(guān)鍵問題.利用OpenSEES軟件，基于已建立的無黏結(jié)部分預(yù)應(yīng)力混凝土梁的截面恢復(fù)力模型，選用基于柔度法的非線性梁柱單元，采用位移控制，計(jì)算了低周反復(fù)荷載作用下無黏結(jié)部分預(yù)應(yīng)力混凝土梁的彎矩-撓度滯回曲線，計(jì)算曲線與試驗(yàn)曲線吻合良好，說明已建立的無黏結(jié)部分預(yù)應(yīng)力混凝土梁的截面恢復(fù)力模型的可靠性，為無黏結(jié)部分預(yù)應(yīng)力混凝土結(jié)構(gòu)的非線性計(jì)算提供了理論依據(jù).

無黏結(jié)部分預(yù)應(yīng)力砼梁；低周反復(fù)荷載；OpenSEES；恢復(fù)力模型；彎矩-撓度滯回曲線

無黏結(jié)預(yù)應(yīng)力砼梁的恢復(fù)力模型研究一直是工程界和學(xué)術(shù)界共同關(guān)注的問題.杜拱辰[1]基于材料的本構(gòu)模型，采用彎矩-曲率分析方法，假定無黏結(jié)筋長(zhǎng)度的增長(zhǎng)與周圍砼長(zhǎng)度變化總和相等，不斷迭代求解無黏結(jié)預(yù)應(yīng)力筋的極限應(yīng)力，該方法在計(jì)算破壞狀態(tài)情況下是可以實(shí)現(xiàn)，但對(duì)于無黏結(jié)預(yù)應(yīng)力砼梁在其他受力狀態(tài)就難以實(shí)現(xiàn)，不能求得無黏結(jié)預(yù)應(yīng)力砼梁的彎矩-撓度滯回曲線；唐昌輝等人[2-3]選用彎矩-曲率分析方法，將梁劃分為更為精細(xì)的N段，基于材料的本構(gòu)關(guān)系，通過迭代梁的撓曲線和無黏結(jié)預(yù)應(yīng)力筋的合力，求解梁的跨中彎矩-撓度曲線，該方法的優(yōu)點(diǎn)是可以計(jì)算無黏結(jié)預(yù)應(yīng)力砼梁的受力全過程彎矩-撓度滯回曲線.上述方法均需進(jìn)行無黏結(jié)預(yù)應(yīng)力筋的應(yīng)力增量的迭代計(jì)算，計(jì)算程序編制較為復(fù)雜.此外，唐昌輝[4]等人根據(jù)試驗(yàn)彎矩-曲率擬合曲線，提出了無黏結(jié)預(yù)應(yīng)力砼梁的截面恢復(fù)力模型，并基于截面恢復(fù)力模型，通過虛梁法，分級(jí)加曲率求解梁的跨中彎矩-撓度曲線，該方法無需迭代計(jì)算無黏結(jié)預(yù)應(yīng)力筋的應(yīng)力增量，且能較準(zhǔn)確求得無黏結(jié)預(yù)應(yīng)力砼梁的受力全過程彎矩-撓度滯回曲線，但計(jì)算程序受截面形式、配筋等參數(shù)的制約，不具有重用性，每次計(jì)算需根據(jù)具體的結(jié)構(gòu)參數(shù)重新編程.

根據(jù)試驗(yàn)結(jié)果和平衡條件，可得到梁的骨架曲線特征點(diǎn)，并利用OpenSEES軟件開源的特點(diǎn)，編制了Visual C++程序，將截面彎矩-曲率滯回模型植入到OpenSEE二次開發(fā)平臺(tái)中，計(jì)算了無黏結(jié)預(yù)應(yīng)力砼梁在低周反復(fù)荷載作用下的彎矩-撓度滯回曲線.此方法無需迭代計(jì)算無黏結(jié)預(yù)應(yīng)力筋的應(yīng)力增量，同時(shí)利用OpenSEES面向?qū)ο蟮某绦蛟O(shè)計(jì)的優(yōu)勢(shì)，實(shí)現(xiàn)程序的重用性.用戶可基于OpenSEES平臺(tái)建立所需的有限元模型，并調(diào)用植入的截面恢復(fù)力模型的命令，進(jìn)行無黏結(jié)預(yù)應(yīng)力砼結(jié)構(gòu)的非線性分析，為OpenSEES解決無黏結(jié)預(yù)應(yīng)力砼結(jié)構(gòu)恢復(fù)力模型的研究提供了一條新的途徑.

1 無黏結(jié)部分預(yù)應(yīng)力砼梁的截面恢復(fù)力模型

已有試驗(yàn)結(jié)果表明，無黏結(jié)部分預(yù)應(yīng)力混凝土梁的彎矩-曲率滯回曲線，除了具有明顯的開裂、屈服、破壞等特征外，還具有壓彎構(gòu)件的滑移、捏縮現(xiàn)象，因此，文獻(xiàn)[4]建立了包含骨架曲線和滯回特性的截面彎矩-曲率恢復(fù)力模型，如圖1所示.

圖1 無黏結(jié)部分預(yù)應(yīng)力砼梁的截面恢復(fù)力模型Fig.1 Restoring force model of unbounded partially prestressed concrete beam

1.1 骨架曲線的確定

基本假定：

1)除預(yù)應(yīng)力筋應(yīng)變外，鋼筋混凝土截面滿足平截面假定；

2)截面開裂后忽略受拉區(qū)混凝土的抗拉作用；

3)非預(yù)應(yīng)力筋屈服前，混凝土與非預(yù)應(yīng)力筋處于彈性階段，其應(yīng)力-應(yīng)變關(guān)系服從虎克定律；

4)無黏結(jié)預(yù)應(yīng)力筋與混凝土之間無摩擦；

5)極限承載力計(jì)算時(shí)，受壓區(qū)混凝土的應(yīng)力圖形可按《混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范》(GB 50010-2010)[5]簡(jiǎn)化為等效的矩形應(yīng)力圖.

開裂點(diǎn)A(D)：

對(duì)照《無黏結(jié)預(yù)應(yīng)力混凝土結(jié)構(gòu)技術(shù)規(guī)程》(JGJ92)[6]，對(duì)稱布置無黏結(jié)預(yù)應(yīng)力混凝土梁的開裂彎矩如式(1)所示，考慮截面塑性變形，取截面剛度為0.85EcI0，開裂點(diǎn)曲率如式(2)所示：

(1)

(2)

式中：σpc為扣除預(yù)應(yīng)力損失后，由預(yù)加力在抗裂驗(yàn)算邊緣產(chǎn)生的混凝土預(yù)壓應(yīng)力；γ為截面塑性矩系數(shù)．

屈服點(diǎn)B(E)：

由平衡條件可得，

(3)

從中求得ξy，再對(duì)受壓區(qū)預(yù)應(yīng)力筋合力作用點(diǎn)取矩得到屈服彎矩的表達(dá)式如下：

(4)

截面的屈服曲率為：

(5)

極限點(diǎn)C(F)：

由平衡條件可得，

(6)

從中求得ξu，再對(duì)受壓區(qū)預(yù)應(yīng)力筋合力作用點(diǎn)取矩得到極限彎矩的表達(dá)式如下：

(7)

截面的極限曲率為：

(8)

如圖1所示，截面加載超過極限點(diǎn)，進(jìn)入負(fù)剛度階段，點(diǎn)T(G)的彎矩和曲率分別為：

(9)

另取kCT=-0.08Mcr/φcr，即

(10)

1.2 滯回規(guī)則特征點(diǎn)的確定

正(反)向定點(diǎn)R(M)的確定：

試驗(yàn)數(shù)據(jù)表明，正(反)向定點(diǎn)彎矩的絕對(duì)值隨綜合配筋指標(biāo)q0的增大而減小，通過回歸分析得到正(反)向定點(diǎn)彎矩的表達(dá)式:

(11)

正(反)向定點(diǎn)對(duì)應(yīng)的曲率大致與屈服點(diǎn)曲率相等，即正(反)向定點(diǎn)的曲率:

(12)

正(反)向捏攏點(diǎn)L(Q)的確定:

同樣對(duì)已有試驗(yàn)數(shù)據(jù)進(jìn)行擬合分析，得到正(反)向捏攏點(diǎn)彎矩的表達(dá)式:

(13)

正(反)向捏攏點(diǎn)曲率的表達(dá)式：

(14)

2 基于OpenSEES截面恢復(fù)力模型的程序?qū)崿F(xiàn)

2.1 程序?qū)崿F(xiàn)

根據(jù)OpenSEES的軟件架構(gòu)[7-10]，若將截面恢復(fù)力模型植入到Open SEES二次平臺(tái)中，只需在Material類基礎(chǔ)上派生出新的材料子類.Material類有3個(gè)派生的子類：Uniaxial Material類、NDMaterial類、Section Force Deformation類.其中，Uniaxial Material類用于定義一維材料模型，包括纖維截面的應(yīng)力-應(yīng)變曲線和單元截面的力-變形曲線；NDMaterial類用于定義多維材料模型，如，實(shí)體單元的應(yīng)力-應(yīng)變響應(yīng)；Section Force Deformation類同樣用于定義多維材料模型，如，壓彎構(gòu)件的力-變形關(guān)系或零單元的彎曲和剪切本構(gòu)關(guān)系的合成應(yīng)力.

本文基于UniaxialMaterial類派生出UPCMaterial類，通過編譯C++程序，植入對(duì)稱布置無黏結(jié)部分預(yù)應(yīng)力混凝土梁的截面恢復(fù)力模型：1)利用派生類的函數(shù)成員覆蓋基類的函數(shù)成員，并定義恢復(fù)力模型的構(gòu)造函數(shù)和私有變量，創(chuàng)建代碼文件UPC Material.h；2)通過編寫恢復(fù)力模型的骨架曲線與滯回規(guī)則，實(shí)現(xiàn)頭文件中的函數(shù)成員，創(chuàng)建代碼文件UPCMaterial.cpp；3)基于VS2010編譯平臺(tái)，將創(chuàng)建的代碼文件和其它類文件載入U(xiǎn)PC Material項(xiàng)目，生成UPC Material工程和動(dòng)態(tài)鏈接庫文件UPCMaterial.dll，并將其與結(jié)構(gòu)分析的Tcl命令流放到一個(gè)子目錄中；4)在OpenSEES中僅需通過Tcl命令流調(diào)用已經(jīng)植入的UPCMaterial材料，進(jìn)行結(jié)構(gòu)的非線性分析計(jì)算.

Tcl命令流具體如下：

uniaxialMaterial UPCmaterial $matTag $mom1p $rot1p $mom2p $rot2p $mom3p $rot3p $mom4p $rot4p $momR $rotR $TmomL $mom1n $rot1n $mom2n $rot2n $mom3n $rot3n $mom4n $rot4n $momM $rotM $TmomQ $q0；

$matTag為材料標(biāo)號(hào)；

$mom1p $rot1p為骨架曲線上正向開裂點(diǎn)的彎矩和曲率；

$mom2p $rot2p為骨架曲線上正向屈服點(diǎn)的彎矩和曲率；

$mom3p $rot3p為骨架曲線上正向極限點(diǎn)的彎矩和曲率；

$mom4p $rot4p為正向骨架曲線考慮截面進(jìn)入負(fù)剛度的彎矩和曲率；

$momR $rotR為正向定點(diǎn)的彎矩和曲率；

$TmomL為正向捏攏點(diǎn)的彎矩；

$mom1n $rot1n為骨架曲線上反向開裂點(diǎn)的彎矩和曲率；

$mom2n $rot2n為骨架曲線上反向屈服點(diǎn)的彎矩和曲率；

$mom3n $rot3n為骨架曲線上反向極限點(diǎn)的彎矩和曲率；

$mom4n $rot4n為反向骨架曲線考慮截面進(jìn)入負(fù)剛度的彎矩和曲率；

$momM $rotM為反向定點(diǎn)的彎矩和曲率；

$TmomQ為反向捏攏點(diǎn)的彎矩；

$q0為截面的綜合配筋指標(biāo) ．

2.2 彎矩-撓度曲線的確定

與文獻(xiàn)[4]的擬梁法不同,本文基于OpenSEES平臺(tái)，選用基于柔度法的非線性梁柱單元，且采用位移控制，模擬試驗(yàn)梁在低周反復(fù)荷載試驗(yàn)的彎矩-撓度滯回曲線.通過分級(jí)增加跨中撓度，得到當(dāng)前結(jié)構(gòu)位移下的單元變形和截面變形，再根據(jù)截面恢復(fù)力模型，得到相應(yīng)的截面柔度矩陣和截面抵抗力；通過單元內(nèi)部迭代，不斷消除截面的不平衡力和單元的殘余變形，可得到單元?jiǎng)偠染仃嚭徒Y(jié)構(gòu)剛度矩陣，從而求得當(dāng)前跨中撓度所對(duì)應(yīng)的荷載大小.具體過程如圖2所示[9-12].

3 基于截面恢復(fù)力模型的無黏結(jié)預(yù)應(yīng)力砼梁數(shù)值模擬

3.1 試驗(yàn)簡(jiǎn)介

本文選取唐昌輝[4]的2根對(duì)稱布置無黏結(jié)部分預(yù)應(yīng)力混凝土梁和MUGURUMA H的無黏結(jié)部分預(yù)應(yīng)力混凝土梁進(jìn)行數(shù)值模擬.試驗(yàn)梁的基本數(shù)據(jù)如表1所示，試驗(yàn)梁的配筋如圖3所示.

圖2 計(jì)算框圖Fig.2 Calculating frame表1 試驗(yàn)梁的基本數(shù)據(jù)Tab.1 Main parameters of tested beam

試件編號(hào)fcu/(N·mm-2)σpcσ'pc/(N·mm-2)ApA'p/mm2fyf'y/(N·mm-2)AsA's/mm2q0q'0AA-3039.68628011431434643802362360.1630.143AA-3135.78188861431434643802362360.1770.169U-2.038.1842404662140.210U-3.538.1862664662140.160

圖3 試驗(yàn)梁配筋圖Fig.3 Reinforcement of tested beam

3.2 數(shù)值模擬結(jié)果

無黏結(jié)部分預(yù)應(yīng)力砼梁選用基于柔度法的非線性梁柱單元，其單元截面特性選用基于上述恢復(fù)力模型的截面定義，非線性分析同樣分為2個(gè)荷載工況，第1種工況為力控制加載，將預(yù)應(yīng)力作用和自重施加到結(jié)構(gòu)上，計(jì)算結(jié)構(gòu)的初始應(yīng)變；第2工況為位移控制加載，以跨中節(jié)點(diǎn)為控制節(jié)點(diǎn)，進(jìn)行往復(fù)加載、卸載、反向加載，以模擬低周反復(fù)荷載試驗(yàn)，模擬結(jié)果如圖4所示.

從圖中可以看出，計(jì)算的彎矩-撓度滯回曲線與試驗(yàn)曲線較吻合.與基于材料本構(gòu)模型相比，基于截面恢復(fù)力模型對(duì)結(jié)構(gòu)進(jìn)行非線性分析無需進(jìn)行截面積分，可直接由恢復(fù)力模型得到截面柔度矩陣，兼顧了非線性計(jì)算的效率和精度.

(a)AA-30

(b)AA-31

(c)U-2.0

(d)U-3.5圖4 基于截面恢復(fù)力模型的彎矩-撓度滯回曲線Fig.4 Moment-deflection hysteretic curves based on restoring force model of section

4 結(jié) 論

1)利用OpenSEES開源的特點(diǎn)，植入了無黏結(jié)預(yù)應(yīng)力筋的截面彎矩-曲率恢復(fù)力模型，利用OpenSEES中基于柔度法的非線性梁柱單元，對(duì)無黏結(jié)預(yù)應(yīng)力砼梁的恢復(fù)力模型進(jìn)行了研究，獲取了基于OpenSEES解決無黏結(jié)預(yù)應(yīng)力砼結(jié)構(gòu)恢復(fù)力模型的一種新方法.

2)計(jì)算了4根無黏結(jié)預(yù)應(yīng)力砼結(jié)構(gòu)梁的彎矩-撓度滯回曲線，與試驗(yàn)曲線比較，結(jié)果吻合較好，驗(yàn)證了二次開發(fā)程序的可靠性.由于截面柔度矩陣可直接求得，無需截面積分，兼顧了非線性計(jì)算的效率和精度.

[1] 杜拱辰.部分預(yù)應(yīng)力混凝土[M].北京:中國建筑工業(yè)出版社,1990:89-94.

DU Gong-chen. Partial prestressed concrete [M]. Beijing: China Architecture and Building Press,1990:89-94. (In Chinese)

[2] 唐昌輝,倪佳,葉林. 無黏結(jié)部分預(yù)應(yīng)力混凝土梁非線性全過程滯回分析[J].湖南大學(xué)學(xué)報(bào)：自然科學(xué)版,2015,42(7):62-68.

TANG Chang-hui, NI Jia, YE Lin. Nonlinear full range hysteretic analysis on unbonded partically prestressed concrete beams[J]. Journal of Hunan University: Natural Sciences,2015,42(7):62-68.(In Chinese)

[3] 唐昌輝,易偉建,沈蒲生. 無黏結(jié)筋對(duì)構(gòu)件抗彎剛度的貢獻(xiàn)分析與研究[J].湖南大學(xué)學(xué)報(bào)：自然科學(xué)版,2001,28(5):75-79．

TANG Chang-hui, YI Wei-jian, SHEN Pu-sheng. The contribution to the flexural stiffness of the unbounded tendons in prestressed concrete beams[J]. Journal of Hunan University: Natural Sciences, 2001,28(5):75-79. (In Chinese)

[4] 唐昌輝,劉冬明,李雲(yún). 對(duì)稱布置無黏結(jié)預(yù)應(yīng)力混凝土梁恢復(fù)力模型研究[J].建筑結(jié)構(gòu),2015,45(4):68-72.

TANG Chang-hui, LIU Dong-min, LI Yun. Research on restoring force model of unbonded partically prestressed concrete beams with symmetrical tendons[J]. Building Structure,2015,45(4):68-72.(In Chinese)

[5] GB50010-2010 混凝土結(jié)構(gòu)設(shè)計(jì)規(guī)范[S]. 北京:中國建筑工業(yè)出版社,2010:37-38.

GB50010-2010 Code for design of concrete structures[S]. Beijing: China Architecture and Building Press, 2010: 37-38. (In Chinese)

[6] JGJ92-2004無黏結(jié)預(yù)應(yīng)力混凝土結(jié)構(gòu)技術(shù)規(guī)程[S]. 北京:中國建筑工業(yè)出版社,1993:13-23．

JGJ92-2004 Technical specification for concrete structures prestressed with unbounded tendons[S]. Beijing: China Architecture and Building Press, 1993: 13-23. (In Chinese)

[7] SCOTT M H, FENVES G L. How to introduce a new material into OpenSees[R]. Berkeley:PEER, University of California, Berkeley, 2001:1-23．

[8] MCKENNA F T. Object-oriented finite element programming: frameworks for analysis, algorithms and parallel computing[D]. Berkeley: University of California, Berkeley, College of Engineering,1997:11-45．

[9] FILIPPOU F C, TAUCER F F, SPACONE E. A fiber beam-column element for seismic response analysis of reinforced concrete structures[R]. Berkeley: Earthquake Engineering Research Center, College of Engineering, University of California, Berkeley, 1991:23-88.

[10]陳學(xué)偉,林哲.結(jié)構(gòu)非線性程序OpenSEES原理與實(shí)例[M].北京:中國建筑工業(yè)出版社,2014:15-36.

CHEN Xue-wei, LIN Zhe. Structural nonlinear analysis program OpenSEES theory and tutorial[M]. Beijing: China Architecture and Building Press, 2014:15-36. (In Chinese)

[11]陳學(xué)偉,韓小雷,孫思為.三種非線性梁柱單元的研究及單元開發(fā)[J].工程力學(xué),2011,28(S1):5-11.

CHEN Xue-wei, HAN Xiao-lei, SUN Si-wei. Research and development of three nonlinear beam-column elements[J]. Engineering Mechanics,2011,28(S1):5-11. (In Chinese)

[12]MUGURUMA H. Study on low-cycle fatigue strength of post-tensioned unbonded prestressed concrete beams[J]. CEB Bulletin D’information,1979,381-384.(In Japanese)

Research on Restoring Force Model of Unbonded Partially Prestressed Concrete Beams

TANG Chang-hui?, NI Jia，CHENG Ming-hui

(College of Civil Engineering, Hunan Univ, Changsha, Hunan 410082, China)

Based on the secondary development platform of OpenSEES software, this paper studied the restoring force model of unbonded partially prestressed concrete beam sections. According to the statistics of experiments and theoretical analysis, the moment-curvature hysteretic model of the section for unbonded partially prestressed concrete beams was established, and planted into Material class in the software of OpenSEES by compiling Visual C++ programs, which is used for the nonlinear analysis of unbonded partially prestressed concrete structure without iterative calculation of stress increment for unbonded prestressed tendon. The crucial problem of stress analysis for unbonded partially prestressed concrete structure was then solved. Moreover, under displacement control, this paper considered the nonlinear beam-column element based on flexibility to calculate the moment-deflection hysteretic curves of unbonded partially prestressed concrete beams under low cyclic reversed loading on the basis of established restoring force model. The calculation results are reasonably identical to experiment results, which indicates the reliability of established restoring force model, and provides theoretical basis for the nonlinear analysis of unbonded partially prestressed concrete structures.

unbonded partially prestressed concrete beams; low cyclic reversed loading; OpenSEES; restoring force model; moment-deflection hysteretic curves

1674-2974(2016)11-0055-06

2015-11-23

國家自然科學(xué)基金資助項(xiàng)目(51278179), National Natural Science Foundation of China(51278179)

唐昌輝(1964-)，男，湖南新寧人，湖南大學(xué)副教授，博士?通訊聯(lián)系人，E-mail: tchui2010@126.com

TU378.2

A