• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    流體推力矢量技術(shù)研究綜述

    2017-09-15 09:09:41肖中云陳作斌
    實(shí)驗(yàn)流體力學(xué) 2017年4期
    關(guān)鍵詞:喉道射流矢量

    肖中云, 江 雄, 牟 斌, 陳作斌

    (中國空氣動力研究與發(fā)展中心 計算空氣動力研究所, 四川 綿陽 621000)

    流體推力矢量技術(shù)研究綜述

    肖中云*, 江 雄, 牟 斌, 陳作斌

    (中國空氣動力研究與發(fā)展中心 計算空氣動力研究所, 四川 綿陽 621000)

    流體推力矢量技術(shù)不采用機(jī)械偏轉(zhuǎn),以流動控制方式實(shí)現(xiàn)推力轉(zhuǎn)向,有望成為一種更加高效的推力矢量控制方法。目前實(shí)現(xiàn)流體推力矢量的主要方法有激波矢量法、雙喉道方法、逆流控制方法和同向流方法等,對以上方法選擇具有共性的計算與試驗(yàn)數(shù)據(jù),對噴管的推力矢量效率、推力損失和流量系數(shù)進(jìn)行了對比分析。結(jié)果表明激波矢量方法、雙喉道方法和逆流方法能夠在大落壓比范圍內(nèi)(NPR=1.89~10)實(shí)現(xiàn)推力矢量控制,并且具有俯仰/偏航耦合甚至多軸控制的潛力。相比激波矢量法和逆流方法,雙喉道和同向流方法在減少推力損失和提高矢量效率上占有優(yōu)勢,不足之處是雙喉道方法對喉道進(jìn)行控制限制了流量系數(shù),而同向流方法的適用落壓比范圍受到嚴(yán)重限制。為尋求更加高效的矢量噴管技術(shù),國內(nèi)外相繼發(fā)展了多種新概念流體推力矢量方法,對每種方法的控制原理、潛在優(yōu)勢和存在的問題挑戰(zhàn)進(jìn)行了探討,新方法著眼于從噴流出口下游進(jìn)行控制,對主流的干擾很小,值得深入研究,同時也為流體推力矢量的下一步研究方向提供了借鑒參考。

    噴管;推力矢量;控制;效率;二次流

    0 引 言

    流體推力矢量技術(shù)是指針對噴氣推進(jìn)方式,在不采用噴口機(jī)械偏轉(zhuǎn)前提下,通過流動控制手段改變和控制發(fā)動機(jī)主噴流方向,以此改變發(fā)動機(jī)的推力方向,獲得飛行控制所需的控制力和控制力矩。推力矢量技術(shù)由發(fā)動機(jī)推力分量產(chǎn)生直接控制力,可以直接對飛行器姿態(tài)進(jìn)行控制或者取代部分原有舵面的控制功能,實(shí)現(xiàn)原有飛機(jī)無法做到的機(jī)動或大迎角飛行能力。大量研究表明,在航空推進(jìn)系統(tǒng)中采用流體推力矢量控制在技術(shù)上是可行的,它通過氣流間的相互作用實(shí)現(xiàn)矢量偏轉(zhuǎn)和喉道調(diào)節(jié),結(jié)構(gòu)簡單、重量和造價大大降低,因此有望取代現(xiàn)有復(fù)雜的機(jī)械推力轉(zhuǎn)向裝置,帶來推力矢量技術(shù)的革命性變化。

    傳統(tǒng)推力矢量技術(shù)的實(shí)現(xiàn)主要采用機(jī)械偏轉(zhuǎn)的方式,主要包括折流板[1]、二元矢量噴管和軸對稱矢量噴管[2](AVEN)幾種方案。推力矢量技術(shù)在高性能飛機(jī)上有著很大的應(yīng)用潛力,包括提高飛機(jī)的機(jī)動性、飛行性能、生存能力以及隱身性能等各個方面。不過,對最新一代飛機(jī)來說這些優(yōu)勢還沒有得到充分的體現(xiàn),主要原因是額外增加的推力轉(zhuǎn)向裝置顯著增加了飛機(jī)重量、機(jī)構(gòu)復(fù)雜程度以及造價。機(jī)械系統(tǒng)的過于復(fù)雜一定程度上阻礙了這項(xiàng)技術(shù)在飛機(jī)上的應(yīng)用。另一方面,由于具有更好的氣動效率,飛翼/無尾布局可能成為未來飛行器的一種主要布局方式,這種布局取消了平尾、垂尾等控制舵面,設(shè)計師希望推力矢量能夠成為一種高效的控制方式,在推力損失不大、控制響應(yīng)不降的情況下,補(bǔ)充一種新的穩(wěn)定和操控手段。流體推力矢量具有達(dá)到這一目標(biāo)的潛在能力,理想設(shè)計狀況下這種方式具有較高的運(yùn)動、控制可靠性,不僅可以為高機(jī)動飛機(jī)提供足夠的俯仰控制力矩,也可以為無尾/飛翼布局飛機(jī)的飛行安定性,提供頻繁工作的三軸穩(wěn)定性控制力矩。

    本文對目前發(fā)展的主要流體推力矢量方法進(jìn)行了回顧,分別介紹了各方法的控制原理和適用范圍,在同一評價指標(biāo)下對推力矢量的矢量效率和推力損失進(jìn)行了比較,探討了實(shí)現(xiàn)工程應(yīng)用還需要解決的問題。在此基礎(chǔ)上,介紹了3種目前還處于實(shí)驗(yàn)室階段的新概念流體推力矢量方法,討論了這些方法的潛在優(yōu)勢和存在的問題挑戰(zhàn),為下一步流體推力矢量研究提供借鑒。

    1 主要流體推力矢量方法

    文獻(xiàn)[3-5]對幾種流體推力矢量的發(fā)展歷程進(jìn)行了很好的綜述,這里主要對比分析幾種方法的優(yōu)劣。20世紀(jì)80年代末,美國航空航天局蘭利研究中心(NASA LaRC)發(fā)起了固定形狀噴管的流體推力矢量技術(shù)研究,前后提出了3種不同的推力矢量控制方法,分別是激波矢量法、喉道偏斜法和逆流方法。3種方式的控制原理如圖1所示,分別利用了激波、喉道聲速線和抽吸負(fù)壓的原理實(shí)現(xiàn)偏轉(zhuǎn)。同向流控制在歐洲國家進(jìn)行了較多研究,是一種利用Coanda效應(yīng)實(shí)現(xiàn)偏轉(zhuǎn)的控制方法,適用于對低速主流進(jìn)行控制。以上控制技術(shù)按控制位置不同又可以分為3類,分別是對喉道施加控制、對噴管擴(kuò)張段的控制和對噴流出口下游進(jìn)行控制。

    1.1 激波矢量法

    激波矢量法[6-10](SVC:Shock Vectoring Control)利用了流體經(jīng)過斜激波后方向發(fā)生偏轉(zhuǎn)的原理,如圖1(a)所示,通過注入二次流動與噴管超聲速氣流干擾形成激波,激波越強(qiáng),偏轉(zhuǎn)角度越大。

    NASA蘭利中心(LaRC)從20世紀(jì)80年代開始開展激波矢量法控制技術(shù)研究,首先是在二維收斂-擴(kuò)張型噴管上進(jìn)行試驗(yàn),實(shí)現(xiàn)了俯仰方向的單軸推力矢量控制。隨后又通過在橫側(cè)方向上增加科恩達(dá)效應(yīng)吹氣控制,實(shí)現(xiàn)俯仰和偏航方向的推力矢量控制[6]。NASA LaRC聯(lián)合普惠公司開展了多軸推力矢量噴管(MATV)研究[8],將二次射流縫布置在噴管擴(kuò)張段的上下壁面和側(cè)壁上,分別對射流產(chǎn)生俯仰和偏航方向的控制。試驗(yàn)采用高壓冷氣源的方式,二次流總壓與主流總壓的比值變化范圍在0~1之間,二次流流量與主流的比值為0~0.12之間,試驗(yàn)獲得了俯仰或偏航單個方向上15°的偏角,俯仰和偏航同時偏轉(zhuǎn)的情況下獲得了10°的偏角。其他嘗試還包括軸對稱收斂-擴(kuò)張噴管[9]的多軸推力矢量控制研究,目的是實(shí)現(xiàn)更好的推力矢量和噴管內(nèi)流性能。

    在激波矢量法中,二次流壓比作為輸入控制參數(shù),理論上二次流壓比越大得到的推力矢量角就越大。在流動形態(tài)上,隨二次流壓比增大,二次流注入位置后側(cè)的流動分離類型從閉式分離向開式分離過渡。當(dāng)二次流壓比增大到一定值以后,所產(chǎn)生的誘導(dǎo)激波會碰到另一側(cè)壁面后發(fā)生反射,起到減小推力矢量角的反向作用。激波矢量法的主要特點(diǎn)是可以產(chǎn)生較大的矢量角和較高的矢量效率(設(shè)計狀態(tài)下引射量為1%時矢量偏角3.3°),但是在超聲速區(qū)使氣流發(fā)生偏轉(zhuǎn),產(chǎn)生大的矢量角時要求噴管工作在過膨脹的狀態(tài),因此SVC不僅有激波損失,還存在過膨脹損失,其推力系數(shù)一般在0.86~0.94之間。文獻(xiàn)[9]提出為解決激波矢量方法的過膨脹流動問題,利用流體注入方法改變喉道面積和噴管膨脹比,減少因過膨脹流動帶來的推力損失,但必然會帶來控制系統(tǒng)的復(fù)雜性。

    1.2 喉道偏斜法

    喉道偏斜法[11]采用在噴管喉道注入二次流的方法控制氣流偏轉(zhuǎn),避免了激波矢量法由于產(chǎn)生激波帶來的流動損失,典型條件下具有更大的推力系數(shù),但獲得的推力矢量角小于激波矢量法。2003年Deere等提出的雙喉道噴管DTN(Dual Throat Nozzle)的概念[12],如圖1(b)所示,NASA隨后開展了比較廣泛的數(shù)值計算和地面試驗(yàn)研究[13-16],噴管形狀涉及到二維平面和軸對稱外形,結(jié)果表明雙喉道方法可以在較大壓比范圍內(nèi)實(shí)現(xiàn)推力矢量。經(jīng)過對雙喉道之間區(qū)域流道的優(yōu)化設(shè)計后,獲得了比原來喉道偏移法更大的推力矢量效率(引射量為1%時矢量偏角約4°), 并且推力系數(shù)也較高(0.92~0.96之間) , 實(shí)現(xiàn)了較好的綜合性能。Deere[12]研究后認(rèn)為,雙喉道方法主流偏轉(zhuǎn)與聲速線偏離并無直接關(guān)系,對上下游喉道之間的分離區(qū)進(jìn)行控制才是產(chǎn)生大的矢量偏轉(zhuǎn)的真正原因。為了產(chǎn)生大的推力矢量角,應(yīng)當(dāng)使二次流的入射角度盡量與主流反向,保持上下游喉道面積相等,同時增大空腔收縮段的角度。為了提高矢量噴管的效率,文獻(xiàn)[17]對二維雙喉道射流推力矢量噴管的靜特性進(jìn)行了參數(shù)化數(shù)值研究,比較了雙喉道噴管的空腔長度、空腔擴(kuò)張角、空腔收斂角、上游喉道高度等設(shè)計參數(shù)對雙喉道射流矢量噴管的影響。文獻(xiàn)[18]采用徑向基神經(jīng)網(wǎng)絡(luò)( RBF)和粒子群算法( PSO)對設(shè)計參數(shù)進(jìn)行尋優(yōu),找出了更好的雙喉道氣動矢量噴管設(shè)計參數(shù)組合。在主流壓比NPR=4,二次流流量3%條件下,經(jīng)過優(yōu)化設(shè)計以后的噴管的矢量效率從4.039提高到5.147,同時推力損失系數(shù)基本保持不變。

    雙喉道方法的成功之處在于在噴管喉道處進(jìn)行二次流控制,減少了激波流動損失;同時去掉了噴管擴(kuò)張段,讓氣流膨脹發(fā)生在管外,減少了過膨脹流動帶來的損失。雙喉道方法在提高推力系數(shù)的同時也帶來了新的問題,由于在噴管喉道引入二次流,減小了喉道面積,造成流動出現(xiàn)壅塞而流量系數(shù)減小。文獻(xiàn)數(shù)據(jù)表明,雙喉道噴管在典型控制狀態(tài)下流量系數(shù)僅為控制前的84%,這就意味著,在施加推力矢量控制以后,由于流量減小導(dǎo)致發(fā)動機(jī)推力損失,同時為了維持發(fā)動機(jī)工作穩(wěn)定,還必須考慮如何處理多余流量的問題。

    1.3 逆流控制方法

    如圖1(c)所示,逆流控制方法[19-24]在噴管出口截面的下游加一個外套,形成反向流動的通道,在需要主流偏轉(zhuǎn)時,啟動真空泵抽吸系統(tǒng)形成負(fù)壓,使主氣流偏轉(zhuǎn)產(chǎn)生側(cè)向力。逆流方法在出口下游增加固壁延伸面,抽吸噴管主流與固壁面之間的氣體形成低壓區(qū),使主流發(fā)生偏轉(zhuǎn)。由于是對射流在切向上進(jìn)行控制,具有推力損失小的特點(diǎn),推力系數(shù)約在0.92~0.97之間,被認(rèn)為是一種有前途的控制方法。Strykowski[21]提出,當(dāng)逆流速度與主流速度的比值超過臨界值時(即U2/U1<-0.14),逆流剪切層由對流不穩(wěn)定轉(zhuǎn)變?yōu)榻^對不穩(wěn)定,剪切層內(nèi)的渦結(jié)構(gòu)和湍動特性增強(qiáng),是射流偏轉(zhuǎn)的內(nèi)在機(jī)制。Hunter[24]在研究中認(rèn)為,不管剪切流是同向還是逆向,只要注入的動量通量相等,所產(chǎn)生的推力矢量效果就應(yīng)該是一樣的,Hunter認(rèn)為套筒壁在吸氣側(cè)上的壓力分布是驅(qū)動推力轉(zhuǎn)向的關(guān)鍵因素,與二次流的方向無關(guān)。Flamm[20]的實(shí)驗(yàn)細(xì)致地測量了矢量角和二次流率的關(guān)系,證明同向和逆向二次流都能獲得推力轉(zhuǎn)向。

    文獻(xiàn)[22]對逆流推力矢量的外流影響進(jìn)行了試驗(yàn),其中噴管落壓比對應(yīng)的完全膨脹馬赫數(shù)為Ma=1.4,來流馬赫數(shù)對應(yīng)了Ma=0.3、0.5和0.7這3種條件。試驗(yàn)表明,在來流條件下,逆流推力矢量方法可以比較容易地實(shí)現(xiàn)矢量角22°的偏轉(zhuǎn),在相同控制條件下,隨來流馬赫數(shù)的升高,噴管的推力矢量有微弱的減小。逆流推力矢量的工程應(yīng)用還有許多問題亟待解決,比如采用什么樣的吸氣裝置、流動遲滯效應(yīng)、大尺寸的后緣板與機(jī)身難以融合等,更為嚴(yán)重的是,主流在偏角較大時容易附著于控制面上,導(dǎo)致推力方向的突然變化,給飛行控制帶來風(fēng)險。

    1.4 同向二次流控制方法

    為了尋找新的控制途徑,英國曼徹斯特大學(xué)Mason等[25]開展了同向二次流控制的射流推力矢量研究。如圖1(d)所示,該方法利用與主流方向平行的二次射流與壁面Coanda效應(yīng),使主流方向發(fā)生偏轉(zhuǎn),在主流馬赫數(shù)Ma=0.33條件下,獲得了推力矢量控制效果,研究認(rèn)為該方案在主流速度為低亞聲速時可以取得很好的控制效率。Banazadeh等[26]在小型渦噴發(fā)動機(jī)上開展試驗(yàn)研究,在發(fā)動機(jī)低轉(zhuǎn)速下(5.8×104~7.8×104r/min)實(shí)現(xiàn)了推力矢量性能,試驗(yàn)噴管為軸對稱外形,具有多軸推力矢量的功能,但隨發(fā)動機(jī)轉(zhuǎn)速升高,射流速度增加,該噴管的矢量控制能力迅速下降。該控制技術(shù)被應(yīng)用到英國無尾飛行器創(chuàng)新研究項(xiàng)目FLAVIIR[27-29](Flapless Air Vehicle Integrated Industrial Research)中,在該項(xiàng)目中Gill[30]探討了從發(fā)動機(jī)中引氣的方式,作為機(jī)翼環(huán)量控制和流體推力矢量控制的二次流氣源。后來這些技術(shù)都集成在驗(yàn)證飛行器Demon上,在2010年進(jìn)行了演示飛行,但是對流體推力矢量的實(shí)際控制能力缺少報道。同向二次流也被用來嘗試對高速射流進(jìn)行控制,Heo[31]通過在噴管下游設(shè)置更大的產(chǎn)生Coanda效應(yīng)的壁面,采用同向二次流的射流控制方法,實(shí)現(xiàn)了對超聲速主流的射流偏轉(zhuǎn)控制。

    2 矢量能力與效率

    評價一種流體推力矢量方法包括了噴管的單軸/多軸控制能力、最大矢量角、引氣量、推力損失大小等多個方面,在控制效率上主要用以下3個指標(biāo)進(jìn)行衡量,一是推力比,即實(shí)際推力與理想等熵推力的比值;二是推力矢量效率,即單位引氣量所產(chǎn)生的推力矢量角大小;三是流量系數(shù),即實(shí)際流量與理想等熵流量的比值。上述指標(biāo)的理論公式如下所示。

    (1) 推力系數(shù)[32]:

    (2) 流量系數(shù)[16]:

    其中wp和ws分別為主流和次流的質(zhì)量流量,wi,p為理想主流流量。理想質(zhì)量流率由噴管喉道面積和總壓、總溫等參數(shù)決定的。

    (3) 推力矢量的效率:

    圖2對比了不同流體推力矢量控制方法的推力系數(shù),其中數(shù)據(jù)來源于公開發(fā)表的試驗(yàn)數(shù)據(jù)[3, 7-9, 15-16, 20]。由于同向二次流控制方法主要適用于控制低速主流,與其他幾種方法的適用范圍不同,這里沒有列入圖中進(jìn)行比較。可以看到,雙喉道噴管在低落壓比的時候推力系數(shù)較高,而激波矢量法和逆流控制方法在落壓比較大的時候效率較高。就雙喉道噴管來說,噴管設(shè)計壓比為1.89,出口速度為聲速,在大落壓比的情況下流動欠膨脹造成推力損失。對于激波矢量法來說,為了使主流能夠偏轉(zhuǎn),噴管的設(shè)計壓比大于實(shí)際流動壓比,噴流處于過膨脹流動,流動壓比越小過膨脹程度就越大,造成較大的推力損失。

    Fig.2 Comparison of thrust coefficients for three control methods

    圖3顯示的是不同推力矢量方法的矢量效率比較,其中逆流矢量法主要利用真空泵抽吸原理,不同于二次流注入方法,此處不便于進(jìn)行橫向比較。這里矢量效率用吹氣量占主流1%時得到的矢量角表示,可以看到雙喉道方法具有較高的矢量效率,在壓比NPR=1.89~10范圍內(nèi),矢量效率在3.4°/1%~5.2°/1%之間變化。激波矢量法的矢量效率要整體低于雙喉道噴管,從文獻(xiàn)收集的數(shù)據(jù)顯示矢量效率覆蓋0.8°/1%~3.3°/1%這個區(qū)間。圖4給出的是雙喉道方法與激波矢量法的流量系數(shù)比較,從圖中可以看到,激波矢量法的流量系數(shù)在0.96以上,遠(yuǎn)遠(yuǎn)高于雙喉道方法。其中原因是激波矢量法在噴管擴(kuò)張段進(jìn)行控制,不改變噴管的喉道面積,對流量的影響較小;相反雙喉道法是在噴管喉道位置注入二次流,減小了噴管的最小流通面面積,導(dǎo)致流量減小。由此可見,雙喉道方法實(shí)現(xiàn)推力矢量是以改變喉道面積為代價,越是要獲得大的矢量偏角,二次流注入越強(qiáng),噴管的實(shí)際流量就越小,發(fā)動機(jī)產(chǎn)生的實(shí)際推力也就越小。

    Fig.3 Comparison of thrust vectoring efficiencies between DTN and SVC

    表1給出了3種流體推力矢量噴管的性能比較,可以看到,3種方法中激波矢量法存在較大的推力損失,矢量效率較低,但是流量系數(shù)較大;雙喉道方法獲得了較大的矢量效率和推力比,唯一存在的不足是流量系數(shù)受到很大的限制,對發(fā)動機(jī)主流的影響很大;逆流方法在3個方面都有較高的指標(biāo),但存在抽吸系統(tǒng)過于復(fù)雜的不足。

    表1 3種流體推力矢量控制方法的比較Table 1 Comparison of three fluidic thrust vectoring methods

    流體推力矢量不同于機(jī)械偏轉(zhuǎn)的另一個特點(diǎn)是,實(shí)際外流條件下的推力矢量角小于靜態(tài)條件下測得的矢量角。文獻(xiàn)[7]開展了激波矢量法的外流計算,自由來流馬赫數(shù)Ma=0.3~1.2,噴管的落壓比NPR=3.6~7.2,二次流的壓比pS/pt=0.6~1.0,結(jié)果表明,來流影響會減弱推力矢量性能,減少量從Ma=0.9(NPR=6)的1.5°到Ma=0.7(NPR=5.2)的2.9°,推力系數(shù)損失則從4%(對應(yīng)壓比3.6)到0.8%(對應(yīng)壓比7.2)。解決從哪里引氣是流體推力矢量面臨的另一難題,一個比較可行的辦法是從發(fā)動機(jī)的壓縮機(jī)段引氣。WANG[10]采用三維CFD方法模擬了激波矢量噴管航空發(fā)動機(jī)耦合的氣動性能,結(jié)果表明采用風(fēng)扇出口引氣的方式優(yōu)于從一級壓力機(jī)引氣,后者對發(fā)動機(jī)性能影響嚴(yán)重,前者從風(fēng)扇出口引氣,當(dāng)推力矢量角達(dá)到19°時,推力減少量在17%以內(nèi)。

    3 新概念流體推力矢量方法

    傳統(tǒng)流體推力矢量方法目前還存在一些難以解決的瓶頸問題,導(dǎo)致該技術(shù)未能在工程實(shí)踐中推廣使用。下面介紹幾種完全不同于傳統(tǒng)方法的流體推力矢量控制技術(shù),作為流體推力矢量的一種新的嘗試,可以為下一步推力矢量發(fā)展提供借鑒。

    3.1 合成射流控制方法

    合成射流技術(shù)是一種基于小尺度擾動引起大的宏觀效應(yīng)的主動流動控制技術(shù)。Smith和Glezer[33]在1998年用合成射流控制了噴口直徑是自身25倍的噴流轉(zhuǎn)向,研究中合成射流方向與噴流方向平行,在合成射流不作用時,噴流沿噴口中心線對稱分布,在合成射流作用時,射流致動器所產(chǎn)生的渦環(huán)引起主噴流上方的壓力降低,導(dǎo)致主噴流向這一側(cè)發(fā)生偏轉(zhuǎn)。本文對合成射流控制主流偏轉(zhuǎn)進(jìn)行了非定常數(shù)值模擬,圖5顯示的是合成射流控制的渦量云圖,其中時間t=0~T/2區(qū)間為射流的噴出階段,可以看到射流在經(jīng)過特殊設(shè)計的噴口噴出以后,右側(cè)渦量分布非常集中,左側(cè)渦量發(fā)展受到抑制,在誘導(dǎo)速度作用下繞右側(cè)旋渦旋轉(zhuǎn)。圖6顯示的是合成射流控制射流偏轉(zhuǎn)的瞬時速度云圖,其中主流速度7m/s,零質(zhì)量射流的峰值速度34m/s,動量比值為1∶10,獲得的射流轉(zhuǎn)向角約為30°。文獻(xiàn)[34-35]對不同出口構(gòu)型合成射流激勵器開展的射流矢量控制進(jìn)行了研究,認(rèn)為低壓、合成射流動量以及合成射流對主流流體的卷吸是控制主流偏轉(zhuǎn)的主要原因。

    為了提高合成射流控制射流偏轉(zhuǎn)的效果,Park等[36]對出口帶有擴(kuò)張段的圓管流進(jìn)行了實(shí)驗(yàn)研究,其中激勵器安裝在圓管出口與擴(kuò)張段的交接位置。實(shí)驗(yàn)結(jié)果表明,出口下游擴(kuò)張段可以起到增強(qiáng)摻混的作用,同時有控制射流轉(zhuǎn)向的效果,但是實(shí)驗(yàn)僅限于對低速主流進(jìn)行控制。合成射流控制射流轉(zhuǎn)向試驗(yàn)向人們展示了合成射流控制的奇妙之處,但是將這一技術(shù)應(yīng)用于推力轉(zhuǎn)向還有很遠(yuǎn)的距離,目前成功的射流轉(zhuǎn)向試驗(yàn)還局限于小速度尺度和低雷諾數(shù)情形,要實(shí)現(xiàn)對更大尺度的發(fā)動機(jī)高速尾噴流的控制還存在很大的難度。

    3.2 Coanda效應(yīng)噴管

    Coanda效應(yīng)噴管在噴管出口下游進(jìn)行矢量偏轉(zhuǎn)控制,通常作法是在噴管出口增加后緣板,利用流體自身的附壁效應(yīng),在較小輸入能量的前提下實(shí)現(xiàn)高速大動量射流的偏轉(zhuǎn),也可以簡稱這種方式為后緣板控制。2013年歐盟資助開展了空氣柯恩達(dá)效應(yīng)高效射流矢量噴管計劃ACHEON(Aerial Coanda High Efficiency Orienting-jet Nozzle project)[37-39],該研究計劃提出利用高速射流在凸表面的Coanda效應(yīng)和等離子體加速流體延遲分離的作用,提高矢量噴管的偏轉(zhuǎn)效率。研究內(nèi)容包括對3種物理效應(yīng)進(jìn)行研究,分別是高速射流的混合效應(yīng)、Coanda效應(yīng)和靜電場控制。文獻(xiàn)[38]設(shè)計了一種雙主流的合成射流噴管,如圖7所示,2股主流在噴管喉道前進(jìn)行匯合并經(jīng)過收縮-擴(kuò)張噴管排出,當(dāng)2股主流存在速度差的時候,合成射流將向速度高的一側(cè)偏轉(zhuǎn)。

    在后緣Coanda壁面采用等離子體控制被認(rèn)為是未來流動控制的一項(xiàng)關(guān)鍵技術(shù),未來有可能用于推力矢量控制。吳云等人[40]對等離子體激勵器種類及流動控制應(yīng)用進(jìn)行了綜述,指出了等離子體氣動激勵在抑制流動分離、激波控制、管道流動方面的諸多應(yīng)用,等離子體產(chǎn)生氣動激勵的物理作用原理包括了動力效應(yīng)、沖擊效應(yīng)和物性改變。Samimy研究組[41]開展了等離子體對高速高雷諾數(shù)管道射流控制的研究,表明等離子體可以在有效促進(jìn)流動摻混,增大管道出口流動的湍動能,縮短主流核心區(qū)長度,降低排氣噪聲等方面發(fā)揮作用。

    3.3 引射效應(yīng)矢量噴管

    對高速射流而言,一個重要的特點(diǎn)是具有很強(qiáng)的引射效應(yīng),射流速度越大,對周圍氣體的引射作用越強(qiáng)?;谏淞鞯倪@種特性,本文提出了一種利用引射效應(yīng)的推力矢量控制方法[42]。如圖8所示,在主噴管出口下游設(shè)置有2塊對稱的后緣板,在主噴流的引射作用下,兩側(cè)氣體被不斷抽走,同時又從下游以回流形式得到補(bǔ)充。此時主流存在一個非常不穩(wěn)定的狀態(tài),既可能保持中立狀態(tài),射流方向保持軸線方向不變,又可能因擾動而發(fā)生向上或者向下的偏轉(zhuǎn),我們把射流的這種行為描述為具有雙穩(wěn)特性[42-44]。在把射流變得十分敏感以后,加入人為擾動可以很容易地實(shí)現(xiàn)射流偏轉(zhuǎn)。文獻(xiàn)[42]采用改變通氣量的方式實(shí)現(xiàn)了射流的上下偏轉(zhuǎn)控制,但是不足是射流容易附著到固體壁面,導(dǎo)致射流方向的突然變化。

    該控制技術(shù)被用于對低速主流進(jìn)行控制[44],通過改變控制縫入口面積實(shí)現(xiàn)了對主射流偏轉(zhuǎn)的比例控制,在測試條件下主射流最大偏轉(zhuǎn)角可達(dá)19°。控制曲線分為了2個區(qū)域,敏感區(qū)和遲鈍區(qū)。敏感區(qū)的控制曲線近似線性,斜率較大,范圍約為±15°;在遲鈍區(qū)射流靠近壁面,射流方向?qū)刂频捻憫?yīng)比較遲緩。由于該方法采用的是被動控制手段,不需要任何能量輸入,相對其他方法而言具有自己獨(dú)特的優(yōu)勢。

    上述3種新概念流體推力矢量方法都試圖從噴管出口下游進(jìn)行控制,對主流或者說發(fā)動機(jī)的干擾較小,代表了流體推力矢量研究的一個新方向。目前這些方法還僅限于在設(shè)計條件下進(jìn)行使用,不能勝任大范圍變化壓比等復(fù)雜流動條件,因此這些方法在走向?qū)嶋H應(yīng)用之前還有許多障礙需要克服。

    4 結(jié)束語

    流體推力矢量技術(shù)由于其潛在的應(yīng)用價值得到了國內(nèi)外研究機(jī)構(gòu)的重視,目前這項(xiàng)技術(shù)還沒有得到工程應(yīng)用,主要原因是制約噴管性能的影響參數(shù)較多,這些參數(shù)之間的相互作用關(guān)系復(fù)雜,控制難度較大。未來隨著微型傳感器技術(shù)、自動控制技術(shù)的發(fā)展,流體推力矢量也一定會走上前臺,在工程實(shí)踐中發(fā)揮作用。本文主要結(jié)論如下:

    (1) 未來流體推力矢量的發(fā)展不僅面臨解決飛行器的機(jī)動問題,實(shí)現(xiàn)單軸(俯仰或偏航方向)或多軸推力矢量控制,還要面臨解決無尾布局飛行器的安定性問題,為飛行穩(wěn)定提供頻繁工作的三軸穩(wěn)定性控制力矩,為實(shí)現(xiàn)真正的高氣動效率飛行器布局設(shè)計提供技術(shù)支撐。

    (2) 流體推力矢量控制受多種因素影響,除形狀參數(shù)外,還包括噴管落壓比,二次流流量,外流速度等流動參數(shù),這些參數(shù)影響噴管的推力比、矢量效率、流量系數(shù)大小。流體推力矢量走向工程應(yīng)用是一個與發(fā)動機(jī)、與飛行控制相結(jié)合的過程,需要解決各項(xiàng)系統(tǒng)參數(shù)之間的匹配與銜接問題,而計算機(jī)技術(shù)、先進(jìn)傳感器與飛行控制技術(shù)的發(fā)展將有助于解決這些問題。

    (3) 機(jī)械偏轉(zhuǎn)噴管在改變內(nèi)流道的同時也改變了噴管周圍外流的流動方向。目前流體推力矢量研究主要集中解決內(nèi)流的偏轉(zhuǎn)問題,對外流影響的研究較少,缺乏減少外流影響的措施。對目前所有流體推力矢量控制來說,都面臨外流速度增加導(dǎo)致噴管矢量效率下降的問題,如何解決這一問題也應(yīng)得到重視。

    [1]Bowers A H, Pahle J W. Thrust vectoring on the NASA F-18 high alpha research vehicle[R]. NASA Technical Memorandum 4771, 1996.

    [2]Kowal H J. Advances in thrust vectoring and the application of flow-control technology[J]. Canadian Aeronautics and Space Journal, 2002, 48(2): 145-151.

    [3]Deere K A. Summary of fluidic thrust vectoring research conducted at NASA langley research center[R]. AIAA-2003-3800, 2003.

    [4]連永久. 射流推力矢量控制技術(shù)研究[J]. 飛機(jī)設(shè)計, 2008, 28(2): 19-24.

    Lian Y J. Fluidic thrust vectoring techniques research[J]. Aircraft Design, 2008, 28(2): 19-24.

    [5]宋亞飛, 高峰, 何至林. 流體推力矢量技術(shù)[J]. 飛航導(dǎo)彈, 2010, (11): 72-75.

    [6]Wing D J. Static investigation of two fluidic thrust-vectoring concepts on a two dimensional convergent-divergent nozzle[R]. NASA Technical Memorandum 4574, 1995.

    [7]Deere K A. Computational investigation of the aerodynamic effects on fluidic thrust vectoring[R]. AIAA-2000-3598, 2000.

    [8]Giuliano V J, Wing D J. Static investigation of a fixed-aperture nozzle employing fluidic injection for multiaxis thrust vector control[R]. AIAA-1997-3149, 1997.

    [9]Wing D J, Giuliano V J. Fluidic thrust vectoring of an axisymmetric exhaust nozzle at static conditions[R]. FEDSM97-3228, 1997.

    [10]Shi J W, Wang Z X, Zhang X B, et al. Performance estimation for fluidic thrust vectoring nozzle coupled with aero-engine[R]. AIAA-2014-3771, 2014.

    [11]Miller D N, Yagle P J, Hamstra J W. Fluidic throat skewing for thrust vectoring in fixed geometry nozzles[R]. AIAA-99-16262, 1999.

    [12]Deere K A, Berrier B L, Flamm J D, et al. Computational study of fluidic thrust vectoring using separation control in a nozzle[R]. AIAA-2003-3803, 2003.

    [13]Deere K A, Flamm J D, Berrier B L, et al. Computational study of an axisymmetric dual throat fluidic thrust vectoring nozzle for a supersonic aircraft application[R]. AIAA-2007-5085, 2007.

    [14]Flamm J D, Deere K A, Berrier B L, et al. Experimental study of a dual-throat fluidic thrust-vectoring nozzle concept[R]. AIAA-2005-3503, 2005.

    [15]Flamm J D, Deere K A, Mason M L, et al. Design enhancements of the two-dimensional, dual throat fluidic thrust vectoring nozzle concept[R]. AIAA-2006-3701, 2006.

    [16]Flamm J D, Deere K A, Mason M L, et al. Experimental study of an axisymmetric dual throat fluidic thrust vectoring nozzle for supersonic aircraft application[R]. AIAA-2007-5084, 2007.

    [17]譚慧俊, 陳智. 二元雙喉道射流推力矢量噴管的數(shù)值模擬研究[J]. 航空動力學(xué)報, 2007, 22(10): 1678-1684.

    Tan H J, Chen Z. A computational study of 2-D dual-throat fluidic thrust-vectoring nozzles[J]. Journal of Aerospace Power, 2007, 22(10): 1678-1684.

    [18]吳正科, 楊青真, 施永強(qiáng), 等. 基于RBF和PSO的雙喉道氣動矢量噴管優(yōu)化設(shè)計[J]. 推進(jìn)技術(shù), 2013, 34(4): 451-456.

    Wu Z K, Yang Q Z, Shi Y Q, et al. Optimization design of the dual throat fluidic thrust vectoring nozzle based on RBF and PSO[J]. Journal of Propulsion Technology, 2013, 34(4): 451-456.

    [19]Hunter C A, Deere K A. Computational investigation of fluidic counterflow thrust vectoring[R]. AIAA-99-2669, 1999.

    [20]Flamm J D. Experimental study of a nozzle using fluidic counterflow for thrust vectoring[R]. AIAA-1998-3255, 1998.

    [21]Strykowski P J, Krothapalli A, Forliti D J. Counterflow thrust vectoring of supersonic jets[J]. AIAA Journal, 1996, 34(11):2306-2314.

    [22]Alvi F S, Strykowski P J. Forward flight effects on counterflow thrust vector control of a supersonic jet[J]. AIAA Journal, 2015, 37(2): 279-281.

    [23]Strykowski P J. An experimentallmodeling study of jet attachment during counterflow thrust vectoring[R]. NASA-CR-204436, 1996.

    [24]Hunter C A. Experimental, theoretical, and computqational investigation of separated nozzle flows[R]. AIAA-98-3107, 1998.

    [25]Mason M S, Crowther W J. Fluidic thrust vectoring for low observable air vehicles[R]. AIAA-2004-2210, 2004.

    [26]Banazadeh A, Saghafi F, Ghoreyshi M, et al. Multi-directional co-flow fluidic thrust vectoring intended for a small gas turbine[R]. AIAA-2007-2940, 2007.

    [27]Fielding J P, Smith H. FLAVIIR, an innovative university/industry research program for collaborative research and demonstration of UAV technologies[C]. 25th International Congress of the Aeronautical Sciences, 2006.

    [28]Abbasi A Y, Clarke A, Lawson C P, et al. Design and development of the eclipse and demon demonstrator UAVs[C]. 26th International Congress of The Aeronautical Sciences, 2008.

    [29]Fielding J P, Lawson C P, Pires R, et al. Development of the demon technology demonstrator UAV[C]. 27th International Congress of The Aeronautical Sciences, 2010.

    [30]Gill K, Wilde P, Gueroult R, et al. Development of an integrated propulsion and pneumatic power supply system for flapless UAVs[R]. AIAA-2007-7726, 2007.

    [31]Heo J Y, Yoo K H, Lee Y, et al. Fluidic thrust vector control of supersonic jet using co-flow injection[R]. AIAA-2009-5174, 2009.

    [32]Lamb M, Taylor J G, Frassinelli M C. Static internal performance of a two dimensional convergent divergent nozzle with external shelf[R]. NASA Technical Memorandum 4719, 1996.

    [33]Smith B L, Glezer A. Vectoring and small-scale motions effected in free shear flows usin synthetic jet actuators[R]. AIAA-97-0213, 1997.

    [34]夏智勛,羅振兵. 合成射流激勵器射流矢量控制的物理因素[J]. 應(yīng)用數(shù)學(xué)和力學(xué), 2007, 28(7): 811-823.

    Xia Z X, Luo Z B. Physical factors of a primary jet vectoring control using synthetic jet actuators[J]. Applied Mathematics and Mechanics, 2007, 28(7): 811-823.

    [35]羅振兵, 夏智勛. 合成射流技術(shù)及其在流動控制中應(yīng)用的進(jìn)展[J]. 力學(xué)進(jìn)展, 2005, 35(2): 221-234.

    Luo Z B, Xia Z X. Advances in synthetic jet technology and applications in flow control[J]. Advances in Mechanics, 2005, 35(2): 221-234.

    [36]Park L G, Seifert A. Periodic excitation for jet vectoring and enhanced spreading[J]. Journal of Aircraft, 2001, 38(3): 486-495.

    [37]Trancossi M, Dumas A, Vucinic D. Mathematical modeling of coanda effect[R]. SAE Technical Paper 2013-01-2195, 2013.

    [38]Trancossi M, Dumas A, Das S S, et al. Design methods of Coanda effect nozzle with two streams[J]. INCAS Bulletin, 2014, 6(1): 83-95.

    [39]Sunol A, Vucinic D. Numerical analysis and UAV application of the ACHEON thrust vectoring nozzle[R]. AIAA-2014-2046, 2014.

    [40]吳云, 李應(yīng)紅. 等離子體流動控制研究進(jìn)展與展望[J]. 航空學(xué)報, 2015, 36(2): 381-405.

    Wu Y, Li Y H. Progress and outlook of plasma flow control[J]. Acta Aeronautica Et Astronautica Sinica, 2015, 36(2): 381-405.

    [41]Sinha A, Alkandry H, Fischer M K, et al. The impulse response of a high-speed jet forced with localized arc filament plasma actuators[J]. Physics of Fluids, 2012, 24(12): 1-20.

    [42]肖中云, 顧蘊(yùn)松, 江雄, 等. 一種基于引射效應(yīng)的流體推力矢量新技術(shù)[J]. 航空學(xué)報, 2012, 33(11): 1967-1974.

    Xiao Z Y, Gu Y S, Jiang X, et al. A new fluidic thrust vectoring technique based on ejecting mixing effects[J]. Acta Aeronautica et Astronautica Sinica, 2012, 33(11): 1967-1974.

    [43]Ozgu M R, Stenning A H. Switching dynamics of bistable fluid amplifiers[R]. AD72383, 1971.

    [44]曹永飛, 顧蘊(yùn)松, 程克明, 等. 基于被動二次流的射流偏轉(zhuǎn)比例控制[J]. 航空學(xué)報, 2015, 36(3): 757-762.

    Cao Y F, Gu Y S, Cheng K M, et al. Proportional control of jet deflection with passive secondary flow[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(3): 757-762.

    (編輯:李金勇)

    Advances influidic thrust vectoring technique research

    Xiao Zhongyun*, Jiang Xiong, Mou Bin, Chen Zuobin

    (Computational Aerodynamics Institute, China Aerodynamics Research and Development Center, Mianyang Sichuan 621000, China)

    In contrast to the mechanical deflecting nozzle, the fluidic thrust vectoring control hires flow control methods to realize the jet vectoring, which is expected to be a more efficient way to manipulate the thrust direction. Among the main fluidic vectoring control methods, including shock vectoring control(SVC), dual throat nozzle(DTN), counter-flow(CC) and co-flow control, performance parameters such as the thrust vectoring efficiency, the thrust ratio and the discharge coefficient are compared based on published experimental and computational data. It shows that SVC, DTN and CC methods produce thrust vectoring in a wide range of Nozzle Pressure Ratio(NPR) from 1.8 to 10, and are extendable to pitch/yaw control or multi-axis control. Comparatively, DTN and co-flow control are superior to SVC and CC in the thrust loss and thrust vectoring efficiency, yet DTN is disadvantageous in the discharge coefficient as a consequence of throat injection, and the working range of co-flow method is highly limited. In pursuit of highly efficient control, some new methods of jet vectoring are introduced, and the principles, potential advantages and challenges of each method are discussed. These methods adopt after-deck-flow control and introduce little disturbance to the main jet, which are desirable for the thrust vectoring control. Such methods show promising prospects and the related experience should be drawn on for further studies.

    nozzle; thrust vectoring; control; efficiency; secondary flow

    1672-9897(2017)04-0008-08

    10.11729/syltlx20160207

    2016-12-22;

    2017-04-20

    國家自然科學(xué)基金項(xiàng)目(11572341)

    XiaoZY,JiangX,MouB,etal.Advancesinfluidicthrustvectoringtechniqueresearch.JournalofExperimentsinFluidMechanics, 2017, 31(4): 8-15. 肖中云, 江 雄, 牟 斌, 等. 流體推力矢量技術(shù)研究綜述. 實(shí)驗(yàn)流體力學(xué), 2017, 31(4): 8-15.

    V211.3

    A

    肖中云(1977-),男,四川大竹人,副研究員。研究方向:流動控制。通信地址:四川省綿陽市二環(huán)路南段6號13信箱08分箱(621000)。 E-mail:scxiaozy@sina.cn

    *通信作者 E-mail: scxiaozy@sina.cn

    猜你喜歡
    喉道射流矢量
    深海逃逸艙射流注水均壓過程仿真分析
    低壓天然氣泄漏射流擴(kuò)散特性研究
    煤氣與熱力(2022年4期)2022-05-23 12:45:00
    矢量三角形法的應(yīng)用
    基于矢量最優(yōu)估計的穩(wěn)健測向方法
    U型渠道無喉道量水槽流動規(guī)律數(shù)值模擬
    三角形法則在動態(tài)平衡問題中的應(yīng)用
    勝利油田致密砂巖油藏微觀孔隙結(jié)構(gòu)特征
    亞聲速二喉道流場不對稱現(xiàn)象研究
    射流齒形噴嘴射流流場與氣動聲學(xué)分析
    地鐵站臺活塞風(fēng)附壁射流起始段的實(shí)測和實(shí)驗(yàn)驗(yàn)證
    国产黄片视频在线免费观看| 涩涩av久久男人的天堂| 国产成人午夜福利电影在线观看| 国产高潮美女av| 七月丁香在线播放| 美女福利国产在线 | 成人美女网站在线观看视频| 精品国产一区二区三区久久久樱花 | 国产精品人妻久久久影院| av在线播放精品| 精品人妻熟女av久视频| 人妻 亚洲 视频| 大话2 男鬼变身卡| 国产精品免费大片| a 毛片基地| 一个人免费看片子| 男女免费视频国产| 夫妻午夜视频| 在线免费十八禁| 又黄又爽又刺激的免费视频.| 女性生殖器流出的白浆| 少妇人妻久久综合中文| 午夜老司机福利剧场| 熟妇人妻不卡中文字幕| 久久久久人妻精品一区果冻| 人妻少妇偷人精品九色| 午夜视频国产福利| 国产永久视频网站| 久久影院123| 亚洲精品第二区| 午夜免费观看性视频| 寂寞人妻少妇视频99o| 高清日韩中文字幕在线| 18禁裸乳无遮挡免费网站照片| 久久久久久久大尺度免费视频| 91精品一卡2卡3卡4卡| 久久精品人妻少妇| 欧美性感艳星| 国产中年淑女户外野战色| 亚洲av中文字字幕乱码综合| 国产精品偷伦视频观看了| 97在线视频观看| 黄色视频在线播放观看不卡| 亚洲av中文字字幕乱码综合| 精品人妻视频免费看| 一级二级三级毛片免费看| 午夜激情福利司机影院| 小蜜桃在线观看免费完整版高清| 国产亚洲5aaaaa淫片| av一本久久久久| 最近2019中文字幕mv第一页| 赤兔流量卡办理| 欧美精品人与动牲交sv欧美| 啦啦啦在线观看免费高清www| 大香蕉97超碰在线| 人妻 亚洲 视频| 日韩av在线免费看完整版不卡| 国产亚洲午夜精品一区二区久久| 午夜激情福利司机影院| 欧美老熟妇乱子伦牲交| 婷婷色综合大香蕉| 欧美xxxx黑人xx丫x性爽| 久久精品国产亚洲av天美| 久久久久国产网址| 色吧在线观看| 乱系列少妇在线播放| 国产精品麻豆人妻色哟哟久久| 中国美白少妇内射xxxbb| 国产精品久久久久久av不卡| 天美传媒精品一区二区| 精品午夜福利在线看| 亚洲国产欧美在线一区| 丰满人妻一区二区三区视频av| 最黄视频免费看| 免费不卡的大黄色大毛片视频在线观看| 国模一区二区三区四区视频| 日本欧美视频一区| 国产极品天堂在线| 欧美zozozo另类| 在线亚洲精品国产二区图片欧美 | 精品人妻一区二区三区麻豆| 你懂的网址亚洲精品在线观看| 国产精品国产三级国产av玫瑰| 一本久久精品| 另类亚洲欧美激情| 精品久久久久久久久亚洲| av国产精品久久久久影院| 免费黄网站久久成人精品| 国产黄片美女视频| 高清av免费在线| 国产黄色视频一区二区在线观看| 99久久精品一区二区三区| 日韩一本色道免费dvd| 日韩强制内射视频| 最近的中文字幕免费完整| 亚洲精华国产精华液的使用体验| 深夜a级毛片| 啦啦啦视频在线资源免费观看| av线在线观看网站| 免费观看av网站的网址| 丰满乱子伦码专区| 99热全是精品| 伦理电影大哥的女人| 九九久久精品国产亚洲av麻豆| 欧美日韩国产mv在线观看视频 | 国产亚洲最大av| 国产乱来视频区| 精品少妇黑人巨大在线播放| 99久久中文字幕三级久久日本| 免费观看的影片在线观看| 亚洲人成网站高清观看| 婷婷色麻豆天堂久久| 综合色丁香网| tube8黄色片| 亚洲av欧美aⅴ国产| av一本久久久久| 亚洲精品中文字幕在线视频 | 纵有疾风起免费观看全集完整版| 黄色怎么调成土黄色| 亚洲国产成人一精品久久久| 寂寞人妻少妇视频99o| 黑人猛操日本美女一级片| 欧美3d第一页| 肉色欧美久久久久久久蜜桃| 中国国产av一级| 香蕉精品网在线| 亚洲成色77777| 久久久久视频综合| 色视频www国产| 久久鲁丝午夜福利片| 菩萨蛮人人尽说江南好唐韦庄| 春色校园在线视频观看| 久久毛片免费看一区二区三区| 国产一区二区三区av在线| 久久影院123| 久久久久久久国产电影| 午夜福利高清视频| 久久 成人 亚洲| 精品亚洲成a人片在线观看 | 久久精品国产亚洲av天美| 午夜福利在线观看免费完整高清在| 欧美另类一区| 伊人久久国产一区二区| 亚洲精品自拍成人| 亚洲成人手机| 精品一区二区三卡| 看非洲黑人一级黄片| 亚洲国产最新在线播放| av专区在线播放| 久久99精品国语久久久| 男男h啪啪无遮挡| 夜夜爽夜夜爽视频| 哪个播放器可以免费观看大片| 好男人视频免费观看在线| 国产色婷婷99| 欧美bdsm另类| 在线 av 中文字幕| 中国美白少妇内射xxxbb| 九九在线视频观看精品| 亚洲内射少妇av| 啦啦啦啦在线视频资源| 久久久久久久久久人人人人人人| 99九九线精品视频在线观看视频| 观看免费一级毛片| 午夜日本视频在线| 18禁在线播放成人免费| 又黄又爽又刺激的免费视频.| 亚洲精品自拍成人| 日韩免费高清中文字幕av| 久久久精品免费免费高清| 只有这里有精品99| 国产高清国产精品国产三级 | 婷婷色综合www| 国产精品久久久久成人av| 又爽又黄a免费视频| 人妻制服诱惑在线中文字幕| 边亲边吃奶的免费视频| 免费观看a级毛片全部| 国产国拍精品亚洲av在线观看| 久久6这里有精品| 久久精品国产a三级三级三级| 成人毛片a级毛片在线播放| 精品一区在线观看国产| videos熟女内射| 男女啪啪激烈高潮av片| 日日摸夜夜添夜夜爱| 亚洲第一av免费看| 最黄视频免费看| 美女cb高潮喷水在线观看| 亚洲精华国产精华液的使用体验| 各种免费的搞黄视频| 久久久久国产精品人妻一区二区| h日本视频在线播放| 午夜激情久久久久久久| 国产爽快片一区二区三区| 永久免费av网站大全| 观看美女的网站| 老女人水多毛片| 免费高清在线观看视频在线观看| 男人和女人高潮做爰伦理| 亚洲,欧美,日韩| 日日撸夜夜添| 国产精品久久久久久久久免| 午夜福利高清视频| 精品久久久精品久久久| 91狼人影院| 18禁裸乳无遮挡动漫免费视频| 51国产日韩欧美| 欧美激情极品国产一区二区三区 | av视频免费观看在线观看| 国产精品无大码| 亚洲,一卡二卡三卡| 精品久久久噜噜| 国产一区二区三区av在线| videossex国产| 亚洲精品一二三| 高清日韩中文字幕在线| 成年免费大片在线观看| 国产精品一区二区性色av| 超碰av人人做人人爽久久| 一个人免费看片子| 日韩不卡一区二区三区视频在线| 美女高潮的动态| 亚洲欧美精品自产自拍| 校园人妻丝袜中文字幕| 成年女人在线观看亚洲视频| 亚洲美女视频黄频| 国产无遮挡羞羞视频在线观看| 欧美xxⅹ黑人| 97热精品久久久久久| 尾随美女入室| 色哟哟·www| 国产在线一区二区三区精| 男人舔奶头视频| 99久久综合免费| 日韩在线高清观看一区二区三区| 午夜激情久久久久久久| 高清毛片免费看| 国产亚洲5aaaaa淫片| 国产亚洲精品久久久com| 国产伦在线观看视频一区| 日韩三级伦理在线观看| 亚洲av在线观看美女高潮| av在线观看视频网站免费| 国产在视频线精品| 国产精品一区二区三区四区免费观看| 中国国产av一级| 99热这里只有是精品50| av免费观看日本| 国产老妇伦熟女老妇高清| 久久精品国产亚洲av涩爱| 97超视频在线观看视频| 亚洲精品乱码久久久v下载方式| 日韩av在线免费看完整版不卡| 欧美日韩国产mv在线观看视频 | 六月丁香七月| 在线精品无人区一区二区三 | 蜜桃在线观看..| 七月丁香在线播放| 直男gayav资源| 只有这里有精品99| 少妇熟女欧美另类| 男女啪啪激烈高潮av片| 久久久色成人| 亚洲最大成人中文| 久久综合国产亚洲精品| 午夜视频国产福利| 国产精品三级大全| 在线观看人妻少妇| 亚洲综合色惰| 久久精品国产a三级三级三级| 只有这里有精品99| 久久人人爽人人片av| 国产免费福利视频在线观看| 毛片一级片免费看久久久久| 黄色配什么色好看| 最近中文字幕高清免费大全6| 人妻 亚洲 视频| 18禁裸乳无遮挡动漫免费视频| 国产一级毛片在线| 亚洲av中文av极速乱| 一本一本综合久久| 肉色欧美久久久久久久蜜桃| 最后的刺客免费高清国语| 国产在线视频一区二区| 啦啦啦在线观看免费高清www| 久久久久久久久久久丰满| 日本与韩国留学比较| 91在线精品国自产拍蜜月| 激情五月婷婷亚洲| 九九久久精品国产亚洲av麻豆| 国精品久久久久久国模美| 亚洲av在线观看美女高潮| 身体一侧抽搐| 我的女老师完整版在线观看| 美女内射精品一级片tv| 国产精品蜜桃在线观看| www.av在线官网国产| 少妇人妻久久综合中文| 欧美高清性xxxxhd video| 99视频精品全部免费 在线| 国国产精品蜜臀av免费| 日韩欧美一区视频在线观看 | 91午夜精品亚洲一区二区三区| 插阴视频在线观看视频| 最近中文字幕2019免费版| 久久青草综合色| 国产精品人妻久久久久久| 国产国拍精品亚洲av在线观看| 日日啪夜夜爽| 亚洲国产精品专区欧美| 人人妻人人澡人人爽人人夜夜| 久久久久人妻精品一区果冻| 久久精品熟女亚洲av麻豆精品| 我的女老师完整版在线观看| 在线观看人妻少妇| 亚洲精品乱久久久久久| 国产成人精品一,二区| www.色视频.com| 中文字幕制服av| 国产精品福利在线免费观看| 免费看av在线观看网站| 久久国产乱子免费精品| 久热久热在线精品观看| 亚洲天堂av无毛| 久久精品久久久久久久性| 亚洲av国产av综合av卡| 亚洲国产精品成人久久小说| 国产爱豆传媒在线观看| 久久精品夜色国产| 亚洲精品自拍成人| 国产精品国产三级国产专区5o| 国产永久视频网站| 国产精品久久久久久精品电影小说 | 久久久久久久久大av| 欧美极品一区二区三区四区| 视频区图区小说| 啦啦啦啦在线视频资源| 久久99蜜桃精品久久| 乱系列少妇在线播放| 免费av不卡在线播放| 欧美国产精品一级二级三级 | 精品人妻偷拍中文字幕| 日本猛色少妇xxxxx猛交久久| 极品教师在线视频| av.在线天堂| 国产成人精品婷婷| 99热6这里只有精品| 国产熟女欧美一区二区| 2018国产大陆天天弄谢| 久久久久精品性色| 91狼人影院| 黑人高潮一二区| 老师上课跳d突然被开到最大视频| 日韩制服骚丝袜av| 丝瓜视频免费看黄片| 日韩伦理黄色片| 国产欧美日韩一区二区三区在线 | 国产黄色免费在线视频| 天堂俺去俺来也www色官网| 婷婷色麻豆天堂久久| 午夜激情福利司机影院| 亚洲精品,欧美精品| 国产视频内射| 亚洲精品日韩av片在线观看| 中文字幕免费在线视频6| 亚洲精品自拍成人| 久久久精品94久久精品| 欧美日韩视频精品一区| 校园人妻丝袜中文字幕| 久久久久久久国产电影| 欧美日韩精品成人综合77777| 美女视频免费永久观看网站| 国内揄拍国产精品人妻在线| 亚洲av在线观看美女高潮| 免费播放大片免费观看视频在线观看| 日日啪夜夜爽| 在线观看三级黄色| 精品一品国产午夜福利视频| 男女边吃奶边做爰视频| 午夜免费男女啪啪视频观看| 18禁在线播放成人免费| 免费大片黄手机在线观看| 免费观看无遮挡的男女| 久久久久久久久久人人人人人人| 国产亚洲精品久久久com| 国产免费福利视频在线观看| 狂野欧美激情性bbbbbb| 在线观看免费日韩欧美大片 | 国产精品久久久久久av不卡| 高清午夜精品一区二区三区| 亚洲精品日本国产第一区| 成人毛片a级毛片在线播放| 国产精品99久久99久久久不卡 | 午夜福利网站1000一区二区三区| 狂野欧美白嫩少妇大欣赏| 91午夜精品亚洲一区二区三区| 一级黄片播放器| 蜜臀久久99精品久久宅男| 久久精品久久久久久久性| 哪个播放器可以免费观看大片| 亚洲精品国产av蜜桃| 国产精品av视频在线免费观看| 久久精品国产鲁丝片午夜精品| 菩萨蛮人人尽说江南好唐韦庄| 国产大屁股一区二区在线视频| 爱豆传媒免费全集在线观看| a级一级毛片免费在线观看| 91久久精品国产一区二区三区| 在线播放无遮挡| 欧美性感艳星| www.av在线官网国产| 嫩草影院入口| 亚洲欧美一区二区三区黑人 | 人人妻人人看人人澡| 欧美一区二区亚洲| 免费大片黄手机在线观看| 六月丁香七月| 高清黄色对白视频在线免费看 | 日韩在线高清观看一区二区三区| 啦啦啦中文免费视频观看日本| h日本视频在线播放| 日韩免费高清中文字幕av| 高清黄色对白视频在线免费看 | 91精品一卡2卡3卡4卡| 天美传媒精品一区二区| 大陆偷拍与自拍| 日本欧美国产在线视频| 男女啪啪激烈高潮av片| 国产v大片淫在线免费观看| 一区二区av电影网| 日韩视频在线欧美| 我要看日韩黄色一级片| 五月开心婷婷网| 欧美xxxx黑人xx丫x性爽| 永久网站在线| 99热全是精品| 精品久久久久久久久亚洲| 尤物成人国产欧美一区二区三区| 久久人人爽人人爽人人片va| 色5月婷婷丁香| 国产成人免费无遮挡视频| 人体艺术视频欧美日本| 国产黄色视频一区二区在线观看| 哪个播放器可以免费观看大片| 亚洲图色成人| 亚洲内射少妇av| 尾随美女入室| 精品久久久久久久末码| 欧美精品一区二区大全| 超碰97精品在线观看| 搡女人真爽免费视频火全软件| 亚洲国产毛片av蜜桃av| 噜噜噜噜噜久久久久久91| 十八禁网站网址无遮挡 | 国产高清三级在线| av一本久久久久| 亚洲欧美日韩东京热| 最后的刺客免费高清国语| 国产精品嫩草影院av在线观看| 国产免费福利视频在线观看| 建设人人有责人人尽责人人享有的 | 偷拍熟女少妇极品色| 国产白丝娇喘喷水9色精品| 2021少妇久久久久久久久久久| 亚洲精品亚洲一区二区| 精品99又大又爽又粗少妇毛片| 国产精品国产三级国产av玫瑰| 人妻 亚洲 视频| 日本爱情动作片www.在线观看| 色网站视频免费| 亚洲欧美中文字幕日韩二区| 91aial.com中文字幕在线观看| 蜜桃久久精品国产亚洲av| 国产精品人妻久久久影院| 蜜桃在线观看..| 国产免费视频播放在线视频| 1000部很黄的大片| 国产69精品久久久久777片| 国内揄拍国产精品人妻在线| 赤兔流量卡办理| 嫩草影院新地址| 在线天堂最新版资源| 最近2019中文字幕mv第一页| 26uuu在线亚洲综合色| 观看免费一级毛片| 伦理电影免费视频| 亚洲经典国产精华液单| 国产成人aa在线观看| 亚洲精品日韩av片在线观看| 久久久久久久亚洲中文字幕| 99九九线精品视频在线观看视频| 大陆偷拍与自拍| 最近的中文字幕免费完整| 在线观看一区二区三区| 中文字幕久久专区| 在线观看美女被高潮喷水网站| 黄色日韩在线| 久久久久国产精品人妻一区二区| 成人国产麻豆网| 免费人成在线观看视频色| 综合色丁香网| 永久网站在线| 久热久热在线精品观看| 亚洲国产欧美在线一区| 十八禁网站网址无遮挡 | 99热这里只有是精品50| 免费大片18禁| 欧美三级亚洲精品| 青春草国产在线视频| www.色视频.com| 另类亚洲欧美激情| 一区二区三区免费毛片| 亚洲av国产av综合av卡| 国产乱人偷精品视频| 成人特级av手机在线观看| www.av在线官网国产| 国产片特级美女逼逼视频| 欧美日韩视频精品一区| 成人美女网站在线观看视频| 天美传媒精品一区二区| 亚洲av日韩在线播放| 国产精品伦人一区二区| av又黄又爽大尺度在线免费看| 亚洲精品一二三| 丰满少妇做爰视频| 婷婷色麻豆天堂久久| 国产在线一区二区三区精| 欧美少妇被猛烈插入视频| 成人免费观看视频高清| 51国产日韩欧美| 亚洲怡红院男人天堂| 亚洲人成网站在线观看播放| 王馨瑶露胸无遮挡在线观看| 男人舔奶头视频| 伊人久久国产一区二区| 成人18禁高潮啪啪吃奶动态图 | 欧美一级a爱片免费观看看| 美女xxoo啪啪120秒动态图| 五月开心婷婷网| 丰满少妇做爰视频| av在线蜜桃| 综合色丁香网| 日韩一区二区三区影片| freevideosex欧美| 国产精品蜜桃在线观看| 亚洲国产av新网站| 日日啪夜夜爽| 亚洲av福利一区| 午夜免费观看性视频| 搡老乐熟女国产| 亚洲久久久国产精品| 久久韩国三级中文字幕| 色网站视频免费| 女性被躁到高潮视频| av.在线天堂| 最近手机中文字幕大全| 性高湖久久久久久久久免费观看| 在线观看免费高清a一片| 波野结衣二区三区在线| 天天躁日日操中文字幕| 全区人妻精品视频| 国产亚洲一区二区精品| xxx大片免费视频| 国产探花极品一区二区| 亚洲av免费高清在线观看| 久久99精品国语久久久| 99热这里只有是精品50| av在线播放精品| 两个人的视频大全免费| 久久久久久久国产电影| 少妇的逼好多水| 久久女婷五月综合色啪小说| 免费黄色在线免费观看| 国内揄拍国产精品人妻在线| 久久久欧美国产精品| 一本久久精品| 亚洲成人中文字幕在线播放| 精品一品国产午夜福利视频| 好男人视频免费观看在线| 国产精品偷伦视频观看了| 欧美一级a爱片免费观看看| 亚洲aⅴ乱码一区二区在线播放| 日韩av在线免费看完整版不卡| 国产日韩欧美亚洲二区| 亚洲精品久久午夜乱码| 国产女主播在线喷水免费视频网站| 在线亚洲精品国产二区图片欧美 | 国产欧美亚洲国产| 国产淫片久久久久久久久| 成人国产麻豆网| 在线观看免费视频网站a站| 视频中文字幕在线观看| av福利片在线观看| 国产精品国产三级国产专区5o| 国产精品人妻久久久影院| 色综合色国产| 青春草国产在线视频| 国产亚洲91精品色在线| 网址你懂的国产日韩在线| 韩国高清视频一区二区三区| 国产亚洲欧美精品永久| 黑丝袜美女国产一区| 男人狂女人下面高潮的视频| 多毛熟女@视频| 夜夜看夜夜爽夜夜摸| 久久久久久久国产电影| 亚洲国产精品一区三区| 亚洲欧美精品自产自拍| 亚洲不卡免费看| 少妇裸体淫交视频免费看高清| 国产免费一级a男人的天堂| 九九爱精品视频在线观看| 久久综合国产亚洲精品| 22中文网久久字幕|