不同材料上層建筑的超級游艇在波浪中結構強度分析*

張水林　裴志勇　蔡　薇　吳衛國

(武漢理工大學交通學院　武漢　430063)

?

不同材料上層建筑的超級游艇在波浪中結構強度分析*

張水林裴志勇蔡薇吳衛國

(武漢理工大學交通學院武漢430063)

摘要：采用基于勢流理論的三維特異點分布法計算指定海況下超級游艇的運動和作用于船體表面的水壓力，將得到的外載荷等效轉換到結構模型上，采用有限元程序分析超級游艇的上層建筑分別采用鋁合金和普通鋼時波浪中的響應特性.計算結果表明,超級游艇的上層建筑對主船體的剖面應力分布有較大的影響，并且不同材料上層建筑影響的程度不同；設置多層甲板會改變剖面應力分布；不同材料上層建筑的超級游艇在相同的外載荷作用下引起的船體結構應力和變形不同.

關鍵詞：超級游艇；上層建筑；結構強度；剖面應力分布

張水林(1990- )：男，碩士生，主要研究領域為船體結構強度

*工信部高技術船舶科研項目資助(批準號：工信部聯裝[2010]337號)

0引言

超級游艇十分重視安全性和舒適性，而在波浪中的結構響應是決定超級游艇的安全性和舒適性的一個重要因素.目前，針對游艇在外載荷作用下的結構響應，國內外已開展諸多研究[1-4].游艇結構響應的研究主要是針對小型游艇，較少涉及超級游艇.此外，當前的研究主要是基于船級社的規范和相關組織的規定，未充分考慮具體船型特點.由于超級游艇尺度較大，作用于船體的波浪載荷十分依賴于船體外形.因此，為將超級游艇船型特點考慮在內，以全面真實地模擬超級游艇在其工作環境中的結構強度，亟須對超級游艇在波浪中的結構強度展開深入研究.

為降低船體結構質量，超級游艇廣泛采用鋁合金，并且采用鋁合金制造的船舶的營運費用和維護費用相比同類型鋼質船舶低[5].超級游艇的上層建筑采用鋁合金，一方面可以降低船舶質量，提高超級游艇的快速性；另一方面也可以降低其重心高度，提高船舶穩性.主船體為普通鋼的超級游艇，由于2種材料物理性質的差異，導致其在波浪的結構強度不同于全鋼船舶.因此，需要對其在波浪中的結構強度展開研究，并與相同尺度的全鋼超級游艇的結構強度進行比較，以探究上層建筑采用鋁合金對船體結構在波浪中結構強度的影響.

本文將載荷/運動計算和結構響應分析結合起來，采用基于勢流理論的三維特異分布點法[6-7]計算指定海況條件下作用于超級游艇的載荷和相應的船體運動，將得到的外載荷施加到整船結構模型上，采用結構有限元程序分析相應的結構響應.在此基礎上，研究上層建筑分別采用不同材料的超級游艇在波浪中結構強度的差異，主要比較了在相同外載荷作用下超級游艇典型剖面應力分布、結構的應力和變形水平等方面的異同，并分析了產生這些差異的原因.

1分析流程

為了比較不同材料上層建筑的超級游艇在波浪中的結構強度，整個分析流程由以下6個步驟組成.

1) 建立用于結構分析的全船有限元模型.

2) 模擬指定裝載狀況下的質量分布，包括船體鋼料質量、舾裝質量、壓載水質量等.

3) 在指定海況和載況下進行載荷/運動計算，得到作用于船體的隨時間而變的外壓力和船體的運動加速度.

4) 將載荷/運動計算得到的作用于船體表面的外壓力以及由于船體運動產生的慣性力轉化為等效節點力.

5) 將等效節點力施加到結構計算模型上，進行結構計算分析.

6) 改變上層建筑材料，采用4)得到的等效外載荷進行結構分析.

