王 雷，萬正權，黃進浩，劉俊杰

(中國船舶科學研究中心，江蘇 無錫214082)

0 引 言

隨著現(xiàn)代計算技術的發(fā)展，基于分析計算、符合第一力學原理的規(guī)范更適應現(xiàn)代船級社規(guī)范的發(fā)展趨勢。2006年4月，國際船級社協(xié)會頒布了油船共同結構規(guī)范和散貨船共同結構規(guī)范。這2 份規(guī)范的發(fā)布，改變了之前各家船級社對于規(guī)范標準各自制定和執(zhí)行的局面，抑制了因市場利益驅動而引發(fā)規(guī)范應用不統(tǒng)一導致的低標準、低安全水平船出現(xiàn)。共同結構規(guī)范與早期的傳統(tǒng)規(guī)范相比具有更高的安全性、合理性和先進性［1］。

但由于共同結構規(guī)范最初的擬定分開進行，某些條款的規(guī)定存在差異。同時，在共同結構規(guī)范的實際應用過程中，工業(yè)界對規(guī)范也提出了進一步完善和更新的建議和要求。為了解決上述問題，IACS針對油船和散貨船共同結構規(guī)范的不協(xié)調問題，通過專門工作組的協(xié)調工作，建立新框架結構下的協(xié)調共同規(guī)范 (HCSR)。由于需要滿足GBS 要求，HCSR 增加了很多以往規(guī)范中沒有的要求，如剩余強度、結構冗余度、人為因素等。在計算船體極限強度時若考慮擱淺、碰撞、腐蝕、磨損和疲勞等損傷因素，即為結構承載能力的剩余強度問題［2］。為了對破損船體的剩余強度進行計算和校核，DNV 提出了余度因子，同時ABS 規(guī)定了舷側和船底的破損尺寸。有學者［3］根據(jù)破損船體極限強度和外載荷的隨機性質，采用概率模型并以可靠性方法進行評估。

國際船級社協(xié)會一直運用結合失效模式的安全系數(shù)來進行結構設計，并在發(fā)布的共同結構規(guī)范中推薦Smith法計算船體極限強度。基于近幾年的船舶運營和結構計算分析，船舶工作者總結經(jīng)驗，在2013年4月推出了以草稿形式發(fā)行的協(xié)調共同結構規(guī)范，面向船舶領域公開征求意見以進一步完善后正式發(fā)行。2013 版協(xié)調共同結構規(guī)范對極限強度和船體梁載荷計算的諸多系數(shù)做出了新的修正，對散貨船和油船的強度校核有重要的影響。尤其在第五章船體梁強度中，新增針對剩余極限強度的計算和校核。本文將以某76 000 t 散貨船為算例，計算完整船體極限強度和碰撞剩余強度，進而對協(xié)調共同結構規(guī)范進行對比分析。

1 基于Smith法的極限強度計算方法

1.1 方法簡介

Smith法是計算船體極限彎矩的常用方法。通過逐漸增加船體縱向彎曲曲率的形式反映剖面構件的破壞過程，結合離散后有效單元的力學性能分析，得到船體梁的彎矩—曲率曲線。當曲線的斜率減小至0 時所對應的彎矩即為船體梁的極限彎矩。Smith法考慮構件的后屈曲特性，具有較好精度，使用簡單方便。

船體梁的失效由縱向結構單元的屈曲、屈服和極限強度綜合影響。計算船體梁極限彎曲能力時需要指定所有主要縱向結構單元的失效模式。每個結構單元的失效模式有彈性－塑性失效、梁柱屈曲、扭轉屈曲、折邊型材的腹板局部屈曲、扁鋼腹板的局部屈曲以及橫向加筋板格的屈曲6 種。結構分析時綜合考慮每種失效模式從而選取剖面結構單元的最弱失效模式。船體梁垂向極限彎曲承載能力的彎矩－曲率曲線由增量迭代法得到，在校核破損船體剩余強度時，破損部分的結構單元需要從橫剖面中去除。

1.2 單元劃分

船體梁橫剖面分為不同種類的可以承受總縱彎曲的單元。構成剖面的結構單元可以分為普通扶強材單元、加筋板單元和硬角單元，如圖1所示。

圖1 硬角單元劃分及尺寸示意Fig.1 Extension of the breadth of elements

圖2 船體剖面的單元構成Fig.2 The configuration of elements on a hull section

表1 單元的失效模式Tab.1 Models of failure of stiffened plate element and stiffener element

