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

冷壓縮法消減7A85鋁合金鍛件淬火殘余應(yīng)力研究

2015-10-10 07:57:48代偉易幼平崔金棟
關(guān)鍵詞:方向有限元變形

代偉,易幼平,崔金棟

?

冷壓縮法消減7A85鋁合金鍛件淬火殘余應(yīng)力研究

代偉,易幼平,崔金棟

(中南大學(xué)機(jī)電工程學(xué)院,湖南長(zhǎng)沙,410083)

為研究冷壓縮法對(duì)7A85鋁合金淬火殘余應(yīng)力消減效果,通過(guò)實(shí)驗(yàn)獲得7A85鋁合金的比熱容、熱傳導(dǎo)系數(shù)、熱膨脹系數(shù)等熱物性參數(shù),利用ABAQUS有限元軟件對(duì)7A85鋁合金試件進(jìn)行淬火及冷壓縮變形工藝仿真,研究淬火及冷變形工藝參數(shù)對(duì)試件殘余應(yīng)力的影響,對(duì)7A85鋁合金試件淬火與冷壓縮變形工藝進(jìn)行實(shí)驗(yàn)研究,采用X線衍射儀分別測(cè)試試件淬火及冷變形后的殘余應(yīng)力。研究結(jié)果表明:當(dāng)冷變形量為1%時(shí),冷壓縮法可降低試件約70%的淬火殘余應(yīng)力,隨著冷變形量的增加,殘余應(yīng)力變化較小,表現(xiàn)出與拉伸法消減殘余應(yīng)力不同的變化特征。實(shí)驗(yàn)與仿真分析獲得的殘余應(yīng)力絕對(duì)誤差在40 MPa以?xún)?nèi)。

7A85鋁合金;淬火;殘余應(yīng)力;冷壓縮

7A85鋁合金是新一代高強(qiáng)韌、高淬透性航空結(jié)構(gòu)材料,被廣泛應(yīng)用于大型飛機(jī)主承力構(gòu)件中[1]。固溶淬火是超高強(qiáng)鋁合金構(gòu)件熱處理強(qiáng)化過(guò)程中重要工序之一,通過(guò)淬火使第二相溶質(zhì)完全固溶于基體后,經(jīng)人工時(shí)效使第二相元素在組織中均勻分布,從而提高其強(qiáng)度、韌度、耐腐蝕、抗疲勞等綜合力學(xué)性能[2]。但固溶淬火過(guò)程形成的殘余應(yīng)力易引起構(gòu)件在機(jī)加工過(guò)程產(chǎn)生較大變形甚至開(kāi)裂,目前已成為我國(guó)航空構(gòu)件制造面臨的重大難題。許多學(xué)者圍繞構(gòu)件的殘余應(yīng)力問(wèn)題進(jìn)行了理論與試驗(yàn)研究,主要針對(duì)預(yù)拉伸厚板淬火工藝及殘余應(yīng)力形成機(jī)制[3?4]、拉伸法消減鋁合金厚板淬火殘余應(yīng)力[5?7]等方面開(kāi)展工作。這些對(duì)殘余應(yīng)力控制及消減方法主要適用于大厚板。雖然預(yù)拉伸板在航空構(gòu)件制造方面具有工藝穩(wěn)定、殘余應(yīng)力小、適應(yīng)性強(qiáng)等優(yōu)勢(shì),但也存在構(gòu)件機(jī)加工過(guò)程流線被切斷、力學(xué)性能各向異性、材料浪費(fèi)嚴(yán)重等不足。同時(shí),飛機(jī)關(guān)鍵承力構(gòu)件通常設(shè)計(jì)為模鍛件,以適應(yīng)飛機(jī)構(gòu)件輕量化、整體化的發(fā)展趨勢(shì),單一的拉伸法消減殘余應(yīng)力不再適用于結(jié)構(gòu)復(fù)雜的模鍛件。國(guó)外航空制造企業(yè)普遍采用冷壓縮變形(Tx52)對(duì)模鍛件進(jìn)行淬火殘余應(yīng)力消減[8?10],國(guó)內(nèi)尚未在此領(lǐng)域開(kāi)展系統(tǒng)深入的研究。本文作者針對(duì)航空模鍛件Tx52消減殘余應(yīng)力工藝,結(jié)合國(guó)內(nèi)最新的7A85超高強(qiáng)鋁合金材料,開(kāi)展材料基礎(chǔ)試驗(yàn),以7A85鋁合金試件為例,通過(guò)有限元仿真方法研究固溶淬火工藝對(duì)材料殘余應(yīng)力的影響規(guī)律,分析不同冷壓縮變形量對(duì)試件淬火殘余應(yīng)力消減效果,在有限元分析的基礎(chǔ)上進(jìn)行7A85鋁合金試件淬火及冷壓縮變形工藝試驗(yàn)與殘余應(yīng)力測(cè)試,為復(fù)雜模鍛件的T7x52工藝研究提供了基礎(chǔ)。

