平板翼型地面效應(yīng)數(shù)值模擬與分析

陳世適, 熊芬芬

(北京理工大學(xué) 飛行器動(dòng)力學(xué)與控制教育部重點(diǎn)實(shí)驗(yàn)室, 北京 100081)

平板翼型地面效應(yīng)數(shù)值模擬與分析

陳世適, 熊芬芬

(北京理工大學(xué) 飛行器動(dòng)力學(xué)與控制教育部重點(diǎn)實(shí)驗(yàn)室, 北京 100081)

目前大多數(shù)有關(guān)地面效應(yīng)的研究僅限于平直地面,并非現(xiàn)實(shí)中更為普遍的彎曲地面,從而極大地影響了研究結(jié)果的有效性。對靠近地面的二維平板翼型升力變化情況進(jìn)行了研究,包括平直地面、山谷和山坡彎曲地面,采用復(fù)變函數(shù)理論計(jì)算流場中的復(fù)速度,同時(shí)運(yùn)用保角變換將平板繞流變換為圓柱繞流,并且采用離散渦法建立了地面數(shù)學(xué)模型,將其表示為一系列離散邊界元。研究結(jié)果表明,當(dāng)平板翼型靠近平直地面時(shí),其升力明顯增大;當(dāng)其飛過山谷或山坡時(shí),其升力波動(dòng)相當(dāng)激烈,故在設(shè)計(jì)地效飛行器時(shí)應(yīng)妥當(dāng)考慮和處理此現(xiàn)象。

地面效應(yīng); 彎曲地面; 復(fù)變函數(shù); 保角變換; 離散渦

0 引言

當(dāng)飛行器從地面起飛或者降落至地面時(shí),由于此時(shí)飛行器非常接近地面,其機(jī)翼周圍氣流的流場有所變化,作用在飛行器上的升力、阻力及俯仰力矩也發(fā)生了變化,該現(xiàn)象稱之為地面效應(yīng)[1-2]。目前,有關(guān)地面效應(yīng)的研究甚多。Wieselsberger[3]首次使用鏡像法對地效飛行器進(jìn)行了研究,Tomotika等人[4]基于單翼翼型的結(jié)果對靠近平直地面的平板翼型升力進(jìn)行了理論計(jì)算研究,Ivanteeva等人[5]研究了靠近地面的特定翼型的氣動(dòng)特性,Terentiev等人[6]使用數(shù)值方法對靠近地面翼型的流場進(jìn)行了仿真研究,張亮等人[7]對二維地效翼氣動(dòng)力進(jìn)行了計(jì)算和性能研究,李盾和王義寧[8]對地效飛行器三維地面效應(yīng)進(jìn)行了數(shù)值模擬。最近,王中清等人[9]對飛行器地面效應(yīng)定常/非定常流場進(jìn)行了數(shù)值模擬。

然而,以上研究多數(shù)僅限于平直地面,鮮有對彎曲地面進(jìn)行研究?，F(xiàn)實(shí)中機(jī)場周圍的地形并非都是平坦的,時(shí)常伴有山谷或者山坡,機(jī)場跑道也偶爾有凸起或者凹坑。當(dāng)飛行器起飛或者降落時(shí),機(jī)翼與跑道地面之間的距離將隨著地面形狀而改變,導(dǎo)致飛行器升力波動(dòng),從而影響飛行器的安全系數(shù)。

本文將對靠近平直和彎曲地面的平板翼型升力變化情況進(jìn)行研究。首先將平板翼型周邊流場描述為平面勢流,通過復(fù)變函數(shù)理論可求得流場的速度勢與流函數(shù);然后運(yùn)用保角變換將平板繞流變換為圓柱繞流,同時(shí)采用離散渦法處理地面對流場的影響,在地表面散布一系列離散點(diǎn)渦,將連續(xù)的地面問題轉(zhuǎn)換為離散的邊界元問題。該方法比其它數(shù)值方法更能節(jié)省計(jì)算時(shí)間,并可運(yùn)用于各種形狀的復(fù)雜彎曲地面問題。文中首先對平直地面情況進(jìn)行研究以驗(yàn)證上述方法的有效性,然后將上述方法應(yīng)用于彎曲地面情況,可求得翼型高度以及地面彎曲半徑對平板翼型升力的影響。

