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

        ?

        異種鋼激光-電弧焊復(fù)合焊接數(shù)值模擬

        2015-06-26 13:03:39許新猴趙小強(qiáng)華鵬李先芬周偉
        精密成形工程 2015年4期
        關(guān)鍵詞:異種應(yīng)力場(chǎng)熱源

        許新猴,趙小強(qiáng),華鵬,李先芬,周偉,2

        (1.合肥工業(yè)大學(xué)材料科學(xué)與工程學(xué)院,合肥230009; 2.新加坡南洋理工大學(xué)機(jī)械與宇航工程學(xué)院,新加坡639798)

        異種鋼激光-電弧焊復(fù)合焊接數(shù)值模擬

        許新猴1,趙小強(qiáng)1,華鵬1,李先芬1,周偉1,2

        (1.合肥工業(yè)大學(xué)材料科學(xué)與工程學(xué)院,合肥230009; 2.新加坡南洋理工大學(xué)機(jī)械與宇航工程學(xué)院,新加坡639798)

        目的 研究異種鋼激光-GMAW復(fù)合焊接溫度場(chǎng)以及應(yīng)力場(chǎng)變化。方法 運(yùn)用ANSYS有限元分析軟件,以5 mm厚D500鋼和A514鋼為研究對(duì)象,采用均勻分布的柱體熱源與橢球熱源組合的方法,建立了激光-GMAW焊接熱源模型,對(duì)異種鋼激光電弧復(fù)合焊接過(guò)程進(jìn)行了模擬計(jì)算,并與實(shí)驗(yàn)所得的焊縫形狀以及焊后殘余應(yīng)力進(jìn)行了對(duì)比。結(jié)果 結(jié)果表明,異種鋼激光電弧復(fù)合焊接過(guò)程焊接變形以及殘余應(yīng)力實(shí)驗(yàn)結(jié)果與數(shù)值計(jì)算結(jié)果吻合較好。結(jié)論 驗(yàn)證了錐體加柱體熱源與橢球熱源的組合熱源模型在異種鋼激光-GMAW復(fù)合焊接溫度場(chǎng)及應(yīng)力場(chǎng)模擬中的適用性,從而為不同焊接工藝條件下異種鋼激光-GMAW復(fù)合焊接的焊縫形狀和尺寸預(yù)測(cè),提供了一種有效的途徑。

        異種鋼;熱源模型;激光-GMAW復(fù)合焊;應(yīng)力場(chǎng);數(shù)值模擬

        激光電弧復(fù)合焊接通過(guò)激光和電弧兩熱源之間的相互作用,彌補(bǔ)了單一熱源焊接工藝的不足,具有焊接熔深大、焊接速度快、工件變形小、熔池搭橋能力強(qiáng)、易于集成等特點(diǎn)[1—4]。隨著社會(huì)的進(jìn)步和科學(xué)技術(shù)的發(fā)展,人們?cè)絹?lái)越重視焊接質(zhì)量和焊接生產(chǎn)效率。有限元技術(shù)和焊接技術(shù)的飛速發(fā)展,為數(shù)值模擬技術(shù)提供了有力的工具,焊接溫度場(chǎng)的研究和殘余應(yīng)力的分布情況可以采用數(shù)值模擬方法進(jìn)行分析,這樣就可以省去大量的實(shí)驗(yàn),從而可以大大節(jié)省人力、物力和時(shí)間,具有很大的經(jīng)濟(jì)效益[5—7]。D500鋼與A514鋼作為一種低合金高強(qiáng)鋼,具有高強(qiáng)度、高韌性的特點(diǎn),被廣泛應(yīng)用于船舶工業(yè)、橋梁鋼構(gòu)、航空航天等領(lǐng)域。低合金鋼在世界范圍內(nèi)需求的不斷提高使得其在焊接結(jié)構(gòu)件上的應(yīng)用越來(lái)越普遍,焊接質(zhì)量和焊接效率要求也越來(lái)越高。目前,針對(duì)異種鋼激光電弧復(fù)合焊接數(shù)值模擬研究較少,因此,異種鋼焊接過(guò)程殘余應(yīng)力及變形的模擬研究對(duì)實(shí)際生產(chǎn)具有一定的指導(dǎo)意義。文中以厚度均為5 mm的D500鋼與A514鋼為研究對(duì)象,借助ANSYS有限元分析軟件,對(duì)異種鋼激光-GMAW復(fù)合焊接溫度場(chǎng)及應(yīng)力場(chǎng)進(jìn)行模擬計(jì)算,并與實(shí)驗(yàn)所得的焊縫形狀以及焊后殘余應(yīng)力進(jìn)行對(duì)比。

        1 激光-電弧焊(GMAW)復(fù)合焊接試驗(yàn)

        試驗(yàn)選用厚度均為5 mm的D500鋼板和A514鋼板,使用激光+GMAW復(fù)合熱源焊機(jī)頭焊接,復(fù)合焊接裝置如圖1所示,其中,激光器為額定功率3.8 kW的Nd:YAG固體激光器,激光波長(zhǎng)為1.06μm,聚焦透鏡焦距為300 mm,焦點(diǎn)直徑為0.3 mm;脈沖GMAW焊接電源為FroniusTPSS000數(shù)字化電源,對(duì)于一定材料和直徑的焊絲,焊接電流、電弧電壓、送絲速度等焊接規(guī)范參數(shù)可實(shí)現(xiàn)一體化調(diào)節(jié)。焊接時(shí),激光前置,GMAW焊槍在后,兩者夾角為10°,光絲間距為3 mm,激光功率為3000 W;電弧功率為3800 W,焊接速度為20 mm/s。焊接試驗(yàn)結(jié)束后,對(duì)試樣進(jìn)行線切割、打磨、拋光、腐蝕,獲取接頭截面宏觀試樣,測(cè)量焊縫表面尺寸,并利用X射線法對(duì)焊接接頭殘余應(yīng)力進(jìn)行測(cè)定。

        圖1 激光-GMAW復(fù)合焊接裝置Fig.1 Set-up for Laser-GMAW welding

        2 激光-GMAW復(fù)合焊接模型

        2.1 控制方程及邊界條件

        由焊接過(guò)程的特點(diǎn)可知,焊接溫度場(chǎng)分析屬于典型的非線性瞬態(tài)熱傳導(dǎo)問(wèn)題,非線性瞬態(tài)熱傳導(dǎo)問(wèn)題的控制方程[8]為:

        式中:Cρ為定壓比熱容;λ是溫度的函數(shù);Hr為熔化潛熱。

        異種鋼激光-GMAW復(fù)合焊接有限元模型如圖2所示。

        圖2 復(fù)合焊接有限元模型Fig.2 The finite elementmodel for hybrid welding

        對(duì)流邊界條件的處理比較復(fù)雜,施加的氣體保護(hù)焊接過(guò)程包括4種對(duì)流現(xiàn)象,即熔池中液態(tài)金屬的對(duì)流,遠(yuǎn)離熔池區(qū)域的自然對(duì)流,噴嘴下方區(qū)域的保護(hù)氣體強(qiáng)制對(duì)流,以及噴嘴鄰近區(qū)域的保護(hù)氣體流散所造成的強(qiáng)制對(duì)流。在數(shù)值模擬中,熔池內(nèi)的對(duì)流損失一般采用人為的提高導(dǎo)熱系數(shù)來(lái)處理。為簡(jiǎn)化整個(gè)計(jì)算過(guò)程,采用一個(gè)總傳熱系數(shù),即20 W/(m2·K),試樣初始溫度設(shè)為環(huán)境溫度(20℃)。

        2.2 激光-GMAW復(fù)合焊接熱源模型

        激光-GMAW復(fù)合熱源焊接過(guò)程中,由于激光和電弧的相互作用,形成一個(gè)溫度極高的熱源(約有11 000 K),對(duì)工件有強(qiáng)烈的加熱作用。對(duì)于激光-GMAW復(fù)合焊接,熱源模型總體上是2種子熱源(電弧、激光)的組合。此外,由于復(fù)合焊接中電弧利用率得以有效提高,等離子體能量在變形熔池內(nèi)部分布較為復(fù)雜。基于激光-電弧能量耦合機(jī)理,復(fù)合熱源模型建立在橢球熱源加錐體熱源加柱體熱源模型基礎(chǔ)之上,通過(guò)調(diào)整相應(yīng)的分布參數(shù),模擬激光和電弧的耦合作用。組合式熱源模型對(duì)于復(fù)合焊而言,相當(dāng)于有一個(gè)橢球分布的體熱源模擬電弧的熱輸入作用,運(yùn)用錐體熱源及柱體熱源來(lái)體現(xiàn)激光焊接典型的“釘頭”和“深熔”的現(xiàn)象。各熱源模型公式見(jiàn)參考文獻(xiàn)[9—12]。

        2.3 材料熱物理性能參數(shù)

        材料的熱物理性能對(duì)激光-電弧復(fù)合焊的溫度場(chǎng)分布及焊縫成形有著非常重要的影響,其物理性能參數(shù)直接影響溫度場(chǎng)和應(yīng)力場(chǎng)的形態(tài)和大小[13]。在對(duì)異種鋼激光-GMAW復(fù)合焊模擬分析時(shí)必須確定以下熱物理性能參數(shù):焊件的初始溫度、焊件的熔點(diǎn)(℃)、材料的密度(kg/m3)、導(dǎo)熱系數(shù)(W/(m·℃))、對(duì)流換熱系數(shù)(W/(m2·℃))和比熱容(J/(kg·℃))。除了材料的熔點(diǎn)、工件的初始溫度(20℃)、材料的密度是常數(shù)外,材料的其他屬性參數(shù)都是溫度的函數(shù)。D500鋼與A514鋼物理性能參數(shù)[14—16]見(jiàn)表1—3。

        表1 A514鋼應(yīng)力場(chǎng)分析熱物理部分參數(shù)Tab le 1 Part therm al param eters of A514 steel

        表2 A514鋼應(yīng)力場(chǎng)分析熱物理部分參數(shù)Tab le 2 Part therm al param eters ofA 514 steel

        表3 D500鋼應(yīng)力場(chǎng)分析熱物理性能參數(shù)Tab le 3 Part therm al param eters of D500 steel

        3 結(jié)果及分析

        3.1 溫度場(chǎng)結(jié)果分析及實(shí)驗(yàn)驗(yàn)證

        圖3所示為2.5 s時(shí)焊件上表面溫度場(chǎng)分布以及1.5 s時(shí)縱截面方向熱循環(huán)曲線圖。

        圖3 溫度場(chǎng)分布Fig.3 Temperature distribution on the surface

        由圖3a可以看出,激光-GMAW復(fù)合焊接過(guò)程中熱影響區(qū)小,焊縫附近的溫度場(chǎng)梯度大;熱源后方的等溫線稀疏,溫度梯度小;在熱源前方的等溫線密,溫度梯度大。隨著熱源的移動(dòng),相鄰位置的金屬經(jīng)歷不同的熱循環(huán)史,熱源后方的熔融金屬開(kāi)始冷卻凝固并逐漸形成焊縫,而恢復(fù)彈性對(duì)后冷卻的焊縫金屬產(chǎn)生拘束作用,這是誘導(dǎo)殘余應(yīng)力及焊接變形的一個(gè)重要原因。圖3b所示為縱截面方向在1.5 s時(shí)刻的熱循環(huán)曲線。由圖3可以看出,焊接時(shí)試樣經(jīng)歷了一個(gè)快速升溫和緩慢降溫的過(guò)程,且最高溫度已經(jīng)遠(yuǎn)遠(yuǎn)高于材料的熔點(diǎn),由于2種材料的熱膨脹系數(shù)不同,容易在焊縫周圍產(chǎn)生過(guò)大的應(yīng)力應(yīng)變,接頭容易產(chǎn)生缺陷。

        圖4為模擬焊縫形狀與實(shí)驗(yàn)焊縫形狀的對(duì)比。實(shí)際焊縫深度為5 mm,上下表面熔寬分別為5.4,1.3 mm,有限元模擬結(jié)果分別為5.0,5.0,1.4 mm。

        由圖4和數(shù)據(jù)對(duì)比結(jié)果可知,組合式熱源模型模擬所得的焊縫形狀與實(shí)驗(yàn)所得的焊縫形狀比較接近,且2種結(jié)果的焊縫尺寸在誤差范圍內(nèi)。

        圖4 模擬焊縫形狀與實(shí)驗(yàn)焊縫形狀的對(duì)比Fig.4 Comparison of simulated and experimental weld seam shape

        3.2 應(yīng)力場(chǎng)結(jié)果分析及實(shí)驗(yàn)驗(yàn)證

        圖5所示為實(shí)驗(yàn)結(jié)果與模擬結(jié)果縱向殘余應(yīng)力對(duì)比。圖5b中實(shí)測(cè)的殘余應(yīng)力分布由X射線法測(cè)得。由圖5可以看出,焊縫處數(shù)值模擬出的焊接殘余應(yīng)力基本上能夠與試驗(yàn)測(cè)量的殘余應(yīng)力結(jié)果相對(duì)應(yīng),曲線變化趨勢(shì)基本一致,實(shí)測(cè)結(jié)果與模擬結(jié)果中的最大殘余拉應(yīng)力均出現(xiàn)在焊縫兩側(cè)的熱影響區(qū),且縱向殘余應(yīng)力主要表現(xiàn)為兩端受較小的拉應(yīng)力,越靠近焊縫中間處拉應(yīng)力越大,而在焊縫中心處應(yīng)力略微下降。因?yàn)楹缚p部分的材料最后冷卻,此時(shí)焊縫縱向收縮時(shí)受到已經(jīng)冷卻材料的阻力,焊縫冷卻受阻,焊縫及其附近區(qū)域的縱向應(yīng)力為拉應(yīng)力。模擬結(jié)果與實(shí)測(cè)結(jié)果中,焊縫左側(cè)的最大殘余應(yīng)力大于焊縫右側(cè)的殘余應(yīng)力,這是由于在同一焊接熱循環(huán)下,A514鋼的熱膨脹系數(shù)比D500鋼的熱膨脹系數(shù)要大,在焊縫冷卻時(shí)產(chǎn)生的收縮較大,致使A514鋼一側(cè)的殘余應(yīng)力較大。

        通過(guò)對(duì)模擬所得應(yīng)力曲線和實(shí)驗(yàn)測(cè)得的應(yīng)力曲線進(jìn)行比較分析可知,模擬的應(yīng)力分布曲線與實(shí)驗(yàn)測(cè)得應(yīng)力分布曲線形狀相近,但是模擬得出的殘余應(yīng)力數(shù)值比其實(shí)驗(yàn)得出的殘余應(yīng)力數(shù)值要小。引起殘余應(yīng)力值偏小的原因可能是在進(jìn)行溫度場(chǎng)及應(yīng)力場(chǎng)模擬時(shí),對(duì)材料的熱物理性能參數(shù)選擇不當(dāng),建模后的網(wǎng)格劃分疏密程度以及模擬計(jì)算時(shí)時(shí)間載荷步的設(shè)置方式等。

        圖5 實(shí)測(cè)殘余應(yīng)力與模擬所得殘余應(yīng)力的對(duì)比Fig.5 Comparison of experimental and simulated residual stress

        4 結(jié)論

        1)采用錐體-橢球-柱體組合式熱源模型,可以對(duì)5 mm厚的A514與D500異種鋼激光-GMAW復(fù)合焊的溫度場(chǎng)進(jìn)行準(zhǔn)確模擬,模擬所得焊縫形狀和實(shí)際焊縫截面形狀比較接近。

        2)通過(guò)應(yīng)力場(chǎng)模擬結(jié)果可以看出,實(shí)測(cè)結(jié)果與模擬結(jié)果中的最大殘余拉應(yīng)力均出現(xiàn)在焊縫兩側(cè)的熱影響區(qū),且縱向殘余應(yīng)力主要表現(xiàn)為兩端受較小的拉應(yīng)力,越靠近焊縫中間處拉應(yīng)力越大,而在焊縫中心處應(yīng)力略微下降,且A514鋼一側(cè)的殘余應(yīng)力大于D500鋼一側(cè)的殘余應(yīng)力。

        3)由試驗(yàn)結(jié)果與模擬結(jié)果對(duì)比分析可知,數(shù)值模擬結(jié)果與測(cè)量的焊接殘余應(yīng)力分布規(guī)律曲線形狀基本一致,但模擬結(jié)果數(shù)值偏小。

        [1]田川,聞義.激光電弧復(fù)合焊接——一種新型焊接方法[J].機(jī)車車輛工藝,2005,4(2):5—10. TIAN Chuan,WEN Yi.A New Welding Process:Laser—Arc Hybrid Welding[J].Locomotive&Rolling Stock Technology,2005,4(2):5—10.

        [2]盧振洋,劉建,黃鵬飛,等.基于ANSYS薄板GTAW焊接溫度場(chǎng)數(shù)值模擬[J].微計(jì)算機(jī)信息,2007,23(4):291—292. LU Zhen-yang,LIU Jian,HUANG Peng-fei,etal.Simulation on Welding Temperature Field of GTAW on ANSYS[J]. Microcomputer Information,2007,23(4):291—292.

        [3]FENG Gao,QIAN Yi-yu.Micromechanical Properties of Heterogeneous Aluminum-silicon Brazed Joint[J].Materials Letters,2004(58):2861—2866.

        [4]KELVIIWei Guo.Influence of In Situ Reaction on the Microstructure of SiCp/AlSi7Mg Welded by Nd:YAG Laser with Ti Filler[J].JMEPEG,2010(19):52—58.

        [5]崔旭明,李劉合,張彥華.激光-電弧復(fù)合熱源焊接[J].焊接技術(shù),2003,32(2):19—21. CUIXu-min,LILiu-he,ZHANG Yan-hua.Laser-arc Hybrid Welding[J].Welding Technology,2003,32(2):19—21.

        [6]武傳松.焊接熱過(guò)程與熔池形狀[M].北京:機(jī)械工業(yè)出版社,2007. WU Chuan-song.Welding Thermal Process and Pool Shape [M].Beijing:Machinery Industry Press,2007.

        [7]LEEW B,YEON Y M,JUNG SB.The ImprovementofMechanical Properties of Friction Stir Welded A356 Al Alloy [J].Materials Science and Engineering,2003(355):154—159.

        [8]李冬林.基于ANSYS軟件焊接溫度場(chǎng)應(yīng)力場(chǎng)模擬研究[J].湖北工業(yè)大學(xué)學(xué)報(bào),2005,20(5):81—83. LIDong-lin.Simulation Research for Welding Temperature and Stress Field Based on ANSYS[J].Journal of Hubei U-niversity of Technology,2005,20(5):81—83.

        [9]張旭東.CO2激光-MIG同軸復(fù)合焊方法及鋁合金焊接的研究[J].應(yīng)用激光,2005,25(1):1—3. ZHANG Xu-dong.Coaxial Hybrid CO2Laser-MIGWelding System and Its Application in Welding of Aluminum[J]. Applied Laser,2005,25(1):1—3.

        [10]AMIRIZADM,KOKABIA H,GHARACHEH MA,etal.E-valuation of Microstructure and Mechanical Properties in Friction Stir Welded A356+15%SiCp Cast Composite[J]. Materials Letters,2006(60):565—568.

        [11]XU C L,JIANG Q C.Morphologies of Primary Silicon in Hypereutectic Al—Sialloyswith MeltOverheating Temperature and Cooling Rate[J].Materials Science and Engineering A,2006(437):451—455.

        [12]張春平,李萌盛,秦琳.基于ANSYS的異種鋼焊接殘余應(yīng)力的數(shù)值模擬[J].合肥工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2003,26(3):460—463. ZHANG Chun-ping,LI Meng-sheng,QIN Lin.Numerical Simulation of Residual Stress in Dissimilar Steel Based on ANSYS[J].Journal of Hefei University of Technology (Natural Science),2003,26(3):460—463.

        [13]CHEN Ming-hua,LI Chen-bin,LIU Li-ming.Research of Discharge Status of Arc Plasma Enhanced by Laser Enhanced Research[J].High Voltage Engineering,2013,39 (7).

        [14]吳圣川,劉建華.鋁合金激光-電弧復(fù)合焊的有限元數(shù)值[J].模擬航空制造技術(shù),2005(12):74—76. WU Shen-chuang,LIU Jian-hua.Finite Element Simulation on Laser-arc Hybrid Welding of Aluminum Alloy[J].Aeronautical Manufacturing Technology,2005(12):74—76.

        [15]GHOSH M,KUMAR K,KAILASS V,et al.Optimization of Friction Stir Welding Parameters for Dissimilar Aluminum Alloys[J].Materials and Design,2010(31):3033—3037.

        [16]MA Zhi-hua,CHEN Dong-gao,TAN Bin.Effect of Laser-MIG Hybrid Welding of5052 Alloy on Weld Forming[J]. Material Science and Engineering,2012,2(35):76—80.

        Simulation on Laser-GMAW Hybrid W elding of Heterogeneous Steel

        XU Xin-hou1,ZHAO Xiao-qiang1,HUA Peng1,LIXian-fen1,ZHOUWei1,2
        (1.School of Material Science and Engineering,Hefei University of Technology,Hefei230009,China; 2.School of Mechanical and Aerospace Engineering,Nan yang Technological University,Singapore 639798,Singapore)

        on the aim of this study was to investigate thewelding thermal field and residual stress field of laser-GMAW hybrid welding.ANSYS finite element softwarewas employed to simulate the welding thermal field and residual stress field of Heterogeneous steel.Laser-GMAW hybrid welding sourcemodelwas built by using Ellipsoid-Cylinder-Cone hybrid heat method,and the results were analyzed.The simulated results were coincided well to the experimental results.The laser-GMAW hybrid welding sourcemodel is applicable in the simulation ofwelding thermal field and residual stress field of Heterogeneous steel,which provided an effectivemethod for predicting weld seam shape and size of heterogeneous steel under differentwelding conditions.

        heterogeneous steel;heat sourcemodel;laser-GMAW hybrid welding;welding stress field;simulation

        10.3969/j.issn.1674-6457.2015.04.015

        TG456.7

        :A

        :1674-6457(2015)04-0071-05

        2015-05-25

        許新猴(1990—),男,安徽,碩士研究生,主要研究方向?yàn)榧す怆娀?fù)合焊接技術(shù)。

        李先芬(1969—),女,安徽人,博士,教授,碩士生導(dǎo)師,主要從事先進(jìn)焊接工藝及焊接材料研究。

        猜你喜歡
        異種應(yīng)力場(chǎng)熱源
        橫流熱源塔換熱性能研究
        煤氣與熱力(2021年3期)2021-06-09 06:16:20
        基于啟發(fā)式動(dòng)態(tài)規(guī)劃的冷熱源優(yōu)化控制
        電子制作(2017年19期)2017-02-02 07:08:31
        異種部門(mén)
        異種部門(mén)
        異種部門(mén)
        異種部門(mén)
        鋁合金多層多道窄間隙TIG焊接頭應(yīng)力場(chǎng)研究
        焊接(2016年9期)2016-02-27 13:05:22
        中部槽激光-MAG復(fù)合熱源打底焊焊接工藝研究
        焊接(2015年8期)2015-07-18 10:59:13
        多類型熱源點(diǎn)共存下的區(qū)域熱力統(tǒng)籌供應(yīng)探討
        考慮斷裂破碎帶的丹江口庫(kù)區(qū)地應(yīng)力場(chǎng)與水壓應(yīng)力場(chǎng)耦合反演及地震預(yù)測(cè)
        国产白嫩护士被弄高潮| 亚洲av免费看一区二区三区 | 欧美亚洲国产片在线播放| 国产美熟女乱又伦av果冻传媒| 亚洲乱码中文字幕综合| 国产一区二区三区色哟哟| 人妻丰满熟妇无码区免费| 国产成人www免费人成看片| 日本一道dvd在线中文字幕| 精品国产免费一区二区久久| 免费人成小说在线观看网站 | 美女高潮无遮挡免费视频| 亚州五十路伊人网| 偷拍韩国美女洗澡一区二区三区 | 粉嫩的18在线观看极品精品| 国产精品人成在线观看免费| 女邻居的大乳中文字幕| 未满十八勿入av网免费| 亚洲成av人片在久久性色av| 欧美熟妇另类久久久久久多毛| 麻豆精品传媒一二三区| 精品熟女少妇免费久久| 国产精品女同一区二区免| 日本动漫瀑乳h动漫啪啪免费| 国产熟妇搡bbbb搡bb七区| a√无码在线观看| 日韩一区二区三区精品视频| 色天使综合婷婷国产日韩av | 亚洲动漫成人一区二区| 麻豆国产精品伦理视频| 久久天堂av综合合色| 欧美老熟妇欲乱高清视频| 玩弄人妻奶水无码AV在线| 丁香婷婷六月综合缴清| 亚欧中文字幕久久精品无码| 欧美巨大精品欧美一区二区| 在线日韩中文字幕乱码视频| 国产精品国产三级国产aⅴ下载| 乱人伦视频中文字幕| 国产亚洲美女精品久久久2020| 久久精品国产亚洲av性瑜伽|