基于參數(shù)化子模型的SWATH船典型節(jié)點抗疲勞設計

任慧龍，張清越，胡雨蒙，江雪云

(1.哈爾濱工程大學 船舶工程學院，黑龍江 哈爾濱150001;2.哈爾濱工程大學 多體船技術國防重點學科實驗室，黑龍江 哈爾濱150001)

0 引 言

小水線面雙體船(small waterplane area twin hull ship，SWATH)是近年來發(fā)展起來的一種新概念船型，因其良好的耐波性及快速性，以及能夠勝任多種工作的優(yōu)點而受到廣泛關注［1-2］。但與此同時，為SWATH 船型帶來諸多優(yōu)勢的細長支柱體結構同樣是整船中結構最為薄弱的部位。特別是連接橋過渡區(qū)域，由于截面發(fā)生突變，剛度不連續(xù)，應力場不連續(xù)［3］，使得該區(qū)域在橫浪和斜浪狀態(tài)下不斷受到較大的交變橫向對開力和片體不同步縱搖扭矩的作用，應力水平很高，應力集中顯著。此外，SWATH船型對于重量的變化非常敏感，船體結構的尺寸不宜過大，因此這些高應力部位通常都會使用高強度鋼來滿足其屈服和屈曲強度要求。綜上所述，SWATH船型連接橋過渡區(qū)域的疲勞強度問題十分突出，在結構設計初期就應該予以足夠的重視。文獻［4］對多種影響船體結構應力的相關因素進行比較分析，提出一種小水線面雙體船典型節(jié)點的抗疲勞設計，為設計者提供參考;文獻［5］針對小水線面雙體船采用規(guī)范法、譜分析法、設計波法進行疲勞強度評估，并對其中的關鍵節(jié)點進行了以滿足疲勞強度為目標的優(yōu)化設計，取得了一定的成果。

通常設計者依據規(guī)范和經驗擬定的結構形式和構件尺寸未必能夠滿足結構強度的要求，而小水線面雙體船的相關規(guī)范尚不完善，且疲勞強度和結構細節(jié)密切相關，對于疲勞載荷和結構應力響應的估算十分困難，疲勞壽命的準確計算還有賴于水動力和有限元分析技術?；谄趶姸确治龅膬?yōu)化設計是給出較好抗疲勞結構方案的有效途徑，但此方法需要對優(yōu)化模型進行大量的迭代計算。針對這些問題和難點，本文提出一種結合參數(shù)化子模型有限元分析方法和基于Miner 線性累積損傷理論及Weibull分布連續(xù)模型的疲勞設計波法的優(yōu)化方法，通過快速生成批量模型的參數(shù)化建模技術，準確描述局部模型外載荷的子模型法及能夠較好反映結構細節(jié)而計算量又相對較小的疲勞設計波分析方法來快速有效地解決結構抗疲勞設計問題。

1 參數(shù)化子模型法

參數(shù)化子模型法是針對局部結構優(yōu)化問題提出的一種快速準確的有限元計算方法，它基于整體結構模型的有限元分析，從計算結果中提取局部模型切割邊界的位移作為約束條件施加在子模型上，從而得到子模型較為精確的外載荷響應計算結果。通過參數(shù)化建模手段建立子模型的參數(shù)化模型，使優(yōu)化模型的大量迭代計算得以實現(xiàn)。

1.1 整船有限元分析

SWATH 船型結構特殊。較小的船長型深比，深入水下的潛體使得SWATH 船的縱彎曲強度有較大的富裕。但較大的側面積，重力浮力沿船寬分布情況和薄弱的支柱體結構使其在遭受橫浪時產生對開彎矩作用及斜浪航行時產生復雜扭轉的情況下非常危險。特別是連接橋過渡結構，主要承受了巨大的交變載荷，很容易發(fā)生疲勞破壞。

本文依據CCS 頒布的《小水線面雙體船指南2005》［6］中的相關要求，使用通用有限元軟件MSC. Patran 建立某小水線面雙體船整船有限元模型，計算應力云圖如圖1 所示。在初期的總強度分析計算中，可以發(fā)現(xiàn)該船的高應力區(qū)域集中在支柱體橫艙壁(或橫框架)上端與舷臺相交的結構過渡處。特別是在橫浪工況下，該區(qū)域的應力水平很高，應力集中嚴重，如圖2 所示。因此，可將此局部結構作為疲勞強度校核和優(yōu)化的典型節(jié)點。

圖1 整船計算應力云圖Fig.1 Stress calculation of whole ship

圖2 舷臺與支柱體橫艙壁應力云圖Fig.2 Stress calculation of typical detail

1.2 子模型法

