水下仿生柔性胸鰭擺動水動力特性試驗研究

張紀(jì)華，劉樂華，楊曉光，馬雪泉

（1.北京機電工程研究所，北京100074；2.上海船舶運輸科學(xué)研究所，上海200135）

0 引 言

水下仿生潛航器具有高效推進、高機動等良好性能，仿生胸鰭推進是模仿魚類左右對稱的柔性胸鰭運動而興起的一種仿生推進模式，典型的海洋生物是俗稱“魔鬼魚”的蝠鲼。針對水下仿生胸鰭推進集成，日本科研工作者Imae采用柔性乙烯樹脂薄膜和淬火鋼帶，以剛性雙四桿機構(gòu)驅(qū)動胸鰭前緣帶動胸鰭鰭面實現(xiàn)擺動運動[1]；德國Festo 公司研制了胸鰭擺動的仿生魚Aqua_ray[2]；新加坡國立大學(xué)選用PVC 片材料作為胸鰭研制了MantaDroid 仿蝠鲼機器魚[3]，并完成了長航時的續(xù)航試驗；新加坡南洋理工大學(xué)與中國國防科學(xué)技術(shù)大學(xué)聯(lián)合開發(fā)了仿蝠鲼機器魚，樣機兩側(cè)胸鰭分別由三組獨立控制的具有一定柔性的鰭條驅(qū)動，并采用類魚漂的機構(gòu)進行深度控制[4]；國防科技大學(xué)采用多鰭條驅(qū)動方式研制了仿胸鰭樣機[5]；哈爾濱工業(yè)大學(xué)研制了SMA 驅(qū)動的仿蝠鲼機器魚[6-7]，實現(xiàn)了高穩(wěn)定性和無噪聲推進；北京航空航天大學(xué)研制了多型仿牛鼻鲼Robo_Ray樣機[8]，并自行搭建了小型測力實驗水槽[9]，完成了樣機水池游動試驗和樣機測力試驗?，F(xiàn)有的研究成果中，仿生胸鰭推進樣機側(cè)重于仿生機構(gòu)設(shè)計、無模型的仿生胸鰭推進游動控制及系統(tǒng)集成[10]等，從研制流程上基本忽略了仿生胸鰭推進水動力設(shè)計這一極為關(guān)鍵的環(huán)節(jié)，這也是造成仿生胸鰭推進效率低下的重要原因。此外，對于仿生胸鰭推進水動力特性的研究，還主要依賴于流場的數(shù)值計算，且數(shù)值仿真計算缺乏柔性材質(zhì)的胸鰭與水流相互作用，造成水動力數(shù)據(jù)失真。為此，本文通過研制仿生柔性胸鰭模型及驅(qū)動控制系統(tǒng)，搭建試驗測試系統(tǒng)，采用拖曳水池測力試驗的方法對仿生柔性胸鰭擺動過程受力進行測試，獲得仿生柔性胸鰭在不同運動參數(shù)時的受力數(shù)據(jù)。

1 試驗裝置與方法

1.1 試驗水池

試驗水池選擇上海船舶運輸科學(xué)研究所船模拖曳試驗水池，如圖1所示。水池主要參數(shù)為：總長192 m，有效試驗長度160 m，寬10 m，水深4 m，拖車運行速度0.01～10 m/s。試驗水池配備NI 公司集成化數(shù)據(jù)實時采集分析系統(tǒng)。

1.2 試驗?zāi)Ｐ?/h3>

1.2.1 物理模型

本文依據(jù)生物胸鰭原型輪廓，研制了仿生柔性胸鰭試驗?zāi)Ｐ?，如圖2 所示。該模型由鰭條、硅橡膠蒙皮經(jīng)過澆筑、固化和脫膜等過程制作而成，硅橡膠確保試驗?zāi)Ｐ彤a(chǎn)生柔性光滑的變形運動。

