吳夢寧, 余建星, 孫震洲, 段晶輝

(1. 天津大學 水利工程仿真與安全國家重點試驗室，天津 300354; 2. 高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心，上海 200240)

雙整體式止屈器結(jié)構(gòu)性能研究及優(yōu)化設(shè)計

吳夢寧1, 2, 余建星1, 2, 孫震洲1, 2, 段晶輝1, 2

(1. 天津大學 水利工程仿真與安全國家重點試驗室，天津 300354; 2. 高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心，上海 200240)

建立了雙整體式止屈器結(jié)構(gòu)的三維有限元模型，模擬其在外部靜水壓力作用下的屈曲穿越過程，總結(jié)不同穿越破壞模式的產(chǎn)生條件及機理。對影響穿越破壞壓力的主要參數(shù)進行敏感性分析，結(jié)果表明止屈器布置間距、兩個止屈器的幾何尺寸均會對整體止屈性能產(chǎn)生影響。在擬合的雙整體式止屈器效率公式的基礎(chǔ)上，依據(jù)非線性有約束優(yōu)化理論，分別采用經(jīng)典算法和智能遺傳算法確定模型最優(yōu)解，提出了一種兼顧結(jié)構(gòu)重量與止屈性能的最優(yōu)雙整體式止屈器結(jié)構(gòu)型式；并經(jīng)計算分析，驗證了該方法的可靠性和適應(yīng)性。

雙整體式止屈器；穿越模式；擬合公式；結(jié)構(gòu)設(shè)計

Abstract: A finite element(FE) model was developed to simulate the crossover of the double-integral buckle arrestors(DIBA) under the ambient externa hydrostatic pressure. In this study, the condition and mechanism of different crossover modes during the process of arresting were summarized. Then, the impact based on series of structure parameters on the DIBA’s crossover pressure was studied. The results indicated that the arresting performance of the model could be affected by many factors, such as, the space between arrestors and the geometrical dimensions of each arrestor. After analyzing the result of sensitivity, a fitting formula on DIBA’s arresting efficiency was put forward in the paper. At last, the structure optimization model of DIBA was established, both arithmetic and genetic algorithms are introduced to solve the problem, which is based on the constrained nonlinear programming theory. The result shows that the optimal double-integral buckle arrestor is an advanced buckling arrest measure, with excellent arresting performance and lower structural weight. And the measure has been verified to have good performance in reliability and adaptability through calculation and analyse.

Keywords: double-integral buckle arrestors; crossover modes; fitting formulas; structural design

近年來，隨著海洋油氣資源勘探開發(fā)的重點開始轉(zhuǎn)向深水及超深水區(qū)域，巨大外部靜水壓力及惡劣復雜的作業(yè)環(huán)境所帶來的海底管道屈曲失效問題愈加嚴重。相較于長距離鋪設(shè)厚壁管道在結(jié)構(gòu)重量、鋪設(shè)難度及維修成本等方面的弊端，沿海底管道長度方向間隔布置止屈器的方法優(yōu)勢明顯。從20世紀80年代開始，各類止屈器就已經(jīng)成為海底管道止屈措施研究的熱點，但目前還處于理論探討和室內(nèi)試驗階段，沒有投入實際應(yīng)用。研究證明，整體式止屈器以其優(yōu)異的止屈性能，是各類止屈器中最適用于(超)深水海底管道鋪設(shè)中的止屈裝置[1-2]。整體式止屈器是焊接在管道指定位置處的厚壁圓環(huán)，通過局部加強管道剛度，從而防止屈曲沿管道軸向傳播的裝置。由于海底管道的鋪設(shè)距離逐年增加，若采用整體式止屈器抑制屈曲傳播的方法，勢必需沿管道布置多個。對于單個整體式止屈器，其止屈性能受多重因素影響，主要包括管道徑厚比[3]、止屈器長度及厚度[4-6]、材料特性[3， 7-8]、外部環(huán)境[9-11]等。相比之下，多個整體式止屈器共同作用的情況就更加復雜，需在此基礎(chǔ)上考慮不同布置間距對整體止屈性能的影響效果、不同尺寸止屈器間搭配所產(chǎn)生的相互作用等問題。近幾年來，對整體式止屈器的研究，總體上局限于單個止屈器的設(shè)計標準及形式方面[12-14]，但僅憑借單個止屈器的試驗研究結(jié)果，不足以有效指導其在實際工程中的應(yīng)用。美國奧斯汀大學曾沿管長布置兩個整體式止屈器進行試驗研究[4-5]，但其目的在于降低試驗次數(shù)，提高試驗效率，因此在試驗中將兩個止屈器遠距離布置，使二者間不足以相互影響。李旭等[15]在對止屈器設(shè)計方法的研究中，僅基于管材采辦費用確定止屈器間距，并未考慮止屈器間的力學性能?？梢姡瑢Χ鄠€止屈器聯(lián)合作用下止屈效果的研究，目前還處于空白階段，對其研究具有極其重要的工程意義。另外，在結(jié)構(gòu)分析方法方面，相較于傳統(tǒng)設(shè)計方法需多次迭代、重復分析，價格昂貴等缺點，遺傳算法、神經(jīng)網(wǎng)絡(luò)等智能優(yōu)化法于20世紀80年代相繼提出，該類方法最優(yōu)解的獲得不需要目標函數(shù)和約束條件的導數(shù)信息，為復雜工程的優(yōu)化設(shè)計提供了新的思路和方法，近年來已成功應(yīng)用于實際工程優(yōu)化設(shè)計中。

