• <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窄帶干擾抑制算法
    亚洲熟女毛片儿| 成人精品一区二区免费| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久久精品吃奶| 天堂动漫精品| 精品一品国产午夜福利视频| 精品少妇一区二区三区视频日本电影| 正在播放国产对白刺激| 大码成人一级视频| 亚洲第一欧美日韩一区二区三区 | 亚洲七黄色美女视频| av在线播放免费不卡| 好男人电影高清在线观看| 一级a爱视频在线免费观看| 多毛熟女@视频| 成人国产av品久久久| 菩萨蛮人人尽说江南好唐韦庄| 国产精品 国内视频| 欧美黄色片欧美黄色片| 久久人人97超碰香蕉20202| 999久久久精品免费观看国产| 啦啦啦视频在线资源免费观看| 一本色道久久久久久精品综合| 中国美女看黄片| 久久久水蜜桃国产精品网| 在线播放国产精品三级| 国产成人av教育| 久久午夜亚洲精品久久| 老司机影院毛片| 国产免费av片在线观看野外av| 精品亚洲乱码少妇综合久久| 国产精品免费视频内射| 中文字幕高清在线视频| av网站在线播放免费| 69精品国产乱码久久久| 国产欧美日韩综合在线一区二区| 国产欧美日韩精品亚洲av| 丁香六月天网| 建设人人有责人人尽责人人享有的| 美女福利国产在线| 黑人猛操日本美女一级片| aaaaa片日本免费| 中亚洲国语对白在线视频| 日韩熟女老妇一区二区性免费视频| 人妻 亚洲 视频| 成人av一区二区三区在线看| 在线观看66精品国产| 丰满少妇做爰视频| 天堂8中文在线网| 亚洲精品乱久久久久久| 亚洲天堂av无毛| 免费观看av网站的网址| 日本黄色日本黄色录像| 午夜福利视频在线观看免费| 成人av一区二区三区在线看| 一边摸一边抽搐一进一出视频| 日日夜夜操网爽| www.精华液| 久9热在线精品视频| 亚洲av美国av| 亚洲国产成人一精品久久久| 下体分泌物呈黄色| 午夜福利视频精品| 欧美大码av| 视频区图区小说| cao死你这个sao货| 狠狠婷婷综合久久久久久88av| 在线观看一区二区三区激情| 日韩欧美国产一区二区入口| 亚洲色图av天堂| 免费在线观看视频国产中文字幕亚洲| 91老司机精品| 国产人伦9x9x在线观看| 91麻豆av在线| 久久午夜亚洲精品久久| 欧美日韩成人在线一区二区| av网站免费在线观看视频| 欧美在线一区亚洲| 十八禁人妻一区二区| 91国产中文字幕| 国产熟女午夜一区二区三区| 男人操女人黄网站| 免费人妻精品一区二区三区视频| 人成视频在线观看免费观看| 99精国产麻豆久久婷婷| 大陆偷拍与自拍| 久久久国产精品麻豆| 91成年电影在线观看| 欧美国产精品va在线观看不卡| 91av网站免费观看| 天堂俺去俺来也www色官网| 高清av免费在线| 丝袜美腿诱惑在线| 热re99久久精品国产66热6| 国产成人免费观看mmmm| 在线天堂中文资源库| 黑人巨大精品欧美一区二区mp4| 欧美日韩一级在线毛片| 一区二区日韩欧美中文字幕| 又黄又粗又硬又大视频| 国产熟女午夜一区二区三区| 一区福利在线观看| 一区二区三区国产精品乱码| 国产色视频综合| 国产高清激情床上av| 国产欧美日韩一区二区三区在线| 欧美成人午夜精品| 欧美性长视频在线观看| 美女福利国产在线| 久久国产精品影院| 免费黄频网站在线观看国产| 婷婷成人精品国产| 国产在视频线精品| 国产欧美日韩一区二区三区在线| 亚洲精品美女久久久久99蜜臀| 天天操日日干夜夜撸| 久久久久视频综合| 久久 成人 亚洲| 国产一区二区三区视频了| 两个人免费观看高清视频| 美女高潮到喷水免费观看| 人成视频在线观看免费观看| 日本wwww免费看| 亚洲精品国产色婷婷电影| 99九九在线精品视频| 少妇粗大呻吟视频| 免费少妇av软件| 女人久久www免费人成看片| 欧美日韩亚洲综合一区二区三区_| 日韩大码丰满熟妇| 极品人妻少妇av视频| 国产日韩一区二区三区精品不卡| 一本久久精品| www日本在线高清视频| 久久久国产成人免费| 欧美乱妇无乱码| 一区在线观看完整版| 99精品久久久久人妻精品| 青草久久国产| 国产黄色免费在线视频| cao死你这个sao货| 日韩有码中文字幕| 欧美日韩一级在线毛片| 中文字幕人妻熟女乱码| 91大片在线观看| 91大片在线观看| 欧美精品av麻豆av| 美女扒开内裤让男人捅视频| 777米奇影视久久| 国产成人啪精品午夜网站| 亚洲avbb在线观看| 精品欧美一区二区三区在线| 亚洲精品乱久久久久久| 色综合欧美亚洲国产小说| 黑人猛操日本美女一级片| 成人18禁在线播放| 国产高清国产精品国产三级| 高清av免费在线| 老司机福利观看| 十分钟在线观看高清视频www| 国产极品粉嫩免费观看在线| 国产av又大| 国产精品偷伦视频观看了| 99香蕉大伊视频| 精品国产一区二区三区久久久樱花| 啦啦啦 在线观看视频| 久久中文看片网| 亚洲av日韩精品久久久久久密| 国产欧美亚洲国产| 露出奶头的视频| 久久中文字幕人妻熟女| 欧美精品一区二区大全| 麻豆av在线久日| 久久99一区二区三区| 久久人妻熟女aⅴ| 丝袜人妻中文字幕| 99riav亚洲国产免费| 男人操女人黄网站| 日本vs欧美在线观看视频| www.熟女人妻精品国产| 国产成人精品久久二区二区免费| 中文字幕精品免费在线观看视频| 亚洲午夜精品一区,二区,三区| 亚洲中文字幕日韩| 精品亚洲成a人片在线观看| 国产av一区二区精品久久| 久久久国产精品麻豆| 亚洲国产欧美一区二区综合| 中国美女看黄片| 性少妇av在线| 久久国产精品人妻蜜桃| 丰满少妇做爰视频| 天堂俺去俺来也www色官网| 日韩免费av在线播放| 日韩成人在线观看一区二区三区| 人人妻人人爽人人添夜夜欢视频| 99国产综合亚洲精品| 夜夜爽天天搞| 国产黄频视频在线观看| 99热网站在线观看| 制服人妻中文乱码| 国产免费视频播放在线视频| 欧美成人午夜精品| 精品一区二区三区视频在线观看免费 | 国产精品一区二区在线观看99| 交换朋友夫妻互换小说| 麻豆国产av国片精品| 18禁黄网站禁片午夜丰满| 日韩中文字幕视频在线看片| av一本久久久久| 一区二区三区国产精品乱码| 一边摸一边抽搐一进一小说 | 国产男女超爽视频在线观看| 欧美日韩成人在线一区二区| 色婷婷av一区二区三区视频| 最新美女视频免费是黄的| 啦啦啦在线免费观看视频4| 9191精品国产免费久久| 一边摸一边做爽爽视频免费| 每晚都被弄得嗷嗷叫到高潮| 精品一区二区三卡| 啦啦啦 在线观看视频| 国产男女超爽视频在线观看| 国产黄色免费在线视频| 国产亚洲欧美精品永久| 丰满迷人的少妇在线观看| 国精品久久久久久国模美| 中文字幕精品免费在线观看视频| 丁香六月欧美| 国产精品一区二区免费欧美| 国产男女内射视频| 18禁黄网站禁片午夜丰满| 亚洲av欧美aⅴ国产| 成人亚洲精品一区在线观看| 五月开心婷婷网| 免费黄频网站在线观看国产| 国产精品98久久久久久宅男小说| 啦啦啦视频在线资源免费观看| 久久ye,这里只有精品| 搡老乐熟女国产| svipshipincom国产片| e午夜精品久久久久久久| 免费观看a级毛片全部| 久久天躁狠狠躁夜夜2o2o| 久久人妻熟女aⅴ| 亚洲av欧美aⅴ国产| 99久久精品国产亚洲精品| 久久99热这里只频精品6学生| 热99国产精品久久久久久7| 不卡一级毛片| 波多野结衣av一区二区av| 捣出白浆h1v1| 国产免费福利视频在线观看| 一进一出抽搐动态| 色视频在线一区二区三区| 亚洲国产欧美日韩在线播放| 啦啦啦 在线观看视频| 日本av免费视频播放| 亚洲avbb在线观看| 在线观看舔阴道视频| 黑人巨大精品欧美一区二区蜜桃| 免费少妇av软件| 国产午夜精品久久久久久| 日韩人妻精品一区2区三区| 黄色片一级片一级黄色片| 成人18禁高潮啪啪吃奶动态图| 在线 av 中文字幕| 亚洲伊人色综图| 国产精品免费一区二区三区在线 | 高清欧美精品videossex| 老汉色∧v一级毛片| 国产91精品成人一区二区三区 | 久久国产精品人妻蜜桃| 超碰成人久久| 在线亚洲精品国产二区图片欧美| avwww免费| 国产一区二区三区综合在线观看| 肉色欧美久久久久久久蜜桃| 精品国产一区二区久久| 久久影院123| 97人妻天天添夜夜摸| 在线观看人妻少妇| 亚洲精品美女久久av网站| 母亲3免费完整高清在线观看| 黄色片一级片一级黄色片| 天堂8中文在线网| 777久久人妻少妇嫩草av网站| 国产亚洲精品第一综合不卡| 热99久久久久精品小说推荐| 久久亚洲精品不卡| 国产真人三级小视频在线观看| 波多野结衣av一区二区av| 亚洲av成人不卡在线观看播放网| 国产视频一区二区在线看| av免费在线观看网站| 国产精品熟女久久久久浪| 99国产精品一区二区三区| 一区二区三区乱码不卡18| 91九色精品人成在线观看| 国产成人精品久久二区二区免费| 中文字幕高清在线视频| 日日爽夜夜爽网站| 精品久久久精品久久久| 国产精品久久久久久精品古装| 纯流量卡能插随身wifi吗| 国产高清视频在线播放一区| 欧美 日韩 精品 国产| 欧美亚洲日本最大视频资源| 超碰97精品在线观看| 在线观看www视频免费| 午夜两性在线视频| 精品一区二区三区视频在线观看免费 | 热99久久久久精品小说推荐| 成人手机av| 亚洲欧洲精品一区二区精品久久久| 黄色成人免费大全| 亚洲中文av在线| 精品国产一区二区久久| 亚洲中文日韩欧美视频| 国产成人欧美| 国产伦理片在线播放av一区| 国产精品久久久久久人妻精品电影 | 午夜成年电影在线免费观看| 欧美激情久久久久久爽电影 | 新久久久久国产一级毛片| 美女高潮到喷水免费观看| 日本wwww免费看| 精品高清国产在线一区| 纯流量卡能插随身wifi吗| 久久久久精品人妻al黑| 免费人妻精品一区二区三区视频| 久久影院123| 国产亚洲av高清不卡| 精品视频人人做人人爽| 国产精品自产拍在线观看55亚洲 | 久久精品国产亚洲av高清一级| 亚洲国产av影院在线观看| 久久99热这里只频精品6学生| 午夜激情av网站| 丁香欧美五月| 叶爱在线成人免费视频播放| 视频区图区小说| 亚洲色图综合在线观看| 国产视频一区二区在线看| 久久精品熟女亚洲av麻豆精品| 精品国产乱码久久久久久男人| 99香蕉大伊视频| 久久人人爽av亚洲精品天堂| 国产精品国产高清国产av | 极品少妇高潮喷水抽搐| 女性生殖器流出的白浆| 在线观看免费日韩欧美大片| 免费av中文字幕在线| 精品福利观看| 日韩精品免费视频一区二区三区| 天堂中文最新版在线下载| 免费不卡黄色视频| 精品一区二区三区av网在线观看 | 精品国产一区二区三区四区第35| 日韩视频一区二区在线观看| 99久久国产精品久久久| 亚洲欧美一区二区三区久久| av天堂在线播放| 国产精品自产拍在线观看55亚洲 | 女警被强在线播放| 男男h啪啪无遮挡| 97在线人人人人妻| 桃红色精品国产亚洲av| 可以免费在线观看a视频的电影网站| 一边摸一边抽搐一进一出视频| 国产亚洲精品第一综合不卡| 最近最新中文字幕大全免费视频| 多毛熟女@视频| av免费在线观看网站| 午夜福利乱码中文字幕| 久久免费观看电影| 亚洲国产中文字幕在线视频| 一级a爱视频在线免费观看| 在线永久观看黄色视频| 免费看十八禁软件| 多毛熟女@视频| 亚洲成人免费电影在线观看| 亚洲av美国av| 日韩有码中文字幕| av电影中文网址| av天堂在线播放| 中文字幕人妻丝袜制服| 成人三级做爰电影| 不卡一级毛片| 啦啦啦 在线观看视频| 精品一区二区三区四区五区乱码| 午夜福利影视在线免费观看| 国产在线观看jvid| 亚洲一码二码三码区别大吗| 成人18禁在线播放| 亚洲三区欧美一区| 成人亚洲精品一区在线观看| 最新在线观看一区二区三区| 夜夜骑夜夜射夜夜干| 丁香欧美五月| 麻豆成人av在线观看| 日日爽夜夜爽网站| 久久久久久久国产电影| 久久精品人人爽人人爽视色| 亚洲伊人久久精品综合| 热re99久久国产66热| 国产欧美日韩一区二区三区在线| 欧美日韩中文字幕国产精品一区二区三区 | 2018国产大陆天天弄谢| 五月天丁香电影| 久久天堂一区二区三区四区| 一区二区av电影网| 久久精品亚洲熟妇少妇任你| 国产高清激情床上av| 成人手机av| 亚洲自偷自拍图片 自拍| 欧美久久黑人一区二区| 自线自在国产av| a在线观看视频网站| 欧美日韩av久久| 叶爱在线成人免费视频播放| 欧美中文综合在线视频| 91成人精品电影| 亚洲一卡2卡3卡4卡5卡精品中文| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日韩亚洲高清精品| 51午夜福利影视在线观看| 夫妻午夜视频| 老熟妇仑乱视频hdxx| 纵有疾风起免费观看全集完整版| 老司机福利观看| 丰满饥渴人妻一区二区三| 国产免费现黄频在线看| 亚洲免费av在线视频| 97人妻天天添夜夜摸| 亚洲熟妇熟女久久| 日韩有码中文字幕| xxxhd国产人妻xxx| 水蜜桃什么品种好| 少妇精品久久久久久久| 性高湖久久久久久久久免费观看| 久久久久久久久久久久大奶| 美女福利国产在线| 99精品久久久久人妻精品| 一本综合久久免费| 日韩欧美三级三区| 色综合欧美亚洲国产小说| 亚洲少妇的诱惑av| 一级,二级,三级黄色视频| 精品卡一卡二卡四卡免费| 久久天堂一区二区三区四区| 亚洲avbb在线观看| 亚洲 国产 在线| 青草久久国产| 午夜福利影视在线免费观看| 国产精品二区激情视频| e午夜精品久久久久久久| 又黄又粗又硬又大视频| 亚洲av日韩精品久久久久久密| 午夜成年电影在线免费观看| 91成人精品电影| 在线观看免费午夜福利视频| 下体分泌物呈黄色| 露出奶头的视频| 欧美中文综合在线视频| 国产男靠女视频免费网站| 久久中文字幕人妻熟女| 可以免费在线观看a视频的电影网站| 亚洲精品成人av观看孕妇| 中文字幕另类日韩欧美亚洲嫩草| 大码成人一级视频| 伊人久久大香线蕉亚洲五| 男女边摸边吃奶| 丰满少妇做爰视频| 中文字幕人妻丝袜一区二区| 法律面前人人平等表现在哪些方面| 亚洲少妇的诱惑av| 国产国语露脸激情在线看| 999久久久精品免费观看国产| 国产精品 国内视频| 欧美日韩国产mv在线观看视频| 91麻豆精品激情在线观看国产 | 日日爽夜夜爽网站| 亚洲男人天堂网一区| 黄色丝袜av网址大全| 午夜精品国产一区二区电影| 一区二区三区精品91| 亚洲午夜精品一区,二区,三区| 精品国产一区二区三区久久久樱花| 午夜免费鲁丝| 国产一卡二卡三卡精品| 成人三级做爰电影| 欧美亚洲日本最大视频资源| 黄片大片在线免费观看| 手机成人av网站| 香蕉丝袜av| 91精品国产国语对白视频| 久热这里只有精品99| 一本—道久久a久久精品蜜桃钙片| 日本黄色视频三级网站网址 | 欧美日韩亚洲国产一区二区在线观看 | 久久久水蜜桃国产精品网| 国内毛片毛片毛片毛片毛片| 亚洲色图 男人天堂 中文字幕| 亚洲av欧美aⅴ国产| 一级片'在线观看视频| 精品国内亚洲2022精品成人 | videos熟女内射| 免费观看av网站的网址| 国内毛片毛片毛片毛片毛片| 这个男人来自地球电影免费观看| 国产高清videossex| 久久 成人 亚洲| 亚洲 国产 在线| www.999成人在线观看| 成人黄色视频免费在线看| 欧美亚洲 丝袜 人妻 在线| 亚洲中文日韩欧美视频| 日韩欧美一区视频在线观看| 精品午夜福利视频在线观看一区 | 欧美日韩黄片免| 国产精品av久久久久免费| 久久久精品94久久精品| 男女午夜视频在线观看| 久久九九热精品免费| 成人18禁在线播放| 在线亚洲精品国产二区图片欧美| 十八禁人妻一区二区| 国产亚洲精品一区二区www | 妹子高潮喷水视频| 久久中文字幕人妻熟女| 精品国产乱子伦一区二区三区| 999久久久国产精品视频| 99国产精品一区二区蜜桃av | 国产熟女午夜一区二区三区| 欧美精品亚洲一区二区| 十八禁网站网址无遮挡| www.精华液| 免费在线观看影片大全网站| 黄片小视频在线播放| 精品国产一区二区久久| 国产精品亚洲一级av第二区| 一级片免费观看大全| 窝窝影院91人妻| h视频一区二区三区| 久久久精品94久久精品| 丝袜喷水一区| 亚洲成国产人片在线观看| 妹子高潮喷水视频| 在线观看免费日韩欧美大片| 欧美性长视频在线观看| 久久中文字幕一级| 久久 成人 亚洲| 国产精品久久久人人做人人爽| 成人18禁在线播放| 亚洲国产成人一精品久久久| 久久精品亚洲av国产电影网| netflix在线观看网站| 啦啦啦免费观看视频1| 亚洲天堂av无毛| 最新美女视频免费是黄的| av线在线观看网站| 国产av精品麻豆| 熟女少妇亚洲综合色aaa.| 欧美+亚洲+日韩+国产| 精品国产一区二区三区久久久樱花| 夫妻午夜视频| 亚洲国产欧美日韩在线播放| 大陆偷拍与自拍| 夫妻午夜视频| 亚洲美女黄片视频| 亚洲av欧美aⅴ国产| 黄色怎么调成土黄色| 亚洲 国产 在线| 久久亚洲真实| 一级毛片女人18水好多| 波多野结衣av一区二区av| 后天国语完整版免费观看| 精品少妇内射三级| 亚洲精品久久成人aⅴ小说| 极品人妻少妇av视频| 国产精品99久久99久久久不卡| 国产欧美亚洲国产| 午夜久久久在线观看| 成人国产一区最新在线观看| 亚洲一区二区三区欧美精品| 超碰成人久久| 亚洲av日韩精品久久久久久密| 最新美女视频免费是黄的| 午夜福利一区二区在线看| 无人区码免费观看不卡 | 久久久久精品人妻al黑| 国产1区2区3区精品| 午夜福利,免费看| 啦啦啦中文免费视频观看日本| 亚洲精品在线观看二区| 欧美亚洲 丝袜 人妻 在线| 精品久久久精品久久久| 国产有黄有色有爽视频| 9色porny在线观看| 国产91精品成人一区二区三区 | 老熟妇乱子伦视频在线观看| 欧美激情极品国产一区二区三区| 亚洲精品自拍成人| 香蕉久久夜色| 大香蕉久久网| 亚洲第一欧美日韩一区二区三区 | 日日摸夜夜添夜夜添小说| 国产在线一区二区三区精| 久久性视频一级片| 9热在线视频观看99| 电影成人av| 国产男女内射视频| 一区二区日韩欧美中文字幕| 亚洲成a人片在线一区二区|