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

    重型貨車車內(nèi)聲壓級(jí)預(yù)測(cè)*

    2016-04-11 08:26:58王彬星連小珉鄭四發(fā)任焱晞
    汽車工程 2016年2期
    關(guān)鍵詞:聲壓級(jí)駕駛室損耗

    王彬星,連小珉,鄭四發(fā),任焱晞,劉 玉

    (1.總裝工程兵科研二所,北京 100093; 2.清華大學(xué),汽車安全與節(jié)能國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100084)

    2016038

    重型貨車車內(nèi)聲壓級(jí)預(yù)測(cè)*

    王彬星1,連小珉2,鄭四發(fā)2,任焱晞1,劉 玉1

    (1.總裝工程兵科研二所,北京 100093; 2.清華大學(xué),汽車安全與節(jié)能國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100084)

    本文中對(duì)一重型貨車在20-1 000Hz的頻率范圍內(nèi)的車內(nèi)聲壓級(jí)進(jìn)行了預(yù)測(cè)。首先根據(jù)各種工況下實(shí)測(cè)車內(nèi)噪聲信號(hào)的能量分布情況,確定了分析的頻率范圍。然后基于經(jīng)過(guò)驗(yàn)證的駕駛室有限元模型,綜合利用有限元和邊界元方法,建立了適用于車內(nèi)低頻聲壓級(jí)預(yù)測(cè)的有限元-邊界元模型;綜合利用有限元方法和統(tǒng)計(jì)能量分析方法,建立了適用于車內(nèi)中頻聲壓級(jí)預(yù)測(cè)的有限元-統(tǒng)計(jì)能量分析混合模型;利用統(tǒng)計(jì)能量分析方法,建立了適用于車內(nèi)高頻聲壓級(jí)預(yù)測(cè)的統(tǒng)計(jì)能量分析模型。最后通過(guò)實(shí)車試驗(yàn)驗(yàn)證了3種模型。結(jié)果表明,在分析頻率范圍內(nèi)的大部分頻率下,駕駛員右耳側(cè)聲壓級(jí)的計(jì)算值與實(shí)測(cè)值的誤差在5dB以內(nèi),所建立的模型可用于下一步車內(nèi)聲場(chǎng)的分析和優(yōu)化。

    重型貨車;噪聲預(yù)測(cè);有限元法;邊界元法;統(tǒng)計(jì)能量分析

    前言

    目前廣泛應(yīng)用的預(yù)測(cè)車內(nèi)噪聲的數(shù)值計(jì)算方法主要有有限元方法(FEM)、邊界元方法(BEM)和統(tǒng)計(jì)能量分析(SEA)方法等。

    有限元方法在工程聲學(xué)的應(yīng)用始于20世紀(jì)60年代。文獻(xiàn)[1]中提出了一個(gè)聲場(chǎng)-結(jié)構(gòu)能量公式,開(kāi)創(chuàng)了有限元法在聲學(xué)領(lǐng)域應(yīng)用的先例。文獻(xiàn)[2]中最早應(yīng)用有限元方法預(yù)測(cè)車內(nèi)噪聲。盡管目前仍有很多研究者喜歡繼續(xù)使用純粹的有限元模型來(lái)研究車內(nèi)噪聲,但越來(lái)越多的人選擇使用混合模型,即結(jié)構(gòu)部分繼續(xù)使用有限元模型,內(nèi)部空腔部分使用邊界元模型。文獻(xiàn)[3]中綜合采用有限元和邊界元方法,建立了包括發(fā)動(dòng)機(jī)和傳動(dòng)系在內(nèi)的轎車聲學(xué)仿真模型,用于預(yù)測(cè)0-200Hz頻率范圍內(nèi)由路面激勵(lì)和發(fā)動(dòng)機(jī)懸置激勵(lì)產(chǎn)生的車內(nèi)結(jié)構(gòu)聲。

    統(tǒng)計(jì)能量分析的基本理論[4-5]于20世紀(jì)80年代末90年代初開(kāi)始應(yīng)用于汽車領(lǐng)域,文獻(xiàn)[6]是應(yīng)用統(tǒng)計(jì)能量分析方法進(jìn)行車內(nèi)噪聲預(yù)測(cè)的最早的文獻(xiàn)之一。

    有限元方法和統(tǒng)計(jì)能量分析方法各有其優(yōu)缺點(diǎn)和不同的適用范圍,將兩種方法有機(jī)結(jié)合、綜合運(yùn)用,以達(dá)到提高計(jì)算精度和速度的目的,是車內(nèi)噪聲預(yù)測(cè)領(lǐng)域一個(gè)新的發(fā)展方向。文獻(xiàn)[7]和文獻(xiàn)[8]中建立了混合FE-SEA模型,對(duì)車內(nèi)的結(jié)構(gòu)聲進(jìn)行了預(yù)測(cè)。文獻(xiàn)[9]中總結(jié)了有限元法和統(tǒng)計(jì)能量分析方法混合建模的理論,并分別從數(shù)值仿真和實(shí)驗(yàn)兩個(gè)方面對(duì)混合建模的計(jì)算結(jié)果進(jìn)行了驗(yàn)證。文獻(xiàn)[10]中對(duì)飛機(jī)的機(jī)艙壁板、地板和內(nèi)飾板等結(jié)構(gòu)進(jìn)行了混合建模,研究了飛機(jī)的結(jié)構(gòu)聲傳遞問(wèn)題。文獻(xiàn)[11]中綜合采用動(dòng)力平衡方程和功率平衡方程,用一種非迭代的確定性技術(shù)和統(tǒng)計(jì)技術(shù)混合方法對(duì)由屬性隨機(jī)或不確定的子系統(tǒng)組成的復(fù)雜聲振系統(tǒng)的集合平均穩(wěn)態(tài)響應(yīng)進(jìn)行了預(yù)測(cè)。數(shù)年前國(guó)內(nèi)也開(kāi)始了這方面的研究[12-13]。

    本文中針對(duì)某在研重型貨車車型,在較寬的頻率范圍內(nèi)建立駕駛室的聲學(xué)仿真模型,以實(shí)現(xiàn)低頻、中頻和高頻段的車內(nèi)聲壓級(jí)預(yù)測(cè),為下一步車內(nèi)聲場(chǎng)的分析和優(yōu)化提供基礎(chǔ),以期提高駕駛室的NVH性能。

    1 分析頻率范圍的確定

    參照《GB/T 18697—2002聲學(xué)汽車車內(nèi)噪聲測(cè)量方法》的規(guī)定,測(cè)量試驗(yàn)樣車的車內(nèi)噪聲,采樣頻率為8 192Hz。將測(cè)量的車內(nèi)噪聲信號(hào)轉(zhuǎn)換到頻域,頻率分辨率取1Hz。

    參考等效聲壓級(jí)的定義[14],計(jì)算頻率范圍f1~f2內(nèi)的總聲壓p∑:

    (1)

    式中pi為第i個(gè)頻率點(diǎn)處的有效聲壓。

    計(jì)算總聲壓級(jí):

    Lp=20lg(p∑/p0)

    (2)

    式中:p0為參考聲壓,p0=2×10-5Pa。

    汽車在各種工況下不同頻率范圍內(nèi)的總聲壓級(jí)見(jiàn)表1。

    表1 各種工況下不同頻率范圍內(nèi)的總聲壓級(jí)

    由表可見(jiàn):在各個(gè)工況下,20-1 000Hz范圍內(nèi)的總聲壓級(jí)僅比20-4 096Hz范圍內(nèi)的總聲壓級(jí)小0.001~0.02dB,該重型貨車車內(nèi)噪聲的主要能量集中在20-1 000Hz這個(gè)頻率范圍內(nèi)。因此,研究的頻率范圍定為20-1 000Hz。

    2 車內(nèi)低頻聲壓級(jí)預(yù)測(cè)模型

    根據(jù)有限元方法和邊界元方法的一般適用頻率范圍以及仿真模型的預(yù)計(jì)規(guī)模,在20-200Hz的低頻范圍內(nèi),綜合運(yùn)用有限元方法和邊界元方法來(lái)預(yù)測(cè)車內(nèi)聲壓級(jí)。

    建立駕駛室的結(jié)構(gòu)有限元模型是進(jìn)行車內(nèi)聲場(chǎng)仿真的基礎(chǔ)。文中研究的駕駛室是一個(gè)由骨架、覆蓋件、車門、風(fēng)窗玻璃及大量?jī)?nèi)飾件、外飾件、功能件等零部件通過(guò)不同方式連接形成的復(fù)雜結(jié)構(gòu),內(nèi)部還有座椅、臥鋪、儲(chǔ)物箱等結(jié)構(gòu)。因此,遵循從簡(jiǎn)到繁、逐步完善、逐步驗(yàn)證的理念,首先建立白車身的有限元模型,并通過(guò)模態(tài)試驗(yàn)進(jìn)行驗(yàn)證,然后向白車身模型添加座椅、臥鋪等結(jié)構(gòu),最終建立整備駕駛室模型,同樣通過(guò)模態(tài)試驗(yàn)進(jìn)行驗(yàn)證。

    最后建成的整備駕駛室結(jié)構(gòu)有限元模型如圖1所示。其中,單元總數(shù)為162 841,節(jié)點(diǎn)總數(shù)為123 848。

    基于駕駛室有限元模型,建立駕駛室的聲學(xué)邊界元模型。計(jì)算駕駛室的結(jié)構(gòu)模態(tài),并將結(jié)果導(dǎo)入駕駛室的聲學(xué)邊界元模型中。然后定義激勵(lì)力,即可根據(jù)模態(tài)疊加法計(jì)算駕駛室結(jié)構(gòu)件的振動(dòng)速度,進(jìn)而計(jì)算車內(nèi)聲壓級(jí)。

    3 車內(nèi)中頻聲壓級(jí)預(yù)測(cè)模型

    3.1 中頻范圍的確定

    首先計(jì)算各結(jié)構(gòu)件在200-1 000Hz各1/3倍頻程內(nèi)的模態(tài)數(shù)。采用Lanczos方法,計(jì)算駕駛室每一個(gè)結(jié)構(gòu)的有限元模型在研究頻帶內(nèi)的模態(tài)數(shù)。表2列出部分骨架和壁板的計(jì)算結(jié)果。

    表2 部分結(jié)構(gòu)在分析頻段內(nèi)的模態(tài)數(shù)

    由表可見(jiàn):在500Hz以上的各1/3倍頻程內(nèi),全部部件的模態(tài)數(shù)均大于5,而在500Hz以下的各1/3倍頻程內(nèi),部分骨架的模態(tài)數(shù)小于5。因此,中頻范圍確定為200-500Hz。

    在200-500Hz頻率范圍內(nèi),采用有限元-統(tǒng)計(jì)能量分析混合仿真建模。對(duì)在這個(gè)頻率范圍的各1/3倍頻程內(nèi)的模態(tài)數(shù)均大于5的結(jié)構(gòu),如壁板,用統(tǒng)計(jì)能量分析方法建模,對(duì)不能滿足這個(gè)條件的結(jié)構(gòu),如絕大多數(shù)骨架,則建立有限元模型。

    3.2 結(jié)構(gòu)有限元模型

    基于建立的整備駕駛室有限元模型,從中逐一選取不滿足在200-500Hz各1/3倍頻程內(nèi)的模態(tài)數(shù)均大于5這個(gè)條件的結(jié)構(gòu)的有限元模型,將其建立為FE子系統(tǒng)。在轉(zhuǎn)換過(guò)程中,原結(jié)構(gòu)的屬性自動(dòng)保存在對(duì)應(yīng)的FE子系統(tǒng)中,不再需要單獨(dú)設(shè)置。建成后的結(jié)構(gòu)有限元模型由37個(gè)FE子系統(tǒng)構(gòu)成,如圖2所示。

    3.3 結(jié)構(gòu)統(tǒng)計(jì)能量分析模型

    3.3.1 結(jié)構(gòu)子系統(tǒng)建模

    對(duì)結(jié)構(gòu)子系統(tǒng)建模時(shí),遵循以下基本原則:

    (1) 子系統(tǒng)在分析頻帶內(nèi)的模態(tài)數(shù)不低于5;

    (2) 子系統(tǒng)之間通過(guò)點(diǎn)、線或面等形式較好地進(jìn)行耦合,保證能量的正確傳遞;

    (3) 外部激勵(lì)作為模型的外部能量輸入,而不視為單獨(dú)的子系統(tǒng);

    (4) 簡(jiǎn)單的板件結(jié)構(gòu)建為平板或曲面板子系統(tǒng),近似梁結(jié)構(gòu)則建為梁子系統(tǒng)。

    根據(jù)以上建模原則,對(duì)駕駛室的各個(gè)壁板建立統(tǒng)計(jì)能量分析模型,共計(jì)33個(gè)子系統(tǒng),如圖3所示。

    3.3.2 參數(shù)的確定

    (1) 模態(tài)密度的確定

    統(tǒng)計(jì)能量分析方法中的模態(tài)密度是指子系統(tǒng)在分析頻率范圍單位帶寬內(nèi)的模態(tài)數(shù)。它描述了子系統(tǒng)貯存能量的能力,同時(shí)也描述了子系統(tǒng)從外界接受能量產(chǎn)生振動(dòng)的能力,對(duì)進(jìn)一步確定子系統(tǒng)響應(yīng)能量與功率流的關(guān)系具有非常重要的作用。

    對(duì)桿、梁、板、圓柱殼等構(gòu)型簡(jiǎn)單的結(jié)構(gòu)和簡(jiǎn)化空腔,其模態(tài)密度的確定基本是個(gè)數(shù)學(xué)問(wèn)題,計(jì)算方法及公式可參閱文獻(xiàn)[15]第22-42頁(yè)。

    對(duì)構(gòu)型復(fù)雜的結(jié)構(gòu),其模態(tài)密度的解析計(jì)算則很困難,一般通過(guò)試驗(yàn)確定。而當(dāng)結(jié)構(gòu)的模態(tài)密度較大或者試驗(yàn)中固有誤差較大時(shí),通過(guò)試驗(yàn)方法測(cè)得的模態(tài)密度的誤差也會(huì)較大。

    本文中基于數(shù)值仿真方法確定各個(gè)子系統(tǒng)的模態(tài)密度,具體做法為:基于子系統(tǒng)對(duì)應(yīng)結(jié)構(gòu)的有限元模型,采用Lanczos方法計(jì)算結(jié)構(gòu)的模態(tài),分析頻率范圍為0-1 500Hz。由于統(tǒng)計(jì)能量分析方法只關(guān)心分析頻帶內(nèi)振型數(shù)目的多少,不關(guān)心與邊界條件有關(guān)的子系統(tǒng)具體模態(tài)參數(shù)的大小,即模態(tài)密度與邊界條件無(wú)關(guān),這也是統(tǒng)計(jì)能量分析中“統(tǒng)計(jì)”的依據(jù)之一[15]57,因此,在計(jì)算模態(tài)時(shí),子系統(tǒng)對(duì)應(yīng)結(jié)構(gòu)的有限元模型不施加約束。

    設(shè)Δω為分析的頻帶帶寬,Nm為頻帶內(nèi)的振型數(shù)目,ω為頻帶中心頻率,則模態(tài)密度n(ω)為

    (3)

    (2) 內(nèi)損耗因子的確定

    在統(tǒng)計(jì)能量分析方法中,內(nèi)損耗因子η用來(lái)表示子系統(tǒng)的阻尼損耗特性[15]68,其定義為子系統(tǒng)在單位頻率內(nèi)單位時(shí)間的損耗能量與子系統(tǒng)平均儲(chǔ)存能量的比值。

    雖然各種材料和結(jié)構(gòu)的阻尼損耗數(shù)據(jù)有表格和計(jì)算公式可查,但由于系統(tǒng)連接方式、幾何形狀等因素對(duì)結(jié)構(gòu)動(dòng)力學(xué)的影響,獲取內(nèi)損耗因子的最常用和可靠的方法還是通過(guò)試驗(yàn)測(cè)量。測(cè)量結(jié)構(gòu)內(nèi)損耗因子的方法主要有模態(tài)法(又稱半功率點(diǎn)帶寬法)、脈沖響應(yīng)衰減法和能量注入法。綜合文獻(xiàn)[16]和文獻(xiàn)[17]可知,采用能量注入法能夠獲取精度滿足工程需要的內(nèi)損耗因子。

    在車輛正向設(shè)計(jì)早期的立項(xiàng)、方案設(shè)計(jì)和設(shè)計(jì)開(kāi)發(fā)等階段,很難甚至不可能得到各個(gè)子系統(tǒng)對(duì)應(yīng)的實(shí)際結(jié)構(gòu)件,從而無(wú)法通過(guò)試驗(yàn)測(cè)量實(shí)際結(jié)構(gòu)的內(nèi)損耗因子,因此本文中根據(jù)能量注入法原理,采用數(shù)值計(jì)算的方法計(jì)算子系統(tǒng)的內(nèi)損耗因子。

    具體做法為:基于結(jié)構(gòu)的有限元模型,建立對(duì)應(yīng)的FE子系統(tǒng)。在轉(zhuǎn)換過(guò)程中,根據(jù)文獻(xiàn)[15]第69頁(yè),對(duì)材料為鋼的結(jié)構(gòu),其結(jié)構(gòu)損耗因子取為4×10-4,對(duì)材料為玻璃的結(jié)構(gòu),其結(jié)構(gòu)損耗因子取為1×10-3。然后,在FE子系統(tǒng)上盡量選擇剛度大、且不是結(jié)構(gòu)各階振型節(jié)點(diǎn)的位置施加激勵(lì)力。在激勵(lì)點(diǎn)附近布置加速度傳感器測(cè)量該點(diǎn)的振動(dòng)速度,同時(shí)在結(jié)構(gòu)子系統(tǒng)的其他位置盡可能多地均勻布置加速度傳感器,對(duì)各個(gè)測(cè)點(diǎn)的振動(dòng)速度求均值得到子系統(tǒng)的平均振動(dòng)速度。計(jì)算的頻率范圍為150-1 500Hz。

    以中間地板為例,通過(guò)數(shù)值計(jì)算獲取其內(nèi)損耗因子的方法如圖4所示。

    內(nèi)損耗因子的計(jì)算公式為

    (4)

    式中:F1為激勵(lì)力;v1為激勵(lì)點(diǎn)附近的振動(dòng)速度;m為結(jié)構(gòu)質(zhì)量;〈v2〉為空間平均的振動(dòng)速度均方值。駕駛室部分結(jié)構(gòu)的內(nèi)損耗因子的計(jì)算結(jié)果見(jiàn)表3。

    表3 部分結(jié)構(gòu)內(nèi)損耗因子的數(shù)值計(jì)算結(jié)果

    (3) 耦合損耗因子的確定

    在統(tǒng)計(jì)能量分析方法中,耦合損耗因子用來(lái)描述兩個(gè)子系統(tǒng)之間耦合程度的強(qiáng)弱。按照子系統(tǒng)的性質(zhì),子系統(tǒng)之間的耦合可分為結(jié)構(gòu)-結(jié)構(gòu)耦合、結(jié)構(gòu)-空腔耦合和空腔-空腔耦合。駕駛室結(jié)構(gòu)之間最常見(jiàn)的耦合形式是結(jié)構(gòu)與結(jié)構(gòu)之間的機(jī)械連接,為便于分析,可簡(jiǎn)化為點(diǎn)、線和面等幾種典型的連接接觸形式。結(jié)構(gòu)與空腔之間和空腔與空腔之間的耦合形式一般都是面連接接觸。

    通過(guò)試驗(yàn)方法測(cè)量耦合損耗因子非常困難,這是因?yàn)楫?dāng)子結(jié)構(gòu)之間是弱耦合時(shí),耦合損耗因子比子系統(tǒng)本身的內(nèi)損耗因子至少要小一個(gè)數(shù)量級(jí),在試驗(yàn)中需要處理數(shù)量級(jí)為1×10-3~1×10-4的小數(shù)及其差數(shù),因此對(duì)測(cè)量技術(shù)和數(shù)據(jù)處理的要求很高。此外,對(duì)由N個(gè)子系統(tǒng)組成的系統(tǒng),需要做N次獨(dú)立試驗(yàn)。對(duì)駕駛室這樣的復(fù)雜系統(tǒng),試驗(yàn)的工作量非常大?;谏鲜鲈颍疚闹兄饕捎美碚摴接?jì)算駕駛室各個(gè)子系統(tǒng)之間的耦合損耗因子。

    駕駛室中的很多結(jié)構(gòu)可以用以點(diǎn)、線、面等形式連接的子系統(tǒng)來(lái)表示。以點(diǎn)、線、面等形式連接的兩結(jié)構(gòu)間的耦合損耗因子可通過(guò)理論公式算得,詳見(jiàn)文獻(xiàn)[15]第80-114頁(yè)。

    但還有一些結(jié)構(gòu),如風(fēng)窗玻璃和周圍壁板,它們之間的耦合損耗因子需要通過(guò)特別的方法進(jìn)行計(jì)算。風(fēng)窗玻璃與周圍壁板之間通過(guò)粘膠連接。在有限元模型中,玻璃、壁板和粘膠都可以選用合適的單元類型進(jìn)行建模。在統(tǒng)計(jì)能量分析模型中,玻璃和壁板可以選用薄板類型的子系統(tǒng)建模,而粘膠則沒(méi)有合適的子系統(tǒng)類型與其對(duì)應(yīng)。針對(duì)這個(gè)問(wèn)題,本文中直接把玻璃和壁板通過(guò)線連接的形式耦合起來(lái)。這樣,就出現(xiàn)了一個(gè)問(wèn)題,即:如何保證這種連接方式的能量傳遞關(guān)系與實(shí)際結(jié)構(gòu)等效,亦即它們的耦合損耗因子等效。本文中采用了能量流方法來(lái)解決這一問(wèn)題。

    能量流法提供了一種利用有限元方法計(jì)算通過(guò)復(fù)雜結(jié)構(gòu)連接的子系統(tǒng)之間的等效耦合損耗因子的途徑,以提高統(tǒng)計(jì)能量分析方法對(duì)結(jié)構(gòu)聲的預(yù)測(cè)精度。該方法可視為能量注入方法的數(shù)值實(shí)現(xiàn)。

    風(fēng)窗玻璃與周圍壁板之間的耦合損耗因子的具體計(jì)算過(guò)程如下。

    基于風(fēng)窗玻璃、壁板和膠條的有限元模型,分別建立其對(duì)應(yīng)的FE子系統(tǒng)。將起連接作用的粘膠FE子系統(tǒng)的狀態(tài)設(shè)為“禁用”,其余兩個(gè)FE子系統(tǒng)的狀態(tài)設(shè)為“啟用”,如圖5所示。

    分析的頻率范圍設(shè)為150-1 500Hz。首先計(jì)算各子系統(tǒng)的模態(tài)、質(zhì)量矩陣和剛度矩陣等,然后在一個(gè)狀態(tài)設(shè)為“啟用”的子系統(tǒng)上施加rain-on-the-roof(雨滴屋頂)類型的激勵(lì)。這種類型的激勵(lì)是空間不相干的寬帶激勵(lì),其在任意位置的幅值與局部質(zhì)量密度成比例。然后,計(jì)算各子系統(tǒng)的頻帶平均輸入功率、系統(tǒng)的頻帶平均能量影響系數(shù)矩陣,得到如下關(guān)系式:

    E(ω0)=MEIC(ω0)Pin(ω0)

    (5)

    式中:E(ω0)為子系統(tǒng)能量矩陣;Pin(ω0)為各子系統(tǒng)的輸入功率矩陣;ω0為圓頻率;MEIC(ω0)為能量影響系數(shù)矩陣,其第m行第n列的元素表示的是向第n個(gè)子系統(tǒng)輸入單位功率時(shí),第m個(gè)子系統(tǒng)中的能量。

    MEIC矩陣與施加給系統(tǒng)的激勵(lì)類型有關(guān)。在很多情況下,對(duì)可能的各種類型的激勵(lì)進(jìn)行集合平均以后的能量流模型會(huì)更有用。施加rain-on-the-roof類型的激勵(lì)的能量流模型與進(jìn)行激勵(lì)平均的能量流模型是等效的,其MEIC矩陣比某種特定激勵(lì)下的模型更能描述系統(tǒng)的能量流屬性。

    由式(5)可得

    Pin(ω0)=MEIC(ω0)-1E(ω0)

    (6)

    MEIC(ω0)逆矩陣的非對(duì)角線元素可視為描述子系統(tǒng)之間能量傳遞關(guān)系的等效耦合損耗因子。使用“等效”一詞是因?yàn)橥ㄟ^(guò)這種方式得到的耦合損耗因子是確定性的,而統(tǒng)計(jì)能量分析方法中使用的耦合損耗因子是集合平均的統(tǒng)計(jì)量。

    3.4 空腔統(tǒng)計(jì)能量分析模型

    3.4.1 建模方法

    駕駛室內(nèi)部的空腔由壁板、骨架、車門、車門玻璃和風(fēng)窗玻璃等結(jié)構(gòu)包圍而成。在車輛處于非穩(wěn)態(tài)工況、車內(nèi)聲場(chǎng)快速變化等情況下,車內(nèi)不同位置的聲場(chǎng)之間可能存在差異,因此不能簡(jiǎn)單把車內(nèi)聲場(chǎng)整體視為一個(gè)混響場(chǎng),需要將其內(nèi)部空間劃分為多個(gè)空腔子系統(tǒng),每個(gè)空腔子系統(tǒng)可用混響場(chǎng)來(lái)近似。劃分空腔子系統(tǒng)的原則為:

    (1) 保證每個(gè)空腔子系統(tǒng)在分析頻帶內(nèi)的模態(tài)數(shù)滿足統(tǒng)計(jì)能量分析方法的要求,即不小于5;

    (2) 根據(jù)所耦合的結(jié)構(gòu)子系統(tǒng)劃分,保證正確的能量傳遞關(guān)系;

    (3) 能比較準(zhǔn)確地反映所要研究位置的聲場(chǎng)情況。

    根據(jù)以上原則,將駕駛室內(nèi)部的空間劃分為9個(gè)空腔子系統(tǒng),如圖6所示。

    3.4.2 參數(shù)確定

    (1) 空腔子系統(tǒng)模態(tài)密度的確定

    在空腔的體積、總表面積、總棱邊長(zhǎng)度和分析頻率確定后,即可計(jì)算空腔的模態(tài)密度,計(jì)算方法詳見(jiàn)文獻(xiàn)[15]第38-39頁(yè)。

    (2) 空腔子系統(tǒng)內(nèi)損耗因子的確定

    根據(jù)空腔內(nèi)部損耗的形成機(jī)理,空腔的內(nèi)損耗因子可以通過(guò)計(jì)算空腔的平均吸聲系數(shù)得到。設(shè)空腔由表面Si(i=1,2,…,J)包圍而成,各表面的吸聲系數(shù)為αi,則空腔的平均吸聲系數(shù)為

    (7)

    由于本文中分析頻率的上限為1 000Hz,且駕駛室內(nèi)空間不大,因此不考慮空氣的聲吸收[18]426。

    已知空腔的平均吸聲系數(shù)后,可根據(jù)以下經(jīng)驗(yàn)公式計(jì)算混響聲場(chǎng)的混響時(shí)間[18]423:

    (8)

    式中:S為空腔的總吸聲表面積;V為空腔體積。進(jìn)而計(jì)算空腔的內(nèi)損耗因子:

    (9)

    3.5 激勵(lì)和約束

    重型貨車駕駛室受到的激勵(lì)包括來(lái)自發(fā)動(dòng)機(jī)艙的聲輻射激勵(lì)、動(dòng)力總成的振動(dòng)激勵(lì)、路面不平度產(chǎn)生的振動(dòng)激勵(lì)和車外的風(fēng)激勵(lì)(風(fēng)噪)。

    由于重型貨車的行駛速度一般都在100km/h以下,風(fēng)噪對(duì)車內(nèi)噪聲的影響并不顯著,因此建模時(shí),暫不考慮車外的風(fēng)激勵(lì)。

    對(duì)來(lái)自動(dòng)力總成的振動(dòng)激勵(lì)和路面不平度產(chǎn)生的振動(dòng)激勵(lì),可用駕駛室前后懸置點(diǎn)附近的4個(gè)力來(lái)等效。由于來(lái)自駕駛室懸置系統(tǒng)的等效激勵(lì)力分別作用于駕駛室的左、右地板上,因此分別實(shí)測(cè)駕駛室左、右地板上前后懸置激勵(lì)點(diǎn)處的振動(dòng)加速度,平均后作為響應(yīng)約束施加在對(duì)應(yīng)左、右地板的SEA子系統(tǒng)上,如圖7所示。

    對(duì)來(lái)自發(fā)動(dòng)機(jī)艙的聲輻射激勵(lì)的建模方法是:在發(fā)動(dòng)機(jī)上方的兩端和中間分別放置傳聲器,將各個(gè)傳聲器實(shí)測(cè)的聲壓值平均后作為擴(kuò)散聲場(chǎng)的聲壓,施加在駕駛室各地板子系統(tǒng)上,見(jiàn)圖7。

    4 車內(nèi)高頻聲壓級(jí)預(yù)測(cè)模型

    由第3.1節(jié)可知,在630-1 000Hz范圍內(nèi)的各1/3倍頻程內(nèi),駕駛室各結(jié)構(gòu)的模態(tài)數(shù)均滿足統(tǒng)計(jì)能量分析方法的要求,因此都建為統(tǒng)計(jì)能量分析的子系統(tǒng)。子系統(tǒng)的劃分原則、建模過(guò)程、子系統(tǒng)的參數(shù)確定方法與激勵(lì)和約束的施加方法均與第3.4節(jié)相同。

    建立的結(jié)構(gòu)統(tǒng)計(jì)能量分析模型共包括70個(gè)子系統(tǒng),如圖8所示。

    車內(nèi)高頻聲壓級(jí)預(yù)測(cè)模型的空腔統(tǒng)計(jì)能量分析模型與第3.4節(jié)中的空腔模型完全相同。

    5 預(yù)測(cè)模型的驗(yàn)證

    采用VA One軟件作為聲學(xué)仿真平臺(tái),有限元建模選擇Altair Hypermesh軟件,動(dòng)力學(xué)計(jì)算選用MSC.Nastran軟件。

    選用具有代表性的怠速(650r/min)、勻速(16擋90km/h)和加速(8擋1 000r/min進(jìn)線、2 200r/min出線)等3種工況,在200-1 000Hz的頻率范圍內(nèi),對(duì)駕駛員右耳側(cè)聲壓級(jí)的計(jì)算值與實(shí)測(cè)值進(jìn)行對(duì)比。圖9~圖11分別為20-200Hz范圍內(nèi)車內(nèi)聲壓級(jí)(線性頻譜)、20-200Hz范圍內(nèi)車內(nèi)聲壓級(jí)(1/3倍頻程)和200-1 000Hz范圍內(nèi)車內(nèi)聲壓級(jí)(1/3倍頻程)計(jì)算值與實(shí)測(cè)值的比較。

    由圖9~圖11可見(jiàn):在20-1 000Hz的范圍內(nèi),在大部分頻率點(diǎn)或頻帶內(nèi),計(jì)算值與實(shí)測(cè)值的誤差在5dB以內(nèi),表明所建立的低頻、中頻和高頻車內(nèi)聲壓級(jí)預(yù)測(cè)模型可用于車內(nèi)聲壓級(jí)的仿真計(jì)算。

    6 結(jié)論

    本文中綜合運(yùn)用各種數(shù)值計(jì)算方法,結(jié)合理論分析和試驗(yàn)等手段,在較寬的頻率范圍內(nèi)建立了某重型貨車車型駕駛室的聲壓級(jí)預(yù)測(cè)模型,著重解決了子系統(tǒng)的科學(xué)劃分、模型參數(shù)(特別是內(nèi)損耗因子和耦合損耗因子)的確定、多種激勵(lì)的施加等問(wèn)題,實(shí)現(xiàn)了低頻、中頻和高頻段的車內(nèi)聲壓級(jí)預(yù)測(cè)。

    對(duì)比結(jié)果表明,所建立的聲壓級(jí)預(yù)測(cè)模型的計(jì)算精度和速度能夠滿足工程需要,可用于下一步車內(nèi)聲場(chǎng)的分析、比較和改進(jìn)。

    [1] GLADWELL G M L, ZIMMERMANN G. On Energy and Complementary Energy Formulations of Acoustic and Structural Vibration Problems[J]. Journal of Sound and Vibration,1966,3(3):233-241.

    [2] NEFSKE D J, WOLF Jr J A, HOWELL L J. Structural-acoustic Finite Element Analysis of the Automobile Passenger Compartment: a Review of Current practice[J]. Journal of Sound and Vibration,1982,80(2):247-266.

    [3] SOL A, Van HERPE F. Numerical Prediction of a Whole Car Vibro-Acoustic Behavior at Low Frequencies[C]. SAE Paper 2001-01-1521.

    [4] LYON R H, MAIDANIK G. Power Flow Between Linearly Coupled Oscillators[J]. The Journal of the Acoustical Society of America,1962,34(5):623-639.

    [5] SMITH Jr P W. Response and Radiation of Structural Modes Excited by Sound[J]. The Journal of the Acoustical Society of America,1962,34(5):640-647.

    [6] DEJONG R G. A Study of Vehicle Interior Noise Using Statistical Energy Analysis[C]. SAE Paper 850960.

    [7] CHARPENTIER A, SREEDHAR P, GARDNER B, et al. Use of a Hybrid FE-SEA Model of a Trimmed Vehicle to Improve the Design for Interior Noise[C]. SAE Paper 2009-01-2199.

    [8] CHARPENTIER A, SREEDHAR P, FUKUI K. Using the Hybrid FE-SEA Method to Predict Structure-borne Noise Transmission in a Trimmed Automotive Vehicle[C]. SAE Paper 2007-01-2181.

    [9] COTONI V, GARDNER B, SHORTER P, et al. Demonstration of Hybrid FE-SEA Analysis of Structure-borne Noise in the Mid Frequency Range[C]. SAE Paper 2005-01-2331.

    [10] COTONI V, GARDNER B, CORDIOLI J A, et al. Advanced Modeling of Aircraft Interior Noise Using the Hybrid FE-SEA Method[C]. SAE Paper 2008-36-0575.

    [11] SHORTER P J, LANGLEY R S. Vibro-Acoustic Analysis of Complex Systems[J]. Journal of Sound and Vibration,2005,288(3):669-699.

    [12] 陳書(shū)明,王登峰,昝建明.基于FE-SEA混合模型的轎車車內(nèi)噪聲預(yù)測(cè)[J].汽車工程,2011,33(3):236-240.

    [13] 楊雄偉,李躍明,耿謙.基于混合FE-SEA法的高溫環(huán)境飛行器寬頻聲振特性分析[J].航空學(xué)報(bào),2011,32(10):1851-1859.

    [14] TORIJA A J, RUIZ D P, Ramos-Ridaongel. Required Stabilization Time, Short-Term Variability and Impulsiveness of the Sound Pressure Level to Characterize the Temporal Composition of Urban Soundscapes[J]. Applied Acoustics,2011,72(2-3):90.

    [15] 姚德源,王其政.統(tǒng)計(jì)能量分析原理及其應(yīng)用[M].北京:北京理工大學(xué)出版社,1995.

    [16] BLOSS B C, RAO M D. Estimation of Frequency-Averaged Loss Factors by the Power Injection and the Impulse Response Decay Methods[J]. The Journal of the Acoustical Society of America,2005,117(1):240-249.

    [17] De LANGHE K. High Frequency Vibrations: Contributions to Experimental and Computational SEA Parameter Identification Techniques[D]. Leuven:Katholieke Universiteit,1996.

    [18] 杜功煥,朱哲民,龔秀芬.聲學(xué)基礎(chǔ)[M].南京:南京大學(xué)出版社,2001.

    Prediction on the Interior Sound Pressure Level of a Heavy-duty Truck

    Wang Binxing1, Lian Xiaomin2, Zheng Sifa2, Ren Yanxi1& Liu Yu1

    1.No.2InstituteofEngineeringCorps,GeneralEquipmentsDepartment,Beijing100093; 2.TsinghuaUniversity,StateKeyLaboratoryofAutomotiveSafetyandEnergy,Beijing100084

    The interior sound pressure level of a heavy duty truck at the frequency range of 20-1000Hz is predicted in this paper. Firstly, the frequency range to be analyzed is determined according to the energy distribution of vehicle interior noise measured under various conditions. Then based on the verified finite element model for cabin and with the integrated utilization of finite element method (FEM) and boundary element method (BEM), a finite element / boundary element model for predicting vehicle interior low-frequency noise is built. With the integrated utilization of FEM and statistical energy analysis (SEA) method, a hybrid finite element / statistical energy analysis model is set up for predicting vehicle interior medium-frequency noise. And a SEA model is also created for predicting vehicle interior high-frequency noise by using SEA method. Finally three models are all verified by real vehicle tests. The results show that for the most part of analysis frequency range, the difference between the calculated value and the measured value is within 5dB, indicating that the models built can be used for the further analysis and optimization of vehicle interior sound field.

    heavy duty truck; noise prediction; FEM; BEM; statistical energy analysis

    *國(guó)家自然科學(xué)基金(51275262)資助。

    原稿收到日期為2014年7月28日,修改稿收到日期為2014年11月17日。

    猜你喜歡
    聲壓級(jí)駕駛室損耗
    機(jī)器噪聲平均聲壓級(jí)計(jì)算方法差異性實(shí)證研究
    淺析8D方法在解決H7車型駕駛室振動(dòng)的應(yīng)用
    基于PID控制的載貨汽車駕駛室半主動(dòng)懸置控制
    一種計(jì)算消聲室聲壓級(jí)的新方法
    全新DXR mkll有源揚(yáng)聲器
    演藝科技(2019年4期)2019-03-30 03:21:46
    自我損耗理論視角下的編輯審讀
    新聞傳播(2016年11期)2016-07-10 12:04:01
    變壓器附加損耗對(duì)負(fù)載損耗的影響
    非隔離型單相光伏并網(wǎng)逆變器的功率損耗研究
    Diodes1.9W D類音頻放大器提供高聲壓級(jí)水平并延長(zhǎng)電池壽命
    輕型卡車駕駛室內(nèi)飾平臺(tái)化設(shè)計(jì)
    中文字幕精品免费在线观看视频| 午夜福利欧美成人| 国产主播在线观看一区二区| 国产亚洲欧美在线一区二区| 久久人妻福利社区极品人妻图片| 国产麻豆成人av免费视频| 欧美黄色淫秽网站| 伦理电影免费视频| 午夜老司机福利片| 日日爽夜夜爽网站| 一边摸一边做爽爽视频免费| 国内少妇人妻偷人精品xxx网站 | 欧美不卡视频在线免费观看 | 国内揄拍国产精品人妻在线 | 成年免费大片在线观看| 欧美中文综合在线视频| 曰老女人黄片| 亚洲精品一区av在线观看| 久热这里只有精品99| 亚洲专区中文字幕在线| 88av欧美| 免费看美女性在线毛片视频| 日本a在线网址| 一本一本综合久久| 中文字幕久久专区| 国产高清视频在线播放一区| 天堂动漫精品| 99在线视频只有这里精品首页| 精品日产1卡2卡| 亚洲 欧美一区二区三区| 亚洲成人免费电影在线观看| 精华霜和精华液先用哪个| 国产欧美日韩精品亚洲av| 叶爱在线成人免费视频播放| 少妇的丰满在线观看| 午夜免费鲁丝| 亚洲国产高清在线一区二区三 | 免费搜索国产男女视频| 免费看日本二区| 国产熟女xx| 亚洲精品在线美女| 91麻豆精品激情在线观看国产| 母亲3免费完整高清在线观看| 在线观看免费午夜福利视频| 欧美大码av| 超碰成人久久| 国产亚洲精品第一综合不卡| 亚洲七黄色美女视频| 99精品久久久久人妻精品| 亚洲国产中文字幕在线视频| 欧美精品亚洲一区二区| 欧美av亚洲av综合av国产av| 黄色视频,在线免费观看| 欧美日韩一级在线毛片| 欧美最黄视频在线播放免费| 精品不卡国产一区二区三区| 久久天堂一区二区三区四区| 久久精品91无色码中文字幕| 欧美中文日本在线观看视频| 亚洲av五月六月丁香网| 日韩精品免费视频一区二区三区| 性色av乱码一区二区三区2| 香蕉国产在线看| 50天的宝宝边吃奶边哭怎么回事| 亚洲一区二区三区色噜噜| 美女免费视频网站| 日本a在线网址| 久久久精品国产亚洲av高清涩受| 免费观看精品视频网站| 国产精品久久久久久亚洲av鲁大| 天堂动漫精品| 中文字幕精品亚洲无线码一区 | 欧美日韩黄片免| 岛国视频午夜一区免费看| 中亚洲国语对白在线视频| 国产精品亚洲av一区麻豆| 好看av亚洲va欧美ⅴa在| 国产一区二区三区视频了| 日韩欧美一区二区三区在线观看| 久久久久国内视频| 欧美大码av| 美女国产高潮福利片在线看| 久久亚洲精品不卡| 免费女性裸体啪啪无遮挡网站| 他把我摸到了高潮在线观看| 精品久久久久久久末码| 中文在线观看免费www的网站 | 久久久久国产精品人妻aⅴ院| www.自偷自拍.com| 亚洲av熟女| xxxwww97欧美| 久久久久久久久中文| 午夜福利在线在线| 黄色视频,在线免费观看| 在线av久久热| 夜夜躁狠狠躁天天躁| 夜夜夜夜夜久久久久| 精品国内亚洲2022精品成人| 亚洲专区中文字幕在线| 久久精品影院6| 午夜免费激情av| 欧美久久黑人一区二区| 不卡一级毛片| 国产av不卡久久| 人人妻人人澡人人看| 久久久久久久久免费视频了| 午夜两性在线视频| 亚洲熟妇熟女久久| 在线观看免费午夜福利视频| 日本一本二区三区精品| 精品电影一区二区在线| 听说在线观看完整版免费高清| 午夜老司机福利片| 久久精品亚洲精品国产色婷小说| 黄色 视频免费看| 美女高潮到喷水免费观看| 久久午夜亚洲精品久久| 国内精品久久久久精免费| 亚洲全国av大片| 亚洲精品国产精品久久久不卡| 很黄的视频免费| 波多野结衣巨乳人妻| 少妇粗大呻吟视频| 视频区欧美日本亚洲| av福利片在线| 久久国产亚洲av麻豆专区| 精品人妻1区二区| 琪琪午夜伦伦电影理论片6080| 久久久国产成人免费| 久久国产亚洲av麻豆专区| 一进一出抽搐gif免费好疼| 亚洲精品色激情综合| av福利片在线| 亚洲国产日韩欧美精品在线观看 | 18禁裸乳无遮挡免费网站照片 | 91字幕亚洲| 亚洲精品美女久久av网站| 国产成人一区二区三区免费视频网站| 亚洲一区中文字幕在线| 亚洲,欧美精品.| 在线观看日韩欧美| 亚洲专区字幕在线| 一a级毛片在线观看| 一边摸一边抽搐一进一小说| 啦啦啦观看免费观看视频高清| 久久国产精品影院| 波多野结衣高清作品| 亚洲五月婷婷丁香| 精品一区二区三区av网在线观看| 神马国产精品三级电影在线观看 | 性色av乱码一区二区三区2| 日本在线视频免费播放| 国产精品免费视频内射| 男女做爰动态图高潮gif福利片| 桃红色精品国产亚洲av| 啦啦啦观看免费观看视频高清| 午夜免费成人在线视频| av视频在线观看入口| 老司机午夜福利在线观看视频| 男男h啪啪无遮挡| 脱女人内裤的视频| 欧美又色又爽又黄视频| 国产一区二区在线av高清观看| 国产精品 国内视频| 国产精品永久免费网站| 亚洲中文字幕日韩| a级毛片a级免费在线| 日本一区二区免费在线视频| 亚洲国产欧美网| 欧美在线黄色| 日韩高清综合在线| 看片在线看免费视频| 韩国精品一区二区三区| 黄色视频,在线免费观看| 色尼玛亚洲综合影院| 侵犯人妻中文字幕一二三四区| 久久婷婷人人爽人人干人人爱| 欧美色视频一区免费| 免费在线观看日本一区| 久久欧美精品欧美久久欧美| 久久久久国内视频| 丁香六月欧美| 搡老岳熟女国产| 欧美在线一区亚洲| 亚洲天堂国产精品一区在线| 999精品在线视频| 长腿黑丝高跟| 精品久久久久久成人av| 不卡一级毛片| 精品国产乱码久久久久久男人| 国内久久婷婷六月综合欲色啪| 黑人巨大精品欧美一区二区mp4| 美女免费视频网站| 一级作爱视频免费观看| 欧美 亚洲 国产 日韩一| 91麻豆精品激情在线观看国产| 最近在线观看免费完整版| 免费在线观看视频国产中文字幕亚洲| 免费在线观看影片大全网站| 免费女性裸体啪啪无遮挡网站| 色av中文字幕| 亚洲电影在线观看av| 听说在线观看完整版免费高清| 欧美成人免费av一区二区三区| 精品国产亚洲在线| 精品久久久久久久毛片微露脸| 99久久精品国产亚洲精品| 变态另类成人亚洲欧美熟女| 19禁男女啪啪无遮挡网站| 91麻豆av在线| 日本成人三级电影网站| 岛国视频午夜一区免费看| 国产欧美日韩一区二区精品| 国产熟女xx| 岛国视频午夜一区免费看| 久久久久久亚洲精品国产蜜桃av| 久久香蕉精品热| 88av欧美| 人妻久久中文字幕网| 亚洲成人免费电影在线观看| 久久精品影院6| 人妻久久中文字幕网| 黑人巨大精品欧美一区二区mp4| 嫩草影院精品99| 99久久综合精品五月天人人| 亚洲欧美精品综合一区二区三区| 成人免费观看视频高清| 99热这里只有精品一区 | 1024手机看黄色片| 色婷婷久久久亚洲欧美| 国产成年人精品一区二区| 成人国语在线视频| 亚洲中文字幕日韩| 国产亚洲欧美在线一区二区| 天堂影院成人在线观看| 999久久久国产精品视频| 欧美在线一区亚洲| 一级毛片精品| 日本 av在线| 亚洲五月婷婷丁香| 又黄又粗又硬又大视频| 亚洲成av片中文字幕在线观看| 亚洲成人久久爱视频| 激情在线观看视频在线高清| 国产午夜精品久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 久久人人精品亚洲av| 亚洲五月婷婷丁香| 一本大道久久a久久精品| 99热只有精品国产| 成年女人毛片免费观看观看9| 欧美精品亚洲一区二区| 香蕉丝袜av| 很黄的视频免费| 亚洲午夜精品一区,二区,三区| 成人免费观看视频高清| 日本三级黄在线观看| 麻豆成人av在线观看| 嫩草影院精品99| 99热这里只有精品一区 | 首页视频小说图片口味搜索| 身体一侧抽搐| 午夜亚洲福利在线播放| √禁漫天堂资源中文www| 日本五十路高清| 99国产精品一区二区蜜桃av| 欧美日本视频| 久久久久久久久久黄片| 精品一区二区三区四区五区乱码| 麻豆国产av国片精品| 黑人操中国人逼视频| 深夜精品福利| 精品午夜福利视频在线观看一区| 波多野结衣高清无吗| 日韩av在线大香蕉| 日日摸夜夜添夜夜添小说| 久久婷婷人人爽人人干人人爱| 久久久国产成人免费| 一级a爱视频在线免费观看| 亚洲精品美女久久久久99蜜臀| av片东京热男人的天堂| 丝袜在线中文字幕| 国产成人精品无人区| 50天的宝宝边吃奶边哭怎么回事| 男女做爰动态图高潮gif福利片| 在线观看日韩欧美| 手机成人av网站| 香蕉丝袜av| 中文在线观看免费www的网站 | 琪琪午夜伦伦电影理论片6080| 老司机午夜福利在线观看视频| 国产高清视频在线播放一区| 色播亚洲综合网| 我的亚洲天堂| 久久精品国产清高在天天线| 精品欧美一区二区三区在线| 久久精品91蜜桃| 一个人观看的视频www高清免费观看 | 欧美性长视频在线观看| 波多野结衣高清作品| 在线av久久热| 精品一区二区三区视频在线观看免费| 日韩大码丰满熟妇| 久久精品亚洲精品国产色婷小说| videosex国产| 国产激情久久老熟女| 欧美日本亚洲视频在线播放| 午夜精品在线福利| 曰老女人黄片| 亚洲黑人精品在线| 中文字幕av电影在线播放| 88av欧美| 久久久久久亚洲精品国产蜜桃av| 又紧又爽又黄一区二区| 午夜亚洲福利在线播放| 精品国产国语对白av| 宅男免费午夜| √禁漫天堂资源中文www| 韩国精品一区二区三区| 久久九九热精品免费| av免费在线观看网站| www.999成人在线观看| 日本撒尿小便嘘嘘汇集6| 韩国av一区二区三区四区| 亚洲美女黄片视频| 国产精品日韩av在线免费观看| 极品教师在线免费播放| 亚洲国产欧美一区二区综合| 观看免费一级毛片| 日本成人三级电影网站| 91字幕亚洲| 一个人免费在线观看的高清视频| 国产亚洲av嫩草精品影院| 特大巨黑吊av在线直播 | 90打野战视频偷拍视频| 欧美中文综合在线视频| 婷婷精品国产亚洲av在线| 国产真实乱freesex| 777久久人妻少妇嫩草av网站| 日本精品一区二区三区蜜桃| 国产激情久久老熟女| 精品久久久久久,| 欧美人与性动交α欧美精品济南到| 欧美国产精品va在线观看不卡| 免费观看人在逋| 精品久久久久久久久久免费视频| 国产高清激情床上av| 亚洲欧美精品综合久久99| 悠悠久久av| 精品欧美国产一区二区三| 曰老女人黄片| 国产一区二区三区在线臀色熟女| 亚洲av片天天在线观看| 久久精品影院6| 怎么达到女性高潮| 丰满的人妻完整版| www国产在线视频色| 每晚都被弄得嗷嗷叫到高潮| 中国美女看黄片| 久热爱精品视频在线9| 亚洲狠狠婷婷综合久久图片| 91在线观看av| 国产亚洲欧美在线一区二区| 色播在线永久视频| 19禁男女啪啪无遮挡网站| 成熟少妇高潮喷水视频| 久久狼人影院| 97超级碰碰碰精品色视频在线观看| 午夜福利在线在线| 操出白浆在线播放| 波多野结衣高清无吗| 亚洲色图 男人天堂 中文字幕| 亚洲自偷自拍图片 自拍| 国产精品日韩av在线免费观看| 国产aⅴ精品一区二区三区波| 国产蜜桃级精品一区二区三区| 国产亚洲欧美98| 又紧又爽又黄一区二区| 给我免费播放毛片高清在线观看| 在线十欧美十亚洲十日本专区| 国产亚洲欧美在线一区二区| 人人澡人人妻人| 久久狼人影院| 日韩有码中文字幕| 精品国产亚洲在线| 国产精华一区二区三区| 老司机在亚洲福利影院| 免费高清在线观看日韩| 日韩国内少妇激情av| 亚洲精品粉嫩美女一区| 国产1区2区3区精品| 国产私拍福利视频在线观看| 精品国产乱子伦一区二区三区| 国产精品久久久人人做人人爽| 欧美中文综合在线视频| 精华霜和精华液先用哪个| 国产精品美女特级片免费视频播放器 | 国产亚洲精品久久久久5区| 无限看片的www在线观看| 亚洲av日韩精品久久久久久密| 男女那种视频在线观看| 久久久久久久久久黄片| 久久伊人香网站| 宅男免费午夜| 国产人伦9x9x在线观看| 99精品欧美一区二区三区四区| 国产精品 欧美亚洲| 亚洲激情在线av| 特大巨黑吊av在线直播 | 精品久久蜜臀av无| 此物有八面人人有两片| 最近在线观看免费完整版| 亚洲熟女毛片儿| 日日干狠狠操夜夜爽| 中文字幕高清在线视频| 国产成人精品久久二区二区免费| av有码第一页| 国产成人av激情在线播放| 丰满的人妻完整版| 91麻豆av在线| 亚洲激情在线av| 人妻久久中文字幕网| 午夜福利一区二区在线看| 91麻豆av在线| 国产av一区二区精品久久| 精品欧美一区二区三区在线| 欧美黑人精品巨大| 日韩欧美在线二视频| 99国产精品一区二区蜜桃av| 一本综合久久免费| 久久久久久久精品吃奶| 黄网站色视频无遮挡免费观看| 日韩高清综合在线| 麻豆久久精品国产亚洲av| 久久精品成人免费网站| 999精品在线视频| 91国产中文字幕| 久久久久国产精品人妻aⅴ院| 色精品久久人妻99蜜桃| 在线观看免费视频日本深夜| 国产熟女午夜一区二区三区| 精品一区二区三区视频在线观看免费| 操出白浆在线播放| 999久久久国产精品视频| 国产亚洲欧美98| 日本五十路高清| 性欧美人与动物交配| 久久久国产欧美日韩av| 妹子高潮喷水视频| 午夜福利在线观看吧| 99久久久亚洲精品蜜臀av| 视频在线观看一区二区三区| 侵犯人妻中文字幕一二三四区| 禁无遮挡网站| 久久香蕉精品热| 久99久视频精品免费| 91成人精品电影| 午夜成年电影在线免费观看| 波多野结衣高清作品| 国产精品永久免费网站| 黄色视频不卡| 黄色毛片三级朝国网站| а√天堂www在线а√下载| 免费在线观看视频国产中文字幕亚洲| 麻豆av在线久日| 婷婷精品国产亚洲av| 亚洲av第一区精品v没综合| 侵犯人妻中文字幕一二三四区| 色精品久久人妻99蜜桃| 亚洲国产欧美日韩在线播放| 麻豆国产av国片精品| 午夜免费成人在线视频| 国产单亲对白刺激| 1024香蕉在线观看| 观看免费一级毛片| 一个人免费在线观看的高清视频| 国产熟女xx| 老鸭窝网址在线观看| 又黄又爽又免费观看的视频| 狠狠狠狠99中文字幕| 亚洲三区欧美一区| 亚洲色图 男人天堂 中文字幕| 久久久久亚洲av毛片大全| 在线看三级毛片| 日韩欧美一区视频在线观看| 级片在线观看| videosex国产| 一本综合久久免费| 麻豆av在线久日| 亚洲三区欧美一区| 亚洲欧美激情综合另类| 天天躁狠狠躁夜夜躁狠狠躁| 人妻久久中文字幕网| 一区二区日韩欧美中文字幕| www国产在线视频色| 宅男免费午夜| xxx96com| 亚洲熟妇熟女久久| 99国产精品99久久久久| 黄色视频,在线免费观看| 黄色成人免费大全| 丝袜人妻中文字幕| www.自偷自拍.com| 欧美黄色淫秽网站| 国产91精品成人一区二区三区| 精品第一国产精品| 99国产综合亚洲精品| 老司机深夜福利视频在线观看| 色精品久久人妻99蜜桃| 国产午夜福利久久久久久| 免费在线观看亚洲国产| 18禁黄网站禁片免费观看直播| 熟女少妇亚洲综合色aaa.| 一二三四社区在线视频社区8| 精品高清国产在线一区| 中亚洲国语对白在线视频| 熟女电影av网| 国产视频内射| 国产一区在线观看成人免费| 久久青草综合色| 视频区欧美日本亚洲| 欧美日韩精品网址| 男人的好看免费观看在线视频 | 在线观看一区二区三区| 中文在线观看免费www的网站 | 国产色视频综合| 色哟哟哟哟哟哟| 久久久久久亚洲精品国产蜜桃av| avwww免费| 成年女人毛片免费观看观看9| 在线观看一区二区三区| 国产激情久久老熟女| 亚洲精品国产精品久久久不卡| 国产1区2区3区精品| 一个人观看的视频www高清免费观看 | 亚洲国产精品999在线| 天堂影院成人在线观看| 国产欧美日韩精品亚洲av| 欧美性长视频在线观看| 最近最新免费中文字幕在线| 久久欧美精品欧美久久欧美| 亚洲第一欧美日韩一区二区三区| 国产欧美日韩精品亚洲av| 可以在线观看毛片的网站| 亚洲aⅴ乱码一区二区在线播放 | 国产精品综合久久久久久久免费| 国产亚洲欧美精品永久| avwww免费| 欧美久久黑人一区二区| 午夜福利一区二区在线看| 最近最新中文字幕大全电影3 | 国产成年人精品一区二区| 国产激情久久老熟女| 嫩草影视91久久| 在线观看免费午夜福利视频| 精品少妇一区二区三区视频日本电影| 亚洲欧美一区二区三区黑人| 99国产精品一区二区三区| 女性生殖器流出的白浆| xxx96com| 亚洲成人国产一区在线观看| 黄片播放在线免费| 两个人免费观看高清视频| 中出人妻视频一区二区| 好男人电影高清在线观看| 久久性视频一级片| 国产精品一区二区精品视频观看| 制服丝袜大香蕉在线| 伦理电影免费视频| 中文字幕人妻丝袜一区二区| 制服人妻中文乱码| 国产精品,欧美在线| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品一区av在线观看| 婷婷丁香在线五月| 巨乳人妻的诱惑在线观看| 国产精品久久电影中文字幕| 亚洲av电影在线进入| 日本黄色视频三级网站网址| 婷婷丁香在线五月| 视频区欧美日本亚洲| 一边摸一边做爽爽视频免费| www国产在线视频色| 国产一区二区在线av高清观看| 一二三四在线观看免费中文在| 两人在一起打扑克的视频| 波多野结衣高清作品| 精品福利观看| av片东京热男人的天堂| 亚洲av成人av| 国产精品1区2区在线观看.| 免费在线观看亚洲国产| 欧美黑人欧美精品刺激| 欧美国产精品va在线观看不卡| 国产极品粉嫩免费观看在线| 1024香蕉在线观看| 欧美国产精品va在线观看不卡| 18禁黄网站禁片午夜丰满| 国产私拍福利视频在线观看| 欧美性长视频在线观看| bbb黄色大片| 免费av毛片视频| 国产成+人综合+亚洲专区| 在线观看日韩欧美| 午夜免费观看网址| 又黄又粗又硬又大视频| 国产精品美女特级片免费视频播放器 | 亚洲 欧美一区二区三区| 1024视频免费在线观看| 亚洲国产精品999在线| 精品久久久久久久人妻蜜臀av| 大香蕉久久成人网| avwww免费| av电影中文网址| 操出白浆在线播放|