石 湘，張易飛，杜 穎，張 偉，黃一峰

(1.中國海洋大學(xué) 工程學(xué)院，山東 青島 266100；2.海洋石油工程股份有限公司 海洋工程技術(shù)服務(wù)分公司，天津 300452)

灌漿卡箍是一種導(dǎo)管架平臺水下構(gòu)件常用的加固方法，其中自應(yīng)力灌漿卡箍以其承載力大、高的缺陷容限和良好的疲勞性能在國外得到了廣泛應(yīng)用[1]。國內(nèi)的工程應(yīng)用很少，陳波等[2]針對HZ21-1A平臺水深101 m 處的節(jié)點裂紋進行了自應(yīng)力灌漿卡箍加固的設(shè)計和施工介紹。程凱旋等[3]對卡箍加固后的損傷導(dǎo)管架平臺結(jié)構(gòu)強度進行了校核和方法驗證?？紤]到傳統(tǒng)自應(yīng)力灌漿卡箍的水下安裝難度大、結(jié)構(gòu)復(fù)雜等問題，一種膨脹式自應(yīng)力灌漿卡箍由楊彬等[4]提出并進行了試驗測試。這種卡箍是在水泥漿中摻加一定比列的膨脹劑，灌漿后可以自動建立卡箍預(yù)應(yīng)力，這樣水下安裝時間和難度相對于傳統(tǒng)自應(yīng)力灌漿卡箍有較大降低。膨脹式自應(yīng)力灌漿卡箍通過防滑校核判斷卡箍的連接尺寸是否足夠，通過防撬校核確定螺栓型號，這些校核尚不能完全滿足卡箍的設(shè)計要求，還需要通過有限元分析進行進一步強度校核[5]。這種灌漿卡箍是一種典型的鋼/灌漿組合結(jié)構(gòu)，灌漿膨脹受限處于三向受壓狀態(tài)，力學(xué)參數(shù)尚無相關(guān)規(guī)范可以查詢，早期設(shè)計僅是按照線彈性模型用于有限元分析[6]，采用合適的灌漿本構(gòu)模型十分重要。下文采用ABAQUS軟件建模，針對膨脹式灌漿卡箍構(gòu)筑了三種灌漿本構(gòu)關(guān)系模型，測試了相關(guān)的灌漿材料參數(shù)，通過一個節(jié)點灌漿卡箍的算例，分析和比較了三種灌漿本構(gòu)模型的計算結(jié)果，從應(yīng)用角度為這種卡箍的有限元分析提供一些指導(dǎo)和建議。

1 灌漿本構(gòu)關(guān)系模型的選取

卡箍內(nèi)的灌漿是一種非常脆的材料，相比于抗壓能力，抗拉能力非常有限。根據(jù)灌漿的材料特性以及組合結(jié)構(gòu)混凝土強度分析的常用本構(gòu)模型，選取了三種灌漿本構(gòu)關(guān)系模型，即線彈性模型、莫爾庫侖模型和塑性損傷模型。

1.1 線彈性模型

線彈性材料的本構(gòu)關(guān)系服從廣義胡克定律，應(yīng)力應(yīng)變在加卸載過程中成線性關(guān)系。在模型中只需要定義彈性模量和泊松比即可。若作用載荷使灌漿內(nèi)各點的應(yīng)力應(yīng)變處在彈性階段，僅用線彈性本構(gòu)模型即可表達材料的力學(xué)行為。

1.2 莫爾庫侖模型

ABAQUS采用的莫爾庫侖模型是經(jīng)典準則的擴展，其流動勢函數(shù)在子午線平面上為雙曲函數(shù)，避免了由于經(jīng)典莫爾庫侖準則在屈服面的尖角而引起的塑性流動方向不唯一的現(xiàn)象[7]。

莫爾庫侖本構(gòu)模型大多應(yīng)用于以顆粒結(jié)構(gòu)為主的材料在單調(diào)荷載下的分析。在彈性模型中定義材料的彈性模量和泊松比，在彈塑性模型中定義材料的內(nèi)摩擦角、剪脹角和黏聚力。其形式簡單且參數(shù)容易測試，常用于巖土材料，也可以作為卡箍灌漿的本構(gòu)關(guān)系模型。

