下载文档原格式
轨道交通列车运行模型与仿真轨道交通是现代化城市交通系统的重要组成部分，它具有高效、安全、环保等特点。
而为了确保轨道交通系统的稳定运行，需要对列车的运行模型进行研究和仿真。
本文将探讨轨道交通列车运行模型与仿真的相关内容。
1. 列车运行模型的概念列车运行模型是指描述列车在轨道上运行的数学模型，它能够模拟列车在不同速度、负载情况下的运行状态和特性。
列车运行模型通常包括列车的机械动力学模型、车辆动力学模型以及轨道特性等。
2. 列车的机械动力学模型列车的机械动力学模型主要研究列车的运动规律和力学特性。
它考虑列车的质量、摩擦力、空气阻力等因素，通过数学公式描述列车的加速度、速度和位移等变化。
机械动力学模型能够帮助我们分析列车的启动、制动和驶入弯道等过程中的力学特性，从而为优化列车运行提供理论依据。
3. 车辆动力学模型车辆动力学模型研究的是列车上各个车辆之间的相互作用和运动特性。
不同类型的车辆在轨道上行驶时，会产生相互作用力，如车辆之间的接触力、车辆与轨道之间的相对位移等。
车辆动力学模型能够模拟不同车辆之间的相互作用，从而通过仿真分析提升列车的运行效能和安全性。
4. 轨道特性模型轨道特性模型考虑轨道的几何形状、轨道弯曲半径、轨道轮廓等因素。
它能够模拟列车在弯道、坡道、过衔接段等复杂轨道条件下的运行。
轨道特性模型能够通过仿真预测列车在特定轨道条件下的运行情况，为轨道交通系统的规划和设计提供可靠依据。
5. 列车运行仿真的意义列车运行仿真是研究轨道交通系统的重要手段之一。
通过运行仿真，可以模拟不同负载条件下列车的运行状态，分析列车的稳定性和安全性，优化列车的运行参数，提升列车的运行效能。
此外，列车运行仿真还可以用于对轨道交通系统的运行策略进行优化，减少拥堵，提高运输能力。
6. 列车运行仿真的方法列车运行仿真可以采用多种方法进行。
一种常用的方法是基于物理模型的仿真。
这种方法通过建立列车运行模型，利用计算机技术模拟列车在不同条件下的运行情况。
车站列车运行数量时间序列模型分析通过给定的数据，判断数据的平稳性和季节性，采取相应的措施，消除这些影响，然后利用B-J方法和P-W方法对所选的数据进行分析，进行模型的识别，阶数的识别，参数的拟合和适应性检验，最后再利用所留数据进行预测。
根据预测结果对求得的模型进行判断。
本文中已知的是某车站1993-1997年各月的列车运行数量60个，单位：千列/千米，具体数据见附录。
一、数据的平稳性判断利用matlab画出数据的二维图，如下图（1），程序（1）见附录图（1）由上图我们可以看到数据是不平稳的，因此需要对数据进行平稳化处理，利用一阶差分的方法得到52个数据，再次利用matlab画图，见图（2），程序（2）见附录图（2）有图（2）我们可以初步判断xt序列是平稳的,再对稳化后的数据进行零均值处理，得到的数据见附录。
二、利用B-J方法建模（1）模型的识别如果一个时间序列是由某一类模型生成，理论上它就应该具有相应的自相关特征，因而我们可以计算出平稳时间序列的样本自相关函数和样本偏自相关函数，将其特性与不同类型的序列的理论自相关函数和偏自相关函数的特性进行比较，进而初步判断序列x t所适合的模型类型。
在求样本的自相关函数和偏自相关函数时，利用了Eviews，得到如下数据下面利用|k i ρ+|≤2121[(12)/N]ml l ρ=+∑判断样本自相关函数的截尾性取]=7，当m=1时，1122221(12)(12*0.679)/0.192ρ+=+=在(k 2,3,4,5,6,7,8)k ρ=中满足|k ρ|<0.192的占6/7=85.7%>68.3%,因此k ρ一步截尾，可初步判断差分后的序列适合MA(1)模型。
（2） 阶数的识别图（4）由图（4）可以看出MA(1)模型的剩余平方和为73855.52图（5）由图（5）可以看出MA(2)模型的剩余平方和为73187.38图（6）由图（6）可以看出MA(3)的剩余平方和为67346.46图（7）由图（7）可知，剩余平方和为：66468.54图（8）由以上的五个图形可以看出对于零均值化平稳后的序列拟合1-5阶的MA 模型剩余平方和：73855.52 73187.38 67346.46 66468.54 68804.69MA(5)的剩余平方和超过了MA （4）的剩余平方和，因此从MA （4）考虑模型阶数是否可以降低对于MA （3）和MA （4）模型=0.647取α=0.05查F 分布表可得，F(1，49)大约为4.03，显然F<F(1，49)所以在α=0.05的显著性水平下MA(3)和MA(4)没有显著性差异。
轨道交通列车运行动力学建模与仿真分析轨道交通是现代城市中不可或缺的一种交通方式，它具有快速、高效和环保等特点。
为了确保轨道交通的安全和稳定运行，轨道交通列车的运行动力学建模和仿真分析显得至关重要。
本文将探讨轨道交通列车运行动力学建模与仿真分析的方法和应用。
一、轨道交通列车运行动力学建模轨道交通列车运行动力学建模是指根据列车的运行规律和机械特性，建立数学模型来描述列车在运行过程中所受到的各种力和运动状态变化。
常见的列车运行动力学模型包括单物体模型和多体模型。
1. 单物体模型单物体模型假设轨道交通列车为一个整体，通过运动学原理和力学方程来描述列车的运动状态。
该模型适用于研究列车的加速度、速度、位移和运动平稳性等基本动力学特性。
2. 多体模型多体模型考虑列车车体、燃油车或电动机、车轮、轨道等多个物体之间的相互作用。
它通过建立列车系统的运动方程，包括车体的平动方程和转动方程，来揭示列车系统的运行机理和特性。
多体模型能够更精确地描述列车的运行过程，适用于研究列车的动力学响应、悬挂系统的特性和车辆稳定性等问题。
二、轨道交通列车运行动力学仿真分析轨道交通列车运行动力学仿真分析是指利用计算机软件模拟列车的运行过程，通过运行结果的模拟和分析，评估列车的运行性能和安全性。
1. 建立仿真模型首先，根据实际列车的参数和运行规律，建立列车的数学模型。
这一步可以采用前述的单物体模型或多体模型，也可以根据实际需求进行模型的适当简化。
2. 数据采集和验证在进行仿真分析之前，需要收集实际运行数据，以验证所建立的模型的准确性。
这些数据可以包括列车的加速度、速度、位移以及与之相关的温度、摩擦等指标。
3. 进行仿真运行根据实际运行情况和仿真模型，输入相应的控制指令和参数，进行仿真运行。
运行过程中，可以观察列车的运动状态变化和各种力的作用情况。
4. 仿真结果分析通过仿真结果的分析，可以评估列车的运行性能和稳定性，并确定是否需要进行相应的优化调整。
第二章列车运行的基本数学模型空气波传播的定性分析在重载列车长大下坡周期性制动策略研究中，首先需要解决的难点就是长大下坡速度控制问题。
由于重力的作用，列车会在长大下坡不断加速，为了保证列车的安全行驶必须施以足够的制动力。
在此过程中，如果一直使用空气制动，长时间制动将使闸瓦过热，其次，由于制动缸漏泄制动力会不断衰减。
所以重载列车长大下坡时必须间断使用空气制动，并且掌握好制动与缓解时机。
由于列车管系中气体是不规则变化的，即无法直接预测整个列车管系中各部分压强的压强变化，在此将通过气体流动理论计算制动特性来可避免各种假定压强带来的计算不正确。
当速度过高时采取制动措施可能会出现运行速度超过限速；速度过低开始制动一方面可能运行速度低，线路通过能力没有完全发挥，另一方面可能出现因为缓解时间过短副风缸没有充风完毕，会出现制动力不足问题。
而缓解时，如果速度过高开始缓解，则可能出现需要制动时副风缸还没有充满，此时制动，制动缸将没有足够的制动力，列车可能会出现失控的危险；如果速度过低开始缓解，则通过线路的速度过低。
2.1.1空气制动系统描述首先对空气制动系统的机械结构和控制过程中的三个过程的认知理解，我们将列车制动机械结构简化抽象，在此基础上对降压制动、充风缓解和控制保压三个过程进行研究。
图1 空气制动系统示意图2.1.1.1制动单元车辆制动单元包括以上多个机械和空气组成部件,由于车辆类型和用途的不同,其具体的组成部件也会有所区别,但是大多数车辆制动单元的基本部件是一致的,在此我们主要简化研究控制阀（120三通阀）、副风缸和制动缸之间气压变化机理。
图2 空气制动单元示意图2.1.1.2制动位列车运行中准备进站停车或者减速时,通常是施行常用制动。
司机施行常用制动减压后,因副风缸压力空气来不及通过滑阀和滑阀座向列车管逆流,于是在主活塞两侧产生了一定的压力差,此压力差产生的向上作用力克服了节制阀与滑阀背面间的摩擦阻力、橡胶膜板的变形阻力和压缩稳定弹簧的阻力以及主活塞重力等向下作用力的总和,使主活塞先带动节制阀上移,然后再带动滑阀上移,此时,由于120阀的动作,阻止了制动管的空气流向副风缸,同时紧急阀上的放风阀也被阻塞。
高速列车运行的数学模型研究随着科技的发展和人们生活水平的提高，高速列车作为高效、便捷的交通工具，已成为人们日常出行和商务活动的首选。
但是高速列车的稳定性和安全性是运营中最为关注的问题之一。
为了保证高速列车的安全性和稳定性，需要建立一套科学的数学模型。
本文将从高速列车的运行原理、数学模型的建立原则和具体模型构建三个方面进行研究。
一、高速列车的运行原理高速列车的运行原理是利用电机发动机提供动力，使列车在铁路轨道上运行。
其中，列车在轨道上运行的方式包括摩擦运动和滚动运动两种。
通过摩擦运动，车轮与轨道之间的摩擦力将列车推动前进；通过滚动运动，车轮在轨道上滚动，实现列车的运行。
此外，高速列车在运行中还会面临一系列问题，如平衡、防止车轮脱轨、减震、防止飞蹬等。
因此，需要建立一套完整的数学模型，对高速列车的运行过程进行模拟和分析，以进一步提高其安全性和稳定性。
二、数学模型的建立原则建立高速列车的数学模型需要遵循一定的原则，即：1.准确性原则：数学模型应能准确地表达高速列车的运动规律。
2.可验证性原则：数学模型应能被实验所验证。
3.实用性原则：数学模型应便于应用和推广。
4.合理性原则：数学模型应依据合理的物理原理和理论基础。
基于以上原则，我们可以建立高速列车的运行模型，并对其进行分析和验证，以确保模型的准确性和实用性。
三、具体模型构建1.高速列车的速度运动模型高速列车的运行速度是其最为关键的参数之一。
因此，我们需要建立一套速度运动模型，以模拟高速列车在运行中的速度变化规律。
其中，高速列车的速度运动模型可以采用泊松分布模型，其表达式为：V(t) = V0 - ΔV(e^(-λt)-1)其中，V(t)表示高速列车在任意时刻t的速度值；V0为高速列车的初始速度；ΔV为高速列车在经过一个区间的速度变化值；λ为泊松分布的参数。
2.高速列车的转弯模型高速列车在转弯时需要保持平衡和稳定。
因此，我们需要建立一套转弯模型，以模拟高速列车在转弯时的运动规律。
第二章列车运行的基本数学模型2.1 空气波传播的定性分析在重载列车长大下坡周期性制动策略研究中，首先需要解决的难点就是长大下坡速度控制问题。
由于重力的作用，列车会在长大下坡不断加速，为了保证列车的安全行驶必须施以足够的制动力。
在此过程中，如果一直使用空气制动，长时间制动将使闸瓦过热，其次，由于制动缸漏泄制动力会不断衰减。
所以重载列车长大下坡时必须间断使用空气制动，并且掌握好制动与缓解时机。
由于列车管系中气体是不规则变化的，即无法直接预测整个列车管系中各部分压强的压强变化，在此将通过气体流动理论计算制动特性来可避免各种假定压强带来的计算不正确。
当速度过高时采取制动措施可能会出现运行速度超过限速；速度过低开始制动一方面可能运行速度低，线路通过能力没有完全发挥，另一方面可能出现因为缓解时间过短副风缸没有充风完毕，会出现制动力不足问题。
而缓解时，如果速度过高开始缓解，则可能出现需要制动时副风缸还没有充满，此时制动，制动缸将没有足够的制动力，列车可能会出现失控的危险；如果速度过低开始缓解，则通过线路的速度过低。
2.1.1空气制动系统描述首先对空气制动系统的机械结构和控制过程中的三个过程的认知理解，我们将列车制动机械结构简化抽象，在此基础上对降压制动、充风缓解和控制保压三个过程进行研究。
图1 空气制动系统示意图2.1.1.1制动单元车辆制动单元包括以上多个机械和空气组成部件,由于车辆类型和用途的不同,其具体的组成部件也会有所区别,但是大多数车辆制动单元的基本部件是一致的,在此我们主要简化研究控制阀（120三通阀）、副风缸和制动缸之间气压变化机理。
图2 空气制动单元示意图2.1.1.2制动位列车运行中准备进站停车或者减速时,通常是施行常用制动。
司机施行常用制动减压后,因副风缸压力空气来不及通过滑阀和滑阀座向列车管逆流,于是在主活塞两侧产生了一定的压力差,此压力差产生的向上作用力克服了节制阀与滑阀背面间的摩擦阻力、橡胶膜板的变形阻力和压缩稳定弹簧的阻力以及主活塞重力等向下作用力的总和,使主活塞先带动节制阀上移,然后再带动滑阀上移,此时,由于120阀的动作,阻止了制动管的空气流向副风缸,同时紧急阀上的放风阀也被阻塞。
空气只能从副风缸流向制动社,推动活塞向左边运动,然后通过基础制动装置作用到车轮上,如图3所示。
（该过程在常用制动和紧急制动时都会发生,它们的区别就在于常用制动时,列车管的空气即不向车辆制动单元流动也不排向大气,但是紧急制动时,120阀除了有上述动作,还会通过紧急放风阀将列车管中的空气直接排向大气。
换言之,紧急制动过程中,120阀通过紧急放风阀的动作,会大大加快列车管的放风动作。
）图3 制动位示意图2.1.1.3缓解位列车管经过减压,使列车进行制动后,当再次向列车管充气时,由于列车管压力持续增加,破坏了主活塞两侧在制动保压的压力平衡状态,当主活塞两侧的压力差产生的向下作用力与主活塞重力之和,超过了滑阀与滑阀座间的摩擦阻力时,主活塞便带着滑阀一起下移至充气及缓解位。
此时,副风缸得到充气,副风缸持续充气直到压力跟列车管压力一致。
一旦副风缸充满了,就可以进行第二次制动操作了。
滑阀的这次动作同时将副风缸和制动缸的通路阻断,而且制动缸排气孔打幵向外排气。
如图4所示。
此时,制动缸内被压缩的缓解弹黃要恢复原状,其复原力将活塞推回,闹瓦也就与车轮分离了。
一旦整列车的副风缸完成充气,列车管停止供气。
图4 缓解位示意图2.1.1.4保压位实施了常用制动后,当压力表显示达到所要求的列车管减压量时,将司机制动手把移到保压位,使列车管停止继续减压,制动保压状态被激活,此时滑阀处于中间位置。
因为副风缸持续向制动缸供气,副风缸的压力逐渐减小直到跟列车管压力一致,此时滑阀移到保压位如图5所示。
三通阀在该位置时,切断了所有列车管、副风缸和制动缸之间的通路,而且本身还处于制动模式,制动的排气孔也是关闭的,因此,制动缸压力保持恒定,达到了一个持续的制动状态。
图5 保压位示意图2.2 长大列车空气管系二维充气特性模型及算法长大列车空气管系二维充气特性模型应用了现代计算流体动力学数值方法，从理论上论证了二维模型的合理性及其处理支管的有效性，即它没有使用类似一维模型引进经验公式帮助求解，并提出了一种改进的算子分裂算法求解压力速度耦合方程。
在该论文的结尾处，论文编写者将计算结果和国内外长大列车充气特性的有关试验数据进行对比分析，验证了所提出的理论和算法的正确性。
研究了列车编组辆数、管系组成、管系泄漏以及支管长度和直径大小等因素对列车管充气压力的影响。
2.2.1二维数学模型的建立一般说来，铁道列车空气制动系统可分为两个子系统，即列车制动系统和车辆制动系统。
其中列车制动系统包括安装在机车上的司机制动阀和沿列车长度方向的列车管系。
为能有效地数值模拟充气过程，假设每节车辆的列车管仅由一根主管和一根支管组成，只需要通过调整管内壁的摩擦系数来等效模拟弯管，软管和折角塞门对列车管内空气流动状态的影响。
列车管简化示意图如下图：图6 列车管简化示意图管内气体随时间变化十分剧烈，是典型的非定常流动。
在假定管壁是刚性的，主管无逆流，并不计空气重力的条件下，二维模型的控制方程可由以下状态方程，连续性方程和动量方程组成。
()()14= 02nL c T c k p U UU p F t f F u u ρρρ⎧⎪⋅⎪⎪∂+∇∙=-∇-⎨∂⎪⎪⎛⎫=⎪ ⎪⎝⎭⎩上式中，,∇∇∙分别表示梯度算子和散度算子，在本模型中可以表示为：,T x y x y ⎛⎫∂∂∂∂∇=∇∙=+ ⎪∂∂∂∂⎝⎭上述方程(),,,,L p t p U u v ρ=分别为密度、压力、时间、泄漏量和速度。
假定管内空气为满足相同状态方程的理想气体，则1=k RT,R T ，分别为气体常数和绝地温度；n 为多变过程指数；4=c f f μ⋅，为管壁摩擦系数与等效阻尼系数的乘积，,c f μ分别指管壁摩擦系数和弯管、软管、支管和折角塞门的等效阻尼。
2.2.2二维数学合理性论证在实际生活中，列车管的空气流动模型应当是一个三维模型，其轴向几何尺寸比径向尺寸大得多，可以假定在主管的横截面上的压力场和速度场相等，将二维模型看成三维方形横截面模型的纵截面模型。
首先所在圆柱型的控制容积上进行积分，再利用散度定理将部分体积分转化为面积分后可得：()+c S S U dV UU ndS p ndS F dV t ρρΩΩ∂⋅=-⋅-∂⎰⎰⎰⎰由于控制容积的长度相同，控制容积所导致的误差由横截面的面积差决定，示意图如下：图7 列车管横截面改变示意图D 表示主管直径，由两个横截面的周长相等可以得到()22D dy π-⋅=，计算可得两个横截面积之差为2414S D π-∆= 。
b D 表示支管的接口宽度，类似的我们可以将支管的圆形横截面等效成方形横截面，在保证宽度一致的条件下修正支管直径为：()0/2b b D D D D π=-⋅-在该论文中取定0.0375 ,D m =00.025 ,b D m ==0.0179 b D m 可以得到2.2.3基于算子分裂法的改进算法该论文所使用的改进算法与传统算子分裂法的不同之处是将压力与速度耦合的双曲型方程变换成关于压力的抛物型方程，保证了算法的稳定性和强健性。
改进的算法中假定了相变指数,01L n p ==泄漏量。
在每个时间步压力场和速度场的求解分两步进行，但不需要进行交替迭代。
第一步：假定密度和质量流已知，并不考虑压力算子，求解速度场如下：()()()()()()1114 02N N N N c T N N c U U U F t f F u u ρρ+++⎧∂⎪+∇∙=-⎪∂⎨⎪⎛⎫= ⎪⎪⎝⎭⎩ 第二步：求解新的压力和速度场，不考虑对流和速度场的源项： 考虑了泄露的连续性方程：()110N N p k U t ρ++∂+∇∙=∂动量方程： ()11N N U p t ρ++∂=-∇∂对动量方程向后欧拉差分离散得到：()()()11*1N N N U U p dt ρρ+++=-∇再将差分处理前后的两个方程进行代入整理可得：()()11*10N N N p k U p dt t ρ+++∂+∇∙-∆=∂计算过程中先使用此方程求解压力，在利用初始状态方程修正密度，最后使用差分处理之后的方程修正质量流和速度。
2.2.4基于有限体积法的离散方程求解方程上述方程时，首先需要对其进行离散，包括时间离散和空间离散，因此该论文选择有限体积法，因为有限体积法推广到任意非规则的区域时比其他几种方法更容易。
此处不再对过程进行赘述，仅给出有用的信息与结论。
列车管初始端的压力可由试验数据按如下的指数函数形式拟合得到：()(),11ct r i f i p p p p e -=+-- 本文利用初始端压力给定，列车管终端的速度为，0并在初始端压力控制容积上对连续性方程积分可得：()1100000N N N N P p p kV u S u S dt ρρ++--+= 由于计算能力和相应数据的欠缺，我们暂时无法做出相应的仿真数据和图像。
因此我们暂时先对论文结果进行分析。
2.3现有论文结果1、随着泄漏量的增加，列车管的增压速度减慢，而且沿列车管长度方向离入口处越远的车辆其增压速度减慢得越多。
图8 泄漏量不同时的充气特性曲线2、随主管长度增加，主管增压速度变慢，在同样充气时间下，充气压力随编组辆数的增加而减少。
图9 编组不同时的充气特性曲线3、管内壁的粗糙度对充气特性的影响较大，随管内壁摩擦系数的增大，主管的增压速度变慢。
图10 管壁粗糙度不同时的充气特性曲线4、越长或支管的直径越大，充气时间会有所延长，但它们的影响不是很大。
图11 支管长度和直径不同时的充气特性曲线2.2 牵引计算系统模型根据牵引计算采用的列车模型,一般分为两类:一类是单质点模型,另一类是多质点模型。
两者的区别在于,单质点把整列车看成一个没有尺寸、大小的质点,所有的受力都作用在一个点上，这样的好处在于受力分析比较简单,牵引计算容易实现。
《牵规》就是以单质点模型为编写依据的，而多质点模型把整列车看成一条有长度的质点链,由于有了长度,一所以每个机车和车厢的受力情况不是时刻处于同一线路段中,由于位置的不同,各个机车和车厢的受力情况也不同,所以进行受力分析就更复杂。
与单质点相比,多质点模型描述的牵引情况 更接近实际,在有计算机辅助的情况下,越来越被人们采用。
单质点模型传统牵引计算模型采用单质点列车模型,将列车简化为单个无尺寸的质点,考虑列车的受力变化时,将列车的受力变化放到质点上进行计算,而不考虑列车内部车辆与车辆之间的受力情况。
单质点的列车模型如图所示:列车单质点模型列车受到外力全部作用在质点上,主要包括:牵引力qy W 基本阻力jb W 由于曲线和坡道等产生的附加阻力ij W 制动力zd W 列车自身的重力mg 以及线路对列车的支持力N 等。
