有初應(yīng)力的FRP約束混凝土圓柱應(yīng)力—應(yīng)變分析型模型

潘　毅，吳曉飛，郭　瑞，蔡聯(lián)亨

(1.西南交通大學(xué) 土木工程學(xué)院，四川 成都　610031；2.西南交通大學(xué) 交通隧道工程教育部重點實驗室，四川 成都　610031；3.成都基準(zhǔn)方中建筑設(shè)計有限公司，四川 成都　610017)

纖維增強復(fù)合材料(Fiber Reinforced Polymer，F(xiàn)RP)外包加固混凝土柱，通過約束混凝土的側(cè)向膨脹，從而提高混凝土柱的抗壓強度和延性。國內(nèi)外學(xué)者對FRP約束混凝土柱的應(yīng)力—應(yīng)變關(guān)系開展了大量研究，提出了很多模型。文獻(xiàn)[1]指出，F(xiàn)RP約束混凝土柱的應(yīng)力—應(yīng)變關(guān)系模型可分為設(shè)計型模型與分析型模型兩類。設(shè)計型模型將FRP 約束混凝土作為一個整體，通過對試驗數(shù)據(jù)的回歸分析，得到應(yīng)力—應(yīng)變曲線上的關(guān)鍵點，從而建立本構(gòu)關(guān)系[2]。而分析型模型則將混凝土和FRP視為2個部分，通過二者的應(yīng)力平衡和應(yīng)變協(xié)調(diào)，考慮FRP對混凝土的約束作用，從而建立應(yīng)力—應(yīng)變關(guān)系，更適合做結(jié)構(gòu)的非線性分析。但相比設(shè)計型模型，分析型模型還比較少。其中，文獻(xiàn)[3]通過修正約束混凝土的體積應(yīng)變計算公式，采用數(shù)值迭代，建立了FRP約束混凝土的應(yīng)力—應(yīng)變關(guān)系；文獻(xiàn)[4]以八面體強度準(zhǔn)則提出了FRP約束混凝土峰值應(yīng)力的計算方法，然后采用增量方法建立了FRP約束混凝土圓柱的分析型模型；文獻(xiàn)[5-6]根據(jù)大量的試驗數(shù)據(jù)對FRP約束混凝土的軸向—側(cè)向應(yīng)變關(guān)系進(jìn)行了修正，并由此建立FRP約束混凝土圓柱的分析型模型。但上述模型都是混凝土柱在無初應(yīng)力的條件下提出的，而實際加固中的混凝土柱都帶有一定的初應(yīng)力。而初應(yīng)力會造成外包FRP的拉應(yīng)變滯后，從而降低FRP的側(cè)向約束強度。目前，盡管對有初應(yīng)力下FRP約束混凝土的研究也已開展[7-13]，但這些研究多數(shù)是設(shè)計型模型，且多數(shù)模型僅與自己的試驗數(shù)據(jù)進(jìn)行了比較，未建立考慮不同初應(yīng)力的分析型模型。

針對上述研究中的不足，在103個有初應(yīng)力的FRP約束混凝土圓柱試驗數(shù)據(jù)的基礎(chǔ)上，本文引入考慮不同初應(yīng)力水平的影響因子，提出有初應(yīng)力的峰值應(yīng)力和峰值應(yīng)變計算公式，根據(jù)FRP約束混凝土軸向—側(cè)向應(yīng)變關(guān)系，構(gòu)建側(cè)向膨脹系數(shù)計算公式，并基于Mander的箍筋約束混凝土本構(gòu)模型，建立考慮初應(yīng)力影響的FRP約束混凝土圓柱應(yīng)力—應(yīng)變分析型模型，并與試驗結(jié)果和其他既有模型進(jìn)行對比。

1　有初應(yīng)力的FRP約束混凝土圓柱峰值應(yīng)力與峰值應(yīng)變

FRP約束混凝土圓柱的峰值應(yīng)力fcc與峰值應(yīng)變εcc主要取決于側(cè)向約束強度fl(即FRP對混凝土提供的側(cè)向約束應(yīng)力)，其表達(dá)式為

