歐陽(yáng)柏添 蔡家富 劉桂雄

摘 ?要：針對(duì)裝配車體撓度調(diào)整過程中存在過程麻煩、容易出錯(cuò)、人力耗費(fèi)大的問題，文章從裝配車體模型入手，研究基于線性規(guī)劃的軌道交通車輛裝配車體撓度調(diào)整方法。通過建立裝配車體有限元模型，求解裝配車體調(diào)整過程各調(diào)整點(diǎn)調(diào)整力與裝配車體撓度變形關(guān)系。通過三次樣條曲線擬合方法對(duì)裝配車體撓度變形曲線進(jìn)行擬合，并利用線性規(guī)劃方法求解最優(yōu)裝配車體撓度調(diào)整方法。結(jié)果表明，利用該方法可快速計(jì)算出各調(diào)整點(diǎn)撓度調(diào)整力最小的最優(yōu)調(diào)整方法，指導(dǎo)技術(shù)人員對(duì)裝配車體進(jìn)行調(diào)整。

關(guān)鍵詞：裝配車體;撓度調(diào)整;三次樣條曲線擬合;線性規(guī)劃

中圖分類號(hào)：TP29 ?文獻(xiàn)標(biāo)志碼：A文章編號(hào)：2096-4706（2021）13-0144-04

Research on Deflection Adjustment Method of Rail Transit Vehicle Assembly Carbody Based on linear Programming

OUYANG Baitian1， CAI Jiafu2， LIU Guixiong2

（1.Guangdong Jiangmen Supervision Testing Institute of Quality & Metrology， Jiangmen ?529030， China;

2. College of Mechanical and Automotive Engineering， South China University of Technology， Guangzhou ?510640， China）

Abstract： In view of the problem of troublesome process， error-proneness， and high labor cost in the assembly car body deflection adjustment process， this paper starts with the assembly car body model and studied the rail transit vehicle assembly car body deflection adjustment method based on the linear programming method. By establishing a finite element model of the assembly car body， the relationship between the adjustment force of each adjustment point and the deflection of the assembly car body during the adjustment process of the assembly car body is solved. The cubic spline curve fitting method is used to fit the deflection deformation curve of the assembly car body， and the linear programming method is used to solve the optimal assembly car body deflection adjustment method. The results show that this method could quickly calculate the optimal adjustment method with the smallest deflection adjustment force at each adjustment point， and guide the technicians to adjust the assembled car body.

Keyword： assembly car body; deflection adjustment; cubic spline curve fitting; linear programming

0 ?引 ?言

