CN107800343B - 异步电机自抗扰控制器的设计方法 - Google Patents
异步电机自抗扰控制器的设计方法 Download PDFInfo
- Publication number
- CN107800343B CN107800343B CN201711187196.XA CN201711187196A CN107800343B CN 107800343 B CN107800343 B CN 107800343B CN 201711187196 A CN201711187196 A CN 201711187196A CN 107800343 B CN107800343 B CN 107800343B
- Authority
- CN
- China
- Prior art keywords
- shaft current
- disturbance rejection
- nonlinear
- rejection controller
- single order
- 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.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/0003—Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
- H02P21/0017—Model reference adaptation, e.g. MRAS or MRAC, useful for control or parameter estimation
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/24—Vector control not involving the use of rotor position or rotor speed sensors
- H02P21/26—Rotor flux based control
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
本发明公开了一种异步电机自抗扰控制器的设计方法,包括步骤:一、构建异步电机在按转子磁链定向的同步旋转正交坐标系下的动态数学模型;二、确定出由转速环一阶自抗扰控制器、磁链环一阶自抗扰控制器、d轴电流环一阶自抗扰控制器和q轴电流环一阶自抗扰控制器四个一阶自抗扰控制器构成的异步电机矢量控制系统的结构;三、分别设计步骤二中所述异步电机矢量控制系统的转速环一阶自抗扰控制器、磁链环一阶自抗扰控制器、d轴电流环一阶自抗扰控制器和q轴电流环一阶自抗扰控制器。本发明方法步骤简单,增强了系统的抗负载扰动和抗电机参数变化的能力,提高了系统的控制精度及鲁棒性,实用性强。
Description
技术领域
本发明属于异步电机控制技术领域,具体涉及一种异步电机自抗扰控制器的设计方法。
背景技术
异步电机具有结构简单、坚固耐用、运行可靠等优点,在工业应用场合得到了非常广泛的应用。矢量控制方式的引入,实现了转矩与磁链的解耦控制,使异步电机变频调速性能达到与直流电机控制性能相媲美的程度。然而,在实际应用过程中,传统的PI矢量控制系统受负载变化及电机参数变化的影响较大,因此,控制性能受到很大影响。滑模变结构控制、自适应控制、预测控制、内模控制等现代控制策略在交流电机变频控制研究领域受到了广泛关注。但是,滑模变结构控制存在“抖振”现象,对于矢量控制系统在低速时会引起较大转矩脉动;自适应控制算法复杂,对处理器性能要求较高;预测控制算法存在模型预测精度不高、滚动优化策略少、反馈校正方法单调等问题;内模控制器的固定时间常数在鲁棒性与快速性之间存在矛盾。
自抗扰控制器是一种新型的非线性控制器,通过对系统内、外扰动的实时估计并补偿,结合非线性控制策略,从而系统获得更好的鲁棒性和动态性能。
然而,就目前的异步电机自抗扰控制器在设计过程中存在没有完全按照自抗扰控制器数学特征来设计,控制器设计中存在通过线性化的方法进行简化,并且仅凭经验对控制器参数的进行取值,使得异步电机变频调速系统的达不到理想的控制性能。此外,对异步电机在同步旋转坐标系下的电流id,iq、转速ωr、磁链ψr在数学模型上的不同特征没有对控制器的设计做出区分,很大程度上影响了控制器的控制性能。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足,提供一种异步电机自抗扰控制器的设计方法,其方法步骤简单,增强了系统的抗负载扰动和抗电机参数变化的能力,提高了系统的控制精度及鲁棒性,实用性强。
为解决上述技术问题,本发明采用的技术方案是:一种异步电机自抗扰控制器的设计方法,其特征在于,该方法包括以下步骤:
步骤一、构建异步电机在按转子磁链定向的同步旋转正交坐标系下的动态数学模型,所述动态数学模型包括转速环动态数学模型、q轴电流环动态数学模型、磁链环动态数学模型和d轴电流环动态数学模型;
步骤二、根据步骤一中构建的动态数学模型,并结合自抗扰控制原理,确定出由转速环一阶自抗扰控制器、磁链环一阶自抗扰控制器、d轴电流环一阶自抗扰控制器和q轴电流环一阶自抗扰控制器四个一阶自抗扰控制器构成的异步电机矢量控制系统的结构;
所述异步电机矢量控制系统包括转速环一阶自抗扰控制器、磁链环一阶自抗扰控制器、d轴电流环一阶自抗扰控制器和q轴电流环一阶自抗扰控制器,以及PARK变换模块、PARK反变换模块、CLARK变换模块、SVPWM模块、逆变器和转子磁链观测器;所述q轴电流环一阶自抗扰控制器与转速环一阶自抗扰控制器的输出端连接,所述d轴电流环一阶自抗扰控制器与磁链环一阶自抗扰控制器的输出端连接,所述q轴电流环一阶自抗扰控制器的输出端和d轴电流环一阶自抗扰控制器的输出端均与PARK反变换模块的输入端连接,所述PARK反变换模块的输出端与SVPWM模块的输入端连接,所述SVPWM模块的输出端与逆变器连接,异步电机与逆变器的A相电流输出端、B相电流输出端和C相电流输出端均连接,所述逆变器的A相电流输出端和B相电流输出端均与CLARK变换模块的输入端连接,所述CLARK变换模块的输出端与PARK变换模块的输入端连接,所述PARK变换模块的异步电机定子侧q轴电流分量输出端与q轴电流环一阶自抗扰控制器连接,所述PARK变换模块的异步电机定子侧d轴电流分量输出端与d轴电流环一阶自抗扰控制器连接,所述PARK变换模块的异步电机定子侧q轴电流分量输出端和异步电机定子侧d轴电流分量输出端均与转子磁链观测器的输入端连接,所述转子磁链观测器的异步电机转子侧d轴磁链信号输出端与磁链环一阶自抗扰控制器连接,所述转子磁链观测器的旋转定向角输出端与PARK变换模块和PARK反变换模块均连接。
步骤三、根据一阶自抗扰控制器由跟踪微分器、扩张状态观测器和非线性状态误差反馈控制律三部分组成的原理,分别设计步骤二中所述异步电机矢量控制系统的转速环一阶自抗扰控制器、磁链环一阶自抗扰控制器、d轴电流环一阶自抗扰控制器和q轴电流环一阶自抗扰控制器。
上述的异步电机自抗扰控制器的设计方法,其特征在于:步骤一中所述转速环动态数学模型为步骤一中所述q轴电流环动态数学模型为步骤一中所述磁链环动态数学模型为步骤一中所述d轴电流环动态数学模型为其中,Tr为异步电机转子时间常数且Tr=Lr/Rr,σ为异步电机漏磁系数且usd为异步电机定子侧d轴电压分量,usq为异步电机定子侧q轴电压分量,isd为异步电机定子侧d轴电流分量,isq为异步电机定子侧q轴电流分量;Rs为异步电机定子侧电阻,Rr为异步电机转子侧电阻,Ls为异步电机定子侧电感,Lr为异步电机转子侧电感;Lm为异步电机定转子间互感;ψrd为异步电机转子侧d轴磁链;ω1为异步电机同步角转速;ωr为异步电机转子转速;np为异步电机极对数;J为异步电机转动惯量;TL为负载转矩。
上述的异步电机自抗扰控制器的设计方法,其特征在于:步骤三中设计所述转速环一阶自抗扰控制器的过程为:
步骤3011、将步骤一中所述转速环动态数学模型变形为其中,b为转速环补偿参数且w2(t)为转速环扰动且
步骤3012、将转速环跟踪微分器过渡过程的离散方程表示为:
其中,v(k)为转速环一阶自抗扰控制器的输入信号,x11(k)为跟踪输入信号v(k)的信号,x11(k+1)为x11(k)的下一时刻的变量,x12(k)为x11(k)的微分,x12(k+1)为x12(k)的下一时刻的变量,T为离散控制周期,h为转速环跟踪微分器的滤波参数,r为决定跟踪输入信号v(k)的速度的参数,fst[x11(k)-v(k),x12(k),r,h]为转速环离散时间系统最优控制函数;
步骤3013、将w2(t)视为扰动,构造转速环扩张状态观测器为:
其中,z11(k)为对ωr的估计,z11(k+1)为z11(k)的下一时刻的变量,y(k)为ωr的反馈值,e10(k)为z11(k)与y(k)的差值,z12(k)为对转速环扩张状态观测器系统的总扰动的估计,z12(k+1)为z12(k)的下一时刻的变量,fal(e10(k),α1,δ1)为转速环扩张状态观测器的非线性组合函数且α1为转速环扩张状态观测器的非线性参数,δ1为转速环扩张状态观测器的滤波参数,β01为扩张状态观测器转速估计值的输出误差校正增益,β02为转速环扩张状态观测器系统扰动估计值的输出误差校正增益,u(k)为转速环一阶自抗扰控制器的输出控制量;
步骤3014、根据跟踪微分器和扩张状态观测器的输出,构造非线性状态反馈控制律为:
其中,u(k)为转速环一阶自抗扰控制器的输出控制量,e11(k)为x11(k)与z11(k)的差值,u10(k)为转速环非线性状态反馈控制律的输出,β1为转速环非线性状态反馈控制律输出误差校正增益,fal(e11(k),α01,δ01)为转速环非线性状态反馈控制律的非线性组合函数且α01为转速环非线性状态反馈控制律的非线性参数,δ01为转速环非线性状态反馈控制律的滤波参数。
上述的异步电机自抗扰控制器的设计方法,其特征在于:步骤3012中根据转速环跟踪微分器的输出波形并根据转速环离散时间系统最优控制函数fst[x11(k)-v(k),x12(k),r,h]对h和r进行调整,其中,转速环离散时间系统最优控制函数fst[x11(k)-v(k),x12(k),r,h]用公式表示为:
其中,a为中间变量且
步骤3013中所述α1的取值范围为0<α1<1,所述δ1的取值为0.01,β01=1/T,β02=1/(5T2);
步骤3014中所述α01的取值范围为0<α01<1,所述δ01的取值为0.01,所述β1的取值为10/b。
上述的异步电机自抗扰控制器的设计方法,其特征在于:步骤三中设计所述q轴电流环一阶自抗扰控制器的过程为:
步骤3021、将步骤一中所述q轴电流环动态数学模型变形为其中,b1为q轴电流环补偿参数且b1=1/σ,w3(t)为q轴电流环扰动且w3(t)=-Lmψrdωr/(σLr)-ω1isd;
步骤3022、将q轴电流环跟踪微分器过渡过程的离散方程表示为:
其中,v1(k)为q轴电流环一阶自抗扰控制器的输入信号,x21(k)为跟踪输入信号v1(k)的信号,x21(k+1)为x21(k)的下一时刻的变量,x22(k)为x21(k)的微分,x22(k+1)为x22(k)的下一时刻的变量,T为离散控制周期,h1为q轴电流环跟踪微分器的滤波参数,r1为决定跟踪输入信号v1(k)的速度的参数,fst[x21(k)-v1(k),x22(k),1r,h1]为q轴电流环离散时间系统最优控制函数;
步骤3023、将视为被控对象的已知部分,构造扩张状态观测器为:
其中,z21(k)为对isq的估计,z21(k+1)为z21(k)的下一时刻的变量,y1(k)为isq的反馈值,e20(k)为z21(k)与y1(k)的差值,z22(k)为对q轴电流环扩张状态观测器系统的总扰动的估计,z22(k+1)为z22(k)的下一时刻的变量,fal(e20(k),α2,δ2)为q轴电流环扩张状态观测器系统的非线性组合函数且α2为q轴电流环扩张状态观测器系统的非线性参数,δ2为q轴电流环扩张状态观测器的滤波参数,β11为扩张状态观测器q轴电流估计值的输出误差校正增益、β12为q轴电流环扩张状态观测器系统扰动估计值的输出误差校正增益,u1(k)为q轴电流环一阶自抗扰控制器的输出控制量;
步骤3024、根据跟踪微分器和扩张状态观测器的输出,构造非线性状态反馈控制律为:
其中,u1(k)为q轴电流环一阶自抗扰控制器的输出控制量,e21(k)为x21(k)与z21(k)的差值,u20(k)为q轴电流环非线性状态反馈控制律的输出,β2为q轴电流环非线性状态反馈控制律输出误差校正增益,fal(e21(k),α11,δ11)为q轴电流环非线性状态反馈控制律的非线性组合函数且α11为q轴电流环非线性状态反馈控制律的非线性参数,δ11为q轴电流环非线性状态反馈控制律的滤波参数。
上述的异步电机自抗扰控制器的设计方法,其特征在于:步骤3022中根据q轴电流环跟踪微分器的输出波形并根据q轴电流环离散时间系统最优控制函数fst[x21(k)-v1(k),x22(k),r1,h1]对h1和r1进行调整,其中,q轴电流环离散时间系统最优控制函数fst[x21(k)-v1(k),x22(k),r1,h1]用公式表示为:
其中,a1为中间变量且
步骤3023中所述α2的取值范围为0<α2<1,所述δ2的取值为0.01,β11=1/T,β12=1/(5T2);
步骤3024中所述α11的取值范围为0<α11<1,所述δ11的取值为0.01,所述β2的取值为10/b1。
上述的异步电机自抗扰控制器的设计方法,其特征在于:步骤三中设计所述磁链环一阶自抗扰控制器的过程为:
步骤3031、将步骤一中所述磁链环动态数学模型变形为其中,b2为磁链环补偿参数且b2=Lm/Tr;
步骤3032、将磁链环跟踪微分器过渡过程的离散方程表示为:
其中,v2(k)为磁链环一阶自抗扰控制器的输入信号,x31(k)为跟踪输入信号v2(k)的信号,x31(k+1)为x31(k)的下一时刻的变量,x32(k)为x31(k)的微分,x32(k+1)为x32(k)的下一时刻的变量,T为离散控制周期,h2为磁链环跟踪微分器的滤波参数,r2为决定跟踪输入信号v2(k)速度的参数,fst[x31(k)-v2(k),x32(k),r2,h2]为磁链环离散时间系统最优控制函数;
步骤3033、将f1(z31(k))=-ψrd/Tr视为被控对象的已知部分,构造扩张状态观测器为:
其中,z31(k)视为对ψrd的估计,z31(k+1)为z31(k)的下一时刻的变量,y2(k)为ψrd的反馈值,e30(k)为z31(k)与y2(k)的差值,z32(k)为对磁链环扩张状态观测器系统的总扰动的估计,z32(k+1)为z32(k)的下一时刻的变量,fal(e30(k),α3,δ3)为磁链环扩张状态观测器系统的非线性组合函数且α3为磁链环扩张状态观测器的非线性参数,δ3为磁链环扩张状态观测器的滤波参数,β21为扩张状态观测器磁链估计值的输出误差校正增益,β22为磁链环扩张状态观测器系统扰动估计值的输出误差校正增益,u2(k)为磁链环一阶自抗扰控制器的输出控制量;
步骤3034、根据跟踪微分器和扩张状态观测器的输出,构造非线性状态反馈控制律为:
其中,u2(k)为磁链环一阶自抗扰控制器的输出控制量,e31(k)为x31(k)与z31(k)的差值,u30(k)为磁链环非线性状态反馈控制律的输出,β3为磁链环非线性状态反馈控制律输出误差校正增益,fal(e31(k),α21,δ21)为磁链环非线性状态反馈控制律的非线性组合函数且α21为磁链环非线性状态反馈控制律的非线性参数,δ21为磁链环非线性状态反馈控制律的滤波参数。
上述的异步电机自抗扰控制器的设计方法,其特征在于:步骤3032中根据磁链环跟踪微分器的输出波形并根据磁链环离散时间系统最优控制函数fst[x31(k)-v2(k),x32(k),r2,h2]对h2和r2进行调整,其中,磁链环离散时间系统最优控制函数fst[x31(k)-v2(k),x32(k),r2,h2]用公式表示为:
其中a2为中间变量且
步骤3033中所述α3的取值范围为0<α3<1,所述δ3的取值为0.01,β21=1/T,β22=1/(5T2);
步骤3034中所述α21的取值范围为0<α21<1,所述δ21的取值为0.01,所述β3的取值为10/b2。
上述的异步电机自抗扰控制器的设计方法,其特征在于:步骤三中设计所述d轴电流环一阶自抗扰控制器的过程为:
步骤3041、将步骤一中所述d轴电流环动态数学模型变形为其中,b3为d轴电流环补偿参数且b3=1/σ,w1(t)为d轴电流环扰动且
步骤3042、将d轴电流环跟踪微分器过渡过程的离散方程表示为:
其中,v3(k)为d轴电流环一阶自抗扰控制器的输入信号,x41(k)为跟踪输入信号v3(k)的信号,x41(k+1)为x41(k)的下一时刻的变量,x42(k)为x41(k)的微分,x42(k+1)为x42(k)的下一时刻的变量,T为离散控制周期,h3为d轴电流环跟踪微分器的滤波参数,r3为决定跟踪输入信号v3(k)的速度的参数,fst[x41(k)-v3(k),x42(k),r3,h3]为d轴电流环离散时间系统最优控制函数;
步骤3043、将视为被控对象的已知部分,构造扩张状态观测器为:
其中,z41(k)为对isd的估计,z41(k+1)为z41(k)的下一时刻的变量,y3(k)为isd的反馈值,e40(k)为z41(k)与y3(k)的差值,z42(k)为对d轴电流环扩张状态观测器系统的总扰动的估计,z42(k+1)为z42(k)的下一时刻的变量,fal(e40(k),α4,δ4)为d轴电流环扩张状态观测器系统的非线性组合函数且α4为d轴电流环扩张状态观测器的非线性参数,δ4为d轴电流环扩张状态观测器的滤波参数,β31为扩张状态观测器d轴电流估计值的输出误差校正增益、β32为d轴电流环扩张状态观测器系统扰动估计值的输出误差校正增益,u3(k)为d轴电流环一阶自抗扰控制器的输出控制量;
步骤3044、根据跟踪微分器和扩张状态观测器的输出,构造非线性状态反馈控制律为:
其中,u3(k)为d轴电流环一阶自抗扰控制器的输出控制量,e41(k)为x41(k)与z41(k)的差值,u40(k)为d轴电流环非线性状态反馈控制律的输出,β4为d轴电流环非线性状态反馈控制律输出误差校正增益,fal(e41(k),α31,δ31)为d轴电流环非线性状态反馈控制律的非线性组合函数且α31为d轴电流环非线性状态反馈控制律的非线性参数,δ31为d轴电流环非线性状态反馈控制律滤波参数。
上述的异步电机自抗扰控制器的设计方法,其特征在于:步骤3042中根据d轴电流环跟踪微分器的输出波形并根据d轴电流环离散时间系统最优控制函数fst对h3和r3进行调整,其中,d轴电流环离散时间系统最优控制函数fst用公式表示为:
其中,a3为中间变量且
步骤3043中所述α4的取值范围为0<α4<1,所述δ4的取值为0.01,β31=1/T,β32=1/(5T2);
步骤3044中所述α31的取值范围为0<α31<1,所述δ31的取值为0.01,所述β4的取值为10/b3。
本发明与现有技术相比具有以下优点:
1、本发明在异步电机自抗扰控制器设计过程中,对未知扰动和已知扰动进行区分,并对自抗扰控制器众多的参数提出了整定方法,根据自抗扰控制器的数学特征和异步电机的数学模型,设计了异步电机矢量控制系统的自抗扰转速环控制器、自抗扰磁链环控制器、自抗扰d轴电流环控制器和自抗扰q轴电流环控制器,对4个自抗扰控制器的设计步骤、设计方法作了详细描述,对跟踪微分器(TD)、扩张状态观测器(ESO)以及非线性状态反馈控制律(NLSEF)的各个参数对自抗扰控制器性能的影响作了详细描述,所设计的异步电机自抗扰控制器的鲁棒性和动态性能好。
2、本发明通过Matlab/Simulink仿真环境,对基于自抗扰控制器的矢量控制系统和基于PI控制器的矢量控制系统进行了对比研究,仿真结果表明自抗扰控制器(ADRC)系统的鲁棒性、动态性能和控制精度均优于PI系统。
3、本发明的方法步骤简单,增强了系统的抗负载扰动和抗电机参数变化的能力,提高了系统的控制精度及鲁棒性,实用性强,使用效果好,便于推广使用。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明的方法流程框图。
图2为本发明异步电机矢量控制系统的结构框图。
图3为本发明实施例中突加/减负载时ADRC系统与PI系统的转速波形对比图。
图4A为本发明实施例中PI控制系统转速波形图。
图4B为本发明实施例中ADRC控制系统转速波形图。
附图标记说明:
1—转速环一阶自抗扰控制器; 2—磁链环一阶自抗扰控制器;
3—d轴电流环一阶自抗扰控制器; 4—q轴电流环一阶自抗扰控制器;
5—PARK变换模块; 6—PARK反变换模块;
7—CLARK变换模块; 8—SVPWM模块;
9—逆变器; 10—转子磁链观测器;
11—异步电机。
具体实施方式
如图1所示,本发明的异步电机自抗扰控制器的设计方法,包括以下步骤:
步骤一、构建异步电机11在按转子磁链定向的同步旋转正交坐标系下的动态数学模型,所述动态数学模型包括转速环动态数学模型、q轴电流环动态数学模型、磁链环动态数学模型和d轴电流环动态数学模型;
步骤二、根据步骤一中构建的动态数学模型,并结合自抗扰控制原理,确定出由转速环一阶自抗扰控制器1、磁链环一阶自抗扰控制器2、d轴电流环一阶自抗扰控制器3和q轴电流环一阶自抗扰控制器4四个一阶自抗扰控制器构成的异步电机矢量控制系统的结构;
如图2所示,所述异步电机矢量控制系统包括转速环一阶自抗扰控制器1、磁链环一阶自抗扰控制器2、d轴电流环一阶自抗扰控制器3和q轴电流环一阶自抗扰控制器4,以及PARK变换模块5、PARK反变换模块6、CLARK变换模块7、SVPWM模块8、逆变器9和转子磁链观测器10;所述q轴电流环一阶自抗扰控制器4与转速环一阶自抗扰控制器1的输出端连接,所述d轴电流环一阶自抗扰控制器3与磁链环一阶自抗扰控制器2的输出端连接,所述q轴电流环一阶自抗扰控制器4的输出端和d轴电流环一阶自抗扰控制器3的输出端均与PARK反变换模块6的输入端连接,所述PARK反变换模块6的输出端与SVPWM模块8的输入端连接,所述SVPWM模块8的输出端与逆变器9连接,异步电机11与逆变器9的A相电流输出端、B相电流输出端和C相电流输出端均连接,所述逆变器9的A相电流输出端和B相电流输出端均与CLARK变换模块7的输入端连接,所述CLARK变换模块7的输出端与PARK变换模块5的输入端连接,所述PARK变换模块5的异步电机定子侧q轴电流分量输出端与q轴电流环一阶自抗扰控制器4连接,所述PARK变换模块5的异步电机定子侧d轴电流分量输出端与d轴电流环一阶自抗扰控制器3连接,所述PARK变换模块5的异步电机定子侧q轴电流分量输出端和异步电机定子侧d轴电流分量输出端均与转子磁链观测器10的输入端连接,所述转子磁链观测器10的异步电机转子侧d轴磁链信号输出端与磁链环一阶自抗扰控制器2连接,所述转子磁链观测器10的旋转定向角输出端与PARK变换模块5和PARK反变换模块6均连接。为定子电压的α轴分量给定;为定子电压的β轴分量给定;isa为A相定子电流;isb为B相定子电流;isα为定子电流的α轴分量;isβ为定子电流的β轴分量;θs为旋转定向角。
步骤三、根据一阶自抗扰控制器由跟踪微分器(TD)、扩张状态观测器(ESO)和非线性状态误差反馈控制律(NLSEF)三部分组成的原理,分别设计步骤二中所述异步电机矢量控制系统的转速环一阶自抗扰控制器1、磁链环一阶自抗扰控制器2、d轴电流环一阶自抗扰控制器3和q轴电流环一阶自抗扰控制器4。
本实施例中,步骤一中所述转速环动态数学模型为步骤一中所述q轴电流环动态数学模型为步骤一中所述磁链环动态数学模型为步骤一中所述d轴电流环动态数学模型为其中,Tr为异步电机转子时间常数且Tr=Lr/Rr,σ为异步电机漏磁系数且usd为异步电机定子侧d轴电压分量,usq为异步电机定子侧q轴电压分量,isd为异步电机定子侧d轴电流分量,isq为异步电机定子侧q轴电流分量;Rs为异步电机定子侧电阻,Rr为异步电机转子侧电阻,Ls为异步电机定子侧电感,Lr为异步电机转子侧电感;Lm为异步电机定转子间互感;ψrd为异步电机转子侧d轴磁链;ω1为异步电机同步角转速;ωr为异步电机转子转速;np为异步电机极对数;J为异步电机转动惯量;TL为负载转矩。
本实施例中,步骤三中设计所述转速环一阶自抗扰控制器1的过程为:
步骤3011、将步骤一中所述转速环动态数学模型变形为其中,b为转速环补偿参数且w2(t)为转速环扰动且
步骤3012、将转速环跟踪微分器(TD)过渡过程的离散方程表示为:
其中,v(k)为转速环一阶自抗扰控制器1的输入信号,对应图2转速环ADRC的转速给定x11(k)为跟踪输入信号v(k)的信号,x11(k+1)为x11(k)的下一时刻的变量,x12(k)为x11(k)的微分,x12(k+1)为x12(k)的下一时刻的变量,T为离散控制周期,h为转速环跟踪微分器(TD)的滤波参数,r为决定跟踪输入信号v(k)的速度的参数,fst[x11(k)-v(k),x12(k),r,h]为转速环离散时间系统最优控制函数;跟踪微分器(TD)中,h和r为可调参数,h越大,对噪声滤波效果越好;r越大跟踪速度越快,但同时也会增加超调量;
步骤3013、将w2(t)视为扰动,构造转速环扩张状态观测器(ESO)为:
其中,z11(k)为对ωr的估计,z11(k+1)为z11(k)的下一时刻的变量,y(k)为ωr的反馈值,对应图2转速环ADRC的转速反馈ωr,e10(k)为z11(k)与y(k)的差值,z12(k)为对转速环扩张状态观测器(ESO)系统的总扰动的估计,z12(k+1)为z12(k)的下一时刻的变量,fal(e10(k),α1,δ1)为转速环扩张状态观测器(ESO)的非线性组合函数且α1为转速环扩张状态观测器(ESO)的非线性参数,δ1为转速环扩张状态观测器(ESO)的滤波参数,β01为扩张状态观测器(ESO)转速估计值的输出误差校正增益,β02为转速环扩张状态观测器(ESO)系统扰动估计值的输出误差校正增益,u(k)为转速环一阶自抗扰控制器1的输出控制量;α1,δ1,β01,β02为可调参数,α1越小,fal函数非线性越强,扩张状态观测器(ESO)对系统扰动的适应能力越强;β01和β02很大程度影响着系统的动态性能,状态变量的估计主要受β01的影响,系统扰动的估计主要受β02的影响,β01和β02越大估计收敛的越快,但若取值过大,扩张状态观测器(ESO)输出会产生振荡发散并产生高频噪声信号;当扩张状态观测器(ESO)的参数整定足够好时,z11(k)就能足够准确地观测出系统的状态变量,z12(k)也能足够准确地估计出扰动w2(t)和未建模的扰动。
步骤3014、根据跟踪微分器(TD)和扩张状态观测器(ESO)的输出,构造非线性状态反馈控制律为:
其中,u(k)为转速环一阶自抗扰控制器1的输出控制量,对应图2 q轴电流环的转矩电流给定e11(k)为x11(k)与z11(k)的差值,u10(k)为转速环非线性状态反馈控制律的输出,β1为转速环非线性状态反馈控制律输出误差校正增益,fal(e11(k),α01,δ01)为转速环非线性状态反馈控制律的非线性组合函数且α01为转速环非线性状态反馈控制律的非线性参数,δ01为转速环非线性状态反馈控制律的滤波参数。α01,δ01,β1为可调参数,α01越小,非线性状态反馈控制律(NLSEF)对系统扰动的适应能力越强;增大β1可加快系统响应速度,但β1取值过大,会引起系统振荡甚至不稳定,因此,需要根据系统整体控制效果调整参数β1。
本实施例中,步骤3012中根据转速环跟踪微分器(TD)的输出波形并根据转速环离散时间系统最优控制函数fst[x11(k)-v(k),x12(k),r,h]对h和r进行调整,其中,转速环离散时间系统最优控制函数fst[x11(k)-v(k),x12(k),r,h]用公式表示为:
其中,a为中间变量且
步骤3013中所述α1的取值范围为0<α1<1,所述δ1的取值为0.01,β01=1/T,β02=1/(5T2);
步骤3014中所述α01的取值范围为0<α01<1,所述δ01的取值为0.01,所述β1的取值为10/b。
本实施例中,步骤三中设计所述q轴电流环一阶自抗扰控制器4的过程为:
步骤3021、将步骤一中所述q轴电流环动态数学模型变形为其中,b1为q轴电流环补偿参数且b1=1/σ,w3(t)为q轴电流环扰动且w3(t)=-Lmψrdωr/(σLr)-ω1isd;
步骤3022、将q轴电流环跟踪微分器(TD)过渡过程的离散方程表示为:
其中,v1(k)为q轴电流环一阶自抗扰控制器4的输入信号,对应图2中q轴电流环的转矩电流给定x21(k)为跟踪输入信号v1(k)的信号,x21(k+1)为x21(k)的下一时刻的变量,x22(k)为x21(k)的微分,x22(k+1)为x22(k)的下一时刻的变量,T为离散控制周期,h1为q轴电流环跟踪微分器(TD)的滤波参数,r1为决定跟踪输入信号v1(k)的速度的参数,fst[x21(k)-v1(k),x22(k),r1,h1]为q轴电流环离散时间系统最优控制函数;跟踪微分器(TD)中,h1和r1为可调参数,h1越大,对噪声滤波效果越好;r1越大跟踪速度越快,但同时也会增加超调量;
步骤3023、将视为被控对象的已知部分,构造扩张状态观测器(ESO)为:
其中,z21(k)为对isq的估计,z21(k+1)为z21(k)的下一时刻的变量,y1(k)为isq的反馈值,对应图2中q轴电流环的转矩电流反馈isq,e20(k)为z21(k)与y1(k)的差值,z22(k)为对q轴电流环扩张状态观测器(ESO)系统的总扰动的估计,z22(k+1)为z22(k)的下一时刻的变量,fal(e20(k),α2,δ2)为q轴电流环扩张状态观测器(ESO)系统的非线性组合函数且α2为q轴电流环扩张状态观测器(ESO)系统的非线性参数,δ2为q轴电流环扩张状态观测器(ESO)的滤波参数,β11为扩张状态观测器(ESO)q轴电流估计值的输出误差校正增益、β12为q轴电流环扩张状态观测器(ESO)系统扰动估计值的输出误差校正增益,u1(k)为q轴电流环一阶自抗扰控制器4的输出控制量;α2,δ2,β11,β12为可调参数,α2越小,fal函数非线性越强,扩张状态观测器(ESO)对系统扰动的适应能力越强;β11和β12很大程度影响着系统的动态性能,状态变量的估计主要受β11的影响,系统扰动的估计主要受β12的影响,β11和β12越大估计收敛的越快,但若取值过大,扩张状态观测器(ESO)输出会产生振荡发散并产生高频噪声信号;当扩张状态观测器(ESO)的参数整定足够好时,z21(k)就能足够准确地观测出系统的状态变量,z22(k)也能足够准确地估计出扰动w3(t)和未建模的扰动。
扰动w3(t)中存在ψrd,isd形成的耦合项,若采用PI调节则耦合项将会对系统控制性能产生较大的影响,可用扩张状态观测器(ESO)对w3(t)实时进行估计和补偿。
步骤3024、根据跟踪微分器(TD)和扩张状态观测器(ESO)的输出,构造非线性状态反馈控制律为:
其中,u1(k)为q轴电流环一阶自抗扰控制器4的输出控制量,对应图2定子电压的q轴分量给定e21(k)为x21(k)与z21(k)的差值,u20(k)为q轴电流环非线性状态反馈控制律的输出,β2为q轴电流环非线性状态反馈控制律输出误差校正增益,fal(e21(k),α11,δ11)为q轴电流环非线性状态反馈控制律的非线性组合函数且α11为q轴电流环非线性状态反馈控制律的非线性参数,δ11为q轴电流环非线性状态反馈控制律的滤波参数。α11,δ11,β2为可调参数,α11越小,非线性状态反馈控制律(NLSEF)对系统扰动的适应能力越强;增大β2可加快系统响应速度,但β2取值过大,会引起系统振荡甚至不稳定,因此,需要根据系统整体控制效果调整参数β2。
本实施例中,步骤3022中根据q轴电流环跟踪微分器(TD)的输出波形并根据q轴电流环离散时间系统最优控制函数fst[x21(k)-v1(k),x22(k),r1,h1]对h1和r1进行调整,其中,q轴电流环离散时间系统最优控制函数fst[x21(k)-v1(k),x22(k),r1,h1]用公式表示为:
其中,a1为中间变量且
步骤3023中所述α2的取值范围为0<α2<1,所述δ2的取值为0.01,β11=1/T,β12=1/(5T2);
步骤3024中所述α11的取值范围为0<α11<1,所述δ11的取值为0.01,所述β2的取值为10/b1。
本实施例中,步骤三中设计所述磁链环一阶自抗扰控制器2的过程为:
步骤3031、将步骤一中所述磁链环动态数学模型变形为其中,b2为磁链环补偿参数且b2=Lm/Tr;
步骤3032、将磁链环跟踪微分器(TD)过渡过程的离散方程表示为:
其中,v2(k)为磁链环一阶自抗扰控制器2的输入信号,对应图2磁链环的磁链给定x31(k)为跟踪输入信号v2(k)的信号,x31(k+1)为x31(k)的下一时刻的变量,x32(k)为x31(k)的微分,x32(k+1)为x32(k)的下一时刻的变量,T为离散控制周期,h2为磁链环跟踪微分器(TD)的滤波参数,r2为决定跟踪输入信号v2(k)速度的参数,fst[x31(k)-v2(k),x32(k),r2,h2]为磁链环离散时间系统最优控制函数;跟踪微分器(TD)中,h2和r2为可调参数,h2越大,对噪声滤波效果越好;r2越大跟踪速度越快,但同时也会增加超调量;
步骤3033、将f1(z31(k))=-ψrd/Tr视为被控对象的已知部分,构造扩张状态观测器(ESO)为:
其中,z31(k)视为对ψrd的估计,z31(k+1)为z31(k)的下一时刻的变量,y2(k)为ψrd的反馈值,对应图2磁链环的磁链反馈ψrd,e30(k)为z31(k)与y2(k)的差值,z32(k)为对磁链环扩张状态观测器(ESO)系统的总扰动的估计,z32(k+1)为z32(k)的下一时刻的变量,fal(e30(k),α3,δ3)为磁链环扩张状态观测器(ESO)系统的非线性组合函数且α3为磁链环扩张状态观测器(ESO)的非线性参数,δ3为磁链环扩张状态观测器(ESO)的滤波参数,β21为扩张状态观测器(ESO)磁链估计值的输出误差校正增益,β22为磁链环扩张状态观测器(ESO)系统扰动估计值的输出误差校正增益,u2(k)为磁链环一阶自抗扰控制器2的输出控制量;α3,δ3,β21,β22为可调参数,α3越小,fal函数非线性越强,扩张状态观测器(ESO)对系统扰动的适应能力越强;β21和β22很大程度影响着系统的动态性能,状态变量的估计主要受β21的影响,系统扰动的估计主要受β22的影响,β21和β22越大估计收敛的越快,但若取值过大,扩张状态观测器(ESO)输出会产生振荡发散并产生高频噪声信号;当扩张状态观测器(ESO)的参数整定足够好时,z31(k)就能足够准确地观测出系统的状态变量,z32(k)也能足够准确地估计出系统的扰动。
步骤3034、根据跟踪微分器(TD)和扩张状态观测器(ESO)的输出,构造非线性状态反馈控制律为:
其中,u2(k)为磁链环一阶自抗扰控制器2的输出控制量,对应图2 d轴电流环的励磁电流给定e31(k)为x31(k)与z31(k)的差值,u30(k)为磁链环非线性状态反馈控制律的输出,β3为磁链环非线性状态反馈控制律输出误差校正增益,fal(e31(k),α21,δ21)为磁链环非线性状态反馈控制律的非线性组合函数且α21为磁链环非线性状态反馈控制律的非线性参数,δ21为磁链环非线性状态反馈控制律的滤波参数。α21,δ21,β3为可调参数,α21越小,非线性状态反馈控制律(NLSEF)对系统扰动的适应能力越强;增大β3可加快系统响应速度,但β3取值过大,会引起系统振荡甚至不稳定,因此,需要根据系统整体控制效果调整参数β3。
本实施例中,步骤3032中根据磁链环跟踪微分器(TD)的输出波形并根据磁链环离散时间系统最优控制函数fst[x31(k)-v2(k),x32(k),r2,h2]对h2和r2进行调整,其中,磁链环离散时间系统最优控制函数fst[x31(k)-v2(k),x32(k),r2,h2]用公式表示为:
其中a2为中间变量且
步骤3033中所述α3的取值范围为0<α3<1,所述δ3的取值为0.01,β21=1/T,β22=1/(5T2);
步骤3034中所述α21的取值范围为0<α21<1,所述δ21的取值为0.01,所述β3的取值为10/b2。
本实施例中,步骤三中设计所述d轴电流环一阶自抗扰控制器3的过程为:
步骤3041、将步骤一中所述d轴电流环动态数学模型变形为其中,b3为d轴电流环补偿参数且b3=1/σ,w1(t)为d轴电流环扰动且
步骤3042、将d轴电流环跟踪微分器(TD)过渡过程的离散方程表示为:
其中,v3(k)为d轴电流环一阶自抗扰控制器3的输入信号,对应图2中d轴电流环的励磁电流给定x41(k)为跟踪输入信号v3(k)的信号,x41(k+1)为x41(k)的下一时刻的变量,x42(k)为x41(k)的微分,x42(k+1)为x42(k)的下一时刻的变量,T为离散控制周期,h3为d轴电流环跟踪微分器(TD)的滤波参数,r3为决定跟踪输入信号v3(k)的速度的参数,fst[x41(k)-v3(k),x42(k),r3,h3]为d轴电流环离散时间系统最优控制函数;
跟踪微分器(TD)中,h3和r3为可调参数,h3越大,对噪声滤波效果越好;r3越大跟踪速度越快,但同时也会增加超调量;
步骤3043、将视为被控对象的已知部分,构造扩张状态观测器(ESO)为:
其中,z41(k)为对isd的估计,z41(k+1)为z41(k)的下一时刻的变量,y3(k)为isd的反馈值,对应图2中d轴电流环的励磁电流反馈isd,e40(k)为z41(k)与y3(k)的差值,z42(k)为对d轴电流环扩张状态观测器(ESO)系统的总扰动的估计,z42(k+1)为z42(k)的下一时刻的变量,fal(e40(k),α4,δ4)为d轴电流环扩张状态观测器(ESO)系统的非线性组合函数且α4为d轴电流环扩张状态观测器(ESO)的非线性参数,δ4为d轴电流环扩张状态观测器(ESO)的滤波参数,β31为扩张状态观测器(ESO)d轴电流估计值的输出误差校正增益、β32为d轴电流环扩张状态观测器(ESO)系统扰动估计值的输出误差校正增益,u3(k)为d轴电流环一阶自抗扰控制器3的输出控制量;α4,δ4,β31,β32为可调参数,α4越小,fal函数非线性越强,扩张状态观测器(ESO)对系统扰动的适应能力越强;β31和β32很大程度影响着系统的动态性能,状态变量的估计主要受β31的影响,系统扰动的估计主要受β32的影响,β31和β32越大估计收敛的越快,但若取值过大,扩张状态观测器(ESO)输出会产生振荡发散并产生高频噪声信号;当扩张状态观测器(ESO)的参数整定足够好时,z41(k)就能足够准确地观测出系统的状态变量,z42(k)也能足够准确地估计出扰动w1(t)和未建模的扰动。
扰动w1(t)中存在ω1和isq形成的耦合项以及易受电机温升、频率变化的转子电阻Rr等,这些耦合项和易变的电机参数将影响磁链闭环子系统的控制性能,进而影响整个系统的控制性能。
步骤3044、根据跟踪微分器(TD)和扩张状态观测器(ESO)的输出,构造非线性状态反馈控制律为:
其中,u3(k)为d轴电流环一阶自抗扰控制器3的输出控制量,对应图2定子电压的d轴分量给定e41(k)为x41(k)与z41(k)的差值,u40(k)为d轴电流环非线性状态反馈控制律的输出,β4为d轴电流环非线性状态反馈控制律输出误差校正增益,fal(e41(k),α31,δ31)为d轴电流环非线性状态反馈控制律的非线性组合函数且α31为d轴电流环非线性状态反馈控制律的非线性参数,δ31为d轴电流环非线性状态反馈控制律滤波参数。
α31,δ31,β4为可调参数,α31越小,非线性状态反馈控制律(NLSEF)对系统扰动的适应能力越强;增大β4可加快系统响应速度,但β4取值过大,会引起系统振荡甚至不稳定,因此,需要根据系统整体控制效果调整参数β4。
本实施例中,步骤3042中根据d轴电流环跟踪微分器(TD)的输出波形并根据d轴电流环离散时间系统最优控制函数fst对h3和r3进行调整,其中,d轴电流环离散时间系统最优控制函数fst用公式表示为:
其中,a3为中间变量且
步骤3043中所述α4的取值范围为0<α4<1,所述δ4的取值为0.01,β31=1/T,β32=1/(5T2);
步骤3044中所述α31的取值范围为0<α31<1,所述δ31的取值为0.01,所述β4的取值为10/b3。
为了验证采用自抗扰控制器(ADRC)的异步电机矢量控制系统性能,在Matlab/Simulink中搭建了系统仿真模型,并进行了仿真研究。仿真时电机的额定参数为:PN=3kW,UN=380V,IN=6.5A,fN=50Hz,Rs=1.798Ω,Rr=1.588Ω,Lm=0.387H,Ls=0.394H,Lr=0.394H,J=0.0067kg·m2,TeN=20N·m,np=2,nN=1430r/min。在Simulink仿真模型中,编写M文件S函数来实现每个自抗扰控制器的功能。
为了验证ADRC系统在抗扰性能上优于PI系统,分别仿真对比研究了抗负载扰动的性能和抗电机参数变化的性能。图3为发生负载扰动时ADRC系统与PI系统的性能对比。仿真中两系统工况相同,转速为额定转速1430r/min,负载于t=0.5s时由空载突增至额定负载20N·m,在t=1.5s时又突变为空载。从图3给出的转速响应波形可以看出:PI控制系统的调节时间ts1=0.45s,ADRC系统调节时间ts2=0.02s;PI系统的起动转速超调量σ=5.6%,ADRC系统起动转速几乎无超调;在t=0.5s两系统分别加额定负载后,PI系统的转速动态降落3.8%且恢复时间tv1=0.5s,ADRC系统的转速动态降落1.5%且恢复时间tv2=0.015s。因此,对比ADRC系统与PI系统,ADRC系统可以实现快速无超调并且对负载扰动具有更好的鲁棒性。
为研究系统的参数鲁棒性,对电机转子电阻不变以及转子电阻变为1.5Rr的两种情况,分别采用ADRC控制器与PI控制器在额定转速下空载起动进行了仿真。图4A和图4B表示了电机转子电阻发生变化的转速响应(实线为转子电阻不变的转速波形,虚线为转子电阻变为1.5Rr的转速波形),图4A为PI控制系统转速波形图,图4A为ADRC控制系统转速波形图。
由图4A和图4B可知,当转子电阻发生变化时,PI控制系统的转速动态性能发生明显变化,ADRC控制系统在转子电阻变化前后转速曲线几乎重合,系统所受影响很小。由于ADRC控制器能够有效地观测到系统参数变化引起的模型变化,并加以实时补偿和控制。因此,当电机参数发生变化时,ADRC系统具有更好的鲁棒性。
综上所述,本发明在自抗扰控制器原理的基础上分析设计了由4个一阶自抗扰控制器构成的异步电机矢量控制系统。通过仿真对比了ADRC控制器和传统PI控制器的抗扰性能。仿真结果表明,ADRC控制器较PI控制器增强了系统的抗负载扰动和抗电机参数变化的能力,提高了系统的控制精度及鲁棒性。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何限制,凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化,均仍属于本发明技术方案的保护范围内。
Claims (9)
1.一种异步电机自抗扰控制器的设计方法,其特征在于,该方法包括以下步骤:
步骤一、构建异步电机(11)在按转子磁链定向的同步旋转正交坐标系下的动态数学模型,所述动态数学模型包括转速环动态数学模型、q轴电流环动态数学模型、磁链环动态数学模型和d轴电流环动态数学模型;
步骤二、根据步骤一中构建的动态数学模型,并结合自抗扰控制原理,确定出由转速环一阶自抗扰控制器(1)、磁链环一阶自抗扰控制器(2)、d轴电流环一阶自抗扰控制器(3)和q轴电流环一阶自抗扰控制器(4)四个一阶自抗扰控制器构成的异步电机矢量控制系统的结构;
所述异步电机矢量控制系统包括转速环一阶自抗扰控制器(1)、磁链环一阶自抗扰控制器(2)、d轴电流环一阶自抗扰控制器(3)和q轴电流环一阶自抗扰控制器(4),以及PARK变换模块(5)、PARK反变换模块(6)、CLARK变换模块(7)、SVPWM模块(8)、逆变器(9)和转子磁链观测器(10);所述q轴电流环一阶自抗扰控制器(4)与转速环一阶自抗扰控制器(1)的输出端连接,所述d轴电流环一阶自抗扰控制器(3)与磁链环一阶自抗扰控制器(2)的输出端连接,所述q轴电流环一阶自抗扰控制器(4)的输出端和d轴电流环一阶自抗扰控制器(3)的输出端均与PARK反变换模块(6)的输入端连接,所述PARK反变换模块(6)的输出端与SVPWM模块(8)的输入端连接,所述SVPWM模块(8)的输出端与逆变器(9)连接,异步电机(11)与逆变器(9)的A相电流输出端、B相电流输出端和C相电流输出端均连接,所述逆变器(9)的A相电流输出端和B相电流输出端均与CLARK变换模块(7)的输入端连接,所述CLARK变换模块(7)的输出端与PARK变换模块(5)的输入端连接,所述PARK变换模块(5)的异步电机定子侧q轴电流分量输出端与q轴电流环一阶自抗扰控制器(4)连接,所述PARK变换模块(5)的异步电机定子侧d轴电流分量输出端与d轴电流环一阶自抗扰控制器(3)连接,所述PARK变换模块(5)的异步电机定子侧q轴电流分量输出端和异步电机定子侧d轴电流分量输出端均与转子磁链观测器(10)的输入端连接,所述转子磁链观测器(10)的异步电机转子侧d轴磁链信号输出端与磁链环一阶自抗扰控制器(2)连接,所述转子磁链观测器(10)的旋转定向角输出端与PARK变换模块(5)和PARK反变换模块(6)均连接;
步骤三、根据一阶自抗扰控制器由跟踪微分器、扩张状态观测器和非线性状态误差反馈控制律三部分组成的原理,分别设计步骤二中所述异步电机矢量控制系统的转速环一阶自抗扰控制器(1)、磁链环一阶自抗扰控制器(2)、d轴电流环一阶自抗扰控制器(3)和q轴电流环一阶自抗扰控制器(4);
步骤一中所述转速环动态数学模型为步骤一中所述q轴电流环动态数学模型为步骤一中所述磁链环动态数学模型为步骤一中所述d轴电流环动态数学模型为其中,Tr为异步电机转子时间常数且Tr=Lr/Rr,σ为异步电机漏磁系数且usd为异步电机定子侧d轴电压分量,usq为异步电机定子侧q轴电压分量,isd为异步电机定子侧d轴电流分量,isq为异步电机定子侧q轴电流分量;Rs为异步电机定子侧电阻,Rr为异步电机转子侧电阻,Ls为异步电机定子侧电感,Lr为异步电机转子侧电感;Lm为异步电机定转子间互感;ψrd为异步电机转子侧d轴磁链;ω1为异步电机同步角转速;ωr为异步电机转子转速;np为异步电机极对数;J为异步电机转动惯量;TL为负载转矩。
2.按照权利要求1所述的异步电机自抗扰控制器的设计方法,其特征在于:步骤三中设计所述转速环一阶自抗扰控制器(1)的过程为:
步骤3011、将步骤一中所述转速环动态数学模型变形为其中,b为转速环补偿参数且w2(t)为转速环扰动且
步骤3012、将转速环跟踪微分器过渡过程的离散方程表示为:
其中,v(k)为转速环一阶自抗扰控制器(1)的输入信号,x11(k)为跟踪输入信号v(k)的信号,x11(k+1)为x11(k)的下一时刻的变量,x12(k)为x11(k)的微分,x12(k+1)为x12(k)的下一时刻的变量,T为离散控制周期,h为转速环跟踪微分器的滤波参数,r为决定跟踪输入信号v(k)的速度的参数,fst[x11(k)-v(k),x12(k),r,h]为转速环离散时间系统最优控制函数;
步骤3013、将w2(t)视为扰动,构造转速环扩张状态观测器为:
其中,z11(k)为对ωr的估计,z11(k+1)为z11(k)的下一时刻的变量,y(k)为ωr的反馈值,e10(k)为z11(k)与y(k)的差值,z12(k)为对转速环扩张状态观测器系统的总扰动的估计,z12(k+1)为z12(k)的下一时刻的变量,fal(e10(k),α1,δ1)为转速环扩张状态观测器的非线性组合函数且α1为转速环扩张状态观测器的非线性参数,δ1为转速环扩张状态观测器的滤波参数,β01为扩张状态观测器转速估计值的输出误差校正增益,β02为转速环扩张状态观测器系统扰动估计值的输出误差校正增益,u(k)为转速环一阶自抗扰控制器(1)的输出控制量;
步骤3014、根据跟踪微分器和扩张状态观测器的输出,构造非线性状态反馈控制律为:
其中,u(k)为转速环一阶自抗扰控制器(1)的输出控制量,e11(k)为x11(k)与z11(k)的差值,u10(k)为转速环非线性状态反馈控制律的输出,β1为转速环非线性状态反馈控制律输出误差校正增益,fal(e11(k),α01,δ01)为转速环非线性状态反馈控制律的非线性组合函数且α01为转速环非线性状态反馈控制律的非线性参数,δ01为转速环非线性状态反馈控制律的滤波参数。
3.按照权利要求2所述的异步电机自抗扰控制器的设计方法,其特征在于:步骤3012中根据转速环跟踪微分器的输出波形并根据转速环离散时间系统最优控制函数对h和r进行调整,其中,转速环离散时间系统最优控制函数fst[x11(k)-v(k),x12(k),r,h]用公式表示为:
其中,a为中间变量且
步骤3013中所述α1的取值范围为0<α1<1,所述δ1的取值为0.01,β01=1/T,β02=1/(5T2);
步骤3014中所述α01的取值范围为0<α01<1,所述δ01的取值为0.01,所述β1的取值为10/b。
4.按照权利要求1所述的异步电机自抗扰控制器的设计方法,其特征在于:步骤三中设计所述q轴电流环一阶自抗扰控制器(4)的过程为:
步骤3021、将步骤一中所述q轴电流环动态数学模型变形为其中,b1为q轴电流环补偿参数且b1=1/σ,w3(t)为q轴电流环扰动且w3(t)=-Lmψrdωr/(σLr)-ω1isd;
步骤3022、将q轴电流环跟踪微分器过渡过程的离散方程表示为:
其中,v1(k)为q轴电流环一阶自抗扰控制器(4)的输入信号,x21(k)为跟踪输入信号v1(k)的信号,x21(k+1)为x21(k)的下一时刻的变量,x22(k)为x21(k)的微分,x22(k+1)为x22(k)的下一时刻的变量,T为离散控制周期,h1为q轴电流环跟踪微分器的滤波参数,r1为决定跟踪输入信号v1(k)的速度的参数,fst[x21(k)-v1(k),x22(k),r1,h1]为q轴电流环离散时间系统最优控制函数;
步骤3023、将视为被控对象的已知部分,构造扩张状态观测器为:
其中,z21(k)为对isq的估计,z21(k+1)为z21(k)的下一时刻的变量,y1(k)为isq的反馈值,e20(k)为z21(k)与y1(k)的差值,z22(k)为对q轴电流环扩张状态观测器系统的总扰动的估计,z22(k+1)为z22(k)的下一时刻的变量,fal(e20(k),α2,δ2)为q轴电流环扩张状态观测器系统的非线性组合函数且α2为q轴电流环扩张状态观测器系统的非线性参数,δ2为q轴电流环扩张状态观测器的滤波参数,β11为扩张状态观测器q轴电流估计值的输出误差校正增益、β12为q轴电流环扩张状态观测器系统扰动估计值的输出误差校正增益,u1(k)为q轴电流环一阶自抗扰控制器(4)的输出控制量;
步骤3024、根据跟踪微分器和扩张状态观测器的输出,构造非线性状态反馈控制律为:
其中,u1(k)为q轴电流环一阶自抗扰控制器(4)的输出控制量,e21(k)为x21(k)与z21(k)的差值,u20(k)为q轴电流环非线性状态反馈控制律的输出,β2为q轴电流环非线性状态反馈控制律输出误差校正增益,fal(e21(k),α11,δ11)为q轴电流环非线性状态反馈控制律的非线性组合函数且α11为q轴电流环非线性状态反馈控制律的非线性参数,δ11为q轴电流环非线性状态反馈控制律的滤波参数。
5.按照权利要求4所述的异步电机自抗扰控制器的设计方法,其特征在于:步骤3022中根据q轴电流环跟踪微分器的输出波形并根据q轴电流环离散时间系统最优控制函数fst[x21(k)-v1(k),x22(k),r1,h1]对h1和r1进行调整,其中,q轴电流环离散时间系统最优控制函数fst[x21(k)-v1(k),x22(k),r1,h1]用公式表示为:
其中,a1为中间变量且
步骤3023中所述α2的取值范围为0<α2<1,所述δ2的取值为0.01,β11=1/T,β12=1/(5T2);
步骤3024中所述α11的取值范围为0<α11<1,所述δ11的取值为0.01,所述β2的取值为10/b1。
6.按照权利要求1所述的异步电机自抗扰控制器的设计方法,其特征在于:步骤三中设计所述磁链环一阶自抗扰控制器(2)的过程为:
步骤3031、将步骤一中所述磁链环动态数学模型变形为其中,b2为磁链环补偿参数且b2=Lm/Tr;
步骤3032、将磁链环跟踪微分器过渡过程的离散方程表示为:
其中,v2(k)为磁链环一阶自抗扰控制器(2)的输入信号,x31(k)为跟踪输入信号v2(k)的信号,x31(k+1)为x31(k)的下一时刻的变量,x32(k)为x31(k)的微分,x32(k+1)为x32(k)的下一时刻的变量,T为离散控制周期,h2为磁链环跟踪微分器的滤波参数,r2为决定跟踪输入信号v2(k)速度的参数,fst[x31(k)-v2(k),x32(k),r2,h2]为磁链环离散时间系统最优控制函数;
步骤3033、将f1(z31(k))=-ψrd/Tr视为被控对象的已知部分,构造扩张状态观测器为:
其中,z31(k)视为对ψrd的估计,z31(k+1)为z31(k)的下一时刻的变量,y2(k)为ψrd的反馈值,e30(k)为z31(k)与y2(k)的差值,z32(k)为对磁链环扩张状态观测器系统的总扰动的估计,z32(k+1)为z32(k)的下一时刻的变量,fal(e30(k),α3,δ3)为磁链环扩张状态观测器系统的非线性组合函数且α3为磁链环扩张状态观测器的非线性参数,δ3为磁链环扩张状态观测器的滤波参数,β21为扩张状态观测器磁链估计值的输出误差校正增益,β22为磁链环扩张状态观测器系统扰动估计值的输出误差校正增益,u2(k)为磁链环一阶自抗扰控制器(2)的输出控制量;
步骤3034、根据跟踪微分器和扩张状态观测器的输出,构造非线性状态反馈控制律为:
其中,u2(k)为磁链环一阶自抗扰控制器(2)的输出控制量,e31(k)为x31(k)与z31(k)的差值,u30(k)为磁链环非线性状态反馈控制律的输出,β3为磁链环非线性状态反馈控制律输出误差校正增益,fal(e31(k),α21,δ21)为磁链环非线性状态反馈控制律的非线性组合函数且α21为磁链环非线性状态反馈控制律的非线性参数,δ21为磁链环非线性状态反馈控制律的滤波参数。
7.按照权利要求6所述的异步电机自抗扰控制器的设计方法,其特征在于:步骤3032中根据磁链环跟踪微分器的输出波形并根据磁链环离散时间系统最优控制函数fst[x31(k)-v2(k),x32(k),r2,h2]对h2和r2进行调整,其中,磁链环离散时间系统最优控制函数fst[x31(k)-v2(k),x32(k),r2,h2]用公式表示为:
其中a2为中间变量且
步骤3033中所述α3的取值范围为0<α3<1,所述δ3的取值为0.01,β21=1/T,β22=1/(5T2);
步骤3034中所述α21的取值范围为0<α21<1,所述δ21的取值为0.01,所述β3的取值为10/b2。
8.按照权利要求1所述的异步电机自抗扰控制器的设计方法,其特征在于:步骤三中设计所述d轴电流环一阶自抗扰控制器(3)的过程为:
步骤3041、将步骤一中所述d轴电流环动态数学模型变形为其中,b3为d轴电流环补偿参数且b3=1/σ,w1(t)为d轴电流环扰动且
步骤3042、将d轴电流环跟踪微分器过渡过程的离散方程表示为:
其中,v3(k)为d轴电流环一阶自抗扰控制器(3)的输入信号,x41(k)为跟踪输入信号v3(k)的信号,x41(k+1)为x41(k)的下一时刻的变量,x42(k)为x41(k)的微分,x42(k+1)为x42(k)的下一时刻的变量,T为离散控制周期,h3为d轴电流环跟踪微分器的滤波参数,r3为决定跟踪输入信号v3(k)的速度的参数,fst[x41(k)-v3(k),x42(k),r3,h3]为d轴电流环离散时间系统最优控制函数;
步骤3043、将视为被控对象的已知部分,构造扩张状态观测器为:
其中,z41(k)为对isd的估计,z41(k+1)为z41(k)的下一时刻的变量,y3(k)为isd的反馈值,e40(k)为z41(k)与y3(k)的差值,z42(k)为对d轴电流环扩张状态观测器系统的总扰动的估计,z42(k+1)为z42(k)的下一时刻的变量,fal(e40(k),α4,δ4)为d轴电流环扩张状态观测器系统的非线性组合函数且α4为d轴电流环扩张状态观测器的非线性参数,δ4为d轴电流环扩张状态观测器的滤波参数,β31为扩张状态观测器d轴电流估计值的输出误差校正增益、β32为d轴电流环扩张状态观测器系统扰动估计值的输出误差校正增益,u3(k)为d轴电流环一阶自抗扰控制器(3)的输出控制量;
步骤3044、根据跟踪微分器和扩张状态观测器的输出,构造非线性状态反馈控制律为:
其中,u3(k)为d轴电流环一阶自抗扰控制器(3)的输出控制量,e41(k)为x41(k)与z41(k)的差值,u40(k)为d轴电流环非线性状态反馈控制律的输出,β4为d轴电流环非线性状态反馈控制律输出误差校正增益,fal(e41(k),α31,δ31)为d轴电流环非线性状态反馈控制律的非线性组合函数且α31为d轴电流环非线性状态反馈控制律的非线性参数,δ31为d轴电流环非线性状态反馈控制律滤波参数。
9.按照权利要求8所述的异步电机自抗扰控制器的设计方法,其特征在于:步骤3042中根据d轴电流环跟踪微分器的输出波形并根据d轴电流环离散时间系统最优控制函数fst对h3和r3进行调整,其中,d轴电流环离散时间系统最优控制函数fst用公式表示为:
其中,a3为中间变量且
步骤3043中所述α4的取值范围为0<α4<1,所述δ4的取值为0.01,β31=1/T,β32=1/(5T2);
步骤3044中所述α31的取值范围为0<α31<1,所述δ31的取值为0.01,所述β4的取值为10/b3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711187196.XA CN107800343B (zh) | 2017-11-24 | 2017-11-24 | 异步电机自抗扰控制器的设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711187196.XA CN107800343B (zh) | 2017-11-24 | 2017-11-24 | 异步电机自抗扰控制器的设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107800343A CN107800343A (zh) | 2018-03-13 |
CN107800343B true CN107800343B (zh) | 2019-07-23 |
Family
ID=61535815
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711187196.XA Active CN107800343B (zh) | 2017-11-24 | 2017-11-24 | 异步电机自抗扰控制器的设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107800343B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109067217B (zh) * | 2018-09-11 | 2020-02-07 | 西安科技大学 | 三相电压型pwm整流器的线性自抗扰控制器的设计方法 |
CN109889113B (zh) * | 2019-03-11 | 2023-11-07 | 中国科学院上海技术物理研究所 | 一种基于自抗扰控制的永磁电机变速扫描控制系统 |
CN110752616A (zh) * | 2019-09-11 | 2020-02-04 | 国电南瑞科技股份有限公司 | 基于自适应自抗扰比例积分的直流输电系统控制方法及系统 |
CN110995090B (zh) * | 2019-12-12 | 2023-02-28 | 大连海事大学 | 一种控制参数未知的异步电动机位置跟踪控制器 |
CN111211718B (zh) * | 2020-01-14 | 2021-06-08 | 浙江大学 | 永磁同步电机矢量控制的自抗扰控制器参数自动调节系统 |
CN111371360A (zh) * | 2020-04-15 | 2020-07-03 | 商飞信息科技(上海)有限公司 | 一种基于抗扰观测器的三相鼠笼式异步电机控制方法 |
CN112039394A (zh) * | 2020-09-24 | 2020-12-04 | 大连海事大学 | 一种基于模糊自抗扰的pmsm伺服控制系统 |
CN112134297A (zh) * | 2020-09-27 | 2020-12-25 | 西安科技大学 | 一种静止无功发生器的自抗扰电流环控制器设计方法 |
CN112953328B (zh) * | 2021-02-06 | 2022-04-22 | 枣庄学院 | 一种电动汽车永磁同步电机自抗扰控制方法 |
CN114629397B (zh) * | 2022-01-18 | 2024-02-20 | 上海理工大学 | 一种双绕组感应电机发电系统负载自适应控制方法 |
CN115483859A (zh) * | 2022-10-12 | 2022-12-16 | 国网福建省电力有限公司 | 一种暖通系统组合电机的控制方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101977010A (zh) * | 2010-10-19 | 2011-02-16 | 天津理工大学 | 基于自抗扰控制技术的无刷双馈风力发电机解耦方法 |
CN101995822A (zh) * | 2010-09-26 | 2011-03-30 | 上海电力学院 | 一种大时滞系统的灰色自抗扰控制方法 |
CN102811015A (zh) * | 2012-08-22 | 2012-12-05 | 电子科技大学 | 一种基于自抗扰控制的交流感应电机控制系统 |
CN102868336A (zh) * | 2012-09-12 | 2013-01-09 | 江苏大学 | 基于模糊二阶自抗扰控制器的三电机同步控制系统 |
CN107359837A (zh) * | 2017-09-06 | 2017-11-17 | 湖南工业大学 | 基于滑模观测器和自抗扰控制的永磁同步电机转矩控制系统及方法 |
-
2017
- 2017-11-24 CN CN201711187196.XA patent/CN107800343B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101995822A (zh) * | 2010-09-26 | 2011-03-30 | 上海电力学院 | 一种大时滞系统的灰色自抗扰控制方法 |
CN101977010A (zh) * | 2010-10-19 | 2011-02-16 | 天津理工大学 | 基于自抗扰控制技术的无刷双馈风力发电机解耦方法 |
CN102811015A (zh) * | 2012-08-22 | 2012-12-05 | 电子科技大学 | 一种基于自抗扰控制的交流感应电机控制系统 |
CN102868336A (zh) * | 2012-09-12 | 2013-01-09 | 江苏大学 | 基于模糊二阶自抗扰控制器的三电机同步控制系统 |
CN107359837A (zh) * | 2017-09-06 | 2017-11-17 | 湖南工业大学 | 基于滑模观测器和自抗扰控制的永磁同步电机转矩控制系统及方法 |
Non-Patent Citations (1)
Title |
---|
《异步电机矢量控制及磁链观测器的研究》;于婷;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20170715;第43-47页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107800343A (zh) | 2018-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107800343B (zh) | 异步电机自抗扰控制器的设计方法 | |
CN105827168B (zh) | 基于滑模观测的永磁同步电机控制方法及系统 | |
CN102969968B (zh) | 一种永磁同步电机控制方法 | |
CN106655938B (zh) | 基于高阶滑模方法的永磁同步电机控制系统及控制方法 | |
CN108649847A (zh) | 基于频率法和模糊控制的电机pi控制器参数整定方法 | |
CN110165962B (zh) | 一种直驱永磁同步风力发电系统及其全自抗扰控制方法 | |
CN108459497B (zh) | 一种基于adrc和nlpid的稳瞄伺服系统的控制方法 | |
CN110764418A (zh) | 基于有限时间收敛扩张状态观测器的自抗扰控制器 | |
CN109066756B (zh) | 一种改善系统暂态稳定的vsc-hvdc线性自抗扰控制方法 | |
CN105577058A (zh) | 基于新型模糊自抗扰控制器的五相容错永磁电机速度控制方法 | |
CN105262393B (zh) | 一种采用新型过渡过程的容错永磁电机速度控制方法 | |
CN103312244A (zh) | 基于分段式滑模变结构的无刷直流电机直接转矩控制方法 | |
CN104639001A (zh) | 融合滑模控制和分数阶神经网络控制的伺服电机控制方法 | |
CN105811826A (zh) | 一种感应电机新型趋近律滑模控制方法 | |
CN108377117A (zh) | 基于预测控制的永磁同步电机复合电流控制系统及方法 | |
CN105790270A (zh) | 通过双馈风机转子侧变流器抑制次同步谐振的方法及装置 | |
CN103762921B (zh) | 一种基于粒子群算法的不平衡电网下dfig的多目标控制方法 | |
CN113193809A (zh) | 一种改进二阶线性自抗扰的永磁同步电机控制方法 | |
Zhu et al. | Second-order sliding-mode control of DFIG-based wind turbines | |
CN109599889A (zh) | 基于模糊自抗扰的不平衡电压下dfig低电压穿越控制方法、系统 | |
CN110266227A (zh) | 一种基于模糊滑膜变结构的永磁同步电机控制系统 | |
CN112448633B (zh) | 一种基于改进adrc的飞轮储能机侧控制系统及方法 | |
Karami et al. | High-order sliding mode control of rotor-side converter in doubly-fed wind power generation system | |
He et al. | Research on active disturbance rejection control of induction motor | |
Abulizi et al. | Research of current control strategies for doubly-fed wind power generation system |
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 |