• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    船舶碰撞機理三維解析法實現(xiàn)及恢復(fù)系數(shù)研究

    2017-01-02 13:56:29劉俊峰胡志強
    中國艦船研究 2017年2期
    關(guān)鍵詞:恢復(fù)系數(shù)動能坐標系

    劉俊峰,胡志強

    1上海交通大學(xué)海洋工程國家重點實驗室,上海200240

    2高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海200240

    船舶碰撞機理三維解析法實現(xiàn)及恢復(fù)系數(shù)研究

    劉俊峰1,2,胡志強1,2

    1上海交通大學(xué)海洋工程國家重點實驗室,上海200240

    2高新船舶與深海開發(fā)裝備協(xié)同創(chuàng)新中心,上海200240

    [目的]采用解析方法分析船舶碰撞動力特性較為快速和準確,其中外部動力學(xué)分析十分重要。[方法]為此,運用MATLAB程序?qū)崿F(xiàn)船舶碰撞外部機理三維簡化解析方法,計算兩艘船舶碰撞的動能損失,并與二維解析方法的計算結(jié)果進行比對。在實現(xiàn)船舶碰撞動能損失快速計算解析方法的基礎(chǔ)上,討論碰撞高度、角度和位置對動能損失的影響。此外,還研究碰撞場景對保守恢復(fù)系數(shù)的影響和保守恢復(fù)系數(shù)對動能損失的影響。[結(jié)果]結(jié)果表明:三維解析方法得到的動能損失小于二維解析方法,碰撞高度對于動能損失有明顯的影響;簡化解析方法中,對于碰撞角度大于90°的場景,恢復(fù)系數(shù)簡單地取0并不安全。[結(jié)論]在今后的外部動力學(xué)分析中,為了使動能損失的計算值更加準確,可以使用三維解析方法代替二維解析方法。

    船舶碰撞;外部動力學(xué);三維解析法;能量耗散;恢復(fù)系數(shù)

    0 引 言

    船舶碰撞是一個高度非線性過程,為分析其碰撞過程,碰撞機理常被分為外部動力學(xué)和內(nèi)部動力學(xué)兩個獨立的過程來開展研究[1]。外部動力學(xué)主要利用剛體運動理論,分析撞擊船和被撞船的運動姿態(tài),以及碰撞過程中的能量損耗。內(nèi)部動力學(xué)主要利用彈塑性力學(xué)理論解決撞擊船和被撞船的結(jié)構(gòu)變形阻力、變形能耗散和結(jié)構(gòu)損傷等問題。

    船舶碰撞會產(chǎn)生災(zāi)難性的后果,因此很多學(xué)者開展了船舶碰撞研究。在外部動力學(xué)的研究中,Minorsky[2]開展了開創(chuàng)性研究,假設(shè)被撞船發(fā)生完全塑性變形,利用動量守恒定律,將周圍附連水的影響視為附加質(zhì)量且碰撞過程中不變化,從而估計船舶碰撞中的動能損失。Pedersen等[3]提出了船舶碰撞的沖擊動力學(xué)模型,通過積分各方向上的接觸力和相對位移,獲得各方向上的能量耗散和沖量變化值。但是該模型僅適用于二維平面空間,并且由于方程建立在總體坐標系下,比較復(fù)雜,較難推廣到三維空間。Stronge[4]提出了三維碰撞解決方案,但是該模型僅涉及碰撞物體的速度和加速度求解。在Stronge[4]的研究基礎(chǔ)之上,Liu等[5]提出了船舶碰撞的三維動力學(xué)模型,求得各方向上的能量損耗和沖量變化值。該研究方法的主要特點是所有的方程都建立在局部坐標系中,并且可以得到該局部坐標系下沿著每個坐標軸能量耗散值的封閉解。另一方面,部分學(xué)者為了真實地模擬船舶碰撞過程,從時域上來研究碰撞能量的耗散過程。Petersen[6]提出了一種時域上解析二維船舶碰撞外部動力學(xué)的方法,基于船舶的水平面運動方程,將被撞船舷側(cè)簡化為4個非線性的彈簧,采用時間上數(shù)值積分,得到船舶碰撞中船舶的運動響應(yīng)、動能損失和破損深度。Tabri等[7-9]針對船舶的碰撞過程,結(jié)合實船實驗和模型實驗數(shù)據(jù),以船舶操縱模型為基礎(chǔ),提出一種解析船舶碰撞中船舶運動的方法。Yu等[10]應(yīng)用顯式非線性有限元軟件LS_DYNA的計算模塊,通過子程序添加水動力以及通過LS_DYNA計算碰撞力,提出了一種耦合外部動力學(xué)和內(nèi)部動力學(xué)的方法。

    在內(nèi)部動力學(xué)的研究中,Lützen等[11]利用超單元法估算船舶碰撞中的變形阻力和能量耗散,考慮了舷側(cè)結(jié)構(gòu)和船艏形狀的影響。Kitamura[12]利用FEM討論了很多簡化解析方法中忽略的問題,包括船體梁橫向彎曲、等效失效應(yīng)變、被撞船的前進速度和碰撞角度。Haris等[13]提出一種簡化解析方法,能夠快速估算船舶碰撞中的損傷和能量耗散,并用LS_DYNA的數(shù)值結(jié)果進行了驗證。Sun等[14]提出的簡化方法則能夠分別估算舷側(cè)外板、骨材和桁材的變形阻力和能量耗散,從而得到整個舷側(cè)結(jié)構(gòu)的變形阻力和能量耗散。

    在外部動力學(xué)簡化解析方法中,恢復(fù)系數(shù)對碰撞能的大小有直接影響,但目前缺少詳細研究成果。Pedersen等[3]提出了快速估算二維船舶碰撞能量損失和碰撞沖量的方法,該方法有封閉解。在其研究中,恢復(fù)系數(shù) e取值范圍為0≤e≤1。對于完全塑性碰撞,恢復(fù)系數(shù)取0;對于完全彈性碰撞,恢復(fù)系數(shù)取1;但是,關(guān)于恢復(fù)系數(shù)的選取方法,并沒有給出詳細解釋。Liu等[5]提出了船舶碰撞外部動力學(xué)的三維模型,并將之用于計算船—冰碰撞的作用力,也沒有對恢復(fù)系數(shù)的選取進行深入討論,只是取0計算?;謴?fù)系數(shù)是外部動力學(xué)模型中很重要的參數(shù),決定了兩船之間的相對速度在碰撞前后的大小關(guān)系。一般認為,如果恢復(fù)系數(shù)取為0,則撞擊船碰撞之后不會被被撞船舷側(cè)反彈,此時應(yīng)變能最大,動能損失最多,偏于保守。如果取恢復(fù)系數(shù)為1,則撞擊船會被彈開,此時動能損失最少,偏于危險。在不了解恢復(fù)系數(shù)的情況下,可以簡單地取恢復(fù)系數(shù)為0。但實際上,取恢復(fù)系數(shù)為0并不是保守的,對于不同的碰撞場景,保守恢復(fù)系數(shù)(這里,定義為使動能損失最大的恢復(fù)系數(shù),且取值范圍0≤e≤1)不同。

    本文將用MATLAB程序?qū)崿F(xiàn)船舶碰撞外部動力學(xué)三維簡化解析方法,討論碰撞高度、角度和位置對碰撞能的影響。此外,還將討論碰撞場景對于恢復(fù)系數(shù)的影響,并將用保守恢復(fù)系數(shù)計算得出的新的碰撞能(即碰撞過程中的動能損失)與恢復(fù)系數(shù)取0所對應(yīng)的碰撞能進行比較。

    1 外部動力學(xué)三維解析方法

    1.1 坐標系統(tǒng)

    總體坐標系(固定坐標系)與局部坐標系的定義如圖1所示,總體坐標系為慣性系,原點位于船舶重心,采用右手坐標系,X軸沿船艏方向,Z軸向上。局部坐標系原點位于碰撞點C處,其3個相互垂直的單位向量分別記為n1,n2和n3。向量n1和n2位于碰撞切平面處,向量n3垂直于該切平面。為了后面闡述的方便,還引入連體坐標系,初始時與總體坐標系重合。

    1.2 船舶碰撞運動方程

    船舶碰撞的過程中,船舶被視為剛體,因此,可以使用牛頓—歐拉方程來描述剛體的運動。剛體一般運動的牛頓—歐拉方程為:

    式中:ω為剛體的在局部坐標系下的角速度向量;v為剛體在局部坐標系下重心處的速度向量;M為局部坐標系下的質(zhì)量矩陣;F為局部坐標系下的力向量;J為局部坐標系下的轉(zhuǎn)動慣量矩陣;G為力矩向量。

    在碰撞中,通過考慮質(zhì)量矩陣和轉(zhuǎn)動慣量矩陣中的附加質(zhì)量和附加慣量,來簡化水動力對碰撞過程的影響。

    式(1)和式(2)具有非線性,不方便求解。假定碰撞時間很短,則 Mω×v和ω×Jω這兩項為小量,得到線性式(3)和式(4):

    式中,r表示從重心到碰撞點處的位置矢量。

    使用符號′區(qū)分撞擊船和被撞船,帶′的物理量表示撞擊船的物理量,否則表示被撞船的物理量。碰撞過程中,假定碰撞點C處的碰撞力不變,兩船的沖量dPi和dP′i表達式為:

    則船舶的運動方程張量形式為:

    式(7)~式(10)中:重復(fù)出現(xiàn)的下標表示求和,εijk表示置換矩陣,當(dāng)其下標為順時針排列時值為+1,逆時針排列值為-1,有重復(fù)下標時值為 0; Mij,M′ij,Iij,I′ij分別表示兩艘船局部坐標系下的質(zhì)量矩陣和慣量矩陣;rj,r′j分別表示兩艘船從各自重心到碰撞點處的位置矢量;Vj,V′j表示重心處的速度。

    上面得到的式(7)~式(10)與Stronge[4]提出的三維碰撞模型一致,可以采取與Liu等[5]同樣的求解過程。

    2 算例分析

    采用MATLAB軟件編制程序,實現(xiàn)船舶碰撞外部動力學(xué)的三維簡化解析方法,得到船舶碰撞中的動能損失。

    2.1 碰撞場景

    本文借用文獻[15]的算例。兩艘供應(yīng)船向前的航行速度均為4.5 m/s,兩艘船的船長82.5 m,船寬18.8 m,吃水7.6 m,排水量4 000 t。X和Y方向上的附加質(zhì)量系數(shù)分別取為0.05和0.85,艏搖的附加轉(zhuǎn)動慣量系數(shù)取為1。慣性半徑取為對應(yīng)船長的0.25。為與算例的結(jié)果進行對比,兩艘船之間的摩擦力系數(shù) μ0取值與算例相同,取0.6。

    為分析不同碰撞位置對碰撞結(jié)果的影響,選取被撞船從船艏到船艉的24個點為碰撞點,在XYZ總體坐標系下位置坐標如表1所示,其中α為水線角,表征被撞船水線的曲率。撞擊船的碰撞位置選取為撞擊船船艏,如圖2所示。并且,為考慮不同碰撞角θ對碰撞結(jié)果的影響,θ值分別取為30°,60°,90°,120°和150°。

    在本文的算例中,選取上述的兩艘供應(yīng)船為分析對象,附加質(zhì)量、慣量系數(shù)和回轉(zhuǎn)半徑均通過經(jīng)驗公式[16]求得。在二維碰撞中,忽略了垂向的影響,即假定被撞船A和撞擊船B碰撞點垂向坐標Zc和Z′c均為0,因此計算結(jié)果中只有垂蕩、橫蕩和艏搖運動。三維碰撞要考慮碰撞點垂向高度的影響,碰撞點在被撞船A處的垂向坐標Zc變化范圍為0~1Rxx,Rxx為被撞船A橫搖運動的回轉(zhuǎn)半徑,該范圍是根據(jù)船舶重心的垂向坐標為吃水的一半這一假定得到的。計算中,假定被撞船A上碰撞點的垂向坐標 Zc分別為0,0.25Rxx,0.5Rxx,0.75Rxx和1Rxx,撞擊船B上碰撞點的垂向坐標 Z′c為0.5R′xx。被撞船A上碰撞點在總體坐標系下的坐標 Xc和Yc及撞擊船B上碰撞點的坐標 X′c和Y′c和二維碰撞場景下相同。兩艘船舶碰撞角度分別為30°,60°,90°,120°和150°。因此,總共有120個碰撞點位置,600個碰撞場景。

    2.2 結(jié)果分析

    計算求得600個碰撞場景下的碰撞能量損耗率(動能損失和初始動能之比),并繪制指定碰撞角度θ和垂向碰撞位置場景下的碰撞能量損耗率隨碰撞點沿船長方向位置變化的曲線(圖3)。

    由圖3(a)~圖3(e)可以看出,碰撞點高度對于碰撞能量損耗率有較大的影響。當(dāng)碰撞點的Xc和Yc坐標確定,碰撞角度給定時,碰撞能量損耗率隨著碰撞點高度Zc的增大而減少,特別是發(fā)生船舯碰撞時,該現(xiàn)象更加明顯。這是由于Zc越大誘導(dǎo)產(chǎn)生的橫搖運動越大,從而動能損失越小,即碰撞能量損耗率越大。相同的碰撞高度,碰撞點靠近船艏或船艉時,會引起強烈的艏搖和橫搖運動;當(dāng)碰撞點在船舯時,由碰撞高度導(dǎo)致的橫搖運動占主要部分,艏搖運動較小。因而,在船舯碰撞時,碰撞高度的影響較為明顯。圖3(a)~圖3(e)還給出了文獻[15]的二維計算結(jié)果。可以看出,三維解析方法得到的碰撞能量小于二維解析方法。當(dāng)碰撞高度等于0時,三維方法得到的結(jié)果與二維方法基本接近。

    從圖3(f)中可以看出,碰撞角度和碰撞點位置對于碰撞能量損耗率的影響也較為明顯。當(dāng)碰撞角度為0°~120°時,碰撞能量損耗率隨碰撞角度的增加而增加;當(dāng)碰撞角度為120°~150°,碰撞能量隨碰撞角度增加而降低。因為,碰撞角度增加到一定數(shù)值時,兩艘船發(fā)生相對運動,使得動能損失減少。此外,碰撞角度為120°和150°時,被撞船A艏部碰撞時碰撞能量明顯大于船艉碰撞所對應(yīng)的能量,并且碰撞角度越大,艏部碰撞時碰撞能量增加得越明顯。原因是被撞船A艏部的曲率不為0,使得碰撞接近于相向碰撞。

    3 恢復(fù)系數(shù)的影響

    恢復(fù)系數(shù)的大小可能與兩艘船各自的特性(舷側(cè)結(jié)構(gòu)、船長、船寬、吃水和船型等)和碰撞場景參數(shù)(被撞船速度、撞擊船速度、碰撞角度和碰撞位置)有關(guān)。本節(jié)將討論碰撞場景和兩艘船各自的特性對保守恢復(fù)系數(shù)的影響,以及保守恢復(fù)系數(shù)對碰撞能的影響。

    3.1 碰撞場景和兩艘船的特性對保守恢復(fù)系數(shù)的影響

    為研究碰撞場景對保守恢復(fù)系數(shù)的影響,保持其他量不變,分別調(diào)整被撞船速度、撞擊船速度、碰撞角度和碰撞位置,得到不同的動能損失—恢復(fù)系數(shù)曲線,如圖4~圖7所示。圖中的標注“x”所在的點表示曲線的最高點,其橫坐標對應(yīng)于保守恢復(fù)系數(shù)。圖的題目中,v_stru代表被撞船的速度,v_stri代表撞擊船的速度。

    從圖4可以看出,對于碰撞角小于90°的情況,在恢復(fù)系數(shù)等于0時,動能損失最大。因為恢復(fù)系數(shù)越小,則被撞船的內(nèi)能(即應(yīng)變能)恢復(fù)為動能的部分就越小,動能損失越大。對于碰撞角大于90°的情況,如碰撞角為120°和150°,在恢復(fù)系數(shù)等于0處,卻沒有取得最大值。這是由于動能損失由摩擦能(切向力做功)和內(nèi)能(法向力做功)構(gòu)成。為驗證這一解釋,分別計算60°和120°時相應(yīng)的摩擦能和內(nèi)能。

    從圖8和圖9可以看出,60°時,摩擦能明顯小于內(nèi)能,因為碰撞角小于90°,撞擊船速度方向與被撞船速度方向成銳角,沿著接觸面的切向力的相對滑動距離短,摩擦力做功小,故動能損失的變化與內(nèi)能的變化一致。120°時,摩擦能與內(nèi)能大小相當(dāng),甚至大于內(nèi)能,因為碰撞角大于90°,沿著接觸面的切向力的相對滑動距離大,摩擦力做功大,此時動能損失變化與內(nèi)能變化不一致。由圖5~圖7可知,碰撞位置、被撞船速度和撞擊船速度對于保守恢復(fù)系數(shù)的選取也有影響。

    除了碰撞場景,兩艘船各自的特性(船長、船寬、吃水、排水量和方形系數(shù)等)對于保守恢復(fù)系數(shù)的影響,可以采用類似的方法得到。經(jīng)過研究發(fā)現(xiàn),兩艘船各自的特性對于保守恢復(fù)系數(shù)影響很小,這是因為恢復(fù)系數(shù)的影響決定于摩擦能和內(nèi)能的相對大小,摩擦能主要與相對滑動距離有關(guān),而兩艘船各自的特性對滑動距離的影響可以忽略。

    3.2 保守恢復(fù)系數(shù)對碰撞能的影響

    為具體分析保守恢復(fù)系數(shù)對碰撞結(jié)果的影響,取Zxx=0.5Rxx,Z′c=0.5R′c,分別求得兩種航速情況下的碰撞能量損耗率。計算過程如圖10所示。計算程序分為主程序和子程序,主程序計算碰撞場景下的動能損失,子程序計算具體碰撞場景下的保守恢復(fù)系數(shù)?;謴?fù)系數(shù)取保守恢復(fù)系數(shù)時,需要調(diào)用子程序,子程序的計算過程如圖10所示,計算得到的結(jié)果見圖11和圖12。

    從圖11和圖12可以看出,對于碰撞角大于90°的碰撞場景,恢復(fù)系數(shù)取0顯得并不保守,或者說并不安全,特別是碰撞角越大,保守恢復(fù)系數(shù)和恢復(fù)系數(shù)取0所對應(yīng)的碰撞能損耗率差距越大。對于碰撞角小于90°的情況,兩種恢復(fù)系數(shù)下的碰撞能損耗率曲線基本接近。圖12中,60°對應(yīng)的曲線在船艏出現(xiàn)上升,因為撞擊船相對于被撞船的相對速度矢量方向正好與船艏水線的切線方向形成小角度,從而有更大的相對滑動,導(dǎo)致摩擦能較大,碰撞能量損耗增加。

    表2和表3分別給出了圖11和圖12中各個碰撞場景下對應(yīng)的保守恢復(fù)系數(shù)??梢钥闯?,當(dāng)碰撞角度接近90°時,大多數(shù)保守恢復(fù)系數(shù)靠近0;150°左右時,大多數(shù)保守恢復(fù)系數(shù)靠近1。

    4 結(jié) 論

    本文給出了一套三維船舶碰撞解析計算方法和程序。相對于二維船舶碰撞解析方法,本文提出的三維船舶碰撞解析方法由于考慮了橫搖、縱搖及垂蕩運動,因此計算得到的碰撞能量小于二維船舶碰撞的結(jié)果。

    通過具體分析不同碰撞場景下的計算結(jié)果,得出結(jié)論如下:

    1)碰撞高度對于碰撞能量有較明顯的影響,發(fā)生船舯碰撞時,影響尤其明顯。

    2)碰撞點位置和碰撞角度對碰撞結(jié)果有較明顯的影響,尤其是發(fā)生相向碰撞時,船艏碰撞會產(chǎn)生更大的動能損失。

    3)對于碰撞角度大于90°的場景,恢復(fù)系數(shù)簡單地取0是不安全的。本文給出的程序,能夠針對每個具體的碰撞場景計算得到對應(yīng)的保守恢復(fù)系數(shù)。

    [1] 王自力,顧永寧.船舶碰撞運動的滯后特性[J].中國造船,2001,42(2):56-62.

    WANG Z L,GU Y N.Motion lag of struck ship in colli?sion[J].Shipbuilding of China,2001,42(2):56-64(in Chinese).

    [2] MINORSKY V U.An analysis of ship collisions with reference to protection of nuclear power plants[R]. New York:Sharp(George G.)Inc.,1958.

    [3] PEDERSEN P T,ZHANG S M.On impact mechanics in ship collisions[J].Marine Structures,1998,11(10):429-449.

    [4] STRONGE W J.Impact mechanics[M].Cambridge:Cambridge University Press,2004.

    [5] LIU Z H,AMDAHL J.A new formulation of the im?pact mechanics of ship collisions and its application to a ship-iceberg collision[J].Marine Structures,2010,23(3):360-384.

    [6] PETERSEN M J.Dynamics of ship collisions[J].Ocean Engineering,1982,9(4):295-329.

    [7] TABRI K,BROEKHUIJSEN J,MATUSIAK J,et al. Analyticalmodelling ofship collision based on full-scale experiments[J].Marine Structures,2009,22(1):42-61.

    [8] TABRI K,M??TT?NEN J,RANTA J.Model-scale experiments of symmetric ship collisions[J].Journal of Marine Science and Technology,2008,13(1):71-84.

    [9] TABRI K,VARSTA P,MATUSIAK J.Numerical and experimental motion simulations of nonsymmetric ship collisions[J].Journal of Marine Science and Technolo?gy,2010,15(1):87-101.

    [10] YU Z L,AMDAHL J R,STORHEIM M.A new ap?proach for coupling external dynamics and internal mechanics in ship collisions[J].Marine Structures,2016,45:110-132.

    [11] LüTZEN M,PEDERSEN P T.Ship collision damage[D].Denmark:Technical University of Denmark,2002.

    [12] KITAMURA O.FEM approach to the simulation of collision and grounding damage[J].Marine Struc?tures,2002,15(4/5):403-428.

    [13] HARIS S,AMDAHL J.Analysis of ship-ship colli?sion damage accounting for bow and side deformation interaction[J].Marine Structures,2013,32:18-48.

    [14] SUN B,HU Z Q,WANG G.An analytical method for predicting the ship side structure response in raked bow collisions[J].Marine Structures,2015,41:288-311.

    [15] ZHANG S M.The mechanics of ship collisions[D]. Denmark:Technical University of Denmark,1999.

    [16] POPOV Y N,F(xiàn)ADDEEV O V,KHEISIN D E,et al. Strength of ships sailing in ice[R].DTIC Document,Washington,DC:Army Foreign Science&Technolo?gy Center,1969.

    3D analytical method for the external dynamics of ship collisions and investigation of the coefficient of restitution

    LIU Junfeng1,2,HU Zhiqiang1,2
    1 State Key Laboratory of Ocean Engineering,Shanghai Jiao Tong University,Shanghai 200240,China
    2 Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration,Shanghai 200240,China

    The analytical method for predicting the dynamic responses of a ship in a collision scenario features speed and accuracy,and the external dynamics constitute an important part.A 3D simplified analytical method is implemented by MATLAB and used to calculate the energy dissipation of ship-ship collisions.The results obtained by the proposed method are then compared with those of a 2D simplified analytical method.The total dissipated energy can be obtained through the proposed analytical method,and the influence of the collision heights,angles and locations on the dissipated energy is discussed on that basis.Furthermore,the effects of restitution on the conservative coefficients and the effects of conservative coefficients on energy dissipation are discussed.It is concluded that the proposed 3D analysis yields a lesser energy dissipation than that of the 2D analysis,and the collision height has a significant influence on the dissipated energy.In using the proposed simplified method,it is not safe to simplify the conservative coefficient as zero when the collision angle is greater than 90 degrees.In the future research, to get more accurate energy dissipation,it is a good way to adopt the 3D simplified analytical method instead of the 2D method.

    ship collisions;external dynamics;3D analytical method;energy dissipation;restitution coefficient

    U661.43

    A

    10.3969/j.issn.1673-3185.2017.02.011

    http://kns.cnki.net/kcms/detail/42.1755.TJ.20170313.1604.018.html

    劉俊峰,胡志強.船舶碰撞機理三維解析法實現(xiàn)及恢復(fù)系數(shù)研究[J].中國艦船研究,2017,12(2):84-91. LIU J F,HU Z Q.3D analytical method for the external dynamics of ship collisions and investigation of the coefficient of restitution[J].Chinese Journal of Ship Research,2017,12(2):84-91.

    2016-07-28 < class="emphasis_bold"> 網(wǎng)絡(luò)出版時間:

    時間:2017-3-13 16:04

    國家自然科學(xué)基金重點項目資助(51239007)

    劉俊峰,男,1993年生,碩士生。研究方向:船舶碰撞與擱淺。E-mail:liujunfeng@sjtu.edu.cn

    胡志強(通信作者),男,1975年生,博士,副教授。研究方向:船舶碰撞與擱淺,海洋可再生能源,海洋工程結(jié)構(gòu)物動力學(xué)性能。E-mail:zhqhu@sjtu.edu.cn

    期刊網(wǎng)址:www.ship-research.com

    猜你喜歡
    恢復(fù)系數(shù)動能坐標系
    剛體彈性碰撞中恢復(fù)系數(shù)的探討
    新動能,源自創(chuàng)新力
    金橋(2022年4期)2022-05-05 06:09:50
    利用恢復(fù)系數(shù)巧解碰撞問題
    為構(gòu)建開放創(chuàng)新新高地增添動能
    “金企對接”轉(zhuǎn)換舊動能
    商周刊(2018年15期)2018-07-27 01:41:16
    解密坐標系中的平移變換
    澎湃新動能
    商周刊(2018年11期)2018-06-13 03:41:52
    坐標系背后的故事
    基于重心坐標系的平面幾何證明的探討
    落石碰撞法向恢復(fù)系數(shù)的模型試驗研究
    日韩大片免费观看网站| 欧美xxxx黑人xx丫x性爽| 大话2 男鬼变身卡| 精品久久久久久久人妻蜜臀av| 国产91av在线免费观看| 最后的刺客免费高清国语| 婷婷色综合大香蕉| 波野结衣二区三区在线| 麻豆成人av视频| 国产又色又爽无遮挡免| 国产色爽女视频免费观看| 成人二区视频| 国产人妻一区二区三区在| 好男人在线观看高清免费视频| 久久国内精品自在自线图片| 噜噜噜噜噜久久久久久91| 高清av免费在线| 男男h啪啪无遮挡| 日韩一区二区三区影片| 另类亚洲欧美激情| 街头女战士在线观看网站| 亚洲最大成人av| 简卡轻食公司| 男女那种视频在线观看| 寂寞人妻少妇视频99o| 国产精品久久久久久精品电影| av卡一久久| 女的被弄到高潮叫床怎么办| 黄色日韩在线| 最新中文字幕久久久久| 嫩草影院入口| 国产精品国产三级国产专区5o| 精品久久久精品久久久| 亚洲性久久影院| 新久久久久国产一级毛片| 午夜精品国产一区二区电影 | 性色av一级| 国产一区有黄有色的免费视频| 国产成人午夜福利电影在线观看| 日日摸夜夜添夜夜添av毛片| 国产午夜精品一二区理论片| 99久久精品一区二区三区| videos熟女内射| 亚洲精品乱码久久久v下载方式| 成人亚洲精品av一区二区| 特级一级黄色大片| 国产精品福利在线免费观看| 91在线精品国自产拍蜜月| 国产精品伦人一区二区| 国产探花极品一区二区| 草草在线视频免费看| 天天躁日日操中文字幕| 丝瓜视频免费看黄片| 欧美激情久久久久久爽电影| 最新中文字幕久久久久| 国产亚洲精品久久久com| 国产综合懂色| 肉色欧美久久久久久久蜜桃 | 国产欧美日韩精品一区二区| 国内精品美女久久久久久| 亚洲av免费在线观看| 插阴视频在线观看视频| 男人添女人高潮全过程视频| 亚洲精品视频女| 国产亚洲最大av| 亚洲成人久久爱视频| 秋霞伦理黄片| 国产亚洲5aaaaa淫片| 舔av片在线| 秋霞在线观看毛片| 18+在线观看网站| 婷婷色麻豆天堂久久| 久久精品综合一区二区三区| 青春草国产在线视频| 亚洲欧美成人精品一区二区| 精品人妻视频免费看| 亚洲精品久久午夜乱码| 男人狂女人下面高潮的视频| 老司机影院毛片| 七月丁香在线播放| 国产精品一区www在线观看| 日韩电影二区| 少妇人妻一区二区三区视频| 亚洲,一卡二卡三卡| 国产欧美日韩一区二区三区在线 | 亚洲成人中文字幕在线播放| 国产欧美另类精品又又久久亚洲欧美| 韩国av在线不卡| 亚洲欧美成人精品一区二区| 成人一区二区视频在线观看| 哪个播放器可以免费观看大片| 嫩草影院新地址| 18禁动态无遮挡网站| 亚洲最大成人中文| 免费不卡的大黄色大毛片视频在线观看| 少妇猛男粗大的猛烈进出视频 | 有码 亚洲区| 欧美xxxx性猛交bbbb| 国产又色又爽无遮挡免| 免费播放大片免费观看视频在线观看| 2021少妇久久久久久久久久久| 青春草视频在线免费观看| 男的添女的下面高潮视频| 国产国拍精品亚洲av在线观看| 久久精品夜色国产| 黄色配什么色好看| 啦啦啦在线观看免费高清www| 麻豆成人午夜福利视频| 亚洲无线观看免费| 久久久久久伊人网av| 热99国产精品久久久久久7| 男人添女人高潮全过程视频| 日本色播在线视频| 亚洲国产最新在线播放| 亚洲三级黄色毛片| 日本午夜av视频| 在线观看免费高清a一片| 亚洲国产精品国产精品| 91午夜精品亚洲一区二区三区| 在线观看国产h片| 欧美 日韩 精品 国产| 亚洲av二区三区四区| 18禁裸乳无遮挡动漫免费视频 | 97人妻精品一区二区三区麻豆| 校园人妻丝袜中文字幕| 在线亚洲精品国产二区图片欧美 | 51国产日韩欧美| 熟女av电影| 日本色播在线视频| 男人爽女人下面视频在线观看| 午夜免费观看性视频| 日日啪夜夜撸| 一个人观看的视频www高清免费观看| 久久精品国产a三级三级三级| 国产熟女欧美一区二区| 日韩强制内射视频| 欧美 日韩 精品 国产| 久久女婷五月综合色啪小说 | 欧美 日韩 精品 国产| 三级经典国产精品| 欧美国产精品一级二级三级 | 久久久久国产精品人妻一区二区| 波野结衣二区三区在线| 我要看日韩黄色一级片| 久久精品国产鲁丝片午夜精品| 国产免费福利视频在线观看| 日本爱情动作片www.在线观看| 直男gayav资源| 两个人的视频大全免费| 亚洲,一卡二卡三卡| 丝袜喷水一区| 国产免费一级a男人的天堂| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲精品乱码久久久v下载方式| 亚洲久久久久久中文字幕| 男人狂女人下面高潮的视频| 别揉我奶头 嗯啊视频| 看黄色毛片网站| 中文天堂在线官网| 亚洲一级一片aⅴ在线观看| 69人妻影院| 亚洲aⅴ乱码一区二区在线播放| 乱系列少妇在线播放| 久久久久网色| 国产成人aa在线观看| 中文字幕制服av| 亚洲天堂国产精品一区在线| 小蜜桃在线观看免费完整版高清| 国产伦在线观看视频一区| 成人免费观看视频高清| 亚洲精品乱码久久久久久按摩| 久久精品国产自在天天线| 午夜激情久久久久久久| 久久6这里有精品| 乱系列少妇在线播放| 老师上课跳d突然被开到最大视频| 精品少妇黑人巨大在线播放| 免费播放大片免费观看视频在线观看| 国产黄色视频一区二区在线观看| 黄色日韩在线| 国产综合精华液| 国产综合懂色| 看非洲黑人一级黄片| 大香蕉97超碰在线| 看免费成人av毛片| 精品一区在线观看国产| 中文字幕人妻熟人妻熟丝袜美| 午夜福利高清视频| 亚洲高清免费不卡视频| 中文资源天堂在线| 七月丁香在线播放| 中文字幕制服av| 亚洲av日韩在线播放| 欧美一区二区亚洲| 国产淫语在线视频| 日韩制服骚丝袜av| 热re99久久精品国产66热6| 亚洲av二区三区四区| 一级a做视频免费观看| 一级av片app| 制服丝袜香蕉在线| 91精品伊人久久大香线蕉| 日韩电影二区| 欧美bdsm另类| 国产精品99久久99久久久不卡 | 国产精品一区www在线观看| 人妻制服诱惑在线中文字幕| 高清av免费在线| 午夜精品国产一区二区电影 | 美女cb高潮喷水在线观看| 免费看不卡的av| 国产精品成人在线| 在线观看美女被高潮喷水网站| 精品亚洲乱码少妇综合久久| 制服丝袜香蕉在线| 久久人人爽人人爽人人片va| av在线app专区| 久久99蜜桃精品久久| 1000部很黄的大片| 在线观看免费高清a一片| 亚洲av中文字字幕乱码综合| 黑人高潮一二区| 干丝袜人妻中文字幕| 九九在线视频观看精品| 18禁在线播放成人免费| 80岁老熟妇乱子伦牲交| 草草在线视频免费看| 亚洲精品成人av观看孕妇| 成年版毛片免费区| 我要看日韩黄色一级片| 日韩一区二区视频免费看| 男女下面进入的视频免费午夜| 国产人妻一区二区三区在| 成人免费观看视频高清| 午夜免费鲁丝| 久久久精品免费免费高清| 国产亚洲5aaaaa淫片| 深爱激情五月婷婷| 日韩欧美精品v在线| av播播在线观看一区| 我要看日韩黄色一级片| 乱系列少妇在线播放| 国产日韩欧美亚洲二区| 精品少妇黑人巨大在线播放| 3wmmmm亚洲av在线观看| 国产黄片视频在线免费观看| 亚洲成人久久爱视频| 国产探花在线观看一区二区| 成年女人看的毛片在线观看| 精品久久久久久久久亚洲| 嫩草影院精品99| 人妻少妇偷人精品九色| 一级爰片在线观看| 人人妻人人澡人人爽人人夜夜| 欧美日韩精品成人综合77777| 日日摸夜夜添夜夜爱| 国产久久久一区二区三区| 国产高清不卡午夜福利| 又粗又硬又长又爽又黄的视频| 97精品久久久久久久久久精品| 亚洲三级黄色毛片| 九九久久精品国产亚洲av麻豆| 少妇熟女欧美另类| 你懂的网址亚洲精品在线观看| 亚洲色图综合在线观看| 内射极品少妇av片p| 又粗又硬又长又爽又黄的视频| 亚洲精品乱码久久久v下载方式| 久久99精品国语久久久| 日韩 亚洲 欧美在线| 插逼视频在线观看| 美女内射精品一级片tv| 少妇高潮的动态图| 一级二级三级毛片免费看| 精品久久久噜噜| 在线观看一区二区三区| 极品教师在线视频| 亚洲天堂av无毛| 麻豆成人午夜福利视频| 亚洲在线观看片| 成人亚洲欧美一区二区av| 看非洲黑人一级黄片| 97热精品久久久久久| 嘟嘟电影网在线观看| 亚洲内射少妇av| av国产久精品久网站免费入址| 老师上课跳d突然被开到最大视频| 亚洲人成网站在线观看播放| 成人国产麻豆网| 亚洲欧美日韩另类电影网站 | 别揉我奶头 嗯啊视频| 纵有疾风起免费观看全集完整版| 国产精品精品国产色婷婷| 国产成人aa在线观看| 最近最新中文字幕大全电影3| 国产爱豆传媒在线观看| 最近中文字幕高清免费大全6| 日韩视频在线欧美| 久久久午夜欧美精品| 国产精品爽爽va在线观看网站| 成人亚洲精品av一区二区| 国产熟女欧美一区二区| 18禁在线播放成人免费| 别揉我奶头 嗯啊视频| 国内精品美女久久久久久| 少妇裸体淫交视频免费看高清| 国产在线一区二区三区精| 午夜福利网站1000一区二区三区| 一级a做视频免费观看| 秋霞在线观看毛片| 看免费成人av毛片| 亚洲综合色惰| 最近的中文字幕免费完整| 十八禁网站网址无遮挡 | 18+在线观看网站| 七月丁香在线播放| 成人鲁丝片一二三区免费| 日本黄大片高清| 久久6这里有精品| 精品一区二区三卡| av播播在线观看一区| 亚洲欧美清纯卡通| 人妻一区二区av| 大码成人一级视频| 嘟嘟电影网在线观看| 一级毛片黄色毛片免费观看视频| 综合色av麻豆| 国产欧美日韩一区二区三区在线 | 国产精品无大码| 在线观看免费高清a一片| 日日啪夜夜爽| 免费av不卡在线播放| 国产成人精品婷婷| 亚洲成人一二三区av| 国产熟女欧美一区二区| 午夜激情久久久久久久| 51国产日韩欧美| 国产精品不卡视频一区二区| 中文字幕av成人在线电影| 久久久久久国产a免费观看| 国产乱来视频区| 日韩,欧美,国产一区二区三区| 七月丁香在线播放| 99热这里只有是精品在线观看| 熟女电影av网| 在线观看国产h片| 涩涩av久久男人的天堂| 男女边摸边吃奶| videossex国产| 精品久久久久久久人妻蜜臀av| 99久久精品热视频| 国产精品一二三区在线看| 午夜爱爱视频在线播放| 亚洲丝袜综合中文字幕| 黄片wwwwww| 一级黄片播放器| 精品视频人人做人人爽| 亚洲欧美日韩无卡精品| 欧美xxxx性猛交bbbb| 99re6热这里在线精品视频| 深爱激情五月婷婷| 又爽又黄无遮挡网站| 免费观看在线日韩| 偷拍熟女少妇极品色| 一级爰片在线观看| 五月伊人婷婷丁香| 日韩亚洲欧美综合| 国产欧美亚洲国产| 国产一区二区亚洲精品在线观看| 丝瓜视频免费看黄片| 午夜福利在线在线| 亚洲欧美清纯卡通| 不卡视频在线观看欧美| 日本猛色少妇xxxxx猛交久久| 99热这里只有精品一区| 午夜日本视频在线| 久久久久九九精品影院| 搞女人的毛片| 亚州av有码| 国产欧美日韩一区二区三区在线 | 日韩不卡一区二区三区视频在线| 久久精品夜色国产| 看十八女毛片水多多多| 一级毛片电影观看| 国产高清三级在线| 3wmmmm亚洲av在线观看| 亚洲欧美一区二区三区国产| 亚洲av成人精品一区久久| 2021天堂中文幕一二区在线观| 亚洲欧美清纯卡通| 免费看光身美女| 九色成人免费人妻av| 国产乱人偷精品视频| 国产精品福利在线免费观看| 在线 av 中文字幕| 成人亚洲精品一区在线观看 | 免费电影在线观看免费观看| 美女高潮的动态| 国内精品美女久久久久久| 亚洲色图综合在线观看| 午夜免费男女啪啪视频观看| 久久久久久伊人网av| 蜜桃亚洲精品一区二区三区| 国产女主播在线喷水免费视频网站| 毛片一级片免费看久久久久| 久久亚洲国产成人精品v| 成人国产av品久久久| www.色视频.com| 一级黄片播放器| 精品熟女少妇av免费看| 搡老乐熟女国产| 日日啪夜夜爽| 自拍偷自拍亚洲精品老妇| 在现免费观看毛片| 欧美高清成人免费视频www| 国产一区二区三区av在线| 免费观看的影片在线观看| 免费大片黄手机在线观看| 欧美精品国产亚洲| 免费看日本二区| 日韩一本色道免费dvd| 亚洲国产最新在线播放| 国产伦在线观看视频一区| 最新中文字幕久久久久| 秋霞在线观看毛片| 九九在线视频观看精品| 18禁在线播放成人免费| 亚洲三级黄色毛片| 街头女战士在线观看网站| 亚洲aⅴ乱码一区二区在线播放| 男女那种视频在线观看| 男女下面进入的视频免费午夜| 国产成人aa在线观看| 99视频精品全部免费 在线| 国产黄片视频在线免费观看| 国产在线男女| www.色视频.com| 晚上一个人看的免费电影| 久久热精品热| 日本一二三区视频观看| 成年版毛片免费区| 免费大片18禁| 亚洲av中文字字幕乱码综合| 一本色道久久久久久精品综合| 久久久久久久国产电影| 免费高清在线观看视频在线观看| 久久国内精品自在自线图片| 亚洲色图av天堂| 色视频在线一区二区三区| 波多野结衣巨乳人妻| av又黄又爽大尺度在线免费看| 特大巨黑吊av在线直播| 国产亚洲5aaaaa淫片| 午夜爱爱视频在线播放| av.在线天堂| 午夜福利视频精品| 久久精品国产亚洲av涩爱| 中文精品一卡2卡3卡4更新| 日韩一区二区三区影片| 亚洲av日韩在线播放| 能在线免费看毛片的网站| 深夜a级毛片| 嫩草影院入口| 人妻一区二区av| 大码成人一级视频| 又黄又爽又刺激的免费视频.| 免费av毛片视频| 午夜精品国产一区二区电影 | 亚洲aⅴ乱码一区二区在线播放| 91精品一卡2卡3卡4卡| 国产亚洲午夜精品一区二区久久 | 又粗又硬又长又爽又黄的视频| 国产高清三级在线| 99热网站在线观看| 少妇熟女欧美另类| 亚洲av日韩在线播放| 日韩免费高清中文字幕av| 免费大片黄手机在线观看| 欧美性感艳星| 欧美高清成人免费视频www| 国产亚洲午夜精品一区二区久久 | 免费av毛片视频| 直男gayav资源| 97人妻精品一区二区三区麻豆| 久久久久九九精品影院| 欧美高清成人免费视频www| 精品久久久久久电影网| 亚洲欧美中文字幕日韩二区| 97超碰精品成人国产| 亚洲美女视频黄频| 久热久热在线精品观看| 成人黄色视频免费在线看| 免费在线观看成人毛片| 精品人妻偷拍中文字幕| 男人爽女人下面视频在线观看| 欧美精品一区二区大全| 能在线免费看毛片的网站| 波野结衣二区三区在线| 久久6这里有精品| 国产成人freesex在线| 我要看日韩黄色一级片| 男人添女人高潮全过程视频| av在线天堂中文字幕| 国产成人精品福利久久| 一区二区三区乱码不卡18| 交换朋友夫妻互换小说| 国产亚洲av嫩草精品影院| 三级国产精品欧美在线观看| 久久99热6这里只有精品| 精华霜和精华液先用哪个| 新久久久久国产一级毛片| 国产av不卡久久| 久久久久久久久久成人| 97热精品久久久久久| 天天躁夜夜躁狠狠久久av| 欧美成人精品欧美一级黄| 国产亚洲精品久久久com| 欧美三级亚洲精品| 精品久久久精品久久久| 久久精品国产亚洲av涩爱| 肉色欧美久久久久久久蜜桃 | 日韩成人伦理影院| 日日撸夜夜添| 一区二区三区免费毛片| h日本视频在线播放| 国产黄片视频在线免费观看| 内射极品少妇av片p| 少妇人妻久久综合中文| 蜜桃亚洲精品一区二区三区| 久久99精品国语久久久| 日韩 亚洲 欧美在线| 国内少妇人妻偷人精品xxx网站| 男人爽女人下面视频在线观看| 边亲边吃奶的免费视频| 久久久色成人| 国产成人精品久久久久久| 性色avwww在线观看| 看十八女毛片水多多多| 一个人观看的视频www高清免费观看| 国产精品一区二区在线观看99| 欧美xxxx黑人xx丫x性爽| 亚洲欧美一区二区三区黑人 | 国产成人精品一,二区| 少妇人妻 视频| 亚洲精品乱码久久久v下载方式| 草草在线视频免费看| 亚洲精品乱码久久久v下载方式| 夫妻性生交免费视频一级片| 国产伦精品一区二区三区视频9| 免费看日本二区| 黄片无遮挡物在线观看| 免费看av在线观看网站| 亚洲欧美日韩无卡精品| 国产69精品久久久久777片| 亚洲欧美精品专区久久| 热re99久久精品国产66热6| 看非洲黑人一级黄片| 直男gayav资源| 亚洲av免费在线观看| 亚洲三级黄色毛片| av一本久久久久| 麻豆成人av视频| 亚洲国产精品成人综合色| 国产精品一区二区在线观看99| 亚洲精品国产av成人精品| 身体一侧抽搐| 51国产日韩欧美| 欧美日韩在线观看h| 免费黄网站久久成人精品| www.色视频.com| 日韩视频在线欧美| 综合色av麻豆| 亚洲在久久综合| 一个人观看的视频www高清免费观看| av天堂中文字幕网| 免费观看的影片在线观看| 亚洲欧美精品自产自拍| 亚洲成人久久爱视频| 亚洲精品,欧美精品| 久久久久久久大尺度免费视频| 免费av毛片视频| 成人一区二区视频在线观看| 国产精品爽爽va在线观看网站| 国产毛片a区久久久久| 在线免费十八禁| 黄片wwwwww| 久久久精品免费免费高清| 王馨瑶露胸无遮挡在线观看| 一级毛片黄色毛片免费观看视频| 免费观看的影片在线观看| 中文字幕av成人在线电影| 99热这里只有是精品在线观看| av播播在线观看一区| 九九爱精品视频在线观看| 免费不卡的大黄色大毛片视频在线观看| 联通29元200g的流量卡| 我要看日韩黄色一级片| 少妇高潮的动态图| 国产在视频线精品| 2022亚洲国产成人精品| 特大巨黑吊av在线直播| 亚洲av男天堂| 欧美潮喷喷水| 18禁裸乳无遮挡动漫免费视频 | 精品久久久久久久人妻蜜臀av| 精品一区二区免费观看| 免费人成在线观看视频色| 国产一区二区亚洲精品在线观看| 亚洲最大成人中文| 欧美成人一区二区免费高清观看| 中文字幕av成人在线电影| 精品99又大又爽又粗少妇毛片| 久久久国产一区二区| 成人综合一区亚洲| 亚洲精品国产成人久久av|