板材多點漸進成形的破裂預測模型的建立

肖洪力　蔡改貧　祁步春　周小磊

摘要：針對方金屬板材在成形過程中易破裂問題，以方錐臺為試驗制件，設計響應曲面試驗方案，對多點漸進成形過程進行數(shù)值模擬，借助DesignExpert，將成形工具頭直徑、進給量、成形角、板厚作為分析因子，厚度減薄率φ與成形深度h作為響應因子，建立多元回歸破裂預測模型，對預測模型得到的最優(yōu)與最差工藝參數(shù)進行數(shù)值模擬與試驗，并對得到的結(jié)果進行對比分析。結(jié)果表明：預測結(jié)果、數(shù)值模擬結(jié)果、試驗結(jié)果在數(shù)值上非常接近，以方錐臺為對象構(gòu)建的破裂預測模型，可用于金屬板材破裂的預測，提高金屬制件成形質(zhì)量。

關鍵詞：多點漸進成形；預測模型；響應曲面；數(shù)值模擬

DOI：10.15938/j.jhust.2018.01.005

中圖分類號： TG306

文獻標志碼： A

文章編號： 1007-2683（2018）01-0023-07

Abstract：Aiming at the problem that square sheet metal is easy to be fracturedexisted in the forming process， we put forward an experimental method of response surface， which uses square cone as experimental material， to simulate the multipoint incremental forming. A multivariate regression fracture prediction model is established， regarding tool diameter， feed， forming angle， sheet thickness as influential factors and thickness reduction ratio（φ）， depth of forming（h） as response factors based on DesignExpert software. A comparison of the calculated results with the test results which use the the optimal and worst process parameters obtained from the prediction model indicates that the prediction results， numerical simulation results and test results are very close to each other. The fracture prediction model constructed by square cone can be used to predict the fracture of sheet metal and improve the forming quality of metal parts.

Keywords：mutipoint incremental forming；prediction model； response surface； numerical simulation

0引言

破裂是金屬板材在漸進成形過程中最常見的問題[1]。針對該問題，國內(nèi)外眾多學者對板材破裂進行了大量的研究。李軍超等[2]通過采用均布策略和順逆相間的加工方式，有效的避免了在數(shù)值模擬過程中由于采用G代碼加工軌跡所造成的板材破裂等失穩(wěn)現(xiàn)象。文[3-4]闡述了縮頸在單點漸進成形中的產(chǎn)生，并對比了兩點和單點漸進成形這兩種成形工藝下制件的成形精度和成形極限，指出破裂成形極限曲線相比傳統(tǒng)的成形極限曲線更能準確的描述破裂。DucToan Nguyen等[5-6]結(jié)合混合硬化法則和Oyane韌性破裂準則通過數(shù)值模擬探討了工具頭直徑和下壓量對鎂合金韌性斷裂的影響，發(fā)現(xiàn)混合硬化法則能更加準確的預測成形極限。Jacob Smith等[7]通過對比SPIF和ADSIF兩種成形工藝，發(fā)現(xiàn)ADSIF板材破裂的產(chǎn)生要滯后于單點漸進成形。李磊等[8]通過四組破裂仿真和試驗所得到的應力應變數(shù)據(jù)，有效的預測了LY12（M）硬質(zhì)鋁板的成形極限。

從圖2中可知，板材多點漸進成形過程中，在成形角和進給量一定的情況下（圖 2θ=64.5°，ΔZ=2.5mm），當初始板厚一定時，方錐臺制件的成形深度隨著工具頭直徑的增大而增大；當成形工具頭直徑一定時，方錐臺制件的成形深度隨著初始板厚的增大而增大；在成形深度為26mm的情況下，若工具頭直徑小于10.7mm，則初始板厚必須大于0.8mm，若工具頭直徑大于10.7mm，則初始板厚必須小于1.1mm，表明方錐臺制件的成形深度一定時，工具頭直徑越大，初始板厚必須相應越小。因此，在板材多點漸進成形過程中，當成形角和進給量一定時，為了充分的發(fā)揮金屬板材的成形性能，必須合理的選擇一定直徑的成形工具頭和一定厚度的金屬板材，才能使制件的成形深度達到最大而不至破裂。

對預測模型中的交互項DC作響應曲面圖和等高線圖，如圖 3所示。

從圖3中可知，板材多點漸進成形過程中，在成形工具頭直徑和進給量一定的情況下（圖中D=10mm，ΔZ=2.5mm），當初始板厚一定時，方錐臺制件的成形深度隨著成形角的增大而而減小；當成形角一定時，方錐臺制件的成形深度隨著初始板厚的增大而增大；在成形深度為35mm的情況下，若成形角小于61.2°，則初始板厚必須大于0.88mm，表明當方錐臺制件的成形深度一定時，成形角隨著初始板厚的增大而增大。因此，在板材多點漸進成形過程中，當成形工具頭直徑和進給量一定時，為了充分的發(fā)揮金屬板材的成形性能，同樣必須合理的選擇一定大小的成形角和一定厚度的金屬板材，才能使制件的成形深度達到最大而不至破裂。

3預測模型的驗證及優(yōu)化

3.1預測模型的驗證

在其他工藝參數(shù)一定的情況下，將預測模型中的成形角作為變量，建立有限元模型[12]。并利用Marc軟件按照表3試驗方案進行數(shù)值模擬。

由表4可知，模擬值與預測值基本吻合，二者的誤差均在±5%以內(nèi)，表明以成形深度h建立的破裂預測模型可靠。

3.2預測模型的優(yōu)化

利用Design Expert中的預測優(yōu)化功能Optimization對預測模型進行最優(yōu)化設計，綜合考慮成形深度和厚度減薄率，對各成形工藝參數(shù)進行尋優(yōu)重要度分配，且因素的重要性隨著度數(shù)的增大而增大，優(yōu)化方案表如表5所示。

4預測模型試驗結(jié)果分析

4.1制件成形深度分析

利用軟件Marc，建立有板材限元模型[13-14]，對兩組工藝參數(shù)進行數(shù)值模擬。將模擬結(jié)果與制件實物及預測結(jié)果進行對比，分析兩組工藝參數(shù)對破裂的影響。數(shù)值模擬與成形結(jié)果如圖4和圖5所示。

由圖4和圖5可知，由試驗和數(shù)值模擬得到的方錐臺制件在外觀形狀和破裂位置上基本一致；在最優(yōu)工藝參數(shù)下，當制件破裂時，通過數(shù)值模擬得到的制件的成形深度為40.9774mm，實際加工得到的制件的成形深度為39.15mm，預測模型得到的制件成形深度為40.7805mm，三者在制件破裂前的成形深度數(shù)值上基本一致，表明預測模型的預測效果良好。

在最差工藝參數(shù)下，當制件破裂時，預測模型計算得到的成形深度為23.8897mm，實際加工得到的成形深度為23.98mm，與預測模型接近，但比最優(yōu)工藝參數(shù)下的成形深度在數(shù)值上要小15.17mm，在板材漸進成形過程中要避免此種參數(shù)組合。

4.2制件厚度減薄率分析

將方錐臺制件對角線各點沿軸線方向自對角線起始點每間隔3.6mm均勻選取12個點利用厚度測量儀進行厚度測量，如圖6（a）所示，并將測量結(jié)果根據(jù)相應公式換算成厚度減薄率，作出各點的厚度減薄率圖，如圖6（b）所示。

由圖6（b）可知，對于厚度減薄率，最優(yōu)工藝參數(shù)在軸向單元8處達到最大，最差工藝參數(shù)在軸向單元5處達到最大，表明在最優(yōu)工藝參數(shù)下制件的破裂得到了延緩；破裂時，最優(yōu)工藝參數(shù)的最大厚度減薄率為64.8%，最差工藝參數(shù)的最大厚度減薄率為56.7%，表明在最優(yōu)參數(shù)組合下得到的制件厚度分布更加均勻，充分發(fā)揮了金屬板材的成形性能。

5結(jié)論

