鋁鋰合金薄板微結構輥壓微成形缺陷及其電致塑性效應數值模擬研究

汪鑫偉，高鵬永，石勝龍，劉琛，徐杰，王春舉，李建偉，郭斌，單德彬

輕合金成形

鋁鋰合金薄板微結構輥壓微成形缺陷及其電致塑性效應數值模擬研究

汪鑫偉1a,1b，高鵬永1b，石勝龍1b，劉琛1a，徐杰1b，王春舉2，李建偉3，郭斌1a,1b，單德彬1b

（1.哈爾濱工業(yè)大學 a.空間環(huán)境與物質科學研究院 b.微系統(tǒng)與微結構制造教育部重點實驗室，哈爾濱 150001；2.蘇州大學 機電工程學院，江蘇 蘇州 215131； 3.北京航星機器制造有限公司，北京 100013）

研究微結構形狀及尺寸、壓下量、溫度、電流強度和電流路徑對電流輔助輥壓成形的影響規(guī)律，分析輥壓過程中多物理場耦合行為，探究電致塑性效應對輥壓成形的作用機理，并以此來解決薄板微結構形性協(xié)同控制關鍵問題。采用ABAQUS建立微結構輥壓有限元模型，通過輥壓成形數值模擬方法進行相關研究。首先設計了三因素三水平正交試驗，探究了溝槽形狀、溝槽深度、壓下量對成形效果的影響。其次通過改變電流大小、對比輥對輥（R2R）和輥對板（R2P）2種電流路徑、對比電流輔助輥壓和等溫輥壓模擬試驗來依次探究電流強度、電流路徑、電流的非熱效應對鋁鋰合金薄板微結構輥壓成形過程中微結構填充率與板材翹曲高度等輥壓成形缺陷的影響規(guī)律。在R2R電流路徑下，電流分布于整個板厚空間，而在R2P電流路徑下，電流主要分布于材料表面，2種路徑下的溫度與應力分布相似，但R2R下的峰值溫度略高，約高20 ℃；峰值應力略小，約小35 MPa。當壓下量與微溝槽寬度越大、溝槽形狀越接近矩形時，輥壓填充高度和翹曲高度均越大。在相同溫度下，與等溫輥壓相比，電流輔助輥壓微結構填充高度更高。三角形微溝槽在100 A、30%壓下量時的輥壓填充高度提升率高達7%。在等溫輥壓下，微結構填充時的塑性流動主要分布于表層；而在電流輥壓下，微結構填充時的塑性流動在整個板厚范圍內分布得較為均勻，且應力集中不明顯。溫度對輥壓微結構薄板翹曲程度的改善作用較小，電致塑性效應大幅改善了微結構薄板輥壓翹曲程度，翹曲高度減小約70%。從焦耳熱均勻性角度來看，電流路徑選擇R2R更優(yōu)。通過改變壓下量、微溝槽寬度與形狀可以改變微結構輥壓成形效果。相較于等溫輥壓，電流輔助輥壓可以提高輥壓填充高度、減小板材翹曲高度，故其成形效果更好，說明電流非熱效應可以促進輥壓成形的進行。

鋁鋰合金；電流輔助；成形缺陷；微結構；電致塑性效應

近年來，越來越多的學者開始關注微結構輥壓成形工藝，該工藝結合了傳統(tǒng)的微壓印工藝和輥壓成形工藝，能夠高效、低成本地實現大面積微結構制備[1]。微結構輥壓成形主要有2種類型，即輥對輥與輥對平板成形工藝[2]。在輥對輥輥壓成形方面，Hirt等[3-4]采用70～400 μm的鋼絲纏繞工藝制備了表面陣列微結構壓輥，并對寬20 mm、厚2 mm的Al 99.5和AlMg3板材進行了輥壓成形。研究發(fā)現，當板材厚度減小55%時，可以輥壓出高質量的半圓形陣列微溝槽結構，該微結構可用于飛行器壁面減阻，但局部鋼絲斷裂損壞、鋼絲松動等問題會影響微溝槽的輥壓成形質量。Zhou等[5]開發(fā)了桌面型輥對輥輥壓成形裝置，其中上壓輥表面加工微溝槽陣列可用來輥壓500 μm厚的AA5052鋁合金薄板。其團隊在薄板表面成形出寬100 μm、深20 μm的微溝槽陣列，同時發(fā)現上、下輥的相對速度對微溝槽的質量有顯著影響。在輥對平板輥壓成形方面，Gao等[6]通過純銅表面的微溝槽成形實驗，分析了輥表面微結構特征尺寸如槽寬（0.4～ 0.8 mm）、槽間距（0.8～1.2 mm）、圓角半徑（0.1～ 0.3 mm）與材料晶粒尺寸對成形質量的影響，并從微溝槽尺寸與晶粒尺寸交互影響的角度，揭示了微結構輥壓成形機理。研究發(fā)現，槽寬尺寸越小，成形越困難，微結構輪廓的不規(guī)則程度以及表面粗糙度隨著晶粒尺寸和輥壓深度的增大而增大；當微溝槽結構尺寸接近甚至達到金屬材料晶粒尺寸數量級時，由于參與變形的各個晶粒具備各向異性，局部塑性流動表現出不均勻性，即材料在塑性變形過程中產生尺寸效應[7-8]。尺寸效應對材料塑性變形的影響在過去的20年里獲得了廣泛的共識，該效應成為微成形領域的研究難點和熱點。總體上，尺寸效應不利于微結構的塑性成形，主要表現在它會降低成形性、增大工藝不穩(wěn)定性、增加成形缺陷、降低成形精度以及加劇摩擦磨損等[9-10]。同時，微結構的成形通常要求材料在微小模具型腔中完全填充，由于受微小模具型腔約束、接觸界面摩擦等影響，微結構的塑性成形行為更加復雜、變形抗力更大[11-12]。因此，在大面積陣列微結構輥壓成形工藝中引入新的物理場，從而提高成形溫度、減輕尺寸效應、改善材料成形極限是目前最有效的解決途徑[13-14]。

