(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 202210641952.6 (22)申请日 2022.06.07 (71)申请人 武汉科技大 学 地址 430000 湖北省武汉市和平大道 947号 (72)发明人 鲁志文　蒋阳　党章　袁锐　吕勇　 (74)专利代理 机构 北京哌智科创知识产权代理 事务所(普通 合伙) 11745 专利代理师 张元媛 (51)Int.Cl. G06K 9/00(2022.01) G06F 7/78(2006.01) G06F 17/16(2006.01) G06F 30/23(2020.01) G01D 21/02(2006.01) G06F 119/02(2020.01) (54)发明名称 一种超谐动力学模式驱动的转轴疲劳裂纹 定位方法 (57)摘要 本发明公开了一种超谐动力学模式驱动 的 转轴疲劳裂纹定位方法， 通过动力学模式分解提 取因疲劳裂纹呼吸作用非线性引起的超谐动力 学模式， 并将其转换到不同的尺度空间， 利用多 尺度超谐动力学模式信息融合， 形成了一种新的 转轴裂纹损伤指标， 通过损伤指标的峰值分布判 定裂纹的位置。 本发明仅需测试运行状态下的转 轴系统多测点振动响应， 无需健康 转轴的模型参 考， 也不需要测试健康状态下的振动信息， 且可 以避免转轴上的阶梯、 沟槽以及转子不对中等因 素的干扰， 具有噪声鲁棒性好、 快速可靠的优点， 适用于转轴 系统的在线监测。 权利要求书2页 说明书6页 附图4页 CN 114912494 A 2022.08.16 CN 114912494 A 1.一种超谐动力学模式驱动的转轴疲劳裂纹定位方法， 其特征在于， 所述方法包括如 下步骤： S1： 非接触测试平稳工况 下的转轴多测点横向振动响应； S2： 构造响应矩阵； S3： 通过动力学模式分解 提取超谐动力学模式 φ(x)； S4： 利用多尺度超谐模式证据融合构造损伤指标； S5： 判断各测点损伤指数分布， 若分布较均匀则不存在裂纹； S6： 否则， 转轴存在裂纹， 裂纹位置对应损伤指标峰值 位置。 2.如权利要求1所述的一种超谐动力学模式驱动的转轴疲劳裂纹定位方法， 其特征在 于， 步骤S2具体为： 将n个测点的振动时域响应组成响应矩阵Y 其中yi为测点i所测得的响应,m为响应数据中的点数。 3.如权利要求1所述的一种超谐动力学模式驱动的转轴疲劳裂纹定位方法， 其特征在 于， 步骤S3具体为： S3.1:取响应矩阵中的前l个数据点形成响应矩阵Y1及其时移 矩阵Y2； S3.2:假设Y1及Y2间存在线性变换， 使得： Y2＝AY1 则转换矩阵 其中 表示Y1的伪逆， 转换矩阵A的特征向量及特征值即为系统 响应矩阵对应的动力学分解模式及特征值， 由于A的维数很高， 以及系统噪声的干扰， 有必 要对转换矩阵A进行降维； S3.3:对Y1进行奇异值分解， 有： Y1＝UΣV* 式中， *表示矩阵的转置， 则有 因此， A＝Y2VΣ‑1U* S3.4:采用本征模式分解对A进行降维， 对其进行r阶秩截断， 则降维后的近似矩阵 即： S3.5:求解 的特征值及特 征向量即可 得系统的特 征值和动力学模式， 即， φ＝Uw 式中， μ为动力学模式分解特 征值， φ为对应的动力学模式； S3.6:根据转轴运行转速， 提取转频的高阶成分对应的超谐动力学模式即φ(x)。 4.如权利要求1所述的一种超谐动力学模式驱动的转轴疲劳裂纹定位方法， 其特征在 于， 步骤S4具体为：权　利　要　求　书 1/2 页 2 CN 114912494 A 2S4.1:将提取的超谐动力学模式φ(x)， 利用高斯多尺度分解方法将其分解到多个尺 度： 其中 为高斯核函数 S4.2:对各尺度下的超谐动力学模式计算 其Teager能量： T(x, σ )＝ L2(x, σ )‑L(x‑1, σ )L(x+1, σ ) S4.3:对各尺度Tea ger能量超谐动力学模式进行 数据融合， 具体操作如下： 在第k个尺度下， 裂纹出现在位置xi(i＝1,2…,n)处的概 率表示为： 因此， 第k个尺度下， 裂纹出现在x1～xn位置的概 率组成向量 利用D‑S证据融合规则， 对各尺度下的基本概率分配函数进行融合， 即可得到各个位置 出现裂纹的概率， 各尺度Teager能量超谐动力学模式对应数据源， 用S1到SM表示； 根据融合 规则， 前两个尺度融合后的结果表示 为： 其中K1表示证据1和证据2之间的冲突权 重： 依此规则， 对剩余的尺度按顺序进行融合， 最终可 得： 式中， 最后所得 即为损伤因子 。权　利　要　求　书 2/2 页 3 CN 114912494 A 3