基于COMSOL的城軌弓網電弧多物理場仿真研究

宋夢杰 于曉英 王澤 蘇洋 楊雋睿

DOI：10.16783/j.cnki.nwnuz.2024.01.008

收稿日期：20230330;修改稿收到日期：20231018

基金項目：甘肅省科技專員專項（23CXGA0023）；中國國家鐵路集團有限公司科技研究開發計劃重點課題（N2022X009）

作者簡介：宋夢杰（1996—），男，河南周口人，碩士研究生.主要研究方向為城市軌道交通弓網電弧溫度特性.

Email：smj02242023@163.com

*通訊聯系人，女，副教授，博士.主要研究方向為電氣化鐵路接觸網安全檢測監測.

Email：yxy1627@163.com

摘要：弓網電弧產生的高溫是燒蝕弓網系統材料的主要原因，研究弓網電弧發生機制，掌握弓網電弧溫度分布規律，有助于減少電弧對弓網系統材料的燒蝕損傷，增加弓網系統服役壽命.文中基于磁流體動力學理論，根據弓網電弧數學模型設置相關初始條件，采用有限元軟件COMSOL對弓網電弧數學模型進行求解，分析弓網電弧溫度分布及電弧兩極表面的溫度分布，并探究弓網電弧持續時間與弓網電弧溫度的關系.結果表明，弓網電弧中心區域溫度最高，接觸線表面溫度遠遠超過接觸線材料的熔點；弓網電弧持續時間與弓網電弧溫度呈正相關.

關鍵詞：城市軌道交通；弓網電弧；磁流體動力學；電弧溫度

中圖分類號：U 225.3??? 文獻標志碼：A??? 文章編號：1001-988Ⅹ（2024）01-0046-06

Research on multiphysical field simulation of pantograph arc

in urban rail transit based on COMSOL

SONG Meng-jie，YU Xiao-ying，WANG Ze，SU Yang，YANG Jun-rui

（School of Automation & Electrical Engineering，Lanzhou Jiaotong University，Lanzhou 730070，Gansu，China）

Abstract：The high temperature produced by the pantograph arc is the main reason for the ablation of the pantograph catenary system materials.Studying the mechanism of the pantograph arc and mastering the temperature distribution of the pantograph arc will help reduce the ablation of the arc on the pantograph catenary system materials and increase the service life of the pantograph catenary system.Based on the theory of magnetohydrodynamics，this paper establishes a mathematical model suitable for pantograph arc plasma by setting relevant initial conditions according to the mathematical model of pantograph arc and using the finite element software COMSOL to solve the mathematical model of pantograph arc，analyzes the temperature distribution of pantograph arc and the surface of the two poles of arc，and explores the relationship between the duration of pantograph arc and the temperature of pantograph arc.The results show that the temperature in the central area of the arc is the highest，the surface temperature of the contact wire far exceeds the melting point of the contact wire material，and the duration of the pantograph arc is positively correlated with the temperature of the pantograph arc.This work can provide a basis for reducing the number of arc occurrences in the pantograph and catenary，as well as the erosion effect of the arc on the pantograph and catenary materials.

Key words：urban rail transit;pantograph arc;magnetohydrodynamics;arc temperature

近年來，隨著我國交通運輸業的迅速發展，城市軌道交通因其具有載人流量大、速度快、安全穩定、綠色環保、方便快捷等優點在各中大城市得到了快速發展［1］.

弓網系統的良好接觸是保證城軌列車安全穩定運行的重要因素之一［2-4］.隨著科技的發展，城軌列車速度不斷提升，受電弓與接觸線之間振動加劇，導致弓網離線次數劇增，城軌弓網電弧頻發，嚴重影響城軌列車安全運行［5-6］.弓網電弧其實是氣體放電現象，在弓網離線瞬間，如果受電弓與接觸線之間的電壓大于弓網起弧電壓，弓網空氣間隙將會被擊穿，產生弓網電弧.發生弓網電弧時，會產生大量的熱量，反復作用于弓網系統材料表面，將會嚴重燒蝕弓網系統材料，影響弓網系統受流質量，縮短弓網系統服役壽命，可能導致接觸線斷裂，造成城軌列車運行事故［5-8］.因此有必要研究城軌弓網電弧溫度特性及電弧溫度分布情況，了解電弧產生機制，掌握弓網電弧溫度分布規律，可為延長弓網設備壽命、提高受流質量、減少燃弧次數等提供關鍵數據［9］.

