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

    具有約束層阻尼旋轉(zhuǎn)復合材料軸的動態(tài)特性研究

    2017-04-21 00:50:34時玉艷任勇生張玉環(huán)
    振動與沖擊 2017年7期
    關鍵詞:鋪層阻尼比固有頻率

    時玉艷, 任勇生, 張玉環(huán)

    (山東科技大學 機械電子工程學院,山東 青島 266590)

    具有約束層阻尼旋轉(zhuǎn)復合材料軸的動態(tài)特性研究

    時玉艷, 任勇生, 張玉環(huán)

    (山東科技大學 機械電子工程學院,山東 青島 266590)

    復合材料驅(qū)動軸在高速旋轉(zhuǎn)下的振動抑制,是先進直升機和汽車傳動系統(tǒng)結(jié)構(gòu)設計需要考慮的重要問題。約束層阻尼技術是工程領域內(nèi)普遍采用的實用有效的阻尼減振設計方法;然而通過采用約束層阻尼處置方式增加復合材料傳動軸的阻尼能力的研究報道,目前國內(nèi)外很少見到。從應力-應變關系、應變-位移關系出發(fā),推導出復合材料Timoshenko軸、約束層與黏彈性層的動能及勢能的數(shù)學表達式,采用 Hamilton原理建立具有約束層阻尼的旋轉(zhuǎn)復合材料軸的運動學方程。采用廣義Galerkin法對動力學方程進行了離散化,建立了廣義坐標表示的自由振動方程組,通過特征值求解得到系統(tǒng)的固有頻率和阻尼比;基于比例阻尼假設和四階Runge-Kutta數(shù)值積分方法求解上述方程,得到系統(tǒng)的自由振動響應曲線。通過數(shù)值分析揭示了約束層材料、黏彈性層材料、鋪層方式、長徑比和轉(zhuǎn)速對具有約束層阻尼的旋轉(zhuǎn)復合材料軸的固有振動特性以及自由振動響應特性的影響。

    被動約束層阻尼;復合材料軸;轉(zhuǎn)子系統(tǒng);阻尼自由振動

    驅(qū)動軸是用于直升機和汽車動力傳輸系統(tǒng)中的一類重要的結(jié)構(gòu)部件。與金屬材料軸相比,復合材料由于具有重量輕、比強度和比剛度高等優(yōu)點,因此,正在逐步取代金屬材料,在現(xiàn)代航空和汽車驅(qū)動軸結(jié)構(gòu)設計中獲得越來越廣泛的應用[1-5]。高性能飛機和汽車的研制對復合材料驅(qū)動軸在高速和超高速跨臨界旋轉(zhuǎn)狀態(tài)下的振穩(wěn)定性和工作可靠性,提出了更高的要求,因此,采用實用有效的方法對旋轉(zhuǎn)復合材料驅(qū)動軸實現(xiàn)振動抑制,就顯得十分重要。

    被動約束層阻尼是通過利用粘貼于受控結(jié)構(gòu)與約束層之間的高阻尼黏彈性材料耗散結(jié)構(gòu)的振動能量,從而達到減振的目的。約束層阻尼的概念自20世紀50年代提出以來,已經(jīng)發(fā)展成為一個成熟的結(jié)構(gòu)振動控制技術。采用被動約束層阻尼處置能夠在很寬的溫度和頻率范圍內(nèi)獲得較高的阻尼比,具有設計簡單、減振性能可靠等優(yōu)點,是工程界普遍認可的結(jié)構(gòu)振動控制措施之一[6]。

    近年來,約束層阻尼結(jié)構(gòu)的研究成果主要集中在含約束層阻尼的非旋轉(zhuǎn)各向同性材料梁、板和殼類結(jié)構(gòu)的阻尼特性研究方面[7-13],而對于旋轉(zhuǎn)軸,尤其是針對旋轉(zhuǎn)各向異性復合材料軸的研究報道較少。NAPOLITANO等[14]等對具有約束層阻尼和拉-扭耦合的復合材料旋轉(zhuǎn)軸進行了阻尼性能研究。他們通過有限元仿真和實驗分析,發(fā)現(xiàn)扭轉(zhuǎn)阻尼可以通過設計來實現(xiàn)。VENKATACHALAM等[15]對具有不同約束層阻尼的軸-盤系統(tǒng)進行了數(shù)值分析和試驗研究,采用基于板殼理論的半解析有限元法,分析了該系統(tǒng)的固有頻率和損耗因子。GHONEIM等[16]基于動力平衡方法建立具有部分約束層的Timoshenko復合材料軸的運動方程,并結(jié)合有限元法和假設振型法,對軸的阻尼能力進行研究。

    REN等[17]從復合材料單層的三維應力應變關系和Timoshenko梁理論出發(fā),在導出復合材料夾層軸的動能和勢能表達式的基礎上,采用Hamilton原理建立了具有約束層阻尼的旋轉(zhuǎn)復合材料軸的拉-彎-扭耦合運動方程,采用廣義Galerkin法對運動方程進行數(shù)值求解,通過對約束層阻尼復合材料軸的自由振動方程進行特征值分析,得到約束層阻尼復合材料軸的前三階固有頻率和阻尼比,研究了約束層厚度、阻尼層厚度的影響。

    本文在REN等工作的基礎上,進一步考察約束層材料、纖維鋪層方式、長徑比和黏彈性材料剪切模量,對旋轉(zhuǎn)復合材料軸的固有頻率及阻尼特性的影響?;诒壤枘峒俣?建立了具有約束層阻尼的復合材料軸自由振動時間響應分析模型,采用四階龍格-庫塔法進行數(shù)值積分,研究了不同參數(shù)對軸的自由振動響應特性的影響。

    1 理論建模

    如圖1所示,長度為L的復合材料軸,是由復合材料軸、黏彈性層和約束層組成。其中r為復合材料軸的任意層的半徑,rs是軸的內(nèi)徑。hs,hv和hc分別為軸的厚度、黏彈性層的厚度和約束層的厚度。

    為了建立具有約束層阻尼的旋轉(zhuǎn)復合材料軸的動力學方程,利用Hamilton原理

    (1)

    式中,δT=δTs+δTc+δTv,δU=δUs+δUc+δUv,下標“s”、“c”和“v”分別為軸、約束層和黏彈性層。

    復合材料軸、約束層和黏彈性層的動能分別表示為

    (a)

    (b)圖1 具有約束層阻尼的復合材料軸Fig. 1 Composite shaft with a constrained layer damping

    (2)

    (3)

    (4)

    復合材料軸、約束層及黏彈性層的勢能分別為

    (5)

    (6)

    (7)

    在柱坐標系 (x,r,α)下,復合材料軸和約束層的任意一層的本構(gòu)關系為

    (8)

    根據(jù)梁理論,σr,σα和ταx可以不予考慮。但在高速旋轉(zhuǎn)狀態(tài)下,應力σr,σα一般不為0,應當予以保留,但σr和σα的具體形式暫不確定,可以假設它們?yōu)檩S對稱。由式(8)的第一、第二和第六個方程可知,εr,εα和γrα能夠用εx,γαx,γxr,σr和σα表示。

    引入Timoshenko剪切修正系數(shù)κ[18],得到表達式

    (9)

    根據(jù)Timoshenko梁理論,復合材料軸和約束層內(nèi)任意一點的應變表示為

    (10)

    黏彈性層的本構(gòu)關系為

    (11)

    式中,Gv為黏彈性層的復數(shù)剪切模量。

    (12)

    將式(9)和式(10)代入式(5)和式(6),將式(11)和式(12)代入式(7),分別進行變分并代入式(1),得位移表示的彎曲-橫向剪切耦合運動方程如下

    (13)

    (14)

    (15)

    (16)

    (17)

    (18)

    設軸的兩端具有簡支邊界條件

    uy(0)=uz(0)=0,uy(L)=uz(L)=0

    (19)

    (20)

    為了得到復合材料軸的近似解,設位移和轉(zhuǎn)角的表達式為

    (21)

    (22)

    (23)

    式中,[M]、[G]和[K]分別為質(zhì)量矩陣、陀螺矩陣和剛度矩陣。

    振動方程式(23)能夠進一步化為廣義特征值問題

    Ay=λy

    (24)

    其中,

    (25)

    對振動系統(tǒng)式(24)進行特征值問題求解,可以得到系統(tǒng)的各階固有頻率和阻尼比。

    為了簡化對自由振動響應問題的分析過程,具有約束層阻尼的旋轉(zhuǎn)復合材料軸的結(jié)構(gòu)阻尼對振動響應的影響,基于Raleigh假設,通過引入如下比例阻尼矩陣進行分析

    [C]=a0[M]+a1[K]

    (26)

    式中:[M]和[K]分別為軸的模態(tài)質(zhì)量矩陣和模態(tài)剛度矩陣;a0,a1為比例系數(shù),由式(27)確定

    (27)

    式中:ω1和ω2為軸的前兩階固有頻率;η1和η2為相應前兩階阻尼比,這些固有頻率和阻尼比參數(shù)可以通過自由振動分析得到。

    將比例阻尼矩陣式(26)與運動方程式(23)中的陀螺矩陣[G]相疊加,得到系統(tǒng)的阻尼矩陣,對此阻尼自由運動方程進行數(shù)值積分,可獲得系統(tǒng)的自由振動時間響應曲線。

    2 數(shù)值結(jié)果與分析

    以具有約束層阻尼的旋轉(zhuǎn)復合材料軸為研究對象,軸的材料為石墨環(huán)氧樹脂(簡稱石墨/環(huán)氧),材料常數(shù)為:E11=139 GPa,E22=11 GPa,G12=G13=6.05 GPa,G23=3.78 GPa,v12=0.313,ρs=1 578 kg/m3;軸的基本幾何尺寸為:軸的截面平均半徑50 mm,總厚度為4 mm,約束層總厚度為2 mm。黏彈性層厚度為2 mm。對振動方程式(24)采用特征值技術進行固有振動特性分析,可得到固有頻率和阻尼比。

    2.1 固有頻率與阻尼比

    2.1.1 約束層材料的影響

    軸長L=1 m,共10層等厚,鋪層方式[θ/-θ]5。約束層分別采用高模量碳纖維環(huán)氧樹脂(簡稱高模碳/環(huán)氧)、石墨/環(huán)氧、鋼和鋁材料。其中,高模碳/環(huán)氧和石墨/環(huán)氧復合材料的鋪層方式[β/-β],鋪層角分別取β=0°、β=45°、β=90°,高模碳/環(huán)氧的材料常數(shù)為:E11=449 GPa,E22=9.3 GPa,G12=G13=3.24 GPa,G23=3.31 GPa,ν12=0.31,ν23=0.39,ρ=1 732 kg/m3;鋼材料常數(shù)為:E=207 GPa,G=79.6 GPa,ν=0.3,ρ=7 680 kg/m3;鋁材料常數(shù)為:E=69 GPa,G=26 GPa,ν=0.33,ρ=2 800 kg/m3;黏彈性材料Gv=69(1+i)MPa,ρv=1 575 kg/m3。

    圖2~圖4表示轉(zhuǎn)速Ω=0 ,具有不同約束層材料,約束層的鋪層角分別取0°、45°和90°時,固有頻率和阻尼比隨軸的鋪層角的變化曲線。結(jié)果表明,隨著約束層復合材料鋪層角的增加,固有頻率和阻尼比表現(xiàn)出明顯下降的趨勢。當β=0°時,采用復合材料約束層軸相對于采用金屬材料約束層軸的固有頻率和阻尼比而言,總體上有較大的增加;當β=90°時,具有復合材料約束層軸的固有頻率和阻尼卻低于具有金屬材料約束層軸的固有頻率和阻尼。原因是復合材料約束層彈性主方向1上的彈性模量與金屬材料彈性模量相比,數(shù)值相當或者較大,而復合材料約束層彈性主方向2上的彈性模量明顯小于金屬材料彈性模量。由此可見,為了能夠最大限度地提高軸的阻尼能力,同時不會以犧牲軸的固有頻率為代價,宜采用復合材料特別是高模量復合材料作為約束層。此外,從圖2~圖4中還可知,復合材料軸的鋪層角θ對于固有頻率和阻尼比也會產(chǎn)生明顯的影響。

    圖5(a)~圖5(b)表示具有不同約束層材料的軸固有頻率和阻尼比隨轉(zhuǎn)速的變化曲線。結(jié)果表明,由于陀螺效應的存在,隨著轉(zhuǎn)速的增加,固有頻率和阻尼比曲線分叉為上下兩支,上支隨著轉(zhuǎn)速的增加而增加,下支隨著轉(zhuǎn)速的增加而減小。由圖5(a)可知,對于給定的轉(zhuǎn)速Ω,固有頻率從大到小排列相對應的約束層材料依次是高模碳/環(huán)氧、石墨/環(huán)氧、鋼和鋁;由圖5(b)可知,阻尼比從大到小排列相對應的約束層材料依次是高模碳/環(huán)氧、石墨/環(huán)氧、鋼和鋁。這與圖2(a)~圖2(b)是一致的。

    (a)固有頻率

    (b)阻尼比圖2 不同約束層材料復合材料軸固有頻率和阻尼比隨鋪層角的變化曲線( β=0°) Fig. 2 The natural frequency and damping ratio versus fiber angle θ for various material of constraining layer( β=0°)

    (a)固有頻率

    (b)阻尼比圖3 不同約束層材料復合材料軸固有頻率和阻尼比隨鋪層角的變化曲線( β=45°) Fig. 3 The natural frequency and damping ratio versus fiber angle θ for various material of constraining layer( β=45°)

    (a)固有頻率

    (b)阻尼比圖4 不同約束層材料復合材料軸固有頻率和阻尼比隨鋪層角變化曲線( β=90°) Fig. 4 The natural frequency and damping ratio versus fiber angle θ for various material of constraining layer( β=90°)

    (a)固有頻率

    (b)阻尼比圖5 不同約束層材料復合材料軸固有頻率和阻尼比隨轉(zhuǎn)速變化曲線( β=0°,θ=60°)Fig. 5 The natural frequency and damping ratio versus rotating speed Ω for various material of constraining layer( β=0°,θ=60°)

    2.1.2 軸的鋪層方式的影響

    軸長L=1 m,共10層等厚,分別采用[θ/-θ]5和[θ]10兩種鋪層方式;約束層的材料為高模碳/環(huán)氧, 共2層等厚,鋪層方式[β/-β],β=0°;黏彈性材料剪切模量為Gv=69(1+i)MPa,ρv=1 575 kg/m3。

    圖6(a)~圖6(b)表示轉(zhuǎn)速Ω=0時,在具有不同鋪層方式的情況下,固有頻率和阻尼比隨鋪層角的變化曲線。結(jié)果表明,與鋪層方式[θ/-θ]5對應的固有頻率,高于與鋪層方式[θ]10對應的固有頻率,對于阻尼比而言,其值的大小與鋪層方式的關系與上述結(jié)論相反。當鋪層角大于65°時,兩種鋪層方式下的固有頻率和阻尼比曲線基本重合。

    (a)固有頻率

    (b)阻尼比圖6 不同鋪層方式下復合材料軸固有頻率和阻尼比隨鋪層角變化曲線( β=0°)Fig. 6 The natural frequency and damping ratio versus fiber angle θ for various stacking sequences( β=0°)

    圖7(a)~圖7(b)表示鋪層角θ=30°時,對應于不同鋪層方式的固有頻率和阻尼比隨轉(zhuǎn)速的變化曲線。結(jié)果表明,固有頻率和阻尼比大小與鋪層方式的對應關系,分別與圖6(a)和圖6(b)是一致的。

    2.1.3 長徑比的影響

    軸長分別采用L=0.5 m、L=1 m、L=1.5 m,L=2.5 m(即軸的長度與軸的截面半徑之比分別為10、20、30、50),共10層等厚,鋪層方式[θ/-θ]5;約束層的材料為高模碳/環(huán)氧, 共2層等厚,鋪層方式[β/-β],β=0°;黏彈性材料剪切模量為Gv=69(1+i)MPa,ρv=1 575 kg/m3。

    (a)固有頻率

    (b)阻尼比圖7 不同鋪層方式下復合材料軸固有頻率和阻尼比隨轉(zhuǎn)速變化曲線( β=0°)Fig. 7 The natural frequency and damping ratio versus rotating speed Ω for various stacking sequences( β=0°)

    圖8(a)~圖8(b)表示轉(zhuǎn)速Ω=0時,對應于不同長徑比固有頻率和阻尼比隨鋪層角的變化曲線。結(jié)果表明,隨著長徑比的增大,固有頻率和阻尼比隨著長徑比的增加逐漸減小,大長徑比下的固有頻率和阻尼比有彼此靠近的趨勢,在這種情況下,鋪層角的影響減弱。

    圖9(a)~圖9(b)表示鋪層角θ=60°時,不同長徑比軸的固有頻率和阻尼比隨轉(zhuǎn)速的變化曲線。結(jié)果表明,當長徑比增至30和50時,上下兩支的分叉現(xiàn)象相對不再明顯。此外,隨著長徑比的增大,復合材料軸的固有頻率值逐漸減小,且相互之間也在逐步靠近。這與圖8(a)~圖8(b)所得出的結(jié)論是相符合的。

    2.1.4 黏彈性層剪切模量的影響

    軸長L=1 m,共10層等厚,鋪層方式[θ/-θ]5;約束層的的材料為高模碳/環(huán)氧, 共2層等厚,鋪層方式[β/-β],β=0°。黏彈性材料的剪切模量分別為Gv=(69×107)*(1+i)Pa、Gv=(69×106)*(1+i)Pa、Gv=(69×105)*(1+i)Pa、Gv=(69×104)*(1+i)Pa,ρv=1 575 kg/m3。

    圖10(a)~圖10(b)表示轉(zhuǎn)速Ω=0時,具有不同黏彈性材料剪切模量軸的固有頻率和阻尼比隨鋪層角的變化曲線。結(jié)果表明,黏彈性材料剪切模量對復合材料軸的阻尼比有重要的影響,隨著黏彈性材料剪切模量的增大,復合材料軸的固有頻率也隨之增大;但阻尼比似乎并非

    (b)阻尼比圖8 不同長徑比復合材料軸固有頻率和阻尼比隨鋪層角變化曲線( β=0°)Fig.8 The natural frequency and damping ratio versus fiber angle θ for various ratio of the length and radius( β=0°)

    (a)固有頻率

    (b)阻尼比圖9 不同長徑比復合材料軸固有頻率和阻尼比隨轉(zhuǎn)速變化曲線( β=0°)Fig. 9 The natural frequency and damping ratio versus rotating speed Ω for various ratio of the length and radius( β=0°)

    隨著黏彈性材料剪切模量的增大而單調(diào)地增大,當黏彈性材料剪切模量增加超過一定的限度,如圖10(b) 中的Gv=(69×105)*(1+i)Pa,隨著黏彈性材料剪切模量的繼續(xù)增加,阻尼比反而會產(chǎn)生下降的趨勢,這與文獻[20]的研究發(fā)現(xiàn)是一致的。

    (a)固有頻率

    (b)阻尼比圖10 不同黏彈性材料剪切模量的復合材料軸固有頻率和阻尼比隨鋪層角θ變化曲線( β=0°)Fig. 10 The natural frequency and damping ratio versus fiber angle θ for various shear modulus of the viscoelastic material( β=0°)

    圖11(a)~圖11(b)表示鋪層角θ=60°時,具有不同黏彈性材料剪切模量軸的固有頻率和阻尼比隨轉(zhuǎn)速的變化曲線。結(jié)果表明,黏彈性材料剪切模量大小的變化,對固有頻率值影響不大,這與圖14(a)~圖14(b)的結(jié)果相一致。

    圖11(a)~圖11(b)表示鋪層角θ=60°時,具有不同黏彈性材料剪切模量軸的固有頻率和阻尼比隨轉(zhuǎn)速的變化曲線。結(jié)果表明,當剪切模量超過Gv=(69×105)*(1+i)Pa時,阻尼比-轉(zhuǎn)速曲線的上下兩支的分叉相對不再明顯。

    2.2 自由振動時間響應

    (a)固有頻率

    (b)阻尼比圖11 不同黏彈性材料剪切模量的復合材料軸固有頻率和阻尼比隨轉(zhuǎn)速Ω變化曲線( β=0°)Fig. 11 The natural frequency and damping ratio versus rotating speed Ω for various shear modulus of the viscoelastic material( β=0°)

    采用四階龍格-庫塔法對振動方程進行數(shù)值積分,得到復合材料軸的跨截面中點位置分別沿y、z方向的時間響應曲線,研究約束層材料、軸的鋪層方式、長徑比以及黏彈性層剪切模量對自由振動響應特性的影響。

    2.2.1 約束層材料的影響

    復合材料軸、約束層和黏彈性層的參數(shù)和鋪層與2.1.1中的取值相同,而軸的鋪層角θ=60°,旋轉(zhuǎn)速度Ω=2 000 r/min。

    圖12~圖15分別表示約束層采用高模碳/環(huán)氧、石墨/環(huán)氧、鋼和鋁時,軸中點(x=L/2)的彎曲位移時間響應uy(L/2,t)和uz(L/2,t)。結(jié)果表明,在軸和約束層的鋪層角一定的條件下,采用復合材料約束層,特別是高模量的碳纖維復合材料高模碳/環(huán)氧,能夠獲得比采用金屬材料約束層更為優(yōu)越的阻尼性能,這與圖2的結(jié)論是一致的。

    (a)y向時間響應 (b)z向時間響應圖12 具有高模碳/環(huán)氧約束層的復合材料軸中點位移時間響應曲線Fig.12 Displacement response curves with time the composite shaft using HMCE constraining layer

    (a)y向時間響應 (b)z向時間響應圖13 具有石墨/環(huán)氧約束層的復合材料軸中點位移時間響應曲線Fig.13 Displacement response curves with time the composite shaft using graphite-epoxy constraining layer

    (a)y向時間響應 (b)z向時間響應圖14 具有鋼約束層的復合材料軸中點位移時間響應曲線Fig.14 Displacement response curves with time the composite shaft using steel constraining layer

    (a)y向時間響應 (b)z向時間響應圖15 具有鋁約束層的復合材料軸中點位移時間響應曲線Fig.15 Displacement response curves with time the composite shaft using aluminum constraining layer

    2.2.2 軸的鋪層方式的影響

    復合材料軸、約束層和黏彈性層的參數(shù)和鋪層與“2.1.2”中的取值相同,旋轉(zhuǎn)速度Ω=2 000 r/min。圖16~圖17分別對應于采用兩種不同的鋪層方式, [40°/-40°]5和[40°]10,軸中點的彎曲位移時間響應。結(jié)果表明,復合材料軸在鋪層方式[40°]10下的響應曲線比鋪層方式[40°/-40°]5下的響應曲線,衰減速度要更快,這一結(jié)論與圖6的結(jié)論是相對應的。

    (a)y向時間響應 (b)z向時間響應圖16 具有鋪層方式[40°/-40°]5的復合材料軸中點位移時間響應曲線Fig.16 Displacement response curves with time the composite shaft with stacking sequence [40°/-40°]5

    (a)y向時間響應 (b)z向時間響應圖17 具有鋪層方式[40°]10的復合材料軸中點位移時間響應曲線Fig.17 Displacement response curves with time the composite shaft with stacking sequence [40°]10

    2.2.3 長徑比的影響

    復合材料軸、約束層和黏彈性層的參數(shù)和鋪層與2.1.3中的取值相同,而軸的鋪層角θ=60°,旋轉(zhuǎn)速度Ω=2 000 r/min。圖18~圖20分別表示軸的長徑比取10、20、50時,軸的彎曲位移時間響應。結(jié)果表明,短粗軸比細長軸的時間響應曲線衰減更快,因此,阻尼能力也更強,這也印證了圖8的結(jié)果。

    (a)y向時間響應 (b)z向時間響應圖18 長徑比為10的復合材料軸中點位移時間響應曲線Fig.18 Displacement response curves with time the composite shaft with the ratio of the length and radius 10

    (a)y向時間響應 (b)z向時間響應圖19 長徑比為20的復合材料軸中點位移時間響應曲線Fig.19 Displacement response curves with time the composite shaft with the ratio of the length and radius 20

    (a)y向時間響應 (b)z向時間響應圖20 長徑比為50的復合材料軸中點位移時間響應曲線Fig.20 Displacement response curves with time the composite shaft with the ratio of the length and radius 50

    2.2.4 黏彈性層剪切模量的影響

    復合材料軸、約束層和黏彈性層的參數(shù)和鋪層與2.1.4中的取值相同,而軸的鋪層角θ=60°,旋轉(zhuǎn)速度Ω=2 000 r/min。圖21~圖23分別表示黏彈性材料的剪切模量分別為Gv=69×103(1+i)Pa,Gv=69×104(1+i)Pa和Gv=69×105(1+i)Pa時,軸的彎曲位移時間響應。結(jié)果表明,在一定的范圍內(nèi),時間響應曲線的衰減速度隨著黏彈性材料的剪切模量的增加而增加,這與圖10的規(guī)律是一致的。

    (a)y向時間響應 (b)z向時間響應圖21 黏彈性層剪切模量Gv=69×103(1+i)Pa的復合材料軸中點位移時間響應曲線Fig.21 Displacement response curves with time the composite shaft with the shear modulus Gv=69×103(1+i)Pa

    (a)y向時間響應 (b)z向時間響應圖22 黏彈性層剪切模量Gv=69×104(1+i)Pa的復合材料軸中點位移時間響應曲線Fig.22 Displacement response curves with time the composite shaft with the shear modulus Gv=69×104(1+i)Pa

    (a)y向時間響應 (b)z向時間響應圖23 黏彈性層剪切模量Gv=69×105(1+i)Pa的復合材料軸中點位移時間響應曲線Fig.23 Displacement response curves with time the composite shaft with the shear modulus Gv=69×105(1+i)Pa

    最后,值得注意的是,在上述不同參數(shù)影響下的時間響應曲線均出現(xiàn)了拍振現(xiàn)象。原因是在旋轉(zhuǎn)狀態(tài)下,由于陀螺矩陣的存在,使得軸在y和z方向上的振動產(chǎn)生相互耦合,從而產(chǎn)生拍振。顯然,拍振頻率隨著轉(zhuǎn)速的增加而增加。

    3 結(jié) 論

    本文采用復合材料軸、約束層和黏彈性層的應力應變本構(gòu)關系、Timoshenko梁理論和Hamilton原理,建立具有約束層阻尼的復合材料軸的動力學方程?;趶V義Galerkin法得到系統(tǒng)的離散化方程。采用特征值技術對上述方程進行求解,分析系統(tǒng)的固有振動特性(固有頻率與阻尼比)。為了分析系統(tǒng)的自由振動時間響應特性,基于Raleigh假設,將比例阻尼矩陣引入系統(tǒng)的離散化方程,采用四階龍格-庫塔法進行數(shù)值積分。根據(jù)本文提出的模型與計算方法,針對不同的約束層材料、纖維鋪層方式、長徑比和黏彈性層剪切模量的復合材料軸,進行了固有頻率、阻尼比和時間響應的計算。主要研究結(jié)論如下:

    (1)復合材料約束層,特別是高模量碳纖維約束層與金屬材料約束層相比,能夠顯著增加約束層阻尼復合材料傳動軸的阻尼能力,同時不會以減低結(jié)構(gòu)本身的固有頻率為代價。

    (2)復合材料軸鋪層方式能夠?qū)S的固有頻率與阻尼比產(chǎn)生不容忽視的影響。

    (3)增大復合材料軸的長徑比,其固有頻率值和阻尼比均隨之減小。

    (4)復合材料軸的固有頻率隨著黏彈性材料剪切模量的增加而增加,在一定的范圍內(nèi),黏彈性材料剪切模量的增加將導致復合材料軸的阻尼比增加,當黏彈性材料剪切模量增加超出一定限度,復合材料軸的阻尼比隨著黏彈性材料剪切模量的增加而減小。

    (5)軸的旋轉(zhuǎn)效應對復合材料軸的振動特性能夠帶來明顯的影響,例如,轉(zhuǎn)速的增加會導致固有頻率和阻尼比出現(xiàn)分叉現(xiàn)象,旋轉(zhuǎn)狀態(tài)下的自由振動時間響應會發(fā)生拍振現(xiàn)象。

    本文的研究可進一步擴展,研究部分粘貼約束阻尼層的旋轉(zhuǎn)復合材料軸的動力學特性、考慮黏彈性材料頻變特性,以及考慮復合材料內(nèi)阻的旋轉(zhuǎn)復合材料軸的動力穩(wěn)定性問題。

    [ 1 ] HETHERINGTON E L, KRAUS R E, DARLOW M S. Demonstration of a super critical composite helicopter power transmission shaft [J]. Journal of American Helicopter Society, 1990, 35(1): 23-28.

    [ 2 ] DARLOW M S, CREONTE J. Optimal design of composite helicopter power transmission shafts with axially varying fiber lay-up[J]. Journal of the American Helicopter Society,1995,40(2):50-56.

    [ 3 ] LEE D G,KIM H S,KIM J W,et al. Design and manufacture of an automotive hybrid aluminum/composite drive shaft[J].Composite Structures,2004,63(1):87-99.

    [ 4 ] TALIB A R A,ALI A, BADIE M A,et al. Developing a hybrid,carbon/glass fiber-reinforced,epoxy composite automotive drive shaft[J].Materials and Design,2010,31(1):514-521.

    [ 5 ] BADIE M A, MAHDI E, HAMOUDA A M S.An investigation into hybrid carbon/glass fiber reinforced epoxy composite automotive drive shaft[J].Materials and Design, 2011,32(3):1485-1500.

    [ 6 ] 戴德沛.阻尼技術的工程應用[M]. 北京:清華大學出版社,1991.

    [ 7 ] ROSS D, UNGAR E E, KERWIN E M. Damping of plate flexural vibrations by means of viscoelastic laminate[M].Oxford:Pergamon Press, 1959.

    [ 8 ] DITARANTO R A. Theory of vibratory bending for elastic and viscoelastic layered finite-length beams[J]. Journal of Applied Mechanics,1965, 32(4):881-886.

    [ 9 ] MEAD D J, MARKUS S. The forced vibrations of a three-layer damped sandwich beam with arbitrary boundary conditions[J]. Journal of Sound & Vibration,1969, 10(2):163-175.

    [10] KERWIN E M. Damping of flexural waves by a constrained viscoelastic layer [J]. Journal of the Acoustical Society of America, 1959, 31(7): 952-962.

    [11] CHEN Y C, HUANG S C. An optimal placement of CLD treatment for vibration suppression of plates[J]. International Journal of Mechanical Sciences, 2002, 44(8):1801-1821.

    [12] DITARANTO R A. Theory of vibratory bending for elastic and viscoelastic layered finite-length beams[J]. Journal of Applied Mechanics,1965,32(4):881-886.

    [13] MEAD, D J, MARKUS S. The forced vibration of a three-layer damping sandwich beam with arbitrary boundary conditions [J]. Journal of Sound and Vibration, 1969, 10(2): 163-175.

    [14] NAPOLITANO K L, GRIPPO W, KOSMATKA J B, et al. A comparison of two cocured damped composite torsion shafts [J]. Composite Structures, 1998, 43(2): 115-125.

    [15] VENKATACHALAM R, PRABU S B. Vibration characteristics of orthotropic shaft-disk system with different constrained damping layers: experimental and numerical study [J].International Journal of Advanced Manufacturing Technology, 2013, 65: 601-610.

    [16] GHONEIM H, LAWRIE D J. Analysis of the flexural vibration of a composite drive shaft with partial cylindrical constrained layer damping treatment[J].Journal of Vibration and Control,2006,12(1):25-55.

    [17] REN Yongsheng,ZHANG Yuhuan. Free vibration and damping of rotating composite shaft with a constrained layer damping[J]. Shock and Vibration, 2016(4):1-20

    [18] DHARMARAJAN S, MCCUTCHEN H. Shear coefficients for orthotropic beams[J].Journal of Composite Materials,1973,7(4):530-535.

    [19] KIM W, ARGENTO A, SCOTT R A. Free vibration of a rotating tapered composite Timoshenko shaft[J]. Journal of Sound and Vibration,1999,226(1):125-147.

    [20] JIN G Y, YANG C M, LIU Z G, et al. A unified method for the vibration and damping analysis of constrained layer damping cylindrical shells with arbitrary boundary conditions[J]. Composite Structures,2015, 130:124-142.

    Vibration characteristics of a rotating composite shaft with constrained layer damping

    SHI Yuyan, REN Yongsheng, ZHANG Yuhuan

    (College of Mechanical and Electronic Engineering, Shandong University of Science and Technology, Qindao 266590, China)

    The vibration suppression of a composite shaft operating at high speeds is an important issue in structural design of advance power transmission of helicopters and automobiles. Constrained layer damping technology is an effective method in engineering field and it is widely used for reducing structural vibration. However, the studies on strengthening damping of a composite shaft by using constrained layer damping treatment are relatively few. Here, a mechanical model was developed for a rotating composite shaft with passive constrained layer damping (PCLD). Based on the constitutive relation of materials and the strain-displacement relation, the kinetic energy and strain energy of the composite shaft, constrained layer and viscoelastic layer were derived. Hamilton’s principle was used to derive the motion equations of the shaft with PCLD. The motion equations were discretized by using the general Galerkin method. The natural frequencies and modal damping ratios were obtained by computing a matrix eigenvalue problem. Using the assumption of proportional damping and Runge-Kutta integration algorithm, the time history responses of the system’s damped free vibration were gained. The effects of constrained layer material, viscoelastic layer material, stacking sequences, ratio of length to diameter and rotating speed on the natural vibration characteristics and damped free vibration responses of the composite shaft with PCLD were analyzed using numerical simulation.

    passive constrained layer damping (PCLD); composite shaft; rotor system; damped free vibration

    國家自然科學基金(11272190)

    2016-06-29 修改稿收到日期:2016-08-31

    時玉艷 女,碩士生,1989年生

    任勇生 男,博士,教授,1956年生

    TB33

    A

    10.13465/j.cnki.jvs.2017.07.008

    猜你喜歡
    鋪層阻尼比固有頻率
    現(xiàn)場測定大型水輪發(fā)電機組軸系的固有頻率
    大電機技術(2021年2期)2021-07-21 07:28:38
    基于細觀結(jié)構(gòu)的原狀黃土動彈性模量和阻尼比試驗研究
    地震研究(2021年1期)2021-04-13 01:05:24
    黏滯阻尼器在時程分析下的附加有效阻尼比研究
    振動與沖擊(2019年4期)2019-02-22 02:33:34
    波形分析法求解公路橋梁阻尼比的探討
    上海公路(2018年3期)2018-03-21 05:55:48
    衛(wèi)星天線復合材料框架的鋪層優(yōu)化設計
    結(jié)構(gòu)構(gòu)件阻尼比對大跨度懸索橋地震響應的影響
    復合材料軸結(jié)構(gòu)力學性能預測及鋪層方案設計
    CPD鋪層自動創(chuàng)建技術
    FIBERSIM軟件在工字梁鋪層設計上的應用
    教練機(2015年2期)2015-04-03 06:51:48
    總溫總壓測頭模態(tài)振型變化規(guī)律研究
    а√天堂www在线а√下载| 国产综合懂色| 黄色配什么色好看| 久久人人精品亚洲av| a级毛片免费高清观看在线播放| 男人的好看免费观看在线视频| 极品教师在线视频| 少妇熟女aⅴ在线视频| 日韩精品青青久久久久久| 精品午夜福利在线看| 在线观看午夜福利视频| 乱码一卡2卡4卡精品| 亚洲经典国产精华液单 | 国产美女午夜福利| 精品国内亚洲2022精品成人| 12—13女人毛片做爰片一| 在线免费观看的www视频| 欧美在线一区亚洲| 别揉我奶头 嗯啊视频| 欧美+日韩+精品| 亚洲真实伦在线观看| 久久久久久国产a免费观看| 久久久色成人| 精品久久久久久久人妻蜜臀av| 又紧又爽又黄一区二区| 精品乱码久久久久久99久播| 成年人黄色毛片网站| 久久久久精品国产欧美久久久| 黄色配什么色好看| 久久99热这里只有精品18| 两个人视频免费观看高清| 97超级碰碰碰精品色视频在线观看| 国产伦精品一区二区三区四那| 亚洲一区二区三区色噜噜| 免费电影在线观看免费观看| 国产精品嫩草影院av在线观看 | 可以在线观看毛片的网站| 国产一区二区三区视频了| 国产私拍福利视频在线观看| 精品欧美国产一区二区三| 最近最新中文字幕大全电影3| 免费av不卡在线播放| 亚洲美女视频黄频| 日韩欧美精品免费久久 | 哪里可以看免费的av片| 美女 人体艺术 gogo| 日日夜夜操网爽| 亚洲一区二区三区色噜噜| 天堂动漫精品| 非洲黑人性xxxx精品又粗又长| 小蜜桃在线观看免费完整版高清| 亚洲自拍偷在线| 欧美性猛交╳xxx乱大交人| 国产中年淑女户外野战色| 麻豆一二三区av精品| 一区二区三区高清视频在线| 俄罗斯特黄特色一大片| 欧美丝袜亚洲另类 | 神马国产精品三级电影在线观看| 一个人免费在线观看电影| 国产精品亚洲av一区麻豆| 韩国av一区二区三区四区| 国产精品影院久久| 亚洲av第一区精品v没综合| 久久久国产成人精品二区| 国产真实乱freesex| 国产大屁股一区二区在线视频| 国产亚洲精品久久久久久毛片| 91麻豆精品激情在线观看国产| 岛国在线免费视频观看| av黄色大香蕉| 成人av一区二区三区在线看| 欧美一区二区亚洲| 国产久久久一区二区三区| 国产精品1区2区在线观看.| 国产伦精品一区二区三区视频9| 中亚洲国语对白在线视频| 欧美国产日韩亚洲一区| 日本黄大片高清| 国产精品一区二区三区四区久久| 1024手机看黄色片| 亚洲中文日韩欧美视频| 婷婷色综合大香蕉| 国模一区二区三区四区视频| 成人欧美大片| 午夜福利在线在线| 天美传媒精品一区二区| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 在线观看av片永久免费下载| 亚洲av美国av| 99国产精品一区二区蜜桃av| 级片在线观看| 国产免费一级a男人的天堂| 国产一区二区激情短视频| 久99久视频精品免费| 天堂av国产一区二区熟女人妻| 又粗又爽又猛毛片免费看| 成人特级av手机在线观看| 中文字幕精品亚洲无线码一区| 99久久99久久久精品蜜桃| 久久国产乱子伦精品免费另类| 国产黄色小视频在线观看| 亚洲av不卡在线观看| 国产白丝娇喘喷水9色精品| 亚洲片人在线观看| av女优亚洲男人天堂| 午夜免费男女啪啪视频观看 | 一区二区三区高清视频在线| 国产精品自产拍在线观看55亚洲| 欧美高清性xxxxhd video| 亚洲人成网站在线播放欧美日韩| 成人特级av手机在线观看| 99久久精品国产亚洲精品| 欧美+日韩+精品| av欧美777| 男插女下体视频免费在线播放| 精品国内亚洲2022精品成人| 97人妻精品一区二区三区麻豆| 国产高清有码在线观看视频| 丰满人妻熟妇乱又伦精品不卡| 国产一级毛片七仙女欲春2| 国产欧美日韩精品亚洲av| 亚洲人成电影免费在线| 欧美+日韩+精品| 精品国产三级普通话版| 亚洲人成电影免费在线| 在线十欧美十亚洲十日本专区| 国内揄拍国产精品人妻在线| 国产成人aa在线观看| 99久久99久久久精品蜜桃| 18禁黄网站禁片免费观看直播| 亚洲黑人精品在线| 可以在线观看毛片的网站| 欧美日本亚洲视频在线播放| 亚洲五月婷婷丁香| 亚洲一区二区三区不卡视频| 高清日韩中文字幕在线| 在线看三级毛片| 亚洲最大成人中文| 九色成人免费人妻av| 免费在线观看亚洲国产| 欧美最新免费一区二区三区 | 黄色丝袜av网址大全| 成年女人看的毛片在线观看| 免费看美女性在线毛片视频| 亚洲天堂国产精品一区在线| 国产精品综合久久久久久久免费| 免费人成视频x8x8入口观看| 天堂网av新在线| 男人舔奶头视频| 成人午夜高清在线视频| 真实男女啪啪啪动态图| 色在线成人网| 嫩草影院入口| 在线观看舔阴道视频| 午夜a级毛片| 日韩亚洲欧美综合| 成人高潮视频无遮挡免费网站| 国产私拍福利视频在线观看| 色视频www国产| 亚洲经典国产精华液单 | 内地一区二区视频在线| 久久精品91蜜桃| 高清日韩中文字幕在线| 俺也久久电影网| 看片在线看免费视频| 91av网一区二区| 午夜老司机福利剧场| 国产精品电影一区二区三区| 在线观看一区二区三区| 国产真实乱freesex| 内射极品少妇av片p| 亚洲综合色惰| 性欧美人与动物交配| 黄色日韩在线| 国产一区二区三区在线臀色熟女| 国产 一区 欧美 日韩| 91字幕亚洲| 亚洲av熟女| 精品一区二区三区av网在线观看| 美女大奶头视频| a在线观看视频网站| a级一级毛片免费在线观看| 精华霜和精华液先用哪个| www.www免费av| 99久久精品一区二区三区| 亚洲精品日韩av片在线观看| 国产真实乱freesex| 成人高潮视频无遮挡免费网站| 亚洲精品色激情综合| 免费在线观看成人毛片| 国产精品久久久久久久电影| 国产黄色小视频在线观看| 真实男女啪啪啪动态图| www日本黄色视频网| 亚洲最大成人手机在线| 日韩精品中文字幕看吧| 九九热线精品视视频播放| 成人亚洲精品av一区二区| 亚洲精品在线美女| 91在线精品国自产拍蜜月| 国产综合懂色| 国产不卡一卡二| 久久久国产成人免费| 熟妇人妻久久中文字幕3abv| 国产一区二区在线av高清观看| 成熟少妇高潮喷水视频| 国产三级中文精品| 久久久久国内视频| 美女被艹到高潮喷水动态| 国产一区二区在线观看日韩| 非洲黑人性xxxx精品又粗又长| h日本视频在线播放| 九色成人免费人妻av| 欧美另类亚洲清纯唯美| 国产av在哪里看| 午夜福利在线观看吧| 国产亚洲精品综合一区在线观看| 99热只有精品国产| 女人被狂操c到高潮| 日本黄大片高清| 最近中文字幕高清免费大全6 | 91在线观看av| 久久精品影院6| 成年女人毛片免费观看观看9| 国产一级毛片七仙女欲春2| 一夜夜www| 久久久久久久亚洲中文字幕 | 最近最新中文字幕大全电影3| 最好的美女福利视频网| 精品熟女少妇八av免费久了| 欧美区成人在线视频| 综合色av麻豆| 久久久久久久亚洲中文字幕 | 精品一区二区三区av网在线观看| eeuss影院久久| 久9热在线精品视频| 少妇熟女aⅴ在线视频| 天堂动漫精品| 亚洲最大成人手机在线| 国产精品美女特级片免费视频播放器| 亚洲综合色惰| 最近最新免费中文字幕在线| 麻豆av噜噜一区二区三区| 日日摸夜夜添夜夜添小说| 人妻夜夜爽99麻豆av| 黄片小视频在线播放| 真人一进一出gif抽搐免费| 日韩大尺度精品在线看网址| 日本免费一区二区三区高清不卡| 精品久久久久久久久av| x7x7x7水蜜桃| 在线观看舔阴道视频| 村上凉子中文字幕在线| 国产伦人伦偷精品视频| 人妻久久中文字幕网| 婷婷精品国产亚洲av| 亚洲av.av天堂| 国产精华一区二区三区| 亚洲av中文字字幕乱码综合| 嫩草影视91久久| 国产精品久久视频播放| 国产毛片a区久久久久| 亚洲国产欧美人成| 国产探花在线观看一区二区| 美女高潮喷水抽搐中文字幕| 国产精品永久免费网站| 精品久久久久久久久亚洲 | 欧美高清性xxxxhd video| 午夜免费成人在线视频| 51国产日韩欧美| 在线观看免费视频日本深夜| 91久久精品国产一区二区成人| 色综合婷婷激情| 精品午夜福利在线看| 久久伊人香网站| 欧美日本亚洲视频在线播放| 国产成人aa在线观看| 熟女电影av网| 中国美女看黄片| 搞女人的毛片| 欧美性猛交╳xxx乱大交人| 啦啦啦观看免费观看视频高清| 亚洲片人在线观看| 精品久久久久久成人av| 亚洲欧美日韩无卡精品| 99热6这里只有精品| 国产一级毛片七仙女欲春2| 欧美黑人巨大hd| 欧美激情国产日韩精品一区| 国产精品女同一区二区软件 | 国产爱豆传媒在线观看| 在现免费观看毛片| 亚洲中文字幕日韩| 嫁个100分男人电影在线观看| 亚洲国产精品合色在线| 精品免费久久久久久久清纯| 国产免费男女视频| 免费看美女性在线毛片视频| 午夜日韩欧美国产| 欧美日韩中文字幕国产精品一区二区三区| 国产精品电影一区二区三区| АⅤ资源中文在线天堂| 国产成人a区在线观看| 国产精品av视频在线免费观看| 欧美3d第一页| 精品人妻1区二区| 99在线视频只有这里精品首页| 国产视频一区二区在线看| 超碰av人人做人人爽久久| 日日摸夜夜添夜夜添小说| 国产成年人精品一区二区| 亚洲不卡免费看| 欧美日本亚洲视频在线播放| 亚洲人与动物交配视频| 日韩人妻高清精品专区| 一进一出抽搐动态| 12—13女人毛片做爰片一| 禁无遮挡网站| 999久久久精品免费观看国产| 国模一区二区三区四区视频| 国产老妇女一区| 夜夜爽天天搞| 老熟妇乱子伦视频在线观看| 岛国在线免费视频观看| 国产精品一区二区三区四区免费观看 | 久久国产精品影院| 国产精品久久视频播放| 国产老妇女一区| 亚洲五月天丁香| 成人国产综合亚洲| 欧美激情在线99| 国产探花极品一区二区| 69人妻影院| 欧美乱色亚洲激情| 性色av乱码一区二区三区2| 在线播放无遮挡| 九九热线精品视视频播放| 成人欧美大片| 丰满乱子伦码专区| 91在线观看av| 免费观看人在逋| 日韩欧美在线乱码| 欧美性感艳星| 国产午夜精品论理片| 舔av片在线| 在线国产一区二区在线| 18禁在线播放成人免费| 久久久久久久久久成人| 日日干狠狠操夜夜爽| 国产午夜精品论理片| 亚洲性夜色夜夜综合| 国产亚洲精品综合一区在线观看| 欧美性猛交╳xxx乱大交人| 午夜日韩欧美国产| 18禁在线播放成人免费| 好看av亚洲va欧美ⅴa在| 黄色视频,在线免费观看| 日韩欧美在线乱码| 一级毛片久久久久久久久女| 日韩国内少妇激情av| 男女之事视频高清在线观看| 激情在线观看视频在线高清| 中亚洲国语对白在线视频| 国产精品98久久久久久宅男小说| 精品久久久久久久末码| 又紧又爽又黄一区二区| 欧美乱妇无乱码| 熟女电影av网| 亚洲欧美日韩高清专用| 91字幕亚洲| 51国产日韩欧美| 男人舔奶头视频| 很黄的视频免费| 免费看a级黄色片| 国产69精品久久久久777片| 久久人人精品亚洲av| 女人被狂操c到高潮| 国产真实伦视频高清在线观看 | 日日干狠狠操夜夜爽| 亚洲最大成人中文| 亚洲五月天丁香| 一级黄色大片毛片| 欧美成狂野欧美在线观看| 亚洲国产高清在线一区二区三| 日韩高清综合在线| 午夜福利在线观看免费完整高清在 | 亚洲va日本ⅴa欧美va伊人久久| 天美传媒精品一区二区| 亚洲va日本ⅴa欧美va伊人久久| 日日干狠狠操夜夜爽| 亚洲欧美日韩东京热| 性色av乱码一区二区三区2| 国产精品乱码一区二三区的特点| 91久久精品国产一区二区成人| h日本视频在线播放| 国产精品久久视频播放| 永久网站在线| 最近在线观看免费完整版| 国产精品一区二区性色av| 午夜精品久久久久久毛片777| 日本免费一区二区三区高清不卡| 国产亚洲欧美98| 成人三级黄色视频| 色尼玛亚洲综合影院| 免费电影在线观看免费观看| 欧美潮喷喷水| 美女xxoo啪啪120秒动态图 | 久久久色成人| 亚洲 国产 在线| 99热这里只有是精品在线观看 | 观看免费一级毛片| 精品午夜福利在线看| 欧美潮喷喷水| 欧美一区二区亚洲| 精品福利观看| 男人和女人高潮做爰伦理| 欧美一区二区精品小视频在线| 可以在线观看的亚洲视频| 午夜福利在线在线| 色av中文字幕| 国产高清激情床上av| 中文在线观看免费www的网站| 五月玫瑰六月丁香| 在线国产一区二区在线| 一区二区三区四区激情视频 | 别揉我奶头~嗯~啊~动态视频| 大型黄色视频在线免费观看| 首页视频小说图片口味搜索| 九九久久精品国产亚洲av麻豆| av在线老鸭窝| 国产三级黄色录像| 特级一级黄色大片| 欧美成人一区二区免费高清观看| 12—13女人毛片做爰片一| 少妇人妻精品综合一区二区 | 18禁黄网站禁片午夜丰满| 国产精品一及| 国产精品一区二区免费欧美| 成人特级av手机在线观看| 国产精品亚洲美女久久久| 久久6这里有精品| 精品不卡国产一区二区三区| 亚洲欧美日韩高清专用| 自拍偷自拍亚洲精品老妇| 91麻豆av在线| 丰满人妻一区二区三区视频av| 欧美一级a爱片免费观看看| 色综合站精品国产| 欧美乱色亚洲激情| 成人无遮挡网站| 身体一侧抽搐| 夜夜夜夜夜久久久久| 亚洲成人中文字幕在线播放| 高清在线国产一区| 亚洲av成人av| 日日摸夜夜添夜夜添av毛片 | 久久精品国产亚洲av天美| 黄片小视频在线播放| a级一级毛片免费在线观看| 亚洲av一区综合| 亚洲av二区三区四区| 青草久久国产| 亚洲精品在线观看二区| 亚洲av不卡在线观看| 国产精品,欧美在线| 国产欧美日韩精品一区二区| 男女之事视频高清在线观看| 国产精品久久久久久久电影| 非洲黑人性xxxx精品又粗又长| 婷婷色综合大香蕉| 精品久久久久久久久久免费视频| 最近视频中文字幕2019在线8| 黄色一级大片看看| 在现免费观看毛片| 色在线成人网| 亚洲天堂国产精品一区在线| a级毛片免费高清观看在线播放| 亚洲va日本ⅴa欧美va伊人久久| 日韩欧美精品免费久久 | 成熟少妇高潮喷水视频| 亚洲自偷自拍三级| 国产熟女xx| 亚洲精华国产精华精| 成人av一区二区三区在线看| 老司机午夜福利在线观看视频| 他把我摸到了高潮在线观看| 少妇的逼水好多| 一区二区三区四区激情视频 | 久久久久性生活片| 久久亚洲精品不卡| 免费在线观看日本一区| 一夜夜www| 全区人妻精品视频| 99热这里只有是精品在线观看 | 99久久精品一区二区三区| 性色av乱码一区二区三区2| 女同久久另类99精品国产91| 欧美xxxx性猛交bbbb| 色噜噜av男人的天堂激情| 最近最新中文字幕大全电影3| 综合色av麻豆| 黄色日韩在线| 国产麻豆成人av免费视频| 婷婷色综合大香蕉| 怎么达到女性高潮| 观看美女的网站| 日本黄色视频三级网站网址| 麻豆av噜噜一区二区三区| 他把我摸到了高潮在线观看| 日韩欧美在线乱码| 国模一区二区三区四区视频| 国产精品久久久久久亚洲av鲁大| 日本三级黄在线观看| 亚洲精品粉嫩美女一区| 黄色配什么色好看| 国产真实乱freesex| 婷婷色综合大香蕉| 欧美最黄视频在线播放免费| 成人午夜高清在线视频| 国产高清视频在线播放一区| 欧美日韩亚洲国产一区二区在线观看| 长腿黑丝高跟| 欧美激情久久久久久爽电影| 男人和女人高潮做爰伦理| 日韩欧美精品v在线| 国产色爽女视频免费观看| 99国产综合亚洲精品| 欧美区成人在线视频| aaaaa片日本免费| 久久久成人免费电影| 中文字幕av在线有码专区| 久久亚洲真实| 久久精品人妻少妇| 99久久成人亚洲精品观看| or卡值多少钱| 99视频精品全部免费 在线| aaaaa片日本免费| 欧美午夜高清在线| 欧美日本视频| 免费观看的影片在线观看| 日本与韩国留学比较| 真实男女啪啪啪动态图| 欧美日韩乱码在线| 人妻夜夜爽99麻豆av| 亚洲av不卡在线观看| 少妇的逼好多水| 欧美激情久久久久久爽电影| 国产精品美女特级片免费视频播放器| 搡女人真爽免费视频火全软件 | 狂野欧美白嫩少妇大欣赏| 搡女人真爽免费视频火全软件 | 99久久精品国产亚洲精品| 我的女老师完整版在线观看| 国产伦一二天堂av在线观看| 一本一本综合久久| 欧美bdsm另类| 中文字幕高清在线视频| 色播亚洲综合网| 97碰自拍视频| 国产大屁股一区二区在线视频| 欧美性感艳星| 精品一区二区三区人妻视频| 国产精品久久久久久久久免 | 最近最新中文字幕大全电影3| 深爱激情五月婷婷| 永久网站在线| 欧美另类亚洲清纯唯美| 久久久久久大精品| 一级黄色大片毛片| 性欧美人与动物交配| 长腿黑丝高跟| 国语自产精品视频在线第100页| 非洲黑人性xxxx精品又粗又长| 亚洲av.av天堂| 国产精品99久久久久久久久| 亚洲欧美日韩高清专用| 在线观看一区二区三区| 久久久久久久久中文| 性色avwww在线观看| 深夜a级毛片| 国产av麻豆久久久久久久| 欧美xxxx黑人xx丫x性爽| 午夜福利18| 欧美xxxx性猛交bbbb| 国产欧美日韩一区二区三| 欧美乱色亚洲激情| 真人一进一出gif抽搐免费| 美女免费视频网站| 国产欧美日韩精品一区二区| 97热精品久久久久久| 国产午夜精品久久久久久一区二区三区 | 色哟哟·www| 国语自产精品视频在线第100页| 国产日本99.免费观看| 深夜精品福利| 国产白丝娇喘喷水9色精品| 欧美色欧美亚洲另类二区| 亚洲熟妇中文字幕五十中出| 99久久九九国产精品国产免费| 中文字幕精品亚洲无线码一区| 精品乱码久久久久久99久播| 黄色配什么色好看| 国内精品久久久久久久电影| 五月玫瑰六月丁香| 亚洲av二区三区四区| 欧美黄色淫秽网站| 91麻豆精品激情在线观看国产| 日本熟妇午夜| 在线免费观看不下载黄p国产 | 欧美色视频一区免费| 国产伦精品一区二区三区视频9| 国产成+人综合+亚洲专区| 成人av在线播放网站| 精品一区二区免费观看|