王云飛，肖 洪，楊 廣，蔣華毅，劉榮強(qiáng)，丁江民 ，郭宏偉

(1.大連交通大學(xué) 機(jī)械工程學(xué)院，大連 116028；2.機(jī)器人技術(shù)與系統(tǒng)國(guó)家重點(diǎn)試驗(yàn)室(哈爾濱工業(yè)大學(xué))，哈爾濱 150001)

長(zhǎng)期以來(lái)傳統(tǒng)機(jī)翼是一種折衷方案，雖然飛機(jī)在多種工況下飛行，但每個(gè)工況的性能都是次優(yōu)的。如機(jī)翼能夠在起飛和降落時(shí)采用小后掠角、大展相比和高相對(duì)厚度的高升力外形，同時(shí)在超聲速巡航和機(jī)動(dòng)狀態(tài)采用大后掠角、小展弦比和低相對(duì)厚度的小阻力外形，使得氣動(dòng)性能最優(yōu)，從而擴(kuò)大飛行包線(xiàn)，適應(yīng)不同的飛行環(huán)境，滿(mǎn)足不同的任務(wù)需求，減少能耗，提高飛行效率[1-3]。

機(jī)翼變形方式最早可追溯到萊特兄弟利用繩索滑輪機(jī)構(gòu)控制和調(diào)節(jié)飛機(jī)飛行。隨著飛行速度與任務(wù)適應(yīng)性的不斷提高，設(shè)計(jì)出升降舵、副翼以及襟翼等小尺度局部變形方式來(lái)改善氣流[4]。為解決軍用飛機(jī)不同設(shè)計(jì)點(diǎn)氣動(dòng)布局矛盾，改善多功能性，提高作戰(zhàn)功能，出現(xiàn)了伸縮、折疊、后掠等這種對(duì)機(jī)翼整體改變的大尺度變形，如F-14和B-1等戰(zhàn)斗機(jī)采用變后掠機(jī)翼技術(shù)。然而這些改變遠(yuǎn)沒(méi)有達(dá)到高效率飛行的目的，一方面剛性部件會(huì)破壞機(jī)翼表面光滑性，產(chǎn)生額外阻力，另一方面，剛性固定機(jī)翼從根本上難以滿(mǎn)足人類(lèi)對(duì)于飛機(jī)性能的要求。研究人員嘗試將仿生結(jié)構(gòu)應(yīng)用到機(jī)翼變形上，如旺超等[5]研究了翼型厚度和翼型彎度對(duì)飛撲翼氣動(dòng)性能的影響，為變翼型的氣動(dòng)設(shè)計(jì)提供了理論指導(dǎo)和技術(shù)參考。Luca等[6]通過(guò)模擬翅膀改變機(jī)翼形狀，提高了低速機(jī)動(dòng)、高速抗風(fēng)性能，并采用不對(duì)稱(chēng)變形提高了滾動(dòng)控制，然而存在機(jī)構(gòu)鉸鏈脆弱，難以等比例縮小等問(wèn)題。Jenett等[7]研制了基于超輕彈性體單元的離散分布、可逆裝配的主動(dòng)扭轉(zhuǎn)變形翼機(jī)構(gòu)，這種離散式多胞元結(jié)構(gòu)實(shí)現(xiàn)了機(jī)翼的輕量化。Andersen等[8]提出了滑動(dòng)蒙皮方案，能實(shí)現(xiàn)面積與后掠角的單獨(dú)變化，創(chuàng)新性的為變形機(jī)翼提供了新思路；然而，其變面積與變后掠的翼肋不連續(xù)，不適用較大樣機(jī)，并且翼根強(qiáng)度較小，機(jī)翼剛度受氣動(dòng)影響較大。變后掠機(jī)翼結(jié)構(gòu)的分布式驅(qū)動(dòng)的位置與數(shù)目?jī)?yōu)化對(duì)飛行效率具有很大影響，最具有代表性的為Joo等[9-10]研究了變形結(jié)構(gòu)中驅(qū)動(dòng)器的優(yōu)化布置，得出驅(qū)動(dòng)器在結(jié)構(gòu)中的最佳位置。Johnson等[11-12]采用遺傳算法和基于梯度的優(yōu)化方法對(duì)不同工況下驅(qū)動(dòng)器的位置進(jìn)行優(yōu)化。朱華等[13]提出了利用超聲電機(jī)來(lái)驅(qū)動(dòng)小型變體飛行器變形所要研究的關(guān)鍵問(wèn)題，目前只適用于機(jī)翼局部變形。于曉至[14]提出了新型分布式結(jié)構(gòu)變形翼，并對(duì)其進(jìn)行了非線(xiàn)性動(dòng)力學(xué)建模與控制器研究，研究基于二、三單元，尚缺少整翼的多單元以及實(shí)驗(yàn)驗(yàn)證，沒(méi)有應(yīng)用電機(jī)實(shí)現(xiàn)結(jié)構(gòu)驅(qū)動(dòng)。在機(jī)構(gòu)變形方面，Xi等[15]提出了一種冗余驅(qū)動(dòng)的單自由度系統(tǒng)來(lái)解決空間和負(fù)載限制的問(wèn)題，Moosavian等[16]使用參數(shù)化曲線(xiàn)，通過(guò)單輸入的多環(huán)連桿機(jī)構(gòu)來(lái)改變翼型，在一定程度上減緩了機(jī)翼增重，是變構(gòu)型在機(jī)翼上的一次創(chuàng)新應(yīng)用，然而并未討論在翼展方向如何布置機(jī)構(gòu)。上述研究的變形翼多數(shù)結(jié)構(gòu)復(fù)雜，質(zhì)量大，承載能力弱，無(wú)法實(shí)現(xiàn)多維度的大變形，并且變弦長(zhǎng)時(shí)受到剛性翼肋的限制，大多是采用翼肋自身變形這種復(fù)雜結(jié)構(gòu)，同時(shí)驅(qū)動(dòng)略顯笨重單一，影響飛機(jī)驅(qū)動(dòng)效率，難以適應(yīng)大變形需求。

