静态下的无阻尼系统振动平衡方程.docVIP

第4章9-500型无极捻筒体振动分析 本次测试是针对大型的刚性旋转部件，振动成因较为复杂，因此在测试之初对设备进行了模型的建立，通过仿真得出的参考值以提升诊断的信心。同时通过研究设备动态特性，可以更好的解读设备运转振动模式（振型），为以后的设计工作提供依据。 4.1求解固有频率 静态下的无阻尼系统振动平衡方程： （4-1） M：等效质量 K：等效刚度 X：位移向量 则方程的解： （4-2） ， ：第i阶振型的固有频率 ：第i阶振型的频率相位角 将（4-2）带入（4-1）得 （4-3） 令 因振动系振幅不全为0，固有 （4-4） 由（4-3）得 （4-5） 令则 （4-6） 是矩阵D的特征向量 初始迭代向量 （4-7） 第一次迭代 （4-8） 将式（4-6）带入（4-8）得 （4-9） 第二次迭代 第m次迭代 （4-10） 由（4-10）可知 至此可得 对于第二阶，第三阶以至于更高阶次的固有频率，可以令，，……这里不再赘述。迭代计算的特征是重复且有规则的计算，因此可以采用计算机计算结果，本文所用的是FFEPLUS解算器计算结果如下表4.1所示： 频率数 弧度/秒 赫兹 秒 1 156.79 24.954 0.040073 2 158.59 25.24 0.039619 3 386.08 61.446 0.016274 4 512.67 81.594 0.012256 5 878.85 139.87 0.0071493 6 892.93 142.11 0.0070366 7 1135.2 180.68 0.0055348 8 1226.1 195.14 0.0051245 9 1306.8 207.98 0.0048081 10 1421 226.16 0.0044216 11 1450.2 230.81 0.0043325 12 1526.2 242.91 0.0041168 13 1528.3 243.24 0.0041111 14 1629.2 259.3 0.0038566 15 1643.2 261.52 0.0038238 16 1688.4 268.71 0.0037214 17 1711 272.31 0.0036723 18 1774.4 282.41 0.003541 19 1838 292.52 0.0034185 20 1853.2 294.94 0.0033905 表4.1 选取了部分振型的图像如下图所示，详见附件1 名称 类型 最小 最大 位移1 URES:合位移 图解 对于模式形状: 1(数值 = 24.9545 Hz) 0 mm 节: 874 44.798 mm 节: 8981 复件 筒体-算例 1-位移-位移1 名称 类型 最小 最大 位移3 URES:合位移 图解 对于模式形状: 3(数值 = 61.4462 Hz) 0 mm 节: 874 49.3169 mm 节: 72622 复件 筒体-算例 1-位移-位移3 名称 类型 最小 最大 位移5 URES:合位移 图解 对于模式形状: 5(数值 = 139.873 Hz) 0 mm 节: 874 47.4884 mm 节: 2080 复件 筒体-算例 1-位移-位移5 4.2求解周期激振的响应 电机输出的磁力矩是周期性变化的，假设其周期函数为，则其输出力为 （4-11） 将周期函数展开成傅里叶级数： （4-12） 为傅氏系数。 将振动方程进行正则化后可得出： （4-13） 则其第i阶有阻尼稳态响应为： （4-14） 为放大因子： 为相位角: 由式（4-14）可知，任意阶正则坐标的响应是多个具有不同频率的激振力引起响应的叠加，因而周期性激振函数产生共振的可能性要比简谐激振高的多。所以很难预料各个振型中哪一个振型将受到激振力的强烈影响。但，当激振力函数展开成傅里叶级数之后，每个激振频率可以和每个固有频率比较，从而预测出强烈振动的形式。 4.3测试数据 本次测试是针对管式捻股机筒体部分的测试。同样采用了AMSsuite的数据库路径分析，简化了测试流程，同时采集了2种不同情况下的数据，对诊断