CN115495831A - 地铁线路缓和曲线段的加宽曲线计算方法 - Google Patents

地铁线路缓和曲线段的加宽曲线计算方法 Download PDF

Info

Publication number
CN115495831A
CN115495831A CN202110673402.8A CN202110673402A CN115495831A CN 115495831 A CN115495831 A CN 115495831A CN 202110673402 A CN202110673402 A CN 202110673402A CN 115495831 A CN115495831 A CN 115495831A
Authority
CN
China
Prior art keywords
curve
point
train
widening
value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110673402.8A
Other languages
English (en)
Inventor
林凤涛
陈孟禹
钱鑫
陈武
方琴
江超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
East China Jiaotong University
Original Assignee
East China Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by East China Jiaotong University filed Critical East China Jiaotong University
Priority to CN202110673402.8A priority Critical patent/CN115495831A/zh
Publication of CN115495831A publication Critical patent/CN115495831A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Train Traffic Observation, Control, And Security (AREA)

Abstract

本发明涉及一种地铁线路缓和曲线段的加宽曲线计算方法,可以精确求解列车过缓和曲线时的加宽曲线。其具体计算步骤为:(1)对铁路缓和曲线线路进行函数拟合(2)基于数学函数的应用,联立曲率公式、正矢公式,求解列车过缓和曲线时的内外侧加宽点(3)确定列车在各个状态下加宽点,利用最小二乘法拟合各点求出列车加宽曲线。本发明在坚实的理论研究基础上开发了一种精确、系统、使用方便的缓和曲线地段限界加宽计算方法。本发明专利能用于任意工况下铁路限界的设计、建设。特别对于应对城市轨道交通小半径、多曲线的复杂工况能够保证良好的设计精度,针对列车过缓和曲线归纳出内侧和外侧的加宽计算方法能大大节省设计者时间和精力。

Description

