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

    縫翼流場結構與低頻窄帶噪聲形成特性研究

    2021-07-12 12:26:00王紅建張巧
    西北工業(yè)大學學報 2021年3期
    關鍵詞:尾緣噪聲源窄帶

    王紅建, 張巧

    (西北工業(yè)大學 航空學院, 陜西 西安 710072)

    飛機起降過程中,前緣縫翼是機體噪聲的重要噪聲源之一[1]??p翼噪聲頻譜主要分為:寬帶、低頻窄帶和高頻窄帶等。普遍認為,高頻窄帶噪聲是縫翼尾緣渦脫落引起的;寬帶噪聲則與剪切層周期性撞擊壁面有關;Khorrami等[2]認為加速通過縫道區(qū)域的對流渦是產生低頻窄帶噪聲的主要原因;Hein等[3]認為低頻窄帶噪聲與聲共振有關;而Imamura等[4]將低頻窄帶噪聲與剪切層的非定常振蕩和縫道區(qū)域的不穩(wěn)定性聯(lián)系起來;由此可知,各學者對低頻窄帶噪聲的產生機理仍有爭議。

    為明確縫翼低頻窄帶噪聲的產生機理,國內外學者已進行過相關研究。K?nig等[5]結合大渦模擬(LES)和聲擾動方程(APE)求解遠場氣動噪聲,發(fā)現(xiàn)低頻窄帶噪聲(1~3 kHz頻率范圍內)集中在縫道區(qū)域;且其噪聲由縫翼尾緣與主翼之間的壓力脈動產生;Terracol等[6]采用RANS和LES方法模擬兩段翼中縫翼凹腔區(qū)域的非定常流動和噪聲輻射,發(fā)現(xiàn)凹腔中的反饋回路是產生離散低頻噪聲的原因;并在Rossiter等[7]研究基礎上完善了低頻窄帶噪聲的預測模型,但其渦結構在剪切層上勻速運動的假設較為理想化。

    POD和DMD方法的提出對縫翼低頻噪聲產生機理的研究具有重要意義。Souza等[8-9]使用POD方法分析低頻窄帶噪聲對應頻率的主要流動特征,發(fā)現(xiàn)窄帶噪聲的幅值由剪切層再附著區(qū)及縫道加速流控制;而縫翼凹腔區(qū)域的流動主要與剪切層非定常脈動有關。進一步采用湍動能和壓力脈動范數(shù)進行POD分析,發(fā)現(xiàn)大尺度渦結構與噪聲頻譜主要峰值頻率相聯(lián)系。渦結構通常與剪切層非定常振蕩有關,其壓力脈動與遠場噪聲的相關性很高。而湍動能與遠場噪聲的相關性不大,但在使用壓力內積求解POD模態(tài)時,無法消除模擬產生的背景壓力脈動,這會顯著影響POD模態(tài),使其對應的流場特性存在一定差異。隨后,陳鵬等[10]采用DDES方法獲得兩段翼流場特性,利用POD方法研究縫翼周圍的壓力脈動與噪聲輻射的關系,發(fā)現(xiàn)再附著區(qū)與主翼前緣是遠場噪聲輻射的關鍵區(qū)域,且主模態(tài)由加速流穿過縫道區(qū)域的對流渦引起,表明主要噪聲源位于縫道區(qū)域;此外,魏佳云等[11]采用LES方法研究了兩段翼中的縫翼噪聲,并通過DMD發(fā)現(xiàn)低頻噪聲源于剪切層中大尺度渦結構與縫翼尾緣下壁面的周期性撞擊效應。至今,各學者對低頻窄帶噪聲機理的研究仍未達成一致,因此,有必要進一步研究縫翼低頻噪聲的產生與傳播機理。

    本文基于30P-30N翼型,采用DDES方法(瞬態(tài)延時脫體渦模擬)和FW-H聲類比方程分析獲得縫翼周圍流場和噪聲的基本分布特性;然后利用POD和DMD方法,抽取瞬態(tài)流場結構特征,并結合傅里葉變換(FFT),研究流場結構的頻譜特性;最后利用EMD和互相關分析技術,建立低頻窄帶噪聲與流場模態(tài)結構之間的內在聯(lián)系,以揭示其產生機理。

    1 仿真模型與結果分析

    1.1 網(wǎng)格劃分

    本文采用30P-30N高升力幾何構型進行研究??p翼結構參數(shù)見表1,仿真模擬的自由來流條件與試驗條件一致,如表2所示。

    表1 縫翼結構參數(shù)

    表2 流場模擬工況

    1.2 網(wǎng)格劃分及計算方法

    多段翼網(wǎng)格(圖1)節(jié)點總數(shù)約1 000萬,第一層網(wǎng)格到壁面距離滿足黏性無量綱尺度與DDES方法對邊界層網(wǎng)格的要求[12]。遠場邊界位于翼型收緊前緣點50倍弦長處,其單位展向長度為Δz=1.98×10-3C,計算域展向尺度滿足縫翼流場的非定常模擬計算要求[13]。

    圖1 模型CFD計算網(wǎng)格

    1.3 仿真模型驗證與結果分析

    基于上述模型仿真結果,對比JAXA[14]實驗結果(見圖2),可知二者的翼型表面壓力系數(shù)吻合較好,初步驗證本文流場仿真方法的準確性。

    圖2 壓力系數(shù)對比

    瞬態(tài)展向渦量見圖3。對比文獻[15]可知,本文仿真結果與粒子圖像測速法(PIV)實驗結果[15]基本一致,較準確地捕捉到剪切層的發(fā)展過程,進一步驗證本文瞬態(tài)流場分析的正確性。

    圖3 瞬態(tài)展向渦ωzc/U∞

    相鄰渦結構及其與壁面之間的相互作用是縫翼噪聲產生的主要原因。本文采用Q準則云圖進行渦結構的可視化分析。圖4顯示了縫翼渦結構的整體發(fā)展過程。由圖可知,隨著自由剪切層的向后發(fā)展,渦結構發(fā)生了較大改變,其發(fā)展規(guī)律與劉超群等[16]的關于分離氣流中渦結構的形成與發(fā)展規(guī)律類似,即剪切層上渦結構從Λ渦演變?yōu)榘l(fā)卡渦,且卡渦含有較大能量,這是縫翼剪切層發(fā)展的重要特征,并會對噪聲源的形成起到關鍵作用。

    圖4 Q準則瞬態(tài)展向渦量圖

    回流區(qū)含有大量低能量的離散渦,這些渦與剪切層相互耦合,一方面對再循環(huán)區(qū)域中渦結構產生影響;另一方面也增強了剪切層的不穩(wěn)定性,使得剪切層中下游區(qū)域中的渦呈現(xiàn)大量混亂的結構形式,因此該區(qū)域是瞬態(tài)流場行為特征的主要集中區(qū)。剪切層與尾緣壁面撞擊區(qū)的流場結構中,含有大量撞擊形成的離散渦,其中一些渦結構在縫道加速氣流的作用下,發(fā)生了強烈變形,這種渦結構短時內的形狀變化會導致局部流場結構能量的劇變,也是縫翼低頻窄帶噪聲形成的重要方式。

    圖5為縫翼典型位置監(jiān)測點的壓力脈動分布圖。觀察發(fā)現(xiàn)中低頻范圍內,縫翼尾緣處點(P4)和主翼前緣點(P5)壓力脈動水平遠高于其他監(jiān)測點,這也反映出縫翼尾緣不穩(wěn)定發(fā)卡渦的能量脈動及其與主翼前緣的相互作用;圖3顯示P4、P5展向渦量較小,說明這兩處的壓力脈動并非來自展向渦,進一步說明發(fā)卡渦的能量脈動相對較大,也是窄帶噪聲輻射的重要原因之一;而P2和P3點主要在凹腔回流區(qū),其壓力脈動則相對較低。

    圖5 縫翼凹腔內部近場監(jiān)測點壓力脈動頻譜

    分析遠場噪聲時,選擇以收緊翼型前緣點為圓心,翼型10倍巡航弦長為半徑,于290°方向設置觀測點。采用Welch方法分析頻譜,遠場噪聲頻譜(見圖6)中出現(xiàn)3個主要的窄帶峰值(1 343,1 984,2 653) Hz,這與BANC Ⅲ[17]中所預測的結果基本一致,進一步驗證了本文噪聲預測方法的準確性。

    圖6 遠場噪聲頻譜

    2 低頻窄帶噪聲機理研究

    POD和DMD[18]方法可抽取瞬態(tài)流場的主要特征,對探究大尺度渦與低頻窄帶噪聲輻射的內在機理具有重要意義。研究表明[10],前緣縫翼周圍的壓力脈動水平與其遠場噪聲輻射密切相關。POD和DMD方法可將瞬態(tài)流場的壓力脈動序列分解為流場特征模態(tài)。選取主要流動模態(tài)進行分析,研究其特征模態(tài)的能量構成及其頻譜特性。

    2.1 基于POD方法分析窄帶噪聲產生機理

    2.1.1 POD方法

    POD方法可識別流場中含能最大的模態(tài),對高階動力學系統(tǒng)進行降維處理。具體操作方法為:提取各時刻流場的壓力值,整個壓力場p(x,t)分解為基本流動和壓力脈動量p′(x,t)疊加。將離散時刻流場提取的N個快照,排列成快照矩陣A

    A=[p′(x,t1),…,p′(x,tN)]

    (1)

    由于該模型中,切片上的網(wǎng)格節(jié)點總數(shù)(m=304 989)遠大于快照個數(shù)(n=1 000),因此本文采用SPOD[19]方法進行后續(xù)求解。首先求解中間矩陣C=ATA;可根據(jù)中間矩陣的特征值(λ)和特征向量(V),得POD的模態(tài)(Φ)和模態(tài)系數(shù)(a)為

    Φ=AV

    (2)

    a=ATΦ

    (3)

    根據(jù)公式(4)可得各階模態(tài)占總模態(tài)能量的累計百分比。

    (4)

    式中:n代表特征值數(shù),EM表示前M階模態(tài)的能量占比之和。

    2.1.2 POD分解結果分析

    圖7為POD前50階模態(tài)特征值,其中第一階模態(tài)約占總模態(tài)能量的98%,第二階模態(tài)約2%,而第四階之后其模態(tài)能量低于0.1%。選取總模態(tài)能量占比99.9%的前四階模態(tài)進行分析。

    圖7 POD特征值

    圖8為具有代表性第一、三階模態(tài)的壓力云圖。圖中看到,第一階模態(tài)捕捉到壓力動態(tài)主要在剪切層下游及碰撞區(qū)。顯示出這些區(qū)域是該模態(tài)能量的集中區(qū),也是主要噪聲源。相比于第一階模態(tài),隨后各階模態(tài)的正負壓力脈動區(qū)域明顯增大,而其脈動間隔尺度有較大程度地減小;模態(tài)壓力分布的主要區(qū)域在靠近尾緣的剪切層沿線上。

    圖8 一、三階POD模態(tài)

    與此同時,第三階模態(tài)云圖中縫翼凹腔和剪切層路徑上也存在不同強度的壓力脈動區(qū),但其模態(tài)能量較弱,所對應的噪聲源強度也較弱。因POD模態(tài)(具有寬頻特性)不能反映窄帶頻率特征,故下面利用DMD方法對流場特性進一步分析。

    2.2 基于DMD方法分析低頻窄帶噪聲機理

    2.2.1 DMD方法

    DMD方法可求解特定頻率下,模態(tài)對應的流場結構特征。具體方法為:將上述流動快照分為兩組X和Y:

    X=[p′(x,t1),…,p′(x,tN-1)]

    (5)

    Y=[p′(x,t2),…,p′(x,tN)]

    (6)

    則滿足Y≈A*X,其中X可分解為

    X≈UΣV*

    (7)

    式中,U,V均為酉矩陣,即U*U=I,V*V=I。

    A=YVΣ-1U*

    (8)

    (9)

    Ψ=YVΣ-1W=UW

    (10)

    本文采用Ψ=UW表示DMD的模態(tài)。在確定特征值后,可根據(jù)DMD模態(tài)求解特征頻率和模態(tài)增長率

    μi=lg(Wi)/Δt

    (11)

    μi實部表示DMD模態(tài)增長率,虛部表示特征頻率。

    2.2.2 DMD分解結果分析

    DMD模態(tài)特征值分布見圖9??芍獔D中特征值絕大部分在單位圓上,部分在單位圓內,該分布對應于系統(tǒng)臨界穩(wěn)定狀態(tài)。

    圖9 DMD模態(tài)特征值

    按照模態(tài)振幅αi=W-1UHp′(x,t1)對不同頻率模態(tài)排序。第一階模態(tài)為靜態(tài)模態(tài),之后模態(tài)成對出現(xiàn)(共軛模態(tài))。選取DMD前九階模態(tài)進行分析,并對2~3、6~7階模態(tài)云圖展開介紹,見圖10,其中低模態(tài)對應于低頻。DMD第二~三階模態(tài)的壓力脈動主要集中在尾緣渦碰撞區(qū);而第六~七階模態(tài)壓力脈動區(qū)域沿剪切層擴散。從而證明低頻窄帶噪聲源主要在縫翼尾緣附近,而較高頻率聲源則更多地與縫翼凹腔和剪切層渦結構相聯(lián)系。

    圖10 DMD模態(tài)(2對共軛模態(tài))

    圖8a)POD第一階模態(tài)(具有寬頻特性)也表現(xiàn)出尾緣區(qū)域壓力脈動的分布特征,而該模態(tài)的寬頻特性似乎說明尾緣渦碰撞區(qū)并非只包含低頻噪聲,該區(qū)域在一定程度上也包含高頻噪聲分量,但其強度相對較弱。

    POD與DMD壓力模態(tài)云圖是流場壓力瞬態(tài)變化特征的解構,在一定程度上反映了縫翼聲源的分布。然而,必須注意到,流場壓力脈動與聲壓傳播是有一定差別的;其原因就是有相當大的壓力脈動能量并沒有輻射到遠場。為了進一步探索縫翼低頻窄帶噪聲源的分布特征,本文對原始瞬態(tài)壓力脈動數(shù)據(jù)進行處理,濾除高頻和對低頻窄帶貢獻較小的成分,并利用相關性分析,探討縫翼窄帶噪聲源的分布特征。

    2.3 低頻窄帶噪聲源分布

    本文在剪切層與縫道區(qū)域分別設置17和5個監(jiān)測點(見圖11),首先從流場分析中獲得這些點的瞬態(tài)壓力幅值;然后利用相關性與傅里葉變換分析,建立剪切層與縫道區(qū)域聲源分布規(guī)律。剪切層總長S=2.20C,縫道總長為T=0.035C。

    圖11 各監(jiān)測點位置及排列順序

    圖12為剪切層(奇數(shù)監(jiān)測點)和縫道監(jiān)測點的壓力脈動頻譜,其中1 343,1 984與2 653 Hz峰值頻率對應遠場的前3個峰值。S17在整個頻率范圍內均具有最大的頻譜幅值,而與其相近的S15、S13和S11等點的頻譜幅值逐漸降低。S13頻譜具有特殊性,其前2個峰值與相鄰點的頻譜幅值變換規(guī)律呈現(xiàn)明顯差異,而S13正是處于剪切層渦結構的驟變區(qū)域;同時,縫道區(qū)域中的T1-T5與S15、S17點則具有類似的窄帶頻譜特性。

    圖13為各監(jiān)測點對應縫翼遠場前2個窄帶峰值的頻譜值。該圖橫坐標為縫翼凹腔與縫道典型位置。從左向右為:P2,P1,S1~S17,P3,P4,T1~T5,P5等??梢钥闯?在縫翼尾源的P4點,剪切層后段S16,S17點,以及主翼前緣P5點,其2個峰值譜值均較大,說明窄帶噪聲源與這些位置密切相關;據(jù)此可推斷出噪聲源主要集中在縫道區(qū)域。

    圖13 前2個峰值各監(jiān)測點的壓力頻譜幅值

    為了對窄帶噪聲源特性進一步分析,利用EMD方法,對監(jiān)測點(P4, P5,S16,S17,T3-T5)的源噪聲信號進行分解,獲得其本征模態(tài)函數(shù)(IMF)分量;然后采用低通重構方法去除對窄帶噪聲影響較小的高頻IMF分量,并將IMF分量幅值較低的濾除掉。相關步驟如下:

    1) 對各監(jiān)測點噪聲進行EMD分解,獲得各點IMF分量;

    2) 計算各監(jiān)測點中每個IMF分量與其原始信號的相關系數(shù)ρimfi(t)(i表示第i個IMF分量);

    4) 基于步驟3)中保留的IMF分量,進一步選取那些幅值較大的分量(選擇Ai>2,Ai表示第i各IMF分量的振幅);

    5) 將最終保留的IMF進行信號重構,即將選擇的IMF分量和殘差(rn)累加獲得新信號x′(t)

    (12)

    新生成的信號x′(t),主要含有與低頻窄帶噪聲相關的信息。下面對各個監(jiān)測點的x′(t)進行相關性分析,以此揭示窄帶噪聲的分布特性。

    為了確定低頻噪聲源的分布形式,將上述各點重構信號分別與特征點P4(頻譜幅值較大點)進行互相關分析,以獲得各點的互相關函數(shù)與特征延遲時間。結合圖17和壓力脈動頻譜(圖5和圖12),可發(fā)現(xiàn):

    1) P3與P4相關函數(shù)值較大,但其曲線較平直,這是P3頻譜幅值較小所致。P4較P3有延遲,顯示P3處于剪切層渦流與縫翼尾緣碰撞形成的湍流微團核心區(qū)域中,其脈動較小,非聲源點。而P4為大尺度渦碰撞后的壓力脈沖及離散渦結構,其中渦結構的驟變存在能量轉換及釋放,對流場影響較大,并會導致噪聲輻射;

    2) T4與P4呈弱相關,顯示出該點與尾緣區(qū)域具有不同的流動結構。T5與P4相關性較強,但流場結構卻存在差異,P4主要為大渦結構,而T5則為小渦或層流混合結構,同時T5還受到與P4之間狹縫共振的影響;

    3) P5雖與P4相關性較低,而其相關函數(shù)規(guī)律則與之相似,同時具有正延遲,反映出主翼前緣點的流動結構雖與縫翼尾緣不同,但卻受到尾緣脈動的直接影響;

    4) S16,S17與P4均表現(xiàn)出強相關性。S16對P4延遲為負,而S17對P4延遲為正,同時S16與P3延遲相近,因而在模態(tài)密集區(qū)各點脈動規(guī)律呈現(xiàn)較為復雜的特性。但由頻譜及相關性分析可知,S16,S17及P4的流場結構特性是低頻窄帶噪聲產生的主要原因。同時S16區(qū)域部分渦卷入低速回流區(qū),也將P3區(qū)域的離散渦脈動能量帶入回流區(qū)。

    上述特征點是根據(jù)縫翼流場典型位置確定的,未必能完全表示低頻噪聲源區(qū)域。下面將探索性地利用特征點延遲時間進一步分析聲源位置的分布特征。

    假設S(x,y,z)為聲源點,Ri為監(jiān)測點,τij為Ri和Rj的延遲時間(τij=τi-τj),τi代表S到Ri的聲傳播時間。因此有關系式:(cτi)2=(xi-x)2+(yi-y)2+(zi-z)2,其中,c表示聲速。建立方程組,求解可得噪聲源位置(見圖14)。

    圖14 主要噪聲源位置示意圖

    由圖14可知,噪聲源主要分布在碰撞區(qū)、剪切層下游、縫道以及部分回流區(qū)。其中,碰撞點附近的剪切層區(qū)域聲源較密集。對比DMD模態(tài)云圖,發(fā)現(xiàn)計算所得聲源點主要集中在模態(tài)云圖中壓力脈動較大的區(qū)域;同時,在縫道以及主翼前緣均有聲源分布,而該結果與KONIG等[5]的研究結論基本一致。這也反映出低頻窄帶噪聲與剪切層上渦結構周期性撞擊壁面引起縫翼尾緣與主翼邊界層之間的壓力脈沖有密切關系。

    根據(jù)相關性分析及聲源位置分布規(guī)律可知,剪切層下游與壓力脈動頻譜峰值點呈現(xiàn)強相關,對低頻窄帶噪聲輻射水平影響較大;因此可基于剪切層干擾的思想提出縫翼尖端變形等技術進行縫翼降噪方法的研究。課題組后續(xù)將展開縫翼前緣結構調整進行縫翼降噪的研究。

    3 結 論

    本文利用POD和DMD方法分析縫翼瞬態(tài)流場結構特征,并基于EMD法和相關性分析研究了低頻窄帶噪聲的分布規(guī)律。主要結論有:

    1) POD第一階模態(tài)區(qū)域顯示,約98%的模態(tài)能量主要集中在縫翼尾緣碰撞區(qū)附近,該區(qū)域是噪聲產生的重要區(qū)域,其頻率組成具有寬頻特性;

    2) DMD各階模態(tài)特征表明,低頻窄帶聲源不僅集中在縫翼尾緣附近區(qū)域,在剪切層下游,以及回流區(qū)也存在分布較廣的低頻窄帶聲源;

    3) 結合EMD分解與互相關分析,低頻窄帶噪聲源不僅存在于特征模態(tài)脈動較大區(qū)域,縫道及主翼前部區(qū)域也存在噪聲源分布。該結論表明,剪切層撞擊縫翼尾緣壁面所形成的離散渦與主翼前緣邊界層相互耦合是低頻窄帶聲源產生的重要原因。在此基礎上,可基于剪切層干擾的思想提出新的縫翼降噪方法。

    猜你喜歡
    尾緣噪聲源窄帶
    基于強化換熱的偏斜尾緣設計
    能源工程(2021年1期)2021-04-13 02:05:46
    汽車后視鏡-A柱區(qū)域氣動噪聲源特征識別
    汽車工程(2021年12期)2021-03-08 02:34:18
    熱軋窄帶鋼Q345B微合金化生產實踐
    山東冶金(2019年1期)2019-03-30 01:34:54
    翼型湍流尾緣噪聲半經驗預測公式改進
    具有尾緣襟翼的風力機動力學建模與恒功率控制
    無線通信中頻線路窄帶臨界調試法及其應用
    電子制作(2017年19期)2017-02-02 07:08:38
    一種基于相位增量隨機化的寬帶噪聲源產生技術
    雷達與對抗(2015年3期)2015-12-09 02:38:55
    鈍尾緣葉片三維建模方法的研究
    基于壓縮感知的窄帶干擾重構與消除
    電視技術(2014年19期)2014-03-11 15:38:13
    利用源強聲輻射模態(tài)識別噪聲源
    纵有疾风起免费观看全集完整版| 久久影院123| 波多野结衣一区麻豆| 欧美xxⅹ黑人| videosex国产| 18禁国产床啪视频网站| 男女无遮挡免费网站观看| 亚洲国产精品国产精品| 韩国高清视频一区二区三区| av不卡在线播放| 一本大道久久a久久精品| 一本大道久久a久久精品| 老司机亚洲免费影院| 亚洲精品av麻豆狂野| 国产精品一区二区在线不卡| 黑人猛操日本美女一级片| 777米奇影视久久| 成人国产av品久久久| 日韩一区二区视频免费看| 亚洲精品一区蜜桃| 丝袜在线中文字幕| 久久久久久人人人人人| 日韩中文字幕欧美一区二区 | 最近手机中文字幕大全| 国精品久久久久久国模美| 精品少妇久久久久久888优播| 精品少妇久久久久久888优播| 中文字幕av电影在线播放| 久久97久久精品| 在线观看国产h片| 欧美另类一区| 国产激情久久老熟女| 亚洲精品国产av成人精品| 精品福利永久在线观看| 在线天堂最新版资源| 大陆偷拍与自拍| 一级毛片黄色毛片免费观看视频| 狠狠婷婷综合久久久久久88av| 免费观看人在逋| 男女免费视频国产| 日韩人妻精品一区2区三区| 你懂的网址亚洲精品在线观看| 亚洲久久久国产精品| 国产一区二区三区综合在线观看| 人妻一区二区av| 亚洲在久久综合| 人体艺术视频欧美日本| 丝袜美腿诱惑在线| 免费高清在线观看日韩| 国产成人精品久久二区二区91 | 大片电影免费在线观看免费| 久久午夜综合久久蜜桃| 久久热在线av| 久久人人爽av亚洲精品天堂| 天堂8中文在线网| 人人妻人人添人人爽欧美一区卜| 国产一区二区 视频在线| www.熟女人妻精品国产| 激情五月婷婷亚洲| 久久影院123| 亚洲国产精品一区三区| 青春草亚洲视频在线观看| 成人18禁高潮啪啪吃奶动态图| 国产精品三级大全| 国产精品 欧美亚洲| 国产男女内射视频| 欧美另类一区| 人人妻人人添人人爽欧美一区卜| 久久国产亚洲av麻豆专区| 91精品伊人久久大香线蕉| 天天躁日日躁夜夜躁夜夜| 日韩免费高清中文字幕av| 在线 av 中文字幕| 免费黄色在线免费观看| 亚洲精品国产区一区二| 99久久人妻综合| 欧美97在线视频| 亚洲av成人精品一二三区| av网站免费在线观看视频| 美女大奶头黄色视频| 十八禁网站网址无遮挡| 在线亚洲精品国产二区图片欧美| 午夜日韩欧美国产| 午夜激情久久久久久久| 18在线观看网站| 男女床上黄色一级片免费看| 中文字幕制服av| 在现免费观看毛片| 久久久精品94久久精品| 国产熟女午夜一区二区三区| 久久国产精品大桥未久av| 久久婷婷青草| 亚洲精品中文字幕在线视频| 王馨瑶露胸无遮挡在线观看| 99久久综合免费| 2021少妇久久久久久久久久久| 精品少妇一区二区三区视频日本电影 | 又大又爽又粗| 欧美亚洲 丝袜 人妻 在线| 精品国产乱码久久久久久小说| 亚洲免费av在线视频| 男女边吃奶边做爰视频| 成人18禁高潮啪啪吃奶动态图| 9191精品国产免费久久| 欧美日韩av久久| 日韩 亚洲 欧美在线| 欧美成人精品欧美一级黄| 看免费av毛片| 免费日韩欧美在线观看| 欧美成人精品欧美一级黄| 国产精品麻豆人妻色哟哟久久| 蜜桃在线观看..| 国产精品.久久久| 国产成人精品久久二区二区91 | 人人妻人人添人人爽欧美一区卜| 超碰成人久久| 建设人人有责人人尽责人人享有的| 天天躁夜夜躁狠狠躁躁| 亚洲av欧美aⅴ国产| av网站免费在线观看视频| 制服诱惑二区| 久久久国产一区二区| 午夜免费鲁丝| 国产精品久久久久久精品电影小说| 在线观看www视频免费| 天天躁狠狠躁夜夜躁狠狠躁| 日韩,欧美,国产一区二区三区| 成人亚洲欧美一区二区av| 久久久久久久久久久久大奶| 99精品久久久久人妻精品| 可以免费在线观看a视频的电影网站 | 精品第一国产精品| xxx大片免费视频| 少妇被粗大猛烈的视频| 亚洲一码二码三码区别大吗| www.av在线官网国产| 青春草视频在线免费观看| 999久久久国产精品视频| 亚洲欧美日韩另类电影网站| 欧美日韩国产mv在线观看视频| 欧美精品av麻豆av| 亚洲一级一片aⅴ在线观看| 国产国语露脸激情在线看| 一级片免费观看大全| 啦啦啦 在线观看视频| 午夜福利一区二区在线看| 高清在线视频一区二区三区| 亚洲五月色婷婷综合| 人妻人人澡人人爽人人| 综合色丁香网| 高清av免费在线| 日本vs欧美在线观看视频| 多毛熟女@视频| 亚洲 欧美一区二区三区| 夜夜骑夜夜射夜夜干| 欧美亚洲日本最大视频资源| 在线观看一区二区三区激情| 男女床上黄色一级片免费看| 欧美在线一区亚洲| 看非洲黑人一级黄片| 国产精品免费大片| 在线观看免费高清a一片| 在线免费观看不下载黄p国产| 老熟女久久久| 国产伦理片在线播放av一区| 精品免费久久久久久久清纯 | 欧美黑人欧美精品刺激| 狠狠婷婷综合久久久久久88av| 黄频高清免费视频| 亚洲美女黄色视频免费看| 亚洲国产中文字幕在线视频| 无限看片的www在线观看| 男人舔女人的私密视频| 亚洲,欧美精品.| 男的添女的下面高潮视频| 久久久久国产一级毛片高清牌| av免费观看日本| 久久久久国产一级毛片高清牌| 亚洲欧美一区二区三区久久| 丝袜在线中文字幕| 免费在线观看完整版高清| 男女下面插进去视频免费观看| 亚洲欧洲国产日韩| 久久久久久久久久久免费av| 国产免费一区二区三区四区乱码| 2021少妇久久久久久久久久久| 国产精品久久久av美女十八| 九草在线视频观看| 日韩大码丰满熟妇| 精品久久蜜臀av无| 日日爽夜夜爽网站| 亚洲欧美精品综合一区二区三区| 午夜福利影视在线免费观看| 亚洲熟女毛片儿| 日本91视频免费播放| 色94色欧美一区二区| 国产不卡av网站在线观看| 久久人妻熟女aⅴ| 国产一区二区在线观看av| 久久青草综合色| 水蜜桃什么品种好| 国产成人系列免费观看| 午夜福利在线免费观看网站| 亚洲av中文av极速乱| 成人国语在线视频| 亚洲 欧美一区二区三区| 久久久精品区二区三区| 欧美变态另类bdsm刘玥| av.在线天堂| 国产麻豆69| 卡戴珊不雅视频在线播放| 丝袜在线中文字幕| 天天躁夜夜躁狠狠躁躁| 秋霞伦理黄片| 免费观看av网站的网址| 国产麻豆69| 卡戴珊不雅视频在线播放| 午夜精品国产一区二区电影| 国产在线视频一区二区| 中文字幕人妻熟女乱码| 成人漫画全彩无遮挡| 国产精品一区二区精品视频观看| av在线播放精品| 另类精品久久| 亚洲精品一二三| 国产精品 国内视频| 中文精品一卡2卡3卡4更新| 夜夜骑夜夜射夜夜干| 午夜福利视频在线观看免费| 18禁动态无遮挡网站| 亚洲精品一二三| 国产精品国产三级专区第一集| 2018国产大陆天天弄谢| av不卡在线播放| 国产日韩欧美亚洲二区| videosex国产| 亚洲av成人精品一二三区| 亚洲,一卡二卡三卡| 久久久久久免费高清国产稀缺| 亚洲久久久国产精品| 亚洲人成电影观看| 伦理电影免费视频| 97精品久久久久久久久久精品| 久久精品熟女亚洲av麻豆精品| 亚洲成人av在线免费| 操美女的视频在线观看| 久久97久久精品| 晚上一个人看的免费电影| 九草在线视频观看| 69精品国产乱码久久久| 午夜免费鲁丝| 热re99久久精品国产66热6| 人妻一区二区av| 久久99精品国语久久久| 亚洲精品久久成人aⅴ小说| 丰满乱子伦码专区| 中文乱码字字幕精品一区二区三区| 精品亚洲成国产av| 欧美日韩综合久久久久久| 精品久久久久久电影网| 国精品久久久久久国模美| 亚洲国产成人一精品久久久| 国产一区二区三区综合在线观看| 午夜日本视频在线| 国产人伦9x9x在线观看| 韩国av在线不卡| 午夜福利,免费看| 超色免费av| 久久人人爽av亚洲精品天堂| 国产精品国产三级国产专区5o| 国产不卡av网站在线观看| 欧美日韩成人在线一区二区| 9色porny在线观看| 日韩一区二区三区影片| 80岁老熟妇乱子伦牲交| 成人毛片60女人毛片免费| www.精华液| 少妇人妻精品综合一区二区| 亚洲视频免费观看视频| 十八禁人妻一区二区| 综合色丁香网| 999精品在线视频| 亚洲成av片中文字幕在线观看| 69精品国产乱码久久久| 丝袜美足系列| 日本91视频免费播放| 国产成人精品无人区| 自拍欧美九色日韩亚洲蝌蚪91| 曰老女人黄片| 欧美变态另类bdsm刘玥| 又大又爽又粗| 国产一区二区 视频在线| 丁香六月欧美| 免费黄网站久久成人精品| 成人亚洲欧美一区二区av| 久久久久视频综合| 最近2019中文字幕mv第一页| 男女无遮挡免费网站观看| 欧美精品一区二区大全| 丰满乱子伦码专区| 啦啦啦啦在线视频资源| 深夜精品福利| 街头女战士在线观看网站| 热re99久久国产66热| 国产一卡二卡三卡精品 | 女性被躁到高潮视频| 亚洲精品一区蜜桃| 一级黄片播放器| 亚洲精品久久成人aⅴ小说| 欧美精品人与动牲交sv欧美| 激情五月婷婷亚洲| 青春草国产在线视频| 岛国毛片在线播放| 亚洲熟女精品中文字幕| 最近中文字幕高清免费大全6| 久久精品久久久久久久性| 人妻一区二区av| 在线观看免费午夜福利视频| 黄色怎么调成土黄色| 18禁动态无遮挡网站| 99国产综合亚洲精品| 视频区图区小说| av卡一久久| 新久久久久国产一级毛片| 免费观看av网站的网址| 自拍欧美九色日韩亚洲蝌蚪91| 午夜激情av网站| 熟女av电影| 人人妻人人澡人人看| 日韩免费高清中文字幕av| h视频一区二区三区| 十八禁人妻一区二区| 午夜激情久久久久久久| 久久韩国三级中文字幕| 日本欧美视频一区| 亚洲熟女精品中文字幕| 中文字幕av电影在线播放| 一个人免费看片子| 日韩人妻精品一区2区三区| 久久99热这里只频精品6学生| www.自偷自拍.com| 男女之事视频高清在线观看 | 最近中文字幕2019免费版| 中文字幕色久视频| 亚洲欧美精品综合一区二区三区| 国产1区2区3区精品| av线在线观看网站| 久久女婷五月综合色啪小说| 成年女人毛片免费观看观看9 | 黄色视频不卡| 最近最新中文字幕免费大全7| 侵犯人妻中文字幕一二三四区| 亚洲,欧美精品.| 交换朋友夫妻互换小说| 免费av中文字幕在线| 亚洲国产成人一精品久久久| 叶爱在线成人免费视频播放| 国产一区亚洲一区在线观看| 亚洲三区欧美一区| 久久久国产精品麻豆| 精品午夜福利在线看| 在线观看www视频免费| 老汉色av国产亚洲站长工具| 日本vs欧美在线观看视频| 国产精品嫩草影院av在线观看| 国产伦理片在线播放av一区| 久久久精品国产亚洲av高清涩受| 日本欧美国产在线视频| 18禁国产床啪视频网站| 亚洲国产av影院在线观看| 水蜜桃什么品种好| 成年人免费黄色播放视频| 爱豆传媒免费全集在线观看| 青春草国产在线视频| 波野结衣二区三区在线| www.精华液| 美女视频免费永久观看网站| 中文字幕人妻丝袜制服| 久久久久久久久免费视频了| 丁香六月天网| 丝袜人妻中文字幕| 久久99热这里只频精品6学生| 精品亚洲成a人片在线观看| 亚洲激情五月婷婷啪啪| www.熟女人妻精品国产| 婷婷成人精品国产| 欧美日韩综合久久久久久| 久久午夜综合久久蜜桃| 精品亚洲成国产av| 亚洲人成电影观看| 欧美日韩国产mv在线观看视频| 人妻 亚洲 视频| 操出白浆在线播放| 丝袜脚勾引网站| 大香蕉久久网| 国产精品 欧美亚洲| 丝袜在线中文字幕| 亚洲av日韩在线播放| 啦啦啦在线免费观看视频4| 成年动漫av网址| bbb黄色大片| 国产精品免费视频内射| 90打野战视频偷拍视频| √禁漫天堂资源中文www| 成年av动漫网址| 久久免费观看电影| 国产av国产精品国产| avwww免费| 国产精品无大码| 国产片特级美女逼逼视频| 99国产精品免费福利视频| 亚洲精品乱久久久久久| 国产精品久久久久久精品古装| 国产伦人伦偷精品视频| 制服人妻中文乱码| 韩国精品一区二区三区| 一级,二级,三级黄色视频| 美女福利国产在线| 秋霞伦理黄片| 日韩不卡一区二区三区视频在线| 曰老女人黄片| 国产日韩欧美亚洲二区| 欧美日韩成人在线一区二区| 午夜福利影视在线免费观看| 精品少妇黑人巨大在线播放| 最新在线观看一区二区三区 | 久久 成人 亚洲| 久热这里只有精品99| 十分钟在线观看高清视频www| 一区二区三区激情视频| 精品免费久久久久久久清纯 | 亚洲av国产av综合av卡| 一二三四中文在线观看免费高清| 一区福利在线观看| 精品国产一区二区久久| 激情视频va一区二区三区| av网站免费在线观看视频| av卡一久久| 熟女av电影| 熟女少妇亚洲综合色aaa.| av线在线观看网站| 中文欧美无线码| 一本大道久久a久久精品| 欧美激情极品国产一区二区三区| 丰满乱子伦码专区| av在线app专区| 久久精品久久久久久噜噜老黄| 国产av一区二区精品久久| 日本vs欧美在线观看视频| 久久久久久免费高清国产稀缺| 国产成人精品久久二区二区91 | 91精品国产国语对白视频| 亚洲欧美清纯卡通| 国产欧美亚洲国产| 国产成人免费观看mmmm| 超碰97精品在线观看| 色精品久久人妻99蜜桃| 夜夜骑夜夜射夜夜干| 国产成人精品福利久久| 精品第一国产精品| 男人添女人高潮全过程视频| 亚洲av男天堂| 自线自在国产av| 成人手机av| 水蜜桃什么品种好| 丝袜在线中文字幕| 久久99精品国语久久久| 建设人人有责人人尽责人人享有的| 男的添女的下面高潮视频| 一二三四中文在线观看免费高清| 精品少妇黑人巨大在线播放| 搡老乐熟女国产| 久久精品aⅴ一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看 | 一边摸一边抽搐一进一出视频| 黄色一级大片看看| 亚洲国产中文字幕在线视频| 国产无遮挡羞羞视频在线观看| 国产伦人伦偷精品视频| 一区二区三区乱码不卡18| 国产欧美亚洲国产| 看免费av毛片| 午夜福利,免费看| 青春草亚洲视频在线观看| 宅男免费午夜| 国产精品三级大全| 成人午夜精彩视频在线观看| 国产精品 欧美亚洲| 亚洲成人国产一区在线观看 | 天天躁日日躁夜夜躁夜夜| 99国产精品免费福利视频| 校园人妻丝袜中文字幕| 国产精品av久久久久免费| 侵犯人妻中文字幕一二三四区| 看非洲黑人一级黄片| 一级爰片在线观看| 成人漫画全彩无遮挡| 夫妻性生交免费视频一级片| 777久久人妻少妇嫩草av网站| 日韩av在线免费看完整版不卡| 黄色 视频免费看| 亚洲成国产人片在线观看| 久久久欧美国产精品| 国产成人免费无遮挡视频| 亚洲国产精品一区二区三区在线| 久久久久国产精品人妻一区二区| 免费在线观看黄色视频的| 精品一区二区三卡| 日本av手机在线免费观看| 国产精品av久久久久免费| 少妇人妻精品综合一区二区| 亚洲国产av新网站| 国产精品久久久人人做人人爽| av国产精品久久久久影院| 欧美日韩一区二区视频在线观看视频在线| 青春草国产在线视频| 日日撸夜夜添| 黄片小视频在线播放| 嫩草影院入口| 麻豆av在线久日| 日韩制服丝袜自拍偷拍| 欧美另类一区| 国产精品.久久久| 欧美日韩精品网址| 99热国产这里只有精品6| 国产淫语在线视频| 国产成人免费无遮挡视频| 天天操日日干夜夜撸| 久久久久久久大尺度免费视频| 成年女人毛片免费观看观看9 | 亚洲欧美一区二区三区黑人| 婷婷成人精品国产| 成人三级做爰电影| 国产探花极品一区二区| 亚洲精品日本国产第一区| 日本色播在线视频| 国产免费又黄又爽又色| 国精品久久久久久国模美| 亚洲精品aⅴ在线观看| 亚洲精品国产av蜜桃| 大片电影免费在线观看免费| 啦啦啦 在线观看视频| 国产老妇伦熟女老妇高清| 日本午夜av视频| 免费少妇av软件| 一区二区三区四区激情视频| 亚洲av欧美aⅴ国产| 免费黄网站久久成人精品| 国产精品 欧美亚洲| 人人妻人人澡人人爽人人夜夜| 国产精品99久久99久久久不卡 | 美女午夜性视频免费| 一区二区三区四区激情视频| 伊人亚洲综合成人网| 久久人人97超碰香蕉20202| 丰满乱子伦码专区| 99精国产麻豆久久婷婷| 久久精品久久久久久久性| 精品久久蜜臀av无| 国产欧美亚洲国产| 日韩 亚洲 欧美在线| 色婷婷av一区二区三区视频| 亚洲第一青青草原| 不卡av一区二区三区| 欧美人与性动交α欧美精品济南到| 美女中出高潮动态图| 一二三四中文在线观看免费高清| 日韩av在线免费看完整版不卡| 十八禁高潮呻吟视频| 国产精品女同一区二区软件| 大话2 男鬼变身卡| 18禁观看日本| www.熟女人妻精品国产| 国产熟女欧美一区二区| 久久久久久免费高清国产稀缺| 自拍欧美九色日韩亚洲蝌蚪91| 国产一区二区激情短视频 | 免费日韩欧美在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 精品人妻在线不人妻| 欧美乱码精品一区二区三区| 欧美日韩视频精品一区| 在线观看免费视频网站a站| 色播在线永久视频| av在线观看视频网站免费| 日本欧美视频一区| 久久韩国三级中文字幕| 超色免费av| 欧美在线黄色| 91精品伊人久久大香线蕉| 只有这里有精品99| 精品一区二区三卡| 大片免费播放器 马上看| 国产黄频视频在线观看| 国产色婷婷99| 亚洲婷婷狠狠爱综合网| 国产精品偷伦视频观看了| 国产xxxxx性猛交| 嫩草影院入口| 91aial.com中文字幕在线观看| 一区福利在线观看| 日韩欧美精品免费久久| 香蕉国产在线看| 秋霞伦理黄片| av电影中文网址| 如日韩欧美国产精品一区二区三区| 亚洲精品久久成人aⅴ小说| av国产久精品久网站免费入址| 制服诱惑二区| 美女大奶头黄色视频| 极品少妇高潮喷水抽搐| 久久久久久久精品精品| 日韩av在线免费看完整版不卡| 国产成人系列免费观看| 午夜福利视频精品|