為了比較相同結構形式下，超級游艇的上層建筑分別采用鋼和鋁時結構強度的異同，假定2種情況下作用于船體的外載荷相同.當然，實際由于材料的不同會引起作用于船體的載荷有所不同，文中暫未考慮.

2計算實例

2.1模型概況

本文研究的對象是一艘總長為64.6 m的超級游艇，其主尺度要素見表1.該超級游艇的上層建筑采用鋁合金，主船體結構采用普通鋼.

表1　超級游艇的主尺度　m

載荷/運動計算和結構響應分析階段使用相同的模型有利于避免因計算模型的差異而引起的不平衡載荷而影響結構分析的準確性.整船彈性有限元分析模型及邊界條件見圖1.計算模型中船體外板、甲板和艙壁板等板件采用板殼單元模擬，縱桁、肋骨等骨材采用彈性梁單元模擬，整個模型共有21 189個節點，22 761個板殼單元，23 044個彈性梁單元.

圖1　全船有限元計算模型和邊界條件

2.2載荷/運動計算

超級游艇的船體通常為細長體結構，在迎浪時船體所受的垂向彎矩是最基本的載荷.因此選取迎浪狀態對船體結構強度進行分析.波浪取規則正弦波，波長取艇長，波高的選取主要是考慮游艇在營運過程中常見的波浪大小結合結構強度分析需要而設置.滿載狀況下波浪參數見表2.

表2　滿載狀況下波浪參數

由于波浪和船體的運動，不同時刻作用于船體的壓力分布是隨時間變化的.本文選取大于一個波浪周期的時間段，并將其均分為20個時刻，將每一時刻作用于船體的外力積分即可得到各個時刻作用于船體的垂向彎矩，不同時刻作用于船體的垂向彎矩沿船長的分布見圖2.在t1時刻，波峰位于船中，作用于船體的彎矩為中拱彎矩.隨著波浪和船體的運動，波峰逐漸向船首移動，作用于船體的彎矩從中拱逐漸變為中垂.在t6時刻至t7時刻之間波谷移至船中，此時作用于船體的中垂彎矩達到最大值.然后波峰向船中移動，作用于船體的彎矩逐漸從中垂變為中拱，在t16時刻波峰移動至船中，作用于船體的中拱彎矩達到最大.

圖2　各時刻垂向彎矩分布

2.3結構響應計算結果

將各個不同時刻作用于船體表面的外壓力和慣性力等效轉換到結構模型上，采用通用商業有限元程序進行線性結構強度分析.先對上層建筑為鋁合金的鋼鋁混合模型進行結構強度分析，然后分析上層建筑采用普通鋼的全鋼模型在相同外載荷作用下的結構強度，分別得到典型剖面應力分布、結構的應力和變形的計算結果.

2.3.1典型剖面應力分布

為了提供更寬敞的內部空間，以滿足布置、居住、飲食和娛樂等要求，超級游艇往往會設置多層甲板.由于設置了多層甲板，各層甲板的局部作用會對剖面的應力分布特性產生影響.為了研究典型剖面應力分布特性，選取作用于船體的中拱彎矩最大的時刻(t16時刻)作為典型時刻，取靠近船中的某一剖面.該剖面在t16時刻的舷外水壓力分布見圖3，圖中箭頭方向表示壓力的方向，與船體外表面垂直；箭頭長短表示作用于艇體表面的壓力幅值. 舷外水壓力作用于下甲板以下的船體外殼，隨著型深高度的增加逐漸減小.

圖3　t16時刻橫剖面外壓力分布

采用有限元分析得到該剖面上各點的縱向應力，以各點的縱向應力值為橫坐標，各點的垂向位置為縱坐標，分別繪出鋼鋁混合模型和全鋼模型中該剖面上各點在有無上層建筑的縱向應力沿型深方向的分布，見圖4.

圖4　典型橫剖面縱向應力沿型深分布