(1)

其中，

εh.rup=kεεfrp

式中：Ef為FRP的彈性模量；tf為FRP的實際厚度；D為混凝土圓柱的直徑；εh.rup為FRP的有效極限拉應(yīng)變；εfrp為FRP的名義極限拉應(yīng)變；kε為由于混凝土開裂處應(yīng)力集中等原因?qū)е翭RP提前斷裂的折減系數(shù)。

不同的學(xué)者對kε取值有所不同。其中，Carey and Harries[14]建議取0.57～0.61，而Lam and Teng[6]建議取0.586，敬登虎[15]建議FRP約束混凝土圓柱取0.65。而Ozbakkaloglu 和Akin[16]則認(rèn)為kε不應(yīng)為一個定值，其取值與未約束混凝土的強度、FRP的彈性模量有關(guān)，其表達(dá)式為

kε=0.9-2.3fco×10-3-0.75Ef×10-6

(2)

式中：fco為未約束混凝土的峰值應(yīng)力。

試驗研究發(fā)現(xiàn)，初應(yīng)力會降低FRP約束混凝土柱的峰值應(yīng)力與峰值應(yīng)變，且降低幅度是初應(yīng)力的函數(shù)。本文收集到有初應(yīng)力的FRP約束混凝土圓柱軸壓試驗數(shù)據(jù)103個，初應(yīng)力水平大多在0.0～0.7之間，其具體信息見表1。

表1　試驗數(shù)據(jù)概況

通過對收集到的103個試驗數(shù)據(jù)進(jìn)行回歸分析，并引入峰值應(yīng)力影響因子kf，建立有初應(yīng)力的峰值應(yīng)力計算公式，其表達(dá)式為

(3)

其中，

式中：n為初應(yīng)力水平，即初應(yīng)力與未約束混凝土峰值應(yīng)力的比值。

同峰值應(yīng)力計算公式的建立方法，引入峰值應(yīng)變影響因子kε，建立有初應(yīng)力的峰值應(yīng)變計算公式，其表達(dá)式為

(4)

其中，

kε=(1-n8.4)2.67

式中：El為FRP側(cè)向約束剛度；εco為未約束混凝土的峰值應(yīng)變；Eco為未約束混凝土的彈性模量。

由表1中的構(gòu)件基本參數(shù)和式(3)和式(4)計算得到峰值應(yīng)力、峰值應(yīng)變，并與試驗結(jié)果比較，二者的誤差基本都在20%以內(nèi)，如圖1所示。其中，峰值應(yīng)力的預(yù)測效果整體優(yōu)于峰值應(yīng)變。這主要有2個原因，其一是應(yīng)變測量主要通過應(yīng)變儀或電阻應(yīng)變片獲得，測量精度較應(yīng)力差一些；其二是峰值應(yīng)變不僅與側(cè)向約束應(yīng)力有關(guān)，而且與FRP的種類有關(guān)，而峰值應(yīng)力基本與FRP的種類無關(guān)[17]。因此，多數(shù)既有計算模型對峰值應(yīng)力的預(yù)測均好于峰值應(yīng)變。

圖1　本文計算值與試驗值的比較

為了驗證峰值應(yīng)力和峰值應(yīng)變計算公式的準(zhǔn)確性，根據(jù)收集到的試驗數(shù)據(jù)，本文通過平均絕對誤差、(計算值/試驗值)的平均值和變異系數(shù)，對本文計算公式和有代表性的現(xiàn)有計算公式進(jìn)行評估，其中，平均絕對誤差、(計算值/試驗值)的變異系數(shù)定義如下

(5)

(6)

對峰值應(yīng)力和峰值應(yīng)變計算公式的評估結(jié)果見表2。由表2可知，本文計算公式的預(yù)測效果整體優(yōu)于其他的現(xiàn)有計算公式。

表2　峰值應(yīng)力和峰值應(yīng)變計算公式的比較