本文設(shè)計(jì)的變形翼要滿(mǎn)足大變形、高剛度，兼顧高、低速飛行工況，并需要結(jié)合所給定變形結(jié)構(gòu)的各項(xiàng)參數(shù)，對(duì)放置其內(nèi)的驅(qū)動(dòng)器位置做出合適優(yōu)化，試圖找到驅(qū)動(dòng)器最優(yōu)位置以及此時(shí)蒙皮初始狀態(tài)，同時(shí)探索分布式驅(qū)動(dòng)與機(jī)翼結(jié)構(gòu)剛度之間的關(guān)系，使得驅(qū)動(dòng)效率能夠達(dá)到最大。

1 平行連桿變形翼設(shè)計(jì)

1.1 變形原理

機(jī)翼是基于二維翼型建立的三維幾何體，機(jī)翼變形最大難點(diǎn)是需要考慮機(jī)翼平面與翼型兩個(gè)維度的協(xié)調(diào)。變形翼由主動(dòng)變形的骨架與柔性剪切變形蒙皮組成，可實(shí)現(xiàn)低速高升力狀態(tài)和高速低阻力狀態(tài)，如圖1所示。

(a)變形翼組成

剪切變形蒙皮[17]由鑲嵌在硅橡膠中的碳纖維棒和凱夫拉線(xiàn)固化而成，在骨架支撐下，維持氣動(dòng)外形。取x軸正方向?yàn)橐碚狗较颍瑈軸正方向在機(jī)翼平面內(nèi)垂直于x軸方向指向飛行方向，xy平面內(nèi)平行連桿構(gòu)成的多單元變后掠翼以O(shè)點(diǎn)為轉(zhuǎn)軸，通過(guò)內(nèi)置在四邊形的線(xiàn)性驅(qū)動(dòng)器改變夾角或者直線(xiàn)驅(qū)動(dòng)O3，使得翼根的O1、O2、O3上下滑動(dòng)來(lái)改變后掠角、機(jī)翼面積、展長(zhǎng)、弦長(zhǎng)等，虛線(xiàn)為變形后的高速狀態(tài)。yz平面內(nèi)，初始翼型為低速翼型，在變形時(shí)，沿任意展向翼剖面翼型均在同一剖面上變化，如，翼肋上的B1、A1、C1一直在一條線(xiàn)上，變形后翼型只在弦長(zhǎng)方向伸長(zhǎng)。

1.2 機(jī)翼幾何參數(shù)建模與分析

