转子动力学基础.ppt

1、2020/10/9,1,第一章 转子动力学基础,本章主要内容： 涡动分析、临界转速 重力影响 弹性支承影响 非轴对称转子影响、稳定性问题 初始弯曲影响 等加速过临界的特点,2020/10/9,2,第一节 转子的涡动,旋转的转子是具有质量和弹性的振动系统，这与其他振动系统相同。 区别：转子是旋转的 涡动：既有自转，又有公转，是一种复合运动。 不平衡力引起的同步正进动分析,2020/10/9,3,第二节 Jeffcott转子涡动分析,Jeffcott转子：垂直安装等截面对称转子、不计重力影响。 一、Jeffcott转子运动微分方程 Jeffcott转子示意图 薄盘：h/D0.1；偏心矩：e 定坐标

2、系：oxyz；基点： 设自转为常数，确定 的运动： x(t)、y(t) 或 r(t)、(t) 假设：扭转刚度无限大(不计扭振) 忽略轴向位移、刚性支承 轴的弯曲刚度为EJ E：弹性模量 J：截面惯性矩,2020/10/9,4,轴的弹性恢复力在坐标轴上投影为： k轴的刚度系数 对称简支梁中点刚度为： 粘性外阻尼力在坐标轴上投影为： c粘性阻尼系数 由牛顿定律可得： 由几何关系可知：,2020/10/9,5,两边对时间求两次导数得： 代入牛顿方程得 点的运动微分方程 根据动量矩定理，可得圆盘绕重心c转动的微分方程： 对于稳态涡动，,2020/10/9,6,代入牛顿方程得 点的运动微分方程 化为标准

3、形式为： 式中： 弹性轴无阻尼横向振动固有频率 相对阻尼系数 运动微分方程与线性阻尼系统强迫振动相同，可设解为,2020/10/9,7,代入运动微分方程解得: 点作圆周运动，参照极坐标几何关系: 故运动半径为轴的动挠度r，为动挠度r与偏心矩e间的相位差，且有:,2020/10/9,8,2020/10/9,9,低转速区,共振区,高转速区,圆盘重边飞出,圆盘轻边飞出； 自动定心或质心转向,2020/10/9,10,临界转速定义(ISO)：系统(位移)共振时主响应的特征转速。 主响应：轴颈运动或转子挠曲 对于Jeffcott转子，临界转速对应 常以cr或c表示，若以转/分或转/秒为单位，则有 或 将

4、转子挠度表达式代入临界转速条件得 解得 可见，阻尼总使临界转速大于横向振动固有频率，与机械振动中的阻尼使固有频率降低作用相反。 当转子系统阻尼很小时，可近似认为： 此时有,2020/10/9,11,p时，/2，与阻尼系数大小无关，利用这一特点可测取转子系统的p，在小阻尼情况下可近似为临界转速。 当0时，p时，0， 三点在一条直线上 p时， 三点在一条直线上 p时，/2，r,不同转速下圆盘偏心位置见图114,2020/10/9,12,=，同步正涡动，或正协调进动； =-，同步反涡动，或反协调进动； ，同方向，正涡动，或非协调正进动； ，反方向，反涡动，或非协调反进动。 当转子圆盘不在中间时，即使

5、是无阻尼系统，其临界转速p，主要是陀螺力矩影响。,同步正进动轴的受力,2020/10/9,13,例：已知：轴长l=57cm，直径d=1.5cm，轴材料弹性模量 ，圆盘厚度h=2cm，直径D=16cm,材料密度 ，不计阻尼。 求：1）临界转速cr 2）e=0.1cm，0.6cr；0.8cr时的动挠度r及支反力幅值F。 解：弹性轴质量： 圆盘质量： 弹性轴中点刚度： 不计轴质量时临界转速：,2020/10/9,14,计入弹性轴等效质量，按照振动理论，梁在中点的等效质量为原质量的17/35，则临界转速为：,0.6cr时挠度为：,支反力幅为：F=kr=74.562N,轴承力与重力之比为：,2020/1

6、0/9,15,0.8cr时挠度为： 支反力幅为：F=kr=235.68N 轴承力与重力之比为：,2020/10/9,16,第二节 刚体绕定点的转动,力学模型：连续质量模型弹性体 集中质量模型盘轴系统 本章以盘轴系统为分析模型 刚体在空间有六个自由度：沿三个垂直轴方向的平移和绕这三个轴的转动。 理论力学：刚体运动可分解成随基点的平动和绕基点的转动。 平动运动规律与基点选择有关； 转动运动规律与基点选择无关。 1.2.1 描述定点刚体位置的欧拉角 刚体球铰定点约束：约束三个平动自由度； 只有三个转动自由度。,2020/10/9,17,定坐标系oxyz与动坐标系的关系 见表11和图16 关系式为：,

