楊 康，黃焯麒，于 龍，陳唐龍

（西南交通大學(xué) 電氣工程學(xué)院，四川 成都610031）

0 引 言

隨著高速鐵路的迅速發(fā)展，電氣化列車的運行速度越來越快，高速氣流產(chǎn)生的空氣動力對列車運行的影響越來越大。目前研究主要集中在列車空氣動力特性及其與列車和周圍環(huán)境相互影響的方面［1］，高速列車受電弓空氣動力學(xué)效應(yīng)對弓網(wǎng)動力學(xué)性能的影響研究并不多。所謂弓網(wǎng)受流，是指電氣化列車在運行過程中，通過安裝在車頂?shù)氖茈姽c道路上空的接觸網(wǎng)懸掛系統(tǒng)之間的滑動接觸，獲取電能來驅(qū)使列車運動的過程。受電弓良好的空氣動力性能是保證弓網(wǎng)受流質(zhì)量的重要因素。當(dāng)前對弓網(wǎng)動力學(xué)的研究很多，但是大多沒有考慮空氣動力的影響［2－4］，這在列車低速運行時是可以的，但在高速情況下受電弓的空氣動力特性，尤其是受電弓受到的空氣抬升力對弓網(wǎng)關(guān)系的影響至關(guān)重要。

對受電弓空氣動力性能的研究主要通過2種方法：風(fēng)洞 實 驗［5－7］和 數(shù) 值 計 算［8，9］。風(fēng) 洞 實 驗 所 需 的 成 本 高，數(shù) 值計算方法沒有考慮受電弓的真實結(jié)構(gòu)，很難得到準(zhǔn)確的結(jié)果。隨著計算機(jī)技術(shù)的發(fā)展，通過計算機(jī)仿真的方法研究受電弓的空氣動力特性成為一種快速有效的方法。本文采用有限體積法結(jié)合N－S方程和連續(xù)方程，建立了受電弓的實體模型，分析了受電弓在不同車速下受到的空氣動力特性；建立了弓網(wǎng)動態(tài)耦合的有限元模型，分析了受電弓空氣抬升力對弓網(wǎng)動態(tài)特性的影響。

1 模型的建立

1.1 受電弓空氣動力學(xué)模型

高速列車在以空氣為介質(zhì)的空間中運行，研究對象可歸納為研究空氣流體以及在空氣流體中的剛體間相互作用力的問題。在空氣動力學(xué)研究中，忽略流體系的份子結(jié)構(gòu)，將流體看作是連續(xù)介質(zhì)，其中不存在真空、分子間無間隙和分子運動。由于流體的連續(xù)性，將流體的運動學(xué)得各個基本量 （速度、壓力等）看作是某些物理量 （時間、空間坐標(biāo)）的連續(xù)函數(shù)，用連續(xù)函數(shù)的數(shù)學(xué)模型求解空氣動力學(xué)問題。

建立高速受電弓的簡化模型，如圖1 所示。受電弓運行的計算區(qū)域的長為22m，寬為12m，高為5m。由于受電弓結(jié)構(gòu)復(fù)雜，計算區(qū)域用非結(jié)構(gòu)網(wǎng)格進(jìn)行離散，受電弓表面為三角形網(wǎng)格，空氣區(qū)域為四面體網(wǎng)格?？紤]到距離受電弓越遠(yuǎn)的區(qū)域?qū)τ嬎憬Y(jié)果的影響越小，采用漸變方式進(jìn)行網(wǎng)格劃分。受電弓表面的網(wǎng)格劃分比較密集，網(wǎng)格大小為10 mm；離受電弓越遠(yuǎn)的區(qū)域，網(wǎng)格越稀疏，流場最外層網(wǎng)格大小為500～600 mm，如圖2所示。

圖1 受電弓簡化仿真模型

圖2 受電弓計算區(qū)域網(wǎng)格劃分

