謝素明,薛寧鑫,馬夢琳( .大連交通大學(xué)交通運輸工程學(xué)院,遼寧大連608; .中國北車集團長春軌道客車股份有限公司技術(shù)中心,吉林長春3006)*
高速動車組車體模態(tài)分析建模方法及試驗驗證
謝素明1,薛寧鑫1,馬夢琳2
( 1.大連交通大學(xué)交通運輸工程學(xué)院,遼寧大連116028; 2.中國北車集團長春軌道客車股份有限公司技術(shù)中心,吉林長春130062)*
基于“剛度等效”原則,建立某自主研發(fā)的高速動車組車體剛結(jié)構(gòu)有限元模型,車體一階垂向彎曲和一階扭轉(zhuǎn)的振動頻率的數(shù)值解與模態(tài)試驗值的誤差分別為2.75%和7.90%;在此基礎(chǔ)上,以整備狀態(tài)車體有限元模型重心與實際重心一致為質(zhì)量分布原則,創(chuàng)建五種整備狀態(tài)車體模態(tài)分析模型.模態(tài)分析結(jié)果表明:與實際重心最接近的模型5的計算結(jié)果與試驗結(jié)果最為接近,誤差分別為0.03%和6.85%;建議車體方案設(shè)計階段采用模型2的建模方法估算整備狀態(tài)車體模態(tài).
高速動車組;整備狀態(tài)車體;建模方法;模態(tài)分析
近年來,關(guān)于軌道客車模態(tài)問題的研究較多[1],郝魯波建立了某25G型空調(diào)客車車體分析模型,詳細(xì)分析了附件質(zhì)量和剛度對車體模態(tài)的影響[2];曲天威等分析了國內(nèi)外機車車體和轉(zhuǎn)向架模態(tài)分析現(xiàn)狀,對其中存在的問題和評價方法存在的差異進(jìn)行了討論[3];趙陽陽等建立了高速列車CRH3承載結(jié)構(gòu)和整備狀態(tài)有限元模型,對比分析了車體下吊分別采用彈性連接和剛性連接時車體模態(tài)頻率的變化[4];宮島等建立高速列車整備狀態(tài)車體有限元模型及包含車體彈性的高速列車車輛剛?cè)狁詈蟿恿W(xué)模型,分析了車下設(shè)備吊掛方式對車體整備狀態(tài)模態(tài)參數(shù)的影響,提出車下設(shè)備隔振橡膠件參數(shù)設(shè)計方法[5-6].
目前,我國自主研發(fā)的高速動車組車體設(shè)計遵循的《200 km/h及以上速度級鐵道車輛強度設(shè)計及試驗鑒定暫行規(guī)定》中要求:整備狀態(tài)下車體一階垂向彎曲振動頻率不得低于10 Hz.同時,TB1335-1996中也規(guī)定了車體垂向彎曲剛度與扭轉(zhuǎn)剛度的限定值.由于車體整體剛度和質(zhì)量分布決定著整備狀態(tài)車體的振動品質(zhì),局部結(jié)構(gòu)改變的影響可以忽略不計,因此車體設(shè)計階段如何較為準(zhǔn)確地分析出各種質(zhì)量狀態(tài)下車體模態(tài)至關(guān)重要.
本文以某自主研發(fā)的動力分散型高速動車組車體為研究對象,采用數(shù)值分析與試驗相結(jié)合的方法,建立了各種質(zhì)量分布狀態(tài)下車體模態(tài)分析模型,并將其計算結(jié)果與試驗結(jié)果進(jìn)行對比,總結(jié)出車體設(shè)計不同階段的模態(tài)分析模型的建模方法.
高速動車組車體主要由EN AW-6005A和EN AW-6082的鋁合金中空擠壓型材、板材通過插接、焊接組裝成底架、側(cè)墻、端墻、車頂?shù)却蟛考唇映蔀檐圀w,車體長度為24.5 m,寬度為3.3 m,高度為3.9 m.
車體結(jié)構(gòu)質(zhì)量定義為車體承載結(jié)構(gòu)的質(zhì)量,為9.4 t.整備狀態(tài)車體質(zhì)量,亦即空車質(zhì)量,是指動車組車輛按出廠技術(shù)條件將所有的設(shè)備、內(nèi)飾等裝備完整之后的質(zhì)量,為37.8 t.車體整備質(zhì)量通過進(jìn)行預(yù)組裝和總組裝完成,預(yù)組裝主要進(jìn)行塞拉門、車窗、隔熱材安裝,地板、衛(wèi)生間預(yù)置,鋁地板以上布線,車底線排布置,空調(diào)機組、車內(nèi)回風(fēng)道、水箱安裝;車頂安裝受電弓、主斷等設(shè)施;總組裝主要進(jìn)行車內(nèi)外設(shè)備裝配,包括內(nèi)裝(側(cè)墻板、行李架、客室座椅等),車內(nèi)設(shè)備,車下吊裝等,圖1為整備狀態(tài)車體示意圖.
圖1 高速動車組整備狀態(tài)車體示意圖
車體剛結(jié)構(gòu)低階整體模態(tài)主要取決于結(jié)構(gòu)整體剛度.建立高速動車組車體剛結(jié)構(gòu)模態(tài)數(shù)值分析模型的關(guān)鍵是如何模擬車體中空擠壓鋁型材之間的焊接關(guān)系,參見圖2中的A和B及C處.
圖2中所示的中空擠壓鋁型材之間的焊接形式A搭接和B與C對接,均是沿車長方向的縱向長焊縫,應(yīng)基于“剛度等效”原則進(jìn)行有限元建模.以對接焊接形式B與C為例,解釋“剛度等效”原則.首先,分別創(chuàng)建兩個模型,模型1將對型材對接處離散為三維實體單元,并在型材可能接觸部位建立接觸關(guān)系;模型2將型材對接處離散為薄殼單元,設(shè)單元厚度為t;其次,分別對兩個模型施加四邊簡支約束、型材上表面施加相同垂向載荷;最后,根據(jù)垂向最大位移相等,確定對接處的等效厚度.高速動車組車體結(jié)構(gòu)有限元模型中單元總數(shù)為2 827 557,結(jié)點總數(shù)為2 345 639,車體的有限元模型如圖3所示.
圖2 高速動車組斷面型材結(jié)構(gòu)示意圖
圖3 車體結(jié)構(gòu)有限元模型
經(jīng)自由模態(tài)數(shù)值分析,高速動車組車體結(jié)構(gòu)的一階垂向彎曲振動頻率為16.307 Hz;一階扭轉(zhuǎn)振動頻率為24.05 Hz,振型云圖如圖4所示.車體一階垂向彎曲和一階扭轉(zhuǎn)的振動頻率試驗值分別為15.87和22.29 Hz[7],計算值與試驗值的誤差分別為2.75%與7.90%.可以說,車體結(jié)構(gòu)有限元模型能夠準(zhǔn)確地模擬車體的整體垂向和扭轉(zhuǎn)剛度.
動車組車體設(shè)計詳細(xì)階段才能給出整備質(zhì)量的詳細(xì)分布和重心位置,在方案設(shè)計階段僅知道車體整備質(zhì)量的大小.因此,在車體設(shè)計階段準(zhǔn)確地模擬出整備狀態(tài)車體模態(tài)極為困難.
在上述車體結(jié)構(gòu)模態(tài)分析模型的基礎(chǔ)上,建立五種整備狀態(tài)車體模態(tài)分析模型.模型1將200 kg以上的空調(diào)、受電弓和車下設(shè)備以集中質(zhì)量的形式布置在相應(yīng)安裝位置,再將其余質(zhì)量( 19.567 t)均布在地板上表面;模型2將模型1中的其余質(zhì)量( 19.567 t)以調(diào)整材料密度的形式均布在整車結(jié)構(gòu)上;模型3將車體質(zhì)量進(jìn)一步細(xì)化,將質(zhì)量大于50 kg的設(shè)備以集中質(zhì)量的形式布置在相應(yīng)安裝位置,再將其余質(zhì)量( 5.927 t)以質(zhì)量單元的形式均布在地板上表面;模型4在模型3的基礎(chǔ)上,將5.927 t質(zhì)量以調(diào)整材料密度的形式均布在整車結(jié)構(gòu)上;模型5增加由上、下兩層鋁板和中間的PVC泡沫構(gòu)成的類似三明治形式的復(fù)合地板的有限元模型,以考慮復(fù)合地板與底架螺栓聯(lián)接后對整車剛度的貢獻(xiàn).并進(jìn)一步細(xì)化車體質(zhì)量,全部考慮車內(nèi)外設(shè)備及內(nèi)裝飾件的質(zhì)量,將它們布置在安裝位置.模型5的質(zhì)量與車體實際設(shè)計質(zhì)量僅差0.561 t,并將這部分質(zhì)量以調(diào)整材料密度的形式均布在整車結(jié)構(gòu)上.
圖4 車體結(jié)構(gòu)自由模態(tài)振型云圖
五種整備狀態(tài)車體模態(tài)分析模型的集中質(zhì)量均是施加到重心上,以保證分析模型的重心與實際車體的重心更為接近.實際車體重心縱向坐標(biāo)為47.5 mm,垂向坐標(biāo)為1 837.9 mm,橫向坐標(biāo)為-13.9 mm;車體坐標(biāo)原點的橫向坐標(biāo)和縱向坐標(biāo)位于車體對稱中心,垂向坐標(biāo)原點位于軌面.表1給出了五種模型的重心與實際重心的對比結(jié)果,其中:△定義為與實際重心差值的絕對值.
表1 整備狀態(tài)車體重心位置
五種整備狀態(tài)車體模態(tài)分析模型的計算結(jié)果如圖5~圖9,圖中給出了車體一階垂向彎曲和一階扭轉(zhuǎn)振型云圖.整備狀態(tài)車體一階垂向彎曲和一階扭轉(zhuǎn)的振動頻率試驗值分別為10.98與11. 37[7].五種分析模型的分析結(jié)果與實驗值的誤差分析如表2所示,表中頻率單位為Hz,誤差百分?jǐn)?shù)定義為︱(計算值-試驗值) /試驗值︱%.
圖5 模型1的振型云圖
圖6 模型2的振型云圖
圖7 模型3的振型云圖
圖8 模型4的振型云圖
圖9 模型5的振型云圖
表2 整備狀態(tài)下車體模態(tài)分析結(jié)果
結(jié)合表1和表2可得:模型1中質(zhì)量分布遠(yuǎn)離實際,重心位置與實際值誤差大,尤其是垂向坐標(biāo)的△為445 mm,所以,模型1的計算頻率值與實際試驗值誤差高達(dá)40.73%,不能采用模型1計算整備狀態(tài)車體模態(tài);模型2的車體一階垂向彎曲振動頻率的計算值為10.258 Hz,接近實際試驗值10.258 Hz,誤差為6.58%.扭轉(zhuǎn)振動頻率的誤差為3.80%.建議在車體方案設(shè)計階段,整備狀態(tài)車體模態(tài)分析時可以采用模型2.
模型3和模型4的區(qū)別僅是5.927 t質(zhì)量的分布方式不同.與模型3相比,模型4的重心位置更接近實際重心,其計算頻率與試驗值誤差相對較小.所以,分析模型中的質(zhì)量分布應(yīng)盡可能細(xì)分質(zhì)量,若不具備條件時,應(yīng)采用調(diào)整材料密度的方法分布質(zhì)量,不能像靜強度分析模型中將質(zhì)量均布在底架上.
模型5的重心位置幾乎與實際值一致,已將內(nèi)裝飾中可能對結(jié)構(gòu)剛度有貢獻(xiàn)的結(jié)構(gòu)納入了計算模型中.車體一階垂向彎曲振動頻〗率與試驗值的誤差僅為0.03%,扭轉(zhuǎn)振動頻率的誤差為6. 85%.在詳細(xì)設(shè)計階段,采用模型5分析整備狀態(tài)車體模態(tài).
高速動車組車體整體剛度和質(zhì)量分布決定著整備狀態(tài)車體的整體低階振動品質(zhì).車體剛結(jié)構(gòu)模態(tài)分析結(jié)果的準(zhǔn)確性能夠反映模型模擬車體剛度的準(zhǔn)確性;整備狀態(tài)車體模態(tài)分析結(jié)果的準(zhǔn)確性能夠反映模型模擬車體質(zhì)量分布的準(zhǔn)確性.
( 1)車體剛結(jié)構(gòu)一階垂向彎曲振動頻率的計算值與試驗值的誤差為2.75%;一階扭轉(zhuǎn)振動頻率的計算值與試驗值的誤差為7.9%.這意味著基于“剛度等效”原則建立的高速動車組車體剛結(jié)構(gòu)模態(tài)分析模型較為準(zhǔn)確地模擬了車體的整體垂向剛度以及扭轉(zhuǎn)剛度;
( 2)與實際車體重心最為接近的整備狀態(tài)車體仿真模型5的計算結(jié)果:一階垂向彎曲振動頻率的計算值與試驗值的誤差為0.03%;一階扭轉(zhuǎn)振動頻率的計算值與試驗值的誤差為6.85%.
( 3)通過改變材料密度的方法將非承載結(jié)構(gòu)質(zhì)量施加到整車結(jié)構(gòu)上的快速且較為精確的建模方法(模型2采用的方法)適用于車體方案設(shè)計階段估算整備狀態(tài)車體的低階整體振動模態(tài).
[1]李芾,安琪.國內(nèi)外高速動車組的發(fā)展[J].電力機車與城軌車輛,2007,30( 5) : 1-5.
[2]郝魯波,胡青泥,李剛.整備狀態(tài)下客車模態(tài)的有限元分析探討[J].鐵道車輛,2004,42( 11) : 4-7.
[3]曲天威,王惠玉.機車車體和轉(zhuǎn)向架模態(tài)分析研究[J].鐵道機車車輛,2012,32( 3) : 5-8.
[4]趙陽陽,劉士煜,張偉.車體下吊掛方式和彈性懸置質(zhì)量對車體模態(tài)頻率的影響[J].佳木斯大學(xué)學(xué)報(自然科學(xué)版),2012,30( 3) : 326-329.
[5]宮島,周勁松,孫文靜,等.高速列車彈性車體與轉(zhuǎn)向架耦合振動分析[J].交通運輸工程學(xué)報,2011,11 ( 4) : 41-47.
[6]宮島,周勁松,孫文靜,等.高速列車車下設(shè)備模態(tài)匹配研究[J].振動與沖擊,2014,33( 8) : 180-185.
[7]青島四方車輛研究所.某動車組車體模態(tài)試驗報告[R].青島四方車輛研究所,2015.
Modeling Method of High-Speed EMU Car Body Modal Analysis and Experimental Verification
XIE Suming1,XUE Ningxin1,MA Menglin2
( 1.School of Traffic and Transportation Engineering,Dalian Jiaotong University,Dalian 116028,China; 2.CNR Changchun Railway Vehicles Co.,Ltd,Changchun 130062,China)
Finite element model of a car body for independent developed high-speed EMU is set up based on stiffness equivalence principle.Errors of numerical and modal experimental results are 2.75% and 7.90% respectively,for first-order vertical bending frequency and torsional frequency.On this basis,five modal analysis models of equipped car body are established based on distributing mass principle in which gravity center of the model agrees with physical prototype.The frequencies from the fifth model are most closed to the experimental results,and the errors are 0.03% and 6.85%.The suggestion is that the second modeling method is feasible to simulate the equipped car body in design stage.
high-speed EMU; equipped car body; modeling method; modal analysis
A
1673-9590( 2016) 01-0011-04
2015-03-25
中國鐵路總公司科技研究開發(fā)計劃資助項目( 2014J004-N) ;國家科技支撐計劃資助項目( 2013BAG21Q01)
謝素明( 1965-),女,教授,博士,主要從事車輛工程CAE關(guān)鍵技術(shù)的研究
E-mail: smx@ djtu.edu.cn.