仲 鵬
(江陰賽福咨詢有限公司,江蘇 無錫 214400)
在船舶種類中,工程駁船是一種類型繁多,用途廣泛的船型。目前駁船在系泊或錨泊狀態下纜索的強度,通常是通過規范的經驗公式進行校核的。由于船型及工作環境的特殊性,這些經驗公式的計算結果往往不夠精確。本文將運用大型有限元軟件SESAM建立駁船的有限元模型,運用三維勢流理論分析其水動力特征,從而可以計算駁船所受到的環境載荷,并在此基礎上,分析駁船在不同的系泊和錨泊工況下,纜索受到的張力,校核系泊纜索的強度。
SESAM軟件是目前計算船舶與海洋結構物的環境載荷較為先進的軟件之一,它采用的理論為三維勢流理論,計算作用于駁船上的波浪載荷。
船體的外載荷矩陣可以表達成如下形式[1~4]:

式中:FV(t)為船體的外載荷,Ffk為佛汝德-克雷洛夫力,Fd為波浪繞射力,Fw為風力,FC為流力,Fsd為二階漂移力。

式(2)~式(6)中:Φ0和Φ7分別代表波浪的入射勢和繞射勢,P為風壓,Vc為流速,Aw和Ac分別為受風和受流面積。Tsd為二階傳遞函數,η(ωm)和 η(ωn)分別為2種頻率的波的幅值,i為次數,ω為頻率,ρ為密度,Rn為雷諾數。
采用 Wilson - theta[5]或者 Runge - Kutta[6]數值積分方法對其進行數值求解,可以得到作用于駁船上的外環境載荷。
駁船主要參數見表1。

表1 駁船主要參數
(1)碼頭頂靠滾裝的系泊布置
碼頭頂靠滾裝的系泊布置如圖1所示。

圖1 碼頭頂靠滾裝時系泊布置圖
(2)碼頭系泊布置
鋼絲繩碼頭系泊布置如圖2所示。

圖2 碼頭系泊布置
(3)航行中避風錨泊布置
航行中避風錨布置如圖3所示。

圖3 航行中避風拋錨布置
(4)大節段拋錨定位
除上述主要系泊布置外,在橋面的運輸和安裝過程中,亦有專門的錨泊布置方案。圖1~圖3可視為沿著來流(首向或者尾向來流),而大節段拋錨定位時則是側向來流。因此,在以下工況分析中,則分別計算駁船受到的首尾向和側向的環境載荷。圖4、圖5分別為大節段拋錨定位方案1、方案2。

圖4 大節段拋錨定位(方案1)

圖5 大節段拋錨定位(方案2)
計算過程中,取流速為3 m/s,風速為15 m/s,波浪取規則波,波高取0.5 m,周期4 s。由于創建工況時選擇的是規則波工況,而風載荷為隨機載荷,故先將風載荷與浪、流載荷分離開來分別計算。先計算流和浪的組合工況,最后加上風載荷。計算其首尾向載荷時采用流和浪與X軸成0°,以及側向載荷時流和浪與X軸成90°的情況。工況選擇時,主要考慮2種工況:2.5 m吃水以及3.5 m吃水。由于裝載了橋梁,故這2種受風面積相對較大,視為危險工況進行計算。
本文主要計算駁船的環境載荷,故只需建立駁船的外表面模型即可。另外全船左右對稱,只需建立一半的模型。取網格長度為2 m。駁船有限元模型、水動力模型分別如圖6、圖7所示。

圖6 駁船有限元模型

圖7 駁船水動力模型
由于軟件的限制,而不能獨立顯示駁船所受的首尾向和側向的環境載荷,但是可以得出作用于錨鏈上的拉力。本文將采用一種特殊方法進行處理:求駁船的首尾向波流載荷的時候,假設駁船通過一平行于X軸的系泊錨鏈連接于河底,流和浪與X軸成0°(沿首尾向),此時求得的系泊錨鏈靠近導纜孔處的拉力沿首尾方向的分力,即為駁船所受到的首尾向波流載荷。同理,求駁船的側向波流載荷的時候,假設駁船通過一平行于Y軸的系泊錨鏈連接于河底,流和浪與X軸成-90°(沿側向),此時求得的系泊錨鏈靠近導纜孔處的拉力沿側方向的分力,即為駁船所受到的側向波流載荷。計算時,系泊錨鏈的材料、尺寸、位置做簡單處理,因為它與所受的拉力無關。2種假設的分析模型分別見圖8和圖9。
圖10為工況1(3.5 m吃水)時,纜索在側向波流載荷作用下的張力時間歷程,圖11為工況1(3.5 m吃水)時,纜索在首尾向波流載荷作用下的張力時間歷程。圖12為工況2(2.5m吃水)時,纜索在側向波流載荷作用下的張力時間歷程;圖13為工況2(2.5 m吃水)時,纜索在首尾向波流載荷作用下的張力時間歷程。

圖8 首尾向受力分析模型

圖9 側向受力分析模型

圖10 工況1(3.5 m吃水)側向波流環境載荷
本文取7級風,平均風速v=15.69 m/s,風速v與風壓ω有如下關系:


圖11 工況1(3.5 m吃水)首尾向波流環境載荷

圖12 工況2(2.5 m吃水)側向波流環境載荷

圖13 工況2(2.5 m吃水)首尾向波流環境載荷
式中:γ為空氣容重,kN/m3;g為重力加速度,m/s2。在標準大氣壓情況下,設γ=0.012018 kN/m3,g=9.81 m/s2,則有

