鋁電解槽鋼棒加高型陰極對鋁液中水平電流的優化

楊 帥，李 劼，徐宇杰，張紅亮，張翮輝，鄒 忠，賴延清

(中南大學 冶金科學與工程學院，長沙 410083)

鋁電解槽鋼棒加高型陰極對鋁液中水平電流的優化

楊 帥，李 劼，徐宇杰，張紅亮，張翮輝，鄒 忠，賴延清

(中南大學 冶金科學與工程學院，長沙 410083)

為減小傳統鋁電解槽內鋁液中存在的較大水平電流，提出一種陰極鋼棒加高型陰極。采用含電接觸的有限元模型計算鋼棒加高型陰極結構的電場，考察此種陰極對于鋁液中水平電流的優化作用。計算結果表明：選擇合適的加高位置及加高高度可使鋁液中X向水平電流密度的最大值和平均值有效減小，并且可以優化X向水平電流密度的分布，而對于減小Z向水平電流密度的最大值也有一定作用。以300 kA鋁電解槽的陰極為例，在距陰極炭塊邊緣800 mm處鋼棒上加高70 mm至90 mm時為最佳加高設計；同時鋼棒加高型陰極結構的電壓降小于傳統結構的，鋁液中的等勢線分布更平緩。

鋁電解；陰極壓降；加高鋼棒；水平電流

工業鋁電解槽內的電解質和鋁液兩種高溫熔體在電磁力、重力等作用下在槽膛內運動。熔體的運動促進了槽內的傳質與傳熱，對電解槽的穩定運行起著重要作用，然而熔體過大的水平流動及垂直波動都不利于電解槽的穩定運行，嚴重時甚至引起滾鋁等極端現象[1?3]。熔體的運動對電解的經濟技術指標如電耗、效率等也有顯著影響，因此，追求槽內熔體合理的流動及盡可能小的波動，一直是鋁電解槽研究與設計的核心之一。

電磁力是引起槽內熔體運動的主要驅動力，由于熔體流動及波動產生的鋁液?電解質界面變形會引起額外的擾動電磁力，擾動電磁力進一步改變熔體的流動及鋁液?電解質界面的波動，在此過程可能產生波動的不穩定分量，若這些不穩定分量得不到有效的抑制就會引起槽內磁流體的不穩定現象[4]。LINDSAY和DAVIDSON[5]以及BOJAREVICS和ROMERIO[6]研究指出，由鋁液中水平電流與垂直磁場作用產生的擾動電磁力是唯一有重要影響的擾動電磁力組分，因此，從穩定磁流體波動的角度考慮，優化鋁液中的垂直磁場以及水平電流都能直接改善槽內磁流體的穩定性。

為改善鋁電解槽中的磁流體穩定性，國外與國內的學者都有許多研究，這部分的工作主要集中在優化電解槽的母線結構及進電制度以及提出更有利于鋁液中減小垂直磁場的新型陰極結構。GUSEV等[7]以及ZIEGLER和RUAN[8]運用數值仿真的方法，設計出用于電解系列最末端槽的磁場補償母線結構，并將其結果運用于工業生產中。周萍等[9]、李茂和周孑民[10]以及LI等[11]分別運用數值仿真的方法研究不同的母線結構及進電制度對鋁液中垂直磁場的影響及對熔體流動的影響，優化了特定槽型的母線配置及進電制度。劉業翔等[12]設計的底部出電鋁電解槽，不僅具有較為理想的垂直磁場分布，而且鋁液中的水平電流也相對較小，為實現鋁電解槽的大幅節能創造了條件。與以上研究相比，對于減小鋁液中水平電流的研究相對要少，而MOREAU和ZIEGLER[13]在很早就指出了減小鋁液中水平電流對于增強磁流體穩定性的意義。QI等[14]計算了陰極結構中的電場，研究鋁液高度及伸腿長度對鋁液中的水平電流密度的影響，指出鋁液中水平電流密度的分布狀況和鋁液高度及伸腿對鋁液中水平電流密度的影響規律。張艷芳等[15]對陰極炭塊類型和鋼棒類型以及鋼棒高度及寬度變化對于鋁液中水平電流的影響進行了分析，指出鋼棒的電導率越高或增加鋼棒的寬度都有助于減小水平電流的分量，并使水平電流密度的分布更均勻。但由于上述針對鋁液中水平電流的兩項研究在陰極結構上都保持傳統形狀，而鋁液高度及伸腿長度受工藝條件制約并不能隨意改變，因此對于優化鋁液中的水平電流作用有限。