軌道交通裝配車體是軌道交通的重要一環(huán)，軌道交通裝配車體可用于運(yùn)輸乘客往返城市與郊區(qū)。裝配車體底架撓度是裝配車體的一項(xiàng)重要結(jié)構(gòu)參數(shù)，其中撓度調(diào)整主要是在裝配車體總成過程中進(jìn)行。目前對(duì)裝配車體底架撓度調(diào)整的方法主要是使用拉頂裝置或者采用火焰調(diào)修實(shí)現(xiàn)裝配車體撓度調(diào)整?；鹧嬲{(diào)修雖然能較好控制裝配車體撓度，但是容易對(duì)裝配車體底架性能產(chǎn)生影響[1，2]。在使用拉頂裝置對(duì)裝配車體撓度進(jìn)行調(diào)整研究中，文獻(xiàn)[3]（2017）通過對(duì)各個(gè)調(diào)整位置上機(jī)械支撐裝置設(shè)置不同高度，實(shí)現(xiàn)撓度對(duì)應(yīng)位置反變形;文獻(xiàn)[4]（2009）通過采用增加液壓拉緊裝置數(shù)量及提高裝置壓力的方法，克服裝配車體制造過程中撓度偏低的問題;文獻(xiàn)[5]在考慮車體底架撓度回彈基礎(chǔ)上，利用F形鋼卡和千斤頂上頂車體中部位置及下拉車體端部位置，調(diào)整裝配車體撓度值到達(dá)設(shè)計(jì)值。國(guó)內(nèi)裝配車體撓度調(diào)整方法多基于實(shí)踐經(jīng)驗(yàn)總結(jié)，缺乏數(shù)學(xué)模型支持，尚未有較好的裝配車體撓度調(diào)整方法。目前國(guó)內(nèi)外許多研究學(xué)者對(duì)車體撓度變形進(jìn)行模擬，文獻(xiàn)[6]（2018）利用有限元方法對(duì)不銹鋼點(diǎn)焊車體進(jìn)行靜強(qiáng)度分析，通過分析車體在不同整備狀態(tài)及不同工況下車體撓度變形及垂向彎曲振型頻率，驗(yàn)證車體質(zhì)量合格狀況。文獻(xiàn)[7]（2017）以殼單元為主體建立客車車身有限元模型，計(jì)算車身在緊急制動(dòng)工況下車身、車架、地板骨架橫向及縱向撓度變形，分析車身?yè)隙茸冃屋^大的危險(xiǎn)位置，用于車身優(yōu)化設(shè)計(jì)。文獻(xiàn)[8]研究鋼桁架平臺(tái)車架體添加勁板結(jié)構(gòu)性能影響，利用有限元方法分析三種不利荷載工況下車架體大梁端部撓度變形并進(jìn)行撓度曲線擬合，研究結(jié)果表明添加勁板可有效減小車架體撓度變形。文獻(xiàn)[9]（2020）提出一種應(yīng)變感知和有限元分析車體敞口梁撓度計(jì)算方法，通過有限元分析方法構(gòu)建車體底部撓度與應(yīng)變間數(shù)學(xué)模型，對(duì)車體敞口梁底部撓度進(jìn)行求解，并采用多項(xiàng)式曲線擬合求解車體撓度曲線方程。本文利用有限元分析思想，研究裝配車體調(diào)整力與裝配車體撓度變形關(guān)系，建立撓度調(diào)整約束條件與撓度調(diào)整目標(biāo)函數(shù)，對(duì)裝配車體調(diào)整過程中各撓度控制點(diǎn)調(diào)整力進(jìn)行求解，指導(dǎo)技術(shù)人員對(duì)裝配車體底架進(jìn)行調(diào)整，提高裝配車體撓度調(diào)整效率，縮短裝配車體總成制造時(shí)間。

1 ?軌道交通車輛裝配車體撓度調(diào)整問題描述及整體方法研究

1.1 ?問題描述

裝配車體底架調(diào)整位置如圖1所示，調(diào)整過程控制點(diǎn)在車體兩側(cè)端部、枕中區(qū)域、過渡區(qū)域、車體中部共14個(gè)位置。

目前裝配車體底架撓度調(diào)整方法為：（1）以4處枕中位置撓度值為基準(zhǔn)點(diǎn)進(jìn)行調(diào)整，使4處枕中位置撓度相等，對(duì)裝配車體底架4處枕中位置進(jìn)行固定。（2）依次對(duì)左側(cè)及右側(cè)車中、端部、過渡區(qū)域位置進(jìn)行撓度調(diào)整，手動(dòng)調(diào)整對(duì)應(yīng)調(diào)整點(diǎn)，使裝配車體底架邊梁撓度達(dá)到要求范圍，使車中、端部、過渡區(qū)域位置撓度ωi∈[ωmin[i]， ωmax[i]]（i∈[1，10]），其中ωmin[i]表示第i個(gè)調(diào)整位置撓度調(diào)整合格最小值，ωmax[i]表示第i個(gè)調(diào)整位置撓度調(diào)整合格最大值。（3）調(diào)整完成后對(duì)各位置撓度進(jìn)行檢查，保證底架各測(cè)量點(diǎn)撓度在合格范圍內(nèi)。

目前軌道交通車輛裝配車體撓度調(diào)整方法尚不完善，主要存在的問題有：（1）過程麻煩。車體底架4處枕中位置撓度固定后，需對(duì)車中、端部、過渡區(qū)域共10個(gè)位置分別進(jìn)行撓度調(diào)整，每次調(diào)整需要將調(diào)整位置外撓度調(diào)整點(diǎn)進(jìn)行固定，調(diào)整效率低下。（2）容易出錯(cuò)。撓度調(diào)整過程中，需保證千斤頂、手拉葫蘆裝置已起到固定約束作用，且調(diào)整量難以控制在合適范圍，調(diào)整容易反復(fù)。（3）人力耗費(fèi)大。撓度調(diào)整過程中，需要2名技術(shù)人員，其中1名控制千斤頂、手拉葫蘆，另1名進(jìn)行撓度值監(jiān)測(cè)。

