999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

纖維體積分數對三維編織復合材料T型梁模態性能的影響

2015-03-16 09:19:28李嘉祿
材料工程 2015年9期
關鍵詞:編織模態復合材料

王 歡,李嘉祿,樊 威

(天津工業大學 復合材料研究所 天津市和教育部共建先進紡織復合材料重點實驗室,天津 300387)

?

纖維體積分數對三維編織復合材料T型梁模態性能的影響

王 歡,李嘉祿,樊 威

(天津工業大學 復合材料研究所 天津市和教育部共建先進紡織復合材料重點實驗室,天津 300387)

采用懸臂梁自由衰減振動的實驗方法研究了碳纖維三維四向編織復合材料T型梁的模態性能,分析了纖維體積分數對其固有頻率和阻尼的影響。將實驗得到的固有頻率結果與有限元得到的固有頻率結果作了對比。結果表明:三維四向編織復合材料T型梁的模態性能與激振點位置無關。且其固有頻率隨纖維體積分數的增大而提高,阻尼性能隨纖維體積分數的增大而降低。通過陣型可以準確判斷T型梁固有頻率的階數。實驗結果與計算結果相一致。

三維四向編織;復合材料T型梁;模態分析

三維編織復合材料是由具有多軸纖維取向的高度整體化的連續纖維集合體增強的復合材料,由于增強體的整體性,顯著提高了復合材料整體的剛度和強度, 并且使復合材料具有減振性好、抗疲勞性佳、材料性能可設計等特點。因此可以在許多部件上替代金屬材料以滿足航空航天領域減重和承載性能的要求,也為其應用于主承力結構件提供了深遠而廣闊的應用前景[1,2]。隨著它在航空、航天等領域的廣泛應用,研究三維編織復合材料的模態性能,提高其在振動和噪聲環境下的可靠性,從而避免在使用中因共振造成不必要的損失是十分必要的。

目前,關于這類材料的模態性能已有科研人員進行研究。蔡敢為等[3]研究了三維編織纖維增強復合材料連桿機構的模態阻尼,導出了此類復合材料的細觀結構參數、組分材料性能參數、纖維體積比等與機構的固有頻率、振型和模態阻尼等參數之間的關系。李典森等[4]對不同編織角、不同體積分數的三維四向和五向玻璃纖維編織復合材料進行了振動阻尼實驗,分析了編織角、纖維體積分數和編織結構對三維編織復合材料阻尼特性的影響。孫雨果等[5]在建立有限元模型基礎上,對三維碳/環氧編織復合材料等網格承力筒結構的靜力學性能和模態進行了分析,得到了該結構的固有頻率和模態振型。裴曉園等[6,7]研究了編織角及纖維體積分數對碳纖維三維四向編織復合材料的模態性能影響,實驗表明,隨纖維體積分數的增大,碳纖維三維四向編織復合材料的固有頻率增大,阻尼性能下降。Gao等[8]采用實驗方法研究了編織角對三維五向編織復合材料模態性能的影響。Ooijevaar等[9]將振動實驗與模態應變能損傷算法相結合,對鋪層T型梁的損傷位置做出了準確的定位。研究得出:彎曲模態固有頻率可表明分層的存在,部分陣型可表明分層的存在及其位置。Herman等[10]利用振動模態分析方法檢測了T型加筋板的損傷。Lee等[11]基于經典的層合理論研究了工字形截面復合材料薄壁梁的振動。Thinh等[12]用實驗方法和有限元方法研究了帶有不同形狀加強筋的玻璃纖維增強聚酯基復合材料的自由振動。Alnefaie[13]利用三維有限元模型計算了分層纖維增強復合材料板的固有頻率和模態位移。本工作對三維編織復合材料T型梁制件的模態性能進行了研究。采用單輸入/單輸出和單輸入/多輸出的頻域模態參數識別方法,對三維四向碳纖維編織復合材料T型梁的模態性能進行了研究,確定了纖維體積分數對復合材料T型梁模態參數( 阻尼比、固有頻率)的影響,并將實驗獲得的固有頻率結果和有限元計算結果作了對比。

1 實驗描述

1.1 測試試件

三維四向編織復合材料T型梁試件采用四步法編織工藝,再經RTM工藝復合固化制得。編織紗為 T700-12K 碳纖維,密度1.76g/cm3,線密度0.8g/m,基體材料為TDE-86 型環氧樹脂。圖1為試件的實物圖,實驗中T型梁試件的長度為280mm,其中實驗的有效尺寸為255mm。試件的詳細尺寸如圖2所示,其具體參數見表1。

