基于神經(jīng)網(wǎng)絡(luò)的短電弧覆蓋次數(shù)最優(yōu)化模型研究

李 杰， 孫 煒， 范凱旋， 商慶清

(南京林業(yè)大學(xué)機(jī)械電子工程學(xué)院，江蘇 南京 210037)

1 研究背景

短電弧加工是指在具有一定氣、液比例且?guī)毫Φ幕旌衔锕ぷ鹘橘|(zhì)的作用下，利用兩個(gè)電極之間產(chǎn)生的受激發(fā)短電弧放電群組或混合有火花放電的放電群組，來(lái)蝕除金屬或非金屬導(dǎo)電材料的一種電加工方法，是一種新型的強(qiáng)焰流、電子流、離子流、弧流混合放電加工方法，屬于特種加工行業(yè)電加工技術(shù)范疇。

短電弧加工技術(shù)主要用于加工各類(lèi)水泥磨輥、立磨輥、渣漿泵葉輪、大型冷軋工作輥、高速線(xiàn)材輥及其他大型工件上的高強(qiáng)度、高硬度的難加工金屬材料，如電焊、等離子堆焊后的表面金屬材料，可以高效率地對(duì)此類(lèi)大型磨輥、軋輥外圓表面進(jìn)行修復(fù)、再制造加工[1]。

作為一種大功率高效電加工方法，短電弧加工過(guò)程中峰值電流可達(dá)上千安培，消耗大量熱量與能量，其間伴隨的能量損耗不容忽視，這些能耗也會(huì)引起電源整機(jī)溫度上升，降低電源內(nèi)部元器件的使用壽命，影響電源整機(jī)穩(wěn)定性，會(huì)對(duì)整個(gè)短電弧加工過(guò)程產(chǎn)生很大影響。因此，研究最少短電弧放電次數(shù)模型是發(fā)展短電弧技術(shù)的關(guān)鍵[2]。

1.1 問(wèn)題引出

當(dāng)使用半徑為60.5 mm的球形短電弧對(duì)半徑為1 000 mm的圓形材料實(shí)現(xiàn)完全覆蓋，且短電弧加工所形成的球形電弧的球心落在待蝕除金屬或非金屬導(dǎo)電材料的平面上時(shí)，此問(wèn)題可以轉(zhuǎn)化為求解平面問(wèn)題，即使用半徑為r的小圓對(duì)半徑為R的大圓區(qū)域?qū)崿F(xiàn)完全覆蓋時(shí)，采用何種方式能夠使半徑為r的小圓個(gè)數(shù)最少，此處r=60.5 mm，R=1 000 mm。

多個(gè)圓形之間的密排方式有相切與相交兩種，其中，以相切的方式進(jìn)行密排必然會(huì)導(dǎo)致出現(xiàn)間隙；因此，本文研究的問(wèn)題轉(zhuǎn)化為小圓之間以何種方式相交時(shí)，能夠令小圓使用個(gè)數(shù)最小，同時(shí)實(shí)現(xiàn)對(duì)大圓的完全覆蓋。多個(gè)圓形相交時(shí)，產(chǎn)生的交點(diǎn)之間的連線(xiàn)能夠形成圓的內(nèi)接封閉多邊形；而當(dāng)多個(gè)大小相同的圓形相交時(shí)，其交點(diǎn)之間的連線(xiàn)所形成的多邊形為正多邊形時(shí)，能夠使圓形之間因相交產(chǎn)生的重復(fù)覆蓋面積達(dá)到最小。

1.2 六邊形覆蓋模型

在已知二維平面的覆蓋問(wèn)題上，利用正n邊形覆蓋特定區(qū)域且要求使用正n邊形的個(gè)數(shù)最少，需要滿(mǎn)足x個(gè)正n邊形的頂點(diǎn)都位于同一個(gè)點(diǎn)上，使得恰好能拼接成360°的一個(gè)圓周。由于正n邊的內(nèi)角度數(shù)為：

(1)

所以在一個(gè)頂點(diǎn)處的x個(gè)正n邊形應(yīng)滿(mǎn)足方程：

(2)

解得：

(3)

其中x必須為整數(shù)，則其解只有以下幾種情況：

n=3，x=6
n=4，x=4
n=6，x=3

(4)

由以上解可以得出，當(dāng)小圓交點(diǎn)連線(xiàn)形成的封閉多邊形為正三角形、正方形、正六邊形時(shí)，頂點(diǎn)都位于同一個(gè)點(diǎn)上，使得恰好能拼接成360°的一個(gè)圓周，3種情況示意圖如圖1所示。

圖1 3種情況示意圖

使用相同大小的圓相交，使每個(gè)小圓內(nèi)的交點(diǎn)分別連接成小圓的內(nèi)接正三角形、正方形和正六邊形，此時(shí)相鄰兩個(gè)小圓的重疊面積占整個(gè)圓面積的占比為：

