0 引 言

排气歧管是发动机主要受热零件，与高温燃气直接接触，承受较高的温度，工作环境恶劣。如果排气歧管受热时引起的变形过大，会造成漏气，使内燃机无法正常工作。因此对排气歧管热固耦合分析很必要。

 

1 介 绍

进行热固耦合分析，必须要获得包括排气歧管内外流场的温度与对流换热系数边界，再使用 FEA 软件进行热应力、密封性和残余应力应变分析等。首先通过 BOOST 计算出排气歧管的进出口边界条件；然后将边界条件应用于 CFD计算中，通过 FIRE 计算排气歧管的内流?。ㄋ蔡┮约巴饬鞒。ㄎ忍?，得到歧管内外壁面的温度以及对流换热系数分布，进行时域平均后，将结果映射到 FEM 的网格上；最后通过 Abaqus 对歧管总成进行了热固耦合分析，对设计方案进行评估。排气歧管耦合分析考虑了排气歧管螺栓预紧力与热应力，主要分析以下几项内容：

l排气歧管温度场分布；

l排气歧管垫片密封性分析；

l热应力分析与残余应力应变分析。

 

其主要分析流程见下图 1

 

 

图 1 分析流程图

 

2 热固耦合分析流程

2.1 Boost 性能仿真

下图 2 是发动机的 BOOST 的计算模型。环境气体从边界 SB1 吸入，通过管 1 到空滤器 CI1，再通过管 2 进入涡轮增压器，出来再通过管 3 导向中冷器 CO1，再经连接管4进入进气歧管的谐振腔，用PL1代替，考虑到 EGR 的存在，容腔的一部分容积转化成连接管 16，17。管子 5—8 代表进气歧管和气道，通过气道将气引到气缸 C1—C4。

 

管子 9—12 代表排气歧管和气道，通过气道将气引到容腔 PL2，再通过管 13 连接到增压器 TC1，经管 14 连接到消声器 PL3，然后经管 15，通过 SB2 流向大气。

 

 

图 2 Boost 计算模型图

 

随着发动机工况的改变，燃烧放热规律也是随着变化的，并且燃烧特性的好坏对发动机性能有很重要的影响。热力学循环计算需要输入燃烧热量释放率曲线。在数学上，该曲线可以通过 VIBE 函数来模拟，VIBE函数的主要参数为燃烧区间和形状参数。燃烧热量释放率曲线影响缸内压力和温度，从而就影响了循环效率、爆发压力和排气温度。计算获得相应位置的瞬态流量、压力和温度分别作为歧管的进出口边界条件。

 

2.2 FEA 模型创建

分析模型包括模拟缸盖、排气歧管螺栓和排气歧管。使用 Hypermesh 对模型进行网格划分。在传热分析中实体单元类型选择DC3D4，面单元选择 DS3。强度分析时单元类型选择 C3D10M。需要注意 FEA 模型的坐标和 CFD 分析的完全一致。排气歧管耦合分析网格没有特别要求，只要保证接触边界的节点一一对应即可。

 

温度场计算要对模型的相关部件施加CFD 计算结果热边界，因此，需要将排气歧管内腔和外壁面的面网格提取出来，并将这些面网格单独写成 inp 文件，以供 CFD 软件提取热边界。对于在前处理中提取的部件面网格单元，定义一个很小的厚度，一般定义0.001mm 即可。

 

2.3 CFD 分析

2.3.1 内流场 CFD 分析

（1）网格化分

利用 Fire 软件的自带的 Fame 工具划分网格，网格的尺寸为 2mm。为了有利于计算的收敛，进出口边界均沿法线方向延长了20 层（每层高度为 2mm）。最终的网格数目约为 16.5 万，98%以上为六面体，其余为四面体，五面体等的混合网格。

 

（2）瞬态计算

瞬态计算共计算了五个循环，前四个循环为了计算收敛。取最后一个循环中气体的温度和换热系数平均值作为热应力计算的边界，如图 3 和 4 所示 。

 

 

图 3 热应力计算温度边界示意图

 

 

图 4 热应力计算换热系数示意图

 

