CN110661463B - 适用于磁悬浮球形电机的分数阶pid滑模观测器设计方法 - Google Patents

适用于磁悬浮球形电机的分数阶pid滑模观测器设计方法 Download PDF

Info

Publication number
CN110661463B
CN110661463B CN201910903030.6A CN201910903030A CN110661463B CN 110661463 B CN110661463 B CN 110661463B CN 201910903030 A CN201910903030 A CN 201910903030A CN 110661463 B CN110661463 B CN 110661463B
Authority
CN
China
Prior art keywords
mode observer
fractional order
sliding mode
magnetic suspension
sliding
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
Application number
CN201910903030.6A
Other languages
English (en)
Other versions
CN110661463A (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.)
Northeastern University China
Original Assignee
Northeastern University China
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 Northeastern University China filed Critical Northeastern University China
Priority to CN201910903030.6A priority Critical patent/CN110661463B/zh
Publication of CN110661463A publication Critical patent/CN110661463A/zh
Application granted granted Critical
Publication of CN110661463B publication Critical patent/CN110661463B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/13—Observer control, e.g. using Luenberger observers or Kalman filters
    • 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/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/18—Estimation of position or speed
    • 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
    • H02P25/00—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02—Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor

Landscapes

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

Abstract

本发明提供一种适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法，涉及无位置传感器技术领域。该方法首先采集磁悬浮球形电机三相电流值和三相电压值，进而构建滑模观测器的数学模型；利用构建的滑模观测器数学模型计算磁悬浮球形电机在两相静止坐标系下的反电动势和定子电流估计值；同时以定子电流估计值的误差作为滑模面，构造滑模观测器的切换函数，并采用Prelu函数作为滑模观测器开关函数；然后利用分数阶PID对滑模观测器的切换函数中的滑模增益进行参数调节和整定；并由计算得到的反电动势得到真实转子电角度和电角速度的估计值。本发明方法有效的抑制了磁悬浮球形电机滑模观测器的抖振，能够设计出稳定的滑模观测器。

Description

