U型鋼冷彎成形過程的有限元仿真

寇長海 高世陽 楊樹寶 許淑芹

摘 要：為了提高U型鋼的抗扭性能和承載能力，改善外彎角的曲率厚度分布均勻性，并減少尖點在圓角部位的出現，本文使用軟件模擬板帶連續冷彎成形過程，以揭示U 型鋼的應力和應變的分布規律。首先，在傳統設計方法的基礎上進行了成型工藝參數和軋輥的設計，確定了成型道次和每次彎曲的角度；然后，利用有限元軟件MSC.Marc，構建了冷彎成型仿真模型，并進行了仿真計算。仿真結果表明：在冷彎過程中，等效應力及等效應變分布呈現彎角處大，立板處次之，而底板處最小的特點，各道次間應力及應變的大小變化較均勻，成型力分布也比較合理。

關鍵詞：冷彎；有限元模擬；變形；成型力

中圖分類號：TH124 文獻標志碼：A? ? 文章編號：1005-2550（2024）02-0064-06

Finite Element Simulation on the Cold Bending Process of U-shaped Steel

KOU Chang-hai1， GAO Shi-yang2， YANG Shu-bao3， XU Shu-qin2

（1.Changyi Agriculture and Rural Bureau Development Service Center.， Weifang 261300， China; 2. Shandong Haoxin Co.， Ltd.， Weifang 261300， China; 3.School of Mechanical Engineering， Anhui University of Technology.， Maanshan 243032， China）

Abstract： In order to improve the torsional performance and support bearing capacity of U-shaped steel， improve the uniformity of bending angle thickness distribution， and reduce the occurrence of sharp points in the rounded corners， software is applied to simulate the continuous cold bending process of plate and strip， revealing the distribution law of stress and strain of U-shaped steel. Firstly， Based on traditional design methods， the forming process parameters and rolling rolls were designed to determine the number of forming passes and each bending angle. Then， with the help of the finite element software MSC. Marc， a cold forming simulation model was established and calculated. The calculation results show that in the cold bending process， the equivalent force and equivalent effect distribution on the workpiece are characterized by large forming angles， followed by vertical plates and smallest bottom plates， and the magnitude changes of stress and strain between passes are more uniform， and the forming force distribution is also reasonable.

Key Words： Cold Bending; Finite Element Simulation; Deformation

引? ? 言

冷彎型鋼作為一種價格比較實惠的斷面型材，因其均勻的斷面形狀、較好的產品質量、制造過程消耗的能源較少和效益較高等優點[1，2]，而被廣泛采用。但由于其成形工藝和工序復雜，成形過程中材料、幾何和接觸邊界等具有多重非線性特點，無法用單純的數學模型描述和精確的計算冷彎過程[3-5]。隨著現代計算技術的發展，有限元仿真技術愈加成熟和完善，利用有限元技術對冷彎成形過程進行模擬越來越受到關注。根據計算結果，能夠較好地把控冷彎成形過程中金屬材料的成形規律，通過模擬過程修正可以得到最優的工序，從而能夠預先了解成形特征、避免在設計軋輥和確定工藝參數時產生不必要的成本浪費，減少了盲目嘗試的風險。目前，國內外學者對此進行了相關研究，為冷彎成型模擬工作提供了參考[6-9]。本文利用MSC.MARC有限元軟件對典型U型件的成型過程進行了模擬分析實踐。

1? ? 分析參數計算

1.1? ?道次參數計算

成型道次的數量級是冷彎設計至關重要的參數之一。隨著其參數量值的增多，所制造的部件精度會提高；圓角曲率的半徑過渡均勻，其力學性能會更好。然而數據過多會影響生產成本和生產效率，型材本身的優勢難以發揮[8]，因此，分析計算出道次參數尤為重要。早期日本教授進行了大量成型型材的斷面結構樣本調研，根據企業實際生產的型材斷面結構與道次參數量分析及整理確定了數據傾向[10]。根據整理出的各個斷面結構數據曲線，進而推斷出了對稱斷面結構的形狀因子與道次數的關系圖，如圖1所示。

