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

1護面塊體位姿變換測量原理

使用IMU慣性傳感器內部的加速度計和陀螺儀分別測量出護面塊體在受到波浪沖擊瞬間產生的加速度值和角速度值⑷。由于IMU慣性傳感器布置在護面塊體上，在進行試驗時所輸出的運動量為載體坐標系下的運動信息。載體坐標系⑷會隨波浪沖擊護面塊體的位置改變而發生變化，通過位姿變換可將塊體相對載體坐標系下的運動信息，經過姿態矩陣轉換到地理坐標系下，同時對地理坐標系下的信號進行處理，可使護面塊體的運動信息不會隨著本身位置的變化而發生改變。

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

通過姿態角的變化來描述空間中的轉換關系包含三次基本旋轉，每次旋轉的角度稱分別為航向角、俯仰角和橫滾角。

各個旋轉對應的變換矩陣表達式為

因此可以得到姿態矩陣表達式為：

化簡后得

2試驗布置及信號處理

試驗是在交通運輸部天津水運工程科學研究院大比尺波浪水槽的消能區進行的。水槽長456 m、寬5 m、高12 m。大比尺波浪水槽能模擬產生3. 5m高的波浪，造波周期為2〜10 s。該水槽能很好模擬出真實的波浪沖擊海岸護面塊體的過程，如圖2、圖3所示。

IMU慣性傳感器的集成模塊為高精度的陀螺儀、加速度計、地磁傳感器，采用高性能的微處理器和先進的動力學解算與卡爾曼動態濾波算法。同時它具有降低測量噪聲，提高測量精度、抗干擾能力強等優點。IMU的特性指標見表1。

試驗前，將SD存儲卡插入IMU慣性傳感器中，并將IMU慣性傳感器放入到已經固定在護面塊體上防水的盒子中。并將需要試驗的護面塊體放置在大比尺波浪水槽的消能區水線附近，如圖4、圖5所示。

圖3消能區

表1特性指標

圖4 IMU慣性傳感器和SD卡

試驗的不規則波具體工況見表2。

表2波高與水線位置

試驗結束后，將IMU慣性傳感器取出并將內置的SD卡插入讀卡器中與深圳維特公司開發的 JY901上位機軟件相關聯后可獲取加速度和角速度信號。該軟件能夠對在試驗過程中的信號進行實時動態顯示，也能顯示出加速度和角速度信號曲線圖，如圖6、圖7所25。

對于加速度信號縱坐標單位為g，角速度信號縱坐標單位為(°)/s,橫坐標為采集時間間隔，每一個正方形代表1s。

加速度信號在受到碰撞時會受到重力加速度影響導致結果不準確，利用現有的技術手段來處理產生的信號難以被修正，所以選取角速度信號作為處理信號，其優點是不受重力加速度影響。絕對角速度信號可由z、y、z3個軸產生的角速度分量合成購。

假設護面塊體繞某個點做純旋轉運動，塊體的公稱直徑為，因此可以得到沖擊速度。

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

3實驗結果分析

針對不規則波的工況，在不同水線位置可以得到穩定性數與碰撞數之間的對應關系。同時，分析了護面塊體在不同水線位置處的最大沖擊速度變化趨勢。

3.1穩定性數與碰撞數

波高Hs與護面塊體公稱直徑Dn的比值記穩定性系數H"/Dn,碰撞數NgjN由塊體受到的碰撞次數N®與波個數N的比值。其中碰撞次數 N罰可以定義為：在受波浪沖擊時所產生的沖擊速度中,統計出大于沖擊速度的平均值個數，如圖10所示。

圖10穩定性數與碰撞數對應關系

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

3.2不同水線位置產生的最大沖擊速度

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

圖11不同波高、不同水線處最大沖擊速度趨勢

從圖11中可以看出，在不同的水線位置處最大沖擊速度的趨勢均是先增加后減少，在z/Da = 0處有最大沖擊速度。同時，波高越大對護面塊體產生的最大沖擊速度越大。如果護面塊體長期在最大沖擊速度的沖擊作用下，會加快護面塊體的失穩以及破壞。

4結論

使用該慣性傳感器在交通運輸部天津水運工程科學研究院大比尺波浪水槽的消能區進行試驗，在接近真實環境的波浪運動條件下，研究了穩定性數與碰撞數之間的關系、不同水線位置處產生的最大沖擊速度。在對穩定性數與碰撞數研究中發現, 處于水線位置z/Da = 0時試驗塊體與周圍塊體發生的碰撞的劇烈程度要比在z/Dn =-2處小，最大沖擊速度在水線z/Dn = 0出現。通過研究可以為防波堤防護提供參考。

參考文獻

[1] 胡原野，王收軍，陳松貴.基于機器學習的斜坡堤越浪量預測方法研究[J].科技和產業,2021,21(2):218-224.

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

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

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

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

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

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

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

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

[10]王晨陽.基于慣性測量元件的波浪測量方法研究［D］. 天津:天津理工大學,2019.

[11]叮蔣雙.基于捷聯慣導的運動軌跡跟蹤技術研究［D］.哈爾濱：哈爾濱工業大學,2015.