7、2020/10/9,18,各方向余弦存在关系： 因此，九个方向余弦中只有三个是独立的（自由度数）。 方向余弦求解复杂，采用夹角欧拉角表示，多种定义。 1、第一种定义（图17）： 1)动坐标与静坐标重合，先绕oz轴转动角进动角； 到达oNN1z，oN称为节线，右手法则 2)绕oN轴转角方位或挠曲角； 到达 3)绕 转角自转角； 到达 引入坐标轴矢量 、,2020/10/9,19,再引入oN、oN1及 的单位矢量 ，则有： 由于： 得到：,2020/10/9,20,2、第二种定义（图18） 1)动坐标与静坐标重合，先绕oy轴转动角,到达ox1yz1； 右手法则 2)绕ox1轴转角,到达 3)绕 转

8、角自转角， 到达 、结合体现进动与方位角。 令ox1、oy1、oz1单位矢量为 则有,2020/10/9,21,由此可导出欧拉角的三角函数表示的方向余弦：,2020/10/9,22,欧拉角表示的刚体绕定点转动的运动为 或 1.2.2 刚体绕定点运动的角速度及速度分布 刚体的角速度为 或 所在的位置称为刚体绕定点转动的瞬时转动轴，瞬时转动轴时刻不同，但总通过定点。 第一种定义法得到矢量 向定坐标系投影得,2020/10/9,23,利用方向余弦关系得 向动坐标系投影得 类似，由第二种定义可得 向定坐标系和动坐标系的投影 刚体上任一点瞬时速度矢量为,2020/10/9,24,将速度向定坐标系和动坐标

9、系投影得 刚体上各点角加速度和加速度为 1.2.3 刚体作定点转动时的动量矩定理 动量矩定理：刚体对定点o的动量矩 对时间t的导数，等于外力系对该点的主矩 则有 对有集中质量的刚体，动量矩为 刚体在绝对运动中对质心的动量矩 ，等于刚体随质心平移动坐标系中运动的相对于质心的动量矩 。,2020/10/9,25,因为由速度合成定理： 则刚体相对质心的绝对运动动量矩为 由于刚体对质心的质量矩等于零，即 因此 若将固定点取在质心o上，则有 在相对随质心平移的动坐标系中，刚体对质心动量矩对时间的导数等于外力系对质心的主矩刚体相对质心的动量矩定理。 因此，对质心动量矩的计算只需考虑相对转动。 刚体作定点转

10、动时，有 刚体动量矩为,2020/10/9,26,向动坐标系投影得 式中： 为刚体对 轴的惯性矩 为刚体对 、 轴的惯性积 为刚体对 、 轴的惯性积 对一般具有圆截面的均质轴对称转子有 对均质薄圆盘有 式中：m圆盘质量 R圆盘半径 类似可得 于是,2020/10/9,27,如果 为刚体对o点的主惯性轴，则各惯性积为零，即 于是有 一般情况下的矢量关系如图19。 若刚体对动坐标系的惯性矩为常数 则有 式中： 欧拉动力学方程,2020/10/9,28,1.2.4 刚体运动的动能 能量定理、拉个朗日方程运动微分方程 设刚体质量为m，基点运动方程为x(t)、y(t)、z(t)，以基点为原点的动坐标系

11、是刚体的惯性主轴，惯性矩分别是 ，则刚体的动能为 通常转子沿oz轴方向的运动为二阶小量，可忽略不计，即有 z(t)=0 故转子的动能计算公式为,2020/10/9,29,第三节 单盘偏置转子的涡动、回转效应,转动惯量：反应刚体质量分布的力学参数。 中心极转动惯量： 绕通过执行的对称轴的转动惯量。 中心直径转动惯量： 绕通过质心的任一直径的转动惯量 均值等厚度圆盘，其转动惯量为： 圆盘的回转效应：转动的刚体有力图保持转轴方向不变的特性。转动物体的惯性的体现。,2020/10/9,30,三个圆盘的动量矩：,的方向沿轴线的切线方向。若转子以角速度,绕z轴转动，则动量矩的变化率：,2020/10/9,

