单位文秘网 2021-07-04 08:07:38 点击: 次
摘要: 以某型大功率船用柴油机为研究对象,介绍Abaqus在船用柴油机开发过程中整机和零部件(固定件、运动件和受热部件等)设计仿真分析方面的应用.结果表明,基于有限元仿真分析的现代设计方法在船用柴油机开发中起到非常重要的作用.
关键词: 船用柴油机; 仿真分析; Abaqus
中图分类号: U664.121; TB115.1文献标志码: B
引言
CAE作为现代设计手段和核心技术,一方面能够从结构、热、流体和电磁等多个物理场综合反映柴油机的内在工作机理,比较分析对标机型的性能,预测产品设计的合理性;另一方面,克服传统设计方式的不足,缩短研发周期,降低产品开发成本,极大提升产品研发效率[1],可以满足船用柴油机设计开发仿真分析的需求.
Abaqus作为一款功能强大的有限元分析软件,不但能分析复杂的固体力学和结构力学系统,还能够驾驭庞大复杂问题和高度非线性等问题,并可以进行单一零部件的力学和多物理场的分析,以及系统级的分析和研究.
本文介绍在大功率船用中速柴油机开发设计阶段,利用Abaqus进行主要零部件结构仿真分析的情况,分固定件、运动件和受热零部件等3种类型,涉及的零部件主要包括机体、连杆、缸盖和缸套等.
1整机变形分析
本文基于新开发的大功率船用中速柴油机,其主要参数见表1.
1.1整机有限元分析
建立整机模型,见图1,主要包括机架、主轴承、缸套、机架侧盖板、凸轮轴盖板、缸盖、缸盖螺栓和油底壳等,单元总数为854 000个,节点总数为1 370 000个,总质量为58 000 kg.各部分的有限元模型及其参数见表2.
1.2机体静变形
考察机体在自重作用下的变形情况,可以得到在不同支撑位置的机架纵向相对变形,基准点为曲轴中心线,测量点为轴承座(曲轴中心线投影)位置.通过结果比较,选取合适的主机支撑方案.整机支撑位置分析见图2. 该模型还可以用于进行整机的振动分析,模拟真实工作条件下整机的振动特性,以证明所开发的柴油机符合船级社的规范要求.由于本文主要介绍柴油机研发过程中Abaqus的使用情况,在此不再赘述.图 2整机支撑位置分析2固定件静力学分析
以机架静力学分析为例,考虑到建立全部的模型需要的计算代价较大,因此选取2个缸段建立有限元模型,见图3,其中包含缸盖(简化)、垫圈、缸盖螺栓、缸套、机架、轴瓦、主轴承盖和主轴承联接螺栓等,进行变形和应力分析、关键部位的高周疲劳分析、接触分析(接触分离和相对滑动)和螺栓动态应力分析等4个方面的分析.
图 3机架应力分析有限元模型
2.1边界条件
一端约束x,y和z等3个方向,另一端用EQUATION约束,支座底平面全部固定,见图4.除了施加位移约束边界条件外,还考虑各零部件之间的间隙和过盈配合等接触边界条件.
图 4机架约束边界条件
2.2载荷
主要有气体压力、螺栓预紧力和轴承载荷等.
2.3计算工况
考虑80%额定转速,100%额定转速,110%超转速和110%超负荷4种工况.
2.4计算结果
2.4.1机架应力结果
选取关键部位的最大主应力和最小主应力进行应力结果评价.机架应力云图见图5.
图 5机架应力云图
2.4.2接触评估
评估所有载荷步下机架各部件间的最大接触压力.主轴承盖接触压力分布云图见图6.
图 6主轴承盖接触压力分布云图
2.4.3接触间隙
在全部载荷步中,最大的接触开口分布见图7.对于初始有间隙的部分,可允许接触开口的发生,但在主轴承盖联接螺栓与主轴承盖的接触部位不允许有接触滑移发生.图 7主轴承盖接触开口分布
2.4.4螺栓动态应力评估
螺栓动态应力的评估一般采用参考规范的方法,分别对轴承盖横向螺栓、轴承盖垂向螺栓和底脚螺栓进行动态疲劳强度分析.经计算,螺栓的安全因数分别为3.6,4.11和13.5(一般要求该值大于2.0即可).螺栓动态应力安全因数见图8.
图 8螺栓动态应力安全因数
3运动件力学分析
以连杆为例进行分析.连杆是柴油机最重要的运动件之一.连杆的长度变形将直接影响燃烧室压缩比,弯曲变形会导致摩擦加剧,引起相连接部件的损坏.对于船用连杆,其主要组成部件有衬套、杆身、连杆大端、连杆大端轴承盖、轴瓦和螺栓等.
3.1计算目的
连杆主要承受的载荷有气体作用力、运动件的惯性力、重力、摩擦力、弯曲载荷和装配静载等.其中,气体力和惯性力导致连杆承受反复拉压作用,易对连杆造成疲劳损坏,气体力还使杆身受压,甚至导致屈曲失稳.
通过有限元计算,研究连杆在不同工况下的应力及变形,并以此为基础分析连杆的疲劳性能,评估轴瓦与连杆以及连杆各部分之间的接触特性,此外,还要进行连杆的屈曲分析和螺栓在动态受力情况下的应力估算.根据计算结果,对连杆的设计进行校核.
3.2计算内容
针对连杆主要开展下几方面的仿真分析:变形分析、应力分析、屈曲分析、接触分析、螺栓动态应力分析和HCF高周疲劳应力分析等.连杆模型见图9.
图 9连杆模型
3.3载荷工况
为全面分析连杆组在不同工况下的应力特性,并为后续的疲劳分析提供必要的数据,共定义5个载荷步,每个载荷步对应一种工况.
工况1100%螺栓预紧力,定义2个轴瓦之间接触时的过盈量;定义轴瓦和曲柄销座的过盈量.
工况2工况1+570 r/min下运动组件惯性力(包括连杆和活塞组件),此转速下柴油机的输出扭矩最大.
工况3工况1+600 r/min下运动组件惯性力(包括连杆和活塞组件),此转速为柴油机的额定转速.
工况4工况1+635 r/min下运动组件惯性力(包括连杆和活塞组件),此转速为柴油机110%超负荷.
工况5工况1+600 r/min下运动组件惯性力(包括连杆和活塞组件)+气体最大爆发压力.
3.4计算结果
与机架静力学分析结果类似,连杆分析结果主要包括最大主应力、最小主应力、von Mises应力、接触面压应力、接触面分离和连杆的变形等.依据相应规范对计算结果进行评估,确定设计方案的可行性.由于篇幅所限,不再展开.
4受热部件
受热部件主要包括气缸盖和缸套及气缸盖、缸套和活塞构成燃烧室空间.在气缸盖内一般有进排气道、冷却水腔、起动阀、安全阀和燃烧室等,并装有配气机构和喷油器等零部件.在内燃机工作过程中,气缸盖和缸套承受很高的交变机械负荷与热负荷,是内燃机工作条件最为恶劣的零部件之一.
对气缸盖和缸套进行预紧力、预紧力热应力、预紧力热应力爆发压力等3种工况下的热固耦合仿真分析,对气缸盖和缸套的应力和变形进行评估,判断设计是否满足要求.
4.1计算目的
评估气缸盖和缸套在热负荷和机械负荷作用下的变形和强度.
4.2计算内容
(1)热分析.(2)结构应力分析.装配载荷,热负荷,热分析和气体压力等.(3)疲劳分析.
4.3有限元模型
计算模型取一个单缸,见图10.模型包含缸盖、缸套、垫片、刮油环、水套、机架(取部分结构)、喷油器组件、起动空气阀组件、进排气阀和进排气阀座等部件.图 10缸盖和缸套有限元模型
4.4热分析
4.4.1初始热边界条件
确定相关零部件各部位的温度和对流换热系数,主要有缸盖底部燃烧面、缸套表面和冷却水腔的热边界条件,缸盖、缸套冷却水侧的换热系数,气阀和阀座处及排气和进气通道处的热边界条件.
4.4.2热分析结果
分别对不同零件进行热分析,得到热固耦合分析的热边界条件,即缸盖底部、排气阀、缸套上部和冷却水腔壁面的温度分布见图11~14.
图 11缸盖底部温度分布
图 12排气阀温度分布图 13缸套温度分布
图 14排气道区域的水侧温度分布4.5应力分析
4.5.1边界条件
4.5.1.1约束
对气缸体的其中一侧约束x方向,另一侧施加EQUATION对称性边界条件;约束气缸体侧面的y方向(坐标轴采用局部坐标轴,y轴垂直于气缸体表面);对气缸体底部约束z方向,见图15.
图 15位移约束边界条件定义
4.5.1.2间隙和过盈
由于涉及的零部件较多,且相互间配合关系不一,定义接触关系需要考虑各零部件之间的间隙或过盈配合.各个零部件的接触应力在计算中采用小滑移接触,摩擦因数取0.19.
4.5.1.3载荷
(1)螺栓预紧力.在缸盖和缸套的有限元计算中,除考虑缸盖螺栓预紧力的影响,还需考虑喷油器夹具螺栓、起动空气夹具螺栓和压力传感器夹具螺栓的预紧力.
(2)燃气压力.分别在缸盖底部的燃烧面、缸套和刮油环的上部施加P=2.50×107 Pa的爆发压力.对于缸套和刮油环的内表面,由于气体压力在竖直向下的方向上递减,所以在结构强度计算中施加线性方程的压力曲线.
进气阀座、排气阀座和起动空气阀座上施加燃气压力,取最大爆发压力.
4.5.2计算结果
主要考察各工况下,缸盖应力分布、缸套变形和垫片受力情况.
4.5.2.1缸盖应力
应力分析是评价缸盖和缸套安全性的基本内容,局部应力集中会引起缸盖的热裂.缸盖受到的应力比较复杂,一般考察缸盖在装配状态、热负荷状态、热负荷+爆压状态(见图16和17)和卸载状态下的最大或最小主应力.图 16热负荷+爆压状态下的最大主应力
图 17热负荷+爆发压力下的最小主应力
4.5.2.2缸套变形
(1)径向变形.在缸套变形影响因素以及控制缸套变形方面,研究发现引起缸套变形的各种影响因素中,温度的影响起主要作用,其次是螺栓预紧力,所以在对缸套变形的分析中,主要考察冷装配状态和热负荷状态下缸套经过傅里叶变换后的变形量,见图18.
(a)缸套示意(b)冷装配 (c)热负荷
图 18缸套变形
由图18可知,缸套的最大径向变形发生在第一活塞环位置处(h=89.5),是因为受到外侧冷却水孔的结构影响;最小径向变形发生在缸套与机体接触部位,变形呈圆形状态;在热负荷工况时,缸套变形较大,最大膨胀的部位发生在缸套顶部h=0;在缸套的上部h为0~340缸套变形呈圆形向外扩张;在缸套下部h为697~1038缸套的变形呈椭圆状态.
(2)径向变形沿轴向分布.在气缸套的轴线方向布置节点集,计算完成后取出这些节点的坐标,然后将其绘制在以缸套内壁处于顶平面上的点作为原点的局部坐标系下,得到缸套径向变形沿轴向分布图,基本可以反映缸套轴向变形的情况,见图19.
图 19傅里叶变换后各阶次幅值下的缸套变形
缸套变形经过傅里叶变换后可以看作是所有阶次下变形的叠加,不同阶次对缸套变形的贡献不同,通过上述傅里叶变换得到各阶幅值,找到对缸套变形影响较大的阶次.对冷装配工况和热负荷工况进行分析.图19中阴影区域为相似机型的缸套在热负荷工况下变形量,各曲线分别表示缸套在冷装配工况下和热负荷工况下的径向变形沿轴向分布,可以看出第2阶和第4阶的幅值较大.
4.5.2.3接触分析
分别考察缸盖垫片在冷装配、热负荷、热负荷+气体压力和卸载工况下的压力分布,见图20.
(a)冷装配工况 (b)热负荷工况(c)热负荷+气体压力工况 (d)卸载工况图 20各工况下的垫片压力
5结束语
基于仿真分析驱动的产品设计开发正逐渐成为普遍现象,有限元软件将会发挥越来越大的作用.随着CAE应用水平的提高,仿真结果也逐步接近实际状态,参考经验数据,许多计算结果无须经过试验验证.示例结果表明,Abaqus以其强大的非线性分析功能,完全能够满足船用柴油机自主研发过程中对结构仿真分析的需求.参考文献:
[1]张坚. 船用柴油机研发与CAE应用[J]. 船舶物资与市场, 2008(5): 2428.
(编辑武晓英)
(责任编辑:单位文秘网) )地址:https://www.kgf8887.com/show-139-57378-1.html
下一篇:汽车装配工艺模块化设计研究
版权声明:
本站由单位文秘网原创策划制作,欢迎订阅或转载,但请注明出处。违者必究。单位文秘网独家运营 版权所有 未经许可不得转载使用