吳惠松，林麒，柳汀, 劉震, 師璐，王曉光

廈門大學(xué) 航空航天學(xué)院，廈門 361102

風(fēng)洞虛擬飛行試驗(Wind Tunnel Based Virtual Flight Testing, WTBVFT)是近年發(fā)展起來的一種連接傳統(tǒng)風(fēng)洞試驗和飛行試驗之間紐帶的風(fēng)洞試驗技術(shù)，它采用某種支撐裝置將飛行器模型支撐在風(fēng)洞試驗段中進行吹風(fēng)試驗，約束模型的3個線位移，釋放模型的3個轉(zhuǎn)動角位移，按照模擬相似準(zhǔn)則設(shè)計飛行控制律，操縱舵面來控制模型偏轉(zhuǎn)，檢驗其操縱響應(yīng)特性，研究氣動/運動耦合機理，達到建立氣動/運動/控制一體化評估試驗平臺的目的。同時，支撐裝置還必須滿足轉(zhuǎn)角范圍大、支撐干擾小等要求，以及具備3個轉(zhuǎn)動自由度的自由耦合和解耦運動功能。

為此，眾多學(xué)者對風(fēng)洞虛擬飛行試驗技術(shù)開展研究。美國Gebet、Wehef、Lawrence等采用多根鋼繩通過模型外部軸承環(huán)將導(dǎo)彈模型懸掛支撐在風(fēng)洞中開展虛擬飛行試驗，實現(xiàn)了導(dǎo)彈的俯仰、滾轉(zhuǎn)和偏航運動。英國Lowenberg和Kyle設(shè)計了鐘擺式支撐裝置，用于研究飛行器的俯仰單自由度及俯仰和升沉兩個自由度的被動驅(qū)動運動；Davison、Gatto等設(shè)計了硬式三自由度支撐機構(gòu)，并進行系統(tǒng)建模和穩(wěn)定性導(dǎo)數(shù)試驗；Pattinson等研發(fā)了五自由度支撐機構(gòu)，對機構(gòu)進行動力學(xué)建模，開展俯仰極限振蕩研究等試驗。俄羅斯中央空氣流體動力研究院研發(fā)了三自由度背撐式虛擬飛行試驗支撐機構(gòu)，開展了大迎角失速/偏離、飛行控制律算法驗證等研究。國內(nèi)趙忠良等在2.4 m跨聲速風(fēng)洞中采用吊架式模型支撐系統(tǒng)，實現(xiàn)模型俯仰、滾轉(zhuǎn)自由運動及偏航驅(qū)動控制運動，開展氣動/運動耦合機理研究；李浩開展了風(fēng)洞虛擬飛行相似準(zhǔn)則和模擬方法研究；郭林亮等采用三自由度接觸型球鉸支撐機構(gòu)、兩自由度萬向鉸和旋轉(zhuǎn)曲桿組合三自由度支撐機構(gòu)開展動力學(xué)特性研究；岑飛等開展了模型自由飛試驗和虛擬飛行試驗相關(guān)技術(shù)研究；南京航空航天大學(xué)、中航工業(yè)空氣動力研究院、中國航天空氣動力技術(shù)研究院等單位也開展了相關(guān)研究。

本文充分利用繩系機器人用于風(fēng)洞試驗的模型支撐機構(gòu)簡單，且繩索鉸接點具備球副的功能，以及繩索對氣流干擾小、可智能控制等優(yōu)勢，基于螺旋理論的機構(gòu)構(gòu)型設(shè)計方法和自由度計算原理，設(shè)計了風(fēng)洞虛擬飛行試驗?zāi)Ｐ屠K系支撐機構(gòu)(Wire Driven Parallel Robot suspension for Virtual Flight Testing model,WDPR-VFT)，進行穩(wěn)定性分析和結(jié)構(gòu)優(yōu)化設(shè)計；再對其機構(gòu)進行運動學(xué)和動力學(xué)建模，并進行相應(yīng)的仿真分析和試驗驗證。

1 構(gòu)型設(shè)計與分析

1.1 構(gòu)型設(shè)計