評(píng)價(jià)機(jī)翼變形帶來(lái)的氣動(dòng)效益，最主要看衡量機(jī)翼氣動(dòng)外形的幾何參數(shù)如何變化。B1B2=l1,B1O=l2，定義K=l1/l2為順氣流優(yōu)勢(shì)，O1B2與翼展方向的夾角為后掠角，記為θ，翼肋與來(lái)流方向的夾角為α=sin-1(Kcosθ)，翼根弦長(zhǎng)為c，半翼展為l/2，機(jī)翼面積為S(單位面積，下同)。展弦比為λ=l2/S=6.45Kcos2θ/cos(α-θ)，相對(duì)厚度為tmax/c。

相比旋轉(zhuǎn)式變后掠的翼肋與氣流方向夾角增大，平行連桿結(jié)構(gòu)變后掠時(shí)最大特點(diǎn)是翼肋做趨向于“順氣流擺動(dòng)”，即：隨著后掠角的增大，翼肋與氣流方向的夾角變小，如圖2所示。

(a)順氣流優(yōu)勢(shì)對(duì)翼肋擺動(dòng)的影響

K越小，翼肋擺動(dòng)程度越小，使得翼型趨向于初始最優(yōu)翼型，這對(duì)需要在高亞音速巡航的飛機(jī)來(lái)說(shuō)具有一定的工程意義。然而，“順氣流擺動(dòng)”的夾角不能無(wú)限變小，其一，在結(jié)構(gòu)上翼肋會(huì)無(wú)限增長(zhǎng)，其二，展弦比變化率會(huì)減小，為滿(mǎn)足翼肋趨向于順氣流方向以及大展弦比變化的寬速域，使得機(jī)翼在低、高速時(shí)均有較好的效率，文中折中選取K=0.75。

后掠角θ在20.0°～70.0°變化時(shí)，展弦比由4.8減小到1.0，相對(duì)厚度由13.4%減小到7.76%;后掠角θ在20.0°～53.4°時(shí)，機(jī)翼面積整體均要大于θ在53.4°～70.0°時(shí)的值;當(dāng)θ≈36.7°時(shí)，機(jī)翼面積達(dá)到最大值，面積變化率達(dá)到175%；小后掠角對(duì)應(yīng)大展弦比和高相對(duì)厚度，帶來(lái)高升力，適合爬升；大后掠角對(duì)應(yīng)小展弦比和低相對(duì)厚度，降低阻力，適合高速機(jī)動(dòng)。

2 變形翼分布式驅(qū)動(dòng)配置

2.1 驅(qū)動(dòng)器在單元中的位置優(yōu)化

2.1.1 單元參數(shù)建模

放置在變形翼內(nèi)部的線(xiàn)性驅(qū)動(dòng)器的尺寸、位置會(huì)影響到機(jī)翼驅(qū)動(dòng)效率，因此需要對(duì)驅(qū)動(dòng)器在單元中的位置合理優(yōu)化[9-10]。

平行連桿式變形翼為單自由度結(jié)構(gòu)，只對(duì)單個(gè)單元進(jìn)行研究，如圖3所示。驅(qū)動(dòng)器主要受外載荷和柔性剪切蒙皮的作用，蒙皮跟隨骨架變形時(shí)，碳纖維棒與凱夫拉線(xiàn)擠壓硅橡膠，與彈簧的作用相似，因此將蒙皮對(duì)骨架的作用等效為OB2、AB1兩個(gè)方向的彈簧，彈性系數(shù)分別為：Kd、Ky，彈簧OB2與A1B2夾角為β。該單元結(jié)構(gòu)所受的外載荷在x軸和y軸方向的投影分別為Fex、Fey；驅(qū)動(dòng)器的位置在單元結(jié)構(gòu)中由兩個(gè)向量Z1和Z2來(lái)確定。其中驅(qū)動(dòng)器的長(zhǎng)度和角度用向量Z3表示，驅(qū)動(dòng)器伸長(zhǎng)方向與x軸的夾角為ψ。

圖3 單元結(jié)構(gòu)在外載荷和驅(qū)動(dòng)力下的模型簡(jiǎn)圖Fig.3 Schematic diagram of a basic structure unit of morphing wing under external load and actuator force

驅(qū)動(dòng)器的長(zhǎng)度為

(1)

因此驅(qū)動(dòng)器伸長(zhǎng)的方向與x軸正方向的夾角為

(2)

(3)