本文以雙整體式止屈器為研究對象，通過數(shù)值仿真技術(shù)完成對不同結(jié)構(gòu)型式的雙止屈器在外部靜水壓力作用下的屈曲穿越過程模擬，研究止屈性能隨各設(shè)計參數(shù)的變化趨勢。在此基礎(chǔ)上，建立雙止屈器結(jié)構(gòu)的優(yōu)化設(shè)計模型，通過結(jié)構(gòu)參數(shù)優(yōu)化設(shè)計，得到滿足止屈性能、制造成本等要求的最優(yōu)結(jié)構(gòu)。

1 屈曲穿越雙整體式止屈器有限元模型

1.1模型建立

首先建立具有0.2%端部橢圓度缺陷的管道帶雙止屈器的整體幾何模型。根據(jù)外部靜水壓力作用下模型壓潰變形呈對稱的特點，將x=0、y=0、z=0三個平面均設(shè)置為對稱面，僅建立1/4模型以進行簡化。以止屈器為分界，將模型分為上游管段LU(包括屈曲誘發(fā)段LUO、純圓段LU1與LUc)、第一止屈器段(LAR1與兩側(cè)過渡段LT1)、中間管段LG、第二止屈器段(LAR2與兩側(cè)過渡段LT2)及下游管段LD(包括LD1和LDc)五部分。具體幾何模型如圖1所示。

圖1 雙整體式止屈器幾何模型Fig. 1 Model of double-integral buckle arrestors

圖2 雙止屈器附近的網(wǎng)格劃分Fig. 2 Meshes near double-integral buckle arrestors

管道及止屈器模型均采用8節(jié)點六面體線性非協(xié)調(diào)模式單元C3D8I。為模擬后屈曲階段管壁的接觸效應(yīng)，引入兩個分別位于XZ平面與YZ平面的解析剛性面作為主面，各自與作為從面的管道內(nèi)壁面構(gòu)成接觸對，并設(shè)置剛體參考點施加六個自由度全約束。在模型外部區(qū)域建立流體加載艙，采用流體質(zhì)量加載方法，通過手工修改原始模型的INP代碼施加外部靜水壓力。網(wǎng)格劃分上，沿管道壁厚方向劃分2個單元、圓周方向劃分24個單元，沿軸線方向分成不同長度的區(qū)段，根據(jù)不同區(qū)段的長度及所關(guān)心的重要性決定劃分單元的個數(shù)，如在雙止屈器兩側(cè)設(shè)置網(wǎng)格加密區(qū)以保證模擬效果，如圖2所示。

1.2模擬結(jié)果

以外徑51 mm、壁厚2 mm的管件，配兩個長度40.8 mm、壁厚4.4 mm的整體式止屈器所組成的雙止屈器模型進行止屈建模計算，兩個止屈器的間距取為4.0倍管道外徑，材質(zhì)均為SS304鋼材。通過ABAQUS軟件進行雙止屈器的屈曲穿越模擬，輸出每一個載荷子步的靜水流體壓力和靜水流體腔體體積，通過換算得到模型所受外壓和流體腔內(nèi)體積的壓縮量，繪出準靜態(tài)加載壓力曲線，如圖3所示。

圖3 準靜態(tài)加載壓力曲線Fig. 3 Quasi-static pressure curve of double-integral buckle arrestors

圖4 模型屈曲穿越變形過程云圖Fig. 4 Mises stress nephogram of buckling and crossover

