推搖組合式造波機參數(shù)研究與結構設計

張 楠,尹 闊,董志奎,姜文通,梅 林

(燕山大學 a.機械工程學院;b.電氣工程學院,河北 秦皇島 066004)

0 引 言

我國海域遼闊，海岸線長達1.5萬km，管轄海域約300萬km2，同時還有6 536個500 m2以上的島嶼，是一個名副其實的海洋大國。我國海洋事業(yè)依然面臨嚴峻挑戰(zhàn)，從海洋大國到海洋強國的征程任重而道遠。我國長期“重陸輕海”，這種傳統(tǒng)觀念使我們嚴重缺乏對海洋的重視，充分利用海洋資源已經(jīng)成為我們發(fā)展的迫切需要。

在海洋工程和港口與近海工程實驗研究領域，由于地理因素以及海洋環(huán)境的復雜性與不確定性，戶外實地實驗往往具有一定的風險[1]。研究人員通常采用實驗室水槽造波技術來模擬海洋波浪，波浪的模擬是一種重要的實驗手段，而造波機是進行水槽實驗的必備實驗裝置[2-3]。根據(jù)造波板運動形式的不同，可以分為平推式和搖擺式造波機。平推式造波機通過造波板沿水平方向的平推運動實現(xiàn)波浪的模擬；搖擺式造波機通過造波板繞底部鉸鏈搖擺運動實現(xiàn)波浪的模擬[4-6]。平推式和搖擺式造波機由于造波原理不同，波形也各有特點，它們分別應用于不同的實驗環(huán)境，由于實驗要求的復雜性和多樣性，在實驗需要不同形式的波浪時，單一運動形式的造波機已不能滿足要求，只能采取更換造波機的方法，提高了實驗成本。

本文提出了一種推搖組合式造波機的新方案，根據(jù)秦皇島近海海域7、8月份海浪浪高特點，對平推造波和搖擺造波兩種工作模式的主要設計及工作性能參數(shù)進行理論分析與推導，并利用MATLAB 2016 M語言建立推導結果的數(shù)學模型，計算得到平推造波和搖擺造波兩種工作模式的主要設計及工作性能參數(shù)。確定了推搖組合式造波機的設計方案，根據(jù)受力分析和仿真計算結果，驗證了此方案的可行性。

1 造波實驗水槽水深的確定

秦皇島面朝渤海，屬開放性海域，對北戴河老虎石外設立的波浪浮標站全年觀測數(shù)據(jù)進行統(tǒng)計分析[7]，發(fā)現(xiàn)小于0.5 m的波浪出現(xiàn)頻率為89.15%。其中7、8月份有效波高不超過0.5 m[8-9]，波浪的平均周期為1.6 s，最大周期為3.2 s[10]。

根據(jù)實際觀測數(shù)據(jù)，確定了推搖組合式造波機造波能力最大波高為0.5 m?；谖⒎ɡ碚摬⒏鶕?jù)破碎波高公式[11]和波浪的色散方程[12]討論水深和周期的關系，破碎波高和波的色散方程如下：

式中：Hmax為最大波高；hmin為最小水深；k為波數(shù)，2π/λ=k;T為波浪周期，2π/T=ω。通過MATLAB中的M語言進行編程，建立了最小水深和周期的關系，如圖1所示。

圖1 最小水深與周期關系

從圖1可以看到：隨著水深的增加，波浪周期無限接近于1.5 s，這說明在最大波高0.5 m的條件下，當水深大于2.29 m時，即使進一步增大水深，也不能形成周期小于1.5 s的波浪。從實驗室水槽建造的經(jīng)濟角度考慮，選擇水槽的水深為1.03 m，在此條件下，可以制造出周期1.6～3.2 s，最大波高為0.5 m的波浪。

2 推搖組合式造波機參數(shù)

設坐標原點位于靜水面與造波板初始位置的交點處，x軸與靜水面重合，y軸垂直向上，d為造波板浸深，坐標示意圖如圖2所示。

圖2 坐標示意圖

