賀媛媛，楊炫，韓慧，王琦琛，張航

北京理工大學(xué) 宇航學(xué)院，北京 100081

微型飛行器現(xiàn)如今處于高速發(fā)展的階段，在軍事和民用上被運用得越來越廣泛。按照飛行器的布局形式與飛行原理可以將微型飛行器分為固定翼微型飛行器、旋翼微型飛行器和撲翼微型飛行器［1］。近年來，通過對鳥類及昆蟲等飛行生物的觀察、試驗以及數(shù)值仿真研究，仿生撲翼飛行器受到越來越多學(xué)者的關(guān)注［2-3］。相對于固定翼飛行器與旋翼飛行器，大量撲翼飛行器具有氣動效率高、機動性好、能垂直起降并保持懸停等優(yōu)點，但同樣存在機械系統(tǒng)復(fù)雜等缺點。針對撲翼飛行器的弊端，英國克蘭菲爾德大學(xué)的Guo 等［4］提出了撲旋翼的概念，結(jié)合撲翼與旋翼的技術(shù)特點，通過安裝一對中心對稱的翼使得撲翼在上下拍動的同時在兩側(cè)翼面上產(chǎn)生一對力偶，帶動撲翼實現(xiàn)旋轉(zhuǎn)，在保留撲翼優(yōu)點的同時使得結(jié)構(gòu)較為簡單，能耗降低。

目前，世界上微型撲翼飛行器與撲旋翼飛行器的研究主要集中在實現(xiàn)仿生拍動的驅(qū)動機構(gòu)設(shè)計方面［5-6］和對飛行器氣動特性進行的數(shù)值模擬和試驗驗證方面［7-13］。Keennon 等［14］對撲翼構(gòu)型進行研究，分析翼膜材料、撲翼展弦比等對于升力產(chǎn)生的影響。Chen 等［15］通過對于撲翼形態(tài)和慣性參數(shù)的改變探究其對于飛行器飛行性能的影響。周超等［16］基于撲旋翼飛行器的發(fā)展，提出了撲旋翼增升設(shè)計。而在對昆蟲翅翼的研究方面，蜻蜓翅翼由于薄且輕但在拍動、滑翔和懸停等飛行狀態(tài)下的穩(wěn)定性和承載能力都極其優(yōu)異因而受到研究者們的廣泛關(guān)注［17-19］。宏觀上，蜻蜓翅翼由翅脈和翅膜組成［20］。翅膜輕薄，為透明狀，使蜻蜓飛行時產(chǎn)生空氣動力，翅脈縱橫交錯，主要起到支撐作用。從微觀角度來看，翅膜的成分為結(jié)構(gòu)性蛋白，包括背部層、中間層和腹部層［21］。表面不同部位微觀結(jié)構(gòu)的不同有效減少了蜻蜓飛行過程中的阻力。此外，部分學(xué)者還對翅翼的材料和力學(xué)性能進行了研究。Kempf［22］使用納米壓痕技術(shù)發(fā)現(xiàn)翅脈的彈性模量為（2.9±0.8） GPa，翅膜的彈性模量為（1.5±0.5） GPa，且前翅的彈性模量和硬度大于后翅。張金［23］和潘春祥等［24］分別對蜻蜓翅翼的前緣脈、徑脈和后緣脈進行了納米力學(xué)性能研究，得到了蜻蜓翅翼翅脈彈性模量和納米硬度沿翅翼展向的變化關(guān)系。

在目前的電機電池等部件尺寸和重量的條件限制下，可實現(xiàn)飛行仿生撲翼微型飛行器的尺度和重量都比昆蟲大一個數(shù)量級，而到目前為止，仿生研究實際上都是在模仿昆蟲翅的拍動模式，對昆蟲翅本身而言，無論是按照幾何相似仿制昆蟲翅，或是制作簡化外形和結(jié)構(gòu)的撲翼，尚無統(tǒng)一而成熟的設(shè)計制作方法。本文基于自然界中不同尺度下昆蟲翅翼抗彎剛度與幾何形態(tài)參數(shù)的關(guān)系，并依據(jù)蜻蜓翅剛度實驗結(jié)果得到其展向抗彎剛度與展長的關(guān)系，建立了適用于仿蜻蜓翅的撲翼設(shè)計剛度相似性準則。首先采用拓撲優(yōu)化方法建立仿生撲翼構(gòu)形分布，進而根據(jù)剛度相似性準則確定具體結(jié)構(gòu)的參數(shù)，設(shè)計制作出仿蜻蜓翅撲翼，最后通過撲旋翼模型試驗驗證了基于剛度相似性研制的撲翼能夠明顯增升并提高效率。

1 基于主成分分析法的昆蟲翅翼剛度

1.1 主成分分析法

