列车动力学的非线性数值分析方法.pdfVIP

、 ＼ 维普资讯 ＼ 总 第 均 期 西 南 交 通 大 学 学 报 1991年第l期 SUM 抽 ．79 J0RNALOFSOUTH~EST~[AOTONOUNIVERSITY 列车动力学的非线性数值分析方法 翟婉明 黄志辉 (机车车辆研究所) 【摘 要】 列车动力学方程是个复杂非线性方程 ，必须采用数值积分方法求解。本 文从工程实用的观点出发，提 供 7： 妾 简单有效的数值积分方法—— Newmark显式积分法及其显一臆式预测一校正积分法。这类积分方法 积 分过程 简捷 ，程序鳊制容易，由于积分过程 中无须求解高阶线性代数 方 程组，这就显著地提高7积分速度，并有效地节约7计算 内存 ，从而在 普通微机上实现 7对太型列车动力学问题的模拟分析。 文中建议 7一种判定工程实际问题之非线性数值积分稳定性 的 实 用 办 法 ，并成功地应 用于重裁列车动力学分析。在此基础上进一步实现 7变 步长积分计算，太大提高了计算速度。通过算例校验厦与试验结果的对 照 ，证实 7方法的非线性数值分析有效性。 【关键词】 列车动力学，非线性分析，数值积分，数值试验，数值稳定性 引 言 重载运输将成为我国铁路运输发展的一种重要模式。为使重载列车能够安全可靠运行， 开展列车动力学的研究意义十分重大 … 。 列车动力学包括列车纵向、垂向和横向动力学三个方面，主要研究列车在运行过程中的 纵向动力作用 以及由此而产生的列车垂向与横向动力作用。可见 ，列车纵向动力学是列车动 力学的核心。 ． 两车动力学方程可 以表述成如下统一形式： [M]{A}+[C]{y}+[K]{X}={P} (1) ● ●● 式中，{x}、{}(或{x})、 {A}(或 {x})分别是广义位移、速度和加速度矢量，CM]、[c]、 [K)依次代表列车动力系统 的质量、阻尼与刚度矩阵，{P}为广义载荷矢量 。 列车动力学方程 (1)是个复杂非线性方程 ，包含有许多非线性因素。其中比较典型的 如：车钩间隙，缓冲器的刚度}机车车辆的阻力，机车牵引特性 以及车辆制动特性等。 目前 ，求解非线性列车动力学方程只能采用直接数值积分方法 。较为常用的有： 四 阶Runge—Kutta法，Newmark一口法 ，Wilson-0法等。然而 ，这些方法都是在线性条件下 本 文于 1989卑 1呵7日收 到 维普资讯 获得的，用于非线性分析迄今为止在数学力学界尚无完全成熟理论。习惯的做法是，将 线 性分析的结果外推到一些非线性 问题中去。这显然存在一定弊端。倒如 ，对于线性问题是稳 定的方法 ，用于非线性分析时却不一定再保持其原有的数值稳定性、 。此外，这些方法用 于分析重载列车动力学这样的大型非线性工程问题 ，由于计算内存大、许用积分步长小，必 须要求大容量高速计算机的支持 ，这将付出昂贵的代 价。 例如 ，Gogineni和 Gsrg 曾在 IBM一370／158中型机上采用四阶Runge-Kutta法对 6O辆 车进 行 过运行模拟，最大积分步 长仅能取 1O～ s，此时模拟列车运行 1s竟需要 i00s机时，而当步长缩小为 8×10 s时 ， 这一数值便增