試驗?zāi)Ｐ驼归L為1.29 m，根部最大弦長為1.1 m，根部剖面最大厚度為0.24 m。仿生胸鰭試驗?zāi)Ｐ偷那?、中、后三根鰭條外接三個伺服電動缸，試驗?zāi)Ｐ驮陔妱痈淄评z杠同步往復(fù)運動的帶動下產(chǎn)生左右擺動，在電動缸推拉絲杠異步等時差往復(fù)運動的帶動下產(chǎn)生左右擺動和鰭面扭轉(zhuǎn)，調(diào)節(jié)電動缸絲杠線位移、頻率以及異步時間差實現(xiàn)試驗?zāi)Ｐ瓦\動參數(shù)和運動形態(tài)的變化。

1.2.2 模型運動

仿生胸鰭試驗?zāi)Ｐ瓦\動由三根鰭條帶動，本文設(shè)計了鰭條的正弦周期性運動規(guī)律如下：

圖1 船模拖曳試驗水池Fig.1 Ship model towing tank

圖2 仿生胸鰭試驗?zāi)Ｐ虵ig.2 Experimental model of bionic flexible pectoral

式中，θ11，θ12，θ13分別為前中后三根鰭條與鉛垂面的夾角；A 為角度幅值，單位為°；f為擺動頻率，單位為Hz；Δ?為相鄰鰭條間的相位差，單位為°。圖3為前中后鰭條運動定義，F(xiàn)x和Fy分別為模型受到的推力和側(cè)向力。

1.3 試驗?zāi)Ｐ桶惭b

試驗?zāi)Ｐ桶惭b如圖4 所示。仿生柔性胸鰭試驗?zāi)Ｐ偷那?、中、后鰭條與對應(yīng)的電動缸相連，形成仿生胸鰭運動系統(tǒng)。整個運動系統(tǒng)以豎直垂吊的方式通過測力天平傳感器、支撐板等連接，被固定在水池拖車的模型安裝測試機構(gòu)上。鰭面轉(zhuǎn)動軸心與水面齊平，試驗時隨著拖車沿導(dǎo)軌在水池內(nèi)前后運動。此時模型全部浸沒在水池中，且其擺動過程中不會露出水面產(chǎn)生額外的興波阻力干擾。支撐和驅(qū)動作用的電動缸直線推桿會在水面以上滑行，但其產(chǎn)生的興波阻力相比于模型全沾濕受到的力而言比重較小，暫時可忽略其影響。

仿生胸鰭試驗?zāi)Ｐ屯弦匪販y力試驗包含兩部分運動，一部分是外部運動，一部分是自身運動。外部運動的產(chǎn)生源于水池拖車的帶動，沿導(dǎo)軌方向在水池內(nèi)前后運動。試驗?zāi)Ｐ妥陨淼倪\動依靠前、中、后鰭條的運動牽連產(chǎn)生，具體實現(xiàn)是：各鰭條在電動缸的帶動下，以鰭條根部為支撐點，左右擺動，實現(xiàn)沿展長方向的擺動；前、中、后鰭條在電動缸異步相位（通過設(shè)置啟動時間延遲實現(xiàn)）驅(qū)動的帶動下，實現(xiàn)鰭面的扭轉(zhuǎn)。通過改變電動缸的頻率、施加到鰭條上的驅(qū)動位移、異步相位差等參數(shù)，可實現(xiàn)仿生胸鰭試驗?zāi)Ｐ瓦\動頻率、擺動幅值和相位差的變化。測力天平傳感器感受仿生胸鰭試驗?zāi)Ｐ褪艿降碾妷鹤兓盘?，并通過電線、信號采集器輸送到計算機，完成試驗數(shù)據(jù)的采集與處理。

圖3 前中后鰭條運動定義Fig.3 Locomotion definition of front fin,middle fin and back fins

圖4 試驗?zāi)Ｐ桶惭bFig.4 Installation of experimental model

1.4 測力天平傳感器

試驗天平傳感器（如圖5 所示）采集頻率為200 Hz，三通道最大量程分別為：推力Fx測量通道1.5 kN、側(cè)向力Fy測量通道1.5 kN、扭矩Mx測量通道0.16 kN?m，單向測力精度為0.5%，串?dāng)_度誤差小于2%。

