CN105171778B - 一种工业机器人关节减速器扭振提取与分析方法 - Google Patents
一种工业机器人关节减速器扭振提取与分析方法 Download PDFInfo
- Publication number
- CN105171778B CN105171778B CN201510658081.9A CN201510658081A CN105171778B CN 105171778 B CN105171778 B CN 105171778B CN 201510658081 A CN201510658081 A CN 201510658081A CN 105171778 B CN105171778 B CN 105171778B
- Authority
- CN
- China
- Prior art keywords
- counter
- numerical value
- instantaneous
- reading
- turn
- 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.)
- Expired - Fee Related
Links
Landscapes
- Manipulator (AREA)
Abstract
一种工业机器人关节减速器扭振提取与分析方法,先通过光栅角度编码器获得瞬时转角,然后通过瞬时转角计算瞬时转速,再观测转速是否稳定,不稳定则继续观测,稳定则开始数据采集,读取的计数器的数值和计时器的时间值,通过分析获得瞬时扭转角,然后绘制时间‑瞬时扭转角曲线,再进行整转截断,最后进行时频转换,获得扭振的频域特性,进而得到扭振幅值,本发明能够对低转速、小振幅的扭振进行测试,同时能够基本保证等采样频率采样,降低变采样频率产生的影响。
Description
技术领域
本发明属于工业机器人技术领域,具体涉及一种工业机器人关节减速器扭振提取与分析方法。
背景技术
在各类回转机械中,轴的扭转振动是重要的振动形式之一。轴系的扭转振动的产生原因有很多种,如轴系的负载扭矩的变化或者动力***工况的突然变化以及***零部件的一些损伤和故障等都会激起轴系的扭振。扭振产生的动应力会导致轴的疲劳甚至损坏。扭振对轴及轴上零件的危害在振动的初期表现得并不明显,但是随着扭振的交变应力的长时间作用,不断累积,就可能造成轴的疲劳损坏,从而产生突发的破坏事故。对于汽轮发电机组这类大型回转机械,扭振的危害是非常严重的,其表现之一是机组的次同步振荡;对于机床,其传动链轴系发生严重扭振时,会破坏加工精度;对于发动机,曲轴受到周期性的扭振作用,可能产生疲劳断裂,造成不可估量的破坏。
对于工业机器人关节减速器,扭转振动不仅可以反映出其内部零部件的运行状态,同时也是工业机器人机械臂进行正确且平稳地轮廓控制的重要特性。
传统的扭转振动测试方法有贴应变片法、安装加速度计法、齿测法。贴应变片和安装加速度计属于接触式测量,需要对轴系结构进行一定的改造,如为保证加速度计切向布置而进行的结构调整、安装集流环等,使轴系的结构发生了变化,不仅结构复杂,还对测试结果产生了一定的影响。对于齿测法,被测轴上的等分齿形结构一般齿数较少,轴转动一周获得的脉冲数也较少,不能测得小幅值的扭振,而且在低转速时,这种缺陷尤为明显。
工业机器人关节减速器传动比大(可达200左右),输出转速非常低(15rpm以下),同时其扭转刚度大,扭振幅值低。目前的测试方法不适应工业机器人关节减速器低转频、小扭振幅值的测试需要。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种工业机器人关节减速器扭振提取与分析方法,能够对低转速、小振幅的扭振进行测试,同时能够基本保证等采样频率采样,降低变采样频率产生的影响。
为了达到上述目的,本发明采取的技术方案为:
一种工业机器人关节减速器扭振提取与分析方法,包括以下步骤:
1)、通过光栅角度编码器获得瞬时转角:用计数器对光栅角度编码器的方波信号进行计数,以采样频率f读取计数器中的数值,
则读数周期ΔT为:
第i次读取计数器的数值时获得的瞬时转角为:
式中:Ki——第i次读取计数器中的数值,
p——光栅角度编码器的分辨率,
2)、通过瞬时转角计算瞬时转速:第i次读取计数器的数值时获得的瞬时转速ωi为:
式中:Ki-1——第i-1次读取计数器中的时间值,
ΔT——计数器的读数周期,
f——采样频率,
3)、观测转速是否稳定,不稳定则继续观测,直至稳定则开始数据采集;
4)、数据采集:
4.1)以光栅角度编码器z相脉冲触发计数器和计时器,使计数器置零并对光栅角度编码器的信号开始计数,使计时器置零并开始计时;
4.2)以固定的读数周期ΔT读取计数器中的数值;
4.3)等光栅角度编码器下一个z相脉冲触发计数器和计时器时,读取计数器中的数值并置零,重新开始计数;读取计时器中的时间值并置零重新开始计时,
则有一转内的平均转速为:
式中:K——一转时,计数器中的数值,
t——一转时,计时器中的时间值,
p——光栅角度编码器的分辨率,
5)、根据步骤4)中读取的计数器的数值和计时器的时间值,通过分析获得瞬时扭转角,
在第i-1次读数与第i次读数之间的扭转角:
式中:θi——在第i-1次读数与第i次读数之间的扭转角,
ωi——第i次读取计数器的数值时获得的瞬时转速,
——一转内的平均转速,
ΔT——计数器的读数周期,
f——采样频率,
Ki——第i次读取计数器中的数值,
Ki-1——第i-1次读取计数器中的数值,
则有第i次读数时的瞬时扭转角θ(i)为:
式中:θ(i)——第i-1次读数与第i次读数之间的扭转角,
Ki——第i次读取计数器中的数值,
K——一转时,计数器中的数值,
p——光栅角度编码器的分辨率,
t——转过一转的时间,
i——读数次数,
6)、根据步骤5)中获得的瞬时扭转角,绘制时间-瞬时扭转角曲线;
7)、对步骤6)中获得时间-瞬时扭转角曲线进行整转截断;
8)、对步骤7)获得的整转截断的时间-瞬时扭转角曲线进行时频转换,获得扭振的频域特性,进而得到扭振幅值。
所述的光栅角度编码器每转360000个脉冲,通过四倍频计数,角度分辨率可达0.9角秒。
本发明的有益效果为:
本发明通过光栅角度编码器直接读取计数器和计时器的值,通过分析获得瞬时扭转角曲线,进而得到扭振的扭角幅值,不需要改造轴系结构,能够对低转速、小振幅的扭振进行测试,同时能够基本保证等采样频率采样,降低变采样频率产生的影响。
附图说明
附图为本发明方法的流程示意图。
具体实施方式
下面结合附图对本发明做进一步描述。
参照附图,一种工业机器人关节减速器扭振提取与分析方法,包括以下步骤:
1)、通过光栅角度编码器获得瞬时转角:用计数器对光栅角度编码器的方波信号进行计数,以采样频率f读取计数器中的数值,
则读数周期ΔT为:
第i次读取计数器的数值时获得的瞬时转角为:
式中:Ki——第i次读取计数器中的数值,
p——光栅角度编码器的分辨率,
2)、通过瞬时转角计算瞬时转速:第i次读取计数器的数值时获得的瞬时转速ωi为:
式中:Ki-1——第i-1次读取计数器中的时间值,
ΔT——计数器的读数周期,
f——采样频率,
3)、观测转速是否稳定,不稳定则继续观测,直至稳定则开始数据采集;
4)、数据采集:
4.1)以光栅角度编码器z相脉冲触发计数器和计时器,使计数器置零并对光栅角度编码器的信号开始计数,使计时器置零并开始计时;
4.2)以固定的读数周期ΔT读取计数器中的数值;
4.3)等光栅角度编码器下一个z相脉冲触发计数器和计时器时,读取计数器中的数值并置零,重新开始计数;读取计时器中的时间值并置零重新开始计时,
则有一转内的平均转速为:
式中:K——一转时,计数器中的数值,
t——一转时,计时器中的时间值,
p——光栅角度编码器的分辨率,
5)、根据步骤4)中读取的计数器的数值和计时器的时间值,通过分析获得瞬时扭转角,
在第i-1次读数与第i次读数之间的扭转角:
式中:θi——在第i-1次读数与第i次读数之间的扭转角,
ωi——第i次读取计数器的数值时获得的瞬时转速,
——一转内的平均转速,
ΔT——计数器的读数周期,
f——采样频率,
Ki——第i次读取计数器中的数值,
Ki-1——第i-1次读取计数器中的数值,
则有第i次读数时的瞬时扭转角θ(i)为:
式中:θ(i)——第i-1次读数与第i次读数之间的扭转角,
Ki——第i次读取计数器中的数值,
K——一转时,计数器中的数值,
p——光栅角度编码器的分辨率,
t——转过一转的时间,
i——读数次数,
6)、根据步骤5)中获得的瞬时扭转角,绘制时间-瞬时扭转角曲线;
7)、对步骤6)中获得时间-瞬时扭转角曲线进行整转截断;
8)、对步骤7)获得的整转截断的时间-瞬时扭转角曲线进行时频转换,获得扭振的频域特性,进而得到扭振幅值。
所述的光栅角度编码器每转360000个脉冲,通过四倍频计数,角度分辨率可达0.9角秒。
Claims (2)
1.一种工业机器人关节减速器扭振提取与分析方法,其特征在于,包括以下步骤:
1)、通过光栅角度编码器获得瞬时转角:用计数器对光栅角度编码器的方波信号进行计数,以采样频率f读取计数器中的数值,
则读数周期ΔT为:
第i次读取计数器的数值时获得的瞬时转角为:
式中:Ki——第i次读取计数器中的数值,
p——光栅角度编码器的分辨率,
2)、通过瞬时转角计算瞬时转速:第i次读取计数器的数值时获得的瞬时转速ωi为:
式中:Ki-1——第i-1次读取计数器中的数值,
ΔT——计数器的读数周期,
f——采样频率,
3)、观测转速是否稳定,不稳定则继续观测,直至稳定则开始数据采集;
4)、数据采集:
4.1)以光栅角度编码器z相脉冲触发计数器和计时器,使计数器置零并对光栅角度编码器的信号开始计数,使计时器置零并开始计时;
4.2)以固定的读数周期ΔT读取计数器中的数值;
4.3)等光栅角度编码器下一个z相脉冲触发计数器和计时器时,读取计数器中的数值并置零,重新开始计数;读取计时器中的时间值并置零重新开始计时,
则有一转内的平均转速为:
式中:K——一转时,计数器中的数值,
t——一转时,计时器中的时间值,
p——光栅角度编码器的分辨率,
5)、根据步骤4)中读取的计数器的数值和计时器的时间值,通过分析获得瞬时扭转角,
在第i-1次读数与第i次读数之间的扭转角:
式中:θi——在第i-1次读数与第i次读数之间的扭转角,
ωi——第i次读取计数器的数值时获得的瞬时转速,
——一转内的平均转速,
ΔT——计数器的读数周期,
f——采样频率,
Ki——第i次读取计数器中的数值,
Ki-1——第i-1次读取计数器中的数值,
则有第i次读数时的瞬时扭转角θ(i)为:
式中:θ(i)——第i-1次读数与第i次读数之间的扭转角,
Ki——第i次读取计数器中的数值,
K——一转时,计数器中的数值,
p——光栅角度编码器的分辨率,
t——转过一转的时间,
i——读数次数,
6)、根据步骤5)中获得的瞬时扭转角,绘制时间-瞬时扭转角曲线;
7)、对步骤6)中获得时间-瞬时扭转角曲线进行整转截断;
8)、对步骤7)获得的整转截断的时间-瞬时扭转角曲线进行时频转换,获得扭振的频域特性,进而得到扭振幅值。
2.根据权利要求1所述的一种工业机器人关节减速器扭振提取与分析方法,其特征在于:所述的光栅角度编码器每转360000个脉冲,通过四倍频计数,角度分辨率可达0.9角秒。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510658081.9A CN105171778B (zh) | 2015-10-13 | 2015-10-13 | 一种工业机器人关节减速器扭振提取与分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510658081.9A CN105171778B (zh) | 2015-10-13 | 2015-10-13 | 一种工业机器人关节减速器扭振提取与分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105171778A CN105171778A (zh) | 2015-12-23 |
CN105171778B true CN105171778B (zh) | 2017-01-18 |
Family
ID=54894445
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510658081.9A Expired - Fee Related CN105171778B (zh) | 2015-10-13 | 2015-10-13 | 一种工业机器人关节减速器扭振提取与分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105171778B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108507664A (zh) * | 2017-02-27 | 2018-09-07 | 华晨汽车集团控股有限公司 | 一种汽车传动轴扭振的快速判定方法 |
CN113532891A (zh) * | 2021-06-16 | 2021-10-22 | 武汉理工大学 | 轮轴两端的转速、扭振和扭矩信号同步测量方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1453929A (zh) * | 2002-04-18 | 2003-11-05 | 佳能株式会社 | 振动波驱动装置、振动体以及振动波驱动装置的驱动*** |
CN1740751A (zh) * | 2004-08-23 | 2006-03-01 | 索尼株式会社 | 角度检测信号处理装置 |
WO2014039846A2 (en) * | 2012-09-07 | 2014-03-13 | Dana Limited | Active torque ripple rectifier |
CN104034407A (zh) * | 2014-06-10 | 2014-09-10 | 国家电网公司 | 减小旋转机械扭振信号脉冲测量法中周期性误差的方法 |
CN102183290B (zh) * | 2011-01-24 | 2015-09-30 | 哈尔滨工程大学 | 一种提高轴系扭振信号时频转换精度的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002139374A (ja) * | 2000-10-30 | 2002-05-17 | Mitsubishi Heavy Ind Ltd | 回転捩り振動計測方法及び装置 |
-
2015
- 2015-10-13 CN CN201510658081.9A patent/CN105171778B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1453929A (zh) * | 2002-04-18 | 2003-11-05 | 佳能株式会社 | 振动波驱动装置、振动体以及振动波驱动装置的驱动*** |
CN1740751A (zh) * | 2004-08-23 | 2006-03-01 | 索尼株式会社 | 角度检测信号处理装置 |
CN102183290B (zh) * | 2011-01-24 | 2015-09-30 | 哈尔滨工程大学 | 一种提高轴系扭振信号时频转换精度的方法 |
WO2014039846A2 (en) * | 2012-09-07 | 2014-03-13 | Dana Limited | Active torque ripple rectifier |
CN104034407A (zh) * | 2014-06-10 | 2014-09-10 | 国家电网公司 | 减小旋转机械扭振信号脉冲测量法中周期性误差的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105171778A (zh) | 2015-12-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Reddy et al. | Detection and monitoring of coupling misalignment in rotors using torque measurements | |
Yoon et al. | On the use of a single piezoelectric strain sensor for wind turbine planetary gearbox fault diagnosis | |
EP2259033B1 (en) | Angular and torsional vibration monitoring of rotordynamic systems | |
CN102507205B (zh) | 一种检测航空发动机风扇叶片颤振故障的方法 | |
CN102066695B (zh) | 用于确定涡轮喷气发动机转子角度位置的方法和*** | |
CN105675113B (zh) | 基于微型传感器的旋转机械角域振动信号获取装置和方法 | |
CN101487756A (zh) | 旋转机械振动分析中的谐分量转速平衡方法 | |
CN104819841B (zh) | 基于内置编码信息的单传感柔性角度域平均方法 | |
Yu et al. | A data processing method for determining instantaneous angular speed and acceleration of crankshaft in an aircraft engine–propeller system using a magnetic encoder | |
CN103363921A (zh) | 一种改进型三点法回转误差、圆度误差计算方法 | |
CN104006966A (zh) | 载荷对齿轮轴系振动和齿轮箱体振动相互耦合影响规律的实验装置及验证方法 | |
CN110319922A (zh) | 基于叶尖定时原理的叶片扭转振动位移测量方法 | |
CN105171778B (zh) | 一种工业机器人关节减速器扭振提取与分析方法 | |
CN105628176A (zh) | 一种旋转机械扭转振动信号采集分析方法 | |
André et al. | Comparison between angular sampling and angular resampling methods applied on the vibration monitoring of a gear meshing in non stationary conditions | |
Wang et al. | A two-stage method using spline-kernelled chirplet transform and angle synchronous averaging to detect faults at variable speed | |
Liang et al. | Application of instantaneous rotational speed to detect gearbox faults based on double encoders | |
CN109540482A (zh) | 一种涡轮机叶片无键相高倍频同步振动参数分析方法及分析装置 | |
CN106441840A (zh) | 一种基于模型的旋转机械伪次同步故障全息诊断方法 | |
Bulatović et al. | Measurement and analysis of angular velocity variations of twelve-cylinder diesel engine crankshaft | |
JPH04270935A (ja) | 機械の回転部分の運動状態をモニターする方法及び装置 | |
CN203349773U (zh) | 精密离心机主轴径向回转误差与倾角回转误差测试装置 | |
CN106092567A (zh) | 旋转机械故障振动同步测量方法 | |
CN117686232A (zh) | 一种燃气轮机振动基频实时提取方法、装置及存储介质 | |
CN110219816A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170118 Termination date: 20191013 |