排气歧管流固耦合计算.docxVIP

排气歧管流固耦合计算

张长岭王福志赵铮

（长城汽车股份有限公司动力事业部，河北省保定市）

摘要：本文采用流固耦合的方法（Fire和Abaqus软件）计算了排气歧管的热应力，首先利用Fire计算出流场边界的对流

换热系数和温度。再用Abaqus计算排气歧管的温度场和热应力。确定热应力最大位置。

关键词：流固耦合；排气歧管；热应力；温度场

1前言

随着内燃机性能的不断提高，其零部件的热负荷也随之增加。排气歧管作为发动机的重要受热件之一，具有

工作热负荷大，热应力高特点，容易产生热疲劳失效、断裂等问题。因此在发动机的设计阶段，对其排气歧管进行热应力分析是非常必要的。本文采用CFD与FEA结合的方法对发动机排气歧管进行热应力分析，并根据计算结

果来初步判断此排气歧管是否存在由于热疲劳而产生的失效、断裂等隐患。

2排气歧管CFD计算

2.1内流场计算与结果

首先考虑热负荷最恶劣的情况，即在发动机额定功率下的工况，通过Boost软件计算得到发动机在额定转速

6000rpm时一个工作循环下排气歧管各个出口的流量、温度随曲轴转角变化的曲线（如图2.1所示），作为Fire软件的瞬态边界条件来计算排气歧管的内流场，得到排气歧管内壁面的瞬态对流换热系数和环境温度，再用时间

平均的方法计算出内壁面稳态对流换热系数和环境温度。排气歧管内流场瞬态边界条件进口给流量和温度，出口给静压，其余为壁面边界条件，壁面温度给550degC。排气歧管内壁CFD结果如图2.2、2.3所示，这个结果将要

作为FEA模型的内壁面边界条件来计算排气歧管的温度场。

图2.1Fire瞬态计算的边界条件

图2.2气体的壁面温度图2.3壁面的换热系数分布

2.2外流场的计算与结果

外流场采用稳态计算方法，计算模型是模拟在实验室情况下，发动机前方布置一个风扇，以2.5kg/s的空气

流量吹向发动机（试验数据），对发动机进行冷却，出口采用静压即一个标准大气压。同样考虑热负荷最恶劣的情况，即在发动机额定功率下的工况，根据性能计算得到排气歧管外部壁面温度作为边界条件进行计算。当计算达

到收敛后将计算后的空间壁面温度与换热系数映射到有限元网格，进行热应力计算。

图2.4外流场网格模型

图2.5歧管外壁面温度图2.6歧管外壁面换热系数

3排气歧管有限元模型

排气歧管计算的FE模型如图3.1所示，有限元网格用Hypermesh完成，边界条件用Patran软件处理，最后

用Abauqs来求解温度场和热应力。其中采用十节点四面体网格。共有节点59787；单元32157个。

图3.1有限元模型

4温度场计算及结果

排气歧管温度场计算控制方程为：

????????????

ttt

?k???k???k???q

??

?x??x??x?y?z??z?

??

?t

热边界条件为第一类边界条件t=g(x,y)和第三类边界条件??

?k?ht?t，和缸盖接触壁面为第一类边界条件，

?f

n其余面给第三类边界条件。排气歧管内外壁的热边界条件由CFD计算得到。图4.1为排气歧管温度场的计算结果；

最高温度为661.9℃，出现在排气支管交汇处；最低温度为283.6℃，出现在与气缸盖接触的法兰盘外侧边缘，由

于此处间断地受到高温废气的加热，所以温度较低。由云图可知，排气歧管温度分布不均匀，存在较大的温度差

异，这就必然产生较大的热应力。

图4.1排气歧管的温度分布

5热应力计算结果

为了分析热负荷对排气歧管热应力、热变形影响，本文计算了排气歧管在只有温度场作用下的热应力、热变

形。计算出排气歧管温度场后，以温度为边界条件计算热应力以及塑变。其基本方程为：

?????k???

R（1）

t

?????D?????????

B???(2)

0

其中[K]为总体刚度矩阵，[D]为弹性矩阵，[B]为应变矩阵，[Rt]为总体载荷矩阵，???节点应力矩阵。如图5.1

为s.mises应力云图，图5.2为X向应力分布云图，图5.3Y向应力分布云图

图5.1排气歧管Von.mises应力云图图5.2排气歧管X向拉伸应力分布图

图5.3排气歧管Y向拉伸应力分布图

由上图可知排气歧管内表面大部分区域受到压应力作用，外表面大部分区域主要受到拉应力作用，排气歧管

外表面受到拉应力最大值达到80.7MPa，主要出现在圆角处。

6热变形计算结果

整个排气歧管的温度较高，膨胀后热变形较大，最大变形出现在第四缸法兰