伍 劍，鄒 宏，徐志純，王坤云

(成都理工大學(xué)地質(zhì)災(zāi)害防治與地質(zhì)環(huán)境保護(hù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，四川 成都 610059)

CFRP鋼管混凝土是在鋼管混凝土外部纏繞碳纖維布的一種組合結(jié)構(gòu)，它具有承載力高、韌性好、制作和施工方便、經(jīng)濟(jì)效果好等優(yōu)點(diǎn)。目前，國(guó)內(nèi)外一些學(xué)者已開(kāi)展了相關(guān)研究。陶忠等[1]分別在圓形和矩形鋼管混凝土上測(cè)試了碳纖維增強(qiáng)復(fù)合材料的加固效果，并提出了一種計(jì)算CFRP包裹的圓形鋼管混凝土短柱極限強(qiáng)度的簡(jiǎn)單模型。王慶利等[2]對(duì)32個(gè)CFRP鋼管混凝土進(jìn)行了軸壓試驗(yàn)研究，將試件的載荷-中截面撓度曲線分為彈性段、彈塑性段和軟化段3個(gè)階段，并用纖維模型法模擬了圓CFRP鋼管混凝土軸壓構(gòu)件的荷載-中截面撓度曲線，計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合良好，且偏于安全。顧威等[3]采用極限平衡法，推導(dǎo)出鋼管屈服時(shí)構(gòu)件承載力和碳纖維片拉斷時(shí)構(gòu)件極限承載力的解析計(jì)算式，并通過(guò)對(duì)8根CFRP-鋼管混凝土軸壓短柱的極限承載力進(jìn)行試驗(yàn)研究，驗(yàn)證解析計(jì)算式的正確性。Li Na等[4]進(jìn)行了長(zhǎng)細(xì)比對(duì)碳纖維布約束鋼管混凝土柱性能的影響的試驗(yàn)研究。Jai Woo Park等[5]以D/t(B/t)比、CFRP層數(shù)為試驗(yàn)參數(shù)，分別對(duì)CFRP加固圓形截面和矩形截面的鋼管混凝土結(jié)構(gòu)的性能進(jìn)行了比較。王靜峰等[6]對(duì)軸壓作用下CFRP部分包裹圓鋼管混凝土短柱的力學(xué)性能進(jìn)行了參數(shù)分析。

上述研究表明：使用CFRP包裹鋼管混凝土可以顯著提高鋼管混凝土的承載力。然而，CFRP鋼管混凝土的極限承載力與很多因素有關(guān)，如鋼管厚度、混凝土強(qiáng)度、CFRP厚度等。傳統(tǒng)的極限承載力是通過(guò)理論計(jì)算、實(shí)驗(yàn)和有限元模擬得到。實(shí)際上，由于材料制造、使用等原因，CFRP鋼管混凝土存在大量的不確定因素，這些不確定因素導(dǎo)致CFRP極限承載力實(shí)際并不是“一個(gè)值”，而是“一組值”。在判斷不同的不確定參數(shù)對(duì)CFRP鋼管混凝土極限承載力的影響程度時(shí)，采用靈敏度指標(biāo)更具有實(shí)際意義。

本文以CFRP鋼管混凝土為研究對(duì)象，使用ABAQUS分析確定參量系統(tǒng)的極限承載力，并對(duì)比試驗(yàn)得到的極限承載力，驗(yàn)證模型的可行性。然后考慮試件參數(shù)的不確定性，進(jìn)一步探討不同隨機(jī)參量對(duì)極限承載力影響的局部靈敏度指標(biāo)和全局靈敏度指標(biāo)。研究結(jié)果可為CFRP鋼管混凝土的設(shè)計(jì)、優(yōu)化以及后期的加固提供參考和依據(jù)。

1 模型的建立及分析

1.1 材料本構(gòu)關(guān)系

某江流鋼管的本構(gòu)關(guān)系采用彈塑性模型中的五段式二次塑流模型[7]。其中鋼材的彈性模量為2.06×105MPa，泊松比為0.3。

核心混凝土的本構(gòu)關(guān)系采用文獻(xiàn)[7]提出的鋼管約束核心混凝土本構(gòu)關(guān)系模型。核心混凝土的抗壓應(yīng)力應(yīng)變曲線和抗拉應(yīng)力應(yīng)變曲線見(jiàn)圖1。其中混凝土軸心抗壓強(qiáng)度標(biāo)準(zhǔn)值為36.85 MPa，泊松比為0.2。

a)混凝土受壓應(yīng)力-應(yīng)變關(guān)系

