• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    激光-電弧復(fù)合熱源及等離子體模擬研究現(xiàn)狀

    2016-12-13 06:53:32丁玲芳張英喬劉盛耀
    電焊機(jī) 2016年3期
    關(guān)鍵詞:熔池熱源電弧

    丁玲芳,張英喬,劉盛耀

    (中北大學(xué)材料科學(xué)與工程學(xué)院,山西太原030051)

    激光-電弧復(fù)合熱源及等離子體模擬研究現(xiàn)狀

    丁玲芳,張英喬,劉盛耀

    (中北大學(xué)材料科學(xué)與工程學(xué)院,山西太原030051)

    激光-電弧復(fù)合焊作為一種新型高效復(fù)合焊接工藝,受到國(guó)內(nèi)外學(xué)者的廣泛重視。綜述了激光-電弧復(fù)合焊熱源耦合機(jī)理及數(shù)值模擬熱源模型,詳細(xì)介紹了復(fù)合熱源模型的應(yīng)用及等離子體的模擬現(xiàn)狀,指出了建立激光-電弧復(fù)合焊統(tǒng)一的數(shù)學(xué)模型是未來(lái)研究的新方向。

    激光-電弧復(fù)合;熱源模型;等離子體;數(shù)值模擬

    0 前言

    激光-電弧復(fù)合熱源焊接技術(shù)由Steen于1979年首次提出[1],20世紀(jì)90年代以來(lái)得到迅速發(fā)展,其兼具快速、經(jīng)濟(jì)及焊接缺陷少等優(yōu)點(diǎn),已廣泛應(yīng)用于汽車、造船、焊管等各個(gè)領(lǐng)域,但對(duì)于激光-電弧復(fù)合熱源的研究并不充分,特別是在數(shù)值模擬方面[2-3]。

    為了更好地認(rèn)識(shí)激光-電弧復(fù)合焊接過(guò)程中激光與電弧兩種不同熱源間作用機(jī)理,需借助數(shù)值模擬技術(shù),提高對(duì)激光-電弧復(fù)合焊復(fù)雜物理過(guò)程的認(rèn)識(shí),從而進(jìn)一步優(yōu)化工藝參數(shù),提高焊接接頭的性能。在此分析了激光-電弧復(fù)合熱源及等離子體數(shù)值模擬研究現(xiàn)狀,具有一定的理論意義。

    1 激光-電弧復(fù)合作用機(jī)理

    激光-電弧復(fù)合焊是由物理性質(zhì)與能量傳輸機(jī)制截然不同的兩種熱源復(fù)合進(jìn)行焊接[4]。復(fù)合熱源相互耦合并作用于同一加熱部位,既能發(fā)揮各自優(yōu)勢(shì),又可彌補(bǔ)自身不足,其能量大于兩熱源的簡(jiǎn)單疊加[5-6]。

    激光-電弧復(fù)合焊中,電弧屬低溫、低密度等離子體,而激光等離子體則有高溫、高密度的特點(diǎn)。激光與電弧應(yīng)保持合適的間距,并產(chǎn)生一個(gè)傳輸帶電粒子的導(dǎo)電通道,從而發(fā)生強(qiáng)烈的交互作用[7]。激光和電弧等離子體相互作用機(jī)制為:①激光能穩(wěn)定電弧、防止電弧漂移。激光與電弧作用于焊接部位,激光束作用處產(chǎn)生金屬蒸汽,電弧電壓趨于穩(wěn)定,場(chǎng)強(qiáng)降低,弧柱電阻減小,激光對(duì)電弧產(chǎn)生明顯的吸引、壓縮作用,增加了電弧的穩(wěn)定性[8-10]。Albright等人研究低能量激光輔助電弧焊時(shí),發(fā)現(xiàn)保護(hù)氣分子吸收特定波長(zhǎng)激光的能量,相互碰撞并伴有微弱的氣體電離,激發(fā)激光等離子?。↙IP),擴(kuò)展電弧導(dǎo)電通路,從而提高電弧穩(wěn)定性[11]。②電弧能強(qiáng)化激光,稀釋等離子體,提高激光吸收率;但電弧等離子體的“負(fù)透鏡效應(yīng)”使激光束產(chǎn)生偏折和散焦,影響激光吸收[4]。電弧稀釋激光等離子體,降低等離子體對(duì)激光能量的吸收、反射,增大激光穿透能力,提升激光能量傳輸率;然而,激光束通過(guò)等離子體時(shí),會(huì)產(chǎn)生逆韌致輻射吸收,導(dǎo)致激光能量衰減,衰減值與激光波長(zhǎng)大小有關(guān)。波長(zhǎng)較長(zhǎng)的CO2激光(10.6 μm)的能量衰減大于波長(zhǎng)較短的Nd:YAG激光(1.06 μm)[12-13],Hu等[14]借助功率計(jì)測(cè)量YAG激光垂直穿過(guò)TIG電弧后的功率衰減,證實(shí)YAG激光穿過(guò)電弧等離子體的能量損耗很低,可以忽略。Bibik等人[15]研究CO2激光垂直穿過(guò)TIG電弧時(shí),發(fā)現(xiàn)激光的傳輸系數(shù)隨電弧電流的增大而減小,當(dāng)電弧電流大于200 A時(shí),電弧對(duì)激光能量的吸收接近30%,此時(shí)能量損耗不可忽略。

    2 激光-電弧復(fù)合焊熱源模型及應(yīng)用

    焊接熱源模型是焊接熱過(guò)程數(shù)值模擬的基礎(chǔ),經(jīng)歷了從單一熱源到組合熱源的演變過(guò)程。激光-電弧復(fù)合焊中,組合熱源模型分別考慮激光與電弧的熱輸入,同時(shí)考慮輻射、表面張力、電弧壓力等因素而建立,模型較復(fù)雜。隨著熱源模型的不斷改進(jìn)和完善,激光-電弧復(fù)合焊的數(shù)值模擬體系也逐漸完善,不僅能計(jì)算焊接溫度場(chǎng)分布情況,而且能有效預(yù)測(cè)焊后殘余應(yīng)力分布及變形情況,還能觀察熔池流動(dòng)及合金元素分布,有助于焊接參數(shù)的優(yōu)化設(shè)計(jì)。以下分別分析激光-電弧復(fù)合焊中電弧及電弧復(fù)合激光后的熱源模型與應(yīng)用。

    2.1電弧部分熱源模型

    復(fù)合焊中電弧部分熱源模型主要討論分析鎢極惰性氣體保護(hù)焊TIG和熔化極氣體保護(hù)焊GMAW。TIG焊適于薄板焊接,且不消耗電極材料,焊接過(guò)程相對(duì)簡(jiǎn)單。因此,TIG電弧熱源作用的表面熱流通常采用平面熱源,如高斯熱源、雙橢圓平面熱源[16-18],其中高斯熱源模型為

    式中r0為熱源有效半徑;Q為熱源有效功率;η為焊接熱效率;U為電弧電壓;Iarc為焊接電流。

    雙橢圓平面熱源模型考慮了熱源移動(dòng)時(shí)焊縫前后能量分布的差異,其函數(shù)表達(dá)式為

    式中qf、qr分別為前、后半橢圓內(nèi)的熱流分布;Qf、Qr分別為電弧前、后的總熱量;a1、a2、b為雙橢圓熱源參數(shù)。

    GMAW焊接過(guò)程不僅受電弧壓力的影響,還存在熔滴過(guò)渡現(xiàn)象,情況復(fù)雜,簡(jiǎn)化處理方式也存在差異。當(dāng)僅考慮電弧等離子體時(shí),由Goldak提出的“雙橢球體熱源模型”[19]被廣泛應(yīng)用,此模型從宏觀傳熱學(xué)角度,以焊縫橫斷面的形貌輪廓為參照,不僅考慮工件厚度方向的熱流分布,還考慮熱源移動(dòng)時(shí)焊接前后熱流分布的不對(duì)稱性,將前、后半部分各分成兩個(gè)不同的1/4橢球,兩半部分橢球內(nèi)熱源分布為

    式中qf、qr分別為前、后半橢球體內(nèi)的熱流分布;ar、af、bh、ch為前后半橢球體的半軸;ff、fr分別為前后半橢球體內(nèi)熱輸入的比份,且ff+fr=1。

    當(dāng)將電弧熱源分為電弧等離子體和熔滴兩部分時(shí),電弧等離子體多采用平面熱源模型,作用于焊件表面,熔滴則選用形狀各異的體熱源模型,表示在熔滴沖擊作用下導(dǎo)致的焊件表面下陷。如有學(xué)者研究激光-GMAW復(fù)合焊接工藝時(shí),電弧等離子體部分采用雙橢圓平面熱源,熔滴部分采用雙橢球體熱源[18];或選用表面高斯熱源模擬氣-金屬界面電弧等離子體熱流分布,圓柱狀熱源用于模擬熔滴沖擊力[20],充分考慮了焊件表面所受的電弧作用力,與實(shí)際焊接過(guò)程中焊件表面的受力情況接近。

    2.2激光-TIG/GMAW復(fù)合熱源模型及應(yīng)用

    激光-電弧復(fù)合熱源模型是根據(jù)焊接結(jié)構(gòu)、焊接工藝參數(shù)、母材等的差異,合理地將激光和電弧熱源模型組合起來(lái),通過(guò)數(shù)值模擬進(jìn)行多物理場(chǎng)耦合,實(shí)現(xiàn)兩者的交互作用。復(fù)合焊中激光的加入主要考慮通過(guò)“小孔”結(jié)構(gòu)來(lái)實(shí)現(xiàn)能量的轉(zhuǎn)換機(jī)制,增加焊件的熔深。以激光熱源模型發(fā)展為主線,介紹復(fù)合熱源模型及其應(yīng)用。

    激光-電弧復(fù)合焊中,針對(duì)激光產(chǎn)生的小孔結(jié)構(gòu),多采用體熱源模型,如較早提出的高斯圓柱體熱源模型,激光能量沿徑向呈高斯分布。Kong等人[21]結(jié)合實(shí)驗(yàn)所得的焊縫橫截面形狀研究激光-GMAW復(fù)合焊接A36低碳鋼時(shí),建立雙橢球體熱源+高斯圓柱體熱源模型分別模擬GMAW熱輸入和激光熱輸入,得到焊后殘余應(yīng)力的分布云圖,發(fā)現(xiàn)在試樣處存在明顯的殘余應(yīng)力。計(jì)算過(guò)程中,雖然通過(guò)提高焊接熔點(diǎn)處的熱傳導(dǎo)率來(lái)彌補(bǔ)熔池內(nèi)因流體流動(dòng)引起的對(duì)流因素,但仍不能代替對(duì)流散熱的復(fù)雜性,而且模型中并未考慮激光能量在小孔內(nèi)傳輸過(guò)程中的損耗,存在一定的局限性。

    針對(duì)激光熱源產(chǎn)生的典型的“釘頭”形狀及“深熔”現(xiàn)象,由吳甦提出的旋轉(zhuǎn)高斯體熱源模型[22],綜合考慮由激光引發(fā)的小孔機(jī)制及匙孔壁的能量傳遞的物理特點(diǎn),其作為一種新型高能束焊接熱源,得到廣泛應(yīng)用,如圖1a所示,此熱源模型的數(shù)學(xué)表達(dá)式為

    Liang等人[23]對(duì)AZ31B鎂合金進(jìn)行激光-TIG復(fù)合焊接模擬時(shí),采用高斯+旋轉(zhuǎn)高斯體組合熱源分別模擬TIG電弧表面熱流和激光熱源,能得到呈“螺釘”狀的焊縫橫斷面,但忽略了焊件表面受到的表面張力、電弧壓力、磁懸浮力等作用,獲得的熔池表面半徑的模擬結(jié)果小于實(shí)際結(jié)果。

    李迅波等人[17]對(duì)鎂合金/不銹鋼板進(jìn)行異種材料對(duì)接焊時(shí)的激光-TIG復(fù)合焊研究,不僅用雙橢圓分布的面熱流密度和旋轉(zhuǎn)高斯體熱源分別模擬電弧和激光熱輸入,還采用雙橢球體積分布的熱流密度模擬因電弧沖擊等作用導(dǎo)致的變形熔池內(nèi)部激光能量分布,彌補(bǔ)了先前學(xué)者僅單純考慮激光和電弧熱源的不足,并采用回歸分析法確定三者之間的能量分配關(guān)系,對(duì)比數(shù)值計(jì)算的焊縫溫度及殘余應(yīng)力分布與實(shí)測(cè)值,發(fā)現(xiàn)焊接熱循環(huán)曲線及殘余應(yīng)力與實(shí)測(cè)變化趨勢(shì)一致。

    此外,還有改進(jìn)的高斯體熱源模型(CIN)[24-25]、旋轉(zhuǎn)曲面體熱源模型[18-26]及呈高斯“杯狀”的新型體熱源模型[27],分別如圖1b、圖1c和圖1d所示。圖1所示四種體熱源模型既考慮到激光能量沿徑向呈高斯分布的特點(diǎn),又考慮到高能量密度的激光形成的“小孔效應(yīng)”,還涉及小孔內(nèi)激光能量損耗,并能調(diào)節(jié)熱源形狀參數(shù)以適應(yīng)不同的焊接條件。

    圖1 熱源模型Fig.1Heat source models

    CIN熱源模型指出小孔內(nèi)能量分布隨熔深增加呈指數(shù)遞減,根據(jù)焊件厚度及焊縫橫斷面輪廓來(lái)設(shè)置熱源的形狀參數(shù),獲得圓錐體狀或平截錐體狀的焊縫形狀。Piekarska和Kubiak[24-25]研究激光-GMAW復(fù)合焊熱源模擬時(shí),為避免過(guò)熱熔滴的復(fù)雜性,僅考慮電弧和激光作用,以S355鋼板為基體材料,雙橢球熱源模型用于電弧熱輸入,激光熱源則是基于實(shí)驗(yàn)研究而改進(jìn)的高斯體熱源模型(CIN),計(jì)算過(guò)程中考慮了熔池內(nèi)由懸浮力驅(qū)動(dòng)的流體流動(dòng)和熔融區(qū)材料的相變潛熱,得到焊件溫度場(chǎng)分布、焊縫橫斷面輪廓和熔池內(nèi)液態(tài)金屬的流動(dòng)速度,表明溫度場(chǎng)分布對(duì)熔池形狀有重要影響,熔池內(nèi)流動(dòng)呈渦旋狀。

    旋轉(zhuǎn)曲面體熱源模型在處理小孔內(nèi)激光的多次反射和菲涅爾吸收時(shí),采用卡普蘭法計(jì)算匙孔壁輪廓尺寸,指出小孔底部的能量密度較高,故分配到底部的峰值能量密度高于頂部。Xu等人[18]研究激光-GMAW復(fù)合焊接熱源模型的適用性時(shí),激光熱輸入選用峰值遞增的旋轉(zhuǎn)曲面體熱源,并探討了小孔內(nèi)的反射、吸收機(jī)制,電弧表面熱流采用雙橢圓平面模型,熔滴熱焓采用雙橢球體熱源。對(duì)比實(shí)驗(yàn)和數(shù)值計(jì)算的焊后6061-T6鋁合金的橫截面形狀及尺寸,發(fā)現(xiàn)預(yù)測(cè)得到的輪廓曲線較實(shí)驗(yàn)結(jié)果平滑,但輪廓形狀和尺寸大小非常接近,表明復(fù)合熱源模型的可靠性,并得到焊接溫度場(chǎng)的分布云圖。文獻(xiàn)[26]中進(jìn)一步研究T型接頭6061-T6鋁合金的焊后殘余應(yīng)力及變形情況,分析焊縫的殘余應(yīng)力云圖及變形曲線圖可知,在焊縫處及附近沿焊接方向存在較大的殘余張應(yīng)力,T型接頭前端存在小焊接角變形,中部及后端則有較大的焊接角變形。此復(fù)合熱源模型的參數(shù)與小孔的形狀和尺寸關(guān)聯(lián)性較大,而非依據(jù)試驗(yàn)數(shù)據(jù)預(yù)測(cè)得出,用于復(fù)合焊的模擬研究較精確、有效,但未考慮熔池內(nèi)復(fù)雜的液態(tài)金屬流動(dòng),與實(shí)際情況存在偏差,有待改進(jìn)和完善。

    高斯“杯狀”熱源模型則是以激光熱源在焊件表面形成的“倒喇叭”狀小孔結(jié)構(gòu)為原型,其上部是圓錐狀體熱源,下部是圓柱狀體熱源,通過(guò)調(diào)節(jié)熱源模型的形狀參數(shù)來(lái)確定上、下部的熱源形狀分界點(diǎn)。Faraji等人[27]研究激光-TIG復(fù)合焊接AA6082鋁合金時(shí),TIG熱流采用表面高斯熱源模型,激光采用高斯“杯狀”體熱源模型,還考慮了熔化潛熱,建立電磁模型表示熔池內(nèi)電弧電流和電磁力的影響,熔池內(nèi)流體流動(dòng)綜合了湍流和層流狀態(tài),涵蓋了熔池內(nèi)小孔效應(yīng)的熱傳輸和流體流動(dòng),計(jì)算所得的橫斷面輪廓與試驗(yàn)焊縫對(duì)比,發(fā)現(xiàn)考慮流體流動(dòng)后的模擬結(jié)果更接近實(shí)際焊縫輪廓,因此建立復(fù)合熱源模型時(shí)有必要同時(shí)考慮傳導(dǎo)和對(duì)流的熱傳輸機(jī)制。

    為提高計(jì)算速度,以上熱源模型均假定焊件表面為平坦界面,忽略匙孔動(dòng)力學(xué)行為。Cho等人[28]研究激光-GMAW復(fù)合焊后熔池內(nèi)合金元素分布時(shí),除計(jì)算質(zhì)量、動(dòng)量、能量守恒方程外,還有描述熔池內(nèi)液態(tài)金屬流動(dòng)的流體體積模型(VOF)方程及計(jì)算合金元素分布的附加方程,并假定熔融金屬具有牛頓粘度和不可壓縮的層流特性。此外,電弧熱流采用表面高斯分布,重力、電磁力和自感應(yīng)磁力作為體積力考慮,激光能量?jī)H考慮菲涅爾吸收,小孔內(nèi)每一束激光照射到孔壁上的能量均呈高斯分布,小孔變形主要考慮反沖壓力作用,得到熔池渦流及Cr元素分布情況,熔池內(nèi)渦流合金元素含量高,并沿同一方向旋轉(zhuǎn),使得小孔后部的合金元素均勻分布在熔池頂部。模型中雖有匙孔動(dòng)力學(xué)行為,但忽略了逆軔致輻射吸收、瑞利散射,熔融金屬和氣泡間的熱量及質(zhì)量傳輸,今后研究中還需綜合考慮這些影響因素。

    由于復(fù)合焊涉及的物理問(wèn)題較復(fù)雜,并考慮到運(yùn)算速度,這些數(shù)學(xué)模型是在簡(jiǎn)化的假設(shè)條件下建立的,有一定的局限性,但在復(fù)合焊接模擬研究中仍起著重要的指導(dǎo)作用。目前焊接熱過(guò)程的模擬研究中,盡可能少的假定條件建立三維熱源模型,從而得到較精確的模擬結(jié)果,仍是熱源模型的發(fā)展方向。

    3 激光-電弧復(fù)合焊等離子體模擬

    等離子體被稱為物質(zhì)第四態(tài),具有很強(qiáng)的導(dǎo)電性,呈現(xiàn)電準(zhǔn)中性,并受磁場(chǎng)控制,其運(yùn)動(dòng)的結(jié)果又能形成電磁場(chǎng)[29]。根據(jù)等離子體的產(chǎn)生途徑,激光-電弧復(fù)合焊中以電弧等離子體和激光等離子體為主。焊接等離子體的研究是從電弧等離子體開始,復(fù)雜的GMAW電弧物理過(guò)程,除存在電場(chǎng)、磁場(chǎng)和流場(chǎng)作用外,還有熔滴現(xiàn)象,因此GMAW電弧等離子體的研究非常少。目前就焊接電弧的數(shù)值分析大都集中于無(wú)熔滴過(guò)渡現(xiàn)象的TIG電弧,是研究激光-電弧復(fù)合焊等離子體的基礎(chǔ)。

    激光-電弧復(fù)合焊中,需考慮等離子體對(duì)激光能量的吸收、LIP及金屬蒸汽的影響[30]。有學(xué)者在TIG焊電弧模型的基礎(chǔ)上建立激光-TIG同軸復(fù)合焊數(shù)學(xué)模型[31-32],如圖2所示,并作出以下假設(shè):a.電弧等離子體和金屬蒸汽都呈軸對(duì)稱分布,滿足局部熱平衡;b.等離子體呈光學(xué)薄,僅考慮激光的逆韌致輻射;a.流體為層流流動(dòng);d.不考慮陽(yáng)極變形和熔池內(nèi)流體流動(dòng)對(duì)焊件溫度場(chǎng)的影響。在激光作用下,母材表面產(chǎn)生大量金屬蒸汽,根據(jù)金屬蒸汽的求解結(jié)果確定計(jì)算區(qū)域陽(yáng)極表面的溫度邊界條件,并求解其質(zhì)量、電量、動(dòng)量和能量四大守恒方程,經(jīng)反復(fù)迭代至收斂,得到復(fù)合焊等離子體的熱場(chǎng)、電場(chǎng)和流場(chǎng)。

    Cho等人[31]研究Nd:YAG激光-TIG復(fù)合焊等離子體時(shí),忽略激光能量損失,得到復(fù)合焊等離子體的溫度場(chǎng)、電場(chǎng)和流場(chǎng)分布,分析表明陽(yáng)極激光輻射區(qū)域附近等離子體的局域溫度隨激光功率的增大而增大,熱流密度高度集中,流動(dòng)速度增大。

    Cho和Na[32]研究CO2激光-TIG復(fù)合焊等離子體時(shí),保護(hù)氣體為Ar,考慮激光能量的逆韌致輻射吸收,結(jié)果表明整個(gè)電弧軸線區(qū)域因等離子體對(duì)激光能量的大量吸收,溫度很高,電弧弧根部分因金屬蒸汽的影響而收縮。

    圖2 激光-TIG同軸復(fù)合焊模型Fig.2Hybrid model of laser-coaxial arc welding

    此模型針對(duì)靜止時(shí)的激光-TIG同軸復(fù)合焊提出,是目前基于激光-電弧復(fù)合等離子體模擬研究所建立的僅有的數(shù)學(xué)模型,Cho和Na也是在復(fù)合等離子體模擬研究領(lǐng)域的杰出學(xué)者代表。而實(shí)際復(fù)合焊接過(guò)程多采用高焊速的激光-電弧旁軸復(fù)合,此時(shí)電弧形態(tài)及電弧熱流密度分布均會(huì)發(fā)生變化。因此,研究高焊速下激光-電弧復(fù)合等離子體機(jī)制,還需進(jìn)一步改進(jìn)和完善模型。

    4 結(jié)論

    目前,激光-電弧復(fù)合熱源模型的研究應(yīng)用已引起眾多學(xué)者的廣泛關(guān)注,其研究體系也日益完善。相對(duì)而言,對(duì)于復(fù)合等離子體復(fù)雜的物理機(jī)制的研究,僅有少數(shù)研究者涉足,且受試驗(yàn)條件和數(shù)學(xué)模型的限制及材料的高溫?zé)嵛锢硇阅軘?shù)據(jù)的匱乏,人們對(duì)電弧物理行為的認(rèn)識(shí)尚不完善。

    激光-電弧復(fù)合等離子體與熔池是不可分割的統(tǒng)一體,兩者存在雙向耦合作用:復(fù)合等離子體作用于焊件表面,導(dǎo)致熔化金屬自由表面的變形;熔池表面變形不僅會(huì)改變?nèi)鄢貎?nèi)熱量及質(zhì)量傳輸,還會(huì)影響復(fù)合等離子體形狀及焊件表面的熱流密度和電流密度分布。將激光-電弧復(fù)合等離子體與熔池統(tǒng)一進(jìn)行分析,避免交互界面條件的假定,能更精確地揭示復(fù)合等離子體與熔池的交互作用機(jī)制,更能反映復(fù)合等離子體和熔池?cái)?shù)值模擬的實(shí)際情況。

    為實(shí)現(xiàn)激光-電弧復(fù)合等離子體與熔池?cái)?shù)值模擬的統(tǒng)一,根據(jù)磁流體動(dòng)力學(xué)理論構(gòu)建激光-電弧復(fù)合等離子體和熔池的三維數(shù)學(xué)模型,通過(guò)應(yīng)用相關(guān)有限元軟件,對(duì)復(fù)合等離子體與熔池的耦合進(jìn)行迭代分析,以此更新兩者間因界面的動(dòng)態(tài)變化而改變的計(jì)算區(qū)域,但考慮到金屬熔化的相變過(guò)程,故采用特定方法解決焊件的熔化區(qū)與非熔化區(qū)的移動(dòng)邊界問(wèn)題。

    研究激光-電弧復(fù)合區(qū)溫度場(chǎng)分布、電流密度分布及等離子體流場(chǎng)分布,討論分析等離子體流動(dòng)機(jī)理過(guò)程,由于焊接過(guò)程中復(fù)雜的多物理現(xiàn)象,其過(guò)程難以用數(shù)學(xué)方法準(zhǔn)確描述,尚待研究。

    現(xiàn)階段,焊接電弧與熔池交互作用的統(tǒng)一模型已經(jīng)建立,但激光-電弧復(fù)合等離子體與熔池的統(tǒng)一模型的建立還尚未報(bào)導(dǎo),為簡(jiǎn)化分析過(guò)程,通常將兩者作為獨(dú)立的研究對(duì)象,并假定適當(dāng)?shù)倪吔鐥l件,其將成為未來(lái)激光-電弧復(fù)合焊接數(shù)值模擬的研究動(dòng)態(tài)及發(fā)展趨勢(shì)。

    [1]Steen W M.Arc augmented laser welding[J].Matal Construction,1979,11(7):332–333.

    [2]吳圣川,劉建華.鋁合金激光-電弧復(fù)合焊的有限元數(shù)值模擬[J].航空制造技術(shù),2005(12):74-76.

    [3]Mahmoud Moradi,Majid Ghoreishi,Jan Frostevarg,et al. An investigation on stability of laser hybrid arc welding[J]. Optics and Lasers in Engineering,2013(51)481-487.

    [4]肖榮詩(shī),吳世凱.激光-電弧復(fù)合焊接的研究進(jìn)展[J].中國(guó)激光,2008,35(11):1680-1685.

    [5]夏源,宋永倫,胡坤平,等.激光-TIG復(fù)合焊接熱源機(jī)理研究現(xiàn)狀與進(jìn)展[J].焊接,2008(12):21-24.

    [6]Steen W M.Arc augmented laser processing materials[J]. Joural of Applied Physics,1980,51(11):5636-4641.

    [7]高明,曾曉雁,嚴(yán)軍,等.激光-電弧復(fù)合焊接的熱源相互作用[J].激光技術(shù),2007,31(5):465-468.

    [8]MoriakiO,YukioS,AkihideetY,et al.Development of laserarc hybrid welding[J].NKK Technical Review,2002(86): 8-12.

    [9]陳彥賓,陳杰,李俐群,等.激光與電弧相互作用時(shí)的電弧形態(tài)及焊縫特征[J].焊接學(xué)報(bào),2003,24(1):55-60.

    [10]張紹彬.激光-電弧雙熱源焊接新工藝的研究[D].天津:天津大學(xué),1989.

    [11]Albright C E,Eastman J,Lempert W.Low-power lasers assist arc welding[J].Welding Journal,2001,80(4):55-59.

    [12]Matsunawa A,Ohnawa T.Beam-plume interaction in plused YAG laser processing[A].International Congress on Applications of Laser&Electro-Optics’84.Boston,U.S.A,1984:34-36.

    [13]Tix C,Gratzke U,Simon G.Absorption of the laser beam by the plasma in deep laser beam welding of metals[J].J.Appl. Phys.,1995,78(11):6448-6453.

    [14]Hu B,Ouden G D.Laser induced stabilisation of the welding arc[J].Science and Technology of Welding and Joining,2005,10(1):76-80.

    [15]Bibik O B,Brodyaqin V N.Special features of inter-action of laser radiation with the electric welding arc in the combined laser-arc welding[J].Physics and Chemistry of Materials Treatment,1990,24(2):176-178.

    [16]劉黎明,遲鳴聲,宋剛.鎂合金激光-TIG復(fù)合熱源焊接熱源模型的建立及其數(shù)值模擬[J].機(jī)械工程學(xué)報(bào),2006,42(2):82-85.

    [17]李迅波,吳剛,曾志.鎂合金/不銹鋼激光-TIG復(fù)合焊的數(shù)值分析[J].焊接學(xué)報(bào),2011,32(8):49-52.

    [18]Xu G X,ChuanSongWu,Qin G L,et al.Adaptive volumetric heat source models for laser beam and laser+pulsed GMAW hybrid welding processes[J].Int J Adv Manuf Technol,2011(57):245-255.

    [19]Goldak J A.Computational Welding Mechanics,Springer,NY,2005.

    [20]Olivier Desmaison,Michel Bellet,Gildas Guillemot.A level set approach for the simulation of the multipass hybrid laser/ GMA welding process[J].Computational Materials Science 2014(91):240-250.

    [21]Fanrong Kong,Junjie Ma.Numerical and experimental study of thermally induced residual stress in the hybrid laser-GMA welding process[J].Journal of Materials Processing Technology,2011,211:1102–1111.

    [22]吳甦,趙海燕,王煜.高能束焊接數(shù)值模擬中的新型熱源模型[J].焊接學(xué)報(bào),2004,25(1):91-94.

    [23]Liang G L,Zhou G,Yuan S Q.Study on hybrid heat source overlap welding of magnesium alloy AZ31B[J].Materials Science and Engineering A,2009(499):93-96.

    [24]Piekarska W,Kubiak M.Modeling of thermal phenomena in single laser beam and laser-arc Hybrid welding processes using projection method[J].Applied Mathematical Modelling,2013(37):2051-2062.

    [25]Piekarska W,Kubiak M.Three-dimensional model for numerical analysis of thermal phenomena in laser-arc hybrid welding process[J].International Journal of Heat and Mass Transfer,2011,(54)4966-4974.

    [26]Guoxiang X U,Chuansong W U,Xuezhou MA.Numerical AnalysisofWeldingResidualStressandDistortioninLaser+ GMAW Hybrid Welding of Aluminum Alloy T-Joint[J]. Acta Metall.Sin.(Engl.Lett.),2013,26(3):352-360.

    [27]Amir Hossein Faraji,Massoud Goodarzi.Numerical modeling of heat transfer and fluid flow in hybrid laser-TIG welding of aluminum alloy AA6082[J].Int J Adv Manuf Technol,2015(77):2067-2082.

    [28]Won-Ik Cho,Suck-Joo Na,Min-Hyun Cho,et al.Numerical study of alloying element distribution in CO2laser-GMA hybrid welding[J].Computational Materials Science,2010(49):792-800.

    [29]過(guò)增元,趙文華.電弧和熱等離子體[M].北京:科學(xué)出版社,1986.

    [30]胥國(guó)祥,胡慶賢.激光+GMAW復(fù)合焊工藝及數(shù)值模擬[M].江蘇:江蘇大學(xué)出版社,2013.

    [31]Cho Y T,Cho W I,Na S J.Numerical analysis of hybrid plasma generated by Nd:YAG laser and gas tungsten arc[J]. Optics&Laser Technology,2011(43):711-720.

    [32]Young Tae Cho,Suck Joo Na.Numerical Analysis of Plasma in CO2Laser and Arc Hybrid Welding[J].International Journal of Precision engineering and manufacturing,2015,16(4):787-795.

    Numerical analysis of thermal model and plasma in laser-arc hybrid welding

    DING Lingfang,ZHANG Yingqiao,LIU Shengyao
    (Department of Materials Science and Engineering,North University of China,Taiyuan 030051,China)

    As a new type of efficient composite welding technology,laser-arc hybrid welding has attracted attention of scholars at home and abroad widely.In this paper,the coupling mechanism of heat source and the required heat source model in numerical simulation of laser-arc hybrid welding are outlined.And then,the application of coupling heat source and the hybrid plasma are detailed.And it points out that the establishment of unified mathematical model in laser-arc hybrid welding would be the future research trend.

    laser-arc hybrid;heat source model;plasma;numerical simulation

    TG456.7

    C

    1001-2303(2016)03-0137-06

    10.7512/j.issn.1001-2303.2016.03.29

    2015-07-06;

    2015-08-18

    山西省自然科學(xué)基金資助項(xiàng)目(2012011021-1)

    丁玲芳(1991—),河南濮陽(yáng)人,碩士,主要從事焊接電弧物理及材料焊接性方面的研究工作。

    猜你喜歡
    熔池熱源電弧
    故障電弧探測(cè)器與故障電弧保護(hù)裝置在工程中的應(yīng)用分析
    橫流熱源塔換熱性能研究
    煤氣與熱力(2021年3期)2021-06-09 06:16:20
    2219鋁合金激光電弧復(fù)合焊接及其溫度場(chǎng)的模擬
    航空電氣系統(tǒng)中故障電弧的分析
    電子制作(2017年22期)2017-02-02 07:10:15
    基于啟發(fā)式動(dòng)態(tài)規(guī)劃的冷熱源優(yōu)化控制
    電子制作(2017年19期)2017-02-02 07:08:31
    電弧焊熔池表征與熔透狀態(tài)映射研究
    焊接(2016年10期)2016-02-27 13:05:32
    一種基于頻域變換的熔池運(yùn)動(dòng)模糊圖像恢復(fù)算法
    焊接(2015年11期)2015-07-18 11:12:46
    MIG焊接熔池形成與凝固過(guò)程數(shù)值模擬
    焊接(2015年9期)2015-07-18 11:03:51
    中部槽激光-MAG復(fù)合熱源打底焊焊接工藝研究
    焊接(2015年8期)2015-07-18 10:59:13
    多類型熱源點(diǎn)共存下的區(qū)域熱力統(tǒng)籌供應(yīng)探討
    亚洲av美国av| 亚洲国产色片| 成人亚洲精品av一区二区| 麻豆一二三区av精品| 国产精品久久久久久久电影 | 村上凉子中文字幕在线| 婷婷六月久久综合丁香| 在线视频色国产色| 免费无遮挡裸体视频| 精品久久久久久久人妻蜜臀av| 淫妇啪啪啪对白视频| 91老司机精品| 久久热在线av| 亚洲国产中文字幕在线视频| 精品久久久久久久久久久久久| 国产aⅴ精品一区二区三区波| 一级毛片女人18水好多| 欧美高清成人免费视频www| 国内精品一区二区在线观看| 国产精品一区二区精品视频观看| 九九在线视频观看精品| 悠悠久久av| 国产精品影院久久| 国内揄拍国产精品人妻在线| 欧美成人性av电影在线观看| 99久久精品国产亚洲精品| 精品国内亚洲2022精品成人| 久久久久久九九精品二区国产| 国产私拍福利视频在线观看| 久久性视频一级片| 高潮久久久久久久久久久不卡| 日本成人三级电影网站| 欧美黄色片欧美黄色片| 男人和女人高潮做爰伦理| av黄色大香蕉| 麻豆成人av在线观看| av欧美777| 最近最新中文字幕大全免费视频| 法律面前人人平等表现在哪些方面| 国产免费男女视频| 中文字幕人成人乱码亚洲影| 国产精品 国内视频| 99久国产av精品| 99久久99久久久精品蜜桃| 久久久久久九九精品二区国产| 长腿黑丝高跟| 国产成人av教育| 日日夜夜操网爽| 亚洲一区二区三区色噜噜| 午夜免费观看网址| 麻豆av在线久日| 欧美中文综合在线视频| 在线a可以看的网站| 国产主播在线观看一区二区| 18禁裸乳无遮挡免费网站照片| 亚洲国产中文字幕在线视频| 欧美不卡视频在线免费观看| 丰满的人妻完整版| 久久久色成人| 日韩欧美在线乱码| 午夜免费成人在线视频| 日日干狠狠操夜夜爽| 久久草成人影院| 欧美绝顶高潮抽搐喷水| 草草在线视频免费看| 男女视频在线观看网站免费| 久久热在线av| 欧美+亚洲+日韩+国产| 成人特级av手机在线观看| bbb黄色大片| 国产1区2区3区精品| 男人舔女人的私密视频| 精品国产乱码久久久久久男人| 在线看三级毛片| 亚洲无线观看免费| 国产又色又爽无遮挡免费看| 欧美成人一区二区免费高清观看 | 岛国视频午夜一区免费看| 丝袜人妻中文字幕| 黑人欧美特级aaaaaa片| 亚洲av成人不卡在线观看播放网| 午夜亚洲福利在线播放| 色av中文字幕| 亚洲国产高清在线一区二区三| 亚洲,欧美精品.| 禁无遮挡网站| 韩国av一区二区三区四区| 免费看a级黄色片| 91av网站免费观看| 校园春色视频在线观看| av片东京热男人的天堂| 精品一区二区三区视频在线 | 99久久99久久久精品蜜桃| 男女午夜视频在线观看| 精品免费久久久久久久清纯| 免费在线观看成人毛片| 亚洲欧美日韩东京热| 国产一区二区在线观看日韩 | 成人av在线播放网站| 天堂av国产一区二区熟女人妻| 欧美+亚洲+日韩+国产| 国产精品九九99| 国产成人一区二区三区免费视频网站| www日本在线高清视频| 亚洲国产看品久久| 嫁个100分男人电影在线观看| 淫妇啪啪啪对白视频| ponron亚洲| 淫秽高清视频在线观看| 精品99又大又爽又粗少妇毛片 | 久99久视频精品免费| 一区二区三区国产精品乱码| 一级a爱片免费观看的视频| 99国产精品一区二区蜜桃av| 亚洲第一电影网av| 好男人电影高清在线观看| 嫩草影院精品99| 首页视频小说图片口味搜索| 精品不卡国产一区二区三区| 久久久久久国产a免费观看| 两人在一起打扑克的视频| 久久久久国产精品人妻aⅴ院| 亚洲片人在线观看| 嫁个100分男人电影在线观看| 人妻丰满熟妇av一区二区三区| 91字幕亚洲| 黄色视频,在线免费观看| 欧美一级毛片孕妇| 一区二区三区激情视频| 最近视频中文字幕2019在线8| 免费在线观看成人毛片| 一个人看的www免费观看视频| 18禁国产床啪视频网站| 亚洲人成伊人成综合网2020| h日本视频在线播放| 宅男免费午夜| 久久久久性生活片| 成年版毛片免费区| 欧美黑人欧美精品刺激| 首页视频小说图片口味搜索| 99视频精品全部免费 在线 | 日韩精品中文字幕看吧| 久久久水蜜桃国产精品网| 久久久久性生活片| 亚洲人成网站在线播放欧美日韩| 欧美性猛交╳xxx乱大交人| 欧美一级毛片孕妇| 一区二区三区高清视频在线| 亚洲aⅴ乱码一区二区在线播放| 高清在线国产一区| 嫩草影视91久久| 男女之事视频高清在线观看| 怎么达到女性高潮| 啦啦啦免费观看视频1| 欧美xxxx黑人xx丫x性爽| 精品国产乱子伦一区二区三区| 91字幕亚洲| 两个人视频免费观看高清| 高清在线国产一区| 嫩草影视91久久| 国产成人精品无人区| 国产aⅴ精品一区二区三区波| 国内少妇人妻偷人精品xxx网站 | 午夜福利高清视频| 男女视频在线观看网站免费| 亚洲中文日韩欧美视频| 久久精品国产清高在天天线| xxxwww97欧美| 在线观看66精品国产| 国产欧美日韩精品亚洲av| 日本撒尿小便嘘嘘汇集6| 国产精品国产高清国产av| 亚洲,欧美精品.| 久久久久性生活片| 国产精品久久久久久久电影 | 日本黄色片子视频| 欧美三级亚洲精品| 亚洲精品中文字幕一二三四区| 级片在线观看| 亚洲欧美日韩无卡精品| 美女cb高潮喷水在线观看 | 亚洲七黄色美女视频| 亚洲欧洲精品一区二区精品久久久| ponron亚洲| 成人午夜高清在线视频| 99久久精品一区二区三区| 亚洲国产精品999在线| 黄片大片在线免费观看| 国产精品乱码一区二三区的特点| 一二三四社区在线视频社区8| 欧美3d第一页| 国产亚洲精品久久久久久毛片| 九九热线精品视视频播放| 校园春色视频在线观看| 日本五十路高清| 欧美xxxx黑人xx丫x性爽| 成人精品一区二区免费| 国产私拍福利视频在线观看| 成人三级做爰电影| 久久欧美精品欧美久久欧美| 国产亚洲av嫩草精品影院| 久久精品亚洲精品国产色婷小说| 国产精品九九99| 精品99又大又爽又粗少妇毛片 | 亚洲欧美精品综合久久99| 欧美三级亚洲精品| 男女视频在线观看网站免费| 国产午夜精品论理片| 亚洲,欧美精品.| 最新在线观看一区二区三区| 三级国产精品欧美在线观看 | 国产精品久久久久久久电影 | 老司机午夜福利在线观看视频| 亚洲电影在线观看av| 俄罗斯特黄特色一大片| 国产av麻豆久久久久久久| 国内揄拍国产精品人妻在线| 午夜精品久久久久久毛片777| 亚洲天堂国产精品一区在线| 国产av在哪里看| 一个人免费在线观看的高清视频| 亚洲第一欧美日韩一区二区三区| 午夜福利高清视频| 亚洲成av人片免费观看| 亚洲va日本ⅴa欧美va伊人久久| 十八禁网站免费在线| 91麻豆av在线| 日本免费一区二区三区高清不卡| 黄片大片在线免费观看| 欧美日韩精品网址| 怎么达到女性高潮| 欧美午夜高清在线| 校园春色视频在线观看| 午夜福利18| 精品欧美国产一区二区三| 在线十欧美十亚洲十日本专区| 九色成人免费人妻av| 久久久国产成人精品二区| 国产一区二区三区视频了| 最近最新中文字幕大全免费视频| 黑人操中国人逼视频| 国产私拍福利视频在线观看| 99热这里只有精品一区 | 亚洲五月婷婷丁香| 老司机在亚洲福利影院| 欧美日韩乱码在线| 少妇的逼水好多| 国产欧美日韩一区二区三| 国产视频内射| 一进一出抽搐gif免费好疼| a级毛片a级免费在线| 很黄的视频免费| 日日摸夜夜添夜夜添小说| 两性午夜刺激爽爽歪歪视频在线观看| 小说图片视频综合网站| 色噜噜av男人的天堂激情| 男人的好看免费观看在线视频| 丁香六月欧美| 国产野战对白在线观看| 18禁美女被吸乳视频| 国产麻豆成人av免费视频| 国产av麻豆久久久久久久| 精品久久久久久,| 久久久国产成人免费| 国产精品1区2区在线观看.| 国产精品精品国产色婷婷| 久久久久久久久免费视频了| 99久久精品热视频| 亚洲一区高清亚洲精品| 久久久久亚洲av毛片大全| 午夜福利视频1000在线观看| 精品不卡国产一区二区三区| 日本 欧美在线| 欧美乱码精品一区二区三区| 国产欧美日韩精品亚洲av| 色尼玛亚洲综合影院| 亚洲欧美日韩无卡精品| 成年免费大片在线观看| 欧美+亚洲+日韩+国产| 毛片女人毛片| 亚洲成人免费电影在线观看| 亚洲国产欧美一区二区综合| 欧美乱色亚洲激情| 听说在线观看完整版免费高清| 国内精品久久久久精免费| 国内揄拍国产精品人妻在线| www国产在线视频色| 日韩免费av在线播放| 十八禁人妻一区二区| 亚洲欧美日韩东京热| 精品久久久久久久毛片微露脸| 夜夜躁狠狠躁天天躁| 久久精品aⅴ一区二区三区四区| 性欧美人与动物交配| 老熟妇仑乱视频hdxx| 日韩高清综合在线| 国产精品1区2区在线观看.| 男女之事视频高清在线观看| 国产免费男女视频| 国产精品自产拍在线观看55亚洲| 欧美高清成人免费视频www| 免费高清视频大片| 欧美一级a爱片免费观看看| 国产欧美日韩一区二区三| 老司机福利观看| 人人妻人人澡欧美一区二区| 日韩有码中文字幕| 悠悠久久av| 日本一二三区视频观看| 精品不卡国产一区二区三区| 18美女黄网站色大片免费观看| 超碰成人久久| 一个人免费在线观看电影 | 51午夜福利影视在线观看| 中文资源天堂在线| 亚洲性夜色夜夜综合| 国产伦精品一区二区三区四那| 免费电影在线观看免费观看| 亚洲欧洲精品一区二区精品久久久| 国产精品香港三级国产av潘金莲| 成熟少妇高潮喷水视频| 亚洲国产看品久久| 99热这里只有精品一区 | 免费在线观看成人毛片| 老熟妇仑乱视频hdxx| 香蕉国产在线看| 女同久久另类99精品国产91| 欧美成人一区二区免费高清观看 | 国产真实乱freesex| 熟女电影av网| 久久久久九九精品影院| 日本黄色视频三级网站网址| 男人的好看免费观看在线视频| 亚洲专区国产一区二区| 中亚洲国语对白在线视频| 国产精品乱码一区二三区的特点| 两个人视频免费观看高清| 欧美绝顶高潮抽搐喷水| 亚洲av第一区精品v没综合| netflix在线观看网站| 国产精品日韩av在线免费观看| 天堂网av新在线| 午夜免费观看网址| 两个人看的免费小视频| 国产欧美日韩一区二区三| 国产久久久一区二区三区| 欧美在线黄色| 欧美另类亚洲清纯唯美| 精品国产三级普通话版| 国产精品永久免费网站| 热99在线观看视频| 亚洲av电影不卡..在线观看| 国产精品香港三级国产av潘金莲| 欧洲精品卡2卡3卡4卡5卡区| 色综合亚洲欧美另类图片| www.精华液| 国产又色又爽无遮挡免费看| 国产久久久一区二区三区| 亚洲片人在线观看| 级片在线观看| 五月玫瑰六月丁香| 欧美最黄视频在线播放免费| 成年女人永久免费观看视频| 少妇的逼水好多| 九九热线精品视视频播放| 免费一级毛片在线播放高清视频| 久久精品综合一区二区三区| 亚洲成av人片在线播放无| 噜噜噜噜噜久久久久久91| 麻豆一二三区av精品| 国产亚洲欧美在线一区二区| 在线国产一区二区在线| 久久久久久久午夜电影| netflix在线观看网站| 天天添夜夜摸| 免费看光身美女| 免费观看人在逋| 熟妇人妻久久中文字幕3abv| а√天堂www在线а√下载| 久久精品亚洲精品国产色婷小说| 岛国在线免费视频观看| 免费av不卡在线播放| 中亚洲国语对白在线视频| av福利片在线观看| 桃红色精品国产亚洲av| 三级国产精品欧美在线观看 | 亚洲成人中文字幕在线播放| 欧洲精品卡2卡3卡4卡5卡区| 美女黄网站色视频| 欧美最黄视频在线播放免费| 露出奶头的视频| 国产亚洲精品av在线| 最近视频中文字幕2019在线8| 搡老熟女国产l中国老女人| 黄色视频,在线免费观看| 一本综合久久免费| 欧美在线一区亚洲| 手机成人av网站| 精品久久久久久久末码| 国产99白浆流出| av女优亚洲男人天堂 | 99久国产av精品| 国产精品香港三级国产av潘金莲| 日韩欧美国产在线观看| 亚洲精品美女久久av网站| 看片在线看免费视频| a级毛片a级免费在线| 亚洲一区二区三区不卡视频| 日本免费a在线| 国产精品一区二区三区四区免费观看 | 国内久久婷婷六月综合欲色啪| 欧美午夜高清在线| 国产综合懂色| 国产成人福利小说| 亚洲午夜理论影院| 性色av乱码一区二区三区2| 一个人看视频在线观看www免费 | 精品久久久久久久久久久久久| 亚洲,欧美精品.| 免费在线观看日本一区| 90打野战视频偷拍视频| 老司机深夜福利视频在线观看| 欧美成人性av电影在线观看| 亚洲欧美日韩东京热| 男女视频在线观看网站免费| 午夜福利在线观看免费完整高清在 | 久久国产精品影院| 免费搜索国产男女视频| 日韩国内少妇激情av| 九九热线精品视视频播放| 丁香六月欧美| 亚洲av成人av| 高潮久久久久久久久久久不卡| 一区福利在线观看| 97超级碰碰碰精品色视频在线观看| 国产午夜福利久久久久久| 麻豆av在线久日| 国产精品一区二区三区四区免费观看 | 亚洲国产日韩欧美精品在线观看 | 亚洲精品在线观看二区| 一本久久中文字幕| 欧美乱色亚洲激情| 少妇裸体淫交视频免费看高清| 日日夜夜操网爽| 久久久久性生活片| 麻豆成人av在线观看| 日本免费a在线| 中文字幕高清在线视频| 久久精品影院6| 观看免费一级毛片| 国产高清激情床上av| 女人被狂操c到高潮| 日韩免费av在线播放| tocl精华| 亚洲 欧美 日韩 在线 免费| 国产精品久久久久久人妻精品电影| 日韩国内少妇激情av| 欧美性猛交╳xxx乱大交人| 窝窝影院91人妻| 男女床上黄色一级片免费看| 国产av一区在线观看免费| 日韩 欧美 亚洲 中文字幕| 成人鲁丝片一二三区免费| 日韩欧美国产一区二区入口| 成人av在线播放网站| 黑人巨大精品欧美一区二区mp4| tocl精华| av中文乱码字幕在线| 亚洲男人的天堂狠狠| 黄片小视频在线播放| 一级毛片高清免费大全| а√天堂www在线а√下载| 黑人巨大精品欧美一区二区mp4| 久9热在线精品视频| 午夜福利在线在线| 精品日产1卡2卡| 岛国视频午夜一区免费看| 国产成人精品久久二区二区免费| 欧美成人性av电影在线观看| 国产一级毛片七仙女欲春2| 无人区码免费观看不卡| 欧美在线黄色| 天堂网av新在线| 国产精品精品国产色婷婷| 国产黄片美女视频| 好男人在线观看高清免费视频| 男女午夜视频在线观看| 国产高清videossex| 国产精品九九99| 亚洲欧美日韩东京热| 一本久久中文字幕| av中文乱码字幕在线| 日本三级黄在线观看| www.999成人在线观看| 韩国av一区二区三区四区| 天天躁日日操中文字幕| 久久久久久人人人人人| 老司机福利观看| 天堂动漫精品| 一个人看视频在线观看www免费 | 成人鲁丝片一二三区免费| 色综合亚洲欧美另类图片| 日日摸夜夜添夜夜添小说| а√天堂www在线а√下载| 成人特级黄色片久久久久久久| 精品无人区乱码1区二区| 热99在线观看视频| 天堂av国产一区二区熟女人妻| 国产成人精品无人区| 麻豆久久精品国产亚洲av| 国产单亲对白刺激| 国产成年人精品一区二区| 香蕉久久夜色| 男人的好看免费观看在线视频| 日本五十路高清| 久久久久精品国产欧美久久久| 视频区欧美日本亚洲| 老熟妇乱子伦视频在线观看| 操出白浆在线播放| 精品国产三级普通话版| 99国产精品99久久久久| 91久久精品国产一区二区成人 | 少妇人妻一区二区三区视频| 国产在线精品亚洲第一网站| 观看美女的网站| 国产精品98久久久久久宅男小说| 亚洲精品久久国产高清桃花| 九九在线视频观看精品| 香蕉国产在线看| 欧美又色又爽又黄视频| 久久草成人影院| 日韩欧美一区二区三区在线观看| 夜夜躁狠狠躁天天躁| 国产一区二区三区在线臀色熟女| 天天一区二区日本电影三级| 亚洲五月天丁香| 色哟哟哟哟哟哟| 在线免费观看不下载黄p国产 | 又大又爽又粗| 两个人的视频大全免费| 午夜激情欧美在线| 日韩成人在线观看一区二区三区| 熟女少妇亚洲综合色aaa.| 午夜免费观看网址| 国产激情欧美一区二区| 偷拍熟女少妇极品色| 国语自产精品视频在线第100页| 啪啪无遮挡十八禁网站| 亚洲av五月六月丁香网| 久久久久国内视频| 国产99白浆流出| 非洲黑人性xxxx精品又粗又长| 最好的美女福利视频网| 色综合欧美亚洲国产小说| 色噜噜av男人的天堂激情| 一卡2卡三卡四卡精品乱码亚洲| 热99在线观看视频| 国产综合懂色| 床上黄色一级片| 一级毛片高清免费大全| 俄罗斯特黄特色一大片| 日韩精品中文字幕看吧| 18禁裸乳无遮挡免费网站照片| 999精品在线视频| 国产成人一区二区三区免费视频网站| 国产男靠女视频免费网站| 中出人妻视频一区二区| 他把我摸到了高潮在线观看| 91老司机精品| 国产三级黄色录像| 中文资源天堂在线| 一个人免费在线观看电影 | 男女之事视频高清在线观看| h日本视频在线播放| 久久性视频一级片| 精品电影一区二区在线| 97超级碰碰碰精品色视频在线观看| 五月伊人婷婷丁香| 亚洲欧洲精品一区二区精品久久久| 亚洲av成人av| 亚洲成av人片在线播放无| 69av精品久久久久久| 国产一级毛片七仙女欲春2| 国产视频内射| 精品国内亚洲2022精品成人| 亚洲狠狠婷婷综合久久图片| 少妇熟女aⅴ在线视频| 精品久久久久久久末码| 亚洲电影在线观看av| 操出白浆在线播放| 97人妻精品一区二区三区麻豆| 国产 一区 欧美 日韩| 视频区欧美日本亚洲| 欧美不卡视频在线免费观看| 麻豆成人午夜福利视频| 亚洲成人精品中文字幕电影| 日本三级黄在线观看| 国产一区二区在线av高清观看| 久久久久免费精品人妻一区二区| 国模一区二区三区四区视频 | 久久久久九九精品影院| 91字幕亚洲| 特大巨黑吊av在线直播| 日韩欧美三级三区| 亚洲av成人一区二区三| 69av精品久久久久久| 老汉色∧v一级毛片| 久久精品国产99精品国产亚洲性色| 国产精品一及| 国产亚洲精品综合一区在线观看| 精品免费久久久久久久清纯| 久久香蕉精品热| 精品无人区乱码1区二区| 国内揄拍国产精品人妻在线| 怎么达到女性高潮|