1.2 ?整體調(diào)整方法框架設(shè)計(jì)

本文針對(duì)裝配車體撓度調(diào)整過程存在問題，提出基于線性規(guī)劃的軌道交通車輛裝配車體撓度調(diào)整方法如圖2所示。方法步驟如下：（1）采集車體調(diào)整信息，根據(jù)車體類型建立車體有限元模型，進(jìn)行裝配車體調(diào)整模擬仿真，求解對(duì)應(yīng)調(diào)整點(diǎn)作用力下車體撓度的變化情況。（2）通過建立車體作用力與撓度變形關(guān)系方程，根據(jù)撓度調(diào)整合格情況建立車體撓度調(diào)整方程，基于MATLAB實(shí)現(xiàn)車體各點(diǎn)調(diào)整力的計(jì)算。（3）對(duì)求解得到車體各點(diǎn)撓度數(shù)值符合撓度標(biāo)準(zhǔn)時(shí)的各點(diǎn)撓度調(diào)整力及各點(diǎn)預(yù)期撓度調(diào)整值進(jìn)行顯示，并繪制車體撓度調(diào)整曲線，指導(dǎo)技術(shù)人員對(duì)裝配車體底架進(jìn)行撓度調(diào)整。技術(shù)人員只需通過計(jì)算出來各點(diǎn)撓度調(diào)整力對(duì)裝配車體底架撓度進(jìn)行調(diào)整，即可快速實(shí)現(xiàn)裝配車體撓度合格調(diào)整。

本文軌道交通車輛裝配車體撓度調(diào)整方法以裝配車體車身為研究對(duì)象，該方法也可應(yīng)用于裝配車體車頭的撓度調(diào)整中。

2 ?基于線性規(guī)劃的軌道交通車輛裝配車體撓度調(diào)整方法設(shè)計(jì)

2.1 ?裝配車體有限元模型建立

裝配車體結(jié)構(gòu)由車頂、側(cè)墻、端墻、底架組成，其中底架由地板、底架邊梁、枕梁、緩沖梁、牽引梁五個(gè)部分組成，采用SolidWorks對(duì)車體整體結(jié)構(gòu)進(jìn)行建模，模型側(cè)壁采用中空擠壓鋁型材結(jié)構(gòu)，通過將裝配車體各部分進(jìn)行裝配固定，得到裝配車體模型如圖3所示。

本文車體模型有限元分析軟件為Hypermesh，可實(shí)現(xiàn)對(duì)車體模型快速可靠有限元分析[10]。裝配車體是由大型中空擠壓鋁型材焊接而成，且各組件大部分由薄壁結(jié)構(gòu)所構(gòu)成，因此使用2D板殼單元來進(jìn)行模擬。裝配車體有限元模型如圖4所示。

對(duì)裝配車體枕中支撐點(diǎn)處施加垂向位移約束，對(duì)裝配車體對(duì)角線撐桿拉頂裝置施加縱向位移約束，對(duì)裝配車體側(cè)墻門口、車窗處施加橫向位移約束，通過在車體調(diào)整位置施加壓力模擬機(jī)械支撐裝置及手拉葫蘆裝置調(diào)整過程，并對(duì)模型進(jìn)行求解，圖5為裝配車體右側(cè)中部施加1 MPa作用力下裝配車體底架邊梁垂向撓度變形圖。

2.2 ?基于三次樣條曲線擬合的裝配車體撓度變形曲線求解方法

基于有限元分析法得到的軌道交通車輛裝配車體邊梁撓度變形數(shù)據(jù)較多，利用傳統(tǒng)多項(xiàng)式擬合得到撓度變形曲線各點(diǎn)誤差較大，曲線擬合效果不好。本文采用三次樣條曲線擬合方法[11]，對(duì)裝配車體撓度變形曲線進(jìn)行求解，建立車體作用力與撓度變形關(guān)系方程。

設(shè)車體撓度變形在點(diǎn)集[x1，x2，…，xn+1]上的撓度值為[f（x1），f（x2），…，f（xn+1）]，設(shè)函數(shù)（x）g（x）為裝配車體撓度三次樣條曲線擬合函數(shù)：

