亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        加筋矩形薄板的平均聲輻射效率

        2016-11-24 06:36:48任惠娟盛美萍
        振動(dòng)與沖擊 2016年20期
        關(guān)鍵詞:模態(tài)振動(dòng)研究

        任惠娟, 盛美萍

        (1. 咸陽(yáng)師范學(xué)院 物理與電子工程學(xué)院,陜西 咸陽(yáng) 712000; 2.西北工業(yè)大學(xué) 航海學(xué)院,西安 710072)

        ?

        加筋矩形薄板的平均聲輻射效率

        任惠娟1, 盛美萍2

        (1. 咸陽(yáng)師范學(xué)院 物理與電子工程學(xué)院,陜西 咸陽(yáng) 712000; 2.西北工業(yè)大學(xué) 航海學(xué)院,西安 710072)

        采取工程等效的方法研究了四邊簡(jiǎn)支矩形加筋薄板的振動(dòng)與聲輻射,推導(dǎo)了點(diǎn)力激勵(lì)下的平均輻射效率計(jì)算公式,并進(jìn)行了數(shù)值仿真和實(shí)驗(yàn)驗(yàn)證。理論分析與實(shí)驗(yàn)的一致性表明:四邊簡(jiǎn)支情況下矩形加筋薄板與勻質(zhì)矩形薄板和圓薄板的平均輻射效率公式相似,可以采用模態(tài)輻射效率加權(quán)平均的方法等效。在中高頻段,縱橫雙向加筋薄板的平均聲輻射效率最高,單向加筋薄板次之,勻質(zhì)平薄板最低;而在低頻段及臨界頻率以上,三種結(jié)構(gòu)平均聲輻射效率相差不大。

        加筋薄板;正交各向異性板;平均聲輻射效率

        在船舶、飛機(jī)等大型結(jié)構(gòu)中,為了增加結(jié)構(gòu)的剛度,通常在勻質(zhì)板殼結(jié)構(gòu)上附加加強(qiáng)筋。加強(qiáng)筋的存在改變了勻質(zhì)結(jié)構(gòu)的振動(dòng)特性,也使得其聲振分析更為復(fù)雜。

        對(duì)加筋薄板的聲輻射研究,早期主要是以無(wú)限大具有相同加強(qiáng)筋的板為研究對(duì)象[1-2]。其中,MACE[1]借助波數(shù)變換和穩(wěn)相法對(duì)于無(wú)限流場(chǎng)作用下有平行加強(qiáng)筋和垂直加強(qiáng)筋的板結(jié)構(gòu)在點(diǎn)激勵(lì)下的遠(yuǎn)場(chǎng)聲壓進(jìn)行了研究。文獻(xiàn)[3-4]研究加筋板的聲輻射問(wèn)題時(shí)僅考慮了加強(qiáng)筋的振動(dòng)效應(yīng)。MOYNE等[5]的研究則認(rèn)為加強(qiáng)筋對(duì)加筋板的聲輻射有振動(dòng)效應(yīng)、聲源效應(yīng)及聲障效應(yīng)三種影響,并指出其聲源效應(yīng)和聲障效應(yīng)可分別用E準(zhǔn)則和λ/6準(zhǔn)則來(lái)評(píng)估。文獻(xiàn)[6]則以四邊簡(jiǎn)支單向?qū)ΨQ加筋矩形板為研究對(duì)象,將薄板沿加強(qiáng)筋劃分成兩個(gè)簡(jiǎn)單單元,運(yùn)用反力法將加強(qiáng)筋的作用等效在板上,研究了加筋板的聲輻射性能。文獻(xiàn)[7]采用有限元/邊界元方法,對(duì)加筋板結(jié)構(gòu)的聲輻射進(jìn)行了計(jì)算分析,研究了加筋板的結(jié)構(gòu)參數(shù)、邊界條件、流體介質(zhì)等因素對(duì)于板的輻射聲功率的影響。文獻(xiàn)[8]在波數(shù)域上求解了流體介質(zhì)中具有相同加強(qiáng)筋平板在集中力作用下的聲輻射。文獻(xiàn)[9]采用有限元/邊界元法研究了手性覆蓋層對(duì)不均勻結(jié)構(gòu)加筋梁聲輻射的影響, 結(jié)果表明手性覆蓋層不僅能夠有效抑制不均勻結(jié)構(gòu)及其傳遞到覆蓋層輻射面的振動(dòng), 而且對(duì)不均勻結(jié)構(gòu)的輻射聲功率有很好的抑制作用。以上研究從不同方面揭示了加筋結(jié)構(gòu)的聲振特性,并為進(jìn)一步深入研究加筋板的聲輻射效率提供了幫助。

        振動(dòng)結(jié)構(gòu)的平均聲輻射效率是其聲輻射性能的一個(gè)重要表征參量,反映了振動(dòng)結(jié)構(gòu)向周圍介質(zhì)中輻射聲能量的能力,是工程技術(shù)人員在進(jìn)行聲學(xué)結(jié)構(gòu)設(shè)計(jì)時(shí)需要考慮的重要參數(shù)。XIE等[10]曾就勻質(zhì)矩形薄板和條形結(jié)構(gòu)的平均輻射效率進(jìn)行了詳細(xì)研究,給出了其平均聲輻射效率的計(jì)算公式。本文作者曾在該文獻(xiàn)的研究基礎(chǔ)上,對(duì)勻質(zhì)圓形薄板結(jié)構(gòu)的平均輻射效率進(jìn)行過(guò)詳細(xì)研究,提出了勻質(zhì)薄圓板平均輻射效率的計(jì)算公式[11]。本文將在以上研究的基礎(chǔ)上結(jié)合加筋薄板的正交各向異性等效理論,進(jìn)一步研究矩形加筋薄板的平均輻射效率計(jì)算方法。

        1 加筋薄板與正交各向異性板的等效條件

        周期加筋的矩形薄板通??梢缘刃檎桓飨虍愋园?。HECKEL[12]指出:當(dāng)筋的間距小于1/4彎曲波長(zhǎng)時(shí),周期加筋薄板便可等效為正交各向異性平板。由于平板上彎曲波的波速為

        (1)

        (2)

        或者其振動(dòng)頻率滿足

        (3)

        時(shí),加筋薄板便可以等效為正交各向異性板。

        SUNDARA K T研究了滿足Heckel等效條件的加筋薄板和正交各向異性板之間的等效關(guān)系,并給出了下列等效公式[13]。

        (4)

        (5)

        Dxy=

        (6)

        D1=σD

        (7)

        (8)

        2 正交各向異性板的受迫振動(dòng)響應(yīng)

        為簡(jiǎn)單起見(jiàn),分析時(shí)忽略板的阻尼,則正交各向異性板的受迫振動(dòng)方程為

        (9)

        對(duì)于四邊簡(jiǎn)支矩形各向異性板,設(shè)其撓度解的形式為

        (10)

        當(dāng)載荷為幅值為F0并作用于板面(x0,y0)上的點(diǎn)簡(jiǎn)諧力qt=F0δ(x-x0)δ(y-y0)ejωt時(shí),將該載荷展開(kāi)成振型函數(shù)的級(jí)數(shù)形式為

        (11)

        將式(10)和式(11)代入式(9)并比較系數(shù)得

        (12)

        式中:ωmn為正交各向異性板的自然頻率,表示為

        (13)

        于是正交各向異性板的振速為

        u(x,y,t)=jωw(x,y,t)

        (14)

        其振速幅值表示為

        Umn=jωL(ω)φmn(x0,y0)φmn(x,y)

        (15)

        3 加筋薄板的平均聲輻射效率

        為了研究加筋板的平均聲輻射效率問(wèn)題,將符合Heckel等效條件的加筋薄板等效為正交各向異性板,令其邊界條件為四邊簡(jiǎn)支,并將之嵌于無(wú)限大剛性障板中。

        3.1 加筋薄板的模態(tài)聲功率及其平均值

        建立圖1的坐標(biāo)系,將該正交各向異性板表面分割成無(wú)限多個(gè)小面元,面元ds′相對(duì)于ds的位置用極坐標(biāo)(h,θ)表示,其中h為兩面元之間的距離,θ為h與x軸反方向之間的夾角。當(dāng)正交各向異性板受到點(diǎn)簡(jiǎn)諧力qt=F0δ(x-x0)δ(y-y0)ejωt激勵(lì)激勵(lì)時(shí),ds′面元的振動(dòng)在ds處產(chǎn)生的聲壓幅值為

        (16)

        式中:k為聲波的波數(shù),ρa(bǔ)為空氣密度,ca為聲波在空氣中的傳播速度,ω為振動(dòng)的角頻率。對(duì)上式積分得到所有面元的振動(dòng)在面元ds處產(chǎn)生的聲壓為

        Pnm(x,y)=∫sdPmn(x,y)

        (17)

        面元ds的聲強(qiáng)為

        (18)

        Wnm=∫sInm(x,y)ds

        (19)

        對(duì)激勵(lì)點(diǎn)進(jìn)行面平均,得到模態(tài)輻射聲功率的平均值為

        (20)

        式中,S為板的面積。

        圖1 正交各向異性矩形板表面積分示意圖Fig.1 Integral sketch in the surface of othotropic plate

        3.2 加筋薄板的模態(tài)輻射效率

        汽車和發(fā)動(dòng)機(jī)在整個(gè)壽命期間如何減少故障、延長(zhǎng)使用壽命、提高在有限時(shí)間內(nèi)的工作質(zhì)量始終是個(gè)難題。隨著時(shí)間的發(fā)展,狀態(tài)監(jiān)控技術(shù)應(yīng)運(yùn)而生,它成為解決這個(gè)難題的重要手段和方法。狀態(tài)監(jiān)控不是判定結(jié)果(故障診斷),而是養(yǎng)護(hù)和維修“注意”和“預(yù)測(cè)”。綜合后的SAE狀態(tài)監(jiān)測(cè)概念列于表2,信息流程如圖1所示。

        等效正交各向異性矩形板的模態(tài)均方振速為

        (21)

        對(duì)上式中的激勵(lì)點(diǎn)作空間平均,得到模態(tài)均方振速的平均值為

        (22)

        等效正交各向異性矩形板的模態(tài)輻射效率為

        (23)

        將式(20)和式(22)代入式(23),則等效正交各向異性矩形板的模態(tài)輻射效率為

        φmn(x,y)·φmn(x′,y′)dsds′

        (24)

        將上式離散化為

        (25)

        3.3 加筋薄板的平均輻射效率

        等效正交各向異性矩形板振動(dòng)時(shí)的平均輻射效率為

        (26)

        (27)

        (28)

        將式(27)和式(28)代入式(26)得到嵌于無(wú)限大剛性障板中的加筋矩形薄板在點(diǎn)激勵(lì)下的平均輻射效率為

        (29)

        式中,lmn定義為

        (30)

        由式(29)可以看出,四邊簡(jiǎn)支加筋薄板的平均輻射效率是其各階模態(tài)輻射效率的加權(quán)平均值,其權(quán)系數(shù)由式(30)決定。將 式(29) 及 式(30) 與勻質(zhì)薄板的平均輻射效率公式[10]對(duì)比發(fā)現(xiàn),二者平均聲輻射效率表達(dá)式及權(quán)系數(shù)的表達(dá)式均非常類似,即平均輻射效率都等于其模態(tài)輻射效率的加權(quán)平均值,而不同之處在于二者權(quán)系數(shù)表達(dá)式中對(duì)自然頻率及板的面密度定義的不同。

        4 數(shù)值分析與討論

        對(duì)于一縱橫雙向加筋的方薄板,a=1 m,b=1 m,hP=0.003 m,雙向筋的參數(shù)相同:筋間距dr=0.02 m,筋條寬度br=0.005 m,筋高h(yuǎn)r=0.02 m,筋條總數(shù)nr=50;筋與板的材質(zhì)均為鋼(楊氏模量E=21.6×1010N/m2,體密度ρV=7.8×103kg/m3,泊松比0.28)。根據(jù)加筋薄板的正交各向異性等效條件,當(dāng)頻率滿足f<4 623 Hz時(shí),該加筋薄板可以等效為正交各向異性板。圖2為該加筋薄板的輻射效率隨頻率的變化關(guān)系曲線。

        由圖2可以看出:低頻段,隨著頻率的增大,平均輻射效率也隨之增大,并且平均輻射效率等于(1,1)階的模態(tài)輻射效率;中頻段,隨著頻率的增大,平均輻射效率則圍繞某一恒值在小范圍內(nèi)波動(dòng)變化;高頻段,隨著頻率的增大平均輻射效率也逐漸增大;臨界頻率之后,平均輻射效率將趨近于1。以上特點(diǎn)與勻質(zhì)矩形薄板[10]及勻質(zhì)圓薄板[11]的平均輻射效率特點(diǎn)一致。

        圖2 雙向加筋薄板的輻射效率……模態(tài)輻射效率 ━━平均輻射效率 ——(1,1)階模態(tài)輻射效率Fig 2 Radiation efficiency of a crisscross thin stiffened plate

        圖3給出了雙向加筋薄板、單向加筋薄板(筋條僅沿x方向,筋及板的所有參數(shù)同雙向加筋方板)及勻質(zhì)薄板的平均輻射效率對(duì)比。由圖3可以看出:在中高頻段,加筋薄板的平均聲輻射效率遠(yuǎn)遠(yuǎn)高于勻質(zhì)板,并且縱橫雙向加筋薄板的平均聲輻射效率高于單向加筋薄板,這是由于加筋薄板的剛度較無(wú)筋板大,而縱橫雙向加筋薄板的剛度較單向加筋薄板大的緣故;而在低頻段及臨界頻率以上,三種結(jié)構(gòu)平均聲輻射效率相差不大。

        圖3 雙向加筋、單向加筋及勻質(zhì)薄板平均輻射效率對(duì)比Fig.3 Average radiation efficiency comparison of the crisscross stiffened plate, the unidirectional stiffened plate and the uniform plate

        5 實(shí)驗(yàn)測(cè)試

        為了驗(yàn)證點(diǎn)簡(jiǎn)諧力激勵(lì)下加筋矩形薄板平均輻射效率的理論分析結(jié)果,作者在西北工業(yè)大學(xué)消聲室對(duì)單向加筋薄板的平均聲輻射效率進(jìn)行了實(shí)驗(yàn)驗(yàn)證。

        加筋板的輻射指數(shù)10lgσ根據(jù)下式進(jìn)行測(cè)試:

        (31)

        實(shí)驗(yàn)中,聲功率級(jí)是通過(guò)測(cè)試加筋薄板附近半球表面的聲壓獲得的,其中聲壓的測(cè)試依照《GB/T 6882—1986》附錄B及《GB/T 6882—2008》標(biāo)準(zhǔn)執(zhí)行,平均振速級(jí)則通過(guò)測(cè)試其表面的加速度獲得。

        實(shí)驗(yàn)時(shí),加筋薄板的參數(shù)為:長(zhǎng)為0.5 m、寬為0.5 m、厚為0.005 m、筋高為0.02 m、筋寬為0.005 m、筋條間隔為0.05 m、筋條數(shù)為10,筋和板的材質(zhì)均為鋼。當(dāng)振動(dòng)頻率滿足f<1 232 Hz時(shí),將其等效為正交各向異性板。由于板面距離地面較近,地面在此充當(dāng)了矩形板周圍的無(wú)限大剛性障板,實(shí)驗(yàn)條件和理論分析條件基本符合。在板面、地面形成的空腔中填充吸聲海綿以減少空腔的聲振耦合,激振器懸掛在矩形板上方,激振點(diǎn)偏離板中心。測(cè)試安裝現(xiàn)場(chǎng)如圖4所示。

        圖4 簡(jiǎn)支邊界加筋方薄板的安裝測(cè)試現(xiàn)場(chǎng)Fig.4 The test scene of a stiffened plate with simply supported boundary condiation

        為了實(shí)現(xiàn)加筋薄板的簡(jiǎn)支邊界條件,本文設(shè)計(jì)并制造了簡(jiǎn)支邊界夾具。夾具的材質(zhì)為鋼,與板接觸的地方做了削尖處理。其剖面示意圖如圖5所示。

        圖5 簡(jiǎn)支夾具剖面示意圖Fig.5 A section through the fixture with simply supported boundary condiation

        圖6 簡(jiǎn)支邊界加筋方薄板輻射指數(shù)測(cè)量值與理論值對(duì)比Fig.6 Radiation index comparison of a stiffened plate between the test value with the theroy value

        圖6給出了加筋薄板輻射指數(shù)(10lgσ)測(cè)量值與理論值對(duì)比結(jié)果,由二者的對(duì)比情況看,實(shí)驗(yàn)值與理論值一致性較好。其中實(shí)驗(yàn)值略大是由于實(shí)驗(yàn)過(guò)程中空腔內(nèi)聲振耦合不能完全消除所導(dǎo)致的,但總體來(lái)講理論值與實(shí)驗(yàn)值一致性良好。

        6 結(jié) 論

        本文采用工程等效方法研究了矩形加筋薄板嵌于無(wú)限大障板中的振動(dòng)與聲輻射,得到了模態(tài)輻射效率及平均輻射效率的計(jì)算公式,并通過(guò)實(shí)驗(yàn)驗(yàn)證了理論分析的正確性,理論與實(shí)驗(yàn)的一致性表明:

        (1) 四邊簡(jiǎn)支情況下矩形加筋薄板與勻質(zhì)矩形薄板和圓薄板的平均輻射效率公式相似,可以采用模態(tài)輻射效率加權(quán)平均的方法等效。

        (2) 與勻質(zhì)矩形薄板及勻質(zhì)圓薄板的平均聲輻射效率特點(diǎn)一致,臨界頻率之后,矩形加筋薄板平均輻射效率將趨近于1。低頻段和高頻段,矩形加筋薄板平均輻射效率隨著頻率的增大而增大;中頻段,平均輻射效率隨著頻率的變化在某一恒值內(nèi)小范圍波動(dòng)。

        (3) 在中高頻段,縱橫雙向加筋薄板的平均聲輻射效率最高,單向加筋薄板次之,勻質(zhì)平薄板最低;而在低頻段及臨界頻率以上,三種結(jié)構(gòu)平均聲輻射效率相差不大。

        [1] MACE B R. Periodically stiffened fluid-load platesⅡ: response to line and point forces [J]. J Sound and Vibration,1980,73(4):487-504.

        [2] MEAD D J. Plates with regular stiffening in acoustic media: vibration and radiation[J]. J Acoust Soc Am, 1990,88:391-401.

        [3] FAHY F. Sound and structural vibration: radiation, transmission and response[M]. London: Academic Press,1985.

        [4] MAIDANIK G, DICKEY J. Acoustic behavior of ribbed planels[J]. J Sound and Vibration, 1988: 293-307.

        [5] MOYNE S L, TEBEC J L, KRAEMER J C. Source effect of ribs in siund radiation of stiffened plates: experimental and calculation investigation[J]. Acustica, 2000,86: 457-464.

        [6] 姚本炎,黃其柏. 加筋薄板結(jié)構(gòu)振動(dòng)與聲輻射特性研究[J].華中科技大學(xué)學(xué)報(bào),2001,29(2):90-92.

        YAO Benyan, HUANG Qibai. Research on the vibration and sound radiation characteristics of stiffened thin plate[J]. Journal of Huazhong University of Science and Technology, 2001,29(2):90-92.

        [7] 黎勝,趙德有. 流體加載下加筋板結(jié)構(gòu)的聲輻射特性研究[J].應(yīng)用聲學(xué),2000,19(6):28-32.

        LI Sheng, ZHAO Deyou. Study on the acoustic radiation characteristics of a fluid-loaded stiffened plate[J]. Applied Acoustic,2000,19(6):28-32.

        [8] 吳文偉,沈順根.具有等間距相同加強(qiáng)筋板的聲輻射[J]. 船舶力學(xué),1998, 2(5): 63-72.

        WU Wenwei, SHEN Shungen. The sound radiation of a plate with equidistant reinforced ribs[J]. Journal of Ship Mechanics,1998, 2(5): 63-72.

        [9] 朱大巍,黃修長(zhǎng),華宏星,等.敷設(shè)手性覆蓋層加筋梁低頻振動(dòng)和聲輻射特性[J].振動(dòng)與沖擊,2014,33(11):178-183.

        ZHU Dawei, HUANG Xiuchang, HUA Hongxing, et al.Vibration and acoustic radiation characteristics of a stiffened beam with a chiral covering layer[J].Journal of Vibration and Shock, 2014,33(11):178-183.

        [10] XIE G, THOMPSON D J, JONES C J C. The radiation efficiency of baffled plates and strips[J]. Journal of Sound and Vibration, 2005, 280:181-209.

        [11] 任惠娟,盛美萍. 彎曲振動(dòng)薄圓板的輻射效率[J]. 振動(dòng)與沖擊, 2012, 31(20):121-125.

        REN Huijuan, SHENG Meiping. Radiation efficiency of a thin circular plate vibrating flexurally[J]. Journal of Vibration and Shock, 2012,31(20):121-125.

        [12] HECKL M. Wave propagation on beam-plate systems[J].Journal of the Acoustical Society of America, 1961, 33:640-651.

        [13] GRAFF K F. Wave motion in elastic solids[M]. New York: Dover Publications INC,1975.

        The average radiation efficiency of a rectangular stiffened thin plate

        REN Huijuan1, SHENG Meiping2

        (1. Department of Physics and Electronic Engineering, Xianyang Normal University, Xianyang 712000, China;2. Department of Marine, Northwestern Polytechnical University, Xi’an 710072, China)

        The vibration and sound radiation of a rectangular stiffened thin plate were investigated by using the method of the engineering equivalence. The average radiation efficiency formula of the rectangular stiffened plate driven by a simple harmonic point force was derived. The numerical simulation and the experiment were implemented. The consistency of the theory analysis and the experiment shows that the average radiation efficiency formula of the rectangular stiffened thin plate is similar to that of the uniform rectangular thin plate and circular thin plate. They are all the weighted mean of the modal radiation efficiency. In the middle and high frequency band, the average radiation efficiency of the crisscross stiffened thin plate is maximal. The unidirectional stiffened thin plate takes the second place, and the uniform thin plate is the minimum. While in the low frequency band and above the critical frequency, there is no significant difference between them.

        rectangular stiffened thin plate; orthotropic plate; average radiation efficiency

        陜西省自然科學(xué)基金 (2011JM1017);陜西省教育廳2013年科學(xué)研究計(jì)劃項(xiàng)目(2013JK0618)

        2015-12-24 修改稿收到日期:2016-03-17

        任惠娟 女,博士,副教授,1972年生

        O42

        A

        10.13465/j.cnki.jvs.2016.20.027

        猜你喜歡
        模態(tài)振動(dòng)研究
        振動(dòng)的思考
        FMS與YBT相關(guān)性的實(shí)證研究
        遼代千人邑研究述論
        振動(dòng)與頻率
        視錯(cuò)覺(jué)在平面設(shè)計(jì)中的應(yīng)用與研究
        科技傳播(2019年22期)2020-01-14 03:06:54
        EMA伺服控制系統(tǒng)研究
        中立型Emden-Fowler微分方程的振動(dòng)性
        國(guó)內(nèi)多模態(tài)教學(xué)研究回顧與展望
        基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識(shí)別
        UF6振動(dòng)激發(fā)態(tài)分子的振動(dòng)-振動(dòng)馳豫
        无码AV高潮喷水无码专区线| 久久久久亚洲精品无码系列| 大地资源网高清在线播放| 欧美性xxxx狂欢老少配| 久久久午夜毛片免费| 日本一区二区在线播放| 人禽杂交18禁网站免费| 午夜精品久久久久久| 在线观看网址你懂的| 国产一级r片内射视频播放| 美女下蹲露大唇无遮挡| 成人免费777777被爆出| 无码日韩AⅤ一区二区三区| 国产黑色丝袜在线观看网站91 | 欧美另类高清zo欧美| 亚洲熟妇色xxxxx欧美老妇| 少妇勾引视频网站在线观看| 亚洲精品国产av成人精品| 亚洲亚洲人成综合网络| 国产免费一级在线观看| 亚洲中文字幕在线精品2021| 国产美女高潮流白浆免费视频| 黄瓜视频在线观看| 亚洲AV无码永久在线观看| 日韩精品视频中文字幕播放| 亚洲a无码综合a国产av中文| 国产香蕉尹人在线观看视频| 久久老子午夜精品无码| 亚洲国产国语对白在线观看 | 日韩欧美在线播放视频| 国成成人av一区二区三区| 久久精品中文字幕无码绿巨人| 女女女女bbbbbb毛片在线| 欧美人与动牲交片免费| 国产三级精品三级在线专区2| 亚洲热妇无码av在线播放| 亚洲美女啪啪| 热综合一本伊人久久精品| 成 人片 黄 色 大 片| 久久精品成人欧美大片| baoyu网址国产最新|