丁瑩

（北京北方華創新能源鋰電裝備技術有限公司　北京　110000）

The Application Simulation Analysis on Lithium-ion Battery Pole Piece Rolling

Ding　ying

(Naura New Energy Technology System Co., Ltd, Beijing 110000, China)

摘要：隨著國內外對電池極片軋制設備的高速、高精度發展需求，促使仿真技術研究極片軋制的課題越來越多。介紹了非對稱涂布的電池極片軋制仿真分析，仿真結果說明液壓伺服加壓式極片軋機能夠軋制非對稱涂布電池極片，克服了氣液增壓泵加壓式電池極片軋機對此類極片無法正常軋制的缺點。綜述了仿真技術在液壓領域和熱輥壓機方面的研究進展，并進行了總結及其所存在的不足。

引言

近年來新能源行業蓬勃發展，給鋰離子電池行業帶來了巨大的發展空間^[1,2]，同時對鋰離子電池的質量和壽命等也提出了更高的要求。在鋰離子電池的生產過程中，極片的制造完成決定了電池80％以上的性能^[3,4]。因此，正負極片質量的好壞對鋰離子電池的好壞有重大的影響，極片的制造在電池制造工藝中占有重要的地位。

國內外大量研究表明，正負極片的壓實密度對鋰離子電池的容量等電池性能有較大的影響^[5-10]，合適的壓實密度可以增大電池的容量，減少內阻，減少極化損失等。軋制后的電池極片密度的均勻性和厚度的一致性直接決定著電池使用壽命的長短和儲能的多少^[11]。可見，用軋機對極片進行軋制是生產極片過程中極其重要的一道工序，由于極片軋機對電池生產的重要性，也必然要求軋制設備向高精度，自動化發展^[12]。

仿真技術以數學模型的建立、驗證、實驗為核心，綜合計算機技術、自動控制及系統工程、信息處理技術等多學科^[13]。用計算機對設備系統進行仿真的一般過程如下：

（1）建立所研究的控制系統數學模型。

（2）通過仿真軟件把數學模型轉化為計算機仿真模型。

（3）采用合適的算法編寫仿真程序或者直接采用他人程序。

（4）通過仿真，得到系統的動態響應曲線和數據。

（5）分析系統仿真結果，提出提高系統動態性能的方法和措施。

隨著計算機技術和科技的不斷發展，人們可以利用仿真技術方便的分析軋機設備系統動態性能，確定最優參數，獲得最佳的設備控制系統。電池極片負載特性是非線性的，軋制過程負載也隨時間變化，電池極片軋機的伺服控制系統具有非線性和時變性，進而有必要對電池極片的軋制過程進行仿真研究，以了解各個因素對電池極片軋制質量和軋制效率的影響規律，對后續電池極片軋機的設計具有重要的理論指導意義^[14]。因此，仿真技術已成為鋰電軋機設備系統設計的重要輔助手段之一^[15]。

國內外通過仿真技術研究軋機的課題越來越多，非對稱涂布電池極片軋制仿真分析是應用仿真技術的典型課題，與此同時，提出了許多仿真技術在液壓領域的應用課題。此外，仿真技術在熱輥壓機中的應用課題也備受關注。因此，本文分別從以下三方面進行介紹及綜述。

1　非對稱涂布電池極片軋制仿真分析

電池極片的涂布方式有對稱涂布與非對稱涂布。因為實際使用氣液增壓泵加壓式電池極片軋機對非對稱涂布的電池極片的軋制效果差，由此，北京北方華創新能源鋰電裝備技術有限公司與王益群等人^[14]合作以該司軋機為主要研究對象，分別對氣液增壓泵加壓式電池極片軋機和液壓伺服控制加壓式極片軋機進行非對稱涂布電池極片軋制仿真分析。

1.1　氣液增壓泵加壓式電池極片軋機

氣液增壓泵加壓方式電池極片軋機^[16]成本較低，能夠軋制對稱涂布的電池極片，但是它采用楔鐵和絲杠離線調節輥縫，不能對軋輥間隙和軋制力進行實時在線調節，使其軋制速度、極片的質量及應用范圍受到一定程度的制約。

研究者們基于氣液增壓泵加壓式電池極片結構與軋制特性的分析，建立的電池極片模型及軋機軋制極片過程的仿真模型如圖1所示。由于實際使用氣液增壓泵加壓式電池極片軋機對非對稱涂布的電池極片（其示意圖如圖2所示）的軋制效果差，下面對此種情況進行仿真分析，仿真結果如圖3所示。

1-上軋輥；2-電池極片；3-下軋輥；4-楔鐵；5-液壓缸作用力

圖1　液壓缸調定壓力與輥縫預調節值匹配的規律

通過仿真結果可以看出在極片從雙層滾壓到單層時，作用在極片上的軋制力由346kN 減小到96kN，單層部分的密實度不能滿足要求。油缸壓力與輥縫間隙是根據雙層部分的軋制需要進行離線設定的，一旦設定，就不能改變。作用在單層上的軋制力是由軋輥和楔鐵的剛度決定的，是不可調節的。軋輥的剛度越大，作用在單層上的力就越小。在軋輥滾壓到極片間隙時，軋輥與極片脫開，作用在極片上的力為零。

圖2　非對稱涂布電池極片示意圖

1-作用在極片上的軋制力；2-消耗在楔鐵上的力；3-極片入口厚度；4-極片出軋出厚度

圖3　軋制極片仿真結果

1.2　液壓伺服控制加壓式極片軋機

液壓伺服控制加壓式極片軋機不再使用楔鐵調節輥縫值，液壓缸壓力能夠完全作用在電池極片上，為了能夠實時控制作用在電池極片上壓力和液壓缸活塞位置，加壓系統采用閥控缸的液壓伺服控制系統。由于沒有楔鐵，當使用位置環軋制時，可以實現負的預輥縫軋制，不再像氣液增壓泵加壓式極片軋機那樣使用塞尺檢測輥縫，只能設定正值輥縫，克服了由于楔鐵的作用而使有效軋制力不能不斷提高的缺點。

氣液增壓泵加壓式電池極片軋機無法軋制非對稱涂布電池極片，由于在做軋制實驗時沒有非對稱涂布的電池極片，研究人員對非對稱涂布情況進行恒軋制力的仿真分析，仿真結果如圖4所示。通過仿真結果可以看出，軋輥從雙側涂布漿料部分滾壓到單側涂布部分時，軋制力有一個減小的波動，但是能夠迅速恢復到設定的軋制力值，能夠保證軋制力基本恒定。說明可以采用本課題所設計的軋機軋制此類電池極片，克服了氣液增壓泵加壓式電池極片軋機的缺點。

伺服加壓式電池極片軋機壓力設定值可以在線實時調節，在軋制非對稱涂布的電池極片時，還可以根據工藝要求，使雙層涂布與單層涂布的地方使用不同的軋制力或不同的輥縫軋制，進一步加大了軋機的應用范圍。

圖4　恒軋制力（單側400kN）軋制非對稱涂布電池極片仿真結果

傳遞的力和功率大的液壓伺服控制系統的引入使得極片軋機能夠實現壓力和輥縫的在線實時調節，軋制單雙層交替涂布的極片時，單層部分也能得到比較好的軋制效果，使得軋制極片的質量大大提高。但是由于極片涂布形式的多樣，極片間隙處厚度突然變薄，極片負載的突變，還是影響了軋制速度的進一步提高。極片軋機具有非線性和時變性，電池極片負載也是非線性的，也都對軋機的設計與控制提出了挑戰。

2　仿真技術在液壓領域的應用

所謂液壓系統仿真就是建立所研究液壓系統的數學模型，并轉化為計算機上的仿真模型，然后進行求解運算來分析所研究系統的動靜態特性的過程。

液壓系統仿真不僅可以找出現有系統的缺陷所在，提出對系統的改進措施，對系統進行實質上的改進和修正；而且可以對系統設計過程中的某些參數進行檢驗分析，確保系統具有良好的動態性能。在實際項目中，所研究內容要在系統最脆弱的地方展開分析。這種情況下要忽略一些次要矛盾，突出主要矛盾，做到有的放矢，解決問題的關鍵。

液壓系統計算機仿真的關鍵和難點在于：一建立描述液壓系統的準確數學模型，二是合理的編寫仿真程序。一個準確的數學模型的建立是仿真的基礎和前提，目前建立數學模型的常用方法主要有傳遞函數法、解析建模法和功率鍵合圖法。

仿真技術的應用在液壓領域主要包括下面幾個方面

（1）設計液壓系統時，通過理論推導建立所設計液壓系統的數學模型，經過仿真，把仿真數據和實驗數據進行對比來驗證所建立數學模型的準確性，并把經過驗證的數學模型作為以后改進和設計類似系統的依據。

（2）利用建立的數學模型和仿真模型，通過仿真實驗來確定已知系統參數的調整范圍，從而縮短系統的調整時間，提高效率。

（3）通過仿真實驗，研究所建液壓系統的可行性，以及各個液壓元件參數對系統的影響，最后獲得最佳的元件參數和控制方案^[15]。

