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

0Cr16Ni5Mo低碳馬氏體不銹鋼的熱變形行為及其熱加工圖

2016-08-16 02:06:56袁武華龔雪輝孫永慶梁劍雄
材料工程 2016年5期
關鍵詞:不銹鋼變形

袁武華,龔雪輝,,孫永慶,梁劍雄

(1湖南大學 材料科學與工程學院,長沙 410082;2 鋼鐵研究總院 特殊鋼研究所,北京 100081)

?

0Cr16Ni5Mo低碳馬氏體不銹鋼的熱變形行為及其熱加工圖

袁武華1,龔雪輝1,2,孫永慶2,梁劍雄2

(1湖南大學 材料科學與工程學院,長沙 410082;2 鋼鐵研究總院 特殊鋼研究所,北京 100081)

在Gleeble-3800 熱模擬試驗機上進行高溫壓縮實驗, 研究 0Cr16Ni5Mo低碳馬氏體不銹鋼在變形溫度為900~1150℃、應變速率為0.01~10s-1條件下的熱變形行為。采用雙曲正弦模型確定了該材料的熱變形參數隨應變量的變化規律,建立了相應的熱變形本構方程。根據動態材料模型建立并分析了其熱加工圖,同時觀察了變形組織。結果表明:在熱壓縮過程中, 流變應力隨變形溫度的升高而降低,隨應變速率的升高而增加,變形條件對材料的組織結構有較大影響。材料熱變形參數與應變量之間可采用四次函數關系式表示,并且具有很好的相關性,獲得了該材料的最佳熱變形工藝參數范圍為:變形溫度980~1150℃,應變速率0.01~0.2s-1。

0Cr16Ni5Mo馬氏體不銹鋼;流變應力;本構模型;動態再結晶;熱加工圖

低碳馬氏體不銹鋼是在傳統馬氏體不銹鋼的基礎上通過減少碳、添加鎳和鉬來提高其綜合性能的新型不銹鋼[1,2]。0Cr16Ni5Mo不銹鋼屬于此類低碳馬氏體不銹鋼,因其具有良好的綜合力學性能以及耐沖蝕、耐腐蝕性能,在水輪機組件、石油輸送管線和航空航天等領域得到了廣泛的應用。由于該鋼中含有較高的鉻含量,因此其組織中常含有一定量的δ-鐵素體[3]。 在高溫鍛造時該鋼為奧氏體-鐵素體的兩相組織,如果鍛造不當,極易在兩相界面上出現微裂紋,造成鍛件開裂,因此研究0Cr16Ni5Mo不銹鋼的熱加工性能可為實際生產提供重要依據,而目前對該鋼的熱變形行為及熱加工圖的研究報道較少。同時,當前在建立合金的變形本構關系時大多是以峰值應力或者穩態應力為基礎建立的,并未考慮塑性應變量的影響[4-7]。事實上,由此建立的模型在描述流變應力隨應變的變化特征時有一定的局限性,這是因為累積塑性應變也是決定材料顯微組織演變的主要參數。一些研究者的研究表明[8-10]采用考慮應變量相關的雙曲正弦模型能更好地反映合金在熱變形過程中流變應力的變化規律?;诖耍竟ぷ髟贕leeble-3800熱模擬試驗機上進行單道次熱壓縮實驗,研究變形溫度和變形速率對0Cr16Ni5Mo低碳馬氏體不銹鋼流變應力的影響,計算了不同應變量下的熱變形參數,并建立了相應的本構模型,同時根據動態材料模型建立熱加工圖。

1 實驗材料與方法

