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

    船舶螺旋槳空化噪聲非均勻調(diào)制特性及其應(yīng)用

    2017-07-01 19:04:10劉啟軍邱家興程玉勝
    艦船科學(xué)技術(shù) 2017年6期
    關(guān)鍵詞:線譜頻帶空化

    劉啟軍,邱家興,程玉勝

    (海軍潛艇學(xué)院,山東 青島 266199)

    船舶螺旋槳空化噪聲非均勻調(diào)制特性及其應(yīng)用

    劉啟軍,邱家興,程玉勝

    (海軍潛艇學(xué)院,山東 青島 266199)

    從螺旋槳空化噪聲譜統(tǒng)計(jì)模型出發(fā),分析研究螺旋槳空化噪聲非均勻調(diào)制特性的機(jī)理和規(guī)律,并利用大量海上實(shí)錄船舶噪聲進(jìn)行驗(yàn)證。從形成機(jī)理看,認(rèn)為非均勻調(diào)制特性主要由船尾形狀、螺旋槳結(jié)構(gòu)以及船舶航行工況共同決定的尾流分布所導(dǎo)致,不同型號(hào)船舶的非均勻調(diào)制特性具有個(gè)體性;從艦船噪聲信號(hào)能量譜看,調(diào)制的非均勻特性,主要受能量譜峰值頻率的影響。最后基于船舶噪聲的非均勻調(diào)制特性,提出一種基于多子帶自適應(yīng)加權(quán)的DEMON增強(qiáng)算法,并利用實(shí)際數(shù)據(jù)進(jìn)行驗(yàn)證。

    螺旋槳空化噪聲;非均勻調(diào)制;調(diào)制度;多子帶;線譜增強(qiáng)

    0 引 言

    船舶螺旋槳一旦空化,其輻射噪聲就會(huì)成為船舶航行時(shí)最主要的噪聲,螺旋槳在船尾非均勻尾流中旋轉(zhuǎn)并產(chǎn)生空化的過程中,空化噪聲能量又會(huì)被周期性調(diào)制,表現(xiàn)為在空化調(diào)制譜中有明顯的具有倍數(shù)關(guān)系的多根線譜,利用這些調(diào)制線譜可以判斷出物理意義明確的螺旋槳轉(zhuǎn)速和槳葉數(shù)等特征,是目前船舶噪聲識(shí)別中非常重要的識(shí)別特征[1–3]。

    螺旋槳空化噪聲在不同頻帶上的調(diào)制程度不均勻,但這種非均勻調(diào)制特性的規(guī)律如何,調(diào)制強(qiáng)弱隨頻帶分布關(guān)系,乃至這種非均勻調(diào)制特性的利用等問題,目前已公開發(fā)表的文獻(xiàn)中未見有詳細(xì)研究和結(jié)論。本文將從螺旋槳空化噪聲譜統(tǒng)計(jì)模型出發(fā),將螺旋槳空化噪聲非均勻調(diào)制特性的機(jī)理和規(guī)律作為研究重點(diǎn)并進(jìn)一步考慮其運(yùn)用。

    1 螺旋槳空化噪聲非均勻調(diào)制特性

    螺旋槳在非均勻尾流中轉(zhuǎn)動(dòng),每個(gè)時(shí)刻都有大量隨機(jī)產(chǎn)生和崩潰的空化氣泡。對(duì)于螺旋槳空化噪聲譜,本文綜合采用文獻(xiàn)[4]中的群空泡崩潰統(tǒng)計(jì)模型和文獻(xiàn)[3]中假設(shè)空泡群平均最大半徑服從正態(tài)分布的特點(diǎn),進(jìn)行推導(dǎo)[5–7],得到螺旋槳空化噪聲能量譜的表達(dá)式為[4]:

    式中:C為一比例常數(shù);t為時(shí)間;ω為頻率;KR為一取決于氣體含量和環(huán)境壓力等因素的常數(shù)因子,海況和目標(biāo)航行深度等的變化會(huì)引起其變化;N為某一時(shí)刻發(fā)生崩潰的空泡個(gè)數(shù);Rmax為氣泡崩潰開始時(shí)最大半徑;ρ為水的密度;P∞為無窮遠(yuǎn)處水的壓力(靜壓力);為同一時(shí)刻大量空泡崩潰最大半徑Rmax的均值;σ為方差。

    為了表征能量譜幅度被調(diào)制的程度,定義mp為能量譜在一個(gè)時(shí)間周期T內(nèi),能量譜幅值變化的深度,即幅度調(diào)制深度。螺旋槳空化噪聲可近似視為周期性局部平穩(wěn)隨機(jī)過程,短時(shí)間內(nèi)(幾個(gè)螺旋槳轉(zhuǎn)動(dòng)周期內(nèi))螺旋槳轉(zhuǎn)速、海洋環(huán)境、尾流等相對(duì)穩(wěn)定,mp可視為隨時(shí)間穩(wěn)定不變,因此mp可僅看作是頻率ω的函數(shù),記為mp(ω);在任一時(shí)間周期T內(nèi),中心時(shí)刻為t,能量譜的時(shí)間平均值可視為相對(duì)穩(wěn)定,記為G-(ω,t),能量譜最大值和最小值分別出現(xiàn)在t1和t2時(shí)刻,則

    根據(jù)信號(hào)調(diào)制理論,能量譜被周期性幅度調(diào)制,可表示為:

    推得能量譜幅度調(diào)制深度為:

    從式(8)可看出,能量譜幅度調(diào)制深度完全可以由能量譜最大值與最小值之比(最大相對(duì)起伏)來表示,用分貝數(shù)表示并記為:

    通過計(jì)算,結(jié)合圖1容易得出:

    1)4 個(gè)參數(shù)中,僅N(t)變化時(shí),mG(ω)在不同頻率ω處為常數(shù),即空泡個(gè)數(shù)的變化引起了能量譜在不同頻帶間的均勻調(diào)制;而其他 3 個(gè)參數(shù)(t),KR,σ的變化均會(huì)引起能量譜在不同頻帶間的非均勻調(diào)制;

    2)引起能量譜的調(diào)制深度mG(ω)較大主要因素為N(t),(t),而KR,σ的變化引起的調(diào)制深度相對(duì)較小,尤其是在能量譜峰值頻率以上的中高頻段上更是如此;

    從以上分析可知,螺旋槳空化在不同頻帶的非均勻調(diào)制特性主要是受空泡平均最大半徑的影響,而影響空泡平均最大半徑的因素是螺旋槳槳葉和海水之間的相對(duì)來流速度[4,9-10],因此,最根本的原因是由船尾形狀、螺旋槳結(jié)構(gòu)以及船舶航行工況共同決定的尾流分布所導(dǎo)致的[3]。由于不同型號(hào)船舶的船尾形狀、螺旋槳結(jié)構(gòu)具有個(gè)體性,因此其也具有特有的螺旋槳空化噪聲非均勻調(diào)制特性。

    此外,據(jù)圖1的仿真分析還可得出,從螺旋槳空化輻射噪聲能量譜看,調(diào)制的非均勻特性,主要受能量譜峰值頻率的影響。在峰值頻率的左側(cè)低頻段,調(diào)制深度最高,在峰值頻率的右側(cè)一小段頻帶內(nèi),調(diào)制深度最低,在高頻段調(diào)制深度基本保持不變,略高于最低調(diào)制深度。分析這種現(xiàn)象的原因:一是螺旋槳空化噪聲能量譜本身受空泡半徑的影響最大,其次是當(dāng)空泡半徑變大時(shí),不僅空化噪聲譜線整體上升,而且峰值頻率還要向低頻移動(dòng),勢(shì)必造成在空化噪聲能量譜在低頻段調(diào)制深度高于高頻段,并在緊靠峰值頻率而又高于峰值頻率一側(cè)具有最低的調(diào)制深度[1]。

    2 實(shí)錄船舶噪聲非均勻調(diào)制特性分析

    下面通過海上實(shí)錄的船舶噪聲信號(hào)進(jìn)行調(diào)制譜分析,驗(yàn)證以上理論是否符合實(shí)際情況。為了定量描述噪聲調(diào)制度大小,采用文獻(xiàn)[4]提出的噪聲調(diào)制度σs計(jì)算方法。雖然螺旋槳空化噪聲不是船舶輻射噪聲的全部,但一旦螺旋槳發(fā)生空化,就將成為船舶輻射噪聲的最主要部分,因此可以用海上實(shí)錄的高工況航行船舶噪聲近似研究螺旋槳空化噪聲的調(diào)制特性。對(duì)船舶噪聲進(jìn)行DEMON分析,通過分析計(jì)算軸頻及其諧波線譜、葉頻及其諧波線譜、干擾線譜之間的大小關(guān)系,判定噪聲的調(diào)制度大小[11]。在分析噪聲不同頻帶的非均勻調(diào)制特性時(shí),首先分頻段解調(diào),然后對(duì)各頻段解調(diào)獲得的DEMON譜分別計(jì)算調(diào)制度大小,最后組成調(diào)制度矢量,得到非均勻調(diào)制特性,計(jì)算過程如圖2所示。

    從實(shí)錄船舶噪聲庫中隨機(jī)選取工況較高的目標(biāo)513個(gè),每個(gè)目標(biāo)噪聲長(zhǎng)度40 s~3 min不等,對(duì)每個(gè)目標(biāo)按照?qǐng)D2所示步驟分別計(jì)算得到一個(gè)非均勻調(diào)制特性矢量圖,從統(tǒng)計(jì)結(jié)果看,本文認(rèn)為船舶噪聲的非均勻調(diào)制特性雖然看似無章可循,但其實(shí)是有一定規(guī)律的,大多數(shù)噪聲基本滿足本文上節(jié)所述的“從螺旋槳空化輻射噪聲能量譜看…在峰值頻率的左側(cè)低頻段,調(diào)制深度最高,在峰值頻率的右側(cè)一小段頻帶內(nèi),調(diào)制深度最低”的特點(diǎn)。

    圖3為某大型船舶噪聲非均勻調(diào)制特性圖。圖中:①為目標(biāo)功率譜圖縱坐標(biāo)為信號(hào)頻率0~10 kHz;②為目標(biāo)全頻帶噪聲DEMON譜圖,橫坐標(biāo)為調(diào)制頻率,范圍為0~50 Hz;③為目標(biāo)DEMON譜歷程圖,橫坐標(biāo)為調(diào)制頻率,縱坐標(biāo)為信號(hào)頻率;④為目標(biāo)非均勻調(diào)制度矢量圖,縱坐標(biāo)為信號(hào)頻率,橫坐標(biāo)為各頻率子帶的調(diào)制度大小,dB。

    該信號(hào)的功率譜圖較好滿足標(biāo)準(zhǔn)船舶噪聲功率譜特點(diǎn):峰值頻率很低,在峰值頻率以下的頻段,由于傳感器濾波的原因,看不到幅值的衰減,在峰值頻率以上的中高頻段大約以–6 dB/倍頻程的速率下降;其非均勻性調(diào)制度矢量圖基本符合上節(jié)分析的結(jié)論特點(diǎn),因信號(hào)頻率基本位于峰值頻率以上,其非均勻調(diào)制特點(diǎn)則只具有峰值以上頻率的調(diào)制特性,在整個(gè)信號(hào)頻率范圍內(nèi),噪聲調(diào)制度逐漸遞增,在高頻部分遞增的幅度較小,而在靠近且低于1 kHz的信號(hào)頻段,由于恰巧處于峰值頻率的后沿,調(diào)制度最弱。

    3 基于多子帶自適應(yīng)加權(quán)的DEMON線譜增強(qiáng)

    傳統(tǒng)DEMON譜分析,一般是依據(jù)噪聲接收的工作頻段,設(shè)置一個(gè)固定的較為通用的解調(diào)頻段,所有噪聲信號(hào)均在該頻段內(nèi)濾波、檢波、DEMON譜分析,并據(jù)此進(jìn)行線譜提取和特征分析。而根據(jù)本文第1節(jié)結(jié)論,船舶噪聲的調(diào)制特性在不同頻帶上是非均勻的,那么在DEMON譜分析中,如果僅僅在調(diào)制度高的信號(hào)頻段進(jìn)行解調(diào),獲得的DEMON譜線譜清晰度是否會(huì)更高呢,這是本節(jié)要研究的問題。

    根據(jù)第1節(jié)的分析,對(duì)艦船輻射噪聲的解調(diào),實(shí)際上為對(duì)某個(gè)頻段ω1噪聲能量Gi(t)的頻譜分析,即DEMON為:

    則全頻段的DEMON譜可以近似表示為N個(gè)子頻帶DEMON的譜值積累:

    考慮到將調(diào)制度大小差異較大的子帶DEMON譜值積累一起,與僅將調(diào)制度高的子帶DEMON譜值積累一起相比,調(diào)制特性應(yīng)該會(huì)有所減弱,因此本文在自動(dòng)計(jì)算各子帶調(diào)制度的基礎(chǔ)上,僅將調(diào)制度高的子帶DEMON譜積累作為最終的DEMON譜,如圖4所示。根據(jù)大量數(shù)據(jù)統(tǒng)計(jì)結(jié)果給定一個(gè)調(diào)制度可用門限?,若某個(gè)子帶的調(diào)制度σs(ωi)大于門限?,則判定該子帶權(quán)值為1,否則為0,最終DEMON譜D(f)為:

    表1 多子帶自適應(yīng)加權(quán)的DEMON譜線譜增強(qiáng)效果Tab.1 The effect of DEMON line spectrum strength based onmulti sub-band adaptive weighted

    仍然采用第2節(jié)使用的513個(gè)海上實(shí)錄船舶噪聲數(shù)據(jù)進(jìn)行分析計(jì)算,分別進(jìn)行全頻段解調(diào)和采用不同加權(quán)門限?時(shí)的多子帶自適應(yīng)加權(quán)解調(diào)。統(tǒng)計(jì)513個(gè)目標(biāo)的平均調(diào)制度和線譜平均增強(qiáng)幅度,如表1所示。2個(gè)典型目標(biāo)噪聲的DEMON譜線譜增強(qiáng)前后對(duì)比如圖5所示。

    從實(shí)際噪聲數(shù)據(jù)的分析結(jié)果看,頻帶的非均勻調(diào)制特性對(duì)DEMON譜線譜幅值的計(jì)算結(jié)果影響還是比較明顯的,尤其是當(dāng)噪聲的調(diào)制特性僅分布在比較窄的頻帶內(nèi),影響更甚。通過多子帶自適應(yīng)加權(quán),可以獲得較優(yōu)的DEMON譜線譜清晰度。另外,權(quán)限?取較低值(5 dB)時(shí)就已經(jīng)取得的較好的增強(qiáng)效果,而繼續(xù)把權(quán)限值?升高,增強(qiáng)的幅度減緩,這是因?yàn)闄?quán)限?給定一個(gè)較低值時(shí),已經(jīng)將沒有調(diào)制的頻帶濾掉,基本上消除了無調(diào)制頻帶對(duì)調(diào)制特性的影響。

    4 結(jié) 語

    本文主要研究了螺旋槳空化非均勻調(diào)制特性的機(jī)理、規(guī)律和應(yīng)用。首先從船舶螺旋槳空化噪聲譜模型出發(fā),仿真分析了螺旋槳空化非均勻調(diào)制特性;然后利用大量海上實(shí)錄船舶噪聲,統(tǒng)計(jì)了實(shí)際船舶噪聲非均勻調(diào)制特性,驗(yàn)證了以上結(jié)論的合理性;最后基于船舶噪聲的非均勻調(diào)制特性,提出了一種基于多子帶自適應(yīng)加權(quán)的DEMON增強(qiáng)算法,試驗(yàn)表明其效果良好。主要形成以下結(jié)論:

    1)從形成機(jī)理看,螺旋槳空化在不同頻帶的非均勻調(diào)制特性主要是由船尾形狀、螺旋槳結(jié)構(gòu)以及船舶航行工況共同決定的尾流分布所導(dǎo)致的,不同型號(hào)船舶的非均勻調(diào)制特性具有個(gè)體性。

    2)從艦船噪聲信號(hào)能量譜看,調(diào)制的非均勻特性,主要受能量譜峰值頻率的影響。在峰值頻率的前沿低頻段,調(diào)制深度最高,在峰值頻率的后沿一小段頻帶內(nèi),調(diào)制深度最低,在高頻段調(diào)制深度基本保持不變,略高于最低調(diào)制深度。

    3)船舶噪聲的非均勻調(diào)制特性對(duì)DEMON譜線譜幅值影響明顯,尤其是當(dāng)噪聲的調(diào)制特性僅分布在比較窄的頻帶內(nèi)時(shí),影響更甚。通過多子帶自適應(yīng)加權(quán),對(duì)DEMON譜線譜增強(qiáng)有明顯效果,可以獲得更優(yōu)的線譜清晰度與譜圖結(jié)構(gòu)特征。

    [ 1 ]史廣智. 艦船噪聲調(diào)制譜諧波族結(jié)構(gòu)特性理論分析[J]. 聲學(xué)學(xué)報(bào), 2007, 32(1): 19–25.

    SHI Guang-zhi. Theoretical analysis of the structure law of ship radiated-noise demodulation spectrum harmonic clan feature[J]. Acta Acustica, 2007, 32(1): 19–25.

    [ 2 ]史廣智. 基于模板匹配的艦船螺旋槳噪聲特征精細(xì)分析[J].聲學(xué)學(xué)報(bào), 2009, 34(5): 401–407.

    SHI Guang-zhi. A method for ship propeller blade-number recognition based on template matching[J]. Acta Acustica, 2009, 34(5): 401–407.

    [ 3 ]陶篤純. 螺旋槳空化噪聲譜[J]. 聲學(xué)學(xué)報(bào), 1982, 7(6): 344–351.

    TAO Du-chun. Spectrum of propeller cavitation noise[J]. Acta Acustica, 1982, 7(6): 344–351.

    [ 4 ]邱家興, 程玉勝, 等. 船舶大側(cè)斜螺旋槳空化噪聲調(diào)制譜諧波結(jié)構(gòu)[J]. 聲學(xué)技術(shù)增刊-中國聲學(xué)學(xué)會(huì)水聲學(xué)分會(huì)2015年學(xué)術(shù)會(huì)議論文集, 2015, 34(2): 364–368.

    QIU Jia-xing, CHENG Yu-sheng, et al. Demodulation spectrum harmonic line structure of high-skewed propeller cavitation noise[J]. Technical Acoustics, 2015, 34(2): 364–368.

    [ 5 ]戚定滿, 魯傳敬, 何友聲. 空化噪聲譜特性研究[J]. 振動(dòng)與沖擊, 1999, 18(3): 32–36.

    QI Ding-man, LU Chuang-jing, HE You-sheng. An investigation on the spectra of bubble noise[J]. Journal of Vibration and Shock, 1999, 18(3): 32–36.

    [ 6 ]ROSS D. Mechanica of Underwater Noise[M]. 北京: 《水下噪聲原理翻譯組》譯. 1983: 277–394.

    [ 7 ]陶篤純. 艦船噪聲節(jié)奏的研究(I)—數(shù)學(xué)模型及功率譜密度[J]. 聲學(xué)學(xué)報(bào), 1983, 8(2): 65–76.

    TAO Du-chun. A study on ship radiated noise rhythms(i)-mathematical model and power spectrum densidy[J]. Acta Acustica, 1983, 8(2): 65–76.

    [ 8 ]陶篤純. 艦船噪聲節(jié)奏的研究(II)—自相關(guān)函數(shù)及節(jié)奏信息的提取[J]. 聲學(xué)學(xué)報(bào), 1983, 8(5): 280–289.

    TAO Du-chun. A study on ship radiated noise rhythms(ii)-autocorrelation function and extraction of rhythm information[J]. Acta Acustica, 1983, 8(5): 280–289.

    [ 9 ]魏以邁. 均流與非均流中螺旋槳空泡噪聲的研究[J]. 聲學(xué)技術(shù), 1983, (4): 39–50.

    WEI Yi-mai. Study on propeller cavitation noise under the steady and unsteady hydrodynamic[J]. Ship Science and Technology, 1983, (4): 39–50.

    [10]何友聲, 王國強(qiáng). 螺旋槳激振力[M]. 上海: 上海交通大學(xué)出版社, 1987: 13–67.

    HE You-sheng, Wang Guo-qiang. Propeller exciting force[M]. Shanghai: Shanghai Jiao Tong University Press, 1987: 13–67.

    [11]邱家興, 程玉勝, 等. 船舶噪聲DEMON譜質(zhì)量評(píng)估方法研究[J]. 艦船科學(xué)技術(shù), 2014, 36(9): 46–49.

    QIU Jia-xing, CHENG Yu-sheng, et al. Study of ship-radiated noise demon quality assessment methods[J]. Ship Science and Technology, 2014, 36(9): 46–49.

    Non-uniform modulation feature of ship propeller cavitation noiseand its application

    LIU Qi-jun, QIU Jia-xing, CHENG Yu-sheng
    (Navy Submarine Academy, Qingdao 266199, China)

    Based on the statistical model of ship propeller cavitation noise spectrum, the mechanism and laws of nonuniform modulation feature of propeller cavitation noise are analyzed and studied, and also be verified by a large number of ship radiated noise. From thepoint of formed mechanism, it shows that non-uniform modulation feature is mainly formed by the distribution of wake flow, which is formed by the combined action of stern shape, propeller structure and working condition. Non-uniform modulation feature vary in different class ship. Then from the point of ship radiated noise spectrum, nonuniform modulation feature is mainly decided by the energy spectrum peak frequency. Finally, an application method of nonuniformity envelope is proposed, that is the technology of DEMON line spectrum strength based on multi sub-band adaptive weighted, and it is also verified by real date.

    propeller cavitation noise;non-uniform modulation;modulation degree;multi sub-band;line spectrum strength

    TB533

    A

    1672 – 7619(2017)06 – 0018 – 05

    10.3404/j.issn.1672 – 7619.2017.06.004

    2017 – 03 – 14

    劉啟軍(1975 – ),男,碩士,講師,研究方向?yàn)樗暷繕?biāo)識(shí)別。

    猜你喜歡
    線譜頻帶空化
    功率超聲作用下鋼液中空化泡尺寸的演變特性
    鋼鐵釩鈦(2023年5期)2023-11-17 08:48:34
    UUV慣導(dǎo)系統(tǒng)多線譜振動(dòng)抑制研究
    Wi-Fi網(wǎng)絡(luò)中5G和2.4G是什么?有何區(qū)別?
    單音及部分頻帶干擾下DSSS系統(tǒng)性能分析
    雙頻帶隔板極化器
    三維扭曲水翼空化現(xiàn)象CFD模擬
    不同運(yùn)動(dòng)形式下水物相互作用空化數(shù)值模擬
    調(diào)諧放大器通頻帶的計(jì)算及應(yīng)用
    “”維譜在艦船輻射噪聲線譜提取中的應(yīng)用
    基于隱馬爾可夫模型的線譜跟蹤技術(shù)
    欧美日韩中文字幕国产精品一区二区三区 | 窝窝影院91人妻| 成年版毛片免费区| 午夜福利视频在线观看免费| 国产成人精品在线电影| 免费观看a级毛片全部| 欧洲精品卡2卡3卡4卡5卡区| 在线十欧美十亚洲十日本专区| 每晚都被弄得嗷嗷叫到高潮| 欧美黑人精品巨大| 老鸭窝网址在线观看| 两人在一起打扑克的视频| 久久青草综合色| 国产精品亚洲一级av第二区| 在线永久观看黄色视频| 久久天躁狠狠躁夜夜2o2o| 亚洲第一欧美日韩一区二区三区| www日本在线高清视频| 国产高清视频在线播放一区| a级片在线免费高清观看视频| 午夜亚洲福利在线播放| 国产成人欧美在线观看 | 国产一区二区三区在线臀色熟女 | 亚洲色图综合在线观看| 亚洲熟女毛片儿| 午夜福利欧美成人| 一级片免费观看大全| 中出人妻视频一区二区| 欧美不卡视频在线免费观看 | 国产精品自产拍在线观看55亚洲 | 国产在线观看jvid| 亚洲av熟女| 国产精品亚洲av一区麻豆| 热re99久久国产66热| 少妇粗大呻吟视频| 国产一区二区激情短视频| 国产精品成人在线| 免费在线观看日本一区| 国产精品久久久久久精品古装| 一进一出好大好爽视频| 久久九九热精品免费| 成人18禁在线播放| 成人特级黄色片久久久久久久| 一进一出好大好爽视频| 国产av又大| 亚洲欧洲精品一区二区精品久久久| 亚洲一区中文字幕在线| 久久狼人影院| 最近最新免费中文字幕在线| 国产片内射在线| 精品福利观看| 国产淫语在线视频| 久久精品国产综合久久久| 日韩欧美国产一区二区入口| 午夜亚洲福利在线播放| 色综合欧美亚洲国产小说| 99国产精品一区二区三区| videos熟女内射| 成人亚洲精品一区在线观看| 三级毛片av免费| 亚洲中文日韩欧美视频| 人人妻,人人澡人人爽秒播| 国产亚洲欧美在线一区二区| av片东京热男人的天堂| a级毛片黄视频| 丁香欧美五月| 久久精品91无色码中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 久热爱精品视频在线9| 超碰97精品在线观看| 夜夜夜夜夜久久久久| 99热网站在线观看| 久久久国产一区二区| 国产人伦9x9x在线观看| 欧美大码av| 国产成人免费无遮挡视频| 欧美精品亚洲一区二区| 免费黄频网站在线观看国产| 欧美黑人精品巨大| 久久久久精品人妻al黑| 亚洲精品中文字幕一二三四区| 久久精品亚洲av国产电影网| 人人妻人人澡人人看| 一级片'在线观看视频| 一区二区三区国产精品乱码| 男人舔女人的私密视频| 色94色欧美一区二区| 99久久国产精品久久久| 免费在线观看视频国产中文字幕亚洲| 国产主播在线观看一区二区| 麻豆av在线久日| 多毛熟女@视频| 国产精品免费视频内射| 精品人妻熟女毛片av久久网站| 国产av又大| 亚洲 欧美一区二区三区| 国产精品综合久久久久久久免费 | 亚洲色图av天堂| 午夜福利一区二区在线看| 最新的欧美精品一区二区| 成人影院久久| 久久久久久久国产电影| 丝袜美腿诱惑在线| 法律面前人人平等表现在哪些方面| 9热在线视频观看99| 黄色女人牲交| 在线观看免费日韩欧美大片| 日日摸夜夜添夜夜添小说| 久久精品国产a三级三级三级| 亚洲性夜色夜夜综合| 国产成人免费观看mmmm| 97人妻天天添夜夜摸| 在线永久观看黄色视频| 久久九九热精品免费| 国产激情欧美一区二区| 国产在线观看jvid| 母亲3免费完整高清在线观看| 多毛熟女@视频| 亚洲精品成人av观看孕妇| 欧美激情 高清一区二区三区| 国产av又大| 在线观看免费高清a一片| 日本五十路高清| 香蕉丝袜av| av网站在线播放免费| 亚洲专区国产一区二区| 日韩熟女老妇一区二区性免费视频| 最新在线观看一区二区三区| 国产精品久久久久成人av| 精品久久久久久,| 69av精品久久久久久| x7x7x7水蜜桃| 国产精品一区二区精品视频观看| 美女高潮到喷水免费观看| 亚洲成a人片在线一区二区| 99国产精品一区二区三区| 久久精品国产亚洲av香蕉五月 | 国产在线一区二区三区精| 欧美激情高清一区二区三区| 黑人猛操日本美女一级片| 久久久精品区二区三区| 俄罗斯特黄特色一大片| 久久人人爽av亚洲精品天堂| 亚洲全国av大片| 久久中文字幕一级| 黄色丝袜av网址大全| 亚洲国产欧美日韩在线播放| а√天堂www在线а√下载 | 狠狠狠狠99中文字幕| 中文字幕最新亚洲高清| 99国产极品粉嫩在线观看| 亚洲欧美一区二区三区黑人| 午夜免费鲁丝| 一级a爱片免费观看的视频| a级毛片在线看网站| 一级毛片女人18水好多| 精品视频人人做人人爽| 99久久国产精品久久久| 亚洲欧美日韩另类电影网站| 成在线人永久免费视频| 夫妻午夜视频| 老熟女久久久| 亚洲av熟女| 亚洲在线自拍视频| 久久国产精品男人的天堂亚洲| 18禁裸乳无遮挡免费网站照片 | 99精品欧美一区二区三区四区| bbb黄色大片| 欧美国产精品va在线观看不卡| 黄色女人牲交| av片东京热男人的天堂| 婷婷成人精品国产| 18禁美女被吸乳视频| 久久午夜综合久久蜜桃| 久久精品国产a三级三级三级| 精品福利永久在线观看| 啦啦啦在线免费观看视频4| 欧美精品高潮呻吟av久久| 午夜福利在线观看吧| 黄色a级毛片大全视频| 午夜免费成人在线视频| 中文欧美无线码| 亚洲午夜理论影院| 日韩欧美三级三区| 精品久久蜜臀av无| 亚洲专区中文字幕在线| 这个男人来自地球电影免费观看| 欧美人与性动交α欧美软件| av在线播放免费不卡| 成人手机av| 99久久99久久久精品蜜桃| 人妻 亚洲 视频| 国产国语露脸激情在线看| 色老头精品视频在线观看| 深夜精品福利| 性少妇av在线| 色综合婷婷激情| 欧美 日韩 精品 国产| 另类亚洲欧美激情| 91成年电影在线观看| 99riav亚洲国产免费| 人人妻人人澡人人爽人人夜夜| av网站免费在线观看视频| 久久ye,这里只有精品| 亚洲欧美激情在线| 国产深夜福利视频在线观看| 午夜福利免费观看在线| e午夜精品久久久久久久| 亚洲国产欧美网| 999精品在线视频| 日韩欧美三级三区| 中文字幕人妻丝袜制服| 国产成人免费观看mmmm| 亚洲专区字幕在线| 村上凉子中文字幕在线| 日本a在线网址| 成年女人毛片免费观看观看9 | 久久久久国产精品人妻aⅴ院 | 国产区一区二久久| 黄色毛片三级朝国网站| 成年版毛片免费区| 日韩欧美国产一区二区入口| 美女高潮喷水抽搐中文字幕| 天天操日日干夜夜撸| 久久久久久久午夜电影 | 在线观看www视频免费| 日韩 欧美 亚洲 中文字幕| 免费在线观看影片大全网站| 每晚都被弄得嗷嗷叫到高潮| 777久久人妻少妇嫩草av网站| 久久国产精品男人的天堂亚洲| 国产精品免费一区二区三区在线 | 嫩草影视91久久| a级片在线免费高清观看视频| 久久人妻福利社区极品人妻图片| 在线永久观看黄色视频| 亚洲av美国av| 丝袜美腿诱惑在线| 国产精品成人在线| 国产男女内射视频| 在线观看免费高清a一片| 午夜两性在线视频| 午夜免费成人在线视频| 欧美日韩中文字幕国产精品一区二区三区 | 国产欧美日韩一区二区三区在线| 亚洲精品美女久久久久99蜜臀| 国产伦人伦偷精品视频| 啪啪无遮挡十八禁网站| 久久九九热精品免费| 欧美丝袜亚洲另类 | 免费在线观看亚洲国产| 两人在一起打扑克的视频| 大码成人一级视频| 亚洲午夜理论影院| 国产精品乱码一区二三区的特点 | 久久婷婷成人综合色麻豆| 悠悠久久av| 亚洲全国av大片| 国产欧美亚洲国产| 久久精品成人免费网站| 午夜老司机福利片| √禁漫天堂资源中文www| 亚洲av成人一区二区三| 婷婷成人精品国产| 欧美黑人欧美精品刺激| 亚洲av日韩在线播放| 国产成人免费观看mmmm| 咕卡用的链子| 国产在线一区二区三区精| 在线免费观看的www视频| 久久久久久久久久久久大奶| 国产欧美日韩综合在线一区二区| 多毛熟女@视频| 一区二区三区国产精品乱码| 亚洲欧美激情在线| 村上凉子中文字幕在线| 免费在线观看完整版高清| 国产男女超爽视频在线观看| 久久久久久亚洲精品国产蜜桃av| 国产99久久九九免费精品| 国产精品乱码一区二三区的特点 | 国产黄色免费在线视频| 婷婷丁香在线五月| 精品久久久精品久久久| 悠悠久久av| 欧美午夜高清在线| 久久久国产精品麻豆| 热99久久久久精品小说推荐| 99精国产麻豆久久婷婷| 中文欧美无线码| 欧美成狂野欧美在线观看| 夫妻午夜视频| 亚洲色图综合在线观看| 亚洲五月色婷婷综合| 国产xxxxx性猛交| 69av精品久久久久久| 99riav亚洲国产免费| 啦啦啦在线免费观看视频4| 亚洲中文字幕日韩| 国产淫语在线视频| 国产精品久久久久久人妻精品电影| 18禁裸乳无遮挡动漫免费视频| 99国产精品99久久久久| 麻豆乱淫一区二区| 天堂俺去俺来也www色官网| 精品一品国产午夜福利视频| 国产精品 国内视频| 精品欧美一区二区三区在线| 极品人妻少妇av视频| 欧美国产精品va在线观看不卡| 涩涩av久久男人的天堂| 亚洲av成人av| 国产极品粉嫩免费观看在线| 精品一区二区三区av网在线观看| 美女视频免费永久观看网站| 免费在线观看完整版高清| 最新的欧美精品一区二区| 黄色怎么调成土黄色| 国产亚洲精品一区二区www | 久久中文看片网| 男女之事视频高清在线观看| 很黄的视频免费| 亚洲av成人一区二区三| 免费少妇av软件| 少妇裸体淫交视频免费看高清 | 两个人看的免费小视频| 啦啦啦视频在线资源免费观看| 精品亚洲成a人片在线观看| 国产精品九九99| 99精品久久久久人妻精品| 香蕉久久夜色| av中文乱码字幕在线| 男女床上黄色一级片免费看| 久久精品亚洲av国产电影网| 久久久国产一区二区| 精品国产亚洲在线| 另类亚洲欧美激情| 久久性视频一级片| 国产97色在线日韩免费| 欧美 亚洲 国产 日韩一| 午夜福利乱码中文字幕| 欧美精品人与动牲交sv欧美| 最新在线观看一区二区三区| 怎么达到女性高潮| 国产精品免费视频内射| 在线观看66精品国产| 国产精品二区激情视频| 亚洲一区二区三区欧美精品| 丰满人妻熟妇乱又伦精品不卡| 亚洲熟女毛片儿| 亚洲片人在线观看| 欧美黄色片欧美黄色片| 久久久精品国产亚洲av高清涩受| 国产人伦9x9x在线观看| 久久久久久免费高清国产稀缺| 91麻豆精品激情在线观看国产 | 久久精品人人爽人人爽视色| 大片电影免费在线观看免费| 美女国产高潮福利片在线看| 成人av一区二区三区在线看| 极品教师在线免费播放| 久久久国产成人精品二区 | 夜夜爽天天搞| www.熟女人妻精品国产| 中文字幕av电影在线播放| 亚洲精品成人av观看孕妇| 欧美黑人欧美精品刺激| 亚洲人成电影免费在线| 一本综合久久免费| 国产xxxxx性猛交| 亚洲九九香蕉| 久久九九热精品免费| 亚洲一区二区三区不卡视频| 午夜激情av网站| 岛国毛片在线播放| 视频区图区小说| 日本vs欧美在线观看视频| 丁香六月欧美| 一级毛片精品| 亚洲色图av天堂| 999久久久精品免费观看国产| 日韩欧美免费精品| 欧美黄色淫秽网站| 国产精品99久久99久久久不卡| 宅男免费午夜| 免费看十八禁软件| 午夜福利影视在线免费观看| 精品电影一区二区在线| 岛国毛片在线播放| 免费av中文字幕在线| 高清在线国产一区| 看黄色毛片网站| 久久精品aⅴ一区二区三区四区| 久久亚洲精品不卡| 欧美激情极品国产一区二区三区| 黄色片一级片一级黄色片| 午夜福利免费观看在线| 精品福利永久在线观看| 91成年电影在线观看| 久久久国产成人免费| 手机成人av网站| 一级a爱片免费观看的视频| 别揉我奶头~嗯~啊~动态视频| 久久国产精品人妻蜜桃| 久久久精品区二区三区| 大陆偷拍与自拍| 欧美成狂野欧美在线观看| 国产熟女午夜一区二区三区| 一区二区三区激情视频| 久久人妻熟女aⅴ| 在线观看免费高清a一片| 久久精品国产亚洲av香蕉五月 | 美女高潮喷水抽搐中文字幕| 免费看十八禁软件| 91九色精品人成在线观看| 亚洲美女黄片视频| 麻豆国产av国片精品| 精品福利观看| 黄片大片在线免费观看| 国产三级黄色录像| 国产片内射在线| 久久国产精品人妻蜜桃| 亚洲熟女毛片儿| 国产亚洲精品久久久久5区| 大香蕉久久成人网| 在线国产一区二区在线| 成年动漫av网址| 欧美日韩亚洲高清精品| 亚洲av成人av| 最新美女视频免费是黄的| tube8黄色片| bbb黄色大片| 亚洲精品一卡2卡三卡4卡5卡| 80岁老熟妇乱子伦牲交| 一级片'在线观看视频| 叶爱在线成人免费视频播放| av一本久久久久| 极品少妇高潮喷水抽搐| 麻豆国产av国片精品| 国产精品一区二区精品视频观看| 国产男女超爽视频在线观看| 精品国产乱子伦一区二区三区| 一区二区日韩欧美中文字幕| 三级毛片av免费| 悠悠久久av| 亚洲人成电影观看| 不卡av一区二区三区| 日韩三级视频一区二区三区| 国产亚洲欧美98| 嫁个100分男人电影在线观看| 成年人黄色毛片网站| 免费不卡黄色视频| 国产亚洲精品一区二区www | 日韩欧美免费精品| 免费一级毛片在线播放高清视频 | 久久热在线av| 国产主播在线观看一区二区| 老司机午夜福利在线观看视频| 老熟妇乱子伦视频在线观看| 成人免费观看视频高清| 欧美精品av麻豆av| 国产xxxxx性猛交| 99久久精品国产亚洲精品| 免费在线观看影片大全网站| 久久久久国内视频| 亚洲国产精品合色在线| 丝袜美腿诱惑在线| 看片在线看免费视频| 亚洲av美国av| 欧美国产精品一级二级三级| 免费黄频网站在线观看国产| 亚洲熟女精品中文字幕| 亚洲欧美精品综合一区二区三区| 性色av乱码一区二区三区2| 精品卡一卡二卡四卡免费| 女人久久www免费人成看片| 丁香六月欧美| 天堂√8在线中文| 亚洲精品国产色婷婷电影| 亚洲精品久久成人aⅴ小说| 国产亚洲欧美精品永久| 少妇粗大呻吟视频| 精品免费久久久久久久清纯 | 五月开心婷婷网| 美女福利国产在线| 国产一区二区激情短视频| 免费一级毛片在线播放高清视频 | 一边摸一边抽搐一进一出视频| 久久人妻福利社区极品人妻图片| 亚洲 欧美一区二区三区| 国产一区二区三区综合在线观看| 久久人妻av系列| 91麻豆精品激情在线观看国产 | 中国美女看黄片| 国产乱人伦免费视频| 男女免费视频国产| av电影中文网址| 19禁男女啪啪无遮挡网站| 中文字幕最新亚洲高清| 久久狼人影院| 精品久久久久久电影网| 精品电影一区二区在线| 69精品国产乱码久久久| 999久久久精品免费观看国产| 99riav亚洲国产免费| 变态另类成人亚洲欧美熟女 | 9热在线视频观看99| 免费在线观看影片大全网站| 777久久人妻少妇嫩草av网站| 99精国产麻豆久久婷婷| 成人亚洲精品一区在线观看| av中文乱码字幕在线| 两个人免费观看高清视频| 亚洲第一av免费看| 18禁国产床啪视频网站| 国产一区有黄有色的免费视频| 99国产极品粉嫩在线观看| 久久中文字幕一级| 国产精品自产拍在线观看55亚洲 | 久久久国产欧美日韩av| 国产日韩一区二区三区精品不卡| 免费在线观看日本一区| 色老头精品视频在线观看| 欧美日韩黄片免| 成人特级黄色片久久久久久久| 精品一区二区三区av网在线观看| 久久性视频一级片| 亚洲久久久国产精品| 久久精品熟女亚洲av麻豆精品| 久久中文看片网| 日本黄色日本黄色录像| 国产国语露脸激情在线看| 99在线人妻在线中文字幕 | 亚洲熟女精品中文字幕| 中文欧美无线码| 国产成人精品无人区| 老鸭窝网址在线观看| 国产欧美日韩一区二区精品| 人妻久久中文字幕网| 欧美一级毛片孕妇| 香蕉国产在线看| 亚洲精品久久午夜乱码| 亚洲专区中文字幕在线| 亚洲精品美女久久久久99蜜臀| 国产成人精品无人区| 最新在线观看一区二区三区| 高清av免费在线| 精品第一国产精品| 日韩免费高清中文字幕av| 日韩制服丝袜自拍偷拍| 人成视频在线观看免费观看| 中文字幕人妻熟女乱码| 久久ye,这里只有精品| 香蕉久久夜色| 另类亚洲欧美激情| 亚洲aⅴ乱码一区二区在线播放 | 热99国产精品久久久久久7| 两个人免费观看高清视频| 色精品久久人妻99蜜桃| 18禁黄网站禁片午夜丰满| 不卡一级毛片| 日韩中文字幕欧美一区二区| 动漫黄色视频在线观看| 免费在线观看完整版高清| 免费女性裸体啪啪无遮挡网站| 最近最新免费中文字幕在线| 91老司机精品| 久久人人97超碰香蕉20202| 十八禁网站免费在线| 操美女的视频在线观看| 免费在线观看视频国产中文字幕亚洲| 丁香欧美五月| 国产精品 国内视频| 一区二区三区国产精品乱码| 无遮挡黄片免费观看| 久久精品亚洲av国产电影网| 午夜91福利影院| 亚洲久久久国产精品| 高清av免费在线| av中文乱码字幕在线| 99riav亚洲国产免费| 国产亚洲精品一区二区www | 在线看a的网站| 亚洲 欧美一区二区三区| 另类亚洲欧美激情| 国产精品亚洲av一区麻豆| 少妇 在线观看| 亚洲av熟女| 国产亚洲精品一区二区www | 国产精品美女特级片免费视频播放器 | 免费久久久久久久精品成人欧美视频| 12—13女人毛片做爰片一| 电影成人av| 国产成+人综合+亚洲专区| 99在线人妻在线中文字幕 | 久久人人97超碰香蕉20202| 亚洲av日韩在线播放| 午夜日韩欧美国产| 久久精品aⅴ一区二区三区四区| 亚洲色图 男人天堂 中文字幕| 色老头精品视频在线观看| 色94色欧美一区二区| 黄色a级毛片大全视频| 国产高清激情床上av| 一级作爱视频免费观看| 久久久久国产精品人妻aⅴ院 | 国产区一区二久久| 老汉色∧v一级毛片| 叶爱在线成人免费视频播放| 国产精品电影一区二区三区 | 十分钟在线观看高清视频www| 黑人巨大精品欧美一区二区蜜桃| 高清毛片免费观看视频网站 | www.999成人在线观看| 精品午夜福利视频在线观看一区|