弓網電弧現象是由熱場、流場、電磁場等多物理場相互耦合變化的過程，與開關電弧相比，由于弓網電弧維持著回路電流，更不容易熄滅［10-13］.對于弓網電弧，國內外學者開展了相關研究.胡怡等［14］通過光譜診斷法和紅外相機對降弓過程中弓網電弧溫度進行研究，得出在降弓過程中，接觸線表面溫度高于碳滑板表面溫度，碳滑板與接觸線溫度上升趨勢變緩.王君鵬［15］通過直流弓網電弧實驗，采集弓網電弧相關電氣數據，使用攝像機捕捉電弧的燃弧圖像，研究了弓網電弧相關特性，得出弓網系統在負向電源供電下容易產生弓網電弧，弓網之間電流越大燃弧就越劇烈.Wei等［16］提出弓網系統磨損率的計算模型，可近似預測弓網系統磨損情況，并通過實驗驗證該計算模型.預測結果與實際磨損相吻合，可為延長弓網系統服役壽命提供數據參考.Xu等［17］搭建了弓網電弧實驗平臺，研究了不同氣壓和氣流速度情況下電弧弧根的位置與時間，出現在低壓下電弧的起弧時間與弧根滯留時間明顯比常壓下的時間長，導致弓網系統燒蝕更加嚴重.

文中在前述研究弓網電弧的基礎上，借鑒開關電弧研究方法，基于磁流體動力學理論，結合麥克斯韋方程組與流體方程組，建立適用于城軌弓網電弧等離子體的數學模型.根據弓網電弧數學模型設置相關初始條件，并利用有限元軟件COMSOL對電弧數學模型進行仿真求解，分析仿真結果并探究弓網電弧持續時間與弓網電弧溫度的關系.

1? 城軌弓網電弧仿真建模

城軌弓網電弧內部物理化學過程十分復雜，為了簡化弓網電弧等離子體的計算過程，根據弓網電弧實際情況提出以下假設［18］：

① 假設弓網電弧穩定存在；

② 弓網電弧存在時，系統物性參數變化較小；

③ 弓網電弧等離子體磁導率為常數；

④ 弓網電弧存在時，電弧等離子體流場符合局部熱力學平衡要求.

1.1? 弓網電弧物理發生過程

弓網電弧發生過程是由熱場、流場和電磁場相互作用耦合形成的結果，相互作用過程如圖1所示.首先設置初始條件對模型初始化，通過設置初始壓力與溫度來確定弓網電弧等離子體的電導率，利用電磁方程得出電位分布和電場分布.

弓網系統中電流會產生焦耳熱，使得弓網電弧溫度快速升高，會導致弓網系統物性參數發生改變.同時弓網電弧受電磁力的影響，使得電弧等離子體流場不斷變化.弓網電弧整個系統會以熱傳導、熱對流和熱輻射的形式與外界流場不斷發生能量傳遞.同時列車高速行駛及弓網系統狀態的變化不僅會引起電弧等離子體流場的變化，也會改變弓網電弧系統物性參數，而物性參數的改變反過來也會影響流場的變化.整個物理過程不斷變化，相互影響，最終形成穩定的弓網電弧等離子體流場分布.

1.2? 弓網電弧數學模型

磁流體動力學理論是基于傳統流體力學與電磁學相結合發展而產生的理論，它是研究在電磁場中，帶電流體運動變化的一門學科［19］.磁流體傳熱理論包括流體動力學方程組和麥克斯韋電磁方程組.

1.2.1? 流體力學方程組

根據流體傳熱理論，流體動力學基本守恒方程包括質量守恒方程、動量守恒方程和能量守恒方程，具體形式如下.