圖1中的中心線左側代表帽形截面形狀曲線，中心線右側代表C型截面結構曲線。橫坐標表示成型道次數，中心線右側曲線的截面結構道次數大于左側曲線。圖中相應的縱坐標表示形狀因子函數。根據圖中的形狀因子函數的坐標可以確定出道次數，其因子公式為（1） [11]。

式中：F為左右側兩立邊長度之和、n為彎曲角度數、t為板材厚度。

本文采用典型的U型件進行模擬成型分析，所選零件的具體尺寸和結構如圖2所示。零件所用材料為Q235，壁厚為1.2mm，彎角半徑為3.0mm。

由圖2可知，在相同的形狀因子下不同截面型鋼所需要的成形道次也不相同。該件的形狀因子φ=Fnt=60×2×1.2=144。將計算出的值對照圖1中心線右側的曲線，由此確定成型道次數為N=6。

1.2? ?彎曲角度參數計算

彎曲角度是根據成型道次數進行計算[12]，成型道次不合理，會影響彎曲角度[13]。當道次數過少時，會導致產品的圓角部位應力分布不均勻，產生制造缺陷，還會導致應力分布不均勻。如果道次數量太多，會導致生產工序過多。成型道次數量確定后，需使用余弦計算公式計算每個成型道次的折彎角度，具體參考公式（2）：

式中，θ0折彎成型后最終角度，N表示折彎道次數，i表示每道次數對應的角度。

考慮成型后要抵消材料的回彈，一般要進行過彎處理，過彎角度由公式（3）進行計算：

根據上述公式可求得每道次的理論彎曲角度和實際折彎考慮回彈角度，相關參數如表1所示：

根據相關參數，確定U型件斷面的成形工序如圖3所示。然后運用三維建模軟件進行相關參數建模工作。

2? ? 有限元仿真模型的建立

MSC.MARC有限元軟件適用于復雜的非線性仿真分析問題，而冷彎則是對板帶進行彎曲成形，該工藝所需成型力相對較小 [11，12]。所以本研究中將成形輥、推板及托板等定義為剛性體。軋輥各道次的間距與現生產一致，設定為300mm，輥縫與材料厚度相同。

2.1? ?模型的導入及基本參數的設置

將六道次成形軋輥三維裝配模型以“.IGES”格式導出，在MSC.MARC軟件導入模型文件并進行相關設置，用于模擬的冷彎有限元模型如圖4所示[13-15]。主要部件包括六對軋輥、零部件推板及支撐托板等。其中板帶尺寸400mm× 136.9mm×1.2mm，設定為變形體，網格節點9922個，單元數4800個，均設定為實體單元。工件與軋輥間的摩擦系數設置為0.3，板帶速度為400mm/s，仿真時間設定為5s，總共仿真步數為1000步[16-17]。

2.2? ?仿真結果分析

2.2.1 變形分析

圖5為仿真分析計算結果。該結果顯示了工件在冷彎過程中經過1-6道次軋機時的等效應變分布。在經過6道次軋制之后，工件已經基本成型。從各個道次的等效塑性應變云圖中可以看出，彎角處的等效應變數值最大，說明此處的塑形變形也比較大，冷彎件的立板處應變值次之，底板等部位的等效塑性應變值最小。因此，塑性變形主要集中在冷彎各個工序的彎角位置。同時發現，塑性應變值隨著成形彎曲角度的增加而累計增加，如圖6所示。第一道次最大塑形應變僅為0.07，通過變形積累第六道次時塑形應變達到了0.29。相比較而言，板料其它部分的塑性應變則小得多，這與彎角處變形較大、其它部分變形較小相符。從圖6可以看出，各道次應變大小分布相對平均，沒有出現明顯的急劇增加現象，這表明設計較為合理。

從圖7（a）可見，第1道次剛咬入時，彎角附近等效應變最大，立板區較小而底板中心區最小。彎曲的部分產生了明顯的塑性變形，而立板的兩側則經歷了輕微的塑性變形，其他區域則繼續保持著彈性變形狀態。第2-5道次變形規律類似：兩側立板變形范圍逐漸擴大，鋼帶邊部地區也進入塑性變形狀態，但直到6道次即成型后，在底板及立板頂端仍存在未發生塑形變形的區域。