子模型法基于圣維南原理，即如果實際分布載荷被等效載荷代替后，應力和應變只在載荷施加的位置附近有改變，因此只要子模型切割邊界避開載荷集中及應力集中位置，子模型內部就可以得到較精確的解［7］。子模型法包含以下步驟:

1)整體模型分析，對整體粗網格模型進行計算分析，并算出局部模型切割邊界的位移響應;

2)根據結構實際構造建立局部精細網格模型;

3)根據整體模型分析切邊的位移響應通過線性插值法加載至局部精細子模型的邊界，作為其邊界條件;

4)局部精細模型自身的載荷及邊界條件不變，進行有限元分析。本文從MSC. Patran 整船有限元計算結果中提取切割邊界位移響應，采用APDL 參數(shù)化編譯語言編寫了Ansys 子模型自動加載接口，為后續(xù)用以優(yōu)化迭代計算的參數(shù)化模型提供準確的外載荷模擬。

1.3 參數(shù)化模型

參數(shù)化模型是一種通過變量參數(shù)及建模命令建模，可以快速生成系列模型的技術手段，它包含幾何模型參數(shù)化生成，有限元網格參數(shù)化劃分，材料屬性參數(shù)化定義，載荷和邊界條件自動添加及參數(shù)化后處理等幾塊內容。本文采用Ansys 參數(shù)化編譯語言APDL 建立小水線面雙體船連接橋過渡區(qū)域典型節(jié)點的參數(shù)化模型，選擇影響該區(qū)域疲勞強度的幾何參數(shù)和尺寸參數(shù)作為優(yōu)化變量，為后續(xù)的優(yōu)化迭代計算打下基礎。本文建立的局部結構參數(shù)化模型和子模型法加載計算應力結果如圖3 和圖4 所示。

圖3 參數(shù)化子模型效果圖Fig.3 Parametric sub-model

圖4 參數(shù)化子模型計算應力云圖Fig.4 Stress calculation of parametric sub-model

2 典型節(jié)點疲勞強度分析

基于譜分析的疲勞評估直接計算雖然具有較高精度，但是需要對波浪譜有效波浪頻率范圍內的多個規(guī)則波和浪向角分別進行計算，計算量過于龐大，優(yōu)化模型的迭代計算難以實現(xiàn);而基于經驗公式的規(guī)范簡化計算方法雖然計算簡便，但是不涉及有限元分析計算，難以體現(xiàn)結構的細節(jié)特征;CCS 《小水線面雙體船指南2005》中的載荷計算簡化公式及有限元加載方式雖然已經較充分地考慮到了SWATH 船的結構及載荷的特點，但針對尺度相近的具體船型卻難以體現(xiàn)不同SWATH 船的差異性。綜上所述，研究認為基于設計波法的疲勞評估直接計算方法通過勢流理論及船舶航行過程的實際海況長期分析得到設計波載荷并通過有限元分析計算得到典型節(jié)點的疲勞應力范圍，能夠反映結構的細節(jié)特征，疲勞載荷更加真實，而計算量又相對較小，適合進行優(yōu)化迭代計算。

2.1 設計波系統(tǒng)的確定

根據CCS《小水線面雙體船指南2005》中規(guī)定的載荷工況及總強度計算結果，認為SWATH 船主要考慮的載荷參數(shù)為:橫向對開彎矩、不同步縱搖扭矩、縱向彎矩、水平扭矩、縱中剖面剪力。由文獻［5］及文獻［8］的計算結果可知在橫浪產生的橫向對開力作用下產生的累積損傷在SWATH 船型的整個疲勞壽命周期內占主要成分，其余載荷的累積損傷相對較小，可以忽略?？紤]到本文的最終目的是針對SWATH 船連接橋過渡區(qū)域的結構進行抗疲勞設計，僅以橫浪狀態(tài)下的橫向對開彎矩作為主要載荷參數(shù)的設計波計算結果已經能夠表征結構的疲勞強度。指南中規(guī)定了橫浪狀態(tài)的時間分配系數(shù)為0.25。

根據本研究所自主開發(fā)的基于三維勢流理論的波浪載荷計算程序計算船體在單位波幅規(guī)則波下的波浪載荷傳遞函數(shù)，根據全球波浪散布圖及數(shù)理統(tǒng)計理論進行主要控制載荷的長期預報，取10-8概率水平的預報值作為疲勞載荷設計值。主要控制載荷的幅頻響應最大值對應的浪向即為設計波的浪向β和ω，而設計波的波長為:

設計波的波幅為載荷設計值與主要載荷參數(shù)幅頻響應最大值之比:

2.2 S － N 曲線選取

結構在某一應力范圍水平為S 的單一循環(huán)載荷作用下達到疲勞破壞所需循環(huán)次數(shù)為N，則稱該結構在應力范圍水平為S 時的疲勞壽命為N，一般根據疲勞試驗結果擬合得到S － N 曲線，記為:

式中，m 和A 為疲勞試驗所得參數(shù)。

對上式等號兩邊取對數(shù)，得

即為常用的S-N 曲線的雙對數(shù)線性模型。CCS《船舶結構疲勞強度指南》［9］中給出了8 條S-N 曲線對應的參數(shù)值及焊接和結構形式。但SWATH 船型連接橋過渡處的結構在規(guī)范中無法找到對應的結構形式，S-N 曲線需由疲勞試驗確定。在此取文獻［10］中的SWATH 船連接橋過渡結構疲勞試驗結果作為疲勞評估的S-N 曲線。

2.3 應力范圍的Weibull 分布

在疲勞評估的設計波法中，設計壽命期應力范圍的長期分布定義為Weibull 分布，分布函數(shù)見下式:

式中:h 為形狀參數(shù);q 為尺度參數(shù)。根據CCS 《小水線面雙體船指南》，h =1 比較有意義。q 則可以用回復期內超越概率為P(S ≥SL)=1/NL的應力范圍SL來表示:

式中:NL為回復期內的應力范圍循環(huán)次數(shù)，依據CCS 《船舶結構疲勞強度指南》，對于設計壽命為20 年的艦船，應取0.6 ×108。同時該指南中規(guī)定應力范圍應在t × t 的精細網格有限元模型下按Lagrange 插值計算得到的熱點應力，插值方法如下:

其中στ/2和σ3τ/2都是按式(9)計算得到的，式(9)中的σ1、σ3、σ4是靠近焊趾的4 個有限元節(jié)點的應力值，而系數(shù)C1、C2、C3、C4按下式計算:

式中:x 為插值點到焊趾的距離;xi為有限元節(jié)點到焊趾的距離。

2.4 累積損傷度的計算

Weibull 應力范圍分布模型對應的線性累積損傷度可以表示為:

代入式(3)和式(5)可得:

對應的GAMMA 函數(shù)值可通過查表得到。

3 優(yōu)化模型

在進行結構優(yōu)化分析中，必須要建立合適的數(shù)學優(yōu)化模型。一個完整的數(shù)學優(yōu)化模型包括所定義的目標函數(shù)、優(yōu)化變量及約束條件等。通常數(shù)學優(yōu)化模型可以表示為如下形式:

式中，X=(x1，x2，…，xm)T為設計變量組成的列向量空間;f(X)為定義的目標函數(shù);gi(X)為定義的約束函數(shù)，包括等式約束和不等式約束;和分別為設計變量xi的上下限值，其大小通常可以根據實際設計情況而定。

針對本文中的抗疲勞設計問題，目標函數(shù)取為結構的累積損傷度D。連接橋過渡區(qū)域的圓弧半徑越大，結構過渡就越平滑，顯然疲勞強度更高，但是過大的圓弧半徑對于SWATH 船的重量控制和阻力性能都不利;過渡區(qū)域的構件尺寸設計未必最合理，部分區(qū)域的構件尺寸對于結構強度的貢獻并不明顯。因此，優(yōu)化變量取連接橋過渡區(qū)域的圓弧半徑R 及該區(qū)域各板構件的厚度ti。同時，對于SWATH 船這類吃水對結構重量非常敏感的船型，盡量地減輕船體重量非常必要，在優(yōu)化時還應以結構重量作為優(yōu)化模型的約束條件。數(shù)學優(yōu)化模型表示如下:

式中:VP 為體積參數(shù)，是優(yōu)化后體積與初始體積之比;R 為過渡圓弧的半徑;Ti為各板構件厚度。

4 優(yōu)化結果與分析

本文選取的算例為某5 000 t 級SWATH 船，根據前文所述的方法及理論建立目標船的參數(shù)化模型及優(yōu)化模型，在Ansys 的優(yōu)化設計平臺中進行優(yōu)化迭代計算。優(yōu)化算法選擇為零階子問題逼近法，此方法迭代效率較高，且不容易陷入局部最優(yōu)解，能夠滿足本次優(yōu)化迭代的要求。目標函數(shù)、優(yōu)化變量和約束條件的選取和初始值見表1。

表1 各變量初始值Tab.1 The initial value of the variables

經過28 步優(yōu)化迭代，計算收斂。目標函數(shù)值隨迭代序列對變化如圖5 所示，部分迭代序列的計算結果見表2。

圖5 目標函數(shù)隨迭代序列變化Fig.5 Optimization iteration of objective function

表中的優(yōu)化序列28 即為本次優(yōu)化的最優(yōu)設計方案。累積損傷度降低為0.62，大大提高了結構的疲勞壽命，同時局部結構重量還有一定的降低。

