朱 虎 王河坡 邢宇飛

沈陽(yáng)航空航天大學(xué)機(jī)電工程學(xué)院，沈陽(yáng)，110136

數(shù)控漸進(jìn)成形不等層間距等高線軌跡生成

朱 虎 王河坡 邢宇飛

沈陽(yáng)航空航天大學(xué)機(jī)電工程學(xué)院，沈陽(yáng)，110136

考慮數(shù)控漸進(jìn)成形中層間距對(duì)成形質(zhì)量的影響，提出了一種基于等殘余高度的、并對(duì)成形角不同曲面區(qū)域生成不等層間距的等高線軌跡生成方法。首先生成等殘余高度的等高線軌跡，然后在大成形角、大層間距區(qū)域中插入等高線來(lái)優(yōu)化層間距。最終生成能滿足等殘余高度要求、并能通過(guò)減小成形力來(lái)減小反彈變形的軌跡。數(shù)值模擬分析與試驗(yàn)結(jié)果表明，該方法生成的不等層間距成形軌跡與等殘余高度軌跡相比，成形件所受的等效應(yīng)變均勻性好，成形件輪廓精度高，能夠獲得較好的工藝質(zhì)量。

板材成形；數(shù)控漸進(jìn)成形；殘余高度；不等層間距

0 引言

金屬板材數(shù)控漸進(jìn)成形是使擠壓工具頭在數(shù)控機(jī)床的控制下按預(yù)先編制的軌跡運(yùn)動(dòng)，對(duì)金屬板材進(jìn)行逐點(diǎn)擠壓使其發(fā)生塑性變形，進(jìn)而漸進(jìn)地完成板材件的成形[1]。由于該技術(shù)不需要專用模具和沖壓設(shè)備，因此具有較高的柔性，可用于小批量多品種生產(chǎn)和樣品試制，在航空、汽車等行業(yè)具有廣泛的應(yīng)用前景[2]。在數(shù)控漸進(jìn)成形中用來(lái)控制擠壓工具運(yùn)動(dòng)的成形加工軌跡對(duì)成形件的尺寸精度、表面質(zhì)量和成形時(shí)間具有巨大的影響[3]。如何充分考慮金屬板材數(shù)控漸進(jìn)成形特性，生成一個(gè)高效、合理的成形軌跡是數(shù)控漸進(jìn)成形技術(shù)的重要研究?jī)?nèi)容。

目前，在金屬板材數(shù)控漸進(jìn)成形中，等高線軌跡由于生成過(guò)程簡(jiǎn)便、適用范圍大的特點(diǎn)得到廣泛應(yīng)用[4]。等高線軌跡的一個(gè)重要參數(shù)是層間距。研究表明，層間距對(duì)成形力[5-6]、成形精度[7-8]、表面質(zhì)量[9-10]、可成形性[11-12]有較大的影響。然而，在金屬板材數(shù)控漸進(jìn)成形中大多采用在Z軸方向等間距的等高線軌跡[13-14]。ATTANASIO[14]和ZHU等[15]分別研究了基于等殘余高度的等高線成形軌跡和螺旋線成形軌跡生成方法，雖然客觀上生成出不等層間距的成形軌跡，但主觀上沒(méi)有考慮層間距對(duì)可成形性的影響，實(shí)際起不到提高成形性的效果。

綜上所述，眾多的研究已表明，在數(shù)控漸進(jìn)成形中成形軌跡層間距對(duì)可成形性和成形質(zhì)量具有較大的影響，但目前依然采用等層間距的等高線軌跡，沒(méi)有給出適應(yīng)不同成形性和成形質(zhì)量要求的不等層間距成形軌跡生成方法，而本文將對(duì)此展開(kāi)研究。

1 不等層間距等高線軌跡生成

1.1算法思路