适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法
技术领域
本发明涉及无位置传感器技术领域，尤其涉及一种适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法。
背景技术
近几年，随着科学技术的发展，社会工业化程度越来越高，电机作为一种在机械能和电能之间相互转换的支撑型机械，受到越来越严格的要求。传统的机械轴承电机，在一些高转速高精度低损耗要求的场合下，已经不再适用。传统电机在高速运行时，其机械式轴承与轴的摩擦会有较大阻力出现，这将引起振动、磨损以及噪音，从而导致电机发烫，寿命大大减小，同时会增加维护费用。
磁悬浮轴承技术通过产生悬浮磁力支承电机转子悬浮，使电机与轴承分离。磁轴承具备不需要润滑、没有摩擦、使用寿命长等优点，在航空航天、机械人和机械加工等领域得以应用。其与传统的机械轴承相比有一定的优势，但是也有一些问题。在电机里使用磁轴承来进行悬浮需要三个轴承，两个径向轴承和一个轴向轴承，对于电机的空间使用进行了一定的限制，并影响了电机的临界转速。针对这个问题，在电机的转矩绕组上加上径向轴承的绕组，实现电机同时旋转和悬浮的功能，这种电机又可称为无轴承电机。
长期以来，人们对单自由度电动机进行了广泛的研究，随着机械臂等多自由度复杂运动精密装置的发展，多自由度电机系统越来越受到人们的重视。磁悬浮球形电机是基于磁悬浮技术和电机技术，使其产生悬浮磁力支承电机转子悬浮，并产生磁转矩直接驱动转子转动。由于转子是球形的，磁悬浮球形电机将转矩绕组和悬浮绕组结合，也称为无轴承电机。
为了对磁悬浮球形电机进行高精度的控制，通常用位置传感器对转子的位置进行检测，加大了系统的体积和转子的转动惯量，增加了系统成本，降低了系统可靠性。因此有必要对无位置传感器技术进行研究。
目前应用于电机的转子位置和速度估计方法大致有两种，一是利用反电动势或者磁链估算转子位置，常用的有龙贝格观测器、滑模观测器和扩展卡尔曼滤波观测器。二是利用电机凸极效应的高频信号注入法。滑模观测器法(SMO)以其强鲁棒性，计算简单，便于工程应用和数字实现而备受青睐，但是存在低速抖振问题。
发明内容
本发明要解决的技术问题是针对上述现有技术的不足，提供一种适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法，消除了磁悬浮球形电机系统的低速抖振，增强了系统的鲁棒性，使传统控制系统的抗干扰能力和动态性能得到了改善。
为解决上述技术问题，本发明所采取的技术方案是：适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法，包括以下步骤：
步骤1：采集磁悬浮球形电机三相电流值iα、ib、ic和三相电压值uα、ub、uc，对三相电压值和电流值进行Clark变换得到两相静止坐标系下的电压uα、uβ和电流iα和iβ；构造在两相静止坐标系下的电压方程，并对电压方程进行变形得到电流iα和iβ的平衡方程，由此得到滑模观测器的数学模型；
所述对三相电压值和电流值进行Clark变换得到两相静止坐标系下的电压uα、uβ和电流iα和iβ，如下公式所示：
Figure BDA0002212426060000021
Figure BDA0002212426060000022
Figure BDA0002212426060000023
所述构造的两相静止坐标系下的电压方程，如下公式所示：
Figure BDA0002212426060000024
其中，Rs为定子电阻，Ls为定子电感，p为微分算子，ω为转子的电角速度，ψf为永磁体与定子交链的磁链，θ为转子电角度；
对电压方程进行变形得到两相坐标系下电流iα和iβ的平衡方程，如下公式所示：
Figure BDA0002212426060000025
对平衡方程进行简化得：
pis＝A·is+B·vs+ψf·vi
其中，is为
Figure BDA0002212426060000026
a为
Figure BDA0002212426060000027
B为
Figure BDA0002212426060000028
vs为
Figure BDA0002212426060000029
vi为
Figure BDA00022124260600000210
E为2×2的单位矩阵；
进一步简化得到滑模观测器的数学模型，如下公式所示：
Figure BDA0002212426060000031
其中，
Figure BDA0002212426060000032
为微分定子电流估计值，
Figure BDA0002212426060000033
为定子电流估计值，ksw为滑模增益，
Figure BDA0002212426060000034
为开关函数；
步骤2：利用步骤1构建的滑模观测器数学模型计算磁悬浮球形电机在两相静止坐标系下的反电动势
Figure BDA0002212426060000035
和定子电流估计值
Figure BDA0002212426060000036
同时以定子电流估计值的误差作为滑模面，构造滑模观测器的切换函数，并采用Prelu函数作为滑模观测器开关函数；
利用两相坐标系下电流iα和iβ的平衡方程得到在两相静止坐标系下的反电动势
Figure BDA0002212426060000037
如下公式所示：
Figure BDA0002212426060000038
其中，t为时间；
利用滑模观测器的数学模型构造滑模观测器的切换函数，如下公式所示：
Figure BDA0002212426060000039
其中，
Figure BDA00022124260600000310
Figure BDA00022124260600000311
为转子电角度估计值，Prelu是参数化修正线性单元，其公式为
Figure BDA00022124260600000312
其中，ai是(1，+∞)的固定参数；
步骤3：构造分数阶PID的传递函数，利用改进Oustaloup近似法将分数阶PID控制器近似成连续的高阶整数阶PID控制器，再构造分数阶微分积分算子模型，并对步骤2中得到的滑模观测器的切换函数中的滑模增益ksw进行参数调节和整定；
步骤3.1：构造分数阶PID控制器的传递函数，如下公式所示：
Figure BDA00022124260600000313
其中，G1(s)为分数阶PID控制器的传递函数，kp、ki、kd为分数阶PID控制器的比例、积分、微分环节，λ、μ为积分阶次和微分阶次，s为频域变量；
步骤3.2：构造分数阶PID控制器的微分方程，如下公式所示：
u(t)＝kp·e(t)+ki·D-λ·e(t)+kd·Dμ·e(t)
其中，u(t)为分数阶PID控制器的输出，e(t)为分数阶PID控制器的误差输入，D为微积分算子；
步骤3.3：利用改进Oustaloup近似法作为滤波器，将分数阶PID控制器近似为连续的高阶整数阶PID控制器；
构造改进Oustaloup近似法的传递函数如下公式所示：
Figure BDA0002212426060000041
其中，α、d、b为常数，ωb、ωh分别为滤波器上限截止频率和下限截止频率；
将k(s)泰勒展开，舍去高阶项，化简构造的改进Oustaloup近似法连续滤波器的传递函数，如下公式所示：
Figure BDA0002212426060000042
其中，G2(s)为改进Oustaloup近似法连续滤波器的传递函数，
Figure BDA0002212426060000043
Figure BDA0002212426060000044
N为滤波器阶次；
步骤3.4：将改进Oustaloup近似法拟合出的传递函数与一个新的低通滤波器
Figure BDA0002212426060000045
封装为分数阶微分积分算子模型；
步骤3.5：在数值仿真软件Matlab/simulink中搭建步骤3.4中的分数阶微分积分算子模型，将滑模观测器的切换函数中的滑模增益ksw作为分数阶微分积分算子模型的输入，对其进行参数整定；
步骤4：由步骤2计算得到的反电动势
Figure BDA0002212426060000046
构造反电动势模型，由反电动势模型计算得到转子电角度的估计值
Figure BDA0002212426060000047
再对转子的电角度进行补偿，得到真实转子电角度θ；并由反电动势模型得到转子磁链
Figure BDA0002212426060000048
步骤4.1：由步骤2计算得到的反电动势
Figure BDA0002212426060000049
构造反电动势模型，如下公式所示：
Figure BDA0002212426060000051
由反电动势模型计算得到的转子电角度的估计值
Figure BDA0002212426060000052
如下公式所示：
Figure BDA0002212426060000053
步骤4.2：对步骤4.1得到的转子电角度进行补偿并变形，得到真实转子电角度，如下公式所示：
Figure BDA0002212426060000054
其中，θef为采集磁悬浮球形电机的三相电流和电压经过低通滤波器时产生的相位延迟误差；
步骤4.3：由反电动势模型得到转子磁链
Figure BDA0002212426060000055
如下公式所示：
Figure BDA0002212426060000056
步骤5：通过对步骤4中得到的反电动势
Figure BDA0002212426060000057
和转子磁链
Figure BDA0002212426060000058
计算得到转子的电角速度的估计值
Figure BDA0002212426060000059
如下公式所示：
Figure BDA00022124260600000510
步骤6：利用步骤2中磁悬浮球形电机的电流估计值
Figure BDA00022124260600000511
和电流实际值构造Lyapunov函数，对所设计的滑模观测器进行稳定性分析；
构造Lyapunov函数，如下公式所示：
Figure BDA00022124260600000512
其中，
Figure BDA00022124260600000513
Figure BDA00022124260600000514
ω分别为转子电角速度的估计值和实际值；
根据Lyapunov稳定性定理，当
Figure BDA00022124260600000515
时，判定滑模观测器是稳定的。
采用上述技术方案所产生的有益效果在于：本发明提供的适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法：
(1)采用了参数化修正线性单元Prelu代替了符号函数，减弱抖振现象，其计算速度和收敛速度相比于常用的sigmod函数和tanh函数更快，具有更强的鲁棒性。
(2)采用改进Oustaloup近似法，将分数阶控制器近似成连续的高阶整数阶系统，克服了Oustaloup近似法在频段端点处的边界拟合效果不理想的缺点。利用分数阶PID对滑模观测器的切换函数中的滑模增益进行整定，有效的抑制了磁悬浮球形电机滑模观测器的抖振。
(3)利用前馈控制对转子电角度进行相位补偿，减小了系统误差对观测结果所造成的影响。
附图说明
图1为本发明实施例提供的适用的磁悬浮球形电机示意图；
图2为本发明实施例提供的适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法的流程图；
图3为本发明实施例提供的设计滑模观测器的系统图；
图4为本发明实施例提供的分数阶PID控制器的原理框图；
图5为本发明实施例提供的分数阶微分积分算子模型的示意图；
图6为本发明实施例提供的分数阶PID控制器的实现图；
图7为本发明实施例提供的将本采用本发明方法设计的滑模观测器对磁悬浮球形电机进行整体控制的系统图。
图中：1、转子转轴；2、永磁体；3、球形转子；4、定子；5、转矩绕组；6、径向绕组；7、托片；8、轴向绕组；9、底座。
具体实施方式
下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明，但不用来限制本发明的范围。
磁悬浮球形电机如图1所示，其中，1为转子转轴，由球形转子3带动，径向绕组6和转矩绕组5共同镶嵌在定子上，同时提供转矩和径向悬浮力，定子底部增加有保护装置托片7当电机停转时，托片有缓冲转子掉落的作用。轴向绕组8提供轴向悬浮力，使转子与定子4之间没有接触，没有摩擦。
本实施例中，采集磁悬浮球形电机三相电流电压构造滑模观测器数学模型，利用电流偏差作为滑模面，构造切换函数，同时用分数阶PID整定滑模观测器的滑模增益，同时对转子位置进行补偿，得到转子的电角度和电角速度，实现适用于磁悬浮球形电机的分数阶PID滑模观测器的设计。适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法，如图2和3所示，包括以下步骤：
步骤1：采集磁悬浮球形电机三相电流值iα、ib、ic和三相电压值uα、ub、uc，对三相电压值和电流值进行Clark变换得到两相静止坐标系下的电压uα、uβ和电流iα和iβ；构造在两相静止坐标系下的电压方程，并对电压方程进行变形得到电流iα和iβ的平衡方程，由此得到滑模观测器的数学模型；
所述对三相电压值和电流值进行Clark变换得到两相静止坐标系下的电压uα、uβ和电流iα和iβ，如下公式所示：
Figure BDA0002212426060000071
Figure BDA0002212426060000072
所述构造的两相静止坐标系下的电压方程，如下公式所示：
Figure BDA0002212426060000073
其中，Rs为定子电阻，Ls为定子电感，p为微分算子，ω为转子的电角速度，ψf为永磁体与定子交链的磁链，θ为转子电角度；
对电压方程进行变形得到两相坐标系下电流iα和iβ的平衡方程，如下公式所示：
Figure BDA0002212426060000074
对平衡方程进行简化得：
pis＝A·is+B·vs+ψf·vi
其中，is为
Figure BDA0002212426060000075
A为
Figure BDA0002212426060000076
B为
Figure BDA0002212426060000077
vs为
Figure BDA0002212426060000078
vi为
Figure BDA0002212426060000079
E为2×2的单位矩阵；
进一步简化得到滑模观测器的数学模型，如下公式所示：
Figure BDA00022124260600000710
其中，
Figure BDA00022124260600000711
为微分定子电流估计值，
Figure BDA00022124260600000712
为定子电流估计值，ksw为滑模增益，
Figure BDA00022124260600000713
为开关函数；
步骤2：利用步骤1构建的滑模观测器数学模型计算磁悬浮球形电机在两相静止坐标系下的反电动势
Figure BDA00022124260600000714
和定子电流估计值
Figure BDA00022124260600000715
同时以定子电流估计值的误差作为滑模面，构造滑模观测器的切换函数，并采用Prelu函数作为滑模观测器开关函数；
利用两相坐标系下电流iα和iβ的平衡方程得到在两相静止坐标系下的反电动势
Figure BDA0002212426060000081
如下公式所示：
Figure BDA0002212426060000082
其中，t为时间；
利用滑模观测器的数学模型构造滑模观测器的切换函数，如下公式所示：
Figure BDA0002212426060000083
其中，
Figure BDA0002212426060000084
Figure BDA0002212426060000085
为转子电角度估计值，Prelu是参数化修正线性单元，其公式为
Figure BDA0002212426060000086
其中，ai是(1，+∞)的固定参数；
步骤3：构造分数阶PID的传递函数，如图4所示，利用改进Oustaloup近似法将分数阶PID控制器近似成连续的高阶整数阶PID控制器，再构造分数阶微分积分算子模型，并对步骤2中得到的滑模观测器的切换函数中的滑模增益ksw进行参数调节和整定；
步骤3.1：构造分数阶PID控制器的传递函数，如下公式所示：
Figure BDA0002212426060000087
其中，G1(s)为分数阶PID控制器的传递函数，kp、ki、kd为分数阶PID控制器的比例、积分、微分环节，λ、μ为积分阶次和微分阶次，s为频域变量；
步骤3.2：构造分数阶PID控制器的微分方程，如下公式所示：
u(t)＝kp·e(t)+ki·D-λ·e(t)+kd·Dμ·e(t)
其中，u(t)为分数阶PID控制器的输出，e(t)为分数阶PID控制器的误差输入，D为微积分算子；
步骤3.3：利用改进Oustaloup近似法作为滤波器，将分数阶PID控制器近似为连续的高阶整数阶PID控制器；
构造改进Oustaloup近似法的传递函数如下公式所示：
Figure BDA0002212426060000091
其中，α、d、b为常数，ωb、ωh分别为滤波器上限截止频率和下限截止频率；
将k(s)泰勒展开，舍去高阶项，化简构造的改进Oustaloup近似法连续滤波器的传递函数，如下公式所示：
Figure BDA0002212426060000092
其中，G2(s)为改进Oustaloup近似法连续滤波器的传递函数，
Figure BDA0002212426060000093
Figure BDA0002212426060000094
N为滤波器阶次；
步骤3.4：将改进Oustaloup近似法拟合出的传递函数与一个新的低通滤波器
Figure BDA0002212426060000095
封装为分数阶微分积分算子模型，如图5所示；
步骤3.5：在数值仿真软件Matlab/simulink中搭建步骤3.4中的分数阶微分积分算子模型，如图6所示，将滑模观测器的切换函数中的滑模增益ksw作为分数阶微分积分算子模型的输入，对其进行参数整定；
步骤4：由步骤2计算得到的反电动势
Figure BDA0002212426060000096
构造反电动势模型，由反电动势模型计算得到转子电角度的估计值
Figure BDA0002212426060000097
再对转子的电角度进行补偿，得到真实转子电角度θ；并由反电动势模型得到转子磁链
Figure BDA0002212426060000098
步骤4.1：由步骤2计算得到的反电动势
Figure BDA0002212426060000099
构造反电动势模型，如下公式所示：
Figure BDA00022124260600000910
由反电动势模型计算得到的转子电角度的估计值
Figure BDA00022124260600000911
如下公式所示：
Figure BDA00022124260600000912
步骤4.2：对步骤4.1得到的转子电角度进行补偿并变形，得到真实转子电角度，如下公式所示：
Figure BDA0002212426060000101
其中，θef为采集磁悬浮球形电机的三相电流和电压经过低通滤波器时产生的相位延迟误差；
步骤4.3：由反电动势模型得到转子磁链
Figure BDA0002212426060000102
如下公式所示：
Figure BDA0002212426060000103
步骤5：通过对步骤4中得到的反电动势
Figure BDA0002212426060000104
和转子磁链
Figure BDA0002212426060000105
计算得到转子的电角速度的估计值
Figure BDA0002212426060000106
如下公式所示：
Figure BDA0002212426060000107
步骤6：利用步骤2中磁悬浮球形电机的电流估计值
Figure BDA0002212426060000108
和电流实际值构造Lyapunov函数，对所设计的滑模观测器进行稳定性分析；
构造Lyapunov函数，如下公式所示：
Figure BDA0002212426060000109
其中，
Figure BDA00022124260600001010
Figure BDA00022124260600001011
ω分别为转子电角速度的估计值和实际值；
根据Lyapunov稳定性定理，当
Figure BDA00022124260600001012
时，判定滑模观测器是稳定的。
本实施例中，将构造的稳定的滑模观测器加入磁悬浮球形电机控制系统中，如图7所示，将两相静止坐标系下的电流iα、iβ和电压uα、uβ作为滑模观测器的输入，输出位置角
Figure BDA00022124260600001013
给Park变换器，转速
Figure BDA00022124260600001014
给PI控制，将PI控制后的输出送给SVPWM控制器，输出脉冲波给三相逆变器，进而实现控制磁悬浮球形电机。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明权利要求所限定的范围。

