吳毅彬，金國芳

(同濟大學 結(jié)構(gòu)工程與防災(zāi)研究所,上海，200092)

增強纖維復(fù)合材料(FRP)的力學性能指標是 FRP應(yīng)用于土木工程領(lǐng)域中的重要參數(shù)之一。各項指標選擇的合理以及對其離散型的評價都會直接影響FRP材料工程應(yīng)用的安全性與經(jīng)濟性。由于鋼筋、混凝土等材料本身的離散性較小，性能較穩(wěn)定，對材料概率統(tǒng)計分布已有一定的研究，因此已廣泛應(yīng)用于設(shè)計規(guī)范中。目前關(guān)于FRP混凝土結(jié)構(gòu)的可靠度研究主要是在對FRP材料概率分布假設(shè)的基礎(chǔ)上進行，研究表明：非正態(tài)分布簡化為正態(tài)分布帶來的可靠度計算誤差達20%～30%[1]。FRP材料概率統(tǒng)計分布研究和討論較少的主要原因在于一方面FRP材性受生產(chǎn)工藝、施工質(zhì)量、測試方案、破壞形態(tài)等因素的影響較為明顯，另一方面是FRP生產(chǎn)廠商習慣于將這些數(shù)據(jù)視為專有，不對外公開，而這些數(shù)據(jù)又是至關(guān)重要的，數(shù)據(jù)不足可能導(dǎo)致假定的性能參數(shù)累積概率分布函數(shù)形式不同，這對FRP材料在土木工程中的應(yīng)用有較大影響。對比已有研究中所采用的分布類型[2-11]可知：FRP材料參數(shù)的概率分布類型差異性較大，大多分布類型采用正態(tài)分布作為概率分布，這主要是基于經(jīng)驗假設(shè)或是建立在鋼筋混凝土材料統(tǒng)計分布研究的基礎(chǔ)上，由于FRP材料的力學性能與土木工程中長期使用的具有明顯彈塑性受力特征的普通低碳鋼和混凝土有較大的差別，且FRP材料本身離散性較強，假設(shè)的正確性有待驗證。同時已有統(tǒng)計參數(shù)主要針對材料的拉伸強度，由于纖維復(fù)合材的破壞形態(tài)有脆性斷裂、層間剝離等，因此有必要對其他指標，如彈性模量、厚度等參數(shù)進行分析。在此，本文作者通過收集整理近1 000份的FRP材性試驗樣本，運用參數(shù)估計和假設(shè)檢驗方法對材料性能各指標的統(tǒng)計參數(shù)和概率分布進行研究，并分析各參數(shù)的相關(guān)性，在此基礎(chǔ)上，重新對指標標準值進行定義、計算，為進一步研究基于可靠度的FRP加固混凝土結(jié)構(gòu)設(shè)計方法提供一定理論依據(jù)。

1 FRP材性試驗研究成果

纖維復(fù)合材應(yīng)用于土木工程階段，國內(nèi)外對材料的力學性能指標展開較多的研究[12-17]，但研究結(jié)果表明，F(xiàn)PR材料強度與破壞呈現(xiàn)離散性、多層次性、各向開性等特點，影響FRP材料的力學性能不僅與纖維布和膠結(jié)樹脂緊密相關(guān)，且受到施工工藝、工藝水平、施工環(huán)境等因素的影響，材料性能離散性較大，各指標統(tǒng)計參數(shù)、分布類型較鋼材等材料復(fù)雜。

本文所選用FRP材料性能測試數(shù)據(jù)均為根據(jù)國內(nèi)外標準試驗方法所進行的試驗，部分試驗數(shù)據(jù)由于試驗條件如試驗環(huán)境、儀器等的不同雖存在一定測試偏差，但對影響試驗結(jié)果的主要條件上基本保持一致性。收集到的國內(nèi)外近 1 000份的 FRP材性試驗樣本見表1。

由各指標的統(tǒng)計描述(表1)可知：各組樣本的性能指標變化幅度較大，個別數(shù)組中最大值為最小值的 2倍多。材料極限抗拉強度隨材料體系的不同有明顯變化，且與黏貼層數(shù)并非呈正比關(guān)系，層數(shù)越多，片材的拉伸強度反而降低。極限抗拉強度變異系數(shù)在0.03～0.23范圍浮動。彈性模量對不同材料與層數(shù)也呈現(xiàn)出較大差異，彈性模量隨著層數(shù)增加略有增長，變異系數(shù)在 0.03～0.28范圍內(nèi)。試件厚度為平均厚度，厚度小于單層厚度與層數(shù)的乘積，表明多層粘貼時存在明顯壓縮現(xiàn)象，厚度變異系數(shù)為0.03～0.13。

