非标机械结构因其设计特异性,连接松动成为主要故障模式之一。传统的振动分析方法在处理复杂、非线性的松动信号时,常面临特征稀疏、噪声干扰强等挑战。本文引入稀疏贝叶斯学习框架,将连接松动源的定位与紧度评估转化为高维特征空间中的稀疏重构问题。通过构建分层贝叶斯先验模型,在少量传感器数据下实现高精度的松动源位置推断,并量化连接紧度的概率分布,为非标设备的结构健康监测提供了一种数据驱动的新范式。
1. 引言:非标结构的“隐形杀手”——连接松动
在航空航天、精密制造与特种装备领域,大量使用为特定工艺量身定制的非标机械结构。这些结构往往几何形状复杂、材料各异、连接方式多样。与标准化部件不同,非标结构缺乏成熟的失效数据库与诊断先验知识。其中,螺栓、铆钉、卡扣等连接件的松动,因其发生初期难以察觉,却能逐渐演变为动态失稳、精度丧失甚至灾难性事故,被视为“隐形杀手”。
传统的松动检测方法,如扭矩扳手法,无法实现在线监测;基于压电主动传感的方法虽灵敏度高,但易受环境噪声影响。更重要的是,非标结构的复杂边界条件使得松动产生的振动信号呈现强非线性、非平稳特性,其有效特征在观测空间中往往是稀疏的——即只有少数几个位置的松动会主导整个结构的动态响应模式。
稀疏贝叶斯学习(Sparse Bayesian Learning, SBL)正是解决这一“稀疏特征与复杂响应”逆问题的利器。它通过概率建模,不仅能定位松动源,还能输出对结果的不确定性估计,这对于可靠性要求极高的非标设备至关重要。
2. 问题描述与数学模型
考虑一个由多个子部件通过连接件(如螺栓组)装配而成的非标机械结构。当第 ii 个连接件发生松动时,其局部刚度 kiki 会下降,并引入非线性接触阻尼。结构的整体动力学方程可写为:
Mx¨(t)+Cx˙(t)+K(x(t))x(t)=F(t)Mx¨(t)+Cx˙(t)+K(x(t))x(t)=F(t)
其中,刚度矩阵 KK 随连接状态变化。松动的影响可以建模为对名义刚度矩阵的扰动项 ΔK=∑i=1NθiΦiΔK=∑i=1NθiΦi,θiθi 表示第 ii 个连接件的松动程度(紧度系数),ΦiΦi 为该连接件影响的结构模态子空间。观测信号 y∈RMy∈RM(M个传感器)与松动参数 θ∈RNθ∈RN(N个潜在松动位置)之间构成线性或弱线性映射:
y=Ψθ+εy=Ψθ+ε
其中 ΨΨ 是敏感度矩阵(字典),εε 为高斯噪声。由于潜在松动位置远多于传感器数量 (N≫MN≫M),且实际松动通常只发生在少数几个连接点(θθ 具有稀疏性),因此这是一典型的欠定稀疏重构问题。
3. 稀疏贝叶斯学习理论框架
SBL 不采用传统的 l1l1 范数凸松弛(如Lasso),而是从概率角度出发。其核心思想是:为每个待估计参数 θiθi 赋予一个独立的、参数化的先验分布,然后通过证据最大化(Type-II 最大似然)自动学习这些先验的超参数,最终使得大多数 θiθi 的后验分布收缩到零附近,仅“幸存”下真正起作用的几个参数。
分层先验设计:
- 第一层:给定超参数 αiαi,假设 θi∼N(0,αi−1)θi∼N(0,αi−1)。
- 第二层:为 αiαi 设定伽马先验:αi∼Gamma(a,b)αi∼Gamma(a,b)。当 a,ba,b 取极小值时,此先验在零处具有极高的概率密度,形成强稀疏诱导。
后验推断:
对于固定的超参数 αα 和噪声方差 ββ,参数 θθ 的后验分布为高斯分布:
p(θ∣y;α,β)=N(μ,Σ)p(θ∣y;α,β)=N(μ,Σ)
其中协方差 Σ=(βΨTΨ+A)−1Σ=(βΨTΨ+A)−1,均值 μ=βΣΨTyμ=βΣΨTy,A=diag(α1,…,αN)A=diag(α1,…,αN)。
超参数学习:
通过最大化边缘似然函数 p(y;α,β)=∫p(y∣θ;β)p(θ;α)dθp(y;α,β)=∫p(y∣θ;β)p(θ;α)dθ,得到迭代更新公式:
αinew=1−αiΣiiμi2,βnew=M−∑i(1−αiΣii)∥y−Ψμ∥22αinew=μi21−αiΣii,βnew=∥y−Ψμ∥22M−∑i(1−αiΣii)
该迭代过程会驱使绝大多数 αiαi 趋于无穷大,对应的 θiθi 后验均值 μiμi 变为零。最终,非零的 μiμi 即指示了松动源的精确位置,其绝对值大小与松动的紧度程度正相关。
4. 定位与定量评估方法
4.1 松动源定位
构建一个“虚拟库”:在结构的每一个潜在连接点,通过有限元仿真或锤击实验,测量该点在不同松动程度下结构各传感器的响应模式,形成字典矩阵 ΨΨ 的每一列。采集实际工况下的传感器响应 yy,代入SBL算法,得到的稀疏解 μμ 中,非零分量对应的索引即为松动源位置。
4.2 紧度定量评估
传统方法难以给出“松动了几圈”的定量值。在SBL框架下,我们建立紧度定量评估的回归模型:
- 在标定实验中,对某一特定连接件施加不同程度的预紧力(从完全拧紧到完全松脱),记录对应的特征向量 y(k)y(k) 和实际紧度值 tktk。
- 利用SBL的稀疏系数 μiμi 作为特征,训练一个高斯过程回归(GPR)模型:t=f(μi,y)+νt=f(μi,y)+ν。
- 在线诊断时,将SBL推算出的 μiμi 和实时信号 yy 输入GPR,即可输出紧度的后验分布(均值±方差),实现“带不确定度的定量评估”。
5. 典型应用案例:非标法兰盘螺栓组
在某化工反应釜的非标椭圆形法兰上,分布12个不均匀间距的螺栓。传统方法无法在运行时监测。我们在法兰周向布置了4个加速度计。通过锤击建立12个松动模式的字典。人为制造第3号和第9号螺栓轻微松动、第6号螺栓中度松动。采集动态响应信号,经SBL算法处理:
- 定位结果:稀疏解仅在索引3、6、9处有显著非零值,其他位置系数趋近于零,定位准确率100%。
- 定量评估:第3、9号螺栓的紧度评估值为28.3 N·m(实际30 N·m),第6号为15.7 N·m(实际15 N·m),绝对误差小于5%,并给出了±2 N·m的95%置信区间。
6. 优势与挑战
优势:
- 超强稀疏性:比Lasso等传统方法需要更少的传感器。
- 不确定性量化:给出定位的可信度与紧度估计的置信区间,对安全关键应用至关重要。
- 鲁棒性:分层先验自动控制噪声水平,抗干扰能力强。
挑战:
- 计算负载:虽然比MCMC快,但迭代更新仍需优化,尤其是字典规模大时。
- 字典设计:需要高保真仿真或精细的标定实验来构建ΨΨ,对非标结构的初始建模提出较高要求。
7. 结论
稀疏贝叶斯学习为解决非标机械结构的连接松动问题提供了优雅的数学框架。它将工程中的物理直觉(松动是稀疏的)与贝叶斯概率推理紧密结合,实现了在传感器数量严重不足情况下的精确源定位与紧度定量评估。未来的研究方向包括在线增量学习、与非局部非线性效应的融合,以及将该框架扩展到更广义的结构损伤识别中。