鋁鋰合金材料是近年來航空航天材料中發(fā)展非常迅速的一種先進輕量化結構材料[15-16]，具有密度低、彈性模量高、比強度和比剛度高、疲勞性能好、耐腐蝕及焊接性能好等諸多優(yōu)異性能[17]。用它代替常規(guī)的高強度鋁合金可使結構質量減輕10%～20%，剛度提高15%～20%，燃油節(jié)約5.4%，飛行費用降低2.1%[18]。因此，近年來，鋁鋰合金在航空航天領域的輕量化進程中扮演著越來越重要的角色，并成為最熱門的結構材料[17-19]。然而，由于鋁鋰合金具備室溫塑性差、屈強比高、各向異性明顯[20]、冷加工時容易開裂等缺點，其成形難度較大。目前鋁鋰合金只能成形較簡單的零件，難以制造復雜的零部件及微結構，因此其應用受到一定程度的限制。

近年來，相關學者通過實驗證明了電致塑性效應[21-22]，即當電流作用于微成形坯料后，電流與材料的相互作用可以抑制尺寸效應的發(fā)生并提高材料的塑性成形性能，具體表現為在材料成形過程中成形力降低、塑性提高、回彈減小等[23]。此外，電流還會產生非熱的電致塑性效應[24]，可以促進材料流動。電流輔助塑性微成形技術具有可實現極端尺寸和優(yōu)異性能微型構件或微結構批量制造的潛力[25-26]，在難變形鋁鋰合金薄板表面復雜微結構的成形及應用方面具有重要的研究前景。因此，本文建立了多參數交互作用下的鋁鋰合金薄板微結構電流輔助輥壓有限元數值模型，分析了不同參數下的電流輔助輥壓過程多物理場耦合行為，闡明了電致塑性效應對微結構填充率與板材翹曲等輥壓成形缺陷影響的作用機理，以期對解決電流輔助輥壓薄板微結構形性協(xié)同控制關鍵問題提供參考。

1 研究方案

1.1 材料及其性能

本文選用的材料為由西南鋁業(yè)生產的O態(tài)2195鋁鋰合金板，板厚為2 mm，其化學成分如表1所示。

利用線切割技術切取10 mm×10 mm×2 mm的試樣，經水砂紙打磨和拋光布機械拋光后，采用HF+HNO3+H2O（體積比為1︰4︰20）的混合溶液腐蝕3～5 s，用酒精沖洗、吹干后，獲得2195鋁鋰合金薄板①、②、③面的金相組織，3個面的具體位置如圖1a所示。采用OLYMPUS-SOL3550光學顯微鏡對鋁鋰合金薄板原始組織進行表征觀察，結果如圖1b～d所示。觀察發(fā)現，①面（見圖1b）的晶粒較為粗大，呈等軸狀，無明顯的拉長現象，而②面（見圖1c）和③面（見圖1d）的晶粒被拉長，呈纖維狀分布，尺寸最大達到數百微米。原始的2195鋁鋰合金板金相呈現出典型的軋制態(tài)特點，晶界及晶粒內部彌散分布著大量的強化相。

表1 2195鋁鋰合金成分

Tab.1 Composition of 2195 Al-Li alloy wt.%

圖1 2195鋁鋰合金原始板材金相組織

采用島津AGXplus萬能電子試驗機對2195鋁鋰合金分別在常溫、100 ℃、200 ℃和300 ℃下進行單向拉伸實驗。拉伸試樣尺寸如圖2所示，標距長度為5 mm，拉伸應變率為10?3s?1，每組拉伸實驗重復3次以減小測試誤差。

圖2 鋁鋰合金單向拉伸試樣尺寸

室溫下，2195鋁鋰合金試樣拉伸方向與軋向呈0o、45o和90o時的真實應力-應變曲線如圖3a所示?？梢钥闯觯牧系母飨虍愋圆⒉伙@著。常溫、100 ℃、200 ℃和300 ℃下鋁鋰合金軋向試樣的拉伸應力-應變曲線如圖3b所示?？梢钥吹?，當溫度在200 ℃以下時，材料呈現出加工硬化變形特點，說明此時溫度升高導致的位錯湮滅速率還不足以抵消塑性變形導致的位錯增殖速率。當溫度達到300 ℃時，材料呈現出明顯的軟化特征，說明此時材料發(fā)生了動態(tài)再結晶，材料的延伸率顯著增大，塑性變形能力大大提高。不同拉伸條件下材料的屈服強度、抗拉強度和延伸率等性能指標如表2所示。