主成分分析法是一種多元統(tǒng)計方法，用于實現(xiàn)多個變量之間相關(guān)性的考察，基于對各個變量構(gòu)成的矩陣進行研究，分析其內(nèi)部關(guān)系，將多個存在相關(guān)性、闡述不同信息的指標轉(zhuǎn)換為盡可能少的不相關(guān)的綜合指標［25］。

主成分分析法的基本步驟如下：

1） 假定有n個樣本，每個樣本共有p個變量，構(gòu)成一個n×p階的數(shù)據(jù)矩陣：

2） 計算相關(guān)系數(shù)矩陣：

式中：rst為原變量xs與xt的相關(guān)系數(shù)（s，t=1，2，…，p），且rst=rts，其計算公式為

3） 計算相關(guān)系數(shù)矩陣R的特征根與特征向量，確定主成分。由式（4）所示的特征根方程式可求得p個特征根λs：

式中：I為單位矩陣。

按其大小順序排列為λ1≥λ2≥…≥λp≥0。由式（4）可知，每一個特征根對應(yīng)一個特征向量es。

4） 計算主成分貢獻率，根據(jù)累計貢獻率確定主成分個數(shù)。

主成分zs的貢獻率為

累積貢獻率為

主成分分析中在確保損失的信息量盡可能少的前提下選取盡量少的m個主成分進行綜合評價，通常選擇累計貢獻率的特征值對應(yīng)的第1，2，…，m(m≤p)個主成分。

5） 計算主成分載荷。

原變量xt在諸主成分z1，z2，…，zm(m≤p)上載荷的計算公式為

6） 記x1，x2，…，xp為原變量指標，各主成分z1，z2，…，zm為新變量指標。各主成分的得分為

1.2 昆蟲翅翼剛度分布

昆蟲翅翼在飛行中由于氣動力、慣性力和翼結(jié)構(gòu)彈性力的耦合作用產(chǎn)生被動的氣動彈性變形，雖然不同尺寸和類型的昆蟲翅外形和脈結(jié)構(gòu)有所差異，但整體具有相似的剛度分布特性。

部分昆蟲的翅長度、弦長、面積、昆蟲質(zhì)量和翅70% 展長處的展向截面抗彎剛度EI 數(shù)據(jù)如表1 所示［26］。從表1 中可以看出，隨著翅翼長度、翅翼弦長、翅翼面積和昆蟲質(zhì)量的增加，昆蟲翅翼展向抗彎剛度逐漸增大。為了得到各個參數(shù)的綜合指標對撲翼展向抗彎剛度的影響，采用主成分分析法進行分析，通過不同尺度下飛行生物翅翼剛度與幾何形態(tài)參數(shù)之間的關(guān)系確定剛度相似性準則。

表1 部分昆蟲翅尺寸和重量及翼展向抗彎剛度數(shù)據(jù)［26］Table 1 Size， weight and bending stiffness of some insect wings［26］

以昆蟲翅翼翅長、翅翼弦長、翅翼面積和昆蟲質(zhì)量為參數(shù)，探究4 者的綜合指標對撲翼展向抗彎剛度的影響。首先，取表1 前4 列的數(shù)據(jù)（翅長、弦長、翼面積、重量）構(gòu)成如式（1）所示的10×4 的數(shù)據(jù)矩陣，標準化處理翅長、弦長、翅翼面積和質(zhì)量的相關(guān)數(shù)據(jù)，通過主成分分析法計算得到相關(guān)系數(shù)矩陣為

通過計算得到特征值及各個主成分的方差貢獻率如表2 所示。由表2 可知，第1 主成分的累計貢獻率高于85%，故只需求出第1 主成分即可。對于特征值λ1=3.884 6，求出其特征向量，再計算得到各變量在主成分上的載荷，則主成分z1得分表達式為

表2 特征值及主成分貢獻Table 2 Eigenvalues and principal component contribution

采用對數(shù)函數(shù)對由翅長、弦長、翅翼面積和質(zhì)量組成的第1 主成分與昆蟲展向抗彎剛度的數(shù)據(jù)進行擬合，得到由翅長、弦長、翅翼面積和質(zhì)量組成的綜合指標z1與昆蟲展向抗彎剛度EI 之間的關(guān)系式為

2 蜻蜓翅剛度實驗

蜻蜓具有前飛、滑翔、倒飛、垂直飛行和懸停等多種飛行性能，翅翼拍動頻率低，飛行速度快，在有翼昆蟲中的綜合飛行性能最優(yōu)，但卻是仿生飛行器實現(xiàn)飛行最具有挑戰(zhàn)的模式。蜻蜓翅展弦比相對其他昆蟲大，在飛行中會發(fā)生大幅度的彎曲變形和扭轉(zhuǎn)變形，剛度分布和大小對飛行性能的影響大。以圖1 所示的黃蜻為例對蜻蜓翅翼的展向抗彎剛度進行進一步分析。

