列車弓網接觸壓力仿真與受流性分析

0高速弓網關系的研究方法當列車在低速運行時，由于發射弓高度的變化，速度減慢，并且受電弓和接觸線之間的接觸壓力幾乎可以被視為一個常量。然而，隨著列車速度的推移，隨著列車速度的推移，弓網之間的接觸壓力非常復雜。因此,高速受流問題(弓網關系)就變得非常突出。法國、德國和日本等高速鐵路發達國家將弓網關系置于輪軌關系同等重要的位置,一直在從事這方面的研究工作,并已取得了豐碩的成果。在指導弓網結構類型,系統參數合理匹配以及保證列車高速安全運行方面起到了重要作用。而至今我國對高速弓網關系的研究還不夠深入,因此,從技術和經濟的角度出發,都體現了深入研究高速弓網關系的重要性和緊迫性。本文針對接觸網不平順分析的需求,通過利用Newmark算法仿真弓網耦合動力學模型,得到弓網間的接觸壓力樣本數據。論文提出采用譜分析方法研究接觸網的不平順性,具體采用功率譜密度來分析影響弓網受流的波長成分。仿真結果表明,該方法具有較好的直觀性和有效性。1射箭耦合系統模型1.1振動微分方程的選取接觸網懸掛是一個復雜的機械系統,除了有接觸線和承力索外,還有吊弦、線夾、定位器、腕臂、電連接器、分段絕緣器等。本文接觸網選擇簡單鏈形懸掛方式,將承力索和接觸線的質量看成是多個集中質點的分布。對于接觸線,在吊弦點之間為多個集中質點;對于承力索,僅考慮吊弦點和吊弦中間點。由于接觸網的橫向振動對弓網接觸的影響很小,可忽略,因此,只討論接觸網的垂向振動。承力索的振動微分方程為接觸線的振動微分方程為式中,ρA,ρB分別為承力索和接觸線的線密度;EIA,EIB分別為承力索和接觸線的彈性模量;SA,SB分別為承力索和接觸線的張力;m為振動階數,Am,Bm分別為承力索和接觸線第m階幅值;mTA,mTB,mD分別為腕臂、定位器和吊弦的質量;p,w分別為吊弦的個數和支撐點的個數;kD,kTA分別為吊弦和腕臂的剛度;L為跨距長度。1.2受電弓的線性化模型受電弓是由框架(下臂桿、推桿、上臂桿、平衡桿)、弓頭、底架和傳動機構等4部分組成的,結構較復雜。受電弓原本是一個非線性系統,但在其工作點位置附近可作線性化處理。計算受電弓與接觸線之間的接觸壓力時,不考慮機車的振動,把受電弓的基座看成固定不變,則受電弓的線性化模型的微分方程如下:受電弓弓頭運動微分方程受電弓框架運動微分方程為式中,M1,K1,C1分別為受電弓弓頭歸算質量、剛度系數和阻尼系數;M2,K2,C2分別為受電弓框架歸算質量、剛度系數和阻尼系數;y1,y2分別為受電弓弓頭和框架垂直位移;F0為受電弓抬升力。1.3弓網耦合系統弓網耦合系統模型為一個二元受電弓線性模型與簡單鏈形懸掛接觸網通過KS耦合在一起,形成弓網耦合系統模型(圖略)。聯立承力索、接觸線和受電弓弓頭、框架的運動微分方程就可得到弓網耦合系統的運動微分方程組。對接觸網來說,如果位移級數取N項,則整個接觸網的廣義自由度將是2N個,加上受電弓的歸算質量M1、M2的自由度y1、y2,整個弓網系統就有2N+2個自由度,其2N+2階矩陣表達式為受電弓與接觸線間的接觸壓力FC通過接觸剛度KS在接觸點處滑板與接觸線上的位移求得,即式中,y1為受電弓弓頭位移;yB(xc,t)為受電弓運行到xc處時,接觸點處接觸線的位移,即至此就完成了建立接觸網、受電弓及其耦合系統的數學模型。2解析法的精確解求解振動微分方程的方法可分為解析法和數值法2大類。由于數學上的困難,通常只有某些簡單的問題才能通過解析法得到精確的解答,而對于多數復雜結構問題,還不能得到精確的解。由于弓網系統的復雜性,進行弓網耦合振動分析已遠遠超出了理論解析的范疇,因此,本文采用Newmark數值方法。2.1節點位移矩陣的確定Newmark法是由加速度法引伸而來的,其假設的速度、位移與加速度的關系如下:式中,α,β為參數,由積分精度和穩定性決定。Newmark法通過t+Δt的平衡方法求解節點位移矩陣Yt+Δt,因此,式(5)可寫成:由式(9)得:將式(8)代入式(11),然后再一并代入式(10),則可得到Yt+Δt:式中,Yt+Δt為t+Δt時刻節點位移矩陣,Ft+Δt為t+Δt時刻節點外激振力列陣,Δt為時間步長。則式(12)可寫為2.2剛度矩陣的計算由上述算法分析,得出Newmark計算弓網耦合接觸壓力和抬升量的步驟如下:(1)初始化。a.形成弓網耦合系統的剛度矩陣F、質量矩陣M和阻尼矩陣;b.給定初始位移、初始速度和初始加速度;c.選擇時間步長Δt和參數α,β。(2)積分計算。對于每一步長Δt:a.形成剛度矩陣K,計算有效剛度矩陣進行三角分解;c.計算t+Δt時刻的有效荷載計算t+Δt時刻的位移矩陣Yt+Δt;e.計算t+Δt時刻的速度和加速度矩陣;f.計算t+Δt時刻的接觸壓力和接觸抬升量。具體流程如圖1所示。3受電弓升降控制本文利用Newmark法編寫弓網耦合仿真程序,分別對車速為80km/h和120km/h的弓網受流情況進行仿真。接觸線跨距為65m,受電弓的抬升力為70N時,接觸壓力和接觸線動態抬升量隨受電弓運動位置的變化情況如圖2、圖3所示。從仿真結果可以看出,接觸壓力和接觸線抬升量以跨距為周期變化,在跨距內接觸壓力和接觸線抬升量在定位器和吊弦處出現明顯變化,接觸壓力的峰值出現在定位器處,而接觸線抬升量出現在跨中。這是因為:受電弓抬升力不變,且接觸網在定位器處的彈性剛度大,在跨中處的彈性剛度小。4接觸網不順流特性接觸網的不平順是指受電弓取流的接觸面沿接觸線延長與理論上平順的接觸面的偏差。造成接觸網不平順的原因有:接觸線卷繞過程中形成周期性的折痕;定位器在長期運營過程中產生的累積變形等。這些不平順激發了受電弓和接觸網之間的有害振動,惡化受流質量,對列車的安全運行產生不利的影響。從前面的仿真結果可以看出,接觸網的不平順直接客觀地反應在弓網接觸壓力上,因此,可以利用弓網接觸壓力研究接觸網的不平順。接觸網不平順與軌道不平順類似,是與線路里程有關的隨機過程,是由不同波長、不同相位、不同幅值的隨機不平順波疊加而成的。接觸網不平順的隨機特性決定了對接觸網不平順的描述不能用一個明確的數學表達式表示,而只能用隨機振動理論中描述隨機數據的“均方差”、“方差”、和“功率譜密度函數”等統計函數來表示接觸網不平順的特征,從時域、頻域對接觸網不平順的幅值特性、波長結構以及是否包含周期性波形作全面的描述。功率譜密度函數是描述作為平穩隨機過程的接觸表面不平順的最重要和最常用的統計函數。工程中常用功率譜描述譜密度對頻率的函數變化。接觸網不平順功率譜圖是以譜密度為縱坐標和以頻率或波長為橫坐標的連續變化曲線,它可清楚地表明不平順的大小隨頻率或波長的變化關系。本文通過功率譜方法,獲得列車速度為80km/h和120km/h時的接觸網不平順功率譜圖,如圖4、圖5所示。從功率譜圖可以看出,圖5顯示出了1,5m波長的特征不平順,圖6顯示出了0.5,1.2,1.5m波長的特征不平順。如果當接觸網不平順的波長為5m時,列車運行速度為80km/h,振動頻率約為4.45Hz。當接觸網本身的固有振動頻率為4.45Hz時,就會發生弓網共振,受電弓與接觸網保持接觸的振幅最低,受流質量惡化,離線的幾率增加。因此,必須通過檢測消除波長為5m的接觸網不平順特征,保證列車安全運行。5吊弦長度的確定接觸網彈性懸掛的不平順是研究受電弓與接觸網系統振動的基礎,接觸網的不平順也將影響弓網受流質量。本文從弓網耦合系統模型出發,首先仿真了列車運行速度為80km/h和120km/h的弓網動態受流情況,得到了接觸壓力和接觸線抬升量變化趨勢。然后,通過建立接觸網不平順功率譜圖,分析了接觸網不平順情況,從功率譜圖可以清晰地看出接觸網不平順波長。接觸網不平順的研究在我國還需要加大力度,爭取早日建立我國的接觸網不平順功率密度函數,這將對我國電氣化鐵道交通的發展起到積極的推動作用。那么A、B處的實際張力為2吊弦計算修正設全補償簡單鏈形懸掛結構高度Hj,懸掛自重g,接觸線自重gj,跨距l,第1吊弦點至懸掛點距離D,承力索張力TC,接觸線張力TJ,接觸線預留弛度F0,則x點處吊弦長度計算公式為其中,承力索張力TC在中間柱處需要增加半個錨段張力差ΔTR,錨段關節處需根據式(16)、式(11)和式(17)增加張力差。3工程實例達成線南充東—南充區間電氣化擴能改造大部分是利舊工程,小曲線半徑線路較多。其懸掛類型為THJ-95+CTHA-120直鏈形懸掛,懸掛自重為22.78kN/m,錨段長度為1400m,跨距為45m,半錨段跨距數15個,曲線半徑600m,平腕臂長2.8m,承力索膨脹系數17×10-6,額定張力15kN,溫差20℃,根據式(8)得張力差為0.645kN;錨段關節內由式(16