（4）利用仿真技術研究新設計元件各結構參數對系統動態特性的影響，確定元件合適的結構參數。

電池極片軋機液伺服控制系統的設計，主要解決傳統極片軋機依靠楔鐵絲杠調節輥縫所帶來的缺陷，實現恒輥縫和恒軋制力制，提高軋制速度。利用仿真驗證所設計液壓伺服系統的準確性，并研究提高極片軋機伺服系統性能的措施^[17,18]。

董敏等人^[19]為了全面研究軋機液壓厚控系統動態特性，建立直觀真實的虛擬軋機模型，提出基于AMESim和ADAMS聯合建模的仿真方法。在ADAMS中構造軋機實體剛柔耦合動力學模型，實現了軋機的負載特性研究；在AMESim中建立液壓系統物理模型，實現了液壓伺服系統精確建模和分析，兩者通過接口實現數據交換，保證了液壓系統模擬的準確性和負載系統模擬的真實性。通過聯合仿真模型得到了系統實時響應以及出口板厚實時數據，將模型仿真輸出數據與實測數據進行比較，證明仿真模型能準確體現系統動態響應，并能體現機械部件在載荷下彈性變形和板厚實時輸出情況。

Liu等人^[20]為了提高極片軋機軋制效率和軋板的質量，建立了液壓系統的動態仿真模型，給出負載特性方程，模擬了恒軋制力、恒輥縫等軋制過程，仿真與實驗結果基本一致。該課題組以北京北方華創新能源鋰電裝備技術有限公司液壓伺服加壓式電池極片軋機為依托，分析了液壓伺服加壓式極片軋機液壓系統，建立液壓泵、蓄能器和減壓閥等的動態仿真模型，對模型進行仿真，并作了相應的實驗，對元件模型進行了實驗驗證。在元件模型的基礎上對整個液壓系統進行建模，使得建立的系統模型考慮了液壓泵和蓄能器對伺服閥前壓力的影響，減壓閥對伺服缸背壓的影響。在模型上對閉環軋制力階躍響應和軋制過程進行動態仿真分析，并在所研究極片軋機上作了相應的實驗，仿真結果與實驗結果基本一致，驗證了模型的正確性，所建立的仿真模型具有實際參考價值，為液壓控制系統參數優化提供了重要依據^[14]。

3　仿真技術在熱輥壓機中的應用

熱輥壓機是指能加熱軋輥的輥壓機，當利用熱輥壓機輥壓極片時，可提高極片的壓實密度、質量和生產速度。運用不同的加熱工藝對軋輥進行加熱的加熱過程、輥面溫度分布均勻性和能耗進行分析，并在此基礎上對不同結構軋輥的加熱過程、輥面溫度分布均勻性和熱應力進行分析，對指導現場加熱工藝，改進軋輥結構，提高極片軋制質量有著重大的意義。國內外對軋輥的研究主要分為兩個方面：1）加熱軋輥的方式及其相應的控制系統，對軋輥的加熱方式主要包括電加熱和油加熱；2）軋輥的溫度分布和內部應力等。

王文成等人研究了感應加熱下圓柱體的電磁場和溫度場，指出電磁加熱具備精度高，能耗低等特點，并開發了軋輥感應加熱系統^[21,22]；李徽通過實際生產論述了邊部加熱器在熱軋寬帶鋼生產中的應用，通過對軋輥邊部進行加熱，有利于提高熱軋鋼板凸度精度^[23]。

申世杰等人利用有限元軟件對萬能軋機進行了溫度分析、熱力耦合分析和疲勞壽命研究，指出輥頸與輥身之間的圓角處是熱應力集中的位置，此處最危險，增大輥頸的直徑和倒角，可使應力集中得到明顯改善^[24]；胡仕成等人建立了軋輥輥套和鑄壞之間的接觸導熱模型和傳熱數學模型，對軋輥和鑄坯進行了仿真分析，軋輥輥套的粗糙度越低，其導熱能力越強^[25]；王永洲等人對軋機軋輥進行了應力分析和模態分析^[26]；羅麗萍等人通過有限元軟件對軋輥表面電渣加熱和渣池熱電場進行了模擬研究，分析渣池內的溫度分布^[27]。

王興東等人^[28]運用傳熱學知識及有限元軟件進行仿真分析的方法，通過更改橫向油道與輥面之間的距離δ，按照同種加熱工藝對三種不同結構的軋輥進行加熱。給出了軋輥內部應變圖如圖5所示，δ=70mm輥面的應變為0.53mm/m，δ=50mm輥面的應變為0.61mm/m，說明δ越小，輥面的變形越大；并利用有限元軟件對軋輥進行溫度和應力分析。分析結果表明δ值越小，加熱時間越短，輥面的溫差越小；三種不同結構軋輥的最大應力和應變差別不大；δ越小，油道處的應力越大，為改進軋輥內部油道結構提供理論依據。

