沈文君，趙志娟，劉利琴，金瑞佳，陳漢寶

（1.交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所，港口水工建筑技術(shù)國(guó)家工程研究中心，天津 300456；2.中國(guó)直升機(jī)設(shè)計(jì)研究所，江西 景德鎮(zhèn)， 333011；3.天津大學(xué)水利工程仿真與安全國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津 300072）

0 引 言

船舶安全系靠泊并高效地進(jìn)行裝卸作業(yè)是港口設(shè)計(jì)和運(yùn)營(yíng)中的重要問題，泊位的泊穩(wěn)條件不僅通過限定允許作業(yè)波高來保證，同時(shí)也考慮船舶的六自由度運(yùn)動(dòng)量作為安全作業(yè)標(biāo)準(zhǔn)[1]。受到防波堤掩護(hù)的港口，泊位處的波浪方向多為順浪和斜浪，在系泊船舶六個(gè)自由度的運(yùn)動(dòng)分量中，橫蕩、縱蕩和橫搖是港內(nèi)系泊船舶運(yùn)動(dòng)的主要形式，容易超出規(guī)范[2]給出的限制。系泊船舶對(duì)波浪周期非常敏感，小噸位船舶對(duì)應(yīng)的固有周期也偏小，與波浪周期更易接近，系泊船舶產(chǎn)生大幅劇烈運(yùn)動(dòng)的可能性更大。而船舶運(yùn)動(dòng)量過大不僅影響裝卸作業(yè)的效率，還易導(dǎo)致纜繩斷纜等安全問題。因此，在港口設(shè)計(jì)和運(yùn)營(yíng)管理中，針對(duì)波浪周期對(duì)小型系泊船舶運(yùn)動(dòng)的影響開展系統(tǒng)研究是非常必要的，明確系泊船舶在不同波浪周期下的運(yùn)動(dòng)響應(yīng)特性，可為設(shè)計(jì)及運(yùn)營(yíng)提供技術(shù)支持。

國(guó)內(nèi)外已有不少學(xué)者針對(duì)波浪作用下系泊船舶的運(yùn)動(dòng)特性開展了研究，研究重點(diǎn)集中在波浪條件、船舶裝載度、系纜方式和碼頭結(jié)構(gòu)等因素對(duì)系泊船舶運(yùn)動(dòng)的影響，然而，針對(duì)小型船舶運(yùn)動(dòng)特性開展系列研究的還非常少。

Ligteringen 等[3]在分析長(zhǎng)周期波浪作用下系泊船舶物理模型和原型測(cè)量數(shù)據(jù)的基礎(chǔ)上，提出了長(zhǎng)周期波浪作用下系泊船舶縱蕩、橫蕩和橫搖三個(gè)自由度運(yùn)動(dòng)量的估算公式。Van der Molen[4]等采用六自由度模型對(duì)澳大利亞西北部Withnell 海灣系泊LNG船舶的運(yùn)動(dòng)響應(yīng)和系泊荷載進(jìn)行了數(shù)值模擬研究，研究結(jié)果給出了涌浪作用下系泊船舶回轉(zhuǎn)運(yùn)動(dòng)的估算公式；Van der Molen等[5-6]還采用亞長(zhǎng)周期重力波模型，聯(lián)合切片理論對(duì)日本Tomakomai 港內(nèi)長(zhǎng)周期波浪作用下的系泊船舶運(yùn)動(dòng)響應(yīng)進(jìn)行了研究，計(jì)算表明，數(shù)值計(jì)算結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù)吻合較好，縱蕩運(yùn)動(dòng)的固有周期與港池振蕩周期接近，會(huì)導(dǎo)致系泊船舶的大幅運(yùn)動(dòng)響應(yīng)；Prashant Kumar 等[7]建立了一種基于港口的幾何結(jié)構(gòu)以及相關(guān)波浪的繞射和折射等特性下分析系泊船舶運(yùn)動(dòng)的水動(dòng)力模型，研究了系泊船舶在不同共振條件下不同方向的小振幅長(zhǎng)波的運(yùn)動(dòng)；結(jié)果表明，入射波的方向和共振模態(tài)是分析港內(nèi)固定位置系泊船舶運(yùn)動(dòng)的關(guān)鍵參數(shù)，根據(jù)波場(chǎng)信息，可在港口或港口站找到安全位置，指導(dǎo)船舶更好地進(jìn)行裝卸作業(yè)。

