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

    諧振面力作用下的無(wú)窮大雙周期加筋板的振動(dòng)響應(yīng)分析

    2015-04-26 08:07:47周海安王曉明宋鳳敏梅玉林
    船舶力學(xué) 2015年11期
    關(guān)鍵詞:加強(qiáng)筋薄板傅里葉

    周海安,王曉明,宋鳳敏,梅玉林

    (1山東理工大學(xué)機(jī)械工程學(xué)院,山東淄博255049;2大連理工大學(xué)a.機(jī)械工程學(xué)院;b.汽車工程學(xué)院,遼寧大連116024;3山東水利職業(yè)學(xué)院,山東日照276826)

    諧振面力作用下的無(wú)窮大雙周期加筋板的振動(dòng)響應(yīng)分析

    周海安1,王曉明2a,宋鳳敏3,梅玉林2b

    (1山東理工大學(xué)機(jī)械工程學(xué)院,山東淄博255049;2大連理工大學(xué)a.機(jī)械工程學(xué)院;b.汽車工程學(xué)院,遼寧大連116024;3山東水利職業(yè)學(xué)院,山東日照276826)

    文章主要研究了流體負(fù)載下的無(wú)窮大雙周期加筋薄板,在周期諧振力作用下的振動(dòng)響應(yīng)。假設(shè)薄板和加強(qiáng)筋只存在法向力的作用,通過(guò)薄板和加強(qiáng)筋的位移連續(xù)條件,以及格林函數(shù)及傅立葉變換公式的應(yīng)用,建立了頻域內(nèi)雙周期加筋薄板的振動(dòng)位移方程。利用空間波數(shù)法,將加筋薄板的振動(dòng)位移表達(dá)為波數(shù)分量的迭加形式。通過(guò)數(shù)值計(jì)算的方法對(duì)波數(shù)分量進(jìn)行求解,最后通過(guò)方程的傅立葉逆變換,得到了加筋薄板的振動(dòng)方程。數(shù)值計(jì)算的結(jié)果表明文中的計(jì)算方法收斂速度快。

    振動(dòng)響應(yīng);加筋板;格林函數(shù);空間波數(shù)法

    0 引言

    加筋板是許多工程結(jié)構(gòu)的組成部分,在工程界得到了廣泛的應(yīng)用,尤其在機(jī)械,船舶和航天領(lǐng)域,加筋板的應(yīng)用更為常見(jiàn)。到目前為止,許多學(xué)者對(duì)加筋板的振動(dòng)特性進(jìn)行過(guò)研究,并且提出了許多不同的求解方法[1-8]。然而大部分的科研工作主要集中在有限尺寸或單周期無(wú)窮大加筋板的振動(dòng)研究上,而對(duì)于水下無(wú)窮大雙周期加筋薄板振動(dòng)響應(yīng)的研究則要少得多。本文的主要目的是尋求一種合適的數(shù)學(xué)方法以便對(duì)雙周期加筋薄板結(jié)構(gòu)進(jìn)行振動(dòng)分析。但是由于加筋薄板的周期性以及流固耦合問(wèn)題的復(fù)雜性,幾種重要的數(shù)學(xué)方法如Rayleigh-Ritz能量法﹑有限差分法﹑傳遞矩陣法﹑模態(tài)法等并不能很好地解決無(wú)窮大周期結(jié)構(gòu)的振動(dòng)問(wèn)題。由Mead和Pujara[9]提出的空間波數(shù)法在處理流固耦合的周期結(jié)構(gòu)振動(dòng)問(wèn)題上體現(xiàn)出了明顯的優(yōu)勢(shì):Mace[1]在空間波數(shù)法的基礎(chǔ)上,推導(dǎo)出流體作用下雙周期加筋薄板的振動(dòng)和輻射聲壓方程,并且分析了點(diǎn)力作用下的水中遠(yuǎn)場(chǎng)輻射聲壓。Lee和Kim[10]等亦通過(guò)空間波數(shù)法對(duì)單周期加筋薄板的振動(dòng)響應(yīng)及聲透射系數(shù)進(jìn)行了理論分析與計(jì)算。此外,傅里葉變換的方法也是研究結(jié)構(gòu)振動(dòng)和聲輻射的一種有力工具。吳文偉等[11]曾經(jīng)采用傅里葉變換的方法對(duì)單周期加筋薄板的聲輻射進(jìn)行過(guò)研究。

    本文對(duì)水下無(wú)窮大雙周期加筋薄板的振動(dòng)響應(yīng)進(jìn)行了分析,并且根據(jù)格林函數(shù)與空間波數(shù)法以及傅里葉變換法,建立了加筋薄板的振動(dòng)方程,并最終通過(guò)數(shù)值計(jì)算的方法得到了方程的解。在計(jì)算結(jié)果中,本文重點(diǎn)分析了加強(qiáng)筋橫截面尺寸以及周期間距等幾個(gè)主要物理參數(shù)對(duì)薄板振動(dòng)響應(yīng)帶來(lái)的影響。

    1 加筋薄板的振動(dòng)方程

    1.1 數(shù)學(xué)模型

    無(wú)窮大雙周期加筋薄板的上表面空間充滿水流體,薄板的下表面空間是空氣。在x軸方向的加強(qiáng)筋沿x=mLx周期排列,在y軸方向的加強(qiáng)筋沿y=nLy(m,n=0,±1,…,±∞)周期排列,周期加筋薄板的局部結(jié)構(gòu)示意圖如圖1所示。

    圖1 雙周期加筋薄板Fig.1 Diagram of the two-dimensional beam-stiffened plate

    雙周期加筋薄板受到隨時(shí)間t周期變化的諧振面力f(x,y)exp(jωt)的作用,其方程表達(dá)式是[1]

    pe是諧振力的幅值,ω是圓頻率,k0是諧振力波數(shù),θ是諧振力與z軸方向的夾角,φ是諧振力在水平面上的分量與x軸方向的夾角。

    對(duì)雙周期加筋薄板進(jìn)行振動(dòng)響應(yīng)分析,假設(shè)加筋薄板的振動(dòng)位移是w(x,y)exp(jωt),加強(qiáng)筋與板的相互作用力僅考慮法向力的作用,則周期加筋板的振動(dòng)方程是(為了公式的簡(jiǎn)潔,在以下所有方程中都將省略時(shí)間因子exp(jωt))

    式中:G(x-ξ,y-η)是格林函數(shù),其物理含義是:薄板由于受到作用在(ξ,η)點(diǎn)的單位激勵(lì)力而在(x,y)位置處的振動(dòng)位移。F1(ξ,η)δ(ξ-mLx)和F2(ξ,η)δ(η-nLy)是加強(qiáng)筋沿直線ξ=mLx和η=nLy作用于薄板的法向力。這里假設(shè)兩組加強(qiáng)筋是完全相同的,利用薄板和加強(qiáng)筋的位移連續(xù)條件可以得到

    式中:E1,I,ρb,Sb分別是加強(qiáng)筋的楊氏模量,慣性矩,質(zhì)量密度和橫截面面積。

    1.2 振動(dòng)方程的傅里葉變換

    要求解方程(3),需要對(duì)方程進(jìn)行傅里葉變換。對(duì)任意滿足要求的函數(shù)g(x,y),其傅里葉變換是

    傅里葉逆變換是

    對(duì)于滿足傅里葉積分定理的任意函數(shù)g1(x,y)和g2(x,y),根據(jù)函數(shù)卷積的傅里葉變換等于函數(shù)的傅里葉變換的乘積的性質(zhì)可知

    將(4)式代入(3)式中,對(duì)(3)式進(jìn)行傅里葉變換,并根據(jù)(7)式的積分性質(zhì)可得

    式中:w*(α,nLy),w*(mLx,β)分別是w(x,nLy),w(mLx,y)關(guān)于x軸和y軸的單重傅里葉變換,并且

    流體負(fù)載下薄板的格林函數(shù)滿足方程

    式中:c是水流體中的聲速。

    薄板和水流體接觸面的位移邊界條件是

    對(duì)(12)式和(13)式進(jìn)行傅里葉變換,求解得到

    對(duì)(11)式進(jìn)行傅里葉變換,結(jié)合(14)式可以得到頻域內(nèi)流固耦合的薄板阻抗

    2 振動(dòng)方程的求解

    按照空間波數(shù)法理論[9]可知,二維周期加筋薄板的振動(dòng)位移可以表達(dá)為如下形式[1]

    式中:Wmn是(m,n)th諧振分量振幅,且kx,mx=kx+2mπ/Lx,ky,n=ky+2nπ/Ly。

    由(16)式可以得到振動(dòng)位移的周期條件

    根據(jù)(17)式加筋薄板的周期條件,可知

    將(18)式代入(8)式中,可得

    雙周期加筋薄板在原點(diǎn)處的振動(dòng)位移分別可以表達(dá)為

    根據(jù)(16)式和(20)式可知,加筋薄板沿x=0和y=0直線上的振動(dòng)位移分別是

    對(duì)(21a)式和(21b)式分別進(jìn)行關(guān)于y軸和x軸的單重傅立葉變換,得到:

    由泊松迭加公式可知

    同理可得

    將(9)、(22a)、(22b)式和(23a)、(23b)式代入(19)式中,可以得到頻域內(nèi)加筋薄板關(guān)于未知量Vn和Nm的振動(dòng)方程

    對(duì)(24)式進(jìn)行關(guān)于α的單重傅立葉逆變換,并且令變量x=0以及β=ky,n,根據(jù)(22a)式得到方程[1]

    同理,對(duì)(24)式進(jìn)行關(guān)于β的單重傅立葉逆變換,并令變量y=0以及α=kx,m,根據(jù)(22b)式得到方程?

    對(duì)(24)式進(jìn)行傅立葉逆變換,最終得到時(shí)域內(nèi)加筋薄板的振動(dòng)方程

    本文提出了無(wú)窮大雙周期加筋板振動(dòng)響應(yīng)的格林函數(shù)法,通過(guò)比較可以發(fā)現(xiàn),雖然本文與文獻(xiàn)[1]采用的方法不同,但是本文最終求解的振動(dòng)位移表達(dá)式(方程(28))和文獻(xiàn)[1]中的表達(dá)式(方程(34))和實(shí)際物理含義是相同的,從而證明了本文方法的正確性。相對(duì)于文獻(xiàn)[1]本文方法的優(yōu)點(diǎn)在于:格林函數(shù)的方法不僅可以用來(lái)計(jì)算諧振面力的作用,而且可以計(jì)算包括點(diǎn)力激勵(lì)在內(nèi)的任何作用力形式;通過(guò)格林函數(shù)的應(yīng)用,得到了頻域內(nèi)流固耦合的薄板阻抗,揭示了加強(qiáng)筋作用力與薄板振動(dòng)位移的周期關(guān)系((8)式);本文中將加強(qiáng)筋的作用力表達(dá)為薄板振動(dòng)位移的分量形式,其求解方法也與文獻(xiàn)[1]不同。

    3 計(jì)算分析

    下面將利用推導(dǎo)的公式來(lái)計(jì)算薄板和加筋薄板的振動(dòng),并且分析加強(qiáng)筋對(duì)薄板振動(dòng)的影響。計(jì)算中用到的部分物理參數(shù)如下:E=E1=2.01×1011×(1+0.02j),v=0.2,h=0.005,ρ=ρb=7 850,ρ0=1 000,c=1 500,pe=1 000,θ=φ=π/4,k0=ω/c。加強(qiáng)筋的橫截面面積是Sb=0.005×0.05,短邊與薄板相接。所有的單位均采用國(guó)際標(biāo)準(zhǔn)單位。本文在下面的算例分析中,通過(guò)計(jì)算發(fā)現(xiàn)當(dāng)Q=10時(shí),已經(jīng)得到滿足收斂性要求的計(jì)算精度。

    圖2是薄板在頻率為2 000 Hz諧振力作用下,角度分別為φ=π/4和φ=π/6,在-1≤x,y≤1且時(shí)間t=0的振動(dòng)響應(yīng)。

    圖2 薄板的振動(dòng)響應(yīng)Fig.2 The vibration responses of the unstiffened plate

    薄板的振動(dòng)特性主要取決于它的振動(dòng)相位角和幅值大小。薄板的振動(dòng)相位角變化方向由φ決定。從圖中可以看出,當(dāng)φ=π/4時(shí),振動(dòng)相位沿與x軸成45°的直線方向變化;如果改變?chǔ)盏拇笮?,如圖所示,振動(dòng)相位角變化方向?qū)?huì)隨之改變。與φ不同,θ并不改變相位角變化方向,但是可以改變薄板的振動(dòng)幅值大小。

    圖3是不同周期間距的加筋薄板在頻率為2 000 Hz諧振力作用下,在-1≤x,y≤1且時(shí)間t=0的振動(dòng)響應(yīng)。

    圖3 加筋薄板的振動(dòng)響應(yīng)Fig.3 The vibration responses of the beam-stiffened plate

    對(duì)比圖2和圖3可以發(fā)現(xiàn),相對(duì)于薄板結(jié)構(gòu),雙周期加筋薄板的振動(dòng)具有二維周期性質(zhì),而且加筋薄板的振幅變化比較明顯。從圖3(a)中可以看出,加筋薄板的振動(dòng)沿x軸和y軸方向上是完全一致的,這是由于沿兩個(gè)軸方向的周期間距相等而且壓力角φ=π/4。如果改變兩軸的周期間距或壓力角φ,如圖3(b)所示,加筋薄板沿x軸和y軸方向的振動(dòng)將不再具有對(duì)稱性。

    圖4給出了薄板和不同周期間距加筋薄板在諧振力頻率為50~2 500 Hz范圍內(nèi)分別在(0,0)點(diǎn)和(0.4,0.4)點(diǎn)處的振幅。周期間距分別是Lx=Ly=0.2和Lx=Ly=0.4。

    從圖4中可以發(fā)現(xiàn),薄板在(0,0)點(diǎn)的振幅幾乎為0,且隨頻率的增加不發(fā)生改變,表明在薄板中傳播的振動(dòng)彎曲波在通過(guò)原點(diǎn)位置時(shí)受到阻抑作用;在(0.4,0.4)點(diǎn)薄板的振幅隨頻率的增加而變化,振動(dòng)波的傳播可以順利通過(guò)。加強(qiáng)筋對(duì)薄板振動(dòng)的影響在100~1 300 Hz赫茲范圍內(nèi)作用明顯,相對(duì)于薄板,加筋薄板的振幅振蕩較大。隨著頻率的進(jìn)一步增加,從整體分析,加強(qiáng)筋的作用減弱,加筋薄板振幅的振蕩峰值減小,振幅的變化隨頻率的增加趨于平緩。當(dāng)周期間距增大后,加筋薄板在更多的頻率位置出現(xiàn)共振現(xiàn)象,表明適當(dāng)增大周期間距可以增強(qiáng)加筋薄板的振蕩效果,從而在較多的頻率范圍內(nèi)影響薄板的振動(dòng)。而且由于加強(qiáng)筋周期間距的增大,使得加筋薄板的共振頻率向低頻范圍偏移,因此可以通過(guò)調(diào)節(jié)周期間距的方法改變加筋薄板的共振頻率。

    圖4 在不同位置點(diǎn)的振動(dòng)響應(yīng)Fig.4 The vibration responses at different points

    圖5 不同加強(qiáng)筋對(duì)振幅的影響Fig.5 The influence of different beam-stiffeners on the vibration amplitudes of the stiffened plate

    為了分析不同尺寸的加強(qiáng)筋對(duì)薄板的影響,圖5給出了加強(qiáng)筋橫截面面積分別為0.01×0.1和0.01× 0.15的周期加筋薄板在(0.4,0.4)點(diǎn)隨頻率變化的振幅曲線。

    對(duì)比圖5(a)和5(b)發(fā)現(xiàn),增加加強(qiáng)筋的橫截面高度對(duì)加筋薄板在200 Hz頻率處的振幅影響比較明顯,但對(duì)其它頻率振幅的影響相對(duì)較小,而且振幅曲線在400 Hz之后隨頻率的增加其變化趨勢(shì)是相似的。當(dāng)同時(shí)改變加強(qiáng)筋橫截面的寬度和高度尺寸后,通過(guò)對(duì)比圖4(b)和5(a)可以發(fā)現(xiàn),周期間距為L(zhǎng)x=Ly=0.4的加筋薄板的振幅曲線發(fā)生變化。從總體上分析,加強(qiáng)筋橫截面尺寸的變化對(duì)1 000 Hz以內(nèi)振動(dòng)響應(yīng)的作用較為顯著,而對(duì)較高頻率振動(dòng)響應(yīng)的影響較弱。

    4 結(jié)論

    根據(jù)格林函數(shù)和空間波數(shù)法以及傅立葉變換的應(yīng)用,本文提出了一種求解受諧振力作用的流體加載下的無(wú)窮大雙周期加筋薄板振動(dòng)方程的方法。通過(guò)將加筋薄板的振動(dòng)位移展開(kāi)為波數(shù)分量迭加的形式,并采用數(shù)值計(jì)算的方法求解分量的值,從而得到了加筋薄板的振動(dòng)方程表達(dá)式。在計(jì)算結(jié)果中,本文詳細(xì)分析了加強(qiáng)筋對(duì)薄板振動(dòng)響應(yīng)的作用,結(jié)果表明周期加強(qiáng)筋對(duì)薄板的低頻振動(dòng)影響較大。

    參考文獻(xiàn):

    [1]Mace B R.Sound radiation from fluid loaded orthogonally stiffened plates[J].Journal of Sound and Vibration,1981,71(3): 435-441.

    [2]Alain Berry.A new formulation for the vibrations and sound radiation of fluid-loaded plates with elastic boundary conditions[J].Journal of the Acoustical Society of America,1994,96(2):889-901.

    [3]洪明,孫盈盈.船上含敷料板與加筋結(jié)構(gòu)的固有振動(dòng)特性研究[J].船舶力學(xué),2008,12(2):311-318. Hong Ming,Sun Yingying.Natural vibration characteristics of stiffened plates and plates with covering layer in ships[J]. Journal of Ship Mechanics,2008,12(2):311-318.

    [4]Lin T R,Pan J.A closed form solution for the dynamic response of finite ribbed plates[J].Journal of the Acoustical Society of America,2006,119(2):917-925.

    [5]黎勝,楊婧媛.水下加筋板振動(dòng)聲輻射的代理模型研究[J].聲學(xué)學(xué)報(bào),2010,35(6):659-664. Li Sheng,Yang Jingyuan.Research on surrogate models for structural vibration and acoustic radiation of underwater stiffened plates[J].Acta Acustica,2010,35(6):659-664.

    [6]Maxit L.Wavenumber space and physical space responses of a periodically ribbed plate to a point drive:A discrete approach[J].Applied Acoustics,2009,70(4):563-578.

    [7]Mace B R.Sound radiation from a plate reinforced by two sets of parallel stiffeners[J].Journal of Sound and Vibration, 1980,71(3):435-441.

    [8]陳海龍,金葉青,龐福振.一種分析單向加筋板隔聲性能的解析—數(shù)值耦合法[J].船舶力學(xué),2013,17(8):952-959. Chen Hailong,Jin Yeqing,pang Fuzhen.Sound insulation performance analysis of single stiffened plate by a combined analytical-numerical method[J].Journal of Ship Mechanics,2013,17(8):952-959.

    [9]Mead D J,Pujara K K.Space harmonic analysis of periodically supported beams:Response to convected random loading [J].Journal of Sound and Vibration,1971,14(2):525-532.

    [10]Lee J H,Kim J.Analysis of sound transmission through periodically stiffened panels by space-harmonic expansion method [J].Journal of Sound and Vibration,2002,251(3):349-366.

    [11]吳文偉,冷文浩,沈順根.具有等間距相同加強(qiáng)筋板的聲輻射[J].中國(guó)造船,1999,146(3):72-81. Wu Wenwei,Leng Wenhao,Shen Shungen.Sound radiation of plate with equidistant reinforced ribs[J].Shipbuilding of China,1999,146(3):72-81.

    Investigation of the vibration response of an infinite two-dimensional rib-stiffened plate under a harmonic plane pressure

    ZHOU Hai-an1,WANG Xiao-ming2a,SONG Feng-min3,MEI Yu-lin2b
    (1.School of Mechanical Engineering,Shandong University of Technology,Zibo 255049,China;2a.School of Mechanical Engineering;b.School of Automotive Engineering,Dalian University of Technology,Dalian 116024,China;3.Shandong Water Polytechnic,Rizhao 276826,China)

    The vibration response of an infinite fluid-loaded plate,reinforced with two sets of orthogonal rib stiffeners,is investigated.The stiffeners are assumed to exert only forces on the plate.With the displacement continuity condition of the plate and the stiffeners,the vibration response of the rib-stiffened plate to a harmonic pressure in the wavenumber space is established by using the Green function and Fourier transforms.Applying the space harmonic method,the vibration displacement of the stiffened plate is then expressed as infinite sets of space harmonic amplitudes.After having a numerical solution for these amplitudes,the Fourier inverse transform is employed to find an expression for the vibration response of the stiffened plate in the physical space.Numerical results show that this method converges rapidly.

    vibration response;rib-stiffened plate;Green function;space harmonic method

    U661.4

    A

    10.3969/j.issn.1007-7294.2015.11.008

    1007-7294(2015)11-1352-08

    2015-07-15

    國(guó)家自然科學(xué)基金資助(51405276,50875030)

    周海安(1981-),男,博士,講師,通訊作者,E-mail:zhouhaian1981@163.com;

    王曉明(1965-),男,教授,博士生導(dǎo)師,E-mail:xiaoming@dlut.edu.cn;

    宋鳳敏(1978-),女,講師,E-mail:stj08@163.com;

    梅玉林(1969-),女,副教授,E-mail:meiyulin@dlut.edu.cn。

    猜你喜歡
    加強(qiáng)筋薄板傅里葉
    一角點(diǎn)支撐另一對(duì)邊固支正交各向異性矩形薄板彎曲的辛疊加解
    10MN鋁合金薄板拉伸機(jī)組的研制
    汽車零件加強(qiáng)筋結(jié)構(gòu)設(shè)計(jì)
    汽車工程師(2019年7期)2019-08-12 07:25:04
    越野車側(cè)壁加強(qiáng)筋板工藝改進(jìn)
    模具制造(2019年3期)2019-06-06 02:10:52
    雙線性傅里葉乘子算子的量化加權(quán)估計(jì)
    基于小波降噪的稀疏傅里葉變換時(shí)延估計(jì)
    加強(qiáng)筋的截面尺寸對(duì)注塑制品翹曲變形的影響
    塑料制造(2016年5期)2016-06-15 20:27:39
    基于傅里葉變換的快速TAMVDR算法
    鋁薄板高速DP-GMAW焊接性能的研究
    焊接(2016年5期)2016-02-27 13:04:42
    基于PLC的薄板激光焊機(jī)控制系統(tǒng)
    焊接(2015年10期)2015-07-18 11:04:46
    少妇的逼好多水| 国产亚洲av片在线观看秒播厂| 久久av网站| 欧美xxⅹ黑人| 亚洲一级一片aⅴ在线观看| 久久国产亚洲av麻豆专区| 建设人人有责人人尽责人人享有的| 日本欧美视频一区| 久久国产精品大桥未久av | 成年av动漫网址| 午夜视频国产福利| 午夜福利在线观看免费完整高清在| 国产精品秋霞免费鲁丝片| 日韩强制内射视频| 亚洲va在线va天堂va国产| 国产69精品久久久久777片| 美女中出高潮动态图| 最近中文字幕2019免费版| av福利片在线| 丰满迷人的少妇在线观看| 下体分泌物呈黄色| a级一级毛片免费在线观看| 日本-黄色视频高清免费观看| 日本wwww免费看| 一区二区三区四区激情视频| 国产精品蜜桃在线观看| 精品酒店卫生间| 亚洲精品日韩av片在线观看| 久久 成人 亚洲| 国产精品久久久久久精品古装| 久久久欧美国产精品| 午夜福利网站1000一区二区三区| 精品久久久久久久久av| 五月天丁香电影| 十八禁高潮呻吟视频 | 亚洲欧美清纯卡通| 国产日韩欧美在线精品| 欧美三级亚洲精品| 少妇的逼好多水| 亚洲精品日本国产第一区| 免费黄频网站在线观看国产| 国产成人精品久久久久久| 久久久久久伊人网av| 国产成人a∨麻豆精品| 国产精品久久久久久精品古装| 国产毛片在线视频| 69精品国产乱码久久久| 99久久综合免费| 中文字幕亚洲精品专区| 亚洲经典国产精华液单| 一本—道久久a久久精品蜜桃钙片| 亚洲婷婷狠狠爱综合网| 久久久a久久爽久久v久久| 黑丝袜美女国产一区| 人人妻人人澡人人看| 高清不卡的av网站| 亚洲国产精品国产精品| 黑丝袜美女国产一区| 少妇人妻久久综合中文| 国内精品宾馆在线| 成人午夜精彩视频在线观看| 只有这里有精品99| 91精品伊人久久大香线蕉| 久久久久久久国产电影| 黑丝袜美女国产一区| 亚洲av男天堂| 国产淫片久久久久久久久| 在线观看一区二区三区激情| 国内少妇人妻偷人精品xxx网站| 又粗又硬又长又爽又黄的视频| 2021少妇久久久久久久久久久| 久久久a久久爽久久v久久| 日本免费在线观看一区| 九色成人免费人妻av| 男人爽女人下面视频在线观看| 夜夜爽夜夜爽视频| 国产精品一区二区性色av| 在线观看三级黄色| 男人和女人高潮做爰伦理| 亚洲av成人精品一区久久| 在线看a的网站| 国产男人的电影天堂91| 亚洲精品日本国产第一区| 亚洲丝袜综合中文字幕| 汤姆久久久久久久影院中文字幕| 免费观看性生交大片5| 特大巨黑吊av在线直播| 熟女电影av网| 亚洲欧洲日产国产| 在线观看www视频免费| 午夜福利,免费看| 国产精品成人在线| 男男h啪啪无遮挡| 老熟女久久久| 边亲边吃奶的免费视频| a级毛片免费高清观看在线播放| 亚洲欧美一区二区三区国产| 国产在线男女| av免费观看日本| 最近中文字幕高清免费大全6| 晚上一个人看的免费电影| 久久鲁丝午夜福利片| 最新中文字幕久久久久| av不卡在线播放| 欧美97在线视频| 国产女主播在线喷水免费视频网站| 一边亲一边摸免费视频| 免费看av在线观看网站| 日韩熟女老妇一区二区性免费视频| 精品国产一区二区久久| 亚洲国产精品成人久久小说| 欧美高清成人免费视频www| 亚洲真实伦在线观看| 久久影院123| 韩国高清视频一区二区三区| 十八禁网站网址无遮挡 | 久久久久久久亚洲中文字幕| freevideosex欧美| 在线观看免费高清a一片| 久久久久久久久久久丰满| 国产av码专区亚洲av| 国产亚洲欧美精品永久| 欧美丝袜亚洲另类| 国产高清不卡午夜福利| √禁漫天堂资源中文www| 麻豆乱淫一区二区| 成人亚洲精品一区在线观看| 久久久久久久久大av| 国产老妇伦熟女老妇高清| a级一级毛片免费在线观看| 亚洲欧美精品专区久久| 99九九在线精品视频 | 亚洲国产精品一区三区| 亚洲不卡免费看| 插逼视频在线观看| 美女xxoo啪啪120秒动态图| 国产爽快片一区二区三区| 亚洲成人av在线免费| 大片电影免费在线观看免费| 成人二区视频| 中文字幕人妻丝袜制服| 久久午夜综合久久蜜桃| 丁香六月天网| 2022亚洲国产成人精品| 国产免费一区二区三区四区乱码| 成人午夜精彩视频在线观看| 菩萨蛮人人尽说江南好唐韦庄| 色婷婷久久久亚洲欧美| 亚洲久久久国产精品| 国产欧美亚洲国产| 各种免费的搞黄视频| 久久久久国产精品人妻一区二区| 国产精品99久久久久久久久| 少妇被粗大猛烈的视频| 夜夜骑夜夜射夜夜干| 一级毛片aaaaaa免费看小| 久久久午夜欧美精品| 久久久久久久精品精品| 免费人成在线观看视频色| av播播在线观看一区| 免费观看无遮挡的男女| av视频免费观看在线观看| 深夜a级毛片| 美女内射精品一级片tv| 人人妻人人看人人澡| 日韩av在线免费看完整版不卡| 亚洲情色 制服丝袜| av线在线观看网站| 在线观看av片永久免费下载| 日日爽夜夜爽网站| 国产欧美日韩精品一区二区| 菩萨蛮人人尽说江南好唐韦庄| 在线 av 中文字幕| 国产在线免费精品| 男女边吃奶边做爰视频| 又爽又黄a免费视频| 欧美精品人与动牲交sv欧美| 国产av国产精品国产| 国产亚洲av片在线观看秒播厂| 91午夜精品亚洲一区二区三区| 日韩成人av中文字幕在线观看| 久久精品国产a三级三级三级| 人妻夜夜爽99麻豆av| 亚洲av电影在线观看一区二区三区| 亚洲国产成人一精品久久久| 我的老师免费观看完整版| 成人18禁高潮啪啪吃奶动态图 | 久久久国产欧美日韩av| 亚洲真实伦在线观看| 亚洲婷婷狠狠爱综合网| 亚洲,欧美,日韩| 极品人妻少妇av视频| 曰老女人黄片| 国产黄片美女视频| 夫妻午夜视频| 极品少妇高潮喷水抽搐| 日本色播在线视频| 黄色视频在线播放观看不卡| av福利片在线观看| 18禁在线播放成人免费| av黄色大香蕉| 欧美日韩在线观看h| 能在线免费看毛片的网站| h视频一区二区三区| 欧美日韩视频精品一区| 最近的中文字幕免费完整| 免费大片黄手机在线观看| 久久人人爽人人爽人人片va| 黑人巨大精品欧美一区二区蜜桃 | 熟女电影av网| 日本av手机在线免费观看| 国产亚洲av片在线观看秒播厂| 欧美亚洲 丝袜 人妻 在线| 精品少妇内射三级| 久久久精品免费免费高清| 少妇丰满av| 国产老妇伦熟女老妇高清| 亚洲国产日韩一区二区| 久久午夜综合久久蜜桃| 大又大粗又爽又黄少妇毛片口| 如日韩欧美国产精品一区二区三区 | 多毛熟女@视频| 中文字幕免费在线视频6| 久久久久久久久久久久大奶| 最近手机中文字幕大全| 亚洲,欧美,日韩| 国产精品人妻久久久影院| 天堂俺去俺来也www色官网| 你懂的网址亚洲精品在线观看| 亚洲情色 制服丝袜| 国产成人精品一,二区| 中文欧美无线码| 美女cb高潮喷水在线观看| 久久av网站| 五月玫瑰六月丁香| 十八禁高潮呻吟视频 | 国产在线一区二区三区精| 91午夜精品亚洲一区二区三区| 18禁在线无遮挡免费观看视频| a级一级毛片免费在线观看| 菩萨蛮人人尽说江南好唐韦庄| a级毛色黄片| 国产熟女欧美一区二区| 美女xxoo啪啪120秒动态图| 成人美女网站在线观看视频| 美女国产视频在线观看| 久久久久久久国产电影| 一区二区三区精品91| 亚洲精品国产av成人精品| 亚洲精华国产精华液的使用体验| 久久久久久久久大av| 热re99久久精品国产66热6| 免费观看在线日韩| 男男h啪啪无遮挡| kizo精华| 一级毛片电影观看| 精品亚洲成国产av| 精品国产一区二区久久| 国产精品蜜桃在线观看| 国产黄频视频在线观看| 国产 精品1| 三级经典国产精品| 亚洲精品国产色婷婷电影| 极品少妇高潮喷水抽搐| 亚洲欧洲国产日韩| 天堂俺去俺来也www色官网| 亚洲av日韩在线播放| 亚洲欧美日韩卡通动漫| 中国美白少妇内射xxxbb| 免费不卡的大黄色大毛片视频在线观看| 欧美日韩亚洲高清精品| 免费人妻精品一区二区三区视频| 亚洲国产色片| 男女免费视频国产| 一区二区三区乱码不卡18| 亚洲性久久影院| 自拍欧美九色日韩亚洲蝌蚪91 | 大码成人一级视频| 老女人水多毛片| 国产一区二区在线观看日韩| 国产美女午夜福利| 亚洲国产精品一区二区三区在线| 99九九在线精品视频 | 中文字幕精品免费在线观看视频 | 亚洲欧美成人精品一区二区| 各种免费的搞黄视频| 免费黄网站久久成人精品| 国产黄色免费在线视频| 久久精品夜色国产| 边亲边吃奶的免费视频| 女性被躁到高潮视频| 国产男女超爽视频在线观看| 亚洲欧美日韩卡通动漫| 欧美日本中文国产一区发布| 国产高清国产精品国产三级| 一级爰片在线观看| 日韩制服骚丝袜av| 纯流量卡能插随身wifi吗| 免费人妻精品一区二区三区视频| 免费黄色在线免费观看| 中文天堂在线官网| 99九九在线精品视频 | 男人和女人高潮做爰伦理| 欧美最新免费一区二区三区| 啦啦啦视频在线资源免费观看| 国产成人午夜福利电影在线观看| 日韩av免费高清视频| 又爽又黄a免费视频| 国产一区二区在线观看日韩| 特大巨黑吊av在线直播| av播播在线观看一区| 五月天丁香电影| 人体艺术视频欧美日本| 久久 成人 亚洲| 亚洲久久久国产精品| 久久ye,这里只有精品| 丝袜喷水一区| 亚洲伊人久久精品综合| 99热这里只有是精品在线观看| 中文天堂在线官网| 亚洲欧美中文字幕日韩二区| 亚洲国产精品一区二区三区在线| 国产精品嫩草影院av在线观看| 美女脱内裤让男人舔精品视频| 在线观看免费高清a一片| 免费看光身美女| 欧美另类一区| 国产欧美日韩精品一区二区| 多毛熟女@视频| 国产色婷婷99| 亚洲自偷自拍三级| 国产色婷婷99| 精品久久久久久久久av| 午夜免费观看性视频| 各种免费的搞黄视频| av黄色大香蕉| 日韩大片免费观看网站| 国产深夜福利视频在线观看| 亚洲欧洲国产日韩| 亚洲人成网站在线播| 亚洲美女黄色视频免费看| 国产乱来视频区| 22中文网久久字幕| 熟女人妻精品中文字幕| 伊人亚洲综合成人网| 国产淫片久久久久久久久| 亚洲欧美精品自产自拍| 99久国产av精品国产电影| 黑丝袜美女国产一区| 热re99久久国产66热| 国产在线视频一区二区| 一二三四中文在线观看免费高清| 人妻制服诱惑在线中文字幕| 卡戴珊不雅视频在线播放| 黄色配什么色好看| 国产精品久久久久久av不卡| 夫妻性生交免费视频一级片| 一本久久精品| 亚洲电影在线观看av| 一区二区三区免费毛片| 黄色日韩在线| 亚洲精品国产av成人精品| 99视频精品全部免费 在线| av国产久精品久网站免费入址| 国产成人精品福利久久| 黄色怎么调成土黄色| 99国产精品免费福利视频| 在线观看国产h片| a级片在线免费高清观看视频| 欧美少妇被猛烈插入视频| 免费人成在线观看视频色| 九色成人免费人妻av| 日韩精品免费视频一区二区三区 | 日韩欧美 国产精品| 男人爽女人下面视频在线观看| 国产免费一级a男人的天堂| 国产亚洲精品久久久com| 成年美女黄网站色视频大全免费 | 成人综合一区亚洲| 欧美国产精品一级二级三级 | 久久精品久久精品一区二区三区| 欧美日本中文国产一区发布| 久久99热这里只频精品6学生| 麻豆成人av视频| 亚洲欧洲精品一区二区精品久久久 | 三级国产精品片| 这个男人来自地球电影免费观看 | 午夜福利,免费看| 久久综合国产亚洲精品| 中文字幕久久专区| 精品国产乱码久久久久久小说| 国精品久久久久久国模美| 国产男人的电影天堂91| 好男人视频免费观看在线| 亚洲第一av免费看| 国产极品天堂在线| 国产综合精华液| 毛片一级片免费看久久久久| 久久精品国产亚洲网站| 91久久精品电影网| 亚洲高清免费不卡视频| 在线精品无人区一区二区三| 久久人人爽人人爽人人片va| 高清毛片免费看| 亚洲精品自拍成人| 2021少妇久久久久久久久久久| 成人亚洲精品一区在线观看| 久久人人爽av亚洲精品天堂| 精品久久久噜噜| 久久久午夜欧美精品| av不卡在线播放| 欧美日韩综合久久久久久| 91久久精品国产一区二区三区| 黑人巨大精品欧美一区二区蜜桃 | 少妇人妻 视频| 国产亚洲精品久久久com| 80岁老熟妇乱子伦牲交| 26uuu在线亚洲综合色| 夜夜骑夜夜射夜夜干| 另类亚洲欧美激情| 久久久久视频综合| 免费看日本二区| 日韩av在线免费看完整版不卡| 亚洲高清免费不卡视频| 午夜免费男女啪啪视频观看| 亚洲av.av天堂| 欧美丝袜亚洲另类| 日韩欧美 国产精品| 91精品伊人久久大香线蕉| 午夜福利,免费看| 亚洲av福利一区| 又爽又黄a免费视频| 亚洲精品日韩av片在线观看| 国产黄色视频一区二区在线观看| 久久人妻熟女aⅴ| 国产视频首页在线观看| 91久久精品国产一区二区三区| 中文字幕制服av| 日本av免费视频播放| 老司机影院成人| 国产精品一区二区三区四区免费观看| av一本久久久久| 美女中出高潮动态图| 99热网站在线观看| 亚洲精品日本国产第一区| .国产精品久久| 色哟哟·www| 国产成人精品一,二区| 亚洲伊人久久精品综合| 欧美日韩亚洲高清精品| 一区二区三区精品91| 人妻少妇偷人精品九色| 插逼视频在线观看| 少妇猛男粗大的猛烈进出视频| 久久久久久久久久成人| 国产成人精品婷婷| 国产色爽女视频免费观看| 精品国产一区二区久久| 婷婷色av中文字幕| 永久免费av网站大全| 免费观看性生交大片5| 国产精品国产三级专区第一集| 99久国产av精品国产电影| 男的添女的下面高潮视频| 在线观看www视频免费| 国内揄拍国产精品人妻在线| 狠狠精品人妻久久久久久综合| 中国三级夫妇交换| tube8黄色片| 97在线视频观看| 内射极品少妇av片p| 久久国产亚洲av麻豆专区| 欧美日韩一区二区视频在线观看视频在线| 丝袜脚勾引网站| 久久99精品国语久久久| 久久精品国产a三级三级三级| 777米奇影视久久| 亚洲人成网站在线观看播放| 国产成人免费无遮挡视频| 精品国产露脸久久av麻豆| 亚洲欧美一区二区三区国产| 国产一区二区在线观看日韩| 精品久久久噜噜| 女人精品久久久久毛片| 久久久欧美国产精品| 噜噜噜噜噜久久久久久91| 少妇高潮的动态图| 2021少妇久久久久久久久久久| 久久久午夜欧美精品| 国产爽快片一区二区三区| 在线观看三级黄色| 亚洲,欧美,日韩| 成人无遮挡网站| 国产永久视频网站| 久久久久久久国产电影| 2018国产大陆天天弄谢| 日本黄色日本黄色录像| 国产在视频线精品| 亚洲av欧美aⅴ国产| 免费在线观看成人毛片| 精品久久久久久久久亚洲| 涩涩av久久男人的天堂| 欧美另类一区| 看免费成人av毛片| 亚洲精品日本国产第一区| 国产精品无大码| 中文在线观看免费www的网站| 有码 亚洲区| 久久久久国产网址| 亚洲精品日本国产第一区| 国产真实伦视频高清在线观看| videos熟女内射| av天堂久久9| 国产男女超爽视频在线观看| 亚洲欧美日韩另类电影网站| 国产欧美日韩综合在线一区二区 | 久久久久国产网址| 欧美日韩国产mv在线观看视频| 日韩欧美 国产精品| 国产精品久久久久久久久免| 日本wwww免费看| 成人午夜精彩视频在线观看| 毛片一级片免费看久久久久| 日韩av不卡免费在线播放| 国产爽快片一区二区三区| 九色成人免费人妻av| 成年人免费黄色播放视频 | .国产精品久久| 国产美女午夜福利| 亚洲av成人精品一二三区| 日韩欧美一区视频在线观看 | 久久国产精品大桥未久av | 亚洲av成人精品一二三区| 日本91视频免费播放| 一级片'在线观看视频| 看免费成人av毛片| 国产精品嫩草影院av在线观看| 亚洲真实伦在线观看| videos熟女内射| 一级二级三级毛片免费看| 熟女电影av网| 婷婷色av中文字幕| 国产男人的电影天堂91| 成年人免费黄色播放视频 | 最近中文字幕2019免费版| 久久精品国产亚洲网站| 色吧在线观看| 亚洲在久久综合| 老司机影院成人| 最黄视频免费看| 欧美97在线视频| 日日摸夜夜添夜夜爱| 久久久久久久久久久久大奶| 久久精品国产亚洲av天美| 国产国拍精品亚洲av在线观看| 99国产精品免费福利视频| 在线 av 中文字幕| 国产极品粉嫩免费观看在线 | 日韩精品免费视频一区二区三区 | 99久久精品热视频| 日韩在线高清观看一区二区三区| 性色avwww在线观看| av播播在线观看一区| 国产日韩欧美在线精品| 久久6这里有精品| 久久免费观看电影| 九九爱精品视频在线观看| 九色成人免费人妻av| 国产黄片美女视频| 日韩一区二区三区影片| 亚洲精品色激情综合| 99国产精品免费福利视频| 老司机影院成人| 久久久久精品久久久久真实原创| 午夜免费观看性视频| 一区二区三区免费毛片| 在线观看免费日韩欧美大片 | 久久久国产精品麻豆| 欧美xxxx性猛交bbbb| 亚洲怡红院男人天堂| 亚洲无线观看免费| 日韩亚洲欧美综合| 精华霜和精华液先用哪个| 黄色欧美视频在线观看| 男男h啪啪无遮挡| 免费高清在线观看视频在线观看| 欧美精品高潮呻吟av久久| 亚洲欧美日韩卡通动漫| 男人舔奶头视频| 日本91视频免费播放| 久久精品久久精品一区二区三区| 久久热精品热| 精品亚洲成国产av| 最近中文字幕高清免费大全6| 22中文网久久字幕| videos熟女内射| 啦啦啦中文免费视频观看日本| 男女啪啪激烈高潮av片| 国产亚洲av片在线观看秒播厂| 久久99热这里只频精品6学生| 日韩人妻高清精品专区| 91成人精品电影| 亚洲,一卡二卡三卡| 在线观看av片永久免费下载| 精品亚洲成a人片在线观看| 日本爱情动作片www.在线观看| 天天操日日干夜夜撸| 又爽又黄a免费视频| 观看美女的网站| 永久网站在线| 三上悠亚av全集在线观看 | .国产精品久久| 少妇 在线观看| 日韩大片免费观看网站| 国产在线免费精品|