2.3.2 外流场 CFD 分析

（1）网格化分

外流场主要是稳态计算，模型包括排气歧管，隔热罩，增压器，弯管等。为了模拟实验室的环境，外流场的计算域取 3 倍发动机的尺寸，参考风扇的鼓风速度，设定外流场入口的风速为 8m/s[2]。下图 5 为外流场计算模型。

 

 

图 5 排气歧管外流场计算模型图

 

使用 FAME 生成以六面体为主计算网格?？悸堑奖诿娓浇谋呓绮阌跋欤诒诿嫔仙梢徊惚呓绮阃?。由于隔热罩厚度仅为 3mm 左右，因此需要进行局部加密。网格总数约为 112 万。

 

（2）瞬态计算

采用迎风离散格式，一阶隐式格式离散时间项，压力与速度耦合算法选择SILMPLE。设定管内空气流动为可压缩粘性湍流流动，空气为理想气体，湍流模型 k-z-f方程，使用混合壁面函数描述壁面附近边界层流体速度、压力等的分布，且要求贴近壁面的网格的 y+值在 11～200 之间[3]。残差小于 0.0001。外流场的计算为稳态计算，为了保证计算收敛，一般计算 3000 个迭代步。具体边界条件如下所示：

 

进口边界条件：设为速度边界，速度大小为 8m/s；出口边界条件：设为梯度为零。排气歧管外壁面的换热系数和温度分别如图 6、7 所示，平均温度为 476.6K，平均换热系数为 76.5W/（m2K）。

 

 

图 6 排气歧管外壁面换热系数分布示意图

 

 

图 7 排气歧管外壁面气体温度分布示意图

 

2.3.3 CFD 瞬态计算结果的时域平均

瞬态计算（时间步或曲轴转角）在每个时刻都会有一个计算结果，而 FEM 计算只需要一个稳态的热边界结果，因此，瞬态计算的结果不能直接应用于 FEM 计算，必须首先对其进行时域平均。

 

 

 

2.3.4 CFD 结果与 Abaqus 有限元模型映射

映射是指 CFD 计算与 FEM 计算之间，在流固交界面上热数据的交换，它实现了FIRE 和 Abaqus 之间的流固耦合。热边界来源与 CFD 模拟结果，其获取方法是将排气歧管网格模型的面网格导出inp 格式（包含部件节点、单元信息即可）提供给 CFD 分析人员，由 CFD 分析人员提取出部件面单元的换热系数与温度文件，文件后缀名为.abaqinp。其中值得注意的两点：1，提出的换热系数与温度信息要附在单元上。2，所有涉及到数值的信息要标明单位，以免产生数量级上的错误。图 8、图 9 分别给出了两类文件的示例。

 

 

图 8 面网格单元 inp 文件示例

 

 

图 9 CFD 提供的热边界文件示例

 

2.4 温度场分析

2.4.1 温度场分析 FEA 设定

排气歧管垫片单元类型采用 DC3D15，实体单元类型选择 DC3D4，面网格单元类型采用 DS3。输入模型部件材料实测的随温度变化的弹性模量、柏松比和导热系数值，考虑到保密事宜，只展示排气歧管GGGSiMo51 的随温度变化弹性模量曲线，如下图 10 所示。

 

 

图 10 排气歧管 GGGSiMo51 的随温度变化弹性模量曲线

 

由于加载热边界时需要使用 include 语句，因此在进行模型搭建时一般进行手工编辑，而不是在 ABAQUS/CAE 中进行，首先在 Hypermesh 中对各接触对进行绑定的定义。

 

定义完绑定后导出 inp 格式，查看单元类型是否设定正确，在 inp 文件结尾添加分析步，如图 11 所示。

 

模型搭建完成后提交计算，利用批处理方式提交，使用内部命令 abaqus job=filename.inp。

 

2.4.2 温度场分析结果后处理

计算完成后，在 Abaqus CAE Visualization ?？橹醒≡?NT11 查看温度场分布云图，评价最高温度是否超过材料的温度极限值，并保留结果文件以供热应力计算。

 