由圖可知，雙止屈器屈曲穿越過程與單止屈器類似。加載過程中壓力上升至①處(約7.2 MPa)時管道局部缺陷處開始壓潰，壓力迅速下降，直至②處管道內(nèi)壁首次相互接觸(約3.3 MPa)時壓力回升，此時壓潰變形停止，屈曲開始向下游傳播。隨后由于第一個止屈器對管道的強化作用，壓力繼續(xù)上升直到屈曲傳播至該止屈器。該過程的壓力峰值③即第一個止屈器的穿越破壞壓力為13.54 MPa。當屈曲穿越第一個止屈器后，管道外壓由于中間管道的壓潰會迅速下降，隨屈曲傳播逐漸接近第二個止屈器又逐步上升，達到的第二個壓力峰值④即第二個止屈器的穿越破壞壓力，為13.02 MPa。模型屈曲穿越過程的應(yīng)力云圖如圖4所示。

1.3穿越破壞模式

雙整體式止屈器結(jié)構(gòu)在屈曲穿越過程中，會呈現(xiàn)出較單止屈器更為多樣的穿越破壞模式。為避免產(chǎn)生歧義，以“未穿越”、“0°穿越”、“90°穿越”對其中每一止屈器的穿越模式進行描述。具體來說，“未穿越模式”主要針對第一個止屈器連同中間管段均不會壓潰而保持完整性的情況，此時稱第一個止屈器發(fā)生“未穿越模式”；“0°穿越模式”表示止屈器下游管段壓潰方向與上游屈曲傳播方向相同；“90°穿越”表示止屈器下游管段壓潰方向與上游屈曲傳播方向呈90°。

近間距雙止屈器結(jié)構(gòu)在外部靜水壓力下，會出現(xiàn)第一個止屈器不穿越、中間管段不壓潰的現(xiàn)象，這種穿越破壞模式統(tǒng)稱為單穿越模式。根據(jù)兩止屈器相對剛度大小及間距的不同，細分為“未穿越-90°穿越”模式和“未穿越-0°穿越”模式兩種。當止屈器布置間距較近時，兩止屈器連同中間管段可看作整體、共同抵抗屈曲傳播，下游管段由于反向橢圓度的作用發(fā)生90°穿越，如圖5(a)所示。當止屈器布置間距擴大到需考慮兩止屈器對中間管段的影響效果、且第一個止屈器剛度較大時，屈曲波形仍可被第一個止屈器有效阻止，中間管段首先產(chǎn)生由上游屈曲帶來的反向橢圓度，隨后使下游管段產(chǎn)生“反向-反向”橢圓度，下游管段的變形方向由該橢圓度方向決定，產(chǎn)生了與原屈曲方向一致的壓潰，同時雙止屈器結(jié)構(gòu)產(chǎn)生“未穿越-0°穿越”模式，如圖5(b)所示。

遠間距雙止屈器結(jié)構(gòu)在外部靜水壓力下會出現(xiàn)兩個止屈器各自發(fā)生穿越、中間管段被壓潰的現(xiàn)象，這種穿越破壞模式統(tǒng)稱為雙穿越模式。具體穿越模式與其中包含的兩個止屈器剛度均相關(guān)，表現(xiàn)為整體穿越模式與其中止屈器單個作用時的穿越模式對應(yīng)，根據(jù)不同剛度組合的雙止屈器模型分為“0°-0°穿越”模式、“0°-90°穿越”模式、“90°-0°穿越”模式和“90°-90°穿越”模式四種。將不同穿越模式的發(fā)生條件總結(jié)至表1，四種穿越模式分別如圖6(a)～6(d)所示。

圖5 單穿越模式Fig. 5 Single corssover mode

雙止屈器穿越模式單個止屈器單獨作用時穿越模式第一個止屈器第二個止屈器0°-0°平行平行0°-90°平行正交90°-0°正交平行90°-90°正交正交

圖6 雙穿越模式Fig. 6 Double-crossover mode

2 雙整體式止屈器結(jié)構(gòu)性能分析

雙整體式止屈器止屈效果的評價，以求得整體結(jié)構(gòu)的穿越破壞壓力為核心目標。主要關(guān)注兩個止屈器的布置間距，第一個止屈器的壁厚、長度，第二個止屈器的壁厚、長度等參數(shù)對結(jié)構(gòu)穿越壓力的影響效果。分析中參數(shù)D表示管道外徑，t表示管道壁厚，La1表示第一個止屈器長度，h1表示第一個止屈器厚度，La2表示第二個止屈器長度，h2表示第二個止屈器厚度，L表示兩止屈器布置間距。

2.1穿越破壞壓力對止屈器布置間距的敏感性