本文在空間平行四桿機構(gòu)和平面三平行五桿機構(gòu)的基礎(chǔ)上，結(jié)合風(fēng)洞虛擬飛行試驗?zāi)Ｐ椭蔚男枨?，設(shè)計了兩個空間五桿平行四邊形機構(gòu)經(jīng)正交疊加組合，并分別與上下兩個偏航轉(zhuǎn)動副構(gòu)成虛擬飛行試驗?zāi)Ｐ椭螜C構(gòu)，如圖1(a)所示。圖1(b)給出的支撐平臺可實現(xiàn)釋放或鎖定模型的轉(zhuǎn)動自由度，及調(diào)整模型不同的轉(zhuǎn)動幅度的功能。即該機構(gòu)可在風(fēng)洞里支撐著飛機模型的同時通過打舵實現(xiàn)飛機模型做俯仰、偏航和滾轉(zhuǎn)運動，也就是做虛擬飛行試驗。

圖1 虛擬飛行試驗?zāi)Ｐ椭螜C構(gòu)Fig.1 Suspension mechanism of virtual flight test model

為便于計算自由度數(shù)目，先將圖1中繩索視為連桿進行分析，且用于張緊繩索的移動副不參與自由度計算。該機構(gòu)可看成由上下2個串并聯(lián)分支(RR[RRPP])組成的并聯(lián)機構(gòu)，每個分支中的空間五桿平行四邊形機構(gòu)又可看成由與Delta機構(gòu)子鏈類似的4S機構(gòu)和四邊形平面垂直的轉(zhuǎn)動副構(gòu)成的串并聯(lián)分支(R[RRPP])，其中R、P、S分別表示轉(zhuǎn)動副、移動副和球副，如圖2(a)所示，紅色十字架為飛機的抽象模型(即動平臺)。機構(gòu)的等效結(jié)構(gòu)簡圖如圖2(b)所示，不失一般性，建立圖中所示的坐標(biāo)系，以便于自由度分析。

(1)

式中：4個運動螺旋分別表示單獨的平行四邊形機構(gòu)的桿可繞、兩個軸線轉(zhuǎn)動，及沿、兩軸線方向的移動。其與Delta機構(gòu)4S閉環(huán)子鏈類似，可將該四桿機構(gòu)等效為2轉(zhuǎn)動2平移的[RRPP]機構(gòu)。

表1 兩個平行四邊形機構(gòu)分析Table 1 Analysis of two parallelogram mechanisms

式(1)與點的運動螺旋、偏航轉(zhuǎn)動螺旋串聯(lián)組成虛擬飛行試驗?zāi)Ｐ椭螜C構(gòu)第1分支(俯仰和偏航)的運動螺旋系：

圖2 WDPR-VFT機構(gòu)簡化分析Fig.2 Simplified analysis of WDPR-VFT mechanism

(2)

運動螺旋和均線性相關(guān)，說明第1分支存在局部自由度。

對式(2)求約束反螺旋，得

(3)

式(3)表明第1個分支對十字架模型產(chǎn)生了一個約束力，約束其沿軸線方向的平移。

(4)

式(4)表明第2個分支對十字架模型也產(chǎn)生了一個約束力，約束其沿軸線方向的平移。

式(3)和式(4)共同構(gòu)成虛擬飛行試驗?zāi)Ｐ椭螜C構(gòu)的反螺旋系，它們共同約束了十字架模型沿軸線方向的平移，即十字架模型此時有5個自由度?？梢钥闯?，兩個反螺旋線性相關(guān)，存在一個過約束。

顯然，若對于剛性硬式桿件構(gòu)成的支撐，十字架模型無法繞、軸線方向旋轉(zhuǎn)，也無法沿、、方向移動。十字架模型只剩下繞軸線方向旋轉(zhuǎn)，以及繞模型質(zhì)心所在相應(yīng)平行四邊形機構(gòu)平面上的轉(zhuǎn)動中心轉(zhuǎn)動的3個自由度。

1.2 自由度數(shù)目驗算

