列车运行的基本数学模型

空气波传播的定性分析

在重载列车长大下坡周期性制动策略研究中，首先需要解决的难点就是长大下坡速度控制问题。由于重力的作用，列车会在长大下坡不断加速，为了保证列车的安全行驶必须施以足够的制动力。在此过程中，如果一直使用空气制动，长时间制动将使闸瓦过热，其次，由于制动缸漏泄制动力会不断衰减。所以重载列车长大下坡时必须间断使用空气制动，并且掌握好制动与缓解时机。由于列车管系中气体是不规则变化的，即无法直接预测整个列车管系中各部分压强的压强变化，在此将通过气体流动理论计算制动特性来可避免各种假定压强带来的计算不正确。

当速度过高时采取制动措施可能会出现运行速度超过限速；速度过低开始制动一方面可能运行速度低，线路通过能力没有完全发挥，另一方面可能出现因为缓解时间过短副风缸没有充风完毕，会出现制动力不足问题。而缓解时，如果速度过高开始缓解，则可能出现需要制动时副风缸还没有充满，此时制动，制动缸将没有足够的制动力，列车可能会出现失控的危险；如果速度过低开始缓解，则通过线路的速度过低。

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.1二维数学模型的建立

一般说来，铁道列车空气制动系统可分为两个子系统，即列车制动系统和车辆制动系统。其中列车制动系统包括安装在机车上的司机制动阀和沿列车长度方向的列车管系。为能有效地数值模拟充气过程，假设每节车辆的列车管仅由一根主管和一根支管组成，只需要通过调整管内壁的摩擦系数来等效模拟弯管，软管和折角塞门对列车管内空气流动状态的影响。列车管简化示意图如下图：

图6 列车管简化示意图

管内气体随时间变化十分剧烈，是典型的非定常流动。在假定管壁是刚性的，主管无逆流，并不计空气重力的条件下，二维模型的控制方程可由以下状态方程，连续性方程和动量方程组成。 ()()14= 02n

L 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 ρρΩΩ∂⋅=-⋅-∂⎰⎰⎰⎰r r

由于控制容积的长度相同，控制容积所导致的误差由横截面的面积差决定，示意图如下：