1.2 研究方案

微結構輥壓成形數值模擬所要研究的變量及其參數如表3所示。首先設計了不同溝槽形狀、溝槽深度、壓下量的三因素三水平正交試驗進行常規(guī)輥壓成形數值模擬研究，選取缺陷顯著的參數組進行電流輔助輥壓成形與等溫輥壓成形的仿真研究。其次，為了研究電流路徑的影響，分別對微結構輥壓過程中輥對輥（R2R）和輥對板（R2P）[26]2種電流路徑下的多物理場耦合行為進行仿真分析。最后，為了研究電流非熱效應對鋁鋰合金薄板微結構輥壓成形的影響規(guī)律，設計了微結構等溫輥壓數值模擬方案，其中將等溫輥壓的溫度設置為接近電流輔助輥壓時的焦耳熱溫度?？紤]到2195鋁鋰合金在300 ℃時已發(fā)生再結晶軟化（見圖3），本文分別在100 ℃、200 ℃溫度下進行了微結構的等溫輥壓模擬。電流輔助輥壓時的電流強度參數分別為60、80、100 A，對應的平均焦耳熱溫度分別為140、200、240 ℃。

圖3 2195鋁鋰合金不同狀態(tài)下的真實應力-應變曲線

表2 2195鋁鋰合金拉伸力學性能

Tab.2 Tensile mechanical properties of 2195 Al-Li alloy

表3 微結構輥壓數值模擬參數表

Tab.3 Numerical simulation parameters of microstructure roll forming

2 有限元模型

采用ABAQUS建立電流輔助輥壓有限元模型，主要步驟包括幾何創(chuàng)建、網格劃分、材料屬性定義、相互作用與邊界條件設置、分析步創(chuàng)建等。

1）幾何創(chuàng)建與網格劃分。分別對電流輔助和非電流輔助輥壓成形進行數值模擬分析，其中后者是前者的對照組。薄板微結構輥壓成形有限元模型及電流路徑示意圖如圖4所示。由圖4a和圖4b可知，模型分為微結構輥、平輥以及板料三部分，板料厚度為0.5 mm，平輥直徑與微結構輥直徑相同，均為40 mm。對于非電流輔助輥壓模型（見圖4a），由于不涉及傳熱和導電等因素，所以為節(jié)約計算時間和提高模型收斂性，采用動力學顯式算法，且將軋輥均設定為剛體。對于電流輔助輥壓模型（見圖4b），由于需要考慮導電、傳熱、變形等多物理場耦合因素，所以采用隱式算法。采用離散剛體網格對圖4a中的軋輥進行劃分，而板料網格單元類型選為C3D8R；采用考慮了熱-電-結構耦合的Q3D8R單元對圖4b中的板料、平輥、微結構輥進行網格劃分，如圖4b中的放大圖所示，對微結構齒進行網格加密處理。2種不同的電流路徑分別如圖4c和圖4d所示。

2）材料屬性定義。在本文的有限元模型中，針對板材和軋輥采用的熱-電-結構多物理性能參數有溫度（）、密度（）、彈性模量（）、熱導率（）、電導率（）、比熱容（）、熱膨脹系數（），其數值分別如表4、表5所示。鋁鋰合金板材的塑性本構關系采用圖3中的試驗數據。

3）邊界與接觸定義。板材與軋輥接觸面間相互作用的定義包括摩擦、接觸電導率、接觸熱導率以及與環(huán)境間的輻射換熱系數，它們之間的接觸方式為面面接觸，接觸模型為罰函數摩擦模型，摩擦因數為0.2。接觸面的接觸電導率和熱導率分別為33 172 Ω?1·mm?1和44 W/(m·K)，與環(huán)境間的輻射換熱系數為1.134×10?9W/(m2·K4)。由于微結構尺寸遠遠小于板料的寬度，因此，為減小計算量，對板料進行對稱邊界建模，在對稱面上施加平面應變約束。

圖4 薄板微結構輥壓成形有限元模型及電流路徑示意圖

表4 鋁鋰合金熱-電-結構多物理性能參數

表5 Cr2Mo冷軋輥用鋼熱-電-結構多物理性能參數

Tab.5 Thermal-electrical-structural multi-physical performance parameters of Cr2Mo cold roll steel

3 結果與分析

3.1 電流輔助輥壓過程多物理場耦合行為

在相同電流強度（100 A）和相同時刻（0.04 s）下，R2R和R2P 2種電流加載路徑下2195鋁鋰合金薄板微結構的多物理場如圖5所示。可以發(fā)現，R2R與R2P的電流密度場有很大不同，R2R的電流密度主要集中分布在輥、板接觸區(qū)域，而R2P的電流密度則主要從輥、板接觸區(qū)向板材接電端呈梯度分布，且R2P的電流密度峰值比R2R的更大，峰值區(qū)域則更小。不同于電流密度場，不同電流路徑下溫度場的分布規(guī)律較相似，均是輥、板接觸區(qū)的溫度最高，且該區(qū)與電流密度集中分布區(qū)完全重合。該區(qū)域由于電流密度集中而形成了局部的焦耳熱熱源區(qū)，而其他區(qū)域則是由焦耳熱熱源區(qū)溫度沿板材軋向呈梯度遞減擴散而形成的，但在整個輥壓過程中，R2R的溫度略高于R2P的，約高了25 ℃。由溫度云圖還可以看出，R2R的溫度峰值區(qū)域占比更大且溫度梯度更均勻，因此從焦耳熱均勻性角度考慮，R2R更優(yōu)。由圖5c可知，R2R和R2P的應力分布基本一致，輥、板接觸區(qū)因受微結構齒擠壓而變形劇烈，應力峰值較大且較集中，而且R2R的應力峰值明顯低于R2P的，這主要是因為在相同壓下量下，在輥壓過程中，R2R的溫度始終高于R2P的。