由圖4可見，由于所遭受的外載荷和船體結構都相同，鋼鋁混合模型和全鋼模型剖面應力分布非常吻合.船底至舭部區域應力分布波動較為劇烈，舭部至下層甲板，以及下甲板至主甲板之間的船體結構的縱向應力與垂向位置之間均呈線性變化，在各層甲板上的縱向應力與其橫向位置相關.

由圖4可見，船底至舭部剖面應力分布仍存在波動，但波動程度較圖4a)和緩；舭部至下甲板，以及各層甲板之間的船體結構的縱向應力與垂向位置之間均近似呈線性變化，在各層甲板上的縱向應力與橫向位置相關.由于上層建筑的影響，圖4b)中主船體以下剖面應力分布不同于圖4a)，并且，由于上層建筑材料的不同，導致2個模型中典型剖面的縱向應力分布有所不同，主要表現在船底部應力值、應力分布曲線的斜率以及各層甲板處應力的變化范圍等方面.

2.3.2結構應力與變形

由于普通鋼與鋁合金的物理性質不同，因而上層建筑分別采用普通鋼和鋁合金時，其參與抵抗外載荷的程度不一.相比而言，采用普通鋼時的船體結構的剛度大于采用鋁合金，因此在相同的波浪載荷作用下，其船體結構的應力水平和變形程度有所不同.選取作用于船體垂向彎矩最大的時刻(t16時刻)作為一種典型狀態，從上層建筑的應力分布和變形、主船體的應力分布和變形這4個方面，比較在相同外載荷作用下超級游艇上層建筑分別采用鋁合金和普通鋼時不同的結構強度.

上層建筑分別采用鋁合金和普通鋼時，上層建筑的Mises應力分布見圖5.其中鋼鋁混合模型的最大應力為29.7 MPa，全鋼模型的最大應力為54.8 MPa.變形見圖6.其中鋼鋁混合模型的最大變形值為10.3 mm，全鋼模型的最大變形值為5.06 mm.

圖5　上層建筑應力分布

圖6　上層建筑變形

上層建筑分別采用鋁合金和普通鋼時，主船體Mises應力分布見圖7.其中鋼鋁混合模型的最大應力為56.2 MPa，全鋼模型的最大應力為57.7 MPa；變形見圖8.其中鋼鋁混合模型的最大變形值為10.9 mm，全鋼模型的最大變形值為9.29 mm.

圖7　主船體應力分布

圖8　主船體變形

3計算結果分析比較

3.1剖面應力分布分析比較

未設置多層甲板的船體剖面的縱向應力沿斷面呈線性分布，即不同位置處的應力與其距中和軸的位置成正比.在中拱彎矩最大的時刻(t16時刻)，設置多層甲板的超級游艇的典型剖面應力分布見圖4.

沒有上層建筑的主船體在相同外載荷作用時，從其剖面應力分布(圖4a))中可以看出其分布特性.從船底至舭部區域：船體由于受中拱彎矩的作用，使得船底構件受壓.船體外板上各點的總縱彎曲應力隨型深的增加逐漸增大，并按線性變化.而由圖3知此區域的船體外板所受到的舷外水壓力較大，使這部分船體構件產生橫向彎曲，船體外板中產生較大的局部應力.總縱彎曲應力與局部應力疊加得到船體構件的縱向應力.由于局部應力的影響使得縱向應力不是嚴格的線性分布，部分觀測點波動較為劇烈.

舭部到主甲板之間的船體外板在中拱彎矩作用下產生的總縱彎曲應力隨型深的增加線性增加；舷外水壓力隨著型深增加急劇減小，船體外板產生的局部應力迅速減小；而總縱彎曲應力和局部應力疊加得到的縱向應力線性增加.因此，船體外板的縱向應力主要由總縱彎曲應力組成，局部作用力的影響較小.所以，在沒有上層建筑的主船體剖面應力分布(圖4a))中此區域船體外板的縱向應力分布表現出良好的線性.由于甲板上各處局部變形不同導致各點的局部應力不同，致使下甲板和主甲板上各點的縱向應力依賴其橫向位置.