12、31,动量矩定理：,圆盘在轴上的反力矩：,圆盘的回转力矩：,2020/10/9,32,回转效应：由于高速旋转圆盘的偏摆运动而使临界转速变化的现象（见图115）。 1.3.1 单盘偏置转子运动微分方程 假设：无阻尼、无偏心 不计轴质量 如图115，圆盘的轴线在空间 画出的轨迹是个锥面。 为分析方便，建立如下坐标系： (图116、图117） 1）定坐标系：oxyz 2）随 点平移坐标系： 3）固联于 动坐标系：,2020/10/9,33,其中： 是轴挠度曲线的切线 、 为两正交直径,2020/10/9,34,薄盘运动可以用xoz、yoz平面投影x(t)、y(t)表示。 采用第二种欧拉角定义有 故可

13、以用x(t)、y(t)、(t)、(t)、(t)确定圆盘空间位置，描述运动状态。 如图118， 点的挠度x和转角为 解出盘对轴的作用力Fx和力矩Mx为：,2020/10/9,35,式中： 式中和Mx的转向如图118所示。 在yoz平面也有类似公式； 为了使力矩矢量都沿坐标轴正方向， My与 Mx的转向规定相反 于是有 根据质心运动定理： 代入力关系式得 点横向运动微分方程为：,2020/10/9,36,化成标准形式： 式中： 由动量矩定理可建立圆盘绕 点转动的运动微分方程。 由于、都是小量，故有： 根据图117，三个轴的角速度为： 显然， 、 、 为圆盘的三个中心主 惯性轴。 令圆盘对 轴转动惯

14、量为Ip，对 、 轴转动惯量为Id 则有：,2020/10/9,37,对 轴的角速度就是自转角速度，即： 故对三个轴的动量矩为： 分别向 、 、 轴投影得： 根据对质心的动量矩定理，对圆盘有 由于 是质心，所以重力对 、 的矩为零。 假设作用圆盘上的所有外力对 的矩为零，则由上式得： 常数,2020/10/9,38,在此条件下，可得盘的偏摆运动微分方程： 与方程 联立求解 1.3.2 单盘转子涡动分析 设联立方程的解为： 代入联立方程，得到A、B、C、D的一次齐次方程组，根据非零解的条件，方程系数行列式的值应等于零，由此得到关于自然频率的高次方程，将解得的代回联立方程，,2020/10/9,3

15、9,可得相应的一组A、B、C、D之间的比值。 转子运动稳定时，动挠度曲线在动坐标系中是不变的；只是绕着oz轴进动，进动角(t)一般从ox轴量起。 令总的动挠度为r，挠曲角为，由图119几何关系得 显然，r、常数，且 常数 将以上关系式代入联立方程得,2020/10/9,40,由此可见：弹性轴发生弯曲不仅有离心力( )，还有回 转力矩 的影响；回转力矩改变了轴的弯曲刚 度。 上式是关于r、的齐次方程，由非零解条件得： 展开后得： 即 上式可求得得四个根，且随而变，令： 为同步正进动，则方程为,2020/10/9,41,存在一个正根（负根舍去）： 讨论：1）令 频率方程为： 解得： 式中： 为弹性

16、轴 点的横向刚度 此时得到的频率数值上等于转子不旋转时的横向固有频率，即不计回转效应时转子临界转速。 频率为： 2）令 为同步反涡动，频率方程为,2020/10/9,42,有两个正根(仍然假设为薄盘，即Ip=2Id)为 3）如果出现IpId的情况(如地面串联式离心压气机)，可能使 得 在正同步涡动情况下为负值，此时陀螺力 矩将降低临界转速。 例：已知：轴长l=57cm，直径d=1.5cm，轴材料弹性模量 ，圆盘厚度h=2cm，直径D=16cm,材料密度 ，a=l/4，b=l(3/4)，不计阻尼。 求：临界转速cr 解：Id=50.192 a=l/4=14.25cm b=l-a=42.75cm

17、3lEJ= k11=5498.589N/cm k12=k21=-67161.07N k22= 1)考虑轮盘回转效应的临界转速为,2020/10/9,43,1420208.1 则 cr=30/=2663.02 转/分 不考虑轮盘回转效应的临界转速为： 回转效应提高临界转速百分比为 对悬臂转子有类似结论，频率特性曲线如图121所示,2020/10/9,44,第四节 重力对临界转速影响、副临界,对水平安装的Jeffcott转子，重力影响： 1）重力产生静挠曲；航空发动机忽略、汽轮机扬度 2）质量偏心：交变力矩mgesint，如图122 薄盘极转动惯量为：Ip= 式中：为回转半径 轮盘角加速度为： 产

