• <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在线播放网站| 三级男女做爰猛烈吃奶摸视频| 一进一出抽搐gif免费好疼| 国产一区二区在线观看日韩| 国产精品女同一区二区软件| 日韩亚洲欧美综合| 国产在线精品亚洲第一网站| 亚洲电影在线观看av| 亚洲国产精品sss在线观看| 色哟哟·www| 一本久久中文字幕| 免费看日本二区| 欧美最黄视频在线播放免费| 亚洲中文字幕一区二区三区有码在线看| www日本黄色视频网| 麻豆国产97在线/欧美| 最近手机中文字幕大全| 校园春色视频在线观看| 亚洲性久久影院| 老熟妇仑乱视频hdxx| 免费电影在线观看免费观看| 免费观看的影片在线观看| 在线播放国产精品三级| 欧美日韩一区二区视频在线观看视频在线 | 中文字幕人妻熟人妻熟丝袜美| 欧美人与善性xxx| 22中文网久久字幕| 久久久久久国产a免费观看| 国产精华一区二区三区| 亚洲在线自拍视频| 欧美三级亚洲精品| 亚洲18禁久久av| 六月丁香七月| 在线播放国产精品三级| 成人特级黄色片久久久久久久| 国产精品av视频在线免费观看| 观看美女的网站| 久久久久性生活片| 国产伦一二天堂av在线观看| av中文乱码字幕在线| 一级毛片电影观看 | 波多野结衣高清作品| 久久韩国三级中文字幕| 国产精品永久免费网站| 在线观看免费视频日本深夜| av免费在线看不卡| 亚洲婷婷狠狠爱综合网| 久久久久久久久久久丰满| 特级一级黄色大片| 97人妻精品一区二区三区麻豆| 丰满乱子伦码专区| 久久人人爽人人片av| 国产激情偷乱视频一区二区| 亚洲激情五月婷婷啪啪| 久久久精品94久久精品| 精品无人区乱码1区二区| 免费观看精品视频网站| 欧美三级亚洲精品| 97在线视频观看| 久久九九热精品免费| 日韩欧美国产在线观看| av福利片在线观看| 99热网站在线观看| 一本一本综合久久| 午夜视频国产福利| ponron亚洲| 天堂√8在线中文| 在现免费观看毛片| av专区在线播放| av.在线天堂| 久久婷婷人人爽人人干人人爱| 91在线精品国自产拍蜜月| aaaaa片日本免费| 免费观看精品视频网站| 国内少妇人妻偷人精品xxx网站| 亚洲自偷自拍三级| 老司机影院成人| 亚洲在线自拍视频| 天天躁夜夜躁狠狠久久av| 欧美精品国产亚洲| 欧美xxxx性猛交bbbb| 午夜日韩欧美国产| 日本免费在线观看一区| 一级毛片 在线播放| videossex国产| 乱人伦中国视频| 国产精品一区二区三区四区免费观看| 少妇人妻 视频| 国产高清有码在线观看视频| 狠狠精品人妻久久久久久综合| 国产精品.久久久| 亚洲欧洲国产日韩| 国产亚洲精品久久久com| 亚洲av在线观看美女高潮| 老司机影院毛片| 亚洲国产av新网站| 亚洲四区av| 久久久久久久亚洲中文字幕| 在线观看av片永久免费下载| 亚洲av电影在线观看一区二区三区| 在线免费观看不下载黄p国产| 人妻 亚洲 视频| 美女大奶头黄色视频| 大片电影免费在线观看免费| av女优亚洲男人天堂| 免费少妇av软件| 亚洲欧美精品自产自拍| 日韩av免费高清视频| 精品一区二区三区视频在线| 国产精品人妻久久久久久| 亚洲精品乱码久久久久久按摩| av又黄又爽大尺度在线免费看| 国语对白做爰xxxⅹ性视频网站| 久久久久久久久久久丰满| .国产精品久久| 日韩不卡一区二区三区视频在线| 亚洲va在线va天堂va国产| 精品亚洲成a人片在线观看| 国产av精品麻豆| 亚洲综合色惰| 亚洲国产色片| 亚洲精品一二三| 久久久午夜欧美精品| 我要看日韩黄色一级片| 亚洲国产欧美日韩在线播放 | 最近手机中文字幕大全| 国产日韩欧美在线精品| 亚洲真实伦在线观看| 亚洲国产精品一区三区| 亚洲人成网站在线播| 欧美日韩综合久久久久久| 久久久久视频综合| 狂野欧美白嫩少妇大欣赏| 国产在线男女| 一区在线观看完整版| 久久 成人 亚洲| 精品99又大又爽又粗少妇毛片| 夫妻性生交免费视频一级片| 中文精品一卡2卡3卡4更新| 简卡轻食公司| 免费看不卡的av| 国产精品成人在线| 欧美日韩国产mv在线观看视频| 丝袜喷水一区| 精品人妻偷拍中文字幕| 91在线精品国自产拍蜜月| 99久久人妻综合| 熟女电影av网| 一级二级三级毛片免费看| 特大巨黑吊av在线直播| 国产精品一区二区在线观看99| 七月丁香在线播放| 午夜激情福利司机影院| 亚洲中文av在线| 少妇人妻 视频| 免费观看的影片在线观看| av专区在线播放| 亚洲美女黄色视频免费看| 91精品伊人久久大香线蕉| 久久精品国产鲁丝片午夜精品| 国产精品久久久久久久电影| 色哟哟·www| 亚洲av综合色区一区| 精品酒店卫生间| 成人国产av品久久久| 色视频www国产| 黄色配什么色好看| 少妇人妻久久综合中文| 亚洲国产精品成人久久小说| 久久av网站| 日韩伦理黄色片| 插阴视频在线观看视频| 国产成人午夜福利电影在线观看| 极品人妻少妇av视频| 国产亚洲精品久久久com| 国产欧美日韩精品一区二区| av免费观看日本| 亚洲成人手机| 成人综合一区亚洲| 日本色播在线视频| 日韩亚洲欧美综合| 亚洲中文av在线| 日韩三级伦理在线观看| 婷婷色av中文字幕| 天堂中文最新版在线下载| 97在线视频观看| 如何舔出高潮| 纵有疾风起免费观看全集完整版| av在线播放精品| 国产亚洲av片在线观看秒播厂| 亚洲精品aⅴ在线观看| 欧美性感艳星| 国产成人精品一,二区| 日日啪夜夜撸| 亚洲va在线va天堂va国产| 国产av一区二区精品久久| 免费观看a级毛片全部| 能在线免费看毛片的网站| 久久久久久久亚洲中文字幕| 精品久久久久久久久亚洲| 一本久久精品| 亚洲欧美日韩东京热| 女人精品久久久久毛片| 国产无遮挡羞羞视频在线观看| 久久女婷五月综合色啪小说| 又粗又硬又长又爽又黄的视频| 欧美xxxx性猛交bbbb| 国产欧美亚洲国产| 国产国拍精品亚洲av在线观看| 国产中年淑女户外野战色| 一个人看视频在线观看www免费| 亚洲综合色惰| 另类精品久久| 大话2 男鬼变身卡| 嫩草影院新地址| 日本91视频免费播放| 精品人妻偷拍中文字幕| 热re99久久国产66热| 久久毛片免费看一区二区三区| 亚洲av成人精品一二三区| 国产片特级美女逼逼视频| 欧美人与善性xxx| 伦理电影免费视频| 免费看日本二区| 免费观看在线日韩| av一本久久久久| 久久久久久人妻| 纵有疾风起免费观看全集完整版| av天堂久久9| 高清在线视频一区二区三区| 国产精品偷伦视频观看了| 插阴视频在线观看视频| 亚洲在久久综合| 桃花免费在线播放| 免费看av在线观看网站| 边亲边吃奶的免费视频| 色视频在线一区二区三区| 国产色爽女视频免费观看| 大话2 男鬼变身卡| 精品少妇久久久久久888优播| 中文字幕人妻熟人妻熟丝袜美| 国产精品国产三级专区第一集| 性色avwww在线观看| 久久久欧美国产精品| 亚洲国产精品专区欧美| 三级国产精品片| 夫妻性生交免费视频一级片| 伦理电影免费视频| tube8黄色片| 国产精品无大码| 国产精品久久久久成人av| 久久av网站| 我的女老师完整版在线观看| 亚洲精品中文字幕在线视频 | 天美传媒精品一区二区| 秋霞伦理黄片| 爱豆传媒免费全集在线观看| 老熟女久久久| 哪个播放器可以免费观看大片| 一本一本综合久久| .国产精品久久| 亚洲精品国产成人久久av| 久久久久久久久久久丰满| 日韩av在线免费看完整版不卡| 在线观看免费高清a一片| 久久婷婷青草| 少妇丰满av| 精品一区二区三区视频在线| 黑丝袜美女国产一区| 18禁在线播放成人免费| 黄色欧美视频在线观看| 视频区图区小说| 狂野欧美激情性bbbbbb| 高清黄色对白视频在线免费看 | 亚洲欧美一区二区三区黑人 | 男女免费视频国产| 成人影院久久| 青春草国产在线视频| 看免费成人av毛片| av又黄又爽大尺度在线免费看| 99热这里只有是精品在线观看| 亚洲成色77777| 99久久中文字幕三级久久日本| 亚洲欧美一区二区三区国产| 欧美三级亚洲精品| 欧美日韩国产mv在线观看视频| 欧美+日韩+精品| 又粗又硬又长又爽又黄的视频| 自线自在国产av| 亚洲一区二区三区欧美精品| 国语对白做爰xxxⅹ性视频网站| 国产毛片在线视频| √禁漫天堂资源中文www| 九九爱精品视频在线观看| 只有这里有精品99| 精品人妻熟女av久视频| 寂寞人妻少妇视频99o| 免费观看无遮挡的男女| 精品少妇黑人巨大在线播放| 91久久精品国产一区二区三区| 国产精品麻豆人妻色哟哟久久| 少妇丰满av| 国产亚洲一区二区精品| 男人舔奶头视频| 春色校园在线视频观看| 国产精品一区二区在线观看99| 亚洲情色 制服丝袜| 又黄又爽又刺激的免费视频.| 国产一区二区三区av在线| 国产老妇伦熟女老妇高清| 欧美激情极品国产一区二区三区 | 国产一区亚洲一区在线观看| 亚洲av在线观看美女高潮| 亚洲美女黄色视频免费看| 国产美女午夜福利| 美女大奶头黄色视频| 国产精品久久久久成人av| 岛国毛片在线播放| 在线观看国产h片| 丝袜喷水一区| 日日撸夜夜添| 三上悠亚av全集在线观看 | 丝瓜视频免费看黄片| freevideosex欧美| 最近2019中文字幕mv第一页| 夜夜看夜夜爽夜夜摸| 国产免费视频播放在线视频| 精品少妇黑人巨大在线播放| 观看美女的网站| 日韩在线高清观看一区二区三区| 日本猛色少妇xxxxx猛交久久| 一个人免费看片子| 国产伦精品一区二区三区视频9| 久热久热在线精品观看| 亚洲真实伦在线观看| 久久久久精品性色| 天堂中文最新版在线下载| 久久韩国三级中文字幕| 欧美激情极品国产一区二区三区 | 又粗又硬又长又爽又黄的视频| 亚洲图色成人| 久久国产乱子免费精品| 日本黄色片子视频| 国产午夜精品久久久久久一区二区三区| 国产精品久久久久久精品电影小说| 欧美区成人在线视频| 免费观看无遮挡的男女| 亚洲久久久国产精品| 亚洲精品国产色婷婷电影| 一级毛片aaaaaa免费看小| 久久久久久久久大av| 亚洲不卡免费看| 一区在线观看完整版| 久久久久国产精品人妻一区二区| 久久热精品热| 女性被躁到高潮视频| 人体艺术视频欧美日本| 精品卡一卡二卡四卡免费| 人妻制服诱惑在线中文字幕| 免费观看无遮挡的男女| 国产白丝娇喘喷水9色精品| av在线播放精品| 国产精品99久久久久久久久| 久久人人爽av亚洲精品天堂| 亚洲精华国产精华液的使用体验| 亚洲欧洲日产国产| 久久婷婷青草| 午夜久久久在线观看| 国产一区有黄有色的免费视频| 亚洲国产精品专区欧美| 十八禁高潮呻吟视频 | 久久久a久久爽久久v久久| 国产极品粉嫩免费观看在线 | 国产一级毛片在线| 亚洲熟女精品中文字幕| 蜜桃久久精品国产亚洲av| 亚洲美女视频黄频| 十八禁高潮呻吟视频 | 亚洲精华国产精华液的使用体验| 免费大片18禁| 亚洲中文av在线| 热re99久久国产66热| 波野结衣二区三区在线| 性色avwww在线观看| 在线天堂最新版资源| 免费观看的影片在线观看| 日本黄大片高清| 一本一本综合久久| 新久久久久国产一级毛片| 国产一区二区在线观看av| 99久久人妻综合| 国产成人免费观看mmmm| 在线观看一区二区三区激情| 亚洲欧洲日产国产| 色5月婷婷丁香| 亚洲情色 制服丝袜| 欧美人与善性xxx| 久久精品久久精品一区二区三区| 国产精品一区二区在线观看99| 欧美日韩国产mv在线观看视频| 国产av码专区亚洲av| 春色校园在线视频观看| 久久久久久久久久久久大奶| 97精品久久久久久久久久精品| 亚洲精品视频女| 一级黄片播放器| 久久99蜜桃精品久久| 免费看av在线观看网站| 亚洲美女视频黄频| 国产日韩欧美在线精品| 精品一区在线观看国产| 国产成人freesex在线| 卡戴珊不雅视频在线播放| 看十八女毛片水多多多| 我要看黄色一级片免费的| 天堂中文最新版在线下载| 国产精品久久久久久精品电影小说| 男女啪啪激烈高潮av片| 麻豆成人午夜福利视频| 亚洲精品,欧美精品| 国产精品蜜桃在线观看| 日韩精品有码人妻一区| 观看免费一级毛片| 国国产精品蜜臀av免费| 91精品一卡2卡3卡4卡| 久久久久久久精品精品| 黑丝袜美女国产一区| 99久久精品国产国产毛片| 天天操日日干夜夜撸| 久久精品国产亚洲网站| 欧美亚洲 丝袜 人妻 在线| 妹子高潮喷水视频| 亚洲伊人久久精品综合| 久久精品熟女亚洲av麻豆精品| 亚洲欧洲精品一区二区精品久久久 | av在线app专区| 97超视频在线观看视频| 欧美变态另类bdsm刘玥| 免费观看性生交大片5| 久久久久国产网址| 一级毛片我不卡| 边亲边吃奶的免费视频| 另类精品久久| 日本vs欧美在线观看视频 | 综合色丁香网| 国产在线免费精品| av不卡在线播放| 三级国产精品欧美在线观看| 亚洲精品成人av观看孕妇| 欧美三级亚洲精品| 99久久精品国产国产毛片| 黄色配什么色好看| 亚洲va在线va天堂va国产| 99精国产麻豆久久婷婷| 免费观看的影片在线观看| 大又大粗又爽又黄少妇毛片口| 中国三级夫妇交换| av专区在线播放| 亚洲av电影在线观看一区二区三区| 免费人成在线观看视频色| 日韩人妻高清精品专区| 色婷婷久久久亚洲欧美| av在线观看视频网站免费| 精品一区二区三区视频在线| 国产亚洲5aaaaa淫片| 中文乱码字字幕精品一区二区三区| 日日啪夜夜撸| 99热这里只有是精品在线观看| 日韩伦理黄色片| 国产一级毛片在线| 女人精品久久久久毛片| 老熟女久久久| 亚洲欧洲日产国产| 国产欧美日韩精品一区二区| 欧美日韩视频高清一区二区三区二| 午夜免费男女啪啪视频观看| 国精品久久久久久国模美| 日日摸夜夜添夜夜添av毛片| 汤姆久久久久久久影院中文字幕| 18禁在线播放成人免费| 伊人久久精品亚洲午夜| 91在线精品国自产拍蜜月| 全区人妻精品视频| 精品国产国语对白av| 中文天堂在线官网| 亚洲欧美中文字幕日韩二区| 久久99精品国语久久久| 国产视频内射| www.色视频.com| 中文字幕制服av| 五月开心婷婷网| 大片免费播放器 马上看| 寂寞人妻少妇视频99o| 日韩欧美精品免费久久| 久久热精品热| 欧美xxxx性猛交bbbb| 少妇被粗大猛烈的视频| 亚洲精品国产av蜜桃| 亚洲av福利一区| 国产日韩欧美在线精品| 一级毛片久久久久久久久女| 日韩成人伦理影院| 国产免费一区二区三区四区乱码| 黄色日韩在线| 一级片'在线观看视频| 精品人妻熟女毛片av久久网站| 久久ye,这里只有精品| 亚洲经典国产精华液单| 六月丁香七月| 国产高清不卡午夜福利| 精品久久久精品久久久| 91午夜精品亚洲一区二区三区| 亚洲国产欧美日韩在线播放 | 国产成人aa在线观看| 日韩伦理黄色片| 亚洲人与动物交配视频| 精品国产一区二区久久| 另类亚洲欧美激情| 久久ye,这里只有精品| 看免费成人av毛片| 久久国产亚洲av麻豆专区| 人妻 亚洲 视频| 久久久久久伊人网av| 最新的欧美精品一区二区| 午夜免费鲁丝| 乱人伦中国视频| 午夜av观看不卡| 狂野欧美激情性bbbbbb| 国产精品偷伦视频观看了| 国产精品久久久久久av不卡| 色5月婷婷丁香| 精品一区二区免费观看| 美女视频免费永久观看网站| 亚洲欧美精品自产自拍| 亚洲国产色片| 成年美女黄网站色视频大全免费 | 人人妻人人澡人人爽人人夜夜| 亚洲一级一片aⅴ在线观看| 18禁动态无遮挡网站| 人妻制服诱惑在线中文字幕| 搡女人真爽免费视频火全软件| 亚洲天堂av无毛| 欧美日韩视频高清一区二区三区二| 亚洲欧美清纯卡通| 国产精品一区二区三区四区免费观看| 亚洲成人手机| 国产男女超爽视频在线观看| 亚洲精品国产av成人精品| 日本免费在线观看一区| 18禁在线播放成人免费| 青春草视频在线免费观看| 亚洲真实伦在线观看| 国产精品久久久久久精品电影小说| 午夜福利在线观看免费完整高清在| 成人漫画全彩无遮挡| 精品亚洲乱码少妇综合久久| 欧美亚洲 丝袜 人妻 在线| 国产乱来视频区| 人体艺术视频欧美日本| 在线观看美女被高潮喷水网站| av卡一久久| 午夜免费男女啪啪视频观看| 极品人妻少妇av视频| 人妻系列 视频| freevideosex欧美| 精品久久久精品久久久| 国产在线免费精品| 少妇猛男粗大的猛烈进出视频| 亚洲三级黄色毛片| 99久久精品热视频| 中文在线观看免费www的网站| 日本与韩国留学比较| 精品久久久精品久久久| 99国产精品免费福利视频| 少妇被粗大猛烈的视频| 嘟嘟电影网在线观看| 成人亚洲欧美一区二区av| 亚洲精品久久午夜乱码| 精品酒店卫生间| 观看美女的网站| 不卡视频在线观看欧美| 大香蕉97超碰在线| 免费观看在线日韩| 欧美成人精品欧美一级黄| 一区在线观看完整版| 伊人久久国产一区二区| 观看美女的网站| 国产成人精品无人区| 乱系列少妇在线播放| 日日撸夜夜添| 日韩,欧美,国产一区二区三区| 日韩 亚洲 欧美在线| 色94色欧美一区二区| 久久99热6这里只有精品| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲自偷自拍三级| 三级经典国产精品| 黄色毛片三级朝国网站 | 极品教师在线视频| 成人国产麻豆网| 久久av网站| 下体分泌物呈黄色| 亚洲精品日韩在线中文字幕| 又黄又爽又刺激的免费视频.| 国产成人免费无遮挡视频| 国产av码专区亚洲av| 亚洲国产成人一精品久久久| www.色视频.com| 亚洲精品一二三| 中文乱码字字幕精品一区二区三区| 美女大奶头黄色视频| 国产亚洲av片在线观看秒播厂| 一区二区三区四区激情视频| 三级经典国产精品| 婷婷色av中文字幕|