在100 A電流及40%壓下量下，輥壓時R2R與R2P鋁鋰合金板截面上的等效應變場（PEEQ）、溫度場（NT）及電流密度場（ECD）如圖6所示。可知，峰值電流密度均位于微結構齒形接觸區(qū)，電流密度以沿厚度方向梯度遞減的方式貫穿分布于板材截面上，但R2R截面上的電流密度分布得更加集中，峰值電流密度區(qū)更小，且R2P的峰值電流密度要高于R2R的。對比圖6b和圖6e可以發(fā)現，在2種加載路徑下，溫度場的分布無明顯差別，但R2R的峰值溫度更高，這可能是因為R2R在微結構齒形接觸區(qū)的電流密度更集中。此外，還可以看出，輥壓變形區(qū)的金屬流動規(guī)律一致，兩者的等效應變場分布也無明顯差別，唯一的區(qū)別是R2R的最大應變略高于R2P的，這與兩者在峰值溫度上的差別相對應。

鋁鋰合金薄板微結構R2P輥壓成形的峰值溫度與電流強度的關系如圖7所示。可知，峰值溫度隨著電流強度的增大而增大，且在60～100 A下，焦耳熱溫度為100～300 ℃，低于鋁鋰合金形成大范圍動態(tài)再結晶時的溫度（見圖3）。還可以發(fā)現，峰值溫度隨時間的延長而非線性增大，即在初始0.03 s間隔內，溫度大幅上升，而在后續(xù)的0.03 s間隔內，溫度上升幅度下降，這符合焦耳加熱特征。此外，模擬結果表明，在輥壓過程中，溫度高低及溫差大小會隨溝槽陣列數目的增加而減小，這可能是因為接觸面積增大導致散熱加快和電流密度降低。

3.2 微結構填充率及其電致塑性效應

溝槽形狀（三角形、橢圓形、矩形）、溝槽深度（100、200、300 μm）以及壓下量（30%、40%、50%）3個因素對鋁鋰合金薄板微結構填充率的影響規(guī)律如圖8所示。不同形狀微溝槽輥壓成形過程中的法向應變分量及反作用力如圖9所示。由圖8可知，隨著壓下量和溝槽深度的增大，微結構溝槽填充率明顯增大，且在相同壓下量和溝槽深度條件下，三角形、橢圓形、矩形的溝槽填充率依次增大。這可以通過圖9中板材法向塑性應變分量（PE33）的差別進行解釋：在相同壓下量及溝槽深度下，三角形溝槽底部的法向壓應變最小，而矩形溝槽底部的壓應變最大；底部材料在擠壓作用下沿著微溝槽側壁向微結構型腔流動，產生法向向上的拉應變，因此，沿填充方向的正應變按照三角形、橢圓形、矩形的順序依次增大，說明材料在三角形溝槽中受到沿溝槽側壁的流動阻力最大，而在矩形中的最小，故填充率遵循圖8中的規(guī)律。3種形狀溝槽在輥壓成形時軋輥上的反作用力如圖9d所示?？梢钥闯?，在獲得相同深度的微結構時，三角形軋輥受到的力最大，因此相應的材料流動阻力最大。

圖5 不同電流路徑下薄板微結構輥壓成形的多物理場

圖6 不同電流路徑下薄板微結構輥壓成形截面上的多物理場

圖7 電流輔助輥壓過程峰值溫度與電流強度的關系

深度為200 μm的三角形微溝槽在等溫輥壓與電流輔助輥壓下的微結構填充輪廓圖及填充率與壓下量之間的關系如圖10所示。當壓下量為40%時，不同輥壓條件下薄板微結構法向塑性應變分量PE33分布云圖如圖11所示。由于微溝槽填充輪廓相差不大，僅選取中間2個溝槽的填充輪廓進行對比。觀察圖10可知，隨著溫度的升高，填充高度也在增大，表明溫度的升高導致鋁鋰合金薄板塑性成形性能提高，從而改善了微溝槽的填充性能。在相同壓下量下，100 A電流輔助輥壓下的三角形微溝槽填充高度大于等溫輥壓的結果，在壓下量為30%時最明顯，提升率約為8.8%，但隨著壓下量的增大，電流輔助導致的填充率增大量逐漸減小，這可能與在微結構輥壓成形填充極限下摩擦阻力顯著增大有關。不同溫度及電流強度下200 μm三角形微溝槽的填充率與壓下量之間的關系如圖10d所示?？梢园l(fā)現，在20、100、200 ℃的等溫條件下，微溝槽輥壓填充率隨著壓下量的增大而增大，但溫度越高，填充率隨壓下量增大的幅度越不明顯。此外，在大壓下量下，填充率隨溫度的變化也同樣遵循上述規(guī)律，該規(guī)律進一步說明，在填充極限下，摩擦阻力的作用主導了材料的填充行為。由圖10d還可以看出，60 A與80 A電流輔助下的填充率分別高于100 ℃及200 ℃下等溫輥壓的結果，100 A時的填充率明顯高于其他條件下的填充率。