本文作者提出一種陰極鋼棒加高型的新型陰極結構，并且運用此前開發的采用電接觸模型的有限元分析方法[16]計算整個陰極結構的電場，探索此種陰極結構對于減小鋁液中水平電流的優化作用，初步獲得了合理的結構參數，結果可為改善鋁電解槽中磁流體穩定性提供一條思路。

1 鋼棒加高型陰極結構簡介

鋁電解槽中，電流從陽極導桿流入槽內，經過陽極炭塊、電解質、鋁液及陰極炭塊與鋼棒后從槽側部流出，由于電流的流出方向與流入方向成90°角，并且鋁液的電阻率遠遠小于電解質以及陰極炭塊的，因此電流的轉向大部分發生在鋁液層中，形成了鋁液中的水平電流，水平電流與垂直磁場共同作用引起電解質?鋁液界面的波動，從而影響到槽內磁流體的穩定性。傳統的鋁電解槽剖面結構及電流流向示意圖如圖1所示，半陰極結構剖面圖如圖2所示。

圖1 預焙陽極鋁電解槽結構示意圖Fig. 1 Structure diagram of aluminum reduction cell: 1—Aluminum rod; 2—Steel draw; 3—Anode carbon block; 4—Side ledges; 5—Electrolyte; 6—Aluminum melt; 7—Cathode carbon block; 8—Collector bar; 9—Side carbon block; 10—Surrounding paste draw; 11—Refractory brick; 12—Insulating brick; 13—Steel shell

本文作者所提出的鋼棒加高型陰極結構的半陰極剖面如圖3所示，其實質是在傳統陰極結構的基礎上對嵌入陰極炭塊內的鋼棒上對稱的某兩點從兩端進行線性加高(圖中所示其中的一點)，而伸出陰極炭塊部分的鋼棒高度保持不變。由于鋼棒的導電性優于陰極炭塊的，因此，加高型鋼棒對于陰極部分的電場分布會有一定影響，從而改變鋁液中的電流分布，為優化鋁液中水平電流提供了可能性。在此以〈a,b〉表示不同的加高方案，其中a代表加高點距離陰極炭塊出電端邊緣的水平距離，b表示加高點與初始結構的高度差，單位均為mm，其具體所表示的位置標于圖3中。

圖3 鋼棒加高型陰極結構示意圖Fig.3 Schematic diagram of new cathode structure with heightened collector bar

2 理論基礎及模型

2.1 控制方程

鋁電解槽中導電部分的導電微分方程可表示為拉普拉斯方程，即

式中：V為標量電位；I為電流；R為電阻；σ為電導率。

2.2 接觸壓降

陰極鋼棒與鋼棒糊間的接觸壓降可由下式計算：

式中：cφ為接觸電壓；I為電流；Rc為接觸電阻率；S為接觸面積。

2.3 模型及邊界條件

以某300 kA預焙鋁電解槽陰極作為研究對象，其初始結構及工藝參數列于表1。此外，加高型鋼棒的模型圖如圖4所示。

考慮到陰極結構的對稱性，截取其1/2為計算模型，同時為建模的便捷及計算的效率，不考慮電解槽中不同陰極上通過的電流在大小上的差異，即假定每塊陰極上通過電流為總電流在所有陰極上的平均值；不考慮各材料電阻率隨溫度的變化；不考慮鋁液界面的變形。

表1 陰極初始結構及工藝參數Table 1 Initial structure parameters of cathode

圖4 加高型鋼棒模型示意圖Fig. 4 Schematic diagram of heightened collector bars model

