CN111262496A - 一种伺服控制系统谐振在线抑制的方法及伺服控制系统 - Google Patents

一种伺服控制系统谐振在线抑制的方法及伺服控制系统 Download PDF

Info

Publication number
CN111262496A
CN111262496A CN202010361439.2A CN202010361439A CN111262496A CN 111262496 A CN111262496 A CN 111262496A CN 202010361439 A CN202010361439 A CN 202010361439A CN 111262496 A CN111262496 A CN 111262496A
Authority
CN
China
Prior art keywords
trap
frequency
control system
servo control
notch
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
Application number
CN202010361439.2A
Other languages
English (en)
Inventor
陈小波
张东卫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Dafeng Cnc Technology Co Ltd
Original Assignee
Nanjing Dafeng Cnc Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing Dafeng Cnc Technology Co Ltd filed Critical Nanjing Dafeng Cnc Technology Co Ltd
Priority to CN202010361439.2A priority Critical patent/CN111262496A/zh
Publication of CN111262496A publication Critical patent/CN111262496A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P23/00Arrangements or methods for the control of AC motors characterised by a control method other than vector control
    • H02P23/04Arrangements or methods for the control of AC motors characterised by a control method other than vector control specially adapted for damping motor oscillations, e.g. for reducing hunting
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P23/00Arrangements or methods for the control of AC motors characterised by a control method other than vector control
    • H02P23/14Estimation or adaptation of motor parameters, e.g. rotor time constant, flux, speed, current or voltage

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Electric Motors In General (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种伺服控制系统谐振在线抑制的方法及伺服控制系统,属于控制技术领域。本发明首先设置伺服控制系统响应幅值的阈值;然后采集速度环调节器输出的电流指令,对采集的电流指令进行频谱分析,并根据分析结果,确定陷波器的陷波参数;然后根据所述陷波器的陷波参数确定陷波器的具体传递函数;最后将该陷波器的传递函数应用于所述伺服控制系统的控制环路。本发明通过直接对速度环调节器输出的电流指令进行FFT分析,确定陷波器的陷波参数,从而实现了陷波器参数的在线调节,并且得到的陷波器的传递函数更加符合系统特性,进而能够对伺服控制进行更好的陷波,使得系统具有更好的稳定性。

Description

一种伺服控制系统谐振在线抑制的方法及伺服控制系统
技术领域
本发明属于控制技术领域,具体涉及一种伺服控制系统谐振在线抑制的方法及伺服控制系统。
背景技术
谐振频率是机械结构的固有特性,对伺服控制系统的控制带宽、动态跟踪性能都有重大影响。在不对谐振采取抑制措施的情况下,为了保证系统工作的稳定性,系统的工作带宽就会被抑制在系统的最小谐振频率点以下,这就大大限制了系统的工作带宽,进而影响了系统的动态性能。为了解决机械系统中的谐振问题,科研人员发展了很多抑制谐振的方法,例如常用的方法是在速度环串入低通滤波器来实现机械谐振的抑制。这种结构设计起来比较简单,只需要选择合适的阻尼比和截止频率即可。其中阻尼比的选择是根据幅值裕度和相位裕度来进行合理的选择;而截止频率是根据谐振频率的大小来选择。但是抑振效果不太明显,还会引起相位滞后,不能满足现代伺服控制系统高精的要求。另一种常用方法是基于陷波器的谐振抑制,通过在伺服控制系统的速度环与电流环之间串入一个陷波器,从而降低共振频率处的增益,达到抑制机械共振的效果。陷波器根据参数的不同,分为二参数陷波器和三参数陷波器,其中,最常用的是二参数陷波器,二参数陷波器传递函数的形式为
Figure 308493DEST_PATH_IMAGE001
ω n 表示陷波中心频率,ζ表示陷波器的阻尼比,然而,采用二参数陷波器时,陷波器的参数虽然数量少,但是由于其中心频率点的陷波大小变化较大,使得伺服控制系统在性能上受到极大的影响,当伺服控制系统工作时,调节起来会比较麻烦,因此,现有技术中提出了很多方法,来解决传统的二参数陷波器进行伺服控制系统谐振抑制时,难以满足需求的问题。
申请号为201910327030.6的中国专利中公开了一种用于伺服控制系统谐振抑制、陷波参数优化方法,该专利中提出采用三参数的陷波器进行伺服控制系统的谐振抑制,并且通过FFT分析,在线对陷波器的参数进行优化,通过FFT的分析结果,确定伺服控制系统工作时陷波器的参数。然而,通过对该专利的方法进行分析,发现该专利在采用FFT进行陷波器的参数优化时,先对转速差集进行FFT分析,得到谐振频率,然后根据谐振频率,确定陷波器的参数。但在实际使用时,发现系统仍然存在相对较大的谐振。经过分析,主要原因是现有技术中对陷波器的参数确定的准确性相对较差,从而使得陷波器不能很好地进行陷波,影响伺服控制系统的稳定性。
发明内容
技术问题:本发明提供一种伺服控制系统谐振在线抑制的方法及伺服控制系统,该抑制方法通过直接对速度环调节器的输出作频谱分析,在线调整陷波器的参数,使得陷波器具有更好的陷波性能,并且将该方法用于伺服控制系统,使得控制系统能够更好地抑制谐振,使得伺服控制系统具有更好的稳定性。
技术方案:本发明的伺服控制系统谐振在线抑制的方法,其特征在于,包括:
设置伺服控制系统响应幅值的阈值;
采集速度环调节器输出的电流指令,对采集的电流指令进行频谱分析,并根据分析结果确定陷波器的陷波参数;
根据所述陷波器的陷波参数确定陷波器的具体传递函数;
将该传递函数应用于所述伺服控制系统的控制环路。
进一步地,采用FFT进行频谱分析。
进一步地,所述陷波器的陷波参数包括陷波中心频率、陷波深度和陷波宽度,传递函数的形式为:
Figure 61685DEST_PATH_IMAGE002
式中,ω n 表示陷波器的陷波中心频率,ξ表示陷波深度,k表示陷波宽度,s表示复变量。
进一步地,所述陷波中心频率的确定方法为:采集FFT分析曲线幅值最高点的处的频率作为陷波器的陷波中心频率。
进一步地,所述陷波深度的确定方法为:
Figure DEST_PATH_IMAGE003
式中,H th 表示伺服控制系统响应幅值的阈值,H max表示FFT分析曲线幅值最高点的处幅值。
进一步地,所述陷波宽度的确定方法为:提取FFT分析曲线幅值最高点两侧幅值均为所述阈值的两个频率点,所述两个频率点对应的频率分别为ω 1ω 2,并根据如下公式确定陷波宽度:
Figure 897048DEST_PATH_IMAGE004
式中,ω n 表示陷波中心频率,
Figure 518523DEST_PATH_IMAGE005
表示中心频率ω n 与频率ω 1间的频带宽度,
Figure 775804DEST_PATH_IMAGE006
表示中心频率ω n 与频率ω 2间的频带宽度,m为正数,表示将
Figure 168739DEST_PATH_IMAGE005
Figure 396458DEST_PATH_IMAGE006
中较大值放大m倍,m的范围为1~20。
进一步地,所述将
Figure 306777DEST_PATH_IMAGE007
Figure 342866DEST_PATH_IMAGE006
中较大值放大m倍时,m的范围为1~2。
进一步地,确定所述陷波器的陷波参数的步骤为:
采集速度调节环输出的电流信号,并统计采集到的数据量,如果数据量超过设定的最少数据量,则执行下一步骤,否则,继续采集;
对采集的数据进行FFT分析,并生成FFT分析曲线,根据FFT分析曲线得到陷波中心频率及相应的幅值;
判断陷波中心频率对应的幅值是否大于伺服控制系统响应幅值的阈值,如果大于所述阈值,则根据FFT分析曲线,计算陷波宽度和陷波深度,并开启陷波器,否则,则不计算陷波宽度和陷波深度。
本发明的伺服控制系统,包括:
伺服电机;
编码器,用于采集伺服电机的实际转速,并将采集到的伺服电机的实际转速进行反馈;
速度环调节器,用于将实际输入的转速与所述编码器反馈回的实际转速做差,并输出电流指令;
自适应陷波器,用于对所述速度环调节器输出的电流指令进行陷波,滤除谐振分量,并输出陷波后的电流指令;
FFT模块,用于采集所述速度环调节器输出的所述电流指令,进行FFT分析,并根据FFT分析曲线,得到陷波器的陷波参数,并将所述陷波参数传递给所述陷波器;
电流环调节器,将陷波器输出的电流信号转化为电压信号,所述电压信号驱动所述伺服电机。
进一步地,所述FFT模块包括:
数据采集器,用于采集速度环调节器输出的电流信号;
数据记录器,用于统计所述数据采集器采集的数据量,并判断是否满足需要的最少数据量;
数据处理器,采用FFT作频谱分析,得到FFT分析曲线;
数据提取器,提取FFT分析曲线中最高点处的频率作为陷波中心频率ω n ,以及对应的幅值,提取最高点两侧幅值均为伺服控制系统响应幅值的阈值的两个频率点对应的频率ω 1ω 2
第一计算器:根据公式
Figure 31336DEST_PATH_IMAGE008
计算陷波宽度k,其中m为正数,表示将
Figure DEST_PATH_IMAGE009
Figure 700346DEST_PATH_IMAGE010
中较大值放大m倍,
Figure 398043DEST_PATH_IMAGE011
表示陷波中心频率ω n 与频率ω 1间的频带宽度,
Figure 757481DEST_PATH_IMAGE010
表示陷波中心频率ω n 与频率ω 2间的频带宽度,其中m为正数,范围为1~20;;
第二计算器:根据公式
Figure 164322DEST_PATH_IMAGE012
计算陷波深度ξ,其中,H th 表示伺服控制系统响应幅值的阈值,H max表示FFT分析曲线幅值最高点的处幅值;
第三计算器:根据所述的中心频率、陷波深度和陷波宽度确定传递函数,传递函数的形式为:
Figure 569896DEST_PATH_IMAGE013
,其中s为复变量。
有益效果:本发明与现有技术相比,具有以下优点:
(1)在传统的伺服系统谐波抑制方法中,大多对速度差进行频谱分析,由于准确性较差导致其并不能很好的抑制伺服系统的谐振。本发明摒弃了传统技术中直接利用速度差进行频谱分析的方法,而是利用电流信号进行频谱分析获得陷波参数;并且是在速度环调节器的输出端进行信号采集,对该输出端的电流信号进行频谱分析,利用分析结果确定陷波器的陷波参数,更符合陷波器的陷波特性,更能有效地对伺服控制系统的谐振进行抑制。
(2)本发明的陷波器的传递函数,引入陷波中心频率、陷波深度及陷波宽度乘积作为传递函数分子中一次项的系数,把陷波中心频率与陷波宽度乘积作为传递函数分母一次项的系数,该传递函数所确定的陷波器能够中和谐振引起的幅值,从而达到更好的谐振抑制效果。
(3)本发明对传统的陷波器进行了改进,并采集速度环调节器输出端的电流信号进行FFT分析,根据FFT分析曲线,得到陷波器的陷波参数,从而实现陷波器的陷波参数的在线调节。通过对陷波器传递函数形式的改进,能够直接利用陷波器的输入端的信号获取陷波参数,该参数更符合陷波器的陷波特性,使陷波器对伺服系统具有更好的谐振抑制能力,同时也使得伺服系统动态性能得以提升。
附图说明
图1为本发明的伺服控制系统结构示意图;
图2为本发明的伺服控制系统谐振在线抑制方法的流程图;
图3为本发明中确定陷波器的陷波参数的流程图;
图4为本发明中对速度环调节器输出的电流信号进行FFT频谱分析的示意图;
图5为本发明中伺服控制系统未开启陷波器时系统波形图;
图6为本发明中伺服控制系统未开启陷波器时FFT频谱图;
图7为本发明中伺服控制系统开启陷波器时系统波形图;
图8为本发明中伺服控制系统开启陷波器时FFT频谱图。
具体实施方式
下面结合实施例和说明书附图对本发明作进一步的说明。
需要说明的是,本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
结合图1所示,本发明的伺服控制系统,包括伺服电机;
编码器,用于采集伺服电机的转速ω,并将采集到的伺服电机的转速反馈;
速度环调节器,用于将实际输入的转速ω ref 与反馈回的伺服电机的转速ω做差,得到转速差集,即ω ref -ω,然后将所述转速差集作为输入,并输出电流指令iq ref
陷波器,用于对所述速度环调节器输出的电流指令iq ref 进行陷波,滤除谐振分量,并输出陷波后的电流指令;
FFT模块,用于采集所述速度环调节器输出的所述电流指令iq ref ,进行FFT(FastFourier Transform,快速傅里叶变换)分析,并得到陷波器的陷波参数,并传递给所述自适应陷波器,从而得到陷波器的具体传递函数,以便将陷波器的传递函数应用于伺服控制系统的控制环路;
电流环调节器,将陷波器输出的电流信号转化为电压信号,所述电压信号驱动所述伺服电机。
陷波器设置在速度环调节器和电流环调节器之间,并通过FFT模块在线确定陷波器的参数,确定出更适合于伺服控制系统的陷波器的具体传递函数,从而更好地对伺服控制系统进行陷波器,滤去谐波振动,使得伺服控制系统具有更好的稳定性。
结合伺服控制系统,如图2所示,本发明提供了一种伺服控制系统谐振抑制的方法,该方法包括以下步骤:
步骤S1:设置伺服控制系统响应幅值的阈值;
伺服控制系统通常根据其具体情况,有一个可以接受的响应幅值,并且为了避免伺服控制系统发生谐振,通常该接受的系统响应幅值小于伺服谐振频率,因此,将伺服控制系统可接受的响应幅值作为伺服控制系统响应幅值的阈值,并记为H th ,将系统的谐振频率点处的幅值记为H max 。因此对于伺服控制系统,系统响应幅值超过阈值H th 的频率点构成伺服控制系统的谐振频率带,对于谐振频率带中的频率点都要进行谐振抑制。现有技术中,将无机械共振时谐振频率最大值所对应的幅值作为阈值,对于伺服控制系统而言,当未达到谐振频率时,系统可能已经产生较为严重的振动,因此,在确定该阈值时考虑伺服控制系统可以接受的响应幅值,并将该幅值作为阈值,更利于伺服系统保持稳定性。
步骤S2:采集速度环调节器输出的电流指令,对采集的电流指令进行频谱分析,并根据分析结果确定陷波器的陷波参数,在本发明的实施例中,所述陷波参数包括将陷波中心频率ω n 、陷波深度ξ和陷波宽度k。在本发明的优选实施例中,采用的是FFT方法对采集的电流指令进行频谱分析,得到FFT分析曲线,然后根据FFT(Fast Fourier Transform,快速傅里叶变换)分析曲线,确定陷波器的陷波参数。
在本发明中,采集的是速度环调节器输出的电流信号iq ref ,在伺服控制系统中,伺服电机转矩T e 计算公式为:
Figure DEST_PATH_IMAGE014
(1)
公式(1)中,p表示伺服电机转子极对数,ψ f 表示伺服电机转子磁链,i d 表示伺服电机d轴电流,i q 表示伺服电机q轴电流,L d 表示伺服电机d轴的电感,L q 表示伺服电机q轴的电感。在采用i d =0的控制方法时,电机的输出转矩只与伺服电机q轴电流i q 呈线性关系,当系统发生谐振时,电机受到谐振影响转速和转矩的输出信号都会发生振荡,因为转矩只与q轴电流有关,所以此时q轴电流也会发生明显的振荡现象,由于电流调节器的作用q轴电流i q 是紧跟着电流指令iq ref 变化,因此,可以看出,只需对电流指令iq ref 进行频谱分析就可以得到系统的谐振信息。而本发明摒弃了现有技术中直接利用速度差进行频谱分析的方法,而是利用电流信号进行频谱分析获得陷波参数;并且是在速度环调节器的输出端进行信号采集,对该输出端的电流信号进行频谱分析,利用分析结果确定陷波器的陷波参数,更能反映陷波器的陷波特性,从而使得根据所得参数确定的陷波器的传递函数,更符合伺服控制系统特性,更能有效地对伺服控制系统的谐振进行抑制。在确定陷波器的陷波参数时,如图3所示,可按照如下步骤进行:
采集速度调节环输出的电流信号,并统计采集到的数据量,如果数据量超过设定的最少数据量,则执行下一步骤,否则,继续采集;在优选的实施例中,最少数据量可设定为1024.
对采集的数据进行FFT(Fast Fourier Fransform,快速傅里叶变换)分析,并生成FFT分析曲线,根据FFT分析曲线得到陷波中心频率ω n 及相应的幅值H max。判断陷波中心频率对应的幅值是否大于伺服控制系统响应幅值的阈值,如果大于所述阈值,即H max > H th ,则根据FFT分析曲线,计算陷波宽度k和陷波深度ξ,并开启陷波器,否则,则不计算陷波宽度和陷波深度,因为此时没有发生谐振,因此也不会开启陷波器进行陷波。
陷波器的各陷波参数采用如下方式确定:
陷波器中心频率ω n :陷波器的陷波中心频率对应的是FFT分析曲线幅值最高点P 0处的频率,参见图4所示,即伺服控制系统的谐振频率,对应的振幅即为H max
陷波宽度k:表示以陷波器的陷波中心频率在陷波器作用下复制增益与初值相比小于3dB的频率带宽,具体计算时,本发明的实施例在计算陷波宽度k时,提取FFT分析曲线最高点两侧幅值均为所述阈值的两个频率点,所述两个频率点对应的频率分别为ω 1ω 2,即谐振频率点两侧幅值为H th 的两个频率点,如图4所示,P 1P 2,其中P 1点对应的频率为ω 1P 2点对应的频率为ω 2,并根据如下公式确定陷波宽度k
Figure 287929DEST_PATH_IMAGE015
(2)
式(2)中,ω n 表示陷波中心频率,
Figure 892085DEST_PATH_IMAGE016
表示中心频率ω n 与频率ω 1间的频带宽度,
Figure 594462DEST_PATH_IMAGE017
表示中心频率ω n 与频率ω 2间的频带宽度,m为正数表示将
Figure 159436DEST_PATH_IMAGE016
Figure 11985DEST_PATH_IMAGE018
中较大值放大m倍。其中n根据伺服控制系统的特性进行设置,因为伺服控制系统中,陷波中心频率ω n 和频率ω 1的之间的宽度与陷波中心频率ω n 和频率ω 2之间的宽度可能并不相等,因此,为了保证陷波器有足够的陷波宽度来覆盖谐振时响应幅度超出阈值的频带反馈,将陷波宽度k设置为
Figure 955801DEST_PATH_IMAGE016
Figure 563500DEST_PATH_IMAGE018
中较大值的m倍。根据系统情况将m的范围设定在1~20,在本发明的优选实施例中,m的范围取1~2。
陷波深度ξ:陷波深度ξ影响着陷波中心频率点处所能够提供的最大的幅值衰减,陷波中心频率处幅值衰减为20lg|ξ|dB,在本发明中,陷波深度
Figure 740404DEST_PATH_IMAGE020
采用下式进行计算:
Figure 803169DEST_PATH_IMAGE021
(3)
根据公式(3)可以看出,陷波深度
Figure 116338DEST_PATH_IMAGE020
的取值是
Figure 247935DEST_PATH_IMAGE023
的。
步骤S3:根据所述陷波器的陷波参数确定陷波器的传递函数。
陷波器的参数包括陷波中心频率、陷波深度和陷波宽度,陷波器的传递函数形式采用如下形式:
Figure 787501DEST_PATH_IMAGE024
(4)
式(4)中,ω n 表示陷波器的陷波中心频率,ξ表示陷波深度,k为陷波宽度,s表示复变量。
本发明采用陷波器传递函数,引入陷波中心频率、陷波深度及陷波宽度乘积作为传递函数分子中一次项的系数,把陷波中心频率与陷波宽度乘积作为传递函数分母一次项的系数,该传递函数所确定的陷波器能够中和谐振引起的幅值,从而达到更好的谐振抑制效果。,根据FFT分析曲线,在线确定了陷波器的陷波参数后,即可得到用于该伺服控制系统的陷波器的传递函数。因此相对于传统的陷波器,本发明的陷波器能够在线调节陷波参数,使得伺服控制系统具有更好的动态性能,同时是伺服系统具有更好的谐振抑制能力,从而在伺服系统的工作过程中,更好地抑制伺服系统的谐振。
应该说明的是,采用本发明的方法得到陷波器的陷波参数,因可以应用到其他形式的陷波器的传递函数中,例如专利201910327030.6中所提到的陷波器的传递函数的形式,但相对与现有技术中的陷波器的传递函数形式,采用本发明中所提的传递函数形式的陷波器,陷波性能更好。
通过对陷波器的传递函数(4)进行离散化,得到差分方程,即可实现陷波器的数字化,为得到传递函数对应的差分方程,通过Z变换的方法,具体如下:
令公式(4)中
Figure 575329DEST_PATH_IMAGE025
Figure 8584DEST_PATH_IMAGE026
,其中z表示复变量,f n 表示采样频率,T表示采样周期,n为正整数,因此有,
Figure 364610DEST_PATH_IMAGE027
(5)
将公式(5)化简并给出差分方程的形式有:
Figure 125893DEST_PATH_IMAGE028
(6)
式(6)中,x(n)表示n时刻输入的信号的值,y(n)表示n时刻输出的信号的值,ɑ 0ɑ 1b 0b 1b 2表示系数,其中ɑ 0ɑ 1b 0b 1b 2的具体表达式如下:
Figure DEST_PATH_IMAGE029
(7)
Figure 796040DEST_PATH_IMAGE030
(8)
Figure DEST_PATH_IMAGE031
(9)
Figure 818222DEST_PATH_IMAGE032
(10)
Figure 935695DEST_PATH_IMAGE033
(11)
步骤S4:将该传递函数应用于所述伺服控制系统的控制环路。
将确定的陷波器的传递函数应用于伺服控制系统的控制环路,具体步骤是:编码器采集伺服电机的转速,并将采集到的伺服电机的转速反馈到速度环调节器;速度环调节器将实际输入的转速与反馈回的伺服电机的实际转速做差,得到转速差集,然后将所述转速差集作为输入,并输出电流指令;FFT模块采集所述速度环调节器输出的所述电流指令,进行FFT分析,并得到陷波器的参数,并传递给所述陷波器,从而得到陷波器的传递函数;陷波器根据所述传递函数,对速度环调节器输出的电流信号进行陷波,然后输出,并通过电流环调节器,将电流信号转换为电压信号,所述电压信号驱动所述伺服电机。
将本发明的方法,应用于所述的伺服控制系统,从而使得伺服控制系统能够进行谐振抑制,在本发明的实施例中,伺服控制系统为了能够实施所述方法,FFT模块包括:
数据采集器,用于采集速度环调节器输出的电流信号iq ref
采集判断器,用于统计数据采集器采集的数据量,并判断是否满足设定的最少数据量;因为如果通过分析伺服控制系统的频率特性,必须采集足够多的信号频率点,在本发明的实施例中,将最少数据量设定为1024,即当采集判断器判断已经采集了1024个点后,FFT模块才可以继续执行其后续功能操作;
数据处理器,采用FFT(Fast Fourier Transform,FFT,快速傅里叶变换)作频谱分析,得到FFT分析曲线;
数据提取器,提取FFT分析曲线中最高点处的幅值和频率、提取最高点两侧幅值均为伺服控制系统响应幅值的阈值的两个频率点;即图4中P 0点的幅值和频率,以及P 1点和P 2点处的频率,其中P 1点和P 2点处的幅值与阈值相同;
第一计算器:根据公式
Figure 715432DEST_PATH_IMAGE034
计算陷波宽度k,其中m为正数,表示将
Figure 579483DEST_PATH_IMAGE035
Figure 987331DEST_PATH_IMAGE036
中较大值放大m倍,
Figure 292017DEST_PATH_IMAGE035
表示陷波中心频率ω n 与频率ω 1间的频带宽度,
Figure 637679DEST_PATH_IMAGE037
表示陷波中心频率ω n 与频率ω 2间的频带宽度,其中m为正数,范围为1~20;
第二计算器:根据公式
Figure 305420DEST_PATH_IMAGE038
计算陷波深度ξ,其中,H th 表示伺服控制系统响应幅值的阈值,H max表示FFT分析曲线幅值最高点的处幅值;
第三计算器:根据所述的中心频率、陷波深度和陷波宽度确定传递函数,传递函数形式为:
Figure 36616DEST_PATH_IMAGE039
,其中s为复变量。
对应用本发明的方法的伺服控制系统的进行试验得到图5~8所示结果。
图5给出了伺服控制系统未开启陷波器时系统波形图,从图5中可以看出,在未开启陷波器前,系统的转动波动较大。
图6给出了伺服控制系统未开启陷波器时FFT频谱图,从图6可以看出,此时系统中有一个频率为786Hz的谐振,因此系统的谐振较严重。
图7给出了伺服控制系统开启陷波器时系统波形图,从图7可以看出,开启陷波器后,系统的转动波动明显减小。
图8给出了伺服控制系统开启陷波器时FFT频谱图,从图8可以看出,开启陷波器后,系统没有产生明显的谐振,说明基于本发明的方法对伺服控制系统的谐振具有明显的抑制作用。
根据与专利201910327030.6对比,对比专利中,陷波器的输入是速度环的输出,而通过对速度环的输入做FFT分析,不能准确地确定陷波器的参数;而本发明通过对速度环的输出做频谱分析,由于速度环的输出陷波器的参数起最直接影响,因此相比现有技术,本发明更能准确地识别出符合伺服系统特性的陷波参数。此外,陷波器的陷波参数是根据频谱分析结果进行计算的,该专利在进行参数确定时,是以谐振频率最大值所对应的幅值作为阈值的,但是对于伺服控制系统,当未达到谐振频率时,系统已经会产生较为严重的振动,而本发明,考虑伺服控制系统可以接受的响应幅值,从而确定出适合伺服控制系统的陷波器参数,采用该参数值确定的陷波器的传递函数,更符合系统特性,更能对伺服系统进行更好地陷波。
本实施所提供的伺服控制系统谐振在线抑制的方法,通过采集速度环调节器输出的电流信号,并对该电流信号进行频谱分析,根据频谱分析结果确定陷波器的陷波参数,实现了陷波器的陷波参数在线调整,并根据所得陷波参数,确定陷波器的传递函数,使得陷波器具有更好的动态特性,从而能够更好的抑制伺服控制系统的谐振。
上述实施例仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和等同替换,这些对本发明权利要求进行改进和等同替换后的技术方案,均落入本发明的保护范围。

Claims (7)

1.一种伺服控制系统谐振在线抑制的方法,其特征在于,包括:
设置伺服控制系统响应幅值的阈值;
采集速度环调节器输出的电流指令,采用FFT对采集的电流指令进行频谱分析,并根据分析结果确定陷波器的陷波参数;
根据所述陷波器的陷波参数获取陷波器的具体传递函数;所述陷波器的陷波参数包括陷波中心频率、陷波深度和陷波宽度,传递函数形式为:
Figure DEST_PATH_IMAGE001
式中,ω n 表示陷波器的陷波中心频率,ξ表示陷波深度,k表示陷波宽度,s表示复变量;所述陷波宽度的确定方法为:提取FFT分析曲线幅值最高点两侧幅值均为所述阈值的两个频率点,所述两个频率点对应的频率分别为ω 1ω 2 ,并根据如下公式确定陷波宽度:
Figure 462994DEST_PATH_IMAGE002
式中,ω n 表示陷波中心频率,
Figure 862358DEST_PATH_IMAGE003
表示中心频率ω n 与频率ω 1间的频带宽度,
Figure 110936DEST_PATH_IMAGE004
表示中心频率ω n 与频率ω 2间的频带宽度,m为正数,表示将
Figure 974987DEST_PATH_IMAGE005
Figure 710731DEST_PATH_IMAGE004
中较大值放大m倍,m的范围为1~20;
将该传递函数应用于所述伺服控制系统的控制环路。
2.根据权利要求1所述的一种伺服控制系统谐振在线抑制的方法,其特征在于,所述陷波中心频率的确定方法为:采集FFT分析曲线幅值最高点的处的频率作为陷波器的陷波中心频率。
3.根据权利要求1所述的一种伺服控制系统谐振在线抑制的方法,其特征在于,所述陷波深度的确定方法为:
Figure 267614DEST_PATH_IMAGE006
式中,H th 表示伺服控制系统响应幅值的阈值,H max 表示FFT分析曲线幅值最高点的处幅值。
4.根据权利要求1所述的一种伺服控制系统谐振在线抑制的方法,其特征在于,将所述
Figure 206751DEST_PATH_IMAGE005
Figure 625226DEST_PATH_IMAGE004
中较大值放大m倍时,m的范围为1~2。
5.根据权利要求1-4任一项所述的一种伺服控制系统谐振在线抑制的方法,其特征在于,确定所述陷波器的陷波参数的步骤为:
采集速度调节环输出的电流信号,并统计采集到的数据量,如果数据量超过设定的最少数据量,则执行下一步骤,否则,继续采集;
对采集的数据进行FFT分析,并生成FFT分析曲线,根据FFT分析曲线得到陷波中心频率及相应的幅值;
判断陷波中心频率对应的幅值是否大于伺服控制系统响应幅值的阈值,如果大于所述阈值,则根据FFT分析曲线,计算陷波宽度和陷波深度,并开启陷波器,否则,则不计算陷波宽度和陷波深度。
6.一种伺服控制系统,其特征在于,包括:
伺服电机;
编码器,用于采集伺服电机的实际转速,并将采集到的伺服电机的实际转速进行反馈;
速度环调节器,用于将实际输入的转速与所述编码器反馈回的实际转速做差,并输出电流指令;
自适应陷波器,用于对所述速度环调节器输出的电流指令进行陷波,滤除谐振分量,并输出陷波后的电流指令;
FFT模块,用于采集所述速度环调节器输出的所述电流指令,进行FFT分析,并根据FFT分析曲线,得到陷波器的陷波参数,并将所述陷波参数传递给所述陷波器;
电流环调节器,将陷波器输出的电流信号转化为电压信号,所述电压信号驱动所述伺服电机。
7.根据权利要求6所述的一种伺服控制系统,其特征在于,所述FFT模块包括:
数据采集器,用于采集速度环调节器输出的电流信号;
数据记录器,用于统计所述数据采集器采集的数据量,并判断是否满足需要的最少数据量;
数据处理器,采用FFT作频谱分析,得到FFT分析曲线;
数据提取器,提取FFT分析曲线中最高点处的频率作为陷波中心频率ω n ,以及对应的幅值,提取最高点两侧幅值均为伺服控制系统响应幅值的阈值的两个频率点对应的频率ω 1ω 2
第一计算器:根据公式
Figure 762946DEST_PATH_IMAGE007
计算陷波宽度k,其中m为正数,表示将
Figure 225151DEST_PATH_IMAGE005
Figure 448322DEST_PATH_IMAGE008
中较大值放大m倍,
Figure 903443DEST_PATH_IMAGE005
表示陷波中心频率ω n 与频率ω 1间的频带宽度,
Figure 161249DEST_PATH_IMAGE008
表示陷波中心频率ω n 与频率ω 2间的频带宽度,其中m为正数,范围为1~20;
第二计算器:根据公式
Figure 59935DEST_PATH_IMAGE009
计算陷波深度ξ,其中,H th 表示伺服控制系统响应幅值的阈值,H max表示FFT分析曲线幅值最高点的处幅值;
第三计算器:根据所述的陷波中心频率、陷波深度和陷波宽度确定传递函数,传递函数的形式为:
Figure 770402DEST_PATH_IMAGE010
,其中s为复变量。
CN202010361439.2A 2020-04-30 2020-04-30 一种伺服控制系统谐振在线抑制的方法及伺服控制系统 Pending CN111262496A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010361439.2A CN111262496A (zh) 2020-04-30 2020-04-30 一种伺服控制系统谐振在线抑制的方法及伺服控制系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010361439.2A CN111262496A (zh) 2020-04-30 2020-04-30 一种伺服控制系统谐振在线抑制的方法及伺服控制系统

Publications (1)

Publication Number Publication Date
CN111262496A true CN111262496A (zh) 2020-06-09

Family

ID=70955200

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010361439.2A Pending CN111262496A (zh) 2020-04-30 2020-04-30 一种伺服控制系统谐振在线抑制的方法及伺服控制系统

Country Status (1)

Country Link
CN (1) CN111262496A (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111650889A (zh) * 2020-06-23 2020-09-11 佛山市普拉迪数控科技有限公司 一种五轴龙门数控加工中心的切削颤振抑制智能控制方法
CN111896039A (zh) * 2020-07-24 2020-11-06 中国水利水电科学研究院 一种便携式安全监测数据采集和诊断系统
CN111934591A (zh) * 2020-07-16 2020-11-13 武汉久同智能科技有限公司 一种基于自适应陷波的机械谐振抑制方法
CN112886896A (zh) * 2021-03-10 2021-06-01 无锡信捷电气股份有限公司 一种伺服系统用在线自适应陷波滤波器实现方法
CN114083336A (zh) * 2021-12-17 2022-02-25 广东技术师范大学 一种驱动控制器的机械共振消除方法、装置及系统
CN114329728A (zh) * 2021-12-30 2022-04-12 苏州汇川控制技术有限公司 伺服系统共振优化方法、设备及介质
CN114825415A (zh) * 2022-05-17 2022-07-29 南京南瑞继保电气有限公司 一种柔性直流高频谐振的自适应抑制方法及系统
CN116505810A (zh) * 2023-06-30 2023-07-28 苏州铼钠克信息技术有限公司 陷波器使用及伺服系统振荡抑制方法、装置、设备及介质

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1078958A (zh) * 1992-05-13 1993-12-01 金星产电株式会社 用于控制电梯速度的装置和方法
US20080309279A1 (en) * 2007-06-14 2008-12-18 Toru Tazawa Motor controller
CN107024617A (zh) * 2017-04-21 2017-08-08 上海电气集团股份有限公司 一种电机转矩振荡频率辨识方法
CN107394784A (zh) * 2017-08-30 2017-11-24 无锡信捷电气股份有限公司 一种陷波滤波器在线谐振抑制参数优化方法
CN110350886A (zh) * 2019-07-12 2019-10-18 广东工业大学 一种永磁伺服系统的谐振频率检测方法、系统及装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1078958A (zh) * 1992-05-13 1993-12-01 金星产电株式会社 用于控制电梯速度的装置和方法
US20080309279A1 (en) * 2007-06-14 2008-12-18 Toru Tazawa Motor controller
CN107024617A (zh) * 2017-04-21 2017-08-08 上海电气集团股份有限公司 一种电机转矩振荡频率辨识方法
CN107394784A (zh) * 2017-08-30 2017-11-24 无锡信捷电气股份有限公司 一种陷波滤波器在线谐振抑制参数优化方法
CN110350886A (zh) * 2019-07-12 2019-10-18 广东工业大学 一种永磁伺服系统的谐振频率检测方法、系统及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
龙丁: "交流伺服系统谐振陷波器参数自整定研究", 《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111650889A (zh) * 2020-06-23 2020-09-11 佛山市普拉迪数控科技有限公司 一种五轴龙门数控加工中心的切削颤振抑制智能控制方法
CN111934591A (zh) * 2020-07-16 2020-11-13 武汉久同智能科技有限公司 一种基于自适应陷波的机械谐振抑制方法
CN111934591B (zh) * 2020-07-16 2022-08-02 武汉久同智能科技有限公司 一种基于自适应陷波的机械谐振抑制方法
CN111896039A (zh) * 2020-07-24 2020-11-06 中国水利水电科学研究院 一种便携式安全监测数据采集和诊断系统
CN112886896A (zh) * 2021-03-10 2021-06-01 无锡信捷电气股份有限公司 一种伺服系统用在线自适应陷波滤波器实现方法
CN114083336A (zh) * 2021-12-17 2022-02-25 广东技术师范大学 一种驱动控制器的机械共振消除方法、装置及系统
CN114329728A (zh) * 2021-12-30 2022-04-12 苏州汇川控制技术有限公司 伺服系统共振优化方法、设备及介质
CN114825415A (zh) * 2022-05-17 2022-07-29 南京南瑞继保电气有限公司 一种柔性直流高频谐振的自适应抑制方法及系统
CN114825415B (zh) * 2022-05-17 2023-08-08 南京南瑞继保电气有限公司 一种柔性直流高频谐振的自适应抑制方法及系统
CN116505810A (zh) * 2023-06-30 2023-07-28 苏州铼钠克信息技术有限公司 陷波器使用及伺服系统振荡抑制方法、装置、设备及介质
CN116505810B (zh) * 2023-06-30 2023-08-22 苏州铼钠克信息技术有限公司 陷波器使用及伺服系统振荡抑制方法、装置、设备及介质

Similar Documents

Publication Publication Date Title
CN111262496A (zh) 一种伺服控制系统谐振在线抑制的方法及伺服控制系统
CN105547698B (zh) 滚动轴承的故障诊断方法及装置
CN107167306B (zh) 基于阶次提取的旋转机械转子运行状态模态分析方法
CN111044814B (zh) 一种对变压器直流偏磁异常进行辨识的方法及系统
US6014617A (en) Method and apparatus for extracting a fundamental frequency based on a logarithmic stability index
CN110907826B (zh) 一种基于卷积神经网络滤波的电机故障诊断方法及系统
CN102353500B (zh) 一种用于动平衡测量的不平衡信号提取方法
CN114061678B (zh) 一种科氏流量计数字驱动方法
Blasco et al. Speed measurement of inverter fed induction motors using the FFT and the rotor slot harmonics
CN108388839A (zh) 一种基于二阶同步提取变换的强转速波动特征提取方法
CN113901379A (zh) 一种边缘端的实时数据动态在线快速处理方法
CN112782421B (zh) 一种基于音频的转速识别方法
CN114710081A (zh) 基于扩张状态观测器与改进型陷波器的在线谐振抑制方法
CN110716534A (zh) 一种基于自整定变分模态分解的工业过程振荡检测方法
CN105044459B (zh) 一种谐波分析方法
CN113643679A (zh) 基于级联滤波器的旋翼和尾桨气动噪声分离方法
CN117686232A (zh) 一种燃气轮机振动基频实时提取方法、装置及存储介质
CN115356108B (zh) 一种调制高阶水平挤压变换机械故障诊断方法与装置
CN113252929B (zh) 转速确定方法、装置、电子设备及计算机可读存储介质
CN106485073A (zh) 一种磨床故障诊断方法
CN109774389B (zh) 一种基于频谱分析的间接式胎压监测方法
CN110779723B (zh) 一种基于霍尔信号的变速工况电机轴承精确故障诊断方法
JP3506449B2 (ja) 騒音制御装置
CN111181528A (zh) 一种同步基频提取滤波器
CN103344990A (zh) 基于matlab电磁式可控震源输出信号的分析方法

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200609