從包括上層建筑的剖面應力分布(圖4b))中可以看出有無上層建筑以及不同上層建筑時其剖面應力分布的特性和差異.從船底至舭部區域，是鋼鋁混合模型與全鋼模型其應力分布都存在波動，但較沒有上層建筑的應力分布的波動程度小. 鋼鋁模型此區域船體構件的縱向應力變化范圍較全鋼模型大.原因在于上層建筑承擔了一部分載荷，并且采用普通鋼的上層建筑承擔的載荷更多，因而主船體承擔的載荷更小.

舭部至下甲板區域船體外板的縱向應力分布呈現出近似線性，個別觀測點存在偏離.鋼鋁混合模型與全鋼模型的縱向應力分布的斜率不同，全鋼模型的分布斜率稍大.在下甲板處上的各點的縱向應力也與其橫向位置有關，但全鋼模型的縱向應力變化范圍小于鋼鋁混合模型.因為此時舷外水壓較小，引起的局部作用力較小，船體外板的縱向應力以彎曲應力為主.

下甲板至主甲板間船體外板的縱向應力主要是總縱彎曲應力，鋼鋁混合模型與全鋼模型的縱向應力分布均表現出良好的線性，兩者的分布斜率較為接近.但是2個模型在此區域的船體外板的縱向應力之間相差較大，并且主甲板上鋼鋁混合模型的縱向應力變化范圍較大.原因是由于鋼鋁混合模型主船體承擔的外載荷較大，因而其在船體外板中引起的總縱彎曲應力較大.

主甲板至陽光甲板間船體外板的縱向應力近似按線性分布，鋼鋁混合模型的縱向應力分布斜率大于全鋼模型，但縱向應力變化范圍小于全鋼模型.原因是上層建筑采用普通鋼時其剛度較大，承擔的載荷較大.因而在駕駛甲板和陽光甲板上，全鋼模型的縱向應力變化范圍較鋼鋁模型大.

全鋼模型和鋼鋁混合模型的陽光甲板至羅經甲板，以及羅經甲板上各構件的縱向應力分布較接近，甲板間也呈現出近似線性，但是在羅經甲板上為壓應力.原因在從陽光甲板至羅經甲板之間的上層建筑結構較為密集，結構較強，因而鋼鋁混合模型和全鋼模型的縱向應力水平較低且相差不大.并且由于羅經甲板與陽光甲板之間的上層建筑長度較小，羅經甲板與陽光甲板及之下的船體之間的相互作用致使羅經甲板出現反向彎曲，在甲板處出現壓應力.

3.2結構應力與變形分析比較

由圖5可見，全鋼模型中無論是單元最大應力還是整體應力水平都較鋼鋁混合模型大.原因是全鋼上層建筑的剛度較鋁合金上層建筑大，因而其承擔了更多的載荷，使得其整體應力水平都高于鋼鋁混合模型.由圖7可見，雖然鋼鋁混合模型的最大應力值比全鋼模型的最大應力值稍小，但是由于鋼鋁混合模型中上層建筑承擔的載荷較全鋼模型小，因此全鋼模型的主船體整體應力分布優于鋼鋁混合模型.

由圖6可見，雖然鋼鋁混合模型的上層建筑承擔的載荷較全鋼模型小，其整體應力水平小于全鋼模型上層建筑相應應力的一半，但是由于鋁合金的彈性模量約為普通鋼的1/3，因此，鋼鋁混合模型的最大變形值和整體變形水平都較全鋼模型大.由圖8可見，鋼鋁混合模型的最大變形值和整體變形均較全鋼模型大.主要原因是全鋼模型的上層建筑承擔的外載荷較鋼鋁混合模型大，在相同外載荷作用下鋼鋁混合模型的主船體遭受的載荷更大，因而變形更大.