不同剛度止屈器所組成的雙止屈器模型，屈曲穿越機理有明顯區(qū)別，因此以止屈器厚度作為剛度評價指標，對雙止屈器模型進行分類，研究不同剛度組合雙止屈器模型中穿越壓力對間距的敏感性。

模型中選用厚度h為3.0t和2.2t的兩種止屈器，分別定義為“厚”止屈器和“薄”止屈器，長度La固定為0.8D，以兩個止屈器之間的距離L作為變化參數(shù)，得到不同模型穿越破壞壓力隨止屈器布置間距的敏感性曲線，如圖7所示。圖中紅色的水平線為相同材料及幾何尺寸的單止屈器穿越破壞壓力。“單穿越壓力”曲線對應(yīng)單穿越模式中整體結(jié)構(gòu)穿越壓力結(jié)果，“第一、第二穿越壓力”曲線分別對應(yīng)雙穿越模式中的第一個、第二個止屈器穿越破壞壓力結(jié)果。

圖7 穿越破壞壓力的間距敏感性曲線Fig. 7 Crossover pressure plotted against the space between arrestors

由敏感性曲線可知，止屈器布置間距對穿越破壞壓力的影響作用顯著：

1)隨止屈器布置間距的增加，各雙止屈器模型均表現(xiàn)出從單穿越模式到雙穿越模式的變化趨勢；

2)單穿越模式中，穿越破壞壓力隨布置間距的增加先升高后降低，穿越破壞壓力水平總體較高；雙穿越模式中，第一穿越壓力隨布置間距的增加先大幅降低、后基本穩(wěn)定在對應(yīng)的單止屈器穿越壓力水平，第二穿越壓力隨布置間距的增加先小幅提升、后基本穩(wěn)定不變，穩(wěn)定值根據(jù)第二止屈器剛度而有所區(qū)別。

綜上，近間距雙止屈器結(jié)構(gòu)適用性更高。一方面，近間距布置能保證中間管段的結(jié)構(gòu)完整性，且整體結(jié)構(gòu)的止屈效果顯著優(yōu)于單個止屈器及間距較遠雙止屈器；另一方面，中間管段越短，整個雙止屈器結(jié)構(gòu)的用料越少，在經(jīng)濟方面具有明顯的優(yōu)勢。

2.2穿越破壞壓力對止屈器厚度的敏感性

鑒于近間距雙止屈器模型兼具止屈效率與結(jié)構(gòu)重量上的優(yōu)勢，針對間距為3.0D以內(nèi)的雙止屈器模型，進行單止屈器厚度對穿越破壞壓力的敏感性計算。

圖8(a)為固定兩個止屈器長度均為0.80D，第二個止屈器厚度h2為1.8t，模型在不同布置間距下，穿越破壞壓力隨第一個止屈器厚度的變化趨勢。圖8(b)中止屈器長度及布置間距與上述相同，固定第一個止屈器厚度h1為3.0t，僅改變第二個止屈器厚度進行計算。兩圖中紅色水平線代表管道的壓潰壓力，不同顏色曲線代表不同布置間距結(jié)果。另外，為比較兩厚度對雙止屈器模型穿越破壞壓力影響效果的差異，計算固定間距(L=1.0D)下穿越破壞壓力隨單止屈器厚度的變化趨勢，如圖8(c)所示。

圖8 穿越破壞壓力的止屈器厚度敏感性曲線Fig. 8 Crossover pressure plotted against thickness of arrestor

由敏感性曲線可知：

1)止屈器間距在1.0D范圍內(nèi)，穿越破壞壓力隨每一止屈器厚度的增加而逐漸增加；當間距擴大到1.0D至3.0D范圍內(nèi)，穿越過程復雜，穿越壓力的大小要考慮反向橢圓度等影響因素，與每一止屈器厚度并無明顯關(guān)系，但各組結(jié)果相差不大，此時止屈器厚度對模型止屈性能的影響效果降低。

2)固定布置間距下，穿越破壞壓力對第一個止屈器厚度的敏感性更強，隨第一個止屈器厚度的增加顯著提升，隨第二個止屈器厚度的增加而變化不明顯。

綜上，第一個止屈器厚度是整體結(jié)構(gòu)穿越破壞壓力的強敏感性因素，第二個止屈器厚度是弱敏感性因素。實際設(shè)計中，可適當減薄第二個止屈器，以降低整體結(jié)構(gòu)的重量，減少安裝難度及成本。

2.3穿越破壞壓力對止屈器長度的敏感性

