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

SAE9310鋼動(dòng)態(tài)再結(jié)晶臨界條件的研究

2014-03-13 10:33:12黃順喆厲勇王春旭劉憲民田志凌
航空材料學(xué)報(bào) 2014年3期
關(guān)鍵詞:變形

黃順喆, 厲勇, 王春旭, 韓 順, 劉憲民, 田志凌

(1.鋼鐵研究總院特殊鋼研究所,北京100081;2.中國鋼研科技集團(tuán)有限公司,北京100081)

動(dòng)態(tài)再結(jié)晶臨界應(yīng)變通常用作判斷材料在熱變形過程中是否發(fā)生動(dòng)態(tài)再結(jié)晶,對(duì)材料熱變形過程中工藝控制有重要的指導(dǎo)意義,確定材料在熱變形中的臨界應(yīng)變量是建立臨界應(yīng)變預(yù)測模型的關(guān)鍵條件[1]。如果材料的應(yīng)力-應(yīng)變曲線上出現(xiàn)峰值σp,那么可以認(rèn)為發(fā)生了動(dòng)態(tài)再結(jié)晶行為,但研究中把峰值應(yīng)力對(duì)應(yīng)的應(yīng)變?chǔ)舙看做臨界應(yīng)變?chǔ)與是不科學(xué)的,事實(shí)上動(dòng)態(tài)再結(jié)晶在應(yīng)變量未達(dá)到峰值應(yīng)變?chǔ)舙時(shí)就已經(jīng)發(fā)生[2,3]。材料的加工硬化率 (θ= dσ/dε)表征流變應(yīng)力隨應(yīng)變的變化速率,流變曲線是熱變形時(shí)材料微觀組織的外在表現(xiàn),一般很難直接從曲線上來確定發(fā)生動(dòng)態(tài)再結(jié)晶的臨界應(yīng)變,而加工硬化率曲線則能反映材料內(nèi)部組織的變化特征[4,5]。Poliak[2,3]認(rèn)為,材料發(fā)生動(dòng)態(tài)再結(jié)晶軟化時(shí),不管應(yīng)力-應(yīng)變曲線上是否出現(xiàn)應(yīng)力峰,其θ-σ曲線上均呈現(xiàn)拐點(diǎn),并認(rèn)為臨界條件與dθ/dσ-σ曲線上的最大值以及θ-σ曲線上的拐點(diǎn)相對(duì)應(yīng)。很多學(xué)者[1,6~9]應(yīng)用加工硬化理論,在判斷材料發(fā)生動(dòng)態(tài)再結(jié)晶的臨界應(yīng)變方面取得較好的成果。

SAE9310鋼是國內(nèi)外廣泛使用的航空用低成本、高強(qiáng)度滲碳鋼,其強(qiáng)度高韌性好,且淬透性和可焊性良好,多用于制造截面尺寸較大且受重載荷的傳動(dòng)構(gòu)件,如航空領(lǐng)域中的關(guān)鍵齒輪、齒輪軸、主旋翼軸、尾槳軸等傳動(dòng)裝置[10-13]。本文基于加工硬化理論,通過熱模擬實(shí)驗(yàn)獲得的應(yīng)力-應(yīng)變曲線數(shù)據(jù),研究SAE9310鋼在大應(yīng)變條件下的加工硬化率曲線、動(dòng)態(tài)再結(jié)晶的拐點(diǎn)判據(jù)以及動(dòng)態(tài)再結(jié)晶的臨界條件。

1 實(shí)驗(yàn)材料及方法

