带非线性阻尼的转子动力学有限元的研究.pdfVIP

带非线性阻尼的转子动力学有限元研究 洪杰赵险峰李其汉 (北京航空航天大学动力系北京100083) 育多十非缋性挤压油膜的转子系统的计算分析寰明．在通用有限元程序基础上该计算方法具有通用性强．计算精度高 和数值稳定性好的优点． i 关t词航空发动机转子动力学挤压油蔗阻尼器SASl_P．Am 1引言 由于转子动力学的特殊性．以往的转子动力特性分析程序大多是用某种算法语盲依据转子动力学 理论独立开发而成。计算程序的兼容性差使用效率低．本文在现有通用有限元程序NASTRAN平台上， 通过加入陀螺矩阵和挤压油膜阻尼器油膜力阵．重新组织求解序列以此宴现对转子动力特性的计算分 析。这种分析方法的优点表现为(1)充分利用各种成熟稳定的求解模块．从而提高了程序的稳定性和 计算结果的可靠性：(2)利用优秀的前后处理功能提高计算数据的可视性：(3)利用通用的数据格式， 提高了鼓据的可移植性：(4)改善计算模型的精度。但是，在目前的通用有限元软件上进行转子动力 特性分析存在着两个难题：一是如何在转子系统的动力分析中考虑陀螺力的影响，因为现有有限元软 件都无法计算与结构转动和进动相关的反对称陀螺矩阵；=是由于挤压油膜阻尼器油膜力的高度非线 形特性，在剐度和阻尼矩阵中出现非线性项，现有有限元程序无法进行稳定高效求解计算。本文正是 针对这两点利用NASTRAN提供的二次开发程序DMAP实现了对转子系统的动力学分析。 2转子系统的有限元模型 转子系统的有限元建模分为两种方案：～是进行发动机的整机建模，再应用GUYAN【静凝聚、 自由度减缩技术把质量、刚度等模型参数凝聚到沿轴心线分布的各点上进行计算求解：．二是利用较传 统的有限元粱轴模型考虑陀螺惯性的影响，对转子进行建模、分析。这种模型较为简单，讨。算量少， 对模型描述的精确程度虽不及全机建模但计算结果还是令人满意的。 2．1转子系统的粱轴离散模型特征 对于转子系统有限元离散模型划分时，以铁木辛科粱理论为基础．轴作为等直径圆截面粱处理， 考虑粱的剪切变形、转动惯性、陀螺力矩和连续质量：机匣简化为转速为零的转子．机塑上的附件当 作集中质量或分布质量加在机匣的相应位置；叶盘及封严盘等作为集中质量点处理，考虑其质量和转 动惯量；轴承及一些联接结掏(联轴器、腹扳等)简化为线性的刚度和阻尼单元：航空发动机的转子 结构通常十分复杂，一些较为特殊的结构如联轴器、惟形腹板等，很难用一种有效的有限元模型来表 达。在此．可以把这些结构简化为广义的弹簧单元和阻尼单元．单元刚度和阻尼的大小可虬由实验测 定或依据经验数据而定。因此。在对转子系统的有限元分析中．考虑到各部件的具体结构单元可划分 成粱单元、集中质量单元，轴承单元和广义单元。 矩阵，只是额外考虑了陀螺力矩的影响，从而保证了与NAsl限^N单元库的～致性。 2．2挤压油奠阻尼嚣的敷学模型 在对挤压油膜阻尼器建立数学模型中可将油膜力视情况作为非线性外力．或者非线性刚度和阻尼 加入系统运动方程。 挤压油膜阻尼器油膜力示意圈见图1所示．按照挤压油膜阻尼器的短轴承半油膜解，挤压油膜阻 尼器油膜力在固定坐标系下瞬态和稳态时计算公式见公式(1)、(2)所示为： ——302— ，．．。，。。，fjiiij。；!；：iiij+；i；；2：；}：i誊+]‘1a’ ．芦了I下一万i薯磊蔫． 、＼ ¨b) 一+—!苎-=二4—1 2 ：fc‘一、一：’】： ； 形心f ．————生■L三1——一l “：一17一：!I：，1．+：：I’：f r=一un。[；i；：却+；：?!÷；j!筹]c2a) 在公式中．“为滑油动粘度．R为轴颈半径．L为有效油膜长度。C为半径间橡．e为油膜轴颈值 心率．从油膜力计算公式中可以看出．等效刚度系教和等效阻尼系敷与挤压油膜阻尼嚣内外环的相对 位移(或速度)相关，表现出非线形的性质．在本文进行的研究中，将采用固定坐标系中的瞬寿和穗 态油膜力计算公式，并不做线