楊憲章[8]針對(duì)受涌浪和長(zhǎng)周期波影響嚴(yán)重的毛里塔尼亞友誼港，對(duì)波浪和系泊船舶的運(yùn)動(dòng)進(jìn)行了研究，研究指出，系泊船舶對(duì)低頻長(zhǎng)周期波浪響應(yīng)敏感，波列中的長(zhǎng)周期波浪成分可引起船舶的劇烈運(yùn)動(dòng)，從而給碼頭系泊及裝卸作業(yè)造成危害；劉必勁[9]對(duì)20萬(wàn)噸級(jí)散貨船和聶鴻鵬[10]對(duì)40萬(wàn)噸級(jí)油輪的系泊試驗(yàn)研究均表明船舶的運(yùn)動(dòng)量隨著波浪周期的增大而增大；Shi 等[11]研究了橫浪作用下26.6×104m3LNG 船舶的系泊運(yùn)動(dòng)響應(yīng)特性，研究結(jié)果表明船舶的橫蕩運(yùn)動(dòng)為周期性運(yùn)動(dòng)，橫蕩自振周期與船舶橫搖固有周期比值在1.11～1.23 之間，橫蕩和橫搖運(yùn)動(dòng)峰值隨著波浪周期的增大而逐漸增大；李越[12]利用OPTIMOOR 軟件對(duì)20 萬(wàn)噸系泊散貨船的運(yùn)動(dòng)進(jìn)行了研究，計(jì)算結(jié)果表明，不同運(yùn)動(dòng)量對(duì)于波浪周期有著各自的響應(yīng)特征，橫蕩和升沉運(yùn)動(dòng)隨著波浪周期的增大先增大然后逐漸減??；Shen等[13]研究了波浪周期和波高對(duì)大型系泊船舶的影響，詳細(xì)分析了運(yùn)動(dòng)隨波浪周期的變化，并對(duì)危險(xiǎn)波浪周期進(jìn)行了預(yù)測(cè)，計(jì)算結(jié)果表明在橫浪作用下，橫搖和橫蕩運(yùn)動(dòng)占主導(dǎo)作用，在計(jì)算周期范圍內(nèi)橫搖運(yùn)動(dòng)隨著波浪周期的增大先增大后減小，橫蕩運(yùn)動(dòng)隨著波浪周期的增大而增大，當(dāng)波浪周期達(dá)到一定數(shù)值后，數(shù)值變化緩慢；史憲瑩等[14]以一艘26.6×104m3LNG 船舶為對(duì)象，對(duì)混合波浪作用下系泊船舶的運(yùn)動(dòng)響應(yīng)規(guī)律進(jìn)行了試驗(yàn)研究，結(jié)果表明當(dāng)混合浪能量一定時(shí)，系泊運(yùn)動(dòng)量大多隨著低頻部分波浪能量的增大而增大；沈雨生等[15]采用物理模型試驗(yàn)對(duì)一艘10 萬(wàn)噸級(jí)散貨船進(jìn)行了不同周期下的系列試驗(yàn)，研究結(jié)果表明，系泊船舶橫蕩運(yùn)動(dòng)量總體上均隨著波浪周期的增大而增大，橫搖和升沉運(yùn)動(dòng)量隨著波周期的增大先增大后減小。以上工作表明，以往工程建設(shè)多以大型船舶系泊的泊穩(wěn)情況為關(guān)注點(diǎn)，因此學(xué)者多選擇大噸位船舶與波浪周期的影響關(guān)系進(jìn)行研究，同時(shí)針對(duì)該耦合系統(tǒng)的運(yùn)動(dòng)頻率特性開展討論的文獻(xiàn)較少。近幾年某些小型碼頭的建設(shè)逐漸增多，如國(guó)內(nèi)LNG運(yùn)輸船、國(guó)外電廠碼頭小型駁船或小型集裝箱船等，這些船的固有周期偏小，對(duì)常見的波浪周期敏感性更高。

