CN113722812A - 一种地铁圆缓段加宽的计算方法 - Google Patents

一种地铁圆缓段加宽的计算方法 Download PDF

Info

Publication number
CN113722812A
CN113722812A CN202110805534.1A CN202110805534A CN113722812A CN 113722812 A CN113722812 A CN 113722812A CN 202110805534 A CN202110805534 A CN 202110805534A CN 113722812 A CN113722812 A CN 113722812A
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.)
Granted
Application number
CN202110805534.1A
Other languages
English (en)
Other versions
CN113722812B (zh
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 CN202110805534.1A priority Critical patent/CN113722812B/zh
Publication of CN113722812A publication Critical patent/CN113722812A/zh
Application granted granted Critical
Publication of CN113722812B publication Critical patent/CN113722812B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

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

Abstract

一种地铁圆缓段加宽的计算方法,所述方法对铁路缓和曲线、目标圆曲线线路进行函数拟合;求解列车过缓和曲线时的内外侧加宽点;利用最小二乘法对各点坐标进行拟合求出列车加宽曲线。所述方法包括以下步骤:(1)求解列车在缓和曲线段因几何偏移造成的内侧曲线加宽值;(1)求解列车在缓和曲线段因几何偏移造成的内侧曲线加宽值;(3)求解列车在缓和曲线段因几何偏移造成的内侧曲线加宽值;(4)求解列车在缓和曲线段的外侧曲线加宽值;(5)求解缓和曲线段列车车体函数内外侧加宽函数曲线。本发明适用于任意工况下铁路限界的设计、建设;特别对于应对城市轨道交通小半径、多曲线的复杂工况能够保证良好的精度。

Description

一种地铁圆缓段加宽的计算方法
技术领域
本发明涉及一种地铁圆缓段加宽的计算方法,属地铁线路设计技术领域。
背景技术
随着我国城市规模和经济建设的快速发展、城市化进程逐步加快,地铁车辆运营速度的提高成为了城市轨道交通的发展趋势;与此同时,城市轨道交通小半径曲线也越来越多这对车辆的限界安全有直接的影响;面对变化越来越复杂的城市轨道线路和运营速度的提高。
目前在建设施工中常用到的设计方法已经无法满足城市轨道交通对限界设计的需求,现有的圆缓段限界加宽计算方法基本采用三种较为常规方法:第一种是先计算出圆曲线范围内的加宽值,然后将这个加宽值沿缓和曲线有规则的分配,这种计算方法的缺陷是精度不高会使得限界设计造成过度的冗余存在着工程浪费现象;第二种是采用图解法使用CAD软件绘制出缓和曲线、圆曲线的图像再模拟列车过曲线的过程;第三种是采用几何分析法,针对线路圆缓段的过渡区段,提出一种基于线路中心线的精确加宽量计算方法。
上述计算方法中有些含有大量三角函数,增加了算法的复杂性;亦或者没有建立车体函数,不便于动态实时的监测列车的运动变化且在一些特殊点上不适用;依靠CAD软件的方法实际运算也是非常复杂的,且对软件设备具有较高的要求。
发明内容
本发明的目的是,针对现有技术对圆缓段限界加宽计算方法存在的问题,提出一种地铁圆缓段加宽的计算方法,
本发明实现的技术方案如下,一种地铁圆缓段加宽的计算方法,所述方法对铁路缓和曲线、目标圆曲线线路进行函数拟合;求解列车过缓和曲线时的内外侧加宽点;利用最小二乘法对各点坐标进行拟合求出列车加宽曲线;
所述方法包括以下步骤:
(1)考虑车体因几何偏移产生的内侧限界加宽,计算车体几何偏移产生的内侧曲线加宽值。
任意选取一种地铁车辆车型;根据车型的转向架定距i以及车长l设定车体函数表达式;
任意选取一种地铁缓和曲线的工况;根据缓和曲线的总长、线性、最大半径拟合出缓和曲线的函数表达式
Figure BDA0003166378530000021
式中,R表示目标圆曲线半径,单位m;L表示缓和曲线的总长度,单位m;根据工况拟合圆曲线表达式,设圆曲线表达式为(x-x1)2+(y-y1)2=R2;用于确定列车前后转向架的函数关系。
求车体函数Lab的函数表达式,动点b(xb,yb)代表列车前转向架中心点;动点a(xa,ya)代表列车后转向架中心点。
设动点a落在缓和曲线上、动点b落在圆曲线上,由拟合的缓和曲线和圆曲线表达式可设动点a的坐标为:xa,
Figure BDA0003166378530000022
动点b的坐标为:xb,
Figure BDA0003166378530000023
列车转向架定距为im,因此线段lab恒等于im,根据两点距离公式,可列出动点a和动点b之间的距离方程式:
Figure BDA0003166378530000031
其中,动点b(xb,yb)代表列车前转向架中心点,动点a(xa,ya)代表列车后转向架中心点;xa为动点a的x座标值;ya为动点a的y座标值;xb为动点b的x座标值;yb动点b的y座标值;R表示目标圆曲线半径;L表示缓和曲线的总长度;x1为圆曲线圆心的横坐标;y1为圆曲线圆心的纵坐标。
已知曲线函数和列车后转向架中心点a的坐标后,可以通过距离方程式确定列车前转向架中心点b点的坐标值以及动点a、b之间的函数关系,在求出列车前后转向架中心点b、a的坐标后,可列出车体函数Lab的表达式:
Figure BDA0003166378530000032
通过车体函数Lab斜率并确定辅助线Lcd的斜率;从车体函数Lab的表达式中可知斜率
Figure BDA0003166378530000033
所构造的辅助线Lcd与车体函数Lab恒垂直,可求得Lcd的斜率:
Figure BDA0003166378530000034
Lcd代表计算内侧曲线加宽值所拟定的辅助线,垂直于车体函数Lab通过车体函数Lab的斜率,以及两者的垂直关系可以得到辅助线Lcd的斜率。
求lcd的最大值,求取最大值的过程关键在于确定动点d的坐标,d点为落在缓和曲线或圆曲线上的一个动点,对缓和曲线或圆曲线求导找出满足y′d=Kab的动点d,满足条件的动点d即是lcd取到最大值的位置,此时会出现三种情况:
情况1:Kab<y′圆缓点,说明车体内侧最大偏移点未过圆缓点处于缓和曲线上,因此,求解动点d时应对缓和曲线进行求导,
情况2:Kab=y′圆缓点,说明车体内侧最大偏移点落在圆缓(YH)点上,缓和曲线和圆曲线在圆缓点处曲率几乎相等,因此,带入任意一种函数都可以进行计算;
情况3:Kab>y′圆缓点,说明车体内侧最大偏移点越过圆缓点处于圆曲线上,因此,求解动点d时应对圆曲线进行求导即:
Figure BDA0003166378530000041
式中,y′圆曲线为圆曲线的导数,y′缓和曲线为缓和曲线的导数,y′圆缓点为圆缓点的导数值,Kab为车体函数Lab的斜率。
此处只针对Kab<y′圆缓点的情况1做分析,情况3只需把缓和曲线函数换成圆曲线函数重新求导计算即可。
Figure BDA0003166378530000042
可得xd用a、b坐标表示的解析式:
Figure BDA0003166378530000043
xd取正值,可用(xa,ya)、(xb,yb)来表示d点坐标:
Figure BDA0003166378530000044
算得动点d的坐标和Lcd的斜率Kcd后,可求出线段Lcd的表达式:
Figure BDA0003166378530000051
求解车辆内侧偏移点c点的坐标,联立Lcd与Lab的函数表达式,求解方程式:
Figure BDA0003166378530000052
得:
Figure BDA0003166378530000053
又因内侧偏移点c点处在车体函数Lab上,把xc带入Lab的方程式即可算得:
Figure BDA0003166378530000054
最后对c点和d点坐标使用一次两点距离方程即可求出Lcd的值:
Figure BDA0003166378530000055
式中,lcd为车体几何偏移产生的内侧曲线加宽值,单位mm。
(2)考虑缓和曲线线路超高产生的垂向内侧限界加宽值,计算圆缓段处因线路超高产生的加宽值。
列车处于圆缓段时,前转向架中心点b点已经落在了圆曲线上,因此在圆缓段的超高应使用圆曲线的超高:
Figure BDA0003166378530000061
式中,Hyh为计算圆缓段曲线加宽用到的曲线超高值,单位mm;Hmax为目标圆曲线的超高值,单位mm,V为过圆曲线时列车的设计速度,单位m/s;R为目标圆曲线半径,单位m;力代表限界控制点自轨面的高度,单位mm;1500为钢轨中心距,单位mm;Wyhv为圆缓段处因线路超高产生的加宽值,单位mm。
(3)考虑缓和曲线上列车内侧限界总加宽,计算列车行驶在缓和曲线处车头和车尾的内侧曲线加宽值。
Figure BDA0003166378530000062
式中,Wi为圆缓段处车辆总的内侧加宽值,单位mm。
(4)考虑缓和曲线外侧限界加宽,计算列车行驶在缓和曲线处车头和车尾的外侧曲线加宽值。
对列车过缓和曲线求解外侧加宽值的方法,定义g(xg,yg)、e(xe,ye)、h(xh,yh)、f(xf,yf)四点;
设点a(xa,ya)为列车后转向架中心点,b(xb,yb)为列车前转向架中心点,地铁A型车总车长lef=Im、转向架定距lab=im,可知
Figure BDA0003166378530000063
设动点e(xe,ye)、f(xf,yf)。
求列车车尾外侧曲线加宽:
Figure BDA0003166378530000071
可列长度方程
Figure BDA0003166378530000072
在已知a点坐标的情况下,很容易算出e点坐标(xe,ye),再做垂直与Lef的直线Leg
Figure BDA0003166378530000073
联立直线Leg和缓和曲线方程即可求出g点坐标(xg,yg),最后计算出来的
Figure BDA0003166378530000074
就是列车车尾外侧曲线加宽。
求列车车头外侧曲线加宽:
Figure BDA0003166378530000075
可列长度方程
Figure BDA0003166378530000076
在已知b点坐标的情况下,很容易算出f点坐标(xf,yf),做垂直与Lef的直线Lhf
Figure BDA0003166378530000077
联立Lhf和圆曲线方程即可求出h点坐标(xh,yh),最后计算出来的
Figure BDA0003166378530000078
就为列车车尾外侧曲线加宽。
Figure BDA0003166378530000079
Figure BDA00031663785300000710
式中,W车头为列车行驶在缓和曲线处车头的外侧曲线加宽值,单位mm;W车尾为列车行驶在缓和曲线处车尾的外侧曲线加宽值,单位mm。
(5)用最小二乘法拟合列车过缓和曲线时的加宽函数,步骤如下:
1)确定列车后转向架xa的取值范围,设后转向架中心点为点a(xa,ya),当x圆缓点≤xa时列车就完全进入了圆曲线,因此,(x圆缓点-i≤xa≤x圆缓点);
2)对xa在(x圆缓点-i≤xa≤x圆缓点)平均取点,以这些点为基准,求解列车相对应的内侧最大偏移点c(xc,yc)的坐标;
3)对xa在(x圆缓点-i≤xa≤x圆缓点)平均取点,以这些点为基准,求出相对应的车尾最大偏移点e(xe,ye)、车头最大偏移点f(xf,yf)的坐标;
4)用最小二乘法对步骤(1)、(2)、(3)、(4)中所求的点进行函数拟合,最终求出在圆缓段处车体函数内外侧加宽函数曲线。
本发明的有益效果在于,本发明在坚实的理论研究基础上开发了一种精确、系统、使用方便的缓和曲线地段限界加宽计算方法,本发明适用于任意工况下铁路限界的设计、建设;特别对于应对城市轨道交通小半径、多曲线的复杂工况能够保证良好的精度。本发明使用微积分、最小二乘法拟合曲线以及导函数的应用通过建立车体函数,研究车体过圆缓段的运行过程,并得出在此过程中,列车需加宽限界的变化规律,针对列车过缓和曲线归纳出内侧和外侧的加宽计算方法能大大节省设计者时间和精力。
根据经验,地铁限界尺寸毫米级的优化或者改动,往往会引起土建工程费用千万元级的变化。所以,确定一个既能保证车辆运行安全,又不增大隧道空间的经济、合理的断面尺寸,是地铁限界设计的任务和目的,本发明对后期编写程序带入不同列车过缓和曲线工况数据,利用计算机进自动计算,可以节省大量的人力、物力,缩短设计周期,提高设计、计算精度。
附图说明
图1为车体在圆缓段上的拟合示意图;
图2为车体在圆缓段上的简化图;
图3为车体在圆缓段的拟合简化图;
图4为列车在圆缓段曲线处的外侧加宽曲线;
图5为列车在圆缓段曲线处的内侧加宽曲线图;
图6为本实施例一种地铁圆缓段加宽的计算方法流程图。
具体实施方式
本实施例一种地铁圆缓段加宽的计算方法,如图6所示。
本实施例选用总车长度为22.1m、转向架定距为15.7m的地铁A型车;选用总长为50m、目标圆曲线半径为200m的缓和曲线工况。为实施例对本发明进行详细的阐述说明,此工况下的Hmax=120mm、通过确定A型车危险点h值取3231mm,包括以下实施步骤:
S1、本实施例表述的是列车前转向架中心点b落在圆曲线上、后转向架中心点a落在缓和曲线的工况,在此工况下,求解列车因几何偏移照成的内侧曲线加宽值:
为便于分析,构造坐标系;如图1所示,定义坐标系原点为直缓(ZH)点(直线和缓和曲线的切点,标注时用ZH表示);直缓(ZH)点到圆缓(YH)点(圆缓点是缓和曲线的终点和下一段圆曲线的起点相交的点)之间的曲线为缓和曲线;圆缓(YH)点之后为圆曲线;长方形ABCD表示地铁A型车车体;定义虚线Lab为车体函数;Lab与曲线相交的动点b、a分别代表列车的前后转向架中心点。
根据《地铁设计规范》规定:“线路平面圆曲线与直线之间应根据曲线半径、超高设置及设计速度等因素设置缓和曲线,常用的缓和曲线为三次方抛物线型起方程为:
Figure BDA0003166378530000101
其中,C=RL
式中,L是缓和曲线总长度,单位m;R是圆曲线半径,单位m。
假设现已知车辆后转向架中心a点坐标(xa,ya),以此为条件求解车辆对应的最大内侧加宽值lcd,求解思路如下:
首先,确定列车处于缓和曲线上后转向架中心点a的横坐标xa的取值范围,当34≤xa列车前转向架进入圆曲线,当50≤xa列车则完全进入圆曲线段,因此xa∈[34,50]。
接着,对前转向架中心点a落在圆曲线上,后转向架中心点b落在缓和曲线的圆缓段曲线加宽进行分析前,需要确定各点、曲线以及车体函数的表达式,对R=200m、L=50m的工况进行函数拟合,拟合缓和曲线的表达式为三次方函数:
Figure BDA0003166378530000102
拟合圆曲线表达式为函数:(x-21.21)2+(y-200)2=40000;动点a为列车后转向架、动点b为列车前转向架,依据地铁A型车的标准尺寸设线段lab=15.7m。
求线段lcd的最大长度即是车体最大的内侧偏移量,直线Lcd是计算内侧曲线加宽值所拟定的辅助线,车辆内侧偏移点c落在车体函数Lab上,动点d落在曲线上,辅助线Lcd垂直于车体函数Lab;如图2所示,Lcd在Lab与曲线形成的封闭平面内任意的移动,Lcd的长度随着列车后转向架中心点a的变化而变化,此过程中可取到无数条Lcd,对曲线的某一个点进行求导(此点为上述定义的动点d),当动点d的曲线导数值与车体函数Lab的斜率Kab相等时,就是Lcd取到最大值的位置。
最后,再考虑超高对车辆限界加宽带来的影响,具体的求解数理推导过程如下:
第一步,求出车体函数Lab的函数表达式。
由拟合曲线的表达式可知,动点a的坐标为
Figure BDA0003166378530000111
动点b的坐标为(xb
Figure BDA0003166378530000112
),已知地铁A型车转向架定距为15.7m,因此线段lab恒等于15.7m,根据两点距离公式,可列出动点a和动点b之间的距离方程式:
Figure BDA0003166378530000113
在已知曲线函数和列车后转向架中心点a的坐标后,可以用此公式确定列车前转向架b点的坐标值以及动点a、b之间的函数关系,为简化方程推导过程用yb代替较为复杂的下式,
Figure BDA0003166378530000114
函数关系如下:
Figure BDA0003166378530000115
化简可得:
Figure BDA0003166378530000116
把xb看作是未知量对其求解二元一次方程:
Figure BDA0003166378530000117
因xa<xb,可确定xa与xb之间的关系:
Figure BDA0003166378530000118
根据动点a、b的坐标,列出车体函数Lab的表达式:
Figure BDA0003166378530000121
第二步,通过车体函数Lab斜率,确定辅助线Lcd的斜率:
Figure BDA0003166378530000122
所构造的辅助线Lcd与车体函数Lab恒垂直,可求得Lcd的斜率:
Figure BDA0003166378530000123
第三步,求lcd的最大值,求取最大值的过程关键在于确定动点d的坐标,对缓和曲线或圆曲线求导找出满足y′d=Kab的动点d,满足条件的动点d即是lcd取到最大值的位置,此时会出现三种情况。
情况一:Kab<y′圆缓点,说明车体内侧最大偏移点未过圆缓点处于缓和曲线上,因此,求解动点d时应对缓和曲线进行求导;
情况二:Kab=y′圆缓点,说明车体内侧最大偏移点落在圆缓(YH)点上,缓和曲线和圆曲线在圆缓点处曲率几乎相等,因此,带入任意一种函数都可以进行计算;
情况三:Kab>y′圆缓点,说明车体内侧最大偏移点越过圆缓点处于圆曲线上,因此,求解动点d时应对圆曲线进行求导即:
Figure BDA0003166378530000124
Figure BDA0003166378530000125
式中,y′圆曲线为圆曲线的导数;y′缓和曲线为缓和曲线的导数;y′圆缓点为圆缓点的导数值;Kab为车体函数的斜率。
此处只针对Kab<y′圆缓点的情况一做分析,情况二只需把缓和曲线函数换成圆曲线重新求导计算即可,
Figure BDA0003166378530000131
可得xd用a、b坐标表示的解析式:
Figure BDA0003166378530000132
xd取正值,可用(xa,ya)、(xb,yb)来表示d点坐标:
Figure BDA0003166378530000133
算得动点d的坐标和Lcd的斜率Kcd后,可求出线段Lcd的表达式:
Figure BDA0003166378530000134
第四步,求解车辆内侧偏移点c点的坐标,Lcd与Lab的交点为c点,联立两直线方程式可求出c点的坐标:
Figure BDA0003166378530000135
解方程式可得c点坐标(xc,yc):
Figure BDA0003166378530000141
Figure BDA0003166378530000142
最后对c点和d点坐标使用一次两点距离方程即可求出Lcd的值:
Figure BDA0003166378530000143
式中,lcd为车体几何偏移产生的内侧曲线加宽值,单位mm。
S2、本实施例表述的是列车前转向架中心点落在圆曲线上、后转向架中心点落在缓和曲线的工况,在此工况下求解列车因线路超高照成的内侧曲线加宽值:
考虑超高对圆缓段车辆曲线内加宽带来的影响,列车处圆缓段时,前转向架中心点b点已经落在圆曲线上,因此在圆缓段的超高应使用圆曲线的超高。
Hyh=Hmax=120mm (15)
Figure BDA0003166378530000144
式中,Hyh为计算圆缓段曲线加宽用到的曲线外轨超高值,单位mm;Hmax为目标圆曲线的超高值,单位mm;V为过圆曲线时列车的设计速度,单位m/s;R为圆曲线的半径,单位m;h代表限界控制点自轨面的高度,单位mm;1500为钢轨中心距,单位mm;Wyhv为圆缓段处因线路超高产生的加宽值,单位mm。
S3、本实施例表述的是列车前转向架中心点落在圆曲线上、后转向架中心点落在缓和曲线的工况,在此工况下求解列车总的内侧曲线加宽值:
Figure BDA0003166378530000151
式中,Wi为圆缓段处,车辆总的内侧加宽值,单位mm。
S4、本实施例表述的是列车前转向架中心点落在圆曲线上、后转向架中心点落在缓和曲线的工况,在此工况下求解列车外侧曲线加宽值。
对列车前转向架落中心点a落在在圆曲线上、后转向架中心点b落在缓和曲线上的情况,求解外侧加宽值的方法,如图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 BDA0003166378530000152
在已知a点坐标的情况下,很容易算出e点坐标(xe,ye),再做垂直与Lef的直线Leg
Figure BDA0003166378530000153
联立直线Leg和缓和曲线方程即可求出g点坐标(xg,yg),最后计算出来的
Figure BDA0003166378530000154
就是列车车尾外侧曲线加宽。
第三步,求列车车头外侧曲线加宽:
由lbf=3.2m可列长度方程
Figure BDA0003166378530000161
在已知b点坐标的情况下,很容易算出f点坐标(xf,yf),做垂直与Lef的直线Lhf
Figure BDA0003166378530000166
联立Lhf和圆曲线方程即可求出h点坐标(xh,yh),最后计算出来的
Figure BDA0003166378530000163
就是列车车尾外侧曲线加宽。
Figure BDA0003166378530000164
Figure BDA0003166378530000165
式中,W车头为列车运行在圆缓段处车头处的外侧曲线加宽值,单位mm;W车尾为列车运行在圆缓段处车尾处的外侧曲线加宽值,单位mm。
S5、收集S1、S2、S3、S4中得到的车尾最大偏移点e(xe,ye)、车头最大偏移点f(xf,yf)的坐标、内侧最大偏移点c(xm,ym)的坐标。最终求出在缓和曲线段车体函数内外侧加宽函数曲线,如图4和图5所示。

