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

    跨聲速風(fēng)洞槽壁干擾評估與修正技術(shù)的應(yīng)用

    2018-03-15 09:51:28劉光遠(yuǎn)魏志彭鑫陳德華1賈智亮
    航空學(xué)報(bào) 2018年2期
    關(guān)鍵詞:洞壁馬赫數(shù)試驗(yàn)段

    劉光遠(yuǎn),魏志,彭鑫,陳德華1,,賈智亮

    1.中國空氣動(dòng)力研究與發(fā)展中心 空氣動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室, 綿陽 621000 2.中國空氣動(dòng)力研究與發(fā)展中心 高速空氣動(dòng)力研究所,綿陽 621000

    跨聲速風(fēng)洞試驗(yàn)段通常采用開孔或開槽壁,高亞聲速范圍內(nèi),槽壁相比孔壁以干擾小、噪聲低、加工調(diào)整方便等特性逐漸成為跨聲速風(fēng)洞的發(fā)展趨勢[1],晚些時(shí)候建成的主要設(shè)備,如美國蘭利中心的國家跨聲速設(shè)備(National Transonic Facility,NTF)、德、英、荷、法四國合建的歐洲跨聲速風(fēng)洞(European Transonic Wind tunnel,ETW)等均采用了開槽形式。但開槽導(dǎo)致的透氣流動(dòng)使洞壁邊界處流動(dòng)異常復(fù)雜,增大了修正的難度,國外自20世紀(jì)40年代開始采用槽壁至今,開展了大量的干擾修正研究工作。

    早期研究主要以開槽參數(shù)調(diào)整的試驗(yàn)試湊法[2]為主,20世紀(jì)80年代后,逐漸出現(xiàn)了基于穿槽流動(dòng)特征和槽壁邊界干擾分析的經(jīng)典線性法、有限基本解法、壁壓信息法等工程修正技術(shù)[3]。2000年后,隨著數(shù)值模擬技術(shù)的發(fā)展,開始出現(xiàn)采用高精度CFD手段的研究成果[4-6],但其重點(diǎn)多為激波較強(qiáng)的跨聲速范圍[7],而且計(jì)算耗時(shí)長,無法滿足工程上實(shí)時(shí)修正的需求。目前大展弦比布局的軍/民用運(yùn)輸機(jī)試驗(yàn)范圍在馬赫數(shù)為0.60~0.90的高亞聲速段,該范圍內(nèi)經(jīng)典線性法仍然具有足夠的準(zhǔn)確度和可靠性[7],而且計(jì)算效率高、成本低,具有試驗(yàn)前評估和試驗(yàn)中實(shí)時(shí)修正的能力,是生產(chǎn)型風(fēng)洞中應(yīng)用最廣泛的修正方法,并在實(shí)踐中不斷改進(jìn)提高。

    歐美主要試驗(yàn)機(jī)構(gòu)均具備基于經(jīng)典方法的槽壁干擾快速評估和修正能力。NASA Ames研究中心在2000年開發(fā)了ANTARES程序[8],目前已應(yīng)用于11 ft(1 ft=0.304 8 m)跨聲速風(fēng)洞(Transonic Wind Tunnel, TWT)和12 ft增壓風(fēng)洞(Pressure Wind Tunnel, PWT)中[9-10],但該程序計(jì)算洞壁擾動(dòng)速度場時(shí),需要?jiǎng)澐侄幢诰W(wǎng)格,當(dāng)試驗(yàn)?zāi)P臀恢冒l(fā)生變化時(shí),網(wǎng)格需要進(jìn)行相應(yīng)的調(diào)整。Langley研究中心NTF采用基于高階面元法的PANCOR程序[11]評估和修正槽壁干擾,在模型表面和尾流區(qū)域布置面元,能夠提高模型擾流模擬的準(zhǔn)確性,因此該程序除修正洞壁干擾外,還能夠計(jì)算飛行器的氣動(dòng)特性,但前期外形建模和面元?jiǎng)澐止ぷ鬏^繁瑣。美國波音跨聲速風(fēng)洞(Boeing Transonic Wind Tunnel, BTWT)[12]開發(fā)了基于空間笛卡兒網(wǎng)格的TRANAIR程序,并利用翼身組合體標(biāo)模開展了模型翼展對洞壁干擾的影響分析工作[13],該程序利用理想均勻邊界條件求解全速勢方程,并添加了附面層的黏性影響修正,能夠用于馬赫數(shù)接近1.0的跨聲速范圍內(nèi)的模型氣動(dòng)特性評估和洞壁干擾修正工作,但與高階面元法類似,評估計(jì)算的前期準(zhǔn)備工作繁瑣。

    國內(nèi)方面,近年來中國空氣動(dòng)力研究與發(fā)展中心(CARDC)采用CFD方法初步研究了槽壁試驗(yàn)段的流動(dòng)特征[14-15],西北工業(yè)大學(xué)開展了開槽幾何參數(shù)對跨聲速流場品質(zhì)的研究工作[16],但內(nèi)容暫未涉及槽壁邊界對模型繞流的干擾修正。高速透氣壁干擾修正工作均針對孔壁,主要方法為試驗(yàn)和計(jì)算相結(jié)合的壁壓信息法[3,17-20],不具備試驗(yàn)前評估的能力。

    2013年CARDC完成了國內(nèi)第一座2 m量級槽壁試驗(yàn)段的流場調(diào)試工作[21],正在建造的先進(jìn)跨聲速設(shè)備也將采用槽壁,但目前國內(nèi)對槽壁干擾特性認(rèn)識尚不成熟,缺乏工程修正經(jīng)驗(yàn)和方法,而槽壁干擾修正技術(shù)的不足又制約了槽壁的大規(guī)模工程應(yīng)用。因此,研究槽壁干擾特性,評估并修正其干擾就顯得尤為迫切。本文基于洞壁處的小擾動(dòng)均勻邊界條件和經(jīng)典方法對槽壁干擾特性進(jìn)行評估和修正,并利用民機(jī)標(biāo)模進(jìn)行準(zhǔn)確性驗(yàn)證,以期為將來建立工程化的槽壁干擾修正體系提供借鑒。

    1 經(jīng)典線性法

    1.1 控制方程

    對于亞聲速可壓縮流動(dòng),擾動(dòng)速度勢滿足以下線性方程:

    (1)

    式中:Ma∞為來流馬赫數(shù);φ為擾動(dòng)速度勢,可以分解為模型在自由流場中的速度勢φm和洞壁誘導(dǎo)的擾動(dòng)速度勢φi。φm為已知量,可以利用φm+φi在洞壁邊界處滿足特定的透氣壁邊界條件來求解φi。洞壁在模型區(qū)域誘導(dǎo)的擾動(dòng)速度為

    (2)

    (3)

    流向擾動(dòng)速度ui給出了模型實(shí)體及其尾流的堵塞干擾影響,法向擾動(dòng)速度vi給出了升力干擾影響,二者沿流向的變化(梯度)分別給出了浮阻影響和流線彎曲影響。

    1.2 邊界條件

    透氣壁邊界條件的描述方法是影響修正準(zhǔn)度的關(guān)鍵因素。對于孔壁邊界,可采用縮比模型流量測量法[7,22]獲得邊界條件;但對于槽壁邊界,目前尚缺少有效的測試方法[23],工程應(yīng)用中以理想槽壁均勻邊界條件對穿槽流動(dòng)進(jìn)行簡化,即

    (4)

    式中:Cps和θs分別為槽壁邊界處壓力系數(shù)和氣流偏角。理想槽壁均勻邊界條件認(rèn)為二者之間存在線性關(guān)系,比例系數(shù)K稱為槽壁流動(dòng)彎曲系數(shù),其量值僅為槽外形參數(shù)的函數(shù)[24],即

    (5)

    式中:a為槽間距;d為槽寬;如圖1所示。

    小擾動(dòng)條件下,式(4)中的壓力系數(shù)Cps和氣流偏角θs可分別寫為流向和法向擾動(dòng)速度φx、φy的函數(shù)形式,見圖1右下角所示,U∞為來流速度,因此式(4)可以寫為

    φx+Kφxy=0

    (6)

    圖1 槽壁理想均勻邊界條件Fig.1 Ideal homogeneous boundary condition for slotted walls

    目前生產(chǎn)型風(fēng)洞槽壁干擾修正均采用式(5)和式(6)定義的理想均勻邊界條件。使用該方法的前提條件是等寬度槽沿?zé)o限長的展向均勻分布,近壁區(qū)域流動(dòng)為線性位勢流,以及流動(dòng)的小擾動(dòng)條件。因此要求模型尺度相對于試驗(yàn)段為小量,忽略洞壁附近流動(dòng)的黏性影響,超臨界范圍內(nèi)模型表面波系未發(fā)展至洞壁區(qū)域,即試驗(yàn)段流場中的非線性流動(dòng)區(qū)域僅存在于近模型區(qū)域,見圖2。

    式(5)和式(6)中的系數(shù)K具有長度單位,由于常規(guī)縱向試驗(yàn)中側(cè)壁干擾相對于上下壁干擾為小量,因此利用試驗(yàn)段高度H對K進(jìn)行歸一化,引入無量綱系數(shù)P:

    (7)

    P=0即K/H→∞時(shí)為實(shí)壁,P=1即K/H=0時(shí)為開口壁。

    圖2 槽壁試驗(yàn)段位勢流動(dòng)Fig.2 Potential flow in test section with slotted walls

    1.3 干擾影響因子

    通過Prandtl-Glauert轉(zhuǎn)換將偏微分方程式(1)變?yōu)長aplace方程,結(jié)合槽壁邊界條件式(6),采用有限基本解方法[25]求解洞壁誘導(dǎo)速度。為能夠在試驗(yàn)中實(shí)時(shí)修正,不劃分模型表面網(wǎng)格,而利用已知解的偶極子、點(diǎn)源分別模擬模型的實(shí)體堵塞和尾流干擾效應(yīng),利用馬蹄渦模擬機(jī)翼的升力效應(yīng),各奇點(diǎn)的強(qiáng)度由模型體積以及升、阻力系數(shù)確定。也就是說僅考慮模型繞流的積分結(jié)果,控制點(diǎn)坐標(biāo)r、θ,以及試驗(yàn)段高度H、寬度B的意義見圖3。用于模擬模型擾流的奇點(diǎn)位于試驗(yàn)段中心的原點(diǎn)位置,擾動(dòng)速度勢φm及其對應(yīng)的強(qiáng)度計(jì)算方法見表1,其中V、S、s分別為模型體積、機(jī)翼參考面積和翼展長度,μm、m、Г分別為偶極子、點(diǎn)源、馬蹄渦的強(qiáng)度,β為Prandtl-Glauert壓縮性因子,是馬赫數(shù)的函數(shù)。

    而不關(guān)注其繞流細(xì)節(jié)??紤]到本方法僅用于評估槽壁干擾量而不計(jì)算模型氣動(dòng)特性,φm僅應(yīng)用于洞壁邊界,因此模型區(qū)域φm的誤差不會(huì)顯著影響φi的計(jì)算結(jié)果。

    根據(jù)以上方法計(jì)算得到的實(shí)壁(P=0)和槽壁(P=0.25)誘導(dǎo)的流向和法向擾動(dòng)速度沿x軸的分布見圖4(ui/U∞,vi/U∞),計(jì)算模型的實(shí)體堵塞度為1.0%,升力系數(shù)為-0.20~0.20。可以看出,與實(shí)壁狀態(tài)相比,槽壁能夠有效降低壁板誘導(dǎo)的擾動(dòng)速度:模型中心處的流向擾動(dòng)速度(與來流速度的比值)從實(shí)壁狀態(tài)的0.015降低至0.004,法向擾動(dòng)速度從0.002~0.005降低至0.001~0.003;同時(shí),由于流向擾動(dòng)速度分布關(guān)于模型中心對稱,因此,槽壁中的實(shí)體堵塞不會(huì)產(chǎn)生由于流向擾動(dòng)速度梯度導(dǎo)致的阻力變化(浮阻)。

    圖3 計(jì)算洞壁擾動(dòng)速度使用的坐標(biāo)系Fig.3 Coordinate system for calculating interference velocities induced by tunnel walls

    表1 試驗(yàn)?zāi)P蛿_動(dòng)速度勢及其強(qiáng)度Table 1 Potentials and strengths of interference velocity of testing model

    干擾項(xiàng)奇點(diǎn)類型擾動(dòng)速度勢φm強(qiáng)度說明實(shí)體堵塞偶極子μm4π·xx2+β2r2()3/2μm=U∞Vμm與模型體積V相關(guān)尾流干擾點(diǎn)源-m4π1x2+β2r2()3/2m=0.5U∞CDSm與模型阻力系數(shù)CD相關(guān)升力干擾馬蹄渦Γs2π1+xx2+β2r2()1/2sinθrΓs=0.5U∞CLSΓs與模型升力系數(shù)CL相關(guān),升力L=ρU∞Γs=0.5ρU2∞CLS

    圖5給出了利用本方法計(jì)算得到的高速客機(jī)CHNT-1標(biāo)模干擾因子隨無量綱槽壁參數(shù)P的變化曲線,升、阻特性由該標(biāo)模在2.4 m跨聲速風(fēng)洞槽壁試驗(yàn)段結(jié)果給出。為驗(yàn)證準(zhǔn)確性,圖5還給出了ETW提供的評估結(jié)果,狀態(tài)為Ma=0.78。圖6給出了升力干擾因子δ0隨馬赫數(shù)的變化關(guān)系??梢钥闯觯?/p>

    1) 本文方法的計(jì)算結(jié)果規(guī)律正常,干擾因子量值隨開槽系數(shù)P的增大向負(fù)值(開口趨勢)發(fā)展。兩座設(shè)備的結(jié)果對比表明,槽壁特性除與開槽參數(shù)有關(guān)外,還與試驗(yàn)段高度有關(guān),相同槽參數(shù)下,高度越大,槽壁干擾越趨于開口特性。2.4 m跨聲速試驗(yàn)段比ETW高0.40 m,計(jì)算結(jié)果表明堵塞因子比ETW量值略小。

    2) 計(jì)算結(jié)果與ETW的評估結(jié)果一致性較好,測試點(diǎn)處量值差異小于1.0%。

    3) 考慮流動(dòng)壓縮性時(shí),升力干擾因子隨馬赫數(shù)增大而向負(fù)值增大,與ETW參考結(jié)果規(guī)律一致。

    因此,本方法計(jì)算結(jié)果與國外參考數(shù)據(jù)一致,可以用于槽壁干擾的評估和分析。

    圖4 槽壁誘導(dǎo)擾動(dòng)速度分布Fig.4 Distributions of interference velocities induced by slotted walls

    圖5 槽壁系數(shù)對洞壁干擾因子的影響Fig.5 Influence of slot parameter on wall interference factors

    圖6 馬赫數(shù)對升力干擾因子的影響Fig.6 Influence of Mach number on lift interference factor

    1.4 修正流程

    獲得槽壁干擾影響因子后,逐項(xiàng)對風(fēng)軸系氣動(dòng)力(矩)系數(shù)進(jìn)行修正,程序流程見圖7,其中ε為堵塞因子,下標(biāo)w、f、wk分別代表機(jī)翼、機(jī)身和尾流堵塞效應(yīng);ΩS為開槽壁的堵塞干擾因子。

    最終結(jié)果包括馬赫數(shù)修正(堵塞干擾),迎角修正(升力干擾、流線彎曲干擾),升、阻力修正(堵塞干擾、升力干擾、尾流浮阻),以及俯仰力矩修正(堵塞干擾、流線彎曲干擾),即

    圖7 槽壁干擾評估/修正程序流程圖Fig.7 Flow chart of assessment/correction sequence for slotted wall interference

    (8)

    式中:下標(biāo)u、c分別表示修正前、后的數(shù)據(jù);CQ為馬赫數(shù)修正變化引起的速壓修正系數(shù),即

    (9)

    尾流浮阻ΔCD,wk采用考慮流動(dòng)分離影響的Maskell方法[26]計(jì)算,即

    (10)

    CDS為扣除升力產(chǎn)生的誘導(dǎo)阻力影響的模型廢阻力,由模型極曲線求解。當(dāng)用于試驗(yàn)前評估時(shí),以附著流模型代替Maskell方法,式(10)中的CDS替換為估算的最小阻力或零阻。

    Δα1和Δαsc分別為升力干擾、流線彎曲干擾誘導(dǎo)的迎角增量,即

    (11)

    (12)

    迎角變化引起的升、阻力系數(shù)修正量通過修正前后的風(fēng)軸系矢量轉(zhuǎn)換得到,即

    ΔCLc,α={CLu[cos(Δα)-1]-CDusin(Δα)}CQ

    ΔCDc,α={CDu[cos(Δα)-1]+CDucos(Δα)}CQ

    (13)

    流線彎曲誘導(dǎo)的俯仰力矩增量ΔCm,sc為

    (14)

    從以上修正方法可以看出,試驗(yàn)段尺寸除影響槽壁特性系數(shù)P以及相應(yīng)的干擾因子外,還影響堵塞度和迎角修正量。槽壁干擾量正比于干擾因子與試驗(yàn)段橫截面積的比值。

    2 模型與槽壁試驗(yàn)介紹

    CHNT-1標(biāo)模是中國最新的大飛機(jī)標(biāo)模,采用大展弦比,超臨界機(jī)翼,大收縮后機(jī)身,是先進(jìn)大飛機(jī)布局的典型代表,外形見圖8,模型尺寸參數(shù)見表2。

    圖8 CHNT-1標(biāo)模外形Fig.8 Standard model CHNT-1 shape

    ETW與2.4 m跨聲速風(fēng)洞均采用上下開槽、左右實(shí)壁的形式,開槽數(shù)均為6條,不同的是2.4 m跨聲速風(fēng)洞包括2條半槽,而ETW為6條全槽,兩座風(fēng)洞的槽寬、槽間距及分布有較大差異,開槽參數(shù)對比見表3,壁板外形見圖9。

    由于理想均勻邊界條件的前提條件之一是等寬度槽沿?zé)o限長展向(即z方向,見圖1和圖3)均勻分布,因此側(cè)邊半槽不會(huì)對槽壁的K系數(shù)產(chǎn)生影響,其量值僅與槽寬d和槽間距a有關(guān)。而在計(jì)算槽壁P系數(shù)時(shí),以試驗(yàn)段高度H進(jìn)行歸一化,因此雖然兩座風(fēng)洞的K系數(shù)差異較大,但P系數(shù)仍然較接近。

    表2 CHNT-1標(biāo)模幾何參數(shù)Table 2 Geometrical parameters of standard model CHNT-1

    表3 槽壁幾何參數(shù)Table 3 Geometrical parameters of slotted walls

    圖9 兩座風(fēng)洞的開槽形式Fig.9 Slotted walls in two wind tunnels

    3 槽壁干擾修正與分析

    3.1 數(shù)據(jù)修正前后對比

    利用本文方法對原始數(shù)據(jù)進(jìn)行槽壁干擾修正,修正前后的對比見圖10,圖中還給出了ETW的參考結(jié)果??梢钥闯?,槽壁對CHNT-1標(biāo)模的影響干擾主要體現(xiàn)為升力干擾,即迎角修正。修正后升力線斜率略增大,與ETW參考結(jié)果一致性較好。

    圖10 CHNT-1標(biāo)模洞壁干擾修正前后對比(Ma=0.78)Fig.10 Comparison between corrected and uncorrected data for standard model CHNT-1 tunnel wall interference (Ma=0.78)

    各馬赫數(shù)下迎角修正量Δαwi隨模型迎角的變化見圖11,在不同馬赫數(shù)范圍內(nèi),本文方法的修正結(jié)果與ETW評估結(jié)果一致。CHNT-1標(biāo)模迎角修正量與馬赫數(shù)呈正比,隨馬赫數(shù)的增大而增大,量值從Ma=0.20時(shí)的0.05°增大至Ma=0.90時(shí)的0.20°。結(jié)合圖6,主要原因是高馬赫數(shù)流動(dòng)壓縮性使升力干擾因子增大。本文計(jì)算結(jié)果與ETW評估結(jié)果在升力線性段內(nèi)差異小于0.01°,說明采用的計(jì)算方法是準(zhǔn)確、可靠的。

    圖11 CHNT-1標(biāo)模迎角修正量對比Fig.11 Comparison of angle of attack correction for standard model CHNT-1

    3.2 兩座風(fēng)洞修正量對比

    Ma=0.70,0.78,0.85時(shí),在設(shè)計(jì)升力系數(shù)CL=0.48狀態(tài)下計(jì)算了ETW和2.4 m跨聲速風(fēng)洞的槽壁干擾量,槽壁的升力干擾Δαwi、堵塞干擾ΔMawi、浮阻干擾ΔCD,wk及流線彎曲干擾ΔCm,SC隨馬赫數(shù)的變化,對比結(jié)果見圖12。

    1) 雖然兩座風(fēng)洞的P系數(shù)量值較接近,但由于高度不同,干擾量存在明顯差異:堵塞度的差異使ETW中的實(shí)體堵塞干擾更嚴(yán)重,導(dǎo)致馬赫數(shù)略低;ETW升力干擾因子比2.4 m跨聲速風(fēng)洞大0.025左右,加之試驗(yàn)段截面積小,因此由式(11)和式(12)計(jì)算的迎角修正量更大。

    圖12 2.4 m風(fēng)洞與ETW槽壁干擾量對比(CL=0.48)Fig.12 Comparison of slotted wall interference for 2.4 m wind tunnel and ETW (CL=0.48)

    2) 對于CHNT-1標(biāo)模,槽壁干擾主要表現(xiàn)為升力干擾,干擾量值隨馬赫數(shù)的增大而增大。2.4 m跨聲速風(fēng)洞中,設(shè)計(jì)點(diǎn)處槽壁干擾使模型迎角降低0.06°~0.08°,ETW中使模型迎角降低0.08°~0.11°。

    3) 2.4 m跨聲速風(fēng)洞中,堵塞干擾使馬赫數(shù)略增大,尾流浮阻使模型阻力略增大,ETW趨勢相反。對比圖5(a)可以發(fā)現(xiàn),雖然兩風(fēng)洞的堵塞干擾因子量值均較小,但符號相反,2.4 m跨聲速風(fēng)洞采用大開閉比的寬槽,干擾影響趨于開口邊界,而ETW采用中等開閉比的窄槽,干擾影響趨于實(shí)壁邊界。

    4) 對于CHNT-1標(biāo)模,堵塞干擾、尾流浮阻以及流線彎曲干擾等量值均較小,試驗(yàn)數(shù)據(jù)修正中,該量級可忽略不計(jì)。這是由于該標(biāo)模堵塞度較小,在2 m量級口徑風(fēng)洞中堵塞度小于或接近1.0%的要求,因此堵塞影響可不修正;同時(shí),CL=0.48時(shí)機(jī)翼繞流均為附著流,尾流范圍和強(qiáng)度小,故尾流浮阻影響可不修正。

    5) 考慮到CHNT-1標(biāo)模類構(gòu)型模型升力線斜率為0.10~0.15,由于迎角修正導(dǎo)致的升力系數(shù)差異將接近0.006~0.015,遠(yuǎn)超過國軍標(biāo)的精度指標(biāo),對模型失速性能預(yù)測和升阻比評估等都將產(chǎn)生明顯影響。因此,該類構(gòu)型模型即使符合堵塞度小于1.0%的試驗(yàn)要求,數(shù)據(jù)也必須進(jìn)行干擾修正,至少進(jìn)行升力干擾修正后才能使用。

    4 結(jié) 論

    基于經(jīng)典方法進(jìn)行了槽壁干擾的評估與修正工作,利用CHNT-1標(biāo)模進(jìn)行了修正準(zhǔn)度驗(yàn)證,在高亞聲速范圍內(nèi),對比分析了2.4 m跨聲速風(fēng)洞和ETW槽壁對大展弦比構(gòu)型模型試驗(yàn)數(shù)據(jù)的干擾規(guī)律和差異。

    1) 采用的計(jì)算方法準(zhǔn)確可行,兼具預(yù)先評估的能力,既可用于國內(nèi)高速風(fēng)洞的槽壁干擾評估、修正,也可用于槽壁外形參數(shù)的設(shè)計(jì)與影響分析。

    2) 開槽參數(shù)與試驗(yàn)段尺寸對槽壁干擾效應(yīng)和修正量的影響較大,合理設(shè)計(jì)試驗(yàn)段和壁板外形是提高數(shù)據(jù)準(zhǔn)度的重要手段。

    3) 大展弦比高升力構(gòu)型模型常規(guī)試驗(yàn)中,堵塞、尾流和流線彎曲干擾可忽略,但升力干擾引起的迎角增量必須修正。

    [1] HIEBERT D, STEINLE F. A historical perspective of design requirements for AEDC’s propulsion wind tunnel and von Kármán Facilities: AIAA-2010-0138[R]. Reston, VA: AIAA, 2010.

    [2] LITTLE J B, CUBBAGE J J. The development of an 8-inch by 8-inch slotted tunnel for Mach numbers up to 1.28: NACA-RM-L55B08[R]. Washington, D.C.: NACA, 1955.

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

    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 (in Chinese).

    [4] NORBERT G, JURGEN Q. The ETW wall interference assessment for full and half models: AIAA-2004-0769 [R]. Reston, VA: AIAA, 2004.

    [5] MASELAND J E J, LABAN M, VEN H V D, et al. Development of CFD-based interference models for the DNW HST transonic wind tunnel: AIAA-2006-3639[R]. Reston, VA: AIAA, 2006.

    [6] KOHZAI M, UENO M, KOGA S, et al. Wall and support interference corrections of NASA common research model wind tunnel tests in JAXA: AIAA-2013-0963[R]. Reston, VA: AIAA, 2013.

    [7] GLAZKOV S A, GORBUSHIN A R, IVANOV A I, et al. Recent experience in improving the accuracy of wall interference corrections in TsAGI T-128 wind tunnel[J]. Progress in Aerospace Sciences, 2001, 37(3): 263-298.

    [8] ULBRICH N. The application of panel method code ANTARES to wind tunnel wall interference problems[C]∥AIAA Aerospace Sciences Meeting & Exhibit. Reston, VA: AIAA, 2013.

    [9] ULBRICH N, BOONE A. Direct validation of the wall interference correction system of the Ames 11-foot transonic wind tunnel: NASA/TM-20030212268[R]. Washington, D.C.: NASA, 2003.

    [10] ULBRICH N. The real-time wall interference correction system of the NASA Ames 12-foot pressure tunnel: NASA/CR-1998-208537[R]. Washington, D.C.: NASA,1998.

    [11] WILLIAM B. User’s guide to PANCOR: A panel method program for interference assessment in slotted wall wind tunnels: NASA CR 187497 [R].Washington, D.C.: NASA, 1990.

    [12] KRYNYTZKY A, HERGERT D. Boeing transonic wind tunnel upgrade assessment: AIAA-2001-0452[R]. Reston, VA: AIAA, 2001.

    [13] KRYNYTZKY A J. Parametric model size study of wall interference in the boeing transonic wind tunnel using TRANAIR: AIAA-2004-2310[R]. Reston, VA: AIAA, 2004.

    [14] 叢成華, 彭強(qiáng), 王海鋒. 槽壁試驗(yàn)段低超聲速流場特性數(shù)值模擬[J]. 航空學(xué)報(bào), 2010, 31(12): 2302-2308.

    CONG C H, PENG Q, WANG H F. Numerical simulation on characteristics of slotted test section of low supersonic field[J]. Acta Aeronautica et Astronautica Sinica, 2010, 31(12): 2302-2308 (in Chinese).

    [15] 汪伏波. 跨聲速風(fēng)洞槽壁試驗(yàn)段流動(dòng)特性初步研究[D]. 綿陽: 中國空氣動(dòng)力研究與發(fā)展中心, 2010.

    WANG F B.Preliminary research on flow characteristic of transonic wind tunnel slotted test section[D]. Mianyang: China Aerodynamics Research and Development Center, 2010 (in Chinese).

    [16] 鞠煉, 白俊強(qiáng), 郭斌, 等. 槽壁幾何參數(shù)對跨聲速風(fēng)洞流場品質(zhì)的影響[J]. 航空學(xué)報(bào), 2016, 37(5): 1440-1453.

    JU L, BAI J Q, GUO B, et al. Effect of geometry parameters on flow field quality in a transonic wind tunnel[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(5):1440-1453 (in Chinese).

    [17] 尹陸平, 賀中, 于志松, 等. 亞聲速大迎角模型試驗(yàn)洞壁干擾修正方法研究[J]. 流體力學(xué)實(shí)驗(yàn)與測量, 2000, 14(3): 37-41.

    YIN L P, HE Z, YU Z S, et al. 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).

    [18] 張其威. 壁壓法用于高速風(fēng)洞洞壁干擾修正[J]. 南京航空學(xué)院學(xué)報(bào), 1992(5): 496-505.

    ZHANG Q W. Wall-signature methods for high speed wind tunnel wall interference corrections[J]. Journal of Nanjing Aeronautical Institute, 1992(5): 496-505 (in Chinese).

    [19] 范召林, 張玉倫, 賀中. 跨音速飛機(jī)模型試驗(yàn)洞壁干擾數(shù)值模擬的初步研究[J]. 航空學(xué)報(bào), 1997, 18(2): 210-214.

    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).

    [20] 范召林, 惲起麟, 賀中. 跨聲速洞壁干擾研究述評[J]. 實(shí)驗(yàn)流體力學(xué), 1996, 10(4): 1-11.

    FAN Z L, YUN Q L, HE Z. Areview on transonic wall interference researches[J]. Journal of Experiments in Fluid Mechanics, 1996, 10(4): 1-11 (in Chinese).

    [21] 陳德華, 王瑞波, 劉光遠(yuǎn), 等. 2.4 m跨聲速風(fēng)洞槽壁試驗(yàn)段調(diào)試及流場校測[J]. 實(shí)驗(yàn)流體力學(xué), 2013, 27(4): 66-70.

    CHEN D H, WANG R B, LIU G Y, et al. The test and flow field calibration of 2.4 m transonic wind tunnel slotted test section[J]. Journal of Experiments in Fluid Mechanics, 2013, 27(4): 66-70 (in Chinese).

    [22] MACMANUS D G, EATON J A. Measurement and analysis of the flowfields induced by suction perforations[J]. AIAA Journal, 1998, 36(9): 1553-1561.

    [23] GLAZKOV S, GORBUSHIN A, IVANOV A, et al. Numerical and experimental investigations of slot flow with respect to wind tunnel wall interference assessment: AIAA-2004-2308[R]. Reston, VA: AIAA, 2004.

    [24] EWALD B F R. Wind tunnel wall correction[R]. Paris: AGARD,1998.

    [25] 陳作斌. 計(jì)算流體力學(xué)及應(yīng)用[M]. 北京: 國防工業(yè)出版社, 2002.

    CHEN Z B. Computational fluid mechanics and applications[M]. Beijing: National Defense Industry Press, 2002 (in Chinese).

    [26] HACKETT J E. Tunnel-induced gradients and their effect on drag[J]. AIAA Journal, 1996, 34(12):2575-2581.

    猜你喜歡
    洞壁馬赫數(shù)試驗(yàn)段
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    “洞壁”聽“琴音”(3)
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    跨聲速風(fēng)洞槽壁試驗(yàn)段流場品質(zhì)提升措施研究
    載荷分布對可控?cái)U(kuò)散葉型性能的影響
    側(cè)壓影響下圓形洞室?guī)r爆雙軸物理模擬試驗(yàn)研究
    合肥軌道交通1號線試驗(yàn)段Ⅰ標(biāo)基坑監(jiān)測及結(jié)果分析研究
    HL-2M真空室試驗(yàn)段制造工藝技術(shù)
    中國核電(2017年2期)2017-08-11 08:00:54
    小展弦比飛翼標(biāo)模FL-2風(fēng)洞跨聲速開孔壁干擾特性修正研究
    三種跨聲速洞壁干擾修正方法及其在小展弦比飛翼標(biāo)模試驗(yàn)中的應(yīng)用
    亚洲av免费高清在线观看| 十八禁高潮呻吟视频 | 观看av在线不卡| 国产黄频视频在线观看| 国产一区亚洲一区在线观看| 丰满乱子伦码专区| 青春草亚洲视频在线观看| 国产在视频线精品| 国产精品麻豆人妻色哟哟久久| 高清视频免费观看一区二区| 日本色播在线视频| 9色porny在线观看| 热re99久久精品国产66热6| 黄色欧美视频在线观看| 免费观看av网站的网址| 久久久久久伊人网av| 在线观看美女被高潮喷水网站| 国产亚洲精品久久久com| 欧美成人精品欧美一级黄| 国产在线男女| 只有这里有精品99| 一级毛片久久久久久久久女| 嘟嘟电影网在线观看| 国产伦在线观看视频一区| 午夜日本视频在线| 视频中文字幕在线观看| 国产无遮挡羞羞视频在线观看| 欧美高清成人免费视频www| 22中文网久久字幕| 亚洲怡红院男人天堂| 久久久精品94久久精品| 最近中文字幕高清免费大全6| 国产成人a∨麻豆精品| 国产免费一区二区三区四区乱码| 少妇猛男粗大的猛烈进出视频| videossex国产| 久久国内精品自在自线图片| 熟女人妻精品中文字幕| 国产精品99久久99久久久不卡 | a级毛片免费高清观看在线播放| 亚洲欧洲日产国产| 色视频在线一区二区三区| 99久久综合免费| 一本一本综合久久| 伊人久久国产一区二区| 成年美女黄网站色视频大全免费 | 免费观看的影片在线观看| 高清在线视频一区二区三区| 亚洲电影在线观看av| 桃花免费在线播放| 美女cb高潮喷水在线观看| 乱人伦中国视频| 人人妻人人澡人人爽人人夜夜| 亚洲天堂av无毛| 国产综合精华液| 国产视频内射| 亚洲国产精品999| 免费人妻精品一区二区三区视频| 黄色毛片三级朝国网站 | 五月玫瑰六月丁香| 精品人妻熟女av久视频| 国产av一区二区精品久久| 大话2 男鬼变身卡| 免费播放大片免费观看视频在线观看| 国模一区二区三区四区视频| 国产日韩欧美在线精品| 免费观看av网站的网址| 丰满人妻一区二区三区视频av| 国精品久久久久久国模美| 香蕉精品网在线| 在线精品无人区一区二区三| 国内揄拍国产精品人妻在线| 99九九在线精品视频 | 国产亚洲午夜精品一区二区久久| 麻豆精品久久久久久蜜桃| 少妇被粗大的猛进出69影院 | 午夜免费男女啪啪视频观看| 欧美精品一区二区免费开放| 久久久久久久精品精品| av线在线观看网站| 国产高清不卡午夜福利| 日本猛色少妇xxxxx猛交久久| 简卡轻食公司| 国产在线视频一区二区| 亚洲一级一片aⅴ在线观看| 免费观看的影片在线观看| 日韩视频在线欧美| 欧美精品一区二区免费开放| 亚洲高清免费不卡视频| 日日啪夜夜爽| 成人免费观看视频高清| 大码成人一级视频| 国产片特级美女逼逼视频| 亚洲三级黄色毛片| 国产精品偷伦视频观看了| 国产精品女同一区二区软件| 亚洲精品乱码久久久v下载方式| 成人亚洲精品一区在线观看| 各种免费的搞黄视频| 最近手机中文字幕大全| 久久久欧美国产精品| 汤姆久久久久久久影院中文字幕| 少妇裸体淫交视频免费看高清| 亚洲国产欧美日韩在线播放 | 国产老妇伦熟女老妇高清| 国产一区有黄有色的免费视频| 精品熟女少妇av免费看| 狂野欧美激情性xxxx在线观看| 热99国产精品久久久久久7| 又黄又爽又刺激的免费视频.| 高清黄色对白视频在线免费看 | 女人久久www免费人成看片| 99久久精品热视频| 亚洲欧美一区二区三区国产| 熟女av电影| 亚洲精品日韩av片在线观看| 一本久久精品| 久久久欧美国产精品| 纵有疾风起免费观看全集完整版| 综合色丁香网| 久久久久久人妻| 黑丝袜美女国产一区| 亚洲精品乱久久久久久| 国产精品久久久久久精品电影小说| 在线观看一区二区三区激情| 国产亚洲5aaaaa淫片| 日韩成人av中文字幕在线观看| 亚洲人成网站在线播| 最近中文字幕2019免费版| xxx大片免费视频| 熟女人妻精品中文字幕| 亚洲经典国产精华液单| 日韩av免费高清视频| 少妇熟女欧美另类| 欧美日本中文国产一区发布| 亚洲av男天堂| 久久精品久久久久久久性| 少妇高潮的动态图| 亚洲欧洲日产国产| 内地一区二区视频在线| 亚洲欧美一区二区三区黑人 | 国产在线视频一区二区| 国产亚洲欧美精品永久| 青春草国产在线视频| 成年人免费黄色播放视频 | 极品教师在线视频| 精品熟女少妇av免费看| 啦啦啦中文免费视频观看日本| 免费黄网站久久成人精品| 肉色欧美久久久久久久蜜桃| 久久午夜综合久久蜜桃| 国产高清三级在线| 久久国产乱子免费精品| 大香蕉97超碰在线| kizo精华| 最近最新中文字幕免费大全7| 熟女人妻精品中文字幕| 国产精品一区二区在线观看99| 大码成人一级视频| 国产乱来视频区| 少妇高潮的动态图| a级毛片免费高清观看在线播放| 亚洲国产欧美日韩在线播放 | 日韩一区二区三区影片| 久久久久久久久久久免费av| 自拍欧美九色日韩亚洲蝌蚪91 | 99久久精品热视频| 99热全是精品| 人体艺术视频欧美日本| 美女xxoo啪啪120秒动态图| 99热这里只有精品一区| 高清黄色对白视频在线免费看 | 80岁老熟妇乱子伦牲交| 国内少妇人妻偷人精品xxx网站| av免费在线看不卡| 嫩草影院入口| 超碰97精品在线观看| 久久这里有精品视频免费| 人妻制服诱惑在线中文字幕| 91久久精品电影网| 最近最新中文字幕免费大全7| 黄色怎么调成土黄色| 午夜激情福利司机影院| 一个人免费看片子| 久久国产精品大桥未久av | 麻豆精品久久久久久蜜桃| 精品一区二区三卡| 亚洲av二区三区四区| 国产精品无大码| 2018国产大陆天天弄谢| 王馨瑶露胸无遮挡在线观看| 免费黄网站久久成人精品| 我的老师免费观看完整版| 国产综合精华液| 97在线人人人人妻| av免费观看日本| 国产亚洲精品久久久com| 久久av网站| 一级毛片 在线播放| 99热6这里只有精品| av天堂中文字幕网| 午夜免费鲁丝| 久久这里有精品视频免费| 久久精品国产自在天天线| 在线 av 中文字幕| 国产黄片美女视频| 视频中文字幕在线观看| 国产av一区二区精品久久| 久久精品国产自在天天线| 一本一本综合久久| 老司机影院毛片| 久久久精品94久久精品| 国产乱人偷精品视频| 日韩人妻高清精品专区| 亚洲综合精品二区| 欧美日本中文国产一区发布| √禁漫天堂资源中文www| 久久久久久人妻| 日韩人妻高清精品专区| 2021少妇久久久久久久久久久| 欧美+日韩+精品| 国产欧美另类精品又又久久亚洲欧美| 内地一区二区视频在线| 欧美日韩视频精品一区| 国产片特级美女逼逼视频| 欧美激情极品国产一区二区三区 | 我的女老师完整版在线观看| 亚洲av不卡在线观看| 交换朋友夫妻互换小说| 久久国产亚洲av麻豆专区| 99久久精品一区二区三区| 亚洲美女视频黄频| 一级av片app| 美女脱内裤让男人舔精品视频| 久久综合国产亚洲精品| 欧美变态另类bdsm刘玥| 国产精品人妻久久久影院| 国产一级毛片在线| 国内少妇人妻偷人精品xxx网站| 18禁在线无遮挡免费观看视频| 欧美高清成人免费视频www| 成人国产av品久久久| 亚洲三级黄色毛片| 七月丁香在线播放| 成年av动漫网址| 日韩一本色道免费dvd| 亚洲国产最新在线播放| 九草在线视频观看| 观看美女的网站| 国产欧美日韩精品一区二区| 日日撸夜夜添| 三上悠亚av全集在线观看 | 久久久久久人妻| 亚洲av国产av综合av卡| 成人特级av手机在线观看| 三上悠亚av全集在线观看 | 高清毛片免费看| 这个男人来自地球电影免费观看 | 国产成人精品一,二区| 免费观看a级毛片全部| 水蜜桃什么品种好| 人妻夜夜爽99麻豆av| 99久久综合免费| 人人妻人人爽人人添夜夜欢视频 | 中文在线观看免费www的网站| 亚洲欧美精品专区久久| 精品一区二区三区视频在线| 波野结衣二区三区在线| 中文字幕av电影在线播放| 少妇的逼水好多| 各种免费的搞黄视频| 狠狠精品人妻久久久久久综合| 亚洲av国产av综合av卡| 国产黄色视频一区二区在线观看| 精品国产一区二区三区久久久樱花| 国产精品女同一区二区软件| 97超碰精品成人国产| 中文精品一卡2卡3卡4更新| 亚洲人成网站在线播| 亚洲精品,欧美精品| 男人爽女人下面视频在线观看| .国产精品久久| 日本av手机在线免费观看| av视频免费观看在线观看| 一本—道久久a久久精品蜜桃钙片| 亚洲av不卡在线观看| 亚洲真实伦在线观看| 中文乱码字字幕精品一区二区三区| 国产黄片美女视频| 国产成人免费无遮挡视频| 亚洲精品一二三| 欧美精品国产亚洲| 久久国产精品男人的天堂亚洲 | av福利片在线| 国产乱来视频区| 中文天堂在线官网| 精品一区二区三区视频在线| 插逼视频在线观看| 最后的刺客免费高清国语| www.av在线官网国产| 中国国产av一级| 久久这里有精品视频免费| 亚州av有码| 蜜桃久久精品国产亚洲av| 国产真实伦视频高清在线观看| 免费黄色在线免费观看| 久久久亚洲精品成人影院| 男人和女人高潮做爰伦理| 人体艺术视频欧美日本| 久久韩国三级中文字幕| 国产黄片美女视频| 嫩草影院入口| 成人亚洲欧美一区二区av| 一区二区三区四区激情视频| 国产精品99久久久久久久久| 观看美女的网站| kizo精华| 一级二级三级毛片免费看| 国产伦精品一区二区三区四那| 亚洲自偷自拍三级| 免费少妇av软件| 菩萨蛮人人尽说江南好唐韦庄| 一本一本综合久久| 人妻夜夜爽99麻豆av| 一本大道久久a久久精品| 国产高清不卡午夜福利| 欧美日韩国产mv在线观看视频| 国产精品.久久久| 天美传媒精品一区二区| 免费看光身美女| 欧美性感艳星| 欧美激情国产日韩精品一区| 天堂8中文在线网| 久久亚洲国产成人精品v| 丰满乱子伦码专区| kizo精华| 在线观看美女被高潮喷水网站| 人妻 亚洲 视频| 亚洲精品,欧美精品| 大香蕉97超碰在线| 80岁老熟妇乱子伦牲交| 精品久久久久久电影网| 亚洲综合精品二区| 熟女av电影| 99热国产这里只有精品6| 亚洲va在线va天堂va国产| 赤兔流量卡办理| 只有这里有精品99| 最后的刺客免费高清国语| 亚洲精华国产精华液的使用体验| 最黄视频免费看| 欧美日韩国产mv在线观看视频| 青春草视频在线免费观看| a级毛片在线看网站| 9色porny在线观看| 精品99又大又爽又粗少妇毛片| 欧美精品一区二区大全| 成人漫画全彩无遮挡| 乱人伦中国视频| 久久久亚洲精品成人影院| 国产成人免费无遮挡视频| 丝袜喷水一区| 免费人妻精品一区二区三区视频| 插阴视频在线观看视频| 久久毛片免费看一区二区三区| 国产精品秋霞免费鲁丝片| 狂野欧美白嫩少妇大欣赏| 黄色日韩在线| 精品久久国产蜜桃| 久久久久久人妻| 日韩av在线免费看完整版不卡| 熟女人妻精品中文字幕| 蜜臀久久99精品久久宅男| 人妻系列 视频| 久久99精品国语久久久| 亚洲国产欧美日韩在线播放 | 精品人妻一区二区三区麻豆| 国产一区亚洲一区在线观看| 日韩一区二区视频免费看| 亚洲欧洲日产国产| 熟女av电影| 一级毛片久久久久久久久女| 国产在线免费精品| 热re99久久国产66热| 亚洲成色77777| 午夜福利视频精品| 夜夜看夜夜爽夜夜摸| 高清欧美精品videossex| 制服丝袜香蕉在线| 国产黄色视频一区二区在线观看| 国产精品偷伦视频观看了| 亚洲av综合色区一区| 日本色播在线视频| 乱人伦中国视频| 十分钟在线观看高清视频www | 成人国产麻豆网| av线在线观看网站| 亚洲精品色激情综合| 99热这里只有是精品50| 免费av中文字幕在线| 97超碰精品成人国产| 一级毛片久久久久久久久女| 女性生殖器流出的白浆| 我要看黄色一级片免费的| 三级国产精品欧美在线观看| 国产成人aa在线观看| 国产精品国产三级专区第一集| 国产无遮挡羞羞视频在线观看| 精品久久久久久久久亚洲| 成年av动漫网址| 日本vs欧美在线观看视频 | 9色porny在线观看| 女人久久www免费人成看片| 王馨瑶露胸无遮挡在线观看| 亚洲av日韩在线播放| 国产精品99久久99久久久不卡 | 久久综合国产亚洲精品| 国产精品.久久久| 一区二区三区乱码不卡18| 免费观看av网站的网址| 51国产日韩欧美| 在线观看免费日韩欧美大片 | 黄色一级大片看看| 国产午夜精品久久久久久一区二区三区| 另类精品久久| 亚洲欧洲国产日韩| 日本91视频免费播放| 国产成人精品婷婷| 中国美白少妇内射xxxbb| 嫩草影院新地址| 成人免费观看视频高清| 人妻人人澡人人爽人人| 在线观看av片永久免费下载| 成人二区视频| 亚洲伊人久久精品综合| 精品亚洲成国产av| 色哟哟·www| a级片在线免费高清观看视频| 日本av手机在线免费观看| 人人妻人人澡人人爽人人夜夜| 一本大道久久a久久精品| 91久久精品电影网| 婷婷色av中文字幕| 国产一区二区在线观看av| 99热这里只有是精品50| 自拍欧美九色日韩亚洲蝌蚪91 | 美女福利国产在线| 国产免费福利视频在线观看| 国产视频内射| 日韩不卡一区二区三区视频在线| 国产精品麻豆人妻色哟哟久久| 国产精品欧美亚洲77777| 亚洲国产精品999| 国产免费视频播放在线视频| 香蕉精品网在线| 久久国产乱子免费精品| 涩涩av久久男人的天堂| 18禁裸乳无遮挡动漫免费视频| 青春草视频在线免费观看| 在线观看国产h片| 国产在线免费精品| 国内精品宾馆在线| 亚洲精品视频女| 嘟嘟电影网在线观看| 久久久久久伊人网av| 日韩精品免费视频一区二区三区 | 日本与韩国留学比较| 99热这里只有精品一区| 欧美精品国产亚洲| 亚洲久久久国产精品| 久久 成人 亚洲| 肉色欧美久久久久久久蜜桃| 麻豆乱淫一区二区| 精品国产一区二区久久| 精品一品国产午夜福利视频| 亚洲精品久久午夜乱码| 精品国产一区二区三区久久久樱花| 简卡轻食公司| 婷婷色综合大香蕉| 欧美精品高潮呻吟av久久| 黑人猛操日本美女一级片| 欧美日本中文国产一区发布| 欧美 亚洲 国产 日韩一| 亚洲精品亚洲一区二区| 蜜桃久久精品国产亚洲av| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲成人手机| av网站免费在线观看视频| 亚洲精品456在线播放app| 99久久精品热视频| 成人午夜精彩视频在线观看| 人人澡人人妻人| av福利片在线观看| 水蜜桃什么品种好| 久久精品久久久久久噜噜老黄| 成人美女网站在线观看视频| 在线观看美女被高潮喷水网站| 99精国产麻豆久久婷婷| 国产中年淑女户外野战色| av免费观看日本| 99热这里只有是精品50| 青青草视频在线视频观看| 精品国产一区二区三区久久久樱花| 熟女人妻精品中文字幕| 国产真实伦视频高清在线观看| 涩涩av久久男人的天堂| 我要看日韩黄色一级片| 久久国产精品大桥未久av | 99热国产这里只有精品6| 日韩视频在线欧美| 久久99热这里只频精品6学生| 在线天堂最新版资源| 午夜激情福利司机影院| 亚洲国产av新网站| 少妇被粗大的猛进出69影院 | 一区二区av电影网| a级一级毛片免费在线观看| 一区二区三区四区激情视频| 两个人的视频大全免费| 18禁动态无遮挡网站| 自拍欧美九色日韩亚洲蝌蚪91 | 大片免费播放器 马上看| 免费大片黄手机在线观看| 大片免费播放器 马上看| 成人亚洲欧美一区二区av| 国产老妇伦熟女老妇高清| 汤姆久久久久久久影院中文字幕| 欧美成人精品欧美一级黄| 美女国产视频在线观看| 一个人免费看片子| 国产亚洲欧美精品永久| 午夜免费男女啪啪视频观看| 久久久久久久精品精品| 看十八女毛片水多多多| 嫩草影院入口| 国产熟女欧美一区二区| 中文字幕精品免费在线观看视频 | 伊人久久国产一区二区| √禁漫天堂资源中文www| 看十八女毛片水多多多| 久久国内精品自在自线图片| 熟女电影av网| 久久国内精品自在自线图片| 亚洲欧美一区二区三区黑人 | 狂野欧美白嫩少妇大欣赏| 一区二区三区精品91| 国产日韩欧美亚洲二区| 久久人人爽av亚洲精品天堂| 在线天堂最新版资源| 另类精品久久| 午夜激情久久久久久久| 少妇高潮的动态图| 少妇人妻精品综合一区二区| 99热网站在线观看| 伦精品一区二区三区| 一本一本综合久久| 少妇丰满av| 国产在线男女| 午夜老司机福利剧场| 免费人妻精品一区二区三区视频| 少妇精品久久久久久久| 精品少妇久久久久久888优播| 国产成人精品无人区| 纯流量卡能插随身wifi吗| 中文精品一卡2卡3卡4更新| 免费看光身美女| a 毛片基地| 日本wwww免费看| a级片在线免费高清观看视频| 亚洲精品国产av成人精品| 国产 一区精品| 成年女人在线观看亚洲视频| 国产成人免费无遮挡视频| 日韩av在线免费看完整版不卡| av在线app专区| 观看美女的网站| 交换朋友夫妻互换小说| 亚洲av在线观看美女高潮| 99九九线精品视频在线观看视频| 永久网站在线| 免费黄网站久久成人精品| 亚洲欧洲精品一区二区精品久久久 | 色视频www国产| 亚洲av福利一区| 亚洲中文av在线| 免费黄色在线免费观看| 亚洲一级一片aⅴ在线观看| 国产在线视频一区二区| 狂野欧美激情性bbbbbb| 日韩大片免费观看网站| 看免费成人av毛片| 中文乱码字字幕精品一区二区三区| 亚洲欧美成人综合另类久久久| 欧美老熟妇乱子伦牲交| 老司机影院成人| 99久久精品热视频| 久久久久久久久久成人| a 毛片基地| 日韩大片免费观看网站| 黄色配什么色好看| 免费在线观看成人毛片| 少妇精品久久久久久久| 亚洲色图综合在线观看| 亚洲精品日韩在线中文字幕| tube8黄色片| 国产精品女同一区二区软件| 亚洲精品aⅴ在线观看| 国产在线一区二区三区精| 欧美日韩视频高清一区二区三区二| 91久久精品电影网| 新久久久久国产一级毛片| 啦啦啦中文免费视频观看日本| 亚洲一区二区三区欧美精品| 午夜影院在线不卡| 精品少妇内射三级|