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

    梁板組合結(jié)構(gòu)的中頻振動(dòng)試驗(yàn)與計(jì)算分析

    2016-08-04 06:15:43張永杰韋冰峰秦朝紅
    振動(dòng)與沖擊 2016年13期
    關(guān)鍵詞:梁板聲壓級(jí)聲場(chǎng)

    張永杰, 陳 璐, 韋冰峰, 王 喆, 秦朝紅

    (北京強(qiáng)度環(huán)境研究所,北京 100076)

    梁板組合結(jié)構(gòu)的中頻振動(dòng)試驗(yàn)與計(jì)算分析

    張永杰, 陳璐, 韋冰峰, 王喆, 秦朝紅

    (北京強(qiáng)度環(huán)境研究所,北京100076)

    針對(duì)梁板組合結(jié)構(gòu)的中頻振動(dòng)問(wèn)題,在混響實(shí)驗(yàn)室中進(jìn)行某一典型梁板組合結(jié)構(gòu)的聲振試驗(yàn),對(duì)試驗(yàn)獲取的梁和板的加速度響應(yīng)進(jìn)行了分析,基于H-混合模型法建立了有限元-統(tǒng)計(jì)能量(FE-SEA)混合模型,通過(guò)模態(tài)試驗(yàn)對(duì)模型進(jìn)行了驗(yàn)證。對(duì)混合模型計(jì)算的結(jié)果進(jìn)行了分析,并與試驗(yàn)結(jié)果進(jìn)行了對(duì)比。得到了內(nèi)外聲場(chǎng)之間通過(guò)板結(jié)構(gòu)進(jìn)行非共振傳遞和共振傳遞的聲傳遞規(guī)律。FE-SEA混合模型計(jì)算與試驗(yàn)結(jié)果吻合,驗(yàn)證了使用H-混合模型方法進(jìn)行中頻振動(dòng)分析的有效性和可靠性。

    中頻振動(dòng);梁板組合結(jié)構(gòu);FE-SEA建模;聲振分析

    梁板組合結(jié)構(gòu)在航空航天領(lǐng)域中被廣泛使用。由于梁板子結(jié)構(gòu)之間存在較大的動(dòng)態(tài)特性差異,使得在比較寬的頻率范圍內(nèi),整個(gè)結(jié)構(gòu)的振動(dòng)表現(xiàn)為梁的局部低頻振動(dòng)和板的局部高頻振動(dòng)同時(shí)存在的復(fù)雜振動(dòng)情況,這就是典型的結(jié)構(gòu)中頻振動(dòng)[1]。

    中頻振動(dòng)分析相對(duì)于低頻振動(dòng)和高頻振動(dòng)分析而言要復(fù)雜得多,自20世紀(jì)90年代以來(lái),中頻振動(dòng)研究已引起越來(lái)越多的關(guān)注,逐漸成為系統(tǒng)結(jié)構(gòu)動(dòng)態(tài)特性研究領(lǐng)域的熱點(diǎn)之一[2-6]。對(duì)于中頻振動(dòng)分析的方法最早采用的主要是基于Soize[3]提出的結(jié)構(gòu)模糊理論(Structural Fuzzy Theory),此后針對(duì)中頻振動(dòng)問(wèn)題的研究大部分都是在該原理基礎(chǔ)上形成的。目前對(duì)中頻振動(dòng)的計(jì)算方法本質(zhì)上有三類[2]:傳統(tǒng)有限元法、統(tǒng)計(jì)能量法、混合模型法。其中,H-混合模型法與其它計(jì)算方法相比具有更好的理論完整性和更廣泛的工程適用性[7-9],已成為聲振分析軟件VA-One在中頻計(jì)算模塊的理論基礎(chǔ)[10]。然而,現(xiàn)有的對(duì)于中頻結(jié)構(gòu)的研究主要集中在理論計(jì)算方法的拓展上,對(duì)于試驗(yàn)與計(jì)算的完整研究還存在不足,依然需要進(jìn)一步的研究和發(fā)展。

    本文針對(duì)某一典型梁板組合結(jié)構(gòu)的中頻振動(dòng)問(wèn)題,在混響實(shí)驗(yàn)室中進(jìn)行了聲振試驗(yàn),并基于H-混合模型法建立了該梁板組合結(jié)構(gòu)的FE-SEA混合模型。對(duì)試驗(yàn)測(cè)量與模型計(jì)算得到的梁和板的加速度響應(yīng)進(jìn)行了對(duì)比分析,通過(guò)試驗(yàn)與計(jì)算的方法揭示了內(nèi)外聲場(chǎng)之間通過(guò)板結(jié)構(gòu)進(jìn)行非共振傳遞和共振傳遞的聲傳遞規(guī)律,驗(yàn)證了使用H-混合模型方法進(jìn)行中頻振動(dòng)分析的有效性和可靠性。

    1 梁板組合結(jié)構(gòu)中頻振動(dòng)聲振試驗(yàn)研究

    1.1試驗(yàn)系統(tǒng)及過(guò)程

    某一典型的梁板組合結(jié)構(gòu)如圖1所示,整個(gè)結(jié)構(gòu)由19根矩形截面箱梁和7塊厚度不同的板組成。梁之間通過(guò)焊接形成框架結(jié)構(gòu),板通過(guò)螺釘和鉚釘連接在梁上。1~7號(hào)梁長(zhǎng)度為1 m,8~11號(hào)梁長(zhǎng)度為0.8 m,12~13號(hào)梁長(zhǎng)度為0.4 m,16~17號(hào)梁長(zhǎng)度為2 m,18~19號(hào)梁長(zhǎng)度為1.4 m。矩形截面箱梁橫截面尺寸為3 cm×3 cm,厚度為3 mm。7塊板的厚度如表1所示。整個(gè)梁板組合結(jié)構(gòu)長(zhǎng)為2 m,寬為1 m,高為0.8 m,材料均為Q235鋼。

    將該梁板組合系統(tǒng)放入混響室內(nèi)進(jìn)行聲振試驗(yàn)。在每塊板上安裝若干振動(dòng)傳感器,在某一不與板直接連接的梁中部布置振動(dòng)傳感器,使用聲壓計(jì)對(duì)混響室內(nèi)混響聲場(chǎng)和梁板組合結(jié)構(gòu)內(nèi)的封閉空間內(nèi)聲場(chǎng)進(jìn)行聲場(chǎng)測(cè)量。

    圖1 Fig.1 (a) The setup in the reverberation room;(b) The complex built-up structure

    Tab.1 The thickness of each plate

    板編號(hào)板1板2板3板4板5板6板7厚度/mm111.511.51.51

    1.2試驗(yàn)結(jié)果分析

    1.2.1內(nèi)外聲場(chǎng)試驗(yàn)結(jié)果分析

    對(duì)處于混響室內(nèi)的結(jié)構(gòu)施加不同總聲壓級(jí)的混響聲載荷,當(dāng)混響室聲場(chǎng)及結(jié)構(gòu)振動(dòng)達(dá)到穩(wěn)定后,采集結(jié)構(gòu)的內(nèi)外聲場(chǎng)和振動(dòng)加速度響應(yīng)。圖2所示為試驗(yàn)中施加的6個(gè)不同總聲壓級(jí)的混響聲場(chǎng),總聲壓級(jí)處于133 dB~147.6 dB之間。

    圖3所示為不同總聲壓級(jí)的混響聲場(chǎng)作用下結(jié)構(gòu)的內(nèi)聲場(chǎng)。內(nèi)聲場(chǎng)的總聲壓級(jí)處于124 dB~137.4 dB之間。不同混響聲場(chǎng)載荷下的結(jié)構(gòu)內(nèi)聲場(chǎng)聲譜具有相似性,并且總聲壓級(jí)隨著外部混響聲場(chǎng)總聲壓級(jí)增加而增加,六種不同混響聲場(chǎng)下內(nèi)外聲場(chǎng)總聲壓級(jí)差值處于8.5 dB~10.2 dB之間。在所有的混響聲場(chǎng)條件中,4 000 Hz以上無(wú)明顯峰值,而結(jié)構(gòu)內(nèi)聲場(chǎng)聲譜在頻率為9 200 Hz和13 700 Hz處均產(chǎn)生了明顯的峰值。這說(shuō)明在以這兩個(gè)頻率為中心頻率的頻帶內(nèi),結(jié)構(gòu)內(nèi)聲場(chǎng)對(duì)外部混響聲場(chǎng)的吸收效率較高。

    圖2 試驗(yàn)中施加的不同總聲壓級(jí)的混響聲場(chǎng)Fig.2 Different diffuseacoustic field loading conditions in the testing

    圖3 不同總聲壓級(jí)的混響聲場(chǎng)作用下結(jié)構(gòu)的內(nèi)聲場(chǎng)Fig.3 The internal acoustic fields with different external diffuse acoustic fields

    1.2.2梁板振動(dòng)試驗(yàn)結(jié)果分析

    對(duì)混響聲場(chǎng)總聲壓級(jí)為147.6 dB試驗(yàn)條件下的梁板組合結(jié)構(gòu)振動(dòng)試驗(yàn)結(jié)果進(jìn)行分析,圖4所示為該試驗(yàn)條件下的梁板振動(dòng)功率譜密度,圖4(a)~圖4(g)分別為1~7號(hào)板上各測(cè)點(diǎn)的振動(dòng)功率譜密度。同一塊板上不同測(cè)點(diǎn)的功率譜密度曲線表現(xiàn)出相似性,并且總均方根值一致,這是由于薄板在高頻振動(dòng)時(shí),波長(zhǎng)較短,邊界條件對(duì)板的能量分布影響較小,從而使得板上除邊界附近的其它各處振動(dòng)均方根值一致,因此可以使用能量來(lái)描述整個(gè)板的振動(dòng)狀態(tài)。在統(tǒng)計(jì)能量分析中,子系統(tǒng)的響應(yīng)是指在時(shí)間、頻率和空間上的平均響應(yīng),因此對(duì)于梁板組合結(jié)構(gòu)中的板件振動(dòng),適合使用統(tǒng)計(jì)能量分析方法。

    圖4(h)所示為圖1(b)中梁上的加速度測(cè)點(diǎn)的功率譜密度曲線。可以看出頻率峰值非常明顯,在200 Hz以內(nèi)有四個(gè)明顯的峰值,分別是56 Hz、74 Hz、84 Hz、134 Hz。這說(shuō)明梁板混合結(jié)構(gòu)中梁的模態(tài)密度較小,梁在頻率和空間平均意義上的平均響應(yīng)并不能代表其真實(shí)的振動(dòng)響應(yīng)情況,因此對(duì)于模態(tài)密度較小的梁,統(tǒng)計(jì)能量分析則不再適用,一般使用有限元法進(jìn)行分析。

    (a)(b)(c)

    (d)(e)(f)

    (g)(h)圖4 混響聲場(chǎng)總聲壓級(jí)為147.6dB試驗(yàn)條件下的梁板振動(dòng)功率譜密度曲線;(a)~(g)分別為1~7號(hào)板上各測(cè)點(diǎn)的振動(dòng)功率譜密度曲線;(h)為梁上測(cè)點(diǎn)的功率譜密度曲線Fig.4Theaccelerationpowerdensityspectrumsofthe

    2 梁板組合結(jié)構(gòu)中頻振動(dòng)的FE-SEA混合模型計(jì)算

    2.1混合模型的建立與驗(yàn)證

    2.1.1混合模型建立

    通過(guò)梁板組合結(jié)構(gòu)中的子系統(tǒng)間的波長(zhǎng)比,可以判斷梁子結(jié)構(gòu)和板子結(jié)構(gòu)之間存在的動(dòng)態(tài)特性差異,從而能夠確定梁板組合結(jié)構(gòu)的中頻振動(dòng)頻率范圍。對(duì)于梁子結(jié)構(gòu),其橫向彎曲振動(dòng)的波長(zhǎng)與材料屬性和橫截面形狀有關(guān),對(duì)于板子結(jié)構(gòu),其橫向彎曲振動(dòng)的波長(zhǎng)與材料屬性和厚度有關(guān),因此,這里選取厚度不同的6號(hào)板和7號(hào)板子結(jié)構(gòu)與梁子結(jié)構(gòu)的波長(zhǎng)進(jìn)行比較,如圖5所示,可以看出相同頻率下1 mm板的彎曲振動(dòng)波長(zhǎng)比1.5 mm板的彎曲波長(zhǎng)小,梁的彎曲波長(zhǎng)與兩個(gè)厚度的板的彎曲波長(zhǎng)之比分別為5和6。說(shuō)明梁子結(jié)構(gòu)與各板子結(jié)構(gòu)動(dòng)態(tài)差異非常大,適合采用中頻分析方法。

    圖5 梁與6號(hào)板和7號(hào)板的波長(zhǎng)及波長(zhǎng)比曲線Fig.5 Structural wavelengths andwavelength ratios

    使用H-混合模型法(Hybrid FE-SEA Method)對(duì)該系統(tǒng)建模。該梁板組合系統(tǒng)整體結(jié)構(gòu)可以分為確定性子結(jié)構(gòu)和非確定性子結(jié)構(gòu)。其中梁結(jié)構(gòu)較為簡(jiǎn)單,動(dòng)態(tài)分析中模態(tài)密度較小,認(rèn)為是確定性子結(jié)構(gòu),用有限元子系統(tǒng)描述。板結(jié)構(gòu)由于在動(dòng)態(tài)分析中包含大量高階模態(tài),同時(shí)板連接邊界參數(shù)具有一定程度的隨機(jī)不確定性,因此認(rèn)為是非確定性子結(jié)構(gòu),用統(tǒng)計(jì)能量子系統(tǒng)描述。梁板組合形成一個(gè)封閉的聲腔,板的振動(dòng)會(huì)改變聲腔的聲壓,而聲腔聲壓變化又會(huì)激勵(lì)板產(chǎn)生振動(dòng),從而形成了一個(gè)結(jié)構(gòu)與空氣相互作用的聲固耦合系統(tǒng)。

    圖6所示為該梁板組合結(jié)構(gòu)的有限元-統(tǒng)計(jì)能量(FE-SEA)混合模型。該模型中7塊板使用結(jié)構(gòu)統(tǒng)計(jì)能量子系統(tǒng)描述,梁使用結(jié)構(gòu)有限元子系統(tǒng)描述,內(nèi)聲腔使用聲統(tǒng)計(jì)能量子系統(tǒng)描述。每塊板上的混響聲載荷用擴(kuò)散聲場(chǎng)(DAF)載荷描述。使用半無(wú)限流體(SIF)子系統(tǒng)來(lái)描述結(jié)構(gòu)表面對(duì)空氣的聲輻射效應(yīng)。求解在混響聲載荷作用下的結(jié)構(gòu)響應(yīng)和內(nèi)聲場(chǎng)響應(yīng)。

    混合模型中梁和板的材料密度為7 870 kg/m3,彈性模量為 2.11×1011Pa,泊松比為0.29。使用基于脈沖激勵(lì)的內(nèi)損耗因子獲取方法[11]得到梁板組合狀態(tài)下1~2號(hào)板、5~7號(hào)板以及梁主結(jié)構(gòu)的內(nèi)損耗因子為0.45%?;祉懯以囼?yàn)中,3號(hào)和4號(hào)板與支撐結(jié)構(gòu)接觸,內(nèi)損耗因子為1%。內(nèi)聲場(chǎng)平均聲吸收系數(shù)為2%。

    圖6 FE-SEA混合模型Fig.6 The hybrid FE-SEA model

    2.1.2模型驗(yàn)證

    對(duì)梁板組合結(jié)構(gòu)進(jìn)行模態(tài)試驗(yàn),試驗(yàn)過(guò)程中采用橡皮繩豎直懸吊模擬自由邊界條件,如圖7所示。使用力錘法進(jìn)行主結(jié)構(gòu)的模態(tài)分析,不考慮結(jié)構(gòu)中板上的局部模態(tài),模態(tài)試驗(yàn)獲取的前四階模態(tài)如表2所示。

    將建立的FE-SEA混合模型中的統(tǒng)計(jì)能量子系統(tǒng)劃分網(wǎng)格,進(jìn)行模態(tài)求解,并與模態(tài)試驗(yàn)分析結(jié)果進(jìn)行對(duì)比,如圖8所示為前四階模態(tài)振型求解結(jié)果,模態(tài)頻率分別為60 Hz、72 Hz、85 Hz、140 Hz。與模態(tài)試驗(yàn)結(jié)果相符,說(shuō)明建立的FE-SEA混合模型能夠反映梁板組合結(jié)構(gòu)的主要?jiǎng)恿W(xué)特性。

    圖7 梁板組合結(jié)構(gòu)模態(tài)試驗(yàn)Fig.7 Experimental modal test of the complex built-up structure

    模態(tài)階頻率/Hz阻尼/%振型描述1600.5主結(jié)構(gòu)彎2720.3主結(jié)構(gòu)彎3850.2主結(jié)構(gòu)扭41400.1局部扭

    圖8 前四階模態(tài)振型Fig.8 The mode shapes for

    2.2計(jì)算結(jié)果分析

    在總聲壓級(jí)為147.6 dB的混響室聲場(chǎng)條件下對(duì)模型進(jìn)行計(jì)算,模型中的擴(kuò)散聲場(chǎng)載荷指的是作用在結(jié)構(gòu)子系統(tǒng)表面上的壓力譜,在典型混響室中,結(jié)構(gòu)表面壓力比混響室內(nèi)聲場(chǎng)高3 dB[12],如圖9所示。在該擴(kuò)散聲場(chǎng)載荷下,模型計(jì)算的結(jié)構(gòu)內(nèi)聲場(chǎng)總聲壓級(jí)和各子系統(tǒng)的振動(dòng)加速度均方根值分布如圖10所示。內(nèi)聲場(chǎng)總聲壓級(jí)為139.4 dB,結(jié)構(gòu)振動(dòng)的最大加速度均方根值為34.41 g。

    圖9 模型中施加的擴(kuò)散聲場(chǎng)載荷譜Fig.9 The diffuse acoustic field loading

    圖10 模型計(jì)算結(jié)果Fig.10 Calculated results from the hybrid FE-SEA model

    2.2.1內(nèi)聲場(chǎng)計(jì)算結(jié)果分析

    結(jié)構(gòu)內(nèi)聲場(chǎng)的試驗(yàn)測(cè)量與計(jì)算結(jié)果圖11所示。計(jì)算結(jié)果與試驗(yàn)值在整個(gè)頻段內(nèi)都吻合的較好,總聲壓級(jí)與試驗(yàn)測(cè)量結(jié)果相差2 dB,能夠反應(yīng)出結(jié)構(gòu)與聲場(chǎng)響應(yīng)的整體變化趨勢(shì)。計(jì)算的結(jié)構(gòu)內(nèi)聲場(chǎng)聲譜在8 000 Hz中心頻率和11 986.5 Hz中心頻率處產(chǎn)生明顯峰值。

    圖11 結(jié)構(gòu)內(nèi)聲場(chǎng)的試驗(yàn)與計(jì)算結(jié)果Fig.11 Comparison of internal acoustic fields obtained by simulation and experiment

    對(duì)于外聲場(chǎng)激勵(lì)板結(jié)構(gòu)振動(dòng)產(chǎn)生內(nèi)聲場(chǎng)問(wèn)題,內(nèi)外聲場(chǎng)之間的聲傳遞可以分為兩方面,一方面是非共振傳遞,另一方面是共振傳遞。圖12所示為該模型中通過(guò)非共振傳遞方式和共振傳遞方式輸入到內(nèi)聲場(chǎng)的功率曲線,圖中DAF表示通過(guò)各板的非共振傳遞功率,Plate表示通過(guò)各板的共振傳遞功率,Total表示總的傳入內(nèi)聲場(chǎng)的功率。在6 350 Hz范圍以內(nèi),內(nèi)聲場(chǎng)主要來(lái)源于混響聲載荷的非共振傳遞,而在6 350 Hz以上頻率范圍,內(nèi)聲場(chǎng)則主要來(lái)源于共振傳遞通方式。通過(guò)1.5 mm厚的3、5和6號(hào)板的共振傳遞峰值位于8 000 Hz,而通過(guò)1 mm厚的1、2、4和7號(hào)板的共振傳遞峰值位于11 986.5 Hz,這說(shuō)明共振傳遞與板的厚度相關(guān)。

    當(dāng)受到擴(kuò)散聲場(chǎng)激勵(lì)時(shí),強(qiáng)迫彎曲波控制著板的聲傳遞特性。板子系統(tǒng)存在一個(gè)臨界頻率,在臨界頻率以下,聲傳遞以非共振傳遞為主。臨界頻率及臨界頻率以上,聲驅(qū)動(dòng)板產(chǎn)生共振,此時(shí)聲傳遞以共振傳遞為主。臨界頻率是指當(dāng)板的彎曲波速cb等于聲波波速ca時(shí)所對(duì)應(yīng)的頻率。

    二維板的彎曲剛度B及臨界頻率為:

    (1)

    (2)

    式(2)中E,h,ν,m分別為板的彈性模量、厚度、泊松比和面密度(kg/m2),代入材料參數(shù)和厚度計(jì)算得到1 mm和1.5 mm厚板的臨界頻率分別為8 000 Hz和118 000 Hz,分別處于1/12倍頻程中以8 000 Hz和11 986.5 Hz為中心頻率的頻帶范圍內(nèi)。試驗(yàn)和計(jì)算的內(nèi)聲場(chǎng)聲譜在高頻范圍內(nèi)的兩個(gè)突出峰值對(duì)應(yīng)的頻率分別是兩種厚度板的臨界頻率。在臨界頻率處,板內(nèi)彎曲波速與聲波波速相等,發(fā)生共振,板對(duì)內(nèi)聲場(chǎng)的聲輻射效率極高,從而形成內(nèi)聲場(chǎng)輸入功率的峰值。因此,內(nèi)聲場(chǎng)試驗(yàn)值在9 200 Hz和13 700 Hz處的兩個(gè)峰值與計(jì)算結(jié)果中的8 000 Hz和11 986.5 Hz處的峰值是對(duì)應(yīng)的,是由于內(nèi)外聲場(chǎng)在兩種厚度板的臨界頻率上的共振傳遞所引起的。由于結(jié)構(gòu)實(shí)際的幾何參數(shù)、尺寸公差等與其名義參數(shù)存在一定誤差,而系統(tǒng)在高頻范圍內(nèi)響應(yīng)對(duì)這些參數(shù)的微小變化都非常敏感,從而造成了兩個(gè)峰值的試驗(yàn)與計(jì)算結(jié)果存在一定的差異。

    圖12 通過(guò)非共振傳遞方式和共振傳遞方式輸入到內(nèi)聲場(chǎng)的功率,DAF示非共振傳遞,Plate表示共振傳遞,Total表示總輸入功率Fig.12 Power input to the internal acoustic field (represented with Total) due to the non-resonant transmission (represented with DAF) and resonant transmission (represented with Plate)

    2.2.2梁板計(jì)算結(jié)果分析

    通過(guò)FE-SEA混合模型計(jì)算的各板加速度響應(yīng)與試驗(yàn)測(cè)量結(jié)果比較如圖13所示。圖13(a)~(g)分別為七塊板加速度響應(yīng)的功率譜密度曲線比較,各板計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合,計(jì)算的功率譜密度曲線能夠反映試驗(yàn)測(cè)量結(jié)果的主要峰值及變化趨勢(shì)。圖13(h)為各板響應(yīng)總均方根值的模型計(jì)算與試驗(yàn)結(jié)果比較,最大誤差不超過(guò)2 dB。

    在梁板組合結(jié)構(gòu)的聲振問(wèn)題中,由于薄板結(jié)構(gòu)模態(tài)密集,使用傳統(tǒng)有限元法存在模態(tài)參數(shù)具有不確定性以及計(jì)算量大的問(wèn)題。而使用統(tǒng)計(jì)能量子系統(tǒng)表征薄板,建立梁板組合結(jié)構(gòu)中頻振動(dòng)的FE-SEA混合模型,能夠有效可靠地得到結(jié)構(gòu)上的響應(yīng),與傳統(tǒng)有限元法相比,減少了計(jì)算量,提高了計(jì)算效率。

    梁板組合結(jié)構(gòu)中梁上加速度傳感器的測(cè)量結(jié)果與計(jì)算結(jié)果對(duì)比如圖14所示。計(jì)算與試驗(yàn)測(cè)量的RMS值基本一致,計(jì)算結(jié)果能夠體現(xiàn)試驗(yàn)測(cè)量結(jié)果的主要峰值和變化趨勢(shì)。結(jié)構(gòu)振動(dòng)頻率越高,主結(jié)構(gòu)對(duì)打孔、螺栓和裝配等因素越敏感,造成了計(jì)算結(jié)果中峰值位置相對(duì)試驗(yàn)測(cè)量結(jié)果在頻率上有一些平移。

    3 結(jié) 論

    針對(duì)某一典型梁板組合結(jié)構(gòu)的中頻振動(dòng)問(wèn)題,在混響實(shí)驗(yàn)室中進(jìn)行了聲振試驗(yàn),并基于H-混合模型法建立了該梁板組合結(jié)構(gòu)的FE-SEA混合模型。對(duì)試驗(yàn)測(cè)量與模型計(jì)算得到的梁和板的加速度響應(yīng)進(jìn)行了對(duì)比分析, 通過(guò)試驗(yàn)與計(jì)算的方法揭示了內(nèi)外聲場(chǎng)之間通過(guò)板結(jié)構(gòu)進(jìn)行非共振傳遞和共振傳遞的聲傳遞規(guī)律,由試驗(yàn)與模型獲取的內(nèi)外聲場(chǎng)聲傳遞規(guī)律的一致性,驗(yàn)證了使用H-混合模型方法進(jìn)行中頻振動(dòng)分析的有效性和可靠性。得到了如下結(jié)論:

    (a)(b)(c)

    (d)(e)(f)

    (g)(h)圖13 FE-SEA混合模型計(jì)算的板響應(yīng)結(jié)果與試驗(yàn)結(jié)果比較Fig.13ComparisonofaccelerationresponsesobtainedbyhybridFE-SEAmethodandexperiment

    圖14 梁響應(yīng)的計(jì)算與試驗(yàn)結(jié)果比較Fig.14 Comparison of beam responses obtained by simulation and experiment

    (1) 對(duì)于外聲場(chǎng)激勵(lì)平板振動(dòng)產(chǎn)生的內(nèi)聲場(chǎng),內(nèi)外聲場(chǎng)之間的聲傳遞方式可以分為非共振傳遞和共振傳遞兩種方式。在臨界頻率處,板內(nèi)彎曲波速與聲波波速相等,發(fā)生共振,板對(duì)內(nèi)聲場(chǎng)的聲輻射效率極高。

    (2) 在梁板組合結(jié)構(gòu)中,梁子結(jié)構(gòu)與板子結(jié)構(gòu)的波長(zhǎng)比往往比較大,說(shuō)明梁子結(jié)構(gòu)與板子結(jié)構(gòu)存在著較大的動(dòng)態(tài)特性差異,使用H-混合模型法建模求解,能夠擴(kuò)大求解的頻率范圍,減少計(jì)算量,提高計(jì)算效率。

    (3) FE-SEA混合模型計(jì)算與試驗(yàn)測(cè)量的各板響應(yīng)均方根值吻合,內(nèi)聲場(chǎng)預(yù)示結(jié)果與試驗(yàn)值在整個(gè)頻段內(nèi)一致,能夠反應(yīng)出結(jié)構(gòu)與聲場(chǎng)響應(yīng)整體的變化趨勢(shì),試驗(yàn)與模型獲取的內(nèi)外聲場(chǎng)聲傳遞規(guī)律的一致性驗(yàn)證了使用H-混合模型方法進(jìn)行中頻振動(dòng)響應(yīng)預(yù)示的有效性和可靠性。

    [1] 紀(jì)琳.中頻振動(dòng)分析方法[M]. 北京:機(jī)械工業(yè)出版社,2013.

    [2] Desmet W. Mid-frequency vibro-acoustic modelling: challenges and potential solutions[C]// Proceedings of ISMA, 2002, 2:835-862.

    [3] Soize C.A model and numerical method in the medium frequency range for vibroacoustic predictions using the theory of structural fuzzy[J].The Journal of the Acoustical Society of America, 1993, 94(2):849-865.

    [4] Zhao X, Vlahopoulos N. A basic hybrid finite element formulation for mid-frequency analysis of beams connected at an arbitrary angle[J]. Journal of Sound and Vibration, 2004, 269(6): 135-164.

    [5] Shorter P J, Langley R S.Vibro-acoustic analysis of complex systems[J].Journal of Sound and Vibration, 2005, 288(3): 669-699.

    [6] Langley R S.Recent advances and remaining challenges in the statistical energy analysis of dynamics systems[C]//Proceedings of the 7th European Conference on Structural Dynamics, Southampton:Southampton University, 2008.

    [7] Langley R S, Bremner P. A hybrid method for the vibration analysis of complex structural-axoustic systems[J]. Journal of the Acoustical Society of America,1999,105(3):1657-1671.

    [8] Shorter P J,Gardner B K,Bremner P G. A hybrid method for full spectrum noise and vibration prediction[J]. Journal of Computational Acoustics,2003,11(2):323-338.

    [9] Shorter P J, Langley R S. On the reciprocity relationship between direct field radiation and diffuse reverberant loading[J]. Journal of the Acoustical Society of America, 2005, 117(1):85-95.

    [10] VA One 2010.5 User’s Guide[M]. ESI Group, 2011.

    [11] 張永杰, 肖健, 韋冰峰,等. 一種新的基于脈沖激勵(lì)的內(nèi)損耗因子獲取方法[J]. 振動(dòng)與沖擊,2014,33(12):161-164.

    ZHANG Yong-jie,XIAO Jian,WEI Bing-feng,et al. A new method for estimating damping loss factor based on multi-pulse excitation[J]. Jounal of Vibration and Shock, 2014, 33(12): 161-164.

    [12] Waterhouse R V. Interference patterns in reverberant sound fields[J]. The Journal of the Acoustical Society of America, 1955, 27(2): 247-258.

    Tests and simulation for vibration in mid-frequency domain of a beam-plate built-up structure

    ZHANG Yong-jie, CHEN Lu, WEI Bing-feng, WANG Zhe, QIN Zhao-hong

    (Beijing Institute of Structure and Environment Engineering, Beijing 100076, China)

    To investigate vibration problems of beam-plate built-up structure in mid-frequency domain, vibro-acoustic tests were performed in a reverberation lab. The measured acceleration responses of the structure were analyzed to establish a hybrid FE-SEA model. The hybrid FE-SEA model’s simulation results agreed well with test ones. The non-resonant and resonant sound transmission rules between external acoustic field and internal one through the structure were then deduced. The results showed that the hybrid FE-SEA model is reliable and effective for vibration analysis in mid-frequency domain of a beam-plate built-up structure.

    beam-plate built-up structure; mid-frequency vibration; hybrid FE-SEA model; vibro-acoustic analysis

    10.13465/j.cnki.jvs.2016.13.031

    2015-04-20修改稿收到日期:2015-07-04

    張永杰 男,博士,高級(jí)工程師,1983年3月生

    TB535

    A

    猜你喜歡
    梁板聲壓級(jí)聲場(chǎng)
    機(jī)器噪聲平均聲壓級(jí)計(jì)算方法差異性實(shí)證研究
    淺談房建施工中現(xiàn)澆梁板模板施工技術(shù)
    淺析房建施工中現(xiàn)澆梁板模板的要求、質(zhì)量控制及其意義
    一種計(jì)算消聲室聲壓級(jí)的新方法
    基于BIM的鐵路車站聲場(chǎng)仿真分析研究
    全新DXR mkll有源揚(yáng)聲器
    演藝科技(2019年4期)2019-03-30 03:21:46
    探尋360°全聲場(chǎng)發(fā)聲門道
    房建施工中現(xiàn)澆梁板模板的要求及質(zhì)量控制
    江西建材(2018年4期)2018-04-10 12:37:08
    Diodes1.9W D類音頻放大器提供高聲壓級(jí)水平并延長(zhǎng)電池壽命
    板結(jié)構(gòu)-聲場(chǎng)耦合分析的FE-LSPIM/FE法
    日本与韩国留学比较| 午夜福利乱码中文字幕| 女性被躁到高潮视频| 一级片'在线观看视频| 久久av网站| 18禁裸乳无遮挡动漫免费视频| 国产极品天堂在线| 成人免费观看视频高清| 亚洲国产精品国产精品| 婷婷色综合www| 国产成人精品无人区| 中文字幕亚洲精品专区| 热99国产精品久久久久久7| 久久久久久久精品精品| av福利片在线| 亚洲成人手机| 又粗又硬又长又爽又黄的视频| 欧美精品一区二区免费开放| 九草在线视频观看| 欧美3d第一页| 丝瓜视频免费看黄片| 又黄又爽又刺激的免费视频.| 日韩精品免费视频一区二区三区 | 人妻一区二区av| 精品人妻偷拍中文字幕| av又黄又爽大尺度在线免费看| 美女内射精品一级片tv| 中文字幕免费在线视频6| 欧美+日韩+精品| 丝袜人妻中文字幕| 免费av中文字幕在线| 免费日韩欧美在线观看| 午夜福利视频在线观看免费| 制服丝袜香蕉在线| 18禁裸乳无遮挡动漫免费视频| 咕卡用的链子| 男女午夜视频在线观看 | 久久青草综合色| 高清av免费在线| 成人毛片60女人毛片免费| av国产久精品久网站免费入址| 欧美日韩精品成人综合77777| 精品国产一区二区久久| 高清黄色对白视频在线免费看| 免费av不卡在线播放| 亚洲国产毛片av蜜桃av| 激情视频va一区二区三区| 飞空精品影院首页| 亚洲久久久国产精品| 蜜臀久久99精品久久宅男| 国国产精品蜜臀av免费| 黄色视频在线播放观看不卡| 国产色婷婷99| 母亲3免费完整高清在线观看 | 桃花免费在线播放| 观看美女的网站| 久久久精品区二区三区| 久久av网站| 两性夫妻黄色片 | 国产一级毛片在线| 91成人精品电影| 精品99又大又爽又粗少妇毛片| 欧美精品av麻豆av| 精品亚洲成a人片在线观看| 伊人久久国产一区二区| 成年av动漫网址| 国产乱人偷精品视频| 日日啪夜夜爽| 午夜福利视频在线观看免费| 午夜视频国产福利| 久热这里只有精品99| 欧美老熟妇乱子伦牲交| 精品人妻一区二区三区麻豆| 免费看av在线观看网站| 午夜精品国产一区二区电影| 国产精品.久久久| 少妇被粗大猛烈的视频| 国产免费视频播放在线视频| 涩涩av久久男人的天堂| 99国产综合亚洲精品| 国产成人91sexporn| 丝袜美足系列| 国产日韩欧美在线精品| 成人毛片60女人毛片免费| 亚洲欧洲日产国产| 亚洲五月色婷婷综合| 另类亚洲欧美激情| 午夜日本视频在线| 免费av中文字幕在线| 99热全是精品| www.av在线官网国产| 国产极品粉嫩免费观看在线| 日韩免费高清中文字幕av| 久久久国产欧美日韩av| 久久人妻熟女aⅴ| 大片电影免费在线观看免费| 欧美成人午夜精品| 国产伦理片在线播放av一区| 又大又黄又爽视频免费| 国产免费一区二区三区四区乱码| 国产成人精品一,二区| 日本黄大片高清| 久久久国产一区二区| 亚洲第一av免费看| 精品熟女少妇av免费看| 色吧在线观看| 狂野欧美激情性bbbbbb| 尾随美女入室| 精品一区二区三区四区五区乱码 | 久久精品国产鲁丝片午夜精品| 国产熟女欧美一区二区| 777米奇影视久久| a级毛片黄视频| 亚洲国产看品久久| 亚洲精品乱码久久久久久按摩| 中文字幕最新亚洲高清| 如何舔出高潮| 日韩,欧美,国产一区二区三区| 少妇被粗大猛烈的视频| 性色avwww在线观看| 一级爰片在线观看| 国产精品 国内视频| 制服人妻中文乱码| 亚洲精品国产av蜜桃| 波野结衣二区三区在线| 亚洲成国产人片在线观看| 久久久久网色| 欧美人与善性xxx| 亚洲欧洲国产日韩| 一本久久精品| 少妇被粗大的猛进出69影院 | 在线 av 中文字幕| 精品视频人人做人人爽| 2022亚洲国产成人精品| 国产成人免费无遮挡视频| videos熟女内射| 日本黄大片高清| 久久ye,这里只有精品| 妹子高潮喷水视频| 捣出白浆h1v1| 五月玫瑰六月丁香| 久久精品夜色国产| 亚洲在久久综合| 极品人妻少妇av视频| 国产精品熟女久久久久浪| 精品人妻在线不人妻| 一区二区三区精品91| 亚洲色图综合在线观看| 乱码一卡2卡4卡精品| 一级片免费观看大全| 久久热在线av| 午夜福利视频在线观看免费| 国产精品无大码| 综合色丁香网| 欧美97在线视频| 婷婷色麻豆天堂久久| 校园人妻丝袜中文字幕| 91成人精品电影| 男女啪啪激烈高潮av片| 欧美精品av麻豆av| 99九九在线精品视频| 一区二区三区乱码不卡18| 日本av手机在线免费观看| 精品一品国产午夜福利视频| 日韩成人伦理影院| 国产免费现黄频在线看| 婷婷色麻豆天堂久久| 七月丁香在线播放| 国产片内射在线| 精品久久久精品久久久| 成人无遮挡网站| 最近手机中文字幕大全| 97人妻天天添夜夜摸| 国产色爽女视频免费观看| 精品酒店卫生间| 日本与韩国留学比较| av一本久久久久| 一二三四在线观看免费中文在 | 两性夫妻黄色片 | 国产在线免费精品| 日韩三级伦理在线观看| 欧美成人午夜精品| 大陆偷拍与自拍| 午夜久久久在线观看| 久久国内精品自在自线图片| 另类精品久久| 国产无遮挡羞羞视频在线观看| av女优亚洲男人天堂| 国产熟女午夜一区二区三区| 波多野结衣一区麻豆| 久久精品国产自在天天线| 一级片'在线观看视频| 国精品久久久久久国模美| 又黄又粗又硬又大视频| 激情视频va一区二区三区| 亚洲国产毛片av蜜桃av| 日韩,欧美,国产一区二区三区| 久久久久久伊人网av| 国产 精品1| 国产在线视频一区二区| 国产国拍精品亚洲av在线观看| 啦啦啦啦在线视频资源| 精品国产一区二区三区久久久樱花| 夫妻午夜视频| 精品亚洲成a人片在线观看| 亚洲av国产av综合av卡| av卡一久久| 亚洲丝袜综合中文字幕| 看非洲黑人一级黄片| 国产精品嫩草影院av在线观看| 精品少妇黑人巨大在线播放| 各种免费的搞黄视频| 妹子高潮喷水视频| 男女免费视频国产| 久久99热这里只频精品6学生| 日本91视频免费播放| 男女边摸边吃奶| 国产成人免费无遮挡视频| 午夜福利乱码中文字幕| 日本wwww免费看| 亚洲欧美中文字幕日韩二区| 久久鲁丝午夜福利片| 高清av免费在线| 久久久久久人人人人人| 永久免费av网站大全| 91在线精品国自产拍蜜月| 精品国产一区二区三区四区第35| 久久久久精品久久久久真实原创| 人妻系列 视频| 欧美日韩综合久久久久久| 久久热在线av| 赤兔流量卡办理| 草草在线视频免费看| 一级黄片播放器| 99视频精品全部免费 在线| www.熟女人妻精品国产 | 美女视频免费永久观看网站| 欧美xxⅹ黑人| 日韩欧美精品免费久久| 在线精品无人区一区二区三| 多毛熟女@视频| 日韩制服丝袜自拍偷拍| 免费人妻精品一区二区三区视频| 丝袜脚勾引网站| 日韩,欧美,国产一区二区三区| 我的女老师完整版在线观看| 啦啦啦啦在线视频资源| 天天操日日干夜夜撸| 亚洲成人手机| 欧美激情 高清一区二区三区| 爱豆传媒免费全集在线观看| 九九爱精品视频在线观看| 亚洲精品自拍成人| 欧美人与性动交α欧美软件 | 大片免费播放器 马上看| 久热这里只有精品99| 亚洲综合色惰| 午夜精品国产一区二区电影| 国产精品一国产av| 女性被躁到高潮视频| 成人亚洲精品一区在线观看| 国产av国产精品国产| 日韩一区二区视频免费看| 国产精品秋霞免费鲁丝片| 99热6这里只有精品| 久久人妻熟女aⅴ| 国精品久久久久久国模美| 热99久久久久精品小说推荐| 各种免费的搞黄视频| 三上悠亚av全集在线观看| 男女啪啪激烈高潮av片| 日本欧美国产在线视频| 亚洲人成77777在线视频| 成人18禁高潮啪啪吃奶动态图| 国产色婷婷99| 爱豆传媒免费全集在线观看| 九九爱精品视频在线观看| 日本与韩国留学比较| 9热在线视频观看99| 免费av不卡在线播放| 黄色怎么调成土黄色| 赤兔流量卡办理| 精品午夜福利在线看| 纯流量卡能插随身wifi吗| 日产精品乱码卡一卡2卡三| 老司机影院成人| 色94色欧美一区二区| 久久影院123| 国产又爽黄色视频| 精品久久久久久电影网| 纯流量卡能插随身wifi吗| 久久99精品国语久久久| 51国产日韩欧美| 国产免费福利视频在线观看| 99久久综合免费| 国产毛片在线视频| av一本久久久久| 国产色婷婷99| 啦啦啦在线观看免费高清www| 内地一区二区视频在线| 亚洲美女黄色视频免费看| 最近2019中文字幕mv第一页| 边亲边吃奶的免费视频| 人妻人人澡人人爽人人| 国产一区二区激情短视频 | 三级国产精品片| 最近最新中文字幕大全免费视频 | 久久久久久久久久久久大奶| 欧美 亚洲 国产 日韩一| 五月开心婷婷网| 成人漫画全彩无遮挡| 精品视频人人做人人爽| 亚洲av福利一区| av线在线观看网站| 五月伊人婷婷丁香| 成人二区视频| 欧美人与性动交α欧美精品济南到 | 最新的欧美精品一区二区| 2021少妇久久久久久久久久久| 岛国毛片在线播放| 人人澡人人妻人| 国产av一区二区精品久久| 在线观看三级黄色| 国产亚洲精品第一综合不卡 | 成年动漫av网址| 亚洲国产精品专区欧美| 中文字幕制服av| 国产黄频视频在线观看| 麻豆精品久久久久久蜜桃| 欧美精品一区二区大全| 热99国产精品久久久久久7| 国产av码专区亚洲av| 在线观看免费高清a一片| 亚洲,欧美,日韩| 国产亚洲最大av| 青春草视频在线免费观看| 99热网站在线观看| 久久精品国产亚洲av天美| 99热国产这里只有精品6| 日韩 亚洲 欧美在线| av视频免费观看在线观看| 国产精品偷伦视频观看了| 免费观看性生交大片5| 免费黄网站久久成人精品| 国产一区二区激情短视频 | 欧美变态另类bdsm刘玥| 免费人妻精品一区二区三区视频| 91精品伊人久久大香线蕉| av在线app专区| 久久99一区二区三区| 欧美性感艳星| 亚洲精品成人av观看孕妇| 国产不卡av网站在线观看| av在线观看视频网站免费| 一区二区三区四区激情视频| 国产免费又黄又爽又色| 91国产中文字幕| 春色校园在线视频观看| 久久av网站| 国产一区亚洲一区在线观看| 97精品久久久久久久久久精品| 免费av不卡在线播放| 国产精品一区二区在线不卡| 在线观看人妻少妇| 国产精品成人在线| 亚洲综合精品二区| 宅男免费午夜| 亚洲综合色惰| 另类精品久久| 天天操日日干夜夜撸| 在线观看美女被高潮喷水网站| 一区二区日韩欧美中文字幕 | 亚洲精品国产av蜜桃| 最黄视频免费看| 久久久久久人人人人人| 亚洲欧美成人综合另类久久久| 九九爱精品视频在线观看| 搡女人真爽免费视频火全软件| 久久久久久久大尺度免费视频| 婷婷色综合大香蕉| 久久久久久久大尺度免费视频| 欧美人与性动交α欧美精品济南到 | 大香蕉97超碰在线| 亚洲第一区二区三区不卡| 亚洲综合精品二区| 亚洲av日韩在线播放| 国产精品久久久久久精品电影小说| 欧美精品一区二区大全| 亚洲欧美日韩另类电影网站| 一区二区三区精品91| 欧美3d第一页| 最新中文字幕久久久久| av在线app专区| 18禁裸乳无遮挡动漫免费视频| 国产一区有黄有色的免费视频| 少妇猛男粗大的猛烈进出视频| 十八禁网站网址无遮挡| 91精品国产国语对白视频| 日韩,欧美,国产一区二区三区| 中文字幕av电影在线播放| 亚洲国产日韩一区二区| 卡戴珊不雅视频在线播放| 免费高清在线观看日韩| 大香蕉97超碰在线| 22中文网久久字幕| 秋霞在线观看毛片| 91aial.com中文字幕在线观看| 欧美3d第一页| 一区二区av电影网| 乱码一卡2卡4卡精品| 精品久久蜜臀av无| 久久精品熟女亚洲av麻豆精品| 欧美另类一区| 不卡视频在线观看欧美| 精品少妇内射三级| 欧美xxⅹ黑人| 亚洲欧美日韩另类电影网站| 一区在线观看完整版| 熟妇人妻不卡中文字幕| 亚洲av.av天堂| 国产色爽女视频免费观看| 久久人人爽人人爽人人片va| 国产精品国产av在线观看| 欧美精品av麻豆av| 成年人免费黄色播放视频| 亚洲高清免费不卡视频| av.在线天堂| 丝袜脚勾引网站| 日产精品乱码卡一卡2卡三| 999精品在线视频| 观看av在线不卡| 内地一区二区视频在线| 日韩中文字幕视频在线看片| 亚洲精品久久成人aⅴ小说| 丰满乱子伦码专区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久韩国三级中文字幕| 精品熟女少妇av免费看| 一级黄片播放器| 制服诱惑二区| 国产乱人偷精品视频| 亚洲国产欧美在线一区| 亚洲精品一二三| 黄色怎么调成土黄色| 色94色欧美一区二区| 欧美人与性动交α欧美软件 | 女人被躁到高潮嗷嗷叫费观| 久久精品久久久久久噜噜老黄| av线在线观看网站| 国产国拍精品亚洲av在线观看| 国产不卡av网站在线观看| 最近中文字幕高清免费大全6| 夫妻性生交免费视频一级片| 亚洲av免费高清在线观看| 亚洲精品自拍成人| 在线亚洲精品国产二区图片欧美| 18禁裸乳无遮挡动漫免费视频| 久久这里有精品视频免费| 久久人人爽av亚洲精品天堂| 性高湖久久久久久久久免费观看| 一本大道久久a久久精品| 午夜福利网站1000一区二区三区| 麻豆乱淫一区二区| 亚洲熟女精品中文字幕| 91成人精品电影| 夫妻午夜视频| 国产欧美亚洲国产| 国内精品宾馆在线| 国产精品偷伦视频观看了| 狠狠精品人妻久久久久久综合| 99热网站在线观看| 黄网站色视频无遮挡免费观看| 日本与韩国留学比较| 内地一区二区视频在线| 亚洲国产欧美在线一区| 免费在线观看完整版高清| 美女福利国产在线| 国产女主播在线喷水免费视频网站| 99热这里只有是精品在线观看| 日韩中字成人| 免费少妇av软件| 91久久精品国产一区二区三区| 看免费成人av毛片| 免费日韩欧美在线观看| 91aial.com中文字幕在线观看| 亚洲国产精品国产精品| 欧美另类一区| 成人免费观看视频高清| a级毛片黄视频| 国产精品一国产av| 少妇精品久久久久久久| 两个人看的免费小视频| 色94色欧美一区二区| 美女国产高潮福利片在线看| 哪个播放器可以免费观看大片| 97在线人人人人妻| 国产亚洲最大av| 看十八女毛片水多多多| 免费日韩欧美在线观看| freevideosex欧美| 亚洲欧美日韩另类电影网站| 99热国产这里只有精品6| 夜夜爽夜夜爽视频| 人人妻人人添人人爽欧美一区卜| 欧美日韩综合久久久久久| 最近的中文字幕免费完整| av免费观看日本| videos熟女内射| 一本色道久久久久久精品综合| 亚洲少妇的诱惑av| 一级毛片 在线播放| 亚洲av男天堂| 欧美 日韩 精品 国产| 在线观看免费日韩欧美大片| 天堂俺去俺来也www色官网| 18在线观看网站| 一级黄片播放器| 香蕉国产在线看| 狠狠精品人妻久久久久久综合| 亚洲av.av天堂| 男人添女人高潮全过程视频| 国产精品一区二区在线不卡| 欧美激情国产日韩精品一区| 啦啦啦中文免费视频观看日本| 国产成人精品在线电影| xxxhd国产人妻xxx| 国产成人午夜福利电影在线观看| 久久ye,这里只有精品| 国产免费又黄又爽又色| 99久久中文字幕三级久久日本| 精品亚洲成国产av| 我的女老师完整版在线观看| 高清在线视频一区二区三区| 亚洲av成人精品一二三区| 精品亚洲乱码少妇综合久久| 亚洲欧洲国产日韩| 99久久人妻综合| 亚洲av福利一区| 国产日韩欧美在线精品| 99热国产这里只有精品6| 久久人人97超碰香蕉20202| 美女国产视频在线观看| 国产乱人偷精品视频| 一区二区av电影网| 丝瓜视频免费看黄片| 五月伊人婷婷丁香| 欧美 日韩 精品 国产| 精品一区二区免费观看| 伊人久久国产一区二区| 日本wwww免费看| 免费黄网站久久成人精品| 国产精品一国产av| 国产精品99久久99久久久不卡 | 欧美日韩国产mv在线观看视频| 免费av不卡在线播放| 春色校园在线视频观看| 免费观看在线日韩| av电影中文网址| 国产熟女欧美一区二区| 热99国产精品久久久久久7| 一级,二级,三级黄色视频| 香蕉精品网在线| 菩萨蛮人人尽说江南好唐韦庄| 日日啪夜夜爽| 国产精品国产三级国产av玫瑰| 免费黄频网站在线观看国产| 一级毛片电影观看| 亚洲国产精品999| 国产精品久久久久久av不卡| 欧美人与善性xxx| 国精品久久久久久国模美| 久久久久国产网址| 青春草亚洲视频在线观看| 蜜桃在线观看..| 制服人妻中文乱码| 久久亚洲国产成人精品v| 国产日韩一区二区三区精品不卡| 捣出白浆h1v1| 最后的刺客免费高清国语| 人成视频在线观看免费观看| 欧美精品亚洲一区二区| 欧美亚洲 丝袜 人妻 在线| 国产69精品久久久久777片| 欧美人与性动交α欧美精品济南到 | 两性夫妻黄色片 | 成人亚洲精品一区在线观看| 18禁在线无遮挡免费观看视频| 欧美xxⅹ黑人| 久久国内精品自在自线图片| 中文天堂在线官网| 亚洲人与动物交配视频| 26uuu在线亚洲综合色| 在线观看人妻少妇| 亚洲av电影在线观看一区二区三区| 波野结衣二区三区在线| 国产黄色视频一区二区在线观看| 中文字幕人妻熟女乱码| 日本欧美视频一区| 色视频在线一区二区三区| 免费女性裸体啪啪无遮挡网站| 少妇人妻精品综合一区二区| 亚洲精品日韩在线中文字幕| 国产一级毛片在线| 亚洲伊人久久精品综合| 青春草视频在线免费观看| 成人18禁高潮啪啪吃奶动态图| 人人妻人人澡人人爽人人夜夜| 久久国内精品自在自线图片| 精品熟女少妇av免费看| 一二三四在线观看免费中文在 | 亚洲,一卡二卡三卡| 九九在线视频观看精品| 免费黄色在线免费观看| 一区在线观看完整版| 丁香六月天网|