表2 部分迭代序列計算結果Tab.2 Partial results of optimization iteration

從優(yōu)化結果中可以發(fā)現(xiàn)，過渡圓弧半徑的增大對于降低該區(qū)域應力集中的作用非常明顯。同時，靠近應力集中點的橫艙壁板1、橫艙壁板2、外板加厚板的尺寸對累積損傷度的影響也很明顯。相對的，遠離應力集中點的構件尺寸對于疲勞強度的貢獻比較低，可以適當減小尺寸。從理論上分析，SWATH船的支柱體結構在橫浪作用下類似于一個深入水下的懸臂梁，一端受到與梁相垂直的力的作用。決定支柱體上端與舷臺連接處應力集中水平的主要因素就是該區(qū)域的幾何連續(xù)程度和水平剖面的抗彎剛度。因此，提升該剖面上的構件尺寸能夠提高抗彎剛度;而增大過渡圓弧的半徑不僅能夠使幾何過渡更加平滑，同時也增大了剖面尺寸。

5 結 語

本文基于參數(shù)化子模型法及疲勞強度評估的設計波法，提出了一套結構抗疲勞設計的優(yōu)化方法，并結合某SWATH 船進行優(yōu)化設計，得出以下結論:

1)參數(shù)化子模型法為局部結構應力響應的準確計算和優(yōu)化迭代問題提供了有效途徑;

2)基于設計波法的疲勞強度評估方法能夠真實反映結構細節(jié)，同時計算量較小，適合作為結構設計初期的抗疲勞設計衡準;

3)SWATH 船連接橋過渡區(qū)域的疲勞強度問題嚴重，增加過渡圓弧的半徑和應力集中點附近的板材厚度對提高該結構疲勞強度幫助較大，本文提出的抗疲勞設計方案為SWATH 船連接橋過渡區(qū)結構設計提供了一定參考。

［1］林偉國，朱云翔，范井峰，等.小水線面雙體船的發(fā)展及在海軍艦船領域中的應用前景［J］.船舶，2007，18(3):1-5.LIN Wei-guo，ZHU Yun-xiang，F(xiàn)AN Jing-feng，et al. Development of small waterplane area twin hull and its application prospect in navy［J］.Ship and Boad，2007，18(3):1-5.

［2］ANANTHA S V.Numeric design and evaluation of swath form［J］.International Shipbuilding Progress，2002，49(2):95-125.

［3］徐志亭，柳存根，王曉宇.小水線面雙體船疲勞強度的規(guī)范計算法研究［J］.中國造船，2011，52(4):66 -73.XU Zhi-ting，LIU Cun-gen，WANG Xiao-yu. Research on rules′ method in fatigue assessment of SWATH［J］.Shipbuilding of China，2011，52(4):66 -73.

［4］陳穎，岳亞霖，崔維成.小水線面雙體船典型節(jié)點抗疲勞設計［J］.艦船科學技術，2010，32(6):8 -13.CHEN Ying，YUE Ya-lin，CUI Wei-cheng. Anti-fatigue design of SWATH typical local details［J］.Ship Science and Technology，2010，32(6):8 -13.

［5］徐志亭.小水線面雙體船結構疲勞分析及典型節(jié)點優(yōu)化研究［D］.上海:上海交通大學，2011.

［6］中國船級社.小水線面雙體船指南［S］.北京:人民交通出版社，2005.ChinaClassification Society.Guid elines of Small Waterplane Area Twin Hull Craft［S］.Beijing:China Communication Press，2005.

［7］史戰(zhàn)新.基于Ansys 子模型法的肘板結構優(yōu)化［J］.艦船科學技術，2014，36(8):19 -26.SHI Zhan-xin. The bracket structure optimization based on Ansyssub-model method［J］.Ship Science and Technology，2014，36(8):19 -26.

［8］甄春博.小水線面雙體船結構疲勞評估方法研究［D］.哈爾濱:哈爾濱工程大學，2009.ZHEN Chun-bo.Study on the methodofthe structural fatigue assessment of SWATH［D］. Harbin:Harbin Engineering University，2009.

［9］中國船級社.船舶結構疲勞強度指南［S］.北京:人民交通出版社，2007.China Classification Society.Guidelines for fatigue strength of ship structure［S］.Beijing:China Communication Press，2005.

［10］甄春博，任慧龍，馮國慶，等.小水線面雙體船連接橋結構疲勞強度試驗［J］.哈爾濱工程大學學報，2012，33(4):414-418.ZHEN Chun-bo，REN Hui-long，F(xiàn)ENG Guo-qing，et al.Fatigue strength test of a SWATH ship′ s cross-deck structure［J］. Journal of Harbin Engineering University，2012，33(4):414 -418.