全海深沉積物保壓容器結(jié)構(gòu)多目標優(yōu)化設計

劉廣平，金永平，彭佑多，萬步炎

（湖南科技大學 海洋礦產(chǎn)資源探采裝備與安全技術國家地方聯(lián)合工程實驗室，湘潭 411201）

0 引言

海洋是人類巨大的共同資源寶庫，蘊藏著豐富的礦產(chǎn)、生物、油氣以及其他多種資源。隨著世界經(jīng)濟的快速發(fā)展，世界各國開發(fā)利用海洋資源的需求與競爭將愈演愈烈[1,2]。深淵海底沉積物中蘊藏著大量的微生物等生命群落，這些微生物群落是人類認識和研究深淵生命演化及深淵環(huán)境變化的重要途徑[3,4]。為了獲取深淵海底沉積物中的嗜壓微生物，必須開展全海深沉積物氣密取樣器的研制工作，而全海深沉積物保壓容器作為全海深沉積物氣密取樣器中轉(zhuǎn)移裝置的主要部件，對轉(zhuǎn)移裝置的保壓性能有著重要的影響。鑒于全海深沉積物轉(zhuǎn)移裝置特殊的作業(yè)環(huán)境，在全海深沉積物保壓容器的設計過程中，應盡可能地減少保壓容器整體重量；然而，保壓容器又屬于高壓容器，在高壓容器的設計過程中，為了安全需要，又應盡可能地降低保壓容器所受的最大應力。針對這一矛盾問題，在進行保壓容器保壓性能優(yōu)化設計中，應同時考慮保壓容器重量和最大應力這兩個不同的目標。因此，該優(yōu)化又屬于一個多目標優(yōu)化問題。在多目標優(yōu)化過程中，各個目標之間既聯(lián)系又相互沖突，即過度減少保壓容器的重量可能會導致保壓容器最大應力值增加，反之亦然。要使多個目標同時達到最優(yōu)值幾乎是不可能的，因此必須對各個目標的組織和協(xié)調(diào)進行深入的研究。

本文以全海深沉積物保壓容器為研究對象，應用有限元軟件ANSYS中的Workbench模塊對全海深保壓容器進行有限元分析研究。采用實驗設計方法，研究全海深沉積物保壓容器結(jié)構(gòu)設計參數(shù)與保壓容器重量及最大應力之間的關系，并以保壓容器重量和最大應力為優(yōu)化目標，對保壓容器的幾何結(jié)構(gòu)設計參數(shù)進行多目標優(yōu)化設計。

1 保壓容器結(jié)構(gòu)參數(shù)與數(shù)值模擬方法

1.1 計算模型

為研究保壓容器結(jié)構(gòu)對保壓容器保壓性能的影響，同時也便于分析比較，本文以國家重點研發(fā)計劃“全海深沉積物氣密取樣器研制”課題所需的保壓容器為研究對象，根據(jù)課題技術指標，單個保壓容器所需的儲存容積為V=150ml。在本文的初步設計中，取保壓容器的高度（筒長）為h=212mm，保壓容器的內(nèi)徑為D1=30mm，外徑為D2=52mm。根據(jù)設計要求，選取保壓容器的筒體材料為TC4鈦合金，TC4鈦合金的力學性能為：抗拉強度σb≥950MPa，為確保保壓容器的使用安全性，選取安全系數(shù)為nb=2.7，則許用應力為保壓容器端蓋的連接方式采用螺栓連接，螺栓對稱分布在保壓容器的端蓋上。保壓容器的結(jié)構(gòu)圖和三維圖如圖1所示，保壓容器的初步設計參數(shù)如表1所示。

圖1 保壓容器結(jié)構(gòu)圖和三維幾何模型

表1 保壓容器參數(shù)

1.2 保壓容器有限元分析