(5)

當(dāng)n=3、n=4或n=6時(shí)，λn結(jié)果如下：

(6)

相鄰兩個(gè)圓的重疊面積占整個(gè)圓面積的百分比越小，所需要的小圓個(gè)數(shù)越少，式(6)表明正六邊形能夠通過(guò)使用最少結(jié)點(diǎn)覆蓋整個(gè)大圓范圍的面積。因此，當(dāng)短電弧小圓半徑為r時(shí)，以它的內(nèi)接正六邊形對(duì)大圓區(qū)域進(jìn)行覆蓋，可實(shí)現(xiàn)重疊覆蓋率最小的完全覆蓋，即正六邊節(jié)點(diǎn)覆蓋模型。每一個(gè)短電弧的覆蓋率ηs為短電弧小圓內(nèi)接正六邊形面積As與短電弧小圓Al之比，即：

(7)

使用正六邊形密排方式代替小圓對(duì)大圓區(qū)域進(jìn)行完全覆蓋，即要求大圓邊界要在最外側(cè)正六邊形最遠(yuǎn)端頂點(diǎn)內(nèi)。正六邊形的密排方式有以下3種方式：

(1)正六邊形中心位于大圓圓心處，其他正六邊形依次向外擴(kuò)展；

(2)3個(gè)正六邊形頂點(diǎn)位于大圓圓心處，其他正六邊形依次向外擴(kuò)展；

(3)正六邊形頂點(diǎn)位于大圓內(nèi)任意位置且不與大圓圓心重合，其他正六邊形依次向外擴(kuò)展。

由于圓和正六邊形都是中心對(duì)稱(chēng)圖形，當(dāng)大圓圓心在任意位置時(shí)，必有一象限需要額外六邊形覆蓋，因此方式3用來(lái)完全覆蓋大圓的正六邊形個(gè)數(shù)無(wú)法達(dá)到最小值，因此本節(jié)著重討論方式1與方式2的六邊形覆蓋情況。方式1和方式2的簡(jiǎn)要示意圖如圖2所示。

圖2 兩種方式示意圖

1.3 模型求解

使用AutoCAD軟件繪圖驗(yàn)證示意如圖3所示，對(duì)半徑為1 000 mm的大圓用外接圓為60.5 mm的正六邊形進(jìn)行填充，此時(shí)正六邊形中心與大圓圓心重合，實(shí)現(xiàn)對(duì)大圓的完全覆蓋。對(duì)圖4中的六邊形的個(gè)數(shù)進(jìn)行計(jì)數(shù)，得367個(gè)，可以證明經(jīng)過(guò)MATLAB計(jì)算所獲得的程序能夠得到與實(shí)際值相同的結(jié)果，程序可靠。

圖3 正六邊形中心與大圓圓心重合

當(dāng)正六邊形的頂點(diǎn)位于大圓圓心處時(shí)，經(jīng)MATLAB計(jì)算得出最少需要378次半徑為60.5 mm的短電弧加工才能夠?qū)Π霃綖? 000 mm的圓形材料實(shí)現(xiàn)完全覆蓋。

使用AutoCAD軟件繪圖驗(yàn)證如圖4所示，對(duì)半徑為1 000 mm的大圓用外接圓為60.5 mm的正六邊形進(jìn)行填充，此時(shí)3個(gè)正六邊形頂點(diǎn)與大圓圓心重合，實(shí)現(xiàn)對(duì)大圓的完全覆蓋。對(duì)圖4中的六邊形的個(gè)數(shù)按行進(jìn)行計(jì)數(shù)，得到共需378個(gè)短電弧才能夠?qū)崿F(xiàn)完全覆蓋。

圖4 三個(gè)正六邊形頂點(diǎn)與大圓圓心重合

1.4 模型對(duì)比

兩種模型結(jié)果見(jiàn)表1。

表1 兩種模型結(jié)果

由表1可以看出，相比模型二，模型一不僅正六邊形使用個(gè)數(shù)更少，且對(duì)正六邊形的利用率更高。因此，當(dāng)使用半徑為60.5 mm的短電弧加工方式加工半徑為1 000 mm的圓形材料時(shí)，模型一能夠使加工次數(shù)達(dá)到最少，最少的加工次數(shù)為367次。

2 問(wèn)題拓展

同理，當(dāng)短電弧半徑改為65.5 mm、70.5 mm和75.5 mm，圓形材料半徑不變時(shí)，依然可以使用正六邊形模型。

2.1 短電弧半徑為65.5 mm