2　軸向—側(cè)向應(yīng)變關(guān)系

FRP約束混凝土研究的基礎(chǔ)問題之一是如何確定其軸向—側(cè)向應(yīng)變關(guān)系。國內(nèi)外學(xué)者對此問題進(jìn)行了大量的研究[5,9,18]。本文定義側(cè)向膨脹系數(shù)μ為側(cè)向應(yīng)變和軸向應(yīng)變之比，其表達(dá)式為

(7)

式中：εl為混凝土的側(cè)向應(yīng)變，其與混凝土環(huán)向應(yīng)變相等；εc為混凝土的軸向應(yīng)變。

在FRP與混凝土之間粘接牢固且無任何相對滑移的條件下，混凝土的側(cè)向應(yīng)變與FRP的環(huán)向拉應(yīng)變相等，即εh=εl。

已有試驗研究表明，加載初期由于軸向壓力較小，F(xiàn)RP基本不發(fā)揮作用，其側(cè)向膨脹系數(shù)μ類似于無約束混凝土的泊松比ν，約為0.2[19]；隨著σc和εc的增加，混凝土側(cè)向不斷膨脹，μ隨著增大，直至FRP側(cè)向約束作用被充分發(fā)揮出來，μ將逐漸趨于穩(wěn)定，其穩(wěn)定值主要由側(cè)向約束強度決定[20]。多數(shù)軸向—側(cè)向應(yīng)變關(guān)系模型未考慮初應(yīng)力的影響，而文獻(xiàn)[21]認(rèn)為側(cè)向膨脹系數(shù)μ的穩(wěn)定值μu不僅與FRP的側(cè)向約束強度有關(guān)，還與初應(yīng)力水平有關(guān)，其表達(dá)式為

(8)

根據(jù)公式(8)，并結(jié)合FRP約束混凝土柱軸向—側(cè)向應(yīng)變關(guān)系的發(fā)展規(guī)律，本文建立了有初應(yīng)力的側(cè)向膨脹系數(shù)μ計算公式，其表達(dá)式為

(9)

將式(9)與文獻(xiàn)[5]的Teng模型、文獻(xiàn)[11]的試驗數(shù)據(jù)進(jìn)行比較，如圖2所示。由圖2比較可知，本文給出的計算公式不僅比Teng模型的公式簡單，而且有較好的計算精度。

圖2　軸向—側(cè)向應(yīng)變計算模型的比較

3　有初應(yīng)力的分析型模型

Mirmiran認(rèn)為FRP約束混凝土在每個給定的側(cè)向約束作用下等同于1個主動約束狀態(tài)，即可先給1個側(cè)向約束應(yīng)力，基于主動約束的本構(gòu)關(guān)系，便可得該側(cè)向約束應(yīng)力下的軸向應(yīng)力和應(yīng)變點，然后采用增量迭代可得到更多的點，最后將不同側(cè)向約束應(yīng)力下的應(yīng)力和應(yīng)變點連起來[22]，便可得完整的FRP約束混凝土應(yīng)力—應(yīng)變曲線，如圖3所示。由于箍筋屈服后所提供的側(cè)向約束應(yīng)力不變，等同于主動約束狀態(tài)，因此本文采用Mander的箍筋約束混凝土本構(gòu)模型作為FRP約束混凝土的主動約束方程，其表達(dá)式為

(10)

式中：σc為混凝土的軸向應(yīng)力；r為混凝土的脆性系數(shù)，r=Ec/(Ec-fcc/εcc)；其余符號與前文相同。

圖3　分析模型計算示意

根據(jù)上述提出的峰值應(yīng)力和峰值應(yīng)變計算公式及側(cè)向膨脹系數(shù)μ的計算公式，基于Mander主動約束方程，采用增量迭代方法，建立有初應(yīng)力的應(yīng)力—應(yīng)變分析型模型，具體計算步驟如下。

(1)給定1個約束混凝土圓柱的軸向應(yīng)變εc，i(i=1，2，…)。