在確定實驗水槽水深為1.03 m后，分別討論平推模式和搖擺模式不同周期條件下的行程與造波板浸深的關系，波高與行程比值如下[13]：

(3)

式中：

L0={sinh(k0h)+(σ-1)sinh[k0(h-d)]+

σ[cosh(k0h-k0d)-cosh(k0h)]/k0d}

H為波高;h為水深;d為浸深;k0為波數(shù);S為造波機運動結構的行程，S=2Xa。當σ=0時，為平推造波模式;σ=1時，為搖擺造波模式。在波高和水深已知情況下，只有行程、浸深和周期未知，通過MATLAB M語言編程，建立了三者之間的關系。平推和搖擺造波模式下，行程和浸深的關系分別如圖3、4所示。

圖3 平推模式行程與浸深的關系圖4 搖擺模式行程與浸深的關系

分析計算結果可知，在最大波高和水深不變的情況下，隨著周期的減小，浸深的增大，造波機運動結構的行程相應減小。對于實驗設備來說，行程越小則絲杠和導軌的價格越低，可有效節(jié)約實驗設備加工和制造費用?；诖嗽瓌t，在最大波高0.5 m前提下，選擇周期1.6 s為對應于此波高的造波機設計周期。

由圖3、4對比可知：同一周期和浸深條件下，搖擺模式的行程大于平推模式的行程，搖擺模式行程接近于平推模式行程的2倍。

平推與搖擺模式下，不同浸深條件，造波機運動結構最大運動速度與周期的關系計算結果分別如圖5、6所示。在相同周期條件下，隨著浸深的增大，造波機運動結構的最大運動速度v相應減小，從絲杠的實際最大安全運動速度的角度出發(fā)，并結合增大浸深亦可減小行程的計算結果，最終確定造波板浸深為1 m。

圖5 平推模式速度與周期的關系圖6 搖擺模式速度與周期的關系

由圖5、6對比可知：同一浸深和周期條件下，搖擺模式的速度大于平推模式的速度，搖擺模式的速度接近于平推模式速度的兩倍。

在確定推搖組合式造波機最大造波波浪高度0.5 m(設計波高)、造波實驗水槽水深、造波機設計周期及造波板浸深，并相應確定推搖組合式造波機運動結構的最大行程和最大速度后，根據(jù)式(3)計算出本實驗設備推搖組合式造波機造波波高能力曲線圖，如圖7所示。由圖7可知，本實驗設備能造出周期1.6～3.2 s，對應波高0.5～0.22 m(平推模式)，0.5～0.19 m(搖擺模式)的波浪。

圖7 造波能力曲線

3 推搖組合式造波機結構設計與作用力分析

3.1 推搖組合式造波機結構設計

根據(jù)造波板浸深和最大波高，確定造波板凈高度為1.28 m，并進一步根據(jù)推搖組合式造波機運動結構的行程及工作特點，設計推搖組合式造波機機械結構如圖8所示。

造波機結構材料采用45號鋼，造波板質(zhì)量207 kg，造波架質(zhì)量122 kg。即平推模式下，造波機運動結構質(zhì)量為329 kg；搖擺模式下，造波機運動結構質(zhì)量為207 kg。

3.2 推搖組合式造波機作用力分析

根據(jù)微幅波理論和伯努利方程[12]，推導得到造波板的壓力場公式：

圖8 推搖組合式造波機機械結構示意圖

(4)

速度勢φ(x,y,t)函數(shù)為：

φ(x,y,t)=cosωt·C0cosh[k0(y+h)]sink0x+

(5)

(6)

(7)

式中：

Ln=sin(knh)+(σ-1)sin[kn(h-d)]+

σ[cos(knh)-cos(knh-knd)]/(knd)

造波板在運動過程中會受到規(guī)則波力、擾動駐波力以及造波板慣性力的作用，已知推出的速度勢φ(x,y,t)函數(shù)中包括規(guī)則波和擾動駐波兩部分，分別對規(guī)則波和擾動駐波沿造波板面板積分得到平推模式的規(guī)則波力和擾動駐波力。