4結 束 語

超級游艇的上層建筑對主船體的結構響應有較大的影響，并且不同材料上層建筑的影響程度不同. 超級游艇設置的多層甲板產生了較強的局部作用力，導致其斷面應力分布特性不同于未設置多層甲板的常規船舶.甲板之間的船體外板的縱向應力基本上都是近似線性分布，符合船體梁理論.不同甲板間的船體外板由于局部作用力、材料不同和所處位置的綜合影響，導致其線性分布的斜率不同.并且，在各層甲板上由于局部作用力的影響，致使同一甲板上的點的縱向應力依賴于其橫向位置.由于普通鋼和鋁合金的物理屬性不同，導致超級游艇的主船體和上層建筑承擔的載荷不同，因此，鋼鋁混合模型與全鋼模型的典型剖面縱向應力分布也存在應力幅值和變化范圍不同以及分布斜率的差異.

當上層建筑分別采用鋁合金和普通鋼時，由于兩者的物理性質不同導致上層建筑抵抗外載荷的水平不同，其中全鋼模型的上層建筑擔更大的載荷.因此，全鋼模型的主船體和上層建筑的最大變形值和變形狀況都較鋼鋁混合模型小；應力最大值都較鋼鋁混合模型大.

參 考 文 獻

[1]SHIN J G, LEE J Y, LEE J H, et al. A Study on the structural design and structural analysis for small yacht [J]. Journal of the Society of Naval Architects of Korea,2006,43(1):75-86.

[2]ALAIMO A, MILAZZO A, TUMINO D. Modal and structural FEM analysis of a 50 ft pleasure yacht[J]. Applied Mechanics and Materials,2012,215:692-697.

[3]BOOTE D, PAIS T,CAMPORSE E.Numerical and experimental analysis of the dynamic behavior of large yacht superstructures[C].Design and Construction of Super & Mega Yachts, Genoa, Italy,2013.

[4]李堯.玻璃鋼雙體游艇結構強度的計算方法研究[D].廈門：集美大學,2014.

[5]LAMB T, BEAVERS N. The all-aluminum naval ship-The way to affordable naval ships[C]∥Proc. 10th International Naval Engineering Conference and Exhibition: HM Naval Base, Portsmouth, United Kingdom，2010.

[6]PEI Z, IIJIMA K, FUJIKUBO M, et al. Simulation on progressive collapse behavior of whole ship model under extreme waves using idealized structural unit method[J]. Marine Structures,2015,40:104-133.

[7]LARSSON L, ELIASSON R, ORYCH M. Principles of yacht design[M]. Wikipedia:A&C Black,2014.

Structural Strength Analysis of Super Yacht in the Waves

with Different Materials Applied in the Superstructure

ZHANG ShuilinPEI ZhiyongCAI WeiWU Weiguo

(SchoolofTransportation,WuhanUniversityofTechnology,Wuhan430063,China)

Abstract：The three dimensional singularity distribution method based on potential flow theory is adopted to calculate the motion of the super yacht and the corresponding external water pressure firstly under specified sea condition. Then, obtained loads are transformed to the whole structural model and structural strength analysis of the super yacht in waves is performed by FEM software with aluminum alloy and ordinary carbon steel used in the superstructures respectively. According to the results, the superstructure of super yacht has a great influence on the characteristics of stress distribution of a typical cross-section, and different superstructures have different effect. Besides, the multi-decks of super yacht will change the stress distribution. Furthermore, stress and deformation of super yacht’s structure are not the same when aluminum alloy and steel are taken for the superstructure respectively.

Key words：super yacht; superstructure; structural response; cross-section stress distribution

收稿日期：2015-11-18

doi:10.3963/j.issn.2095-3844.2016.01.035

中圖法分類號：U663.2