羊 衛(wèi)，趙 園，陸響輝

(中國船舶及海洋工程設(shè)計(jì)研究院，上海 200011)

0 引言

海上風(fēng)力發(fā)電作為可再生能源開發(fā)的重要方向之一，已經(jīng)成為深遠(yuǎn)海智能海洋牧場、超大型海上浮式基地建設(shè)和海上石油鉆探等海上重大項(xiàng)目重點(diǎn)研究的能源供給方式之一。目前，海上風(fēng)電已成為全球關(guān)注的焦點(diǎn)。我國沿海風(fēng)能資源開發(fā)潛力較大，為加快海上風(fēng)電的發(fā)展，部分沿海地區(qū)海上風(fēng)電場已投入運(yùn)營[1]。隨著海上風(fēng)電場運(yùn)營數(shù)量的不斷增加，專門用于其日常維護(hù)工作的運(yùn)維船需求量將出現(xiàn)逐年遞增的趨勢。

目前，國內(nèi)海上風(fēng)電運(yùn)維船以專業(yè)雙體運(yùn)維船最為適用。該型船具有甲板面積大、布置寬敞、穩(wěn)性好、吃水淺、航速較高，能承受較大風(fēng)浪能力等特點(diǎn)，可確保運(yùn)維工作順利實(shí)施和保障海上作業(yè)人員的安全，同時(shí)雙體船還具備操控性良好、使用可靠、維修方便等特點(diǎn)，可保證海上風(fēng)電運(yùn)維工作能高效完成并降低海上風(fēng)電廠的日常運(yùn)營成本[2]。本文研究的37 m雙體交通運(yùn)維船主要承擔(dān)海上風(fēng)電日常巡查、維護(hù)、交通等工作。該船結(jié)構(gòu)主體部分包括：兩個(gè)瘦長型的平行片體(單體船)以及抗扭箱和連接橋，后者布置在兩片體上部中間，將兩個(gè)片體通過焊接等方式牢固地連接在一起，起到“扁擔(dān)”效應(yīng)。相比于單體船，雙體船的船長型深比很小，船體具有足夠的總縱強(qiáng)度和剛度，因而其縱向承載能力具備一定的儲(chǔ)備。一般對(duì)于船長小于50 m的雙體船可以不進(jìn)行總縱強(qiáng)度校核[3]，但當(dāng)雙體船遭受90°浪向角和45°或135°斜浪時(shí)，兩片體之間的連接橋?qū)?huì)承受較大的總橫力矩和扭轉(zhuǎn)力矩載荷。因此，在結(jié)構(gòu)設(shè)計(jì)時(shí)，必須充分考慮整船尤其是連接橋結(jié)構(gòu)的總橫強(qiáng)度和扭轉(zhuǎn)強(qiáng)度，以保證雙體船總強(qiáng)度達(dá)到規(guī)范及設(shè)計(jì)的要求。

本文基于有限元分析方法，對(duì)雙體交通運(yùn)維船進(jìn)行總強(qiáng)度分析計(jì)算。根據(jù)規(guī)范[4]中相關(guān)要求，通過設(shè)計(jì)波法，采用SESAM水動(dòng)力軟件計(jì)算出雙體船總橫彎矩和扭轉(zhuǎn)力矩的長期預(yù)報(bào)值，計(jì)及其他相關(guān)載荷，對(duì)本船總橫強(qiáng)度和扭轉(zhuǎn)強(qiáng)度開展分析。通過分析本船整船的應(yīng)力分布情況，對(duì)雙體船結(jié)構(gòu)優(yōu)化提出建議，為船體輕量化和局部結(jié)構(gòu)加強(qiáng)提供有益參考。

1 主要參數(shù)及有限元模型建立

1.1 主要參數(shù)介紹

37 m雙體交通運(yùn)維船為雙機(jī)、雙槳、雙舵型雙體船，其參數(shù)為：垂線間長37.2 m，型寬10.4 m，片體寬3.2 m，型深3.8 m，設(shè)計(jì)吃水1.8 m，肋距0.5 m，梁拱0.1 m，設(shè)計(jì)航速15 kn，航區(qū)為近海航區(qū)。本船主要建造材料為CCSA鋼，采用橫骨架式建造，兩單船體之間的抗扭箱和連接橋結(jié)構(gòu)也采用橫骨架式。