不考慮環(huán)境風(fēng)對列車的影響，在列車運行速在160 km／h到350km／h時，對應(yīng)的馬赫數(shù)小于0.3，空氣的壓縮性可以不考慮。受電弓周圍空氣場一般可以看做為定常、等溫、不可壓縮的三維流場，通常具有較高的雷諾數(shù)，應(yīng)按湍流處理［1］。采用標(biāo)準(zhǔn)k－ε 模型來模擬湍流模型，其方程為［10］

其中

式中：ρ——空氣密度；k——湍流動能；ε——湍流耗散率；xi或xj——坐標(biāo)的3個分量，分別代表x，y，z 這3個方向坐標(biāo)；ui或uj——列車周圍流場速度，分別代表ux，uy，uz這3個坐標(biāo)方向的速度分量；μ——流體的黏性系數(shù)；C1，C2——經(jīng)驗常數(shù)；σk，σε——分別為湍流動能和湍流動耗散率對應(yīng)的普朗特數(shù)。

根據(jù)氣體黏性應(yīng)力與氣體運動速度之間的關(guān)系，運用牛頓第二定律，忽略空氣的質(zhì)量力，得到直角坐標(biāo)系下不可壓縮黏性流體的3 個運動方程，即不可壓縮流的N－S方程

受電弓周圍流體運動遵循質(zhì)量守恒定律，可得到連續(xù)方程

1.2 受電弓—接觸網(wǎng)系統(tǒng)模型

接觸網(wǎng)采用簡單鏈形接觸網(wǎng)，由承力索、接觸線、吊弦、線夾以及定位器等組成，是一個復(fù)雜的機(jī)械系統(tǒng)。受電弓由上框架、下框架、弓頭等結(jié)構(gòu)組成。這樣的系統(tǒng)很難做出其真實的模型，通常對其模型進(jìn)行簡化。承力索和接觸線的有限元模型采用歐拉－伯努利梁，吊弦采用彈簧單元。受電弓采用三元質(zhì)量塊單元。弓網(wǎng)耦合系統(tǒng)的模型［11］如圖3所示。

圖3 弓網(wǎng)仿真模型

圖3中，Tc是承力索的張力；Tj是接觸線的張力；m1、m2、m3分別是受電弓弓頭、上框架和下框架的等效質(zhì)量；k1、k2、k3分別是弓頭與上框架、上框架與下框架、下框架與機(jī)車頂部之間的等效剛度；c1、c2、c3分別是弓頭與上框架、上框架與下框架、下框架與機(jī)車頂部之間的等效阻尼；y1、y2、y3分別是弓頭、上框架和下框架的位移；F0是受電弓的靜態(tài)抬升力。

弓網(wǎng)耦合系統(tǒng)的運動方程可簡寫為矩陣形式［2，12］

式中：［M］——弓網(wǎng)系統(tǒng)的質(zhì)量矩陣；［C］——弓網(wǎng)系統(tǒng)的阻尼矩陣；［K］——弓網(wǎng)系統(tǒng)的剛度矩陣；［F］——外載荷矩陣；y——位移矩陣。受電弓和接觸網(wǎng)的仿真參數(shù)見表1和表2。

表1 受電弓模型參數(shù)

表2 接觸網(wǎng)參數(shù)

2 仿真結(jié)果與分析

2.1 受電弓空氣動力特性仿真

將受電弓計算區(qū)域的頂面和側(cè)面設(shè)置為滑移壁面；受電弓表面和地面設(shè)置為無滑移壁面；出口截面設(shè)置為壓力出口邊界；入口截面設(shè)置為速度入口邊界，速度值為列車的運行速度。

