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

    直板葉片振動響應(yīng)試驗的相似關(guān)系確定方法

    2019-12-13 08:33:52張文笛于清文
    航空發(fā)動機(jī) 2019年5期
    關(guān)鍵詞:原型撓度載荷

    張文笛,羅 忠,于清文,王 菲

    (東北大學(xué)機(jī)械工程與自動化學(xué)院,沈陽110819)

    0 引言

    航空發(fā)動機(jī)壓氣機(jī)工作葉片在設(shè)計中可簡化為懸臂薄板,基于薄板類模型,通常進(jìn)行的振動特性分析主要包括系統(tǒng)固有特性(固有頻率、振型)分析以及動力響應(yīng)分析。其中動力響應(yīng)為結(jié)構(gòu)在特定或隨機(jī)載荷作用下?lián)隙入S時間變化的時程曲線。葉片通常會受到離心力、氣動力、熱載荷、振動載荷等不均勻脈沖載荷的作用,因此進(jìn)行葉片在脈沖載荷下的振動響應(yīng)分析具有重要意義。由于航空發(fā)動機(jī)壓氣機(jī)工作葉片具有復(fù)雜的工作環(huán)境、加工復(fù)雜性以及成本、場地等的局限性,常采用相似理論中的模型對原型固有特性進(jìn)行預(yù)測,而對于振動響應(yīng)與載荷的相似關(guān)系的研究較少。

    在相似模型動力響應(yīng)研究中,Quercetti等[1]以燃料桶下落為對象,以沖擊載荷下動力響應(yīng)為主要參數(shù)驗證了原型與模型的相似關(guān)系;Ramu等[2]采用量綱分析法確定結(jié)構(gòu)模型參數(shù)的相似關(guān)系,并采用實(shí)例驗證懸臂梁分別受橫向和縱向階躍載荷的撓度、應(yīng)力等相似關(guān)系;Qian等[3-4]對層合板進(jìn)行沖擊試驗,分析了彈性變形階段與塑性變形階段完全相似關(guān)系的預(yù)測情況,得到彈性變形階段可以使用完全相似關(guān)系,塑性變形階段不能使用完全相似關(guān)系的結(jié)論;陳喆等[5]在相似理論和量綱分析的基礎(chǔ)上結(jié)合模態(tài)分析理論推導(dǎo)了結(jié)構(gòu)固有特性、頻域響應(yīng)與時域響應(yīng)的完全相似關(guān)系,并對機(jī)翼進(jìn)行有限元仿真,得到加速度與位移的時程曲線;Wu[6-7]采用完全相似模型預(yù)測彈性支撐板在圓環(huán)動載荷作用下的彈性振動特性,根據(jù)動力學(xué)方程,采用方程分析法確定各物理量之間的相似關(guān)系。針對不完全相似的研究,羅忠等[8]針對板、殼等模型進(jìn)行固有特性的不完全幾何相似關(guān)系的確定;James等[9-10]將渦輪發(fā)動機(jī)葉片簡化為簡支矩形板,通過方程分析法預(yù)測其完全幾何相似與不完全幾何相似下的固有特性;又將其簡化為簡支薄壁復(fù)合工字梁,通過控制方程推導(dǎo)出靜態(tài)撓度的相似律,預(yù)測原型與模型最大靜撓度的相似性與準(zhǔn)確性;Simites等[11]對線性載荷下簡支板相似關(guān)系進(jìn)行分析,通過對不同方程下的預(yù)測撓度進(jìn)行比較,得到預(yù)測誤差最小的相似關(guān)系;張瑋等[12]根據(jù)有限元預(yù)測系數(shù)法對單參數(shù)畸變相似模型進(jìn)行沖擊響應(yīng)研究,實(shí)現(xiàn)預(yù)測原模型的振動特性。而目前關(guān)于直板葉片在脈沖載荷下振動響應(yīng)相似關(guān)系的研究較少。

    本文針對直板葉片的振動響應(yīng)相似模型設(shè)計問題,結(jié)合相似理論中的方程分析法和敏感性分析法,建立葉片試驗?zāi)P偷恼駝禹憫?yīng)分析的相似關(guān)系。分析振動響應(yīng)機(jī)理并通過有限元預(yù)測系數(shù)法,提高預(yù)測精度,消除尺寸效應(yīng)帶來的影響。通過數(shù)值仿真結(jié)果對振動響應(yīng)進(jìn)行驗證,得到了理想的預(yù)測效果。

    1 動力學(xué)模型的確定

    將航空發(fā)動機(jī)壓氣機(jī)工作葉片簡化為懸臂薄板結(jié)構(gòu),如圖1所示。左端為固定端,右端為自由端,以固定端前端點(diǎn)為原點(diǎn)建立直角坐標(biāo)系,定義懸臂薄板結(jié)構(gòu)長、寬、厚分別為l、b、h,薄板可受集中載荷F或均布載荷q。

    圖1 懸臂薄板結(jié)構(gòu)

    根據(jù)板殼振動理論[13]可得懸臂薄板的振動微分方程

    當(dāng)載荷作用時間非常小時,視其為沖擊情況下的振型分量

    2 脈沖載荷下振動響應(yīng)的相似關(guān)系

    在評價脈沖響應(yīng)特性時,由于頻響特性只提供響應(yīng)幅度隨頻率的分布,而結(jié)構(gòu)的位移、速度、加速度以及彈性力學(xué)參數(shù)可作為表征線彈性振動問題中的主要指標(biāo),因此時域響應(yīng)特性也尤為重要。且由于碰撞、爆炸、地震等工況產(chǎn)生的脈沖載荷加載時間十分短暫,振動響應(yīng)只受到相對少數(shù)頻率支配,因此本文采用時域響應(yīng)作為主要指標(biāo)進(jìn)行相似模型的預(yù)測。

    2.1 完全幾何相似下振動響應(yīng)相似關(guān)系

    將式(1)中的Laplace算子展開,并運(yùn)用積分類比法可得

    下標(biāo)p代表原型,m代表模型,根據(jù)相似理論中的方程分析法,原型與模型的表達(dá)式分別為

    對原型的幾何尺寸進(jìn)行縮放,假定長度、寬度、厚度尺寸縮放比例相同,即為完全相似條件

    式中:λ 為幾何相似比;λl、λb、λh分別為葉片長度、寬度、厚度相似比;λx、λy為葉片上任一點(diǎn)橫、縱坐標(biāo)的相似比。

    對抗彎剛度方程進(jìn)行分析可得抗彎剛度的相似比

    式中:λE為葉片彈性模量的相似比。

    對式(4)、(5)進(jìn)行比較可得

    式中:λw為葉片撓度相似比;λρ為密度相似比;λω0為靜固有頻率相似比;λq為均布載荷相似比;λt為時間相似比。

    根據(jù)相似原理,式(8)相似關(guān)系為

    從而可得速度與加速度的相似關(guān)系λv、λa分別為

    對于集中載荷與均布載荷,有

    式中:λF為集中載荷相似比;λq為均布載荷相似比。

    2.2 不完全幾何相似下振動響應(yīng)相似關(guān)系

    壓氣機(jī)工作葉片在完全幾何相似下,可以很好地預(yù)測動力學(xué)特性,但在實(shí)際工程及試驗中,由于葉片的幾何特性,厚度一般小于長度和寬度1個量級以上,縮放后會產(chǎn)生生產(chǎn)制造等諸多問題,此時需要改變模型的局部幾何尺寸縮放比例,這種模型可稱為畸變模型或幾何不完全相似模型。畸變主要為對葉片厚度、長度或?qū)挾鹊幕儭?/p>

    2.2.1 對葉片的厚度進(jìn)行畸變

    假定 λl=λb=λ≠λh,則式(8)為

    從而得到相關(guān)參數(shù)的相似比

    由此可知,葉片的厚度相似比并不與長度、寬度的相似比發(fā)生沖突,λh可作為獨(dú)立性條件的相似因子。

    2.2.2 對葉片的長度或?qū)挾冗M(jìn)行畸變

    與對厚度進(jìn)行畸變的結(jié)果不同,當(dāng)λl≠λb時,通過代入到振動微分方程中得到

    通過上述的完全幾何相似關(guān)系的推導(dǎo),得到關(guān)于固有頻率的待定相似關(guān)系的統(tǒng)一格式表達(dá)式為

    式中:α為長度相似比加權(quán)冪數(shù);β為寬度相似比加權(quán)冪數(shù)。

    從而可知時間的相似比統(tǒng)一表達(dá)式為

    位移的待定相似比關(guān)系為統(tǒng)一的表達(dá)格式為

    當(dāng)載荷為集中力時,聯(lián)立式(12)、(18)可寫為

    根據(jù)文獻(xiàn)[14]可知,對長度和寬度的敏感性分析建立精確相似關(guān)系,即可得到α與β的值。葉片的長、寬對葉片不同階固有頻率有不同影響,不同階固有頻率對應(yīng)的相似關(guān)系也不同,對于時間的相似關(guān)系很難確定。通過對2種載荷下振動響應(yīng)分析可知,載荷強(qiáng)度相同時,產(chǎn)生最大幅值的時間是固定的,通過對不同相似比模型仿真得到各相似比下的離散點(diǎn)進(jìn)行擬合,再對擬合函數(shù)求導(dǎo),即可得到在小范圍相似比內(nèi)長度對撓度產(chǎn)生的敏感率。由于敏感性分析只針對結(jié)構(gòu)尺寸發(fā)生小范圍變化的情況,因此對于大范圍變化模型進(jìn)行相似設(shè)計時,可以采用校正模型,原型、畸變模型與校正模型的關(guān)系如圖2所示。

    圖2 不完全幾何相似過程

    在實(shí)際工程中,往往無法針對原型載荷匹配完全滿足時間與強(qiáng)度幅值要求的模型載荷,會出現(xiàn)原型載荷作用時間較長,通過縮比后模型載荷加載時間無法滿足模型設(shè)計要求等問題。由式(14)、(17)可知,時間的相似比由材料參數(shù)以及幾何參數(shù)的相似獲得,針對厚度較小的結(jié)構(gòu),在對長度、寬度縮放的同時可以對厚度進(jìn)行縮小,得到更大的時間相似比,即可達(dá)到滿足模型試驗的設(shè)計要求。

    3 脈沖載荷下振動響應(yīng)預(yù)測方法驗證

    3.1 完全幾何相似下振動響應(yīng)預(yù)測

    以Ti-6Al-4V鈦合金作為原型材料,40Cr鋼作為模型材料參數(shù)為例,原型與模型的幾何參數(shù)、材料參數(shù)分別見表1、2,采用ANSYS有限元仿真軟件對葉片進(jìn)行動力響應(yīng)做相似比的驗證。

    表1 原型與模型幾何參數(shù)

    表2 原型與模型材料參數(shù)

    為了分析相似比λw、λq與時間的關(guān)系,采用2種不同類型的脈沖載荷,1種為矩形載荷,載荷強(qiáng)度與時間無關(guān);另1種為三角形載荷,載荷強(qiáng)度與時間呈線性關(guān)系,最大載荷強(qiáng)度與矩形載荷強(qiáng)度相同,根據(jù)給出的相似關(guān)系和材料幾何參數(shù)、物理參數(shù),得到原型與模型的載荷形式,如圖3所示。

    圖3 不同形式載荷

    采用表1中原型幾何尺寸在原型葉片上施加100 Pa均布載荷,通過有限元仿真得到直板葉片不考慮阻尼的瞬態(tài)振動響應(yīng),拾取自由端中點(diǎn)的位移曲線。2種不同形式載荷在不同加載時間下的時程曲線對比如圖4所示。通過對彈性力學(xué)及振動理論的分析可知,載荷峰值的強(qiáng)度影響振動響應(yīng)的幅值,但不影響最大峰值產(chǎn)生的時間與軌跡;當(dāng)載荷加載時間小于葉片產(chǎn)生最大幅值的時間時,振動響應(yīng)近似于諧波函數(shù);對于矩形載荷,當(dāng)加載時間大于葉片產(chǎn)生最大幅值的時間時,由于相同時間段產(chǎn)生沖量相同,不同加載時間的振動響應(yīng)在到達(dá)位移峰值前的軌跡相同,但加載時間越長振動響應(yīng)的自由振動幅值越小,這是由于載荷方向與葉片回彈方向相反,阻礙葉片運(yùn)動并且此刻載荷對葉片作負(fù)功消耗葉片能量,因此葉片在自由振動時的幅值減??;當(dāng)載荷加載時間過長時可視為階躍載荷的瞬態(tài)響應(yīng)部分(本文不予考慮)。對于三角形載荷,產(chǎn)生沖量速度與加載時間有關(guān),由圖3(b)分析可知,在相同峰值情況下,加載時間越長,產(chǎn)生沖量速度越小,到達(dá)葉片振動峰值時間越長,由于葉片隨時間向下彎曲,時間越長彎曲產(chǎn)生的抵抗載荷力越大,因此加載時間加長可能無法到達(dá)該載荷峰值下的最大位移,符合圖4(b)中的軌跡變化規(guī)律。由于不同加載時間產(chǎn)生的時程曲線軌跡規(guī)律不同,無法單獨(dú)使用有限元預(yù)測系數(shù)法,因此必須通過相似比得到加載時間的近似值作為導(dǎo)向,再運(yùn)用有限元預(yù)測系數(shù)法加以修正來保證模型的預(yù)測精度。

    圖4 不同加載時間下的時程曲線

    根據(jù)得到的理論相似比,原型葉片施加不同加載時間時100 Pa、-2 N的2種載荷下的時程曲線預(yù)測結(jié)果如圖5、6所示。其中模型載荷的相似比為0.1。從圖中可見,由于幾何尺寸的縮小導(dǎo)致尺寸效應(yīng)的產(chǎn)生,通過選取原型與預(yù)測后時程曲線上幅值最大點(diǎn)或遠(yuǎn)端相位差較大的自由振動幅值來得到時間比值的修正系數(shù),將修正系數(shù)代入到模型的加載時間中,最后得到修正后的時程曲線,保證不同階段振動軌跡皆可得到很好的預(yù)測效果。

    圖5 矩形載荷預(yù)測修正曲線

    圖6 三角形載荷預(yù)測修正曲線

    3.2 不完全幾何相似驗證

    葉片的參數(shù)和由ANSYS仿真得到在施加均布載荷為50 Pa的加載時間及最大撓度分別見表3、4,得到擬合曲線如圖7~10所示。

    表3 長度變化對撓度產(chǎn)生的影響

    表4 寬度變化對撓度產(chǎn)生的影響

    圖7 λl變化對時間擬合曲線

    圖8 λb變化對時間擬合曲線

    圖9 λl變化對位移擬合曲線

    圖10 λb變化對位移擬合曲線

    得到多項式擬合曲線方程

    計算得到其敏感性

    得到多項式擬合曲線方程

    計算得到其敏感性

    根據(jù)精確相似關(guān)系計算方法

    解得

    因此得到葉片加載時間的相似關(guān)系

    得到多項式擬合曲線方程

    計算得到其敏感性

    得到多項式擬合曲線方程

    計算得到其敏感性

    解得

    因此得到葉片位移的相似關(guān)系

    聯(lián)立式(25)與(30)即可得到關(guān)于葉片振動響應(yīng)中速度與加速的相似比。同樣可求得葉片前2階固有頻率的相似比

    本文只針對不完全幾何相似有阻尼情況給出位移與加速度的時程預(yù)測圖,因為完全幾何相似下不考慮阻尼可以很好地看出預(yù)測曲線與原型曲線的相位差。當(dāng)考慮不同材料模型預(yù)測原型的動力響應(yīng)時,必須考慮阻尼帶來的影響[15-16],采用瑞利阻尼,通過任意2階固有頻率及其阻尼比即可得到質(zhì)量阻尼與剛度阻尼。由于阻尼比為無量綱數(shù),可假定理想模型的各階阻尼比與原型的各階阻尼比相同,將原型參數(shù)帶入式(31)得到校正模型的相對應(yīng)階次的固有頻率,再代入到式(32)中即可求出校正模型的質(zhì)量與剛度阻尼。

    式中:i、j為任意階數(shù);ω 為固有頻率;ξ、Ψ 分別為質(zhì)量阻尼與剛度阻尼;η為阻尼比。

    假定前2階固有頻率阻尼比η1=0.17,η2=0.14[17],以 λl=0.9,λb=1.1為例,施加產(chǎn)生最大位移時間0.0118 s和大于產(chǎn)生最大位移時間0.025 s的2種矩形載荷加載,分別對以上相似關(guān)系進(jìn)行驗證,結(jié)果如圖11所示。從圖中可見,模型的位移、加速度曲線可以很好地預(yù)測原型曲線,并且無明顯相位差,無漏峰、錯峰現(xiàn)象,數(shù)值仿真選取時間增量的大小是導(dǎo)致峰值預(yù)測結(jié)果出現(xiàn)偏差的原因之一,同時可以運(yùn)用在時間相似比較大的工作情況下。針對完全幾何相似情況,方程分析法結(jié)合有限元預(yù)測系數(shù)法提高了模型預(yù)測精確度;針對不完全幾何相似情況,方程分析法結(jié)合敏感性分析法可以準(zhǔn)確預(yù)測出原型的振動響應(yīng)。

    圖11 考慮阻尼情況不完全幾何相似的預(yù)測曲線

    4 結(jié)論

    (1)針對壓氣機(jī)工作葉片的振動響應(yīng)相似模型設(shè)計問題,基于板殼振動理論,結(jié)合方程分析法與有限元預(yù)測系數(shù)法,完善了葉片振動響應(yīng)問題的動力學(xué)相似關(guān)系,放寬了振動響應(yīng)中對載荷相似比的約束條件。

    (2)通過對不同模型施加不同時間的載荷,分析了振動響應(yīng)軌跡的變化機(jī)理,得到具有明顯的特征點(diǎn)作為幾何相似關(guān)系求解的參數(shù),可以作為敏感性分析的重要指標(biāo)。

    (3)針對幾何不完全相似,解決了大時間相似比的工作情況下振動響應(yīng)模型設(shè)計問題;采用數(shù)值仿真的方法對動力學(xué)相似進(jìn)行了驗證,對不同形式載荷有無阻尼的彈性時程位移曲線進(jìn)行了對比,皆得到了很好的預(yù)測效果。

    (4)通過完全幾何相似與不完全幾何相似關(guān)系的對比,可以得到結(jié)合敏感性分析法與方程分析法可以有效地解決尺寸效應(yīng)問題的結(jié)論。

    (5)線彈性振動響應(yīng)相似關(guān)系方法的確定,為得到?jīng)_擊載荷下應(yīng)力波對振動響應(yīng)的影響規(guī)律以及相似模型設(shè)計奠定基礎(chǔ)。

    猜你喜歡
    原型撓度載荷
    交通運(yùn)輸部海事局“新一代衛(wèi)星AIS驗證載荷”成功發(fā)射
    水上消防(2022年2期)2022-07-22 08:45:00
    包裹的一切
    Spontaneous multivessel coronary artery spasm diagnosed with intravascular ultrasound imaging:A case report
    《哈姆雷特》的《圣經(jīng)》敘事原型考證
    滾轉(zhuǎn)機(jī)動載荷減緩風(fēng)洞試驗
    論《西藏隱秘歲月》的原型復(fù)現(xiàn)
    原型理論分析“門”
    人間(2015年20期)2016-01-04 12:47:08
    一種基于白噪聲響應(yīng)的隨機(jī)載荷譜識別方法
    底排藥受力載荷及其分布規(guī)律
    懸高測量在橋梁撓度快速檢測中的應(yīng)用
    老汉色av国产亚洲站长工具| 国产精品亚洲av一区麻豆| 黄色视频,在线免费观看| 91字幕亚洲| 中文字幕高清在线视频| 欧美精品一区二区免费开放| 久久中文看片网| 天堂俺去俺来也www色官网| 人人妻人人添人人爽欧美一区卜| 欧美日韩一级在线毛片| 久久久久精品人妻al黑| 成人国语在线视频| 黄色视频不卡| 久久人人97超碰香蕉20202| 黑人欧美特级aaaaaa片| 欧美亚洲 丝袜 人妻 在线| avwww免费| aaaaa片日本免费| 日韩视频一区二区在线观看| 久久av网站| 欧美黑人欧美精品刺激| 久久久国产精品麻豆| 电影成人av| 亚洲色图 男人天堂 中文字幕| 99久久精品国产亚洲精品| 国产精品一区二区免费欧美| 国产精品98久久久久久宅男小说| 日本欧美视频一区| 国产欧美亚洲国产| bbb黄色大片| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久久久精品电影小说| 极品少妇高潮喷水抽搐| 肉色欧美久久久久久久蜜桃| 国产精品秋霞免费鲁丝片| 淫妇啪啪啪对白视频| 亚洲精品一卡2卡三卡4卡5卡| 啦啦啦免费观看视频1| 国产精品自产拍在线观看55亚洲 | 亚洲伊人久久精品综合| 亚洲av成人一区二区三| 蜜桃国产av成人99| 亚洲精品粉嫩美女一区| 午夜老司机福利片| 国产亚洲精品一区二区www | 极品教师在线免费播放| 国产成人免费观看mmmm| 在线 av 中文字幕| 国产欧美日韩一区二区三| 纵有疾风起免费观看全集完整版| 久久久精品94久久精品| 亚洲欧美激情在线| 男女床上黄色一级片免费看| 在线观看免费日韩欧美大片| 久久热在线av| 麻豆乱淫一区二区| 欧美日本中文国产一区发布| 在线观看免费午夜福利视频| 亚洲精品久久午夜乱码| 欧美亚洲日本最大视频资源| 日韩熟女老妇一区二区性免费视频| 成年人黄色毛片网站| 搡老乐熟女国产| 少妇 在线观看| 久久久精品免费免费高清| 满18在线观看网站| 一区在线观看完整版| 99热国产这里只有精品6| 黄色片一级片一级黄色片| 性少妇av在线| 91大片在线观看| 日本av手机在线免费观看| 午夜福利视频精品| 亚洲专区国产一区二区| 另类精品久久| 亚洲成人免费av在线播放| 亚洲欧美色中文字幕在线| 欧美黑人精品巨大| 91成年电影在线观看| 亚洲av第一区精品v没综合| cao死你这个sao货| 桃红色精品国产亚洲av| 国产精品美女特级片免费视频播放器 | 午夜精品国产一区二区电影| 日韩视频一区二区在线观看| 777米奇影视久久| 亚洲午夜精品一区,二区,三区| 成年动漫av网址| 夜夜爽天天搞| 50天的宝宝边吃奶边哭怎么回事| 亚洲五月色婷婷综合| 黑人巨大精品欧美一区二区蜜桃| 久久久精品免费免费高清| av有码第一页| av网站免费在线观看视频| 男人操女人黄网站| av又黄又爽大尺度在线免费看| 久久久久精品国产欧美久久久| 在线播放国产精品三级| 久久青草综合色| 国产老妇伦熟女老妇高清| 欧美日韩黄片免| 18禁观看日本| 精品人妻1区二区| 国产成人系列免费观看| av在线播放免费不卡| 欧美日本中文国产一区发布| 男女免费视频国产| 麻豆av在线久日| 亚洲熟女精品中文字幕| 国产精品亚洲av一区麻豆| 久久人妻av系列| 美国免费a级毛片| 12—13女人毛片做爰片一| 久久狼人影院| 男女无遮挡免费网站观看| 亚洲av电影在线进入| 又大又爽又粗| 国产精品秋霞免费鲁丝片| 成人18禁高潮啪啪吃奶动态图| 国产av国产精品国产| 国产精品99久久99久久久不卡| av天堂在线播放| 两个人免费观看高清视频| 99精国产麻豆久久婷婷| 超色免费av| 狂野欧美激情性xxxx| 国产成人av激情在线播放| 一本一本久久a久久精品综合妖精| 久久天躁狠狠躁夜夜2o2o| 菩萨蛮人人尽说江南好唐韦庄| 最新美女视频免费是黄的| 欧美人与性动交α欧美软件| 久久婷婷成人综合色麻豆| 美女午夜性视频免费| 国产一区二区激情短视频| 国产黄色免费在线视频| 精品久久久久久电影网| 最新美女视频免费是黄的| 激情在线观看视频在线高清 | 精品福利观看| 欧美激情 高清一区二区三区| 久9热在线精品视频| 成人精品一区二区免费| 欧美日韩亚洲高清精品| 日韩欧美一区二区三区在线观看 | 国产精品久久久av美女十八| 欧美av亚洲av综合av国产av| 亚洲久久久国产精品| 国产有黄有色有爽视频| 精品一区二区三区四区五区乱码| 人人妻人人澡人人爽人人夜夜| 最近最新中文字幕大全免费视频| 嫩草影视91久久| 一区在线观看完整版| 国产精品九九99| 丝袜喷水一区| 精品久久久久久久毛片微露脸| 精品国内亚洲2022精品成人 | 成年动漫av网址| 亚洲国产欧美网| 人妻一区二区av| 中文字幕色久视频| 国产一区二区三区视频了| 亚洲一码二码三码区别大吗| 亚洲av成人一区二区三| 两个人免费观看高清视频| 亚洲精品自拍成人| 麻豆成人av在线观看| 亚洲avbb在线观看| 亚洲男人天堂网一区| 久久影院123| 欧美黄色片欧美黄色片| 精品一区二区三区av网在线观看 | 日韩一区二区三区影片| 日韩大码丰满熟妇| 两性午夜刺激爽爽歪歪视频在线观看 | 看免费av毛片| 波多野结衣一区麻豆| 日韩欧美三级三区| 婷婷成人精品国产| 99国产精品99久久久久| 最新在线观看一区二区三区| 欧美日韩一级在线毛片| 日韩一卡2卡3卡4卡2021年| 精品一区二区三区四区五区乱码| 国产av又大| 汤姆久久久久久久影院中文字幕| 国产成人精品无人区| 女性生殖器流出的白浆| 欧美老熟妇乱子伦牲交| 在线观看免费视频网站a站| 999精品在线视频| 欧美国产精品va在线观看不卡| 久久久国产一区二区| 桃红色精品国产亚洲av| 国产有黄有色有爽视频| 少妇精品久久久久久久| 一边摸一边做爽爽视频免费| 最新美女视频免费是黄的| 在线观看免费午夜福利视频| 黑丝袜美女国产一区| 久久久国产欧美日韩av| 在线观看66精品国产| 精品一区二区三区av网在线观看 | 高清av免费在线| 操出白浆在线播放| 在线观看免费日韩欧美大片| 亚洲av第一区精品v没综合| 国产欧美日韩一区二区三| 黑人巨大精品欧美一区二区蜜桃| 老司机在亚洲福利影院| 亚洲av片天天在线观看| 精品一区二区三区av网在线观看 | 最新美女视频免费是黄的| 亚洲自偷自拍图片 自拍| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品av久久久久免费| 欧美人与性动交α欧美精品济南到| 制服诱惑二区| 成人手机av| 欧美在线一区亚洲| 亚洲欧美一区二区三区黑人| 天堂动漫精品| 亚洲精品一二三| 亚洲第一av免费看| 乱人伦中国视频| 极品少妇高潮喷水抽搐| 99久久人妻综合| 免费在线观看影片大全网站| 黄网站色视频无遮挡免费观看| 岛国在线观看网站| 国产在视频线精品| 99久久99久久久精品蜜桃| 99国产精品一区二区三区| 在线 av 中文字幕| 老司机在亚洲福利影院| 欧美日韩一级在线毛片| 香蕉丝袜av| 国产伦人伦偷精品视频| 一级,二级,三级黄色视频| 国产国语露脸激情在线看| 777米奇影视久久| 成人三级做爰电影| 亚洲欧洲日产国产| 国产精品欧美亚洲77777| 国产成人啪精品午夜网站| 最新在线观看一区二区三区| 精品久久久久久电影网| 免费不卡黄色视频| 亚洲综合色网址| 久久久久久久久免费视频了| 男男h啪啪无遮挡| 免费日韩欧美在线观看| 亚洲午夜精品一区,二区,三区| 777久久人妻少妇嫩草av网站| 国产精品免费视频内射| 一区二区三区激情视频| 亚洲精品中文字幕在线视频| 国产一区二区三区视频了| 亚洲伊人色综图| 国产男靠女视频免费网站| 在线观看一区二区三区激情| 日韩欧美一区视频在线观看| 免费在线观看完整版高清| 成年人免费黄色播放视频| 国内毛片毛片毛片毛片毛片| 最近最新中文字幕大全免费视频| 视频在线观看一区二区三区| 欧美黄色淫秽网站| 亚洲av片天天在线观看| 三上悠亚av全集在线观看| 国产在线视频一区二区| 欧美亚洲 丝袜 人妻 在线| 国产亚洲午夜精品一区二区久久| 久久精品亚洲精品国产色婷小说| 伊人久久大香线蕉亚洲五| 久久久国产精品麻豆| 99香蕉大伊视频| 50天的宝宝边吃奶边哭怎么回事| 黄色视频在线播放观看不卡| 国产精品久久久av美女十八| 久久中文字幕一级| 欧美黄色淫秽网站| 欧美性长视频在线观看| 日本欧美视频一区| 国产亚洲精品久久久久5区| 菩萨蛮人人尽说江南好唐韦庄| 丁香六月欧美| 少妇粗大呻吟视频| 美女高潮喷水抽搐中文字幕| 久久天躁狠狠躁夜夜2o2o| svipshipincom国产片| 欧美性长视频在线观看| 亚洲少妇的诱惑av| 成人免费观看视频高清| 欧美乱妇无乱码| 国产人伦9x9x在线观看| 九色亚洲精品在线播放| 一区福利在线观看| 天天操日日干夜夜撸| 亚洲第一欧美日韩一区二区三区 | 中文亚洲av片在线观看爽 | 精品一品国产午夜福利视频| 51午夜福利影视在线观看| 十分钟在线观看高清视频www| 咕卡用的链子| 亚洲少妇的诱惑av| 亚洲九九香蕉| 桃花免费在线播放| 一区二区三区国产精品乱码| 国产精品久久久久久精品古装| 久久毛片免费看一区二区三区| a级毛片在线看网站| 大片电影免费在线观看免费| 久久久久网色| 久久狼人影院| 国产福利在线免费观看视频| 亚洲第一欧美日韩一区二区三区 | 亚洲少妇的诱惑av| 亚洲性夜色夜夜综合| 日本撒尿小便嘘嘘汇集6| 9191精品国产免费久久| 中文字幕人妻丝袜制服| 亚洲伊人久久精品综合| 精品少妇内射三级| 精品人妻1区二区| 老司机影院毛片| 无人区码免费观看不卡 | 欧美老熟妇乱子伦牲交| 国产黄色免费在线视频| 法律面前人人平等表现在哪些方面| 亚洲午夜精品一区,二区,三区| 欧美变态另类bdsm刘玥| 视频在线观看一区二区三区| av一本久久久久| 久热爱精品视频在线9| 亚洲精品一二三| 两性夫妻黄色片| 怎么达到女性高潮| 欧美亚洲日本最大视频资源| 成年人免费黄色播放视频| 精品国产一区二区久久| 午夜91福利影院| 在线观看一区二区三区激情| 亚洲av欧美aⅴ国产| 丰满人妻熟妇乱又伦精品不卡| 亚洲三区欧美一区| 亚洲精品美女久久av网站| 国产精品秋霞免费鲁丝片| 99九九在线精品视频| 免费少妇av软件| 免费看十八禁软件| 国产成人av激情在线播放| 中文字幕制服av| 欧美日韩亚洲综合一区二区三区_| av又黄又爽大尺度在线免费看| 一边摸一边抽搐一进一出视频| 岛国在线观看网站| 80岁老熟妇乱子伦牲交| 大香蕉久久成人网| 日本av手机在线免费观看| 亚洲精品在线观看二区| 天天操日日干夜夜撸| 国产成人免费无遮挡视频| 悠悠久久av| 午夜福利在线免费观看网站| 亚洲成国产人片在线观看| 亚洲自偷自拍图片 自拍| av欧美777| 超色免费av| 亚洲精品久久午夜乱码| 91字幕亚洲| 丝袜美腿诱惑在线| 一级片'在线观看视频| 精品福利观看| 99精品欧美一区二区三区四区| 亚洲第一欧美日韩一区二区三区 | 国产主播在线观看一区二区| 黄片大片在线免费观看| 日日爽夜夜爽网站| 91成人精品电影| cao死你这个sao货| 精品国产乱子伦一区二区三区| 狠狠精品人妻久久久久久综合| 成在线人永久免费视频| 免费久久久久久久精品成人欧美视频| 在线观看免费视频日本深夜| www.999成人在线观看| 老司机午夜十八禁免费视频| 成年人午夜在线观看视频| 男女床上黄色一级片免费看| 99热网站在线观看| 黄色成人免费大全| 国产午夜精品久久久久久| 国产97色在线日韩免费| 九色亚洲精品在线播放| 岛国毛片在线播放| av又黄又爽大尺度在线免费看| 丰满饥渴人妻一区二区三| 在线天堂中文资源库| 18禁裸乳无遮挡动漫免费视频| 一二三四在线观看免费中文在| 亚洲精品美女久久久久99蜜臀| 99久久精品国产亚洲精品| 激情在线观看视频在线高清 | 久久99热这里只有精品18| 美女cb高潮喷水在线观看 | 亚洲国产欧洲综合997久久,| 国产日本99.免费观看| www日本在线高清视频| 国产精品免费一区二区三区在线| 欧美日韩福利视频一区二区| 免费看光身美女| 美女高潮喷水抽搐中文字幕| 很黄的视频免费| 久久国产精品影院| 成人特级av手机在线观看| 午夜福利成人在线免费观看| 身体一侧抽搐| 99国产综合亚洲精品| 欧美大码av| 婷婷精品国产亚洲av在线| cao死你这个sao货| 精品国产乱子伦一区二区三区| 亚洲美女黄片视频| 中国美女看黄片| 一级毛片精品| 久久精品aⅴ一区二区三区四区| 亚洲人成电影免费在线| 久久久久国产精品人妻aⅴ院| 亚洲熟妇熟女久久| 校园春色视频在线观看| 日韩免费av在线播放| 69av精品久久久久久| 窝窝影院91人妻| 国产精品香港三级国产av潘金莲| 成人av在线播放网站| 老熟妇乱子伦视频在线观看| 琪琪午夜伦伦电影理论片6080| 精品久久久久久成人av| 久久久久久人人人人人| 麻豆av在线久日| 亚洲午夜理论影院| 天堂√8在线中文| 国产单亲对白刺激| 欧美一区二区精品小视频在线| 日韩av在线大香蕉| 小蜜桃在线观看免费完整版高清| 久久精品aⅴ一区二区三区四区| 中文字幕精品亚洲无线码一区| 九九在线视频观看精品| 亚洲av第一区精品v没综合| 日日夜夜操网爽| 在线永久观看黄色视频| 国产视频内射| 国产极品精品免费视频能看的| 国产高清有码在线观看视频| 国产欧美日韩精品一区二区| 最新中文字幕久久久久 | 看免费av毛片| 国产伦精品一区二区三区视频9 | 国产精品av久久久久免费| 夜夜躁狠狠躁天天躁| 美女大奶头视频| 国产美女午夜福利| 日本成人三级电影网站| 亚洲性夜色夜夜综合| 国产高清三级在线| 97人妻精品一区二区三区麻豆| 亚洲人成电影免费在线| 国产精品1区2区在线观看.| 亚洲狠狠婷婷综合久久图片| 国产探花在线观看一区二区| 他把我摸到了高潮在线观看| 欧美中文日本在线观看视频| 国内精品美女久久久久久| 色播亚洲综合网| 欧美乱色亚洲激情| АⅤ资源中文在线天堂| 日本在线视频免费播放| 精品久久久久久,| 日韩有码中文字幕| 国产亚洲精品久久久久久毛片| 51午夜福利影视在线观看| 小蜜桃在线观看免费完整版高清| 久久天堂一区二区三区四区| 老司机午夜福利在线观看视频| 久久久久久九九精品二区国产| 亚洲五月天丁香| 99国产精品一区二区三区| 久久香蕉国产精品| 最近最新中文字幕大全电影3| 18禁黄网站禁片免费观看直播| 狂野欧美白嫩少妇大欣赏| 999久久久精品免费观看国产| 亚洲av第一区精品v没综合| 少妇的逼水好多| 搡老熟女国产l中国老女人| 天堂动漫精品| 伦理电影免费视频| av国产免费在线观看| 亚洲精品456在线播放app | 身体一侧抽搐| 午夜免费激情av| 免费av不卡在线播放| 最近最新中文字幕大全免费视频| 欧美乱色亚洲激情| 午夜福利免费观看在线| 啦啦啦观看免费观看视频高清| 他把我摸到了高潮在线观看| 国产单亲对白刺激| 日韩免费av在线播放| 禁无遮挡网站| 黄色成人免费大全| 亚洲av熟女| 欧美乱色亚洲激情| 露出奶头的视频| 怎么达到女性高潮| 久久中文看片网| 级片在线观看| 午夜福利欧美成人| 美女高潮的动态| 国产男靠女视频免费网站| 午夜激情欧美在线| 成人av一区二区三区在线看| 日本一二三区视频观看| 三级男女做爰猛烈吃奶摸视频| 国产精品久久久久久人妻精品电影| 超碰成人久久| 叶爱在线成人免费视频播放| 久久精品国产综合久久久| 99在线人妻在线中文字幕| 国产高清三级在线| 黄色女人牲交| 美女扒开内裤让男人捅视频| 男女午夜视频在线观看| 亚洲人与动物交配视频| 老汉色∧v一级毛片| 久久久久九九精品影院| av天堂中文字幕网| 国产精品久久久久久人妻精品电影| 99国产精品99久久久久| 久久久国产欧美日韩av| 成人性生交大片免费视频hd| 老熟妇乱子伦视频在线观看| 国产高清视频在线观看网站| 欧美日韩亚洲国产一区二区在线观看| 无遮挡黄片免费观看| aaaaa片日本免费| 亚洲精品在线观看二区| 在线免费观看不下载黄p国产 | 精品人妻1区二区| 国内精品久久久久久久电影| av中文乱码字幕在线| 黄色片一级片一级黄色片| ponron亚洲| 日韩有码中文字幕| 日韩人妻高清精品专区| 99精品久久久久人妻精品| 午夜视频精品福利| 国产淫片久久久久久久久 | 听说在线观看完整版免费高清| 亚洲成av人片在线播放无| 成在线人永久免费视频| 国产熟女xx| 亚洲狠狠婷婷综合久久图片| 国产成人精品久久二区二区免费| 琪琪午夜伦伦电影理论片6080| 一级毛片高清免费大全| 视频区欧美日本亚洲| 在线国产一区二区在线| 淫秽高清视频在线观看| 波多野结衣巨乳人妻| 天天添夜夜摸| 91在线观看av| 偷拍熟女少妇极品色| 国产午夜精品论理片| 听说在线观看完整版免费高清| 搡老熟女国产l中国老女人| 淫秽高清视频在线观看| 亚洲av片天天在线观看| or卡值多少钱| 国产精品亚洲美女久久久| 国产伦人伦偷精品视频| 757午夜福利合集在线观看| 国产精品久久久久久人妻精品电影| 欧美不卡视频在线免费观看| 天天躁日日操中文字幕| 亚洲专区国产一区二区| 免费看十八禁软件| 成人特级黄色片久久久久久久| 国产精华一区二区三区| 悠悠久久av| 国产麻豆成人av免费视频| 国产精品久久久久久亚洲av鲁大| 亚洲人成伊人成综合网2020| 亚洲自拍偷在线| 免费看日本二区| 久久精品国产亚洲av香蕉五月| 日本三级黄在线观看| 18禁国产床啪视频网站| 亚洲在线观看片| 亚洲国产欧美人成| 久久精品国产综合久久久| 日日干狠狠操夜夜爽| 国语自产精品视频在线第100页| 欧美成人一区二区免费高清观看 | 中亚洲国语对白在线视频| 亚洲电影在线观看av| 在线观看免费午夜福利视频| 免费看a级黄色片| 人妻久久中文字幕网| 男女下面进入的视频免费午夜| www.熟女人妻精品国产|