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

    環(huán)狀粘彈層合懸臂薄壁圓柱殼的振動(dòng)特性研究

    2017-11-28 07:00:12翟敬宇
    船舶力學(xué) 2017年11期
    關(guān)鍵詞:模態(tài)振動(dòng)

    王 宇,翟敬宇,李 昌,宋 華

    (1.遼寧科技大學(xué) 機(jī)械工程學(xué)院,遼寧 鞍山 114051;2.大連理工大學(xué) 機(jī)械工程學(xué)院,遼寧 大連 116024)

    環(huán)狀粘彈層合懸臂薄壁圓柱殼的振動(dòng)特性研究

    王 宇1,翟敬宇2,李 昌1,宋 華1

    (1.遼寧科技大學(xué) 機(jī)械工程學(xué)院,遼寧 鞍山 114051;2.大連理工大學(xué) 機(jī)械工程學(xué)院,遼寧 大連 116024)

    文章采用傳遞矩陣法分析了懸臂邊界條件下環(huán)狀局部覆蓋粘彈層合薄壁圓柱殼的振動(dòng)特性?;跇?lè)甫薄殼理論,結(jié)合粘彈性阻尼的變形協(xié)調(diào)關(guān)系和層間作用力關(guān)系,建立了基層和約束層薄壁圓柱殼的一階狀態(tài)微分方程,利用傳遞矩陣法推導(dǎo)了結(jié)構(gòu)的整體傳遞矩陣,并通過(guò)高精度的精細(xì)積分方法進(jìn)行求解,得到了固有頻率、損耗因子和三維振型,探討了約束阻尼層位置變化時(shí)對(duì)振動(dòng)特性的影響,并通過(guò)有限元法進(jìn)行了比較,通過(guò)算例驗(yàn)證了傳遞矩陣法對(duì)模態(tài)特性分析的有效性,前25階模態(tài)以周向振動(dòng)為主,最低階固有頻率對(duì)應(yīng)的三維模態(tài)振型為(1,5),并且在懸臂端的振動(dòng)位移最大,約束阻尼層覆蓋位置對(duì)薄壁圓柱殼振動(dòng)特性的影響較大。

    粘彈層合薄壁圓柱殼;懸臂邊界條件;環(huán)狀覆蓋;三維振型

    0 引 言

    薄壁圓柱殼構(gòu)件在船舶、潛艇和航空航天等機(jī)械領(lǐng)域中應(yīng)用廣泛,在外部復(fù)雜載荷等實(shí)際工況條件下,容易發(fā)生共振、失穩(wěn)和疲勞損傷等故障,在服役過(guò)程通常使用敷設(shè)阻尼層對(duì)其進(jìn)行減振降噪處理,由于實(shí)際管系布置和結(jié)構(gòu)優(yōu)化等因素,有時(shí)采用局部覆蓋約束阻尼層的方式,近年來(lái)其振動(dòng)控制一直是研究的熱點(diǎn)問(wèn)題之一[1-2]。對(duì)于懸臂邊界條件(即固支—自由約束條件)的薄壁圓柱殼構(gòu)件,在基頻激勵(lì)作用下,作者對(duì)發(fā)生共振時(shí)的振動(dòng)響應(yīng)特征進(jìn)行了研究,說(shuō)明了對(duì)殼體結(jié)構(gòu)進(jìn)行模態(tài)振型、模態(tài)數(shù)量和振動(dòng)特征分析的關(guān)聯(lián)和重要性[3-4]。然而,在懸臂邊界條件下,對(duì)帶有局部覆蓋粘彈層合薄壁圓柱殼的高階振動(dòng)問(wèn)題,例如約束阻尼層位于殼體不同位置時(shí)的三維模態(tài)振型等方面,尚缺少相應(yīng)的研究成果,因此,對(duì)其振動(dòng)特性進(jìn)行深入研究,具有重要的理論意義和工程應(yīng)用價(jià)值[5-8]。

    利用粘彈性材料的高阻尼性能對(duì)薄壁圓柱殼構(gòu)件進(jìn)行振動(dòng)抑制處理,改善薄壁圓柱殼構(gòu)件的動(dòng)力學(xué)性能,很多學(xué)者進(jìn)行了這一領(lǐng)域的研究,并取得了眾多成果,具有代表性的文獻(xiàn)如下:Wang和Chen[9]利用模態(tài)假設(shè)法研究了局部覆蓋約束阻尼層圓柱殼的動(dòng)力學(xué)問(wèn)題,結(jié)果表明結(jié)構(gòu)阻尼隨著約束層厚度和剛度的增大而增加,當(dāng)約束層厚度達(dá)到一定值時(shí)阻尼效果并不明顯。Chen和Huang[10-11]研究了固定大小的帶狀被動(dòng)約束阻尼層個(gè)數(shù)對(duì)圓柱殼振動(dòng)衰減的影響。Sainsbury[12]和Mastib[13]采用有限元法研究了部分覆蓋約束阻尼層圓柱殼應(yīng)變能強(qiáng)度的分布規(guī)律和減振方法。Mohammadi等[14]根據(jù)粘彈性層厚度方向上的位移變化關(guān)系,推導(dǎo)了約束阻尼層合殼的有限元模型,研究得到了其振動(dòng)特性。李恩奇等[15]基于唐奈爾薄殼理論和哈密頓原理,推導(dǎo)出局部覆蓋約束層阻尼圓柱殼的運(yùn)動(dòng)方程,并通過(guò)Laplace變換和傳遞函數(shù)法得到了對(duì)邊固支邊界條件下的固有頻率、損耗因子和頻響曲線(xiàn)。楊少紅和王安穩(wěn)[16]應(yīng)用混合分層理論,推導(dǎo)出兩端簡(jiǎn)支粘彈性層合圓柱殼的動(dòng)力學(xué)方程,得到了結(jié)構(gòu)的振動(dòng)頻率和損耗因子。王祖華等[17]利用模態(tài)展開(kāi)和Green函數(shù)法,分析了兩端簡(jiǎn)支邊界條件下敷設(shè)阻尼層環(huán)肋圓柱殼在流場(chǎng)中的聲輻射特性。曹雄濤等[18]利用桑德斯薄殼理論和波傳播法,求解了兩端簡(jiǎn)支被動(dòng)約束層阻尼圓柱殼的自由振動(dòng)特性。李超等[19]建立了圓柱殼阻尼材料的優(yōu)化模型,對(duì)約束阻尼和自由阻尼材料布局進(jìn)行了優(yōu)化研究。王淼等[20]基于哈密頓變分原理,獲得了懸臂夾層圓柱殼的運(yùn)動(dòng)微分方程,采用譜有限元法研究了長(zhǎng)徑比和徑厚比對(duì)固有頻率和損耗因子的影響。

    本文針對(duì)懸臂邊界條件下的環(huán)狀局部覆蓋粘彈層合薄壁圓柱殼,基于薄殼理論和粘彈性理論,采用傳遞矩陣法進(jìn)行了振動(dòng)特性研究,求解得到了結(jié)構(gòu)的高階固有頻率、三維模態(tài)振型和模態(tài)損耗因子,并通過(guò)有限元方法進(jìn)行了比較,最后分析了約束阻尼層位置對(duì)整體結(jié)構(gòu)振動(dòng)特性的影響。

    1 整體結(jié)構(gòu)的模態(tài)特性分析

    1.1 環(huán)狀覆蓋粘彈層合懸臂薄壁圓柱殼模型

    如圖1所示為局部環(huán)狀覆蓋粘彈性層合懸臂薄壁圓柱殼的模型,在柱坐標(biāo)系Oxθz中,基層長(zhǎng)度、中面半徑和壁厚表示為L(zhǎng)、R1和H1,粘彈性阻尼層和約束層的長(zhǎng)度均為L(zhǎng)2-L1,阻尼層的中面半徑和壁厚表示為R2和H2,約束層的中面半徑和和厚度表示為R3和H3,uj、vj和wj表示基層和約束層在x、y和z方向上的振動(dòng)位移。

    圖1 局部覆蓋粘彈層合懸臂薄壁圓柱殼模型Fig.1 Model of a partially covered viscoelastic laminated cantilever cylindrical shell

    1.2 變形協(xié)調(diào)關(guān)系和層間作用力

    粘彈性阻尼層通過(guò)基層和約束層的彎曲變形,起到振動(dòng)耗能減振,由基爾霍夫平面假定和樂(lè)甫殼體理論[21],假設(shè)分別表示基層和約束層的中面沿x、θ和z方向上的振動(dòng)位移,阻尼層內(nèi)沿x和θ方向的剪應(yīng)力為:

    圖2 層間徑向相互作用力和阻尼層內(nèi)的剪切力Fig.2 Radial interactions between layers and shear force in damping layer

    如圖2所示,對(duì)于粘彈阻尼層,當(dāng)僅考慮法向慣性力時(shí)的平衡方程為

    1.3 振動(dòng)特性求解

    首先,對(duì)基層(即j=1)進(jìn)行分析,忽略任一區(qū)段的下角標(biāo)i,振動(dòng)微分方程表示為

    內(nèi)力與內(nèi)力矩和中面位移間的關(guān)系為:

    式中:E1為楊氏模量,μ1為泊松比,K1為薄膜剛度,D1為彎曲剛度。

    對(duì)于面內(nèi)剪力、橫向剪力和中法線(xiàn)繞θ軸的轉(zhuǎn)角,表達(dá)式分別為:

    令m表示軸向半波數(shù),n表示周向波數(shù),ωmn表示固有頻率,定義位移解的形式為:

    定義任一截面ξi上狀態(tài)向量中的8個(gè)元素為:

    則位移、轉(zhuǎn)角、內(nèi)力和內(nèi)力矩表示為:

    整理得到一階狀態(tài)向量微分方程組為:

    同理,對(duì)于約束層薄壁圓柱殼,任一截面上狀態(tài)向量為:

    一階狀態(tài)向量微分方程組表示為:

    對(duì)于基層和約束層,傳遞矩陣的關(guān)系為:

    消除變量后,定義狀態(tài)向量中的12個(gè)元素為:

    將變形協(xié)調(diào)關(guān)系式(1)和法向平衡方程式(2),代入基層與約束層的振動(dòng)微分方程中,對(duì)變量進(jìn)行整理,推導(dǎo)得到一階微分方程組為:

    式中:U是12×12階的矩陣。

    如圖1所示,對(duì)整體結(jié)構(gòu)進(jìn)行求解時(shí),可沿縱向分成a、b、c三段,將狀態(tài)向量式(14)進(jìn)行處理,令

    對(duì)于a段和b段的薄壁圓柱殼采用(10)式和(12)式描述,對(duì)于覆蓋區(qū)域采用(14)式描述。在軸向截面x=0、L1、L2和L位置,狀態(tài)向量之間存在關(guān)系為:

    式中:Ta、Tb和 Tc分別為 a、b和 c段的傳遞矩陣。

    消去 Y1(La)和 Y1(Lb),整理得

    對(duì)于某一階模態(tài),整合各個(gè)變量后可以得到ξ0層的徑向位移,然后利用每個(gè)區(qū)段的子傳遞矩陣關(guān)系,可以得到任一截面處的振動(dòng)位移關(guān)系,最后求出整體結(jié)構(gòu)的模態(tài)振型。

    則固有頻率和模態(tài)損耗因子分別為:

    2 算例求解與分析

    在懸臂邊界條件下,通過(guò)算例對(duì)模態(tài)特性進(jìn)行求解,選取構(gòu)件的幾何參數(shù)與材料參數(shù)為:L=0.3 m,L1=0.1m,L2=0.2 m,R1=0.3 m,H1=0.003 m,H2=0.001 m,H3=0.002 m,基層和約束層的彈性模量為 E1=E3=70 GPa,密度為 ρ1=ρ3=2 700 kg/m3,泊松比為 μ1=μ3=0.3,粘彈性層的復(fù)剪切模量為 G2=0.896泊松比 μ2=0.498,ρ2=999 kg/m3。

    (1)中部覆蓋阻尼層薄壁圓柱殼的振動(dòng)特性

    對(duì)于中部覆蓋約束阻尼層的薄壁圓柱殼,求得的固有頻率和模態(tài)損耗因子如表1所示,三維模態(tài)振型圖如表2所示,并通過(guò)有限元程序ANSYS的結(jié)果進(jìn)行比較,其中單元類(lèi)型取三維高階實(shí)體單元SOLID186,采用VROATA命令從截面繞軸旋轉(zhuǎn)生成模型,從相鄰面掃掠體的方法生成網(wǎng)格,包括2 400個(gè)單元和14 580個(gè)節(jié)點(diǎn)。

    由表1可以看出,隨著階次的增加,利用傳遞矩陣法和有限元法得到的前25階模態(tài)振型變化關(guān)系一致,以周向模態(tài)振型變化為主,軸向模態(tài)振型出現(xiàn)在m=1、2和3,周向振型最高階出現(xiàn)在n=14,其中,前4階三維模態(tài)振型依次出現(xiàn)在m=1和n=5、6、4、7;第1階模態(tài)振型發(fā)生在軸向半波數(shù)為1和周向波數(shù)為5,對(duì)應(yīng)的固有頻率分別為350.36 Hz和347.8 Hz,模態(tài)損耗因子分別為0.013 44和0.016 5,通過(guò)誤差比較可知,固有頻率的最大誤差值小于2%,模態(tài)損耗因子的最大誤差值小于7%,說(shuō)明了采用傳遞矩陣法對(duì)固有頻率和模態(tài)損耗因子求解的有效性。

    表1 局部覆蓋粘彈性層合薄壁圓柱殼的固有頻率和模態(tài)損耗因子Tab.1 Natural frequencies and modal loss factors of a partially covered viscoelastic laminated cylindrical shell

    由表1可知,前4階模態(tài)發(fā)生在軸向半波數(shù)為m=1和周向波數(shù)為n=4、5、6、7,則由兩種方法求解m=1、2和n=4、5、6、7得到的三維模態(tài)振型如表2所示,整體結(jié)構(gòu)的振動(dòng)位移變化關(guān)系一致,反映了構(gòu)件在懸臂端的振動(dòng)狀態(tài),即在懸臂端處節(jié)點(diǎn)的振動(dòng)位移最大,其中最低階(1,5)對(duì)應(yīng)的三維模態(tài)振型表現(xiàn)為懸臂端周向波數(shù)為5的花瓣形狀,說(shuō)明了使用傳遞矩陣法對(duì)三維模態(tài)振型求解的有效性。

    (2)約束阻尼層位置對(duì)振動(dòng)特性的影響

    當(dāng)基層、阻尼層和約束層的尺寸參數(shù)和材料參數(shù)不變時(shí),僅改變阻尼層的軸向覆蓋位置,當(dāng)軸向半波數(shù)為1時(shí)振動(dòng)特性如圖3~5所示。由圖3和圖4可知,隨著周向波數(shù)的增加,當(dāng)約束阻尼層位于薄壁圓柱殼下部時(shí)固有頻率值最高,固有頻率值減小的順序依次為約束阻尼層位于中部和上部,全覆蓋阻尼層時(shí)的固有頻率值最低,同時(shí)上部覆蓋阻尼層和全覆蓋阻尼層的固有頻率值比較接近,固有頻率最低的一階均出現(xiàn)在第(1,5)階;對(duì)于模態(tài)損耗因子,當(dāng)阻尼層位于薄壁圓柱殼下部時(shí)模態(tài)損耗因子最小,模態(tài)損耗因子增加的順序依次為約束阻尼層位于中部和上部,全覆蓋約束阻尼層時(shí)的模態(tài)損耗因子最高,損耗因子越高時(shí)減振效果越好。這是由于約束阻尼層的影響,對(duì)基層薄壁圓柱殼起到了約束和耗能作用,當(dāng)約束阻尼層在下側(cè)時(shí)對(duì)于整體結(jié)構(gòu)剛度和阻尼影響不大,在上側(cè)時(shí)導(dǎo)致了固有頻率降低和損耗因子升高,也就更接近于全覆蓋約束阻尼層時(shí)的情況。

    表2 局部覆蓋粘彈性層合薄壁圓柱殼的模態(tài)振型Tab.2 Three-dimensional modal shapes of a partially covered viscoelastic laminated cylindrical shell

    圖3 約束阻尼層位置對(duì)固有頻率的影響Fig.3 The influence of constrained damping layer’s position to natural frequencies

    圖4 約束阻尼層位置對(duì)模態(tài)損耗因子的影響Fig.4 The influence of constrained damping layer’s position to modal loss factors

    圖5 約束阻尼層位置不同時(shí)整體結(jié)構(gòu)的三維模態(tài)振型Fig.5 Three-dimensional modal shapes with different position of constrained damping layer

    當(dāng)約束阻尼層分別位于薄壁圓柱殼的上部、中部、下部和全覆蓋時(shí),以最低階固有頻率對(duì)應(yīng)的振型(1,5)為例,如圖5和表2所示,三維模態(tài)振型均表現(xiàn)為懸臂端周向波數(shù)為5的花瓣形狀,這是由于下端是固定的,對(duì)于結(jié)構(gòu)有較強(qiáng)的約束,而對(duì)于上部約束的影響較小,引起懸臂端的振動(dòng)位移最大,這一階也是整體結(jié)構(gòu)最容易發(fā)生共振的一階。

    3 結(jié) 論

    在懸臂邊界條件下,通過(guò)算例求解得到了環(huán)狀局部覆蓋粘彈層合薄壁圓柱殼的模態(tài)振動(dòng)特性,并通過(guò)有限元法進(jìn)行了比較,主要結(jié)論如下:

    (1)對(duì)于中部覆蓋約束阻尼層的薄壁圓柱殼,隨著階次的增加,采用兩種方法得到的固有頻率的誤差值小于2%,模態(tài)損耗因子的誤差值小于7%,三維模態(tài)振型的變化關(guān)系一致,并且以周向模態(tài)的振動(dòng)為主,前4階三維模態(tài)振型依次出現(xiàn)在m=1和n=5、6、4、7,固有頻率最低的(1,5)階三維振型表現(xiàn)為周向波數(shù)為5的花瓣形狀,在懸臂端的振動(dòng)位移最大,說(shuō)明了采用傳遞矩陣法對(duì)于此類(lèi)結(jié)構(gòu)模態(tài)振動(dòng)特性分析的有效性。

    (2)當(dāng)約束阻尼層的覆蓋位置不同時(shí),當(dāng)位于下部時(shí)固有頻率值最高,固有頻率依次減小的順序?yàn)榧s束阻尼層位于中部和上部,全覆蓋時(shí)固有頻率最低,固有頻率最低階模態(tài)振型均出現(xiàn)在(1,5)階;對(duì)于模態(tài)損耗因子,當(dāng)約束阻尼層位于下部時(shí)損耗因子最小,數(shù)值結(jié)果依次增加的順序?yàn)榧s束阻尼層位于中部和上部,全覆蓋時(shí)模態(tài)損耗因子最高。因此,從結(jié)構(gòu)優(yōu)化設(shè)計(jì)和減振效果的角度出發(fā),當(dāng)約束阻尼層位于上部時(shí)的減振效果比較好,應(yīng)根據(jù)實(shí)際情況對(duì)構(gòu)件進(jìn)行減振處理。

    [1]黃志誠(chéng),秦朝燁,褚福磊.附加粘彈阻尼層的薄壁構(gòu)件振動(dòng)問(wèn)題研究綜述[J].振動(dòng)與沖擊,2014,33(07):105-113.Huang Zhicheng,Qin Zhaoye,Chu Fulei.A review about vibration problems of thin-walled structures with viscoelastic damping layer[J].Journal of Vibration and Shock,2014,33(07):105-113.

    [2]俞孟薩,劉延利,廖彬彬.腔室內(nèi)部聲場(chǎng)與結(jié)構(gòu)振動(dòng)耦合特性及噪聲控制研究綜述[J].船舶力學(xué),2012,16(1-2):191-201.Yu Mengsa,Liu Yanli,Liao Binbin.Summarization on the research of coupling character between cavity acoustics and structure vibration together with noise control[J].Journal of Ship Mechanics,2012,16(1-2):191-201.

    [3]王 宇,羅 忠.薄壁圓柱殼構(gòu)件受迫振動(dòng)的響應(yīng)特征研究[J].振動(dòng)與沖擊,2015,34(07):103-108.Wang Yu,Luo Zhong.Study on response characteristics for forced vibrations of thin cylindrical shell[J].Journal of Vibration and Shock,2015,34(07):103-108.

    [4]王 宇,翟敬宇,李 暉,等.模態(tài)數(shù)量對(duì)薄壁短圓柱殼振動(dòng)響應(yīng)分析的影響[J].噪聲與振動(dòng)控制,2014,34(02):50-55.Wang Yu,Zhai Jingyu,Li Hui,et al.Influence of modal numbers of short thin cylindrical shell on its forced vibration response[J].Noise and Vibration Control,2014,34(02):50-55.

    [5]韓清凱,王 宇,李學(xué)軍.旋轉(zhuǎn)薄壁圓柱殼的高節(jié)徑振動(dòng)特性以及篦齒結(jié)構(gòu)的影響[J].中國(guó)科學(xué):物理學(xué) 力學(xué) 天文學(xué),2013,43(04):436-458.Han Qingkai,Wang Yu,Li Xuejun.High nodal diameter vibration characteristics of rotating shell and the effects of its sealing teeth[J].Scientia Sinica Physica,Mechanica&Astronomica,2013,43(04):436-458.

    [6]解妙霞,陳花玲,吳九匯.圓柱殼高頻彎曲振動(dòng)的能量有限元分析[J].西安交通大學(xué)學(xué)報(bào),2008,42(09):1113-1116.Xie Miaoxia,Chen Hualing,Wu Jiuhui.Energy finite element analysis to high-frequency bending vibration in cylindrical shell[J].Journal of Xi’an Jiaotong University,2008,42(09):1113-1116.

    [7]漆萬(wàn)鵬,侯 磊.正交各向異性圓柱殼的軸對(duì)稱(chēng)和梁式振動(dòng)分析[J].船舶力學(xué),2010,14(8):908-914.Qi Wanpeng,Hou Lei.Axisymmetric and beamlike vibration of orthotropic circular cylindrical shell[J].Journal of Ship Mechanics,2010,14(8):908-914.

    [8]李 暉,孫 偉,許 卓,等.基于激光旋轉(zhuǎn)掃描的約束態(tài)薄壁圓柱殼模態(tài)振型測(cè)試新方法[J].振動(dòng)與沖擊,2014,33(16):155-159.Li Hui,Sun Wei,Xu Zhuo,et al.Experimental method of laser rotating sanning to measure mode shapes of constrained thin cylindrical shell[J].Journal of Vibration and Shock,2014,33(16):155-159.

    [9]Wang H J,Chen L W.Finite element dynamic analysis of orthotropic cylindrical shells with a constrained damping layer[J].Finite Elements in Analysis and Design,2004,40(7):737-755.

    [10]Chen L H,Huang S C.Vibrations of a cylindrical shell with partially constrained layer damping(CLD)treatment[J].International Journal of Mechanical Sciences,1999,41(12):1485-1498.

    [11]Chen L H,Huang S C.Vibration attenuation of a cylindrical shell with constrained layer damping strip treatment[J].Computers and Structures,2001,79(14):1355-1362.

    [12]Sainsbury M G,Mastib R S.Vibration damping of cylindrical shells using strain-energy-based distribution of an add-on viscoelastic treatment[J].Finite Elements in Analysis and Design,2007,43:175-192.

    [13]Mastib R S,Sainsbury M G.Vibration damping of cylindrical shells partially coated with a constrained viscoelastic treatment having a standoff layer[J].Thin-Walled Structures,2005,43:1355-1379.

    [14]Mohammadi F,Sedaghati R.Linear and nonlinear vibration analysis of sandwich cylindrical shell with constrained viscoelastic core layer[J].International Journal of Mechanical Sciences,2012,54:156-171.

    [15]李恩奇,李道奎,唐國(guó)金,等.基于傳遞函數(shù)方法的局部覆蓋環(huán)狀CLD圓柱殼動(dòng)力學(xué)分析[J].航空學(xué)報(bào),2007,28(6):1487-1493.Li Enqi,Li Daokui,Tang Guojin,et al.Dynamic analysis of cylindrical shell with partially covered ring-shape constrained layer damping by the transfer function method[J].Acta Aeronautica et Astronautica Sinica,2007,28(6):1487-1493

    [16]楊少紅,王安穩(wěn).基于混合分層理論的粘彈層合圓柱殼的阻尼振動(dòng)和層間應(yīng)力研究[J].工程力學(xué),2008,25(3):38-42.Yang Shaohong,Wang Anwen.Damping vibration and interlaminar stresses atudy using mixed layerwise theories for viscoelastic laminated cylindrical shell[J].Engineering Mechanics,2008,25(3):38-42.

    [17]王祖華,彭 旭,謝官模.部分敷設(shè)阻尼層的有限長(zhǎng)環(huán)肋圓柱殼聲輻射計(jì)算方法研究[J].船舶力學(xué),2010,14(4):430-438.Wang Zuhua,Peng Xu,Xie Guanmo.Sound radiation from submerged finite ringstiffened cy-lindrical shells,partially covered with compliant layers[J].Journal of Ship Mechanics,2010,14(4):430-438.

    [18]曹雄濤,張志誼,華宏星.約束層阻尼圓柱殼的自由振動(dòng)[J].應(yīng)用數(shù)學(xué)和力學(xué),2011,32(4):470-482.Cao Xiongtao,Zhang Zhiyi,Hua Hongxing.Free vibration of circular cylindrical shell with constrained layer damping[J].Applied Mathematics and Mechanics,2011,32(4):470-482.

    [19]李 超,李以農(nóng),施 磊,等.圓柱殼體阻尼材料布局拓?fù)鋬?yōu)化研究[J].振動(dòng)與沖擊,2012,31(4):48-52.Li Chao,Li Yinong,Shi Lei,et al.Topological optimization for placement of damping material on cylindrical shells[J].Journal of Vibration and Shock,2012,31(4):48-52.

    [20]王 淼,孟 光,方之楚,等.夾層圓柱殼振動(dòng)的譜有限元分析[J],應(yīng)用力學(xué)學(xué)報(bào),2009,26(4):715-720.Wang Miao,Meng Guang,Fang Zhichu,et al.Spectral finite element analysis for vibration of sandwiched shell[J].Chinese Journal of Applied Mechanics,2009,26(4):715-720.

    [21]Soedel W.Vibrations of shells and plates[M].CRC Press,2004.

    [22]向 宇,黃玉盈,黃健強(qiáng).一種新型齊次擴(kuò)容精細(xì)積分法[J].華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2002,3(11):74-76.Xiang Yu,Huang Yuying,Huang Jianqiang.A method of homogenization of high precision direct integration[J].Journal Huazhong University of Science and Technology(Nature Science Edition),2002,30(11):74-76.

    Study on vibration characteristics of a ring covered viscoelastic laminated cantilever thin cylindrical shell

    WANG Yu1,ZHAI Jing-yu2,LI Chang1,SONG Hua1
    (1.School of Mechanical Engineering,University of Science and Technology Liaoning,Anshan 114051,China;2.School of Mechanical Engineering,Dalian University of Technology,Dalian 116024,China)

    Vibration characteristics of a ring covered viscoelastic laminated cantilever thin cylindrical shell are studied by transfer matrix method.A first-order ordinary differential equations for the host shell and the constraining layer are established based on Love’s shell theory combining with deformation coordination of viscoelastic damping layer and interactional forces,and an integrated transfer matrix of the whole component is derived by transfer matrix method.Then,the natural frequencies,loss factors and three-dimensional shapes are derived by high-precise integration method,and its vibration characteristics are discussed when the covered position is different.Finally,the results are compared with results of finite element method.The example results show that transfer matrix method is effective to solve the vibration characteristics of the ring covered viscoelastic laminated cantilever thin cylindrical shell,and the main modal vibration occurs in the circumferential direction.The three-dimensional modal shape corresponding to the lowest order natural frequency is(1,5),and the cantilever’s vibration displacement is maximum.At the same time,the covered location of constrained damping layer has obvious impact on modal characteristics of the thin cylindricalshell.

    viscoelastic laminated thin cylindrical shell;cantilever boundary condition;ring covered;three-dimensional modal shape

    O326 TH113.1

    A

    10.3969/j.issn.1007-7294.2017.11.009

    1007-7294(2017)11-1393-11

    2017-05-18

    遼寧省科技廳博士科研啟動(dòng)基金項(xiàng)目(201601295);遼寧省教育廳自由探索重點(diǎn)項(xiàng)目(2016TSZD07);遼寧科技大學(xué)優(yōu)秀人才項(xiàng)目

    王 宇(1979-),男,博士,講師,E-mail:wangyu435@126.com;

    翟敬宇(1983-),男,博士,講師。

    猜你喜歡
    模態(tài)振動(dòng)
    振動(dòng)的思考
    噴水推進(jìn)高速艇尾部振動(dòng)響應(yīng)分析
    This “Singing Highway”plays music
    振動(dòng)攪拌 震動(dòng)創(chuàng)新
    中立型Emden-Fowler微分方程的振動(dòng)性
    車(chē)輛CAE分析中自由模態(tài)和約束模態(tài)的應(yīng)用與對(duì)比
    國(guó)內(nèi)多模態(tài)教學(xué)研究回顧與展望
    高速顫振模型設(shè)計(jì)中顫振主要模態(tài)的判斷
    基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識(shí)別
    UF6振動(dòng)激發(fā)態(tài)分子的振動(dòng)-振動(dòng)馳豫
    免费黄频网站在线观看国产| 国产午夜精品一二区理论片| 丝袜美腿诱惑在线| 少妇人妻 视频| 亚洲精品国产色婷婷电影| 国产1区2区3区精品| 美女国产高潮福利片在线看| 无遮挡黄片免费观看| 国产精品.久久久| 日韩一本色道免费dvd| 国产欧美日韩综合在线一区二区| 日本爱情动作片www.在线观看| 亚洲五月色婷婷综合| 欧美97在线视频| 亚洲人成77777在线视频| 青草久久国产| 丰满迷人的少妇在线观看| 美女国产高潮福利片在线看| 无遮挡黄片免费观看| 久久精品人人爽人人爽视色| 在线观看免费午夜福利视频| 国产免费视频播放在线视频| 亚洲av成人不卡在线观看播放网 | 黄色视频在线播放观看不卡| 精品一品国产午夜福利视频| 在线看a的网站| 亚洲一级一片aⅴ在线观看| 91精品国产国语对白视频| 亚洲成人一二三区av| 国产又色又爽无遮挡免| 只有这里有精品99| 最近中文字幕高清免费大全6| 国产乱来视频区| 国产av一区二区精品久久| 久久精品国产亚洲av高清一级| 免费不卡黄色视频| 日韩欧美精品免费久久| tube8黄色片| 大陆偷拍与自拍| bbb黄色大片| 免费观看性生交大片5| 久久久精品区二区三区| 99re6热这里在线精品视频| 欧美日韩成人在线一区二区| 黄片播放在线免费| 亚洲国产成人一精品久久久| 成人亚洲精品一区在线观看| 操美女的视频在线观看| 老司机深夜福利视频在线观看 | 亚洲精品国产av蜜桃| 欧美人与善性xxx| 国产黄色免费在线视频| 黄片播放在线免费| av线在线观看网站| 国产熟女欧美一区二区| 一级片'在线观看视频| 超碰成人久久| 两个人免费观看高清视频| 国产精品 国内视频| 午夜激情久久久久久久| 久久久久国产精品人妻一区二区| 久久精品aⅴ一区二区三区四区| 满18在线观看网站| 亚洲精品国产一区二区精华液| av电影中文网址| 在线观看人妻少妇| 国产女主播在线喷水免费视频网站| 成年女人毛片免费观看观看9 | 精品少妇一区二区三区视频日本电影 | 高清视频免费观看一区二区| 欧美日韩福利视频一区二区| 黄色毛片三级朝国网站| 亚洲伊人色综图| 欧美最新免费一区二区三区| 成人国语在线视频| 一边摸一边抽搐一进一出视频| 日韩电影二区| 免费观看av网站的网址| 国产男女内射视频| 天堂俺去俺来也www色官网| 国产不卡av网站在线观看| 亚洲精品视频女| 精品一区在线观看国产| 久久精品久久久久久久性| 亚洲国产精品999| 黄片无遮挡物在线观看| 久久精品久久精品一区二区三区| 亚洲av在线观看美女高潮| 午夜福利免费观看在线| 亚洲欧美成人综合另类久久久| 精品少妇内射三级| 亚洲欧美色中文字幕在线| 成人毛片60女人毛片免费| 国产男人的电影天堂91| 国产乱来视频区| 91精品三级在线观看| 亚洲精品视频女| 黄片播放在线免费| 女的被弄到高潮叫床怎么办| 久久久久久久久久久免费av| 性少妇av在线| 国产 精品1| 国产一卡二卡三卡精品 | av女优亚洲男人天堂| 国产在线一区二区三区精| videosex国产| 久久精品久久久久久噜噜老黄| 久久久久久久精品精品| 日韩中文字幕欧美一区二区 | 街头女战士在线观看网站| 看免费av毛片| 精品人妻一区二区三区麻豆| 国产精品av久久久久免费| 免费不卡黄色视频| 热99国产精品久久久久久7| 美女中出高潮动态图| 电影成人av| 搡老乐熟女国产| 免费看av在线观看网站| 最近中文字幕2019免费版| 性少妇av在线| 久久人人爽人人片av| 人人妻人人爽人人添夜夜欢视频| 一区二区av电影网| a 毛片基地| 国产成人午夜福利电影在线观看| 亚洲天堂av无毛| 久久国产亚洲av麻豆专区| 亚洲美女搞黄在线观看| 人体艺术视频欧美日本| av一本久久久久| 国产精品一二三区在线看| 国产在线视频一区二区| 精品久久久久久电影网| 韩国av在线不卡| 亚洲精华国产精华液的使用体验| 观看av在线不卡| 好男人视频免费观看在线| 国产精品麻豆人妻色哟哟久久| 色网站视频免费| 可以免费在线观看a视频的电影网站 | 人体艺术视频欧美日本| 成人漫画全彩无遮挡| 亚洲国产成人一精品久久久| 亚洲av欧美aⅴ国产| 久久久久精品国产欧美久久久 | 亚洲精品久久午夜乱码| 国产成人系列免费观看| 日韩大码丰满熟妇| 日日摸夜夜添夜夜爱| 高清在线视频一区二区三区| av网站免费在线观看视频| 国产人伦9x9x在线观看| 精品亚洲乱码少妇综合久久| 麻豆av在线久日| 我的亚洲天堂| 亚洲精品一区蜜桃| 亚洲av欧美aⅴ国产| 天美传媒精品一区二区| 国产精品麻豆人妻色哟哟久久| 日韩av在线免费看完整版不卡| 日韩中文字幕视频在线看片| 欧美日韩精品网址| 成人亚洲精品一区在线观看| 搡老乐熟女国产| 天天躁夜夜躁狠狠躁躁| 天天躁夜夜躁狠狠久久av| 多毛熟女@视频| 亚洲成人一二三区av| 亚洲成人av在线免费| av网站在线播放免费| 激情视频va一区二区三区| 精品酒店卫生间| 日韩中文字幕视频在线看片| 在线亚洲精品国产二区图片欧美| 中文天堂在线官网| 欧美日韩国产mv在线观看视频| 美女午夜性视频免费| 汤姆久久久久久久影院中文字幕| 精品久久久精品久久久| 亚洲精品一二三| 欧美日韩福利视频一区二区| 亚洲五月色婷婷综合| 日韩熟女老妇一区二区性免费视频| 午夜激情久久久久久久| 99久久人妻综合| 亚洲一码二码三码区别大吗| 国产99久久九九免费精品| 欧美精品高潮呻吟av久久| 亚洲精品日韩在线中文字幕| 在线观看免费高清a一片| 精品视频人人做人人爽| 最黄视频免费看| 国产老妇伦熟女老妇高清| 欧美激情 高清一区二区三区| 极品人妻少妇av视频| 亚洲精品日本国产第一区| 亚洲视频免费观看视频| 国产精品亚洲av一区麻豆 | 日韩不卡一区二区三区视频在线| 曰老女人黄片| 国产男女超爽视频在线观看| 久久热在线av| 国产又爽黄色视频| 观看av在线不卡| 五月天丁香电影| 亚洲男人天堂网一区| 日韩一卡2卡3卡4卡2021年| 一级a爱视频在线免费观看| 岛国毛片在线播放| 欧美中文综合在线视频| 精品亚洲乱码少妇综合久久| 哪个播放器可以免费观看大片| 午夜福利影视在线免费观看| 老司机影院毛片| 国产精品三级大全| 亚洲国产精品国产精品| 国产精品无大码| 午夜91福利影院| 午夜福利乱码中文字幕| 亚洲成国产人片在线观看| 日韩中文字幕欧美一区二区 | 久久久久国产一级毛片高清牌| 另类精品久久| 免费女性裸体啪啪无遮挡网站| 黑人欧美特级aaaaaa片| 欧美在线黄色| 操美女的视频在线观看| 亚洲精品国产av蜜桃| 亚洲国产欧美一区二区综合| 亚洲婷婷狠狠爱综合网| 欧美黑人欧美精品刺激| 国产精品三级大全| 亚洲欧美色中文字幕在线| 一区在线观看完整版| 亚洲欧美一区二区三区久久| 久久亚洲国产成人精品v| 99久久精品国产亚洲精品| 国产亚洲一区二区精品| 黄频高清免费视频| 精品国产国语对白av| 亚洲人成网站在线观看播放| 国产一区亚洲一区在线观看| 女的被弄到高潮叫床怎么办| 久久精品亚洲av国产电影网| 一边摸一边抽搐一进一出视频| 亚洲av日韩精品久久久久久密 | 久久 成人 亚洲| 亚洲精品久久成人aⅴ小说| 少妇被粗大猛烈的视频| 久久精品久久久久久久性| 大话2 男鬼变身卡| 一二三四中文在线观看免费高清| 免费看不卡的av| 国产精品亚洲av一区麻豆 | av又黄又爽大尺度在线免费看| 久久久久久人人人人人| 国产极品天堂在线| 亚洲国产精品成人久久小说| av.在线天堂| 亚洲色图综合在线观看| h视频一区二区三区| 日韩不卡一区二区三区视频在线| 一本大道久久a久久精品| 日韩电影二区| 精品国产乱码久久久久久男人| 久久天躁狠狠躁夜夜2o2o | 免费久久久久久久精品成人欧美视频| 侵犯人妻中文字幕一二三四区| 黄色视频不卡| 久久精品久久久久久噜噜老黄| 九色亚洲精品在线播放| 欧美黑人欧美精品刺激| 一本—道久久a久久精品蜜桃钙片| 欧美日韩精品网址| 看免费成人av毛片| 中国三级夫妇交换| 黄网站色视频无遮挡免费观看| 99re6热这里在线精品视频| 欧美人与性动交α欧美精品济南到| 亚洲色图综合在线观看| 国产免费又黄又爽又色| 青草久久国产| 又黄又粗又硬又大视频| 日韩一区二区视频免费看| 只有这里有精品99| 大片免费播放器 马上看| h视频一区二区三区| 两性夫妻黄色片| 好男人视频免费观看在线| 亚洲欧美色中文字幕在线| 男女免费视频国产| 精品第一国产精品| 999久久久国产精品视频| 纵有疾风起免费观看全集完整版| 亚洲情色 制服丝袜| 国产乱人偷精品视频| 欧美精品一区二区大全| 久久久久久久精品精品| 男女午夜视频在线观看| 电影成人av| 国产99久久九九免费精品| 国产一区有黄有色的免费视频| 亚洲熟女精品中文字幕| a 毛片基地| 成人亚洲欧美一区二区av| 亚洲欧洲日产国产| 亚洲美女搞黄在线观看| 这个男人来自地球电影免费观看 | 丝袜美腿诱惑在线| 青青草视频在线视频观看| 高清在线视频一区二区三区| 又大又爽又粗| 亚洲四区av| 国产精品.久久久| 国产一区二区 视频在线| 一级爰片在线观看| 亚洲欧洲国产日韩| 国产精品一二三区在线看| 大话2 男鬼变身卡| 嫩草影院入口| avwww免费| 日韩不卡一区二区三区视频在线| 国产亚洲av高清不卡| 欧美激情高清一区二区三区 | videosex国产| 欧美另类一区| 亚洲成人免费av在线播放| 哪个播放器可以免费观看大片| 下体分泌物呈黄色| 国产精品av久久久久免费| 制服人妻中文乱码| 韩国精品一区二区三区| 两性夫妻黄色片| 狂野欧美激情性bbbbbb| 99国产综合亚洲精品| 亚洲,一卡二卡三卡| 考比视频在线观看| 可以免费在线观看a视频的电影网站 | 老汉色av国产亚洲站长工具| 亚洲欧洲日产国产| 少妇精品久久久久久久| 美女主播在线视频| 99热全是精品| 男女下面插进去视频免费观看| 中文字幕高清在线视频| 亚洲,欧美精品.| 夫妻午夜视频| 亚洲成人免费av在线播放| 免费在线观看黄色视频的| 丝袜脚勾引网站| 国产欧美亚洲国产| 男女边吃奶边做爰视频| 亚洲美女搞黄在线观看| 亚洲精品久久久久久婷婷小说| 午夜91福利影院| 99国产精品免费福利视频| 精品午夜福利在线看| 老汉色∧v一级毛片| 午夜福利网站1000一区二区三区| 日本av手机在线免费观看| 国产免费又黄又爽又色| 99热国产这里只有精品6| 韩国av在线不卡| 免费观看av网站的网址| 精品少妇久久久久久888优播| 一级爰片在线观看| 久久久久精品久久久久真实原创| 777久久人妻少妇嫩草av网站| 精品人妻熟女毛片av久久网站| av又黄又爽大尺度在线免费看| 中文字幕亚洲精品专区| 成人毛片60女人毛片免费| 91精品伊人久久大香线蕉| 亚洲成人国产一区在线观看 | √禁漫天堂资源中文www| 男女午夜视频在线观看| 只有这里有精品99| av国产久精品久网站免费入址| 一区福利在线观看| 亚洲伊人久久精品综合| 国产亚洲一区二区精品| 国产福利在线免费观看视频| 丁香六月欧美| 久久国产精品男人的天堂亚洲| 最黄视频免费看| 欧美精品亚洲一区二区| 精品国产露脸久久av麻豆| svipshipincom国产片| 熟妇人妻不卡中文字幕| av在线观看视频网站免费| 国产老妇伦熟女老妇高清| 亚洲成色77777| 中文字幕最新亚洲高清| 亚洲第一av免费看| 成人黄色视频免费在线看| 无限看片的www在线观看| 男女床上黄色一级片免费看| 在线观看免费午夜福利视频| 一二三四在线观看免费中文在| 19禁男女啪啪无遮挡网站| 国产男女内射视频| 三上悠亚av全集在线观看| 亚洲七黄色美女视频| 国产成人av激情在线播放| 天天躁狠狠躁夜夜躁狠狠躁| a级毛片在线看网站| 国语对白做爰xxxⅹ性视频网站| 亚洲精华国产精华液的使用体验| 色吧在线观看| 婷婷成人精品国产| 美女脱内裤让男人舔精品视频| 日本爱情动作片www.在线观看| 欧美黑人精品巨大| 国产97色在线日韩免费| 中文字幕亚洲精品专区| 国产黄色免费在线视频| 99久国产av精品国产电影| 女人久久www免费人成看片| av一本久久久久| 一级爰片在线观看| 午夜老司机福利片| 国产精品一二三区在线看| 中文天堂在线官网| 一级毛片我不卡| 在线天堂中文资源库| 在线观看国产h片| 日韩视频在线欧美| 精品一区二区三区四区五区乱码 | 男男h啪啪无遮挡| 亚洲国产毛片av蜜桃av| 欧美日韩成人在线一区二区| 五月天丁香电影| 国产成人精品在线电影| 精品午夜福利在线看| 日日啪夜夜爽| 观看av在线不卡| 夜夜骑夜夜射夜夜干| 高清视频免费观看一区二区| 国产精品 欧美亚洲| 91成人精品电影| 人人妻人人澡人人看| 国产无遮挡羞羞视频在线观看| videosex国产| 人人妻人人添人人爽欧美一区卜| 久久狼人影院| 成人国语在线视频| 综合色丁香网| 久久亚洲国产成人精品v| www.精华液| 亚洲国产看品久久| 18禁裸乳无遮挡动漫免费视频| 伊人久久大香线蕉亚洲五| 日本午夜av视频| 麻豆乱淫一区二区| 天天躁夜夜躁狠狠久久av| 免费日韩欧美在线观看| 久久久久精品人妻al黑| 9色porny在线观看| 国产又色又爽无遮挡免| svipshipincom国产片| 午夜免费观看性视频| 久久毛片免费看一区二区三区| 大码成人一级视频| 亚洲图色成人| 少妇人妻 视频| 三上悠亚av全集在线观看| av有码第一页| 久久青草综合色| 日本黄色日本黄色录像| 日韩成人av中文字幕在线观看| 国产欧美亚洲国产| av在线播放精品| 成人免费观看视频高清| 久久99热这里只频精品6学生| 另类亚洲欧美激情| 亚洲精品国产色婷婷电影| 国产成人午夜福利电影在线观看| av网站在线播放免费| 精品人妻一区二区三区麻豆| 亚洲av男天堂| 午夜影院在线不卡| 亚洲免费av在线视频| 亚洲精品一区蜜桃| 又大又黄又爽视频免费| 婷婷成人精品国产| 只有这里有精品99| 一级a爱视频在线免费观看| 交换朋友夫妻互换小说| 欧美 亚洲 国产 日韩一| 欧美激情高清一区二区三区 | 赤兔流量卡办理| 免费久久久久久久精品成人欧美视频| av在线老鸭窝| 2018国产大陆天天弄谢| 成年动漫av网址| 日本猛色少妇xxxxx猛交久久| av在线老鸭窝| 一本—道久久a久久精品蜜桃钙片| 国产野战对白在线观看| av片东京热男人的天堂| 精品人妻一区二区三区麻豆| 操美女的视频在线观看| 亚洲图色成人| 少妇的丰满在线观看| 午夜福利视频在线观看免费| 成人国产麻豆网| 免费黄网站久久成人精品| 性少妇av在线| 精品国产一区二区三区四区第35| 成人国语在线视频| 波多野结衣av一区二区av| 亚洲免费av在线视频| 亚洲图色成人| 精品一区二区免费观看| 涩涩av久久男人的天堂| 爱豆传媒免费全集在线观看| 久久婷婷青草| 中文字幕高清在线视频| 日韩,欧美,国产一区二区三区| 99久久综合免费| 免费女性裸体啪啪无遮挡网站| 超碰成人久久| 亚洲国产欧美一区二区综合| 午夜激情av网站| 国产极品天堂在线| 婷婷成人精品国产| 最黄视频免费看| 天天躁夜夜躁狠狠躁躁| 久久性视频一级片| 国产日韩欧美视频二区| 男人操女人黄网站| 欧美成人精品欧美一级黄| av福利片在线| 飞空精品影院首页| 久久精品久久久久久噜噜老黄| 狠狠精品人妻久久久久久综合| 午夜老司机福利片| 极品少妇高潮喷水抽搐| 亚洲成国产人片在线观看| 天天影视国产精品| 日本猛色少妇xxxxx猛交久久| 午夜福利乱码中文字幕| 国产视频首页在线观看| 国产伦理片在线播放av一区| 最近中文字幕高清免费大全6| 国产欧美亚洲国产| 精品国产一区二区久久| 日本爱情动作片www.在线观看| 免费看不卡的av| 男女高潮啪啪啪动态图| 国产不卡av网站在线观看| 成人漫画全彩无遮挡| 亚洲男人天堂网一区| 丁香六月欧美| 丝袜美腿诱惑在线| a 毛片基地| 十八禁人妻一区二区| 国产亚洲av高清不卡| 国产黄频视频在线观看| 黄网站色视频无遮挡免费观看| a级毛片在线看网站| 青春草国产在线视频| 亚洲美女黄色视频免费看| 亚洲精品成人av观看孕妇| 亚洲第一av免费看| 视频在线观看一区二区三区| 亚洲欧美日韩另类电影网站| 午夜日本视频在线| 国产精品久久久久久精品古装| 黄片无遮挡物在线观看| 国产爽快片一区二区三区| 国产一区亚洲一区在线观看| 国产成人午夜福利电影在线观看| 午夜久久久在线观看| 18禁国产床啪视频网站| 国产日韩一区二区三区精品不卡| 久久鲁丝午夜福利片| 色播在线永久视频| 狠狠婷婷综合久久久久久88av| 日韩免费高清中文字幕av| 男女免费视频国产| 99热全是精品| 欧美日韩福利视频一区二区| 精品亚洲成国产av| 国产精品免费大片| 亚洲成人国产一区在线观看 | 亚洲欧洲日产国产| 深夜精品福利| 午夜福利视频精品| 国产成人精品久久久久久| 侵犯人妻中文字幕一二三四区| 午夜av观看不卡| 女人高潮潮喷娇喘18禁视频| 日韩av免费高清视频| 18禁国产床啪视频网站| 亚洲精品av麻豆狂野| 男女免费视频国产| 国产亚洲一区二区精品| 一区二区av电影网| 在线观看一区二区三区激情| 国产精品欧美亚洲77777| 一本一本久久a久久精品综合妖精| 国产老妇伦熟女老妇高清| 菩萨蛮人人尽说江南好唐韦庄| 国产毛片在线视频| 黄色视频在线播放观看不卡| 桃花免费在线播放| 国产在线免费精品| 久久人人爽人人片av| 色综合欧美亚洲国产小说| 久久国产精品大桥未久av| 在线观看www视频免费| 最近中文字幕高清免费大全6| 色精品久久人妻99蜜桃|