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

    基于流固耦合的硬涂層整體葉盤振動(dòng)特性研究

    2023-08-17 01:46:42楊錚鑫孫榮城鄭偉黨鵬飛
    機(jī)床與液壓 2023年14期
    關(guān)鍵詞:葉盤共振固有頻率

    楊錚鑫,孫榮城,鄭偉,黨鵬飛

    (沈陽化工大學(xué)機(jī)械與動(dòng)力工程學(xué)院,遼寧沈陽 110142)

    0 前言

    航空發(fā)動(dòng)機(jī)是一種結(jié)構(gòu)復(fù)雜的高速旋轉(zhuǎn)機(jī)械,在復(fù)雜多變的工況下運(yùn)作,導(dǎo)致內(nèi)部結(jié)構(gòu)出現(xiàn)了不同方式的振動(dòng)損壞[1]。整體葉盤作為航空發(fā)動(dòng)機(jī)中核心元件之一,在轉(zhuǎn)速快、溫度高和氣壓較大的工況下運(yùn)作,受到離心載荷、溫度載荷、氣動(dòng)載荷等多載荷的共同作用[2]。隨著航空航天技術(shù)對(duì)整體葉盤性能的不斷優(yōu)化,整體葉盤向著質(zhì)量輕、推重比大、體積小等方向發(fā)展。相對(duì)于傳統(tǒng)的榫槽結(jié)構(gòu)葉盤,整體葉盤缺少了榫槽結(jié)構(gòu)的摩擦效應(yīng),導(dǎo)致結(jié)構(gòu)的阻尼性能較低[3]。隨著硬涂層阻尼技術(shù)的發(fā)展,由金屬基和陶瓷基合成的硬涂層不但優(yōu)化了整體葉盤的阻尼性能,還提高了結(jié)構(gòu)的阻尼能力[4]。在復(fù)雜多變的載荷作用下,為了避免整體葉盤出現(xiàn)振動(dòng)導(dǎo)致結(jié)構(gòu)損壞的問題,以及提高結(jié)構(gòu)系統(tǒng)的可靠性,有必要利用有限元軟件對(duì)整體葉盤的振動(dòng)特性進(jìn)行分析。

    魏武國等[5]采用分塊Lanczos法找出了高壓壓氣機(jī)轉(zhuǎn)子葉片應(yīng)力較大的薄弱區(qū)域以及氣流尾跡造成的共振。張露、卞祥德[6]通過坎貝爾圖找出發(fā)生雙扭復(fù)合共振的位置,采取削角的方法降低了振動(dòng)應(yīng)力。金業(yè)壯等[7]采用優(yōu)化的八節(jié)點(diǎn)超參數(shù)單元法分析研究了可能發(fā)生的危險(xiǎn)共振情況,為轉(zhuǎn)子的故障檢測(cè)以及預(yù)報(bào)提供了依據(jù)。寇海軍等[8]利用流體力學(xué)求解葉片表面的氣動(dòng)載荷,在轉(zhuǎn)子正常轉(zhuǎn)速范圍內(nèi)找到10個(gè)共振轉(zhuǎn)速點(diǎn),表明了共振應(yīng)力峰值與臨界轉(zhuǎn)速并無正相關(guān)性。楊錚鑫等[9]研究了改變NiCrAlY硬涂層的涂敷方式對(duì)整體葉盤減振效果的影響,發(fā)現(xiàn)硬涂層技術(shù)可以降低整體葉盤共振點(diǎn)的數(shù)量。吳正人等[10]利用CFX和ANSYS對(duì)離心風(fēng)機(jī)葉輪進(jìn)行了流固耦合模擬,發(fā)現(xiàn)葉輪的固有頻率部分落入局部的共振區(qū)域。

    本文作者考慮某航空噴氣發(fā)動(dòng)機(jī)一級(jí)轉(zhuǎn)子工作的實(shí)際工況,對(duì)整體葉盤的工況進(jìn)行設(shè)計(jì)。選用鎳基高溫合金GH4169作為整體葉盤材料,對(duì)整體葉盤葉片的壓力面涂敷NiCrAlY硬涂層。利用ANSYS Workbench中Fluent組件求解整體葉盤的氣動(dòng)載荷,采用流固耦合方法將氣動(dòng)載荷導(dǎo)入預(yù)應(yīng)力分析模塊,再采用預(yù)應(yīng)力模態(tài)分析法進(jìn)行求解。分析在氣動(dòng)載荷情況下,硬涂層對(duì)整體葉盤固有頻率的影響;分析在離心載荷、氣動(dòng)載荷以及2種載荷共同作用下對(duì)涂敷硬涂層的整體葉盤固有頻率和等效應(yīng)力的影響;選取具有特征性的振型以及葉片的變形情況進(jìn)行分析總結(jié);最后通過繪制工況下無硬涂層整體葉盤和有硬涂層整體葉盤的坎貝爾圖,分析共振點(diǎn)對(duì)應(yīng)的共振轉(zhuǎn)速。

    1 建模及參數(shù)

    1.1 整體葉盤建模

    基于Creo軟件建立整體葉盤有限元模型,并導(dǎo)入ANSYS Geometry中進(jìn)行模型處理。圖1所示為整體葉盤有限元模型,圖2所示為葉片壓力面涂敷NiCrAlY硬涂層的葉片有限元模型。表1為硬涂層整體葉盤的幾何參數(shù)。

    表1 硬涂層整體葉盤的幾何參數(shù)

    圖1 整體葉盤示意 圖2 硬涂層葉片

    整體葉盤的材料選用沉淀強(qiáng)化鎳基高溫合金——GH4169合金,硬涂層材料選用NiCrAlY。以上2種材料參數(shù)通過查閱中國航空材料手冊(cè)獲得,表2為硬涂層整體葉盤的材料參數(shù)。

    表2 硬涂層整體葉盤的材料參數(shù)

    1.2 整體葉盤流場(chǎng)模型

    通過查閱文獻(xiàn)[11]可知,渦輪噴氣發(fā)動(dòng)機(jī)一級(jí)轉(zhuǎn)子的實(shí)際工況,整體葉盤的轉(zhuǎn)速范圍為30 000~39 600 r/min,文中設(shè)計(jì)整體葉盤工況見表3。

    表3 整體葉盤工況設(shè)計(jì)

    整體葉盤流場(chǎng)物理模型如圖3所示,利用Boolean對(duì)整體葉盤、內(nèi)流域和外流域進(jìn)行布爾操作,再利用Fluent前處理網(wǎng)格工具對(duì)模型整體采用四面體單元進(jìn)行網(wǎng)格劃分。為了提高計(jì)算結(jié)果的精度,對(duì)整體葉盤流場(chǎng)物理模型和流域接觸邊界進(jìn)行網(wǎng)格局部加密。利用Fluent對(duì)5種工況下的整體葉盤表面氣動(dòng)載荷進(jìn)行求解。模型邊界條件采用壓力進(jìn)口和壓力出口,選擇標(biāo)準(zhǔn)的κ-ε湍流模型和標(biāo)準(zhǔn)壁面函數(shù),求解時(shí)選用SIMPLE算法和二階迎風(fēng)格式。為了驗(yàn)證網(wǎng)格數(shù)量對(duì)整體計(jì)算結(jié)果的影響,在481萬網(wǎng)格基礎(chǔ)上增加30%(625萬網(wǎng)格),求解整體葉盤表面氣動(dòng)載荷,相差1.3%以內(nèi),所以認(rèn)為481萬網(wǎng)格滿足求解要求。根據(jù)上述方法劃分模型為847 884個(gè)節(jié)點(diǎn)、4 819 430個(gè)單元。

    圖3 流場(chǎng)模型

    2 求解模型理論

    2.1 單向流固耦合法

    單向流固耦合是流固耦合分析中分離解法之一,其不需要耦合控制方程,是通過流固交界面(FS Interface)把計(jì)算結(jié)果單向傳遞到固體。發(fā)動(dòng)機(jī)整體葉盤在高速的氣流中運(yùn)轉(zhuǎn),高速氣流會(huì)對(duì)整體葉盤的振動(dòng)性能產(chǎn)生影響。根據(jù)整體葉盤在流場(chǎng)中的工作狀態(tài),文中采用標(biāo)準(zhǔn)的κ-ε湍流模型,如下:

    YM+Sκ

    (1)

    (2)

    其中:

    (3)

    (4)

    (5)

    (6)

    (7)

    式中:湍流黏度系數(shù)μt可表示成κ與ε的函數(shù);Gκ是速度梯度產(chǎn)生的湍動(dòng)能;Gb是浮力影響產(chǎn)生的湍動(dòng)能;YM是可壓縮湍流脈動(dòng)膨脹對(duì)總的耗散率的影響;經(jīng)驗(yàn)系數(shù)默認(rèn)C1ε=1.44、C2ε=1.92、C3ε=0.09;默認(rèn)的湍動(dòng)能和湍動(dòng)耗散率對(duì)應(yīng)的普朗特系數(shù)分別是σκ=1.0、σε=1.3;湍動(dòng)普朗特系數(shù)默認(rèn)Prt=0.85;gi是在i方向上的重力加速度;β是熱膨脹系數(shù);Mt是湍動(dòng)馬赫數(shù);a是聲速。

    2.2 預(yù)應(yīng)力模態(tài)分析法

    文中將整體葉盤在流場(chǎng)中受到的氣動(dòng)載荷導(dǎo)入靜力分析模塊,通過對(duì)整體葉盤進(jìn)行預(yù)應(yīng)力分析后,在預(yù)應(yīng)力基礎(chǔ)上進(jìn)行模態(tài)分析。

    結(jié)構(gòu)系統(tǒng)的動(dòng)力學(xué)方程:

    (8)

    式中:M是系統(tǒng)質(zhì)量矩陣;C是系統(tǒng)阻尼矩陣;K是系統(tǒng)剛度矩陣。

    (9)

    x=Xsin(ωt+α)

    (10)

    將式(10)代入式(9)可得到:

    (K-ω2M)X=0

    (11)

    用A表示K-ω2M,令detA=0,得:

    (12)

    求解得到整體葉盤的固有頻率:

    (13)

    3 特性分析

    3.1 固有頻率分析

    由于整體葉盤具有周期對(duì)稱性,對(duì)整體葉盤采用循環(huán)對(duì)稱分析方式進(jìn)行有限元分析。分別計(jì)算了整體葉盤和硬涂層整體葉盤在氣動(dòng)載荷、離心載荷以及2種載荷共同作用下的前72階頻率,其中1—18階為葉片一階頻率,19—36階為葉片二階頻率,37—54階為葉片三階頻率,55—72階為葉片四階頻率。

    分別求解了整體葉盤和硬涂層整體葉盤在氣動(dòng)載荷作用下的固有頻率,見表4。對(duì)整體葉盤和硬涂層整體葉盤的固有頻率進(jìn)行比較,結(jié)果表明氣動(dòng)載荷對(duì)整體葉盤的固有頻率影響較小,變化率在1.2%以內(nèi);硬涂層對(duì)一階固有頻率影響較小,而使二、三階固有頻率有明顯的提高,并且會(huì)降低四階固有頻率,變化率在0.9%~2.3%內(nèi)。

    表4 氣動(dòng)載荷下整體葉盤的固有頻率(靜頻)

    考慮到整體葉盤在轉(zhuǎn)動(dòng)時(shí)受到離心載荷和氣動(dòng)載荷的影響,分別計(jì)算了整體葉盤在離心載荷、離心載荷和氣動(dòng)載荷共同作用下的固有頻率,見表5。對(duì)不同轉(zhuǎn)速下的整體葉盤固有頻率進(jìn)行比較,離心載荷會(huì)提高一、三階固有頻率,對(duì)二、四階固有頻率影響較小,變化率在1.1%以內(nèi)。在離心載荷和氣動(dòng)載荷共同作用下,與離心載荷作用下的固有頻率相比較,會(huì)降低整體葉盤的固有頻率,且四階固有頻率降低較明顯,變化率在2.6%以內(nèi)。

    表5 多載荷下整體葉盤的固有頻率(動(dòng)頻)

    分別計(jì)算了硬涂層整體葉盤在離心載荷、離心載荷和氣動(dòng)載荷共同作用下的固有頻率,見表6。對(duì)不同轉(zhuǎn)速下的硬涂層整體葉盤固有頻率進(jìn)行比較,離心載荷會(huì)提高一、三階固有頻率,對(duì)二、四階固有頻率影響較小,變化率在1.1%以內(nèi)。在離心載荷和氣動(dòng)載荷共同作用下,與離心載荷作用下的固有頻率相比較,會(huì)降低硬涂層整體葉盤的固有頻率,且四階固有頻率降低較明顯,變化率在2.7%以內(nèi)。

    表6 多載荷下硬涂層整體葉盤的固有頻率(動(dòng)頻)

    由以上分析可說明,在靜頻時(shí),氣動(dòng)載荷對(duì)整體葉盤固有頻率的影響較??;在動(dòng)頻時(shí),整體葉盤在高速旋轉(zhuǎn)狀態(tài)下,離心載荷會(huì)導(dǎo)致葉盤的伸張且增強(qiáng)了葉片的彎曲效應(yīng);而作用在葉片上氣動(dòng)載荷會(huì)減緩彎曲效應(yīng);綜合考慮2種載荷同時(shí)作用時(shí),可以降低葉片的彎曲頻率。

    3.2 振型分析

    在預(yù)應(yīng)力模態(tài)分析下得到72階模態(tài)振型圖,根據(jù)循環(huán)對(duì)稱性,得到葉片第一、二、三、四階模態(tài)振型圖。其中第37階和第55階以節(jié)圓振動(dòng)方式呈現(xiàn),其余階次以節(jié)徑振動(dòng)方式呈現(xiàn),文中取特征性振型圖進(jìn)行分析,如圖4所示。

    圖4 模態(tài)振型

    1—36階、41—54階、62—72階主要是葉片振動(dòng)占主導(dǎo)。其中1—18階主要是葉片的彎曲振動(dòng),19—36階主要是葉片的扭轉(zhuǎn)振動(dòng),其余階次主要是葉片和輪盤形成的耦合振動(dòng)。由圖4可以看出:葉盤的形變量從中心到葉端逐漸增加,在葉端達(dá)到最大。

    3.3 應(yīng)力分析

    在設(shè)計(jì)工況的情況下,對(duì)整體葉盤進(jìn)行了應(yīng)力分析,圖5所示是工況為30 000 r/min的整體葉盤在不同載荷下的等效應(yīng)力分布圖??芍涸诓煌d荷下整體葉盤的等效應(yīng)力變化趨勢(shì)大致相同,最小的等效應(yīng)力分布在輪盤和葉端區(qū)域,最大的等效應(yīng)力分布在靠近葉根葉片上,伴隨載荷數(shù)目的增加,最大等效應(yīng)力也隨之增加。由圖5還可知,與離心載荷相比,氣動(dòng)載荷對(duì)整體葉盤等效應(yīng)力的影響更大,但是整體葉盤的等效應(yīng)力分布與離心載荷有關(guān)。

    圖5 不同載荷下整體葉盤的等效應(yīng)力云圖

    對(duì)其他工況下的整體葉盤的等效應(yīng)力進(jìn)行研究,發(fā)現(xiàn)整體葉盤等效應(yīng)力分布大致相同,因此不再給出整體葉盤的等效應(yīng)力云圖。

    為了解不同載荷下硬涂層對(duì)整體葉盤等效應(yīng)力的影響,對(duì)不同載荷下的整體葉盤和硬涂層等效應(yīng)力進(jìn)行采集并分析。圖6所示為不同載荷下最大等效應(yīng)力的變化曲線。

    圖6 不同載荷下的最大等效應(yīng)力

    由圖6(a)可發(fā)現(xiàn):在離心載荷作用下硬涂層對(duì)整體葉盤的等效應(yīng)力的影響較小,可忽略不計(jì)。由圖6(b)可發(fā)現(xiàn):在離心載荷和氣動(dòng)載荷共同作用下硬涂層對(duì)整體葉盤的最大等效應(yīng)力的影響較大,變化范圍在11%以內(nèi);隨著轉(zhuǎn)速和氣動(dòng)載荷的增加,硬涂層對(duì)整體葉盤最大等效應(yīng)力增長(zhǎng)速率有減緩的趨勢(shì),但在實(shí)際運(yùn)轉(zhuǎn)中,還應(yīng)考慮其對(duì)溫度等因素的影響。

    3.4 共振分析

    整體葉盤在實(shí)際運(yùn)轉(zhuǎn)下要完全避免共振難以實(shí)現(xiàn),借助設(shè)計(jì)工況以無硬涂層整體葉盤和硬涂層整體葉盤模型為對(duì)象進(jìn)行分析研究。根據(jù)坎貝爾圖(Campbell),分析整體葉盤發(fā)生共振的條件,可以避免具有破壞性的振動(dòng)。根據(jù)表5和表6繪制坎貝爾圖,如圖7所示。圖中曲線代表轉(zhuǎn)速和固有頻率的關(guān)系,射線為激振頻率射線,兩者的交點(diǎn)為共振點(diǎn),共振點(diǎn)對(duì)應(yīng)的轉(zhuǎn)速為共振轉(zhuǎn)速。

    圖7 整體葉盤的坎貝爾圖

    對(duì)于整體葉盤,當(dāng)激振頻率與固有頻率滿足公式(14)時(shí),葉盤將發(fā)生共振現(xiàn)象。

    Fi=αFe

    (14)

    式中:Fi為葉盤的固有頻率;Fe為激振力頻率;當(dāng)α=1時(shí)成為共振。

    由圖7(a)可知,第三階固有頻率與K=7倍激振力的交點(diǎn)距30 000 r/min工況點(diǎn)比較接近,對(duì)應(yīng)轉(zhuǎn)速為30 458 r/min,共振裕度為1.5%。第四、二階固有頻率分別與K=9、K=5倍激振力的交點(diǎn)距32 400 r/min工況點(diǎn)比較接近,對(duì)應(yīng)轉(zhuǎn)速分別為32 148、31 652 r/min,共振裕度分別為0.8%、2.3%。第三階固有頻率與K=6倍激振力的交點(diǎn)距34 800 r/min工況點(diǎn)比較接近,對(duì)應(yīng)轉(zhuǎn)速為35 551 r/min,共振裕度為2.1%。第四階固有頻率與K=8倍激振力的交點(diǎn)距37 200 r/min工況點(diǎn)比較接近,對(duì)應(yīng)轉(zhuǎn)速為36 191 r/min,共振裕度為2.7%。第二、一階固有頻率分別與K=4、K=1倍激振力的交點(diǎn)距39 600 r/min工況點(diǎn)比較接近,對(duì)應(yīng)轉(zhuǎn)速分別為39 538、39 131 r/min,共振裕度分別為0.2%、1.1%。結(jié)果表明:整體葉盤在30 458、32 148、31 652、35 551、36 191、39 538、39 131 r/min轉(zhuǎn)速下較容易產(chǎn)生共振,為避免出現(xiàn)振動(dòng)損壞現(xiàn)象需要改變激振力和整體葉盤固有頻率。

    由圖7(b)可知,第三階固有頻率與K=7倍激振力的交點(diǎn)距30 000 r/min工況點(diǎn)比較接近,對(duì)應(yīng)轉(zhuǎn)速為30 554 r/min,共振裕度為1.8%。第四、二階固有頻率分別與K=9、K=5倍激振力的交點(diǎn)距32 400 r/min工況點(diǎn)比較接近,對(duì)應(yīng)轉(zhuǎn)速分別為32 079、32 453 r/min,共振裕度分別為1.0%、0.2%。第三階固有頻率與K=6倍激振力的交點(diǎn)距34 800 r/min工況點(diǎn)比較接近,對(duì)應(yīng)轉(zhuǎn)速為35 767 r/min,共振裕度為2.8%。第四階固有頻率與K=8倍激振力的交點(diǎn)距37 200 r/min工況點(diǎn)比較接近,對(duì)應(yīng)轉(zhuǎn)速為36 140 r/min,共振裕度為2.8%。第一階固有頻率與K=1倍激振力的交點(diǎn)距39 600 r/min工況點(diǎn)比較接近,對(duì)應(yīng)轉(zhuǎn)速為39 118 r/min,共振裕度為1.2%。結(jié)果表明:硬涂層整體葉盤比無硬涂層的整體葉盤共振點(diǎn)數(shù)由7個(gè)降為6個(gè)。

    4 總結(jié)

    基于流固耦合方法,在設(shè)計(jì)工況下對(duì)有無硬涂層的整體葉盤進(jìn)行葉盤表面氣動(dòng)載荷模擬,再利用預(yù)應(yīng)力模態(tài)分析方法對(duì)葉盤進(jìn)行靜力和模態(tài)分析,求解葉盤的表面氣動(dòng)載荷、固有頻率以及等效應(yīng)力的分布,繪制葉盤的坎貝爾圖,結(jié)果如下:

    (1)葉片壓力面涂敷硬涂層后對(duì)整體葉盤的固有頻率影響較小。離心載荷會(huì)提高整體葉盤的固有頻率,在離心載荷和氣動(dòng)載荷同時(shí)作用下會(huì)降低整體葉盤的固有頻率。

    (2)整體葉盤葉片在第一、二階分別以彎曲、扭轉(zhuǎn)振動(dòng)為主導(dǎo),第三、四階以葉片和輪盤的耦合振動(dòng)為主。整體葉盤前72階模態(tài)振型中,輪盤形變量最小,葉片頂端的形變量最大。

    (3)與離心載荷相比,氣動(dòng)載荷對(duì)整體葉盤等效應(yīng)力的影響更大。在離心載荷和氣動(dòng)載荷共同作用下硬涂層會(huì)降低整體葉盤的最大等效應(yīng)力,變化范圍在11%以內(nèi)。

    (4)比較分析葉片壓力面涂敷硬涂層整體葉盤與無硬涂層整體葉盤的動(dòng)頻特性,可以發(fā)現(xiàn)葉片壓力面涂敷硬涂層會(huì)使整體葉盤共振點(diǎn)數(shù)量從7個(gè)降為6個(gè),對(duì)減少共振現(xiàn)象有積極作用。

    猜你喜歡
    葉盤共振固有頻率
    現(xiàn)場(chǎng)測(cè)定大型水輪發(fā)電機(jī)組軸系的固有頻率
    安然 與時(shí)代同頻共振
    選硬人打硬仗——紫陽縣黨建與脫貧同頻共振
    CTA 中紡院+ 化纖聯(lián)盟 強(qiáng)強(qiáng)聯(lián)合 科技共振
    某型航空發(fā)動(dòng)機(jī)整體葉盤強(qiáng)度分析*
    總溫總壓測(cè)頭模態(tài)振型變化規(guī)律研究
    改革是決心和動(dòng)力的共振
    A novel functional electrical stimulation-control system for restoring motor function of post-stroke hemiplegic patients
    轉(zhuǎn)向系統(tǒng)固有頻率設(shè)計(jì)研究
    基于FE-ERSM航空發(fā)動(dòng)機(jī)葉盤結(jié)構(gòu)可靠性研究
    我的女老师完整版在线观看| 757午夜福利合集在线观看| 美女xxoo啪啪120秒动态图 | 夜夜躁狠狠躁天天躁| 国产黄a三级三级三级人| 91久久精品电影网| 欧美黑人欧美精品刺激| 一本综合久久免费| 麻豆成人av在线观看| 亚洲熟妇熟女久久| 变态另类成人亚洲欧美熟女| 青草久久国产| 亚洲自偷自拍三级| 免费av不卡在线播放| 国产亚洲欧美98| 亚洲av免费在线观看| 变态另类成人亚洲欧美熟女| 日韩欧美在线二视频| 亚洲黑人精品在线| 午夜影院日韩av| 在线观看舔阴道视频| 午夜激情欧美在线| 变态另类丝袜制服| 欧美不卡视频在线免费观看| 99热这里只有是精品在线观看 | 成人亚洲精品av一区二区| 欧美日本亚洲视频在线播放| 日日摸夜夜添夜夜添小说| 国产精华一区二区三区| 村上凉子中文字幕在线| 可以在线观看毛片的网站| 久久久国产成人免费| 久久精品影院6| 精品一区二区三区人妻视频| 黄色丝袜av网址大全| 91久久精品国产一区二区成人| 午夜福利高清视频| 特级一级黄色大片| 国内精品久久久久久久电影| 色哟哟·www| 97碰自拍视频| 91狼人影院| 久久久国产成人精品二区| 九九热线精品视视频播放| 亚洲国产精品999在线| 欧美+日韩+精品| 免费在线观看成人毛片| 久久欧美精品欧美久久欧美| 男人狂女人下面高潮的视频| av视频在线观看入口| 又粗又爽又猛毛片免费看| 中文字幕精品亚洲无线码一区| 91久久精品电影网| 国产精品电影一区二区三区| 精品欧美国产一区二区三| 床上黄色一级片| av福利片在线观看| 免费在线观看成人毛片| 最近视频中文字幕2019在线8| 午夜两性在线视频| 国产成人啪精品午夜网站| 亚洲七黄色美女视频| 国产私拍福利视频在线观看| 午夜福利在线观看免费完整高清在 | 一个人观看的视频www高清免费观看| 变态另类成人亚洲欧美熟女| 高清在线国产一区| 97热精品久久久久久| 欧美激情在线99| 国产白丝娇喘喷水9色精品| 91在线观看av| 国产精品精品国产色婷婷| 欧美区成人在线视频| 国产精品久久电影中文字幕| 最近中文字幕高清免费大全6 | 五月伊人婷婷丁香| 免费一级毛片在线播放高清视频| 男人的好看免费观看在线视频| 色哟哟·www| 亚洲乱码一区二区免费版| 日韩欧美在线乱码| 精品国产三级普通话版| 精品久久久久久久末码| 99久久久亚洲精品蜜臀av| 精华霜和精华液先用哪个| 尤物成人国产欧美一区二区三区| 99riav亚洲国产免费| 国产亚洲欧美在线一区二区| 色在线成人网| 国产老妇女一区| 精品久久久久久久人妻蜜臀av| 美女 人体艺术 gogo| 舔av片在线| 免费电影在线观看免费观看| 国产三级中文精品| 国产精品乱码一区二三区的特点| 夜夜躁狠狠躁天天躁| 亚洲一区高清亚洲精品| 国产黄a三级三级三级人| 亚洲美女视频黄频| 老熟妇仑乱视频hdxx| 日韩国内少妇激情av| 尤物成人国产欧美一区二区三区| 久久久久性生活片| 在线观看午夜福利视频| 综合色av麻豆| 久久精品国产清高在天天线| 日本免费一区二区三区高清不卡| 欧美高清成人免费视频www| 中文字幕免费在线视频6| 美女cb高潮喷水在线观看| 深夜精品福利| 国产精品影院久久| 少妇高潮的动态图| 校园春色视频在线观看| 成年人黄色毛片网站| 噜噜噜噜噜久久久久久91| 波多野结衣高清无吗| 亚洲最大成人中文| 亚洲不卡免费看| 在线a可以看的网站| 啦啦啦韩国在线观看视频| 国产精品亚洲av一区麻豆| 黄片小视频在线播放| 亚洲国产高清在线一区二区三| 天堂网av新在线| 日本黄色片子视频| 久久久久久久精品吃奶| 久久久久久国产a免费观看| 免费av毛片视频| 日本与韩国留学比较| 国产精品久久久久久久电影| 婷婷亚洲欧美| 天堂av国产一区二区熟女人妻| 欧美3d第一页| 女人被狂操c到高潮| 国产精品av视频在线免费观看| 欧美区成人在线视频| 日日摸夜夜添夜夜添小说| www.熟女人妻精品国产| 简卡轻食公司| 日日干狠狠操夜夜爽| 成人特级av手机在线观看| 男人舔女人下体高潮全视频| 亚洲一区二区三区色噜噜| 在线十欧美十亚洲十日本专区| 乱码一卡2卡4卡精品| 国产69精品久久久久777片| 亚洲欧美清纯卡通| 岛国在线免费视频观看| 国产成人福利小说| 亚洲国产精品999在线| 欧美又色又爽又黄视频| 日韩中文字幕欧美一区二区| 国产精品人妻久久久久久| 极品教师在线视频| 亚洲最大成人手机在线| 美女xxoo啪啪120秒动态图 | www.色视频.com| 好男人电影高清在线观看| 麻豆成人av在线观看| 成人一区二区视频在线观看| 亚洲成人久久爱视频| 久久久久久久午夜电影| 亚洲片人在线观看| 91字幕亚洲| 国产黄片美女视频| 国产午夜精品论理片| 亚洲欧美日韩高清在线视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲国产精品成人综合色| 久久精品综合一区二区三区| 91麻豆精品激情在线观看国产| 国内精品一区二区在线观看| 丁香六月欧美| 高清在线国产一区| 欧美另类亚洲清纯唯美| 高清毛片免费观看视频网站| 69人妻影院| 亚洲内射少妇av| 国内少妇人妻偷人精品xxx网站| 99久久99久久久精品蜜桃| 尤物成人国产欧美一区二区三区| 国内毛片毛片毛片毛片毛片| 少妇丰满av| 国产精品乱码一区二三区的特点| 色综合欧美亚洲国产小说| 嫩草影院精品99| 99久久精品热视频| 日韩免费av在线播放| 中文字幕av在线有码专区| 亚洲av第一区精品v没综合| 久久这里只有精品中国| 国产一区二区三区视频了| 人人妻人人看人人澡| 国产又黄又爽又无遮挡在线| 老司机福利观看| x7x7x7水蜜桃| 午夜精品在线福利| 精品一区二区三区人妻视频| 给我免费播放毛片高清在线观看| 精品免费久久久久久久清纯| 色噜噜av男人的天堂激情| 久久精品人妻少妇| 女同久久另类99精品国产91| 熟妇人妻久久中文字幕3abv| 精品福利观看| 怎么达到女性高潮| 一级av片app| 别揉我奶头~嗯~啊~动态视频| 搡老岳熟女国产| 在线天堂最新版资源| 成人高潮视频无遮挡免费网站| 日日夜夜操网爽| 国产精品伦人一区二区| 最后的刺客免费高清国语| 欧美一区二区国产精品久久精品| 高清日韩中文字幕在线| 久久久久久久久大av| 一个人看视频在线观看www免费| 精品福利观看| 国产亚洲精品综合一区在线观看| 级片在线观看| 三级毛片av免费| 1024手机看黄色片| 亚洲欧美清纯卡通| 亚洲欧美日韩卡通动漫| 少妇裸体淫交视频免费看高清| 国产精品免费一区二区三区在线| 欧美在线黄色| 1000部很黄的大片| 一区二区三区激情视频| 午夜日韩欧美国产| 亚洲精品久久国产高清桃花| 一级毛片久久久久久久久女| 国产精品爽爽va在线观看网站| 亚洲av免费在线观看| 两个人的视频大全免费| 亚洲国产高清在线一区二区三| 久久国产精品人妻蜜桃| 亚洲av电影不卡..在线观看| 欧美3d第一页| 午夜激情福利司机影院| 别揉我奶头~嗯~啊~动态视频| 亚洲精品久久国产高清桃花| 精品久久久久久久末码| 淫秽高清视频在线观看| av福利片在线观看| 国产白丝娇喘喷水9色精品| .国产精品久久| 激情在线观看视频在线高清| 在线十欧美十亚洲十日本专区| 91狼人影院| 精品国内亚洲2022精品成人| 我要搜黄色片| 乱人视频在线观看| 乱人视频在线观看| 中文字幕av成人在线电影| 国产老妇女一区| 全区人妻精品视频| h日本视频在线播放| 国产精品99久久久久久久久| 中文字幕精品亚洲无线码一区| 高清毛片免费观看视频网站| 波多野结衣高清无吗| 国产成+人综合+亚洲专区| 免费在线观看成人毛片| 亚洲成人久久性| 亚洲三级黄色毛片| 国产成人a区在线观看| 亚洲美女搞黄在线观看 | 美女xxoo啪啪120秒动态图 | 国产精品久久久久久精品电影| 日韩中文字幕欧美一区二区| 两人在一起打扑克的视频| 午夜福利在线观看吧| 午夜精品久久久久久毛片777| 日韩 亚洲 欧美在线| 热99re8久久精品国产| 国产人妻一区二区三区在| 久久99热6这里只有精品| 免费高清视频大片| 国产精品av视频在线免费观看| 久久这里只有精品中国| 欧美日韩综合久久久久久 | 国产精品一区二区三区四区免费观看 | 性色avwww在线观看| 欧美zozozo另类| 色播亚洲综合网| 亚洲最大成人手机在线| 18禁黄网站禁片午夜丰满| 毛片女人毛片| 婷婷六月久久综合丁香| 国产高清视频在线播放一区| 97人妻精品一区二区三区麻豆| 国产男靠女视频免费网站| eeuss影院久久| 男女床上黄色一级片免费看| 精品人妻偷拍中文字幕| 婷婷精品国产亚洲av在线| 国产探花极品一区二区| 久久久精品大字幕| 午夜激情福利司机影院| 国产高潮美女av| 无遮挡黄片免费观看| 在线观看一区二区三区| www.色视频.com| 麻豆国产97在线/欧美| 制服丝袜大香蕉在线| 欧美激情在线99| 免费观看的影片在线观看| 国产成人av教育| 少妇人妻精品综合一区二区 | 久久精品国产99精品国产亚洲性色| 国产人妻一区二区三区在| 精品无人区乱码1区二区| 国产成人av教育| 高清毛片免费观看视频网站| 真人一进一出gif抽搐免费| 桃红色精品国产亚洲av| 久久天躁狠狠躁夜夜2o2o| 在线观看午夜福利视频| 有码 亚洲区| 日韩中文字幕欧美一区二区| 免费看日本二区| 一级作爱视频免费观看| 国产中年淑女户外野战色| 禁无遮挡网站| 亚洲成人久久爱视频| 成人性生交大片免费视频hd| 桃色一区二区三区在线观看| 亚洲国产色片| 在现免费观看毛片| 国产精品嫩草影院av在线观看 | 国产精品人妻久久久久久| 午夜精品一区二区三区免费看| 亚洲成av人片免费观看| 男女床上黄色一级片免费看| 少妇丰满av| 久久国产乱子伦精品免费另类| 999久久久精品免费观看国产| 黄色配什么色好看| 看十八女毛片水多多多| 国产精品国产高清国产av| 看黄色毛片网站| 在线观看舔阴道视频| 又黄又爽又免费观看的视频| 欧美日韩黄片免| 18禁黄网站禁片免费观看直播| 午夜精品一区二区三区免费看| 51国产日韩欧美| 又爽又黄无遮挡网站| 在线播放国产精品三级| 久久久久久久久中文| 欧美乱色亚洲激情| 久久国产精品人妻蜜桃| 看黄色毛片网站| 国产成人欧美在线观看| 成人毛片a级毛片在线播放| 国产亚洲欧美98| 日本三级黄在线观看| 免费人成在线观看视频色| 国产伦精品一区二区三区四那| 精品人妻一区二区三区麻豆 | АⅤ资源中文在线天堂| 精品一区二区免费观看| 色噜噜av男人的天堂激情| 亚洲18禁久久av| 岛国在线免费视频观看| 日韩有码中文字幕| 午夜老司机福利剧场| 亚洲美女黄片视频| 成人国产综合亚洲| 欧美日韩瑟瑟在线播放| 在线观看免费视频日本深夜| av福利片在线观看| 日日摸夜夜添夜夜添av毛片 | 国产一区二区三区在线臀色熟女| 久久6这里有精品| 国产精品永久免费网站| 99久久精品一区二区三区| 又紧又爽又黄一区二区| 婷婷精品国产亚洲av在线| 欧美性猛交黑人性爽| 日韩欧美精品免费久久 | 久久久久国产精品人妻aⅴ院| 全区人妻精品视频| 深夜a级毛片| 成人特级av手机在线观看| 中文字幕熟女人妻在线| 夜夜爽天天搞| 午夜亚洲福利在线播放| 亚洲人成电影免费在线| 一本精品99久久精品77| 人妻丰满熟妇av一区二区三区| 国产亚洲精品av在线| 在线观看免费视频日本深夜| 免费看a级黄色片| 久久久久久久午夜电影| 给我免费播放毛片高清在线观看| 如何舔出高潮| 极品教师在线视频| 白带黄色成豆腐渣| 日韩欧美精品v在线| 九色成人免费人妻av| 观看美女的网站| 麻豆国产av国片精品| 久久久久精品国产欧美久久久| 欧美激情国产日韩精品一区| 一区二区三区激情视频| 亚洲美女视频黄频| 久久久久久久精品吃奶| 99热这里只有精品一区| 老司机深夜福利视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 我的女老师完整版在线观看| 亚洲欧美日韩卡通动漫| 免费在线观看亚洲国产| 免费人成在线观看视频色| 大型黄色视频在线免费观看| 亚洲熟妇中文字幕五十中出| 国产麻豆成人av免费视频| 亚洲第一电影网av| 内射极品少妇av片p| 热99在线观看视频| 免费搜索国产男女视频| 欧美日韩亚洲国产一区二区在线观看| 欧美绝顶高潮抽搐喷水| 久久热精品热| 亚洲av.av天堂| 精品免费久久久久久久清纯| 老司机午夜十八禁免费视频| 国产激情偷乱视频一区二区| 国产精品久久久久久久电影| 1000部很黄的大片| 久久这里只有精品中国| 性色avwww在线观看| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 俺也久久电影网| 国内精品久久久久久久电影| 亚洲在线观看片| 亚洲精品一卡2卡三卡4卡5卡| av在线观看视频网站免费| 久久精品综合一区二区三区| 很黄的视频免费| av专区在线播放| 国内少妇人妻偷人精品xxx网站| 国产一区二区三区在线臀色熟女| 国产欧美日韩精品亚洲av| 亚洲自偷自拍三级| 日本 欧美在线| 欧美国产日韩亚洲一区| 婷婷亚洲欧美| 婷婷色综合大香蕉| 99久久无色码亚洲精品果冻| 老司机福利观看| 老熟妇乱子伦视频在线观看| 亚洲国产欧洲综合997久久,| 最近最新免费中文字幕在线| 日本三级黄在线观看| 好看av亚洲va欧美ⅴa在| 成年版毛片免费区| 99国产极品粉嫩在线观看| 欧美一区二区国产精品久久精品| 一个人免费在线观看电影| www.999成人在线观看| 男人舔女人下体高潮全视频| 三级国产精品欧美在线观看| 日日摸夜夜添夜夜添小说| 波野结衣二区三区在线| 亚洲真实伦在线观看| 神马国产精品三级电影在线观看| 国产一区二区亚洲精品在线观看| 美女cb高潮喷水在线观看| 很黄的视频免费| 久久久久久久久中文| 神马国产精品三级电影在线观看| 美女高潮喷水抽搐中文字幕| 亚洲欧美日韩卡通动漫| 久久草成人影院| 成人一区二区视频在线观看| 亚洲成人久久性| 国产一区二区在线av高清观看| 1024手机看黄色片| 窝窝影院91人妻| 免费在线观看影片大全网站| 精品一区二区免费观看| 亚洲av熟女| 亚洲av电影不卡..在线观看| 国产精品98久久久久久宅男小说| 精华霜和精华液先用哪个| 欧美一区二区国产精品久久精品| 99久久无色码亚洲精品果冻| 内地一区二区视频在线| 看十八女毛片水多多多| 18禁黄网站禁片免费观看直播| 国产精品女同一区二区软件 | 18禁黄网站禁片免费观看直播| 国产av不卡久久| 好男人在线观看高清免费视频| 日韩欧美在线乱码| 宅男免费午夜| 给我免费播放毛片高清在线观看| 男女那种视频在线观看| ponron亚洲| 一个人免费在线观看的高清视频| 精品人妻熟女av久视频| 久久久久免费精品人妻一区二区| 老熟妇仑乱视频hdxx| 精华霜和精华液先用哪个| 日韩高清综合在线| 欧美最黄视频在线播放免费| 91字幕亚洲| 亚洲av第一区精品v没综合| 桃色一区二区三区在线观看| 搞女人的毛片| 亚洲精品亚洲一区二区| 国产精品爽爽va在线观看网站| 国产精品美女特级片免费视频播放器| 尤物成人国产欧美一区二区三区| 亚洲精品在线美女| 97超视频在线观看视频| 好男人在线观看高清免费视频| 国内精品久久久久精免费| 国产精品国产高清国产av| 久久精品夜夜夜夜夜久久蜜豆| 天美传媒精品一区二区| 欧美丝袜亚洲另类 | 宅男免费午夜| 禁无遮挡网站| 精品久久久久久久人妻蜜臀av| 亚洲精品在线观看二区| 免费看美女性在线毛片视频| 精品日产1卡2卡| 黄色一级大片看看| 精品久久国产蜜桃| 国产日本99.免费观看| 国产极品精品免费视频能看的| 亚洲无线在线观看| 男女那种视频在线观看| 国内精品美女久久久久久| 色综合站精品国产| 久久草成人影院| 精品国产亚洲在线| 国产亚洲欧美在线一区二区| 禁无遮挡网站| 国产精品国产高清国产av| 免费高清视频大片| 久久伊人香网站| 国产乱人视频| 偷拍熟女少妇极品色| 精品久久久久久久久久久久久| 欧美日韩综合久久久久久 | 男人和女人高潮做爰伦理| 日韩国内少妇激情av| 免费电影在线观看免费观看| 18禁黄网站禁片午夜丰满| 人妻制服诱惑在线中文字幕| 亚州av有码| 日韩有码中文字幕| 身体一侧抽搐| 精品免费久久久久久久清纯| 亚洲人成网站高清观看| 波多野结衣高清无吗| 日韩高清综合在线| 国产av麻豆久久久久久久| 3wmmmm亚洲av在线观看| 久久久久九九精品影院| 精品不卡国产一区二区三区| АⅤ资源中文在线天堂| 男人狂女人下面高潮的视频| 亚洲中文字幕一区二区三区有码在线看| 亚洲av日韩精品久久久久久密| 国产成人啪精品午夜网站| 久久精品人妻少妇| 欧美成人a在线观看| 成人一区二区视频在线观看| 高清日韩中文字幕在线| 色视频www国产| 亚洲 欧美 日韩 在线 免费| 超碰av人人做人人爽久久| 国产 一区 欧美 日韩| 午夜福利高清视频| 91在线精品国自产拍蜜月| 久久人妻av系列| 一级黄色大片毛片| 亚洲av电影在线进入| 中文字幕人妻熟人妻熟丝袜美| 久久久久亚洲av毛片大全| 天天躁日日操中文字幕| 99久久精品一区二区三区| 草草在线视频免费看| 一级作爱视频免费观看| 亚洲av二区三区四区| 一个人免费在线观看的高清视频| 性欧美人与动物交配| av天堂中文字幕网| 淫秽高清视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 成人国产综合亚洲| 男女国产视频网站| videos熟女内射| 国产真实伦视频高清在线观看| 国模一区二区三区四区视频| 小蜜桃在线观看免费完整版高清| 99热国产这里只有精品6| 久久97久久精品| 日韩视频在线欧美| 又大又黄又爽视频免费| 插阴视频在线观看视频| 特级一级黄色大片| av又黄又爽大尺度在线免费看| 国产v大片淫在线免费观看| 91午夜精品亚洲一区二区三区| 一个人观看的视频www高清免费观看|