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

    非線性能量阱剛度優(yōu)化計算與振動臺試驗

    2016-11-24 06:17:43劉中坡呂西林烏建中
    振動與沖擊 2016年20期
    關(guān)鍵詞:振動臺振子線性

    劉中坡, 呂西林, 王 棟, 烏建中

    (1.同濟大學 機械與能源工程學院,上海 201804;2.同濟大學 土木工程防災(zāi)國家重點實驗室, 上海 200092)

    ?

    非線性能量阱剛度優(yōu)化計算與振動臺試驗

    劉中坡1, 呂西林2, 王 棟2, 烏建中1

    (1.同濟大學 機械與能源工程學院,上海 201804;2.同濟大學 土木工程防災(zāi)國家重點實驗室, 上海 200092)

    介紹了一種通過幾何非線性產(chǎn)生非線性回復力的NES,并對它的寬頻控制特性進行了研究。首先通過解析方式求得1∶1共振下頻率與能量之間的關(guān)系并繪制頻率能量圖,頻率能量圖直觀的表達了當NES與線性振子連接到一起的時候,其相對于線性振子的振動頻率與系統(tǒng)中能量水平有著直接關(guān)系。然后,在多自由度振動系統(tǒng)的模態(tài)空間中使用解析的方法將NES剛度與線性振子的能量、頻率與振型建立起相關(guān)的表達式。由表達式得到振動系統(tǒng)在一定能量水平下,NES達到較優(yōu)控制效果所需要的非線性剛度。根據(jù)較優(yōu)剛度計算方法,設(shè)計了NES振動控制試驗并進行振動臺試驗。為了試驗NES寬頻控制特性,振動臺臺面輸入采用了針對頻域測試的Chirp信號,試驗過程中通過改變NES的彈簧剛度與增加被控框架質(zhì)量來改變框架動力學特性的方法以檢驗NES寬頻控制的效果。試驗結(jié)果表明NES有較好的寬頻控制效果,即使其剛度偏離優(yōu)化值或被控對象動力學特性發(fā)生一定改變,NES依然能發(fā)揮較好的振動控制作用。

    NES;振動控制;振動臺試驗;非線性振動

    振動控制裝置被廣泛用于航天、汽車、機械與建筑工程中,目前大多數(shù)振動控制裝置都是基于線性吸振器。隨著對非線性振動研究的深入,非線性吸振器由于其寬頻控制的特性引起了相關(guān)研究者們的廣泛關(guān)注。ROBERSON等[1]最早提出了采用非線性振子作為吸振器進行振動控制,VAKAKIS等[2]將一個非線性振子與一個線性振子相連接,研究了系統(tǒng)在沖擊載荷下1∶1的共振響應(yīng),并提出非線性靶向能量傳遞現(xiàn)象,他將該類剛度近似為立方剛度用于振動控制的非線性振子命名為“非線性能量阱”(Nonlinear Energy Sink,NES)[3]。GOURDON等[4]通過數(shù)值方法驗證了該現(xiàn)象并指出NES優(yōu)于傳統(tǒng)的線性調(diào)頻類吸振裝置。PANAGOPOULOS等[5]研究了具有多個自由度的線性振動系統(tǒng)與NES振子相連并發(fā)生1∶1共振時的共振頻率作用范圍的問題。GENDELMAN等[6]利用非線性分叉理論對NES的動力學特性進行分析,從而揭示了單自由度主振動系統(tǒng)連接NES后的各個分叉枝,為NES的減振機理研究提供了分析基礎(chǔ);LEE等[7]采用經(jīng)驗?zāi)B(tài)分解法對動力響應(yīng)時程進行分解,繪制了頻率能量圖,這兩位研究者的成果為研究NES的振動控制提供了有效的工具。國內(nèi)學者張也弛等[8-9]對引發(fā)最優(yōu)靶向能量傳遞的條件與NES的參數(shù)優(yōu)化做了研究;孔憲仁等[10]對兩自由度NES在簡諧荷載下的力學特性與振動控制效果使用數(shù)值方法進行了研究;陳勇等[11]使用數(shù)值仿真方法研究了NES在高聳結(jié)構(gòu)振動控制中的應(yīng)用。本文通過解析方式來對試驗用框架在一定外界載荷作用下NES的剛度進行優(yōu)化計算,求得較優(yōu)的NES剛度后進行試驗驗證。在頻域中考查NES在較優(yōu)與偏離優(yōu)化值的剛度參數(shù)下振動控制效果,并通過改變試驗框架的質(zhì)量來改變框架的動力學特性來考查NES的寬頻控制效果。

    1 NES及其特性

    1.1 NES結(jié)構(gòu)與動力方程

    NES由非線性回復力裝置、質(zhì)量塊和阻尼器構(gòu)成。在本文中利用線性彈簧,通過構(gòu)造幾何非線性使其提供非線性的回復力,如圖1所示。

    圖1 非線性振子實現(xiàn)方式Fig.1 Structure of nonlinear oscillator

    質(zhì)量塊m在平衡位置時由兩根剛度相同且等長的彈簧拉住,當彈簧受力沿中垂線移動到平衡位置,可以得到回復力與位移關(guān)系的表達式:

    (1)

    式中:F為沿中垂線方向作用到質(zhì)量塊的力;k為線性彈簧的剛度;un為質(zhì)量塊相對絕對坐標系沿中線垂直方向的位移;l為彈簧被拉伸的長度。使用麥克勞林級數(shù)將式(1)展開得到下式:

    (2)

    式中:l0為彈簧未受拉力下的原始長度,在本文中質(zhì)量塊m處于水平平衡位置的時候彈簧處于不受力的狀態(tài)l0=l,代入式(2)中得到質(zhì)量塊m回復力近似表達式:

    (3)

    式中:kn=k/l2,當質(zhì)量塊位移un的絕對值在l長度40%以內(nèi)時,式(3)的誤差在5%以內(nèi)[12]。

    圖2 單自由度線性振子與NESFig.2 SDOF linear oscillator and NES

    把NES與單自由度線性振子連接到一起,系統(tǒng)的動力學方程為:

    (4)

    (5)

    1.2 NES連接到線性振動系統(tǒng)的動力學特性

    NES在與線性系統(tǒng)連接到一起時,對線性系統(tǒng)的振動控制主要是通過與線性振子發(fā)生內(nèi)共振的方式來實現(xiàn)的,內(nèi)共振形式有:主共振、亞諧波共振和非線性拍振,這些共振形式并不是同時出現(xiàn),而是隨系統(tǒng)能量水平的不同被激發(fā),但其中起主要作用的是主共振[13]。在目前的研究中大多數(shù)學者一般通過頻率能量圖來描述NES的振動控制特性,本文對這種相對于以往非線性振動分析中常用的幅頻、相頻圖較新的工具使用解析的方法進行介紹,并使用頻率能量圖來描述NES振動控制中起主要作用的1∶1主共振,由此來說明使用NES振動控制的時候,NES振動頻率與被控對象能量相關(guān)的聯(lián)系。這里的1∶1主共振指的是式(4)所定義的非線性微分方程在式(6)所定義的初始狀態(tài)下:

    x(-T/4)=u(-T/4)=0

    (6)

    滿足式(7)條件且半周期一致的周期解[13],其中T為非線性微分方程式(4)周期解的周期值。

    (7)

    使用復化平均法來求解單自由度線性振子與NES連接到一起時候的1∶1共振[15],定義復變量:

    (8)

    由式(8)可以得到使用定義的復變量來表示的NES與線性振子的位移與加速度:

    (9)

    式中:星號表示復變量的共軛關(guān)系。定義的復變量可以使用式(10)進行近似表達:

    ψ1(t)=φ1(t)ejωt,ψ2(t)=φ2(t)ejωt

    (10)

    將式(9)與(10)代入式(5)中可以得到:

    (11)

    根據(jù)復化平均法對式(11)消除快變頻率高于ω的項(e2jωt,e3jωt,…),整理可以得到:

    (12)

    引入極坐標變換:

    φ1=Aejα,φ2=Bjβ

    (13)

    式中:A與B為實幅值,α與β為實相位,將式(13)代入式(12)之后,兩個方程的實部與虛部各自為零,得到4個方程。定義相位α=β與對時間t的微分為0,將4個方程整理得到以下方程組:

    (14)

    式(14)中的第一個方程組整理并積分可以得到:

    A2+εB2=N2

    (15)

    式中:N為積分常數(shù),式(15)為保守系統(tǒng)中的全部動能。對式(14)第二、三個方程求解可以得到A,B的表達式,將其代入式(13)、(10)與式(9)中的x與u的表達式,根據(jù)歐拉公式可以得到:

    (16)

    根據(jù)式(16)可以得到保守系統(tǒng)的勢能與頻率之間的關(guān)系為:

    (17)

    由式(17)得到的關(guān)系式,即可以在FEP圖中繪制S11±脊線的曲線。

    1.3 頻率能量圖算例

    考慮由式(4)定義并去掉阻尼項的一個保守系統(tǒng)的例子,令參數(shù)ε=0.05,C=1,ω∈[0.01~2],并根據(jù)表達式(17)與(16)在Matlab中進行求解,由求解結(jié)果繪制NES的頻率與能量相關(guān)的曲線,如圖3所示。

    圖3 S11±脊線的頻率能量圖Fig.3 FEP of backbone branch S11±

    圖4 線性振子與NES的位移時程Fig.4 Displacement time history of linear oscillator and NES

    在時程數(shù)據(jù)中計算u-v,將得到的相對位移時程數(shù)據(jù)在Matlab中進行小波變換,得到分解的各個振動頻率、頻率的幅值和時間序列,將這些數(shù)據(jù)通過時間序列與振動時程中系統(tǒng)總體的能量建立關(guān)系,得到圖5:圖中x軸為系統(tǒng)的動能與勢能之和,y軸為小波變換后分解得到的振動頻率,z軸為各個分解頻率能量。對照圖5可以看出各個頻率的能量峰值對應(yīng)系統(tǒng)能量的坐標值與S11+脊線所描述的關(guān)系基本一致,可以從圖中觀察出NES這種寬頻振動的特性。由此算例說明FEP圖,相對傳統(tǒng)的幅頻響應(yīng)圖和相頻響應(yīng)圖可以更為形象的描述非線性振動中頻率變化與系統(tǒng)能量變化的關(guān)系。另外圖5中分析的結(jié)果,使用NES對被控目標特定頻率范圍的振動進行控制需要根據(jù)振動系統(tǒng)的能量水平來設(shè)計NES的參數(shù)。

    圖5 小波變換后各個頻率幅值與系統(tǒng)振動能量關(guān)系Fig.5 Variation of frequency amplitude with the variation of the system’s energy

    2 NES剛度優(yōu)化與仿真計算

    2.1 NES剛度優(yōu)化理論

    從上述分析可以得出結(jié)論在使用NES進行振動控制時,需要針對線性振子的能量水平來設(shè)計NES的參數(shù),文獻[8-10,17-19]都有涉及到NES參數(shù)優(yōu)化的內(nèi)容,在這里對NES的剛度優(yōu)化推導過程做簡要介紹。為了優(yōu)化單NES對多自由度線性振動系統(tǒng)的振動控制效果,需要建立起線性振子頻率、能量與NES剛度之間的聯(lián)系。除此之外,與其他被動振動控制裝置類似,為了獲得較好的振動控制效果,還需要結(jié)合被控系統(tǒng)的振型設(shè)計NES的安裝位置。

    為了方便進行分析,將式(4)表達的動力學模型轉(zhuǎn)到集中質(zhì)量多自由度振動系統(tǒng)的模態(tài)空間中。定義X=φq,φ為多自由度線性系統(tǒng)的模態(tài)振型矩陣,q為模態(tài)坐標向量。在本文中關(guān)注的是NES對試驗對象一階模態(tài)振動的控制,所以通過該式將方程轉(zhuǎn)換到一階模態(tài)空間之中, 如下:

    (18)

    (19)

    得到式(18)經(jīng)過變換之后的近似表達式:

    (20)

    使用式(8)與(10)定義的復變量對s,vn進行替換,并使用式(21)定義的多尺度進行展開:

    φ1=φ10+εφ11,φ2=φ20+εφ21

    (21)

    展開之后消除久期項,并使用式(22)定義的極坐標變量替換復變量:

    φ10=R1(τ1)eiδ1(τ1)

    φ20=R2(τ1)eiδ2(τ1)

    (22)

    得到方程:

    (23)

    (24)

    得到式(24),通過該式將被控的線性振子能量水平與NES的參數(shù):質(zhì)量、剛度與阻尼聯(lián)系到一起,方便了對NES振動控制性能的研究。將式(24)的第二式兩側(cè)對En求偏導,得到E1和En的解析表達式:

    (25)

    圖6 E1與En多重解關(guān)系Fig.6 Multiplicity solution between E1 and En

    為了簡化極值點的表達式,定義以下變量:

    (26)

    E1極值與對應(yīng)En的值表達式為:

    (27)

    (28)

    式(27)與式(28)建立關(guān)系,求得當滿足這個條件時的NES剛度表達式為:

    (29)

    通過以上推導過程,得到了NES對一定能量水平的線性振子進行較優(yōu)振動控制所需的剛度,為后面試驗設(shè)計提供了計算依據(jù)。

    2.2 基于試驗對象的較優(yōu)剛度參數(shù)仿真計算

    由于在試驗過程中無法直接測得框架頂層的速度,而NES剛度優(yōu)化計算又與線性振子的速度有關(guān),所以采用數(shù)值仿真計算的方法來求得被控框架的速度。仿真計算的對象為試驗中使用的小型兩層框架,每層高度0.3 m、寬度0.25 m、質(zhì)量4 kg,立柱為高0.6 m、寬0.15 m、厚1.2 mm的65Mn彈簧鋼板。試驗過程中使用振動測試系統(tǒng)中常用的Chirp信號對振動臺臺面進行基于位移的激勵,設(shè)定信號的幅值2 mm,頻率范圍1~15 Hz。從振動臺控制系統(tǒng)中獲取到如圖7所示Chirp信號激勵下振動臺臺面實際的位移時程作為計算用的外部激勵,使用有限元軟件計算得到框架弱軸向的剛度陣和質(zhì)量陣作為仿真計算用模型的參數(shù),將這些數(shù)值帶入兩自由度線性振子動力學方程并在Matlab中編寫方程積分計算的程序,計算程序中動力學微分方程的積分求解使用四階龍格庫塔方法。框架在下的空框架速度時程如圖8所示。

    圖7 Chirp信號Fig.7 Chirp signal

    圖8 仿真得到框架頂層的速度時程Fig.8 Velocity time history at roof of the frame

    3 NES振動控制的振動臺試驗

    3.1 試驗對象與設(shè)備

    試驗對象為一個小型二層鋼框架,框架詳細尺寸信息與2.2節(jié)中仿真對象的尺寸一致。經(jīng)過掃頻測試,框架一階頻率為3.8 Hz、二階頻率為10.6 Hz。振動臺使用Quanser公司出品的Shake Table Ⅱ型振動臺,該設(shè)備在全球各大高校振動相關(guān)研究課題中被廣泛使用,其質(zhì)量與振動模擬的能力都達到國際標準,是國際上機械、土木工程試驗中采用的一種標準小型振動試驗裝置。振動臺與框架的安裝狀態(tài)如圖9(a)所示。框架每層安裝一個IEPE型加速度計用來記錄框架動力反應(yīng)信息,加速度信號經(jīng)過信號調(diào)理后輸入數(shù)據(jù)采集卡。數(shù)據(jù)采集設(shè)備使用NI公司制造的9125型16位四通道同步采集卡,試驗過程中采樣頻率200 Hz。

    NES質(zhì)量塊為一通過直線軸承在光軸上滑動的鋼塊,鋼塊兩側(cè)使用線性彈簧與框架的頂板連接,結(jié)構(gòu)如圖9(b)所示。NES的質(zhì)量與2.2節(jié)中計算取值一致,為130 g。在這里需要指出的是張也馳等[9]與KERSCHEN 等[13]指出NES相對被控結(jié)構(gòu)質(zhì)量比大于等于5.6%是觸發(fā)最優(yōu)靶向能量傳遞的條件,在試驗方案設(shè)計的時候考慮到在應(yīng)用中NES的質(zhì)量較輕的話利于實施,同時參考了結(jié)構(gòu)振動中線性調(diào)頻質(zhì)量阻尼器的質(zhì)量比設(shè)定范圍,將其設(shè)定為一階模態(tài)質(zhì)量的2.5%。該質(zhì)量比雖然不能滿足動力學模型中最優(yōu)的完全能量傳遞,但在試驗中還是發(fā)揮了較好的振動控制作用。

    圖9 振動臺、試驗框架與NESFig.9 Shake table, steel frame and NES

    3.2 試驗過程

    振動臺臺面的輸入使用本文2.2節(jié)介紹過的Chirp信號,起始與截至頻率為1~15 Hz,幅值2 mm。在試驗中除了驗證仿真計算得到NES的彈簧剛度對試驗用框架加速度的控制,還改變NES彈簧的剛度試驗不同剛度對框架的振動控制效果的影響,所采用的彈簧如表1所示。另外為了考察NES對振動控制的魯棒性,在試驗過程中保持NES參數(shù)不變,改變框架的動力特性以考察其控制的寬頻特性。由于改變框架的剛度比較難,試驗中通過給框架頂板添加附加質(zhì)量塊來改變框架的振動頻率。共有4個附加質(zhì)量塊,每塊質(zhì)量為0.4 kg,試驗過程中每次遞加一個,附加質(zhì)量塊與框架的安裝如圖10所示。

    表1 試驗中使用的彈簧

    圖10 附加質(zhì)量塊Fig.10 Attached mass of the frame

    3.3 試驗結(jié)果

    NES彈簧使用表1中編號1的彈簧來近似得到仿真計算中求得的較優(yōu)剛度216 N/m,無附加質(zhì)量空框架和NES控制下框架頂層加速度的時程與加速度響應(yīng)譜對比分別如圖11(a)、(b)所示。

    圖11 框架頂層加速度時程對比Fig.11 Comparison of acceleration time history and acceleration response spectrum

    由圖11(a)所示框架頂層加速度的時程曲線可以看出:在Chirp信號作用的起始階段的5~13 s范圍內(nèi),振動臺臺面振動的頻率接近框架一階模態(tài)頻率。在該范圍內(nèi)框架頂層的加速度得到了較好的控制,最大加速度由1.9 g降低至0.79 g。在13 s之后,振動臺加載的頻率接近框架的二階頻率,框架振動以二階頻率為主。此時從圖中加速度時程曲線可以觀察到:相對無控狀態(tài),NES控制下框架頂層加速度雖然有一定量的削弱但是從加速度減小的比例上不如13 s之前。對兩個加速度時程分別計算其均方根,空框架的為0.913,NES控制下的為0.761相對空框架減少17%。

    將頂層加速度時程數(shù)據(jù)進行FFT變換,轉(zhuǎn)到頻域范圍內(nèi)進行分析,得到頂層加速度的反應(yīng)譜如圖11(b)所示。從圖11(b)中可以觀察到一階頻率3.8 Hz下加速的響應(yīng)譜幅值從0.39減小到0.11,降低了70%。二階頻率10.8 Hz下響應(yīng)譜的幅值并無變化,這個與剛度優(yōu)化計算中對一階模態(tài)振動控制的目標是一致的。另外在頂層加速度的時程曲線中,在振動臺激勵頻率接近框架二階模態(tài)頻率范圍內(nèi)的時候,NES控制下的頂層加速度相對于無控下的也有減小,這種現(xiàn)象由頻域分析的結(jié)果可以得出結(jié)論:這個階段NES還是削弱振動中一階振動的部分。

    圖12 不同剛度彈簧控制效果對比Fig.12 Comparison of control effect between different stiffness of NES spring

    表1中其他不同剛度的彈簧對試驗用框架振動控制的頻域效果對比如圖12所示,從圖中結(jié)果可以得到:① 當NES使用比優(yōu)化計算值大50%,剛度為300 N/m的彈簧時,響應(yīng)譜中一階模態(tài)頻率的峰值為0.15 g。與空框架的峰值0.39 g相比還是降低了61%,這個試驗工況通過將彈簧剛度偏移優(yōu)化值,對比空框架的動力響應(yīng),說明了NES對振動控制的寬頻特性;② 隨著剛度的增加,NES對一階模態(tài)的振動控制開始減弱,對二階模態(tài)的振動控制開始加強。這個與剛度優(yōu)化計算公式中所體現(xiàn)的隨著被控頻率的增加,需要同時增加NES剛度的關(guān)系是一致的。需要指出的是,雖然振動控制的效果在趨向二階模態(tài)振動控制的同時,一階模態(tài)的振動還是在一定程度上被削弱,這體現(xiàn)了NES的寬頻振動控制的優(yōu)點。另外,由于Chirp信號激勵的特點,在接近二階頻率時候振動臺輸入的能量相對一階頻率時候要高,使得試驗框架頂板振動的速度值相對較大,同時在NES剛度優(yōu)化計算中隨著被控結(jié)構(gòu)能量水平越高反而需要相對低能量水平時更低的彈簧剛度。這樣造成了試驗中隨著彈簧剛度的增加,剛度越來越趨近較高速度條件下對頂板二階振動控制所需的較優(yōu)的NES剛度,試驗中出現(xiàn)了對頂板二階振動控制效果開始明顯提升的現(xiàn)象。

    圖13 附加額外質(zhì)量的框架頂層加速度時程與響應(yīng)譜對比Fig.13 Comparison of acceleration time history and acceleration response spectrum of the frame with attached mass

    在通過改變試驗框架的動力特性來考察NES寬頻控制效果試驗中,將4個質(zhì)量塊同時附加在框架頂板上的時候,同時保持NES的質(zhì)量130 g與剛度參數(shù)210 N/m不變。對附加質(zhì)量的框架在無控與NES控制下分別進行試驗,外界激勵依然使用與之前試驗相同的Chirp信號,無控與有控下頂層加速度時程、響應(yīng)譜對比分別如圖13(a)、(b)所示。對圖13(a)所示的兩個加速度時程分別統(tǒng)計其均方根,無控下為0.507,NES控制下為0.411相對無控下降低了19%。從圖13(b)可以觀察到:在給框架增加了1.6 kg的附加質(zhì)量后一階頻率變?yōu)?.8 Hz。盡管一階頻率改變較大,但是從響應(yīng)譜對比可以看出NES對框架的一階振動還是起到了較好的控制效果,將框架的一階頻率對應(yīng)的加速度響應(yīng)譜最大幅值從0.27降至0.15,減少44%。由于附加質(zhì)量安裝在框架的頂部,框架的二階頻率改變并不大僅降低為10 Hz,但是由于附加質(zhì)量對框架振型的改變,NES對框架的二階振動也有了一定的控制效果,響應(yīng)譜幅值從0.23降至0.17,減少26%。除此之外,從框架頂層的加速度時程可以看出,由于附加質(zhì)量對框架的動力學特性做了較大的改變,頂板的加速度值也低于未附加的。在振動臺臺面激勵不變的前提下,可以判斷頂板的速度也低于未附加質(zhì)量的框架。在該試驗工況中,相對剛度優(yōu)化計算中設(shè)定的條件,兩方面較為顯著的不同:① 框架動力學特性不同;② 優(yōu)化對象的速度不同。在這兩個實際條件顯著偏移優(yōu)化計算設(shè)定的參數(shù)下,對比空框架與附加質(zhì)量的框架分別在無控與有控下加速度的時域與頻域的對比結(jié)果,NES依然對一階振動有著較好的控制作用,可以說明NES在振動控制中有較高的魯棒性。

    4 結(jié) 論

    本文對被動式NES進行了研究,使用頻率能量圖說明了NES與線性振子連接到一起時候的動力學特性,介紹了一種解析的方式將單個NES剛度與多自由度線性振子能量水平、頻率、振型建立起關(guān)系,并依此關(guān)系對剛度進行優(yōu)化計算的方法。

    根據(jù)優(yōu)化計算表達式與試驗對象的動力學特征,設(shè)計并開展了NES振動控制試驗,著重從頻域中考查NES的寬頻控制特性,從試驗結(jié)果得到如下結(jié)論:

    (1) 剛度優(yōu)化計算公式是有效的,根據(jù)精確的參數(shù)進行優(yōu)化計算得到的剛度來設(shè)計NES,可以有效控制優(yōu)化計算中指定的目標模態(tài)振動,但是對其他模態(tài)振動并無明顯的控制作用。由于公式中體現(xiàn)了多自由度線性振子的振型、頻率和NES安裝位置,優(yōu)化公式本身并不能體現(xiàn)NES寬頻控制的特性。

    (2) 對彈簧剛度改變試驗結(jié)果證明了NES的寬頻控制特性,在剛度超出優(yōu)化計算值50%得時候,對目標模態(tài)振動的控制效果還較為明顯。這說明,在使用NES時候因其自身剛度參數(shù)精確度對控制效果造成負面影響的敏感性并不高,可以減少制作、實施的難度與工作量。

    (3) 改變試驗框架頂板質(zhì)量的試驗結(jié)果同樣說明了NES的寬頻控制特性與魯棒性,在被控目標實際動力參數(shù)偏移計算參數(shù)范圍不大的情況下,NES依舊可以發(fā)揮振動控制作用。這使NES在類似精密設(shè)備機柜、運輸載具這類工作中質(zhì)量常常不能精準確定,而又希望有消振措施的設(shè)備中有著良好應(yīng)用場所。

    目前關(guān)于非線性能量阱在振動控制中的研究和應(yīng)用還較少,考慮到其寬頻與較為魯棒的控制特性,這一類結(jié)構(gòu)振動控制設(shè)備具有較大的研究與應(yīng)用空間。

    [1] ROBERSON R E. Synthesis of a nonlinear dynamic vibration absorber [J]. Journal of the Franklin Institute, 1952, 254(3): 205-220.

    [2] VAKAKIS A F, GENDELMAN O. Energy pumping in nonlinear mechanical oscillators: part Ⅱ: resonance capture [J]. Journal of Applied Mechanics, 2001, 68(1): 42-48.

    [3] VAKAKIS A F. Inducing passive nonlinear energy sinks in vibrating systems[J]. Transactions-American Society of Mechanical Engineers Journal of Vibration and Acoustics, 2001, 123(3): 324-332.

    [4] GOURDON E, LAMARQUE C H. Energy pumping with various nonlinear structures: numerical evidences [J].Nonlinear Dynamics, 2005, 40(3): 281-307.

    [5] PANAGOPOULOS P N, VAKAKIS A F, TSAKIRTZIS S. Transient resonant interactions of finite linear chains with essentially nonlinear end attachments leading to passive energy pumping [J]. International Journal of Solids and Structures, 2004, 41(22): 6505-6528.

    [6] GENDELMAN O, MANEVITCH L, VAKAKIS A F, et al. A degenerate bifurcation structure in the dynamics of coupled oscillators with essential stiffness nonlinearities[J]. Nonlinear Dynamics, 2003, 33(1): 1-10.

    [7] LEE Y S, KERSCHEN G, VAKAKIS A F, et al. Complicated dynamics of a linear oscillator with a light, essentially nonlinear attachment[J]. Physica D: Nonlinear Phenomena, 2005, 204(1):41-69.

    [8] 張也弛,孔憲仁,楊正賢,等.非線性吸振器的靶能量傳遞及參數(shù)設(shè)計[J].振動工程學報,2011, 34(2):111-117.

    ZHANG Yechi, KONG Xianren, YANG Zhengxian, et al. Targeted energy transfer and parameter design of a nonlinear vibration absorber[J]. Journal of Vibration Engineering, 2011, 34(2):111-117.

    [9] 張也弛,孔憲仁,張紅亮. 非線性耦合振子間的靶能量傳遞研究:保守系統(tǒng)中的完全能量傳遞[J].振動與沖擊,2012, 31(1):150-155.

    ZHANG Yechi, KONG Xianren, ZHANG Hongliang. Targeted energy transfer among coupled nonlinear oscillators:complete energy exchange in a conservative system[J]. Journal of Vibration and Shock, 2012, 31(1):150-155.

    [10] 孔憲仁,張也弛.兩自由度非線性吸振器在簡諧激勵下的振動抑制[J].航空學報, 2012, 33(6):1020-1029.

    KONG Xianren, ZHANG Yechi. Vibration suppression of a two-degree-of-freedom nonlinear energy sink under harmonic excitation[J]. Acta Aeronautica et Astronautica Sinica, 2012, 33(6):1020-1029.

    [11] 陳勇,徐羿. 基于非線性能量吸振器的高聳結(jié)構(gòu)減振分析[J].振動與沖擊,2014, 33(9):27-32.

    CHEN Yong, XU Yi. Vibration suppression analysis for a tall structure attached with a nonlinear energy sink absorber[J]. Journal of Vibration and Shock, 2014, 33(9):27-32.

    [12] CARRELLA A, BRENNAN M J, WATERS T P. Static analysis of a passive vibration isolator with quasi-zero-stiffness characteristic[J]. Journal of Sound and Vibration, 2007, 301(3/4/5):678-689.

    [13] KERSCHEN G, LEE Y S, VAKAKIS A F, et al. Irreversible passive energy transfer in coupled oscillators with essential nonlinearity[J]. SIAM Journal on Applied Mathematics, 2005, 66(2):648-679.[14] KERSCHEN G, PEETERS M, GOLINVAL J C, et al. Nonlinear normal modes, part Ⅰ: A useful framework for the structural dynamicist[J]. Mechanical Systems and Signal Processing, 2009, 23(1):170-194.

    [15] MANEVITCH L, GOURDON E, LAMARQUE C H. Towards the design of an optimal energetic sink in a strongly inhomogeneous two degree of freedom system[J]. Journal of Applied Mechanics, 2007, 74(6):1078-1086.

    [16] PILIPCHUN V N. Transient mode localization in coupled strongly nonlinear exactly solvable oscillators[J]. Nonlinear Dynamics, 2008, 51(1): 245-258.

    [17] GENDELMAN O, STAROSVETSKY Y, FELDMAN M. Attractors of harmonically forced linear oscillator with attached nonlinear energy sink I: description of response regimes[J]. Nonlinear Dynamics, 2008, 51(1/2):31-46.

    [18] VAURIGAUD B, SAVADKOOHI A T, LAMARQUE C H. Efficient targeted energy transfer with parallel nonlinear energy sinks: theory and experiment[J]. Journal of Computational and Nonlinear Dynamics, 2011, 6(4):1-10.

    [19] 熊懷,孔憲仁,劉源.阻尼對耦合非線性能量阱系統(tǒng)影響研究[J].振動與沖擊,2015, 34(11):116-121.

    XIONG Huai,KONG Xianren,LIU Yuan. Influence of structural damping on a system with nonlinear energy sinks[J]. Journal of Vibration and Shock, 2015, 34(11):116-121.

    Stiffness optimization of nonlinear energy sink and shaking table test

    LIU Zhongpo1, Lü Xilin2, WANG Dong2, WU Jianzhong1

    (1. School of Mechanical Engineering, Tongji University, Shanghai 201804, China; 2. State Key Laboratory for Disaster Reduction in Civil Engineering, Tongji University, Shanghai 200092, China)

    Nonlinear energy sink (NES) was proposed and its main feature of mitigating vibrations in broad frequency range was explored in this paper. First, the expressions of backbone branch S11± that represent 1∶1 internal resonance was derived through an analytical method. According to the expressions, the frequency-energy plot was obtained. The vibration frequency of the NES that attached to the linear oscillator and the energy levels in the linear oscillator had a direct relationship, which was clearly illustrated in the frequency-energy plot. Then, the formulas defining the relationship between the stiffness of NES and the energy of MDOF linear oscillators were derived in the modal space. So, the optimal solution of NES stiffness could be calculated with regard to the energy level of the target to be controlled. Based on the optimal computations of experimental target, shaking table experiments were designed and conducted to investigate the performance of NES vibration mitigation. To validate the wide frequency vibration controlling attributes of NES, a set of springs with different stiffness were employed and dynamic properties of the experimental frame were modified in the tests. Experimental results demonstrate that NES has broad frequency band vibration controlling feature. It works well even its stiffness deviate from the optimal value or there are some differences between the real and the designed dynamic parameters of experimental targets.

    nonlinear energy sink(NES); vibration control; shaking table experiment; nonlinear vibration

    國家自然科學基金重大研究計劃重點國際合作項目(51261120377);國家自然科學基金重大研究計劃集成項目(91315301-4)

    2015-10-28 修改稿收到日期:2016-02-28

    劉中坡 男,博士生,1980年生

    呂西林 男,教授,博士生導師,1955年生

    E-mail:lxlst@#edu.cn

    TB123;TU317.1

    A

    10.13465/j.cnki.jvs.2016.20.013

    猜你喜歡
    振動臺振子線性
    漸近線性Klein-Gordon-Maxwell系統(tǒng)正解的存在性
    基于振動臺試驗的通信機柜地震易損性分析
    彈簧振子問題的分析與求解
    線性回歸方程的求解與應(yīng)用
    二階線性微分方程的解法
    非線性Duffing擾動振子共振機制的研究
    基于近似熵和混沌振子的電力諧波檢測與估計
    大型液壓離心振動臺控制策略的仿真研究
    電磁彈簧振子實驗裝置的改進
    物理實驗(2015年10期)2015-02-28 17:36:55
    420 kV避雷器振動臺抗震試驗
    两人在一起打扑克的视频| 久久久久久亚洲精品国产蜜桃av| 国产片内射在线| 看黄色毛片网站| 亚洲三区欧美一区| 国产91精品成人一区二区三区| 亚洲av熟女| bbb黄色大片| 精品久久久久久久久久免费视频 | 成人av一区二区三区在线看| 久久国产精品男人的天堂亚洲| 嫩草影院精品99| 中文字幕人妻熟女乱码| 人人妻人人添人人爽欧美一区卜| 丝袜在线中文字幕| 国产精品九九99| 如日韩欧美国产精品一区二区三区| 色尼玛亚洲综合影院| 99国产精品一区二区三区| 亚洲成人久久性| 久久中文看片网| 啪啪无遮挡十八禁网站| 日韩欧美在线二视频| 亚洲va日本ⅴa欧美va伊人久久| 好男人电影高清在线观看| 亚洲av成人av| 国产精华一区二区三区| 国产在线精品亚洲第一网站| 亚洲国产欧美日韩在线播放| 亚洲中文av在线| 99在线视频只有这里精品首页| 丰满迷人的少妇在线观看| 国产精品一区二区精品视频观看| 亚洲视频免费观看视频| 欧美日韩视频精品一区| 久久 成人 亚洲| 亚洲精品中文字幕在线视频| 国产一区二区在线av高清观看| 嫁个100分男人电影在线观看| 午夜成年电影在线免费观看| 中文亚洲av片在线观看爽| 久热这里只有精品99| 性色av乱码一区二区三区2| 97超级碰碰碰精品色视频在线观看| 国产成人免费无遮挡视频| 亚洲精品中文字幕在线视频| 午夜福利免费观看在线| 首页视频小说图片口味搜索| 韩国av一区二区三区四区| 好男人电影高清在线观看| 五月开心婷婷网| 69精品国产乱码久久久| 成年女人毛片免费观看观看9| 韩国av一区二区三区四区| 亚洲av成人不卡在线观看播放网| 91av网站免费观看| 一级片'在线观看视频| 不卡av一区二区三区| 在线视频色国产色| 久久久久久久午夜电影 | 村上凉子中文字幕在线| 久久中文字幕人妻熟女| 欧美一级毛片孕妇| 变态另类成人亚洲欧美熟女 | 欧美成人性av电影在线观看| 大码成人一级视频| 一级毛片精品| 亚洲 国产 在线| 亚洲中文av在线| 人人妻人人澡人人看| 看片在线看免费视频| 亚洲 国产 在线| 国产黄a三级三级三级人| 久久久国产成人免费| 国产99久久九九免费精品| 高清欧美精品videossex| 又黄又粗又硬又大视频| 国产亚洲精品综合一区在线观看 | 亚洲欧美一区二区三区久久| 又黄又爽又免费观看的视频| 国产伦一二天堂av在线观看| 宅男免费午夜| 亚洲性夜色夜夜综合| 国产av精品麻豆| 啪啪无遮挡十八禁网站| 视频区图区小说| 免费在线观看黄色视频的| 日韩人妻精品一区2区三区| 亚洲国产欧美一区二区综合| 久久久久九九精品影院| 成年人黄色毛片网站| 欧美日韩福利视频一区二区| 中亚洲国语对白在线视频| 日韩欧美一区二区三区在线观看| 国产一区在线观看成人免费| 久久精品国产亚洲av香蕉五月| 嫩草影院精品99| 精品免费久久久久久久清纯| 欧美在线一区亚洲| 久久99一区二区三区| 母亲3免费完整高清在线观看| 99在线视频只有这里精品首页| 国产精品成人在线| 日韩一卡2卡3卡4卡2021年| 美女 人体艺术 gogo| 国产人伦9x9x在线观看| 久久影院123| 欧美最黄视频在线播放免费 | netflix在线观看网站| 亚洲专区国产一区二区| 国产野战对白在线观看| 啦啦啦 在线观看视频| 免费搜索国产男女视频| 黄色片一级片一级黄色片| 久久伊人香网站| 国产精品野战在线观看 | 一进一出好大好爽视频| 黑丝袜美女国产一区| 18美女黄网站色大片免费观看| 亚洲一区二区三区欧美精品| 免费不卡黄色视频| 成年人黄色毛片网站| 日日爽夜夜爽网站| 免费少妇av软件| 一级片免费观看大全| 伦理电影免费视频| 久久久久久免费高清国产稀缺| 国产精品香港三级国产av潘金莲| 丰满的人妻完整版| 男人舔女人下体高潮全视频| 99热只有精品国产| 高清黄色对白视频在线免费看| 久久久久国内视频| 在线观看免费午夜福利视频| 自拍欧美九色日韩亚洲蝌蚪91| 高清在线国产一区| 不卡av一区二区三区| 国产精品98久久久久久宅男小说| 欧美+亚洲+日韩+国产| 久久国产精品男人的天堂亚洲| 日本黄色日本黄色录像| 成年人免费黄色播放视频| 久久亚洲精品不卡| 久久香蕉国产精品| 国产成人影院久久av| 大码成人一级视频| 午夜日韩欧美国产| 国产一区二区三区视频了| 日韩欧美一区视频在线观看| 午夜亚洲福利在线播放| 日韩一卡2卡3卡4卡2021年| 人妻丰满熟妇av一区二区三区| 老司机在亚洲福利影院| 天天躁狠狠躁夜夜躁狠狠躁| 成人亚洲精品一区在线观看| 欧美乱妇无乱码| 亚洲情色 制服丝袜| 黑人巨大精品欧美一区二区mp4| av在线播放免费不卡| 性色av乱码一区二区三区2| 99riav亚洲国产免费| 亚洲人成电影观看| 丝袜美腿诱惑在线| 天堂中文最新版在线下载| tocl精华| 成人18禁在线播放| 亚洲自偷自拍图片 自拍| 久热爱精品视频在线9| 亚洲少妇的诱惑av| 十八禁网站免费在线| 久久精品91无色码中文字幕| 操出白浆在线播放| 久久久久精品国产欧美久久久| 午夜免费鲁丝| 久久精品亚洲av国产电影网| 久久精品国产99精品国产亚洲性色 | 一本综合久久免费| 亚洲精品国产区一区二| 国产欧美日韩一区二区三| 成人亚洲精品一区在线观看| 男人操女人黄网站| 久久人妻av系列| 亚洲欧美日韩另类电影网站| 精品国产一区二区三区四区第35| 中文字幕另类日韩欧美亚洲嫩草| 国产极品粉嫩免费观看在线| 精品国产乱子伦一区二区三区| 一边摸一边做爽爽视频免费| 一区二区日韩欧美中文字幕| 亚洲一码二码三码区别大吗| 亚洲一区高清亚洲精品| 中文字幕高清在线视频| 亚洲专区中文字幕在线| 成人免费观看视频高清| 久久这里只有精品19| 精品卡一卡二卡四卡免费| 天天躁夜夜躁狠狠躁躁| www.www免费av| 真人做人爱边吃奶动态| 99热只有精品国产| 亚洲avbb在线观看| 免费在线观看影片大全网站| 亚洲黑人精品在线| 搡老熟女国产l中国老女人| 他把我摸到了高潮在线观看| 成人国语在线视频| 国产成人欧美| 中亚洲国语对白在线视频| 美女福利国产在线| 国产一卡二卡三卡精品| 亚洲欧美精品综合久久99| 国产成人欧美| 国产成人精品久久二区二区91| 黑人欧美特级aaaaaa片| 国产又色又爽无遮挡免费看| aaaaa片日本免费| 老司机在亚洲福利影院| 免费观看人在逋| 成人影院久久| 国产精品久久久av美女十八| 久99久视频精品免费| 日韩免费高清中文字幕av| 亚洲伊人色综图| 免费搜索国产男女视频| 国产成人精品久久二区二区免费| 91精品三级在线观看| 成人三级做爰电影| 国产精品影院久久| 国产av又大| 国产免费男女视频| 亚洲av电影在线进入| 国产高清videossex| 国产三级在线视频| 久9热在线精品视频| 美女国产高潮福利片在线看| 中文字幕另类日韩欧美亚洲嫩草| 国产成人精品久久二区二区免费| 亚洲激情在线av| 国产单亲对白刺激| www.自偷自拍.com| 久久久国产成人免费| 精品福利观看| 最近最新中文字幕大全免费视频| 亚洲欧美一区二区三区久久| 久久香蕉国产精品| 精品国产乱子伦一区二区三区| av在线播放免费不卡| 国产精品秋霞免费鲁丝片| 少妇 在线观看| 天天添夜夜摸| 中文字幕av电影在线播放| av视频免费观看在线观看| 69av精品久久久久久| 国产激情久久老熟女| 好男人电影高清在线观看| 十八禁人妻一区二区| 亚洲一区二区三区不卡视频| 一进一出抽搐gif免费好疼 | 成人国产一区最新在线观看| 色在线成人网| 99在线人妻在线中文字幕| 欧美日韩视频精品一区| 国产成人免费无遮挡视频| 国产男靠女视频免费网站| 日韩欧美一区二区三区在线观看| 午夜福利免费观看在线| 国产成人一区二区三区免费视频网站| 男男h啪啪无遮挡| 精品人妻1区二区| 一夜夜www| 自拍欧美九色日韩亚洲蝌蚪91| 日韩精品中文字幕看吧| 国产精品久久久人人做人人爽| 天堂动漫精品| 一级,二级,三级黄色视频| 欧洲精品卡2卡3卡4卡5卡区| 亚洲专区国产一区二区| 91麻豆精品激情在线观看国产 | 少妇 在线观看| 国产欧美日韩精品亚洲av| 97超级碰碰碰精品色视频在线观看| 热99国产精品久久久久久7| 国产高清激情床上av| 桃红色精品国产亚洲av| 黄色毛片三级朝国网站| 国产男靠女视频免费网站| 亚洲aⅴ乱码一区二区在线播放 | 在线观看66精品国产| 成年人免费黄色播放视频| 精品久久蜜臀av无| 九色亚洲精品在线播放| 91精品三级在线观看| 色尼玛亚洲综合影院| 国产极品粉嫩免费观看在线| 纯流量卡能插随身wifi吗| 丝袜人妻中文字幕| 在线观看66精品国产| 欧美成人性av电影在线观看| 热99国产精品久久久久久7| 亚洲狠狠婷婷综合久久图片| 亚洲男人的天堂狠狠| 国产伦一二天堂av在线观看| 女人精品久久久久毛片| 亚洲精品av麻豆狂野| 日韩免费av在线播放| 女生性感内裤真人,穿戴方法视频| 国产精品美女特级片免费视频播放器 | 午夜两性在线视频| 成年人免费黄色播放视频| 国产精品亚洲av一区麻豆| 国产精品一区二区精品视频观看| 十八禁网站免费在线| 桃色一区二区三区在线观看| 成年版毛片免费区| 日韩三级视频一区二区三区| 最新在线观看一区二区三区| 久久国产精品影院| 日本黄色日本黄色录像| 精品久久蜜臀av无| 热re99久久精品国产66热6| 99国产综合亚洲精品| 制服人妻中文乱码| 女人高潮潮喷娇喘18禁视频| 久久精品成人免费网站| 久久精品亚洲精品国产色婷小说| 宅男免费午夜| 亚洲精品国产一区二区精华液| 在线视频色国产色| 精品欧美一区二区三区在线| 午夜福利欧美成人| 亚洲avbb在线观看| 午夜福利欧美成人| 午夜老司机福利片| 丝袜美足系列| 久久伊人香网站| 亚洲色图av天堂| 99精品欧美一区二区三区四区| 一级片免费观看大全| 国产亚洲精品综合一区在线观看 | 午夜福利欧美成人| 免费在线观看黄色视频的| 欧美人与性动交α欧美软件| 在线永久观看黄色视频| 国产精品美女特级片免费视频播放器 | 国产三级黄色录像| 亚洲欧美精品综合久久99| 国产精品1区2区在线观看.| 国产一卡二卡三卡精品| 免费看十八禁软件| 免费一级毛片在线播放高清视频 | 国产精品乱码一区二三区的特点 | 亚洲自拍偷在线| 久久中文看片网| 这个男人来自地球电影免费观看| 久久亚洲真实| 国产三级黄色录像| 女人精品久久久久毛片| 九色亚洲精品在线播放| 五月开心婷婷网| 热99re8久久精品国产| 自拍欧美九色日韩亚洲蝌蚪91| 成熟少妇高潮喷水视频| 国产激情欧美一区二区| 国产精品影院久久| 免费女性裸体啪啪无遮挡网站| 精品午夜福利视频在线观看一区| 久久香蕉精品热| 国产黄a三级三级三级人| 99久久国产精品久久久| 日韩大尺度精品在线看网址 | 老司机深夜福利视频在线观看| 丝袜在线中文字幕| 性色av乱码一区二区三区2| 99国产精品99久久久久| 久热这里只有精品99| 十八禁人妻一区二区| 亚洲av片天天在线观看| 伊人久久大香线蕉亚洲五| 男女下面插进去视频免费观看| 伦理电影免费视频| 亚洲欧美精品综合一区二区三区| 国产一区在线观看成人免费| a级毛片在线看网站| 精品久久久久久久久久免费视频 | 欧美日本中文国产一区发布| 亚洲欧美激情综合另类| 少妇 在线观看| 天天添夜夜摸| 欧美大码av| 欧美日本亚洲视频在线播放| 午夜免费激情av| 国产熟女xx| 1024香蕉在线观看| 美女 人体艺术 gogo| 亚洲成人免费电影在线观看| 精品国产乱码久久久久久男人| 我的亚洲天堂| 美女大奶头视频| 免费在线观看黄色视频的| 97人妻天天添夜夜摸| 亚洲性夜色夜夜综合| 人妻丰满熟妇av一区二区三区| 少妇被粗大的猛进出69影院| 精品久久蜜臀av无| 麻豆久久精品国产亚洲av | 免费看a级黄色片| 操美女的视频在线观看| 久久狼人影院| 国产高清国产精品国产三级| 亚洲成av片中文字幕在线观看| 精品国内亚洲2022精品成人| 黑丝袜美女国产一区| 亚洲av成人一区二区三| 18禁美女被吸乳视频| 在线十欧美十亚洲十日本专区| 精品福利永久在线观看| 国产亚洲av高清不卡| 香蕉国产在线看| 悠悠久久av| 日本一区二区免费在线视频| 精品久久久久久,| 欧美日韩亚洲国产一区二区在线观看| 两个人免费观看高清视频| 91大片在线观看| 长腿黑丝高跟| 搡老熟女国产l中国老女人| 国产有黄有色有爽视频| 一边摸一边抽搐一进一出视频| 一级,二级,三级黄色视频| 久久精品91蜜桃| 久99久视频精品免费| 亚洲少妇的诱惑av| 熟女少妇亚洲综合色aaa.| 日韩欧美在线二视频| 久久天堂一区二区三区四区| 脱女人内裤的视频| 叶爱在线成人免费视频播放| 水蜜桃什么品种好| 两性夫妻黄色片| 成人av一区二区三区在线看| 国产亚洲欧美98| 一a级毛片在线观看| 搡老乐熟女国产| 99国产精品免费福利视频| 精品高清国产在线一区| 日日干狠狠操夜夜爽| av免费在线观看网站| 大香蕉久久成人网| 午夜精品国产一区二区电影| 亚洲一区二区三区不卡视频| 久久热在线av| 国产男靠女视频免费网站| 99在线视频只有这里精品首页| 亚洲五月色婷婷综合| 51午夜福利影视在线观看| 久久 成人 亚洲| 日本精品一区二区三区蜜桃| 黑人猛操日本美女一级片| 亚洲成人免费av在线播放| av福利片在线| 黄色 视频免费看| 大香蕉久久成人网| 欧美 亚洲 国产 日韩一| 国产午夜精品久久久久久| 欧美久久黑人一区二区| 久久精品国产亚洲av香蕉五月| av网站在线播放免费| 免费看a级黄色片| 久久这里只有精品19| 97碰自拍视频| 午夜老司机福利片| 12—13女人毛片做爰片一| 欧美日韩视频精品一区| 午夜两性在线视频| 丝袜美腿诱惑在线| 69精品国产乱码久久久| 亚洲精品中文字幕在线视频| 真人一进一出gif抽搐免费| 久久香蕉精品热| 日韩高清综合在线| 亚洲人成网站在线播放欧美日韩| 久久精品亚洲av国产电影网| 国产一区二区三区视频了| 亚洲一卡2卡3卡4卡5卡精品中文| 一区二区三区激情视频| 久久久久国内视频| 亚洲熟女毛片儿| 日韩av在线大香蕉| 岛国在线观看网站| 99国产综合亚洲精品| 国产亚洲精品一区二区www| 国产av精品麻豆| 亚洲自偷自拍图片 自拍| 18禁美女被吸乳视频| 精品福利观看| 色精品久久人妻99蜜桃| 久热爱精品视频在线9| 久久精品成人免费网站| 在线观看免费视频日本深夜| 一边摸一边抽搐一进一小说| 亚洲欧美日韩另类电影网站| 久久99一区二区三区| 国产精品国产高清国产av| 黄频高清免费视频| 亚洲一区二区三区不卡视频| 亚洲七黄色美女视频| 亚洲av五月六月丁香网| 免费女性裸体啪啪无遮挡网站| 50天的宝宝边吃奶边哭怎么回事| 一区福利在线观看| 色综合欧美亚洲国产小说| 亚洲国产看品久久| 亚洲国产中文字幕在线视频| av欧美777| 国产激情久久老熟女| av视频免费观看在线观看| 久久国产精品人妻蜜桃| 国产欧美日韩精品亚洲av| 婷婷丁香在线五月| 九色亚洲精品在线播放| 日韩大码丰满熟妇| 动漫黄色视频在线观看| 午夜精品久久久久久毛片777| 乱人伦中国视频| 可以免费在线观看a视频的电影网站| 日韩 欧美 亚洲 中文字幕| 久久久国产一区二区| 他把我摸到了高潮在线观看| 久久久久久久久免费视频了| 久久久水蜜桃国产精品网| 国产aⅴ精品一区二区三区波| 免费在线观看完整版高清| 亚洲欧美日韩高清在线视频| 亚洲久久久国产精品| 日日夜夜操网爽| 国产精品自产拍在线观看55亚洲| 欧美乱色亚洲激情| 日韩av在线大香蕉| 国产亚洲精品久久久久久毛片| 国产av又大| 波多野结衣高清无吗| 中出人妻视频一区二区| 久久精品aⅴ一区二区三区四区| 亚洲情色 制服丝袜| 日韩三级视频一区二区三区| 亚洲av成人一区二区三| 满18在线观看网站| 亚洲五月色婷婷综合| 欧美日韩视频精品一区| 午夜免费成人在线视频| 校园春色视频在线观看| 黄片小视频在线播放| 脱女人内裤的视频| 成在线人永久免费视频| 色老头精品视频在线观看| 十八禁人妻一区二区| 精品电影一区二区在线| 首页视频小说图片口味搜索| 80岁老熟妇乱子伦牲交| 人人澡人人妻人| 淫秽高清视频在线观看| 日韩国内少妇激情av| 欧美成人性av电影在线观看| 侵犯人妻中文字幕一二三四区| 青草久久国产| 精品电影一区二区在线| 国产97色在线日韩免费| 免费少妇av软件| 国产精品永久免费网站| 中文字幕人妻丝袜一区二区| 久久精品国产综合久久久| 精品福利永久在线观看| 动漫黄色视频在线观看| 99久久精品国产亚洲精品| 日韩精品免费视频一区二区三区| 最新美女视频免费是黄的| 久久香蕉国产精品| 老司机深夜福利视频在线观看| 一个人观看的视频www高清免费观看 | 电影成人av| 一边摸一边做爽爽视频免费| 午夜精品久久久久久毛片777| 久久人人97超碰香蕉20202| 欧美中文日本在线观看视频| 少妇被粗大的猛进出69影院| 人妻丰满熟妇av一区二区三区| 成在线人永久免费视频| 欧美激情久久久久久爽电影 | 桃色一区二区三区在线观看| 少妇被粗大的猛进出69影院| 高潮久久久久久久久久久不卡| 视频区欧美日本亚洲| 国产精品乱码一区二三区的特点 | 两性夫妻黄色片| 神马国产精品三级电影在线观看 | 纯流量卡能插随身wifi吗| 日本vs欧美在线观看视频| 美女午夜性视频免费| 丰满人妻熟妇乱又伦精品不卡| 成年人免费黄色播放视频| 淫妇啪啪啪对白视频| 成人免费观看视频高清| 免费日韩欧美在线观看| 中国美女看黄片| 久久久久久免费高清国产稀缺| av片东京热男人的天堂| 国产精品久久久人人做人人爽| 看黄色毛片网站| 久久久久久久精品吃奶| 欧美乱色亚洲激情| 免费搜索国产男女视频| 丝袜在线中文字幕| 可以免费在线观看a视频的电影网站| 国产一卡二卡三卡精品| 亚洲av片天天在线观看|