1）借助 Design Expert ，利用響應曲面法建立了方錐臺制件破裂前所能達到的成形深度的破裂預測模型，且利用該破裂預測模型，可得最優(yōu)工藝參數(shù)組合為工具頭直徑D=6mm，進給量ΔZ=2.61mm，成形角θ=62.21°，初始板厚t0=1.15mm，最差工藝參數(shù)組合為工具頭直徑D=6mm，進給量ΔZ=3mm，成形角θ=69°，初始板厚t0=0.5mm，在金屬板材在成形過程中，為工藝參數(shù)的選擇提供了參考。

2）對預測模型所得兩組工藝參數(shù)組合方案進行了試驗驗證，結(jié)合軟件Marc，對比分析了成形深度和厚度減薄率的試驗結(jié)果、數(shù)值模擬結(jié)果、預測結(jié)果，得到三者在數(shù)值上非常接近，達到成形深度時破裂區(qū)位置一致，進一步說明建立的破裂預測模型可靠，且利用該該模型，能夠避免金屬板材成形時破裂的情況，進一步提高了在板材生產(chǎn)中制件的成形質(zhì)量。

參 考 文 獻：

[1]姜志宏， 朱成堯， 黃信建，等. 板材兩點對稱式漸進成形軌跡的數(shù)值模擬研究[J]. 熱加工工藝， 2016，45（5）：167-170.

[2]李軍超， 毛鋒， 周杰.板材單點漸進成形工藝數(shù)值模擬與成形缺陷研究[J].金屬鑄鍛焊技術(shù).2010，39（5）：71-75.

[3]MARIA B.Silva， PETER S. Nielsen， NIELS Bay， et al. Failure mechanisms in singlepoint incremental forming of metals[J].Int J Adv Manuf Technol.2011，56（9）：893-903.

[4]SILVA M.B， MARTINS P.A.F. TwoPoint Incremental Forming with Partial Die： Theory and Experimentation[J].JMEPEG.2012，22（4）：1018-1027.

[5]NGUYEN DucToan， PARK JinGee， KIM YoungSuk. Ductile Fracture Prediction in Rotational Incremental Forming for Magnesium Alloy Sheets Using Combined Kinematic/Isotropic Hardening Model[J]. Metallurgical and Materials Transactions A.2010，41（8）：1983-1994.

[6]NGUYEN DucToan， PARK JinGee， KIM YoungSuk. Combined Kinematic/Isotropic Hardening Behavior Study for Magnesium Alloy Sheets to Predict Ductile Fracture of Rotational Incremental Forming[J].Int J Mater Form . 2010，3（S1）：939-942.

[7]JACOB Smith， RAJIV Malhotra， LIU W. K.， et al. Deformation mechanics in singlepoint and accumulative doublesided incremental forming[J].Int J Adv Manuf Technol.2013，69（5）：1185-1201.

[8]李磊， 周晚林， HUSSAING.金屬板料單點漸進成形極限的數(shù)值模擬預測[J].機械工程學報.2010，46（18）：102-107.

[9]陳蘊弛， 李芳， 王秋成. 基于響應面法的漸進成形參數(shù)優(yōu)化[J]. 鍛壓技術(shù)， 2014， 39（11）：11-15.

[10]王海， 王進， 王廷和，等. 鎂合金板料工具轉(zhuǎn)動漸進成形工藝參數(shù)實驗研究[J]. 鍛壓技術(shù)， 2015， 40（7）：48-52.

[11]周小磊. 金屬板材多點復合漸進成形破裂缺陷研究[D]. 贛州：江西理工大學， 2015.

[12]蔡改貧， 周小磊， 熊洋. 成形角對多點復合漸進成形過程影響的有限元分析[J]. 熱加工工藝， 2016（7）：126-130.

[13]蔡改貧， 劉志剛， 曾艷祥，等. 金屬板材多點復合漸進成形破裂機理分析[J]. 熱加工工藝， 2016，45（15）：118-121.

[14]姜志宏， 曾艷祥， 劉志剛，等. 基于MSC.Marc軟件的板材多點復合漸進成形缺陷研究[J]. 熱加工工藝， 2016，45（15）：156-159.

（編輯：關毅）