CN107707168A - 一种基于双电流观测器的永磁同步电机控制方法 - Google Patents

一种基于双电流观测器的永磁同步电机控制方法 Download PDF

Info

Publication number
CN107707168A
CN107707168A CN201711066158.9A CN201711066158A CN107707168A CN 107707168 A CN107707168 A CN 107707168A CN 201711066158 A CN201711066158 A CN 201711066158A CN 107707168 A CN107707168 A CN 107707168A
Authority
CN
China
Prior art keywords
msub
mrow
mover
mfrac
beta
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
CN201711066158.9A
Other languages
English (en)
Other versions
CN107707168B (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.)
Ningbo University of Technology
Original Assignee
Ningbo University 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 Ningbo University of Technology filed Critical Ningbo University of Technology
Priority to CN201711066158.9A priority Critical patent/CN107707168B/zh
Publication of CN107707168A publication Critical patent/CN107707168A/zh
Application granted granted Critical
Publication of CN107707168B publication Critical patent/CN107707168B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/22Current control, e.g. using a current control loop

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Motors That Do Not Use Commutators (AREA)
  • Control Of Ac Motors In General (AREA)
  • External Artificial Organs (AREA)

Abstract

本发明提供一种基于双电流观测器的永磁同步电机控制方法,实现电机无编码器控制技术,首先建立永磁同步电机数学模型,接着设计两种电流观测器,从这两种电流观测器进而得到反电动势值,然后利用Lyapunov函数验证系统稳定性,最后通过仿真结果确定系统参数。设计电流观测器过程中采取了滑模变结构方法,有效提高了电流估计精度,双电流观测器大大降低了外在干扰的影响,实现电机速度和位置精确估计。由于机械式编码器在永磁同步电机中的使用,不仅增加成本,而且降低可靠性,所以本发明有效克服了现有技术中的种种缺点而具高度应用价值。

Description