(2)通過描述軸向—側(cè)向應(yīng)變關(guān)系的式(9)，得到軸向應(yīng)變εc，i對應(yīng)的側(cè)向應(yīng)變εl，i。

(3)通過混凝土與FRP的應(yīng)變協(xié)調(diào)條件，即εh，i=εl，i，得到FRP的應(yīng)變εh，i。

(4)由式(1)計算應(yīng)變εh，i下FRP對混凝土的側(cè)向約束力fl，i。

(5)將fl，i和初應(yīng)力水平n代入式(3)—式(4)，計算當(dāng)前側(cè)向約束力fl，i下對應(yīng)的峰值應(yīng)力fcc和峰值應(yīng)變εcc。

(6)將fcc和εcc代入式(10)，計算本次εc，i對應(yīng)的σc，i，得到應(yīng)力—應(yīng)變曲線上的1個點。

(7)逐步加大εc，i，重復(fù)步驟(1)—步驟(6)，即可得到應(yīng)力—應(yīng)變關(guān)系曲線上更多的點。

(8)直至εh，i大于εh.rup，結(jié)束迭代過程，得到完整的應(yīng)力—應(yīng)變關(guān)系曲線。

本文在Visual Basic 6.0編譯器下開發(fā)了相應(yīng)的計算程序，實現(xiàn)了該方法。

4　分析型模型的驗證與比較

4.1　分析型模型的驗證

為驗證本文分析型模型的正確性，取與初應(yīng)力水平n=0.00～0.71對應(yīng)的試驗數(shù)據(jù)進(jìn)行計算，并將計算得到的曲線與試驗得到的曲線進(jìn)行比較，限于篇幅，本文僅列出部分比較結(jié)果，如圖4所示。

圖4　本文模型計算曲線與試驗曲線的比較

由圖4可知，除個別構(gòu)件的計算結(jié)果與試驗結(jié)果的誤差稍大外，大多數(shù)分析型模型的計算曲線與試驗曲線吻合較好，且計算值整體略小于試驗值，偏于安全。這說明本文建立的分析型模型可以較好地考慮不同初應(yīng)力對應(yīng)力—應(yīng)變曲線的影響。

4.2　分析型模型的比較

為了檢查本文分析型模型的適用性，將本文的分析型模型與文獻(xiàn)[2,3,18,21]中分析型或設(shè)計型模型和收集的試驗數(shù)據(jù)進(jìn)行比較，限于篇幅，本文僅列出初應(yīng)力水平n較高的結(jié)果，如圖5所示。相比既有模型，由于考慮了初應(yīng)力水平的影響，本文的分析型模型計算曲線與試驗曲線整體吻合較好，峰值應(yīng)力和峰值應(yīng)變與試驗值的誤差也較小，優(yōu)于部分既有模型。

圖5　本文模型與既有模型的比較

5　結(jié)　語

(1)考慮不同的初應(yīng)力水平，引入峰值應(yīng)力和峰值應(yīng)變影響因子，所建立的峰值應(yīng)力與峰值應(yīng)變的計算公式能較好考慮初應(yīng)力的影響，與試驗結(jié)果的誤差小于20%。

(2)建立了1個考慮初應(yīng)力影響的軸向—側(cè)向應(yīng)變關(guān)系的三折線計算模型，該模型能較好地預(yù)測FRP約束混凝土圓柱側(cè)向膨脹的變化規(guī)律。

(3)建立的分析型模型計算曲線與試驗曲線吻合較好，比其他的既有模型有更好的預(yù)測精度。

[1]TOGAY O，JIAN C L，THOMAS V. FRP-Confined Concrete in Circular Sections：Review and Assessment of Stress-Strain Models [J]. Engineering Structures，2013，49(4)：1068-1088.

[2]XIAO Y，WU H. Compressive Behavior of Concrete Confined by Carbon Fiber Composite Jackets [J]. Journal of Material Civil Engineering，2000，12(2)：139-146.

