图2.9偏心轴质心加速度随仿真时间变化曲线在运动学分析中,采用的驱动约束使偏心轴按照给定的转速匀速旋转,因此偏心轴质心位移在x、y两个坐标方向上的分量随仿真时间的变化曲线均为正弦曲线,如图2.7所示,偏心轴的质心距离偏心轴旋转轴线的距离为0.961mm。这部分内容会在后面偏心轴的动力学输出部分详细讲解。......
2023-06-23
小型转子发动机的偏心轴扭转振动系统由偏心轴及与其相连的轴承、转子和相位齿轮组成。在周期性变化的燃气压力和运动质量的惯性力作用下,轴系会发生振动。振动的形态有扭转振动、弯曲振动、纵向振动和滚振,这些形态的振动不仅有各自的固有振动特性,还存在一定程度的耦合[8]。
以上提到的振动将给小型转子发动机带来的危害主要是:破坏偏心轴的平衡状态,使基体的振动和噪声急剧增大,使相位齿轮齿圈的磨损加剧,严重时甚至断裂。因此,精确仿真小型转子发动机动力学时,必须考虑到偏心轴的柔性效应。
1.模态分析的有限元原理
结构模态分析是为了获得其固有振动特性,它是由结构形状、所用的材料、边界约束条件决定的,是结构动力学分析之前的一个重要环节。分析系统的模态可以用很多种方法,常用的有有限元方法[9],其主要步骤是:
(1)将结构离散成单个的有限元个体,然后用合适的方法描述它的位移,再对单元进行动力学的特点分析。可以用多种方法建立刚度矩阵、质量矩阵、阻尼矩阵。
式中,V为单元体积;B为单元应变-位移矩阵;上标T为矩阵转置运算符;D为应变-应力矩阵;N为形状函数;ρ为结构材料的密度;ν为材料阻尼系数。
(2)质量矩阵和刚度矩阵组合,建立平衡方程,也就是整体结构动力方程:
式中,M为总质量矩阵;C为总阻尼矩阵;K为总刚度矩阵;X为节点位移列阵;f为施加于整体结构的外力。
(3)找到合适的方法得到相关的导出量,求得相应的解。
模态分析主要过程包括:①建立模型;②加载并求解;③扩展模态;④结果观察及分析。
振动模态分析可根据边界的约束情况分为自由模态分析和约束模态分析。小型转子发动机的自由模态分析非常方便和简单,它在一定程度上反映系统的固有属性。但它与实际情况不符。研究发现,在进行模态分析时可以忽略阻尼的因素。固有频率和振动形式是结构的固有特性,与是否受到外部的冲击无关。
结构的无阻尼自由振动方程可以表达为
式中,M为结构的质量矩阵;K为结构的刚度矩阵;为节点的加速度矢量矩阵;U为节点的位移矢量矩阵。
对于线性系统,简谐振动形式为
式中,φi为各阶固有频率相对应的特征向量;ωi为结构的各阶固有频率;t为振动周期。
把式(2.17)代入式(2.16)得到
得到特征方程
由特征方程可以计算得到特征值ω2i,再将特征值ω2i代入结构的无阻尼自由振动方程中,可以得到与其对应的特征向量φi。
2.偏心轴的自由模态分析
本节利用有限元软件进行子结构和模态分析,图2.23所示为偏心轴的有限元网格模型。
图2.23 偏心轴的有限元网格模型
对偏心轴进行自由模态分析,获得偏心轴的固有频率和各阶振型,其固有频率如表2.2所示,部分振型如图2.24~图2.29所示。
表2.2 偏心轴自由模态频率
图2.24 偏心轴自由模态第7阶振型
图2.25 偏心轴自由模态第8阶振型
图2.26 偏心轴自由模态第9阶振型
图2.27 偏心轴自由模态第10阶振型
图2.28 偏心轴自由模态第11阶振型
图2.29 偏心轴自由模态第12阶振型
以上仿真分析结果表明,偏心轴的前6阶模态频率为零,第7阶、第9阶、第11阶、第12阶模态呈现y 向弯曲振型,第8阶、第10阶模态呈现x向弯曲振型。自由模态的偏心轴振动主要表现为偏心轴径向的变形和扭曲。偏心轴模态以高频为主,对机体低频振动特性响应较小。
3.偏心轴的约束模态分析
约束模态分析指在进行模态分析之前对零部件间的接触部位加上相应的约束条件,然后再进行模态计算。偏心轴在小型转子发动机内做旋转运动,主要与前轴承、主轴承、后轴承接触,绕前后轴承中心旋转。求解分析得到偏心轴约束模态的固有频率如表2.3所示,部分振型如图2.30~图2.34所示。
表2.3 偏心轴约束模态固有频率
图2.30 偏心轴约束模态第2阶振型
图2.31 偏心轴约束模态第3阶振型
图2.32 偏心轴约束模态第4阶振型
图2.33 偏心轴约束模态第5阶振型
图2.34 偏心轴约束模态第6阶振型
偏心轴约束模态仿真结果表明,偏心轴的约束模态和自由模态振动基本一致,主要表现为径向的弯曲和扭曲。与自由模态不同的是,约束模态表现了个别地方单独振动的模态。约束模态的第1阶频率为0,第2阶、第3阶为中频模态,其他的都是高频模态。第2阶振型与第3阶振型相位不同,振型相似。第4阶振型与第5阶振型相位不同,振型相似。结构模态以高频为主,对机体低频振动特性响应较小。
有关小型转子发动机数值建模与仿真分析的文章
图2.9偏心轴质心加速度随仿真时间变化曲线在运动学分析中,采用的驱动约束使偏心轴按照给定的转速匀速旋转,因此偏心轴质心位移在x、y两个坐标方向上的分量随仿真时间的变化曲线均为正弦曲线,如图2.7所示,偏心轴的质心距离偏心轴旋转轴线的距离为0.961mm。这部分内容会在后面偏心轴的动力学输出部分详细讲解。......
2023-06-23
双击桌面上的SAMCEF Field图标或相应目录中的SAMCEF Field以打开程序。在主菜单中选择File→Open命令,浏览并选择模型文件crank-stress.sfield。输入文件名称crank-mode。如图9-2所示,在弹出的Solver Driver Settings对话框中的Analysis Type下拉列表框中选择Modal选项。图9-1 转换分析类型图9-2 Solver Driver Settings对话框单击OK按钮。单击Solver模块图标,进入求解模块。输入要计算的特征值数,例如输入10。图9-3 Eigen Values选项卡图9-4 显示执行状态......
2023-10-27
在工作室气体合力的作用下推动转子旋转,从而带动偏心轴旋转对外输出扭矩。从图2.38可以看出,偏心轴轴颈的受力也是周期变化的,这是由于离心惯性力在转速不变的情况下是不变的,偏心轴轴颈的受力变化情况只与气体作用力相关。图2.37偏心轴输出扭矩图2.37偏心轴输出扭矩图2.38偏心轴轴颈受力图图2.38偏心轴轴颈受力图......
2023-06-23
向心轴承在承受径向载荷Fr后,上半圈的滚动体不承受载荷,下半圈的滚动体承受载荷。图2 - 19 向心轴承中的弹性变形图2 - 20 向心轴承中的径向变位在图2-19中:ψ——各滚动体中心与最大载荷滚动体之间的夹角,ψ0=0°,,ψ2’,…对球轴承:对滚子轴承 t=1.1图2-21为轴承载荷分布图。每个滚动体载荷可以分解为两个分量,由力的平衡可得将式代入可得由上式可看出,Qmax与作用于轴承上的载荷Fr和滚动体数量z有关。......
2023-06-26
传统的方法都是通过附加质量法来计算结构在水中的模态,即湿模态,本文利用APDL语言命令流和Workbench平台耦合的方法来实现对转子部件在预应力情况下湿模态的计算。插入命令流如下求解之后,前三阶振型图如图5.4-4所示图5.4-4 湿模态前三阶振型图则前六阶固有频率如图5.4-5所示。另外由图5.4-4可以看出,在水中各阶振型对应的振幅相比于在真空中也相应降低,这是水介质阻尼作用的结果。......
2023-06-26
表2-10 向心轴承JrJr与载荷分布参数T的关系也可由图2-22中查出。再由径向游隙Gr的数值计算出载荷分布范围角,最后计算出作用于每个滚动体上的载荷数值。Gr=0 Qmax=1.43kNGr=0.025mm Qmax=1.647kNGr=0.050mm Qmax=1.833kN而球轴承寿命与最大滚动体载荷的三次方成反比,因此,轴承中径向游隙的大小对轴承寿命的影响是很大的。......
2023-06-26
【主要内容】设二元函数u=u(x,y),v=v(x,y)在点(x,y)处偏导数存在,函数z=f(u,v)在对应点(u,v)处可微,则复合函数z=f(u(x,y),v(x,y))在点(x,y)处的偏导数存在,且二元以上复合函数的偏导数也有类似的计算公式.注 在计算二元或二元以上复合函数的偏导数时,应先画出复合函数的关系图,然后按此图用复合函数求偏导数公式计算.例如,二元复合函数z=f(u(x,y),......
2023-10-27
【主要内容】1.二元隐函数偏导数的计算设二元函数z=z(x,y)由方程F(x,y,z)=0确定,则称z=z(x,y)是二元隐函数.它的偏导数可按以下步骤计算:(1)所给方程两边对x(或y)求偏导数(此时应注意z是关于x与y的二元函数)得;(2)解上述以为未知数的方程,即得注 关于方程F(x,y,z)=0确定隐函数z=z(x,y)有以下的隐函数定理:设F(x,y,z)在点(x0,y0,z0)的某个邻......
2023-10-27
相关推荐