宋宇豪，孫麗萍，王玉艷

(大連交通大學機車車輛工程學院，遼寧 大連 116028)

目前，我國城市呈現(xiàn)多中心或者放射狀發(fā)展趨勢，交通擁堵問題成為制約城鎮(zhèn)化進一步發(fā)展的瓶頸。市域快軌列車以打造一小時內(nèi)都市圈中程出行為主旨，將成為跨區(qū)間客流運營的中流砥柱[1]。在保證靜態(tài)性能(如車輛剛度和強度)和動態(tài)特性(如固有頻率)的基礎(chǔ)上，優(yōu)化市域快軌車體結(jié)構(gòu)并減少質(zhì)量，是目前市域快軌車研究中最重要的課題之一[2]。同時，現(xiàn)有市域快軌車體內(nèi)的筋板布局方案多由設(shè)計人員憑經(jīng)驗設(shè)計，在方案的試驗與改進中耗費了大量成本，而且往往無法設(shè)計出經(jīng)濟和性能最優(yōu)的方案。優(yōu)化設(shè)計可以有效提高車輛性能、避免材料浪費、節(jié)約實驗時間和降低制造成本，為設(shè)計提供參考。因此，有必要對列車車體內(nèi)的筋板布局進行優(yōu)化設(shè)計。

本文以某市域快軌車為研究對象，對其車頂中空鋁合金型材斷面進行拓撲優(yōu)化。車頂采用鋁合金擠壓型材，以及車體結(jié)構(gòu)的剛度是本次優(yōu)化設(shè)計過程中需要考慮的重要因素。根據(jù)優(yōu)化計算結(jié)果設(shè)計出新的車頂筋板布局方案并驗證其可行性，實現(xiàn)車頂筋板合理布局，達到輕量化設(shè)計的目的。

1 車體結(jié)構(gòu)簡介

該市域快軌車輛車體結(jié)構(gòu)采用全長的大型中空鋁合金擠壓型材組焊成筒型整體承載結(jié)構(gòu)，車體由底架、側(cè)墻、車頂、端墻等焊接而成。車輛最高運行速度為160 km/h，車體長度為21880 mm，車體寬度為3000 mm，車體高度(距軌面)為3500 mm，車輛定距為15700 mm。

2 車體有限元模型

本文利用HyperMesh有限元分析軟件建立車體有限元模型，采用20 mm左右的殼單元離散車體結(jié)構(gòu)。本車主要采用的是四邊形單元，同時采用三角形單元進行局部過渡，沖擊座、抬車座和牽引變壓器安裝座采用實體單元。螺栓采用RBE2剛性單元和CBEAM單元。車體鋁合金的彈性模量為69 GPa，泊松比為0.3。車體有限元模型包括272.3萬個單元和201.5萬個節(jié)點。原車頂由5段3種擠壓型材焊接而成，長21596 mm，寬2175.7 mm，質(zhì)量為926.75 kg，車頂及其橫斷面筋板布局如圖1所示。

2.1 計算工況

根據(jù)《BS EN12663-1：2010+A1：2014鐵路應用-鐵路車輛車體的結(jié)構(gòu)要求》和該市域快軌車輛車體設(shè)計任務書，選取了四種比較危險的計算工況如下：

(1)超載工況：車體自重+超員重量。

(2)最大垂載工況：1.3×(車體自重+超員重量)。

(3)超載壓縮工況：1200 kN縱向壓縮載荷+車體自重+超員重量。

(4)超載拉伸工況：960 kN縱向拉伸載荷+車體自重+超員重量。