實驗材料為0Cr16Ni5Mo低碳馬氏體不銹鋼,采用真空感應+電渣重熔冶煉,鑄錠在1150℃開坯鍛造成φ100mm的圓棒。材料的化學成分(質量分數/%)為:C 0.040,Cr 15.75,Ni 5.19,Mo 0.66,Si 0.30,Mn 0.77,P 0.0064,S 0.009,其余為Fe。為使初始晶粒大小均勻,首先將棒材進行1100℃/1h固溶處理,然后從固溶后棒材上取樣,試樣尺寸為φ8mm×12mm, 在Gleeble-3800熱模擬試驗機上進行壓縮實驗。圖1為試樣固溶處理后的晶粒形貌。進行熱壓縮實驗時,先以20℃/s的升溫速率將試樣加熱到預定變形溫度后,保溫3min,使試樣溫度均勻, 再進行等溫壓縮變形, 變形結束后立即水冷,以保留試樣的高溫變形組織。用線切割機將變形后的試樣沿軸向切開,經研磨、拋光后采用KMnO4+H2SO4+H2O腐蝕顯示其奧氏體晶粒。實驗變形溫度分別為900,950,1000,1050,1100,1150℃,應變速率分別為0.01,0.1,1,10s-1,總壓下量為60%(真應變為0.9)。

圖1 試樣固溶熱處理后顯微組織Fig.1 Microstructure of the specimen after solution treatment

2 實驗結果與分析

2.1真應力-應變曲線分析

圖2為0Cr16Ni5Mo馬氏體不銹鋼在不同變形條件下的真應力-應變曲線。由圖2可知,變形條件對0Cr16Ni5Mo馬氏體不銹鋼的流變應力的影響較大:在同一應變速率時,隨著變形溫度的升高,材料的流變應力下降;在同一變形溫度下, 流變應力隨著應變速率的增大而增加。當變形速率和變形溫度固定的條件下,在壓縮初期,流變應力隨著應變的增大而迅速增加,此時為加工硬化階段;隨后出現兩種情況,一種是應變速率較小(0.01s-1)時,流變應力到達峰值后下降或者逐漸趨于穩定,此時表現為動態再結晶;另一種是應變速率較大(>0.01s-1)時,流變應力隨應變的增加趨勢放緩,達到一定值后,隨著應變的繼續增加,流變應力變化不大,沒有明顯的峰值應力,這種情況單從應力-應變曲線上很難判斷是發生了動態回復還是動態再結晶, 需要借助微觀組織的觀察來進一步確定其狀態。

圖2 0Cr16Ni5Mo馬氏體不銹鋼應力-應變曲線Fig.2 Stress-strain curves of 0Cr16Ni5Mo martensitic stainless =10s-1

2.2流變應力模型的建立

(1)

(2)

式中:A為常數;σ為曲線的峰值應力或者流變應力;α為應力因子;n為與應變速率敏感性相關的指數。

在低應力(ασ<0.8)時,式(2)可簡化為

(3)

在高應力(ασ>1.2)時,式(2)可簡化為

(4)

其中α,β,n1為常數,且滿足

α=β/n1

(5)

由圖2可知,應變量對該材料在較高變形速率時的流變應力σ有較大的影響,但式(2)中并未考慮到應變量的影響,因此建立與應變量相關的本構關系,對于更準確地預測材料的流變應力σ具有重要的意義。本工作假定材料熱變形參數A,α,n和Q等均與應變量有關,根據式(2)雙曲正弦函數模型來建立0Cr16Ni5Mo不銹鋼在不同應變量下的本構關系。依次在0Cr16Ni5Mo不銹鋼的真應力-真應變曲線上選取ε=0.1~0.9(每隔0.1取一數值)時對應的流變應力σ,然后求各個應變量下材料的熱變形參數,最后對各個參數進行擬合得到與應變量相關的多項式函數。

(6)

(7)

將以上求得的數值帶入公式(1),則Zener-Hollomon參數為:

(8)

根據Z=A[sinh(ασ)]n對其兩邊取對數,則有

lnZ=lnA+nln[sinh(ασ)]

(9)