實(shí)驗(yàn)采用真空感應(yīng)+真空自耗重熔冶煉的SAE9310鋼(以下簡稱實(shí)驗(yàn)鋼),其化學(xué)成分為(質(zhì)量分?jǐn)?shù)/%):C 0.11,Cr 1.25,Ni3.22,Mo 0.13,Mn 0.68,Si0.26,其余為Fe。鋼錠經(jīng)高溫?cái)U(kuò)散退火后,鍛造成材,將試樣加工成尺寸為φ10mm×15mm的圓柱形壓縮試樣后,在Gleeble-3800熱模擬實(shí)驗(yàn)機(jī)上進(jìn)行等溫恒應(yīng)變速率的軸向壓縮實(shí)驗(yàn),測得該鋼的高溫流變曲線。變形溫度為900℃,1000℃,1100℃,1200℃,變形速率為0.01s-1,0.1s-1,1s-1,10s-1,最大真應(yīng)變?yōu)?.9,變形后立刻噴水淬火,以分析其高溫組織。用線切割機(jī)把熱壓縮試樣沿壓縮方向從中間剖開,磨平、拋光后,用飽和苦味酸水溶液+少量洗劑靈溶液腐蝕奧氏體晶界,用光學(xué)顯微鏡觀察顯微組織。

2 實(shí)驗(yàn)結(jié)果與分析

2.1 加工硬化理論

材料在加工過程中的流變應(yīng)力與位錯(cuò)密度的變化存在如下關(guān)系[14]:

式中:σ為流變應(yīng)力(MPa),ε為真應(yīng)變,ρ為位錯(cuò)密度(m-3),M為Taylor因子,μ為剪切模量,θ為加工硬化率,其他為材料常數(shù)。式(1),(2),(3)表明:材料的加工硬化率θ與位錯(cuò)密度ρ有關(guān)。根據(jù)加工硬化理論[15-21],材料的加工硬化率隨應(yīng)力變化的規(guī)律可分為五個(gè)階段,分別是Ⅰ-易滑移階段、Ⅱ-線性硬化階段、Ⅲ-動(dòng)態(tài)回復(fù)硬化階段、Ⅳ-大應(yīng)變硬化階段、Ⅴ-動(dòng)態(tài)再結(jié)晶軟化階段,如圖1所示。

圖1 加工硬化率θ隨應(yīng)力σ變化示意圖[1]Fig.1 Sketchmap of strain hardening rateθand stressσ

本工作重點(diǎn)研究第Ⅴ階段-動(dòng)態(tài)再結(jié)晶軟化階段。在這一階段里,變形進(jìn)行到一定程度后,變形組織內(nèi)部的位錯(cuò)密度達(dá)到一定的臨界值,開始發(fā)生動(dòng)態(tài)再結(jié)晶,此時(shí)材料的θ-σ曲線上呈現(xiàn)拐點(diǎn)特征,拐點(diǎn)處對(duì)應(yīng)的應(yīng)力即是動(dòng)態(tài)再結(jié)晶臨界應(yīng)力σc[20,22]原字符串,且lnθ-lnσ曲線與lnθ-ε曲線上均呈現(xiàn)拐點(diǎn)特征[22],由此通過分析應(yīng)變硬化速率θ與流變應(yīng)力σ曲線上的拐點(diǎn)或者-?θ/?σ-σ曲線上的最小值,即可得到動(dòng)態(tài)再結(jié)晶的臨界條件。

2.2 臨界應(yīng)變的確定

