通常在海港處的防波堤上都會(huì)鋪一層護(hù)面塊體 ,這層護(hù)面塊體可以保護(hù)防波堤不會(huì)受到波浪沖擊而損壞，同時(shí)，也能起到消除波浪的作用。護(hù)面塊體受波浪沖擊會(huì)發(fā)生瞬時(shí)搖擺，長(zhǎng)期沖擊會(huì)發(fā)生位移變化頁(yè)。在實(shí)際環(huán)境下研究護(hù)面塊體的搖擺需要花費(fèi)很多人力、物力、財(cái)力，成本需求較大。所以，在實(shí)驗(yàn)水槽中進(jìn)行物理模型試驗(yàn)是科學(xué)研究的常用選擇。目前，護(hù)面塊體受波浪沖擊瞬間發(fā)生的角度變化與周?chē)�物體發(fā)生碰撞只有定性的分析并沒(méi)有定量的評(píng)價(jià)。王鐵凝等研究了不同波高下塊體的穩(wěn)定性 并且定性評(píng)價(jià)了塊體的位姿變化旳。李賀青等給出了塊體受波浪影響瞬間的穩(wěn)定性系數(shù)⑷�？梢钥闯觯� 對(duì)于塊體位姿的測(cè)量研究不足同時(shí)沒(méi)有具體的數(shù)值 來(lái)刻畫(huà)出塊體受波浪作用的變化規(guī)律。本文使用IMU慣性傳感器直接對(duì)護(hù)面塊體進(jìn)行研究能反映出塊體在受到?jīng)_擊時(shí)瞬間產(chǎn)生的加速度值與角速度值。在護(hù)面塊體外部安裝IMU慣性傳感器,當(dāng)波浪沖擊瞬間會(huì)隨護(hù)面塊體發(fā)生搖擺，并輸出角速度信號(hào)。通過(guò)對(duì)角速度信號(hào)處理得到?jīng)_擊速度和碰撞數(shù)。 本文研究了在不同水線(xiàn)處穩(wěn)定性數(shù)與碰撞數(shù)的關(guān)系、 不同水線(xiàn)處產(chǎn)生最大沖擊速度的變化趨勢(shì)，與前人研究相比本文給出塊體受波浪作用的變化規(guī)律并借助IMU慣性傳感器與護(hù)面塊體相結(jié)合的方式在試驗(yàn)過(guò)程中記錄塊體的位姿變化，試驗(yàn)結(jié)束后可以通過(guò)傳感 器內(nèi)數(shù)據(jù)查看塊體受波浪作用時(shí)的具體數(shù)值。

1護(hù)面塊體位姿變換測(cè)量原理

使用IMU慣性傳感器內(nèi)部的加速度計(jì)和陀螺儀分別測(cè)量出護(hù)面塊體在受到波浪沖擊瞬間產(chǎn)生的加速度值和角速度值⑷。由于IMU慣性傳感器布置在護(hù)面塊體上，在進(jìn)行試驗(yàn)時(shí)所輸出的運(yùn)動(dòng)量為載體坐標(biāo)系下的運(yùn)動(dòng)信息。載體坐標(biāo)系⑷會(huì)隨波浪沖擊護(hù)面塊體的位置改變而發(fā)生變化，通過(guò)位姿變換可將塊體相對(duì)載體坐標(biāo)系下的運(yùn)動(dòng)信息，經(jīng)過(guò)姿態(tài)矩陣轉(zhuǎn)換到地理坐標(biāo)系下，同時(shí)對(duì)地理坐標(biāo)系下的信號(hào)進(jìn)行處理，可使護(hù)面塊體的運(yùn)動(dòng)信息不會(huì)隨著本身位置的變化而發(fā)生改變。

載體坐標(biāo)系與地理坐標(biāo)系之間的關(guān)系可以拆分成三次簡(jiǎn)單的坐標(biāo)變換，首先固定Z軸并繞軸旋 轉(zhuǎn)角度0，然后固定X軸并繞軸旋轉(zhuǎn)0，最后繞Y軸旋轉(zhuǎn)y ,如圖1所示。

