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

    基于回歸方法和神經(jīng)網(wǎng)絡(luò)的電弧增材制造單道成形形貌預(yù)測

    2023-02-24 09:25:08李峰光曹宜發(fā)郭睿姜淑馨劉建永胡勝波戎博川
    精密成形工程 2023年2期
    關(guān)鍵詞:熔深增材電弧

    李峰光,曹宜發(fā),郭睿,姜淑馨,劉建永,胡勝波,戎博川

    基于回歸方法和神經(jīng)網(wǎng)絡(luò)的電弧增材制造單道成形形貌預(yù)測

    李峰光1,曹宜發(fā)1,郭睿1,姜淑馨1,劉建永1,胡勝波1,戎博川2

    (1.湖北汽車工業(yè)學(xué)院 材料科學(xué)與工程學(xué)院,湖北 十堰 442002;2.北京市房山區(qū)特種設(shè)備檢測所機(jī)電一室,北京 102401)

    針對電弧增材制造技術(shù)實(shí)際應(yīng)用中工藝參數(shù)選取困難和成形結(jié)果難預(yù)測的問題,確定高效、準(zhǔn)確的電弧增材制造單道成形形貌預(yù)測的數(shù)學(xué)方法,以快速、方便地選取絲材電弧增材制造工藝參數(shù)并指導(dǎo)成形質(zhì)量控制。在單道單層絲材電弧增材制造實(shí)驗(yàn)的基礎(chǔ)上,采用多種回歸方法和神經(jīng)網(wǎng)絡(luò)方法分別建立焊接電流、電壓和焊槍移動速度等多個(gè)工藝參數(shù)與增材層寬度、增材層高度及熔池深度等成形形貌參數(shù)之間的數(shù)學(xué)關(guān)系模型。電弧增材制造單道成形形貌與焊接電流、電壓和焊槍移動速度顯著相關(guān),且各參數(shù)間存在非線性交互作用;采用多元線性回歸法可較準(zhǔn)確地預(yù)測單道增材層寬度,但對于增材層高度和熔深的預(yù)測效果較差;神經(jīng)網(wǎng)絡(luò)可良好地處理各工藝參數(shù)間復(fù)雜的非線性關(guān)系,其對增材層寬度、增材層高度和熔深的預(yù)測平均誤差率分別為4.17%、6.60%和7.01%,顯著優(yōu)于多元線性回歸法。采用神經(jīng)網(wǎng)絡(luò)法可以準(zhǔn)確預(yù)測電弧增材制造單道成形的形貌參數(shù),進(jìn)而指導(dǎo)增材制造工藝參數(shù)的選取和成形質(zhì)量的控制。

    電弧增材制造;成形形貌;神經(jīng)網(wǎng)絡(luò);H13鋼;回歸方法

    絲材電弧增材制造(Wire and Arc Additive Manufacturing,WAAM)是一種快速近凈成形制造工藝,因其材料利用率和成形效率高、制造成本低且增材件組織結(jié)構(gòu)致密、力學(xué)性能好等突出優(yōu)勢而廣泛應(yīng)用于汽車、航空航天等領(lǐng)域的零部件制造和修復(fù)[1-4]。但目前絲材電弧增材制造技術(shù)仍存在工藝參數(shù)難選取、成形件性能難調(diào)控和成形結(jié)果難預(yù)測等問題[5-8]。根據(jù)不同的成形件結(jié)構(gòu)和性能要求,研究者一般需要進(jìn)行大量的單道單層電弧增材實(shí)驗(yàn),通過單道成形形貌和性能的分析來確定合理的工藝參數(shù),進(jìn)而制定完整的增材制造工藝方案[9-13],所需物力、人力和時(shí)間成本均較高。部分學(xué)者在進(jìn)行激光增材制造研究時(shí)提出,可在實(shí)驗(yàn)的基礎(chǔ)上通過數(shù)學(xué)方法建立工藝參數(shù)與熔覆層形貌之間的關(guān)系模型,對熔覆層形貌進(jìn)行預(yù)測以快速確定工藝參數(shù),進(jìn)而降低制造成本。相關(guān)學(xué)者開展了大量研究工作,證明了隨機(jī)森林回歸[14]、遺傳算法[15-16]、線性回歸[17]和神經(jīng)網(wǎng)絡(luò)[18-21]等方法均能較好地預(yù)測激光增材制造熔覆層形貌,其中線性或非線性回歸方法和神經(jīng)網(wǎng)絡(luò)應(yīng)用最為廣泛[22-24],但相關(guān)數(shù)學(xué)方法在電弧增材制造形貌預(yù)測領(lǐng)域應(yīng)用仍較少。石玗等[1]對金屬電弧增材成形控制研究現(xiàn)狀進(jìn)行了較為詳盡的文獻(xiàn)調(diào)研研究,認(rèn)為增材層寬度和高度等關(guān)鍵形狀參數(shù)的精確控制是電弧增材制造成型控制最為重要的技術(shù)。目前,相關(guān)研究者已經(jīng)建立了基于模糊控制方法的焊縫成形尺寸控制模型[25]、基于激光視覺傳感器的增材層形貌控制模型[26]和基于單變量階躍響應(yīng)傳遞函數(shù)的增材層尺寸動態(tài)控制模型[27]等多種增材層形貌預(yù)測和控制模型,但相關(guān)模型均或多或少地依賴昂貴的硬件設(shè)備支持,尚未能大規(guī)模推廣應(yīng)用。因此,文中在單道單層電弧增材制造實(shí)驗(yàn)的基礎(chǔ)上,利用多種回歸方法和神經(jīng)網(wǎng)絡(luò)建立電弧增材制造工藝參數(shù)與其增材層形貌參數(shù)之間的數(shù)學(xué)關(guān)系模型,比較各方法預(yù)測的精確度,以期為后續(xù)研究中工藝參數(shù)的快速選取和成形質(zhì)量控制提供參考。

    1 實(shí)驗(yàn)方法

    1.1 實(shí)驗(yàn)設(shè)備及方案

    實(shí)驗(yàn)在自行搭建的數(shù)控絲材電弧增材制造系統(tǒng)上進(jìn)行,該電弧增材制造系統(tǒng)主要由保護(hù)氣系統(tǒng)、焊機(jī)(松下YD-350GR5型,熔化極活性氣體保護(hù)焊)、送絲機(jī)構(gòu)和數(shù)控機(jī)床(西門子802C型)4個(gè)部分組成,關(guān)鍵設(shè)備如圖1所示。其中,焊機(jī)的焊接電流和電壓可以分別獨(dú)立調(diào)節(jié),以控制熱輸入;通過數(shù)控編程控制機(jī)床帶動焊槍以一定速度和路徑運(yùn)動而實(shí)現(xiàn)增材制造。增材實(shí)驗(yàn)所用基板為調(diào)質(zhì)處理后的H13鋼,焊絲名為GLG-H13,其成分如表1所示,直徑1.2 mm。由于本文主要研究的是焊接電流()、電壓()和焊槍移動速度()及焊道長度()4個(gè)參數(shù)對單道成形形貌的影響,結(jié)合筆者前期研究,選擇保護(hù)氣為85%Ar+15%CO2;焊接電流變化范圍為80~160 A(送絲速度在1.7~3.7 m/min范圍內(nèi)與焊接電流自動匹配),焊接電壓變化范圍為18~22 V,焊槍移動速度變化范圍為2~6 mm/s,焊接長度變化范圍為50~150 mm;采用L25(56)的4因素5水平正交實(shí)驗(yàn)方案安排實(shí)驗(yàn),如表2所示。

    表1 GLG-H13焊絲成分

    Tab.1 Composition of GLG-H13 welding wire wt.%

    1.2 實(shí)驗(yàn)結(jié)果測量方法

    電弧增材制造實(shí)驗(yàn)結(jié)束后,在每個(gè)單道增材層中部沿垂直于焊槍移動方向采用線切割方式取下試樣,采用400#~1500#砂紙進(jìn)行逐級打磨拋光,經(jīng)無水乙醇清洗后,浸入(HF):(HNO3):(H2O)=1:3:7的腐蝕液中腐蝕10 s,磨樣并機(jī)械拋光后在體視顯微鏡下觀察增材層及熔池形貌。典型的熔池和增材層形貌如圖2所示,其中,為增材層寬度,為增材層高度,為熔池深度,對每一試樣均測量上述3個(gè)參數(shù)并填入實(shí)驗(yàn)結(jié)果表中。

    圖1 電弧增材制造系統(tǒng)關(guān)鍵設(shè)備

    表2 增材制造實(shí)驗(yàn)方案

    Tab.2 Experimental scheme of additive manufacturing

    圖2 增材層形貌參數(shù)測量示意圖

    2 結(jié)果與分析

    2.1 實(shí)驗(yàn)結(jié)果

    各參數(shù)組合下電弧增材制造所得增材層宏觀形貌如圖3所示,部分參數(shù)組合下無法正常成形,未在圖中列出。

    圖3 增材層宏觀形貌

    筆者團(tuán)隊(duì)在激光增材制造研究中發(fā)現(xiàn),由增材熱輸入功率除以焊接速度計(jì)算得到的線能量L(即/,單位J/mm)與增材制造形貌間存在一定的關(guān)系[28],因此,文中針對每個(gè)參數(shù)組合均計(jì)算線能量并測量、和后得到表3所示實(shí)驗(yàn)結(jié)果。

    表3 增材制造實(shí)驗(yàn)方案及結(jié)果

    Tab.3 Experimental scheme and results of additive manufacturing

    Note: "×" indicates that the specimen cannot be formed normally under this parameter combination.

    根據(jù)表3中數(shù)據(jù)可知,各工藝參數(shù)會對增材制造單道成形的增材層寬度、增材層高度和熔深產(chǎn)生影響,在工藝參數(shù)選擇不合理時(shí),甚至無法成形。當(dāng)線能量較低或較高時(shí),均有較大可能出現(xiàn)無法單道成形的情況,但線能量的高低并不能直接決定增材是否成形。例如,在線能量較低時(shí),第3組和第8組實(shí)驗(yàn)中線能量同樣為400.0 J/mm,第8組參數(shù)條件下可以正常成形,但第3組條件下無法成形;同樣,在線能量較高時(shí),第23組(線能量1 066.7 J/mm)無法成形,但與之線能量接近的第10組(線能量1 100.0 J/mm)和線能量更高的第14組(線能量1 260.0 J/mm)均能正常成形。為進(jìn)一步分析線能量與增材層寬度、增材層高度和熔深之間的關(guān)系,將無法成形的7組實(shí)驗(yàn)組合剔除,剩余的18組數(shù)據(jù)按照線能量排序后,繪制線能量L與增材層寬度、增材層高度及熔深的關(guān)系曲線,如圖4所示。

    圖4 線能量對增材層寬度、增材層高度及熔深的影響

    由圖4可知,增材層寬度隨著線能量的增大有逐漸增大的趨勢,但在線能量較低(約350~500 J/mm)時(shí),增材層寬度隨著線能量的增大而有較明顯的波動,在線能量較高(約880 J/mm以上)時(shí),增材層寬度基本保持穩(wěn)定。增材層高度隨著線能量的增大呈現(xiàn)出較明顯的上升趨勢,但上升幅度隨著線能量的增大而逐漸減弱。熔深隨線能量的增大而產(chǎn)生較大的波動,整體上無明顯上升或下降的趨勢,僅在線能量較高(約880 J/mm以上)時(shí)才基本趨于穩(wěn)定。

    2.2 增材層形貌回歸預(yù)測

    2.2.1 基于線能量的單變量回歸預(yù)測

    分別對增材層寬度、增材層高度及熔深3個(gè)結(jié)果數(shù)據(jù)與線能量L之間的關(guān)系進(jìn)行線性回歸分析,所得回歸方程及相應(yīng)2和adj2如表4所示。

    表4 增材層寬度、高度及熔深與線能量的線性回歸關(guān)系

    Tab.4 Linear regression relationship between width, height of additive layer, and depth of weld penetration

    表4中,2和adj2分別表示回歸決定系數(shù)和調(diào)整后回歸決定系數(shù),二者取值均為[0, 1],2和adj2越接近1且二者越接近時(shí)表示回歸效果越好,對線性回歸而言,即表示回歸方程中因變量和自變量之間的線性關(guān)系越明顯。顯然,由表4可知,、和與線能量之間的線性關(guān)系并不明顯,尤其與線能量之間幾乎沒有任何線性關(guān)系。

    同樣,考慮到增材層寬度、高度隨線能量的增大均呈現(xiàn)出逐漸增大且增速逐漸變緩的趨勢,分別回歸了、及與L之間的二次函數(shù)關(guān)系,如圖5所示。由圖中所給出的二次函數(shù)關(guān)系式及回歸質(zhì)量信息可知,、和L之間二次回歸的2和adj2均只有0.6左右,即增材層寬度、高度與線能量之間的二次函數(shù)關(guān)系非常弱;和L之間二次回歸的2和adj2分別為0.173和0.062 7,即熔深和線能量之間幾乎不存在二次函數(shù)關(guān)系。

    圖5 增材制造形貌參數(shù)與線能量的二次函數(shù)關(guān)系

    通過以上分析可知,增材層寬度、增材層高度和熔深與線能量之間的線性和二次回歸結(jié)果質(zhì)量均較差,難以直接用線能量來表征電弧增材制造單道成形的關(guān)鍵質(zhì)量信息。

    2.2.2 多元線性回歸預(yù)測

    利用多元線性回歸分析可以建立增材制造工藝參數(shù)與單道成形質(zhì)量參數(shù)之間的定量函數(shù)關(guān)系。多元回歸方程建立的基本過程為:將電流、電壓、焊槍移動速度、焊道長度作為自變量參數(shù),將增材層寬度、增材層高度和熔深分別作為目標(biāo),建立各目標(biāo)與自變量的定量化函數(shù)方程;而后檢查各自變量對方程影響的顯著性,如某個(gè)自變量的影響不顯著(取置信度95%),則剔除該參數(shù)后重新回歸直至得到最優(yōu)方程。

    在進(jìn)行多元線性回歸預(yù)測時(shí),將18組數(shù)據(jù)劃分為兩部分,一部分進(jìn)行回歸計(jì)算建立預(yù)測方程,另一部分用于驗(yàn)證回歸方程預(yù)測的準(zhǔn)確性。根據(jù)線能量分布均勻的原則,選取第8、11、24這3組數(shù)據(jù)作為驗(yàn)證數(shù)據(jù),其余15組用于建立回歸方程。經(jīng)逐步回歸后可得到3個(gè)目標(biāo)值與各工藝參數(shù)之間的最佳回歸方程及其回歸決定系數(shù)和調(diào)整后的回歸決定系數(shù),如表5所示。

    表5 多元線性回歸結(jié)果

    Tab.5 Results of multiple linear regression

    由表5中的回歸方程及其回歸質(zhì)量參數(shù)可以看出,多元線性回歸法對于增材層寬度的擬合效果較好,對于增材層高度的擬合效果次之,而對于熔深的擬合效果則較差。將第8、11、24這3個(gè)驗(yàn)證組的增材制造工藝參數(shù)數(shù)據(jù)代入上述回歸方程中,可以得出多元回歸預(yù)測的增材層寬度'、增材層高度'和熔深',并與實(shí)驗(yàn)實(shí)測結(jié)果進(jìn)行對比,如表6所示。表中、分別為多元回歸法預(yù)測的'、'和'的誤差率,負(fù)值表示預(yù)測值小于實(shí)驗(yàn)值。誤差率計(jì)算方法如式(1)所示。

    式中:為誤差率;p為預(yù)測值;為實(shí)驗(yàn)值。誤差率平均值為各誤差率絕對值的均值。

    由表6可知,多元回歸模型對于增材層寬度的預(yù)測值較精確,且不同工藝參數(shù)下預(yù)測結(jié)果誤差率波動不大,平均預(yù)測誤差率僅為4.66%。但對于不同工藝參數(shù)下增材層高度和熔深的預(yù)測結(jié)果誤差波動較大,例如第8組增材層高度預(yù)測誤差率僅?1.59%,但第11組達(dá)到了16.52%。同時(shí),增材層高度和熔深預(yù)測的平均誤差率均較高,分別達(dá)到了8.07%和11.79%。此外,由各擬合方程可以看出,電流、電壓和焊槍移動速度3個(gè)參數(shù)均出現(xiàn)在增材層寬度和高度的預(yù)測方程中,即這兩個(gè)參數(shù)對增材層寬度及高度產(chǎn)生顯著的影響作用;而焊道長度在3個(gè)方程中均未出現(xiàn),即焊道長度對增材層形貌的影響并不顯著;對于熔深而言,僅有電流這一個(gè)參數(shù)的影響是顯著的,但考慮到該方程擬合度不高,并結(jié)合焊接基本原理,認(rèn)為電壓和焊接速度對熔深應(yīng)該有較明顯的影響,但各參數(shù)間存在非線性的交互作用,導(dǎo)致其多元線性回歸效果不佳。

    表6 實(shí)驗(yàn)結(jié)果與多元線性回歸預(yù)測結(jié)果對比

    Tab.6 Comparison of experimental results and predicted results by multiple linear regression

    以上分析說明,無論是基于線能量的單變量回歸模型還是基于焊接電流、電壓及焊槍移動速度的多元線性回歸模型,均難以實(shí)現(xiàn)對3個(gè)增材制造形貌關(guān)鍵參數(shù)的準(zhǔn)確預(yù)測。這是由于增材制造成形時(shí),電流、電壓和焊槍移動速度等工藝參數(shù)共同對焊接熱源形狀、熱輸入、焊絲和基材熔化及凝固、增材制造體系散熱等產(chǎn)生相互耦合的影響,這些參數(shù)相互間的影響是高度非線性的,難以用常規(guī)的線性回歸或高階多項(xiàng)式回歸方法得到增材制造工藝參數(shù)和最終成形質(zhì)量參數(shù)間的關(guān)系。

    2.3 增材層形貌神經(jīng)網(wǎng)絡(luò)預(yù)測

    2.3.1 神經(jīng)網(wǎng)絡(luò)基本原理及預(yù)測方法

    圖6顯示了一種典型的人工神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)。人工神經(jīng)網(wǎng)絡(luò)主要由輸入層、隱含層和輸出層構(gòu)成,每層中的圓圈稱為人工神經(jīng)元,其中,輸入層用于接收變化的各影響因素信息,隱藏層(可以有多層)將輸入層傳遞來的信息轉(zhuǎn)換為輸出層可以使用的信息,輸出層則根據(jù)上一個(gè)隱含層的信息輸出研究者關(guān)心的指標(biāo)結(jié)果;通過不斷調(diào)整各神經(jīng)元之間的連接強(qiáng)度(即權(quán)重),最終可建立輸入信息和輸出指標(biāo)間高精度的對應(yīng)關(guān)系,實(shí)現(xiàn)結(jié)果預(yù)測。神經(jīng)網(wǎng)絡(luò)具有高度的非線性函數(shù)逼近能力、自適應(yīng)學(xué)習(xí)能力和容錯(cuò)能力,因而在非線性系統(tǒng)建模中得到了廣泛應(yīng)用。

    采用神經(jīng)網(wǎng)絡(luò)法進(jìn)行預(yù)測的首要任務(wù)是確定神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu),即確定神經(jīng)網(wǎng)絡(luò)的輸入、輸出和隱含層層數(shù)及每個(gè)隱含層神經(jīng)元的個(gè)數(shù)。由于文中主要研究電流、電壓、焊槍移動速度和焊道長度4個(gè)工藝對增材層寬度、增材層高度和熔深3個(gè)結(jié)果指標(biāo)的影響,其中,通過前文的分析可知,焊道長度對3個(gè)指標(biāo)的影響不顯著,故可確定輸入為3個(gè)參數(shù),輸出為3個(gè)結(jié)果指標(biāo),即可以建立一個(gè)3輸入3輸出的神經(jīng)網(wǎng)絡(luò)來一次性完成增材層寬度、增材層高度和熔深的預(yù)測。但根據(jù)神經(jīng)網(wǎng)絡(luò)原理和筆者前期試算可知,多輸出的神經(jīng)網(wǎng)絡(luò)預(yù)測精度會降低,因而建立了3個(gè)3輸入1輸出的神經(jīng)網(wǎng)絡(luò)分別對增材層寬度、增材層高度和熔深3個(gè)指標(biāo)進(jìn)行預(yù)測。在神經(jīng)網(wǎng)絡(luò)中,隱含層的層數(shù)和每個(gè)隱含層的神經(jīng)元個(gè)數(shù)決定了神經(jīng)網(wǎng)絡(luò)的預(yù)測精度。由于本文預(yù)測問題的復(fù)雜程度較低,因而只使用1個(gè)隱含層。隱含層神經(jīng)元個(gè)數(shù)可由經(jīng)驗(yàn)公式(式(2))計(jì)算得到[17]。

    圖6 神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)示意圖

    式中:為輸入神經(jīng)元個(gè)數(shù);為輸出神經(jīng)元個(gè)數(shù);為[1, 10]之間的常數(shù)。通常而言,的取值需要通過多次神經(jīng)網(wǎng)絡(luò)訓(xùn)練并對比結(jié)果優(yōu)劣來確定,但在SPSS中,軟件可自動完成這一過程,根據(jù)預(yù)測指標(biāo)的不同,的最終取值不同。

    在確定神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)后,使用SPSS軟件搭建神經(jīng)網(wǎng)絡(luò),因輸入、輸出數(shù)據(jù)的單位和數(shù)值量級不同,為保證各變量之間的可比性,需要對輸入、輸出數(shù)據(jù)進(jìn)行歸一化處理,此外,還需要選定隱藏層激活函數(shù)和輸出層激活函數(shù),這些過程均可在SPSS軟件中自動完成。最終采用的歸一化公式如式(3)所示,隱藏層激活函數(shù)和輸出層激活函數(shù)分別為雙曲正切函數(shù)和恒等式。

    式中:'為歸一化結(jié)果;為樣本值;min和max分別為最小樣本值和最大樣本值。

    與多元線性回歸預(yù)測一樣,第8、11、24這3組數(shù)據(jù)不參與神經(jīng)網(wǎng)絡(luò)訓(xùn)練,而用于檢驗(yàn)?zāi)P妥罱K預(yù)測效果,其余15組用于神經(jīng)網(wǎng)絡(luò)訓(xùn)練(在訓(xùn)練過程中,神經(jīng)網(wǎng)絡(luò)自動按照8:2的比例將這15組數(shù)據(jù)劃分為訓(xùn)練集和測試集)。在訓(xùn)練完成后,保存網(wǎng)絡(luò)并將第8、11、24組數(shù)據(jù)代入神經(jīng)網(wǎng)絡(luò)即可得出增材層寬度、增材層高度和熔深的神經(jīng)網(wǎng)絡(luò)預(yù)測值。由于神經(jīng)網(wǎng)絡(luò)中各神經(jīng)元初始權(quán)重值是隨機(jī)賦予的,因而每次神經(jīng)網(wǎng)絡(luò)訓(xùn)練結(jié)果均不一樣,本文針對每個(gè)預(yù)測指標(biāo)均進(jìn)行了兩次神經(jīng)網(wǎng)絡(luò)訓(xùn)練。

    2.3.2 神經(jīng)網(wǎng)絡(luò)預(yù)測精度分析

    每個(gè)指標(biāo)的2次神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果及其誤差率如表7所示。由表7可知,神經(jīng)網(wǎng)絡(luò)的預(yù)測精度與其訓(xùn)練過程有關(guān),不同的訓(xùn)練過程得出的預(yù)測結(jié)果精度差異較大,但若網(wǎng)絡(luò)訓(xùn)練合理,可以得到較高精度的預(yù)測結(jié)果。在實(shí)際應(yīng)用中,可多次訓(xùn)練神經(jīng)網(wǎng)絡(luò),通過驗(yàn)證集的平均誤差等指標(biāo)選擇預(yù)測效果最佳的網(wǎng)絡(luò),文中對于增材層寬度的預(yù)測選第2次訓(xùn)練的網(wǎng)絡(luò),對于增材層高度和熔深則均取第1次訓(xùn)練的網(wǎng)絡(luò),3個(gè)參數(shù)的預(yù)測平均誤差率分別為4.17%、6.60%和7.01%。

    表7 實(shí)驗(yàn)結(jié)果與神經(jīng)網(wǎng)絡(luò)預(yù)測結(jié)果對比

    Tab.7 Comparison of experimental results and predicted results by neural network

    Note: Subscript "1" represents the first neural network prediction training; subscript "2" represents the second neural network prediction training.

    為了進(jìn)一步對比神經(jīng)網(wǎng)絡(luò)與多元回歸法對于增材制造單道成形形貌預(yù)測效果的好壞,將兩種方法對于驗(yàn)證組預(yù)測結(jié)果的誤差率繪制在同一幅圖中,如圖7所示。由圖7可知,對于增材層寬度而言,2種方法的預(yù)測精度均較高,平均誤差率均在5%以內(nèi),神經(jīng)網(wǎng)絡(luò)針對不同組別的預(yù)測結(jié)果誤差波動較大,對于部分組別的預(yù)測誤差大于多元線性回歸(如第8組),但整體而言,神經(jīng)網(wǎng)絡(luò)預(yù)測精度更高。這是由于增材層寬度的變化與焊接電流、電壓和焊槍移動速度之間有較好的多元線性關(guān)系(由表5中的回歸方程可以看出),而神經(jīng)網(wǎng)絡(luò)對于多元線性關(guān)系同樣能進(jìn)行很好的預(yù)測,且其精度更高。

    圖7 預(yù)測結(jié)果誤差率對比

    對于熔深而言,神經(jīng)網(wǎng)絡(luò)對于第8組和第24組的預(yù)測誤差率均明顯低于多元線性回歸,神經(jīng)網(wǎng)絡(luò)的預(yù)測平均誤差率為7.01%,而多元線性回歸預(yù)測的平均誤差率達(dá)到11.79%,神經(jīng)網(wǎng)絡(luò)的預(yù)測精度顯著高于多元線性回歸。這是由于熔深的變化與除電流之外的其他參數(shù)之間的線性關(guān)系不明顯,而是受到各參數(shù)間較復(fù)雜的非線性交互作用的影響,因而多元線性回歸法不能很好地進(jìn)行預(yù)測,但神經(jīng)網(wǎng)絡(luò)法卻可以準(zhǔn)確、高效地處理非線性交互作用,預(yù)測結(jié)果更加精確。同樣,對于增材層高度而言,由于焊接電流、電壓和速度的交互作用會對熔滴的流動性、鋪展性等產(chǎn)生較復(fù)雜的影響,因而多元線性回歸法對于不同組別的預(yù)測結(jié)果精度差異很大,而神經(jīng)網(wǎng)絡(luò)法對于不同組別的預(yù)測結(jié)果精度穩(wěn)定性較好,且平均預(yù)測誤差明顯低于多元線性回歸。

    總體而言,神經(jīng)網(wǎng)絡(luò)法對于增材制造單道成形形貌預(yù)測具有更高的精度,且由于神經(jīng)網(wǎng)絡(luò)的特性,當(dāng)參與訓(xùn)練的數(shù)據(jù)記錄越多時(shí),神經(jīng)網(wǎng)絡(luò)的預(yù)測精度將更高。此外,神經(jīng)網(wǎng)絡(luò)相較于多元線性回歸還可以接受更多的輸入,即考慮更多工藝參數(shù)的影響。在實(shí)際生產(chǎn)中,可以通過部分實(shí)驗(yàn)結(jié)果訓(xùn)練建立工藝參數(shù)與單道成形形貌指標(biāo)間的神經(jīng)網(wǎng)絡(luò),隨著生產(chǎn)的進(jìn)行,可將得到的結(jié)果進(jìn)一步送入神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練以優(yōu)化其性能,最終實(shí)現(xiàn)工藝參數(shù)的選取、優(yōu)化和生產(chǎn)結(jié)果的精確預(yù)測。

    3 結(jié)論

    1)電弧增材制造單道成形形貌與焊接電流、電壓和焊槍移動速度顯著相關(guān),單道焊接長度對成形形貌的影響不顯著;各參數(shù)間存在非線性交互作用,焊接線能量與成形形貌各指標(biāo)間沒有明顯的線性或二次相關(guān)關(guān)系。

    2)單道增材層寬度與焊接電流、電壓和焊槍移動速度間存在明顯的多元線性關(guān)系,采用多元線性回歸法可較準(zhǔn)確地預(yù)測單道增材層寬度(平均預(yù)測誤差率僅為4.66%),但對于增材層高度和熔深的預(yù)測效果較差(平均誤差率分別為8.07%和11.79%)。

    3)神經(jīng)網(wǎng)絡(luò)可良好地處理各工藝參數(shù)間復(fù)雜的非線性關(guān)系,其對于增材層寬度、增材層高度和熔深的預(yù)測平均誤差率分別為4.17%、6.60%和7.01%,顯著高于多元線性回歸法;采用神經(jīng)網(wǎng)絡(luò)法可以準(zhǔn)確預(yù)測電弧增材制造單道成形的形貌參數(shù),進(jìn)而指導(dǎo)增材制造工藝參數(shù)的選取和成形質(zhì)量的控制。

    [1] 石玗, 朱珍文, 張剛, 等. 金屬電弧增材成形控制關(guān)鍵技術(shù)及研究現(xiàn)狀[J]. 材料導(dǎo)報(bào), 2022, 36(12): 135-142.

    SHI Yu, ZHU Zhen-wen, ZHANG Gang, et al. Key Tec-h-n-ology and Status of Metal Arc Additive Morphology Control[J]. Materials Reports, 2022, 36(12): 135-142.

    [2] WANG Ting-ting, ZHANG Yuan-bin, WU Zhi-hong, et al. Microstructure and Properties of Die Steel Fabricated by WAAM Using H13 Wire[J]. Vacuum, 2018, 149: 185-189.

    [3] LIU Zhi-yuan, ZHAO Dan-dan, WANG Pei, et al. Additive Manufacturing of Metals: Microstructure Evolution and Multistage Control[J]. Journal of Materials Science & Technology, 2022, 100: 224-236.

    [4] 姜淑馨, 李峰光. 絲材電弧增材制造技術(shù)的研究與應(yīng)用[J]. 鑄造技術(shù), 2022, 43(5): 369-374.

    JIANG Shu-xin, LI Feng-guang. Research and Application of Wire and Arc Additive Manufacturing Technology[J]. Foundry Technology, 2022, 43(5): 369-374.

    [5] 王華明. 高性能大型金屬構(gòu)件激光增材制造: 若干材料基礎(chǔ)問題[J]. 航空學(xué)報(bào), 2014, 35(10): 2690-2698.

    WANG Hua-ming. Materials' Fundamental Issues of Laser Additive Manufacturing for High-Performance Large Metallic Components[J]. Acta Aeronautica et Astronautica Sinica, 2014, 35(10): 2690-2698.

    [6] SINGH S, SHARMA S K, RATHOD D W. A Review on Process Planning Strategies and Challenges of WAAM[J]. Materials Today: Proceedings, 2021, 47: 6564-6575.

    [7] 盧秉恒. 增材制造技術(shù)——現(xiàn)狀與未來[J]. 中國機(jī)械工程, 2020, 31(1): 19-23.

    LU Bing-heng. Additive Manufacturing—Current Situation and Future[J]. China Mechanical Engineering, 2020, 31(1): 19-23.

    [8] 張金田, 王杏華, 王濤, 等. 電弧增材制造單道單層工藝特性研究[J]. 材料導(dǎo)報(bào), 2020, 34(24): 24132- 24137.

    ZHANG Jin-tian, WANG Xing-hua, WANG Tao, et al. Study on the Processing Characteristics of Single-Bead and Single-Layer in the WAAM[J]. Materials Reports, 2020, 34(24): 24132-24137.

    [9] DEREKAR K S. A Review of Wire Arc Additive Manufacturing and Advances in Wire Arc Additive Manufacturing of Aluminium[J]. Materials Science and Technology, 2018, 34(8): 895-916.

    [10] 許劍. 藥芯絲材脈沖TIG電弧增材制造工藝及數(shù)值模擬研究[D]. 天津: 天津理工大學(xué), 2021: 21-28.

    XU Jian. Research on Process and Numerical Simulation of Flux Cored Wire and Pulsed TIG Arc Additive Manufacturing[D]. Tianjin: Tianjin University of Technology, 2021: 21-28.

    [11] LI Kun, KLECKA M A, CHEN Shu-ying, et al. Wire-Arc Additive Manufacturing and Post-Heat Treatment Optimization on Microstructure and Mechanical Properties of Grade 91 Steel[J]. Additive Manufacturing, 2021, 37: 101734.

    [12] KUMAR A, MAJI K. Numerical and Experimental Investigations on Deposition of Stainless Steel in Wire Arc Additive Manufacturing[J]. International Journal of Manufacturing, Materials, and Mechanical Engineering, 2021, 11(4): 40-56.

    [13] 姜淑馨, 李峰光. 工藝參數(shù)對H13鋼電弧增材制造成形品質(zhì)的影響[J]. 特種鑄造及有色合金, 2022, 42(5): 642-646.

    JIANG Shu-xin, LI Feng-guang. Influence of Process Parameters on Forming Quality of H13 Steel by WAAM[J]. Special Casting & Nonferrous Alloys, 2022, 42(5): 642-646.

    [14] 梁旭東, 王煒, 趙凱, 等. 隨機(jī)森林回歸分析在激光熔覆形貌預(yù)測中的應(yīng)用[J]. 中國有色金屬學(xué)報(bào), 2020, 30(7): 1644-1652.

    LIANG Xu-dong, WANG Wei, ZHAO Kai, et al. Application of Random Forest Regression Analysis in Trace Geometry Prediction of Laser Cladding[J]. The Chinese Journal of Nonferrous Metals, 2020, 30(7): 1644-1652.

    [15] YANG Yang, CAO Long-chao, ZHOU Qi, et al. Multi-Objective Process Parameters Optimization of Laser-Magnetic Hybrid Welding Combining Kriging and NSGA-II[J]. Robotics and Computer-Integrated Manufacturing, 2018, 49: 253-262.

    [16] 溫海駿, 孟小玲, 許向川, 等. 基于神經(jīng)網(wǎng)絡(luò)和遺傳算法的激光熔覆工藝參數(shù)多目標(biāo)優(yōu)化[J]. 應(yīng)用激光, 2019, 39(5): 734-740.

    WEN Hai-jun, MENG Xiao-ling, XU Xiang-chuan, et al. Multi-Objective Optimization of Laser Cladding Process Parameters Based on Neural Network and Genetic Algorithm[J]. Applied Laser, 2019, 39(5): 734-740.

    [17] 范鵬飛, 張冠. 基于線性回歸和神經(jīng)網(wǎng)絡(luò)的金屬陶瓷激光熔覆層形貌預(yù)測[J]. 表面技術(shù), 2019, 48(12): 353-359, 368.

    FAN Peng-fei, ZHANG Guan. Prediction on Geometrical Characteristics of Cermet Laser Cladding Based on Linear Regression and Neural Network[J]. Surface Technology, 2019, 48(12): 353-359, 368.

    [18] LE Xia-wei, HAO Jing-bin, MENG Qing-dong, et al. Prediction of 18Ni300 Laser Cladding Topography Based on Back-Propagation Neural Network and Particle Swarm Optimization[J]. Journal of Physics: Conference Series, 2021, 1775(1): 012009.

    [19] LIU J C, NI L B. Prediction of Laser Clad Parameters Based on Neural Network[J]. Materials Technology, 2012, 27(1): 11-14.

    [20] 劉干成, 黃博. 基于GA-BP神經(jīng)網(wǎng)絡(luò)的鎳基合金熔覆涂層形貌預(yù)測[J]. 應(yīng)用激光, 2018, 38(4): 527-535.

    LIU Gan-cheng, HUANG Bo. Prediction of Ni-Based Alloy Cladding Coatings Topography Based on GA-BP Neural Network[J]. Applied Laser, 2018, 38(4): 527- 535.

    [21] 吳許祥, 王成, 薛華軍, 等. 基于BP神經(jīng)網(wǎng)絡(luò)的激光焊接工藝參數(shù)優(yōu)化及組織性能研究[J]. 鍛壓裝備與制造技術(shù), 2021, 56(6): 77-82.

    WU Xu-xiang, WANG Cheng, XUE Hua-jun, et al. Research on Optimization of Laser Welding Process Parameters and Microstructure Performance Based on BP Neural Network[J]. China Metalforming Equipment & Manufacturing Technology, 2021, 56(6): 77-82.

    [22] 康靚, 米曉希, 王海蓮, 等. 人工神經(jīng)網(wǎng)絡(luò)在材料科學(xué)中的研究進(jìn)展[J]. 材料導(dǎo)報(bào), 2020, 34(21): 21172- 21179.

    KANG Jing, MI Xiao-xi, WANG Hai-lian, et al. Research Progress of Artificial Neural Networks in Material Science[J]. Materials Reports, 2020, 34(21): 21172-21179.

    [23] 米曉希, 湯愛濤, 朱雨晨, 等. 機(jī)器學(xué)習(xí)技術(shù)在材料科學(xué)領(lǐng)域中的應(yīng)用進(jìn)展[J]. 材料導(dǎo)報(bào), 2021, 35(15): 15115-15124.

    MI Xiao-xi, TANG Ai-tao, ZHU Yu-chen, et al. Research Progress of Machine Learning in Material Science[J]. Materials Reports, 2021, 35(15): 15115-15124.

    [24] RAMPRASAD R, BATRA R, PILANIA G, et al. Machine Learning in Materials Informatics: Recent Applications and Prospects[J]. Npj Computational Materials, 2017, 3: 54.

    [25] XU Jian-ning, ZHANG Hua, ZHANG Guang-yun, et al. The Choice of Welding Parameters and Prediction of Weld Seam Dimensions for Welding Rapid Prototyping[C]//2008 IEEE Conference on Cybernetics and Intelligent Systems. Chengdu, China: IEEE, 2008: 976-979.

    [26] XIONG Jun, LI Yan-jiang, YIN Zi-qiu, et al. Determination of Surface Roughness in Wire and Arc Additive Manufacturing Based on Laser Vision Sensing[J]. Chinese Journal of Mechanical Engineering, 2018, 31(1): 74.

    [27] XIONG Jun, YIN Zi-qiu, ZHANG Wei-hua. Forming Appearance Control of Arc Striking and Extinguishing Area in Multi-Layer Single-Pass GMAW-Based Additive Manufacturing[J]. The International Journal of Advanced Manufacturing Technology, 2016, 87(1): 579-586.

    [28] 劉洋. H13鋼汽車熱鍛模具激光仿生強(qiáng)化工藝研究[D]. 十堰: 湖北汽車工業(yè)學(xué)院, 2021: 21-23.

    LIU Yang. Research on the Laser Bionic Strengthening Process of H13 Steel Automobile Hot Forging Die[D]. Shiyan: Hubei University of Automotive Technology, 2021: 21-23.

    Prediction of Single-pass Arc Additive Manufacturing Forming Morphology Based on Regression and Neural Network

    LI Feng-guang1, CAO Yi-fa1, GUO Rui1, JIANG Shu-xin1, LIU Jian-yong1, HU Sheng-bo1, RONG Bo-chuan2

    (1. School of Materials Science and Engineering, Hubei University of Automotive Technology, Hubei Shiyan 442002, China; 2. No.1 Mechanical and Electrical Office of Beijing Fangshan District Special Equipment Testing Institute, Beijing 102401, China)

    The work aims to determine an efficient and accurate mathematical method for predicting the forming morphology of single pass in arc additive manufacturing to solve the difficulties in selecting process parameters and predicting forming results in practical application of arc additive manufacturing technology, so as to select the process parameters of wire arc additive manufacturing quickly and conveniently and guide the forming quality control. A variety of regression methods and neural network methods were used to establish the mathematical relationship model between the multiple process parameters (welding current, voltage, welding torch moving speed, et al) and forming morphology parameters (width of additive layer, height of additive layer and depth of weld penetration) based on the single-pass single-layer wire arc additive manufacturing experiment. The results showed that the single-pass forming morphology of arc additive manufacturing was significantly related to welding current, voltage and moving speed of welding torch, and nonlinear interaction existed between the parameters. The multiple linear regression method could accurately predict the width of single additive layer, but the prediction effect of additive layer height and weld penetration depth was poor. The neural network could handle the complex nonlinear relationship among the process parameters well, its average prediction error rates on the width, height and penetration of the additive layer were respectively 4.17%, 6.60% and 7.01%, which were significantly lower than those of the multiple linear regression method. The neural network method can accurately predict the morphology parameters of arc additive manufacturing single-pass forming, and then guide the selection of additive manufacturing process parameters and the control of forming quality.

    arc additive manufacturing; forming morphology; neural network; H13 steel;regression method

    10.3969/j.issn.1674-6457.2023.02.020

    TG44

    A

    1674-6457(2023)02-0171-09

    2022?07?28

    2022-07-28

    國家自然科學(xué)基金(51604103);湖北省教育廳科研計(jì)劃重點(diǎn)項(xiàng)目(D20221801);湖北省教育廳科研計(jì)劃(Q20211804);湖北汽車工業(yè)學(xué)院“增材制造及表面強(qiáng)化”校級創(chuàng)新團(tuán)隊(duì)(B05)

    National Natural Science Foundation of China (51604103); Key Scientific Research Project of Education Department of Hubei Province (D20221801); Scientific Research Project of Education Department of Hubei Province (Q20211804); Innovation Team Found of Hubei University of Automotive Technology: Additive Manufacturing and Surface Strengthening (B05)

    李峰光(1987—),男,博士,教授,主要研究方向?yàn)樵霾闹圃?、激光表面?qiáng)化。

    LI Feng-guang (1987-), Male, Doctor, Professor, Research focus: additive manufacturing and laser surface strengthening.

    郭睿(1989—),男,碩士,講師,主要研究方向?yàn)樵霾闹圃臁㈣T造工藝優(yōu)化。

    GUO Rui (1989-), Male, Master, Lecturer, Research focus: additive manufacturing and casting process optimization.

    李峰光, 曹宜發(fā), 郭睿, 等. 基于回歸方法和神經(jīng)網(wǎng)絡(luò)的電弧增材制造單道成形形貌預(yù)測[J]. 精密成形工程, 2023, 15(2): 171-179.

    LI Feng-guang, CAO Yi-fa, GUO Rui, et al. Prediction of Single-pass Arc Additive Manufacturing Forming Morphology Based on Regression and Neural Network[J]. Journal of Netshape Forming Engineering, 2023, 15(2): 171-179.

    猜你喜歡
    熔深增材電弧
    動力電池焊接熔深在線監(jiān)測與閉環(huán)控制
    厚板單邊V型坡口對接焊縫熔深控制方法
    石材增材制造技術(shù)研究
    石材(2022年4期)2022-06-15 08:55:02
    故障電弧探測器與故障電弧保護(hù)裝置在工程中的應(yīng)用分析
    2219鋁合金激光電弧復(fù)合焊接及其溫度場的模擬
    激光增材制造仿真過程分析
    下降特性焊接電源中焊接熔深的控制研究
    電子制作(2019年10期)2019-06-17 11:45:00
    我國增材制造技術(shù)的應(yīng)用方向及未來發(fā)展趨勢
    航空電氣系統(tǒng)中故障電弧的分析
    電子制作(2017年22期)2017-02-02 07:10:15
    關(guān)于焊接熔深試驗(yàn)方法的探討
    青春草视频在线免费观看| 超色免费av| 秋霞伦理黄片| 2021少妇久久久久久久久久久| 最近最新中文字幕免费大全7| 如日韩欧美国产精品一区二区三区 | 亚洲av国产av综合av卡| 丝袜喷水一区| 国产片特级美女逼逼视频| 欧美bdsm另类| 天堂中文最新版在线下载| 蜜桃在线观看..| 99热6这里只有精品| av天堂久久9| 最新中文字幕久久久久| 免费观看av网站的网址| 中国美白少妇内射xxxbb| 久久精品国产亚洲网站| 男男h啪啪无遮挡| 国产亚洲欧美精品永久| 精品久久久久久久久亚洲| 美女大奶头黄色视频| 亚洲人成77777在线视频| av视频免费观看在线观看| 特大巨黑吊av在线直播| 欧美精品亚洲一区二区| 国产成人免费无遮挡视频| 黑人猛操日本美女一级片| 国产一区二区在线观看日韩| 亚洲精品视频女| 成年人免费黄色播放视频| 成人午夜精彩视频在线观看| 九九爱精品视频在线观看| 精品一区二区三区视频在线| 热re99久久国产66热| 亚洲情色 制服丝袜| 日本猛色少妇xxxxx猛交久久| 日本vs欧美在线观看视频| 亚洲天堂av无毛| 日本色播在线视频| av播播在线观看一区| 久久久久久久久久久久大奶| 国产av码专区亚洲av| 99re6热这里在线精品视频| 久久久精品免费免费高清| 国产淫语在线视频| 日日啪夜夜爽| videos熟女内射| 在线观看人妻少妇| 十八禁高潮呻吟视频| 中文字幕亚洲精品专区| 婷婷成人精品国产| 少妇人妻精品综合一区二区| 亚洲精品第二区| 国产成人aa在线观看| 免费黄色在线免费观看| 亚洲国产精品成人久久小说| 99九九在线精品视频| 王馨瑶露胸无遮挡在线观看| 插阴视频在线观看视频| 成年人午夜在线观看视频| 免费日韩欧美在线观看| 99久久精品国产国产毛片| av电影中文网址| 国产在视频线精品| 九色亚洲精品在线播放| 婷婷色综合www| 啦啦啦视频在线资源免费观看| 日本av手机在线免费观看| 久久精品久久精品一区二区三区| 王馨瑶露胸无遮挡在线观看| 五月天丁香电影| 国模一区二区三区四区视频| 能在线免费看毛片的网站| 亚洲婷婷狠狠爱综合网| 国模一区二区三区四区视频| 一本久久精品| 国产老妇伦熟女老妇高清| 国产视频内射| 欧美亚洲日本最大视频资源| 高清毛片免费看| 中文字幕免费在线视频6| 亚洲少妇的诱惑av| 一个人免费看片子| 中文字幕久久专区| 国产精品人妻久久久久久| 日韩在线高清观看一区二区三区| 亚洲精品自拍成人| 少妇被粗大的猛进出69影院 | 亚洲欧洲精品一区二区精品久久久 | 久久久久久久国产电影| 亚洲精品国产av蜜桃| 纯流量卡能插随身wifi吗| 黑人巨大精品欧美一区二区蜜桃 | 飞空精品影院首页| 亚洲图色成人| 亚洲国产毛片av蜜桃av| 久久久精品免费免费高清| 伦精品一区二区三区| 国产白丝娇喘喷水9色精品| 久久精品人人爽人人爽视色| 精品国产一区二区三区久久久樱花| 欧美激情极品国产一区二区三区 | 久久久久久久久久成人| 欧美人与善性xxx| 女人久久www免费人成看片| 国产高清不卡午夜福利| 黑丝袜美女国产一区| 你懂的网址亚洲精品在线观看| 人成视频在线观看免费观看| 欧美精品人与动牲交sv欧美| 一区二区三区精品91| av播播在线观看一区| 99热网站在线观看| 午夜激情久久久久久久| 97在线视频观看| 久久综合国产亚洲精品| 午夜福利视频精品| 我要看黄色一级片免费的| 亚洲精品第二区| 狠狠精品人妻久久久久久综合| 精品亚洲成a人片在线观看| 国产亚洲精品第一综合不卡 | 在线观看www视频免费| 久久久久久久大尺度免费视频| 欧美成人午夜免费资源| 人妻 亚洲 视频| 国产在线免费精品| 天天操日日干夜夜撸| 中文欧美无线码| 国产伦精品一区二区三区视频9| 国产成人精品婷婷| 亚洲情色 制服丝袜| 午夜福利视频精品| 国产精品99久久99久久久不卡 | 丝袜脚勾引网站| 精品人妻熟女av久视频| 欧美亚洲日本最大视频资源| 黄色欧美视频在线观看| 精品酒店卫生间| 国产av码专区亚洲av| 一区在线观看完整版| 又粗又硬又长又爽又黄的视频| 18禁在线播放成人免费| 卡戴珊不雅视频在线播放| 曰老女人黄片| 欧美亚洲日本最大视频资源| 久久狼人影院| 欧美激情国产日韩精品一区| 成人亚洲欧美一区二区av| 亚洲精品视频女| 欧美bdsm另类| 丰满乱子伦码专区| 国模一区二区三区四区视频| 热re99久久精品国产66热6| 蜜桃久久精品国产亚洲av| 久久久精品免费免费高清| 国产有黄有色有爽视频| 亚洲熟女精品中文字幕| 一本色道久久久久久精品综合| 制服人妻中文乱码| 一区二区三区乱码不卡18| 久久99热6这里只有精品| 啦啦啦啦在线视频资源| 欧美97在线视频| 免费大片黄手机在线观看| 久久女婷五月综合色啪小说| 美女大奶头黄色视频| 国产熟女午夜一区二区三区 | 亚洲av不卡在线观看| 亚洲人成网站在线观看播放| av有码第一页| 久久久亚洲精品成人影院| 日本vs欧美在线观看视频| 国产精品一区二区在线不卡| 亚洲精品日本国产第一区| 亚洲不卡免费看| 亚洲丝袜综合中文字幕| 亚洲经典国产精华液单| 成人手机av| 一级a做视频免费观看| 午夜福利视频精品| 韩国av在线不卡| 国产乱人偷精品视频| 麻豆成人av视频| 人成视频在线观看免费观看| 成年人免费黄色播放视频| 午夜福利网站1000一区二区三区| 日本vs欧美在线观看视频| 国产成人免费无遮挡视频| 妹子高潮喷水视频| 免费看不卡的av| 人成视频在线观看免费观看| 国产成人免费观看mmmm| 国产成人免费无遮挡视频| 欧美亚洲日本最大视频资源| 久久毛片免费看一区二区三区| 你懂的网址亚洲精品在线观看| 国产精品偷伦视频观看了| 精品久久久久久久久av| 一本一本综合久久| 蜜桃久久精品国产亚洲av| 国产精品一区二区三区四区免费观看| av播播在线观看一区| 亚洲av电影在线观看一区二区三区| 久久久精品区二区三区| 久久久久国产网址| 一区在线观看完整版| 性色avwww在线观看| 少妇的逼水好多| 欧美日韩精品成人综合77777| 亚洲综合精品二区| 国产精品99久久久久久久久| 一边亲一边摸免费视频| 国产亚洲精品第一综合不卡 | 男女国产视频网站| 新久久久久国产一级毛片| 日产精品乱码卡一卡2卡三| 成人无遮挡网站| 另类亚洲欧美激情| xxx大片免费视频| 另类精品久久| 国产永久视频网站| 插阴视频在线观看视频| 国产精品一区二区三区四区免费观看| 在线观看国产h片| 乱码一卡2卡4卡精品| 亚洲av中文av极速乱| 亚洲国产最新在线播放| 18+在线观看网站| 亚洲伊人久久精品综合| 亚洲欧美中文字幕日韩二区| 在线天堂最新版资源| 欧美性感艳星| 一个人免费看片子| 一区二区av电影网| 91精品国产九色| 亚洲av日韩在线播放| 日本av手机在线免费观看| 久久狼人影院| 成人国语在线视频| 亚洲图色成人| 在线观看免费日韩欧美大片 | 2022亚洲国产成人精品| 美女国产高潮福利片在线看| 你懂的网址亚洲精品在线观看| 观看美女的网站| 国产成人精品一,二区| 最黄视频免费看| a 毛片基地| 看十八女毛片水多多多| 91久久精品国产一区二区三区| 亚洲av在线观看美女高潮| 视频在线观看一区二区三区| 亚洲国产精品成人久久小说| 免费观看在线日韩| 热re99久久国产66热| 日韩成人伦理影院| 精品久久久久久久久av| 亚洲国产精品国产精品| 一区在线观看完整版| 午夜激情av网站| 亚洲成色77777| 美女xxoo啪啪120秒动态图| 国产免费又黄又爽又色| 亚洲综合色惰| 青春草国产在线视频| 满18在线观看网站| 国产成人免费无遮挡视频| 亚洲精华国产精华液的使用体验| 午夜激情久久久久久久| 国产精品麻豆人妻色哟哟久久| 国产欧美日韩综合在线一区二区| 嘟嘟电影网在线观看| 亚洲精品aⅴ在线观看| 国产成人a∨麻豆精品| 亚洲无线观看免费| 亚洲内射少妇av| 国产乱人偷精品视频| 亚洲精品日韩av片在线观看| 国产日韩欧美视频二区| 国产成人精品久久久久久| 免费观看无遮挡的男女| 欧美日本中文国产一区发布| 成人漫画全彩无遮挡| 日本欧美视频一区| 一级毛片 在线播放| 午夜精品国产一区二区电影| 日本av免费视频播放| 亚洲欧美一区二区三区国产| 三级国产精品欧美在线观看| 日韩欧美精品免费久久| 夫妻性生交免费视频一级片| 久久 成人 亚洲| 国产毛片在线视频| 亚洲国产精品一区三区| 亚洲人成网站在线观看播放| 欧美少妇被猛烈插入视频| av在线老鸭窝| 久久久久久伊人网av| 中文字幕亚洲精品专区| 五月伊人婷婷丁香| 嘟嘟电影网在线观看| 久久久午夜欧美精品| 国产精品熟女久久久久浪| 亚洲人成77777在线视频| 人人澡人人妻人| 久久这里有精品视频免费| 七月丁香在线播放| 亚洲av在线观看美女高潮| 欧美人与性动交α欧美精品济南到 | 九草在线视频观看| 免费观看性生交大片5| 亚洲欧美成人精品一区二区| 久久久久人妻精品一区果冻| 97超碰精品成人国产| 亚洲综合色网址| 久久久久久久久久久丰满| 国产一区二区在线观看日韩| 丝袜在线中文字幕| 欧美精品一区二区免费开放| 人妻系列 视频| 永久免费av网站大全| 夜夜骑夜夜射夜夜干| 亚洲成色77777| 一区二区三区四区激情视频| 精品人妻熟女毛片av久久网站| 青春草亚洲视频在线观看| 在线观看免费视频网站a站| 亚洲av.av天堂| 国产深夜福利视频在线观看| 精品久久久久久久久亚洲| 麻豆乱淫一区二区| 男女边摸边吃奶| 美女国产高潮福利片在线看| 欧美精品一区二区大全| 午夜激情av网站| 满18在线观看网站| 亚洲国产精品成人久久小说| 欧美精品一区二区大全| 国产成人aa在线观看| 国产白丝娇喘喷水9色精品| 亚洲av福利一区| 老熟女久久久| 国产熟女午夜一区二区三区 | 国产免费一级a男人的天堂| 夜夜骑夜夜射夜夜干| 成人亚洲欧美一区二区av| 嫩草影院入口| 日韩强制内射视频| 两个人免费观看高清视频| 大片电影免费在线观看免费| 欧美97在线视频| 亚洲精品乱久久久久久| 国产色婷婷99| 欧美精品国产亚洲| av专区在线播放| 亚洲内射少妇av| 青春草视频在线免费观看| 在线观看三级黄色| 国产免费又黄又爽又色| 国产男女内射视频| 两个人免费观看高清视频| 亚洲欧洲日产国产| 蜜桃在线观看..| 国产精品欧美亚洲77777| 一级毛片aaaaaa免费看小| 亚洲av福利一区| 欧美一级a爱片免费观看看| 精品久久蜜臀av无| 精品视频人人做人人爽| 男女高潮啪啪啪动态图| 亚洲国产精品999| 亚洲国产精品一区三区| 91精品三级在线观看| 99九九线精品视频在线观看视频| 99国产精品免费福利视频| 欧美性感艳星| 成年人免费黄色播放视频| 最近中文字幕高清免费大全6| 国产探花极品一区二区| 免费观看性生交大片5| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲第一区二区三区不卡| 最后的刺客免费高清国语| 久久女婷五月综合色啪小说| 一本久久精品| 永久网站在线| 日韩 亚洲 欧美在线| 国产女主播在线喷水免费视频网站| 精品一区二区三区视频在线| 亚洲精品视频女| 日韩一区二区三区影片| 免费观看a级毛片全部| 亚洲av不卡在线观看| av福利片在线| 国产深夜福利视频在线观看| 一本—道久久a久久精品蜜桃钙片| 涩涩av久久男人的天堂| 又粗又硬又长又爽又黄的视频| 久久精品熟女亚洲av麻豆精品| 日韩av免费高清视频| 久久狼人影院| 搡老乐熟女国产| kizo精华| 亚洲色图综合在线观看| 久久久久久久久久人人人人人人| 免费日韩欧美在线观看| 成人毛片a级毛片在线播放| 简卡轻食公司| 亚洲精品久久久久久婷婷小说| 丰满饥渴人妻一区二区三| 日本与韩国留学比较| 欧美丝袜亚洲另类| 国产极品粉嫩免费观看在线 | 如日韩欧美国产精品一区二区三区 | 97超碰精品成人国产| 一级二级三级毛片免费看| 久久精品国产亚洲av涩爱| 日本爱情动作片www.在线观看| av有码第一页| 中文字幕精品免费在线观看视频 | 精品久久蜜臀av无| 成人午夜精彩视频在线观看| 久久人人爽av亚洲精品天堂| 性高湖久久久久久久久免费观看| 最后的刺客免费高清国语| 永久网站在线| 一级黄片播放器| 只有这里有精品99| 如何舔出高潮| 七月丁香在线播放| 秋霞伦理黄片| 午夜免费鲁丝| 亚洲av二区三区四区| 最近手机中文字幕大全| 国产成人精品无人区| 亚洲欧美中文字幕日韩二区| 日本与韩国留学比较| 18禁动态无遮挡网站| 精品一区二区免费观看| 一级爰片在线观看| 成人影院久久| 精品酒店卫生间| av电影中文网址| 永久免费av网站大全| 少妇熟女欧美另类| 免费观看a级毛片全部| 黑人巨大精品欧美一区二区蜜桃 | 美女cb高潮喷水在线观看| 纵有疾风起免费观看全集完整版| 少妇人妻精品综合一区二区| 欧美日韩精品成人综合77777| 精品亚洲乱码少妇综合久久| 久久久久视频综合| 国产不卡av网站在线观看| 黄片播放在线免费| 亚洲经典国产精华液单| 国产片特级美女逼逼视频| 国产极品天堂在线| 另类精品久久| 亚洲第一av免费看| 国产精品欧美亚洲77777| 麻豆成人av视频| 卡戴珊不雅视频在线播放| 高清午夜精品一区二区三区| 亚洲人与动物交配视频| 国产精品无大码| av线在线观看网站| 亚洲综合色网址| 91精品一卡2卡3卡4卡| 中文字幕av电影在线播放| 夜夜看夜夜爽夜夜摸| 简卡轻食公司| 免费观看a级毛片全部| 新久久久久国产一级毛片| 一边亲一边摸免费视频| 91精品三级在线观看| 成人无遮挡网站| 欧美精品一区二区大全| 丝袜在线中文字幕| 国产精品国产三级国产av玫瑰| 久久狼人影院| 久久久久久久精品精品| 大陆偷拍与自拍| 欧美精品亚洲一区二区| 欧美日韩精品成人综合77777| 国产亚洲av片在线观看秒播厂| 欧美国产精品一级二级三级| 欧美激情极品国产一区二区三区 | 国产精品麻豆人妻色哟哟久久| 激情五月婷婷亚洲| 性高湖久久久久久久久免费观看| 在线 av 中文字幕| 欧美精品一区二区免费开放| av在线app专区| 2018国产大陆天天弄谢| 伦精品一区二区三区| 亚洲国产最新在线播放| 国产黄色免费在线视频| av线在线观看网站| 欧美国产精品一级二级三级| 女人久久www免费人成看片| 国产精品嫩草影院av在线观看| 这个男人来自地球电影免费观看 | 成人国产av品久久久| 超色免费av| 国产欧美日韩一区二区三区在线 | 中国美白少妇内射xxxbb| 丝袜美足系列| 国产永久视频网站| 亚洲人成网站在线播| 国产男人的电影天堂91| 在线精品无人区一区二区三| 人妻人人澡人人爽人人| 日韩av在线免费看完整版不卡| 飞空精品影院首页| 少妇丰满av| 久久久久久久大尺度免费视频| 精品久久久久久电影网| 国产精品国产av在线观看| 伦理电影免费视频| 国产精品人妻久久久久久| 免费大片18禁| 日本-黄色视频高清免费观看| 韩国高清视频一区二区三区| 日产精品乱码卡一卡2卡三| 亚洲不卡免费看| √禁漫天堂资源中文www| 国产探花极品一区二区| 国产精品熟女久久久久浪| 国产精品一国产av| 在线观看人妻少妇| 欧美性感艳星| 精品熟女少妇av免费看| 成人二区视频| 久久久久视频综合| 亚洲在久久综合| 99久久精品国产国产毛片| 成人手机av| 国产精品麻豆人妻色哟哟久久| 久久综合国产亚洲精品| 午夜福利在线观看免费完整高清在| 日韩一区二区三区影片| 十分钟在线观看高清视频www| 亚洲无线观看免费| 我要看黄色一级片免费的| 中文欧美无线码| 久久久久久久亚洲中文字幕| 人人妻人人添人人爽欧美一区卜| 午夜福利在线观看免费完整高清在| 亚洲成人一二三区av| 国产精品人妻久久久影院| 人妻一区二区av| 国产精品久久久久成人av| av播播在线观看一区| 国产精品麻豆人妻色哟哟久久| 天美传媒精品一区二区| 午夜福利网站1000一区二区三区| 三上悠亚av全集在线观看| 国产亚洲午夜精品一区二区久久| 久久女婷五月综合色啪小说| 麻豆乱淫一区二区| 亚洲精品乱久久久久久| 国产一区二区三区综合在线观看 | 精品亚洲成国产av| 国产精品女同一区二区软件| 国产成人精品久久久久久| 国产精品蜜桃在线观看| 男女边吃奶边做爰视频| 久久久久久久久久人人人人人人| 国产精品一二三区在线看| 七月丁香在线播放| 99国产精品免费福利视频| 国产老妇伦熟女老妇高清| 最新中文字幕久久久久| 久久久欧美国产精品| 一级毛片我不卡| 只有这里有精品99| 日本欧美视频一区| 黄色毛片三级朝国网站| 亚洲一区二区三区欧美精品| 在线观看人妻少妇| 精品少妇久久久久久888优播| 欧美精品一区二区大全| 国产在线一区二区三区精| 国产亚洲最大av| 亚洲成人av在线免费| 另类精品久久| 欧美精品国产亚洲| 久久精品人人爽人人爽视色| 青春草国产在线视频| 精品午夜福利在线看| 少妇高潮的动态图| 少妇熟女欧美另类| 亚洲精品色激情综合| 欧美日本中文国产一区发布| 亚洲人与动物交配视频| 中国三级夫妇交换| 免费大片黄手机在线观看| 97精品久久久久久久久久精品| 成人国产av品久久久| 国产不卡av网站在线观看| av免费在线看不卡| tube8黄色片| 久久99热这里只频精品6学生| 久久久久久久久久久久大奶| 新久久久久国产一级毛片| 久久99热这里只频精品6学生| 男女啪啪激烈高潮av片| 新久久久久国产一级毛片| 久久久久久久精品精品| 日韩一本色道免费dvd| 丁香六月天网| av在线播放精品|