图 12 为排气歧管外壁面的温度分布，可以看出外壁面最高温度 543.2℃。

 

 

图 12 排气歧管外壁面温度分布

 

图 13 为排气歧管内壁面的温度分布，可以看出外壁面最高温度 546.6℃。

 

 

图 13 排气歧管内壁面温度分布

 

为提高计算精度，可以进行多轮反复迭代，将 Abaqus 中的结算结果文件提取出内外壁面的温度结果.rpt 文件，提交给 CFD 分析工程师，进行多轮重复迭代：在 Abaqus后处理 Visualization ?？橹械慊?Create Display Group，分别选出并显示排气歧管内外壁网格；再点击 Report Filed Output，在Variable 中选择结果类型；在 Setup 中定义文件名称和存储地址，确认后即可生成温度结果.rpt 文件，如下图 14 所示。

 

 

图 14 排气歧管内壁面温度分布

 

蠕墨铸铁 TL047-SiMo4.5 排气歧管材料的最高温度限值为 800℃，排气歧管的外壁最高分析温度远小于此材料的许用温度限值,满足设计要求。

 

2.5 热应力分析

2.5.1 热应力分析 FEA 设定

排气歧管垫片单元类型采用C3D12MN，实体单元类型选择 C3D10M，在热应力分析中不需要壳体，删除面单元。热应力计算需要输入材料随温度变化的弹性模量、柏松比、线膨胀系数、密度、比热和塑性材料的应力应变曲线。同时添加排气歧管垫片的压缩卸载曲线，如下图 15 所示；对缸盖（或模拟缸盖）的一端施加固定约束，即约束 1、2、3 自由度；施加排气歧管螺栓的预紧力和热载。

 

 

图 15 垫片属性加载示例图

 

进行三个分析步热应力分析设定：

Step1：施加螺栓预紧力；

Step2：Fix boltlength，加载温度??；

Step3：卸载。

 

2.5.2 热应力分析结果后处理

计算完成后，在 Abaqus 后处理Visualization ?？橹薪胁僮?。图 16 为排气歧管热应力分布云图，可以看出，排气歧管的热应力集中区域在 EGR 连接管道过渡区域，最大热应力约在 129.4－200MPa，其它区域都低于 150MPa，都远低于材料的屈服极限（≥400MPa），满足静强度设计要求。

 

 

图 16 排气歧管热应力分布

 

2.5.3 密封性后处理

垫片上的面压应大于 15MPa，如图 17所示，满足设计要求。同时可以提取螺栓的振荡响应应力曲线。

 

 

图 17 垫片面压分布

 

2.5.4 残余热应力及残余应变结果后处理

图 18 为排气歧管残余热应力结果，可以看出排气歧管在卸载工况下残余热应力较小，在 61.18MPa 以下（螺栓预紧力加载区域不考察），静强度和疲劳都满足设计要求。

 

 

图 18 排气歧管残余热应力分布

 

图 19 为排气歧管残余应变结果，可以看出排气歧管在卸载工况下最大残余应变为 1.775%。根据 AVL 的标准，该材料的累积塑性应变最大限值为 2％，可以判断最大残余应变满足设计要求。

 

 

图 19 排气歧管残余应变分布

 

3 结 语

根据以上的分析可知，排气歧管热应力、密封性等都满足设计要求，建议采用原设计方案，并进行排气歧管支架有限元分析，重点考虑排气歧管搭载各附件下的子系统在外界加速度激励的应力状况和疲劳安全系数。

 

江淮汽车动力总成研究院 CAE 部通过使用 AVL.Boost、Fire 以及 Abaqus 软件，已经建立起较为成熟的固耦合分析能力，能对排气歧管的温度场、热应力、密封性、螺栓响应力、排气歧管的残余热应力及残余应变进行评估。

 

为了提升热固耦合分析的置信度，JAC通过近 2 年的时间，完成了发动机总成近 30种受温度影响部件材料的随温度变化性能测试，包括随温度应力应变曲线和弹性模量等力学性能，以及包括随温度变化的热膨胀系数、比热、热传到率和密度等物热性能。

 

资料来源：达索官方