• <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ù)法
    亚洲美女视频黄频| 男女边吃奶边做爰视频| 9191精品国产免费久久| 女的被弄到高潮叫床怎么办| 宅男免费午夜| 国产极品天堂在线| 亚洲内射少妇av| 99九九在线精品视频| 久久久久久久大尺度免费视频| 国产不卡av网站在线观看| 亚洲国产精品一区二区三区在线| 国产亚洲欧美精品永久| 国产在视频线精品| av黄色大香蕉| 人妻人人澡人人爽人人| 女人被躁到高潮嗷嗷叫费观| 国产精品不卡视频一区二区| 国产永久视频网站| 亚洲av福利一区| 在线观看免费高清a一片| 精品一区在线观看国产| 亚洲内射少妇av| av有码第一页| 国产成人91sexporn| 全区人妻精品视频| 久久久久久久亚洲中文字幕| 久久久欧美国产精品| 国产高清三级在线| 亚洲,欧美,日韩| 狠狠婷婷综合久久久久久88av| 亚洲中文av在线| 亚洲第一区二区三区不卡| 两性夫妻黄色片 | 免费高清在线观看视频在线观看| 99热国产这里只有精品6| 国产精品成人在线| 亚洲综合色网址| 啦啦啦中文免费视频观看日本| 老司机亚洲免费影院| 国产 精品1| 亚洲国产精品国产精品| 国产女主播在线喷水免费视频网站| 黄色毛片三级朝国网站| 丝袜人妻中文字幕| 黄色配什么色好看| 精品亚洲乱码少妇综合久久| 夜夜骑夜夜射夜夜干| 精品久久久久久电影网| 人成视频在线观看免费观看| 性高湖久久久久久久久免费观看| 国产精品久久久av美女十八| 亚洲精品国产av蜜桃| 午夜激情av网站| 美女国产高潮福利片在线看| av黄色大香蕉| 亚洲精品一二三| 色94色欧美一区二区| 成年人免费黄色播放视频| 国产免费福利视频在线观看| 精品一区二区三区四区五区乱码 | 宅男免费午夜| 大陆偷拍与自拍| 亚洲精品久久久久久婷婷小说| 少妇人妻久久综合中文| 亚洲色图 男人天堂 中文字幕 | 成人黄色视频免费在线看| 99re6热这里在线精品视频| 国产午夜精品一二区理论片| 欧美3d第一页| 国产一级毛片在线| 纯流量卡能插随身wifi吗| 久久久久久人人人人人| 一级爰片在线观看| 秋霞伦理黄片| 九草在线视频观看| 深夜精品福利| 欧美激情国产日韩精品一区| 看十八女毛片水多多多| 色5月婷婷丁香| 丝袜喷水一区| 欧美成人精品欧美一级黄| 99热全是精品| 久久久久久久国产电影| 伊人久久国产一区二区| 国产亚洲av片在线观看秒播厂| 女人被躁到高潮嗷嗷叫费观| 黄色一级大片看看| 国产熟女午夜一区二区三区| 久久久久久久精品精品| 国产无遮挡羞羞视频在线观看| 午夜福利乱码中文字幕| av网站免费在线观看视频| 中国三级夫妇交换| 一本大道久久a久久精品| 天天躁夜夜躁狠狠躁躁| 国产一区二区激情短视频 | 久久久欧美国产精品| av.在线天堂| 亚洲美女视频黄频| 免费观看无遮挡的男女| 亚洲精品国产av蜜桃| 亚洲精品美女久久久久99蜜臀 | 成人黄色视频免费在线看| 国产片内射在线| 男人爽女人下面视频在线观看| 纯流量卡能插随身wifi吗| 尾随美女入室| 免费高清在线观看视频在线观看| 成人影院久久| 丝袜脚勾引网站| av在线老鸭窝| 欧美最新免费一区二区三区| 亚洲av.av天堂| 国产精品国产三级国产av玫瑰| 在线观看免费视频网站a站| 国产一区二区三区av在线| 免费大片黄手机在线观看| 极品少妇高潮喷水抽搐| 热re99久久精品国产66热6| 99热全是精品| 亚洲伊人色综图| 国产亚洲av片在线观看秒播厂| 99久久人妻综合| 国产高清三级在线| 国产高清三级在线| 国产精品蜜桃在线观看| 美女国产高潮福利片在线看| 国产成人免费无遮挡视频| 久久人人爽av亚洲精品天堂| av不卡在线播放| 欧美丝袜亚洲另类| 成年女人在线观看亚洲视频| 国产一区二区三区综合在线观看 | 少妇人妻久久综合中文| 伦理电影免费视频| 久久99精品国语久久久| 18禁在线无遮挡免费观看视频| 黑人巨大精品欧美一区二区蜜桃 | 午夜免费观看性视频| 亚洲欧美成人精品一区二区| 精品熟女少妇av免费看| 日本vs欧美在线观看视频| 我要看黄色一级片免费的| 在线观看免费视频网站a站| 日本与韩国留学比较| 日日摸夜夜添夜夜爱| 丁香六月天网| 丝袜人妻中文字幕| 中文字幕av电影在线播放| 国产日韩欧美亚洲二区| 成人毛片60女人毛片免费| 男女国产视频网站| 国产在线一区二区三区精| 大话2 男鬼变身卡| 欧美日韩视频高清一区二区三区二| 午夜精品国产一区二区电影| 亚洲av日韩在线播放| 中文字幕最新亚洲高清| 九九爱精品视频在线观看| 成人亚洲欧美一区二区av| 少妇的逼好多水| 99热这里只有是精品在线观看| 精品酒店卫生间| 亚洲国产欧美在线一区| 妹子高潮喷水视频| 久久人妻熟女aⅴ| 免费av中文字幕在线| 国产精品久久久久久久久免| 大话2 男鬼变身卡| 91在线精品国自产拍蜜月| 97在线人人人人妻| 亚洲激情五月婷婷啪啪| 日韩欧美精品免费久久| 精品卡一卡二卡四卡免费| 99久久人妻综合| 寂寞人妻少妇视频99o| 国产亚洲最大av| 丝袜在线中文字幕| 国产免费一级a男人的天堂| 两个人看的免费小视频| 乱码一卡2卡4卡精品| 日韩av在线免费看完整版不卡| 久久青草综合色| 精品亚洲成国产av| 免费人妻精品一区二区三区视频| 亚洲精品国产av成人精品| 久久久精品免费免费高清| 制服诱惑二区| a 毛片基地| 国产欧美亚洲国产| 多毛熟女@视频| 中文字幕亚洲精品专区| 交换朋友夫妻互换小说| 巨乳人妻的诱惑在线观看| 80岁老熟妇乱子伦牲交| 中文欧美无线码| 精品卡一卡二卡四卡免费| 日韩欧美一区视频在线观看| 国产欧美日韩综合在线一区二区| 中文乱码字字幕精品一区二区三区| 宅男免费午夜| 中文精品一卡2卡3卡4更新| 韩国高清视频一区二区三区| 全区人妻精品视频| 日本爱情动作片www.在线观看| 乱码一卡2卡4卡精品| 汤姆久久久久久久影院中文字幕| 欧美 亚洲 国产 日韩一| 欧美老熟妇乱子伦牲交| 日韩三级伦理在线观看| 成年av动漫网址| 成人毛片a级毛片在线播放| 日日爽夜夜爽网站| 日韩欧美一区视频在线观看| 各种免费的搞黄视频| 国产亚洲一区二区精品| 男女免费视频国产| 亚洲成人手机| 欧美老熟妇乱子伦牲交| 久久人人爽人人爽人人片va| 久久精品国产a三级三级三级| 中文天堂在线官网| 大片免费播放器 马上看| 久久精品aⅴ一区二区三区四区 | 久久久国产一区二区| 香蕉精品网在线| 亚洲精品aⅴ在线观看| 亚洲精品av麻豆狂野| 美女中出高潮动态图| 欧美日韩视频精品一区| a级毛片黄视频| 国产在线视频一区二区| 国产国拍精品亚洲av在线观看| 三级国产精品片| 欧美激情国产日韩精品一区| 满18在线观看网站| 欧美成人精品欧美一级黄| 亚洲精品久久成人aⅴ小说| 久久精品aⅴ一区二区三区四区 | 亚洲综合精品二区| 国产欧美亚洲国产| 夫妻午夜视频| av在线app专区| 国产精品无大码| 国产精品嫩草影院av在线观看| www.av在线官网国产| 午夜久久久在线观看| 免费黄网站久久成人精品| 亚洲一区二区三区欧美精品| a级毛片黄视频| 亚洲经典国产精华液单| 天美传媒精品一区二区| 久久久欧美国产精品| 黑人高潮一二区| 久久久久久久大尺度免费视频| 日韩av不卡免费在线播放| videossex国产| 制服丝袜香蕉在线| 久久精品国产综合久久久 | 性色avwww在线观看| 亚洲性久久影院| 亚洲精品色激情综合| 人妻少妇偷人精品九色| 国产又色又爽无遮挡免| 妹子高潮喷水视频| 日韩一区二区三区影片| h视频一区二区三区| 国产毛片在线视频| 国产黄频视频在线观看| 9色porny在线观看| 婷婷色av中文字幕| 尾随美女入室| 如日韩欧美国产精品一区二区三区| 视频中文字幕在线观看| 国产精品国产三级国产专区5o| 在线观看美女被高潮喷水网站| 女性被躁到高潮视频| 波野结衣二区三区在线| 久久韩国三级中文字幕| 免费看光身美女| 黑人高潮一二区| 国产亚洲精品久久久com| 少妇人妻 视频| 久久国产亚洲av麻豆专区| 美女视频免费永久观看网站| 国产一区有黄有色的免费视频| 永久免费av网站大全| 夫妻午夜视频| 三级国产精品片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 美国免费a级毛片| 亚洲图色成人| 女人久久www免费人成看片| 一本色道久久久久久精品综合| 日韩成人伦理影院| 90打野战视频偷拍视频| 亚洲内射少妇av| 久久精品久久久久久久性| 黑丝袜美女国产一区| 日本wwww免费看| 国产精品久久久久成人av| 亚洲人成网站在线观看播放| 国产精品无大码| 亚洲高清免费不卡视频| 免费黄频网站在线观看国产| 天堂中文最新版在线下载| 日韩中文字幕视频在线看片| 中文字幕免费在线视频6| av在线观看视频网站免费| 黄色视频在线播放观看不卡| 一区二区av电影网| 中文字幕免费在线视频6| 香蕉国产在线看| 97精品久久久久久久久久精品| 国产免费一级a男人的天堂| 咕卡用的链子| 香蕉丝袜av| 国产免费一区二区三区四区乱码| 女性生殖器流出的白浆| 美女福利国产在线| 亚洲国产日韩一区二区| 亚洲精品av麻豆狂野| 欧美激情国产日韩精品一区| 欧美性感艳星| 亚洲久久久国产精品| 狂野欧美激情性xxxx在线观看| 精品一品国产午夜福利视频| 亚洲内射少妇av| 亚洲第一av免费看| 久久久久视频综合| 国产精品久久久av美女十八| 日日撸夜夜添| 国产一区有黄有色的免费视频| 日本wwww免费看| 午夜精品国产一区二区电影| 国产国语露脸激情在线看| 女人被躁到高潮嗷嗷叫费观| 国产欧美日韩综合在线一区二区| 免费黄网站久久成人精品| 美女主播在线视频| 国产成人欧美| 99热国产这里只有精品6| 香蕉精品网在线| 午夜av观看不卡| 国产精品国产三级国产专区5o| 91午夜精品亚洲一区二区三区| 久久久久网色| 边亲边吃奶的免费视频| 一二三四在线观看免费中文在 | 日日爽夜夜爽网站| 国产白丝娇喘喷水9色精品| 亚洲精品久久午夜乱码| 不卡视频在线观看欧美| 国产av国产精品国产| 亚洲精品美女久久久久99蜜臀 | av.在线天堂| 最近手机中文字幕大全| 在线看a的网站| 色婷婷久久久亚洲欧美| 国产永久视频网站| 国产精品国产三级国产专区5o| 丝袜人妻中文字幕| 超碰97精品在线观看| 欧美亚洲 丝袜 人妻 在线| 国产日韩欧美视频二区| 欧美日韩成人在线一区二区| 欧美3d第一页| 成人手机av| 免费看av在线观看网站| 国产高清不卡午夜福利| 国产综合精华液| 侵犯人妻中文字幕一二三四区| 欧美精品高潮呻吟av久久| 人成视频在线观看免费观看| 精品久久国产蜜桃| 亚洲精品成人av观看孕妇| 国产精品人妻久久久久久| 日产精品乱码卡一卡2卡三| 久久国产亚洲av麻豆专区| 亚洲欧美日韩另类电影网站| 日本欧美视频一区| 欧美日韩视频精品一区| 18禁在线无遮挡免费观看视频| 女的被弄到高潮叫床怎么办| 免费观看av网站的网址| 精品久久国产蜜桃| 三上悠亚av全集在线观看| 欧美最新免费一区二区三区| 国产亚洲午夜精品一区二区久久| 麻豆精品久久久久久蜜桃| 国产永久视频网站| 丰满饥渴人妻一区二区三| 天天躁夜夜躁狠狠久久av| 日韩 亚洲 欧美在线| 亚洲精品日韩在线中文字幕| 国产精品久久久久久精品古装| 欧美人与性动交α欧美精品济南到 | 制服诱惑二区| 国产爽快片一区二区三区| 精品国产乱码久久久久久小说| 高清毛片免费看| 日本wwww免费看| 国产亚洲精品第一综合不卡 | 两个人免费观看高清视频| 五月伊人婷婷丁香| 精品少妇内射三级| 亚洲精品色激情综合| 国内精品宾馆在线| 久久精品熟女亚洲av麻豆精品| 亚洲av中文av极速乱| 欧美性感艳星| 亚洲美女搞黄在线观看| 亚洲第一区二区三区不卡| 免费久久久久久久精品成人欧美视频 | 久久国内精品自在自线图片| 少妇的逼水好多| 亚洲国产精品专区欧美| 91在线精品国自产拍蜜月| 日日啪夜夜爽| 国产一区二区三区av在线| 国产老妇伦熟女老妇高清| 国产成人精品无人区| 欧美xxxx性猛交bbbb| 欧美日韩成人在线一区二区| 18禁观看日本| 日韩中文字幕视频在线看片| 免费看不卡的av| 日本av手机在线免费观看| 九草在线视频观看| 一区二区三区乱码不卡18| 高清黄色对白视频在线免费看| 亚洲人与动物交配视频| 哪个播放器可以免费观看大片| 国产欧美日韩一区二区三区在线| 七月丁香在线播放| 久热这里只有精品99| 久久午夜综合久久蜜桃| 精品酒店卫生间| 精品国产国语对白av| 你懂的网址亚洲精品在线观看| 国产老妇伦熟女老妇高清| freevideosex欧美| 国产一级毛片在线| av免费在线看不卡| 久久婷婷青草| 日韩成人av中文字幕在线观看| 国产在线免费精品| 亚洲四区av| 这个男人来自地球电影免费观看 | 久久人人爽人人片av| 欧美日韩视频精品一区| 男女午夜视频在线观看 | 欧美日本中文国产一区发布| av电影中文网址| av网站免费在线观看视频| 激情视频va一区二区三区| 观看av在线不卡| 亚洲伊人色综图| 大片电影免费在线观看免费| 久久久久久久久久成人| 亚洲高清免费不卡视频| 美女大奶头黄色视频| 激情视频va一区二区三区| 搡女人真爽免费视频火全软件| 国产精品麻豆人妻色哟哟久久| 一级a做视频免费观看| a级片在线免费高清观看视频| 黑人欧美特级aaaaaa片| 国产免费现黄频在线看| 亚洲欧美中文字幕日韩二区| 国产黄频视频在线观看| 欧美成人午夜免费资源| 在线观看三级黄色| 日韩成人伦理影院| 精品视频人人做人人爽| 91午夜精品亚洲一区二区三区| 国产亚洲一区二区精品| 国产成人91sexporn| 日韩一区二区视频免费看| 国产精品女同一区二区软件| 99精国产麻豆久久婷婷| 在线免费观看不下载黄p国产| 少妇被粗大猛烈的视频| 欧美97在线视频| 免费大片18禁| 99国产综合亚洲精品| 中文乱码字字幕精品一区二区三区| 久久99精品国语久久久| av在线app专区| 丰满乱子伦码专区| 国产白丝娇喘喷水9色精品| 久久久久网色| 五月伊人婷婷丁香| 午夜影院在线不卡| 国产成人av激情在线播放| 精品国产露脸久久av麻豆| 少妇人妻精品综合一区二区| 一区二区av电影网| 成人午夜精彩视频在线观看| 日韩成人av中文字幕在线观看| 在线天堂最新版资源| 亚洲精品成人av观看孕妇| 亚洲欧洲国产日韩| 亚洲第一av免费看| 亚洲av电影在线观看一区二区三区| 少妇高潮的动态图| 女性被躁到高潮视频| 制服诱惑二区| 亚洲第一av免费看| 日本爱情动作片www.在线观看| av.在线天堂| a级毛色黄片| 日本wwww免费看| 亚洲精品久久成人aⅴ小说| av国产精品久久久久影院| 国产免费福利视频在线观看| 国产成人欧美| 中文精品一卡2卡3卡4更新| 91精品三级在线观看| 国产在视频线精品| 黄片无遮挡物在线观看| 国产亚洲精品第一综合不卡 | 性高湖久久久久久久久免费观看| 2021少妇久久久久久久久久久| 国产精品国产三级专区第一集| 伊人久久国产一区二区| 成人18禁高潮啪啪吃奶动态图| 久久久国产精品麻豆| 色婷婷久久久亚洲欧美| 中文字幕人妻熟女乱码| 国产又爽黄色视频| 观看美女的网站| 国产免费福利视频在线观看| 黄色 视频免费看| 欧美精品国产亚洲| 国产在线免费精品| 国产黄色视频一区二区在线观看| 人人妻人人澡人人看| 国产av一区二区精品久久| 啦啦啦视频在线资源免费观看| 成人毛片a级毛片在线播放| 国产一级毛片在线| 高清av免费在线| 少妇高潮的动态图| 51国产日韩欧美| 美女主播在线视频| 精品亚洲成a人片在线观看| 我的女老师完整版在线观看| 一级爰片在线观看| 日韩av免费高清视频| 免费在线观看完整版高清| 精品一区二区三区四区五区乱码 | 日日撸夜夜添| 亚洲精品国产av蜜桃| 国产成人91sexporn| 亚洲精品久久成人aⅴ小说| 天堂8中文在线网| 久久久久国产网址| 亚洲欧美清纯卡通| 亚洲精品视频女| 少妇 在线观看| 日韩伦理黄色片| 国产精品久久久av美女十八| 国产亚洲最大av| 日韩成人伦理影院| av天堂久久9| 亚洲精华国产精华液的使用体验| 黑人欧美特级aaaaaa片| 亚洲成色77777| 久久影院123| 久久久久久久久久久免费av| 欧美日韩视频高清一区二区三区二| 黄色一级大片看看| 咕卡用的链子| 九九在线视频观看精品| 久久青草综合色| videosex国产| 五月天丁香电影| 国产极品天堂在线| 亚洲欧美一区二区三区国产| 精品熟女少妇av免费看| 又黄又粗又硬又大视频| 99国产综合亚洲精品| 熟女人妻精品中文字幕| 五月伊人婷婷丁香| 国产无遮挡羞羞视频在线观看| 欧美精品人与动牲交sv欧美| 麻豆精品久久久久久蜜桃| 69精品国产乱码久久久| 伊人久久国产一区二区| 国产欧美另类精品又又久久亚洲欧美| 日韩电影二区| 欧美日韩视频高清一区二区三区二| 成人午夜精彩视频在线观看| 免费观看无遮挡的男女| 亚洲精品aⅴ在线观看| 激情视频va一区二区三区| 日韩中字成人| 乱码一卡2卡4卡精品| 成人影院久久| 久久久久久久久久人人人人人人| 国产老妇伦熟女老妇高清| 蜜臀久久99精品久久宅男| 又黄又粗又硬又大视频| 精品熟女少妇av免费看| 日本欧美国产在线视频| 国产极品天堂在线| 内地一区二区视频在线| 在线天堂最新版资源| 女人精品久久久久毛片| 国产深夜福利视频在线观看| 亚洲精品,欧美精品| 啦啦啦视频在线资源免费观看| 99国产精品免费福利视频| 美女福利国产在线|