在上述所有參數(shù)中，后掠角θ、Fex、Fey和K在求解過(guò)程中已知，以Z1為變量進(jìn)行優(yōu)化求解。

2.1.2 蒙皮等效剛度測(cè)定

定義初始后掠角θ0時(shí)蒙皮處于初始狀態(tài)，此時(shí)蒙皮不對(duì)骨架單元產(chǎn)生張拉作用。如圖4所示，當(dāng)蒙皮初始狀態(tài)為碳纖維棒與凱夫拉線(xiàn)的夾角ε=90°時(shí)，初始后掠角θ0≈36.7°，此時(shí)單元結(jié)構(gòu)剛好為矩形，以此狀態(tài)進(jìn)行蒙皮剛度拉伸試驗(yàn)。在靜加載試驗(yàn)機(jī)上通過(guò)拉伸帶有蒙皮的單元框架，得到加載載荷Fload與對(duì)角線(xiàn)位移x之間的關(guān)系，為了簡(jiǎn)化，等效彈簧的Kd(Ky)值應(yīng)該滿(mǎn)足：kminx≥Fload，得到擬合直線(xiàn)的斜率即為等效彈簧的剛度值，由于是單對(duì)角線(xiàn)受拉，故取一半kmin值，即為柔性剪切變形蒙皮等效彈簧剛度Kd=Ky≈kmin/2=0.625 N/mm。

(a)不同蒙皮初始狀態(tài) (b)測(cè)量實(shí)驗(yàn)

2.1.3 驅(qū)動(dòng)器力的確定

機(jī)翼在變形時(shí)，柔性蒙皮除了初始狀態(tài)，整個(gè)過(guò)程均對(duì)骨架產(chǎn)生作用力，單元在蒙皮、驅(qū)動(dòng)器和外載荷的作用下，通過(guò)準(zhǔn)靜態(tài)力學(xué)分析方法得到驅(qū)動(dòng)力的表達(dá)式：

(4)

其中：

式(4)說(shuō)明了驅(qū)動(dòng)力的大小取決于蒙皮初始狀態(tài)、蒙皮剛度、外載荷Fex、Fey的大小、驅(qū)動(dòng)器的位置Z1、驅(qū)動(dòng)器的方向Ψ以及可變后掠角θ。當(dāng)外載荷作用與輸出方向相反時(shí)，驅(qū)動(dòng)器的力大于Fine時(shí)才能實(shí)現(xiàn)該變形結(jié)構(gòu)的運(yùn)動(dòng)。

為此研究出不同后掠角時(shí)驅(qū)動(dòng)力Fine與Z1的關(guān)系，以便在給驅(qū)動(dòng)力賦值時(shí)滿(mǎn)足Fin≥Fine。定義驅(qū)動(dòng)器初始長(zhǎng)度為10(單位長(zhǎng)度，下同)，F(xiàn)ex=Fey=10(單位力，下同)。如圖5所示，驅(qū)動(dòng)力負(fù)值代表驅(qū)動(dòng)力方向與伸長(zhǎng)方向一致。

圖5 不同后掠角時(shí)驅(qū)動(dòng)力與設(shè)計(jì)參數(shù)Z1的關(guān)系Fig.5 Relationship between actuator force and design parameter Z1 under different sweep angles

驅(qū)動(dòng)器的位置Z1距離單元變形結(jié)構(gòu)的鉸接(B1或B2)越近時(shí)，所需要的驅(qū)動(dòng)力就越大，距離鉸接處越遠(yuǎn)驅(qū)動(dòng)力越??；隨著后掠角的增大，驅(qū)動(dòng)器所需最小驅(qū)動(dòng)力(minFine)的Z1值在減小，因此預(yù)估Z1的最優(yōu)值在5～8范圍內(nèi)。

2.1.4 變形效率的描述

在整個(gè)變形范圍內(nèi)，結(jié)構(gòu)在驅(qū)動(dòng)器、蒙皮和外載荷作用下，所做的變形量越大效率越高。定義變形后的角度為θf(wàn)，根據(jù)虛功原理，θf(wàn)滿(mǎn)足

(5)

定義單元結(jié)構(gòu)輸出位移為：

Xout=2l1(cosθ0-cosθf(wàn))

(6)

Yout=l1(sinθf(wàn)-sinθ0)+l2(cosαf-cosα0)

(7)