圖9 不同形狀微溝槽輥壓成形的法向應變分量及反作用力

圖10 等溫與電流輔助條件下三角形微溝槽（200 μm）輥壓成形填充輪廓及其填充率

圖11 相同壓下量（40%）下不同輥壓條件薄板微結構法向塑性應變分量PE33分布云圖

產生上述現象的原因是電致塑性效應中的非熱效應促進了材料在微結構型腔中的流動填充。對比觀察圖11中填充方向的塑性應變分量可知，在同一壓下量下，等溫輥壓和電流輔助輥壓在微結構填充區(qū)域的金屬塑性流動行為存在明顯差異。在等溫輥壓時，如圖10a～c所示，塑性變形最劇烈的位置是微結構齒側壁尖角處，此處材料壓應變最大，且壓應變呈對稱分布，導致相鄰的材料被反向擠壓，產生最大的拉應變，使材料沿填充方向流入微結構型腔。此外還可看出，沿填充方向的材料塑性流動主要發(fā)生在板材表面較小的深度范圍內，板材大部分底部區(qū)域的壓應變大小一致，這說明板材底部區(qū)域只起到傳遞變形的作用，故金屬塑性流動區(qū)域比較有限。在電流輔助輥壓下，如圖10d～f所示，雖然最大壓應變仍然處在微結構齒側壁尖角處，但最大壓應變區(qū)沒有呈對稱分布，且明顯向板厚的深度方向延伸，在溝槽下方的整個板厚范圍內幾乎都存在較大的壓應變。同時，溝槽相鄰處的材料由于受到溝槽材料的反向擠壓作用而產生了最小壓應變，促使材料沿填充方向流入微結構型腔，但與等溫輥壓不同的是，在板厚方向上很大范圍內都存在與溝槽相鄰區(qū)域相近的應變。因此，金屬材料的流動范圍幾乎是整個板材厚度范圍，這使電流輔助條件下的材料填充性能得以提高。

考慮到當電流強度參數為60、80、100 A時，電流輔助輥壓對應的平均焦耳熱溫度分別為140、200、240 ℃，故選擇同等溫度下的等溫輥壓填充高度作為參照數據，計算電流輔助輥壓下電致塑性效應導致的填充高度的增大率，結果如圖12所示。可知，在不同壓下量下，電致塑性效應的非熱作用促進了微結構的填充，在電流為100 A、壓下量為30%時該值高達7%。

圖12 相同溫度下等溫輥壓電流輔助導致的微結構填充高度增大率

3.3 板材翹曲及其電致塑性效應

在薄板微結構輥壓過程中，微結構輥面和平輥面的接觸差異導致板材不可避免地產生了翹曲變形缺陷。由于本研究中輥壓試樣的長度一致，所以可采用翹曲高度定量表征板材翹曲程度的大小。翹曲高度為板材輪廓線上最高點到輪廓線首尾連線的距離。對室溫輥壓過程進行了模擬，獲得了溝槽形狀、溝槽深度、壓下量3個因素對翹曲高度的影響規(guī)律，如圖13所示?？梢钥闯?，隨著壓下量和溝槽深度的增大，板材的翹曲高度顯著增大；在相同壓下量和溝槽深度下，三角形、橢圓形、矩形的翹曲高度依次增大。實際上，在輥壓過程中，微結構輥和平輥作用在板材上的摩擦力不同，進而會導致板材上下面的金屬流動速率不同，微溝槽深度、形狀、壓下量的變化均會導致金屬塑性流動中的摩擦力發(fā)生變化，進而導致板材局部區(qū)域金屬流動程度產生差異而影響翹曲高度，該設想可在圖14中得到證明。在同一壓下量下，三角形、橢圓形、矩形微結構溝槽的等效塑性應變PEEQ云圖如圖14a～c所示?？芍?，峰值應變均在溝槽底部圓角處，且兩側呈對稱分布，矩形和橢圓形的PEEQ峰值相差不大，但三角形的明顯低于前兩者的，證明材料在三角形微結構溝槽中的局部流動程度不如在橢圓形和矩形中的劇烈。因此，等效塑性應變的局部分布差異是導致翹曲高度不同的直接原因。圖14d為不同形狀微結構輥在相同壓下量及溝槽深度下的軋向反作用力，即摩擦力大小，可知，三角形微結構輥摩擦力明顯低于矩形和橢圓形的，這進一步證明了三角形溝槽輥下材料塑性變形程度最輕，翹曲高度最小。

