周俊霖， 夏小浩， 邵漢東

(1.揚(yáng)帆集團(tuán)股份有限公司， 浙江 舟山 316100； 2.浙江嘉藍(lán)海洋電子有限公司， 浙江 舟山 316100)

0 引 言

由于穩(wěn)性好、舒適性好等眾多優(yōu)點(diǎn)，雙體船成為海上休閑觀光、海釣的理想船型。本文結(jié)合中國船級社《海上高速船入級與建造規(guī)范(2015)》[1]及《材料與焊接規(guī)范(2015)》[2]的要求，采用MSC.Patran/Nastran軟件，對34.7 m雙體休閑船應(yīng)用直接計(jì)算方法進(jìn)行船體總強(qiáng)度(包括總縱強(qiáng)度、橫向強(qiáng)度和扭轉(zhuǎn)強(qiáng)度) 計(jì)算，采用全船整體三維模型進(jìn)行總體結(jié)構(gòu)分析，計(jì)算載荷以及總縱彎矩、總橫彎矩和總扭矩。

1 雙體休閑船結(jié)構(gòu)強(qiáng)度計(jì)算

1.1 雙體休閑船相關(guān)參數(shù)和結(jié)構(gòu)型式

目標(biāo)船為鋼鋁復(fù)合結(jié)構(gòu)、單底、單甲板、前傾艏、方艉、艉縱傾、雙體折角線型、雙機(jī)、雙槳、雙舵、艉機(jī)型雙體船。主船體由連接橋結(jié)構(gòu)連接左、右兩個(gè)剛性水密片體組成的雙體結(jié)構(gòu)以及艏附加體組成。主甲板以上設(shè)置兩層甲板室。甲板室全部采用橫骨架式、鋁合金全焊接結(jié)構(gòu)。該雙體休閑船主要參數(shù)如表1所示。

表1 雙體休閑船主要參數(shù)

1.2 水動力載荷計(jì)算

采用SESAM軟件的Patran pre模塊建立目標(biāo)船外表面模型，原點(diǎn)位于艉封板與船體中線交點(diǎn)處，x軸指向船首為正方向，y軸指向左舷為正方向，z軸垂直向上為正方向。各裝載狀態(tài)下的實(shí)船質(zhì)量分布用沿船長方向分布的21根質(zhì)量棒及其兩端的質(zhì)量點(diǎn)模擬。在HydroD模塊中濕表面模型與質(zhì)量模型組合成水動力模型，如圖1所示。

圖1 水動力模型

全船設(shè)置21個(gè)橫向計(jì)算剖面和7個(gè)縱向計(jì)算剖面。橫向計(jì)算剖面(T1～T21)沿船長方向均勻分布?？v向計(jì)算剖面沿船寬方向，L1為中縱剖面，L2和L3分別為左右片體與連接橋連接處剖面，L4和L5分別為左右片體中心線位置剖面，L6和L7分別為左右片體與舷臺相交處剖面。波浪頻率按波長與船長比的0.2～3.0確定，步長取0.1。浪向角取0°～180°，間隔為15°。

基于三維勢流理論，應(yīng)用中國沿海波浪散布圖，用P-M波浪譜模擬散布圖中的海況，用二維Weibull分布擬合波浪長期分布，對波浪載荷進(jìn)行長期預(yù)報(bào)，得到船舶在航行壽命期中10-8超越概率下各計(jì)算截面處的波浪垂向彎矩、總橫彎矩和扭矩的長期預(yù)報(bào)值，如圖2～圖4所示。計(jì)算得到滿載出港工況下垂向彎矩長期預(yù)報(bào)值在浪向角為180°(迎浪)時(shí)最大，總橫彎矩長期預(yù)報(bào)值在浪向角為90°(橫浪)時(shí)最大，扭矩長期預(yù)報(bào)值在浪向角為75°(斜浪)時(shí)最大。

船體所受總縱彎矩(中拱彎矩和中垂彎矩)用靜水彎矩(由完整穩(wěn)性裝載計(jì)算書得到)加上波浪彎矩確定。

圖2 波浪垂向彎矩長期預(yù)報(bào)值

圖3 總橫彎矩長期預(yù)報(bào)值

圖4 扭矩長期預(yù)報(bào)值

1.3 全船有限元模型

建立雙體休閑船整船結(jié)構(gòu)三維有限元模型。采用板殼單元模擬船體外板、甲板、艙壁、肋板、隔板等板材以及強(qiáng)橫梁、強(qiáng)肋骨、桁材的腹板，采用梁單元模擬肋骨、縱骨、加強(qiáng)筋、扶強(qiáng)材、桁材的面板等[3-4]。網(wǎng)格大小按肋距和縱骨間距劃分。網(wǎng)格單元的邊長比須小于3。整船模型的單元數(shù)量為49 712個(gè)，節(jié)點(diǎn)數(shù)為34 112個(gè)。全船有限元模型如圖5所示。

