有軌電車動(dòng)態(tài)限界計(jì)算軟件設(shè)計(jì)與實(shí)現(xiàn)

楊 宇 張濟(jì)民 耿 躍

1(同濟(jì)大學(xué)鐵道與城市軌道交通研究院 上海 201804) 2(上海匯煦交通科技有限公司 上海 201108)

有軌電車動(dòng)態(tài)限界計(jì)算軟件設(shè)計(jì)與實(shí)現(xiàn)

楊 宇1張濟(jì)民1耿 躍2

1(同濟(jì)大學(xué)鐵道與城市軌道交通研究院 上海 201804)2(上海匯煦交通科技有限公司 上海 201108)

針對(duì)目前沒有形成有軌電車結(jié)構(gòu)和模塊編組特點(diǎn)的限界計(jì)算方法的現(xiàn)狀，提出有軌電車動(dòng)態(tài)限界(水平限界)算法。該算法基于解析幾何方法，結(jié)合列車動(dòng)力學(xué)分別推導(dǎo)出有軌電車在直線和曲線運(yùn)行情況下以及有、無強(qiáng)側(cè)向風(fēng)情況下車體橫向偏移量，計(jì)算出車輛動(dòng)態(tài)限界。同時(shí)通過Matlab軟件開發(fā)平臺(tái)設(shè)計(jì)了有軌電車動(dòng)態(tài)限界計(jì)算軟件。該軟件具有導(dǎo)入計(jì)算參數(shù)與導(dǎo)出計(jì)算結(jié)果的功能，通過軟件計(jì)算，可以得到不同運(yùn)行情況下有軌電車車輛橫向限界值。最后給出計(jì)算實(shí)例來對(duì)算法進(jìn)行驗(yàn)證，實(shí)驗(yàn)結(jié)果表明,該算法滿足實(shí)踐工程需求。該軟件可供有軌電車車輛限界設(shè)計(jì)參考使用。

有軌電車 動(dòng)態(tài)限界 解析幾何 Matlab GUI

0 引 言

自1899年以來，有軌電車在我國(guó)的發(fā)展歷經(jīng)了引進(jìn)—衰落—再發(fā)展的階段[1]，目前，現(xiàn)代有軌電車以諸多優(yōu)勢(shì)取代了舊式有軌電車，在國(guó)內(nèi)得到較為廣泛的應(yīng)用。據(jù)《2014年中國(guó)城市軌道交通市場(chǎng)發(fā)展報(bào)告》預(yù)測(cè)，2015年中國(guó)現(xiàn)代有軌電車新增運(yùn)營(yíng)線路里程達(dá)133.43公里，建設(shè)規(guī)模大幅增長(zhǎng)。其中，北京、上海、廣州、成都等諸多省市積極規(guī)劃建設(shè)有軌電車系統(tǒng)[2]。但總體來看，現(xiàn)代有軌電車在我國(guó)尚處于初步發(fā)展階段，尚未建立起系統(tǒng)的技術(shù)標(biāo)準(zhǔn)，車輛設(shè)計(jì)理念和方法不能簡(jiǎn)單照搬地鐵的相關(guān)規(guī)范，在起步階段需要建立起結(jié)合我國(guó)國(guó)情的現(xiàn)代有軌電車系統(tǒng)的國(guó)家標(biāo)準(zhǔn)，用以指導(dǎo)有軌電車系統(tǒng)的發(fā)展。

截至2016年11月，有軌電車系統(tǒng)的相關(guān)規(guī)范《現(xiàn)代有軌電車交通工程技術(shù)標(biāo)準(zhǔn)》正在接受審查，但是技術(shù)標(biāo)準(zhǔn)中關(guān)于有軌電車車輛限界的相關(guān)規(guī)定只列舉了原則性的要求，未明確給出限界的計(jì)算方法。目前，用于計(jì)算有軌電車車輛限界的方法大多參考《城市軌道交通建設(shè)和運(yùn)營(yíng)規(guī)則》(BOStrab)和《地鐵限界標(biāo)準(zhǔn)》(CJJ 96—2003)，其設(shè)計(jì)理念基本與地鐵限界設(shè)計(jì)標(biāo)準(zhǔn)一致[3-5]。由于有軌電車的運(yùn)行環(huán)境與地鐵有很大差異，因此借用地鐵限界設(shè)計(jì)理念計(jì)算有軌電車限界的做法是否準(zhǔn)確適用，仍需實(shí)踐檢驗(yàn)。

