发明内容
本发明的目的在于提供一种伺服系统电流环控制参数自整定方法及伺服系统,旨在在伺服系统中未知参量对电流环动态特性存在影响的情况下,实现电流环控制参数的自整定。
本发明所述伺服系统电流环控制参数自整定方法包括以下步骤:
电机转子堵转,输入q轴电流周期方波指令信号;在t=0时刻,对所有粒子初始化,在允许取值范围内随机设置粒子的初始化位置x,将第i个粒子的自身个体极值设置成当前位置,全局极值设置成粒子群中的最优粒子位置;
随机给定电流环初始控制参数P=x,通过ADC采样,坐标变换后得到q轴电流跟踪响应信号,更新粒子位置,计算粒子i的适应度;
如果粒子i的适应度优于自身个体极值的适应度;如果当前进化代数中,粒子i的适应度优于全局极值的适应度;根据公式计算群体适应度方差;
判断算法是否满足收敛条件,如果满足就执行根据公式计算群体适应度方差,否则就对全局最优解按照公式执行变异操作并转回对所有粒子初始化;
求出全局最优解的目标函数值,并输出全局最优解,算法结束;通过伺服系统校验最优值等于全局极值,如果满足响应要求则整定成功,否则继续整定;相同整定结构,在确定伺服系统最优P值之后,整定系统I、D值;
最终校验伺服系统整体电流闭环响应特性。
进一步,所述伺服系统电流环控制参数自整定方法进一步包括以下步骤:
电机转子堵转,消除d轴电流反向电动势的影响;
在方波信号作用下,分析q轴一个周期T电流环PID参数;
分析时将方波信号一个周期分为高电平区间[0,T/2]、低电平区间[T/2,T];
电流方波信号函数用e(t)表示,高电平区间电流响应函数用e1(t)表示,低电平区间电流响应函数用e2(t)表示;
ITAE整定准则表达式为t表示时间,|e(t)|表示实际输出与期望输出的偏差值绝对值,ITAE准则控制系统瞬态响应振荡性小,对系统参数具有良好的选择性;对于伺服系统,通过ADC采样得到反馈相电流,然后进行坐标变换得到电流环跟踪响应电流;
对P值进行整定,初值P(0)对应ITAE指标为E(0);P(i)对应ITAE指标为E(i);i∈[1,∞),i∈n;
按照粒子群优化算法对P值进行动态赋值,变量P(i)值所对应的适应度函数用fi表示,当fi<2%时,此时得到最优伺服整定P(i)值,粒子群优化算法公式如下:
x(t+1)=wx(t)+c1r1(pbest-x(t))+c2r2(gbest-x(t));
w=(wmax-wmin)×exp(-β(t/Tmax)2)+wmin;
式中w为惯性权重,初始值取0.8,c1、c2为常数2,r1、r2为分布于[0,1]范围内的随机数,pbest为粒子本身找到的最优解,全局极值gbest为整个粒子群当前最优解;式中β取值由经验决定,为β∈[15,20];
根据群体适应度方差δ2判别局部极值是否是全局极值,群体适应度方差定义为下式:
式中n为粒子数,fi为第i个粒子适应度,favg为粒子群目前平均适应度,f为归一化定标因子,f的取值为下式:
f=max{1,max|f1-favg|},i∈[1,n];
如果出现粒子群过早收敛,则执行变异操作:
gbest=gbest×(1+τ×0.5)
τ为服从标准正态分布的随机变量,对gbest执行随机变异操作用来提高离子群算法跳出局部最优解的能力;
在确定最优伺服系统控制参数P值后,分别使D值取0,整定I值,I值取0,整定D值;
对得到的整定参数进行校验,若作用下的电流闭环阶跃响应满足快速、稳态误差小等特征,则认为参数整定结果满足电流环控制整定要求,整定过程结束,否则重新进行整定。
进一步,所述ITAE准则表达为:
进一步,所述惯性权重调整策略为:进化前期使用较大惯性权重,保证全局性;进化后期使用较小惯性权重,保证局部最优性。
一种利用上述的伺服系统电流环控制参数自整定方法的伺服控制系统,所述伺服控制系统包括:
上位机,与伺服驱动器通讯,用于发送或接收系统参数,发送速度信号;
伺服驱动器,根据速度值调节伺服系统的转速,将伺服系统实际速度值反馈给上位机;
上位机,用于通过人机界面输出实际速度值;
执行机构,用于按照伺服驱动器指令将电能转化为机械能。
本发明提供的交流伺服系统在许多高科技领域得到了非常广泛的应用,如激光加工、机器人、数控机床、大规模集成电路制造、办公自动化设备、雷达和各种军用武器随动系统、以及柔性制造系统等。交流伺服系统一般由机械执行机构、交流伺服驱动器、交流伺服系统和位置反馈检测装置构成。交流伺服驱动器是交流伺服系统的控制核心,其内部一般采用固定参数的PID控制结构来完成系统的调节过程。虽然PID控制结构具有算法简单、鲁棒性强和可靠性高等优点,但其控制性能与设置的控制参数密切相关,当控制参数设置不恰当时,无法使PID控制结构获得满意的控制性能。由于实际的伺服系统中存在非线性摩擦、放大器饱和、执行机构的死区与饱和等非线性特征,以及存在负载对象的不确定性,因此,在交流伺服系统的应用过程中,常常无法获得进行控制参数整定的准确依据,以致难以得到与实际应用过程相匹配的最优控制参数。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
下面结合附图对本发明的应用原理作详细的描述。如图1所示电流环参数自整定试验平台,固定电机转子后系统输入方波信号,结合伺服系统理论分析得到的电流环控制参数,作为粒子群算法初始参数值,在寻优规则下找寻适应度fi<2%时所对应的P值作为电流环自整定最优参数。校验寻优规则下的最优P控制参数,以实际的跟踪曲线为标准,观察响应电流闭环阶跃上升时间、超调量、调节时间。如果校验结果满意,则可认为自整定成功,将参数应用于控制器;如果校验结果不满意,可将粒子群算法复位,直至找寻到校验满意结果。
本发明提供的伺服系统电流环控制参数自整定方法及伺服系统,在实验平台上予以验证。图1为电流环控制参数自整定试验平台,锁死电机转子,输入周期方波指令信号,由粒子群算法随机给定P值,通过上位机采样电流环输出Iq跟踪响应曲线。图2为伺服系统在P1、P2两种电流环控制参数P作用下的q轴电流响应曲线,由图可见,P1上升时间小于P2,P2超调量大于P1,P1调节时间小于P2,使用ITAE公式运算结果亦可知E(p1)<E(p2),综上可见P1参数值优于P2。通过ITAE公式作为标准,结合粒子群算法,最终选定当fi<2%时,此时得到最优伺服整定P(i)值。
如图3所示,本发明实施例的伺服系统电流环控制参数自整定方法及伺服系统包括以下步骤:
S101:电机转子堵转,输入q轴电流周期方波指令信号;
S102:t=0时刻,对所有粒子初始化,在允许取值范围内随机设置粒子的初始化位置x,将第i个粒子的自身个体极值设置成当前位置,全局极值设置成粒子群中的最优粒子位置;
S103:随机给定电流环初始控制参数P=x,通过ADC采样,坐标变换后得到q轴电流跟踪响应信号,更新粒子位置,计算粒子i的适应度;
S104:如果粒子i的适应度优于自身个体极值的适应度;如果当前进化代数中,粒子i的适应度优于全局极值的适应度;根据公式计算群体适应度方差;
S105:判断算法是否满足收敛条件,如果满足就执行根据公式计算群体适应度方差,否则就对全局最优解按照公式执行变异操作并转回对所有粒子初始化;
S106:求出全局最优解的目标函数值,并输出全局最优解,算法结束;通过伺服系统校验最优值等于全局极值,如果满足响应要求则整定成功,否则继续整定;相同整定结构,在确定伺服系统最优P值之后,整定系统I、D值;
S107:最终校验伺服系统整体电流闭环响应特性。
本发明的具体步骤如下:
Step1:电机转子堵转,输入周期方波指令信号。
Step2:t=0时刻,对所有粒子初始化,在允许取值范围内随机设置粒子的初始化位置x,将第i个粒子的pbest设置成当前位置,gbest设置成粒子群中的最优粒子位置,i∈[1,n],n为粒子数。
Step3:随机给定电流环初始P=x,通过ADC采样,得到离散反馈跟踪响应信号。
Step4:按公式(2)、(3)更新粒子位置;
Step5:计算粒子i的适应度fi;
Step6:如果粒子i的适应度fi优于自身个体极值pbest的适应度f(Pbest),就用粒子当前的位置xi替换pbest;
Step7:如果当前进化代数中,粒子i的适应度fi优于全局极值gbest的适应度f(gbest),就用粒子当前位置xi替换全局极值gbest;
Step8:根据公式(4)、(5)计算群体适应度方差δ2;
Step9:判断算法是否满足收敛条件,如果满足就执行Step8,否则就对gbest按照公式(6)执行变异操作并转回Step2;
Step10:求出全局最优解gbest的目标函数值,并输出全局最优解gbest,算法结束。
Step11:通过伺服系统校验最优值P=gbest,如果满足响应要求则整定成功,否则继续整定。
Step12:相同整定结构,在确定伺服系统最优P值之后,整定系统I、D值。
Step13:最终校验伺服系统整体电流闭环响应特性。
下面结合附图对本发明的应用原理作进一步的描述。
本发明伺服系统电流环整定选用q轴和d轴电流控制参数作为整定分析对象。在电流周期方波信号作用下,通过粒子群优化算法调整PID控制参数,将电流阶跃响应ITAE指标作为粒子群算法适应函数,最终根据所需ITAE指标选取对应最优PID整定参数值。分析计算q轴和d轴电流动态响应特性,电流环q轴和d轴电流控制结构相同,选取q轴电流作为分析对象,将整定控制参数应用于d轴电流控制过程。
1.电机转子堵转,消除d轴电流反向电动势的影响。
2.在电流周期方波指令信号作用下,分析q轴一个周期T电流环PID参数。
3.分析时将电流周期方波指令信号周期分为高电平区间[0,T/2]、低电平区间[T/2,T]。
4.电流误差函数用e(t)表示,高电平区间电流响应函数用e1(t)表示,低电平区间电流响应函数用e2(t)表示。
5.ITAE整定准则表达式为t表示时间,|e(t)|表示实际输出与期望输出的偏差值绝对值。ITAE准则控制系统瞬态响应振荡性小,对系统参数具有良好的选择性。对于伺服系统,电流环通过ADC采样得到跟踪响应电流,通过坐标变换和公式计算可得到电流反馈量,则ITAE准则可表达为:
6.为保证动态响应特性,首先对P值进行整定。设初值P(0)对应ITAE指标为E(0);P(i)对应ITAE指标为E(i);i∈[1,∞),i∈n。
7.按照粒子群优化算法对P值进行动态赋值,变量P(i)值所对应的适应度函数用fi表示。当fi<2%时,此时得到最优伺服整定P(i)值。粒子群优化算法公式如下:
x(t+1)=wx(t)+c1r1(pbest-x(t))+c2r2(gbest-x(t))(2)
w=(wmax-wmin)×exp(-β(t/Tmax)2)+wmin(3)
式(2)中w为惯性权重,初始值取0.8,c1、c2为常数2,r1、r2为分布于[0,1]范围内的随机数,pbest为粒子本身找到的最优解,全局极值gbest为整个粒子群当前最优解。
式(3)中β取值由经验决定,一般为β∈[15,20]。惯性权重调整策略为:进化前期使用较大惯性权重,保证全局性;进化后期使用较小惯性权重,保证局部最优性。
8.为防止粒子群算法早熟收敛,可根据群体适应度方差δ2判别局部极值是否是全局极值。群体适应度方差可定义为下式:
式中n为粒子数,fi为第i个粒子适应度,favg为粒子群目前平均适应度,f为归一化定标因子,f的取值为下式:
f=max{1,max|f1-favg|},i∈[1,n](5)
如果出现粒子群过早收敛,则执行变异操作:
gbest=gbest×(1+τ×0.5)(6)
τ为服从标准正态分布的随机变量,对gbest执行随机变异操作用来提高离子群算法跳出局部最优解的能力。
9.在确定最优伺服系统P值后,分别整定I值(D值取0)和D值(I值取0)。伺服系统电流环自整定结构图如图4所示:
10.对得到的整定参数进行校验,若其作用下的电流闭环阶跃响应满足快速、稳态误差小等特征,则认为参数整定结果满足电流环控制整定要求,整定过程结束,否则重新进行整定。
图2为伺服控制系统结构图,根据此图对伺服做简要介绍:
上位机与伺服驱动器通讯,发送或接收系统参数,如发送速度信号,伺服驱动器根据速度值调节伺服系统的转速;如伺服驱动器将伺服系统实际速度值反馈给上位机,上位机通过人机界面呈现给我们。伺服系统作为能量转换机构,为执行机构提供机械能,按照伺服驱动器指令将电能绝大部分转化为机械能,辅助执行机构完成指定任务。执行机构由伺服系统提供的机械能,完成设定的速度、位置等指定。由最简伺服系统结构图可见,上位机通过人机界面与我们进行交流,我们将需要执行机构完成的指令通过上位机下达给伺服驱动器,驱动器调节伺服系统完成电能转化为机械能,伺服系统将机械能通过机械结构传递给执行机构最终实现我们下达指令。伺服驱动器在执行机构完成任务后,将系统运行情况反馈给上位机,上位机通过人机界面反馈给我们,完成伺服系统任务。电流环作为速度环和位置环的基础,电流环控制参数直接影响系统的整体性能。为此可见,伺服系统电流环控制参数的选择显得尤为重要。本发明专利提出一种伺服系统电流环控制参数的自整定方法,能有效提高伺服系统控制精度。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。