圖1 三維四向編織復合材料T型梁試件實物圖Fig.1 The picture of specimen of three-dimension and four-direction braided composites T-beam

圖2 T型梁試件的尺寸Fig.2 The size of T-beam specimen

NumberofsampleSurfacebraidingangle/(°)Volumefractionoffiber/%QuantityV129383V229513V329613

1.2 測試方法

本實驗采用一端固支懸臂梁的自由振動衰減法測量T型梁模態性能,振動模式為y向彎曲振動,實驗環境溫度為室溫。實驗裝置示意圖如圖3所示。用力錘在激振點沿y方向激振,力錘上的力傳感器感受的力信號經電荷放大器放大后輸入到數據采集儀的1 號通道,同時加速度傳感器測量的加速度響應信號經電荷放大器放大后輸入到數據采集儀的 2 號通道,再由Vib′SYS模態分析軟件進行分析計算,得到頻率響應函數。本次實驗還選取了不同的激振點進行測試,A點(圖2)和B點為所選取的不同的激振點,其中A點位于遠離腹板一側的肋板平面上,B點位于腹板外側沿Z軸的中心線上,且它們都距T型梁的頂端245mm,如圖3所示。在實驗中分別采用單點激勵單點響應和單點激勵多點響應的方法。

2 實驗結果與討論

以B點為激振點,分別對不同纖維體積分數的三維四向編織碳纖維/環氧樹脂復合材料T型梁進行懸臂梁振動實驗模態分析(每一種纖維體積分數,有3個試件進行實驗),獲得了各個模態參數的平均值( 阻尼比、固有頻率),見表 2。

圖3 模態實驗測試裝置圖Fig.3 The device for modal experimental measurement

NumberofsampleNaturalfrequency/HzDampingratio/%FirstorderSecondorderFirstorderSecondorderV1213.871408.281.300.80V2261.281616.061.100.53V3290.151773.000.970.37

2.1 不同激振點的測試

圖4 T型梁在A,B點激勵后頻率響應函數圖的對比Fig.4 The comparison of frequency responses of T-beam excited at points A and B

本實驗選擇A,B兩個不同的激振點對三維四向編織復合材料T型梁的模態性能進行測試。圖4是分別對纖維體積分數為61%的T型梁的A,B兩點進行激振,得到的頻率響應函數圖,這兩條曲線的吻合性較好。獲得的峰值中最大的相對標準偏差為0.12%。由于所有的相對標準偏差都不超過1%,所以實驗中得到的固有頻率值具有很好的重復性,這也與Ooijevaar等[9]得到的結果相一致。因此在之后的實驗中就選取B點為激振點。

2.2 纖維體積分數對三維四向編織復合材料T型梁固有頻率的影響

圖5是相同表面編織角、不同纖維體積分數的三維四向編織復合材料T型梁試件的前兩階固有頻率的對比,可以看出前兩階固有頻率都隨著纖維體積分數的增加而增大。根據梁的振動理論,懸臂梁的一階振動頻率服從以下關系[14]:

(1)

圖5 不同纖維體積分數T型梁固有頻率的對比Fig.5 The comparison of natural frequencies of T-beam with different fiber volume fractions

式中:k1為對應一階頻率的系數(取1.875);l為梁的長度, mm;E為梁的彈性模量,GPa;h為平行于梁振動方向的高度,mm;ρ為梁的質量密度,g/mm3。由公式(1)可見,當梁的外部尺寸相同時,梁的一階振動固有頻率正比于梁的比模量的平方根。由于三維四向編織復合材料T型梁的纖維剛度遠遠大于基體的剛度,所以纖維體積分數的增大就意味著整體剛度的增大,即梁的彈性模量增大。因此三維四向編織復合材料T型梁的固有頻率隨纖維體積分數的增大而升高,故實驗結果是合理的。

圖6是不同纖維體積分數的三維四向編織復合材料T型梁在激振力作用下的傳遞函數幅值圖。可以看出隨纖維體積分數的增大,各階峰值所對應的固有頻率增大。固有頻率是由系統本身的質量、形狀和剛度所決定[4],相同尺寸的碳纖維三維四向編織復合材料T型梁,其材料的彈性模量越大,則剛度越高,因此梁的固有頻率也越大。固有頻率越高表明其在脈沖載荷作用下的抵抗變形能力越大,即其抗彎剛度大。碳纖維三維四向編織復合材料T型梁的固有頻率越高,在實際使用中就越不容易達到其共振頻率,即避免發生共振的可能性越大。所以為提高三維四向編織復合材料T型梁的固有頻率,在實際應用中應適當地增加其纖維體積分數。