固定第一個止屈器厚度為1.8t，第二個止屈器厚度為3.0t，止屈器原始長度均為0.80D，分別改變單個止屈器的長度，研究近間距雙止屈器模型在不同布置間距下，穿越破壞壓力隨單止屈器長度的變化趨勢，結(jié)果如圖9(a)、9(b)所示。同時，建立固定間距L為1.0D的雙“薄”止屈器及雙“厚”止屈器模型，僅改變單一止屈器長度，對比分析兩止屈器長度對穿越破壞壓力的影響效果，結(jié)果如圖9(c)、9(d)所示。圖9(c)中藍色曲線在0.4D、0.5D的不連續(xù)點，代表模型在該幾何尺寸下發(fā)生了雙穿越模式。

由敏感性曲線可知：

與止屈器厚度的影響效果類似，間距在1.0D范圍內(nèi)穿越破壞壓力隨單止屈器長度的增加而近似線性增加，間距在1.0D至 3.0D范圍內(nèi)長度的影響效果降低。

固定間距下，第一個止屈器的長度對雙止屈器模型的止屈性能影響更為明顯，但影響效果并不具有很好的規(guī)律性，需針對具體情況分析；第二個止屈器長度對雙止屈器模型止屈性能的影響總體不大。

綜上，與第二個止屈器長度相比，雙止屈器結(jié)構(gòu)穿越破壞壓力對第一個止屈器長度的敏感性更強。在第一個止屈器長度較短情況下，厚度較薄時會造成整體剛度不夠而使兩個止屈器均發(fā)生穿越；厚度較厚時可能使穿越模式發(fā)生變化同時穿越壓力發(fā)生波動。因此，在設(shè)計中可適當加長第一個止屈器，以保證雙止屈器結(jié)構(gòu)的止屈性能穩(wěn)定性。

圖9 穿越破壞壓力的止屈器長度敏感性曲線Fig. 9 Crossover pressure plotted against length of arrestor

3 雙整體式止屈器結(jié)構(gòu)優(yōu)化設(shè)計

在雙整體式止屈器的設(shè)計過程中，理論上存在一個最佳尺寸。結(jié)構(gòu)尺寸過大，會造成材料的浪費；但尺寸過小又起不到抑制屈曲傳播的目的，最佳結(jié)構(gòu)需在止屈效率與結(jié)構(gòu)重量間達到平衡。

3.1數(shù)學模型

3.1.1 設(shè)計變量

以無量綱參數(shù)止屈器間距(L/D)、第一個止屈器長度(La1/D)、第一個止屈器厚度(h1/t)、第二個止屈器長度(La2/D)、第二個止屈器厚度(h2/t)作為設(shè)計變量，以式(1)所示的向量X代表雙整體式止屈器的設(shè)計方案。

3.1.2 目標函數(shù)

由于設(shè)計變量僅有尺寸參數(shù)，運用結(jié)構(gòu)最小重量的設(shè)計方法，將雙整體式止屈器結(jié)構(gòu)重量最低作為目標函數(shù)，如式(2)所示。式中一、二項分別為兩個止屈器的重量，第三項為中間管段的重量。

3.1.3 約束條件

由于止屈器最重要的功能是抑制屈曲形式的傳播，以止屈效率約束其結(jié)構(gòu)設(shè)計。止屈效率η[16]計算公式如式(3)所示。在單止屈器的設(shè)計中，一般認為止屈效率達到78%即為高效率止屈器[17]。為體現(xiàn)雙整體式止屈器結(jié)構(gòu)在止屈性能上的優(yōu)勢，結(jié)構(gòu)設(shè)計中規(guī)定止屈效率不低于1.0，即如式(4)：

鑒于雙整體式止屈器的止屈效率沒有現(xiàn)有理論及規(guī)范參考，利用有限元計算結(jié)果進行止屈效率公式的擬合。之前研究成果證明，單止屈器穿越壓力與管道直徑D、壁厚t、屈曲傳播壓力Pp均相關(guān)。加上雙止屈器模型穿越壓力大小依賴于止屈器間距及兩個止屈器尺寸，因此將雙穿越破壞壓力按式(5)表示：

根據(jù)白金漢(Buckingham)π定律，將D和t作為基本量綱，對式(5)進行無量綱化，得：

由于海底管道的PX≥PP，因此A0=1。為降低擬合的難度，將n>1的各情況省略后，得到近似關(guān)系為：

再通過將止屈效率公式(3)帶入式(8)，得到止屈效率的近似公式(9)：

式中：αi(i=1～5)為待定相關(guān)系數(shù)。