1 有限元模型及材料參數(shù)

以7A85鋁合金長(zhǎng)×寬×高為100 mm×60 mm× 40 mm試件為例研究7A85-Tx52工藝,借助有限元仿真方法分析試件在固溶淬火、冷壓縮變形過(guò)程中殘余應(yīng)力演變規(guī)律。依據(jù)目前國(guó)內(nèi)企業(yè)已有的熱處理工藝制度,試件加熱至470 ℃,保溫4 h;試件入水轉(zhuǎn)移時(shí)間10 s,淬火水溫25℃,試件厚度40 mm。由于7A85鋁合金試件在冷壓縮變形后有一定回彈,根據(jù)仿真及實(shí)驗(yàn)測(cè)試結(jié)果設(shè)計(jì)了補(bǔ)償量,冷壓縮方案如表1所示。

表1 冷壓縮方案

材料物性參數(shù)是進(jìn)行淬火與冷變形工藝有限元仿真的基礎(chǔ),對(duì)7A85鋁合金材料主要參數(shù)進(jìn)行了測(cè)試。

1) 采用排水法測(cè)試密度,測(cè)得密度為2 820 kg/m3。

2) 采用差熱法(DSC)測(cè)量比熱容。實(shí)驗(yàn)結(jié)果如圖1所示,7A85鋁合金比熱容在410 ℃較高,達(dá)到1 253.2 J/(kg·℃),存在明顯的峰值。

圖1 7A85鋁合金比熱容

3) 采用JR?3激光導(dǎo)熱儀測(cè)量導(dǎo)熱系數(shù)。先通過(guò)熱脈沖法測(cè)量材料熱擴(kuò)散率,再計(jì)算出導(dǎo)熱系 數(shù)[11?12]。在獲得材料比熱容及熱擴(kuò)散率后,導(dǎo)熱系數(shù)可由式(1)求出[13]:

7A85鋁合金熱膨脹系數(shù)采用德國(guó)NETZSCH DIL402C熱膨脹儀在Ar氣保護(hù)下測(cè)量。7A85鋁合金導(dǎo)熱系數(shù)及熱膨脹系數(shù)如表2所示。

表2 7A85鋁合金導(dǎo)熱系數(shù)及熱膨脹系數(shù)

4) 7A85材料力學(xué)性能實(shí)驗(yàn)數(shù)據(jù)參見(jiàn)文獻(xiàn)[14]。在應(yīng)力?應(yīng)變曲線中,將應(yīng)變分解為彈性應(yīng)變與塑性應(yīng)變,并輸入ABAQUS材料特性。對(duì)于7A85鋁合金材料應(yīng)力?應(yīng)變的非線性關(guān)系采用Ramberg?Osgood方程描述,對(duì)于塑性段,總應(yīng)變由彈性應(yīng)變和塑性應(yīng)變2部分構(gòu)成[15?16],

其中:為應(yīng)變;為應(yīng)力;為彈性模量;和為與材料相關(guān)的常數(shù)。