根據公式(9)作lnZ-ln[sinh(ασ)]之間的關系曲線,如圖4所示,可以看出,實驗值基本落在回歸直線的附近,線性相關系數為0.982。通過確定回歸直線的截距即可確定式(9)中的lnA=36.64。因此可知應變量在0.4時材料的熱變形參數為:α=0.007,n=6.342,Q=426.324kJ/mol,lnA=36.64。其他應變量下的熱變形參數均按上述方法計算。

圖3 0Cr16Ni5Mo不銹鋼峰值應力與變形速率 (a)和變形溫度 (b)的關系(ε=0.4)Fig.3 Relationship between peak stress and deformation temperature (a) and strain rate (b) at ε=0.4 for 0Cr16Ni5Mo stainless steel

圖4 0Cr16Ni5Mo不銹鋼峰值應力與Z參數的關系(ε=0.4)Fig.4 Relationship between ln[sinh(ασ)] and lnZ (ε=0.4)

圖5為0Cr16Ni5Mo不銹鋼的熱變形參數α,n,Q和lnA隨應變量的變化關系。由圖(5)可知,其變形參數均隨著應變量的增大而發生變化,這主要是因為隨著變形的不斷深入,其動態再結晶程度不斷加大而使晶粒細化。以熱激活能Q的變化為例,應變量從0.1增大至0.7時,熱激活能降低,這主要是由于在變形初期,大多滑移系和晶粒并未開動,因而此時材料的熱變

形激活能Q相對較高,但是隨著應變量的增大,使更多的滑移系被激活啟動,導致激活能Q開始隨之下降;而當應變量高于0.7后,其熱激活能隨應變量的增加略有上升,這是由于材料在應變量進一步增大后,動態再結晶作用逐漸增大,晶內大量的位錯消失,晶粒也不斷細化,導致變形激活能也開始隨變形量的增大而略有增加。

根據計算值的變化趨勢用Origin軟件擬合材料的熱變形參數α,n,Q和lnA隨應變量變化的函數關系,可得函數關系式(10)~(13),可以看出,各擬合曲線均與實驗數據的計算值吻合較好。

α=0.0111ε4-0.0347ε3+0.0369ε2-0.0163ε+0.00957

(10)

n=26.64ε4-59.33ε3+53.81ε2-25.39ε+11.01

(11)

Q=2563.72ε4-5793.32ε3+5059.49ε2-

2160.19ε+785.86

(12)

lnA=-103.47ε4+211.67ε3-96.88ε2-

31.16ε+54.14

(13)

圖5  0Cr16Ni5Mo不銹鋼材料參數與應變量之間的變化關系 (a)α-ε;(b)n-ε;(c)Q-ε;(d)lnA-εFig.5 Relations between material parameters and strain for 0Cr16Ni5Mo stainless steel (a)α-ε;(b)n-ε;(c)Q-ε;(d)lnA-ε

將式(10),(11),(12),(13)代入式(2)可得新的本構方程式:

(14)

經過優化的本構方程式(14)更好地反映了應變量對本構關系的影響,同時在實際生產中也具有更重要的意義。

2.30Cr16Ni5Mo不銹鋼熱加工圖的建立

(15)

(16)

應變速率敏感性指數m為G與J之間的分配系數:

(17)

(18)

在功率耗散圖中,不是所有的高能量耗散區域都適合熱加工,同時還必須考慮材料的加工失穩區。Prasad等[16-18]根據不可逆力學極值原理推導出材料的流變失穩條件:

(19)

材料在真應變分別為0.6和0.8的熱加工圖如圖6所示,圖6中的數字為能量耗散率η,陰影部分則代表失穩區。在熱加工圖中,能量耗散效率值高于0.3時即可認為是該區域具有較佳加工性能[19]。由圖6可知,變形溫度的增加或應變速率的降低,都會使材料的能量耗散效率增加,在高功率耗散區域更有利于材料的加工。當應變量為0.6時,在低溫變形區域或高變形速率區域,能量耗散功率都不高。當變形溫度在983~1150℃,應變速率為0.01~0.4s-1時,該區域的能量耗散系數可以達到0.3~0.46,尤其變形溫度為1150℃,應變速率為0.01s-1時,其能量耗散效率達到46%,該區域具有良好的加工性能。當應變量為0.8時,變形溫度在980~1150℃,應變速率為0.01~0.2s-1,該區域的能量耗散系數可以達到0.3~0.51,而在整個低溫變形區域(900~980℃)內,其能量耗散效率均較低,因此該區域容易出現組織不均勻等問題。