采用最小二乘曲線擬合的方法，得到有限元結(jié)果與擬合結(jié)果關(guān)系散點圖，如圖10所示。圖中隨止屈器模型的幾何尺寸增加，止屈效率結(jié)果的離散性逐漸顯著，特別是在止屈效率大于1.1范圍內(nèi)，止屈效率與擬合公式的關(guān)系不再呈現(xiàn)出遞增的態(tài)勢?；陔p止屈器結(jié)構(gòu)優(yōu)化設(shè)計的目的是得到止屈效率在1.0附近、能使整體結(jié)構(gòu)重量最輕的結(jié)構(gòu)型式，因此為提高擬合曲線的精度，在數(shù)據(jù)處理階段去掉止屈效率超過1.1的各計算結(jié)果，重新擬合止屈效率公式，結(jié)果如式(10)所示。

圖10 雙止屈器結(jié)構(gòu)止屈效率的擬合Fig. 10 Fitting result of arresting efficiency

經(jīng)計算，該公式的相關(guān)系數(shù)r=0.904 4，可決系數(shù)R2=0.891 1，擬合優(yōu)度理想。另外，隨機抽取10組擬合公式計算結(jié)果同有限元結(jié)果對比，誤差均在5%以內(nèi)，計算結(jié)果較為可靠。將式(10)帶入至約束條件(4)中，作為雙止屈器結(jié)構(gòu)優(yōu)化設(shè)計的止屈效率取值依據(jù)。

3.1.4 最終模型

由于雙止屈器結(jié)構(gòu)中設(shè)計參數(shù)眾多，采用控制單一變量的方法逐一確定設(shè)計變量范圍困難較大。考慮采用遺傳算法與BP神經(jīng)網(wǎng)絡(luò)結(jié)合的方法，不預設(shè)止屈效率的顯性表達式，以雙止屈器結(jié)構(gòu)重量作為適應(yīng)度值，利用神經(jīng)網(wǎng)絡(luò)的強大非線性擬合能力和遺傳算法的非線性尋優(yōu)能力，尋找雙止屈器模型的最優(yōu)結(jié)構(gòu)尺寸。在BP神經(jīng)網(wǎng)絡(luò)中，雖然通過遺傳算法可選取最優(yōu)初始權(quán)值及閥值，但訓練及測試數(shù)據(jù)是在已有數(shù)據(jù)中隨機選取的。該方法會導致每次的最優(yōu)結(jié)果有所區(qū)別，因而不能將某次的優(yōu)化結(jié)果作為最優(yōu)雙止屈器結(jié)構(gòu)尺寸。鑒于通過BP神經(jīng)網(wǎng)絡(luò)方法得到的結(jié)果理論上均為所選變量域內(nèi)的最優(yōu)解，因此可通過統(tǒng)計BP神經(jīng)網(wǎng)絡(luò)計算得到的一系列最優(yōu)結(jié)果的范圍，作為設(shè)計變量的上下限約束值。

依據(jù)該思路，最終確定設(shè)計變量的取值范圍如表2所示。

表2 設(shè)計變量范圍Tab. 2 Scope of design variables

通過前述分析，最終確定止屈器優(yōu)化設(shè)計的數(shù)學模型如式(11)所示：

3.2優(yōu)化結(jié)果

依據(jù)非線性有約束優(yōu)化理論，采用經(jīng)典算法和遺傳算法分別對優(yōu)化模型(11)進行求解，結(jié)果如表3所示。

表3 最優(yōu)結(jié)果對比表Tab. 3 Comparison of optimal solution

經(jīng)比較，兩種算法結(jié)果基本一致。最優(yōu)結(jié)構(gòu)中第一止屈器厚度較厚、兩止屈器間距較近，符合敏感性研究結(jié)果。

3.3最優(yōu)型式的進一步分析

利用數(shù)值仿真模擬最優(yōu)結(jié)構(gòu)的屈曲穿越過程。在Abaqus中根據(jù)優(yōu)化得到的參數(shù)建立最優(yōu)雙整體式止屈器結(jié)構(gòu)的有限元模型，如圖11所示。

圖11 最優(yōu)雙止屈器結(jié)構(gòu)模型Fig. 11 Optimal model of double-integral buckle arrestors

通過Abaqus模擬最優(yōu)雙止屈器在靜水壓力中的穿越過程，得到的止屈效率為99.15%，對比擬合公式計算值100%，二者間差異僅為0.86%，認為通過優(yōu)化設(shè)計得到的最優(yōu)雙止屈器結(jié)構(gòu)型式滿足要求。

