转子系统非线性动力学DEM建模研究简介.pdfVIP

转子系统非线性动力学DEM 建模研究的简介1) 闰 民 陈予恕 曹树谦 (天津大学机械学院力学系，天漳∞∞72) 麓要将DEM方法引入转子系统非线性动力学的分析计算，对其建模过程和有关本构关系，作 了比较详尽的舟绍。该模型可用于碰摩、弯曲振动和扭转振动等转子系统的大部分动力学和运 动学问暮分析计算．不仅适用于单跨转子系统．也适用于多跨转子系统，尤其是在碰摩计算时． 可考虑其摩擦特性齄时问历程而变的特点，从而使分析计算更加接近实际情况。 关键词转子系磊讳线性动力学．DEM方法．转子碰摩，本构关系 人 引言 转子连同它的轴承、油膜、支座以及在运动中可能与之发生接触的部件被称为转子系统， 典型的大型高速转子系统有汽轮发电机组、燃气轮机、离心压缩机、以及各类大型旋转机械 等。转子系统中各种异常振动的存在，常常引发灾难性的事故f1】。从20世纪初，H．H．Je位on 以单圆盘挠性转子为模型伫】对转子系统的动力学行为研究以来，建立转子系统的动力学模型 一直是人们最关注的问题之一。在转子系统中，存在大量的非线性因素，这些非线性因素对 转子系统的局部和全局动态特性具有重大影响。一现在已是不争的事实。如从转子系统的工 作过程看．若转子的挠曲变形或涡动幅度大于某一数值时，它的某些部件就可能与其相邻的 零部件(如机壳或轴承等)发生接触和碰撞，习惯上把这一现象称为“碰摩”，不断的碰摩 导致巨大的转动能量向弹性变形积聚，最后可能引起转子的突然断裂。这种转子与其它零部 件发生碰摩的决定因素既是多方面的池是极其复杂的。碰摩是否发生以及碰摩发生的初始状 态既第一次碰摩的位置、碰摩力与摩擦力的大小以及它们之间的数值关系等一系列初始条 件，决定着转子系统以后的工作性态。此外在碰摩过程中，转子与其相邻的零部件有时接触 有时分开的离散特点，使碰摩的发生具有了随机特点。上述原因为转子系统的碰摩理论研究 和数值模拟带来一定的困难．尤其是在建立其动力学模型时表现的更为明显。尽管目前在这 方面的研究甚多【·】，但是无论是非线性研究方法还是研究内容本身，都还处于发展和完善之 Me廿坶d)简称DEMMI，用于转子 中。在这种背景下，本文将离散单元法(Dis∞咖EI啪em 系统．试图运用DEM建立转子系统较为完善的动力学模型，以便对转子系统的复杂动力学 行为进行研究。 DEM方法最早是由PA．cund“l于1971年提出的№卅，是研究不连续体力学行为的一种 新的数值方法，它可以细致地模拟备离散单元的相互作用。这种数值方法明显的优点是适用 于模拟离散组合体的接触或碰撞过程。DljM不同于基于最小势能变分原理的有限单元法，而 是建立在最基本的牛顿第二定律之上，它的基本思想，可以追溯到古老的受力分析方法上， 任何一个物体作为分离体来分析，总会受到相邻物体的力和力矩的作用，正是在其台力和合 力矩的作用下，产生变形和运动，这种思维方式在计算机上很容易实现，以每个单元为受力 体，建立其显式方程组之后．根据牛顿第二定律井结台不同的本构关系，以差分 法进行求解和迭代计算，结合动画技术，可形象直观地反映系统中每个单元运动变化的力场、 位移场、速度场等力学参量的变化。运用DEM建立的转子系统动力学模型及其单元间的本 构关系如图l所示。这里认为轴颈与轴承、转子与机壳之间的作用力具有相同的本构关系， 即都被简化为一对法向和切向弹簧阻尼器 辟C和野e以及一个切向滑动摩擦器 p。值得说明的是，弹簧阻尼器丘．C和滑 动摩擦器Ⅳ是联合作用的，当两单元之间 的切向力大于摩擦力时．两单元之间即发 生相对滑动，这时滑动摩擦器起作用，弹 簧阻尼器不起作用；否则，弹簧阻尼器起 作用，滑动摩擦器不起作用；此外轴颈和 图1力学模型及本构关系 轴承单元的作用，当最小油膜间隙大于等 嗣萨 于零或大于某一数值时，即认为油膜未被破坏时，疋．G和砰e由油膜力确定，且F=0； 反之以上参数则按碰摩情况确定，具体计算方法在下面的介绍中将会提及。轴颈与转子单元 的