CCSA鋼參數(shù)為：密度7 850 kg/m3，泊松比0.3，屈服強(qiáng)度235 MPa，彈性模量2.06×105MPa。

1.2 船體濕表面及整體結(jié)構(gòu)有限元模型

采用MSC.Patran完成整船有限元模型的構(gòu)建，其中采用三角形單元、四邊形單元來模擬甲板、舷側(cè)板、船底板、船體縱桁腹板、艙壁板、連接橋甲板、連接橋抗扭箱底板、連接橋橫隔板等殼板結(jié)構(gòu)，采用梁單元來模擬縱向骨材、加強(qiáng)筋、桁材面板等。37 m雙體交通運(yùn)維船有限元模型包括主甲板以下整船結(jié)構(gòu)。按右手坐標(biāo)系法則建立坐標(biāo)系，其原點(diǎn)位于兩片體艉封板連線中點(diǎn)基線處。整船有限元模型總節(jié)點(diǎn)個(gè)數(shù)為32 659，單元個(gè)數(shù)為46 719。整船有限元模型如圖1所示。

圖1 整船有限元模型

2 邊界條件及載荷工況

2.1 邊界條件

根據(jù)規(guī)范[4]，采用的邊界條件如圖2所示。在連接橋上甲板中縱剖面上取艏、艉部點(diǎn)A和點(diǎn)B，片體中部點(diǎn)C，對(duì)A點(diǎn)約束X、Y、Z三個(gè)位移分量，對(duì)B點(diǎn)的Y、Z位移分量進(jìn)行限定，對(duì)C點(diǎn)的Z向位移進(jìn)行約束。

圖2 總體模型施加的約束條件

2.2 計(jì)算工況及載荷施加

按照規(guī)范要求，本船計(jì)算工況主要分為兩個(gè)：滿載出水時(shí)總橫強(qiáng)度計(jì)算工況(波浪的浪向角為90°)和滿載出水時(shí)扭轉(zhuǎn)強(qiáng)度計(jì)算工況(波浪的浪向角為45°或135°)。每個(gè)計(jì)算工況包括的載荷有：船員及行李重量、舷外靜水壓力、艙內(nèi)淡水、燃油等液體壓力及總橫彎矩或扭轉(zhuǎn)彎矩。

2.2.1船員及行李重量

船員及行李重量按船員主要活動(dòng)區(qū)域，對(duì)該重量進(jìn)行均布近似處理。

2.2.2靜水壓力

舷外靜水壓力在船長方向以梯形分布，可根據(jù)艏、艉吃水位置來確定液面位置；在吃水垂直方向靜水壓力以三角形梯度變化。艙內(nèi)淡水、燃油等靜水壓力在吃水方向也呈線性變化規(guī)律，在有限元模型上加載時(shí)以壓力場形式完成液壓載荷的施加。

2.2.3總橫彎矩及橫向?qū)﹂_力計(jì)算

當(dāng)波浪為橫浪(浪向角為90°)時(shí)，雙體船將承受橫向彎曲載荷作用，在連接橋上將產(chǎn)生較大的橫向彎矩。根據(jù)規(guī)范[4]，連接橋的總橫彎矩可采取等效載荷施加的方式進(jìn)行，即通過沿片體橫向強(qiáng)框架的中縱剖線上的節(jié)點(diǎn)來施加等效后的橫向?qū)﹂_力，其中橫向?qū)﹂_力fy按式(1)進(jìn)行等效計(jì)算：

(1)

式中：Mbx為總橫彎矩；Z為橫向?qū)﹂_力作用點(diǎn)到連接橋中縱剖面中和軸的垂向距離；n1為包括橫艙壁在內(nèi)的橫向強(qiáng)框架個(gè)數(shù)之和。

2.2.4扭轉(zhuǎn)彎矩及反對(duì)稱垂向均布力計(jì)算

