CN110119084A - 一种分数阶pid的无轴承电机驱动***参数整定方法 - Google Patents
一种分数阶pid的无轴承电机驱动***参数整定方法 Download PDFInfo
- Publication number
- CN110119084A CN110119084A CN201810143975.8A CN201810143975A CN110119084A CN 110119084 A CN110119084 A CN 110119084A CN 201810143975 A CN201810143975 A CN 201810143975A CN 110119084 A CN110119084 A CN 110119084A
- Authority
- CN
- China
- Prior art keywords
- bearing
- free motor
- fractional order
- order pid
- displacement
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000006073 displacement reaction Methods 0.000 claims abstract description 24
- 230000033228 biological regulation Effects 0.000 claims abstract description 5
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000013178 mathematical model Methods 0.000 claims description 6
- 230000035945 sensitivity Effects 0.000 claims description 5
- 230000004044 response Effects 0.000 claims description 4
- 230000003321 amplification Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 230000001934 delay Effects 0.000 claims description 2
- 238000005339 levitation Methods 0.000 claims description 2
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 2
- 235000013399 edible fruits Nutrition 0.000 claims 1
- 238000000465 moulding Methods 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 12
- 239000000725 suspension Substances 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 2
- 238000004146 energy storage Methods 0.000 description 2
- 238000005461 lubrication Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000005299 abrasion Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B11/00—Automatic controllers
- G05B11/01—Automatic controllers electric
- G05B11/36—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential
- G05B11/42—Automatic controllers electric with provision for obtaining particular characteristics, e.g. proportional, integral, differential for obtaining a characteristic which is both proportional and time-dependent, e.g. P. I., P. I. D.
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
本发明公开了一种分数阶PID的无轴承电机驱动***参数整定方法,该方法通过分数阶PID结合D‑分割法,选取加性不确定函数抑制外界扰动,提高稳定性,主要包括D‑分割技术模块、分数阶PID控制器、径向力控制模块、电流调节脉冲宽度调制模块、无轴承电机转子和位移感应器检测电路组成,无轴承电机转子的径向位移由位移感应器进行检测,位移感应器输出位移输出信号z,位移输出信号z与给定的参考位置信号zref进行比较,得到两者的偏差e。将得到的偏差结果e通过由D‑分割技术模块处理后的分数阶PID控制器调节后输出力信号Fz至径向力控制模块,进而输出电流参考信号iref,该信号iref经由电流调节脉冲宽度调制模块输出控制电流iz驱动无轴承电机转子的控制线圈,从而实现无轴承电机转子的闭环控制。
Description
技术领域
本发明涉及一种无轴承电机驱动***参数整定方法,可用于人造卫星、导弹等军事化装备,以及高速高精电机传动领域的飞轮储能、发电机、人工心脏泵的高可靠性稳定控制***等,属于电力传动控制设备以及机械润滑密封技术领域。
背景技术
无轴承电机是一种结构新颖的磁交流悬浮电机。相对一般的高速电机而言,具有高转速、无磨损、寿命长、无须润滑等优势,其控制***较为复杂,转子稳定悬浮性较难控制。磁轴承在飞轮储能,高速高精电机,人工心脏泵等电力传动领域,具有广泛的应用,特别适用于高温等环境恶劣场合。无轴承电机转子稳定悬浮***,本质上是一个开环不稳定的非线性***,目前无轴承悬浮***中仍然多数采用PID经典控制实现闭环控制***对其进行调节。但PID经典线性控制方式很难满足高速高精无轴承驱动***这一非线性数学模型高性能的要求,基于这一点,由分数阶微积分所提出的分数阶PID控制,可以更好的将PID控制方法进一步拓展应用到无轴承电机驱动***中。
相比PID经典线性控制,分数阶PID控制器具有时间记忆功能,能够根据跟随误差的历史信息合适地调节控制器输出,从而获得更好的控制性能,改善了动态响应性能和鲁棒性;具有分数阶微积分特性可以更好地说明无轴承电机驱动***,能够获得更好的实现其控制性能。
经检索国内外相关文献和专利,由D-分割法确定的分数阶PID的D稳定域未曾应用在无轴承电机驱动***控制中。
发明内容
本发明的目的是在于克服现有技术的不足,对于存在模型化误差的无轴承电机非线性数学模型高精度控制中,针对外界未知干扰,进行在线观测并能抑制扰动的分数阶PID稳定域转子的悬浮稳定性方法,实现对外界干扰的有效抑制,提高磁轴承的悬浮精度。
为实现上述目的,本发明的技术方案是采用以下步骤:
无轴承电机转子5的径向位移由位移感应器6进行检测,位移感应器6输出位移输出信号z,位移输出信号z与给定的参考位置信号zref进行比较,得到两者的偏差e。将得到的偏差结果e通过由D-分割技术模块1处理后的分数阶PID控制器2调节后输出力信号Fz至径向力控制模块3,进而输出电流参考信号iref,该信号iref经由电流调节脉冲宽度调制模块4输出控制电流iz驱动无轴承电机转子5的控制线圈,从而实现无轴承电机转子的闭环控制。
本发明与现有技术相比,有益效果是:
(1)针对无轴承电机控制***存在时滞现象,选取合适的加性不确定的权函数,通过边界穿越定理,将D稳定域的边界分为奇异边界(v=0,v=∞)和非奇异边界(0<v<∞),设计出满足不确定时滞***的H∞加权灵敏度约束分数阶PID的D稳定域。
(2)依靠传递函数的频率响应,采用H∞期望指标能够补偿外界未知干扰、无轴承电机控制***数学模型和实际数学模型的差异,提高了无轴承电机驱动控制***的精度。
(3)本发明对于具有时延的无轴承电机驱动控制***参数整定方法有着一定的实用性。
(4)本发明将D-分割方法和H∞期望指标运用到分数阶PID控制中,相比PID整定方法具有更好的精确性和稳定性。
下面结合附图和具体实施方式对本发明作进一步详细说明,但本发明的实施方式并不仅限于此。
附图说明
图1为本发明所述无轴承电机驱动控制的分数阶PID稳定域总体框图;
图中:1为D-分割技术模块;2为分数阶PID控制器;3为径向力控制模块;4为电流调节脉冲宽度调制模块;5为无轴承电机转子;6为位移感应器。
具体实施方式
参见图1所示,是本发明所述分数阶PID的稳定域方法的控制***总体框图,无轴承电机转子5的径向位移由位移感应器6进行检测,位移感应器6输出位移输出信号z,位移输出信号z与给定的参考位置信号zref进行比较,得到两者的偏差e。将得到的偏差结果e通过由D-分割技术模块1处理后的分数阶PID控制器2调节后输出力信号Fz至径向力控制模块3,进而输出电流参考信号iref,该信号iref经由电流调节脉冲宽度调制模块4输出控制电流iz驱动无轴承电机转子5的控制线圈,从而实现无轴承电机转子的闭环控制。
无轴承电机位移和悬浮力模型的传递函数可以表示为:
G(s)=(ms2-kz)-1e-τs (1)
其中,kz为单边拉力与该方向的放大系数,kw为功放系数,m为转子质量,s为拉普拉斯算子,τ为数字控制滞后时间常数。
H∞期望指标为:
其中,θ=∠(WA(s)K(s)/(1+GΔ(s)K(s)))。分数阶PID控制器K(s)、加性不确定的灵敏度函数WA(s)以及GΔ(s)的频率响应表达式分别为:
GΔ(jω)=R(ω)+jI(ω) (5)
kp,k1,kd分别代表比例值,积分值,微分值;X(s)和Y(s)均为加性不确定的灵敏度函数有理多项式;A(s)和B(s)分别为加性不确定的灵敏度函数的实数和虚数,R(ω)和I(ω)分别为轴向磁轴承传递函数的实数和虚数。λ和μ均为分数阶阶数。
将式(3)和式(4)代入到式(2)转换为驱动控制***的闭环特征多项式为:
式(6)中,R,I,A和B分别是R(ω),I(ω),A(s)和B(s)的简化表达。直接在式(6)中求解Hurwitz稳定域,计算量会非常大,这里用D分割法计算稳定域,将进一步精确缩小稳定域,令连续函数
s=ρ(v)eJθ(v) (7)
将式(7)带入到式(6)中,由此Hurwitz稳定域转化为D稳定域的求解。标量参数v为广义频率。
广义频率v将D稳定域的边界划分为奇异边界(v=0,v=∞)和非奇异边界(0<v<∞),令驱动控制***的闭环特征多项式(6)的实部和虚部分别等于零可以得到上述边界。如图1所示,本发明的分数阶PID的高速磁浮电机驱动***参数整定方法的具体实现步骤如下:
步骤一:D稳定域可以分为复数根边界、实数根边界、无限根边界。
步骤二:复数根边界的求解方法是令kd为(kp,ki)平面D稳定域的固定值,将式(7)带入式(6),并使其实部和虚部等于零,可得:
式(7)中,
a11=A cosθρλcos(θvλ)-A sinθρλsin(θvλ)-B cosθρλsin(θvλ)-B sinθρλcos(θvλ)-(R+A)ρλcos(θvλ)+(I+B)ρλsin(θvλ)
a12=A cosθ-B sinθρλ-(R+A)
a21=A cosθρλsin(θvλ)+A sinθρλcos(θvλ)+B cosθρλcos(θvλ)-B sinθρλsin(θvλ)-(R+A)ρλsin(θvλ)-(I+B)ρλcos(θvλ)
a22=A sinθ+B cosθ-(I+B)
b1=kd{A cosθρλ+μcos(θv(λ+μ))-A sinθρλ+μsin(θv(λ+μ))-B cosθρλ+μsin(θv(λ+μ))-(R+A)ρλ+μcos(θv(λ+μ))+(I+B)ρλ+μsin(θv(λ+μ))}-B sinθρλ+μcos(θv(λ+μ))-ρλcosθvλ
b2=kd{A cosθρλ+μsin(θv(λ+μ))+A sinθρλ+μcos(θv(λ+μ))+B cosθρλ+μsin(θv(λ+μ))-(R+A)ρλ+μsin(θv(λ+μ))-(I+B)ρλ+μcos(θv(λ+μ))}-B sinθρλ+μsin(θv(λ+μ))-ρλsinθvλ
由D稳定域的边界分析法。当v=0时,由闭环特征多项式(6)可知ki=0,kp是任意的,说明此时是奇异边界。当v=∞时的奇异边界,此时闭环特征多项式方程无解。综合以上分析可得(kp,ki)平面的参数非奇异稳定域。当0<v<∞时,此时闭环特征多项式方程有唯一连续解曲线,令闭环特征多项式方程的实部和虚部等于零可以直接求得解,分为两种情况,当存在标量函数v使得a11a22-a12a21≠0,由此可以通过kd,λ和μ,求解出kp,ki
当存在标量函数v使得a11a22-a12a21=0,式(7)可以求解得(kp,ki)平面中的若干条奇异曲线,具体表示为,
A cosθρλcos(θvλ)kp-(I+B)ρλsin(θvλ)ki+((R+A)ρλcos(θvλ)+(I+B)ρλsin(θvλ)kd=0 (9)
步骤三:实数根边界的求解方法是,设定ki是(kp,kd)平面的参数鲁棒稳定域的任意固定值,令式(6)的实部和虚部等于零,可得到
式(10)中,a31=a11,a41=a21,
a32=A cosθρλ+μcos(θv(λ+μ))-A sinθρλ+μsin(θv(λ+μ))-B cosθρλ+μsin(θv(λ+μ))-(R+A)ρλ+μcos(θv(λ+μ))+(I+B)ρλ+μsin(θv(λ+μ))
a42=A cosθρλ+μsin(θv(λ+μ))+A sinθρλ+μcos(θv(λ+μ))+B cosθρλ+μsin(θv(λ+μ))-(R+A)ρλ+μsin(θv(λ+μ))-(I+B)ρλ+μcos(θv(λ+μ))
b3=-A cosθ+B sinθρλ+(R+A)-B sinθρλ+μcos(θv(λ+μ))-ρλcosθvλ,
b4=-A sinθ-B cosθ+(I+B)-B sinθρλ+μsin(θv(λ+μ))-ρλsinθvλ
当存在标量函数v使得θ(v)=0,可知s(v)=ρ(v),求解得(kp,ki)平面的D稳定域由以下曲线得出:
步骤四:当ρ(v)→∞时,由闭环特征多项式(6)可知,(ki,kd)平面D稳定域的无限根边界是一组直线相交的凸多边形,这组直线表示为kd=aki+b,a和b为直线方程系数。
此处针对无轴承电机悬浮***的控制模型采用分数阶PID参数整定方法,可以有效控制转子的悬浮稳定性,并且通过H∞期望指标能够补偿外界未知干扰、无轴承电机控制***数学模型和实际数学模型的差异,本发明具有显著提高无轴承电机悬浮***精度的优点。
本发明说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (5)
1.一种分数阶PID的无轴承电机驱动***参数整定方法,其特征在于,本发明技术方案是:无轴承电机转子5的径向位移由位移感应器6进行检测,位移感应器6输出位移输出信号z,位移输出信号z与给定的参考位置信号zref进行比较,得到两者的偏差e。将得到的偏差结果e通过由D-分割技术模块1处理后的分数阶PID控制器2调节后输出力信号Fz至径向力控制模块3,进而输出电流参考信号iref,该信号iref经由电流调节脉冲宽度调制模块4输出控制电流iz驱动无轴承电机转子5的控制线圈,从而实现无轴承电机转子的闭环控制。
2.根据权利要求1所述的一种分数阶PID的无轴承电机驱动***参数整定方法,其特征在于,无轴承电机位移和悬浮力模型的传递函数可以表示为:
G(s)=(ms2-kz)-1e-τs
其中,kz为单边拉力与该方向的放大系数,kw为功放系数,m为转子质量,s为拉普拉斯算子,τ为数字控制滞后时间常数。
3.根据权利要求1所述的一种分数阶PID的无轴承电机驱动***参数整定方法,其特征在于,依靠传递函数的频率响应,采用H∞期望指标能够补偿外界未知干扰、无轴承电机控制***数学模型和实际数学模型的差异,提高了无轴承电机驱动控制***的精度。
4.根据权利要求1所述的一种分数阶PID的无轴承电机驱动***参数整定方法,其特征在于,针对无轴承电机控制***存在时滞现象,选取合适的加性不确定的权函数,通过边界穿越定理,将D稳定域的边界分为奇异边界(v=0,v=∞)和非奇异边界(0<v<∞),设计出满足不确定时滞***的H∞加权灵敏度约束分数阶PID的D稳定域。
5.根据权利要求1所述的一种分数阶PID的无轴承电机驱动***参数整定方法,其特征在于,本发明将D-分割方法和H∞期望指标运用到分数阶PID控制中,相比PID整定方法具有更好的精确性和稳定性。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810143975.8A CN110119084A (zh) | 2018-02-03 | 2018-02-03 | 一种分数阶pid的无轴承电机驱动***参数整定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810143975.8A CN110119084A (zh) | 2018-02-03 | 2018-02-03 | 一种分数阶pid的无轴承电机驱动***参数整定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110119084A true CN110119084A (zh) | 2019-08-13 |
Family
ID=67519638
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810143975.8A Pending CN110119084A (zh) | 2018-02-03 | 2018-02-03 | 一种分数阶pid的无轴承电机驱动***参数整定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110119084A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110824910A (zh) * | 2019-12-05 | 2020-02-21 | 苏州市职业大学 | 一种确定磁轴承的乘性pid稳定域的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073270A (zh) * | 2011-01-27 | 2011-05-25 | 浙江工业大学 | 单输入单输出时滞***的分数阶pid控制方法 |
CN103034127A (zh) * | 2012-12-24 | 2013-04-10 | 江苏大学 | 一种轴向磁轴承控制*** |
CN103513574A (zh) * | 2013-05-06 | 2014-01-15 | 江苏大学 | 轴向混合磁轴承鱼群算法分数阶pid控制器的构造方法 |
CN103529697A (zh) * | 2013-10-16 | 2014-01-22 | 江苏大学 | 一种确定轴向磁轴承pid参数鲁棒稳定域的方法 |
CN104199293A (zh) * | 2014-08-12 | 2014-12-10 | 上海交通大学 | 基于最大灵敏度指标的稳定分数阶pid参数优化方法 |
WO2018022942A1 (en) * | 2016-07-27 | 2018-02-01 | Smigelski Jeffrey R | Fractional scaling digital signal processing |
-
2018
- 2018-02-03 CN CN201810143975.8A patent/CN110119084A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073270A (zh) * | 2011-01-27 | 2011-05-25 | 浙江工业大学 | 单输入单输出时滞***的分数阶pid控制方法 |
CN103034127A (zh) * | 2012-12-24 | 2013-04-10 | 江苏大学 | 一种轴向磁轴承控制*** |
CN103513574A (zh) * | 2013-05-06 | 2014-01-15 | 江苏大学 | 轴向混合磁轴承鱼群算法分数阶pid控制器的构造方法 |
CN103529697A (zh) * | 2013-10-16 | 2014-01-22 | 江苏大学 | 一种确定轴向磁轴承pid参数鲁棒稳定域的方法 |
CN104199293A (zh) * | 2014-08-12 | 2014-12-10 | 上海交通大学 | 基于最大灵敏度指标的稳定分数阶pid参数优化方法 |
CN105045093A (zh) * | 2014-08-12 | 2015-11-11 | 上海交通大学 | 基于最大灵敏度指标的稳定分数阶pid参数优化方法 |
WO2018022942A1 (en) * | 2016-07-27 | 2018-02-01 | Smigelski Jeffrey R | Fractional scaling digital signal processing |
Non-Patent Citations (2)
Title |
---|
杨益飞: "D-分割技术确定磁轴承的PID参数鲁棒稳定域", 《电机与控制学报》 * |
杨益飞: "永磁同步发电机混沌运动分析及最优输出反馈H∞控制", 《物理学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110824910A (zh) * | 2019-12-05 | 2020-02-21 | 苏州市职业大学 | 一种确定磁轴承的乘性pid稳定域的方法 |
CN110824910B (zh) * | 2019-12-05 | 2023-03-24 | 苏州市职业大学 | 一种确定磁轴承的乘性pid稳定域的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Watil et al. | Multi-objective output feedback control strategy for a variable speed wind energy conversion system | |
Chen et al. | Design of robust MPPT controller for grid-connected PMSG-Based wind turbine via perturbation observation based nonlinear adaptive control | |
CN101938246B (zh) | 无速度传感器电机转速的模糊融合辨识方法 | |
Saarinen et al. | Field measurements and system identification of three frequency controlling hydropower plants | |
CN101581272B (zh) | 定桨距变速风力发电机组在失速区的功率控制方法 | |
Li et al. | Sensorless-based active disturbance rejection control for a wind energy conversion system with permanent magnet synchronous generator | |
CN104579090B (zh) | 一种永磁同步电机功率补偿控制***及方法 | |
CN105743107B (zh) | 一种电力孤网***频率调节的控制方法 | |
CN105871282A (zh) | 一种基于电机转动惯量的控制器pi参数整定方法 | |
CN103762921B (zh) | 一种基于粒子群算法的不平衡电网下dfig的多目标控制方法 | |
CN112542855B (zh) | 一种双馈风力发电***相量模型建模及仿真方法 | |
CN105262393A (zh) | 一种采用新型过渡过程的容错永磁电机速度控制方法 | |
CN109507869A (zh) | 一种适用于永磁同步电机的电机控制pi参数的优化方法 | |
CN105843321A (zh) | 一种开关磁阻风力发电***追踪最大功率控制方法 | |
Xiao et al. | Experimental evaluation of extremum seeking based region-2 controller for CART3 wind turbine | |
CN110784144B (zh) | 内置式永磁同步电机的改进控制方法 | |
CN110119084A (zh) | 一种分数阶pid的无轴承电机驱动***参数整定方法 | |
CN103485978A (zh) | 一种电磁转矩补偿实现最大风能快速平稳跟踪的控制方法 | |
Chakib et al. | Comparative study of active disturbance rejection control with RST control for variable wind speed turbine based on doubly fed induction generator connected to the grid | |
CN103529697A (zh) | 一种确定轴向磁轴承pid参数鲁棒稳定域的方法 | |
CN204376797U (zh) | 一种永磁同步电机功率补偿控制*** | |
Hailemariam et al. | Real-Time speed control of a PMSM for wind turbine application | |
CN115102180B (zh) | 一种基于构网型变流器的svg主动支撑控制方法及*** | |
CN104819098B (zh) | 一种无速度传感器的风力发电最大功率跟踪方法 | |
CN112072693B (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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190813 |
|
WD01 | Invention patent application deemed withdrawn after publication |