1.3 塑性損傷模型

在ABAQUS中塑性損傷模型的定義是以塑性為基礎(chǔ)，假定導(dǎo)致材料損傷破壞的是各向受壓和受拉。其模型在彈性段仍用彈性模量和泊松比來定義。在塑性段除了定義材料的受壓和受拉的應(yīng)力及對應(yīng)的非彈性應(yīng)變，還可以通過設(shè)置損傷因子來定義材料受壓和受拉時的損傷度[8]。如圖1所示，當(dāng)材料壓應(yīng)變超過極限應(yīng)變εc,e0后，材料由彈性段進入受壓強化段，再超過εc0后進入軟化段。當(dāng)材料拉應(yīng)變超過極限應(yīng)變εt0后由彈性段進入受拉軟化段。

圖1 塑性損傷模型的應(yīng)力應(yīng)變關(guān)系

塑性損傷模型設(shè)置包括了應(yīng)變硬化和軟化，從而考慮灌漿的壓碎和開裂，可以充分描述卡箍內(nèi)灌漿的各種響應(yīng)行為。

2 三種灌漿本構(gòu)關(guān)系模型的材料參數(shù)測定

2.1 灌漿材料介紹

對于膨脹式自應(yīng)力灌漿卡箍，灌漿使用的水泥型號是P.O42.5，膨脹劑類型是FEA100，膨脹劑摻量不超過18%。選取最常用的10%與15%兩種摻量進行測試，水灰比為0.45，則分為以下兩種灌漿型號：

1)膨脹劑摻量10%灌漿，膨脹劑∶水泥∶水=1∶10∶4.95；

2)膨脹劑摻量15%灌漿，膨脹劑∶水泥∶水=1∶6.67∶3.45。

按照GB/T 17671—1999，在無圍壓下測試灌漿的28天時期立方體抗壓強度，膨脹劑摻量10%和15%灌漿分別為48.3 MPa和40.5 MPa。由于灌漿在卡箍內(nèi)是三向受壓狀態(tài)，所以相關(guān)力學(xué)性能需要考慮預(yù)應(yīng)力(或稱圍壓)的影響。表1是實測的實際尺寸灌漿卡箍膨脹壓力的數(shù)據(jù)[9]，參考表中數(shù)據(jù)膨脹劑摻量10%灌漿取1.5～1.6 MPa、膨脹劑摻量15%灌漿取2.2 MPa作為圍壓值進行相關(guān)測試。

表1 膨脹式自應(yīng)力灌漿卡箍的膨脹壓力

2.2 彈性模量的測試

灌漿的彈性模量測試參照測試塑性混凝土彈性模量的方法，首先測試灌漿的棱柱體抗壓強度fc，再分別測試0.7倍和0.3倍抗壓強度下試件相應(yīng)的應(yīng)變ε0.7fc與ε0.3fc，根據(jù)應(yīng)力應(yīng)變關(guān)系得出彈性模量：

(1)

制作兩個灌漿棱柱體試件(截面40 mm×40 mm，長160 mm)，在相對的兩面各貼上一個應(yīng)變片，采用全橋方法測試試件在相應(yīng)抗壓強度下的應(yīng)變，如圖2所示。則膨脹劑摻量10%和15%灌漿的彈性模量測試值分別為13 892 MPa和12 601 MPa。

圖2 灌漿材料彈性模量的測試

2.3 莫爾庫侖模型的黏聚力和內(nèi)摩擦角測試

灌漿材料的黏聚力和內(nèi)摩擦角采用常規(guī)三軸壓縮試驗進行測試，根據(jù)莫爾庫侖失效理論通過不同圍壓下測試的試件抗壓強度求解。將水泥漿澆筑在直徑50 mm、高100 mm的圓柱模具中，并與模具一起在標準養(yǎng)護箱養(yǎng)護28天，然后取出試模。由于灌漿膨脹與模具結(jié)合緊密，使用萬能試驗機采用推出法進行試件的啟模，如圖3所示。

圖3 采用推出法進行Φ50×100 mm試件的啟模

