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

基于灰色關(guān)聯(lián)與最弱鏈理論的鑄鋁氣缸蓋組織與力學(xué)性能評(píng)估

2025-05-23 00:00:00蒲博聞周海濤王根全毛郭靈胡定云何曉東
車用發(fā)動(dòng)機(jī) 2025年2期
關(guān)鍵詞:力學(xué)性能

氣缸蓋是車用發(fā)動(dòng)機(jī)的核心部件之一,在發(fā)動(dòng)機(jī)工作時(shí)所受工況條件嚴(yán)苛,往往承受強(qiáng)烈的熱沖擊作用[1]。伴隨車輛輕量化需求,鋁合金氣缸蓋獲得廣泛應(yīng)用[-3]。氣缸蓋的力學(xué)性能、失效概率很大程度上取決于其微觀組織,準(zhǔn)確表征氣缸蓋材料組織與力學(xué)性能、失效概率之間的相關(guān)性具有重要意義。為了滿足發(fā)動(dòng)機(jī)服役性能要求,鋁合金氣缸蓋鑄件內(nèi)部應(yīng)盡量避免產(chǎn)生縮孔、縮松、裂紋等缺陷,因此,對(duì)鋁合金氣缸蓋鑄造工藝的控制和鑄件質(zhì)量的評(píng)估極其必要[4-6]

對(duì)氣缸蓋本體材料的組織與性能進(jìn)行評(píng)估時(shí),需要考慮的特征參量較多,鑄鋁氣缸蓋本體材料的微觀組織包括硅相、孔隙缺陷等,同時(shí)涉及其各自的分布情況與尺寸大小等因素,而對(duì)材料力學(xué)性能評(píng)價(jià)通常要從抗拉強(qiáng)度、硬度、彈性模量等方面進(jìn)行綜合考核。因此,考慮微觀組織特征對(duì)材料力學(xué)性能進(jìn)行預(yù)測(cè)與評(píng)估,實(shí)際上是多要素分析和多目標(biāo)優(yōu)化問題。基于灰色關(guān)聯(lián)理論的分析方法是一種分析灰色系統(tǒng)中各要素關(guān)聯(lián)程度的量化方法,適用于在數(shù)據(jù)量少、信息貧乏的條件下,有效規(guī)避數(shù)據(jù)欠完整性和主觀因素的影響,厘清多要素多目標(biāo)之間的關(guān)聯(lián)性[7-9]。灰色關(guān)聯(lián)分析方法在碳鋼和低合金鋼大氣腐蝕主要環(huán)境影響因素分析、碳纖維增強(qiáng)樹脂復(fù)合材料控制臂最優(yōu)鋪層方案制定、基于材料組成的半剛性材料疲勞性能模型構(gòu)建等領(lǐng)域,均已展現(xiàn)出較好的應(yīng)用效果[10-12]。

鋁合金氣缸蓋鑄件在不同位置存在各種隨機(jī)分布的缺陷,缺陷位置往往存在應(yīng)力集中和非均勻局部應(yīng)力,是疲勞裂紋萌生和擴(kuò)展的敏感區(qū)域,容易誘發(fā)產(chǎn)生疲勞破壞[13]。缺陷對(duì)于疲勞失效的這種敏感性特征常伴隨著明顯的尺寸效應(yīng),即疲勞強(qiáng)度隨著試樣尺寸的增長而減小[14]。由于氣缸蓋內(nèi)部有較多薄壁結(jié)構(gòu),很多部位無法取得標(biāo)準(zhǔn)試樣進(jìn)行力學(xué)性能測(cè)試,常采用小型試樣進(jìn)行替代。然而,小型試樣不可避免地會(huì)出現(xiàn)尺寸效應(yīng),使材料性能受試樣尺寸及結(jié)構(gòu)影響而表現(xiàn)出較大差異性。最弱鏈理論、臨界缺陷方法等是考慮尺寸效應(yīng)的常用方法,其中最弱鏈理論的應(yīng)用較為廣泛,該理論最初由W.

