CN108279571A - 一种有限转角机电伺服系统的模型参数辨识方法 - Google Patents

一种有限转角机电伺服系统的模型参数辨识方法 Download PDF

Info

Publication number
CN108279571A
CN108279571A CN201810108306.7A CN201810108306A CN108279571A CN 108279571 A CN108279571 A CN 108279571A CN 201810108306 A CN201810108306 A CN 201810108306A CN 108279571 A CN108279571 A CN 108279571A
Authority
CN
China
Prior art keywords
servo system
signal amplitude
electromechanical servo
amplitude
model
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.)
Granted
Application number
CN201810108306.7A
Other languages
English (en)
Other versions
CN108279571B (zh
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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201810108306.7A priority Critical patent/CN108279571B/zh
Publication of CN108279571A publication Critical patent/CN108279571A/zh
Application granted granted Critical
Publication of CN108279571B publication Critical patent/CN108279571B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

一种有限转角机电伺服系统的模型参数辨识方法,本发明涉及有限转角机电伺服系统的模型参数辨识方法。本发明的目的是为了解决现有开环扫频测定法使用时,机电伺服系统运动存在位置漂移、方向不确定,可能会使机电伺服系统位置转动角度过大,进而超过位置限制,引起机械结构的损坏的问题。过程为:一、建立机电伺服系统机理模型;二、获得实际控制信号幅值;三、计算补偿信号幅值;四、改变输入信号幅值;五、重复二到四,直到得到理想控制信号幅值,获得频率点的幅相特性;六、改变输入信号频率,重复二到五,得到所需频带内各频率点的幅相特性;七、得到机电伺服系统的模型参数。本发明用于机电伺服系统的模型参数辨识领域。

Description

一种有限转角机电伺服系统的模型参数辨识方法
技术领域
本发明涉及有限转角机电伺服系统的模型参数辨识方法。
背景技术
随着电机种类的增多和制造水平的提高,电机已不仅仅作为简单的动力输出设备为生产服务,在电力电子和控制科学水平不断提高的今天,对电机的高性能控制已经可以为精密加工、自动控制、动态跟踪等高精尖领域服务。这不仅促进了电机制造技术的发展,同时也为机电一体化控制的研究提供了广阔舞台。对于高性能的机电控制系统,通常需要建立负反馈回路以使电机输出跟踪参考输入信号,通过这种方法对电机进行精确的控制。这样的机电控制系统被称为机电伺服系统。在设计机电伺服系统的数字控制器时,一般采用PID控制方法或频域校正法,其中频域校正法具有较高的动态性能,能够适应高精度、高动态的控制任务要求。频域校正法的设计需要较为精确的机电伺服系统数学模型参数,使校正后的系统达到频带要求和“双十”指标,获得较好的动静态性能。
在工程上常用的机电伺服系统数学模型参数的确定方法为开环扫频测定法,其主要思想是:给伺服系统加入正弦输入激励信号,则系统的输出响应也是同一频率的正弦信号,通过改变输入信号的频率,对输入的信号和输出端信号进行数据处理,可获得被测系统或对象的频率特性,进而通过最小二乘法拟合得到系统的开环传递函数。这种测定方法虽然能够得到较为准确的机电伺服系统数学模型参数,但是在实际应用时,由于三阶伺服系统输入正弦信号的响应存在稳态分量和系统中存在各种非线性因素的影响,开环扫频存在位置漂移、方向不确定的问题,应用在转角有限的机电伺服系统上时,可能会机电伺服系统的位置转动角度过大,进而超过位置限制,对电机和限位模块造成冲击,严重时可能会引起机械结构的损坏。
发明内容
本发明的目的是为了解决现有的开环扫频测定法由于非线性因素和正弦响应存在的稳态分量的影响导致机电伺服系统运动存在位置漂移、方向不确定,可能会使机电伺服系统位置转动角度过大,进而超过位置限制,引起机械结构的损坏的问题,而提出一种有限转角机电伺服系统的模型参数辨识方法。
一种有限转角机电伺服系统的模型参数辨识方法具体过程为:
步骤一、根据机电伺服系统的结构建立机电伺服系统机理模型,即机电伺服系统开环传递函数模型G(s);
步骤二、在机电伺服系统闭环的状态下,输入信号幅值为Ar、频率为f的正弦信号,使用相关辨识法获得该频率点上的实际控制信号幅值Au
步骤三、建立步骤二中输入信号幅值Ar与对应频率点上的实际控制信号幅值Au的线性关系,根据所需理想控制信号幅值(经验值)计算补偿信号幅值ΔA;
步骤四、根据步骤三中得到的补偿信号幅值ΔA,将输入信号幅值从Ar渐变到实际信号幅值Ar′;
步骤五、令步骤二中的Ar=Ar′,重复步骤二到四,直到实际控制信号幅值Au与理想控制信号幅值相等,对此时的实际控制信号u和位置信号θ进行FFT数据处理,获得机电伺服系统模型在频率点f上的幅相特性;
所述幅相特性包括幅值比和相角差
步骤六、改变输入信号频率f,重复步骤二到五,得到所需频带内各频率点的幅相特性;
步骤七、使用最小二乘法对步骤六中得到的各频率点的幅相特性数据进行拟合,得到步骤一中的机电伺服系统开环传递函数模型的参数τe、τm和K;
τe为电磁时间常数,τm为机电时间常数,K为标称模型等效增益。
本发明的有益效果为:
本发明提出了一种有限转角机电伺服系统的模型参数辨识方法。为了实现上述目的,本发明的一种有限转角机电伺服系统的模型参数辨识方法为:在机电伺服系统闭环的状态下,通过建立输入信号与使用相关辨识法获得的实际控制信号之间的线性关系获得与开环扫频一样的控制信号,进而进行闭环扫频获得机电伺服系统的幅相特性,使用最小二乘法拟合得到机电伺服系统的模型参数。
相较于现有技术,本发明一种有限转角机电伺服系统的模型参数辨识方法在机电伺服系统转角范围有限时,步骤一使用相关辨识法获得该频率点上的实际控制信号幅值,对比现有的理论计算得到控制信号幅值具有实时性高、辨识准确的优点;步骤三将输入信号幅值渐变到实际信号幅值,对比直接改变输入信号幅值,避免了信号突变引起的加速度过大,对系统造成冲击过大的问题,具有安全性高的优点;本发明设计的闭环扫频方法结构简单,能够有效地控制机电伺服系统位置转动角度在初始位置附近,并且获得与开环扫频测定法一样的结果,具有设计方便、实现简单、安全性高的优点,解决了现有的开环扫频测定法在使用时机电伺服系统位置转动角度存在位置漂移、方向不确定,可能会使机电伺服系统位置转动角度过大,进而超过位置限制,引起机械结构的损坏的问题。表1验证了本发明提出的模型参数辨识方法能够准确的获得机电伺服系统在各频率点上的幅相特性,分别用开环扫频法和本发明提出的闭环扫频法对转台这一典型的机电伺服系统进行扫频,获得其1到10Hz的幅相特性,可以看出,本发明提出的闭环扫频法获得的系统在各频率点上的幅相特性与开环扫频法得到的幅相特性有略微的差异,但是考虑到实际系统中存在的不确定性和实验结果的偶然误差,这些差异都在误差允许的范围内;图2说明了本发明提出的闭环扫频法能够有效克服现有开环扫频法存在的位置漂移问题,向转台输入3V,1Hz的正弦控制信号,使用开环扫频法时,转台从初始位置183°运动到了155°(图2中实线),有较大的位置转动,而本发明有效地控制了机电伺服系统位置转动角度在初始位置155°附近(图2中虚线);图3说明了本发明提出的方法能够避免了信号突变引起的加速度过大的问题,在获得2Hz点的频率特性时,若在2.5s时直接将所需的补偿信号叠加到输入信号上,实际控制信号有较大的波动(图3中虚线),对转台伺服系统的冲击较大,而使用本发明提出的在输入信号上叠加渐变的补偿信号的方法(图3中实线),实际控制信号则比较平滑,不会有较大突变。
附图说明
图1为本发明机电伺服系统的结构机理图,其中:u为控制信号,s为拉普拉斯算子,iq为交轴电流,ω为电机机械角速度,θ为机电伺服系统实际位置输出;re为电机等效电阻;kps为驱动器放大系数;kV为电流放大系数;ke为反电势系数;kT为电机力矩系数;L为电枢电感;JΣ为轴系总转动惯量;Te为驱动器输出力矩;Tl为干扰力矩,包含摩擦力矩及波动力矩;
图2为本发明实施例中转台伺服控制系统在使用现有开环扫频法和本发明提出的闭环扫频法获取1Hz点幅相特性时的位置输出对比图;
图3为本发明实施例中转台伺服控制系统在输入信号上直接叠加补偿信号和本发明提出的渐变叠加补偿信号时的实际控制信号对比图。
具体实施方式
具体实施方式一:本实施方式的一种有限转角机电伺服系统的模型参数辨识方法具体过程为:
本发明主要针对转角有限的机电伺服系统在进行开环正弦扫频获取模型参数时,机电伺服系统位置转动角度存在位置漂移、方向不确定,可能超过限位的问题,从而提出的一种有限转角机电伺服系统的模型参数辨识方法。
步骤一、根据机电伺服系统的结构建立机电伺服系统机理模型,即机电伺服系统开环传递函数模型G(s);
步骤二、在机电伺服系统闭环的状态下,输入信号幅值为Ar、频率为f的正弦信号,使用相关辨识法获得该频率点上的实际控制信号幅值Au
步骤三、建立步骤二中输入信号幅值Ar与对应频率点上的实际控制信号幅值Au的线性关系,根据所需理想控制信号幅值(经验值)计算补偿信号幅值ΔA;
步骤四、根据步骤三中得到的补偿信号幅值ΔA,将输入信号幅值从Ar渐变到实际信号幅值Ar′;Ar′=Ar+ΔA;
步骤五、令步骤二中的Ar=Ar′,重复步骤二到四,直到实际控制信号幅值Au与理想控制信号幅值相等,对此时的实际控制信号u和位置信号θ进行FFT数据处理,获得机电伺服系统模型在频率点f上的幅相特性;
所述幅相特性包括幅值比和相角差
步骤六、改变输入信号频率f,重复步骤二到五,得到所需频带内各频率点的幅相特性;
步骤七、使用最小二乘法对步骤六中得到的各频率点的幅相特性数据进行拟合,得到步骤一中的机电伺服系统开环传递函数模型的参数τe、τm和K;
τe为电磁时间常数,τm为机电时间常数,K为标称模型等效增益。
具体实施方式二:本实施方式与具体实施方式一不同的是:所述步骤一中根据机电伺服系统的结构(图1)建立机电伺服系统机理模型,即机电伺服系统开环传递函数模型G(s);具体过程为:
机电伺服系统机理模型为:
式中,id,iq为同步旋转d-q坐标系中的直轴电流和交轴电流,d为直轴,q为交轴;ωr为电机机械角速度;re为电机等效电阻;L为电机电枢电感;JΣ为轴系总转动惯量;D为电机旋转时的摩擦系数;Tl为干扰力矩,包含摩擦力矩及波动力矩;Pm为电机极对数;ψf为转子永磁效应对应的每对磁极磁通;ud,uq为同步旋转d-q坐标系中的直轴电压和交轴电压;为微分算子;
上式(1)中建立了d轴和q轴的动态方程,但实际对电机进行矢量控制时,通常忽略机电伺服系统机理模型中id的影响,通过控制iq来直接控制电机的力矩输出,即
根据式(2),将电机电流反馈放大系数kf设为1,并忽略电机旋转时的摩擦系数D,令反电势系数电机力矩系数机电伺服系统机理模型变为如图1所示;
则由实际控制信号u到位置输出θ(单位为rad)的机电伺服系统开环传递函数模型为:
式中,G(s)为机电伺服系统开环传递函数模型,s为拉普拉斯算子,kps为电机驱动器放大系数,kV为电机电流放大系数;记为电磁时间常数,JΣ为轴系总转动惯量,为机电时间常数,为标称模型等效增益,式(3)简化为:
对于一般的机电伺服系统,总能够满足τm>>τe,所以上式(4)又简化为:
至此,机电伺服系统开环传递函数模型G(s)建立完毕;
在机电伺服系统开环传递函数模型G(s)基础上使用开环扫频法,会存在一个位置漂移,下面对位置漂移产生原理进行解释:
假设系统输入信号为:
r=Arsin(ωt) (6)
其中,Ar和ω分别是系统输入信号r的幅值和角频率,ω=2πf,t为时间,则系统输入信号r的拉普拉斯变换为:
则系统输出信号的拉普拉斯变换为:
将系统输出信号进行反拉普拉斯变换,得到系统输出信号为:
当系统输出进入稳态时,即t→∞时,式(9)变为:
即响应输出中包含一个稳态分量和一个正弦波动,因此在给机电伺服系统输入正弦指令信号时,会存在一个位置漂移,大小为方向与初始运动方向一致。然而在实际系统中,由于摩擦力矩、偏载力矩等非线性因素的影响,这种位置漂移的方向和大小都存在不确定性,因此当机电伺服系统转角范围有限时,这种不确定性可能会导致机电伺服系统位置运动角度过大,超过转角限制范围,进而造成危险,对机械结构造成损坏。
因此,有必要设计一种即能够控制扫频时机电伺服系统的转角范围,又不影响扫频获得的机电伺服系统开环频率特性的方法,解决机电伺服系统在进行开环正弦扫频获取模型参数时,出现的位置漂移、方向不确定,可能超过限位的问题。
其它步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:所述步骤二中在机电伺服系统闭环的状态下,输入信号幅值为Ar、频率为f的正弦信号,使用相关辨识法获得该频率点上的实际控制信号幅值Au
具体过程为:
设计一个简单比例控制器,将机电伺服系统闭环运行,在能够保证系统稳定的前提下,控制器的结构和控制精度不影响本发明后续步骤的进行。
虽然在理论上能够对频率点f上的实际控制信号幅值Au进行计算,但是当系统闭环后,闭环控制器的输入变为输入信号r与实际位置θ之间的偏差,因此实际控制信号u会与理论上有较大的差别,即使只考虑频率点f上的信号分量,由于机电伺服系统中存在的非线性因素影响,也与理论有差别,因此需要使用相关辨识法得到实际的频率点f上的实际控制信号幅值Au。首先介绍一下相关辨识法的原理:
相关辨识是利用一对正交的基准函数计算它们与需滤波数据(函数)之间的相关函数,从而辨识出所需频率分量的幅值及相位。
设所要处理的采样信号S(t)的数学描述表示为
其中,t表示时间;是采样信号中的有用分量,即基波信号,As、ωs分别是基波信号的幅值、角频率和相位;代表采样信号中的中、低频强干扰,其中n为有限正整数,表示这部分干扰等效分布在有限频率点上,Aj、ωj分别是第j项中、低频强干扰的幅值、角频率和相位,j=0,1,2,3,...n;vs(t)表示随机噪声,用于描述采样信号中的高频噪声。
设所要处理的基波信号为:
其中,Re、Im分别表示正弦分量和余弦分量的系数。
取两个基准信号fs01(t)和fs02(t):
fs01(t)=F0sinωst (13)
fs02(t)=F0cosωst (14)
fs(t)与fs01(t)和fs02(t)的互相关函数R1(τ)和R2(τ)用下式计算:
其中,Ns为采样点个数,正整数,Ts为采样时间,τ为积分变量,当τ=0时:F0为基准信号的幅值;
将式(12)和(13)带入式(17),则:
因为则:
同理可求:
由式(20)和式(21)可得:
由式(22)和式(23)可知,只要求出采样信号S(t)与基准信号fs01(t)和fs02(t)之间的相关函数R1(τ)和R2(τ)在τ=0点的值,即求出R1(0)和R2(0)就可计算出所求某个频率信号的幅值和相位。
计算相关函数要用到采样信号的数据,由于采样信号中不仅仅含有它的基波信号,因此要考虑信号中直流分量、谐波分量及异频分量的影响。很容易验证直流分量与基准信号的相关函数为零,谐波分量与基准信号的相关函数亦为零。对于异频分量,只要NsTs是基频周期及异频周期的公倍数,则其相关函数也为零。由上可见,除了基波分量之外,其它信号对所求的R1(0)和R2(0)均无影响。因此可以直接使用采样数据计算信号基波与基准信号之间的相关函数。
由于计算机采样的数据是离散形式的,因此计算相关函数也要用离散形式,因此,式(17)和(18)的离散形式表示为:
其中,NR为取样点数,k=0,1,2,3,...NR-1。
至此,相关辨识法的原理介绍完毕。
利用上述相关辨识法的原理,可以获得实际控制信号u中频率点f上的实际控制信号幅值Au,具体流程为:
步骤二一、根据输入信号的频率f和采样时间Ts,选取及合适的F0(经验值),得到两个基准信号分别为fs01(k)=F0sin(2πf·kTs)和fs02(k)=F0cos(2πf·kTs);
步骤二二、当k=0,1,2,3,...NR-1时,根据式(24)和(25)计算相关函数R1(0)和R2(0),其中fs(k)=u(t),即执行计算时的实际控制信号u,t为时间;
fs(k)为基波信号;
步骤二三、根据式(22)计算实际控制信号幅值
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是:所述步骤三中建立步骤二中输入信号幅值Ar与实际控制信号幅值Au的线性关系,根据所需理想控制信号幅值计算补偿信号幅值ΔA;具体过程为:
根据步骤二得到的实际控制信号幅值Au,由于闭环控制器是线性的,得到输入信号幅值Ar与实际控制信号幅值Au的对应关系,根据所需理想控制信号幅值计算补偿信号幅值ΔA为:
其中,理想控制信号幅值与开环扫频测定法时的输入信号一致。
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是:所述理想控制信号幅值按照经验选取,一般为3V~5V电压信号。
其它步骤及参数与具体实施方式一至四之一相同。
具体实施方式六:本实施方式与具体实施方式一至五之一不同的是:所述步骤四中根据步骤三中得到的补偿信号幅值ΔA,将输入信号幅值从Ar渐变到实际信号幅值Ar′;具体过程为:
通过步骤三得到补偿信号幅值ΔA后,若直接将输入信号幅值修改到Ar′=Ar+ΔA,则会出现输入信号斜率突变的情况,这会引起加速度的突变,对系统造成冲击。因此,将输入信号的幅值在1秒内按式(27)所示从输入信号幅值Ar逐渐增加到实际信号幅值Ar′;
其中rA(t)为补偿的输入信号,tsys为系统运行时间,tA为补偿开始的时间,tp为渐变时间。ΔA+Ar=Ar′,实际信号幅值Ar′对应的信号为rA(t)。
其它步骤及参数与具体实施方式一至五之一相同。
具体实施方式七:本实施方式与具体实施方式一至六之一不同的是:所述步骤五中令步骤二中的Ar=Ar′,重复步骤二到四,直到实际控制信号幅值Au与理想控制信号幅值相等,对此时的实际控制信号u和位置信号θ进行FFT数据处理,获得机电伺服系统模型在频率点f上的幅相特性;具体过程为:
由于实际机电伺服系统中存在摩擦力矩、偏载力矩等非线性因素的影响,当输入信号幅值Ar较小时,受非线性因素影响较大,实际控制信号u中频率点f上的实际控制信号幅值Au不能够与Ar呈完全线性关系,因此要多次修改实际信号幅值Ar′的大小,才能使实际控制信号幅值Au与理想控制信号幅值相等。
通过几次重复步骤二到四即可达到上述目的。
对此时的实际控制信号u和位置信号θ进行FFT数据处理,具体流程为:
步骤五一、对实际控制信号u进行FFT数据处理,得到其在频率点f上的幅值Au′和相角
步骤五二、对位置信号θ进行FFT数据处理,得到其在频率点f上的幅值Aθ′和相角
步骤五三、获得机电伺服系统模型在频率点f上的幅值比为Af=Aθ′/Au′、相角差为
其它步骤及参数与具体实施方式一至六之一相同。
采用以下实施例验证本发明的有益效果:
实施例一:
本实施例一种有限转角机电伺服系统的模型参数辨识方法具体是按照以下步骤制备的:
飞行仿真转台作为典型的机电伺服系统可以用来对多种运动形式进行模拟和跟踪,在航空航天、军工武器等需要高精度、高可靠性的领域具有重要地位。为验证本发明的有效性,以某型两轴飞行仿真转台的俯仰轴为被控对象进行实验验证。分别用开环扫频测定法和本发明提出的一种有限转角机电伺服系统的模型参数辨识方法获得主轴的各频率点幅相特性,验证本发明中闭环扫频方法的有效性。具体步骤如下:
①根据经验设计一超前校正控制器,使系统能够闭环稳定,本实施例中设计的超前校正控制器为:
将控制器进行离散化嵌入到转台的控制程序中(采样周期为0.5ms,采用双线性离散化方法);
②选取闭环扫频所需参数为:待扫频的频率范围为f=1,2...10Hz,基准信号为fs01(t)=2sin(2πf·t)和fs01(t)=2cos(2πf·t),理想控制信号幅值采样时间为Ts=0.0005s,每个频率的采样点数为然后编写扫频程序在实际系统中进行实施;
③令频率f=1,2...10Hz,转台伺服系统进行闭环运动,当实际控制信号u在频率点f上的实际控制信号幅值Au与理想控制信号幅值相等时,对此时的实际控制信号u和位置信号θ进行FFT数据处理,并记录转台模型在频率点f上的幅相特性(幅值比、相角差);
④使用最小二乘法对得到的各频率点的FFT数据进行拟合,得到最终的转台伺服系统的模型参数τe、τm和K。
为验证本发明提出的模型参数辨识方法能够有效克服现有开环扫频法存在的位置漂移问题,分别用开环扫频法和本发明提出的闭环扫频法对转台进行扫频,获得其1到10Hz的幅相特性,得到实验结果如表1所示。可以看出,本发明提出的闭环扫频法获得的转台伺服系统在各频率点上的幅相特性与开环扫频法得到的幅相特性有略微的差异,但是考虑到实际系统中存在的不确定性和实验结果的偶然误差,这些差异都在误差允许的范围内,可以说本发明提出的方法能够准确的获得转台伺服系统在各频率点上的幅相特性。
通过表1得到的闭环扫频数据,使用最小二乘法进行对象拟合,得到最终的转台伺服系统的开环传递函数模型参数τe=0.0011、τm=5.2792和K=1300.31。
以输入3V,1Hz的正弦控制信号为例,将两种方法在实施时的转台实际位置曲线进行对比,如图2所示,可以看出,使用开环扫频法时,转台从初始位置183°运动到了155°(图2中实线),有较大的位置转动,而本发明有效地控制了转台位置转动角度在初始位置155°附近(图2中虚线)。以2Hz实验为例,若在2.5s时直接将所需的补偿信号叠加到输入信号上,画出此时实际控制信号u(图3中虚线)可以看出,控制信号有较大的波动,对转台伺服系统的冲击较大,而使用本发明提出的在输入信号上叠加渐变的补偿信号的方法(图3中实线),实际控制信号则比较平滑,不会有较大突变。
表1开环扫频法与本发明提出的闭环扫频法实验对比
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (7)

1.一种有限转角机电伺服系统的模型参数辨识方法,其特征在于:所述方法具体过程为:
步骤一、根据机电伺服系统的结构建立机电伺服系统机理模型,即机电伺服系统开环传递函数模型G(s);
步骤二、在机电伺服系统闭环的状态下,输入信号幅值为Ar、频率为f的正弦信号,使用相关辨识法获得该频率点上的实际控制信号幅值Au
步骤三、建立步骤二中输入信号幅值Ar与对应频率点上的实际控制信号幅值Au的线性关系,根据所需理想控制信号幅值计算补偿信号幅值ΔA;
步骤四、根据步骤三中得到的补偿信号幅值ΔA,将输入信号幅值从Ar变到实际信号幅值Ar′;
步骤五、令步骤二中的Ar=Ar′,重复步骤二到四,直到实际控制信号幅值Au与理想控制信号幅值相等,对此时的实际控制信号u和位置信号θ进行FFT数据处理,获得机电伺服系统模型在频率f上的幅相特性;
所述幅相特性包括幅值比和相角差;
步骤六、改变输入信号频率f,重复步骤二到步骤五,得到所需频带内各频率点的幅相特性;
步骤七、使用最小二乘法对步骤六中得到的各频率点的幅相特性数据进行拟合,得到步骤一中的机电伺服系统开环传递函数模型的参数τe、τm和K;
τe为电磁时间常数,τm为机电时间常数,K为标称模型等效增益。
2.根据权利要求1所述一种有限转角机电伺服系统的模型参数辨识方法,其特征在于:所述步骤一中根据机电伺服系统的结构建立机电伺服系统机理模型,即机电伺服系统开环传递函数模型G(s);具体过程为:
机电伺服系统机理模型为:
式中,id,iq为同步旋转d-q坐标系中的直轴电流和交轴电流,d为直轴,q为交轴;ωr为电机机械角速度;re为电机等效电阻;L为电机电枢电感;JΣ为轴系总转动惯量;D为电机旋转时的摩擦系数;Tl为干扰力矩,包含摩擦力矩及波动力矩;Pm为电机极对数;ψf为转子永磁效应对应的每对磁极磁通;ud,uq为同步旋转d-q坐标系中的直轴电压和交轴电压;为微分算子;
忽略机电伺服系统机理模型中id的影响,通过控制iq来直接控制电机的力矩输出,即
根据式(2),将电机电流反馈放大系数kf设为1,并忽略电机旋转时的摩擦系数D,令反电势系数电机力矩系数机电伺服系统机理模型变为
则由实际控制信号u到位置输出θ的机电伺服系统开环传递函数模型为:
式中,G(s)为机电伺服系统开环传递函数模型,s为拉普拉斯算子,kps为电机驱动器放大系数,kV为电机电流放大系数;记为电磁时间常数,JΣ为轴系总转动惯量,为机电时间常数,为标称模型等效增益,式(3)简化为:
对于机电伺服系统,满足τm>>τe,所以上式(4)又简化为:
至此,机电伺服系统开环传递函数模型G(s)建立完毕。
3.根据权利要求2所述一种有限转角机电伺服系统的模型参数辨识方法,其特征在于:所述步骤二中在机电伺服系统闭环的状态下,输入信号幅值为Ar、频率为f的正弦信号,使用相关辨识法获得该频率点上的实际控制信号幅值Au;具体过程为:
步骤二一、根据输入信号的频率f和采样时间Ts,选取及合适的F0,得到两个基准信号分别为fs01(k)=F0sin(2πf·kTs)和fs02(k)=F0cos(2πf·kTs);
式中,Ts为采样时间,F0为基准信号的幅值;k为取样点数;NR为取样点数;
步骤二二、当k=0,1,2,3,...NR-1时,根据计算相关函数R1(0)和R2(0);
其中fs(k)=u(t),即执行计算时的实际控制信号u;t为时间;fs(k)为基波信号;
步骤二三、计算实际控制信号幅值
4.根据权利要求3所述一种有限转角机电伺服系统的模型参数辨识方法,其特征在于:所述步骤三中建立步骤二中输入信号幅值Ar与实际控制信号幅值Au的线性关系,根据所需理想控制信号幅值计算补偿信号幅值ΔA;具体过程为:
根据步骤二得到的实际控制信号幅值Au,由于闭环控制器是线性的,得到输入信号幅值Ar与实际控制信号幅值Au的对应关系,根据所需理想控制信号幅值计算补偿信号幅值ΔA为:
5.根据权利要求4所述一种有限转角机电伺服系统的模型参数辨识方法,其特征在于:所述理想控制信号幅值为3V~5V电压信号。
6.根据权利要求5所述一种有限转角机电伺服系统的模型参数辨识方法,其特征在于:所述步骤四中根据步骤三中得到的补偿信号幅值ΔA,将输入信号幅值从Ar变到实际信号幅值Ar′;具体过程为:
将输入信号的幅值在1秒内按下式从输入信号幅值Ar增加到实际信号幅值Ar′;
其中rA(t)为补偿的输入信号,tsys为系统运行时间,tA为补偿开始的时间,tp为渐变时间。
7.根据权利要求6所述一种有限转角机电伺服系统的模型参数辨识方法,其特征在于:所述步骤五中令步骤二中的Ar=Ar′,重复步骤二到四,直到实际控制信号幅值Au与理想控制信号幅值相等,对此时的实际控制信号u和位置信号θ进行FFT数据处理,获得机电伺服系统模型在频率点f上的幅相特性;具体过程为:
步骤五一、对实际控制信号u进行FFT数据处理,得到其在频率点f上的幅值Au′和相角
步骤五二、对位置信号θ进行FFT数据处理,得到其在频率点f上的幅值Aθ′和相角
步骤五三、获得机电伺服系统模型在频率点f上的幅值比为Af=Aθ′/Au′、相角差为
CN201810108306.7A 2018-02-02 2018-02-02 一种有限转角机电伺服系统的模型参数辨识方法 Active CN108279571B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810108306.7A CN108279571B (zh) 2018-02-02 2018-02-02 一种有限转角机电伺服系统的模型参数辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810108306.7A CN108279571B (zh) 2018-02-02 2018-02-02 一种有限转角机电伺服系统的模型参数辨识方法

Publications (2)

Publication Number Publication Date
CN108279571A true CN108279571A (zh) 2018-07-13
CN108279571B CN108279571B (zh) 2020-09-01

Family

ID=62807462

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810108306.7A Active CN108279571B (zh) 2018-02-02 2018-02-02 一种有限转角机电伺服系统的模型参数辨识方法

Country Status (1)

Country Link
CN (1) CN108279571B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109115533A (zh) * 2018-09-07 2019-01-01 天津福云天翼科技有限公司 一种正弦波实时精确复现方法
CN109491251A (zh) * 2018-12-18 2019-03-19 华中科技大学 考虑数据扰动补偿的交流伺服系统模型辨识方法及设备
CN113037168A (zh) * 2019-12-25 2021-06-25 大族激光科技产业集团股份有限公司 一种永磁同步电机辨识电感电阻参数的方法
CN113110021A (zh) * 2021-03-17 2021-07-13 华南理工大学 一种用于伺服系统辨识与控制器设计的方法
CN113189868A (zh) * 2021-03-26 2021-07-30 哈尔滨工大航博科技有限公司 伺服系统动态误差的精确补偿方法
CN113325709A (zh) * 2021-05-26 2021-08-31 哈尔滨工大航博科技有限公司 转台变负载情况下的自适应控制方法
CN113406495A (zh) * 2021-06-28 2021-09-17 歌尔股份有限公司 振动电机的扫频特性曲线生成方法、装置及存储介质
CN114518709A (zh) * 2022-01-26 2022-05-20 哈尔滨工业大学 一种全姿态四轴转台框架角指令解算方法、设备和介质
CN115686090A (zh) * 2022-12-28 2023-02-03 武汉至驱动力科技有限责任公司 有限转角力矩电机的转角位置控制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080123811A1 (en) * 2006-11-08 2008-05-29 Curtis Steven E System and Method for Improved Collision Detection in an Imaging Device
CN104615130A (zh) * 2015-03-05 2015-05-13 哈尔滨工业大学 基于时间采样的机电伺服系统位置域周期性干扰的抑制方法
CN104101465B (zh) * 2014-08-05 2016-08-17 上海交通大学 转子硬支承动不平衡测试的电路频率特性误差补偿方法
CN106227040A (zh) * 2016-08-26 2016-12-14 哈尔滨工业大学 基于极值控制的伺服系统输入补偿方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080123811A1 (en) * 2006-11-08 2008-05-29 Curtis Steven E System and Method for Improved Collision Detection in an Imaging Device
CN104101465B (zh) * 2014-08-05 2016-08-17 上海交通大学 转子硬支承动不平衡测试的电路频率特性误差补偿方法
CN104615130A (zh) * 2015-03-05 2015-05-13 哈尔滨工业大学 基于时间采样的机电伺服系统位置域周期性干扰的抑制方法
CN106227040A (zh) * 2016-08-26 2016-12-14 哈尔滨工业大学 基于极值控制的伺服系统输入补偿方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
JINSHU WU等: "Dual EKF-Based Model Parameters Identification for Flight Simulation Turntable", 《PROCEEDINGS OF THE 36TH CHINESE CONTROL CONFERENCE》 *
陈宇飞: "机电伺服系统鲁棒控制设计与实现", 《中国优秀硕士学位论文全文数据库(电子期刊) 信息科技辑》 *

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109115533B (zh) * 2018-09-07 2020-10-27 天津福云天翼科技有限公司 一种正弦波实时精确复现方法
CN109115533A (zh) * 2018-09-07 2019-01-01 天津福云天翼科技有限公司 一种正弦波实时精确复现方法
CN109491251A (zh) * 2018-12-18 2019-03-19 华中科技大学 考虑数据扰动补偿的交流伺服系统模型辨识方法及设备
CN109491251B (zh) * 2018-12-18 2020-07-10 华中科技大学 考虑数据扰动补偿的交流伺服系统模型辨识方法及设备
CN113037168A (zh) * 2019-12-25 2021-06-25 大族激光科技产业集团股份有限公司 一种永磁同步电机辨识电感电阻参数的方法
CN113037168B (zh) * 2019-12-25 2023-05-09 大族激光科技产业集团股份有限公司 一种永磁同步电机辨识电感电阻参数的方法
CN113110021A (zh) * 2021-03-17 2021-07-13 华南理工大学 一种用于伺服系统辨识与控制器设计的方法
CN113189868B (zh) * 2021-03-26 2022-07-26 哈尔滨工大航博科技有限公司 伺服系统动态误差的精确补偿方法
CN113189868A (zh) * 2021-03-26 2021-07-30 哈尔滨工大航博科技有限公司 伺服系统动态误差的精确补偿方法
CN113325709B (zh) * 2021-05-26 2022-07-26 哈尔滨工大航博科技有限公司 转台变负载情况下的自适应控制方法
CN113325709A (zh) * 2021-05-26 2021-08-31 哈尔滨工大航博科技有限公司 转台变负载情况下的自适应控制方法
CN113406495B (zh) * 2021-06-28 2022-06-21 歌尔股份有限公司 振动电机的扫频特性曲线生成方法、装置及存储介质
CN113406495A (zh) * 2021-06-28 2021-09-17 歌尔股份有限公司 振动电机的扫频特性曲线生成方法、装置及存储介质
CN114518709A (zh) * 2022-01-26 2022-05-20 哈尔滨工业大学 一种全姿态四轴转台框架角指令解算方法、设备和介质
CN114518709B (zh) * 2022-01-26 2024-05-17 哈尔滨工业大学 一种全姿态四轴转台框架角指令解算方法、设备和介质
CN115686090A (zh) * 2022-12-28 2023-02-03 武汉至驱动力科技有限责任公司 有限转角力矩电机的转角位置控制方法

Also Published As

Publication number Publication date
CN108279571B (zh) 2020-09-01

Similar Documents

Publication Publication Date Title
CN108279571A (zh) 一种有限转角机电伺服系统的模型参数辨识方法
CN107065551A (zh) 一种基于模型参数精确辨识的仿真转台自校正控制方法
Mao et al. Design and implementation of continuous finite-time sliding mode control for 2-DOF inertially stabilized platform subject to multiple disturbances
CN103823379B (zh) 基于迭代学习的高频率角振动转台滑模控制方法
CN108716471A (zh) 一种磁悬浮分子泵转子极小位移主动控制方法
CN107133387B (zh) 转子不平衡系数变步长多边形迭代搜寻的不平衡补偿方法
CN107193211A (zh) 基于自抗扰和反演技术的单臂机械手控制器及其设计方法
CN105763119A (zh) 一种cmg框架永磁同步电机的控制系统及控制方法
CN107992063B (zh) 基于变参数章动阻尼的变速倾侧动量轮进动控制方法
CN111458531B (zh) 一种磁悬浮主轴基于转子位移的转速监测系统
CN103780188A (zh) 基于动态摩擦补偿的永磁球形电机转子自适应控制系统
Deng et al. Discrete-time direct model reference adaptive control application in a high-precision inertially stabilized platform
Zhou et al. A high precision compound control scheme based on non-singular terminal sliding mode and extended state observer for an aerial inertially stabilized platform
CN104009696A (zh) 一种基于滑模控制的交互式模型参考自适应速度与定子电阻的辨识方法
CN106100469A (zh) 基于自适应的电机伺服系统鲁棒位置控制器的实现方法
Zhou et al. A High‐Precision Control Scheme Based on Active Disturbance Rejection Control for a Three‐Axis Inertially Stabilized Platform for Aerial Remote Sensing Applications
Li et al. Dual drive system modeling and analysis for synchronous control of an H-type gantry
CN106208859B (zh) 基于干扰观测器与重复控制器的永磁同步电机控制方法
CN107203130B (zh) 基于极值搜索控制的仿真转台模型参数辨识方法
CN109143861B (zh) 一种基于力矩器的变速倾侧动量轮系统主动振动抑制方法
Mao et al. Output feedback stabilization of inertial stabilized platform with unmatched disturbances using sliding mode approach
CN106208844B (zh) 一种全局稳定的电机伺服系统输出反馈鲁棒控制方法
Yasien et al. Design of Fuzzy-Like Position Controller for Permanent Magnet Stepper Motor
CN116301081B (zh) 一种惯性测试设备的速率控制方法、装置、设备及介质
Yue et al. An observer-based robust adaptive controller for permanent magnet synchronous motor drive with initial rotor angle uncertainty

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