由上述3個指標的統(tǒng)計描述可看出：纖維復(fù)合材料，特別是現(xiàn)場濕貼型FRP，受纖維布本身質(zhì)量、樹脂性能、周邊環(huán)境及施工條件等因素的影響很大，因此，在設(shè)計階段應(yīng)考慮FRP材料的不確定因素，以保證工程安全可靠。

2 概率分布模型與參數(shù)估計

Normal分布、Lognormal分布、二參數(shù) Weibull分布和Gamma分布4種分布類型常用于工程中，且較適合于材料模型。每一種分布均可采用概率分布函數(shù)(PDF)與累積概率分布函數(shù)(CDF)進行描述。

(1) Normal分布

式中：均值-∞＜x，μ＜+∞；標準差σ＞0。

(2) Lognormal分布

式中：均值-∞＜μ＜+∞；標準差σ＞0；x＞0。

(3) 二參數(shù)Weibull分布

式中：β為形狀系數(shù)；η為比例系數(shù)；η，β＞0；x≥0。

(4) Gamma分布

式中：α為形狀系數(shù)；β為比例系數(shù)；α，β＞0；x≥0。

累積概率分布可由函數(shù)CDF圖進行描述，A1和D1組極限抗拉強度CDF圖如圖1所示。由圖1可知：4種分布類型均可適用于該數(shù)組，因此由CDF圖無法直觀選擇最合適的類型分布，其他數(shù)組也存在類似情況。

表1 FRP材料力學性能試驗樣本與參數(shù)統(tǒng)計Table 1 FRP data sets and descriptive statistics for mechanical property

圖1 A1和D1組極限抗拉強度累積概率分布圖Fig.1 Plot of CDF for set A1 and D1 tensile strength

3 參數(shù)分布類型擬合優(yōu)度檢驗

本文采用 3種擬合優(yōu)度檢驗方法，χ2檢驗(Chi-Square)、K-S檢驗(Kolmogorov-Smirnov)和 A-D檢驗(Anderson-Darling)對試驗樣本進行統(tǒng)計分析。由于不同分布類型對不同檢驗方法的接受標準不同，為了對比不同分布類型的擬合程度，以在同一顯著性水平α檢驗數(shù)組接受或拒絕檢驗的程度來衡量分布類型的擬合程度，χ2檢驗和K-S檢驗采用顯著性水平α=0.1進行檢驗；A-D檢驗采用α=0.25進行檢驗。

3.1 極限抗拉強度指標檢驗

由χ2檢驗結(jié)果如表2所示。由表2可知：各數(shù)組對4種分布類型均可接受檢驗，但從測試參數(shù)χq2看，數(shù)值越小表示分布類型接受檢驗的程度越高。表2中黑體為該數(shù)組的最優(yōu)分布類型，可見各組的最優(yōu)分布類型一致性不強，這主要是由于χ2檢驗將數(shù)值進行分組，容易導(dǎo)致差異較大值被忽略，而使假設(shè)通過檢驗。但從總體上看 Weibull分布較其他分布更具適用性。K-S檢驗與A-D檢驗結(jié)果也列于表2。由表2可知：Weibull分布在顯著性水平下檢驗接受程度高于其他類型分布。雖然 Weibull分布并非完全適用于所有數(shù)組，但與其他分布類型相比，Weibull分布的擬合效果最好，且在航空、航天工業(yè)中纖維復(fù)合材也已廣泛采用Weibull概率分布。因此，選取Weibull分布為材料極限抗拉強度的概率分布類型。

表2 極限抗拉強度擬合優(yōu)度檢驗結(jié)果Table 2 Goodness of fit test result for tensile strength

3.2 彈性模量指標檢驗

χ2檢驗結(jié)果表明：Weibull分布與Lognormal分布對彈性模量參數(shù)的擬合程度較好；K-S檢驗表明：Weibull分布的接受程度較Lognormal分布高；而A-D檢驗表明：Weibull分布與Gamma分布較優(yōu)。因此綜合3種檢驗方法結(jié)果，建議選取Weibull分布為材料彈性模量的概率分布類型。

3.3 厚度指標檢驗

χ2檢驗、K-S檢驗及A-D檢驗對材料厚度檢驗結(jié)果均表明：Lognormal分布的接受程度優(yōu)于其他分布類型，因此選取Lognormal分布為材料厚度的概率分布類型。

3.4 參數(shù)相關(guān)性分析

已有研究[2]中均視材料參數(shù)為相互獨立變量，然而若各參數(shù)間存在相關(guān)性，則會對可靠度分析產(chǎn)生影響。以ρ表示各參數(shù)之間的相關(guān)程度，當|ρ|＜0.2時，視各參數(shù)間不相關(guān)；當|ρ|＞0.8時，視各參數(shù)間完全相關(guān)。經(jīng)參數(shù)相關(guān)性分析可知：極限抗拉強度、彈性模量和厚度3個參數(shù)間呈現(xiàn)一定的弱相關(guān)性，相關(guān)系數(shù)在不同數(shù)組間變異性較大。極限抗拉強度與厚度、彈性模量與厚度呈負相關(guān)性，這是由于材料強度、彈模與應(yīng)力相對應(yīng)，厚度減小，材料橫截面積減小，在恒定拉力作用下，應(yīng)力增大。因此，高強度、高彈模與較薄厚度相對應(yīng)。同時，負相關(guān)性與廠商也有一定關(guān)系，越薄的材料代表制造水平越高，材料性能越好。

4 基于Weibull分布的強度標準值

4.1 已有FRP強度標準值

材料強度標準值的選取是進行設(shè)計分析的基礎(chǔ)，一般要求對至少20個樣本統(tǒng)計分析，并考慮環(huán)境、工藝、測試方法等因素，取概率分布某一分位值計算確定強度標準值，其計算表達式如下：

式中：xk為強度標準值；μk為樣本均值；σk為樣本標準差；n為常數(shù)。不同規(guī)范對常數(shù)n有不同。由于Weibull分布更合適描述拉伸強度指標，因此本文重新定義拉伸強度標準值。

4.2 標準值定義

Weibull概率密度函數(shù)中尺度參數(shù)η、形狀參數(shù)β與樣本均值μ和變異系數(shù)Ck有關(guān)。

采用式(7)和(8)計算參數(shù)η和β時，未考慮樣本數(shù)量限制而使參數(shù)η和β產(chǎn)生統(tǒng)計不確定性。因此本文采用最大似然數(shù)估計方法，對服從Weibul概率分布的η和β進行估計。假設(shè)樣本個數(shù)為n，其中x1，x2，…，xn是Weibull概率分布的觀察值，Weibull分布中參數(shù)η和β未知，Weibull分布似然函數(shù)為：

根據(jù)式(9)，采用Matlab編程計算參數(shù)η，β的估計量。由η，β估計量，根據(jù)分位值p確定拉伸強度標準值，本文按照可靠度統(tǒng)一理論取p=5%分位值，即P[X＜xp]=p，拉伸強度標準值xp為：

4.3 材料標準值計算與對比分析

根據(jù)上述推導(dǎo)分析，現(xiàn)以 A1數(shù)組為例計算拉伸強度標準值，已知A1組試樣個數(shù)n＝100；β=10.81；η=4 121；假設(shè)置信水平r＝0.95；p為5%分位值，則由此可得：x0.05為3.131 GPa。

將計算得到各組樣本的等效 Weibull分布拉伸強度標準值與3個規(guī)范[18-20]的規(guī)范值進行對比，結(jié)果如表 3所示，其中，Pd為標準值與規(guī)范值的對比度，Pd=[(xp-xk)/xk]×100%。由表3可知，在相同分位值下，等效 Weibull分布材料標準值較各規(guī)范建議值偏于安全。由于規(guī)范建議值并非基于驗證的概率分布，且系數(shù)的選取基本是以混凝土或鋼筋的統(tǒng)計分析為基礎(chǔ)，因此采用經(jīng)驗證的 Weibull分布及得到的標準值較規(guī)范取值具有較高的可靠性及適用性，同時也可為進一步研究基于可靠度的FRP加固混凝土結(jié)構(gòu)設(shè)計方法提供一定理論依據(jù)。

表3 等效Weibull分布拉伸強度標準值與規(guī)范值對比Table 3 Comparision of equivalent weibull characteristic values to current guideline characteristic values

5 結(jié)論

(1) 材料極限抗拉強度與彈性模量近似服從二參數(shù)Weibull分布；厚度近似服從對數(shù)正態(tài)分布。

(2) 3個指標間呈現(xiàn)一定的弱相關(guān)性，相關(guān)系數(shù)在不同數(shù)組間變異性較大。極限抗拉強度與厚度、彈性模量與厚度呈負相關(guān)性。

(3) 通過相同分位值下等效 Weibull分布拉伸強度標準值與現(xiàn)行規(guī)范取值的對比分析可知，現(xiàn)行規(guī)范對FRP拉伸強度標準值取值差異較大，而采用經(jīng)驗證的Weibull分布函數(shù)計算公式可較好估計材性標準值，計算結(jié)果較規(guī)范取值更為合理、安全。

[1]Alqam M,Bennett R M,Zureick A H. Three-parameter vs.two–parameter Weibull distribution for pultruded composite material properties[J]. Composite Structures,2002,58: 497-503.

[2]Pevris N,Triantafillou T C,Veneziano A. Reliability of RC members strengthened with CFRP laminates[J]. Journal of Structural Engineering,1995,121(7): 1037-1044.