1 數(shù)學(xué)模型及解決方法

圖1 山谷彎曲地面簡化原理圖Fig.1 Simplified schematic diagram of hollow (<0)

不可壓流體的任何一種平面勢流必定具有一個(gè)確定的復(fù)勢W(z)=φ(x,y)+iψ(x,y),其中,φ(x,y)和ψ(x,y)分別為速度勢和流函數(shù)[10]。研究中采用儒可夫斯基變換z=(ζ+c2/ζ)/2(c為實(shí)常數(shù)),其能夠?qū)⒃趜平面中平板繞流變換為在相應(yīng)ζ平面中圓柱繞流[10-11],有利于將已知的圓柱繞流研究結(jié)果運(yùn)用于平板翼型的研究。

假定沿著x和y方向的速度分別為u和v。根據(jù)復(fù)變函數(shù)理論,地面中每個(gè)控制點(diǎn)在直勻流和點(diǎn)渦流中的復(fù)速度為[11]:

(1)

根據(jù)離散渦法,地面可分為n個(gè)邊界單元,同時(shí)每個(gè)邊界單元都含有一個(gè)控制點(diǎn)和點(diǎn)渦。其中,控制點(diǎn)位于邊界單元長度的1/4處,而點(diǎn)渦位于3/4處,jth個(gè)點(diǎn)渦的強(qiáng)度為Γj(見圖1)。受到地面的約束,地面中每個(gè)控制點(diǎn)的法向速度都為零,且其計(jì)算可分為兩部分:受到地面上方直勻流的影響;受到地面中分布離散點(diǎn)渦的影響。

由圓定理可知,在ζ二維平面直勻流和點(diǎn)渦流中,圓柱體無環(huán)量流動(dòng)的復(fù)勢函數(shù)可分別表示為[11]:

(2)

(3)

式中,V∞為無邊界直勻流的復(fù)速度;Rc為平板翼型保角變換至ζ平面的圓半徑;Γ為點(diǎn)渦的強(qiáng)度;ζv為點(diǎn)渦在ζ平面的位置。

式(4)表示ith控制點(diǎn)在jth點(diǎn)渦流中復(fù)速度的簡化形式,其中Ai,j定義為jth點(diǎn)渦對ith控制點(diǎn)復(fù)速度的影響因子。

(4)

為了滿足庫塔-儒可夫斯基假設(shè)(K-J條件),平板翼型的后緣應(yīng)當(dāng)為一個(gè)額外的控制點(diǎn)(n+1),其沿翼面的法向速度由于受到翼型表面的影響也為零,同時(shí)需要在ζ平面的原點(diǎn)處添加一個(gè)額外環(huán)繞翼型的環(huán)量Γn+1。匯總直勻流和所有點(diǎn)渦流,ith控制點(diǎn)在流場中的總復(fù)速度為:

=ui-ivi(i,j=1,2,…,n,n+1)

(5)

由于受到地面和K-J條件的約束,每個(gè)控制點(diǎn)的法向速度都為零。

對于平直地面情況,根據(jù)幾何關(guān)系(見圖2),該區(qū)域控制點(diǎn)在流場中的法向速度可表示為:

Vi,n= -uisinα+vicosα=0

(i=1,2,…,n+1)

(6)

圖2 平直地面上控制點(diǎn)的法向速度Fig.2 Normal speed of control point on flat ground

對于山谷與山坡彎曲地面情況,根據(jù)它們不同的幾何關(guān)系(見圖3和圖4),對應(yīng)區(qū)域控制點(diǎn)的法向速度為:

Vi,n=uicos(θi+α)+visin(θi+α)=0