為驗證上述機構(gòu)的自由度數(shù)目，用修正的G-K自由度公式進行計算。

十字架模型的自由度為

6(10-12-1)+28+1-(4+2+2)=3

(5)

式中：、、、和分別表示包含機架的構(gòu)件數(shù)、運動副數(shù)、第個運動副的自由度數(shù)、過約束數(shù)和局部自由度數(shù)。

這里有8個局部自由度，除、、、桿分別存在一個局部自由度外，還有偏航轉(zhuǎn)動副運動螺旋分別與對應(yīng)的分支組合后產(chǎn)生的局部自由度，以及十字架模型繞自身桿件軸線的旋轉(zhuǎn)相當(dāng)于另一個平行四邊形機構(gòu)的轉(zhuǎn)動，其存在2個局部自由度。

因此，若采用硬式桿件，該虛擬飛行試驗?zāi)Ｐ椭螜C構(gòu)自由度總數(shù)為3，即十字架模型具有俯仰、滾轉(zhuǎn)、偏航自由度，說明該機構(gòu)滿足虛擬飛行試驗?zāi)Ｐ偷闹涡枨蟆?/p>

2 結(jié)構(gòu)優(yōu)化設(shè)計與穩(wěn)定性

2.1 結(jié)構(gòu)優(yōu)化設(shè)計

由于真實飛機模型的氣動外形較為復(fù)雜，若采用繩索直接連接在機身外殼上，則圖1中的支撐桿、繩索和模型結(jié)合起來很難構(gòu)成完整的平行四邊形機構(gòu)。這會導(dǎo)致模型在轉(zhuǎn)動過程中，轉(zhuǎn)動前后模型質(zhì)心與機構(gòu)轉(zhuǎn)動中心不重合，產(chǎn)生不穩(wěn)定現(xiàn)象。

為解決這個問題，在飛機模型內(nèi)部設(shè)計了如圖3所示的連接機構(gòu)，主要是在縱橫軸兩端分別安裝了不同長度的軸承套搖桿用于繩索連接。

圖3 飛機模型內(nèi)部連接機構(gòu)Fig.3 Internal connecting mechanism of model

同時將繞軸承中心轉(zhuǎn)動的支撐桿優(yōu)化設(shè)計為圖4所示的形式，即支撐桿可沿偏航轉(zhuǎn)臺方向移動，移動距離等于模型質(zhì)心到原轉(zhuǎn)動中心的垂直距離。

圖4 優(yōu)化后的支撐平臺Fig.4 Optimized support platform

本文采用如圖5所示的一款動態(tài)試驗標(biāo)模?？紤]到風(fēng)洞大小以及阻塞比要求，在國內(nèi)動態(tài)標(biāo)準(zhǔn)模型的基礎(chǔ)上進行了0.591倍縮尺比設(shè)計加工。模型幾何尺寸見圖5，總質(zhì)量2.476 kg。在滿足幾何相似外，動力學(xué)相似準(zhǔn)則主要是雷諾數(shù)。由于本文的試驗來流速度約10 m/s，特征長度取為風(fēng)洞出口截面邊長450 mm，則約為3.23×10。

圖5 動態(tài)試驗標(biāo)模縮比模型Fig.5 Scaled model of dynamic test standard model

用該動態(tài)模型以圖4的支撐平臺搭建了WDPR-VFT模型原理樣機。其機構(gòu)示意圖如圖6所示。圖中所示的3個坐標(biāo)系，分別是地面靜坐標(biāo)系、機體坐標(biāo)系(機身軸線沿方向，豎直向下指向地心)、桿件坐標(biāo)系、(下標(biāo)u代表上支撐桿，d代表下支撐桿)。

圖6 WDPR-VFT模型原理樣機示意圖Fig.6 Principle prototype of WDPR-VFT

2.2 穩(wěn)定性