本文提出一種新型的基于解析幾何方法的有軌電車車輛限界(橫向部分)計(jì)算算法，相對(duì)于參考地鐵限界計(jì)算的方法，涉及到的參數(shù)較少，更加簡(jiǎn)潔方便。由于不同型號(hào)的有軌電車各項(xiàng)參數(shù)不同，但計(jì)算過程基本相同，為避免重復(fù)計(jì)算過程，利用Matlab軟件在復(fù)雜數(shù)值分析、運(yùn)算方面的明顯優(yōu)勢(shì)，采用Matlab軟件自帶的圖形用戶界面GUI(Graphical User Interface)開發(fā)可視化界面。通過該軟件可以很方便地實(shí)現(xiàn)限界的計(jì)算過程，用戶只需輸入相應(yīng)車型的參數(shù)，選定相應(yīng)運(yùn)行情況，無需掌握代碼具體內(nèi)容，便可以快捷得到動(dòng)態(tài)限界計(jì)算結(jié)果[6]。

1 有軌電車動(dòng)態(tài)限界計(jì)算方法

當(dāng)有軌電車在線路上運(yùn)行時(shí)，計(jì)算車體橫向動(dòng)態(tài)限界需考慮到列車的橫擺和側(cè)滾運(yùn)動(dòng)，經(jīng)分析計(jì)算得到車體最大偏移點(diǎn)處的內(nèi)、外側(cè)限界，最終相加得到總限界值。有軌電車動(dòng)車和拖車橫擺示意圖如圖1所示，側(cè)滾示意圖如圖2所示。

圖1 車體橫擺示意圖

圖2 車體側(cè)滾示意圖

無論列車橫擺和側(cè)滾運(yùn)動(dòng)方向如何，外側(cè)限界總?cè)Q于動(dòng)車的B1點(diǎn)，無需考慮拖車與懸浮車；而計(jì)算內(nèi)側(cè)限界時(shí)，拖車的內(nèi)側(cè)限界需要與動(dòng)車的內(nèi)側(cè)限界比較，找出動(dòng)車和拖車內(nèi)側(cè)限界的較大值[5]，作為車體內(nèi)側(cè)限界。

考慮橫擺運(yùn)動(dòng)的影響，車體外側(cè)限界的偏移最大點(diǎn)位于動(dòng)車的端部截面的一點(diǎn)處，即圖1中B1點(diǎn)。圖1中A1、A2點(diǎn)分別是轉(zhuǎn)向架兩車軸中心線與線路中心線的交點(diǎn)，以轉(zhuǎn)向架中心點(diǎn)為坐標(biāo)原點(diǎn)建立直角坐標(biāo)系。

根據(jù)所給參數(shù)和三角形相似計(jì)算出A1、A2點(diǎn)的坐標(biāo)。當(dāng)車輛在直線運(yùn)行時(shí)，可根據(jù)兩點(diǎn)坐標(biāo)計(jì)算求得偏移角，根據(jù)相似三角形定理可分別得到B1點(diǎn)橫向偏移量以及B1點(diǎn)據(jù)車體中心線的距離，即直線運(yùn)行工況下橫擺作用的外側(cè)動(dòng)態(tài)限界。

當(dāng)車輛在曲線運(yùn)行時(shí)，需要找出曲線的圓心坐標(biāo)及B1點(diǎn)距圓心距離，后者減去曲線半徑(外擺)或曲線半徑減去后者(內(nèi)擺)，即為曲線運(yùn)行工況下橫擺作用的外側(cè)動(dòng)態(tài)限界。

當(dāng)通過大半徑曲線時(shí)，動(dòng)車有可能向曲線內(nèi)側(cè)擺動(dòng)，當(dāng)曲線圓點(diǎn)橫坐標(biāo)大于前部端墻到車體中心縱向距離時(shí)，就需要考慮動(dòng)車車體內(nèi)側(cè)輪廓點(diǎn)造成的限界。找出曲線的圓心坐標(biāo)及B2點(diǎn)距圓心距離，曲線半徑減去后者即為動(dòng)車內(nèi)側(cè)動(dòng)態(tài)限界，并可以確定車體向內(nèi)偏移的距離。利用鉸接點(diǎn)與圓心坐標(biāo)即可求出拖車內(nèi)側(cè)動(dòng)態(tài)限界，以及拖車向內(nèi)偏移的距離。進(jìn)行比較后，較大值作為橫擺作用的內(nèi)側(cè)動(dòng)態(tài)限界。