在數(shù)控漸進(jìn)成形中，層間距越小，板材的塑性變形區(qū)越小，參與變形的板料越少，變形越均勻，成形件表面質(zhì)量越好，但成形效率低。當(dāng)層間距較大時(shí)，雖然能提高板材件的成形效率，但板材的塑性變形和成形力會(huì)比較大，板材容易被擠破，且成形的表面質(zhì)量較差。為適應(yīng)不同的成形性和成形質(zhì)量要求，對(duì)成形角較大的曲面區(qū)域采用較小的層間距(減小成形力，防止大反彈)，而對(duì)成形角較小的曲面區(qū)域采用較大的層間距(節(jié)省時(shí)間)，使得在保證成形質(zhì)量的同時(shí)提高成形效率，即以滿足用戶表面質(zhì)量要求的等殘余高度等高線軌跡為基準(zhǔn)，在板材件容易破裂區(qū)域的成形軌跡中插入等高線軌跡，通過(guò)減小易破裂區(qū)域成形軌跡的層間距來(lái)提高板材件的可成形性。

1.2基于等殘余高度軌跡生成

(1)

(2)

式中，α為STL模型各三角面片處的成形角；zi為等殘余高度軌跡線在成形內(nèi)表面上對(duì)應(yīng)的Z軸方向高度值；i為等殘余高度軌跡的層數(shù)；r為擠壓工具頭半徑；z為STL模型中位于zi處的三角面片法向量的Z坐標(biāo)值。

圖1 刀位點(diǎn)Fig.1 Tool location point

每層等殘余高度軌跡線在成形內(nèi)表面上對(duì)應(yīng)的高度值z(mì)i為上一層等殘余高度軌跡對(duì)應(yīng)的高度值減去要插入的層間距H，即zi=zi+1-H。殘余高度滿足用戶要求的等高線軌跡的層間距H可由下式[16]求出：

(3)

式中，h為滿足用戶要求的殘余高度值。

生成殘余高度滿足用戶要求的等殘余高度軌跡流程如下。

(1)求出待成形件STL模型內(nèi)表面在Z軸方向的最大高度值z(mì)max和最小高度值z(mì)0，設(shè)定步距Δz。

(2)從zmax-Δz高度位置開(kāi)始，以步距Δz沿z軸方向豎直向下搜索STL模型內(nèi)表面上的三角面片，并求出該高度處各三角面片的成形角α。

(3)根據(jù)式(3)求出對(duì)應(yīng)的殘余高度值h*，并與給定的殘余高度值h比較，如果h*

(4)如果h*>h，利用插值法在高度值z(mì)*和z*-Δz區(qū)域之間，找到殘余高度值h*-h≤ε的高度值z(mì)i；如果殘余高度值h*=h，則殘余高度值h*對(duì)應(yīng)的高度值就是zi。高度值z(mì)i-1與zi的差，為殘余高度h對(duì)應(yīng)的等高線軌跡層間距H。以此類推求出各個(gè)層間距H值，進(jìn)而求出各層刀觸點(diǎn)的Z坐標(biāo)值z(mì)i。

(6)將刀位點(diǎn)沿Z軸負(fù)方向偏置成形工具頭半徑r距離后，按照相鄰關(guān)系依次連接得到圖2所示的等殘余高度軌跡。

圖2 等殘余高度軌跡Fig.2 Path with constant scallop height

1.3等殘余高度軌跡層間距優(yōu)化

在數(shù)控漸進(jìn)成形時(shí)，擠壓工具頭對(duì)板材件的成形力隨著層間距的增大而增大。以等高線軌跡為例，在成形角大的區(qū)域，擠壓工具頭每走一層，層間距越大板材減薄區(qū)域就越大，板材就越容易破裂。為了減小板材在大成形角、大層間距區(qū)域破裂的可能性，提高板材件的可成形性，在大成形角、大層間距區(qū)域插入一系列的等高線，通過(guò)減小層間距來(lái)減小成形力，從而提高板材件的可成形性。在生成的一系列等殘余高度軌跡的大成形角、大層間距區(qū)域中插入等高線的流程如下。

(1)設(shè)定極大層間距Hlimt和極大成形角αlimt，即在板材件的成形角大于αlimt，層間距大于Hlimt的等殘余高度軌跡中插入等高線軌跡，從而減小該區(qū)域的層間距，提高板材件的可成形性。