圖7為模擬的道次1及道次6的等效應力分布圖?？梢娍闯觯簭澖亲冃螀^附近等效應力最大，立板區較小而底板中心區則最小，與等效應變分布相似。脫離軋輥后，鋼帶頭部形狀發生變化，有較明顯的彈性恢復現象。各道次對應的最大等效應力值呈現先增加后平穩最后下降的趨勢，最大值為284MPa，遠小于Q235的抗拉強度，如圖8所示，因此不存在產生變形開裂的危險。

2.3? ?成型力分析

成型力是冷彎成型工藝的關鍵參數之一，圖9為工件經過各道次時的最大成型力?？梢钥闯?，隨著冷彎過程的進行，成型力呈現先增加后減小的趨勢。前三道次時成型力逐漸增加，第三道次成型力達到最大至2350N；后三道次則相反，成型力逐漸減小，最小至1420N。這個規律在彎曲角度分配時，首尾小中間大的原則是一致的。從圖中可以看出，成型過程中成型力變化相對平穩，進一步說明了相關設計較合理。

3? ? 結論

本文以冷彎成形90°的U型材為研究對象，進行了成形工藝參數及軋輥的設計，以設計結果為基礎并借助有限元MSC.Marc軟件進行仿真分析。通過模擬仿真分析，發現冷彎過程彎角附近變形最大，立板區較小而底板中心區則最小。U型件成型過程中的等效應力、應變及成型力分布規律表明了設計的可行性。因此，有限元方法可以為冷彎工藝的設計及生產實踐提供理論指導。

參考文獻：

[1]趙月偉，帥美榮，楚志兵，等.不銹鋼彎曲應力中性層偏移模型構建與實驗研究[J]. 塑性工程學報，2023，30（8）：210-217.

[2]喬治. 哈姆斯著，劉繼英譯.冷彎成形技術手冊[M].北京：化學工業出版社，2008：1～4.

[3]劉淼兒，楊亮，范嘉堃，等.深水柔性管道骨架層的多道次冷彎工藝分析[J].模具技術.2021（3）：40-47.

[4]李子軒，C型鋼的冷彎成形工藝參數分析及成形質量研究[D].寧波大學，2017.

[5]王全華，李培勇，黃啟明，等.型材冷彎成形過程中的回彈和殘余應力[J]. 船舶工程，2022.44（02）：116-122.

[6]顧澤中，孟珂，任琨，等. 支架類槽鋼冷彎成形參數設計及優化[J].鍛壓技術，2022.47（10）：90-95.

[7]付中昊. 典型開口型鋼的孔型設計及連續冷彎成形過程仿真[D].燕山大學，2015.

[8]葉文華. 工程中常用截面的冷彎鋼構件的優化和設計 [D].寧波大學，2018.

[9]B. Shirani Bidabadi，H. Moslemi Naeini，M. Salmani Tehrani，H. Barghikar. Experimental and numerical study of bowing defects in cold roll-formed， U-channel sections[J]. Journal of Constructional Steel Research，2016，118.

[10]小奈弘著，劉繼英譯. 冷彎成形技術[M].北京：化學工業出版社，2007：1～11.

[11]張漢威. 航天薄壁件多道次滾彎成形曲率預測方法研究 [D].東華大學，2022.

[12]T. Welo， Frode Paulsen. Predicting Tube Opalization in Cold Bending： An Analytical Approach[J]. Key Engineering Materials，2015，3974（651）.

[13]孫達，楊龍允，蘇春建.基于彎曲角度分配函數的輥彎成形研究[J].塑性工程學報，2019，26 （ 4） ： 107-112．

[14]吳一昊，黃站立.彎曲工藝對金屬 U 型件彎曲角度影響的有限元分析[J].熱加工工藝，2016，45（9）：162-164.

[15]王翠竹，牛衛中.彎曲角度與半徑對 Q345 鋼矩形管繞彎截面畸變的影響[J]．塑性工程學報，2021，28 （ 8） ： 112 -118.

[16]陳林，李曉謙.基于 MSC Marc 軟件平臺的鑄軋輥輥型三維有限元仿真[J]. 計算機輔助工程，2006，15（9）：321-325.

[17]陳麗，蘇建榕，陳鍇.基于MSC.Marc 的平板對接接頭殘余應力分布規律及影響因素[J].熱加工工藝，50（23）：135-145.