挠曲轴系横向振动计算及分析

船舶柴油机的轴系扭转振动的分析与研究【摘要】本文通过一些国内因轴系扭转振动而引起的断轴断桨的事故实例，来分析引起轴系扭转振动的主要原因，分析扭振主要特性，并提取一些减振和防振的基本控制措施。

【关键词】船舶柴油机轴系扭振危害分析措施在现代船舶机械工程中，船舶柴油机轴系扭转振动已经成为一个很普遍的问题，它是引起船舶动力装置故障的一个很常见的原因，国内外因轴系扭转而引起的断轴断桨的事故也屡见不鲜，随着科学水平的提高和航运业的发展，人们越来越重视船舶柴油机组的轴系扭转振动，我国《长江水系钢质船舶建造规范》和《钢质海船入级与建造规范》（简称《钢规》）和也均规定了在设计和制造船舶过程中，必须要向船级社呈报柴油机组的轴系扭转振动测量和计算报告，以此来表明轴系扭转振动的有关测量特性指标均在“规范”的允许范围内。

1 船舶柴油机轴系扭转振动现象简介凡具有弹性与惯性的物体，在外力作用下都能产生振动现象。

它在机械，建筑，电工，土木等工程中非常普遍的存在着。

振动是一种周期性的运动，在许多场合下以谐振的形式出现的，船舶振动按其特点和形式可分为三种，船体振动，机械设备及仪器仪表振动，和轴系振动。

船舶柴油机轴系振动按其形式可分为三种：扭转振动，纵向振动，横向振动。

柴油机扭转振动主要是由气缸内燃气压力周期性变化引起的，它的主要表现是轴系上各质点围绕轴系的旋转方向来回不停的扭摆，各轴段产生不相同的扭角。

纵向振动主要是由螺旋桨周期性的推力所引起的。

横向振动主要是由转抽的不平衡，如螺旋桨的悬重以及伴流不均匀产生的推力不均匀等的力的合成。

船舶由于振动引起的危害不但可以产生噪音，严重影响旅客和船员休息，还会造成仪器和仪表的损害，严重的时候甚至出现船体裂缝断轴断桨等海损事故，直接影响船舶的航行安全。

而在船舶柴油机轴系的三种振动中，产生危害最大的便是扭转振动，因扭转振动而引起的海损事故也最多，因此对扭转振动的研究也最多。

而且当柴油机轴系出现扭转振动时，一般情况下，船上不一定有振动的不适感，因此这种振动也是最容易被忽视的一种振动形式，一旦出现扭转振动被忽视，往往意味着会发生重大的事故。

梁的挠曲与振动文中关于梁的挠曲与振动的内容，可以按照以下方式进行论述：梁是一种常见的结构元件，主要用于支撑和传递载荷。

在工程应用中，梁的挠曲和振动问题是一个重要的研究方向。

本文将从梁的基本理论开始，介绍梁的挠曲和振动原理，以及相关的方法和应用。

一、梁的基本理论梁的基本理论包括梁的结构模型、梁的受力分析和梁的位移方程。

在这一部分中，我们将详细介绍梁的结构模型，如欧拉梁理论和蒙元梁理论，并推导出梁的受力分析和位移方程的表达式。

二、梁的挠曲理论梁的挠曲是指在受力作用下，梁发生的曲度变形。

这部分将介绍梁的弯曲应力和挠曲变形的计算，包括梁的弯矩-曲率关系、梁的挠度和梁的挠曲曲线等内容。

同时，还可以讨论梁的挠曲问题在工程中的应用，如在梁设计中的影响因素和设计原则。

三、梁的振动理论梁的振动是指梁在受到外力激励后产生的自由振动或强迫振动。

这部分可以介绍梁的振动特性，如梁的共振频率、振型和振动响应等内容。

同时，还可以讨论梁的振动问题在工程中的应用，如梁的减振措施和振动测试方法等。

四、梁的挠曲与振动的分析方法在梁的挠曲与振动分析中，有多种数值分析方法可以应用，如有限元方法和模态分析等。

本部分可以介绍这些分析方法的基本原理和步骤，并以实例说明其在梁的挠曲与振动分析中的应用。

五、梁的挠曲与振动的应用梁的挠曲与振动问题在工程中具有广泛的应用背景。

这部分可以以实例的形式介绍梁的挠曲与振动问题在不同领域的应用，如桥梁结构、航空航天和机械工程等，以及相应的安全性评估和优化设计等内容。

六、总结通过对梁的挠曲与振动的论述，我们可以得出结论，总结研究的结果和成果，并思考未来在这一领域的发展方向。

同时，还可以指出该领域的研究挑战和存在的问题，为进一步的研究提供借鉴和启示。

以上所述为梁的挠曲与振动文章的一个可能的论述框架，具体内容需要根据实际情况进行发挥和拓展，以充分满足文章的字数要求和信息表达的完整性。