18、生一个切向惯性力： 垂直分量为： 式中的常值部分产生静挠度；交变 部分在垂直方向投影作简谐变化； 当cr/2时，此时，动挠度达到 极大值副临界。,2020/10/9,45,1.4.1 水平Jeffcott转子运动微分方程,水平安装的Jeffcott转子如图123 无阻尼、盘在中间、无回转效应 三点总在一条直线上， 如图124，弹性恢复力投影为 根据质心运动定理得 绕质心轴线的动量矩为 圆盘质量m对oz轴的动量矩为：,2020/10/9,46,根据动量矩方程有： 即 为了化简方程，作坐标变换，令 变换后联立运动微分方程为: 1.5.2 水平Jeffcott转子涡动分析 假设：Mz=0， ，则 令

19、的初值为零，则（cr/2）t，代入运动微分方程第三式： 联立方程的解为：,2020/10/9,47,式中：s=mg/k 为圆盘重力作用下弹性轴跨中静挠度。 根据图124的几何关系，圆盘质心运动方程为： 进动由个分量组成：基频分量与偏心矩e有关； 倍频分量与静挠度s有关。 例1.6：求第四节例子转子，重力作用下，副临界时两个运动分量幅值。 解：基频分量为：e=3=0.1/3=0.0333cm 倍频分量为：s=mg/k=0.0232cm 可见各幅值较小，只有s较大情况下，才能观察到明显副临界现象，轴横截面非对称是一个更为主要原因，下节讨论。,2020/10/9,48,第五节 弹性支承单盘转子涡动分

20、析,1.5.1 弹性支承单盘转子的运动微分方程 典型鼠笼式弹性支承结构见图125 弹性支承单盘转子计算模型见图126,2020/10/9,49,如图126，假设支承的质量、刚度、阻尼参数分别为： mb、kx、ky、cx、cy，支承转子的作用力为： 由于阻尼存在，应用非保守系统的拉格朗日方程： 式中：L=T-V拉格朗日函数 T系统的动能函数 V系统的势能函数 系统阻尼的耗散函数 Qj作用在系统上的广义力 qj系统广义坐标 n系统自由度数,2020/10/9,50,假设支承对称，自由度减为二个：xb(t)、yb(t)；圆盘位于中间，自由度减为三个：xc(t)、yc(t)、(t)，系统共五自由度 设

21、圆盘角速度为常数： 则系统缩减为四个自由度数。 圆盘动能： 支承动能： 圆盘势能： 支承势能： 圆盘外阻尼耗散函数： 支承粘性阻尼耗散函数： 因假设常数，合力矩为零，若忽略重力影响，则系统广义力为零。,2020/10/9,51,将以上各式代入拉格朗日方程得 1.5.2 弹性支承单盘转子涡动分析 为简化分析，不计阻尼，可设 将其代入运动微分方程得,2020/10/9,52,由第三式解得 代入第一式得 于是 分母为零，则X达到无限大，即x方向的临界转速cr，令 即 或 由求根公式可得 讨论：1）kxk，支承几乎不动，支承质量可视为零，mb=0, 频率方程简化为：,2020/10/9,53,2）kk

22、x，频率方程简化为： 支承动刚度定义：作用在支承上的简谐激振力与力方向上的振动位移之比。与支承刚度、阻尼、参振质量有关，一般为复数，是随而变化的曲线（特性）。也称为位移阻抗，用i除后，得到速度阻抗，也就是常用的机械阻抗(参见六章) 3）同样，若x、y方向互不耦合，则有 频率方程化简为： y方向临界转速为:,2020/10/9,54,一般kxky，所以通常XY，因而盘心 轨迹是一个椭圆。 4）设kxky，且kkx， kky，则有： 表明转子首先发生支承运动，即转子刚体涡动 运动轨迹为一椭圆柱面（图127）； 考虑转子转动惯量，还会出现一椭圆锥面（图128） 为研究涡动随自转速度变化的性质，盘心运

23、动采用复数表示,2020/10/9,55,令： 此时： Xc=X， Xs=0； Yc=0， Ys=-Y 由复数关系式得： 正、反进动两个分量相位角为： 由前面公式可得椭圆长半轴a、短半轴b及长半轴与ox轴夹角 如果XpXr，则为正涡动；如果Xpky，则：cxcy。 a）当0，YsXc 则：XpXr，为正涡动，a=|Ys|，b=Xc， 。,2020/10/9,56,b）当cy时，Xc0，Ys 则：a=，b=Xc， 。 c）当 时，Xc0，Ys0 则：Xp0 则：Xp0，Ys0，且XcYs 则：Xp0 则：a=， b=Xc， g）当cx时，Xc0，且|Xc|Ys 则：XpXr，为正涡动， a=|X

