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

    對轉(zhuǎn)槳和單槳空泡水筒噪聲測量對比試驗(yàn)研究

    2018-07-30 09:19:52杜選民劉竹青
    船舶力學(xué) 2018年7期
    關(guān)鍵詞:空泡空化螺旋槳

    曾 賽,杜選民,范 威 ,劉竹青

    (1.上海船舶電子設(shè)備研究所 水聲對抗技術(shù)重點(diǎn)實(shí)驗(yàn)室,上海201108;2.中國船舶科學(xué)研究中心,江蘇 無錫 214082)

    0 引 言

    對轉(zhuǎn)螺旋槳具有推進(jìn)效率高、扭矩平衡等優(yōu)點(diǎn),被應(yīng)用于水下高速航行體、高速快艇等[1]。單槳螺旋槳?jiǎng)t被廣泛應(yīng)用在水面船舶上,如民船領(lǐng)域。近年來,水下目標(biāo)的分類與識別技術(shù)備受各國聲納和艦船科研工作者的關(guān)注,用于水下監(jiān)測系統(tǒng)、水下作戰(zhàn)武器自導(dǎo)系統(tǒng)、水聲引信的水聲探測的目標(biāo)分類識別技術(shù),是研發(fā)水聲對抗器材、水下探測系統(tǒng)、港口預(yù)警防御系統(tǒng)以及水下攻防設(shè)備的關(guān)鍵技術(shù)[2]。隨著人類海洋活動(dòng)的日益頻繁,水下或者水上航行體的聲隱身指標(biāo)越來越重要,作為重要噪聲源的螺旋槳始終是被關(guān)注的重點(diǎn)。已有的船舶噪聲治理文獻(xiàn)表明,艦船的三大水下噪聲源中螺旋槳噪聲對總噪聲貢獻(xiàn)最大[3]。從螺旋槳輻射噪聲中提取與水下以及水面運(yùn)動(dòng)目標(biāo)相關(guān)信號特征是分類水下目標(biāo)的有效途徑。為了探索螺旋槳空泡以及非空泡輻射噪聲的水聲信號特征以及相關(guān)的信號處理方法,比較對轉(zhuǎn)螺旋槳和單槳在輻射噪聲上的差異性,研究螺旋槳噪聲測試的實(shí)驗(yàn)方法以及獲得對轉(zhuǎn)螺旋槳和單槳的輻射噪聲數(shù)據(jù)是十分必要的。目前,與航空領(lǐng)域?qū)β菪龢肼暤难芯肯啾龋瑢W(xué)術(shù)界對于水下螺旋槳空泡以及非空泡噪聲的產(chǎn)生機(jī)理和模型建立尚不成熟和完備,還不能精確預(yù)報(bào)螺旋槳的空泡和非空泡噪聲譜級,加之海試實(shí)驗(yàn)存在設(shè)備和人力成本高昂的問題,且水聲環(huán)境的影響不易控制,故在研究的初期采用符合ITTC(國際拖曳水池會(huì)議International Towing Tank Conference)國際標(biāo)準(zhǔn)的水筒試驗(yàn)方式是經(jīng)濟(jì)可行的,且不斷有科技工作者在大型空泡水筒或者減壓水池等設(shè)備中來探索螺旋槳的空泡和非空泡噪聲特性,試圖發(fā)現(xiàn)螺旋槳的輻射噪聲規(guī)律[4]。

    眾多國內(nèi)外研究人員在空泡水筒或水洞中應(yīng)用縮比模型針對螺旋槳空化噪聲進(jìn)行了試驗(yàn)研究[5]。ITTC亦給出了進(jìn)行水筒噪聲測試的注意要點(diǎn)和進(jìn)行水洞試驗(yàn)的相似準(zhǔn)則,并且逐步成為此后水洞和水筒噪聲測試的準(zhǔn)則[6]。國內(nèi)科研工作者也利用水筒或水洞設(shè)備對螺旋槳噪聲進(jìn)行了探索,突出的工作來自于中國船舶科學(xué)研究中心的朱錫清[7],其利用空泡水筒測試了不同對轉(zhuǎn)槳的噪聲和水動(dòng)力性能,得到了一系列關(guān)于對轉(zhuǎn)槳噪聲特性的研究成果。曾賽[8-10]針對對轉(zhuǎn)螺旋槳的無空化噪聲進(jìn)行了深入的理論分析,進(jìn)行了數(shù)值仿真驗(yàn)證,并通過空泡水筒試驗(yàn)驗(yàn)證了理論分析和數(shù)值計(jì)算結(jié)果。楊勇等[11]利用中型空泡水筒對單槳的空泡噪聲進(jìn)行了測量,提出了利用橫向和垂向加速度計(jì)進(jìn)行噪聲消除的技術(shù)。馮源等[12]利用水洞研究了螺旋槳在加速狀態(tài)下的輻射噪聲特征。胡健等[13]利用水筒測試了單槳的空泡噪聲以此來驗(yàn)證其空泡數(shù)值模型。

    然而,已有的研究工作都是針對單一槳型展開模型試驗(yàn),對轉(zhuǎn)槳和單槳由于結(jié)構(gòu)上的差異,其輻射噪聲譜在結(jié)構(gòu)上也應(yīng)該存在差異,尚沒有公開文獻(xiàn)表明有人對此展開過對比研究。本文討論了在ITTC準(zhǔn)則下空泡水筒測量螺旋槳噪聲的試驗(yàn)條件,提出了相應(yīng)的試驗(yàn)工況設(shè)計(jì)方法,針對水筒測量信號局限性的問題,提出了多頻段調(diào)制譜處理方法的信號處理手段,提取了對轉(zhuǎn)槳和單槳在空泡以及非空泡工況下的低頻調(diào)制線譜,比較了二者的特征,為進(jìn)一步的試驗(yàn)、信號特征和處理方法奠定了基礎(chǔ)。

    1 試驗(yàn)介紹

    1.1 試驗(yàn)水筒簡介

    為了比較研究對轉(zhuǎn)槳和單槳分別在空化和非空化狀態(tài)下輻射噪聲的水聲信號特征,在中國船舶科學(xué)研究中心的大型空泡水筒中進(jìn)行了對轉(zhuǎn)槳和單槳的空化和非空化輻射噪聲測量試驗(yàn)。該大型空泡水筒實(shí)驗(yàn)段的體積為3.2 m×0.8 m×0.8 m,水筒可以進(jìn)行較大范圍內(nèi)的壓力調(diào)節(jié),水流速度范圍為2~20 m/s,流速不均勻度小于1%,最小空泡數(shù)為0.15??张菟补ぷ鞫问疽鈭D如圖1所示。

    空泡水筒實(shí)驗(yàn)段的側(cè)壁安置有體積為1 m×0.6 m×0.6 m的水聲測量消聲艙,為了盡可能減少水流湍動(dòng)對噪聲測量的影響,在水筒實(shí)驗(yàn)段的側(cè)壁上采用有機(jī)玻璃透聲窗將水聲測量消聲艙與實(shí)驗(yàn)段隔離,在數(shù)據(jù)處理階段,可以忽略有機(jī)透聲窗對聲信號傳播的影響。水聲艙體積較為狹小,不能產(chǎn)生開闊水域的聲場測量條件,為了減小壁面反射對噪聲測量的影響,在水聲艙內(nèi)層四壁布滿了消聲尖劈。在消聲艙與實(shí)驗(yàn)段相結(jié)合的部位安裝有B&k8105標(biāo)準(zhǔn)水聽器。采集的數(shù)據(jù)最終在顯示屏上呈現(xiàn)。以上的性能指標(biāo)使得該大型空泡水筒能夠用于進(jìn)行螺旋槳噪聲測試。

    1.2 試驗(yàn)槳模簡介

    本次實(shí)驗(yàn)擬在空泡水筒中測量對轉(zhuǎn)槳和單槳在均勻進(jìn)流條件下的無空化和空化輻射噪聲信號。選用的對轉(zhuǎn)槳槳模參數(shù)如表1所示,前槳為DTMB3686,后槳為DTMB3849。單槳槳模選用DTMB3686槳模。

    圖1 空泡水筒工作段的示意圖Fig.1 Schematic diagram of cavitation tunnel working section

    表1 實(shí)驗(yàn)?zāi)P突緟?shù)Tab.1 The basic parameters of experimental model

    圖2 加工后槳模示意圖Fig.2 Schematic diagram of propeller model

    圖3 槳模型安裝圖Fig.3 Schematic diagram of mounted propeller model

    加工好的槳模如圖2所示。對轉(zhuǎn)槳的前后槳采用不同的動(dòng)力儀進(jìn)行驅(qū)動(dòng),前槳用直流動(dòng)力儀驅(qū)動(dòng),后槳用斜流動(dòng)力儀驅(qū)動(dòng)。模型安裝如圖3所示。

    2 試驗(yàn)方法和工況設(shè)計(jì)

    2.1 試驗(yàn)準(zhǔn)則

    根據(jù)ITTC的建議和中國船舶科學(xué)研究中心的試驗(yàn)規(guī)程,在該大型空泡水筒中進(jìn)行對轉(zhuǎn)槳槳模試驗(yàn)時(shí),前后槳安裝在空泡水筒的轉(zhuǎn)動(dòng)軸上,前槳使用同軸的長軸動(dòng)力儀驅(qū)動(dòng),后槳使用斜流動(dòng)力儀驅(qū)動(dòng),斜流動(dòng)力儀安裝在筒壁上,前后槳的軸線就是旋轉(zhuǎn)軸的中心線,安裝的過程中要確保對轉(zhuǎn)槳的前后槳相對來流水速無攻角和無漂角,前后槳的間距固定,前后槳在軸上通過鍵槽固定牢固,確保在試驗(yàn)中前后槳在來流的沖擊下不出現(xiàn)抖動(dòng)。為了確保水筒中水流的“純凈”,在槳模安裝好以后,需要進(jìn)行一段時(shí)間的除氣處理,因?yàn)槿苡谒械目諝獾刃в跉夂?,氣核的存在將直接影響到對轉(zhuǎn)槳以及單槳的空泡起始,這對于非空化噪聲的測量試驗(yàn)是不利的;若不進(jìn)行除氣處理,一旦產(chǎn)生一定程度的空泡,這部分氣核空泡會(huì)阻止其他部分空泡的崩潰,影響測量的總噪聲級,還會(huì)直接影響到聲學(xué)介質(zhì)的特性,這對于空泡噪聲的測量試驗(yàn)是不利的。試驗(yàn)測量開始之前進(jìn)行排氣處理,經(jīng)過處理后,水筒中的氣泡含量比為0.85。

    為了模擬螺旋槳所在流場的真實(shí)情形,要求螺旋槳0.7R處雷諾數(shù)超過臨界雷諾數(shù)[6],即,

    式中:va為螺旋槳前進(jìn)速度;C0.75R為0.75R處槳葉切面弦長;R為螺旋槳半徑;n為螺旋槳轉(zhuǎn)速;D為螺旋槳直徑;ν為水的運(yùn)動(dòng)粘滯系數(shù)。在試驗(yàn)開始前的參數(shù)設(shè)計(jì)階段要考慮到這些準(zhǔn)則的要求和限制,使設(shè)計(jì)的參數(shù)滿足準(zhǔn)則要求。

    2.2 試驗(yàn)平臺搭建

    在空泡水筒中進(jìn)行螺旋槳噪聲測試試驗(yàn),包含有兩部分系統(tǒng)的搭建,第一部分是動(dòng)力儀部分,在對轉(zhuǎn)槳噪聲測量時(shí),前槳使用直流長軸動(dòng)力儀驅(qū)動(dòng),后槳使用斜流動(dòng)力儀驅(qū)動(dòng),如圖3所示;第二部分為水聲信號的采集與分析部分,采用注水艙單水聽器測量的方法。試驗(yàn)測量平臺的原理如圖4所示。試驗(yàn)中,為了觀察槳葉表面的空化情形以及附近的流場形態(tài),采用頻閃儀輔助觀察。圖5為顯示螺旋槳轉(zhuǎn)速的裝置,結(jié)合該裝置給出的螺旋槳轉(zhuǎn)速便可以調(diào)節(jié)頻閃儀輔助觀察螺旋槳槳葉表面的流場以及空化情形。

    圖4 試驗(yàn)測量平臺示意圖Fig.4 Schematic diagram of the experimental measurement platform

    圖5 轉(zhuǎn)速顯示與頻閃輔助觀察設(shè)備Fig.5 Schematic diagram of speed display and auxiliary stroboscopic observation equipment

    2.3 試驗(yàn)工況設(shè)計(jì)

    根據(jù)1.1給出的試驗(yàn)準(zhǔn)則以及試驗(yàn)條件,對轉(zhuǎn)槳和單槳非空化和空化噪聲測試的工況設(shè)計(jì)步驟為:

    (1)根據(jù)螺旋槳尺寸以及動(dòng)力儀量程初步估算空泡水筒中進(jìn)行螺旋槳試驗(yàn)時(shí)的極限轉(zhuǎn)速,為了避免給動(dòng)力儀帶來損壞,螺旋槳的設(shè)計(jì)工況轉(zhuǎn)速不能高于該極限轉(zhuǎn)速。

    (2)結(jié)合螺旋槳的水動(dòng)力性能曲線和極限轉(zhuǎn)速的限制,計(jì)算起始試驗(yàn)的參數(shù),即水筒來流速度、水筒壓力調(diào)節(jié)值、空泡數(shù)和模型槳起始轉(zhuǎn)速。確保模型槳能夠滿足臨界雷諾數(shù)的要求。

    (3)只有槳轂時(shí)噪聲的測量。在進(jìn)行有槳噪聲測試之前,針對每一種工況下的槳轂噪聲進(jìn)行單獨(dú)測量,方便與有槳葉時(shí)的噪聲進(jìn)行比較研究。

    (4)保持水筒壓力和來流速度不變,調(diào)節(jié)控制平臺的螺旋槳轉(zhuǎn)速控制按鈕,使轉(zhuǎn)速達(dá)到設(shè)定值,待穩(wěn)定后獲得螺旋槳噪聲水筒試驗(yàn)數(shù)據(jù)。

    3 試驗(yàn)結(jié)果與數(shù)據(jù)分析

    3.1 試驗(yàn)工況

    空泡水筒在設(shè)定水速和固定加壓情形下,通過調(diào)節(jié)轉(zhuǎn)速控制按鈕使得螺旋槳處于試驗(yàn)工況。對轉(zhuǎn)槳和單槳的噪聲試驗(yàn)工況如下表所示。表2中,A表示對轉(zhuǎn)槳的工況1,即只有前后槳轂轉(zhuǎn)動(dòng)時(shí)的噪聲,稱之為對轉(zhuǎn)槳背景噪聲工況;B表示對轉(zhuǎn)槳噪聲測試工況,其中B-5為對轉(zhuǎn)槳空化工況;含氣比是指空泡水筒經(jīng)過排氣處理后,水筒的含氣量與原含氣量的比值;進(jìn)速系數(shù)按照定義計(jì)算而來,此時(shí)的進(jìn)速系數(shù)是指前槳的進(jìn)速系數(shù),當(dāng)試驗(yàn)的水流速度和對轉(zhuǎn)槳轉(zhuǎn)速確定后便可以計(jì)算得到進(jìn)速系數(shù);下表3中的項(xiàng)目意義與表2一樣,其中C表示單槳背景噪聲工況,D為單槳噪聲測試工況。D-5為單槳空化工況。

    表2 對轉(zhuǎn)槳噪聲試驗(yàn)工況表Tab.2 The noise test table of counter-rotation propeller

    表3 單槳噪聲試驗(yàn)工況表Tab.3 The noise test table of single propeller

    續(xù)表3

    3.2 對轉(zhuǎn)槳和單槳聲壓頻譜比較

    在空泡水筒中進(jìn)行螺旋槳噪聲測試試驗(yàn),應(yīng)遵守ITTC給出的規(guī)范準(zhǔn)則,即聲模數(shù)準(zhǔn)則。

    式中:N為聲模數(shù),c為水中聲速,f為噪聲頻率,V為實(shí)驗(yàn)段體積。利用本試驗(yàn)使用的空泡水筒參數(shù)可以計(jì)算得到,可以利用的噪聲頻段為800 Hz以上。本試驗(yàn)利用的噪聲測試頻段為800 Hz~63 kHz。

    按照表2和表3給出的工況情形,分別測量離槳盤面0.642 m處的螺旋槳噪聲,其中對轉(zhuǎn)槳測量為前槳的盤面處。對所測得的噪聲時(shí)域信號進(jìn)行1/3倍頻程分析,得到各類工況下的1/3倍頻程功率譜如下圖所示。圖中直接從水筒的有效頻段開始作圖,去掉了800 Hz以下部分,n為螺旋槳轉(zhuǎn)速,rps表示轉(zhuǎn)每秒。為了作圖美觀,橫坐標(biāo)使用對數(shù)坐標(biāo)。

    圖6 對轉(zhuǎn)槳A工況槳轂1/3倍頻程功率譜Fig.6 Schematic diagram of 1/3 octave spectrum of counter-rotation propeller hub

    圖7 單槳C工況槳轂1/3倍頻程功率譜Fig.7 Schematic diagram of 1/3 octave spectrum of single propeller hub

    如圖6-7所示為對轉(zhuǎn)槳和單槳槳轂在試驗(yàn)空泡水筒中1/3倍頻程功率譜曲線。從圖中可以看出,無論是對轉(zhuǎn)槳還是單槳,其槳轂噪聲在整個(gè)測量頻段上均有分布,動(dòng)力儀槳轂噪聲在不同的頻段上分布有較大差異,隨著軸轉(zhuǎn)速的升高,背景噪聲級逐漸增加,在測量的頻段內(nèi),對轉(zhuǎn)槳典型工況下的總背景噪聲級分別為109.1 dB、112.8 dB、115.1 dB、116.2 dB和120.8 dB;單槳典型工況下的總背景噪聲級為100.7 dB、103.8 dB、104.6 dB、105.6 dB和109.3 dB。從中可以看出,由于對轉(zhuǎn)槳后槳槳轂的轉(zhuǎn)動(dòng),使得對轉(zhuǎn)槳的背景噪聲比單槳大9~11.5 dB,驅(qū)動(dòng)后槳轉(zhuǎn)動(dòng)的斜流動(dòng)力儀對整個(gè)背景噪聲貢獻(xiàn)較大。對轉(zhuǎn)槳和單槳的槳轂背景噪聲存在一些特征,首先,在3.8 kHz和8 kHz附近均存在功率譜峰值。其次,對轉(zhuǎn)槳背景噪聲在800 Hz~10 kHz頻段內(nèi)噪聲功率譜較高,10 kHz~30 kHz頻段內(nèi)噪聲功率譜按照倍頻程10 dB衰減,在30~63 kHz頻段內(nèi)為高頻連續(xù)譜,對轉(zhuǎn)槳噪聲衰減非常緩慢;單槳的背景噪聲頻譜特征與對轉(zhuǎn)槳的基本一致。

    圖8 對轉(zhuǎn)槳非空泡工況1/3倍頻程功率譜Fig.8 Schematic diagram of non-cavitation 1/3 octave spectrum of CRP

    圖9單槳非空泡工況1/3倍頻程功率譜Fig.9 Schematic diagram of non-cavitation 1/3 octave spectrum of single propeller

    圖8 所示為對轉(zhuǎn)槳模在非空化工況下的噪聲1/3倍頻程功率譜,在非空化工況情形下,總噪聲級隨著轉(zhuǎn)速的升高而增加,在測試頻段內(nèi)總聲級分別為116.8 dB、119.2 dB、126.2 dB和131.1d B。在800 Hz~5 kHz頻段內(nèi),噪聲功率譜存在兩個(gè)峰值,這與水筒自身特性有關(guān),主要由水筒的混響造成,詳見3.3節(jié)。在5~40 kHz頻段內(nèi),噪聲功率譜呈現(xiàn)穩(wěn)定的衰減特性,在測量工況情形下以每倍頻程9.25~10.5 dB規(guī)律快速衰減,至40 kHz時(shí),噪聲功率譜趨于穩(wěn)定;對比圖8和圖6所示的對轉(zhuǎn)槳背景輻射噪聲頻譜,可以看出,在800 Hz~63 kHz頻段內(nèi),對轉(zhuǎn)槳噪聲級比背景噪聲高5 dB以上,說明測試信號具有較好的信噪比。在1.8 kHz和3.8 kHz附近,有槳時(shí)的譜峰值明顯高于無槳時(shí)的背景噪聲。

    圖9所示為單槳槳模在非空化工況下的噪聲1/3倍頻程功率譜。其噪聲級隨轉(zhuǎn)速而變化的規(guī)律與對轉(zhuǎn)槳一致,在測試頻段內(nèi),典型工況對應(yīng)的總噪聲級分別為110 dB、113.4 dB、118.6 dB和132.2 dB。在800 Hz~5 kHz頻段內(nèi)存在兩個(gè)明顯的峰值,分別在3.8 kHz附近和8 kHz附近,這兩個(gè)峰值產(chǎn)生的原因與對轉(zhuǎn)槳非空化噪聲功率譜曲線中峰值產(chǎn)生的原因是不同的,詳見3.3節(jié)。單槳轉(zhuǎn)速為14 rps時(shí),在8 kHz附近產(chǎn)生非常劇烈的峰值,高出背景噪聲達(dá)40 dB,原因是此時(shí)單槳發(fā)生了“唱音”效應(yīng),在試驗(yàn)測試中也聽到了單槳螺旋槳發(fā)出的“嘯叫”。在10~40 kHz頻段內(nèi),噪聲功率譜呈現(xiàn)穩(wěn)定的衰減特性,在測量工況情形下以每倍頻程11~14 dB規(guī)律快速衰減,至40 kHz時(shí),噪聲功率譜趨于穩(wěn)定;比較圖9和圖7所示的單槳背景輻射噪聲頻譜,可以看出背景輻射噪聲較高,在800 Hz~3 kHz范圍內(nèi),信噪比在1~3 dB,至3 kHz以后,信噪比大于6 dB,說明測試信號具有較好的信噪比。

    比較圖8和圖9對轉(zhuǎn)槳和單槳的非空化噪聲功率譜,可以發(fā)現(xiàn)其在噪聲級、噪聲衰減特性、信噪比和頻譜峰值上存在差異。在相同工況下,對轉(zhuǎn)槳非空化噪聲級明顯高于單槳非空化噪聲級。這是由于后槳的存在以及流場特征差異造成的;在噪聲衰減特性上,二者基本一致,單槳的非空化噪聲每倍頻程的衰減要大于對轉(zhuǎn)槳非空化噪聲;在800 Hz以上的研究頻段,空泡水筒中對轉(zhuǎn)槳非空化噪聲的測試信噪比要優(yōu)于單槳;頻譜峰值上,對轉(zhuǎn)槳的頻譜峰值集中在5 kHz以下,主要與空泡水筒自身的結(jié)構(gòu)特征有關(guān),單槳的頻譜峰值除了與水筒自身結(jié)構(gòu)特征有關(guān),與槳模自身的特征有關(guān)。

    圖10單槳與對轉(zhuǎn)槳空泡工況1/3倍頻程功率譜Fig.10 Schematic diagram of cavitation 1/3 octave spectrum of CRP and single propeller

    圖10 所示為單槳和對轉(zhuǎn)槳空化工況下噪聲測試結(jié)果,從中可以看出,空泡工況下,對轉(zhuǎn)槳和單槳的1/3倍頻程功率譜曲線變化趨勢基本一致,在3.8 kHz附近均存在頻譜峰值,這主要與水筒自身的結(jié)構(gòu)特征有關(guān)。5~40 kHz頻段內(nèi),噪聲倍頻程衰減值比非空泡工況下要小,對轉(zhuǎn)槳和單槳的空泡工況下倍頻程噪聲級按照6 dB衰減。

    3.3 空泡水筒自身特征自振頻率的影響

    根據(jù)對圖6~10的單槳和對轉(zhuǎn)槳噪聲1/3倍頻程功率譜圖分析可知,在1.8 kHz(圖8)、3.8 kHz(圖6~10)和 8 kHz(圖 7、圖 9)附近會(huì)出現(xiàn)峰值譜線。 根據(jù) Ross的《Mechanics of Underwater Noise》[14],這是由空泡水筒的特征頻率所決定的,聲波在有限體積內(nèi)傳播會(huì)產(chǎn)生振蕩現(xiàn)象??张菟驳奶卣髯哉耦l率為:

    其中:nx,ny和nz為空泡水筒邊壁的方向余弦,lx,ly和lz為水筒三個(gè)方向的特征尺度,C為此時(shí)水中聲速。對于試驗(yàn)使用空泡水筒,三個(gè)方向的特征尺度分別為0.4 m、0.4 m和1.6 m,水中聲速取C=1 450 m/s。當(dāng)取水筒x方向的筒壁方向余弦值為0.5,y方向的筒壁方向余弦值為0,z方向的筒壁方向余弦值為0時(shí),f=1.86 kHz;當(dāng)取水筒x方向的筒壁方向余弦值為1,y方向的筒壁方向余弦值為0,z方向的筒壁方向余弦值為0時(shí), f=3.7 kHz;計(jì)算說明,該空泡水筒的特征自振頻率為1.86 kHz和3.7 kHz,試驗(yàn)值在該值附近會(huì)出現(xiàn)異?,F(xiàn)象,事實(shí)上,由圖6~10中,在1.8 kHz和3.8 kHz附近出現(xiàn)了峰值譜線。而單槳測試中在14 rps時(shí)在8 kHz時(shí)出現(xiàn)非常高的峰值,這是由單槳自身的“唱音”所引起。

    3.4 調(diào)制特征線譜分析

    螺旋槳噪聲中有本質(zhì)上存在的一些特征,這些本征特征對于水下以及水面目標(biāo)的識別具有重要意義,線譜特征又稱為“聲紋”特征,是一些穩(wěn)定的離散線譜分量,其不會(huì)隨著背景噪聲的增加而消失。大量的文獻(xiàn)研究表明,螺旋槳的低頻離散線譜對高頻連續(xù)譜存在調(diào)制作用[15-18]。由于水筒測試噪聲的局限性,使得800 Hz以下的低頻噪聲沒有足夠的信噪比,無法直接使用傳統(tǒng)的傅里葉變換分析螺旋槳低頻離散線譜。3.2節(jié)的分析表明,空泡水筒中對轉(zhuǎn)槳噪聲在800 Hz~40 kHz頻段內(nèi)具有良好的信噪比,單槳噪聲在3~40 kHz頻段內(nèi)具有良好的信噪比,可以使用DEMON(Detection of Envelope Modulation On Noise)分析方法處理具有良好信噪比的中高頻段噪聲,解調(diào)出低頻特征線譜。DEMON特征具有明確的物理意義,可以獲得螺旋槳不變的物理特征,如確定螺旋槳的軸頻、葉頻甚至葉數(shù)等。然而,不同頻段內(nèi)的解調(diào)效果是不同的,本節(jié)采用多頻段調(diào)制譜處理有良好信噪比的實(shí)驗(yàn)數(shù)據(jù)。多頻段調(diào)制譜方法流程圖如圖11所示。

    圖11 多頻段調(diào)制譜理論框圖Fig.11 Schematic diagram of multi-band modulation spectrum theory

    對單槳和對轉(zhuǎn)槳的非空泡工況7 rps和14 rps以及空泡工況18 rps測量數(shù)據(jù)進(jìn)行多頻段調(diào)制譜分析,得到的這三個(gè)工況下的調(diào)制譜如下圖所示,其中,選擇的解調(diào)區(qū)間為3~40 kHz,劃分的子帶帶寬為1 kHz,橫坐標(biāo)為頻率,頻段為0~1 kHz,縱坐標(biāo)為歸一化幅度值,幅度的大小表示該線譜調(diào)制的深度。

    圖12 單槳7 r/s時(shí)全頻帶內(nèi)的調(diào)制譜Fig.12 The modulation in whole band spectrum of single propeller in 7r/s

    圖13 對轉(zhuǎn)槳7 r/s時(shí)全頻帶內(nèi)的調(diào)制譜Fig.13 The modulation in whole band spectrum of CR propeller in 7r/s

    圖14 單槳14 r/s時(shí)全頻帶內(nèi)的調(diào)制譜Fig.14 The modulation in whole band spectrum of single propeller in 14 r/s

    圖15 對轉(zhuǎn)槳14 r/s時(shí)全頻帶內(nèi)的調(diào)制譜Fig.15 The modulation in whole band spectrum of CR propeller in 14 r/s

    圖16 單槳18 r/s時(shí)全頻帶內(nèi)的調(diào)制譜Fig.16 The modulation in whole band spectrum of single propeller in 18 r/s

    圖17對轉(zhuǎn)槳18 r/s時(shí)全頻帶內(nèi)的調(diào)制譜Fig.17 The modulation in whole band spectrum of CR propeller in 18 r/s

    圖12 所示為單槳在7 rps工況下的調(diào)制譜。從中可以看出,單槳一階葉頻的調(diào)制作用最強(qiáng),其低頻調(diào)制線譜的頻段主要集中在200 Hz以下,低頻調(diào)制線譜比較單一,且頻段較窄。圖13為對轉(zhuǎn)槳在7 rps時(shí)的調(diào)制譜,從中可以看出,對轉(zhuǎn)槳在低頻段的調(diào)制線譜要比單槳的低頻調(diào)制線譜豐富得多,對轉(zhuǎn)槳前后槳一階葉頻的組合調(diào)制作用最強(qiáng),對轉(zhuǎn)槳調(diào)制線譜的頻段分布范圍較廣,在0~1 kHz的頻段內(nèi)有豐富的調(diào)制線譜出現(xiàn)。圖14與圖15分別為單槳和對轉(zhuǎn)槳在14 rps工況時(shí)的調(diào)制譜,從中可以看出,其對比規(guī)律與圖12和圖13的對比規(guī)律一致。

    圖16所示為單槳在18 rps空泡工況下的調(diào)制譜。從中可以看出,在空化情形下,單槳一階葉頻的調(diào)制作用最強(qiáng),其調(diào)制線譜的分布頻段仍然較窄且較為單一。圖17為對轉(zhuǎn)槳在18 rps時(shí)的調(diào)制譜,從中可以看出,在空泡工況下,對轉(zhuǎn)槳前后槳一階葉頻的組合的調(diào)制作用最強(qiáng),且其調(diào)制譜分布頻段比單槳的調(diào)制譜分布頻段廣,調(diào)制線譜的數(shù)目也較為豐富。

    3.5 單槳與對轉(zhuǎn)槳空泡發(fā)展情形

    設(shè)計(jì)工況B-4為對轉(zhuǎn)槳發(fā)生初生空化的工況,B-5為對轉(zhuǎn)槳發(fā)生梢渦空化的工況;D-4為單槳發(fā)生初生空化的工況,D-5為單槳發(fā)生梢渦空化時(shí)的工況。借助于高分辨率攝像機(jī),在高速頻閃儀的輔助下,可以觀察到發(fā)生在對轉(zhuǎn)槳和單槳葉面上的初生空化和梢渦空化,如圖18和圖20所示。此時(shí)對轉(zhuǎn)槳的槳葉隨邊處面上出現(xiàn)弱空化,但是空化面積非常小,處于初生空化狀態(tài),隨著轉(zhuǎn)速的增加,對轉(zhuǎn)槳和單槳的梢渦渦核壓力進(jìn)一步降低,空化進(jìn)一步生長,如圖19和圖21所示。由于前后槳的相互干擾,使得對轉(zhuǎn)槳的渦線分布發(fā)生紊亂,如圖19所示。在梢渦工況下,渦尾跡出現(xiàn)明顯的收縮現(xiàn)象。

    圖18 對轉(zhuǎn)槳初生空化Fig.18 The initial cavitation of CR propeller

    圖19 對轉(zhuǎn)槳梢渦空化Fig.19 The tip vortex cavitation of CR propeller

    圖20 單槳初生空化Fig.20 The initial cavitation of single propeller

    圖21 單槳梢渦空化Fig.21 The tip vortex cavitation of single propeller

    4 結(jié) 論

    在大型空泡水筒中測量了單槳和對轉(zhuǎn)槳的非空化以及空化噪聲,通過比較二者噪聲特性的差異性,以期為水下目標(biāo)識別和分類提供判斷依據(jù)。且針對測試使用空泡水筒的特性進(jìn)行了探討。得到如下結(jié)論:

    (1)空泡水筒中的背景噪聲較高,其中驅(qū)動(dòng)后槳的斜流動(dòng)力儀背景噪聲最強(qiáng)。

    (2)空泡水筒由于自身結(jié)構(gòu)上的特殊性,在某些頻率附近易出現(xiàn)噪聲異常。

    (3)對轉(zhuǎn)槳和單槳在非空泡噪聲上存在較大差異,表現(xiàn)在噪聲級、噪聲衰減特性、信噪比和三分之一倍頻程譜峰值方面。在相同工況下,對轉(zhuǎn)槳非空泡噪聲級要比單槳高6 dB以上;對轉(zhuǎn)槳非空泡噪聲的倍頻程衰減要略微小于單槳非空泡噪聲的衰減;在該空泡水筒中,800 Hz以上的測試頻段,對轉(zhuǎn)槳噪聲的測試信噪比要優(yōu)于單槳;對轉(zhuǎn)槳非空化噪聲頻譜在1.8 kHz和3.8 kHz附近出現(xiàn)了頻譜峰值,這主要與空泡水筒的混響有關(guān),而單槳非空化噪聲頻譜在3.8 kHz和8 kHz附近出現(xiàn)頻譜峰值,其主要原因分別是混響和單槳“唱音”作用的結(jié)果。

    (4)對轉(zhuǎn)槳和單槳在空泡噪聲上差異較小,其三分之一倍頻程譜結(jié)構(gòu)變化趨勢基本一致。相同工況下,對轉(zhuǎn)槳空泡噪聲要高于單槳空泡噪聲。

    (5)對轉(zhuǎn)槳和單槳的非空化以及空化多頻段調(diào)制譜在結(jié)構(gòu)上存在較大差異,單槳的多頻段調(diào)制譜較為簡單,其一階葉頻的調(diào)制深度最大,且調(diào)制頻段較窄,一般在0~200 Hz。對轉(zhuǎn)槳的多頻段調(diào)制譜較為復(fù)雜,其前后槳一階葉頻的組合調(diào)制深度最大,調(diào)制頻段較寬。

    參 考 文 獻(xiàn):

    [1]常 煜,洪方文,張志榮,等.對轉(zhuǎn)槳水動(dòng)力性能的數(shù)值分析[C]//2008年船舶水動(dòng)力學(xué)學(xué)術(shù)會(huì)議暨中國船舶學(xué)術(shù)界進(jìn)入ITTC30周年紀(jì)念會(huì)論文集.杭州:中國造船工程學(xué)會(huì),2008.Chang Yu,Hong Fanwen,Zhang Zhirong,et al.Numerical study of contra-rotating propellers hydrodynamic performance[C]//The Ship Hydrodynamics Conference and The Shipbuilding Academic of China Step into The ITTC30 Anniversary Conference.Hangzhou:China Shipbuliding Engineering Society,2008.

    [2]Johnson E.Unmanned under sea vehicles and guided missile submarines:Technological and operational synergies[D].Occasional Paper No.27,2002.

    [3]張永坤.艦船螺旋槳輻射噪聲預(yù)報(bào)研究[D].武漢:海軍工程大學(xué),2009.Zhang Yongkun.Investigation on predication ship propeller radiated noise[D].Wuhan:Naval University of Engineering,2009.

    [4]Atlar M,Takinaci A C,Korkut E,et al.Cavitation tunnel tests for propeller noise of a FRV and comparisions with fullscale measurements[C]//Proc of the 4th International Symposium on Cavitation.Pasadena,2001:1-14.

    [5]Park C,Seol H,Kim K,et al.A study on propeller noise source localization in a cavitation tunnel[J].Ocean Engineering,2009,36:754-762.

    [6]ITTC.Cavitation committee report[C]//Proc of 18th ITTC.Kobe,1987:183-186.

    [7]朱錫清,吳武生.水下高速航行體對轉(zhuǎn)螺旋槳線譜噪聲預(yù)報(bào)研究[J].聲學(xué)學(xué)報(bào),1998(2):124-132.Zhu Xiqing,Wu Wusheng.Prediction of line-spectrum noise induced by high speed vehicle counter-rotation propellers in water[J].Acta Acoustica,1998(2):124-132.(in Chinese)

    [8]曾 賽,杜選民,范 威.水下對轉(zhuǎn)槳非空化線譜噪聲分析與數(shù)值研究[J].兵工學(xué)報(bào),2015,36(6):1052-1060.Zeng Sai,Du Xuanmin,Fan Wei.Numerical simulation and analysis of non-cavitation noise line-spectrum frequency of underwater counter-rotation propeller[J].Acta Armamentarii,2015,36(6):1051-1060.(in Chinese)

    [9]曾 賽,杜選民,范 威.水下對轉(zhuǎn)槳無空化噪聲調(diào)制理論分析與試驗(yàn)研究[J].聲學(xué)學(xué)報(bào),2017,42(6):641-650.Zeng Sai,Du Xuanmin,Fan Wei.Theoretically analysis and experimental research on non-cavitation noise modulation mechanism of underwater counter-rotation propeller[J].Acta Acoustica,2017,42(6):641-650.(in Chinese)

    [10]Zeng Sai,Du Xuanmin.A semi-empirical modulation model on non-cavitation noise of underwater counter-rotation propeller[C]//2nd International Conference on Computational Modeling,Simulation and Applied Mathematics,October 22-23,2017.Beijing,China,2017.

    [11]楊 勇,熊 鷹,時(shí)立攀.螺旋槳空泡噪聲測量與分析[J].中國艦船研究,2013(1):85-90.Yang Yong,Xiong Ying,Shi Lipan.Measurement and analysis of the cavitation noise of propellers[J].Chinese Journal of Ship Research,2013(1):85-90.(in Chinese)

    [12]馮 源,陶 然,卓智海.加速旋轉(zhuǎn)螺旋槳噪聲測試水筒試驗(yàn)研究[J].兵工學(xué)報(bào),2010(12):1611-1616.Feng Yuan,Tao Ran,Zhuo Zhihai.Study on cavitation tunnel experiment of accelerating rotating propeller noise[J].Acta Armamentarii,2010(12):1611-1616.(in Chinese)

    [13]胡 健.螺旋槳空泡性能及低噪聲螺旋槳設(shè)計(jì)研究[D].哈爾濱:哈爾濱工程大學(xué),2006.Hu Jian.Research on propeller cavitation characteristics and low noise propeller design[D].Harbin:Harbin Engineering University,2006.

    [14]Ross D.Mechanics of underwater noise[M].Pergaman Press,New York,1976.

    [15]汪德昭.水聲學(xué)[M].第二版.北京:科學(xué)出版社,2009.

    [16]熊紫英,朱錫清.基于LOFAR譜和DEMON譜特征的艦船輻射噪聲研究[J].船舶力學(xué),2007,11(2):301-308.Xiong Ziying,Zhu Xiqing.Ship radiated noise research based on the LOFAR spectrum and DEMON spectrum characteristics[J].Journal of Ship Mechanics,2007,11(2):301-308.(in Chinese)

    [17]陶篤純.艦船噪聲節(jié)奏的研究(Ⅰ)-數(shù)學(xué)模型及功率譜密度[J].聲學(xué)學(xué)報(bào),1983,8(2):66-72.TAO Duchun.A study on ship radiated noise rhythms(I)-Mathematical model and power spectrum densidy[J].Acta A-coustica,1983,8(2):66-74.

    [18]陶篤純.艦船噪聲節(jié)奏的研究(Ⅱ)-自相關(guān)函數(shù)及節(jié)奏信息的提取[J].聲學(xué)學(xué)報(bào),1983,8(5):281-288.TAO Duchun.A study on ship radiated noise rhythms(II)-Autocorrelation function and extraction of rhythm information[J].Acta Acoustica,1983,8(5):280-289.(in Chinese)

    猜你喜歡
    空泡空化螺旋槳
    功率超聲作用下鋼液中空化泡尺寸的演變特性
    鋼鐵釩鈦(2023年5期)2023-11-17 08:48:34
    水下航行體雙空泡相互作用數(shù)值模擬研究
    基于CFD的螺旋槳拉力確定方法
    三維扭曲水翼空化現(xiàn)象CFD模擬
    不同運(yùn)動(dòng)形式下水物相互作用空化數(shù)值模擬
    基于LPV的超空泡航行體H∞抗飽和控制
    基于CFD的對轉(zhuǎn)槳無空泡噪聲的仿真預(yù)報(bào)
    船海工程(2015年4期)2016-01-05 15:53:28
    3800DWT加油船螺旋槳諧鳴分析及消除方法
    廣東造船(2015年6期)2015-02-27 10:52:46
    螺旋槳轂帽鰭節(jié)能性能的數(shù)值模擬
    SPH在水下高速物體空泡發(fā)展模擬中的應(yīng)用
    亚洲无线观看免费| 在线视频色国产色| 天堂动漫精品| 免费搜索国产男女视频| 久久草成人影院| 人妻丰满熟妇av一区二区三区| 国产精品99久久久久久久久| 日本与韩国留学比较| 麻豆国产97在线/欧美| 日本 欧美在线| 国产麻豆成人av免费视频| 男女那种视频在线观看| 国产精品久久久久久久电影 | 日本与韩国留学比较| 九色成人免费人妻av| 国产精品 国内视频| 成年免费大片在线观看| 国产熟女xx| 18禁观看日本| 亚洲精品粉嫩美女一区| 精品一区二区三区视频在线 | 久久国产精品影院| 色噜噜av男人的天堂激情| 国产精品av久久久久免费| 国产精华一区二区三区| 日本一本二区三区精品| 亚洲av五月六月丁香网| 国产不卡一卡二| 我要搜黄色片| 国产v大片淫在线免费观看| 国内精品久久久久精免费| 中文字幕高清在线视频| 国产精品一区二区三区四区免费观看 | 校园春色视频在线观看| 在线十欧美十亚洲十日本专区| av中文乱码字幕在线| 久久久久久久久久黄片| 国产亚洲av嫩草精品影院| 熟女少妇亚洲综合色aaa.| 在线观看舔阴道视频| ponron亚洲| 波多野结衣高清作品| 午夜免费激情av| 在线播放国产精品三级| 变态另类成人亚洲欧美熟女| 国产精品一及| 两性夫妻黄色片| av黄色大香蕉| 欧美日韩国产亚洲二区| 国产黄片美女视频| 大型黄色视频在线免费观看| 国产一区在线观看成人免费| 亚洲精品美女久久av网站| 麻豆av在线久日| 国产精品 欧美亚洲| 亚洲专区字幕在线| 美女免费视频网站| 色播亚洲综合网| 成年女人毛片免费观看观看9| 韩国av一区二区三区四区| 我的老师免费观看完整版| 亚洲av免费在线观看| 欧美xxxx黑人xx丫x性爽| 欧美一级毛片孕妇| 国产乱人伦免费视频| 国产不卡一卡二| 亚洲av熟女| 欧美一级毛片孕妇| 黄色成人免费大全| 日韩欧美一区二区三区在线观看| 男人的好看免费观看在线视频| www.熟女人妻精品国产| 精品一区二区三区视频在线观看免费| 十八禁网站免费在线| 美女cb高潮喷水在线观看 | 一二三四社区在线视频社区8| 18禁裸乳无遮挡免费网站照片| 免费看日本二区| 国产精品电影一区二区三区| 黄色 视频免费看| 波多野结衣高清作品| 香蕉国产在线看| 淫妇啪啪啪对白视频| 999久久久精品免费观看国产| 岛国视频午夜一区免费看| 久久性视频一级片| 又黄又粗又硬又大视频| 人人妻人人澡欧美一区二区| 国产一区二区三区在线臀色熟女| 国产精品久久久久久人妻精品电影| 亚洲色图 男人天堂 中文字幕| 国产av一区在线观看免费| av福利片在线观看| 91麻豆精品激情在线观看国产| 久久久久国内视频| 欧美一级a爱片免费观看看| 日本黄色视频三级网站网址| 亚洲精品粉嫩美女一区| 国产69精品久久久久777片 | 久久精品91无色码中文字幕| 久久人妻av系列| 欧美三级亚洲精品| 深夜精品福利| 精品无人区乱码1区二区| 亚洲av成人一区二区三| 美女大奶头视频| 亚洲欧美一区二区三区黑人| 高清在线国产一区| 亚洲中文字幕一区二区三区有码在线看 | 国产单亲对白刺激| 激情在线观看视频在线高清| 男女做爰动态图高潮gif福利片| 天天躁日日操中文字幕| 男女那种视频在线观看| 婷婷精品国产亚洲av在线| 免费看美女性在线毛片视频| 亚洲国产精品久久男人天堂| av在线蜜桃| 久久香蕉国产精品| 狂野欧美激情性xxxx| 天天一区二区日本电影三级| 久久精品国产清高在天天线| 日韩有码中文字幕| 亚洲最大成人中文| 久久久久久久久久黄片| 亚洲欧美激情综合另类| 亚洲自偷自拍图片 自拍| 美女午夜性视频免费| 亚洲黑人精品在线| 日韩欧美国产在线观看| 色哟哟哟哟哟哟| 亚洲精品国产精品久久久不卡| 一进一出抽搐动态| 在线十欧美十亚洲十日本专区| 久久精品aⅴ一区二区三区四区| 成人无遮挡网站| 国产aⅴ精品一区二区三区波| 动漫黄色视频在线观看| 亚洲av电影在线进入| 亚洲中文日韩欧美视频| 男人舔女人的私密视频| 精品国产乱码久久久久久男人| 老司机福利观看| 小说图片视频综合网站| 男女床上黄色一级片免费看| 国产亚洲欧美98| 成人高潮视频无遮挡免费网站| 九九久久精品国产亚洲av麻豆 | 特级一级黄色大片| 在线免费观看的www视频| 精品熟女少妇八av免费久了| 一个人看视频在线观看www免费 | 丝袜人妻中文字幕| 免费观看精品视频网站| 亚洲自拍偷在线| 色播亚洲综合网| 九色国产91popny在线| 日韩欧美 国产精品| 亚洲熟妇熟女久久| av片东京热男人的天堂| 在线a可以看的网站| 欧美黑人欧美精品刺激| 欧美黄色淫秽网站| 久久精品91蜜桃| 亚洲国产欧洲综合997久久,| 首页视频小说图片口味搜索| 久久香蕉精品热| 亚洲精品美女久久久久99蜜臀| 国产麻豆成人av免费视频| 亚洲自偷自拍图片 自拍| 久久亚洲精品不卡| 99精品久久久久人妻精品| 欧美最黄视频在线播放免费| 亚洲乱码一区二区免费版| 啪啪无遮挡十八禁网站| 国产伦精品一区二区三区四那| 2021天堂中文幕一二区在线观| 69av精品久久久久久| 真实男女啪啪啪动态图| 午夜福利18| 欧美黑人巨大hd| 免费在线观看日本一区| 午夜福利欧美成人| 操出白浆在线播放| 久久久久国产精品人妻aⅴ院| 1024香蕉在线观看| 欧美黄色淫秽网站| 国内精品美女久久久久久| 久久精品aⅴ一区二区三区四区| 欧美性猛交黑人性爽| 性色avwww在线观看| 成年人黄色毛片网站| 制服人妻中文乱码| 最好的美女福利视频网| 两性夫妻黄色片| 国产熟女xx| 99久久精品热视频| 亚洲欧美激情综合另类| 国产黄片美女视频| 亚洲av成人av| tocl精华| 成人一区二区视频在线观看| 午夜福利在线观看免费完整高清在 | 在线观看66精品国产| 九九久久精品国产亚洲av麻豆 | 熟妇人妻久久中文字幕3abv| 长腿黑丝高跟| www.精华液| 两个人的视频大全免费| 免费人成视频x8x8入口观看| 亚洲av中文字字幕乱码综合| 精华霜和精华液先用哪个| 美女cb高潮喷水在线观看 | 免费搜索国产男女视频| 亚洲av成人一区二区三| 亚洲,欧美精品.| 国产成人av教育| 老司机在亚洲福利影院| 国产成人福利小说| 午夜成年电影在线免费观看| 久久久久免费精品人妻一区二区| 一进一出抽搐gif免费好疼| 久久久精品欧美日韩精品| 蜜桃久久精品国产亚洲av| 一进一出抽搐gif免费好疼| 91久久精品国产一区二区成人 | 日韩欧美在线乱码| 欧美黄色片欧美黄色片| 精品午夜福利视频在线观看一区| 久久久久久九九精品二区国产| 性色av乱码一区二区三区2| 国产又黄又爽又无遮挡在线| 亚洲美女黄片视频| 精品熟女少妇八av免费久了| 91在线观看av| 亚洲精品在线观看二区| 狠狠狠狠99中文字幕| 亚洲欧洲精品一区二区精品久久久| 午夜福利在线在线| 午夜精品久久久久久毛片777| 手机成人av网站| 亚洲中文av在线| 午夜成年电影在线免费观看| 精品久久久久久,| 欧美成狂野欧美在线观看| 特级一级黄色大片| 午夜福利免费观看在线| 真实男女啪啪啪动态图| 久久亚洲精品不卡| 少妇的丰满在线观看| 精品国产超薄肉色丝袜足j| 国内毛片毛片毛片毛片毛片| 国产又黄又爽又无遮挡在线| 午夜a级毛片| 热99re8久久精品国产| 在线a可以看的网站| av在线天堂中文字幕| 最近最新中文字幕大全电影3| 亚洲乱码一区二区免费版| 久久九九热精品免费| 国产综合懂色| 国产视频一区二区在线看| av欧美777| 国产精品久久久久久人妻精品电影| 日韩欧美国产在线观看| 亚洲中文字幕日韩| 免费高清视频大片| 级片在线观看| 男人舔女人下体高潮全视频| 女人被狂操c到高潮| 身体一侧抽搐| 夜夜爽天天搞| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av电影在线进入| 真实男女啪啪啪动态图| 久久精品国产亚洲av香蕉五月| 久久久久性生活片| 国产精品综合久久久久久久免费| 亚洲五月婷婷丁香| av天堂中文字幕网| 最近最新中文字幕大全免费视频| 日韩欧美在线乱码| 国产一区二区在线av高清观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲成人精品中文字幕电影| 中文字幕精品亚洲无线码一区| 国产亚洲av高清不卡| 91av网站免费观看| 一区二区三区国产精品乱码| 久久精品影院6| 日韩欧美在线二视频| 神马国产精品三级电影在线观看| 欧美乱妇无乱码| 免费无遮挡裸体视频| a级毛片a级免费在线| 中文字幕最新亚洲高清| 久久精品夜夜夜夜夜久久蜜豆| 欧美色欧美亚洲另类二区| 亚洲精品在线美女| 欧美成人一区二区免费高清观看 | 最近最新免费中文字幕在线| 欧美激情久久久久久爽电影| 搡老妇女老女人老熟妇| 香蕉久久夜色| 天天添夜夜摸| 久久久久久久久中文| 精品无人区乱码1区二区| 亚洲精品在线美女| 高清毛片免费观看视频网站| 国产精品99久久99久久久不卡| 免费看美女性在线毛片视频| av黄色大香蕉| 99久久精品热视频| 法律面前人人平等表现在哪些方面| 好男人电影高清在线观看| 这个男人来自地球电影免费观看| 国产精品98久久久久久宅男小说| 亚洲av熟女| 搡老熟女国产l中国老女人| 99在线人妻在线中文字幕| av黄色大香蕉| 丰满的人妻完整版| 老司机午夜十八禁免费视频| ponron亚洲| 国产精华一区二区三区| 熟女电影av网| 国产精品久久视频播放| 国产乱人伦免费视频| 精品国产乱码久久久久久男人| 麻豆一二三区av精品| 宅男免费午夜| 色尼玛亚洲综合影院| av国产免费在线观看| 在线免费观看不下载黄p国产 | 成人性生交大片免费视频hd| 999久久久精品免费观看国产| 精品久久久久久久末码| 国产毛片a区久久久久| 免费看美女性在线毛片视频| 人人妻人人看人人澡| 久久精品影院6| 级片在线观看| 国产精品一及| 久久人妻av系列| 精品国产美女av久久久久小说| 精品日产1卡2卡| 精品久久久久久,| 老司机午夜十八禁免费视频| 欧美不卡视频在线免费观看| 国产精品免费一区二区三区在线| 欧美日韩一级在线毛片| 在线观看66精品国产| 亚洲欧美日韩卡通动漫| 国产精品国产高清国产av| 日韩欧美免费精品| 国产亚洲精品综合一区在线观看| 嫩草影院精品99| 国产激情偷乱视频一区二区| 床上黄色一级片| 91老司机精品| 十八禁网站免费在线| xxx96com| 国产成年人精品一区二区| av黄色大香蕉| 丰满人妻一区二区三区视频av | 久久精品国产99精品国产亚洲性色| 成年女人看的毛片在线观看| 色综合站精品国产| 蜜桃久久精品国产亚洲av| 国产真人三级小视频在线观看| 国产精品亚洲av一区麻豆| 午夜影院日韩av| 国产97色在线日韩免费| 无遮挡黄片免费观看| 欧美丝袜亚洲另类 | 免费在线观看亚洲国产| 久久人人精品亚洲av| 欧美色欧美亚洲另类二区| 婷婷精品国产亚洲av在线| 国产日本99.免费观看| 亚洲av成人不卡在线观看播放网| 亚洲欧美激情综合另类| 日韩欧美精品v在线| 少妇的丰满在线观看| 搡老岳熟女国产| 国产v大片淫在线免费观看| 18禁黄网站禁片免费观看直播| 老司机午夜福利在线观看视频| 亚洲成人中文字幕在线播放| 久久这里只有精品19| 女同久久另类99精品国产91| 免费看美女性在线毛片视频| 国产97色在线日韩免费| 国产精品国产高清国产av| 国产av不卡久久| 久久久久国产一级毛片高清牌| 欧美日韩瑟瑟在线播放| 亚洲av美国av| 老司机福利观看| 两个人的视频大全免费| 久久久久国内视频| 国产精品99久久久久久久久| 国产黄a三级三级三级人| 精品久久久久久成人av| 此物有八面人人有两片| 国产高清有码在线观看视频| 黑人欧美特级aaaaaa片| 无遮挡黄片免费观看| 精品熟女少妇八av免费久了| 香蕉久久夜色| 亚洲熟妇中文字幕五十中出| 夜夜夜夜夜久久久久| 搞女人的毛片| 国产91精品成人一区二区三区| 免费在线观看成人毛片| 亚洲精品一卡2卡三卡4卡5卡| 日韩三级视频一区二区三区| 成人无遮挡网站| 精品久久蜜臀av无| 久久久久国产一级毛片高清牌| 三级男女做爰猛烈吃奶摸视频| 一个人免费在线观看的高清视频| 99久久成人亚洲精品观看| 久久久久精品国产欧美久久久| av中文乱码字幕在线| 天堂√8在线中文| 欧洲精品卡2卡3卡4卡5卡区| 我的老师免费观看完整版| 琪琪午夜伦伦电影理论片6080| 欧美丝袜亚洲另类 | 国产又黄又爽又无遮挡在线| 国产又黄又爽又无遮挡在线| 人妻夜夜爽99麻豆av| 国产主播在线观看一区二区| 亚洲精品456在线播放app | 久久国产精品人妻蜜桃| 国产蜜桃级精品一区二区三区| 欧美成狂野欧美在线观看| 亚洲成人免费电影在线观看| 国产三级黄色录像| 校园春色视频在线观看| 97人妻精品一区二区三区麻豆| 性色avwww在线观看| 搡老熟女国产l中国老女人| 中文字幕人成人乱码亚洲影| 日本黄色片子视频| 国产精品亚洲一级av第二区| 久久中文字幕人妻熟女| 黑人操中国人逼视频| 国产一区二区在线观看日韩 | 亚洲五月婷婷丁香| 国产精品久久久久久精品电影| 欧美成人性av电影在线观看| 欧美另类亚洲清纯唯美| 亚洲一区二区三区色噜噜| 亚洲精品国产精品久久久不卡| 后天国语完整版免费观看| 午夜精品一区二区三区免费看| 一二三四在线观看免费中文在| 叶爱在线成人免费视频播放| 91老司机精品| 1000部很黄的大片| 午夜福利在线观看吧| 亚洲天堂国产精品一区在线| av天堂中文字幕网| 少妇裸体淫交视频免费看高清| 亚洲国产精品久久男人天堂| av福利片在线观看| 亚洲18禁久久av| 狂野欧美白嫩少妇大欣赏| 欧美日韩黄片免| 精品久久久久久久毛片微露脸| 性欧美人与动物交配| 亚洲国产精品999在线| 老司机午夜十八禁免费视频| 久久亚洲精品不卡| 女人被狂操c到高潮| 国产免费av片在线观看野外av| 男女那种视频在线观看| 99热只有精品国产| 看黄色毛片网站| 亚洲精品粉嫩美女一区| 日本熟妇午夜| 观看美女的网站| 黄色日韩在线| 亚洲欧美日韩卡通动漫| 欧美日韩黄片免| 国模一区二区三区四区视频 | 亚洲aⅴ乱码一区二区在线播放| 韩国av一区二区三区四区| 夜夜爽天天搞| 亚洲国产精品合色在线| 欧美性猛交╳xxx乱大交人| h日本视频在线播放| 一个人免费在线观看电影 | 又黄又粗又硬又大视频| 国产精品女同一区二区软件 | 精品日产1卡2卡| 神马国产精品三级电影在线观看| 99精品久久久久人妻精品| 亚洲国产精品sss在线观看| x7x7x7水蜜桃| 在线永久观看黄色视频| 国产成人啪精品午夜网站| АⅤ资源中文在线天堂| 亚洲在线自拍视频| 免费看美女性在线毛片视频| 国产伦在线观看视频一区| 中亚洲国语对白在线视频| 不卡av一区二区三区| 亚洲av电影在线进入| 麻豆一二三区av精品| 欧美黄色片欧美黄色片| 欧美日韩国产亚洲二区| 天天添夜夜摸| 国产视频一区二区在线看| 免费人成视频x8x8入口观看| av在线蜜桃| 久久草成人影院| 老熟妇仑乱视频hdxx| 精品日产1卡2卡| 美女高潮喷水抽搐中文字幕| 看免费av毛片| 视频区欧美日本亚洲| 精品福利观看| 噜噜噜噜噜久久久久久91| 黄色成人免费大全| 亚洲七黄色美女视频| 日韩国内少妇激情av| 亚洲成av人片免费观看| 国产精品亚洲一级av第二区| 99久久精品国产亚洲精品| 网址你懂的国产日韩在线| 麻豆成人午夜福利视频| 亚洲自拍偷在线| 91麻豆av在线| 国产成人福利小说| 欧美最黄视频在线播放免费| 国产高清有码在线观看视频| 人妻夜夜爽99麻豆av| 精品一区二区三区视频在线 | av在线蜜桃| 午夜免费激情av| 人人妻,人人澡人人爽秒播| 国产精品久久久久久人妻精品电影| 日韩欧美 国产精品| 亚洲成人久久爱视频| 热99在线观看视频| 亚洲男人的天堂狠狠| 久久久久久久午夜电影| 99久久精品一区二区三区| 母亲3免费完整高清在线观看| 天天一区二区日本电影三级| 精品一区二区三区视频在线 | 日韩精品青青久久久久久| 熟女少妇亚洲综合色aaa.| 久久精品aⅴ一区二区三区四区| 99国产精品99久久久久| 久久久精品欧美日韩精品| 久久香蕉精品热| 首页视频小说图片口味搜索| 久久精品国产清高在天天线| 国产人伦9x9x在线观看| 成在线人永久免费视频| 亚洲人成网站在线播放欧美日韩| 国产亚洲av嫩草精品影院| 丁香欧美五月| 欧美中文日本在线观看视频| 国产欧美日韩精品亚洲av| 免费在线观看日本一区| 岛国在线免费视频观看| 国产精品一区二区精品视频观看| bbb黄色大片| av天堂在线播放| 欧美大码av| 亚洲精品美女久久久久99蜜臀| 日韩精品中文字幕看吧| 18美女黄网站色大片免费观看| 免费看a级黄色片| 很黄的视频免费| 91麻豆精品激情在线观看国产| 亚洲在线自拍视频| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区二区三区四区免费观看 | 国产视频一区二区在线看| 亚洲成av人片免费观看| 一级黄色大片毛片| 给我免费播放毛片高清在线观看| 国产毛片a区久久久久| 国产成人av教育| 色在线成人网| 啦啦啦韩国在线观看视频| 天堂动漫精品| 国产精品一区二区三区四区免费观看 | 国内精品一区二区在线观看| 成人av在线播放网站| 久久久久久九九精品二区国产| 高清在线国产一区| 九九久久精品国产亚洲av麻豆 | 欧美在线黄色| e午夜精品久久久久久久| 美女高潮喷水抽搐中文字幕| 久久国产精品人妻蜜桃| 色综合站精品国产| 亚洲国产高清在线一区二区三| 亚洲一区高清亚洲精品| 91av网一区二区| 美女被艹到高潮喷水动态| 最近最新免费中文字幕在线| 午夜日韩欧美国产| 两性夫妻黄色片| 久久亚洲精品不卡| 久久久久久久久免费视频了| 国产高清有码在线观看视频| 婷婷精品国产亚洲av|