(2)先找出每層等殘余高度軌跡的高度值，求出層間距Hi，判斷每層的層間距Hi與極大層間距Hlimt的大小。如果|zi-zi-1|=Hi>Hlimt，則提取出兩層等高線之間的三角面片，再判斷提取區(qū)域內(nèi)三角面片的最大成形角αi。如果αi>αlimt，則對(duì)該區(qū)域再次進(jìn)行分層處理，即在該區(qū)域的兩層等高線軌跡之間插入等高線軌跡。

(3)根據(jù)下式計(jì)算在兩層等高線軌跡之間需要插入的等高線軌跡的層數(shù)n：

(4)

其中，H0為基準(zhǔn)層間距，由用戶根據(jù)其加工條件、效率和其他需求來(lái)設(shè)定，用于控制分層間距。

(4)如果n為整數(shù)，則在該兩層等高線軌跡中插入m(m=n)層等間距的等高線軌跡；如果n為非整數(shù)，則求出比n大的最小整數(shù)m，并在該兩層等高線軌跡之間插入m層等間距等高線軌跡；此時(shí)，插入等高線軌跡的層間距Hj可根據(jù)下式求出：

(5)

(5)這樣生成的不等層間距軌跡，既能滿足等殘余高度要求，又能提高可成形性，進(jìn)而有效防止成形件的破裂。圖3所示為等殘余高度軌跡中插入等高線的過(guò)程。

(a)成形軌跡

(b)插入等高線軌跡圖3 不等層間距等高線軌跡生成Fig.3 Generation of the contour path with unequal layer distance

1.4算法應(yīng)用實(shí)例

本文在Windows XP環(huán)境下，利用VC++6.0，OpenGL圖形庫(kù)和C++語(yǔ)言完成了上述算法的軟件系統(tǒng)實(shí)現(xiàn)。該系統(tǒng)具有STL模型可視化、等殘余高度軌跡和不等層間距軌跡生成、NC代碼輸出等功能。圖 4顯示了本系統(tǒng)的用戶界面和一個(gè)待成形件的STL模型。該模型的最大直徑為106 mm，高25 mm，最大成形角為60°，其中區(qū)域1和區(qū)域2為成形角較大的難成形區(qū)域，如圖 5所示。

圖4 STL模型Fig.4 STL model

圖5 成形區(qū)域Fig.5 Forming areas

圖6所示為以殘余高度0.03 mm生成的等殘余高度等高線軌跡，共46層。該等高線軌跡的最大層間距為0.988 mm，最小層間距為0.077 mm，平均層間距為0.543 mm。

圖6 等殘余高度等高線軌跡Fig.6 Contour path with constant scallop height

對(duì)于以殘余高度為0.03 mm生成的等殘余高度軌跡，設(shè)極大層間距為0.5 mm，極大成形角為45°，則在層間距大于0.5 mm以及成形角大于45°的等殘余高度軌跡中插入等高線軌跡。設(shè)基準(zhǔn)層間距為0.3 mm，則插入等高線軌跡后生成的不等層間距軌跡共有87層，該不等層間距軌跡的最大層間距為0.476 mm，最小層間距為0.077 mm(為等殘余高度軌跡的原始層間距)，平均層間距為0.288 mm。而在成形角比較大的難成形區(qū)域(圖5所示的區(qū)域1和區(qū)域2)插入等高線軌跡后，最大層間距為0.3 mm，最小層間距為0.206 mm。圖7所示為生成的不等層間距軌跡。

圖7 不等層間距等高線軌跡Fig.7 Contour path with unequal layer distance

2 成形過(guò)程的有限元分析

為研究在數(shù)控漸進(jìn)成形中不等層間距軌跡對(duì)成形質(zhì)量的影響，本文以圖5所示的板材件模型為例，采用ANSYS/LS-DYNA有限元分析軟件，分別對(duì)殘余高度為0.03 mm的等殘余高度軌跡和不等層間距軌跡進(jìn)行數(shù)值模擬分析。分析過(guò)程中采用Belystchko-Wong-Chiang算法，板材選用厚度為0.88 mm的1060鋁板并簡(jiǎn)化為單元類型為shell163的殼單元，采用1.5 mm的映射網(wǎng)格劃分方式對(duì)其進(jìn)行網(wǎng)格劃分；直徑為10mm的擠壓工具頭材料選用W6Mo5Cr4V2高速鋼，支撐模材料選用GCr15軸承鋼，將擠壓工具和支撐簡(jiǎn)化為solid164的實(shí)體單元，并分別采用1.5 mm和4 mm的自由網(wǎng)格劃分方式對(duì)其進(jìn)行網(wǎng)格劃分。各材料的力學(xué)性能參數(shù)如表1所示。