變形量L=Xout+Yout，建立驅(qū)動(dòng)器位置優(yōu)化的目標(biāo)函數(shù)與約束條件為：

(8)

通過(guò)matlab優(yōu)化工具箱fmincon函數(shù)可以計(jì)算出任意初始后掠角時(shí)的驅(qū)動(dòng)器最優(yōu)位置。

2.1.5 最優(yōu)位置求解

在求解最優(yōu)位置之前，進(jìn)行實(shí)驗(yàn)來(lái)驗(yàn)證理論模型的準(zhǔn)確性，并將測(cè)試結(jié)果與理論分析進(jìn)行比較。實(shí)驗(yàn)裝置如圖6所示，采用線(xiàn)性驅(qū)動(dòng)器驅(qū)動(dòng)帶有柔性剪切蒙皮的鉸鏈單元框架來(lái)測(cè)定變形量L，單元框架尺寸為270 mm×360 mm，驅(qū)動(dòng)器初始長(zhǎng)度為270 mm，以后掠角36.7°為初始狀態(tài)，即ε=90°。驅(qū)動(dòng)力Fine=500 N，測(cè)得單元框架的變形量L，重復(fù)4次取平均值，調(diào)整驅(qū)動(dòng)器的位置Z1，在多個(gè)位置進(jìn)行測(cè)試，得到驅(qū)動(dòng)器位置與變形量L的關(guān)系。單元理論分析和實(shí)驗(yàn)測(cè)試結(jié)果如圖6(c)所示，數(shù)值優(yōu)化結(jié)果表明，最大變形量L對(duì)應(yīng)的最優(yōu)位置Z1=190.755 mm。由于選取等效彈簧剛度值遵循了Kminx≥Fload，轉(zhuǎn)動(dòng)副之間有摩擦，實(shí)驗(yàn)測(cè)得數(shù)據(jù)整體大于理論分析值，但實(shí)驗(yàn)數(shù)據(jù)與理論計(jì)算趨勢(shì)一致，因此采用等效彈簧的單元模型與柔性蒙皮實(shí)驗(yàn)?zāi)Ｐ突窘咏顑?yōu)位置相差5.7%。

(a)初始狀態(tài) (b)平衡狀態(tài)

將該研究方法擴(kuò)展到整個(gè)后掠角范圍內(nèi)以確定最佳執(zhí)行器位置，通過(guò)之前分析，不同后掠角時(shí)最佳驅(qū)動(dòng)器位置也不同，但是在飛行時(shí)驅(qū)動(dòng)器不能主動(dòng)改變位置，因此需要一個(gè)在全局變后掠范圍內(nèi)將L最大化的位置。優(yōu)化中，后掠角由20°到70°，每7.5°和特殊位置(36.7°)作為節(jié)點(diǎn)。

首先在整個(gè)變后掠范圍內(nèi)給定一個(gè)后掠角，記作@θ，代表此時(shí)為初始后掠角，蒙皮在此時(shí)為初始狀態(tài)，得出對(duì)應(yīng)的驅(qū)動(dòng)器最佳位置Z1(單位長(zhǎng)度，下同)；其次對(duì)每個(gè)最佳位置Z1給定初始后掠角，計(jì)算變形時(shí)的L，如表1所示。通過(guò)式(9)可以計(jì)算出每一個(gè)位置對(duì)應(yīng)的全局變形量Ltotal：

(9)

表1 不同后掠角的最優(yōu)位置在全局變后掠中的變形量Tab.1 Deformation of the optimal position of different sweep angles in global sweep

采用Lagrange插值對(duì)全局變形量Ltotal進(jìn)行擬合，得到全局最優(yōu)驅(qū)動(dòng)器位置Z1=6.39,如圖7所示，與之前所預(yù)測(cè)的最優(yōu)區(qū)間相符合。

圖7 驅(qū)動(dòng)器不同位置與變形量的關(guān)系Fig.7 Relationship between actuator position and deformation

對(duì)全局變后掠驅(qū)動(dòng)器最佳位置優(yōu)化后，得到當(dāng)Z1=6.39時(shí)，在變形范圍內(nèi)輸出位移之和最大，最優(yōu)位置對(duì)應(yīng)的初始后掠角為45°，因此最優(yōu)的柔性剪切變形蒙皮的初始狀態(tài)為碳纖維棒與凱夫拉線(xiàn)的夾角ε=68.2°。這與變形機(jī)翼在飛行過(guò)程中保持中等后掠角條件相一致。