地铁线路缓和曲线段的加宽曲线计算方法
技术领域
本发明涉及一种地铁线路缓和曲线段的加宽曲线计算方法,具体为一种求解地铁车辆过缓和曲线时内外侧加宽曲线的方法。
背景技术
随着我国城市规模和经济建设的快速发展、城市化进程逐步加快,地铁车辆运营速度的提高成为了城市轨道交通的发展趋势;与此同时,城市轨道交通小半径曲线也越来越多这对车辆的限界安全有直接的影响;面对变化越来越复杂的城市轨道线路和运营速度的提高。目前在建设施工中常用到的设计方法已经无法满足城市轨道交通对限界设计的需求,现有的缓和曲线地段限界加宽计算方法基本采用三种较为常规方法:一是先计算出圆曲线范围内的加宽值,然后将这个加宽值沿缓和曲线有规则的分配,这种计算方法的缺陷是精度不高会使得限界设计造成过度的冗余存在着工程浪费现象;二是采用图解法使用CAD软件绘制出缓和曲线的图像再模拟列车过曲线的过程;三是利用几何分析法针对线路缓和曲线段的过渡区段,提出一种基于线路中心线的精确加宽量计算方法,上述计算方法中有些含有大量三角函数,增加了算法的复杂性;亦或者没有建立车体函数,不便于动态实时的监测列车的运动变化且在一些特殊点上不适用;依靠CAD软件的方法实际运算也是非常复杂的,且对软件设备具有较高的要求。因此,本发明专利针对上述问题以建立数学函数为基础,使用微积分、最小二乘法拟合曲线以及导函数的应用通过建立车体函数,研究车体过缓和曲线段的运行过程;针对列车过缓和曲线归纳出内侧和外侧的加宽计算方法,最后对线路用三次幂函数拟合缓和曲线表达式进行限界加宽的分析。
发明内容
本发明的目的在于提供一种地铁缓和曲线段加宽的计算方法,以解决上述背景技术中提出的问题。
为实现上述目的,本发明提供如下技术方案:一种列车过缓和曲线段内外侧限界加宽的计算方法。
地铁线路缓和曲线段的加宽曲线计算方法,包括以下几个步骤:
S1,考虑车体因几何偏移产生的内侧限界加宽,包括以下几个小的步骤:
a1.任意选取一种地铁车辆车型,根据车型的转向架定距i以及车长I设定车体函数表达式;
b1.任意选取一种地铁缓和曲线的工况,根据《地铁设计规范》规定:线路平面圆曲线与直线之间应根据曲线半径、超高设置及设计速度等因素设置缓和曲线,常用的缓和曲线为三次方抛物线型起方程为:
Figure BDA0003119701010000011
式中,R为目标圆曲线的半径,单位m,L为缓和曲线的总长度,单位m,用于确定列车前后转向架的函数关系;
c1.建立一个直角坐标系,定义坐标系原点为直缓(ZH)点(直线和缓和曲线的切点,标注时用ZH表示);直缓(ZH)点到圆缓(YH)点(圆缓点是缓和曲线的终点和下一段圆曲线的起点相交的点)之间的曲线为缓和曲线;圆缓(YH)点之后为圆曲线;长方形ABCD表示地铁A型车车体;定义虚线Lab为车体函数;,Lab与曲线相交的动点b(xb,yb)、a(xa,ya)分别代表列车的前后转向架中心点;
d1.根据步骤a、b的求解车体函数的表达式,用直线Lab表示简化的车体,此时前后转向架都落在缓和曲线上,因此,可设
Figure BDA0003119701010000012
转向架定距lab=im,可知a、b是满足直线距离恒等于im的两个点,通过这一条件可以确定动点a、b之间的关系式:
Figure BDA0003119701010000013
在已知a点坐标(xa,ya)的情况下,可以通过计算a、b之间的关系式解出动点b的坐标,最后利用b的坐标,列出缓和曲线段列车的车体函数:
Figure BDA0003119701010000014
e1.求解车体最大内侧偏移点c的坐标值,建立动点a(xa,ya)与列车内侧最大偏移点c(xm,ym)之间的关系,已知lab=im,内侧最大偏移点c点同时也是列车中点,因此有:
Figure BDA0003119701010000015
根据距离公式,可列出最大偏移点c点与车辆后转向架中心点a的关系式:
Figure BDA0003119701010000021
在已知a点坐标的基础上就可以求解c(xm,ym)的坐标值;
f1.求解缓和曲线的曲率半径,现已知c(xm,ym)的坐标值,带入c点坐标对缓和曲线三次方抛物线求Kec值,可得:
Figure BDA0003119701010000022
式中,xm为c点的横坐标,Kec为缓和曲线对应c点的曲率,缓和曲线的曲率半径
Figure BDA0003119701010000023
因此,可建立曲率半径Rec随着c点的横坐标值xm的变化而变化的关系式;
以三次抛物线为函数分析其值等于:
Figure BDA0003119701010000024
式中,L是缓和曲线总长度,单位m,R是圆曲线半径,单位m,xm为c点的横坐标,Rec为缓和曲线的曲率半径,单位m;
g1.求解列车在缓和曲线上因几何偏移引起的内侧加宽值,使用正矢公式
Figure BDA0003119701010000025
可以对圆曲线上列车因几何偏移造成的加宽值做计算,但在缓和曲线段则需联立正矢公式和曲率公式计算缓和曲线上任意一点的正矢值fec
以拟合的缓和曲线函数为
Figure BDA0003119701010000026
为例,可到出处在于缓和曲线
Figure BDA0003119701010000027
上任意一点的正矢值求解方法:
Figure BDA0003119701010000028
式中,R为目标圆曲线的半径,单位m,L为缓和曲线长度,单位m,lab为车辆转向架定距,单位m,xm为c点的横坐标,fec为列车在缓和曲线上因几何偏移引起的内侧加宽值,单位mm;
S2,考虑缓和曲线线路超高产生的垂向内侧限界加宽值,包括以下小步骤:
a2.计算缓和曲线的超高Hec,缓和曲线的超高Hec为变量,Hec的值在0到Hmax的范围内均匀的变化,即Hec∈[0,Hmax];
Hec的计算方法:
Figure BDA0003119701010000029
其中:
Figure BDA00031197010100000210
其中:
Figure BDA00031197010100000211
式中,L为从直缓点算起至列车行驶到缓和曲线上任意一点的距离,单位m,L为缓和曲线的总长度,单位m,xb为列车前转向架中心点b的横坐标,xYH为圆缓点的横坐标;
由于积分为不可积积分,因此用幂级数展开成多项式,舍去高次项,再进行积分,可得:
Figure BDA00031197010100000212
Figure BDA00031197010100000213
前转向架中心点b点的超高最大因此需要带入前转向架中心b点的横坐标xb来计算,考虑xb、xYH,再带入原函数
Figure BDA00031197010100000214
做积分,就可以求出L、L、Hec的值;
b2.计算缓和曲线上超高引起的垂向内侧限界加宽值Wecv
Figure BDA00031197010100000215
式中,Wecv为列车在缓和曲线上超高引起的垂向内侧限界加宽值,单位mm,1500为钢轨中心距,单位mm,Hec为缓和曲线外轨超高值,单位mm,h为限界控制点自轨面的高度,单位mm;
S3,考虑缓和曲线上列车内侧限界总加宽,包括以下公式:
缓和曲线段内侧限界总加宽计算公式为:
Figure BDA0003119701010000031
式中,Wi为列车过缓和曲线时总的内侧加宽值,单位mm,fec为列车在缓和曲线上因几何偏移引起的内侧加宽值,单位mm,Wecv为列车在缓和曲线上超高引起的垂向内侧限界加宽值,单位mm,h为限界控制点自轨面的高度,单位mm;
S4,考虑缓和曲线外侧限界加宽,包括以下小步骤:
a4.对列车过缓和曲线求解外侧加宽值的方法,定义g(xg,yg)、e(xe,ye)、h(xh,yh)、f(xf,yf)四点;
b4.设点a(xa,ya)为列车后转向架中心点,b(xb,yb)为列车前转向架中心点,地铁A型车总车长lef=Im、转向架定距lab=im,可知
Figure BDA0003119701010000032
设动点e(xe,ye)、f(xf,yf);
c4.求列车车尾外侧曲线加宽:
Figure BDA0003119701010000033
可列长度方程
Figure BDA0003119701010000034
在已知a点坐标的情况下,很容易算出e点坐标(xe,ye),再做垂直与Lef的直线Leg
Figure BDA0003119701010000035
联立直线Leg和缓和曲线方程即可求出g点坐标(xg,yg),最后计算出来的
Figure BDA0003119701010000036
就是列车车尾外侧曲线加宽;
d4.求列车车头外侧曲线加宽:
Figure BDA0003119701010000037
可列长度方程
Figure BDA0003119701010000038
在已知b点坐标的情况下,很容易算出f点坐标(xf,yf),再做垂直与Lef的直线Lhf
Figure BDA0003119701010000039
联立Lhf和圆曲线方程即可求出h点坐标(xh,yh),最后计算出来的
Figure BDA00031197010100000310
就是列车车尾外侧曲线加宽;
其中:
Figure BDA00031197010100000311
其中:
Figure BDA00031197010100000312
式中,W车头为列车行驶在缓和曲线处车头的外侧曲线加宽值,单位mm,W车尾为列车行驶在缓和曲线处车尾的外侧曲线加宽值,单位mm;
S5,用最小二乘法拟合列车过缓和曲线时的加宽函数,包括以下步骤:
a5.确定列车后转向架xa的取值范围,设后转向架中心点为点a(xa,ya),那么(0≤xa<xYH);
b5.对xa在(0≤xa<xYH)平均取点,以这些点为基准,求解列车相对应的内侧最大偏移点c(xm,ym)的坐标;
c5.对xa在(0≤xa<xYH)平均取点,以这些点为基准,求出相对应的车尾最大偏移点e(xe,ye)、车头最大偏移点f(xf,yf)的坐标;
d5.用最小二乘法对S1、S2、S3、S4中所求的点进行函数拟合,最终求出在缓和曲线段车体函数内外侧加宽函数曲线。
优选的,考虑车体因几何偏移产生的内侧限界加宽的计算方法,S1的步骤a1中i为车辆转向架定距长度,单位m,I为车体车长长度,单位m,此算法适应与任意车型的限界加宽计算,S1的步骤b1中R代表缓和曲线后圆曲线的半径,单位m,L代表缓和曲线总长度,单位m,S1的步骤c1中动点b代表列车前转向架中心点,动点a代表列车后转向架中心点,Lab代表车体函数,此算法适应与任意铁路缓和曲线工况下的计算,其目的在于建立直角坐标系对S1的步骤a1、b1的缓和曲线进行拟合,建立车体函数便于监测车体的运动轨迹,S1的步骤d1其目的在于以列车后转向架中心点a为已知的基准点对车体函数进行分析,S1的步骤e1中lab代表转向架定距,c(xm,ym)代表列车内侧最大偏移点,S1步骤中的xm代表c点的横坐标,Kec代表缓和曲线对应c点的曲率,Rec代表缓和曲线对应c点的曲率半径,m-1,S1的步骤f1其目的在于找出a点与c点的函数关系,S1的步骤g1中fec表示列车在缓和曲线上因几何偏移引起的内侧加宽值,单位mm。
优选的,考虑缓和曲线线路超高产生的垂向内侧限界加宽值的计算方法,其中:S2的步骤a2中Hec代表缓和曲线外轨超高值,单位mm,L为从直缓点算起至列车行驶到缓和曲线上任意一点的距离,单位m,L为缓和曲线的总长度,单位m,xb为列车前转向架中心点b的横坐标,xYH为圆缓点的横坐标,其目的在于利用微积分等知识计算变量Hec,S2的步骤b2中Wecv代表列车在缓和曲线上超高引起的垂向内侧限界加宽值,单位mm,1500为钢轨中心距,单位mm,Hec代表缓和曲线外轨超高值,单位mm,h代表限界控制点自轨面的高度,单位mm。
优选的,考虑缓和曲线上列车内侧限界总加宽值的计算方法,S3中Wi代表列车过缓和曲线时总的内侧加宽值,单位mm。
优选的,考虑缓和曲线外侧限界加宽的计算方法,其中:S4的步骤a4中定义g(xg,yg)、e(xe,ye)、h(xh,yh)、f(xf,yf)四点,S4的步骤b4中的lef代表列车总车长,S4的步骤c4中leg代表列车车尾外侧加宽值,lea代表车为到后转架的距离,lbf代表车头到前转向架的距离,lhf代表列车车头外侧加宽值。
与现有的技术相比,本发明的有益效果是:
此种算法可提高限界加宽计算值的精度且能在保证限界安全的前提下不会有过多的冗余,能够有效的控制地铁施工的成本。
附图说明
图1为本发明实施例中车体在缓和曲线段的拟合示意图;
图2为本发明实施例中车体在缓和曲线段的简化图;
图3为本发明实施例中车体在缓和曲线段的拟合简化图;
图4为本发明实施例中列车在缓和曲线处的内侧加宽曲线;
图5为本发明实施例中列车在缓和曲线处的外侧侧加宽曲线;
具体实施方式
下面将结合本发明实施例,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例一:以总车长22.1m、转向架定距为15.7m的地铁A型车,以总长为50m,目标圆曲线半径为200m的缓和曲线工况为实施例对本发明专利进行详细的阐述说明,此工况下的Hmax=120mm,通过确定A型车危险点,h取3231mm,包括以下步骤:
S1,考虑车体几何偏移产生的内侧限界加宽计算方法:
对列车过缓和曲线时的限界曲线加宽进行计算,构造x、y轴坐标系。如附图1所示,直缓(ZH)点为坐标系原点;曲线为拟合后的缓和曲线;长方形ABCD表示地铁A型车车体定义虚线Lab为车体函数;Lab与曲线相交的动b(xb,yb)、a(xa,ya)分别代表列车的前后转向架中心点。
根据《地铁设计规范》规定:线路平面圆曲线与直线之间应根据曲线半径、超高设置及设计速度等因素设置缓和曲线,常用的缓和曲线为三次方抛物线型起方程为:
Figure BDA0003119701010000051
式中,L是缓和曲线总长度,单位m;R是圆曲线半径,单位m。
又已知曲率公式:
Figure BDA0003119701010000052
如附图2所示,在已知车辆后转向架中心a点坐标(xa,ya)的情况下,求车辆对应的几何偏移产生的限界加宽值fec,求解思路如下:
根据实施例工况可得缓和曲线的函数表达式为
Figure BDA0003119701010000053
前后转向架都落在缓和曲线上因此可设
Figure BDA0003119701010000054
依据地铁A型车的标准尺寸转向架定距lab=15.7m,可知a、b是满足直线距离恒等于15.7m的两个点;通过这一条件可以确定动点a、b之间的关系式
Figure BDA0003119701010000055
在已知a点坐标(xa,ya)的情况下可以用a、b之间的关系式求解出动点b的坐标,解出点b的坐标后,再列出缓和曲线段A型车的车体函数Lab如下:
Figure BDA0003119701010000056
接着建立动点a(xa,ya)于列车内侧最大偏移点c(xm,ym)之间的关系,已知lab=15.7m,内侧最大偏移点c点同时也是列车中点,因此有:
Figure BDA0003119701010000057
根据距离公式,可列出:
Figure BDA0003119701010000058
得c(xm,ym)的坐标值。
在已知c(xm,ym)的坐标值的前提下,带入c点坐标对缓和曲线三次方抛物线求Kec值,可得:
Figure BDA0003119701010000059
式中,xm为c点的横坐标;Kec为缓和曲线对应c点的曲率。
缓和曲线是指设置在直线与圆曲线之间的一种曲率连续变化的曲线。缓和曲线曲率Kec的值从零平滑且连续的增长到目标圆曲线的曲率值,可有效的减少因为曲率突变对列车运行带来的不利影响。圆曲线曲率值
Figure BDA00031197010100000510
因此缓和曲线Kec值的取值范围在
Figure BDA00031197010100000511
Figure BDA00031197010100000512
因此,缓和曲线可视为由不同曲率半径的圆弧组成的平面线型,利用曲率公式可对缓和曲线的任意一点求得相对应的Kec值,缓和曲线的曲率半径
Figure BDA00031197010100000513
可建立曲率半径Rec随着c点的横坐标值xm的变化而变化的关系式,以三次抛物线为函数分析其值等于:
Figure BDA0003119701010000061
式中,xm为c点的横坐标;Rec为缓和曲线的曲率半径。
公式(5)中,E为定值代表圆曲线的半径,圆曲线的K值也是固定不变的,即
Figure BDA0003119701010000062
但所求的缓和曲线曲率Kec则在
Figure BDA0003119701010000063
区间里不断变化,因而缓和曲线相对应的Rec值也在不停的变化,列车运行在圆曲线上时K值不变。因此,可以使用已知的正矢公式
Figure BDA0003119701010000064
对列车处在圆曲线上因几何偏移造成的限界加宽做分析,然而在缓和曲线上则需要额外考虑变量Rec带来的影响,需联立正矢公式和曲率公式才能计算出缓和曲线上任意一点的正矢值fec,以缓和曲线函数为
Figure BDA0003119701010000065
为例,可得出任意点在曲线
Figure BDA0003119701010000066
的正矢值求解方法如下:
Figure BDA0003119701010000067
式中,l为车辆转向架定距,单位m;xm为c点的横坐标;fec为列车在缓和曲线上因几何偏移引起的内侧加宽值,单位mm;
如公式(6)所示,fec的大小与c点位置以及缓和曲线的表达式有关,在缓和曲线变化比较陡峭(即曲率大,曲率半径R小)时计算出的fec值就会变大,在缓和曲线变化较为平缓(即曲率小,曲率半径R大)时fec值就小。
S2,考虑缓和曲线线路超高产生的垂向内侧限界加宽值:
考虑缓和曲线上超高对限界加宽的影响,已知计算圆曲线超高的公式(7),车辆在圆曲线上垂向内侧限界加宽公式:
Figure BDA0003119701010000068
Figure BDA0003119701010000069
式中,Hmax为目标圆曲线外轨超高值此工况为120,单位mm;V为过圆曲线时列车的设计速度,单位m/s;R为圆曲线的半径,单位m;h为限界控制点自轨面的高度,单位mm;1500为钢轨中心距,单位mm;Wv为圆曲线超高引起的垂向内侧限界加宽值,单位mm。
现求解因缓和曲线线路超高产生的车辆垂向内侧加宽值Wecv
Figure BDA00031197010100000610
式中,Wecv为列车在缓和曲线上超高引起的垂向内侧限界加宽值,单位mm;1500为钢轨中心距,单位mm;Hec为缓和曲线外轨超高值,单位mm;3231为限界控制点自轨面的高度,单位mm。
在圆曲线线路上超高为定值,而缓和曲线的超高Hec则为一个变量,Hec的值在0到Hmax的范围内均匀的变化,减少超高的突变对列车运行造成的不利影响,即Hec∈[0,Hmax]。Hec的计算方法如下:
Figure BDA00031197010100000611
Figure BDA00031197010100000612
式中,L为从直缓点算起至列车行驶到缓和曲线上任意一点的距离,单位m;L为缓和曲线的总长度,单位m;xb为列车前转向架中心点b的横坐标;120mm为目标圆曲线外轨超高值。
由于(11)式不可积分,因此用幂级数展开成多项式,舍去高次项,再进行积分,可得:
Figure BDA0003119701010000071
缓和曲线的超高是渐变的,当列车车辆行驶时,前转向架中心点b点的超高最大因此需要带入前转向架中心b点的横坐标xb来计算;根据公式(12),考虑xb,再带入原函数
Figure BDA0003119701010000072
做积分,就可以求出L,最后带入公式(10)即可求出Hec的值。
S3,缓和曲线下总内侧限界加宽计算方法,计算公式如下:
Figure BDA0003119701010000073
式中,Wi为列车过缓和曲线时总的内侧加宽值,单位mm;fec为列车在缓和曲线上因几何偏移引起的内侧加宽值,单位mm;Wecv为列车在缓和曲线上超高引起的垂向内侧限界加宽值;单位mm。
S4,缓和曲线外侧限界加宽计算方法:
考虑的是对列车过缓和曲线求解外侧加宽值的方法,定义g(xg,yg)、e(xe,ye)、h(xh,yh)、f(xf,yf)四点,如附图3所示。
根据车辆后转向架中心点所在坐标a(xa,ya)求车体的最大外侧曲线加宽值,求解的思路如下:
第一步,设点a(xa,ya)为列车后转向架中心点,b(xb,yb)为列车前转向架中心点,地铁A型车总车长lef=22.1m、转向架定距lab=15.7m,可知lea=3.2m,设动点e(xe,ye)、f(xf,yf);
第二步,求列车车尾外侧曲线加宽:
由lea=3.2m可列长度方程
Figure BDA0003119701010000074
在已知a点坐标的情况下,很容易算出e点坐标(xe,ye),再做垂直与Lef的直线Leg
Figure BDA0003119701010000075
联立直线Leg和缓和曲线方程即可求出g点坐标(xg,yg),最后计算出来的
Figure BDA0003119701010000076
就是列车车尾外侧曲线加宽。
第三步,求列车车头外侧曲线加宽:
由lbf=3.2m可列长度方程
Figure BDA0003119701010000077
在已知b点坐标的情况下,很容易算出f点坐标(xf,yf),做垂直与Lef的直线Lhf
Figure BDA0003119701010000078
联立Lhf和圆曲线方程即可求出h点坐标(xh,yh),最后计算出来的
Figure BDA0003119701010000079
就是列车车尾外侧曲线加宽;
Figure BDA00031197010100000710
Figure BDA0003119701010000081
式中,W车头为列车行驶在缓和曲线处车头的外侧曲线加宽值,单位mm;W车尾为列车行驶在缓和曲线处车尾的外侧曲线加宽值,单位mm。
S5,收集S1、S2、S3、S4中得到的车尾最大偏移点e(xe,ye)、车头最大偏移点f(xf,yf)的坐标、内侧最大偏移点c(xm,ym)的坐标。最终求出在缓和曲线段车体函数内外侧加宽函数曲线,如附图4、5所示。
本发明优点是:根据经验,地铁限界尺寸毫米级的优化或者改动,往往会引起土建工程费用千万元级的变化。所以,确定一个既能保证车辆运行安全,又不增大隧道空间的经济、合理的断面尺寸,是地铁限界设计的任务和目的,后期可以编写程序带入不同列车过缓和曲线工况数据,利用计算机进自动计算,可以节省大量的人力、物力,缩短设计周期,提高设计、计算精度。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。