當(dāng)波浪為斜浪(浪向角為45°或135°)時(shí)，雙體船兩片體發(fā)生不同步縱搖，使船產(chǎn)生扭轉(zhuǎn)變形。根據(jù)規(guī)范[4]，雙體船扭轉(zhuǎn)力矩也可采取等效載荷來施加，即將扭轉(zhuǎn)力矩等效為反對(duì)稱分布在片體中縱剖面內(nèi)的垂向均布力，垂向均布等效力關(guān)于雙體船中縱剖面反對(duì)稱，也關(guān)于中橫剖面反對(duì)稱，等效的垂向分布載荷pz由式(2)計(jì)算：

(2)

式中：Mty為扭轉(zhuǎn)彎矩；L為雙體船船長；n2為片體中前部或中后部所需施加的節(jié)點(diǎn)個(gè)數(shù)之和。

2.3 總橫彎矩和扭轉(zhuǎn)彎矩等效載荷施加

總橫彎矩和扭轉(zhuǎn)彎矩等效加載的關(guān)鍵是計(jì)算確定Mbx和Mty的值，可根據(jù)文獻(xiàn)[5]提供的經(jīng)驗(yàn)公式計(jì)算確定，也可根據(jù)設(shè)計(jì)波法計(jì)算得到橫浪或斜浪下總橫彎矩或扭轉(zhuǎn)力矩的長期預(yù)報(bào)值。本文采用第二種方法完成式(1)和式(2)中的彎矩計(jì)算過程。

2.3.1設(shè)計(jì)波參數(shù)確定

根據(jù)規(guī)范[4]，波浪載荷可采用設(shè)計(jì)波法來計(jì)算，應(yīng)用SESAM水動(dòng)力軟件的HydroD(WADAM)模塊分別計(jì)算90°、45°和135°浪向下的波浪彎矩的RAO(傳遞函數(shù))。圖3為本船橫向彎矩和扭轉(zhuǎn)彎矩的傳遞函數(shù)。

獲取了橫向彎矩和扭轉(zhuǎn)彎矩載荷的RAO后，就可以對(duì)相應(yīng)的載荷進(jìn)行長期預(yù)報(bào)。根據(jù)挪威船級(jí)社提供的北大西洋波浪散布圖，利用相應(yīng)的RAO，依據(jù)P-M波浪譜來模擬散布圖，運(yùn)用二維韋布爾分布來擬合長期分布，最終對(duì)本船所受的兩個(gè)彎矩載荷分別進(jìn)行長期預(yù)報(bào)。應(yīng)用SESAM中的Postresp模塊，便可計(jì)算出相應(yīng)波浪載荷的長期預(yù)報(bào)值，其結(jié)果見表1。

圖3 按設(shè)計(jì)波計(jì)算的波浪彎矩的RAO

表1 長期預(yù)報(bào)所得結(jié)果

從表1分析可知，當(dāng)浪向角為45°時(shí)雙體船所受的扭轉(zhuǎn)力矩較大，確定橫向彎矩Mbx=3.08×106N·m，扭轉(zhuǎn)彎矩Mty=5.75×105N·m。

2.3.2等效載荷施加

根據(jù)所建立的整船模型，可確定式(1)中Z=3.95 m，n1=22，將其與橫向彎矩Mbx代入式(1)中，可求得對(duì)開力fy的值；將扭轉(zhuǎn)彎矩Mty、船長L、片體中前部節(jié)點(diǎn)個(gè)數(shù)77或片體中后部節(jié)點(diǎn)個(gè)數(shù)67代入式(2)中，可求得總橫扭矩作用下的等效反對(duì)稱垂向均布力。

3 計(jì)算結(jié)果匯總

將液壓載荷、90°橫浪及45°斜浪動(dòng)壓力產(chǎn)生的波浪彎矩和整船運(yùn)動(dòng)的慣性力施加到整個(gè)雙體船有限元模型上，運(yùn)用MSC.Nastran分析軟件便可求得整船應(yīng)力分布結(jié)果。

