转子动平衡的几个优化算法.pdf

2003年第 1期 东方汽轮机 7 转子动平衡的几个优化算法 胡廷瑞 摘 要 ：介绍了转子动平衡的影响系数法，影响系数余量法，定相位角最佳算法，单平 面加重的优化搜索法及优化影响系数 的算法及其 MATLAB语言的实现。数值算例验证 了各 程序的正确与有效，同时说 明通过计算得到的多种配重方案 ，给平衡工作者提供 了宽广的选 择空间，有利于提高平衡工作的质量、效率和降低成本。 关键词 ：动平衡 影响系数法 影响系数余量法 优化 0 前 言 为影响系数本身，形成非线性约束优化问 题 ，并用罚函数法求解出更准确的影响系 计算机辅助转子动平衡常用的影响系 数。 数法 文【 1】是基于使平衡后的残余振动为 上述算法都有一定的应用业绩，值得 零 ，所以求出的配重往往过大 ，且有时也会 我们在平衡工作 中引用。为此 ，我们采用功 出现各修正面配重 的效应相互抵消 的现 能强大、使用方便的MATLAB 自行编制了 象。在许多试 图对其改进的努力中，文【 2】 这几个程序。经若干算例验证，不但算得 讨论了允许残余振动不为零 的 “影响系数 快，计算结果与上述文献吻合，且计算精度 余量法”，并对此非线性约束规划问题用 更高。 Rokafeller乘子法求解；并指出它能给有残 余弯曲的转子的平衡一个较好的结果。文【 1各算法的MATLAB语言表述 3】及 文【 6】的作者从避免各平面配重的效 应相互抵消出发 ，认为各平面配重的相位 1．1影响系数法 (最小二乘法 。Goodman 角以按单平面平衡计算所得到的相位角为 法) 最佳 ；配重的大小再用最小二乘法得到。本 转子的动平衡方程为 文 中称之为 “定相位角最佳算法”。对于只 KW 4-V0=0 (1) 用一个平衡面平衡多转速多测点的转子 ， 其中：K为影响系数矩阵，W 为待求配重 如燃机转子等，文【4】给了一个单平面最佳 列向量，V0为初始振动列向量，三者均为 配重的优化搜索算法。 复数。 由于各转子 的个性表现使同类转子的 于是W=K ＼(一V0) (2) 影响系数会有较大的分散度，往往使借用 “＼”表示左除，即，a~b=inv(a)*b，当方 同类转子的影响系数所作的平衡计算给出 程数多于未知数个数时，则执行最小二乘 不适用的结果 ，甚至于使得有人对这种作 算法。 法持否定态度。但是，毕竟 “有胜于无 ”：减 加配重后的残余振动Vcy为 ： 少些盲 目性，减少一两次平衡 中的开机成 Vcy=K W 4-V0 (3) 本。所以，文【5】把优化的目标函数直接定 自编程序文件名为 phen1．m。 8 东方汽轮机 2003年第 1期 到的相位角。其平衡方程为： 1．2影响系数余量法 K [Wx