短電弧半徑為65.5 mm時(shí)，首先使用MATLAB進(jìn)行求解，圓心位置不同時(shí)，得到最小加工次數(shù)分別為313個(gè)和321個(gè)；使用AutoCAD繪圖軟件繪圖驗(yàn)證，大圓圓心位于正六邊形中心和大圓圓心位于正六邊形頂點(diǎn)時(shí)，求解結(jié)果對(duì)比見(jiàn)表2。

表2 短電弧半徑為65.5 mm時(shí)的結(jié)果

由上表可以得出，短電弧半徑為65.5 mm的時(shí)候，至少需要313次短電弧加工才可以對(duì)半徑為1 000 mm的圓形材料實(shí)現(xiàn)表面完全覆蓋加工，此時(shí)大圓圓心位置在正六邊形中心位置，即采用模型一方式比模型二方式更好。

2.2 短電弧半徑為70.5 mm

短電弧半徑為70.5 mm時(shí)，首先使用MATLAB進(jìn)行求解，得到最小次數(shù)分別為283個(gè)和282個(gè)，使用AutoCAD繪圖軟件進(jìn)行繪圖驗(yàn)證，大圓圓心位于正六邊形中心和大圓圓心位于正六邊形頂點(diǎn)時(shí)，求解結(jié)果對(duì)比見(jiàn)表3。

表3 短電弧半徑為70.5 mm時(shí)的結(jié)果

由上表可以得出，短電弧半徑為70.5 mm的時(shí)候，至少需要282次短電弧加工才可以對(duì)半徑為1 000 mm的圓形材料實(shí)現(xiàn)表面完全覆蓋加工，此時(shí)大圓圓心位置在正六邊形頂點(diǎn)位置，即采用模型二方式比模型一方式更好。

2.3 短電弧半徑為75.5 mm

短電弧半徑為75.5 mm時(shí)，首先使用MATLAB進(jìn)行求解，得到最小次數(shù)分別為253個(gè)和246個(gè)，使用AutoCAD繪圖軟件進(jìn)行繪圖驗(yàn)證，大圓圓心位于正六邊形中心和大圓圓心位于正六邊形頂點(diǎn)時(shí)，求解結(jié)果對(duì)比見(jiàn)表4。

表4 短電弧半徑為75.5 mm時(shí)的結(jié)果

由上表可以得出，短電弧半徑為75.5 mm的時(shí)候，至少需要246次短電弧加工才可以對(duì)半徑為1 000 mm的圓形材料實(shí)現(xiàn)表面完全覆蓋加工，此時(shí)大圓圓心位置在正六邊形頂點(diǎn)位置。

總結(jié)，當(dāng)短電弧半徑分別為60.5 mm、65.5 mm、70.5 mm和75.5 mm時(shí)，短電弧最少覆蓋次數(shù)見(jiàn)表5。

表5 短電弧最少覆蓋次數(shù)

如何更好的獲得短電弧最少覆蓋次數(shù)，然后可以根據(jù)此方法優(yōu)化改進(jìn)短電弧加工方式，減少加工成本，一直是短電弧加工制造探索的重點(diǎn)。近年來(lái)，神經(jīng)網(wǎng)絡(luò)因其優(yōu)秀的數(shù)據(jù)擬合與預(yù)測(cè)性能，已在許多領(lǐng)域得到了較好的應(yīng)用，如果能將這種技術(shù)應(yīng)用于短電弧最少覆蓋次數(shù)的研究，將會(huì)大大減少短電弧加工成本的支出。

3 神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型

3.1 神經(jīng)網(wǎng)絡(luò)理論

BP(BackPropagation)神經(jīng)網(wǎng)絡(luò)是一種在生物神經(jīng)網(wǎng)絡(luò)的啟發(fā)下建立的數(shù)據(jù)處理模型，并依據(jù)逆向算法訓(xùn)練的多層前饋型神經(jīng)網(wǎng)絡(luò)。BP網(wǎng)絡(luò)結(jié)構(gòu)如圖5所示，其結(jié)構(gòu)通常由輸入層、隱含層和輸出層組成。信息正向傳播，而誤差反向傳播。在正向傳播過(guò)程中，向輸入層中輸入n組數(shù)據(jù)樣本作為輸入信號(hào)，經(jīng)隱含層計(jì)算后傳遞至輸出層，輸入層與隱含層及隱含層與輸出層之間均采用權(quán)重進(jìn)行連接。當(dāng)在輸出層未能得到預(yù)期的訓(xùn)練結(jié)果時(shí)，計(jì)算輸出層的誤差變化值，并將誤差反向傳播。之后，根據(jù)誤差修改各層之間的權(quán)重，直至輸出信號(hào)達(dá)到期望目標(biāo)，停止訓(xùn)練并輸出結(jié)果。其優(yōu)點(diǎn)可以概括為自學(xué)習(xí)和自適應(yīng)性強(qiáng)、容錯(cuò)率高、具有較高的計(jì)算精度和儲(chǔ)存速度[3]，所以本文將使用BP神經(jīng)網(wǎng)絡(luò)作為預(yù)測(cè)模型。

