為了掌握高海況下艦船的縱搖運動規(guī)律并對其進行準確預報,分析了艦船的縱搖運動規(guī)律。運用譜分析原理和頻率響應法,分別在4級、5級和6級海況下,根據(jù)海浪譜密度函數(shù)和艦船的縱搖頻率響應函數(shù),求得其縱搖運動譜密度函數(shù);再運用公式推導,得到其縱搖運動的時域函數(shù);使用MATLAB仿真軟件對艦船進行縱搖運動仿真,得到3種高海況條件下的艦船縱搖運動時域函數(shù)圖像,并驗證仿真結果的準確性。仿真結論:在高海況下,海況等級對艦船縱搖影響顯著,海況等級越高,艦船縱搖的角度值范圍就越大。以某型艦船為例的預報結果:在4級海況下,艦船縱搖的幅度在±6°之間;在5級海況下,艦船縱搖的幅度在±15°之間;在6級海況下,艦船縱搖幅度可以達到±30°之間。 

【文章來源】：數(shù)字海洋與水下攻防. 2020,3(04)

【文章頁數(shù)】：9 頁

【部分圖文】：

輸入與輸出的轉換Fig.1Conversionofinputandoutput

第4期張鈺，等：基于譜分析原理和頻率響應法的艦船縱搖運動分析·317·速度；U為離海面19.5m處的風速。本文采用式（13）作為海浪譜密度公式。4縱搖運動分析4.1運動方程的建立與求解運用切片法[8]對船舶在波中的縱搖和垂蕩的耦合運動進行分析。在分析船體某切片的受力前，首先對其切片的運動進行分析，當船體作為小幅度縱搖和垂蕩運動時，可以認為船體切片作垂向運動。對運動坐標系G-XYZ中X處船體橫剖面的運動進行分析，見圖2所示。切片隨重心上升了z，由于縱搖下降了Xsin，近似為X。此時剖面出的波面坐標為ζ，根據(jù)式（15）求得該剖面與波面的垂向相對位置為RzzX（15）圖2船體橫剖面切片的運動Fig.2Motionofthehullcrosssection當波浪運動導致切片受到外力并引起切片運動時，在對其進行分析時需要考慮流場內史密斯效應的影響因素[9]。本文使用等效波對其進行計算，而不是直接用表面波的升高來進行相關計算。實際應用中，通常將等效波的深度用其切片平均的吃水深度來表示[10]。因此，若水線以下切片面積為S，切片寬度為2b，則平均吃水為2mSTb（16）可知其等效波面方程為*ecos(cossin)mkTAekxkyt（17）對于切片來說，有3種流體動力作用在垂向運動的切片上，其中包括流體靜力、興波阻力以及附加慣性力。首先，考慮流體靜力，由于切片吃水的變化，單位長度的上浮力為*1F2gb(zX)（18）其次，考慮興波阻力，對式（15）進行微分可以得到切片垂向速度**dd()ddRXzzXzXtt（19）當船在以速度V在水中前進時，在半固定坐標系上觀察船的運動就相當于觀察該

·320·數(shù)字海洋與水下攻防第3卷4.4基于MATLAB的縱搖運動仿真根據(jù)上述海況條件，基于MATLAB軟件對某型艦船進行縱搖運動仿真。取采樣頻率f=50Hz，分別對4級，5級和6級海況下UUV的縱搖幅度進行隨機仿真，得到如下結果，如圖3-5所示。其橫坐標為時間t/s，縱坐標為縱搖幅度的角度值Ψ（t）/（°）。圖34級海況下艦船縱搖幅度Fig.3UUVpitchamplitudeatseastate4圖45級海況下艦船縱搖幅度Fig.4UUVpitchamplitudeatseastate5圖56級海況下艦船縱搖幅度Fig.5UUVpitchamplitudeatseastate6為了驗證縱搖運動仿真的準確性，本文利用由圖3-5中縱搖幅度時域函數(shù)求得的縱搖譜密度函數(shù)與根據(jù)相關參數(shù)計算得到的縱搖譜密度函數(shù)1()pS、2()pS、3()pS進行對比。對比

【參考文獻】：
期刊論文
[1]隨機海浪的數(shù)值仿真與頻譜分析[J]. 許景波,邊信黔,付明玉.  計算機工程與應用. 2010(36)
[2]基于逆變換的路面不平度仿真研究[J]. 劉獻棟,鄧志黨,高峰.  中國公路學報. 2005(01)
[3]船舶穩(wěn)性理論研究的方法及進展[J]. 張文斌,姚震球,蔣志勇.  華東船舶工業(yè)學院學報(自然科學版). 2002(01)
[4]高速排水型船的運動性能預報[J]. 段文洋,賀五洲.  清華大學學報(自然科學版). 2001(12)
[5]不規(guī)則波作用下海床面上的波浪壓力計算[J]. 別社安,趙子丹,王光綸.  水利學報. 1998(03)
[6]由JONSWAP譜和PM譜計算的風浪波高之間的關系[J]. 徐德倫.  海洋湖沼通報. 1987(01)

博士論文
[1]基于水動力在線預報的艦船減縱—橫搖方法研究[D]. 韓陽.哈爾濱工程大學 2017



本文編號：3088411