圖5 全船有限元模型

坐標(biāo)系統(tǒng)采用右手笛卡爾坐標(biāo)系，原點(diǎn)O位于FR 0號肋位船底中線處，x軸從船尾指向船首為正方向，y軸從右舷指向左舷為正方向，z軸垂直向上為正方向。

(1) 邊界條件。給有限元模型施加邊界條件，約束6個(gè)位移分量以限制模型的空間剛體運(yùn)動，并且不能影響各部分結(jié)構(gòu)的相對變形[5-9]。在雙體釣魚船中縱剖面上取艏、艉部各1點(diǎn)A和B，中部舷側(cè)1點(diǎn)C。A點(diǎn)約束x、y、z方向的3個(gè)位移分量，B點(diǎn)約束y和z方向的2個(gè)位移分量，C點(diǎn)約束z方向的分量。邊界條件示例如圖6所示。

圖6 邊界條件示例

(2) 材料參數(shù) 。主船體及連接橋材料采用CCSA級船用鋼材。所有甲板室圍壁板材采用5083-H321鋁合金材質(zhì)，所有甲板室型材采用6082-T6鋁合金材質(zhì)。主船體及甲板室之間采用鋁鈦鋼三復(fù)合過渡接頭3A21+TA2+CCSA連接。材料參數(shù)：鋼的彈性模量取2.06×105N/mm2，鋁合金的彈性模量取0.70×105N/mm2；泊松比ν= 0.3。

1.4 載荷施加

1.4.1 總縱彎矩加載

《海上高速船入級與建造規(guī)范》(簡稱《規(guī)范》)假定船體總縱彎矩沿船長方向按正弦函數(shù)曲線形式分布，其表達(dá)式為

(1)

式中：x為從艉垂線量起的橫剖面縱向坐標(biāo)；L為船長。

函數(shù)曲線的幅值為船中橫剖面的總縱彎矩My,B值。

M(x)與施加沿船長分布的垂向力q(x)等效，q(x)(向上為正)計(jì)算式為

(2)

(3)

分別計(jì)算中拱與中垂兩種情況，在有限元模型片體甲板中線上施加沿船長分布的q(x)，施加到有限元模型上的節(jié)點(diǎn)力為q0/2。 在MSC.Patran軟件中模型的加載如圖7和圖8所示。

圖7 總縱彎矩加載(中拱)

圖8 總縱彎矩加載(中垂)

1.4.2 總橫彎矩加載

根據(jù)《規(guī)范》，用總橫彎矩Mx,B計(jì)算等效的橫向?qū)﹂_力Fy的公式為

(4)

式中：z為設(shè)計(jì)水線到連接橋中橫剖面中和軸的距離，m；d為設(shè)計(jì)吃水，m。

將橫向?qū)﹂_力施加于模型吃水的一半高度處，如圖9所示，分別按向外作用和向內(nèi)作用的兩個(gè)工況進(jìn)行計(jì)算。

圖9 橫向?qū)﹂_力加載示例

在實(shí)際計(jì)算時(shí)，將Fy換算為分布在連接橋整個(gè)長度范圍內(nèi)的分布載荷施加于船體模型，分布載荷q的計(jì)算式為

(5)

式中：Lb為連接橋縱向長度，m。

將分布載荷q換算成等效集中力P，作用于船體橫向強(qiáng)框架上。P計(jì)算式為

(6)

式中：S1和S2分別為橫向強(qiáng)框架的前后間距，m。

施加到模型上的節(jié)點(diǎn)力為P/2，在MSC.Patran軟件中模型加載如圖10和圖11所示。

圖10 橫向?qū)﹂_力加載(向外)

圖11 橫向?qū)﹂_力加載(向內(nèi))

1.4.3 總扭矩加載

根據(jù)《規(guī)范》，雙體船繞船寬方向y軸的總扭矩My,t以作用在片體半船長上反對稱分布的均布載荷p來等效。反對稱分布指的是：在同一片體上以中橫剖面為分界，前后載荷的方向相反，左右片體上載荷的方向也相反，如圖12所示。

圖12 等效均布載荷加載示例

等效均布載荷計(jì)算式為

(7)

使用集中力的形式加載時(shí)，其大小等于分布力與加載區(qū)間長度的乘積。對于該雙體船，加載區(qū)間即片體甲板的長度，該區(qū)間節(jié)點(diǎn)的數(shù)量為128個(gè)。在MSC.Patran軟件中模型加載如圖13所示。

圖13 總扭矩加載

2 屈服強(qiáng)度計(jì)算結(jié)果

2.1 許用應(yīng)力

雙體休閑船各構(gòu)件的最大計(jì)算應(yīng)力應(yīng)不大于表2所列的許用應(yīng)力。

表2 總強(qiáng)度許用應(yīng)力