一种基于双电流观测器的永磁同步电机控制方法
技术领域
本发明涉及一种基于双电流观测器的永磁同步电机控制方法。
背景技术
永磁同步电机(PMSM)具备十分优良的低速性能、可以实现弱磁高速控制,调速范围宽广、动态特性和效率都很高,而且无需激磁电流,提高了电机效率和功率密度,永磁同步电机已经成为伺服系统的主流之选,广泛应用于数控机床、工业机器人等领域。
随着微电子技术、微处理器、控制技术的发展,使得很多算法复杂的控制策略可以应用到电机控制中。目前针对交流电机数学模型的控制策略主要有转差频率控制、矢量控制、解耦控制、直接转矩控制等,而其中的矢量控制、直接转矩控制以其优良控制性能受到普遍关注。然而降低成本、提高性能一直是电机控制领域追求的目标,而编码器的使用无疑带来安装、维护、干扰以及增加成本等问题,因此无传感器技术成为重要的研究方向。目前无传感器技术主要存在干外在扰的问题,还处于理论探索阶段,还有一些理论问题未解决。
近年来,无传感器技术已成为一个重要研究方向,本发明对永磁同步电机的无传感器技术进行研究。
发明内容
鉴于以上所述现有技术的缺点,为了提高永磁同步电机性能,降低成本,本发明提供一种基于双电流观测器的永磁同步电机控制方法。
步骤1:针对凸装子结构的永磁同步电机,建立数学模型,其在静止坐标系α,β中模型为
vα=-ψfωrsin(θ) (3)
vβ=-ψfωrcos(θ) (4)
其中iα,iβ分别为定子α,β轴的电流,uα,uβ分别为定子α,β轴的电压,vα,vβ为反电动势,R为定子电阻,L为定子电感,Tα,Tβ分别为定子α,β轴的电磁转矩,ψf为磁链,ωr、θ分别为转速和位置;
步骤2:考虑参数变化,设计电流观测器一为
反电动势势第一个估计值为
其中“^”表示观测量,fα,fβ表示系统输入控制,
由式(5)(6)及式(1)(2)得到误差方程为
其中“~”表示误差量,
滑模面取为:
由式(9)(10)得到等效控制
步骤3:为了进一步消除参数变化的影响,设计电流观测器二为
由式(13)(14)及式(1)(2),可得到误差方程
步骤4:综合步骤2和步骤3的两种电流观测器,可获得反电动势的第二个估计值为
其中Γ12为常数;
步骤5:取速度观测器和扰动观测器
取fα,fβ
其中γ1、γ2、γ3、γ4为常数;
步骤6:利用Lyapunov函数,验证系统稳定性;
步骤7:由MATLAB仿真结果,判断是否需要参数调整,若需要调整,返回步骤5。
综上所述,无论是直接转矩控制还是矢量控制,都需要编码器信号,然而机械式编码器在PMSM中的使用,不仅增加成本,而且降低可靠性。为了解决上述问题,本发明提出了基于双电流观测器的无传感器方法,双电流观测器大大降低了外在干扰的影响。设计电流观测器采取了滑模变结构方法,有效提高了电流估计精度,实现电机速度和位置精确估计。本方法实现简单,具有很好的应用前景。
附图说明
图1为本发明基于变结构滑模观测器的矢量控制原理图。
图2为本发明实施例速度观测器的仿真图。
图3为本发明实施例位置实际值与估计值。
图4为本发明实施例速度实际值与估计值。
图5为本发明实施例三相电流值。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地实施。
步骤1:针对凸装子结构的永磁同步电机,建立数学模型,其在静止坐标系α,β中模型为
vα=-ψfωrsin(θ) (3)
vβ=-ψfωrcos(θ) (4)
其中iα,iβ分别为定子α,β轴的电流,uα,uβ分别为定子α,β轴的电压,vα,vβ为反电动势,R为定子电阻,L为定子电感,Tα,Tβ分别为定子α,β轴的电磁转矩,ψf为磁链,ωr、θ分别为转速和位置;
步骤2:考虑参数变化,设计电流观测器一为
反电动势势第一个估计值为
其中“^”表示观测量,fα,fβ表示系统输入控制,
由式(5)(6)及式(1)(2)得到误差方程为
其中“~”表示误差量,
滑模面取为:
由式(9)(10)得到等效控制
步骤3:为了进一步消除参数变化的影响,设计电流观测器二为
由式(13)(14)及式(1)(2),可得到误差方程
步骤4:综合步骤2和步骤3的两种电流观测器,可获得反电动势的第二个估计值为
其中Γ12为常数;
步骤5:取速度观测器和扰动观测器
取fα,fβ
其中γ1、γ2、γ3、γ4为常数;
步骤6:利用Lyapunov函数,验证系统稳定性;
假设Lyapunov函数为
V=V1+V2 (26)
其中γ12为正常数,对式(27)求导
考虑(9)~(12)、(19)(20)及(29)有
由式(26)(30)有
为了渐进稳定,即
Γ1=Γ2=-1/L (33)
则式(32)变为
定义
假设由式(3)(4)(35)得到
由式(36)(37),取
在(21)-(25)满足式(38)~(40)情况下,(34)可写为
得到系统渐进稳定,简化(39)(40)为
其中γ3、γ4为较小常数;
步骤7:由MATLAB仿真结果,判断是否需要参数调整,若需要调整,返回步骤5。
图1为本发明基于变结构滑模观测器的矢量控制原理图,ASR、ACR分别表示转速和电流控制器,I_PARK、PARK、CLARK表示坐标系变换,I_PARK是dq-αβ变换,PARK是αβ-dq变换,CLARK是abc-αβ变换。图2为本发明实施例速度观测器的仿真图。仿真所用的三相PMSM参数为:定子电阻R=2.875Ω,电感Ld=Lq=8.5mH,永磁体磁通0.175Wb,转动惯量J=0.8×10- 3kg·m2,极对数pn=4,粘滞摩擦系数B=0。转速ωr=50rad/s时的仿真结果见图3~图5。图3、图4分别是ωr=50rad/s时位置和速度实际值与估计值的波形,横坐标是时间,单位秒(s),纵坐标是分别是位置和转速,位置单位弧度(rad),转速单位是弧度/秒(rad/s)。图5是使用速度观测器三相电流ia、ib、ic的波形,横坐标是时间,单位秒(s),纵坐标是电流,单位安培(A)。
综上所述,机械式编码器在PMSM中的使用,不仅增加成本,而且降低可靠性。为了解决上述问题,本发明提出了基于双电流观测器的无传感器方法,双电流观测器大大降低了外在干扰的影响。设计电流观测器采取了滑模变结构方法,有效提高了电流估计精度,实现电机速度和位置精确估计,实现方便,所以本发明有效克服了现有技术中的种种缺点而具高度应用价值。