Claims (5)

1.一种地铁线路缓和曲线段的加宽曲线计算方法,其特征在于,包括以下几个步骤:
S1,考虑车体因几何偏移产生的内侧限界加宽,其特征在于,包括以下几个小的步骤:
a1.任意选取一种地铁车辆车型,根据车型的转向架定距i以及车长I设定车体函数表达式;
b1.任意选取一种地铁缓和曲线的工况,根据《地铁设计规范》规定:线路平面圆曲线与直线之间应根据曲线半径、超高设置及设计速度等因素设置缓和曲线,常用的缓和曲线为三次方抛物线型起方程为:
Figure FDA0003119701000000011
式中,R为目标圆曲线的半径,单位m,L为缓和曲线的总长度,单位m,用于确定列车前后转向架的函数关系;
c1.建立一个直角坐标系,定义坐标系原点为直缓(ZH)点(直线和缓和曲线的切点,标注时用ZH表示);直缓(ZH)点到圆缓(YH)点(圆缓点是缓和曲线的终点和下一段圆曲线的起点相交的点)之间的曲线为缓和曲线;圆缓(YH)点之后为圆曲线;长方形ABCD表示地铁A型车车体;定义虚线Lab为车体函数,Lab与曲线相交的动点b(xb,yb)、a(xa,ya)分别代表列车的前后转向架中心点;
d1.根据步骤a、b的求解车体函数的表达式,用直线Lab表示简化的车体,此时前后转向架都落在缓和曲线上,因此,可设
Figure FDA0003119701000000012
转向架定距lab=im,可知a、b是满足直线距离恒等于im的两个点,通过这一条件可以确定动点a、b之间的关系式:
Figure FDA0003119701000000013
在已知a点坐标(xa,ya)的情况下,可以通过计算a、b之间的关系式解出动点b的坐标,最后利用b的坐标,列出缓和曲线段列车的车体函数:Lab
Figure FDA0003119701000000014
e1.求解车体最大内侧偏移点c的坐标值,建立动点a(xa,ya)与列车内侧最大偏移点c(xm,ym)之间的关系,已知lab=im,内侧最大偏移点c点同时也是列车中点,因此有:
Figure FDA0003119701000000015
根据距离公式,可列出最大偏移点c点与车辆后转向架中心点a的关系式:
Figure FDA0003119701000000016
在已知a点坐标的基础上就可以求解c(xm,ym)的坐标值;
f1.求解缓和曲线的曲率半径,现已知c(xm,ym)的坐标值,带入c点坐标对缓和曲线三次方抛物线求Kec值,可得:
Figure FDA0003119701000000017
式中,xm为c点的横坐标,Kec为缓和曲线对应c点的曲率,缓和曲线的曲率半径
Figure FDA0003119701000000018
因此,可建立曲率半径Rec随着c点的横坐标值xm的变化而变化的关系式;
以三次抛物线为函数分析其值等于:
Figure FDA0003119701000000019
式中,L是缓和曲线总长度,单位m,R是圆曲线半径,单位m,xm为c点的横坐标,Rec为缓和曲线的曲率半径,单位m;
g1.求解列车在缓和曲线上因几何偏移引起的内侧加宽值,使用正矢公式
Figure FDA00031197010000000110
可以对圆曲线上列车因几何偏移造成的加宽值做计算,但在缓和曲线段则需联立正矢公式和曲率公式计算缓和曲线上任意一点的正矢值fec
以拟合的缓和曲线函数为
Figure FDA00031197010000000111
为例,可到出处在于缓和曲线
Figure FDA00031197010000000112
上任意一点的正矢值求解方法:
Figure FDA00031197010000000113
式中,R为目标圆曲线的半径,单位m,L为缓和曲线长度,单位m,lab为车辆转向架定距,单位m,xm为c点的横坐标,fec为列车在缓和曲线上因几何偏移引起的内侧加宽值,单位mm;
S2,考虑缓和曲线线路超高产生的垂向内侧限界加宽值,其特征在于,包括以下小步骤:
a2.计算缓和曲线的超高Hec,缓和曲线的超高Hec为变量,Hec的值在0到Hmax的范围内均匀的变化,即Hec∈[0,Hmax];
Hec的计算方法:
Figure FDA0003119701000000021
其中:
Figure FDA0003119701000000022
其中:
Figure FDA0003119701000000023
式中,L为从直缓点算起至列车行驶到缓和曲线上任意一点的距离,单位m,L为缓和曲线的总长度,单位m,xb为列车前转向架中心点b的横坐标,xYH为圆缓点的横坐标;
由于积分为不可积积分,因此用幂级数展开成多项式,舍去高次项,再进行积分,可得:
Figure FDA0003119701000000024
Figure FDA0003119701000000025
前转向架中心点b点的超高最大因此需要带入前转向架中心b点的横坐标xb来计算,考虑xb、xYH,再带入原函数
Figure FDA0003119701000000026
做积分,就可以求出L、L、Hec的值;
b2.计算缓和曲线上超高引起的垂向内侧限界加宽值Wecv
Figure FDA0003119701000000027
式中,Wecv为列车在缓和曲线上超高引起的垂向内侧限界加宽值,单位mm,1500为钢轨中心距,单位mm,Hec为缓和曲线外轨超高值,单位mm,h为限界控制点自轨面的高度,单位mm;
S3,考虑缓和曲线上列车内侧限界总加宽,其特征在于,包括以下公式:
缓和曲线段内侧限界总加宽计算公式为:
Figure FDA0003119701000000028
式中,Wi为列车过缓和曲线时总的内侧加宽值,单位mm,fec为列车在缓和曲线上因几何偏移引起的内侧加宽值,单位mm,Wecv为列车在缓和曲线上超高引起的垂向内侧限界加宽值,单位mm,h为限界控制点自轨面的高度,单位mm;
S4,考虑缓和曲线外侧限界加宽,其特征在于,包括以下小步骤:
a4.对列车过缓和曲线求解外侧加宽值的方法,定义g(xg,yg)、e(xe,ye)、h(xh,yh)、f(xf,yf)四点;
b4.设点a(xa,ya)为列车后转向架中心点,b(xb,yb)为列车前转向架中心点,地铁A型车总车长lef=Im、转向架定距lab=im,可知
Figure FDA0003119701000000029
设动点e(xe,ye)、f(xf,yf);
c4.求列车车尾外侧曲线加宽:
Figure FDA00031197010000000210
可列长度方程
Figure FDA00031197010000000211
在已知a点坐标的情况下,很容易算出e点坐标(xe,ye),再做垂直与Lef的直线Leg
Leg
Figure FDA00031197010000000212
联立直线Leg和缓和曲线方程即可求出g点坐标(xg,yg),最后计算出来的
Figure FDA0003119701000000031
就是列车车尾外侧曲线加宽;
d4.求列车车头外侧曲线加宽:
Figure FDA0003119701000000032
可列长度方程
Figure FDA0003119701000000033
在已知b点坐标的情况下,很容易算出f点坐标(xf,yf),再做垂直与Lef的直线Lhf
Lhf
Figure FDA0003119701000000034
联立Lhf和圆曲线方程即可求出h点坐标(xh,yh),最后计算出来的
Figure FDA0003119701000000035
就是列车车尾外侧曲线加宽;
其中:
Figure FDA0003119701000000036
其中:
Figure FDA0003119701000000037
式中,W车头为列车行驶在缓和曲线处车头的外侧曲线加宽值,单位mm,W车尾为列车行驶在缓和曲线处车尾的外侧曲线加宽值,单位mm;
S5,用最小二乘法拟合列车过缓和曲线时的加宽函数,其特征在于,包括以下步骤:
a5.确定列车后转向架xa的取值范围,设后转向架中心点为点a(xa,ya),那么(0≤xa<xYH);
b5.对xa在(0≤xa<xYH)平均取点,以这些点为基准,求解列车相对应的内侧最大偏移点c(xm,ym)的坐标;
c5.对xa在(0≤xa<xYH)平均取点,以这些点为基准,求出相对应的车尾最大偏移点e(xe,ye)、车头最大偏移点f(xf,yf)的坐标;
d5.用最小二乘法对S1、S2、S3、S4中所求的点进行函数拟合,最终求出在缓和曲线段车体函数内外侧加宽函数曲线。
2.根据权利要求1所述的考虑车体因几何偏移产生的内侧限界加宽,其特征在于:S1的步骤a1中i为车辆转向架定距长度,单位m,I为车体车长长度,单位m,此算法适应与任意车型的限界加宽计算,S1的步骤b1中R代表缓和曲线后圆曲线的半径,单位m,L代表缓和曲线总长度,单位m,S1的步骤c1中动点b代表列车前转向架中心点,动点a代表列车后转向架中心点,Lab代表车体函数,此算法适应与任意铁路缓和曲线工况下的计算,其特征在于建立直角坐标系对S1的步骤a1、b1的缓和曲线进行拟合,建立车体函数便于监测车体的运动轨迹,S1的步骤d1其特征在于以列车后转向架中心点a为已知的基准点对车体函数进行分析,S1的步骤e1中lab代表转向架定距,c(xm,ym)代表列车内侧最大偏移点,S1步骤中的xm代表c点的横坐标,Kec代表缓和曲线对应c点的曲率,Rec代表缓和曲线对应c点的曲率半径,m-1,S1的步骤f1其特征在于找出a点与c点的函数关系,S1的步骤g1中fec表示列车在缓和曲线上因几何偏移引起的内侧加宽值,单位mm。
3.根据权利要求1所述的考虑缓和曲线线路超高产生的垂向内侧限界加宽值,其特征在于:S2的步骤a2中Hec代表缓和曲线外轨超高值,单位mm,L为从直缓点算起至列车行驶到缓和曲线上任意一点的距离,单位m,L为缓和曲线的总长度,单位m,xb为列车前转向架中心点b的横坐标,xYH为圆缓点的横坐标,其特征在于利用微积分等知识计算变量Hec,S2的步骤b2中Wecv代表列车在缓和曲线上超高引起的垂向内侧限界加宽值,单位mm,1500为钢轨中心距,单位mm,Hec代表缓和曲线外轨超高值,单位mm,h代表限界控制点自轨面的高度,单位mm。
4.根据权利要求1所述的考虑缓和曲线上列车内侧限界总加宽,其特征在于:S3中Wi代表列车过缓和曲线时总的内侧加宽值,单位mm。
5.根据权利要求1所述的考虑缓和曲线外侧限界加宽,其特征在于:S4的步骤a4中定义g(xg,yg)、e(xe,ye)、h(xh,yh)、f(xf,yf)四点,S4的步骤b4中的lef代表列车总车长,S4的步骤c4中leg代表列车车尾外侧加宽值,lea代表车为到后转架的距离,lbf代表车头到前转向架的距离,lhf代表列车车头外侧加宽值。
CN202110673402.8A 2021-06-17 2021-06-17 地铁线路缓和曲线段的加宽曲线计算方法 Pending CN115495831A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110673402.8A CN115495831A (zh) 2021-06-17 2021-06-17 地铁线路缓和曲线段的加宽曲线计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110673402.8A CN115495831A (zh) 2021-06-17 2021-06-17 地铁线路缓和曲线段的加宽曲线计算方法