WEIBULL[15]提出并用于脆性或準(zhǔn)脆性材料的高周疲勞壽命和疲勞極限評(píng)估,具有不存在微觀建模過程、計(jì)算成本低的優(yōu)點(diǎn)。趙承章等[16基于最弱鏈理論分別構(gòu)建了鋁合金氣缸蓋本體材料小型試樣和標(biāo)準(zhǔn)試樣抗拉強(qiáng)度的Weibull分布模型,獲得了小型試樣抗拉強(qiáng)度對(duì)相應(yīng)標(biāo)準(zhǔn)試樣抗拉強(qiáng)度的預(yù)測(cè)關(guān)系。最弱鏈理論在系統(tǒng)故障預(yù)測(cè)分析中也得到了驗(yàn)證,如翁家慶等[1]根據(jù)發(fā)動(dòng)機(jī)大修數(shù)據(jù),結(jié)合回歸分析和優(yōu)化技術(shù),建立了發(fā)動(dòng)機(jī)大修里程的兩參數(shù)Weibull分布模型。

本研究在氣缸蓋頂板、立墻和底板位置分別取得標(biāo)準(zhǔn)試樣和小型試樣,通過試樣力學(xué)性能測(cè)試和微觀組織表征,結(jié)合灰色關(guān)聯(lián)分析方法,發(fā)現(xiàn)了鑄鋁氣缸蓋不同位置材料的抗拉強(qiáng)度受內(nèi)部缺陷影響較大,同時(shí)建立了廣義的雙參數(shù)Weibull分布模型,可以有效降低隨機(jī)缺陷影響強(qiáng)度分布的試樣尺寸效應(yīng),便于對(duì)氣缸蓋材料內(nèi)部缺陷與強(qiáng)度進(jìn)行綜合評(píng)估。

1試驗(yàn)材料及方法

用于研究的氣缸蓋材料為鑄造鋁硅合金,其化學(xué)成分如表1所示。鑄鋁氣缸蓋采用低壓鑄造成型并進(jìn)行T6熱處理,氣缸蓋澆注溫度為 ,T6熱處理包括 三級(jí)固溶處理及 時(shí)效處理。對(duì)同一個(gè)氣缸蓋多個(gè)位置進(jìn)行切割加工取樣,通過拉伸性能測(cè)試,對(duì)氣缸蓋整體結(jié)構(gòu)的力學(xué)性能進(jìn)行分析。氣缸蓋單缸結(jié)構(gòu)及取樣位置如圖1所示。分別從頂板、立墻和底板處取樣加工得到相應(yīng)標(biāo)準(zhǔn)拉伸試樣,為了對(duì)比研究,在臨近標(biāo)準(zhǔn)試樣取樣位置也取樣加工得到相應(yīng)小型拉伸試樣(取樣位置如圖1虛線框所示)。合金拉伸性能根據(jù)GB/T228.1—2010《金屬材料拉伸試驗(yàn)第1部分:室溫試驗(yàn)方法》的規(guī)定執(zhí)行測(cè)試,標(biāo)準(zhǔn)試樣和小型試樣尺寸如圖2所示。切取標(biāo)準(zhǔn)拉伸試樣及小型拉伸試樣斷口附近的橫截面 (?6,?3 尺寸的圓形觀測(cè)平面),進(jìn)行機(jī)械打磨和拋光,并利用金相顯微鏡對(duì)材料金相組織進(jìn)行觀察。

表1合金的化學(xué)成分
圖1氣缸蓋及不同位置的示意
圖2試樣尺寸規(guī)格

2 結(jié)果與分析

2.1氣缸蓋不同位置本體材料的顯微組織

圖3示出氣缸蓋不同位置本體材料的顯微組織。從圖3可見,氣缸蓋不同位置表現(xiàn)出材料顯微組織差異性,這是由不同位置具有不同的凝固冷卻速率導(dǎo)致[18]。如圖3箭頭所示,頂板和底板部位的共晶硅分布均勻,呈細(xì)小纖維狀形態(tài);立墻部位的共晶硅分布不均勻,為狹長的針狀形態(tài)。這種狹長針狀形態(tài)的共晶組織容易引起其周圍的鋁晶粒產(chǎn)生應(yīng)力集中而降低材料的強(qiáng)塑性[19]。此外,從圖中還可以觀察到一些沿枝晶臂邊緣形成的、形狀不規(guī)則的孔隙缺陷(圖3虛線圈所示位置),其形成原因一方面是凝固冷卻時(shí)發(fā)生體積收縮,另一方面是溶解氣體在逸出時(shí)未獲得充分補(bǔ)縮[20]。

a氣缸蓋頂板部位
圖3氣缸蓋不同位置本體材料的顯微組織

2.2氣缸蓋微觀組織特征對(duì)力學(xué)性能影響的灰色關(guān)聯(lián)分析

由于氣缸蓋結(jié)構(gòu)的復(fù)雜性,氣缸蓋的鑄造成型引起不同位置的鋁合金微觀組織表現(xiàn)出分散性。以鋁合金氣缸蓋微觀組織特征中的二維孔隙率、孔隙面積、孔隙最大Feret尺寸、硅相面積、硅相長寬比、硅相最大Feret尺寸和硅相顆粒圓度等7個(gè)特征參數(shù)作為比較序列,利用灰色理論分析了氣缸蓋本體材料微觀組織特征分別與抗拉強(qiáng)度、彈性模量和維氏硬度測(cè)試結(jié)果的關(guān)聯(lián)關(guān)系。

用于建立灰色關(guān)聯(lián)分析模型的樣本量為15組,分別取自氣缸蓋頂板、立墻、底板3個(gè)位置,對(duì)每個(gè)樣本的7個(gè)微觀組織特征和抗拉強(qiáng)度、彈性模量、維氏硬度測(cè)試結(jié)果分別計(jì)算平均值,如表2所示。以氣缸蓋本體材料力學(xué)性能測(cè)試結(jié)果為參考序列 0 ,此處分別將材料的抗拉強(qiáng)度、彈性模量和維氏硬度設(shè)置為參考序列。材料的7個(gè)微觀組織特征設(shè)置為比較序列 。對(duì)

分別求平均值,得到 ,計(jì)算 。對(duì)上述計(jì)算處理之后的比較序列分別求 層次上的最小差和最大差:

再針對(duì)比較序列在 層次上的最小差和最大差分別求 n 層次上的最小差和最大差:

將氣缸蓋本體材料的7個(gè)微觀組織特征與抗拉強(qiáng)度、彈性模量、維氏硬度測(cè)試結(jié)果的關(guān)聯(lián)度進(jìn)行比較,結(jié)果如表3所示。其中,氣缸蓋微觀組織特征對(duì)抗拉強(qiáng)度、彈性模量、維氏硬度影響的關(guān)聯(lián)度強(qiáng)弱順序依次是孔隙最大Feret尺寸、硅相最大Feret尺寸、硅相顆粒圓度、硅相長寬比、硅相面積、孔隙面積、二維孔隙率。分析結(jié)果表明,孔隙最大Feret尺寸對(duì)鑄造鋁合金氣缸蓋力學(xué)性能的影響較大。

表2氣缸蓋本體材料微觀組織特征與力學(xué)性能測(cè)試結(jié)果的平均值
表3氣缸蓋本體材料微觀組織特征與力學(xué)性能測(cè)試結(jié)果的關(guān)聯(lián)度
續(xù)表

2.3氣缸蓋不同位置材料強(qiáng)度的雙參數(shù)最弱鏈統(tǒng)計(jì)模型

氣缸蓋本體材料缺陷最大Feret尺寸與抗拉強(qiáng)度存在較強(qiáng)相關(guān)性,說明氣缸蓋不同位置材料的拉伸性能受缺陷的影響較大,缺陷的隨機(jī)性導(dǎo)致了材料強(qiáng)度的分散性,這一結(jié)果符合Weibull分布的最弱環(huán)節(jié)理論相關(guān)假設(shè)[16,21]。因此,可以用基于最弱鏈統(tǒng)計(jì)模型的雙參數(shù)Weibull分布分析材料強(qiáng)度數(shù)據(jù)與失效概率,形式如下式所示:

式中: P 為試樣在抗拉強(qiáng)度小于或等于 σ 情況下的失效概率; V 為試樣在測(cè)試區(qū)間的體積; 為參考體積,通常按單位體積 取用; 為特征強(qiáng)度,即尺寸參數(shù); m 為Weibull模量,即形狀參數(shù); σ 為進(jìn)行材料失效概率估計(jì)時(shí)的抗拉強(qiáng)度。失效概率 P 通常由Bergman統(tǒng)計(jì)估計(jì)方法計(jì)算[22]:

式中: i 為參與統(tǒng)計(jì)試樣的抗拉強(qiáng)度值按遞增順序進(jìn)行排列的排列數(shù); n 為參與統(tǒng)計(jì)的試樣數(shù)量。為便于對(duì)模型進(jìn)行參數(shù)擬合,式(6)改寫為如下形式:

氣缸蓋不同位置標(biāo)準(zhǔn)試樣基于式(8)計(jì)算的數(shù)據(jù)分布如圖4所示,氣缸蓋不同位置小型試樣基于式(8)計(jì)算的數(shù)據(jù)分布如圖5所示,相關(guān)參數(shù)通過函數(shù)擬合得到,詳見表4。

圖4氣缸蓋不同位置標(biāo)準(zhǔn)試樣抗拉強(qiáng)度的Weibull分布
圖5氣缸蓋不同位置小型試樣抗拉強(qiáng)度的Weibull分布
表4氣缸蓋不同位置本體材料的Weibul參數(shù)

事實(shí)上,孔隙率具有尺寸效應(yīng),其在不同尺寸的試樣中具有明顯差異,由于小尺寸試樣在測(cè)試區(qū)間的體積較小,其孔隙率普遍高于標(biāo)準(zhǔn)尺寸試樣。事實(shí)上,氣缸蓋不同位置的孔隙率也不盡相同。同時(shí),缺陷的空間分布并不是均勻的。研究發(fā)現(xiàn),將描述缺陷分布偏離理想均勻分布程度的 β 系數(shù)[23-24]引人式(6)中,建立廣義的雙參數(shù)Weibull分布統(tǒng)計(jì)模型,一方面能夠削弱孔隙率對(duì)力學(xué)性能影響的尺寸效應(yīng),另一方面能準(zhǔn)確表示不同位置材料缺陷均勻性與材料強(qiáng)度的關(guān)系,該模型表示為

式(9)可以轉(zhuǎn)化為如式(10)所示的形式:

氣缸蓋不同位置不同尺寸試樣基于式(10)計(jì)算的數(shù)據(jù)分布如圖6所示。

圖6氣缸蓋不同位置不同試樣尺寸下廣義的雙參數(shù)Weibull分布

由圖6可知, β 值均小于1,說明氣缸蓋不同位置的孔隙缺陷密度均小于缺陷均勻空間分布時(shí)的密度。其中,底板位置的 β 值為0.9,大于立墻和頂板位置的 β 值 ,說明三個(gè)位置中,底板位置的缺陷偏離理想均勻分布的程度更大,容易在局部位置集中分布,具有較高的鑄造缺陷誘導(dǎo)失效概率,在鑄造工藝制定時(shí)應(yīng)特別關(guān)注。

3結(jié)論

a)鋁合金氣缸蓋鑄件不同位置的微觀組織具有差異化特征,頂板和底板部位的微觀組織相對(duì)較好,共晶硅分布均勻,為細(xì)小纖維狀形態(tài);立墻部位的微觀組織較差,共晶硅分布不均,且表現(xiàn)為狹長的針狀形態(tài);