2.2 多單元驅(qū)動(dòng)器布局

圖8 在外載荷和驅(qū)動(dòng)器作用下三單元結(jié)構(gòu)示意圖Fig.8 Schematic diagram of three-cell model with actuators and external load

為了找出結(jié)構(gòu)剛度與分布執(zhí)行器布置之間的關(guān)系，研究以下6種情況：①在單元1中驅(qū)動(dòng)力1 500 N、②在單元3中驅(qū)動(dòng)力1 500 N、③分布在單元1、2中驅(qū)動(dòng)力各1 000 N、④分布在單元1、3中驅(qū)動(dòng)力各1 000 N、⑤分布在單元2、3中驅(qū)動(dòng)力各1 000 N、⑥分布在單元1、2、3中驅(qū)動(dòng)力各500 N。通過(guò)式(5)可知6種狀態(tài)下的理論平衡位置是一致的，其中當(dāng)驅(qū)動(dòng)器數(shù)量有限制時(shí)，③和⑤中的驅(qū)動(dòng)器集中在結(jié)構(gòu)一端，稱(chēng)為弱分布式驅(qū)動(dòng)布局，④分散在兩側(cè)，為強(qiáng)分布式驅(qū)動(dòng)布局。

以結(jié)構(gòu)桿件的截面積參數(shù)d為變量，對(duì)蒙皮初始狀態(tài)在后掠角θ0=36.7°的情況進(jìn)行研究，在外載荷、彈簧剛度、驅(qū)動(dòng)力和設(shè)計(jì)參數(shù)Z1下達(dá)到平衡角，通過(guò)ANSYS進(jìn)行瞬態(tài)動(dòng)力學(xué)仿真，得到6種驅(qū)動(dòng)器布置在一系列截面積參數(shù)情況的總變形，6種布置下驅(qū)動(dòng)器布局與變形量的關(guān)系如圖9所示?？梢钥闯?，當(dāng)結(jié)構(gòu)剛度較低時(shí)，驅(qū)動(dòng)器分布在1、2、3單元中所獲得的變形量均比其余的大，若驅(qū)動(dòng)器數(shù)量減少，位于分散單元1、3的強(qiáng)分布式驅(qū)動(dòng)方式要比位于相對(duì)集中單元1、2或2、3的弱分布式驅(qū)動(dòng)的變形量大；當(dāng)機(jī)翼結(jié)構(gòu)剛度逐漸增強(qiáng)時(shí)，分布式驅(qū)動(dòng)與單個(gè)驅(qū)動(dòng)的變形量均收斂在理想值。

圖9 驅(qū)動(dòng)器布置和剛度對(duì)變形量的影響Fig.9 Influence of actuator arrangement and stiffness on deformation

當(dāng)驅(qū)動(dòng)器靠近翼根，部分驅(qū)動(dòng)器輸入能量作為應(yīng)變能存儲(chǔ)在結(jié)構(gòu)中，結(jié)構(gòu)剛度越低，能量傳遞越無(wú)法傳遞到翼尖，因此，需要分布式驅(qū)動(dòng)以最大化變形量。當(dāng)截面積逐漸增加或采用分布式驅(qū)動(dòng)器，會(huì)通過(guò)減小機(jī)構(gòu)遠(yuǎn)端上的外力引起的撓度來(lái)最大化的傳遞能量。因此分布式驅(qū)動(dòng)系統(tǒng)對(duì)于柔性結(jié)構(gòu)是必不可少的。

3 結(jié)構(gòu)設(shè)計(jì)與樣機(jī)功能演示