SAE9310鋼在不同變形條件下的高溫流變曲線表現(xiàn)為動(dòng)態(tài)回復(fù)型和動(dòng)態(tài)再結(jié)晶型,如圖2所示。由圖可知,當(dāng)變形溫度為900℃時(shí),實(shí)驗(yàn)鋼在不同應(yīng)變速率下的真應(yīng)力-真應(yīng)變曲線均表現(xiàn)為動(dòng)態(tài)回復(fù)型,隨著變形溫度的提高,實(shí)驗(yàn)鋼的真應(yīng)力-真應(yīng)變曲線可明顯分為兩種類型:動(dòng)態(tài)回復(fù)型和具有單峰特征的動(dòng)態(tài)再結(jié)晶型。材料的應(yīng)力-應(yīng)變曲線上出現(xiàn)峰值意味著材料內(nèi)部組織發(fā)生了動(dòng)態(tài)再結(jié)晶,但應(yīng)力-應(yīng)變曲線并不能直觀的反應(yīng)材料發(fā)生動(dòng)態(tài)再結(jié)晶的臨界條件[1,3]。應(yīng)用加工硬化理論,對(duì)應(yīng)力-應(yīng)變曲線數(shù)據(jù)進(jìn)行處理,得到θ-σ曲線及dθ/dσ-σ曲線,采用加工硬化率拐點(diǎn)判據(jù)來判斷發(fā)生動(dòng)態(tài)再結(jié)晶的臨界應(yīng)力σc和臨界應(yīng)變?chǔ)與。為了確定應(yīng)力-應(yīng)變曲線上各應(yīng)變下的斜率,先對(duì)曲線進(jìn)行擬合,再對(duì)擬合方程求導(dǎo),從而得到各應(yīng)變下的斜率,然后繪制加工硬化率θ與真應(yīng)力σ的關(guān)系曲線,進(jìn)而確定發(fā)生動(dòng)態(tài)再結(jié)晶的臨界條件。

圖2 SAE9310鋼在不同變形條件下的高溫流變曲線Fig.2 True stress-strain curves of SAE9310 steel at different deformation conditions (a)900℃;(b)1000℃;(c)1100℃;(d)1200℃

以變形溫度為1100℃,應(yīng)變速率為0.1s-1時(shí)的真應(yīng)力-真應(yīng)變曲線為例,見圖3所示,對(duì)圖中曲線中應(yīng)力σ從0到峰值應(yīng)力σp段進(jìn)行擬合,得到此變形階段各應(yīng)力對(duì)應(yīng)的加工硬化率θ。求得各應(yīng)變下的加工硬化率θ,得到加工硬化率θ與真應(yīng)力σ關(guān)系曲線,如圖4a所示。對(duì)圖中曲線進(jìn)行3次多項(xiàng)式擬合得到θ-σ方程:

對(duì)方程(4)求導(dǎo)并繪制dθ/dσ-σ曲線,如圖4b所示。當(dāng)d2θ/dσ2=0時(shí)對(duì)應(yīng)的應(yīng)力σ即為臨界應(yīng)力σc,對(duì)應(yīng)的應(yīng)變?chǔ)偶礊榕R界應(yīng)變?chǔ)與,求得在變形條件為 1100℃,應(yīng)變速率為 0.1s-1時(shí) SAE9310鋼的臨界應(yīng)力σc為55.8MPa,對(duì)應(yīng)臨界應(yīng)變?chǔ)與為 0.09。同理得到不同變形條件下SAE9310鋼的加工硬化率θ與應(yīng)力σ的關(guān)系曲線,見圖5所示,為了確定θ-σ曲線中拐點(diǎn)數(shù)值,對(duì)圖中曲線求一階偏導(dǎo),結(jié)果如圖6所示。

圖3 應(yīng)變速率為0.1s-1時(shí)SAE9310鋼的真應(yīng)力-真應(yīng)變曲線Fig.3 True stress-true strain cures of SAE9310 steel deformed at strain rate of 0.1s-1

圖4 變形條件為1100℃,0.1s-1時(shí)的θ-σ關(guān)系曲線(a)及dθ/dσ-σ曲線(b)Fig.4 Curves ofθversusσ(a)and dθ/dσversusσ(b)deformed at1100℃ with 0.1 s-1

根據(jù)圖5、圖6可得到SAE9310鋼在不同應(yīng)變速率和變形溫度下發(fā)生動(dòng)態(tài)再結(jié)晶的臨界應(yīng)力和臨界應(yīng)變,如圖7、圖8所示。由圖可知,當(dāng)變形溫度一定時(shí),隨著應(yīng)變速率的增加,臨界應(yīng)變和臨界應(yīng)力都相應(yīng)提高;而當(dāng)變形速率一定時(shí),提高變形溫度,臨界應(yīng)變和臨界應(yīng)力均降低。