b)揭示了氣缸蓋本體材料的孔隙最大Feret尺寸、硅相最大Feret尺寸、二維孔隙率等7個(gè)微觀組織特征與本體材料的抗拉強(qiáng)度、彈性模量和維氏硬度測(cè)試結(jié)果的關(guān)聯(lián)程度,其中孔隙最大Feret尺寸對(duì)鑄造鋁合金氣缸蓋力學(xué)性能的影響最大;

c)基于最弱鏈統(tǒng)計(jì)模型,引入 β 系數(shù)以消除試樣尺寸對(duì)力學(xué)性能的影響,建立了氣缸蓋不同位置材料強(qiáng)度數(shù)據(jù)與失效概率的歸一化雙參數(shù)Weibull分布模型,模型表征結(jié)果顯示,底板位置容易出現(xiàn)局部缺陷的集中分布,使該部位具有較高的失效概率,在鑄造工藝制定時(shí)應(yīng)予以著重考慮。

參考文獻(xiàn):

[1]佘祖新,張凱,王長朋,等.汽車鑄造鋁合金缸蓋失效分析[J].裝備環(huán)境工程,2018,15(3):81-85.

[2] 陳大輝,賈祥磊,朱秀榮,等.發(fā)動(dòng)機(jī)缸蓋鑄造鋁合金的研究進(jìn)展[J].鑄造技術(shù),2010,31(7):882-887.

[3] 蔡啟舟,駱海賀,魏伯康,等.鋁合金缸蓋鑄造缺陷分析及對(duì)策[J].特種鑄造及有色合金,2007,27(增刊1):70-72.

[4] 廖治東,彭寶斌,李浩如,等.鋁合金氣缸蓋鑄造氣孔缺陷分析及解決方法J].鑄造,2019,68(3):311-314.

[5] 趙葵,何德道.冷鐵在鋁合金缸蓋鑄造工藝中的應(yīng)用[J].特種鑄造及有色合金,2017,37(3):338-339.

[6] 馬欣宇,范云波,陳敏,等.前缸蓋的間接擠壓鑄造技術(shù)研究[J].兵器材料科學(xué)與工程,2010,33(3):70-74.

[7] 錢毅,江超,李楠,等.基于灰色評(píng)估理論的某旋轉(zhuǎn)機(jī)械的可靠性分配J.機(jī)械強(qiáng)度,2021,43(6):1335-1340.

[8] 胡成亮,劉全坤,王強(qiáng),等.基于灰色關(guān)聯(lián)和模糊邏輯的齒輪鍛模多目標(biāo)優(yōu)化設(shè)計(jì)[J].中國機(jī)械工程,2007(14):1739-1742.

[9] 梁強(qiáng),張賢明,杜彥斌,等.基于灰色關(guān)聯(lián)分析的齒環(huán)熱精鍛成形工藝參數(shù)優(yōu)化J.計(jì)算機(jī)集成制造系統(tǒng),2022,28(4):1020-1029.