圖5　δ=70mm（左）與δ=50mm（右）軋輥內部應變圖

此外，通過軋輥輥面及內部溫度圖6可知，軋輥進油端的溫度高于另一端；兩端軸頸表面溫度達到了100℃以上，應選擇耐高溫的軸承；軋輥輥身溫度分布較為均勻，不存在較大的溫差，芯部的溫度高于輥面。該研究還指出，軋輥相當于一個集熱容器，加熱時間越長，軋輥對導熱油的吸熱量越少，進出口處油溫差越小，導熱油損失的能量越少，軋輥內部溫度越均勻，軋輥輥面溫度分布越均勻。

圖6　軋輥輥面溫度（左）及內部溫度（右）分布云圖

陳國樑等人運用有限元軟件建立了油加熱軋輥模型，對熱軋輥三維溫度場進行分析計算，分析結果表明：在軋輥內部，溫度的梯度變化主要是徑向的變化，但還有一些溫度梯度沿軸向變化，提高軋輥端部溫度均勻性可改善軋輥溫度分布均勻性^[29,30]。李立新等人通過測量工作輥表面溫度值，利用有限元軟件建立CSP工作輥溫度場模型，研究了軋制過程中的軋輥的溫度場^[31]。蔡輝發明了一種主要用于加熱軋輥軋輥溫度控制系統，將導熱油控制在一定的溫度范圍內，同時可以向軋輥內部導入常溫導熱油，導熱油被密封在系統中，不結觸鎂合金板帶，不對其產生仍何影響，導熱油也不會污染環境，從而使軋輥的工作溫度控制在合適的范圍，提高了產品成品率和生產效率，軋制出產品質量良好^[32]。

Denis等利用有限元軟件對熱輥壓機的軋輥進行了熱應力分析和疲勞壽命分析^[33]。Luks等人通過實驗的方法分析了軋輥之間的接觸應力和輥面溫度^[34]。

綜上所述，目前對于軋輥的溫度場、應力場和使用壽命研究得比較多，有必要進一步對油加熱的軋輥如軋輥輥面溫度分布、加熱過程和能耗、應力應變等進行深入的分析，如探討油道數量對軋輥的輥面溫差、加熱過程及應力等的影響。

4　結語

經過實驗驗證了的電池極片軋機模型上，進行了非對稱涂布的電池極片的虛擬軋制，仿真結果說明了氣液增壓泵加壓式電池極片軋機無法正常軋制非對稱涂布的電池極片，而液壓伺服加壓式電池極片軋機能夠軋制類似的電池極片。仿真技術的應用為液壓控制策略的改進提供了有效的模型平臺。對軋輥的加熱過程和應力等有待進行深入的仿真分析。

在仿真分析過程中，改變模型參數，可以模擬實際物理系統中各個元件，機械部件等參數的變化，從而可以了解這些實際參數的變化對整個系統的影響，便于尋找影響系統性能的主要因素。這些都非常利于設備的改進設計，提高設備的性能，縮短設計周期，調高效率。對于特別的控制系統，還可以在仿真模型上進行控制算法的仿真實驗，便于探索更好的提高系統性能的控制算法，減小進行多次物理實驗的時間和財力投入。

在理論計算和仿真分析過程中，采用了很多假設條件，運用實驗關聯式求解，存在較大誤差，因此還需進一步修正仿真分析中的邊界條件。

參考文獻

[1]郭永興. 鋰離子動力電池制造關鍵技術基礎及安全性研究[D]. 長沙：中南大學學位論文, 2010：9-10.

[2]胡喜鵬. 電池極片輥軋設備廠商迎機遇[N]. 中國能源報, 2010-06-14(2).

[3]楊洪, 何顯峰, 李峰. 壓實密度對高倍率鋰電池性能的影響[J]. 電源技術, 2009, 33(11): 959-962.

[4]Zheng H H, Li J, Song X Y. A comprehensive understanding of electrode thickness effects on the electrochemical performances of Li-ion battery cathodes[J]. Electrochimica Acta, 2012, 71(2): 258-261.

[5]李華成, 李普良, 盧道煥, 等. 粒度分布對錳酸鋰極片壓實密度的影響[J]. 礦業工程, 2011, 31(5): 94-96.

[6]王麗君. 壓實密度及電解液種類對鋰離子電池的影響[J]. 廣州化工, 2013, 41(15): 76-78.

