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

    基于 Layerwise 理論的共固化粘彈阻尼復(fù)合材料動(dòng)特性分析

    2016-01-12 10:32:51徐超,林松,王立峰
    振動(dòng)與沖擊 2015年1期
    關(guān)鍵詞:粘彈性復(fù)合材料有限元

    基于Layerwise理論的共固化粘彈阻尼復(fù)合材料動(dòng)特性分析

    徐超1,林松1,王立峰2,王建月2

    (1.西北工業(yè)大學(xué)航天學(xué)院,西安710072;2. 航天材料及工藝研究所,北京100076)

    摘要:共固化粘彈性復(fù)合材料兼具結(jié)構(gòu)承載和阻尼減振功能。針對(duì)傳統(tǒng)的混合單元法在應(yīng)用于粘彈性夾層復(fù)合材料結(jié)構(gòu)阻尼性能分析時(shí)存在著前處理困難、計(jì)算規(guī)模大、精度低以及難以考慮正交各向異性鋪層自身損耗能力的缺點(diǎn),推導(dǎo)了一種基于Layerwise離散層理論的四節(jié)點(diǎn)四邊形復(fù)合材料層合板單元,并利用直接復(fù)特征值解法建立了共固化粘彈性復(fù)合材料結(jié)構(gòu)的阻尼性能分析方法。將該方法應(yīng)用于不同的阻尼結(jié)構(gòu),分析結(jié)果與文獻(xiàn)中已公開結(jié)果和混合單元法的計(jì)算結(jié)果進(jìn)行了對(duì)比驗(yàn)證。結(jié)果表明,基于離散層理論的層合板單元具有計(jì)算精度高、前處理建模簡(jiǎn)單和計(jì)算規(guī)模小的優(yōu)點(diǎn),可有效應(yīng)用于復(fù)雜共固化粘彈性復(fù)合材料結(jié)構(gòu)的阻尼性能分析和設(shè)計(jì)。

    關(guān)鍵詞:粘彈性;復(fù)合材料;離散層理論;有限元

    中圖分類號(hào):V414文獻(xiàn)標(biāo)志碼:A

    基金項(xiàng)目:中國(guó)博士后科學(xué)基金面上資助項(xiàng)目(2013M530192); 上海市博士后科研資助項(xiàng)目(12R21413800)

    收稿日期:2013-11-18修改稿收到日期:2014-01-02

    Layerwise dynamic analysis of composite laminates with co-cured viscoelastic damping layers

    XUChao1,LINSong1,WANGLi-feng2,WANGJian-yue2(1. School of Astronautics, Northwestern Polytechnical University,Xi’an 710072, China; 2. Aerospace Research Institute of Material and Processing Technology,Beijing 100076,China)

    Abstract:The Co-curing of thin viscoelastic damping materials in the composite laminates can improve their structural damping with less loss of stiffness and strength. However, classic numerical simulation of the damped structures by using hybrid finite element method demands a cumbersome computational work, especially for complicated layered structures. A new 4-node quadrangular plate finite element based on the discrete layerwise plate theory was developed for dynamics analysis of composite laminates with co-cured viscoelastic layers. The natural frequencies and modal loss factors of different viscoelastic damped laminated beams, cylindrical shells and stiffened plates were derived. The results show that the proposed model can correctly describe the high shear pattern developed inside thin viscoelastic layers and provide results with less computational work.

    Key words:viscoelasticity; damped composite; discrete layerwise theory; finite element analysis

    先進(jìn)纖維增強(qiáng)復(fù)合材料因其突出的比強(qiáng)度、比剛度和可設(shè)計(jì)性強(qiáng)等優(yōu)點(diǎn),已成為航空、航天飛行器的一種主要結(jié)構(gòu)材料。復(fù)合材料自身的損耗因子比金屬材料高約1~2個(gè)數(shù)量級(jí),但其作為結(jié)構(gòu)使用時(shí),多采用整體化成型工藝制造,因而裝配結(jié)合面大為減少,界面摩擦損耗能力減弱,整體結(jié)構(gòu)阻尼值并不高,在實(shí)際使用中仍顯偏低[1]。共固化粘彈性復(fù)合材料是一種新型的結(jié)構(gòu)阻尼復(fù)合材料。所謂共固化,是指在纖維增強(qiáng)復(fù)合材料成型工藝過程中,將高分子粘彈性材料作為特定的鋪層嵌入層合結(jié)構(gòu)中去,然后共同固化成型[2]。它的最大特點(diǎn)是經(jīng)過共固化處理,粘彈性阻尼材料與結(jié)構(gòu)復(fù)合材料合為一體,既充分利用了粘彈性材料的高損耗能力,又不顯著降低復(fù)合材料的剛度和強(qiáng)度性能,還能利用復(fù)合材料的可設(shè)計(jì)性進(jìn)行“事前”阻尼設(shè)計(jì),達(dá)到了結(jié)構(gòu)承載和阻尼減振功能一體化的目的。

    一般說來,粘彈性阻尼材料的模量較復(fù)合材料鋪層低約4~5個(gè)數(shù)量級(jí),作為夾層嵌入層合結(jié)構(gòu)中后,受到相鄰剛度較大的彈性層變形約束,發(fā)生剪切變形,通過分子鏈段間的摩擦耗散振動(dòng)能量。分析共固化粘彈性復(fù)合材料的阻尼性能必須要考慮阻尼層的剪切變形以及沿厚度方向材料的不均勻性。通用有限元軟件中的層合板、殼等結(jié)構(gòu)單元多是基于經(jīng)典層合板理論或一階剪切變形理論,前者忽略了層合板的橫向剪切變形,而后者未能考慮橫向剪切變形沿板厚方向分布的不均勻性,忽略了粘彈性材料層高剪切變形的事實(shí),因而都不能很好的用于共固化粘彈性復(fù)合材料阻尼性能的分析[3]。目前,工程中一般采用“混合單元法”進(jìn)行阻尼結(jié)構(gòu)的分析和設(shè)計(jì)[4-6],即采用三維實(shí)體單元描述阻尼層的剪切變形,采用傳統(tǒng)的層合板殼單元描述復(fù)合材料彈性層,板殼單元和實(shí)體單元之間通過定義節(jié)點(diǎn)偏置或多點(diǎn)約束方程相聯(lián)系。張少輝等[4]對(duì)復(fù)合材料彈性層也采用實(shí)體單元模擬,應(yīng)用全三維實(shí)體模型和模態(tài)應(yīng)變能法分析了共固化粘彈性復(fù)合材料結(jié)構(gòu)的損耗因子。徐超等[5]應(yīng)用混合單元法對(duì)衛(wèi)星飛輪支架進(jìn)行了共固化阻尼減振設(shè)計(jì)。林松等[6]應(yīng)用混合單元法和多目標(biāo)遺傳算法建立了共固化粘彈性復(fù)合材料的阻尼/結(jié)構(gòu)一體化設(shè)計(jì)方法?;旌蠁卧ㄔ趯?shí)際應(yīng)用中主要存在三大不足:一是隨著阻尼層數(shù)增多,有限元模型規(guī)模驟增,建模前處理困難,計(jì)算耗費(fèi)巨大;二是阻尼層厚度較小,體單元長(zhǎng)厚比大,計(jì)算誤差較大;三是修改阻尼層厚度或插入位置時(shí)需要進(jìn)行有限元網(wǎng)格重構(gòu),不易于進(jìn)行優(yōu)化設(shè)計(jì)。

    Reddy[7]提出的Layerwise離散層理論是描述復(fù)合材料層合結(jié)構(gòu)力學(xué)行為的精細(xì)理論?;贚ayerwise理論的結(jié)構(gòu)單元已被成功用于層間應(yīng)力計(jì)算、損傷預(yù)測(cè)等復(fù)合材料力學(xué)問題中。它在統(tǒng)一的位移場(chǎng)描述下,獨(dú)立地考慮各層的變形,并通過引入層間位移連續(xù)性假設(shè),將三維問題退化為二維問題,計(jì)算建模簡(jiǎn)單,易于進(jìn)行優(yōu)化設(shè)計(jì)。文獻(xiàn)[8]基于Layerwise理論推導(dǎo)了多層阻尼梁?jiǎn)卧挠邢拊惺?,分析了共固化粘彈?fù)合材料平面梁的阻尼性能。本文基于Layerwise板理論,推導(dǎo)一種四節(jié)點(diǎn)四邊形等參數(shù)板單元,將其應(yīng)用于共固化粘彈性復(fù)合材料層合結(jié)構(gòu)固有頻率和阻尼性能的分析,并通過與文獻(xiàn)中已有結(jié)果和傳統(tǒng)分析方法的對(duì)比驗(yàn)證方法的有效性。

    1基于Layerwise理論的復(fù)合材料層合板單元

    根據(jù)Reddy的Layerwise理論,首先將復(fù)合材料層合結(jié)構(gòu)沿厚度方向分為若干層。各層既可以是真實(shí)物理材料層,也可以是由多個(gè)材料層構(gòu)成的等效層或者把單層材料再離散的數(shù)學(xué)意義上的數(shù)值層等。為準(zhǔn)確分析復(fù)合材料結(jié)構(gòu)的阻尼性能,要把嵌入的各粘彈性材料層劃分為單獨(dú)層。然后將劃分后的每一層視作考慮面內(nèi)、彎曲和橫向剪切變形的Mindlin板,且各層之間滿足位移連續(xù)條件[9]?;贚ayerwise理論推導(dǎo)有限元單元時(shí),引入如下計(jì)算假設(shè):①忽略層板橫向正應(yīng)變的影響;②各層內(nèi)材料均勻,且符合線彈性或線粘彈性假設(shè);③考慮各層的平動(dòng)質(zhì)量和轉(zhuǎn)動(dòng)慣量。

    1.1位移場(chǎng)

    如圖1所示,將層合板沿厚度方向劃分為n層,各層厚度為hk(k=1,…n)。取k=1層中面所在平面為參考面,建立正交坐標(biāo)系oxyz,其中z軸沿板厚度方向,oxy位于參考平面內(nèi),u,v,w分別為各層沿坐標(biāo)軸方向的位移,θx,θy分別為各層法線繞x,y軸的轉(zhuǎn)角,zk[-hk/2,hk/2]為各層z向的局部坐標(biāo)。

    圖1 Layerwise板的位移場(chǎng) Fig.1 Displacement field of the Layerwise plate

    按照Layerwise理論和計(jì)算假設(shè)(1),復(fù)合材料層合結(jié)構(gòu)中第k(k≥2)層的位移場(chǎng){u}k為

    {u}k=

    (1a)

    第k=1層的位移場(chǎng)與通常的Mindlin板相同,即

    (1b)

    (2)

    這樣,第k層的位移場(chǎng)可記為

    {u}k=[N]kj5i0abt0b

    (3)

    式中,矩陣[N]k為

    [N]k=

    1.2應(yīng)變場(chǎng)

    (5)

    其中,

    式中,M、C、B和S分別表示膜(membrane)、耦合(coupling)、彎曲(bending)和橫向剪切(transverse shear)應(yīng)變項(xiàng)。

    因此,各層的應(yīng)變場(chǎng)可表示為

    (7)

    (8)

    (9)

    即重新排列后的應(yīng)變可表示為

    {ε}k=[B]kj5i0abt0b

    8×18×(2n+3)(2n+3)×1

    (10)

    1.3應(yīng)力場(chǎng)

    仍按式(8)的排列形式,根據(jù)應(yīng)力-應(yīng)變物理方程和假設(shè)(2),應(yīng)力場(chǎng){σ}k為

    {σ}k=[D]k{ε}k

    8×18×88×1

    (11)

    式中,[D]k為物理矩陣。

    對(duì)于各向同性不可壓粘彈性阻尼材料,通常采用復(fù)模量模型描述其本構(gòu)關(guān)系,即材料彈性模量為

    E*=E(1+jβ)

    (12)

    1.4應(yīng)變能和動(dòng)能

    層合結(jié)構(gòu)的應(yīng)變能表示為各層應(yīng)變能的疊加,利用式(10)和式(11),在待求域Ωk上積分可得

    (13)

    層合結(jié)構(gòu)的動(dòng)能可表示為各層動(dòng)能的疊加,即

    (14)

    (15)

    1.5四節(jié)點(diǎn)四邊形等參數(shù)板單元

    由式(3)可知,廣義位移場(chǎng)為二維坐標(biāo)x,y的連續(xù)函數(shù)。為求得廣義位移場(chǎng),按照有限元法基本思想,將求解域離散化為四邊形等參數(shù)單元組成的有限元計(jì)算模型,每個(gè)單元有四個(gè)節(jié)點(diǎn),如圖2所示,其位移插值函數(shù)和等參坐標(biāo)變換函數(shù)均為[10]

    (i=1,...,4)

    (16)

    將每個(gè)單元節(jié)點(diǎn)處參考平面上3個(gè)平動(dòng)位移和各層法線繞坐標(biāo)軸的轉(zhuǎn)角取為節(jié)點(diǎn)位移,則每個(gè)單元節(jié)點(diǎn)自由度的數(shù)目為4×(2n+3)=8n+12。可見,單元層數(shù)分的越多,單元矩陣的維數(shù)就越大。

    圖2 四節(jié)點(diǎn)等參單元局部與參數(shù)坐標(biāo)系 Fig. 2 Local and natural coordinate system of 4 node quadrilateral element

    (17)

    進(jìn)一步,將式(17)代入式(13)和式(15),并在單元子域上計(jì)算,可得單元應(yīng)變能和動(dòng)能的表達(dá)式,再應(yīng)用Hamilton原理,最終可得單元?jiǎng)偠染仃嚭唾|(zhì)量矩陣分別為

    (18)

    (19)

    其中,[B]k是單元第k層的幾何矩陣

    [B]k=[B]k[N]

    (20)

    [J]k是單元第k層的一致慣性矩陣

    [J]k=ρk[N]Tk[N]k

    (21)

    2系統(tǒng)運(yùn)動(dòng)方程

    式(18)、(19)定義的單元?jiǎng)偠染仃嚭唾|(zhì)量矩陣是定義在單元坐標(biāo)系下的,將其進(jìn)行坐標(biāo)變換至結(jié)構(gòu)整

    體坐標(biāo)系下,并組裝可得系統(tǒng)自由振動(dòng)的運(yùn)動(dòng)微分方程為

    (22)

    式中,[M]為總體質(zhì)量矩陣,[K]為總體剛度矩陣,{di}為結(jié)構(gòu)全局節(jié)點(diǎn)自由度向量。

    (23)

    (24)

    本文用Matlab編寫了相應(yīng)有限元求解程序,并采用eigs函數(shù)直接求解式(22)的復(fù)特征值問題[11]。

    3數(shù)值驗(yàn)證和結(jié)果分析

    3.1粘彈性阻尼夾層懸臂梁

    為驗(yàn)證本文的分析方法和程序,以文獻(xiàn)[12]給出的粘彈性阻尼夾層懸臂梁解析模型為對(duì)象,該梁平面尺寸為177.8 mm×25.4 mm,上、下彈性層及阻尼夾層的厚度及材料參數(shù)見表1所示。

    表1 粘彈性阻尼夾層梁材料參數(shù)

    表2 粘彈性阻尼夾層懸臂梁固有頻率和損耗因子

    采用本文的Layerwise板單元和分析方法,將夾層梁用四邊形板單元離散,長(zhǎng)度方向單元數(shù)14,寬度方向?yàn)?,厚度方向離散為3層。表2給出了針對(duì)該模型的六階梁理論解析解[12],文獻(xiàn)[13]采用三維實(shí)體單元建模的有限元解以及本文的分析結(jié)果。

    由表2可知,無論在粘彈性阻尼材料低損耗因子還是高損耗因子的情況下,基于Layerwise理論的板單元都能準(zhǔn)確描述阻尼材料的剪切變形和結(jié)構(gòu)的動(dòng)剛度特性,所給出的固有頻率和模態(tài)損耗因子的計(jì)算結(jié)果都與文獻(xiàn)中已公開結(jié)果符合很好,說明本文單元和分析方法可以用于分析粘彈性阻尼夾層結(jié)構(gòu)的固有頻率和阻尼特性。

    3.2共固化粘彈性復(fù)合材料層合圓柱殼

    分析文獻(xiàn)[14]中的共固化粘彈性復(fù)合材料層合圓柱殼的阻尼性能。圓柱殼體的幾何尺寸如圖3所示。殼體固支ABCD面,鋪層為[(02/902)/V/(902/ 02)],其中V表示粘彈性阻尼層。正交各向異性鋪層單層厚度為0.25 mm,計(jì)算時(shí)同時(shí)考慮復(fù)合材料自身的損耗能力,具體材料參數(shù)為:E1=119 GPa,E2=8.67 GPa,μ12=0.31,G12=G13=5.18 GPa,G23=3.9 GPa,ρ=1 570kg/m3,η11=0.118%,η22=0.620%,η12=η13=0.812%,η23=0.846%。

    粘彈性阻尼材料層厚度為0.25 mm,材料參數(shù)為:E=4.25 MPa,ν=0.49,ρ=968.8 kg/m3,β=1.56。

    圖3 共固化粘彈性復(fù)合材料層合圓柱殼幾何尺寸 Fig.3 Geometry and configuration of cylindrical composite blades with co-cured viscoelastic layer

    表3 共固化粘彈性復(fù)合材料圓柱殼固有頻率和損耗因子

    取Layerwise板單元網(wǎng)格密度為16×8,厚度方向離散為3層,其中粘彈性阻尼材料層為單獨(dú)層,上下層合板分別考慮為等效層。表3給出了無阻尼層和共固化粘彈性復(fù)合材料層合圓柱殼固有頻率和模態(tài)損耗因子的計(jì)算結(jié)果,并與文獻(xiàn)[14],以及傳統(tǒng)的采用混合單元法建模使用MSC.NASTRAN直接復(fù)特征值分析(SOL107)得到的結(jié)果進(jìn)行了對(duì)比。

    需要注意的是,圓柱殼為曲面殼體,采用混合單元法建模時(shí),等效板單元節(jié)點(diǎn)偏置定義需保證其沿當(dāng)?shù)厍娴姆ň€方向,而這對(duì)復(fù)雜曲面殼體而言非常繁瑣。采用layerwise板單元時(shí),只需將參考面按通常的平面單元離散,各層的位置和材料屬性轉(zhuǎn)化為物理參數(shù)輸入,從而大大簡(jiǎn)化了前處理的過程,這對(duì)復(fù)雜形狀結(jié)構(gòu)和嵌入多阻尼層的情況是極方便的。此外,通用有限元軟件中往往不支持考慮復(fù)合材料自身各向異性的損耗能力,而在Layerwise板單元,可將各層的彈性系數(shù)都考慮為復(fù)剛度形式,從而能夠有效的計(jì)算復(fù)合材料自身的損耗能力。

    由表3結(jié)果可知,只考慮復(fù)合材料自身損耗能力的情況下,本文計(jì)算的前四階模態(tài)的固有頻率和損耗因子的結(jié)果與文獻(xiàn)[14]都符合很好,證明了Layerwise平板單元模擬曲面殼體動(dòng)態(tài)特性的能力。共固化粘彈性阻尼材料層后,結(jié)構(gòu)各階模態(tài)固有頻率變化不大,但模態(tài)損耗因子卻增大了約1個(gè)數(shù)量級(jí),表明嵌入阻尼層能顯著增大復(fù)合材料層合結(jié)構(gòu)的阻尼性能。對(duì)比嵌入阻尼層后的損耗因子分析結(jié)果,本文的計(jì)算結(jié)果與SOL107解法獲得的結(jié)果較符合,但與文獻(xiàn)[14]的差別近1倍。這主要是因?yàn)槲墨I(xiàn)[14]是采用的模態(tài)法求解結(jié)構(gòu)阻尼特性,而在此過程中,需要用到無阻尼結(jié)構(gòu)的振型矩陣以進(jìn)行物理空間和模態(tài)空間的坐標(biāo)變換。在材料損耗因子較小的情況下,無阻尼結(jié)構(gòu)的實(shí)特征值解與有阻尼結(jié)構(gòu)的復(fù)特征值解差別較小,而在大損耗因子的情況下,這種近似則可能會(huì)導(dǎo)致較大的誤差。為進(jìn)一步驗(yàn)證上述結(jié)論,使用MSC.NASTRAN的模態(tài)復(fù)特征值解法(SOL110)分析算例問題,結(jié)果也列于表3中??梢?,使用SOL110求解序列獲得結(jié)果與文獻(xiàn)[14]吻合的很好,從而說明本文計(jì)算結(jié)果與文獻(xiàn)[14]的誤差主要是由于所采用的復(fù)特征值解法不同造成的。進(jìn)一步結(jié)合3.1節(jié)中與解析解的比對(duì)可知,模態(tài)解法通常僅適用于材料損耗因子較低的情況,在材料高損耗因子的情況下,直接法獲得的結(jié)果更為可信[15]。

    3.3共固化粘彈性復(fù)合材料加筋板

    加筋板是航空航天飛行器上廣泛采用的薄壁結(jié)構(gòu)形式??紤]一四邊簡(jiǎn)支的復(fù)合材料加筋板,面板尺寸為500 mm×450 mm,面板鋪層為[0/45/-45/90/V/90/-45/45/0]T,其中V表示粘彈性材料阻尼層,單層鋪層厚度為0.2 mm;沿長(zhǎng)度方向有兩對(duì)稱分布筋條,間距為150 mm,高度為20 mm,筋條全為0度鋪層,總厚度為2 mm。各向異性鋪層材料屬性與3.2中相同。分別計(jì)算阻尼層厚度tv=0,0.2 mm, 0.5 mm時(shí)加筋結(jié)構(gòu)的阻尼性能。取Layewise板單元模型網(wǎng)格密度為25×24,厚度方向分為3層。

    需要注意的是,對(duì)于復(fù)雜復(fù)合材料阻尼結(jié)構(gòu),混合單元法中采用三維實(shí)體單元模擬阻尼層,采用基于經(jīng)典層合板理論的普通板單元模擬彈性層以及筋條,兩者之間需定義多點(diǎn)約束方程進(jìn)行連接,前處理過程較為繁瑣;而采用Layerwise板單元時(shí),只涉及四邊形板單元,可按傳統(tǒng)的板單元建模方法構(gòu)造模型。而且,兩者模型形成的計(jì)算規(guī)模也有差別,阻尼層厚度較薄時(shí),為了減小實(shí)體單元的長(zhǎng)厚比(本例中取10:1),網(wǎng)格密度需取得很密,結(jié)構(gòu)矩陣維數(shù)很高,而Layerwise板單元在較疏的網(wǎng)格密度下就可取得滿意的精度。從面向優(yōu)化設(shè)計(jì)應(yīng)用的角度看,使用混合單元時(shí),若修改阻尼層厚度,需要改變實(shí)體單元節(jié)點(diǎn)坐標(biāo),導(dǎo)致有限元網(wǎng)格的重構(gòu);而使用Layerwise板單元時(shí),阻尼層的材料參數(shù)和各層厚度都是物理參量,對(duì)它們進(jìn)行修改時(shí)與有限元網(wǎng)格無關(guān),大大方便了優(yōu)化設(shè)計(jì)。

    表4 復(fù)合材料加筋阻尼板固有頻率和損耗因子

    表4中對(duì)比了本文的分析結(jié)果與采用MSC.NASTRAN混合單元法分析的結(jié)果。由表4可知,當(dāng)tv=0時(shí),本文的分析結(jié)果與NASTRAN的分析結(jié)果非常符合,并且本文的分析還能給出僅考慮正交各向異性鋪層損耗因子時(shí)的復(fù)合材料結(jié)構(gòu)的阻尼性能。當(dāng)tv=0.2 mm和0.5 mm時(shí),隨著阻尼層厚度的增大,各階模態(tài)的固有頻率下降,損耗因子顯著增大,兩類方法預(yù)測(cè)的阻尼值和變化規(guī)律都一致,再次證明了本文方法的有效性。

    4結(jié)論

    本文利用Layerwise離散層理論推導(dǎo)了可用于共固化粘彈性層合結(jié)構(gòu)阻尼性能分析的四節(jié)點(diǎn)四邊形等參數(shù)板單元,應(yīng)用該單元分析了典型粘彈阻尼結(jié)構(gòu)的固有頻率和模態(tài)損耗因子,并與文獻(xiàn)已公開結(jié)果和傳統(tǒng)分析方法進(jìn)行了對(duì)比驗(yàn)證。結(jié)果表明:

    (1)基于Layerwise理論構(gòu)造的復(fù)合材料層合板單元能夠正確描述阻尼層的高剪切變形,計(jì)算精度好且能夠考慮各向異性彈性層自身的損耗能力,可應(yīng)用于共固化粘彈性復(fù)合材料結(jié)構(gòu)性能的分析。

    (2)相比于傳統(tǒng)的混合單元法,應(yīng)用本文的層合板單元對(duì)復(fù)雜阻尼結(jié)構(gòu)進(jìn)行分析,能夠大大簡(jiǎn)化前處理建模過程,計(jì)算規(guī)模也較小。

    另外,采用層合板單元避免了傳統(tǒng)方法在鋪層順序或阻尼夾層厚度改變時(shí)需要重構(gòu)有限元往網(wǎng)格的問題,可方便的推廣于粘彈性阻尼復(fù)合材料結(jié)構(gòu)/阻尼的一體化設(shè)計(jì)。

    參考文獻(xiàn)

    [1]Finegan I C, Gibson R F. Recent research on enhancement of damping in polymer composites [J]. Composite Structures, 1999, 44: 89-98.

    [2]余啟勇,馬玉璞,郭萬濤,等. 結(jié)構(gòu)阻尼復(fù)合材料及其研究進(jìn)展[J]. 材料工程, 2007(S1): 253-258.

    YU Qi-yong,MA Yu-pu,GUO Wan-tao,et al. Structural damping composite and its research progress[J]. Journal of Materiale Engineering,2007, (S1):253-258.

    [3]Carrera E, Brischetto S. A survey with numerical assessment of classical and refined theories for the analysis of sandwich plates [J]. Applied Mechanics Reviews, 2009, 62:301-316.

    [4]張少輝,陳花玲. 共固化復(fù)合材料粘彈阻尼結(jié)構(gòu)的損耗因子研究[J]. 航空材料學(xué)報(bào),2005, 25(1): 53-57.

    ZHANG Shao-hui, CHEN Hua-ling. Damping analysis of cocured composites with interleaved viscoelastic layers [J]. Journal of Aeronautical Materials,2005, 25(1): 53-57.

    [5]徐超,李瑞杰,游少雄. 衛(wèi)星飛輪支架的共固化阻尼減振設(shè)計(jì)[J]. 宇航學(xué)報(bào), 2010, 31(3):907-911.

    XU Chao, LI Rui-jie, YOU Shao-xiong. Passive vibration control design with co-curing damping composite for a satellite flywheel bracket [J], Journal of Astronautics, 2010, 31(3):907-911.

    [6]林松,徐超,吳斌. 共固化粘彈性復(fù)合材料的結(jié)構(gòu)多目標(biāo)優(yōu)化設(shè)計(jì)[J]. 宇航學(xué)報(bào), 2010, 31(8): 1900-1905.

    LIN Song, XU Chao, WU Bin. Multi-objective evolutionary optimization design of composite structure embedded with cocured viscoelastic layers [J]. Journal of Astronautics, 2010, 31 (8):1900-1905.

    [7]Reddy J N. Mechanics oflaminated composite plates: theory and analysis. Boca Raton: CRC Press, 1997.

    [8]徐超,張醒,李瑞杰,等. 多粘彈性膠膜夾層約束阻尼梁損耗因子分析[J]. 宇航學(xué)報(bào), 2009, 30(2): 458-462.

    XU Chao, ZHANG Xing, LI Rui-jie, et al. Modal loss factor analysis for sandwich beams embedded with multiple viscoelastic damping thin layers[J]. Journal of Astronautics, 2009, 30(2): 458-462.

    [9]Moreira R A S, Rodrigues J D. A Layerwise model for thin soft core sandwich plates [J]. Computers and Structures, 2006, 84: 1256-1263.

    [10]Cook R D, et al. 有限元分析的概念與應(yīng)用[M]. 關(guān)正西,強(qiáng)洪夫譯.西安:西安交通大學(xué)出版社,2007:176-179.

    [11]MATLAB Version R2009a. The Math Works Inc: Natick, U.S.A., 2009.

    [12]Rao D K. Frequency and loss factors of sandwich beams under various boundary conditions [J]. Journal of Mechanical Engineering Science,1978, 20: 271-282.

    [13]Soni M L,Bogner F K. Finite element vibration analysis of damped structures [J]. AIAA Journal,1981, 20(5): 700-707.

    [14]Oh Il-Kwon. Dynamic characteristics of cylindrical hybrid panels containing viscoelastic layer based on Layerwise mechanics [J]. Composites Part B: Engineering,2007,38(2):159-171.

    [15]Torvik P J, Runyon B. Modifications to the method of modal strain energy for improved estimates of loss factors for damped structures [J]. Shock and Vibration, 2007, 14:339-353.

    《振動(dòng)與沖擊》入選中國(guó)精品科技期刊

    近日,國(guó)家科技部中國(guó)科學(xué)技術(shù)信息研究所公布第3屆中國(guó)精品科技期刊遴選結(jié)果和2014 年度“領(lǐng)跑者5000——中國(guó)精品科技期刊頂尖論文(F5000)”收錄情況?!墩駝?dòng)與沖擊》作為中國(guó)精品科技期刊頂尖學(xué)術(shù)論文(F5000)項(xiàng)目來源期刊,有20篇論文入選F5000。自2012年以來,《振動(dòng)與沖擊》共有50篇論文被評(píng)為F5000論文。

    中國(guó)精品科技期刊頂尖學(xué)術(shù)論文(F5000)是由中國(guó)科學(xué)技術(shù)信息研究所于2012年底啟動(dòng),每年評(píng)選一次。中國(guó)科學(xué)技術(shù)信息研究所從300多家中國(guó)精品科技期刊近5年發(fā)表的論文中評(píng)選出“領(lǐng)跑者5000——中國(guó)精品科技期刊頂尖學(xué)術(shù)論文”,以集中展示和交流我國(guó)的頂尖學(xué)術(shù)論文。通過“領(lǐng)跑者5000——中國(guó)精品科技期刊頂尖學(xué)術(shù)論文平臺(tái)(F5000)”與國(guó)際和國(guó)內(nèi)重要檢索系統(tǒng)鏈接,并利用WOK國(guó)際檢索平臺(tái)與SCI數(shù)據(jù)庫在同一平臺(tái)內(nèi)實(shí)現(xiàn)鏈接和國(guó)際引文檢索,向世界和科技界同行展示和推廣我國(guó)最重要的科研成果,擴(kuò)大中國(guó)學(xué)術(shù)期刊和出版機(jī)構(gòu)的國(guó)際影響。

    第一作者裴亞魯男,碩士生,1990年1月生

    通信作者張志誼男,博士生導(dǎo)師,1970年11月生

    郵箱:chychang@sjtu.edu.cn

    猜你喜歡
    粘彈性復(fù)合材料有限元
    二維粘彈性棒和板問題ADI有限差分法
    時(shí)變時(shí)滯粘彈性板方程的整體吸引子
    民機(jī)復(fù)合材料的適航鑒定
    復(fù)合材料無損檢測(cè)探討
    不可壓粘彈性流體的Leray-α-Oldroyd模型整體解的存在性
    磨削淬硬殘余應(yīng)力的有限元分析
    TiO2/ACF復(fù)合材料的制備及表征
    基于SolidWorks的吸嘴支撐臂有限元分析
    RGO/C3N4復(fù)合材料的制備及可見光催化性能
    環(huán)境溫度作用下瀝青路面熱粘彈性溫度應(yīng)力分析
    河南科技(2014年15期)2014-02-27 14:12:43
    www日本在线高清视频| 91av网站免费观看| 国产伦理片在线播放av一区| 午夜福利欧美成人| tube8黄色片| 亚洲专区中文字幕在线| 人人妻,人人澡人人爽秒播| 波多野结衣av一区二区av| 中文字幕制服av| 亚洲av成人不卡在线观看播放网| videos熟女内射| 亚洲少妇的诱惑av| 精品人妻1区二区| 99国产精品一区二区三区| 美女国产高潮福利片在线看| 国产精品久久久久久精品电影小说| 岛国在线观看网站| 不卡av一区二区三区| 99国产精品一区二区三区| 国产成+人综合+亚洲专区| av免费在线观看网站| 欧美成人午夜精品| 一区二区三区激情视频| 黄色视频在线播放观看不卡| 麻豆成人av在线观看| 在线 av 中文字幕| 99国产精品一区二区三区| 一二三四在线观看免费中文在| 黄色视频在线播放观看不卡| 免费在线观看影片大全网站| 女人高潮潮喷娇喘18禁视频| av欧美777| 国产欧美日韩一区二区三| 日本精品一区二区三区蜜桃| videosex国产| 欧美在线一区亚洲| 国产成人影院久久av| 国产黄色免费在线视频| 汤姆久久久久久久影院中文字幕| 亚洲天堂av无毛| 国产麻豆69| 一级片'在线观看视频| 午夜久久久在线观看| 久久ye,这里只有精品| 色视频在线一区二区三区| 一区二区三区国产精品乱码| 超碰成人久久| 一级a爱视频在线免费观看| 色老头精品视频在线观看| 午夜两性在线视频| 亚洲全国av大片| 日本欧美视频一区| 欧美黑人欧美精品刺激| videosex国产| 国产成人精品无人区| 最新的欧美精品一区二区| 亚洲天堂av无毛| 在线av久久热| 午夜成年电影在线免费观看| 黄色丝袜av网址大全| 脱女人内裤的视频| 我要看黄色一级片免费的| 老司机在亚洲福利影院| 欧美日韩精品网址| 亚洲专区中文字幕在线| 免费黄频网站在线观看国产| 国产不卡一卡二| 99国产精品免费福利视频| 国产精品自产拍在线观看55亚洲 | 9热在线视频观看99| 亚洲欧美激情在线| 欧美日韩亚洲高清精品| 91成人精品电影| 少妇被粗大的猛进出69影院| 国产欧美亚洲国产| 午夜福利在线免费观看网站| 51午夜福利影视在线观看| 超色免费av| 日本撒尿小便嘘嘘汇集6| 熟女少妇亚洲综合色aaa.| 怎么达到女性高潮| 亚洲精品av麻豆狂野| 欧美乱码精品一区二区三区| av不卡在线播放| 中文亚洲av片在线观看爽 | 老司机靠b影院| 国产xxxxx性猛交| 99re6热这里在线精品视频| 久久国产精品大桥未久av| 亚洲一码二码三码区别大吗| 香蕉国产在线看| 少妇的丰满在线观看| 搡老乐熟女国产| 精品卡一卡二卡四卡免费| a级毛片在线看网站| 亚洲色图 男人天堂 中文字幕| 一边摸一边抽搐一进一小说 | 99久久国产精品久久久| 免费av中文字幕在线| 精品少妇内射三级| 亚洲色图av天堂| 日韩中文字幕欧美一区二区| 91成人精品电影| 欧美黄色片欧美黄色片| 露出奶头的视频| 又大又爽又粗| 日本一区二区免费在线视频| 亚洲国产av新网站| 高清在线国产一区| 精品高清国产在线一区| 少妇 在线观看| 极品少妇高潮喷水抽搐| 国产一卡二卡三卡精品| 中文字幕人妻丝袜制服| 天天操日日干夜夜撸| 日本wwww免费看| 国产在线精品亚洲第一网站| a级毛片黄视频| 纯流量卡能插随身wifi吗| 老熟妇仑乱视频hdxx| 亚洲国产看品久久| 岛国毛片在线播放| 亚洲精华国产精华精| 成人国产av品久久久| 一边摸一边做爽爽视频免费| 丁香六月天网| 亚洲av美国av| 日本a在线网址| 嫁个100分男人电影在线观看| 精品熟女少妇八av免费久了| svipshipincom国产片| 精品少妇内射三级| 一二三四在线观看免费中文在| 97在线人人人人妻| 亚洲av美国av| 男女下面插进去视频免费观看| 最黄视频免费看| 丰满人妻熟妇乱又伦精品不卡| 超碰成人久久| 午夜精品国产一区二区电影| 欧美变态另类bdsm刘玥| 亚洲人成77777在线视频| 久久精品91无色码中文字幕| 免费不卡黄色视频| 国产精品电影一区二区三区 | 在线观看66精品国产| 嫩草影视91久久| 精品乱码久久久久久99久播| 国产福利在线免费观看视频| 亚洲国产av影院在线观看| 人人澡人人妻人| 9热在线视频观看99| 乱人伦中国视频| 精品人妻熟女毛片av久久网站| 视频在线观看一区二区三区| 中文字幕色久视频| av国产精品久久久久影院| 欧美日韩精品网址| 欧美av亚洲av综合av国产av| 免费在线观看完整版高清| 亚洲黑人精品在线| 十分钟在线观看高清视频www| 50天的宝宝边吃奶边哭怎么回事| 手机成人av网站| 91成人精品电影| 一二三四社区在线视频社区8| 最近最新免费中文字幕在线| 91麻豆av在线| 成人精品一区二区免费| 建设人人有责人人尽责人人享有的| 天天影视国产精品| 99精品欧美一区二区三区四区| 亚洲av国产av综合av卡| 国产精品亚洲av一区麻豆| 真人做人爱边吃奶动态| 老司机靠b影院| 美女高潮喷水抽搐中文字幕| 纵有疾风起免费观看全集完整版| e午夜精品久久久久久久| 99精品在免费线老司机午夜| 精品福利观看| 99国产极品粉嫩在线观看| 成年人免费黄色播放视频| 啪啪无遮挡十八禁网站| 十分钟在线观看高清视频www| 日日摸夜夜添夜夜添小说| 日本黄色视频三级网站网址 | 久久久精品免费免费高清| 久久国产精品大桥未久av| 日本一区二区免费在线视频| 日本av手机在线免费观看| 国产精品亚洲一级av第二区| 久久国产精品男人的天堂亚洲| videosex国产| 操出白浆在线播放| 欧美久久黑人一区二区| 水蜜桃什么品种好| 国产一区有黄有色的免费视频| 19禁男女啪啪无遮挡网站| 嫁个100分男人电影在线观看| 国产人伦9x9x在线观看| 老熟妇乱子伦视频在线观看| 丝袜喷水一区| 亚洲av电影在线进入| 国产欧美日韩一区二区精品| 国产精品一区二区在线不卡| 一二三四在线观看免费中文在| 美女午夜性视频免费| 黄色 视频免费看| 午夜福利欧美成人| 亚洲九九香蕉| 成人国产av品久久久| 大型黄色视频在线免费观看| 他把我摸到了高潮在线观看 | 亚洲欧美一区二区三区久久| av线在线观看网站| 国产主播在线观看一区二区| 99国产综合亚洲精品| 99久久国产精品久久久| 一级毛片女人18水好多| 动漫黄色视频在线观看| 久热这里只有精品99| 99国产精品一区二区蜜桃av | 国产高清videossex| 亚洲色图av天堂| 一区福利在线观看| 中文字幕另类日韩欧美亚洲嫩草| 国产麻豆69| 热re99久久国产66热| 午夜精品国产一区二区电影| 岛国毛片在线播放| 午夜精品久久久久久毛片777| 黄色怎么调成土黄色| 亚洲成人免费电影在线观看| 两个人免费观看高清视频| 亚洲一码二码三码区别大吗| 国产亚洲欧美精品永久| 中文字幕精品免费在线观看视频| 免费看a级黄色片| 久久中文字幕一级| 可以免费在线观看a视频的电影网站| 精品久久久精品久久久| 在线观看人妻少妇| 国产精品av久久久久免费| 操美女的视频在线观看| 王馨瑶露胸无遮挡在线观看| 大型黄色视频在线免费观看| 日本黄色日本黄色录像| 制服诱惑二区| 50天的宝宝边吃奶边哭怎么回事| 超碰97精品在线观看| 亚洲第一青青草原| 最新的欧美精品一区二区| 91大片在线观看| 日韩中文字幕欧美一区二区| 亚洲成av片中文字幕在线观看| 欧美精品一区二区免费开放| 国产高清国产精品国产三级| 午夜两性在线视频| 国产伦理片在线播放av一区| 亚洲人成电影观看| 高清在线国产一区| 侵犯人妻中文字幕一二三四区| 考比视频在线观看| 交换朋友夫妻互换小说| √禁漫天堂资源中文www| 一区二区三区乱码不卡18| 涩涩av久久男人的天堂| 国内毛片毛片毛片毛片毛片| 老司机午夜福利在线观看视频 | 无遮挡黄片免费观看| 大香蕉久久网| 精品国产一区二区久久| 精品高清国产在线一区| 国产日韩欧美亚洲二区| 亚洲性夜色夜夜综合| 亚洲免费av在线视频| 久久精品亚洲精品国产色婷小说| 美女午夜性视频免费| 满18在线观看网站| 国产高清国产精品国产三级| 久久久久久免费高清国产稀缺| 亚洲精品粉嫩美女一区| 操出白浆在线播放| 亚洲午夜理论影院| 国产在线一区二区三区精| 亚洲性夜色夜夜综合| 十分钟在线观看高清视频www| 高清欧美精品videossex| 男女无遮挡免费网站观看| 久久久久国内视频| 久久人妻av系列| 一级黄色大片毛片| 国产欧美日韩精品亚洲av| 无人区码免费观看不卡 | 国内毛片毛片毛片毛片毛片| 熟女少妇亚洲综合色aaa.| 在线观看www视频免费| 91精品三级在线观看| 欧美日韩一级在线毛片| 欧美日韩国产mv在线观看视频| 天天影视国产精品| 精品福利永久在线观看| 正在播放国产对白刺激| 不卡一级毛片| 日韩三级视频一区二区三区| 欧美日韩视频精品一区| 一区二区三区国产精品乱码| 精品亚洲成国产av| 亚洲国产欧美一区二区综合| 国产一区二区 视频在线| 老熟妇仑乱视频hdxx| 91av网站免费观看| 免费av中文字幕在线| 999精品在线视频| 99国产综合亚洲精品| 777米奇影视久久| 亚洲欧美日韩另类电影网站| 日韩三级视频一区二区三区| 一进一出抽搐动态| 欧美国产精品va在线观看不卡| 国产精品久久久人人做人人爽| 国产色视频综合| 大片免费播放器 马上看| 国产不卡一卡二| 超色免费av| 久久青草综合色| 欧美精品高潮呻吟av久久| 久久久久久久久久久久大奶| 我的亚洲天堂| 99久久人妻综合| 电影成人av| 人人澡人人妻人| 国产野战对白在线观看| 亚洲黑人精品在线| 久久久久精品国产欧美久久久| 午夜福利影视在线免费观看| 在线观看免费视频网站a站| 国产黄频视频在线观看| 99国产精品免费福利视频| 免费人妻精品一区二区三区视频| 久久 成人 亚洲| 丝袜喷水一区| 超色免费av| 99国产精品99久久久久| 久久狼人影院| 久久午夜亚洲精品久久| 热re99久久国产66热| 一区福利在线观看| 国产国语露脸激情在线看| 午夜激情av网站| 18禁观看日本| cao死你这个sao货| 久久久久久久久免费视频了| 新久久久久国产一级毛片| 丝袜人妻中文字幕| 一区二区三区国产精品乱码| 啦啦啦免费观看视频1| 国产精品成人在线| 亚洲欧美一区二区三区久久| 国产又爽黄色视频| 国产色视频综合| 无限看片的www在线观看| 国产不卡av网站在线观看| 日本wwww免费看| 黄色成人免费大全| 国产欧美亚洲国产| 成人三级做爰电影| 国产一区二区三区视频了| 免费观看人在逋| 午夜免费鲁丝| 亚洲欧美色中文字幕在线| 欧美精品一区二区免费开放| 国产亚洲精品一区二区www | 宅男免费午夜| 亚洲国产av影院在线观看| 黑人巨大精品欧美一区二区mp4| 国产精品自产拍在线观看55亚洲 | 99精品在免费线老司机午夜| 国产精品香港三级国产av潘金莲| 精品国产乱码久久久久久小说| 精品福利观看| 国产一区二区激情短视频| 一级片'在线观看视频| 男人舔女人的私密视频| 99国产精品一区二区蜜桃av | 国产免费视频播放在线视频| 欧美变态另类bdsm刘玥| 亚洲国产av影院在线观看| 麻豆乱淫一区二区| 国产av精品麻豆| 精品少妇久久久久久888优播| 日韩 欧美 亚洲 中文字幕| 日本一区二区免费在线视频| 在线十欧美十亚洲十日本专区| 高清视频免费观看一区二区| 成人三级做爰电影| 嫩草影视91久久| 电影成人av| 亚洲av片天天在线观看| 国产精品久久久人人做人人爽| 欧美乱妇无乱码| 在线十欧美十亚洲十日本专区| 国产免费现黄频在线看| svipshipincom国产片| 精品亚洲乱码少妇综合久久| 91国产中文字幕| 黄色视频在线播放观看不卡| 女人久久www免费人成看片| 伦理电影免费视频| 十八禁高潮呻吟视频| 91国产中文字幕| 五月天丁香电影| 亚洲av电影在线进入| 黑人欧美特级aaaaaa片| 亚洲第一av免费看| av福利片在线| 少妇 在线观看| 亚洲欧美一区二区三区黑人| 国精品久久久久久国模美| 在线永久观看黄色视频| 成人国产av品久久久| 少妇裸体淫交视频免费看高清 | 丰满少妇做爰视频| 一二三四社区在线视频社区8| 一级片'在线观看视频| 久久久久久久国产电影| 又大又爽又粗| 在线 av 中文字幕| 黑丝袜美女国产一区| 精品国内亚洲2022精品成人 | 热re99久久精品国产66热6| 日韩人妻精品一区2区三区| 99国产精品99久久久久| 黄片小视频在线播放| 国产激情久久老熟女| 久久久久久久国产电影| 在线观看免费日韩欧美大片| av不卡在线播放| 成在线人永久免费视频| 午夜免费鲁丝| 一区二区日韩欧美中文字幕| 成人国产av品久久久| 老鸭窝网址在线观看| av电影中文网址| 亚洲中文av在线| 国产人伦9x9x在线观看| 亚洲少妇的诱惑av| 成人黄色视频免费在线看| 侵犯人妻中文字幕一二三四区| 91av网站免费观看| 亚洲午夜理论影院| 日韩三级视频一区二区三区| 欧美大码av| 欧美黑人精品巨大| 亚洲成国产人片在线观看| 91字幕亚洲| 麻豆国产av国片精品| 国产精品久久久久久精品电影小说| 国产伦人伦偷精品视频| 精品国产国语对白av| 黄色视频不卡| 精品第一国产精品| 一级a爱视频在线免费观看| 一级毛片女人18水好多| 亚洲天堂av无毛| 成年人黄色毛片网站| 老熟妇乱子伦视频在线观看| 肉色欧美久久久久久久蜜桃| 黑人操中国人逼视频| 美女高潮到喷水免费观看| 中文字幕色久视频| 精品人妻1区二区| 免费在线观看黄色视频的| 99精品在免费线老司机午夜| 亚洲av美国av| 久久久久精品国产欧美久久久| 欧美变态另类bdsm刘玥| 2018国产大陆天天弄谢| 免费黄频网站在线观看国产| 女警被强在线播放| 国产无遮挡羞羞视频在线观看| 桃红色精品国产亚洲av| 变态另类成人亚洲欧美熟女 | 日韩有码中文字幕| 12—13女人毛片做爰片一| 国产色视频综合| 国产亚洲av高清不卡| 麻豆av在线久日| 曰老女人黄片| 91av网站免费观看| 国产精品久久久久成人av| 亚洲性夜色夜夜综合| 热re99久久国产66热| 国产精品.久久久| 成人影院久久| 宅男免费午夜| 久久久久久亚洲精品国产蜜桃av| 国产高清视频在线播放一区| 男人操女人黄网站| 日韩中文字幕欧美一区二区| 高清毛片免费观看视频网站 | 日本精品一区二区三区蜜桃| 久久影院123| 亚洲精品一二三| 亚洲精品国产色婷婷电影| 久久精品aⅴ一区二区三区四区| 久久久国产一区二区| 美女午夜性视频免费| 亚洲国产欧美在线一区| 搡老岳熟女国产| 精品午夜福利视频在线观看一区 | 亚洲一码二码三码区别大吗| 欧美人与性动交α欧美软件| 香蕉丝袜av| 超碰成人久久| 精品免费久久久久久久清纯 | 妹子高潮喷水视频| 在线观看免费日韩欧美大片| 日本黄色视频三级网站网址 | 精品国产一区二区久久| 成人国产一区最新在线观看| 操出白浆在线播放| 国产欧美日韩一区二区三区在线| 国产熟女午夜一区二区三区| 国产欧美日韩一区二区三区在线| 18禁黄网站禁片午夜丰满| 亚洲专区中文字幕在线| 国产欧美日韩综合在线一区二区| 不卡一级毛片| 午夜精品久久久久久毛片777| 国产欧美日韩一区二区三区在线| 亚洲色图 男人天堂 中文字幕| 美女视频免费永久观看网站| 制服诱惑二区| 亚洲欧美日韩另类电影网站| 国产黄色免费在线视频| 一区二区三区激情视频| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜视频精品福利| bbb黄色大片| 老熟女久久久| 国产精品亚洲av一区麻豆| 搡老岳熟女国产| 2018国产大陆天天弄谢| 久久精品aⅴ一区二区三区四区| svipshipincom国产片| 国产精品久久久人人做人人爽| 正在播放国产对白刺激| 高清视频免费观看一区二区| 黄色a级毛片大全视频| 啦啦啦 在线观看视频| 99国产精品99久久久久| 国产av精品麻豆| 国产精品熟女久久久久浪| 欧美在线一区亚洲| 精品少妇一区二区三区视频日本电影| 亚洲一区二区三区欧美精品| 天堂俺去俺来也www色官网| 午夜精品久久久久久毛片777| 天天躁夜夜躁狠狠躁躁| 精品久久久久久电影网| 涩涩av久久男人的天堂| 国产单亲对白刺激| avwww免费| 久久久精品94久久精品| 一本一本久久a久久精品综合妖精| 看免费av毛片| 黑丝袜美女国产一区| 在线十欧美十亚洲十日本专区| 大型av网站在线播放| 精品少妇一区二区三区视频日本电影| 日本av免费视频播放| 三上悠亚av全集在线观看| 日本vs欧美在线观看视频| 国产真人三级小视频在线观看| 成年女人毛片免费观看观看9 | 九色亚洲精品在线播放| 久久天躁狠狠躁夜夜2o2o| 丰满饥渴人妻一区二区三| 成年人黄色毛片网站| 精品一区二区三区视频在线观看免费 | 女人久久www免费人成看片| 亚洲av片天天在线观看| 精品国产一区二区三区久久久樱花| av不卡在线播放| 久久亚洲真实| 日本黄色视频三级网站网址 | 中文字幕人妻熟女乱码| 91成人精品电影| 午夜久久久在线观看| 亚洲成人国产一区在线观看| 看免费av毛片| 黄色视频,在线免费观看| 成年人午夜在线观看视频| 十八禁人妻一区二区| 久久久久久人人人人人| 日韩大片免费观看网站| 80岁老熟妇乱子伦牲交| 精品少妇久久久久久888优播| 90打野战视频偷拍视频| 不卡av一区二区三区| 最新的欧美精品一区二区| 国产高清视频在线播放一区| 国产97色在线日韩免费| 丁香欧美五月| 黑人欧美特级aaaaaa片| 成人免费观看视频高清| 国产精品久久久久成人av| 一级,二级,三级黄色视频| 欧美 日韩 精品 国产| 18禁裸乳无遮挡动漫免费视频| 国产av一区二区精品久久| 怎么达到女性高潮| 母亲3免费完整高清在线观看|