亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        纖維取向?qū)C4/PEEK/Cf層板抗高速?zèng)_擊性能的影響

        2023-02-27 13:35:40高立新袁瀟灑王夢麟薛鵬博吳焱兵
        振動(dòng)與沖擊 2023年4期
        關(guān)鍵詞:層板鋪層交叉

        滕 威,高立新,袁瀟灑,王夢麟,薛鵬博,吳焱兵,潘 蕾

        (1.南京航空航天大學(xué) 材料科學(xué)與技術(shù)學(xué)院,南京 210016; 2.南京航空航天大學(xué) 馬克思主義學(xué)院,南京 210016;3.中國航發(fā)商用航空發(fā)動(dòng)機(jī)有限責(zé)任公司,上海 200241; 4.中航西安飛機(jī)工業(yè)集團(tuán)股份有限公司,西安 710089)

        纖維金屬層板(fiber metal laminates,F(xiàn)MLs)是由樹脂基復(fù)合材料和金屬薄板交替疊層鋪放成預(yù)設(shè)樣式后,在一定工藝下固化成型的一種層間混雜結(jié)構(gòu)[1-2]。TC4/PEEK/Cf層板不僅將Cf/PEEK復(fù)合材料優(yōu)異的抗疲勞性能、高比強(qiáng)度等與TC4金屬薄板的耐蝕性、耐高溫性能等結(jié)合起來,而且其抗沖擊破壞性能、承載能力也較好[3-5]。當(dāng)其應(yīng)用于飛機(jī)的結(jié)構(gòu)件例如高溫壓氣機(jī)葉片、高超音速飛機(jī)殼體、機(jī)翼等難免會出現(xiàn)碰撞沖擊問題,而高速?zèng)_擊的危害比低速?zèng)_擊更加嚴(yán)重。所以研究TC4/PEEK/Cf層板在高速?zèng)_擊下的響應(yīng)特征和損傷破壞機(jī)理等是保障航空安全的重要問題。

        目前,國內(nèi)外對于樹脂基復(fù)合材料和纖維金屬層板的沖擊問題進(jìn)行了大量的研究。樹脂基復(fù)合材料在航空航天領(lǐng)域已被廣泛應(yīng)用,對其抗沖擊性能的研究也較普遍。Higuchi等[6]對碳纖維增強(qiáng)樹脂基復(fù)合材料層板(carbon fiber reinforced polymer,CFRP)進(jìn)行空氣炮高速?zèng)_擊試驗(yàn),研究了其穿透和損傷破壞行為;Katunin等[7-8]對比分析了幾種復(fù)合材料層板低速反復(fù)沖擊行為;汪洋等[9]模擬冰雹高速撞擊復(fù)合材料層板,研究不同參數(shù)對層合板失效模式和損傷情況的影響。

        基于復(fù)合材料層板沖擊問題的研究方法和成果,眾多學(xué)者在纖維金屬層板抗沖擊問題上進(jìn)行了探索。Sangsefidi等[10]研究了剛性、半剛性和可變形彈丸對Glare層板高速?zèng)_擊能量吸收機(jī)制和損傷破壞形貌的影響;Sharma等[11]研究了Glare層板金屬層厚度和鋪放結(jié)構(gòu)對其抗高速?zèng)_擊性能的影響;Zhu等[12]針對FMLs高速?zèng)_擊問題建立了精確描述層板復(fù)雜損傷模式和破壞失效過程的材料損傷本構(gòu)關(guān)系。

        以往在樹脂基復(fù)合材料及FMLs抗沖擊性能方面的研究主要集中在熱固性樹脂基層板,而TC4/PEEK/Cf層板是一種新型熱塑性樹脂基層板,對該層板在高速?zèng)_擊問題上的研究相對較少。本文通過高速?zèng)_擊試驗(yàn),探究纖維取向?qū)C4/PEEK/Cf抗高速?zèng)_擊性能的影響,并對有限元模型進(jìn)行驗(yàn)證,分析層板高速?zèng)_擊下?lián)p傷破壞機(jī)制和沖擊響應(yīng)特征。

        1 試 驗(yàn)

        1.1 試驗(yàn)材料及試樣制備

        制備TC4/PEEK/Cf層板所用材料是TC4金屬薄板、Cf/PEEK預(yù)浸料和PEEK膠膜。本文所用TC4金屬薄板是由江蘇雙威金屬材料有限公司生產(chǎn),單層厚度為0.3 mm;所用單向Cf/PEEK預(yù)浸料是由英國ICI公司生產(chǎn),單層厚度為0.125 mm;所用PEEK膠膜是由英國VICTREX公司生產(chǎn),單層厚度為0.1 mm。3種材料的基本性能參數(shù)分別如表1~表3所示。

        表1 TC4金屬薄板基本性能參數(shù)

        表2 Cf/PEEK預(yù)浸料基本性能參數(shù)

        表3 PEEK膠膜基本性能參數(shù)

        為了提高TC4/PEEK/Cf層板中TC4金屬薄板和Cf/PEEK預(yù)浸料界面結(jié)合強(qiáng)度,對TC4金屬薄板進(jìn)行了噴砂+陽極氧化表面處理??紤]到TC4/PEEK/Cf層板的成型溫度高,常規(guī)熱壓罐成型工藝并不太適用,故使用如圖1所示的XLB-D400×400×1的平板硫化機(jī)及工藝曲線進(jìn)行熱壓成型。

        圖1 熱壓成型設(shè)備及工藝路線Fig.1 Hot pressing molding equipment and process route

        1.2 高速?zèng)_擊試驗(yàn)裝備及方案

        高速?zèng)_擊試驗(yàn)采用空氣炮裝置進(jìn)行,裝置示意圖及實(shí)物圖如圖2所示,由氣泵、炮管、彈拖分離器、紅外速度感應(yīng)器、防護(hù)艙和夾具組成。夾具由上下模組成,四周設(shè)置螺栓固定孔,下模設(shè)置有凹槽,便于放置試樣,且可避免在試樣四周鉆孔裝夾對其性能的損害,夾緊后,開口部分尺寸為120 mm×70 mm。

        層板經(jīng)熱壓制備成型后,按夾具尺寸裁剪為150 mm×100 mm×1.8 mm沖擊試樣。由于高速?zèng)_擊速度大,沖擊物尺寸較小,所以擬采用直徑為9 mm質(zhì)量為3 g的球狀沖擊物。分別在不同的速度點(diǎn)下沖擊不同纖維取向的TC4/PEEK/Cf層板,具體方案如表4所示。編號OE1~OE3中,“O”為纖維取向;“E”為試驗(yàn)研究。

        圖2 空氣炮高速?zèng)_擊裝置示意圖和實(shí)物圖Fig.2 Schematic diagram and physical drawing of air gun high speed impact device

        2 數(shù)值模擬方案

        2.1 有限元模型的建立

        本文采用ABAQUS軟件數(shù)值模擬軟件模擬層板受高速?zèng)_擊過程,在模型網(wǎng)格劃分上采取沖擊中心區(qū)域網(wǎng)格細(xì)化、四周網(wǎng)格粗化的方式來均衡求解效率和求解精度。TC4金屬層、Cf/PEEK復(fù)合材料層和沖擊物均采用C3D8R單元,并采用COH3D8單元模擬沖擊過程中金屬/復(fù)合材料界面響應(yīng)。限制XYZ方向上的所有自由度為0,層板各層之間設(shè)置綁定約束方式(Tie),沖擊物與層板間采用非線性通用接觸方式(general contact explicit),采用罰剛度接觸算法(Penalty)計(jì)算沖擊過程中的潛在接觸,并設(shè)置摩擦因數(shù)為0.3。建立的有限元模型如圖3所示。

        表4 不同纖維取向?qū)影鍥_擊試驗(yàn)方案

        圖3 TC4/PEEK/Cf層板高速?zèng)_擊非線性有限元模型Fig.3 High-speed impact nonlinear finite element model of TC4/PEEK/ Cf laminates

        2.2 材料本構(gòu)關(guān)系及失效模型

        使用Johnson Cook模型建立TC4金屬層本構(gòu)關(guān)系及失效參數(shù)[13-14]以反映應(yīng)變硬化、應(yīng)變率硬化和溫度軟化對金屬強(qiáng)度極限影響;Cf/PEEK是典型的橫觀各向異性纖維增強(qiáng)復(fù)合材料,故采用橫觀各向異性本構(gòu)關(guān)系[15]并使用能較全面反映其在沖擊過程中4種損傷類型的Hashin準(zhǔn)則;高速?zèng)_擊過程中,金屬/復(fù)合材料界面在混合模式(模式I、模式II)下常發(fā)生分層失效,為此,引入基于牽引-分離描述的內(nèi)聚力模型模擬界面分層失效[16]。

        3 結(jié)果與討論

        3.1 損傷模式和失效機(jī)理

        圖4是不同纖維取向的層板未擊穿時(shí)正面、背面及截面損傷形貌??梢园l(fā)現(xiàn),OE1,OE2和OE3層板在沖擊后正面均出現(xiàn)沖擊物壓入留下的凹坑和周向裂紋,背面均出現(xiàn)金屬塑性變形后留下的鼓包和金屬開裂。這可能是由于沖擊物接觸層板時(shí),壓縮應(yīng)力波沿層板厚度方向傳播,在不同介質(zhì)交界處,部分壓縮應(yīng)力波反射回原介質(zhì),并轉(zhuǎn)化為沿層板面內(nèi)徑向傳播的拉伸應(yīng)力波。金屬層在高值壓縮應(yīng)力和拉伸應(yīng)力作用下,正面形成金屬塑性變形凹坑和周向開裂,背面形成鼓包和裂紋。

        OE1,OE2和OE3層板在稍高速度點(diǎn)未擊穿時(shí)截面破壞區(qū)域均呈“V”形,主要損傷破壞模式是鈦合金塑性變形、金屬/復(fù)合材料界面分層和沖擊物作用中心區(qū)域復(fù)合材料碎裂折斷,還可發(fā)現(xiàn)較小部分的復(fù)合材料層間分層。原因可能是壓縮應(yīng)力波在不同介質(zhì)間沿層板厚度方向傳播時(shí),由于介質(zhì)差異會發(fā)生部分反射,反射的壓縮應(yīng)力波轉(zhuǎn)化為沿層板面內(nèi)徑向的拉伸應(yīng)力波,不同介質(zhì)中徑向拉伸應(yīng)力不同,從而產(chǎn)生金屬/復(fù)合材料界面的分層;又由于層板上半部分先接觸沖擊物,可充分耗散沖擊能量,造成下半部分未能充分發(fā)揮耗散沖擊能量的潛力,所以破壞區(qū)域呈“V”形。

        圖5是不同纖維取向?qū)影鍝舸顟B(tài)下正面、背面及截面形貌。可觀察到層板正面被沖擊后形成圓環(huán)狀柱塞和沿層板徑向擴(kuò)展的裂紋,沖擊斷口處金屬呈壓延斷裂特征,背面出現(xiàn)“C”形金屬掀開斷裂。從圖5(a)和圖5(b)中可以觀察到層板正面金屬裂紋主要沿著與其縱向呈一定角度擴(kuò)展,在層板橫向擴(kuò)展的金屬裂紋尺寸較小,但在圖5 (c)中金屬裂紋呈十字交叉狀,另外,OE3層板沖擊后復(fù)合材料碎裂劈開更嚴(yán)重,金屬背面開裂方向與上述特征類似。

        圖4 不同纖維取向?qū)影逦磽舸顟B(tài)下正面,背面及截面形貌Fig.4 Morphology of front, back and section of laminates with different fiber orientations without breakdown

        圖5 不同纖維取向?qū)影鍝舸顟B(tài)下正面,背面及截面形貌Fig.5 Morphology of front, back and section of laminates with different fiber orientations under breakdown state

        對比圖5中3種纖維取向?qū)影褰孛嫘蚊策€可發(fā)現(xiàn),OE1和OE2層板在沖擊物作用中心附近出現(xiàn)嚴(yán)重的復(fù)合材料彎曲、擠壓碎裂和較輕的復(fù)合材料內(nèi)部層間分層,而OE3層板在此區(qū)域的損傷形貌特征與OE1和OE2層板相反,反而是復(fù)合材料內(nèi)部層間分層更嚴(yán)重,且該類型層板整體分層區(qū)域和前兩者相比要大,說明OE3層板在高速?zèng)_擊載荷下可通過更大范圍的損傷來耗散沖擊能量,從而展現(xiàn)出良好的抗沖擊性能,而其他纖維取向?qū)影鍥_擊破壞區(qū)域集中在較小范圍,不能更好的發(fā)揮此作用。

        3.2 有限元模型的驗(yàn)證

        為了驗(yàn)證所建非線性有限元模型的準(zhǔn)確性以及后續(xù)的數(shù)值仿真研究,本文采用直徑9 mm質(zhì)量3 g的軸承鋼圓球沖擊物進(jìn)行5種纖維取向TC4/PEEK/Cf層板的沖擊仿真,具體方案如表5所示,編號OS1~OS5中:“O”為纖維取向;“S”為模擬研究。

        本文分別對比了0°/30°,0°/45°和0°/90° 3種纖維取向的層板在不同沖擊速度下的初始速度-剩余速度擬合曲線、彈道極限和截面形貌的試驗(yàn)與模擬結(jié)果,論證了所建有限元模型的有效性。

        圖6是0°/30°,0°/45°和0°/90° 3種纖維取向的層板沖擊后初始速度-剩余速度擬合曲線試驗(yàn)和模擬結(jié)果??梢园l(fā)現(xiàn),模擬和試驗(yàn)結(jié)果曲線走勢一致,隨著沖擊速度的增加,兩者之間的差距有縮小的趨勢。

        表5 不同纖維取向?qū)影鍥_擊仿真方案

        圖6 不同纖維取向?qū)影宓某跏?剩余速度擬合關(guān)系曲線試驗(yàn)和模擬結(jié)果Fig.6 Experimental and simulation results of initial-residual velocity fitting curves for laminates with different fiber orientations

        表6所示為0°/30°,0°/45°和0°/90°三種纖維取向鋪層的層板在球形沖擊物高速作用下的彈道極限的具體統(tǒng)計(jì)數(shù)據(jù)。其試驗(yàn)和模擬結(jié)果誤差率分別為11.8%,12.1%和12.7%,這低于Schwab等[17]采用試驗(yàn)和數(shù)值模擬方法研究纖維鋪層交叉角度對層板抗沖擊性能時(shí)的誤差率。

        圖7是0°/30°和0°/45°兩種纖維取向的層板沖擊后截面形貌試驗(yàn)和模擬結(jié)果對比。可以發(fā)現(xiàn),模擬結(jié)果截面形貌能大致反映實(shí)際試驗(yàn)時(shí)層板沖擊后的損傷破壞特征,上述結(jié)果進(jìn)一步說明了針對不同纖維取向?qū)影逅ㄓ邢拊P偷挠行院瓦m用性。

        表6 不同纖維取向?qū)影鍙椀罉O限試驗(yàn)和模擬結(jié)果

        圖7 不同纖維取向?qū)影鍥_擊后截面形貌試驗(yàn)和模擬結(jié)果Fig.7 Cross section morphology test and simulation results of different fiber orientations after impact

        3.3 數(shù)值仿真研究

        為評估纖維取向?qū)影鍥_擊性能的影響,采用驗(yàn)證后的有限元模型對不同纖維取向的層板分析其初始速度-剩余速度擬合曲線、彈道極限和能量耗散率曲線來總結(jié)其沖擊響應(yīng)特征。

        圖8是不同纖維取向?qū)影逶诓煌瑳_擊速度作用后初始速度-剩余速度擬合關(guān)系曲線、彈道極限對比以及能量耗散率曲線。圖8(a)和圖8(b)可以發(fā)現(xiàn),在相同剩余速度時(shí),OS3層板所需沖擊速度最大,原因可能是纖維交叉角度越大,層板在面臨高能沖擊壓縮應(yīng)力和由此引起的面內(nèi)拉伸應(yīng)力作用時(shí)能更好的將應(yīng)力分散,發(fā)揮整體耗散沖擊能量作用;OS3層板具有最大彈道極限為171.1 m/s,OS4層板具有最小彈道極限為122.1 m/s,隨著纖維交叉角度的減小,層板的彈道極限呈下降趨勢;圖8(c)中,各層板能量耗散率總體變化趨勢是隨著沖擊速度增加,能量耗散率降低,降低的速率變小,OS3層板相應(yīng)曲線在最右上方,OS4層板相應(yīng)曲線在最左下方,即在相同沖擊速度下,OS3層板能量耗散率最高。

        以上現(xiàn)象說明,纖維交叉鋪放相對于單向鋪放層板在高速?zèng)_擊時(shí)能更有效的耗散沖擊能量,其中0°/90°鋪層層板耗散沖擊能量效率最高,抗沖擊性能最好。

        圖8 不同纖維取向?qū)影逶诓煌瑳_擊速度作用Fig.8 Laminates with different fiber orientations are subjected to different impact velocities

        圖9是不同纖維取向?qū)影逶?80 m/s左右擊穿時(shí)速度和加速度隨時(shí)間變化特征。從圖9(a)可以發(fā)現(xiàn),在沖擊物與層板作用過程中,其速度先急劇降低,然后較緩和降低,最后趨于不變,OS4層板與交叉鋪放層板相比,沖擊后剩余速度更大,纖維交叉鋪層角度越大,層板沖擊后剩余速度越小。說明單向鋪層層板抵抗沖擊物穿透能力較弱,纖維鋪層交叉角度越大層板抗擊穿性能越好,纖維鋪層角度為0°/90°的層板抗擊穿性能最好。

        從圖9(b)可以發(fā)現(xiàn),在沖擊物與層板作用過程中,加速度具有先急劇增加,然后逐漸降低的趨勢,在0~20 μs內(nèi),層板對沖擊物減速效果達(dá)到最大,這與圖9(a)中曲線變化規(guī)律相對應(yīng);OS4層板對沖擊物的減速效果弱于交叉鋪層層板,表現(xiàn)為OS4層板加速度曲線基本在交叉鋪層層板之上;隨著纖維交叉角度增加,層板減速效果變大,其中OS3層板具有最好的減速效果,其加速度最大為-3 310 km/s2。原因可能是纖維單向鋪層時(shí),層板容易發(fā)生局部破壞,不能很好的發(fā)揮整體耗散沖擊能量作用,而交叉鋪層能更好的調(diào)動(dòng)層板整體耗散沖擊能量作用。

        圖10、圖11分別選取了OS2和OS3層板在相似速度沖擊下7.5 μs時(shí)刻的瞬時(shí)應(yīng)力分布和損傷情況。可以發(fā)現(xiàn),OS2層板TP1,TP2層應(yīng)力分布在與層板橫向呈一定角度的橢圓中,而OS3層板TP層瞬時(shí)應(yīng)力分布在較規(guī)則的圓形區(qū)域內(nèi),說明纖維取向會影響層板高速?zèng)_擊過程中金屬層應(yīng)力波的擴(kuò)展。

        圖9 不同纖維取向?qū)影逶?80 m/s左右擊穿時(shí)Fig.9 When the laminates with different fiber orientations break down at about 180 m/s

        CP層應(yīng)力分布特征同樣受到纖維取向的影響,其分布規(guī)律與上述TP層類似。CP1分層損傷區(qū)域小于CP2層,可能是由于沖擊物先與層板上部作用,沖擊能量較大,沖擊物可瞬時(shí)穿過層板上部,待與層板下部作用時(shí),沖擊物剩余能量不足以瞬時(shí)擊穿層板,而充分與其作用,因此CP2分層更嚴(yán)重。與OS2層板相比,OS3層板在稍高速度點(diǎn)下,無論正面還是背面CP層的分層區(qū)域均比OS2層板相應(yīng)位置要小,說明增大纖維鋪層交叉角度有利于減小金屬/復(fù)合材料界面分層損傷區(qū)域。

        注:TP.TC4金屬層;CP.內(nèi)聚力層;FP.復(fù)合材料層。圖10 OS2層板在196 m/s速度下7.5 μs時(shí)刻各層應(yīng)力云圖Fig.10 Stress nephogram of OS2 laminates at 7.5 μs time at 196 m/s speed

        注:TP.TC4金屬層;CP.內(nèi)聚力層;FP.復(fù)合材料層。圖11 OS3層板在201 m/s速度下7.5 μs時(shí)刻各層應(yīng)力云圖Fig.11 Stress nephogram of OS3 laminates at 7.5 μs time at 196 m/s speed

        從FP層可以發(fā)現(xiàn),沖擊時(shí)復(fù)合材料層中應(yīng)力波主要沿著纖維方向擴(kuò)展,層板上部FP1-FP2單層復(fù)合材料中應(yīng)力分布范圍小于層板下部FP7-FP8中相應(yīng)范圍,原因可能與上述類似。

        4 結(jié) 論

        (1) 不同纖維取向的TC4/PEEK/Cf層板在高速?zèng)_擊未擊穿狀態(tài)下正面形成凹坑和周向裂紋,背面出現(xiàn)鼓包、金屬開裂和復(fù)合材料折斷,損傷破壞模式主要是復(fù)合材料層間分層和金屬塑性變形;擊穿狀態(tài)下,層板正面損傷區(qū)域較小,出現(xiàn)沿面內(nèi)徑向擴(kuò)展的裂紋、周向擴(kuò)展的金屬壓延斷口,背面損傷區(qū)域較大,出現(xiàn)花瓣?duì)罱饘匍_裂和復(fù)合材料斷裂,損傷破壞模式主要是金屬/復(fù)合材料界面分層、復(fù)合材料碎裂折斷及其層間分層。

        (2) TC4/PEEK/Cf層板在球形沖擊物作用下試驗(yàn)和仿真彈道極限誤差在12%左右,試驗(yàn)和模擬結(jié)果存在一定差異,但誤差率均在現(xiàn)有技術(shù)水平之內(nèi)。且數(shù)值仿真結(jié)果和實(shí)際試驗(yàn)時(shí)層板受高速?zèng)_擊的響應(yīng)特征和損傷破壞形貌較為吻合,驗(yàn)證了所建有限元模型的有效性和適用性。

        (3) 數(shù)值模擬結(jié)果表明,TC4/PEEK/Cf層板受沖擊后的的能量耗散率和彈道極限隨著鋪層交叉角度的增加而增大,剩余速度則隨著鋪層交叉角度的增加而減小。同時(shí)纖維鋪層交叉角度會影響層板應(yīng)力波的擴(kuò)展,增大鋪層交叉角度有利于減小金屬/復(fù)合材料界面分層損傷區(qū)域。以上現(xiàn)象說明增加纖維鋪層交叉角度可提高層板抗沖擊性能,0°/90°鋪層層板抗外物擊穿能力最好。

        猜你喜歡
        層板鋪層交叉
        纖維金屬層板成形極限測試的新試樣設(shè)計(jì)
        “六法”巧解分式方程
        復(fù)合材料層板連接受彎釘載分配方法研究
        衛(wèi)星天線復(fù)合材料框架的鋪層優(yōu)化設(shè)計(jì)
        復(fù)合材料軸結(jié)構(gòu)力學(xué)性能預(yù)測及鋪層方案設(shè)計(jì)
        連一連
        玻璃纖維增強(qiáng)鋁合金層板連接孔擠壓性能實(shí)驗(yàn)研究
        CPD鋪層自動(dòng)創(chuàng)建技術(shù)
        基于Fast-ICA的Wigner-Ville分布交叉項(xiàng)消除方法
        FIBERSIM軟件在工字梁鋪層設(shè)計(jì)上的應(yīng)用
        教練機(jī)(2015年2期)2015-04-03 06:51:48
        亚洲国产香蕉视频欧美| 亚洲av成人噜噜无码网站| 国产熟妇人妻精品一区二区动漫| 五月中文字幕| 日本久久精品国产精品| 亚洲一二三区免费视频| 蜜臀色欲av在线播放国产日韩| 久久久国产精品麻豆| 国产三级在线看完整版| 国产实拍日韩精品av在线| 欧美精品黑人粗大免费| 亚洲产国偷v产偷v自拍色戒| 强d乱码中文字幕熟女1000部| 青青草小视频在线观看| 一本色道久久综合无码人妻| 亚洲人成影院在线无码观看| 日韩丝袜人妻中文字幕| 一区二区三区国产在线视频| 末发育娇小性色xxxx| 日韩AV有码无码一区二区三区| 亚洲一区二区在线视频,| 在线观看午夜视频一区二区| 亚洲av成人一区二区三区| 国产精品美女黄色av| 亚洲精品大全中文字幕| 国产日韩av在线播放| 996久久国产精品线观看| 日本精品极品视频在线| 美女用丝袜脚玩我下面| 色橹橹欧美在线观看视频高清| 98在线视频噜噜噜国产| 亚洲国产天堂av成人在线播放| 中文字幕亚洲无线码在线一区| 国产精品jizz在线观看老狼| 欧美亚洲尤物久久综合精品| 国产av熟女一区二区三区密桃| 亚洲一区自拍高清亚洲精品| 国产精品一区二区在线观看99| 日韩精品一级在线视频| 69国产成人精品午夜福中文| 午夜福利电影|