表1 材料力學(xué)性能參數(shù)

在有限元分析中無(wú)法將復(fù)雜的成形軌跡直接加載到擠壓工具頭上，因此需要將各擠壓刀位點(diǎn)的X、Y、Z坐標(biāo)及對(duì)應(yīng)的時(shí)間t分別保存到文本文件中，然后利用ANSYS軟件的數(shù)組讀入功能導(dǎo)入數(shù)據(jù)，從而得到圖 8所示的軌跡加載曲線。

圖9所示為用不等層間距軌跡和等殘余高度軌跡成形時(shí)的等效應(yīng)變分布。圖9a中板材件的等效應(yīng)變?cè)?.000 15～1.153范圍內(nèi)變化，等效應(yīng)變分布較均勻。圖9b中板材件的等效應(yīng)變?cè)?.0006～0.93范圍內(nèi)變化，有效應(yīng)變小且分布不均勻。為分析成形軌跡在數(shù)值模擬中對(duì)輪廓精度的影響，以模型X=0的中間輪廓線上點(diǎn)的坐標(biāo)值，在Excel中繪制出X=0截面輪廓曲線，如圖10所示。為進(jìn)一步比較成形輪廓曲線與理論模型輪廓曲線的差異，在Excel中擬合出有限元分析模型輪廓與理論模型輪廓在Z方向的偏差曲線，如圖11所示。利用不等層間距軌跡模擬成形的最大偏差為1.06 mm，最小偏差為0.28 mm，平均偏差為0.517 mm；等殘余高度軌跡模擬成形的最大偏差為2.04 mm，最小偏差為0.72 mm,平均偏差為1.09 mm。通過(guò)對(duì)比可知，采用不等層間距軌跡擠壓板材時(shí)，板材件的Z方向偏差小，輪廓精度高，等效應(yīng)變大且分布均勻。采用不等層間距軌跡擠壓板材時(shí)，板材件的輪廓精度高，等效應(yīng)變大且分布均勻。

(a)X軸方向

(b)Y軸方向

(c)Z軸方向圖8 軌跡加載曲線Fig.8 Path loading curves

(a)不等層間距軌跡成形應(yīng)變

(b)等殘余高度軌跡成形應(yīng)變圖9 等效應(yīng)變分布Fig.9 Effective strain distribution

圖10 X=0中間截面輪廓曲線Fig.10 Intermediate section profile curve at X=0

圖11 基于有限元分析的Z軸方向偏差Fig.11 Z-axis direction deviation based on finite element analysis

3 成形試驗(yàn)與測(cè)試分析

3.1成形試驗(yàn)

為了進(jìn)一步驗(yàn)證在數(shù)控漸進(jìn)成形時(shí)，不等層間距軌跡對(duì)板材件成形質(zhì)量產(chǎn)生的影響，以圖5所示的板材件模型為例，利用上述不等層間距軌跡和等殘余高度軌跡分別進(jìn)行了數(shù)控漸進(jìn)成形試驗(yàn)。圖12所示為在模型模具機(jī)上以尼龍板為材料銑削加工支撐的過(guò)程。

圖12 加工支撐Fig.12 Support machining

數(shù)控漸進(jìn)成形試驗(yàn)時(shí)選用厚度為0.88 mm的1060鋁板，直徑為10 mm擠壓工具頭，機(jī)油潤(rùn)滑，主軸停止旋轉(zhuǎn)，進(jìn)給速度為1000 mm/min，成形加工試驗(yàn)過(guò)程如圖13所示，成形件如圖14所示。

圖13 成形試驗(yàn)Fig.13 Forming experiment

(a)不等層間距軌跡成形件(b)等殘余高度軌跡成形件圖14 成形件Fig.14 Formed part