CFRP只承受纖維方向的拉應(yīng)力，將其他方向的應(yīng)力值設(shè)為0.001 MPa，在斷裂失效前為線彈性材料，認(rèn)為纖維方向(鋼管環(huán)向)上CFRP的抗拉強(qiáng)度為斷裂強(qiáng)度，并滿足胡克定律：

σcf=Ecfεcf

(1)

式中σcf——CFRP的環(huán)向應(yīng)力；Ecf——CFRP的彈性模量；εcf——CFRP的環(huán)向應(yīng)變。

1.2 計(jì)算模型

CFRP鋼管混凝土模型見(jiàn)圖2，試件長(zhǎng)度為400 mm、鋼管內(nèi)徑(半徑)為62 mm、鋼管壁厚為2.5 mm、CFRP厚度為0.17 mm。核心混凝土、鋼管均采用C3D8R單元，計(jì)算精度稍低，但相對(duì)于高次等參單元來(lái)說(shuō)，它的自由度較少，可以節(jié)省計(jì)算成本的特點(diǎn)；CFRP采用M3D4R單元，它只傳遞面內(nèi)的力，不承受彎矩，即沒(méi)有彎曲剛度。

圖2 CFRP鋼管混凝土模型

采用掃掠網(wǎng)格劃分技術(shù)分別對(duì)核心混凝土、鋼管、CFRP進(jìn)行網(wǎng)格劃分，經(jīng)試算當(dāng)網(wǎng)格尺寸取為模型1/30長(zhǎng)度左右時(shí)，可在較短的時(shí)間內(nèi)取得滿意的計(jì)算效果。在CFRP鋼管混凝土柱模型中，存在多種接觸關(guān)系。核心混凝土與鋼管之間相互作用采用“表面對(duì)表面接觸”，法向行為定為“硬接觸”，切向行為定義為“罰”，摩擦系數(shù)設(shè)為0.6。CFRP與鋼管之間接觸方式設(shè)置為“綁定”。在計(jì)算模型中，底端和頂端各設(shè)一個(gè)參考點(diǎn)，分別與模型的底面和頂面進(jìn)行耦合。通過(guò)參考點(diǎn)將構(gòu)件底端完全固定，頂端施加位移荷載。

1.3 承載力分析

為驗(yàn)證上述建模方法的可行性和準(zhǔn)確性，試件的試驗(yàn)和模擬的荷載-應(yīng)變曲線見(jiàn)圖3，可以看出，有限元模擬所得的荷載-應(yīng)變曲線與文獻(xiàn)中試驗(yàn)得到荷載-位移曲線大致吻合，都具有相同的變化規(guī)律，在初始受力階段，曲線呈線性上升。隨著應(yīng)變的逐漸增加，荷載達(dá)到最大值(即極限承載力)，繼續(xù)施加荷載，曲線發(fā)生驟降，這是因?yàn)镃FRP發(fā)生了斷裂，之后曲線達(dá)到穩(wěn)定狀態(tài)。本次模擬得到的極限承載力為1 276 kN，試驗(yàn)中的極限承載力為1 294 kN，模擬結(jié)果略小于試驗(yàn)結(jié)果，相對(duì)誤差值為1.39%，表明本次有限元模擬過(guò)程較好。

圖3 試驗(yàn)曲線與模擬曲線對(duì)比

2 靈敏度分析的實(shí)現(xiàn)

2.1 參數(shù)化建模

Abaqus腳本接口可以實(shí)現(xiàn)Abaqus/CAE中的所有功能[12]。使用者可以通過(guò)圖形用戶界面(GUI)窗口、命令行接口(CLI)和腳本(SCRIPT)來(lái)執(zhí)行相應(yīng)的命令，不過(guò)所用的命令都必須經(jīng)過(guò)Python解釋器后才能進(jìn)入Abaqus/CAE中執(zhí)行，同時(shí)生成擴(kuò)展名為.rpy的文件。再經(jīng)過(guò)求解器分析，最終得到輸出數(shù)據(jù)庫(kù)(ODB)文件。