圖1 黃蜻Fig.1 Yellow dragonfly

蜻蜓翅樣件幾何形態(tài)參數(shù)測量結(jié)果如表3 所示。由于蜻蜓后翅較多采集樣件破損，因此只基于蜻蜓前翅進行剛度測試實驗。同一蜻蜓的2 個前翅的長度和最大弦長基本相同。

表3 蜻蜓幾何形態(tài)參數(shù)測量結(jié)果Table 3 Results of dragonfly geometric morphological parameters measurement

對蜻蜓前翅進行剛度實驗，采用懸臂梁彎曲試驗方法測量其抗彎剛度。測量黃蜻翅翼抗彎剛度方法如圖2 所示。將蜻蜓翅根粘接在固定支座上，再將質(zhì)量塊通過細線連接在翅脈上，測量加載荷點處的位移。

圖2 蜻蜓翅翼剛度測量Fig.2 Measurement of stiffness of dragonfly wings

根據(jù)測量的位移結(jié)果，結(jié)合式（12）可計算得出蜻蜓翅在展向不同截面的抗彎剛度。為提高精確度，同一位置反復(fù)施加不同大小載荷得到平均值，加載力的大小以線彈性變形為準則，加載力和位置及計算結(jié)果如表4～表9 所示。

表4 1 號蜻蜓右翅剛度實驗測量結(jié)果Table 4 Experimental measurements of right wing stiffness of Dragonfly 1

式中：F為翅所受的加載力；L為翅根到加載截面的距離；δ為加載點處的位移，如圖2（a）所示。

表4 和表5 的測量結(jié)果顯示，1 號蜻蜓左右兩翅在相同展向位置處的抗彎剛度存在較大差異，這是因為在實驗之前將1 號蜻蜓左右兩翅同時從蜻蜓上摘下，導(dǎo)致做完右翅實驗時，左翅隨著時間流逝變干，剛度增大，相同載荷下位移變小，這說明翅翼的剛度與生物體的活性有關(guān)，因此后續(xù)的實驗中沒有同時摘下翅翼，而是在做完一只翅翼的試驗后再摘下另一只翅翼。

表5 1 號蜻蜓左翅剛度實驗測量結(jié)果Table 5 Experimental measurements of left wing stiffness of Dragonfly 1

由于1 號蜻蜓左翅是在翅翼活性降低剛度增大的情況下測得的，誤差較大，所以將忽略表5 測得的1 號蜻蜓左翅剛度，不將其計入后續(xù)數(shù)據(jù)處理。根據(jù)表4 及表6～表9 所示的蜻蜓翅在不同展長處加載的位移實驗結(jié)果，可以得出如圖3（a）所示的蜻蜓翅不同展長處載荷大小與位移的變化關(guān)系。結(jié)果顯示，在相同展長位置施加載荷時，載荷與位移基本呈線性關(guān)系，說明實驗是在蜻蜓翅的線彈性變形范圍內(nèi)，滿足式（12）的條件。因此，根據(jù)式（12）計算可得表4 及表6～表9所示的各蜻蜓翅的剛度結(jié)果。取上述實驗結(jié)果的平均值，得到如圖3（b）所示的各蜻蜓翅在不同展向位置抗彎剛度。

圖3 蜻蜓翅剛度實驗結(jié)果Fig.3 Results of dragonfly wing stiffness experiments

表6 2 號蜻蜓左翅剛度實驗測量結(jié)果Table 6 Experimental measurements of left wing stiffness of Dragonfly 2

表7 2 號蜻蜓右翅剛度實驗測量結(jié)果Table 7 Experimental measurements of right wing stiffness of Dragonfly 2

表8 3 號蜻蜓左翅剛度實驗測量結(jié)果Table 8 Experimental measurements of left wing stiffness of Dragonfly 3

表9 3 號蜻蜓右翅剛度實驗測量結(jié)果Table 9 Experimental measurements of right wing stiffness of Dragonfly 3

由圖3（b）可見，蜻蜓翅翼的抗彎剛度從50%展長至90%展長處逐漸減小。目前的研究發(fā)現(xiàn)蜻蜓翅翼剛度分布從翼根到翼尖大致呈線性關(guān)系變化［27］，結(jié)合圖3（b）可以得出蜻蜓翅剛度與展長關(guān)系的擬合計算式：

式中：E為翅材料的彈性模量；Ii為撲翼展長i×10%處的截面慣性矩（i取0，1，…，10）；m為Ii沿展長變化的線性一次函數(shù)的斜率；n為線性一次函數(shù)的截距。