圖15 三坐標(biāo)測(cè)量?jī)xFig.15 Three coordinate measuring instrument

3.2測(cè)試分析

為了檢驗(yàn)不等層間距軌跡和等殘余高度軌跡的成形效果，本文采用三坐標(biāo)測(cè)量?jī)x對(duì)成形件進(jìn)行測(cè)量，如圖15所示。在X軸方向以2 mm為間隔，測(cè)量成形件的中間截面輪廓(即以模型輪廓中心為坐標(biāo)原點(diǎn)的Y=0截面)坐標(biāo)值，然后擬合出成形件的中間截面輪廓曲線并與理論模型輪廓曲線進(jìn)行對(duì)比，如圖16所示。從圖16可以看出，基于不等層間距軌跡的成形件輪廓與理論模型輪廓基本相吻合，且輪廓誤差更小。為進(jìn)一步分析成形件輪廓與理論模型輪廓截面曲線的差異，比較了兩截面輪廓曲線在Z軸方向的偏差，如圖17所示，其中最大偏差為1.17 mm，最小偏差為0.001 mm，平均偏差為0.297 mm; 等殘余高度軌跡成形板材件的最大偏差為1.97 mm，最小偏差為0.3 mm，平均偏差為0.76 mm。在先成形的大成形角區(qū)域，偏差值較大，這是由于在數(shù)控漸進(jìn)成形時(shí)，板材受到拉力和成形工具頭的壓力，會(huì)對(duì)已成形的板材區(qū)域產(chǎn)生拉伸作用，而且板材產(chǎn)生回彈，所以成形輪廓會(huì)出現(xiàn)向內(nèi)偏移的現(xiàn)象，這是在數(shù)控漸進(jìn)成形時(shí)出現(xiàn)的普遍現(xiàn)象。

圖16 Y=0中間截面輪廓曲線對(duì)比Fig.13 Contours lines comparison at Y=0 intermediate section

圖17 Z軸方向偏差Fig.17 Deviation in Z-axis direction

4 結(jié)論

(1)本文提出了基于等殘余高度的不等層間距等高線軌跡生成算法，通過(guò)數(shù)字模擬和成形試驗(yàn)驗(yàn)證了該算法的可行性。

(2)數(shù)字模擬和成形試驗(yàn)結(jié)果表明，不等層間距軌跡與等殘余高度軌跡相比，成形的最大有效應(yīng)變較大且更均勻，并且基于不等層間距軌跡的成形件輪廓誤差更小。

(3)不等層間距軌跡在滿足等殘余高度要求的同時(shí)能夠有效地提高板材件的成形性能和輪廓精度。

[1] LEACOCKA G. The Future of Sheet Metal Forming Research [J]. Materials and Manufacturing Processes, 2012, 27(4): 366-369.

[2] TISZA M. General Overview of Sheet Incremental Forming [J]. Journal of Achievements in Materials and Manufacturing Engineering, 2012, 55(1): 113-120.

[3] RAUCH M, HASCOET J Y, HAMANN J C. A New Approach for Tool Path Programming in Incremental Sheet Forming[J]. International Journal of Material Forming, 2008,1(1): 1191-1194.

[4] ZHU H, LIN W W，BAI J L. An Overview of the Sheet Metal CNC Incremental Forming Toolpath Generation[J]. Advanced Materials Research, 2012, 503/504: 35-39.

[5] BAGUDANCH I, CENTENO G, VALLELLANO C, et al. Forming Force in Single Point Incremental Forming under Different Bending Conditions [J]. Procedia Engineering, 2013, 63: 354-360.

[6] 姚梓萌, 李言, 楊明順,等. 基于響應(yīng)面法的單點(diǎn)增量成形過(guò)程變形能優(yōu)化[J]. 中國(guó)機(jī)械工程, 2017, 28(7): 862-866. YAO Zimeng, LI Yan, YANG Mingshun, et al. Optimization of Deformation Energy in Single Point Incremental Forming Processes Using Response Surface Methodology[J]. China Mechanical Engineering, 2017, 28(7): 862-866.

