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

    大型結(jié)冰風(fēng)洞氣流場適航符合性驗證

    2021-05-04 03:27:24郭向東柳慶林賴慶仁楊升科
    空氣動力學(xué)學(xué)報 2021年2期
    關(guān)鍵詞:馬赫數(shù)試驗段風(fēng)洞

    郭向東,柳慶林,賴慶仁,楊升科,趙 照

    (中國空氣動力研究與發(fā)展中心 結(jié)冰與防除冰重點實驗室,綿陽 621000)

    0 引 言

    飛機結(jié)冰廣泛存在于飛行實踐中,并嚴(yán)重威脅飛行安全[1-2]。民用航空適航規(guī)章要求制造商通過結(jié)冰防護等手段,確保飛機在結(jié)冰氣象條件下的飛行安全。結(jié)冰風(fēng)洞試驗[3],即利用結(jié)冰風(fēng)洞在地面模擬飛機結(jié)冰過程驗證飛機防除冰系統(tǒng)性能,是目前重要的飛機結(jié)冰適航審定手段,已經(jīng)廣泛應(yīng)用于飛機適航審定實踐中。為滿足結(jié)冰風(fēng)洞適航應(yīng)用要求[4],結(jié)冰風(fēng)洞必須開展全面的流場品質(zhì)評估,以驗證其適航符合性,其中氣流場品質(zhì)是重要的評估驗證內(nèi)容。

    美國NASA Glenn IRT結(jié)冰風(fēng)洞在該領(lǐng)域最早開展研究[5],建立了系統(tǒng)的氣流場校測方法,該風(fēng)洞圍繞1999年和2011年開展的兩次大型設(shè)備升級改造[6-11],相繼開展了多期氣流場校測,采用風(fēng)洞試驗和數(shù)值計算兩種手段,全面評估了改造后的風(fēng)洞穩(wěn)定段和試驗段內(nèi)流場品質(zhì),支撐了該風(fēng)洞流場的適航符合性驗證。意大利CIRA IWT結(jié)冰風(fēng)洞作為目前國際上尺寸最大、性能最完善的結(jié)冰風(fēng)洞之一[12],同樣發(fā)展了配套的氣流場校測設(shè)備和方法,開展了全面的氣流場校測,為該風(fēng)洞適航應(yīng)用奠定了基礎(chǔ)[13]。此外,美國Cox結(jié)冰風(fēng)洞于2007年也開展了氣流場校測[14],評估了氣流場品質(zhì),進而為風(fēng)洞升級改造提供了數(shù)據(jù)支撐。但是,我國在大型結(jié)冰風(fēng)洞氣流場適航符合性驗證方面缺乏系統(tǒng)的研究,僅針對小尺寸風(fēng)洞,開展了部分儀器應(yīng)用研究[15]。近年來,隨著國內(nèi)大型結(jié)冰風(fēng)洞3 m × 2 m結(jié)冰風(fēng)洞的建成,其適航符合性亟待得到驗證,其中氣流場是符合性驗證的關(guān)鍵一環(huán)。

    因此,鑒于目前國內(nèi)研究現(xiàn)狀以及急迫的型號試驗需求,本文首先發(fā)展結(jié)冰風(fēng)洞氣流場適航符合性驗證方法,然后針對3 m × 2 m結(jié)冰風(fēng)洞主試驗段構(gòu)型,開展氣流場適航符合性驗證試驗,考察試驗段氣流速度和噴嘴干空氣射流對流場特征參數(shù)(氣流速度、氣流偏角和氣流湍流度)的影響,最后評估試驗段內(nèi)氣流場空間均勻性和時間穩(wěn)定性,獲得試驗段氣流總壓、靜壓和馬赫數(shù)修正關(guān)系,形成結(jié)冰風(fēng)洞氣流場控制包線,為3 m × 2 m結(jié)冰風(fēng)洞適航應(yīng)用奠定基礎(chǔ)。

    1 3 m × 2 m結(jié)冰風(fēng)洞簡介

    中國空氣動力研究與發(fā)展中心3 m × 2 m結(jié)冰風(fēng)洞是一座閉口回流式高亞聲速風(fēng)洞(圖1),主要包括結(jié)冰噴霧系統(tǒng)、制冷系統(tǒng)、高度模擬系統(tǒng)和風(fēng)機動力系統(tǒng)。結(jié)冰噴霧系統(tǒng)利用噴霧耙和噴嘴產(chǎn)生結(jié)冰云霧,其中噴霧耙安裝于穩(wěn)定段蜂窩器下游、收縮段入口位置,共有20排,每排設(shè)置了50個噴嘴安裝位置,共計1 000個噴嘴。該風(fēng)洞擁有三個可更換的試驗段:主試驗段、次試驗段和高速試驗段,三個試驗段收縮比分別為14.67、5.73和29.33。本文僅針對主試驗段構(gòu)型開展研究,主試驗段尺寸為6.5 m(長) × 3 m(寬) × 2 m(高)。

    圖1 3 m × 2 m結(jié)冰風(fēng)洞Fig. 1 Configuration of the icing wind tunnel at CARDC

    結(jié)冰風(fēng)洞總壓和靜壓探針如圖2所示??傡o壓探針為自研設(shè)備,均具有加熱功能,用于結(jié)冰風(fēng)洞動力系統(tǒng)氣流速度控制,其中總壓探針位于噴霧耙與蜂窩器之間,共計4根,分別安裝于各洞壁水平中心線處,靜壓探針位于試驗段入口,在左右洞壁各布置2根探針。

    2 結(jié)冰風(fēng)洞氣流場適航符合性驗證方法

    結(jié)冰風(fēng)洞氣流場適航符合性驗證方法流程圖如圖3所示,該方法主要包括試驗標(biāo)準(zhǔn)、試驗內(nèi)容、試驗儀器、試驗方法和數(shù)據(jù)處理五部分內(nèi)容,本節(jié)將分別進行詳細(xì)論述。

    圖3 結(jié)冰風(fēng)洞氣流場適航符合性驗證方法流程圖Fig. 3 Flow diagram of the airworthiness compliance verification method for the aerodynamic flowfield of icing wind tunnels

    2.1 試驗標(biāo)準(zhǔn)及內(nèi)容

    《Calibration and Acceptance of Icing Wind Tunnels》(SAE ARP5905)為國際結(jié)冰適航領(lǐng)域公認(rèn)的結(jié)冰風(fēng)洞流場校測和符合性標(biāo)準(zhǔn)[16],是結(jié)冰風(fēng)洞氣流場適航符合性驗證的基礎(chǔ)。該標(biāo)準(zhǔn)給出了結(jié)冰風(fēng)洞氣流場品質(zhì)指標(biāo),如表1所示。表中針對氣流速度、氣流偏角和湍流度三個流場特征參數(shù),分別給出了測試設(shè)備最大不確定度、空間均勻性、風(fēng)洞中心線時間穩(wěn)定性和最大限制范圍四個典型參數(shù)的符合性指標(biāo)。因此,本文依據(jù)SAE ARP5905標(biāo)準(zhǔn),以氣流速度、氣流偏角和湍流度三個流場特征參數(shù)為試驗對象,考察試驗段氣流速度和噴嘴干空氣射流對流場參數(shù)的影響,評估其空間均勻性和時間穩(wěn)定性,進而驗證氣流場適航符合性。

    表1 結(jié)冰風(fēng)洞氣流場品質(zhì)指標(biāo)Table 1 Quality indexes of the air flowfield of icing wind tunnels

    2.2 試驗儀器

    五孔氣流方向探針和五孔探針排管架裝置用于測量試驗段內(nèi)氣流總靜壓和方向角,如圖4所示。

    五孔氣流方向探針前端為半球形,總長370 mm,外徑Φ10 mm。探針共有13個測壓孔(見圖5):總壓測壓孔位于探針前端中心;4個氣流偏角測壓孔位于半球形前端45°位置處,并沿周向等間距分布;8個靜壓測壓孔位于距探針頭部45 mm處,并沿周向等間距分布。試驗前已經(jīng)對五孔探針進行了標(biāo)定[17],獲得了不同馬赫數(shù)下的探針校準(zhǔn)系數(shù)和測量精度,其中速度測量精度為0.5%,氣流偏角測量精度為±0.1°。

    五孔探針排管架裝置尺寸為厚25 mm、寬200 mm、高2 000 mm,其前緣設(shè)置9個探針安裝孔,安裝孔間距為200 mm。排管架豎直安裝于試驗段內(nèi),上下洞壁處設(shè)置了11個橫向安裝位置,位置間隔為230 mm。

    熱線風(fēng)速儀和熱線排管架裝置用于測量試驗段內(nèi)氣流湍流度,如圖6所示。

    圖4 五孔氣流方向探針和五孔探針排管架Fig. 4 Five-hole flow angle probe and five-hole vertical rake

    圖5 熱線風(fēng)速儀和熱線排管架Fig. 5 Hot wire anemometer and hot wire vertical rake

    圖6 五孔探針排管架測量點位置矩陣Fig. 6 Location matrix of measuring points of the five-hole vertical rake

    熱線風(fēng)速儀由熱線探針和數(shù)據(jù)分析機組成,其中熱線探針外徑Φ4 mm、長250 mm,其前部安裝熱線探頭,尾部通過BNC線連接至數(shù)據(jù)分析機。試驗前對熱線風(fēng)速儀進行了校準(zhǔn),校準(zhǔn)結(jié)果顯示其速度測量精度和湍流度測量精度分別為1%和0.1%。應(yīng)該指出的是,熱線探頭可以實現(xiàn)多維度速度矢量的測量,本期試驗僅采用一維氣流速度測量探頭。

    熱線探針排管架裝置尺寸為厚21 mm、寬250 mm、高2 000 mm,其前緣設(shè)置7個探針安裝孔,安裝孔間距為250 mm。排管架豎直安裝于試驗段內(nèi),上下洞壁處設(shè)置11個橫向安裝位置(與五孔探針排管架裝置共用安裝板),位置間隔為230 mm。

    2.3 試驗方法及數(shù)據(jù)處理

    2.3.1 氣流速度和氣流偏角

    表2給出了氣流速度和氣流偏角試驗工況,其中排管架測點位置如圖7所示。試驗段名義氣流速度(VTS)選取40、80、100、120、140、160 m/s,覆蓋主試驗段主要試驗速度范圍。噴嘴噴氣壓力(pa)選取0、0.4、0.9 MPa,其中0 MPa和0.9 MPa分別對應(yīng)噴嘴噴氣關(guān)閉工況和噴嘴最大噴氣壓力工況。

    表2 氣流速度和氣流偏角試驗工況Table 2 Test conditions for airspeeds and airflow angularities

    圖7 五孔氣流方向探針測壓孔位置示意圖Fig. 7 Schematics of pressure tap positions of the five-hole probe

    圖7所示的排管架測量點位置矩陣,其中x-y平面為風(fēng)洞試驗段測試截面,CP為試驗段中心位置,x軸從試驗段左壁(沿流向左側(cè)為左壁)指向右壁,y軸從試驗段下壁面指向上壁面。試驗段內(nèi)共設(shè)置11列測量位置(C1~C11),每列包括9個測量點,共計99個測量位置,測量點橫向(x方向)和縱向(y方向)間距分別為230 mm和200 mm,分別覆蓋77%橫向和80%縱向試驗段區(qū)域。

    試驗時,首先測量試驗段中心處測量點(位置對應(yīng)圖7中C6列),然后沿x軸依次測量其余10個位置測量點(位置對應(yīng)圖7中C1~C5 和 C7 ~C11列),各測量點壓力參數(shù)采樣時間為60 s,采樣頻率為1 Hz。應(yīng)該指出的是,在每個橫向位置(x軸方向)壓力測量前,應(yīng)先測量排管架的安裝俯仰角和偏航角,獲得氣流偏角修正參數(shù)。

    試驗后利用測量的壓力參數(shù),采用五孔探針計算公式,計算當(dāng)?shù)貧饬骺傡o壓、馬赫數(shù)和氣流偏角。根據(jù)五孔探針測壓孔位置示意圖(見圖5),當(dāng)?shù)貧饬骺倝海╬t,Local)、靜壓(ps,Local)和馬赫數(shù)(MaLocal)計算式表示為:

    式中,Maprobe為探針計算馬赫數(shù),該參數(shù)經(jīng)過修正獲得當(dāng)?shù)貧饬黢R赫數(shù);C(Ma)為馬赫數(shù)修正系數(shù),通過探針標(biāo)定獲得。

    當(dāng)?shù)貧饬髌前ǜ┭鼋牵é粒┖推浇牵é拢?,計算公式為?/p>

    式中,θ、Cα和Cβ分別為當(dāng)?shù)貧饬髌墙^對角度、俯仰角壓力系數(shù)和偏航角壓力系數(shù);K0,α、K1,α、K0,β和K2,β分別為俯仰角和偏航角系數(shù),由探針標(biāo)定獲得;Δα和Δβ為探針安裝偏角,試驗前測得。

    氣流速度和氣流偏角空間分布以試驗段中心線處流場參數(shù)為基準(zhǔn),采用馬赫數(shù)比空間偏差(ΔMaR)和氣流偏角空間偏差(Δα,Δβ)表征,表示為:

    其中下標(biāo)C表示試驗段中心線參數(shù),MaR、MaTS和MaWT分別為馬赫數(shù)比、試驗段馬赫數(shù)和風(fēng)洞馬赫數(shù)。進而采用馬赫數(shù)比空間偏差標(biāo)準(zhǔn)差(σ(ΔMaR))和最大絕對值(|ΔMaR|max)評估試驗段內(nèi)氣流速度均勻性;同理,采用氣流偏角空間偏差標(biāo)準(zhǔn)差(σ(Δα),σ(Δβ))和最大絕對值(|Δα|max,|Δβ|max)評估氣流偏角均勻性。

    氣流速度時間分布以試驗段中心線處馬赫數(shù)時間平均值((MaC)mt)為基準(zhǔn),采用馬赫數(shù)時間偏差(ΔMaT)表征,表示為:

    式中下標(biāo)mt表示時間平均。進而采用馬赫數(shù)時間偏差標(biāo)準(zhǔn)差(σ(ΔMaT))和最大絕對值(|ΔMaT|max)評估試驗段中心線處氣流速度時間穩(wěn)定性。

    最后,利用試驗段氣流總壓(pt,TS)、靜壓(ps,TS)和馬赫數(shù)(MaTS)與風(fēng)洞測量的氣流總壓(pt,WT)、靜壓(ps,WT)和馬赫數(shù)(MaWT),獲得試驗段中心線處總壓、靜壓和馬赫數(shù)修正關(guān)系,表示為:

    式中,Kpt、Kps和KMa分別試驗段中心線處總壓、靜壓和馬赫數(shù)修正函數(shù)。

    2.3.2 氣流湍流度

    氣流湍流度試驗工況如表3所示,其中測量點位置如圖8所示;試驗段名義氣流速度選取40、60、80、100、120、140 m/s,其中120 m/s和140 m/s工況僅測量試驗段中心線處(P5位置)的氣流湍流度;噴嘴噴氣壓力(pa)選取0 MPa和0.4 MPa,其中0.4 MPa為常用試驗條件下噴嘴最大噴氣壓力。

    表3 氣流湍流度評估工況Table 3 Test conditions for the turbulence intensity

    圖8 熱線排管架測量點位置矩陣Fig. 8 Location matrix of measuring points of the hot wire probe rake

    熱線探針排管架測量點位置矩陣如圖8所示,圖中坐標(biāo)系統(tǒng)與圖6相同。試驗段內(nèi)共設(shè)置11列測量位置,每列包括7個測量位置(對應(yīng)探針安裝位置),共計77個測量位置(如圖灰色點所示),測量點間橫向(x方向)和縱向(y方向)間距分別為230 mm和250 mm。實際試驗中,考慮到熱線探頭易損壞等因素,僅測量試驗段模型區(qū)內(nèi)P1~P9九個位置的氣流湍流度(如圖紅色點所示),其中P5為測試截面中心點位置。

    試驗前,應(yīng)在試驗段內(nèi)對熱線風(fēng)速儀進行標(biāo)定,標(biāo)定風(fēng)速范圍覆蓋試驗范圍(見表3)。試驗時,為降低熱線探頭損壞的風(fēng)險,僅進行單點測試,首先測量截面中心位置測量點(位置對應(yīng)圖8中P5點),然后依次測量其余8個位置(P1~P4 和 P6~P9),各測量點參數(shù)采樣時間為5 s,采樣頻率為3 000 Hz。

    試驗后利用測量的風(fēng)速參數(shù),計算當(dāng)?shù)貧饬魍牧鞫龋═I),表達(dá)式為

    式中,Ui、Urms和N分別為氣流速度(沿流向一維速度)、平均速度和速度采樣個數(shù)。

    氣流湍流度空間分布和時間分布以試驗段中心線處氣流湍流度(TITS,C)為基準(zhǔn),采用湍流度空間偏差(ΔTIS)和時間偏差(ΔTIT)表征,表示為

    式中TITS為試驗段氣流湍流度。進而采用湍流度空間偏差和時間偏差標(biāo)準(zhǔn)差(σ(ΔTIS)和σ(ΔTIT))和最大絕對值(|ΔTIS|max和|ΔTIT|max)評估試驗段內(nèi)氣流湍流度均勻性和試驗段中心處時間穩(wěn)定性。

    3 試驗結(jié)果與分析討論

    3.1 氣流速度

    首先,為評估試驗段氣流速度空間均勻性,圖9給出了試驗段內(nèi)馬赫數(shù)空間分布云圖,圖中包括0 MPa和0.4 MPa兩種噴嘴氣壓條件下40 m/s、80 m/s和160 m/s三個典型速度對應(yīng)的試驗結(jié)果,紅色虛線框表示模型區(qū),模型區(qū)范圍為-920 mm≤x≤920 mm,-600 mm≤y≤600 mm。從圖中可以看出:80 m/s和160 m/s工況下,馬赫數(shù)分布一致性較好,其中模型區(qū)內(nèi)馬赫數(shù)比均在±1%范圍內(nèi);但是在40 m/s工況下,馬赫數(shù)空間分布差異較大,試驗段內(nèi)存在偏差超過2%的非均勻峰值區(qū)域,并且這些非均勻區(qū)域主要集中在C3(-670 mm)和C9(670 mm)位置,對應(yīng)噴霧耙豎直支撐排架位置(見圖1(b))。由此可見,噴霧耙結(jié)構(gòu)會影響試驗段內(nèi)馬赫數(shù)空間分布形態(tài),進而導(dǎo)致40 m/s工況下形成了非均勻的氣流速度峰值區(qū)。

    圖9 試驗段馬赫數(shù)空間分布云圖Fig. 9 Spatial distribution of the Mach number contour in the test section

    進一步,為定量評估試驗段氣流速度空間均勻性,圖10給出了試驗段模型區(qū)內(nèi)馬赫數(shù)比空間標(biāo)準(zhǔn)差和最大絕對值。從圖中可以看出:隨著氣流速度的增加,馬赫數(shù)比空間偏差標(biāo)準(zhǔn)差和最大絕對值不斷減小,氣流速度空間均勻性不斷增強,尤其當(dāng)氣流速度超過80 m/s后,標(biāo)準(zhǔn)差和最大絕對值分布近似為0.25%和1%,但在40 m/s工況下分別近似為1.2%和4%,此時模型區(qū)內(nèi)出現(xiàn)超過2%的非均勻峰值區(qū);噴嘴干空氣射流對氣流速度均勻性并無顯著影響,參數(shù)變化幅度小于0.5%。

    然后,為評估試驗段中心線處氣流速度時間穩(wěn)定性,圖11給出了試驗段中心線處馬赫數(shù)時間偏差變化曲線,圖中包括0 MPa、0.4 MPa和0.9 MPa三種噴氣狀態(tài)下20 m/s、40 m/s、80 m/s和160 m/s四個速度對應(yīng)的試驗結(jié)果,從圖中可以看出:各試驗工況下,試驗段中心線處馬赫數(shù)時間偏差均在±2%范圍內(nèi),其中20 m/s和40 m/s條件下氣流馬赫數(shù)出現(xiàn)顯著的脈動,噴嘴氣壓越高,脈動越強烈,氣流穩(wěn)定性越弱,而80 m/s以上時,氣流馬赫數(shù)較穩(wěn)定,偏差在±0.2%范圍內(nèi)。

    圖10 試驗段模型區(qū)內(nèi)馬赫數(shù)空間比偏差標(biāo)準(zhǔn)差和最大絕對值Fig. 10 Standard deviation and maximum absolute value of the Mach number ratio spatial deviation in the model area of the test section

    圖11 試驗段中心線處馬赫數(shù)時間偏差變化曲線Fig. 11 Time variation of the Mach number time deviation in the centerline of the test section

    進一步,圖12給出了試驗段中心線處馬赫數(shù)時間偏差標(biāo)準(zhǔn)差和最大絕對值,從圖中可以看出:隨著氣流速度的增加,馬赫數(shù)時間偏差標(biāo)準(zhǔn)差和最大絕對值不斷減小,氣流速度時間穩(wěn)定性不斷增強,整體而言,標(biāo)準(zhǔn)差和最大絕對值分別均小于0.8%和2%,尤其當(dāng)氣流速度超過60 m/s后,標(biāo)準(zhǔn)差和最大絕對值分別均小于0.2%和0.4%;增大噴嘴氣壓會增大馬赫數(shù)時間偏差標(biāo)準(zhǔn)差和最大絕對值,但增大幅度會隨著氣流速度的增加而不斷減小,其中對于20 m/s工況,標(biāo)準(zhǔn)差和最大絕對值的最大變化幅度約為0.4%和1%,但當(dāng)氣流速度超過60 m/s后,變化幅度均小于0.2%,由此可見噴嘴干空氣射流僅對60 m/s以下試驗工況的氣流速度時間穩(wěn)定性產(chǎn)生較顯著影響,并且隨著噴嘴氣壓的增大,氣流速度時間穩(wěn)定性不斷減弱。

    最后,試驗段中心線處總壓、靜壓和馬赫數(shù)修正關(guān)系如圖13所示,圖中包括噴嘴氣壓為0 MPa、0.4 MPa和0.9 MPa對應(yīng)的試驗結(jié)果、擬合曲線和擬合曲線公式。從圖中可以看出,風(fēng)洞測量的氣流總壓、靜壓和馬赫數(shù)與試驗段中心線處氣流總壓、靜壓和馬赫數(shù)具有顯著的線性關(guān)系,擬合公式表示為:

    進一步,圖14給出了各工況下試驗段中心線處氣流馬赫數(shù)比,從圖中可以看出:隨著氣流速度的增加,馬赫數(shù)比逐漸增大,但是隨著噴嘴氣壓的增大,馬赫數(shù)比的增大趨勢不斷減弱,尤其當(dāng)氣壓為0.9 MPa時,各氣流速度下馬赫數(shù)比近似為0.995;增大噴嘴氣壓,會增大馬赫數(shù)比,但增大幅度會隨著氣流速度的增加而不斷減小,其中在20 m/s工況下馬赫數(shù)比最大偏差(對應(yīng)0.9 MPa與0 MPa馬赫數(shù)比之差)接近7%,但當(dāng)氣流速度超過60 m/s后,馬赫數(shù)比最大偏差將小于1%。這主要是因為噴嘴會向試驗氣流中注入高速空氣射流,這些高速空氣射流會將自身攜帶的動能傳遞給試驗氣流,進而增大試驗段氣流速度,并且氣壓越大,則射流流量和動能越大,試驗段氣流速度增大幅度越大;但是隨著試驗段氣流速度的增加,試驗氣流流量和動能不斷增大,此時噴嘴干空氣射流對試驗氣流的影響不斷減弱,速度增加幅度則不斷減小。

    圖12 試驗段中心線處馬赫數(shù)時間偏差標(biāo)準(zhǔn)差和最大絕對值Fig. 12 Standard deviation and maximum absolute value of the Mach number time deviation in the centerline of the test section

    圖13 試驗段中心線處氣流總壓、靜壓和馬赫數(shù)修正關(guān)系Fig. 13 Correction relationships of the total pressure, the static pressure and the Mach number in the centerline of the test section

    圖14 試驗段中心線處馬赫數(shù)比Fig. 14 Mach number ratio in the centerline of the test section

    3.2 氣流偏角

    為評估試驗段內(nèi)氣流偏角空間均勻性,圖15給出了試驗段氣流偏角空間分布云圖,圖中包括噴嘴不噴氣條件下40 m/s、80 m/s和160 m/s三個典型速度對應(yīng)的試驗結(jié)果,分布結(jié)果以氣流偏角絕對角度θ表示,圖中箭頭方向表示當(dāng)?shù)貧饬髌窃跍y量平面內(nèi)的投影,紅色虛線框表示模型區(qū),范圍為-920 mm≤x≤920 mm,-600 mm≤y≤600 mm。從圖中可以看出:各工況下氣流偏角方向分布形態(tài)相似,呈現(xiàn)上下區(qū)域向中心收縮、左右區(qū)域向外側(cè)發(fā)散的形態(tài);各工況下普遍存在偏角超過1°的氣流偏角非均勻峰值區(qū),并且這些區(qū)域主要集中在C3(-670 mm)和C9(670 mm)位置,對應(yīng)噴霧耙豎直支撐排架位置(見圖1(b)),這與氣流速度空間分布特征相一致。由此可見,噴霧耙結(jié)構(gòu)同樣會影響氣流偏角空間分布形態(tài),進而導(dǎo)致試驗段內(nèi)出現(xiàn)非均勻的氣流偏角峰值區(qū)。

    進一步,為定量評估試驗段內(nèi)氣流偏角空間均勻性,圖16給出了試驗段模型區(qū)內(nèi)氣流偏角(俯仰角和偏航角)最大絕對值、空間偏差標(biāo)準(zhǔn)差和空間偏差最大絕對值。從圖中可以看出:試驗段氣流速度對俯仰角和偏航角最大絕對值并無顯著影響,其中俯仰角最大絕對值均小于1.5°,而偏航角最大絕對值在氣流速度80 m/s以上均小于1.5°,但在40 m/s工況下,由于受到噴嘴干空氣射流的顯著影響,偏航角最大絕對值會隨著噴嘴氣壓不斷增大,最大值約為3°;試驗段氣流速度對俯仰角和偏航角空間偏差標(biāo)準(zhǔn)差也無顯著影響,表明氣流速度對氣流偏角空間均勻性的影響不顯著,其中俯仰角空間偏差標(biāo)準(zhǔn)差均小于0.5°,而偏航角空間偏差標(biāo)準(zhǔn)差在80 m/s以上小于0.6°,但對于40 m/s工況,由于受到噴嘴干空氣射流的顯著影響,標(biāo)準(zhǔn)差會隨著噴嘴氣壓不斷增大,最大值約為1.1°;俯仰角空間偏差最大絕對值均小于1.2°,而偏航角空間偏差最大絕對值在80 m/s以上小于2°,但對于40 m/s工況,在噴嘴干空氣射流的顯著影響下,最大絕對值會隨著噴嘴氣壓不斷增大,進而導(dǎo)致模型區(qū)內(nèi)出現(xiàn)氣流偏角超過2°的非均勻峰值區(qū)。

    圖15 試驗段氣流偏角空間分布云圖Fig. 15 Spatial distribution of the airflow angularity in the test section

    圖16 試驗段模型區(qū)內(nèi)氣流偏角(俯仰角和偏航角)最大絕對值、空間偏差標(biāo)準(zhǔn)差和最大絕對值Fig. 16 The maximum absolute value, standard deviation and its maximum absolute value of the airflow angularities(pitch angle and yaw angle) in the model area of the test section

    3.3 氣流湍流度

    為評估試驗段內(nèi)氣流湍流度空間均勻性,圖17給出了試驗段氣流湍流度空間分布曲線,圖中橫軸P1~P9位置如圖8所示,其中P5對應(yīng)試驗段中心線位置,圖中包括0 MPa和0.4 MPa兩種噴氣狀態(tài)下40 、60、80、100 m/s四個速度對應(yīng)的試驗結(jié)果,同時給出了平均氣流湍流度(TIm)。從圖中可以看出:各試驗工況下,試驗段氣流湍流度均小于1%;噴嘴干空氣射流會顯著提高氣流湍流度,尤其對于40 m/s工況,平均氣流湍流度增大幅度約0.44%。

    圖17 試驗段氣流湍流度空間分布曲線Fig. 17 Spatial distribution of the airflow turbulence intensity in the test section

    進一步,為定量評估試驗段湍流度空間均勻性,圖18給出了試驗段內(nèi)湍流度空間偏差標(biāo)準(zhǔn)差和最大絕對值,從圖中可以看出:湍流度空間偏差標(biāo)準(zhǔn)差和最大絕對值隨試驗段氣流速度并無明顯變化,表明試驗段氣流速度對氣流湍流度空間均勻性影響不顯著,整體而言,湍流度空間偏差標(biāo)準(zhǔn)差和最大絕對值分別小于0.2%和0.5%;噴嘴干空氣射流對40 m/s工況下氣流湍流度均勻性影響較大,其中標(biāo)準(zhǔn)差和最大絕對值增幅接近0.15%和0.5%。

    圖18 試驗段湍流度空間偏差標(biāo)準(zhǔn)差和最大絕對值Fig. 18 The standard deviation and maximum absolute value of the turbulence intensity spatial deviation in the test section

    為評估試驗段中心線處氣流湍流度時間穩(wěn)定性,圖19給出了試驗段中心線處湍流度時間偏差變化曲線,圖中包括0 MPa和0.4 MPa兩種噴氣狀態(tài)下60 m/s、80 m/s、100 m/s和140 m/s四個速度對應(yīng)的試驗結(jié)果,從圖中可以看出:各試驗工況下,試驗段中心線處湍流度時間偏差均在±0.2%范圍內(nèi);噴嘴氣壓對湍流度時間偏差影響不顯著,影響幅度小于0.1%。

    圖20給出了試驗段中心線處湍流度時間偏差標(biāo)準(zhǔn)差和最大絕對值,從圖中可以看出:湍流度時間偏差標(biāo)準(zhǔn)差和最大絕對值隨試驗段氣流速度并無明顯變化,表明試驗段氣流速度對氣流湍流度時間均勻性影響不顯著,整體而言,湍流度時間偏差標(biāo)準(zhǔn)差和最大絕對值分別小于0.1%和0.2%;噴嘴干空氣射流會增大標(biāo)準(zhǔn)差和最大絕對值,進而減弱氣流湍流度時間穩(wěn)定性,其中標(biāo)準(zhǔn)差和最大絕對值最大增幅接近0.08%和0.15%。

    3.4 結(jié)冰風(fēng)洞氣流場控制包線

    根據(jù)3 m×2 m結(jié)冰風(fēng)洞主試驗段氣流場品質(zhì)評估結(jié)果,圖21給出了3 m×2 m結(jié)冰風(fēng)洞主試驗段氣流場控制包線,圖中黑色點為試驗工況點,紅線區(qū)域為控制包線,該包線針對氣流速度、氣流偏角和氣流湍流度三個流場特征參數(shù),綜合考慮流場空間均勻性和時間穩(wěn)定性品質(zhì),給出了滿足適航審定要求的試驗氣流速度范圍。從圖中可以看出,3 m×2 m結(jié)冰風(fēng)洞主試驗段氣流場品質(zhì)在主要試驗速度范圍內(nèi)均滿足適航審定要求,其中試驗速度在80 m/s至140 m/s間的氣流場品質(zhì)最優(yōu)(如圖藍(lán)色線所示)。

    圖19 試驗段中心線處湍流度時間偏差變化曲線Fig. 19 Time variation of the turbulence intensity time deviation in the centerline of the test section

    圖20 試驗段中心線處湍流度時間偏差標(biāo)準(zhǔn)差和最大絕對值Fig. 20 Standard deviation and maximum absolute value of the turbulence intensity time deviation in the centerline of the test section

    圖21 3 m×2 m結(jié)冰風(fēng)洞主試驗段氣流場控制包線Fig. 21 Operating envelop of the aerodynamic flowfield in the main test section of the icing wind tunnel at CARDC

    4 結(jié) 論

    本文發(fā)展了結(jié)冰風(fēng)洞氣流場適航符合性驗證方法,開展了3 m × 2 m結(jié)冰風(fēng)洞主試驗段氣流場適航符合性驗證試驗,評估了其氣流場品質(zhì),主要得到以下結(jié)論:

    1)噴霧耙結(jié)構(gòu)會影響試驗段內(nèi)氣流速度和氣流偏角空間分布形態(tài),進而導(dǎo)致了非均勻峰值區(qū)的形成。

    2)增大試驗段氣流速度會增強試驗段氣流速度空間均勻性和時間穩(wěn)定性,但噴嘴干空氣射流會減弱氣流偏角和氣流湍流度空間均勻性,同時降低氣流速度和氣流湍流度時間穩(wěn)定性,尤其對試驗段氣流速度低于60 m/s的氣流場品質(zhì)影響最為顯著。

    3) 3 m × 2 m結(jié)冰風(fēng)洞主試驗段氣流場品質(zhì)在主要試驗速度范圍內(nèi)滿足適航審定要求,其中80 m/s至140 m/s范圍內(nèi)的氣流場品質(zhì)最優(yōu)。

    致謝:感謝中國商飛上海飛機設(shè)計研究院李海星博士對本文試驗方法和數(shù)據(jù)處理部分提出的寶貴建議。

    猜你喜歡
    馬赫數(shù)試驗段風(fēng)洞
    高馬赫數(shù)激波作用下單模界面的Richtmyer-Meshkov不穩(wěn)定性數(shù)值模擬
    爆炸與沖擊(2024年7期)2024-11-01 00:00:00
    一維非等熵可壓縮微極流體的低馬赫數(shù)極限
    跨聲速風(fēng)洞槽壁試驗段流場品質(zhì)提升措施研究
    載荷分布對可控擴散葉型性能的影響
    斑頭雁進風(fēng)洞
    合肥軌道交通1號線試驗段Ⅰ標(biāo)基坑監(jiān)測及結(jié)果分析研究
    黃風(fēng)洞貂鼠精
    基于NI cRIO平臺的脈沖燃燒風(fēng)洞控制系統(tǒng)設(shè)計
    HL-2M真空室試驗段制造工藝技術(shù)
    中國核電(2017年2期)2017-08-11 08:00:54
    能源監(jiān)控管理系統(tǒng)在滬寧高速公路照明試驗段的應(yīng)用
    亚洲精品成人av观看孕妇| 国产精品久久久久久精品电影小说| 成人黄色视频免费在线看| 男女免费视频国产| 女性生殖器流出的白浆| 国产精品久久久av美女十八| 乱人伦中国视频| 亚洲精品国产区一区二| 日日摸夜夜添夜夜添小说| 久久久久精品国产欧美久久久| 国产精品亚洲av一区麻豆| 人人妻人人澡人人看| 一边摸一边抽搐一进一小说 | 亚洲专区中文字幕在线| 在线十欧美十亚洲十日本专区| 色尼玛亚洲综合影院| 午夜福利在线观看吧| 女人久久www免费人成看片| 捣出白浆h1v1| 国产高清激情床上av| 国产精品麻豆人妻色哟哟久久| 十八禁高潮呻吟视频| 曰老女人黄片| 91精品国产国语对白视频| 亚洲精品久久午夜乱码| 久久国产精品人妻蜜桃| av天堂久久9| 国产高清国产精品国产三级| 男女高潮啪啪啪动态图| 成人国语在线视频| 日韩一卡2卡3卡4卡2021年| tocl精华| 欧美精品一区二区免费开放| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品自拍成人| 夜夜爽天天搞| 一区二区三区精品91| 又大又爽又粗| 如日韩欧美国产精品一区二区三区| 黄片大片在线免费观看| 男女免费视频国产| 肉色欧美久久久久久久蜜桃| 国产精品久久久久久人妻精品电影 | 91成人精品电影| 午夜福利欧美成人| 色在线成人网| 美女视频免费永久观看网站| 精品一区二区三区视频在线观看免费 | 老司机亚洲免费影院| 国产精品国产高清国产av | 在线观看免费视频日本深夜| 国产老妇伦熟女老妇高清| 在线 av 中文字幕| 国产精品1区2区在线观看. | 伦理电影免费视频| 国产精品 欧美亚洲| 精品第一国产精品| 黄色毛片三级朝国网站| 成人精品一区二区免费| 精品高清国产在线一区| 中文字幕av电影在线播放| 中文字幕av电影在线播放| 精品高清国产在线一区| 成人精品一区二区免费| 亚洲欧美一区二区三区久久| 日韩免费高清中文字幕av| 欧美日韩中文字幕国产精品一区二区三区 | av一本久久久久| 国产高清videossex| 人人澡人人妻人| 国产成人av激情在线播放| 亚洲精品av麻豆狂野| 香蕉久久夜色| 777久久人妻少妇嫩草av网站| 欧美激情久久久久久爽电影 | 一级a爱视频在线免费观看| 美女国产高潮福利片在线看| 一二三四在线观看免费中文在| 亚洲熟女毛片儿| 中文亚洲av片在线观看爽 | 一级a爱视频在线免费观看| 亚洲欧美一区二区三区久久| 超碰成人久久| 亚洲专区国产一区二区| 天堂俺去俺来也www色官网| 午夜精品国产一区二区电影| 久久久久久久国产电影| 老熟女久久久| 亚洲美女黄片视频| 亚洲国产毛片av蜜桃av| 亚洲欧美色中文字幕在线| 男女无遮挡免费网站观看| 午夜福利一区二区在线看| 国产一卡二卡三卡精品| 日本黄色日本黄色录像| 美女高潮喷水抽搐中文字幕| 欧美黄色片欧美黄色片| 亚洲中文av在线| 亚洲精品国产精品久久久不卡| 久久久水蜜桃国产精品网| 日韩视频一区二区在线观看| 国产精品av久久久久免费| 亚洲国产精品一区二区三区在线| 狠狠婷婷综合久久久久久88av| 在线看a的网站| 国产福利在线免费观看视频| 亚洲国产欧美一区二区综合| 国产一区二区三区在线臀色熟女 | 国产一区二区在线观看av| 成在线人永久免费视频| 国产亚洲欧美在线一区二区| 丁香六月天网| 国产精品影院久久| 人人妻,人人澡人人爽秒播| 日本五十路高清| 精品国产一区二区三区四区第35| 久久狼人影院| 久久久精品免费免费高清| 久久久精品免费免费高清| 一区二区日韩欧美中文字幕| 国产淫语在线视频| 欧美变态另类bdsm刘玥| 午夜福利,免费看| 老司机影院毛片| 肉色欧美久久久久久久蜜桃| 在线看a的网站| 欧美日韩成人在线一区二区| 99热国产这里只有精品6| 久久这里只有精品19| 欧美黄色淫秽网站| 夜夜夜夜夜久久久久| 成人三级做爰电影| 欧美精品av麻豆av| 一级毛片精品| 99九九在线精品视频| 中文字幕av电影在线播放| 精品久久久精品久久久| 色婷婷久久久亚洲欧美| 黑人巨大精品欧美一区二区mp4| 日韩一区二区三区影片| 丁香六月天网| 欧美激情 高清一区二区三区| 肉色欧美久久久久久久蜜桃| 大香蕉久久成人网| 一区二区av电影网| 欧美一级毛片孕妇| 亚洲欧洲日产国产| 高潮久久久久久久久久久不卡| 国产成人av激情在线播放| 女人爽到高潮嗷嗷叫在线视频| av免费在线观看网站| 国产高清视频在线播放一区| 精品亚洲成国产av| 精品一品国产午夜福利视频| 欧美激情极品国产一区二区三区| 交换朋友夫妻互换小说| 国产精品偷伦视频观看了| 日本wwww免费看| 亚洲精品国产色婷婷电影| 日日爽夜夜爽网站| 欧美精品啪啪一区二区三区| 亚洲精品国产区一区二| 免费黄频网站在线观看国产| 搡老岳熟女国产| 深夜精品福利| 久久99一区二区三区| 国产精品成人在线| 日本wwww免费看| 一级毛片电影观看| 日韩免费高清中文字幕av| 另类精品久久| 免费黄频网站在线观看国产| 深夜精品福利| 久久中文字幕人妻熟女| 麻豆av在线久日| 国产精品秋霞免费鲁丝片| 两人在一起打扑克的视频| svipshipincom国产片| 美女福利国产在线| 黄色视频在线播放观看不卡| 91字幕亚洲| 亚洲成人国产一区在线观看| 大陆偷拍与自拍| 91九色精品人成在线观看| 黑人巨大精品欧美一区二区mp4| 日本五十路高清| 99在线人妻在线中文字幕 | 国产亚洲精品久久久久5区| 一区二区三区精品91| e午夜精品久久久久久久| 丰满迷人的少妇在线观看| 免费在线观看视频国产中文字幕亚洲| 婷婷成人精品国产| 国产又爽黄色视频| 日韩三级视频一区二区三区| 最近最新免费中文字幕在线| 在线av久久热| 我的亚洲天堂| 老司机在亚洲福利影院| 十八禁高潮呻吟视频| www.999成人在线观看| cao死你这个sao货| 99久久精品国产亚洲精品| 制服人妻中文乱码| 欧美日韩国产mv在线观看视频| 丝袜在线中文字幕| 国产又爽黄色视频| 黑丝袜美女国产一区| 两性夫妻黄色片| av网站在线播放免费| 色综合婷婷激情| 欧美成人午夜精品| 国产真人三级小视频在线观看| 午夜福利,免费看| 俄罗斯特黄特色一大片| 亚洲av国产av综合av卡| 欧美黄色淫秽网站| 日日夜夜操网爽| 水蜜桃什么品种好| 亚洲第一av免费看| 午夜福利在线免费观看网站| 桃红色精品国产亚洲av| 国产亚洲午夜精品一区二区久久| 丰满少妇做爰视频| 女人被躁到高潮嗷嗷叫费观| 久久精品aⅴ一区二区三区四区| 国产精品av久久久久免费| 99香蕉大伊视频| 亚洲情色 制服丝袜| av免费在线观看网站| 91成人精品电影| 久久久国产精品麻豆| 欧美日韩成人在线一区二区| 国产男女内射视频| 成人永久免费在线观看视频 | 久久久欧美国产精品| 国产又色又爽无遮挡免费看| 亚洲男人天堂网一区| 久久亚洲真实| 婷婷丁香在线五月| 国产亚洲午夜精品一区二区久久| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲精品第一综合不卡| 日本黄色日本黄色录像| 国产精品 国内视频| 成人三级做爰电影| 老司机亚洲免费影院| 无遮挡黄片免费观看| 国产亚洲欧美在线一区二区| 一边摸一边抽搐一进一小说 | 精品久久蜜臀av无| 一级毛片精品| 国产精品久久久久久精品古装| 人妻一区二区av| 亚洲av欧美aⅴ国产| 精品欧美一区二区三区在线| 在线观看免费高清a一片| 午夜视频精品福利| 十八禁网站网址无遮挡| 99riav亚洲国产免费| 国产片内射在线| 热99国产精品久久久久久7| 国产亚洲精品第一综合不卡| 日韩熟女老妇一区二区性免费视频| 欧美日韩一级在线毛片| 黄色a级毛片大全视频| 亚洲精品成人av观看孕妇| av在线播放免费不卡| 成人亚洲精品一区在线观看| 国产欧美日韩综合在线一区二区| 最近最新中文字幕大全电影3 | 久久毛片免费看一区二区三区| 国产欧美日韩精品亚洲av| 欧美亚洲 丝袜 人妻 在线| 亚洲精品国产色婷婷电影| 国产精品99久久99久久久不卡| 丁香欧美五月| 国产亚洲精品久久久久5区| 性少妇av在线| 国产免费视频播放在线视频| 欧美日韩亚洲高清精品| 99re6热这里在线精品视频| 久9热在线精品视频| xxxhd国产人妻xxx| 亚洲精品国产区一区二| 国产在线视频一区二区| 在线观看免费午夜福利视频| 看免费av毛片| 水蜜桃什么品种好| 午夜福利欧美成人| 满18在线观看网站| 亚洲精品乱久久久久久| 一区二区三区乱码不卡18| 免费看a级黄色片| 一夜夜www| 自拍欧美九色日韩亚洲蝌蚪91| www日本在线高清视频| 国产xxxxx性猛交| 国产伦人伦偷精品视频| 精品国产一区二区久久| 丁香六月欧美| 黑人巨大精品欧美一区二区蜜桃| 国产一区二区三区综合在线观看| 亚洲色图 男人天堂 中文字幕| 精品国产乱子伦一区二区三区| 999精品在线视频| 十八禁网站网址无遮挡| 成人黄色视频免费在线看| 看免费av毛片| 久久久久久人人人人人| 久久精品人人爽人人爽视色| 老司机午夜十八禁免费视频| 成人精品一区二区免费| 午夜福利视频在线观看免费| 欧美日韩亚洲综合一区二区三区_| 精品一区二区三区四区五区乱码| 亚洲国产欧美网| 国产精品麻豆人妻色哟哟久久| 男女下面插进去视频免费观看| 男女高潮啪啪啪动态图| 欧美日韩黄片免| 女性生殖器流出的白浆| 99久久精品国产亚洲精品| 丰满少妇做爰视频| 亚洲国产毛片av蜜桃av| 亚洲av第一区精品v没综合| 亚洲成国产人片在线观看| 在线十欧美十亚洲十日本专区| av在线播放免费不卡| 亚洲一区二区三区欧美精品| 国产片内射在线| 国产一区二区三区视频了| 90打野战视频偷拍视频| 欧美日韩福利视频一区二区| 99精品在免费线老司机午夜| 日韩制服丝袜自拍偷拍| 日韩视频一区二区在线观看| 欧美黄色淫秽网站| 老熟妇乱子伦视频在线观看| 1024香蕉在线观看| 中文欧美无线码| 国产1区2区3区精品| 丁香六月天网| 人成视频在线观看免费观看| 成年女人毛片免费观看观看9 | 高清黄色对白视频在线免费看| 成人18禁在线播放| 国产老妇伦熟女老妇高清| 国产亚洲一区二区精品| 国产一区二区三区综合在线观看| 99国产精品一区二区蜜桃av | av一本久久久久| 另类精品久久| 女性生殖器流出的白浆| 日韩中文字幕视频在线看片| 精品福利观看| 十八禁网站免费在线| 无人区码免费观看不卡 | 黑丝袜美女国产一区| aaaaa片日本免费| 久久久久网色| 精品乱码久久久久久99久播| 十八禁人妻一区二区| 美女福利国产在线| 自线自在国产av| 亚洲精品粉嫩美女一区| 一边摸一边抽搐一进一小说 | 淫妇啪啪啪对白视频| 另类精品久久| 黄色视频,在线免费观看| 国产激情久久老熟女| 一区二区三区激情视频| av有码第一页| 成人国语在线视频| 在线观看免费视频日本深夜| 亚洲精品粉嫩美女一区| 亚洲成a人片在线一区二区| 老熟妇乱子伦视频在线观看| 午夜福利,免费看| 欧美黄色淫秽网站| 亚洲专区国产一区二区| 国产午夜精品久久久久久| 777米奇影视久久| 蜜桃国产av成人99| 国产av精品麻豆| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品在线美女| 成在线人永久免费视频| www.精华液| 亚洲精品中文字幕在线视频| 精品国产一区二区三区久久久樱花| av网站在线播放免费| 国产精品1区2区在线观看. | 亚洲中文字幕日韩| 国产三级黄色录像| 可以免费在线观看a视频的电影网站| 激情在线观看视频在线高清 | 日韩视频在线欧美| 久久精品aⅴ一区二区三区四区| 国产一区二区三区综合在线观看| 视频区图区小说| 制服人妻中文乱码| 欧美另类亚洲清纯唯美| e午夜精品久久久久久久| 一进一出抽搐动态| 久久香蕉激情| 纵有疾风起免费观看全集完整版| 国产一区二区三区在线臀色熟女 | 亚洲专区字幕在线| 午夜福利乱码中文字幕| 无人区码免费观看不卡 | 一夜夜www| 精品国内亚洲2022精品成人 | 国产老妇伦熟女老妇高清| 国产欧美亚洲国产| 亚洲国产欧美在线一区| 激情视频va一区二区三区| 性高湖久久久久久久久免费观看| 精品一区二区三区四区五区乱码| 国产精品免费视频内射| 国产亚洲欧美在线一区二区| 久久ye,这里只有精品| 丝瓜视频免费看黄片| 色精品久久人妻99蜜桃| 淫妇啪啪啪对白视频| 午夜福利视频在线观看免费| a在线观看视频网站| 黑人操中国人逼视频| 少妇被粗大的猛进出69影院| 1024视频免费在线观看| 天堂8中文在线网| 十八禁高潮呻吟视频| 国产成人影院久久av| av电影中文网址| 欧美老熟妇乱子伦牲交| 极品人妻少妇av视频| 黄色成人免费大全| 国产免费现黄频在线看| 男女高潮啪啪啪动态图| 久久精品国产a三级三级三级| 大香蕉久久成人网| 国产成人精品久久二区二区免费| av视频免费观看在线观看| 精品国产乱码久久久久久小说| 欧美日韩av久久| 成人黄色视频免费在线看| 777米奇影视久久| 成人av一区二区三区在线看| 国产成人精品久久二区二区91| 一级毛片女人18水好多| 精品国产一区二区三区四区第35| 免费观看a级毛片全部| 欧美人与性动交α欧美软件| 黄色视频不卡| www.熟女人妻精品国产| 午夜福利视频在线观看免费| 久9热在线精品视频| 亚洲精品一卡2卡三卡4卡5卡| 国产精品一区二区在线观看99| 国产片内射在线| 欧美精品高潮呻吟av久久| 操出白浆在线播放| 欧美 亚洲 国产 日韩一| 99精品欧美一区二区三区四区| 欧美一级毛片孕妇| 欧美 日韩 精品 国产| 欧美精品人与动牲交sv欧美| 欧美乱码精品一区二区三区| 精品人妻熟女毛片av久久网站| 熟女少妇亚洲综合色aaa.| 成人国产一区最新在线观看| 欧美乱妇无乱码| 成年人午夜在线观看视频| 91成年电影在线观看| 精品第一国产精品| 精品一区二区三卡| 亚洲国产av新网站| 激情视频va一区二区三区| 母亲3免费完整高清在线观看| 啦啦啦中文免费视频观看日本| 午夜福利一区二区在线看| 国产区一区二久久| 一边摸一边做爽爽视频免费| 狠狠狠狠99中文字幕| 亚洲国产欧美网| 美女高潮到喷水免费观看| 久久国产精品人妻蜜桃| 国产不卡一卡二| 亚洲精品乱久久久久久| 亚洲av成人不卡在线观看播放网| 亚洲国产欧美一区二区综合| netflix在线观看网站| 国产欧美日韩精品亚洲av| 丰满人妻熟妇乱又伦精品不卡| 免费少妇av软件| av片东京热男人的天堂| 国产精品久久久久久精品古装| 视频在线观看一区二区三区| 无限看片的www在线观看| 亚洲精品美女久久久久99蜜臀| 正在播放国产对白刺激| 国产又爽黄色视频| 国产一区二区在线观看av| 12—13女人毛片做爰片一| 波多野结衣一区麻豆| 一区二区三区激情视频| 中国美女看黄片| a级毛片黄视频| 中文字幕人妻丝袜制服| 亚洲av片天天在线观看| 亚洲久久久国产精品| 搡老乐熟女国产| 午夜福利乱码中文字幕| 久久精品国产99精品国产亚洲性色 | 汤姆久久久久久久影院中文字幕| 99国产精品一区二区三区| 免费高清在线观看日韩| 日日摸夜夜添夜夜添小说| a级片在线免费高清观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 十八禁高潮呻吟视频| 99精品久久久久人妻精品| 久久婷婷成人综合色麻豆| 成人手机av| 一区二区三区激情视频| 男女边摸边吃奶| 黄色视频在线播放观看不卡| 免费观看av网站的网址| 女人爽到高潮嗷嗷叫在线视频| 午夜福利在线免费观看网站| av国产精品久久久久影院| 天堂动漫精品| 久久午夜综合久久蜜桃| 日韩欧美一区视频在线观看| 中文欧美无线码| 无限看片的www在线观看| 热99久久久久精品小说推荐| 欧美日韩视频精品一区| 女人爽到高潮嗷嗷叫在线视频| 80岁老熟妇乱子伦牲交| 国产成人一区二区三区免费视频网站| 老熟妇乱子伦视频在线观看| 欧美激情极品国产一区二区三区| av网站在线播放免费| 在线亚洲精品国产二区图片欧美| 亚洲精品中文字幕一二三四区 | 999久久久国产精品视频| 亚洲精品中文字幕一二三四区 | 在线观看人妻少妇| 久久精品91无色码中文字幕| 丁香六月欧美| 日韩欧美一区视频在线观看| 黄色片一级片一级黄色片| 国产精品一区二区精品视频观看| 久久久水蜜桃国产精品网| 777米奇影视久久| 亚洲全国av大片| 1024视频免费在线观看| 久久国产精品影院| 汤姆久久久久久久影院中文字幕| 在线观看免费视频日本深夜| 久久精品成人免费网站| 亚洲精品成人av观看孕妇| 国产精品偷伦视频观看了| 国产无遮挡羞羞视频在线观看| 青草久久国产| 亚洲全国av大片| 免费不卡黄色视频| 美女午夜性视频免费| 多毛熟女@视频| 18禁国产床啪视频网站| 国产人伦9x9x在线观看| 国产91精品成人一区二区三区 | 考比视频在线观看| 成人18禁在线播放| 久久天堂一区二区三区四区| 精品第一国产精品| www.999成人在线观看| 久久青草综合色| 亚洲av电影在线进入| 高清欧美精品videossex| 久久久久久久精品吃奶| 欧美精品亚洲一区二区| 免费人妻精品一区二区三区视频| 热re99久久精品国产66热6| 麻豆国产av国片精品| 亚洲国产成人一精品久久久| av福利片在线| www.自偷自拍.com| 操美女的视频在线观看| 成人av一区二区三区在线看| 久久午夜亚洲精品久久| 精品国产亚洲在线| 精品免费久久久久久久清纯 | 欧美久久黑人一区二区| a级片在线免费高清观看视频| 中文字幕最新亚洲高清| 国产精品国产av在线观看| 女人精品久久久久毛片| 一个人免费看片子| 日本wwww免费看| 午夜两性在线视频| 午夜福利在线观看吧| 欧美日韩国产mv在线观看视频| 极品教师在线免费播放| 亚洲欧美一区二区三区久久| 免费观看av网站的网址| 丝袜美足系列| 99精品欧美一区二区三区四区| 国产精品一区二区免费欧美| 中文字幕人妻丝袜一区二区| 久久热在线av| 中文字幕人妻丝袜一区二区| 国产av国产精品国产| 国产在线视频一区二区|