為驗證上述設(shè)計的WDPR-VFT機構(gòu)的穩(wěn)定性，采用Adams軟件進行運動仿真分析，以考察模型穩(wěn)定性。假設(shè)圖5的動態(tài)標(biāo)模受到舵偏作用，即相對于其抽象模型十字架的質(zhì)心作用一個100 N·m 的俯仰力矩，兩根上拉和兩根下拉繩索初始預(yù)緊力分別為32.13 N和20 N。下面取其中一個平行四邊形機構(gòu)(以俯仰為例)進行分析。為分析方便，將標(biāo)模對稱面簡化為直線和矩形，并利用CAD作圖分析標(biāo)模轉(zhuǎn)動前后質(zhì)心和轉(zhuǎn)動中心是否重合。

當(dāng)模型縱向豎直對稱面簡化為直線時，即繩索連接在縱橫軸上，支撐桿、兩根繩索和模型一同構(gòu)成平行四邊形機構(gòu)，如圖7(a)所示。當(dāng)模型在俯仰力矩(例如，給模型升降舵一個舵偏角)作用下繞軸承中心轉(zhuǎn)過角時，平行四邊形機構(gòu)轉(zhuǎn)動前后的模型質(zhì)心、支撐機構(gòu)的轉(zhuǎn)動中心均重合，不會出現(xiàn)振蕩現(xiàn)象，由于無摩擦耗散，俯仰角呈線性遞增趨勢，如圖7(b)所示。

圖7 模型對稱面簡化為直線Fig.7 Model symmetry plane simplified as straight line

當(dāng)模型對稱面簡化為矩形時，即繩索連接在軸承套搖桿上，支撐桿、兩根繩索、模型和兩個軸承套搖桿不能構(gòu)成完整的平行四邊形機構(gòu)，如圖8(a) 所示。這時若模型在俯仰力矩作用下繞軸承中心轉(zhuǎn)過角，平行四邊形機構(gòu)轉(zhuǎn)動前后的模型質(zhì)心、支撐機構(gòu)的轉(zhuǎn)動中心會不重合，轉(zhuǎn)動過程中右邊的繩索傾斜，存在使模型恢復(fù)到平衡狀態(tài)的反力矩，俯仰角出現(xiàn)來回振蕩現(xiàn)象，如圖8(b)所示。

采用2.1節(jié)的結(jié)構(gòu)優(yōu)化，即將繞軸承中心轉(zhuǎn)動的支撐桿沿偏航轉(zhuǎn)臺方向移動，移動距離等于模型質(zhì)心到原轉(zhuǎn)動中心的垂直距離，繩索長度相應(yīng)加長，如圖9(a)所示。這時機構(gòu)能保持穩(wěn)定狀態(tài)，模型姿態(tài)角變化如圖9(b)所示，不會出現(xiàn)來回振蕩現(xiàn)象。

圖9 模型對稱面簡化為矩形(優(yōu)化后)Fig.9 Model symmetry plane simplified as rectangle (after optimization)

3 運動學(xué)和動力學(xué)建模

本文所設(shè)計的飛行器模型繩系支撐機構(gòu)屬于欠約束機構(gòu)。在不改變初始繩長的條件下，模型在繩索支撐下，受來流和舵面力矩作用，實現(xiàn)模型在3個轉(zhuǎn)軸上的自由運動。飛行器模型做受迫運動時(例如舵偏作用時)，位姿不僅受繩索運動狀態(tài)的影響，還與該機構(gòu)的動力學(xué)特性有關(guān)。

為驗證上述飛行器模型支撐機構(gòu)的可行性，并預(yù)測飛行器模型在虛擬飛行試驗過程中的運動，對該機構(gòu)進行飛行器模型的運動學(xué)和動力學(xué)建模。

3.1 運動學(xué)建模

首先建立牽引繩長度與飛行器模型位姿之間的運動學(xué)關(guān)系。

如圖10所示，在地面坐標(biāo)系下，第根繩的繩長矢量為

=-=+--

(6)

圖10 WDPR-VFT的運動學(xué)建模Fig.10 Kinematic modeling of WDPR-VFT

(7)

3.2 動力學(xué)建模