考慮到雙止屈器模型可能出現(xiàn)由于結(jié)構(gòu)重量大、制造費用高而影響其實際應(yīng)用的問題，以經(jīng)濟指標為判斷依據(jù)，就制造費用對相同止屈效率的最優(yōu)雙止屈器結(jié)構(gòu)型式與單止屈器型式進行對比。

以有限元結(jié)果作為最優(yōu)雙整體式止屈器的止屈效率，即η優(yōu)=99.15%，將單止屈器厚度取為雙止屈器中第一個止屈器的厚度，即h1=3.0t。經(jīng)不斷試算得到當單止屈器長度為1.00D時，其止屈效率與最優(yōu)雙止屈器結(jié)構(gòu)止屈效率相當。分別計算單止屈器和雙止屈器結(jié)構(gòu)的重量，比較在止屈性能相當?shù)那闆r下，二者由重量引起的實際鋪設(shè)難易。h為單止屈器厚度，La為單止屈器長度，鋼材密度ρ為7 850 kg/m3。

1)單止屈器：

2)雙止屈器：

經(jīng)計算，發(fā)現(xiàn)對于同一預期止屈效率，最優(yōu)雙整體式止屈器結(jié)構(gòu)的重量低于單整體式止屈器，具有明顯的重量優(yōu)勢?？梢?，在止屈器深水領(lǐng)域適用性問題的研究上，除了制造成本之外，設(shè)備重量具有更為重要的意義。盡管在制造和鋪設(shè)過程中，雙止屈器方案肯定更為耗時與昂貴，但其重量較低的優(yōu)勢使管道能夠更易鋪設(shè)與應(yīng)用。

4 結(jié) 語

本文對雙整體式止屈器進行了一系列止屈性能的研究，得到如下結(jié)論：

1)雙整體式止屈器根據(jù)止屈器剛度及布置間距差異，在屈曲穿越時會產(chǎn)生六種不同的穿越破壞模式。

2)止屈器布置間距對整體結(jié)構(gòu)的穿越破壞壓力影響顯著。近間距雙整體式止屈器結(jié)構(gòu)具有明顯的止屈性能優(yōu)勢。

3)相較于第二個止屈器，第一個止屈器長度及厚度對整體結(jié)構(gòu)的穿越破壞壓力影響顯著。在雙止屈器結(jié)構(gòu)的實際設(shè)計中，可適當增加第一個止屈器尺寸、減小第二個止屈器尺寸，以在保證止屈性能的前提下減輕結(jié)構(gòu)重量，最大化地降低安裝難度。

4)采用最小二乘曲線方法擬合得到雙整體式止屈器止屈效率擬合公式，其變量均為雙止屈器特征參數(shù)的無量綱化形式，適用性較廣，可對除本文外的雙止屈器模型止屈性能進行有效預測。

5)以整體結(jié)構(gòu)重量最低為目標，止屈效率盡可能接近1.0為約束條件，建立最優(yōu)雙整體式止屈器結(jié)構(gòu)優(yōu)化設(shè)計的數(shù)學模型，得到最優(yōu)設(shè)計參數(shù)。通過對比，認為最優(yōu)雙止屈器結(jié)構(gòu)無論在止屈性能還是經(jīng)濟成本上，均具有明顯優(yōu)勢。

[1] NETTO T A, ESTENFEN S F. Buckle arrestors for deepwater pipelines[J]. Marine Structures, 1996, 9(95): 873-883.

[2] 汪紅艷, 陳曉賢, 林金保, 等. 深水海底管道止屈器設(shè)計研究[J].中國造船, 2011, 52(z1): 123-131. (WANG Hongyan, CHEN Xiaoxian, LIN Jinbao, et al. The design of internal buckle arrestors on deepwater pipelines.[J]. Shipbuilding of China, 2011, 52(z1): 123-131.(in Chinese))

[3] LEE L H, KYRIAKIDES S, NETTO T A. Integral buckle arrestors for offshore pipelines: Enhanced design criteria[J]. International Journal of Mechanical Sciences, 2008, 50(6): 1 058-1 064.

[4] PARK T D, KYRIAKIDES S. On the performance of integral buckle arrestors for offshore pipelines[J]. International Journal of Mechanical Sciences, 1997, 39(6): 643-669.

[5] KYRIAKIDES S, PARK T D, NETTO T A. On the design of integral buckle arrestors for offshore pipelines[J]. Applied Ocean Research, 1998, 20(1): 95-104.