模型在ANSYS中生成并進行網格劃分及計算，具體模型及計算流程見文獻[16]。半陰極有限元模型如圖5所示，X向表示指向電解槽水平短軸方向，Y向指向電解槽高度方向，Z向指向電解槽水平長軸方向。

2.4 材料屬性

計算中所需的材料屬性見表2。此外，鋼棒糊與陰極鋼棒間的接觸電阻率根據文獻[16]取1×10?5?·m。

表2 實驗材料的電阻率Table 2 Electrical resistivity of experimental materials

圖5 半陰極有限元模型Fig. 5 FEM model of half cathode

3 結果與分析

3.1 鋼棒加高位置對鋁液中水平電流的影響

分別計算初始結構、鋼棒加高方案分別為〈300,50〉、〈500,50〉、〈700,50〉、〈800,50〉、〈1 000,50〉的陰極結構的電場，考察鋼棒加高位置對鋁液中水平電流的影響。各方案鋁液中X向水平電流密度分布分別如圖6所示，Z向水平電流密度分布如圖7所示，電流密度最大值與平均值列于表3。

由表3可知，隨著鋼棒加高位置向陰極內部延伸，鋁液中X向水平電流密度的最大值先減小，繼而在距炭塊邊緣800 mm左右的某一位置開始增大，但增大的趨勢比減小的趨勢小，而X向水平電流密度平均值則隨著加高位置向陰極內部延伸而不斷減小，減小的絕對量不太大；陰極鋼棒的加高位置對Z向水平電流密度最大值的影響與對X向水平電流密度最大值的影響有相似的規律，但其影響的絕對量要小得多；此外，鋼棒加高會增大鋁液中的Z向水平電流密度平均值，但這種影響非常小，并且大小基本不隨加高位置變化而變化，即加高位置的選擇對其幾乎沒有影響。

圖6 鋼棒不同加高位置鋁液中X向水平電流密度分布Fig. 6 Horizontal current density alongXaxis in metal pad with various heightened positions of collector bars (A/m2): (a) Initial structure; (b) 〈300,50〉; (c) 〈500,50〉; (d) 〈700,50〉; (e) 〈800,50〉; (f) 〈1 000,50〉

圖7 鋼棒不同加高位置鋁液中Z向水平電流密度分布Fig. 7 Horizontal current density alongZaxis in metal pad with various heightened positions of collector bars (A/m2): (a) Initial structure; (b) 〈300,50〉; (c) 〈500,50〉; (d) 〈700,50〉; (e) 〈800,50〉; (f) 〈1 000,50〉

表3 鋼棒不同加高位置鋁液中水平電流密度Table 3 Horizontal current density in metal pad with various heightened positions of collector bars

由圖6中X向水平電流密度分布也可以看出，不同加高位置的鋁液中X向水平電流密度的分布形態基本一致，電流密度較大的區域分布于偏靠近陰極邊緣的位置，鋼棒加高位置為距炭塊邊緣為800 mm處時X向水平電流密度分布最為均勻，電流密度的變化梯度較小，并且最大值達到最小。圖7所示的不同加高位置方案的鋁液中Z向水平電流密度分布形態基本相同，電流密度較大的區域分布于靠近炭塊邊緣兩端的部位，所占體積較小并且向鋁液內部逐漸減小，僅在電流密度最大值上有很小的差別。

綜上所述，鋼棒加高對于減小鋁液中的水平電流以及優化其分布有積極效果，加高位置的選擇主要影響鋁液中的X向水平電流的分布及大小，加高位置在距炭塊邊緣800 mm處時電流密度最大值達到最小，平均值也有所減小，鋁液中X向水平電流密度的分布更加均勻，因而距炭塊邊緣800 mm處是較為合適的加高位置。

3.2 鋼棒加高高度對鋁液中水平電流的影響

分別計算鋼棒加高方案分別為〈800,30〉、〈800,50〉、〈800,70〉、〈800,90〉、〈800,110〉的陰極結構的電場，考察鋼棒加高高度對于鋁液中水平電流的影響。各方案鋁液中X向水平電流密度分布分別如圖8所示，Z向水平電流密度分布如圖9所示，電流密度最大值與平均值列于表4。