3 基于剛度相似性的撲翼結(jié)構(gòu)優(yōu)化設(shè)計

3.1 撲翼結(jié)構(gòu)拓撲優(yōu)化

蜻蜓翅翼上具有很多的微觀結(jié)構(gòu)，翅脈縱橫交錯，即便采用3D 打印加工方法也難以做到對其翅脈分布完全復(fù)制或仿制，也不能保證其剛度和質(zhì)量分布與蜻蜓翅的相似性，所以首先采用拓撲方法對仿生撲翼結(jié)構(gòu)布局進行優(yōu)化，對仿蜻蜓翅的撲翼結(jié)構(gòu)采用殼單元建立的有限元模型如圖4 所示。蜻蜓由于前緣脈是蜻蜓翅承力的主要部件，翅脈直徑相對大，是不可或缺的懸臂梁，因此將前緣脈位置（即圖中藍色區(qū)域）模型定義為非設(shè)計域，圖中紫色部分為結(jié)構(gòu)設(shè)計域。蜻蜓在實際飛行中的氣動力變化復(fù)雜，此處簡化為垂直于翅翼表面的均布載荷，其合力為蜻蜓體重的一半。另外，在翅翼根部設(shè)置六自由度約束，模擬蜻蜓翅與體的連接。

圖4 蜻蜓翅翼平面結(jié)構(gòu)有限元模型Fig.4 Finite element model of plane structure of dragonfly wing

撲翼需要具有一定的承載能力以承受飛行過程中產(chǎn)生的氣動載荷，這就對撲翼的剛度有所要求。此外，減重是撲翼設(shè)計一直致力的目標。因此，建立如式 （14）所示的結(jié)構(gòu)拓撲優(yōu)化數(shù)學(xué)模型，其中，設(shè)計變量為各個單元的偽密度，優(yōu)化目標為結(jié)構(gòu)節(jié)點力下所做的功，約束條件為材料保留比。

式中：xα為第α個單元的偽密度；C為結(jié)構(gòu)節(jié)點力F所做的功；F為結(jié)構(gòu)節(jié)點力向量；U為結(jié)構(gòu)節(jié)點位移向量；M為結(jié)構(gòu)質(zhì)量；M0為結(jié)構(gòu)原質(zhì)量；f為材料保留比，表示優(yōu)化后結(jié)構(gòu)保留的材料占原結(jié)構(gòu)材料的比例；K為結(jié)構(gòu)整體剛度矩陣；xmin為防止剛度矩陣奇異引入的單元偽密度最小值。其中，C值越小代表結(jié)構(gòu)剛度越大。

對于如圖4 所示的仿蜻蜓翅撲翼平面結(jié)構(gòu)有限元模型，利用ABAQUS 軟件中的Optimization模塊進行結(jié)構(gòu)拓撲優(yōu)化，插值模型罰系數(shù)取p=4，材料保留比分別取f=0.189～0.191，0.199～0.201，0.219～0.221，0.259～0.261，0.299～0.301，得到的結(jié)構(gòu)拓撲優(yōu)化結(jié)果如圖5 所示。

圖5 不同材料保留比設(shè)置下的結(jié)構(gòu)拓撲優(yōu)化結(jié)果Fig.5 Results of structural topology optimization with different settings of material retention ratio

從拓撲優(yōu)化結(jié)果可見，對應(yīng)不同材料保留比的撲翼結(jié)構(gòu)都具有較為清晰且與蜻蜓翅相似的脈絡(luò)分布，翼根處的脈絡(luò)密度高，材料保留量大，是由于該處應(yīng)力集中并承受全部載荷。在撲翼前緣梁處和后緣中部延伸出多條脈絡(luò)，將分布載荷傳遞到翼前梁。

在對撲翼結(jié)構(gòu)進行拓撲優(yōu)化過程中，由于將復(fù)雜的氣動力簡化為垂直于翅翼表面的均布載荷，忽略了前緣渦的增升作用，因此導(dǎo)致拓撲優(yōu)化的模型前緣脈絡(luò)細節(jié)缺失，與蜻蜓翅脈有差別。而蜻蜓翅具有更多細致分布的主脈和從脈，如圖6 所示，基于蜻蜓翅的脈絡(luò)分布，在拓撲優(yōu)化的有限元模型上添亞前緣脈、橫脈和中部主脈等缺失的細節(jié)，得到仿蜻蜓撲翼的基本結(jié)構(gòu)模型如圖7 所示。

圖6 蜻蜓翅翼樣本及脈絡(luò)分布Fig.6 Dragonfly wing samples and vein distribution

圖7 仿蜻蜓簡化翅翼設(shè)計模型Fig.7 Simplified design model of dragonfly wing

3.2 基于幾何相似性的撲翼設(shè)計