本文利用有限元軟件ANSYS中的Workbench模塊對全海深保壓容器進行有限元分析研究[8～10]。在有限元分析之前，首先，需要對保壓容器的材料進行定義，根據(jù)課題的要求，定義保壓容器筒體和保壓容器端蓋的材料均為TC4鈦合金，彈性模量為110GPa，泊松比為0.34，螺栓的材料為45鋼；螺栓和保壓容器端蓋之間的接觸單元設置為Frictional單元，保壓容器端蓋和保壓容器筒體的接觸單元設置為No separation單元。然后分別對保壓容器筒體、保壓容器端蓋及連接螺栓進行網(wǎng)格劃分，如圖2所示，并對保壓容器施加約束，最后，利用有限元分析ANSYS軟件，開展模擬全海深超高壓環(huán)境下的保壓容器性能分析。

圖2 端蓋和筒體網(wǎng)格劃分

圖3和圖4分別給出了在全海深超高壓環(huán)境下（120MPa），保壓容器筒體、保壓容器端蓋及連接螺栓的應力、應變云圖。由圖4和圖5可知，保壓容器端蓋的應力、應變均很小，在初步設計中所選取的8個螺栓的應力均小于螺栓材料的許用應力。保壓容器筒體所受的應力沿軸線分布，保壓容器筒體中間部分所受應力較大，同時，應力沿著筒體徑向由內(nèi)向外逐漸變小，最大應力位于保壓容器筒體的內(nèi)壁，對應該位置的最大應力值為320.66MPa，最大應變?yōu)?.003mm，屬于安全范圍，但還可以進一步優(yōu)化。

圖3 保壓容器筒體、保壓容器端蓋及連接螺栓應力云圖

圖4 保壓容器筒體、保壓容器端蓋及連接螺栓應變云圖

2 響應面建模與分析

2.1 實驗設計

為研究保壓容器結(jié)構(gòu)設計參數(shù)對保壓容器最大應力值及重量的影響，考慮到保壓容器所需的儲存容積為V=150ml。因此，本文選取保壓容器的3個設計參數(shù)（保壓容器的內(nèi)徑D1、保壓容器壁厚δ1、保壓容器底部厚度δ2）為設計變量，應用實驗設計方法來揭示保壓容器的3個設計參數(shù)與保壓容器最大應力值及重量之間的關系，從而為保壓容器的優(yōu)化設計提供參考。綜合考慮實驗設計的各種方法，本文采用中心組合設計來安排實驗方案。在實驗設計之前，通過單設計變量因子實驗確定出各設計變量的變化范圍，單設計變量因子實驗結(jié)果表明：保壓容器的3個設計參數(shù)D1、δ1、δ2對保壓容器最大應力值和重量均有影響。由單設計變量因子實驗確定出實驗中各設計變量的變化范圍（為了便于分析將D1、δ1、δ2分別由X1、X2、X3代替）各因素的水平與編碼如表2所示，實驗方案如表3所示。

表2 各因素水平與編碼

表3 實驗方案

2.2 實驗結(jié)果分析

以全海深沉積物保壓容器內(nèi)徑D1、壁厚δ1、底部厚度δ2為設計變量，分別以保壓容器重量和最大應力為響應量，進行響應面回歸分析，分析結(jié)果如表4、表5所示。

由表4、表5可知，回歸分析檢驗是顯著的，能較好地反映保壓容器重量和最大應力與保壓容器內(nèi)徑D1、壁厚δ1、底部厚度δ2的關系，即通過改變保壓容器內(nèi)徑D1、壁厚δ1、底部厚度δ2可以有效的減少保壓容器重量和最大應力值。保壓容器內(nèi)徑D1、壁厚δ1、底部厚度δ2對保壓容器重量和最大應力的響應曲面及其等高線如圖5所示。

對各因素進行多遠回歸擬合，得到以保壓容器重量F1（X1，X2，X3）及最大應力F2（X1，X2，X3）為目標函數(shù)的回歸模型：

該模型較好地反映了保壓容器重量和最大應力與保壓容器內(nèi)徑D1、壁厚δ1、底部厚度δ2的關系，所得的回歸模型能較好地預測保壓容器重量和最大應力隨各參數(shù)變化規(guī)律，因此可以利用該回歸模型確定保壓容器內(nèi)徑D1、壁厚δ1、底部厚度δ2。

3 保壓容器結(jié)構(gòu)參數(shù)多目標優(yōu)化與分析