由表4可知，鋼棒加高高度對鋁液中X向水平電流密度影響較大，隨著鋼棒加高高度的增加，鋁液中X向水平電流密度的最大值先減小再增大，在加高高度為70 mm左右達到最小值，而X向水平電流密度的平均值則隨著鋼棒加高高度的增加而不斷減小；陰極鋼棒的加高高度對鋁液中的Z向水平電流密度也有一定影響，Z向水平電流密度的最大值隨著加高高度的增加而不斷減小，但平均值隨加高高度的增加而不斷增加。

由圖8可以看出，隨著鋼棒加高高度的增加，鋁液中X向水平電流密度較大的區域向槽中心移動，并且趨向于較為平均的分布，但當加高高度超過90 mm附近的某一值時，開始再次趨向于集中分布。由圖9可知，鋁液中Z向水平電流密度分布形態基本不受加高高度變化的影響，僅在最大值上有所差別，并且這種差別在數值上相對于X向水平電流密度的變化要小得多。

綜上所述，鋼棒的加高高度對鋁液中X向水平電流密度的分布及大小影響較大，同時對Z向水平電流的也有一定影響，雖然更高的鋼棒加高高度能更多的減小鋁液中X向水平電流密度的平均值，但從避免Z向水平電流密度平均值升高及較好的X向水平電流密度分布形態考慮，加高高度為70 mm到90 mm是較為合理的加高高度。在加高方案為〈800,70〉時，X向水平電流密度最大值相對于初始結構下降了28.36%，平均值下降了11.11%。

圖8 鋼棒不同加高位置鋁液中X向水平電流密度分布Fig. 8 Horizontal current density alongXaxis in metal pad with various heightened positions of collector bars (A/m2): (a) 〈800,30〉; (b) 〈800,50〉; (c) 〈800,70〉; (d) 〈800,90〉; (e) 〈800,110〉

圖9 鋼棒不同加高位置鋁液中Z向水平電流密度分布Fig. 9 Horizontal current density alongZaxis in metal pad with various heightened positions of collector bars(A/m2): (a) 〈800,30〉; (b) 〈800,50〉; (c) 〈800,70〉; (d) 〈800,90〉; (e) 〈800,110〉

表4 鋼棒不同加高位置鋁液中水平電流密度Table 4 Horizontal current density in metal pad with various heightened positions of collector bars

3.3 鋼棒加高對陰極壓降的影響

分別計算初始結構、鋼棒加高方案分別為〈300,50〉、〈500,50〉、700,50〉、〈800,50〉、〈1 000,50〉以及〈800,30〉、〈800,50〉、〈800,70〉、〈800,90〉、〈800,110〉的陰極結構壓降，考察不同鋼棒加高方案對于陰極結構電壓降的影響。不同方案的陰極壓降列于表5中，其中初始結構及鋼棒加高方案為〈800,70〉時鋁液中X向截面等電勢線分布如圖10所示。

表5 不同陰極結構的電壓降Table 5 Voltage drop with different cathode structures

圖10 鋁液中X向截面等電勢線分布Fig. 10 Isopotential curves in section alongXaxis of metal pad: (a)Initial structure; (b) 〉〈800,90

由表5可知，由于鋼棒加高使得陰極結構中電阻率較高的陰極炭塊部分被鋼棒取代，所以鋼棒加高型陰極結構的電壓降都小于初始結構的壓降，且加高位置越靠近炭塊邊緣或加高高度越高，則陰極結構的壓降越小；在優化水平電流最優的加高方案〈800,70〉時陰極結構的壓降為283 mV，相對于初始結構降低了19 mV。

此外，由圖10可以看出，鋼棒加高型陰極結構的鋁液中X向截面上等電勢線分布形態與初始結構相似，但分布更加平緩，更趨向于在鋁液中沿垂直的方向分布，而這是對于減小鋁液中水平電流最根本的原因。

4 結論

1) 鋼棒加高對于減小鋁液中的水平電流以及優化其分布有積極效果。隨著鋼棒加高位置向陰極內部延伸，鋁液中X向水平電流密度的最大值先減小，加高位置距炭塊邊緣800 mm處時電流密度最大值達到最小，平均值也有所減小，鋁液中X向水平電流密度的分布最均勻。

