转子动力学 | 模态分析 附转子动力学钟一谔下载
模态分析是一种用于确定结构振动特性的技术:
• 固有频率
- 结构在什么频率下会自然地振动
• 模态形状
- 结构在每种频率下会以什么形状振动
• 模态参与系数
-每种模态参与给定方向的质量
模态分析是所有动力学分析类型中最基本的分析。
• 允许设计避免共振或以指定频率振动(例如扬声器箱)。
• 使工程师了解设计如何响应不同类型的动态负载。
• 帮助计算其他动态分析的解决方案控制(时间步骤等)。
建议:由于结构的振动特性决定了它如何响应任何类型的动态载荷,因此通常建议在尝试任何其他动态分析之前先进行模态分析。
转子动力学假设:
• 线性转子动力学
- 表示响应频率与激励频率相同
- 对两个组合激励的响应是两个单独应用的总和
- 这是转子从平衡位置小位移的合理假设(除非转子轴承系统出现灾难性故障)
• 稳态响应
- 所有瞬态都结束了
单自由度系统的自由振动:
- “自由”意味着没有强迫项
- EOM用于平行连接弹簧和阻尼器的刚性质量块
- 一个自由度系统的EOM的一般形式:
- 固有频率:
- 阻尼比ζ可通过以下方式获得:
假设为阻尼EOM形式的试算解:
特征方程可写为:
这个特征方程的根:
基于ζ值
复平面中不同情况的根位置
• 如果系统中的阻尼量变大,响应将不再振荡。
• 临界阻尼定义为振荡和非振荡行为之间的阈值。
• 阻尼比是系统阻尼与临界阻尼之比,由下式得出:
- 它是一种旋转机械,相当于单弹簧质量阻尼系统,在无质量弹性轴上具有集中质量
- 它是一个单自由度系统,通常用于引入转子动力学特性
- 本章将使用下面显示的Jeffcott转子来演示转子动力学的一些概念
无质量弹性轴轴长米,轴直径= 0.06米
盘直径= 0.5m,盘厚= 0.06m,盘质量m = 91kg
• 我们将讨论两种简单的Jeffcott转子:
• Jeffcott Rotor 1
-转子安装在非常刚的轴承上,形成“轴模式”
-- 与轴承和基础支撑相比,转子轴非常柔
-- 也称为简单的“柔性转子”
• Jeffcott Rotor 2
-转子安装在非常柔的轴承上,形成“轴承模式”
-- 与轴承和基础支撑相比,转子轴更加刚
-- 也称为简单的“刚性转子”
• 找到非旋转Jeffcott转子的前三种模态形状和相应的频率
- 假设转子没有旋转
- 平面运动的结果
- 轴端边界条件
• 情况1- 非常刚的轴承(考虑固定- 固定)
• 案例2- 软轴承,刚度为1* 10 ^ 6 N / m
案例1:考虑非常刚的轴承支撑(或简单支撑的轴端):
惯性矩:
轴刚度:
自然频率:
案例2:考虑软轴承支撑:
总刚度
自然频率
注意:添加轴承时,系统总刚度会降低,软轴承的第一固有频率会降低
注意:轴承刚度与轴刚度比对模态形状有很大影响•由于没有施加旋转,因此在这种情况下没有陀螺运动
• 对于绕轴旋转的结构,如果围绕垂直于Z轴的轴旋转(进动运动)应用于结构:
–反作用力矩称为陀螺力矩,其Z轴垂直于旋转轴和进动轴。
–由此产生的陀螺矩阵[g]对垂直于旋转轴Z的平面上的自由度进行了配对。–由此产生的陀螺矩阵[g]将是斜对称的。
• 考虑使用恒定角速度Ω围绕Z轴旋转的圆盘
- 允许质心在X,Y和Z平移
- 转子也可以在X和Y方向上旋转少量Ѳ和ψ
EOM用于柔性支撑的旋转转子,在x和y方向具有不同的刚度:
其中:简易转子模型的自由度、对称质量矩阵、斜对称陀螺矩阵、对称刚度矩阵
旋转:当旋转结构以其共振频率振动时,自旋轴上的点发生轨道运动,称为旋转。
后旋(BW):当转子在旋转轨道上以与转速Ω相反的方向旋转时
前旋(FW):当转子在旋转轨道上以与转速Ω相同的方向旋转时
轴中心椭圆旋转轨迹(本例显示了后旋(BW) ,由于轴承支承是各向异性的,因此会导致椭圆振荡)
• 坎贝尔图仅在模态分析结果中有效
• 由于陀螺效应,旋转分量的固有频率随旋转速度而变化
• 坎贝尔图用于绘制旋转结构部件在不同旋转速度下的变化动态特性
• 临界速度是与结构的共振频率(或频率)对应的转速
• 当固有频率等于激励频率时,出现临界速度
• 激励可能来自与旋转速度同步的不平衡激励或来自任何异步激励
激励线
旋转频率=斜率*旋转速度
阻尼模态求解器
打开阻尼
对阻尼模态的后处理支持
输出表是工程项(阻尼频率,稳定性等),而不是数学项(实/虚特征值)
• 在Workbench求解结果中,仅列出虚部的正值项以避免混淆
• 在solution information 中,所有8个虚模态频率都以旋转速度列出,该旋转速度是QR阻尼求解器的输出
陀螺矩阵- 无能量耗散
• 在公式中列出了具有阻尼项的陀螺矩阵,因为它取决于旋转速度
• 但它不会导致任何能量耗散
• 因此,除非存在一些阻尼(如轴承阻尼),否则固有频率的实部始终为零
在阻尼模态分析动画期间,打开/关闭时间衰减以进行切换控制循环次数以显示衰变的影响
• 为轴承单元增加阻尼(800 Ns / m)
• 会导致能量耗散
• 提供非零实部
在垂直于自旋轴的平面中,节点的轨道是椭圆形它由三个特征定义:局部坐标系(x,y,z)中的半轴A,B和相位ψ,其中x是旋转轴角度φ是节点相对于主轴A的初始位置。
轨道图只可用于梁模型
打印轨道:PRORB
动画与轨道图
悬臂转子型号4000rpm ,前两种模态
SET,4,3
/show,png
/view,,-10,-2,8
/rgb,index,100,100,100,0 !white
background
/rgb,index,0,0,0,15
PLORB
注意:这是一个梁和点质量悬臂转子模型,点质量不可见
案例
案例研究了悬臂转子的行为。旨在说明如何创建坎贝尔图并理解前旋和后旋。
问题陈述:
• 该模型由一个代表转子和轴的简单装配文件组成。
• 目标是估计0到4000 rpm之间的前六个固有频率和模态形状,创建Campbell图,观察由于模态对的分离而看到的陀螺效应(前旋和后旋)
轴承建模(ROMAC)
ROMAC Tilting Pad Bearing程序生成的ASCII文件-ANSYS轴承单元COMBI214
• COMBI214单元用于模拟轴承
• COMBI214在二维应用中具有纵向和交叉耦合功能
• 它是一个拉伸压缩单元,每个节点最多有两个自由度:
- 任意两个节点方向(x,y或z)的平动
• 轴承宏IMPORTBEARING1.MAC用于从ASCII文件导入与转速相关的轴承特性
ANSYS轴承单元COMBI214
轴承支撑系数:
• ASCII文件由THPAD创建(ROMAC的轴承分析软件)
• 轴承宏将通过Command Snippet在分析中输入
• 轴承特性放入APDL表参数中
• 然后可以在COMBI214轴承元件实际中使用这些表常数定义。
• 临界速度图可用于显示转子临界转速相对于轴承刚度的演变
• 水平刻度表示支撑刚度,垂直刻度表示转子速度,rpm
• 显示轴承刚度对临界速度的影响
• 绘制曲线,计算轴承平均恒定刚度的不同值的前几个横向临界速度
需要将以下输入作为Workbench Mechanical中的宏的参数提供
参数:
ARG1 临界速度数
ARG2 每转激励次数
= 0默认为1.0(同步激励)
轴承刚度
ARG3 最低值
ARG4 最高值
ARG5 步数= 0 默认为10
旋转速度方向(归一化为单位)
ARG6 X分量
ARG7 Y分量
ARG8 Z分量
• 宏CRITSPEEDMAP的特征
- 每步最多修改COMBI214轴承单元的2个实常数
- 仅将实常数K11和K22设置为ARG3至ARG5指定的非零值(K11 = K22)
- 然后通过求解特定的特征值问题直接获得临界速度
- OMEGA ARG6 ARG8命令将用于根据到定义整个结构的旋转速度
- 模型无阻尼
- 确保除了轴承阻尼之外没有定义阻尼,这将自动抑制
- 临界速度没有排序算法
- 它们的打印/显示顺序与计算顺序相同
- 使用Log-Log标度
- 轴承刚度增量在对数刻度上是恒定的
旋转结构中的自激振动导致振动幅度随时间增加,如下所示
如果不加以控制,这种不稳定性可能导致设备损坏。
最常见的不稳定因素是:
i. 轴承特性(特别是当存在非对称交叉项时)
II . 内部旋转阻尼(材料阻尼)
III . 旋转和静止部件之间的接触