計算船體梁極限彎曲能力時需要指定所有主要縱向結構單元的失效模式。結構分析時綜合考慮各種模式，選取剖面結構單元的最弱失效模式。在校核破損船體剩余強度時，破損部分的結構單元需要從橫剖面中去除。

2 算 例

本文以某76 000 t 散貨船作為算例，通過船體梁載荷計算，得出散貨船受到的垂向彎矩載荷?；赟mith法，應用本文所編寫的程序計算船體極限強度，最后進行極限強度校核。本文選取的散貨船為BC－A 級。

2.1 船體梁極限載荷

2013年版的HCSR［4］中船體梁載荷部分規(guī)定:中拱狀態(tài)和中垂狀態(tài)下的垂向靜水彎矩分別由MSW－h(huán)和MSW－s表示，單位kN·m。

1)中拱工況:

2)中垂工況:

在任意船長位置上的垂向波浪彎矩采用下面的公式進行計算，單位kN·m:

1)中拱工況:

2)中垂工況:

對于散貨船，在完整航行情況下，選取HCSR規(guī)范中設計載荷工況A(S+ D)進行校核:

式中:γS為靜水彎矩部分的安全因子，取值1.0;γW為波浪彎矩部分的安全因子，取值1.2。

根據(jù)2013年版協(xié)調共同結構規(guī)范，算例中散貨船所受載荷的計算結果如表2所示。

表2 完整船體梁載荷(2013 版HCSR)Tab.2 Intact ship hull girder loads (HCSR in 2013)

根據(jù)共同結構規(guī)范(CSR)的船體梁載荷規(guī)定，散貨船所受載荷的計算結果如表3所示。

表3 完整船體梁載荷(2006 版CSR)Tab.3 Intact ship hull girder loads (HCSR in 2006)

和CSR 相比，2013 版HCSR的極限彎矩、靜水彎矩和中垂工況下的垂向波浪彎矩數(shù)值都有所減小，并且中拱工況下的極限彎矩大于中垂工況。

2.2 完整船體極限強度

本文應用Smith法，將散貨船的中剖面分為182個單元。經(jīng)過Fortran 語言編寫的程序計算，船體梁剖面總面積3 722 640 mm2，初始中和軸高度7 832.563 mm，剖面靜矩2.915 781 4E+10 mm3，剖面慣性矩2.323 272 3E+14 mm4，最大曲率5.856 804 3E－07。在中垂和中拱條件下的彎矩－曲率關系如圖3和圖4所示。中拱極限彎矩值為720 158 1kNm，中垂極限彎矩值為613 027 6kNm。

圖3 中拱工況下彎矩－曲率曲線Fig.3 Bending moment versus curvature curve in hogging condition

圖4 中垂工況下彎矩－曲率曲線Fig.4 Bending moment versus curvature curve in sagging condition

依據(jù)2013 版HCSR的相關規(guī)定，校核船體梁極限彎曲能力需要保證其滿足如下給出的校核標準:

1)需要在中拱和中垂2 種工況下校核垂向船體梁極限彎矩能力;

2)對于散貨船，設計載荷方案為靜水力載荷(S)+動水力載荷(D);

3)任一橫截面的垂向船體梁極限彎曲能力需要滿足下式:

式中:γR為船體梁垂向極限彎曲能力的安全系數(shù)，取為γR=γMγDB;γM為和材料、幾何、強度預報相關的安全系數(shù)，γM=1.1;γDB為和雙層底影響相關的安全系數(shù)，由于本文散貨船屬于BC-A 型，中拱工況，取γDB=1.25 中垂工況，γDB=1.0。

中拱情況的校核:

中垂情況的校核:

2 種工況均滿足HCSR 規(guī)定的極限強度校核條件。

相比于CSR 中極限強度校核規(guī)定的γR=1.1，2013年版HCSR的γR更為細致，考慮了材料、幾何、雙層底的影響，中拱工況和中垂工況的取值也各不相同，最后導致HCSR 對船體的垂向極限彎曲能力要求更高，船舶安全性更好。

2.3 破損船體剩余強度

船舶發(fā)生碰撞或擱淺事故后，船體具有的極限承載能力減弱，需要重新計算破損狀態(tài)下船體梁的承載能力以滿足剩余強度校核條件。本節(jié)以碰撞為例，考慮船舶碰撞后有效剩余剖面的非對稱性，以及碰撞后剖面的傾斜。應用Fortran 語言編寫程序，計算散貨船在破損狀態(tài)下的剩余強度，并依據(jù)2013版HCSR 新增的剩余強度部分進行校核。

2013 版HCSR 規(guī)定，破損橫剖面的碰撞評估中，破損位于底部的更為不利的橫剖面位置，破損尺寸和范圍如表4和圖5所示。

表4 碰撞破損范圍Tab.4 Damage extents for collision

圖5 2013 版HCSR 規(guī)定的碰撞破損范圍Fig.5 Damage extents for collision ruled by HCSR in 2013

在碰撞狀況下，船體梁的剩余強度校核需要考慮中拱和中垂工況下的垂向彎矩載荷MD依據(jù)下式計算:

式中:γSD=1.1 為碰撞狀態(tài)下靜水彎矩安全系數(shù);γWD=0.67 為碰撞狀態(tài)下垂向波浪彎矩安全系數(shù)。

表5 碰撞狀態(tài)船體梁載荷(HCSR)Tab.5 Hull girder loads for collision (HCSR)

在整船剖面模型的基礎上，按照2013 版HCSR規(guī)范，將破損部分的單元移除，不計入承載構件，得到剖面的有效單元數(shù)為169，跨距為2.6 m，材料的彈性模量E=206 000 N/mm2，屈服應力σs=320 MPa。經(jīng)過程序計算，船體梁剖面總面積3 364 395 mm2，初始中和軸高度7 076.091 mm，剖面靜矩2.380 676 7E10 mm3，剖面慣性矩2.076 5657E14 mm4，最大曲率6.482 927 7E－07。

圖6 中拱工況下彎矩－曲率曲線Fig.6 Bending moment versus curvature curve in hogging condition

圖7 中垂工況下彎矩－曲率曲線Fig.7 Bending moment versus curvature curve in sagging conditon

碰撞狀態(tài)下，船體梁在中拱和中垂條件下的彎矩－曲率關系如圖6～圖7所示。中拱剩余極限彎矩值為5 612 297 kNm，中垂剩余極限彎矩值為4 894 826 kNm。

在破損狀況下，船體梁垂向極限彎曲能力需要在中拱和中垂2 種工況下進行校核，設計載荷方案為靜水力載荷S+ 動水力載荷D，即考慮垂向靜水彎矩和垂向波浪彎矩。

在破損狀況下，船體梁任一橫剖面的垂向極限彎曲能力需要滿足以下判據(jù):

式中:γRD為破損狀況下，船體梁垂向極限彎矩的安全系數(shù)，γRD=1.0 ;CNA為中和軸系數(shù)，計算碰撞時CNA=1.1。

中拱工況下的校核:

中垂工況下的校核:

2 種工況均滿足HCSR 規(guī)定的極限強度校核條件。

3 碰撞破損模型對比

2013 版HCSR 規(guī)范中增加了剩余強度分析部分，其中關于破損模型做出了新的規(guī)定。結合本文完整船體極限強度和破損船體剩余強度計算，對比ABS和DNV的規(guī)定，最后檢驗2013 版HCSR的準確性。

在完整船體極限強度計算中，通過HCSR 給出的公式校核，完整船體的極限承載能力能夠抵抗環(huán)境施加的載荷，安全系數(shù)γM和γDB的選取合適、嚴格。在破損船體的船體梁載荷計算中，相對于其他船級社的規(guī)定，垂向波浪彎矩和靜水彎矩的安全系數(shù)取值較大，加之判據(jù)校核系數(shù)γRD和CNA分別取為1和1.1，致使HCSR 對碰撞后破損船體的剩余強度要求較高，可見2013 版HCSR 規(guī)定嚴格。

將2013 版HCSR 規(guī)范和ABS、DNV 給出的碰撞模型進行對比。

1)2013 版HCSR

HCSR 規(guī)定，破損橫剖面的碰撞評估中，破損位于底部的更為不利的橫向位置，尺寸如表5和圖5所示。

2)ABS

ABS［5］對船舶碰撞時船舷最危險的部分規(guī)定為離首柱后0.15 L 與離尾柱前0.2 L的一段船體。至少考慮2個剖面，一個是中剖面，一個是剪力最大的剖面。受撞部位假定在舷側板的上部，在強力甲板以下。這些結構必須整體或部分地從船體梁剖面模數(shù)計算中扣除。舷側板的垂向受損范圍h 定為4 m 或D/4 中的大者，D 為型深。強力甲板(包含縱向加強筋)的受損范圍為由舷側向內(nèi)h/2的部分，邊水艙的受損范圍如圖8所示。

圖8 ABS 規(guī)定的碰撞破損范圍Fig.8 Damage extents for collision ruled by ABS

3)DNV

DNV［6］定義的損傷比ABS 嚴重，損傷區(qū)域按最不利情況考慮。其中碰撞貫穿舷側的寬度為B/15，破損程度如表6所示。

表6 DNV 規(guī)定的碰撞破損范圍Tab.6 Damage extents for collision ruled by DNV

以76 000 t 散貨船為例，船寬B=32.26 m，3個規(guī)范的碰撞寬度和高度范圍如表7所示。

表7 三個規(guī)范的破損范圍對比Tab.7 Comparison of damage extents of the three rules

經(jīng)過對比可看出，HCSR 給出的碰撞范圍和DNV 相同。在破損范圍內(nèi)，包括甲板邊板、上甲板及其加強筋、舷頂列板、舷側板及其加強筋、邊水艙的斜板和隔板等構件在計算剖面模數(shù)時全部扣除。這種規(guī)定方式對碰撞時的總縱極限彎矩影響極大，校核最嚴格。而ABS 給出的破損高度在數(shù)值上比前兩者大大減小，破損寬度略大于前兩者。結合現(xiàn)實的海損事故案例，隨著航區(qū)船舶通行密度越來越大，船橋碰撞和船船碰撞事故急劇增多。HCSR 在破損范圍的規(guī)定能夠覆蓋絕大多數(shù)的相撞船舶，進而計算出準確的船舶碰撞時的總縱強度，HCSR 給出的船舶碰撞破損范圍適合。

4 結 語

國際船級社協(xié)會在2013年4月以草稿形式發(fā)行協(xié)調共同結構規(guī)范，面向船舶領域公開征求修改、驗證意見。2013 版HCSR 對極限強度和船體梁載荷計算的安全系數(shù)做出新的修正。在第五章船體梁強度中，新增加針對船體梁剩余強度的計算和校核。

通過本文的計算校核和分析，可得出以下結論:

1)2013 版HCSR 中新增的校核碰撞后剩余強度破損范圍選取適合，剩余強度安全性校核更嚴格;2013 版HCSR 中完整船體極限強度校核的中拱工況安全系數(shù)取值更高，極限強度校核公式更嚴格。

2)基于Smith法，應用Fortran 語言編寫極限強度計算程序，可以計及材料屈服、結構單元屈曲以及后屈曲的特性，很好地反映了船體剖面的實際破壞過程，繼續(xù)成為HCSR的推薦方法，計算結果準確，可為今后船體梁載荷和強度計算提供技術參考。

3)2013 版HCSR 打破了原來的規(guī)范體系和基于經(jīng)驗的傳統(tǒng)規(guī)范結構，理論分析型、技術背景充分且透明、內(nèi)容與相應軟件系統(tǒng)結合應用，將是未來規(guī)范發(fā)展模式。原各船級社規(guī)范體系中涉及工業(yè)界層面的標準在內(nèi)容的一致性、完備性等方面存在的問題需進一步研究。

［1］王剛，張道坤.IMO GBS 要求下的油船散貨船共同結構規(guī)范［C］//紀念徐秉漢院士船舶與海洋結構力學學術會議，2011:430－436.

［2］陳鐵云，王德禹，黃震球.船舶結構終極承載能力［M］.上海:上海交通大學出版社，2005.CHEN Tie-yun，WANG De-yu，HUANG Zhen-qiu.Ultimate carrying capacity of ship structures ［M］.Shanghai:Shanghai Jiaotong University Press，2005.

［3］祁恩榮，崔維成，彭興寧，等.破損船體非對稱彎曲及可靠性評估［J］.中國造船，2000，41(2):41－48.QI En-rong，CUI Wei-cheng，PENG Xing-ning，et al.Ultimate strength analysis and reliability assessment of damaged ship hulls in unsymmetrical bending ［J］.Shipbuilding of China，2000，41(2):41－48.

［4］IACS，Common structural rules for bulk carriers and oil tankers ［S］.London: International Association of Classification Society，April 2013.

［5］ABS，Rules for building and classing steel vessels［S］.New York:American Bureau of Shipping，2010.

［6］DNV，Hull structural design ships with length 100 meters and above［S］.Norway:Det Norske Veritas，2011.