[3]陶忠，高獻(xiàn)，于清，等. FRP約束混凝土的應(yīng)力—應(yīng)變關(guān)系 [J]. 工程力學(xué)，2005，22(4)：187-195.

(TAO Zhong，GAO Xian，YU Qing，et al. Stress-Strain Relation of FRP-Confined Concrete [J]. Engineering Mechanics，2005，22(4)：187-195. in Chinese)

[4]王震宇，王代玉，呂大剛. FRP 約束混凝土圓柱應(yīng)力—應(yīng)變分析模型 [J]. 哈爾濱工業(yè)大學(xué)學(xué)報，2010，42(2)：200-206.

(WANG Zhenyu, WANG Daiyu, Lü Dagang. Stress-Strain Analytical Model for FRP-Confined Concrete Circular Column [J]. Journal of Harbin Institute of Technology，2010，42(2)：200-206. in Chinese)

[5]JIANG T，TENG J G. Analysis-Oriented Stress-Strain Models for FRP-Confined Concrete [J]. Engineering Structures，2007，29(11)：2968-2986.

[6]TENG J G，HUANG Y L，LAM L，et al. Theoretical Model for Fiber Reinforced Polymer-Confined Concrete [J]. Journal of Composites for Construction，2007，11(2)：201-210.

[7]潘毅，萬里，吳曉飛. 負(fù)載下CFRP約束混凝土應(yīng)力—應(yīng)變關(guān)系分析型模型 [J]. 西南交通大學(xué)學(xué)報，2015，50(3)：461-465.

(PAN Yi, WAN Li, WU Xiaofei. Analysis-Oriented Stress-Strain Model of CFRP Confined Concrete with Preload [J]. Journal of Southwest Jiaotong University，2015，50(3)：461-465. in Chinese)

[8]潘從建. 考慮二次受力的FRP約束混凝土圓柱軸壓性能的試驗研究 [D]. 北京：中國建筑科學(xué)研究院建筑結(jié)構(gòu)研究所，2007：28-37.

(PAN Congjian. Experiment Study on the Compressive Behavior of Circular Concrete Columns Confined with FRP Considering the Secondary Load [D]. Beijing：Institute of Building Structure，China Academy of Building Research，2007：28-37. in Chinese)

[9]金建平. CFRP帶載約束加固混凝土圓柱軸壓力學(xué)性能試驗研究 [D]. 沈陽：東北大學(xué)資源與土木工程學(xué)院，2008：33-58.

(JIN Jianping. Experimental Study on the Axial Mechanical Behavior of Preloaded Circular Concrete Columns with CFRP Confinement [D]. Shenyang: School of Resources and Civil Engineering，Northeastern University，2008：33-58. in Chinese)

[10]李趙怡. CFRP持載加固混凝土柱的有限元分析 [D]. 大連：大連理工大學(xué)土木水利學(xué)院，2009：40-43.

(LI Zhaoyi. Finite Element Analysis of Concrete Column Confined by CFRP under Sustained Load [D]. Dalian：School of Civil Engineering and Water Conservancy，Dalian University of Technology，2009：40-43. in Chinese)

[11]曹雙寅，潘毅，敬登虎，等. 負(fù)載下碳纖維布約束混凝土圓柱應(yīng)力—應(yīng)變關(guān)系的試驗研究與分析 [J]. 東南大學(xué)學(xué)報，2008，30(2)：788-793.

(CAO Shuangyin，PAN Yi，JING Denghu，et al. Experimental Research on Axially Loaded Circular Concrete Columns Confined by CFRP under Preload [J]. Journal of Southeast University，2008，30(2)：788-793. in Chinese)

[12]潘毅，曹雙寅，敬登虎，等. 負(fù)載下碳纖維布約束混凝土方柱應(yīng)力—應(yīng)變關(guān)系的試驗研究與分析 [J]. 土木工程學(xué)報，2009，42(1)：23-29.