在分析過(guò)程中,材料的各種熱物性參數(shù)隨溫度的變化而變化,可用于計(jì)算7A85鋁合金材料在淬火過(guò)程中殘余應(yīng)力的形成過(guò)程。

2 淬火及冷壓縮變形殘余應(yīng)力仿真分析

利用ABAQUS有限元軟件對(duì)7A85鋁合金試件淬火及冷壓縮變形的殘余應(yīng)力進(jìn)行仿真分析。7A85鋁合金試件淬火結(jié)束后方向應(yīng)力分布如圖2所示,試件淬火后的殘余應(yīng)力呈典型的內(nèi)拉外壓特征。由于試件具有對(duì)稱(chēng)性,取試件芯部點(diǎn)1與表面點(diǎn)2進(jìn)行分析,位置如圖2所示。為方便實(shí)驗(yàn)結(jié)果與仿真值比對(duì)分析,取試件仿真結(jié)果側(cè)面及截面位置的中間點(diǎn)測(cè)量其應(yīng)力,分別表征試件,和方向淬火后的殘余應(yīng)力,其方向及位置如圖2所示。

圖2 7A85鋁合金試件淬火Y方向應(yīng)力場(chǎng)

采取下模固定、上模下壓的方式對(duì)試件進(jìn)行冷壓縮變形,下壓速度為0.1 mm/s,冷壓縮后,試件芯部點(diǎn)1和表面點(diǎn)2的殘余應(yīng)力隨壓縮率變化曲線如圖3所示,其中試件未變形時(shí)的應(yīng)力即為淬火殘余應(yīng)力。由圖3可知:試件淬火后,尺寸方向最大殘余應(yīng)力比其他2個(gè)方向的大,且表面點(diǎn)2的法向殘余應(yīng)力約為0 MPa。模壓1%對(duì)試件殘余應(yīng)力消減效果較好,芯部點(diǎn)1的方向殘余應(yīng)力由265 MPa降至?33.8 MPa,表面點(diǎn)2的殘余應(yīng)力由?187.9 MPa降至51.7 MPa,冷變形可將最大值方向的殘余應(yīng)力消減70%左右。在變形量1%~5%范圍內(nèi),冷壓縮變形時(shí)殘余應(yīng)力隨變形量的增加變化較小,表現(xiàn)出與拉伸法消減殘余應(yīng)力不同的變化規(guī)律。

試件1與2點(diǎn)在壓縮率為1%時(shí)殘余應(yīng)力變化曲線如圖4所示,由于試件淬火后有一定變形,上下模初始位置距試件各為2 mm,上模下壓4.0 mm時(shí)合模。由圖4可知:試件在上模下壓4.0~4.6 mm時(shí)應(yīng)力急劇增大,且承載的方向應(yīng)力增量最大,在上模回退后各方向應(yīng)力回復(fù)并形成殘余應(yīng)力。試件冷變形后,芯部點(diǎn)1的方向殘余應(yīng)力由147.3 MPa減至?22.4 MPa,方向的殘余應(yīng)力由265.1 MPa減至?34.8 MPa,方向殘余應(yīng)力雖然在加載過(guò)程中有較大的增加,但卸載回復(fù)后,應(yīng)力水平無(wú)較大變化;表面點(diǎn)2的方向應(yīng)力為垂直表面的應(yīng)力,其絕對(duì)值在模壓過(guò)程中均無(wú)明顯變化,方向殘余應(yīng)力水平由?175.5 MPa減至51.7 MPa,方向的殘余應(yīng)力由?133.4減至?6.6 MPa。通過(guò)冷變形可消減本文試件70%左右的淬火殘余應(yīng)力,并使試件殘余應(yīng)力分布更均勻。

(a) P1點(diǎn);(b) P2點(diǎn)