不同輥壓條件下200 μm三角形微溝槽薄板翹曲輪廓與翹曲高度如圖15所示。由圖15a可以看到，溫度的變化對薄板等溫輥壓翹曲的影響較小，這在圖15b中也可以得到印證。觀察圖15b可知，與20 ℃的常溫輥壓相比，100 ℃和200 ℃下的翹曲高度增大比例非常小，當壓下量為30%、40%、50%時，翹曲高度增大比例分別為2.5%、1.1%、1.25%，故溫度對翹曲程度的影響幾乎可以忽略。然而，在電流輔助輥壓條件下，微結構薄板翹曲問題得到了明顯的改善，如圖15a所示，但R2P和R2R 2種電流路徑下的板材翹曲輪廓一致。觀察圖15b可知，壓下量的增大會導致翹曲高度增大，對比電流輔助輥壓與等溫輥壓下的板材翹曲高度可知，無論在何種壓下量下，電流輔助導致的翹曲高度都比等溫輥壓的小70%以上，這種大幅度的翹曲高度改善效果已無法通過電流導致的溫度效應來解釋。

選擇同等溫度下的等溫輥壓翹曲高度作為參照數據，計算電流輔助輥壓下電致塑性效應導致的翹曲高度減小率，其結果如圖16所示?？芍诓煌瑝合铝肯?，電致塑性效應的非熱作用改善了微結構薄板輥壓的翹曲程度，翹曲高度的減小幅度約為70%。在同一壓下量（40%）下，等溫與電流輔助時微結構薄板輥壓截面區(qū)的等效塑性應變PEEQ與Mises的應力分布如圖17所示。可以看到，與等溫輥壓相比，在電流輔助輥壓過程中，鋁鋰合金薄板微結構局部處金屬的塑性流動行為發(fā)生了很大改變，具體表現如下：1）雖然兩者最大的塑性變形區(qū)均在微結構溝槽底角處，但電流輔助條件下的最大變形區(qū)明顯大于等溫輥壓條件下的，塑性變形從微溝槽接觸面向板材深度方向延伸的距離遠遠大于等溫條件下的；2）電流作用下輥壓變形區(qū)的最大PEEQ明顯低于等溫條件下的，甚至同一深度處的塑性應變大小均小于等溫輥壓的，這說明電流輔助輥壓下的局部塑性變形更加均勻。為進一步分析導致上述現象的原因，對等溫輥壓與電流輔助輥壓板厚切面上的Mises應力云圖進行對比觀察，發(fā)現等溫輥壓下溝槽底角處出現了應力集中現象，由于溝槽處存在擠壓與反擠壓作用，所以應力集中區(qū)沿板材表面呈周期性對稱分布；而在電流輔助輥壓條件下，最大應力區(qū)不在溝槽底角處，溝槽齒形區(qū)與板厚左右兩側的應力最大，最小應力區(qū)位于板厚切面的中心并呈類X形對稱分布。因此，可以認為在電流輔助輥壓條件下，應力分布更加分散，應力集中不明顯，材料流動均勻性更好，故其翹曲程度更小。

圖13 微結構薄板翹曲高度與溝槽形狀、深度及壓下量之間的關系

圖14 不同形狀微溝槽輥壓成形的等效塑性應變云圖及摩擦力

圖16 相同溫度下等溫輥壓電流輔助導致的薄板翹曲高度減小率

4 結論

采用數值模擬研究了2195鋁鋰合金薄板微結構電流輔助輥壓成形效果，分析了電流輔助輥壓過程中的多物理場耦合行為及薄板微結構輥壓成形缺陷影響規(guī)律，研究了電流路徑、電流強度、微溝槽深度、微溝槽形狀、壓下量等因素對薄板微結構電流輔助輥壓微成形缺陷的影響規(guī)律，闡明了微結構填充率、板材翹曲等輥壓成形缺陷的電致塑性效應作用機理。主要結論如下：

1）在R2R電流路徑下，電流分布于整個板厚空間，而在R2P電流路徑下，電流主要分布于材料表面，2種路徑下的溫度與應力分布相似，但R2R的峰值溫度略高，約高20 ℃；峰值應力略小，約小35 MPa。

2）當壓下量與微溝槽寬度越大、溝槽形狀越接近矩形時，輥壓填充高度和翹曲高度越大。

3）在相同溫度下，電流輔助輥壓微結構填充高度大于等溫輥壓的，三角形微溝槽在100 A與30%壓下量時的輥壓填充高度提升率高達7%。

4）溫度對輥壓微結構薄板翹曲程度的改善作用較小，電致塑性效應大幅改善了微結構薄板輥壓翹曲程度，翹曲高度減小約70%。

5）在等溫輥壓下，微結構填充時的塑性流動區(qū)域主要分布于表層；而在電流輥壓下，微結構填充時的塑性流動區(qū)域在整個板厚范圍內分布得較為均勻，且應力集中不明顯。

[1] LI J, GAO Z Y, HAN P D. Roll-to-plate Micro-imprinting of Large-area Surface Micro Features: Experiments and a Local Simulation Method[J]. Materials Today Communications, 2023, 34: 105011.

[2] 高緒杰, 王回航, 朱光明, 等. 減阻微溝槽成形工藝研究進展[J]. 塑性工程學報, 2021, 28(5): 183-191. GAO Xu-jie, WANG Hui-hang, ZHU Guang-ming, et al. Research Progress in Forming Process of Micro-grooves with Drag Reduction[J]. Journal of Plasticity Engineering, 2021, 28(5): 183-191.