Claims (1)

1.一种基于双电流观测器的永磁同步电机控制方法,其特征在于,其包括如下步骤:
步骤1:建立凸装子结构的永磁同步电机数学模型,其在静止坐标系α,β中模型为
<mrow> <mfrac> <mrow> <msub> <mi>di</mi> <mi>&amp;alpha;</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>R</mi> <mi>L</mi> </mfrac> <msub> <mi>i</mi> <mi>&amp;alpha;</mi> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <msub> <mi>v</mi> <mi>&amp;alpha;</mi> </msub> <mo>+</mo> <mfrac> <msub> <mi>u</mi> <mi>&amp;alpha;</mi> </msub> <mi>L</mi> </mfrac> <mo>+</mo> <msub> <mi>T</mi> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mrow> <msub> <mi>di</mi> <mi>&amp;beta;</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>R</mi> <mi>L</mi> </mfrac> <msub> <mi>i</mi> <mi>&amp;beta;</mi> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <msub> <mi>v</mi> <mi>&amp;beta;</mi> </msub> <mo>+</mo> <mfrac> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> <mi>L</mi> </mfrac> <mo>+</mo> <msub> <mi>T</mi> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
vα=-ψfωrsin(θ) (3)
vβ=-ψfωrcos(θ) (4)
其中iα,iβ分别为定子α,β轴的电流,uα,uβ分别为定子α,β轴的电压,vα,vβ为反电动势,R为定子电阻,L为定子电感,Tα,Tβ分别为定子α,β轴的电磁转矩,ψf为磁链,ωr、θ分别为转速和位置;
步骤2:考虑参数变化,设计电流观测器一为
<mrow> <mfrac> <mrow> <mi>d</mi> <msub> <mover> <mi>i</mi> <mo>^</mo> </mover> <mi>&amp;alpha;</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>R</mi> <mi>L</mi> </mfrac> <msub> <mover> <mi>i</mi> <mo>^</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <msub> <mover> <mi>v</mi> <mo>^</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>+</mo> <mfrac> <msub> <mi>u</mi> <mi>&amp;alpha;</mi> </msub> <mi>L</mi> </mfrac> <mo>+</mo> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <msub> <mi>f</mi> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mrow> <mi>d</mi> <msub> <mover> <mi>i</mi> <mo>^</mo> </mover> <mi>&amp;beta;</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>R</mi> <mi>L</mi> </mfrac> <msub> <mover> <mi>i</mi> <mo>^</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <msub> <mover> <mi>v</mi> <mo>^</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>+</mo> <mfrac> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> <mi>L</mi> </mfrac> <mo>+</mo> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <msub> <mi>f</mi> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
反电动势势第一个估计值为
<mrow> <msub> <mover> <mi>v</mi> <mo>^</mo> </mover> <mrow> <mn>1</mn> <mi>&amp;alpha;</mi> </mrow> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>&amp;psi;</mi> <mi>f</mi> </msub> <msub> <mi>&amp;omega;</mi> <mi>r</mi> </msub> <mi>sin</mi> <mrow> <mo>(</mo> <mover> <mi>&amp;theta;</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mover> <mi>v</mi> <mo>^</mo> </mover> <mrow> <mn>1</mn> <mi>&amp;beta;</mi> </mrow> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>&amp;psi;</mi> <mi>f</mi> </msub> <msub> <mover> <mi>&amp;omega;</mi> <mo>^</mo> </mover> <mi>r</mi> </msub> <mi>cos</mi> <mrow> <mo>(</mo> <mover> <mi>&amp;theta;</mi> <mo>^</mo> </mover> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
其中“^”表示观测量,fα,fβ表示系统输入控制,
由式(5)(6)及式(1)(2)得到误差方程为
<mrow> <mfrac> <mrow> <mi>d</mi> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>R</mi> <mi>L</mi> </mfrac> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <msub> <mover> <mi>v</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>+</mo> <msub> <mover> <mi>T</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <msub> <mi>f</mi> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mrow> <mi>d</mi> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>R</mi> <mi>L</mi> </mfrac> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <msub> <mover> <mi>v</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>+</mo> <msub> <mover> <mi>T</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <msub> <mi>f</mi> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
其中“~”表示误差量,
滑模面取为:
由式(9)(10)得到等效控制
<mrow> <msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>&amp;alpha;</mi> </msub> <mo>)</mo> </mrow> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <msub> <mover> <mi>v</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>+</mo> <msub> <mover> <mi>T</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>&amp;beta;</mi> </msub> <mo>)</mo> </mrow> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <msub> <mover> <mi>v</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>+</mo> <msub> <mover> <mi>T</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>
步骤3:为了进一步消除参数变化的影响,设计电流观测器二为
<mrow> <mfrac> <mrow> <mi>d</mi> <msub> <mover> <mi>i</mi> <mo>^</mo> </mover> <mi>&amp;alpha;</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>R</mi> <mi>L</mi> </mfrac> <msub> <mover> <mi>i</mi> <mo>^</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <msub> <mi>U</mi> <mi>&amp;alpha;</mi> </msub> <mo>+</mo> <mfrac> <msub> <mi>u</mi> <mi>&amp;alpha;</mi> </msub> <mi>L</mi> </mfrac> <mo>+</mo> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mrow> <mi>d</mi> <msub> <mover> <mi>i</mi> <mo>^</mo> </mover> <mi>&amp;beta;</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>R</mi> <mi>L</mi> </mfrac> <msub> <mover> <mi>i</mi> <mo>^</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <msub> <mi>U</mi> <mi>&amp;beta;</mi> </msub> <mo>+</mo> <mfrac> <msub> <mi>u</mi> <mi>&amp;beta;</mi> </msub> <mi>L</mi> </mfrac> <mo>+</mo> <msub> <mover> <mi>T</mi> <mo>^</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>14</mn> <mo>)</mo> </mrow> </mrow>
由式(13)(14)及式(1)(2),可得到误差方程
<mrow> <mfrac> <mrow> <mi>d</mi> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>R</mi> <mi>L</mi> </mfrac> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <mrow> <mo>(</mo> <msub> <mi>U</mi> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>&amp;alpha;</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mover> <mi>T</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <mfrac> <mrow> <mi>d</mi> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>R</mi> <mi>L</mi> </mfrac> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mi>L</mi> </mfrac> <mrow> <mo>(</mo> <msub> <mi>U</mi> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <msub> <mi>v</mi> <mi>&amp;beta;</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mover> <mi>T</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>U</mi> <mrow> <mi>&amp;alpha;</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>v</mi> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <mi>L</mi> <msub> <mover> <mi>T</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>17</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>U</mi> <mrow> <mi>&amp;beta;</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>v</mi> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <mi>L</mi> <msub> <mover> <mi>T</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>18</mn> <mo>)</mo> </mrow> </mrow>
步骤4:综合步骤2和步骤3的两种电流观测器,可获得反电动势的第二个估计值为
<mrow> <msub> <mover> <mi>v</mi> <mo>^</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>=</mo> <mi>L</mi> <mo>&amp;CenterDot;</mo> <mo>&amp;lsqb;</mo> <msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>&amp;alpha;</mi> </msub> <mo>)</mo> </mrow> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>&amp;Gamma;</mi> <mn>1</mn> </msub> <msub> <mi>U</mi> <mrow> <mi>&amp;alpha;</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mover> <mi>v</mi> <mo>^</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>=</mo> <mi>L</mi> <mo>&amp;CenterDot;</mo> <mo>&amp;lsqb;</mo> <msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mi>&amp;beta;</mi> </msub> <mo>)</mo> </mrow> <mrow> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>&amp;Gamma;</mi> <mn>2</mn> </msub> <msub> <mi>U</mi> <mrow> <mi>&amp;beta;</mi> <mi>e</mi> <mi>q</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>
其中Γ12为常数;
步骤5:取速度观测器和扰动观测器
<mrow> <msub> <mover> <mover> <mi>&amp;omega;</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>r</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&amp;psi;</mi> <mi>f</mi> </msub> <msub> <mi>&amp;gamma;</mi> <mn>1</mn> </msub> </mrow> <mi>L</mi> </mfrac> <mrow> <mo>(</mo> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mo>(</mo> <mover> <mi>&amp;theta;</mi> <mo>^</mo> </mover> <mo>)</mo> <mo>-</mo> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>&amp;CenterDot;</mo> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mo>(</mo> <mover> <mi>&amp;theta;</mi> <mo>^</mo> </mover> <mo>)</mo> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mover> <mover> <mi>T</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>=</mo> <msub> <mi>&amp;gamma;</mi> <mn>2</mn> </msub> <mo>&amp;CenterDot;</mo> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>22</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mover> <mover> <mi>T</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>=</mo> <msub> <mi>&amp;gamma;</mi> <mn>2</mn> </msub> <mo>&amp;CenterDot;</mo> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>23</mn> <mo>)</mo> </mrow> </mrow>
取fα,fβ
<mrow> <msub> <mi>f</mi> <mi>&amp;alpha;</mi> </msub> <mo>=</mo> <msub> <mi>&amp;gamma;</mi> <mn>3</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;alpha;</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>24</mn> <mo>)</mo> </mrow> </mrow>
<mrow> <msub> <mi>f</mi> <mi>&amp;beta;</mi> </msub> <mo>=</mo> <msub> <mi>&amp;gamma;</mi> <mn>4</mn> </msub> <mi>s</mi> <mi>i</mi> <mi>g</mi> <mi>n</mi> <mrow> <mo>(</mo> <msub> <mover> <mi>i</mi> <mo>~</mo> </mover> <mi>&amp;beta;</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>25</mn> <mo>)</mo> </mrow> </mrow>
其中γ1、γ2、γ3、γ4为常数;
步骤6:利用Lyapunov函数,验证系统稳定性;
步骤7:由MATLAB仿真结果,判断是否需要参数调整,若需要调整,返回步骤5。
CN201711066158.9A 2017-11-02 2017-11-02 一种基于双电流观测器的永磁同步电机控制方法 Expired - Fee Related CN107707168B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711066158.9A CN107707168B (zh) 2017-11-02 2017-11-02 一种基于双电流观测器的永磁同步电机控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711066158.9A CN107707168B (zh) 2017-11-02 2017-11-02 一种基于双电流观测器的永磁同步电机控制方法