為了提高測力試驗平臺與仿生胸鰭試驗?zāi)Ｐ瓦B接的剛度，采用兩個三分力天平傳感器串聯(lián)。兩個天平傳感器前后共軸線布置，并與仿生胸鰭試驗?zāi)Ｐ娃D(zhuǎn)動軸線平行且共平面，以便于數(shù)據(jù)處理時力與力矩的解算。測力天平傳感器可以測量模型受到的水平軸向力和垂直于試驗?zāi)Ｐ亡捗娴膫?cè)向以及繞試驗?zāi)Ｐ娃D(zhuǎn)軸轉(zhuǎn)動的扭矩。試驗?zāi)Ｐ褪芰εc天平測量數(shù)據(jù)的關(guān)系為

式中，F(xiàn)x1、Fy1、Mx1、Fx2、Fy2、Mx2分別為測力天平1、測力天平2 采集的力和力矩信號；h 表示測力天平傳感器軸心到模型轉(zhuǎn)動軸的距離。

圖5 測力天平傳感器Fig.5 Force balance

1.5 試驗方法

仿生柔性胸鰭擺動拖曳水池試驗過程中，首先將拖車開至船塢，完成試驗?zāi)Ｐ团c拖車的固定連接（胸鰭尖端距離水池底部2 m 左右），檢查無誤后將拖車緩慢開至預(yù)定港位；然后電動缸上電，胸鰭運動控制發(fā)控系統(tǒng)（如圖6 所示）按照預(yù)定工況發(fā)出電動缸作動指令，電動缸單向啟動至絲杠的最大位置，延遲30 s后鰭條在電動缸的帶動下在鉛垂面內(nèi)產(chǎn)生周期性左右擺動，拖車延時60 s后啟動并緩慢加速至試驗速度；最后待拖車速度穩(wěn)定后，測力天平傳感器采集胸鰭受力的電壓信號，通過信號轉(zhuǎn)換獲得模型受的力和力矩，截取幾個周期內(nèi)的數(shù)據(jù)作為最終結(jié)果。通過控制電動缸輸出的頻率、輸出軸位移以及3個電動缸之間作動時間間隔獲得鰭條不同頻率、擺幅和相位差。

試驗過程中，拖車速度保持為0.4 m/s，對應(yīng)流動的雷諾數(shù)Re為5.16×105。

圖6 上位機控制界面Fig.6 PC control interface

2 試驗結(jié)果分析

2.1 仿生胸鰭運動形態(tài)分析

圖7給出了仿生胸鰭典型擺動工況下仿生柔性胸鰭模型的形狀變化圖像。仿生胸鰭在驅(qū)動系統(tǒng)的帶動下產(chǎn)生了左右擺動+鰭面扭轉(zhuǎn)耦合的主動柔性大變形，模型胸鰭尖端的薄壁區(qū)域受到水流的阻滯作用，發(fā)生了較為嚴重的不規(guī)則被動變形，此現(xiàn)象較好地反映了胸鰭的柔性特征。

圖7 水下仿生胸鰭擺動形態(tài)圖像Fig.7 Shape image of underwater bionic pectoral fin oscillation

2.2 流場形態(tài)分析

圖8 給出了仿生胸鰭典型擺動工況下尾流場形態(tài)。在尾流場中產(chǎn)生了明顯的尾渦，相鄰尾渦的旋向相反，并沿仿生胸鰭轉(zhuǎn)軸中心線等間距排布，對比圖8（b）中卡門渦街形態(tài)，仿生胸鰭擺動尾渦呈現(xiàn)與卡門渦街旋轉(zhuǎn)方向相反的結(jié)構(gòu)特征，可稱之為反卡門渦街形態(tài)，這也是仿生擺動推力產(chǎn)生的根源。分析系列的試驗圖像發(fā)現(xiàn)（見圖9～10），擺動頻率越大，尾流中相鄰渦核水平間距L 越??；擺動角度幅值越大，尾流中相鄰渦核縱向間距H越大。