[7]Deiss E, Wokaun A. Average voltage, energy density, and specific energy of lithium-ion batteries[J]. Journal of the Electrochemical Society, 1997, 144(11): 3877-3881. 

[8]Zu C X, Li H. Thermodynamic on analysis on energy densities of batteries[J]. Energy & Environmental Science, 2011, 4(8): 2614-2624.

[9]Oyama N, Hatozaki O. Lithium polymer battery with high energy density[J]. Macromolecular Symposia, 2000, 156(45): 171-178.

[10]Helene K, Wolfgang B, Jean-Pierre B. Industrial awareness of Lithium batteries in the word during the past two years[J]. Journal of Power Sources, 1998, 3: 43-50.  

[11]程建軍. 國內鋰離子動力電池關鍵工藝裝備[J]. 科技傳播, 2012, 7(5): 82-83.

[12]趙東林. 國內外鋰離子動力電池極片制造專用設備的新進展[J]. 新材料產業, 2009(9): 67-68.

[13]吳重光. 仿真技術[M]. 北京: 化學工業出版社, 2000: 9-11.

[14]袁叢林. Φ800電池極片軋機液壓控制系統建模仿真及實驗研究[D]. 西安：燕山大學學位論文, 2013.

[15]李成功, 彥淼. 液壓系統建模與仿真分析[M]. 北京: 航空工業出版社, 2008: 1-30.

[16]王建斌. 配備自動厚度調節的電池極片軋制設備[J]. 電源技術研究與設計, 2008, 32(11): 790-792.

[17]韓虎, 劉印鋒, 孫成通, 等. 基于 MATLAB 液壓系統的仿真技術研究[M]. 液壓氣動與密封, 2007: 20-22.

[18]王勇, 張勇, 李從心. 液壓仿真軟件的研究進展[J]. 系統仿真學報, 1998, 10: 54-57.

[19]董敏, 馬昌飛, 楊利輝, 等. 基于AMESim與ADAMS聯合建模的軋機厚控系統仿真[J]. 機床與液壓, 2017, 45(5): 124-128.

[20]Zhang W, Yuan C, Liu G. Modeling simulation and experimental research of hydraulic pressure system in battery pole piece rolling mill[J]. China Mechanical Engineering, 2014, 25(24): 3267-3270.

[21]王文成, 張國山. 軸形工件感應加熱溫度控制系統開發[J]. 計算機與應用化學,  2012, 29(7): 773-777.

[22]王文成. 軋輥感應加熱系統開發[D]. 天津: 天津大學, 2012.

[23]李微. 邊部加熱器在熱軋寬帶鋼生產中的應用[J]. 冶金設備, 2009, (z1): 33-35.

[24]申世杰, 郭喜平. 熱應力對萬能軋機軋輥變形及其疲勞壽命影響分析研究[D]. 內蒙古: 內蒙古科技大學, 2011.

[25]胡仕成, 湛利華, 曹浚. 連續鑄軋輥套傳熱分析[J]. 中南工業大學學報(自然科學版), 2002, 33(3): 293-296.

[26]王永州, 畢鳳榮, 李晶華. 電池極片軋機軋輥有限元分析[D]. 天津: 天津大學, 2013.

[27]羅麗萍, 耿茂鵬, 徐萬里, 等. 基于ANSYS 軋輥表面電渣加熱渣池熱電場模擬研究[J]. 鑄造, 2009, 4: 353-356．

[28]肖述文. 熱輥壓機軋輥加熱過程熱應力分析和試驗研究[D]. 武漢: 武漢科技大學, 2015.

[29]陳國樑. 有限元法在熱軋輥三維溫度場分析中的應用[J]. 化纖與紡織技術, 2005, 5: 32-35.

[30] 陳國樑. ANSYS在熱軋輥研究中的應用[J]. 廣西紡織科技, 2005, 24(3): 43-46.

[31]李立新, 劉建洋, 袁翔等. CSP生產線高速鋼軋輥溫度場有限元分析[J]. 武漢科技大學(自然科學版). 2014(2): 121-124.

[32]蔡輝. 一種軋輥溫度控制系統:中國, CN202779192U[P]. 2013.03.13.

[33]Denis B. On thermal stress and fatigue life evaluation in work rolls of hot rolling mill[J]. Journal of Strain Analysis, 2012, 45(5): 297-312.

[34]Luks T, Horsky J, Nilsson A, et al. Contact stress distribution and roll surface temperatures in the roll gap analyzed with different sensors[J]. ROLLING, 2013, 7(10): 17-23.