Publications (1)

Publication Number Publication Date
CN115495831A true CN115495831A (zh) 2022-12-20

Family

ID=84464181

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110673402.8A Pending CN115495831A (zh) 2021-06-17 2021-06-17 地铁线路缓和曲线段的加宽曲线计算方法

Country Status (1)

Country Link
CN (1) CN115495831A (zh)

Similar Documents

Publication Publication Date Title
CN105551251B (zh) 一种无信号交叉口机动车冲突概率判断方法
Guo et al. Numerical study for the aerodynamic performance of double unit train under crosswind
CN103693040B (zh) 一种基于双模式协作的车辆避撞系统
CN109214554A (zh) 一种不落轮工况下的铁道车辆轮对整车镟修策略决策算法
Zhang et al. Crosswind stability of high-speed trains in special cuts
CN114332212B (zh) 基于车载移动激光点云的轨道超高和前后高低检测方法
CN106951657A (zh) 一种磨损钢轨打磨目标廓形快速设计方法
CN105512397A (zh) 一种独立车轮踏面外形设计方法及独立车轮
CN113771871B (zh) 一种车辆偏离车道的预警方法及其系统、介质、车辆
CN109407658A (zh) 无人车的行车轨迹规划方法及计算机可读存储介质
CN108763806B (zh) 高速列车过特长双线隧道时列车及隧道表面所受静压的数值模拟分析方法
Cui et al. Multi-objective optimization of electric multiple unit wheel profile from wheel flange wear viewpoint
CN112818612A (zh) 基于隧道口行车安全仿真研究的安全管控措施确定方法
Josefsson et al. Numerical analysis of aerodynamic impact on passenger vehicles during cornering
CN115495831A (zh) 地铁线路缓和曲线段的加宽曲线计算方法
CN106844862B (zh) 一种基于cae分析的铝制车身节点刚度评价方法
CN106169000B (zh) 一种转弯弯道导轨及其设计方法
CN115495832A (zh) 地铁线路直缓段处的加宽曲线计算方法
CN101906758A (zh) 高速铁路中声屏障脉动力大小的确定方法及应用
CN113722812B (zh) 一种地铁圆缓段加宽的计算方法
Lei et al. The Effect of Lateral Offset Distance on the Aerodynamics and Fuel Economy of Vehicle Queues.
CN113722811B (zh) 一种列车明线交会时压力波幅值与车速比关系的估算方法
Wang et al. Research on aerodynamic characteristic for EMU passing by windbreak wall gap under crosswind
Ji et al. Aerodynamic Shape Design of Pantograph Network Monitoring Device on High-Speed Trains
Kamal et al. Effect of Crosswind on Aerodynamic Characteristics of a Generic Train Model using ANSYS

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination