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

    飛機(jī)尾噴流誘導(dǎo)速度建模與仿真

    2021-10-22 03:34:44楊立波
    紅外技術(shù) 2021年10期
    關(guān)鍵詞:尾渦噴流噴口

    楊立波,叢 巖

    飛機(jī)尾噴流誘導(dǎo)速度建模與仿真

    楊立波1,叢 巖2

    (1. 廣東科技學(xué)院,廣東 東莞 523083;2. 空軍航空大學(xué),吉林 長(zhǎng)春 130000)

    本文對(duì)飛機(jī)尾噴流誘導(dǎo)速度進(jìn)行了仿真研究。建立了飛機(jī)尾渦模型以及尾噴口噴流模型,并仿真了飛機(jī)流場(chǎng)中的誘導(dǎo)速度。然后應(yīng)用CFD(computational fluid dynamics)對(duì)飛機(jī)的流場(chǎng)進(jìn)行計(jì)算,并將CFD計(jì)算結(jié)果與尾噴流模型計(jì)算結(jié)果進(jìn)行對(duì)比。仿真對(duì)比結(jié)果表明:飛機(jī)尾后100 m區(qū)域內(nèi),尾渦模型計(jì)算結(jié)果誤差較大,此時(shí)應(yīng)采用CFD對(duì)尾流場(chǎng)進(jìn)行計(jì)算;而在尾后100 m區(qū)域外尾渦模型與CFD計(jì)算結(jié)果較為一致,尾渦模型的計(jì)算結(jié)果能夠達(dá)到精度要求。

    尾噴流;誘導(dǎo)速度;尾渦模型;噴流模型;CFD

    0 引言

    由于空氣的粘性作用,飛機(jī)在飛行過(guò)程中,其機(jī)身和尾后會(huì)產(chǎn)生一系列的脫體渦,主要分為前緣渦、后緣渦以及翼尖渦。前緣渦和后緣渦相比機(jī)翼翼尖產(chǎn)生的翼尖渦強(qiáng)度要小得多,且不穩(wěn)定,脫體后很快就會(huì)削弱。而翼尖渦可以在機(jī)翼后方200倍翼展內(nèi)穩(wěn)定,渦核可以持續(xù)數(shù)分鐘不消散。同時(shí),飛行過(guò)程中發(fā)動(dòng)機(jī)向后噴射出大量高速氣流也會(huì)對(duì)飛機(jī)的尾流場(chǎng)產(chǎn)生一定的影響。由于尾渦和發(fā)動(dòng)機(jī)噴流的影響,對(duì)周圍大氣會(huì)產(chǎn)生變化率大、分布不規(guī)律的誘導(dǎo)速度。因此,飛機(jī)的尾噴流誘導(dǎo)速度可以分為尾渦誘導(dǎo)速度和噴流誘導(dǎo)速度,其中尾渦誘導(dǎo)速度主要由翼尖渦引起,而噴流誘導(dǎo)速度主要指尾噴口噴射的高速氣流產(chǎn)生的誘導(dǎo)速度。

    飛機(jī)尾噴流誘導(dǎo)速度可以對(duì)經(jīng)由尾流場(chǎng)區(qū)域內(nèi)的飛機(jī)產(chǎn)生強(qiáng)烈的震顫、側(cè)翻下沉、甚至是空中停車,嚴(yán)重威脅飛行安全[1-5]。同時(shí),誘導(dǎo)速度還會(huì)影響紅外誘餌彈、箔條彈的運(yùn)動(dòng)擴(kuò)散規(guī)律,影響干擾效果。誘導(dǎo)速度甚至?xí)?dǎo)致加油機(jī)加油椎管在空中來(lái)回?cái)[動(dòng),導(dǎo)致加油失敗。因此,對(duì)飛機(jī)尾噴流誘導(dǎo)速度進(jìn)行研究意義十分重大。

    國(guó)內(nèi)外的學(xué)者對(duì)飛機(jī)尾噴流誘導(dǎo)速度進(jìn)行了大量研究。Loucel[6]建立了尾渦模型,并研究了加快尾渦消亡及避免飛行事故的方法。程學(xué)東[7]對(duì)飛機(jī)的尾渦進(jìn)行了研究,并分析了飛機(jī)編隊(duì)飛行中尾流的規(guī)避方法。魏志強(qiáng)[8]對(duì)尾渦的消散和渦核的運(yùn)動(dòng)進(jìn)行了建模研究,并將計(jì)算結(jié)果與激光雷達(dá)實(shí)測(cè)結(jié)果對(duì)比,驗(yàn)證了模型的可信性。黃爍橋[9]提出了利用發(fā)動(dòng)機(jī)噴流加快尾渦消亡的觀點(diǎn),并對(duì)最終的結(jié)果進(jìn)行了仿真分析。

    由于前緣渦、后緣渦脫體后很快消散,因此文中主要求解翼尖渦及噴流產(chǎn)生的誘導(dǎo)速度。

    1 尾噴流誘導(dǎo)速度模型

    本章主要建立尾渦模型和噴流模型,分別求解翼尖渦和尾噴口噴流產(chǎn)生的誘導(dǎo)速度。

    1.1 飛機(jī)尾渦模型

    飛機(jī)在飛行過(guò)程中,機(jī)翼上下表面會(huì)產(chǎn)生一個(gè)壓力差。而這個(gè)壓力差在機(jī)翼翼尖兩側(cè)匯合,從而形成了一對(duì)向內(nèi)卷曲的翼尖渦流[10],如圖1所示。翼尖渦流渦核強(qiáng)度較大,且相對(duì)穩(wěn)定,因此其產(chǎn)生的誘導(dǎo)速度較大,可以對(duì)飛機(jī)尾后處在該區(qū)域內(nèi)的其他飛機(jī)產(chǎn)生強(qiáng)烈的影響,甚至造成飛行事故,如圖2所示。

    圖1 飛機(jī)翼尖渦

    1.1.1 尾渦的初始強(qiáng)度

    飛機(jī)的升力主要由附著在機(jī)翼剖面上的渦環(huán)量產(chǎn)生[11],設(shè)()為機(jī)翼上的渦環(huán)量,則機(jī)翼產(chǎn)生的升力p可以表示為:

    式中:為大氣密度;為飛機(jī)翼展;p為飛機(jī)的飛行速度。

    對(duì)橢圓形翼和后掠翼,()可以表示為:

    式中:為渦核的初始強(qiáng)度。將公式(2)代入公式(1),則飛機(jī)的升力可以表示為:

    設(shè)飛機(jī)在飛行過(guò)程中所承受的法向過(guò)載為n、為飛機(jī)質(zhì)量,則飛機(jī)的升力p=nmg,代入公式(3)可以得到渦核的初始強(qiáng)度為[12-13]:

    =4nmg/pVp(4)

    設(shè)0為機(jī)翼兩側(cè)翼尾核的初始間距,其值為:

    0=(5)

    式中:為翼形參數(shù),值的大小與機(jī)翼升力的橫向分布有關(guān)。對(duì)于橢圓形機(jī)翼=p/4,對(duì)于后掠翼取值為0.75~0.8。

    由公式(4)可知,尾渦的初始強(qiáng)度主要取決于飛機(jī)質(zhì)量、飛機(jī)法向過(guò)載、飛行速度、大氣密度、飛機(jī)翼展以及翼形參數(shù)。

    則飛機(jī)平飛時(shí)尾渦的初始渦強(qiáng)度為[12-13]:

    設(shè)c尾渦渦核半徑,一般c表示為[10]:

    c=0.050(7)

    1.1.2 尾渦速度模型

    設(shè)渦線微元為d,渦核強(qiáng)度為。則其對(duì)渦線外任意一點(diǎn)產(chǎn)生的切線速度為:

    式中:m為點(diǎn)到渦線的垂直距離;t為點(diǎn)與渦線微元之間的夾角。

    對(duì)公式(8)積分,則渦線在點(diǎn)處產(chǎn)生的切線速度w為:

    式中:t和t為點(diǎn)與渦線兩端之間的夾角。當(dāng)渦線無(wú)限長(zhǎng)時(shí)t=t=0°,則w為:

    由于尾渦渦核初始強(qiáng)度很大,在渦核半徑c內(nèi)部可以認(rèn)為旋轉(zhuǎn)角速度恒定。令公式(11)中的m=c,則渦核內(nèi)部的旋轉(zhuǎn)角速度c為:

    式中:w,c為渦核半徑上的切線速度。

    則當(dāng)點(diǎn)位于渦核內(nèi)部,則切線速度w為:

    Hallock-burnham模型無(wú)論對(duì)渦核內(nèi)還是渦核外的點(diǎn)都適用[14-15],因此文中采用該模型對(duì)尾渦切線速度進(jìn)行計(jì)算,即:

    由于飛機(jī)的尾渦為渦線與飛機(jī)縱軸平行且強(qiáng)度相等的兩個(gè)反向自由渦,而自由渦只存在切線的誘導(dǎo)速度無(wú)徑向速度,因此尾渦對(duì)飛機(jī)尾后流場(chǎng)中任意點(diǎn)的誘導(dǎo)速度都只在軸和軸上存在一定的分量,在軸上無(wú)分量。設(shè)尾渦的兩個(gè)點(diǎn)渦核心坐標(biāo)為(1,1)和(2,2),由公式(14)可知,任意點(diǎn)在翼尖渦中的誘導(dǎo)速度為[8]:

    VV1+V2(20)

    VV1+V2(21)

    式中:(,,)為點(diǎn)坐標(biāo);1和2為尾渦渦核強(qiáng)度,且1=-2[16]。

    1.1.3 尾渦飄降模型

    由于飛機(jī)尾渦左右對(duì)稱且反向旋轉(zhuǎn),在互相誘導(dǎo)的作用下,其渦核高度不斷地下降[17-18]。飛機(jī)尾渦的下降速度可以達(dá)到每分鐘幾百英尺,但當(dāng)下降到一定高度后,尾渦就停止下降,直至消散。

    當(dāng)渦核高度高于地面50m且無(wú)風(fēng)時(shí),尾渦的下降速度近似為定值。下降速度v為:

    則時(shí)刻尾渦渦線的坐標(biāo)為(0-vt,0/2)和(0-vt, -0/2)。

    1.2 尾噴口噴流模型

    飛機(jī)在飛行過(guò)程中,發(fā)動(dòng)機(jī)噴口在飛機(jī)尾后噴射出大量的高速、高溫氣體,因此尾噴口噴流能夠在飛機(jī)尾后空間內(nèi)產(chǎn)生較大的誘導(dǎo)速度。噴流產(chǎn)生的誘導(dǎo)速度雖大,但是氣流在噴射出后迅速減速,相比尾渦,噴流的影響區(qū)域并不大。

    將發(fā)動(dòng)機(jī)的噴流場(chǎng)沿軸向分為穩(wěn)定區(qū)和不穩(wěn)定區(qū),沿徑向分為穩(wěn)速核心區(qū)、減速區(qū)、混流增速區(qū)以及邊界層區(qū),如圖4所示。

    圖4 噴流模型

    設(shè)m為噴流場(chǎng)中心速度,由于初始階段處于穩(wěn)速核心區(qū)內(nèi),因此其m=n0不變。在減速區(qū)內(nèi),中心速度m可以表示為:

    式中:0為發(fā)動(dòng)機(jī)噴口直徑;1為流場(chǎng)軸心線上的點(diǎn)到噴口之間的距離;0為圖4中的轉(zhuǎn)折截面到噴口的距離,其值為:

    式中:n為噴流馬赫數(shù)。聯(lián)立公式(23)和公式(24)則可以得到噴流場(chǎng)軸線上的速度。

    噴流流場(chǎng)徑向速度可以表示為:

    式中:為計(jì)算點(diǎn)到軸線的距離;tp為噴流邊界層到軸線之間的距離。在不穩(wěn)定區(qū)域內(nèi)tp大概為噴口半徑的3.5倍,tp的具體值可以由下式求得,即:

    在穩(wěn)定區(qū)域內(nèi),其值為:

    式中:1為射流擴(kuò)散角。

    2 飛機(jī)流場(chǎng)CFD計(jì)算

    本章對(duì)飛機(jī)的流場(chǎng)進(jìn)行CFD計(jì)算,其中尾流CFD計(jì)算主要計(jì)算尾渦誘導(dǎo)速度,而噴流CFD計(jì)算則主要計(jì)算飛機(jī)噴流誘導(dǎo)速度。

    2.1 飛機(jī)尾流CFD計(jì)算

    由于飛機(jī)尾渦可以在尾后很長(zhǎng)區(qū)域內(nèi)保持穩(wěn)定,因此,本文將計(jì)算域設(shè)置為520m×75m×48m的長(zhǎng)方體,以滿足尾渦誘導(dǎo)速度的計(jì)算要求。如圖5所示,飛機(jī)機(jī)頭距計(jì)算域前端100m,飛機(jī)質(zhì)心位于計(jì)算域原點(diǎn),來(lái)流方向沿軸負(fù)方向與機(jī)身平行。邊界層網(wǎng)格第一層厚度設(shè)置為10-2m,采用四面體非結(jié)構(gòu)網(wǎng)格對(duì)計(jì)算域進(jìn)行剖分,網(wǎng)格數(shù)總計(jì)為5137.6萬(wàn)。初始條件設(shè)置為高度=0 km、來(lái)流馬赫數(shù)=0.8、壓力=101325Pa,大氣溫度=288.15K。尾噴口設(shè)置為質(zhì)量入口,質(zhì)量流量設(shè)置為111kg/s,總溫為774.17K,計(jì)算域邊界設(shè)置為壓強(qiáng)遠(yuǎn)場(chǎng),飛機(jī)壁面設(shè)置為無(wú)滑移絕熱壁面。對(duì)機(jī)翼后緣、機(jī)翼下表面以及飛機(jī)尾后進(jìn)行網(wǎng)格加密處理,重點(diǎn)計(jì)算飛機(jī)機(jī)翼前緣渦、后緣渦以及翼尖渦的流場(chǎng)。

    圖6為=0.8時(shí),CFD計(jì)算得到的飛機(jī)尾后渦線圖。

    由圖6的計(jì)算結(jié)果可知,飛機(jī)在0°迎角下其前緣渦、后緣渦以及翼尖渦較為分散,渦核強(qiáng)度不大。在5°迎角下,前緣渦、后緣渦以及翼尖渦匯合在了一起,在飛機(jī)尾后形成了一個(gè)十分強(qiáng)烈的渦。飛機(jī)在5°迎角飛行時(shí),機(jī)翼產(chǎn)生的升力增大,機(jī)翼上下表面壓力差也隨之增大,最終導(dǎo)致翼尖處形成的誘導(dǎo)渦強(qiáng)度增強(qiáng)。

    2.2 飛機(jī)噴流CFD計(jì)算

    控制方程為Navier-Stokes方程,采用分離隱式求解器,壓力速度耦合使用SIMPLEC算法,各物理量的離散采用二階迎風(fēng)格式,迭代收斂準(zhǔn)則為殘差小于1×10-4。噴管壁面為無(wú)滑移絕熱壁面,尾噴口為質(zhì)量入口,外流場(chǎng)前端為壓強(qiáng)入口,徑向四周邊界設(shè)為壓強(qiáng)遠(yuǎn)場(chǎng),外流場(chǎng)后端設(shè)為壓強(qiáng)出口。計(jì)算模型選用標(biāo)準(zhǔn)-湍流模型。

    圖5 計(jì)算域示意圖

    圖6 飛機(jī)尾后渦線圖

    設(shè)噴口半徑為0.4m,計(jì)算域?yàn)檫呴L(zhǎng)16m×16m×500 m的長(zhǎng)方體。計(jì)算域的前端距噴口1m,噴口軸心位于來(lái)流的中心,如圖7所示。采用六面體結(jié)構(gòu)網(wǎng)格對(duì)流場(chǎng)進(jìn)行計(jì)算,對(duì)噴口附近和噴口軸心線區(qū)域進(jìn)行網(wǎng)格加密,網(wǎng)格數(shù)總計(jì)117.7萬(wàn)。初始條件設(shè)置為高度=0km、來(lái)流馬赫數(shù)=0、壓力=101325Pa,大氣溫度=288.15K。發(fā)動(dòng)機(jī)處于非加力最大狀態(tài),噴管進(jìn)口質(zhì)量流量111kg/s,總溫774.17K。采用CFD專業(yè)求解器FLUENT進(jìn)行流場(chǎng)計(jì)算。計(jì)算得到發(fā)動(dòng)機(jī)尾流場(chǎng)速度分布如圖8所示。

    由圖8的計(jì)算結(jié)果可知,在發(fā)動(dòng)機(jī)噴流的影響下,噴口后的流場(chǎng)速度迅速增加,但影響的區(qū)域僅局限于尾噴口正后方的狹窄區(qū)域,且噴流在尾噴口噴出后速度迅速減小,發(fā)動(dòng)機(jī)噴流的誘導(dǎo)速度在尾噴口后50m內(nèi)逐漸降為0。

    3 模型驗(yàn)證與結(jié)果分析

    應(yīng)用CFD流場(chǎng)計(jì)算的方法求解飛機(jī)尾噴流誘導(dǎo)速度是目前公認(rèn)可信度最高、最為精確的數(shù)值計(jì)算方法,但是CFD流場(chǎng)計(jì)算速度較慢,且由于尾渦可以在飛機(jī)尾后很長(zhǎng)區(qū)域內(nèi)保持穩(wěn)定,因此計(jì)算域的設(shè)置需要足夠長(zhǎng)。同時(shí)每改變初始條件都需要對(duì)結(jié)果進(jìn)行重新計(jì)算。而采用飛機(jī)尾噴流模型計(jì)算誘導(dǎo)速度的方法計(jì)算速度較快,但是結(jié)果不準(zhǔn)確。本節(jié)將尾噴流誘導(dǎo)速度模型的計(jì)算結(jié)果與CFD計(jì)算結(jié)果進(jìn)行對(duì)比,從而研究能夠同時(shí)兼顧計(jì)算精度和計(jì)算效率的誘導(dǎo)速度計(jì)算方法。

    3.1 尾流場(chǎng)誘導(dǎo)速度結(jié)果對(duì)比

    設(shè)飛機(jī)翼展=15m,飛機(jī)質(zhì)量=27273kg,飛行馬赫數(shù)=0.8,0=0.75,飛行高度=0 km。則由2.1節(jié)建立的尾渦模型可以計(jì)算尾流場(chǎng)中任意點(diǎn)的誘導(dǎo)速度。

    CFD流場(chǎng)計(jì)算中,初始條件的設(shè)置與2.1節(jié)一致。通過(guò)流場(chǎng)計(jì)算結(jié)果可知,當(dāng)迎角=0.5°時(shí),飛機(jī)的升力系數(shù)為0.10,此時(shí)飛機(jī)的升力與重力相等,飛機(jī)處于平飛狀態(tài)。圖9~圖12為飛機(jī)尾后不同距離時(shí),尾渦模型與CFD計(jì)算得到的尾流場(chǎng)誘導(dǎo)速度對(duì)比圖。

    圖7 尾噴口計(jì)算域示意圖

    圖8 尾噴流速度分布圖

    由圖9和圖11的對(duì)比可知,飛機(jī)尾后10 m處CFD計(jì)算的誘導(dǎo)速度分布規(guī)律較為雜亂無(wú)章,與尾渦模型的計(jì)算結(jié)果差距較大。這是由于流經(jīng)飛機(jī)的氣流受到機(jī)身、機(jī)腹、進(jìn)氣道外壁等的干擾,產(chǎn)生的誘導(dǎo)速度極為不規(guī)律,同時(shí)飛機(jī)尾后的氣流不但受到機(jī)翼前緣渦、后緣渦、翼尖渦的影響,垂尾以及平尾產(chǎn)生的誘導(dǎo)渦都會(huì)對(duì)流場(chǎng)內(nèi)的誘導(dǎo)速度產(chǎn)生干擾。同時(shí)由于飛機(jī)飛行姿態(tài)、飛行速度的差異,還存在不確定的渦破碎現(xiàn)象。因此,飛機(jī)尾后的誘導(dǎo)速度及其復(fù)雜、不規(guī)律,與CFD計(jì)算結(jié)果差異較大。但是上述干擾的影響較不穩(wěn)定,只停留在飛機(jī)尾后的一小段區(qū)域內(nèi)。

    圖9 飛機(jī)尾后10 m處尾流場(chǎng)誘導(dǎo)速度

    圖10 飛機(jī)尾后100 m處尾流場(chǎng)誘導(dǎo)速度

    圖11 飛機(jī)尾后10 m對(duì)比結(jié)果

    圖12 飛機(jī)尾后100 m對(duì)比結(jié)果

    由圖10和圖12的對(duì)比可知,飛機(jī)尾后100m處,尾渦模型與CFD的計(jì)算結(jié)果差距進(jìn)一步減小,模型的計(jì)算結(jié)果能夠真實(shí)地反應(yīng)流場(chǎng)中誘導(dǎo)速度的變化規(guī)律。因?yàn)榇藭r(shí),飛機(jī)尾流場(chǎng)距離飛機(jī)較遠(yuǎn),飛機(jī)機(jī)翼產(chǎn)生的前緣渦、后緣渦開始逐步消散,渦核強(qiáng)度大大降低,而此時(shí)機(jī)身對(duì)氣流的干擾作用也被消減的很弱。由于翼尖渦十分穩(wěn)定,在飛機(jī)尾后很長(zhǎng)的距離都不會(huì)發(fā)生消散,流場(chǎng)中的誘導(dǎo)速度主要由翼尖渦決定。因此,尾渦模型與CFD結(jié)果較為接近。

    圖9~圖12中的尾流場(chǎng)誘導(dǎo)速度對(duì)比結(jié)果可知,由于機(jī)身、機(jī)腹、進(jìn)氣道外壁等的干擾的存在,如果不依靠CFD,單純采用理論建模的方法計(jì)算尾流場(chǎng)誘導(dǎo)速度,必須要對(duì)以上各種因素深入分析,分別建模最終確定結(jié)果,實(shí)現(xiàn)較為困難,可行性不高,并且效果甚微。因此,單純采用建模的方法來(lái)計(jì)算尾流場(chǎng)誘導(dǎo)速度顯然不能滿足精度的要求。文中采用CFD流場(chǎng)計(jì)算和尾流模型相結(jié)合的方法來(lái)計(jì)算誘導(dǎo)速度。對(duì)于飛機(jī)尾后距離尾噴口100m之內(nèi)的區(qū)域,采用CFD流場(chǎng)計(jì)算的方法計(jì)算誘導(dǎo)速度,對(duì)于距離尾噴口大于100m的區(qū)域,采用尾流模型計(jì)算誘導(dǎo)速度。這樣大大的縮短了飛機(jī)CFD流場(chǎng)計(jì)算中的計(jì)算域長(zhǎng)度,減少網(wǎng)格數(shù)量,在保證計(jì)算精度的前提下加快了計(jì)算速度。

    3.2 噴流場(chǎng)誘導(dǎo)速度結(jié)果對(duì)比

    設(shè)發(fā)動(dòng)機(jī)噴口直徑0=0.8m,噴管出口速度n0=867.5,尾噴流馬赫數(shù)n=1.78,來(lái)流馬赫數(shù)=0.8,其余條件的設(shè)置與2.2節(jié)一致。則發(fā)動(dòng)機(jī)噴流速度CFD計(jì)算結(jié)果和噴流模型計(jì)算結(jié)果如圖13和圖14所示。

    圖13 發(fā)動(dòng)機(jī)噴流速度CFD計(jì)算結(jié)果

    圖14 噴流模型計(jì)算結(jié)果

    由圖13和圖14的計(jì)算結(jié)果可知,噴流場(chǎng)中心軸線速度最大,距中心軸越遠(yuǎn)速度越小,噴流對(duì)于邊界層外的流場(chǎng)區(qū)域幾乎沒(méi)有影響。噴流場(chǎng)核心區(qū)速度大,但核心區(qū)空間體積較小,減速區(qū)噴流速度迅速下降。噴流模型計(jì)算的誘導(dǎo)速度分布規(guī)律與CFD計(jì)算結(jié)果相一致,因此采用噴流模型計(jì)算噴流速度能夠正確地反應(yīng)噴流場(chǎng)的速度分布規(guī)律,不需要對(duì)模型進(jìn)行二次校正。

    4 結(jié)論

    本文對(duì)飛機(jī)尾噴流誘導(dǎo)速度進(jìn)行了仿真分析。通過(guò)建立尾渦模型和尾噴口噴流模型研究了尾噴流誘導(dǎo)速度的數(shù)值計(jì)算方法。然后,應(yīng)用CFD對(duì)飛機(jī)的尾噴流流場(chǎng)進(jìn)行了計(jì)算,并與數(shù)學(xué)模型的計(jì)算結(jié)果進(jìn)行對(duì)比。結(jié)果表明,在飛機(jī)尾后100m區(qū)域內(nèi),CFD計(jì)算的誘導(dǎo)速度與尾渦模型差距較大,而噴流模型與CFD的計(jì)算結(jié)果吻合得較好。為了兼顧計(jì)算精度和計(jì)算效率,文中認(rèn)為在尾后100m區(qū)域內(nèi)應(yīng)當(dāng)采用CFD計(jì)算結(jié)果,而在100m區(qū)域外可以采用尾渦模型計(jì)算流場(chǎng)誘導(dǎo)速度。

    [1] Proctor F H, Hamilton D W, Switzer G F. TASS driven algorithms for wake prediction[C]//41, 2006: 1-20.

    [2] Hinton A, O Connor C J. Development of a wake vortex spacing system for airport capacity enhancement and delay reduction[C]//19, 2000, 1: 3E6/1-3E610.

    [3] Reimer H M, Vicroy D D. A preliminary study of a wake vortex encounter hazard boundary for a B737-100 airplane[R]. NASA-96-TM110223, 1996: 1-18.

    [4] Crichley J B , Foot P B. Analysis of incidents reported between 1972 and 1990 [C]//, 1992: 1-10.

    [5] Murphy B, Callaghan O, Fox M. Overview of the structures investigation for the American airline flight 587 investigation[C]//46, 2005: 1-9.

    [6] Loucel R E, Crouch J D. Flight- simulator study of airplane encounters with perturbed trailing vortices[J]., 2005, 42(4): 924-931.

    [7] 程學(xué)東, 范修宏. 大編隊(duì)飛行尾流規(guī)避研究[J]. 保定學(xué)院學(xué)報(bào), 2012, 23(3): 104-107.

    CHENG Xuedong, FAN Xiuhong. Research on avoiding the wake of the larger formatting[J]., 2012, 23(3):104-107.

    [8] 魏志強(qiáng), 徐肖豪. 飛機(jī)尾渦流場(chǎng)的建模與仿真計(jì)算研究[J]. 交通運(yùn)輸系統(tǒng)工程與信息, 2010, 10(4): 186-191.

    WEI Zhiqiang, XU Xiaohao. Modeling and simulating of flow field for aircraft wake vortex[J]., 2010, 10(4): 186-191.

    [9] 黃爍橋, 申功炘, Robert Konrat, 等. 噴流對(duì)飛機(jī)尾流渦影響的試驗(yàn)研究[J]. 航空學(xué)報(bào), 2010, 31(5): 899-908.

    HUANG Shuoqiao, SHEN Gongxin, ROBERT Konrat, et al. Experimental investigation of influence of jets on aircraft wake vortices[J]., 2010, 31(5): 899-908.

    [10] 周彬, 王雪松, 王濤, 等. 側(cè)向風(fēng)速對(duì)飛機(jī)尾流運(yùn)動(dòng)的影響[J]. 航空學(xué)報(bào), 2009, 30(5): 773-779.

    ZHOU Bin, WANG Xuesong, WANG Tao, et al. Influence of crosswind speeds on aircraft wake vortex movement[J]., 2009, 30(5): 773-779.

    [11] 錢翼稷. 空氣動(dòng)力學(xué)[M]. 北京: 北京航空航天大學(xué)出版社, 2005.

    QIAN Yishe.[M]. Beijing: Beijing University of Aeronautics and Astronautics Press, 2005.

    [12] Laurence H M. ATM decision support tool for wake vortex hazard management combining sensors and modeling[C]//6, 2014: doi:10.2514/6.2014-2332.

    [13] Schwarz C, Hahn K U, Fischenberg D. Wake encounter severity assessment based on validated aerodynamic interaction models[C]//,,, 2010: 1-9.

    [14] 周彬, 王雪松, 王濤, 等. 飛機(jī)尾流的介電常數(shù)分布特性分析[J]. 微波學(xué)報(bào), 2008, 24: 24-32.

    ZHOU Bin, WANG Xuesong, WANG Tao, et al. Analysis of the dielectric constant distributing characteristic of aircraft wake vortices[J]., 2008, 24:24-32.

    [15] 李大偉, 王宏倫. 自動(dòng)空中加油階段加油機(jī)尾渦流場(chǎng)建模與仿真[J]. 北京航空航天大學(xué)學(xué)報(bào), 2010, 36(7): 776-797.

    LI Dawei, WANG Honglun. Wake vortex effectmodeling and simulation in automated aerial refueling[J]., 2010, 36(7): 776-797.

    [16] Loucel R E, Crouch J D. Flight-simulator study of airplane encounters with perturbed trailing vortices[J]., 2005, 42(4): 924-931.

    [17] Proctor F H. Numerical simulation of wake vortices measured during the Idaho falls and Memphis field programs[C]//14,1996: 1-18.

    [18] 周彬. 飛機(jī)尾流的微結(jié)構(gòu)特征及散射特性研究[D]. 長(zhǎng)沙: 國(guó)防科學(xué)技術(shù)大學(xué), 2009.

    ZHOU Bin. Study on the microstructure and scattering characteristics of aircraft wake cortices[D]. Changsha: National University of Defense Technology, 2009.

    Modeling and Simulation of Jet-and Wake-Flow-Induced Velocity of Aircraft

    YANG Libo1,CONG Yan2

    (1.,523083,;2.,130013,)

    The jet-and wake-flow-induced velocities of aircraft were analyzed in this study. Wake vortex and jet flow models were established, and the induced velocity of the aircraft was simultaneously simulated in the wake flow field. Then, the flow field of the aircraft was computed via computational fluid dynamics(CFD), the results of which were compared with those of the wake vortex models. The comparison showed that the errors computed using the wake vortex models were large, which suggests that the wake flow field should be computed via CFD within 100m of the aircraft tail; these results are consistent with the CFD results beyond 100m from the aircraft tail, and the wake vortex models meet the accuracy requirements.

    jet and wake flow, induced velocity, wake vortex models, jet flow models, CFD

    TJ7

    A

    1001-8891(2021)10-0940-09

    2020-03-19;

    2020-04-28.

    楊立波(1981-),男,黑龍江木蘭縣人,漢族,碩士,副教授,研究方向?yàn)橹悄芸刂萍夹g(shù)。E-mail:yanglibo_1981ylb@163.com。

    東莞市社會(huì)科技發(fā)展(一般)項(xiàng)目(2019507154529)。

    猜你喜歡
    尾渦噴流噴口
    不同B-V頻率下的飛機(jī)尾渦數(shù)值模擬研究
    噴口形狀對(duì)噴水推進(jìn)器性能的影響
    飛參數(shù)據(jù)在某型飛機(jī)噴口工作監(jiān)控中的應(yīng)用
    高空巡航階段的飛機(jī)尾渦流場(chǎng)演化特性研究
    “慧眼”發(fā)現(xiàn)迄今距離黑洞最近的高速噴流
    基于激光雷達(dá)回波的動(dòng)態(tài)尾渦特征參數(shù)計(jì)算
    干擾板作用下飛機(jī)尾渦流場(chǎng)近地演變機(jī)理研究
    小蜘蛛
    高速?zèng)_床噴口壓力及位置對(duì)工件出模影響分析
    噴流干擾氣動(dòng)熱數(shù)值模擬的若干影響因素
    亚洲精品美女久久久久99蜜臀| 最近中文字幕2019免费版| 国产成人精品无人区| 最黄视频免费看| 欧美国产精品一级二级三级| 欧美国产精品va在线观看不卡| 悠悠久久av| 欧美变态另类bdsm刘玥| 最近最新中文字幕大全免费视频| 最黄视频免费看| 高潮久久久久久久久久久不卡| 下体分泌物呈黄色| 天天躁狠狠躁夜夜躁狠狠躁| 国产免费现黄频在线看| 如日韩欧美国产精品一区二区三区| 国产精品 欧美亚洲| 欧美 亚洲 国产 日韩一| 国产在线一区二区三区精| 成年人免费黄色播放视频| 欧美亚洲 丝袜 人妻 在线| 美女扒开内裤让男人捅视频| 欧美人与性动交α欧美软件| 国产精品影院久久| 19禁男女啪啪无遮挡网站| 中文欧美无线码| av天堂在线播放| 在线天堂中文资源库| 两性夫妻黄色片| 久久ye,这里只有精品| 亚洲国产精品成人久久小说| 精品亚洲成国产av| 亚洲人成77777在线视频| 99热国产这里只有精品6| 午夜免费鲁丝| 一区二区三区激情视频| h视频一区二区三区| 亚洲七黄色美女视频| e午夜精品久久久久久久| 少妇猛男粗大的猛烈进出视频| 每晚都被弄得嗷嗷叫到高潮| 国产精品二区激情视频| 成人手机av| 男女国产视频网站| 九色亚洲精品在线播放| 国产xxxxx性猛交| 国产精品av久久久久免费| 精品久久久久久电影网| 男男h啪啪无遮挡| 久久久久久人人人人人| 国产欧美亚洲国产| 美女主播在线视频| 啦啦啦在线免费观看视频4| 免费在线观看影片大全网站| 日本精品一区二区三区蜜桃| 中文欧美无线码| 久久影院123| 久久国产精品影院| 久久精品熟女亚洲av麻豆精品| 男女免费视频国产| 精品欧美一区二区三区在线| 老司机亚洲免费影院| 丝袜在线中文字幕| 亚洲精品美女久久久久99蜜臀| 香蕉国产在线看| 色综合欧美亚洲国产小说| 色婷婷久久久亚洲欧美| avwww免费| 久久99一区二区三区| 亚洲综合色网址| 欧美 亚洲 国产 日韩一| 天天影视国产精品| 国产成人精品久久二区二区91| 色综合欧美亚洲国产小说| 色老头精品视频在线观看| 青春草视频在线免费观看| 交换朋友夫妻互换小说| 亚洲精品国产一区二区精华液| 精品福利观看| av视频免费观看在线观看| 麻豆av在线久日| 亚洲欧洲精品一区二区精品久久久| 黄色视频在线播放观看不卡| 考比视频在线观看| 久久久久国内视频| 久久久久精品国产欧美久久久 | 妹子高潮喷水视频| 久久国产亚洲av麻豆专区| 少妇人妻久久综合中文| 成人国产av品久久久| 一二三四社区在线视频社区8| 国产亚洲欧美在线一区二区| 永久免费av网站大全| 中文欧美无线码| 91麻豆av在线| 高潮久久久久久久久久久不卡| 高清欧美精品videossex| 日韩一区二区三区影片| 国产亚洲精品第一综合不卡| 久久亚洲精品不卡| 成人免费观看视频高清| 中亚洲国语对白在线视频| 久久精品久久久久久噜噜老黄| 在线av久久热| 久久国产精品人妻蜜桃| 夜夜骑夜夜射夜夜干| 天天躁狠狠躁夜夜躁狠狠躁| 大型av网站在线播放| 亚洲国产欧美在线一区| 久热爱精品视频在线9| 欧美老熟妇乱子伦牲交| 欧美精品一区二区免费开放| 啦啦啦视频在线资源免费观看| av线在线观看网站| 国产伦人伦偷精品视频| 成人三级做爰电影| 久久人人爽人人片av| 美女高潮喷水抽搐中文字幕| 中文字幕人妻熟女乱码| 91九色精品人成在线观看| 波多野结衣av一区二区av| 高清黄色对白视频在线免费看| 欧美国产精品一级二级三级| 亚洲精品第二区| 美女高潮喷水抽搐中文字幕| 久久精品久久久久久噜噜老黄| 搡老熟女国产l中国老女人| 日韩视频在线欧美| 中国美女看黄片| 天堂8中文在线网| 婷婷成人精品国产| 精品人妻熟女毛片av久久网站| 亚洲欧美成人综合另类久久久| 9191精品国产免费久久| 久久国产精品男人的天堂亚洲| 电影成人av| 国产日韩一区二区三区精品不卡| 老熟妇仑乱视频hdxx| 欧美日韩亚洲综合一区二区三区_| 可以免费在线观看a视频的电影网站| 亚洲午夜精品一区,二区,三区| 在线永久观看黄色视频| 国产一区二区三区av在线| 在线观看免费高清a一片| 97在线人人人人妻| 交换朋友夫妻互换小说| 国产成人av激情在线播放| av线在线观看网站| 天天操日日干夜夜撸| 国产高清videossex| 黄色 视频免费看| 亚洲性夜色夜夜综合| 美女扒开内裤让男人捅视频| 在线观看www视频免费| 亚洲精品美女久久av网站| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久成人av| 久久精品国产亚洲av香蕉五月 | 亚洲七黄色美女视频| 搡老乐熟女国产| 亚洲国产成人一精品久久久| 国内毛片毛片毛片毛片毛片| 自拍欧美九色日韩亚洲蝌蚪91| 一级a爱视频在线免费观看| 亚洲精品美女久久av网站| 久久久久久亚洲精品国产蜜桃av| 一级a爱视频在线免费观看| 久久久水蜜桃国产精品网| av片东京热男人的天堂| 成年人午夜在线观看视频| 少妇被粗大的猛进出69影院| 麻豆国产av国片精品| 9191精品国产免费久久| av视频免费观看在线观看| 99精品久久久久人妻精品| 亚洲专区字幕在线| 国产精品久久久人人做人人爽| 男人操女人黄网站| 一本大道久久a久久精品| 久久精品aⅴ一区二区三区四区| 国产av又大| 夜夜骑夜夜射夜夜干| 18禁国产床啪视频网站| 热99re8久久精品国产| www日本在线高清视频| 国产一卡二卡三卡精品| 国产国语露脸激情在线看| 女性被躁到高潮视频| 亚洲av美国av| 夜夜夜夜夜久久久久| 国产精品免费视频内射| 亚洲精品国产精品久久久不卡| 成年女人毛片免费观看观看9 | 黄片播放在线免费| 人人妻,人人澡人人爽秒播| cao死你这个sao货| av国产精品久久久久影院| 亚洲国产毛片av蜜桃av| 嫁个100分男人电影在线观看| 午夜福利视频精品| 久久久精品免费免费高清| 欧美日韩福利视频一区二区| www日本在线高清视频| www.999成人在线观看| 国产成人免费观看mmmm| 一边摸一边做爽爽视频免费| 美女脱内裤让男人舔精品视频| 一级毛片电影观看| 韩国精品一区二区三区| 国产黄色免费在线视频| 欧美日韩国产mv在线观看视频| 午夜免费观看性视频| 视频区图区小说| 国产精品一区二区精品视频观看| 欧美激情极品国产一区二区三区| 国产亚洲欧美精品永久| 每晚都被弄得嗷嗷叫到高潮| 亚洲黑人精品在线| 日日摸夜夜添夜夜添小说| 美女午夜性视频免费| 久久九九热精品免费| 精品久久久精品久久久| 宅男免费午夜| 女人精品久久久久毛片| 不卡一级毛片| 国产免费福利视频在线观看| 高潮久久久久久久久久久不卡| 国产真人三级小视频在线观看| 午夜精品久久久久久毛片777| 亚洲国产精品成人久久小说| 视频区欧美日本亚洲| 国产成人啪精品午夜网站| 老熟妇仑乱视频hdxx| 老鸭窝网址在线观看| 亚洲avbb在线观看| 欧美性长视频在线观看| 国产日韩欧美亚洲二区| 成年美女黄网站色视频大全免费| 2018国产大陆天天弄谢| 9热在线视频观看99| 欧美激情 高清一区二区三区| 日本a在线网址| 国产精品亚洲av一区麻豆| netflix在线观看网站| 999精品在线视频| 久久毛片免费看一区二区三区| 人成视频在线观看免费观看| 老司机午夜福利在线观看视频 | 亚洲精品国产一区二区精华液| 黑人操中国人逼视频| av线在线观看网站| 国产欧美亚洲国产| 在线观看免费午夜福利视频| 各种免费的搞黄视频| 王馨瑶露胸无遮挡在线观看| 少妇粗大呻吟视频| 久久久久久久精品精品| 国产精品一区二区精品视频观看| 国产男女超爽视频在线观看| 青青草视频在线视频观看| 国产亚洲av片在线观看秒播厂| 成人黄色视频免费在线看| 欧美日韩精品网址| 俄罗斯特黄特色一大片| 91老司机精品| 久久久国产欧美日韩av| 精品熟女少妇八av免费久了| 97精品久久久久久久久久精品| 免费观看人在逋| 亚洲欧美日韩高清在线视频 | 中国美女看黄片| 日韩免费高清中文字幕av| 十八禁网站免费在线| 成人三级做爰电影| 亚洲少妇的诱惑av| 午夜精品久久久久久毛片777| 久久国产精品影院| 成人亚洲精品一区在线观看| 天堂8中文在线网| 欧美黑人精品巨大| 国产97色在线日韩免费| 好男人电影高清在线观看| 亚洲精品日韩在线中文字幕| 老司机亚洲免费影院| 国产精品一区二区在线观看99| 天堂8中文在线网| 人妻 亚洲 视频| 各种免费的搞黄视频| 91字幕亚洲| 汤姆久久久久久久影院中文字幕| 不卡av一区二区三区| 久久中文看片网| 女人精品久久久久毛片| 九色亚洲精品在线播放| 国产一区二区三区综合在线观看| 91麻豆av在线| 亚洲精品中文字幕在线视频| 国产免费福利视频在线观看| 久久青草综合色| 日韩视频在线欧美| 国产成人影院久久av| 亚洲欧美一区二区三区久久| 人妻 亚洲 视频| 欧美日韩视频精品一区| 搡老熟女国产l中国老女人| 考比视频在线观看| 日韩有码中文字幕| 久久精品亚洲熟妇少妇任你| 国产伦理片在线播放av一区| 桃花免费在线播放| 午夜激情av网站| 亚洲中文av在线| 欧美黄色淫秽网站| 欧美激情久久久久久爽电影 | 黄网站色视频无遮挡免费观看| 亚洲,欧美精品.| 亚洲色图 男人天堂 中文字幕| 中文字幕高清在线视频| 日本五十路高清| av网站免费在线观看视频| 亚洲国产精品一区三区| 国产成人欧美| 国产男女超爽视频在线观看| 大码成人一级视频| 色婷婷久久久亚洲欧美| 黄片小视频在线播放| 热99re8久久精品国产| 老熟妇乱子伦视频在线观看 | 黄色a级毛片大全视频| 另类亚洲欧美激情| 亚洲九九香蕉| 国产男女内射视频| 亚洲精品国产av成人精品| 日本欧美视频一区| 成年美女黄网站色视频大全免费| 精品欧美一区二区三区在线| 国产精品国产三级国产专区5o| 亚洲精品美女久久久久99蜜臀| 国产麻豆69| www日本在线高清视频| 久久中文看片网| 三级毛片av免费| 成年人午夜在线观看视频| 精品少妇黑人巨大在线播放| bbb黄色大片| 五月天丁香电影| 老汉色av国产亚洲站长工具| 9色porny在线观看| 免费在线观看完整版高清| 欧美中文综合在线视频| 男人添女人高潮全过程视频| 成年女人毛片免费观看观看9 | 高清在线国产一区| 精品一品国产午夜福利视频| 欧美一级毛片孕妇| 波多野结衣av一区二区av| 十分钟在线观看高清视频www| 国产在线一区二区三区精| 热99re8久久精品国产| 男人舔女人的私密视频| 久久精品久久久久久噜噜老黄| 欧美日韩黄片免| 久久九九热精品免费| 性色av一级| 国产精品成人在线| 999久久久精品免费观看国产| 免费人妻精品一区二区三区视频| 日韩中文字幕欧美一区二区| 久久性视频一级片| 亚洲熟女毛片儿| 亚洲一区二区三区欧美精品| 一二三四在线观看免费中文在| 国产精品久久久久久人妻精品电影 | 亚洲精品av麻豆狂野| 亚洲国产日韩一区二区| 久久 成人 亚洲| av有码第一页| 国产免费一区二区三区四区乱码| 久久国产精品大桥未久av| 国产免费一区二区三区四区乱码| 叶爱在线成人免费视频播放| 无限看片的www在线观看| 狠狠婷婷综合久久久久久88av| 国产免费一区二区三区四区乱码| 丝袜美腿诱惑在线| 欧美国产精品va在线观看不卡| 久久精品国产亚洲av高清一级| 久久精品亚洲熟妇少妇任你| 狠狠狠狠99中文字幕| 欧美黄色片欧美黄色片| 欧美大码av| 国产精品国产av在线观看| √禁漫天堂资源中文www| 国产av一区二区精品久久| www.精华液| 美女中出高潮动态图| www.精华液| 热99久久久久精品小说推荐| 成人18禁高潮啪啪吃奶动态图| kizo精华| 黄色怎么调成土黄色| 亚洲成av片中文字幕在线观看| xxxhd国产人妻xxx| 色94色欧美一区二区| 啦啦啦视频在线资源免费观看| 日韩熟女老妇一区二区性免费视频| 正在播放国产对白刺激| av在线老鸭窝| 国产成人免费观看mmmm| 一区在线观看完整版| 极品人妻少妇av视频| 韩国高清视频一区二区三区| 久久人妻熟女aⅴ| 中文字幕av电影在线播放| 在线永久观看黄色视频| 成人三级做爰电影| 搡老岳熟女国产| 日本一区二区免费在线视频| 不卡av一区二区三区| 亚洲,欧美精品.| 免费av中文字幕在线| 国产又爽黄色视频| 99久久国产精品久久久| 久久久久久免费高清国产稀缺| 纵有疾风起免费观看全集完整版| 午夜91福利影院| 亚洲欧美色中文字幕在线| 免费不卡黄色视频| 韩国高清视频一区二区三区| 日韩视频一区二区在线观看| 国产男人的电影天堂91| 伊人久久大香线蕉亚洲五| 91国产中文字幕| 精品少妇内射三级| 久久久精品免费免费高清| 色综合欧美亚洲国产小说| 男女午夜视频在线观看| 操出白浆在线播放| 亚洲三区欧美一区| 69av精品久久久久久 | www.精华液| 亚洲精品av麻豆狂野| 亚洲av美国av| 亚洲精品国产色婷婷电影| 国产精品久久久久久精品古装| 久久中文字幕一级| 午夜日韩欧美国产| 人人妻,人人澡人人爽秒播| 99精国产麻豆久久婷婷| 久久国产亚洲av麻豆专区| 两个人免费观看高清视频| 在线天堂中文资源库| 狂野欧美激情性xxxx| 飞空精品影院首页| 亚洲国产欧美一区二区综合| 日日夜夜操网爽| 91麻豆精品激情在线观看国产 | 精品一区二区三区四区五区乱码| 啦啦啦啦在线视频资源| 欧美在线一区亚洲| netflix在线观看网站| 国产在线免费精品| 高潮久久久久久久久久久不卡| 蜜桃在线观看..| 国产激情久久老熟女| 日韩 欧美 亚洲 中文字幕| 两人在一起打扑克的视频| 狠狠精品人妻久久久久久综合| 日本av免费视频播放| 日韩中文字幕欧美一区二区| 国产精品 国内视频| 热99国产精品久久久久久7| 俄罗斯特黄特色一大片| 亚洲中文日韩欧美视频| 亚洲精品国产精品久久久不卡| 久久综合国产亚洲精品| 91麻豆精品激情在线观看国产 | 一边摸一边做爽爽视频免费| bbb黄色大片| 国产又爽黄色视频| 欧美久久黑人一区二区| 大陆偷拍与自拍| a级毛片黄视频| 久久九九热精品免费| 国产在线一区二区三区精| 免费久久久久久久精品成人欧美视频| 亚洲精品美女久久久久99蜜臀| 天天躁日日躁夜夜躁夜夜| 久久免费观看电影| 成年av动漫网址| 午夜激情久久久久久久| 男女床上黄色一级片免费看| 国产亚洲精品一区二区www | 交换朋友夫妻互换小说| 悠悠久久av| 亚洲精品乱久久久久久| 亚洲,欧美精品.| 精品卡一卡二卡四卡免费| 满18在线观看网站| 国产一级毛片在线| 精品欧美一区二区三区在线| 99re6热这里在线精品视频| 99九九在线精品视频| 最近最新免费中文字幕在线| 亚洲专区中文字幕在线| 91麻豆av在线| 亚洲精品美女久久av网站| av不卡在线播放| 国产精品偷伦视频观看了| 韩国精品一区二区三区| 黑人猛操日本美女一级片| 国产三级黄色录像| a 毛片基地| 欧美日韩福利视频一区二区| 国产淫语在线视频| 日日摸夜夜添夜夜添小说| 国产成人精品无人区| 老汉色∧v一级毛片| 各种免费的搞黄视频| 亚洲男人天堂网一区| 亚洲国产欧美网| 亚洲av欧美aⅴ国产| 三上悠亚av全集在线观看| 大片电影免费在线观看免费| 亚洲精品美女久久av网站| 在线十欧美十亚洲十日本专区| 制服诱惑二区| 国产xxxxx性猛交| 久久99热这里只频精品6学生| 51午夜福利影视在线观看| 国产精品久久久久久精品古装| 成人国产av品久久久| 中文字幕另类日韩欧美亚洲嫩草| 亚洲中文字幕日韩| 日韩一卡2卡3卡4卡2021年| 欧美午夜高清在线| 在线观看免费高清a一片| 一级片免费观看大全| 欧美+亚洲+日韩+国产| svipshipincom国产片| 日韩,欧美,国产一区二区三区| 热99国产精品久久久久久7| 亚洲天堂av无毛| 国产一区有黄有色的免费视频| 我的亚洲天堂| 精品免费久久久久久久清纯 | 1024视频免费在线观看| 久久久久久人人人人人| 夫妻午夜视频| videosex国产| 亚洲专区字幕在线| 国产成人精品无人区| 亚洲精品第二区| 成人手机av| 亚洲国产成人一精品久久久| 最黄视频免费看| 国产欧美日韩一区二区三区在线| 国产成人av教育| 嫩草影视91久久| 亚洲免费av在线视频| av在线app专区| 中文字幕色久视频| 美女中出高潮动态图| 国产高清videossex| 国产精品二区激情视频| av在线播放精品| 99国产综合亚洲精品| 无限看片的www在线观看| cao死你这个sao货| 亚洲三区欧美一区| 精品熟女少妇八av免费久了| 亚洲少妇的诱惑av| 老司机福利观看| 人妻 亚洲 视频| 人人妻人人添人人爽欧美一区卜| 日韩熟女老妇一区二区性免费视频| 午夜福利影视在线免费观看| 亚洲欧美一区二区三区黑人| 老司机在亚洲福利影院| 国产黄频视频在线观看| 老熟妇仑乱视频hdxx| 成人国语在线视频| 久热这里只有精品99| 国产一区二区 视频在线| 成人免费观看视频高清| 欧美 亚洲 国产 日韩一| 亚洲av国产av综合av卡| 天堂俺去俺来也www色官网| 国产成人欧美在线观看 | 蜜桃在线观看..| 老司机深夜福利视频在线观看 | 美女午夜性视频免费| 亚洲 欧美一区二区三区| 欧美精品人与动牲交sv欧美| 十八禁高潮呻吟视频| 久久综合国产亚洲精品| 交换朋友夫妻互换小说| 精品人妻一区二区三区麻豆| 黑人欧美特级aaaaaa片| cao死你这个sao货| 极品少妇高潮喷水抽搐| 人人澡人人妻人| 美女大奶头黄色视频| 国产野战对白在线观看| 日日夜夜操网爽| 日韩大码丰满熟妇| 亚洲第一av免费看| 91精品伊人久久大香线蕉| 91国产中文字幕| 亚洲精品久久成人aⅴ小说| 人妻人人澡人人爽人人| av福利片在线| 国产免费av片在线观看野外av| 欧美人与性动交α欧美软件| 两性夫妻黄色片| 国产av一区二区精品久久|