基于上述撲翼的結(jié)分布模型，可進一步針對各翼脈細節(jié)進行仿生設(shè)計分析。如圖8 所示，蜻蜓翅脈為環(huán)形空心結(jié)構(gòu)，根據(jù)文獻［28］中對蜻蜓翅的尺寸測量，取翅長度為47.6 mm，將翅脈截面簡化為外圍半徑為0.067 5 mm、厚度為0.025 mm 的圓環(huán)，翅膜厚度為0.004 mm。鑒于加工技術(shù)的限制，將該翅脈截面簡化為相同剛度的梁單元模型。根據(jù)計算得出，在材料參數(shù)不變的情況下，半徑為0.064 7 mm 的圓形截面梁單元模型具有和空心翅脈相同的剛度。

圖8 蜻蜓翅脈的橫截面及其幾何構(gòu)型［28］Fig.8 Cross section of dragonfly wing veins and its geometric configuration ［28］

進而可根據(jù)需要設(shè)計的微型仿生飛行器的撲翼尺度對仿蜻蜓翅撲翼模型按幾何相似比例放大。以半展長160 mm 的撲翼為例，翼展放大比例為3.36（=160/47.6），在原始模型中對每個截面的展弦比按照3.36 倍進行放大，按照原始模型脈絡(luò)走向與趨勢連接新的點形成放大后的翅翼脈絡(luò)，以此保持原始模型的翅翼外形結(jié)構(gòu)、脈絡(luò)分布不變，同時將模型中翅脈的截面半徑放大為0.217 5 mm（=0.064 7×3.36）?？紤]到加工因素，將翅脈截面設(shè)計為矩形，在保持截面積不變的情況下，選定截面尺寸為0.4 mm×0.37 mm，得到按幾何相似性設(shè)計的撲翼模型如圖9 所示。

圖9 基于幾何相似性放大的撲翼模型Fig.9 Flapping wing model based on geometrical similarity amplification

3.3 基于剛度相似性的撲翼設(shè)計

3.3.1 撲翼剛度分布

在設(shè)計和制作時需根據(jù)所采用的材料參數(shù)和尺寸限制按照剛度相似性準則確定撲翼模型中各梁元的剛度，以同樣的撲翼模型為例，其半展長、最大弦長、表面積和設(shè)計飛行器重量分別為：l=160 mm，b=40 mm，S=5 000 mm2，m=30。根據(jù)建立的剛度相似性準則式（10）及式（11），再結(jié)合式（13）得到該撲翼模型展向抗彎剛度隨展長的變化為

根據(jù)式（15）可計算出撲翼在不同展長處的抗彎剛度如表10 所示。

表10 撲翼沿展長的展向抗彎剛度Table 10 Bending stiffness along span of flapping wing

3.3.2 變截面撲翼參數(shù)優(yōu)化

由于蜻蜓翅脈的截面尺寸從翅根到翅尖，從前緣到后緣都逐漸減小，且沿展長的抗彎剛度按線性函數(shù)減小，因此設(shè)撲翼脈絡(luò)尺寸從翼根到翼尖，從前緣到后緣的變化也用線性函數(shù)描述。為了在滿足撲翼剛度要求的同時使結(jié)構(gòu)最輕，對撲翼進行優(yōu)化設(shè)計，優(yōu)化數(shù)學(xué)模型為

式中：lj=aj x+bj為各個梁的線寬函數(shù)，x為脈絡(luò)延伸方向上某點的位置，aj為線性一次函數(shù)的斜率，bj為線性一次函數(shù)的截距；Ii為撲翼展長i×10%處的截面慣性矩（i取0，1，…，10）；j為第j條脈絡(luò)；W為撲翼模型的質(zhì)量；lmin和lmax為表示線寬的上下邊界。

基于式（16）以撲翼質(zhì)量最小為優(yōu)化目標，各梁線寬函數(shù)的參數(shù)為變量，以表10 所示各展長處的剛度相似性和變量的上下邊界為約束條件，采用序列二次規(guī)劃方法和Matlab 中的Fmincon 工具箱對圖9 所示的仿生撲翼各條脈絡(luò)尺寸進行多次優(yōu)化計算，并取每次優(yōu)化結(jié)果平均值獲得最終結(jié)果。圖10 所示為優(yōu)化過程的迭代收斂曲線，在迭代次數(shù)達到約為15 次時，目標函數(shù)值開始收斂到接近一個常數(shù)（撲翼質(zhì)量約為1.1 g）。

圖10 撲翼質(zhì)量在優(yōu)化過程中的迭代收斂結(jié)果Fig.10 Iterative convergence results of flapping wing mass in optimization process

