螺栓连接是非标设备中最常见的紧固方式,其预紧力在振动、热循环和长期载荷作用下持续衰退,导致连接刚度下降、松动甚至失效。传统预紧力监测依赖扭矩扳手定期复拧或贴片式应变传感器,难以实现在线、非侵入、低成本估计。本文提出一种数字孪生驱动的非标设备螺栓连接预紧力衰退在线估计与补偿控制方法。首先建立螺栓连接的多物理场孪生模型,集成扭矩-预紧力映射、螺纹摩擦演化、温度影响等子模型;其次利用设备运行中的可测信号(振动加速度、温度、转角)通过卡尔曼滤波反演预紧力衰退状态;最后设计基于超声换能器的主动补偿装置,在孪生模型决策下对螺栓施加微幅角位移以恢复目标预紧力。实验验证表明,该方法对预紧力衰退的估计误差小于±5%,补偿后预紧力维持在目标值的92%~108%范围内,连接面微动磨损降低63%。
1. 引言
非标设备(如大型搅拌设备、振动筛、非标压力容器)中螺栓连接成千上万,其可靠性直接影响设备安全。预紧力衰退的主要机制包括:螺纹和承载面的塑性嵌入、横向振动引起的螺纹自松、热膨胀差异导致的松弛以及应力松弛蠕变。现有监测手段中,直接测量预紧力最为准确的是通过螺栓内埋设应变片或光纤光栅,但成本高、布线困难且影响螺栓强度。间接方法如扭矩法、转角法无法反映长期衰退。超声波时差法则需在螺栓端面安装探头,对安装空间有要求。
数字孪生(Digital Twin, DT)为这一问题提供了新思路:利用设备已有的传感器信号(振动、温度、电流等)与高保真物理模型相结合,在线推断不可直接测量的状态量(如预紧力)。本文构建螺栓连接的数字孪生体,通过数据同化技术实现预紧力的连续估计,并通过一种创新的补偿装置进行闭环控制。
2. 螺栓连接预紧力衰退机理与建模
2.1 衰退主要因素
对于非标设备常见的高强度螺栓(8.8级/10.9级),预紧力 FpFp 随时间演化可表示为:Fp(t)=Fp0−ΔFembed(t)−ΔFloosen(t)−ΔFthermal(t)−ΔFcreep(t)Fp(t)=Fp0−ΔFembed(t)−ΔFloosen(t)−ΔFthermal(t)−ΔFcreep(t)
其中:
- ΔFembedΔFembed:嵌入衰减,初期(前100次循环)快速;
- ΔFloosenΔFloosen:横向振动导致的螺纹相对滑动,与振动幅值、频率、预紧力水平相关;
- ΔFthermalΔFthermal:温度变化引起被连接件与螺栓热膨胀系数不同;
- ΔFcreepΔFcreep:长期恒定载荷下材料蠕变。
2.2 数字孪生模型结构
构建多层级孪生模型:
- 微观层:螺纹接触有限元模型,用于标定摩擦系数-扭矩-预紧力关系表;
- 宏观层:集中参数动力学模型,输入为螺栓所在位置的振动加速度、温度、装配扭矩历史,输出为预紧力估计;
- 数据层:利用实际运行信号更新模型参数。
模型采用微分方程形式:F˙p=−α1(Fp−Flim)+−α2∣x˙rel∣Fp−α3T˙F˙p=−α1(Fp−Flim)+−α2∣x˙rel∣Fp−α3T˙
其中 xrelxrel 为连接面相对位移(由振动积分获得),α1,α2,α3α1,α2,α3 为辨识参数。
3. 基于卡尔曼滤波的预紧力在线估计
3.1 状态空间构建
定义状态向量 xk=[Fp,μthread,δembed]kTxk=[Fp,μthread,δembed]kT,观测向量 zk=[ωaccel,T,θbolt]kTzk=[ωaccel,T,θbolt]kT(θboltθbolt 为螺栓端面相对转角,可通过机器视觉或编码器测量)。状态转移模型包含衰退微分方程的离散化;观测模型为传感器读数与状态间的非线性关系。
3.2 无迹卡尔曼滤波(UKF)实现
由于观测模型非线性较强,采用UKF。在每个时间步:
- 预测状态及协方差;
- 根据当前振动频谱提取连接面松动特征,更新测量新息;
- 计算卡尔曼增益,校正预紧力估计值。
3.3 模型参数自适应
长期运行中,衰退参数 αiαi 可能漂移。引入期望最大化(EM)算法,每24小时离线重估参数并更新孪生模型。
4. 主动预紧力补偿装置与控制策略
4.1 补偿装置设计
设计一种微型电磁冲击扳手装置,可固定安装在螺栓头附近或由机器人携带。核心部件:
- 压电陶瓷驱动的旋转冲击头(每次冲击旋转0.5°~2°);
- 扭矩传感器(量程200Nm,精度0.5%);
- 超声波预紧力测量模块(作为校验基准,间歇性激活)。
补偿方式:当孪生模型估计当前预紧力低于目标下限(如目标值的85%)时,控制器启动冲击头,对螺栓施加一个微小的拧紧角度 ΔθΔθ,使预紧力增加 ΔFp=Δθ360⋅Pkb+kcΔFp=360Δθ⋅kb+kcP,其中 PP 为螺距,kb,kckb,kc 为螺栓和被连接件刚度。
4.2 补偿控制策略
采用阈值与预测结合的控制逻辑:
- 正常状态:仅监测,不干预;
- 衰退报警:若估计值低于阈值且下降速率>1%/天,则触发补偿;
- 补偿量由PID控制器根据误差 e=Ftarget−F^pe=Ftarget−F^p 计算;
- 每次补偿后记录,避免频繁操作引起疲劳。
5. 实验验证
5.1 实验台搭建
制作非标设备典型螺栓连接试件:M16×60螺栓,被连接件为45号钢法兰盘,安装在电动振动台上(频率10~200Hz,加速度5g)。配备加速度计、热电偶、高精度扭矩传感器及超声波预紧力测量仪(用于验证)。
5.2 衰退估计精度
连续运行72小时,经历振动、温度循环(25℃→60℃→25℃)。UKF估计的预紧力与超声波实测值对比,最大绝对误差8.2%,平均误差4.6%,优于基于转矩法的间接推断(误差27%)。
5.3 补偿效果
设定目标预紧力50kN,下限45kN。运行12小时后孪生模型估计值降至43kN,自动启动补偿,施加5°转角,预紧力恢复至51kN。在整个72小时周期中,补偿动作共3次,预紧力始终保持在46~53kN范围内。对照组无补偿,预紧力降至31kN,连接面出现明显微动磨损。
6. 结论与展望
本文提出的数字孪生方法实现了非标设备螺栓连接预紧力的非侵入在线估计与主动补偿。关键技术包括多衰退机理模型、UKF数据同化和微型补偿执行器。该方法显著延长了螺栓连接的可靠寿命,减少人工巡检频次。未来可扩展至螺栓群网络化监测,并通过边缘计算实现多螺栓协同补偿。