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

    溝槽表面對翼型繞流流場的氣動作用研究

    2021-11-20 14:38:40陳子龍楊愛玲陳二云
    能源研究與信息 2021年3期
    關(guān)鍵詞:升力溝槽脈動

    陳子龍,楊愛玲,陳二云,焦 躍

    (上海理工大學(xué) 能源與動力工程學(xué)院/上海市動力工程多相流動與傳熱重點實驗室,上海 200093)

    葉片表面流動狀態(tài)是決定旋轉(zhuǎn)葉輪機械氣動性能與噪聲特性的關(guān)鍵因素。由于氣體的黏性,氣體流經(jīng)葉片表面時會形成緊貼壁面的流體薄層。該薄層厚度小、黏性影響大,即所謂的流動邊界層。流動邊界層內(nèi)速度梯度大,壁面黏性切應(yīng)力大,而且邊界層內(nèi)渦流在葉片尾緣周期性地脫落,在某些條件下邊界層還可能發(fā)生流動分離。這些正是葉片表面流動阻力增加、流動噪聲源增強的根本原因。

    為了提高葉片的氣動特性及降噪性能,國內(nèi)外學(xué)者從葉片結(jié)構(gòu)優(yōu)化、表面設(shè)計等方面做了大量的研究工作,試圖從設(shè)計源頭改善葉片的氣動降噪性能。非光滑表面技術(shù)作為一門新興的邊緣學(xué)科,涉及仿生學(xué)、流體力學(xué)、聲學(xué)等眾多學(xué)科。目前,較為成熟的研究方法有渦流發(fā)生器[1-2]、疏水表面[3]、鋸齒尾緣[4]、仿生非光滑表面[5-6]等。

    流動邊界層控制技術(shù)被美國NASA 研究中心列為21 世紀(jì)航空關(guān)鍵技術(shù)之一[7]。溝槽非光滑表面技術(shù)是流動邊界層控制技術(shù)的重要分支。該方法通過改變物體表面形狀,控制壁面湍流附面層的應(yīng)力分布與渦流結(jié)構(gòu),達到流動控制的目的。NACA Langley 研究中心Walsh[8]最先開展了將微溝槽非光滑表面應(yīng)用于平板減阻的研究,突破了表面越光滑阻力越小的傳統(tǒng)思維方式,從而引起了各國研究者對非光滑表面減阻研究的興趣。Sundaram 等[9]對來流攻角為0°~6°的NACA012翼型的溝槽減阻研究結(jié)果表明,減阻效果隨來流攻角的增加而增加,最大減阻效果達13%??罩锌蛙嚬驹贏320 試驗機機翼的70%表面貼上溝槽面薄膜,獲得了1%~2%的節(jié)油效果[10]。楊林等[11]基于“凹槽”和“陷窩”技術(shù)對低雷諾數(shù)下渦輪流動損失控制進行了實驗研究與數(shù)值分析,研究表明:凹槽結(jié)構(gòu)的存在能加強流動摻混,促進分離泡轉(zhuǎn)捩,有效降低葉柵總壓損失系數(shù)。

    目前相關(guān)學(xué)者們對溝槽非光滑表面技術(shù)的研究更多關(guān)注其減阻效果與工程應(yīng)用,并取得了大量研究成果[12-14]。而關(guān)于溝槽內(nèi)渦流結(jié)構(gòu)及與下游湍流流動相互影響的研究相對較少。因此,對溝槽非光滑表面技術(shù)的湍流特性進行研究非常有必要。

    本文基于NACA6510 翼型表面的二維圓弧型溝槽為研究對象,通過數(shù)值方法研究溝槽相對位置及相對深度的變化對翼型表面湍流流動的影響,為溝槽非光滑表面技術(shù)在流動噪聲控制應(yīng)用提供理論基礎(chǔ)。

    1 計算模型和數(shù)值計算方法

    1.1 幾何模型和網(wǎng)格劃分

    圖1 為NACA6510 翼型表面二維圓弧型溝槽幾何示意圖,其中:翼型弦長c=100 mm,溝槽位于翼型吸力面;x表示弦向坐標(biāo);λ、h分別表示溝槽的流向長度和深度。本文改變溝槽相對深度h/λ及溝槽在弦長方向的相對位置x/c獲得15 組不同的翼型模型,溝槽幾何參數(shù)如表1 所示。通過二維非定常湍流流場模擬分析上述參數(shù)對溝槽內(nèi)外流動、翼型表面湍流的影響規(guī)律。

    圖1 NACA6510 翼型圓弧型溝槽幾何示意圖Fig.1 Sketch of the arc-shaped groove on NACA6510 airfoil

    表1 溝槽幾何參數(shù)Tab.1 Geometric parameters of the groove

    計算時采用二維結(jié)構(gòu)網(wǎng)格。為節(jié)省計算資源,將計算域劃分為A、B 兩個子域,數(shù)值計算域示意圖如圖2 所示。子域A 為外流場,長20c、寬16c,保證遠場邊界條件。內(nèi)部子域B 長3.5c、寬2c、尾跡2c。子域B 采用O 型拓撲結(jié)構(gòu)進行網(wǎng)格離散,兩個子域交接面的網(wǎng)格一一對應(yīng)。溝槽翼型網(wǎng)格數(shù)為140 萬~160 萬,其中子域B 網(wǎng)格數(shù)為110 萬~125 萬。

    圖2 數(shù)值計算域示意圖Fig.2 Schematic diagram of the computation field

    溝槽幾何尺寸較小,為提高網(wǎng)格精度,對翼型近壁區(qū)、溝槽附近進行局部加密,近壁第一層網(wǎng)格高度為0.01 mm,確保Y+≤1,溝槽內(nèi)部采用C 型拓撲。溝槽表面翼型網(wǎng)格示意圖如圖3 所示。

    圖3 溝槽表面翼型網(wǎng)格示意圖Fig.3 Mesh of the airfoil with groove

    1.2 控制方程及邊界條件

    運動馬赫數(shù)Ma<0.3的空氣繞翼型流動,形成繞翼型的不可壓湍流流場。本文采用二維不可壓縮大渦模擬方法(large eddy simulation,LED)求解該翼型外部湍流場。大渦模擬方法求解瞬時N?S 方程獲得大于計算網(wǎng)格尺度的湍流運動,而小尺度渦對大尺度渦的影響則通過亞格子應(yīng)力模型體現(xiàn)。對于變量 φ,空間濾波后變?yōu)?/p>

    式中:D為流動區(qū)域;G(x,x′)為濾波函數(shù)。

    式中,V為控制體體積。

    將式(2)代入式(1)可得

    利用式(3)處理連續(xù)性方程及瞬時N?S 方程,可得二維不可壓縮大渦模擬運動方程[15],即

    式中:ρ、μ 分別為空氣密度、分子黏度;u為速度分量;p為修正壓力;τij表 示亞格子應(yīng)力,τij=,通過它可以體現(xiàn)小尺度渦對大渦運動方程的影響;帶有上劃線的為濾波后的變量;t為時間。

    為使運動方程封閉,采用Smagorinsky?Lilly亞格子模型,即

    式中:μt為亞格子湍流黏性系數(shù);Δ 為網(wǎng)格尺度;Cs為 亞格子常數(shù);Sij為網(wǎng)格平均應(yīng)變率。

    流動方程采用二階迎風(fēng)格式離散,基于SIMPLEC 算法實現(xiàn)速度場和壓力場分離迭代求解。計算域進口給定速度邊界條件,來流速度V=30m·s?1,基于翼型來流速度與弦長雷諾數(shù)Re=2×105。計算域出口位于翼型下游11倍弦長處,采用自由出流條件,翼型物面及溝槽采用無滑移絕熱條件。時間步長 Δt=1×10?5s,監(jiān)測翼型升力系數(shù)呈明顯周期性變化后,推進10 000步作為非定常流場信息。大渦模擬計算時,采用SSTk?ω湍流模型求解RANS 方程,將時均收斂解作為非定常流場計算的初場。

    2 計算結(jié)果與分析

    2.1 計算方法驗證

    為驗證數(shù)值計算結(jié)果的準(zhǔn)確性,計算驗證時采用Selig等[16]的經(jīng)典實驗,實驗?zāi)P蜑镋387翼型,弦長為0.304 8 m,基于弦長的雷諾數(shù)Re=2×105。驗證時采用與本文相同的網(wǎng)格拓撲結(jié)構(gòu)及計算方法。圖4 為E387 翼型極曲線試驗數(shù)據(jù)與本文模擬值的對比,圖中橫坐標(biāo)Cd為阻力系數(shù),縱坐標(biāo)Cl為升力系數(shù)??梢姡瑪?shù)值計算結(jié)果與實驗結(jié)果保持一致,表明本文數(shù)值計算方法可行。

    圖4 E387 翼型極曲線Fig.4 Polar curve of E387 airfoil

    2.2 翼型氣動性能及繞流流場特征

    極曲線反映了翼型的升、阻力特性,是翼型重要的氣動數(shù)據(jù)之一[16]。圖5為雷諾數(shù)Re=2×105,來流攻角α分別為0°、3°、6°、···、18°時NACA6510 翼型極曲線,其中:圖5(a)給出了相對深度h/λ=0.35 時翼型極曲線隨溝槽相對位置x/c的變化;圖5(b)則為x/c=0.65 時溝槽深度對翼型極曲線的影響。由圖5(a)可知:小攻角下溝槽相對位置對翼型升、阻力幾乎沒有影響;隨來流攻角增加,影響逐漸明顯,在α=18°時,溝槽使翼型阻力系數(shù)和升力系數(shù)減小,相對位置x/c=0.65 時,阻力系數(shù)從光滑翼型的0.226 5 減小為0.162 5,下降了28.26%。圖5(b)表明,大來流攻角下溝槽深度對翼型升、阻力有較明顯的影響。當(dāng)相對深度h/λ=0.35 時,阻力系數(shù)從光滑翼型的0.226 5 減小為0.169 1,減少了25.34%;升力系數(shù)從光滑翼型的1.743 9 減小為1.437 5,下降了17.57%,阻力系數(shù)的減少量大于升力系數(shù)的減少量。

    圖5 NACA6510 翼型極曲線Fig.5 Polar curves of NACA6510 airfoil

    圖6 給出了來流攻角α=6°、t=0.135~0.140 s、溝槽相對位置x/c=0.65 時NACA6510 翼型升力系數(shù)時域曲線隨相對深度h/λ的變化??梢?,升力系數(shù)隨時間呈明顯的周期性波動。根據(jù)流體力學(xué)理論可知,這正是翼型尾緣周期性渦流脫落引起的升力脈動。表2 為α=6°、x/c=0.65 時經(jīng)快速傅里葉變換(FFT)后獲得的翼型升力系數(shù)脈動主頻及對應(yīng)的幅值,在溝槽相對深度小于0.35時,脈動頻率幾乎不變,但對應(yīng)的峰值有變化。隨相對深度進一步增加,當(dāng)h/λ=0.50 時,升力系數(shù)脈動頻率出現(xiàn)突變,頻率明顯減小。

    圖6 升力系數(shù)時間歷程Fig.6 Time history of lift coefficient

    表2 不同相對深度溝槽升力系數(shù)脈動主頻及對應(yīng)的幅值Tab.2 Main fluctuation frequency and the corresponding amplitude of lift coefficient for the groove with various relative depth

    為觀測翼型近壁區(qū)流場壓強脈動,在翼型吸力面?zhèn)仍O(shè)置g、p、f三個監(jiān)測點,如圖7 所示。監(jiān)測點g位于0.3 倍弦長位置,離吸力壁面2 mm處;監(jiān)測點p在凹坑內(nèi)部監(jiān)測點f則位于0.8 倍弦長位置,離吸力壁面2 mm 處。

    圖7 監(jiān)測點示意圖Fig.7 Distribution of monitor points

    圖8 為溝槽相對深度h/λ=0.35 時監(jiān)測點f的壓強脈動在頻域的分布,表3 為α=6°、h/λ=0.35 時的主頻及對應(yīng)幅值。可見,該監(jiān)測點壓強脈動在頻域的分布呈明顯的離散特性,具有溝槽的翼型壓強脈動峰值均小于光滑翼型,且主脈動峰值隨溝槽相對位置增加而減小,當(dāng)x/c=0.75時,峰值由光滑翼型的225.11 Pa 減小為110.43 Pa,降幅為50.9%。從表3 還可知,相對位置x/c=0.55、0.75 時,壓強脈動頻率均有所降低。

    表3 不同相對位置溝槽監(jiān)測點f 壓強脈動主頻及對應(yīng)幅值Tab.3 Main frequency and the corresponding amplitude of pressure fluctuation at point f for the groove at different relative positions

    圖8 不同相對位置溝槽監(jiān)測點f 壓強脈動頻譜Fig.8 Pressure fluctuation spectrum of point f for the groove at different relative positions

    表4 為來流攻角α=6°、相對位置x/c=0.65時不同相對深度溝槽表面翼型及光滑翼型監(jiān)測點f的壓強脈動主頻及對應(yīng)幅值。由表4 可知,同一相對位置下,隨著溝槽深度的增加,監(jiān)測點f脈動幅值呈先降后升的分布趨勢,h/λ=0.35 時脈動幅值達到最小。表4 同時表明,溝槽較淺(h/λ=0.10、h/λ=0.20、h/λ=0.35)時,溝槽的存在對頻率的影響較?。粶喜圻_到一定深度(h/λ=0.50、h/λ=0.80)后,壓強脈動頻率發(fā)生突變,脈動頻率減小了110 Hz。監(jiān)測點f壓強脈動的頻譜特性與翼型升力系數(shù)(參見表2)的變化吻合。通過上文分析說明,翼型吸力面的溝槽結(jié)構(gòu)會改變下游渦流脫落頻率和幅值。

    表4 不同相對深度溝槽監(jiān)測點f 壓強脈動主頻及對應(yīng)幅值Tab.4 Main frequency and the corresponding amplitude of pressure fluctuation at point f for the groove at different relative depth

    2.3 溝槽內(nèi)流動特征分析

    氣體流經(jīng)翼型表面,部分流體進入溝槽,在溝槽內(nèi)部形成復(fù)雜的流動狀態(tài)。圖9 為來流攻角α=6°、相對位置x/c=0.65 時溝槽內(nèi)部時均速度及流線隨相對深度h/λ的變化。當(dāng)h/λ較小時,流體流過溝槽結(jié)構(gòu)區(qū)域,在內(nèi)部產(chǎn)生了一個反向、穩(wěn)定的扁平狀二次渦,流動速度較小。當(dāng)h/λ=0.35 時,有更多較高能量流體流入溝槽內(nèi)部,形成了一大一小兩個反向渦對。反向渦對的大小隨溝槽相對增加呈先增大后減小的趨勢。當(dāng)溝槽深度進一步增加到h/λ=0.8 時,高能量的大渦擠壓溝槽內(nèi)較小的渦,使之變小。

    圖9 溝槽結(jié)構(gòu)內(nèi)部時均速度及流線分布Fig.9 Distribution of mean velocity and streamline in the groove

    圖10(a)為來流攻角 α=6、x/c=0.65 時沿溝槽中心線l的時均速度分布,縱坐標(biāo)ΔY指中心線上的觀測點到光滑翼型表面的長度。相關(guān)定義如圖10(b)所示,l為過凹坑圓心的直線,ΔY<0表示監(jiān)測點位于溝槽內(nèi)部。由圖10 可知,隨壁面相對距離的增加,速度明顯呈先增大后減少再增大的趨勢,說明流體流經(jīng)溝槽結(jié)構(gòu)時在溝槽內(nèi)部形成了漩渦。深徑比對 ΔY>0區(qū)域的速度分布剖面影響較小,但溝槽內(nèi)部的速度分布發(fā)生明顯變化。相對于光滑壁面,溝槽相對深度增加,溝槽內(nèi)部速度增量越大,這與上文對溝槽內(nèi)部速度流線圖的分析相對應(yīng)。

    圖10 溝槽中心線時均速度分布Fig.10 Distribution of mean velocity along the center of the groove

    圖11 為α=6°、x/c=0.65、h/λ=0.35 時溝槽內(nèi)部瞬時流線在一個主脈動周期ΔT(ΔT為監(jiān)測點p壓強脈動主頻對應(yīng)的周期)內(nèi)分布。可以看出,在6°來流攻角下,附面層在溝槽上游已形成小分離渦,溝槽內(nèi)部的二次渦由于該上游分離渦的擠壓、附面層流體的裹挾共同作用經(jīng)歷了產(chǎn)生、發(fā)展、脫落過程,無法穩(wěn)定滯留在溝槽內(nèi)部。

    圖11 1 個周期內(nèi)各相位瞬時流線分布Fig.11 Distribution of instantaneous streamline during each phase of a period

    圖12 為α=6°、x/c=0.65、h/λ=0.35 時監(jiān)測點g、p、f的壓強脈動頻譜分布(監(jiān)測點位置如圖7 所示)。可見,在溝槽上游的監(jiān)測點g,壓強脈動只在1 600 Hz 附近有較小的峰值,監(jiān)測點的壓強脈動在3 200 Hz 處出現(xiàn)了大小約35 Pa的峰值。該頻率應(yīng)為溝槽內(nèi)部二次渦周期性產(chǎn)生、發(fā)展和脫落的頻率,而在溝槽下游監(jiān)測點f處,壓強脈動在1 600、3 200 Hz 處均有明顯的脈動峰,相較于g、p兩點,監(jiān)測點f的脈動幅值大幅增長。氣體流經(jīng)溝槽結(jié)構(gòu),翼型尾緣保持了上游及溝槽內(nèi)部壓強脈動的頻率特性,脈動幅值逐漸增強。

    圖12 溝槽上、下游壓強脈動頻譜Fig.12 Pressure fluctuation spectrum at the upstream and downstream of the groove

    表5 給出了α=6°、h/λ=0.35 時,不同相對位置溝槽上、下游監(jiān)測點g、p、f的壓強脈動主頻及對應(yīng)幅值。表6 則為α=6°、x/c=0.65 時,不同相對深度溝槽上、下游壓強脈動主頻及對應(yīng)幅值??芍煌瑓?shù)溝槽的存在對其翼型表面上、下游壓強脈動的影響,其頻率呈現(xiàn)類似圖12的特性,溝槽內(nèi)部監(jiān)測點p的壓強脈動主頻表現(xiàn)為上、下游主頻的2 倍,說明溝槽內(nèi)部二次渦產(chǎn)生、發(fā)展和脫落過程始終存在。溝槽結(jié)構(gòu)的存在對上游監(jiān)測點g的壓強脈動幅值幾乎沒有影響。在不同相對位置下,隨弦向的增加,溝槽內(nèi)部監(jiān)測點p脈動幅值逐漸增加,但溝槽下游監(jiān)測點f幅值卻呈減小趨勢。在同一位置,隨溝槽相對深度的增加,溝槽內(nèi)部監(jiān)測點p脈動幅值逐漸增加,下游監(jiān)測點f幅值呈先增大后減小趨勢,在h/λ=0.35 時出現(xiàn)最小值。

    表5 不同相對位置溝槽上、下游壓強脈動主頻及對應(yīng)幅值Tab.5 Main frequency and the corresponding amplitude of pressure fluctuationat the upstream and downstream of the groove at different relative positions

    表6 不同相對深度溝槽上、下游壓強脈動主頻及對應(yīng)幅值Tab.6 Main frequency and the corresponding amplitude of pressure fluctuationat the upstream and downstream of the groove at different relative depth

    3 結(jié) 論

    本文以NACA510 翼型溝槽結(jié)構(gòu)為研究對象,通過大渦模擬方法對翼型繞流流場進行了數(shù)值計算,分析了溝槽結(jié)構(gòu)對翼型氣動特性與渦流結(jié)構(gòu)的影響規(guī)律,取得了以下結(jié)論:

    (1)在小攻角下,溝槽表面對翼型的升、阻力系數(shù)影響很小,大攻角下升、阻力系數(shù)明顯減小,溝槽位于x/c=0.65、相對深度h/λ=0.35時,翼型阻力系數(shù)的減小量大于升力系數(shù)的減小量。

    (2)溝槽結(jié)構(gòu)的存在改變了下游渦脫落情況,對上游影響不明顯。6°來流攻角下,隨弦向方向,溝槽下游監(jiān)測點f壓強脈動幅值逐漸減??;隨溝槽深度的增加,監(jiān)測點壓強脈動幅值呈先降后增的分布趨勢,在h/λ=0.35 時達到最小。溝 槽 較 淺(h/λ=0.10、h/λ=0.20、h/λ=0.35)時,溝槽的存在對頻率的影響不明顯;溝槽達到一定深度(h/λ=0.50、h/λ=0.80)后,壓強脈動頻率發(fā)生突變,頻率明顯降低。

    (3)氣體流經(jīng)溝槽結(jié)構(gòu)時會在溝槽內(nèi)部形成反向二次渦,其壓強脈動主頻為上、下游主頻的2 倍。在一定條件下,溝槽內(nèi)部渦流失穩(wěn),存在周期性渦流產(chǎn)生、發(fā)展和脫離。氣體流經(jīng)溝槽結(jié)構(gòu),翼型尾緣保持了上游及溝槽內(nèi)部壓強脈動的頻率特性,脈動幅值逐漸增強。

    猜你喜歡
    升力溝槽脈動
    新學(xué)期,如何“脈動回來”?
    家教世界(2023年25期)2023-10-09 02:11:56
    高速列車車頂–升力翼組合體氣動特性
    RBI在超期服役脈動真空滅菌器定檢中的應(yīng)用
    一種具有多形式鋼片結(jié)構(gòu)的四季胎
    無人機升力測試裝置設(shè)計及誤差因素分析
    基于自適應(yīng)偽譜法的升力式飛行器火星進入段快速軌跡優(yōu)化
    一種低噪音的全路況輪胎
    地球脈動(第一季)
    升力式再入飛行器體襟翼姿態(tài)控制方法
    溝槽爆破參數(shù)優(yōu)化及成本分析
    av福利片在线| 亚洲av日韩精品久久久久久密 | 欧美日韩一级在线毛片| 日本vs欧美在线观看视频| 精品一区在线观看国产| 高清欧美精品videossex| 99热全是精品| 宅男免费午夜| 久久久久久久久久久久大奶| 亚洲欧洲日产国产| 五月天丁香电影| 国产精品久久久久久精品古装| 亚洲国产精品成人久久小说| 丝袜脚勾引网站| av不卡在线播放| 在线观看免费高清a一片| 成人国语在线视频| 中文字幕人妻熟女乱码| 波多野结衣一区麻豆| 日韩伦理黄色片| 桃花免费在线播放| 大片电影免费在线观看免费| 夫妻午夜视频| 一级a爱视频在线免费观看| 999久久久国产精品视频| 高清视频免费观看一区二区| 亚洲天堂av无毛| 国产人伦9x9x在线观看| 一二三四在线观看免费中文在| 国产一区有黄有色的免费视频| 国产在视频线精品| 国产精品99久久99久久久不卡| 91成人精品电影| 久久精品人人爽人人爽视色| 一区二区av电影网| 一级毛片我不卡| 伊人久久大香线蕉亚洲五| 在线天堂中文资源库| 欧美亚洲 丝袜 人妻 在线| 免费av中文字幕在线| 亚洲国产精品成人久久小说| 久久亚洲精品不卡| 男的添女的下面高潮视频| 国产男人的电影天堂91| 曰老女人黄片| 1024香蕉在线观看| 中文字幕人妻熟女乱码| 天天躁日日躁夜夜躁夜夜| 91字幕亚洲| 国产欧美亚洲国产| 日韩一卡2卡3卡4卡2021年| 啦啦啦视频在线资源免费观看| 欧美日韩黄片免| 尾随美女入室| 又大又黄又爽视频免费| av电影中文网址| 国产日韩欧美亚洲二区| 亚洲欧美一区二区三区黑人| 大码成人一级视频| 精品久久蜜臀av无| 免费在线观看黄色视频的| 最新的欧美精品一区二区| 亚洲欧美一区二区三区黑人| 满18在线观看网站| 嫁个100分男人电影在线观看 | 久久久久久久大尺度免费视频| 男女午夜视频在线观看| 精品亚洲成国产av| 80岁老熟妇乱子伦牲交| 国产又色又爽无遮挡免| 捣出白浆h1v1| 高清视频免费观看一区二区| 热99久久久久精品小说推荐| 国产午夜精品一二区理论片| 国产欧美日韩精品亚洲av| 国产又色又爽无遮挡免| 欧美老熟妇乱子伦牲交| 国产免费现黄频在线看| 母亲3免费完整高清在线观看| 成人黄色视频免费在线看| xxx大片免费视频| 多毛熟女@视频| 国产精品 欧美亚洲| xxx大片免费视频| 免费不卡黄色视频| 少妇的丰满在线观看| 亚洲免费av在线视频| 在线观看免费午夜福利视频| 各种免费的搞黄视频| 久久人人爽av亚洲精品天堂| 国产成人a∨麻豆精品| 咕卡用的链子| 欧美av亚洲av综合av国产av| 男女边摸边吃奶| √禁漫天堂资源中文www| 少妇猛男粗大的猛烈进出视频| 99国产精品99久久久久| 国产99久久九九免费精品| 亚洲av片天天在线观看| 国产精品久久久久成人av| 亚洲午夜精品一区,二区,三区| 美女视频免费永久观看网站| 国产精品久久久av美女十八| 国产一区二区三区综合在线观看| av网站免费在线观看视频| 亚洲国产中文字幕在线视频| 999精品在线视频| 黄色视频不卡| 一区二区av电影网| 国产精品99久久99久久久不卡| 久久毛片免费看一区二区三区| 国产精品国产三级专区第一集| 国产精品一区二区免费欧美 | 精品视频人人做人人爽| 在线天堂中文资源库| 久久国产精品大桥未久av| 男女免费视频国产| 老鸭窝网址在线观看| 又大又爽又粗| 少妇人妻 视频| 国产成人欧美在线观看 | 免费不卡黄色视频| 国产激情久久老熟女| 午夜久久久在线观看| 91老司机精品| 99久久综合免费| 婷婷色综合www| 亚洲五月色婷婷综合| 另类精品久久| 久久亚洲精品不卡| 国产精品av久久久久免费| 汤姆久久久久久久影院中文字幕| 首页视频小说图片口味搜索 | 日韩人妻精品一区2区三区| 国产高清videossex| 少妇 在线观看| 一二三四社区在线视频社区8| 日韩伦理黄色片| av电影中文网址| h视频一区二区三区| 亚洲精品日韩在线中文字幕| 亚洲国产欧美网| 国产在线观看jvid| 精品卡一卡二卡四卡免费| 天天躁夜夜躁狠狠躁躁| 日韩人妻精品一区2区三区| 人人澡人人妻人| 国产精品99久久99久久久不卡| 亚洲国产欧美日韩在线播放| tube8黄色片| 亚洲av美国av| 国产三级黄色录像| 性色av乱码一区二区三区2| 精品熟女少妇八av免费久了| 亚洲,一卡二卡三卡| 亚洲精品一卡2卡三卡4卡5卡 | 日韩 亚洲 欧美在线| 国产主播在线观看一区二区 | 国产精品久久久久成人av| 久久国产精品男人的天堂亚洲| 久久精品国产亚洲av涩爱| 免费高清在线观看视频在线观看| svipshipincom国产片| 亚洲成人手机| 欧美日韩av久久| 精品人妻一区二区三区麻豆| 午夜影院在线不卡| a级片在线免费高清观看视频| 人体艺术视频欧美日本| 天天躁夜夜躁狠狠久久av| 亚洲一卡2卡3卡4卡5卡精品中文| 黄色片一级片一级黄色片| 国产片内射在线| 又黄又粗又硬又大视频| 欧美+亚洲+日韩+国产| avwww免费| 黄色视频不卡| 电影成人av| 操美女的视频在线观看| 久久久久精品国产欧美久久久 | 欧美日韩一级在线毛片| 啦啦啦在线观看免费高清www| 91老司机精品| 超碰97精品在线观看| 国产欧美日韩一区二区三 | 国产成人系列免费观看| 成人亚洲欧美一区二区av| 日韩一卡2卡3卡4卡2021年| 日日爽夜夜爽网站| 中文乱码字字幕精品一区二区三区| 男人爽女人下面视频在线观看| 免费少妇av软件| 久久精品人人爽人人爽视色| 国产精品成人在线| 男女国产视频网站| 免费在线观看视频国产中文字幕亚洲 | 视频区欧美日本亚洲| 天天躁夜夜躁狠狠久久av| 少妇人妻久久综合中文| 欧美大码av| 午夜日韩欧美国产| 亚洲精品国产区一区二| 尾随美女入室| 亚洲av欧美aⅴ国产| 又大又黄又爽视频免费| 国产精品欧美亚洲77777| 亚洲成人手机| 久久精品久久久久久久性| 精品一区二区三卡| 免费高清在线观看日韩| 国产黄色免费在线视频| 国产精品 国内视频| 少妇人妻久久综合中文| 男人操女人黄网站| 青春草亚洲视频在线观看| www日本在线高清视频| 精品亚洲乱码少妇综合久久| 午夜福利一区二区在线看| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩视频精品一区| 美女扒开内裤让男人捅视频| 老司机深夜福利视频在线观看 | 手机成人av网站| 亚洲成人免费av在线播放| 久久久久精品人妻al黑| 最近最新中文字幕大全免费视频 | 999久久久国产精品视频| 人人妻人人爽人人添夜夜欢视频| 丝袜人妻中文字幕| 高潮久久久久久久久久久不卡| 九色亚洲精品在线播放| 国产又爽黄色视频| 久久人妻福利社区极品人妻图片 | 一级黄片播放器| 国产真人三级小视频在线观看| 亚洲精品美女久久av网站| 十八禁人妻一区二区| 黄片播放在线免费| 免费看av在线观看网站| 成年动漫av网址| 免费观看av网站的网址| 一级毛片 在线播放| 少妇猛男粗大的猛烈进出视频| 日韩中文字幕欧美一区二区 | 一级毛片 在线播放| 精品少妇内射三级| 韩国高清视频一区二区三区| 久久毛片免费看一区二区三区| 伊人亚洲综合成人网| 国产男人的电影天堂91| 亚洲成人手机| 婷婷成人精品国产| 欧美97在线视频| 午夜福利免费观看在线| 一区在线观看完整版| 久久99精品国语久久久| 一边亲一边摸免费视频| 亚洲国产欧美网| 9热在线视频观看99| 欧美变态另类bdsm刘玥| 视频在线观看一区二区三区| 国产成人一区二区三区免费视频网站 | 久久性视频一级片| 高清av免费在线| 99九九在线精品视频| 妹子高潮喷水视频| 亚洲,欧美,日韩| 亚洲天堂av无毛| av在线播放精品| 欧美亚洲日本最大视频资源| 亚洲中文av在线| 十分钟在线观看高清视频www| 超碰成人久久| 亚洲国产看品久久| 日韩熟女老妇一区二区性免费视频| 欧美人与性动交α欧美软件| 在线亚洲精品国产二区图片欧美| 日本黄色日本黄色录像| 又大又黄又爽视频免费| 我要看黄色一级片免费的| 亚洲少妇的诱惑av| 色94色欧美一区二区| 桃花免费在线播放| 亚洲国产欧美网| tube8黄色片| 亚洲国产精品一区二区三区在线| 欧美日韩亚洲国产一区二区在线观看 | 最近手机中文字幕大全| 亚洲精品美女久久av网站| 多毛熟女@视频| videos熟女内射| 欧美人与性动交α欧美精品济南到| 成年人午夜在线观看视频| 视频区欧美日本亚洲| 在线观看国产h片| 巨乳人妻的诱惑在线观看| 别揉我奶头~嗯~啊~动态视频 | 欧美变态另类bdsm刘玥| av又黄又爽大尺度在线免费看| 晚上一个人看的免费电影| 国产成人精品久久二区二区91| 男男h啪啪无遮挡| 亚洲第一青青草原| 国产亚洲午夜精品一区二区久久| 国产日韩欧美在线精品| 又大又黄又爽视频免费| 亚洲国产欧美日韩在线播放| 日本av免费视频播放| 狂野欧美激情性bbbbbb| 男女下面插进去视频免费观看| 精品卡一卡二卡四卡免费| 亚洲成av片中文字幕在线观看| 国产精品三级大全| 一区福利在线观看| 午夜激情久久久久久久| 嫩草影视91久久| 久久人妻熟女aⅴ| 桃花免费在线播放| 侵犯人妻中文字幕一二三四区| 国产高清不卡午夜福利| 精品高清国产在线一区| 国产欧美亚洲国产| 脱女人内裤的视频| 久久精品国产亚洲av涩爱| 伦理电影免费视频| 丝袜喷水一区| 亚洲五月婷婷丁香| 亚洲欧美清纯卡通| 人人妻人人爽人人添夜夜欢视频| 桃花免费在线播放| 久久人人爽人人片av| 亚洲成av片中文字幕在线观看| 亚洲人成电影免费在线| 国产精品久久久久成人av| 男女无遮挡免费网站观看| 欧美精品人与动牲交sv欧美| 日本av手机在线免费观看| 18禁观看日本| 少妇裸体淫交视频免费看高清 | 一本久久精品| 国产欧美日韩一区二区三 | 亚洲九九香蕉| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美激情在线| 亚洲成色77777| 爱豆传媒免费全集在线观看| 国产精品一二三区在线看| 黄色视频不卡| 七月丁香在线播放| 亚洲三区欧美一区| 最近中文字幕2019免费版| 中国国产av一级| 亚洲国产毛片av蜜桃av| 亚洲av欧美aⅴ国产| 少妇的丰满在线观看| 欧美xxⅹ黑人| 女人精品久久久久毛片| 最近手机中文字幕大全| 新久久久久国产一级毛片| 叶爱在线成人免费视频播放| 性少妇av在线| 咕卡用的链子| videos熟女内射| 老汉色av国产亚洲站长工具| 国产欧美日韩一区二区精品| 窝窝影院91人妻| 黄色 视频免费看| 亚洲第一青青草原| 别揉我奶头~嗯~啊~动态视频| 精品国内亚洲2022精品成人| 日韩欧美 国产精品| 国产精品爽爽va在线观看网站 | 成人三级做爰电影| 国产99白浆流出| 国内毛片毛片毛片毛片毛片| 国产又黄又爽又无遮挡在线| 亚洲av成人不卡在线观看播放网| 亚洲欧美一区二区三区黑人| 国产片内射在线| 草草在线视频免费看| 亚洲第一欧美日韩一区二区三区| 精品人妻1区二区| 国内精品久久久久久久电影| 午夜影院日韩av| e午夜精品久久久久久久| 国产1区2区3区精品| www.999成人在线观看| 欧美一级a爱片免费观看看 | 给我免费播放毛片高清在线观看| 亚洲精品中文字幕在线视频| 99精品久久久久人妻精品| 精品一区二区三区视频在线观看免费| 午夜精品在线福利| 国产精品综合久久久久久久免费| 国产亚洲av嫩草精品影院| 侵犯人妻中文字幕一二三四区| 国产精品久久视频播放| 精品少妇一区二区三区视频日本电影| 亚洲国产精品成人综合色| 精品一区二区三区av网在线观看| 国产又爽黄色视频| a级毛片在线看网站| www国产在线视频色| 一级黄色大片毛片| 国产欧美日韩一区二区三| 长腿黑丝高跟| 满18在线观看网站| 免费在线观看亚洲国产| 成年免费大片在线观看| 午夜免费鲁丝| 欧美日韩精品网址| 亚洲国产看品久久| av视频在线观看入口| 老司机午夜福利在线观看视频| 国产伦一二天堂av在线观看| 老司机靠b影院| 亚洲成av片中文字幕在线观看| 嫁个100分男人电影在线观看| 国产成人影院久久av| 91麻豆av在线| 搞女人的毛片| 此物有八面人人有两片| 深夜精品福利| 人成视频在线观看免费观看| 国产国语露脸激情在线看| 淫妇啪啪啪对白视频| 日韩欧美国产一区二区入口| 亚洲国产精品成人综合色| 亚洲国产中文字幕在线视频| 亚洲电影在线观看av| 美女高潮到喷水免费观看| 亚洲欧美激情综合另类| 少妇熟女aⅴ在线视频| 亚洲成av片中文字幕在线观看| 精品久久久久久成人av| 国产男靠女视频免费网站| 男人操女人黄网站| 亚洲专区国产一区二区| 国产精品自产拍在线观看55亚洲| 成人18禁在线播放| 免费电影在线观看免费观看| 亚洲精品av麻豆狂野| 嫩草影院精品99| 不卡av一区二区三区| 色播亚洲综合网| 亚洲最大成人中文| 国产精品自产拍在线观看55亚洲| 最近最新中文字幕大全电影3 | 久久精品夜夜夜夜夜久久蜜豆 | 精品乱码久久久久久99久播| 亚洲成a人片在线一区二区| 欧美性长视频在线观看| 国产精品二区激情视频| 欧美黄色淫秽网站| 成人特级黄色片久久久久久久| 精品免费久久久久久久清纯| 美国免费a级毛片| 欧美在线一区亚洲| 亚洲一区二区三区不卡视频| 最近最新免费中文字幕在线| 亚洲中文字幕日韩| 国产成人一区二区三区免费视频网站| 精华霜和精华液先用哪个| 校园春色视频在线观看| 国产成年人精品一区二区| 中文字幕最新亚洲高清| 亚洲精品久久成人aⅴ小说| 亚洲免费av在线视频| 美国免费a级毛片| 国产欧美日韩精品亚洲av| 国产成人啪精品午夜网站| 99国产极品粉嫩在线观看| 在线观看免费午夜福利视频| ponron亚洲| 国产真实乱freesex| 丝袜在线中文字幕| 免费看十八禁软件| 亚洲一区二区三区不卡视频| 亚洲国产欧洲综合997久久, | 少妇 在线观看| 午夜亚洲福利在线播放| 男男h啪啪无遮挡| 国产午夜精品久久久久久| 久久国产精品人妻蜜桃| 日本 av在线| 成人亚洲精品av一区二区| 久久国产乱子伦精品免费另类| 欧美成狂野欧美在线观看| 亚洲五月色婷婷综合| 日日干狠狠操夜夜爽| 两性午夜刺激爽爽歪歪视频在线观看 | 制服人妻中文乱码| 国产蜜桃级精品一区二区三区| 欧美一级毛片孕妇| 男人舔奶头视频| 麻豆国产av国片精品| 黄网站色视频无遮挡免费观看| 性色av乱码一区二区三区2| 国内精品久久久久精免费| 久久天堂一区二区三区四区| 中文字幕精品亚洲无线码一区 | 国产一区在线观看成人免费| 中亚洲国语对白在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 18禁黄网站禁片免费观看直播| 国产成人欧美在线观看| 欧美日本亚洲视频在线播放| 亚洲免费av在线视频| 精品久久久久久,| 久久久久久亚洲精品国产蜜桃av| 欧美+亚洲+日韩+国产| 中国美女看黄片| 亚洲国产欧美网| 免费电影在线观看免费观看| 性欧美人与动物交配| 久久久久亚洲av毛片大全| 男女视频在线观看网站免费 | 色婷婷久久久亚洲欧美| 亚洲精品色激情综合| 可以在线观看的亚洲视频| 妹子高潮喷水视频| 又大又爽又粗| 欧美日韩精品网址| 欧美激情高清一区二区三区| 亚洲一区二区三区不卡视频| 国产亚洲精品久久久久5区| 国产亚洲精品综合一区在线观看 | videosex国产| 窝窝影院91人妻| 19禁男女啪啪无遮挡网站| 怎么达到女性高潮| 91麻豆精品激情在线观看国产| 国产高清videossex| 日本 av在线| 午夜久久久在线观看| 91成人精品电影| 狂野欧美激情性xxxx| 国产在线观看jvid| 中文字幕人妻熟女乱码| 精品国产超薄肉色丝袜足j| 视频在线观看一区二区三区| 欧美绝顶高潮抽搐喷水| 亚洲熟妇熟女久久| 精品午夜福利视频在线观看一区| 此物有八面人人有两片| 操出白浆在线播放| 无限看片的www在线观看| 免费av毛片视频| 国产精品永久免费网站| 国产精品香港三级国产av潘金莲| 桃色一区二区三区在线观看| 亚洲真实伦在线观看| 久久久久久亚洲精品国产蜜桃av| 国产亚洲精品第一综合不卡| 老司机午夜十八禁免费视频| 村上凉子中文字幕在线| 欧美一级毛片孕妇| 国产高清激情床上av| 在线播放国产精品三级| 久久久久久久久久黄片| 黄色丝袜av网址大全| 波多野结衣高清作品| 中文字幕av电影在线播放| 亚洲性夜色夜夜综合| 亚洲成人免费电影在线观看| 久久久国产成人精品二区| 久久久久久久久久黄片| 亚洲最大成人中文| 亚洲成人久久爱视频| 日韩欧美一区视频在线观看| 亚洲av中文字字幕乱码综合 | 不卡一级毛片| 欧美黑人精品巨大| 丁香六月欧美| 大型av网站在线播放| 久久久久久人人人人人| 成人18禁在线播放| 俄罗斯特黄特色一大片| 身体一侧抽搐| 亚洲国产中文字幕在线视频| 久久久久久久午夜电影| 又紧又爽又黄一区二区| 十分钟在线观看高清视频www| 国产精品二区激情视频| 亚洲av五月六月丁香网| 热99re8久久精品国产| 欧美在线黄色| 亚洲欧洲精品一区二区精品久久久| 国产精品二区激情视频| 精品国内亚洲2022精品成人| 黄色视频不卡| 91在线观看av| 母亲3免费完整高清在线观看| 亚洲三区欧美一区| 欧美人与性动交α欧美精品济南到| 亚洲精品在线美女| 两性夫妻黄色片| 一区二区三区激情视频| 亚洲精华国产精华精| 少妇 在线观看| 青草久久国产| 每晚都被弄得嗷嗷叫到高潮| 99国产精品99久久久久| 日本免费a在线| 亚洲一区中文字幕在线| 国产精品精品国产色婷婷| 免费一级毛片在线播放高清视频| 久久久国产精品麻豆| 成人永久免费在线观看视频| 免费看美女性在线毛片视频| 看片在线看免费视频| 天堂√8在线中文| 亚洲在线自拍视频| 免费高清在线观看日韩| 色老头精品视频在线观看|