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

    具有非線性能量阱的彈性邊界約束軸向載荷梁結(jié)構(gòu)動(dòng)力學(xué)行為研究

    2023-01-03 04:37:08趙雨皓杜敬濤張樹(shù)奇陳依林
    振動(dòng)與沖擊 2022年24期
    關(guān)鍵詞:端點(diǎn)阻尼軸向

    趙雨皓, 杜敬濤, 張樹(shù)奇, 劉 楊, 陳依林

    (1. 哈爾濱工程大學(xué) 動(dòng)力與能源工程學(xué)院, 哈爾濱 150001; 2. 北京京能國(guó)際能源技術(shù)有限公司, 北京 100041)

    軸向載荷梁作為一種基本單元廣泛應(yīng)用于船舶軸系結(jié)構(gòu)振動(dòng)分析中,其中軸向載荷通常由動(dòng)力設(shè)備引入。在工程中,軸向載荷梁結(jié)構(gòu)在外部激勵(lì)下發(fā)生不利振動(dòng)。為降低軸向載荷梁結(jié)構(gòu)的過(guò)大振動(dòng),研究其動(dòng)力學(xué)行為與振動(dòng)控制方式就顯得尤為重要。

    近年來(lái),眾多學(xué)者采用傅里葉級(jí)數(shù)研究不同邊界條件下梁結(jié)構(gòu)振動(dòng)特性和動(dòng)力學(xué)行為[1-7]。在船舶軸系振動(dòng)分析中,支承軸承通常被簡(jiǎn)化為旋轉(zhuǎn)、平動(dòng)支撐剛度。傳統(tǒng)傅里葉級(jí)數(shù)由于其邊界不連續(xù)性難以應(yīng)用于具有一般邊界軸系結(jié)構(gòu)的振動(dòng)特性分析中。為提高傅里葉級(jí)數(shù)的工程適用性,Li[8]提出邊界光滑傅立葉級(jí)數(shù)用于研究具有一般邊界約束的梁結(jié)構(gòu)振動(dòng)特性。在此基礎(chǔ)上,眾多學(xué)者采用光滑邊界傅里葉級(jí)數(shù)研究了彈性邊界約束梁結(jié)構(gòu)的振動(dòng)特性[9-11]??紤]到工程中軸向載荷的存在,Xu[12-13]等利用邊界光滑傅里葉級(jí)數(shù)和能量原理建立具有一般邊界約束的軸向載荷梁結(jié)構(gòu)的振動(dòng)特性分析模型。

    非線性能量阱作為一種有效的被動(dòng)振動(dòng)控制手段由Vakakis[14]提出。非線性能量阱具有輕質(zhì)、寬頻以及能量靶向傳遞特性,自提出以來(lái)受到了國(guó)內(nèi)外學(xué)者的廣泛關(guān)注。Georgiades等[15]研究了具有非線性能量阱的Euler-Bernoulli梁動(dòng)力學(xué)行為。Samani等[16]通過(guò)最大振動(dòng)幅值和耗能估計(jì)非線性能量阱的吸振性能。Ahmadabadi等[17-18]研究了具有非線性能量阱的懸臂梁結(jié)構(gòu)非線性控制。Kani等[19-20]研究了具有非線性能量阱和經(jīng)典邊界條件梁結(jié)構(gòu)振動(dòng)控制問(wèn)題。Zhang等[21]研究了具有非線性能量阱的軸向運(yùn)動(dòng)梁結(jié)構(gòu)振動(dòng)控制,表明非線性能量阱對(duì)梁結(jié)構(gòu)振動(dòng)抑制具有良好的效果。Fang等[22]設(shè)計(jì)了一種單穩(wěn)/雙穩(wěn)非線性能量阱控制具有阻尼層懸臂梁結(jié)構(gòu)的振動(dòng)。李響[23]、王闖[24]和賀騰[25]分別采用解析與數(shù)值方法研究了非線性能量阱對(duì)梁結(jié)構(gòu)的振動(dòng)控制效果。魯正等[26]結(jié)合現(xiàn)有研究,綜述了非線性能量阱技術(shù)在振動(dòng)控制中的應(yīng)用與優(yōu)勢(shì)。姚永玉等[27]分析了非線性能量阱對(duì)懸臂梁結(jié)構(gòu)的振動(dòng)抑制特性。Zhang等[28]研究了邊界非線性能量阱對(duì)梁結(jié)構(gòu)自由、強(qiáng)迫響應(yīng)的振動(dòng)抑制效果。值得注意的是上述研究多是在特定邊界下開(kāi)展,忽略梁結(jié)構(gòu)邊界處旋轉(zhuǎn)約束,當(dāng)邊界條件改變時(shí),需對(duì)模型進(jìn)行重新建立。此外,上述研究未考慮工程中常見(jiàn)的軸向載荷且未充分研究非線性能量阱對(duì)梁結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)、減振效果的影響規(guī)律。因此,建立具有非線性能量阱的彈性邊界約束軸向載荷梁結(jié)構(gòu)振動(dòng)分析模型具有重要意義。

    本文研究具有非線性能量阱的彈性邊界約束軸向載荷梁結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)與減振特性規(guī)律。采用Galerkin截?cái)喾?Galerkin truncation method,GTM)預(yù)報(bào)具有非線性能量阱的軸向載荷梁結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)。在Galerkin離散化過(guò)程中,具有線彈性邊界的軸向載荷Euler-Bernoulli梁結(jié)構(gòu)的模態(tài)函數(shù)作為Galerkin截?cái)喾ǖ臋?quán)函數(shù)和試函數(shù)。采用諧波平衡法(harmonic balance method,HBM)驗(yàn)證Galerkin截?cái)喾ǖ挠?jì)算結(jié)果。在上述基礎(chǔ)上,研究非線性能量阱對(duì)梁結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)、減振效果的影響。

    1 理論模型

    如圖1所示,考慮具有非線性能量阱的彈性邊界約束軸向載荷梁結(jié)構(gòu)橫向振動(dòng)問(wèn)題。在梁結(jié)構(gòu)兩端引入約束彈簧對(duì)邊界條件進(jìn)行模擬;kL與kR分別為左、右邊界平動(dòng)剛度;KL與KR分別為左、右邊界旋轉(zhuǎn)剛度;CL與CR分別為左、右邊界黏性阻尼。梁結(jié)構(gòu)模型為Euler-Bernoulli梁;E,L,ρ,I與S分別為其楊氏模量、長(zhǎng)度、密度、慣性矩及截面積;P為作用在梁結(jié)構(gòu)上的軸向載荷,其中拉力載荷為負(fù),壓力載荷為正;F(x,t)為作用在梁結(jié)構(gòu)上的外部激勵(lì),在本文研究中,外部激勵(lì)為簡(jiǎn)諧激勵(lì)力,其形式為F(x,t)=δ(x-xF)·F0sin(ωt)。其中xF為外部激勵(lì)作用位置,F(xiàn)0為外部激勵(lì)幅值。非線性能量阱由黏性阻尼CNES,質(zhì)量mNES與非線性剛度knNES組成。本文考慮具有慣性效應(yīng)的非線性能量阱,非線性能量阱的有效慣性質(zhì)量(bN)為慣性系數(shù)(bc)與非線性能量阱質(zhì)量(mNES)的乘積。u(x,t)與uNES分別為梁結(jié)構(gòu)的橫向振動(dòng)位移與非線性能量阱的振動(dòng)位移。

    圖1 具有非線性能量阱的彈性邊界約束軸向載荷梁結(jié)構(gòu)振動(dòng)分析模型Fig.1 Vibration analysis model of an axially loaded beam structure with a nonlinear energy sink and elastic boundary restraints

    根據(jù)牛頓第二定律推導(dǎo)梁結(jié)構(gòu)、非線性能量阱的動(dòng)力學(xué)控制方程與振動(dòng)系統(tǒng)邊界條件。梁結(jié)構(gòu)左端點(diǎn)邊界條件為

    (1)

    (2)

    右端邊界條件為

    (3)

    (4)

    式(5)為非線性能量阱的動(dòng)力學(xué)控制方程;式(6)為梁結(jié)構(gòu)動(dòng)力學(xué)控制方程。

    knNES[uNES-u(xNES)]3=0

    (5)

    本文采用Galerkin截?cái)喾A(yù)報(bào)具有非線性能量阱的彈性邊界約束軸向載荷梁結(jié)構(gòu)的動(dòng)力學(xué)特性。在Galerkin截?cái)喾ㄖ?,選取一組滿足邊界條件的函數(shù)作為權(quán)函數(shù)與試函數(shù),實(shí)現(xiàn)梁結(jié)構(gòu)系統(tǒng)動(dòng)力學(xué)控制方程的離散??紤]到邊界阻尼的存在,通過(guò)狄利克雷函數(shù)將邊界阻尼引入梁結(jié)構(gòu)動(dòng)力學(xué)控制方程中此時(shí)邊界條件與梁結(jié)構(gòu)動(dòng)力學(xué)控制方程被改寫(xiě)為式(7)、式(8)與式(9)。由式(7)與式(8)可知,梁結(jié)構(gòu)振動(dòng)系統(tǒng)具有線彈性邊界條件,因此具有彈性邊界約束的軸向載荷梁結(jié)構(gòu)系統(tǒng)的模態(tài)函數(shù)被選取為Galerkin截?cái)喾ㄖ械臋?quán)函數(shù)與試函數(shù)。

    2 Galerkin截?cái)喾ㄅc求解

    在本文分析中,根據(jù)模態(tài)疊加原理,梁結(jié)構(gòu)振動(dòng)系統(tǒng)的橫向振動(dòng)位移展開(kāi)為如下形式

    (10)

    式中:φi(x)為具有彈性邊界約束的軸向載荷梁結(jié)構(gòu)振動(dòng)系統(tǒng)的第i階模態(tài)函數(shù);N為截?cái)鄶?shù);具有彈性邊界約束的軸向載荷梁結(jié)構(gòu)振動(dòng)系統(tǒng)的模態(tài)函數(shù)通過(guò)邊界光滑傅里葉級(jí)數(shù)與能量法求解得到。

    在上述基礎(chǔ)上,將式(10)代入式(5),得到非線性能量阱動(dòng)力學(xué)方程的具體表達(dá)形式,如式(11)所示。將式(10)代入式(9),并對(duì)其進(jìn)行Galerkin離散,得到梁結(jié)構(gòu)振動(dòng)系統(tǒng)的殘差方程,其中第m階殘差方程的形式如式(12)所示。

    (11)

    (12)

    將梁結(jié)構(gòu)振動(dòng)系統(tǒng)的第m階殘差方程簡(jiǎn)化為如下形式

    Rm1+Rm2+Rm3+Rm4+Rm5+Rm6+Rm7+Rm8+Rm9=0

    (13)

    式(13)中各項(xiàng)的具體形式為

    (14)

    (15)

    (16)

    (17)

    Rm5=ψm(xF)F0sin(ωt)

    (18)

    (19)

    (20)

    (22)

    將式(13)中的加速度項(xiàng)移至等式的另一側(cè),式(13)改寫(xiě)為

    Rm1=-(Rm2+Rm3+Rm4+Rm5+Rm6+Rm7+Rm8+Rm9)

    (23)

    根據(jù)式(14),Rm1可以展開(kāi)為

    (24)

    在上述基礎(chǔ)上,將梁結(jié)構(gòu)振動(dòng)系統(tǒng)的殘差方程整理為矩陣的形式,即

    (25)

    將式(11)與式(25)整理為標(biāo)準(zhǔn)矩陣形式,即

    (26)

    式中,RNES的具體表達(dá)式為

    (27)

    式(26)可通過(guò)4階龍格-庫(kù)塔算法直接求解,將求解結(jié)果代入式(10)即可得到具有非線性能量阱與彈性邊界約束的軸向載荷梁結(jié)構(gòu)振動(dòng)系統(tǒng)任意點(diǎn)的動(dòng)力學(xué)響應(yīng)。

    3 數(shù)值結(jié)果與分析

    按照上述理論推導(dǎo),對(duì)具有非線性能量阱的彈性邊界約束軸向載荷梁結(jié)構(gòu)振動(dòng)系統(tǒng)進(jìn)行編程仿真。首先驗(yàn)證本文模型在預(yù)報(bào)梁結(jié)構(gòu)振動(dòng)系統(tǒng)動(dòng)力學(xué)響應(yīng)時(shí)的準(zhǔn)確性和可靠性。在此基礎(chǔ)上,研究非線性能量阱參數(shù)對(duì)梁結(jié)構(gòu)振動(dòng)系統(tǒng)減振、動(dòng)力學(xué)特性的影響規(guī)律。

    3.1 模型驗(yàn)證

    本節(jié)采用Galerkin截?cái)喾A(yù)報(bào)具有非線性能量阱的彈性邊界約束軸向載荷梁結(jié)構(gòu)的動(dòng)力學(xué)響應(yīng)。采用諧波平衡法所得到的梁結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)與Galerkin截?cái)喾ㄋ媒Y(jié)果進(jìn)行對(duì)比,驗(yàn)證本文模型的正確性。在此基礎(chǔ)上,研究截?cái)鄶?shù)對(duì)Galerkin截?cái)喾ㄋ媒Y(jié)果的影響。表1列出了梁結(jié)構(gòu)材料參數(shù)、幾何參數(shù)與邊界約束彈簧剛度值。表2給出了非線性能量阱的參數(shù)值。外部激勵(lì)幅值F0=10 N,外部激勵(lì)位置xF=0。

    表1 梁結(jié)構(gòu)材料、幾何與邊界參數(shù)值Tab.1 Material, geometric and boundary parameters of beam structure

    表2 非線性能量阱參數(shù)值Tab.2 Parameters of the nonlinear energy sink

    首先驗(yàn)證基于Galerkin截?cái)喾ㄋ鶚?gòu)建模型的正確性。Galerkin截?cái)喾ǖ挠?jì)算參數(shù)取為N=M=2。為消除瞬態(tài)響應(yīng),選取0~420Te為計(jì)算域,其中Te為簡(jiǎn)諧外激勵(lì)作用周期。選取400~420Te的結(jié)果作為穩(wěn)態(tài)周期結(jié)果。4階龍格-庫(kù)塔算法的初始值取為

    (28)

    圖1為采用諧波平衡法與Galerkin截?cái)喾ǖ玫降木哂蟹蔷€性能量阱與梁結(jié)構(gòu)振動(dòng)系統(tǒng)端點(diǎn)處的動(dòng)力學(xué)響應(yīng)。由圖1可知,兩種方法所得結(jié)果吻合良好,驗(yàn)證了Galerkin截?cái)喾ㄔ谇蠼獗疚哪P蜁r(shí)的正確性。

    圖2 梁結(jié)構(gòu)穩(wěn)態(tài)幅頻響應(yīng)曲線Fig.2 Stable steady-state amplitude-frequency response curves of beam structure

    圖3給出當(dāng)截?cái)鄶?shù)分別為2,4和6時(shí)梁結(jié)構(gòu)的穩(wěn)態(tài)幅頻響應(yīng)曲線。由圖3可知,當(dāng)截?cái)鄶?shù)為4和6時(shí),梁結(jié)構(gòu)的動(dòng)力學(xué)響應(yīng)結(jié)果已經(jīng)收斂。因此,在后續(xù)研究中,選取Galerkin截?cái)喾ǖ慕財(cái)鄶?shù)為N=M=4。

    圖3 不同截?cái)鄶?shù)下梁結(jié)構(gòu)穩(wěn)態(tài)幅頻響應(yīng)曲線Fig.3 Stable steady-state amplitude-frequency response curves of beam structure with different truncated terms

    3.2 非線性能量阱與激勵(lì)位置對(duì)梁結(jié)構(gòu)振動(dòng)行為的影響

    本節(jié)研究非線性能量阱以及激勵(lì)位置對(duì)梁結(jié)構(gòu)振動(dòng)行為的影響。梁結(jié)構(gòu)邊界條件取值為kL=50 000 N/m,kR=500 N/m,KL=10 000 Nm/rad,KR=100 Nm/rad,CL=CR=0。梁結(jié)構(gòu)參數(shù)、外部激勵(lì)參數(shù)與3.1節(jié)中相同。

    首先研究非線性能量阱參數(shù)對(duì)梁結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)的影響。圖4研究非線性能量阱對(duì)梁結(jié)構(gòu)端點(diǎn)處穩(wěn)態(tài)響應(yīng)的振動(dòng)抑制效果。非線性能量阱參數(shù)取值為kNES=1×109N/m,CNES=1×10 Ns/m,mNES=0.001 kg,xNES=0.2 m。由圖4可知,非線性能量阱的存在能夠有效抑制梁結(jié)構(gòu)在第2階、第3階主共振區(qū)振動(dòng)。在第2階主共振區(qū)處非線性能量阱對(duì)左、右端點(diǎn)的吸振效果分別為41.9 %與46.3%;在第3階主共振區(qū)處非線性能量阱對(duì)左、右端點(diǎn)的吸振效果分別為66.0%與65.9%。

    圖4 knNES=1×109 N/m3時(shí)梁結(jié)構(gòu)兩端穩(wěn)態(tài)響應(yīng)的減振效果Fig.4 Vibration suppression effect for the stable responses at both ends of beam structure with knNES=1×109 N/m3

    圖5研究了當(dāng)knNES=1×1010N/m時(shí),非線性能量阱對(duì)梁結(jié)構(gòu)端點(diǎn)處穩(wěn)態(tài)響應(yīng)的減振效果。其他參數(shù)與本節(jié)前述部分一致。由圖5可知,此時(shí)非線性能量阱對(duì)梁結(jié)構(gòu)振動(dòng)系統(tǒng)端點(diǎn)處的吸振效果較knNES=1×109N/m時(shí)增強(qiáng),但在第2階主共振區(qū)處出現(xiàn)復(fù)雜動(dòng)力學(xué)行為。上述現(xiàn)象說(shuō)明非線性能量阱對(duì)梁結(jié)構(gòu)端點(diǎn)處的振動(dòng)狀態(tài)影響顯著。為深刻揭示非線性能量阱對(duì)梁結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)的影響規(guī)律,提高梁結(jié)構(gòu)的減振效果,研究其非線性剛度、阻尼以及位置對(duì)梁結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)的影響規(guī)律就顯得尤為重要。此外,根據(jù)工程經(jīng)驗(yàn),外部激勵(lì)位置亦對(duì)梁結(jié)構(gòu)的動(dòng)力學(xué)響應(yīng)具有顯著影響。因此,在本節(jié)后續(xù)研究中,研究非線性能量阱的非線性剛度、阻尼、位置以及外部激勵(lì)位置對(duì)梁結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)的影響規(guī)律。

    圖5 knNES=1×1010 N/m3時(shí)梁結(jié)構(gòu)兩端穩(wěn)態(tài)響應(yīng)的減振效果Fig.5 Vibration suppression effect for the stable responses at both ends of beam structure with knNES=1×1010 N/m3

    圖6研究非線性剛度對(duì)系統(tǒng)端點(diǎn)處動(dòng)力學(xué)響應(yīng)的影響。非線性剛度由1×108N/m3變化至1×1012N/m3。外部簡(jiǎn)諧激勵(lì)頻率為48 Hz。其他參數(shù)與本節(jié)前述部分相同。由圖6可知,幅值-非線性剛度響應(yīng)曲線存在兩個(gè)不穩(wěn)定區(qū),在不穩(wěn)定區(qū)內(nèi)系統(tǒng)呈現(xiàn)復(fù)雜的動(dòng)力學(xué)行為。在一定范圍內(nèi)增加非線性能量阱的非線性剛度對(duì)系統(tǒng)的減振具有有益的影響。

    圖6 梁結(jié)構(gòu)兩端幅值-非線性剛度響應(yīng)曲線Fig.6 Amplitude-nonlinear-stiffness response curves at both ends of beam structure

    圖7分析了非線性能量阱阻尼變化對(duì)系統(tǒng)端點(diǎn)處動(dòng)力學(xué)響應(yīng)的影響。非線性能量阱的剛度取為1.5×1010N/m3。非線性能量阱阻尼由9 Ns/m變化至25 Ns/m。外部簡(jiǎn)諧激勵(lì)頻率為48 Hz。其他參數(shù)與本節(jié)前述部分相同。由圖7可知,當(dāng)非線性能量阱阻尼較小時(shí),幅值-阻尼響應(yīng)曲線存在一個(gè)不穩(wěn)定區(qū),在不穩(wěn)定區(qū)內(nèi)梁結(jié)構(gòu)呈現(xiàn)復(fù)雜的動(dòng)力學(xué)行為。增加非線性能量阱的阻尼能夠抑制梁結(jié)構(gòu)的復(fù)雜動(dòng)力學(xué)行為且對(duì)端點(diǎn)處減振具有有益的效果。

    圖7 梁結(jié)構(gòu)兩端幅值-阻尼響應(yīng)曲線Fig.7 Amplitude-damping response curves at both ends of beam structure

    圖8研究非線性能量阱位置變化對(duì)梁結(jié)構(gòu)端點(diǎn)處動(dòng)力學(xué)響應(yīng)的影響。非線性能量阱的非線性剛度取為1×1010N/m3。非線性能量阱位置由0.1 m變化至0.9 m。外部簡(jiǎn)諧激勵(lì)頻率為48 Hz。其他參數(shù)與本節(jié)前述部分相同。由圖8可知,當(dāng)非線性能量阱靠近較硬側(cè)邊界時(shí),幅值-NES位置響應(yīng)曲線呈現(xiàn)穩(wěn)態(tài)單周期狀態(tài)。當(dāng)非線性能量阱靠近較軟側(cè)邊界時(shí),幅值-NES位置響應(yīng)曲線出現(xiàn)一個(gè)不穩(wěn)定區(qū)。在不穩(wěn)定區(qū)內(nèi)梁結(jié)構(gòu)呈現(xiàn)復(fù)雜的動(dòng)力學(xué)行為。合理的布置非線性能量阱的位置能夠避免梁結(jié)構(gòu)出現(xiàn)復(fù)雜的動(dòng)力學(xué)行為且對(duì)端點(diǎn)處的減振具有有益的效果。

    圖8 梁結(jié)構(gòu)兩端幅值-NES位置響應(yīng)曲線Fig.8 Amplitude-NES-position response curves at both ends of beam structure

    圖9研究激勵(lì)位置變化對(duì)梁結(jié)構(gòu)振動(dòng)系統(tǒng)端點(diǎn)處動(dòng)力學(xué)響應(yīng)的影響。非線性能量阱的非線性剛度取為1×1010N/m3。外部簡(jiǎn)諧激勵(lì)位置由0變化至1 m,外部簡(jiǎn)諧激勵(lì)頻率為48 Hz。其他參數(shù)與本節(jié)前述部分相同。由圖9可知,幅值-激勵(lì)位置響應(yīng)曲線存在3個(gè)不穩(wěn)定區(qū),在不穩(wěn)定區(qū)內(nèi)梁結(jié)構(gòu)呈現(xiàn)復(fù)雜的動(dòng)力學(xué)行為。合理的布置外部激勵(lì)位置能夠有效抑制梁結(jié)構(gòu)的復(fù)雜動(dòng)力學(xué)行為且對(duì)其端點(diǎn)處的減振具有有益的效果。

    綜上所述,合適的非線性能量阱參數(shù)能夠有效抑制梁結(jié)構(gòu)端點(diǎn)處第2階、第3階主共振區(qū)的振動(dòng)。在一定條件下,非線性能量阱的非線性剛度、阻尼、位置以及外部激勵(lì)位置均能夠顯著改變梁結(jié)構(gòu)的動(dòng)力學(xué)響應(yīng)特性。

    3.3 非線性能量阱對(duì)梁結(jié)構(gòu)端點(diǎn)處減振率的影響

    在工程中,梁結(jié)構(gòu)受到的外部簡(jiǎn)諧激勵(lì)通常由動(dòng)力設(shè)備引入。當(dāng)動(dòng)力設(shè)備、梁結(jié)構(gòu)安裝完成后其位置難以變更。梁結(jié)構(gòu)所受的外部簡(jiǎn)諧激勵(lì)位置亦難以變化。特別的,在船舶工程中,船舶軸系往往被簡(jiǎn)化為軸向載荷梁結(jié)構(gòu)進(jìn)行振動(dòng)特性分析,其所受的外部簡(jiǎn)諧激勵(lì)通常作用于軸系結(jié)構(gòu)端部。此外,非線性能量阱的布置通常取決于結(jié)構(gòu)的剩余空間且非線性能量阱的質(zhì)量調(diào)節(jié)范圍有限。綜合考慮工程實(shí)踐,本節(jié)研究非線性能量阱的非線性剛度、阻尼對(duì)梁結(jié)構(gòu)減振率的影響。在本節(jié)計(jì)算中,梁結(jié)構(gòu)參數(shù)、外部激勵(lì)幅值、外部激勵(lì)位置以及非線性能量阱質(zhì)量、位置參數(shù)與3.2節(jié)中相同。外部激勵(lì)頻率選取為12 Hz,46.4 Hz,103.8 Hz,分別對(duì)應(yīng)于線性系統(tǒng)的前3階固有頻率。定義梁結(jié)構(gòu)左、右端點(diǎn)處的減振率為

    (29)

    (30)

    式中,uW(·)為無(wú)非線性能量阱時(shí)具有彈性邊界約束的軸向載荷梁結(jié)構(gòu)振動(dòng)系統(tǒng)的橫向振動(dòng)位移。

    圖10研究當(dāng)外部激勵(lì)頻率為12 Hz時(shí)非線性能量阱非線性剛度、阻尼對(duì)梁結(jié)構(gòu)端點(diǎn)處減振率的影響。非線性能量阱非線性剛度由1×107N/m3變化至1×1010N/m3。非線性能量阱阻尼由10 Ns/m變化至200 Ns/m。由圖10可知,12 Hz下梁結(jié)構(gòu)端點(diǎn)處的減振率隨著非線性能量阱的非線性剛度、阻尼的變化在小范圍內(nèi)波動(dòng)。非線性剛度、阻尼的變化對(duì)梁結(jié)構(gòu)振動(dòng)系統(tǒng)較軟側(cè)邊界處的減振率微弱。對(duì)于外部激勵(lì)為12 Hz的情況,改變非線性能量阱的非線性剛度與阻尼對(duì)梁結(jié)構(gòu)端點(diǎn)處減振率的提升效果有限。

    圖10 12 Hz下NES非線性剛度、阻尼對(duì)梁結(jié)構(gòu)端點(diǎn)處減振率的影響Fig.10 Influence of the NES nonlinear stiffness and damping on the vibration suppression rate at both ends of beam structure under 12 Hz

    圖11研究外部激勵(lì)頻率為46.4 Hz時(shí)非線性能量阱非線性剛度、阻尼對(duì)梁結(jié)構(gòu)端點(diǎn)處減振率的影響。由圖11可知,46.4 Hz下梁結(jié)構(gòu)端點(diǎn)處非線性能量阱非線性剛度、阻尼的增加能夠提高梁結(jié)構(gòu)振動(dòng)系統(tǒng)端點(diǎn)處的減振率。非線性能量阱的非線性剛度對(duì)梁結(jié)構(gòu)端點(diǎn)處減振率的影響存在剛度敏感區(qū)。在剛度敏感區(qū)內(nèi)增加非線性剛度能夠顯著增加梁結(jié)構(gòu)端點(diǎn)處的減振率。對(duì)于本節(jié)所研究的梁結(jié)構(gòu)參數(shù)而言,其非線性剛度敏感區(qū)為1×109N/m3至1×1010N/m3。非線性能量阱的阻尼對(duì)梁結(jié)構(gòu)端點(diǎn)處減振率的影響存在臨界值。當(dāng)阻尼超過(guò)臨界值時(shí),梁結(jié)構(gòu)端點(diǎn)處的減振率顯著提高。

    圖11 46.4 Hz下NES非線性剛度、阻尼對(duì)梁結(jié)構(gòu)端點(diǎn)處減振率的影響Fig.11 Influence of the NES nonlinear stiffness and damping on the vibration suppression rate at both ends of beam structure under 46.4 Hz

    圖12研究外部激勵(lì)頻率為103.8 Hz時(shí)非線性能量阱的非線性剛度、阻尼對(duì)梁結(jié)構(gòu)端點(diǎn)處減振率的影響。由圖12可知,103.8 Hz下梁結(jié)構(gòu)端點(diǎn)處非線性能量阱非線性剛度、阻尼的增加能夠顯著提高梁結(jié)構(gòu)振動(dòng)系統(tǒng)端點(diǎn)處的減振率,其最佳減振效果能達(dá)到90%以上。與圖11的分析類似,非線性能量阱的非線性剛度對(duì)梁結(jié)構(gòu)端點(diǎn)處減振率的影響存在敏感區(qū)。在剛度敏感區(qū)內(nèi)增加非線性剛度能夠顯著增加梁結(jié)構(gòu)端點(diǎn)處的減振率。非線性能量阱的阻尼對(duì)梁結(jié)構(gòu)端點(diǎn)處減振率的影響存在臨界值。當(dāng)阻尼超過(guò)臨界值時(shí),梁結(jié)構(gòu)端點(diǎn)處的減振率顯著提高。值得一提的是,103.8 Hz外激勵(lì)頻率下的非線性剛度敏感區(qū)、阻尼臨界值與46.4 Hz外激勵(lì)頻率下的非線性剛度敏感區(qū)、阻尼臨界值重合。

    圖12 103.8 Hz下NES非線性剛度、阻尼對(duì)梁結(jié)構(gòu)端點(diǎn)處減振率的影響Fig.12 Influence of the NES nonlinear stiffness and damping on the vibration suppression rate at both ends of beam structure under 103.8 Hz

    綜上所述,當(dāng)激勵(lì)頻率為線性系統(tǒng)第1階固有頻率時(shí),非線性能量阱非線性剛度、阻尼對(duì)梁結(jié)構(gòu)端點(diǎn)處隔振率的影響微弱。當(dāng)激勵(lì)頻率為線性系統(tǒng)第2、第3階固有頻率時(shí),非線性能量阱非線性剛度、阻尼對(duì)梁結(jié)構(gòu)端點(diǎn)處隔振率的影響顯著。在其他條件一定時(shí),合適的非線性剛度與阻尼的組合能夠顯著提高梁結(jié)構(gòu)振動(dòng)系統(tǒng)端點(diǎn)處的減振率。在工程中,對(duì)于軸向載荷梁結(jié)構(gòu),在進(jìn)行非線性能量阱的非線性剛度、阻尼進(jìn)行優(yōu)化時(shí),首先應(yīng)確定被控對(duì)象受到的外部激勵(lì)頻率。之后,研究非線性能量阱非線性剛度的剛度敏感區(qū)與阻尼臨界值。在上述基礎(chǔ)上,選取與被控對(duì)象振動(dòng)特性匹配的非線性剛度與阻尼參數(shù),實(shí)現(xiàn)軸向載荷梁結(jié)構(gòu)的振動(dòng)控制。

    4 結(jié) 論

    本文建立了具有非線性能量阱的彈性邊界約束軸向載荷梁結(jié)構(gòu)動(dòng)力學(xué)預(yù)報(bào)分析模型。利用Galerkin截?cái)喾ㄇ蠼饬航Y(jié)構(gòu)動(dòng)力學(xué)響應(yīng),并與諧波平衡法所得結(jié)果進(jìn)行對(duì)比,驗(yàn)證了本文方法的正確性,探討了Galerkin截?cái)鄶?shù)對(duì)結(jié)果穩(wěn)定性的影響。在此基礎(chǔ)上,研究了非線性能量阱參數(shù)與激勵(lì)位置對(duì)梁結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)、減振效果的影響規(guī)律。主要結(jié)論如下:

    (1)基于Galerkin截?cái)喾ㄋ鶚?gòu)建模型能夠準(zhǔn)確預(yù)報(bào)具有非線性能量阱的彈性邊界約束軸向載荷梁結(jié)構(gòu)動(dòng)力學(xué)響應(yīng),當(dāng)截?cái)鄶?shù)為4時(shí)梁結(jié)構(gòu)動(dòng)力學(xué)響應(yīng)趨于收斂。

    (2)非線性能量阱的引入能夠有效抑制梁結(jié)構(gòu)兩端頻率響應(yīng)曲線第2、第3階主共振區(qū)內(nèi)的振動(dòng)。

    (3)非線性能量阱參數(shù)及激勵(lì)位置對(duì)梁結(jié)構(gòu)振動(dòng)狀態(tài)影響顯著。非線性能量阱非線性剛度在一定范圍內(nèi)使得梁結(jié)構(gòu)端點(diǎn)處呈現(xiàn)復(fù)雜動(dòng)力學(xué)行為;非線性能量阱阻尼增加能夠有效抑制梁結(jié)構(gòu)的復(fù)雜動(dòng)力學(xué)行為,并對(duì)梁結(jié)構(gòu)減振具有有益作用;合理的布置非線性能量阱的位置以及外部激勵(lì)位置能夠避免梁結(jié)構(gòu)出現(xiàn)復(fù)雜的動(dòng)力學(xué)行為且對(duì)梁結(jié)構(gòu)端點(diǎn)處的減振具有有益的效果。

    (4)當(dāng)激勵(lì)頻率為線性系統(tǒng)第2階、第3階固有頻率時(shí),非線性能量阱剛度、阻尼參數(shù)的增加可以改善梁結(jié)構(gòu)端點(diǎn)處的減振率;合適的非線性剛度、阻尼參數(shù)能夠顯著提高梁結(jié)構(gòu)端點(diǎn)處的減振率。

    猜你喜歡
    端點(diǎn)阻尼軸向
    非特征端點(diǎn)條件下PM函數(shù)的迭代根
    大型立式單級(jí)引黃離心泵軸向力平衡的研究
    N維不可壓無(wú)阻尼Oldroyd-B模型的最優(yōu)衰減
    關(guān)于具有阻尼項(xiàng)的擴(kuò)散方程
    具有非線性阻尼的Navier-Stokes-Voigt方程的拉回吸引子
    不等式求解過(guò)程中端點(diǎn)的確定
    荒銑加工軸向切深識(shí)別方法
    參數(shù)型Marcinkiewicz積分算子及其交換子的加權(quán)端點(diǎn)估計(jì)
    具阻尼項(xiàng)的Boussinesq型方程的長(zhǎng)時(shí)間行為
    基丁能雖匹配延拓法LMD端點(diǎn)效應(yīng)處理
    少妇高潮的动态图| 人人妻人人看人人澡| 亚洲av电影在线观看一区二区三区 | 日韩一区二区视频免费看| 久久久精品欧美日韩精品| 国产乱人视频| 精品久久久久久成人av| 91精品伊人久久大香线蕉| 亚洲美女视频黄频| 国产精品1区2区在线观看.| 直男gayav资源| 亚洲精品影视一区二区三区av| 99热网站在线观看| 大话2 男鬼变身卡| 麻豆一二三区av精品| 成人无遮挡网站| 久久久久久久久久黄片| 精品国产露脸久久av麻豆 | 特级一级黄色大片| 欧美成人午夜免费资源| 久久精品综合一区二区三区| 最近中文字幕2019免费版| 午夜福利成人在线免费观看| 欧美zozozo另类| 日韩欧美精品v在线| 内射极品少妇av片p| 男人狂女人下面高潮的视频| 久久久精品欧美日韩精品| 人妻夜夜爽99麻豆av| 亚洲人成网站高清观看| 国产精华一区二区三区| 色吧在线观看| 九九在线视频观看精品| av黄色大香蕉| 性色avwww在线观看| 久久久久久伊人网av| 小说图片视频综合网站| 久久精品国产亚洲av天美| 99久久无色码亚洲精品果冻| 亚洲欧美精品专区久久| 特级一级黄色大片| 久久精品国产99精品国产亚洲性色| 噜噜噜噜噜久久久久久91| 国产精品福利在线免费观看| 国产精品三级大全| 三级经典国产精品| 啦啦啦韩国在线观看视频| 亚洲欧美一区二区三区国产| 91在线精品国自产拍蜜月| 国产伦精品一区二区三区视频9| 天堂网av新在线| 一区二区三区免费毛片| 99久久中文字幕三级久久日本| 精品不卡国产一区二区三区| 亚洲在线自拍视频| 老司机影院成人| 欧美性感艳星| 黄片wwwwww| 看免费成人av毛片| 女的被弄到高潮叫床怎么办| 免费搜索国产男女视频| 热99在线观看视频| 大香蕉久久网| 精品无人区乱码1区二区| 亚洲最大成人av| 97人妻精品一区二区三区麻豆| 国产精品一区二区在线观看99 | 国产伦理片在线播放av一区| 晚上一个人看的免费电影| 久久韩国三级中文字幕| 九九久久精品国产亚洲av麻豆| 国产一区二区在线av高清观看| 国产午夜福利久久久久久| 欧美+日韩+精品| 丝袜喷水一区| 久热久热在线精品观看| 99在线视频只有这里精品首页| 长腿黑丝高跟| 久久精品国产自在天天线| 永久网站在线| 国产成人freesex在线| 99久久成人亚洲精品观看| 午夜精品在线福利| 久久午夜福利片| 亚洲av中文av极速乱| 久久人妻av系列| 免费观看的影片在线观看| 舔av片在线| 长腿黑丝高跟| www.色视频.com| 91av网一区二区| 国产精品综合久久久久久久免费| 国产午夜精品论理片| 国产成人精品一,二区| 国产麻豆成人av免费视频| 久久99热这里只频精品6学生 | 精品酒店卫生间| 又爽又黄a免费视频| 99久国产av精品国产电影| 男女国产视频网站| 欧美bdsm另类| 精品一区二区三区人妻视频| 日韩国内少妇激情av| 黄色配什么色好看| 嘟嘟电影网在线观看| 成人午夜高清在线视频| av免费观看日本| 日韩,欧美,国产一区二区三区 | av在线老鸭窝| 久热久热在线精品观看| 男人和女人高潮做爰伦理| 一级毛片久久久久久久久女| 国产精品国产三级国产av玫瑰| 在线天堂最新版资源| 精品国内亚洲2022精品成人| 中文乱码字字幕精品一区二区三区 | 高清日韩中文字幕在线| 日韩av不卡免费在线播放| 国产黄片美女视频| 久久久久免费精品人妻一区二区| 国产精品永久免费网站| 狂野欧美白嫩少妇大欣赏| 亚洲五月天丁香| 亚洲国产欧美在线一区| 自拍偷自拍亚洲精品老妇| 成年女人看的毛片在线观看| 成人午夜精彩视频在线观看| 亚洲人成网站在线播| 1024手机看黄色片| 免费看美女性在线毛片视频| 国产精品,欧美在线| 午夜a级毛片| 中文亚洲av片在线观看爽| 亚洲性久久影院| 亚洲av成人精品一区久久| 麻豆精品久久久久久蜜桃| 欧美xxxx黑人xx丫x性爽| 偷拍熟女少妇极品色| 国产精品一区二区三区四区免费观看| 午夜激情欧美在线| 在现免费观看毛片| 一级黄色大片毛片| 亚洲av成人精品一区久久| 国产色爽女视频免费观看| 秋霞伦理黄片| 尾随美女入室| 伦理电影大哥的女人| 国产三级中文精品| 别揉我奶头 嗯啊视频| 少妇熟女欧美另类| 男人和女人高潮做爰伦理| 18禁裸乳无遮挡免费网站照片| 久久综合国产亚洲精品| 国产精品女同一区二区软件| 亚洲精品成人久久久久久| 日韩一本色道免费dvd| 欧美一区二区精品小视频在线| 高清午夜精品一区二区三区| 国产免费男女视频| 久久精品熟女亚洲av麻豆精品 | 国产一级毛片在线| 久久午夜福利片| 日日干狠狠操夜夜爽| 卡戴珊不雅视频在线播放| 亚洲五月天丁香| 最近中文字幕高清免费大全6| 亚洲欧洲日产国产| 高清av免费在线| 国产成人免费观看mmmm| 成人毛片a级毛片在线播放| 久久久久久久久中文| 尤物成人国产欧美一区二区三区| 久久人人爽人人片av| 91精品伊人久久大香线蕉| 亚洲av电影在线观看一区二区三区 | 国产亚洲91精品色在线| 日韩精品有码人妻一区| 我的老师免费观看完整版| 天天一区二区日本电影三级| 人人妻人人看人人澡| 久久久精品大字幕| 亚洲不卡免费看| 亚洲国产精品成人久久小说| 少妇裸体淫交视频免费看高清| 国产精品一区二区三区四区久久| 黄片wwwwww| 午夜福利在线在线| 国产一区二区三区av在线| 自拍偷自拍亚洲精品老妇| 男的添女的下面高潮视频| 永久免费av网站大全| 精品不卡国产一区二区三区| 九草在线视频观看| 国产精品一区二区三区四区免费观看| 1000部很黄的大片| 边亲边吃奶的免费视频| 夜夜看夜夜爽夜夜摸| 日韩人妻高清精品专区| 亚洲中文字幕日韩| 三级男女做爰猛烈吃奶摸视频| 久久久久免费精品人妻一区二区| 中文字幕av在线有码专区| 可以在线观看毛片的网站| 特级一级黄色大片| 免费看光身美女| 女人久久www免费人成看片 | 国产精品野战在线观看| 久久精品国产亚洲网站| 久久99热6这里只有精品| 男人和女人高潮做爰伦理| 亚洲四区av| 我要搜黄色片| 亚洲欧洲日产国产| 精品久久久久久电影网 | 在线观看66精品国产| 日韩,欧美,国产一区二区三区 | h日本视频在线播放| 日韩一区二区视频免费看| 亚洲在线观看片| 国产又色又爽无遮挡免| 国产大屁股一区二区在线视频| 99久久成人亚洲精品观看| 国产亚洲av片在线观看秒播厂 | 国产亚洲5aaaaa淫片| 2022亚洲国产成人精品| 国产精品熟女久久久久浪| 有码 亚洲区| av国产久精品久网站免费入址| 三级毛片av免费| 日本一本二区三区精品| 啦啦啦观看免费观看视频高清| 日日摸夜夜添夜夜添av毛片| 亚洲欧美精品自产自拍| 久久精品国产自在天天线| 国产成人精品婷婷| 青青草视频在线视频观看| 一本久久精品| 狂野欧美白嫩少妇大欣赏| 中文字幕免费在线视频6| 欧美一级a爱片免费观看看| 亚洲成av人片在线播放无| h日本视频在线播放| 三级男女做爰猛烈吃奶摸视频| 亚洲av日韩在线播放| 久久这里只有精品中国| 最近视频中文字幕2019在线8| 亚洲婷婷狠狠爱综合网| 欧美色视频一区免费| 亚洲精品456在线播放app| 久久这里有精品视频免费| 国产老妇女一区| 亚洲丝袜综合中文字幕| 尾随美女入室| 国产爱豆传媒在线观看| 亚洲av男天堂| 欧美zozozo另类| 狂野欧美白嫩少妇大欣赏| 欧美极品一区二区三区四区| 中文字幕久久专区| 特级一级黄色大片| 1000部很黄的大片| 久久精品夜夜夜夜夜久久蜜豆| 国产成人福利小说| 久久精品久久精品一区二区三区| 丰满少妇做爰视频| 欧美xxxx性猛交bbbb| 国产乱人偷精品视频| 亚洲人成网站在线观看播放| .国产精品久久| 又粗又爽又猛毛片免费看| 中文欧美无线码| 五月玫瑰六月丁香| 青青草视频在线视频观看| 久久久久九九精品影院| 麻豆精品久久久久久蜜桃| 亚洲av成人精品一区久久| 亚洲在久久综合| 国产午夜福利久久久久久| 精品国产三级普通话版| 免费观看人在逋| 亚洲国产精品成人久久小说| 免费黄网站久久成人精品| 国产在视频线在精品| 亚洲欧美日韩东京热| 国产免费福利视频在线观看| 国产精品福利在线免费观看| 免费不卡的大黄色大毛片视频在线观看 | 国产淫片久久久久久久久| 成人亚洲精品av一区二区| 日韩在线高清观看一区二区三区| 成年女人看的毛片在线观看| 成人国产麻豆网| 亚洲av一区综合| 精品人妻熟女av久视频| 久久人妻av系列| 色综合色国产| 亚洲av中文av极速乱| 国产男人的电影天堂91| 99久久无色码亚洲精品果冻| 啦啦啦啦在线视频资源| 久久久a久久爽久久v久久| 国产白丝娇喘喷水9色精品| 免费看a级黄色片| 色哟哟·www| 日韩欧美精品免费久久| 国产精品国产三级国产专区5o | eeuss影院久久| 国产精品精品国产色婷婷| 国产精品熟女久久久久浪| 直男gayav资源| 日韩国内少妇激情av| 国产极品天堂在线| 中文精品一卡2卡3卡4更新| 91久久精品国产一区二区成人| 欧美3d第一页| 人人妻人人看人人澡| 国产免费一级a男人的天堂| 日日撸夜夜添| 69av精品久久久久久| 成年免费大片在线观看| 免费一级毛片在线播放高清视频| 精品国产露脸久久av麻豆 | 国产精品人妻久久久久久| 亚洲中文字幕一区二区三区有码在线看| 久久精品夜夜夜夜夜久久蜜豆| 成人漫画全彩无遮挡| 亚洲成人久久爱视频| 国产成人a区在线观看| 亚洲成人av在线免费| 亚洲综合色惰| 亚洲五月天丁香| 内地一区二区视频在线| 国产精品国产高清国产av| 99热网站在线观看| 亚洲国产精品成人久久小说| 国产精品永久免费网站| 亚洲在久久综合| 国产亚洲5aaaaa淫片| 一级二级三级毛片免费看| 亚洲av不卡在线观看| 免费无遮挡裸体视频| 国产精品久久久久久精品电影| 一区二区三区四区激情视频| 建设人人有责人人尽责人人享有的 | 搡女人真爽免费视频火全软件| 观看美女的网站| 一级av片app| 黄色一级大片看看| 日日撸夜夜添| 97热精品久久久久久| 国产色爽女视频免费观看| 亚洲电影在线观看av| 色吧在线观看| 性色avwww在线观看| 国产成人精品一,二区| 又粗又爽又猛毛片免费看| 狂野欧美激情性xxxx在线观看| 亚洲欧美中文字幕日韩二区| 欧美一级a爱片免费观看看| 永久免费av网站大全| 久久精品影院6| 久99久视频精品免费| 不卡视频在线观看欧美| 国产av在哪里看| 丝袜美腿在线中文| 久久精品国产亚洲av天美| 黄色配什么色好看| 日产精品乱码卡一卡2卡三| 我要看日韩黄色一级片| 最后的刺客免费高清国语| 中文字幕人妻熟人妻熟丝袜美| 一区二区三区高清视频在线| 美女黄网站色视频| 中文乱码字字幕精品一区二区三区 | 久久久精品94久久精品| a级毛色黄片| 国产精品一区二区在线观看99 | 精品一区二区三区人妻视频| www.色视频.com| 亚洲va在线va天堂va国产| 国产大屁股一区二区在线视频| 永久网站在线| 国产精品三级大全| 国产精品久久久久久精品电影小说 | 成人三级黄色视频| 深爱激情五月婷婷| 国产v大片淫在线免费观看| 小说图片视频综合网站| 国产黄a三级三级三级人| 色视频www国产| 少妇熟女aⅴ在线视频| 免费av毛片视频| 成人无遮挡网站| 成人一区二区视频在线观看| 日韩精品青青久久久久久| 偷拍熟女少妇极品色| 九九爱精品视频在线观看| 久久久久久久久久久免费av| 99久久精品热视频| 成人av在线播放网站| a级一级毛片免费在线观看| 少妇人妻一区二区三区视频| 亚洲av熟女| 波多野结衣高清无吗| 成人三级黄色视频| 久久韩国三级中文字幕| 亚洲人成网站高清观看| 婷婷色麻豆天堂久久 | 看片在线看免费视频| 中文欧美无线码| 国产在视频线精品| 亚洲丝袜综合中文字幕| 蜜臀久久99精品久久宅男| 久久韩国三级中文字幕| 国产在视频线在精品| 国产亚洲5aaaaa淫片| 成人无遮挡网站| 欧美一级a爱片免费观看看| 蜜臀久久99精品久久宅男| 亚洲高清免费不卡视频| 亚洲国产精品国产精品| 一区二区三区高清视频在线| 综合色丁香网| 观看免费一级毛片| 纵有疾风起免费观看全集完整版 | 国产成人精品久久久久久| 久久久久久大精品| 精品久久久久久成人av| 三级国产精品欧美在线观看| 青春草亚洲视频在线观看| 国产大屁股一区二区在线视频| 日韩人妻高清精品专区| 丰满少妇做爰视频| 免费av观看视频| 七月丁香在线播放| 三级毛片av免费| 内射极品少妇av片p| 免费播放大片免费观看视频在线观看 | 日韩欧美三级三区| 午夜激情福利司机影院| 麻豆成人av视频| 国产精品国产三级国产av玫瑰| 少妇熟女欧美另类| 夜夜看夜夜爽夜夜摸| 日本av手机在线免费观看| 成人欧美大片| 精品国产三级普通话版| 少妇人妻精品综合一区二区| 网址你懂的国产日韩在线| 久久精品影院6| 女的被弄到高潮叫床怎么办| 精品熟女少妇av免费看| 狂野欧美激情性xxxx在线观看| 成年av动漫网址| 成人毛片a级毛片在线播放| 日本爱情动作片www.在线观看| 国产 一区 欧美 日韩| 只有这里有精品99| 午夜老司机福利剧场| 少妇的逼水好多| 久久亚洲精品不卡| av女优亚洲男人天堂| 日本欧美国产在线视频| 亚洲精品一区蜜桃| 日韩一本色道免费dvd| 国产一区二区亚洲精品在线观看| 国产精品蜜桃在线观看| 日本黄色视频三级网站网址| 乱人视频在线观看| 久久久久久久久久久丰满| 麻豆久久精品国产亚洲av| 一个人免费在线观看电影| 国产大屁股一区二区在线视频| 日韩av在线免费看完整版不卡| 免费观看性生交大片5| 天堂网av新在线| 成人三级黄色视频| 亚洲欧美日韩高清专用| 国产毛片a区久久久久| 狂野欧美激情性xxxx在线观看| 亚洲av成人精品一区久久| 国内揄拍国产精品人妻在线| 韩国av在线不卡| 人妻夜夜爽99麻豆av| 久久午夜福利片| 特级一级黄色大片| 一级毛片电影观看 | 国产毛片a区久久久久| 我要搜黄色片| 美女被艹到高潮喷水动态| 午夜a级毛片| 国产单亲对白刺激| 国产免费一级a男人的天堂| 女人十人毛片免费观看3o分钟| 亚洲精华国产精华液的使用体验| 一边亲一边摸免费视频| videos熟女内射| 免费观看a级毛片全部| 好男人在线观看高清免费视频| 人人妻人人澡人人爽人人夜夜 | 午夜福利在线在线| 精品欧美国产一区二区三| 国产精品一区二区三区四区免费观看| 一级毛片电影观看 | 夫妻性生交免费视频一级片| 免费观看人在逋| 欧美日本视频| 菩萨蛮人人尽说江南好唐韦庄 | 日韩国内少妇激情av| 中国国产av一级| 99热精品在线国产| 亚洲av男天堂| 久久精品国产鲁丝片午夜精品| 亚洲国产精品sss在线观看| 亚洲欧美日韩高清专用| 网址你懂的国产日韩在线| 国产黄色小视频在线观看| 人人妻人人看人人澡| 国产欧美另类精品又又久久亚洲欧美| 少妇熟女欧美另类| 亚洲欧美精品自产自拍| 久久韩国三级中文字幕| 18+在线观看网站| av又黄又爽大尺度在线免费看 | 久久草成人影院| 国产精品国产三级专区第一集| 国产人妻一区二区三区在| 国产极品精品免费视频能看的| 亚洲国产精品成人综合色| 内射极品少妇av片p| 少妇裸体淫交视频免费看高清| 国产在线男女| 国产精品日韩av在线免费观看| 啦啦啦韩国在线观看视频| 亚洲怡红院男人天堂| 如何舔出高潮| 国产亚洲5aaaaa淫片| 最近2019中文字幕mv第一页| 国产精品野战在线观看| 久久久久精品久久久久真实原创| 成人亚洲精品av一区二区| 亚洲成人精品中文字幕电影| 中文字幕精品亚洲无线码一区| 欧美精品一区二区大全| 国内精品宾馆在线| 黄色日韩在线| 一区二区三区免费毛片| 91狼人影院| or卡值多少钱| 欧美潮喷喷水| 99在线人妻在线中文字幕| 黄色欧美视频在线观看| 亚洲国产欧洲综合997久久,| 美女被艹到高潮喷水动态| 亚洲国产欧美在线一区| 亚洲av成人精品一区久久| 亚洲综合色惰| 久久久久久久国产电影| av播播在线观看一区| 日韩一区二区三区影片| 国产免费一级a男人的天堂| 男女视频在线观看网站免费| 国产成人freesex在线| 舔av片在线| 亚洲丝袜综合中文字幕| 亚洲一级一片aⅴ在线观看| 在线播放国产精品三级| 91av网一区二区| 久久精品久久久久久久性| 亚洲激情五月婷婷啪啪| 一本一本综合久久| 免费av不卡在线播放| 国产精品久久久久久av不卡| 国产免费男女视频| 一级毛片久久久久久久久女| 小蜜桃在线观看免费完整版高清| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩三级伦理在线观看| 国产色爽女视频免费观看| 婷婷六月久久综合丁香| 精品99又大又爽又粗少妇毛片| 国产精品国产三级专区第一集| 水蜜桃什么品种好| 国产女主播在线喷水免费视频网站 | 在线观看av片永久免费下载| 亚洲最大成人手机在线| 国产熟女欧美一区二区| 一区二区三区乱码不卡18| 老司机影院毛片| 一个人观看的视频www高清免费观看| 岛国毛片在线播放| 岛国在线免费视频观看| 黄色一级大片看看| 91在线精品国自产拍蜜月| 国产精品乱码一区二三区的特点| 青青草视频在线视频观看| 久久亚洲国产成人精品v| 青青草视频在线视频观看| 国产伦在线观看视频一区| 观看美女的网站| 亚洲最大成人中文| 一个人观看的视频www高清免费观看| 国内精品一区二区在线观看| 久久精品久久精品一区二区三区| 人妻制服诱惑在线中文字幕| 天美传媒精品一区二区| 99久久九九国产精品国产免费| 国产免费视频播放在线视频 | 白带黄色成豆腐渣| 2021少妇久久久久久久久久久| 三级男女做爰猛烈吃奶摸视频| 18禁裸乳无遮挡免费网站照片| 熟女电影av网| 午夜精品在线福利| 午夜福利网站1000一区二区三区|