考慮側(cè)滾運(yùn)動(dòng)的影響，以線路水平方向?yàn)闄M軸，列車垂向中心線為縱軸，建立直角坐標(biāo)系，如圖2所示。在考慮橫擺影響車體產(chǎn)生偏移的情況下，假設(shè)車體側(cè)墻為曲面，只需計(jì)算得到車體側(cè)墻曲面的圓點(diǎn)坐標(biāo)，用曲面半徑減去(外側(cè)側(cè)滾)或加上(內(nèi)側(cè)側(cè)滾)圓點(diǎn)橫坐標(biāo)值即可最終得到最大偏移點(diǎn)處的內(nèi)外側(cè)動(dòng)態(tài)限界。

當(dāng)有軌電車以空車運(yùn)行時(shí)，計(jì)算動(dòng)態(tài)限界只需考慮橫擺運(yùn)動(dòng)的影響；以重車運(yùn)行時(shí)，需要考慮橫擺和側(cè)滾的共同影響。

在有軌電車實(shí)際運(yùn)行時(shí)，還可能會(huì)受到強(qiáng)側(cè)向風(fēng)的影響。此時(shí)，車體產(chǎn)生較大的橫向位移，導(dǎo)致橫向止擋與構(gòu)架接觸，限制車體更大的橫向位移。因此在考慮橫擺的影響時(shí)，A1、A2點(diǎn)的縱坐標(biāo)只需考慮一系懸掛橫向位移和車體與轉(zhuǎn)向架橫向止擋間隙兩項(xiàng)參數(shù)。其余步驟與上述計(jì)算方法相同。

2 限界計(jì)算軟件實(shí)現(xiàn)

2.1 軟件概述

有軌電車動(dòng)態(tài)限界計(jì)算軟件采用Matlab的GUI模塊進(jìn)行編寫。該軟件采用了本文提出的新型的基于解析幾何的動(dòng)態(tài)限界計(jì)算方法，主要功能就是計(jì)算有軌電車動(dòng)態(tài)限界范圍。該軟件很好地避免了人工計(jì)算的重復(fù)繁雜和帶來的誤差，節(jié)省了人力物力。目前軟件可針對(duì)不同車型的有軌電車進(jìn)行動(dòng)態(tài)限界的計(jì)算。圖3給出了有軌電車動(dòng)態(tài)限界計(jì)算軟件的計(jì)算流程。

圖3 軟件計(jì)算方法流程圖

2.2 軟件主界面

Matlab下的GUI圖形用戶開發(fā)環(huán)境GUIDE提供了一系列的開發(fā)工具，設(shè)計(jì)者只需根據(jù)功能需要自由調(diào)用、組合不同控件，布局完成后進(jìn)行各個(gè)控件的程序代碼編寫[7]。有軌電車動(dòng)態(tài)限界計(jì)算軟件中每個(gè)GUI界面設(shè)計(jì)均包含有后綴名為.fig的界面布局文件和.m的運(yùn)行函數(shù)文件。

軟件主界面如圖4所示。主界面主要包括2個(gè)Axes用于放置背景圖；“橫向參數(shù)”、“垂向參數(shù)”、“計(jì)算”3個(gè)按鈕；“曲線半徑”靜態(tài)文本框及其后的可編輯文本框；“強(qiáng)側(cè)向風(fēng)”的靜態(tài)文本框和單選按鈕。

圖4 軟件主界面

對(duì)應(yīng)于每個(gè)按鈕都可以在對(duì)應(yīng)的.m文件中找到回調(diào)函數(shù)，對(duì)其進(jìn)行代碼編寫以實(shí)現(xiàn)相應(yīng)的功能，此處為調(diào)用對(duì)應(yīng)的GUI界面。

線路曲線半徑在可編輯文本框中輸入，其對(duì)應(yīng)的回調(diào)函數(shù)為：

r=str2num(get(handles.edit14,′string′));

單選按鈕選定后將側(cè)向風(fēng)的有關(guān)數(shù)據(jù)賦值給wind并導(dǎo)入到計(jì)算過程中；其回調(diào)函數(shù)為：

wind=get(handles.wind,′value′);