平推模式規(guī)則波力：

sinh[k0(h-d)]}/k0

(8)

(9)

平推模式擾動駐波力：

sin[kn(h-d)]}/kn

(10)

(11)

式中：B為造波板寬度，本設計方案造波板寬度為0.6 m;kn通過ω2=-gkntan(knh)求得[14]。

搖擺模式造波過程中[15]，不同水深對造波板底部鉸鏈存在一個力臂(y+d)，電動缸對造波板底部鉸鏈也存在一個力臂(d1+d)，d1為電動缸與坐標原點的垂直距離，本案例d1=0.3 m，因而推導出的搖擺模式規(guī)則波力和擾動駐波力如下：

搖擺模式規(guī)則波力：

(12)

C3=[k0dsinh(k0h)+cosh(k0h-k0d)-

cosh(k0h)]2

(13)

搖擺模式擾動駐波力：

(14)

C4=[kndsin(knh)-cos(knh-knd)+

cos(knh)]2

(15)

造波板慣性力通過下式求得:

(16)

式中：y1為造波板質(zhì)心相對于坐標系原點的位置，本案例y1=-0.293 m；m為造波機運動結構質(zhì)量。

電動缸作用于造波板的總力：

F=F0cosωt+Fnsinωt+Fasinωt

(17)

最大單元負載力：

(18)

電動缸功率：

(19)

將推出的作用力公式通過MATLAB M語言進行編程，計算結果如圖9～11所示。

圖9 平推和搖擺規(guī)則波力圖10 平推和搖擺擾動駐波力

圖11 平推和搖擺慣性力

通過分析計算結果可知：同一波高、周期條件下，平推模式的規(guī)則波力大于搖擺模式的規(guī)則波力，平推模式的規(guī)則波力接近于搖擺模式規(guī)則波力的2倍；平推模式的擾動駐波力隨著周期增大，逐漸趨于零，搖擺模式的擾動駐波力在給定周期范圍內(nèi)則相對變化不大，比較穩(wěn)定；平推模式的慣性力大于搖擺模式的慣性力。

通過以上推導與計算，最終確定推搖組合式造波機的主要參數(shù)(設計周期條件下)如表1所示。

表1 推搖組合式造波機主要參數(shù)

3.3 推搖組合式造波機受力分析

由式(4)可知，造波板所受到的力是隨時間和空間變化的，因此采用ANSYS中的瞬態(tài)動力學分析模塊對其受力作仿真分析，將推導出的壓力場公式以Function函數(shù)的形式導入到ANSYS中，對應力和總變形進行求解,得到圖12～15結果。

圖12 平推模式總應力云圖圖13 平推模式總變形云圖

圖14 搖擺模式總應力云圖圖15 搖擺模式總變形云圖

工作過程中，造波機機械結構的最大變形為2.3 mm，最大應力為62.432 MPa，最大變形在本造波機機械運動最大允許誤差范圍內(nèi)，最大應力遠小于材料屈服極限335 MPa，從瞬態(tài)動力學分析的角度，認為本設計方案是可行的。

4 結 語

根據(jù)實驗造波需要，本文針對秦皇島近海海域7、8月份海浪特點，依據(jù)微幅波理論，研究了在波高一定的情況下造波實驗水槽水深與周期的關系，平推造波和搖擺造波兩種工作模式的主要設計及工作性能參數(shù)。利用MATLAB 2016 M語言建立平推造波和搖擺造波主要設計及工作性能參數(shù)的數(shù)學模型，計算出兩種工作模式下推搖組合式造波機運動結構的周期、浸深與行程，浸深、周期與速度，作用力與周期等的關系，得到了平推造波和搖擺造波兩種工作模式下的主要設計及工作性能參數(shù)。在此基礎上設計推搖組合式造波機機械結構，并對其做受力及瞬態(tài)動力學分析，仿真結果驗證了本設計方案的可行性，為加工制造推搖組合式造波機和開展實驗造波工作提供了重要理論基礎及設計依據(jù)。