圖8 尾渦流場試驗圖像Fig.8 Experimental image of vortex flow

圖9 不同頻率下渦核水平位置Fig.9 The horizontal distribution of vortex core region under different frequencies

圖10 不同擺動角度下渦核縱向位置Fig.10 The vertical distribution of vortex core region under different swing angle

2.3 水動力特性分析

2.3.1 水動力參數(shù)定義

此外，定義Ct_max、Ct_min 分別為推力系數(shù)周期性震蕩幅值的最大和最小值，Cn_max、Cm_max 分別為側(cè)向力系數(shù)、力矩系數(shù)周期性震蕩幅值的最大值（最小值與最大值的絕對值相等）。

2.3.2 運動頻率影響

圖11給出了仿生胸鰭模型不同擺動頻率下（分別為0.2、0.3、0.4 Hz，此時拖車速度、擺動角度和相位差均一致）的水動力特性變化曲線，圖中曲線的負值代表正推力（與運動方向相反，下同）。圖中結(jié)果表明：頻率越大，推力曲線Ct峰值和均值也越大（由負推力變?yōu)檎屏Γ?、?cè)向力曲線Cn和繞胸鰭轉(zhuǎn)動軸的力矩曲線Cm峰值也越大；仿生胸鰭擺動過程中柔性鰭面受到流體力的作用，胸鰭材料產(chǎn)生明顯的彈性變形（0.4 Hz工況下尤為明顯），導(dǎo)致整個胸鰭發(fā)生不規(guī)則形狀變化影響，推力曲線Ct在均值上下呈現(xiàn)不對稱現(xiàn)象，均值線以上Ct曲線時間行程較長，且頻率越大，這種非對稱現(xiàn)象越顯著；側(cè)向力曲線Cn和力矩曲線Cm受鰭面不規(guī)則形狀變化影響較小，呈現(xiàn)較好的正弦變化特征。

圖11 頻率對水動力參數(shù)影響曲線Fig.11 The hydrodynamic parameter curves under different frequencies

將各力系數(shù)在周期內(nèi)進行平均計算，獲得力系數(shù)周期均值隨柔性胸鰭運動頻率變化關(guān)系，如圖12所示。圖中，無量綱推力均值為正，代表此時試驗?zāi)Ｐ褪艿降乃畡恿樽枇Γㄏ峦?。隨著頻率的增大，柔性胸鰭受到的正向推力均值在減小，表示此時凈推力在增大，即推力均值隨頻率增大而近似線性增大。此外，各力和力矩系數(shù)周期振蕩的最大幅值隨頻率增大而近似線性增大。

圖12 頻率對水動力參數(shù)均值、幅值影響曲線Fig.12 Influence of frequency on the average and the amplitude of hydrodynamic parameters

2.3.3 擺動角度幅值影響

圖13給出了仿生胸鰭模型不同擺動角度下（此時拖車速度、擺動頻率和相位差均一致）水動力特性變化曲線。圖中曲線表明：擺動角度越大，推力曲線Ct峰值和均值也越大（由負推力變?yōu)檎屏Γ?、?cè)向力曲線Cn和繞胸鰭轉(zhuǎn)動軸的力矩曲線Cm峰值也越大；與水動力曲線受頻率影響變化特性類似，柔性鰭面不規(guī)則形狀變化引起了推力曲線Ct均值上下曲線呈現(xiàn)不對稱現(xiàn)象，擺動角度越大，非對稱現(xiàn)象越顯著；側(cè)向力曲線Cn和力矩曲線Cm受鰭面不規(guī)則形狀變化影響較小，呈現(xiàn)較好的正弦變化特征。

圖14給出了力系數(shù)周期均值隨柔性胸鰭擺動角度變化關(guān)系。隨著擺動角度的增大，柔性胸鰭受到的正向推力均值也在減小，表示此時凈推力在增大，即推力均值隨擺動角度增大而近似線性增大。此外，各力和力矩系數(shù)周期振蕩的最大幅值隨擺動角度的增大而近似線性增大。

