基于Abaqus轮轨耦合有限元法和固定界面模态综合法对比下的车-桥耦合振动分析.docxVIP

基于Abaqus轮轨耦合有限元和固定界面模态综合法

对比下的车-桥耦合振动分析

万鹏

江西省交通运输厅

摘要:在传统的车桥耦合仿真分析中，多为采用有限元软件求解的方法，将轮轨力加载在桥

梁之上，在有限元软件中直接积分。他的缺点在于计算时间长，效率低。近来，一种基于子

结构或者模态综合法的求解方法，通过缩减自由度提高了计算效率。本文对这两种方法进行

了对比。

关键词:ABAQUS联合仿真固定界面模态综合法车桥耦合振动仿真分析

1.引言

近年来高速列车行驶引起的桥梁振动安全问题受到了社会的高度重视，成为当今铁路交通研

究领域的一个热点问题。研究车辆-桥梁耦合振动需分析车辆荷载激发的振动波在桥面-桥梁-桥墩

-桥梁复杂系统中的传播。由于有限元方法有较好的适用性，因此在桥梁振动理论研究中得到了

广泛的应用。

然而，当前采用有限元方法模拟车辆-桥梁相互作用体系的动力问题还存在一些困难。当前

学者一般采用二维、三维车辆-桥梁直接耦合的有限元整体模型求解。有限元整体分析法可有效

求解复杂的车辆-桥梁耦合模型以及非线性问题如地震作用、温度作用、风荷载作用等，但其也

有明显缺点。在传统有限元的车辆-桥梁动力学模型中，网格划分尤其是轮胎与桥梁接触的部分

需要划分得很细、模型离散区域需很大，这直接引起有限元模型的单元数量和自由度非常庞大，

多达十几万甚至几十万，这大大的增加了计算工作量，也由于研究者因各种参数设置不合理导致

最终仿真计算结果不准确。

除有限元整体分析法外，另一种分析车辆-桥梁耦合振动的方法为固定界面下的模态综合法。

固定界面下的模态综合法将整体结构划分为两个子结构，再对各子结构通过坐标变换进行自由度

缩减，综合缩减后的各个子结构以形成整体系统的运动方程，最终对降阶的运动方程进行求解，

进而获得原结构体系的动态响应[1]。

各类动态子结构模态综合法中，固定界面模态综合法（本文采用UM软件联合仿真来实现）

由于概念清晰，易于理解，不会出现悬浮子结构[2]，常被使用。UM联合仿真法由Hurty[3]于1964

年正式提出，Craig和Bampton[4]于1968年改进了此方法。

UM联合仿真法基于固定界面模态综合法，常用于外部荷载作用下刚体-结构动力相互作用

分析。本文采用固定界面模态综合法分析车辆引起的桥梁振动问题，并与整体有限元法（简称

FEM法）进行对比分析，探讨了联合仿真法在车辆-桥梁振动分析中的适用性。

2.固定界面模态综合法基本理论及联合仿真关键流程

2.1固定界面模态综合法基本理论

（1）将整个结构分割成若干个子结构，子结构之间以及子结构与剩余结构（非子结构）的

连接处定义为交界面。每个子结构的节点被分为交界面节点b和内部节点i。

（2）固定交界面节点的自由度，由特征方程（1）式求解子结构固定界面主模态?k

K?（1）

ii?Mii

??

式中，k为保留的低阶模态数，Kii为分块刚度矩阵，Mii为分块质量矩阵。

（3）依次给交界面上的节点自由度以单位位移，同时保持其余自由度固定，由式(2)求

出子结构静力反应的约束模态?c

????

1

cKK

iiib

（2）

（4）子结构位移的物理坐标可由固定界面主模态和约束模态表示成模态坐标：

???????q?

u

??Tq

ikck?

u??

??????

u0Iq

??????

bbbb

（3）

式中T为坐标变换矩阵，q为模态坐标。

(5)子结构的在物理坐标下的运动方程为：

M???????t?（4）

qCqKqF

T

T

将式（3）代入式（4），并在方程左侧乘以

，则模态坐标下的子结构运动方程表示为：

M??????（5）

qCqKqF

?t?

式中，M、C、K、F分别为缩减的广义质量、阻尼、刚度和广义荷载向量。

（6）结合子结构交界面位移协调条件，将子结构广义矩阵组装到非子结构上，得到整体结

构的缩减矩阵方程。求得广义坐标下的位移解，按照式（3）即可最终返回体系在物理坐标下的

位移时程响应u?t?。

2.2UM联合仿真技术流程

UM联合仿真技术流程主要有三个步骤：有限元部分、多体系统部分和联合组装部分。

（1）有限元部分主要是通过在Abaqus/CAE软件中建立好需要导入至UM软件的部分，

并提取主自由度将其进行自由度凝聚，提取需导入的模型的动力特性信息，并生成FSS后

缀文件。

（2）多体系统部分是在UM软件中完成的。通过UM软件建立好车辆系统的模型。由

于UM软件的建模可组装性强，模型可套性强，建立多车系统更为简便。

（3）联合组装部分。通过坐标对齐和模型正则化校正，将有