在主界面的.m文件中包含了所有的計(jì)算過程。當(dāng)輸入的參數(shù)保存在對(duì)應(yīng)的mat數(shù)據(jù)文件中，“計(jì)算”按鈕可將各個(gè)數(shù)據(jù)分配至對(duì)應(yīng)的參數(shù)名稱中，通過文件后面的計(jì)算得到相應(yīng)的限界計(jì)算結(jié)果，儲(chǔ)存在Matlab數(shù)據(jù)庫(kù)中。另外定義一個(gè)result_data.m函數(shù)文件用于調(diào)用結(jié)果數(shù)據(jù)，并將數(shù)據(jù)導(dǎo)入進(jìn)另一個(gè)result.m文件，用于結(jié)果的顯示。其中，點(diǎn)擊“計(jì)算”按鈕即可調(diào)用result.m文件對(duì)應(yīng)的GUI界面。計(jì)算過程的流程圖如圖5所示。

圖5 計(jì)算過程流程圖

2.3 參數(shù)輸入與輸出

主界面中的2個(gè)參數(shù)輸入按鈕分別對(duì)應(yīng)2個(gè)橫向和垂向參數(shù)輸入的GUI界面，如圖6-圖7所示。每個(gè)參數(shù)輸入的GUI界面分別包括1個(gè)Axes用于放置各項(xiàng)參數(shù)的示意圖；若干個(gè)顯示參數(shù)名稱的靜態(tài)文本框以及對(duì)應(yīng)的可編輯文本框，用于輸入不同參數(shù)值；1個(gè)“確定”按鈕。

圖6 橫向參數(shù)輸入界面

圖7 垂向參數(shù)輸入界面

每一個(gè)靜態(tài)文本框都可以直接在GUI圖形設(shè)計(jì)編輯界面里操作完成；每一個(gè)可編輯文本框的回調(diào)函數(shù)格式為：

L=str2num(get(handles.edit1,′string′));

右下角的“確定”按鈕用于保存輸入的參數(shù)值并關(guān)閉此界面，返回主界面。

2.4 計(jì)算結(jié)果

主界面的計(jì)算按鈕對(duì)應(yīng)于另外一個(gè)顯示結(jié)果的GUI界面，如圖8所示。界面包括1個(gè)Axes用于放置限界結(jié)果示意圖；若干個(gè)靜態(tài)文本框用于結(jié)果參數(shù)名稱設(shè)置及對(duì)應(yīng)的可編輯文本框用于結(jié)果數(shù)據(jù)的顯示；“導(dǎo)出”和“返回”2個(gè)按鈕。

圖8 計(jì)算結(jié)果界面

結(jié)果輸出界面對(duì)應(yīng)result.m文件，只需調(diào)用存入了結(jié)果數(shù)據(jù)的result_data.m文件，將結(jié)果顯示到對(duì)應(yīng)的可編輯文本框中即可。

“導(dǎo)出”按鈕可以用于計(jì)算結(jié)果數(shù)據(jù)的導(dǎo)出并保存，使用者可以自行定義保存路徑，結(jié)果數(shù)據(jù)以.xls格式保存。“返回”按鈕用于返回主界面。

3 實(shí)例驗(yàn)證

本文采用某型號(hào)有軌電車為例，將其相關(guān)參數(shù)輸入軟件進(jìn)行計(jì)算，得出其動(dòng)態(tài)限界，驗(yàn)證軟件的可行性和準(zhǔn)確性。

某型號(hào)有軌電車的相關(guān)參數(shù)詳見表1。

表1 某型號(hào)有軌電車相關(guān)參數(shù)

根據(jù)表1中的參數(shù)計(jì)算該車型有軌電車的動(dòng)態(tài)限界。設(shè)定曲線半徑為150 m，有強(qiáng)側(cè)向風(fēng)的情況下，計(jì)算結(jié)果為車體外側(cè)動(dòng)態(tài)限界值為1 485.7 mm，內(nèi)側(cè)動(dòng)態(tài)限界值取拖車內(nèi)側(cè)限界，為1 443.9 mm，總限界值為2 929.6 mm。當(dāng)設(shè)定曲線半徑為150 m，無強(qiáng)側(cè)向風(fēng)情況下運(yùn)行，計(jì)算結(jié)果為車體外側(cè)動(dòng)態(tài)限界值為1 509.4 mm，內(nèi)側(cè)動(dòng)態(tài)限界值取拖車內(nèi)側(cè)限界，為1 459 mm，總限界值為2 968.4 mm。由計(jì)算結(jié)果可知，軟件的有軌電車動(dòng)態(tài)限界算法具有一定的可行性，能夠滿足實(shí)際工程的需求。

