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

    板結(jié)構(gòu)多模態(tài)動力吸振設(shè)計方法

    2024-02-29 01:45:52鄧子偉孫啟航張保成
    中國海洋平臺 2024年1期
    關(guān)鍵詞:模態(tài)振動

    彭 磊, 鄧子偉, 劉 鵬, 孫啟航, 張保成

    (中國海洋大學(xué) 工程學(xué)院 機(jī)電工程系, 山東 青島 266404)

    0 引 言

    海洋平臺的振動噪聲會影響平臺的安全性和工作人員的舒適性,對平臺的正常作業(yè)帶來多方面的挑戰(zhàn),如何有效實現(xiàn)平臺的減振降噪是當(dāng)前專家學(xué)者廣泛關(guān)注的問題[1]。海洋平臺是由船體外板、甲板、艙壁等組成的復(fù)雜板架結(jié)構(gòu),進(jìn)行板結(jié)構(gòu)的振動控制研究,有助于化繁為簡地探索平臺復(fù)雜多模態(tài)的振動演化機(jī)制。

    針對平臺艙室、甲板等板結(jié)構(gòu)常用的減振方法有加筋加肋[2]、敷設(shè)阻尼材料[3]、安裝動力吸振器[4]等。加筋加肋可改善結(jié)構(gòu)頻率達(dá)到頻率避讓的目的,但改善結(jié)構(gòu)的處理方式僅能調(diào)節(jié)結(jié)構(gòu)的頻率,結(jié)構(gòu)的共振特性仍然存在,對于寬頻能量引起的振動存在局限性。敷設(shè)黏彈性阻尼材料可耗散振動能量、明顯降低結(jié)構(gòu)共振峰的響應(yīng),但阻尼材料的減振性能與模態(tài)應(yīng)變能相關(guān),對低階模態(tài)振動的抑制效果并不理想[5]。動力吸振器具有結(jié)構(gòu)簡單、可靠性強(qiáng)、無需外界能源等優(yōu)點,應(yīng)用更為廣泛。ZHU等[6]研究抑制薄板振動的分布式動力吸振器的參數(shù)優(yōu)化問題,推導(dǎo)得到分布式動力吸振器的參數(shù)優(yōu)化公式,在此基礎(chǔ)上采用分布式吸振器對薄板前三階模態(tài)進(jìn)行振動控制。王克肖等[7]推導(dǎo)附加動力吸振器的薄板結(jié)構(gòu)聲輻射效率計算方法,研究動力吸振器對薄板聲輻射效率的調(diào)控作用。高聰?shù)萚8]針對海洋平臺上層建筑低頻線譜的振動控制難題,以海洋平臺板架結(jié)構(gòu)為例研究分布式動力吸振器的布放位置、數(shù)量對減振效果的影響規(guī)律。然而海洋平臺中的板結(jié)構(gòu)大多為對稱結(jié)構(gòu),存在頻率相同或相近的密集模態(tài),目前對密集模態(tài)主導(dǎo)的多個模態(tài)耦合振動的吸振研究較少,關(guān)于動力吸振器設(shè)計方法的系統(tǒng)性研究還不夠深入。

    本文基于上述研究中存在的不足,建立薄板與動力吸振器的耦合模型,從理論上研究動力吸振器的質(zhì)量比、阻尼比、頻率比和安裝位置等4個參數(shù)的設(shè)計問題。在此基礎(chǔ)上,形成一套用于板結(jié)構(gòu)多模態(tài)減振的動力吸振器設(shè)計方法,將該方法應(yīng)用于艙室模型的振動控制,獲得的仿真結(jié)果可驗證平臺板結(jié)構(gòu)動力吸振方法的有效性。

    1 理論分析

    1.1 分析模型

    薄板與單個動力吸振器耦合的系統(tǒng)模型如圖1所示。

    圖1 安裝動力吸振器的薄板模型

    薄板受到簡諧激勵f(t),動力吸振器的質(zhì)量、阻尼和剛度分別為m1、c1和k1。根據(jù)結(jié)構(gòu)動力學(xué),該系統(tǒng)的運動方程為

    (1)

    式中:M、C、K分別為系統(tǒng)的質(zhì)量矩陣、阻尼矩陣和剛度矩陣;x(t)為系統(tǒng)的振動位移。激勵力向量F(t)的表達(dá)式為

    (2)

    式中:Φ=[φ1(x,y),φ2(x,y),φ3(x,y),…,φJ(rèn)(x,y)]T為薄板的模態(tài)振型向量,其中J為薄板參與振動的模態(tài)數(shù)量。

    當(dāng)一個結(jié)構(gòu)的振動模態(tài)在目標(biāo)頻率范圍內(nèi)可很好地分離時,振動響應(yīng)通常由單個模態(tài)控制[9]。假設(shè)薄板的振動響應(yīng)由其第i階模態(tài)控制,這時系統(tǒng)的運動方程可簡化為

    (3)

    式中:Mi、Ci、Ki分別為薄板第i階模態(tài)的模態(tài)質(zhì)量、模態(tài)阻尼和模態(tài)剛度;m1、c1、k1分別為動力吸振器的質(zhì)量、阻尼和剛度;x1為動力吸振器的位移;φi為薄板第i階模態(tài)的振型;φil為薄板第i階模態(tài)振型在動力吸振器安裝位置處的幅值。

    將式(3)進(jìn)行拉普拉斯變換并求解得到第i階模態(tài)位移響應(yīng)的復(fù)振幅為

    (4)

    式中:Xi(s)為xi(t)的拉普拉斯變換;Fi,g(s)=φiF(s)為廣義力在第i階模態(tài)的分量,其中F(s)為f(t)的拉普拉斯變換。

    令s=jω,其中ω為激勵力頻率。整理可以得到位移頻響函數(shù)為

    (5)

    對式(5)求導(dǎo)可得到速度響應(yīng)和加速度響應(yīng)的頻響函數(shù)為

    (6)

    為便于研究,引入5個無量綱參數(shù):

    (7)

    式中:μ為動力吸振器質(zhì)量與薄板第i階模態(tài)的模態(tài)質(zhì)量之比;f為動力吸振器頻率與薄板第i階模態(tài)固有頻率的比值;g為激勵力頻率與薄板第i階模態(tài)固有頻率的比值;ζi為薄板第i階模態(tài)的阻尼比;ζk為動力吸振器的阻尼比;ω1為動力吸振器的固有頻率;ωi為薄板第i階模態(tài)的固有頻率。

    將式(7)代入式(5)和式(6),位移、速度和加速度分別用下角標(biāo)di、ve和ac區(qū)分,得到第i階模態(tài)3個響應(yīng)量的歸一化頻響函數(shù)分別為

    (8)

    式中:

    (9)

    1.2 動力吸振器最優(yōu)參數(shù)

    動力吸振器的減振性能與其參數(shù)密切相關(guān),安裝在連續(xù)結(jié)構(gòu)上的動力吸振器共有4個參數(shù)需要設(shè)計,分別是動力吸振器的質(zhì)量、阻尼、剛度和安裝位置。其中,質(zhì)量根據(jù)主系統(tǒng)的質(zhì)量來選定,通常指定為與主系統(tǒng)質(zhì)量的比值,表示為質(zhì)量比μ。質(zhì)量確定后,剛度可由動力吸振器的固有頻率得到,頻率指定為與主系統(tǒng)固有頻率的比值,表示為頻率比f。阻尼則以臨界阻尼比ζk表示。質(zhì)量比、阻尼比和頻率比的表達(dá)式在式(7)中已給出。

    (1)動力吸振器的安裝位置。為便于討論,忽略薄板的阻尼。未安裝動力吸振器時,薄板的位移頻響函數(shù)為

    (10)

    式中:KP和MP分別為薄板的剛度矩陣和質(zhì)量矩陣,表達(dá)式分別為

    (11)

    根據(jù)式(11),式(10)還可寫為

    (12)

    安裝動力吸振器后,薄板的位移頻響函數(shù)為

    (13)

    式中:m1、k1、c1分別為動力吸振器的質(zhì)量、剛度和阻尼;(xa,ya)為動力吸振器安裝位置的坐標(biāo)。

    根據(jù)式(10)和式(12),式(13)可重新寫為

    (14)

    式中:ωa為動力吸振器的頻率。

    將薄板第i階模態(tài)作為目標(biāo)模態(tài),假設(shè)ωi-1、ωi和ωi+1為相鄰模態(tài)的固有頻率。若ωi-1、ωi和ωi+1的值相差較大,當(dāng)激勵力頻率ω在ωi附近時,有:

    (15)

    保持其他參數(shù)不變,薄板的振動響應(yīng)由第i階模態(tài)主導(dǎo),目標(biāo)頻率處的振動不受其他模態(tài)的影響,這時第i階模態(tài)為離散模態(tài)。對于離散模態(tài),動力吸振器安裝在目標(biāo)模態(tài)振型的幅值最大點效果最佳[10]。

    若ωi-1、ωi和ωi+1的值相近,則3個模態(tài)為密集模態(tài),當(dāng)激勵力頻率ω在ωi附近時,有:

    (16)

    這時薄板的振動響應(yīng)由3個密集模態(tài)共同主導(dǎo),僅針對第i階模態(tài)設(shè)置吸振器,吸振器也會與其他兩階模態(tài)耦合振動,又由于振型的不同,動力吸振器的參數(shù)對于其他兩階模態(tài)不是最優(yōu)參數(shù),會產(chǎn)生不利振動,導(dǎo)致達(dá)不到理想的減振效果。

    當(dāng)薄板包含阻尼時,安裝動力吸振器的薄板的位移頻響函數(shù)為

    (17)

    這時薄板的模態(tài)為復(fù)模態(tài),阻尼的引入并不會影響薄板的固有頻率,可用極點的虛部近似固有頻率,上述結(jié)論仍然適用于阻尼薄板[11]。

    (2)質(zhì)量比的選定。已有研究[12]表明質(zhì)量比越大,動力吸振器的減振性能越強(qiáng)。但質(zhì)量比并非越大越好,隨著質(zhì)量比的增大,主系統(tǒng)響應(yīng)降低的幅度有限,在設(shè)計質(zhì)量比時需要綜合考慮吸振性能與經(jīng)濟(jì)性。

    (3)頻率比和阻尼比。這2個參數(shù)可通過優(yōu)化設(shè)計公式得到,位移和加速度是平臺減振常用的目標(biāo)響應(yīng)量,吸振的目標(biāo)函數(shù)為

    (18)

    當(dāng)位移為目標(biāo)響應(yīng)量時,最優(yōu)頻率比和最優(yōu)阻尼比的優(yōu)化公式[13]分別為

    (19)

    (20)

    當(dāng)加速度為目標(biāo)響應(yīng)量時,最優(yōu)頻率比和最優(yōu)阻尼比的優(yōu)化公式[13]分別為

    (21)

    (22)

    2 板結(jié)構(gòu)動力吸振方法

    根據(jù)第1節(jié)的研究對板結(jié)構(gòu)動力吸振方法的流程進(jìn)行總結(jié),該方法可分為振動特性分析、受控模態(tài)與目標(biāo)響應(yīng)量確定、動力吸振器的安裝位置和數(shù)量確定、動力吸振器最優(yōu)參數(shù)確定、吸振效果驗證以及補(bǔ)正設(shè)計共6個步驟。動力吸振器設(shè)計方法的流程圖如圖2所示。

    圖2 板結(jié)構(gòu)動力吸振方法流程圖

    各個步驟的說明如下:

    (1)振動特性分析。振動特性分析是對板結(jié)構(gòu)的自由振動和受迫振動進(jìn)行分析。自由振動分析是為了得到板結(jié)構(gòu)的固有頻率、模態(tài)振型和模態(tài)質(zhì)量等模態(tài)參數(shù)。受迫振動分析是根據(jù)板結(jié)構(gòu)實際工況施加激勵力,確定考核區(qū)域得到板結(jié)構(gòu)的頻響曲線,獲得結(jié)構(gòu)的動態(tài)振動特性。

    (2)受控模態(tài)與目標(biāo)響應(yīng)量確定。在受迫振動特性分析的基礎(chǔ)上,根據(jù)實際振動控制要求,確定待控制的目標(biāo)模態(tài)和響應(yīng)量。響應(yīng)量包括位移、速度和加速度等3種。平臺通常選擇位移和加速度作為目標(biāo)響應(yīng)量。其中:位移響應(yīng)會導(dǎo)致結(jié)構(gòu)的疲勞、變形、碰撞或斷裂;加速度響應(yīng)則與人體舒適度相關(guān)。

    (3)動力吸振器的安裝位置和數(shù)量確定。對于離散模態(tài),應(yīng)選擇目標(biāo)模態(tài)的振型幅值最大點作為動力吸振器的安裝位置。對于密集模態(tài),吸振器的最佳安裝位置應(yīng)使目標(biāo)模態(tài)的振型幅值盡可能大,使其他頻率相近模態(tài)的振型幅值盡可能小。動力吸振器的數(shù)量由其安裝位置決定。

    (4)動力吸振器最優(yōu)參數(shù)確定。動力吸振器的安裝位置確定后,還差動力吸振器的質(zhì)量、頻率和阻尼即可得到所需動力吸振器的全部參數(shù)。吸振器的常用質(zhì)量比范圍為0~0.2[14],頻率比和阻尼比可通過優(yōu)化設(shè)計公式計算得到。

    (5)吸振效果驗證。設(shè)計完動力吸振器后需要通過數(shù)值計算或有限元仿真對吸振器的減振結(jié)果進(jìn)行檢驗,若吸振效果不理想則需要后續(xù)的補(bǔ)正設(shè)計。

    (6)補(bǔ)正設(shè)計。補(bǔ)正設(shè)計是在現(xiàn)有動力吸振器的基礎(chǔ)上對其參數(shù)進(jìn)行調(diào)整,先對比附加吸振器后形成的2個共振峰的相對偏差:若吸振器2個共振峰相對偏差較小,需要增大吸振器的質(zhì)量;若吸振器2個共振峰相對偏差較大,對動力吸振器的頻率比和阻尼比進(jìn)行調(diào)整可提升減振效果。頻率比的增大會使2個共振峰中低階的共振峰增大,使高階的共振峰減小;頻率比的減小會使高階的共振峰增大,使低階的共振峰減小;阻尼比的增大會降低2個共振峰的峰值,但會降低2個共振峰之間頻帶的減振效果[15]??筛鶕?jù)此規(guī)律對2個共振峰進(jìn)行調(diào)整。

    3 艙室動力吸振應(yīng)用

    第2節(jié)對板結(jié)構(gòu)中動力吸振器的設(shè)計流程進(jìn)行總結(jié),形成了板結(jié)構(gòu)動力吸振方法。因此,基于該方法,以一個艙室模型為具體應(yīng)用對象,通過對艙室的艙壁進(jìn)行振動控制來驗證其有效性。

    3.1 艙室模型

    艙室的結(jié)構(gòu)形式參考文獻(xiàn)[3],艙室的整體尺寸為長3.0 m、寬2.7 m、高2.5 m。艙壁的厚度為0.006 m。在艙室的6個艙壁內(nèi)面采用加強(qiáng)肋結(jié)構(gòu)。加強(qiáng)肋的高度為0.1 m,其截面為0.1 m×0.01 m的矩形。加強(qiáng)肋的尺寸和分布如圖3所示。

    圖3 艙室內(nèi)部加強(qiáng)肋的尺寸和分布

    按照上述尺寸進(jìn)行建模,艙室模型如圖4所示。艙室材料均選用結(jié)構(gòu)鋼,彈性模量E=2×1011Pa,密度ρ=7 850 kg/m3,泊松比μ=0.3。艙室的阻尼比設(shè)置為0.01。

    圖4 艙室模型

    3.2 艙室振動特性分析

    對艙室的自由振動進(jìn)行分析。單元尺寸設(shè)置為0.1 m,艙室有限元模型共包含單元50 610個。艙室底部設(shè)置固定支撐。分析得到艙室前6階模態(tài)的固有頻率和振型如圖5所示。

    圖5 艙室前6階自由振動模態(tài)及振型

    將艙室頂部的壁板作為振動考察區(qū)域,對頂部艙壁的受迫振動進(jìn)行分析。在頂部艙壁的法線方向施加大小為1 000 N的簡諧激勵力。由于艙室前6階模態(tài)的頻率都在50 Hz以內(nèi),因此在0~50 Hz范圍內(nèi)對激勵力進(jìn)行掃頻,掃頻間隔為0.2 Hz。激勵點和考核點位置如圖6所示。

    圖6 激勵點和考核點在頂部艙壁的位置

    考核點的頻響曲線如圖7所示,共有5個共振峰,即有5階模態(tài)參與了振動。根據(jù)自由振動的分析結(jié)果,這5階模態(tài)分別為第一、二、四、五、六階模態(tài)。由于考核點和激勵點在第三階模態(tài)振型的節(jié)線上,因此沒有激發(fā)第三階模態(tài)的振動。在參與振動的5階模態(tài)中,第一、五、六階模態(tài)的共振峰相距較遠(yuǎn),可視作離散模態(tài)。第二、四階頻率相近,是密集模態(tài)。下文根據(jù)板結(jié)構(gòu)動力吸振方法設(shè)計動力吸振器對參與振動的5階模態(tài)進(jìn)行振動控制。

    圖7 艙壁考核點頻響曲線

    3.3 動力吸振器設(shè)計及吸振效果驗證

    先確定吸振器的安裝位置。第一、五、六階模態(tài)為離散模態(tài),因此吸振器安裝在各自模態(tài)的振型幅值最大點。第二、四階模態(tài)為密集模態(tài),針對其中一階模態(tài)的吸振器的安裝位置應(yīng)盡可能避開另一階模態(tài)的影響。綜合分析設(shè)置吸振的安裝位置如圖8所示。針對第一、二、四、五、六階模態(tài)的動力吸振器分別為DVA1、DVA2、DVA4、DVA5和DVA6。

    圖8 動力吸振器安裝位置

    安裝位置確定后,計算動力吸振器的最優(yōu)參數(shù)還需要目標(biāo)模態(tài)的模態(tài)質(zhì)量和動力吸振器安裝位置處目標(biāo)模態(tài)振型的幅值大小。模態(tài)質(zhì)量和振型幅值可在自由振動分析中得到。吸振器的質(zhì)量、數(shù)量以及根據(jù)參數(shù)設(shè)計公式得到的吸振器剛度和阻尼如表1所示。

    表1 動力吸振器最優(yōu)參數(shù)

    安裝動力吸振器前后,艙壁振動的頻率響應(yīng)對比曲線如圖9所示。第一階模態(tài)的位移響應(yīng)降低90.6%(從11.24 mm至1.05 mm),加速度響應(yīng)降低92.2%(從214.70 m/s2至16.62 m/s2)。第二階模態(tài)的位移響應(yīng)降低35.9%(從1.03 mm至0.66 mm),加速度響應(yīng)降低32.4%(從31.73 m/s2至21.44 m/s2)。第四階模態(tài)的位移響應(yīng)降低61.2%(從1.70 mm至0.66 mm),加速度響應(yīng)降低63.4%(從58.63 m/s2至21.44 m/s2)。第五階模態(tài)的位移響應(yīng)降低51.5%(從0.66 mm至0.32 mm),加速度響應(yīng)降低48.3%(從34.71 m/s2至17.95 m/s2)。第六階模態(tài)的位移響應(yīng)降低89.7%(從5.61 mm至0.58 mm),加速度響應(yīng)降低88.9%(從488.90 m/s2至54.25 m/s2)。由結(jié)果可知,采用本方法設(shè)計的動力吸振器取得了理想的減振效果,有效控制艙壁多個模態(tài)的振動。尤其是第二、四階模態(tài)所處的頻段,動力吸振器在抑制2個共振峰振動的同時并沒有引發(fā)其他不良振動,這驗證了所總結(jié)的板結(jié)構(gòu)動力吸振方法的有效性。

    圖9 艙壁減振前后頻率響應(yīng)對比曲線

    3.4 補(bǔ)正設(shè)計

    在第3.3節(jié)中,所設(shè)計的動力吸振器取得了良好的減振效果,如果想要提升減振效果,可對動力吸振器進(jìn)行補(bǔ)正設(shè)計,通過調(diào)整動力吸振器的參數(shù)來提升其吸振性能。以第一階和第五階模態(tài)的位移響應(yīng)為例進(jìn)行說明。

    由圖9(a)可知,第一階和第五階模態(tài)新產(chǎn)生的共振峰并非等高,第一階模態(tài)的左側(cè)共振峰偏高,第五階模態(tài)的右側(cè)共振峰偏高。這是由于針對第一階和第五階模態(tài)的動力吸振器安裝位置相同,且安裝位置都為各自模態(tài)振型的幅值最大點。盡管兩者為離散模態(tài),但動力吸振器相對非目標(biāo)模態(tài)相當(dāng)于一個附加質(zhì)量,會影響艙室的模態(tài),從而導(dǎo)致2個共振峰的高度不同。此外,第五階模態(tài)參與振動的程度小,其共振峰較低,右側(cè)非共振峰處的振動也會對其產(chǎn)生影響,這也是其右側(cè)共振峰偏高的原因。

    根據(jù)補(bǔ)正設(shè)計方法,調(diào)節(jié)頻率比可調(diào)整2個共振峰的相對大?。簩τ诘谝浑A模態(tài),動力吸振器的頻率比應(yīng)減??;對于第五階模態(tài),動力吸振器的頻率比應(yīng)增大。補(bǔ)正后動力吸振器的參數(shù)如表2所示。

    表2 補(bǔ)正后的動力吸振器參數(shù)

    補(bǔ)正前后減振效果的對比如圖10所示。第一階和第五階模態(tài)的位移響應(yīng)分別降低91.5%(從11.24 mm至0.95 mm)和57.6%(從0.66 mm至0.28 mm)。對比補(bǔ)正前的動力吸振器,減振效果分別提升0.9個百分點和6.1個百分點,由此可見,通過補(bǔ)正設(shè)計可在一定程度上提高動力吸振器的減振效果。

    圖10 補(bǔ)正設(shè)計前后減振效果對比

    4 結(jié) 論

    根據(jù)板結(jié)構(gòu)的頻響公式從理論上探討多個模態(tài)耦合振動時動力吸振器的最佳安裝位置,分析動力吸振器質(zhì)量比的設(shè)計問題,給出動力吸振器阻尼比和頻率比的優(yōu)化設(shè)計公式,形成一套用于板結(jié)構(gòu)減振的動力吸振器的設(shè)計方法,該方法共有6個步驟,分別是:振動特性分析、受控模態(tài)與目標(biāo)響應(yīng)確定、動力吸振器的安裝位置和數(shù)量確定、吸振器最優(yōu)參數(shù)確定、吸振效果驗證、補(bǔ)正設(shè)計。將該方法應(yīng)用于艙室模型的振動控制,相關(guān)結(jié)論如下:

    (1)動力吸振器對密集模態(tài)的振動抑制效果差是因為這時振動響應(yīng)由多個模態(tài)共同主導(dǎo)。吸振器的最佳安裝位置應(yīng)使目標(biāo)模態(tài)振型的幅值盡可能大,使其他模態(tài)振型的幅值盡可能小。

    (2)在艙壁的減振控制中,第一階模態(tài)減振效果最佳,其位移和加速度響應(yīng)分別降低90.6%和92.2%。這驗證了所總結(jié)的板結(jié)構(gòu)動力吸振方法的有效性。

    (3)對動力吸振器進(jìn)行補(bǔ)正設(shè)計后,針對艙壁第一階和第五階模態(tài)位移響應(yīng)的減振效果分別提升0.9個百分點和6.1個百分點,補(bǔ)正設(shè)計可提升動力吸振器的性能。

    猜你喜歡
    模態(tài)振動
    振動的思考
    噴水推進(jìn)高速艇尾部振動響應(yīng)分析
    This “Singing Highway”plays music
    振動攪拌 震動創(chuàng)新
    中國公路(2017年18期)2018-01-23 03:00:38
    中立型Emden-Fowler微分方程的振動性
    車輛CAE分析中自由模態(tài)和約束模態(tài)的應(yīng)用與對比
    國內(nèi)多模態(tài)教學(xué)研究回顧與展望
    高速顫振模型設(shè)計中顫振主要模態(tài)的判斷
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識別
    UF6振動激發(fā)態(tài)分子的振動-振動馳豫
    計算物理(2014年2期)2014-03-11 17:01:44
    亚洲欧美激情在线| 久久久精品94久久精品| 国产视频一区二区在线看| 欧美乱妇无乱码| 51午夜福利影视在线观看| 国产亚洲欧美在线一区二区| 18禁观看日本| 国产日韩欧美亚洲二区| a级毛片黄视频| 黄片播放在线免费| 午夜老司机福利片| 桃红色精品国产亚洲av| 极品教师在线免费播放| 久久精品熟女亚洲av麻豆精品| 欧美国产精品va在线观看不卡| 91大片在线观看| 亚洲精品国产一区二区精华液| 中文欧美无线码| 99在线人妻在线中文字幕 | 王馨瑶露胸无遮挡在线观看| 久久精品国产99精品国产亚洲性色 | 亚洲欧美日韩高清在线视频 | 亚洲一区中文字幕在线| 正在播放国产对白刺激| 国产97色在线日韩免费| 少妇粗大呻吟视频| 一本久久精品| 中文字幕色久视频| 免费观看av网站的网址| 在线永久观看黄色视频| 妹子高潮喷水视频| 王馨瑶露胸无遮挡在线观看| 黄色丝袜av网址大全| 美女主播在线视频| av超薄肉色丝袜交足视频| 精品人妻在线不人妻| 中文字幕制服av| 手机成人av网站| 午夜福利视频精品| 夫妻午夜视频| 久久久精品国产亚洲av高清涩受| 欧美日韩成人在线一区二区| 久久人妻熟女aⅴ| 在线天堂中文资源库| 啦啦啦免费观看视频1| 日韩制服丝袜自拍偷拍| 日韩中文字幕欧美一区二区| 国精品久久久久久国模美| 免费av中文字幕在线| 成人18禁在线播放| 久久久精品94久久精品| 搡老乐熟女国产| 丰满迷人的少妇在线观看| 中文字幕最新亚洲高清| 啦啦啦中文免费视频观看日本| 国产精品一区二区在线不卡| 国产精品欧美亚洲77777| 婷婷丁香在线五月| 夜夜夜夜夜久久久久| 午夜精品久久久久久毛片777| 日韩大码丰满熟妇| 天天操日日干夜夜撸| 日韩视频一区二区在线观看| 捣出白浆h1v1| kizo精华| 91麻豆精品激情在线观看国产 | av网站在线播放免费| www日本在线高清视频| 久久ye,这里只有精品| 十分钟在线观看高清视频www| 99香蕉大伊视频| 久久性视频一级片| 视频在线观看一区二区三区| 国产成人精品无人区| 极品教师在线免费播放| 欧美精品人与动牲交sv欧美| 一夜夜www| 女性生殖器流出的白浆| www.999成人在线观看| 欧美黑人欧美精品刺激| 丁香六月天网| 亚洲 国产 在线| 少妇粗大呻吟视频| 亚洲中文字幕日韩| 久热爱精品视频在线9| 国产精品国产高清国产av | av又黄又爽大尺度在线免费看| 极品人妻少妇av视频| 久久精品成人免费网站| 黄频高清免费视频| 精品福利永久在线观看| 狠狠狠狠99中文字幕| 两性夫妻黄色片| 精品久久久久久电影网| 熟女少妇亚洲综合色aaa.| 国产成人av激情在线播放| 人妻久久中文字幕网| 亚洲成国产人片在线观看| 老司机午夜十八禁免费视频| 极品人妻少妇av视频| 国产精品一区二区免费欧美| 777米奇影视久久| 夜夜夜夜夜久久久久| 国产国语露脸激情在线看| 免费人妻精品一区二区三区视频| 欧美精品一区二区大全| 男女下面插进去视频免费观看| 国内毛片毛片毛片毛片毛片| 啦啦啦 在线观看视频| 免费在线观看影片大全网站| 侵犯人妻中文字幕一二三四区| 美女福利国产在线| 久久国产亚洲av麻豆专区| 人人澡人人妻人| 夫妻午夜视频| 国产亚洲午夜精品一区二区久久| 国产精品免费大片| 成人国语在线视频| 桃花免费在线播放| 国产麻豆69| 精品国产一区二区久久| 在线观看免费午夜福利视频| 免费人妻精品一区二区三区视频| 欧美成狂野欧美在线观看| 亚洲色图 男人天堂 中文字幕| 嫁个100分男人电影在线观看| 最近最新免费中文字幕在线| av网站免费在线观看视频| a级毛片在线看网站| 电影成人av| 男女边摸边吃奶| 亚洲人成电影观看| 亚洲精品自拍成人| 1024视频免费在线观看| 欧美亚洲 丝袜 人妻 在线| 中国美女看黄片| 国产精品偷伦视频观看了| 免费人妻精品一区二区三区视频| 成年人黄色毛片网站| 国产色视频综合| 中文欧美无线码| 国产熟女午夜一区二区三区| 国产欧美日韩综合在线一区二区| 国产一区二区在线观看av| 性少妇av在线| 欧美黄色淫秽网站| 成人三级做爰电影| √禁漫天堂资源中文www| 黄色丝袜av网址大全| 国产精品一区二区在线不卡| 欧美久久黑人一区二区| 久久久久久久久免费视频了| 黄色a级毛片大全视频| 午夜福利视频精品| 在线观看免费高清a一片| 欧美精品亚洲一区二区| 丝袜喷水一区| 十八禁网站免费在线| 狠狠狠狠99中文字幕| 欧美日韩视频精品一区| 国产激情久久老熟女| 99久久99久久久精品蜜桃| 欧美精品啪啪一区二区三区| 欧美精品亚洲一区二区| 亚洲中文字幕日韩| 午夜精品国产一区二区电影| 日韩视频在线欧美| av有码第一页| 19禁男女啪啪无遮挡网站| 老司机影院毛片| 俄罗斯特黄特色一大片| 欧美日韩亚洲综合一区二区三区_| 国内毛片毛片毛片毛片毛片| 99国产精品一区二区蜜桃av | 性色av乱码一区二区三区2| 久久久久久久精品吃奶| 欧美日韩成人在线一区二区| av福利片在线| 欧美黄色淫秽网站| 亚洲avbb在线观看| 国产成人av激情在线播放| 在线观看免费高清a一片| 欧美日韩国产mv在线观看视频| 美女国产高潮福利片在线看| 成人特级黄色片久久久久久久 | 99re在线观看精品视频| 成年动漫av网址| 美女高潮喷水抽搐中文字幕| 99在线人妻在线中文字幕 | 欧美黑人欧美精品刺激| 精品久久久久久久毛片微露脸| 日韩欧美免费精品| 久久免费观看电影| 欧美成狂野欧美在线观看| 无人区码免费观看不卡 | 精品人妻1区二区| 日本欧美视频一区| kizo精华| 国产一区二区三区在线臀色熟女 | 欧美日韩国产mv在线观看视频| 五月天丁香电影| 久久久久久免费高清国产稀缺| 丝瓜视频免费看黄片| 国产淫语在线视频| 精品乱码久久久久久99久播| 国产成人精品在线电影| 亚洲九九香蕉| a级毛片在线看网站| 亚洲精品一二三| 中文字幕人妻丝袜一区二区| 久久国产精品大桥未久av| 人成视频在线观看免费观看| 亚洲国产精品一区二区三区在线| 在线观看免费日韩欧美大片| 国产成人免费观看mmmm| 黑人猛操日本美女一级片| 一级毛片女人18水好多| 一边摸一边抽搐一进一出视频| 国产真人三级小视频在线观看| 亚洲精品中文字幕一二三四区 | 99久久人妻综合| www.999成人在线观看| 一区在线观看完整版| 两个人看的免费小视频| 激情视频va一区二区三区| 亚洲av美国av| 国产成人影院久久av| 欧美成人免费av一区二区三区 | 亚洲色图 男人天堂 中文字幕| 少妇的丰满在线观看| 精品国内亚洲2022精品成人 | 欧美精品一区二区大全| 99香蕉大伊视频| 欧美成人免费av一区二区三区 | 欧美日韩成人在线一区二区| 亚洲免费av在线视频| 精品第一国产精品| 欧美日韩国产mv在线观看视频| 亚洲精品一卡2卡三卡4卡5卡| 韩国精品一区二区三区| 欧美午夜高清在线| 一个人免费在线观看的高清视频| 欧美精品人与动牲交sv欧美| bbb黄色大片| av一本久久久久| 亚洲av成人一区二区三| 天堂俺去俺来也www色官网| 欧美av亚洲av综合av国产av| 18禁国产床啪视频网站| 精品国产乱码久久久久久男人| 欧美精品一区二区免费开放| 成人国语在线视频| 国产精品电影一区二区三区 | 咕卡用的链子| 欧美日韩黄片免| 在线观看免费视频网站a站| 国产免费现黄频在线看| 嫁个100分男人电影在线观看| 成年女人毛片免费观看观看9 | 天天操日日干夜夜撸| 精品一区二区三区四区五区乱码| 免费av中文字幕在线| 久久精品国产综合久久久| 国产无遮挡羞羞视频在线观看| 国产xxxxx性猛交| 女人久久www免费人成看片| 这个男人来自地球电影免费观看| 操美女的视频在线观看| 人人妻人人添人人爽欧美一区卜| 激情在线观看视频在线高清 | 少妇粗大呻吟视频| 在线永久观看黄色视频| 中文字幕人妻丝袜制服| 一进一出抽搐动态| 黑人欧美特级aaaaaa片| 女性被躁到高潮视频| 狠狠精品人妻久久久久久综合| 国产高清激情床上av| 99国产综合亚洲精品| 两性夫妻黄色片| 久久久精品94久久精品| 女人爽到高潮嗷嗷叫在线视频| 蜜桃国产av成人99| 久久狼人影院| 成在线人永久免费视频| 后天国语完整版免费观看| 一级片'在线观看视频| 咕卡用的链子| 亚洲专区中文字幕在线| 日本精品一区二区三区蜜桃| 亚洲av欧美aⅴ国产| 一级黄色大片毛片| 久久毛片免费看一区二区三区| tocl精华| 亚洲午夜理论影院| 国产精品麻豆人妻色哟哟久久| 一级,二级,三级黄色视频| 人人澡人人妻人| 97在线人人人人妻| cao死你这个sao货| 黑人欧美特级aaaaaa片| 伊人久久大香线蕉亚洲五| 51午夜福利影视在线观看| av天堂久久9| 日韩一区二区三区影片| 1024香蕉在线观看| 黄片大片在线免费观看| 9色porny在线观看| 女人被躁到高潮嗷嗷叫费观| 性高湖久久久久久久久免费观看| 日本撒尿小便嘘嘘汇集6| 亚洲专区字幕在线| 亚洲国产看品久久| 黄色视频不卡| 又大又爽又粗| 精品久久蜜臀av无| 天天躁日日躁夜夜躁夜夜| 欧美+亚洲+日韩+国产| 欧美日本中文国产一区发布| 两性夫妻黄色片| 久久国产精品人妻蜜桃| 国产男女内射视频| 国产av国产精品国产| 麻豆成人av在线观看| 久久精品国产亚洲av香蕉五月 | 国产午夜精品久久久久久| 亚洲精品中文字幕一二三四区 | 成人av一区二区三区在线看| 国产免费av片在线观看野外av| 欧美日韩亚洲综合一区二区三区_| 嫁个100分男人电影在线观看| 国产成人一区二区三区免费视频网站| 国产日韩一区二区三区精品不卡| 下体分泌物呈黄色| 亚洲国产看品久久| 黄色视频,在线免费观看| 国产熟女午夜一区二区三区| 淫妇啪啪啪对白视频| 久热爱精品视频在线9| 国产av一区二区精品久久| 色精品久久人妻99蜜桃| 18禁裸乳无遮挡动漫免费视频| 91九色精品人成在线观看| 亚洲av成人不卡在线观看播放网| 国产午夜精品久久久久久| 别揉我奶头~嗯~啊~动态视频| 人人澡人人妻人| 亚洲第一av免费看| 日本av免费视频播放| 久久人人爽av亚洲精品天堂| 亚洲五月色婷婷综合| 在线亚洲精品国产二区图片欧美| 人人妻人人爽人人添夜夜欢视频| 麻豆av在线久日| 国产精品亚洲av一区麻豆| 99久久99久久久精品蜜桃| 国产在线一区二区三区精| 91成年电影在线观看| bbb黄色大片| 99国产精品99久久久久| 好男人电影高清在线观看| 大香蕉久久网| 美女主播在线视频| 国产精品电影一区二区三区 | 老司机在亚洲福利影院| 精品一品国产午夜福利视频| 久久久久久久久久久久大奶| 亚洲avbb在线观看| 午夜激情av网站| 久久精品国产亚洲av香蕉五月 | 久久99一区二区三区| 亚洲综合色网址| 757午夜福利合集在线观看| 国产日韩欧美在线精品| 亚洲av日韩在线播放| 精品久久蜜臀av无| 色老头精品视频在线观看| 在线观看免费视频网站a站| 久久免费观看电影| 国产精品一区二区在线不卡| 午夜两性在线视频| 亚洲久久久国产精品| 国产精品美女特级片免费视频播放器 | av在线播放免费不卡| 黄片大片在线免费观看| 男男h啪啪无遮挡| 久久狼人影院| 久久九九热精品免费| 精品国产国语对白av| 国产主播在线观看一区二区| 免费看十八禁软件| 国产精品熟女久久久久浪| 女人久久www免费人成看片| 成人手机av| 老司机在亚洲福利影院| 久久中文字幕一级| 国产欧美日韩综合在线一区二区| 色婷婷久久久亚洲欧美| 精品少妇黑人巨大在线播放| 成人三级做爰电影| 日本黄色视频三级网站网址 | 18在线观看网站| 不卡av一区二区三区| 国产精品 欧美亚洲| 真人做人爱边吃奶动态| h视频一区二区三区| 亚洲人成伊人成综合网2020| 亚洲熟女毛片儿| 视频在线观看一区二区三区| 国产精品98久久久久久宅男小说| 高清黄色对白视频在线免费看| 一本综合久久免费| 久久中文字幕人妻熟女| 啦啦啦中文免费视频观看日本| 80岁老熟妇乱子伦牲交| 成人国产av品久久久| 国产不卡av网站在线观看| 一本—道久久a久久精品蜜桃钙片| 老熟妇仑乱视频hdxx| 青草久久国产| 一二三四在线观看免费中文在| 亚洲国产精品一区二区三区在线| 人人澡人人妻人| 亚洲精品一二三| 男女高潮啪啪啪动态图| 欧美精品高潮呻吟av久久| 亚洲国产av新网站| 欧美亚洲日本最大视频资源| 视频在线观看一区二区三区| 亚洲午夜精品一区,二区,三区| 肉色欧美久久久久久久蜜桃| 中文字幕高清在线视频| 日韩三级视频一区二区三区| 久久人人爽av亚洲精品天堂| 热99久久久久精品小说推荐| 欧美日韩亚洲国产一区二区在线观看 | 午夜免费鲁丝| 宅男免费午夜| 伦理电影免费视频| 80岁老熟妇乱子伦牲交| 欧美另类亚洲清纯唯美| 天天影视国产精品| 18禁裸乳无遮挡动漫免费视频| 亚洲伊人色综图| 国产精品美女特级片免费视频播放器 | 18禁裸乳无遮挡动漫免费视频| 亚洲av成人不卡在线观看播放网| 99香蕉大伊视频| 999精品在线视频| 欧美日韩成人在线一区二区| 少妇 在线观看| 中亚洲国语对白在线视频| 亚洲一区中文字幕在线| 亚洲精品中文字幕一二三四区 | 黄频高清免费视频| 一本综合久久免费| 大陆偷拍与自拍| 99久久99久久久精品蜜桃| 免费看a级黄色片| 十八禁网站免费在线| 啦啦啦中文免费视频观看日本| 99国产精品免费福利视频| 久久国产精品人妻蜜桃| 国产伦理片在线播放av一区| 91老司机精品| 深夜精品福利| 午夜免费鲁丝| 欧美黑人欧美精品刺激| 我的亚洲天堂| 18禁国产床啪视频网站| 黑人猛操日本美女一级片| 国产不卡一卡二| 免费久久久久久久精品成人欧美视频| 丁香欧美五月| 久9热在线精品视频| 国产亚洲一区二区精品| 性色av乱码一区二区三区2| 精品少妇内射三级| 亚洲成av片中文字幕在线观看| 国产精品久久电影中文字幕 | 丝袜美足系列| 宅男免费午夜| 亚洲伊人久久精品综合| 天堂俺去俺来也www色官网| 精品免费久久久久久久清纯 | 国产精品一区二区免费欧美| 亚洲第一av免费看| 国产精品美女特级片免费视频播放器 | 欧美黑人欧美精品刺激| 亚洲精品国产一区二区精华液| 亚洲伊人色综图| 亚洲熟女精品中文字幕| 国产精品1区2区在线观看. | 亚洲av成人不卡在线观看播放网| 国产成人影院久久av| 在线 av 中文字幕| 91成年电影在线观看| 人人澡人人妻人| 久热这里只有精品99| 国产97色在线日韩免费| 午夜视频精品福利| 亚洲人成电影观看| 老熟妇乱子伦视频在线观看| 国产精品av久久久久免费| 久久 成人 亚洲| 91精品三级在线观看| 最新美女视频免费是黄的| 国产主播在线观看一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 乱人伦中国视频| 免费人妻精品一区二区三区视频| 成人永久免费在线观看视频 | 多毛熟女@视频| 亚洲美女黄片视频| 欧美日韩福利视频一区二区| 变态另类成人亚洲欧美熟女 | 日韩欧美免费精品| av国产精品久久久久影院| 国产精品国产av在线观看| 少妇的丰满在线观看| 男人舔女人的私密视频| 嫁个100分男人电影在线观看| 欧美激情久久久久久爽电影 | 色94色欧美一区二区| 国产精品一区二区免费欧美| 在线观看免费视频网站a站| 99热网站在线观看| 亚洲国产成人一精品久久久| 国产黄色免费在线视频| 欧美激情久久久久久爽电影 | 热re99久久精品国产66热6| 亚洲av欧美aⅴ国产| 美国免费a级毛片| 亚洲中文av在线| 午夜免费鲁丝| 蜜桃在线观看..| 久久久国产一区二区| 午夜免费成人在线视频| 操出白浆在线播放| 亚洲精品在线美女| 午夜福利在线观看吧| 亚洲色图av天堂| 亚洲欧美一区二区三区黑人| 一本综合久久免费| 亚洲av电影在线进入| 精品少妇内射三级| 露出奶头的视频| 日日爽夜夜爽网站| 美女高潮喷水抽搐中文字幕| 高清在线国产一区| 人妻 亚洲 视频| √禁漫天堂资源中文www| 欧美黑人精品巨大| 黑人猛操日本美女一级片| 两个人免费观看高清视频| 午夜成年电影在线免费观看| 成年动漫av网址| 国产成人啪精品午夜网站| 国产真人三级小视频在线观看| 91麻豆av在线| 免费少妇av软件| 狂野欧美激情性xxxx| 久久精品熟女亚洲av麻豆精品| 欧美大码av| 捣出白浆h1v1| 日日摸夜夜添夜夜添小说| 色94色欧美一区二区| 久久精品亚洲熟妇少妇任你| 少妇裸体淫交视频免费看高清 | 啦啦啦中文免费视频观看日本| 午夜福利在线免费观看网站| 热99国产精品久久久久久7| 久久天躁狠狠躁夜夜2o2o| 久久久久久久大尺度免费视频| 亚洲精品乱久久久久久| 欧美久久黑人一区二区| 国产亚洲欧美在线一区二区| 免费日韩欧美在线观看| 欧美日韩成人在线一区二区| 老汉色∧v一级毛片| 国产精品亚洲av一区麻豆| 欧美变态另类bdsm刘玥| 精品国内亚洲2022精品成人 | 高清视频免费观看一区二区| 久久精品国产综合久久久| av超薄肉色丝袜交足视频| 少妇粗大呻吟视频| 欧美一级毛片孕妇| 美女高潮喷水抽搐中文字幕| 丰满少妇做爰视频| 久久天躁狠狠躁夜夜2o2o| 纯流量卡能插随身wifi吗| 久久久久视频综合| h视频一区二区三区| 亚洲色图综合在线观看| 欧美日韩亚洲综合一区二区三区_| 久热爱精品视频在线9| 日本精品一区二区三区蜜桃| 女人精品久久久久毛片| 男女无遮挡免费网站观看| 黄色怎么调成土黄色| 久久人妻熟女aⅴ| 久久久精品区二区三区| 免费不卡黄色视频| kizo精华| 咕卡用的链子| 亚洲人成电影观看| 女人高潮潮喷娇喘18禁视频| 亚洲国产av新网站| 欧美变态另类bdsm刘玥| 精品高清国产在线一区| 夜夜爽天天搞| 亚洲 国产 在线| 少妇的丰满在线观看| 日本欧美视频一区|