淬火及冷變形后,按圖2所示試件側(cè)面及截面中心點(diǎn)殘余應(yīng)力仿真結(jié)果如表3所示,試件表面中間點(diǎn)淬火殘余應(yīng)力在130~180 MPa,且尺寸較大的長(zhǎng)度方向殘余應(yīng)力絕對(duì)值較大。試件冷變形1%后,各個(gè)方向的殘余應(yīng)力絕對(duì)值均有顯著下降,試件殘余應(yīng)力隨著變形量的增加絕對(duì)值變化比較小,這是由于試件形狀規(guī)則、尺寸小,且7A85鋁合金材料延展性較好。

表3 7A85鋁合金試件淬火及模壓仿真殘余應(yīng)力

3 淬火及冷壓縮變形工藝實(shí)驗(yàn)

7A85鋁合金試件經(jīng)線切割與表面打磨后固溶處理,固溶溫度為470 ℃,保溫4 h后進(jìn)行水淬,采用X線衍射儀測(cè)量淬火后試件的殘余應(yīng)力。測(cè)量時(shí),輻射線CoK,X線管流為4 mA,管壓為20 kV,曝光時(shí)間為2 s,每個(gè)位置曝光次數(shù)為10次,儀器對(duì)試件掃描后計(jì)算機(jī)自動(dòng)采集并處理數(shù)據(jù),給出殘余應(yīng)力。將淬火后的試件進(jìn)行冷壓縮變形,實(shí)驗(yàn)設(shè)備為4 000 t數(shù)控液壓機(jī),上下模具采用平砧。對(duì)各個(gè)試件按標(biāo)記分別施加1%,2%,3%,4%和5%的壓縮率,模壓速度為0.1 mm/s,壓縮量由液壓控制系統(tǒng)精確控制。

壓縮后采用X線衍射儀按圖2所示方向測(cè)量試件側(cè)面和截面中心點(diǎn)殘余應(yīng)力,分別表征試件,和方向模壓后的殘余應(yīng)力。淬火及冷變形后殘余應(yīng)力如表4所示,其中壓縮率這0的應(yīng)力即試件淬火殘余應(yīng)力。實(shí)驗(yàn)結(jié)果表明:對(duì)試件施加冷變形1%可使各個(gè)方向的殘余應(yīng)力絕對(duì)值顯著下降,冷壓縮法削減7A85鋁合金鍛件淬火殘余應(yīng)力具有較理想的效果,且試件殘余應(yīng)力隨著冷壓縮變形量的增加絕對(duì)值變化較小。

表4 7A85鋁合金試件淬火及模壓實(shí)驗(yàn)殘余應(yīng)力

7A85鋁合金試件淬火及冷壓縮變形工藝實(shí)驗(yàn)測(cè)得的殘余應(yīng)力與仿真值對(duì)比如圖5所示。可見(jiàn):試件3個(gè)方向上殘余應(yīng)力趨勢(shì)相同,特別當(dāng)淬火及冷變形量較小時(shí),應(yīng)力變化趨勢(shì)一致,絕對(duì)誤差在20 MPa以?xún)?nèi);當(dāng)變形量較大時(shí),實(shí)驗(yàn)值出現(xiàn)波動(dòng),與仿真曲線出現(xiàn)交叉,絕對(duì)誤差在40 MPa以?xún)?nèi)。實(shí)驗(yàn)結(jié)果表明:本文采用的冷壓縮法削減7A85鋁合金鍛件淬火殘余應(yīng)力有限元模型具有較高的準(zhǔn)確性。

(a) X方向;(b) Y方向;(c) Z方向

4 結(jié)論

1) 建立了7A85鋁合金試件淬火與冷壓縮變形有限元仿真模型。

2) 研究了7A85鋁合金試件淬火及冷變形工藝參數(shù)對(duì)試件殘余應(yīng)力的影響,冷壓縮變形量為1%時(shí),可降低試件約70%的淬火殘余應(yīng)力,冷壓縮變形量對(duì)殘余應(yīng)力影響較小,不同于拉伸法消減殘余應(yīng)力的變化特征。

3) 7A85鋁合金試件仿真結(jié)果與試驗(yàn)獲得的殘余應(yīng)力變化規(guī)律一致,絕對(duì)誤差在40 MPa以?xún)?nèi)。

