一种基于数字孪生的控制阀流致振动疲劳寿命预测方法

1.本发明涉及阀门设备领域，是基于数字孪生技术对控制阀流致振动疲劳寿命进行预测技术，主要应用于电力、石油、化工的控制阀。


背景技术：

2.工业工艺出现了大量需要高压差控制类阀门控制的场合，如高温高压安全阀、疏水阀、止回阀、调节阀、泄放阀及溢流阀等。由于控制阀的流致振动疲劳损伤是影响使用寿命的最常见问题之一，仅通过通用疲劳寿命预测方法难以实现控制阀疲劳寿命的及时有效预测。这就需要及时监测控制类阀门流致振动疲劳损伤，对其疲劳寿命进行有效评估，以保证所设计的控制阀满足长寿命和安全可靠性的要求。
3.目前基于仿真模型和数字孪生技术对机械设备的疲劳寿命的预测方法，是针对机械振动工况下疲劳寿命预测，没有针对流体激振力对结构疲劳寿命影响的预测方法，并且没有有效考虑结构在实际工作过程中，工作环境变化等因素导致结构实测值相对于模拟预测值出现偏差的问题，不能实现控制阀的流致振动疲劳寿命的有效预测。


技术实现要素：

4.本发明的目的是提供一种基于数字孪生的控制阀流致振动疲劳寿命预测方法。
5.本发明是一种基于数字孪生的控制阀流致振动疲劳寿命预测方法，由以下步骤组成：步骤(1)利用压力传感器、加速度传感器及应变片，采集控制阀在工作工况下的脉动压力、振动加速度及应变信息，并对振动信号进行特征提取；步骤(2)将传感器监测获得的信号数据通过链路传输到处理器；步骤(3)建立数字孪生仿真模型，对控制阀在工作工况下的流致振动进行模拟分析，获得振动信息；步骤(4)根据传感器采集获得的振动信息及模拟分析获得的振动信息，建立漂移和扩散模型，生成监测数据与模拟数据的偏差数据集，并采用外推方法对偏差数据进行外推生成外推偏差数据集，修正并更新数字孪生仿真模型；步骤(5)基于miner-palmgren线性累积损伤理论，根据偏差数据预测控制阀的流致振动疲劳寿命。
6.本发明的有益效果是：本发明所述的基于数字孪生的控制阀流致振动疲劳寿命预测方法，解决了现有疲劳寿命预测方法难以实现控制阀疲劳寿命及时有效的预测的问题，尤其是工况及工作环境变化情况下。这种准确的基于数字孪生技术的流致振动疲劳寿命预测方法充分考虑控制阀实际工作过程中，由于工作环境变化等因素导致结构实测值相对于模拟预测值出现偏差的问题，采用漂移-扩散估计技术构建偏差模型，通过外推偏差数据的确定性分量和随机分量来生成外推偏差数据。根据外推偏差数据，通过数字孪生仿真模型及寿命预测模型预测控制阀的流致振动疲劳寿命，根据寿命结果对接近、已达到或超过规
定使用寿命的部件进行维护或更换。基于数字孪生技术的疲劳寿命预测系统通过对结构疲劳寿命及时有效的预测，有效降低了疲劳损伤引起的故障风险，解决了阀门实时健康监测维修管理的问题。
附图说明
7.图1是本发明的流程图，图2是多传感器测试流程图，图3是外推模块流程图。实施方式
8.如图1所示，本发明是一种基于数字孪生的控制阀流致振动疲劳寿命预测方法，由以下步骤组成：步骤(1)利用压力传感器、加速度传感器及应变片，采集控制阀在工作工况下的脉动压力、振动加速度及应变信息，并对振动信号进行特征提取；步骤(2)将传感器监测获得的信号数据通过链路传输到处理器；步骤(3)建立数字孪生仿真模型，对控制阀在工作工况下的流致振动进行模拟分析，获得振动信息；步骤(4)根据传感器采集获得的振动信息及模拟分析获得的振动信息，建立漂移和扩散模型，生成监测数据与模拟数据的偏差数据集，并采用外推方法对偏差数据进行外推生成外推偏差数据集，修正并更新数字孪生仿真模型；步骤(5)基于miner-palmgren线性累积损伤理论，根据外推数据预测控制阀的流致振动疲劳寿命。
9.以上所述方法，步骤(1)中所述，建立控制阀三维模型，采用ansys有限元分析软件对控制阀进行模态分析与静强度分析，根据计算获得的结构模态振型、固有频率及应力结果，确定传感器布置位置，测得测点位置处的脉动压力p
i1
、振动加速度a
i1
及应变ε
i1
信息，结合改进双参数阈值量化函数与小波阈值降噪算法，提取测点位置的振动信息，准确获得测点的振动信息。
10.以上所述方法，步骤(3)中所述，采用isight优化软件调用ansys有限元分析软件，对控制阀工作工况下的流致振动进行流固耦合数值模拟分析，以有限元分析中的输入变量，包括阀前压力p1、阀后压力p2、温度t以及阀门开度(li/l)参数变量，作为机器学习模型的输入，计算获得的振动加速度a
i2
以及应力时间历程σ
i2
作为模型训练的目标，获得训练样本数据；采用10％的数据集被选为训练数据集，5％的数据集被选为验证数据集，5％的数据集被选为测试数据集，采用非线性回归机器学习算法对分析获得的数据集进行训练、验证与测试，建立控制阀流致振动数字孪生仿真模型，对同步骤一相应工况下的控制阀流致振动进行模拟分析，获得控制阀测点位置处的振动加速度与应力时间历程数据。
11.以上所述方法，步骤(4)中所述，采用漂移-扩散估计方法，根据传感器采集获得的振动信息及模拟分析获得的振动加速度信息，对数字孪生模型模拟生成的结果与传感器监测获得的实际结果之间的偏差进行建模，生成漂移-扩散偏差模型；漂移-扩散模型将偏差动力学确定为两个分量，即确定性分量f(x，t)和随机分量g(x，t)，每个变量d(x)是确定性漂移分量及随机分量的组合，遵循微分方程(1)：d(x)＝f(x，t)dt+g(x，t)dw
t
(1)其中，x为控制阀振动信号偏差，w
t
为一个标准布朗运动过程，t是时间；f(x，t)是
palmgren线性累积损伤理论，计算在使用周期内的疲劳损伤总量其中，d为控制阀在使用周期内产生的疲劳损伤，ni为实际应力水平下的循环次数，ni为材料的s-n曲线在相应应力幅值下的循环次数。进而计算获得控制阀流致振动疲劳寿命其中，l为控制阀的流致振动疲劳寿命。根据寿命结果对接近、已达到或超过规定使用寿命的部件进行维护或更换。
13.下面用实施例进一步展开本发明的技术内容。如图1所示，本发明是一种基于数字孪生的控制阀流致振动疲劳寿命预测方法，其步骤为：步骤一、利用压力传感器、加速度传感器及应变片，采集控制阀在工作工况下的脉动压力、振动加速度及应变信息，并对振动信息进行特征提取和去噪处理；步骤二、将传感器监测获得的信号数据通过链路传输到处理器；步骤三、建立数字孪生仿真模型，并对控制阀在工作工况下的流致振动进行模拟分析，获得振动信息；步骤四、根据传感器采集获得的振动信息及模拟分析获得的振动信息，建立漂移和扩散模型，生成监测数据与模拟数据的偏差数据集，并采用外推方法对偏离数据进行外推生成外推偏差数据集；步骤五、基于miner-palmgren线性累积损伤理论，根据外推数据预测控制阀的流致振动疲劳寿命。
14.其中在上述步骤一中，如图2首先建立控制阀三维模型，采用ansys有限元分析软件对控制阀进行模态分析与静强度分析，获得结构的模态特性、模态振型、固有频率及应力结果等，根据分析得到的结构模态结果及应力结果确定传感器布置位置。
15.将多个传感器即压力传感器、加速度传感器及应变传感器固定安装在确定的传感器布置位置，测得测点位置处的脉动压力p
i1
、振动加速度a
i1
及应变ε
i1
信息，包含周期性更新的数据，即正常工作工况下监测的数据，以及突发性更新的数据，如阀门关闭及开启时测得的数据。
16.基于时频谱特征对传感器提取的测点位置的振动信号进行标定。结合改进双参数阈值量化函数与小波阈值降噪算法，基于标定后的特征频段，提取测点的振动信号，准确获得控制阀测点处的振动信息，为数字孪生仿真模型的正确建立与修正提供数据支撑。
17.其中在上述步骤二中，将传感器监测获得的参数通过链路传输到处理器，链路包括有线电缆、无线连接或者光缆技术。
18.其中在上述步骤三中，采用isight优化软件调用ansys有限元分析软件，对控制阀工作工况下的流致振动进行流固耦合数值模拟分析，以有限元分析中的输入变量，包括阀前压力p1、阀后压力p2、温度t以及阀门开度(li/l)参数变量，作为机器学习模型的输入，计算获得的振动加速度a
i2
以及应力时间历程σ
i2
作为模型训练的目标，获得训练样本数据。其中10％的数据集被选为训练数据集，5％的数据集被选为验证数据集，5％的数据集被选为测试数据集，采用非线性回归机器学习算法对分析获得的数据集进行训练、验证与测试，建立控制阀流致振动数字孪生仿真模型，对同步骤一相应工况下的控制阀流致振动进行模拟分析，获得控制阀测点处的振动加速度与应力时间历程数据。
19.其中在上述步骤四中，根据传感器采集以及模拟分析获得的振动加速度信息，采用漂移-扩散估计方法，对数字孪生模型模拟生成的结果与传感器监测获得的实际结果之间的偏差分别进行建模。漂移-扩散模型将偏差动力学确定为两个分量，给定随时间变化的
偏差度量，则偏差动力学的两个分量的函数分别为确定性分量f(x，t)和随机分量g(x，t)，每个变量d(x)看作是确定性漂移分量及随机分量的组合，遵循微分方程(1)：
20.d(x)＝f(x，t)dt+g(x，t)dw
t
(1)
21.其中，x为控制阀振动信号偏差，w
t
为一个标准布朗运动过程，是混沌运动产生的不稳定变量，t是时间；f(x，t)是确定性分量，g(x，t)是随机分量，随机分量服从正态分布dw
t
～n(0，1)。
22.采用基于核函数的回归算法，从偏差数据估计f(x，t)和g(x，t)，假设动力学不会随时间变化，则f(x，t)和g(x，t)简化为f(x)和g(x)。基于核函数的回归算法可以采用nadaraya-watson核回归方法、priestley-chao核回归估计方法或gasser-muller核估计方法来实现。
23.对历史数据包含的数据，即当变量xi∈[x-，x
+
]，x-是左端点的位置，x
+
是右端点的位置，通过三次样条差值拟合函数，如公式(2)所示，预测历史数据范围内的偏差数据：
[0024][0025]
其中，系数ai，bi，ci，di分别为ai＝yi；；hi为步长hi＝x
i+1-xi。
[0026]
历史数据范围之外的数据，采用三种不同的外推函数对历史数据进行外推。控制阀流致振动数据信息满足平稳随机性，假设外推函数可以避免人为地创建不动点，且在预测中可以避免出现过早发散。根据假设，f(x)的外推函数分为三种情况。在第一种情况下，外推继续遵循现有的曲率，但会衰减以避免变形。在第二种情况下，外推遵循沿轴线f(x)＝0的渐近线，第三种情况是沿f(x)＝αx的渐近线，其中系数α≠0。
[0027]
第一种情况通过公式(3)拟合来满足，如公式(3)所示：
[0028][0029]
其中，x
+
是右端点的位置，左端点使用f(-x)和-x
+
进行外推。
[0030]
在第二种情况下，外推遵循沿轴线f(x)＝0的渐近线，通过龚伯兹曲线进行拟合，公式(4)如下所示：
[0031][0032]
其中，k为待定参数，f
xi
″
ex
+
＞1，0＜f
xi
′‑fxi
′
x
+
＜1。
[0033]
第三种情况适用于以下情况：f(x
+
)＜0和f
′
(x
+
)＜0，或f(x
+
)＞0和f
′
(x
+
)＞0表示左端点；f(x
+
)＞0和f
′
(x
+
)＜0，或f(x
+
)＜0和f
′
(x
+
)＞0表示右端点。在这种情况下，可以根据公式(5)进行简单衰减：
[0034][0035]
以上三种情况适用于以下规则，如附图3所示，对于左端点，如果sgn(f)≠sgn(f
′
)且sgn(f)＝sgn(f
″
)，则用公式(5)。不满足时，如果sgn(f)＝sgn(f
′
)，则用公式(4)。对于右端点，如果sgn(f)＝sgn(f
′
)且sgn(f)≠sgn(f
″
)，则用公式(5)，不满足时，若sgn(f)＝sgn(f
′
)，则用公式(3)。否则，在其他情况下使用公式(3)，其中“sgn”为符号函数。
[0036]
根据外推的偏差数据和模拟获得的数据的组合生成更新的偏差数据。利用更新的
偏差数据进而更新改进数字孪生仿真模型，并利用改进的数字孪生仿真模型进行控制阀的流致振动模拟分析，将更新的数值模拟信息传递至疲劳预测模型。
[0037]
其中在上述步骤五中，根据模拟计算获得的控制阀测点应力谱，采用雨流计数法统计控制阀测点位置应力循环的应力幅值σ
fi
、平均应力σm和应力循环次数ni。根据控制阀零部件材料的s-n曲线获得各个应力幅值作用下的循环次数ni，基于miner-palmgren线性累积损伤理论，计算在使用周期内的疲劳损伤总量其中，d为控制阀在使用周期内产生的疲劳损伤，ni为实际应力水平下的循环次数，ni为材料的s-n曲线在相应应力幅值下的循环次数。进而计算获得控制阀流致振动疲劳寿命其中，l为控制阀的流致振动疲劳寿命次数。根据寿命结果对接近、已达到或超过规定使用寿命的部件进行维护或更换。

