• <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电影在线观看| 一二三四社区在线视频社区8| 黄色女人牲交| 亚洲全国av大片| 免费观看精品视频网站| 国产精品久久视频播放| 午夜福利,免费看| 精品一区二区三区视频在线观看免费 | 伦理电影免费视频| 很黄的视频免费| 久久久久国内视频| 久久中文字幕人妻熟女| 亚洲精品美女久久av网站| 国产精品久久久av美女十八| 欧美日本亚洲视频在线播放| 国产精品亚洲一级av第二区| 看免费av毛片| 狂野欧美激情性xxxx| 超色免费av| 亚洲人成77777在线视频| 黄色女人牲交| 乱人伦中国视频| 在线观看66精品国产| 久久久久久久久中文| 在线播放国产精品三级| 久久性视频一级片| 黄色视频不卡| 久久久国产欧美日韩av| 亚洲 欧美 日韩 在线 免费| 午夜日韩欧美国产| 久久久久精品国产欧美久久久| 99热国产这里只有精品6| 欧美在线一区亚洲| 亚洲精品美女久久久久99蜜臀| 波多野结衣高清无吗| 韩国av一区二区三区四区| 制服人妻中文乱码| 欧洲精品卡2卡3卡4卡5卡区| 精品无人区乱码1区二区| 成人手机av| 一个人观看的视频www高清免费观看 | 老司机靠b影院| 国产精品综合久久久久久久免费 | 两人在一起打扑克的视频| 男人的好看免费观看在线视频 | 国产精品99久久99久久久不卡| 免费搜索国产男女视频| av视频免费观看在线观看| 一区二区三区激情视频| 亚洲五月色婷婷综合| 少妇 在线观看| 亚洲欧洲精品一区二区精品久久久| 午夜老司机福利片| 高清毛片免费观看视频网站 | a在线观看视频网站| 日韩人妻精品一区2区三区| 极品教师在线免费播放| 人人妻人人爽人人添夜夜欢视频| 精品久久久久久,| 天天躁夜夜躁狠狠躁躁| 色哟哟哟哟哟哟| 老熟妇仑乱视频hdxx| 在线观看www视频免费| 亚洲欧美日韩另类电影网站| 欧美最黄视频在线播放免费 | 久久人人爽av亚洲精品天堂| 在线播放国产精品三级| 搡老乐熟女国产| 欧美日韩福利视频一区二区| 一级片免费观看大全| 久久久国产精品麻豆| 80岁老熟妇乱子伦牲交| 国产精品久久久人人做人人爽| 成人免费观看视频高清| 18禁国产床啪视频网站| av天堂在线播放| 亚洲三区欧美一区| 亚洲午夜理论影院| 国产日韩一区二区三区精品不卡| www国产在线视频色| bbb黄色大片| 中文欧美无线码| 国产成人精品无人区| 午夜福利在线免费观看网站| 深夜精品福利| 男女做爰动态图高潮gif福利片 | 成人18禁高潮啪啪吃奶动态图| 中出人妻视频一区二区| 女警被强在线播放| 后天国语完整版免费观看| 叶爱在线成人免费视频播放| 嫩草影院精品99| 精品国产乱码久久久久久男人| 啦啦啦 在线观看视频| 精品一品国产午夜福利视频| 天天影视国产精品| 男女午夜视频在线观看| 欧美av亚洲av综合av国产av| 国产主播在线观看一区二区| 亚洲国产精品999在线| 亚洲欧美日韩无卡精品| 搡老岳熟女国产| 国产免费av片在线观看野外av| 欧美在线黄色| 一进一出抽搐gif免费好疼 | 国产av一区二区精品久久| 村上凉子中文字幕在线| 天堂中文最新版在线下载| www.999成人在线观看| 国产欧美日韩精品亚洲av| 99精品欧美一区二区三区四区| 欧美黑人精品巨大| 丰满的人妻完整版| 国产精品二区激情视频| 国产精品九九99| 久久久久久大精品| 午夜亚洲福利在线播放| 免费看十八禁软件| 99riav亚洲国产免费| 香蕉国产在线看| 午夜老司机福利片| 国产极品粉嫩免费观看在线| 欧美黄色淫秽网站| 在线播放国产精品三级| 一区二区三区国产精品乱码| 国产精品国产高清国产av| 国产精品99久久99久久久不卡| 曰老女人黄片| 欧美乱色亚洲激情| 亚洲专区国产一区二区| 妹子高潮喷水视频| 国产真人三级小视频在线观看| 国产高清视频在线播放一区| 成人特级黄色片久久久久久久| 一区二区三区激情视频| 亚洲精品国产区一区二| 悠悠久久av| 自线自在国产av| 日韩 欧美 亚洲 中文字幕| 午夜福利欧美成人| 国产亚洲精品第一综合不卡| 久久人人精品亚洲av| 亚洲伊人色综图| 波多野结衣高清无吗| 人人妻,人人澡人人爽秒播| 97碰自拍视频| www.自偷自拍.com| 天天躁夜夜躁狠狠躁躁| 好男人电影高清在线观看| 日本免费一区二区三区高清不卡 | 亚洲va日本ⅴa欧美va伊人久久| 久久精品亚洲熟妇少妇任你| 在线观看66精品国产| 自线自在国产av| 另类亚洲欧美激情| 午夜a级毛片| 久久精品aⅴ一区二区三区四区| 午夜福利欧美成人| 日本黄色日本黄色录像| 99国产精品一区二区蜜桃av| 欧美色视频一区免费| 午夜老司机福利片| 国产区一区二久久| 久久伊人香网站| 国产1区2区3区精品| 亚洲欧美一区二区三区久久| 久久 成人 亚洲| 欧美最黄视频在线播放免费 | 黑人巨大精品欧美一区二区mp4| 女人被躁到高潮嗷嗷叫费观| 男女做爰动态图高潮gif福利片 | 久久青草综合色| 真人做人爱边吃奶动态| 不卡av一区二区三区| 国产精品国产高清国产av| 亚洲国产欧美日韩在线播放| 亚洲国产精品sss在线观看 | 国产一区二区三区在线臀色熟女 | 亚洲欧美一区二区三区久久| 激情在线观看视频在线高清| 又紧又爽又黄一区二区| 久久人人97超碰香蕉20202| 日韩欧美一区视频在线观看| 久久久久久久午夜电影 | 高潮久久久久久久久久久不卡| 91大片在线观看| 人妻丰满熟妇av一区二区三区| 久久久久久久久中文| 99国产精品一区二区三区| 9色porny在线观看| 真人做人爱边吃奶动态| 两性午夜刺激爽爽歪歪视频在线观看 | 极品人妻少妇av视频| 熟女少妇亚洲综合色aaa.| 日韩大尺度精品在线看网址 | 一级毛片精品| 窝窝影院91人妻| 黑人操中国人逼视频| 欧美性长视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 妹子高潮喷水视频| 久久欧美精品欧美久久欧美| 黑人欧美特级aaaaaa片| 国产精品乱码一区二三区的特点 | 精品高清国产在线一区| 国产精品久久久久成人av| 黑人巨大精品欧美一区二区mp4| 免费在线观看影片大全网站| 国产亚洲精品第一综合不卡| 一区二区日韩欧美中文字幕| 十八禁网站免费在线| 俄罗斯特黄特色一大片| 电影成人av| 欧美日韩一级在线毛片| 久久久久国产一级毛片高清牌| 天堂影院成人在线观看| 免费一级毛片在线播放高清视频 | 久久久久久久久中文| 激情在线观看视频在线高清| 身体一侧抽搐| 久久国产乱子伦精品免费另类| 日本精品一区二区三区蜜桃| 日韩欧美免费精品| 日韩一卡2卡3卡4卡2021年| 搡老乐熟女国产| 国产精品国产高清国产av| 村上凉子中文字幕在线| 久久天堂一区二区三区四区| 久久久久久大精品| 青草久久国产| 成人影院久久| 欧美日韩av久久| 日韩精品免费视频一区二区三区| 欧美激情高清一区二区三区| 69精品国产乱码久久久| 国产激情久久老熟女| 91精品国产国语对白视频| 激情视频va一区二区三区| 午夜福利一区二区在线看| 色综合站精品国产| 国产主播在线观看一区二区| 久久国产精品男人的天堂亚洲| 免费在线观看完整版高清| 美女午夜性视频免费| 国产三级黄色录像| 亚洲人成77777在线视频| 69av精品久久久久久| 精品熟女少妇八av免费久了| 久久久国产成人精品二区 | 一级a爱片免费观看的视频| 成年人免费黄色播放视频| 1024视频免费在线观看| 色综合欧美亚洲国产小说| 国产伦一二天堂av在线观看| 亚洲成人久久性| 国产欧美日韩一区二区三| 黄频高清免费视频| 很黄的视频免费| 国产成人免费无遮挡视频| 亚洲色图av天堂| 亚洲自拍偷在线| 一个人免费在线观看的高清视频| 精品熟女少妇八av免费久了| 亚洲情色 制服丝袜| 国产精品99久久99久久久不卡| 老熟妇仑乱视频hdxx| 一区福利在线观看| 琪琪午夜伦伦电影理论片6080| 丝袜人妻中文字幕| 欧美日韩亚洲国产一区二区在线观看| 国产免费男女视频| 黄网站色视频无遮挡免费观看| 国产三级在线视频| 长腿黑丝高跟| 交换朋友夫妻互换小说| 日本一区二区免费在线视频| xxx96com| 日韩国内少妇激情av| 老熟妇仑乱视频hdxx| 日韩精品免费视频一区二区三区| 亚洲熟妇熟女久久| 欧美色视频一区免费| 亚洲精品国产一区二区精华液| 人人妻人人添人人爽欧美一区卜| 色老头精品视频在线观看| 波多野结衣高清无吗| 十八禁网站免费在线| 欧美在线一区亚洲| 国产精品亚洲av一区麻豆| 男男h啪啪无遮挡| 一进一出抽搐gif免费好疼 | 在线观看午夜福利视频| 精品久久久久久电影网| 久久影院123| 在线观看66精品国产| 免费av毛片视频| 亚洲五月婷婷丁香| 国产在线精品亚洲第一网站| 亚洲五月色婷婷综合| 亚洲精品久久成人aⅴ小说| 少妇裸体淫交视频免费看高清 | 琪琪午夜伦伦电影理论片6080| 国产精品成人在线| 免费久久久久久久精品成人欧美视频| 亚洲色图综合在线观看| 色老头精品视频在线观看| 99re在线观看精品视频| 日本a在线网址| 亚洲成a人片在线一区二区| 淫妇啪啪啪对白视频| 欧美日韩瑟瑟在线播放| 欧美av亚洲av综合av国产av| 一个人免费在线观看的高清视频| 高清毛片免费观看视频网站 | 久久性视频一级片| 亚洲人成电影免费在线| 99精品在免费线老司机午夜| 在线观看免费视频日本深夜| 亚洲激情在线av| 十八禁人妻一区二区| 亚洲精品国产一区二区精华液| 精品少妇一区二区三区视频日本电影| 80岁老熟妇乱子伦牲交| 国产黄a三级三级三级人| 桃红色精品国产亚洲av| xxx96com| 国产欧美日韩一区二区三区在线| 黄色视频不卡| 99久久国产精品久久久| 啦啦啦在线免费观看视频4| 中文亚洲av片在线观看爽| 久久狼人影院| 满18在线观看网站| 成人亚洲精品av一区二区 | 欧美日本中文国产一区发布| 亚洲欧美一区二区三区黑人| 桃色一区二区三区在线观看| 99国产精品99久久久久| 成人国产一区最新在线观看| 琪琪午夜伦伦电影理论片6080| 99热国产这里只有精品6| 日日摸夜夜添夜夜添小说| 午夜精品国产一区二区电影| 欧美日韩亚洲高清精品| 老熟妇仑乱视频hdxx| 日韩高清综合在线| 久久性视频一级片| 涩涩av久久男人的天堂| 久久国产精品影院| 成人av一区二区三区在线看| 欧美乱色亚洲激情| 老司机亚洲免费影院| 曰老女人黄片| 亚洲专区国产一区二区| 欧美亚洲日本最大视频资源| 欧美在线黄色| av超薄肉色丝袜交足视频| 久久久久久久久久久久大奶| 国产成人av教育| 制服人妻中文乱码| 久久精品亚洲熟妇少妇任你| 少妇粗大呻吟视频| 欧美老熟妇乱子伦牲交| 婷婷六月久久综合丁香| 久久亚洲真实| 婷婷丁香在线五月| 激情在线观看视频在线高清| 久久精品国产亚洲av香蕉五月| videosex国产| 欧美日韩福利视频一区二区| 看免费av毛片| e午夜精品久久久久久久| 男人操女人黄网站| xxx96com| 午夜两性在线视频| 91精品三级在线观看| 欧美色视频一区免费| 国产成人av激情在线播放| 亚洲五月色婷婷综合| 久久中文看片网| 国产黄色免费在线视频| 夜夜躁狠狠躁天天躁| 免费在线观看影片大全网站| 十八禁网站免费在线| 久99久视频精品免费| 真人做人爱边吃奶动态| 国产欧美日韩一区二区三区在线| 美女高潮到喷水免费观看| 18禁黄网站禁片午夜丰满| 丰满迷人的少妇在线观看| 日韩免费av在线播放| 五月开心婷婷网| 精品国产一区二区久久| 成人永久免费在线观看视频| 亚洲国产精品sss在线观看 | 国产成人av教育| 日韩国内少妇激情av| 午夜免费观看网址| 久久午夜亚洲精品久久| 正在播放国产对白刺激| 亚洲熟女毛片儿| 天堂动漫精品| 视频区图区小说| 香蕉久久夜色| 精品久久久久久,| 亚洲人成电影观看| 欧美人与性动交α欧美精品济南到| 欧美一级毛片孕妇| 不卡av一区二区三区| 亚洲成人精品中文字幕电影 | 丰满的人妻完整版| 99久久国产精品久久久| 老司机午夜福利在线观看视频| 啦啦啦在线免费观看视频4| 三上悠亚av全集在线观看| 国产一区二区三区在线臀色熟女 | 老汉色∧v一级毛片| 极品教师在线免费播放| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲午夜精品一区,二区,三区| 黄色女人牲交| 亚洲一区二区三区欧美精品| 高清在线国产一区| 在线观看66精品国产| 亚洲久久久国产精品| 亚洲欧美精品综合久久99| 97人妻天天添夜夜摸| 午夜福利一区二区在线看| 亚洲人成电影观看| 日日夜夜操网爽| 久久久久亚洲av毛片大全| 亚洲av日韩精品久久久久久密| 夜夜夜夜夜久久久久| 黑人巨大精品欧美一区二区mp4| 两个人免费观看高清视频| 国产色视频综合| 国产精品野战在线观看 | 后天国语完整版免费观看| 欧美最黄视频在线播放免费 | 成年人免费黄色播放视频| 日韩一卡2卡3卡4卡2021年| 色婷婷久久久亚洲欧美| 亚洲激情在线av| 在线观看日韩欧美| 久久精品亚洲精品国产色婷小说| 国产无遮挡羞羞视频在线观看| 免费高清视频大片| 中出人妻视频一区二区| 国产成人欧美| 精品国产乱码久久久久久男人| 国产高清videossex| 18禁观看日本| 人成视频在线观看免费观看| 怎么达到女性高潮| 国产真人三级小视频在线观看| 久久九九热精品免费| 欧美日韩国产mv在线观看视频| 日韩精品免费视频一区二区三区| 中文欧美无线码| 无限看片的www在线观看| 亚洲成人国产一区在线观看| 这个男人来自地球电影免费观看| 亚洲国产欧美网| 国产主播在线观看一区二区| 后天国语完整版免费观看| 777久久人妻少妇嫩草av网站| 19禁男女啪啪无遮挡网站| 亚洲一区高清亚洲精品| 亚洲欧美激情在线| 视频区欧美日本亚洲| 夫妻午夜视频| 村上凉子中文字幕在线| 超色免费av| 国产亚洲精品一区二区www| 国产免费av片在线观看野外av| 99精品久久久久人妻精品| 国产精品久久久久成人av| 亚洲av美国av| 女人被狂操c到高潮| 一级a爱视频在线免费观看| 五月开心婷婷网| 丰满饥渴人妻一区二区三| 国产成人精品在线电影| 美女扒开内裤让男人捅视频| 丝袜在线中文字幕| 亚洲欧美一区二区三区黑人| 91老司机精品| avwww免费| 日本wwww免费看| 女同久久另类99精品国产91| 午夜亚洲福利在线播放| 亚洲第一欧美日韩一区二区三区| 亚洲人成网站在线播放欧美日韩| 99精国产麻豆久久婷婷| 91大片在线观看| 亚洲国产精品sss在线观看 | 亚洲精品国产色婷婷电影| 色在线成人网| 欧美人与性动交α欧美精品济南到| www日本在线高清视频| 黄色视频不卡| 久久伊人香网站| 精品久久久久久,| 18禁黄网站禁片午夜丰满| 欧美成人午夜精品| av超薄肉色丝袜交足视频| 91九色精品人成在线观看| 亚洲熟妇中文字幕五十中出 | 国产精品影院久久| 男人的好看免费观看在线视频 | 久热爱精品视频在线9| av在线天堂中文字幕 | 欧美精品亚洲一区二区| 露出奶头的视频| 人人妻人人添人人爽欧美一区卜| 亚洲片人在线观看| 亚洲欧美精品综合久久99| 欧美日本亚洲视频在线播放| 日本一区二区免费在线视频| 亚洲精品中文字幕在线视频| 在线观看免费日韩欧美大片| 亚洲一区中文字幕在线| 久久中文字幕一级| 丝袜美足系列| 久久久水蜜桃国产精品网| 丰满迷人的少妇在线观看| 少妇被粗大的猛进出69影院| 色婷婷av一区二区三区视频| 久久久久久人人人人人| a在线观看视频网站| 女性被躁到高潮视频| 精品久久久久久久毛片微露脸| 日日爽夜夜爽网站| a级片在线免费高清观看视频| 精品欧美一区二区三区在线| 嫩草影视91久久| 亚洲色图 男人天堂 中文字幕| 一级a爱视频在线免费观看| 大型av网站在线播放| 亚洲va日本ⅴa欧美va伊人久久| 亚洲狠狠婷婷综合久久图片| 麻豆一二三区av精品| 亚洲精品美女久久久久99蜜臀| 久久午夜综合久久蜜桃| 国产蜜桃级精品一区二区三区| 两人在一起打扑克的视频| 国产野战对白在线观看| 精品久久久久久电影网| 国产欧美日韩一区二区三区在线| 黄片大片在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 国产精品av久久久久免费| 女人爽到高潮嗷嗷叫在线视频| 黄频高清免费视频| 高潮久久久久久久久久久不卡| 欧美黑人欧美精品刺激| 精品国产乱子伦一区二区三区| 99国产精品免费福利视频| 日日夜夜操网爽| 久久久国产精品麻豆| 91九色精品人成在线观看| 日韩欧美免费精品| 美女国产高潮福利片在线看| 久久久久久久精品吃奶| av免费在线观看网站| 久久九九热精品免费| 久久久久久人人人人人| 国产成人一区二区三区免费视频网站| 亚洲av成人不卡在线观看播放网| 两人在一起打扑克的视频| 巨乳人妻的诱惑在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产xxxxx性猛交| 夫妻午夜视频| 99热只有精品国产| 日日爽夜夜爽网站| 露出奶头的视频| 欧美性长视频在线观看| 亚洲一码二码三码区别大吗| 亚洲色图av天堂| 老鸭窝网址在线观看| 亚洲三区欧美一区| 99香蕉大伊视频| 国产亚洲精品第一综合不卡| 亚洲第一欧美日韩一区二区三区| 精品一区二区三卡| 久久伊人香网站| 亚洲男人天堂网一区| 人人妻人人添人人爽欧美一区卜| 黄色视频不卡| av电影中文网址| 一夜夜www| 午夜激情av网站| tocl精华| 日韩欧美国产一区二区入口| 两人在一起打扑克的视频| 国产真人三级小视频在线观看| svipshipincom国产片| 岛国视频午夜一区免费看|