將風速v代入式(8),便可得出風壓。
設結構物的迎風面積為A,則結構的風載荷為F=Aω。
本文所研究的駁船3.5 m吃水和2.5 m吃水工況的受風面積分別為2595 m2和2782 m2,故駁船的風載荷曲線分別為399.3 kN和428 kN。
將上述各種工況下的風、浪、流載荷進行匯總,其結果見表2。

表2 環境總載荷計算結果
表2中,波浪載荷一項表示波浪對錨鏈張力的 作用,而實際的波浪載荷則需要乘上夾角的余弦值。由表2可以發現,工況2(2.5 m吃水)的側向載荷是最大的,為 630.3 kN。
根據圖1碼頭頂靠滾裝的系泊布置情況,認為船舶受到最為惡劣的工況2(2.5m吃水)的側向載荷的作用,此時的最大側向荷載F為630 kN。根據力平衡和彎矩平衡原理,建立力學模型,并對其進行求解,如圖14所示。其中T為鋼絲繩的張力,T1為系泊于5000 t級碼頭系纜墩上的鋼絲繩張力,T2為系泊于涂裝車間鋼絲繩的張力,T3為駁船尾部系泊于碼頭的鋼絲繩張力,T4和T5分別為碼頭作用于駁船尾部的支反力,T6為駁船首部連接于拖船的鋼絲繩的張力,可以得到:與涂裝車間連接的鋼絲繩的張力為376.5 kN;與碼頭相連的鋼絲繩的張力為330 kN;與定位船連接的鋼絲繩的張力為121.6 kN。

圖14 碼頭頂靠滾裝時纜索系泊受力分析模型
根據規范,取鋼絲繩的安全系數為2.5,則與涂裝車間連接的鋼絲繩要求能夠承受的張力不小于941.5 kN,與碼頭相連的鋼絲繩要求能夠承受的張力不小于825 kN,與定位船連接的鋼絲繩要求能夠承受的張力不小于302.5 kN。根據設計方提供的尺寸,T2和T3處為Φ46 mm鋼絲繩,其破斷負荷為1320 kN;T6處為Φ28 mm鋼絲繩,其破斷負荷大于400 kN。因此,這3條鋼絲繩的強度能滿足要求,在外載荷的作用下不會破斷。
根據圖2碼頭系泊布置情況,認為船舶受到工況2(2.5 m吃水)的側向風載荷的作用,此時的最大側向荷載為428.04 kN,首尾向受到波流載荷的作用力66.93 kN。若4根纜索橫向放置,根據力平衡和彎矩平衡原理,可以求得作用于這4根纜索上的張力為108 kN。取鋼絲繩的安全系數為2,則鋼絲繩要求能夠承受的張力不小于226 kN。根據設計方提供的尺寸,這4根纜索中有2根的直徑為Φ46 mm,另外2根的直徑為Φ22 mm,其破斷負荷均大于226 kN。因此,這4根鋼絲繩的強度能滿足要求,在外載荷的作用下不會破斷。
根據圖3航行中避風拋錨的系泊布置情況,認為船舶受到工況1(3.6 m吃水)的首尾向載荷的作用,此時的最大首尾向載荷為75.7 kN,而認為這一拋錨過程,船首始終正對風浪流。根據力平衡和彎矩平衡原理以及懸鏈線方程,可以求得:AM2Φ52 mm電焊有檔錨鏈的張力為40 kN。
錨鏈的安全系數為2,錨鏈要求能夠承受的張力不小于80 kN,而錨鏈的破斷負荷為1479 kN,故滿足強度要求,在外載荷的作用下不會破斷。
根據圖4航行中避風拋錨的系泊布置情況,認為船舶受到最為惡劣的工況2(2.5 m吃水)的側向載荷的作用,此時的最大側向載荷為630 kN。根據力平衡和彎矩平衡原理以及懸鏈線方程,建立力學模型,并對其進行求解,如圖15所示。其中,T1和T2分別為順流和逆流的錨鏈張力,T3為橋墩的支反力??梢郧蟮?AM2Φ52 mm電焊有檔錨鏈的張力為 577.5 kN。

圖15 大節段拋錨定位時纜索系泊受力分析模型(方案1)
錨鏈的安全系數為2,錨鏈要求能夠承受的張力不小于1155 kN,而錨鏈的破斷負荷為1479 kN,Φ46 mm鋼絲繩破斷負荷為1320 kN,故設計方案滿足強度要求。在外載荷的作用下錨鏈不會破斷。
根據圖5航行中避風拋錨的 系泊布置情況,船舶受到最為惡劣的工況2(2.5 m吃水)的側向載荷的作用,此時的最大側向荷載為630 kN。根據力平衡和彎矩平衡原理以及懸鏈線方程,建立力學模型,并對其進行求解,見圖16。其中,T1、T2、T3、T4分別為4根錨鏈的張力??梢郧蟮?AM2Φ52 mm電焊有檔錨鏈的張力為445.5 kN。錨鏈的安全系數為2,錨鏈要求能夠承受的張力不小于891 kN,而錨鏈的破斷負荷為1479 kN,Φ46 mm鋼絲繩破斷負荷為1320 kN,故設計方案滿足強度要求。在外載荷的作用下錨鏈不會破斷。

[1]肖越.系泊系統時域非線性計算分析[D].遼寧:大連理工大學,2005.
[2]戴仰山,沈進威,宋競正.船舶波浪載荷[M].北京:國防工業出版社,2007.
[3]Y M Low,R S Langley.Time and frequency domain coupled analysis of deep -water floating production systems[J].Applied Ocean Research,2006,(28):371 -385.
[4]陳新權.深海半潛式平臺初步設計中的若干關鍵問題研究[D].上海:上海交通大學,2007.