利用FLUENT 流場計算軟件，采用標(biāo)準(zhǔn)k－ε模型來模擬受電弓周圍的空氣流場，流場中的控制方程可用公式（1）～式 （5）來描述，仿真迭代800次，計算得到受電弓開口運行時各部件的氣動阻力和氣動升力。圖4顯示的是車速350km／h情況下受電弓所受的阻力和升力在迭代運行過程中的變化情況。從圖中可以看出，迭代約400 次后，受電弓受到的空氣阻力和升力趨于穩(wěn)定。整弓的阻力穩(wěn)定在902.2N 左右，升力穩(wěn)定在49.3N 左右。受電弓的阻力與列車的運行方向相反，阻礙受電弓的前行，是受電弓產(chǎn)生氣動噪聲的主要來源；受電弓的升力方向垂直于路面，使受電弓框架和弓頭產(chǎn)生垂向的振動，影響受電弓－接觸網(wǎng)之間接觸受流的質(zhì)量。

圖4 受電弓空氣特性仿真

不同速度情況下，受電弓開口運行各部件所受到的空氣動力值如表3所示。由表3可知，受電弓弓頭由于結(jié)構(gòu)復(fù)雜，所受的阻力最大；上下臂桿的空氣升力大小相近，方向相反，基本能夠相互抵消，受電弓受到的空氣升力主要來自于弓頭。由此可知，受電弓的弓頭是空氣動力的主要受力部件。

表3 開口運行時受電弓各部件的氣動力／N

受電弓整弓的氣動阻力和氣動升力隨著速度的增加而增大，與速度的平方成正比，如圖5所示。受電弓的靜態(tài)抬升力為70N，當(dāng)運行速度為160km／h時，整弓的升力約為11.5 N，是靜態(tài)抬升力的16.4%；而速度達(dá)到350 km／h時，整弓的升力達(dá)到49.4 N，是靜態(tài)抬升力的70.7%。因此，在高速條件下的弓網(wǎng)受流，必須考慮受電弓空氣升力的影響。

2.2 弓網(wǎng)動態(tài)仿真

圖5 受電弓空氣動力與速度的關(guān)系

利用ANSYS軟件建立弓網(wǎng)耦合仿真模型，如圖6 所示。計算弓網(wǎng)接觸壓力時，只有垂向的作用力才能對其產(chǎn)生影響，因此只需考慮受電弓空氣抬升力的作用。分別將受電弓各部件的氣動升力加載于受電弓三元質(zhì)量塊的相應(yīng)位置，通過弓網(wǎng)耦合動力學(xué)仿真計算得到車速為350km／h時的弓網(wǎng)接觸壓力。圖7 給出了考慮 （虛線）和未考慮（實線）受電弓氣動升力情況下接觸壓力的比較。由圖可知，考慮受電弓氣動升力的影響后，接觸力明顯增大，且接觸力的波動性也增大。

圖6 弓網(wǎng)仿真模型的建立

表4為考慮和未考慮受電弓氣動升力作用時的接觸壓力數(shù)據(jù)比較。由表可知，考慮受電弓氣動升力后，接觸力的平均值增加了19.5N，最大值增加了42.2N，最小值和標(biāo)準(zhǔn)差均有增加?？梢娛茈姽臍鈩由佑|壓力的影響在高速時十分的明顯。

圖7 空氣升力對弓網(wǎng)接觸壓力的影響

表4 接觸壓力數(shù)據(jù)比較

3 結(jié)束語

本文利用流體分析軟件仿真計算了受電弓在高速運行時的空氣動力特性，并在考慮受電弓空氣動力的情況下，對弓網(wǎng)動態(tài)接觸壓力進(jìn)行了分析，主要得到以下兩點結(jié)論：

（1）仿真計算了受電弓運行時，周圍空氣場對受電弓的空氣動力特性。結(jié)果表明，受電弓所受到的空氣阻力和升力隨速度的增加而增大，與速度的平方成正比，高速情況下的弓網(wǎng)受流，必須考慮氣動升力的影響。其中弓頭是主要的受力部件。