車體整備質(zhì)量為30662 kg，其中空調(diào)質(zhì)量為760 kg、牽引變流器質(zhì)量為1776 kg、變壓器質(zhì)量為3721.9 kg。以重力加速度的形式在車體鋼結(jié)構(gòu)施加其重量，g取9.81 m/s2；超員重量均布在車體底架地板上表面；空調(diào)重量均布在車頂?shù)目照{(diào)安裝座上，其他各吊掛設(shè)備的重量均在對應位置的質(zhì)心處以集中力的形式施加。

四種工況的垂向約束均施加在該市域快軌車底架枕梁的四個空氣彈簧處，選取中心銷處施加橫向約束，在二位端沖擊座處施加縱向約束。

2.2 評價標準

根據(jù)規(guī)范和設(shè)計任務書的要求，確定該車的評價標準如下：

剛度評價標準：在超載工況作用下，車體靜撓度不超過兩轉(zhuǎn)向架支撐點間距離(車輛定距)的1‰。該市域快軌車的車輛定距為15700 mm，則在超載工況作用下，車體靜撓度應不超過15.7 mm。

靜強度評價標準：在各計算載荷工況下車體結(jié)構(gòu)的Von Mises應力應不大于許用應力，許用應力取材料的屈服強度。

3 拓撲優(yōu)化

3.1 拓撲優(yōu)化簡介

優(yōu)化設(shè)計有三要素，即設(shè)計變量、目標函數(shù)和約束條件。設(shè)計變量是發(fā)生改變從而提高性能的一組參數(shù)；目標函數(shù)要求最優(yōu)的設(shè)計性能；約束條件是對設(shè)計的限制。拓撲優(yōu)化是一種根據(jù)給定的負載情況、約束條件和性能指標，在給定的區(qū)域內(nèi)對材料分布進行優(yōu)化的數(shù)學方法，是結(jié)構(gòu)優(yōu)化的一種[3]。拓撲優(yōu)化的研究領(lǐng)域主要分為連續(xù)體拓撲優(yōu)化和離散結(jié)構(gòu)拓撲優(yōu)化，兩個領(lǐng)域都需要依賴于有限元方法。本文采用成熟的連續(xù)體拓撲優(yōu)化進行研究，其方法是將優(yōu)化空間的材料離散成有限個單元，通過給定算法進行計算，確定可以去除的部分，綜合考慮制造工藝以及運行條件等方面的因素，從而設(shè)計出更加合理的材料布局。拓撲優(yōu)化中常用的拓撲表達形式和插值模式有：均勻化方法、密度法(即方法SIMP)、變厚度法和拓撲函數(shù)描述方法[4]。

本文應用的OptiStruct拓撲優(yōu)化的材料模式采用密度法(SIMP方法)，即將有限元模型設(shè)計空間的每個單元的“單元密度”作為設(shè)計變量。單元密度是在0～1之間連續(xù)取值，而取值的大小同該結(jié)構(gòu)的彈性模量E之間存在某種函數(shù)關(guān)系。優(yōu)化求解后單元密度等于或趨近1，表明該單元位置處的材料非常重要，需要保留；如果單元材料密度等于或趨近0，則表明該單元處的材料不重要，可以去除，從而達到材料的高效率利用，實現(xiàn)輕量化設(shè)計[3]。SIMP方法用公式表達為：

E(xi)=Emin+(xi)P(E0-Emin)

(1)

式中，插值之后的彈性模量為E(xi)，實體部分的彈性模量為E0，孔洞部分材料的彈性模量為Emin，單元相對密度為xi(取值為1表示有材料，取值為0表示無材料即孔洞)，懲罰因子為P[5]。

3.2 車頂斷面拓撲優(yōu)化

刪除有限元模型中車頂?shù)慕畎?，并用實體單元對車頂空間進行填充，將該填充部分設(shè)定為拓撲優(yōu)化的設(shè)計區(qū)域。由于車頂為擠壓型材，為保證型材橫截面在縱向上保持一致，對設(shè)計區(qū)域施加沿縱向的擠壓約束。車體結(jié)構(gòu)的剛度是設(shè)計中需要考慮的重要因素，在設(shè)計過程中需要考慮多載荷工況，應變能越小剛度越大，而各工況下車體結(jié)構(gòu)的應變能均不相同。為保證在各工況下車頂結(jié)構(gòu)的剛度盡可能大，本文將剛度最大化的拓撲優(yōu)化問題轉(zhuǎn)為要求各工況加權(quán)應變能最小，可以避免多工況拓撲優(yōu)化問題轉(zhuǎn)化成多目標拓撲優(yōu)化問題，否則求解將變得十分困難[6]?？紤]到車頂筋板體積占車頂部分區(qū)域的比例很小，故體積分數(shù)上限不宜設(shè)置過大，否則優(yōu)化結(jié)果將保留大量冗余材料；同時，由于車頂應力較小，為防止在優(yōu)化過程中設(shè)計區(qū)域的材料被全部去除，需要對體積分數(shù)設(shè)定下限。這樣基于變密度法，在優(yōu)化過程中建立了兩個響應，即體積分數(shù)和加權(quán)應變能。拓撲優(yōu)化的參數(shù)如下：