技术特征：
1.一种基于数字孪生的控制阀流致振动疲劳寿命预测方法，其特征在于，由以下步骤组成：步骤(1)利用压力传感器、加速度传感器及应变片，采集控制阀在工作工况下的脉动压力、振动加速度及应变信息，并对振动信号进行特征提取；步骤(2)将传感器监测获得的信号数据通过链路传输到处理器；步骤(3)建立数字孪生仿真模型，对控制阀在工作工况下的流致振动进行模拟分析，获得振动信息；步骤(4)根据传感器采集获得的振动信息及模拟分析获得的振动信息，建立漂移和扩散模型，生成监测数据与模拟数据的偏差数据集，并采用外推方法对偏差数据进行外推生成外推偏差数据集，修正并更新数字孪生仿真模型；步骤(5)基于miner-palmgren线性累积损伤理论，根据外推数据预测控制阀的流致振动疲劳寿命。2.根据权利要求1所述的基于数字孪生的控制阀流致振动疲劳寿命预测方法，其特征在于：步骤(1)中所述，建立控制阀三维模型，采用ansys有限元分析软件对控制阀进行模态分析与静强度分析，根据计算获得的结构模态振型、固有频率及应力结果，确定传感器布置位置，测得测点位置处的脉动压力p
i1
、振动加速度a
i1
及应变ε
i1
信息，结合改进双参数阈值量化函数与小波阈值降噪算法，提取测点位置的振动信息，准确获得测点的振动信息。3.根据权利要求1所述的基于数字孪生的控制阀流致振动疲劳寿命预测方法，其特征在于：步骤(3)中所述，采用isight优化软件调用ansys有限元分析软件，对控制阀工作工况下的流致振动进行流固耦合数值模拟分析，以有限元分析中的输入变量，包括阀前压力p1、阀后压力p2、温度t以及阀门开度(l
i
/l)参数变量，作为机器学习模型的输入，计算获得的振动加速度a
i2
以及应力时间历程σ
i2
作为模型训练的目标，获得训练样本数据；采用10％的数据集被选为训练数据集，5％的数据集被选为验证数据集，5％的数据集被选为测试数据集，采用非线性回归机器学习算法对分析获得的数据集进行训练、验证与测试，建立控制阀流致振动数字孪生仿真模型，对同步骤一相应工况下的控制阀流致振动进行模拟分析，获得控制阀测点位置处的振动加速度与应力时间历程数据。4.根据权利要求1所述的基于数字孪生的控制阀流致振动疲劳寿命预测方法，其特征在于：步骤(4)中所述，采用漂移-扩散估计方法，根据传感器采集获得的振动信息及模拟分析获得的振动加速度信息，对数字孪生模型模拟生成的结果与传感器监测获得的实际结果之间的偏差进行建模，生成漂移-扩散偏差模型；漂移-扩散模型将偏差动力学确定为两个分量，即确定性分量f(x，t)和随机分量g(x，t)，每个变量d(x)是确定性漂移分量及随机分量的组合，遵循微分方程(1)：d(x)＝f(x，t)dt+g(x，t)dw
t
ꢀꢀꢀ
(1)其中，x为控制阀振动信号偏差，w
t
为一个标准布朗运动过程，t是时间；f(x，t)是确定性分量，g(x，t)是随机分量，随机分量服从正态分布dw
t
～n(0，1)；采用基于核函数的回归算法，从偏差数据估计f(x，t)和g(x，t)，假设动力学不会随时间变化，则f(x，t)和g(x，t)简化为f(x)和g(x)，基于核函数的回归算法采用nadaraya-watson核回归方法、priestley-chao核回归估计方法或gasser-muller核估计方法来实现；对历史数据包含的数据，即当变量x
i
∈[x-，x
+
]，x-是左端点的位置，x
+
是右端点的位置，
通过三次样条差值拟合函数，如公式(2)所示，预测历史数据范围内的偏差数据：其中，系数a
i
，b
i
，c
i
，d
i
分别为a
i
＝y
i
；；h
i
为步长h
i
＝x
i+1-x
i
；历史数据范围之外的数据，采用三种不同的外推函数对历史数据进行外推，假设外推函数能避免人为地创建不动点，且在预测中能避免出现过早发散；根据假设，f(x)的外推函数分为三种情况，在第一种情况下，外推继续遵循现有的曲率，但会衰减以避免变形；在第二种情况下，外推遵循沿轴线f(x)＝0的渐近线，第三种情况是沿f(x)＝αx的渐近线，其中系数a≠0；第一种情况通过公式(3)拟合来满足，如公式(3)所示：其中，x
+
是右端点的位置，左端点使用f(-x)和-x
+
进行外推；在第二种情况下，外推遵循沿轴线f(x)＝0的渐近线，通过龚伯兹曲线进行拟合，公式(4)如下所示：其中，k为待定参数，f
xi
″
ex
+
＞1，0＜f
xi
′‑
f
xi
′
x
+
＜1；第三种情况适用于以下情况：f(x
+
)＜0和f
′
(x
+
)＜0，或f(x
+
)＞0和f
′
(x
+
)＞0表示左端点；f(x
+
)＞0和f
′
(x
+
)＜0，或f(x
+
)＜0和f
′
(x
+
)＞0表示右端点；在这种情况下，根据公式(5)进行简单衰减：以上三种情况适用于以下规则，对于左端点，如果sgn(f)≠sgn(f
′
)且sgn(f)＝sgn(f
″
)，则用公式(5)；不满足时，如果sgn(f)＝sgn(f
′
)，则用公式(4)；对于右端点，如果sgn(f)＝sgn(f
′
)且sgn(f)≠sgn(f
″
)，则用公式(5)；不满足时，若sgn(f)＝sgn(f
′
)，则用公式(4)，否则，在其他情况下使用公式(3)，其中“sgn”为符号函数；根据外推的偏差数据和模拟获得的数据的组合生成更新的偏差数据，利用更新的偏差数据进而更新改进数字孪生仿真模型，并利用改进的数字孪生模型进行控制阀的流致振动模拟分析，将更新的数值模拟信息传递至疲劳预测模型。5.根据权利要求1所述的基于数字孪生的控制阀流致振动疲劳寿命预测方法，其特征在于：步骤(5)中所述，根据模拟计算获得的控制阀测点应力谱，采用雨流计数法统计控制阀测点位置应力循环的应力幅值σ
fi
、平均应力σ
m
和应力循环次数n
i
，根据控制阀零部件材料的s-n曲线获得各个应力幅值作用下的循环次数n
i
，基于miner-palmgren线性累积损伤理论，计算在使用周期内的疲劳损伤总量其中，d为控制阀在使用周期内产生的疲劳损伤，n
i
为实际应力水平下的循环次数，n
i
为材料的s-n曲线在相应应力幅值下的循环次数；进而计算获得控制阀流致振动疲劳寿命其中，l为控制阀的流致振动疲劳寿
命；根据寿命结果对接近、已达到或超过规定使用寿命的部件进行维护或更换。

技术总结
一种基于数字孪生的控制阀流致振动疲劳寿命的预测方法，首先采用多传感器采集控制阀在工作工况下的振动信息，并进行特征提取。其次建立数字孪生仿真模型进行模拟分析，进而建立漂移和扩散模型，生成监测数据与模拟数据的偏差数据，并对偏差数据进行外推生成外推偏差数据集，修正并更新数字孪生仿真模型，最终基于Miner-Palmgren线性累积损伤理论，根据外推数据预测控制阀的流致振动疲劳寿命；本发明提出的基于数字孪生的疲劳寿命预测方法，充分考虑控制阀在实际工作过程中，由于环境变化等因素导致结构实测值相对于模拟预测值出现偏差的问题，实现了对控制阀流致振动疲劳寿命及时有效的预测，有效降低了疲劳损伤引起的故障风险。险。险。


技术研发人员：杨玲霞 李树勋 侯建军 郑铭兴 韦聪 张钰 邓国龙 赵倩 刘丹 雷志军 张建正 张博浩
受保护的技术使用者：兰州理工大学
技术研发日：2023.05.30
技术公布日：2023/8/13