變形翼骨架為翼梁與翼肋組成的平行連桿結(jié)構(gòu)，如圖10所示。變形翼根部與導(dǎo)向系統(tǒng)連接，主翼梁與樞軸座鉸接，變形時(shí)，在翼根處絲杠驅(qū)動(dòng)翼肋或內(nèi)置在四邊形單元的分布式驅(qū)動(dòng)器改變單元的夾角，使得前緣前移，可以減緩變后掠帶來(lái)的氣動(dòng)中心后移[18]。該種翼肋布置方式在變后掠過(guò)程中會(huì)使翼肋產(chǎn)生沿氣流方向的擺動(dòng)，因此在變后掠翼稍處安裝順氣流機(jī)架，使得翼尖沿氣流方向，并在機(jī)架上安裝變彎度翼稍，翼梢通過(guò)改變沿翼展方向分布的載荷來(lái)降低巡航時(shí)誘導(dǎo)阻力。得益于“擺動(dòng)式”翼肋，變形翼的機(jī)翼弦長(zhǎng)與相對(duì)厚度也發(fā)生變化，翼型的變化并不是由翼肋本身變形來(lái)實(shí)現(xiàn)的，解決了剛性翼肋變化難的結(jié)構(gòu)問(wèn)題，實(shí)現(xiàn)了機(jī)翼多維度的大變形。

圖10 變形翼骨架結(jié)構(gòu)Fig.10 Morphing wing skeleton structure

根據(jù)設(shè)計(jì)的變形翼進(jìn)行變后掠部分的加工與蒙皮的制作、安裝，為下一步實(shí)驗(yàn)做準(zhǔn)備。變形翼樣機(jī)如圖11所示，在翼根處的集中式驅(qū)動(dòng)下，變形翼可以實(shí)現(xiàn)后掠角20°到70°的連續(xù)變形。后掠角20°時(shí)，機(jī)翼面積S=0.779 m2，最大半翼展1.41 m，最小弦長(zhǎng)0.92 m；后掠角36.7°，Smax=0.890 m2；后掠角70°，Smin=0.540 m2，最小翼展為0.634 m，最大弦長(zhǎng)1.643 m。

(a)θ=20° (b)θ=70°圖11 變形翼樣機(jī)Fig.11 Prototype of morphing wing

4 結(jié) 論

本文設(shè)計(jì)了一種兼顧高低速飛行的多維度變形翼骨架，確定了分布在變形翼骨架結(jié)構(gòu)中執(zhí)行器的最佳位置和蒙皮初始狀態(tài)，對(duì)多個(gè)單元進(jìn)行機(jī)構(gòu)運(yùn)動(dòng)與彈性變形耦合分析，主要結(jié)論如下：

1)翼肋擺動(dòng)的平行連桿變形翼機(jī)構(gòu)在變后掠的同時(shí)，可實(shí)現(xiàn)變面積、變弦長(zhǎng)、變相對(duì)厚度，解決了機(jī)翼平面維度和翼型維度之間的結(jié)構(gòu)矛盾，從結(jié)構(gòu)上減緩變后掠帶來(lái)的氣動(dòng)中心后移。

2)確定當(dāng)順氣流優(yōu)勢(shì)K=0.75時(shí)，后掠角θ在20°～70°范圍內(nèi)，展弦比由4.8減小到1.0，面積變化率為175%，相對(duì)厚度由13.4%減小到7.76%。當(dāng)θ≈36.7°時(shí)，機(jī)翼面積最大，滿(mǎn)足高升力爬升、小阻力機(jī)動(dòng)狀態(tài)。

3)隨著后掠角的增大，驅(qū)動(dòng)器所需最小驅(qū)動(dòng)力的Z1值逐漸減??；建立的單元數(shù)值模型能夠在一定誤差范圍內(nèi)模擬驅(qū)動(dòng)器位置與變形量的關(guān)系，當(dāng)Z1在63.9%處、初始后掠角為45°時(shí)，變形翼在整個(gè)變形范圍內(nèi)驅(qū)動(dòng)效率最高，此時(shí)柔性剪切蒙皮的最佳初始狀態(tài)ε=68.2°。

4)當(dāng)結(jié)構(gòu)剛度較低時(shí)，分布式驅(qū)動(dòng)布局所獲得的變形量均比其余的大。若驅(qū)動(dòng)器數(shù)量減少，驅(qū)動(dòng)器傳遞的一部分能量會(huì)以應(yīng)變能的形式存儲(chǔ)在結(jié)構(gòu)中，強(qiáng)分布式驅(qū)動(dòng)方式要比弱分布式驅(qū)動(dòng)的變形量大。當(dāng)機(jī)翼結(jié)構(gòu)剛度增強(qiáng)時(shí)，分布式驅(qū)動(dòng)與單個(gè)驅(qū)動(dòng)所做的變形量均收斂在理想值，采用多個(gè)驅(qū)動(dòng)器和單個(gè)驅(qū)動(dòng)器對(duì)變形能力影響不大。