[3] HIRT G, THOME M. Large Area Rolling of Functional Metallic Micro Structures[J]. Production Engineering, 2007, 1(4): 351-356.

[4] HIRT G, THOME M. Rolling of Functional Metallic Surface Structures[J]. CIRP Annals, 2008, 57(1): 317-320.

[5] ZHOU R, CAO J, EHMANN K, et al. An Investigation On Deformation-based Surface Texturing[J]. Journal of Manufacturing Science and Engineering, 2011, 133(6): 061017.

[6] GAO Z Y, PENG L F, YI P Y, et al. Grain and Geometry Size Effects on Plastic Deformation in Roll-to-plate Micro/meso-imprinting Process[J]. Journal of Materials Processing Technology, 2015, 219: 28-41.

[7] GU X, WANG X, MA Y J, et al. Investigation on Grain Size Effect and Forming Mechanism of Laser Shock Hydraulic Microforming of Copper Foil[J]. The International Journal of Advanced Manufacturing Technology, 2021, 114(3): 1049-1064.

[8] XU Z T, PENG L F, LAI X M. Investigation on the Roll-to-plate Microforming of Riblet Features with the Consideration of Grain Size Effect[J]. The International Journal of Advanced Manufacturing Technology, 2020, 109(7): 2055-2064.

[9] VOLLERTSEN F, BIERMANN D, HANSEN H N, et al. Size Effects in Manufacturing of Metallic Components[J]. CIRP Annals, 2009, 58(2): 566-587.

[10] RAJA P C, RAMESH T. Influence of Size Effects and Its Key Issues During Microforming and Its Associated Processes-a Review[J]. Engineering Science and Technology, 2021, 24(2): 556-570.

[11] WANG C J, WANG C J, XU J, et al. Interactive Effect of Microstructure and Cavity Dimension on Filling Behavior in Micro Coining of Pure Nickel[J]. Scientific Reports, 2016, 6(1): 23895.

[12] WANG H H, GAO X J, ZHU G M, et al. Process Analysis and Hole Type Optimization of Micro-groove Multi- pass Rolling[J]. The International Journal of Advanced Manufacturing Technology, 2021, 119(3): 2201-2212.

[13] 孟寶, 潘豐, 杜默, 等. 陣列微流道電場輔助輥壓成形技術研究[J]. 機械工程學報, 2022, 58(20): 231-241. MENG Bao, PAN Feng, DU Mo, et al. Research on Arrayed Micro-channel Electric Field Assisted Roll Forming Technology[J]. Journal of Mechanical Engineering, 2022, 58(20): 231-241.

[14] ZHENG G Y, LUO X, KOU Z D, et al. Microstructural Evolution of Al-Zn-Mg-Cu Alloy during Ultrasonic Surface Rolling Process[J]. Materials Characterization, 2022, 194: 112418.

[15] 黃曉敏, 管奔, 臧勇. 航空鋁鋰合金熱成形研究進展[J]. 稀有金屬材料與工程, 2022, 51(12): 4745-4756. HUANG Xiao-ming, GUAN Ben, ZANG Yong. Research Progress on Thermal Formability of Al-Li Alloys for Aeronautic Industry[J]. Rare Metal Materials and Engineering, 2022, 51(12): 4745-4756.

[16] ARBI E H, KENZA B, MUSTAPHA F, et al. A Review of Manufacturing Processes, Mechanical Properties and Precipitations for Aluminum Lithium Alloys Used in Aeronautic Applications[J]. Heliyon, 2023, 9(3): e12565.

[17] HUANG X M, GUAN B, WANG B Y, et al. Hot Tensile Deformation Behavior, Fracture Mechanism and Microstructural Evolution of 2195 Al-Li Alloy[J]. Journal of Alloys and Compounds, 2023, 946: 169426.

[18] 謝磊. 2A97鋁鋰合金晶粒細化及超塑性研究[D]. 長沙: 中南大學, 2013: 1-6. XIE Lei. Investigation of Grain Refinement and Superplasticity of 2A97 Al-Li Alloy[D]. Changsha: Central South University, 2013: 1-6.

[19] LEQUEU P, LASSINCE P, WARNER T. Aluminum Alloy Development for the Airbus A380-Part 2[J]. Advanced Materials and Processes, 2007, 165(7): 41-45.

[20] LIU T L, LI H R, MA Y L, et al. Investigation of Anisotropy and Structure Variation of Spray-formed 2195 Al-Li Alloy via Final Temperature-controlled Rolling and Cold Rolling[J]. Journal of Alloys and Compounds, 2023, 937: 168414.

[21] LI C Z, XU Z T, PENG L F, et al. An Electric-pulse-assisted Stamping Process towards Springback Suppression and Precision Fabrication of Micro Channels[J]. International Journal of Mechanical Sciences, 2022, 218: 107081.

[22] XU Z T, JIANG T H, HUANG J H, et al. Electroplasticity in Electrically-assisted Forming: Process Phenomena, Performances and Modelling[J]. International Journal of Machine Tools and Manufacture, 2022, 175: 103871.

[23] PERKINS T A, KRONENBERGER T J, ROTH J T. Metallic Forging Using Electrical Flow as an Alternative to Warm/Hot Working[J]. Journal of Manufacturing Science and Engineering-transactions of The Asme, 2007, 129: 84-94.