因此本文選取一艘小型駁船進(jìn)行研究，通過對(duì)不同波浪作用角度、不同波浪周期下船舶與纜繩及護(hù)舷組成的非線性耦合系統(tǒng)響應(yīng)特性進(jìn)行模擬，旨在掌握該船舶的運(yùn)動(dòng)響應(yīng)規(guī)律；然后通過頻譜分析得到系統(tǒng)的頻率特性，并基于小波分析，進(jìn)一步了解該耦合系統(tǒng)在時(shí)域內(nèi)的頻譜變化規(guī)律。

1 數(shù)值計(jì)算

1.1 數(shù)值理論及方程的建立

系泊船舶在環(huán)境載荷、系泊載荷和護(hù)舷荷載的共同作用下處于平衡，本文通過Hydrostar軟件獲得頻域內(nèi)波浪激振力、附加質(zhì)量及阻尼系數(shù)，然后基于時(shí)頻轉(zhuǎn)換理論，進(jìn)行傅立葉變換得到時(shí)域下波浪作用力、附加質(zhì)量和遲滯函數(shù)，進(jìn)而在時(shí)域內(nèi)建立系泊浮體的運(yùn)動(dòng)響應(yīng)方程，如式（1）所示。

式中：[aij]為浮體的慣性質(zhì)量矩陣；[mij(t)]為浮體的附加質(zhì)量矩陣；[Kij(t)]為延遲函數(shù)矩陣；[Cij]為靜水恢復(fù)力矩陣；[xj(t)]為浮體六自由度位移矩陣；Fwave(t)為波浪激勵(lì)力矩陣，包含一階力和二階慢漂力；Fcurrent(t)為流力矩陣，根據(jù)OCIMF 中的流力系數(shù)計(jì)算；Fwind(t)為風(fēng)激勵(lì)力矩陣，根據(jù)受風(fēng)面積和OCIMF 中的風(fēng)力系數(shù)計(jì)算；Ffender(t)為船舶所受的護(hù)舷擠靠力矩陣，根據(jù)圖1 中的護(hù)舷性能曲線計(jì)算；Fmoor(t)為纜繩拉力荷載，根據(jù)圖1中的纜繩性能曲線計(jì)算。

圖1 纜繩與護(hù)舷性能曲線Fig.1 Performance curve of mooring line and fender

延遲函數(shù)矩陣[Kij(t)]為

該函數(shù)可由頻域水動(dòng)力求解處的輻射阻尼Bij(ω)經(jīng)傅里葉逆變換求出。

纜繩和護(hù)舷荷載的計(jì)算方法類似，以纜繩荷載的計(jì)算為例進(jìn)行解釋說明，再計(jì)算模擬的每個(gè)時(shí)間步長(zhǎng)，計(jì)算每根系泊纜繩的系船柱至導(dǎo)纜孔的水平距離和方位角，方位角根據(jù)圖2確定。

圖2 導(dǎo)纜孔與系船柱之間的系泊纜繩Fig.2 Mooring line between fairlead and anchor

式中：xFk、yFk、zFk分別為第k根系泊纜繩導(dǎo)纜孔在船體坐標(biāo)系中的坐標(biāo)；Dh(k)為第k根系泊纜繩系船柱至導(dǎo)纜孔的水平距離；Ak為第k根系泊纜繩導(dǎo)纜孔至系船柱方向與固定坐標(biāo)系N軸的方位角；XAEk為第k根系泊纜繩系船柱在固定坐標(biāo)系中的坐標(biāo)；(XE,XN)為坐標(biāo)原點(diǎn)O在固定坐標(biāo)系中的坐標(biāo)；Ψ為船頭與固定坐標(biāo)系N軸的方位角。

通過式（3）計(jì)算出纜繩長(zhǎng)度，進(jìn)而得到纜繩長(zhǎng)度的變化量，然后根據(jù)纜繩的變形與拉力曲線（圖1）即可得到纜繩的拉力，將所有纜繩的拉力加在一起得到方程（1）中的Fmoor，根據(jù)同樣的方法計(jì)算Ffender。

1.2 算例模型

計(jì)算所用駁船主尺度如表1 所示，水動(dòng)力計(jì)算模型如圖3 所示。該船系泊在碼頭上，碼頭為高樁透空結(jié)構(gòu)，鋼管樁直徑為1.0 m，間距為10 m?？紤]到這種細(xì)長(zhǎng)體結(jié)構(gòu)對(duì)波浪影響較小，因此未在Hydrostar 中建立碼頭的水動(dòng)力模型。