通過(guò)姿態(tài)角的變化來(lái)描述空間中的轉(zhuǎn)換關(guān)系包含三次基本旋轉(zhuǎn)，每次旋轉(zhuǎn)的角度稱(chēng)分別為航向角、俯仰角和橫滾角。

各個(gè)旋轉(zhuǎn)對(duì)應(yīng)的變換矩陣表達(dá)式為

因此可以得到姿態(tài)矩陣表達(dá)式為：

化簡(jiǎn)后得

2試驗(yàn)布置及信號(hào)處理

試驗(yàn)是在交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究 院大比尺波浪水槽的消能區(qū)進(jìn)行的。水槽長(zhǎng)456 m、寬5 m、高12 m。大比尺波浪水槽能模擬產(chǎn)生3. 5m高的波浪，造波周期為2〜10 s。該水槽能 很好模擬出真實(shí)的波浪沖擊海岸護(hù)面塊體的過(guò)程， 如圖2、圖3所示。

IMU慣性傳感器的集成模塊為高精度的陀螺 儀、加速度計(jì)、地磁傳感器，采用高性能的微處理器和先進(jìn)的動(dòng)力學(xué)解算與卡爾曼動(dòng)態(tài)濾波算法。 同時(shí)它具有降低測(cè)量噪聲，提高測(cè)量精度、抗干擾能力強(qiáng)等優(yōu)點(diǎn)。IMU的特性指標(biāo)見(jiàn)表1。

試驗(yàn)前，將SD存儲(chǔ)卡插入IMU慣性傳感器中，并將IMU慣性傳感器放入到已經(jīng)固定在護(hù)面塊體上防水的盒子中。并將需要試驗(yàn)的護(hù)面塊體放置在大比尺波浪水槽的消能區(qū)水線(xiàn)附近，如圖4、圖5所示。

圖3消能區(qū)

表1特性指標(biāo)

圖4 IMU慣性傳感器和SD卡

試驗(yàn)的不規(guī)則波具體工況見(jiàn)表2。

表2波高與水線(xiàn)位置

試驗(yàn)結(jié)束后，將IMU慣性傳感器取出并將內(nèi) 置的SD卡插入讀卡器中與深圳維特公司開(kāi)發(fā)的 JY901上位機(jī)軟件相關(guān)聯(lián)后可獲取加速度和角速度 信號(hào)。該軟件能夠?qū)υ谠囼?yàn)過(guò)程中的信號(hào)進(jìn)行實(shí) 時(shí)動(dòng)態(tài)顯示，也能顯示出加速度和角速度信號(hào)曲線(xiàn) 圖，如圖6、圖7所25。

對(duì)于加速度信號(hào)縱坐標(biāo)單位為g，角速度信號(hào)縱坐標(biāo)單位為(°)/s,橫坐標(biāo)為采集時(shí)間間隔，每一 個(gè)正方形代表1s。

加速度信號(hào)在受到碰撞時(shí)會(huì)受到重力加速度 影響導(dǎo)致結(jié)果不準(zhǔn)確，利用現(xiàn)有的技術(shù)手段來(lái)處理 產(chǎn)生的信號(hào)難以被修正，所以選取角速度信號(hào)作為處理信號(hào)，其優(yōu)點(diǎn)是不受重力加速度影響。絕對(duì)角 速度信號(hào)可由z、y、z3個(gè)軸產(chǎn)生的角速度分量 合成購(gòu)。

假設(shè)護(hù)面塊體繞某個(gè)點(diǎn)做純旋轉(zhuǎn)運(yùn)動(dòng)，塊體的公稱(chēng)直徑為，因此可以得到?jīng)_擊速度。

通過(guò)沖擊速度可以繪制出護(hù)面塊體受到波浪沖擊瞬間的沖擊峰，沖擊峰能直觀地反映出塊體受 到波浪沖擊時(shí)的速度分布以及受波浪沖擊時(shí)的最大、最小沖擊速度，如圖8、圖9所示。

3實(shí)驗(yàn)結(jié)果分析