Python語(yǔ)言是Abaqus軟件的標(biāo)準(zhǔn)設(shè)計(jì)語(yǔ)言，在Abaqus中的所有操作，都以Python語(yǔ)言的形式保存在工作路徑下的.rpy文件中。通過(guò)修改.rpy文件，提取需要參數(shù)化的變量，并編寫(xiě)Python語(yǔ)句進(jìn)行自動(dòng)后處理，提取試件的承載力最大值。最后將修改后的.rpy文件保存名為“Input.py”的文件，即完成了腳本文件的創(chuàng)建。

2.2 靈敏度分析

文件批處理是實(shí)現(xiàn)腳本多次運(yùn)行的基礎(chǔ)。建立批處理文件“run_abaqus.bat”，其內(nèi)容為：abaqus cae noGUI=Input.py。通過(guò)批處理文件自動(dòng)調(diào)用Abaqus，運(yùn)行腳本文件“Input.py”進(jìn)行仿真計(jì)算。建立一個(gè)data.txt文件，用于存放需要計(jì)算的每一組參數(shù)值。提取每次計(jì)算的極限承載力(Max_RF)。分析各參數(shù)變化對(duì)極限承載力的影響程度，即靈敏度分析見(jiàn)圖4。

圖4 靈敏度分析流程

3 靈敏度分析方法

靈敏度分析主要是考察輸入變量對(duì)輸出不確定性的貢獻(xiàn)大小，對(duì)優(yōu)化模型具有重要的工程意義。本文中的靈敏度是指各參數(shù)變化對(duì)CFRP鋼管混凝土極限承載力的影響程度。靈敏度分析主要分為局部靈敏度和全局靈敏度。

3.1 局部靈敏度

局部靈敏度分析(LSA)是指模型響應(yīng)函數(shù)對(duì)輸入變量在名義值點(diǎn)的偏導(dǎo)[8]，研究的是各個(gè)變量參數(shù)的單獨(dú)的變化對(duì)承載力所造成的影響，見(jiàn)式(2)：

(2)

式中SLSAi——局部靈敏度指數(shù)；ΔF——極限承載力的相對(duì)變化大小；ΔXi——初始斷面各個(gè)輸入變量的相對(duì)變化大小；X0——均值向量。

由于不同的輸入?yún)?shù)通常具有不同的維數(shù)和尺度，因此很難對(duì)同一框架下的輸入?yún)?shù)對(duì)輸出響應(yīng)的貢獻(xiàn)進(jìn)行分析。因此定義局部相對(duì)靈敏度，即輸出響應(yīng)的相對(duì)變量率對(duì)輸入?yún)⒘康南鄬?duì)變化率。

(3)

式中SRLSAi——局部相對(duì)靈敏度指數(shù)；F0——輸入變量為X0時(shí)的值，二者不能同時(shí)為0。

3.2 全局靈敏度

相對(duì)局部靈敏度是無(wú)量綱的，只能反映輸入?yún)?shù)在均值處輸出的局部梯度信息。全局靈敏度分析提供了隨機(jī)參數(shù)的變化對(duì)目標(biāo)函數(shù)在整個(gè)參數(shù)空間內(nèi)變化的詳細(xì)影響，以及參數(shù)與目標(biāo)函數(shù)之間的相互作用。

全局靈敏度分析方法主要基于Monte Carlo方法，考慮了參數(shù)的實(shí)際概率分布，所有變量都可以同時(shí)變化，變動(dòng)范圍可以是整個(gè)設(shè)計(jì)空間。通過(guò)響應(yīng)的方差來(lái)衡量輸入?yún)?shù)的相對(duì)重要度[9]。然而，Monte Carlo法的一個(gè)主要缺點(diǎn)是需要對(duì)大量的模型輸出進(jìn)行評(píng)估，從而使靈敏度指數(shù)的最終逼近足夠準(zhǔn)確。考慮到本次模型較為復(fù)雜，如果基于Monte Carlo法進(jìn)行分析，則計(jì)算量非常龐大，可操作性大大受到限制，因此可以考慮借助非統(tǒng)計(jì)類(lèi)方法進(jìn)行研究。

(4)

式中R(x)——CFRP鋼管混凝土的極限承載力；X——輸入變量，二者不能同時(shí)為0；β——多項(xiàng)式的待定系數(shù)；k——輸入變量個(gè)數(shù)；ε——函數(shù)R(x)的誤差。