24、c|，b=Ys， h）当时，Xc=-e，Yc=e 则：XpXr，为正涡动， ab=e，,2020/10/9,57,盘心运动轨迹和方向如图129。 通常轴的两个方向刚度差异不大，两个临界转速靠得很近，一般不允许在临界转速附近停留，故一般只能看到正进动。,2020/10/9,58,第六节 非轴对称单盘转子的涡动分析,实际转子非对称典型结构见图130 假设支承刚度远大于轴的弯曲刚度，由 于轴刚度不对称，重力激励频率为转速 二倍。当转速到达临界转速一半时，会 使动挠度达到极值副临界。是引起 副临界的重要原因（图131）。,2020/10/9,59,1.6.1 非对称轴的单盘转子运动微分方程,如图131

25、，建立与定坐标系平行的坐标系 ，再建一动坐标系 ，令 、 为轴截面两个主惯性轴，转子以绕 轴转动，设截面两个主惯性矩为J、J，相应弯曲刚度为： 圆盘偏心为 ，相对动坐标系的相位角为e 则有： 动坐标系中牵连惯性力为 ，在、 轴上的分量为： 哥氏惯性力的分量为： 轴的弹性反力分量为：,2020/10/9,60,圆盘粘性外阻尼力 分量为： 式中：c圆盘粘性外阻尼系数。 弹性轴内阻尼力 分量为： 式中：ci弹性轴内阻尼系数(分析见后面章节)。 转子受到的重力 分量为： 根据质心运动定理，参考图132，得 令：,2020/10/9,61,质心运动方程为： 为线性非齐次方程 利用线性叠加原理 可分别讨论

26、重力和 偏心的作用。 1.6.2 重力作用下非对称单盘转子的涡动 重力单独作用时运动微分方程为： 运动方程右端采用复数表示，并规定取其实部，有,2020/10/9,62,设运动方程解为： 、 代入运动微分方程得： 求解上式，得 若忽略阻尼，即令c=ci=0，则,2020/10/9,63,上式中都以实部为运动方程 无阻尼情况下，取极值条件，即分母为零，则有 令： 代入上式得： 通常k/k1，k/m ，则有 表明会出现副临界现象。 根据坐标变换关系：,2020/10/9,64,可得固定坐标系下进动方程： 由此可见：重力使双刚度转子 以两倍自转角速度作正进动； 轨迹是以静平衡位置为中心的 圆，半径与

27、k成正比。 副临界峰值比临界峰值小； 利用副临界测临界转速(安全)。 1.6.3 偏心作用下非轴对称单盘转子的涡动 仅偏心作用下运动方程为：,2020/10/9,65,可见偏心引起的离心惯性力在动坐标系中如同一个静力，只能产生相对静挠度e、e，且有： 代入运动方程得： 解得 动坐标系中的动挠度re和相位角e为：,2020/10/9,66,可见：re随动坐标系一起以自转角速度旋转。 在固定坐标系中，圆盘作同步正涡动，盘心轨迹是以re为半径的圆。 讨论：1）响应与内阻尼ci无关； 2）对无阻尼（c=0），根据动挠度极值定义临界转速，有： a)支承不对称，两临界转速之间，转子作稳态 反进动涡动； b

28、)轴不对称，两临界转速之间，转子可能出现 不稳定。,2020/10/9,67,第七节 弹性轴有初始弯曲时的转子涡动,设：一Jeffcott转子，不计重力，无质量偏心，轴有初始弯曲 圆盘处初始弯曲为rs，转子以角速度旋转，进一步弯 曲rd，动挠度为r。见图133 s：初始盘心位置 c: 瞬时盘心位置 定坐标系：oxy 动坐标系：o 建立动坐标系运动微分方程 1）盘心绝对加速度为： 2）弹性轴恢复力：,2020/10/9,68,3）轮盘粘性外阻尼力： 4）根据质心运动定律得： 化为标准形式微分方程： 代入坐标变换，得定坐标系方程： 由于旋转，xs、ys随时间变化，可表示为： 代入定坐标系运动方程得： 与质量偏心运动方程形式相似，但特点不同。,2020/10/9,69,令： 式中： 式中： 代入运动方程得： 由强迫振动可知方程稳态解为： 与偏心盘的比较： 1）动挠度： 偏心盘分子上因子为 轴弯曲分子上因子为 2）动挠度极值频率： 轴弯曲为： 偏心盘为： 无阻尼时两者相等。 3）同时存在轴初弯曲rs和盘偏心e，并假设相位差位，由线性叠加原理得：,2020/10/9,70,第八节 转子在越过临界转速时的行为,临界转速附近挠度最大，主要讨论越过临界转速时的行为。 假设：转子等加速（或等减速）越过临界转速 1.8.1 变转速时转子运动微分方程 如图134，以水平Jeff