圖6 不同纖維體積分數T型梁頻率響應函數圖Fig.6 The frequency response plot of T-beam with different fiber volume fractions

2.3 纖維體積分數對三維四向編織復合材料T型梁阻尼比的影響

阻尼通常是指材料內部在經受振動變形過程中,將機械振動能轉變為熱能耗散掉的能力。它是試件的固有特性,在振動過程中耗散能量,對降低共振區峰值十分重要[15],可以避免結構因振動應力達到極限而造成的破壞。阻尼性能越好,將振動能轉化成為熱能耗散掉的能力越強,這是阻尼減振性能的本質所在。因此,三維四向編織復合材料T型梁的結構設計中同樣要考慮阻尼因素。復合材料制件的剛度由纖維增強體確保,在保證復合材料制件剛度的同時,也需要適當增大其阻尼性能,以獲得良好的減振效果。

阻尼比可以通過半功率帶寬的方法計算得到,如圖7[16]所示,同時根據以下公式可以得到其阻尼比:

(2)

式中:ωn是由極值點對應的固有頻率;Δω是ω2與ω1兩個點的差值。故結合公式(2)可以從圖6中看出試件的固有頻率越小,對應的峰值越低、越寬,所得的阻尼比越大;試件的固有頻率越大,對應的共振峰值越高、越尖銳,所得的阻尼比越小。

圖7 根據半功率帶寬的方法定義ω1,ω2與ωn[15]Fig.7 Definitions of ω1,ω2 and ωn according to the half-band width method[15]

圖8 不同纖維體積分數T型梁阻尼比的對比Fig.8 The comparison of damping ratios of T-beam with different fiber volume fractions

圖8是相同表面編織角、不同纖維體積分數的三維四向編織復合材料T型梁試件的前兩階阻尼比的對比,可以看出前兩階阻尼比都隨著纖維體積分數的增加而降低。這是由于基體良好的黏彈性質,使得樹脂基體具有較高的振動衰減損耗因子,它在復合材料中起到很好的阻尼作用。所以纖維體積分數越低,制件的黏彈性效應也越明顯。反之隨著纖維體積分數的增加,即基體含量的減少,使得基體對復合材料制件的衰減損失因子貢獻減小。因此,纖維體積分數低的試件的振動衰減性能比纖維體積分數高的試件優異,即前者的阻尼比要比后者的大。這也說明了在相同的共振情況下,纖維體積分數低的T型梁比纖維體積分數高的T型梁回到平衡位置的時間要短。從圖8中還可以看出相同纖維體積分數的第二階阻尼比都比其第一階阻尼比低,這可能是由于梁在二階振動時,接近梁的中心存在很大的振動幅值[13]。

2.4 三維四向編織復合材料T型梁的陣型

振型是結構所固有的屬性,它不取決于作用在結構上的力或者載荷。如果結構的材料屬性(質量、剛度、阻尼)或者邊界條件發生變化,陣型就會改變。陣型是指在共振時一個點相對于另一個點的運動情況,它是獨一無二的。對應于頻率而言,一個固有頻率對應于一個振型。按照頻率從低到高排列為一階振型,二階振型等。實際結構的振動形態并不是一個規則的形狀,而是各階振型相疊加的結果,然而在共振頻率或者共振頻率附近的振動形態是以一個單一的陣型為主。在實際應用中,可以用陣型來判定損傷的位置。本實驗采用單點激勵多點響應的方法研究三維四向編織復合材料T型梁的陣型,激振點有1個(B點),響應點有5個,它們將T型梁的有效長度沿著z軸的方向平均分成5份,如圖9所示。圖10(a),(b)分別是纖維體積分數為61%的T型梁的一、二階彎曲陣型圖。通過測量三維四向編織復合材料T型梁的陣型,可以看出它的振動形態,也可以更加準確地斷定其固有頻率的階數。

圖9 T型梁上的激振點和5個測量點Fig.9 Excitation point and 5 measurement points at T-beam

圖10 實驗測量的彎曲陣型 (a)一階彎曲陣型;(b)二階彎曲陣型Fig.10 Experimental measured bending mode shape (a)first bending mode shape;(b)second bending mode shape

3 有限元計算與模擬