其中，g（x）在區(qū)間[xi，xi+1]（i∈[1，n]）內(nèi)是三次多項(xiàng)式。由于有限元模型車體邊梁各網(wǎng)格劃分均勻，因此邊梁上各節(jié)點(diǎn)間距離相等。設(shè)撓度變形曲線各點(diǎn)間距為h=xi+1-xi，函數(shù)在xi上二階導(dǎo)數(shù)為mi=gi”（xi），根據(jù)三次樣條曲線原理，其在區(qū)間[xi，xi+1]內(nèi)需滿足：

式中，Ai、Bi為常數(shù)。聯(lián)立式（1）和式（4）可得：

聯(lián)立式（3）及式（5）可得：

由gi-1’（xi）=gi’（xi）可得：

（mi-1+4mi+mi+1）h2=6[f（xi+1）-2f（xi）+f（xi-1）]

式中，i∈[2，n]。結(jié)合邊界條件m1=0和mn+1=0可求得裝配車體撓度變形曲線方程。

設(shè)裝配車體撓度調(diào)整過程中共需對(duì)u個(gè)點(diǎn)進(jìn)行調(diào)整，設(shè)求得的撓度變形曲線方程為{g（1）（x），g（2）（x），…，g（u）（x）}。

2.3 ?基于線性規(guī)劃的裝配車體撓度調(diào)整力計(jì)算方法

建立裝配車體撓度變形曲線后，可得到各作用力對(duì)撓度調(diào)整點(diǎn)撓度的影響大小，由撓度變形曲線方程有第j（j∈[1，u]）個(gè)調(diào)整作用力Fj對(duì)調(diào)整點(diǎn)xadjk（k∈[1，u]）的撓度變形影響為Fjg（j）（xadjk）。設(shè)撓度調(diào)整至合格值的最小調(diào)整量為，撓度調(diào)整至合格值的最大調(diào)整量為。設(shè)撓度變形矩陣G為：

則撓度調(diào)整需要滿足約束條件：

約束條件式（10）中，F(xiàn)=[F1，F(xiàn)2，…，F(xiàn)u]，。

考慮到裝配車體調(diào)整點(diǎn)調(diào)整力過大會(huì)對(duì)車體整體造成不良影響，各調(diào)整點(diǎn)撓度調(diào)整力應(yīng)較小，取線性規(guī)劃目標(biāo)函數(shù)：

z=min（|F1|+|F2|+…+|Fu|）

目標(biāo)函數(shù)含有絕對(duì)值，需對(duì)目標(biāo)函數(shù)進(jìn)行轉(zhuǎn)化。設(shè) ，約束條件可轉(zhuǎn)化為：

式中U=[U1，U2，…，Uu]T，V=[V1，V2，…，Vu]。目標(biāo)函數(shù)可轉(zhuǎn)為：

3 ?試驗(yàn)與分析

對(duì)HyperView軟件中裝配車體邊梁位置處各節(jié)點(diǎn)編號(hào)進(jìn)行導(dǎo)出，一側(cè)邊梁包含491個(gè)橫向節(jié)點(diǎn)撓度數(shù)據(jù)，利用MATLAB對(duì)配車體撓度變形數(shù)據(jù)進(jìn)行三次樣條擬合，擬合結(jié)果如圖6所示?？梢钥闯隼萌螛訔l曲線擬合方法可以對(duì)裝配車體撓度變形數(shù)據(jù)擬合具有較好效果，設(shè)裝配車體撓度調(diào)整點(diǎn)數(shù)為u=10，可分別求得單位作用力下車體各點(diǎn)撓度變形數(shù)據(jù)g（j）（xadjk）如表1所示。

撓度調(diào)整力計(jì)算基于MATLAB實(shí)現(xiàn)，預(yù)期調(diào)整最小值矩陣和預(yù)期調(diào)整最大值矩陣為=[-3.1，-0.1.12.9. -0.1，-3.1.-0.1.12.9.-0.1.-3.1]，=[-2.9，0.1，13.1，0.1， -2.9，0.1，13.1，0.1，-2.9]，利用COM組件法集成MATLAB和LabVIEW，使裝配車體撓度調(diào)整計(jì)算結(jié)果在LabVIEW中顯示，裝配車體撓度調(diào)整計(jì)算結(jié)果如圖7所示。