參考表1的壓力值，對于膨脹劑摻量10%的灌漿試件，測試試件在圍壓1.2 MPa、1.5 MPa以及2.5 MPa下的抗壓強度；對于膨脹劑摻量15%的灌漿試件，測試試件在圍壓1.8 MPa、2.2 MPa以及3.0 MPa下的抗壓強度。圖4是圍壓下壓縮灌漿試件的情況，采用的三軸試驗機型號為DYS-1000。

圖4 圍壓下壓縮灌漿試件

以主應(yīng)力為水平軸、剪應(yīng)力為豎直軸，以((σ1+σ3)/2,0)為圓心、(σ1-σ3)/2為半徑(σ3為圍壓，σ1為相應(yīng)圍壓下的抗壓強度即主應(yīng)力峰值)，作莫爾應(yīng)力圓，求出這些莫爾應(yīng)力圓的公切線即莫爾包絡(luò)線。如圖5所示在主應(yīng)力和剪應(yīng)力坐標系中，黏聚力就是莫爾包絡(luò)線與剪應(yīng)力軸的截距c，內(nèi)摩擦角就是它與主應(yīng)力軸的夾角φ。表2是膨脹劑摻量10%灌漿的三軸壓縮測試結(jié)果，然后采用專用的分析軟件進行莫爾庫侖參數(shù)分析，可得膨脹劑摻量10%灌漿的黏聚力是9.789 MPa；內(nèi)摩擦角是49.414°。

圖5 莫爾庫侖準則的黏聚力與內(nèi)摩擦角

表2 三軸壓縮實驗膨脹劑摻量10%灌漿測試結(jié)果

表3是膨脹劑摻量15%灌漿的三軸壓縮測試結(jié)果。同樣方法分析可得膨脹劑摻量15%灌漿的黏聚力是9.605 MPa；內(nèi)摩擦角是44.589°。

表3 三軸壓縮實驗膨脹劑摻量15%灌漿測試結(jié)果

2.4 塑性損傷模型的抗拉強度和壓縮段應(yīng)力應(yīng)變曲線測試

工程上一般用巴西劈裂法測定灌漿的抗拉強度。劈裂法的基本原理是基于圓柱試件受相對徑向壓縮后的彈性理論解。作用在試件中心的最大拉應(yīng)力σt為：

(2)

式中：σt即為抗拉強度；P是試件破壞時的極限壓力；d是試件的直徑；t是試件的厚度。

按照GB 50266—2013，水泥漿澆筑在直徑50 mm、高50 mm 的圓柱模具中，試件的養(yǎng)護方法和啟模方法與2.3節(jié)相同。圖6是進行灌漿試件抗拉強度測試的情況，測試結(jié)果如表4所示，可以看出膨脹劑摻量10%的灌漿抗拉強度為1.97 MPa，膨脹劑摻量15%的灌漿抗拉強度為1.55 MPa。

圖6 灌漿試件的抗拉強度測試

表4 灌漿的抗拉強度測試結(jié)果

在2.3節(jié)測試圍壓下的灌漿試件壓縮過程記錄了10%膨脹劑摻量(圍壓1.5 MPa)和15%摻量(圍壓2.2 MPa)灌漿壓縮的應(yīng)力應(yīng)變曲線，如圖7所示。這些應(yīng)力應(yīng)變數(shù)據(jù)可用于灌漿塑性損傷模型受壓塑性段的數(shù)據(jù)輸入和計算。

圖7 圍壓下灌漿壓縮的應(yīng)力應(yīng)變曲線

3 膨脹式自應(yīng)力灌漿卡箍有限元分析算例

選取南海某導(dǎo)管架平臺中的受損K型節(jié)點(圖8)進行膨脹式灌漿卡箍的設(shè)計，該平臺位于水深為38米的海域，受到的載荷以波浪、海風(fēng)與海流為主。平臺極限載荷的選取如下：以南海百年一遇的波浪作為設(shè)計的極限波浪載荷，其波高14.3 m、周期10.4 s；選取與極限波浪載荷對應(yīng)的持續(xù)風(fēng)速作為風(fēng)載荷，風(fēng)速36.2 m/s；流載荷同樣選取極限波浪載荷對應(yīng)的海流流速，上層流速1.319 m/s，中層流速1.113 m/s，底層流速0.792 m/s。求解海洋平臺在X、Y與45°三個方向極限載荷組合工況下(波浪、風(fēng)和流同向)的受損節(jié)點位置未損壞情況下的各單元內(nèi)力作為卡箍的外載荷進行結(jié)構(gòu)設(shè)計?？紤]了波浪相位角對平臺靜力分析的影響。