圖3 駁船水動(dòng)力網(wǎng)格模型Fig.3 Hydrodynamic model of the vessel

表1 浮體主尺度參數(shù)Tab.1 Parameters of principal dimension of the vessel

纜繩為聚酰胺雙編織的纜繩，直徑為38 mm，最小破斷力為210 kN，安全工作荷載為105 kN。選用SANPE1.1 型護(hù)舷，設(shè)計(jì)反力為1081 kN，間距為6.6 m。纜繩和護(hù)舷的性能模擬曲線如圖1 所示。船舶?？看a頭系泊作業(yè)時(shí)的系泊布置如圖4所示。

圖4 作業(yè)時(shí)系泊布置方案圖（3：2）Fig.4 Mooring layout plan during operation（3：2）

1.3 計(jì)算工況

試驗(yàn)船舶裝載為滿載，分別進(jìn)行0°、30°和45°三種浪向下不同周期的試驗(yàn)組次，試驗(yàn)組次見表2。

表2 計(jì)算工況參數(shù)表Tab.2 Calculation cases

2 計(jì)算結(jié)果與分析

2.1 計(jì)算結(jié)果統(tǒng)計(jì)

Mamoun 等[16]開展了一艘液化天然氣船作業(yè)于15 m 水深時(shí)的水動(dòng)力荷載計(jì)算研究，對(duì)比了7個(gè)主流商業(yè)計(jì)算軟件的計(jì)算結(jié)果，結(jié)果顯示Hydrostar軟件展現(xiàn)了良好的計(jì)算性能，因此本文未再進(jìn)行數(shù)值驗(yàn)證。

將頻域的水動(dòng)力結(jié)果經(jīng)過傅里葉變換后，代入到式（1），計(jì)算過程采用全耦合計(jì)算方法，即船舶運(yùn)動(dòng)響應(yīng)與系泊纜的動(dòng)力響應(yīng)完全耦合分析，船舶波頻運(yùn)動(dòng)和低頻運(yùn)動(dòng)在時(shí)域范圍內(nèi)共同求解，進(jìn)而得到船舶的運(yùn)動(dòng)響應(yīng)曲線。圖5 為不同浪向時(shí)船舶橫搖和縱搖運(yùn)動(dòng)響應(yīng)RAO，圖6 為船舶在45°浪向、譜峰周期Tp為16 s時(shí)的橫蕩和橫搖運(yùn)動(dòng)時(shí)歷曲線。

圖5 橫搖和縱搖頻域計(jì)算結(jié)果Fig.5 Calculation results of roll and pitch in frequency domain

將表2中的所有工況進(jìn)行模擬計(jì)算，每組工況重復(fù)3次，統(tǒng)計(jì)分析3次結(jié)果最大值的平均值，觀察船舶運(yùn)動(dòng)響應(yīng)隨波浪周期、浪向等的變化規(guī)律。

2.2 不同波浪因素對(duì)運(yùn)動(dòng)響應(yīng)的影響

2.2.1 波浪周期對(duì)運(yùn)動(dòng)量的影響

按照上述統(tǒng)計(jì)方法，三個(gè)浪向下譜峰周期與六自由度運(yùn)動(dòng)量的關(guān)系如圖7 所示。從圖中可以看出，不同浪向下的規(guī)律基本相同，船舶的橫搖運(yùn)動(dòng)均隨著譜峰周期的增大先增大后減小再增大，三個(gè)浪向橫搖運(yùn)動(dòng)的最大值對(duì)應(yīng)的譜峰周期不同，0°、30°、45°浪向橫搖運(yùn)動(dòng)對(duì)應(yīng)的峰值周期分別在10 s、9 s、8 s附近。其他自由度基本隨著譜峰周期的增大而增大。

圖7 波浪周期對(duì)六自由度運(yùn)動(dòng)量的影響Fig.7 Influence of wave period on motion of six degrees of freedom

2.2.2 浪向?qū)\(yùn)動(dòng)量的影響