ρt+（ρu）x+（ρv）y+（ρw）z=0;

（ρvi）t+div（ρvvi）=

∑2k=1xkηvixk+vkxi-

px+Svi;

（ρT）t+div（ρvT）=

divλCpgradT+ST.

（1）（2）（3）

其中，ρ為流體密度；u，v，w為各個方向的流體速度；vi為不同坐標上的速度分量；η為粘性系數；p為等離子流體壓強；T為電弧溫度；λ為等離子體熱導率；Cp為比熱容；Svi為流體動量守恒方程的源項；ST為流體能量守恒方程的源項.

（2）式中Svi表示動量守恒方程的源項為

Svi=j×（B0+Bi），（4）

其中，j為電流密度；B0為電弧自身電流所產生磁感應強度；Bi為接觸線電流所產生的磁感應強度.

（3）式中ST表示能量守恒方程的源項為

ST=ρt+1σj2+V-QR，（5）

其中，σ為電導率；V為電弧的粘性耗散項；QR為電弧輻射消耗的能量.

1.2.2? 麥克斯韋電磁方程

根據動量守恒方程中的源項方程可知，要求解源項就需要求出電流密度和磁感應強度，求解方程為

div（σgrad）=0，（6）

j=-σgrad.（7）

可通過（6）式和（7）式求解出電流密度j.

磁感應強度的求解包括兩部分，一部分是電弧自身電流產生的；另一部分是接觸線電流產生的.計算接觸線電流產生磁感應強度方法為，將導線等效為一根半無限長電流線，用畢奧-薩伐定律來求解，可通過（8）式得到.計算自身電流產生磁感應強度方法可以利用磁矢位方法進行求解，可通過（9）～（11）式得到，這里

Bi=14πrμieφ，（8）

其中，i為接觸線電流；μ為電弧等離子體磁導率；r為垂直距離.而

2×A=-μ0j，（9）

div（gradA）=-μ0j，（10）

B0=×A，（11）

其中，μ0為真空磁導率，A為磁矢位.

1.2.3? 輻射方程

熱輻射過程對電弧溫度的影響不可忽略，但是輻射散熱與吸收過程十分復雜，為了簡化計算，采用文獻［20］中介紹的簡化公式確定熱輻射方程，即

QR=4αk（T4-T40），（12）

其中，k為吸收系數，在1 atm標準大氣壓下k=13 m-1；α=5.67057 W·（m2K）-1為Stefan-Boltzmann常數.

1.3? 弓網幾何模型及初始條件

1.3.1? 弓網幾何模型

COMSOL有限元軟件計算二維模型和三維模型的方法一樣，將二維模對稱處理后為三維模型，因此，建立弓網系統幾何模型只需要二維模型即可.沿著接觸導線橫向方向建立幾何模型，弓網系統簡化幾何模型如圖2所示.

圖2中，接觸導線半徑為6.5 mm，為了簡化計算將受電弓滑板看作矩形，厚度為20 mm.接觸線材料為純銅導線，受電弓滑板材料為浸銅碳滑板，物性參數如表1所示.

1.3.2? 初始條件與邊界條件

弓網電弧數學模型的求解，需要對電弧仿真模型進行邊界條件與初始條件的設置.將接觸線設置為弓網電弧的陽極，受電弓滑板設置為弓網電弧陰極，陰極作為發射極發射電子，陽極作為接收極接收電子，所以零電勢作為陽極的邊界條件.整個環境初始溫度設置為300 K；系統壓強設置為1 atm；弓網系統電流設置為100 A；弓網間隙設置為4 mm；設置電弧持續時間為100 ms；對流方式選擇強制對流，設置傳熱系數為100 W·（m2·K）-1.

2? 城軌弓網電弧仿真結果分析

2.1? 弓網電弧溫度分布

發生弓網電弧時產生的高溫會對弓網系統表面材料造成燒蝕，嚴重影響機車受流.利用COMSOL仿真軟件，根據數學模型，建立弓網電弧幾何模型，設置邊界條件及弓網系統材料物性參數，加入方程中源項，接觸線為陽極，受電弓滑板為陰極，計算城軌弓網電弧溫度分布，如圖3所示.