(PAN Yi，CAO Shuangyin，JING Denghu，et al. Test and Analysis of the Axial Stress-Strain Relationship of Square-Section Concrete Columns Confined by CFRP under Preload [J]. Civil Engineering Journal，2009，42(1)：23-29. in Chinese)

[13]黃靖. 玻璃纖維聚合持載加固混凝土柱抗壓性能試驗研究及模擬 [D]. 大連：大連理工大學(xué)土木水利學(xué)院，2005：19-24.

(HUANG Jing. Experimental Study and Simulation on Compressive Behavior of Preloaded Strengthened Concrete Columns with GFRP [D]. Dalian：School of Civil Engineering and Water Conservancy，Dalian University of Technology，2005：19-24. in Chinese)

[14]CAREY S A，HARRIES K A. An Axial Behavior and Modeling of Small Medium and Large-Scale Circular Sections Confined with CFRP Jackets [J]. ACI Structural Journal，2005，102(4)：596-604.

[15]敬登虎. FRP約束混凝土的應(yīng)力—應(yīng)變模型及其在加固中的應(yīng)用研究 [D]. 南京：東南大學(xué)土木工程學(xué)院，2006：71-74.

(JING Denghu. Researches on Models of Stress-Strain and Application in Rehabilitation for Concrete Confined by FRP [D]. Nanjing：School of Civil Engineering，Southeast University，2006：71-74. in Chinese)

[16]OZBAKKALOGLU T，AKIN E. Behavior of FRP-Confined Normal- and High-Strength Concrete under Cyclic Axial Compression [J]. Journal of Composites for Construction，2012，16(4)：451-463.

[17]騰錦光，陳建飛，史密斯 S T，等. FRP加固混凝土結(jié)構(gòu) [M]. 李榮，騰錦光，顧磊，譯. 北京：中國建筑工業(yè)出版社，2005.

[18]LAM L，TENG J G. Design-Oriented Stress-Strain Model for FRP-Confined Concrete [J]. Construction and Building Materials, 2009，17(6)：471-489.

[19]MINMIRAN A，SHAHAWY M. Behavior of Concrete Columns Confined by Fiber Composite [J]. Journal of Structural Engineering，1997，123(5)：583-590.

[20]吳剛，呂志濤. FRP約束混凝土圓柱無軟化段時的應(yīng)力—應(yīng)變關(guān)系研究 [J]. 建筑結(jié)構(gòu)學(xué)報，2003，24(5)：1-9.

(WU Gang，Lü Zhitao. Study on the Stress-Strain Relationship of FRP-Confined Concrete Circular Column without a Stress-Softening Response [J]. Journal of Building Structures，2003，24(5)：1-9. in Chinese)

[21]潘毅，吳曉飛，曹雙寅，等. 長期軸壓下有初應(yīng)力的CFRP約束混凝土柱應(yīng)力—應(yīng)變關(guān)系分析 [J]. 土木工程學(xué)報，2016，49(9)：9-19.

(PAN Yi，WU Xiaofei，CAO Shuangying，et al. Analysis on the Stress-Strain Relationship of CERP Confined Concrete with Initial Stress under Long-Term Sustained Load [J]. China Civil Engineering Journal，2016，49(9)：9-19. in Chinese)

[22]MANDER J B，PRIESTLEY M J N. Theoretical Stress-Strain Model for Confined Concrete [J]. Journal of Structural Engineering, 1988，114(8)：1805-1826.

[23]WU Y，WANG L. Unified Strength Model for Square and Circular Concrete Columns Confined by External Jacket [J]. Journal of Structural Engineering，2009，135(3)：253-261.

[24]趙海波. CFRP加固高強混凝土柱子軸心受壓性能試驗研究 [D]. 大連：大連理工大學(xué)土木水利學(xué)院，2006：52-54.

(ZHAO Haibo. Test Study on Mechanical Performance of High-Strength Concrete Column Confined by CFRP under Axial Load [D]. Dalian：School of Civil Engineering and Water Conservancy，Dalian University of Technology，2006：52-54. in Chinese)