本文在飛行器模型無約束六自由度運動方程基礎(chǔ)上添加由4根繩索施加在模型機體上的約束力和約束力矩，同時考慮繩索的彈性阻尼和支撐桿轉(zhuǎn)動軸承的摩擦力矩，建立支撐桿繞桿件質(zhì)心轉(zhuǎn)動的運動方程和有繩支撐約束下的飛行器模型動力學(xué)方程。

3.2.1 飛行器模型無約束六自由度運動方程

飛行器模型屬于六自由度運動體，基于動量定理、動量矩定理和哥氏定理可推導(dǎo)出模型無約束六自由度運動方程，鑒于篇幅原因，不在此詳細列出。

3.2.2 支撐桿繞其質(zhì)心轉(zhuǎn)動的運動學(xué)/動力學(xué)方程

建立如圖11所示的坐標(biāo)系，為建模方便，將機體重心初始位置設(shè)置在地面坐標(biāo)系下的[0,0,-0.775] m，上桿坐標(biāo)系的軸和機體坐標(biāo)系的軸相互垂直，下桿坐標(biāo)系的軸和機體坐標(biāo)系的軸同向。

1) 繩索張力和軸承摩擦力矩建模

由于繩索受拉會產(chǎn)生彈性變形，十字架模型在運動過程中可能出現(xiàn)微小振蕩，對機構(gòu)的穩(wěn)定性產(chǎn)生一定的影響，因此，動力學(xué)建?？紤]了繩索的彈性阻尼，其張力模型為

(8)

式中：為繩拉力矢量；為繩索剛度；為繩阻尼系數(shù)；為實時繩長矢量；為未變形的繩長矢量。本文采用直徑1.2 mm的Kevlar繩作為牽引繩，其彈性模量=43.9 GPa,=/=1.323× 10N/m，為繩索橫截面積。

上支撐桿轉(zhuǎn)動軸承的摩擦力矩模型為

=sgn()(+(-)e-||)+

(9)

式中：為庫侖摩擦力矩；為靜摩擦力矩；參數(shù)的引入使模型符合Stribeck效應(yīng)；為黏性系數(shù)。

2) 繩索對模型的合力/力矩

4根繩索對機體產(chǎn)生的合力、合力矩為

(10)

式中：g是第根繩索在地面坐標(biāo)系下的繩拉力；是地面坐標(biāo)系到機體坐標(biāo)系的變換矩陣；b為第根繩索對應(yīng)飛行器模型上的點在體軸系下的位置向量，=1,2,3,4分別表示圖11中上左繩、上右繩、下前繩、下后繩。

圖11 WDPR-VFT的動力學(xué)建模坐標(biāo)系定義Fig.11 Dynamic modeling coordinate system definition of WDPR-VFT

3) 支撐桿所受合力矩

單根繩索對支撐桿在桿系坐標(biāo)下的作用力矩為

u(d)=u(d)×u(d)

(11)

式中：u(d)是第根繩索對應(yīng)上桿(下桿)上的點在桿系下的位置向量；u(d)是第根繩索在上桿(下桿)坐標(biāo)系下的繩拉力。

上、下支撐桿所受合力矩和為

(12)

4) 支撐桿繞其質(zhì)心轉(zhuǎn)動的運動學(xué)/動力學(xué)方程

上(下)支撐桿繞自身質(zhì)心轉(zhuǎn)動的運動學(xué)和動力學(xué)微分方程為

(13)

(14)

式中：、、表示上支撐桿的三軸姿態(tài)角；u、u、u表示上支撐桿的角速率沿3個坐標(biāo)軸的分量；u、u、u表示上支撐桿所受合力矩在3個坐標(biāo)軸上的分量；為歐拉轉(zhuǎn)換矩陣；為上支撐桿的轉(zhuǎn)動慣量矩陣。

3.2.3 繩索支撐約束下飛行器模型的運動方程

本文開展虛擬飛行試驗?zāi)Ｐ屠K系支撐機構(gòu)驗證試驗，采用圖5所示的飛行器模型，其氣動模型(氣動力和氣動力矩)從風(fēng)洞試驗氣動數(shù)據(jù)庫和相關(guān)資料中收集并通過插值計算得到。

飛行器模型所受合力和合力矩為

(15)

