所在位置: > 合乐888在线客服 >

合乐888在线客服
联系方式
电话:0319 7588019
传真:0319 7588019
邮编:055151
地址:河北省任县 邢家湾镇西黄庄工业区
使用 COMSOL® 软件预测转子系统的稳定性
发布时间:2024-03-06 点击: 次 编辑:admin

  在使用旋转零件时,稳定性分析至关重要,因为不稳定性会导致灾难性故障。旋转系统会因为圆盘的不对称惯性、轴的不对称刚度或轴承引起的交叉耦合效应而导致不稳定的响应。从设计者的角度来看,重要的是要确保潜在的不稳定模式不在机器的工作范围内。本文让我们来探索如何使用 COMSOL Multiphysics® 软件预测转子系统的不稳定性。

  静稳定性是基于系统从平衡状态受到干扰时的初始响应,可进一步分为以下几种情况:

  动态稳定性是基于系统从其平衡位置受到干扰后的运动时间历程,可进一步分为以下几种情况:

  现在我们已经熟悉了基本术语,接下来我们来了解稳定性的数学定义和系统变得不稳定的条件。

  如果所有的系数都是正的,那么控制方程的解也是有界限的并且表现良好。然而,在某些情况下,系数变为负值,解变得无界限。为了理解这个问题,让我们来看看弹簧-质量-阻尼器系统在各种情况下的响应,如下图所示。

  从系统的预测响应中,我们可以说,当系统的刚度为负值时,系统是静态不稳定的;如果阻尼为负值,则系统是动态不稳定的。我们也可以从上图中发现,系统静态稳定才能实现动态稳定。然而,情况并非总是如此,尤其是当陀螺效应出现时。一个常见的例子是自行车,它是静态不稳定,但在一定速度范围内动态稳定鞍形弹性垫圈。

  线性稳定性分析可以用来预测系统的不稳定性。线性稳定性分析框架包括以下步骤:

  转子动力学问题不同于传统的振动分析。附加项的出现是因为运动方程中存在附加的加速力项,这取决于所选的物理性质。

  在转子动力学模块有两种方法可以解决这个问题。实心转子在旋转坐标系中形成,而梁转子在固定坐标系中形成。基于所选的物理场,运动方程中出现的附加项和效应如下所示:

  实心转子在旋转坐标系中建模,因此在该坐标系中不会观察到转子的物理旋转。在动量平衡方程中,旋转的影响被解释为坐标系加速力。因为这个运动方程包含附加的项,所以旋转坐标系中的运动方程是:

  其中,\textbf M是质量矩阵;\textbf C_c是科里奥利力效应的贡献(取决于转子的速度);\textbf C是外部(固定)和内部(旋转)阻尼的贡献;\textbf K_c是旋转坐标系中外部阻尼的贡献(取决于转子的速度);\textbf K是由内部刚度、外部刚度、离心软化和应力硬化引起的刚度矩阵;\textbf F是旋转坐标系中的力。

  对于实心转子,没有陀螺项;相反,会得到离心力(与刚度矩阵的对称部分结合在一起)、科里奥利力(作为阻尼矩阵的反对称部分出现)、欧拉力(作为反对称刚度矩阵出现),以及作为刚度矩阵反对称项中的外部阻尼的贡献。

  其中,\textbf M是质量矩阵;\textbf C_g是陀螺效应的贡献(取决于转子的速度);\textbf C是外部(静止)和内部(旋转)阻尼的贡献;\textbf K_c是固定坐标系中内部阻尼的贡献(取决于转子的速度);\textbf K 是刚度矩阵;\textbf F是力矢量,包含力和力矩(梁上的弯矩和扭矩)。

  在梁转子中,陀螺项构成了阻尼矩阵的反对称部分,内部阻尼在刚度矩阵中表现为反对称项。

  考虑一个与全局x-轴对齐的轴。转子上的局部y和z方向与全局y和z方向相同方向。

  对于转子动力学问题加湿器,横向运动方程是耦合的;因此,它们比轴向和扭转振动的方程更复杂。在本节中,我们将讨论直接刚度和交叉耦合刚度对转子横向运动的影响合乐888在线客服。考虑一个对称的转子,两边都有对称的轴承。轴承包含直接和交叉耦合项。

  让我们讨论每个项对转子系统横向运动的单独影响。根据输入,有两种类型的旋转模式:向前和向后。在向前旋转的情况下,转子的轨道运动和旋转方向是相同的,而向后旋转则相反。

  直接刚度项(\textbf K_{yy}, \textbf K_{zz})影响偏转矢量方向的力。如果系数为正,那么力的方向与偏转方向相反。如果系数值为负,则力与挠度方向相同。

  正直接阻尼项(\textbf C_{yy}, \textbf C_{zz})产生一个垂直于偏转矢量的力,方向与涡旋速度方向相反。如果数值为负,力的方向是旋转速度方向。

  如果\textbf K_{yz}, \textbf K_{zy}有相反的符号,则它们会产生一个垂直于偏转向量的力。方向取决于系数的符号。

  简而言之,产生径向力的项(直接刚度和交叉耦合阻尼)不会导致系统不稳定,而产生与涡旋轨道相切的力的项会导致系统不稳定。该系数的影响也根据涡动模式而变化。例如,如果\textbf K_{yz}(+ve), \textbf K_{zy}(-ve),它会产生一个与涡动轨道相切的力,这个力会使前向涡动不稳定,使后向涡动稳定。基于旋转模式的效果总结如下:

  作为一名设计师,几乎不可能创造出完美的转子,因为模型中总是存在一些不对称。不对称可能是由于质量分布对于无阻尼的转子,也可能是由于形状造成的刚度。此外,转子可能会受到交叉耦合力、内部摩擦、转子摩擦和汽流涡旋的影响。这些条件的存在会影响转子系统的稳定性。我们来看几个这样的案例…

  在轴承中,交叉耦合力是不稳定的关键原因。交叉耦合力导致系统阻尼的快速损失,并导致转子次同步振动。在某些情况下,这可能导致零阻尼或负阻尼,从而使系统不稳定。

  通常,在液体动压轴承中,交叉耦合力在转子中起到负阻尼的作用。当接近临界速度时,这些力会导致不受控制的振动,甚至导致轴承故障。

  COMSOL案例库中的模型讨论了交叉耦合力引起的稳定性。该模型包括一个涡轮增压器转子,由两个轴承支撑:一个靠近压缩机,另一个靠近涡轮机,使压缩机和涡轮机都悬挂在轴上。

  针对两种情况进行分析,一种是没有交叉耦合力,另一种是有交叉耦合力。在 COMSOL Multiphysics 中,有三种不同的方式来包含轴承的交叉耦合效应:

  在上面的图中,在没有交叉耦合刚度的情况下,对数衰减大部分是正的,表明固有模式是稳定的。只有在高转速下,才会有一些模式变得不稳定。在存在交叉耦合刚度的情况下,许多模态具有负对数衰减,即使在低转子速度下也是如此。

  结果表明,交叉耦合刚度的存在使振动模式不稳定额定转速,因此在这些速度下运行涡轮增压器转子是危险的。

  机械耗散的元素会导致不稳定,这听起来很奇怪,但在旋转组件的情况下,这是真的去钝化。由于旋转阻尼附着在物体上,所以它位于旋转坐标系中啮合。在静止坐标系中进行坐标变换后,它不再仅仅依赖于速度。这种转变导致耗散项和循环项,它们取决于转子速度和位移。循环项是一个不稳定的原因。

  考虑一个 Laval-Jeffcott 转子,它由一个无质量的轴和一个由两端具有对称阻尼和刚度的两个相同轴承支撑的圆盘组成。现在,如果旋转阻尼效应也包含在系统中,循环项出现在控制方程中会导致不稳定。

  为了进一步研究这种影响,我们来画一个系统的力图。从力图中可以清楚地看出,循环项在向前旋转的运动方向上产生力。这意味着循环项是在给系统增加能量,由于外部阻尼而产生的耗散项是在不断地从系统中移除能量。

  从力图中,我们可以说条件\Omega = \omega \frac{C_e+C_i}{C_i}在异步涡动条件下触发不稳定性。当旋转速度和旋转频率不同时缓冲接合过程,就会发生异步旋转。转子的内部阻尼是由材料的滞后阻尼或冷缩配合零件界面摩擦产生的库仑阻尼产生的。

  上图比较了有内部阻尼和没有内部阻尼的转子模型。交叉耦合项产生的力本质上是非保守的。内部阻尼作为交叉耦合项对刚度矩阵有贡献,它取决于速度,所以对数衰减不是常数,随速度而变化。对数衰减是针对两种初始模式(向前和向后模式)绘制的。在零自旋时,两种模式具有相同的对数衰减。随着速度的增加,后向涡动的对数衰减增加,前向涡动的对数衰减减小,最终变为负值。这证实了交叉耦合项使后向涡动稳定,并使前向涡动不稳定。

  如果我们进行频域分析,结果与没有阻尼的情况相同,并且涡动是圆形的。这意味着模型中只有向前的涡动。

  由于旋转阻尼引起的交叉耦合项仅影响阻尼特性,而不改变系统的响应。现在,如果我们看坎贝尔(Campbell)图,自然频率不再是一个恒定值,或者如果有内部阻尼,它会随着速度而增加。在内部阻尼的情况下,或者如果我们想用坎贝尔图找到不稳定性阈值,这种效应会增加阻尼临界速度。然后,我们可以画一条线\omega=\Omega(\frac{C_i}{C_i+C_e})\,线与自然频率的交点给出了阈值。

  当沿圆周发生渗漏导致径向位移时,在轴流式机器(如燃气轮机)中会产生蒸汽涡动。这会产生一个垂直于位移的合力。激发出来的力在某些工况下产生不稳定性,在另一些工况下产生稳定性。蒸汽旋转力以反对称的方式增加了刚度项,就像内部阻尼一样,但不依赖于旋转速度。

  为了更容易理解这种影响,我们来讨论在蒸汽涡动作用下滚道,具有外部阻尼和刚度的Laval-Jeffcott转子的力图。

  从力图中,我们可以看到什么时候达到稳定阈值K_{sw}=C_e\cdot\omega_n,其中w_{n}是自然频率。这意味着可以通过增加外部阻尼或固有频率来提高稳定性阈值。

  让考虑拉 Laval-Jeffcott 转子轴流叶轮,由于反对称项,系统最初是不稳定的。如果我们在系统中引入外部阻尼,对数衰减将趋近于零,最终变为正值。通过增加外部阻尼的值,可以增加系统不稳定阈值的频率值。

  本文我们讨论了稳定性的类型、阻尼和刚度矩阵系数对系统响应的影响、实心转子和梁转子物理场的控制方程、使用特征频率分析预测系统不稳定性的方法以及导致转子系统不稳定性的一些应用。

Copyright 2017 合乐HL8客服 All Rights Reserved