• <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)模式的探討
    国产97色在线日韩免费| 手机成人av网站| 精品人妻1区二区| 熟女av电影| 女性生殖器流出的白浆| 国产精品偷伦视频观看了| 可以免费在线观看a视频的电影网站| 欧美日韩视频高清一区二区三区二| 2021少妇久久久久久久久久久| 欧美日韩黄片免| 亚洲国产看品久久| 欧美日韩综合久久久久久| 中文字幕人妻熟女乱码| 人人澡人人妻人| 9热在线视频观看99| 丝袜喷水一区| 捣出白浆h1v1| 十八禁人妻一区二区| 亚洲精品国产一区二区精华液| 婷婷色综合www| 成人国产一区最新在线观看 | 黄色怎么调成土黄色| 侵犯人妻中文字幕一二三四区| 搡老岳熟女国产| 国产高清不卡午夜福利| 日本一区二区免费在线视频| 高潮久久久久久久久久久不卡| 精品少妇一区二区三区视频日本电影| 欧美黑人精品巨大| 成年美女黄网站色视频大全免费| 99热网站在线观看| 99国产精品99久久久久| 国产高清视频在线播放一区 | 嫁个100分男人电影在线观看 | 啦啦啦 在线观看视频| 久久性视频一级片| 亚洲精品在线美女| 免费少妇av软件| 亚洲九九香蕉| 精品国产乱码久久久久久小说| 黄色a级毛片大全视频| 一边摸一边做爽爽视频免费| h视频一区二区三区| 51午夜福利影视在线观看| av片东京热男人的天堂| 亚洲av片天天在线观看| www.自偷自拍.com| 国产一区二区激情短视频 | xxxhd国产人妻xxx| 日韩中文字幕欧美一区二区 | 国产成人精品久久二区二区免费| 少妇粗大呻吟视频| 十八禁网站网址无遮挡| 日本色播在线视频| 校园人妻丝袜中文字幕| 丝袜在线中文字幕| 亚洲,欧美,日韩| 国产黄色免费在线视频| 丰满人妻熟妇乱又伦精品不卡| a级毛片黄视频| 欧美97在线视频| 久久99一区二区三区| 国产一区亚洲一区在线观看| 亚洲成人手机| 亚洲欧美激情在线| 首页视频小说图片口味搜索 | 婷婷色综合www| 国产黄色视频一区二区在线观看| 欧美亚洲日本最大视频资源| 亚洲精品日韩在线中文字幕| 国产欧美日韩综合在线一区二区| 欧美大码av| 久久ye,这里只有精品| 中文字幕av电影在线播放| 亚洲欧美清纯卡通| 男人舔女人的私密视频| av欧美777| 91精品三级在线观看| 超色免费av| 国产成人91sexporn| 中文字幕最新亚洲高清| 国产免费视频播放在线视频| 99九九在线精品视频| 999精品在线视频| 日本91视频免费播放| 人妻 亚洲 视频| 国产免费福利视频在线观看| 国产女主播在线喷水免费视频网站| 99热网站在线观看| 欧美大码av| 一区二区三区乱码不卡18| 亚洲国产最新在线播放| 999久久久国产精品视频| tube8黄色片| 亚洲男人天堂网一区| 叶爱在线成人免费视频播放| 国产免费福利视频在线观看| 国产熟女午夜一区二区三区| 色综合欧美亚洲国产小说| 午夜91福利影院| 自线自在国产av| 精品人妻在线不人妻| 在线天堂中文资源库| 超色免费av| 香蕉国产在线看| 国产97色在线日韩免费| 欧美乱码精品一区二区三区| 中文字幕亚洲精品专区| 久久国产精品男人的天堂亚洲| 亚洲久久久国产精品| 狠狠婷婷综合久久久久久88av| 国产成人av教育| 欧美少妇被猛烈插入视频| 99国产精品一区二区三区| 国产一级毛片在线| 欧美av亚洲av综合av国产av| 成人亚洲欧美一区二区av| 香蕉国产在线看| 日韩伦理黄色片| 亚洲 欧美一区二区三区| 亚洲中文日韩欧美视频| 欧美日韩av久久| 男人舔女人的私密视频| 亚洲第一青青草原| 香蕉丝袜av| 后天国语完整版免费观看| 国产精品三级大全| 一级毛片我不卡| 婷婷成人精品国产| 欧美黑人精品巨大| 免费观看av网站的网址| 国产极品粉嫩免费观看在线| 国产高清国产精品国产三级| 超碰成人久久| 午夜免费观看性视频| 日韩大码丰满熟妇| 国产高清国产精品国产三级| 国产高清视频在线播放一区 | 黄色片一级片一级黄色片| 一区二区三区乱码不卡18| 成年动漫av网址| 最新的欧美精品一区二区| 亚洲欧美日韩高清在线视频 | 午夜两性在线视频| 亚洲成av片中文字幕在线观看| 国产欧美日韩综合在线一区二区| 夜夜骑夜夜射夜夜干| 在线观看免费视频网站a站| 亚洲精品久久久久久婷婷小说| 中文字幕av电影在线播放| 丝袜喷水一区| 亚洲激情五月婷婷啪啪| 久久精品久久精品一区二区三区| 丝瓜视频免费看黄片| 男女之事视频高清在线观看 | tube8黄色片| 久久天躁狠狠躁夜夜2o2o | 丝袜美腿诱惑在线| 亚洲综合色网址| 精品少妇一区二区三区视频日本电影| 亚洲精品美女久久av网站| 亚洲伊人久久精品综合| a级毛片黄视频| 中文字幕人妻丝袜一区二区| 久久久精品国产亚洲av高清涩受| 亚洲国产成人一精品久久久| 涩涩av久久男人的天堂| 亚洲国产欧美一区二区综合| 精品少妇久久久久久888优播| 女人高潮潮喷娇喘18禁视频| 国产女主播在线喷水免费视频网站| 9色porny在线观看| 女人高潮潮喷娇喘18禁视频| 久久天躁狠狠躁夜夜2o2o | 51午夜福利影视在线观看| 成人影院久久| 国产欧美日韩一区二区三区在线| 中文字幕另类日韩欧美亚洲嫩草| 男女国产视频网站| 纯流量卡能插随身wifi吗| 又大又黄又爽视频免费| 国产亚洲欧美精品永久| xxxhd国产人妻xxx| 亚洲五月婷婷丁香| 各种免费的搞黄视频| 欧美日韩国产mv在线观看视频| 丝袜在线中文字幕| 亚洲精品国产av蜜桃| 久久人人97超碰香蕉20202| 久久久国产精品麻豆| 国产欧美日韩综合在线一区二区| 国产一区二区 视频在线| av不卡在线播放| 在线亚洲精品国产二区图片欧美| 亚洲,欧美,日韩| 两性夫妻黄色片| 色视频在线一区二区三区| 热re99久久精品国产66热6| 国产一卡二卡三卡精品| 男女免费视频国产| 色婷婷av一区二区三区视频| av国产精品久久久久影院| 久久久久久人人人人人| 国产精品国产av在线观看| 肉色欧美久久久久久久蜜桃| 99国产精品一区二区蜜桃av | 免费在线观看视频国产中文字幕亚洲 | 极品少妇高潮喷水抽搐| 啦啦啦 在线观看视频| 嫁个100分男人电影在线观看 | 三上悠亚av全集在线观看| 亚洲,欧美精品.| 亚洲七黄色美女视频| 七月丁香在线播放| 亚洲色图综合在线观看| 亚洲欧美精品综合一区二区三区| 女人高潮潮喷娇喘18禁视频| 国产熟女欧美一区二区| 国产精品一区二区免费欧美 | 精品免费久久久久久久清纯 | 国产精品99久久99久久久不卡| 亚洲男人天堂网一区| 日韩伦理黄色片| 热99国产精品久久久久久7| 成年动漫av网址| 久久99精品国语久久久| 久久精品国产综合久久久| 美女中出高潮动态图| 汤姆久久久久久久影院中文字幕| 亚洲精品久久午夜乱码| 久久久国产欧美日韩av| 母亲3免费完整高清在线观看| 人人澡人人妻人| 国产日韩欧美亚洲二区| 久久天躁狠狠躁夜夜2o2o | 国产xxxxx性猛交| 在线亚洲精品国产二区图片欧美| 亚洲精品一卡2卡三卡4卡5卡 | 男女床上黄色一级片免费看| 欧美黄色片欧美黄色片| 1024视频免费在线观看| 麻豆av在线久日| 国产亚洲欧美精品永久| 91九色精品人成在线观看| 黄网站色视频无遮挡免费观看| 亚洲视频免费观看视频| 久久性视频一级片| 国产精品熟女久久久久浪| 欧美日本中文国产一区发布| 国产成人啪精品午夜网站| 亚洲专区中文字幕在线| 汤姆久久久久久久影院中文字幕| 国产精品久久久久成人av| www.999成人在线观看| 高清av免费在线| 99久久人妻综合| 国产成人影院久久av| 国产在视频线精品| 黄色视频不卡| 精品一区在线观看国产| 午夜福利视频在线观看免费| 波多野结衣一区麻豆| 97人妻天天添夜夜摸| 美女脱内裤让男人舔精品视频| 色视频在线一区二区三区| 日韩精品免费视频一区二区三区| 性少妇av在线| 每晚都被弄得嗷嗷叫到高潮| 亚洲一区二区三区欧美精品| 不卡av一区二区三区| 女人久久www免费人成看片| 久久久精品94久久精品| 丁香六月天网| 国产欧美日韩综合在线一区二区| 一边摸一边做爽爽视频免费| 国产成人精品无人区| 在线观看免费高清a一片| 国产不卡av网站在线观看| 一本—道久久a久久精品蜜桃钙片| 国产无遮挡羞羞视频在线观看| 国产有黄有色有爽视频| 亚洲精品久久午夜乱码| 午夜免费男女啪啪视频观看| xxx大片免费视频| 好男人电影高清在线观看| 国产精品国产三级专区第一集| 老司机深夜福利视频在线观看 | 亚洲精品日韩在线中文字幕| 欧美国产精品va在线观看不卡| 五月开心婷婷网| 欧美大码av| avwww免费| 欧美黄色淫秽网站| 狠狠精品人妻久久久久久综合| 婷婷色麻豆天堂久久| 日本欧美国产在线视频| 深夜精品福利| 精品国产乱码久久久久久小说| 久久人人爽人人片av| 久久久国产精品麻豆| 亚洲成国产人片在线观看| 久久久久久人人人人人| 国产成人精品久久二区二区91| 搡老乐熟女国产| 成年动漫av网址| 国产精品亚洲av一区麻豆| 久久久久久人人人人人| 中文欧美无线码| 不卡av一区二区三区| 在线天堂中文资源库| 国产一区二区 视频在线| 欧美黑人精品巨大| 国产男人的电影天堂91| 国产亚洲av片在线观看秒播厂| 18禁黄网站禁片午夜丰满| 欧美亚洲日本最大视频资源| 超色免费av| 日韩 亚洲 欧美在线| 久久精品国产亚洲av高清一级| 精品少妇内射三级| 日本黄色日本黄色录像| www.精华液| 一个人免费看片子| 极品人妻少妇av视频| 狠狠精品人妻久久久久久综合| 超色免费av| www日本在线高清视频| 亚洲,欧美精品.| 亚洲欧洲日产国产| 美女中出高潮动态图| 国产成人欧美在线观看 | 18禁裸乳无遮挡动漫免费视频| 亚洲国产欧美日韩在线播放| 永久免费av网站大全| 咕卡用的链子| 亚洲人成77777在线视频| 久久国产精品人妻蜜桃| 午夜福利一区二区在线看| 人妻久久中文字幕网| 亚洲真实伦在线观看| 成年人黄色毛片网站| 精品福利观看| 一个人免费在线观看的高清视频| 久久 成人 亚洲| 亚洲熟妇中文字幕五十中出| 国产aⅴ精品一区二区三区波| 黄色女人牲交| 丰满人妻熟妇乱又伦精品不卡| 国内揄拍国产精品人妻在线 | 亚洲专区中文字幕在线| 麻豆国产av国片精品| 欧美三级亚洲精品| 国产成人影院久久av| 麻豆成人av在线观看| 欧美在线一区亚洲| 757午夜福利合集在线观看| 欧美日韩乱码在线| 欧美绝顶高潮抽搐喷水| 黄频高清免费视频| 不卡av一区二区三区| 男女之事视频高清在线观看| 老司机深夜福利视频在线观看| 黑丝袜美女国产一区| 99国产精品99久久久久| 久久香蕉国产精品| 精品久久久久久久人妻蜜臀av| 亚洲av成人av| 99精品在免费线老司机午夜| 一二三四社区在线视频社区8| 我的亚洲天堂| 久久精品影院6| 免费观看精品视频网站| 精品国产美女av久久久久小说| 午夜成年电影在线免费观看| 美女免费视频网站| 久久久国产成人精品二区| 男女午夜视频在线观看| 非洲黑人性xxxx精品又粗又长| ponron亚洲| bbb黄色大片| 日韩精品青青久久久久久| 美女 人体艺术 gogo| 天天躁夜夜躁狠狠躁躁| 国产麻豆成人av免费视频| 18禁观看日本| 日韩一卡2卡3卡4卡2021年| 国产精品影院久久| 国内精品久久久久久久电影| 人人妻,人人澡人人爽秒播| 久久这里只有精品19| 韩国精品一区二区三区| 精品无人区乱码1区二区| 久久欧美精品欧美久久欧美| 老司机午夜十八禁免费视频| 国产亚洲精品第一综合不卡| 成人国语在线视频| 午夜福利欧美成人| 国产精品久久电影中文字幕| 久久精品aⅴ一区二区三区四区| 法律面前人人平等表现在哪些方面| 亚洲片人在线观看| 亚洲天堂国产精品一区在线| 国产精品二区激情视频| 中文字幕久久专区| 成人亚洲精品一区在线观看| 欧美日本视频| 少妇裸体淫交视频免费看高清 | 动漫黄色视频在线观看| 欧美人与性动交α欧美精品济南到| 悠悠久久av| 国产成人av教育| 人人妻人人澡欧美一区二区| 久久久久国产一级毛片高清牌| 巨乳人妻的诱惑在线观看| 成人三级做爰电影| 国产av一区二区精品久久| 久久久久精品国产欧美久久久| 欧美午夜高清在线| 亚洲人成网站高清观看| 久久久水蜜桃国产精品网| 国产激情偷乱视频一区二区| 俄罗斯特黄特色一大片| 日韩国内少妇激情av| 国产精品久久久久久人妻精品电影| 欧美精品亚洲一区二区| 青草久久国产| 成人国产一区最新在线观看| 性欧美人与动物交配| 亚洲国产精品999在线| 国产精品久久电影中文字幕| 黄片小视频在线播放| 亚洲 欧美一区二区三区| 91成年电影在线观看| 欧美乱码精品一区二区三区| 大型av网站在线播放| 国产精品一区二区精品视频观看| 成人手机av| 亚洲专区中文字幕在线| 久久久国产成人精品二区| 男女下面进入的视频免费午夜 | 美女午夜性视频免费| 丰满人妻熟妇乱又伦精品不卡| 亚洲男人天堂网一区| 免费一级毛片在线播放高清视频| 久久人妻福利社区极品人妻图片| 不卡av一区二区三区| 国产精品久久视频播放| 国产午夜福利久久久久久| bbb黄色大片| 亚洲欧美日韩无卡精品| 高清在线国产一区| 变态另类丝袜制服| 国产熟女xx| 日日干狠狠操夜夜爽| 侵犯人妻中文字幕一二三四区| 精品不卡国产一区二区三区| 亚洲五月天丁香| 熟女电影av网| 人人妻人人澡人人看| 国产亚洲精品第一综合不卡| 国内久久婷婷六月综合欲色啪| 国产精品一区二区三区四区久久 | 午夜福利欧美成人| 亚洲天堂国产精品一区在线| 色综合站精品国产| 中文字幕精品免费在线观看视频| 亚洲精品久久国产高清桃花| 大型黄色视频在线免费观看| 成熟少妇高潮喷水视频| 亚洲成av片中文字幕在线观看| 免费在线观看日本一区| 国产精品自产拍在线观看55亚洲| 国产精品国产高清国产av| 18禁裸乳无遮挡免费网站照片 | 亚洲七黄色美女视频| 黑人巨大精品欧美一区二区mp4| 国产精华一区二区三区| av天堂在线播放| 亚洲av电影不卡..在线观看| 日韩欧美三级三区| 在线观看舔阴道视频| 日日干狠狠操夜夜爽| 侵犯人妻中文字幕一二三四区| 99精品久久久久人妻精品| 亚洲熟妇中文字幕五十中出| 中文字幕高清在线视频| 免费搜索国产男女视频| 日韩欧美 国产精品| 可以在线观看的亚洲视频| 亚洲自拍偷在线| 精品久久久久久久久久久久久 | 日韩欧美在线二视频| 国产欧美日韩精品亚洲av| 国产伦一二天堂av在线观看| 99riav亚洲国产免费| 日韩精品青青久久久久久| 欧美一级毛片孕妇| 亚洲国产欧洲综合997久久, | 国产精品自产拍在线观看55亚洲| 欧美激情久久久久久爽电影| 午夜免费激情av| 久久精品人妻少妇| 欧洲精品卡2卡3卡4卡5卡区| 久久香蕉国产精品| 激情在线观看视频在线高清| 久久九九热精品免费| 国产男靠女视频免费网站| 国内久久婷婷六月综合欲色啪| 欧美激情极品国产一区二区三区| 一级毛片女人18水好多| 亚洲成人国产一区在线观看| 亚洲无线在线观看| 很黄的视频免费| 日本免费一区二区三区高清不卡| 久久九九热精品免费| 熟女电影av网| 这个男人来自地球电影免费观看| 婷婷精品国产亚洲av| 亚洲av熟女| 搞女人的毛片| 69av精品久久久久久| 变态另类成人亚洲欧美熟女| 久久亚洲真实| 国产精品免费视频内射| 在线观看免费日韩欧美大片| 国产精品亚洲一级av第二区| 亚洲国产精品成人综合色| 黑人操中国人逼视频| 精品乱码久久久久久99久播| 国产99久久九九免费精品| 欧美日韩乱码在线| 日本黄色视频三级网站网址| 最近在线观看免费完整版| 欧美av亚洲av综合av国产av| 99热6这里只有精品| 99国产极品粉嫩在线观看| a级毛片在线看网站| 亚洲一区二区三区不卡视频| 国产精品免费一区二区三区在线| 在线看三级毛片| 亚洲美女黄片视频| 日韩大码丰满熟妇| 99国产精品一区二区三区| 久久中文字幕一级| 韩国av一区二区三区四区| 亚洲七黄色美女视频| 国产乱人伦免费视频| 老鸭窝网址在线观看| 国产视频一区二区在线看| 国产精品久久视频播放| 又黄又粗又硬又大视频| 变态另类丝袜制服| 一区二区三区国产精品乱码| 亚洲精华国产精华精| 嫩草影视91久久| 正在播放国产对白刺激| 丝袜美腿诱惑在线| 亚洲精品中文字幕在线视频| 成人亚洲精品av一区二区| 女生性感内裤真人,穿戴方法视频| 欧美日韩精品网址| 黄色片一级片一级黄色片| 在线永久观看黄色视频| 女人爽到高潮嗷嗷叫在线视频| 91国产中文字幕| 97人妻精品一区二区三区麻豆 | 日韩中文字幕欧美一区二区| 国产单亲对白刺激| 久久久水蜜桃国产精品网| 老司机在亚洲福利影院| 在线观看免费午夜福利视频| 国产免费av片在线观看野外av| 亚洲国产日韩欧美精品在线观看 | 99国产极品粉嫩在线观看| 国产aⅴ精品一区二区三区波| 免费女性裸体啪啪无遮挡网站| 丝袜在线中文字幕| 精品欧美一区二区三区在线| 久久久水蜜桃国产精品网| 国产精品99久久99久久久不卡| 日韩一卡2卡3卡4卡2021年| 日韩欧美三级三区| 欧美日韩瑟瑟在线播放| 无人区码免费观看不卡| 免费av毛片视频| 成人永久免费在线观看视频| 欧美成狂野欧美在线观看| 在线观看一区二区三区| 亚洲国产欧美网| 午夜免费激情av| 久久国产精品男人的天堂亚洲| 欧美最黄视频在线播放免费| 午夜福利在线观看吧| 最好的美女福利视频网| 人人澡人人妻人| 色播亚洲综合网| 天天躁夜夜躁狠狠躁躁| 制服人妻中文乱码| 国产一区二区在线av高清观看| 国产精品98久久久久久宅男小说| 午夜激情福利司机影院| 亚洲男人天堂网一区| 亚洲av五月六月丁香网| 日本撒尿小便嘘嘘汇集6| 中文字幕人妻丝袜一区二区| 午夜激情av网站| 久久久久久亚洲精品国产蜜桃av| 在线天堂中文资源库| 精品少妇一区二区三区视频日本电影| 夜夜夜夜夜久久久久| 日韩三级视频一区二区三区| 又紧又爽又黄一区二区| 国产v大片淫在线免费观看| 亚洲无线在线观看| 国产一级毛片七仙女欲春2 |