本文根據(jù)規(guī)范[4]確定本船所用CCSA鋼的許用應(yīng)力值，即等效應(yīng)力許用值為200 MPa，剪切應(yīng)力許用值為90 MPa。雙體交通運(yùn)維船結(jié)構(gòu)應(yīng)力結(jié)果見表2。本船整船結(jié)構(gòu)應(yīng)力云圖、整體變形圖分別如圖4、圖5所示。

圖4 整船等效應(yīng)力云圖

圖5 整船變形云圖

表2結(jié)構(gòu)應(yīng)力結(jié)果匯總表MPa

計(jì)算工況橫彎扭轉(zhuǎn)等效應(yīng)力剪切應(yīng)力等效應(yīng)力剪切應(yīng)力校核結(jié)果片體甲板90.918.3滿足片體外板65.921.1滿足片體縱向構(gòu)件62.732.514.86.95滿足片體橫艙壁184.077.136.814.10滿足片體橫框架161.072.324.810.10滿足連接橋甲板106.026.3滿足連接橋抗扭箱底板184.037.3滿足連接橋縱向構(gòu)件42.923.291.523.00滿足連接橋橫向構(gòu)件183.030.687.212.50滿足

4 結(jié)論

(1)根據(jù)表2中的應(yīng)力結(jié)果數(shù)據(jù)可知，本船的總橫強(qiáng)度和扭轉(zhuǎn)強(qiáng)度滿足規(guī)范要求。

(2)在斜浪45°浪向角作用下，雙體船主要受扭轉(zhuǎn)力矩的作用。根據(jù)整船應(yīng)力云圖和變形云圖可知，連接橋縱向構(gòu)件和橫向構(gòu)件是主要受力構(gòu)件。此外，連接橋與兩片體相交接附近區(qū)域產(chǎn)生的應(yīng)力響應(yīng)也較其他部位大。在橫浪即90°浪向角作用下，連接橋橫向構(gòu)件與片體橫向強(qiáng)構(gòu)件如橫艙壁和橫框架產(chǎn)生的應(yīng)力響應(yīng)較其他區(qū)域大，最大應(yīng)力達(dá)到了184 MPa。

(3)對(duì)比橫浪和45°斜浪作用下的應(yīng)力響應(yīng)結(jié)果可知，雙體船在橫浪工況下比較

危險(xiǎn)，即雙體船總橫強(qiáng)度在船體設(shè)計(jì)時(shí)應(yīng)作充分的考慮；而在斜浪作用下，雙體船受扭轉(zhuǎn)力矩的影響較小，說明本船具有較富余的扭轉(zhuǎn)強(qiáng)度。

(4)綜上分析可知：連接橋與片體相銜接附近區(qū)域的構(gòu)件如橫梁、梁肘板、抗扭箱橫構(gòu)件及靠近片體艏、艉部區(qū)域均受到較大應(yīng)力，且變形也較大，因此對(duì)于該區(qū)域的構(gòu)件應(yīng)采取加厚加強(qiáng)處理。而片體中靠近船中部分的構(gòu)件的應(yīng)力和變形較小，因此在船體結(jié)構(gòu)設(shè)計(jì)時(shí)可恰當(dāng)協(xié)調(diào)片體中部區(qū)域構(gòu)件的結(jié)構(gòu)尺寸，以滿足結(jié)構(gòu)輕量化的設(shè)計(jì)要求。

參考文獻(xiàn)：

[1] 楊程. 淺析海上風(fēng)電運(yùn)維船的發(fā)展[J].海峽科學(xué), 2016(12):78-80.

[2] 周華. 海上風(fēng)電運(yùn)維之風(fēng)電運(yùn)維船[J].風(fēng)能產(chǎn)業(yè)，2017(9)：36-39.

[3] 胡犇, 許晟, 梅國輝,等. 基于Patran的高速小水線面雙體船有限元結(jié)構(gòu)強(qiáng)度分析[J]. 艦船科學(xué)技術(shù), 2011, 33(10):46-49.

[4] 中國船級(jí)社. 國內(nèi)航行海船建造規(guī)范[M]. 北京: 人民交通出版社, 2015.

[5] 中國船級(jí)社．海上高速船入級(jí)與建造規(guī)范[M]．北京：人民交通出版社，2012.