將飛行器模型所受合力和合力矩代入飛行器模型的無約束六自由度運動方程，可得繩索支撐約束下飛行器模型的運動學(xué)和動力學(xué)方程。

式(15)建立了有繩支撐約束下飛行器模型狀態(tài)量與舵面控制量之間的函數(shù)關(guān)系，以及上、下支撐桿的轉(zhuǎn)動角度、角速度與繩索拉力、軸承摩擦力矩之間的關(guān)系。

飛行器模型的位姿與繩長、繩拉力存在內(nèi)在聯(lián)系，即對于本文設(shè)計的欠約束、可重構(gòu)(牽引點時變)、被動驅(qū)動的繩系支撐機構(gòu)，幾何運動方程和力平衡方程相互耦合，可通過聯(lián)合求解確定其相關(guān)參數(shù)。

4 仿真分析

為驗證上述飛行動力學(xué)建模的正確性和支撐機構(gòu)的可行性，利用MATLAB中的Simulink模塊開展仿真分析。選取3種不同工況進行對比研究：模型無約束3自由度轉(zhuǎn)動運動(3DOF Unconstrained)；模型在有繩支撐約束下的無摩擦力矩的運動(WDPRVFT)；模型在有繩支撐約束下的有摩擦力矩運動(WDPRVFT+ Friction)。

4.1 數(shù)學(xué)模型穩(wěn)定性

首先，初始化配平(配平迎角為2.1°)，在無舵偏角和風(fēng)速情況下，給飛行器模型1°的抬頭俯仰角，繩索初始預(yù)緊力為50 N，仿真考察模型和俯仰支撐桿的變化情況。

如圖12所示，初始時刻飛行器模型在受到負的俯仰力矩作用下俯仰角減小，模型低頭向下，俯仰支撐桿(下桿)受正的俯仰力矩作用使其俯仰角增大，兩者在平衡位置來回振蕩，振幅不斷衰減，模型俯仰角最終穩(wěn)定在0.79°，說明曲線趨勢是合理的，機構(gòu)是穩(wěn)定的，其中俯仰支撐桿振蕩幅值較大，與繩索預(yù)緊力和轉(zhuǎn)動慣量有關(guān)。

上述結(jié)果表明，由于繩索的彈性阻尼，機構(gòu)運動過程中會產(chǎn)生微小的振蕩。但是很快就穩(wěn)定下來。若通過機構(gòu)絲桿組件中的移動副拉緊繩索，可以增大繩索的預(yù)緊力以減弱這種微小的振蕩。

圖12 WDPR-VFT機構(gòu)振蕩仿真Fig.12 WDPR-VFT mechanism oscillation simulation

4.2 典型縱向操縱響應(yīng)

本文就WDPR-VFT的典型縱向操縱響應(yīng)進行研究。選取典型縱向操縱輸入信號，即在來流風(fēng)速為25 m/s條件下，待風(fēng)速穩(wěn)定后，輸入1°的升降舵階躍，階躍點在10 s處，升降舵的操縱響應(yīng)結(jié)果如圖13和圖14所示。仿真結(jié)果說明飛行器模型對1°的升降舵階躍會產(chǎn)生負的純俯仰力矩，模型的俯仰角及其變化率在幾個來回振蕩后進入穩(wěn)定狀態(tài)。

圖13 俯仰角變化Fig.13 Change of pitch angle

圖14 俯仰角速率變化Fig.14 Change of pitch angle rate

摩擦力矩會對模型的俯仰運動產(chǎn)生阻尼作用，阻尼使俯仰角變化速率峰值變小，在階躍后的第1個波峰處減小了0.12(°)/s，第2個波峰處減小了0.2(°)/s，并產(chǎn)生了死區(qū)現(xiàn)象，即在第2個波峰過零時(約在階躍后的1.92 s)不再振蕩，而是直接進入穩(wěn)定狀態(tài)，比無摩擦力矩的結(jié)果提前約1.4 s。