[6] 翟宇軒, 余建星, 孫震洲, 等. 海底管道整體式止屈器的止屈性能研究[J]. 天津理工大學學報, 2014(5): 46-51. (ZHAI Yuxuan, YU Jianxing, SUN Zhenzhou, et al. The arresting perfbrmance of integral buckle arrestors for offshore pipelines[J]. Journal of Tianjin Unversity of Technology, 2014(5): 46-51.(in Chinese))

[7] MANTOVANO L O, AMENTA P. Finite element modeling and experimental validation of buckle arrestors for deepwater pipelines.[J]. Mecnica Computacional, 2006, 6(1): 687-706.

[8] TOSCANO R G, MANTOVANO L O, AMENT P M, et al. Collapse arrestors for deepwater pipelines. Cross-over mechanisms[J]. Computers & Structures, 2008, 86(78): 728-743.

[9] NETTO T A, KYRIAKIDES S. Dynamic performance of integral buckle arrestors for offshore pipelines. Part I: Experiments[J]. International Journal of Mechanical Sciences, 2000, 42(7): 1 405-1 423.

[10] NETTO T A, KYRIAKIDES S. Dynamic performance of integral buckle arrestors for offshore pipelines. Part II: Analysis[J]. International Journal of Mechanical Sciences, 2000, 42(7): 1 425-1 452.

[11] 余建星, 林曉龍, 楊源, 等. 深水海底管道整體式止屈器屈曲動態(tài)傳播及止屈研究[J]. 海洋技術(shù)學報, 2013, 32(2): 60-65. (YU Jianxing, LIN Xiaolong, YANG Yuan, et al. Dynamic performance of propagation buckle and integral buckle arrestors in offshore pipelines[J]. Journal of Ocean Technology, 2013, 32(2): 60-65.(in Chinese))

[12] LEE L H, KYRIAKIDES S, NETTO T A. Integral buckle arrestors for offshore pipelines: Enhanced design criteria[J]. International Journal of Mechanical Sciences, 2008, 50(6):1 058-1 064.

[13] WANG H, SUN J, JUKES P. FEA of a laminate internal buckle arrestor for deep water pipe-in-pipe flowlines[C]// ASME 2009 28th International Conference on Ocean, Offshore and Arctic Engineering. America: American Society of Mechanical Engineers, 2009: 447-452.

[14] MANSOUR G N, TASSOULAS J L. Crossover of integral-ring buckle arrestor: Computational results[J]. Journal of Engineering Mechanics, 2014, 123(4): 359-366.

[15] 李旭, 李慶, 張英. 深水海管止屈器設(shè)計方法研究[J]. 海洋工程裝備與技術(shù), 2014, 1(3): 227-229. (LI Xu, LI Qing, ZHANG Ying. Deepwater pipeline arrestor design study[J]. Ocean Engineering Equipment and Technology, 2014, 1(3): 227-229.(in Chinese))

[16] KYRIAKIDES S, BABCOCK C D. On the “Slip-On” buckle arrestor for offshore pipelines[J]. Journal of Pressure Vessel Technology, 1980, 102(2):188-193.

[17] 羅超, 王琮, 任翠青. 海底管道整體式止屈器結(jié)構(gòu)設(shè)計及數(shù)值模擬[C]// 第十五屆中國海洋(岸)工程學術(shù)討論會論文集(上). 太原, 中國: 海洋出版社. 2011: 363-367. (LUO Chao, WANG Zong, REN Cuiqing. Structure design and numerical simulation of intergral buckle arrestors on deepwater pipelines[C]// Proceedings of the Ffifteenth China Ocean (Coastal) Engineering Symposium (Volume 1). Taiyuan, China: China Ocean Press, 2011:363-367.(in Chinese))

Arresting performance analysis and optimum design for the structure of double-integral buckle arrestors

WU Mengning1, 2，YU Jianxing1, 2, SUN Zhenzhou1, 2, DUAN Jinghui1, 2

(1. State Key Laboratory of Hydraulic Engineering Simulation and Safety ,Tianjin University, Tianjin 300072, China; 2. Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration, Shanghai 200240, China)

1005-9865(2017)02-0056-11

P751

A

10.16483/j.issn.1005-9865.2017.02.009

2016-03-04

國家自然科學基金資助項目(51239008)；國家重點基礎(chǔ)研究發(fā)展計劃資助項目(973)(2014CB046805)；國家自然科學基金創(chuàng)新研究群體科學基金資助項目(51321065)

吳夢寧(1991-)，女，天津人，碩士研究生，從事海洋工程專業(yè)研究。E-mail：Lv_yico@163.com