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

    三種跨聲速洞壁干擾修正方法及其在小展弦比飛翼標(biāo)模試驗(yàn)中的應(yīng)用

    2016-04-06 03:02:54鐘世東蘇繼川李永紅
    關(guān)鍵詞:洞壁飛翼迎角

    鐘世東,李 巍,蘇繼川,李永紅,賀 中

    (中國(guó)空氣動(dòng)力研究與發(fā)展中心高速空氣動(dòng)力研究所,四川綿陽(yáng) 621000)

    三種跨聲速洞壁干擾修正方法及其在小展弦比飛翼標(biāo)模試驗(yàn)中的應(yīng)用

    鐘世東*,李 巍,蘇繼川,李永紅,賀 中

    (中國(guó)空氣動(dòng)力研究與發(fā)展中心高速空氣動(dòng)力研究所,四川綿陽(yáng) 621000)

    在FL-24風(fēng)洞進(jìn)行了帶壁壓信息測(cè)量的小展弦比飛翼標(biāo)模測(cè)力試驗(yàn),并在FL-26風(fēng)洞進(jìn)行了洞壁干擾驗(yàn)證試驗(yàn)。本文利用小擾動(dòng)位流壁壓信息法、全速勢(shì)位流方法、基于RANS的壁壓信息法三種方法開(kāi)展飛翼標(biāo)模的洞壁干擾修正研究,并與試驗(yàn)結(jié)果進(jìn)行了比較。結(jié)果表明,飛翼標(biāo)模洞壁干擾呈現(xiàn)與翼身尾常規(guī)布局不同的規(guī)律,三種方法在飛翼標(biāo)模洞壁干擾修正中有各自的適用性。

    飛翼標(biāo)模;洞壁干擾;數(shù)值模擬

    0 引 言

    高速洞壁干擾修正技術(shù)歷經(jīng)發(fā)展,概括起來(lái)有以下修正方法:純?cè)囼?yàn)修正(單模型大小風(fēng)洞對(duì)比法、多模型單風(fēng)洞外推法)、經(jīng)驗(yàn)修正(統(tǒng)計(jì)數(shù)據(jù)類(lèi)推法)、映象法和Maskwell法、有限基本解法(渦格法、面元法)、壁壓信息法(雙參數(shù)壁壓信息積分法、單參數(shù)壁壓信息法)和基于CFD的非線性修正法等。

    國(guó)內(nèi)從20世紀(jì)80年代開(kāi)始,采用壁壓法對(duì)高速試驗(yàn)洞壁干擾進(jìn)行修正,范召林等在FL-21風(fēng)洞中先后進(jìn)行了一系列二維和三維的風(fēng)洞試驗(yàn)[1],研制出了適用于不同風(fēng)洞的壁壓測(cè)量裝置,發(fā)展了一系列適用于各種高亞聲速風(fēng)洞試驗(yàn)的線性洞壁干擾修正方法[23]。90年代,中國(guó)空氣動(dòng)力研究與發(fā)展中心高速所與總體技術(shù)部針對(duì)Ma=1附近洞壁干擾修正這一難題,開(kāi)展了跨聲速三維非線性洞壁干擾修正的研究[4-6],采用Euler方程和N-S方程模擬模型繞流場(chǎng),以實(shí)測(cè)的透氣壁附近的壓力分布作為風(fēng)洞流場(chǎng)的邊界條件,數(shù)值求解風(fēng)洞流場(chǎng)和自由流場(chǎng),由兩者之差得到洞壁干擾對(duì)模型氣動(dòng)力的影響,取得了令人滿(mǎn)意的結(jié)果。此后,對(duì)數(shù)值模擬方法和網(wǎng)格生成技術(shù)的不斷完善,將其應(yīng)用范圍不斷擴(kuò)展[7]。

    2000年以后,國(guó)內(nèi)外對(duì)洞壁干擾的研究主要圍繞基于CFD的非線性修正方法。GLASGOW大學(xué)計(jì)算流體動(dòng)力試驗(yàn)室利用RANS方程求解器研究洞壁對(duì)三角翼渦產(chǎn)生、發(fā)展與破裂的影響[8]。日本空間探索局(JAAX)用含壁壓信息的CFD非線性模擬手段研究ONERA M5翼身組合體洞壁干擾[9]。DNW采用CFD方法完全模擬了HST風(fēng)洞的試驗(yàn)段、駐室和擴(kuò)散段,并研究了不同尺度的DLR-F4機(jī)翼洞壁干擾特性[10]。

    此外,波音公司一直致力于洞壁干擾的工程化修正方法,從早期的面元法程序PANAIR[11]到全速勢(shì)方程求解器TRANAIR。利用TRANAIR開(kāi)展了翼身組合體標(biāo)模不同翼展模型洞壁干擾擾動(dòng)速度和擾動(dòng)迎角[12],試圖建立洞壁干擾數(shù)據(jù)庫(kù),并用CFD++驗(yàn)證其不確定度[13]。

    飛翼布局因高效、隱身等優(yōu)點(diǎn)廣泛地應(yīng)用于現(xiàn)代飛行器設(shè)計(jì)中,而針對(duì)飛翼布局洞壁干擾修正方法的研究尚未開(kāi)展。飛翼布局洞壁干擾與翼-身-尾常規(guī)布局是否有相同的規(guī)律性,傳統(tǒng)方法的適用性如何,都值得深入分析。裝備預(yù)研聯(lián)合攻關(guān)項(xiàng)目“風(fēng)洞試驗(yàn)技術(shù)”以小展弦比飛翼布局標(biāo)模為研究對(duì)象,廣泛開(kāi)展相關(guān)試驗(yàn)技術(shù)的研究[14-16],其中包括洞壁干擾修正方法的研究。

    為摸清飛翼標(biāo)模洞壁干擾影響規(guī)律并研究洞壁干擾修正方法對(duì)飛翼布局的適用性,該項(xiàng)目在FL-24風(fēng)洞進(jìn)行了帶壁壓測(cè)量的測(cè)力試驗(yàn),并在FL-26風(fēng)洞進(jìn)行了驗(yàn)證試驗(yàn)。本文除了利用小擾動(dòng)位流壁壓信息法和基于RANS的壁壓信息法計(jì)算飛翼標(biāo)模在FL-24風(fēng)洞中的洞壁干擾修正量,還首次在國(guó)內(nèi)建立了與波音公司類(lèi)似的全速勢(shì)位流方法并實(shí)現(xiàn)了洞壁干擾修正。通過(guò)對(duì)三種方法修正結(jié)果與試驗(yàn)結(jié)果的比較,綜合評(píng)價(jià)三種方法的優(yōu)缺點(diǎn)及在飛翼標(biāo)模洞壁干擾修正中的適用性,為類(lèi)似飛翼布局風(fēng)洞試驗(yàn)洞壁干擾修正提供參考。

    1 飛翼標(biāo)模風(fēng)洞試驗(yàn)

    1.1 試驗(yàn)設(shè)備與方法

    飛翼標(biāo)模采用單前緣融合體飛翼布局,前緣后掠角65°,后緣后掠角47°。試驗(yàn)?zāi)P蜑榭s比1∶19的全金屬模型,機(jī)身長(zhǎng)度為0.806m,翼展0.602m,平均氣動(dòng)弦長(zhǎng)為0.5032m。

    中國(guó)空氣動(dòng)力研究與發(fā)展中心高速所FL-24風(fēng)洞試驗(yàn)段尺寸1.2m×1.2m,跨聲速試驗(yàn)段上下壁面為60°斜孔壁,開(kāi)閉比為4.3%;左右為直孔壁,開(kāi)閉比為21.4%;孔徑和壁厚為12mm。飛翼標(biāo)模采用尾支撐,前緣5%當(dāng)?shù)叵议L(zhǎng)處粘貼高度為0.12mm的柱狀轉(zhuǎn)捩帶,零迎角時(shí)在FL-24風(fēng)洞的堵塞度約為0.97%。試驗(yàn)分為無(wú)壁壓管測(cè)力試驗(yàn)、安裝壁壓管空風(fēng)洞試驗(yàn)和安裝壁壓管帶模型測(cè)力試驗(yàn)。修正氣流偏角后的無(wú)壁壓條測(cè)力結(jié)果作為飛翼標(biāo)模在FL-24風(fēng)洞中的基準(zhǔn)數(shù)據(jù),帶模型壁壓減去相同Ma數(shù)空風(fēng)洞壁壓作為洞壁干擾修正用的壁壓數(shù)據(jù)。共安裝有12根壁壓條,上下洞壁各4根、左右洞壁各2根,長(zhǎng)度覆蓋了從試驗(yàn)段入口到支架段區(qū)域。圖1給出了試驗(yàn)照片。

    圖1 FL-24風(fēng)洞試驗(yàn)照片F(xiàn)ig.1 FL-24wind tunnel test picture

    圖2 FL-26風(fēng)洞試驗(yàn)照片F(xiàn)ig.2 FL-26wind tunnel test picture

    中國(guó)空氣動(dòng)力研究與發(fā)展中心高速所FL-26風(fēng)洞試驗(yàn)段尺寸為2.4m×2.4m,四壁都為斜孔壁,綜合開(kāi)閉比為4.3%。飛翼標(biāo)模采用FL-24風(fēng)洞相同的天平、支桿通過(guò)過(guò)渡接頭安裝在中部支架上,零迎角時(shí)堵塞度約為0.24%。修正氣流偏角后的測(cè)力結(jié)果近似認(rèn)為無(wú)洞壁干擾數(shù)據(jù),與FL-24風(fēng)洞數(shù)據(jù)之差作為FL-24風(fēng)洞洞壁干擾試驗(yàn)結(jié)果。

    1.2 壁壓信息

    扣除空風(fēng)洞后的壁壓數(shù)據(jù)由于測(cè)壓孔加工質(zhì)量等原因仍會(huì)存在一定程度的跳動(dòng)散布,可采用三次樣條曲線進(jìn)行光順,然后使用樣條插值法將壁壓數(shù)據(jù)插值到風(fēng)洞壁面處。圖3給出了Ma=0.8,壁壓信息插值到整個(gè)風(fēng)洞洞壁位置的壓力系數(shù)分布圖??梢钥吹?,迎角4°時(shí),模型上方為氣流加速區(qū),模型下方為減速區(qū);迎角16°時(shí),加速區(qū)和減速區(qū)范圍增加、強(qiáng)度增強(qiáng),并在側(cè)壁上出現(xiàn)近似與模型平行的分界區(qū)域。

    圖3 FL-24風(fēng)洞壁壓信息Fig.3 Wall pressure signature in FL-24wind tunnel

    1.3 試驗(yàn)結(jié)果分析

    飛翼標(biāo)模氣動(dòng)特性是以渦的產(chǎn)生、發(fā)展與渦破裂為主導(dǎo)的,在跨聲速階段還伴隨著渦與激波的交叉干擾。飛翼標(biāo)模氣動(dòng)特性詳細(xì)描述見(jiàn)參考文獻(xiàn)[16],這里只給出FL-24風(fēng)洞和FL-26風(fēng)洞試驗(yàn)結(jié)果升力特性對(duì)比曲線。升力曲線可分為線性段、非線性渦升力段、渦破裂段和渦破裂后段。小迎角附著流狀態(tài)為線性段,迎角4°左右開(kāi)始產(chǎn)生前緣渦升力,升力系數(shù)呈非線性增加,12°左右渦破裂,升力曲線出現(xiàn)拐點(diǎn),之后升力系數(shù)繼續(xù)增加。

    圖5給出了FL-26、FL-24風(fēng)洞試驗(yàn)結(jié)果升力系數(shù)差量曲線。兩風(fēng)洞試驗(yàn)采用相同的天平與支桿,保證近場(chǎng)支撐干擾相似,試驗(yàn)結(jié)果均進(jìn)行了氣流偏角修正,可以認(rèn)為兩風(fēng)洞數(shù)據(jù)差異是由洞壁干擾引起的??梢钥闯?,飛翼標(biāo)模洞壁干擾量與翼-身-尾組合體常規(guī)布局不同,而表現(xiàn)為隨迎角變化有不同的發(fā)展趨勢(shì):在小迎角階段,洞壁干擾量基本線性且隨迎角增加向負(fù)向增加,在中等迎角開(kāi)始,趨勢(shì)反向,隨迎角增加洞壁干擾量向正向增加。這與飛翼標(biāo)模氣動(dòng)特性為渦主導(dǎo)的特點(diǎn)是密切相關(guān)的,洞壁干擾在線性段、渦升力非線性段、渦破裂不同階段有不同的影響。需要指出的是,兩座風(fēng)洞升力曲線渦破裂迎角有差異,而渦破裂影響因素不僅與洞壁干擾相關(guān),還與流場(chǎng)湍流度、不均勻度和遠(yuǎn)場(chǎng)支架干擾[15]等因素都有關(guān)。渦破裂后的氣動(dòng)特性是與渦破裂起始迎角密切相關(guān)的,因此渦破裂后兩風(fēng)洞試驗(yàn)結(jié)果的差量不能完全代表洞壁干擾影響,而是摻雜了其他因素的影響。

    圖4 FL-24與FL-26風(fēng)洞試驗(yàn)結(jié)果比較(M=0.8)Fig.4 Comparison of test results in FL-24and FL-26

    2 洞壁干擾修正方法

    2.1 小擾動(dòng)位流壁壓信息法

    小擾動(dòng)位流壁壓信息法是基于線化小擾動(dòng)位勢(shì)方程,根據(jù)模型、支桿等物面條件得到無(wú)洞壁存在情況下模型、支桿等網(wǎng)格的基本解強(qiáng)度,然后計(jì)算模型基本解與洞壁基本解在洞壁網(wǎng)格控制點(diǎn)出的誘導(dǎo)速度,根據(jù)壁面壓力信息代替的洞壁邊界條件,求解洞壁上基本解強(qiáng)度,從而得到洞壁對(duì)模型的干擾壓力分布和干擾阻力,迭代該過(guò)程得到洞壁干擾的影響量。壁壓條測(cè)量壓力系數(shù)Cp,洞壁邊界條件為:

    式中φm為模型引起的攏動(dòng)速勢(shì)。該邊界條件與控制方程構(gòu)成一個(gè)三維Dirichlet問(wèn)題,可以采用偶極子面元法和有限差分法求解。

    該方法經(jīng)過(guò)多年工程化應(yīng)用發(fā)展,形成了一套成熟的洞壁干擾修正軟件,給定模型的參數(shù)信息、模型的氣動(dòng)力系數(shù)測(cè)量值及壁面壓力系數(shù)測(cè)量值,可得到各氣動(dòng)力系數(shù)的洞壁干擾修正結(jié)果和模型上的擾動(dòng)速度分布。其適用范圍為Ma<0.9,試驗(yàn)迎角范圍可達(dá)±90°,風(fēng)洞壁可以是實(shí)壁或通氣壁。模型參數(shù)默認(rèn)為翼身尾組合體模型,飛翼布局只能近似設(shè)置。

    2.2 全速勢(shì)位流方法

    全速勢(shì)位流方法的控制方程為:

    根據(jù)給定的模型、支桿等物面網(wǎng)格,自動(dòng)生成空間笛卡爾網(wǎng)格,數(shù)值求解控制方程與邊界層方程得到無(wú)洞壁條件下的流場(chǎng)信息。壁面是通過(guò)面源和面偶極子奇點(diǎn)分布的模擬,由無(wú)洞壁條件下求取的壁面附近的速度場(chǎng)和特定壁面邊界條件得到奇點(diǎn)強(qiáng)度,進(jìn)而得到洞壁對(duì)模型的氣動(dòng)特性的影響。由于其控制方程為全速勢(shì)方程,因此可應(yīng)用于更高速度的跨聲速流動(dòng),其Ma數(shù)應(yīng)用范圍為Ma<1.0。

    壁面區(qū)域的邊界條件使用擾動(dòng)速勢(shì) 表示,實(shí)壁時(shí):

    開(kāi)孔壁邊界條件為:

    其中R為透氣壁參數(shù),不同的透氣壁參數(shù)R模擬不同的開(kāi)閉比影響,R=0為實(shí)壁,無(wú)窮大時(shí)為自由流。R可通過(guò)由該方法計(jì)算得到的壁面壓力與壁壓條所測(cè)的壁面壓力結(jié)果的吻合程度來(lái)估計(jì)。

    圖6為飛翼標(biāo)模全速勢(shì)位流方法模擬用網(wǎng)格,圖7給出了不同的R參數(shù)壁壓計(jì)算結(jié)果與FL-24風(fēng)洞壁壓測(cè)量結(jié)果的對(duì)比曲線。可以看到,R=1時(shí)計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合最好,這與FL-24風(fēng)洞歷史壁特性參數(shù)測(cè)量結(jié)果是一致的,可用來(lái)模擬FL-24風(fēng)洞開(kāi)孔壁。

    圖6 全速勢(shì)位流方程方法用網(wǎng)格Fig.6 Gird for full potential equation method

    圖7 不同R參數(shù)壁壓計(jì)算結(jié)果與測(cè)量結(jié)果比較Fig.7 Comparison of variant Rparameter calculation results and measured result

    2.3 基于RANS的壁壓信息法

    基于RANS的壁壓信息法在自主研發(fā)的多塊結(jié)構(gòu)網(wǎng)格求解器基礎(chǔ)上融合了壁壓信息邊界條件處理方法,采用相同的控制方程和模型網(wǎng)格,分別以遠(yuǎn)場(chǎng)和試驗(yàn)測(cè)得壁壓為邊界條件,計(jì)算出模型自由流和洞壁存在時(shí)的表面力(力矩)系數(shù),從而得到洞壁干擾量。求解器采用有限體積法,控制方程為雷諾平均N-S方程組,湍流模型為兩方程k-wSST模型,為了加快收斂速度,采用了局部時(shí)間步長(zhǎng)和多重網(wǎng)格技術(shù)。圖8給出了數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果的對(duì)比曲線,可以看出,RANS求解器在升力曲線拐點(diǎn)前與試驗(yàn)結(jié)果吻合良好,拐點(diǎn)迎角比試驗(yàn)結(jié)果推后。

    圖8 RANS CFD結(jié)果與試驗(yàn)結(jié)果比較Fig.8 Comparison of RANS CFD and test results

    壁壓信息處理需要將壁壓條測(cè)得的壁壓插值到洞壁邊界計(jì)算網(wǎng)格點(diǎn)上,并將前后遠(yuǎn)場(chǎng)外插到0。洞壁邊界條件為:

    其中,下標(biāo)B表示最外層邊界網(wǎng)格,下標(biāo)“內(nèi)”表示內(nèi)層網(wǎng)格。

    從理論上講,由于這種方法的控制方程是非線性的,可以突破線性理論的限制,對(duì)于各種復(fù)雜流動(dòng)均可以進(jìn)行洞壁干擾修正。

    2.4 洞壁干擾修正結(jié)果分析

    采用上述三種方法分別對(duì)小展弦比飛翼標(biāo)模FL-24風(fēng)洞試驗(yàn)結(jié)果進(jìn)行洞壁干擾修正。圖9給出了三種修正方法洞壁干擾修正量與FL-26風(fēng)洞驗(yàn)證試驗(yàn)結(jié)果的對(duì)比曲線。

    可以看到,三種方法中,基于RANS的壁壓信息法與全速勢(shì)位流方法在線性段吻合良好,與試驗(yàn)結(jié)果趨勢(shì)一致,量值略有差異,原因是試驗(yàn)結(jié)果同樣存在誤差,從規(guī)律上看,試驗(yàn)結(jié)果曲線偏下。試驗(yàn)結(jié)果中等迎角趨勢(shì)反向問(wèn)題在基于RANS的壁壓信息法結(jié)果中同樣有所反映,迎角12°后量值與試驗(yàn)結(jié)果差異較大,這與RANS求解器在渦破裂迎角的捕捉上與試驗(yàn)結(jié)果差異較大有關(guān)(見(jiàn)圖8)。傳統(tǒng)的小擾動(dòng)位流壁壓信息法結(jié)果在趨勢(shì)上能夠保持一致,但量值差異稍大。這可能與該方法默認(rèn)的模型為翼身尾組合體,與飛翼布局存在差異導(dǎo)致。

    圖9 不同洞壁干擾修正方法結(jié)果比較Fig.9 Comparison of variant wall interference correction method results

    洞壁干擾在不同迎角范圍內(nèi)影響方式不同,整體看,存在三段不同影響區(qū)間。在迎角-2°~4°附著流為主要流態(tài)階段,洞壁干擾的影響表現(xiàn)為上(下)洗流的影響和Ma數(shù)增量的影響。在4°~12°左右,模型上表面逐漸形成穩(wěn)定的渦系,渦升力表現(xiàn)突出,升力呈現(xiàn)非線性增加,該階段洞壁影響不再是洗流和Ma數(shù)增量影響,而是影響了渦的空間位置和強(qiáng)度,使得洞壁干擾量呈現(xiàn)反向增加趨勢(shì)。在渦發(fā)生破裂后,洞壁影響更加復(fù)雜,試驗(yàn)結(jié)果存在諸多偽洞壁干擾因素,而計(jì)算結(jié)果難以捕捉渦破裂角度而存在較大誤差。

    圖10 不同迎角有無(wú)洞壁模擬表面壓力分布對(duì)比Fig.10 Comparison of surface pressure distribution on different angles of attack with and without wind tunnel wall simulation

    為了分析洞壁干擾影響在不同角度出現(xiàn)趨勢(shì)相反的問(wèn)題,圖10給出了基于RANS方法得到的迎角8°與12°、有無(wú)洞壁模擬表面壓力分布的對(duì)比云圖??梢钥闯?,在迎角8°,洞壁的存在使上下表面壓力系數(shù)負(fù)值增加,12°時(shí)洞壁的存在卻是使其負(fù)值減小。從8°與12°表面壓力差異看,8°僅部分展向位置受渦的影響出現(xiàn)吸力鼓包。而在12°時(shí),大部分展向位置出現(xiàn)壓力鼓包,也就是說(shuō),12°迎角渦升力顯著,此時(shí),洞壁干擾修正量出現(xiàn)趨勢(shì)反向,這在一定程度論證了洞壁干擾在不同迎角區(qū)間影響方式不同的說(shuō)法。

    2.5 洞壁干擾修正方法比較

    不同風(fēng)洞試驗(yàn)結(jié)果差異不僅僅表現(xiàn)在洞壁干擾上,來(lái)流湍流度、支撐干擾等都會(huì)有一定影響,因此,試驗(yàn)得到的洞壁干擾結(jié)果也并不能認(rèn)為是真值所在。綜合比較三種方法,可以得出:小擾動(dòng)位流壁壓信息法在量值上與其他結(jié)果相差較大,模型模擬方法需要針對(duì)飛翼布局進(jìn)行改進(jìn);全速勢(shì)位流方法在線性段能夠捕捉規(guī)律,且量值與基于RANS的壁壓信息法較為接近,計(jì)算模型構(gòu)造簡(jiǎn)單,計(jì)算速度較快,而且不需要壁壓信息;基于RANS的壁壓信息法在壁面無(wú)強(qiáng)激波存、氣流無(wú)大面積分離或者無(wú)渦破裂條件下,都能得到較為準(zhǔn)確的結(jié)果,但該方法需要壁壓信息,且空間網(wǎng)格質(zhì)量要求較高,網(wǎng)格準(zhǔn)備復(fù)雜,計(jì)算時(shí)間長(zhǎng)。

    表1 飛翼布局洞壁干擾修正方法對(duì)比Table 1 Comparison of variant wall interference correction method for flying wing configuration

    3 結(jié) 論

    本文采用了小擾動(dòng)壁壓信息法、全速勢(shì)位流方法和基于RANS的壁壓信息法針對(duì)小展弦比飛翼標(biāo)模FL-24風(fēng)洞試驗(yàn)結(jié)果開(kāi)展洞壁干擾修正計(jì)算,通過(guò)與試驗(yàn)結(jié)果的比較,綜合分析了三種方法的優(yōu)缺點(diǎn)。結(jié)果表明:

    (1)飛翼標(biāo)模洞壁干擾不同于翼-身-尾組合體布局,小擾動(dòng)壁壓信息法中模型的模擬方法需針對(duì)飛翼布局進(jìn)行改進(jìn)。

    (2)引入全速勢(shì)位流方法評(píng)估洞壁干擾,該方法具備較強(qiáng)的工程實(shí)用性。

    (3)基于RANS的壁壓信息法能在較大的范圍內(nèi)獲得良好結(jié)果。

    [1] 范召林,惲起麟,崔乃明.跨聲速洞壁干擾修正的壁壓法及其應(yīng)用[J],氣動(dòng)實(shí)驗(yàn)與測(cè)量控制,1990,4(4):45-51

    [2] 范召林,崔乃明,惲起麟.壁壓法應(yīng)用于高速風(fēng)洞大迎角試驗(yàn)洞壁干擾的修正[J],空氣動(dòng)力學(xué)學(xué)報(bào),1991,9(2):243-250.

    [3] 范召林,尹陸平.近聲速測(cè)壁壓試驗(yàn)研究[R].CN-930262,CARDC-2,1993.

    [4] 張玉倫,陳作斌,等.實(shí)用外形亞跨超流場(chǎng)數(shù)值模擬研究[J].空氣動(dòng)力學(xué)學(xué)報(bào),1995,13(1):57-62.

    [5] Fan Z L,Zhang Y L,He Z.Preliminary research on numerical simulation of transonic wall interference for aircraft model tests[J].Acta Aeronautica et Astronautica Sinica,1997,18(2):210-214.(in Chinese)范召林,張玉倫,賀中.跨聲速飛機(jī)模型試驗(yàn)洞壁干擾數(shù)值模擬的初步研究[J].航空學(xué)報(bào),1997,18(2):210-214.

    [6] Fan Z L,Zhang Y L,He Z,et al.Numerical calculation for three dimensional non-linear transonic wall interference[J].Acta Aerodynamica Sinica,1997,15(2):145-152.(in Chinese)范召林,張玉倫,賀中,等.跨聲速三維非線性洞壁干擾的數(shù)值計(jì)算[J].空氣動(dòng)力學(xué)學(xué)報(bào),1997,15(2):145-152.

    [7] Yin L P,He Z,Yu Z S,F(xiàn)ang Z L,Chen Z B.Research on subsonic wall interference correction for model tests at high angle of attack[J].Experiments and Measurements in Fluid Mechanics,2000,14(3):37-41.(in Chinese)尹陸平,賀中,于志松,范召林,陳作斌.亞聲速大迎角模型試驗(yàn)洞壁干擾修正方法研究[J].流體力學(xué)實(shí)驗(yàn)與測(cè)量,2000,14(3):37-41.

    [8] Allan M R,Badcock K J.A RANS investigation of wind tunnel interference effects on delta wing aerodynamics[R].AIAA 2003-4214.

    [9] Masataka Kohzai,Makoto Ueno.Wall and Support Interference Corrections of NASA Common Research Model Wind Tunnel Tests in JAXA[R].AIAA 2013-0963.

    [10]Maseland J E J,Laban M.Development of CFD-based interference models for DNW-HST transonic wind tunnel[R].AIAA 2006-3639

    [11]Robert C Griffiths,Kenneth D Visser.Panel code estimates of wind tunnel wall and support interference effects on the HSCT[R].AIAA 2000-4020.

    [12]Alexander J Krynytzky.Parametric model size study of wall interference in the boeing transonic wind tunnel using TRANAIR[R].AIAA 2004-2310.

    [13]Alexander J Krynytzky.Uncertainty evaluation of wall interference in a large transonic wind tunnel[R].AIAA 2010-4341.

    [14]Li Y H,Huang Y,Wang Y Q,et al.Wing-body thickness effects on aerodynamic and vortex flow characteristics of a common low aspect ratio flying-wing configuration at transonic flow[J].Acta Aerodynamica Sinica 2015,33(3):302-306.(in Chinese)李永紅,黃勇,王義慶,等.翼身相對(duì)厚度對(duì)小展弦比飛翼布局跨聲速氣動(dòng)特性及流動(dòng)機(jī)理的影響研究[J].空氣動(dòng)力學(xué)學(xué)報(bào),2015,33(3):302-306.

    [15]Su J C,Huang Y,Li Y H,et al.Support interference of low aspect ratio flying-wing from subsonic to supersonic speed[J].Acta Aerodynamica Sinica,2015,33(3):289-295.(in Chinese)蘇繼川,黃勇,李永紅,等.小展弦比飛翼亞、跨、超聲速支撐干擾研究[J].空氣動(dòng)力學(xué)學(xué)報(bào),2015,33(3):289-295.

    [16]Su J C,Huang Y,Zhong S D,et al.Research on flow characteristics of low-aspect-ratio flying-wing at transonic speed[J].Acta Aerodynamica Sinica,2015,33(3):307-312.(in Chinese)蘇繼川,黃勇,鐘世東,等.小展弦比飛翼跨聲速典型流動(dòng)特性研究[J].空氣動(dòng)力學(xué)學(xué)報(bào),2015,33(3):307-312.

    Three types of transonic wind tunnel wall interference correction method and their application on low aspect ratio flying wing calibration model

    Zhong Shidong*,Li Wei,Su Jichuan,Li Yonghong,He Zhong

    (High Speed Aerodynamic Institute of China Aerodynamics Research and Development Center,Mianyang Sichuan 621000,China)

    Force test on wall pressure measurement has been conducted in CARDC FL-24 wind tunnel,and the corresponding verification test has been carried out in CARDC FL-26wind Tunnel for a flying wing calibration model with low aspect ratio.Three types of transonic wind tunnel wall interference correction method,i.e.,small perturbation wall pressure signature correction method,full potential equation method and RANS based wall pressure signature correction method,are adopted to calculate wall interference.The calculation results are compared with test result.It shows that the wall interference tendency of low aspect ratio flying wing calibration model is different from that of a wing-body-tail configuration.The three methods have their individual applicability when they are used for flying wing configuration wall interference correction.

    flying wing calibration model;wind tunnel wall interference;numerical simulation

    V211.71

    A

    10.7638/kqdlxxb-2015.0095

    0258-1825(2016)01-0113-06

    2015-07-21;

    2015-09-22

    鐘世東*(1983-),男,山東青州人,工程師,研究方向:氣動(dòng)布局設(shè)計(jì).E-mail:zhongshidong1983@163.com

    鐘世東,李巍,蘇繼川,等.三種跨聲速洞壁干擾修正方法及其在小展弦比飛翼標(biāo)模試驗(yàn)中的應(yīng)用[J].空氣動(dòng)力學(xué)學(xué)報(bào),2016,34(1):113-118.

    10.7638/kqdlxxb-2015.0095 Zhong S D,Li W,Su J C,et al.Three types of transonic wind tunnel wall interference correction method and their application on low aspect ratio flying wing calibration model[J].Acta Aerodynamica Sinica,2016,34(1):113-118.

    猜你喜歡
    洞壁飛翼迎角
    “洞壁”聽(tīng)“琴音”(3)
    連續(xù)變迎角試驗(yàn)數(shù)據(jù)自適應(yīng)分段擬合濾波方法
    側(cè)壓影響下圓形洞室?guī)r爆雙軸物理模擬試驗(yàn)研究
    創(chuàng)新
    飛翼無(wú)人機(jī)機(jī)動(dòng)飛行非線性魯棒控制方法
    飛翼布局飛機(jī)側(cè)風(fēng)起降特性
    小展弦比飛翼標(biāo)模FL-2風(fēng)洞跨聲速開(kāi)孔壁干擾特性修正研究
    飛翼無(wú)人機(jī)嗡鳴氣動(dòng)彈性響應(yīng)分析
    飛翼隱身特性數(shù)值模擬
    失速保護(hù)系統(tǒng)迎角零向跳變研究
    科技傳播(2014年4期)2014-12-02 01:59:42
    午夜福利乱码中文字幕| 成人亚洲精品一区在线观看| 最新的欧美精品一区二区| 一本—道久久a久久精品蜜桃钙片| 免费一级毛片在线播放高清视频 | 亚洲欧美色中文字幕在线| 久久久久精品人妻al黑| 一区二区三区精品91| 午夜成年电影在线免费观看| 久久精品国产综合久久久| 久久久久久久大尺度免费视频| 国产主播在线观看一区二区| 亚洲精品自拍成人| av一本久久久久| 色婷婷久久久亚洲欧美| 国产日韩欧美亚洲二区| 国产三级黄色录像| 国产国语露脸激情在线看| 久久99热这里只频精品6学生| 一边摸一边做爽爽视频免费| 人人妻人人爽人人添夜夜欢视频| 十八禁人妻一区二区| 免费av中文字幕在线| netflix在线观看网站| 久久久久久久大尺度免费视频| 中文字幕人妻丝袜一区二区| 亚洲天堂av无毛| 两个人看的免费小视频| 日韩中文字幕欧美一区二区| 99国产精品免费福利视频| 大香蕉久久成人网| 黄色a级毛片大全视频| 九色亚洲精品在线播放| 他把我摸到了高潮在线观看 | 日本猛色少妇xxxxx猛交久久| 国产精品久久久久成人av| 97在线人人人人妻| 欧美xxⅹ黑人| 亚洲欧洲日产国产| 日韩 亚洲 欧美在线| 精品一品国产午夜福利视频| a在线观看视频网站| 国产高清videossex| 狂野欧美激情性bbbbbb| 深夜精品福利| 日本一区二区免费在线视频| 人妻人人澡人人爽人人| 日韩电影二区| 亚洲人成电影观看| 中文字幕av电影在线播放| 免费在线观看影片大全网站| 国内毛片毛片毛片毛片毛片| 国产熟女午夜一区二区三区| 黑人操中国人逼视频| 999久久久精品免费观看国产| 亚洲,欧美精品.| 国产黄色免费在线视频| 精品少妇久久久久久888优播| 美女视频免费永久观看网站| 久久国产精品大桥未久av| 日本vs欧美在线观看视频| 国产成人精品久久二区二区免费| 国产在线观看jvid| a在线观看视频网站| 久久热在线av| 国产成人欧美| 国产男女超爽视频在线观看| 欧美一级毛片孕妇| 女性被躁到高潮视频| 久久国产精品男人的天堂亚洲| 久久国产精品大桥未久av| 99热网站在线观看| 九色亚洲精品在线播放| 黄片小视频在线播放| 在线观看免费日韩欧美大片| 99国产精品99久久久久| 波多野结衣av一区二区av| 国产精品av久久久久免费| 日韩大片免费观看网站| 欧美国产精品va在线观看不卡| 日本欧美视频一区| 亚洲黑人精品在线| av天堂久久9| 丰满饥渴人妻一区二区三| 99精品欧美一区二区三区四区| 免费高清在线观看视频在线观看| 中文字幕av电影在线播放| 日韩制服丝袜自拍偷拍| 男女下面插进去视频免费观看| 多毛熟女@视频| 少妇 在线观看| 亚洲男人天堂网一区| 一区在线观看完整版| 日日夜夜操网爽| 999久久久国产精品视频| 五月开心婷婷网| 亚洲国产av新网站| 午夜福利影视在线免费观看| 午夜免费观看性视频| 丝袜人妻中文字幕| 热re99久久国产66热| 十八禁高潮呻吟视频| 中文字幕人妻丝袜制服| 欧美国产精品一级二级三级| 亚洲精品成人av观看孕妇| a级片在线免费高清观看视频| 青草久久国产| 国产成人免费无遮挡视频| 国产精品熟女久久久久浪| av在线app专区| 美女脱内裤让男人舔精品视频| e午夜精品久久久久久久| 国产精品 欧美亚洲| 国产精品熟女久久久久浪| 自线自在国产av| 亚洲成人国产一区在线观看| 极品人妻少妇av视频| 欧美 亚洲 国产 日韩一| 精品国产一区二区三区四区第35| 亚洲人成电影观看| 精品国产国语对白av| 亚洲欧美色中文字幕在线| 电影成人av| 91大片在线观看| 黄色 视频免费看| 国产成人精品在线电影| 婷婷丁香在线五月| 亚洲精品中文字幕在线视频| 肉色欧美久久久久久久蜜桃| 制服人妻中文乱码| av线在线观看网站| 久久免费观看电影| 一级毛片精品| 亚洲精品久久成人aⅴ小说| 最新的欧美精品一区二区| 一级毛片女人18水好多| 午夜福利视频在线观看免费| 最新的欧美精品一区二区| 啦啦啦中文免费视频观看日本| av线在线观看网站| 这个男人来自地球电影免费观看| 精品高清国产在线一区| 亚洲精品久久成人aⅴ小说| 久久性视频一级片| 美国免费a级毛片| 精品一品国产午夜福利视频| 国产欧美日韩一区二区精品| av不卡在线播放| 久久久精品免费免费高清| 大香蕉久久成人网| 老司机深夜福利视频在线观看 | 久久 成人 亚洲| 日本91视频免费播放| 精品福利永久在线观看| 亚洲精品av麻豆狂野| av在线app专区| 伊人亚洲综合成人网| 两个人看的免费小视频| 91成人精品电影| 男女边摸边吃奶| 久久99热这里只频精品6学生| 极品少妇高潮喷水抽搐| 免费人妻精品一区二区三区视频| 亚洲自偷自拍图片 自拍| 国产精品国产av在线观看| 亚洲天堂av无毛| 一区二区三区乱码不卡18| 午夜免费观看性视频| 69精品国产乱码久久久| 美女国产高潮福利片在线看| 性高湖久久久久久久久免费观看| 多毛熟女@视频| 欧美 日韩 精品 国产| 中文字幕人妻丝袜制服| 窝窝影院91人妻| 国产精品久久久av美女十八| 波多野结衣一区麻豆| 在线观看舔阴道视频| 国产精品香港三级国产av潘金莲| 一进一出抽搐动态| 精品少妇一区二区三区视频日本电影| 国产男女超爽视频在线观看| 欧美久久黑人一区二区| 午夜影院在线不卡| 日本一区二区免费在线视频| 美女国产高潮福利片在线看| 色94色欧美一区二区| 午夜久久久在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 永久免费av网站大全| 久久精品国产亚洲av香蕉五月 | 国产精品一二三区在线看| 十八禁高潮呻吟视频| 国产在线视频一区二区| 男女边摸边吃奶| 老司机福利观看| 99久久99久久久精品蜜桃| 国产成人欧美| 在线观看一区二区三区激情| 国产不卡av网站在线观看| 99国产精品一区二区蜜桃av | 十八禁网站免费在线| 这个男人来自地球电影免费观看| 国产精品欧美亚洲77777| 色老头精品视频在线观看| 午夜影院在线不卡| 51午夜福利影视在线观看| 日韩中文字幕视频在线看片| 中文字幕精品免费在线观看视频| 日本91视频免费播放| 51午夜福利影视在线观看| 大码成人一级视频| 91九色精品人成在线观看| 91麻豆精品激情在线观看国产 | 80岁老熟妇乱子伦牲交| 18禁裸乳无遮挡动漫免费视频| 97人妻天天添夜夜摸| 黄色视频不卡| 搡老乐熟女国产| 久久久国产欧美日韩av| 99精品久久久久人妻精品| 欧美国产精品va在线观看不卡| www.999成人在线观看| 欧美一级毛片孕妇| 日本五十路高清| 桃红色精品国产亚洲av| 亚洲国产精品成人久久小说| 成在线人永久免费视频| 免费在线观看日本一区| 国产欧美日韩一区二区三 | 老司机靠b影院| 男女下面插进去视频免费观看| 窝窝影院91人妻| 日韩欧美一区视频在线观看| 亚洲第一欧美日韩一区二区三区 | 搡老岳熟女国产| 国产欧美日韩一区二区精品| 99热全是精品| 69av精品久久久久久 | 成人国产一区最新在线观看| 制服诱惑二区| 精品福利观看| 亚洲九九香蕉| 看免费av毛片| 免费黄频网站在线观看国产| 国产成人精品久久二区二区91| 国产精品一区二区免费欧美 | 天堂8中文在线网| 男人添女人高潮全过程视频| 成人国语在线视频| 成人av一区二区三区在线看 | 天天躁夜夜躁狠狠躁躁| 中文欧美无线码| 国产精品熟女久久久久浪| 国产三级黄色录像| 亚洲情色 制服丝袜| 国产一区有黄有色的免费视频| 亚洲熟女精品中文字幕| 久久人人爽人人片av| 国产有黄有色有爽视频| 日日摸夜夜添夜夜添小说| 他把我摸到了高潮在线观看 | 岛国在线观看网站| 国产免费现黄频在线看| 国产99久久九九免费精品| av视频免费观看在线观看| 亚洲av男天堂| 亚洲av成人不卡在线观看播放网 | 涩涩av久久男人的天堂| 久久这里只有精品19| 亚洲中文字幕日韩| 精品免费久久久久久久清纯 | 黄色a级毛片大全视频| 国产淫语在线视频| 91麻豆av在线| 欧美av亚洲av综合av国产av| 少妇人妻久久综合中文| 久久人人97超碰香蕉20202| 黄片大片在线免费观看| 一边摸一边做爽爽视频免费| 亚洲激情五月婷婷啪啪| 亚洲精品一卡2卡三卡4卡5卡 | 99久久综合免费| 亚洲第一青青草原| 亚洲伊人久久精品综合| 欧美激情久久久久久爽电影 | 一个人免费在线观看的高清视频 | 亚洲一区中文字幕在线| 老汉色av国产亚洲站长工具| 日本wwww免费看| 亚洲欧美清纯卡通| 欧美激情 高清一区二区三区| 老司机午夜十八禁免费视频| 天堂俺去俺来也www色官网| 一本一本久久a久久精品综合妖精| 亚洲色图 男人天堂 中文字幕| 一本综合久久免费| 亚洲视频免费观看视频| 亚洲国产av新网站| www.熟女人妻精品国产| 建设人人有责人人尽责人人享有的| 精品国产一区二区三区久久久樱花| av一本久久久久| 老汉色∧v一级毛片| 亚洲精品国产精品久久久不卡| 国产97色在线日韩免费| 色老头精品视频在线观看| 久久av网站| 国产精品 国内视频| 黑人猛操日本美女一级片| 亚洲精品国产av蜜桃| 国产精品九九99| 久久久久国产一级毛片高清牌| 黄网站色视频无遮挡免费观看| a级片在线免费高清观看视频| 久久这里只有精品19| 美女国产高潮福利片在线看| 亚洲国产毛片av蜜桃av| 满18在线观看网站| 美女扒开内裤让男人捅视频| 午夜福利影视在线免费观看| 12—13女人毛片做爰片一| 欧美午夜高清在线| 伊人亚洲综合成人网| 777米奇影视久久| 国产无遮挡羞羞视频在线观看| 在线精品无人区一区二区三| 国产野战对白在线观看| 久久久国产精品麻豆| 狂野欧美激情性xxxx| 婷婷色av中文字幕| 国产免费福利视频在线观看| 久久人人爽av亚洲精品天堂| 国产精品国产三级国产专区5o| 岛国毛片在线播放| 大片电影免费在线观看免费| 激情视频va一区二区三区| h视频一区二区三区| 亚洲美女黄色视频免费看| 午夜福利视频精品| 午夜两性在线视频| 99精国产麻豆久久婷婷| 99久久国产精品久久久| 精品一区在线观看国产| 精品国产国语对白av| 国产精品国产三级国产专区5o| 午夜福利视频在线观看免费| www.999成人在线观看| 久久精品国产亚洲av高清一级| 亚洲av男天堂| 老司机靠b影院| 久久精品国产亚洲av高清一级| 人人澡人人妻人| 美女午夜性视频免费| 波多野结衣一区麻豆| 午夜激情av网站| 精品国产乱码久久久久久男人| 天堂8中文在线网| avwww免费| 麻豆国产av国片精品| 丝袜脚勾引网站| 国产精品自产拍在线观看55亚洲 | 国产日韩欧美在线精品| 国产av又大| 丰满少妇做爰视频| 国产欧美亚洲国产| 久久久精品94久久精品| 国产精品久久久久久精品古装| 考比视频在线观看| 操出白浆在线播放| 亚洲欧美清纯卡通| 国产精品久久久人人做人人爽| 操出白浆在线播放| 国产视频一区二区在线看| av有码第一页| 天堂8中文在线网| 亚洲国产毛片av蜜桃av| 亚洲精品自拍成人| 欧美日韩国产mv在线观看视频| 亚洲欧美色中文字幕在线| 久久99一区二区三区| 久久ye,这里只有精品| 热re99久久精品国产66热6| 国产男女超爽视频在线观看| 18禁黄网站禁片午夜丰满| 欧美另类一区| 国产深夜福利视频在线观看| 男女之事视频高清在线观看| 少妇的丰满在线观看| 免费日韩欧美在线观看| 成年动漫av网址| 老鸭窝网址在线观看| 成人影院久久| 啦啦啦啦在线视频资源| 黄色视频不卡| 99久久精品国产亚洲精品| 亚洲精品国产av成人精品| 久久国产精品男人的天堂亚洲| 69av精品久久久久久 | 久久久久国产一级毛片高清牌| 亚洲va日本ⅴa欧美va伊人久久 | 日韩欧美免费精品| 国产成人啪精品午夜网站| 欧美 日韩 精品 国产| 色精品久久人妻99蜜桃| 人人妻人人添人人爽欧美一区卜| 国产精品麻豆人妻色哟哟久久| 汤姆久久久久久久影院中文字幕| 男人爽女人下面视频在线观看| 成人免费观看视频高清| 高清av免费在线| 国产成人欧美| 69av精品久久久久久 | 午夜福利影视在线免费观看| 热99久久久久精品小说推荐| 国产成人av教育| 久久人人爽人人片av| 90打野战视频偷拍视频| 成人影院久久| 天天躁狠狠躁夜夜躁狠狠躁| 免费观看a级毛片全部| 午夜激情av网站| 精品国产一区二区三区久久久樱花| 我的亚洲天堂| 国产主播在线观看一区二区| 国产一区二区在线观看av| 动漫黄色视频在线观看| 永久免费av网站大全| 亚洲欧美精品综合一区二区三区| 后天国语完整版免费观看| 国产欧美日韩一区二区三区在线| 欧美黄色淫秽网站| 欧美变态另类bdsm刘玥| 91老司机精品| videosex国产| 一区二区三区乱码不卡18| 国产成人精品久久二区二区91| 侵犯人妻中文字幕一二三四区| 国产精品偷伦视频观看了| 国产欧美日韩一区二区三区在线| 亚洲精品乱久久久久久| 久久天堂一区二区三区四区| 啦啦啦视频在线资源免费观看| 美女脱内裤让男人舔精品视频| 欧美 亚洲 国产 日韩一| 欧美精品人与动牲交sv欧美| 精品国内亚洲2022精品成人 | 欧美 亚洲 国产 日韩一| 母亲3免费完整高清在线观看| 天堂俺去俺来也www色官网| 日韩有码中文字幕| 日本欧美视频一区| 制服人妻中文乱码| av有码第一页| 老司机影院成人| 免费看十八禁软件| 男人舔女人的私密视频| 爱豆传媒免费全集在线观看| 国产精品国产av在线观看| 久久99热这里只频精品6学生| 一个人免费在线观看的高清视频 | 久久国产精品大桥未久av| 免费日韩欧美在线观看| √禁漫天堂资源中文www| 欧美老熟妇乱子伦牲交| 国产1区2区3区精品| 国内毛片毛片毛片毛片毛片| 一级黄色大片毛片| 中文字幕人妻熟女乱码| 老熟女久久久| 人妻久久中文字幕网| 久久天堂一区二区三区四区| 国产成人免费无遮挡视频| 最近中文字幕2019免费版| 欧美激情久久久久久爽电影 | 中文字幕最新亚洲高清| av有码第一页| 一边摸一边抽搐一进一出视频| 91国产中文字幕| 99久久人妻综合| 欧美另类亚洲清纯唯美| 99精品欧美一区二区三区四区| 亚洲av成人一区二区三| 国产精品久久久久久人妻精品电影 | 法律面前人人平等表现在哪些方面 | 午夜激情久久久久久久| 免费在线观看完整版高清| 男男h啪啪无遮挡| 亚洲欧美一区二区三区黑人| 韩国高清视频一区二区三区| 操美女的视频在线观看| 久久天堂一区二区三区四区| 大片免费播放器 马上看| 精品福利永久在线观看| 亚洲国产成人一精品久久久| 久久久久视频综合| 午夜福利,免费看| 久久精品国产综合久久久| 一级片免费观看大全| 欧美亚洲日本最大视频资源| 亚洲精品乱久久久久久| 亚洲av美国av| 亚洲欧美清纯卡通| 91国产中文字幕| 国产日韩欧美亚洲二区| av天堂久久9| 久久精品人人爽人人爽视色| 老司机在亚洲福利影院| 在线观看免费午夜福利视频| 欧美97在线视频| 久久久久久久精品精品| 日本a在线网址| 一级黄色大片毛片| 国产一区二区在线观看av| 久久久精品免费免费高清| 丁香六月欧美| 电影成人av| 天天躁夜夜躁狠狠躁躁| 国产亚洲午夜精品一区二区久久| 国产精品.久久久| 老熟妇乱子伦视频在线观看 | 又紧又爽又黄一区二区| 人人妻人人澡人人看| 美女中出高潮动态图| 亚洲精品久久午夜乱码| 国产伦人伦偷精品视频| 动漫黄色视频在线观看| 可以免费在线观看a视频的电影网站| 欧美黄色片欧美黄色片| 精品卡一卡二卡四卡免费| 中国美女看黄片| 国产精品影院久久| 国产高清国产精品国产三级| 精品国产一区二区三区久久久樱花| 久久人人爽人人片av| 99热网站在线观看| 久久av网站| 欧美 亚洲 国产 日韩一| 亚洲美女黄色视频免费看| 国产一区有黄有色的免费视频| 国产成人av激情在线播放| 国产亚洲精品久久久久5区| 韩国高清视频一区二区三区| 各种免费的搞黄视频| 在线精品无人区一区二区三| 女性被躁到高潮视频| 亚洲熟女精品中文字幕| 美女中出高潮动态图| 日日摸夜夜添夜夜添小说| 一级片免费观看大全| 日韩欧美国产一区二区入口| 搡老乐熟女国产| 叶爱在线成人免费视频播放| 另类精品久久| 久久久国产一区二区| 女性生殖器流出的白浆| 亚洲精品一卡2卡三卡4卡5卡 | 12—13女人毛片做爰片一| 久久久国产一区二区| 国产精品.久久久| 久热爱精品视频在线9| 操美女的视频在线观看| 美女高潮喷水抽搐中文字幕| 老司机影院成人| 18禁黄网站禁片午夜丰满| 国产精品久久久人人做人人爽| 在线观看舔阴道视频| 欧美精品一区二区免费开放| 日本五十路高清| 天堂8中文在线网| 国产成人精品无人区| 亚洲精品国产av蜜桃| www.熟女人妻精品国产| 香蕉丝袜av| 欧美亚洲日本最大视频资源| 日本猛色少妇xxxxx猛交久久| 久久精品人人爽人人爽视色| 亚洲欧美激情在线| 一本色道久久久久久精品综合| 99re6热这里在线精品视频| 岛国在线观看网站| 在线观看免费午夜福利视频| 国产又色又爽无遮挡免| 亚洲自偷自拍图片 自拍| 国产熟女午夜一区二区三区| 久久久精品免费免费高清| 一区二区三区乱码不卡18| 亚洲天堂av无毛| 成年动漫av网址| 久久精品国产综合久久久| 如日韩欧美国产精品一区二区三区| 精品视频人人做人人爽| 午夜日韩欧美国产| 香蕉国产在线看| 国产主播在线观看一区二区| 91九色精品人成在线观看| 成人黄色视频免费在线看| 精品人妻在线不人妻| 国产成+人综合+亚洲专区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品久久久人人做人人爽| 又黄又粗又硬又大视频| 免费在线观看黄色视频的| 久久中文看片网| 亚洲一区中文字幕在线| 少妇人妻久久综合中文| 亚洲av美国av| 欧美日韩精品网址| 成人国产av品久久久| 叶爱在线成人免费视频播放| 飞空精品影院首页| 国产人伦9x9x在线观看| 在线看a的网站| 乱人伦中国视频| 一级毛片女人18水好多|