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

    聲-固耦合系統(tǒng)高頻局部能量響應(yīng)的FEM-SmEdA分析

    2019-10-21 06:17王攀陳強(qiáng)李彥斌
    振動(dòng)工程學(xué)報(bào) 2019年4期
    關(guān)鍵詞:有限元法

    王攀 陳強(qiáng) 李彥斌

    摘要: 針對(duì)高頻激勵(lì)下聲-固耦合系統(tǒng)的局部能量響應(yīng)預(yù)示問(wèn)題,采用有限元(FEM)和統(tǒng)計(jì)模態(tài)能量分布分析(SmEdA)相結(jié)合的方法預(yù)示復(fù)雜聲-固耦合系統(tǒng)局部能量響應(yīng)。首先,以平板/聲腔耦合系統(tǒng)為研究對(duì)象,驗(yàn)證FEM-SmEdA方法的準(zhǔn)確性;隨后,為了驗(yàn)證該方法的適用性,分別將該方法應(yīng)用于湍流邊界層(TBL)激勵(lì)下平板/聲腔耦合系統(tǒng)局部響應(yīng)預(yù)示及加筋艙段/聲腔耦合系統(tǒng)局部響應(yīng)預(yù)示。利用雙模態(tài)方程法(DMF)對(duì)基于FEM-SmEdA法的TBL激勵(lì)作用下結(jié)構(gòu)局部響應(yīng)預(yù)示結(jié)果進(jìn)行驗(yàn)證,基于FEM-SmEdA方法揭示TBL激勵(lì)作用下艙段結(jié)構(gòu)和聲腔的能量分布隨頻率的變化規(guī)律,并研究加筋對(duì)聲-固耦合系統(tǒng)能量分布的影響。結(jié)果表明:FEM-SmEdA方法能夠準(zhǔn)確地預(yù)示高頻環(huán)境下聲-固耦合系統(tǒng)的局部能量響應(yīng);該方法適用于TBL激勵(lì)下結(jié)構(gòu)局部響應(yīng)預(yù)示;隨著分析頻率的升高,艙段和聲腔的能量分布更加均勻;加筋后結(jié)構(gòu)的能量分布更加均勻,加筋對(duì)聲腔能量響應(yīng)的抑制更加明顯。

    關(guān)鍵詞: 聲-固耦合系統(tǒng); 統(tǒng)計(jì)模態(tài)能量分布分析; 模態(tài)能量; 能量分布; 有限元法

    中圖分類(lèi)號(hào): O421+.6; V414.4 文獻(xiàn)標(biāo)志碼: A 文章編號(hào): 1004-4523(2019)04-0590-12

    DOI:10.16385/j.cnki.issn.1004-4523.2019.04.005

    引 言

    聲-固耦合問(wèn)題廣泛存在于航空航天結(jié)構(gòu)中,尤其是對(duì)于高頻激勵(lì)下的薄壁結(jié)構(gòu),結(jié)構(gòu)和聲場(chǎng)之間極易產(chǎn)生聲-固耦合效應(yīng),引起結(jié)構(gòu)的振動(dòng)并改變聲場(chǎng)的分布,進(jìn)而影響結(jié)構(gòu)的可靠性和儀器設(shè)備的功能。因此,結(jié)構(gòu)的振動(dòng)水平及聲場(chǎng)的噪聲水平是結(jié)構(gòu)設(shè)計(jì)階段必須考慮的重要指標(biāo),研究薄壁結(jié)構(gòu)在高頻激勵(lì)下的聲-固耦合問(wèn)題具有重要的理論意義和應(yīng)用價(jià)值。

    目前,獲取聲-固耦合響應(yīng)的方法主要有試驗(yàn)方法、理論方法和數(shù)值方法。試驗(yàn)方法的結(jié)果真實(shí)可信,其不足之處在于耗費(fèi)巨大、試驗(yàn)周期較長(zhǎng)且僅能實(shí)現(xiàn)有限實(shí)驗(yàn)條件和工況。理論分析方法較難適用于復(fù)雜工程結(jié)構(gòu)。近年來(lái),隨著計(jì)算機(jī)性能的提高,數(shù)值方法成為一種有效的分析手段[1]。對(duì)于低頻聲-固耦合問(wèn)題,通常采用有限元法(Finite Element Method, FEM)[2]和邊界元法等離散化方法預(yù)示結(jié)構(gòu)響應(yīng)。在結(jié)構(gòu)單元尺寸滿(mǎn)足同一波長(zhǎng)內(nèi)包含6個(gè)單元的條件后[3],離散化方法能較為精確地預(yù)示聲-固耦合系統(tǒng)在低頻段的局部響應(yīng)。但當(dāng)分析頻率升高時(shí),若要保證FEM有足夠的精度,須對(duì)單元進(jìn)行細(xì)化,這將導(dǎo)致計(jì)算規(guī)模呈指數(shù)級(jí)增長(zhǎng)。對(duì)于高頻聲-固耦合問(wèn)題,通常采用統(tǒng)計(jì)能量分析(Statistical Energy Analysis, SEA)[4-5]預(yù)示結(jié)構(gòu)的高頻能量響應(yīng)。SEA將結(jié)構(gòu)劃分為若干子系統(tǒng),通過(guò)引入模態(tài)能量均一化和雨流(Rain On the Roof, ROR)載荷等基本假設(shè),構(gòu)建子系統(tǒng)間的功率流平衡方程,進(jìn)而預(yù)示子系統(tǒng)在空間和分析頻帶內(nèi)的平均能量響應(yīng)[6]。但SEA的各項(xiàng)假設(shè)在工程應(yīng)用中較難完全滿(mǎn)足,而且僅能得到空間平均能量響應(yīng),難以得到結(jié)構(gòu)局部能量響應(yīng)[7]。

    為了克服SEA中的模態(tài)能量分析均分假設(shè),Maxit等基于雙模態(tài)方程(Dual Modal Formulation, DMF)[8-9],推導(dǎo)了子系統(tǒng)模態(tài)能量平衡方程,提出統(tǒng)計(jì)模態(tài)能量分布分析(Statistical Modal Energy Distribution Analysis, SmEdA)[10-11]。SmEdA理論以分析頻帶內(nèi)子系統(tǒng)每階模態(tài)的能量為基本未知量,考慮了同一子系統(tǒng)中不同模態(tài)能量差異較大的情形。基于SmEdA理論預(yù)示子系統(tǒng)的模態(tài)能量,僅需非耦合子系統(tǒng)在耦合邊界的模態(tài)信息即可獲取子系統(tǒng)的模態(tài)能量。對(duì)于聲-固耦合子系統(tǒng)而言,提高了計(jì)算效率,尤其是聲腔子系統(tǒng)在高頻段模態(tài)數(shù)密集的情況下。此外,SmEdA理論以每階模態(tài)為基本輸入?yún)?shù),可考慮同一子系統(tǒng)內(nèi)模態(tài)受到不均勻激勵(lì)的情況,從而可考慮局部激勵(lì),克服SEA理論中雨流載荷的假設(shè)。Aragonès等[12]基于SmEdA理論對(duì)典型聲腔/平板/聲腔耦合系統(tǒng)進(jìn)行了聲功率傳輸路徑分析,識(shí)別出與主要傳遞路徑相關(guān)的平板模態(tài),并在此基礎(chǔ)上對(duì)平板加筋的降噪效果進(jìn)行研究。Totaro等[13]基于SmEdA理論預(yù)示局部激勵(lì)作用下平板/聲腔耦合系統(tǒng)的動(dòng)能和勢(shì)能密度的空間分布。Qiao等[14]利用耦合系統(tǒng)中耦合損耗因子與輻射效率的關(guān)系,基于SmEdA理論研究了圓柱殼和圓錐殼的輻射效率。各子系統(tǒng)的模態(tài)信息為SmEdA理論的重要輸入條件之一,用于計(jì)算各模態(tài)間的耦合損耗因子、各模態(tài)的阻尼耗散功率及各模態(tài)上的載荷輸入功率。對(duì)于簡(jiǎn)單聲-固耦合系統(tǒng)的模態(tài)信息可由解析法給出,而對(duì)于復(fù)雜聲-固耦合系統(tǒng),可借助FEM獲取各子系統(tǒng)的模態(tài)信息[15],進(jìn)而結(jié)合SmEdA理論,預(yù)示各子系統(tǒng)的高頻局部能量響應(yīng)。目前,統(tǒng)計(jì)模態(tài)能量分布分析法的應(yīng)用尚限于ROR載荷和單點(diǎn)力載荷,ROR載荷即空間上處處不相關(guān)的面壓載荷。工程上常把噪聲載荷處理成完全不相關(guān)的面壓載荷[16],而航天器面臨的實(shí)際噪聲環(huán)境存在空間相關(guān)性。因此,有必要開(kāi)展考慮載荷空間相關(guān)性的統(tǒng)計(jì)模態(tài)能量分布分析法研究。

    本文采用FEM和SmEdA結(jié)合的方法,對(duì)于復(fù)雜聲-固耦合系統(tǒng)局部能量響應(yīng)進(jìn)行預(yù)示。首先,以平板/聲腔耦合系統(tǒng)為研究對(duì)象驗(yàn)證FEM-SmEdA方法的準(zhǔn)確性;隨后,分別將該方法應(yīng)用于湍流邊界層(Turbulent Boundary Layer,TBL)激勵(lì)下平板/聲腔耦合系統(tǒng)局部能量響應(yīng)預(yù)示及艙段/聲腔耦合系統(tǒng)的高頻局部響應(yīng)預(yù)示,利用DMF法對(duì)基于FEM-SmEdA法TBL激勵(lì)下結(jié)構(gòu)局部響應(yīng)預(yù)示結(jié)果進(jìn)行驗(yàn)證,并基于FEM-SmEdA法揭示結(jié)構(gòu)和聲腔能量隨著頻率的變化規(guī)律,研究加筋對(duì)聲-固耦合系統(tǒng)能量分布的影響。

    由式(24)和(25)可知,基于SmEdA法預(yù)示結(jié)構(gòu)的能量響應(yīng),不需要整體結(jié)構(gòu)的剛度矩陣,僅需節(jié)點(diǎn)質(zhì)量及節(jié)點(diǎn)剛度。因此,基于SmEdA可以預(yù)示結(jié)構(gòu)的局部能量響應(yīng),克服了傳統(tǒng)的SEA方法僅能預(yù)示結(jié)構(gòu)平均能量響應(yīng)的不足。此外,由于結(jié)構(gòu)的特定模態(tài)信息在分析頻帶內(nèi)是不變的,根據(jù)方程(24)和(25)可知,能量分布的準(zhǔn)確性取決于頻帶內(nèi)基于SmEdA獲得的模態(tài)能量。

    2 方法驗(yàn)證

    以平板/聲腔耦合系統(tǒng)為對(duì)象,基于FEM-SmEdA開(kāi)展結(jié)構(gòu)的局部能量預(yù)示研究,通過(guò)與有限元方法的對(duì)比驗(yàn)證本文方法的準(zhǔn)確性。平板和聲腔耦合示意圖如圖 2所示,平板和聲腔的幾何尺寸及材料屬性分別如表1和2所示。參考文獻(xiàn)[12,15,18],平板和聲腔的結(jié)構(gòu)阻尼均取0.01。平板四周簡(jiǎn)支,聲腔除與平板耦合面,其余各面為固定邊界。

    2.1 模態(tài)能量分析

    為保證有限元計(jì)算結(jié)果的準(zhǔn)確性,以中心頻率為630 Hz的1/3倍頻程頻帶[561, 707] Hz作為分析頻段Δω開(kāi)展研究。基于商用有限元軟件獲得平板及聲腔的模態(tài)信息,頻帶Δω內(nèi)平板和聲腔分別有11階和20階模態(tài)振型。將平板和聲腔的模態(tài)計(jì)算結(jié)果代入式(2),可以求得如圖 3所示的平板與聲腔之間模態(tài)耦合損耗因子。

    在平板坐標(biāo)為(0.7, 0.2)的一點(diǎn)施加垂直于板面方向的單位力激勵(lì)F,由式(6)得到平板模態(tài)上載荷的輸入功率,由式(7)得到平板和聲腔在分析頻帶Δω內(nèi)的模態(tài)能量,如圖 4所示。圖中dB表示兩個(gè)量的比值大小,對(duì)于輸入功率dB_P=10lg(P/ref_P),P表示輸入功率,ref_P表示參考值,且ref_P=1 W。對(duì)于振動(dòng)能量dB_E=10lg(E/ref_E),E表示振動(dòng)能量,ref_E表示參考值,且ref_E=1 J。由圖 4可知,在集中力載荷作用下,即子系統(tǒng)受到局部激勵(lì)的情況下,各階模態(tài)上的載荷輸入功率不同;平板和聲腔的各階模態(tài)能量均具有明顯的差異,平板的各階模態(tài)能量差值最大為20 dB,聲腔差值最大為25 dB,不滿(mǎn)足傳統(tǒng)SEA中的模態(tài)能量均一化假設(shè);平板的各階模態(tài)能量變化趨勢(shì)與模態(tài)上的載荷輸入功率一致。

    圖5為分別基于FEM-SmEdA和FEM預(yù)示得到分析頻帶Δω內(nèi)平板的動(dòng)能分布,圖 6為分別基于FEM-SmEdA和FEM預(yù)示得到分析頻帶Δω內(nèi)平板的勢(shì)能分布。

    從圖5及6對(duì)比可知,兩者存在一定差異的。分析主要存在以下兩點(diǎn)原因:(1)SmEdA理論基于僅考慮了共振傳輸,由于本文研究對(duì)象結(jié)構(gòu)阻尼為0.01,對(duì)于小阻尼系統(tǒng),子系統(tǒng)間的非共振傳輸可忽略,因此本文研究中只考慮分析頻帶內(nèi)的共振模態(tài)。(2)該方法強(qiáng)調(diào)“系統(tǒng)高頻局部能量響應(yīng)預(yù)示”,適用于子系統(tǒng)有較多模態(tài)數(shù)的分析頻段,對(duì)于足夠大的頻帶,滿(mǎn)足模態(tài)數(shù)足夠多,可以忽略同一子系統(tǒng)中不同模態(tài)間的能量,即忽略Tinp(Q,ω)和Vinp(Q,ω)的影響。隨著分析頻段中心頻率的升高,分析頻段的模態(tài)數(shù)增多,忽略同一子系統(tǒng)中不同模態(tài)間的能量所造成的誤差減小,基于FEM-SmEdA預(yù)示的局部響應(yīng)結(jié)果更準(zhǔn)確。

    3 湍流邊界層激勵(lì)下結(jié)構(gòu)局部能量響應(yīng)預(yù)示 ?正如已經(jīng)在文獻(xiàn)[19]指出的,SmEdA理論以每階模態(tài)為基本輸入?yún)?shù),可考慮同一子系統(tǒng)內(nèi)模態(tài)受到不均勻激勵(lì)的情況,克服SEA理論中雨流載荷的假設(shè)。為了證明FEM-SmEdA也適用于預(yù)測(cè)復(fù)雜環(huán)境激勵(lì)下結(jié)構(gòu)的局部響應(yīng)預(yù)示,在平板/聲腔耦合系統(tǒng)中考慮TBL激勵(lì)。如1.2節(jié)已經(jīng)指出的,平板和聲腔的特定模態(tài)信息是不變的,因此能量分布的精度由基于FEM-SmEdA法在頻帶中獲得的模態(tài)能量決定。由于TBL激勵(lì)下系統(tǒng)的局部響應(yīng)難以用有限元方法預(yù)測(cè),因此在第3.2節(jié)中用DMF方法驗(yàn)證了TBL激勵(lì)下子系統(tǒng)模態(tài)能量預(yù)示的準(zhǔn)確性。

    3.1 湍流邊界層激勵(lì)

    3.2 模態(tài)能量分析

    在第2節(jié)所示的平板/聲腔耦合系統(tǒng)中施加TBL激勵(lì),給出了U=200 m/s的來(lái)流速度,分別基于FEM-SmEdA和DMF預(yù)示中心頻率為630 Hz的1/3倍頻程頻帶內(nèi)板和聲腔的模態(tài)能量,如圖 7和8所示。平板和聲腔的模態(tài)能量具有明顯的離散性,平板模態(tài)能量隨模態(tài)變化的能量差值可達(dá)20 dB,聲腔模態(tài)能量隨模態(tài)變化的能量差值可達(dá)30 dB,表明在考慮TBL激勵(lì)的情況下,能量均分的假設(shè)是不合理的。結(jié)果表明,基于FEM-SmEdA和DMF方法預(yù)示的子系統(tǒng)各階模態(tài)的能量值基本一致。如圖 8所示,聲腔的模態(tài)能量略有差異,且兩個(gè)結(jié)果之間的差值小于1 dB。因此,基于FEM-SmEdA得到的TBL激勵(lì)下子系統(tǒng)的模態(tài)能量是準(zhǔn)確合理的。

    3.3 結(jié)構(gòu)局部能量響應(yīng)預(yù)示

    基于FEM-SmEdA預(yù)示TBL激勵(lì)下平板的動(dòng)能和勢(shì)能分布,如圖9和10所示。通過(guò)與圖5和6對(duì)比可知,在中心頻率均為630 Hz的1/3倍頻程頻帶內(nèi),TBL激勵(lì)下平板的局部響應(yīng)與單點(diǎn)力作用不同。

    4 艙段/聲腔耦合系統(tǒng)局部能量響應(yīng)預(yù)示

    4.1 艙段/聲腔耦合系統(tǒng)能量分布預(yù)示 ?以廣泛用于航空航天器的艙段/聲腔耦合系統(tǒng)為研究對(duì)象,開(kāi)展湍流邊界層激勵(lì)下耦合結(jié)構(gòu)的局部能量響應(yīng)預(yù)示。艙段和聲腔的有限元模型如圖11所示,其幾何尺寸及材料屬性分別如表3和4所示。

    基于FEM獲得艙段及聲腔的模態(tài)結(jié)果,由式(2)可得分析頻段Δω1(561.2-2828.4 Hz)內(nèi)艙段與聲腔之間的模態(tài)耦合損耗因子,如圖12所示。由圖12可知,較大的模態(tài)耦合損耗因子多存在于共振模態(tài)間,非共振模態(tài)間的耦合損耗因子相對(duì)較小。由文獻(xiàn)[18]可知:對(duì)于小阻尼系統(tǒng),子系統(tǒng)間的非共振傳輸可忽略,因此式(7)中可只考慮分析頻帶內(nèi)的共振模態(tài)。本文研究對(duì)象屬于小阻尼系統(tǒng),為提高計(jì)算效率,分析中僅考慮共振模態(tài)間的功率傳輸。

    基于FEM-SmEdA計(jì)算艙段和聲腔在中心頻率分別為630,1000,2500,4000 Hz的1/3倍頻程頻帶內(nèi)的能量空間分布,計(jì)算結(jié)果如表 5所示。由表 5可知:(1)在不同頻段內(nèi),結(jié)構(gòu)及聲腔的動(dòng)能和勢(shì)能分布基本一致,艙段動(dòng)能的最大值高于勢(shì)能的最大值,而聲腔勢(shì)能的最大值高于動(dòng)能的最大值。(2)結(jié)構(gòu)的能量響應(yīng)遠(yuǎn)大于聲腔。(3)隨著分析頻率的升高,艙段和聲腔的能量分布更加均勻。當(dāng)分析頻帶中心頻率分別為630,1000,2500 Hz時(shí),艙段和聲腔的能量分布不均勻程度較大;當(dāng)分析頻帶中心頻率為4000 Hz時(shí),艙段和聲腔的能量分布較為均勻。隨著分析頻率的升高,艙段和聲腔的能量分布更接近SEA中子系統(tǒng)能量平均假設(shè),此時(shí)使用SEA方法計(jì)算艙段的平均能量即可較為準(zhǔn)確地表征子系統(tǒng)的能量大小。

    4.2 加筋艙段/聲腔耦合系統(tǒng)能量分布預(yù)示

    為了驗(yàn)證FEM-SmEdA的適用性,以復(fù)雜加筋艙段/聲腔耦合系統(tǒng)為研究對(duì)象,基于FEM-SmEdA預(yù)示聲-固耦合系統(tǒng)的能量分布。加筋艙段的幾何模型如圖13所示。其幾何尺寸及材料屬性如表3所示。筋條包括8根縱向筋,3根環(huán)向筋。筋條材料屬性與艙段一致,筋條高度為0.02 m、厚度為4 mm。

    基于FEM獲得加筋艙段及聲腔的模態(tài)結(jié)果,在中心頻率為1000 Hz的1/3倍頻程頻段Δω2(890.9-1122.5 Hz)內(nèi),加筋艙段和聲腔分別有14階和12階模態(tài)振型。將加筋艙段和聲腔的模態(tài)結(jié)果代入式(2),可得如圖 14所示的Δω2頻段內(nèi)加筋艙段與聲腔之間模態(tài)耦合損耗因子。由圖 14可知,Δω2頻段內(nèi)的共振模態(tài)間的模態(tài)耦合損耗因子只有少數(shù)模態(tài)耦合損耗因子較大,絕大多數(shù)模態(tài)耦合損耗因子較小,不滿(mǎn)足均一化假設(shè)。

    在加筋艙段表面施加湍流邊界層激勵(lì),基于FEM-SmEdA分別計(jì)算加筋艙段和聲腔在Δω2內(nèi)的能量分布,如表 6所示。

    通過(guò)對(duì)比表 5和6可知,加筋后艙段結(jié)構(gòu)和聲腔的能量分布更加均勻;加筋后艙段和聲腔能量分布最大值變小。由4.1節(jié)圖 12可得分析頻段Δω2(890.9-1122.5 Hz)內(nèi)艙段與聲腔之間的模態(tài)耦合損耗因子,如圖 15所示。對(duì)比圖 14和15可知,加筋后艙段與聲腔的耦合強(qiáng)度減弱。這主要是由于筋條改變了艙段的固有頻率,并同時(shí)抑制了部分模態(tài)振型,使得艙段結(jié)構(gòu)向聲腔內(nèi)傳遞的能量減少。

    為更加深入地說(shuō)明加筋對(duì)艙段與聲腔間的耦合強(qiáng)度影響,對(duì)加筋前后的艙段在分析頻帶內(nèi)施加相同的輸入功率,對(duì)比加筋前后艙段和聲腔的總能量,如表 7所示。由表 7可知,加筋后艙段的能量減少了21.4%,聲腔能量減少了82.0%,加筋對(duì)聲腔響應(yīng)抑制更加明顯。

    5 結(jié) 論

    本文結(jié)合FEM和SmEdA理論,形成了適用于聲-固耦合系統(tǒng)的高頻局部能量響應(yīng)預(yù)示方法。首先,以平板/聲腔耦合系統(tǒng)為研究對(duì)象,通過(guò)與FEM的對(duì)比驗(yàn)證了FEM-SmEdA方法的準(zhǔn)確性;隨后,分別將該方法應(yīng)用于TBL激勵(lì)下結(jié)構(gòu)局部響應(yīng)預(yù)示及艙段/聲腔耦合系統(tǒng)的高頻局部能量響應(yīng)預(yù)示,利用DMF法對(duì)基于FEM-SmEdA法的TBL激勵(lì)下結(jié)構(gòu)局部響應(yīng)預(yù)示結(jié)果進(jìn)行驗(yàn)證,并將該方法應(yīng)用于艙段/聲腔耦合系統(tǒng),揭示艙段結(jié)構(gòu)和聲腔的能量分布隨頻率的變化規(guī)律,并研究加筋對(duì)聲-固耦合系統(tǒng)能量分布的影響。結(jié)果表明:

    1)FEM-SmEdA預(yù)示得到的能量響應(yīng)大小和空間分布與FEM預(yù)示結(jié)果基本一致,這說(shuō)明FEM-SmEdA方法能夠準(zhǔn)確地預(yù)示聲-固耦合系統(tǒng)的局部能量響應(yīng);

    2)基于FEM-SmEdA和DMF方法預(yù)示的子系統(tǒng)各階模態(tài)能量結(jié)果的基本一致。因此,基于FEM-SmEdA得到的TBL激勵(lì)下子系統(tǒng)的模態(tài)能量是準(zhǔn)確合理的,該方法適于TBL激勵(lì)下結(jié)構(gòu)的局部響應(yīng)預(yù)示。

    3)對(duì)于TBL激勵(lì)下的艙段/聲腔耦合系統(tǒng),艙段及聲腔的動(dòng)能和勢(shì)能的分布基本一致,艙段的能量響應(yīng)遠(yuǎn)大于聲腔;隨著分析頻率的升高,艙段和聲腔的能量分布更加均勻;

    4)對(duì)于TBL激勵(lì)下的艙段/聲腔耦合系統(tǒng),加筋后艙段結(jié)構(gòu)和聲腔的能量分布更加均勻;加筋后艙段和聲腔能量分布最大值變小;與加筋后艙段的總能量減少值相比,聲腔總能量減少較大,故加筋對(duì)聲腔響應(yīng)抑制更加明顯。

    參考文獻(xiàn):

    [1] 李彥斌, 張 鵬, 吳邵慶,等. 復(fù)合材料加筋板計(jì)及熱效應(yīng)的聲-固耦合分析[J]. 振動(dòng)工程學(xué)報(bào), 2015,28(4):531-540.

    Li Yanbin, Zhang Peng, Wu Shaoqing, et al. Structural- acoustic coupling analysis of a composite stiffened panel in a thermal environment[J]. Journal of Vibration Engineering, 2015,28(4):531-540.

    [2] Yao L Y, Yu D J, Cui X Y, et al. Numerical treatment of acoustic problems with the smoothed finite element method[J]. Applied Acoustics, 2010, 71(8):743-753.

    [3] 李彥斌, 張 鵬, 鄒元杰,等. 隨機(jī)基礎(chǔ)激勵(lì)下承力筒-蒙皮結(jié)構(gòu)的聲-固耦合分析[J]. 宇航學(xué)報(bào), 2015, 36(2):236-242.

    Li Yanbin, Zhang Peng, Zou Yuanjie, et al. Structural- acoustic coupling analysis of a bearing cylinder-skin structure under random base excitation[J]. Journal of Astronautics, 2015, 36(2):236-242.

    [4] Lyon R H. Theory and Application of Statistical Energy Analysis[M]. Elsevier, 2014.

    [5] Chen Qiang, Fei Qingguo, Wu Shaoqing, et al. An efficient transient analysis method for time-varying structures based on statistical energy analysis[J]. Mechanics Research Communications, 2018, 91:93-99.

    [6] Chen Qiang, Fei Qingguo, Li Yanbin, et al. Prediction of the transient energy response for complex vibro-acoustic systems[J]. Journal of Mechanical Science and Technology, 2019, 33(2): 495-504.

    [7] Mace B R. Statistical energy analysis: Coupling loss factors, indirect coupling and system modes[J]. Journal of Sound & Vibration, 2005, 279(1):141-170.

    [8] Karnopp D. Coupled vibratory-system analysis, using the dual formulation[J]. Journal of the Acoustical Society of America, 1966, 40 (2):380-384.

    [9] Tabarrok B. Dual formulations for acoustic-structural vibrations[J]. International Journal for Numerical Methods in Engineering, 1978, 13(1):197-201.

    [10] Maxit L, Guyader J L. Estimation of SEA coupling loss factors using a dual formulation and FEM modal information, part I: Theory[J]. Journal of Sound and Vibration, 2001, 238 (5):907-930.

    [11] Maxit L, Guyader J L. Estimation of SEA coupling loss factors using a dual formulation and FEM modal information, part II: Numerical applications[J]. Journal of Sound and Vibration, 2001, 238 (5):931-948.

    [12] ngels Aragonès, Maxit L, Guasch O. A graph theory approach to identify resonant and non-resonant transmission paths in statistical modal energy distribution analysis[J]. Journal of Sound & Vibration, 2015, 350(8):91-110.

    [13] Totaro N, Guyader J L. Extension of the statistical modal energy distribution analysis for estimating energy density in coupled subsystems[J]. Journal of Sound & Vibration, 2012,331(13):3114-3129.

    [14] Qiao Y, Chen H B, Luo J L. Estimation of shell radiation efficiency using a FEM-SmEdA algorithm[J]. Journal of Vibroengineering, 2013,15(3):1130-1146.

    [15] 于士甲, 張 鵬, 李彥斌,等. 基于FEM-MODENA的加筋板聲-固耦合分析[J]. 振動(dòng)工程學(xué)報(bào), 2017,30(2):262-269.

    Yu Shijia, Zhang Peng, Li Yanbin, et al. Structural-acoustic coupling analysis of a stiffened plate based on FEM-MODENA[J]. Journal of Vibration Engineering, 2017,30(2):262-269.

    [16] Ichchou M N, Hiverniau B, Troclet B. Equivalent ‘rain on the roof loads for random spatially correlated excitations in the mid-high frequency range[J]. Journal of Sound and Vibration, 2009,322(4):926-940.

    [17] Karnopp D. Coupled vibratory-system analysis, using the dual formulation[J]. Journal of the Acoustical Society of America, 2005,40(2):380-384.

    [18] Maxit L, Ege K, Totaro N, et al. Non resonant transmission modelling with statistical modal energy distribution analysis[J]. Journal of Sound & Vibration, 2013, 333(2):499-519.

    [19] Maxit L, Guyader J L. Extension of SEA model to subsystems with non-uniform modal energy distribution[J]. Journal of Sound and Vibration, 2003, 265(2): 337-358.

    [20] Corcos G M. The structure of the turbulent pressure field in boundary-layer flows[J]. Journal of Fluid Mechanics, 1964, 18(3):353-378.

    [21] Graham W R. A comparison of models for the wavenumber frequency spectrum of turbulent boundary layer pressures[J]. Journal of Sound and Vibration, 1997, 206(4): 541-565.

    Abstract: By combining the finite element method (FEM) and the statistical modal energy distribution analysis (SmEdA), a FEM-SmEdA approach is presented to predict the local energy response of structural-acoustic coupling system under high frequency excitation. First, the accuracy of the FEM-SmEdA is verified by a plate/cavity coupling system. Then the approach is applied to a plate/cavity coupling system and a cabin/cavity coupling system under Turbulent Boundary Layer (TBL) excitation to verify its applicability respectively. The energy responses predicted by FEM-SmEdA are verified by Dual Modal Formulation (DMF) method. The influence of analysis frequency band and stiffener on the energy distribution of cabin/cavity coupling system under TBL excitation are investigated by the FEM-SmEdA. Results show that the FEM-SmEdA is capable of predicting the local energy response of structural-acoustic coupling system under the high frequency excitation with sufficient accuracy. This approach is well adapted for predicting the local response of systems under TBL excitation. With the increase of analysis frequency band, the energy distributions of cabin and cavity become more uniform. The energy distribution of stiffened cabin/cavity coupling system is more uniform than that of cabin/cavity coupling system. The suppression of stiffener to response of cavity is more obvious than structure.

    Key words: structural-acoustic coupling system; SmEdA; modal energy; energy distribution; FEM

    作者簡(jiǎn)介: 王 攀(1993-),男,碩士研究生。電話(huà): (025)83790168; E-mail: 220160935@seu.edu.cn

    通訊作者: 費(fèi)慶國(guó)(1977-),男,博士,教授,博士生導(dǎo)師。電話(huà): (025)83790168; E-mail: qgfei@seu.edu.cn

    猜你喜歡
    有限元法
    帶式輸送機(jī)卸料小車(chē)車(chē)架結(jié)構(gòu)靜力分析與結(jié)構(gòu)改進(jìn)設(shè)計(jì)研究
    壩基面受力狀態(tài)的有限元計(jì)算誤差影響分析
    轉(zhuǎn)向驅(qū)動(dòng)橋轉(zhuǎn)向節(jié)臂的分析與改進(jìn)
    機(jī)械有限元課程在本科教學(xué)中的建設(shè)與實(shí)踐
    機(jī)械類(lèi)碩士生有限元法課程教學(xué)方法研究
    隧洞圍巖錨桿支護(hù)模擬方法對(duì)比分析
    CFRP補(bǔ)強(qiáng)混凝土板彎矩作用下應(yīng)力問(wèn)題研究
    基于非線(xiàn)性有限元的空氣彈簧垂向剛度分析
    二階傳輸條件撕裂互連法在電大輻射中的應(yīng)用
    有限元法模擬GFRP筋肋深與其拉伸力學(xué)性能關(guān)系研究
    女性生殖器流出的白浆| 欧美激情 高清一区二区三区| 久久国产精品男人的天堂亚洲| 欧美 日韩 精品 国产| 午夜91福利影院| 成年人午夜在线观看视频| 侵犯人妻中文字幕一二三四区| 一级a爱视频在线免费观看| 久热爱精品视频在线9| 999精品在线视频| 纯流量卡能插随身wifi吗| 新久久久久国产一级毛片| 欧美97在线视频| 国产欧美日韩一区二区三 | 日本五十路高清| 精品亚洲成国产av| 亚洲成av片中文字幕在线观看| 国产男女超爽视频在线观看| 亚洲一区二区三区欧美精品| 日本wwww免费看| 欧美另类亚洲清纯唯美| 正在播放国产对白刺激| 老司机午夜福利在线观看视频 | 色婷婷av一区二区三区视频| 国产精品一区二区在线观看99| 欧美老熟妇乱子伦牲交| 亚洲 国产 在线| 一边摸一边做爽爽视频免费| 欧美国产精品va在线观看不卡| 脱女人内裤的视频| 大码成人一级视频| 一本—道久久a久久精品蜜桃钙片| 国产黄色免费在线视频| 午夜视频精品福利| 一二三四在线观看免费中文在| 国产高清视频在线播放一区 | 成在线人永久免费视频| 国产精品一区二区在线观看99| 51午夜福利影视在线观看| 老鸭窝网址在线观看| 69精品国产乱码久久久| 欧美另类亚洲清纯唯美| 亚洲成人手机| 亚洲va日本ⅴa欧美va伊人久久 | 免费日韩欧美在线观看| 高潮久久久久久久久久久不卡| av在线播放精品| 91精品国产国语对白视频| av天堂久久9| 欧美性长视频在线观看| 久久亚洲精品不卡| 国产精品一区二区免费欧美 | 制服诱惑二区| 夫妻午夜视频| av免费在线观看网站| 欧美97在线视频| 99国产精品一区二区蜜桃av | 久久香蕉激情| 久久久久久久久免费视频了| 日韩一区二区三区影片| 十八禁网站免费在线| 两性夫妻黄色片| 欧美久久黑人一区二区| 国产欧美日韩一区二区三区在线| 久久久久精品人妻al黑| 亚洲国产av影院在线观看| 午夜福利视频精品| av在线老鸭窝| 国产精品 国内视频| 一级毛片电影观看| 亚洲精品中文字幕在线视频| 色94色欧美一区二区| 18禁黄网站禁片午夜丰满| 一级毛片精品| 纯流量卡能插随身wifi吗| 国产精品.久久久| 国产国语露脸激情在线看| 啦啦啦视频在线资源免费观看| 成人av一区二区三区在线看 | 国产精品熟女久久久久浪| 午夜福利视频在线观看免费| 香蕉丝袜av| 18禁裸乳无遮挡动漫免费视频| 国产成+人综合+亚洲专区| videos熟女内射| 国产有黄有色有爽视频| 亚洲成国产人片在线观看| 欧美 亚洲 国产 日韩一| 午夜影院在线不卡| 亚洲av国产av综合av卡| 国精品久久久久久国模美| 亚洲精品国产av蜜桃| 一个人免费在线观看的高清视频 | 婷婷色av中文字幕| www.999成人在线观看| 免费观看a级毛片全部| 久久国产亚洲av麻豆专区| 黑人巨大精品欧美一区二区蜜桃| videosex国产| 自拍欧美九色日韩亚洲蝌蚪91| www日本在线高清视频| 日韩 欧美 亚洲 中文字幕| 国产成人精品久久二区二区91| 2018国产大陆天天弄谢| 亚洲精品成人av观看孕妇| 宅男免费午夜| 成人黄色视频免费在线看| 一区二区三区四区激情视频| 亚洲国产毛片av蜜桃av| 亚洲激情五月婷婷啪啪| 中文字幕av电影在线播放| 自线自在国产av| 99精国产麻豆久久婷婷| 国产人伦9x9x在线观看| 男人舔女人的私密视频| 高清黄色对白视频在线免费看| 高清在线国产一区| 久久久久久久国产电影| 黑丝袜美女国产一区| 大码成人一级视频| 黄色 视频免费看| 大香蕉久久网| 大码成人一级视频| 一级黄色大片毛片| videosex国产| 久久久久精品国产欧美久久久 | 男女下面插进去视频免费观看| 亚洲欧洲精品一区二区精品久久久| 超碰97精品在线观看| 亚洲精品国产精品久久久不卡| 精品国产一区二区久久| 黄色毛片三级朝国网站| 国产精品一区二区精品视频观看| 欧美日韩av久久| 精品人妻一区二区三区麻豆| 精品一品国产午夜福利视频| 欧美黑人欧美精品刺激| 乱人伦中国视频| 他把我摸到了高潮在线观看 | 男女下面插进去视频免费观看| 操美女的视频在线观看| 精品少妇久久久久久888优播| 又紧又爽又黄一区二区| 国产成人系列免费观看| 五月开心婷婷网| 免费高清在线观看视频在线观看| 国产亚洲午夜精品一区二区久久| 美女福利国产在线| 国产精品香港三级国产av潘金莲| 国产欧美日韩一区二区精品| 一本大道久久a久久精品| 考比视频在线观看| 人人妻人人澡人人看| 国产精品熟女久久久久浪| kizo精华| av片东京热男人的天堂| 精品少妇一区二区三区视频日本电影| 国产免费现黄频在线看| 日本猛色少妇xxxxx猛交久久| 另类精品久久| 天堂中文最新版在线下载| 午夜成年电影在线免费观看| 999久久久国产精品视频| 99国产综合亚洲精品| 2018国产大陆天天弄谢| 婷婷色av中文字幕| 99国产精品一区二区蜜桃av | svipshipincom国产片| 亚洲国产欧美日韩在线播放| 黄频高清免费视频| 久久人人爽av亚洲精品天堂| 国产成人系列免费观看| 99香蕉大伊视频| 午夜福利影视在线免费观看| 一本一本久久a久久精品综合妖精| 视频区图区小说| 不卡av一区二区三区| 亚洲,欧美精品.| 热99re8久久精品国产| 老熟妇乱子伦视频在线观看 | 日韩中文字幕视频在线看片| 亚洲精品国产区一区二| 久久99一区二区三区| 丰满迷人的少妇在线观看| 自线自在国产av| 国产成人精品无人区| 亚洲欧美日韩高清在线视频 | 亚洲av电影在线进入| 爱豆传媒免费全集在线观看| 亚洲国产成人一精品久久久| 最近最新中文字幕大全免费视频| 精品高清国产在线一区| 国产又色又爽无遮挡免| 一进一出抽搐动态| 色老头精品视频在线观看| 午夜福利,免费看| 亚洲人成电影免费在线| 欧美精品啪啪一区二区三区 | 欧美日韩亚洲国产一区二区在线观看 | 欧美日本中文国产一区发布| 丰满人妻熟妇乱又伦精品不卡| 欧美 日韩 精品 国产| 狠狠狠狠99中文字幕| 久久综合国产亚洲精品| 美女午夜性视频免费| 两性夫妻黄色片| 国产亚洲精品一区二区www | 国产精品国产av在线观看| 99精品久久久久人妻精品| 岛国在线观看网站| 久久精品人人爽人人爽视色| 免费av中文字幕在线| 韩国精品一区二区三区| 国产精品偷伦视频观看了| 黑人操中国人逼视频| 考比视频在线观看| 久久久久久亚洲精品国产蜜桃av| 中文字幕制服av| 久久99一区二区三区| 日日爽夜夜爽网站| 亚洲国产精品999| 国产区一区二久久| 成人av一区二区三区在线看 | 欧美亚洲日本最大视频资源| 国产精品 国内视频| 丝袜美足系列| 日本精品一区二区三区蜜桃| 午夜久久久在线观看| 久久久久久人人人人人| 亚洲精品久久久久久婷婷小说| av电影中文网址| 黄片小视频在线播放| e午夜精品久久久久久久| 精品少妇一区二区三区视频日本电影| 最近最新中文字幕大全免费视频| 丰满少妇做爰视频| av有码第一页| 久久久久久久精品精品| 国产91精品成人一区二区三区 | 免费在线观看视频国产中文字幕亚洲 | 久久久久视频综合| 午夜精品国产一区二区电影| 免费高清在线观看视频在线观看| 777久久人妻少妇嫩草av网站| 日韩熟女老妇一区二区性免费视频| 伦理电影免费视频| 精品一品国产午夜福利视频| 成年女人毛片免费观看观看9 | 国产成人免费无遮挡视频| avwww免费| 国产欧美日韩综合在线一区二区| 国产亚洲一区二区精品| 久久久久久久久久久久大奶| 欧美精品啪啪一区二区三区 | 人人妻人人添人人爽欧美一区卜| 午夜福利一区二区在线看| 日韩制服骚丝袜av| 少妇粗大呻吟视频| a在线观看视频网站| 午夜福利影视在线免费观看| 亚洲人成77777在线视频| 久久人人爽av亚洲精品天堂| 久久狼人影院| 首页视频小说图片口味搜索| av天堂在线播放| 免费人妻精品一区二区三区视频| a 毛片基地| 叶爱在线成人免费视频播放| 亚洲国产欧美网| 1024香蕉在线观看| 一本色道久久久久久精品综合| 亚洲熟女毛片儿| 亚洲中文av在线| 亚洲一卡2卡3卡4卡5卡精品中文| 中国国产av一级| 每晚都被弄得嗷嗷叫到高潮| 99re6热这里在线精品视频| 最近最新免费中文字幕在线| 一区二区三区激情视频| 久久国产精品男人的天堂亚洲| 国产男女超爽视频在线观看| 中国美女看黄片| 久久天躁狠狠躁夜夜2o2o| 国产区一区二久久| 久久 成人 亚洲| 97人妻天天添夜夜摸| 欧美日韩亚洲高清精品| 欧美老熟妇乱子伦牲交| 99国产综合亚洲精品| 熟女少妇亚洲综合色aaa.| 老司机影院成人| avwww免费| 精品久久久久久久毛片微露脸 | 日韩制服丝袜自拍偷拍| 亚洲国产精品999| 久久香蕉激情| 国产老妇伦熟女老妇高清| 搡老岳熟女国产| 啪啪无遮挡十八禁网站| 高清av免费在线| 亚洲国产欧美日韩在线播放| 大型av网站在线播放| 国产一区二区 视频在线| 成年av动漫网址| 在线看a的网站| av天堂久久9| 少妇精品久久久久久久| 免费看十八禁软件| 一区二区三区四区激情视频| 欧美成人午夜精品| 十八禁网站网址无遮挡| 日韩制服骚丝袜av| 色播在线永久视频| 精品人妻一区二区三区麻豆| 国产成人影院久久av| 久久99热这里只频精品6学生| 亚洲第一av免费看| 久久精品成人免费网站| 色视频在线一区二区三区| 国产精品麻豆人妻色哟哟久久| 国产精品 国内视频| 国产成人a∨麻豆精品| 日本欧美视频一区| 亚洲一卡2卡3卡4卡5卡精品中文| 高清欧美精品videossex| 日韩视频在线欧美| 国产一区二区三区av在线| 久久 成人 亚洲| tube8黄色片| 免费人妻精品一区二区三区视频| 国产成人免费无遮挡视频| 国产亚洲精品第一综合不卡| 国产精品 国内视频| 日日爽夜夜爽网站| 亚洲精品久久成人aⅴ小说| 99精品久久久久人妻精品| 国产精品熟女久久久久浪| 亚洲男人天堂网一区| 国产97色在线日韩免费| 久久国产精品影院| 一个人免费在线观看的高清视频 | 亚洲午夜精品一区,二区,三区| 老汉色av国产亚洲站长工具| 亚洲第一欧美日韩一区二区三区 | 亚洲自偷自拍图片 自拍| 午夜福利,免费看| 大码成人一级视频| 在线亚洲精品国产二区图片欧美| 精品国产乱码久久久久久男人| 黄色毛片三级朝国网站| 国产成人精品在线电影| 亚洲成人手机| 欧美日韩亚洲综合一区二区三区_| 精品高清国产在线一区| 一个人免费在线观看的高清视频 | 深夜精品福利| 国产成人av激情在线播放| 欧美老熟妇乱子伦牲交| 午夜免费鲁丝| 国产精品.久久久| 丁香六月欧美| 国产男女超爽视频在线观看| 黑人操中国人逼视频| 国产精品自产拍在线观看55亚洲 | 亚洲 欧美一区二区三区| 99re6热这里在线精品视频| 精品一区二区三区av网在线观看 | 中文字幕另类日韩欧美亚洲嫩草| 亚洲专区字幕在线| 大型av网站在线播放| 国产老妇伦熟女老妇高清| 99国产精品一区二区蜜桃av | 中文字幕精品免费在线观看视频| 高清黄色对白视频在线免费看| 国产精品熟女久久久久浪| 男女免费视频国产| 老司机在亚洲福利影院| 日韩大码丰满熟妇| 久久久精品区二区三区| 精品熟女少妇八av免费久了| 99国产极品粉嫩在线观看| 精品国产一区二区久久| 极品少妇高潮喷水抽搐| 各种免费的搞黄视频| 亚洲成人免费av在线播放| 国产成人欧美在线观看 | 久久人人97超碰香蕉20202| 精品熟女少妇八av免费久了| av天堂在线播放| 亚洲免费av在线视频| 欧美av亚洲av综合av国产av| 成人三级做爰电影| 国产欧美日韩一区二区三 | 亚洲精品久久成人aⅴ小说| √禁漫天堂资源中文www| 一区福利在线观看| 91大片在线观看| 国产精品二区激情视频| 久久青草综合色| 成年人免费黄色播放视频| 久久 成人 亚洲| 人人妻人人澡人人爽人人夜夜| 精品少妇内射三级| 精品国产乱码久久久久久男人| 久久精品亚洲熟妇少妇任你| 亚洲激情五月婷婷啪啪| 老汉色∧v一级毛片| 69精品国产乱码久久久| 亚洲少妇的诱惑av| 中文字幕色久视频| 日韩欧美一区二区三区在线观看 | 亚洲五月婷婷丁香| 黄色 视频免费看| 人人妻人人澡人人看| 国产日韩欧美亚洲二区| 亚洲avbb在线观看| 亚洲国产精品999| 久久青草综合色| 久久久国产精品麻豆| 精品卡一卡二卡四卡免费| 免费观看人在逋| a级毛片黄视频| 视频在线观看一区二区三区| 曰老女人黄片| 人人妻人人澡人人爽人人夜夜| 免费av中文字幕在线| 欧美日韩黄片免| 满18在线观看网站| 午夜免费鲁丝| 精品人妻一区二区三区麻豆| 欧美日韩视频精品一区| 亚洲精品一卡2卡三卡4卡5卡 | 国产成+人综合+亚洲专区| 亚洲精品美女久久av网站| 日日爽夜夜爽网站| 十八禁网站免费在线| 精品少妇久久久久久888优播| 女性生殖器流出的白浆| 两性午夜刺激爽爽歪歪视频在线观看 | 国产极品粉嫩免费观看在线| 丝袜喷水一区| 午夜成年电影在线免费观看| 91精品三级在线观看| 精品高清国产在线一区| 男女免费视频国产| 建设人人有责人人尽责人人享有的| 国产av又大| 亚洲专区中文字幕在线| 欧美一级毛片孕妇| 另类亚洲欧美激情| 久久99热这里只频精品6学生| 一个人免费在线观看的高清视频 | 国产精品九九99| 午夜日韩欧美国产| 欧美精品一区二区大全| 日韩 亚洲 欧美在线| 日本av免费视频播放| 999久久久国产精品视频| 一本—道久久a久久精品蜜桃钙片| 亚洲精品久久成人aⅴ小说| 黑人欧美特级aaaaaa片| videosex国产| 亚洲午夜精品一区,二区,三区| 欧美亚洲 丝袜 人妻 在线| 老司机靠b影院| 99久久国产精品久久久| 亚洲国产毛片av蜜桃av| 国产日韩一区二区三区精品不卡| 视频区欧美日本亚洲| av不卡在线播放| 日韩人妻精品一区2区三区| 蜜桃国产av成人99| 日本撒尿小便嘘嘘汇集6| 国产精品香港三级国产av潘金莲| 我的亚洲天堂| 日韩视频一区二区在线观看| 亚洲午夜精品一区,二区,三区| 久久精品亚洲av国产电影网| 精品亚洲乱码少妇综合久久| √禁漫天堂资源中文www| 男女床上黄色一级片免费看| 搡老熟女国产l中国老女人| 男女下面插进去视频免费观看| 一区二区三区激情视频| 精品久久久久久电影网| 亚洲国产看品久久| 国产一区二区 视频在线| 90打野战视频偷拍视频| 黄色视频在线播放观看不卡| 国产av又大| 男女午夜视频在线观看| 亚洲欧洲精品一区二区精品久久久| 欧美激情极品国产一区二区三区| 免费日韩欧美在线观看| 国产av又大| videos熟女内射| av一本久久久久| 国产精品.久久久| 亚洲一区中文字幕在线| 麻豆国产av国片精品| av在线app专区| www.熟女人妻精品国产| 久久久国产精品麻豆| 曰老女人黄片| 日韩有码中文字幕| 精品视频人人做人人爽| 日日夜夜操网爽| 亚洲七黄色美女视频| 精品国产超薄肉色丝袜足j| 国产精品秋霞免费鲁丝片| 精品一区二区三卡| 久久综合国产亚洲精品| 久久人妻福利社区极品人妻图片| 亚洲成人免费电影在线观看| 亚洲av男天堂| 欧美激情 高清一区二区三区| 国产免费一区二区三区四区乱码| 人成视频在线观看免费观看| 999精品在线视频| 国产成人欧美在线观看 | 日韩欧美一区视频在线观看| 免费在线观看日本一区| 满18在线观看网站| 久久久久久人人人人人| 美女扒开内裤让男人捅视频| 国产在线免费精品| 水蜜桃什么品种好| 国产精品秋霞免费鲁丝片| 久久久精品免费免费高清| 在线观看免费高清a一片| 在线看a的网站| 悠悠久久av| 一级片'在线观看视频| 男女床上黄色一级片免费看| 老司机深夜福利视频在线观看 | 精品人妻1区二区| 多毛熟女@视频| 久久久国产精品麻豆| 一级a爱视频在线免费观看| 午夜老司机福利片| 午夜福利乱码中文字幕| 黄色毛片三级朝国网站| a级毛片在线看网站| 久久久久久久大尺度免费视频| 精品久久蜜臀av无| 在线天堂中文资源库| 下体分泌物呈黄色| 精品国产一区二区久久| 国产不卡av网站在线观看| 搡老熟女国产l中国老女人| 桃红色精品国产亚洲av| 天天躁狠狠躁夜夜躁狠狠躁| 黄色怎么调成土黄色| 国产福利在线免费观看视频| 老汉色∧v一级毛片| 一级片'在线观看视频| 国产精品九九99| 午夜激情av网站| 日本wwww免费看| 91麻豆精品激情在线观看国产 | 欧美日韩成人在线一区二区| 免费黄频网站在线观看国产| 久久久久久久精品精品| 人成视频在线观看免费观看| 嫩草影视91久久| svipshipincom国产片| 天天添夜夜摸| 热re99久久精品国产66热6| 首页视频小说图片口味搜索| 十八禁网站网址无遮挡| 欧美日韩av久久| 菩萨蛮人人尽说江南好唐韦庄| 午夜福利一区二区在线看| 在线观看舔阴道视频| 两性夫妻黄色片| 亚洲国产看品久久| 精品少妇内射三级| 黄片大片在线免费观看| 国产又爽黄色视频| 美女福利国产在线| 超色免费av| 欧美av亚洲av综合av国产av| av在线播放精品| 中亚洲国语对白在线视频| 成人国产一区最新在线观看| 中国国产av一级| 亚洲国产精品一区二区三区在线| 18禁观看日本| 青春草亚洲视频在线观看| 亚洲欧美精品自产自拍| 婷婷成人精品国产| 多毛熟女@视频| 欧美日韩亚洲综合一区二区三区_| 亚洲精品中文字幕在线视频| 国产成人啪精品午夜网站| 最近中文字幕2019免费版| 老汉色∧v一级毛片| 夫妻午夜视频| 欧美激情高清一区二区三区| 少妇粗大呻吟视频| 正在播放国产对白刺激| 欧美xxⅹ黑人| 婷婷成人精品国产| 黑人猛操日本美女一级片| 国产男女超爽视频在线观看| 不卡一级毛片| 国产精品 欧美亚洲| 亚洲av电影在线进入| 亚洲精品国产区一区二| √禁漫天堂资源中文www| 国产精品av久久久久免费| 女人被躁到高潮嗷嗷叫费观| 国产高清videossex| 亚洲七黄色美女视频| 老熟女久久久|