根據(jù)優(yōu)化結(jié)果得到7 條脈絡(luò)梁的線寬函數(shù)中的參數(shù)aj和bj，如表11 所示。根據(jù)表11 中的線寬函數(shù)參數(shù)aj和bj的優(yōu)化結(jié)果可以確定撲翼模型中7 條脈絡(luò)的線寬，線寬沿展向變化的結(jié)果如圖11 所示，結(jié)合圖12 所示的仿生撲翼模型的脈絡(luò)分布可見，撲翼各脈絡(luò)梁的線寬由翼根到末端、從前緣到后緣逐漸減小，脈絡(luò)l1和脈絡(luò)l3所處位置對應(yīng)蜻蜓翅翼的前緣脈和徑脈，由于前緣脈和徑脈是蜻蜓翅翼中承力的主要構(gòu)件，因此對于設(shè)計的撲翼這2 處梁的尺寸較大，符合自然界中飛行生物翅翼的結(jié)構(gòu)特點。

圖11 撲翼各脈絡(luò)線寬優(yōu)化結(jié)果Fig.11 Optimal results of flapping wing choroids width

圖12 基于幾何和剛度相似性的撲翼優(yōu)化模型Fig.12 Optimized flapping wing model based on geometry and stiffness similarity

表11 撲翼模型部分梁線寬函數(shù)的參數(shù)Table 11 Parameters of partial beam line-width function of flapping wing model

雖然上述剛度相似性是基于第2 節(jié)中昆蟲翅的展向剛度測量值獲得的，但基于蜻蜓翅脈的分布對仿生撲翼開展的結(jié)構(gòu)拓撲和優(yōu)化結(jié)果也部分包括了弦向剛度的相似性。前期對相近尺寸的撲旋翼的展向和弦向剛度與升力關(guān)系的試驗研究結(jié)果表明［29］，在展向剛度不變的情況下，撲旋翼的平均升力隨著弦向剛度的增加而增加，例如展向剛度與弦向剛度比值為30，對應(yīng)的升力為35.7 gf；當展向與弦向剛度比值減小到12 時，平均升力提高36%。隨著弦向剛度增加或隨著展向剛度與弦向剛度的比值減小，平均升力呈繼續(xù)增大的趨勢。表11 和圖11 所示的結(jié)果表明，優(yōu)化的仿生撲翼的展向和弦向剛度比小于12，例如在70% 展長與弦長處的展向剛度（2.418×103N·mm2）和弦向剛度（6.786×102N·mm2）比值大致為4，符合升力隨弦向剛度增加的趨勢。由于撲旋翼的最優(yōu)弦向剛度與拍動升力和旋轉(zhuǎn)推力的比例相關(guān)，換言之，最優(yōu)弦向剛度需根據(jù)一個撲旋翼的設(shè)計目標是最大升力還是最大效率而決定。

4 實驗與結(jié)果

4.1 升力測試實驗

4.1.1 撲翼升力測試系統(tǒng)

采用課題組研制出的微型撲旋翼飛行器作為實物測試模型，軸對稱安裝在模型兩側(cè)的撲翼通過齒輪和連桿傳動機構(gòu)在電機的驅(qū)動下實現(xiàn)上下拍動，在拍動過程中產(chǎn)生方向相反的推力使撲翼同時繞中心軸產(chǎn)生旋轉(zhuǎn)運動。模型試驗測試平臺由壓力傳感器、信號放大器、直流電源、NI數(shù)字采集系統(tǒng)、固定夾具、計算機組成，測試平臺原理示意圖及實物圖分別如圖13 和圖14所示。

圖13 升力測試系統(tǒng)平臺示意圖Fig.13 Schematic diagram of lift test platform

圖14 撲旋翼升力測試系統(tǒng)實驗平臺Fig.14 Experimental platform of flapping wing rotor lift test system

由夾具固定在壓力傳感器上方的撲旋翼模型在試驗中將沿豎直方向上產(chǎn)生的力傳輸至壓力傳感器（SKU-314990000），壓力傳感器將力信號轉(zhuǎn)換為電信號且傳輸至信號放大器（RWST01A）放大并由NI 數(shù)字采集卡（NI USB-6009，采樣率1 000 Hz）采集數(shù)據(jù)，最后傳輸至計算機對數(shù)據(jù)進行處理與存儲。

升力測試實驗流程包括：① 撲翼稱重；② 傳感器標定；③ 撲旋翼實驗樣機組裝；④ 升力實驗數(shù)據(jù)測試；⑤ 升力實驗數(shù)據(jù)與結(jié)果分析。

4.1.2 撲翼模型制作