[10]王海濤,韓恩厚,柯偉.碳鋼、低合金鋼大氣腐蝕的灰色模型預(yù)測(cè)及灰色關(guān)聯(lián)分析[J].腐蝕科學(xué)與防護(hù)技術(shù),2006(4):278-280.

[11]蔣榮超,慈樹坤,劉大維,等.基于灰色關(guān)聯(lián)分析的碳纖維增強(qiáng)樹脂復(fù)合材料控制臂鋪層優(yōu)化[J].復(fù)合材料學(xué)報(bào),2022,39(1):390-398.

[12]周浩,沙愛民.半剛性材料疲勞性能材料組成影響因素分析[J].武漢理工大學(xué)學(xué)報(bào),2012,34(1):41-45.

[13]LUL,CHEN H,REN M,et al.Study on fatigue lifeof aluminum alloy 6o61-T6 based on random defectcharacteristics[J].Materials,2024,17(5):1133.

[14]AIY,ZHU SP,LIAO D.Probabilistic modelling ofnotch and size effect of components under fatigueloadings[J].Procedia structural integrity,2019,22:70-77.

[15]WEIBULL W.A statistical distribution function ofwide applicability[J].Journal of applied mechanics,1951:25-30.

[16]趙承章,黃渭清,左正興,等.鑄造鋁硅合金缸蓋拉伸性能分析及威布爾建模[J].內(nèi)燃機(jī)學(xué)報(bào),2022,40(5):463-470.

[17]翁家慶,張志剛,沈穎剛,等.中型載貨汽車發(fā)動(dòng)機(jī)大修需求預(yù)測(cè)[J].車用發(fā)動(dòng)機(jī),2008(2):20-22.

[18]蒲博聞,王根全,周海濤,等.凝固速率對(duì)鑄造鋁硅合金組織及力學(xué)性能的影響規(guī)律[J].鑄造,2024,73(1):76-80.

[19]VANDERSLUIS E,LOMBARDI A,RAVINDRANC,et al.Factors influencing thermal conductivity andmechanical properties in 319 Al alloy cylinder heads[J].Materials science and engineering A,2015,648:401-411.

[20]JAVIDANI M,LAROUCHE D.Application of castAl-Si alloys in internal combustion engine components[J].International materials reviews,2ol4,59(3):132-158.

[21]WANG W,ZHONG Y,LU K,et al.Size effects andstrength fluctuation in nanoscale plasticity[J].Actamaterialia,2012,60(8):3302-3309.

[22]BERGMAN B.On the estimation of the Weibull mod-ulus[J].Journal of materials science letters,1984,3:689-692.

[23]LEI W S.A generalized weakest-link model for sizeeffect on strength of quasi-brittlematerials[J].Journalof materials science,2018,53(2):1227-1245.

[24]竺宇洋,翟建廣,高春,等.基于廣義雙參數(shù)最弱鏈統(tǒng)計(jì)模型分析纖維斷裂強(qiáng)度與尺寸的相關(guān)性[J].上海工程技術(shù)大學(xué)學(xué)報(bào),2021,35(1):28-32.

Evaluation of Microstructure and Mechanical Properties of Cast Aluminum Cylinder Head Based on Grey Relation and Weakest Chain Theory

PU Bowen1,ZHOU Haitao1,2,WANG Genquan1,2,MAO Guoling1,HU Dingyun,HE Xiaodong2 (1.National Key Laboratory of Vehicle Power System,China North Engine Research Institute, Tianjin,China;2.China North Engine Research Institute(Tianjin),Tianjin ,China)

Abstract:Anevaluationmethodcollctivelyreflectingtheinfluenceof materialtensile strengthand specimen sizeonthefail ureprobabilitywas introduced toauratelycharacterizethecorelationbetweenmicrostructureandfailureprobabilityofcylinderheadmaterial.Firstly,basedontheresultsofmicrostructuralcharacterizationandmechanicalpropertytestingofdifferent positionsincylinderhead,thegreyrelationalanalysismethodwasemployedtorevealtherelationshipbetweenmicrostructural characteristicsandmechanicalpropertiesbyusingtwo-dimensionalporosity,porearea,pore maximumFeretsize,siliconphase area,silicon phase length-widthratio,siliconphase maximumFeret size,silconphaseparticle roundness asrelevant sequences andusingtensile strength,elastic modulus,Vickers hardness asreference sequences.Basedonthe weakest chain theory,a generaliedtwo-parameter Weibulldistributionmodelwithweakenedspecimensizeefectwasbuilttodescribetherelationship between failureprobabilityandmaterialstrengthatdifferentlocationsofcylinderheadbyconsideringthedegreeofdefectdistributiondeviationfromtheidealuniformdistribution.TheresultsshowedthattheporemaximumFeretsizehadarelatively largeroleininfluencingthemechanicalpropertiesofcastaluminumcylinderheadbodymaterial,thebotomplatepositionwas pronetothecentralized distributionof local defects,andtheprobabilityoffailurecaused bycastingdefects was higher.

Key Words:cylinder head;Al-Si aloy;grey relational analysis;defect distribution;weakest chain theory