由圖3可知，弓網電弧弧柱最高溫度可達到7 000 K以上，弓網電弧中心區域溫度最高，依次向四周逐漸遞減.電弧向兩極收縮，越靠近兩極收縮越明顯，兩極溫度分布情況也不相同，陽極附近溫度要高于陰極附近溫度.

利用仿真軟件定義一條二維截線為弧柱中心軸線，采集弧柱中心軸線上溫度數據并繪制溫度曲線圖，如圖4所示.

由圖4可知，電弧溫度分布并不對稱，電弧弧柱中心向陽極靠近時溫度下降速度比向陰極靠近下降速度快，在0～2 mm之間溫度變化最快.

城軌弓網電弧作用下，接觸線與受電弓滑板表面溫度分布如圖5所示.

由圖5可知，接觸線與受電弓滑板溫度分布情況都是表面溫度最高，溫度依次向內部減小.接觸線表面最高溫度可達7 730 K，接觸線材料一般為銅合金材料，銅的熔點為1 358 K，一般銅合金材料的熔點都比銅的熔點低.因此，在發生弓網電弧時，有可能會發生接觸線熔化現象.

根據傅里葉定律可知，熱流密度是與材料的導熱系數有關，接觸線使用純銅導線，導熱系數為400 W·（m·K-1），高于受電弓滑板的導熱系數，所以接觸線和滑板溫度分布相差較大.

2.2? 弓網電弧電流密度分布

由圖3可知，弓網電弧溫度從滑板到接觸線逐漸上升，這是因為發生弓網電弧現象時，滑板陰極電子向接觸線陽極移動，陽極表面電流密度升高，使得陽極處溫度比陰極處溫度高.具體弓網電弧電流密度分布如圖6所示.

從圖6可以看出，弓網電弧電流密度越靠近陽極區域，電流密度越高，電弧弧柱收縮越明顯，電流密度由中心區域逐漸向四周降低，在陽極最高可達2.69×107 A·m-2.

3? 弓網電弧持續時間對電弧溫度的影響

列車在運營時，電力機車輸送電力為DC1.5 kV，弓網系統將會通過大電流，容易產生電弧現象.弓網電弧產生的高溫會對弓網系統表面材料造成燒蝕，嚴重影響機車受流.研究弓網電弧持續時間與電弧溫度場的關系，掌握電弧溫度場變化規律，減少弓網電弧發生次數，提高機車受流質量.因此，文中使用仿真模型探究電弧溫度場隨時間的變化情況，設置弓網間隙為4 mm，電弧持續時間分別為60 ms，100 ms，140 ms和180 ms，如圖7所示.

從圖7可以看出，隨著電弧持續時間的增加，電弧覆蓋面積變大，電弧溫度呈上升趨勢.弓網電弧弧柱中心區域溫度最高，電弧弧柱中心區域溫度向四周逐漸降低.

為了更好地觀察弓網電弧持續時間與電弧溫度的關系，從弓網電弧溫度分布云圖中采集10次弓網電弧最高溫度隨時間的變化數據，如表2所示.弓網電弧最高溫度與電弧持續時間的關系曲線，如圖8所示.從圖8可以看出，弓網電弧最高溫度與弓網電弧持續時間有相關性，電弧持續時間越長，

電弧溫度就越高.因此，電弧持續時間與電弧溫度呈正相關，但并不是線性關系.

4? 結論

以磁流體動力學理論為基礎，建立與實際情況相符合的弓網電弧數學模型，使用有限元軟件COMSOL對數學模型進行求解，計算弓網電弧溫度場，分析仿真結果及探究弓網電弧溫度與持續時間的關系，得出結論如下：

1）弓網電弧弧柱中心區域溫度最高，依次向四周逐漸降低，電弧向兩極收縮，越靠近兩極收縮越明顯；由于兩極材料參數不同，溫度分布也不相同；接觸導線表面溫度遠超過接觸線材料的熔點，容易導致接觸線材料發生熔化現象，縮短弓網使用壽命.

