陳爐云 侯國華 張?jiān)7?/p>
1 上海交通大學(xué)船舶海洋與建筑工程學(xué)院,上海 200240 2 總裝備部陸軍裝備科訂部,北京 100034
船舶碰撞是一種在巨大沖擊載荷作用下的復(fù)雜非線性瞬態(tài)響應(yīng)過程,它具有非常明顯的動(dòng)力特性,在碰撞區(qū)內(nèi)的構(gòu)件一般都會(huì)迅速超越彈性階段而進(jìn)入塑性流動(dòng)狀態(tài),并出現(xiàn)撕裂、屈曲等形式的破壞和失效。船舶碰撞形式受很多因素影響,如碰撞時(shí)的環(huán)境條件,風(fēng)、浪、流等的影響,兩船相撞時(shí)的相對(duì)位置和速度,碰撞部位等等。
對(duì)于船舶結(jié)構(gòu)碰撞問題,國內(nèi)外開展了一系列的研究工作:Yagi等[1]對(duì)船舶尖艏結(jié)構(gòu)的碰撞特性進(jìn)行研究,對(duì)比分析不同形狀的艏部碰撞力,提出了緩沖艏部結(jié)構(gòu)設(shè)計(jì);Wang等[2]對(duì)船舶與橋墩碰撞的結(jié)構(gòu)變形進(jìn)行研究,分析了不同碰撞形式結(jié)構(gòu)失效模式,并對(duì)比了不同結(jié)構(gòu)形式的耐撞性;Tabri等[3]對(duì)船舶碰撞問題進(jìn)行了模擬實(shí)驗(yàn),并與數(shù)值分析進(jìn)行對(duì)比。胡志強(qiáng)等[4]綜述了船舶碰撞機(jī)理與耐撞性結(jié)構(gòu)設(shè)計(jì),對(duì)耐撞性船舶的研究方向進(jìn)行展望;Hu等[5]對(duì)海洋平臺(tái)的碰撞特性進(jìn)行系統(tǒng)的研究,提出了一種新的耐撞結(jié)構(gòu)形式;劉敬喜等[6]對(duì)單殼船舷側(cè)結(jié)構(gòu)的碰撞進(jìn)行分析,分析了耐撞性結(jié)構(gòu)形式;Hou等[7]在對(duì)結(jié)構(gòu)碰撞分析的基礎(chǔ)上,開展了基于多目標(biāo)優(yōu)化的結(jié)構(gòu)耐撞性分析;Klanac等[8]在對(duì)潛艇結(jié)構(gòu)碰撞分析的基礎(chǔ)上,開展了潛艇結(jié)構(gòu)的耐撞性分析;劉峰[9]基于耐撞性的理論,提出了一種新型船舶結(jié)構(gòu)形式,并進(jìn)行了數(shù)值驗(yàn)證;郭君等[10]采用理論和數(shù)值分析相結(jié)合的方式,研究了被撞船舷側(cè)典型帶板梁結(jié)構(gòu)的吸能和破壞機(jī)理,分析了殼單元和體單元及其網(wǎng)格尺寸對(duì)碰撞力和吸能結(jié)果的影響及帶板梁的吸能成分等問題;莊科挺等[11]提出了計(jì)及摩擦力影響后,船舶舷側(cè)加筋板耐撞性能分析的一種簡化分析方法,討論了球鼻艏撞擊作用下舷側(cè)加筋板的漸進(jìn)破壞過程。
船舶碰撞研究大多針對(duì)大型船舶,對(duì)于諸如巡邏艇、交通艇這樣噸位小、速度快的結(jié)構(gòu)碰撞研究較少。小型快艇碰撞特性與大型船舶有一定區(qū)別,開展小型快艇碰撞研究和耐撞性結(jié)構(gòu)設(shè)計(jì)具有重要意義。在這方面,Toyama[12]對(duì)小型快艇艏結(jié)構(gòu)碰撞特性進(jìn)行了研究;高嵩等[13]對(duì)小型艇舷側(cè)結(jié)構(gòu)碰撞特性進(jìn)行研究,指出了小型快艇與大型船舶碰撞特性的一些重要區(qū)別。
本文以小型快艇結(jié)構(gòu)耐撞性問題為研究對(duì)象,在分析一些與大型船舶結(jié)構(gòu)碰撞區(qū)別的基礎(chǔ)上,開展小型快艇結(jié)構(gòu)耐撞性設(shè)計(jì)研究,探索評(píng)價(jià)小型快艇結(jié)構(gòu)耐撞性衡準(zhǔn)問題,為中、小型快艇滿足特殊用途的設(shè)計(jì)建造提供技術(shù)支持。
耐撞性是衡量結(jié)構(gòu)抗撞性能好壞的標(biāo)準(zhǔn)。耐撞性設(shè)計(jì)是指采用合理的結(jié)構(gòu)設(shè)計(jì)來改善結(jié)構(gòu)的耐撞性,可用耐撞指標(biāo)(Crashworthiness Index)βc來表達(dá)。結(jié)構(gòu)耐撞性設(shè)計(jì)方法包括:采用新型耐撞結(jié)構(gòu)形式、應(yīng)用新型耐撞性材料。在結(jié)構(gòu)耐撞性優(yōu)化中,需建立相應(yīng)的耐撞性衡準(zhǔn),為結(jié)構(gòu)設(shè)計(jì)和優(yōu)化提供理論基礎(chǔ)和考核指標(biāo)。
在大型船舶結(jié)構(gòu)耐撞性分析中,采用結(jié)構(gòu)材料有效利用率的結(jié)構(gòu)質(zhì)量吸能能力作為結(jié)構(gòu)耐撞性評(píng)價(jià)指標(biāo),該指標(biāo)在表征大型船舶結(jié)構(gòu)的耐撞性時(shí)具有一定的說服力。
然而,不同于大型船舶,小型艇具有如下的碰撞特性:
1)被撞艇的損傷變形區(qū)域占全艇的比例較大;
2)被撞艇的碰撞損傷主要是艇體的總體彎曲變形,其艏艉結(jié)構(gòu)在船體結(jié)構(gòu)發(fā)生總彎曲變形時(shí)可以相互自由趨近;
3)碰撞后航態(tài)變化復(fù)雜且與撞擊速度密切有關(guān);
4)相互碰撞的小型艇都是變形體。
針對(duì)小型艇的碰撞特點(diǎn),如果僅將結(jié)構(gòu)吸能能力定義為結(jié)構(gòu)耐撞性指標(biāo),不能全面描述小型艇的碰撞特性。
將艇體結(jié)構(gòu)撞深(Collision Penetration)定義為小型艇結(jié)構(gòu)耐撞性的一個(gè)衡量標(biāo)準(zhǔn),寫為:
式中,D(x)為在確定的結(jié)構(gòu)破壞狀態(tài)下艇體結(jié)構(gòu)的撞深。用結(jié)構(gòu)撞深表征碰撞過程中結(jié)構(gòu)被壓潰的位移,它直觀地反映了被撞艇的結(jié)構(gòu)損傷程度。在結(jié)構(gòu)耐撞性設(shè)計(jì)中,追求結(jié)構(gòu)比較小的最大結(jié)構(gòu)撞深。
將被撞艇艇體結(jié)構(gòu)的最大塑性應(yīng)變(Plastic Strain)εp作為小型艇結(jié)構(gòu)耐撞性的衡量標(biāo)準(zhǔn),寫成:
式中,ε(x)為被撞艇艇體破壞時(shí)的結(jié)構(gòu)塑性應(yīng)變,表示艇體結(jié)構(gòu)材料的失效形式。小型快艇耐撞性設(shè)計(jì)目標(biāo)之一是追求結(jié)構(gòu)比較小的最大塑性應(yīng)變。
根據(jù)結(jié)構(gòu)撞深衡準(zhǔn)和結(jié)構(gòu)最大塑性應(yīng)變衡準(zhǔn),綜合考慮兩個(gè)衡準(zhǔn)參數(shù),建立小型艇結(jié)構(gòu)耐撞性歸一化評(píng)價(jià)模型:
式中,X=(x1,x2,…,xi)T為設(shè)計(jì)變量空間,包括結(jié)構(gòu)尺寸、拓?fù)湫螤畹葏?shù);為初始設(shè)計(jì)結(jié)構(gòu)最大撞深;為改進(jìn)設(shè)計(jì)后結(jié)構(gòu)最大撞深;為初始設(shè)計(jì)結(jié)構(gòu)最大塑性應(yīng)變;為改進(jìn)設(shè)計(jì)后結(jié)構(gòu)的最大塑性應(yīng)變;gi(x)為快艇結(jié)構(gòu)其它約束函數(shù),如結(jié)構(gòu)質(zhì)量、結(jié)構(gòu)強(qiáng)度等;和分別是設(shè)計(jì)變量xi的設(shè)計(jì)下限和設(shè)計(jì)上限;參數(shù)s為目標(biāo)函數(shù)的權(quán)重因子,并滿足s∈[0,1],當(dāng)s=1時(shí),表示只將結(jié)構(gòu)最大撞深指標(biāo)定義為優(yōu)化評(píng)價(jià)目標(biāo),s=0時(shí)表示只將結(jié)構(gòu)最大塑性應(yīng)變定義為優(yōu)化評(píng)價(jià)目標(biāo)。在文中,取s=0.5,可得
在小型快艇碰撞分析中,分析撞擊艇(Striking Ship)船艏與被撞艇(Stuck Ship)舷側(cè)結(jié)構(gòu)的正向垂直的碰撞情景。由于撞擊艇艇艏區(qū)域結(jié)構(gòu)剛度遠(yuǎn)大于被撞艇舷側(cè)結(jié)構(gòu),被撞艇的舷側(cè)被船艏撞擊是最危險(xiǎn)的形式,撞擊艇損失的動(dòng)能大部分由被撞艇的舷側(cè)結(jié)構(gòu)吸收,因而通常是被撞艇的舷側(cè)結(jié)構(gòu)發(fā)生破損。
選取兩艘型號(hào)相同、排水量約為6.7 t的快艇為研究對(duì)象。利用MSC.Dytran程序模擬一艘快艇撞擊另一艘快艇舷側(cè)的動(dòng)態(tài)過程,碰撞時(shí)的航速為36 kn。在碰撞計(jì)算中,撞擊艇和被撞艇同時(shí)當(dāng)作變形體來分析,被撞艇的大部分結(jié)構(gòu)均可能在碰撞過程中變形吸能。同時(shí),撞擊艇的速度較高、碰撞接觸面積所占比例較大。快艇碰撞有限元模型如圖1所示。
圖1 快艇碰撞有限元模型Fig.1 The collision finite element model of fast boats
對(duì)船體周圍水動(dòng)力的影響,采用附連水質(zhì)量法處理碰撞過程中的流體耦合作用,附連水質(zhì)量法是將相撞船舶周圍水的影響以船體附加質(zhì)量的形式加以考慮[14]。當(dāng)兩艘艇發(fā)生對(duì)中正撞時(shí),撞擊艇主要發(fā)生進(jìn)退運(yùn)動(dòng),被撞艇主要發(fā)生橫漂運(yùn)動(dòng),因此,只需計(jì)算橫漂運(yùn)動(dòng)、進(jìn)退運(yùn)動(dòng)的附連水質(zhì)量。附連水質(zhì)量采用Motora經(jīng)驗(yàn)公式確定,撞擊船附連水質(zhì)量為:
質(zhì)量賦給水線以下的舷側(cè)外板及船底板上。考慮到船舶的縱移或橫移等自由度的附連水質(zhì)量,采用改變結(jié)構(gòu)密度法來處理,將船殼板材料的密度加大到使船舶結(jié)構(gòu)的質(zhì)量等于船舶初始結(jié)構(gòu)與附連水質(zhì)量之和為止。
撞擊艇和被撞艇材料均為船用結(jié)構(gòu)鋼,其本構(gòu)模型為線性強(qiáng)化彈塑性材料,材料屈服模型選用雙線性模型。其材料力學(xué)參數(shù)如下:材料密度ρ=7 800 kg/m3,楊氏模量 E=2.1×1011Pa,泊松比μ=0.3,硬化模量 Eh=2.1×1011Pa,屈服應(yīng)力 [δ]=255 MPa。材料失效準(zhǔn)則采用最大塑性應(yīng)變作為材料的破壞準(zhǔn)則,結(jié)合本模型碰撞區(qū)域的網(wǎng)格大小,最大塑性應(yīng)變 εmax=0.2。材料應(yīng)變率敏感性按Cowper symonds本構(gòu)方程進(jìn)行處理。
當(dāng)兩小型艇發(fā)生側(cè)向?qū)χ写怪迸鲎矔r(shí),被撞艇的撞擊區(qū)域外板向內(nèi)凹陷比較嚴(yán)重,撞擊區(qū)域的外板和骨架應(yīng)力遠(yuǎn)高于屈服應(yīng)力,從提高結(jié)構(gòu)耐撞性角度考慮,應(yīng)加強(qiáng)撞擊區(qū)舷側(cè)結(jié)構(gòu)的橫向剛度。內(nèi)部骨架如底肋骨和舷側(cè)肋骨的連接處的應(yīng)力比較大,加強(qiáng)骨架橫向剛度是提高小艇耐撞性能的關(guān)鍵。
小艇對(duì)結(jié)構(gòu)重量比較敏感,且內(nèi)部空間有限,在進(jìn)行艇結(jié)構(gòu)耐撞性的優(yōu)化時(shí),要求在滿足船舶結(jié)構(gòu)強(qiáng)度的條件下,對(duì)遠(yuǎn)離碰撞區(qū)的結(jié)構(gòu)進(jìn)行適當(dāng)?shù)膬?yōu)化以控制重量;對(duì)碰撞區(qū)域的舷側(cè)結(jié)構(gòu)形式進(jìn)行優(yōu)化設(shè)計(jì),以加強(qiáng)橫向剛度。在此基礎(chǔ)上,修改3種改進(jìn)型舷側(cè)結(jié)構(gòu),運(yùn)用結(jié)構(gòu)耐撞性優(yōu)化理論,對(duì)這4種結(jié)構(gòu)進(jìn)行評(píng)估,分析最優(yōu)耐撞結(jié)構(gòu)形式。
在艇體結(jié)構(gòu)中,原型艇的結(jié)構(gòu)形式如圖2(a)所示;具有改進(jìn)結(jié)構(gòu)I的小型艇結(jié)構(gòu)形式如圖2(b)所示,其在#3~#10中的每個(gè)肋位新增一組豎直撐材;具有改進(jìn)結(jié)構(gòu)Ⅱ的小型艇的結(jié)構(gòu)形式如圖2(c)所示,其在#3~#10中的每個(gè)肋位新增一組上斜撐材;具有改進(jìn)結(jié)構(gòu)III的小型艇結(jié)構(gòu)形式如圖2(d)所示,其在#3~#10中的每個(gè)肋位新增上、下斜撐材各一組。
圖2 幾種舷側(cè)結(jié)構(gòu)形式Fig.2 Several types of side structures
在3種改進(jìn)型舷側(cè)結(jié)構(gòu)形式的設(shè)計(jì)中,同時(shí)對(duì)遠(yuǎn)離碰撞區(qū)船側(cè)結(jié)構(gòu)的參數(shù)進(jìn)行一定的修改,其目的是滿足被撞艇的在結(jié)構(gòu)改進(jìn)后總質(zhì)量不超過原型艇的要求。表1為原型艇和各個(gè)改進(jìn)艇的主尺度參數(shù)。
如表1中所示,修改后,各改進(jìn)舷側(cè)結(jié)構(gòu)的艇的質(zhì)量都比原型艇有所減少,滿足結(jié)構(gòu)重量不增加的約束條件。
表1 幾種改進(jìn)艇的主尺度參數(shù)Tab.1 Main dimensions of several types of new boats
在數(shù)值仿真計(jì)算中,假設(shè)碰撞時(shí)艇的吃水狀態(tài)是相同的,都處于正浮狀態(tài)。撞擊艇以36 kn的速度撞擊被撞艇的船舯部分。撞擊艇為原型艇,而被撞艇則分別用原型艇、改進(jìn)結(jié)構(gòu)I型艇、改進(jìn)結(jié)構(gòu)II型艇和改進(jìn)結(jié)構(gòu)III型艇。在利用有限元分析程序MSC/Dytran進(jìn)行數(shù)值仿真分析中,總共碰撞計(jì)算時(shí)間為 0.2 s,計(jì)算的步長為2.0×10-9s。
下面分別分析被撞艇結(jié)構(gòu)應(yīng)力應(yīng)變、最大塑性應(yīng)變以及撞深的計(jì)算結(jié)果。
圖3所示為撞擊艇以36 kn的撞擊速度撞向原型被撞艇的被撞艇結(jié)構(gòu)變形圖。由圖3可知,對(duì)于小型快艇結(jié)構(gòu),在碰撞后的結(jié)構(gòu)變形為整體性變形,船體結(jié)構(gòu)發(fā)生總體彎曲變形。結(jié)構(gòu)碰撞區(qū)占整個(gè)艇的比例比較大。
圖3 快艇碰撞變形云圖Fig.3 The collision strain of a fast boat
同時(shí),撞擊艇在與被撞艇碰撞后,由于碰撞力達(dá)到一定峰值后會(huì)出現(xiàn)卸載現(xiàn)象,在碰撞的后期階段,撞擊艇的艇艏上移,撞擊艇的艇艏和被撞艇的舷側(cè)已經(jīng)沒有接觸,向上飛起。此時(shí),撞擊艇的一部分已經(jīng)騎在被撞艇的結(jié)構(gòu)上,對(duì)被撞艇形成向下方向的壓力。
圖4給出了撞擊艇的撞擊速度為36 kn時(shí),原型結(jié)構(gòu)、改進(jìn)結(jié)構(gòu)I、改進(jìn)結(jié)構(gòu)II和改進(jìn)結(jié)構(gòu)III的被撞艇在最大撞深時(shí)舷側(cè)結(jié)構(gòu)的損傷變形情況。
圖4 幾種舷側(cè)結(jié)構(gòu)最大應(yīng)變狀態(tài)Fig.4 The maximum strain of several types of side structures
撞深表示船舶碰撞過程中被撞艇被壓潰的位移,反映了被撞艇的損傷程度,選擇最大撞深作為小型艇結(jié)構(gòu)耐撞性的衡量標(biāo)準(zhǔn)。
圖5所示為在碰撞過程中,不同舷側(cè)結(jié)構(gòu)的被撞艇在發(fā)生碰撞時(shí)的結(jié)構(gòu)撞深—時(shí)間曲線。由圖5可知,初始結(jié)構(gòu)具有最大的結(jié)構(gòu)撞深,而改進(jìn)結(jié)構(gòu)II具有最小的結(jié)構(gòu)撞深。在結(jié)構(gòu)最大撞深方面,改進(jìn)結(jié)構(gòu)I和改進(jìn)結(jié)構(gòu)III都比初始結(jié)構(gòu)有一定的改善。
圖5 不同結(jié)構(gòu)形式艇的最大撞深Fig.5 The maximum penetration of different side structures
根據(jù)對(duì)不同舷側(cè)結(jié)構(gòu)的被撞艇進(jìn)行數(shù)值分析,在碰撞區(qū)域內(nèi),部分結(jié)構(gòu)發(fā)生大的變形并超越了彈性范圍,結(jié)構(gòu)發(fā)生了破壞。
不同舷側(cè)結(jié)構(gòu)的快艇在同一碰撞速度的作用下,將被撞艇的最大結(jié)構(gòu)碰撞力和最大結(jié)構(gòu)撞深進(jìn)行歸一化處理,結(jié)合式(1)、式(2)和式(3),可得到不同舷側(cè)結(jié)構(gòu)形式的結(jié)構(gòu)耐撞性參數(shù),見表2。
表2 結(jié)構(gòu)耐撞性對(duì)比Tab.2 The comparison of boat crashworthiness
對(duì)表2進(jìn)行分析后可知,相比原結(jié)構(gòu),改進(jìn)結(jié)構(gòu)在以下方面均有提高:
1)撞深指標(biāo)方面,改進(jìn)結(jié)構(gòu)I,II和III的小型艇的耐撞性有所提高;
2)塑性應(yīng)變方面,具有改進(jìn)結(jié)構(gòu)I,II和III小型艇的耐撞性均有所提高,其中改進(jìn)結(jié)構(gòu)I的小型艇耐撞性提高最多。
綜合以上分析數(shù)據(jù),改進(jìn)結(jié)構(gòu)I型的艇耐撞性提高最多,是最佳的結(jié)構(gòu)形式,這主要是加撐材后船側(cè)板架結(jié)構(gòu)能有效地將碰撞外力傳遞給整個(gè)船體結(jié)構(gòu)進(jìn)而減少舷側(cè)結(jié)構(gòu)的損傷。
本文以小型快艇相互碰撞的舷側(cè)結(jié)構(gòu)損傷為研究對(duì)象,開展小型快艇的耐撞性評(píng)價(jià)模型研究,得出關(guān)于小型艇耐撞結(jié)構(gòu)設(shè)計(jì)的結(jié)論:
1)文中以行駛速度為36 kn的小艇為研究對(duì)象進(jìn)行碰撞分析,實(shí)際上對(duì)于不同的撞擊速度,其耐撞性評(píng)價(jià)模型的評(píng)價(jià)結(jié)果是有所不同的,因此在工程實(shí)際中,要綜合考慮各種典型的碰撞速度進(jìn)行相應(yīng)的處理。
2)在結(jié)構(gòu)最大撞深指標(biāo)和結(jié)構(gòu)塑性應(yīng)變指標(biāo)的歸一化處理過程中,對(duì)于各個(gè)指標(biāo)權(quán)重的處理還需進(jìn)一步研究,這與兩個(gè)艇的相對(duì)大小、撞擊速度和撞擊艇艏部的強(qiáng)度定義有關(guān),在這方面還有深入研究的空間。
[1]YAGI S,KUMAMOTO H,MURAGISHI O,et al.A study on collision buffer characteristic of sharp entrance angle bow structure[J].Marine Structures,2009,22(1):12-23.
[2]WANG L L,YANG L M,HUANG D J,et al.An impact dynamics analysis on a new crashworthy device against ship-bridge collision[J].International Journal of Impact Engineering,2008,35(8):895-904.
[3]TABRI K,MATUSIAK J,VARSTA P.Sloshing interaction in ship collisions-an experimental and numerical study[J].Ocean Engineering,2009,36(17/18):1366-1376.
[4]胡志強(qiáng),崔維成.船舶碰撞機(jī)理與耐撞性結(jié)構(gòu)設(shè)計(jì)研究綜述[J].船舶力學(xué),2005,9(2):131-142.HU Z Q,CUI W C.Review of the researches on the ship collision mechanisms and the structural designs against collision[J].Journal of Ship Mechanics,2005,9(2):131-142.
[5]HU Z Q,CUI W C.Study on collision characters for SPAR platform[J].China Ocean Engineering,2008,22(1):131-139.
[6]劉敬喜,葉文兵,胡紫劍.單殼船舷側(cè)結(jié)構(gòu)的碰撞分析[J].中國造船,2008,49(增刊):124-133.LIU J X,YE W B,HU Z J.The analysis of single-hull structure’s crashworthiness[J].Shipbuilding of China,2008,49(Supplement):124-133.
[7]HOU S J,LI Q,LONG S Y,et al.Multiobjective optimization of multi-cell sections for the crashworthiness design[J].International Journal of Impact Engineering,2008,35(11):1355-1367.
[8]KLANAC A,EHLERS S,JELOVICA J.Optimization of crashworthy marine structures[J].Marine Structures,2009,22(4):670-690.
[9]劉峰.基于耐撞性的新型船舶結(jié)構(gòu)形式研究[D].上海:上海交通大學(xué),2007.LIU F.Research on new types of ship structures based on crashworthiness[D].Shanghai:Shanghai Jiao Tong University,2007.
[10]郭君,朱楓,姚熊亮,等.被撞船舷側(cè)典型結(jié)構(gòu)單元的吸能機(jī)理分析[J].中國艦船研究,2011,6(1):12-17.GUO J,ZHU F,YAO X L,et al.Energy-absorbing analysis of typical structural elements on a collided ship broadside[J].Chinese Journal of Ship Research,2011,6(1):12-17.
[11]莊科挺,劉敬喜,劉元丹,等.船舶加筋板結(jié)構(gòu)耐撞性能分析[J].中國艦船研究,2011,6(3):16-20.ZHUANG K T,LIU J X,LIU Y D,et al.Analysis on the crashworthiness behavior of ship stiffened plate[J].Chinese Journal of Ship Research,2011,6(3):16-20.
[12]TOYAMA Y.Drift-wood collision load on bow structure of high-speed vessels[J].Marine Structures,2009,22(1):24-41.
[13]高嵩,陳爐云,張?jiān)7迹?小型艇舷側(cè)結(jié)構(gòu)碰撞損傷性能研究[J].中國造船,2011,52(1):112-118.GAO S,CHEN L Y,ZHANG Y F,et al.Study on damage performance of small ship's collision[J].Shipbuilding of China,2011,52(1):112-118.
[14]PEDERSEN P T,ZHANG S M.On impact mechanics in ship collisions[J].Marine Structures,1998,11(10):429-449.