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

    飛機機體表面聲壓及艙內(nèi)降噪優(yōu)化設(shè)計

    2019-11-30 05:48:16李晨曦徐俊偉趙華勇何立燕
    應(yīng)用聲學(xué) 2019年5期
    關(guān)鍵詞:客艙壁板聲壓級

    李晨曦 胡 瑩 韓 峰 徐俊偉 趙華勇 何立燕

    (中國商飛上海飛機設(shè)計研究院 上海 201210)

    0 引言

    隨著乘坐飛機出行頻率的增加和乘坐時間的延長,乘客對民用客機客艙內(nèi)的噪聲要求越來越高??团搩?nèi)聲環(huán)境的舒適性越來越成為各大民航客機制造商關(guān)注的要點,也成為客機商業(yè)競爭能力的重要指標(biāo)[1?2]。商用客機的艙內(nèi)聲學(xué)設(shè)計也成為聲學(xué)領(lǐng)域研究的熱點之一。飛機飛行時,艙內(nèi)的噪聲是由多個不同的噪聲源以及不同特征頻譜共同產(chǎn)生的[3]。這些噪聲源從產(chǎn)生的位置可以分為外部噪聲和內(nèi)部噪聲。外部噪聲包括附面層噪聲、分離湍流層噪聲、發(fā)動機風(fēng)扇噪聲、發(fā)動機噴流噪聲、發(fā)動機振動噪聲等[4]。外部噪聲主要通過機體表面結(jié)構(gòu)以結(jié)構(gòu)聲或振動的形式傳遞到艙內(nèi)。外部噪聲產(chǎn)生的機理比較復(fù)雜,常常是多種噪聲和振動的綜合效應(yīng)[5]。因此,雖然商用客機的外部噪聲可以用計算流體力學(xué)(Computational fluid dynamics,CFD)或經(jīng)驗公式來模擬,試驗仍然是準(zhǔn)確獲得機體表面聲壓的重要方法。而外部噪聲的主要處理方法是優(yōu)化傳遞路徑,即在機體表面噪聲傳遞到艙內(nèi)的路徑上[6],對機身壁板和降噪聲學(xué)包進行合理的聲學(xué)設(shè)計,以降低最終傳遞到客艙內(nèi)部的噪聲。

    商用客機結(jié)構(gòu)復(fù)雜,對商用飛機對機身壁板和降噪聲學(xué)包進行聲學(xué)研究時,難以使用純理論模型。而商用飛機體積巨大,如果使用有限元算法(Finite element analysis,FEA),中高頻計算量龐大,費時費力。在民用航空和船舶工業(yè)等領(lǐng)域,統(tǒng)計能量分析(Statistical energy analysis,SEA)是解決大型結(jié)構(gòu)中高頻噪聲問題的常用方法之一[7]。該方法基于能量守恒的原理,對結(jié)構(gòu)細節(jié)可模糊處理,計算速度快[8?9],但其精確度取決于子系統(tǒng)的劃分、敏感參數(shù)的獲取以及外部聲源的輸入等因素[10]。

    本文以某型號客機為研究對象,從試驗數(shù)據(jù)分析和聲學(xué)建模兩方面研究機體表面聲壓分布,及其對艙內(nèi)壁板近場輻射聲壓的影響。首先根據(jù)試飛數(shù)據(jù)分析了在巡航狀態(tài)四種工況下的機體表面聲壓分布,然后利用統(tǒng)計能量法建立飛機客艙中后段的聲學(xué)模型,以試飛數(shù)據(jù)作為聲源輸入,研究機體表面聲壓分布對客艙內(nèi)部壁板附近聲壓分布的影響,并在分析基礎(chǔ)上提出優(yōu)化設(shè)計方案。優(yōu)化設(shè)計方案的有效性通過聲學(xué)模型進行驗證。本文建立了實際飛行數(shù)據(jù)與艙內(nèi)聲場仿真結(jié)果之間的關(guān)系,用實際試飛數(shù)據(jù)模擬聲源,增強統(tǒng)計能量法仿真的準(zhǔn)確度;而在試飛數(shù)據(jù)的基礎(chǔ)上用仿真方法研究并優(yōu)化客艙內(nèi)部壁板附近聲場,可以驗證降噪聲學(xué)包設(shè)計方案的有效性,降低優(yōu)化設(shè)計的試飛成本,提高飛機設(shè)計的經(jīng)濟性。

    1 飛機機體表面聲壓分布

    1.1 飛行試驗和數(shù)據(jù)處理

    某型號客機按表1所示的工況在不同巡航高度和巡航速度下進行飛行試驗。為采集機體外表面聲壓數(shù)據(jù),在飛機順航向右側(cè)機體表面布設(shè)聲級計。傳聲器采集到的原始數(shù)據(jù)為時域數(shù)據(jù)。本文通過軟件對原始數(shù)據(jù)進行處理,將其轉(zhuǎn)換為頻域數(shù)據(jù),頻率范圍為50~10000 Hz,利用三角網(wǎng)格差值算法繪制了外表面總聲壓級聲載荷分布云圖和發(fā)動機N1、N2 頻率對應(yīng)的1/3 倍頻程頻帶上的聲壓級聲載荷分布云圖。

    表1 試飛工況Table1 Operation conditions of the flight test

    1.2 線性總聲壓級聲載荷分布

    根據(jù)處理后的頻域數(shù)據(jù),本文繪制了飛行試驗中四種工況對應(yīng)的機體表面總聲壓級分布云圖,如圖1~圖4所示。對比發(fā)現(xiàn),在整個有效測量范圍內(nèi),線性總聲壓級在后應(yīng)急門前方、靠近地板處最大,并從此處向四周蔓延;而在機頭及前機身區(qū)域,線性總聲壓級在登機門后側(cè)到客艙第一舷窗之間較大。這一趨勢與波音在B777-300ER 上的測試數(shù)據(jù)相符[6]。

    圖1 工況1 的線性總聲壓級云圖Fig.1 Distribution of the overall sound pressure on the aircraft surface of operation condition 1

    圖2 工況2 的線性總聲壓級云圖Fig.2 Distribution of the overall sound pressure on the aircraft surface of operation condition 2

    圖3 工況3 的線性總聲壓級云圖Fig.3 Distribution of the overall sound pressure on the aircraft surface of operation condition 3

    圖4 工況4 的線性總聲壓級云圖Fig.4 Distribution of the overall sound pressure on the aircraft surface of operation condition 4

    圖1、圖2顯示在巡航高度相同時(35000 ft),巡航速度從0.72 Ma 增加到0.78 Ma,后應(yīng)急門前方、靠近地板處的線性總聲壓級較大區(qū)域面積增加,高總聲壓級范圍向航向前方蔓延。而在機頭及前機身區(qū)域的線性總聲壓級云圖無明顯變化。在巡航高度均為30000 ft,圖3、圖4呈現(xiàn)類似趨勢,證明了該趨勢的可重復(fù)性。

    當(dāng)巡航速度均為0.72 Ma 時,圖1、圖3中不同巡航高度的云圖無明顯區(qū)別。當(dāng)巡航高度均為0.78 Ma時,圖2和圖4也基本無差別??梢娧埠礁叨葘C體表面聲壓分布無明顯影響。

    1.3 發(fā)動機N1 頻率所在1/3 倍頻程聲壓級聲載荷分布

    發(fā)動機噪聲是客機機體表面聲壓的主要影響因素之一。在不同的巡航飛行工況下,發(fā)動機N1頻率在86~90 Hz 之間,因此可能影響中心頻率為80 Hz和100 Hz兩個1/3倍頻程頻帶。圖5~圖8是發(fā)動機N1 頻率所在的1/3 倍頻程頻帶在工況1 到工況4 對應(yīng)的聲壓級云圖,其中(a)為中心頻率為80 Hz 的1/3 倍頻程頻帶對應(yīng)的云圖,(b)為中心頻率為100 Hz的1/3倍頻程頻帶對應(yīng)的云圖。由云圖可以看出,中心頻率為80 Hz 的云圖與中心頻率為100 Hz 的云圖無明顯差異;在發(fā)動機N1 頻率所在的1/3 倍頻程頻帶內(nèi),機體表面聲壓在后應(yīng)急門前方、靠近地板處最大;巡航速度加大會使此處的聲載荷上升并向四周蔓延;巡航高度對機體表面聲壓無明顯影響。此結(jié)論與總聲壓級云圖結(jié)論基本一致。

    圖5 發(fā)動機N1 頻率在工況1 的聲壓級云圖Fig.5 Distribution of the sound pressure at the N1 frequency of the engine on the aircraft surface of operation condition 1

    圖6 發(fā)動機N1 頻率在工況2 的聲壓級云圖Fig.6 Distribution of the sound pressure at the N1 frequency of the engine on the aircraft surface of operation condition 2

    圖7 發(fā)動機N1 頻率在工況3 的聲壓級云圖Fig.7 Distribution of the sound pressure at the N1 frequency of the engine on the aircraft surface of operation condition 3

    圖8 發(fā)動機N1 頻率在工況4 的聲壓級云圖Fig.8 Distribution of the sound pressure at the N1 frequency of the engine on the aircraft surface of operation condition 4

    1.4 發(fā)動機N2 頻率所在1/3 倍頻程聲壓級聲載荷分布

    在不同的巡航飛行工況下,發(fā)動機N2 頻率在273±10 Hz,影響中心頻率為250 Hz 的1/3 倍頻程頻帶。圖9~圖12是發(fā)動機N2頻率所在的1/3倍頻程頻帶在工況1 到工況4 對應(yīng)的聲壓級云圖。由云圖可以看出,在發(fā)動機N2頻率所在的1/3倍頻程頻帶內(nèi),機體表面聲壓在后應(yīng)急門前方、靠近地板處最大;巡航速度加大會使此處的聲載荷上升并向四周蔓延;巡航高度對機體表面聲壓無明顯影響。此結(jié)論與總聲壓級云圖、發(fā)動機N1頻率所在1/3倍頻程云圖的結(jié)論基本一致,說明發(fā)動機在N2 頻率上的噪聲對機體表面聲壓分布趨勢的影響與N1 頻率的影響基本一致。

    圖9 發(fā)動機N2 頻率在工況1 的聲壓級云圖,無計權(quán),中心頻率為250 HzFig.9 Distribution of the sound pressure on the aircraft surface in the frequency band which center frequency is 250 Hz and includes the N2 frequency of the engine of operation condition 1

    圖11 發(fā)動機N2 頻率在工況3 的聲壓級云圖,無計權(quán),中心頻率為250 HzFig.11 Distribution of the sound pressure on the aircraft surface in the frequency band which center frequency is 250 Hz and includes the N2 frequency of the engine of operation condition 3

    圖12 發(fā)動機N2 頻率在工況4 的聲壓級云圖,無計權(quán),中心頻率為250 HzFig.12 Distribution of the sound pressure on the aircraft surface in the frequency band which center frequency is 250 Hz and includes the N2 frequency of the engine of operation condition 4

    2 飛機客艙中后段聲學(xué)建模

    2.1 聲學(xué)模型

    飛行試驗數(shù)據(jù)表明機體表面聲壓在后應(yīng)急門前方、靠近地板處最大,并從此處向四周蔓延。因此,為滿足客艙內(nèi)的舒適性,應(yīng)對該型號客機客艙中后段進行聲學(xué)設(shè)計。為研究機體表面聲壓的分布趨勢對客艙內(nèi)部噪聲分布的影響,本文用統(tǒng)計能量法建立了飛機客艙中后段的聲學(xué)模型如圖13所示。

    圖13 客艙中后段聲學(xué)模型Fig.13 SEA model of the middle and after fuselage

    聲學(xué)模型參考飛行試驗表面聲壓測點布置位置進行飛機蒙皮部分的子系統(tǒng)劃分,并根據(jù)統(tǒng)計能量法的算法原理,盡量劃分較大的壁板、聲腔等子系統(tǒng),以增加模型的準(zhǔn)確性。本聲學(xué)模型中飛機蒙皮的肋板結(jié)構(gòu)、客艙玻璃的參數(shù)均為試驗測量結(jié)果,具體參數(shù)見表2。

    表2 材料參數(shù)Table2 Parameters of the materials

    本模型主要研究測量得到的外場噪聲對客艙壁板近場聲輻射的影響以及聲學(xué)降噪包的降噪效果,因此未考慮行李架、客艙座椅、空調(diào)噪聲等對客艙聲場的影響。統(tǒng)計能量法中,聲空腔的聲壓可表示為

    圖13的模型中,飛行試驗測得的機體表面聲壓數(shù)據(jù)作為聲源從機艙壁板外表面輸入模型。由于飛行試驗測得的四種工況下的機體表面聲壓分布特征具有高度一致性,在仿真建模中僅選取工況1 對應(yīng)的聲載荷作為輸入聲源。為便于研究優(yōu)化設(shè)計方案,對機艙壁板上的降噪聲學(xué)包安裝區(qū)域進行編號,如圖13所示。

    為研究外場噪聲對客艙壁板近場聲輻射,在客艙內(nèi)壁板附近設(shè)置半無限流體(Semi-infinite fluid,SIF)測點。這些SIF測點監(jiān)測對應(yīng)壁板的近場輻射聲壓,只與對應(yīng)壁板子系統(tǒng)連接,未考慮客艙聲場混響和其他壁板、舷窗等子系統(tǒng)對該測點聲壓的影響。因此,SIF測點測量的聲壓值與客艙內(nèi)的混響聲壓有一定差異,兩者間不能直接對比。監(jiān)測客艙內(nèi)壁板附近SIF 測點的聲壓級可以排除其他子系統(tǒng)對特定壁板的聲輻射影響,直接研究聲學(xué)降噪包鋪設(shè)方案在壁板局部的有效性,有助于評估和優(yōu)化聲學(xué)降噪包鋪設(shè)方案。

    統(tǒng)計能量法中頻帶內(nèi)的模態(tài)數(shù)高于5才能保證統(tǒng)計能量算法的準(zhǔn)確性。本文統(tǒng)計模型中每個子系統(tǒng)各個頻帶內(nèi)的模態(tài)數(shù)后發(fā)現(xiàn),只有在315 Hz以上頻段才能保證所有子系統(tǒng)模態(tài)數(shù)大于5。因此,仿真結(jié)果的有效頻帶范圍是315~10000 Hz。仿真結(jié)果中的線性總聲壓級是315~10000 Hz 頻段的總聲壓級,與飛行試驗測得的機體表面總聲壓級頻率范圍不同,仿真結(jié)果也無法與中低頻段內(nèi)的實測數(shù)據(jù)進行對比。由于本文研究聲學(xué)降噪包鋪設(shè)的優(yōu)化方案,而以纖維材料為主的降噪聲學(xué)包的有效范圍也在中高頻段,統(tǒng)計能量法模型的有效頻率范圍是315~10000 Hz,可以滿足研究要求。如果需要用仿真方法研究中低頻段的艙內(nèi)噪聲,建議使用有限元法,但有限元仿真超出本文的研究范圍。

    2.2 壁板仿真結(jié)果

    以在工況1測量的機體表面聲壓作為聲源輸入的情況下,本文應(yīng)用建立的聲學(xué)模型計算出僅考慮蒙皮結(jié)構(gòu)、無降噪聲學(xué)包和內(nèi)飾板的客艙內(nèi)壁板近場輻射聲壓級云圖,如圖14 所示。由圖14 可以看出,試驗測量得到的機體表面聲壓通過肋板蒙皮結(jié)構(gòu)傳入客艙后,在蒙皮壁板附近形成的近場聲壓的分布趨勢是在機艙后部應(yīng)急門到球面框區(qū)和應(yīng)急門下方靠近地板處聲壓最大,即在8 區(qū)和9 區(qū)聲壓最大;高聲壓區(qū)在3 區(qū)中后方和4 區(qū)下方也有分布。對比圖1、圖14 發(fā)現(xiàn),機體表面聲壓在6 區(qū)和7 區(qū)下方較大。當(dāng)機體表面聲壓通過肋板結(jié)構(gòu)傳入機艙后,艙內(nèi)壁板附近的近場聲壓分布與機體表面聲壓分布趨勢類似。但艙內(nèi)壁板近場聲壓的高聲壓區(qū)除了6 區(qū)和7 區(qū)下方以外,還增加了8 區(qū)、9 區(qū)、3 區(qū)中部以及4區(qū)下方。

    圖14 工況1 的客艙中后段的聲學(xué)模型線性總聲壓級仿真云圖僅考慮蒙皮結(jié)構(gòu),無降噪聲學(xué)包和內(nèi)飾板Fig.14 Distribution of the simulated overall sound pressure on the aircraft surface in the mid-after cabin of operation condition 1,skin included only

    圖15 3 區(qū)、4 區(qū)、8 區(qū)、9 區(qū)的輸入聲功率Fig.15 The simulated results of the power inputs on the Areas 3,4,8,and 9

    為分析3 區(qū)、4 區(qū)、8 區(qū)、9 區(qū)的高聲壓區(qū)產(chǎn)生的原因,本文用聲學(xué)模型分析了對應(yīng)區(qū)域壁板的輸入聲功率,如圖15所示。圖15中,3區(qū)的輸入聲功率圖中,聲源貢獻量曲線在1250 Hz 以上頻段內(nèi)幾乎與總輸入聲功率曲線重合,而在1250 Hz以下頻段內(nèi)4區(qū)壁板的貢獻量稍高于聲源貢獻量??梢妼? 區(qū)的壁板結(jié)構(gòu)而言,在1250 Hz 以上頻段主要受聲源影響,在1250 Hz 以下頻段內(nèi)主要受4 區(qū)壁板傳遞的能量影響。而4 區(qū)、8 區(qū)、9 區(qū)的輸入聲功率圖中,總輸入聲功率與聲源貢獻量重合程度高,顯示出主要受聲源影響的特性。因此,3 區(qū)、4區(qū)、8 區(qū)、9 區(qū)中噪聲的主要傳遞路徑是從機體表面通過壁板結(jié)構(gòu)傳遞到艙內(nèi)??梢灶A(yù)見,在3 區(qū)、4 區(qū)、8 區(qū)、9區(qū)使用聲學(xué)降噪包可以消耗噪聲主要傳遞路徑上的聲能量,降低傳遞到艙內(nèi)的聲能量,從而降低艙內(nèi)聲壓級。

    3 艙內(nèi)降噪優(yōu)化設(shè)計

    3.1 設(shè)計方案

    由飛行試驗測量得到的機體表面聲壓和僅考慮蒙皮結(jié)構(gòu)的艙內(nèi)壁板近場聲壓分布仿真可以發(fā)現(xiàn),機體表面聲壓傳入機艙后,在蒙皮附近形成數(shù)個高聲壓區(qū),包括6 區(qū)和7 區(qū)下方,8 區(qū)、9 區(qū)和3 區(qū)中部以及4 區(qū)下方。因此,可在降噪聲學(xué)包的原始設(shè)計方案基礎(chǔ)上,在以上高聲壓區(qū)用高密度隔音棉代替普通隔音棉,以達到在厚度相同的條件下優(yōu)化艙內(nèi)聲場的目的。兩種隔音棉的材料參數(shù)見表3。

    表3 隔音棉材料參數(shù)Table3 Parameters of the fibers used in the thermal acoustic insulation blankets

    該型號飛機的降噪聲學(xué)包原始設(shè)計方案是在客艙壁板和天花板區(qū)鋪設(shè)普通隔音棉和內(nèi)飾板,其中隔音棉使用表3中的普通隔音棉。本文提出的優(yōu)化設(shè)計方案中,優(yōu)化方案1只在5 區(qū)到9區(qū)使用高密度隔音棉,優(yōu)化方案2在3 區(qū)到9區(qū)全部使用高密度隔音棉,如圖16所示。三種方案的詳細對比見表4。

    圖16 優(yōu)化設(shè)計方案示意圖Fig.16 Proposed designs for the thermal acoustic insulation blankets.

    表4 統(tǒng)計能量法仿真的降噪聲學(xué)包設(shè)計方案Table4 Designs of the thermal acoustic insulation blankets used in the SEA model

    3.2 設(shè)計方案仿真結(jié)果

    圖17~圖19 顯示使用降噪聲學(xué)包的原始設(shè)計方案、優(yōu)化設(shè)計方案1 和優(yōu)化設(shè)計方案2 后的艙內(nèi)壁板近場聲壓分布仿真結(jié)果。對比圖14 和圖17 發(fā)現(xiàn),在客艙壁板上按照原始設(shè)計方案添加聲學(xué)降噪包和內(nèi)飾板后,機艙壁板近場聲壓整體有所降低,但仍在6 區(qū)和7 區(qū)下方,8 區(qū)、9 區(qū)和3 區(qū)中部以及4 區(qū)下方有較高的聲壓分布。

    圖17 工況1 的客艙中后段的聲學(xué)模型線性總聲壓級仿真云圖,使用原始降噪聲學(xué)包設(shè)計Fig.17 Distribution of the simulated overall sound pressure on the aircraft surface in the midafter cabin of operation condition 1,using the original design of the thermal acoustic blankets

    圖18 工況1 的客艙中后段的聲學(xué)模型線性總聲壓級仿真云圖,優(yōu)化設(shè)計方案1Fig.18 Distribution of the simulated overall sound pressure on the aircraft surface in the midafter cabin of operation condition 1,using the proposed design 1 of the thermal acoustic blankets

    圖19 工況1 的客艙中后段的聲學(xué)模型線性總聲壓級仿真云圖,優(yōu)化設(shè)計方案2Fig.19 Distribution of the simulated overall sound pressure on the aircraft surface in the midafter cabin of operation condition 1,using the proposed design 2 of the thermal acoustic blankets

    對比圖17 和圖18 發(fā)現(xiàn),使用優(yōu)化方案1、在5 區(qū)~9 區(qū)用高密度隔音棉替換普通隔音棉后,5 區(qū)~9 區(qū)的壁板近場聲壓降低,但3 區(qū)和4 區(qū)的高聲壓區(qū)域依然存在。因此可以將高密度隔音棉區(qū)從5 區(qū)~9 區(qū)向3 區(qū)和4 區(qū)擴展,即為優(yōu)化方案2。對比優(yōu)化方案2 和優(yōu)化方案1 的仿真結(jié)果發(fā)現(xiàn),在3 區(qū)和4 區(qū)用高密度隔音棉替換低密度隔音棉不僅能降低此處的壁板近場聲壓,相鄰的1 區(qū)、2 區(qū)、5 區(qū)、6 區(qū)、10 區(qū)聲壓也有所降低,而這種降低效果與壁板的耦合關(guān)系和聲能量傳遞路徑有關(guān)。

    3.3 艙內(nèi)聲壓級仿真結(jié)果

    由艙內(nèi)壁板的近場輻射聲壓仿真結(jié)果可知,在降噪聲學(xué)包的優(yōu)化設(shè)計方案中,用高密度隔音棉替換低密度隔音棉能降低對應(yīng)的壁板近場聲壓和鄰近區(qū)域壁板的近場輻射聲壓。而本文仿真模型中的近場輻射聲壓只與對應(yīng)壁板有關(guān),可以驗證降噪聲學(xué)包的有效性,未考慮客艙聲場混響和其他壁板、舷窗等的影響。圖20 是仿真模型計算出的客艙內(nèi)聲壓級,體現(xiàn)了客艙內(nèi)的混響以及壁板、舷窗的整體影響。其聲源激勵如圖13 所示,分為只考慮蒙皮和蒙皮+降噪聲學(xué)包等構(gòu)型,降噪聲學(xué)包設(shè)計方案見表4。

    由圖20可以看出,僅考慮蒙皮結(jié)構(gòu)的艙內(nèi)聲壓級較高,說明蒙皮結(jié)構(gòu)無法有效降低艙內(nèi)噪聲,有必要使用降噪聲學(xué)包和內(nèi)飾板。比較圖20 中的實線和虛線可以看出,添加原始設(shè)計方案的降噪聲學(xué)包和內(nèi)飾板能有效降低艙內(nèi)聲壓級,在2500 Hz 最大可降低10.5 dB。三種隔音棉設(shè)計方案的降噪效果在1600 Hz 以上基本一致,而在315~1600 Hz 頻段上,使用三種隔音棉設(shè)計方案得到的艙內(nèi)聲壓級依次降低,證明增加高密度隔音棉的鋪設(shè)范圍可以降低艙內(nèi)聲壓級。此結(jié)論與艙內(nèi)壁板近場輻射聲壓仿真結(jié)果所得到的結(jié)論一致,也符合2.2 小節(jié)中對噪聲傳遞路徑的分析和對降噪方法的預(yù)期。因此,在實際工程設(shè)計中,應(yīng)綜合考慮壁板近場輻射聲壓、艙內(nèi)聲壓、設(shè)計方案重量等因素對設(shè)計方案進行取舍和優(yōu)化,以達到在保證客機經(jīng)濟性的同時優(yōu)化客艙聲環(huán)境的目的。

    為驗證模型,本文將仿真得到的A 計權(quán)艙內(nèi)聲壓級與在艙內(nèi)不同位置測量得到的A 計權(quán)聲壓級進行對比,如圖21 所示。圖中細線代表在艙內(nèi)不同測點位置測量得到的A 計權(quán)聲壓級,黑色粗實線代表仿真得到的A 計權(quán)艙內(nèi)聲壓級。由公式(1)可知,仿真得到的聲壓級是聲空腔內(nèi)對體積求平均得到的聲壓級。由圖21 可以看出,仿真數(shù)據(jù)與實測數(shù)據(jù)呈現(xiàn)相同的趨勢,但仿真數(shù)據(jù)在整個有效頻段內(nèi)比實測數(shù)據(jù)高1~8 dB。這是由于本文使用的仿真模型是一個簡化的客機艙段模型,未考慮行李架、客艙座椅等對客艙聲場的影響。而實測數(shù)據(jù)則包括所有客艙組件對客艙聲空腔的影響。圖22是實際測量得到的艙內(nèi)座椅吸聲系數(shù)。對比圖21和圖22發(fā)現(xiàn),艙內(nèi)座椅的吸聲系數(shù)在315 Hz 以上能達到0.6 以上。由于艙內(nèi)座椅鋪設(shè)面積大,總表面積大,可以推斷仿真數(shù)據(jù)與實測數(shù)據(jù)之間存在差距主要是因為仿真中沒有考慮客艙的艙內(nèi)座椅吸聲等因素。

    圖20 工況1 的客艙中后段的聲學(xué)模型計算出的艙內(nèi)聲壓級Fig.20 The simulated sound pressure levels in the mid-after cabin under the excitation of the measured sound pressure levels on the outer surface of the aircraft with operation condition 1

    圖21 工況1 的客聲學(xué)模型仿真結(jié)果與在不同艙內(nèi)位置實測的A 計權(quán)聲壓級對比,使用原始降噪聲學(xué)包設(shè)計Fig.21 The comparison of the simulated sound pressure levels and the measured sound pressure levels in the mid-after cabin under the excitation of the measured sound pressure levels with operation condition 1

    圖22 實際測量的艙內(nèi)座椅吸聲系數(shù)Fig.22 The measured sound absorption coefficient of the passenger seats in the cabin

    4 結(jié)論與討論

    本文從試驗數(shù)據(jù)分析和聲學(xué)建模兩方面研究機體表面聲壓及其對艙內(nèi)壁板近場輻射聲壓的影響,提出了優(yōu)化設(shè)計方案,并用聲學(xué)模型驗證了優(yōu)化設(shè)計方案的有效性。其中,聲學(xué)建模用試驗數(shù)據(jù)作為聲源激勵,其他結(jié)構(gòu)參數(shù)均采用該型客機的實際參數(shù),材料參數(shù)均來自實際測量。

    通過分析飛行試驗獲得的線性總聲壓級云圖、發(fā)動機N1、N2 頻率對應(yīng)的1/3 倍頻程云圖發(fā)現(xiàn),總聲壓級和發(fā)動機N1、N2 頻率對應(yīng)的1/3 倍頻程聲壓在機體表面的聲載荷分布趨勢基本一致。該型號客機在巡航狀態(tài)下的機體表面聲壓分布特性可總結(jié)為

    (1)在整個測量范圍內(nèi),聲壓級在后應(yīng)急門前方、靠近地板處最大,并從此處向四周蔓延;而在機頭及前機身區(qū)域,聲壓級在登機門后側(cè)到客艙第一舷窗之間較大。

    (2)在相同巡航高度,隨著巡航速度的增加,后應(yīng)急門前方、靠近地板處的聲壓級較大區(qū)域面積增加,高聲壓級范圍向航向前方蔓延;而機頭及前機身區(qū)域的聲壓級云圖無明顯變化。

    (3)在相同巡航速度,不同巡航高度對機體表面聲壓分布無明顯影響。

    (4)發(fā)動機N1、N2 頻率對應(yīng)的1/3 倍頻程云圖中,機體表面聲壓的分布趨勢無明顯差異。

    通過建立聲學(xué)模型、分區(qū)加載飛行試驗測量得到的機體表面聲壓作為聲源進行統(tǒng)計能量法聲學(xué)仿真,并結(jié)合仿真結(jié)果提出降噪聲學(xué)包的優(yōu)化方案1 和優(yōu)化方案2。由仿真模型得到的壁板近場聲輻射仿真結(jié)果和艙內(nèi)聲壓仿真結(jié)果一致表明:

    (1)機體表面聲壓通過肋板結(jié)構(gòu)傳入機艙后,艙內(nèi)壁板附近的近場聲壓分布與機體表面聲壓分布趨勢類似;

    (2)僅蒙皮結(jié)構(gòu)無法有效降低客艙噪聲,有必要使用降噪聲學(xué)包和內(nèi)飾板;

    (3)添加原始設(shè)計方案的降噪聲學(xué)包和內(nèi)飾板能有效降低艙內(nèi)聲壓級,但仍需要在特定區(qū)域?qū)υ荚O(shè)計方案進行優(yōu)化設(shè)計;

    (4)本文降噪聲學(xué)包優(yōu)化方案2 的降噪效果好于優(yōu)化方案1,在實際工程設(shè)計中,應(yīng)綜合考慮設(shè)計方案聲學(xué)仿真結(jié)果和重量等因素對設(shè)計方案進行取舍和優(yōu)化。

    本文用實際試飛數(shù)據(jù)模擬聲源進行統(tǒng)計能量法聲學(xué)仿真,增強了統(tǒng)計能量法仿真的準(zhǔn)確度;而在試飛數(shù)據(jù)的基礎(chǔ)上用仿真方法研究客艙內(nèi)部壁板附近聲場和艙內(nèi)聲壓,以仿真結(jié)果為基礎(chǔ)提出了優(yōu)化設(shè)計方案,并用仿真建模驗證了降噪聲學(xué)包設(shè)計方案的有效性。這種試驗數(shù)據(jù)分析與聲學(xué)仿真相結(jié)合的工程設(shè)計方案可以有效降低優(yōu)化設(shè)計的試飛成本,提高飛機設(shè)計的經(jīng)濟性,值得繼續(xù)研究和推廣。

    猜你喜歡
    客艙壁板聲壓級
    機器噪聲平均聲壓級計算方法差異性實證研究
    電動工具(2024年1期)2024-02-29 01:40:24
    客艙內(nèi)的空氣從何而來?
    國內(nèi)主流機型客艙聲品質(zhì)表現(xiàn)分析
    某大型飛機復(fù)合材料壁板工藝仿真及驗證技術(shù)
    一種計算消聲室聲壓級的新方法
    民用飛機客艙干燥系統(tǒng)的性能分析
    淺談航空安全員應(yīng)變能力在客艙擒敵術(shù)教學(xué)中的培養(yǎng)
    全新DXR mkll有源揚聲器
    演藝科技(2019年4期)2019-03-30 03:21:46
    航天器復(fù)雜整體壁板加工精度控制
    機翼下壁板裂紋擴展分析
    智富時代(2018年5期)2018-07-18 17:52:04
    日韩国内少妇激情av| 在线天堂最新版资源| 亚洲第一区二区三区不卡| 97人妻精品一区二区三区麻豆| 高清av免费在线| 色播亚洲综合网| 久久精品久久久久久久性| 久久精品熟女亚洲av麻豆精品| 国产黄片美女视频| 少妇丰满av| 久久综合国产亚洲精品| 少妇丰满av| 校园人妻丝袜中文字幕| 欧美成人精品欧美一级黄| 日韩一区二区视频免费看| 免费观看av网站的网址| 18禁裸乳无遮挡动漫免费视频 | 中文字幕制服av| 国产精品人妻久久久久久| 午夜福利网站1000一区二区三区| 夫妻午夜视频| 亚洲精品日本国产第一区| 国产 一区 欧美 日韩| 白带黄色成豆腐渣| 亚洲欧美日韩无卡精品| 欧美zozozo另类| 久热久热在线精品观看| 黄色配什么色好看| 亚洲国产av新网站| 我要看日韩黄色一级片| 一区二区三区四区激情视频| 校园人妻丝袜中文字幕| 一级毛片我不卡| 国语对白做爰xxxⅹ性视频网站| 免费观看的影片在线观看| 少妇熟女欧美另类| 久久精品国产亚洲av天美| 99热这里只有是精品50| 特级一级黄色大片| 婷婷色综合大香蕉| 亚洲精品影视一区二区三区av| 中文字幕制服av| 九草在线视频观看| 国产伦精品一区二区三区视频9| 亚洲色图综合在线观看| 国产乱人偷精品视频| 亚洲经典国产精华液单| 爱豆传媒免费全集在线观看| 国产人妻一区二区三区在| 国产免费又黄又爽又色| 亚洲精品亚洲一区二区| 久久精品综合一区二区三区| 色视频www国产| 亚洲人成网站在线观看播放| 免费高清在线观看视频在线观看| 免费观看性生交大片5| 超碰97精品在线观看| 国产成人精品一,二区| 欧美97在线视频| 大又大粗又爽又黄少妇毛片口| 久久精品国产a三级三级三级| 日日啪夜夜爽| 国产伦精品一区二区三区四那| 亚洲在久久综合| 高清午夜精品一区二区三区| 亚洲真实伦在线观看| 日韩不卡一区二区三区视频在线| 国产精品国产三级国产av玫瑰| 黄色一级大片看看| 久久久久国产网址| 久久女婷五月综合色啪小说 | 菩萨蛮人人尽说江南好唐韦庄| 大码成人一级视频| 搡老乐熟女国产| 久久人人爽av亚洲精品天堂 | 欧美成人精品欧美一级黄| 免费av观看视频| 成人免费观看视频高清| 22中文网久久字幕| 日韩精品有码人妻一区| 一级片'在线观看视频| 日本av手机在线免费观看| 亚洲精品日本国产第一区| 精品亚洲乱码少妇综合久久| 1000部很黄的大片| 少妇裸体淫交视频免费看高清| 亚洲精品第二区| 亚洲精华国产精华液的使用体验| 看非洲黑人一级黄片| 国产免费视频播放在线视频| 一级爰片在线观看| 神马国产精品三级电影在线观看| 国产高清有码在线观看视频| 国产高清三级在线| 亚洲av中文av极速乱| av在线天堂中文字幕| 亚洲精品456在线播放app| 成人高潮视频无遮挡免费网站| 日韩伦理黄色片| 欧美成人a在线观看| 国产乱人偷精品视频| 亚洲第一区二区三区不卡| 国产国拍精品亚洲av在线观看| 国产 精品1| 免费看不卡的av| 街头女战士在线观看网站| 色视频在线一区二区三区| 极品教师在线视频| av.在线天堂| 久久人人爽人人爽人人片va| 亚洲精品456在线播放app| h日本视频在线播放| 亚洲综合色惰| 亚洲欧美中文字幕日韩二区| 99re6热这里在线精品视频| 欧美bdsm另类| 在现免费观看毛片| 亚洲,欧美,日韩| 日韩伦理黄色片| 一级毛片我不卡| 国产淫语在线视频| 深夜a级毛片| av卡一久久| 真实男女啪啪啪动态图| 久久99精品国语久久久| 国模一区二区三区四区视频| 欧美97在线视频| 纵有疾风起免费观看全集完整版| 国产精品久久久久久精品电影| 嫩草影院精品99| 亚洲成人久久爱视频| 国产亚洲av片在线观看秒播厂| 国产成人一区二区在线| 亚洲熟女精品中文字幕| 大香蕉97超碰在线| 色播亚洲综合网| 美女高潮的动态| 国产精品一区二区三区四区免费观看| 久久久久久久久大av| 亚洲精品日本国产第一区| 人人妻人人看人人澡| 亚洲精品日韩av片在线观看| 色婷婷久久久亚洲欧美| 成人高潮视频无遮挡免费网站| 高清日韩中文字幕在线| 国产 精品1| 久久精品国产亚洲av涩爱| 少妇人妻 视频| 热99国产精品久久久久久7| 亚洲人成网站在线观看播放| 精品一区二区免费观看| 一级黄片播放器| 黄色一级大片看看| 综合色av麻豆| 国产大屁股一区二区在线视频| 天天躁夜夜躁狠狠久久av| 少妇裸体淫交视频免费看高清| 国产在线一区二区三区精| 欧美精品人与动牲交sv欧美| 91久久精品国产一区二区成人| 在线观看人妻少妇| 国产 一区精品| av国产久精品久网站免费入址| 国产成人精品一,二区| 蜜臀久久99精品久久宅男| 国产成人免费无遮挡视频| videos熟女内射| 日韩人妻高清精品专区| 王馨瑶露胸无遮挡在线观看| kizo精华| 国产成年人精品一区二区| 久久精品国产亚洲网站| 精品熟女少妇av免费看| 嫩草影院新地址| 日本一本二区三区精品| 97人妻精品一区二区三区麻豆| 日韩免费高清中文字幕av| 久热这里只有精品99| 日本与韩国留学比较| 久久午夜福利片| 成人一区二区视频在线观看| 亚洲怡红院男人天堂| 乱系列少妇在线播放| 免费大片黄手机在线观看| 亚洲av免费在线观看| 青春草国产在线视频| 黄色日韩在线| 偷拍熟女少妇极品色| 九九爱精品视频在线观看| 人妻夜夜爽99麻豆av| 色视频在线一区二区三区| 真实男女啪啪啪动态图| 久久久久精品久久久久真实原创| 亚洲欧美日韩另类电影网站 | 美女脱内裤让男人舔精品视频| 国产一区二区亚洲精品在线观看| 久久久久久久久久久丰满| 午夜福利高清视频| 日韩欧美精品免费久久| 高清视频免费观看一区二区| 少妇的逼好多水| 免费看日本二区| 永久网站在线| 少妇猛男粗大的猛烈进出视频 | 国产精品女同一区二区软件| 精品久久久噜噜| 亚洲熟女精品中文字幕| 少妇人妻精品综合一区二区| 亚洲电影在线观看av| 久久影院123| 久久精品国产a三级三级三级| 亚洲精品乱码久久久久久按摩| 亚洲高清免费不卡视频| 欧美+日韩+精品| 久久99热这里只有精品18| 麻豆成人午夜福利视频| 国产午夜精品一二区理论片| 亚洲欧美一区二区三区黑人 | 国产一级毛片在线| 国产精品99久久久久久久久| 国产精品无大码| 亚洲精品久久午夜乱码| 免费看av在线观看网站| 日韩,欧美,国产一区二区三区| 91aial.com中文字幕在线观看| 青青草视频在线视频观看| 亚洲成人一二三区av| 免费看光身美女| 欧美亚洲 丝袜 人妻 在线| 日本av手机在线免费观看| 成人美女网站在线观看视频| 97人妻精品一区二区三区麻豆| 国产片特级美女逼逼视频| 白带黄色成豆腐渣| 校园人妻丝袜中文字幕| 亚洲精品一二三| 国产成人免费观看mmmm| 欧美成人精品欧美一级黄| 国产爱豆传媒在线观看| 禁无遮挡网站| 亚洲成人av在线免费| 偷拍熟女少妇极品色| 午夜精品一区二区三区免费看| 日本色播在线视频| 国产美女午夜福利| 精品国产露脸久久av麻豆| 视频中文字幕在线观看| 精品久久久久久久人妻蜜臀av| 少妇人妻久久综合中文| 亚洲av不卡在线观看| 国内少妇人妻偷人精品xxx网站| 国产免费视频播放在线视频| 国产探花极品一区二区| 噜噜噜噜噜久久久久久91| 最近最新中文字幕大全电影3| 日本与韩国留学比较| 好男人在线观看高清免费视频| 大香蕉久久网| av在线播放精品| 久久久精品免费免费高清| 亚洲国产av新网站| 狠狠精品人妻久久久久久综合| 日韩人妻高清精品专区| 中文乱码字字幕精品一区二区三区| 大香蕉久久网| 91在线精品国自产拍蜜月| 国产欧美日韩精品一区二区| 日日啪夜夜爽| 亚洲熟女精品中文字幕| 国产乱来视频区| 永久网站在线| 午夜福利视频1000在线观看| 成人亚洲欧美一区二区av| 国产精品女同一区二区软件| 欧美日本视频| 国产黄a三级三级三级人| 国产精品久久久久久精品古装| 国产成人精品福利久久| 国产精品国产三级专区第一集| 国产精品女同一区二区软件| 日韩一区二区三区影片| 91狼人影院| 日韩精品有码人妻一区| 色播亚洲综合网| 丰满少妇做爰视频| 色哟哟·www| 免费黄色在线免费观看| 国产探花极品一区二区| 国产伦精品一区二区三区视频9| 99久久精品国产国产毛片| 久久久精品免费免费高清| 国产一区二区三区av在线| 国模一区二区三区四区视频| 看十八女毛片水多多多| 一级毛片久久久久久久久女| 色5月婷婷丁香| 日韩av免费高清视频| 又爽又黄a免费视频| 国产一区二区在线观看日韩| 18禁动态无遮挡网站| 免费看a级黄色片| 乱系列少妇在线播放| 亚洲第一区二区三区不卡| 精品久久久久久久末码| 白带黄色成豆腐渣| 久久韩国三级中文字幕| 国产毛片在线视频| 成人漫画全彩无遮挡| 午夜视频国产福利| 一级片'在线观看视频| 国产综合懂色| 七月丁香在线播放| 久久精品国产亚洲av天美| av国产免费在线观看| 日本黄色片子视频| 汤姆久久久久久久影院中文字幕| 秋霞在线观看毛片| 国产亚洲5aaaaa淫片| 丝袜脚勾引网站| 日本免费在线观看一区| 插阴视频在线观看视频| 夜夜看夜夜爽夜夜摸| 六月丁香七月| 在线a可以看的网站| 亚洲色图综合在线观看| 国产成人freesex在线| 欧美性感艳星| 精品少妇黑人巨大在线播放| 一级毛片久久久久久久久女| 久久精品国产亚洲av涩爱| av在线天堂中文字幕| 小蜜桃在线观看免费完整版高清| 国国产精品蜜臀av免费| 国产伦精品一区二区三区四那| 又粗又硬又长又爽又黄的视频| 最近最新中文字幕大全电影3| 久久久成人免费电影| 亚洲美女搞黄在线观看| 精品久久久久久久久亚洲| 欧美3d第一页| 深夜a级毛片| 性色av一级| 男人舔奶头视频| 亚洲人与动物交配视频| 神马国产精品三级电影在线观看| 在线免费观看不下载黄p国产| 纵有疾风起免费观看全集完整版| 国产人妻一区二区三区在| 边亲边吃奶的免费视频| 中文乱码字字幕精品一区二区三区| 午夜福利在线观看免费完整高清在| 一区二区三区乱码不卡18| 国产黄片美女视频| 欧美成人a在线观看| 亚洲色图av天堂| 国产精品精品国产色婷婷| 99久久人妻综合| 国产精品偷伦视频观看了| 亚洲欧美日韩另类电影网站 | 少妇高潮的动态图| 毛片一级片免费看久久久久| 黄色日韩在线| 午夜福利视频1000在线观看| 国产免费一级a男人的天堂| 国产精品熟女久久久久浪| 肉色欧美久久久久久久蜜桃 | 天天躁夜夜躁狠狠久久av| 熟女人妻精品中文字幕| 欧美日韩综合久久久久久| 美女国产视频在线观看| freevideosex欧美| 国产精品精品国产色婷婷| 久久久久久久久久久丰满| 亚洲欧美精品专区久久| 国产精品99久久99久久久不卡 | 久久精品久久久久久久性| 黄色日韩在线| 综合色av麻豆| 日日啪夜夜爽| 亚洲成人中文字幕在线播放| 日韩 亚洲 欧美在线| 色综合色国产| 精品午夜福利在线看| 国产乱来视频区| 亚洲av日韩在线播放| 好男人在线观看高清免费视频| 国产高潮美女av| 亚洲欧美一区二区三区国产| 国产黄频视频在线观看| 国产精品久久久久久久久免| 亚洲美女视频黄频| 99热这里只有是精品50| 欧美国产精品一级二级三级 | 一级爰片在线观看| 亚洲欧洲日产国产| 91狼人影院| 久久99热这里只频精品6学生| 免费看光身美女| 一二三四中文在线观看免费高清| 三级男女做爰猛烈吃奶摸视频| 欧美xxxx黑人xx丫x性爽| 国产有黄有色有爽视频| 国产淫片久久久久久久久| 男女啪啪激烈高潮av片| 成年免费大片在线观看| 国内揄拍国产精品人妻在线| av在线老鸭窝| 久久久久久久午夜电影| 久久精品国产亚洲av涩爱| 69av精品久久久久久| 成人亚洲精品一区在线观看 | 午夜激情福利司机影院| 国产 一区精品| 天堂俺去俺来也www色官网| 成年女人看的毛片在线观看| 舔av片在线| 国产精品国产三级国产专区5o| 成人高潮视频无遮挡免费网站| 国产精品一及| 免费看日本二区| 午夜日本视频在线| 亚洲综合精品二区| 特级一级黄色大片| 亚洲国产欧美人成| 国产日韩欧美亚洲二区| 欧美日韩视频精品一区| 亚洲欧美一区二区三区国产| 美女脱内裤让男人舔精品视频| 视频中文字幕在线观看| 日韩一区二区视频免费看| 99久久精品热视频| 一级av片app| 久久久午夜欧美精品| 午夜福利视频1000在线观看| 大香蕉久久网| 观看美女的网站| 人妻一区二区av| 国产高清有码在线观看视频| 插逼视频在线观看| 国产人妻一区二区三区在| 欧美亚洲 丝袜 人妻 在线| 香蕉精品网在线| 亚洲国产精品专区欧美| 久久久久久久亚洲中文字幕| 精品久久国产蜜桃| 欧美变态另类bdsm刘玥| 亚洲av福利一区| 国产乱人偷精品视频| 亚洲精品视频女| 亚洲精品乱码久久久久久按摩| 王馨瑶露胸无遮挡在线观看| 欧美+日韩+精品| 成人漫画全彩无遮挡| 久久这里有精品视频免费| 成人国产av品久久久| 成年女人在线观看亚洲视频 | 又爽又黄a免费视频| 26uuu在线亚洲综合色| 国产69精品久久久久777片| 亚洲精品亚洲一区二区| 中文乱码字字幕精品一区二区三区| 夫妻性生交免费视频一级片| 国精品久久久久久国模美| 黄色配什么色好看| 日本一二三区视频观看| 久久久午夜欧美精品| 成人亚洲精品av一区二区| av免费在线看不卡| av专区在线播放| 街头女战士在线观看网站| 亚洲国产成人一精品久久久| 国产精品人妻久久久影院| 国产又色又爽无遮挡免| 国产免费视频播放在线视频| 1000部很黄的大片| a级毛色黄片| 草草在线视频免费看| 亚洲三级黄色毛片| 久久人人爽av亚洲精品天堂 | 男女下面进入的视频免费午夜| 亚洲天堂国产精品一区在线| 91久久精品国产一区二区三区| 免费大片黄手机在线观看| 亚洲欧美清纯卡通| 汤姆久久久久久久影院中文字幕| 国产久久久一区二区三区| 久久久国产一区二区| av在线亚洲专区| 国模一区二区三区四区视频| 亚洲精品亚洲一区二区| 九九爱精品视频在线观看| 大片电影免费在线观看免费| 日韩国内少妇激情av| 能在线免费看毛片的网站| 国产永久视频网站| 国产精品蜜桃在线观看| 中文精品一卡2卡3卡4更新| 亚洲欧美成人精品一区二区| av卡一久久| 国产爱豆传媒在线观看| 狠狠精品人妻久久久久久综合| 自拍欧美九色日韩亚洲蝌蚪91 | 久久这里有精品视频免费| 最近中文字幕高清免费大全6| 夜夜看夜夜爽夜夜摸| 三级国产精品片| 99久久精品热视频| 小蜜桃在线观看免费完整版高清| 天天躁日日操中文字幕| 中文天堂在线官网| 免费看av在线观看网站| 国产一区亚洲一区在线观看| 欧美激情国产日韩精品一区| 2021少妇久久久久久久久久久| 亚洲精品日韩在线中文字幕| 九九爱精品视频在线观看| 一级a做视频免费观看| 精品少妇黑人巨大在线播放| 一本一本综合久久| 亚洲精品成人久久久久久| 国产高清三级在线| 日韩成人av中文字幕在线观看| av.在线天堂| 美女高潮的动态| 亚洲一级一片aⅴ在线观看| 国产伦理片在线播放av一区| 婷婷色综合www| av播播在线观看一区| 老师上课跳d突然被开到最大视频| 激情 狠狠 欧美| 国产高清不卡午夜福利| 国产精品国产三级国产专区5o| 久久久久久伊人网av| 国产精品不卡视频一区二区| 国模一区二区三区四区视频| 蜜臀久久99精品久久宅男| videos熟女内射| 18禁在线播放成人免费| 在线天堂最新版资源| 成年免费大片在线观看| 日韩三级伦理在线观看| 日韩,欧美,国产一区二区三区| 亚洲精品日本国产第一区| 欧美日韩精品成人综合77777| 国产成人精品一,二区| 国产黄片视频在线免费观看| 永久网站在线| 国产v大片淫在线免费观看| 一级片'在线观看视频| 九草在线视频观看| 国产黄色视频一区二区在线观看| 日本色播在线视频| 极品少妇高潮喷水抽搐| 777米奇影视久久| 午夜视频国产福利| 在线观看美女被高潮喷水网站| 亚洲国产精品999| 亚洲av一区综合| 成人免费观看视频高清| 18禁裸乳无遮挡免费网站照片| 九九在线视频观看精品| 一区二区三区乱码不卡18| 老司机影院毛片| 美女cb高潮喷水在线观看| 日本黄大片高清| 下体分泌物呈黄色| 亚洲av国产av综合av卡| 99热6这里只有精品| 少妇人妻一区二区三区视频| 女人被狂操c到高潮| 久久国产乱子免费精品| 内地一区二区视频在线| 亚洲av.av天堂| 精品久久久噜噜| 综合色av麻豆| 免费大片黄手机在线观看| 国模一区二区三区四区视频| 大话2 男鬼变身卡| 午夜激情福利司机影院| 免费观看a级毛片全部| 日日啪夜夜爽| 少妇的逼好多水| 亚洲性久久影院| 天天躁日日操中文字幕| 国产免费福利视频在线观看| 亚洲欧美日韩另类电影网站 | 欧美三级亚洲精品| 狂野欧美激情性bbbbbb| 国产精品无大码| av黄色大香蕉| 男女边吃奶边做爰视频| 高清视频免费观看一区二区| av卡一久久| 一级毛片黄色毛片免费观看视频| 中文乱码字字幕精品一区二区三区| 自拍偷自拍亚洲精品老妇| 免费看日本二区| 日韩伦理黄色片| 直男gayav资源| 一级片'在线观看视频| 特级一级黄色大片| 午夜老司机福利剧场| 三级国产精品欧美在线观看| 秋霞在线观看毛片| 国产一区二区在线观看日韩| 一级爰片在线观看| 在线观看国产h片| .国产精品久久| 午夜爱爱视频在线播放| 97人妻精品一区二区三区麻豆| 极品少妇高潮喷水抽搐| 男女边吃奶边做爰视频| 日韩,欧美,国产一区二区三区| 亚洲欧美清纯卡通| 亚洲真实伦在线观看|