本工作運用ANSYS軟件建立了T型梁的實體有限元模型,計算與模擬了T型梁的固有頻率和陣型。按照復合材料的理論,三維四向編織復合材料是橫觀各向同性的,所以在模擬中選取的材料屬性為orthotropic。采用六面體劃分網格,共劃分了52480個單元,129283個節點。對T型梁的端部固定約束,然后計算。計算中采用的彈性性能數據來自于李典森等[17]的有限元計算,三維四向編織復合材料的彈性性能見表3。T型梁的前兩階彎曲陣型位移云圖如圖11所示,計算的頻率結果與實驗的頻率結果對比如表4所示。

表3 三維四向編織復合材料的彈性性能[17]

圖11 陣型位移云圖 (a)一階彎曲陣型;(b)二階彎曲陣型Fig.11 Cloud pictures of mode shape displacements (a)first bending mode shape;(b)second bending mode shape

No.Naturalfrequency/HzFirstorderSecondorderV3-278.431646.00V3290.151773.00V3+301.331781.20

在有限元計算時不僅可以得到彎曲陣型,還可以得到扭轉陣型。但因本工作只研究彎曲振動,所以只選取了有限元計算中彎曲振動的結果來研究。由于理論計算中采用的T型梁細觀結構參數與實驗中的并不相同,不過表面編織角大體一致,所以在計算中選取了兩種纖維體積分數,以更準確地估計實際中T型梁的固有頻率。圖11為有限元模態計算得到的T型梁的第一階和第二階彎曲陣型。與實驗模態分析得到的第一階和第二階陣型相比,有限元模擬得到的彎曲陣型更加生動形象、直觀,并可以觀察到T型梁上任意點的位移變化。但實驗得到的彎曲陣型和有限元計算得到的彎曲陣型的振動方式是一致的。此外,與陣型對應的實驗測量和有限元計算的固有頻率也具有較好的一致性。

表4中V3是實驗測量的固有頻率,V3-,V3+是采用文獻[17]的數據計算的結果,其纖維體積分數的關系是V3-< V3< V3+。通過對比,可以看出任意階實驗得到的固有頻率都介于計算得到的兩個固有頻率之間,且與它們之間的差值百分比最大不超過8%,所以若理論計算中采用的T型梁細觀結構參數與實驗中的相同,其差值百分比也必然不超過8%,故實驗得到的固有頻率結果是合理的。此外,從實驗與計算結果的一致性還可以得出對于三維四向編織的復雜結構,可以通過有限元計算來大致地預測其振動的固有頻率。

4 結論

(1)表面編織角相同時,三維四向編織結構的碳纖維/環氧樹脂復合材料T型梁的固有頻率隨纖維體積分數的增加而增大,阻尼性能隨纖維體積分數的增加而降低。

(2)由于三維四向編織結構的碳纖維/環氧樹脂復合材料T型梁的抗激振性能隨著纖維體積分數的增加而增大,減振性能隨著纖維體積分數的增加而減小,所以在實際應用中應綜合考慮,合理地選擇其纖維體積分數。

(3)通過觀察三維四向編織結構的碳纖維/環氧樹脂復合材料T型梁的陣型可以準確地判斷其固有頻率的階數。

(4)可以利用有限元計算來大致地預測三維四向編織結構的碳纖維/環氧樹脂復合材料T型梁或者其他復雜結構的振動固有頻率。

[1] 陳光偉,陳利,李嘉祿,等. 三維多向編織復合材料T型梁抗彎應力分析[J]. 紡織學報,2009,30(8):54-58.

CHEN Guang-wei, CHEN Li, LI Jia-lu, et al. Study on flexural stress of T beam of 3-D multidirectional braided composites[J].Journal of Textile Research,2009,30(8):54-58. [2] 李典森,盧子興,李嘉祿,等. 三維編織T300/環氧復合材料的彎曲性能及破壞機理[J]. 航空材料學報,2009,29(5):82-87.

LI Dian-sen, LU Zi-xing, LI Jia-lu,et al. Bending properties and failure mechanism of three dimensional T300/epoxy braided composites[J].Journal of Aeronautical Materials, 2009,29(5):82-87.

[3] 蔡敢為,鐘掘,廖道訓. 三維編織復合材料構件的機構模態阻尼[J]. 中國機械工程,2000,11(5):481-484.

CAI Gan-wei, ZHONG Jue, LIAO Dao-xun. Modal damping of mechanisms fabricated from three-dimensional braided composite[J].China Mechanical Engineering,2000,11(5): 481-484.