從圖6還可以看出,應變量對該材料的失穩圖的影響較大。該材料的失穩區主要出現在低變形溫度和高應變速率區域。應變量為0.6時,出現了兩個失穩區,一個是在變形溫度947~1081℃,應變速率0.88~10s-1區域;另一個失穩區比較小,在變形溫度900~910℃,應變速率為0.71~10s-1之間。當應變量增加到0.8時,失穩區的面積增大,此時材料的失穩區主要集中在低溫高應變速率區域,其變形溫度為900~1086℃,應變速率在1s-1以上。實際生產中,在選擇熱加工工藝參數時,應盡量避免在失穩區域進行塑性變形,以免出現對微觀組織不利的各種缺陷,影響加工后材料的性能。綜上所述,可以獲得0Cr16Ni5Mo低碳馬氏體不銹鋼的最佳熱變形工藝參數范圍:變形溫度為980~1150℃,應變速率為0.01~0.2s-1。這個區域的能量耗散效率較高同時又避開了失穩區,有利于材料的加工,并可以獲得均勻無缺陷的組織。

圖6 0Cr16Ni5Mo不銹鋼的熱加工圖 (a)ε=0.6;(b)ε=0.8Fig.6 Processing maps of 0Cr16Ni5Mo stainless steel at different strains  (a)ε=0.6;(b)ε=0.8

2.40Cr16Ni5Mo不銹鋼組織的演變

以真應變為0.8的熱加工圖為例對不同區域的熱變形試樣的組織進行觀察,如圖7所示。熱變形組織可以分為以下4類:1)形變組織(圖7(b)),此時對應的區域變形溫度為900℃,應變速率為10s-1, 材料處于失穩區,這時在變形過程中只發生動態回復,組織為變形的大晶粒,未觀察到動態再結晶跡象;2)部分再結晶組織(圖7(a),(c)),這時的熱變形組織發生了不同程度的動態再結晶,但是動態再結晶進行得不夠充分,組織為變形的長條狀晶粒和一些細小的動態再結晶晶粒組成的混晶,在此區域進行熱加工會導致材料的性能不均勻;3)完全再結晶組織(圖7(d),(e)),對應的區域變形溫度為1100℃,應變速率分別為0.1,1s-1,此時在熱變形過程中發生了完全動態再結晶,獲得大小均勻的等軸狀再結晶晶粒;4)粗大晶粒組織(圖7(f)),對應的區域變形溫度為1150℃,應變速率0.01s-1,其η值達到51%,由于η的提高,此時發生動態再結晶的程度增大。然而,由于變形溫度較高并且應變速率低,再結晶晶粒出現長大的情況。

圖7 不同變形條件下0Cr16Ni5Mo不銹鋼的典型熱變形組織 (a)T=900℃,℃, (c) T=950℃,℃,℃,℃,Fig.7 Typical microstructures of 0Cr16Ni5Mo stainless steel after hot deformation at different conditions (a)T=900℃, ℃,℃,℃,(e)T=1100℃,℃,

3 結論

(1) 0Cr16Ni5Mo馬氏體不銹鋼在高溫壓縮變形時,流變應力隨變形溫度的升高而降低,隨應變速率的升高而增加。當應變速率較低時,有比較明顯的峰值應力,而隨著應變速率升高,流變曲線變得平緩,峰值應力不顯著。