在三個(gè)浪向下將六個(gè)自由度運(yùn)動(dòng)規(guī)律分別進(jìn)行對(duì)比，結(jié)果如圖8所示。從圖中可以看出，橫蕩和艏搖運(yùn)動(dòng)隨著波浪入射角度的增大而增大，在0°、30°和45°浪向下，波浪周期為16 s 時(shí)橫蕩運(yùn)動(dòng)分別為0.7 m、1.21 m 和2.16 m，艏搖運(yùn)動(dòng)分別為1.19°、1.98°和2.87°；縱蕩運(yùn)動(dòng)在三個(gè)波浪入射角度下的結(jié)果比較接近，0°波浪對(duì)應(yīng)的縱蕩略大；在波浪周期為4～8 s 時(shí)，船舶的垂蕩運(yùn)動(dòng)非常小，三個(gè)波浪入射角度下基本一致，當(dāng)周期大于8 s 后，垂蕩運(yùn)動(dòng)基本隨著波浪入射角度的增大而增大，波浪周期為16 s時(shí)三個(gè)角度對(duì)應(yīng)的垂蕩運(yùn)動(dòng)分別為0.38 m、0.45 m 和0.66 m；在波浪周期為4～8 s時(shí)，船舶的橫搖運(yùn)動(dòng)隨著波浪入射角度的增大而增大，例如波浪周期為8 s時(shí)三個(gè)角度對(duì)應(yīng)的橫搖運(yùn)動(dòng)分別為0.48°、0.78°和0.97°；波浪周期為8～10 s時(shí)，30°對(duì)應(yīng)的橫搖角度最大，45°次之，0°對(duì)應(yīng)的橫搖角度最小，例如波浪周期為9 s時(shí)三個(gè)角度對(duì)應(yīng)的橫搖運(yùn)動(dòng)分別為0.67°、1.05°和0.93°，而波浪周期大于10 s后，橫搖運(yùn)動(dòng)又隨著波浪入射角度的增大而增大，波浪周期為16 s 時(shí)三個(gè)角度對(duì)應(yīng)的橫搖運(yùn)動(dòng)分別為0.66°、0.89°和1.36°。

圖8 波浪入射角度對(duì)六自由度運(yùn)動(dòng)量的影響Fig.8 Influence of wave incident angle on motion of six degrees of freedom

2.3 頻率特性分析

為了掌握船舶運(yùn)動(dòng)的頻率特性，進(jìn)一步進(jìn)行頻譜分析，對(duì)比了代表性周期4 s 和16 s 在0°、30°和45°波浪入射角度時(shí)對(duì)應(yīng)的結(jié)果。從圖9 可以看出，當(dāng)受入射角度45°作用，波浪周期為4 s（0.25 Hz）時(shí)，除縱蕩運(yùn)動(dòng)表現(xiàn)出明顯的低頻特性（縱蕩固有頻率）外，其他自由度的運(yùn)動(dòng)均以波頻運(yùn)動(dòng)為主；波浪周期為16 s（0.062 5 Hz）時(shí)，橫蕩、垂蕩、縱搖和艏搖的主特征頻率為0.062 5 Hz，也還是以波頻為主，但縱蕩和橫搖兩個(gè)自由度的運(yùn)動(dòng)除了含有波頻成分（0.062 5 Hz）外，還出現(xiàn)了明顯的倍頻，縱蕩運(yùn)動(dòng)呈現(xiàn)出低頻（縱蕩固有頻率）和1倍波頻的特點(diǎn)，橫搖運(yùn)動(dòng)呈現(xiàn)出2個(gè)特征頻率，一個(gè)為1倍波頻，另一個(gè)在0.12 Hz 附近。從圖5 橫搖運(yùn)動(dòng)的RAO 曲線可知，橫搖運(yùn)動(dòng)的第一個(gè)峰值對(duì)應(yīng)的周期在8.5 s左右，從而判斷較長(zhǎng)周期的波浪激勵(lì)力不僅使系統(tǒng)產(chǎn)生了波頻運(yùn)動(dòng)，同時(shí)，也在固有頻率附近發(fā)生了振動(dòng)。對(duì)0°和30°對(duì)應(yīng)的結(jié)果也進(jìn)行了頻譜分析，與45°時(shí)的規(guī)律類似。

圖9 船舶運(yùn)動(dòng)的頻譜特性分析Fig.9 Analysis of frequency spectrum characteristics of ship motion