圖13 擺動角度對水動力參數(shù)影響曲線Fig.13 Influence of swing angle on hydrodynamic parameters

圖14 擺動角度對水動力參數(shù)均值、幅值影響曲線Fig.14 Influence of swing angle on the average and the amplitude of hydrodynamic parameters

2.3.4 相位差影響

圖15給出了仿生胸鰭模型不同相位差下（此時拖車速度、擺動頻率和擺動角度一致）水動力特性變化曲線。圖中曲線表明：相位差越?。ù藭r對應(yīng)鰭面整體扭轉(zhuǎn)較小），推力曲線Ct峰值越大、側(cè)向力曲線Cn和繞胸鰭轉(zhuǎn)動軸的力矩曲線Cm峰值也越大；與水動力曲線受頻率影響變化特性類似，柔性鰭面不規(guī)則形狀變化引起了推力曲線Ct在均值上下呈現(xiàn)不對稱現(xiàn)象，均值線以上Ct曲線時間行程較長，相位差越小，非對稱現(xiàn)象越顯著；側(cè)向力曲線Cn和力矩曲線Cm受鰭面不規(guī)則形狀變化影響較小，呈現(xiàn)較好的正弦變化特征。

圖15 相位差對水動力參數(shù)影響曲線Fig.15 Influence of phase difference on hydrodynamic parameters

圖16給出了力系數(shù)周期均值隨柔性胸鰭運動相位差變化關(guān)系。隨著相位差角度的增大，柔性胸鰭受到的正向推力均值呈現(xiàn)先減小后增大現(xiàn)象，表示此時凈推力先增大后減小，即存在一個最佳相位差，使柔性胸鰭推力達到最大值。各力和力矩系數(shù)周期振蕩的最大幅值隨相位差變化減小。

圖16 相位差對水動力參數(shù)均值、幅值影響曲線Fig.16 Influence of phase difference on the average and the amplitude of hydrodynamic parameters

3 結(jié) 論

本文針對仿生柔性胸鰭推進模式下水動力特性研究，研制了仿生柔性胸鰭測力試驗?zāi)Ｐ图膀?qū)動控制系統(tǒng)，首次在大型拖曳水池環(huán)境下完成了仿生胸鰭擺動過程測力試驗，試驗測試了頻率、擺動角度、相位差等因素對仿生胸鰭模型水動力的影響，通過試驗結(jié)果的分析獲得如下結(jié)論：

（1）仿生柔性胸鰭模型擺動過程中受到水流的阻滯作用，模型胸鰭尖端的薄壁區(qū)域發(fā)生了較為嚴重的不規(guī)則被動變形，此現(xiàn)象較好地反映了胸鰭的柔性特征。

（2）仿生胸鰭擺動過程中，在尾流場中產(chǎn)生明顯的反卡門渦街，且擺動頻率越大，相鄰渦核水平間距就越??；擺動角度幅值越大，相鄰渦核縱向間距就越大。

（3）仿生胸鰭擺動頻率越大、擺動角度越大、胸鰭鰭條間的相位差越小，胸鰭產(chǎn)生的推力、側(cè)向力和轉(zhuǎn)動力矩峰值就越大；受仿生胸鰭擺動過程柔性鰭面不規(guī)則形狀變化的影響，推力曲線Ct在均值上下呈現(xiàn)不對稱現(xiàn)象，頻率越大、擺動角度越大、相位差越小，非對稱現(xiàn)象就越顯著；側(cè)向力曲線Cn、力矩曲線Cm受鰭面不規(guī)則形狀變化影響較小，呈現(xiàn)較好的正弦變化特征。

（4）頻率越大、擺動角度越大，柔性胸鰭凈推力就越增大，且推力均值、各力和力矩系數(shù)周期振蕩的最大幅值隨頻率、擺動角度增大而近似線性增大。

（5）相位差越大，柔性胸鰭凈推力先增大后減小，各力和力矩系數(shù)周期振蕩的最大幅值也減小?？梢酝茢?，存在一個最佳相位差，使柔性胸鰭推力達到最大值。