[1] 賈逢博, 易幼平, 黃施全,等. 7A85鋁合金熱壓縮流變行為與本構(gòu)方程研究[J]. 材料熱處理技術(shù), 2010, 39(16): 19?20. JIA Fengbo, YI Youping, HUANG Shiquan, et al. Study on flow behavior and constitutive equation of 7A85 aluminum alloy during hot compression[J]. Material & Heat Treatment, 2010, 39(16): 19?20.

[2] 廖凱. 鋁合金厚板淬火?預(yù)拉伸內(nèi)應(yīng)力形成機(jī)理及其測(cè)試方法研究[D]. 長(zhǎng)沙: 中南大學(xué)機(jī)電工程學(xué)院, 2010: 1, 21?26. LIAO Kai. Researchon mechanism and measurement method of quenching-prestretching stress in aluminum alloy thick plate[D]. Changsha: Central South University. School of Mechanical and Electrical Engineering, 2010: 1, 21?26.

[3] 龔海, 吳運(yùn)新, 廖凱. 不同淬火工藝對(duì)7075鋁合金厚板殘余應(yīng)力的影響[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2010, 41(4): 1354?1359. GONG Hai, WU Yunxin, LIAO Kai.Influence of different quenching techniques on residual stress of 7075 aluminum alloy thick-plate[J].Journal of Central South University (Science and Technology) ,2010, 41(4): 1354?1359.

[4] 胡少虬, 張輝, 楊立斌,等. 7075鋁合金厚板淬火溫度場(chǎng)及熱應(yīng)力場(chǎng)的數(shù)值模擬[J]. 湘潭大學(xué)自然科學(xué)學(xué)報(bào), 2004, 26(2): 66?71. Hu Shaoqiu, ZHANG Hui, YANG Libin, et al.Numerical analysis of temperature and thermal stresses fields in quenching 7075 aluminum alloy thick plates[J].Natural Science Journal of Xiangtan University, 2004, 26(2): 66?71.

[5] 朱偉, 彭大暑, 張輝,等. 7075鋁合金厚板淬火殘余應(yīng)力消除工藝的研究[J]. 加工工藝, 2002, 25(2): 12?14. ZHU Wei, PEN Dashu, ZHANG Hui, et al. Study on process of elimination residual stress in quenching 7075 aluminium alloy plates[J].Processing Technic, 2002, 25(2): 12?14.

[6] 張園園, 吳運(yùn)新, 李麗敏,等. 7075鋁合金預(yù)拉伸板淬火后殘余應(yīng)力的有限元模擬[J]. 材料熱處理技術(shù), 2008, 37(14): 88?91. ZHANG Yuanyuan, WU Yunxin, LI Limin, et al.Finite element simulation of residual stress in prestretching thick-plates of 7075 aluminum alloy after quenching[J].Material & Heat Treatment, 2008, 37(14): 88?91.

[7] 王秋成, 柯映林. 深冷處理消除7050鋁合金淬火殘余應(yīng)力的研究[J]. 浙江大學(xué)學(xué)報(bào)(工學(xué)版), 2003, 37(6): 748?751. WANG Qiucheng, KE Yinglin.Relief of residual stresses in 7050 aluminum alloy by deep cryogenic treatment[J].Journal of Zhejiang University (Engineering Science), 2003, 37(6): 748?751.

[8] 林高用. 高性能7X75系鋁合金厚板加工技術(shù)相關(guān)基礎(chǔ)[D]. 長(zhǎng)沙: 中南大學(xué)材料科學(xué)與工程學(xué)院, 2006: 92?93. LIN Gaoyong.Fundamental research related to the fabrication technology for high quality thick plates of 7X75 series aluminum alloys[D]. Changsha:Central South University. School of Materials Science and Engineering, 2006: 92?93.

[9] 姜海峰, 陸政, 楊洪濤,等. 超高強(qiáng)鋁合金等溫模鍛工藝[J]. 新工藝、新技術(shù)、新設(shè)備, 2002(1): 56?58. JIANG Haifeng, LU Zheng, YANG Hongtao, et al.Ultra high strength aluminum alloy isothermal stamping process[J].New Process, New Technology, New Equipment, 2002(1): 56?58.

[10] 趙懷瑞, 王生武, 溫愛(ài)玲. 模壓強(qiáng)化工藝及其數(shù)值仿真[J]. 計(jì)算機(jī)仿真, 2007, 24(10): 240?244. ZHAO Huairui, WANG Shengwu, WEN Ailing.Die formed strengthen process and its F.E simulation[J].Computer Simulation, 2007, 24(10): 240?244.

[11] Kabayabaya T, YU Fan, ZHANG Xinxin. Thermal diffusivity measurement of glass at high temperature by using flash method[J]. Thermal Science, 2003,13(1): 91?96.

[12] 陳躍飛, 張金濤, 于帆,等. 激光閃光法測(cè)量材料熱擴(kuò)散率國(guó)際比對(duì)的實(shí)驗(yàn)研究與分析[J]. 測(cè)量與設(shè)備, 2008(12): 29?32. CHEN Yuefei, ZHANG Jintoa, YU Fan, et al.Experimental study and analysis of international comparison for the material thermal diffusion rate measuring with laser flash method[J].Measurement and Equipment, 2008(12): 29?32.

[13] 姚山, 陳寶慶, 戴利欣,等. 基于激光快速加熱的金屬材料熱擴(kuò)散率及導(dǎo)熱系數(shù)快速測(cè)量方法的研究[J]. 熱科學(xué)與技術(shù), 2005, 4(1): 87?89. YAO Shan, CHEN Baoqing, DAI Lixin, et al. Study on rapid measure of metallic material heat conductivity based on laser rapid heating[J]. Journal of Thermal Science and Technology, 2005, 4(1): 87?89.

[14] 王少輝. 7A85鋁合金整體結(jié)構(gòu)件淬火殘余應(yīng)力分析及其削減工藝研究[D]. 長(zhǎng)沙: 中南大學(xué)機(jī)電工程學(xué)院, 2011: 14?18. WANG Shaohui.Analysis of quenching residual stress and its reduction process of 7A85 aluminum alloy component[D]. Changsha:Central South University. School of Mechanical and Electrical Engineering, 2011: 14?18.

[15] 謝奇, 呂運(yùn)冰, 徐兆雄. Ramberg?Osgood形冪硬化材料的硬化系數(shù)及硬化指數(shù)的確定[J]. 武漢交通科技大學(xué)學(xué)報(bào), 1996, 20(3): 319?322. XIE Qi, Lü Yunbing, XU Zhaoxiong.The determination of the hardening coefficient and exponent for Ramberg?Osgoods power hardening material[J].Journal of Wuhan Transportation University, 1996, 20(3): 319?322.

[16] 胡勤, 戚承志. Ramberg?Osgood土動(dòng)力非線性模型在ABAQUS 軟件上的開(kāi)發(fā)及應(yīng)用[J]. 巖土力學(xué), 2012, 33(4): 1268?1274. HU Qin, QI Chengzhi.Development and application of Ramberg?Osgood soil dynamicnonlinear constitutive model on ABAQUS code[J].Rock and Soil Mechanics, 2012, 33(4): 1268?1274.

Reduction of quenching residual stress for 7A85 aluminum alloy forging by cold compression

DAI Wei, YI Youping, CUI jindong

(School of Mechanical and Electrical Engineering, Central South University, Changsha 410083, China)

In order to study the quenching residual stress reducing effect of 7A85 aluminum alloy by cold compression method, the experiments on 7A85 aluminum alloy were carried out to get thermal properties parameters such as the specific heat, thermal conductivity and thermal expansion coefficient. The quenching and cold compression processes for 7A85 aluminum alloy specimen were simulated by ABAQUS finite element software. The influences of process parameters on the quenching and cold compressed residual stress were studied. The experiments of quenching and cold compressed processes for 7A85 aluminum alloy specimen were performed and the residual stresses of specimens were tested by X-ray diffraction. The results show that the quenching residual stress decreases by 70% on the condition of 1% cold compression. The residual stress changes a little with the increase of cold compression, which displays different characters from the method of residual stress reduction by tension. The absolute error between the tested results of residual stress and the simulated results is less than 40 MPa.

7A85 aluminum alloy; quenching; residual stress; cold compression

10.11817/j.issn.1672-7207.2015.05.006

V414.3+2

A

1672?7207(2015)05?1609?06

2015?02?02;

2015?05?15

國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃(973計(jì)劃)項(xiàng)目(2012CB619504) (Project(2012CB619504) supported by the National Basic Research Program of China (973 Program))

易幼平,博士,教授,從事航空鍛件制造工藝與組織性能控制技術(shù)研究;E-mail: yyp@mail.csu.edu.cn

(編輯 趙俊)

猜你喜歡
方向有限元變形
2022年組稿方向
2021年組稿方向
2021年組稿方向
談詩(shī)的變形
“我”的變形計(jì)
例談拼圖與整式變形
會(huì)變形的餅
磨削淬硬殘余應(yīng)力的有限元分析
位置與方向
基于SolidWorks的吸嘴支撐臂有限元分析
主站蜘蛛池模板: 久久亚洲中文字幕精品一区| 91探花在线观看国产最新| 亚洲av综合网| 一级毛片免费观看久| 亚洲欧美综合在线观看| 色婷婷天天综合在线| 美女免费黄网站| 国产欧美亚洲精品第3页在线| 色国产视频| 国产福利微拍精品一区二区| 成人午夜免费观看| 婷婷久久综合九色综合88| 人妻无码中文字幕一区二区三区| 午夜无码一区二区三区| 欧美性猛交xxxx乱大交极品| 性视频一区| 国产精品v欧美| 亚洲天堂久久久| 91福利国产成人精品导航| 少妇人妻无码首页| 久久久久久久蜜桃| 久久永久精品免费视频| 中文字幕久久亚洲一区| 不卡午夜视频| 免费一看一级毛片| 无码福利视频| 天天躁夜夜躁狠狠躁躁88| 无码综合天天久久综合网| 自慰网址在线观看| 97人人做人人爽香蕉精品| 91欧美在线| 亚洲av无码久久无遮挡| 婷婷丁香色| 国产精女同一区二区三区久| 国产精品亚洲а∨天堂免下载| 国产福利2021最新在线观看| 国产另类视频| 国产swag在线观看| 精品一区二区三区视频免费观看| 精品久久久久成人码免费动漫| 亚洲无码高清一区| 玖玖精品在线| 国产素人在线| 国产精品美人久久久久久AV| 九九这里只有精品视频| 人妻少妇久久久久久97人妻| Jizz国产色系免费| 色妞www精品视频一级下载| 中文国产成人精品久久一| 国产97视频在线| 就去色综合| 91在线无码精品秘九色APP| 91麻豆精品国产91久久久久| 免费人成在线观看成人片| 亚洲无码不卡网| 国产亚洲精| 亚洲综合极品香蕉久久网| 专干老肥熟女视频网站| 精品剧情v国产在线观看| 久久永久免费人妻精品| 无码专区国产精品一区| 99久久精品视香蕉蕉| 日韩精品欧美国产在线| 成年午夜精品久久精品| 国产aⅴ无码专区亚洲av综合网| 日韩二区三区无| 漂亮人妻被中出中文字幕久久| 中国黄色一级视频| 亚洲综合精品香蕉久久网| 青青青伊人色综合久久| 久久永久视频| 国产视频欧美| 国产91色| 青青青国产精品国产精品美女| 九色最新网址| 园内精品自拍视频在线播放| 538国产在线| 日本尹人综合香蕉在线观看| 国产在线视频导航| 97国产精品视频人人做人人爱| 91毛片网| 91网在线|