基于Abaqus某柴油机排气歧管热固耦合分析.docVIP

基于Abaqus某柴油机排气歧管热固耦合分析 　　摘要： 排气歧管在工作状态下会承受很高的环境温度，而在较高的环境温度下金属材料的力学性能和热学性能通常会发生较大变化.通过添加实测的各种部件材料随温度变化的性能参数，介绍某柴油机排气歧管热固耦合分析过程，展现使用AVL系列软件与Abaqus软件进行的排气歧管热固耦合分析流程. 　　关键词： 排气歧管； 热固耦合； 温度； 力学性能； 热学性能 　　中图分类号： U464.138；TB115.1文献标志码： B 　　引言 　　排气歧管是发动机主要受热零件之一，与高温燃气直接接触，承受较高温度，工作环境恶劣.如果排气歧管受热时引起的变形过大，会造成漏气，使内燃机无法正常工作.因此，很有必要对排气歧管热固耦合分析. 　　进行热固耦合分析，必须获得包括排气歧管内外流场的温度和对流换热系数边界等，再使用有限元软件进行热应力、密封性和残余应力应变分析等工作. 　　首先通过BOOST计算出排气歧管的进出口边界条件；然后将边界条件应用于CFD计算中，通过Fire计算排气歧管的内流场（瞬态）以及外流场（稳态），得到歧管内外壁面的温度以及对流换热系数分布，进行时域平均后，将结果映射到有限元模型的网格上；最后通过Abaqus对歧管总成进行热固耦合分析，评估设计方案. 　　排气歧管耦合分析考虑排气歧管螺栓预紧力和热应力，主要分析以下3项内容：（1）排气歧管温度场分布；（2）排气歧管垫片密封性分析；（3）热应力分析与残余应力应变分析. 　　主要分析流程见图1. 　　图 1分析流程 　　1热固耦合分析流程 　　1.1BOOST性能仿真 　　发动机BOOST计算模型见图2. 　　图 2BOOST计算模型 　　环境气体从边界SB1吸入，通过管1到空滤器CI1，然后通过管2进入涡轮增压器，再通过管3导向中冷器CO1，最后经连接管4进入进气歧管的谐振腔，用PL1代替.考虑到EGR的存在，容腔的一部分容积转化成连接管16和17.管子5～8代表进气歧管和气道，通过气道将气引到气缸C1～C4.管子9～12代表排气歧管和气道，通过气道将气引到容腔PL2，再通过管13连接到增压器TC1，经管14连接到消声器PL3，然后经管15，通过SB2流向大气. 　　随着发动机工况的改变，燃烧放热规律也随着变化，并且燃烧特性的好坏对发动机性能有很重要的影响.热力学循环计算需要输入燃烧热量释放率曲线.在数学上，该曲线可以通过VIBE函数模拟，VIBE函数的主要参数为燃烧区间和形状参数.燃烧热量释放率曲线影响缸内压力和温度，从而影响循环效率、爆发压力和排气温度等.计算获得相应位置的瞬态流量、压力和温度等，分别作为歧管的进出口边界条件. 　　1.2有限元模型创建 　　分析模型包括模拟缸盖、排气歧管螺栓和排气歧管.使用HyperMesh对模型进行网格划分.在传热分析中实体单元类型选择DC3D4，面单元选择DS3.强度分析时单元类型选择C3D10M. 　　需要注意有限元模型的坐标须与CFD分析完全一致.排气歧管耦合分析对网格没有特别要求，只要保证接触边界的节点一一对应即可. 　　温度场计算要对模型的相关部件施加CFD计算结果热边界，因此需要将排气歧管内腔和外壁面的面网格提取出来，并将这些面网格单独写成inp文件，以供CFD软件提取热边界.对于在前处理中提取的部件面网格单元，定义一个很小的厚度，一般定义0.001 mm即可. 　　1.3CFD分析 　　1.3.1内流场CFD分析 　　（1）网格划分.利用Fire软件自带的Fame工具划分网格，网格的尺寸为2 mm.为有利于计算收敛，进出口边界均沿法线方向延长20层（每层高度为2 mm）.最终的网格数目约为16.5万个，98%以上为六面体，其余为四面体和五面体等的混合网格. 　　（2）瞬态计算.瞬态计算共计算5个循环，前4个循环的目的是为计算收敛.取最后一个循环中气体的温度和换热系数平均值，作为热应力计算的边界，见图3和4. 　　图 3热应力计算温度边界 　　图 4热应力计算换热系数 　　1.3.2外流场CFD分析 　　（1）网格划分.外流场主要是稳态计算，模型包括排气歧管、隔热罩，增压器和弯管等.为模拟实验室的环境，外流场计算域取3倍发动机的尺寸，参考风扇的鼓风速度，设定外流场入口的风速为8 m/s.[1]外流场计算模型见图5. 　　图 5排气歧管外流场计算模型 　　使用Fame生成以六面体为主的计算网格.考虑到壁面附近的边界层影响，在壁面上生成一层边界层网格.由于隔热罩厚度仅为3 mm左右，因此需要进行局部加密.网格总数约为112 万个. 　　（2）瞬态计算.采用迎风离散格式，1阶隐式格式离散时间项，压力与速度耦合算法选择SIMPLE