Claims (6)

1.一种适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法，其特征在于：包括以下步骤：
步骤1：采集磁悬浮球形电机三相电流值ia、ib、ic和三相电压值ua、ub、uc，对三相电压值和电流值进行Clark变换得到两相静止坐标系下的电压uα、uβ和电流iα和iβ；构造在两相静止坐标系下的电压方程，并对电压方程进行变形得到电流iα和iβ的平衡方程，由此得到滑模观测器的数学模型；
步骤2：利用步骤1构建的滑模观测器数学模型计算磁悬浮球形电机在两相静止坐标系下的反电动势
Figure FDA0002212426050000011
和定子电流估计值
Figure FDA0002212426050000012
同时以定子电流估计值的误差作为滑模面，构造滑模观测器的切换函数，并采用Prelu函数作为滑模观测器开关函数；
步骤3：构造分数阶PID的传递函数，利用改进Oustaloup近似法将分数阶PID控制器近似成连续的高阶整数阶PID控制器，再构造分数阶微分积分算子模型，并对滑模观测器的切换函数中的滑模增益进行参数调节和整定；
步骤4：由步骤2计算得到的反电动势
Figure FDA0002212426050000013
构造反电动势模型，由反电动势模型计算得到转子电角度的估计值
Figure FDA0002212426050000014
再对转子的电角度进行补偿，得到真实转子电角度θ；并由反电动势模型得到转子磁链
Figure FDA0002212426050000015
步骤5：通过对步骤4中得到的反电动势
Figure FDA0002212426050000016
和转子磁链
Figure FDA0002212426050000017
计算得到转子的电角速度的估计值
Figure FDA0002212426050000018
如下公式所示：
Figure FDA0002212426050000019
步骤6：利用步骤2中磁悬浮球形电机的电流估计值
Figure FDA00022124260500000110
和电流实际值构造Lyapunov函数，对所设计的滑模观测器进行稳定性分析。
2.根据权利要求1所述的适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法，其特征在于：所述步骤1的具体方法为：
对采集的磁悬浮球形电机的三相电压值和电流值进行Clark变换得到两相静止坐标系下的电压uα、uβ和电流iα和iβ，如下公式所示：
Figure FDA0002212426050000021
Figure FDA0002212426050000022
所述构造的两相静止坐标系下的电压方程，如下公式所示：
Figure FDA0002212426050000023
其中，Rs为定子电阻，Ls为定子电感，p为微分算子，ω为转子的电角速度，ψf为永磁体与定子交链的磁链，θ为转子电角度；
对电压方程进行变形得到两相坐标系下电流iα和iβ的平衡方程，如下公式所示：
Figure FDA0002212426050000024
对平衡方程进行简化得：
pis＝A·is+B·vs+ψf·vi
其中，is为
Figure FDA0002212426050000025
A为
Figure FDA0002212426050000026
B为
Figure FDA0002212426050000027
vs为
Figure FDA0002212426050000028
vi为
Figure FDA0002212426050000029
E为2×2的单位矩阵；
进一步简化得到滑模观测器的数学模型，如下公式所示：
Figure FDA00022124260500000210
其中，
Figure FDA00022124260500000211
为微分定子电流估计值，
Figure FDA00022124260500000212
为定子电流估计值，ksw为滑模增益，
Figure FDA00022124260500000213
为开关函数。
3.根据权利要求2所述的适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法，其特征在于：所述步骤2的具体方法为：
利用两相坐标系下电流iα和iβ的平衡方程得到在两相静止坐标系下的反电动势
Figure FDA00022124260500000214
如下公式所示：
Figure FDA00022124260500000215
其中，t为时间；
利用滑模观测器的数学模型构造滑模观测器的切换函数，如下公式所示：
Figure FDA0002212426050000031
其中，
Figure FDA0002212426050000032
Figure FDA0002212426050000033
为转子电角度估计值，Prelu是参数化修正线性单元，其公式为
Figure FDA0002212426050000034
其中，ai是(1，+∞)的固定参数。
4.根据权利要求3所述的适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法，其特征在于：所述步骤3的具体方法为：
步骤3.1：构造分数阶PID控制器的传递函数，如下公式所示：
Figure FDA0002212426050000035
其中，G1(s)为分数阶PID控制器的传递函数，kp、ki、kd为分数阶PID控制器的比例、积分、微分环节，λ、μ为积分阶次和微分阶次，s为频域变量；
步骤3.2：构造分数阶PID控制器的微分方程，如下公式所示：
u(t)＝kp·e(t)+ki·D-λ·e(t)+kd·Dμ·e(t)
其中，u(t)为分数阶PID控制器的输出，e(t)为分数阶PID控制器的误差输入，D为微积分算子；
步骤3.3：利用改进Oustaloup近似法作为滤波器，将分数阶PID控制器近似为连续的高阶整数阶PID控制器；
构造改进Oustaloup近似法的传递函数如下公式所示：
Figure FDA0002212426050000036
其中，α、d、b为常数，ωb、ωh分别为滤波器上限截止频率和下限截止频率；
将k(s)泰勒展开，舍去高阶项，化简构造的改进Oustaloup近似法连续滤波器的传递函数，如下公式所示：
Figure FDA0002212426050000037
其中，G2(s)为改进Oustaloup近似法连续滤波器的传递函数，
Figure FDA0002212426050000041
Figure FDA0002212426050000042
N为滤波器阶次；
步骤3.4：将改进Oustaloup近似法拟合出的传递函数与一个新的低通滤波器
Figure FDA0002212426050000043
封装为分数阶微分积分算子模型；
步骤3.5：在数值仿真软件Matlab/simulink中搭建步骤3.4中的分数阶微分积分算子模型，将滑模观测器的切换函数中的滑模增益ksw作为分数阶微分积分算子模型的输入，对其进行参数整定。
5.根据权利要求4所述的适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法，其特征在于：所述步骤4的具体方法为：
步骤4.1：由步骤2计算得到的反电动势
Figure FDA0002212426050000044
构造反电动势模型，如下公式所示：
Figure FDA0002212426050000045
由反电动势模型计算得到的转子电角度的估计值
Figure FDA0002212426050000046
如下公式所示：
Figure FDA0002212426050000047
步骤4.2：对步骤4.1得到的转子电角度进行补偿并变形，得到真实转子电角度，如下公式所示：
Figure FDA0002212426050000048
其中，θef为采集磁悬浮球形电机的三相电流和电压经过低通滤波器时产生的相位延迟误差；
步骤4.3：由反电动势模型得到转子磁链
Figure FDA0002212426050000049
如下公式所示：
Figure FDA00022124260500000410
6.根据权利要求5所述的适用于磁悬浮球形电机的分数阶PID滑模观测器设计方法，其特征在于：所述步骤6的具体方法为：
构造Lyapunov函数，如下公式所示：
Figure FDA00022124260500000411
其中，
Figure FDA0002212426050000051
Figure FDA0002212426050000052
ω分别为转子电角速度的估计值和实际值；
根据Lyapunov稳定性定理，当
Figure FDA0002212426050000053
时，判定滑模观测器是稳定的。
CN201910903030.6A 2019-09-24 2019-09-24 适用于磁悬浮球形电机的分数阶pid滑模观测器设计方法 Active CN110661463B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910903030.6A CN110661463B (zh) 2019-09-24 2019-09-24 适用于磁悬浮球形电机的分数阶pid滑模观测器设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910903030.6A CN110661463B (zh) 2019-09-24 2019-09-24 适用于磁悬浮球形电机的分数阶pid滑模观测器设计方法

Publications (2)

Publication Number Publication Date
CN110661463A CN110661463A (zh) 2020-01-07
CN110661463B true CN110661463B (zh) 2021-03-09

Family

ID=69038361

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910903030.6A Active CN110661463B (zh) 2019-09-24 2019-09-24 适用于磁悬浮球形电机的分数阶pid滑模观测器设计方法

Country Status (1)

Country Link
CN (1) CN110661463B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111473049B (zh) * 2020-04-17 2021-08-20 河海大学 一种实心定子磁悬浮励磁电流的控制方法
CN114859715B (zh) * 2022-04-21 2024-04-26 武汉大学 一种磁悬浮转台运动控制系统的分数阶滑模运动控制方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103034126A (zh) * 2012-12-24 2013-04-10 江苏大学 恒流源偏置外转子轴向磁轴承的控制系统及其控制方法
CN206712688U (zh) * 2017-04-13 2017-12-05 刘若珂 基于滑模观测器的电机驱动系统

Also Published As

Publication number Publication date
CN110661463A (zh) 2020-01-07

Similar Documents

Publication Publication Date Title
CN108631672B (zh) 计及最优占空比调制的永磁同步电机预测磁链控制方法
CN108880358A (zh) 基于角位移无时滞观测器的永磁同步电机控制方法及装置
CN108512473B (zh) 三相四开关永磁同步电机调速系统的直接转矩控制方法
CN105406784B (zh) 单绕组无轴承电机转矩和悬浮力直接控制器及构造方法
CN105227025A (zh) 一种永磁同步电机低载波比无位置传感器控制系统及其控制方法
CN110661463B (zh) 适用于磁悬浮球形电机的分数阶pid滑模观测器设计方法
CN111181458A (zh) 基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法
CN108183651A (zh) 一种永磁直线同步电机无位置传感器的控制方法
CN111211717B (zh) 非奇异滑模结构的ipmsm无位置传感器电机闭环控制方法
CN103414428B (zh) 无轴承同步磁阻电机转子偏心位移控制器及其构造方法
Chen et al. Radial displacement sensorless control of bearingless flux-switching permanent magnet machines based on difference of symmetric-winding flux linkages
CN112713831A (zh) 一种基于模型预测的三相四开关逆变器永磁同步电机系统的电压控制方法
CN103427755B (zh) 一种无轴承永磁薄片电机转子径向位移控制器的构造方法
CN111756288A (zh) 一种提升永磁同步电机无位置传感器估算性能的方法
CN113659904A (zh) 一种基于非奇异快速终端滑模观测器的spmsm无传感器矢量控制方法
CN114567220A (zh) 永磁同步电机全速域无位置传感器控制方法
CN112688599B (zh) 一种提高无刷直流电机抗扰能力的复合控制方法
CN102013870A (zh) 五自由度无轴承同步磁阻电机逆系统解耦控制器
Wu et al. Sensorless model reference adaptive control of permanent magnet synchronous motor for industrial robots
CN108540031B (zh) 无轴承同步磁阻电机的转速估计方法及转矩控制系统
CN115940719A (zh) 一种新型锁相环的永磁同步电机无位置传感器控制方法
CN115001334A (zh) 基于自抗扰的无位置传感器超高速永磁同步电机转速控制方法及系统
CN113890438A (zh) 一种基于内置式永磁同步电机的无速度传感器控制方法
CN114465543A (zh) 一种永磁同步电机无位置传感器控制方法
CN113315430A (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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20200107

Assignee: Shenyang aibeike Technology Co.,Ltd.

Assignor: Northeastern University

Contract record no.: X2021210000051

Denomination of invention: Design method of Fractional Order PID sliding mode observer for magnetic levitation spherical motor

Granted publication date: 20210309

License type: Common License

Record date: 20211125

EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20200107

Assignee: Dongneng (Shenyang) Energy Engineering Technology Co.,Ltd.

Assignor: Northeastern University

Contract record no.: X2023210000291

Denomination of invention: Design method of fractional order PID sliding mode observer for magnetic levitation spherical motors

Granted publication date: 20210309

License type: Common License

Record date: 20231207

EE01 Entry into force of recordation of patent licensing contract