與俯仰角變化速率對應(yīng)，阻尼使俯仰角峰值在階躍后的0.96 s時減小0.03°，并在1.9 s后進入穩(wěn)定狀態(tài)，與無摩擦的結(jié)果相比減小了0.04°。說明摩擦力矩對支撐機構(gòu)的功能影響很小。

上述仿真結(jié)果說明本文設(shè)計的虛擬飛行試驗?zāi)Ｐ屠K系支撐機構(gòu)與模型無約束3自由度轉(zhuǎn)動仿真結(jié)果趨勢是一致的。由于摩擦力矩和繩索彈性的影響，WDPR-VFT模型并聯(lián)機構(gòu)與無約束3自由度轉(zhuǎn)動相比，滯后約0.26 s。

在上述縱向操縱輸入下，模型的質(zhì)心位移變化如圖15所示。圖15說明飛行器模型質(zhì)心沿三軸的位移變化很小，沿機頭負方向的位移最大為0.016 mm，振蕩幅值越來越小并最終進入穩(wěn)定狀態(tài)，說明該虛擬飛行試驗?zāi)Ｐ屠K系支撐機構(gòu)基本約束了飛行器模型的3個平動自由度，釋放了3個轉(zhuǎn)動自由度。

圖15 模型質(zhì)心位移變化Fig.15 Displacement of model center of mass

上述分析結(jié)果表明，對于縱向操縱，WDPR-VFT模型支撐機構(gòu)對飛行器模型運動影響很小，具有良好的自穩(wěn)定性，驗證了動力學(xué)建模仿真結(jié)果的準(zhǔn)確性，可以滿足風(fēng)洞虛擬飛行試驗的需求。

5 試驗驗證

為驗證所設(shè)計的虛擬飛行試驗?zāi)Ｐ屠K系支撐機構(gòu)的可行性，在風(fēng)洞里開展了俯仰單自由度試驗，即進行升降舵開環(huán)控制試驗，研究模型運動響應(yīng)過程。

5.1 試驗設(shè)備

本文根據(jù)動力學(xué)相似準(zhǔn)則，設(shè)計并加工了如圖16所示的帶舵面的試驗?zāi)Ｐ?具體尺寸參見圖5)，配置了開源飛行控制器Pixhawk、金屬齒數(shù)碼舵機等零件，用于模型姿態(tài)角測量和舵面控制，并將其安裝于實驗室的低速直流式風(fēng)洞出口處，見圖17。

圖16 試驗?zāi)Ｐ图皟?nèi)部配置Fig.16 Experimental model and internal parts configuration

圖17 虛擬飛行試驗繩系支撐機構(gòu)樣機Fig.17 Prototype of WDPR-VFT

該風(fēng)洞試驗段橫截面為300 mm×300 mm，但對于圖16的模型來說顯得太小，故將模型放在風(fēng)洞擴散段出口處，此處橫截面積為450 mm×450 mm，風(fēng)速可達約30 m/s。但即便如此，由于風(fēng)洞試驗條件限制，擴散段出口尺寸仍小于模型翼展。而且風(fēng)洞擴散段后的流場也不太好，不是均勻流場，對試驗結(jié)果必然造成影響。

5.2 試驗結(jié)果

5.2.1 零舵偏時模型的縱向動態(tài)特性

驗證試驗在風(fēng)洞來流速度10 m/s下進行。當(dāng)對模型釋放俯仰單自由度后，在無舵偏的情況下，模型呈負的俯仰角，并出現(xiàn)小幅的俯仰阻尼振蕩，俯仰角繼續(xù)減小，最終穩(wěn)定在平衡迎角處，如圖18所示。

圖18 俯仰單自由度釋放后模型運動響應(yīng)結(jié)果Fig.18 Model motion response after release of single pitch degree of freedom

課題組曾在在空氣動力研究院FL-5低速風(fēng)洞以圖5的模型采用WDPR繩系機器人支撐和彎刀尾支撐進行吹風(fēng)，得到過升降舵偏角為0°時的俯仰力矩系數(shù)對比試驗結(jié)果。另外，國內(nèi)某氣動單位也對相同模型進行過吹風(fēng)試驗。由以上風(fēng)洞試驗結(jié)果可知，當(dāng)=0°時，在=0°處，俯仰力矩系數(shù)<0，存在微小的低頭力矩。