第２０卷第６期２０２２年１２月动力学与控制学报ＪＯＵＲＮＡＬＯＦＤＹＮＡＭＩＣＳＡＮＤＣＯＮＴＲＯＬＶｏｌ．２０Ｎｏ．６Ｄｅｃ．２０２２文章编号：１６７２ ６５５３ ２０２２ ２０（６） １０１ ０５ＤＯＩ：１０．６０５２／１６７２ ６５５３ ２０２２ ０４７　２０２２ ０８ ２０收到第１稿，２０２２ ０９ ２８收到修改稿．国家自然科学基金资助项目（１２１７２２８１，１１９７２２８４），基础加强１７３计划基金项目（２０２１ ＪＣＪＱ ＪＪ ０５６５），陕西省科技创新团队资助（２０２２ＴＤ ６１）和陕西高校青年教师创新团队资助 通信作者Ｅ ｍａｉｌ：ｗｐｈｕ＠ｎｗｐｕ．ｅｄｕ．ｃｎ轴向运动功能梯度梁横向振动问题的保结构分析刘涛１　周洋忻２　胡伟鹏２（１．榆林市城市投资经营集团有限公司，榆林　７１９０００）（２．西安理工大学土木建筑工程学院，西安　７１００４８）摘要　轴向运动速度和材料的非均匀性对轴向运动功能梯度梁振动问题分析提出了严峻挑战．本文在简要回顾轴向运动功能梯度梁横向振动动力学模型基础上，基于无限维动力学系统的对称破缺理论和广义多辛分析方法，构造了横向振动模型的保结构数值格式，并在给定材料参数时给出了数值格式具有良好保结构性能的条件．分别采用微分求积法、复模态法和保结构方法分析横向振动模型的前六阶频率，发现保结构方法得到的频率结果与复模态法得到的结果吻合较好，在此基础上分析了微分求积法的主要误差来源，以指导微分求积法的改进，并为复杂动力学系统的数值求解提供了新途径．关键词　保结构，　轴向运动功能梯度梁，　对称破缺，　广义多辛，　横向振动中图分类号：Ｏ３０２文献标志码：Ａ引言功能梯度材料由于控制界面的成分和组织连续变化，使材料的热应力大为缓和，而在航空航天、机械工程、生物医药等领域应用广泛［１ ３］．智慧建造［４］这一全新概念的提出，使得传统单一均匀材料无法满足建筑设计工程的需求，因此，功能梯度材料将是未来实现很多智慧建造特殊功能的不二选择．作为智慧建造中的基本力学构件，功能梯度梁的动力学行为分析尤为重要．特别是在装配式智慧建造过程中，功能梯度梁运输及吊装过程的横向振动特性对运输和吊装过程的稳定性影响显著．Ｓａｎｋａｒ［５］基于Ｅｕｌｅｒ Ｂｅｒｎｏｕｌｌｉ梁理论，得到了横向载荷作用下功能梯度梁弹性范围内的解．Ｒｅｄｄｙ［６］基于ｖｏｎＫａｒｍａｎ几何非线性理论，建立了功能梯度梁的非线性Ｅｕｌｅｒ Ｂｅｒｎｏｕｌｌｉ梁模型和Ｔｉｍｏｓｈｅｎｋｏ梁模型．丁虎［７］、王忠民等［８］轴向运动功能梯度梁振动模型，并分别采用伽辽金法和微分求积法分析其振动特性，为本文分析功能梯度梁横向振动过程奠定了基础．刘金建等［９］基于Ｅｕｌｅｒ梁理论研究了轴向运动功能梯度粘弹性梁横向振动的稳定性问题．Ｂａｌｉｒｅｄｄｙ和Ｐｉｔｃｈａｉｍａｎｉ［１０］分析了时变轴向载荷作用下功能梯度梁振动特性及稳定性．从本质上讲，功能梯度梁的材料非均匀性和梁式结构的轴向运动均属于动力学对称破缺［１１］因素．对于含有对称破缺因素的动力学系统，本课题组基于多辛分析方法，建立了广义多辛分析方法［１２］这一保结构理论框架，并解决了一系列复杂动力学问题［１３ １６］．因此，本文将基于保结构思想，分析轴向运动功能梯度梁的横向振动频率特性，为功能梯度梁的横向振动控制提供参考．１　轴向运动功能梯度梁横向振动模型本节参考文献［８，９］，简要回顾轴向运动功能梯度梁横向振动动力学模型的建立过程．考虑一轴向运动的简支功能梯度矩形截面梁（图１），梁长度为Ｌ，横截面宽度为ｂ、高为ｈ，轴向运动速度为定常速度，大小为η．为了刻画材料特性沿界面高度方向的梯度，假定功能梯度材料有效杨氏模量和有效Copyright ©博看网. All Rights Reserved.动　力　学　与　控　制　学　报２０２２年第２０卷密度均为ｚ坐标的函数，即Ｅ（ｚ）和ρ（ｚ）．图１　轴向运动功能梯度梁的物理模型Ｆｉｇ．１　Ｐｈｙｓｉｃａｌｍｏｄｅｌｏｆｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｂｅａｍｗｉｔｈａｎａｘｉａｌｖｅｌｏｃｉｔｙ以含两种组分（如金属材料和陶瓷材料）的功能梯度材料为例，其有效材料参数可表述为：Ｅ（ｚ）＝（Ｅｃ－Ｅｍ）（ｚ／ｈ＋１／２）ｋ＋Ｅｍ ＝Ｅｍ［（βＥ－１）（ｚ／ｈ＋１／２）ｋ＋１］ρ（ｚ）＝（ρｃ－ρｍ）（ｚ／ｈ＋１／２）ｋ＋ρｍ ＝ρｍ［（βρ－１）（ｚ／ｈ＋１／２）ｋ＋１］（１）其中Ｅｃ，Ｅｍ，ρｃ，ρｍ分别为两种材料组分的物理参数，βＥ＝Ｅｃ／Ｅｍ，βρ＝ρｃ／ρｍ，ｋ为梯度指标．需要说明的是，从式（１）即可推导出功能梯度梁的中性层与几何对称中心重合．基于Ｅｕｌｅｒ Ｂｅｒｎｏｕｌｌｉ梁基本假设，依据文献［８］，功能梯度梁上任意点的位移可由梁轴线上任意点的轴向位移ｕ（ｘ，ｔ）和横向位移ｗ（ｘ，ｔ）表述：ｕｘ（ｘ，ｚ，ｔ）＝ｕ（ｘ，ｔ）－ｚ ｘｗ（ｘ，ｔ）＋ηｔｕｚ（ｘ，ｚ，ｔ）＝ｗ（ｘ，ｔ）（２）功能梯度梁上任意点的正应变分量和正应力分量分别为：εｘ＝ ｘｕ－ｚ ｘｘｗσｘ＝Ｅ（ｚ）εｘ＝Ｅ（ｚ）（ ｘｕ－ｚ ｘｘｗ）（３）由其描述的梁的应变能可表述为：Ｕ＝１２∫Ｌ０∫ＡσｘεｘｄＡｄｘ ＝１２∫Ｌ０［Ｄ１（ ｘｕ）２－２Ｄ２ ｘｕ ｘｘｗ＋ Ｄ３（ ｘｘｗ）２］ｄｘ（４）其中，Ａ为梁的横截面面积，并且：（Ｄ１，Ｄ２，Ｄ３）＝∫ＡＥ（ｚ）（１，ｚ，ｚ２）ｄＡ功能梯度梁上任意点两个方向的速度分量分别为：ｖｘ＝ ｔｕｘ（ｘ，ｚ，ｔ）　＝ ｔｕ（ｘ，ｔ）－ｚ ｔｘｗ（ｘ，ｔ）＋ηｖｚ＝ ｔｗ（ｘ，ｔ）＋η ｘｗ（ｘ，ｔ）（５）由此描述的梁的动能可表述为：Ｋ＝１２∫Ｌ０∫Ａρ（ｚ）（ｖ２ｘ＋ｖ２ｚ）ｄＡｄｘ＝　１２∫Ｌ０｛Ｉ１［（ ｔｕ）２＋η２＋２η ｔｕ＋（ ｔｗ）２＋　η２（ ｘｗ）２＋２η ｔｗ ｘｗ］－２Ｉ２η ｔｘｗ－　２Ｉ２ ｔｕ ｔｘｗ＋Ｉ３（ ｔｘｗ）２｝ｄｘ（６）其中（Ｉ１，Ｉ２，Ｉ３）＝∫Ａρ（ｚ）（１，ｚ，ｚ２）ｄＡ由哈密顿原理，忽略梁的轴向惯性力及其由轴向惯性力诱导的横向分布载荷项，并消去轴向位移项，得到轴向运动功能梯度梁横向振动方程：（Ｄ３－Ｄ２２Ｄ１） ｘｘｘｘｗ－（Ｉ３－Ｉ２Ｄ２Ｄ１） ｔｔｘｘｗ＋　Ｉ１（ ｔｔｗ＋η２ｘｘｗ＋２η ｔｘｗ）＝０ｗｉｔｈｗ（０，ｔ）＝０ｗ（Ｌ，ｔ）＝０ｘｘｗ（０，ｔ）＝０ ｘｘｗ（Ｌ，ｔ）＝{０（７）２　振动模型的近似对称形式及保结构离散引入如下中间变量： ｔｗ＝ ｘψ＝Ｄ１φ－Ｄ１Ｉ１χＤ２Ｉ２－Ｄ１Ｉ３， ｘｗ＝φ， ｘχ＝φ，并定义状态向量：ｚ＝（ｗ，φ，χ，ψ，φ）Ｔ，轴向运动功能梯度梁横向振动方程（不含边界条件）可以写成如下近似一阶对称形式：Ｍ ｔｚ＋Ｋ ｘｚ＝ ｚＳ（ｚ）＋τ（ｚ）（８）其中，Ｍ，Ｋ∈Ｒ５×５为反对称矩阵：Ｍ＝００００１０００００００００００００００－１００００，Ｋ＝　０Ｄ３－Ｄ２２Ｄ１００１Ｄ２２Ｄ１－Ｄ３０００００００Ｉ２Ｄ２Ｄ１－Ｉ３０００Ｉ３－Ｉ２Ｄ２Ｄ１０００００１０拟哈密顿函数为：Ｓ（ｚ）＝－１２［Ｉ１η２ｗ２＋（Ｄ３－Ｄ２２Ｄ１）φ２－２０１Copyright ©博看网. All Rights Reserved.第６期刘涛等：轴向运动功能梯度梁横向振动问题的保结构分析 Ｉ１χ２＋（Ｉ３－Ｉ２Ｄ２Ｄ１）ψ２］余项为：τ（ｚ）＝［－２Ｉ１ηψ，０，φ，０，０］Ｔ．与标准的多辛形式不同，近似对称形式含有如下对称破缺因素［１１］：①系数矩阵Ｍ，Ｋ及哈密顿函数Ｓ（ｚ）显含空间变量；②哈密顿函数梯度存在余项τ（ｚ）；③系数矩阵Ｍ非严格地反对称，因此将其分解Ｋ＝Ｋ＋Ｋ⌒０Ｄ３－Ｄ２２Ｄ１００１／２Ｄ２２Ｄ１－Ｄ３０００００００Ｉ２Ｄ２Ｄ１－Ｉ３０００Ｉ３－Ｉ２Ｄ２Ｄ１０－１／２－１／２００１／２０＋００００１／２００００００００００００００１／２１／２００１／２０ ①和③两个对称破缺因素引起的横向振动模型多辛结构残差和局部能量耗散均可以参照文献［１１］显式给出，第②个对称破缺因素在模拟仿真中的处理方式可参照参考文献［１７］进行．为避免与已有工作重复，在此不给出详细表达式和具体处理步骤，只在模拟结果中给出离散的多辛结构残差，以间接证明后续构造算法的有效性和保结构性能．在梁长度方位内（０≤ｘ≤Ｌ）采用空间步长进行均匀划分单元，并对系统采用时间步长进行Ｐｒｅｉｓｓｍａｎｎ离散，得到保结构差分格式：Ｍδ＋ｔｚｊｉ＋１／２＋Ｋδ＋ｘｚｊ＋１／２ｉ＝ ｚＳ（ｚｊ＋１／２ｉ＋１／２）＋τ（ｚｊ＋１／２ｉ＋１／２）（９）其中：ｚｊ＋１／２ｉ＋１／２＝１４（ｚｊｉ＋ｚｊｉ＋１＋ｚｊ＋１ｉ＋ｚｊ＋１ｉ＋１），δ＋ｘ，δ＋ｔ均为一阶前向差分．限于篇幅，格式的展开形式和消参后的形式不再给出，同时，离散的多辛结构残差和离散的局部能量耗散项也不再列出．需要强调的是，多辛结构残差是衡量格式保结构性能的重要依据，后续在数值结果中会详细讨论．３　数值算例为了将结果与文献［８，９］的部分结果进行对比，材料参数取值如下：Ｅｃ＝３９０ＧＰａ，Ｅｍ＝２１０ＧＰａ，ρｃ＝３９６０ｋｇ／ｍ３，ρｍ＝７８００ｋｇ／ｍ３．为保证数值格式的保结构性能，依照广义多辛理论［１２］，需要选取合适的时间步长使得在每一时间步内，离散的多辛结构绝对残差不超过差分格式的数值截断误差，即Δｉ≤ｏ（Δｔ，Δｘ），其中ｏ（Δｔ，Δｘ）为格式的数值截断误差．为了计算方便，忽略高阶项并取Δｔ／Δｘ＝０．５后，可以将数值截断误差上限估计值近似取为：ｏ（Δｔ，Δｘ）≤［ｏ］＝７Δｔ２（１０）在考虑梯度指标取值较大的情形下，确定容许的最大时间步长．取ｋ＝１０５，将时间步长取值从Δｔ＝０．００１ｓ逐渐增大，当式（１０）刚好严格满足时，得到最大允许时间步长为Δｔ＝０．０６４ｓ，此时的多辛结构残差与数值截断误差上限估计值之间的关系如图２所示．因此，在后续模拟过程中，取时间步长为Δｔ＝０．０５ｓ，空间步长为Δｘ＝０．１ｍ，就能保证所构造的格式具有良好的保结构性能．分别取ｋ＝０．００１，１００两种梯度指标，分别采用微分求积法（ＤＱＭ）［８］、复模态法（ＣＭＭ）［９］和保结构方法（ＳＰＭ）模拟轴向运动功能梯度梁的横向振动过程，得到梁的前六阶频率值如表１所示．从表１中不难发现，采用保结构分析方法得到的结果与复模态法得到的结果整体吻合较好．随着频率阶次升高，复模态法和保结构方法得到的频率结果明显低于微分求积法得到的结果．考察微分求积法的求解过程，可知微分求积法得到的结果产生以上偏差的主要原因在于以下两个方面：①在进行微分求积运算之前，将偏微分方程化为常微分方程过程中，只考虑了方程解的一阶频率分３０１Copyright ©博看网. All Rights Reserved.动　力　学　与　控　制　学　报２０２２年第２０卷量而忽略了高阶频率分量；②微分求积法采用非均匀网格离散，无法判断每一时间步内不等式（Δｉ≤ｏ（Δｔ，Δｘ））的满足情况，不具有评价其保结构性能的条件．复模态法在一定程度上克服了上述两方面的问题，故得到的结果与本文保结构方法得到的结果吻合较好．上述结果表明，复模态法和保结构方法在分析轴向运动功能梯度梁横向振动问题中均具有较好的数值精度．图２　轴向运动功能梯度梁的物理模型Ｆｉｇ．２　Ｅｖｏｌｕｔｉｏｎｏｆｔｈｅａｂｓｏｌｕｔｅｒｅｓｉｄｕａｌｏｆｔｈｅｍｕｌｔｉ ｓｙｍｐｌｅｃｔｉｃｓｔｒｕｃｔｕｒｅ表１　前六阶频率结果对比（Ｈｚ）Ｔａｂｌｅ１　Ｃｏｍｐａｒｉｓｉｏｎｏｆｔｈｅｆｉｒｓｔｓｉｘｆｒｅｑｕｅｎｃｉｅｓ（Ｈｚ）ｋＭｏｄｅＮｏ．ＤＱＭＣＭＭＳＰＭ１ｓｔ１８．０３８５１８．０３８５１８．０３８５２ｎｄ７２．５８０１７２．５３３９７２．５３３９０．００１３ｒｄ１６１．１９７５１６０．００１８１６０．００１８４ｔｈ２８９．８８０６２８６．２１４７２８６．２１４２５ｔｈ４５８．９３８０４５２．７３１１４５２．７０９６６ｔｈ６６６．２０３９６５９．９０１８６５９．３０８９１ｓｔ９．７８４９９．７８４９９．７８４８２ｎｄ３２．９０９１３２．２２８６３２．２２５９１００３ｒｄ８０．３６１０７８．４３９２７８．４２９８４ｔｈ１４８．１３１５１４４．３６１８１４３．８１００５ｔｈ２３７．２０２７２３１．２１５６２３０．９０３５６ｔｈ３４６．７７３８３３８．８０３３３３８．３２７１４　结论基于动力学系统的对称破缺理论和广义多辛分析方法，本文针对轴向运动功能梯度梁横向振动的动力学模型，发展了保结构分析方法，并用于分析轴向运动功能梯度梁横向振动的频率分布情况．研究结果表明：本文构造的数值求解算法在求解步长满足给定条件时具有良好的保结构性能，得到的前六阶频率值与复模态法得到的结果吻合较好，同时分析了微分求积法得到的结果与保结构方法和复模态法得到的结果有明显差距的原因，为微分求积法的进一步改进指明了方向，也为轴向运动功能梯度梁横向振动这类复杂动力学问题的求解提供了新途径．参　考　文　献１ＲｅｄｄｙＪＮ，ＣｈｉｎＣＤ．Ｔｈｅｒｍｏｍｅｃｈａｎｉｃａｌａｎａｌｙｓｉｓｏｆｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｃｙｌｉｎｄｅｒｓａｎｄｐｌａｔｅｓ．ＪｏｕｒｎａｌｏｆＴｈｅｒ ｍａｌＳｔｒｅｓｓｅｓ，１９９８，２１（６）：５９３～６２６２ＮａｅｂｅＭ，ＳｈｉｒｖａｎｉｍｏｇｈａｄｄａｍＫ．Ｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｍａｔｅｒｉａｌｓ：ａｒｅｖｉｅｗｏｆｆａｂｒｉｃａｔｉｏｎａｎｄｐｒｏｐｅｒｔｉｅｓ．ＡｐｐｌｉｅｄＭａｔｅｒｉａｌｓＴｏｄａｙ，２０１６，５：２２３～２４５３ＢａｒｔｌｅｔｔＮＷ，ＴｏｌｌｅｙＭＴ，ＯｖｅｒｖｅｌｄｅＪＴＢ，ｅｔａｌ．Ａ３Ｄ ｐｒｉｎｔｅｄ，ｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｓｏｆｔｒｏｂｏｔｐｏｗｅｒｅｄｂｙｃｏｍ ｂｕｓｔｉｏｎ．Ｓｃｉｅｎｃｅ，２０１５，３４９（６２４４）：１６１～１６５４ＴｕａｎＡＮ，ＡｉｅｌｌｏＭ．Ｅｎｅｒｇｙｉｎｔｅｌｌｉｇｅｎｔｂｕｉｌｄｉｎｇｓｂａｓｅｄｏｎｕｓｅｒａｃｔｉｖｉｔｙ：ａｓｕｒｖｅｙ．ＥｎｅｒｇｙａｎｄＢｕｉｌｄｉｎｇｓ，２０１３，５６：２４４～２５７５ＳａｎｋａｒＢＶ．Ａｎｅｌａｓｔｉｃｉｔｙｓｏｌｕｔｉｏｎｆｏｒｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｂｅａｍｓ．ＣｏｍｐｏｓｉｔｅｓＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ，２００１，６１（５）：６８９～６９６６ＲｅｄｄｙＪＮ．Ｍｉｃｒｏｓｔｒｕｃｔｕｒｅ ｄｅｐｅｎｄｅｎｔｃｏｕｐｌｅｓｔｒｅｓｓｔｈｅｏ ｒｉｅｓｏｆｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｂｅａｍｓ．ＪｏｕｒｎａｌｏｆｔｈｅＭｅｃｈａｎ ｉｃｓａｎｄＰｈｙｓｉｃｓｏｆＳｏｌｉｄｓ，２０１１，５９（１１）：２３８２～２３９９７ＤｉｎｇＨ，ＣｈｅｎＬＱ．Ｇａｌｅｒｋｉｎｍｅｔｈｏｄｓｆｏｒｎａｔｕｒａｌｆｒｅ ｑｕｅｎｃｉｅｓｏｆｈｉｇｈ ｓｐｅｅｄａｘｉａｌｌｙｍｏｖｉｎｇｂｅａｍｓ．ＪｏｕｒｎａｌｏｆＳｏｕｎｄａｎｄＶｉｂｒａｔｉｏｎ，２０１０，３２９（１７）：３４８４～３４９４８姚晓莎，王忠民，赵凤群．轴向运动功能梯度梁的横向振动．机械工程学报，２０１３，４９（２３）：１１７～１２２（ＹａｏＸＳ，ＷａｎｇＺＭ，ＺｈａｏＦＱ．Ｔｒａｎｓｖｅｒｓｅｖｉｂｒａｔｉｏｎｏｆａｘｉａｌｌｙｍｏｖｉｎｇｂｅａｍｍａｄｅｏｆｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｍａｔｅｒｉ ａｌｓ．ＪｏｕｒｎａｌｏｆＭｅｃｈａｎｉｃａｌＥｎｇｉｎｅｅｒｉｎｇ，２０１３，４９（２３）：１１７～１２２（ｉｎＣｈｉｎｅｓｅ））９刘金建，蔡改改，谢锋，等．轴向运动功能梯度粘弹性梁横向振动的稳定性分析．动力学与控制学报，２０１６，１４（６）：５３３～５４１（ＬｉｕＪＪ，ＣａｉＧＧ，ＸｉｅＦ，ｅｔａｌ．Ｓｔａｂｉｌｉｔｙａｎａｌｙｓｉｓｏｎｔｒａｎｓｖｅｒｓｅｖｉｂｒａｔｉｏｎｏｆａｘｉａｌｌｙｍｏｖｉｎｇｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｖｉｓｃｏｅｌａｓｔｉｃｂｅａｍｓ．ＪｏｕｒｎａｌｏｆＤｙｎａｍｉｃｓａｎｄＣｏｎｔｒｏｌ，２０１６，１４（６）：５３３～５４１（ｉｎＣｈｉ ｎｅｓｅ））１０ＢａｌｉｒｅｄｄｙＳＮ，ＰｉｔｃｈａｉｍａｎｉＪ．Ｓｔａｂｉｌｉｔｙａｎｄｄｙｎａｍｉｃｂｅｈａｖｉｏｕｒｏｆｂｉ ｄｉｒｅｃｔｉｏｎａｌｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｂｅａｍｓｕｂｊｅｃ ｔｅｄｔｏｖａｒｉａｂｌｅａｘｉａｌｌｏａｄ．ＭａｔｅｒｉａｌｓＴｏｄａｙＣｏｍｍｕｎｉｃａ ｔｉｏｎｓ，２０２２，３２：１０４０４３１１ＨｕＷ，ＷａｎｇＺ，ＺｈａｏＹ，ｅｔａｌ．Ｓｙｍｍｅｔｒｙｂｒｅａｋｉｎｇｏｆｉｎｆｉｎｉｔｅ ｄｉｍｅｎｓｉｏｎａｌｄｙｎａｍｉｃｓｙｓｔｅｍ．ＡｐｐｌｉｅｄＭａｔｈｅｍａｔｉｃｓＬｅｔｔｅｒｓ，２０２０，１０３：１０６２０７４０１Copyright ©博看网. All Rights Reserved.第６期刘涛等：轴向运动功能梯度梁横向振动问题的保结构分析１２ＨｕＷＰ，ＤｅｎｇＺＣ，ＨａｎＳＭ，ｅｔａｌ．Ｇｅｎｅｒａｌｉｚｅｄｍｕｌｔｉ ｓｙｍｐｌｅｃｔｉｃｉｎｔｅｇｒａｔｏｒｓｆｏｒａｃｌａｓｓｏｆｈａｍｉｌｔｏｎｉａｎｎｏｎｌｉｎｅａｒｗａｖｅＰＤＥｓ．ＪｏｕｒｎａｌｏｆＣｏｍｐｕｔａｔｉｏｎａｌＰｈｙｓｉｃｓ，２０１３，２３５：３９４～４０６１３宋明哲，邓子辰，赵云平，等．含弱阻尼空间结构的耦合动力学保结构分析．动力学与控制学报，２０１９，１７（５）：４１９～４２４（ＳｏｎｇＭＺ，ＤｅｎｇＺＣ，ＺｈａｏＹＰ，ｅｔａｌ．Ｃｏｕｐｌｉｎｇｄｙｎａｍｉｃｓｔｒｕｃｔｕｒｅ ｐｅｒｓｅｖｅｒｉｎｇａｎａｌｙｓｉｓｏｆｓｐａｔｉａｌｓｔｒｕｃｔｕｒｅｗｉｔｈｗｅａｋｄａｍｐｉｎｇ．ＪｏｕｒｎａｌｏｆＤｙｎａｍｉｃｓａｎｄＣｏｎｔｒｏｌ，２０１９，１７（５）：４１９～４２４（ｉｎＣｈｉｎｅｓｅ））１４ＨｕＷ，ＨｕａｉＹ，ＸｕＭ，ｅｔａｌ．Ｍｅｃｈａｎｏｅｌｅｃｔｒｉｃａｌｆｌｅｘｉｂｌｅｈｕｂ ｂｅａｍｍｏｄｅｌｏｆｌｏｎｉｃ ｔｙｐｅｓｏｌｖｅｎｔ ｆｒｅｅｎａｎｏｆｌｕｉｄｓ．ＭｅｃｈａｎｉｃａｌＳｙｓｔｅｍｓａｎｄＳｉｇｎａｌＰｒｏｃｅｓｓｉｎｇ，２０２１，１５９：１０７８３３１５ＨｕＷ，ＸｕＭ，ＳｏｎｇＪ，ｅｔａｌ．Ｃｏｕｐｌｉｎｇｄｙｎａｍｉｃｂｅｈａｖｉｏｒｓｏｆｆｌｅｘｉｂｌｅｓｔｒｅｔｃｈｉｎｇｈｕｂ ｂｅａｍｓｙｓｔｅｍ．ＭｅｃｈａｎｉｃａｌＳｙｓｔｅｍｓａｎｄＳｉｇｎａｌＰｒｏｃｅｓｓｉｎｇ，２０２１，１５１：１０７３８９１６ＨｕＷ，ＸｕＭ，ＺｈａｎｇＦ，ｅｔａｌ．Ｄｙｎａｍｉｃａｎａｌｙｓｉｓｏｎｆｌｅｘｉｂｌｅｈｕｂ ｂｅａｍｗｉｔｈｓｔｅｐ ｖａｒｉａｂｌｅｃｒｏｓｓ ｓｅｃｔｉｏｎ．Ｍｅｃｈａｎｉ ｃａｌＳｙｓｔｅｍｓａｎｄＳｉｇｎａｌＰｒｏｃｅｓｓｉｎｇ，２０２２，１８０：１０９４２３１７ＨｕＷＰ，ＤｅｎｇＺＣ，ＷａｎｇＢ，ｅｔａｌ．Ｃｈａｏｓｉｎａｎｅｍｂｅｄｄｅｄｓｉｎｇｌｅ ｗａｌｌｅｄｃａｒｂｏｎｎａｎｏｔｕｂｅ．ＮｏｎｌｉｎｅａｒＤｙｎａｍｉｃｓ，２０１３，７２（１ ２）：３８９～３９８ＳＴＲＵＣＴＵＲＥ ＰＲＥＳＥＲＶＩＮＧＡＮＡＬＹＳＩＳＯＮＴＲＡＮＳＶＥＲＳＥＶＩＢＲＡＴＩＯＮＯＦＦＵＮＣＴＩＯＮＡＬＬＹＧＲＡＤＥＤＢＥＡＭＷＩＴＨＡＮＡＸＩＡＬＶＥＬＯＣＩＴＹＬｉｕＴａｏ１　ＺｈｏｕＹａｎｇｘｉｎ２　ＨｕＷｅｉｐｅｎｇ２（１．ＹｕｌｉｎＣｉｔｙＩｎｖｅｓｔｍｅｎｔＣｏｎｓｔｒｕｃｔｉｏｎＤｅｖｅｌｏｐｍｅｎｔＣｏ．，Ｌｔｄ．，Ｙｕｌｉｎ　７１９０００，Ｃｈｉｎａ）（２．ＳｃｈｏｏｌｏｆＣｉｖｉｌＥｎｇｉｎｅｅｒｉｎｇａｎｄＡｒｃｈｉｔｅｃｔｕｒｅ，Ｘｉ’ａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ，Ｘｉ’ａｎ　７１００４８，Ｃｈｉｎａ）Ａｂｓｔｒａｃｔ　Ｔｈｅａｘｉａｌｖｅｌｏｃｉｔｙａｎｄｔｈｅｍａｔｅｒｉａｌ’ｓｈｅｔｅｒｏｇｅｎｅｉｔｙｉｎｔｒｏｄｕｃｅｔｈｅｇｒｅａｔｃｈａｌｌｅｎｇｅｏｎｔｈｅｖｉｂｒａｔｉｏｎａｎａｌｙｓｉｓｏｆｔｈｅｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｂｅａｍｗｉｔｈａｎａｘｉａｌｖｅｌｏｃｉｔｙ．Ｉｎｔｈｉｓｗｏｒｋ，ｔｈｅｄｙｎａｍｉｃｍｏｄｅｌｏｆｔｈｅｔｒａｎｓｖｅｒｓｅｖｉｂｒａｔｉｏｎｏｆｔｈｅｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｂｅａｍｗｉｔｈａｎａｘｉａｌｖｅｌｏｃｉｔｙｉｓｒｅｖｉｅｗｅｄｉｎｂｒｉｅｆｆｉｒｓｔｌｙ．Ｂａｓｅｄｏｎｔｈｅｄｙｎａｍｉｃｓｙｍｍｅｔｒｙｂｒｅａｋｉｎｇｔｈｅｏｒｙａｎｄｔｈｅｇｅｎｅｒａｌｉｚｅｄｍｕｌｔｉ ｓｙｍｐｌｅｃｔｉｃｍｅｔｈｏｄｆｏｒｔｈｅｉｎｆｉｎｉｔｅ ｄｉｍｅｎｓｉｏｎａｌｓｙｓｔｅｍ，ａｓｔｒｕｃｔｕｒｅ ｐｒｅｓｅｒｖｉｎｇｎｕｍｅｒｉｃａｌｓｃｈｅｍｅｆｏｒｔｈｅｄｙｎａｍｉｃｍｏｄｅｌｉｓｄｅｖｅｌｏｐｅｄ．Ｉｎｔｈｅｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎ，ｔｈｅｃｒｉｔｉ ｃａｌｓｔｅｐｌｅｎｇｔｈｓａｔｉｓｆｙｉｎｇｔｈｅｇｅｎｅｒａｌｉｚｅｄｍｕｌｔｉ ｓｙｍｐｌｅｃｔｉｃｃｏｎｄｉｔｉｏｎｉｓｏｂｔａｉｎｅｄｗｉｔｈｔｈｅｇｉｖｅｎｍａｔｅｒｉａｌｐａｒａｍｅ ｔｅｒｓ．Ｔｈｅｆｉｒｓｔｓｉｘｆｒｅｑｕｅｎｃｉｅｓｏｆｔｈｅｔｒａｎｓｖｅｒｓｅｖｉｂｒａｔｉｏｎｍｏｄｅｌａｒｅｐｒｅｓｅｎｔｅｄｅｍｐｌｏｙｉｎｇｔｈｅｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄ，ｔｈｅｃｏｍｐｌｅｘｍｏｄａｌｍｅｔｈｏｄａｎｄｔｈｅｓｔｒｕｃｔｕｒｅ ｐｒｅｓｅｒｖｉｎｇｍｅｔｈｏｄｒｅｓｐｅｃｔｉｖｅｌｙ．Ｆｒｏｍｔｈｅｎｕｍｅｒｉｃａｌｒｅ ｓｕｌｔｓ，ｉｔｃａｎｂｅｆｏｕｎｄｔｈａｔｔｈｅｆｉｒｓｔｓｉｘｆｒｅｑｕｅｎｃｉｅｓｏｂｔａｉｎｅｄｂｙｕｓｉｎｇｔｈｅｓｔｒｕｃｔｕｒｅ ｐｒｅｓｅｒｖｉｎｇｍｅｔｈｏｄａｒｅｈｉｇｈｌｙｃｏｎｓｉｓｔｅｎｔｗｉｔｈｔｈｏｓｅｏｂｔａｉｎｅｄｂｙｕｓｉｎｇｔｈｅｃｏｍｐｌｅｘｍｏｄａｌｍｅｔｈｏｄ．Ｔｏｉｍｐｒｏｖｅｔｈｅｐｒｅｃｉｓｉｏｎｏｆｔｈｅｄｉｆｆｅｒｅｎｔｉａｌｑｕａｄｒａｔｕｒｅｍｅｔｈｏｄ，ｔｈｅｍａｉｎｆａｃｔｏｒｓｒｅｓｕｌｔｉｎｇｉｎｔｈｅｅｒｒｏｒａｒｅｉｎｖｅｓｔｉｇａｔｅｄ．Ｔｈｅｍａｉｎｃｏｎｔｒｉｂｕｔｉｏｎｏｆｔｈｉｓｗｏｒｋｉｓｐｒｏｐｏｓｉｎｇａｎｅｗａｐｐｒｏａｃｈｔｏａｎａｌｙｚｅｔｈｅｃｏｍｐｌｅｘｄｙｎａｍｉｃｐｒｏｂｌｅｍｌｉｋｅｔｈｅｔｒａｎｓｖｅｒｓｅｖｉｂｒａｔｉｏｎｏｆｔｈｅｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｂｅａｍｗｉｔｈａｎａｘｉａｌｖｅｌｏｃｉｔｙｃｏｎｓｉｄｅｒｅｄｉｎｔｈｉｓｐａｐｅｒ．Ｋｅｙｗｏｒｄｓ　ｓｔｒｕｃｔｕｒｅ ｐｒｅｓｅｒｖｉｎｇ，　ｆｕｎｃｔｉｏｎａｌｌｙｇｒａｄｅｄｂｅａｍｗｉｔｈａｎａｘｉａｌｖｅｌｏｃｉｔｙ，　ｓｙｍｍｅｔｒｙｂｒｅａｋｉｎｇ，　ｇｅｎｅｒａｌｉｚｅｄｍｕｌｔｉ ｓｙｍｐｌｅｃｔｉｃ，　ｔｒａｎｓｖｅｒｓｅｖｉｂｒａｔｉｏｎＲｅｃｅｉｖｅｄ２０Ａｕｇｕｓｔ２０２２，ｒｅｖｉｓｅｄ２８Ｓｅｐｔｅｍｂｅｒ２０２２．ＴｈｅｐｒｏｊｅｃｔｓｕｐｐｏｒｔｅｄｂｙｔｈｅＮａｔｉｏｎａｌＮａｔｕｒａｌＳｃｉｅｎｃｅＦｏｕｎｄａｔｉｏｎｏｆＣｈｉｎａ（１２１７２２８１，１１９７２２８４），ＦｏｕｎｄａｔｉｏｎＳｔｒｅｎｇｔｈｅｎｉｎｇＰｒｏｇｒａｍｍｅＴｅｃｈｎｉｃａｌＡｒｅａＦｕｎｄ（２０２１ ＪＣＪＱ ＪＪ ０５６５），ｔｈｅＦｕｎｄｏｆｔｈｅＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙＩｎｎｏｖａｔｉｏｎＴｅａｍｏｆＳｈａａｎｘｉ（２０２２ＴＤ ６１）ａｎｄＦｕｎｄｏｆｔｈｅＹｏｕｔｈＩｎｎｏ ｖａｔｉｏｎＴｅａｍｏｆＳｈａａｎｘｉＵｎｉｖｅｒｓｉｔｉｅｓ ＣｏｒｒｅｓｐｏｎｄｉｎｇａｕｔｈｏｒＥ ｍａｉｌ：ｗｐｈｕ＠ｎｗｐｕ．ｅｄｕ．ｃｎ５０１Copyright ©博看网. All Rights Reserved.。

编号：曲轴系动力学计算分析指南（第I 版）目录目录 (1)1 FE模型 (2)1.1 FE模型建立 (2)1.1.1 曲轴系坐标系的定义 (2)1.1.2 曲轴系FE模型建立 (2)1.1.3 简易主轴承壁FE模型建立 (3)1.1.4边界条件 (3)1.2 FE模型压缩 (4)2 EXCITE模型建立 (4)2.1 FE子结构模型 (4)2.2 EXCITE模型建立 (5)2.2.1连接单元参数的确定 (6)2.2.2 曲轴参考点的定义 (6)3 EXCITE仿真计算 (8)3.1 轴承负荷计算 (8)3.2 轴承性能计算 (8)3.2.1 最大油膜压力 (8)3.2.2 最小油膜厚度 (9)3.2.3 轴心轨迹 (10)4 动应力计算 (10)5 曲轴系平衡计算 (12)6 曲轴扭振及其阶次分析 (13)7 疲劳分析 (14)1 FE模型1.1 FE模型建立几何模型的完整是FE模型建立的前提条件，一个完整的曲轴系几何模型主要由曲轴，与其相连的飞轮和正时齿轮以及皮带轮组成，其中各个零部件之间由非线性连接体连接。

FE模型可以用前处理软件PATRAN、HYPERMESH等来建立。

采用手动划分网格的方法，各部件间通过合并接触面上分布一致的节点来构成一个整体。

1.1.1 曲轴系坐标系的定义整体坐标系采用右手法则的直角坐标系，如图1，坐标系的中心在曲轴第三段主轴颈的中心，X轴为曲轴的轴线方向，Y轴的方向为曲轴的侧向，Z轴与气缸同向，同时要求曲轴的第一拐朝上放置.为了利用A VL-EXCITE软件进行曲轴系的动力学计算，需要曲轴系FE模型以及一个简易的主轴承壁FE模型。

1.1.2 曲轴系FE模型建立曲轴系FE模型采用的六面体网格如图1～图3所示：图1 曲轴系的有限元模型图2 单拐的有限元模型 图3 主轴颈和拐处圆角的有限元模型在进行曲轴系的动态分析时，为了保证在圆角处有足够的网格精度，通常需要6层以上的单元。

杆件的挠曲与弯曲分析杆件是工程中常见的结构元件，通常用于承受压力或拉力。

在实际应用中，杆件常常会受到挠曲与弯曲的力作用，因此对杆件的挠曲与弯曲进行分析十分重要。

本文将分析杆件在挠曲与弯曲力作用下的应力与变形，并基于此进行讨论。

1. 挠曲力学分析挠曲是指杆件在纵轴方向上受到的弯曲力作用，导致杆件产生曲率。

挠曲过程中，杆件内部产生应力分布。

挠曲现象在悬臂梁、梁柱等结构中常常出现。

杆件的挠曲行为可以通过弯曲矩和截面惯性矩来描述。

弯曲矩代表了杆件上点的弯曲力矩大小，而截面惯性矩则反映了杆件截面形状对挠曲的抵抗能力。

挠曲时，杆件上任意一点的曲率与弯曲矩和截面惯性矩之间存在一定的关系。

2. 挠曲的应力分析杆件在挠曲过程中产生应力，这些应力主要分为压应力和拉应力。

在杆件挠曲点的外侧产生压应力，而在内侧产生拉应力。

应力分布沿着截面的纵轴方向呈现三角形分布。

挠曲引起的应力主要由弯曲应力和剪切应力组成。

弯曲应力与弯曲矩成正比，而剪切应力则是由弯曲力矩在杆件横截面上产生的。

在挠曲过程中，应力的最大值通常出现在杆件截面的最外侧纤维上，该处被称为受压纤维。

对于圆形截面的杆件，受压纤维位于截面的最上方。

3. 弯曲力学分析弯曲是指杆件受到横向力作用而产生的曲率。

与挠曲不同，弯曲是杆件整体弯曲，而不是仅在纵轴方向上发生曲率变化。

在弯曲过程中，杆件的上部受到压力，而下部受到拉力。

类似于挠曲，弯曲引起的应力也主要分为弯曲应力和剪切应力。

弯曲应力与杆件受到的弯矩成正比，而剪切应力则是由弯矩在杆件截面上产生的。

弯曲应力的最大值通常出现在杆件截面的最外侧纤维上，受压纤维和受拉纤维分别位于截面的上下方。

4. 挠曲与弯曲的变形分析挠曲和弯曲力作用下，杆件会发生不可忽视的变形。

挠曲引起的变形主要是杆件曲率的改变，而弯曲引起的变形则是杆件整体形状的改变。

杆件挠曲的变形会导致杆件长度的增加，同时产生横向位移。

而杆件弯曲的变形主要表现为杆件整体形状的改变，类似于一个弯曲的弧线。

汽车发动机曲轴扭转振动分析及控制社会经济在进行着快速的发展中，人们对于汽车的使用量也在逐渐的增加，我国对于汽车建设中是要求也越加严苛。

在汽车公司进行汽车设计的过程中，对于发动机及行驶中的稳定程度越加重视。

汽车发动机曲轴扭转振动是汽车公司在对于发动机研究中的热点课题。

本为对于发动机的曲轴扭转技术进行较为全面的分析。

标签：曲轴系；扭转振动；优化设计0 前言增加对于汽车发动机的振动分析与控制，在一定程度上面可以将汽车的内部结构进行优化，增加发动机的使用时间与汽车行驶过程中的稳定性能。

曲轴扭转是发动机在工作过程中的主要部件，性能的好坏将直接对于汽车的整体性能进行影响。

本文主要对于汽车中的曲轴扭转振动进行分析研究，这项研究是十分具有实际意义的。

1 汽车发动机曲轴扭转振动系统理论分析1.1 ADAMS多刚体动力学理论ADAMS动力学理论主要使用坐标方程式进行汽车在行驶中的发动机系统的分析。

在ADAMS动力学理论中，将动力系统内的关性参考系中的坐标与方位坐标进行标注，并使用相对应的数学方程式进行多余坐标的约束，进而将已经标注的坐标进行变量。

在对于动力学的分析过程中，使用数学方程式可以将计算的效率进行大幅度提升。

1.2 ADAMS多柔体动力学理论在进行汽车生产建设中，在机械系统中已经广泛使用柔性材料，是生产设备运行中速度较快，但是运行的精度也在不断的提升，设备内的动力学性能变得更加繁琐。

刚性研究体系已经不能满足对于动力学的研究，因此柔体动力学理论就在这种情况下产生。

这种研究体系一般情况下是以刚性动力学体系作为参照依据，在对于柔体的研究中进常采用不同的处理形式。

在一定程度上面刚性与柔性的个、动力学体系进行共同使用，可以对系统中的动力学进行更加全面的认识[1]。

2 曲轴动力学研究模型2.1 三维几何模型三维几何模型可以将曲轴系统的中每个零件间的关系进行清晰的展示。

按照零件的规格与参数，利用相对应的三维软件就可以建立相对应的三维几何模型。

发动机曲轴系扭振分析许涛;阮仁宇【摘要】采用EXCITE软件,对某汽油机的曲轴系进行了扭振分析,首先进行了减震皮带轮的频率选定,然后针对选定频率的曲轴系进行扭振分析,分别计算了转速波动、耗散功以及扭振角.计算结果显示:以上指标均在可接受的范围内,满足使用要求.【期刊名称】《汽车零部件》【年(卷),期】2016(000)007【总页数】3页(P56-58)【关键词】发动机;曲轴系;扭振分析【作者】许涛;阮仁宇【作者单位】安徽江淮汽车股份有限公司技术中心,安徽合肥230601;安徽江淮汽车股份有限公司技术中心,安徽合肥230601【正文语种】中文【中图分类】U464发动机曲轴系包括曲轴以及与曲轴相连的各运动件(如飞轮、减震皮带轮等)，曲轴系是发动机最为关键的零部件，其作用是把活塞、连杆传来的气体力转化成力矩，用来驱动汽车的传动系统、发动机的配气机构以及其他辅助装置。

曲轴系是一个弹性体，存在一系列的扭振固有频率，在发动机运行过程中，曲轴系将产生明显的扭转振动。

而这个扭转振动是发动机上危害最大的振动问题，不仅消耗发动机的有用功，严重时甚至会导致曲轴断裂。

因此，在设计开发时必须对曲轴的扭振情况进行分析。

借助于EXCITE软件，以某四缸汽油机为研究对象，对其曲轴系进行扭振分析，优选皮带轮的固有频率，评价曲轴系的扭振特性，为曲轴系的设计提供理论参考。

1.1 主要参数发动机的主要参数如表1所示。

计算用的爆发压力(额定工况)来源于热力学的计算结果，如图1所示。

1.2 建立计算模型在EXCITE中搭建曲轴系的等效计算模型，如图2所示，定义好各个部件的设计参数、润滑油特性、仿真控制参数，接下来就可以进行计算。

其中曲轴使用Shaftmodeler模块和Autoshaft模块完成，连杆使用Conrod Modeler完成。

2.1 皮带轮优选由于曲轴是一个柔性轴，以及爆发压力的变动幅度很大，曲轴系在运转过程中会产生扭转振动，这对曲轴的寿命和发动机的工作稳定性都有影响。