撲翼由骨架和翼膜組成，其中骨架按照圖9和圖12 的設(shè)計結(jié)果采用碳纖維板機加切割制成，翼膜選擇厚度為0.015 mm 的聚酰亞胺薄膜，通過膠水粘連到撲翼骨架上。將圖15（a）所示基于幾何相似性設(shè)計制作的撲翼命名為JH 翼，圖15（b）所示兼顧幾何和剛度相似性制作的撲翼命名為GD 翼。

圖15 仿蜻蜓撲翼實驗樣件Fig.15 Test samples of dragonfly flapping wing test

為了將仿生撲翼與目前撲翼飛行器研究中通常采用的矩形撲翼進行對比，設(shè)計制作了一對矩形撲翼。矩形撲翼的面積和展弦比與仿生撲翼一致，骨架由主梁、次梁和斜次梁構(gòu)成，如圖16（a）所示。由表10 可知，剛度沿展向線性變化，于是設(shè)置翼根處與70%展長處的剛度與仿生撲翼相同并滿足表10 中剛度要求。各梁具體尺寸見表12。矩形撲翼的骨架采用碳纖維桿制成，翼膜采用與仿生撲翼相同的薄膜，通過膠水粘接得到如圖16（b）所示的試驗件實物，將該矩形翼試驗件命名為JX 翼。

圖16 矩形撲翼實驗?zāi)Ｐ虵ig.16 Test sample of rectangular flapping wing

表12 矩形翼模型的各梁尺寸Table 12 Dimensions of each beam of rectangular flapping wing model

4.2 升力實驗結(jié)果與分析

4.2.1 GD 翼升力隨電壓的變化

為明確GD 翼在不同電壓下的升力表現(xiàn)，進行了不同電壓下GD 翼的升力測試研究，表13和圖17 為實驗得到的GD 翼的升力測試數(shù)據(jù)及變化規(guī)律。實驗結(jié)果表明，GD 翼的升力隨著電壓增加而提高。當電壓較低時，最大升力和最大負升力的數(shù)值比較接近，隨著電壓的升高，最大升力的增長趨勢較快，最大負升力的增長趨勢相對較緩，導(dǎo)致平均升力增加。這是由于GD 翼剛度沿展向減小，隨著電壓和頻率的提高，撲翼在氣動力作用下發(fā)生扭轉(zhuǎn)，減小了上拍過程中產(chǎn)生的負升力，符合自然界中昆蟲的運動方式。在電壓為3.4 V 時升力達到了14.7 gf，接近文中用于試驗的仿生撲旋翼模型的重量，可實現(xiàn)垂直起飛。

圖17 GD 翼平均升力隨電壓變化結(jié)果Fig.17 Increase of average lift with input voltage

表13 不同電壓下GD 翼升力測試結(jié)果Table 13 Results of lift of GD wing at different voltages

4.2.2 GD 翼與JH 翼的升力對比實驗

為了探究撲翼剛度對升力的影響，將JH 翼的升力測試結(jié)果與表13 的GD 翼升力結(jié)果進行對比。由于JH 翼僅根據(jù)幾何相似性設(shè)計制作而成，忽略了剛度相似性準則，且剛度相對小，為獲得穩(wěn)定的測試數(shù)據(jù)，將試驗的拍動頻率定義在較低范圍進行。JH 翼在不同拍動頻率下的升力測試結(jié)果如表14 所示，二者在不同拍動頻率下的升力變化情況如圖18 所示。

圖18 GD 和JH 翼在不同拍動頻率時產(chǎn)生的平均升力Fig.18 Average lift of GD wing and JH wing at different flapping frequency

表14 JH 翼升力測試結(jié)果Table 14 Test results of lift from JH wing

由圖18 可知，GD 翼的平均升力在4.5～9.5 Hz 的低頻拍動范圍內(nèi)隨拍動頻率的提高呈線性增加，與之不同，JH 翼的平均升力在4.5～6.6 Hz 低頻拍動時隨拍動頻率呈線性提高，且與GD 翼的結(jié)果接近，但在高于6.6 Hz 的拍動頻率時升力急劇下降，這是因為JH 翼的剛度偏小，其彈性變形隨著拍動頻率和相應(yīng)的慣性力和氣動力的提高而超出了承受范圍。由表14 可見，盡管JH 翼在9.3 Hz 拍動頻率時產(chǎn)生的正負力峰值均大于6.6 Hz 時的結(jié)果，但得到的平均升力卻明顯減小，而基于幾何和剛度相似性設(shè)計的GD 撲翼所產(chǎn)生的升力比僅基于幾何相似性設(shè)計的JH 撲翼高出2 倍以上（表15）。

表15 不同拍動頻率下GD 翼與JH 翼的平均升力及提升比Table 15 Average lift and increasing ratio of GD wing and JH wing at different flapping frequency

4.2.3 GD 翼與JX 翼的升力對比實驗