[4] 李典森,盧子興,李嘉祿,等. 三維編織復合材料振動阻尼特性的實驗研究[J]. 機械強度,2009,31(2):211-214.

LI Dian-sen, LU Zi-xing, LI Jia-lu, et al. Experimental research on the vibration damping properties of 3D braided composites[J].Journal of Mechanical Strength,2009,31(2):211-214.

[5] 孫雨果,嚴實,梁紅. 新型三維編織復合材料等網格結構力學性能分析[J]. 宇航學報,2007,28(4): 827-830.

SUN Yu-guo, YAN Shi, LIANG Hong. The mechanical behavior study of a new type of 3D braided composites iso-lattice truss[J]. Journal of Astronautics, 2007,28(4): 827-830.

[6] 裴曉園,李嘉祿,何玉強. 纖維體積分數對炭纖維三維四向編織復合材料實驗模態性能的影響[J]. 固體火箭技術,2012,35(2):267-270.

PEI Xiao-yuan, LI Jia-lu, HE Yu-qiang. Effects of fiber volume fraction on experiment modal behavior of three-dimension and four-direction braiding composites[J].Journal of Solid Rocket Technology,2012,35(2):267-270.

[7] 裴曉園,李嘉祿,何玉強. 纖維取向對三維四向編織復合材料與層合復合材料振動性能的影響[J]. 材料工程,2013,(7):16-23.

PEI Xiao-yuan, LI Jia-lu, HE Yu-qiang. Effects of fiber orientation on modal experimental analysis of three-dimension and four-direction braided composites and laminated composites[J]. Journal of Materials Engineering,2013,(7):16-23.

[8] GAO Y, LI J L. Effects of braiding angle on modal experimental analysis of three-dimensional and five-directional braided composites[J]. Composites: Part B,2012,43(5):2423-2428.

[9] OOIJEVAAR T H, LOENDERSLOOT R, WARNET L L,et al. Vibration based structural health monitoring of a composite T-beam[J].Composite Structures,2010,92(9):2007-2015.

[10] HERMAN A P, ORIFICI A C, MOURITZ A P. Vibration modal analysis of defects in composite T-stiffened panels[J].Composite Structures,2013,104(10):34-42.

[11] LEE J, KIM S E. Free vibration of thin-walled composite beams with I-shaped cross-sections[J].Composite Structures,2002,55(2):205-215.

[12] THINH T I, QUOC T H. Finite element modeling and experimental study on bending and vibration of laminated stiffened glass fiber/polyester composite plates[J].Computational Materials Science, 2010,49(4):S383-S389.

[13] ALNEFAIE K. Finite element modeling of composite plates with internal delamination[J].Composite Structures,2009,90(1):21-27.

[14] 鄭麗娟,趙玉濤,金明江,等. 纖維/樹脂/鋁合金疊層復合材料的顯微組織與阻尼性能[J]. 中國有色金屬學報,2003,13(3):667-670.

ZHENG Li-juan, ZHAO Yu-tao, JIN Ming-jiang, et al. Microstructure and damping property of fiber/resin/aluminum alloy laminated composite[J].The Chinese Journal of Nonferrous Metals,2003,13(3):667-670.

[15] 余建新,陳立平. 復合材料結構阻尼測試分析[J]. 哈爾濱商業大學學報:自然科學版,2011,27(5): 729-732.

YU Jian-xin, CHEN Li-ping. Analysis on composite structural damping measurement[J].Journal of Harbin University of Commerce: Natural Sciences Edition,2011,27(5): 729-732.

[16] KHAN S U, LI C Y, SIDDIQUI N A. Vibration damping characteristics of carbon fiber-reinforced composites containing multi-walled carbon nanotubes[J].Composites Science and Technology,2011,71(12):1486-1494.

[17] 李典森,盧子興,藺曉明,等. 三維四向編織復合材料彈性性能的有限元預報[J]. 北京航空航天大學學報,2006,32(7):828-832.

LI Dian-sen, LU Zi-xing, LIN Xiao-ming, et al. Finite element prediction of elastic properties of three dimensional and four directional braided composites[J].Journal of Beijing University of Aeronautics and Astronautics,2006,32(7):828-832.

Effect of Fiber Volume Fraction on Modal Properties of Three-dimension Braided Composite T-beams

WANG Huan,LI Jia-lu,FAN Wei

(Tianjin and Ministry of Education Key Laboratory of Advanced Textile Composite Materials,Institute of Composite Materials, Tianjin Polytechnic University,Tianjin 300387,China)