使用上述裝配車體撓度調(diào)整方法，可快速計(jì)算出車體各點(diǎn)撓度數(shù)值符合撓度標(biāo)準(zhǔn)時(shí)的各點(diǎn)撓度調(diào)整力及各點(diǎn)預(yù)期撓度調(diào)整值，指導(dǎo)技術(shù)人員對(duì)裝配車體底架進(jìn)行撓度調(diào)整。

3 ?結(jié) ?論

本文通過對(duì)裝配車體撓度調(diào)整過程進(jìn)行問題分析與方法框架設(shè)計(jì)，提出基于線性規(guī)劃的軌道交通車輛裝配車體撓度調(diào)整方法，建立裝配車體有限元模型并對(duì)撓度調(diào)整過程進(jìn)行仿真模擬。

利用三次樣條曲線擬合方法對(duì)裝配車體撓度變形曲線進(jìn)行擬合，并基于線性規(guī)劃方法求解最優(yōu)裝配車體撓度調(diào)整方法。計(jì)算表明，本方法可快速計(jì)算出各調(diào)整點(diǎn)撓度調(diào)整力最小的最優(yōu)調(diào)整方法，指導(dǎo)技術(shù)人員對(duì)裝配車體進(jìn)行調(diào)整。

在實(shí)際裝配車體撓度調(diào)整過程中，除撓度調(diào)整力盡量小外還需考慮其他調(diào)整因素，如何在考慮其他調(diào)整因素的情況下找到撓度調(diào)整最佳方法，并應(yīng)用于撓度調(diào)整實(shí)際過程中，這是下一步研究主要內(nèi)容。

參考文獻(xiàn)：

[1] 樓強(qiáng)天，顧亞桃.機(jī)車車架撓度控制技術(shù)及應(yīng)用 [J].機(jī)車車輛工藝，2014（6）：4-6.

[2] 王立夫，王金金，劉東軍，等.B型鋁合金地鐵車輛車體制造工藝 [J].軌道交通裝備與技術(shù)，2013（1）：1-3.

[3] 黃軍才.淺析城軌鋁合金車體扭曲變形 [J].金屬加工（熱加工），2017（16）：35-39.

[4] 朱向東.上海軌道交通11號(hào)線車輛車體扭曲變形和撓度控制工藝改進(jìn) [J].電力機(jī)車與城軌車輛，2009，32（3）：32-33.

[5] 馮亮，楊佳奇，劉宇，等.鉸接型車輛車體撓度研究 [J].機(jī)車車輛工藝，2018（4）：1-3.

[6] 蒼松.某型號(hào)時(shí)速250 km的不銹鋼點(diǎn)焊車體鋼結(jié)構(gòu)強(qiáng)度分析 [J].鐵路計(jì)算機(jī)應(yīng)用，2018，27（6）：40-42.

[7] 蔣冬清，李三雁，代春香.基于有限元法的客車制動(dòng)工況整車情況分析 [J].中國(guó)測(cè)試，2017，43（11）：140-144.

[8] 張緒進(jìn)，張祎.重大件碼頭鋼桁架平臺(tái)車有限元分析 [J].武漢理工大學(xué)學(xué)報(bào)（交通科學(xué)與工程版），2012，36（1）：145-148.

[9] WU Y，LIANG X，ZHOU W，et al. Load Identification and Deflection Monitoring of Opening Beam on Well-Hole Freight Trains [J]. IEEE Access，2020（99）：1.

[10] 張航，鄔平波.基于ABAQUS的車體強(qiáng)度分析 [J].機(jī)車電傳動(dòng)，2017（4）：32-35.

[11] 王冰冰，李淮江.基于三次樣條插值的硅壓阻式壓力傳感器的溫度補(bǔ)償 [J].傳感技術(shù)學(xué)報(bào)，2015，28（7）：1003-1007.

作者簡(jiǎn)介：歐陽(yáng)柏添（1976—），男，漢族，廣東江門人，總工程師，高級(jí)工程師，本科，主要研究方向：計(jì)量測(cè)試。