[24] HUANG T, XING B H, SONG K X, et al. Thermal and Non-thermal Effects of Cu/Al Laminated Composite during Electrically Assisted Tension[J]. Materials Science and Engineering: A, 2023, 878: 145237.

[25] 孫宏霖. 2024鋁合金陣列微結構電流輔助輥壓成形工藝研究[D]. 哈爾濱: 哈爾濱工業(yè)大學, 2022: 7-10. SUN Hong-lin. Research on Current-assisted Roll Forming Process of Microchannel Array of 2024 Aluminum[D]. Harbin: Harbin Institute of Technology, 2022: 7-10.

[26] 席瑞. 2195鋁鋰合金電流輔助熱沖壓工藝研究[D]. 哈爾濱: 哈爾濱工業(yè)大學, 2020: 6-12. XI Rui. Research on Electrically-assisted Hot Forming-die Quenching Integrated Process of 2195 Al-Li Alloy[D]. Harbin: Harbin Institute of Technology, 2020: 6-12.

Numerical Simulation on Electrically-assisted Roll Micro-forming Defects of Micro-structures and Electroplastic Effects on Al-Li Alloy Sheets

WANG Xin-wei1a,1b, GAO Peng-yong1b, SHI Sheng-long1b, LIU Chen1a, XU Jie1b, WANG Chun-ju2, LI Jian-wei3, GUO Bin1a,1b, SHAN De-bin1b

(1.a. Laboratory for Space Environment and Physical Sciences, b. Key Laboratory of Micro-Systems and Micro-Structures Manufacturing, Ministry of Education, Harbin Institute of Technology, Harbin 150001, China; 2. School of Mechanical and Electrical Engineering, Soochow University, Jiangsu Suzhou 215131, China; 3. Beijing Hangxing Machinery Manufacture Limited Corporation, Beijing 100013, China)

The work aims to study the effect of microstructure shape and size, reduction, temperature, current intensity and current path on electrically-assisted roll forming, analyze the coupling behavior of multiple physical fields in the roll forming process, and clarify the mechanism of electroplastic effect in roll forming, so as to solve the key problem of the synergistic control of the accuracy and property of micro-structures on sheets. ABAQUS was used to establish a finite element model, and the related research was carried out by numerical simulation of roll forming. Firstly, an orthogonal experiment with three factors and three levels was designed to explore the effect of groove shape, groove depth and reduction on forming. Then, the effect of current intensity, current path and non-thermal effect on microstructure rolling defects of Al-Li alloy sheet such as microstructure filling rate and sheet warping height was explored, by changing the magnitude of current, comparing roll-to-roll (R2R) connection and roll-to-plate (R2P) connection, and comparing electrically-assisted and isothermal roll forming simulation tests. The current density for the R2R was distributed across the whole thickness of the sheet, while it tended to be distributed on the surface in the R2P. The distributions of temperature and stress for the two connection modes were similar, but compared with R2P, R2R had a slightly higher peak temperature of about 20 ℃ and a slightly lower peak stress of about 35 MPa. When the roll forming reduction and the width of micro-groove were larger, and the groove shape was closer to rectangle, the groove filling height and the sheet warping height were larger. The filling height of microstructure by electrically-assisted roll forming at the same temperature was higher than that by isothermal case, and the filling height of triangular microgrooves could be increased by up to 7% at 100 A and 30% reduction. The plastic flow tended to be located at the surface layer in the isothermal roll forming of micro-structures, while the microstructure filling was uniform throughout the thickness direction, and the stress concentration was not obvious by electrically-assisted case. Temperature had little effect on the warping degree of the rolled sheet, while the electroplastic effect greatly reduced the warping height of the rolled sheet, causing about 70% reduction of the warping height. According to the above results, the following conclusions can be drawn. R2R is better current path from the view of Joule heating uniformity. The roll forming effect on microstructure can be changed by varying the forming reduction, the width and shape of micro-grooves. Compared with the isothermal case, the electrically-assisted roll forming can improve the filling height of micro-structures and reduce the warping height of sheets, which shows better formability since the non-thermal effect of current plays a promoting role in roll forming defect control.

Al-Li alloy; electrically-assisted; forming defects; micro-structures; electroplastic effect

10.3969/j.issn.1674-6457.2023.010.010

TG335.5+5

A

1674-6457(2023)010-0082-13

2023-07-31

2023-07-31

黑龍江省自然科學基金優(yōu)秀青年項目（YQ2020E014）

Heilongjiang Provincial Natural Science Foundation of China (YQ2020E014)

汪鑫偉, 高鵬永, 石勝龍, 等. 鋁鋰合金薄板微結構輥壓微成形缺陷及其電致塑性效應數值模擬研究[J]. 精密成形工程, 2023, 15(10): 82-94.

WANG Xin-wei, GAO Peng-yong, SHI Sheng-long, et al. Numerical Simulation on Electrically-assisted Roll Micro-forming Defects of Micro-structures and Electroplastic Effects on Al-Li Alloy Sheets[J]. Journal of Netshape Forming Engineering, 2023, 15(10): 82-94.

責任編輯：蔣紅晨