Z參數(shù)(Zener-Hollomon參數(shù))被廣泛用以表征變形溫度及應(yīng)變速率對(duì)熱變形過程的綜合作用。在熱變形過程中有[23]:

臨界應(yīng)變?chǔ)與和峰值應(yīng)變?chǔ)舙與Z參數(shù)的關(guān)系如圖9所示。由圖可知,Z參數(shù)能反應(yīng)臨界應(yīng)變及峰值應(yīng)變與變形條件的變化規(guī)律;臨界應(yīng)變與熱加工參數(shù)Z是密切相關(guān)的,臨界應(yīng)變與峰值應(yīng)變均隨Z參數(shù)的增大而增加,且臨界應(yīng)變與峰值應(yīng)變之間存在一定相關(guān)性,比值約為εc/εp=0.30~0.42。

2.3 組織演變與動(dòng)態(tài)再結(jié)晶狀態(tài)圖

圖10為SAE9310鋼在變形溫度為1100℃、應(yīng)變速率為0.1s-1條件下的加工硬化率與真應(yīng)變的關(guān)系曲線,曲線中所標(biāo)A,B,C,D對(duì)應(yīng)應(yīng)變分別為臨界應(yīng)變?chǔ)與、峰值應(yīng)變?chǔ)舙、最大軟化率處應(yīng)變?chǔ)?和穩(wěn)態(tài)流變應(yīng)力對(duì)應(yīng)應(yīng)變 εs,圖11為對(duì)應(yīng)四點(diǎn)的金相組織。

圖5 不同變形溫度下SAE9310鋼的θ-σ關(guān)系曲線Fig.5 Curves ofθversusσfor SAE9310 steel under different deformation temperatures (a)900℃;(b)1000℃;(c)1100℃;(d)1200℃

圖6 不同變形溫度下SAE9310鋼的dθ/dσ-σ關(guān)系曲線Fig.6 Curves of dθ/dσversusσfor SAE9310 steel under different deformation temperatures (a)900℃;(b)1000℃;(c)1100℃;(d)1200℃

圖7 不同變形條件下SAE9310鋼的σc-˙ε曲線與σc-T曲線Fig.7 Curves ofσc versus˙εandσc versus T for SAE9310 steel under different deformation conditions

圖8 不同變形條件下SAE9310鋼的εc-˙ε曲線與εc-T曲線Fig.8 Curves ofεc versus˙εandεc versus T for SAE9310 steel under different deformation conditions

圖9 Z參數(shù)與εc,εp的關(guān)系曲線Fig.9 Relationship between Z andεc,εp

在變形初期,實(shí)驗(yàn)鋼的奧氏體晶粒被壓扁、拉長,晶界發(fā)生局部遷移,當(dāng)應(yīng)變量達(dá)到臨界應(yīng)變?chǔ)與時(shí),少量新鮮等軸狀動(dòng)態(tài)再結(jié)晶晶粒在變形的原奧氏體晶界處形成(圖11a);當(dāng)應(yīng)變量為峰值應(yīng)變?chǔ)舙時(shí),加工硬化率為0,實(shí)驗(yàn)鋼發(fā)生部分動(dòng)態(tài)再結(jié)晶,變形晶粒與等軸狀再結(jié)晶晶粒共存(圖11b),此后動(dòng)態(tài)再結(jié)晶軟化作用將大于加工硬化作用;當(dāng)應(yīng)變量達(dá)到最大軟化率處應(yīng)變?chǔ)?時(shí),加工硬化降到最小值,動(dòng)態(tài)再結(jié)晶基本完成,金相組織變現(xiàn)為不均勻的再結(jié)晶晶粒(圖11c);隨著變形量的進(jìn)一步增加,當(dāng)位錯(cuò)密度增大產(chǎn)生的加工硬化和再結(jié)晶過程中軟化作用趨于平衡時(shí),流變曲線進(jìn)入穩(wěn)態(tài)階段,此時(shí)組織為充分動(dòng)態(tài)再結(jié)晶并發(fā)生長大的等軸狀?yuàn)W氏體晶粒(圖11d)。