從上述分析可知，在較長(zhǎng)周期時(shí)，橫搖和縱蕩自由度的運(yùn)動(dòng)不僅含有與波浪激勵(lì)頻率一致的成分，還含有自身固有特性的頻率成分。為了進(jìn)一步分析縱蕩和橫搖運(yùn)動(dòng)頻率特性在時(shí)域內(nèi)的變化，引入Morlet 小波變換，一般來說，任意級(jí)數(shù)ξ（t）的連續(xù)小波變換可以通過在a尺度上展開（與頻率成反比）并在τ中平移母小波ψ（t）得到

式中，星號(hào)（*）表示復(fù)共軛。計(jì)算中，選取Morlet小波作為母小波基函數(shù)，

按照上述方法，縱蕩和橫搖運(yùn)動(dòng)的小波分析結(jié)果見圖10～11。從圖中可以看出，系泊船舶的運(yùn)動(dòng)頻率在較高周期時(shí)具有明顯的多頻和時(shí)變特性。

圖10 Tp=4 s時(shí)，縱蕩和橫搖運(yùn)動(dòng)時(shí)頻分析（45°浪向）Fig.10 Time frequency analysis of surge and roll motion with a wave period of 4 s(45°wave direction)

圖10是4 s波浪周期對(duì)應(yīng)的結(jié)果，縱蕩運(yùn)動(dòng)的主要特征頻率在0.03 Hz左右（縱蕩固有頻率），橫搖運(yùn)動(dòng)的主要特征頻率范圍稍大，主要在0.15～0.35 Hz 之間（波浪頻率）。圖11 是16 s 波浪周期對(duì)應(yīng)的結(jié)果，縱蕩運(yùn)動(dòng)的主要特征頻率有兩個(gè)范圍，能量較多的部分還是集中在0.03 Hz左右，能量少部分集中在0.05～0.1 Hz之間（波浪頻率）；橫搖運(yùn)動(dòng)的主要特征頻率也有兩個(gè)范圍，兩個(gè)范圍對(duì)應(yīng)的能量相差不太大，一部分集中在0.05～0.1 Hz之間（波浪頻率），一部分集中在0.1～0.15 Hz之間（橫搖固有頻率），更加直觀地顯示了系泊船舶的時(shí)頻變化特性。

圖11 Tp=16 s時(shí)，縱蕩和橫搖運(yùn)動(dòng)時(shí)頻分析（45°浪向）Fig.11 Time frequency analysis of surge and roll motion with a wave period of 16 s（45°wave direction）

3 結(jié) 論

為了掌握小型船舶運(yùn)動(dòng)響應(yīng)特性與波浪周期之間的關(guān)系，本文基于頻域勢(shì)流理論獲得了該船的水動(dòng)力荷載，然后基于時(shí)頻轉(zhuǎn)換理論，建立了系泊船舶與護(hù)舷、系泊纜繩組成的非線性耦合系統(tǒng)的動(dòng)力學(xué)方程，針對(duì)不同波浪作用角度、波浪周期，對(duì)該系統(tǒng)進(jìn)行了模擬分析。主要結(jié)論如下：

（1）船舶的橫搖運(yùn)動(dòng)均隨著譜峰周期的增大先增大后減小再增大；隨著浪向角的增大，橫搖運(yùn)動(dòng)的最大值對(duì)應(yīng)的譜峰周期逐漸變??；不同波浪入射角度下的規(guī)律大體相同。

（2）波浪周期為4～8 s 時(shí)，船舶的橫搖運(yùn)動(dòng)隨著波浪入射角度的增大而增大；波浪周期為8～10 s時(shí)，30°對(duì)應(yīng)的橫搖角度最大，45°次之，0°對(duì)應(yīng)的橫搖角度最??；波浪周期大于10 s 后，橫搖運(yùn)動(dòng)又隨著波浪入射角度的增大而增大。

（3）短周期波浪作用時(shí)，除縱蕩外，船舶的其它五個(gè)自由度的運(yùn)動(dòng)均以波頻運(yùn)動(dòng)為主，主要是波浪強(qiáng)迫激勵(lì)運(yùn)動(dòng)；隨著波浪周期的增大，縱蕩和橫搖兩個(gè)自由度的運(yùn)動(dòng)含有雙頻率特征，除波浪強(qiáng)迫激勵(lì)運(yùn)動(dòng)外，還存在自身的固有運(yùn)動(dòng)頻率成分。