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

    一種振動能級梳狀濾波精確計算方法

    2015-10-24 05:45:51高慧中
    水下無人系統(tǒng)學報 2015年2期
    關鍵詞:梳狀倍頻程能級

    孫 濤,高慧中

    (中國船舶重工集團公司 第705研究所,陜西 西安,710075)

    一種振動能級梳狀濾波精確計算方法

    孫濤,高慧中

    (中國船舶重工集團公司 第705研究所,陜西 西安,710075)

    現(xiàn)有的1/3倍頻程處理方法都是以測試信號不包含明顯干擾為前提,因而在處理帶噪聲信號時存在一定先天不足。針對現(xiàn)場測量信號中經常存在規(guī)律性噪聲(噪聲線譜呈梳狀等間隔分布)的現(xiàn)象,提出了以Fourier核函數(shù)為基函數(shù)的梳狀帶阻濾波器與1/3倍頻程有機結合來消除規(guī)律性噪聲對結果影響的方法。該方法能夠在頻域下直接消除干擾并計算濾波后信號的振動能級,計算量小、精度高。通過對動力裝置功率試驗中出現(xiàn)規(guī)律性噪聲的振動信號進行合理修正,重新評估其振級水平,避免了多次重復試驗,節(jié)約了試驗費用。

    動力裝置; 振動噪聲; 信號處理; 1/3倍頻程

    0 引言

    噪聲作為人們不希望出現(xiàn)的聲音,卻大量的出現(xiàn)在人們的日常生活與生產活動中。在工業(yè)集中、人口密集的大中型城市,噪聲嚴重危害人民群眾的身心健康; 而在激烈的軍事對抗中,潛艇、戰(zhàn)機與軍艦的噪聲大小更直接與隱身性能直接相關,正是基于這樣的原因,越來越多的學者對其加以研究以實現(xiàn)降噪。聲音來源于振動,二者之間有很好的相關性,因而針對在工業(yè)現(xiàn)場難以使用測量儀器直接獲取的噪聲,通常采取“以振代噪”的方法進行評價。

    倍頻程方法是根據人耳對聲音的響應特性構造的一種頻域分析方法,具有譜線少、頻帶寬的特點,廣泛應用于聲學和振動的評估,并得到普遍認可[1]。一些學者從提高計算精度[2]、濾波器數(shù)字化[3]和增強處理功能[4]等方面對標準倍頻程算法進行改進,大大推動了該方法的工程應用。在1/3倍頻程方法的相關研究中,都是以測試信號為正常信號,不包含明顯干擾為前提,因而在處理帶噪聲信號時存在一定程度的先天不足。

    實際生產中,現(xiàn)場獲得的振動信號往往包含由諸如齒輪、軸承等外部系統(tǒng)引起的等間隔頻率成分,甚至不明干擾源引入的干擾譜線。如果不對這些干擾頻率成分進行處理就直接計算振動能級,必然導致計算結果偏大,影響評判的準確性。盡管目前已有多種成熟的梳狀濾波器構造方法[5-6],但不同的構造方法各具特點[7],尚無學者研究哪種梳狀濾波器適合與1/3倍頻程方法結合并實現(xiàn)它們的有機整合。

    該文提出以Fourier核函數(shù)為基函數(shù)的梳狀帶阻濾波器與1/3倍頻程有機結合來消除規(guī)律性噪聲對結果影響的方法,并通過模擬數(shù)據與實測數(shù)據驗證了該方法的可行性。

    1 算法流程

    目前的梳狀濾波器通常采用Morlet小波作為基礎單元來構造,這是因為Morlet小波作為高斯包絡下的一個復指數(shù)函數(shù),時域下收斂性好,可以通過改變其函數(shù)參數(shù)來改善邊線下降速度,得到較為理想的窄帶濾波器,獲取某些特征頻段,最終通過包絡檢波完成分析[8-9]。但是,Morlet小波在頻域下呈現(xiàn)“鐘”型曲線,頂部不平坦,且幅值一般不為1,并不能較好的滿足梳狀濾波需求[6]。因而,文中選擇以Fourier核函數(shù)構建梳狀濾波器,利用時域卷積等效于頻域相乘的原理,將原始信號的功率譜和構造出的梳狀帶阻濾波器相乘,濾除原始信號中的周期性干擾成分,然后再通過濾波后的功率譜計算1/3倍頻程譜,其簡要流程如圖1所示。該方法在頻域下直接消除干擾并計算濾波后信號的振動能級,計算量小、精度高,可用于修正因干擾導致的振動能級誤報。

    Fourier核函數(shù)的定義如下

    圖1  消除梳狀噪聲影響的1/3倍頻程方法流程圖Fig. 1 Flow chart of 1/3 octave spectrum method for eliminating comb-noise

    由于s(t)為自變量過零點的對稱函數(shù),若直接取t=0計算,無法獲得正確的數(shù)值解,文中提出可使用極限方法求解,即

    s(t)對應的頻域分布特性為

    其中ω0≥0,表示所構造帶通濾波器的截止頻率,其頻域分布特性曲線如圖2所示,取ω0=2π。

    圖2 Fourier核函數(shù)頻域波形Fig. 2 Waveform of Fourier kernel function in frequency domain

    可以看出,F(xiàn)ourier的頻域特性曲線相當于一個幅值為1的平頂帶通濾波器。根據Fourier變換的性質,通過設置頻移因子可構造出由多個帶通濾波器組成的基于Fourier核函數(shù)的梳狀濾波器

    構造出的hb(t)是以fb為等間隔,通帶帶寬為2ω0的梳狀濾波器,在頻域下將只能通過頻率為fb及其倍頻成分,還無法滿足濾除梳狀干擾頻率的要求,因此可在分析頻率fa范圍內構造一通頻帶通濾波器

    ha(t)的帶寬將覆蓋整個分析頻率范圍。

    在此基礎上,將構造得到的ha(t)與hb(t)進行Fourier變換得到Ha(ω)與Hb(ω)。根據Fourier變換的特性可知,梳狀濾波器Hc(ω)=Ha(ω)-Hb(ω)的頻域分布特性如圖3所示。

    圖3 梳狀帶阻濾波器頻域分布Fig. 3 Frequency domain distribution of comb bandstop filter

    Hc(ω)將僅對fb及其倍頻產生帶阻特性,而對其余頻率成分具有帶通特性。

    對原始信號x(t)的頻域值X(ω)進行梳狀濾波,得到對應的頻譜W(ω)后計算功率譜密度G(ω),其求解過程如下

    由此得到在以fc為中心頻率的第h個頻段[fl,fu]內功率譜的平均值為

    其中,n為頻段內譜線個數(shù)。若以a0為零分貝基準值,則ah的能級可表示為

    相比于原有濾波方法,此方法的濾波過程在頻域下直接進行,形象直觀,且得到的頻域結果可直接計算功率譜,無需再次進行傅里葉變換,避免了誤差的引入。

    2 工程應用與仿真

    2.1仿真測試

    構造含有間隔為10 Hz的梳狀譜噪聲的模擬信號,其數(shù)學表達式為

    圖4 模擬信號頻譜分布Fig. 4 Frequency spectrum distribution of simulated signal

    觀察模擬信號的頻譜X(ω)可以發(fā)現(xiàn),10,20,30,40 Hz處分布有等間隔的譜線,且幅值遠大于其余頻段的數(shù)值。以Fourier核函數(shù)為基礎,構造幅值為1中心頻率fc=10,20,30,40 Hz,帶寬B=2 Hz 的平頂帶阻濾波器Hc(ω)。

    在得到了梳狀帶阻濾波器的頻譜Hc(ω)后,根據時域卷積等效于頻域相乘的原理,濾波過程可直接在頻域下完成,濾波后的結果如圖5。

    圖5 濾波后的頻譜分布Fig. 5 Frequency spectrum distribution of filtered signal

    通過與圖4對比可看出周期性干擾譜線得到了明顯抑制,且其余頻段信息得到了完整保留。根據國際電工委員會對1/3倍頻程各頻帶范圍的規(guī)定[10]可以確定,4根干擾譜線位于中心頻率為10 Hz,20 Hz,31.5 Hz和40 Hz這4個1/3倍頻程頻帶內。原始信號的振動能級L1與濾波后信號的振動能級L2對比見圖6。為了便于觀察,該圖的1/3倍頻程中心頻率范圍取為8~40 Hz。可以清晰地看到,干擾譜線所在頻程的振動能級得到了顯著抑制。需要注意的是,由于干擾信號的30 Hz諧波與有用信號的頻率(29.5 Hz)十分接近,導致在濾波過程中有用信號也被濾除。針對這一問題,需要在密集分布區(qū)域單獨設計濾波器帶寬,避免有用信號被濾除。而對于更復雜的重合問題,可在預處理階段選擇更加有效的特征提取與分析技術對其加以分離,再完成濾波。

    2.2工程驗證

    動力裝置的功率試驗是在陸上全面檢驗動力裝置性能最有效的試驗手段。振動是功率試驗最重要的考核指標之一。受現(xiàn)場測試條件的限制,振動測量數(shù)據中不可避免地包含各種干擾。

    圖6 仿真信號振動能級對比圖Fig. 6 Contrast of vibration levels for simulated signal

    在某階段現(xiàn)場試驗中,連續(xù)幾次功率試驗評估點振動信號都包含間隔約100 Hz的周期性干擾,如圖7所示,并且干擾頻率的時間間隔和幅值不隨動力裝置速制的變化而變化。動力裝置機械運動激勵頻率等于發(fā)動機轉頻的分數(shù)倍或整數(shù)倍,當發(fā)動機的轉頻變化時,動力裝置的機械運動激勵頻率必然隨之變化。因此,判斷評估點振動頻譜的尖峰簇不是由動力裝置組件的機械運動激勵產生的,而是來自于外部干擾源。

    圖7 實測信號的時域與頻域曲線Fig. 7 Measured signal curves in time and frequency domains

    采用消除梳狀噪聲影響的1/3倍頻程方法,取帶寬B=10 Hz,間隔fb=100 Hz的梳狀濾波器對異常信號進行處理,濾波前后的振動能級曲線對比如圖8所示(濾波僅對規(guī)律譜線集中出現(xiàn)的頻段),可以看出,濾波后分貝值有所降低,符合振動要求(隨后的湖試測噪也間接證明消除干擾噪聲前振級超標的動力裝置振動處于正常范圍)。通過對功率試驗中出現(xiàn)類似特征的振動信號進行合理修正重新評估其振級水平,避免了多次重復試驗,節(jié)約了試驗費用。

    圖8 實測信號振動能級對比圖Fig. 8 Contrast of vibration levels for measured signal

    3 結束語

    隨著人們生活品質追求的日益提高和軍事對抗技術的不斷進步,噪聲的危害引起了主要工業(yè)部門的普遍關注。噪聲測試與評估是噪聲綜合治理的前提,評估方法的先進程度直接影響評估的準確性和效率。本文針對現(xiàn)場測量信號中經常存在規(guī)律性噪聲(噪聲線譜呈梳狀等間隔分布)的現(xiàn)象,將1/3倍頻程技術與梳狀帶阻濾波構造技術有機結合,并在此基礎上對求解過程加以優(yōu)化,提高了求解效率。該方法已廣泛應用于魚雷動力裝置及主要組件在復雜現(xiàn)場環(huán)境下的振動評估,大大提高了評估的準確性。

    [1]孫濤,高愛軍,王煒,等. 一種多測點變工況倍頻程快速處理方法[J]. 魚雷技術,2011,19(6): 455-458.

    Sun Tao,Gao Ai-jun,Wang Yi,et al. A Fast CPB Processing Method of Off-design Vibration and Noise Data from Multiple Channels[J]. Torpedo Technology,2011,19(6): 455-458.

    [2]成峰,周寶焜,程利青,等. 基于MATLAB的1/3倍頻程FIR數(shù)字濾波器設計[J]. 福州大學學報(自然科學版),2003,31(2): 160-165.

    Cheng Feng,Zhou Bao-kun,Cheng Li-qing,et al. Design of Third-Octave Band FIR Digital Filters Based on MATLAB[J]. Journal of Fuzhou University(Natural Science),2003,31(2): 160-165.

    [3]Flink J H,Bertrand J. Specrrum Analyzer Using Digital Filters: US,747053[P]. 1976-12-3.

    [4]孫濤,高愛軍,李留成. 直接獲得同工況比鄰的多測點變工況倍頻程譜陣的方法: 中國,ZL201110012751.1[P]. 2014-3-14.

    [5]李志勇,危韌勇,張濤. 基于Morlet組合小波的梳狀濾波與包絡檢波方法[J]. 中南大學學報(自然科學版),2006,37(2): 336-340.

    Li Zhi-yong,Wei Ren-yong,Zhang Tao. A Comb-filter and Envelope-demodulation Method Based on Combined Morlet Wavelets[J]. Journal of Central South University(Science and Technology),2006,37(2): 336-340.

    [6]趙俊龍,郭正剛,張志新,等. 梳狀濾波器在滾動軸承早期故障診斷中的應用[J]. 振動與沖擊,2008,27(12): 171-174.

    Zhao Jun-long,Guo Zheng-gang,Zhang Zhi-xin,et al. Application of Comb Fil the in Fault Early Diagnosis of Rolling Element Bearing[J]. Journal of Vibration and Shock,2008,27(12): 171-174.

    [7]侯者非,楊杰,王克成,等. 改進梳狀包絡檢波方法在滾動軸承故障診斷中的應用[J]. 煤礦機械,2010,31(2): 229-232.

    Hou Zhe-fei,Yang Jie,Wang Ke-cheng,et al. Novel Comb-filter and Envelope-demodulation Method Applied in Fault Diagnosis of Rolling Element Bearing[J]. Coal Mine Machinery,2010,31(2): 229-232.

    [8]Bruce A,Donoho D,Gao Hongye. Wavlete analysis[C]// IEEE Spectrum,1996: 26-35.

    [9]Grossmann A,Morlet J. Decomposition of Hardy Functions into Square Integrable Wavelets of Constant Shape[J]. SIAM Journal on Mathematical Analysis,1984,15(4): 14.

    [10]王濟. MATLAB在振動信號處理中的應用[M].北京: 中國水利水電出版社,2006: 134-135.

    (責任編輯: 陳曦)

    An Accurate Computation Method of Vibration Level Based on Comb Band-stop Filter

    SUN Tao,GAO Hui-zhong
    (The 705 Research Institute,China Shipbuilding Industry Corporation,Xi′an 710075,China)

    The existing 1/3 octave spectrum approaches are usually supposed to test signals without interference. Obviously,this strategy is no more effective when confronted with the signals containing band noise. In this paper,aiming at the regular noise interference,such as comb noise,the Fourier kernel function is employed to construct a comb bandstop filter,and further the filter is combined with the 1/3 octave spectrum to remove the effects exerted by noise. The proposed strategy can remove the interference in frequency domain and calculate the filtered vibration level more efficiently and conveniently,compared with the existing approaches. For the vibration signals with regular noise interference in the power test of a propulsion system,proper revision and consequent reevaluation can avoid repeated tests and save test costs.

    propulsion system; vibration noise; signal process; 1/3 octave spectrum

    TJ630.32; TB523

    A

    1673-1948(2015)02-0119-05

    2014-11-18;

    2015-01-08.

    海軍裝備重大專項預研(40110),船舶預研支撐技術基金(14J4.4.1)

    孫濤(1977-),男,博士,高級工程師,主要研究方向為測試測量技術、振源識別與故障診斷、多物理場耦合仿真、動力試驗技術及其系統(tǒng)設計.

    猜你喜歡
    梳狀倍頻程能級
    超聲導波目標模式的激勵方法研究
    振動與沖擊(2024年4期)2024-03-02 02:02:14
    “拼、搶、快、優(yōu)”,展現(xiàn)錢塘“高能級”擔當
    杭州(2023年3期)2023-04-03 07:22:04
    一種抗干擾變電站1/3倍頻程噪聲測量方法*
    應用聲學(2022年6期)2022-11-23 10:51:14
    噪聲聲譜控制算法的研究
    1500萬年前的南酸棗化石長啥樣等
    提升醫(yī)學教育能級 培養(yǎng)拔尖創(chuàng)新人才
    糾纏Ξ-型三能級原子與糾纏腔場相互作用熵的糾纏演化
    常規(guī)倍頻程純音測聽聽閾無異常的耳鳴患者的半倍頻程頻率測試結果分析
    幾種三分之一倍頻程中心頻率定義方法的比較
    PVDF梳狀換能器接收非線性蘭姆波的實驗研究?
    應用聲學(2019年3期)2019-07-25 07:19:58
    欧美最黄视频在线播放免费| 精品久久蜜臀av无| 亚洲成人久久爱视频| 日本三级黄在线观看| 欧美日韩国产亚洲二区| 18禁观看日本| 黄色视频不卡| 禁无遮挡网站| 高清在线国产一区| 亚洲熟女毛片儿| avwww免费| 嫩草影院精品99| 在线观看日韩欧美| 黑人欧美特级aaaaaa片| 亚洲 欧美 日韩 在线 免费| 国产又色又爽无遮挡免费看| 亚洲欧美一区二区三区黑人| 校园春色视频在线观看| 黄色成人免费大全| 国产私拍福利视频在线观看| tocl精华| 久久久久久久久中文| 少妇粗大呻吟视频| 亚洲精品美女久久av网站| 高清毛片免费观看视频网站| 欧美一级a爱片免费观看看 | 欧美性猛交╳xxx乱大交人| 男人的好看免费观看在线视频 | 国产欧美日韩精品亚洲av| 在线观看免费日韩欧美大片| 亚洲在线自拍视频| 老熟妇仑乱视频hdxx| 99精品欧美一区二区三区四区| 久久精品成人免费网站| 91老司机精品| 国产爱豆传媒在线观看 | 国产午夜精品论理片| 国产av在哪里看| 视频区欧美日本亚洲| 制服诱惑二区| 天堂av国产一区二区熟女人妻 | 国模一区二区三区四区视频 | 老司机在亚洲福利影院| 免费电影在线观看免费观看| 最近在线观看免费完整版| 精品第一国产精品| 12—13女人毛片做爰片一| 色在线成人网| 亚洲专区中文字幕在线| 久久中文看片网| 一个人观看的视频www高清免费观看 | 国产精品九九99| 麻豆成人午夜福利视频| 成人一区二区视频在线观看| 老司机靠b影院| 国产精品一区二区三区四区久久| 美女 人体艺术 gogo| 丁香六月欧美| 久久久久久人人人人人| 久久久精品欧美日韩精品| 热99re8久久精品国产| 国产精品野战在线观看| 欧美av亚洲av综合av国产av| 亚洲aⅴ乱码一区二区在线播放 | 亚洲精品国产精品久久久不卡| 欧美色视频一区免费| 欧洲精品卡2卡3卡4卡5卡区| 国产亚洲av嫩草精品影院| 麻豆成人av在线观看| 老司机深夜福利视频在线观看| 在线国产一区二区在线| 成人一区二区视频在线观看| 国产伦一二天堂av在线观看| 国产精品99久久99久久久不卡| 国产一区二区三区视频了| 免费电影在线观看免费观看| 黄色片一级片一级黄色片| 国产激情偷乱视频一区二区| 亚洲自拍偷在线| 久久人妻福利社区极品人妻图片| 国产精品永久免费网站| 999精品在线视频| 长腿黑丝高跟| 啦啦啦韩国在线观看视频| 午夜两性在线视频| 九色成人免费人妻av| bbb黄色大片| 波多野结衣巨乳人妻| 神马国产精品三级电影在线观看 | 久久久精品大字幕| 两个人的视频大全免费| 女生性感内裤真人,穿戴方法视频| 99久久国产精品久久久| www日本在线高清视频| 一a级毛片在线观看| 母亲3免费完整高清在线观看| √禁漫天堂资源中文www| 久久久久国产精品人妻aⅴ院| 国产不卡一卡二| 1024视频免费在线观看| 亚洲精品av麻豆狂野| 久久久久久人人人人人| 舔av片在线| 欧美色欧美亚洲另类二区| 亚洲专区字幕在线| 91麻豆精品激情在线观看国产| 黄片大片在线免费观看| 欧美日本视频| 国产精品久久久人人做人人爽| 这个男人来自地球电影免费观看| 一区二区三区国产精品乱码| 村上凉子中文字幕在线| 黄频高清免费视频| 国产精品一区二区免费欧美| 亚洲成人精品中文字幕电影| 十八禁人妻一区二区| 免费在线观看日本一区| 午夜影院日韩av| 亚洲av第一区精品v没综合| 美女扒开内裤让男人捅视频| 亚洲av中文字字幕乱码综合| 麻豆一二三区av精品| 免费看美女性在线毛片视频| 欧美日韩中文字幕国产精品一区二区三区| 亚洲国产精品sss在线观看| 亚洲片人在线观看| 国产亚洲精品久久久久久毛片| 国产精品亚洲一级av第二区| 亚洲美女黄片视频| 久久久国产欧美日韩av| 嫩草影视91久久| 国产精品影院久久| 黄色 视频免费看| 亚洲专区中文字幕在线| 日本黄色视频三级网站网址| 久久精品91蜜桃| 激情在线观看视频在线高清| 日本熟妇午夜| 女人爽到高潮嗷嗷叫在线视频| 午夜激情av网站| 久久国产精品影院| 变态另类成人亚洲欧美熟女| 黑人巨大精品欧美一区二区mp4| 在线a可以看的网站| 日本黄色视频三级网站网址| 欧美又色又爽又黄视频| 日韩av在线大香蕉| 天堂av国产一区二区熟女人妻 | 国产亚洲欧美在线一区二区| 欧美日韩亚洲国产一区二区在线观看| 国产精品综合久久久久久久免费| 久久精品国产亚洲av高清一级| 伦理电影免费视频| 国产三级在线视频| 久久精品91无色码中文字幕| 国产成人精品久久二区二区免费| 国产不卡一卡二| 成人手机av| 亚洲午夜理论影院| 久久精品aⅴ一区二区三区四区| 真人一进一出gif抽搐免费| 中文字幕av在线有码专区| 国产成人系列免费观看| 国产熟女午夜一区二区三区| 老司机在亚洲福利影院| 精品久久久久久久毛片微露脸| 久久久久国产一级毛片高清牌| 一级黄色大片毛片| 亚洲狠狠婷婷综合久久图片| 日本一本二区三区精品| 在线视频色国产色| 两性夫妻黄色片| 免费一级毛片在线播放高清视频| 欧美黑人精品巨大| 在线看三级毛片| 亚洲精品一卡2卡三卡4卡5卡| 床上黄色一级片| 女人被狂操c到高潮| 国产精品免费视频内射| 精品国产乱码久久久久久男人| 日本在线视频免费播放| 久久精品人妻少妇| 中国美女看黄片| 精品高清国产在线一区| 亚洲,欧美精品.| 老司机午夜十八禁免费视频| 一边摸一边做爽爽视频免费| 亚洲精品中文字幕一二三四区| 黑人欧美特级aaaaaa片| 在线视频色国产色| 女警被强在线播放| 国产伦在线观看视频一区| 99久久精品热视频| xxxwww97欧美| 国产亚洲av高清不卡| 亚洲成人国产一区在线观看| 亚洲激情在线av| 日本熟妇午夜| 欧美三级亚洲精品| 村上凉子中文字幕在线| 色综合站精品国产| 亚洲一卡2卡3卡4卡5卡精品中文| 久久热在线av| 精品不卡国产一区二区三区| 久久精品夜夜夜夜夜久久蜜豆 | 午夜精品久久久久久毛片777| 久久人妻av系列| 夜夜看夜夜爽夜夜摸| 丰满的人妻完整版| 久久久久国内视频| 制服人妻中文乱码| 一区二区三区激情视频| 99国产综合亚洲精品| 免费在线观看成人毛片| 国产午夜福利久久久久久| 国产精品av久久久久免费| 成人18禁高潮啪啪吃奶动态图| 男女午夜视频在线观看| 最近在线观看免费完整版| www日本在线高清视频| 变态另类丝袜制服| 欧美国产日韩亚洲一区| 国产成人系列免费观看| 很黄的视频免费| 国产精品一区二区精品视频观看| 国产精品自产拍在线观看55亚洲| 日日夜夜操网爽| 一级a爱片免费观看的视频| 免费看十八禁软件| 国内精品一区二区在线观看| 国产黄a三级三级三级人| 国产精品久久电影中文字幕| 亚洲真实伦在线观看| 久久久久国内视频| 老司机在亚洲福利影院| 午夜亚洲福利在线播放| 日韩欧美在线乱码| 老汉色av国产亚洲站长工具| 亚洲第一电影网av| 啦啦啦免费观看视频1| 欧美日韩一级在线毛片| 国产精品1区2区在线观看.| 亚洲中文av在线| 国产片内射在线| 亚洲成人精品中文字幕电影| 一进一出抽搐动态| 国产精品1区2区在线观看.| 天天添夜夜摸| 天堂av国产一区二区熟女人妻 | 免费看日本二区| 大型av网站在线播放| 九色国产91popny在线| 老熟妇乱子伦视频在线观看| 久久久久久久久久黄片| 国内精品久久久久久久电影| 国产av麻豆久久久久久久| 一区二区三区高清视频在线| 视频区欧美日本亚洲| aaaaa片日本免费| 亚洲最大成人中文| 草草在线视频免费看| 毛片女人毛片| 看片在线看免费视频| 久久久久亚洲av毛片大全| 午夜亚洲福利在线播放| svipshipincom国产片| 日本黄色视频三级网站网址| АⅤ资源中文在线天堂| 很黄的视频免费| 少妇熟女aⅴ在线视频| 国内揄拍国产精品人妻在线| 亚洲五月婷婷丁香| 国产成人精品久久二区二区免费| 999久久久国产精品视频| 欧美黄色片欧美黄色片| 亚洲九九香蕉| 可以在线观看的亚洲视频| 欧美国产日韩亚洲一区| 一级毛片精品| 一进一出抽搐gif免费好疼| 久久久久久大精品| 国产精华一区二区三区| 九色国产91popny在线| 夜夜看夜夜爽夜夜摸| 久久久国产欧美日韩av| 亚洲黑人精品在线| 国产午夜精品论理片| 国产一区二区在线av高清观看| 成人18禁在线播放| 亚洲美女视频黄频| 欧洲精品卡2卡3卡4卡5卡区| 午夜a级毛片| 国产精品亚洲av一区麻豆| 国产爱豆传媒在线观看 | 国产精品免费视频内射| 99久久无色码亚洲精品果冻| 99国产精品一区二区三区| 91成年电影在线观看| 中文字幕久久专区| 日韩欧美一区二区三区在线观看| 亚洲avbb在线观看| 男女下面进入的视频免费午夜| x7x7x7水蜜桃| 好男人电影高清在线观看| 在线观看免费日韩欧美大片| 亚洲熟女毛片儿| 母亲3免费完整高清在线观看| 免费在线观看视频国产中文字幕亚洲| 欧美日韩乱码在线| 国产亚洲精品综合一区在线观看 | 国产亚洲精品综合一区在线观看 | 免费在线观看视频国产中文字幕亚洲| 少妇裸体淫交视频免费看高清 | 亚洲男人天堂网一区| 亚洲精品国产一区二区精华液| 色在线成人网| 狂野欧美激情性xxxx| 日本五十路高清| 一进一出抽搐动态| 又粗又爽又猛毛片免费看| 国产欧美日韩一区二区精品| 12—13女人毛片做爰片一| 青草久久国产| 我的老师免费观看完整版| 亚洲五月婷婷丁香| 在线a可以看的网站| 亚洲aⅴ乱码一区二区在线播放 | 国产一区二区在线av高清观看| 给我免费播放毛片高清在线观看| 亚洲中文av在线| 欧美色视频一区免费| 丁香六月欧美| 亚洲自拍偷在线| 搡老妇女老女人老熟妇| 日韩欧美国产在线观看| 亚洲精品一区av在线观看| 日韩大码丰满熟妇| 精品一区二区三区av网在线观看| 国产亚洲av高清不卡| 国产av麻豆久久久久久久| 男人舔奶头视频| 亚洲无线在线观看| 99热这里只有精品一区 | 国产人伦9x9x在线观看| 久久久国产欧美日韩av| 18禁黄网站禁片午夜丰满| 小说图片视频综合网站| 19禁男女啪啪无遮挡网站| 国产69精品久久久久777片 | 韩国av一区二区三区四区| 国产午夜精品久久久久久| 久久国产精品人妻蜜桃| 99热只有精品国产| 国产人伦9x9x在线观看| 亚洲av美国av| 亚洲 欧美一区二区三区| 亚洲精品色激情综合| 九色国产91popny在线| 国产精品,欧美在线| 国产高清videossex| 亚洲精品久久成人aⅴ小说| 欧美成人免费av一区二区三区| 五月伊人婷婷丁香| 法律面前人人平等表现在哪些方面| av在线播放免费不卡| 校园春色视频在线观看| 在线永久观看黄色视频| 男女那种视频在线观看| 久久久水蜜桃国产精品网| 高潮久久久久久久久久久不卡| 国产高清视频在线观看网站| 国产精品野战在线观看| 久久久精品欧美日韩精品| 成人精品一区二区免费| 亚洲九九香蕉| 国内精品一区二区在线观看| 中亚洲国语对白在线视频| 欧美精品亚洲一区二区| 黄片大片在线免费观看| 一区二区三区高清视频在线| 99热6这里只有精品| 国产精品影院久久| 久久精品国产99精品国产亚洲性色| 亚洲欧美日韩无卡精品| 国产激情久久老熟女| 午夜a级毛片| 久久伊人香网站| 亚洲国产看品久久| 高清在线国产一区| 亚洲18禁久久av| 精品一区二区三区av网在线观看| 久久香蕉国产精品| 亚洲一区二区三区不卡视频| cao死你这个sao货| 亚洲成人免费电影在线观看| 一级作爱视频免费观看| 国产av在哪里看| 女人被狂操c到高潮| 国产在线观看jvid| www.熟女人妻精品国产| 88av欧美| 精品高清国产在线一区| 少妇裸体淫交视频免费看高清 | 嫩草影院精品99| 亚洲一区二区三区不卡视频| 亚洲黑人精品在线| 亚洲 欧美 日韩 在线 免费| 欧美人与性动交α欧美精品济南到| 欧美一区二区精品小视频在线| 无人区码免费观看不卡| 亚洲美女黄片视频| 日韩欧美在线二视频| 午夜免费观看网址| 国产精品电影一区二区三区| 成在线人永久免费视频| 亚洲一区中文字幕在线| 嫁个100分男人电影在线观看| 一区二区三区激情视频| 亚洲人与动物交配视频| 亚洲国产精品999在线| 亚洲男人的天堂狠狠| 久久精品国产清高在天天线| 国产精品电影一区二区三区| 亚洲九九香蕉| 欧美精品亚洲一区二区| 成年人黄色毛片网站| 每晚都被弄得嗷嗷叫到高潮| 亚洲成a人片在线一区二区| 丁香欧美五月| 国产片内射在线| 亚洲18禁久久av| 波多野结衣巨乳人妻| 我要搜黄色片| 午夜激情av网站| 黄片大片在线免费观看| 好男人电影高清在线观看| 国产69精品久久久久777片 | 日韩大码丰满熟妇| 亚洲av五月六月丁香网| 国产伦人伦偷精品视频| 国产午夜精品久久久久久| e午夜精品久久久久久久| 亚洲av成人一区二区三| 欧美黑人欧美精品刺激| 女人被狂操c到高潮| 大型黄色视频在线免费观看| 亚洲成av人片免费观看| 欧美激情久久久久久爽电影| 色综合站精品国产| 久久 成人 亚洲| 丝袜美腿诱惑在线| 757午夜福利合集在线观看| 欧美成人一区二区免费高清观看 | 丁香欧美五月| 国产aⅴ精品一区二区三区波| 亚洲熟女毛片儿| 毛片女人毛片| 亚洲全国av大片| 亚洲免费av在线视频| 日本 av在线| 国产成人av激情在线播放| 久久久国产精品麻豆| 不卡av一区二区三区| 最近最新免费中文字幕在线| 99久久久亚洲精品蜜臀av| 亚洲成a人片在线一区二区| 国产精品香港三级国产av潘金莲| 精品第一国产精品| 亚洲狠狠婷婷综合久久图片| 国产精品99久久99久久久不卡| 制服诱惑二区| 亚洲自拍偷在线| 亚洲第一电影网av| 久久久久久国产a免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 国产高清videossex| 九色成人免费人妻av| 女人爽到高潮嗷嗷叫在线视频| 色播亚洲综合网| 亚洲精品中文字幕在线视频| 国产精品亚洲av一区麻豆| 少妇被粗大的猛进出69影院| 成人一区二区视频在线观看| 精品国产乱码久久久久久男人| 国产乱人伦免费视频| 亚洲色图 男人天堂 中文字幕| 久久欧美精品欧美久久欧美| 日韩欧美在线二视频| 亚洲精品美女久久av网站| 久久精品亚洲精品国产色婷小说| 日韩高清综合在线| 男人舔奶头视频| 欧美另类亚洲清纯唯美| 99在线人妻在线中文字幕| 又黄又粗又硬又大视频| 又大又爽又粗| av中文乱码字幕在线| svipshipincom国产片| 黄色视频不卡| 色老头精品视频在线观看| 一本一本综合久久| 少妇人妻一区二区三区视频| 成人18禁高潮啪啪吃奶动态图| 国产精品爽爽va在线观看网站| 国内精品一区二区在线观看| 麻豆国产av国片精品| 少妇人妻一区二区三区视频| xxx96com| 国产精品精品国产色婷婷| 在线观看66精品国产| 99久久国产精品久久久| 久久婷婷人人爽人人干人人爱| 黄色毛片三级朝国网站| 啪啪无遮挡十八禁网站| 夜夜看夜夜爽夜夜摸| 青草久久国产| 黄色 视频免费看| 99国产极品粉嫩在线观看| 亚洲欧美日韩高清专用| 欧美日韩黄片免| 日韩大码丰满熟妇| 小说图片视频综合网站| 制服诱惑二区| 国产精品久久久av美女十八| 午夜精品一区二区三区免费看| 91av网站免费观看| 精品国产超薄肉色丝袜足j| 久久久久久亚洲精品国产蜜桃av| 精品乱码久久久久久99久播| 成人永久免费在线观看视频| a在线观看视频网站| www.自偷自拍.com| 国产av一区在线观看免费| 90打野战视频偷拍视频| 麻豆成人av在线观看| 国内精品久久久久久久电影| 欧美午夜高清在线| 欧美乱码精品一区二区三区| 午夜福利视频1000在线观看| 欧美av亚洲av综合av国产av| 一级毛片女人18水好多| 岛国在线免费视频观看| 日韩欧美 国产精品| 老司机午夜福利在线观看视频| 天堂动漫精品| 成人国产综合亚洲| 午夜精品久久久久久毛片777| 国产又色又爽无遮挡免费看| 日日摸夜夜添夜夜添小说| 看片在线看免费视频| 大型av网站在线播放| 日本免费一区二区三区高清不卡| 欧美成人免费av一区二区三区| av欧美777| 久久久久精品国产欧美久久久| 黄片小视频在线播放| 日韩欧美国产在线观看| 欧美黄色片欧美黄色片| 日日夜夜操网爽| 黄频高清免费视频| 黑人巨大精品欧美一区二区mp4| 1024手机看黄色片| 欧美+亚洲+日韩+国产| 老司机福利观看| 国产成年人精品一区二区| 91在线观看av| a级毛片在线看网站| 一区二区三区高清视频在线| 国产日本99.免费观看| 午夜a级毛片| 怎么达到女性高潮| 久久久久久久午夜电影| 国产精品乱码一区二三区的特点| 日本 欧美在线| 久久久久久免费高清国产稀缺| 一进一出抽搐动态| 丁香欧美五月| 亚洲天堂国产精品一区在线| 国产91精品成人一区二区三区| 精品第一国产精品| 老汉色av国产亚洲站长工具| 在线十欧美十亚洲十日本专区| 国产精品美女特级片免费视频播放器 | 黄频高清免费视频| 国产亚洲精品av在线| 亚洲人成77777在线视频| 亚洲成av人片在线播放无| 桃红色精品国产亚洲av| 国产野战对白在线观看| 免费观看精品视频网站| 国产精品久久久久久人妻精品电影| 久久久久亚洲av毛片大全| xxx96com| 亚洲国产精品成人综合色| 亚洲精品在线观看二区| 50天的宝宝边吃奶边哭怎么回事| 欧美成人性av电影在线观看| 日本免费一区二区三区高清不卡| 级片在线观看| 亚洲精品中文字幕一二三四区| 好男人电影高清在线观看| 久久国产精品人妻蜜桃| 人成视频在线观看免费观看| 精品人妻1区二区| 怎么达到女性高潮| 1024视频免费在线观看| 国产在线观看jvid| 中文字幕高清在线视频| 午夜精品在线福利| 国产又色又爽无遮挡免费看| 黑人操中国人逼视频| 久9热在线精品视频| 一进一出好大好爽视频| 在线十欧美十亚洲十日本专区| 高清毛片免费观看视频网站|