2）弓網電弧電流密度越靠近陽極區域，電流密度越高，電弧弧柱收縮越明顯，電流密度由中心區域逐漸向四周降低；弓網電弧持續時間與弓網電弧溫度呈正相關，但不是線性關系.

文中研究結論有助減少弓網電弧發生次數，改進弓網系統材料，增強弓網系統材料抗燒蝕能力，對提高弓網系統服役壽命都有現實的工程意義.

參考文獻：

［1］? 于曉英.基于日盲區的城市軌道交通弓網電弧檢測系統的研制與應用［D］.蘭州：蘭州交通大學，2021.

［2］? 朱光亞，吳廣寧，韓偉峰，等.高速列車靜態升降弓時弓網電弧穩態特性仿真與分析［J］.鐵道學報，2016，38（2）：42.

［3］? 陳忠華，張忠傲，時光，等.弓網系統滑動電接觸區域溫度場仿真［J］.高壓電器，2018，54（3）：193.

［4］? JUNG K，SUH Y.Medium voltage power supply with enhanced ignition characteristics for plasma torches ［J］.Journal of Power Electronics，2011，11（4）：591.

［5］? 徐旻.弓網實驗平臺的研制及弓網電弧特性的研究［D］.北京：北京交通大學，2020.

［6］? 許潘.弓網電弧特性及其與弓網材料相互作用的建模仿真研究［D］.成都：西南交通大學，2019.

［7］? SIGNORINO D，GIORDANO D，MARISCOTTI A，et al.Dataset of measured and commented pantograph electric arcs in DC railways ［J］.Data in Brief，2020，31：105978.

［8］? CHO Y，LIM J，SEO H，et al.A series arc fault detection strategy for single-phase boost PFC rectifiers ［J］.Journal of Power Electronics，2015，15（6）：1664.

［9］? 徐建忠.弓網電弧的物理場仿真分析［D］.北京：北京交通大學，2012.

［10］? FAN F，WANK A，SEFERI Y，et al.Pantograph arc location estimation using resonant frequencies in DC railway power systems ［J］.in IEEE Transactions on Transportation Electrification，2021，7（4）：3083.

［11］? 郝靜，高國強，許潘，等.考慮電極熔化的弓網電弧與電極耦合模型［J］.高電壓技術，2018，44（5）：1668.

［12］? 于曉英，蘇宏升.基于PMT電壓一次積分值的城軌弓網電弧檢測系統［J］.鐵道學報，2019，41（9）：51.

［13］? 郝靜，高國強.受電弓降弓過程弓網電弧等離子體動態特性研究［J］.鐵道學報，2018，40（9）：65.

［14］? 胡怡，高國強，陳旭坤，等.基于光譜的降弓過程弓網電弧研究［J］.高電壓技術，2018，44（12）：3980.

［15］? 王君鵬.直流弓網電弧關鍵特性研究［D］.阜新：遼寧工程技術大學，2021.

［16］? WEI X，MENG H，HE J，et al.Wear analysis and prediction of rigid catenary contact wire and pantograph strip for railway system［J］.Wear，2020，442-443：203118.

［17］? XU Z，GAO G，WEI W，et al.Characteristics of pantograph-catenary arc under low air pressure and strong airflow［J］.High Volt，2022，7（2）：369.

［18］? 雷棟，張婷婷，段緒偉，等.列車運行速度對弓網電弧電氣特性的影響研究［J］.鐵道學報，2019，41（7）：50.

［19］? 李鑫，朱峰，邱日強，等.地鐵弓網電弧對機場儀表著陸系統影響研究［J］.鐵道學報，2018，40（5）：97.

［20］? KARETTA F，LINDMAYER M.Simulation of the gasdynamic and electromagnetic processes in low voltage switching arcs［J］.IEEE Transactions on Components Packaging and Manufacturing Technology，1998，21（1）：96.

（責任編輯? 孫對兄）