由實(shí)驗(yàn)數(shù)據(jù)繪得SAE9310鋼的動(dòng)態(tài)再結(jié)晶狀態(tài)圖,如圖12所示,其中A區(qū)—加工硬化區(qū),B區(qū)—部分動(dòng)態(tài)再結(jié)晶區(qū),C區(qū)—完全動(dòng)態(tài)再結(jié)晶區(qū),從動(dòng)態(tài)再結(jié)晶狀態(tài)圖可以看出,隨著Z參數(shù)的增加,即隨著變形溫度的降低和應(yīng)變速率的增加,材料發(fā)生動(dòng)態(tài)再結(jié)晶的臨界變形量逐漸增加,也就是說發(fā)生動(dòng)態(tài)再結(jié)晶變得較為困難。

3 結(jié)論

(1)SAE9310鋼的應(yīng)力-應(yīng)變曲線隨變形條件的不同出現(xiàn)兩種類型的曲線特征,即動(dòng)態(tài)再結(jié)晶型和動(dòng)態(tài)回復(fù)型,兩種類型曲線均可采用θ-σ曲線拐點(diǎn)判據(jù)來確定動(dòng)態(tài)再結(jié)晶臨界條件。

圖10 變形條件為1100℃,0.1s-1時(shí)的θ-σ關(guān)系曲線Fig.10 Relationship betweenθandσ deformed at1100℃ with 0.1 s-1

圖11 對(duì)應(yīng)圖10中點(diǎn)A,B,C,D處的金相組織Fig.11 Metallograph corresponding to point A,B,C,D in fig.10(a)A;(b)B;(c)C;(d)D

圖12 SAE9310鋼高溫?zé)嶙冃蔚膭?dòng)態(tài)再結(jié)晶狀態(tài)圖Fig.12 State diagram of DRX for SAE9310 steel under hot temperature deformation

(2)SAE9310鋼發(fā)生動(dòng)態(tài)再結(jié)晶時(shí),θ-σ曲線存在拐點(diǎn),dθ/dσ-σ曲線上出現(xiàn)最大值,利用拐點(diǎn)判據(jù)可以確定9310鋼發(fā)生動(dòng)態(tài)再結(jié)晶的臨界應(yīng)變。

(3)SAE9310鋼發(fā)生動(dòng)態(tài)再結(jié)晶的臨界應(yīng)變?chǔ)與和臨界應(yīng)力σc均隨應(yīng)變速率的增大和變形溫度的降低而增加,臨界應(yīng)變與峰值應(yīng)變之間滿足εc/εp= 0.30~0.42。隨著Z參數(shù)的增加,即隨著變形溫度的降低和應(yīng)變速率的增加,實(shí)驗(yàn)鋼發(fā)生動(dòng)態(tài)再結(jié)晶的臨界變形量逐漸增加,發(fā)生動(dòng)態(tài)再結(jié)晶變得更加困難。

[1]歐陽德來,魯世強(qiáng),崔霞,等.應(yīng)用加工硬化率研究TA15鈦合金β區(qū)變形的動(dòng)態(tài)再結(jié)晶臨界條件[J].航空材料學(xué)報(bào),2010,30(2):17-23. (OUYANG D L,LU SQ,CUIX,et al.Study on critical strains of dynamic recrystallization duringβprocess in TA15 Titanium alloy using working hardening rate[J]. Journal of Aeronautical Materials,2010,30(2):17 -23.)

[2]POLIAK E I,JONASJJ.A one-parameter approach to determining the critical conditions for the initiation of dynamic recrystallization[J].Acta Materialia,1996,44(1):127 -136.