The modal properties of carbon fiber three-dimension and four-direction braided composite T-beams with different fiber volume fractions were studied by the free vibration method of cantilever beams. The effect of fiber volume fraction on natural frequency and damping was analyzed. The measured and finite element calculated results of natural frequency were compared. The results show that the modal properties of T-beams have nothing to do with the position of excitation point. The natural frequency of the composite T-beams increases and damping properties of the composite T-beams decrease with the increasing of fiber volume fraction. The order of the natural frequency of T-beam could be estimated accurately by the mode shape. The measured results are in agreement with the calculated results.

three-dimension and four-direction braiding;composite T-beam;modal analysis

10.11868/j.issn.1001-4381.2015.09.013

TB332

A

1001-4381(2015)09-0080-07

天津市科技計劃項目(11ZCKFSF00500,10SYSYJC27800)

2014-05-15;

2015-03-21

李嘉祿(1947-),男,教授,主要從事紡織復合材料結構及性能的研究,聯系地址:天津市西青區賓水西道399號天津工業大學復合材料研究所(300387),E-mail:lijialu@tjpu.edu.cn

猜你喜歡
編織模態復合材料
編織一座橋
竹自清高,編織美好
民機復合材料的適航鑒定
復合材料無損檢測探討
電子測試(2017年11期)2017-12-15 08:57:13
度假吧!帶上你的編織鞋包
Coco薇(2017年7期)2017-07-21 07:16:04
國內多模態教學研究回顧與展望
基于HHT和Prony算法的電力系統低頻振蕩模態識別
TiO2/ACF復合材料的制備及表征
應用化工(2014年10期)2014-08-16 13:11:29
由單個模態構造對稱簡支梁的抗彎剛度
計算物理(2014年2期)2014-03-11 17:01:39
RGO/C3N4復合材料的制備及可見光催化性能
主站蜘蛛池模板: 国产午夜小视频| 婷婷综合在线观看丁香| 国产另类视频| 午夜国产大片免费观看| 欧美成人精品在线| 亚洲久悠悠色悠在线播放| 一本色道久久88综合日韩精品| 一个色综合久久| 在线毛片免费| 九九久久99精品| 国产精品妖精视频| 男女男精品视频| 亚洲成人在线网| 国产欧美日韩va| 亚洲第一综合天堂另类专| 性欧美在线| 成人综合久久综合| 亚洲精品天堂自在久久77| 中文字幕啪啪| 亚洲欧美日韩中文字幕一区二区三区| 久久semm亚洲国产| 久久久久国产精品免费免费不卡| 尤物精品视频一区二区三区| 色综合久久综合网| 久久精品最新免费国产成人| 国产精品无码翘臀在线看纯欲| 国产成人免费视频精品一区二区| 精品国产自| 91精品亚洲| 在线视频亚洲欧美| 亚洲AV成人一区国产精品| 在线观看亚洲天堂| 欧美第二区| 亚洲第一在线播放| 欧美激情视频二区| 激情爆乳一区二区| 国产一级毛片yw| 国产sm重味一区二区三区| 无码精油按摩潮喷在线播放 | 亚洲免费三区| 一级全免费视频播放| 欧美成人第一页| 狠狠综合久久| 国产在线八区| 黄色免费在线网址| 波多野结衣一区二区三区88| 亚洲三级电影在线播放| 丝袜久久剧情精品国产| 88国产经典欧美一区二区三区| 亚洲人成亚洲精品| 中文字幕1区2区| 一级全黄毛片| 亚洲国产成人精品青青草原| 1024你懂的国产精品| 秘书高跟黑色丝袜国产91在线| 欧美高清日韩| 亚洲无码视频图片| 亚洲美女久久| 亚洲成人黄色网址| 真实国产精品vr专区| 国产原创演绎剧情有字幕的| 国产成人乱无码视频| 亚洲天堂免费| 欧美激情成人网| 19国产精品麻豆免费观看| 国产91小视频| 亚洲丝袜第一页| 国产精品一区在线观看你懂的| 日本黄色a视频| 亚洲天堂在线视频| 国产精品乱偷免费视频| 亚洲最大福利视频网| 国产视频自拍一区| 国产超薄肉色丝袜网站| 国产日韩久久久久无码精品| 国产18在线播放| 欧美成人二区| 欧美国产日韩另类| 亚洲国产日韩欧美在线| 亚洲乱码在线视频| 麻豆国产在线观看一区二区| 日韩高清中文字幕|