Publications (2)

Publication Number Publication Date
CN107707168A true CN107707168A (zh) 2018-02-16
CN107707168B CN107707168B (zh) 2020-06-30

Family

ID=61176610

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711066158.9A Expired - Fee Related CN107707168B (zh) 2017-11-02 2017-11-02 一种基于双电流观测器的永磁同步电机控制方法

Country Status (1)

Country Link
CN (1) CN107707168B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111025038A (zh) * 2019-08-28 2020-04-17 湖南大学 一种电驱系统辐射电磁干扰预测方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130229135A1 (en) * 2012-03-02 2013-09-05 University Of Nebraska-Lincoln Drive systems including sliding mode observers and methods of controlling the same
CN104135197A (zh) * 2014-08-18 2014-11-05 济南大学 一种无传感器永磁同步电机调速控制策略
CN104868814A (zh) * 2015-05-06 2015-08-26 北京动力机械研究所 凸极式永磁同步电机无位置传感器的控制方法及装置
US10008966B2 (en) * 2012-03-02 2018-06-26 Deere & Company Drive systems including sliding mode observers and methods of controlling the same

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20130229135A1 (en) * 2012-03-02 2013-09-05 University Of Nebraska-Lincoln Drive systems including sliding mode observers and methods of controlling the same
US10008966B2 (en) * 2012-03-02 2018-06-26 Deere & Company Drive systems including sliding mode observers and methods of controlling the same
CN104135197A (zh) * 2014-08-18 2014-11-05 济南大学 一种无传感器永磁同步电机调速控制策略
CN104868814A (zh) * 2015-05-06 2015-08-26 北京动力机械研究所 凸极式永磁同步电机无位置传感器的控制方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
GUO SHUAI等: "Adaptive Dynamic Terminal Sliding Mode Control Method", 《2009 SECOND INTERNATIONAL CONFERENCE ON INTELLIGENT COMPUTATION TECHNOLOGY AND AUTOMATION》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111025038A (zh) * 2019-08-28 2020-04-17 湖南大学 一种电驱系统辐射电磁干扰预测方法及系统

