武 闖,范晨光,2,郭春艷,代成棟,徐劍鋒
(1.西南交通大學(xué) 力學(xué)與航空航天學(xué)院,四川 成都 610031;2.西南交通大學(xué) 應(yīng)用力學(xué)與結(jié)構(gòu)安全四川省重點(diǎn)實(shí)驗(yàn)室,四川 成都 610031;3.西南交通大學(xué) 電氣工程學(xué)院,四川 成都 610031)
隨著我國經(jīng)濟(jì)的快速發(fā)展,公共交通己經(jīng)無法滿足城市的日??瓦\(yùn)要求,大多城市開始建設(shè)地鐵來緩解常規(guī)公共交通的壓力[1]。接觸網(wǎng)是城市軌道交通工程中的重要系統(tǒng)設(shè)備之一,它擔(dān)負(fù)著為地鐵輸送電能的重要作用,目前接觸網(wǎng)系統(tǒng)分為兩種:架空接觸網(wǎng)系統(tǒng)(柔性接觸網(wǎng))和供電軌系統(tǒng)。供電軌系統(tǒng)又分為架空供電軌(剛性接觸網(wǎng))和第三軌兩種。其中,架空接觸網(wǎng)和架空供電軌由弓網(wǎng)系統(tǒng)供電;第三軌通過靴軌系統(tǒng)供電。接觸線高度是指軌道平面與接觸線下表面的垂直距離,是表征接觸網(wǎng)狀態(tài)的重要靜態(tài)幾何參數(shù)之一[2]。
由于剛性接觸網(wǎng)具有空間占用率低、免維護(hù)、安全可靠性高等特性,廣泛應(yīng)用于地鐵線路段,比如北京、上海、深圳等多個(gè)城市的地鐵線路正線區(qū)段。雖然柔性接觸網(wǎng)架設(shè)需要極大的空間,其結(jié)構(gòu)復(fù)雜,較剛性接觸網(wǎng)更易斷裂,但由于其具有良好的彈性,可以支持列車高速度運(yùn)行,因此,地鐵線路段往往也會出現(xiàn)大量的柔性接觸網(wǎng)區(qū)段,比如上海地鐵2、3、11號線,深圳地鐵11號線和成都地鐵18號線等。對于地鐵線路接觸網(wǎng)系統(tǒng),柔性接觸網(wǎng)段和剛性接觸網(wǎng)段是共同存在的[1]。受電弓作為動車組從接觸網(wǎng)獲取電能的重要電氣設(shè)備,工作環(huán)境非常惡劣,長期受到外部氣動載荷、受電弓與接觸網(wǎng)之間的振動及摩擦等隨機(jī)因素的影響,容易發(fā)生疲勞破壞。
近年來國內(nèi)外諸多學(xué)者著力于探究地鐵受電弓的疲勞壽命預(yù)測問題。文獻(xiàn)[1—7]建立了受電弓三維有限元模型,通過仿真計(jì)算得到應(yīng)力譜,并對受電弓部件進(jìn)行壽命預(yù)測;陳花麗[6]考慮了弓網(wǎng)耦合作用下受電弓參數(shù)改變對接觸網(wǎng)壽命的影響;黃思俊[7]等利用HyperMesh建立受電弓的有限元模型,將疲勞壽命分析方法運(yùn)用到受電弓系統(tǒng)的隨機(jī)振動中,評估隨機(jī)振動載荷作用下高速受電弓的疲勞強(qiáng)度。以上研究多是基于仿真分析得到受電弓的應(yīng)力譜,然后進(jìn)行疲勞壽命預(yù)測,對于受電弓經(jīng)過剛?cè)峤佑|網(wǎng)過渡段時(shí)(即導(dǎo)線高度突變時(shí))部分部件的應(yīng)力發(fā)生跳變的現(xiàn)象關(guān)注較少。本文采用現(xiàn)場實(shí)測的應(yīng)力時(shí)程數(shù)據(jù)對受電弓的關(guān)鍵部件進(jìn)行疲勞可靠性分析,并且著重關(guān)注導(dǎo)線高度突變部分對受電弓疲勞壽命的貢獻(xiàn),為后續(xù)對地鐵受電弓的安全性評估提供參考。
受電弓屬于軌道交通車輛受流裝置,安裝在車輛頂部,與供電網(wǎng)接觸,集取接觸網(wǎng)電流,供地鐵車輛運(yùn)行使用。受電弓在正常運(yùn)行過程中,不斷與接觸網(wǎng)產(chǎn)生摩擦,直接受到弓網(wǎng)接觸交變載荷作用;尤其是經(jīng)過接觸網(wǎng)剛?cè)徇^渡段時(shí),由于導(dǎo)線高度的劇烈變化,會產(chǎn)生瞬間的沖擊作用,引起結(jié)構(gòu)較強(qiáng)的振動響應(yīng),應(yīng)力幅值也相應(yīng)地產(chǎn)生較大的變化。某型號受電弓[8]運(yùn)行里程在50.4~65.5萬km時(shí),上框架處產(chǎn)生了裂紋,如圖1所示。
圖1 某型號受電弓上框架焊趾處裂紋[8]
本文基于國內(nèi)某地鐵線路的弓網(wǎng)線路試驗(yàn),分析兩種不同受電弓的載荷特性。通過線路試驗(yàn)測試接觸網(wǎng)線路導(dǎo)線高度變化、各部件的應(yīng)力變化等。在受電弓關(guān)鍵部位粘貼應(yīng)變片,如圖2所示。受電弓升起后由于受到各種外力,張貼在其表面的應(yīng)變片阻值發(fā)生改變,通過應(yīng)變采集系統(tǒng)計(jì)算出相應(yīng)的應(yīng)變數(shù)值。應(yīng)變信號采集使用美國NI-9235C系列應(yīng)變/橋輸入模塊,通過高壓采集箱將信號引入車內(nèi)。在受電弓開口方向前端安裝滑板動態(tài)高度檢測相機(jī),如圖3所示,用來檢測導(dǎo)線高度變化。地鐵車內(nèi)配電柜提供220 V交流電,此外,為了防止車輛意外斷電的情況,通過穩(wěn)壓器和UPS給檢測設(shè)備供電。
圖2 局部應(yīng)變片布置
圖3 受電弓導(dǎo)線高度參數(shù)測量裝置
試驗(yàn)工況為地鐵正線ATO模式(列車自動駕駛系統(tǒng),能夠?qū)崿F(xiàn)列車自動行駛、精確停車、站臺自動化作業(yè)、無人折返、列車自動運(yùn)行調(diào)整等功能),線路全長51 km,弓網(wǎng)靜態(tài)接觸力為120 N。為對比兩種受電弓性能差異,兩種型號受電弓均采用了相同的碳滑板。在測試完成后,將碳滑板拆除,安裝在下一次測試的受電弓弓頭上,以確保試驗(yàn)遵循單一變量原則。分析碳滑板高度變化曲線,地鐵線路由柔性接觸網(wǎng)和剛性接觸網(wǎng)混合編制而成,剛性網(wǎng)局部導(dǎo)線高度出現(xiàn)較小的跳變,剛性網(wǎng)與柔性網(wǎng)交界處導(dǎo)線高度出現(xiàn)了較大跳變。其中剛性接觸網(wǎng)區(qū)段占比較大。如圖4所示,剛性接觸網(wǎng)區(qū)段,接觸網(wǎng)導(dǎo)線高度有4 040 mm和4 200 mm兩種,而柔性接觸網(wǎng)區(qū)段,接觸網(wǎng)導(dǎo)線高度為4 640 mm。
圖4 全程導(dǎo)線高度時(shí)程曲線
首先對受電弓進(jìn)行有限元靜力學(xué)分析得到各部分的應(yīng)力分布,找到應(yīng)力較大的部位,由此來確定應(yīng)變測點(diǎn)的位置。根據(jù)有限元靜力學(xué)仿真結(jié)果,應(yīng)力較大的部位主要集中在受電弓上框架、下臂桿和底架處,在這些位置的測點(diǎn)布置了應(yīng)變片;同時(shí)還在弓頭、拉桿、平衡桿處布置了測點(diǎn),測點(diǎn)布置如圖5、圖6所示。由此可以全面準(zhǔn)確得到整個(gè)受電弓的應(yīng)力分布狀態(tài)。材料參數(shù)如表1所示。試驗(yàn)表明,當(dāng)受電弓經(jīng)過接觸網(wǎng)導(dǎo)線高度變化區(qū)段時(shí),由于接觸網(wǎng)導(dǎo)線高度的變化引起的應(yīng)力突變部位主要位于底架、上框架及其附屬連接部分。
表1 受電弓各部件參數(shù)統(tǒng)計(jì)
圖5 A型受電弓測點(diǎn)布置圖
圖6 B型受電弓測點(diǎn)布置圖
選取A型受電弓的測點(diǎn)5(拉桿連接件)處,工況選取為前弓ATO模式;B型受電弓的測點(diǎn)3(底架)處,工況選取為后弓ATO模式。兩測點(diǎn)的應(yīng)力時(shí)程曲線如圖7、圖8所示,從圖中可以看出,應(yīng)力突變發(fā)生時(shí)間與全程導(dǎo)線高度時(shí)程曲線的突變位置基本吻合,由此判斷應(yīng)力突變是由導(dǎo)線高度變化引起的;并且在導(dǎo)線高度不變的穩(wěn)定段,應(yīng)力波動不大。
圖7 A型受電弓測點(diǎn)5應(yīng)力時(shí)程曲線
圖8 B型受電弓測點(diǎn)3應(yīng)力時(shí)程曲線
地鐵正線運(yùn)行時(shí),受電弓與接觸網(wǎng)直接接觸,由于接觸網(wǎng)性質(zhì)的改變引起較大應(yīng)力幅值的突變并且這種突變是周期性的,這種周期性的交變載荷可能會導(dǎo)致受電弓各部件產(chǎn)生疲勞破壞,危及行車安全,嚴(yán)重情況下可能導(dǎo)致線路停運(yùn)。因此有必要對受電弓進(jìn)行疲勞壽命預(yù)測以確保弓網(wǎng)系統(tǒng)正常運(yùn)行。
受電弓所受到的外界載荷為隨機(jī)交變載荷。疲勞載荷的描述中使用應(yīng)力幅值σa和應(yīng)力范圍Δσ以及平均應(yīng)力σm表示[9]:
σa=1/2(σmax-σmin)
(1)
Δσ=σmax-σmin
(2)
σm=1/2(σmax+σmin)
(3)
本文對受電弓的疲勞壽命進(jìn)行預(yù)測,具體方法如下:
(1) 對采集到的應(yīng)變數(shù)據(jù)進(jìn)行處理(σ=Eε),得到受電弓各部件動應(yīng)力時(shí)程曲線;
(2) 通過雙參數(shù)雨流計(jì)數(shù)法實(shí)現(xiàn)對應(yīng)力幅值和均值的統(tǒng)計(jì),得到非對稱二維應(yīng)力譜;
(3) 考慮雨流計(jì)數(shù)統(tǒng)計(jì)的具體情況,對材料的S-N曲線修正,得到構(gòu)件的P-S-N曲線;
(4) 采用“波動中心法”得到非對稱一維應(yīng)力譜并經(jīng)過Goodman直線修正,得到其8級對稱應(yīng)力譜;
(5) 基于結(jié)構(gòu)8級應(yīng)力譜統(tǒng)計(jì)結(jié)果,結(jié)合材料P-S-N曲線和Miner線性累積損傷理論,對受電弓關(guān)鍵部件進(jìn)行疲勞壽命預(yù)測。
總體技術(shù)路線如圖9所示。
圖9 技術(shù)路線
雨流計(jì)數(shù)法是20世紀(jì)50年代由英國的兩位工程師MATSUISHI和ENDO提出來的[10]。為保證載荷譜塊的完整性,將前弓與后弓的應(yīng)力測試數(shù)據(jù)拼接起來作為完整的應(yīng)力時(shí)程曲線進(jìn)行統(tǒng)計(jì)。
文章采用雨流計(jì)數(shù)法對動應(yīng)力時(shí)程數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析,基于軟件編程實(shí)現(xiàn)[11-13]。具體步驟如下:
(1) 對應(yīng)力時(shí)程信號再造,使其只包含峰谷值。原始信號應(yīng)定義為Ci,同時(shí)令A(yù)i=Ci,Bi=Ci;依次判斷Ai數(shù)據(jù),如果不滿足峰谷值判斷條件,記相應(yīng)下標(biāo)Bi為NaN并刪除元素,同時(shí)實(shí)現(xiàn)等數(shù)值壓縮,得到處理后新載荷時(shí)間歷程記作Bi。
(2) 用四峰谷值法對Bi進(jìn)行一階段雨流計(jì)數(shù)。定義兩個(gè)空數(shù)列F、J用以存儲應(yīng)力循環(huán)的幅值和均值。當(dāng)length(Bi)>3,令前兩個(gè)點(diǎn)的差值c=0(首次計(jì)數(shù)人為給定),計(jì)算a1=abs(Bj+1-Bj),a2=abs(Bj+1-Bj+2),m3=(Bj+1+Bj)/2。
如果滿足應(yīng)力循環(huán)判斷條件a1≤a2且a1≤c,則提取循環(huán)幅值a1/2,均值m3,分別記錄在數(shù)列F、J中,同時(shí)刪去Bj和Bj+1兩點(diǎn)。為了存儲提取的載荷點(diǎn),將定義數(shù)列D,依次將Bj、Bj+1、Bj填入,用NaN與下一個(gè)循環(huán)隔開。此時(shí)再令c=a1,繼續(xù)循環(huán)。
(3) 一階段雨流計(jì)數(shù)結(jié)束后,剩下的發(fā)散-收斂波形無法構(gòu)成循環(huán),尋找最高波峰(最低波谷)點(diǎn),從該點(diǎn)處截?cái)?左右兩應(yīng)力波段對換位置首尾相接,形成收斂-發(fā)散波形,返回步驟(2),繼續(xù)循環(huán)。
(4) 直到length(Bi)<3,跳出循環(huán)。記錄最后一個(gè)應(yīng)力循環(huán),此時(shí)幅值為(Bmax-Bmin)/2,均值為(Bmax+Bmin)/2,程序結(jié)束。
數(shù)據(jù)處理步驟如圖10所示。
圖10 雨流計(jì)數(shù)程序流程圖
采用如下公式估算受電弓材料在指定概率P下的疲勞壽命[1]:
(lgN)P=AP+BPlgS
(4)
式中:AP和BP分別為截距參數(shù)和斜率參數(shù);N為疲勞壽命;S為應(yīng)力幅值;P為存活概率。
傳統(tǒng)的疲勞壽命理論認(rèn)為當(dāng)應(yīng)力幅值低于材料疲勞極限時(shí),結(jié)構(gòu)不會出現(xiàn)疲勞結(jié)構(gòu)破壞,認(rèn)為此時(shí)結(jié)構(gòu)為無限壽命。但是在工程實(shí)際中,受電弓結(jié)構(gòu)在低應(yīng)力幅值工作下,其結(jié)構(gòu)依然會出現(xiàn)疲勞損傷問題。因此,采用2m-1指數(shù)法外推得到材料的P-S-N曲線。計(jì)算公式如下:
(lgN)P=EP+FPlgS
(5)
不同存活概率下的外推P-S-N曲線參數(shù)取值如圖11、圖12所示。
圖11 鋁合金材料P-S-N曲線參數(shù)
圖12 不銹鋼材料P-S-N曲線參數(shù)
應(yīng)力時(shí)間歷程經(jīng)雨流循環(huán)計(jì)數(shù)處理之后得到的輸出數(shù)據(jù)文件為整個(gè)歷程中各次應(yīng)力循環(huán)的均值和幅值。由于數(shù)據(jù)量大,不便于統(tǒng)計(jì),因此有必要對數(shù)據(jù)進(jìn)行整理,編制出一維應(yīng)力譜。
采用“波動中心法”,只考慮應(yīng)力循環(huán)的幅值而不考慮應(yīng)力循環(huán)的均值,將應(yīng)力譜簡化為一元隨機(jī)變量,具體方法如下[14]:
(1) 檢索均值數(shù)據(jù),計(jì)算所有應(yīng)力循環(huán)的總平均值作為波動中心。
(6)
式中:M為波動中心;Mi為第i級應(yīng)力均值;ni為第i級應(yīng)力均值對應(yīng)的頻次。
求得波動中心后,將幅值作為應(yīng)力循環(huán)的動應(yīng)力分量疊加于波動中心之上。
(2) 檢索幅值數(shù)據(jù),求得應(yīng)力幅值最大值和最小值的差值Δ。由于實(shí)測應(yīng)力波數(shù)據(jù)的特殊性,幅值統(tǒng)計(jì)數(shù)據(jù)呈現(xiàn)兩端多,中間少的特點(diǎn),因此采用非等距數(shù)值進(jìn)行分組[15]。
(7)
式中:Ri為組間距;N為應(yīng)力幅值級數(shù),i=1,2…N-1。
(3) 依次讀取幅值數(shù)據(jù),判斷應(yīng)力幅值屬于哪一級應(yīng)力范圍之內(nèi),存儲于新建變量之中。為統(tǒng)計(jì)方便,各級應(yīng)力幅值組采用其組中值來表示。
由于測點(diǎn)數(shù)較多,選取受電弓關(guān)鍵部位的一維應(yīng)力譜,如表2、表3所示。從表中可以看出,低應(yīng)力幅值占總計(jì)數(shù)的絕大多數(shù),這是由于受電弓在恒定工作高度運(yùn)行穩(wěn)定的結(jié)果;高應(yīng)力幅值占比較小,是由于導(dǎo)線高度突變次數(shù)較少產(chǎn)生的結(jié)果。
表2 A型受電弓各部件8級應(yīng)力譜
表3 B型受電弓各部件8級應(yīng)力譜
統(tǒng)計(jì)分析得到8級非對稱應(yīng)力譜,為保守估計(jì),基于Goodman方程對應(yīng)力譜進(jìn)行修正,將其轉(zhuǎn)化為對稱循環(huán)應(yīng)力譜,修正公式如下[9]:
(8)
式中:σN為對稱循環(huán)應(yīng)力幅值;σa,σm為應(yīng)力循環(huán)的幅值和均值;σu為材料極限強(qiáng)度。
文獻(xiàn)[16]詳細(xì)介紹了線性、非線性、兩級線性疲勞累積損傷理論。Miner線性累積損傷理論是工程上廣泛采用的一種疲勞壽命計(jì)算方法。
經(jīng)典Miner線性累積損傷理論認(rèn)為,各級應(yīng)力造成的損傷可線性疊加,當(dāng)總損傷量D達(dá)到1時(shí)結(jié)構(gòu)發(fā)生破壞。對于一個(gè)8級載荷譜而言,累積的總損傷為:
(9)
式中:ni為第i級應(yīng)力水平下的循環(huán)次數(shù);Ni為第i級應(yīng)力水平下的壽命。實(shí)際運(yùn)營時(shí)考慮安全與可靠性,D一般取0.7[17]。
由實(shí)測數(shù)據(jù),該8級應(yīng)力譜對應(yīng)的測試?yán)锍蘈1為102 km,根據(jù)Miner線性累積損傷理論,可知結(jié)構(gòu)在存活概率P下運(yùn)行的總里程L為:
L=L1D/D8
(10)
統(tǒng)計(jì)結(jié)果如表4、表5所示,可以看出,隨著可靠度的增加,受電弓各部件的疲勞壽命逐漸減小;對于A型受電弓,上框架處壽命最短,容易產(chǎn)生裂紋,應(yīng)關(guān)注;下臂桿和弓頭連接件壽命相近,底架、拉桿處均為長壽命;對于B型受電弓,同樣是上框架的壽命最短,應(yīng)關(guān)注;底架次之,下臂桿、平衡桿、拉桿處均為長壽命。
表4 A型受電弓ATO工況各部件壽命預(yù)測 萬km
表5 B型受電弓ATO工況各部件壽命預(yù)測 萬km
參照文獻(xiàn)[17]的5.3小節(jié),受電弓的設(shè)計(jì)壽命里程不應(yīng)少于1.2×103萬km。A型受電弓的預(yù)測運(yùn)行里程為7.72×103萬km,B型受電弓的預(yù)測運(yùn)行里程為3.22×103萬km,由此可見,兩種型號的受電弓在ATO工況運(yùn)行狀態(tài)下,各部件均滿足壽命設(shè)計(jì)要求,但B型受電弓的壽命較短。
基于A、B兩種不同型號的受電弓現(xiàn)場試驗(yàn)測試結(jié)果,得到受電弓全程ATO工況的動應(yīng)力載荷譜;應(yīng)用雨流計(jì)數(shù)遞歸算法實(shí)現(xiàn)雨流計(jì)數(shù)程序設(shè)計(jì)并對受電弓部件進(jìn)行疲勞可靠性分析,得到以下結(jié)論:
(1) 在導(dǎo)線高度不變的情況下,速度等級變化產(chǎn)生的交變應(yīng)力幅值很小,對受電弓壽命的貢獻(xiàn)幾乎可以忽略不計(jì);導(dǎo)線高度變化時(shí),會對受電弓產(chǎn)生沖擊因而出現(xiàn)較大的應(yīng)力幅值;
(2) 對于改進(jìn)后兩種型號的受電弓來說,最短壽命預(yù)測結(jié)果均為上框架部位,且B型受電弓的壽命較短,但都滿足壽命設(shè)計(jì)要求;
(3) 在搭建接觸網(wǎng)系統(tǒng)時(shí),應(yīng)盡量減小接觸網(wǎng)之間落差以降低應(yīng)力突變,提高受電弓的疲勞壽命。