(θi=0～π 且i∈山谷地面)

(7)

Vi,n= -uicos(θi-α)+visin(θi-α)=0

(θi=0～π 且i∈山坡地面)

(8)

圖3 山谷部分控制點(diǎn)的法向速度Fig.3 Normal speed of control point in hollow

圖4 山坡部分控制點(diǎn)的法向速度Fig.4 Normal speed of control point on hill

故可通過解式(6)～式(8)獲得每個(gè)點(diǎn)渦的強(qiáng)度Γj。根據(jù)流場中的疊加性,在靠近地面和滿足K-J條件的情況下,平板翼型周圍流場的復(fù)勢函數(shù)可表示為:

圖5 平板翼型靠近平直地面、山谷與山坡 彎曲地面時(shí)的流線圖 Fig.5 Streamlines around plate airfoil near flat ground, hollow and hill

環(huán)繞平板翼型的全部渦環(huán)量可以通過閉環(huán)積分計(jì)算可得。

根據(jù)K-J理論,單位翼展長度的升力系數(shù)可表示為:

式中,ρ和V∞分別為流體的密度和直勻流的速度。

為了能清晰地顯示地面對翼型升力系數(shù)大小的影響,可額外定義一個(gè)新的無量綱參數(shù),即翼型在有邊界流和無邊界流中升力系數(shù)的比值:

2 仿真結(jié)果及分析

2.1 平直地面

(1)迎角對翼型升力大小的影響

圖6 迎角對升力系數(shù)的影響Fig.6 Influence of angle of attack on lift coefficient

(2)平板翼型離平直地面高度對升力大小的影響

圖7 翼型高度對升力系數(shù)的影響Fig.7 Influence of airfoil height on lift coefficient

2.2 彎曲地面

(1)地面彎曲半徑對翼型升力系數(shù)的影響

圖8 山谷彎曲半徑對升力系數(shù)的影響Fig.8 Influence of hollow radius on lift coefficient

由圖5(b)可見,當(dāng)流體流過山谷彎曲地面的上表面時(shí),流場的速度方向沿著地面的形狀,先降低后爬升。因此,當(dāng)平板翼型位于山谷的右邊時(shí),其迎角由于流場速度的上升而變大,從而導(dǎo)致升力的增加。相反地,當(dāng)平板翼型位于山谷的左邊時(shí),其迎角由于流場速度的下降而變小,從而導(dǎo)致升力減小。而且,由于山谷會(huì)增大周圍流場空間,將導(dǎo)致山谷周圍的流場速度降低,進(jìn)而也引起升力減小。

(2)平板翼型離彎曲地面高度對升力大小的影響

圖9 山坡彎曲半徑對升力系數(shù)的影響Fig. 9 Influence of hill radius on lift coefficient

圖10 翼型高度在山谷情況下對升力系數(shù)的影響Fig.10 Influence of airfoil height on lift coefficient for hollow

圖11 翼型高度在山坡情況下對升力系數(shù)的影響Fig.11 Influence of airfoil height on lift coefficient for hill

3 結(jié)束語

本文采用復(fù)變函數(shù)理論、保角變換法和離散渦法,對平板翼型靠近平直和彎曲地面時(shí)引起升力變化的情況進(jìn)行了研究。由研究結(jié)果可知,當(dāng)平板翼型飛過平直地面時(shí),其升力由于受到地面影響會(huì)相對增大,并且升力增量隨著迎角的增大而減小,但是當(dāng)翼型的迎角增大至平衡迎角后,翼型升力反而比其在對應(yīng)無邊界流中的升力還要小。平板翼型越接近彎曲地面,其升力波動(dòng)越嚴(yán)重,而且,當(dāng)山谷或山坡地面的彎曲半徑越大時(shí),其對平板翼型的升力擾動(dòng)也越大。如果平板翼型以較小迎角飛過彎曲地面的上表面時(shí),其可能會(huì)產(chǎn)生負(fù)升力現(xiàn)象,加大了飛行器起飛或降落時(shí)墜落的危險(xiǎn)性。雖然本文僅對實(shí)際工程中使用較少的平板翼型進(jìn)行了研究,但研究結(jié)果對將來地效飛行器的研究具有一定的參考價(jià)值。