因此，圖18中模型的俯仰角變化是合理的，符合該模型的氣動特性，也表明WDPR-VFT的縱向動態(tài)特性良好。雖然上述試驗是在風(fēng)洞來流條件欠佳的情況下獲得，但從一個側(cè)面反映本文提出的WDPR-VFT支撐機構(gòu)滿足風(fēng)洞虛擬飛行試驗的基本要求。

5.2.2 俯仰短周期模態(tài)運動

圖19所示為用WDPR-VFT支撐模型吹風(fēng)時的升降舵開環(huán)控制試驗曲線。模型在配平后進行升降舵階躍方波機動，激勵模型做俯仰短周期模態(tài)運動。試驗中，模型配平初始俯仰角為-4.5°，在第一個階躍方波機動中，升降舵從0°先往舵面負方向(即圖中Δ的正值方向)偏轉(zhuǎn)至10°，模型獲得抬頭的俯仰力矩，俯仰角升至22°左右并伴有小幅振蕩；升降舵再往舵面正方向(向下偏轉(zhuǎn)為正)偏轉(zhuǎn)至14°，模型獲得低頭的俯仰力矩，俯仰角降至-12°左右并伴有小幅振蕩后趨于穩(wěn)定；最后再將升降舵回零，模型獲得抬頭的俯仰力矩，俯仰角回至0°附近。

圖19 升降舵開環(huán)控制試驗曲線Fig.19 Open-loop control test curves in elevator maneuvering flight model test

階躍方波機動重復(fù)了3次，模型運動響應(yīng)狀態(tài)基本一致。結(jié)果表明，圖5所示的動態(tài)試驗標(biāo)模的升力特性明顯，升降舵往舵面負方向偏轉(zhuǎn)獲得的模型俯仰角變化比升降舵往舵面正方向偏轉(zhuǎn)獲得的模型俯仰角變化大得多。每次方波階躍機動后模型俯仰角都有一個小幅的振蕩并趨于穩(wěn)定。這是由于模型的慣性與空氣阻尼相互作用的結(jié)果，也與典型縱向操縱仿真結(jié)果(見圖13)總體趨勢類似。

6 結(jié) 論

1) 運用螺旋理論對WDPR-VFT進行了構(gòu)型設(shè)計和自由度分析，證明該機構(gòu)可以實現(xiàn)飛行器模型3個轉(zhuǎn)動自由度的自由耦合和解耦運動，滿足風(fēng)洞虛擬飛行試驗?zāi)Ｐ椭蔚男枨蟆?/p>

2) 對WDPR-VFT進行結(jié)構(gòu)優(yōu)化和穩(wěn)定性分析，證明優(yōu)化后的機構(gòu)具有良好的穩(wěn)定性。

3) 對WDPR-VFT建立了繩系欠約束、可重構(gòu)、被動驅(qū)動機構(gòu)的運動學(xué)和動力學(xué)的數(shù)學(xué)模型，并采用仿真方法分析在虛擬飛行試驗中的模型對典型縱向操縱的響應(yīng)，結(jié)果表明WDPR-VFT模型支撐機構(gòu)對飛行器模型運動影響很小，驗證了所建數(shù)學(xué)模型的可行性和有效性。

4) 在低速直流式風(fēng)洞對WDPR-VFT機構(gòu)的可行性進行了初步的試驗驗證，結(jié)果表明機構(gòu)穩(wěn)定性好，與典型縱向操縱仿真結(jié)果趨勢一致。

本文設(shè)計的用于風(fēng)洞虛擬飛行試驗?zāi)Ｐ偷睦K系支撐欠約束機構(gòu)WDPR-VFT，可為開展飛行器模型氣動/運動/控制一體化研究，探索飛行器模型氣動/運動耦合機理、開展模型氣動參數(shù)辨識提供條件，為風(fēng)洞虛擬飛行試驗?zāi)Ｐ驮O(shè)計繩系支撐機構(gòu)提供參考。