[7] AMBROGIO G, COSTANTINO I, De NAPOLI L, et al. Influence of Some Relevant Process Parameters on the Dimensional Accuracy in Incremental Forming: a Numerical and Experimental Investigation [J]. Journal of Materials Processing Technology, 2004, 153/154(32): 501-507.

[8] HAM M, JESWIET J. Dimensional Accuracy of Single Point Incremental Forming[J]. International Journal of Material Forming, 2008, 1(1): 1171-1174.

[9] ECHRIF S B M, HRAIRI M. Significant Parameters for the Surface Roughness in Incremental Forming Process [J]. Materials and Manufacturing Processes, 2014,29(6): 697-703.

[10] DURANTEM, FORMISANO A, LANGELLA A. Comparison between Analytical and Experimental Roughness Values of Components Created by Incremental Forming[J]. Journal of Materials Processing Technology, 2010, 210: 1934-1941.

[11] SHANMUGANATANS P, KUMAR V S S. Modeling of Incremental Forming Process Parameters of Al3003(O) by Response Surface Methodology [J]. Procedia Engineering, 2014, 97: 346-356.

[12] MCANULTYT, JESWIET J, DOOLAN M. Formability in Single Point Incremental Forming: a Comparative Analysis of the State of the Art[J]. CIRP Journal of Manufacturing Science & Technology, 2016: 1-12.

[13] JESWIETJ, MICARI F, HIRT G. Asymmetric Single Point Incremental Forming of Sheet Metal[J]. Annals of the CIRP—Manufacturing Technology, 2005,54(2): 623-649.

[14] ATTANASIO A, CERETTI E, GIARDINI C. Optimization of Tool Path in Two Points Incremental Forming [J]. Journal of Materials Processing Technology, 2006, 177(1/3): 409-412.

[15] ZHU H, LIU Z J, FU J H. Spiral Tool-path Generation with Constant Scallop Height for Sheet Metal CNC Incremental Forming [J]. The International Journal of Advanced Manufacturing Technology, 2011, 54(6): 911-919.

[16] 宋修成, 陸彬, 陳軍, 等. 板材漸進(jìn)成形表面質(zhì)量的影響因素分析[J]. 機(jī)械工程學(xué)報(bào), 2013, 49(8): 84-90. SONG Xiucheng, LU Bin, CHEN Jun, et al. Influencing Factor Analysis on the Surface Quality of Incremental Forming Parts[J]. Journal of Mechanical Engineering, 2013, 49(8): 84-90.

(編輯王旻玥)

GenerationofContourPathwithUnequalLayerDistancesinCNCIncrementalForming

ZHU Hu WANG Hepo XING Yufei

College of Mechanical and Electrical Engineering,Shenyang Aerospace University,Shenyang,110136

Considering the influences of layer distance on forming quality in CNC incremental forming, a method to generate the contour path with unequal layer distance was proposed based on the constant scallop height and surface region with different forming angles. First, the forming paths with constant scallop height were generated. Then, the contour paths were inserted in to the region with large forming angle and large layer distance to optimize the layer distances. Finally, the forming paths were generated, which might meet the requirements of scallop height and reduce the forming forces and the springback deformations. The results of numerical simulations and forming experiments show that the unequal layer distance forming path generated by this method has better uniformity of effective strain, higher forming accuracy and better process quality than that of constant scallop height path.

sheet metal forming; CNC incremental forming; scallop height; unequal layer distance

2016-11-08

國(guó)家自然科學(xué)基金資助項(xiàng)目(51175351)；遼寧省自然科學(xué)基金資助項(xiàng)目(2014024011)

TG335.5

10.3969/j.issn.1004-132X.2017.17.019

朱虎，男，1964年生。沈陽(yáng)航空航天大學(xué)機(jī)電工程學(xué)院教授、博士。主要研究方向?yàn)閿?shù)字化設(shè)計(jì)與制造。E-mail：zhuhu10@163.com。王河坡，男，1989年生。沈陽(yáng)航空航天大學(xué)機(jī)電工程學(xué)院碩士研究生。邢宇飛，男，1983年生。沈陽(yáng)航空航天大學(xué)機(jī)電工程學(xué)院講師、博士。