• <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避雷器振動臺抗震試驗
    亚洲国产欧美人成| 国产私拍福利视频在线观看| 欧美精品国产亚洲| 午夜精品在线福利| 精品日产1卡2卡| 中文资源天堂在线| 亚洲精品自拍成人| 91久久精品国产一区二区三区| 成人无遮挡网站| 久久久国产成人精品二区| 一级毛片aaaaaa免费看小| 成年女人看的毛片在线观看| 国产亚洲91精品色在线| 国产麻豆成人av免费视频| 99久国产av精品国产电影| 搡女人真爽免费视频火全软件| 国产一区二区三区av在线 | 我要搜黄色片| 日日摸夜夜添夜夜爱| 日韩强制内射视频| 人体艺术视频欧美日本| 亚洲精华国产精华液的使用体验 | 日本免费一区二区三区高清不卡| 亚洲人成网站在线观看播放| 国产高清三级在线| 亚洲av中文字字幕乱码综合| 尾随美女入室| 国内精品一区二区在线观看| 看十八女毛片水多多多| 欧美日韩一区二区视频在线观看视频在线 | 亚洲国产精品合色在线| 亚洲激情五月婷婷啪啪| 久久精品国产自在天天线| 男女做爰动态图高潮gif福利片| 99热只有精品国产| 国产一区二区亚洲精品在线观看| 久久精品国产清高在天天线| 成人国产麻豆网| 搡女人真爽免费视频火全软件| 国产成人91sexporn| 99在线人妻在线中文字幕| 麻豆国产97在线/欧美| 校园人妻丝袜中文字幕| 中文字幕制服av| 人妻制服诱惑在线中文字幕| 搡女人真爽免费视频火全软件| 91久久精品国产一区二区成人| 免费看日本二区| 中文字幕av在线有码专区| 成年版毛片免费区| 寂寞人妻少妇视频99o| 91午夜精品亚洲一区二区三区| 此物有八面人人有两片| 亚洲乱码一区二区免费版| 日韩欧美 国产精品| 国产精品福利在线免费观看| 精品国产三级普通话版| 国产一区亚洲一区在线观看| 一区二区三区四区激情视频 | 一个人看的www免费观看视频| 看非洲黑人一级黄片| 亚洲欧洲国产日韩| 青青草视频在线视频观看| 18禁在线无遮挡免费观看视频| 精品无人区乱码1区二区| 小蜜桃在线观看免费完整版高清| 国产三级在线视频| 91久久精品国产一区二区成人| 欧美另类亚洲清纯唯美| 亚洲中文字幕日韩| 变态另类丝袜制服| 国产视频内射| 伦精品一区二区三区| 91aial.com中文字幕在线观看| 黑人高潮一二区| 免费不卡的大黄色大毛片视频在线观看 | 九九久久精品国产亚洲av麻豆| 国产精品无大码| 亚洲欧洲国产日韩| 精品午夜福利在线看| 青春草亚洲视频在线观看| 久久6这里有精品| 免费观看的影片在线观看| 老女人水多毛片| 欧美日韩一区二区视频在线观看视频在线 | 嫩草影院入口| 国产 一区 欧美 日韩| 全区人妻精品视频| 国产一区二区激情短视频| 成人性生交大片免费视频hd| 国产伦一二天堂av在线观看| 男女做爰动态图高潮gif福利片| 一区二区三区免费毛片| 成人综合一区亚洲| 日本五十路高清| 在线观看美女被高潮喷水网站| 桃色一区二区三区在线观看| 欧美不卡视频在线免费观看| 久久久久久久久久黄片| 一级黄片播放器| 麻豆成人av视频| www日本黄色视频网| 少妇猛男粗大的猛烈进出视频 | 亚洲欧洲国产日韩| 久久久久久伊人网av| 欧美最新免费一区二区三区| 亚洲成人久久性| 一区福利在线观看| 国产探花极品一区二区| 欧美成人精品欧美一级黄| 亚洲不卡免费看| 99热这里只有是精品50| 91精品一卡2卡3卡4卡| 一个人看的www免费观看视频| 边亲边吃奶的免费视频| 九九久久精品国产亚洲av麻豆| 久久精品综合一区二区三区| 色视频www国产| www.av在线官网国产| 国产高清激情床上av| 日韩欧美精品免费久久| 99热精品在线国产| 国产三级中文精品| 毛片一级片免费看久久久久| 精品国内亚洲2022精品成人| 国产精品福利在线免费观看| 美女cb高潮喷水在线观看| 成人一区二区视频在线观看| 国产av在哪里看| 国产黄片视频在线免费观看| 国产色爽女视频免费观看| 悠悠久久av| 国产一区二区三区在线臀色熟女| 别揉我奶头 嗯啊视频| 熟女电影av网| 我要看日韩黄色一级片| 日韩欧美国产在线观看| 日产精品乱码卡一卡2卡三| 一夜夜www| 亚洲美女搞黄在线观看| 国产精品综合久久久久久久免费| av视频在线观看入口| 尤物成人国产欧美一区二区三区| 亚洲精品亚洲一区二区| 午夜福利在线观看吧| 亚洲人成网站高清观看| 91aial.com中文字幕在线观看| 丝袜美腿在线中文| 99热这里只有是精品50| 18禁黄网站禁片免费观看直播| 亚洲高清免费不卡视频| 久久久久久久久中文| 亚洲精品日韩av片在线观看| 亚洲欧美日韩高清专用| 国产免费男女视频| 国产男人的电影天堂91| 嘟嘟电影网在线观看| 三级男女做爰猛烈吃奶摸视频| 国产一区亚洲一区在线观看| 极品教师在线视频| 12—13女人毛片做爰片一| 久久综合国产亚洲精品| av在线老鸭窝| 日产精品乱码卡一卡2卡三| 99久久无色码亚洲精品果冻| 国产在线男女| 久久久久久伊人网av| 国产成人aa在线观看| 国产精品99久久久久久久久| 亚洲熟妇中文字幕五十中出| 中文字幕免费在线视频6| 亚洲国产高清在线一区二区三| 淫秽高清视频在线观看| 午夜久久久久精精品| 22中文网久久字幕| 一个人观看的视频www高清免费观看| 国产成人精品一,二区 | 可以在线观看的亚洲视频| 欧美人与善性xxx| 久久久久久伊人网av| 狂野欧美白嫩少妇大欣赏| 国产毛片a区久久久久| 国产精品一区二区三区四区久久| 亚洲av不卡在线观看| 亚洲电影在线观看av| 久久久精品94久久精品| 天堂影院成人在线观看| 国产亚洲5aaaaa淫片| 国产精品久久久久久av不卡| 久久久久免费精品人妻一区二区| 99久久精品热视频| 韩国av在线不卡| 三级经典国产精品| 嫩草影院精品99| 国产av一区在线观看免费| 国产爱豆传媒在线观看| 能在线免费看毛片的网站| 久久久欧美国产精品| 亚洲欧洲国产日韩| 亚洲一区高清亚洲精品| 国产高潮美女av| 日本av手机在线免费观看| 中文在线观看免费www的网站| 日韩精品青青久久久久久| 99久久无色码亚洲精品果冻| 日韩在线高清观看一区二区三区| 国产三级中文精品| 我的女老师完整版在线观看| 日本撒尿小便嘘嘘汇集6| 日韩大尺度精品在线看网址| 九草在线视频观看| 黄片无遮挡物在线观看| 精品久久国产蜜桃| 卡戴珊不雅视频在线播放| 色视频www国产| 免费av观看视频| 国产极品精品免费视频能看的| 久久婷婷人人爽人人干人人爱| 69av精品久久久久久| 午夜精品国产一区二区电影 | 国产精品三级大全| 国产在线精品亚洲第一网站| 男女下面进入的视频免费午夜| 人妻少妇偷人精品九色| 亚洲国产高清在线一区二区三| 偷拍熟女少妇极品色| 国产亚洲av片在线观看秒播厂 | 小蜜桃在线观看免费完整版高清| 国产成年人精品一区二区| 欧美3d第一页| 国产伦一二天堂av在线观看| 久久久欧美国产精品| 精品国内亚洲2022精品成人| 久久欧美精品欧美久久欧美| 综合色av麻豆| 少妇熟女欧美另类| a级一级毛片免费在线观看| 亚洲激情五月婷婷啪啪| 综合色av麻豆| 可以在线观看的亚洲视频| 久久久久久久亚洲中文字幕| 国产白丝娇喘喷水9色精品| 亚洲精品成人久久久久久| 我要搜黄色片| 色视频www国产| 天天一区二区日本电影三级| 久久人人爽人人爽人人片va| 中文在线观看免费www的网站| 黄片无遮挡物在线观看| 中文字幕精品亚洲无线码一区| .国产精品久久| 美女xxoo啪啪120秒动态图| 亚洲天堂国产精品一区在线| 少妇丰满av| 国产黄色视频一区二区在线观看 | av黄色大香蕉| 可以在线观看的亚洲视频| 高清毛片免费看| 精品国内亚洲2022精品成人| 变态另类成人亚洲欧美熟女| 国产不卡一卡二| 激情 狠狠 欧美| 啦啦啦啦在线视频资源| 伊人久久精品亚洲午夜| av又黄又爽大尺度在线免费看 | 久久精品国产亚洲网站| 一进一出抽搐gif免费好疼| 男女视频在线观看网站免费| 日本熟妇午夜| 久久综合国产亚洲精品| 亚洲欧美成人精品一区二区| 男女边吃奶边做爰视频| 精品国内亚洲2022精品成人| 亚洲av男天堂| 欧洲精品卡2卡3卡4卡5卡区| 91av网一区二区| 一进一出抽搐动态| 国产高潮美女av| 中国美女看黄片| 伊人久久精品亚洲午夜| 我要搜黄色片| 国产伦在线观看视频一区| 哪里可以看免费的av片| 国产在线精品亚洲第一网站| 国产精品免费一区二区三区在线| 亚洲国产精品久久男人天堂| 少妇人妻一区二区三区视频| 亚洲欧美中文字幕日韩二区| 欧美三级亚洲精品| 国产成人影院久久av| 熟妇人妻久久中文字幕3abv| 成人漫画全彩无遮挡| 一级毛片aaaaaa免费看小| av免费在线看不卡| 天堂影院成人在线观看| 少妇丰满av| 性欧美人与动物交配| 天天躁夜夜躁狠狠久久av| 国产亚洲av嫩草精品影院| 亚洲经典国产精华液单| 免费搜索国产男女视频| 免费无遮挡裸体视频| 欧美性感艳星| 国产成人一区二区在线| av在线天堂中文字幕| 久久这里有精品视频免费| 亚洲av免费在线观看| 亚洲久久久久久中文字幕| 亚洲成人精品中文字幕电影| 久久99热6这里只有精品| 日韩人妻高清精品专区| 熟女人妻精品中文字幕| 国产女主播在线喷水免费视频网站 | 在线观看一区二区三区| 女的被弄到高潮叫床怎么办| 3wmmmm亚洲av在线观看| 女的被弄到高潮叫床怎么办| 少妇被粗大猛烈的视频| 极品教师在线视频| eeuss影院久久| 国产视频内射| 黄片无遮挡物在线观看| 成人二区视频| 五月伊人婷婷丁香| 精品日产1卡2卡| 亚洲国产精品sss在线观看| 久久精品国产亚洲网站| 只有这里有精品99| 欧美色欧美亚洲另类二区| 久久久久久九九精品二区国产| 给我免费播放毛片高清在线观看| 国产一级毛片七仙女欲春2| 亚洲精品影视一区二区三区av| 人妻制服诱惑在线中文字幕| 69人妻影院| 色5月婷婷丁香| 中文字幕熟女人妻在线| 久久亚洲国产成人精品v| a级毛片免费高清观看在线播放| 成人美女网站在线观看视频| 国产一区亚洲一区在线观看| 中文字幕久久专区| 天堂√8在线中文| 国产在线精品亚洲第一网站| 亚洲熟妇中文字幕五十中出| 一级毛片aaaaaa免费看小| 欧美精品国产亚洲| av专区在线播放| 国产亚洲91精品色在线| 久久这里有精品视频免费| 一边亲一边摸免费视频| 99热精品在线国产| 国产精品日韩av在线免费观看| 欧美一级a爱片免费观看看| 久久久久久久久久黄片| 大又大粗又爽又黄少妇毛片口| 国产高清视频在线观看网站| 欧美潮喷喷水| a级毛片免费高清观看在线播放| 国产欧美日韩精品一区二区| a级毛片免费高清观看在线播放| a级一级毛片免费在线观看| 床上黄色一级片| 啦啦啦韩国在线观看视频| АⅤ资源中文在线天堂| 欧美成人a在线观看| 色5月婷婷丁香| 欧美性感艳星| 久久鲁丝午夜福利片| 免费观看人在逋| 久久亚洲国产成人精品v| 青春草亚洲视频在线观看| 中国美白少妇内射xxxbb| 一级二级三级毛片免费看| 亚洲国产日韩欧美精品在线观看| 深爱激情五月婷婷| 两性午夜刺激爽爽歪歪视频在线观看| 日本爱情动作片www.在线观看| 欧美日韩国产亚洲二区| 狠狠狠狠99中文字幕| 天天一区二区日本电影三级| 国产精品久久久久久av不卡| 国产精品一区二区三区四区免费观看| 日日摸夜夜添夜夜爱| 欧美激情久久久久久爽电影| 好男人视频免费观看在线| 夫妻性生交免费视频一级片| 毛片女人毛片| 久久精品91蜜桃| 在现免费观看毛片| 午夜福利高清视频| 久久精品国产亚洲网站| 日本色播在线视频| 亚洲国产精品久久男人天堂| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 97人妻精品一区二区三区麻豆| 成人鲁丝片一二三区免费| 亚洲精品日韩av片在线观看| 午夜福利成人在线免费观看| www.色视频.com| 久久久久久久久久久丰满| 久久久成人免费电影| 尤物成人国产欧美一区二区三区| 亚洲18禁久久av| 国产探花在线观看一区二区| 久久久精品欧美日韩精品| 伦精品一区二区三区| 狠狠狠狠99中文字幕| 三级经典国产精品| 免费观看的影片在线观看| 国产精品一区二区性色av| 久99久视频精品免费| 国产三级中文精品| 国产单亲对白刺激| 日韩制服骚丝袜av| 日韩视频在线欧美| 最近视频中文字幕2019在线8| 日韩成人伦理影院| av在线播放精品| 三级经典国产精品| 亚洲av免费高清在线观看| 一进一出抽搐gif免费好疼| 三级毛片av免费| 村上凉子中文字幕在线| 高清毛片免费看| 美女黄网站色视频| 午夜视频国产福利| 嫩草影院新地址| 婷婷精品国产亚洲av| 午夜免费激情av| 我的老师免费观看完整版| 精品人妻熟女av久视频| 日韩欧美精品v在线| 精品99又大又爽又粗少妇毛片| 色综合亚洲欧美另类图片| 国产综合懂色| 色综合亚洲欧美另类图片| 亚洲av不卡在线观看| 亚洲精品久久久久久婷婷小说 | 亚洲综合色惰| 日本撒尿小便嘘嘘汇集6| 色哟哟哟哟哟哟| 一区二区三区免费毛片| 久久久久国产网址| 成人漫画全彩无遮挡| 最后的刺客免费高清国语| 免费观看精品视频网站| 91麻豆精品激情在线观看国产| 亚洲精品乱码久久久v下载方式| 成人国产麻豆网| 性欧美人与动物交配| 黄片无遮挡物在线观看| 欧美3d第一页| 国产午夜精品久久久久久一区二区三区| 好男人在线观看高清免费视频| 在线a可以看的网站| 狠狠狠狠99中文字幕| 国产伦一二天堂av在线观看| 久久久久久大精品| 亚洲欧美日韩高清在线视频| 久久99热这里只有精品18| 中文字幕制服av| 精品人妻视频免费看| 99热6这里只有精品| 日本免费一区二区三区高清不卡| 99久国产av精品国产电影| 插逼视频在线观看| 亚洲精品影视一区二区三区av| 九草在线视频观看| 日韩成人av中文字幕在线观看| 欧美高清成人免费视频www| 有码 亚洲区| 青春草视频在线免费观看| 国产精品一及| 日本在线视频免费播放| 亚洲欧美成人综合另类久久久 | 男女做爰动态图高潮gif福利片| 欧美丝袜亚洲另类| 欧美日韩一区二区视频在线观看视频在线 | 国产亚洲5aaaaa淫片| 亚洲国产精品成人久久小说 | 国模一区二区三区四区视频| 国产黄片美女视频| 免费人成视频x8x8入口观看| 两性午夜刺激爽爽歪歪视频在线观看| 有码 亚洲区| av视频在线观看入口| 51国产日韩欧美| 亚洲中文字幕一区二区三区有码在线看| 一区二区三区免费毛片| 黄色配什么色好看| 噜噜噜噜噜久久久久久91| 国产精品人妻久久久久久| 欧美性猛交╳xxx乱大交人| kizo精华| 欧美一区二区亚洲| 久久这里有精品视频免费| 一级毛片电影观看 | 高清日韩中文字幕在线| 1024手机看黄色片| 白带黄色成豆腐渣| 最近中文字幕高清免费大全6| 亚洲,欧美,日韩| 久久综合国产亚洲精品| 亚洲成av人片在线播放无| 亚洲精品色激情综合| 变态另类丝袜制服| 男女那种视频在线观看| 久久精品国产亚洲av香蕉五月| 蜜桃久久精品国产亚洲av| 18+在线观看网站| 草草在线视频免费看| 少妇人妻精品综合一区二区 | 你懂的网址亚洲精品在线观看 | 天天躁夜夜躁狠狠久久av| 亚洲第一区二区三区不卡| 日日干狠狠操夜夜爽| 欧美潮喷喷水| 精品人妻视频免费看| 亚洲精品国产av成人精品| 在线天堂最新版资源| 精品国产三级普通话版| 国产v大片淫在线免费观看| 亚洲成人精品中文字幕电影| 精品欧美国产一区二区三| 波野结衣二区三区在线| 亚洲精品影视一区二区三区av| 亚洲在久久综合| 我的女老师完整版在线观看| 性欧美人与动物交配| 欧美日韩综合久久久久久| 成年女人永久免费观看视频| 国产精品久久电影中文字幕| av在线观看视频网站免费| 禁无遮挡网站| 亚洲最大成人中文| 91久久精品国产一区二区三区| 免费人成在线观看视频色| 一边亲一边摸免费视频| 十八禁国产超污无遮挡网站| 国产真实伦视频高清在线观看| 日本欧美国产在线视频| h日本视频在线播放| 乱人视频在线观看| 人人妻人人看人人澡| 亚洲av免费在线观看| 国内久久婷婷六月综合欲色啪| 成人三级黄色视频| 人妻系列 视频| 成人毛片60女人毛片免费| 久久午夜亚洲精品久久| 国产精品.久久久| 淫秽高清视频在线观看| 国产精品野战在线观看| 丰满乱子伦码专区| 久久这里有精品视频免费| 精品久久国产蜜桃| 国产在线精品亚洲第一网站| 深夜精品福利| 精品久久久久久久久久免费视频| 亚洲自偷自拍三级| 中文字幕久久专区| 亚洲美女视频黄频| 99热这里只有精品一区| 自拍偷自拍亚洲精品老妇| 麻豆成人午夜福利视频| 精品不卡国产一区二区三区| 春色校园在线视频观看| 精品熟女少妇av免费看| 男人的好看免费观看在线视频| 亚洲中文字幕日韩| 亚洲国产色片| 岛国毛片在线播放| 免费看美女性在线毛片视频| 在线国产一区二区在线| 久久久久久国产a免费观看| 亚洲三级黄色毛片| 亚洲自偷自拍三级| 国产在线精品亚洲第一网站| 久久久久久久久大av| 亚洲美女搞黄在线观看| www日本黄色视频网| 精品久久久久久成人av| 欧美bdsm另类| 免费黄网站久久成人精品| 五月伊人婷婷丁香| 精品久久久久久成人av| 搡女人真爽免费视频火全软件| 国产 一区 欧美 日韩| 免费不卡的大黄色大毛片视频在线观看 | 国产一级毛片在线| 黄片无遮挡物在线观看| 久久久久性生活片| 亚洲真实伦在线观看| 嘟嘟电影网在线观看| 亚洲乱码一区二区免费版| 国产视频内射| 国产三级在线视频| 欧美成人精品欧美一级黄| 欧美日韩在线观看h| 久久精品国产99精品国产亚洲性色| 亚洲乱码一区二区免费版| 国产精品蜜桃在线观看 | 欧美性猛交黑人性爽| 18+在线观看网站| av在线天堂中文字幕| 久久久久久久久久成人| 岛国在线免费视频观看| 看十八女毛片水多多多| 亚洲五月天丁香| 亚洲一区二区三区色噜噜| 亚洲av.av天堂| 男人狂女人下面高潮的视频| 精华霜和精华液先用哪个| 精品日产1卡2卡|