Claims (7)

1.一种地铁圆缓段加宽的计算方法,其特征在于,所述方法对铁路缓和曲线、目标圆曲线线路进行函数拟合;求解列车过缓和曲线时的内外侧加宽点;利用最小二乘法对各点坐标进行拟合求出列车加宽曲线;
所述方法包括以下步骤:
(1)求解列车在缓和曲线段因几何偏移造成的内侧曲线加宽值;
(2)求解列车在缓和曲线段因几何偏移造成的内侧曲线加宽值;
(3)求解列车在缓和曲线段因几何偏移造成的内侧曲线加宽值:
(4)求解列车在缓和曲线段的外侧曲线加宽值:
(5)求解缓和曲线段列车车体函数内外侧加宽函数曲线。
2.根据权利要求1所述的一种地铁圆缓段加宽的计算方法,其特征在于,所述列车在缓和曲线段因几何偏移造成的内侧曲线加宽值求解如下;
求车体函数Lab的函数表达式,设动点a落在缓和曲线上、动点b落在圆曲线上,由拟合的缓和曲线和圆曲线表达式可设动点a的坐标为:xa
Figure RE-FDA0003283666800000011
动点b的坐标为:xb
Figure RE-FDA0003283666800000012
列车转向架定距为im,因此线段lab恒等于im;根据两点距离公式,可列出动点a和动点b之间的距离方程式:
Figure RE-FDA0003283666800000013
其中,动点b(xb,yb)代表列车前转向架中心点,动点a(xa,ya)代表列车后转向架中心点;xa为动点a的x座标值;ya为动点a的y座标值;xb为动点b的x座标值;yb动点b的y座标值;R表示目标圆曲线半径;L表示缓和曲线的总长度;x1为圆曲线圆心的横坐标;y1为圆曲线圆心的纵坐标;
已知曲线函数和列车后转向架中心点a的坐标后,可以通过距离方程式确定列车前转向架中心点b点的坐标值以及动点a、b之间的函数关系,在求出列车前后转向架中心点b、a的坐标后,可列出车体函数Lab的表达式:
Figure RE-FDA0003283666800000021
通过车体函数Lab斜率,确定辅助线Lcd的斜率;从车体函数Lab的表达式中可知斜率
Figure RE-FDA0003283666800000022
所构造的辅助线Lcd与车体函数Lab恒垂直,因此,可得Lcd的斜率:
Figure RE-FDA0003283666800000023
求lcd的最大值:求取最大值的过程关键在于确定动点d的坐标,d点为落在缓和曲线或圆曲线上的一个动点,对缓和曲线或圆曲线求导找出满足y′d=Kab的动点d,满足条件的动点d即是lcd取到最大值的位置;点d(xd,yd)表示运动在缓和曲线上或圆曲线上的动点,c(xc,yc)表示车体因几何偏移产生的最大偏移点;
当Kab<y′圆缓点,说明车体内侧最大偏移点未过圆缓点,处于缓和曲线上,因此,求解动点d时应对缓和曲线进行求导;
当Kab=y′圆缓点,说明车体内侧最大偏移点落在圆缓YH点上,缓和曲线和圆曲线在圆缓点处曲率几乎相等,因此,带入任意一种函数都可以进行计算;
当Kab>y′圆缓点,说明车体内侧最大偏移点越过圆缓点处于圆曲线上,因此,求解动点d时应对圆曲线进行求导,即:
Figure RE-FDA0003283666800000031
式中,y′圆曲线为圆曲线的导数;y′缓和曲线为缓和曲线的导数;y′圆缓点为圆缓点的导数值;Kab为车体函数的斜率;R表示目标圆曲线半径;L表示缓和曲线的总长度;
3.根据权利要求1所述的一种地铁圆缓段加宽的计算方法,其特征在于,所述列车在缓和曲线段因几何偏移造成的内侧曲线加宽值求解如下:
列车处于圆缓段时,前转向架中心点b点已经落在了圆曲线上,因此在圆缓段的超高应使用圆曲线的超高:
Figure RE-FDA0003283666800000032
式中,Hyh为计算圆缓段曲线加宽用到的曲线超高值,单位mm;Hmax为目标圆曲线的超高值,单位mm,V为过圆曲线时列车的设计速度,单位m/s;R为目标圆曲线半径,单位m;h代表限界控制点自轨面的高度,单位mm;1500为钢轨中心距,单位mm;Wyhv为圆缓段处因线路超高产生的加宽值,单位mm;
4.根据权利要求1所述的一种地铁圆缓段加宽的计算方法,其特征在于,所述列车在缓和曲线段因几何偏移造成的内侧曲线加宽值求解如下:
Figure RE-FDA0003283666800000041
其中,Wi为圆缓段处,车辆总的内侧加宽值,单位mm;
5.根据权利要求1所述的一种地铁圆缓段加宽的计算方法,其特征在于,所述列车在缓和曲线段的外侧曲线加宽值求解如下:
对列车过圆缓段求解外侧加宽值的方法,定义g(xg,yg)、e(xe,ye)、h(xh,yh)、f(xf,yf)四点;
设点a(xa,ya)为列车后转向架中心点;b(xb,yb)为列车前转向架中心点;地铁A型车总车长lef=Im、转向架定距lab=im;可知
Figure RE-FDA0003283666800000042
设动点e(xe,ye)、f(xf,yf);
求列车车尾外侧曲线加宽:
Figure RE-FDA0003283666800000043
可列长度方程:
Figure RE-FDA0003283666800000044
lea代表车尾到后转架中心点的距离;
在已知a点坐标的情况下,很容易算出e点坐标(xe,ye),再做垂直与Lef的直线Leg;Leg
Figure RE-FDA0003283666800000045
联立直线Leg和缓和曲线方程即可求出g点坐标xg,yg
计算出来的
Figure RE-FDA0003283666800000051
就是列车车尾外侧曲线加宽;
求列车车头外侧曲线加宽:
Figure RE-FDA0003283666800000052
可列长度方程
Figure RE-FDA0003283666800000053
在已知b点坐标的情况下,很容易算出f点坐标xf,yf,做垂直与Lef的直线Lhf;代表列车车头外侧加宽值:
Figure RE-FDA0003283666800000054
联立Lhf和圆曲线方程即可求出h点坐标xh,yh
计算出来的
Figure RE-FDA0003283666800000055
就是列车车尾外侧曲线加宽;
Figure RE-FDA0003283666800000056
Figure RE-FDA0003283666800000057
式中,W车头为列车行驶在缓和曲线处车头的外侧曲线加宽值,单位mm;W车尾为列车行驶在缓和曲线处车尾的外侧曲线加宽值,单位mm;
6.根据权利要求1所述的一种地铁圆缓段加宽的计算方法,其特征在于,所述缓和曲线段列车车体函数内外侧加宽函数曲线求解如下:
确定列车后转向架xa的取值范围,设后转向架中心点为点a(xa,ya),当x圆缓点≤xa时列车就完全进入了圆曲线,因此,x圆缓点-i≤xa≤x圆缓点
对xa在x圆缓点-i≤xa≤x圆缓点平均取点,以这些点为基准,求解列车相对应的内侧最大偏移点c(xc,yc)的坐标;
对xa在x圆缓点-i≤xa≤x圆缓点平均取点,以这些点为基准,求出相对应的车尾最大偏移点e(xe,ye)、车头最大偏移点f(xf,yf)的坐标;
用最小二乘法对步骤(1)、(2)、(3)(4)中所求的点进行函数拟合,最终求出在缓和曲线段车体函数内外侧加宽函数曲线。
7.根据权利要求2所述的一种地铁圆缓段加宽的计算方法,其特征在于,所述Kab>y′圆缓点时,只需把缓和曲线函数换成圆曲线函数重新求导计算即可:
Figure RE-FDA0003283666800000061
可得xd用a、b坐标表示的解析式:
Figure RE-FDA0003283666800000062
xd取正值,可用(xa,ya)、(xb,yb)来表示d点坐标:
Figure RE-FDA0003283666800000063
算得动点d的坐标和Lcd的斜率Kcd后可求出线段Lcd的表达式:
Figure RE-FDA0003283666800000071
接着,求解车辆内侧偏移点c点的坐标,联立Lcd与Lab的函数表达式,求解方程式:
Figure RE-FDA0003283666800000072
得:
Figure RE-FDA0003283666800000073
又因内侧偏移点c点处在车体函数Lab上,把xc带入Lab的方程式即可算得:
Figure RE-FDA0003283666800000074
最后对c点和d点坐标使用一次两点距离方程即可求出Lcd的值:
Figure RE-FDA0003283666800000075
式中,lcd为车体几何偏移产生的内侧曲线加宽值,单位mm。
CN202110805534.1A 2021-07-16 2021-07-16 一种地铁圆缓段加宽的计算方法 Active CN113722812B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110805534.1A CN113722812B (zh) 2021-07-16 2021-07-16 一种地铁圆缓段加宽的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110805534.1A CN113722812B (zh) 2021-07-16 2021-07-16 一种地铁圆缓段加宽的计算方法

