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

        ?

        擠壓態(tài)AZ81鎂合金的熱壓縮變形行為

        2011-11-23 03:03:40曹鳳紅龍思遠COLE
        中國有色金屬學報 2011年6期
        關(guān)鍵詞:再結(jié)晶鎂合金軟化

        曹鳳紅,龍思遠,G. S. COLE

        (1. 重慶大學 材料科學與工程學院,重慶 400030;2. 成都理工大學 工程技術(shù)學院,樂山 614007;3. 美國汽車通用公司, 美國,48090)

        擠壓態(tài)AZ81鎂合金的熱壓縮變形行為

        曹鳳紅1,2,龍思遠1,G. S. COLE3

        (1. 重慶大學 材料科學與工程學院,重慶 400030;2. 成都理工大學 工程技術(shù)學院,樂山 614007;3. 美國汽車通用公司, 美國,48090)

        在溫度為320~440 ℃、應(yīng)變率為0.001~1 s?1的變形條件下,采用Gleeble?1500熱模擬機對擠壓態(tài)AZ81鎂合金的熱壓縮變形行為進行研究。結(jié)果表明:擠壓態(tài) AZ81鎂合金的流變應(yīng)力隨變形溫度的升高而降低,隨應(yīng)變率的升高而升高,且隨應(yīng)變的增加,流動應(yīng)力很快達到峰值,然后逐漸降低并趨于穩(wěn)定。為評價擠壓態(tài)AZ81鎂合金在熱模壓成形過程中流動應(yīng)力,結(jié)合Arrhenius方程并引入Zener-Hollomon參數(shù),對流動應(yīng)力做出相應(yīng)的修正,根據(jù)修正后的流動應(yīng)力構(gòu)建擠壓態(tài)AZ81鎂合金流變應(yīng)力高溫變形本構(gòu)模型。模擬結(jié)果表明:該模型的應(yīng)力預(yù)測值與試驗值吻合較好,計算精度較高,為后續(xù)的模壓近/凈終成形工藝參數(shù)的制定提供一定的理論參考。

        擠壓態(tài)AZ81鎂合金;熱壓縮模擬;流變應(yīng)力;本構(gòu)模型

        鎂的密度小,是最輕的金屬結(jié)構(gòu)材料,與鑄造鎂合金相比,變形鎂合金在汽車、航空航天以及電子通訊等許多領(lǐng)域有更廣闊的應(yīng)用更前景[1]。目前,工業(yè)發(fā)達國家已大力著手高性能變形鎂合金加工型材的研究與開發(fā)[2]。但鎂合金具有密排六方結(jié)構(gòu),滑移系較少,低溫塑性較差,加工難度大,且隨著溫度的升高,滑移系增多,塑性大大提高[4?5],龍思遠和曹韓學[6]指出鑄鍛復合成形方法獲得產(chǎn)品其力學性能得到了一定程度的改善,但因變形程度有限,不能滿足高強韌鎂合金承力結(jié)構(gòu)件對用大變形加工提升材料機械性能的要求。為提高鎂合金的塑性變形性能,提出了鑄?擠?鍛復合成形方法[7],通過擠壓預(yù)成形塑性變形機制,使鎂合金具有更大的塑性變形能力,通過模壓一次成形工藝,使鎂合金的綜合力學性能得到明顯改善。研究認為[8],擠壓預(yù)成形坯Mg-Al-Zn系合金及其模壓坯料Al含量為8.0%時可獲得令人滿意的強度和塑性指標,因此擠壓態(tài)AZ81鎂合金具有制備高性能結(jié)構(gòu)件的潛力。目前,對板材AZ81鎂合金的變形性能,尤其高溫變形性能的研究未見報道。在模壓成形工藝的研究中,流變應(yīng)力常作為已知的基本參數(shù),是合理進行數(shù)值模擬、制定成形工藝及選擇變形設(shè)備的理論依據(jù)。

        鑒于此,本文作者在溫度320~440 ℃及應(yīng)變速率為 0.001~1 s?1的變形條件下,采用 Gleeble?1500D熱模擬試驗機對 AZ81鎂合金進行單向壓縮試驗。利用Arrheniius方程并引入Zener-Hollomon參數(shù)建立 AZ81鎂合金的流動應(yīng)力模型,同時進行修正,為模壓成形工藝的制定提供一定的理論實驗參考。

        1 實驗

        試驗用合金成分為,其合金成分為Mg-Al(8.15%)-Zn(0.82%)-Mn(0.27%)(質(zhì)量分數(shù)),其中 w(Fe)≤0.005,w(Si)≤0.019,w(Ni)≤0.000 7,w(Cu)≤0.001 6,所用原料均采用I級單質(zhì)態(tài),在坩堝電阻爐中熔煉合金,同時采用氣體及少量熔劑進行保護,連鑄成d 112(±0.2)mm×250 mm半連續(xù)鑄錠,連鑄過程中采用SO2(1%,體積分數(shù))和CO2(99%,體積分數(shù))混合氣體進行保護,澆注后空冷至室溫,然后在XJ?800臥室擠壓機上進行擠壓,擠壓比16,擠壓溫度為400 ℃,擠壓速度為13 mm/s,得到所需擠壓型材,沿垂直型材擠壓方向切割制備成12 mm×6 mm×8 mm壓縮試樣,如圖1所示。

        圖1 壓縮試樣取樣示意圖Fig.1 Sampling diagram of compression specimen

        高溫壓縮試驗在Gleeble?1500材料熱模擬試驗機上進行恒應(yīng)變速率等溫熱壓縮實驗,壓縮時在試樣兩端與壓頭接觸的ICI墊上鈕片,保證試樣在加熱和壓縮過程中溫度均勻減少摩擦。試驗溫度為320、360、400和440 ℃,熱壓縮實驗開始前,試樣以10 ℃/s的加熱速度升溫至預(yù)定溫度,并保溫 3 min;應(yīng)變速率為0.001、0.01、0.1和1 s?1;變形程度均為60%,壓縮變形后的試樣立即進行水淬處理,以保留形變后的組織;由Geeble?1500計算機系統(tǒng)自動采集應(yīng)力、應(yīng)變、壓力、溫度、位移及時間等數(shù)據(jù),繪制真應(yīng)力—真應(yīng)變曲線。

        2 結(jié)果與分析

        擠壓態(tài)AZ81鎂合金熱壓縮變形的真應(yīng)力—真應(yīng)變實驗曲線如圖2所示,其曲線屬于典型的單峰動態(tài)再結(jié)晶,同時也說明AZ81鎂合金在熱變形過程中易于發(fā)生動態(tài)再結(jié)晶。

        熱壓縮變形過程是加工硬化與動態(tài)再結(jié)晶軟化相互競爭的一個動態(tài)平衡過程。擠壓態(tài)AZ81鎂合金在熱壓縮時流變應(yīng)力的變化規(guī)律是(見圖2):微應(yīng)變階段曲線斜率較大,應(yīng)力上升快,主要是由于加工硬化占主導作用,合金中發(fā)生了部分動態(tài)回復或動態(tài)再結(jié)晶,流變應(yīng)力隨真應(yīng)變的增加迅速上升;且隨著變形量的增加,由于位錯密度不斷提高,使得動態(tài)回復和動態(tài)再結(jié)晶加快,軟化作用增強,加工硬化逐漸被抵消,表現(xiàn)為曲線斜率逐漸減小,應(yīng)力上升減緩,當流變應(yīng)力達到峰值時,加工硬化和動態(tài)再結(jié)晶軟化達到平衡;出現(xiàn)峰值后,隨著變形的繼續(xù)進行,動態(tài)再結(jié)晶繼續(xù)發(fā)展,軟化作用繼續(xù)增強,使得流變應(yīng)力逐漸下降,當達到一定的真應(yīng)變量后,其流變應(yīng)力基本不變。同時還可以看出,峰值應(yīng)力與所對應(yīng)的穩(wěn)態(tài)流動應(yīng)力之間的差值隨溫度的升高而減小,在440 ℃時,峰值流變應(yīng)力與穩(wěn)態(tài)流動應(yīng)力基本變化不明顯;當應(yīng)變率一定時,其峰值流變應(yīng)力隨變形溫度的升高而逐漸下降,此時非基面滑移的臨界切應(yīng)力下降,導致合金變形抗力降低,動態(tài)回復和動態(tài)再結(jié)晶更易發(fā)生,因此,峰值應(yīng)力和峰值應(yīng)變均減小[9?12]。

        當變形溫度一定時,峰值應(yīng)力隨應(yīng)變率的增加,加工硬化作用也增大,其硬化速率超過了動態(tài)再結(jié)晶的軟化速率,再結(jié)晶軟化與加工硬化兩者作用效果相互平衡時的峰值應(yīng)力和峰值應(yīng)變均增大;隨應(yīng)變速率的增大和溫度的降低,進入穩(wěn)態(tài)變形時對應(yīng)的真應(yīng)變值增大。

        圖2 不同變形條件下擠壓態(tài)AZ81鎂合金的流動應(yīng)力—應(yīng)變曲線Fig.2 Flow stress—strain curves of AZ81 alloy under different deformation conditions: (a) 320 ℃; (b) 360 ℃; (c) 400 ℃; (d) 440 ℃

        3 本構(gòu)模型的建立

        由于材料在塑性加工過程中的動態(tài)響應(yīng)是材料內(nèi)部組織演化過程引起的硬化和軟化過程中綜合作用的結(jié)果,故本構(gòu)關(guān)系是高度非線性的,不存在普遍使用的構(gòu)造方法。

        熱加工變形是一個熱激活過程,對外加應(yīng)力、應(yīng)變速率和溫度極其敏感,是受熱激活的控制。通過對不同熱加工數(shù)據(jù)的研究得出在任何應(yīng)變或穩(wěn)態(tài)下的高溫流變應(yīng)力σ強烈地取決于變形溫度T和應(yīng)變速率ε˙,SELLARS 和 MCTEGRAT[13]和 WHITTENBERGER[14]提出一種包含變形激活能Q和溫度T的雙曲線正弦形式修正Arrhenius關(guān)系來描述這種熱激活行為,通常可用下列3種形式進行描述[15]:1) 在低應(yīng)力水平下,流變應(yīng)力峰值σp和Z之間的關(guān)系可用指數(shù)關(guān)系描述,見式(1);2) 而在高應(yīng)力水平下可用冪指數(shù)關(guān)系描述,見式(2);3) 在整個應(yīng)力水平下可用雙曲函數(shù)關(guān)系描述,見式(3)。

        式中:A1、A2、A、n1、α和β均為常數(shù),1/nβα=;Q為變形激活能,J/mol,又稱動態(tài)軟化激活能,反映材料變形材料的難易程度,也是材料在熱變形過程中重要的力學性能參數(shù);R為摩爾氣體常數(shù),8.314 J/(mol·K);T 為絕對溫度。

        目前所報道的關(guān)于鎂合金應(yīng)力方程的描述有多種。TAKUDA等[16]采用指數(shù)關(guān)系對AZ31和AZ91鎂合金熱變形鎂合金流變應(yīng)力進行描述,GALIYE等[17?18]采用冪指數(shù)關(guān)系對ZK60鎂合金熱變形鎂合金流變應(yīng)力進行描述,BARNETT[19]采用雙曲函數(shù)關(guān)系對AZ31鎂合金熱變形鎂合金流變應(yīng)力進行描述,結(jié)合本實驗情況將式(1)~(3)分別兩邊取自然對數(shù),

        對式(4)~(6)兩邊分別求偏導,得

        其中:

        圖3所示為擠壓態(tài)AZ81鎂合金不同變形條件下的峰值應(yīng)力、應(yīng)變率和溫度的關(guān)系。由曲線 l n ε ˙—lnσp和 ln ε˙—σp看出,峰值應(yīng)力與應(yīng)變率基本呈線性關(guān)系,并由分析可知,n1=6.303 1,β=0.027 1,經(jīng)優(yōu)化處理α值為0.004 3。

        通過對lnσp—T?1和 ln [sinh(ασp)]—T?1及圖4 ln ε ˙ 與 ln [sinh(ασp)]的曲線分析可知,n=4.486 5,b=4.883 8。且從圖4可以看出,不同變形條件下的峰值應(yīng)力與應(yīng)變率及溫度均為較為近似的平行直線,相關(guān)系數(shù)均在0.984 3以上,說明擠壓態(tài)AZ81鎂合金在熱壓縮實驗條件下的變形激活能幾乎保持不變,其變形激活能為Q=182.17 kJ/mol,A=8.905 7×1012。

        將上述分析得到的參數(shù)帶入式(3),則擠壓態(tài)

        圖3 不同條件下AZ81鎂合金峰值應(yīng)力與應(yīng)變速率和溫度的關(guān)系:Fig.3 Relationship between peak stress, strain-rate and temperature under different conditions: (a) l n ε ˙ — ln σp; (b) l nε˙ —σp; (c)lnσp —T ?1; (d) l n[sinh(ασp)]—T?1

        圖4 l n ε ˙ 與 ln [sinh(ασ) ]關(guān)系Fig.4 Relationship between l n ε ˙ and l n[sinh(ασ)]

        圖6 計算值與實測值之間的比較Fig.6 Comparison between predicted value and measured value

        AZ81合金高溫壓縮下的流變應(yīng)力行為的本構(gòu)方程為

        在式(2)中通過引入?yún)?shù)Z,并采用最小系數(shù)二乘法進行線性擬合得到擠壓態(tài)AZ81鎂合金流動應(yīng)力模型,

        圖5所示為擠壓態(tài)AZ81鎂合金壓縮變形參數(shù)Z與峰值應(yīng)力σp的關(guān)系。由圖5的擬合結(jié)果分析看出,參數(shù)lnZ與峰值應(yīng)力σp之間線性關(guān)系很明顯,其線性相關(guān)系數(shù)為0.976 8,實測值與預(yù)測值之間的數(shù)據(jù)吻合的較好(見圖6),其線性相關(guān)系數(shù)均大于0.987 6,因此本回歸模型的相關(guān)性高度顯著。

        圖 5 擠壓態(tài)AZ81鎂合金壓縮變形 Z參數(shù)與峰值應(yīng)力σp的關(guān)系Fig.5 Relationship between compression deformation parameter Z and flow stress σp of as-extruded AZ81 alloy

        4 結(jié)論

        1) 擠壓態(tài)AZ81鎂合金高溫壓縮變形時流變應(yīng)力取決于變形溫度與變形速率,當應(yīng)變率一定時,流變應(yīng)力隨變形溫度的升高而降低,當變形溫度一定時,隨應(yīng)變率的升高而升高。且在變形初期,由于加工硬化占主導,流變應(yīng)力迅速達到峰值;此外在高溫和低應(yīng)變率下,流變應(yīng)力相對比較穩(wěn)定,沒有明顯的峰值;低溫和高溫應(yīng)變速率下,流變應(yīng)力持續(xù)增加到峰值后呈現(xiàn)流變軟化趨勢。

        2) 擠壓態(tài)AZ81鎂合金高溫塑性變形存在熱激活過程。在實驗條件下,AZ81鎂合金熱壓縮變形時的應(yīng)力指數(shù)為4.486 5,其變形激活能為182.17 kJ/mol,其 高 溫 壓 縮 流 變 應(yīng) 力 模 型 為 ε˙= 8 .905 7×1012[sinh(0.0043σ)]4.4865exp[?182170/(RT)]。

        3) 在應(yīng)變速率為 0.001~1 s?1和變形溫度為320~440 ℃條件下,擠壓態(tài)AZ81鎂合金熱壓縮變形的流變應(yīng)力模型,可以用Zener-Hollomon參數(shù)的雙曲函數(shù)形式進行描述,該模型為鎂合金進一步的模壓塑性成形工藝參數(shù)的制定提供較為科學的實驗理論依據(jù)。

        REFERENCES

        [1] 耿義海, 李保成, 張曉芳, 張治民. AZ80變形鎂合金高溫變形流變應(yīng)力分析[J]. 熱加工工藝, 2008, 37(2): 46?49.GENG Yi-hai, LI Bao-cheng, ZHANG Xiao-fang, ZHANG Zhi-min. Analysis on flow stress of AZ80 magnesium alloy during high temperature deformation[J]. Hot Working Technology, 2008, 37(2): 46?49.

        [2] 欒 娜, 李落星, 李光耀, 鐘志華. AZ80鎂合金的高溫熱壓縮變形行為[J]. 中國有色金屬學報, 2007, 17(10): 1678?1684.LUAN Na, LI Luo-xing, LI Guang-yao, ZHONG Zhi-hua. Hot compression deformation behaviors of AZ80 magnesium alloy at elevated temperature[J]. The Chinese Journal of Nonferrous Metals, 2007, 17(10): 1678?1684.

        [3] 吳立鴻, 關(guān)紹康, 王利國, 劉 俊. 鍛造鎂合金及影響鍛造成形的幾個關(guān)鍵因素[J]. 鍛壓技術(shù), 2006, 31(4): 7?11.WU Li-hong, GUAN Shao-kang, WANG Li-guo, LIU Jun.Wrought magnesium alloys and several key factors affecting the forging forming[J]. Forging & Stamping Technology, 2006,31(4): 7?11.

        [4] ZHANG P, LINDEMANN J. Influence of shot peening on high cycle fatigue properties of the high-strength wrought magnesium alloy AZ80[J]. Scripta Materialia, 2005, 52: 485?490.

        [5] AGHION E, ELIEZER D. Magnesium and magnesium alloys[M]. Israel: Avi Moshe Bregman Ltd, 2004.

        [6] 龍思遠, 曹韓學. 鎂合金復合成形方法. 中國,200410040086.7[P]. 2005?03?16.LONG Si-yuan, CAO Han-xue. Magnesium alloy composite forming method. CN 200410040086.7[P]. 2005?03?16.

        [7] 龍思遠, 曹鳳紅, 廖慧敏. 鎂合金復合成形方法. 中國,200810069225.7[P]. 2008?01?10.LONG Si-yuan, CAO Feng-hong, LIAO Hui-min. Magnesium alloy composite forming method. CN 200810069225.7[P].2008?01?10.

        [8] 杜 勇. 熱擠壓和Al含量對AZ系鎂合金組織和性能影響研究[D]. 重慶: 重慶大學,2009: 38?43.DU Yong. Influence of hot extrusion and Al on microstructure and properties of AZ based magnesium alloys[D]. Chongqing:Chongqing University, 2009: 38?43.

        [9] 郭 強, 嚴紅革, 陳振華, 張 輝. AZ31鎂合金高溫熱壓縮變形特性[J]. 中國有色金屬學報, 2005, 15(6): 900?906.GUO Qiang, YAN Hong-ge, CHEN Zhen-hua, ZHANG Hui.Hot compression deformation behavior of AZ31 magnesium alloy at elevated temperature[J]. The Chinese Journal of Nonferrous Metals, 2005, 15(6): 900?906.

        [10] 劉曉霏, 嚴 巍, 陳國學. AZ31B鎂合金塑性變形動態(tài)再結(jié)晶的實驗研究[J]. 塑性工程學報, 2005, 12(3): 10?13.LIU Xiao-fei, YAN Wei, CHEN Guo-xue. The experimental study on recrystallization of magnesium alloy AZ31B during plastic defromation[J]. Journal of Plasticity Engineering, 2005,12(3): 10?13.

        [11] 楊亞琴, 李保成. 變形鎂合金熱壓縮變形流變應(yīng)力研究[J].中北大學學報, 2007, 28: 145?148.YANG Ya-qin, LI Bao-cheng. Flow stress of wrought magnesium alloys during hot compression deformation[J].Journal of North University of China, 2007, 28: 145?148.

        [12] 黃光杰, 趙國丹. AZ31鎂合金熱變形規(guī)律的研究[J]. 重慶工學院學報, 2006, 20(2): 61?64.HUANG Guang-jie, ZHAO Guo-dan. Deformation behavior for AZ31 magnesium alloy in thermomechanical processes[J].Journal of Chongqing Institute of Technology, 2006, 20(2):61?64.

        [13] SELLARS C M, MCTEGART W J. On the mechanism of hot deformation[J]. Acta Metallurgica, 1966, 14(9): 1136?1138.

        [14] WHITTENBERGER D J. The influence of grain size and composition on slow plastic flow in FeAl between 1100 and 1 400 K[J]. Mater Sci Eng A, 1986, 77: 103?113.

        [15] MCQUEEN H J. Metal forming: Industrial, mechanical,computational and microstructure[J]. Journal of Materials Processing Technology, 1993, 37(1/4): 3?36.

        [16] TAKUDA H, FUJIMOTO H, HATTA N. Modeling on flow stress of Mg-Al-Zn alloys at elevated temperatures[J]. Journal of Materials Processing Technology, 1998, 80/81(4): 513?516.

        [17] GALIYEV A, KAIBYSHEV R, SAKAI T. Continuous dynamic recrystallization in magnesium alloy[J]. Materials Science Forum, 2003, 419/422(1): 509?514.

        [18] GALIYEV A, KAIBYSHEV R, GOTTSTEIN G. Correlation of plastic deformation and dynamic recrystallization in magnesium alloy ZK60[J]. Acta Materialia, 2001, 49(7): 1199?1207.

        [19] BARNETT M R. Influence of deformation conditions and texture on the high temperature flow stress of magnesium AZ31[J]. Journal of Light Metals, 2001, 1(3): 167?177.

        Hot compression deformation behaviors of as-extruded AZ81 magnesium alloy at elevated temperature

        CAO Feng-hong1,2, LONG Si-yuan1, G. S. COLE3
        (1. College of Material Science and Engineering, Chongqing University, Chongqing 400030, China;2. College of Engineering and Technical, Chengdu University of Technology, Leshan 614007, China;3. General Motors Corporation, USA, 48090)

        The hot compression tests of AZ81 magnesium alloy were performed on Gleeble?1500 at deformation temperature of 320?440 ℃ and strain rate of 0.01?10 s?1. The results show that the relationship between stress and strain is affected obviously by the strain rate and deformation temperature. The flow stress becomes higher with increasing strain rate at constant temperature and lower with increasing deformation temperature at constant strain rate. In order to evaluate the hot-pressing process, a flow stress model was constituted based on Arrhenius equation and temperature-compensated strain rate factor, the Zener-Hollomon parameter, and modified further by considering the effect of strain on stress. The results show that the modified model agrees well with the experimental stress value under the experimental condition what offer a experimental theoretical basis for close-pressing net-shape.

        as-extruded AZ81 alloy; hot compression simulation; flow stress; constitutive model

        TG146.2; TG316.1

        A

        1004-0609(2011)06-1214-06

        重慶市科技攻關(guān)項目(CSTC2007AA4008)

        2010-06-28;

        2010-09-28

        曹鳳紅,講師,博士;電話:15348137697;E-mail: caofh2004@yahoo.com.cn

        (編輯 李艷紅)

        猜你喜歡
        再結(jié)晶鎂合金軟化
        牡丹皮軟化切制工藝的優(yōu)化
        中成藥(2018年10期)2018-10-26 03:41:30
        先進鎂合金技術(shù)與應(yīng)用高峰論壇在滬召開
        AZ91D鎂合金磷酸鹽轉(zhuǎn)化膜的制備與表征
        常化工藝對低溫Hi-B鋼初次及二次再結(jié)晶的影響
        上海金屬(2016年3期)2016-11-23 05:19:38
        鎂合金的化學表面處理探討
        軟骨延遲增強磁共振成像診斷早期髕骨軟化癥
        鑄態(tài)30Cr2Ni4MoV鋼動態(tài)再結(jié)晶行為研究
        大型鑄鍛件(2015年1期)2016-01-12 06:32:58
        髕骨軟化癥的研究進展
        我國鎂及鎂合金相關(guān)項目獲國際標準立項
        Cu元素對7XXX 系列鋁合金再結(jié)晶的影響
        上海金屬(2014年3期)2014-12-19 13:09:04
        不卡日韩av在线播放| 96精品免费视频大全| 亚洲精品二区三区在线观看| 日本免费观看视频一区二区| 亚洲精品第一国产综合精品| 日韩精品无码免费专区网站| 久久国产成人亚洲精品影院老金| 免费av网址一区二区| 国产精品专区第一页天堂2019| 国产精品久久777777| 六月婷婷国产精品综合| 久久精品韩国日本国产| 日韩精品视频免费在线观看网站| 久久99国产精品久久99果冻传媒| 精品国产精品久久一区免费式| 丰满人妻AV无码一区二区三区| 国产免费成人自拍视频| 亚洲人成影院在线无码按摩店| 三年片在线观看免费大全电影| 久久亚洲精品成人综合| 极品美女调教喷水网站| 人妻体内射精一区二区三四| 欧美情侣性视频| 中国老太老肥熟女视频| 亚洲成人av在线蜜桃| 少妇aaa级久久久无码精品片| 国产在线91观看免费观看| 亚洲国产精品二区三区| 国产亚洲精品综合一区| 日本熟妇色xxxxx欧美老妇| 中文字幕无码免费久久99| 中文字幕熟女激情50路| 亚洲国产精彩中文乱码av| 久久久久久久99精品国产片| 日本经典中文字幕人妻| 日本一级特黄aa大片| 亚洲国产精品日韩av专区| avtt一区| 国产三级精品av在线| 欧美内射深喉中文字幕| 1精品啪国产在线观看免费牛牛|