根據(jù)回歸模型的二次多項式，以保壓容器重量和最大應力值最小化為優(yōu)化目標，以保壓容器內(nèi)筒體積和許用應力為約束條件，經(jīng)過優(yōu)化后得到保壓容器重量和最大應力值最小的條件：X1=-1，X2=-0.711，X3=-1，此時，保壓容器重量為2.066kg，最大應力為317.243MPa，轉(zhuǎn)化為實際參數(shù)，即D1=25mm，δ1=8.87mm，δ2=15mm，與原始保壓容器重量減少了4.79%，最大應力減少了1.07%。

表4 保壓容器重量回歸分析結(jié)果

表5 保壓容器最大應力回歸分析結(jié)果

圖5 三維響應面曲面和等高線圖

應用有限元軟件ANSYS中的Workbench模塊對優(yōu)化后的保壓容器進行有限元分析。其分析結(jié)構(gòu)與優(yōu)化結(jié)果相吻合。圖6、圖7為對應優(yōu)化后保壓容器應力應變圖。

圖6 保壓容器筒體、保壓容器端蓋及連接螺栓應力云圖

圖7 保壓容器筒體、保壓容器端蓋及連接螺栓應變云圖

4 結(jié)論

針對全海深沉積物保壓容器工作環(huán)境與性能要求，本文應用響應面法和有限元分析方法開展了保壓容器幾何結(jié)構(gòu)設計參數(shù)多目標優(yōu)化設計，通過有限元軟件ANSYS中的Workbench模塊對保壓容器進行了有限元分析，并利用BOX-BEHNKEN試驗設計和響應面法建立了全海深沉積物保壓容器重量和最大應力與內(nèi)徑D1、壁厚δ1、底部厚度δ2之間的數(shù)學模型。其擬合的回歸模型能很好地反映保壓容器重量和最大應力與內(nèi)徑D1、壁厚δ1、底部厚度δ2的關系。通過對該模型進行優(yōu)化，結(jié)果表明，優(yōu)化后的全海深沉積物保壓容器重量減少了4.79%，最大應力減少了1.07%。

[1]王元文.高壓容器殼體的優(yōu)化設計[J].機械設計與制造,2005(12)：16-17.

[2]夏鋒社,朱哲,淡勇.高壓容器筒體結(jié)構(gòu)的最優(yōu)化設計[J].西安石油大學學報(自然科學版),2010,25(1)：81-83.

[3]汲方林,王有杰,陳萬勝,等.基于SolidWorks Simulation對高壓容器外形結(jié)構(gòu)尺寸的優(yōu)化分析[J].機械,2009,36(8)：33-36.

[4]Tauchert T R.Optimum Design of a Reinforced Cylindrical Pressure Vessel[J].Journal of Composite Materials,1981,15(5)：390-402.

[5]Mian H H, Wang G, Dar U A, et al. Optimization of Composite Material System and Lay-up to Achieve Minimum Weight Pressure Vessel[J].Applied Composite Materials,2013,20(5)：873-889.

[6]Joung T H, Lee J H, Nho I S, et al. A study on the pressure vessel design, structural analysis and pressure test of a 6000 m depth-rated unmanned underwater vehicle[J].Ships & Offshore Structures,2005,3(3)：205-214.

[7]Ng R K H, Yousefpour A, Uyema M, et al. Design, Analysis,Manufacture, and Test of Shallow Water Pressure Vessels Using E-Glass/Epoxy Woven Composite Material for a Semi-Autonomous Underwater Vehicle[J].Journal of Composite Materials,2002,36(21)：2443-2478.

[8]Joung T H, Lee J H, Nho I S, et al. A study on the design and manufacturing of a deep-sea unmanned underwater vehicle based on structural reliability analysis[J].Ships & Offshore Structures,2009,4(1)：19-29.

[9]Bjeletich J G, Morton T M. Fracture mechanics technology for optimum pressure vessel design[J].Astm Special Technical Publication,1973.

[10]Paknahad A,Nourani R. Mix Model of FE Method and IPSO Algorithm for Dome Shape Optimization of Articulated Pressure Vessels Considering the Effect of Non-geodesic Trajectories[J].Journal of The Institution of Engineers (India)：Series C,2014,95(2)：151-158.