CN105241589A - 一种机器人手臂应变测试数据处理方法 - Google Patents
一种机器人手臂应变测试数据处理方法 Download PDFInfo
- Publication number
- CN105241589A CN105241589A CN201510555314.2A CN201510555314A CN105241589A CN 105241589 A CN105241589 A CN 105241589A CN 201510555314 A CN201510555314 A CN 201510555314A CN 105241589 A CN105241589 A CN 105241589A
- Authority
- CN
- China
- Prior art keywords
- stress
- value
- curve
- damage
- robot arm
- 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.)
- Pending
Links
Abstract
本发明公开了一种机器人手臂应变测试数据处理方法,它由应变片(1),应变采集器(2),含有峰谷预处理器(4)和改进的雨流计数器(5)的数据处理系统(3),损伤与寿命计算器(6),报警器(8)和显示器(7)并依次连接组成的测试装置,其数据处理方法为:(1)贴应变片;(2)峰谷预处理;(3)改进的雨流计数统计;(4)损伤计算;(5)损伤累积计算。该方法中监测装置结构简单,体积小,数据处理方法先进,精度高,能自动显示数据处理结果,操作简单,直接安装在机器人手臂上,实时自动监控和自动报警。
Description
技术领域
本发明涉及机器人手臂的测试和数据处理,具体涉及一种机器人手臂的应变测试数据处理方法。
背景技术
伺服电机驱动机器人手臂进行水平方向和竖直方向的运动,实现对目标物的抓放作业。机器人手臂在操作的过程中,受变幅交变应力作用而产生疲劳损伤,随着工作时间的增加,机器人手臂(机器人手臂零部件)的疲劳损伤也在增加,其使用寿命逐渐减少。机器人手臂的疲劳损伤程度和剩余寿命,是由实验计算和使用经验积累获得并在此基础上编制维护资料,它规定了维护的时间间隔并提供了机器人手臂的寿命值。虽然这些数据对保证机器人的连续运转起到了重要的作用,但是该资料无通用性。对于不同用途、不同型号、不同工作载荷的机器人,随着运作时间的增加,机器人手臂的疲劳损伤程度也不同,仅靠维护资料来预估每次运作后的疲劳损伤程度及剩余寿命,不能及时准确地提供机器人手臂的疲劳损伤程度及剩余寿命。由于难以预测机器人手臂寿命值,也就难以预防某一台机器人在运作过程中突然失效,进而影响机器人运作的可靠稳定性。为了提高机器人运作的可靠稳定性,对各机器人手臂的损伤程度及剩余寿命进行实时监控尤为重要。
发明内容
针对已有技术存在的不足,本发明的目的是提供一种结构简单,体积小,安装方便,操作简单,能实时自动监测机器人手臂损伤和寿命的机器人手臂应变测试数据处理方法。
本发明是通过如下技术方案实现以上目的的。
本发明的一种机器人手臂应变测试数据处理方法,它由应变片1,应变采集器2,数据处理系统3、装在数据处理系统3内的峰谷预处理器4和改进的雨流计数器5,损伤与寿命计算器6,显示器7和报警器8并依次连接组成的测试装置,其数据处理方法的步骤如下:
(1)、贴应变片
在被测试的机器人手臂的应力(或疲劳损伤)集中部位上贴上应变片1,应变采集器2通过应变片1以一定的时间间隔连续获得贴片部位的应变大小,并将其转换成应力数值信号,生成机器人手臂各应力数值信号的数据文件,输入到含有峰谷预处理器4和改进的雨流计数器5的数据处理系统3;
(2)、峰谷预处理
以一定的时间间隔采集机器人手臂各应力数值信号的数据,将采集到的随时间变化的应力数据点连成应力-时间历程曲线,去掉相邻峰谷点间的中间点,保留峰谷点,得到包含峰谷点的随时间变化的应力-时间历程曲线;
(3)、采用改进的雨流计数器5对封闭历程曲线进行应力循环统计
在应力-时间坐标系中,将保留峰谷点的应力-时间历程曲线以过曲线终止点的竖直线为对称轴,镜像原应力-时间历程曲线,得到对称应力-时间历程曲线,并以对称应力-时间历程曲线最小应力值点为分割点,将对称应力-时间历程曲线分为左、右两段应力-时间历程曲线,将左、右两段应力-时间历程曲线的位置互换,得到起始点应力值与终止点应力值,其应力值相等且均为最小应力值的封闭历程曲线,
采用改进的雨流计数器5对封闭历程曲线进行应力循环统计,得到每一应力循环的应力均值σm、应力幅值σa和循环次数n,根据得到的应力均值σm、应力幅值σa和循环次数n计算出等效对称循环应力σ和有效应力循环次数n',其计算式如下:
式中:σi为第i次等效对称循环应力值,σa'为计算应力幅值,σm'为计算应力均值,σb'为计算强度极限值,[n]为设计安全系数;
(4)、损伤计算
根据材料S-N曲线表达式,采用损伤与寿命计算器6,计算第i次的损伤值:
式中:a和b为待定系数,Ni为第i次的疲劳寿命值,Di为第i次的损伤值;
(5)、损伤累积计算
采用损伤与寿命计算器6计算损伤累积值,其计算式为:
D=∑Di
式中:D为损伤累积值,
当D=1时,剩余寿命减少到零,更换或修复机器人手臂。
本发明的一种机器人手臂应变测试数据处理方法的优点是:该方法中监测装置结构简单,体积小,数据处理方法先进,精度高,能自动显示数据处理结果,操作简单,直接安装在机器人手臂上,实时自动监控和自动报警。
附图说明
图1是本发明的一种机器人手臂应变测试数据处理方法的流程图。
图2是一段应力随时间变化的应力-时间历程曲线图。
图3是峰谷预处理后的应力-时间历程曲线图。
图4是镜像处理后的应力-时间历程曲线图。
图5是分割处理后的封闭时间-历程历程曲线图。
图6是分解封闭时间-历程历程曲线图成单峰循环区间示意图。
图中,1、应变片,2、应变采集器,3、数据处理系统,4、峰谷预处理器,5、改进的雨流计数器,6、损伤与寿命计算器,7、显示器,8、报警器。
具体实施方式
由图1可见,本发明的一种机器人手臂应变测试数据处理方法,由应变片1,应变采集器2,数据处理系统3,装在数据处理系统3内的峰谷预处理器4和改进的雨流计数器5,损伤与寿命计算器6,显示器7,报警器8并依次连接组成测试装置,其所述方法的步骤如下:
(1)、贴应变片
应变片1贴在被检机器人手臂的应力(或疲劳损伤)集中部位上,应变采集器2通过应变片1以一定的时间间隔(如每过0.01秒将测到的值保存起来,10秒钟内将有1000个按时间顺序排列的测试的数据)连续获得贴片部位的应变大小,并将其转换成应力数值信号,生成机器人手臂各应力数值信号的数据文件,输入到含有峰谷预处理器4和改进的雨流计数器5的数据处理系统3;
(2)、峰谷预处理
将以一定的时间间隔采集的应力数据点(如10秒内测得并转换而得的1000个随时间变化的应力数据点)连成应力-时间历程曲线,采用峰谷预处理器4将随时间变化的应力-时间历程曲线图(如图2所示)中去掉相邻峰谷点间的中间点,保留峰谷点,得到包含峰谷点的随时间变化的应力-时间历程曲线图(如图3所示);
(3)、采用改进的雨流计数器5对封闭历程曲线进行应力循环统计
在应力-时间坐标系中,将保留峰谷点的应力-时间历程曲线以过曲线终止点的竖直线为对称轴,镜像原应力-时间历程曲线,得到对称的应力-时间历程曲线图(如图4所示),并以对称的应力-时间历程曲线最小应力值点为分割点,将对称应力-时间历程曲线分割为左、右两段应力-时间历程曲线,将左、右两段应力-时间历程曲线互换位置,得到起始点应力值与终止点应力值,其应力值相等且均为最小应力值的封闭历程曲线(如图5所示)。
采用改进的雨流计数器5对封闭历程曲线进行应力循环统计,得到每一应力循环的应力均值σm、应力幅值σa和循环次数n,根据得到的应力均值σm、应力幅值σa和循环次数n计算出等效对称循环应力σ和有效应力循环次数n',其计算式如下:
式中:σi为第i次等效对称循环应力值,σa'为计算应力幅值,σm'为计算应力均值,σb'为计算强度极限值,[n]为设计安全系数;
设经统计得第i次应力循环的应力均值σm为80Mpa、应力幅值σa为40Mpa和循环次数n为1,令材料计算强度极限值σb'为720Mpa,安全系数[n]为2.4,则经过上式计算可得第i次等效对称循环应力值σi为130.91Mpa和有效应力循环次数n'为0.5;
(4)、损伤计算
根据材料S-N曲线双坐标表达式,采用损伤与寿命计算器6,计算第i次的损伤值:
式中,a和b为待定系数,Ni为第i次的疲劳寿命值,Di为第i次的损伤值;
(5)、损伤累积计算
采用损伤与寿命计算器6计算损伤累积值,其计算式为:
D=∑Di
式中:D为损伤累积值
当D大于等于0.99时,报警器8自动报警,并在显示器7中实时监控,提示机器人手臂损伤累积值,当机器人手臂的损伤累积值达到寿命极限时,及时对机器人手臂进行维护检查或更换,预防机器人的突然失效,提高生产效率。
所述改进的雨流计数器5内置有计算机程序语言,将上述封闭应力-时间历程曲线分区,分成为一个个单峰循环区间局部应力-时间历程曲线(如图6所示),设置全局静态数组S[n](Nstart≤n≤Nend)是经过处理的从起始点Nstart到终止点Nend的应力-时间历程记录,且S[Nstart]=S[Nend]=Smin,该改进的雨流计数器5用计算机程序语言描述如下:
VoidRain_flow(intNstart,intNend,doubleSmin)
{
(1)如果(Nend-Nstart<2),则(return;)
(2)如果(Nend-Nstart=2),则
(①σmax=S[Nstart+1],σmin=Smin;
②σm=(σmax+σmin)/2,σa=(σmax-σmin)/2;
③输出σm和σa;
④return;)
(3)如果(Nend-Nstart>2),则
(①在数组S[n],Nstart﹤n﹤Nend中寻找最大值Smax,记下Npeak及
Smax=S[Npeak];
②如果Nstart+1=Npeak,则令NLV=Nstart,否则,在数组S[n],Nstart
﹤n﹤Npeak中寻找最小值SLmin,记下NLV及SLmin=S[NLV];
③如果Npeak+1=Nend,则令NRV=Nend,否则,在数组S[n],Npeak
﹤n﹤Nend中寻找最小值SRmin,记下NRV及SRmin=S[NRV];
④Rain_flow(Nstart,NLV,SLmin);
⑤Rain_flow(NLV,NRV,Smin);
⑥Rain_flow(NRV,Nend,SRmin);)
}。
Claims (1)
1.一种机器人手臂应变测试数据处理方法,其特征在于,它由应变片(1),应变采集器(2),数据处理系统(3)、装在数据处理系统(3)内的峰谷预处理器(4)和改进的雨流计数器(5),损伤与寿命计算器(6),显示器(7)和报警器(8)并依次连接组成的测试装置,其数据处理方法的步骤如下:
(1)、贴应变片
在被测试的机器人手臂的应力(或疲劳损伤)集中部位上贴上应变片(1),应变采集器(2)通过应变片(1)以一定的时间间隔连续获得贴片部位的应变大小,并将其转换成应力数值信号,生成机器人手臂各应力数值信号的数据文件,输入到含有峰谷预处理器(4)和改进的雨流计数器(5)的数据处理系统(3);
(2)、峰谷预处理
以一定的时间间隔采集机器人手臂各应力数值信号的数据,将采集到的随时间变化的应力数据点连成应力-时间历程曲线,去掉相邻峰谷点间的中间点,保留峰谷点,得到包含峰谷点的随时间变化的应力-时间历程曲线;
(3)、采用改进的雨流计数器(5)对封闭历程曲线进行应力循环统计
在应力-时间坐标系中,将保留峰谷点的应力-时间历程曲线以过曲线终止点的竖直线为对称轴,镜像原应力-时间历程曲线,得到对称应力-时间历程曲线,并以对称应力-时间历程曲线最小应力值点为分割点,将对称应力-时间历程曲线分为左、右两段应力-时间历程曲线,将左、右两段应力-时间历程曲线的位置互换,得到起始点应力值与终止点应力值,其应力值相等且均为最小应力值的封闭历程曲线,
采用改进的雨流计数器(5)对封闭历程曲线进行应力循环统计,得到每一应力循环的应力均值σm、应力幅值σa和循环次数n,根据得到的应力均值σm、应力幅值σa和循环次数n计算出等效对称循环应力σ和有效应力循环次数n',其计算式如下:
式中:σi为第i次等效对称循环应力值,σa'为计算应力幅值,σm'为计算应力均值,σb'为计算强度极限值,[n]为设计安全系数;
(4)、损伤计算
根据材料S-N曲线表达式,采用损伤与寿命计算器(6),计算第i次的损伤值:
式中:a和b为待定系数,Ni为第i次的疲劳寿命值,Di为第i次的损伤值;
(5)、损伤累积计算
采用损伤与寿命计算器(6)计算损伤累积值,其计算式为:
D=∑Di
式中:D为损伤累积值,
当D=1时,剩余寿命减少到零,更换或修复机器人手臂。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510555314.2A CN105241589A (zh) | 2015-09-02 | 2015-09-02 | 一种机器人手臂应变测试数据处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510555314.2A CN105241589A (zh) | 2015-09-02 | 2015-09-02 | 一种机器人手臂应变测试数据处理方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105241589A true CN105241589A (zh) | 2016-01-13 |
Family
ID=55039313
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510555314.2A Pending CN105241589A (zh) | 2015-09-02 | 2015-09-02 | 一种机器人手臂应变测试数据处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105241589A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106043739A (zh) * | 2016-07-30 | 2016-10-26 | 董超超 | 一种具有智能监测功能的飞机 |
CN106290559A (zh) * | 2016-07-30 | 2017-01-04 | 董超超 | 发动机剩余寿命预测系统 |
CN107220713A (zh) * | 2017-06-06 | 2017-09-29 | 上海理工大学 | 基于健康状态的机器人手臂实时保养方法 |
CN108709672A (zh) * | 2018-05-31 | 2018-10-26 | 哈尔滨理工大学 | 一种手臂压力信号采集方法 |
CN110926785A (zh) * | 2019-09-02 | 2020-03-27 | 北京源清慧虹信息科技有限公司 | 一种监测钢结构疲劳损伤的方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101158570A (zh) * | 2007-11-20 | 2008-04-09 | 中国航空动力机械研究所 | 减速器应变测试装置及数据处理方法 |
CN101498575A (zh) * | 2007-12-11 | 2009-08-05 | 汉密尔顿毕克布兰德斯有限公司 | 用于测量转动叶片上的应力的装置及方法 |
CN102767133A (zh) * | 2011-05-06 | 2012-11-07 | 东南大学 | 一种钢箱梁桥疲劳应力监测的传感器布置方法 |
CN104101500A (zh) * | 2014-07-25 | 2014-10-15 | 北车风电有限公司 | 一种监测风电机组疲劳状态的方法 |
-
2015
- 2015-09-02 CN CN201510555314.2A patent/CN105241589A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101158570A (zh) * | 2007-11-20 | 2008-04-09 | 中国航空动力机械研究所 | 减速器应变测试装置及数据处理方法 |
CN101498575A (zh) * | 2007-12-11 | 2009-08-05 | 汉密尔顿毕克布兰德斯有限公司 | 用于测量转动叶片上的应力的装置及方法 |
CN102767133A (zh) * | 2011-05-06 | 2012-11-07 | 东南大学 | 一种钢箱梁桥疲劳应力监测的传感器布置方法 |
CN104101500A (zh) * | 2014-07-25 | 2014-10-15 | 北车风电有限公司 | 一种监测风电机组疲劳状态的方法 |
Non-Patent Citations (1)
Title |
---|
毕继红等: "基于简化载荷谱法与雨流计数法的接触网疲劳寿命对比分析", 《铁道标准设计》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106043739A (zh) * | 2016-07-30 | 2016-10-26 | 董超超 | 一种具有智能监测功能的飞机 |
CN106290559A (zh) * | 2016-07-30 | 2017-01-04 | 董超超 | 发动机剩余寿命预测系统 |
CN107220713A (zh) * | 2017-06-06 | 2017-09-29 | 上海理工大学 | 基于健康状态的机器人手臂实时保养方法 |
CN107220713B (zh) * | 2017-06-06 | 2020-10-09 | 上海理工大学 | 基于健康状态的机器人手臂实时保养方法 |
CN108709672A (zh) * | 2018-05-31 | 2018-10-26 | 哈尔滨理工大学 | 一种手臂压力信号采集方法 |
CN110926785A (zh) * | 2019-09-02 | 2020-03-27 | 北京源清慧虹信息科技有限公司 | 一种监测钢结构疲劳损伤的方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105241589A (zh) | 一种机器人手臂应变测试数据处理方法 | |
CN107144430B (zh) | 一种基于增量学习的轴承故障诊断方法 | |
CN103971171B (zh) | 一种输电设备状态评估方法 | |
CN112629647B (zh) | 大跨悬索桥涡振事件的实时识别和监测预警方法 | |
CN105300692B (zh) | 一种基于扩展卡尔曼滤波算法的轴承故障诊断及预测方法 | |
JP6557110B2 (ja) | 状態診断装置及びプログラム | |
CN107014668A (zh) | 一种基于压电和智能涂层传感器的疲劳裂纹综合监测方法 | |
CN108376184A (zh) | 一种桥梁健康监控的方法及系统 | |
CN109858104A (zh) | 一种滚动轴承健康评估与故障诊断方法及监测系统 | |
CN106289388A (zh) | 一种船闸健康状态的远程监测系统及其监测方法 | |
CN105716814A (zh) | 一种评估桁架结构损伤的实时监测系统及其方法 | |
CN109236587B (zh) | 一种用于检测风力发电机异常工作的报警系统 | |
CN110990978B (zh) | 一种螺栓状态监测方法和装置 | |
CN105045983A (zh) | 一种基于轴温数据的高速列车车轴老化分析方法 | |
CN116739384A (zh) | 基于5g无线通讯的矿用设备运行管理系统 | |
CN112000654B (zh) | 一种高速铁路大风监测数据质量控制方法和装置 | |
CN110398384A (zh) | 旋转机械的机器健康监测 | |
CN107543581A (zh) | 多功能变电构架健康监测和损伤识别方法 | |
CN103487250B (zh) | 基于二维投射的煤矿设备预知维护方法 | |
CN113177646A (zh) | 一种基于自适应边缘代理的配电设备在线监测方法及系统 | |
CN114738205B (zh) | 一种浮式风机基础的状态监测方法、装置、设备和介质 | |
CN116893643A (zh) | 一种基于数据分析的智能机器人行驶轨迹安全管控系统 | |
CN110987396A (zh) | 一种用于采煤机摇臂的智能故障诊断及寿命预测方法 | |
CN113343541B (zh) | 长大跨桥梁的涡激振动预警方法、装置及终端 | |
CN116877358B (zh) | 一种风力发电机组振动监测分析系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160113 |
|
WD01 | Invention patent application deemed withdrawn after publication |