明確了兼顧幾何和剛度相似性的GD 翼升力特性優(yōu)于僅考慮幾何相似性的JH 翼，為了進一步驗證仿生撲翼形狀和結(jié)構(gòu)對提高升力的影響，采用前文設(shè)計制作的JX 翼進一步開展了實驗研究，對比GD 翼與JX 翼在試驗條件下的升力表現(xiàn)。表16 為JX 翼的升力測試結(jié)果，GD 翼和JX翼平均升力的變化情況如圖19 所示。

圖19 GD 翼和JX 翼在不同拍動頻率時產(chǎn)生的平均升力Fig.19 Average lift of GD wing and JX wing at different flapping frequency

表16 JX 翼升力測試結(jié)果Table 16 Test results of lift from JX wing

由圖19 可知，GD 翼和JX 翼的平均升力隨著拍動頻率的增加都以近似線性變化增加。GD翼的平均升力顯著高于JX 翼，且由表17 可知，在12.5 Hz 的拍動頻率時比JX 翼的平均升力提高了25.2%，說明仿生撲翼形狀與結(jié)構(gòu)對撲旋翼的氣動力起到了有效的增升作用。但本研究限于氣動模擬方法對幾何外形進行優(yōu)化，僅限于與傳統(tǒng)矩形翼為初始設(shè)計開展剛度相似性的研究，所以不能對該仿蜻蜓翅幾何外形即為最優(yōu)外形下結(jié)論。

表17 不同拍動頻率下GD 翼與JX 翼的平均升力及提升比Table 17 Average lift and increasing ratio of GD wing and JX wing at different flapping frequency

5 結(jié) 論

基于主成分分析法獲得的昆蟲翅幾何參數(shù)與抗彎剛度的關(guān)系，通過蜻蜓翅的剛度實驗數(shù)據(jù)可建立用于仿生撲翼設(shè)計的剛度相似性準則。基于拓撲方法和幾何相似性設(shè)計制作了仿蜻蜓JH 撲翼，進而采用剛度相似性準則確定仿蜻蜓JH 翼的結(jié)構(gòu)參數(shù)，設(shè)計制作了GD 撲翼。為了深入對比分析，僅基于剛度相似性設(shè)計制作了普遍采用的由主次梁構(gòu)成的矩形JX 撲翼。最后通過撲旋翼模型試驗測試了3 種撲翼的升力表現(xiàn)，結(jié)果對比驗證了仿蜻蜓GD 翼具有顯著的氣動力增升效果，并得到如下結(jié)論：

1） 昆蟲翅脈結(jié)構(gòu)分布和幾何相似性可作為仿生撲翼設(shè)計的必要條件，但不是充分條件。因為結(jié)構(gòu)剛度不足的撲翼產(chǎn)生的升力僅在低頻拍動時隨頻率提高而線性增加，但升力也偏?。划斉膭宇l率提高到一定值后（JH 翼大于6.6 Hz），升力反而隨頻率的提高而明顯下降。換言之，撲翼的剛度應(yīng)根據(jù)飛行器設(shè)計的升力要求而定。

2） 基于昆蟲翅的剛度相似性也可作為仿生撲翼設(shè)計的必要條件，但仍不充分。如采用通常的矩形撲翼簡化結(jié)構(gòu)，基于剛度相似性確定其主次梁參數(shù)后制作的JX 翼，在拍動頻率為9.3 Hz時產(chǎn)生的升力是僅基于幾何相似性設(shè)計的JH 翼的4 倍。雖然JX 翼相比JH 翼在拍動頻率和升力上都可成倍提高，但仍不夠理想，這也是迄今為止仿生撲翼優(yōu)化面臨的問題。

3） 將昆蟲翅脈分布幾何相似性和結(jié)構(gòu)剛度相似性相結(jié)合的準則可作為仿生撲翼設(shè)計的必要和充分條件。試驗結(jié)果表明，基于這一準則設(shè)計制作的GD 撲翼在5～15 Hz 拍動頻率范圍內(nèi)測試的升力都明顯優(yōu)于僅基于幾何相似性或剛度相似性的JH 翼和JX 翼。例如，在拍動頻率12.5 Hz 時GD 翼產(chǎn)生的平均升力比矩形JX 翼提高了25.2%，即在滿足仿生幾何外形和仿生翅剛度相似性的條件下GD 翼是最優(yōu)設(shè)計。

4） 實驗驗證了撲旋翼結(jié)構(gòu)和剛度以及仿生撲翼形狀對升力有顯著的影響，根據(jù)幾何和剛度相似性制作GD 翼產(chǎn)生的升力隨輸入電壓的提高而增加，在電壓達到3.4 V 時產(chǎn)生的升力達到14.7 gf，接近用于試驗的仿生撲旋翼模型的重量，可實現(xiàn)垂直起飛。