内燃机曲轴系统的扭转振动——XiaoPeng内因:轴系本身就有自由扭振特性(惯性、弹性)外因:周期性变化的激励载荷作用在轴系上(扭振的能量来源)曲轴系统(轴系)扭振的产生原因轴系扭振的危害使曲轴间的夹角产生变化,破坏曲轴的原有平衡状态使机体的振动和噪声显著增大导致配气定时和喷油定时失去最佳状态,使内燃机工作性能变坏使传动齿轮间的撞击、摩擦加剧扭振附加应力的增加,有可能使曲轴及其传动齿轮断裂轴系扭转振动的计算①建立物理模型(将复杂的轴系简化换算成振动特性与之相同的当量系统②计算当量系统的自振特性(固有频率、相应频率下的振型)③对作用在各曲拐上的激励载荷进行简谐分析④进行轴系的强迫振动计算,求出共振时的实际振幅与各轴段的扭转振动附加应力⑤针对上述计算结果,全面评定整个轴系工作是否可靠轴系的当量换算原则:振动特性相同惯量较大且较集中的部件惯量较小且较分散的部件阻尼非弹性的惯量元件无惯量的弹性元件弹性元件的轴段阻尼和惯性元件的质量阻尼激励载荷只作用在惯性元件上轴系的当量系统图单位曲柄转动惯量的计算单位曲柄=主轴颈+曲柄销+曲柄臂主轴颈(mainjournal)对曲轴中心线的转动惯量)(3244mdDLI曲柄销(crankpin)对曲轴中心线的转动惯量2p'ppRmII——曲柄销对自身中心线的转动换量'pI曲柄臂(crankweb)和平衡重(blanceweight)对曲轴中心线的转动惯量方法:分割后,求各部分的转动惯量,再求和对应于圆心角的圆弧带的转动惯量)-(2360414'iiiiiRRLIi综上,单位曲柄(crank)的转动惯量为整个曲柄臂的转动惯量用同样的方法可求得平衡重的转动惯量bwpmcIIIII22上述转动惯量可在三维CAD软件中求得niiiiiRRLI1414w)(2360活塞、连杆当量转动惯量的换算往复运动质量()的运动动能原则:运动动能不变c1pjmmm222j2jK)2sin2sin(ω2121RmvmE)821(ω2121222j20KKmRmdEE曲柄转动一周,往复运动质量的平均动能据动能相等则的值很小,可略去不计,因此,往复运动质量的当量转动惯量近似表示为设往复运动质量的当量转动惯量为,则其动能为rcI2rcrc21IErcKmEE)821(22jRmIrc2j2j)21(21RmRmIrc82连杆组换算在曲柄销中心的作旋转运动的质量对曲轴中心线的转动惯量为所以,往复运动质量的对曲轴中心线的转动惯量,相当于将质量的1/2加在曲柄销中心上所产生的转动惯量2c22RmIc综上,内燃机单缸的当量转动惯量为22)22(crcbwpmcrccsIIIIIIIIII齿轮传动系统的当量转动惯量的换算通常把系统简化成仍然按原主动轴的转速(即曲轴转速)回转的系统原则:换算前后,动能不变22''221I21'''221I所以,换算后的当量转动惯量为:令传动比则同理:21'''222''22121II''22122'''2II12i''22'''2IiI'323IiI齿轮传动系统的柔度的换算原则:弹性位能(U)不变2232323)(121eU2'23'23'23)(121eU由得两齿轮在啮合处的周向弧长位移相等'2323UU2'2323'2323)(ee231'23222DDinnDD12112'23232'2323iee所以,简化成单枝系统后,柔度是原来的1/i2倍弹性参数的当量转换通常把不同直径的轴段,转换为一标准直径的轴段,但对其长度做出相应的变化原则:换算前后,柔度不变转换前转换后其中,截面惯性矩:由得:同理,可得:411111321DGLGJLe42'12'1'1321DGLGJLe4321dJ'11ee)(41421'1DDLL)(43423'3DDLL换算后,整根轴的当量长度为:'32'1'LLLL单位曲柄的柔度计算曲柄形状复杂且形式多样,既受弯曲又受扭转,又受过渡部分及轴承间隙的影响等目前均采用半经验公式这些公式只有对与之类型相同的曲轴才适用轴系的自由振动内燃机的阻尼比较小,在计算自振频率时,一般都忽略阻尼)()(111kkkkkkkkKKII—当量转动惯量K—扭转刚度—振动扭转角据达达朗贝尔原理,得:依次得各质量的运动微分方程:0)()(111kkkkkkkkKKI0)(0)()(0)()(0)(111113222112221111nnnnkkkkkkkkKIKKIKKIKI短阵形式:0KIknIIII21Tn21111122-111111KKKKKKKKKKKKKnnkkkTn210KIk轴系作自由振动,时各质量作简谐振动,所以,上式的特解为:其中代入原式,得:tAsin02AIAK0KIk解此方程,可得n个特征值及其对应的特征矩阵[A]矩阵[A]的第i列矢量{A}i就是轴系振动的第i阶固有圆频率Ωi的振形矢量2iTnAAAA21轴系自由扭转振动振形图振形图:各质量在每阶固有圆频率Ωi下的相对振幅相对振幅:将振形矢量{A}i的第一个元素进行归1化,但不改变各质量间的相对振幅比例关系不同的自振频率有不同的振形图工程上只计算一、二、三节点频率节点处扭转应力最大根据自由振动的计算结果,可以定性地了解内燃机的扭振特性轴系扭转振动的激励载荷①缸内燃气压力变化产生的激振力矩②运动部件的惯性载荷引起的激振力矩③部件吸收功率不均匀而产生的激振力矩(次要)(主要)①缸内燃气压力变化引起的干扰力矩得力矩:)(sincos142gPDT气体力沿连杆传递到曲柄销时,将它分解成垂直于曲柄的切向力T和沿曲柄半径的径向力Z,其中RPDRTMgg)(sincos142将干扰力矩按傅里叶级数展开:)sin(3sinb2sinbsinbcos3acos2acosa'''1'03213210MMMMg令:t—干扰力矩的圆频率)sin('''1'0tMMMg'其中:——v’简谐力矩的初相位对于二冲程内燃机,有对于四冲程内燃机,有2)sin('''1'0tMMMg)2sin('''1'0tMMMg写成统一行式:)sin(0tMMMg对于二冲程内燃机,对于四冲程内燃机,,3,2,1,3,25,2,23,1,21)sin(0tMMMg其中,0M使轴系以角速度ω做匀速旋转运动使轴系以不同频率νω做简谐振动(激发轴系扭转振动的内在根源))sin(tM高谐次的简谐力矩的振幅很小一般只考虑到ν=12②曲柄连杆机构惯性力引起的干扰力矩往复惯性力)2coscos(ω--2jjjRmamP往复惯性力引起的干扰力矩比较小实际中,往往只考虑燃气压力引起的干扰力矩往复惯力性力矩)4sin43sin432sin21sin4(ω)2sin2(sin)(sincos1222jRmRPRPRTMjjj干扰力矩的相位图对于单列多缸内燃机,各缸作用在曲轴上的干扰力矩的大小是相同的但各干扰力矩间有一定的相位差——第i缸与第1缸的发火间隔角——第i缸的ν次干扰力矩——第i缸的ν次干扰力矩与第一缸的ν次干扰力矩相位差因此,知道发火顺序,即可画出各缸的ν次干扰力矩的相位图干扰力矩相位图(例)四冲程六缸内燃机发火顺序1-5-3-6-2-4可以看出:简谐次数ν相差曲柄可见数的任意整数倍(即qm)时,干扰力矩矢量图相同ν为曲柄可见数的整数倍(即ν=qm)时,简谐力矩的矢量方向相同,这些谐次的干扰力矩称为主简谐ν=3/2,9/2,15/2…时,简谐力矩的矢量在同一直线上,方向相反,这些谐次的干扰力矩称为强简谐对于主简谐、强简谐引起的共振最危险,进行扭振计算时都要计算消减扭振的措施调整频率减少输出能量装设减振器改变转动惯量改变轴段的柔度插入弹性联轴器动力减振器阻尼减振器复合式减振器