圖8 ANSYS平臺模型及受損節(jié)點

受損K型節(jié)點位于水面以下7.8 m，弦桿尺寸為Φ609.6×19 mm，兩個撐桿尺寸均為Φ304.8×9.5 mm，假設(shè)3 628號單元因疲勞裂紋與弦桿斷開來設(shè)計卡箍。

3.1 節(jié)點灌漿卡箍的結(jié)構(gòu)設(shè)計

結(jié)構(gòu)設(shè)計包括尺寸初步設(shè)計、防滑校核、防撬校核等，參照王冬冬等[5]介紹的設(shè)計方法做簡要說明。卡箍鞍板的強度要保證法蘭板不被螺栓預(yù)緊力碾碎，鞍板材料選取Q315鋼，螺栓材料選取42CrMo合金鋼，灌漿厚度取為30 mm，材料為膨脹劑摻量10%灌漿。圖9是初步設(shè)計的節(jié)點卡箍結(jié)構(gòu)。

圖9 節(jié)點卡箍初步設(shè)計圖

防滑校核只需校核3 628單元，根據(jù)計算外載荷引起的3 628單元的最大滑動應(yīng)力為0.30 MPa(45°極限載荷工況)，而根據(jù)滑動應(yīng)力經(jīng)驗公式卡箍的許用滑動應(yīng)力為0.79 MPa，滿足防滑校核要求。

防撬校核需要校核節(jié)點所有單元處的螺栓，根據(jù)計算在外載荷和膨脹壓力聯(lián)合作用下，在3 598單元端部需提供最大螺栓預(yù)緊力134.045 kN，再經(jīng)計算螺栓有效直徑d>25.6 mm，則螺栓的型號選用M30，其有效直徑為26.7 mm，共40個。

3.2 卡箍有限元模型的建立

模型分為三部分：卡箍兩瓣鞍板、灌漿環(huán)和受損節(jié)點(圖10)。由于研究重點為灌漿環(huán)，為了簡化并未建立螺栓模型，而是將螺栓預(yù)緊力等效地加在法蘭板的印記面上。模型中受損撐桿端部與弦桿完全斷開?？ü俊⒐酀{環(huán)和受損節(jié)點均采用六面體實體單元C3D8R，網(wǎng)格尺寸設(shè)置為30 mm，模型共劃分為54 766個單元，88 710個節(jié)點。

圖10 ABAQUS節(jié)點灌漿卡箍有限元模型

3.2.1 三種灌漿本構(gòu)關(guān)系的設(shè)置

卡箍和受損節(jié)點鋼材的彈性模量設(shè)定為2.06×105MPa，泊松比設(shè)定為0.3。對于本設(shè)計中的10%膨脹劑摻量灌漿材料，分別設(shè)置如下：

1)線彈性模型：彈性模量和泊松比分別為1.39×104MPa和0.2。

2)莫爾庫侖模型：黏聚力為9.789 MPa，內(nèi)摩擦角為49.414°，剪脹角采用缺省值36.31°。

3)塑性損傷模型：由于抗拉強度測試結(jié)果是無圍壓下的測試值，在1.5 MPa 圍壓下考慮增加30%[10]，則10%膨脹劑摻量灌漿抗拉強度ft設(shè)為1.97×1.3=2.56 MPa。再根據(jù)GB 50010—2002混凝土結(jié)構(gòu)設(shè)計規(guī)范[11]計算得到受拉軟化段應(yīng)力應(yīng)變作為灌漿受拉軟化段數(shù)據(jù)(圖1)。受壓部分的強化和軟化段采用1.5 MPa圍壓下10%膨脹劑摻量灌漿測試得到的應(yīng)力應(yīng)變，選取圖7所示的εc,e0—εc0和εc0—εcu兩段的應(yīng)力應(yīng)變作為受壓強化和軟化段數(shù)據(jù)。然后采用以下兩式分別計算灌漿的非彈性應(yīng)變和損傷因子[8]：