目標函數(shù)：選定工況下的加權(quán)應變能最小；

設(shè)計變量：設(shè)計區(qū)域的密度；

約束條件：設(shè)計區(qū)域的體積分數(shù)大于10%，小于30%。

為使優(yōu)化更加充分，設(shè)置收斂容差為10-5；為防止優(yōu)化由于步數(shù)限制提前終止，設(shè)置優(yōu)化步數(shù)上限為100步；由于設(shè)計區(qū)域為體單元，為使材料去除更加充分，設(shè)置離散度參數(shù)(Discrete)為3。

在最大垂載工況、超載壓縮工況和超載拉伸工況下進行拓撲優(yōu)化，經(jīng)88步迭代后計算收斂，得到迭代曲線如圖2所示，拓撲優(yōu)化計算結(jié)果如圖3所示。

從圖2迭代曲線可以看出，最終拓撲優(yōu)化結(jié)果的加權(quán)應變能相比于最初的加權(quán)應變能變小，說明車頂結(jié)構(gòu)的剛度有所提升。上述利用Optistruct得出的優(yōu)化計算結(jié)果可以作為車頂筋板布局設(shè)計的基礎(chǔ)，車頂兩側(cè)筋板布局根據(jù)圖3中的優(yōu)化結(jié)果確定；密度云圖3顯示車頂中間部分無材料，為加強車頂剛度，結(jié)合車頂受力情況，需要在中間部分添加筋板。最后設(shè)計得到的車頂斷面筋板布局如圖4所示。優(yōu)化設(shè)計后的車頂結(jié)構(gòu)質(zhì)量為839.45 kg，較原車頂結(jié)構(gòu)減重87.3 kg，減重率為9.42%。

4 優(yōu)化結(jié)果校核

根據(jù)優(yōu)化設(shè)計得到的筋板布局建立新的有限元模型，在各工況下對模型進行計算分析，根據(jù)計算結(jié)果數(shù)據(jù)檢驗優(yōu)化設(shè)計后的車體結(jié)構(gòu)力學性能是否滿足標準要求。

剛度校核：超載工況下車體底架邊梁中部的垂向位移為12 mm，小于剛度評價標準所得的15.7 mm，滿足剛度要求。

靜強度校核：

(1)最大垂載工況下的最大VonMises應力為168.2 MPa，出現(xiàn)在側(cè)墻門下角處，小于材料的許用應力215 MPa，應力云圖如圖5所示。

(2)超載壓縮工況下的最大VonMises應力為154.2 MPa，出現(xiàn)在一位端牽引梁加強板處，小于材料的許用應力240 MPa，應力云圖如圖6所示。

(3)超載拉伸工況下的最大VonMises應力為189.5 MPa，出現(xiàn)在牽引變壓器安裝座處，小于材料的許用應力415 MPa，應力云圖如圖7所示。

各計算工況下車體結(jié)構(gòu)VonMises應力計算結(jié)果及其發(fā)生部位如表1所示。

表1 各計算工況下車體VonMises應力計算結(jié)果

由計算數(shù)據(jù)可知，優(yōu)化后各結(jié)構(gòu)的應力均小于對應材料的許用應力，滿足靜強度要求。

5 結(jié)論

本文利用Optistruct求解器對某市域快軌車車頂中空鋁合金型材斷面進行拓撲優(yōu)化，根據(jù)優(yōu)化結(jié)果顯示的車頂橫斷面的材料密度云圖，結(jié)合車頂受力情況，設(shè)計出新的車頂斷面筋板布局。新方案車頂筋板布局主要采用直筋，車頂結(jié)構(gòu)剛度較原車頂結(jié)構(gòu)有所提高。在根據(jù)規(guī)范和設(shè)計任務書的要求選取出的工況下，對新方案車體結(jié)構(gòu)的位移和VonMises應力進行計算分析，結(jié)果顯示優(yōu)化方案的車體剛度和靜強度均滿足評價標準要求。新方案的車頂結(jié)構(gòu)較原車頂結(jié)構(gòu)采用了更少的材料，實現(xiàn)減重87.3 kg，減重率為9.42%，輕量化效果明顯，同時降低了制造成本，為車頂中空鋁合金型材斷面布局設(shè)計提供參考。