• <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在线观看视频网站| 美女扒开内裤让男人捅视频| 两性午夜刺激爽爽歪歪视频在线观看 | 999精品在线视频| 丰满的人妻完整版| 交换朋友夫妻互换小说| 欧美日韩亚洲国产一区二区在线观看| 女性被躁到高潮视频| 国产av在哪里看| 亚洲中文av在线| 一边摸一边抽搐一进一小说| 亚洲精品国产一区二区精华液| 国产免费男女视频| 国产真人三级小视频在线观看| 美女高潮到喷水免费观看| 夜夜夜夜夜久久久久| 视频区图区小说| 亚洲欧美激情综合另类| 久久精品亚洲av国产电影网| 中国美女看黄片| 国产黄色免费在线视频| 亚洲五月天丁香| 久久天躁狠狠躁夜夜2o2o| 久久久久国内视频| 99热国产这里只有精品6| 一区二区三区精品91| 美女国产高潮福利片在线看| 欧美大码av| 久久伊人香网站| 极品人妻少妇av视频| 免费在线观看黄色视频的| 在线观看日韩欧美| bbb黄色大片| 夜夜夜夜夜久久久久| 自拍欧美九色日韩亚洲蝌蚪91| 校园春色视频在线观看| 最近最新中文字幕大全电影3 | 欧美亚洲日本最大视频资源| 又紧又爽又黄一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 日日干狠狠操夜夜爽| 少妇的丰满在线观看| 亚洲av片天天在线观看| 国产在线观看jvid| a级毛片在线看网站| 美女扒开内裤让男人捅视频| 亚洲成人免费电影在线观看| 久久久水蜜桃国产精品网| 成人手机av| 国产日韩一区二区三区精品不卡| 动漫黄色视频在线观看| 99热只有精品国产| 精品熟女少妇八av免费久了| 91在线观看av| 国产欧美日韩一区二区三区在线| av有码第一页| 国产97色在线日韩免费| 色老头精品视频在线观看| 男人舔女人下体高潮全视频| 麻豆av在线久日| 国产精品永久免费网站| 一边摸一边做爽爽视频免费| 欧美激情极品国产一区二区三区| 欧美日韩福利视频一区二区| 色播在线永久视频| 美女午夜性视频免费| 亚洲国产精品一区二区三区在线| 天堂影院成人在线观看| 久久久久久人人人人人| 国产精品av久久久久免费| 国产区一区二久久| 成年版毛片免费区| 国产精品一区二区三区四区久久 | 亚洲第一青青草原| 国产精品美女特级片免费视频播放器 | 黄色a级毛片大全视频| 国产无遮挡羞羞视频在线观看| 国产精品二区激情视频| 亚洲人成伊人成综合网2020| 激情在线观看视频在线高清| 18禁国产床啪视频网站| 午夜精品在线福利| 香蕉久久夜色| 丰满的人妻完整版| 国产国语露脸激情在线看| 久久久久久久久中文| 人人妻,人人澡人人爽秒播| 极品教师在线免费播放| 日韩欧美在线二视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲激情在线av| 久久性视频一级片| 午夜91福利影院| 村上凉子中文字幕在线| 男女下面插进去视频免费观看| 亚洲精品国产一区二区精华液| 久久国产精品影院| 亚洲精品国产精品久久久不卡| 国产成人一区二区三区免费视频网站| 色老头精品视频在线观看| av视频免费观看在线观看| 日韩大码丰满熟妇| 热re99久久国产66热| av中文乱码字幕在线| 国产伦人伦偷精品视频| 日韩欧美在线二视频| 黄网站色视频无遮挡免费观看| 中文字幕人妻丝袜一区二区| 欧美丝袜亚洲另类 | 午夜视频精品福利| 久久人妻熟女aⅴ| 欧美成人免费av一区二区三区| 亚洲色图av天堂| 国产蜜桃级精品一区二区三区| 电影成人av| 久久午夜综合久久蜜桃| 亚洲成人久久性| 午夜福利在线免费观看网站| 国内久久婷婷六月综合欲色啪| 午夜两性在线视频| 久久久精品欧美日韩精品| 久久精品国产99精品国产亚洲性色 | av国产精品久久久久影院| 97碰自拍视频| 老鸭窝网址在线观看| 国产人伦9x9x在线观看| 999久久久国产精品视频| 天天添夜夜摸| 国产av一区在线观看免费| 999久久久精品免费观看国产| 国产三级黄色录像| av网站在线播放免费| 国产亚洲av高清不卡| 黄色视频不卡| 在线av久久热| 少妇 在线观看| 国产一区二区三区综合在线观看| 精品午夜福利视频在线观看一区| 国产无遮挡羞羞视频在线观看| 90打野战视频偷拍视频| 成年女人毛片免费观看观看9| av免费在线观看网站| 欧美黄色片欧美黄色片| 久久草成人影院| 亚洲va日本ⅴa欧美va伊人久久| 黄色女人牲交| 免费在线观看亚洲国产| 亚洲欧美激情在线| 国产免费现黄频在线看| 久久久久精品国产欧美久久久| 午夜视频精品福利| 大型av网站在线播放| 久久人妻av系列| 久久香蕉精品热| av视频免费观看在线观看| 欧美激情高清一区二区三区| 黄频高清免费视频| 丝袜人妻中文字幕| 久久久国产成人精品二区 | 91成年电影在线观看| 亚洲aⅴ乱码一区二区在线播放 | 国产精品亚洲一级av第二区| 19禁男女啪啪无遮挡网站| aaaaa片日本免费| 老司机亚洲免费影院| 国产三级在线视频| 亚洲av第一区精品v没综合| www.熟女人妻精品国产| 91成人精品电影| 人成视频在线观看免费观看| 真人做人爱边吃奶动态| 亚洲精品一二三| 亚洲人成伊人成综合网2020| 亚洲精品在线美女| 男女之事视频高清在线观看| 脱女人内裤的视频| 男女做爰动态图高潮gif福利片 | 99久久国产精品久久久| 国产精品98久久久久久宅男小说| 19禁男女啪啪无遮挡网站| 久久草成人影院| 热re99久久精品国产66热6| 日韩视频一区二区在线观看| 级片在线观看| 精品人妻1区二区| 久久久久久久久久久久大奶| 麻豆一二三区av精品| 午夜福利欧美成人| 性少妇av在线| 国产亚洲精品久久久久久毛片| 亚洲久久久国产精品| 国产亚洲欧美98| 欧美 亚洲 国产 日韩一| 天天躁狠狠躁夜夜躁狠狠躁| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲精品久久午夜乱码| 美女国产高潮福利片在线看| 岛国在线观看网站| 成人18禁高潮啪啪吃奶动态图| 又紧又爽又黄一区二区| av中文乱码字幕在线| 天堂中文最新版在线下载| 国产欧美日韩综合在线一区二区| 美女福利国产在线| 国产视频一区二区在线看| 国产精品久久久人人做人人爽| av超薄肉色丝袜交足视频| 成人18禁在线播放| 水蜜桃什么品种好| 天堂动漫精品| 一区在线观看完整版| 亚洲精品国产色婷婷电影| 水蜜桃什么品种好| 国产主播在线观看一区二区| 久久精品91无色码中文字幕| 法律面前人人平等表现在哪些方面| 欧美激情久久久久久爽电影 | 激情在线观看视频在线高清| 欧美日本亚洲视频在线播放| av在线天堂中文字幕 | 无限看片的www在线观看| 日本五十路高清| 国产一区二区三区在线臀色熟女 | 久久人人精品亚洲av| 国产一卡二卡三卡精品| 这个男人来自地球电影免费观看| 淫妇啪啪啪对白视频| 国产欧美日韩一区二区三区在线| 亚洲九九香蕉| 美女扒开内裤让男人捅视频| 不卡一级毛片| 亚洲精品美女久久久久99蜜臀| 香蕉国产在线看| 欧美乱码精品一区二区三区| 国产又爽黄色视频| 免费搜索国产男女视频| 多毛熟女@视频| 水蜜桃什么品种好| 长腿黑丝高跟| 黄色片一级片一级黄色片| 又黄又爽又免费观看的视频| 人人妻人人爽人人添夜夜欢视频| 久久性视频一级片| 日日摸夜夜添夜夜添小说| 久久国产乱子伦精品免费另类| 交换朋友夫妻互换小说| 久久影院123| 激情视频va一区二区三区| 欧美日韩中文字幕国产精品一区二区三区 | 好男人电影高清在线观看| 少妇裸体淫交视频免费看高清 | 丰满人妻熟妇乱又伦精品不卡| 夜夜看夜夜爽夜夜摸 | 99在线视频只有这里精品首页| 亚洲中文日韩欧美视频| 涩涩av久久男人的天堂| 色精品久久人妻99蜜桃| 一级毛片高清免费大全| 黑人巨大精品欧美一区二区蜜桃| 亚洲精品中文字幕一二三四区| 老司机福利观看| 日韩免费高清中文字幕av| 乱人伦中国视频| 国产亚洲精品久久久久5区| 免费av中文字幕在线| 亚洲成人免费电影在线观看| 19禁男女啪啪无遮挡网站| 在线av久久热| 亚洲成人免费电影在线观看| 国产黄a三级三级三级人| 国产成人av教育| 女性被躁到高潮视频| av视频免费观看在线观看| 母亲3免费完整高清在线观看| 成人影院久久| 老司机靠b影院| 国产黄a三级三级三级人| 欧美成人免费av一区二区三区| 亚洲一区二区三区色噜噜 | 一级,二级,三级黄色视频| 中文字幕最新亚洲高清| 久久国产精品影院| 久久天堂一区二区三区四区| 久久九九热精品免费| 黑人操中国人逼视频| 一区二区日韩欧美中文字幕| 视频在线观看一区二区三区| 久久久久久久久中文| 丁香欧美五月| 黄片大片在线免费观看| 精品高清国产在线一区| 国产又色又爽无遮挡免费看| 嫁个100分男人电影在线观看| 丰满人妻熟妇乱又伦精品不卡| 99久久久亚洲精品蜜臀av| 国产黄色免费在线视频| 精品一区二区三区av网在线观看| 国产成人欧美| 亚洲九九香蕉| 夫妻午夜视频| 色哟哟哟哟哟哟| 成年人免费黄色播放视频| aaaaa片日本免费| 国产成人免费无遮挡视频| 亚洲国产欧美日韩在线播放| 88av欧美| 极品人妻少妇av视频| 老司机亚洲免费影院| 在线免费观看的www视频| 在线观看免费午夜福利视频| 韩国av一区二区三区四区| 久久狼人影院| 国产精品偷伦视频观看了| 好看av亚洲va欧美ⅴa在| 在线永久观看黄色视频| 麻豆国产av国片精品| 超碰成人久久| 日韩欧美三级三区| 很黄的视频免费| 黄色成人免费大全| 国产亚洲欧美在线一区二区| 精品欧美一区二区三区在线| 999久久久精品免费观看国产| 天天添夜夜摸| 高清在线国产一区| 在线观看66精品国产| 国产精品一区二区三区四区久久 | 一本综合久久免费| 亚洲精品一二三| 精品国产一区二区久久| 脱女人内裤的视频| 999久久久国产精品视频| 午夜福利免费观看在线| 法律面前人人平等表现在哪些方面| 女同久久另类99精品国产91| 久久久水蜜桃国产精品网| 国产国语露脸激情在线看| 757午夜福利合集在线观看| 宅男免费午夜| 国产成人一区二区三区免费视频网站| 色老头精品视频在线观看| 十八禁网站免费在线| 狂野欧美激情性xxxx| 丰满饥渴人妻一区二区三| 精品国产乱码久久久久久男人| 级片在线观看| av天堂在线播放| 丝袜人妻中文字幕| 精品久久久精品久久久| 亚洲,欧美精品.| 欧美丝袜亚洲另类 | 日韩人妻精品一区2区三区| 新久久久久国产一级毛片| 亚洲欧美一区二区三区黑人| av天堂久久9| 女人高潮潮喷娇喘18禁视频| 十八禁人妻一区二区| 精品卡一卡二卡四卡免费| 精品久久久久久,| 少妇裸体淫交视频免费看高清 | 老司机亚洲免费影院| 涩涩av久久男人的天堂| 高清黄色对白视频在线免费看| 怎么达到女性高潮| 亚洲第一av免费看| 久久九九热精品免费| 色播在线永久视频| 亚洲精品一二三| 最新在线观看一区二区三区| 一级黄色大片毛片| 久久精品人人爽人人爽视色| 三级毛片av免费| 国产精品一区二区精品视频观看| 久久精品国产清高在天天线| 性少妇av在线| 久久香蕉精品热| 99国产综合亚洲精品| 午夜老司机福利片| 亚洲国产精品合色在线| 亚洲精品国产精品久久久不卡| 久久精品国产99精品国产亚洲性色 | 午夜福利一区二区在线看| 欧美日韩一级在线毛片| 欧美 亚洲 国产 日韩一| 久久精品国产亚洲av香蕉五月| 日本一区二区免费在线视频| www日本在线高清视频| 欧美日本中文国产一区发布| 人妻久久中文字幕网| 三上悠亚av全集在线观看| 欧美在线一区亚洲| 亚洲国产欧美网| 深夜精品福利| 老司机福利观看| 免费看a级黄色片| 国产欧美日韩一区二区精品| 精品福利观看| 黑丝袜美女国产一区| 正在播放国产对白刺激| 搡老岳熟女国产| 三级毛片av免费| 久久久久久免费高清国产稀缺| 亚洲一码二码三码区别大吗| 中文字幕高清在线视频| 久久中文看片网| 成年女人毛片免费观看观看9| 老熟妇仑乱视频hdxx| 十八禁人妻一区二区| 日韩精品中文字幕看吧| 91成年电影在线观看| 黄色a级毛片大全视频| 啦啦啦免费观看视频1| 人妻久久中文字幕网| 精品人妻在线不人妻| 欧美大码av| www.999成人在线观看| 精品国内亚洲2022精品成人| 国产一区二区三区在线臀色熟女 | 国产午夜精品久久久久久| 久久久久精品国产欧美久久久| 18禁观看日本| 麻豆av在线久日| 香蕉久久夜色| 国产免费av片在线观看野外av| 国产成人一区二区三区免费视频网站| 国产精品1区2区在线观看.| 国产亚洲欧美在线一区二区| 丰满人妻熟妇乱又伦精品不卡| 91字幕亚洲| 久久亚洲真实| 久久精品国产综合久久久| 久久精品国产亚洲av香蕉五月| 国产精品秋霞免费鲁丝片| 久热爱精品视频在线9| 窝窝影院91人妻| 久久精品aⅴ一区二区三区四区| 日本wwww免费看| 麻豆成人av在线观看| 国产一区二区三区在线臀色熟女 | 亚洲少妇的诱惑av| 久9热在线精品视频| 女同久久另类99精品国产91| 亚洲av熟女| 色综合婷婷激情| 午夜福利一区二区在线看| 黄色片一级片一级黄色片| 色哟哟哟哟哟哟| 两性夫妻黄色片| 欧美人与性动交α欧美软件| 在线天堂中文资源库| 欧美日韩av久久| 日韩免费av在线播放| 男人的好看免费观看在线视频 | www.自偷自拍.com| 最近最新中文字幕大全免费视频| 91精品国产国语对白视频| 日本精品一区二区三区蜜桃| 高清在线国产一区| 一a级毛片在线观看| 欧美中文日本在线观看视频| www.自偷自拍.com| 色综合欧美亚洲国产小说| 91麻豆精品激情在线观看国产 | 久久中文字幕一级| 亚洲专区国产一区二区| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日本亚洲视频在线播放| 免费av中文字幕在线| 国产欧美日韩精品亚洲av| 亚洲av片天天在线观看| 人成视频在线观看免费观看| 国产一区二区三区视频了| 老司机午夜十八禁免费视频| 国产又色又爽无遮挡免费看| 淫妇啪啪啪对白视频| 久久久精品欧美日韩精品| 欧美激情久久久久久爽电影 | 99热只有精品国产| 亚洲av成人一区二区三| 一级片免费观看大全| 激情视频va一区二区三区| 久久久精品国产亚洲av高清涩受| 精品一品国产午夜福利视频| 免费在线观看日本一区| 少妇裸体淫交视频免费看高清 | 国内毛片毛片毛片毛片毛片| 麻豆成人av在线观看| 免费在线观看黄色视频的| 精品久久久久久,| 国产高清videossex| 久久香蕉国产精品| 757午夜福利合集在线观看| 黄片大片在线免费观看| 人人妻,人人澡人人爽秒播| 亚洲国产中文字幕在线视频| 啦啦啦 在线观看视频| 中文字幕精品免费在线观看视频| 超色免费av| 神马国产精品三级电影在线观看 | 男女床上黄色一级片免费看| 美国免费a级毛片| 麻豆av在线久日| 亚洲全国av大片| 精品午夜福利视频在线观看一区| 男女床上黄色一级片免费看| 一本大道久久a久久精品| 国产成人影院久久av| 欧美日韩亚洲高清精品| 老司机福利观看| www.999成人在线观看| 无限看片的www在线观看| 国产成人影院久久av| 久久久久亚洲av毛片大全| 国产成人啪精品午夜网站| 99在线人妻在线中文字幕| 一本大道久久a久久精品| 国产成人影院久久av| 国产精品影院久久| 国产av一区在线观看免费| 亚洲av日韩精品久久久久久密| 欧美激情久久久久久爽电影 | av国产精品久久久久影院| 日本免费a在线| 午夜两性在线视频| 国产亚洲欧美98| 久久中文看片网| 香蕉久久夜色| 国产色视频综合| 在线视频色国产色| 国产在线精品亚洲第一网站| 日本免费a在线| 精品一区二区三区视频在线观看免费 | 91老司机精品| 国产一区二区在线av高清观看| 久久久久久久久中文| 99国产综合亚洲精品| 午夜福利影视在线免费观看| 日韩视频一区二区在线观看| 午夜91福利影院| 日韩欧美免费精品| 成人精品一区二区免费| 久久久久久久午夜电影 | 国产成年人精品一区二区 | 亚洲国产看品久久| 可以免费在线观看a视频的电影网站| 国产av在哪里看| 老司机午夜十八禁免费视频| 18禁观看日本| e午夜精品久久久久久久| av天堂在线播放| 亚洲欧美激情综合另类| 最好的美女福利视频网| 久久久久国产一级毛片高清牌| 搡老岳熟女国产| 免费av中文字幕在线| 欧美激情 高清一区二区三区| 免费观看精品视频网站| 日韩三级视频一区二区三区| 国产激情欧美一区二区| 国产免费现黄频在线看| 亚洲av成人一区二区三| 国产有黄有色有爽视频| 咕卡用的链子| 亚洲精品中文字幕一二三四区| 成年版毛片免费区| 日韩大码丰满熟妇| 12—13女人毛片做爰片一| 黄色片一级片一级黄色片| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲第一欧美日韩一区二区三区| 国产蜜桃级精品一区二区三区| 最近最新中文字幕大全电影3 | 在线观看免费日韩欧美大片| 亚洲国产欧美日韩在线播放| 在线观看日韩欧美| 亚洲一区二区三区色噜噜 | 老司机午夜福利在线观看视频| 男人舔女人下体高潮全视频| 免费久久久久久久精品成人欧美视频| 亚洲片人在线观看| 成人18禁在线播放| 国产一卡二卡三卡精品| 黄色a级毛片大全视频| 欧美激情极品国产一区二区三区| 国产精品野战在线观看 | 美女福利国产在线| 精品国产乱码久久久久久男人| 亚洲av五月六月丁香网| 一区二区三区精品91| 婷婷精品国产亚洲av在线| √禁漫天堂资源中文www| 1024香蕉在线观看| 男男h啪啪无遮挡| 亚洲av第一区精品v没综合| 国产视频一区二区在线看| 日本一区二区免费在线视频| 亚洲国产精品999在线| 夫妻午夜视频| 欧美黄色淫秽网站| 一a级毛片在线观看| 久久人妻av系列| 琪琪午夜伦伦电影理论片6080| 亚洲一区二区三区不卡视频| 久久中文字幕人妻熟女| 成年版毛片免费区| 国产精华一区二区三区| 一进一出抽搐gif免费好疼 | 精品一区二区三区视频在线观看免费 | 男人操女人黄网站| 欧美一区二区精品小视频在线| 美女 人体艺术 gogo|