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

    基于迭代RMSE 法的約束阻尼板動(dòng)力特性分析

    2023-12-21 09:43:34劉全民葉孝意宋立忠孫逸飛劉林芽
    關(guān)鍵詞:固有頻率阻尼損耗

    劉全民,葉孝意,宋立忠,孫逸飛,劉林芽

    (華東交通大學(xué)軌道交通基礎(chǔ)設(shè)施性能監(jiān)測與保障國家重點(diǎn)實(shí)驗(yàn)室,江西 南昌 330013)

    黏彈性阻尼結(jié)構(gòu)具有高效的減振降噪能力及質(zhì)量輕等優(yōu)點(diǎn),廣泛應(yīng)用在航空航天、土木工程、交通運(yùn)輸?shù)阮I(lǐng)域[1-2],如在車輪表面粘貼約束阻尼層以增大車輪阻尼[3],在飛機(jī)壁板敷設(shè)約束阻尼層控制中低頻結(jié)構(gòu)振動(dòng)和噪聲[4]等.黏彈性阻尼材料具有黏性流體和彈性固體2 種特性[5]:黏性流體在受到外力時(shí)耗散能量,而彈性固體則是儲(chǔ)存能量;彈性材料的應(yīng)力和應(yīng)變幾乎是同時(shí)增加和減小的,基本不存在相位差[6];黏彈性材料的應(yīng)變滯后于應(yīng)力,應(yīng)力應(yīng)變曲線呈橢圓形,橢圓形滯回曲線所包含的面積表示黏彈性材料耗散的能量.

    根據(jù)實(shí)際所需,表面敷設(shè)黏彈性阻尼結(jié)構(gòu)的形式主要有以下4 種:1)自由阻尼層結(jié)構(gòu)(FLD),直接將黏彈性材料敷設(shè)在需要減振的結(jié)構(gòu)表面上;2)被動(dòng)約束阻尼層結(jié)構(gòu)(PCLD),在自由阻尼層表面再敷設(shè)一層彈性約束層;3)主動(dòng)約束阻尼層結(jié)構(gòu)(ACLD),采用可控壓電材料約束層代替PCLD 中不可控約束層;4)智能約束阻尼層結(jié)構(gòu)(SCLD),將傳統(tǒng)的被動(dòng)約束阻尼層與主動(dòng)控制技術(shù)相結(jié)合的新興阻尼減振降噪結(jié)構(gòu)[7].被動(dòng)約束阻尼結(jié)構(gòu)是一種將黏彈性阻尼材料鋪設(shè)在結(jié)構(gòu)層和具有較大剛度的約束層之間的結(jié)構(gòu)[8-9],要求阻尼層與結(jié)構(gòu)層之間、約束層與阻尼層之間粘結(jié)牢固,阻尼層隨結(jié)構(gòu)層彎曲振動(dòng)產(chǎn)生拉壓變形,同時(shí)約束層能夠約束住阻尼層的拉壓變形,從而使阻尼層產(chǎn)生剪切變形,耗散更多的振動(dòng)能量實(shí)現(xiàn)減振降噪[10],具有減振效果好、可靠性高、成本低等優(yōu)點(diǎn),適用于鋼結(jié)構(gòu)等的振動(dòng)控制.

    衡量黏彈性材料耗能性能的主要指標(biāo)是材料的損耗因子,約束阻尼結(jié)構(gòu)模態(tài)損耗因子則是衡量阻尼處理是否合理、減振效果是否顯著的關(guān)鍵,不少學(xué)者對此進(jìn)行了積極的探索:郭中澤等[11]利用Reissner夾層板理論計(jì)算了約束阻尼結(jié)構(gòu)的固有頻率和模態(tài)損耗因子;Johnson 和Kienholz[12]采用模態(tài)應(yīng)變能法(MSE 法)計(jì)算了復(fù)合夾層梁前5 階模態(tài)損耗因子;Ren 和Zhao[13]采用有限元法計(jì)算了不考慮阻尼材料頻變特性時(shí)約束阻尼板前5 階模態(tài)損耗因子;Rao等[14]用直接頻率響應(yīng)技術(shù)和模態(tài)應(yīng)變能方法求解黏彈性約束阻尼梁損耗因子;Ravi 等[15-16]將黏彈性材料的力學(xué)性能參數(shù)(剪切模量和損耗因子)考慮成常數(shù),采用有限元法建立了彈性-黏彈性復(fù)合梁動(dòng)力學(xué)模型,分析了結(jié)構(gòu)前2 階固有頻率和模態(tài)損耗因子.

    已有的約束阻尼結(jié)構(gòu)動(dòng)力分析大多將黏彈性材料的剪切模量和損耗因子看作常數(shù),而黏彈性阻尼材料特性隨頻率變化顯著,這將導(dǎo)致在對約束阻尼結(jié)構(gòu)動(dòng)力特性進(jìn)行計(jì)算時(shí)對計(jì)算結(jié)果帶來誤差.MSE法計(jì)算效率高,較多地應(yīng)用在約束阻尼結(jié)構(gòu)模態(tài)損耗因子計(jì)算中,但已有研究往往忽略了黏彈性材料虛剛度的影響.為了減小阻尼層黏彈性材料虛剛度和參數(shù)頻變特性對約束阻尼結(jié)構(gòu)模態(tài)損耗因子計(jì)算造成的誤差,計(jì)算出的約束阻尼結(jié)構(gòu)更接近真實(shí)的模態(tài)損耗因子,本文結(jié)合修正模態(tài)應(yīng)變能法(RMSE法)和迭代算法,探討?zhàn)椥圆牧咸搫偠燃皡?shù)頻變特性對約束阻尼板的振型、固有頻率和模態(tài)損耗因子的影響.

    1 模態(tài)應(yīng)變能法

    模態(tài)應(yīng)變能法被認(rèn)為是黏彈性阻尼結(jié)構(gòu)模態(tài)損耗因子計(jì)算中最實(shí)用的一種方法,反映了整個(gè)結(jié)構(gòu)在敷設(shè)約束阻尼層后,每一階振型對于諧振峰的抑制作用.其主要做法是假定黏彈性阻尼結(jié)構(gòu)的模態(tài)近似于對應(yīng)結(jié)構(gòu)的無阻尼模態(tài),有限元法計(jì)算時(shí),只需將黏彈性阻尼材料當(dāng)作具有實(shí)剛度模量的純彈性體,從而避免了復(fù)雜的復(fù)特征值計(jì)算.該方法同時(shí)具有較高的計(jì)算效率和計(jì)算精度[17].

    約束阻尼結(jié)構(gòu)形式作以下假設(shè)[18-20]:1)各層材料為線彈性、均勻、各向同性;2)結(jié)構(gòu)能正常工作,各層之間粘結(jié)非常牢固,在彎曲過程中各層之間不會(huì)產(chǎn)生滑動(dòng),將約束阻尼層整體作為一個(gè)耗能元件,不考慮能量在約束層和阻尼層之間的流動(dòng);3)垂直于板中間的各層無擠壓;4)黏彈性材料對結(jié)構(gòu)阻尼貢獻(xiàn)是主要的.

    一般地,無阻尼自由振動(dòng)運(yùn)動(dòng)方程[21]為

    式(1)的特征方程為

    其中:λ 為特征值;Φ 為特征向量.

    當(dāng)結(jié)構(gòu)為含有黏彈性材料的復(fù)合結(jié)構(gòu),結(jié)構(gòu)剛度矩陣KG為復(fù)剛度矩陣,如式(3).

    式中:Ke為復(fù)合結(jié)構(gòu)彈性部分剛度矩陣;KvR為復(fù)合結(jié)構(gòu)黏彈性部分剛度矩陣實(shí)部;KvI為復(fù)合結(jié)構(gòu)黏彈性部分剛度矩陣虛部;ηv為黏彈性阻尼材料損耗因子.

    黏彈性復(fù)合結(jié)構(gòu)運(yùn)動(dòng)方程為

    復(fù)合結(jié)構(gòu)特征方程為

    式中:λG為復(fù)合結(jié)構(gòu)復(fù)特征值,如式(6);ΦG為復(fù)合結(jié)構(gòu)復(fù)特征向量.

    式中:λR、λI分別為復(fù)特征值的實(shí)部和虛部;η 為復(fù)合結(jié)構(gòu)模態(tài)損耗因子.

    將式(6)代入式(5)得

    將 ΦG由實(shí)特征向量 ΦR近似,則

    根據(jù)式(9)兩邊實(shí)部與虛部分別相等,可得

    將式(10)代入式(11),可得復(fù)合結(jié)構(gòu)模態(tài)損耗因子的計(jì)算表達(dá)式為

    聯(lián)合式(13)、(14)可得

    2 基于迭代的RMSE 算法

    本文采用基于RMSE 的迭代算法求解復(fù)合結(jié)構(gòu)的頻率和損耗因子有以下2 個(gè)原因:1)一般而言,約束阻尼結(jié)構(gòu)不同階次下模態(tài)損耗因子不同,不同階次采用相同修正系數(shù)可能導(dǎo)致修正不足或過量;2)黏彈性材料的剪切模量及損耗因子隨頻率變化較大,這使得不同模態(tài)階次下剪切模量、材料損耗因子具有顯著差異,這對模態(tài)損耗因子和固有頻率的計(jì)算有很大影響.RMSE 方法僅在計(jì)算特征向量時(shí)使用復(fù)剛度的絕對值來考慮虛剛度的影響,從而避免特征向量和剛度雙重修正引起的過度修正.黏彈性材料參數(shù)是頻變的,而約束阻尼結(jié)構(gòu)的固有頻率又與材料參數(shù)有關(guān),而這可以通過迭代算法來實(shí)現(xiàn).通過迭代計(jì)算考慮黏彈性阻尼材料的頻變特性,可更準(zhǔn)確地分析模態(tài)損耗因子和固有頻率.

    該迭代法的具體實(shí)現(xiàn)步驟:

    步驟1約束阻尼結(jié)構(gòu)進(jìn)行實(shí)模態(tài)分析得各階固有頻率.

    步驟2選取該頻率下阻尼層材料損耗因子和剪切模量,進(jìn)行模態(tài)分析得到新的固有頻率.

    步驟3將前、后2 組的固有頻率進(jìn)行比較,若滿足誤差公式(其中,fi,j為第j次模態(tài)分析后約束阻尼結(jié)構(gòu)的第i階頻率),則輸出固有頻率,并采用RMSE 方法計(jì)算該階模態(tài)損耗因子,否則,重復(fù)步驟1 和步驟2.

    3 計(jì)算方法驗(yàn)證

    為驗(yàn)證本文提出的基于RMSE 的迭代法較MSE 迭代法的準(zhǔn)確性,建立與文獻(xiàn)[23]相同的模型,計(jì)算結(jié)構(gòu)固有頻率和模態(tài)損耗因子,并與MSE 迭代法及文獻(xiàn)[23]的模態(tài)試驗(yàn)測試數(shù)據(jù)對比分析.本文RMSE 迭代法考慮了阻尼層虛剛度的貢獻(xiàn),這與已有MSE 迭代法顯著不同.

    文獻(xiàn)[23]試驗(yàn)材料為ZN-3 黏彈性阻尼材料,試驗(yàn)的頻率范圍為0.1~200.0 Hz,且在該范圍內(nèi)材料的儲(chǔ)能模量E和損耗因子 ηv、頻率f的關(guān)系為

    文獻(xiàn)[23]中試驗(yàn)采用的結(jié)構(gòu)形式是一邊固定的約束阻尼板.本文均采用實(shí)體單元模擬結(jié)構(gòu)層、約束層和阻尼層,計(jì)算模型邊界條件為懸臂約束,如圖1所示.基層和約束層彈性模量均為70 GPa,模型其他參數(shù)如表1 所示,計(jì)算結(jié)果對比如表2 所示.

    表1 文獻(xiàn)[23]模型參數(shù)Tab.1 Model parameters in reference [23]

    表2 試驗(yàn)與計(jì)算結(jié)果比較Tab.2 Comparison of experimental and calculated results

    圖1 結(jié)構(gòu)模型Fig.1 Structural model

    從表2 中可以看出,對比文獻(xiàn)[23]的實(shí)測數(shù)據(jù),本文算法考慮阻尼層虛剛度貢獻(xiàn)和材料頻變特性后得到的固有頻率和模態(tài)損耗因子誤差更小,表明本文提出的算法較傳統(tǒng)MSE 迭代法的計(jì)算精度更高.由于本文算法涉及的仍然是與傳統(tǒng)MSE 法相同的實(shí)特征求解,通常迭代3~5 次計(jì)算結(jié)果即可收斂,計(jì)算效率雖低于傳統(tǒng)MSE 法,但仍然明顯高于復(fù)特征值法.

    4 約束阻尼板動(dòng)力特性分析

    黏彈性阻尼材料的損耗因子和儲(chǔ)能剪切模量隨頻率變化顯著.本文選用ISD-112 型黏彈性阻尼材料,頻變儲(chǔ)能剪切模量和材料損耗因子如圖2.在0~1 000.0 Hz 內(nèi),儲(chǔ)能剪切模量G和損耗因子ηv與頻率f的關(guān)系[1]為

    圖2 ISD-112 阻尼材料頻變特性Fig.2 Frequency-dependent characteristics of damping material ISD-112

    在采用迭代算法計(jì)算約束阻尼板模態(tài)損耗因子時(shí),利用式(18)、(19)可得到任意頻率下阻尼材料的儲(chǔ)能剪切模量和材料損耗因子.

    計(jì)算模型邊界條件為四邊簡支,基層和約束層彈性模量為68.95 GPa,其他具體物理參數(shù)如表3 所示.

    表3 模型參數(shù)Tab.3 Model parameters

    4.1 模態(tài)振型對比

    由于本文關(guān)心的頻率范圍是0~1.0 kHz,所以選取中間頻率500.0 Hz 對應(yīng)的儲(chǔ)能剪切模量和材料損耗因子作為約束阻尼板非頻變參數(shù)進(jìn)行模態(tài)分析.迭代法則可以選擇任意頻率下的儲(chǔ)能剪切模量和材料損耗因子作為初值,按照第2 節(jié)描述的步驟進(jìn)行迭代模態(tài)分析.

    迭代前后計(jì)算所得的前20 階模態(tài)振型相差不大,現(xiàn)取其中的4 階計(jì)算結(jié)果進(jìn)行對比分析.約束阻尼板結(jié)構(gòu)層第1、12、15、20 階模態(tài)振型結(jié)果如圖3 所示.

    圖3 結(jié)構(gòu)層模態(tài)振型Fig.3 Modal shapes of structural layer

    從圖3 可以看出:迭代前后,各階模態(tài)振型形狀基本相同,僅部分模態(tài)振型相位相差180°.這是由于阻尼層對整個(gè)結(jié)構(gòu)的剛度貢獻(xiàn)不大,所以剪切模量的改變對振型形狀的影響可忽略,但會(huì)改變振型的相位.

    4.2 固有頻率對比

    迭代前后黏彈性約束阻尼板前20 階固有頻率計(jì)算結(jié)果如圖4 所示.由圖可見:對于前10 階,非迭代計(jì)算的固有頻率均高于迭代后的,且差值逐漸減??;對于后10 階,非迭代計(jì)算的固有頻率又低于迭代后的結(jié)果,且差值逐漸增大;在第10 階時(shí),由于該頻率下的阻尼層剪切模量與實(shí)際值最為接近,故迭代前后的差值最小.

    圖4 迭代前后各階固有頻率對比Fig.4 Comparison of intrinsic frequencies before and after iteration

    迭代前后,約束阻尼結(jié)構(gòu)的質(zhì)量矩陣不發(fā)生改變,對于影響固有頻率的因素而言,唯一的變化在于結(jié)構(gòu)的剛度.迭代前后剛度的不同是阻尼層剪切模量變化引起的,阻尼材料的剪切模量隨頻率變化函數(shù)如式(18)所示,阻尼層剪切模量隨頻率增大而增大,所以迭代法計(jì)算時(shí),結(jié)構(gòu)剛度隨階數(shù)增大而增大.非迭代計(jì)算,阻尼層剪切模量只能取定值(本文取500.0 Hz 對應(yīng)的剪切模量),結(jié)構(gòu)剛度矩陣不變,顯然低階模態(tài)高估了阻尼層剪切模量,高階模態(tài)又低估了阻尼層剪切模量.

    以上分析表明,忽略阻尼層材料虛剛度貢獻(xiàn)及其頻變特性,阻尼層模量取定值會(huì)造成固有頻率計(jì)算出現(xiàn)明顯偏差.

    4.3 模態(tài)損耗因子對比

    為了修正由于未考慮黏彈性材料頻變特性而引起的誤差,Johnson 等[12]提出了利用頻率函數(shù)f(ω)對損耗因子進(jìn)行修正,損耗因子修正計(jì)算公式為

    式中:ω為圓頻率;ωi為第i階模態(tài)頻率;Gd(ωi) 為第i階模態(tài)對應(yīng)的黏彈性材料剪切模量;Gd,ref為實(shí)模態(tài)分析中采用的剪切模量;ηi為第i階非迭代模態(tài)損耗因子;ηr為修正模態(tài)損耗因子.

    將采用頻率修正法計(jì)算得到的模態(tài)損耗因子與迭代前后的模態(tài)損耗因子進(jìn)行對比,如圖5 所示.可以看出,非迭代計(jì)算的模態(tài)損耗因子在低頻時(shí)偏差較大,頻率修正法計(jì)算的模態(tài)損耗因子在高頻時(shí)偏差較大.

    圖5 模態(tài)損耗因子對比Fig.5 Comparison of modal loss factors

    將迭代法計(jì)算結(jié)果作為參照,對比分析不進(jìn)行迭代和頻率修正對模態(tài)損耗因子的影響.不進(jìn)行迭代的計(jì)算誤差范圍為-56.00%~3.72%,對于頻率修正后的計(jì)算誤差范圍為-10.20%~3.62%.所以通過頻率修正能在一定程度減小損耗因子的計(jì)算誤差,但在高階時(shí)會(huì)導(dǎo)致修正過度,說明頻率修正得到的損耗因子也存在較大誤差,不能較好地解決不考慮阻尼材料頻變特性帶來的影響.

    3 種方式下求得的模態(tài)損耗因子在第10 階最為接近,原因是第10 階時(shí)結(jié)構(gòu)頻率在500.0 Hz 左右,而不進(jìn)行迭代采用的材料參數(shù)為頻率500.0 Hz對應(yīng)的材料屬性值,所以在該頻段附近的計(jì)算結(jié)果會(huì)比較接近真實(shí)值,才會(huì)出現(xiàn)三者的計(jì)算結(jié)果相近.

    4.4 約束阻尼層參數(shù)分析

    阻尼層厚度和約束層厚度對約束阻尼板減振效果至關(guān)重要,所以本節(jié)主要分析這2 個(gè)參數(shù)對約束阻尼結(jié)構(gòu)的模態(tài)損耗因子的影響規(guī)律.

    4.4.1 阻尼層厚度對比分析

    將約束層厚度取0.762 mm,阻尼層的厚度依次取0.154、0.254、0.354、0.454、0.554 mm,討論在不同厚度下約束阻尼板模態(tài)損耗因子的變化,計(jì)算結(jié)果如圖6 所示.由此可見,隨著阻尼層厚度的增加,損耗因子不斷增大,但增加相同的厚度,損耗因子增量不斷減小.

    圖6 不同阻尼層厚度下CLD 板模態(tài)損耗因子Fig.6 Modal loss factors of CLD plate under various thicknesses of damping layer

    4.4.2 約束層厚度對比分析

    阻尼層厚度取0.254 mm,約束層的厚度依次取0.262、0.562、0.762、1.262、1.762 mm 和2.262 mm,討論在不同厚度下約束阻尼結(jié)構(gòu)模態(tài)損耗因子變化,計(jì)算結(jié)果如圖7 所示.由此可見,隨著約束層厚度增加,結(jié)構(gòu)模態(tài)損耗因子不斷增大,但增大到一定值時(shí),開始減小.這說明并不是約束層厚度越大越好,而是存在最佳值使模態(tài)損耗因子最大.約束層厚度與基層厚度相等時(shí),模態(tài)損耗因子最大.

    圖7 不同約束層厚度下CLD 板模態(tài)損耗因子Fig.7 Modal loss factors of CLD plate under various thicknesses of constraint layer

    5 結(jié)論

    1)通過與相關(guān)文獻(xiàn)試驗(yàn)實(shí)測數(shù)據(jù)對比驗(yàn)證了本文所提出算法的準(zhǔn)確性,并且較已有MSE 迭代法具有更高的計(jì)算精度,說明有必要考慮阻尼層虛剛度的貢獻(xiàn).

    2)不考慮阻尼材料參數(shù)頻變特性,結(jié)構(gòu)層模態(tài)振型形狀基本不變,但部分模態(tài)的振型相位相反.

    3)對于本文選取的阻尼材料而言,阻尼層剪切模量取中間頻段值會(huì)造成固有頻率在低階模態(tài)時(shí)計(jì)算結(jié)果偏大,在高階模態(tài)時(shí)偏小.

    4)不考慮阻尼材料參數(shù)頻變特性計(jì)算得到的模態(tài)損耗因子只有在中間頻段處較為準(zhǔn)確,其他頻段處有較大誤差,特別是低階模態(tài)對應(yīng)頻段.

    5)在一定范圍內(nèi),結(jié)構(gòu)損耗因子與阻尼層厚度成正比.對約束層而言,其厚度與基層厚度相等時(shí),結(jié)構(gòu)模態(tài)損耗因子最大.

    猜你喜歡
    固有頻率阻尼損耗
    現(xiàn)場測定大型水輪發(fā)電機(jī)組軸系的固有頻率
    N維不可壓無阻尼Oldroyd-B模型的最優(yōu)衰減
    關(guān)于具有阻尼項(xiàng)的擴(kuò)散方程
    具有非線性阻尼的Navier-Stokes-Voigt方程的拉回吸引子
    自我損耗理論視角下的編輯審讀
    新聞傳播(2016年11期)2016-07-10 12:04:01
    具阻尼項(xiàng)的Boussinesq型方程的長時(shí)間行為
    變壓器附加損耗對負(fù)載損耗的影響
    非隔離型單相光伏并網(wǎng)逆變器的功率損耗研究
    總溫總壓測頭模態(tài)振型變化規(guī)律研究
    A novel functional electrical stimulation-control system for restoring motor function of post-stroke hemiplegic patients
    亚洲av一区综合| 国产爱豆传媒在线观看| 国产爽快片一区二区三区| 亚洲最大成人av| 国产成人精品一,二区| 嫩草影院精品99| 日本色播在线视频| 免费看a级黄色片| 亚洲精品亚洲一区二区| 国内精品宾馆在线| av在线播放精品| 日本av手机在线免费观看| 网址你懂的国产日韩在线| 国产成人freesex在线| 男女啪啪激烈高潮av片| 有码 亚洲区| 国产精品熟女久久久久浪| 极品少妇高潮喷水抽搐| 亚洲成色77777| 日日啪夜夜撸| 日韩av免费高清视频| 国产成人免费无遮挡视频| 99热全是精品| 欧美变态另类bdsm刘玥| 午夜福利在线观看免费完整高清在| 精品国产三级普通话版| 嫩草影院新地址| 久久精品国产亚洲网站| 国产精品人妻久久久久久| 天天躁夜夜躁狠狠久久av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国精品久久久久久国模美| 色视频www国产| 久久久精品免费免费高清| 成人毛片60女人毛片免费| 国产精品一及| 丰满少妇做爰视频| videos熟女内射| 亚洲欧美日韩东京热| 精品人妻视频免费看| 久久人人爽av亚洲精品天堂 | 日韩一区二区视频免费看| 日本免费在线观看一区| 亚洲av免费在线观看| 亚洲高清免费不卡视频| 日本三级黄在线观看| 九草在线视频观看| 亚洲av成人精品一二三区| 久久99热这里只频精品6学生| 久久人人爽人人片av| 秋霞伦理黄片| av网站免费在线观看视频| 一本色道久久久久久精品综合| 亚洲精品影视一区二区三区av| 久久久久久久国产电影| 看十八女毛片水多多多| 又大又黄又爽视频免费| 久久6这里有精品| 干丝袜人妻中文字幕| 亚洲欧美成人精品一区二区| 日产精品乱码卡一卡2卡三| 少妇的逼好多水| 春色校园在线视频观看| 国产成人精品一,二区| 成年女人看的毛片在线观看| 18禁裸乳无遮挡动漫免费视频 | 国产永久视频网站| 欧美国产精品一级二级三级 | 免费观看的影片在线观看| 久久精品国产鲁丝片午夜精品| 人人妻人人澡人人爽人人夜夜| 黄片wwwwww| 国产高清有码在线观看视频| 极品教师在线视频| 久久精品国产亚洲av涩爱| 国产午夜精品久久久久久一区二区三区| 丰满乱子伦码专区| 久久精品夜色国产| 国产精品国产av在线观看| 2022亚洲国产成人精品| 国产一区二区亚洲精品在线观看| 狂野欧美激情性bbbbbb| 人人妻人人看人人澡| 久久精品人妻少妇| 在线观看国产h片| 国产精品国产三级国产专区5o| 国产精品不卡视频一区二区| 国产精品无大码| 最近中文字幕高清免费大全6| 日韩一区二区三区影片| 伦精品一区二区三区| 日日啪夜夜撸| 深爱激情五月婷婷| 欧美成人精品欧美一级黄| 高清视频免费观看一区二区| 欧美日韩在线观看h| 少妇裸体淫交视频免费看高清| 日本黄色片子视频| 人妻制服诱惑在线中文字幕| 亚洲美女视频黄频| 一级毛片我不卡| 日本一本二区三区精品| 日韩欧美 国产精品| 精品午夜福利在线看| 午夜福利视频1000在线观看| 欧美一级a爱片免费观看看| 欧美xxⅹ黑人| 国产成人福利小说| 午夜激情福利司机影院| 欧美老熟妇乱子伦牲交| 久久99蜜桃精品久久| 日韩强制内射视频| 干丝袜人妻中文字幕| 内射极品少妇av片p| 久久99精品国语久久久| 免费大片18禁| 欧美极品一区二区三区四区| 大陆偷拍与自拍| 亚洲一级一片aⅴ在线观看| 日韩免费高清中文字幕av| 真实男女啪啪啪动态图| 寂寞人妻少妇视频99o| 啦啦啦在线观看免费高清www| 丝袜脚勾引网站| 久久久欧美国产精品| 深夜a级毛片| 亚洲国产精品成人综合色| 久久精品国产亚洲av天美| 在线天堂最新版资源| 国产成人精品婷婷| 欧美xxxx性猛交bbbb| 久久久精品免费免费高清| 婷婷色麻豆天堂久久| 亚洲av福利一区| 亚洲精品久久午夜乱码| 激情五月婷婷亚洲| 麻豆国产97在线/欧美| 在线免费十八禁| 精品久久久久久久人妻蜜臀av| 神马国产精品三级电影在线观看| 欧美成人一区二区免费高清观看| 国产久久久一区二区三区| 全区人妻精品视频| 啦啦啦在线观看免费高清www| 亚洲人与动物交配视频| 伦精品一区二区三区| 搡女人真爽免费视频火全软件| 综合色丁香网| 亚洲精品一区蜜桃| 欧美三级亚洲精品| 九色成人免费人妻av| 97精品久久久久久久久久精品| 亚洲精华国产精华液的使用体验| 免费观看a级毛片全部| 卡戴珊不雅视频在线播放| 久久久久久久久久久丰满| 最近中文字幕2019免费版| 99热网站在线观看| 国产精品99久久99久久久不卡 | 青青草视频在线视频观看| 晚上一个人看的免费电影| 亚洲成人av在线免费| 国产成人免费观看mmmm| 男人爽女人下面视频在线观看| 成人午夜精彩视频在线观看| 精华霜和精华液先用哪个| 97热精品久久久久久| 可以在线观看毛片的网站| 免费高清在线观看视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 伊人久久国产一区二区| 亚洲av日韩在线播放| 国产成人福利小说| 晚上一个人看的免费电影| 99视频精品全部免费 在线| 九色成人免费人妻av| 人妻系列 视频| 色综合色国产| 精品午夜福利在线看| 97精品久久久久久久久久精品| 人人妻人人爽人人添夜夜欢视频 | 久久精品国产亚洲av涩爱| 亚洲欧洲国产日韩| 亚洲最大成人手机在线| 国产精品不卡视频一区二区| 精品国产乱码久久久久久小说| 久久99热这里只频精品6学生| 亚洲精品成人av观看孕妇| 香蕉精品网在线| 男女下面进入的视频免费午夜| 日韩人妻高清精品专区| 国产精品人妻久久久影院| 久久精品久久精品一区二区三区| 五月天丁香电影| 香蕉精品网在线| 一级毛片aaaaaa免费看小| 欧美激情久久久久久爽电影| 国产老妇女一区| 成人毛片60女人毛片免费| 亚洲精品色激情综合| 国产成人免费无遮挡视频| 欧美成人精品欧美一级黄| 国产一级毛片在线| 中文字幕人妻熟人妻熟丝袜美| av专区在线播放| 久久久欧美国产精品| 乱码一卡2卡4卡精品| 国产乱人偷精品视频| 日韩欧美一区视频在线观看 | 99热网站在线观看| 国产成人免费观看mmmm| 久久精品综合一区二区三区| 久久午夜福利片| 亚洲最大成人手机在线| 国产精品久久久久久精品古装| 久久人人爽人人片av| 久久久久久久大尺度免费视频| videos熟女内射| 国产毛片a区久久久久| 我的老师免费观看完整版| 欧美亚洲 丝袜 人妻 在线| 七月丁香在线播放| 国产精品一及| 久久久精品94久久精品| 卡戴珊不雅视频在线播放| 久久久久久久精品精品| 国产亚洲精品久久久com| 午夜老司机福利剧场| 一本久久精品| 国产视频内射| 精品一区二区三卡| 精品少妇黑人巨大在线播放| 午夜视频国产福利| videos熟女内射| 欧美高清性xxxxhd video| 男人添女人高潮全过程视频| 国产精品一区二区三区四区免费观看| 高清午夜精品一区二区三区| 日韩欧美精品v在线| 久久久精品94久久精品| 国产精品国产三级国产专区5o| 成年版毛片免费区| 97在线视频观看| 中文字幕制服av| 寂寞人妻少妇视频99o| 夫妻午夜视频| 少妇 在线观看| 丝袜喷水一区| 国产伦理片在线播放av一区| 国内精品美女久久久久久| 人妻夜夜爽99麻豆av| 在线观看三级黄色| av国产久精品久网站免费入址| 尤物成人国产欧美一区二区三区| 精品酒店卫生间| 少妇丰满av| 午夜免费鲁丝| 中文字幕制服av| 成人亚洲欧美一区二区av| 国产精品久久久久久av不卡| 免费观看的影片在线观看| 新久久久久国产一级毛片| 天天一区二区日本电影三级| 97人妻精品一区二区三区麻豆| 国产精品精品国产色婷婷| 在线 av 中文字幕| 真实男女啪啪啪动态图| 日本熟妇午夜| 日本与韩国留学比较| 偷拍熟女少妇极品色| 亚洲成人一二三区av| 看免费成人av毛片| 亚洲精品久久久久久婷婷小说| 免费播放大片免费观看视频在线观看| 七月丁香在线播放| 爱豆传媒免费全集在线观看| 中国三级夫妇交换| 亚洲自偷自拍三级| 简卡轻食公司| 在线精品无人区一区二区三 | 免费观看a级毛片全部| 久久久国产一区二区| 午夜激情福利司机影院| 欧美+日韩+精品| 一级毛片 在线播放| 久久久久网色| 国产精品麻豆人妻色哟哟久久| 99热这里只有是精品50| 久久久久久久亚洲中文字幕| 天堂网av新在线| 日本-黄色视频高清免费观看| 日韩不卡一区二区三区视频在线| 久久ye,这里只有精品| 丰满少妇做爰视频| 丝袜喷水一区| 欧美xxⅹ黑人| 少妇的逼好多水| 插阴视频在线观看视频| 久久久久久久久久久丰满| 91精品伊人久久大香线蕉| 五月玫瑰六月丁香| 我的女老师完整版在线观看| 99热网站在线观看| 亚洲av福利一区| 亚洲av福利一区| 欧美日韩国产mv在线观看视频 | 国产精品国产三级国产专区5o| 国产探花在线观看一区二区| 亚洲成色77777| 成人国产av品久久久| 看十八女毛片水多多多| 黄色怎么调成土黄色| 在线观看一区二区三区激情| 亚洲国产精品专区欧美| 乱系列少妇在线播放| 你懂的网址亚洲精品在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品不卡视频一区二区| 日韩强制内射视频| 人妻 亚洲 视频| 大码成人一级视频| 真实男女啪啪啪动态图| 免费观看a级毛片全部| 亚洲精品成人av观看孕妇| 在线看a的网站| 性插视频无遮挡在线免费观看| 精品少妇黑人巨大在线播放| 美女xxoo啪啪120秒动态图| 成人欧美大片| 亚洲国产欧美在线一区| 乱码一卡2卡4卡精品| 精品国产三级普通话版| 高清毛片免费看| 少妇人妻一区二区三区视频| 国产精品一二三区在线看| 交换朋友夫妻互换小说| 国产淫片久久久久久久久| av.在线天堂| 六月丁香七月| 网址你懂的国产日韩在线| 久久久久久久精品精品| 99re6热这里在线精品视频| 搡老乐熟女国产| 国产伦精品一区二区三区四那| 亚洲国产欧美在线一区| 欧美 日韩 精品 国产| 精品熟女少妇av免费看| 青青草视频在线视频观看| 99热这里只有是精品50| 久久国产乱子免费精品| 亚洲欧美日韩另类电影网站 | 少妇人妻一区二区三区视频| 热re99久久精品国产66热6| 国产精品国产三级专区第一集| 国产淫片久久久久久久久| 精品国产露脸久久av麻豆| 亚洲国产精品成人综合色| 天天一区二区日本电影三级| 午夜免费鲁丝| 亚洲,欧美,日韩| 国产视频首页在线观看| 可以在线观看毛片的网站| 国产欧美另类精品又又久久亚洲欧美| 国产大屁股一区二区在线视频| 久久精品久久久久久噜噜老黄| 女人被狂操c到高潮| 日日啪夜夜撸| 秋霞在线观看毛片| av在线观看视频网站免费| 亚洲欧美一区二区三区国产| 亚洲精品自拍成人| 亚洲精品第二区| 久久精品人妻少妇| 麻豆乱淫一区二区| 汤姆久久久久久久影院中文字幕| 国产中年淑女户外野战色| 亚洲av欧美aⅴ国产| 春色校园在线视频观看| 肉色欧美久久久久久久蜜桃 | 日韩一本色道免费dvd| 一本一本综合久久| 午夜福利在线观看免费完整高清在| 国产免费一区二区三区四区乱码| 九草在线视频观看| 亚洲欧美日韩卡通动漫| 亚洲av二区三区四区| 校园人妻丝袜中文字幕| 人妻制服诱惑在线中文字幕| 91精品伊人久久大香线蕉| 亚洲色图av天堂| 啦啦啦中文免费视频观看日本| 麻豆成人av视频| av在线天堂中文字幕| 欧美高清成人免费视频www| 亚洲欧洲日产国产| 真实男女啪啪啪动态图| 欧美日韩在线观看h| 2022亚洲国产成人精品| 亚洲精品自拍成人| 天美传媒精品一区二区| 亚洲精品国产av蜜桃| 国产成年人精品一区二区| 久久午夜福利片| 80岁老熟妇乱子伦牲交| 亚洲精品,欧美精品| 久久精品久久精品一区二区三区| 亚洲图色成人| 在线天堂最新版资源| 国产国拍精品亚洲av在线观看| 国产精品熟女久久久久浪| 大又大粗又爽又黄少妇毛片口| h日本视频在线播放| 午夜激情久久久久久久| 男人狂女人下面高潮的视频| 色网站视频免费| 观看美女的网站| 国产淫语在线视频| 汤姆久久久久久久影院中文字幕| 亚洲欧美清纯卡通| 麻豆成人av视频| 亚洲av电影在线观看一区二区三区 | 亚洲成色77777| 亚洲精品影视一区二区三区av| .国产精品久久| 一区二区三区免费毛片| 国产精品久久久久久久电影| 亚洲四区av| 成年女人看的毛片在线观看| 人妻少妇偷人精品九色| 91久久精品电影网| 久久久久国产精品人妻一区二区| av播播在线观看一区| 麻豆精品久久久久久蜜桃| 国产精品无大码| 内地一区二区视频在线| 日韩大片免费观看网站| 成人一区二区视频在线观看| 欧美bdsm另类| 日韩一区二区三区影片| 内射极品少妇av片p| 亚洲国产精品专区欧美| 一二三四中文在线观看免费高清| 亚洲精品国产av蜜桃| 黄色一级大片看看| 少妇人妻久久综合中文| 久久精品国产亚洲网站| 国产精品国产av在线观看| 久久99热这里只有精品18| 成人无遮挡网站| 久久久久久久久久成人| 高清av免费在线| 乱码一卡2卡4卡精品| 黄片wwwwww| 99热这里只有是精品在线观看| 纵有疾风起免费观看全集完整版| 亚洲国产最新在线播放| 成人高潮视频无遮挡免费网站| 汤姆久久久久久久影院中文字幕| 久久久久久九九精品二区国产| 国产成人精品婷婷| 亚洲精品中文字幕在线视频 | 亚洲精品乱码久久久久久按摩| 69av精品久久久久久| 香蕉精品网在线| 精品国产露脸久久av麻豆| 一级片'在线观看视频| 国产亚洲一区二区精品| 亚洲欧美一区二区三区黑人 | 女人十人毛片免费观看3o分钟| 欧美97在线视频| 精品久久久久久电影网| 久久精品久久久久久久性| 亚洲电影在线观看av| 美女xxoo啪啪120秒动态图| 麻豆精品久久久久久蜜桃| 看十八女毛片水多多多| 亚洲自偷自拍三级| 国产精品久久久久久精品古装| 全区人妻精品视频| 搡老乐熟女国产| videos熟女内射| 精品午夜福利在线看| 国产伦精品一区二区三区四那| 中国国产av一级| 精品视频人人做人人爽| 成人国产麻豆网| 成年免费大片在线观看| 黄色配什么色好看| 七月丁香在线播放| 午夜老司机福利剧场| 久久精品国产鲁丝片午夜精品| 久久久久久九九精品二区国产| av在线播放精品| 中文字幕久久专区| 亚洲国产av新网站| 国产精品国产三级国产专区5o| 天天一区二区日本电影三级| 国产伦精品一区二区三区视频9| av国产久精品久网站免费入址| 成人美女网站在线观看视频| 91aial.com中文字幕在线观看| 亚洲精品日本国产第一区| 岛国毛片在线播放| 亚洲精品国产色婷婷电影| 国产淫语在线视频| 亚洲欧美成人综合另类久久久| 51国产日韩欧美| 夫妻性生交免费视频一级片| 国产视频首页在线观看| 欧美成人午夜免费资源| 3wmmmm亚洲av在线观看| 国产乱来视频区| 寂寞人妻少妇视频99o| 插阴视频在线观看视频| 欧美三级亚洲精品| av卡一久久| 午夜福利网站1000一区二区三区| 婷婷色av中文字幕| 少妇被粗大猛烈的视频| 18禁动态无遮挡网站| 亚洲自拍偷在线| 热99国产精品久久久久久7| 一个人看视频在线观看www免费| 最近最新中文字幕免费大全7| 国产精品国产三级国产专区5o| 中文天堂在线官网| 王馨瑶露胸无遮挡在线观看| 秋霞伦理黄片| 97在线视频观看| 亚洲精品一二三| 偷拍熟女少妇极品色| 国产伦在线观看视频一区| 色哟哟·www| 国产在线一区二区三区精| av免费观看日本| 免费在线观看成人毛片| 成年女人看的毛片在线观看| 毛片一级片免费看久久久久| 91狼人影院| 深爱激情五月婷婷| 久久亚洲国产成人精品v| 综合色av麻豆| 欧美一级a爱片免费观看看| 久久久久久伊人网av| 老女人水多毛片| 五月玫瑰六月丁香| 日韩国内少妇激情av| 色婷婷久久久亚洲欧美| 黄片wwwwww| 在线亚洲精品国产二区图片欧美 | 国产乱人偷精品视频| 少妇高潮的动态图| 免费大片18禁| 91狼人影院| 超碰av人人做人人爽久久| 麻豆成人av视频| 嫩草影院精品99| 十八禁网站网址无遮挡 | 久久人人爽人人爽人人片va| 大香蕉久久网| 久久久久国产精品人妻一区二区| 国产精品久久久久久精品电影| 又爽又黄a免费视频| 各种免费的搞黄视频| 国产高清三级在线| 高清日韩中文字幕在线| 亚洲精品乱码久久久v下载方式| 日本色播在线视频| 成人鲁丝片一二三区免费| 日产精品乱码卡一卡2卡三| 成年免费大片在线观看| 国产午夜精品久久久久久一区二区三区| 欧美xxxx性猛交bbbb| 亚洲精品乱码久久久久久按摩| 麻豆乱淫一区二区| 美女国产视频在线观看| 久久久久久久久久成人| 久久久久精品性色| 男女下面进入的视频免费午夜| 久久久久久久亚洲中文字幕| 在线观看国产h片| 精品国产一区二区三区久久久樱花 | 大香蕉久久网| 国产老妇伦熟女老妇高清| 蜜臀久久99精品久久宅男| 极品教师在线视频| 免费观看a级毛片全部| 国产爽快片一区二区三区| 99久久精品热视频| av黄色大香蕉| 黑人高潮一二区| 又爽又黄a免费视频| 少妇人妻一区二区三区视频| 国产精品蜜桃在线观看| 日韩中字成人| 一级av片app| 亚洲怡红院男人天堂| tube8黄色片| 亚洲国产精品成人综合色| 汤姆久久久久久久影院中文字幕| 99热全是精品| 国产亚洲91精品色在线| 欧美最新免费一区二区三区| 在现免费观看毛片| 亚洲欧美成人综合另类久久久| 亚洲国产精品999| 午夜福利视频精品| 久久99蜜桃精品久久| 亚洲丝袜综合中文字幕| 观看美女的网站| 22中文网久久字幕| 女人被狂操c到高潮| 97人妻精品一区二区三区麻豆| 高清在线视频一区二区三区| 日韩不卡一区二区三区视频在线|