Also Published As

Publication number Publication date
CN107707168B (zh) 2020-06-30

Similar Documents

Publication Publication Date Title
CN110224648B (zh) 永磁同步电机参数辨识和无位置传感器控制方法及系统
Bodson et al. High-performance induction motor control via input-output linearization
CN101567652B (zh) 永磁电动机启动
CN103715962B (zh) 双级矩阵变换器驱动的永磁同步电机滑模速度观测器
CN103036496B (zh) 自适应反推控制的永磁同步电机dtc系统及其控制方法
CN105871282A (zh) 一种基于电机转动惯量的控制器pi参数整定方法
CN106059419B (zh) 一种永磁同步电机并联矢量控制方案
CN103853891A (zh) 一种基于有限元分析的变参数永磁同步电动机建模方法
CN104393814B (zh) 一种永磁同步电机控制方法
Suman et al. Speed control of permanent magnet synchronous motor drive system using PI, PID, SMC and SMC plus PID controller
CN105227030A (zh) 电动机温度变化控制装置及方法
CN103812412B (zh) 用于估算感应电动机中参数的装置
CN103199787A (zh) 基于混合调节器的抗负载扰动方法及其装置
CN107086836A (zh) 一种改进的永磁同步电机弱磁调速方法
CN104467597A (zh) 一种抑制感应电动机电流振荡的v/f控制方法
CN106788049A (zh) 基于级联滑模观测器的无速度传感器转矩控制系统及方法
CN104660133A (zh) 电机控制装置及方法
CN109600091A (zh) 一种高精度速度要求电机的无速度传感器控制方法
CN106655937A (zh) 双级矩阵变换器驱动的同步磁阻电机模型预测控制方法
CN107395080A (zh) 基于级联非奇异终端滑模观测器的无速度传感器转矩控制系统及方法
CN107707168A (zh) 一种基于双电流观测器的永磁同步电机控制方法
CN109687781B (zh) 一种通过相位超前策略提高无刷直流电机最大转速的方法
CN105322848A (zh) 一种永磁同步电机的启动方法和系统
CN205725545U (zh) 一种感应电动机反馈型间接矢量控制系统
CN204539029U (zh) 永磁同步电机控制装置

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
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: 20200630

Termination date: 20201102