2) 隨著鋼棒加高高度的增加，鋁液中X向水平電流密度的最大值先減小再增大，在加高高度為70 mm左右達到最小值，平均值則不斷減小；Z向水平電流密度的最大值隨加高高度的增加而不斷減小，但平均值不斷增加，最優加高高度為70 mm至90 mm。

3) 在最優鋼棒加高方案〈800,70〉時，X向水平電流密度最大值相對于初始結構下降了28.36%，平均值下降了11.11%，并且陰極結構的電壓降相對于初始結構降低了19 mV，鋁液中的等勢線分布更加平緩。

REFERENCES

[1]LI Xiang-peng, LI Jie, LAI Yan-qing, ZHAO Heng-qin, LIU Ye-xiang. Mathematical simulation of gas induced bath flow in drained aluminum reduction cell[J]. Transactions of Nonferrous Metals Society of China, 2004, 14(6): 1221?1226.

[2]劉業翔, 李 劼. 現代鋁電解[M]. 北京: 冶金工業出版社, 2008: 222. LIU Ye-xiang, LI Jie. Modern aluminum electrolysis[M]. Beijing: Metallurgical Industry Press, 2008: 222.

[3]李 劼, 徐宇杰, 劉 偉, 賴延清. 基于波動模態耦合的鋁電解槽磁流體穩定性傅立葉級數法分析[J]. 計算力學學報, 2010, 27(2): 213?217. LI Jie, XU Yu-jie, LIU Wei, LAI Yan-qing. Aspect ratio research based on magneto-hydrodynamic stability in aluminum reduction cells[J]. Chinese Journal of Computational Mechanics, 2010, 27(2): 213?217.

[4]徐宇杰, 李 劼, 張紅亮, 賴延清. 基于非線性淺水模型的鋁電解磁流體動力學計算[J]. 中國有色金屬學報, 2011, 21(1): 191?197. XU Yu-jie, LI Jie, ZHANG Hong-liang, LAI Yan-qing. MHD calculation for aluminium electrolysis based on nonlinear shallow water model[J]. The Chinese Journal of Nonferrous Metals, 2011, 21(1): 191?197.

[5]LINDSAY R I, DAVIDSON P A. Stability of interfacial waves in aluminum reduction cell[J]. J Fluid Mech, 1998, 362: 327?331.

[6]BOJAREVICS V, ROMERIO M V. Long waves instability of liquid metal-electrolyte interface in aluminum electrolysis cells: A generalization of Sele's criteria[J]. European Journal of Mechanics B, 1994, 13: 33?56.

[7]GUSEV A, KRYLOV L, PLATONOV V, VABISHCHEVICH P. Compensation of the magnetic field of the extreme reduction cells in potlines[C]//KVANDE H. Light Metals 2005. SanFrancisco, CA: TMS, 2005: 475?48.

[8]ZIEGLER D P, RUAN Y M. Busbar arrangement optimization for end cells[C]//BEARNE G. Light Metals 2009. Warrendale, PA: TMS, 2009: 535?538.

[9]周 萍, 周乃君, 梅 熾, 姜昌偉, 蔡祺風. 鋁電解槽內鋁液電磁攪拌流動的數值模擬[J]. 過程工程學報, 2003, 3(4): 295?301. ZHOU Ping, ZHOU Nai-jun, MEI Chi, JIANG Chang-wei, CAI Qi-feng. Numerical simulation of electromagnetic stirring flow of metal pad in aluminum reduction cells[J]. The Chinese Journal of Process Engineering, 2003, 3(4): 295?301.

[10]李 茂, 周孑民. 大型鋁電解槽母線配置的數值仿真與優化[J]. 中南大學學報: 自然科學版, 2009, 40(3): 562?567. LI Mao, ZHOU Jie-min. Numerical simulation and optimization of busbar configuration in large aluminum electrolysis cell[J]. Journal of Central South University: Science and Technology, 2009, 40(3): 562?567.

[11]LI Jie, LIU Wei, LAI Yan-qing, Li Qing-yu, LIU Ye-xiang. Coupled simulation of 3D electro-magneto-flow field in Hall-Heroult cells using finite element method[J]. Acta Metallurgica Sinica: English Letters, 2006, 19(2): 105?116.

[12]劉業翔, 梁學民, 李 劼, 張紅亮, 徐宇杰, 丁鳳其, 鄒 忠.底部出電型鋁電解槽母線結構與電磁流場仿真優化[J]. 中國有色金屬學報, 2011, 21(7): 1688?1695. LIU Ye-xiang, LIANG Xue-min, LI Jie, ZHANG Hong-liang, XU Yu-jie, DING Feng-qi, ZOU Zhong. Simulation and optimization of bus structure and electro-magneto-flow field of aluminum reduction cells with vertical bottom bars[J]. The Chinese Journal of Nonferrous Metals, 2011, 21(7): 1688?1695.

[13]MOREAU R J, ZIEGLER D. Stability of aluminium cells-a new approach[C]//MILLER R E. Light Metals 1986. New Orleans, LA: TMS, 1986: 359?364.

[14]QI Xi-quan, Lü Ding-xiong, MA Shao-xian, WANG De-quan, MAO Ji-hong. Study of current distribution in the metal pad of aluminum reduction cells[C]//BEARNE G. Light Metals 2009. Warrendale, PA: TMS, 2009: 575?580.

[15]張艷芳, 趙慶云, 劉吉波, 柴登鵬, 唐新平, 周益文, 王躍勇.鋁電解槽陰極對水平電流影響的模擬與分析[C]//2007中國國際鋁冶金技術論壇論文集. 北京: 冶金工業出版社, 2007: 540?543. ZHANG Yan-fang, ZHAO Qing-yun, LIU Ji-bo, CHAI Deng-peng, TANG Xin-ping, ZHOU Yi-wen, WANG Yue-yong. The simulation and analysis of the affection of the cathode to the horizontal current in the aluminum reduction cell[C]// Proceedings of 2007 China International Conference on Aluminum Metallurgy. Beijing: Metallurgical Industry Press, 2007: 540?543.

[16]LI Jie, LIU Wei, LAI Yan-qing, WANG Zhi-gang, LIU Ye-xiang. Analysis of cathode voltage drop in aluminum electrolysis cells with an electric contact model[C]//SORLIE M. Light Metals 2007. Orlando, FL: TMS, 2007: 465?470.

(編輯 李艷紅)

Optimization of horizontal current in metal pad by using cathode with heightened collector bars in aluminum reduction cell

YANG Shuai, LI Jie, XU Yu-jie, ZHANG Hong-liang, ZHANG He-hui, ZOU Zhong, LAI Yan-qing
(School of Metallurgical Science and Engineering, Central South University, Changsha 410083, China)

Based on the target of reducing the large horizontal current in metal pad of aluminum reduction cells, a kind of cathode with heightened collector bars was improved. The electric field of the cathode was calculated by using finite element method with an electric contact model, and the optimization effect to horizontal current in metal pad was examined. The results show that an appropriate design scheme of heightened collector bars can effectively reduce the maximum and average values of horizontal current density alongXaxis in the metal pad, and optimize the current density. The maximum value of horizontal current density alongZaxis in metal pad can be reduced at the same time. For the case of cathodes in 300 kA aluminum reduction cells, the best design of heightened position is at a distance of about 800 mm from the edge of cathode carbon block, and the beneficial heightened height ranges from 70 to 90 mm, the cathode with heightened collector bars has a lower voltage drop than the traditional one, and the iso-potential curves in the metal pad shows slightly.

aluminum electrolysis; cathode voltage drop; heightened collector bar; horizontal current

TF821

A

國家自然科學基金青年基金資助項目(51104187)；教育部高等學校博士學科點專項科研基金資助項目(20100162120008)

2011-09-19；

2012-03-20

張紅亮，講師，博士；電話：0731-88830474；E-mail: net_hotang@163.com

1004-0609(2012)10-2951-09