[发明专利]一种滑动轴承-转子系统不平衡量识别方法

权利要求书

1.一种滑动轴承-转子系统不平衡量识别方法，其特征在于，包括以下步骤：

步骤S1、获取滑动轴承-转子系统不平衡量先验分布空间、初始样本及滑动轴承-转子系统的测试不平衡响应；

步骤S2，利用贝叶斯理论求解代价函数，并获得不平衡量的后验联合概率密度分布；

步骤S3，计算上述代价函数的最小值，当所述最小值大于收敛精度时，更新所述滑动轴承-转子系统不平衡量先验分布空间，并利用遗传智能采样技术获取新样本，进入步骤S2，否则，进入步骤S4；所述遗传智能采样技术是将旧样本遗传到新样本中；

步骤S4，采用MCMC方法识别得到滑动轴承-转子系统不平衡量的近似边缘后验概率密度分布，进而确定滑动轴承-转子系统不平衡量的均值、置信区间。

2.根据权利要求1所述的方法，其特征在于，所述利用遗传智能采样技术获取新样本，为通过遗传拉丁超立方实验设计更新先验分布空间中的新样本，包括：

根据极大极小距离准则，对落入到下代先验分布空间的样本进行筛选；

采用模拟火优化求解器使新样本在各个设计变量上投影均匀且产生的新样本到遗传样本的距离最大。

3.根据权利要求1或2所述的方法，其特征在于，所述求解代价函数，包括：

采用传递矩阵法将不平衡量作为输入变量通过滑动轴承-转子系统不平衡响应仿真分析，求取滑动轴承-转子系统不平衡响应；

将上述滑动轴承-转子系统的测试不平衡响应与求取的不平衡响应之间的误差作为代价函数。

4.根据权利要求3所述的方法，其特征在于，所述更新先验分布空间，为根据信赖域更新方法获取下代先验分布空间。

5.根据权利要求4所述的方法，其特征在于，所述信赖域更新方法为：通过黄金分割法自适应调整不平衡量先验分布空间的上下边界X_L和X_R，进而获得下代采样点的先验分布空间。

6.根据权利要求1所述的方法，其特征在于，所述获取不平衡量先验分布包括：

设定误差常数ε，给定滑动轴承-转子系统不平衡量先验均匀分布的初始上下边界和给定每一迭代步样本点数量N_sⁱ。

7.根据权利要求6所述的方法，其特征在于，所述不平衡量的后验联合概率密度分布p(X|q)为：

式中，X为滑动轴承-转子系统不平衡量；q为滑动轴承-转子系统的测试不平衡响应；a为常数；g(X)为代价函数；σ²为滑动轴承-转子系统随机噪声的方差。

8.根据权利要求7所述的方法，其特征在于，代价函数计算公式为：

式中，X为滑动轴承-转子系统不平衡量，h为测点个数，q为滑动轴承-转子系统的测试不平衡响应，q_j(X)为不平衡量作为输入变量通过滑动轴承-转子系统不平衡响应仿真分析获得的不平衡响应。

9.根据权利要求5所述的方法，其特征在于，所述获得下代采样点的先验分布空间为：

式中，1-α为黄金分割点，α＝0.382；和为第i代采样先验分布空间的上、下边界；和为第i+1代采样先验分布空间的上、下边界。

10.根据权利要求9所述的方法，其特征在于，在MCMC方法中，链长取10⁵量级。