[3]Okeil A M,Tawil E L,Shahawy M. Flexural reliability of reinforced concrete bridge girders strengthened with carbon fiber-reinforced polymer laminates[J]. Journal of Bridge Engineering,2002,7(5): 290-299.

[4]Neocleous K,Pilakoutas K,Waldron P. Structural reliability for fiber reinforced polymer reinforced concrete structures[C]//Proceedings of the International Symposium on Fiber Reinforced Polymer Reinforced for Reinforced Concrete Structures(FRPRCS-4). Baltimore,Maryland,Detroit: American Concrete Institute,1999: 65-74.

[5]Pham H B,Al-Mahaidi R. Reliability analysis of bridge beams retrofitted with fibre reinforced polymers[J]. Composite Structures,2008,82: 177-184.

[6]Atadero R,Lee L,Karbhari V M. Consideration of material variability in reliability analysis of FRP strengthened bridge decks[J]. Composite Structures,2005,70: 430-443.

[7]Pliakoutas K,Neocleous K,Guaandagnini M. Design philosophy issue of fiber reinforced polymer reinforced concrete[J]. Journal of Composites for Construction,2002,6(3): 154-161.

[8]何政,李光. 基于可靠度的FRP筋材料分項系數(shù)的確定[J].工程力學,2008,25(9): 214-221.HE Zheng,LI Guang. Reliability-based calibration of material partial factor of fiber reinforced polymer rod[J]. Engineering Mechanics,2008,25(9): 214-221.

[9]楊勇新. 碳纖維布加固混凝土梁的可靠度分析[J]. 建筑結(jié)構(gòu)學報,2008(S1): 88-92.YANG Yong-xin. Reliability analysis of RC concrete beams strengthened with CFRP sheets[J]. Journal of Building Structures,2008(S1): 88-92.

[10]李業(yè)績,岳清瑞,楊勇新,等. 加固修復(fù)用碳纖維布力學性能試驗值的研究[J]. 玻璃鋼/復(fù)合材料,2008,25(9): 9-12.LI Ye-ji,YUE Qing-rui,YANG Yong-xin,et al. Stduy on mechanical property testing values of carbon fiber sheets uesd in strengthening and rehabilitation[J]. FRP/CM,2008,25(9): 9-12.

[11]于明偉. 單一纖維FRP筋基本力學性能與混雜FRP筋的研制[D]. 哈爾濱: 哈爾濱工業(yè)大學土木工程學院,2005: 35-42.YU Ming-wei. Mechanics properties study on FRP rebars and development of hybrid FRP rebars[D]. Harbin: Harbin Institute of Technology. School of Civil Engineering,2005: 35-42.

[12]岳清瑞,曹勁松,楊勇新,等. 碳纖維布標準檢測方法中試件尺寸影響的試驗研究[J]. 工業(yè)建筑,2005,35(8): 1-4.YUE Qing-rui,CAO Jin-song,YANG Yong-xin,et al. Normative specimen of carbon fiber reinforced polymer sheets for checking property[J]. Industrial Construction,2005,35(8): 1-4.

[13]孫秀紅,徐向東,徐茂波. 影響 FRP基本力學性能的因素[J].山東建筑工程學院學報,2005,20(2): 18-23.SUN Xiu-hong,XU Xiang-dong,XU Mao-bo. Factors to influence fundamental mechanical behaviors of FRP[J]. Journal of Shandong University of Architecture and Engineering,2005,20(2): 18-23.

[14]朱紅,張繼文. 結(jié)構(gòu)加固過程中 CFRP布基本力學性能的檢測[C]//第二屆全國土木工程用纖維增強復(fù)合材料(FRP)應(yīng)用技術(shù)學術(shù)交流會論文集. 昆明: 2002: 324-329.ZHU Hong,ZHANG Ji-wen. Test on the main mechanical properties of FRP[C]//The 2nd National Conference on FRP Concrete Structures. Kunming,2002: 324-329.

[15]Lee L S. Monitoring and service life estimation of reinforced concrete bridge decks rehabilitated with externally bonded carbon fiber reinforced polymer (CFRP) composites[D]. San Doego: University of California. Department of Structural Engineering,2005: 126-128.

[16]Atadero R A,Karbhari V M. Source of uncertainty and design value for field-manufactured FRP[J]. Composite Structure,2009,89(5): 83-93.

[17]Zureick A,Scott D. Short-team behavior and design of fiber-reinforced polymer slender members under axial compression[J]. Journal of Composite for Construction,1997,2(1): 140-150.

[18]ACI440.12R－02. Guide for the design and construction of externally bonded FRP systems for strengthening concrete structures[S].

[19]TR-55. Design guidance for strengthening concrete structures using fiber composite materials[S].

[20]GB 50367—2006,混凝土結(jié)構(gòu)加固設(shè)計規(guī)范[S].GB 50367—2006,Design code for strengthening concrete structure[S].