圖5 BP網(wǎng)絡(luò)結(jié)構(gòu)

3.2 生成訓(xùn)練數(shù)據(jù)

當(dāng)圓形材料R與短電弧小圓半徑r變化時(shí)，最少短電弧覆蓋次數(shù)就是一個(gè)二元方程，其最少次數(shù)與圓形材料R與短電弧小圓半徑r都有關(guān)系，即：

n=f(r，R)

(8)

且圓形材料R與短電弧小圓半徑r是基本參數(shù)，兩者沒(méi)有相關(guān)性。

利用MATLAB生成當(dāng)圓形材料半徑R=1 000 mm時(shí)，短電弧的半徑r從40.5～239.5 mm，計(jì)算出最少短電弧覆蓋次數(shù)，共200個(gè)數(shù)據(jù)，然后當(dāng)圓形材料半徑R=1 001 mm時(shí)，短電弧的半徑r從40.5～239.5 mm生成200個(gè)數(shù)據(jù)，計(jì)算出最少短電弧覆蓋次數(shù)，以此類(lèi)推，直到圓形材料半徑R=1 199 mm時(shí)，總共生成40 000(200*200)組輸入為大圓半徑R與小圓半徑r，輸出為短電弧最少覆蓋次數(shù)的數(shù)據(jù)。部分?jǐn)?shù)據(jù)見(jiàn)表6。

表6 部分訓(xùn)練數(shù)據(jù)

3.3 模型的建立與預(yù)測(cè)

模型建立前，需要確定輸入層、輸出層、隱含層節(jié)點(diǎn)數(shù)等參數(shù)[6]。

對(duì)于輸入層，需要各個(gè)變量之間沒(méi)有相關(guān)性或較低相關(guān)性，圓形材料R與短電弧小圓半徑r是基本參數(shù)，兩者沒(méi)有相關(guān)性，同時(shí)直接影響最少覆蓋次數(shù)[5]。

BP神經(jīng)網(wǎng)絡(luò)的一個(gè)隱含層就可以以任意精度逼近任意的非線(xiàn)性函數(shù)，其隱含層節(jié)點(diǎn)數(shù)目可以根據(jù)式(9)確定，經(jīng)多次實(shí)驗(yàn)，最終確定隱含層節(jié)點(diǎn)數(shù)為5。

(9)

式中：m為隱含層節(jié)點(diǎn)數(shù)目；n為輸入層神經(jīng)元數(shù)目；l為輸出層神經(jīng)元數(shù)目；α為1～10之間的常數(shù)。

將數(shù)據(jù)導(dǎo)入Neural Fitting Tool工具箱，輸入為大圓半徑R與小圓半徑r，輸出為短電弧最少覆蓋次數(shù)，并進(jìn)行隨機(jī)排列，將排列后的數(shù)據(jù)集按70∶15∶15的比例分給訓(xùn)練集、驗(yàn)證集與測(cè)試集。

隨后將這些數(shù)據(jù)集代入一個(gè)具有1層5節(jié)點(diǎn)數(shù)的網(wǎng)絡(luò)，進(jìn)行訓(xùn)練，其網(wǎng)絡(luò)模型結(jié)構(gòu)如圖6所示，模型中各結(jié)構(gòu)參數(shù)見(jiàn)表7，訓(xùn)練過(guò)程如圖7所示，各組訓(xùn)練結(jié)果如圖8所示。

圖6 網(wǎng)絡(luò)模型結(jié)構(gòu)

表7 模型結(jié)構(gòu)參數(shù)取值

圖7 訓(xùn)練過(guò)程

圖8 各組訓(xùn)練結(jié)果

最后可以得到訓(xùn)練組、驗(yàn)證組與測(cè)試組的相關(guān)系數(shù)分別為0.999 82、0.999 83、0.999 82，從而可以得出BP網(wǎng)絡(luò)對(duì)于短電弧最少覆蓋次數(shù)具有精確預(yù)測(cè)能力。

4 結(jié)論

驗(yàn)證并選用正六邊形模型對(duì)短電弧最少覆蓋次數(shù)進(jìn)行求解，利用MATLAB計(jì)算出不同短電弧半徑下最少覆蓋次數(shù)，并生成40 000組訓(xùn)練數(shù)據(jù)，利用BP神經(jīng)網(wǎng)絡(luò)，對(duì)數(shù)據(jù)集進(jìn)行訓(xùn)練與預(yù)測(cè)，結(jié)果預(yù)測(cè)較好，因此可以使用BP神經(jīng)網(wǎng)絡(luò)對(duì)短電弧最少覆蓋次數(shù)進(jìn)行預(yù)測(cè)。