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

    一種機(jī)械系統(tǒng)非線性類別辨識(shí)方法

    2015-03-07 02:22:35權(quán)雙璐劉杰韓羅峰何安琦李兵
    關(guān)鍵詞:測(cè)系統(tǒng)機(jī)械系統(tǒng)類別

    權(quán)雙璐,劉杰,韓羅峰,何安琦,李兵

    (西安交通大學(xué)機(jī)械工程學(xué)院, 710049, 西安)

    ?

    一種機(jī)械系統(tǒng)非線性類別辨識(shí)方法

    權(quán)雙璐,劉杰,韓羅峰,何安琦,李兵

    (西安交通大學(xué)機(jī)械工程學(xué)院, 710049, 西安)

    針對(duì)含剛度非線性連續(xù)體系統(tǒng)非線性類別辨識(shí)問(wèn)題,提出了一種基于瞬時(shí)頻率峰值檢測(cè)的機(jī)械系統(tǒng)非線性類別辨識(shí)方法。該方法首先根據(jù)被測(cè)系統(tǒng)結(jié)構(gòu)特點(diǎn),分別建立含間隙、立方剛度等5種機(jī)械系統(tǒng)常見(jiàn)的剛度非線性系統(tǒng)簡(jiǎn)化模型;通過(guò)研究各模型自由衰減振動(dòng)瞬時(shí)頻率分布特點(diǎn),建立不同非線性模型瞬時(shí)頻率曲線庫(kù);結(jié)合圖形匹配算法,匹配被測(cè)系統(tǒng)實(shí)測(cè)瞬時(shí)頻率與曲線庫(kù)曲線,以匹配結(jié)果作為辨識(shí)指標(biāo),實(shí)現(xiàn)非線性類別辨識(shí)。將基于瞬時(shí)頻率檢測(cè)的非線性類別辨識(shí)方法應(yīng)用于含間隙非線性懸臂梁系統(tǒng),仿真分析與實(shí)驗(yàn)計(jì)算表明:被測(cè)系統(tǒng)瞬時(shí)頻率曲線與間隙非線性曲線庫(kù)匹配較好,匹配值明顯小于與其他非線性曲線庫(kù)匹配結(jié)果,驗(yàn)證了該非線性類別辨識(shí)方法的有效性。

    機(jī)械系統(tǒng);非線性類別;瞬時(shí)頻率;圖形匹配算法

    在對(duì)非線性機(jī)械系統(tǒng)進(jìn)行動(dòng)力學(xué)分析時(shí),系統(tǒng)特性描述是最重要也是最困難的工作,主要包括非線性作用位置定位和非線性類別描述兩大問(wèn)題[1]。精確的系統(tǒng)特性描述對(duì)系統(tǒng)模型的建立及非線性參數(shù)的辨識(shí)至關(guān)重要,一旦系統(tǒng)特性描述不準(zhǔn)確將直接導(dǎo)致參數(shù)辨識(shí)失敗[2]。然而,在機(jī)械系統(tǒng)非線性類別辨識(shí)方面,由于導(dǎo)致非線性振動(dòng)產(chǎn)生的因素眾多,如幾何因素,慣性作用、阻尼、邊界條件等,使得非線性類別辨識(shí)充滿挑戰(zhàn)[3]。文獻(xiàn)[4]通過(guò)波特圖及奈奎斯特圖區(qū)分了阻尼非線性與剛度非線性;文獻(xiàn)[5]利用伏爾加特級(jí)數(shù)實(shí)現(xiàn)了多項(xiàng)式非線性類別辨識(shí);文獻(xiàn)[6]用恢復(fù)力曲線圖(RFS)直觀地區(qū)分了機(jī)械系統(tǒng)中常見(jiàn)的幾種單自由度非線性;文獻(xiàn)[2]利用頻譜分析辨識(shí)了含立方剛度與平方剛度非線性懸臂板系統(tǒng);文獻(xiàn)[7]建立了非線性頻域恢復(fù)力曲線庫(kù),實(shí)現(xiàn)了單自由度和多自由度系統(tǒng)非線性類別辨識(shí)。雖然上述研究均取得了令人矚目的成果,但對(duì)于機(jī)械系統(tǒng)故障分析中常見(jiàn)的含剛度非線性連續(xù)體系統(tǒng),由于不同剛度非線性振動(dòng)響應(yīng)間無(wú)明顯差別,加之連續(xù)體系統(tǒng)各單元間運(yùn)動(dòng)相互影響,使得其剛度非線性類別辨識(shí)更加困難,很少有人在這方面進(jìn)行研究[4],目前還沒(méi)有一種方法能有效地解決上述問(wèn)題[8]。

    本文針對(duì)含剛度非線性連續(xù)體系統(tǒng),利用被測(cè)系統(tǒng)瞬時(shí)頻率分布特點(diǎn),提出了一種基于瞬時(shí)頻率峰值檢測(cè)的系統(tǒng)非線性類別辨識(shí)方法。該方法首先根據(jù)被測(cè)系統(tǒng)結(jié)構(gòu)特點(diǎn)建立含不同剛度非線性的系統(tǒng)簡(jiǎn)化模型,仿真求解各模型自由衰減振動(dòng)響應(yīng)及其瞬時(shí)頻率分布曲線;通過(guò)研究不同非線性模型瞬時(shí)頻率分布特點(diǎn),建立各模型瞬時(shí)頻率曲線庫(kù);結(jié)合圖形匹配算法,匹配被測(cè)系統(tǒng)瞬時(shí)頻率與曲線庫(kù)曲線,以匹配結(jié)果作為辨識(shí)指標(biāo),最終實(shí)現(xiàn)系統(tǒng)剛度非線性類別辨識(shí)。與傳統(tǒng)的非線性類別辨識(shí)方法相比,本文方法可有效應(yīng)用于連續(xù)體系統(tǒng)非線性類別辨識(shí),并經(jīng)仿真與實(shí)驗(yàn)驗(yàn)證了方法的有效性。

    1 剛度非線性系統(tǒng)

    在機(jī)械系統(tǒng)故障診斷中,由于剛度非線性引起的機(jī)械故障頻發(fā),常常受到人們關(guān)注[9]。本文以剛度非線性系統(tǒng)為主要研究對(duì)象,表1列出了幾種常見(jiàn)的剛度非線性系統(tǒng)及其恢復(fù)力曲線。

    根據(jù)牛頓第二定律,上述剛度非線性系統(tǒng)自由衰減振動(dòng)動(dòng)力學(xué)方程為

    (1)

    式中:u(t)為系統(tǒng)響應(yīng);M、C分別為系統(tǒng)質(zhì)量矩陣、阻尼矩陣;F(u)為系統(tǒng)非線性恢復(fù)力。

    當(dāng)F(u)=Ku(t)時(shí),系統(tǒng)為線性系統(tǒng)(K為系統(tǒng)剛度),系統(tǒng)固有頻率與振動(dòng)幅值無(wú)關(guān)

    (2)

    當(dāng)F(u)≠Ku(t)時(shí),系統(tǒng)呈非線性振動(dòng)。系統(tǒng)振動(dòng)頻率不再固定不變,而是由系統(tǒng)基頻和一系列高階諧波疊加而成,隨振動(dòng)幅值變化不斷變化。

    表1 剛度非線性系統(tǒng)及其恢復(fù)力曲線

    2 剛度非線性類別辨識(shí)

    2.1 瞬時(shí)頻率檢測(cè)法

    根據(jù)式(1)及正則化攝動(dòng)法[10],表1所列的剛度非線性系統(tǒng)自由衰減振動(dòng)的瞬時(shí)頻率可估計(jì)為系統(tǒng)振動(dòng)幅值及其非線性剛度參數(shù)的函數(shù),即

    (3)

    式中:a(t)為系統(tǒng)振動(dòng)幅值;k1,k2,…,kn為系統(tǒng)非線性剛度。

    根據(jù)式(3)可知,由于剛度非線性函數(shù)不同,不同剛度非線性系統(tǒng)自由衰減振動(dòng)的瞬時(shí)頻率將呈不同的分布規(guī)律,因此可根據(jù)瞬時(shí)頻率分布實(shí)現(xiàn)剛度非線性類別辨識(shí)。

    對(duì)非線性時(shí)變系統(tǒng),常用時(shí)頻分布來(lái)估計(jì)系統(tǒng)瞬時(shí)頻率[11]。本文采用時(shí)頻聚集性能較好且抗干擾強(qiáng)的SPWVD(smoothed pseudo-Winger distri-butions)分布[12],其分布規(guī)律為

    (4)

    式中:h(t)和g(t)為平滑函數(shù)且滿足h(0)=g(0)=1;u為系統(tǒng)時(shí)域響應(yīng);ξ為系統(tǒng)頻域響應(yīng);f(u)為系統(tǒng)響應(yīng)函數(shù);f*(u)為f(u)的復(fù)共軛函數(shù)。

    其時(shí)頻

    寬度與u(t)相同。記fa(t)=a(t)ejφ(t),則系統(tǒng)瞬時(shí)頻率為

    (5)

    式中:a(t)為幅值響應(yīng);φ(t)為系統(tǒng)角速度。

    本文運(yùn)用峰值檢測(cè)法進(jìn)行系統(tǒng)振動(dòng)瞬時(shí)頻率估計(jì),即以每個(gè)時(shí)間點(diǎn)SPWVD最大值對(duì)應(yīng)的頻率作為信號(hào)瞬時(shí)頻率的估計(jì)結(jié)果。以圖1單自由度非線性系統(tǒng)為例,圖2給出了表1所列剛度非線性系統(tǒng)瞬時(shí)頻率分布曲線,系統(tǒng)仿真參數(shù)見(jiàn)表2。

    圖1 單自由度剛度非線性系統(tǒng)

    系統(tǒng)k/N·m-1a0/mk1/N·m-1k2/N·m-1初速度/m·s-1A00.16320.3B00.16320.3C00.16321580.3D00.11586320.3E15830000.3

    (a)系統(tǒng)A (b)系統(tǒng)B (c)系統(tǒng)C

    (d)系統(tǒng)D (e)系統(tǒng)E(硬彈簧) (f)系統(tǒng)E(軟彈簧)圖2 含剛度弱非線性單自由度系統(tǒng)瞬時(shí)頻率

    從圖2中可得,立方剛度系統(tǒng)E瞬時(shí)頻率分布曲線與攝動(dòng)法估計(jì)的瞬時(shí)頻率變化趨勢(shì)相同,即

    (6)

    為二次曲線,可直接從其他幾類非線性系統(tǒng)中區(qū)分出來(lái)。系統(tǒng)A、D,系統(tǒng)B、C瞬時(shí)頻率分布曲線相似,這是由于系統(tǒng)A為系統(tǒng)D的一種特殊情況即k1=0,系統(tǒng)B為系統(tǒng)C的一種特殊情況即k2=0。

    2.2 非線性類別辨識(shí)

    為有效地辨識(shí)A、B、C、D 4類分段剛度非線性系統(tǒng),本文結(jié)合被測(cè)系統(tǒng)瞬時(shí)頻率分布特征和圖形匹配算法實(shí)現(xiàn)非線性類別辨識(shí),其技術(shù)路線如圖3所示。

    圖3 系統(tǒng)非線性分類方法技術(shù)路線圖

    該類別辨識(shí)方法基本分為以下3個(gè)步驟。

    (1)根據(jù)被測(cè)系統(tǒng)結(jié)構(gòu)特點(diǎn)及剛度非線性系統(tǒng)恢復(fù)力曲線,建立含不同剛度非線性的被測(cè)系統(tǒng)簡(jiǎn)化模型,求解相同初始位移下各模型自由衰減振動(dòng)。

    (2)采用SPWVD峰值檢測(cè)法獲得不同模型的瞬時(shí)頻率變化曲線,建立各非線性模型瞬時(shí)頻率曲線庫(kù)。

    (3)結(jié)合圖形匹配算法,匹配被測(cè)系統(tǒng)實(shí)測(cè)自由衰減振動(dòng)瞬時(shí)頻率曲線與曲線庫(kù)曲線,匹配程度最好的系統(tǒng)即為該被測(cè)系統(tǒng)主要非線性類別。

    以圖1所示的單自由度非線性系統(tǒng)為被測(cè)系統(tǒng),根據(jù)系統(tǒng)結(jié)構(gòu)特點(diǎn),建立含表1所列A~D4種剛度非線性系統(tǒng)簡(jiǎn)化模型,并求解系統(tǒng)自由衰減振動(dòng)響應(yīng),仿真參數(shù)見(jiàn)表3,非線性項(xiàng)參數(shù)與表2相同。

    表3 簡(jiǎn)化模型參數(shù)

    (a)系統(tǒng)A (b)系統(tǒng)B

    (c)系統(tǒng)C (d)系統(tǒng)D圖4 單自由度剛度非線性系統(tǒng)瞬時(shí)頻率數(shù)據(jù)庫(kù)

    根據(jù)峰值檢測(cè)法估計(jì)各模型自由衰減振動(dòng)及其瞬時(shí)頻率分布曲線。為使建立的曲線庫(kù)具有一定的可比性,瞬時(shí)頻率分布曲線需進(jìn)行歸一化處理,圖4給出了歸一化瞬時(shí)頻率曲線庫(kù),縱坐標(biāo)以系統(tǒng)基頻為基準(zhǔn),橫坐標(biāo)以仿真總時(shí)長(zhǎng)為基準(zhǔn),即

    (7)

    式中:ω(t)為系統(tǒng)瞬時(shí)頻率;ω0為系統(tǒng)基頻;t為仿真時(shí)間;T為仿真總時(shí)長(zhǎng)。

    (a)與曲線庫(kù)A對(duì)比 (b)與曲線庫(kù)B對(duì)比

    圖5給出了被測(cè)系統(tǒng)為間隙非線性系統(tǒng)A時(shí),系統(tǒng)實(shí)測(cè)瞬時(shí)頻率分布與曲線庫(kù)曲線對(duì)比圖。從圖中可明顯看出,被測(cè)系統(tǒng)瞬時(shí)頻率曲線與曲線庫(kù)中A類間隙非線性系統(tǒng)瞬時(shí)頻率曲線接近??紤]結(jié)合圖形匹配算法,通過(guò)匹配被測(cè)系統(tǒng)實(shí)測(cè)瞬時(shí)頻率與曲線庫(kù)曲線,實(shí)現(xiàn)非線性類別辨識(shí)。本文以兩曲線間距離作為匹配參數(shù),即

    (c)與曲線庫(kù)C對(duì)比 (d)與曲線庫(kù)D對(duì)比圖5 被測(cè)系統(tǒng)與曲線庫(kù)瞬時(shí)頻率對(duì)比結(jié)果

    (8)

    式中:YE為被測(cè)系統(tǒng)瞬時(shí)頻率分布曲線縱坐標(biāo);YB為曲線數(shù)據(jù)庫(kù)內(nèi)各模型瞬時(shí)頻率分布曲線縱坐標(biāo);N為采樣點(diǎn)數(shù)。

    系統(tǒng)瞬時(shí)頻率分布曲線與曲線庫(kù)瞬時(shí)頻率曲線的匹配結(jié)果見(jiàn)表4。從表中可以看出,當(dāng)被測(cè)系統(tǒng)與曲線庫(kù)對(duì)應(yīng)類別一致時(shí),對(duì)應(yīng)的DS值最小,即兩曲線匹配程度最好,從而實(shí)現(xiàn)了非線性類別的定量辨識(shí)。

    表4 單自由度非線性系統(tǒng)類別辨識(shí)DS值

    隨機(jī)選取50組不同間隙與間隙剛度的非線性單自由度系統(tǒng)作為被測(cè)系統(tǒng)。圖6給出了這50組樣本自由衰減振動(dòng)瞬時(shí)頻率分布曲線與曲線庫(kù)匹配結(jié)果。從圖6中可以明顯看出,被測(cè)系統(tǒng)瞬時(shí)頻率分布與曲線庫(kù)內(nèi)間隙非線性匹配程度最好,從而確定被測(cè)系統(tǒng)為間隙非線性系統(tǒng)。

    圖6 含間隙單自由度非線性系統(tǒng)類別的辨識(shí)結(jié)果

    3 連續(xù)體系統(tǒng)非線性類別辨識(shí)

    3.1 傳統(tǒng)方法非線性類別辨識(shí)

    由于連續(xù)體系統(tǒng)各單元間振動(dòng)的相互耦合,使得傳統(tǒng)的非線性類別辨識(shí)方法很難實(shí)現(xiàn)連續(xù)體系統(tǒng)非線性類別的辨識(shí)。以圖7所示的剛度非線性懸臂梁為例,懸臂梁任意點(diǎn)的位移響應(yīng)用絕對(duì)位移u(x,t)表示,該懸臂梁系統(tǒng)的參數(shù)如下:長(zhǎng)度l、橫截面積A、截面慣性矩I、密度ρ、彈性模量E、非線性位置距離固定端的距離lc、表1所列的剛度非線性項(xiàng)f(u)。

    圖7 非線性懸臂梁系統(tǒng)

    圖8給出了立方剛度與間隙非線性系統(tǒng)波特圖。從圖中可以看出:同一種非線性系統(tǒng)在一定輸入輸出下,系統(tǒng)幅頻響應(yīng)往往呈現(xiàn)不同的分布。此外,分段非線性系統(tǒng)波特圖受非線性影響嚴(yán)重,因此頻響函數(shù)的形狀不能作為確定系統(tǒng)剛度非線性類別的依據(jù)[2]。

    (a)立方剛度非線性系統(tǒng) (b)含間隙非線性系統(tǒng)圖8 剛度非線性系統(tǒng)波特圖

    圖9對(duì)比了含間隙單自由度系統(tǒng)及含間隙懸臂梁系統(tǒng)RFS[2]曲線圖。從圖中可以看出,傳統(tǒng)的RFS恢復(fù)力曲線法可有效解決單自由度系統(tǒng)非線性類別辨識(shí)問(wèn)題。但是,在連續(xù)體系統(tǒng)中,由于系統(tǒng)各單元間會(huì)發(fā)生進(jìn)一步相對(duì)位移,曲線呈嚴(yán)重的混疊現(xiàn)象,從原始的恢復(fù)力-位移曲線分布圖中無(wú)法準(zhǔn)確辨識(shí)出非線性類別。

    (a) 含間隙非線性單自由度 (b) 含間隙非線性懸臂梁系統(tǒng)RFS圖 系統(tǒng)RFS圖圖9 含間隙非線性系統(tǒng)RFS曲線

    3.2 懸臂梁系統(tǒng)非線性類別辨識(shí)

    根據(jù)哈密爾頓原理,圖7所示的非線性懸臂梁系統(tǒng)自由衰減振動(dòng)微分方程為

    (9)

    式中:utt是u(x,t)關(guān)于時(shí)間t的二階導(dǎo)數(shù);uxxxx是u(x,t)關(guān)于位置x的四階導(dǎo)數(shù);fs(x,t)為非線性力。

    根據(jù)系統(tǒng)邊界條件及振形疊加原理,系統(tǒng)固有頻率可寫(xiě)為

    (10)

    其中基頻為

    (11)

    利用本文方法,首先根據(jù)系統(tǒng)結(jié)構(gòu)特點(diǎn)及其動(dòng)力學(xué)方程,建立含剛度非線性系統(tǒng)簡(jiǎn)化模型,并求解模型自由衰減振動(dòng)響應(yīng),仿真參數(shù)見(jiàn)表5。

    采用峰值檢測(cè)法獲得上述不同剛度非線性系統(tǒng)瞬時(shí)頻率分布曲線。以基頻為基準(zhǔn)對(duì)各系統(tǒng)瞬時(shí)頻率分布進(jìn)行歸一化處理,并建立被測(cè)系統(tǒng)瞬時(shí)頻率曲線庫(kù),如圖10所示。

    從圖10中可以看出,上述剛度非線性系統(tǒng)瞬時(shí)頻率均呈不同的分布規(guī)律,尤其是立方剛度非線性系統(tǒng)E,其瞬時(shí)頻率呈光滑二次曲線,可直接從其他幾種非線性類別中區(qū)分出來(lái)。

    表5 懸臂梁非線性系統(tǒng)數(shù)據(jù)庫(kù)仿真參數(shù)

    (a)系統(tǒng)A (b)系統(tǒng)B (c)系統(tǒng)C

    (d)系統(tǒng)D (e)系統(tǒng)E(硬彈簧) (f)系統(tǒng)E(軟彈簧)圖10 剛度非線性懸臂梁系統(tǒng)瞬時(shí)頻率分布圖

    圖11 含間隙懸臂梁系統(tǒng)類別辨識(shí)結(jié)果

    以系統(tǒng)A為被測(cè)系統(tǒng),隨機(jī)選取50組不同的間隙與間隙剛度進(jìn)行仿真,圖11給出了被測(cè)系統(tǒng)實(shí)測(cè)瞬時(shí)頻率與曲線庫(kù)圖形匹配結(jié)果。從圖中可以看出,這50組不同間隙非線性系統(tǒng)均與曲線庫(kù)A的匹配程度最高,從而驗(yàn)證了該類別辨識(shí)方法的有效性。

    3.3 懸臂梁系統(tǒng)非線性類別辨識(shí)實(shí)驗(yàn)

    本文在假設(shè)系統(tǒng)剛度非線性未知的情況下,以系統(tǒng)A為實(shí)際被測(cè)系統(tǒng),通過(guò)改變系統(tǒng)間隙和間隙剛度驗(yàn)證該類別辨識(shí)方法的有效性。實(shí)驗(yàn)臺(tái)如圖12所示,間隙位于小懸臂梁與懸臂梁之間,通過(guò)小懸臂梁厚度調(diào)整間隙剛度,通過(guò)間隙頭與間隙梁之間的墊片厚度調(diào)整間隙大小。本實(shí)驗(yàn)用敲擊法獲得初始位移,并通過(guò)電渦流傳感器及數(shù)據(jù)采集儀進(jìn)行數(shù)據(jù)采集,具體實(shí)驗(yàn)參數(shù)見(jiàn)表6。

    (a)局部放大圖

    (b)主視圖圖12 非線性類別辨識(shí)實(shí)驗(yàn)臺(tái)

    實(shí)驗(yàn)采樣點(diǎn)數(shù)小懸臂梁厚度/mm間隙/mm1204880.12204880.23204880.34204840.15204860.162048100.1

    圖13給出了間隙為0.1 mm、小懸臂梁厚度為8 mm的衰減響應(yīng)曲線及其瞬時(shí)頻率分布曲線,從圖中可以看出該被測(cè)系統(tǒng)更接近系統(tǒng)A。

    分別將上述實(shí)驗(yàn)分析結(jié)果與懸臂梁非線性曲線庫(kù)進(jìn)行瞬時(shí)頻率圖形匹配,匹配結(jié)果見(jiàn)表7。從表7中可以看出,被測(cè)系統(tǒng)瞬時(shí)頻率與曲線庫(kù)A匹配的DS值最小,確定被測(cè)系統(tǒng)為間隙非線性系統(tǒng),從而驗(yàn)證了該方法的有效性。

    (a)非線性位置自由衰減振動(dòng) (b)瞬時(shí)頻率分布曲線 圖13 實(shí)驗(yàn)1結(jié)果

    實(shí)驗(yàn)DS系統(tǒng)A系統(tǒng)B系統(tǒng)C系統(tǒng)D110.488115.199115.32747.879230.826130.025130.34257.321316.916123.627123.89651.500430.787122.134122.64050.322521.246117.226117.66345.337624.074120.125120.54647.758

    4 結(jié) 論

    本文基于瞬時(shí)頻率峰值檢測(cè),利用被測(cè)系統(tǒng)瞬時(shí)頻率分布特點(diǎn),提出了一種可用于連續(xù)體系統(tǒng)的剛度非線性類別辨識(shí)方法。該方法根據(jù)被測(cè)系統(tǒng)結(jié)構(gòu)特點(diǎn),建立了含不同剛度非線性系統(tǒng)簡(jiǎn)化模型及其瞬時(shí)頻率曲線庫(kù),并結(jié)合圖形匹配算法實(shí)現(xiàn)了剛度非線性系統(tǒng)類別辨識(shí)。本文通過(guò)剛度非線性懸臂梁系統(tǒng)仿真分析與實(shí)驗(yàn)計(jì)算,驗(yàn)證了該方法的有效性。同時(shí),本文方法對(duì)其他非線性具有普適性,可推廣應(yīng)用于機(jī)械系統(tǒng)常見(jiàn)非線性的類別辨識(shí)中。

    [1] ADAMS D E, ALLEMANG R J. Characterization of nonlinear vibrating systems using internal feedback and frequency response modulation [J]. ASME Journal of Vibration and Acoustics, 1999, 121(4): 495-500.

    [2] MALATKAR P, NAYFEH A H. A plethora of nonlinear dynamics phenomena observed in a simple cantilever plate [C]∥ASME 2003 International Design Engineering Technical Conferences and Computers and Information in Engineering Conference. New York, USA: ASME, 2003: 2603-2610.

    [3] GAЁTAN KERSCHEN, WORDEN K, VAKAKIS A F, et al. Past, present and future of nonlinear system identification in structural dynamics [J]. Mechanical Systems & Signal Processing, 2006, 20(3): 505-592.

    [4] VAKAKIS A F, EWINS D J. Effects of weak non-linearities on modal analysis [J]. Mechanical Systems & Signal Processing, 1994, 8(2): 175-198.

    [5] CHATTERJEE A, VYAS N S. Stiffness non-linearity classification through structured response component analysis using Volterra series [J]. Mechanical Systems & Signal Processing, 2001, 15(2): 323-336.

    [6] BENDAT J S, PALO P A, COPPOLINO R N. A general identification technique for nonlinear differential equations of motion [J]. Probabilistic Engineering Mechanics, 1992, 92(7): 43-61.

    [7] GONDHALEKAR A C, PETROV E P, IMREGUN M. Parameters identification for nonlinear dynamic systems via genetic algorithm optimization [J]. Journal of Computational & Nonlinear Dynamics, 2009, 4(4): 1724-1732.

    [8] BALESTRINO A, CAITI A, CRISOSTOMI E. A classification of nonlinear systems: an entropy based approach [EB/OL].[2014-09-10]. http: ∥www. nt.ntnu.no/users/skoge/prost/proceedings/icheap8-pres07/icheap8webpapers/193%20Balestrino.pdf.

    [9] GALLEANI L, PRESTI L L, STEFANO A D. A method for nonlinear system classification in the time-frequency plane [J]. Signal Processing, 1998, 65(1): 147-153.

    [10]STROGATS Z, FRIEDMAN M, MALLINCKRODT A J. Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering [M]. New York, USA: Perseus Books Publishing, 1994: 216-227.

    [11]COHEN L. Time-frequency analysis [M]. Englewood Cliffs, NJ, USA: Prentice Hall, 1995: 98-153.

    [12]張賢達(dá), 保錚. 平穩(wěn)信號(hào)分析與處理 [M]. 北京: 國(guó)防工業(yè)出版社, 1998: 53-68, 84-92, 392-422.

    (編輯 杜秀杰)

    Type Identification for Mechanical Nonlinearity

    QUAN Shuanglu,LIU Jie,HAN Luofeng,HE Anqi,LI Bing

    (School of Mechanical Engineering, Xi’an Jiaotong University, Xi’an 710049, China)

    For type identification of nonlinear continuum systems with stiffness nonlinearity, a classification strategy with peak detection algorithm is put forward. According to the structure of under test system, simplified models with 5 kinds of stiffness nonlinearity are constructed. The nonlinearity types are particularly interested in structural diagnostic, such as clearance nonlinearity and cubic nonlinearity. An instantaneous frequency library is then set up by studying the free oscillations and instantaneous frequencies of different models. The shape-matching algorithm is used to compare the real instantaneous frequency with the curves in the library. And the results are taken up to the index to identify the nonlinearity types. For a cantilever beam system with clearance, the matching value between the real instantaneous frequency and the curve of clearance nonlinearity measured in this strategy gets obviously smaller than the others, verifying the effectiveness.

    mechanical system; type of nonlinear system; instantaneous frequency; shape-matching algorithm

    2015-04-17。

    權(quán)雙璐(1991—),女,碩士生;李兵(通信作者),男,教授,博士生導(dǎo)師。

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51475356)。

    時(shí)間:2015-09-13

    10.7652/xjtuxb201512017

    TH17

    A

    0253-987X(2015)12-0104-08

    網(wǎng)絡(luò)出版地址:http:∥www.cnki.net/kcms/detail/61.1069.T.20150913.1824.006.html

    猜你喜歡
    測(cè)系統(tǒng)機(jī)械系統(tǒng)類別
    油液混合動(dòng)力工程機(jī)械系統(tǒng)及控制策略研究綜述
    基于隱馬爾可夫模型的航空機(jī)械系統(tǒng)故障診斷算法設(shè)計(jì)
    防洪非工程措施設(shè)計(jì)實(shí)例——嘉興市水文巡測(cè)系統(tǒng)項(xiàng)目設(shè)計(jì)
    電快速瞬變脈沖群對(duì)核測(cè)系統(tǒng)的影響及對(duì)策
    服務(wù)類別
    運(yùn)行中的機(jī)械系統(tǒng)的動(dòng)態(tài)熱分析
    基于廣域量測(cè)系統(tǒng)的電力系統(tǒng)綜合負(fù)荷辨識(shí)模型的研究
    三維標(biāo)測(cè)系統(tǒng)指導(dǎo)下射頻消融治療房顫的護(hù)理觀察
    論類別股東會(huì)
    商事法論集(2014年1期)2014-06-27 01:20:42
    中醫(yī)類別全科醫(yī)師培養(yǎng)模式的探討
    亚洲av五月六月丁香网| 女生性感内裤真人,穿戴方法视频| 日韩大尺度精品在线看网址 | 久久久久国内视频| 两个人免费观看高清视频| 伊人久久大香线蕉亚洲五| 国产亚洲av高清不卡| 99国产精品99久久久久| 人妻丰满熟妇av一区二区三区| 天堂影院成人在线观看| 欧美日韩福利视频一区二区| 国产精品久久久久久人妻精品电影| 一个人免费在线观看的高清视频| 88av欧美| 精品国产一区二区三区四区第35| 两性午夜刺激爽爽歪歪视频在线观看 | 超碰成人久久| 看黄色毛片网站| 国产在线精品亚洲第一网站| 欧美日韩一级在线毛片| 国内毛片毛片毛片毛片毛片| 欧美日本中文国产一区发布| 国产主播在线观看一区二区| 91字幕亚洲| 美女国产高潮福利片在线看| 久久久久精品国产欧美久久久| 可以在线观看的亚洲视频| 久久九九热精品免费| tocl精华| 9191精品国产免费久久| 一边摸一边抽搐一进一小说| 免费av毛片视频| 在线观看一区二区三区| 好男人电影高清在线观看| 久久香蕉精品热| 岛国在线观看网站| 午夜福利影视在线免费观看| 老司机深夜福利视频在线观看| 亚洲成av人片免费观看| 两个人免费观看高清视频| 欧美色欧美亚洲另类二区 | 亚洲一区二区三区不卡视频| 午夜精品国产一区二区电影| 国产一区二区在线av高清观看| 18美女黄网站色大片免费观看| 亚洲成a人片在线一区二区| 可以在线观看毛片的网站| 国产精品免费视频内射| 亚洲精品在线观看二区| 亚洲精品粉嫩美女一区| 亚洲精华国产精华精| 人妻久久中文字幕网| 搡老熟女国产l中国老女人| 91字幕亚洲| 99在线人妻在线中文字幕| 亚洲五月婷婷丁香| 脱女人内裤的视频| 亚洲 国产 在线| 免费高清视频大片| 国产欧美日韩综合在线一区二区| 国产aⅴ精品一区二区三区波| 手机成人av网站| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品国产亚洲av高清一级| 亚洲人成伊人成综合网2020| 91在线观看av| 久久国产精品男人的天堂亚洲| 99国产综合亚洲精品| av福利片在线| 波多野结衣巨乳人妻| 不卡一级毛片| 免费高清视频大片| 亚洲一卡2卡3卡4卡5卡精品中文| 免费一级毛片在线播放高清视频 | 99精品久久久久人妻精品| 久久影院123| 在线观看66精品国产| 精品一品国产午夜福利视频| 免费看十八禁软件| 亚洲情色 制服丝袜| 亚洲一区高清亚洲精品| 久久久久久国产a免费观看| 亚洲午夜精品一区,二区,三区| 高清毛片免费观看视频网站| 超碰成人久久| 久久国产精品人妻蜜桃| 欧美一级毛片孕妇| 精品无人区乱码1区二区| 青草久久国产| 久久精品成人免费网站| 国产成人精品久久二区二区91| 国产一区二区三区综合在线观看| 一级a爱片免费观看的视频| 大陆偷拍与自拍| aaaaa片日本免费| 一级毛片高清免费大全| 久久欧美精品欧美久久欧美| 俄罗斯特黄特色一大片| 熟女少妇亚洲综合色aaa.| 少妇熟女aⅴ在线视频| 波多野结衣av一区二区av| 国产一区二区三区在线臀色熟女| 国产成人欧美在线观看| 岛国在线观看网站| 色综合亚洲欧美另类图片| 久久精品91无色码中文字幕| 亚洲情色 制服丝袜| 精品国产亚洲在线| 女生性感内裤真人,穿戴方法视频| 真人一进一出gif抽搐免费| 国内精品久久久久精免费| 国产99白浆流出| 深夜精品福利| 黄色片一级片一级黄色片| 免费一级毛片在线播放高清视频 | 两人在一起打扑克的视频| 黑人操中国人逼视频| 国产av在哪里看| 亚洲狠狠婷婷综合久久图片| 久久久久国产一级毛片高清牌| 正在播放国产对白刺激| 好男人在线观看高清免费视频 | 黄片大片在线免费观看| 久久久水蜜桃国产精品网| 宅男免费午夜| 波多野结衣av一区二区av| 麻豆成人av在线观看| 在线av久久热| 美女高潮到喷水免费观看| av在线播放免费不卡| 日本五十路高清| 日本免费一区二区三区高清不卡 | 亚洲精品国产精品久久久不卡| 亚洲精品久久成人aⅴ小说| 国产在线观看jvid| 亚洲国产精品999在线| 亚洲精品在线观看二区| 色在线成人网| 国产野战对白在线观看| 美女大奶头视频| cao死你这个sao货| 最近最新中文字幕大全电影3 | 国产亚洲av高清不卡| 欧美激情 高清一区二区三区| 一区二区三区精品91| 欧美最黄视频在线播放免费| 免费看a级黄色片| 免费观看人在逋| 人妻久久中文字幕网| 国语自产精品视频在线第100页| 国产高清激情床上av| 99riav亚洲国产免费| 欧美久久黑人一区二区| 女人高潮潮喷娇喘18禁视频| 最近最新免费中文字幕在线| 日日爽夜夜爽网站| 亚洲专区中文字幕在线| 亚洲七黄色美女视频| 亚洲男人的天堂狠狠| 精品不卡国产一区二区三区| 国产三级黄色录像| 国产亚洲欧美精品永久| 精品一区二区三区视频在线观看免费| 女人高潮潮喷娇喘18禁视频| 色哟哟哟哟哟哟| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕另类日韩欧美亚洲嫩草| 精品人妻1区二区| 久久欧美精品欧美久久欧美| 黑人巨大精品欧美一区二区蜜桃| 国产精品影院久久| 97碰自拍视频| 免费少妇av软件| 男人舔女人下体高潮全视频| 亚洲成人国产一区在线观看| 欧美黄色片欧美黄色片| 久久久久久久午夜电影| 国产精品精品国产色婷婷| 国产亚洲精品久久久久5区| 国内毛片毛片毛片毛片毛片| 欧美另类亚洲清纯唯美| 999久久久国产精品视频| 制服丝袜大香蕉在线| 亚洲精品国产一区二区精华液| 国产亚洲精品久久久久久毛片| 日韩三级视频一区二区三区| 欧美丝袜亚洲另类 | 久久精品91蜜桃| 久久久国产成人精品二区| 国产成人av激情在线播放| 欧美激情久久久久久爽电影 | 精品欧美一区二区三区在线| 91精品三级在线观看| 99国产精品一区二区蜜桃av| 欧美黄色淫秽网站| 精品国产国语对白av| 男人舔女人的私密视频| 亚洲av成人av| 亚洲精品久久国产高清桃花| 这个男人来自地球电影免费观看| 亚洲人成电影免费在线| 欧美大码av| 一本综合久久免费| 国产三级在线视频| 亚洲av电影在线进入| 精品熟女少妇八av免费久了| 久久国产精品人妻蜜桃| 一边摸一边做爽爽视频免费| 午夜福利,免费看| 国产亚洲av嫩草精品影院| 91麻豆av在线| 淫妇啪啪啪对白视频| 色播在线永久视频| 午夜精品国产一区二区电影| 中文字幕人妻丝袜一区二区| 大香蕉久久成人网| 中文字幕最新亚洲高清| 一区二区三区激情视频| 亚洲一码二码三码区别大吗| 91在线观看av| 深夜精品福利| 色综合欧美亚洲国产小说| 欧美亚洲日本最大视频资源| 欧美成人性av电影在线观看| netflix在线观看网站| 91麻豆精品激情在线观看国产| 国产精品秋霞免费鲁丝片| 国产精品98久久久久久宅男小说| 亚洲av日韩精品久久久久久密| 国产精品秋霞免费鲁丝片| 国产精品亚洲美女久久久| 亚洲情色 制服丝袜| 99香蕉大伊视频| 在线国产一区二区在线| 日本五十路高清| 国产一区二区激情短视频| 麻豆国产av国片精品| 成人欧美大片| 日本五十路高清| 亚洲狠狠婷婷综合久久图片| 日本一区二区免费在线视频| 超碰成人久久| 亚洲自拍偷在线| 老鸭窝网址在线观看| 国产av一区在线观看免费| 一边摸一边抽搐一进一出视频| 亚洲国产毛片av蜜桃av| 久久天躁狠狠躁夜夜2o2o| 色综合欧美亚洲国产小说| 不卡av一区二区三区| 一区二区日韩欧美中文字幕| 午夜福利视频1000在线观看 | 在线观看舔阴道视频| 人人妻,人人澡人人爽秒播| 99国产精品免费福利视频| 国产精华一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 欧美不卡视频在线免费观看 | 99精品在免费线老司机午夜| 99久久久亚洲精品蜜臀av| 两个人视频免费观看高清| 我的亚洲天堂| 男女做爰动态图高潮gif福利片 | 亚洲中文字幕一区二区三区有码在线看 | 性色av乱码一区二区三区2| 黄色a级毛片大全视频| 在线观看一区二区三区| 99久久精品国产亚洲精品| 国产一级毛片七仙女欲春2 | www.精华液| 国产精品久久久久久亚洲av鲁大| 99精品久久久久人妻精品| 少妇 在线观看| 国内毛片毛片毛片毛片毛片| 国产精品野战在线观看| 美女高潮喷水抽搐中文字幕| 多毛熟女@视频| 欧美日本亚洲视频在线播放| 美女 人体艺术 gogo| 夜夜看夜夜爽夜夜摸| 满18在线观看网站| 精品高清国产在线一区| 中文字幕人妻丝袜一区二区| 国产精品一区二区免费欧美| 精品人妻在线不人妻| 最新在线观看一区二区三区| 日韩免费av在线播放| 久久久久九九精品影院| 变态另类成人亚洲欧美熟女 | 亚洲欧美日韩另类电影网站| 18禁美女被吸乳视频| 真人做人爱边吃奶动态| 亚洲精品美女久久久久99蜜臀| www.熟女人妻精品国产| 久久国产乱子伦精品免费另类| 波多野结衣一区麻豆| 一级a爱视频在线免费观看| 日本免费a在线| 欧美最黄视频在线播放免费| 亚洲激情在线av| 成人三级做爰电影| 国产成人系列免费观看| 亚洲人成伊人成综合网2020| 亚洲av成人不卡在线观看播放网| 日本一区二区免费在线视频| 久久久久久人人人人人| 69精品国产乱码久久久| 老司机靠b影院| 18禁裸乳无遮挡免费网站照片 | 久久香蕉精品热| 男女下面插进去视频免费观看| 午夜日韩欧美国产| 99在线人妻在线中文字幕| 国产aⅴ精品一区二区三区波| 在线观看一区二区三区| 欧美中文日本在线观看视频| 国产男靠女视频免费网站| 69精品国产乱码久久久| 一级a爱视频在线免费观看| 99精品欧美一区二区三区四区| 日本 欧美在线| 国产一区二区在线av高清观看| tocl精华| 亚洲午夜理论影院| 在线视频色国产色| 亚洲色图综合在线观看| av天堂在线播放| 日本 欧美在线| 99久久99久久久精品蜜桃| 一二三四社区在线视频社区8| 中文字幕最新亚洲高清| 亚洲激情在线av| 亚洲国产精品久久男人天堂| 一级a爱片免费观看的视频| av片东京热男人的天堂| 国产亚洲精品第一综合不卡| 97超级碰碰碰精品色视频在线观看| 99久久久亚洲精品蜜臀av| 91国产中文字幕| 欧美成人午夜精品| 亚洲精华国产精华精| 欧美不卡视频在线免费观看 | 中文字幕人成人乱码亚洲影| 人妻久久中文字幕网| 午夜福利18| 国产成人av教育| 久久久久久人人人人人| 18禁黄网站禁片午夜丰满| 每晚都被弄得嗷嗷叫到高潮| 国产精品美女特级片免费视频播放器 | 亚洲av日韩精品久久久久久密| 精品国内亚洲2022精品成人| 长腿黑丝高跟| 国产成人av激情在线播放| 精品乱码久久久久久99久播| 夜夜躁狠狠躁天天躁| 欧美黄色片欧美黄色片| 高潮久久久久久久久久久不卡| 日韩欧美三级三区| 人妻丰满熟妇av一区二区三区| 岛国视频午夜一区免费看| 欧美大码av| 亚洲激情在线av| 国产高清videossex| 国产av一区二区精品久久| www.精华液| 亚洲精品久久成人aⅴ小说| 色综合欧美亚洲国产小说| 精品久久久久久,| 欧美黄色淫秽网站| 免费少妇av软件| 狂野欧美激情性xxxx| 在线视频色国产色| 一卡2卡三卡四卡精品乱码亚洲| 精品国产一区二区三区四区第35| 日韩欧美在线二视频| 亚洲片人在线观看| 国产精品一区二区在线不卡| 1024香蕉在线观看| 久久国产精品人妻蜜桃| 亚洲av美国av| 91老司机精品| 亚洲熟妇中文字幕五十中出| 日韩大尺度精品在线看网址 | 女人爽到高潮嗷嗷叫在线视频| 日韩视频一区二区在线观看| 无限看片的www在线观看| 黄色视频,在线免费观看| av电影中文网址| 中文字幕另类日韩欧美亚洲嫩草| 熟妇人妻久久中文字幕3abv| 日日干狠狠操夜夜爽| 岛国在线观看网站| 亚洲国产欧美一区二区综合| 九色亚洲精品在线播放| 一本久久中文字幕| 久久久久久免费高清国产稀缺| 香蕉国产在线看| 丝袜在线中文字幕| 久99久视频精品免费| 91字幕亚洲| 国产精品秋霞免费鲁丝片| 99国产极品粉嫩在线观看| 十分钟在线观看高清视频www| 在线视频色国产色| 午夜久久久在线观看| 久久精品aⅴ一区二区三区四区| 校园春色视频在线观看| 国产91精品成人一区二区三区| 最近最新中文字幕大全电影3 | 亚洲精品粉嫩美女一区| 99久久99久久久精品蜜桃| 久久中文字幕一级| 欧美性长视频在线观看| av超薄肉色丝袜交足视频| 欧美成人午夜精品| 波多野结衣一区麻豆| 精品电影一区二区在线| 夜夜爽天天搞| 老司机在亚洲福利影院| 午夜免费激情av| 亚洲精品中文字幕在线视频| 99久久国产精品久久久| 色在线成人网| 中文字幕色久视频| 国产蜜桃级精品一区二区三区| 欧美大码av| 亚洲久久久国产精品| 亚洲色图综合在线观看| 亚洲九九香蕉| 成人国产综合亚洲| 男人的好看免费观看在线视频 | 无限看片的www在线观看| 亚洲av熟女| 国产精品免费视频内射| 国产av又大| 嫩草影视91久久| 国产精品亚洲av一区麻豆| 中文字幕人成人乱码亚洲影| 亚洲精品国产色婷婷电影| 国产精品二区激情视频| 亚洲av片天天在线观看| 亚洲av熟女| 国产视频一区二区在线看| 国产亚洲av嫩草精品影院| 老熟妇仑乱视频hdxx| aaaaa片日本免费| 免费少妇av软件| 国产精品永久免费网站| www.999成人在线观看| 日日夜夜操网爽| 叶爱在线成人免费视频播放| 9热在线视频观看99| 久久久久九九精品影院| 亚洲全国av大片| 欧美中文日本在线观看视频| 国产主播在线观看一区二区| 亚洲欧美日韩高清在线视频| 成人免费观看视频高清| 老司机深夜福利视频在线观看| 免费人成视频x8x8入口观看| 在线观看www视频免费| 久久香蕉激情| 国产成人系列免费观看| 天堂影院成人在线观看| av有码第一页| 不卡av一区二区三区| 在线观看午夜福利视频| 天天躁夜夜躁狠狠躁躁| 91麻豆av在线| 88av欧美| 在线观看舔阴道视频| 色播亚洲综合网| 激情视频va一区二区三区| 亚洲人成电影免费在线| 亚洲国产精品sss在线观看| 多毛熟女@视频| 91精品三级在线观看| 久久国产亚洲av麻豆专区| 久久人妻熟女aⅴ| 午夜激情av网站| 真人一进一出gif抽搐免费| 99国产极品粉嫩在线观看| 成熟少妇高潮喷水视频| 国产精品一区二区三区四区久久 | 熟妇人妻久久中文字幕3abv| a级毛片在线看网站| 国产99白浆流出| 亚洲人成伊人成综合网2020| 男女下面进入的视频免费午夜 | 免费一级毛片在线播放高清视频 | 欧美成人午夜精品| 久久久精品国产亚洲av高清涩受| 美女午夜性视频免费| 人人澡人人妻人| 天天添夜夜摸| 变态另类成人亚洲欧美熟女 | 一级a爱片免费观看的视频| 亚洲天堂国产精品一区在线| 亚洲人成77777在线视频| 一边摸一边抽搐一进一出视频| 香蕉国产在线看| 村上凉子中文字幕在线| 国产日韩一区二区三区精品不卡| 大陆偷拍与自拍| 亚洲专区国产一区二区| 国产极品粉嫩免费观看在线| 啦啦啦免费观看视频1| 黄色a级毛片大全视频| 91麻豆av在线| 亚洲成av片中文字幕在线观看| 香蕉久久夜色| 亚洲人成电影观看| 亚洲伊人色综图| 国产精品久久电影中文字幕| 国产欧美日韩精品亚洲av| 成人永久免费在线观看视频| a在线观看视频网站| 神马国产精品三级电影在线观看 | 人人妻人人澡人人看| 亚洲国产精品999在线| 九色国产91popny在线| 国产精品自产拍在线观看55亚洲| 亚洲专区国产一区二区| 久99久视频精品免费| av在线播放免费不卡| 久久青草综合色| 亚洲欧美一区二区三区黑人| 激情视频va一区二区三区| 十分钟在线观看高清视频www| 91麻豆精品激情在线观看国产| 丝袜美足系列| 久久久久国内视频| 亚洲少妇的诱惑av| 久久午夜综合久久蜜桃| 国产精品av久久久久免费| 亚洲精品一卡2卡三卡4卡5卡| 一本久久中文字幕| 欧美日韩黄片免| 欧美日韩一级在线毛片| 亚洲最大成人中文| 亚洲国产看品久久| 久久久久久久久久久久大奶| 欧美精品亚洲一区二区| www.www免费av| 国产精品国产高清国产av| 久久久久久大精品| 在线国产一区二区在线| 久久人妻福利社区极品人妻图片| 色婷婷久久久亚洲欧美| 亚洲精华国产精华精| 久久香蕉精品热| 久久性视频一级片| 午夜两性在线视频| 一级作爱视频免费观看| 最好的美女福利视频网| 国产熟女午夜一区二区三区| 午夜福利高清视频| 久久久精品国产亚洲av高清涩受| 99久久久亚洲精品蜜臀av| 好男人在线观看高清免费视频 | 夜夜看夜夜爽夜夜摸| 嫁个100分男人电影在线观看| 日本黄色视频三级网站网址| 日本一区二区免费在线视频| 久久久久久人人人人人| 亚洲色图 男人天堂 中文字幕| 亚洲全国av大片| 欧美中文日本在线观看视频| 搡老岳熟女国产| 成人三级做爰电影| 黄色女人牲交| 日韩大尺度精品在线看网址 | 免费在线观看日本一区| 精品乱码久久久久久99久播| 波多野结衣巨乳人妻| 丝袜美腿诱惑在线| 欧美日韩瑟瑟在线播放| 免费女性裸体啪啪无遮挡网站| 国产精品国产高清国产av| 在线观看午夜福利视频| 亚洲精品国产区一区二| 午夜免费激情av| 黑人巨大精品欧美一区二区蜜桃| 最新在线观看一区二区三区| 无遮挡黄片免费观看| 国产激情久久老熟女| 亚洲久久久国产精品| 真人做人爱边吃奶动态| 搡老妇女老女人老熟妇| 老熟妇乱子伦视频在线观看| 亚洲五月天丁香| 久久久久久人人人人人| 老司机福利观看| 国产高清激情床上av| 国产精品一区二区在线不卡| 天堂√8在线中文| 精品高清国产在线一区| 99国产精品一区二区蜜桃av| 亚洲 欧美 日韩 在线 免费| 热99re8久久精品国产| 午夜福利18| 黄色丝袜av网址大全| 免费高清视频大片| 国产激情欧美一区二区| 日本a在线网址| 在线观看www视频免费| 一进一出抽搐动态| 亚洲美女黄片视频| 国内久久婷婷六月综合欲色啪| 18禁国产床啪视频网站| 脱女人内裤的视频| 99国产综合亚洲精品| 免费无遮挡裸体视频| 国产精品久久视频播放|