4 結(jié) 語

本文基于解析幾何的方法，結(jié)合有軌電車動(dòng)力學(xué)提出一種新型的計(jì)算有軌電車動(dòng)態(tài)限界的算法，具有一定的現(xiàn)實(shí)指導(dǎo)意義。在此基礎(chǔ)上，通過Matlab GUI模塊將算法整合出一個(gè)圖形用戶界面，設(shè)計(jì)成有軌電車動(dòng)態(tài)限界計(jì)算軟件，并通過實(shí)例計(jì)算驗(yàn)證了該軟件的有效性，可供有軌電車車輛限界設(shè)計(jì)參考使用。

[1] 李元坤,苗彥英.國(guó)外現(xiàn)代有軌電車建設(shè)發(fā)展的啟示[J].城市軌道交通研究,2013,16(6):29-32.

[2] 中規(guī)院交通院綜合所.2015年中國(guó)現(xiàn)代有軌電車現(xiàn)狀一覽[J].中國(guó)軌道交通雜志,2015-6-29.

[3] 葉芹祿.有軌電車車輛限界的計(jì)算方法[J].鐵道勘測(cè)與設(shè)計(jì),2014(6):1-3.

[4] 王嘉鑫.浮動(dòng)車體型有軌電車車輛限界(橫向部分)計(jì)算方法分析[J].地下工程與隧道,2015(1):44-46,50.

[5] 呂鳳梅,趙建秋,閆曉明,等.100%低地板現(xiàn)代城市有軌電車限界和小曲線通過能力分析[J].鐵道車輛,2013,51(9):5-7.

[6] 秦輝,席裕庚.基于Matlab GUI的預(yù)測(cè)控制仿真平臺(tái)設(shè)計(jì)[J].系統(tǒng)仿真學(xué)報(bào),2006,18(10):2778-2781.

[7] 王華喬.基于MATLAB GUI的自動(dòng)控制原理仿真軟件設(shè)計(jì)[J].計(jì)量與測(cè)試技術(shù),2011,38(2):12-14.

[8] 李連峰,樊金鳳,周菁.鐵路機(jī)車車輛限界規(guī)校準(zhǔn)方法研究[J].鐵道技術(shù)監(jiān)督,2016,44(11):21-23,26.

[9] 李磊,張茂帆,李芾,等.懸掛單軌交通限界計(jì)算方法研究[J].鐵道機(jī)車車輛,2016,36(1):119-126.

[10] 張愷,張玉鈞.基于Matlab-GUI的衍射圖像分析[J].計(jì)算機(jī)應(yīng)用與軟件,2011,28(3):243-245,281.

DESIGNANDIMPLEMENTATIONOFDYNAMICGAUGECALCULATIONSOFTWAREFORTRAMS

Yang Yu1Zhang Jimin1Geng Yue2

1(InstituteofRailTransit,TongjiUniversity,Shanghai201804,China)2(ShanghaiHuixuTransportationScienceandTechnologyCo.,Ltd.,Shanghai201108,China)

Because there is no calculation method of dynamic gauge for trams with its structure and module marshalling characteristics, this paper proposes a calculation method of dynamic gauge (horizontal gauge) for trams. Based on analytic geometry method and train dynamics, this algorithm derived the lateral offset of car-body when trams operated under different circumstances including in line and curve and with a strong side wind or without strong side wind, then the dynamic gauge of trams was calculated. At the same time, the dynamic gauge calculation software for trams was developed through the Matlab software platform. This software has functions of import calculation parameters and export calculation results. Through this software, the lateral gauge values of trams under different situations can be calculated. Experimental results show that this calculation method can meet the needs of practical engineering. Therefore this software can be used as references in vehicle gauge design for trams.

Trams Dynamic gauge Analytic geometry Matlab GUI

2017-02-22。楊宇，碩士生，主研領(lǐng)域：軌道車輛結(jié)構(gòu)及動(dòng)力學(xué)，列車空氣動(dòng)力學(xué)。張濟(jì)民，教授。耿躍，高工。

TP3

A

10.3969/j.issn.1000-386x.2017.12.026