本船主船體采用CCSA級鋼，取σSW=σS=235 MPa；上層建筑結(jié)構(gòu)全部采用鋁合金材料，板材用5083-H321，取σSW=215 MPa，型材用6082-T6，取σSW=250 MPa。

2.2 計(jì)算結(jié)果與分析

在雙體休閑船結(jié)構(gòu)總強(qiáng)度分析中，計(jì)算的載荷組合工況如表3所示。

各工況的整船等效應(yīng)力分布如圖14所示。

由圖14各工況下整船von Mises應(yīng)力計(jì)算分布可得主船體應(yīng)力分析結(jié)果如表4所示。主船體板單元最大等效應(yīng)力和最大剪應(yīng)力出現(xiàn)在工況LC 01和LC 02，即橫向?qū)﹂_力作用時(shí)；梁單元最大正應(yīng)力出現(xiàn)在工況LC 07和LC 08，工況LC 07和LC 08的應(yīng)力大于工況LC 09和LC 10的應(yīng)力。這說明：主船體最大應(yīng)力受總橫彎矩影響最大，其次為扭矩。在LC 03～LC 06工況時(shí)，主船體的應(yīng)力較小，說明總縱彎矩對該雙體休閑船的屈服總強(qiáng)度影響不大。

表3 計(jì)算工況

表4 主船體應(yīng)力分析結(jié)果匯總

關(guān)鍵構(gòu)件在最危險(xiǎn)工況下的應(yīng)力結(jié)果如表5所示。

表5 關(guān)鍵構(gòu)件應(yīng)力結(jié)果

圖14 整船等效應(yīng)力分布

圖15～圖20為最大應(yīng)力工況對應(yīng)的關(guān)鍵構(gòu)件應(yīng)力云圖。

圖15 上層建筑等效應(yīng)力(LC 01)

圖16 片體主甲板等效應(yīng)力(LC 07)

圖17 橫艙壁等效應(yīng)力(LC 01)

圖18 橫框架等效應(yīng)力(LC 01)

圖19 連接橋結(jié)構(gòu)等效應(yīng)力(LC 07)

圖20 船底龍骨等效應(yīng)力(LC 07)

3 結(jié) 論

對34.7 m雙體休閑船進(jìn)行了整船結(jié)構(gòu)總強(qiáng)度分析，計(jì)算結(jié)果表明：

(1) 對于雙體船的結(jié)構(gòu)總強(qiáng)度來說，總橫彎矩的影響最大，其次為扭矩，總縱彎矩影響較小。因此，在進(jìn)行雙體船船體結(jié)構(gòu)設(shè)計(jì)時(shí)，應(yīng)注意加強(qiáng)連接橋抗橫彎和扭轉(zhuǎn)的強(qiáng)度，可采取適當(dāng)增加連接橋上封閉式箱型結(jié)構(gòu)等措施。

(2) 由總強(qiáng)度計(jì)算結(jié)果發(fā)現(xiàn)，橫向構(gòu)件應(yīng)力較大，尤其是連接橋部位、強(qiáng)框和橫艙壁圓弧處。連接橋結(jié)構(gòu)的應(yīng)力分布與片體結(jié)構(gòu)的布置密切相關(guān), 連接橋中與片體橫艙壁相接的橫隔板的計(jì)算應(yīng)力相對較大，而與片體橫艙壁錯(cuò)開布置的橫隔板的計(jì)算應(yīng)力相對較小，說明連接橋結(jié)構(gòu)的應(yīng)力大小與在對應(yīng)處與其相接的片體結(jié)構(gòu)的剛度有很大關(guān)系。連接橋甲板的應(yīng)力相對較小而連接橋底板的應(yīng)力相對較大，這是因?yàn)檫B接橋甲板與片體甲板保持連續(xù)，而連接橋底板在橫向是間斷的。連接橋與片體橫艙壁連接處的甲板、底板以及橫隔板等有明顯的應(yīng)力集中現(xiàn)象，所以在結(jié)構(gòu)設(shè)計(jì)時(shí)應(yīng)特別考慮連接橋與片體相接處的結(jié)構(gòu)和連接方式。

(3) 扭矩作用下雙體船連接橋首尾部的應(yīng)力相對較大, 中間部分應(yīng)力相對較小。由于雙體船在扭矩作用下連接橋首尾部的主要是抗扭結(jié)構(gòu)，因此應(yīng)保證連接橋首尾部結(jié)構(gòu)有足夠的強(qiáng)度。

(4) 該雙體船的上層建筑尤其是第一層甲板室較大程度地參與了橫向和扭轉(zhuǎn)強(qiáng)度，能有效改善主船體的應(yīng)力狀態(tài)，對總縱強(qiáng)度的影響程度很小。

(5) 在橫向組合應(yīng)力作用下，機(jī)艙外龍骨也存在局部較大應(yīng)力，特別是在機(jī)艙區(qū)域，在設(shè)計(jì)建造時(shí)應(yīng)注意。