針對(duì)不規(guī)則波的工況，在不同水線(xiàn)位置可以得 到穩(wěn)定性數(shù)與碰撞數(shù)之間的對(duì)應(yīng)關(guān)系。同時(shí)，分析了護(hù)面塊體在不同水線(xiàn)位置處的最大沖擊速度變化趨勢(shì)。

3.1穩(wěn)定性數(shù)與碰撞數(shù)

波高Hs與護(hù)面塊體公稱(chēng)直徑Dn的比值記穩(wěn)定性系數(shù)H"/Dn,碰撞數(shù)NgjN由塊體受到的 碰撞次數(shù)N®與波個(gè)數(shù)N的比值。其中碰撞次數(shù) N罰可以定義為：在受波浪沖擊時(shí)所產(chǎn)生的沖擊速 度中,統(tǒng)計(jì)出大于沖擊速度的平均值個(gè)數(shù)，如圖10所示。

圖10穩(wěn)定性數(shù)與碰撞數(shù)對(duì)應(yīng)關(guān)系

在保持塊體的公稱(chēng)直徑D"不變的條件下，穩(wěn)定性數(shù)會(huì)隨著波高的增加而增加，同時(shí)碰撞數(shù)也伴隨著增大。由圖10可知，碰撞數(shù)增加的幅度不大, 反映在圖中曲線(xiàn)的變化趨勢(shì)很小。對(duì)比表2中9組工況不同水線(xiàn)與波高位置處，穩(wěn)定性系數(shù)與碰撞數(shù)之間的關(guān)系會(huì)隨著穩(wěn)定性系數(shù)的增加，碰撞數(shù)也增加。而在水線(xiàn)z/Dn =- 2處的碰撞數(shù)明顯大于水 線(xiàn)z/D” = 0和z/Dn = 2，說(shuō)明水線(xiàn)z/Dn =2處塊體發(fā)生的搖擺比z/Dn = 0和z/Dn = 2處要?jiǎng)?烈，同時(shí)與周?chē)�的塊體互相碰撞的頻率也更高。在 穩(wěn)定性系數(shù)不變的條件下，z/Dn =-2處的碰撞數(shù) 大于z/Dn = 0和z/Dn = 2處的碰撞數(shù)。說(shuō)明水線(xiàn) 位置的選取會(huì)影響塊體在受到波浪沖擊時(shí)搖擺的劇烈程度。

3.2不同水線(xiàn)位置產(chǎn)生的最大沖擊速度

護(hù)面塊體安放在不同水線(xiàn)位置時(shí)對(duì)防波堤的 防護(hù)起著關(guān)鍵作用。同時(shí)，在不同水線(xiàn)位置處護(hù)面塊體受到的最大沖擊速度均不同。護(hù)面塊體安放過(guò)程中選擇合適的水線(xiàn)位置有助于提高護(hù)面塊體的壽命。為研究塊體在不同水線(xiàn)位置的最大沖擊速度進(jìn)行9組試驗(yàn)工況，在不同波高、不同水線(xiàn)位置 處受到的最大沖擊速度的變化趨勢(shì)，如圖11所示。

圖11不同波高、不同水線(xiàn)處最大沖擊速度趨勢(shì)

從圖11中可以看出，在不同的水線(xiàn)位置處最大沖擊速度的趨勢(shì)均是先增加后減少，在z/Da = 0處 有最大沖擊速度。同時(shí)，波高越大對(duì)護(hù)面塊體產(chǎn)生 的最大沖擊速度越大。如果護(hù)面塊體長(zhǎng)期在最大 沖擊速度的沖擊作用下，會(huì)加快護(hù)面塊體的失穩(wěn)以 及破壞。

4結(jié)論

