SAMCEF有限元转子动力学分析实例.docVIP

第章 15.1转子动力学理论背景 15.1.1 概述 带有旋转部件的整体结构，比如说飞机引擎，它的动力学行为分析需要旋转部件、定子部件和不同的连接设备的模型。模型的处理过程采用的是有限元方法。15.1.2 转子 有三种有限元模型可用来描述转子系统。 1D模型：转子用梁、弹簧和集中质量单元来模拟，这个模型计算速度快，它适用于有大量参数需要调整时的初期设计分析。但是，比如前文提到的飞机引擎，这个模型的细节可能需要许多专业的经验还有可能耗费很多时间。 图2.1 梁-弹簧-集中质量模型 傅立叶模型：转子采用2D傅立叶多谐波单元模拟。对于转子动力学方面的应用，为了描述轴向变形、扭矩变形和弯曲变形与陀螺力矩耦合，必须考虑0和1节径的谐波。这种模型允许旋转设备的更精细模型的开发，它非常适合于带有多数目叶片的旋转机械的建模，也适合圆锥杆的更好的建模 图2.2 2D傅立叶转子模型 为了促进轴对称转子模型和与它耦合的固定部件或其它旋转部件的使用，还要考虑陀螺效应和阻尼的影响，更新了多谐波实体单元和壳单元库。 为了考虑转子的离心刚度，同时采用不同谐波是一个简单的方法，在瞬态分析中，它可以在同样的运行中做到。另一方面，单元库包含一种专门的连接单元，这种单元可以考虑3D模型和傅立叶模型之间或轴上的一个3D节点和傅立叶模型之间的耦合。 位移场采用傅立叶级数形式展开，有限元建模在参考子无午面上进行。 3D模型：用体单元（六面体、棱柱或者四面体）或者壳单元描述转子系统。 图2.3 3D转子模型 当需要加速计算时，就会想到对转子模型进行简化，这可以通过超单元来实现。 位移场可以通过自由度矢量q描述，在惯性坐标系中，惯性力通过动能的离散表达式得到的。 这里： M是质量矩阵 G是陀螺矩阵 Ω是相应转子的角速度 由于转子和支撑设备的弹性，陀螺效应随着转子向量的方向而变化。按照前面的规定，旋转速度是一个将要解决的已知问题。 弹性力是通过惯性坐标系里的应变能的离散方程得到的。 这里K是刚度矩阵。由于初始装配、离心荷载或者任一其他的荷载存在预应力，弹力就有如下修改： 这里K?是几何刚度矩阵，这里的刚度不一定非要各向同性，但是为了保持常数矩阵在圆周的方向，各向同性是需要的。 如果假设有粘性阻尼，阻尼力通过分散函数来表达。 这里qr 是旋转坐标系里的广义位移。在惯性坐标系里，得到的粘性阻尼力有如下形式： 这里： B是粘性阻尼矩阵； BCirc 是补充阻尼或者循环力矩阵。 通常，旋转组件中的阻尼和装配时摩擦力有关。粘性模型就不适合了，滞后模型是比较适合的。滞后阻尼的数学表述仅在频域有效。在时域里描述阻尼行为的一种方法是用可变的粘性阻尼。在这种情况下，粘性阻尼力就由特征频率给分开了。 注意这种方法是近似的，更好的方法就是在装配阶段包含摩擦力。 15.1.3 定子 从单自由度的等效弹簧到完全的3D模型，定子模型是比较典型的，并且非转动部分可以被有限元库来描述。对于转子动力学，在定子部分也要引入阻尼。 qS的自由度向量描述了定子的位移域，惯性力由惯性坐标系下的动能离散方程来表达。 这里，MS是定子质量矩阵，弹力来源于惯性坐标系中的应变能的离散方程 这里KS是刚度矩阵。预加应力的二阶影响可以考虑进来。 如果阻尼是粘性阻尼，阻尼力源自于分散方程。 粘性阻尼力有如下形式 这里BS是粘性阻尼矩阵。 通常，旋转组件中的阻尼和装配时摩擦力有关。粘性模型就不适合了，滞后模型是比较适合的。因为这个模型没有遵从因果关系，所以滞后阻尼的数学表述仅在频域有效。在时域里描述阻尼行为的一种方法是用可变的粘性阻尼。在这种情况下，粘性阻尼力就由特征频率给分开了。 图3.1 定子模型 15.1.4 简化模型 在整个引擎中使用的简化模型主要是指超单元，超单元通过简化的刚度、质量、阻尼矩阵和荷载向量描述了子结构。如果子结构是转子的一部分，那么就会产生回转和循环力矩阵。 产生这些简化矩阵有很多方法，当前最常用的是GUYAN和CRAIG和BAMPTON简化方法。第一个是第二个的特殊情况，CRAIG和BAMPTON方法是组合模态方法的一部分，作为模态的两种，组合模态方法被用来产生简化模型。 注意GUYAN方法对线性静力问题是准确的，在动力学分析中，和CRAIG和BAMPTON方法一样，它是近似的并且可以通过使用附加模态来改良。 当部件模态方法在动力学分析中使用时，每一个子结构的行为通过两种类型的部件模态的联合来表述：约束模态（静力变形）和子结构的正则振动模态（独动态变形）。这些部件模态由自由子结构的无阻尼动力学方程决定，自由度的分类是通过几个子结构共享的边界自由度，或者是仅考虑子结构的内部自由度来划分的。强迫模态由赋值、每一个边界自由度的单元位移和其它固定的边界自由度来决定。正则模态对应的自由模态是通过固定子结构边界来