基于流固耦合法对空气弹簧的有限元分析

孔繁晓 株洲市职工大学（工业学校）

基于流固耦合法对空气弹簧的有限元分析

孔繁晓 株洲市职工大学（工业学校）

流固耦合数值计算在目前许多工程领域需求越来越强烈，本文简要介绍了流固耦合的特点和目前可行的求解方法，采用Star-CD和Abaqus直接双向耦合的方法对空气弹簧的垂向动刚度做了数值计算，证明其方法的可行性且计算结果与试验较为吻合。

流固耦合；空气弹簧；垂向动刚度；Star-CD；Abaqus

1、概述

流固耦合是一种物理现象，如果流体与固体的相互作用很强烈（单独求解固体或者流体不能满足工程需要），那么我们有必要计算这样的物理现象。流固耦合力学是流体力学与固体力学交叉而形成的一门力学分支，它是研究变形固体在流场作用下的各种行为以及固体位形对流场的影响，这二者交互作用的一门科学。流固耦合力学的重要特征是两相介质之间的交互作用：变形固体在流体载荷作用下会产生变形或运动，而变形或运动又反过来影响流场，从而改变流体载荷的分布和大小。流固耦合问题从控制方程解法上可以直接求解的强耦合和分区迭代求解的弱耦合，所谓强耦合，是将流体域、固体域和耦合作用构造在同一控制方程中，在同一时间步内同时求解所有变量，商用软件ADINA、COMSOL等均是采用强耦合方法求解流固耦合问题；弱耦合是在每每一时间步内分别依次对CFD方程和CSD方程求解，通过中介交换固体域和流体域的计算结果数据，从而实现耦合求解，弱耦合是现在工程问题中较为普遍的方法，商用软件MPCCI作为中介交换数据的商用软件，可以实现固体软件与流体软件之间的数据交换，而Star-CD与Abaqus在最新的版本之间也实现了双向直接耦合，且无需中介软件做数据传递，本文采用Star-CD与Abaqus作为流固耦合求解器完成对空气弹簧动态垂向刚度工况的计算模拟。

空气弹簧是一种新型的弹性元件，它是在柔性的密闭气室内充入压缩空气，利用空气的可压缩性实现弹性作用的一种非金属弹簧，具有减振、缓冲等功能，由于其变刚度、低自振频率、高度可控以及良好的高频振动吸收和隔声性能，具有相当广泛的应用，主要领域有轨道车辆、汽车悬挂、动力机械等方面。

2、空气弹簧有限元模型的建立

2.1 计算模型

本文以某地铁转向架上所用的空气弹簧为计算对象，由于空气弹簧结构和垂向载荷的对称性，故而在本次计算中采用平面对称模型，采用CAD软件Pro/e建立1/2的空气弹簧三维模型，如图1所示，包括了上下盖板、橡胶气囊、空气，再通过专业有限元前处理软件Hypermesh对其进行网格划分，为了提高计算的速度与精度，在模型中全部采用六面体网格，图2所示网格模型， 其中包括了流体软件Star-CD所需要的网格和固体软件Abaqus需要的网格，分别把各自需要的网格倒入到求解器中。

图1 空气弹簧三维几何模型

图2 网格模型

2.2 材料参数

由于橡胶材料是体积模量远远大于其弹性模量的各向同性的不可压缩的超弹性材料，采用二参数的Mooney-Rivlin理论描述橡胶类不可压缩超弹性材料应变能：

式（1）中：Ii是第i方向的应变不变量；k是体积弹性模量；C1和C2是待定的Mooney-Rivlin参数，本文采用Abaqus有限元软件中的Rebar单元技术建立空气弹簧中的增强帘线。表1所示为Mooney-Rivlin具体参数值，而表2为帘子布力学特性。

表1 Mooney-Rivlin参数

表2 帘子布力学特性

3、 流固耦合计算结果与分析

垂向动态刚度是空气弹簧的重要动态力学性能参数，根据试验的要求，以1Hz的频率在垂向下压10mm位移，观察橡胶气囊的应力状态，图3（a）、（b）为0.1 s时刻和0.4s时刻气囊的应力云图，可以观察到在底部的应力受到较大的压强，出现了局部应力较大的情况。

图3 垂向动刚度各时刻应力云图

而图4（a）、（b）相对应时刻内的压力场，从云图上可以看出从开始时刻压力场较小，随着下压，囊内体积变小，压力会随之升高，符合实际的过程，再次验证了该方法的可行性。

图4 对应时刻的囊内压力场

不考虑橡胶材料本身的阻尼，得到空簧的动刚度曲线，如图5所示，随着位移量的增加，气囊的中的压力的变化对刚度产生了一定的振荡，本模拟是基于1Hz下的计算，可以预见随着加载频率改变的同时，对刚度的影响会随着频率的变化而不同。

图5 垂向动刚度曲线（1Hz）

4 总结

流固耦合数值计算在目前许多工程领域需求越来越强烈，本文简要介绍了流固耦合的特点和分类，主要分为强耦合和弱耦合，以及最终采用弱耦合方法中的Star-CD与Abaqus直接耦合的方法对空气弹簧进行流固耦合分析，分析了在垂向以1Hz频率下下压10mm的作用，橡胶的应力分布及内部气体压强的变化，在不考虑橡胶本身材料的阻尼情况下，得到空气弹簧气囊的垂向动刚度曲线，发现气体的流动对动刚度的影响较为明显。

[1]袁春元,周孔亢,吴琳琪，等 车用空气弹簧有限元分析方法[J].机械工程学报.2009，45（6）：262～266

[2]A.Homeyer，等.空气弹簧系统的现代计算方法[J]. 电力牵引快报.1999,4:21～26

[3]叶珍霞，朱海潮，鲁克明，等.囊式空气弹簧刚度特性的非线性有限元法研究[J].振动与冲击.2006,25（4）：94～97

10.3969/j.issn.1001-8972.2011.14.012