[3]POLIAK E I,JONAS J J.Initiation of dynamic recrystallization in constantstrain rate hotdeformation[J].ISIJInternational,2003,43(5):684-691.

[4]PRASAD G V,GOERDELER M,GOTTSTEIN G.Work hardening model based on multiple dislocation densities[J].Materials Science and Engineering(A),2005,400 -401;231-233.

[5]ROLLETT A D,KOCKS U F.A review of the stages of work hardening[J].Solid State Phenomena,1993,35-36;1-18.

[6]GOTTSTEIN G,BRUNGER E,F(xiàn)ROMMERT M,et al. Prediction of the critical conditions for dynamic recrystallization inmetals[J].Zeitschrift for Metallkunde,2003,94 (5):628-635.

[7]POLIAK E I,JONAS J J.Critical strain for dynamic recrystallization in variable strain rate hot deformation[J]. ISIJ International,2003,43(5):692-700.

[8]NAJAFIZADEH A,JONAS JJ.Predicting the critical stress for initiation of dynamic recrystallization[J].ISIJ International,2006,46(11):1679-1684.

[9]黃光杰,錢寶華,汪凌云,等.AZ31鎂合金初始動(dòng)態(tài)再結(jié)晶的臨界條件研究[J].稀有金屬材料與工程,2007,36 (12):2080-2083. (HUANG G J,QIAN B H,WANG L Y,et al.Study on the critical conditions for initial dynamic recrystallization of AZ31 magnesium alloy[J].Rare Metal Materials and Engineering,2007,36(12):2080-2083.

[10]周鳳云,李熙章.高強(qiáng)度滲碳鋼制構(gòu)件的斷裂分析[J].理化檢驗(yàn)-物理分冊(cè),2003,39(4):206-209. (ZHOU F Y,LIX Z.Fracture analysis of the components made of high strength and carburizing steel[J].Physical Testing and Chemical Analysis Part A:Physical Testing,2003,39(4):206-209.)

[11]吳秋平,王春旭,劉憲民,等.回火溫度對(duì)9310鋼力學(xué)性能及組織的影響[J].熱加工工藝,2012,41(6): 179-183. (WU Q P,WANG C X,LIU X M,et al.Influence of tempering temperature on mechanical properties and microstructure of 9310 steel[J].Hot Working Technology,2012,41(6):179-183.)

[12]黃順喆,厲勇,王春旭,等.9310鋼的奧氏體晶粒長大規(guī)律研究[J].熱加工工藝,2010,39(18):31-33. (HUANG SZ,LIY,WANG C X,et al.Study on growth law of austenite grain in 9310 steel[J].Hot Working Technology,2010,39(18):31-33.)

[13]厲勇,王春旭,劉憲民,等.SAE9310鋼奧氏體的冷卻轉(zhuǎn)變行為[J].機(jī)械工程材料,2010,34(5):12-15. (LI Y,WANG C X,LIU X M,et al.Microstructure transformation behavior of austenite cooling for SAE9310 steel[J].Mater Mech.Eng,2010,34(5):12-15.)

[14]ESTRIN Y.Dislocation theory based constitutive modeling:foundations and applications[J].Journal of Materials Processing Technology,1998,80/81;33-39.

[15]UUN T,ZEHETBAUER M.Stage IV work hardening in cell forming materials,part II:a new mechanism[J]. Scripta Materialia,1996,35(12):1467-1473.

[16]GOTTSTEIN G,F(xiàn)ROMMERTM,GOERDELER M,etal. Prediction of the critical conditions for dynamic recrystallization in the austenitic steel 800H[J].Materials Science and Engineering(A),2004,387/389;604-608.

[17]NABARRO F,BASINSKIZ S,HOLTD B.The plasticity of pure single crystals[J].Advances in Physics,1964,13 (50):193-323.

[18]MECKING H,KOCKS U F.Kinetics of flow and strain hardening[J].Acta Metall,1981,29(11):1865 -1875.

[19]KOCKSU F.Laws ofwork hardening and low-temperature creep[J].JEng Mater Technol,1976,98(1):76-85.

[20]GIL S J,VAN H P,AERNOUDT E.Large strain work hardening and textures[J].Progressing Materials Science,1980,25(2-4):69-134.

[21]MECKING H.Dislocation modeling of physical systems[M].Oxford:Pergamon Press,1981:197-211.

[22]POLIAk E I,JONAS JJ.Initiation of dynamic recrystallization in constant strain rate hot deformation[J].ISIJInternational,2003,43(5):684-691.

[23]MEDINA SF,HERNANDEZ C A.General expression of Zener-Hollomon parameter as a function of the chemical composition of low alloy and micro-ally steels[J].Acta Materialia,1996,44(1):137-1548.

猜你喜歡
變形
變形記
談詩的變形
中華詩詞(2020年1期)2020-09-21 09:24:52
柯西不等式的變形及應(yīng)用
“變形記”教你變形
不會(huì)變形的云
“我”的變形計(jì)
會(huì)變形的折紙
童話世界(2018年14期)2018-05-29 00:48:08
變形巧算
例談拼圖與整式變形
會(huì)變形的餅
主站蜘蛛池模板: 亚洲第一中文字幕| 亚国产欧美在线人成| 欧美一区二区三区香蕉视| 1级黄色毛片| 国产欧美专区在线观看| 日韩高清在线观看不卡一区二区 | 国产爽爽视频| 国产成人AV大片大片在线播放 | 一级毛片不卡片免费观看| 福利视频一区| 久久综合色天堂av| 国产欧美日韩综合一区在线播放| 日本日韩欧美| 在线欧美国产| 亚洲VA中文字幕| 97国产成人无码精品久久久| 免费jjzz在在线播放国产| 欧美成人免费| 精品福利国产| 国产视频入口| 日本一区二区不卡视频| 亚洲日韩在线满18点击进入| 三上悠亚精品二区在线观看| 无码免费视频| 网友自拍视频精品区| 亚洲制服丝袜第一页| 国产美女主播一级成人毛片| 东京热av无码电影一区二区| 亚洲成人一区二区| 91免费国产高清观看| 日韩欧美网址| 在线看片中文字幕| 国产91视频免费观看| 亚洲天堂网在线观看视频| 国产在线专区| 一本大道香蕉高清久久| 日韩精品毛片人妻AV不卡| 国产成人精品三级| 日韩免费毛片| 亚洲一区波多野结衣二区三区| 日本午夜影院| 性69交片免费看| 77777亚洲午夜久久多人| 朝桐光一区二区| www精品久久| 国产精品区网红主播在线观看| 大乳丰满人妻中文字幕日本| 国产亚洲精品自在久久不卡| 免费视频在线2021入口| 久热这里只有精品6| 久久久成年黄色视频| a在线观看免费| 亚洲人成在线精品| 99re这里只有国产中文精品国产精品| 亚洲大尺码专区影院| 性视频久久| 国产欧美专区在线观看| 国产波多野结衣中文在线播放 | 日a本亚洲中文在线观看| 2020国产精品视频| 成人久久精品一区二区三区| 色偷偷男人的天堂亚洲av| 人妻无码中文字幕第一区| 国产91熟女高潮一区二区| 亚洲国语自产一区第二页| 很黄的网站在线观看| 成人精品视频一区二区在线| 999在线免费视频| 亚洲资源在线视频| 无码日韩人妻精品久久蜜桃| 亚洲A∨无码精品午夜在线观看| 一级毛片在线免费看| 亚洲第一综合天堂另类专| 亚洲福利视频一区二区| 久久久无码人妻精品无码| 久久综合伊人 六十路| 亚洲国产日韩一区| 凹凸精品免费精品视频| 精品无码国产自产野外拍在线| 伊人久综合| 国产美女91视频| 亚洲人成在线精品|