JournalofMechanicalStrength
2021ꎬ43(5):1184 ̄1190
DOI:10 16579/j.issn.1001 9669 2021 05 023
∗20200907收到初稿ꎬ20201107收到修改稿ꎮ国防科技重点实验室基金项目(6142212190307)资助ꎮ
∗∗邹萍萍ꎬ女ꎬ1987年8月ꎬ江西抚州人ꎬ汉族ꎬ江铃汽车股份有限公司工程师ꎬ硕士ꎬ主要研究方向为发动机缸盖疲劳强度与寿命研究ꎮ
∗∗∗景国玺(通信作者)ꎬ男ꎬ1984年1月ꎬ甘肃会宁人ꎬ河北工业大学机械工程学院教授ꎬ博士ꎬ主要从事发动机结构疲劳强度研究工作ꎮ
基于Sehitoglu模型的发动机气缸盖热机疲劳寿命预测∗
THERMAL ̄MECHANICALFATIGUEPREDICTIONOFENGINE
CYLINDERHEADBASEDONTHESEHITOGLUMODEL
邹萍萍∗∗1㊀景国玺∗∗∗2ꎬ3
㊀曾小春1㊀
骆旭薇1㊀袁晓军1㊀魏㊀涛1㊀林宇星1
(1.江铃汽车股份有限公司动力总成开发部ꎬ南昌330001)
(2.河北工业大学机械工程学院ꎬ天津300400)
(3.天津市新能源汽车动力传动与安全技术重点实验室ꎬ天津300400)
ZOUPingPing1㊀
JINGGuoXi2ꎬ3㊀
ZENGXiaoChun1㊀
LUOXuWei1㊀
YUANXiaoJun1㊀
WEITao1㊀
LINYuXing1
(1.PowertrainEngineeringDepartmentꎬJianglingMotorsCorporationꎬLimitedꎬNanchang330001ꎬChina)
(2.SchoolofMechanicalEngineeringꎬHebeiUniversityofTechnologyꎬTianjin300400ꎬChina)(3.TianjinKeylaboratoryofPowerTransmissionandSafetyTechnologyforNewEnergyVehiclesꎬ
Tianjin300400ꎬChina)
摘要㊀以某发动机气缸盖为研究对象ꎬ开展了气缸盖本体材料力学性能测试与表征㊁气缸盖热机疲劳寿命预测研究工作ꎮ采用流固耦合方法ꎬ获得了准确的气缸盖热边界和温度场结果ꎬ温度场计算值与实测结果相符合ꎮ依照整机热冲击试验规范ꎬ采用Sehitoglu模型对缸盖的热机疲劳寿命进行了预测ꎬ仿真结果显示ꎬ缸盖最低寿命出现在第二缸火力面排气侧鼻梁区ꎬ主要由环境损伤引起ꎬ最低寿命为6860次ꎬ满足设计要求ꎬ该气缸盖顺利通过整机热机疲劳台架试验考核ꎮ
关键词㊀发动机㊀气缸盖㊀热机疲劳㊀寿命预测中图分类号㊀TK422
Abstract㊀Theresearchworkwascarriedoutonthemechanicalperformancetestingandcharacterizationofanengine
cylinderheadmaterialꎬandthepredictionofthecylinderheadthermalenginefatiguelife.Thefluid ̄solidcouplingmethodwas
usedtoobtainaccurateresultsofthethermalboundaryandtemperaturefieldofthecylinderhead.Thecalculatedtemperaturefieldisconsistentwiththem
easuredresults.AccordingtothethermalshocktestspecificationsofthewholemachineꎬtheSehitoglumodelwasusedtopredictthethermalfatiguelifeofthecylinderhead.Thesimulationresultsshowthattheminimum
lifeofthecylinderheadoccursinthenosebridgeareaontheexhaustsideofthesecondcylinderfiresurfaceꎬwhichismainlycausedbyenvironmentaldamage.Theservicelifeis6860ꎬwhichmeetsthedesignrequirements.Atthesametimeꎬthecylinderheadhassuccessfullypassedtheheatenginefatiguebenchtestofthewholemachine.
Keywords㊀InternalcombustionengineꎻCylinderheadꎻThermalmechanicalfatigueꎻLifepredictionCorrespondingauthor:JINGGuoXiꎬE ̄mail:okjgx@163.comꎬTel:+86 ̄22 ̄60435125ꎬFax:+86 ̄22 ̄60435125Thepro
jectsupportedbytheFundofNationalDefenseScienceandTechnologyKeyLaboratory(No.6142212190307).Manuscriptreceived20200907ꎬinrevisedform20201107.
㊀㊀引言
近年来ꎬ随着排放法规日益严格和对轻量化设计的需求ꎬ发动机强化程度越来越高ꎮ随着发动机功率密度增加ꎬ气缸盖等结构件热机械负荷突出ꎬ疲劳寿命
问题已成为研究的热点[1]1 ̄12[2 ̄3]和发动机产品研制过程中的拦路虎[4]ꎮ
发动机气缸盖设计开发必须满足高周机械疲劳和
低周热机械疲劳(ThermalMechanicalFatigueꎬTMFꎬ简称热机疲劳)设计要求[5]ꎮ针对气缸盖高周疲劳评估
㊀第43卷第5期邹萍萍等:基于Sehitoglu模型的发动机气缸盖热机疲劳寿命预测1185㊀㊀
问题ꎬ国内外已开展较为广泛的研究[6 ̄8]ꎬ为缸盖疲劳强度评估提供了重要依据ꎮ由于低周热机疲劳损伤机理复杂ꎬ热机疲劳寿命预测问题是当前研究的难点和焦点[1]1 ̄12[9]ꎮ
国内采用总应变法㊁修正的能量法等方法对气缸
盖低周疲劳寿命进行了预测ꎬ获得了热冲击载荷下气缸盖火力面损伤结果ꎮ由于气缸盖鼻梁区呈现典型的反相位热机疲劳载荷特性ꎬ常规低周疲劳寿命预测模型精度较低ꎬ目前ꎬSehitoglu模型[10]1769 ̄1783被行业认为精度高最权威的TMF模型ꎬ本论文将以该理论模型为依据进行气缸盖材料热机疲劳性能测试与表征及热机疲劳寿命预测研究ꎮ
1㊀热机疲劳分析流程汽车发动机盖
气缸盖热机疲劳寿命预测流程见图1所示ꎬ主要包括材料力学性能测试与表征㊁低周循环载荷条件下温度和应力应变行为分析和热机疲劳寿命评估等方面ꎮ
测试气缸盖铸铝材料不同温度下导热率ꎬ为温度
场预测提供材料边界ꎻ开展材料低周疲劳试验ꎬ获取材料弹塑性本构模型参数ꎬ为气缸盖热机耦合应力有限元分析提供边界ꎻ开展不同温度和相位条件下材料热机疲劳试验和蠕变试验等ꎬ获得材料应力应变行为与热机疲劳寿命之间的关系ꎬ进而构建Sehitoglu热机疲劳寿命预测模型和完成模型参
数识别ꎮ
采用AVL ̄Boost软件建立整机一维热力学仿真模型ꎬ得到发动机在怠速㊁额定功率负荷下进排气道气体温度㊁质量流量等参数ꎬ为三维缸内燃烧分析提供计算边界ꎮ然后ꎬ通过AVL ̄Fire软件进行缸内燃烧分析和缸体缸盖水套分析计算ꎬ并将燃烧和水侧换热边界映射至有限元模型ꎬ采用Abaqus软件对气缸盖进行温度场分析ꎬ通过多轮次耦合计算使得温度结果趋于稳定ꎬ并与实测结果相符合ꎮ
依照整机热冲击试验规范ꎬ建立基于材料弹塑性非线性的有限元分析模型ꎬ获得怠速 ̄额定转速交变工况条件下气缸盖应力应变行为ꎬ进而基于Sehitoglu模型实现气缸盖热机疲劳损伤和寿命预测
图1㊀气缸盖热机疲劳寿命评估流程Fig.1㊀ThecylinderheadTMFlifeevaluationprocess
2㊀热机疲劳损伤模型
本文采用Sehitoglu模型[10]1769 ̄1783ꎬ该模型考虑机
械疲劳㊁环境效应和蠕变三种不同的损伤机制ꎮ疲劳(机械)损伤主要指由室温条件下应变幅主导的疲劳损伤ꎬ高温条件下的损伤主要包括环境损伤和蠕变损伤ꎮ总损伤与机械损伤㊁环境损伤和蠕变损伤之间的关系见公式(1)ꎬ寿命关系式描述见公式(2)ꎮ
DTMF=Dfat+Denv+Dcreep
(1)
式中ꎬDTMF为总损伤ꎬDfat为机械疲劳损伤ꎬDenv为环境损伤ꎬDcreep为蠕变损伤ꎮ
1/Ntotalf=1/Nfatf+1/Nenvf+1/Ncreepf(2)
式中ꎬNtotalf为总寿命ꎬNfatf为机械寿命ꎬNenv
为环境寿命ꎬNcreepf
为蠕变寿命ꎮ在热机疲劳循环中ꎬ采用应变 ̄寿命方法来描述疲
劳行为ꎬ疲劳寿命采用Manson ̄Coffin方程来描述
㊀1186㊀机㊀㊀械㊀㊀强㊀㊀度2021年㊀
εmech2
εemech2
εplmech2
σᶄfE
(2Nfat
)
+εᶄf(2Nfatf)
(3)
式中ꎬΔεmech的机械应变幅ꎬσᶄf为疲劳强度系数ꎬE为弹性模量ꎬb为疲劳强度指数ꎬεᶄf为疲劳延展性系数ꎬ
c为疲劳延展性指数ꎮ这些参数通过室温条件下应变寿命试验来确定ꎮ
环境损伤主要认为在环境温度影响下氧化等现象
诱导裂纹扩展ꎬ不断将新鲜的金属材料暴露在环境中ꎬ环境损伤寿命可描述为
1Nenvf=hcrδoBΦoxKeffpéëêêùû
úú-1/β
2Δεmech()
2/β()+1
ε mech1-α/β()
(4)
式中ꎬhcr为临界氧化层厚度ꎬδ0㊁α㊁β和B均为材料常
数ꎬε mech为机械应变速率ꎬKeffp为有效氧化系数ꎬΦ
ox
为环境损伤相位调整系数ꎮ
考虑载荷作用历程ꎬ蠕变寿命描述如下
1Ncreepf
=Φcreep
ʏ
Ae-ΔH/RT()α1σ+α2σHKæèçöø÷mdt(5)
式中ꎬσ-是有效应力ꎬσH
是流体静力学应力ꎬK是拖拉压力ꎬα1和α2是应力作用系数ꎬ表示在拉伸和压缩过程中发生损伤的相对量ꎬΦ
creep
是蠕变损伤相位调整
系数ꎬA和m是材料常数ꎬt是时间ꎮ
3㊀气缸盖材料力学性能测试与表征
为获得热机疲劳模型材料参数ꎬ设计了一系列低周应变寿命曲线试验和热机疲劳试验ꎬ试验过程分别见图2和图3ꎬ典型试验载荷见图4所示ꎮ
低周疲劳试验基于Instron8801疲劳试验机开展ꎬ应变比为-1ꎬ测试温度包括室温㊁150ħ㊁200ħ㊁250ħ和300ħꎬ测试应变速率为1%/sꎮ另外ꎬ在250ħ和
300ħ下设计低应变速率0 1%/s下应变寿命试验ꎬ以考虑蠕变对疲劳寿命的影响ꎮ通过低周应变疲劳寿命试验得到的塑性应变与应力关系见图5ꎬ低周疲劳试验结果见图6ꎮ
材料TMF测试基于MTS810热机疲劳试验机进行ꎬ测试最高上限温度包括200ħ㊁250ħ和300ħꎬ每种温度考虑两种不同相位和约束系数ꎬ采用电感器加热方式加热ꎬ加热速率约为8ħ/sꎬ冷却使用压缩空气ꎬ最大温度保温60sꎬ截止寿命选定为3000次ꎮ图7为温度载荷100ħ~200ħꎬ约束比为1ꎬ下缸盖材料同相位TMF试验结果ꎬ图7a为不同循环下总应变㊁机械应变和热应变变化曲线ꎬ图7b为第二个㊁中值寿命和最后一个循环对应的应力应变滞后回环ꎬ图7b和图7c分别为中
值寿命对应循环内温度和应变与时间的变化曲线
图2㊀材料低周疲劳测试装置Fig.2㊀Lowcyclefatiguetest
devicet
图3㊀材料TMF测试装置Fig.3㊀TMFtestdevice
基于上述低周疲劳试验和热机疲劳实验结果ꎬ以中值寿命循环应力应变为载荷输入和Sehitoglu模型为依托ꎬ以试验寿命为目标ꎬ采用多目标优化算法进行材料参数多目标优化和识别ꎬ模型计算结果和实测对比见图8所示ꎬ在69个试验结果中ꎬ90%以上预测结果在3倍误差带范围之内ꎮ
4㊀热机疲劳寿命预测
4 1㊀缸内燃烧换热边界确定
本文缸盖热分析工作中缸内燃烧侧边界采用缸内燃烧数值模拟方法获得ꎮ缸内燃烧分析的湍流模型采用k ̄ζ ̄f模型ꎬ壁面传热采用StandardWallFunction模型ꎬPoisson
方程求解采用Simple数值方法ꎮ对于三维缸内燃烧分析ꎬ喷雾模型中ꎬ破碎模型采用Wave模型ꎬ蒸发模型采用Dukowicz模型ꎬ壁面油膜模型采用Walljet1模型ꎬ燃烧模型采用Ecfm ̄3Z模型ꎮCFD ̄FEA燃气侧耦合交接面包括缸盖火力面㊁进排气门㊁进排气道和缸孔内侧ꎮ通过缸内3D燃烧分析获得燃气温度和换热系数ꎬ并将内流场温度和换热系数边界映射到缸盖三维FEA模型上ꎮ
额定工况下CFD ̄FEA耦合单缸的燃气侧在一个发动机循环周期内平均温度和换热系数如图9㊁图10所示ꎬ其他三缸燃气测的温度和换热系数采用与第一缸相同的CFD结果ꎮ从图9可以看出ꎬ缸盖火力面温
㊀第43卷第5期邹萍萍等:基于Sehitoglu模型的发动机气缸盖热机疲劳寿命预测1187㊀
图4㊀材料疲劳试验载荷形式示意图
Fig.4㊀Schematicdiagramofmaterialfatiguetestload
form
图5㊀塑性应变 ̄应力关系曲线Fig.5㊀Plasticstrain ̄stress
curves
图6㊀应变 ̄寿命曲线Fig.6㊀Strain ̄lifecurves
度分布很不均匀ꎬ与燃烧过程密切相关ꎮ
4 2㊀水套流动换热边界确定
水套的主要作用是吸收气缸盖燃烧产生的热量ꎬ确保气缸盖不会因为温度过高而出现故障ꎮ同时ꎬ气缸盖的热负荷也取决于温度分布的均匀性ꎬ温度梯度越大ꎬ气缸盖热负荷越高ꎮCFD ̄FEA耦合水冷边界主要是指缸体和缸盖的水套壁面ꎮ与三维燃烧CFD分
析相似ꎬ三维水套分析发动机额定工况和怠速工况两个点ꎬ为FEA温度场分析的水冷边界提供壁面温度和换热系数ꎮ经三个轮次CFD ̄FEA耦合计算ꎬ缸体缸盖水套近壁面的温度和对流换热系数趋于稳定ꎮ缸体缸盖水套最终速度和对流换热系数分布如图11㊁图12所示ꎬ流速越高的区域对应的换热系数也越大ꎬ带走的热量越多ꎮ
4 3㊀温度场分析与试验对标
缸盖温度场分析CAE模型见图13ꎬ模型网格总数约410万ꎬ缸盖材料导热率结果采用实测值ꎮ缸内燃烧和冷却水腔换热边界由前述CFD结果提供ꎮ额定工况下机油温度为140ħꎬ怠速工况为100ħꎬ换热系数为l30W/(m2 K)ꎮ
额定工况下缸盖温度分布如图14所示ꎬ结果显示
缸盖火力面温度明显高于其他区域ꎮ受进气和排气温
㊀1188㊀机㊀㊀械㊀㊀强㊀㊀度2021年
图7㊀TMF试验结果Fig.7㊀TMFtest
results
图8㊀实测值和模型预测值对比Fig.8㊀Comparisonoftestandpredic
ted
results
图9㊀额定工况燃气侧温度分布
Fig.9㊀Gastemperaturedistributionatratedcondition
度的影响ꎬ缸盖排气侧结构温度高于进气侧ꎬ整个模型最高温度出现在第二缸排气侧鼻梁区ꎬ最高温度为
203ħꎮ
为验证仿真分析的准确性ꎬ本文采用硬度塞法对
发动机额定工况下各缸温度进行了试验测试ꎮ
每缸布
图10㊀额定工况燃气侧换热系数分布Fig.10㊀Gasheattransfercoefficientdistribution
atrated
condition
图11㊀额定工况下水套速度分布
Fig.11㊀Velocitydistributionofcoolantwateratra
tedcondition
置了12个测试点ꎬ四缸测点均采用相同的布置形式ꎬ各点分布位置如图15所示ꎮ四缸各测点温度测试平均值和仿真对比果如图16所示ꎬ结果表明仿真计算温度与测试温度趋势一致ꎬ相对误差基本在5%范围以内ꎮ
4 4㊀气缸盖热机疲劳寿命预测
气缸盖热机疲劳损伤发生部位主要位于火力面气