猜你喜歡
力學(xué)性能
反擠壓Zn-Mn二元合金的微觀組織與力學(xué)性能
Pr對(duì)20MnSi力學(xué)性能的影響
云南化工(2021年11期)2022-01-12 06:06:14
Mn-Si對(duì)ZG1Cr11Ni2WMoV鋼力學(xué)性能的影響
山東冶金(2019年3期)2019-07-10 00:54:00
采用稀土-B復(fù)合變質(zhì)劑提高ZG30MnSi力學(xué)性能
碳纖維增強(qiáng)PBT/ABS—g—MAH復(fù)合材料的力學(xué)性能和流變行為
中國塑料(2016年6期)2016-06-27 06:34:16
紡織纖維彎曲力學(xué)性能及其應(yīng)用
MG—MUF包覆阻燃EPS泡沫及力學(xué)性能研究
中國塑料(2015年12期)2015-10-16 00:57:14
EHA/PE復(fù)合薄膜的力學(xué)性能和阻透性能
中國塑料(2015年9期)2015-10-14 01:12:26
PA6/GF/SP三元復(fù)合材料的制備及其力學(xué)性能研究
中國塑料(2015年4期)2015-10-14 01:09:18
INCONEL625+X65復(fù)合管的焊接組織與力學(xué)性能
焊接(2015年9期)2015-07-18 11:03:53
主站蜘蛛池模板: 久久香蕉国产线看观看亚洲片| 99激情网| 无码一区中文字幕| AV不卡在线永久免费观看| 亚洲精品黄| 国产视频欧美| 欧美日韩一区二区三区在线视频| 亚洲综合第一页| 狠狠色婷婷丁香综合久久韩国| 亚洲性色永久网址| 欧美激情综合| 国产欧美精品午夜在线播放| 欧美19综合中文字幕| 欲色天天综合网| 香蕉精品在线| 亚洲一区二区精品无码久久久| 亚洲第一香蕉视频| 女人av社区男人的天堂| 欧美日韩中文字幕二区三区| 成人国产免费| 久久福利网| 欧美曰批视频免费播放免费| 亚洲欧美成人综合| 亚洲乱伦视频| 人妻丰满熟妇av五码区| 91麻豆精品国产91久久久久| 国产成人乱无码视频| 久久婷婷六月| 日韩欧美国产三级| 亚洲动漫h| 国产精品视频3p| 夜夜拍夜夜爽| 亚洲男人的天堂在线观看| 99re在线观看视频| 国产精品女熟高潮视频| 欧美人在线一区二区三区| 曰AV在线无码| 毛片手机在线看| 久青草网站| 一本久道热中字伊人| 久久精品无码国产一区二区三区| 找国产毛片看| 中文精品久久久久国产网址| 狠狠躁天天躁夜夜躁婷婷| 自慰高潮喷白浆在线观看| 综合久久五月天| 国产精品漂亮美女在线观看| 99热最新网址| 国产精品福利导航| 亚洲天堂成人在线观看| 9久久伊人精品综合| 国禁国产you女视频网站| 国产主播在线一区| 日韩最新中文字幕| 一级毛片在线免费视频| 亚洲全网成人资源在线观看| 精品小视频在线观看| 亚洲精品第一页不卡| 国产一区二区三区免费观看| 亚洲无码高清一区| 18禁黄无遮挡免费动漫网站| 中文字幕日韩欧美| 玖玖精品在线| 日本午夜精品一本在线观看 | 亚洲欧美日本国产综合在线| 国产黄色爱视频| 天天综合网色| 网友自拍视频精品区| 欧美综合区自拍亚洲综合天堂| 免费人成在线观看成人片 | 欧美特级AAAAAA视频免费观看| 午夜天堂视频| 国内精品小视频在线| 国产97色在线| 欧美日韩一区二区在线播放| 国产激情在线视频| 天天综合色天天综合网| 毛片久久网站小视频| 婷婷综合在线观看丁香| 日韩成人午夜| 四虎亚洲精品| 亚洲高清中文字幕在线看不卡|