Publications (2)

Publication Number Publication Date
CN113722812A true CN113722812A (zh) 2021-11-30
CN113722812B CN113722812B (zh) 2023-12-08

Family

ID=78673480

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110805534.1A Active CN113722812B (zh) 2021-07-16 2021-07-16 一种地铁圆缓段加宽的计算方法

Country Status (1)

Country Link
CN (1) CN113722812B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100474261B1 (ko) * 2004-02-05 2005-03-14 주식회사삼안건설기술공사 컴퓨터를 이용한 지하철 및 철도 선형설계 시스템 및설계방법
CN102161342A (zh) * 2010-12-08 2011-08-24 中铁四局集团电气化工程有限公司 地铁隧道设备限界检测装置
CN106126834A (zh) * 2016-06-28 2016-11-16 西安交通大学 一种基于惯导角度量测的轨道平面最佳线形确定算法
CN107024196A (zh) * 2017-05-11 2017-08-08 中南大学 一种高速铁路缓和曲线轨道横向偏差检测方法
CN107832518A (zh) * 2017-11-02 2018-03-23 北京交通大学 一种基于运动病诱发率的铁路线路设计方法
CN107901946A (zh) * 2017-12-12 2018-04-13 中铁第四勘察设计院集团有限公司 一种地铁列车限界检测系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100474261B1 (ko) * 2004-02-05 2005-03-14 주식회사삼안건설기술공사 컴퓨터를 이용한 지하철 및 철도 선형설계 시스템 및설계방법
CN102161342A (zh) * 2010-12-08 2011-08-24 中铁四局集团电气化工程有限公司 地铁隧道设备限界检测装置
CN106126834A (zh) * 2016-06-28 2016-11-16 西安交通大学 一种基于惯导角度量测的轨道平面最佳线形确定算法
CN107024196A (zh) * 2017-05-11 2017-08-08 中南大学 一种高速铁路缓和曲线轨道横向偏差检测方法
CN107832518A (zh) * 2017-11-02 2018-03-23 北京交通大学 一种基于运动病诱发率的铁路线路设计方法
CN107901946A (zh) * 2017-12-12 2018-04-13 中铁第四勘察设计院集团有限公司 一种地铁列车限界检测系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
KE WU 等: "Construction Mechanical Mechanism of Shallow Subway Tunnel Group with Large-Span Variable Cross Section", SPRINGERLINK *
林凤涛 等: "多边形磨耗对轮轨力和轴箱加速度的影响研究", 机电工程 *