為了求解式中的待定系數(shù)，使用最優(yōu)拉丁超立方設(shè)計(jì)方法(Optimal Latin Hypercube Design)生成多個(gè)采樣點(diǎn)，最優(yōu)拉丁超立方設(shè)計(jì)改進(jìn)了隨機(jī)拉丁超立方設(shè)計(jì)的均勻性，使因子和響應(yīng)的擬合更加精確真實(shí)。并用確定性分析方法計(jì)算配點(diǎn)處的函數(shù)值，然后通過(guò)采樣點(diǎn)及其函數(shù)值擬合多項(xiàng)式中的待定系數(shù)，即可得到函數(shù)的具體表達(dá)式，最后通過(guò)表達(dá)式可進(jìn)一步計(jì)算系統(tǒng)輸出的統(tǒng)計(jì)特征值。

二階響應(yīng)面誤差見(jiàn)圖5，并結(jié)合響應(yīng)面的相關(guān)系數(shù)R2來(lái)評(píng)判擬合的精度，當(dāng)R2值越接近1，說(shuō)明擬合和預(yù)測(cè)效果越好[10]。本次通過(guò)二階響應(yīng)面擬合響應(yīng)的R2值為0.97。表明所建立的近似模型擬合程度較好，可用于代替實(shí)際模型并進(jìn)行全局靈敏度分析。其中，R2的表達(dá)式為：

圖5 二階響應(yīng)面近似模型誤差分析

(5)

在靈敏度分析方法方面，采用計(jì)算量較小，精度較高的基于方差分解的FAST法[11]計(jì)算CFRP鋼管混凝土的軸壓極限承載力的全局靈敏度指標(biāo)，分析各個(gè)參量對(duì)極限承載力的貢獻(xiàn)大小。參量Pi的FAST靈敏度指標(biāo)Si計(jì)算表達(dá)式為：

(6)

式中Si——全局靈敏度指數(shù)；E——均值；V——方差。

4 承載力的靈敏度分析

4.1 局部靈敏度分析

在局部靈敏度分析中，將CFRP鋼管混凝土長(zhǎng)度、鋼管內(nèi)半徑、鋼管壁厚、CFRP厚度、混凝土抗壓強(qiáng)度、鋼材強(qiáng)度、CFRP強(qiáng)度、混凝土彈性模量、鋼管彈性模量、CFRP彈性模量作為輸入?yún)?shù)，使其在各自均值的±15%范圍內(nèi)變化，參數(shù)均值見(jiàn)表1。

表1 參數(shù)的統(tǒng)計(jì)特性

各參數(shù)對(duì)承載力的局部相對(duì)靈敏度見(jiàn)圖6。圖中，斜率越大表明參數(shù)對(duì)響應(yīng)的影響程度越大。斜率為正說(shuō)明參數(shù)對(duì)響應(yīng)呈正相關(guān)，即隨著參數(shù)相對(duì)變化率的增加，響應(yīng)的相對(duì)變化率隨之增大，反之亦然。本次分析中，鋼管內(nèi)徑R1對(duì)響應(yīng)的影響最大并呈正相關(guān)，試件長(zhǎng)度L呈負(fù)相關(guān)。

圖6 極限承載力的相對(duì)局部靈敏度

表2給出局部靈敏度指標(biāo)和全局靈敏度指標(biāo)。其中第二列為各參數(shù)的局部相對(duì)靈敏度指數(shù)，靈敏度指數(shù)較大的前5個(gè)參數(shù)排序?yàn)椋轰摴軆?nèi)徑(R1)>鋼管強(qiáng)度(fy)>鋼管壁厚(ts)>混凝土立方體抗壓強(qiáng)度(fcu)>CFRP極限抗拉強(qiáng)度(fcf)。其中R1的靈敏度指數(shù)最大，表明R1對(duì)響應(yīng)影響最大；L的靈敏度指數(shù)為負(fù)，表明L對(duì)響應(yīng)的影響呈負(fù)相關(guān)。這與圖5得出的結(jié)論一致。

4.2 全局靈敏度分析

考慮到各參數(shù)的實(shí)際分布情況[13-15]，均服從正態(tài)分布，其均值和變異系數(shù)見(jiàn)表2，對(duì)CFRP鋼管混凝土進(jìn)行全局靈敏度分析。

各參數(shù)對(duì)承載力的FAST靈敏度直方圖見(jiàn)圖7，利用式(4)共計(jì)計(jì)算了3 221個(gè)數(shù)據(jù)點(diǎn)計(jì)算得到結(jié)果。從圖中可以看出，鋼管內(nèi)徑R1對(duì)響應(yīng)的影響最大。表2中的第三列給出了各參數(shù)的全局(FAST)靈敏度指數(shù)，其中靈敏度指數(shù)較大的前5個(gè)參數(shù)排序?yàn)椋轰摴軆?nèi)徑(R1)>混凝土立方體抗壓強(qiáng)度(fcu)> 鋼管強(qiáng)度(fy)>鋼管壁厚(ts)>CFRP極限抗拉強(qiáng)度(fcf)。

表2 局部-全局靈敏度指數(shù)對(duì)比

圖7 極限承載力的FAST靈敏度指標(biāo)的直方

對(duì)比局部靈敏度分析和全局靈敏度分析，可以發(fā)現(xiàn)，2種分析方法中都是鋼管內(nèi)徑對(duì)極限承載力的影響最大。但影響較大的前5個(gè)參數(shù)的重要性排序并不完全一致，在局部靈敏度分析中，混凝土立方體抗壓強(qiáng)度的靈敏度指數(shù)低于鋼管強(qiáng)度和鋼管壁厚的靈敏度指數(shù)，而在全局靈敏度分析中，混凝土立方體抗壓強(qiáng)度的靈敏度指數(shù)高于鋼管強(qiáng)度和鋼管壁厚的靈敏度指數(shù)。這是混凝土立方體抗壓強(qiáng)度的變異系數(shù)遠(yuǎn)遠(yuǎn)大于其他參數(shù)的變異系數(shù)引起的。因此，基于參數(shù)的變異性和實(shí)際概率分布，且多個(gè)參數(shù)同時(shí)發(fā)生變化的全局靈敏度分析可以更準(zhǔn)確地找出對(duì)目標(biāo)函數(shù)最為敏感的設(shè)計(jì)變量。但全局靈敏度分析的計(jì)算量十分龐大，若各參數(shù)的變異性相差不大或參數(shù)都在小范圍內(nèi)變化，使用局部靈敏度分析更能快速得到理想的結(jié)果。

5 結(jié)語(yǔ)

a)使用ABAQUS有限元分析軟件建立CFRP鋼管混凝土軸壓短柱數(shù)值分析模型，提取荷載-應(yīng)變曲線并于試驗(yàn)曲線對(duì)比，驗(yàn)證模型的準(zhǔn)確性。通過(guò)ABAQUS的腳本語(yǔ)言Python，考慮了9個(gè)參數(shù)并生成參數(shù)的樣本點(diǎn)，對(duì)CFRP鋼管混凝土進(jìn)行參數(shù)化建模，使用ABAQUS計(jì)算并提取CFRP鋼管混凝土的極限承載力。分析了極限承載力的局部靈敏度和全局靈敏度。在全局靈敏度分析中，針對(duì)Monte Carlo法分析復(fù)雜模型耗時(shí)長(zhǎng)的問(wèn)題，通過(guò)二階響應(yīng)面法擬合多項(xiàng)式函數(shù)，從而進(jìn)行全局靈敏度分析。

b)對(duì)比局部靈敏度分析和全局靈敏度分析，鋼管內(nèi)徑對(duì)極限承載力的影響最大，但是靈敏度指數(shù)大小排序并不一致。局部靈敏度考慮的是各參數(shù)在相同變化區(qū)間內(nèi)的梯度變化，每次只有一個(gè)參數(shù)變化；而全局靈敏度分析考慮了參數(shù)的變異性和實(shí)際概率分布，且多個(gè)參數(shù)同時(shí)發(fā)生變化。在各參數(shù)的變異性相差較小或參數(shù)都在小范圍內(nèi)變化時(shí)，使用局部靈敏度分析能快速得到理想的結(jié)果；若各參數(shù)變異性相差較大，使用考慮了參數(shù)變異性和實(shí)際概率分布的全局靈敏度分析更能得到準(zhǔn)確的結(jié)果。