使用該慣性傳感器在交通運(yùn)輸部天津水運(yùn)工 程科學(xué)研究院大比尺波浪水槽的消能區(qū)進(jìn)行試驗(yàn)， 在接近真實(shí)環(huán)境的波浪運(yùn)動(dòng)條件下，研究了穩(wěn)定性 數(shù)與碰撞數(shù)之間的關(guān)系、不同水線(xiàn)位置處產(chǎn)生的最大沖擊速度。在對(duì)穩(wěn)定性數(shù)與碰撞數(shù)研究中發(fā)現(xiàn), 處于水線(xiàn)位置z/Da = 0時(shí)試驗(yàn)塊體與周?chē)鷫K體發(fā) 生的碰撞的劇烈程度要比在z/Dn =-2處小，最大 沖擊速度在水線(xiàn)z/Dn = 0出現(xiàn)。通過(guò)研究可以為 防波堤防護(hù)提供參考。

參考文獻(xiàn)

[1] 胡原野，王收軍，陳松貴.基于機(jī)器學(xué)習(xí)的斜坡堤越浪量 預(yù)測(cè)方法研究[J].科技和產(chǎn)業(yè),2021,21(2):218-224. 

[2] 龔照龍，張曉宇.防波堤護(hù)面塊裂縫產(chǎn)生的原因與控制分 析[J].工程技術(shù)研究,2020,5(11):187-188. 

[3] 郭棟，柳淑學(xué),李金宣，等.波陡對(duì)于扭王字塊穩(wěn)定性影響 的試驗(yàn)研究[J].水道港口 ,2021,42(4):421-426.

[4]戈龍仔，高峰，陳漢寶.深水防波堤新型消浪護(hù)面塊體開(kāi) 發(fā)及其性能參數(shù)測(cè)試研究［J］.海岸工程,2020,39(3)： 204-212. 

[5]王鐵凝，蘇莉源，陳國(guó)平，等.波浪周期對(duì)扭王字塊體穩(wěn)定 性的影響［J1水運(yùn)工程,2019(4):45-50. 

[6]李賀青，柳玉良.從護(hù)面塊體穩(wěn)定性論深水防波堤［J］.水 運(yùn)工程,2014(6):1-5. 

[7]張洪，馬振浩，丁詩(shī)泳.基于IMU的機(jī)器人滑移軌跡優(yōu)化 研究口］.傳感器與微系統(tǒng),2021,40(4):37-40,44. M 

[8]王晨陽(yáng)，王收軍，陳松貴，等.基于加速度計(jì)和陀螺儀的波 浪測(cè)量方法口］・科學(xué)技術(shù)與工程,2019,19(6):44-49. 

[9]廖杰華.激光雷達(dá)/IMU組合導(dǎo)航定位方法研究［D］.南 昌：南昌大學(xué)，2020. 爾濱：哈爾濱工業(yè)大學(xué),2015.

[10]王晨陽(yáng).基于慣性測(cè)量元件的波浪測(cè)量方法研究［D］. 天津:天津理工大學(xué),2019. 

[11]叮蔣雙.基于捷聯(lián)慣導(dǎo)的運(yùn)動(dòng)軌跡跟蹤技術(shù)研究［D］.哈爾濱：哈爾濱工業(yè)大學(xué),2015.

[12]劉路，王收軍,陳松貴，等.基于MEMS加速度計(jì)的波高 測(cè)量裝置口]・科學(xué)技術(shù)與工程,2019,19(30):70-77. 

[13]陳鴻躍，陳雨，趙曉偉，等.一種單軸旋轉(zhuǎn)的車(chē)載IMU/DTU 組合導(dǎo)航方法口]・現(xiàn)代防御技術(shù),2021,49(1):6^74. 

[14]崔海路•基于LiDAR和IMU融合的智能車(chē)組合定位導(dǎo) 航技術(shù)研究[D].濟(jì)南:齊魯工業(yè)大學(xué),2021. 

[15]HOFLAND B,AREFIN S S,LEM C V D. Smart rocking armour units [C]// Selected papers Coastal Research, Netherlands: Coastlab ,2018： 5-6.

[16]VAN D, HEYDRA G. Impact velocities of rocking armour units[CJ// Concrete Armor Units. Netherlands： Coastal Engineering,1990：29-31. 

[17]LE T N. Rocking Revisited 1: Rocking of a singe cube on a breakwater slope[J]. Coastal Engineering, 2016,56 (4):389-399.