(3)

(4)

然后把應(yīng)力、非彈性應(yīng)變和損傷因子輸入塑性損傷模型的材料設(shè)置。其他材料參數(shù)采用缺省值。

3.2.2 接觸的設(shè)置

模型各部件之間的接觸設(shè)置如下：在受損撐桿與灌漿環(huán)的接觸面上設(shè)置法向硬接觸、切向摩擦接觸的接觸對，切向摩擦系數(shù)設(shè)為0.4[12]；受損節(jié)點其它桿件與灌漿環(huán)接觸面的接觸設(shè)定為綁定。灌漿環(huán)與兩瓣卡箍鞍板的接觸也分別設(shè)定為綁定。在兩瓣卡箍法蘭板接觸面上設(shè)置法向硬接觸、切向摩擦接觸的接觸對，切向摩擦系數(shù)也設(shè)為0.4[12]。

3.2.3 約束和載荷的設(shè)置

參見圖11的參考點位置，約束參考點4，在其它3個參考點施加受損結(jié)構(gòu)桿件端部未損情況下的內(nèi)力載荷，表5和表6是施加的兩個卡箍危險工況載荷。工況1是引起卡箍最大螺栓撬力時管節(jié)點所受載荷，對應(yīng)平臺環(huán)境載荷是Y方向、波浪相位角為335°；工況2是引起受損撐桿最大滑動應(yīng)力時管節(jié)點所受載荷，對應(yīng)平臺環(huán)境載荷是45°方向、波浪相位角為9°。這些載荷由ANSYS平臺模型提取，坐標系與ANSYS的總體坐標系一致。

圖11 施加在模型上的約束和載荷

表5 引起卡箍最大螺栓撬力的環(huán)境載荷對應(yīng)各參考點內(nèi)力(工況1)

表6 引起卡箍最大滑動應(yīng)力的環(huán)境載荷對應(yīng)各參考點內(nèi)力(工況2)

另外在卡箍法蘭板添加螺栓作用的印記面，將螺栓預(yù)緊力134.05 kN等效施加在印記面上，并在灌漿環(huán)外表面與卡箍鞍板的內(nèi)表面施加膨脹壓力1.55 MPa。施加在模型上的約束和載荷的情況如圖11所示。

3.3 灌漿的強度校核準則

卡箍灌漿的強度校核一般使用Tresca應(yīng)力判斷，采用DNV在風(fēng)機灌漿連接段的強度校核準則[13]：

(5)

式中：fs是最大剪應(yīng)力即最大Tresca應(yīng)力；fcck是圍壓下圓柱體抗壓強度；材料系數(shù)γm=3.0。根據(jù)2.3節(jié)測試的10%膨脹劑摻量灌漿在1.5～1.6 MPa的平均fcck為62.7 MPa，則fcck/3=20.9 MPa。

3.4 不同本構(gòu)關(guān)系模型下的灌漿環(huán)強度校核分析

3.4.1 引起卡箍最大螺栓撬力的工況1的灌漿強度校核

由圖12可以看出，三種本構(gòu)關(guān)系模型的灌漿環(huán)最大Tresca應(yīng)力均小于20.9 MPa，符合要求。最大應(yīng)力位置都在與受損撐桿接觸灌漿環(huán)的端部。

圖12 工況1卡箍灌漿環(huán)的Tresca應(yīng)力

3.4.2 引起卡箍最大滑動應(yīng)力的工況2的灌漿強度校核

同樣由圖13可以看出，三種本構(gòu)關(guān)系模型的灌漿環(huán)最大Tresca應(yīng)力均小于20.9 MPa，符合要求。最大應(yīng)力位置也都在與受損撐桿接觸灌漿環(huán)的端部。

圖13 工況2卡箍灌漿環(huán)的Tresca應(yīng)力

3.4.3 增加工況1載荷下各本構(gòu)模型的灌漿強度校核分析