[1] 林海.飛機(jī)地面效應(yīng)及其試飛方法[J].飛行力學(xué),1996,14(4):52-58.

[2] 秦緒國,劉沛清,屈秋林,等.三維多段機(jī)翼地面效應(yīng)數(shù)值模擬[J].航空學(xué)報(bào),2011,32(2):257-264.

[3] Wieselsberger C.Wing resistance near the ground[R].NASA TM-77,1922.

[4] Tomotika S,Nagamiya T,Takemouti Y.The lift on a flat plate placed near a plane wall with special reference to the effect of the ground upon the lift of a monoplane aerofoil[J].Proceedings of Aeronautical Research Institute,1933,97(8):1-60.

[5] Ivanteeva L G,Konovalov S I,Pavlovets G A.Calculation of aerodynamic characteristics of the foil near the ground for given magnitudes of its geometric parameters[J].Proceedings of Central Aerohydrodynamics Institute,1980,11(2):1-9.

[6] Terentiev A G,Kartuzova T V.Numerical investigation of flow over airfoil near the ground[J].Proceedings of the National Academy of Sciences and Arts of the Chuvash Republic,1996,(6):94-104.

[7] 張亮,陳萌萌,吳德銘.二維地效翼氣動(dòng)力的計(jì)算和性能研究[J].船舶力學(xué),2000,4(2):1-5.

[8] 李盾, 王義寧.地效飛行器三維地面效應(yīng)的數(shù)值模擬[J].飛行力學(xué),2006,24(2):18-21.

[9] 王中清,白俊強(qiáng),王元元.飛行器地面效應(yīng)定常/非定常流場數(shù)值模擬[J].科學(xué)技術(shù)與工程,2012,12(5):1073-1075.

[10] Houghton E L,Carpenter P W.Aerodynamics for engineering students (5th edition)[M].Great Britain,Oxford:Butterworth-Heinemann,2003:160-171.

[11] 周云龍,郭婷婷.高等流體力學(xué)[M].北京:中國電力出版社,2008:26-32;67-75.

Numericalsimulationandanalysisonaerodynamicsofflatplateairfoilingroundeffect

CHEN Shi-shi, XIONG Fen-fen

(Key Laboratory of Dynamics and Control of Flight Vehicle, Ministry of Education, BIT, Beijing 100081, China)

At present, only the flat ground surface problem is investigated in most studies of ground effect, which is not effective for practical applications due to the universal existence of curved ground surface. In this paper, the lift calculation of a two-dimensional flat plate airfoil placed near the ground surface is investigated, including flat ground surface, hollow and hill curved ground surface. In order to calculate the complex velocity in the flow field, the complex variable function theory is employed. The conformal mapping method is applied to transform the flat plate flow to a flow around a circle cylinder. Moreover, the discrete vortices method is modified to convert the curved ground surface to a series of discrete boundary elements. It is shown that, when the flat plate airfoil flies close to the flat ground, its lift increases obviously. In addition, the lift fluctuation of the flat plate airfoil located near a hollow or a hill is considerably serious, which should be considered during the design of ground effect vehicle.

ground effect; curved ground surface; complex variable function; conformal mapping; discrete vortices

V211.4

A

1002-0853(2013)06-0486-05

2013-03-22;

2013-06-16; < class="emphasis_bold">網(wǎng)絡(luò)出版時(shí)間

時(shí)間:2013-10-22 14:13

國家自然科學(xué)基金資助 (51105034;10972034)

陳世適(1987-)，男，廣東湛江人，博士研究生，研究方向?yàn)轱w行器總體設(shè)計(jì)；

(編輯：姚妙慧)