基于ABAQUS CEL方法的轮胎滑水仿真v.docx

基于ABAQUSCEL方法的轮胎滑水仿真

李兵,李炜,吴福麒

佳通轮胎（中国）研发中心

摘要:轮胎的滑水性能与车辆的行驶安全性密切相关。本文基于AbaqusCEL方法建立轮胎

滑水仿真技术，并使用若干组计及复杂胎面花纹的轮胎有限元模型，进行轮胎滑水案例分

析，考察路面积水深度、车辆行驶速度等工况参数对轮胎滑水性能的影响。结果表明，积

水深度和行驶速度对轮胎滑水性能影响显著，轮胎滑水性能随着水深的增加而降低，随着

行驶速度的增加而显著降低。还针对三种花纹样式进行轮胎滑水性能对比分析，并从设计原理上给出解释。仿真结果与实际情况或测试结果相符，这表明，使用ABAQUS软件进行轮胎滑水仿真能够得到合理的分析结果。

关键词:轮胎滑水有限元分析流固耦合ABAQUSCEL

1.引言

轮胎的滑水性能与车辆的行驶安全性密切相关。当以较快速度行驶的车辆驶过覆有较

深积水的路面时，因水的流体动力学作用使轮胎与地面之间形成水膜，导致轮胎不能充分

与地面接触甚至完全脱离地面，这种轮胎被水托起的现象称为轮胎的滑水（Hydroplaning）

现象；相应的轮胎使用性能被称为轮胎的滑水性能。轮胎研发工程师努力改进轮胎设计，

尤其是胎面花纹设计，以提高轮胎的滑水性能。然而长期以来，工程师们主要凭借经验进

行相关的优化设计，缺乏方便有效的评价手段。在试验方面，轮胎滑水性能测试对试验场

地、水深、车辆、试车手、试验工况、安全保障措施等多方面的试验条件要求苛刻；在仿

真方面，由于轮胎滑水现象涉及到复杂的力学问题，包括高度的非线性、复杂的流固耦合

问题等，长期缺乏有效可靠的仿真工具，也难以建立有效可靠的仿真方法。

2008年，Abaqus在其新版本中正式推出CEL（CouplingEuler-Lagrange，耦合的欧拉-

拉格朗日）方法，大幅拓展了Abaqus软件的流固耦合仿真能力。这也是Abaqus技术发展的

一个重要的里程碑。本文企图基于AbaqusCEL方法，建立一套完整的轮胎滑水仿真建模和

分析技术；与此同时，使用若干个计及复杂胎面花纹的轮胎有限元模型进行轮胎滑水案例

分析，考察路面积水深度、轮胎行驶速度、胎面花纹设计等因素对轮胎滑水性能的影响。

2.轮胎滑水仿真模型

轮胎滑水仿真模型由轮胎模型、水/空气复合模型、轮辋模型和路面模型四部分组成，

叙述如下。

2011SIMULIA中国区用户大会1

轮胎模型均为周期性的组合花纹轮胎模型[1]，总体建模方法如图1所示。限于篇幅，

本文不再对轮胎模型的建模方法进行详细介绍，感兴趣的读者可参阅笔者曾发表的相关文献[1-5]。本文共建立了三种组合花纹轮胎模型PT-A、PT-B和PT-C（如图2所示），分别对应三种花纹样式A、B和C；其中样式A为对称花纹（中心对称），样式B为单导向花纹（镜面对称），样式C为一款非对称花纹。

图1.周期性的组合花纹轮胎模型及其总体建模策略.

(a)模型PT-A(b)模型PT-B(c)模型PT-C

图2.用于滑水仿真的3个组合花纹轮胎模型.

(a)优化前(b)优化后

图3.用于滑水仿真的水/空气复合模型.

水/空气复合模型直接使用Abaqus/CAE建模得到，并将离散为Euler单元，用于描述

液体（水）可能填充的空间。经过多次试算后，根据水流的流动方向、流速分布以及水可

能填充区域的范围，对模型的几何形状进行了优化。优化前和优化后的水/空气复合模型的

22011SIMULIA中国区用户大会

几何形状及网格划分分别如图3(a)、(b)所示。与优化前的模型（图3(a)）相比，优化后的

模型（图3(b)）在水流可能流向的方向尤其是在高度方向上，增加了Euler网格的空间范

围，在水流不可能流向的方向删除了不必要的网格单元，并且在流体运动比较复杂的区域

进行了更加精细的网格划分。借助这些措施，优化后的模型将得到更加精确的计算结果，

同时避免了计算时间的大幅增加。后文进行滑水仿真时均使用优化后的水/空气复合模型。

轮辋和路面均使用解析刚体进行模拟，非常简单，不予详述。

3.轮胎滑水仿真方法

在进行有限元仿真分析时，采用“先隐式后显式”[6]的求解策略。所谓“先隐式”是

指先使用Abaqus隐式求解器进行轮胎装配、充气、加载和滚动工况的模拟，即将轮胎、轮

辋和路面组装在一起，定义合理的位移/力边界条件和接触边界条件，使用Abaqus/Standard进行轮胎装配、充气、加载和自由滚动工况的求解，并得到给定滚动速度下的稳定的应力

应变场。关于此过程的模拟方法，可参考相关论文[1,7]

获悉详细描述。所谓“后显式”是指

Abaqus显式求解器进行轮胎滑水工况的模拟，