王培良, 吳曉芳, 張 婷
(1. 上海海事大學(xué),航運(yùn)仿真技術(shù)教育部工程研究中心,上海 201306; 2. 濰坊科技學(xué)院 汽車工程學(xué)院,山東 濰坊 262700; 3. 交通運(yùn)輸部水運(yùn)科學(xué)研究院,北京 100089;4. 山東交通職業(yè)學(xué)院 航海學(xué)院,山東 濰坊 262700)
集裝箱船舶關(guān)鍵位置的搖擺概率模型
王培良1,2, 吳曉芳3, 張 婷4
(1. 上海海事大學(xué),航運(yùn)仿真技術(shù)教育部工程研究中心,上海 201306; 2. 濰坊科技學(xué)院 汽車工程學(xué)院,山東 濰坊 262700; 3. 交通運(yùn)輸部水運(yùn)科學(xué)研究院,北京 100089;4. 山東交通職業(yè)學(xué)院 航海學(xué)院,山東 濰坊 262700)
為揭示船舶航行安全與船舶搖擺(尤其是橫搖)之間的內(nèi)在關(guān)系,針對集裝箱船舶關(guān)鍵位置在船舶正常航行過程中產(chǎn)生的搖擺構(gòu)建搖擺概率模型,討論和確定與之相關(guān)的參數(shù)。研究結(jié)果表明,集裝箱船舶不同位置的搖擺角度及加速度的概率模型均主要服從高斯分布,但其相關(guān)參數(shù)不同。該方法從集裝箱船舶搖擺的角度為船舶結(jié)構(gòu)安全和航行安全研究提供一定的理論依據(jù)。
集裝箱船舶;船舶搖擺;高斯分布;航行安全
為滿足世界航運(yùn)貿(mào)易發(fā)展的需求,集裝箱船舶正朝著高速化、大型化和多用途方向發(fā)展,隨之產(chǎn)生的船舶安全性問題備受關(guān)注。由船舶搖擺(尤其是橫搖)引發(fā)的船舶傾覆事件時有發(fā)生,造成的生命和財產(chǎn)損失巨大。[1]因此,對船舶安全航行中關(guān)鍵位置(如艏部、舯部等)的搖擺進(jìn)行建模,并求解關(guān)鍵位置在平衡狀態(tài)下的搖擺概率模型,對于檢測船舶航行的安全性和預(yù)測船舶橫搖而言是十分必要的。
將船舶在航行過程中某關(guān)鍵位置處的搖擺記為隨機(jī)變量V,包括縱搖角度、橫搖角度和垂直加速度等3個維度,分別記為VZ,VH和VA。VZ,VH和VA在船舶航行過程中的采樣值分別記為集合Z={xz1,…,xzi,…,xzn},H={xh1,…,xhi,…,xhn},A={xa1,…,xai,…,xan},其中采樣點(diǎn)個數(shù)為n。
高斯分布也稱為正態(tài)分布,自然界中有很多隨機(jī)現(xiàn)象都服從該分布。相關(guān)研究[2-3]證明,海浪波面起伏、船舶搖擺角度和波浪中船體的彎曲應(yīng)力等都服從高斯分布。因此,假定集裝箱船舶正常航行過程中不同位置處的搖擺(包括縱搖角度、橫搖角度及垂直加速度)變化值能被描述成高斯分布,且所服從分布的參數(shù)是不同的。[4]
設(shè)隨機(jī)變量X的密度函數(shù)為
(1)
式(1)中:μ和σ為常數(shù)(待估),且σ>0。稱X為服從參數(shù)μ和σ的正太分布,或高斯分布,記為X~N(μ,σ2),其相應(yīng)的分布函數(shù)為
(2)
首先假定船舶航行過程中的搖擺變化值符合高斯分布或威布爾分布,為驗證假定的可接受程度,需對樣本觀測數(shù)據(jù)進(jìn)行總體分布的粗略推斷。采用直方圖技術(shù)和q-q圖技術(shù)對其進(jìn)行驗證:首選采用q-q圖技術(shù)對船舶2個數(shù)據(jù)觀測點(diǎn)進(jìn)行驗證,判定其是否服從相同的分布,若服從相同的分布,則后續(xù)采用直方圖技術(shù)繪制其中一組樣本數(shù)據(jù)即可初步判定服從哪類分布,否則需分別判定2個觀測點(diǎn)的分布特征。
q-q圖是樣本分位數(shù)與指定分布的分位數(shù)的關(guān)系曲線圖,可用q-q圖技術(shù)來檢驗樣本數(shù)據(jù)是否服從指定的分布。利用MATLAB統(tǒng)計工具箱中的qqplot函數(shù)既可繪制1個樣本的q-q圖,檢驗樣本是否服從指定分布,還可繪制2個樣本的q-q圖,檢驗2個樣本是否服從相同的分布。[5]
在非參數(shù)統(tǒng)計領(lǐng)域,直方圖技術(shù)可在密度估計及數(shù)據(jù)分析方面發(fā)揮重要作用,能對總體分布進(jìn)行初步粗略估計,為研究總體的概率密度提供重要依據(jù)。[6]本文采用MATLAB自帶的hist函數(shù)對觀測點(diǎn)樣本值進(jìn)行縱向、橫向和垂向的數(shù)據(jù)統(tǒng)計處理,獲取各方向角度及加速度大小變化的概率密度函數(shù)的初步估計。[7]
初步確定概率密度函數(shù)之后,需對其參數(shù)進(jìn)行確定,本文采用最大似然估計法確定高斯分布或威布爾分布的位置參數(shù)。該方法的理論基礎(chǔ)是實際推斷原理——概率最大的事件最可能發(fā)生。
設(shè)總體X的密度函數(shù)為f(x;θ),其中θ=(θ1,θ2,…,θk),且(X1,X2,…,Xn)為總體X的一個樣本,則樣本的聯(lián)合密度函數(shù)為
(3)
若一次抽取得到的樣本觀察值為(y1,y2,…,yn),則樣本(X1,X2,…,Xn)在觀察值(y1,y2,…,yn)鄰域內(nèi)的概率為θ的似然函數(shù)可表示為
(4)
式(4)中:L(θ)為θ的函數(shù)。因此,要找到一個θ*使L(θ)達(dá)到最大,以此作為未知數(shù)θ的一個估計,該方法稱為最大似然估計法。通常利用導(dǎo)數(shù)求極限法來求L(θ)的最大點(diǎn),本文可采用該方法求得式(1)及式(2)中μ和σ在較高置信度下的最大似然估計,從而確定概率模型中的待估參數(shù)。
在確定模型參數(shù)之后,還需對船舶搖擺概率模型進(jìn)行檢驗,主要校驗樣本值與假定的概率模型之間的匹配程度。本文采用MATLAB自帶工具箱中的normplot和wblplot實現(xiàn),若樣本值符合假定的概率模型,則樣本值與某條概率密度曲線是線性重合的。[8]
船舶搖擺概率模型估計算法可描述為:
1)選定船舶結(jié)構(gòu)中某2處位置為搖擺測試點(diǎn),記錄這2個點(diǎn)在船舶航行過程中各方向上的搖擺變化值,分別記為
Z={xz1,…,xzi,…,xzn}
(5)
H={xh1,…,xhi,…,xhn}
(6)
A={xa1,…,xai,…,xan}
(7)
2)采用直方圖技術(shù)和q-q圖分布對式(5)~式(7)進(jìn)行統(tǒng)計,獲取隨機(jī)變量Z,H和A各自符合概率密度分布類型的初步估計。
3)對步驟2)推斷的概率模型進(jìn)行檢驗,確保所推斷的模型與理想模型之間的誤差在可接受范圍之內(nèi)。
4)基于步驟3)驗證所得的假定概率模型,采用最大似然估計法對概率密度分布模型中的參數(shù)進(jìn)行計算,最終確定具體的概率分布模型。
5)執(zhí)行以上步驟之后即可獲取航行中集裝箱船舶某2個搖擺位置處搖擺變化值的概率模型。
研究數(shù)據(jù)來自于上海海事大學(xué)萬噸級遠(yuǎn)洋教學(xué)實習(xí)船“育鋒”輪。該船的具體參數(shù)見表1。
為測定船舶的搖擺值,在艏部位置V1及舯部位置V2處安置角度傳感器(見圖1)。[9]
“育鋒”輪在正常航行過程中,其角度傳感器采集V1及V2處的縱向搖擺角度、橫向搖擺角度和垂向加速度,共持續(xù)3 h。在該時間段內(nèi),傳感器按照時間序列共觸發(fā)47次,每次觸發(fā)之后每個傳感器可連續(xù)采集并返回100個角度和加速度的樣本值。因此,每個測試點(diǎn)均有4 700個隨機(jī)采樣值,圖2和圖3分別為每次觸發(fā)時的角度及加速度大小的平均值序列。[8]
表1 “育鋒”輪具體參數(shù)
圖1 “育鋒”輪傳感器測點(diǎn)布置
2.2.1分布一致性判定
首先使用q-q圖判定V1和V2各方向上的變化值是否符合相同分布,結(jié)果見圖4。
由圖4可初步推斷V1和V2處縱向角度變化值、橫向角度變化值及垂直方向加速度變化值均服從相同分布。后續(xù)只對V1處進(jìn)行概率模型判定,V2省略。
2.2.2概率模型判定
利用直方圖技術(shù)對“育鋒”輪上2處傳感器采集的各方向上的搖擺角度及加速度變化值進(jìn)行初步統(tǒng)計分析和概率驗證,結(jié)果見圖5~圖7。
a)縱向搖擺角度均值
b)橫向搖擺角度均值
c)垂向加速度均值
圖2V1點(diǎn)在傳感器依此觸發(fā)時的角度及加速度變化值的平均序列
a)縱向搖擺角度均值
b)橫向搖擺角度均值
c)垂向加速度均值
a)縱向搖擺角度均值
b)橫向搖擺角度均值
c)垂向加速度均值
圖4V1及V2各方向?qū)?yīng)的q-q圖
a)縱搖角度變化頻率分布
b)縱搖角度變化正太分布概率驗證圖
a)橫搖角度變化頻率分布
b)橫搖角度變化正太分布概率驗證圖
a)垂向加速度變化頻率分布
b)V1垂向加速度正太分布概率驗證圖
由圖5可知,其縱搖變化值近似服從高斯分布。利用normplot函數(shù)進(jìn)行驗證,結(jié)果顯示‘+’號均在參考線附近,說明樣本觀測數(shù)據(jù)近似服從正太分布,與直方圖初步估計的結(jié)果一致,同時圖6和圖7中均有相同的結(jié)果。因此,該處的搖擺角度變化值及加速度變化值近似服從高斯分布是可信的。
2.2.3模型參數(shù)估計及分析
根據(jù)上述校驗的各方向上的搖擺角度概率模型結(jié)果,采用最大似然估計法確定相應(yīng)的高斯模型的參數(shù)(主要是總體均值和方差),具體見表2。
表2 搖擺角度及加速度概率模型參數(shù)
由表2可知:測量位置V1及V2處的縱向搖擺角度變化值趨向相同;橫向搖擺角度變化值相差較大,說明船舶在航行過程中受橫風(fēng)、橫浪影響較大(與“育鋒”輪航行過程監(jiān)控視頻所得結(jié)果一致);垂向的加速度變化值V1處較V2處波動大,與實際航行過程中艏部垂向變化幅度比舯部位置大相一致。
探索一種新型集裝箱船舶安全航行過程中不同位置處的搖擺變化概率模型,該模型包括縱向及橫向的角度變化概率模型和垂向的加速度變化概率模型。研究發(fā)現(xiàn),集裝箱船舶在安全航行過程中,其關(guān)鍵位置處的搖擺角度及加速度模型均服從不同參數(shù)的高斯分布。該研究從集裝箱船舶搖擺的角度及加速度2方面為船舶結(jié)構(gòu)安全和航行安全研究提供了一定的理論依據(jù)。
[1] 高儒. 集裝箱船綜合安全評估(FSA)研究[D]. 大連:大連海事大學(xué), 2013.
[2] 萬培峰. 砰擊和波浪誘導(dǎo)響應(yīng)的組合方法及考慮維修的船體可靠性分析[D]. 武漢:武漢理工大學(xué), 2001.
[3] 孫德建, 胡雄, 金永興,等. 集裝箱船舶結(jié)構(gòu)狀態(tài)的在線監(jiān)測技術(shù)研究[C]// “工程與振動”科技論壇, 2009.
[4] 金永興, 張世斌. 集裝箱船體波浪載荷應(yīng)力短期分布的統(tǒng)計分析[J]. 上海交通大學(xué)學(xué)報, 2012(8):1243-1247.
[5] 趙嶷飛, 王晨, 王紅勇. 基于ADS-B的航跡誤差分布規(guī)律研究[J]. 中國民航大學(xué)學(xué)報, 2012, 30(6):48-52.
[6] CHEN X R, ZHAO L C. Almost Sure L 1-Norm Convergence for Data-Based Histogram Density Estimates[J]. Journal of Multivariate Analysis, 1987, 21(1):179-188.
[7] 肖筱南. 新編概率論與數(shù)理統(tǒng)計[M]. 北京:北京大學(xué)出版社, 2013.
[8] 金永興, 胡雄, 施朝健. 集裝箱船舶結(jié)構(gòu)狀態(tài)監(jiān)測與評估系統(tǒng)[J]. 上海海事大學(xué)學(xué)報, 2008, 29(3):1-4.
[9] 金永興, 伍生春. 船舶結(jié)構(gòu)與設(shè)備[M]. 北京:人民交通出版社, 2012.
[10] 王倪傳,石玉虎,覃聞銘.船舶結(jié)構(gòu)不同應(yīng)力點(diǎn)的應(yīng)力概率模型[J]. 上海海事大學(xué)學(xué)報,2015,36(4):62-67.
ProbabilisticModelforSwayingatKeyPositionsonContainerShip
WANGPeiliang1,2,WUXiaofang3,ZHANGTing4
(1. Engineering Research Center of Simulation Technology of Ministry of Education, Shanghai Maritime University, Shanghai 201306, China; 2. Mechanical Engineering Department, Weifang University of Science and Technology, Weifang 262700, China; 3. China Waterborne Transport Research Institute, Beijing 100089, China; 4. Marine Academy, Shandong Transport Vocational College, Weifang 262700, China)
To reveal the inherent relationship between the operation safety of the ship and the swaying of the ship, especially the rolling, a probabilistic model for swaying at the key positions on container ships in normal navigation is constructed. The corresponding parameters of the probabilistic model are studied and defined. It is concluded that the probabilistic model regarding ship swaying follows Gaussian distribution, but the parameters can vary widely, dependent on ship design. The proposed method is of value for the study of structure safety and navigation safety of container ships at sea.
container ship; ship swaying; Gaussian distribution; navigation safety
2017-10-25
上海市科學(xué)技術(shù)委員會項目(14170501600);山東省教育科學(xué)研究課題(16SB057)
王培良(1987—),男,山東濰坊人,博士生,研究方向為載運(yùn)工具運(yùn)用工程。E-mail:gfy5216@126.com
1000-4653(2018)01-0070-04
U661.32+1
A