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

    用于機(jī)械系統(tǒng)固有頻率及阻尼比計(jì)算的改進(jìn)頻域方法

    2015-12-27 05:49:15孟慶虎朱永生刁瑞朋
    關(guān)鍵詞:阻尼比頻域時(shí)域

    孟慶虎,,朱永生,刁瑞朋

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

    ?

    用于機(jī)械系統(tǒng)固有頻率及阻尼比計(jì)算的改進(jìn)頻域方法

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

    針對(duì)傳統(tǒng)的機(jī)械系統(tǒng)固有頻率及阻尼比計(jì)算方法中,時(shí)域方法易造成主觀誤差,頻域方法存在一定精度限制的問題,提出了一種利用指數(shù)衰減正弦信號(hào)頻譜校正算法計(jì)算系統(tǒng)阻尼比和固有頻率的新方法。該方法首先利用基于漢寧窗的插值頻譜校正方法獲得頻率偏移量和衰減因子,再結(jié)合有阻尼自由振動(dòng)響應(yīng)與指數(shù)衰減正弦信號(hào)的相似性關(guān)系求取固有頻率,進(jìn)而獲得系統(tǒng)阻尼比和有阻尼固有頻率,最終計(jì)算出無阻尼固有頻率。分別使用傳統(tǒng)時(shí)域方法和本文提出的頻域方法對(duì)壓電懸臂梁的固有頻率及阻尼比進(jìn)行計(jì)算,結(jié)果表明新方法擁有更高的計(jì)算精度和更好的信噪比適應(yīng)性。與傳統(tǒng)的頻率擬合方法相比,本文所提方法也更簡單,計(jì)算效率更高。

    衰減信號(hào);頻譜校正;固有頻率;阻尼比;壓電懸臂梁

    系統(tǒng)自由振動(dòng)幅值逐步減少的過程稱為阻尼現(xiàn)象[1],其本質(zhì)是振動(dòng)系統(tǒng)能量的轉(zhuǎn)換,即從機(jī)械能轉(zhuǎn)換為熱能或其他形式能量的過程。阻尼比是指阻尼系數(shù)與臨界阻尼系數(shù)之比,用來表達(dá)結(jié)構(gòu)阻尼的大小,描述結(jié)構(gòu)在振動(dòng)過程中能量的耗散情況。較為精確地獲得系統(tǒng)的固有頻率和阻尼比,在結(jié)構(gòu)動(dòng)力學(xué)計(jì)算、結(jié)構(gòu)故障診斷、振動(dòng)實(shí)時(shí)監(jiān)控、響應(yīng)預(yù)測(cè)、荷載識(shí)別、能量轉(zhuǎn)換等課題研究和工程實(shí)際應(yīng)用方面,有著很重要的作用。

    目前常見的計(jì)算系統(tǒng)固有頻率和阻尼比的方法可分為時(shí)域法和頻域法兩種。時(shí)域法有對(duì)數(shù)衰減率法、易卜拉欣時(shí)域法(Ibrahim time domain, ITD)[2]、少時(shí)時(shí)域法(spare time domain, STD)[3]、隨機(jī)減量法[4]等;頻域法有半功率帶寬法[5]、峰值法、導(dǎo)納圓法等[6]。其中,時(shí)域方法具有操作簡單、易于實(shí)現(xiàn)等特點(diǎn),但是該方法需要人工進(jìn)行選點(diǎn),會(huì)帶來一定的主觀輸入誤差。此外,在信噪比較低的信號(hào)時(shí)域波形上也難以準(zhǔn)確地估計(jì)信號(hào)周期,這也為正確地選點(diǎn)造成了困難。半功率點(diǎn)法應(yīng)用較廣,但存在柵欄效應(yīng)、能量泄漏的離散譜線上很難精確找到半功率點(diǎn)。傳統(tǒng)的頻域識(shí)別方法在進(jìn)行快速傅里葉變換(fast Fourier transform, FFT)的時(shí)候,亦會(huì)導(dǎo)致一定的能量泄漏,從而在精度上受到一定的限制。如何提高阻尼比識(shí)別精度,使之便于工程實(shí)際應(yīng)用,便成為了一項(xiàng)難題。

    受到指數(shù)衰減信號(hào)頻域插值校正理論的啟發(fā),本文提出了一種新型的識(shí)別系統(tǒng)固有頻率和阻尼比的指數(shù)衰減正弦信號(hào)頻譜校正方法。利用離散傅里葉變換后兩相鄰譜線的比值關(guān)系可以精確獲得衰減因子,而這個(gè)重要的頻譜校正參數(shù)與阻尼比存在確定的關(guān)系。與傳統(tǒng)的時(shí)域方法和頻域方法相比,本文方法利用校正的信號(hào)參數(shù)進(jìn)行計(jì)算避免了人工選點(diǎn)造成的主觀誤差,而通過插值方法校正出的固有頻率和阻尼比基本不受采樣率的影響,具有較高的估計(jì)精度。同時(shí),本文方法又具有較強(qiáng)的適應(yīng)性和更簡單的計(jì)算過程,并適用于多模態(tài)的阻尼比的計(jì)算。使用本文方法計(jì)算了壓電懸臂梁固有頻率和阻尼比,并與傳統(tǒng)時(shí)域方法進(jìn)行了對(duì)比,獲得了較為準(zhǔn)確的結(jié)果,證明本文所提方法可以作為計(jì)算系統(tǒng)固有頻率和阻尼比的一種可選方法。

    1 方法描述

    1.1 基于頻域差值的非平穩(wěn)信號(hào)參數(shù)識(shí)別方法

    由M個(gè)頻率分量組成的指數(shù)衰減正弦信號(hào)一般形式為

    (1)

    式中:Am、fm、φm為第m個(gè)頻率分量的幅值、頻率、初相位;αm為第m個(gè)頻率分量的衰減因子,該參數(shù)表征指數(shù)衰減正弦信號(hào)的衰減速率,衰減因子越大,信號(hào)幅值的衰減速率越快;N為采樣點(diǎn)數(shù)。其中一單頻指數(shù)衰減正弦信號(hào)可表示為

    (2)

    式(2)的離散傅里葉變換為

    X[k]=A0ejφ0W[k-f0;α0],k=0,…,N-1

    (3)

    式中:W[k-fm;αm]由下式給出

    (4)

    假設(shè)頻率f0的理論位置為θ0,其兩側(cè)的兩根離散譜線的位置分別為k0-1和k0,則存在如下關(guān)系

    (5)

    由式(3)可知,f0兩側(cè)的兩根譜線對(duì)應(yīng)值為

    (6)

    由式(5)和式(6)可知,f0兩側(cè)兩根譜線的離散傅里葉變換比值為

    (7)

    定義變量

    (8)

    將式(4)和式(8)代入式(7),可解出τ的值為

    (9)

    由式(8)可知,τ為一復(fù)數(shù),根據(jù)復(fù)數(shù)模和幅角的定義可知

    (10)

    式中:arg(τ)為τ的幅角。

    因此,指數(shù)衰減正弦信號(hào)頻率的校正值為

    (11)

    由式(3)可知

    A0=|X[k]|/|W[-δ0;α0]|

    φ0=arg(X[k])-arg(W[-δ0;α0])+π/2

    (12)

    基于矩形窗的頻域差值方法存在頻譜泄漏大、衰減因子偏差較大的不足。為了有效地抑制頻譜泄漏的影響,提高頻譜參數(shù)估計(jì)的精度,刁瑞朋等提出了基于漢寧窗的指數(shù)衰減正弦信號(hào)矢量插值頻譜校正方法并推導(dǎo)出了具體的實(shí)現(xiàn)公式,可獲得衰減因子和頻率偏移量的校正值[7]

    (13)

    (14)

    頻率的校正值為

    (15)

    幅值和相位的校正值為

    (16)

    (17)

    至此,我們獲得了指數(shù)衰減正弦信號(hào)基于漢寧窗的插值頻譜校正公式。為了方便起見,將此插值頻譜校正方法簡稱為HDIDFT(Hanning damping interpolated discrete Fourier transform)。

    本文擬采用基于漢寧窗的指數(shù)衰減正弦信號(hào)矢量插值頻譜校正方法,并根據(jù)振動(dòng)力學(xué)的公式,提出一種計(jì)算系統(tǒng)阻尼比和固有頻率的方法。

    1.2 基于頻域方法的系統(tǒng)固有頻率和阻尼比計(jì)算

    由振動(dòng)力學(xué)知識(shí)可知,系統(tǒng)欠阻尼情形單自由度系統(tǒng)的自由振動(dòng)響應(yīng)[8]為

    (18)

    根據(jù)三角函數(shù)和差化積公式,式(18)可表示為

    (19)

    其中

    式(19)表明,欠阻尼情形的單自由度系統(tǒng)對(duì)初始條件的響應(yīng)是一種振幅按指數(shù)規(guī)律逐漸衰減的簡諧振動(dòng),與指數(shù)衰減信號(hào)式(1)具有相同的表示形式。對(duì)比兩個(gè)公式可以看出,式(19)中的-ζωn即是式(1)中的衰減因子項(xiàng)2πα,而有阻尼系統(tǒng)的固有頻率ωd與式(1)中2πf意義相同,即2πf為系統(tǒng)的有阻尼固有頻率,而衰減因子2πα和頻率f剛好是本文提出的頻譜校正方法所校正出的兩個(gè)參數(shù),見式(13)和式(15)。由此可見,指數(shù)衰減信號(hào)的頻域插值頻譜校正方法可以用來計(jì)算欠阻尼情形下的阻尼比及固有頻率。

    下面給出本文提出的利用指數(shù)衰減信號(hào)頻譜校正方法計(jì)算阻尼比的公式。由式(1)和式(19)可知

    α=-ζωn/2π

    (20)

    則有阻尼固有頻率ωd和衰減因子α的比值為

    (21)

    根據(jù)式(20)即可求出阻尼比

    (22)

    式中:fd=ωd/2π為有阻尼固有頻率,單位為Hz。

    具體計(jì)算阻尼比時(shí),首先用指數(shù)衰減正弦信號(hào)頻域插值校正方法獲得衰減因子α和有阻尼固有頻率fd兩個(gè)值,代入式(21)即可求出阻尼比ζ和固有頻率ωd。然后根據(jù)ωd與ωn的關(guān)系即可求出系統(tǒng)的無阻尼固有頻率為

    (23)

    2 壓電懸臂梁固有頻率及阻尼比實(shí)驗(yàn)

    隨著無線傳感設(shè)備及微機(jī)電系統(tǒng)的日益發(fā)展,以化學(xué)電池為主的供能方式的弊端日漸顯露。為實(shí)現(xiàn)微機(jī)電系統(tǒng)的長時(shí)間供能,需要大力發(fā)展能量回收技術(shù),將外界能量轉(zhuǎn)換成系統(tǒng)可利用的電能。由于振動(dòng)現(xiàn)象幾乎無處不在且具有較高的能量密度,因此對(duì)振動(dòng)能量的轉(zhuǎn)化及回收利用的研究在近年逐步興起。目前,振動(dòng)能量回收技術(shù)主要分為靜電式、電磁式和壓電式3類。其中,壓電式方法以其具有的結(jié)構(gòu)簡單、不發(fā)熱、無電磁干擾、清潔環(huán)保和易于微型化等諸多優(yōu)點(diǎn)而得到了極大重視[9-10]。

    壓電發(fā)電裝置利用環(huán)境中的振動(dòng)使壓電材料發(fā)生變形來產(chǎn)生電能,其核心部件是換能器,它決定了系統(tǒng)的能量轉(zhuǎn)換效率[10]。目前,最常用的換能器結(jié)構(gòu)為懸臂梁式結(jié)構(gòu),如圖1所示。換能器設(shè)計(jì)的一個(gè)基本要求是其諧振頻率與外界的振動(dòng)頻率一致,使電能輸出最大,因此需要對(duì)其形狀和結(jié)構(gòu)進(jìn)行優(yōu)化設(shè)計(jì)。理論分析壓電懸臂梁的電能輸出時(shí),將其阻尼設(shè)為瑞利阻尼,忽略質(zhì)量阻尼,而只考慮材料的結(jié)構(gòu)阻尼,此時(shí)壓電懸臂梁的機(jī)電耦合微分方程中的阻尼系數(shù)可由阻尼比和固有頻率計(jì)算得到。

    圖1 壓電懸臂梁結(jié)構(gòu)圖

    具體實(shí)驗(yàn)中,將壓電懸臂梁固定在Z型支架上,如圖2所示。給其自由端一個(gè)微小的瞬時(shí)位移激勵(lì),然后讓其自由振動(dòng),利用電渦流位移傳感器測(cè)量質(zhì)量塊在垂直方向發(fā)生的位移量,進(jìn)而利用該位移量計(jì)算求取其固有頻率及結(jié)構(gòu)阻尼比。采集信號(hào)時(shí),設(shè)置采樣率為1 000 Hz,從采集信號(hào)中截取采樣時(shí)長為601 ms的數(shù)據(jù)。圖3給出了測(cè)得的質(zhì)量塊自由衰減時(shí)域曲線和FFT波形圖。

    圖2 壓電懸臂梁實(shí)驗(yàn)圖

    3 對(duì)比分析

    下面分別用時(shí)域方法和基于頻譜校正的頻域方法計(jì)算壓電懸臂梁的固有頻率和阻尼比,并將兩種方法進(jìn)行對(duì)比。

    3.1 時(shí)域方法

    時(shí)域方法[11-13]是一種簡便有效的計(jì)算自由振動(dòng)阻尼比和固有頻率的方法。具體來說,分別記錄時(shí)域曲線上相隔固定周期的幾個(gè)時(shí)刻的位移值,然后利用式(24)計(jì)算出系統(tǒng)的阻尼比,根據(jù)采集到的信號(hào)周期估計(jì)出固有頻率[10]

    (24)

    式中:N為計(jì)算阻尼比需要的數(shù)據(jù)點(diǎn)數(shù);τ為相鄰兩個(gè)數(shù)據(jù)點(diǎn)間所包含的信號(hào)周期數(shù)。

    (a)時(shí)域波形圖

    (b)FFT頻譜圖圖3 壓電換能器自由衰減時(shí)域波形及FFT頻譜圖

    本例中每隔2個(gè)周期選擇一個(gè)數(shù)據(jù)點(diǎn),共選擇7個(gè)數(shù)據(jù)點(diǎn),即N=7,τ=2。數(shù)據(jù)點(diǎn)的位置和對(duì)應(yīng)的位移值分別見圖4和表1。

    圖4 計(jì)算阻尼比的7個(gè)數(shù)據(jù)點(diǎn)

    使用時(shí)域方法,利用7個(gè)數(shù)據(jù)點(diǎn)計(jì)算出了6組阻尼比和固有頻率。表2給出了使用不同區(qū)間的數(shù)據(jù)的計(jì)算結(jié)果。

    表1 不同數(shù)據(jù)點(diǎn)的位移

    表2 使用不同區(qū)間計(jì)算出的阻尼比

    3.2 基于衰減信號(hào)參數(shù)估計(jì)的頻域方法

    為了用頻域方法計(jì)算懸臂梁的固有頻率和阻尼比,首先需要用頻域插值校正出固有頻率和衰減因子的準(zhǔn)確值。表3給出了FFT和HDIDFT計(jì)算出的信號(hào)參數(shù),由于FFT不能計(jì)算出衰減因子,因此該項(xiàng)的值為空。

    將表3中由HDIDFT計(jì)算出的有阻尼固有頻率和衰減因子代入式(22),則阻尼比為

    表3 信號(hào)的FFT和HDIDFT計(jì)算結(jié)果

    (25)

    將有阻尼固有頻率的值代入式(23),則系統(tǒng)的無阻尼固有頻率為

    fn=[1-0.01982]-1/2×28.5=28.506 (Hz)

    (26)

    以上對(duì)阻尼比和固有頻率的計(jì)算結(jié)果表明,時(shí)域方法和本文提出的利用指數(shù)衰減正弦信號(hào)頻譜校正算法計(jì)算系統(tǒng)固有頻率和阻尼比的方法具有基本相似的計(jì)算結(jié)果。本文方法還可應(yīng)用于多模態(tài)的阻尼比的計(jì)算中。由于本文方法校正出的固有頻率值和衰減因子相對(duì)準(zhǔn)確,基于二者計(jì)算出的參數(shù)也較時(shí)域方法要準(zhǔn)確些。

    4 結(jié) 論

    指數(shù)衰減正弦信號(hào)廣泛存在于機(jī)械系統(tǒng)中,因此對(duì)該類信號(hào)進(jìn)行分析和參數(shù)識(shí)別具有重要的理論和現(xiàn)實(shí)意義。在分析有阻尼自由振動(dòng)響應(yīng)方程的基礎(chǔ)上,結(jié)合指數(shù)衰減正弦信號(hào)插值頻譜校正技術(shù),本文提出一種利用指數(shù)衰減正弦信號(hào)頻譜校正算法來計(jì)算阻尼比和固有頻率的新方法,并給出了計(jì)算公式。該方法利用校正的信號(hào)參數(shù)進(jìn)行計(jì)算,從而避免了人工選點(diǎn)造成的主觀誤差,而通過插值方法校正出的固有頻率和衰減因子基本不受采樣率的影響,具有較高的估計(jì)精度。為了驗(yàn)證該方法的性能,使用其對(duì)壓電懸臂梁的固有頻率和阻尼比進(jìn)行計(jì)算。與傳統(tǒng)時(shí)域方法相比,本文所提方法適應(yīng)性更好,精度更高,計(jì)算過程也更簡單。與傳統(tǒng)的頻率擬合方法相比,本文所提方法也更簡單,計(jì)算效率更高。

    [1] CHOPRA A K. Dynamic of structures [M]. 2nd ed. Englewood Cliffs, New Jersey, USA: Prentice Hall, 2000: 12-13.

    [2] IBRAHIM S R, MIKULCIK E C. Method for the direct identification of vibration parameters from the free response [J]. Shocked Vibration Bulletin, 1977, 47: 183-198.

    [3] IBRAHIM S R. An approach for reducing computational requirements in modal identification [J]. AIAA Journal, 1986, 24(10): 1725-1727.

    [4] IBRAHIM S R. Random decrement technique for modal identification of structures [J]. Journal of Spacecraft and Rockets, 1997, 14(11): 696-700.

    [5] 黃文虎, 夏松波, 劉巖. 設(shè)備故障診斷原理、技術(shù)及應(yīng)用 [M]. 北京: 科學(xué)出版社, 1997: 53-62.

    [6] 周傳榮, 趙淳生. 機(jī)械振動(dòng)參數(shù)識(shí)別及其應(yīng)用 [M]. 北京: 科學(xué)出版社, 1989: 173-176.

    [7] 刁瑞朋, 孟慶豐. 一種衰減信號(hào)加窗頻域插值算法 [J]. 西安交通大學(xué)學(xué)報(bào), 2013, 47(7): 85-90. DIAO Ruipeng, MENG Qingfeng. Interpolation algorithm for discrete Fourier transforms of weighted damped sinusoidal signals [J]. Journal of Xi’an Jiaotong University, 2013, 47(7): 85-90.

    [8] 倪振華. 振動(dòng)力學(xué) [M]. 西安: 西安交通大學(xué)出版社, 1989: 61-69.

    [9] 劉祥建, 陳仁文. 壓電振動(dòng)能量收集裝置研究現(xiàn)狀及發(fā)展趨勢(shì) [J]. 振動(dòng)與沖擊, 2012, 31(16): 169-176. LIU Xiangjian, CHEN Renwen. Current situation and developing trend of piezoelectric vibration energy harvesters [J]. Journal of Vibration and Shock, 2012, 31(16): 169-176.

    [10]王宏金. 基于壓電材料的振動(dòng)能量收集技術(shù)及其應(yīng)用研究 [D]. 西安: 西安交通大學(xué)機(jī)械工程學(xué)院, 2013.

    [11]祁泉泉. 基于振動(dòng)信號(hào)的結(jié)構(gòu)參數(shù)識(shí)別系統(tǒng)方法研究 [D]. 北京: 清華大學(xué)土木工程系, 2011.

    [12]段志平, 張亞. 結(jié)構(gòu)阻尼識(shí)別的方法及比較 [J]. 福州大學(xué)學(xué)報(bào): 自然科學(xué)版, 2005, 33(S1): 208-212. DUAN Zhiping, ZHANG Ya. Structural damping identification methodology and comparison [J]. Journal of Fuzhou University: Natural Science. 2005, 33(S1): 208-212.

    [13]李中付, 華宏星, 宋漢文, 等. 用時(shí)域峰值法計(jì)算頻率和阻尼 [J]. 振動(dòng)與沖擊, 2001, 20(3): 7-8. LI Zhongfu, HUA Hongxing, SONG Hanwen, et al. Frequency and damping calculation based on the peak time domain method [J]. Journal of Vibration and Shock, 2001, 20(3): 7-8.

    (編輯 武紅江)

    A New Frequency-Domain Algorithm to Calculate Natural Frequency and Damping Ratio of Mechanical Systems

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

    A new algorithm to calculate the system damping ratio and natural frequency using frequency spectrum correction of exponentially decaying sinusoidal signals is proposed to solve the problem that time domain methods are easy to cause subjective errors and there exists the accuracy limit problem in frequency domain methods in traditional methods to calculate natural frequencies and damping ratios of mechanical systems. The interpolation spectrum correcting method based on Hanning-window is used to obtain frequency deviations and damped factors. Then the natural frequency is calculated by combining similarity relations between damping free vibration responses and exponentially decaying sinusoidal signals, and the system damping ratio and damped natural frequency are obtained. The undamped natural frequency is finally calculated. Both a traditional time domain method and the proposed frequency domain method are respectively used to calculate the natural frequency and damping ratio of a piezoelectric cantilever beam. The results show that the proposed algorithm has a higher calculation accuracy and better adaption of signal to noise ratio. Comparisons with the traditional frequency fitting method show that the new algorithm is simpler and its computational efficiency is higher.

    damped signal; frequency spectrum correction; natural frequency; damping ratio; piezoelectric cantilever

    2015-03-16 。 作者簡介:孟慶虎(1983—),男,博士生;朱永生(通信作者),男,教授,博士生導(dǎo)師。 基金項(xiàng)目:國家自然科學(xué)基金資助項(xiàng)目(51275380)。

    時(shí)間:2015-05-21

    10.7652/xjtuxb201508001

    TN911.72

    A

    0253-987X(2015)08-0001-05

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

    猜你喜歡
    阻尼比頻域時(shí)域
    基于細(xì)觀結(jié)構(gòu)的原狀黃土動(dòng)彈性模量和阻尼比試驗(yàn)研究
    地震研究(2021年1期)2021-04-13 01:05:24
    黏滯阻尼器在時(shí)程分析下的附加有效阻尼比研究
    基于時(shí)域信號(hào)的三電平逆變器復(fù)合故障診斷
    頻域稀疏毫米波人體安檢成像處理和快速成像稀疏陣列設(shè)計(jì)
    波形分析法求解公路橋梁阻尼比的探討
    上海公路(2018年3期)2018-03-21 05:55:48
    結(jié)構(gòu)構(gòu)件阻尼比對(duì)大跨度懸索橋地震響應(yīng)的影響
    基于極大似然準(zhǔn)則與滾動(dòng)時(shí)域估計(jì)的自適應(yīng)UKF算法
    基于改進(jìn)Radon-Wigner變換的目標(biāo)和拖曳式誘餌頻域分離
    基于時(shí)域逆濾波的寬帶脈沖聲生成技術(shù)
    一種基于頻域的QPSK窄帶干擾抑制算法
    久久伊人香网站| 久久热在线av| 久久天堂一区二区三区四区| 亚洲国产欧洲综合997久久, | 首页视频小说图片口味搜索| 性欧美人与动物交配| 国产成人啪精品午夜网站| 可以在线观看毛片的网站| 黑人操中国人逼视频| 亚洲国产精品999在线| 97人妻精品一区二区三区麻豆 | 天堂动漫精品| 国产在线观看jvid| av片东京热男人的天堂| 操出白浆在线播放| 美女 人体艺术 gogo| 夜夜夜夜夜久久久久| 国产激情久久老熟女| 成年女人毛片免费观看观看9| 国产日本99.免费观看| 久久久精品国产亚洲av高清涩受| 亚洲欧美激情综合另类| 在线av久久热| 久久伊人香网站| 精品久久久久久久久久免费视频| 在线国产一区二区在线| 好看av亚洲va欧美ⅴa在| 免费无遮挡裸体视频| 亚洲三区欧美一区| 欧美一区二区精品小视频在线| 国产男靠女视频免费网站| 一进一出好大好爽视频| 免费高清在线观看日韩| 色播亚洲综合网| 在线观看免费午夜福利视频| 91老司机精品| 女性被躁到高潮视频| 色婷婷久久久亚洲欧美| 俺也久久电影网| 久久久国产成人免费| 免费看日本二区| 成人一区二区视频在线观看| 不卡av一区二区三区| 国产99白浆流出| 巨乳人妻的诱惑在线观看| 他把我摸到了高潮在线观看| 黄色视频,在线免费观看| 亚洲国产欧美一区二区综合| 午夜激情福利司机影院| 久久久久精品国产欧美久久久| 在线观看舔阴道视频| 黄色毛片三级朝国网站| 老司机在亚洲福利影院| 悠悠久久av| 亚洲,欧美精品.| 免费看美女性在线毛片视频| 真人一进一出gif抽搐免费| 欧美一级毛片孕妇| 国产精品野战在线观看| 亚洲电影在线观看av| 999久久久精品免费观看国产| 国产午夜福利久久久久久| 免费无遮挡裸体视频| 国产精品野战在线观看| 欧美大码av| www.精华液| 国产97色在线日韩免费| 99久久久亚洲精品蜜臀av| а√天堂www在线а√下载| 欧美性猛交╳xxx乱大交人| av有码第一页| 丝袜美腿诱惑在线| 日韩视频一区二区在线观看| 男人操女人黄网站| 色综合欧美亚洲国产小说| 校园春色视频在线观看| 法律面前人人平等表现在哪些方面| 人妻久久中文字幕网| 91老司机精品| 国产视频一区二区在线看| 悠悠久久av| 日韩成人在线观看一区二区三区| 午夜福利18| 国产免费av片在线观看野外av| 国产精品二区激情视频| 精品久久久久久久毛片微露脸| 亚洲精品久久国产高清桃花| 国产单亲对白刺激| 正在播放国产对白刺激| 亚洲国产欧美网| 中文亚洲av片在线观看爽| 波多野结衣av一区二区av| 欧美在线一区亚洲| 成人午夜高清在线视频 | 亚洲 欧美一区二区三区| 又黄又爽又免费观看的视频| 每晚都被弄得嗷嗷叫到高潮| 亚洲片人在线观看| 国产精品日韩av在线免费观看| 一个人观看的视频www高清免费观看 | 国产黄色小视频在线观看| 高清毛片免费观看视频网站| 少妇的丰满在线观看| 国产激情偷乱视频一区二区| av免费在线观看网站| 欧美+亚洲+日韩+国产| 日韩欧美三级三区| 亚洲av第一区精品v没综合| 国产极品粉嫩免费观看在线| 国产aⅴ精品一区二区三区波| 老司机深夜福利视频在线观看| 18禁观看日本| 88av欧美| 欧美成人一区二区免费高清观看 | 成人亚洲精品一区在线观看| 成人手机av| 免费看日本二区| 99久久99久久久精品蜜桃| 亚洲精品国产区一区二| 欧美黑人精品巨大| 国内毛片毛片毛片毛片毛片| 精品熟女少妇八av免费久了| 国产av一区在线观看免费| av电影中文网址| 国语自产精品视频在线第100页| 特大巨黑吊av在线直播 | 此物有八面人人有两片| 亚洲国产欧洲综合997久久, | 亚洲国产高清在线一区二区三 | 亚洲五月婷婷丁香| 成人亚洲精品av一区二区| 精品国产国语对白av| 老熟妇仑乱视频hdxx| 国产精品久久久人人做人人爽| 国产伦人伦偷精品视频| 无人区码免费观看不卡| 成人18禁高潮啪啪吃奶动态图| 亚洲av片天天在线观看| 亚洲一区二区三区色噜噜| 日本 av在线| 久久久久久久精品吃奶| 国内精品久久久久久久电影| 久久久久久久久免费视频了| 激情在线观看视频在线高清| 此物有八面人人有两片| 亚洲男人天堂网一区| 久久精品国产亚洲av高清一级| 美女 人体艺术 gogo| 最近最新中文字幕大全电影3 | 亚洲国产精品sss在线观看| 久久久久久久午夜电影| 日本 欧美在线| 首页视频小说图片口味搜索| 啦啦啦观看免费观看视频高清| 色婷婷久久久亚洲欧美| 日韩免费av在线播放| 香蕉av资源在线| 亚洲无线在线观看| 黄片大片在线免费观看| 成人国语在线视频| 久久久久久久久免费视频了| 国产精品野战在线观看| 日本黄色视频三级网站网址| 中文字幕人成人乱码亚洲影| 国产精品自产拍在线观看55亚洲| 好男人在线观看高清免费视频 | 国语自产精品视频在线第100页| 黄色 视频免费看| 一进一出抽搐gif免费好疼| 一a级毛片在线观看| 欧美性猛交╳xxx乱大交人| 国产成人一区二区三区免费视频网站| 欧美最黄视频在线播放免费| 精品国产超薄肉色丝袜足j| 国产一卡二卡三卡精品| 99国产精品一区二区蜜桃av| 女人被狂操c到高潮| 成人精品一区二区免费| 亚洲成人久久性| 亚洲精品粉嫩美女一区| 欧美成狂野欧美在线观看| 中文字幕最新亚洲高清| 国产精品98久久久久久宅男小说| 久热爱精品视频在线9| 波多野结衣av一区二区av| 国产一区二区三区在线臀色熟女| 禁无遮挡网站| 久久久久精品国产欧美久久久| 欧美国产精品va在线观看不卡| 99久久久亚洲精品蜜臀av| 欧美激情 高清一区二区三区| av超薄肉色丝袜交足视频| 国产久久久一区二区三区| 婷婷精品国产亚洲av| 久久99热这里只有精品18| 久久精品国产亚洲av高清一级| 久久精品国产99精品国产亚洲性色| 美女国产高潮福利片在线看| 宅男免费午夜| 国产真实乱freesex| 亚洲av片天天在线观看| 亚洲av成人不卡在线观看播放网| 国产91精品成人一区二区三区| 男人的好看免费观看在线视频 | 亚洲国产精品成人综合色| 欧美成人一区二区免费高清观看 | 少妇被粗大的猛进出69影院| 18禁裸乳无遮挡免费网站照片 | 天天躁夜夜躁狠狠躁躁| 国产三级黄色录像| 日本一区二区免费在线视频| 黑人操中国人逼视频| 最好的美女福利视频网| 少妇裸体淫交视频免费看高清 | 久久久久国内视频| 亚洲av电影不卡..在线观看| 99国产精品一区二区三区| 日韩精品青青久久久久久| 啦啦啦免费观看视频1| 91成年电影在线观看| 久9热在线精品视频| 久久久久久久久免费视频了| 国产乱人伦免费视频| 中文字幕人成人乱码亚洲影| 国产成人精品无人区| 亚洲av电影不卡..在线观看| 在线十欧美十亚洲十日本专区| 国产成人系列免费观看| 亚洲黑人精品在线| 成年女人毛片免费观看观看9| 岛国在线观看网站| 在线观看一区二区三区| 免费搜索国产男女视频| 日本三级黄在线观看| 国产乱人伦免费视频| 久久亚洲精品不卡| 人成视频在线观看免费观看| 亚洲免费av在线视频| netflix在线观看网站| 午夜福利18| 中文字幕另类日韩欧美亚洲嫩草| 国产视频内射| 又黄又粗又硬又大视频| 久久人人精品亚洲av| 亚洲欧美一区二区三区黑人| 国产精品亚洲av一区麻豆| www日本在线高清视频| 亚洲第一av免费看| 免费观看人在逋| 一本一本综合久久| 色综合婷婷激情| 搡老岳熟女国产| 1024视频免费在线观看| 国产精品一区二区免费欧美| 成人av一区二区三区在线看| 久久久国产精品麻豆| 成年女人毛片免费观看观看9| 婷婷亚洲欧美| 久久天堂一区二区三区四区| 男人操女人黄网站| 免费一级毛片在线播放高清视频| 欧美又色又爽又黄视频| 一个人免费在线观看的高清视频| 免费在线观看成人毛片| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美性长视频在线观看| av在线播放免费不卡| 丝袜人妻中文字幕| 国产在线精品亚洲第一网站| 久久久久久人人人人人| 露出奶头的视频| 久久久国产成人精品二区| 国产精品免费视频内射| 亚洲男人天堂网一区| 女警被强在线播放| 久久精品国产综合久久久| 亚洲国产精品999在线| 啪啪无遮挡十八禁网站| 中文在线观看免费www的网站 | 波多野结衣高清无吗| 香蕉丝袜av| 成人18禁在线播放| 一区二区三区激情视频| www.999成人在线观看| 欧美人与性动交α欧美精品济南到| 亚洲 欧美一区二区三区| 好男人电影高清在线观看| 男女之事视频高清在线观看| 麻豆av在线久日| 久久久久久久久免费视频了| 久久国产精品影院| 在线免费观看的www视频| 精品一区二区三区四区五区乱码| 欧美精品啪啪一区二区三区| 日本 av在线| 好男人在线观看高清免费视频 | 19禁男女啪啪无遮挡网站| 欧美日韩一级在线毛片| 一本综合久久免费| 午夜a级毛片| 一区二区三区精品91| 免费观看精品视频网站| 亚洲精品在线美女| 99在线人妻在线中文字幕| 免费在线观看亚洲国产| 亚洲午夜精品一区,二区,三区| 精品国产美女av久久久久小说| 真人一进一出gif抽搐免费| 黄色毛片三级朝国网站| 欧美又色又爽又黄视频| 91av网站免费观看| 在线看三级毛片| 国产精品亚洲美女久久久| 一区福利在线观看| 18禁裸乳无遮挡免费网站照片 | 好男人电影高清在线观看| 国产高清有码在线观看视频 | 狂野欧美激情性xxxx| 亚洲五月色婷婷综合| 国产精品九九99| 非洲黑人性xxxx精品又粗又长| 757午夜福利合集在线观看| 色综合亚洲欧美另类图片| 久9热在线精品视频| 久久精品亚洲精品国产色婷小说| 午夜a级毛片| 成人手机av| 欧美日本视频| 一区二区日韩欧美中文字幕| 亚洲av中文字字幕乱码综合 | 亚洲成a人片在线一区二区| 日韩大尺度精品在线看网址| 侵犯人妻中文字幕一二三四区| 欧美激情 高清一区二区三区| 欧美日韩一级在线毛片| 欧美性长视频在线观看| 国产av一区二区精品久久| 欧美日韩福利视频一区二区| 在线国产一区二区在线| 草草在线视频免费看| 精品一区二区三区四区五区乱码| 色老头精品视频在线观看| 一二三四社区在线视频社区8| 亚洲人成伊人成综合网2020| 曰老女人黄片| 18禁国产床啪视频网站| 国产成人一区二区三区免费视频网站| 亚洲国产精品合色在线| 1024香蕉在线观看| 热99re8久久精品国产| 精品福利观看| 琪琪午夜伦伦电影理论片6080| 午夜久久久在线观看| 日韩精品青青久久久久久| 国产成人影院久久av| 夜夜躁狠狠躁天天躁| 亚洲国产精品sss在线观看| 啪啪无遮挡十八禁网站| 黄色毛片三级朝国网站| 久久狼人影院| 人人妻人人澡人人看| 老熟妇乱子伦视频在线观看| 哪里可以看免费的av片| 首页视频小说图片口味搜索| 亚洲国产高清在线一区二区三 | 欧美中文日本在线观看视频| 亚洲电影在线观看av| 一夜夜www| 国产精品一区二区免费欧美| 波多野结衣av一区二区av| 亚洲成av人片免费观看| 香蕉久久夜色| 中文字幕人妻丝袜一区二区| 亚洲欧美激情综合另类| 久热爱精品视频在线9| 日韩视频一区二区在线观看| 制服人妻中文乱码| 波多野结衣av一区二区av| 1024手机看黄色片| 国产伦一二天堂av在线观看| 国产高清激情床上av| 热99re8久久精品国产| 亚洲国产高清在线一区二区三 | 精品久久久久久成人av| 嫩草影视91久久| 人成视频在线观看免费观看| 波多野结衣巨乳人妻| 欧美激情 高清一区二区三区| www国产在线视频色| 国产精品影院久久| 18禁裸乳无遮挡免费网站照片 | 国内揄拍国产精品人妻在线 | 白带黄色成豆腐渣| 亚洲精品美女久久久久99蜜臀| 特大巨黑吊av在线直播 | 久久久久久亚洲精品国产蜜桃av| 精品一区二区三区视频在线观看免费| 亚洲七黄色美女视频| 中文字幕高清在线视频| 亚洲片人在线观看| 免费在线观看亚洲国产| 妹子高潮喷水视频| av视频在线观看入口| 99久久99久久久精品蜜桃| 亚洲成av片中文字幕在线观看| 成熟少妇高潮喷水视频| 国产一级毛片七仙女欲春2 | 欧美激情久久久久久爽电影| 国产午夜精品久久久久久| 国产激情久久老熟女| 欧美日本亚洲视频在线播放| 亚洲色图 男人天堂 中文字幕| 99国产精品一区二区三区| 国产av不卡久久| www.精华液| 精品人妻1区二区| 久久性视频一级片| 搞女人的毛片| 色综合站精品国产| 日韩大尺度精品在线看网址| 草草在线视频免费看| 中出人妻视频一区二区| 国产精品久久久久久亚洲av鲁大| 亚洲色图 男人天堂 中文字幕| 亚洲欧美精品综合一区二区三区| 制服丝袜大香蕉在线| 99精品在免费线老司机午夜| 欧美午夜高清在线| 成人三级做爰电影| 亚洲精品在线美女| 男女视频在线观看网站免费 | 一卡2卡三卡四卡精品乱码亚洲| 欧美乱妇无乱码| 精品福利观看| 嫩草影院精品99| 亚洲在线自拍视频| 婷婷六月久久综合丁香| 久久久国产成人精品二区| 日韩欧美国产一区二区入口| 亚洲色图 男人天堂 中文字幕| 国产精品1区2区在线观看.| 国产亚洲精品久久久久5区| 国产高清videossex| 精品国产一区二区三区四区第35| av福利片在线| 妹子高潮喷水视频| 最近最新免费中文字幕在线| 两性午夜刺激爽爽歪歪视频在线观看 | 国产主播在线观看一区二区| 给我免费播放毛片高清在线观看| 久久久久国产精品人妻aⅴ院| 欧美激情极品国产一区二区三区| 伦理电影免费视频| 久久国产乱子伦精品免费另类| 麻豆成人午夜福利视频| 国产乱人伦免费视频| 90打野战视频偷拍视频| 精品国产乱子伦一区二区三区| 免费av毛片视频| 99久久综合精品五月天人人| 2021天堂中文幕一二区在线观 | 麻豆成人av在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产1区2区3区精品| 久久精品影院6| 婷婷丁香在线五月| 亚洲真实伦在线观看| 午夜激情福利司机影院| 亚洲中文字幕一区二区三区有码在线看 | 免费在线观看日本一区| 国产精品久久久人人做人人爽| 精品久久久久久成人av| 无遮挡黄片免费观看| 亚洲欧美激情综合另类| 亚洲成人免费电影在线观看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲人成网站在线播放欧美日韩| 免费人成视频x8x8入口观看| 欧美日本视频| 国产精品精品国产色婷婷| 黄片播放在线免费| 精品乱码久久久久久99久播| 国产激情久久老熟女| 国产不卡一卡二| 99久久精品国产亚洲精品| 97碰自拍视频| 中文字幕精品亚洲无线码一区 | 国产精品爽爽va在线观看网站 | 亚洲精品国产区一区二| www.www免费av| 黑人巨大精品欧美一区二区mp4| 国产精品久久视频播放| 曰老女人黄片| 日本a在线网址| 搞女人的毛片| 亚洲av电影在线进入| 国产激情偷乱视频一区二区| 亚洲精品国产精品久久久不卡| 黄色a级毛片大全视频| 国产精品二区激情视频| 欧美中文综合在线视频| 巨乳人妻的诱惑在线观看| 97超级碰碰碰精品色视频在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲人成电影免费在线| 成人av一区二区三区在线看| 久久午夜亚洲精品久久| 亚洲专区国产一区二区| 亚洲精华国产精华精| 老熟妇乱子伦视频在线观看| 国产成人精品久久二区二区91| 又黄又粗又硬又大视频| 香蕉丝袜av| 亚洲国产日韩欧美精品在线观看 | 午夜激情福利司机影院| 神马国产精品三级电影在线观看 | 欧美一级毛片孕妇| 午夜福利18| 首页视频小说图片口味搜索| 精品电影一区二区在线| 亚洲人成电影免费在线| 91字幕亚洲| 一进一出抽搐动态| 在线观看www视频免费| 久久性视频一级片| 欧美日韩黄片免| 禁无遮挡网站| 精品国产超薄肉色丝袜足j| 日韩高清综合在线| 精品国产亚洲在线| tocl精华| 999久久久精品免费观看国产| 香蕉国产在线看| 欧美三级亚洲精品| 后天国语完整版免费观看| 国产一级毛片七仙女欲春2 | 亚洲成人免费电影在线观看| 国产成人影院久久av| 亚洲精品在线美女| 一本一本综合久久| 老汉色∧v一级毛片| 久久性视频一级片| 亚洲一区二区三区色噜噜| 成年人黄色毛片网站| 亚洲精品国产区一区二| 91国产中文字幕| 99久久久亚洲精品蜜臀av| 国产成人一区二区三区免费视频网站| 男人舔奶头视频| www国产在线视频色| 久久国产精品男人的天堂亚洲| 老司机午夜福利在线观看视频| 97人妻精品一区二区三区麻豆 | 国产黄片美女视频| 禁无遮挡网站| 亚洲精品中文字幕在线视频| 中亚洲国语对白在线视频| 999精品在线视频| 国产色视频综合| 91字幕亚洲| 一进一出抽搐动态| av中文乱码字幕在线| 制服诱惑二区| 老汉色∧v一级毛片| 欧美成人一区二区免费高清观看 | 国产伦人伦偷精品视频| 国产熟女xx| 亚洲av熟女| 亚洲一区二区三区不卡视频| 在线天堂中文资源库| 日本三级黄在线观看| 美女高潮喷水抽搐中文字幕| 精品久久久久久久久久免费视频| 日韩三级视频一区二区三区| 看片在线看免费视频| 国产亚洲欧美在线一区二区| 亚洲中文日韩欧美视频| 国产精品精品国产色婷婷| 国产免费男女视频| 午夜福利在线观看吧| 午夜两性在线视频| www.www免费av| 国内少妇人妻偷人精品xxx网站 | 黄色女人牲交| 久久人妻av系列| 天堂√8在线中文| 丁香六月欧美| 久久中文字幕一级| 又黄又爽又免费观看的视频| 91国产中文字幕| 亚洲国产欧美日韩在线播放| 99久久精品国产亚洲精品| 欧美av亚洲av综合av国产av| 亚洲avbb在线观看| 国产高清激情床上av| 国产精品一区二区免费欧美| 久久久国产精品麻豆| 国产伦在线观看视频一区| 国产精品 国内视频| 亚洲三区欧美一区| 午夜精品久久久久久毛片777| 日韩免费av在线播放| 狂野欧美激情性xxxx| 他把我摸到了高潮在线观看| 两个人看的免费小视频| 亚洲熟妇中文字幕五十中出| 人人妻人人看人人澡| 搡老妇女老女人老熟妇| 激情在线观看视频在线高清| 国产精品久久久久久人妻精品电影| 亚洲午夜精品一区,二区,三区| 久久久久九九精品影院| 91字幕亚洲| 国产亚洲欧美精品永久| 男女做爰动态图高潮gif福利片|