CN110954349B - 一种基于残余应力畸变率的起重机结构健康状态监测方法 - Google Patents
一种基于残余应力畸变率的起重机结构健康状态监测方法 Download PDFInfo
- Publication number
- CN110954349B CN110954349B CN201911196181.9A CN201911196181A CN110954349B CN 110954349 B CN110954349 B CN 110954349B CN 201911196181 A CN201911196181 A CN 201911196181A CN 110954349 B CN110954349 B CN 110954349B
- Authority
- CN
- China
- Prior art keywords
- residual stress
- value
- ith
- area
- fixed time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000036541 health Effects 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000012544 monitoring process Methods 0.000 title claims abstract description 16
- 239000013078 crystal Substances 0.000 claims description 30
- 230000006978 adaptation Effects 0.000 claims description 29
- 230000002194 synthesizing effect Effects 0.000 claims description 3
- 238000003466 welding Methods 0.000 claims description 2
- 239000000956 alloy Substances 0.000 claims 1
- 229910045601 alloy Inorganic materials 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 4
- 238000011156 evaluation Methods 0.000 abstract 1
- 230000008859 change Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000003862 health status Effects 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
- G01M99/007—Subject matter not provided for in other groups of this subclass by applying a load, e.g. for resistance or wear testing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L1/00—Measuring force or stress, in general
- G01L1/25—Measuring force or stress, in general using wave or particle radiation, e.g. X-rays, microwaves, neutrons
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
- G01L5/0047—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes measuring forces due to residual stresses
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Toxicology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明公开了一种基于残余应力畸变率的起重机结构健康状态监测方法,包括如下步骤:S1、确定起重机结构易损伤的区域;S2、基于X射线的关键区域残余应力实时采集;S3、固定时间节点的残余应力密集度计算;S4、固定时间节点的残余应力畸变率计算;S5、起重机结构健康状态判量的计算。该方法检测精度高,可实时、高效地实现起重机结构健康状态的评估。
Description
技术领域
本发明涉及起重机结构状态监测方法,特别涉及一种基于残余应力畸变率的起重机结构健康状态监测方法。
背景技术
当前工业生产中安全生产问题的越来越被强化,起重机的运行安全性问题越来越受到社会的关注。起重机因其吨位大,工作环境恶劣,微小的故障就有可能会造成巨大的伤害,所以对起重机结构的安全状态进行实时监测是必要且紧迫的。起重机的健康状态监测方法有很多种,但大多集中在通过单一应力与损伤力学模型结合来实现整体起重机健康状态的评估,但由于实际工况下的起重机损伤机理比较复杂,单一通过宏观指标来进行安全监测会很多局限性,准确性不高,检测效率低。如何从更细尺度的内部残余应力出发来对起重机健康状态进行监测是目前这一研究领域存在的共性问题。
发明内容
发明目的:本发明目的是提供一种检测精度高的基于残余应力畸变率的起重机结构健康状态监测方法。
技术方案:本发明提供一种基于残余应力畸变率的起重机结构健康状态监测方法,包括如下步骤:
S1、确定起重机结构易损伤的区域;
S2、基于X射线的关键区域残余应力实时采集;
S3、固定时间节点的残余应力密集度计算;
S4、固定时间节点的残余应力畸变率计算;
S5、健康状态判量的计算。
进一步地,所述步骤S1中确定起重机结构易损伤的区域后,对相应的焊缝结构区域进行标记。
进一步地,所述步骤S2中在易损伤的区域布置X射线残余应力测量装置,分别实时采集每个损伤区域的X、Y、Z三个方向的残余应力值,即:σriqx、σriqy、σriqz,
然后按式(1)对三向残余应力进行合成,即可得到每个损伤区域的实时综合残余应力值σriq,
其中,σriq为在q时刻第i个损伤区域所对应的综合残余应力值;σriqx为在q时刻第i个损伤区域所对应的X方向残余应力值;σriqy为在q时刻第i个损伤区域所对应的Y方向残余应力值;σriqz为在q时刻第i个损伤区域所对应的Z方向残余应力值;α为起重机结构所对应的X向残余应力晶向适配系数,取值区间为1.035-1.15;β为起重机结构所对应的Y向残余应力晶向适配系数,取值区间为0.855-0.95;γ为起重机结构所对应的Z向残余应力晶向适配系数,取值区间为0.81-0.9;i为损伤区域的编号,i=1.2.3....N,N为损伤区域的全部个数。
其中,tm为固定时间节点的时间值;m为时间节点编号,取值区间为1,2,3…;为固定时间节点tm时第i个损伤区域的残余应力密集度;为在固定时间节点tm时第i个损伤区域所对应的综合残余应力值;q为时间;α为起重机结构所对应的X向残余应力晶向适配系数,取值区间为1.035-1.15;β为起重机结构所对应的Y向残余应力晶向适配系数,取值区间为0.855-0.95;γ为起重机结构所对应的Z向残余应力晶向适配系数,取值区间为0.81-0.9;S(σriq)max为0-tm时刻内第i个损伤区域采集到的所有综合残余应力值中的最大值;S(σriq)min为0-tm时刻内第i个损伤区域采集到的所有综合残余应力值中的最小值;A为残余应力采集面积系数,取值为2。
进一步地,所述步骤S4中根据各损伤区域的残余应力密集度,代入下式对固定时间节点tm的残余应力畸变率进行计算,
其中,为固定时间节点tm时第i个损伤区域的残余应力畸变率;为固定时间节点tm时第i个损伤区域的的残余应力密集度值;为在固定时间节点tm时第i个损伤区域所对应的综合残余应力值;α为起重机结构所对应的X向残余应力晶向适配系数,取值区间为1.035-1.15;β为起重机结构所对应的Y向残余应力晶向适配系数,取值区间为0.855-0.95;γ为起重机结构所对应的Z向残余应力晶向适配系数,取值区间为0.81-0.9;S(σriq)max为0-tm时刻内第i个损伤区域采集到的所有综合残余应力值中的最大值;S(σriq)min为0-tm时刻内第i个损伤区域采集到的所有综合残余应力值中的最小值;为0-tm时刻内第i个损伤区域采集到的所有综合残余应力值中的平均值,
其中,H为固定时间节点tm时整体起重机结构的健康状态判量;Hi为固定时间节点tm时起重机第i个损伤区域的健康状态判量;S(Hi)max为固定时间节点tm时所有损伤区域的健康状态判量中的最大值;S(Hi)min为固定时间节点tm时所有损伤区域的健康状态判量中的最小值;N为损伤区域的个数;为固定时间节点tm时第i个损伤区域的残余应力畸变率;为固定时间节点tm时第i个损伤区域的的残余应力密集度值;为0-tm时刻内第i个损伤区域的残余应力畸变率的最小值;为0-tm时刻内第i个损伤区域的残余应力畸变率的最大值;A为残余应力采集面积系数,取值为2。
以上α为起重机结构所对应的X向残余应力晶向适配系数,取值优选为1.15;β为起重机结构所对应的Y向残余应力晶向适配系数,取值优选为0.95;γ为起重机结构所对应的Z向残余应力晶向适配系数,取值优选为0.9。
有益效果:本发明方法中的残余应力作为能够体现起重机内部能量的一项重要指标,其实时畸变率更能够体现起重机结构在外载作用下的内部能量变化过程,也更能够体现实时工况条件下起重机的健康状态。本发明方法可以实现工况条件下起重机结构健康状态的监测,并从微观残余应力实时递进效应的角度出发来对起重机内部能量的变化过程进行预测,实现起重机健康状态的参数监测,从而避免了通过传统阈值判据或断裂力学模型来评判整体起重机健康状态的局限性,更有效地提高起重机健康状态检测的准确性与效率。
附图说明
图1为本发明方法的流程图。
具体实施方式
实施例1
如图1所示,本实施例的基于残余应力畸变率的起重机结构健康状态监测方法,包括以下步骤:
S1.确定起重机结构易损伤的区域:
根据经验确定起重机结构的常见易损伤区域,并在相应的焊缝结构区域进行标记,以便后续残余应力分布的实时采集;
S2.基于X射线的关键区域残余应力实时采集:
在S1的基础上,对常见损伤区域布置X射线残余应力测量装置来实时分别采集每个损伤区域的X、Y、Z三个方向的残余应力值,即:σriqx、σriqy、σriqz。
然后按式(1)对三向残余应力进行合成,即可得到每个损伤区域的实时综合残余应力值σriq。
其中,σriq为在q时刻第i个损伤区域所对应的综合残余应力值;σriqx为在q时刻第i个损伤区域所对应的X方向残余应力值;σriqy为在q时刻第i个损伤区域所对应的Y方向残余应力值;σriqz为在q时刻第i个损伤区域所对应的Z方向残余应力值;α为起重机结构所对应的X向残余应力晶向适配系数,取值区间为1.035-1.15;β为起重机结构所对应的Y向残余应力晶向适配系数,取值区间为0.855-0.95;γ为起重机结构所对应的Z向残余应力晶向适配系数,取值区间为0.81-0.9;i为损伤区域的编号,i=1.2.3....N,N为损伤区域的全部个数。
S3.固定时间节点的残余应力密集度计算:
其中,tm为固定时间节点的时间值;m为时间节点编号,取值区间为1,2,3…;为固定时间节点tm时第i个损伤区域的残余应力密集度;为在固定时间节点tm时第i个损伤区域所对应的综合残余应力值;q为时间;α为起重机结构所对应的X向残余应力晶向适配系数,取值区间为1.035-1.15;β为起重机结构所对应的Y向残余应力晶向适配系数,取值区间为0.855-0.95;γ为起重机结构所对应的Z向残余应力晶向适配系数,取值区间为0.81-0.9;S(σriq)max为0-tm时刻内第i个损伤区域采集到的所有综合残余应力值中的最大值;S(σriq)min为0-tm时刻内第i个损伤区域采集到的所有综合残余应力值中的最小值;A为残余应力采集面积系数,取值为2。
S4.固定时间节点的残余应力畸变率计算:
根据S3中各损伤区域的残余应力密集度,代入下式对固定时间节点tm的残余应力畸变率进行计算。
其中,为固定时间节点tm时第i个损伤区域的残余应力畸变率;为固定时间节点tm时第i个损伤区域的的残余应力密集度值;为在固定时间节点tm时第i个损伤区域所对应的综合残余应力值;α为起重机结构所对应的X向残余应力晶向适配系数,取值区间为1.035-1.15;β为起重机结构所对应的Y向残余应力晶向适配系数,取值区间为0.855-0.95;γ为起重机结构所对应的Z向残余应力晶向适配系数,取值区间为0.81-0.9;S(σriq)max为0-tm时刻内第i个损伤区域采集到的所有综合残余应力值中的最大值;S(σriq)min为0-tm时刻内第i个损伤区域采集到的所有综合残余应力值中的最小值;为0-tm时刻内第i个损伤区域采集到的所有综合残余应力值中的平均值,
S5.健康状态判量的计算:
Claims (3)
1.一种基于残余应力畸变率的起重机结构健康状态监测方法,其特征在于:包括如下步骤:
S1、确定起重机结构易损伤的区域;
S2、基于X射线的易损伤的区域残余应力实时采集;
S3、固定时间节点的残余应力密集度计算;
S4、固定时间节点的残余应力畸变率计算;
S5、健康状态判量的计算,
所述步骤S1中确定起重机结构易损伤的区域后,对相应的焊缝结构区域进行标记;
所述步骤S2中在易损伤的区域布置X射线残余应力测量装置,分别实时采集每个易损伤的区域的X、Y、Z三个方向的残余应力值,即:σriqx、σriqy、σriqz,
然后按式(1)对三向残余应力进行合成,即可得到每个易损伤的区域的实时综合残余应力值σriq,
其中,σriq为在q时刻第i个易损伤的区域所对应的综合残余应力值;σriqx为在q时刻第i个易损伤的区域所对应的X方向残余应力值;σriqy为在q时刻第i个易损伤的区域所对应的Y方向残余应力值;σriqz为在q时刻第i个易损伤的区域所对应的Z方向残余应力值;α为起重机结构所对应的X向残余应力晶向适配系数,取值区间为1.035-1.15;β为起重机结构所对应的Y向残余应力晶向适配系数,取值区间为0.855-0.95;γ为起重机结构所对应的Z向残余应力晶向适配系数,取值区间为0.81-0.9;i为易损伤的区域的编号,i=1.2.3…N,N为易损伤的区域的全部个数;
其中,tm为固定时间节点的时间值;m为时间节点编号,取值区间为1,2,3…;为固定时间节点tm时第i个易损伤的区域的残余应力密集度;为在固定时间节点tm时第i个易损伤的区域所对应的综合残余应力值;q为时间;α为起重机结构所对应的X向残余应力晶向适配系数,取值区间为1.035-1.15;β为起重机结构所对应的Y向残余应力晶向适配系数,取值区间为0.855-0.95;γ为起重机结构所对应的Z向残余应力晶向适配系数,取值区间为0.81-0.9;S(σriq)max为0-tm时刻内第i个易损伤的区域采集到的所有综合残余应力值中的最大值;S(σriq)min为0-tm时刻内第i个易损伤的区域采集到的所有综合残余应力值中的最小值;A为残余应力采集面积系数,取值为2。
2.根据权利要求1所述的基于残余应力畸变率的起重机结构健康状态监测方法,其特征在于:所述步骤S4中根据各易损伤的区域的残余应力密集度,代入下式对固定时间节点tm的残余应力畸变率进行计算,
其中,为固定时间节点tm时第i个易损伤的区域的残余应力畸变率;为固定时间节点tm时第i个易损伤的区域的残余应力密集度值;为在固定时间节点tm时第i个易损伤的区域所对应的综合残余应力值;α为起重机结构所对应的X向残余应力晶向适配系数,取值区间为1.035-1.15;β为起重机结构所对应的Y向残余应力晶向适配系数,取值区间为0.855-0.95;γ为起重机结构所对应的Z向残余应力晶向适配系数,取值区间为0.81-0.9;S(σriq)max为0-tm时刻内第i个易损伤的区域采集到的所有综合残余应力值中的最大值;S(σriq)min为0-tm时刻内第i个易损伤的区域采集到的所有综合残余应力值中的最小值;为0-tm时刻内第i个易损伤的区域采集到的所有综合残余应力值中的平均值,
3.根据权利要求1所述的基于残余应力畸变率的起重机结构健康状态监测方法,其特征在于:所述步骤S5中根据固定时间节点tm时第i个易损伤的区域的残余应力畸变率对在固定时间节点tm时整体起重机结构的健康状态判量H进行计算,
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911196181.9A CN110954349B (zh) | 2019-11-28 | 2019-11-28 | 一种基于残余应力畸变率的起重机结构健康状态监测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911196181.9A CN110954349B (zh) | 2019-11-28 | 2019-11-28 | 一种基于残余应力畸变率的起重机结构健康状态监测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110954349A CN110954349A (zh) | 2020-04-03 |
CN110954349B true CN110954349B (zh) | 2021-06-01 |
Family
ID=69978843
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911196181.9A Active CN110954349B (zh) | 2019-11-28 | 2019-11-28 | 一种基于残余应力畸变率的起重机结构健康状态监测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110954349B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113111564B (zh) * | 2021-05-24 | 2023-05-23 | 扬州大学 | 基于自适应预测区间内置式芦苇收割机健康状态评估方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102426068A (zh) * | 2011-09-14 | 2012-04-25 | 华东理工大学 | 一种表面薄膜残余应力的预测方法 |
CN102967397A (zh) * | 2012-11-19 | 2013-03-13 | 辽宁省电力有限公司电力科学研究院 | 焊接结构焊缝超声波残余应力测量装置及方法 |
CN104850691A (zh) * | 2015-05-05 | 2015-08-19 | 南京市特种设备安全监督检验研究院 | 一种基于多因素融合修正的结构件裂纹扩展预测方法 |
CN104848969A (zh) * | 2015-05-22 | 2015-08-19 | 华中科技大学 | 一种基于有限测试点的构件残余应力场的预测方法 |
CN106815419A (zh) * | 2017-01-03 | 2017-06-09 | 东南大学 | 一种基于裂纹信息预测的起重机运行状态在线评估方法 |
CN107186342A (zh) * | 2017-06-18 | 2017-09-22 | 西安优耐特容器制造有限公司 | 一种薄壁管件焊接成型控制系统及方法 |
CN107292035A (zh) * | 2017-06-23 | 2017-10-24 | 江苏科技大学 | 含残余应力支架随机振动疲劳寿命的预测方法 |
CN107506544A (zh) * | 2017-08-18 | 2017-12-22 | 中国航发北京航空材料研究院 | 一种金属材料疲劳裂纹扩展预测方法 |
CN109388878A (zh) * | 2018-09-30 | 2019-02-26 | 北京航空航天大学 | 一种综合考虑喷丸强化效果的疲劳寿命预测方法 |
CN110133016A (zh) * | 2019-04-16 | 2019-08-16 | 平高集团有限公司 | 一种焊接数值模拟辅助x射线衍射检测残余应力的方法 |
-
2019
- 2019-11-28 CN CN201911196181.9A patent/CN110954349B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102426068A (zh) * | 2011-09-14 | 2012-04-25 | 华东理工大学 | 一种表面薄膜残余应力的预测方法 |
CN102967397A (zh) * | 2012-11-19 | 2013-03-13 | 辽宁省电力有限公司电力科学研究院 | 焊接结构焊缝超声波残余应力测量装置及方法 |
CN104850691A (zh) * | 2015-05-05 | 2015-08-19 | 南京市特种设备安全监督检验研究院 | 一种基于多因素融合修正的结构件裂纹扩展预测方法 |
CN104848969A (zh) * | 2015-05-22 | 2015-08-19 | 华中科技大学 | 一种基于有限测试点的构件残余应力场的预测方法 |
CN106815419A (zh) * | 2017-01-03 | 2017-06-09 | 东南大学 | 一种基于裂纹信息预测的起重机运行状态在线评估方法 |
CN107186342A (zh) * | 2017-06-18 | 2017-09-22 | 西安优耐特容器制造有限公司 | 一种薄壁管件焊接成型控制系统及方法 |
CN107292035A (zh) * | 2017-06-23 | 2017-10-24 | 江苏科技大学 | 含残余应力支架随机振动疲劳寿命的预测方法 |
CN107506544A (zh) * | 2017-08-18 | 2017-12-22 | 中国航发北京航空材料研究院 | 一种金属材料疲劳裂纹扩展预测方法 |
CN109388878A (zh) * | 2018-09-30 | 2019-02-26 | 北京航空航天大学 | 一种综合考虑喷丸强化效果的疲劳寿命预测方法 |
CN110133016A (zh) * | 2019-04-16 | 2019-08-16 | 平高集团有限公司 | 一种焊接数值模拟辅助x射线衍射检测残余应力的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110954349A (zh) | 2020-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102221503B (zh) | 单轴拉伸全程真应力-真应变曲线测试技术 | |
CN101870075B (zh) | 一种基于性能退化模型的数控机床丝杠副寿命预测方法 | |
CN101870076B (zh) | 一种基于性能退化模型的数控机床导轨副寿命预测方法 | |
CN114619292B (zh) | 一种基于小波降噪和注意力机制融合gru网络的铣削加工刀具磨损监测方法 | |
CN103551709B (zh) | 一种a-tig焊焊缝成型的检测装置及方法 | |
CN1701877A (zh) | 一种基于振动监测的连铸下渣检测方法及装置 | |
CN108709744B (zh) | 一种变负载工况下电机轴承故障诊断方法 | |
CN110954349B (zh) | 一种基于残余应力畸变率的起重机结构健康状态监测方法 | |
CN111460702B (zh) | 一种基于正逆向损伤特征融合的结构件损伤识别方法 | |
CN110887899B (zh) | 一种汽轮机叶片水蚀缺陷监测与识别方法 | |
CN110544034A (zh) | 一种输电线路交叉跨越风险评估方法 | |
CN106202937B (zh) | M50钢锻造组织中碳化物尺寸预测方法 | |
CN109408998B (zh) | 基于样本持续增量快速获取应力谱进行疲劳寿命评估方法 | |
CN102682212B (zh) | 机电类产品的一种可靠性度量方法 | |
CN118195153A (zh) | 一种基于大数据的园区电碳排放预测方法 | |
CN112287586A (zh) | 基于中小桥损伤识别结果的快速可靠度评估方法 | |
CN107679630A (zh) | 一种基于比例维修模型的维修作业时间估计方法 | |
CN115600695A (zh) | 一种计量设备的故障诊断方法 | |
CN111079267B (zh) | 一种基于微观残余应力突度的装载机底盘可靠性检测方法 | |
CN114596296A (zh) | 一种高灵敏度的热轧钢卷端面缺陷识别系统及方法 | |
CN102518012B (zh) | 融合多种设备检测数据进行轨道检测的方法 | |
CN112275985A (zh) | 一种模锻压机行程改进补偿方法 | |
CN110980527B (zh) | 基于顺位竞争系数修正残余应力的起重机健康监测方法 | |
CN112116003B (zh) | 一种复杂传感器网络自适应动态故障诊断方法及系统 | |
Bai et al. | A probabilistic approach for long-term fatigue analysis of onshore wind turbine tower |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |