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

    水中振動(dòng)矩形板的附加質(zhì)量研究

    2016-12-12 05:34:43任惠娟薛小慶
    關(guān)鍵詞:面密度波數(shù)高階

    任惠娟, 姚 展, 薛小慶, 劉 婷, 雷 燁

    ?

    水中振動(dòng)矩形板的附加質(zhì)量研究

    任惠娟1, 姚 展1, 薛小慶1, 劉 婷1, 雷 燁2

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

    從水下振動(dòng)矩形板的振動(dòng)響應(yīng)公式出發(fā),根據(jù)附加質(zhì)量的定義,推導(dǎo)了水下矩形板的附加質(zhì)量表達(dá)式,并進(jìn)行了數(shù)值仿真.研究表明:(1)模態(tài)階數(shù)越低,矩形板附加面密度越大.(2)當(dāng)波數(shù)比kw/kmnw<10-1時(shí),矩形板的附加面密度與波數(shù)比無(wú)關(guān)而幾乎保持恒值;當(dāng)kw/kmnw≈1時(shí),除過(guò)(1,1)階模態(tài),各階模態(tài)附加面密度均出現(xiàn)了峰值;之后隨著波數(shù)比的增大,模態(tài)附加面密度迅速衰減.(3)水下矩形板的附加質(zhì)量不可忽略,研究水下結(jié)構(gòu)的振動(dòng)與聲輻射必須考慮附加質(zhì)量的影響.

    矩形板; 附加質(zhì)量; 附加阻尼

    0 引言

    當(dāng)結(jié)構(gòu)在重質(zhì)流體中振動(dòng)時(shí), 由于流體與結(jié)構(gòu)的相互耦合作用[1-6]而使其振動(dòng)和聲輻射變得復(fù)雜.流固耦合問(wèn)題的弱耦合研究法[7]認(rèn)為流體對(duì)結(jié)構(gòu)的影響是以附加質(zhì)量和附加阻尼的形式體現(xiàn)出來(lái)的, 因此圍繞各種結(jié)構(gòu)在流體中附加質(zhì)量的研究很受人們的關(guān)注.王基盛等[7]采用弱耦合法分析了流體與結(jié)構(gòu)的耦合作用,探討了流體環(huán)境中影響結(jié)構(gòu)附加質(zhì)量的因素,并給出了流體動(dòng)能的計(jì)算方法.錢(qián)勤等[8]借助Green公式,將結(jié)構(gòu)附連水質(zhì)量的計(jì)算化為一個(gè)邊界積分,計(jì)算了水中圓柱體和球體的附連水質(zhì)量.于丹竹等[9]通過(guò)建立帶有附加質(zhì)量主結(jié)構(gòu)的有限元模型,將其與聲流體介質(zhì)的瑞利積分相互耦合以獲得耦合系統(tǒng)的響應(yīng).馬燁等[10]采用商用計(jì)算軟件,運(yùn)用Fluent非常定運(yùn)動(dòng)和動(dòng)網(wǎng)格技術(shù)數(shù)值計(jì)算得到了物體六自由度運(yùn)動(dòng)的附加質(zhì)量.劉瑩等[11]在勢(shì)流理論的基礎(chǔ)上,求解了結(jié)構(gòu)剛彈耦合的附加質(zhì)量.李明等[12]通過(guò)建立水下航行體的梁模型,研究了其運(yùn)動(dòng)過(guò)程中附加質(zhì)量的變化情況.以上研究針對(duì)不同水下結(jié)構(gòu)采用不同方法研究了其附加質(zhì)量問(wèn)題,也從不同方面剖析了水下結(jié)構(gòu)附加質(zhì)量的產(chǎn)生機(jī)理,為求解結(jié)構(gòu)在流體中的附加質(zhì)量提供了很好的幫助.板結(jié)構(gòu)是構(gòu)成各種復(fù)雜結(jié)構(gòu)的基本結(jié)構(gòu)之一,本文則以水中點(diǎn)簡(jiǎn)諧力激勵(lì)下的振動(dòng)矩形板為研究對(duì)象,根據(jù)弱耦合研究理論, 對(duì)其振動(dòng)響應(yīng)公式進(jìn)行修正,結(jié)合振動(dòng)結(jié)構(gòu)附加質(zhì)量的定義,推導(dǎo)其附加質(zhì)量的計(jì)算方法.本研究為進(jìn)一步研究水介質(zhì)對(duì)水下振動(dòng)結(jié)構(gòu)聲輻射的影響奠定基礎(chǔ).

    1 水中振動(dòng)矩形板的振動(dòng)響應(yīng)

    當(dāng)結(jié)構(gòu)在水中振動(dòng)時(shí),結(jié)構(gòu)會(huì)向水中輻射聲場(chǎng),該聲場(chǎng)將反作用于結(jié)構(gòu)本身,從而對(duì)結(jié)構(gòu)自身的振動(dòng)和聲輻射產(chǎn)生影響.聲場(chǎng)對(duì)振動(dòng)結(jié)構(gòu)的影響體現(xiàn)在兩個(gè)方面[7,13]:一方面由于結(jié)構(gòu)的振動(dòng)帶動(dòng)了周?chē)黧w的運(yùn)動(dòng),這相當(dāng)于振動(dòng)結(jié)構(gòu)自身的質(zhì)量產(chǎn)生了增量,我們稱(chēng)之為附加質(zhì)量; 另一方面,由于結(jié)構(gòu)在振動(dòng)的同時(shí)向周?chē)慕橘|(zhì)中輻射出聲能量從而損耗了機(jī)械能,這部分能量的損耗可以用附加阻尼(又稱(chēng)為輻射阻尼)來(lái)進(jìn)行描述.

    設(shè)長(zhǎng)為a、寬為b的矩形板置于水中,板的邊界條件為四邊簡(jiǎn)支并且被嵌于無(wú)限大障板中(如圖1所示).當(dāng)矩形板受到作用于點(diǎn)(x0,y0)處的簡(jiǎn)諧力Fejωt激勵(lì)而發(fā)生振動(dòng)時(shí),考慮到振動(dòng)過(guò)程中矩形板附加質(zhì)量和附加阻尼的產(chǎn)生,將空氣介質(zhì)中點(diǎn)(x′,y′)的振速幅值表達(dá)式[14]中板的密度和阻尼分別進(jìn)行修正,得到水下矩形板在點(diǎn)簡(jiǎn)諧激勵(lì)下的振速幅值表達(dá)式為:

    (1)

    D為板的彎曲剛度.

    圖1 矩形板積分示意圖

    2 振動(dòng)矩形板的附加質(zhì)量

    當(dāng)矩形板受到作用于點(diǎn)(x0,y0) 處的簡(jiǎn)諧力Fejωt激勵(lì)而發(fā)生振動(dòng)時(shí),面圓ds′的振動(dòng)在ds處產(chǎn)生的聲壓幅值為

    (2)

    式(2)中:kw為水中聲波波數(shù),ρw為水的密度,cw為水中的聲速.將(1)式代入(2)式并對(duì)整個(gè)輻射面進(jìn)行積分,得到所有面元的振動(dòng)在面元ds處產(chǎn)生的聲壓為

    (3)

    于是面圓ds受到聲場(chǎng)的反作用力為

    Fmnw(x,y)=Pmnw(x,y)·ds

    (4)

    面元ds的輻射阻抗為

    (5)

    面元ds的附加質(zhì)量[13]為

    (6)

    由于板的振動(dòng),面圓的附加質(zhì)量也隨之一起振動(dòng),于是面元ds的附加質(zhì)量對(duì)應(yīng)的附加動(dòng)能為

    (7)

    將(6)式代入(7)式并對(duì)(7)式積分,得到整個(gè)板面的附加動(dòng)能為

    (8)

    從另一方面看,若令整個(gè)板面的(m,n)模態(tài)的附加質(zhì)量為Δmmn,考慮到整個(gè)板面的振速均方值為

    (9)

    則整個(gè)板面的附加質(zhì)量引起的整個(gè)板面的附加動(dòng)能又可寫(xiě)為

    (10)

    于是由T′=T得

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

    (11)

    (11)式即為水中振動(dòng)矩形板(m,n)模態(tài)的附加質(zhì)量表達(dá)式.

    3 數(shù)值仿真

    圖2給出了邊長(zhǎng)為1 m的四邊簡(jiǎn)支方板前9階模態(tài)的附加面密度隨著波數(shù)比的變化情況.圖中橫坐標(biāo)為輻射聲波的波數(shù)kw和矩形板(m,n)階模態(tài)彎曲波波數(shù)kmnw的比值,即波數(shù)比.

    圖2表明,在雙對(duì)數(shù)坐標(biāo)系下,水中簡(jiǎn)支邊界方板的附加面密度隨波數(shù)比的變化關(guān)系呈現(xiàn)出以下特點(diǎn):

    (1)模態(tài)階數(shù)越低,方板附加面密度越大.

    (2)當(dāng)kw?kmnw時(shí),方板各階模態(tài)的附加面密度曲線(xiàn)近似為與橫軸平行的直線(xiàn),即方板的模態(tài)附加面密度幾乎與波數(shù)比無(wú)關(guān).

    (3)當(dāng)kw

    (4)當(dāng)kw≈kmnw時(shí),除(1,1)模態(tài)外,各模態(tài)的附加面密度均出現(xiàn)峰值,且模態(tài)階數(shù)越高,峰值越明顯.

    (5)當(dāng)kw?kmnw時(shí),所有模態(tài)的附加面密度均迅速減小.

    圖2 簡(jiǎn)支邊界矩形板前9階模態(tài)的附加面密度

    圖3給出了邊長(zhǎng)為1 m的方板在簡(jiǎn)支邊界條件下單面臨水及在空氣中的附加面密度對(duì)比情況.

    由圖3可以看出,當(dāng)方板在水中振動(dòng)時(shí),其附加面密度遠(yuǎn)遠(yuǎn)大于其在空氣中的情形.另外,對(duì)于邊長(zhǎng)為1 m、厚度為6 mm的鋼質(zhì)方板在空氣中振動(dòng)時(shí),計(jì)算表明,其附加面密度遠(yuǎn)在1 kg/m2數(shù)量級(jí)以下,遠(yuǎn)遠(yuǎn)小于方板自身的面密度(約為47 kg/m2).這說(shuō)明在空氣介質(zhì)中計(jì)算板的振動(dòng)及聲輻射時(shí)由于振動(dòng)板的附加質(zhì)量遠(yuǎn)遠(yuǎn)小于自身質(zhì)量,因而可以完全不考慮其附加質(zhì)量的影響.而當(dāng)該板在水中振動(dòng)時(shí),其附加面密度在102kg/ m2數(shù)量級(jí)左右,遠(yuǎn)大于或接近于其自身的面密度,所以當(dāng)板在水下振動(dòng)時(shí),必須考慮其附加質(zhì)量的影響.

    (a) (1,1)模態(tài) (b) (6,3)模態(tài)

    (c) (4,4)模態(tài) (d) (3,2)模態(tài)圖3 簡(jiǎn)支邊界方板在空氣及水中的附加面密度對(duì)比

    4 本文附加質(zhì)量與波數(shù)變換法附加質(zhì)量計(jì)算結(jié)果對(duì)比

    文獻(xiàn)[15]用波數(shù)變換法推導(dǎo)出了臨界頻率以下頻段簡(jiǎn)支邊界矩形板的附加質(zhì)量,其表達(dá)式為:

    (12)

    文獻(xiàn)[15]同時(shí)指出,這一表達(dá)式適用于計(jì)算較高階模態(tài)矩形板的附加質(zhì)量,當(dāng)其用于計(jì)算較低階模態(tài)的附加質(zhì)量時(shí)將會(huì)有較大的誤差.這里以邊長(zhǎng)為0.5 m的方板為例,將本文方法的低階模態(tài)計(jì)算結(jié)果(如圖4所示)及高階模態(tài)計(jì)算結(jié)果(如圖5所示)分別與上式計(jì)算結(jié)果進(jìn)行對(duì)比.

    對(duì)于圖4的各低階模態(tài),本文的計(jì)算結(jié)果均高于(12)式的計(jì)算結(jié)果,尤其是對(duì)于(1,1)階模態(tài),本文的結(jié)果遠(yuǎn)遠(yuǎn)高于(12)式的結(jié)果.但是圖4也表明隨著模態(tài)階數(shù)的逐步升高,本文的計(jì)算結(jié)果越來(lái)越接近于(12)式的計(jì)算結(jié)果.對(duì)于圖5的各高階模態(tài),本文的計(jì)算結(jié)果與(12)式的結(jié)果一致性良好.

    (a) (1,1)模態(tài) (b) (1,2)模態(tài)

    (c) (1,3)模態(tài) (d) (2,2)模態(tài)圖4 低階模態(tài)附加質(zhì)量對(duì)比

    (a) (4,5)模態(tài) (b) (5,6)模態(tài)

    (c) (6,7)模態(tài) (d) (7,7)模態(tài)圖5 高階模態(tài)附加質(zhì)量對(duì)比

    以上對(duì)比結(jié)果與文獻(xiàn)[15]中指出的(12)式適用范圍相吻合.由于(12)式在推導(dǎo)的過(guò)程中做了一定的近似處理,這種近似對(duì)于高階模態(tài)是適合的,但對(duì)于低階模態(tài)則不適用[15],因此(12)式僅適用于計(jì)算高階模態(tài)的附加質(zhì)量.本文給出的矩形板附加質(zhì)量計(jì)算公式(11)式是根據(jù)弱耦合理論嚴(yán)格推導(dǎo)得到的.(11)式和(12)式是在不同的物理模型下推導(dǎo)的結(jié)果,其高階模態(tài)計(jì)算結(jié)果良好的一致性則表明了本文物理模型的有效性及推導(dǎo)結(jié)果正確性.另外,從仿真結(jié)果上看,對(duì)于低階模態(tài),隨著模態(tài)階數(shù)的逐步升高,兩種計(jì)算結(jié)果的差別有逐漸縮小的趨勢(shì),到了高階模態(tài)時(shí)則達(dá)到了完全一致的地步,而這恰恰是(12)式推導(dǎo)過(guò)程中所做的近似處理更適用于高階模態(tài)的具體反映.因此,不僅由該模型所推得的高階模態(tài)附加質(zhì)量計(jì)算結(jié)果是正確的,而且其低階模態(tài)附加質(zhì)量計(jì)算結(jié)果也是可信的.

    本文的附加質(zhì)量計(jì)算方法由于既可用于計(jì)算高階模態(tài),又可用于計(jì)算低階模態(tài),因此較(12)式適用范圍更廣.

    5 結(jié)論

    本文從修正后的水下矩形板振動(dòng)響應(yīng)公式出發(fā),推導(dǎo)了其附加質(zhì)量表達(dá)式,并根據(jù)該表達(dá)式進(jìn)行了數(shù)值仿真,結(jié)果表明:(1)模態(tài)階數(shù)越低,矩形板附加面密度越大.(2)當(dāng)波數(shù)比kw/kmnw<10-1時(shí),矩形板的附加面密度幾乎保持恒值而與波數(shù)比無(wú)關(guān);隨著波數(shù)比的增大,附加面密度逐漸波動(dòng)變化;當(dāng)kw/kmnw≈1時(shí),除過(guò)(1,1)階模態(tài),其他模態(tài)的附加面密度均出現(xiàn)了峰值;之后隨著波數(shù)比的增大,附加面密度迅速衰減.(3)水下矩形板的附加質(zhì)量不可忽略,研究水下結(jié)構(gòu)的振動(dòng)與聲輻射必須考慮附加質(zhì)量的影響.

    [1] 沈惠明,趙德有.流固耦合振動(dòng)問(wèn)題的特征值解法[J].大連理工大學(xué)學(xué)報(bào),1990,30(3):369-371.

    [2] 吳 芳,趙德有.水對(duì)船舶與海洋建筑物結(jié)構(gòu)振動(dòng)影響的研究[J].中國(guó)海洋平臺(tái),2007,22(3):22-26.

    [3] K.H.Jeong.Free vibration of two identical circular plates coupled with bounded fluid[J].Journal of Sound and Vibration,2003,260:653-670.

    [4] K.H.Jeong,K.J.Kim.Hydroelastic vibration of acircular plate submerged in a bounded compressible fluid[J].Journal of Sound and Vibration,2005,283:153-172.

    [5] Kyeong Hoon Jeong,Gyu Mahn Lee,Tae Wan Kim. Free vibration analysis of a circular plate partially in contact with a liquid[J]. Journal of Sound and Vibration,2009,324:194-208.

    [6] Kenneth D.Frampton.Radiation efficiency of convected fluid-loaded plates[J].J.Acoust.Soc Am,2003,113(5):2 663-2 673.[7] 王盛基,楊慶山.流體環(huán)境中結(jié)構(gòu)附加質(zhì)量的計(jì)算[J].北方交通大學(xué)學(xué)報(bào),2003,27(3):40-43.

    [8] 錢(qián) 勤,黃玉盈,劉忠族.求附連水質(zhì)量的一種直接方法[J].力學(xué)與實(shí)踐,1996,18(5):19-21.

    [9] 于丹竹,黎 勝.附加質(zhì)量對(duì)于水中平面板振動(dòng)噪聲特性的影響[J].艦船科學(xué)技術(shù),2015,37(6):26-29.

    [10] 馬 燁,單雪雄.數(shù)值計(jì)算復(fù)雜外形物體附加質(zhì)量的新方法[J].計(jì)算機(jī)仿真,2007,24(5):75-78,113.

    [11] 劉 瑩,楊 衡,王毅娜,等.三維結(jié)構(gòu)水中剛彈耦合運(yùn)動(dòng)附加質(zhì)量特性分析[J].江蘇船舶,2014,31(5):5-9.

    [12] 李 明,尹云玉.水下航行體動(dòng)態(tài)響應(yīng)計(jì)算的附加質(zhì)量探討[J].導(dǎo)彈與航天運(yùn)載技術(shù),2008(4):16-18.

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

    [14] G.Xie,D.J.Thompson,C.J.C.Jones.The radiation efficiency of baffled plates and strips[J].Journal of Sound and Vibration,2005,280:181-209.

    [15] Fahy F.Sound and structural vibration:Radiation,transmission and response[M].London:Academic Press,1985.

    【責(zé)任編輯:蔣亞儒】

    Study on the attached mass of the rectangular plate in water

    REN Hui-juan1, YAO Zhan1, XUE Xiao-qing1, LIU Ting1, LEI Ye2

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

    Starting from the vibration response formula of the rectangular plate in water and according to the definition of the attached mass,the expression of attached mass is derived for a rectangular plate in water,and the numerical simulations are implemented.The results show:First, the lower the modal order,the bigger the attached mass.Second, whenkw/kmnw<10-1,the attached mass almost keeps a constant value.Whenkw/kmnw≈1,the attached mass appears a peak value except the (1,1) mode,then it decays fastly.Third,the attached mass of the rectangular plate in water can not be neglected,and it should be taken into account when vestigating its vibration and sound radiation.

    rectangular plate; attached mass; attached damping

    2016-09-27

    陜西省科技廳自然科學(xué)基金項(xiàng)目 (2011JM1017); 陜西省教育廳科學(xué)研究計(jì)劃項(xiàng)目(2013JK0618); 咸陽(yáng)師范學(xué)院專(zhuān)項(xiàng)科研基金項(xiàng)目(13XSYK016)

    任惠娟(1972-),女,陜西周至人,副教授,博士,研究方向:結(jié)構(gòu)振動(dòng)與聲輻射

    1000-5811(2016)06-0183-04

    O422.6

    A

    猜你喜歡
    面密度波數(shù)高階
    聲場(chǎng)波數(shù)積分截?cái)嗖〝?shù)自適應(yīng)選取方法
    低面密度PE/PP雙組分紡粘非織造布的制備及其性能
    毛紡科技(2022年10期)2022-10-31 03:56:48
    一種基于SOM神經(jīng)網(wǎng)絡(luò)中藥材分類(lèi)識(shí)別系統(tǒng)
    有限圖上高階Yamabe型方程的非平凡解
    高階各向異性Cahn-Hilliard-Navier-Stokes系統(tǒng)的弱解
    滾動(dòng)軸承壽命高階計(jì)算與應(yīng)用
    哈爾濱軸承(2020年1期)2020-11-03 09:16:02
    纖維層面密度的光學(xué)新算法
    水刺產(chǎn)品面密度均勻度波動(dòng)因素分析與自動(dòng)控制
    基于Bernstein多項(xiàng)式的配點(diǎn)法解高階常微分方程
    重磁異常解釋的歸一化局部波數(shù)法
    亚洲精品国产av蜜桃| 又大又黄又爽视频免费| 久久国产精品影院| 国产成人啪精品午夜网站| 久久久久久亚洲精品国产蜜桃av| 男人操女人黄网站| 在线看a的网站| 女性生殖器流出的白浆| 亚洲精品一二三| 赤兔流量卡办理| tube8黄色片| 色婷婷久久久亚洲欧美| 欧美国产精品va在线观看不卡| 高清欧美精品videossex| 国产成人a∨麻豆精品| 亚洲欧美一区二区三区久久| 国产免费一区二区三区四区乱码| 男女高潮啪啪啪动态图| e午夜精品久久久久久久| 国产成人一区二区三区免费视频网站 | av在线播放精品| 伦理电影免费视频| 日本猛色少妇xxxxx猛交久久| 首页视频小说图片口味搜索 | 亚洲专区中文字幕在线| 天天躁夜夜躁狠狠躁躁| 久久久久久久国产电影| 999久久久国产精品视频| 美国免费a级毛片| 亚洲成人手机| 久久国产亚洲av麻豆专区| 成人18禁高潮啪啪吃奶动态图| 可以免费在线观看a视频的电影网站| 十分钟在线观看高清视频www| 国产熟女午夜一区二区三区| 亚洲天堂av无毛| 亚洲精品日韩在线中文字幕| 欧美精品一区二区大全| 永久免费av网站大全| 国产高清不卡午夜福利| 亚洲av日韩精品久久久久久密 | 中文欧美无线码| 这个男人来自地球电影免费观看| 天天躁夜夜躁狠狠久久av| 如日韩欧美国产精品一区二区三区| 如日韩欧美国产精品一区二区三区| 久久性视频一级片| 久久九九热精品免费| 亚洲七黄色美女视频| 黄色一级大片看看| 亚洲国产日韩一区二区| av有码第一页| 免费观看av网站的网址| 满18在线观看网站| 老鸭窝网址在线观看| 天天躁夜夜躁狠狠久久av| 汤姆久久久久久久影院中文字幕| 只有这里有精品99| 亚洲人成电影免费在线| 亚洲av男天堂| 最新在线观看一区二区三区 | 精品第一国产精品| 国产片内射在线| 熟女少妇亚洲综合色aaa.| 亚洲国产中文字幕在线视频| 午夜久久久在线观看| 欧美精品av麻豆av| 久久99热这里只频精品6学生| 日本av手机在线免费观看| 黄色毛片三级朝国网站| 午夜91福利影院| 丝袜脚勾引网站| 国产精品久久久久久精品电影小说| 手机成人av网站| 老司机深夜福利视频在线观看 | 午夜免费鲁丝| 丰满迷人的少妇在线观看| 欧美久久黑人一区二区| 久久国产亚洲av麻豆专区| 婷婷色综合www| 丰满少妇做爰视频| 色婷婷久久久亚洲欧美| 狂野欧美激情性bbbbbb| 国产精品一二三区在线看| 少妇 在线观看| 免费在线观看影片大全网站 | 无遮挡黄片免费观看| 欧美成人精品欧美一级黄| 蜜桃国产av成人99| 久久久久久久国产电影| 97在线人人人人妻| 男女无遮挡免费网站观看| 亚洲av美国av| 十八禁网站网址无遮挡| 看免费成人av毛片| 久久女婷五月综合色啪小说| 黑丝袜美女国产一区| 欧美激情高清一区二区三区| 精品国产乱码久久久久久男人| 黄色一级大片看看| 九草在线视频观看| 国产男女超爽视频在线观看| 你懂的网址亚洲精品在线观看| 久久99一区二区三区| 天天操日日干夜夜撸| 一区二区三区精品91| 午夜老司机福利片| 欧美另类一区| 少妇猛男粗大的猛烈进出视频| 久久精品国产综合久久久| 91九色精品人成在线观看| 男人爽女人下面视频在线观看| 欧美精品亚洲一区二区| 老汉色∧v一级毛片| 亚洲少妇的诱惑av| 午夜免费男女啪啪视频观看| kizo精华| 人人妻人人爽人人添夜夜欢视频| 黑丝袜美女国产一区| 久久人妻熟女aⅴ| 久久天躁狠狠躁夜夜2o2o | 精品国产国语对白av| 亚洲自偷自拍图片 自拍| 手机成人av网站| 国产精品99久久99久久久不卡| 一区二区av电影网| 精品福利永久在线观看| 免费人妻精品一区二区三区视频| 日韩伦理黄色片| 免费不卡黄色视频| 欧美日韩亚洲国产一区二区在线观看 | 国产片特级美女逼逼视频| 一二三四在线观看免费中文在| 亚洲五月色婷婷综合| 嫩草影视91久久| 久热这里只有精品99| 尾随美女入室| 叶爱在线成人免费视频播放| 亚洲七黄色美女视频| 婷婷色综合大香蕉| 777米奇影视久久| 一级毛片女人18水好多 | 亚洲第一青青草原| 看免费成人av毛片| 国产精品久久久久久人妻精品电影 | 亚洲激情五月婷婷啪啪| 欧美xxⅹ黑人| 91国产中文字幕| 最新在线观看一区二区三区 | 热99国产精品久久久久久7| 国产一级毛片在线| 亚洲精品一区蜜桃| 国产精品二区激情视频| 午夜福利在线免费观看网站| 久久久久久久久久久久大奶| 婷婷色麻豆天堂久久| 婷婷色综合大香蕉| 成年女人毛片免费观看观看9 | 亚洲免费av在线视频| 在线观看国产h片| 欧美乱码精品一区二区三区| 秋霞在线观看毛片| 亚洲专区中文字幕在线| 男女午夜视频在线观看| 宅男免费午夜| 男女免费视频国产| 交换朋友夫妻互换小说| 国产成人av教育| 2018国产大陆天天弄谢| 精品卡一卡二卡四卡免费| 国产精品一区二区精品视频观看| 精品人妻1区二区| kizo精华| 男的添女的下面高潮视频| 人人妻人人澡人人看| 亚洲精品久久午夜乱码| 九色亚洲精品在线播放| 久久久久久久大尺度免费视频| 久久亚洲精品不卡| 精品国产一区二区三区久久久樱花| 久久九九热精品免费| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区久久| 建设人人有责人人尽责人人享有的| 天天影视国产精品| 国产成人精品在线电影| 成年人黄色毛片网站| 免费久久久久久久精品成人欧美视频| 国产成人精品久久二区二区免费| 一边摸一边做爽爽视频免费| 丁香六月天网| 丝袜美腿诱惑在线| 男人操女人黄网站| 久久久久精品国产欧美久久久 | e午夜精品久久久久久久| 9色porny在线观看| 99国产精品一区二区三区| 丝瓜视频免费看黄片| 精品人妻熟女毛片av久久网站| 色精品久久人妻99蜜桃| 国产成人免费观看mmmm| av国产精品久久久久影院| 国产亚洲一区二区精品| 国产真人三级小视频在线观看| 日韩熟女老妇一区二区性免费视频| 日韩大码丰满熟妇| 亚洲国产中文字幕在线视频| 天天躁日日躁夜夜躁夜夜| 亚洲国产欧美日韩在线播放| 亚洲第一青青草原| 男女边摸边吃奶| 成人国产av品久久久| 日韩电影二区| 捣出白浆h1v1| 国产精品国产三级专区第一集| svipshipincom国产片| av不卡在线播放| www.自偷自拍.com| 黄色视频不卡| 只有这里有精品99| 多毛熟女@视频| 一本大道久久a久久精品| 999精品在线视频| 精品少妇黑人巨大在线播放| 纵有疾风起免费观看全集完整版| 欧美黄色淫秽网站| 亚洲中文日韩欧美视频| 亚洲 国产 在线| 久久中文字幕一级| 免费高清在线观看视频在线观看| 又粗又硬又长又爽又黄的视频| 国产精品久久久久成人av| 久热爱精品视频在线9| 巨乳人妻的诱惑在线观看| 你懂的网址亚洲精品在线观看| 极品人妻少妇av视频| 国产成人精品久久二区二区免费| 亚洲av综合色区一区| 人妻一区二区av| 天天操日日干夜夜撸| 一本—道久久a久久精品蜜桃钙片| 一区二区三区乱码不卡18| 人人妻,人人澡人人爽秒播 | 麻豆国产av国片精品| 久久国产精品人妻蜜桃| 水蜜桃什么品种好| 日韩欧美一区视频在线观看| 99国产综合亚洲精品| 人人妻人人澡人人看| 久久国产亚洲av麻豆专区| 国产高清视频在线播放一区 | 国产欧美日韩一区二区三区在线| 91精品三级在线观看| 美女大奶头黄色视频| av福利片在线| 日韩中文字幕视频在线看片| 又黄又粗又硬又大视频| 一个人免费看片子| 大陆偷拍与自拍| 女人精品久久久久毛片| 亚洲,一卡二卡三卡| 激情视频va一区二区三区| 亚洲视频免费观看视频| 一区二区三区精品91| 久久影院123| 日韩免费高清中文字幕av| 精品亚洲乱码少妇综合久久| 日韩,欧美,国产一区二区三区| 亚洲色图 男人天堂 中文字幕| 脱女人内裤的视频| 一级毛片我不卡| a级毛片在线看网站| 欧美日韩一级在线毛片| 精品少妇内射三级| 女人爽到高潮嗷嗷叫在线视频| 亚洲精品乱久久久久久| 天堂8中文在线网| 国产精品免费视频内射| 国产精品久久久av美女十八| 大话2 男鬼变身卡| 狂野欧美激情性xxxx| 婷婷色综合www| 国产精品香港三级国产av潘金莲 | 亚洲,一卡二卡三卡| 韩国精品一区二区三区| 国产成人精品在线电影| 亚洲国产av新网站| 99久久99久久久精品蜜桃| 国产在线观看jvid| 久久久久久久久久久久大奶| 一本色道久久久久久精品综合| 国产精品久久久久成人av| 精品视频人人做人人爽| 日韩av在线免费看完整版不卡| 亚洲欧洲精品一区二区精品久久久| 久久精品久久久久久久性| 欧美av亚洲av综合av国产av| 亚洲精品久久久久久婷婷小说| 成人18禁高潮啪啪吃奶动态图| 飞空精品影院首页| 日本黄色日本黄色录像| 大香蕉久久成人网| 黑人欧美特级aaaaaa片| 国产黄色视频一区二区在线观看| av视频免费观看在线观看| 最黄视频免费看| 久久精品国产综合久久久| 日本91视频免费播放| 男的添女的下面高潮视频| 各种免费的搞黄视频| 老司机影院成人| 国产91精品成人一区二区三区 | 亚洲av在线观看美女高潮| 欧美日韩综合久久久久久| 成人三级做爰电影| 一边摸一边抽搐一进一出视频| 欧美黑人精品巨大| 麻豆av在线久日| 国产亚洲欧美精品永久| 国产精品国产三级国产专区5o| 啦啦啦 在线观看视频| 国产精品秋霞免费鲁丝片| netflix在线观看网站| 美女中出高潮动态图| 国产福利在线免费观看视频| 一级片'在线观看视频| 久久国产精品影院| 亚洲,欧美,日韩| 99精品久久久久人妻精品| av电影中文网址| 欧美在线黄色| 亚洲激情五月婷婷啪啪| 日本av手机在线免费观看| 亚洲 国产 在线| 亚洲av日韩在线播放| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品在线美女| 欧美激情极品国产一区二区三区| 欧美日韩福利视频一区二区| 精品久久久精品久久久| 日韩大码丰满熟妇| 人人澡人人妻人| 亚洲少妇的诱惑av| 国产成人精品久久二区二区91| 一级毛片电影观看| 侵犯人妻中文字幕一二三四区| 国产伦理片在线播放av一区| 久久精品久久久久久噜噜老黄| 在线看a的网站| 国产亚洲欧美精品永久| 一级黄片播放器| 国产一区二区激情短视频 | 日本午夜av视频| 精品国产国语对白av| 国产主播在线观看一区二区 | 水蜜桃什么品种好| 欧美黄色片欧美黄色片| 国产成人影院久久av| 国精品久久久久久国模美| 涩涩av久久男人的天堂| 少妇裸体淫交视频免费看高清 | 曰老女人黄片| 亚洲精品一二三| 日本一区二区免费在线视频| av不卡在线播放| 国产午夜精品一二区理论片| 国产精品成人在线| 精品卡一卡二卡四卡免费| h视频一区二区三区| 婷婷色麻豆天堂久久| 热re99久久国产66热| 免费黄频网站在线观看国产| 丝袜人妻中文字幕| 成人影院久久| 女人被躁到高潮嗷嗷叫费观| 国产精品.久久久| 中文精品一卡2卡3卡4更新| 国产精品久久久久久精品古装| 老司机影院毛片| 久久久精品94久久精品| 黄色 视频免费看| 男女之事视频高清在线观看 | 亚洲专区国产一区二区| av有码第一页| 国产激情久久老熟女| 黄网站色视频无遮挡免费观看| 国产伦人伦偷精品视频| 脱女人内裤的视频| 日韩电影二区| 午夜av观看不卡| 一本久久精品| 人人妻人人爽人人添夜夜欢视频| 国产成人免费观看mmmm| 一级a爱视频在线免费观看| 成在线人永久免费视频| 欧美精品av麻豆av| 午夜精品国产一区二区电影| 亚洲欧美精品自产自拍| 精品第一国产精品| 精品人妻1区二区| 亚洲精品久久久久久婷婷小说| 亚洲av片天天在线观看| 人妻 亚洲 视频| 国产精品秋霞免费鲁丝片| 人妻人人澡人人爽人人| 午夜福利,免费看| 国产精品人妻久久久影院| 国产亚洲一区二区精品| 一级毛片黄色毛片免费观看视频| 国产人伦9x9x在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲精品日韩在线中文字幕| 精品高清国产在线一区| 精品一品国产午夜福利视频| 精品国产乱码久久久久久小说| 下体分泌物呈黄色| 90打野战视频偷拍视频| 高清av免费在线| 91字幕亚洲| 亚洲一区二区三区欧美精品| 国产成人精品久久二区二区免费| 成人午夜精彩视频在线观看| 午夜老司机福利片| 2021少妇久久久久久久久久久| 国产在线视频一区二区| 婷婷丁香在线五月| 亚洲国产成人一精品久久久| 人成视频在线观看免费观看| 久久99热这里只频精品6学生| 精品一区二区三卡| 中文欧美无线码| 2021少妇久久久久久久久久久| 欧美日本中文国产一区发布| 桃花免费在线播放| 免费在线观看影片大全网站 | 国产成人免费观看mmmm| 国产精品国产三级国产专区5o| 色网站视频免费| 国产一区二区 视频在线| 国产精品久久久av美女十八| 亚洲一区中文字幕在线| 亚洲熟女毛片儿| 51午夜福利影视在线观看| 91老司机精品| 精品视频人人做人人爽| 国产精品秋霞免费鲁丝片| 一级黄色大片毛片| 国产有黄有色有爽视频| 亚洲国产精品一区三区| 日本欧美视频一区| 精品国产乱码久久久久久小说| 国产片特级美女逼逼视频| 黄色一级大片看看| 国产高清videossex| 涩涩av久久男人的天堂| 久久天躁狠狠躁夜夜2o2o | 久久精品国产亚洲av涩爱| 免费久久久久久久精品成人欧美视频| 18禁观看日本| 亚洲国产精品一区二区三区在线| 少妇猛男粗大的猛烈进出视频| av视频免费观看在线观看| 精品一区二区三区四区五区乱码 | 午夜免费男女啪啪视频观看| 十分钟在线观看高清视频www| 国产视频一区二区在线看| 性色av乱码一区二区三区2| 高潮久久久久久久久久久不卡| avwww免费| 各种免费的搞黄视频| 好男人电影高清在线观看| 看免费成人av毛片| av一本久久久久| 亚洲伊人色综图| 成人黄色视频免费在线看| 9热在线视频观看99| 你懂的网址亚洲精品在线观看| 91精品国产国语对白视频| 在线精品无人区一区二区三| 亚洲精品国产色婷婷电影| h视频一区二区三区| 国产97色在线日韩免费| 欧美+亚洲+日韩+国产| 亚洲天堂av无毛| 精品国产超薄肉色丝袜足j| 欧美精品一区二区大全| 欧美乱码精品一区二区三区| 五月开心婷婷网| 欧美日韩视频高清一区二区三区二| 日韩欧美一区视频在线观看| 18禁裸乳无遮挡动漫免费视频| 丝瓜视频免费看黄片| 视频区图区小说| 中文字幕制服av| 超色免费av| 国产亚洲av高清不卡| 亚洲一码二码三码区别大吗| 咕卡用的链子| 亚洲av电影在线进入| 日日摸夜夜添夜夜爱| 中文字幕人妻丝袜一区二区| 国产亚洲午夜精品一区二区久久| videos熟女内射| 精品一区二区三区四区五区乱码 | 欧美人与善性xxx| 成人18禁高潮啪啪吃奶动态图| 久久国产精品男人的天堂亚洲| 日本一区二区免费在线视频| 欧美成人精品欧美一级黄| 大香蕉久久成人网| 久久国产精品影院| 日韩av不卡免费在线播放| 黄色视频在线播放观看不卡| 日本a在线网址| 欧美国产精品va在线观看不卡| 五月开心婷婷网| 国产免费视频播放在线视频| 男女床上黄色一级片免费看| 久久ye,这里只有精品| 另类亚洲欧美激情| 国产亚洲午夜精品一区二区久久| 欧美97在线视频| 日韩av免费高清视频| 亚洲国产看品久久| 精品欧美一区二区三区在线| 国产一区有黄有色的免费视频| 99久久精品国产亚洲精品| 伊人久久大香线蕉亚洲五| 观看av在线不卡| 久久久久久久大尺度免费视频| 成人黄色视频免费在线看| 热re99久久精品国产66热6| 日韩大码丰满熟妇| 一边摸一边做爽爽视频免费| 麻豆av在线久日| 一二三四在线观看免费中文在| 视频在线观看一区二区三区| 久久国产精品影院| 视频在线观看一区二区三区| 国产精品99久久99久久久不卡| 交换朋友夫妻互换小说| 777久久人妻少妇嫩草av网站| cao死你这个sao货| 日韩 亚洲 欧美在线| 1024视频免费在线观看| 侵犯人妻中文字幕一二三四区| 汤姆久久久久久久影院中文字幕| 国产亚洲av高清不卡| 超色免费av| 啦啦啦啦在线视频资源| 国产精品成人在线| 美女视频免费永久观看网站| 最黄视频免费看| 成年av动漫网址| 女警被强在线播放| 亚洲精品一区蜜桃| 两人在一起打扑克的视频| 日本91视频免费播放| 一本色道久久久久久精品综合| 天天躁狠狠躁夜夜躁狠狠躁| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品香港三级国产av潘金莲 | 一个人免费看片子| 国产成人欧美在线观看 | 老汉色av国产亚洲站长工具| 天堂8中文在线网| 99国产精品免费福利视频| 日本av手机在线免费观看| 国产三级黄色录像| 美女福利国产在线| 最近中文字幕2019免费版| 一本大道久久a久久精品| 欧美在线一区亚洲| 亚洲欧美日韩另类电影网站| 大香蕉久久成人网| 亚洲国产成人一精品久久久| 亚洲国产av新网站| 欧美av亚洲av综合av国产av| 日本猛色少妇xxxxx猛交久久| 18禁国产床啪视频网站| 久久青草综合色| 看十八女毛片水多多多| 99九九在线精品视频| 肉色欧美久久久久久久蜜桃| 国产老妇伦熟女老妇高清| 电影成人av| 久久精品久久久久久久性| 国产熟女午夜一区二区三区| 亚洲av成人不卡在线观看播放网 | 成人国产av品久久久| 精品一区二区三区四区五区乱码 | 久久国产精品大桥未久av| 最近手机中文字幕大全| 亚洲欧美一区二区三区黑人| 五月天丁香电影| 国产成人影院久久av| 热re99久久精品国产66热6| 超色免费av| 黄色视频在线播放观看不卡| 满18在线观看网站| 久热爱精品视频在线9| 久久性视频一级片| 男人舔女人的私密视频| 国产男人的电影天堂91| 十八禁高潮呻吟视频| 中文精品一卡2卡3卡4更新| 两个人免费观看高清视频| 黄片播放在线免费| 午夜福利免费观看在线| 男男h啪啪无遮挡| 亚洲av日韩在线播放| 午夜激情av网站| 国产熟女午夜一区二区三区| 少妇粗大呻吟视频| 成人国语在线视频|