Also Published As

Publication number Publication date
CN113722812B (zh) 2023-12-08

Similar Documents

Publication Publication Date Title
CN111442736B (zh) 一种基于激光扫描仪的铁路隧道变形检测方法及其装置
CN104573397B (zh) 一种计算铁道车辆动态包络线的方法
CN109214554A (zh) 一种不落轮工况下的铁道车辆轮对整车镟修策略决策算法
CN107685748B (zh) 基于激光位移传感器的列车车轮尺寸参数在线检测方法
CN106951657A (zh) 一种磨损钢轨打磨目标廓形快速设计方法
CN110414113B (zh) 自动生成铁路站场平面总体布置图的方法
CN104878667A (zh) 道岔转辙区的钢轨打磨方法
CN108763806B (zh) 高速列车过特长双线隧道时列车及隧道表面所受静压的数值模拟分析方法
CN203037249U (zh) 一种高铁隧道断面积检测系统
CN107024196B (zh) 一种高速铁路缓和曲线轨道横向偏差检测方法
CN114332212A (zh) 基于车载移动激光点云的轨道超高和前后高低检测方法
CN111877108A (zh) 一种基于线激光的路面平整度测量方法
CN204944427U (zh) 一种接触网几何参数动态检测小车
CN113722812A (zh) 一种地铁圆缓段加宽的计算方法
CN108020162B (zh) 基于二维激光扫描与三角原理的轨距仪及其使用方法
CN110580364A (zh) 多模块铰接式低地板城轨车辆曲线几何偏移量的计算方法
CN113010969A (zh) 一种高速列车进隧道压力波最大值的计算校正方法
CN101906758A (zh) 高速铁路中声屏障脉动力大小的确定方法及应用
CN115495832A (zh) 地铁线路直缓段处的加宽曲线计算方法
CN115495831A (zh) 地铁线路缓和曲线段的加宽曲线计算方法
CN108896000A (zh) 钢轨廓形对齐方法及装置、钢轨磨耗计算方法及装置
CN209162532U (zh) 一种时速120km/h地铁车辆段试车线
CN110345930B (zh) 一种基于改进的布尔莎模型的不整平自由设站方法
Zheng et al. Evaluation on braking stability of autonomous vehicles running along curved sections based on asphalt pavement adhesion properties
CN113722811B (zh) 一种列车明线交会时压力波幅值与车速比关系的估算方法

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
GR01 Patent grant
GR01 Patent grant