(2) 基于雙曲正弦模型,建立了與應變量相關的本構關系模型,計算結果表明材料熱變形參數與應變量之間可采用四次函數關系式表示,并且具有很好的相關性。

(3) 通過熱加工圖分析可知,0Cr16Ni5Mo低碳馬氏體不銹鋼的最佳熱變形工藝參數范圍:變形溫度為980~1150℃,應變速率為0.01~0.2s-1。

[1]QIN B,WANG Z Y, SUN Q S. Effect of tempering temperature on properties of 00Cr16Ni5Mo stainless steel[J]. Materials Characterization,2008,59(8):1096-1100.

[2]MA X P, WANG L J, QIN B, et al. Effect of N on microstructure and mechanical properties of 16Cr5Ni1Mo martensitic stainless steel[J]. Materials and Design, 2012,34: 74-81.

[3]DAWOOD M A, MAHALLAWI IS E, AZIM ME A, et al. Thermal aging of 16Cr-5Ni-1Mo stainless steel: part 1-microstructural analysis[J]. Materials Science and Technology, 2004, 20(3): 363-369.

[4]WANG Z H, FU W T, WANG B Z, et al. Study on hot deformation characteristics of 12%Cr ultra-super-critical rotor steel using processing maps and Zener-Hollomon parameter[J]. Materials Characterization, 2010, 61(1):25-30.

[5]俞秋景, 張偉紅, 于連旭, 等. 鑄態Inconel 625合金熱加工圖的建立及熱變形機制分析[J].材料工程,2014,(1):30-34.

YU Qiu-jing, ZHANG Wei-hong, YU Lian-xu, et al. Development of thermal processing map and analysis of hot deformation mechanism of cast alloy Inconel 625[J]. Journal of Materials Engineering, 2014, (1): 30-34.

[6]熊毅,熊良銀,張凌峰,等. GH4199合金的熱變形行為與微觀組織演變[J]. 中國有色金屬學報,2010,20(4):655-661.

XIONG Yi, XIONG Liang-yin, ZHANG Ling-feng, et al. Hot deformation behavior and microstructure evolution of superalloy GH4199[J]. The Chinese Journal of Nonferrous Metals, 2010, 20(4):655-661.

[7]張威,閆東娜,鄒德寧,等.超低碳13Cr-5Ni-2Mo馬氏體不銹鋼熱變形行為及本構關系[J].鋼鐵,2012,47(5):69-74.

ZHANG Wei, YAN Dong-na, ZOU De-ning, et al. Hot deformation behavior and constitutive relationship for super-low 13Cr-5Ni-2Mo martensitic stainless steel[J].Iron and Steel, 2012,47(5):69-74.

[8]吳凱,劉國權,胡本芙,等. 新型鎳基粉末高溫合金的高溫變形行為[J].航空材料學報,2010,30(4):1-7.

WU Kai, LIU Guo-quan, HU Ben-fu, et al. Research on high temperature deformation behavior of new type nickel-based P/M superalloy[J]. Journal of Aeronautical Materials, 2010,30(4):1-7.

[9]孫朝陽,劉金榕,李瑞,等.Incoloy 800H高溫變形流動應力預測模型[J].金屬學報,2011,47(2):191-196.

SUN Chao-yang, LIU Jin-rong, LI Rui,et al. Constitutive modeling for elevated temperature flow behavior of Incoloy 800H superalloy[J]. Acta Metallurgica Sinica,2011,47(2):191-196.

[10]王智祥, 劉雪峰,謝建新.AZ91鎂合金高溫變形本構關系[J].金屬學報,2008,44(11):1378-1383.

WANG Zhi-xiang, LIU Xue-feng, XIE Jian-xin. Constitutive relationship of hot deformation of AZ91 magnesium alloy[J].Acta Metallurgica Sinica,2008,44(11):1378-1383.

[11]李波, 潘清林, 張志野, 等. 含鈧Al-Zn-Mg合金的熱變形行為和顯微組織[J]. 材料工程, 2013, (11):6-11.

LI Bo, PAN Qing-lin, ZHANG Zhi-ye, et al. Hot deformation behavior and microstructure of Al-Zn-Mg alloy containing Sc during hot compression at elevated temperature[J]. Journal of Materials Engineering, 2013, (11): 6-11.

[12]劉娜, 李周, 袁華, 等. 粉末冶金TiAl合金的熱變形行為研究[J]. 航空材料學報, 2013,33(5):1-5.

LIU Na, LI Zhou, YUAN Hua, et al. Hot deformation behavior of powder metallurgical TiAl alloy[J]. Journal of Aeronautical Materials, 2013, 33(5): 1-5.

[13]黃由林,王建波,凌學士,等.熱加工圖理論的研究進展[J].材料導報,2008,22(增刊3): 173-176.

HUANG You-lin, WANG Jian-bo, LING Xue-shi, et al. Research development of hot processing map theory[J]. Materials Review, 2008,22(Suppl 3): 173-176.

[14]康婭雪,蔡大勇,張春玲,等. 微碳鋼的熱變形方程及熱加工圖[J].材料熱處理學報,2012,33(6):74-79.

KANG Ya-xue,CAI Da-yong,ZHANG Chun-ling, et al. Hot deformation equation and processing map of a micro-carbon steel[J]. Transactions of Materials and Heat Treatment,2012,33(6):74-79.

[15]馬龍騰,王立民,胡勁,等. AISI403馬氏體不銹鋼的熱變形特性研究[J]. 材料工程,2013,(5):38-43.

MA Long-teng, WANG Li-min, HU Jin, et al. Hot deformation features of AISI403 martensitic stainless steel[J]. Journal of Materials Engineering, 2013, (5):38-43.

[16]SRINIVASAN N, PRASAD Y V R K. Hot working characteristics of nimonic75, 80A and 90 superalloys: a comparison using processing maps[J]. Journal of Materials Processing Technology,1995, 51(1-4):171-192.

[17]PRASAD Y V R K. Processing maps:a status report[J]. Journal of Materials Engineering and Performance,2003,12(6):638-645.

[18]PRASAD Y V R K, GEGEL H L, DORAIVELU S M. Modeling of dynamic material behavior in hot deformation: forging of Ti-6242[J]. Metallurgical and Materials Transactions A,1984,15(10):1883-1892.

[19]于洋,熊柏青,惠松驍,等.鑄態Ti-6Al-4V-0.1B合金的熱變形行為及加工圖[J].材料熱處理學報,2012,33(6):142-146.

YU Yang, XIONG Bai-qing, HUI Song-xiao,et al. Hot deformation behavior and processing map of as-cast Ti-6Al-4V-0.1B alloy[J].Transactions of Materials and Heat Treatment,2012,33(6):142-146.

Hot Deformation Behavior and Processing Map of 0Cr16Ni5Mo Low Carbon Martensitic Stainless Steel

YUAN Wu-hua1,GONG Xue-hui1,2,SUN Yong-qing2,LIANG Jian-xiong2

(1 College of Materials Science and Engineering,Hunan University,Changsha 410082,China;2 Institute for Special Steel,Central Iron and Steel Research Institute,Beijing 100081,China)

The hot deformation behavior of 0Cr16Ni5Mo low carbon martensitic stainless steel was studied by the isothermal compression of cylindrical specimens at 900-1150℃ with the strain rate of 0.01-10s-1on a Gleeble-3800 simulated machine. The relations of the thermomechanical parameters with strain were obtained using the hyperbolic-sine mathematics model and the hot deformation constitutive relationship was established. Processing map was also established based on the dynamic materials model. The microstructure evolution at different conditions was analyzed. The results show that the flow stress decreases with the increase of deformation temperature and increases with the increase of strain rate, the deformation condition has a great influence on the material microstructure. The relationship between deformation parameters and strain have a good relativity, which can be expressed using four polynomial fitting. An optimum processing parameters of hot deformation for this steel can also be obtained by the maps, in which the hot temperature is 980-1150℃ and the strain rate is 0.01-0.2s-1.

0Cr16Ni5Mo martensitic stainless steel;flow stress;constitutive model;dynamic recrystallization;processing map

袁武華(1973-),男,教授,博士,主要從事金屬與金屬基復合材料相關方面研究,聯系地址:湖南省長沙市岳麓區湖南大學材料科學與工程學院(410082),E-mail:yuan46302@163.com

10.11868/j.issn.1001-4381.2016.05.002

TG142.71

A

1001-4381(2016)05-0008-07

2014-09-16;

2015-11-21

猜你喜歡
不銹鋼變形
談詩的變形
中華詩詞(2020年1期)2020-09-21 09:24:52
80t不銹鋼GOR底吹轉爐工藝設備改造
山東冶金(2019年1期)2019-03-30 01:35:32
“我”的變形計
你知道不銹鋼中“304”的含義嗎
不銹鋼微鉆削的切屑形成與仿真分析
變形巧算
FV520(B)不銹鋼焊接接頭的斷裂性能
例談拼圖與整式變形
會變形的餅
關于不銹鋼厚壁管的焊接
主站蜘蛛池模板: a毛片免费在线观看| 国产亚洲精品无码专| 毛片最新网址| 成人日韩欧美| 国产不卡在线看| 亚洲成人动漫在线观看| 欧美一区精品| 国产美女无遮挡免费视频| 一级做a爰片久久毛片毛片| 一级一毛片a级毛片| 亚洲一区二区三区麻豆| 4虎影视国产在线观看精品| 在线视频亚洲色图| 日本精品视频一区二区| 2021国产v亚洲v天堂无码| 天天爽免费视频| 在线精品亚洲一区二区古装| 精品国产www| 国产理论最新国产精品视频| 四虎亚洲国产成人久久精品| 国产女人18毛片水真多1| 欧美成人精品在线| 特级毛片8级毛片免费观看| 国产精品大白天新婚身材| 免费观看精品视频999| 国产精品林美惠子在线观看| 婷婷六月激情综合一区| 18禁黄无遮挡免费动漫网站| 91福利在线看| 综合天天色| 色噜噜久久| 婷婷久久综合九色综合88| 四虎AV麻豆| 亚洲区视频在线观看| 女人18一级毛片免费观看| 日韩欧美成人高清在线观看| 色婷婷色丁香| 欧美区在线播放| 亚洲首页在线观看| av大片在线无码免费| 婷婷色中文| 99热这里只有精品在线播放| 日本一区二区三区精品视频| 尤物国产在线| 日韩在线播放欧美字幕| 亚洲无码免费黄色网址| 国产日韩久久久久无码精品| 欧美成人一级| 中文成人在线视频| 91色在线观看| 欧美一级色视频| 久久动漫精品| 99无码中文字幕视频| 国产成年女人特黄特色毛片免| 国产又黄又硬又粗| 久久国产高清视频| 午夜电影在线观看国产1区| 九九久久精品免费观看| 国产手机在线小视频免费观看| 91麻豆国产视频| 69综合网| 亚洲午夜福利在线| 欧美www在线观看| a级毛片视频免费观看| 中文字幕在线免费看| 亚洲熟妇AV日韩熟妇在线| 青草91视频免费观看| 91成人在线观看视频| 亚洲天堂首页| 成人国产精品网站在线看| 國產尤物AV尤物在線觀看| 免费又黄又爽又猛大片午夜| 精品福利视频网| 免费国产高清视频| 国产第八页| 无码免费视频| 免费国产高清视频| 2020国产免费久久精品99| 又黄又湿又爽的视频| 久久亚洲精少妇毛片午夜无码| a级毛片一区二区免费视频| 亚洲中文字幕av无码区|