• <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ù)值模擬的若干影響因素
    一本久久精品| 免费黄色在线免费观看| 9热在线视频观看99| 大码成人一级视频| 久久亚洲国产成人精品v| 欧美日韩av久久| 国产精品一区二区在线不卡| 亚洲av欧美aⅴ国产| 在线观看免费视频网站a站| 国产视频首页在线观看| 在线天堂最新版资源| 国产乱来视频区| 最近最新中文字幕免费大全7| 国产伦理片在线播放av一区| 久久久久精品性色| 亚洲内射少妇av| 免费高清在线观看视频在线观看| 91精品三级在线观看| 欧美成人精品欧美一级黄| 考比视频在线观看| 亚洲成人手机| 啦啦啦在线观看免费高清www| 五月开心婷婷网| 波野结衣二区三区在线| 成人免费观看视频高清| 成人毛片60女人毛片免费| 国产免费福利视频在线观看| 波野结衣二区三区在线| 久久精品国产鲁丝片午夜精品| 9色porny在线观看| 国产精品国产三级国产av玫瑰| 色哟哟·www| 9色porny在线观看| 亚洲综合精品二区| 99热这里只有是精品在线观看| 韩国精品一区二区三区 | 国产成人精品久久久久久| 国产伦理片在线播放av一区| 精品国产国语对白av| 搡女人真爽免费视频火全软件| 免费av不卡在线播放| 久久精品久久久久久久性| 中文字幕精品免费在线观看视频 | 亚洲五月色婷婷综合| 亚洲伊人久久精品综合| 欧美人与善性xxx| 最近的中文字幕免费完整| 久久影院123| 免费高清在线观看视频在线观看| 亚洲经典国产精华液单| 97在线人人人人妻| 国产老妇伦熟女老妇高清| 黄片播放在线免费| 大香蕉久久成人网| 亚洲性久久影院| 午夜福利视频在线观看免费| 五月天丁香电影| videossex国产| 亚洲欧美清纯卡通| 我要看黄色一级片免费的| 欧美日韩国产mv在线观看视频| 亚洲欧美日韩卡通动漫| 一级毛片电影观看| 人人妻人人爽人人添夜夜欢视频| 久久人人97超碰香蕉20202| 久久影院123| 大香蕉久久成人网| 国产精品人妻久久久久久| 性色av一级| 久久午夜综合久久蜜桃| 午夜福利视频精品| 欧美xxxx性猛交bbbb| 丝瓜视频免费看黄片| 18禁国产床啪视频网站| 日本-黄色视频高清免费观看| 蜜桃在线观看..| 欧美亚洲日本最大视频资源| 人妻人人澡人人爽人人| 人人妻人人澡人人看| 欧美激情极品国产一区二区三区 | 男女免费视频国产| 18禁观看日本| 国产日韩欧美视频二区| 国产男人的电影天堂91| 香蕉丝袜av| 少妇精品久久久久久久| 春色校园在线视频观看| 91aial.com中文字幕在线观看| 五月开心婷婷网| 男人添女人高潮全过程视频| 午夜视频国产福利| 午夜福利,免费看| 伦精品一区二区三区| 久久av网站| 成人二区视频| 亚洲国产最新在线播放| 国产成人精品婷婷| 女性生殖器流出的白浆| 丝袜喷水一区| av电影中文网址| 日韩av免费高清视频| 97在线视频观看| 美国免费a级毛片| 亚洲国产成人一精品久久久| 18+在线观看网站| 最近最新中文字幕大全免费视频 | 最新中文字幕久久久久| 国产av精品麻豆| 亚洲国产看品久久| 久久精品国产亚洲av天美| 夜夜骑夜夜射夜夜干| 男女边吃奶边做爰视频| 超色免费av| 国产男女超爽视频在线观看| 色5月婷婷丁香| 国产色婷婷99| 在线亚洲精品国产二区图片欧美| 老司机影院成人| av在线播放精品| 成人影院久久| 精品久久久久久电影网| 欧美精品一区二区免费开放| 69精品国产乱码久久久| 欧美少妇被猛烈插入视频| 国产片特级美女逼逼视频| av在线老鸭窝| 国产成人一区二区在线| 全区人妻精品视频| 最新中文字幕久久久久| 97在线视频观看| 美女福利国产在线| 韩国av在线不卡| 97在线人人人人妻| 高清视频免费观看一区二区| 精品国产露脸久久av麻豆| 国产一区二区在线观看日韩| 午夜久久久在线观看| 亚洲精品色激情综合| 少妇人妻久久综合中文| 久久99蜜桃精品久久| 1024视频免费在线观看| 国产国语露脸激情在线看| 2018国产大陆天天弄谢| 日韩成人av中文字幕在线观看| 少妇被粗大猛烈的视频| 大陆偷拍与自拍| 欧美变态另类bdsm刘玥| 欧美 日韩 精品 国产| 咕卡用的链子| 五月玫瑰六月丁香| 免费不卡的大黄色大毛片视频在线观看| 免费女性裸体啪啪无遮挡网站| 丁香六月天网| 国语对白做爰xxxⅹ性视频网站| 又粗又硬又长又爽又黄的视频| 国产极品粉嫩免费观看在线| av播播在线观看一区| 欧美国产精品一级二级三级| 国产老妇伦熟女老妇高清| 亚洲欧美色中文字幕在线| 美女脱内裤让男人舔精品视频| 精品一区二区三区四区五区乱码 | 欧美人与性动交α欧美软件 | 亚洲一级一片aⅴ在线观看| 久久午夜福利片| 亚洲欧美中文字幕日韩二区| 精品人妻一区二区三区麻豆| 久久久国产精品麻豆| 国产男人的电影天堂91| 国产成人aa在线观看| 日韩视频在线欧美| 十八禁网站网址无遮挡| 丝袜在线中文字幕| 亚洲高清免费不卡视频| 国产亚洲午夜精品一区二区久久| 最新中文字幕久久久久| 久久免费观看电影| 波多野结衣一区麻豆| 国产毛片在线视频| 久久久欧美国产精品| 欧美老熟妇乱子伦牲交| 在线精品无人区一区二区三| 99国产精品免费福利视频| 亚洲成av片中文字幕在线观看 | 好男人视频免费观看在线| 在线免费观看不下载黄p国产| a级毛色黄片| 免费少妇av软件| 国产日韩欧美在线精品| 国产在线视频一区二区| 男女下面插进去视频免费观看 | 亚洲av福利一区| 日日撸夜夜添| 国产国语露脸激情在线看| 哪个播放器可以免费观看大片| 色视频在线一区二区三区| 蜜桃国产av成人99| 亚洲精品美女久久av网站| 五月伊人婷婷丁香| 一二三四中文在线观看免费高清| 国产 精品1| 18禁动态无遮挡网站| 岛国毛片在线播放| 亚洲av成人精品一二三区| 亚洲高清免费不卡视频| 91精品伊人久久大香线蕉| 国产国拍精品亚洲av在线观看| 在线观看免费日韩欧美大片| 狠狠婷婷综合久久久久久88av| 国产成人一区二区在线| 18禁国产床啪视频网站| 22中文网久久字幕| 亚洲成人av在线免费| 久久精品国产亚洲av涩爱| 男人添女人高潮全过程视频| 免费黄频网站在线观看国产| 自线自在国产av| 国产深夜福利视频在线观看| a 毛片基地| 午夜福利影视在线免费观看| 在线观看免费高清a一片| 久久国产亚洲av麻豆专区| 亚洲色图综合在线观看| 午夜免费鲁丝| av国产精品久久久久影院| 久久 成人 亚洲| 国产免费视频播放在线视频| 黄色毛片三级朝国网站| 欧美日韩视频高清一区二区三区二| 色婷婷久久久亚洲欧美| 欧美丝袜亚洲另类| 巨乳人妻的诱惑在线观看| 免费大片黄手机在线观看| 看非洲黑人一级黄片| 丰满迷人的少妇在线观看| 久久精品国产综合久久久 | 国产不卡av网站在线观看| 人妻少妇偷人精品九色| 毛片一级片免费看久久久久| 人人澡人人妻人| 亚洲精品国产色婷婷电影| 国产亚洲欧美精品永久| 精品午夜福利在线看| 我的女老师完整版在线观看| 日本vs欧美在线观看视频| 搡老乐熟女国产| 精品一区在线观看国产| 国产国拍精品亚洲av在线观看| 国产精品三级大全| 国产精品三级大全| 自拍欧美九色日韩亚洲蝌蚪91| 又粗又硬又长又爽又黄的视频| 午夜久久久在线观看| √禁漫天堂资源中文www| 欧美精品国产亚洲| 国产欧美另类精品又又久久亚洲欧美| 五月玫瑰六月丁香| 多毛熟女@视频| 在线观看美女被高潮喷水网站| 18禁在线无遮挡免费观看视频| 精品久久久精品久久久| 亚洲精品日本国产第一区| 51国产日韩欧美| av女优亚洲男人天堂| 国产精品麻豆人妻色哟哟久久| 肉色欧美久久久久久久蜜桃| 日韩中文字幕视频在线看片| 亚洲综合精品二区| 丁香六月天网| 日本-黄色视频高清免费观看| 亚洲精品乱久久久久久| 国产男女内射视频| 国国产精品蜜臀av免费| 啦啦啦啦在线视频资源| 国产精品.久久久| 一本色道久久久久久精品综合| 97人妻天天添夜夜摸| 亚洲成av片中文字幕在线观看 | 亚洲欧美色中文字幕在线| 精品一区二区三卡| 久热这里只有精品99| 国产一区有黄有色的免费视频| 精品一区二区三卡| 国产av一区二区精品久久| 99re6热这里在线精品视频| 国产免费现黄频在线看| 一本色道久久久久久精品综合| 一二三四在线观看免费中文在 | 日日撸夜夜添| 国产精品99久久99久久久不卡 | 婷婷成人精品国产| 99精国产麻豆久久婷婷| 免费av中文字幕在线| 纯流量卡能插随身wifi吗| 永久网站在线| 欧美性感艳星| 精品99又大又爽又粗少妇毛片| 亚洲国产日韩一区二区| 我的女老师完整版在线观看| 亚洲精品中文字幕在线视频| 亚洲精品视频女| www日本在线高清视频| videosex国产| 啦啦啦在线观看免费高清www| 老司机亚洲免费影院| 美女中出高潮动态图| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 美女内射精品一级片tv| 一二三四中文在线观看免费高清| 一本久久精品| 亚洲欧美一区二区三区国产| 18禁动态无遮挡网站| 亚洲国产看品久久| 岛国毛片在线播放| 日本黄大片高清| 在线观看免费日韩欧美大片| 99久久中文字幕三级久久日本| 精品久久国产蜜桃| 卡戴珊不雅视频在线播放| 午夜福利影视在线免费观看| 少妇人妻精品综合一区二区| 成年av动漫网址| 久久精品国产亚洲av天美| 国产色婷婷99| 亚洲中文av在线| 久久99蜜桃精品久久| 中国三级夫妇交换| 视频中文字幕在线观看| 亚洲四区av| 国产69精品久久久久777片| 韩国精品一区二区三区 | 欧美精品亚洲一区二区| 国产一区二区激情短视频 | 免费少妇av软件| 香蕉丝袜av| 日日撸夜夜添| 夫妻午夜视频| 黄色毛片三级朝国网站| 国产老妇伦熟女老妇高清| 亚洲精品456在线播放app| 一级毛片黄色毛片免费观看视频| 国产视频首页在线观看| 国产日韩欧美亚洲二区| 国产黄频视频在线观看| 最近手机中文字幕大全| 99re6热这里在线精品视频| 777米奇影视久久| 亚洲国产成人一精品久久久| 交换朋友夫妻互换小说| 毛片一级片免费看久久久久| 赤兔流量卡办理| 亚洲av中文av极速乱| 日韩不卡一区二区三区视频在线| 亚洲内射少妇av| 又黄又粗又硬又大视频| 韩国精品一区二区三区 | 亚洲综合色网址| av国产精品久久久久影院| videos熟女内射| 不卡视频在线观看欧美| 亚洲精品视频女| 日日撸夜夜添| 黑人欧美特级aaaaaa片| 国产成人a∨麻豆精品| 水蜜桃什么品种好| 蜜桃国产av成人99| 日韩成人伦理影院| 亚洲激情五月婷婷啪啪| 美女国产视频在线观看| 久久韩国三级中文字幕| 51国产日韩欧美| tube8黄色片| 久久精品人人爽人人爽视色| av播播在线观看一区| 国产亚洲最大av| 五月天丁香电影| 午夜影院在线不卡| 欧美日韩综合久久久久久| 五月伊人婷婷丁香| 91久久精品国产一区二区三区| 一二三四在线观看免费中文在 | 国产精品国产三级国产av玫瑰| 天堂俺去俺来也www色官网| 五月开心婷婷网| 国产精品人妻久久久影院| 亚洲欧美日韩另类电影网站| 亚洲国产精品999| 99热网站在线观看| 精品一区二区三区视频在线| 成人午夜精彩视频在线观看| 美女主播在线视频| 一级片'在线观看视频| 性色avwww在线观看| 久久人人97超碰香蕉20202| 国产一区二区激情短视频 | 日本色播在线视频| 最近2019中文字幕mv第一页| 人体艺术视频欧美日本| 久久青草综合色| 国产av一区二区精品久久| 亚洲成国产人片在线观看| 男人舔女人的私密视频| av播播在线观看一区| a 毛片基地| 两个人免费观看高清视频| 国产免费一级a男人的天堂| 十分钟在线观看高清视频www| 久久97久久精品| 一二三四在线观看免费中文在 | av在线老鸭窝| 国产免费现黄频在线看| 午夜福利网站1000一区二区三区| 少妇人妻 视频| 日韩中文字幕视频在线看片| 18禁裸乳无遮挡动漫免费视频| 观看美女的网站| 少妇 在线观看| 久久精品国产亚洲av天美| 少妇被粗大的猛进出69影院 | 日本黄大片高清| 亚洲欧美成人综合另类久久久| 在线观看三级黄色| videossex国产| 精品亚洲成国产av| a级片在线免费高清观看视频| 亚洲精品视频女| 黄色怎么调成土黄色| 日韩制服丝袜自拍偷拍| 寂寞人妻少妇视频99o| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品嫩草影院av在线观看| 国产精品久久久av美女十八| 日韩伦理黄色片| 在线亚洲精品国产二区图片欧美| 成人手机av| 久久精品夜色国产| 色网站视频免费| 国产精品免费大片| 国产在线一区二区三区精| 久久免费观看电影| av国产久精品久网站免费入址| 国产又色又爽无遮挡免| 亚洲欧洲国产日韩| 国产精品99久久99久久久不卡 | 免费久久久久久久精品成人欧美视频 | 全区人妻精品视频| 美女国产视频在线观看| 久久午夜福利片| 国产精品麻豆人妻色哟哟久久| 22中文网久久字幕| 两个人看的免费小视频| 麻豆乱淫一区二区| av.在线天堂| 99九九在线精品视频| 成人免费观看视频高清| 啦啦啦中文免费视频观看日本| 亚洲,欧美精品.| 国产高清国产精品国产三级| 亚洲人成网站在线观看播放| 一级片'在线观看视频| 麻豆精品久久久久久蜜桃| a级毛片在线看网站| 丰满迷人的少妇在线观看| 自线自在国产av| 亚洲av电影在线进入| 国产成人免费观看mmmm| 精品少妇内射三级| 精品人妻一区二区三区麻豆| 一级片免费观看大全| 久久久欧美国产精品| 美女大奶头黄色视频| 久久综合国产亚洲精品| 亚洲人成网站在线观看播放| av一本久久久久| 亚洲人成网站在线观看播放| 高清av免费在线| 乱码一卡2卡4卡精品| 国产熟女午夜一区二区三区| a级毛色黄片| 亚洲精品乱码久久久久久按摩| 美国免费a级毛片| 亚洲av电影在线进入| 久久99热这里只频精品6学生| 人妻 亚洲 视频| 日本色播在线视频| 国产精品99久久99久久久不卡 | 久久国产亚洲av麻豆专区| tube8黄色片| 如日韩欧美国产精品一区二区三区| 免费日韩欧美在线观看| 大话2 男鬼变身卡| av不卡在线播放| 成年美女黄网站色视频大全免费| 夜夜骑夜夜射夜夜干| 午夜福利在线观看免费完整高清在| 99久国产av精品国产电影| 在线看a的网站| 久久久久精品性色| 美女大奶头黄色视频| 久久久久国产精品人妻一区二区| 亚洲av中文av极速乱| 国产乱来视频区| 两性夫妻黄色片 | 老司机影院毛片| 丰满少妇做爰视频| 自线自在国产av| 高清黄色对白视频在线免费看| 黄片播放在线免费| 成人综合一区亚洲| 精品国产一区二区久久| 51国产日韩欧美| 亚洲第一区二区三区不卡| 最近中文字幕高清免费大全6| 午夜免费观看性视频| 菩萨蛮人人尽说江南好唐韦庄| 少妇 在线观看| 色网站视频免费| 丰满迷人的少妇在线观看| 日本91视频免费播放| 香蕉国产在线看| 欧美国产精品va在线观看不卡| 亚洲国产日韩一区二区| 有码 亚洲区| 丁香六月天网| 最黄视频免费看| 9热在线视频观看99| 最近中文字幕2019免费版| 亚洲欧洲精品一区二区精品久久久 | 高清毛片免费看| 91午夜精品亚洲一区二区三区| 欧美日韩国产mv在线观看视频| 91精品三级在线观看| 伦理电影免费视频| 男人添女人高潮全过程视频| 精品一区二区免费观看| 亚洲国产欧美日韩在线播放| 一级,二级,三级黄色视频| 日本欧美国产在线视频| 成人免费观看视频高清| 人人妻人人添人人爽欧美一区卜| kizo精华| 亚洲,一卡二卡三卡| 99国产综合亚洲精品| 精品一区二区三区视频在线| 国产精品 国内视频| 99热全是精品| 一级,二级,三级黄色视频| 丝袜脚勾引网站| 蜜臀久久99精品久久宅男| a级毛色黄片| 大香蕉久久成人网| 欧美3d第一页| 日本与韩国留学比较| 国产69精品久久久久777片| freevideosex欧美| 曰老女人黄片| 国产精品 国内视频| 热99国产精品久久久久久7| 大陆偷拍与自拍| 伦理电影免费视频| 另类精品久久| 亚洲一区二区三区欧美精品| 亚洲婷婷狠狠爱综合网| 亚洲一级一片aⅴ在线观看| 午夜影院在线不卡| 丰满乱子伦码专区| 捣出白浆h1v1| 人妻人人澡人人爽人人| 男的添女的下面高潮视频| a级毛色黄片| 人人妻人人澡人人看| 久久久久久久亚洲中文字幕| 一级毛片电影观看| 欧美日本中文国产一区发布| 亚洲国产最新在线播放| 天天操日日干夜夜撸| 欧美丝袜亚洲另类| 精品久久久久久电影网| 青春草视频在线免费观看| 久久热在线av| 久久女婷五月综合色啪小说| 中文乱码字字幕精品一区二区三区| 精品人妻偷拍中文字幕| 丝袜脚勾引网站| 日韩人妻精品一区2区三区| 国产精品久久久久久精品电影小说| 亚洲欧美一区二区三区黑人 | 日韩欧美精品免费久久| 全区人妻精品视频| 一区二区三区乱码不卡18| 成人亚洲欧美一区二区av| 18禁在线无遮挡免费观看视频| 久久精品人人爽人人爽视色| 美国免费a级毛片| 大陆偷拍与自拍| 丰满少妇做爰视频| 99精国产麻豆久久婷婷| 午夜老司机福利剧场| 插逼视频在线观看| 亚洲欧美清纯卡通| 中文字幕另类日韩欧美亚洲嫩草| 国产精品熟女久久久久浪| 交换朋友夫妻互换小说| 黄色视频在线播放观看不卡| 视频在线观看一区二区三区| 国产亚洲一区二区精品| 男人操女人黄网站| 自线自在国产av| 少妇的丰满在线观看| 久久久精品免费免费高清| 91久久精品国产一区二区三区| 人人妻人人爽人人添夜夜欢视频| 国精品久久久久久国模美| 老司机影院毛片| 2021少妇久久久久久久久久久|