為了分析不同本構(gòu)模型對灌漿環(huán)最大Tresca應(yīng)力的影響，把工況1參考點2(受損桿件端)的載荷按倍數(shù)進行逐步增大，得到的最大Tresca應(yīng)力如表7所示。

表7 增加工況1參考點2的載荷下灌漿環(huán)最大Tresca應(yīng)力

可以發(fā)現(xiàn)在5倍載荷下線彈性和莫爾庫侖本構(gòu)模型下的灌漿環(huán)最大Tresca應(yīng)力都超過了20.9 MPa，不符合強度要求，而塑性損傷本構(gòu)模型下的最大Tresca應(yīng)力17.71 MPa，符合強度要求。圖14是5倍載荷下灌漿環(huán)的最大Tresca應(yīng)力位置，仍處于與受損撐桿接觸灌漿環(huán)的端部，位置基本相同。

圖14 工況1參考點2加5倍載荷下灌漿環(huán)的Tresca應(yīng)力

把表7中的參考點2載荷倍數(shù)與灌漿最大Tresca應(yīng)力進行做圖分析(圖15)，發(fā)現(xiàn)線彈性模型下的最大Tresca應(yīng)力與載荷倍數(shù)成正比且隨著載荷增加它增加最快，莫爾庫侖模型下的次之，塑性損傷模型下的最大Tresca應(yīng)力增加最慢?？梢耘袛喑鰧τ诠酀{的強度校核，線彈性模型最保守，3倍載荷之后最大Tresca應(yīng)力就超過20.9 MPa，而塑性損傷模型最不保守，5倍載荷下最大Tresca應(yīng)力仍小于20.9 MPa。

圖15 載荷增加下卡箍灌漿環(huán)的最大Tresca應(yīng)力

考慮整個灌漿卡箍是一個組合結(jié)構(gòu)，灌漿即使有一些塑性變形也會得到周圍鋼結(jié)構(gòu)的強度補償，只要灌漿與管件界面不產(chǎn)生滑移、與鞍板界面不分離，而這兩點分別由防滑和防撬校核保證。由于線彈性本構(gòu)高估了灌漿中的最大Tresca應(yīng)力，而塑性損傷模型可以精細地反映灌漿的壓碎和開裂等響應(yīng)行為，因此在卡箍強度校核的最終設(shè)計環(huán)節(jié)建議采用塑性損傷模型作為膨脹式灌漿卡箍灌漿的本構(gòu)模型。

另外從結(jié)構(gòu)非線性的計算量考慮，莫爾庫侖模型可將線彈性模型預(yù)測和塑性損傷模型預(yù)測之間的差異減半，計算量也少于塑性損傷模型，設(shè)計初期灌漿用莫爾庫侖模型也是一種選擇。本算例中線彈性、莫爾庫侖、塑性損傷三種模型的求解計算時間之比約為1∶1.6∶2.3，除了本構(gòu)模型的非線性外，主要是各部件之間的接觸非線性。

4 結(jié) 語

針對膨脹式自應(yīng)力灌漿卡箍的強度校核分析構(gòu)筑了三種灌漿本構(gòu)關(guān)系模型，測試了相關(guān)灌漿材料力學(xué)參數(shù)，使用ABAQUS軟件通過一個K型節(jié)點灌漿卡箍的算例，比較了三種灌漿本構(gòu)模型的分析結(jié)果和保守程度。主要研究結(jié)論如下：

1)測試了兩種灌漿型號的黏聚力和內(nèi)摩擦角，確立了莫爾庫侖模型主要參數(shù)；測試了兩種灌漿型號的抗拉強度和壓縮段應(yīng)力應(yīng)變曲線，確立了塑性損傷模型主要參數(shù)。

2)對于卡箍的灌漿強度校核，線彈性模型最保守，容易校核不合格，而塑性損傷模型最不保守。

3)考慮到灌漿卡箍是一個組合結(jié)構(gòu)，建議采用能精細地反映灌漿的壓碎和開裂行為的塑性損傷模型作為膨脹式灌漿卡箍灌漿的本構(gòu)關(guān)系模型。

研究成果可以為自應(yīng)力灌漿卡箍的有限元強度校核方法提供有益的借鑒和指導(dǎo)。