（2）分析了高速情況下受電弓空氣升力對弓網(wǎng)接觸壓力的影響。結(jié)果表明空氣升力增大了弓網(wǎng)接觸壓力，且接觸力的波動性也增大，使得弓網(wǎng)受流質(zhì)量變差。

［1］TIAN Hongqi.Train aerodynamic［M］.Beijing：China Railway Press，2007 （in Chinese）. ［田紅旗.列車空氣動力學(xué)［M］.北京：中國鐵道出版社，2007.］

［2］ZHOU Ning，ZHANG Weihua.Dynamical performance simulation of the pantograph－catenary coupled system based on direct intergration method ［J］.China Railway Science，2008，29（6）：71－76 （in Chinese）.［周寧，張衛(wèi)華.基于直接積分法的弓網(wǎng)耦合系統(tǒng)動態(tài)性能仿真分析 ［J］.中國鐵道科學(xué)，2008，29 （6）：71－76.］

［3］ZHAO Fei，LIU Zhigang，ZHANG Xiaoxiao.Simulation of high－speed pantograph－catenary system dynamic performance based on finite element model［J］.Journal of the China Railway Society，2012，34 （8）：33－38 （in Chinese）. ［趙飛，劉志剛，張曉曉.基于有限元的高速弓網(wǎng)系統(tǒng)動態(tài)性能仿真研究［J］.鐵道學(xué)報，2012，34 （8）：33－38.］

［4］Cho YH，Lee K.Influence of contact wire pre－sag on the dynamics of pantograph－railway catenary ［J］.International Journal of Mechanical Sciences，2010，52 （11）：1471－1490.

［5］CAI Guohua.The experimental technique of pantograph of super express train in low speed wind tunnel［J］.Journal of Railway Engineering Society，2006，4 （94）：67－70 （in Chinese）.［蔡國華.高速列車受電弓低速風(fēng)洞試驗技術(shù) ［J］.鐵道工程學(xué)報，2006，4 （94）：67－70.］

［6］ZHANG Jianhui，YANG Jiong，YAO Yong.Wind tunnel test investigation on pantograph drag reduction of high speed train［J］.Journal of Railway Science and Engineering，2010，7（6）：116－121 （in Chinese）.［張建輝，楊炯，姚勇.高速列車受電弓減阻的風(fēng)洞試驗研究 ［J］.鐵道科學(xué)與工程學(xué)報，2010，7 （6）：116－121.］

［7］Bocciolone M，Resta F.Pantograph aerodynamic effects on the pantograph－catenary interaction ［J］.Vehicle System Dynamics，2006，44 （S1）：560－570.

［8］YANG Zhen.Study on current collection of high－speed pantograph based on aerodynamic ［J］.Electric Railway，2009，3：17－20 （in Chinese）.［楊楨.基于空氣動力學(xué)的受電弓高速受流研究 ［J］.電氣化鐵道，2009，3：17－20.］

［9］TIAN Hongqi.Study on the characteristics of train air resistance under wind environment ［J］.China Railway Science，2008，29 （5）：108－112 （in Chinese）.［田紅旗.風(fēng)環(huán)境下的列車空氣阻力特性研究 ［J］.中國鐵道科學(xué)，2008，29 （5）：108－112.］

［10］ZHU Hongjun，LIN Yuanhua，XIE Longhan.Fluent12fluid analysis and engineering simulation ［M］.Beijing：Tsinghua University Press，2011 （in Chinese）.［朱紅鈞，林元華，謝龍漢.Fluent12流體分析及工程仿真 ［M］.北京：清華大學(xué)出版社，2011.］

［11］Zhou Ning，Zhang Weihua.Investigation on dynamic performance and parameter optimization design of pantograph and catenary system ［J］.Finite Elements in Analysis and Design，2011，47 （3）：288－295.

［12］Simon Walters.Simulation of fuzzy control applied to a railway pantograph－catenary system ［C］／／Proceedings of the 14th International Conference on Knowledge－based and Intelligent Information and Engineering Systems：Part II，2010：322－330.