CN107370426A - 一种永磁球形电机四元数反馈线性化的运动控制方法 - Google Patents
一种永磁球形电机四元数反馈线性化的运动控制方法 Download PDFInfo
- Publication number
- CN107370426A CN107370426A CN201710470172.9A CN201710470172A CN107370426A CN 107370426 A CN107370426 A CN 107370426A CN 201710470172 A CN201710470172 A CN 201710470172A CN 107370426 A CN107370426 A CN 107370426A
- Authority
- CN
- China
- Prior art keywords
- msub
- mtd
- mrow
- mtr
- rotor
- 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
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
- H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
- H02P6/34—Modelling or simulation for control purposes
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
- External Artificial Organs (AREA)
Abstract
本发明涉及一种永磁球形电机四元数反馈线性化的运动控制方法,所适用的电机的定子线圈位于定子球壳内接正八面体的顶点、棱中点和面心处；永磁磁极的N、S极交替分布于转子球壳内接正六面体的顶点处，包括：根据PMSM定转子结构及洛伦兹法则，建立转子电磁转矩模型；根据PMSM定转子结构及楞次定律，建立定子反电动势模型；根据四元数特性，建立转子运动学模型；根据刚体牛顿‑欧拉旋转公式，建立转子动力学模型；根据运动学模型及动力学模型，建立PMSM整体运动模型；对整体运动模型进行坐标变换，得到解耦后的整体运动模型；对解耦后的整体运动模型进行反馈控制，得到闭环控制系统。
Description
技术领域
本发明属于永磁球形电机控制技术领域，尤其涉及一种将四元数及非线性系统解耦与电机运动特性相结合的控制算法。
背景技术
传统多自由度执行器由多个单自由度直驱电机和复杂机械传动装置组合而成，体积大、控制精度差且易受到外界扰动影响。与传统多自由度执行器相比，永磁球形电机(PMSM)具有体积小、精度高、响应速度快等优点，在航空航天、工业自动化和智能机器人等领域具有广泛的应用前景。由于PMSM结构特殊且发展时间较短，目前尚无统一的控制理论，因此在实际运用之前建立科学合理的动力模型及控制方案至关重要。
由于PMSM可以实现空间自由转动，为实施精确控制，需要对其转子位置进行准确描述。现有控制方法通常采用欧拉角作为球形电机的姿态描述量，但是当PMSM转动角度较大时，这种方法在控制过程中会周期性地出现奇异，增加了控制难度并且恶化了控制效果。此外，与传统电机定转子的平面排布不同，PMSM定转子为空间排布，转子磁场分布函数复杂，定子反电动势模型及转子电磁转矩模型与球形电机结构密切相关，使得PMSM数学模型具有强耦合、非线性等特点。
发明内容
本发明的目的是克服传统控制方法中由于采用欧拉角而产生的奇异，以及球形电机固有的轴间耦合，提出一种线性化的运动控制方法，技术方案如下：
一种永磁球形电机四元数反馈线性化的运动控制方法，所适用的电机的定子线圈位于定子球壳内接正八面体的顶点、棱中点和面心处，共计26个；永磁磁极的N、S极交替分布于转子球壳内接正六面体的顶点处，共计8个，控制方法如下：
(1)根据PMSM定转子结构及洛伦兹法则，建立转子电磁转矩模型：
式中Ti,j为线圈j、磁极i作用下转子电磁转矩，T为所有线圈磁极共同作用时的转子整体电磁转矩，I为定子线圈电流矢量，G为描述转矩-电流关系的系数矩阵，ri和sj分别为磁极i和线圈j的单位位置矢量，θi,j为线圈j、磁极i间的夹角，f(θi,j)为转子磁场分布的拟合函数。
(2)根据PMSM定转子结构及楞次定律，建立定子反电动势模型：
式中Ei,j为单线圈单磁极下线圈反电动势，Ej和E分别为定子各线圈反电动势及定子整体反电动势，ω为线圈j、磁极i的相对角速度，g(θi,j)为定子线圈反电动势的拟合函数。
(3)根据四元数特性，建立转子运动学模型：
式中为q基于时间的导数，q＝[q0 q1 q2 q3]T为姿态四元数，其中qi，i＝1,2,3为转子单位方向向量在XYZ轴上的投影，为姿态四元数补充值；ω＝[ωx ωyωz]T为转子角速度矢量，其中ωx、ωy和ωz为转子角速度在XYZ轴上分量。此外，Q(q)为转子运动学模型关系矩阵；此外，根据刚体牛顿-欧拉旋转公式，建立转子动力学模型：
式中J＝diag(Jx,Jy,Jz)为转子转动惯量矩阵，其中Jx、Jy和Jz分别为以X、Y和Z轴为旋转轴的转子转动惯量；T＝[Tx Ty Tz]T为施加在转子上的转矩矢量，其中Tx、Ty和Tz分别为转矩矢量在XYZ轴上分量。
(4)根据运动学模型及动力学模型，建立PMSM整体运动模型：
式中a1＝(Jy-Jz)/Jx，a2＝(Jz-Jx)/Jy，a3＝(Jx-Jy)/Jz。
(5)对整体运动模型进行坐标变换，得到解耦后的整体运动模型：
式中yi＝qi,i＝1,2,3为输出量，vi为虚拟输入量。
(6)对解耦后的整体运动模型实施反馈控制，得到闭环控制系统：
式中Kp和Kd分别为比例系数矩阵及微分系数矩阵；y＝[y1 y2 y3]T和分别为实际输出矢量与期望矢量。
本发明的技术效果如下：
1、以四元数替代欧拉角作为球形电机姿态描述量，避免了欧拉角带来的转动奇异问题，降低了电机控制难度，优化了电机控制效果。
2、采用反馈线性化方式对球形电机非线性的运动模型进行解耦处理，消除轴间耦合，最终得出简单二阶线性模型，便于后续增加优化算法。
附图说明
图1是PMSM结构示意图,(a)为定子，(b)为转子
图中：1、永磁球形电机定子部分；101、定子线圈；102、定子球壳；2、永磁球形电机转子部分；201、永磁磁极；202、转子球壳。
图2是电磁转矩ANSOFT仿真模型
图3是电磁转矩拟合结果
图4是反电动势拟合结果
图5是PMSM整体模型
图6是闭环控制框图
图7是无扰动定点运动仿真结果
图8是带扰动定点运动仿真结果
图9是连续轨迹跟踪仿真结果
图10是不连续轨迹跟踪仿真结果
具体实施方式
本发明提供了一种永磁球形电机四元数反馈线性化的运动控制方法，下面结合附图和仿真实例对本发明进行详细的描述。本发明具体实施步骤如下：
1.本发明选择如图1所示永磁球形电机作为示例，电机包括定子部分1、转子部分2和支撑部件。定子部分1由定子线圈101和定子球壳102组成，其中定子球壳102所用材质为聚乳酸(PLA)，定子线圈101由PLA圆柱上绕制铜漆包线构成，且定子线圈101位于定子球壳102内接正八面体的顶点、棱中点和面心处，共计26个；转子部分2由永磁磁极201和转子球壳202组成，转子球壳202材质为PLA，永磁磁极201是圆柱型表面粘贴式稀土钕铁硼永磁体，N、S极交替分布于转子球壳202内接正六面体的顶点处，共计8个。此外，定转子间由支撑部件固定，以确保两者之间无摩擦。
根据定转子结构参数，利用有限元软件ANSOFT建立单线圈单磁极模型，如图2所示。对定子线圈通1A电流条件下磁极所受电磁转矩进行仿真，并利用MATLAB对所得的电磁转矩进行多项式拟合，拟合结果如图3所示，所得拟合函数为：
其中，f(θi,j)为电磁转矩拟合函数，θi,j为磁极线圈间夹角，pm为拟合多项式的各阶系数。根据洛伦兹力特性，可知电磁转矩与线圈电流大小成正比例关系，因此建立单线圈单磁极下电磁转矩模型为：
其中，Ti,j为线圈j和磁极i作用时转子电磁转矩，ri和sj分别为磁极i和线圈j的单位方向矢量，ij为线圈j中所通电流。将上式转换为矩阵形式并采用叠加定理，可以得到所有线圈和磁极共同作用时转子整体电磁转矩模型为：
其中，T为转子整体电磁转矩，I为定子线圈电流矢量，G为转矩-电流系数矩阵。若已知球形电机转矩，可通过：
I＝G+T (4)
求解所需定子电流，其中G+＝GT(GGT)-1为转矩矩阵G的广义逆。
2.与步骤1类似，利用有限元软件ANSOFT建立单线圈单磁极模型，如图2所示。对定子线圈通恒定电流1A，转子旋转角速度为1rad/s条件下线圈反电动势进行仿真，并利用MATLAB对所得的反电动势进行多项式拟合，拟合结果如图4所示，所得拟合函数为：
其中，g(θi,j)为反电动势拟合函数，lm为拟合多项式的各阶系数。根据楞次定律可知，线圈反电动势大小与线圈、磁极相对角速度成正比例关系，因此建立单线圈单磁极下线圈反电动势模型为：
Ei,j＝g(θi,j)ω (6)
其中，Ei,j为线圈j和磁极i作用时线圈中反电动势，ω为线圈j和磁极i的相对角速度。将上式转换为矩阵形式并采用叠加定理，可以得到定子整体反电动势模型为：
其中，Ej和E分别为所有线圈及磁极作用时线圈j反电动势及定子反电动势矢量。
此外，参考传统直流电机，建立PMSM电压方程：
其中，U为定子电压矢量，R为定子电阻矢量，L＝diag(L1,L2,…L26)为定子电感矩阵。
3.(1)建立PMSM运动学模型
四元数最早由Hamilton提出，现在常用于表述三维空间内的坐标旋转变换。本发明中定义四元数为q＝[q0 q1 q2 q3]T，其中q1、q2和q3分别为转子单位方向向量在XYZ轴上的投影，并且有根据四元数特性，可以得到旋转公式
其中，qs为初始时刻球形电机转子位置，qr为旋转变换后转子位置，q为将qs转化成qr的旋转四元数。左右两边同时求导可以得到
其中，为q的共轭四元数，且假定ω＝[ωx ωy ωz]T为转子角速度矢量，由于和化简(10)后可以得到PMSM运动学模型：
(2)建立PMSM动力学模型
由于PMSM转子磁极与转子球体紧密连接，运动时可以将转子整体视为刚体。根据刚体牛顿-欧拉旋转公式，建立PMSM动力学模型为：
其中，J＝diag(Jx,Jy,Jz)为转子转动惯量矩阵，Jx、Jy和Jz分别为以X、Y和Z轴为旋转轴的转子转动惯量；T＝[Tx Ty Tz]T为施加在转子上的转矩矢量，在实际控制中T为电磁转矩矢量，其中Tx、Ty和Tz分别为转矩矢量在XYZ轴上分量。
综合上述，根据公式(3)、(8)、(11)及(12)，可以得到永磁球形电机整体模型如图5所示。PMSM模型包含电气模型、转矩模型、动力学模型、运动学模型和反电动势模型等五部分，其中动力学模型和运动学模型具有非线性、强耦合特性，在下文中将对这两部分进行反馈线性化解耦。
4.根据PMSM动力学模型及运动学模型，建立电机整体运动模型为：
其中，a1＝(Jy-Jz)/Jx，a2＝(Jz-Jx)/Jy，a3＝(Jx-Jy)/Jz。定义输入矢量u＝T，输出矢量y＝[q1 q2 q3]T，状态变量x＝[ωx ωy ωz q0 q1 q2 q3]T，则上式可写成：
其中，
(1)坐标变换
对于(14)，可以计算得到：
其中，L为李导数。为保证对均有||LgLfh||≠0，定义：
其中，ε为极小值。此时有相对阶r＝6。定义映射z＝Φ(x)为：
其中，ξ为外状态，η为内状态。因此有：
其中，C＝diag(C1,C2,C3)，C1＝C2＝C3＝[1 0]T。此外，
b(ξ,η)＝[b1 b2 b3 b4 b5 b6]T (19)
A(ξ,η)＝[A1 A2 A3 A4 A5 A6]T (22)
A1＝A3＝A5＝[0 0 0] (23)
进一步变换(18)中外状态部分可得：
其中，
A0＝diag(A01,A02,A03) (26)
B0＝diag(B01,B02,B03) (27)
u＝A(ξ,η)+(A0ξ+B0v-b(ξ,η)) (29)
其中，v是虚拟输入，和A(ξ,η)+分别为B0和A(ξ,η)的广义逆，且
将(30)代入(18)可得
即原PMSM运动模型解耦为三个互不关联的二阶线性系统。
(2)反馈控制
针对(31)中的二阶线性系统，采用PD算法进行闭环控制。定义误差信号为
其中，y和yd分别为实际输出与期望输出，e和分别为误差信号和差分误差信号。设定反馈控制率为：
其中，Kp和Kd分别为比例系数矩阵以及微分系数矩阵，y＝[y1 y2 y3]T和分别为实际输出矢量和期望矢量。结合公式(4)，可以得到解耦后整体闭环控制框图如图6所示。
5.为证实本方法的有效性，在MATLAB/Simulink环境下构建模块，分别进行定点运动和轨迹跟踪仿真验证。定点运动实验分为无扰动和带外部扰动两种，在无扰动仿真，初始位置定为[0 0 1]T，目标轨迹为[1 0 0]T、[-1 0 0]T、[0 0 1]T和[1 0 0]T，比例及微分系数矩阵分别设为Kp＝diag(9.5,9,10)和Kd＝diag(5,5,6)，仿真结果如图5所示。在带外部扰动仿真中，初始位置定为[0 0 1]T，目标位置为[1 0 0]T。在4s≤t≤4.5s内，施加扰动TE＝[0.03 0.03 0.03]T N·m,结果如图6所示；轨迹跟踪实验中，连续轨迹跟踪结果及不连续轨迹跟踪结果分别如图7和图8所示。较好的仿真结果验证了本方法的正确性和有效性。
Claims (1)
1.一种永磁球形电机四元数反馈线性化的运动控制方法,所适用的电机的定子线圈位于定子球壳内接正八面体的顶点、棱中点和面心处，共计26个；永磁磁极202的N、S极交替分布于转子球壳内接正六面体的顶点处，共计8个。包括下列步骤：
(1)根据PMSM定转子结构及洛伦兹法则，建立转子电磁转矩模型：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>T</mi>
<mo>=</mo>
<mstyle>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>8</mn>
</munderover>
</mstyle>
<mstyle>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>26</mn>
</munderover>
</mstyle>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>G</mi>
<mi>I</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>f</mi>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
<mfrac>
<mrow>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>s</mi>
<mi>j</mi>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>s</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
</mrow>
</mfrac>
<msub>
<mi>i</mi>
<mi>j</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msup>
<mi>cos</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfrac>
<mrow>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<msub>
<mi>s</mi>
<mrow>
<mo>.</mo>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>|</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>s</mi>
<mi>j</mi>
</msub>
<mo>|</mo>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中Ti,j为线圈j、磁极i作用下转子电磁转矩，T为所有线圈磁极共同作用时的转子整体电磁转矩，I为定子线圈电流矢量，G为描述转矩-电流关系的系数矩阵，ri和sj分别为磁极i和线圈j的单位位置矢量，θi,j为线圈j、磁极i间的夹角，f(θi,j)为转子磁场分布的拟合函数；
(2)根据PMSM定转子结构及楞次定律，建立定子反电动势模型：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>E</mi>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msub>
<mi>E</mi>
<mn>26</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>8</mn>
</munderover>
<msub>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>,</mo>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mn>26</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>g</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mi>&omega;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中Ei,j为单线圈单磁极下线圈反电动势，Ej和E分别为定子各线圈反电动势及定子整体反电动势，ω为线圈j、磁极i的相对角速度，g(θi,j)为定子线圈反电动势的拟合函数；
(3)根据四元数特性，建立转子运动学模型：
<mrow>
<mover>
<mi>q</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>3</mn>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>q</mi>
<mn>3</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&omega;</mi>
<mi>x</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&omega;</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&omega;</mi>
<mi>z</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mi>Q</mi>
<mrow>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
</mrow>
<mi>&omega;</mi>
</mrow>
式中为q基于时间的导数，q＝[q0 q1 q2 q3]T为姿态四元数，其中qi，i＝1,2,3为转子单位方向向量在XYZ轴上的投影，为姿态四元数补充值；ω＝[ωx ωyωz]T为转子角速度矢量，其中ωx、ωy和ωz为转子角速度在XYZ轴上分量；此外，Q(q)为转子运动学模型关系矩阵；此外，根据刚体牛顿-欧拉旋转公式，建立转子动力学模型：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>J</mi>
<mover>
<mi>&omega;</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>+</mo>
<mi>&Omega;</mi>
<mi>J</mi>
<mi>&omega;</mi>
<mo>=</mo>
<mi>T</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Omega;</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>&omega;</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&omega;</mi>
<mi>z</mi>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>x</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&omega;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>&omega;</mi>
<mi>x</mi>
</msub>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中J＝diag(Jx,Jy,Jz)为转子转动惯量矩阵，其中Jx、Jy和Jz分别为以X、Y和Z轴为旋转轴的转子转动惯量；T＝[Tx Ty Tz]T为施加在转子上的转矩矢量，其中Tx、Ty和Tz分别为转矩矢量在XYZ轴上分量；
(4)根据运动学模型及动力学模型，建立PMSM整体运动模型：
<mrow>
<mfrac>
<mi>d</mi>
<mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&omega;</mi>
<mi>x</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&omega;</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&omega;</mi>
<mi>z</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>q</mi>
<mn>0</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>q</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>q</mi>
<mn>2</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>q</mi>
<mn>3</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&omega;</mi>
<mi>y</mi>
</msub>
<msub>
<mi>&omega;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>&omega;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>&omega;</mi>
<mi>z</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<msub>
<mi>&omega;</mi>
<mi>x</mi>
</msub>
<msub>
<mi>&omega;</mi>
<mi>y</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mi>Q</mi>
<mrow>
<mo>(</mo>
<mi>q</mi>
<mo>)</mo>
</mrow>
<mi>&omega;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mfrac>
<mn>1</mn>
<msub>
<mi>J</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msub>
<mi>J</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msub>
<mi>J</mi>
<mi>z</mi>
</msub>
</mfrac>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mn>0</mn>
<mrow>
<mn>4</mn>
<mo>&times;</mo>
<mn>3</mn>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mi>x</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>T</mi>
<mi>z</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
式中a1＝(Jy-Jz)/Jx，a2＝(Jz-Jx)/Jy，a3＝(Jx-Jy)/Jz；
(5)对整体运动模型进行坐标变换，得到解耦后的整体运动模型：
<mrow>
<msub>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>i</mi>
</msub>
<mo>=</mo>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
</mrow>
式中yi＝qi,i＝1,2,3为输出量，vi为虚拟输入量；
(6)对解耦后的整体运动模型进行反馈控制，得到闭环控制系统：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>v</mi>
<mo>=</mo>
<msub>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mo>-</mo>
<msub>
<mi>K</mi>
<mi>p</mi>
</msub>
<mo>(</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>d</mi>
</msub>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mo>(</mo>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>-</mo>
<msub>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>K</mi>
<mi>p</mi>
</msub>
<mo>=</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>p</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>p</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>p</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>K</mi>
<mi>d</mi>
</msub>
<mo>=</mo>
<mi>d</mi>
<mi>i</mi>
<mi>a</mi>
<mi>g</mi>
<mo>(</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>d</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>d</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>d</mi>
<mn>3</mn>
</mrow>
</msub>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
式中Kp和Kd分别为比例系数矩阵及微分系数矩阵；y＝[y1 y2 y3]T和分别为实际输出矢量与期望矢量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710470172.9A CN107370426B (zh) | 2017-06-20 | 2017-06-20 | 一种永磁球形电机四元数反馈线性化的运动控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710470172.9A CN107370426B (zh) | 2017-06-20 | 2017-06-20 | 一种永磁球形电机四元数反馈线性化的运动控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107370426A true CN107370426A (zh) | 2017-11-21 |
| CN107370426B CN107370426B (zh) | 2020-04-17 |
Family
ID=60305492
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710470172.9A Expired - Fee Related CN107370426B (zh) | 2017-06-20 | 2017-06-20 | 一种永磁球形电机四元数反馈线性化的运动控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107370426B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107911060A (zh) * | 2017-12-19 | 2018-04-13 | 天津大学 | 一种永磁球形电机定子线圈激励方法 |
| CN108061855A (zh) * | 2017-11-30 | 2018-05-22 | 天津大学 | 一种基于mems传感器的球形电机转子位置检测方法 |
| CN109284569A (zh) * | 2018-10-15 | 2019-01-29 | 天津大学 | 一种永磁电机周期定子的组合不稳定分析及预测方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101527491A (zh) * | 2009-04-21 | 2009-09-09 | 安徽大学 | 三自由度运动永磁球形步进电动机 |
| US20130307450A1 (en) * | 2012-05-15 | 2013-11-21 | Robert Lyle Fuller | Electrical Control System |
| CN104724301A (zh) * | 2015-03-19 | 2015-06-24 | 天津大学 | 一种基于球形电机三自由度动量交换的卫星姿态控制方法 |
| CN105629732A (zh) * | 2016-01-29 | 2016-06-01 | 北京航空航天大学 | 一种考虑控制受限的航天器姿态输出反馈跟踪控制方法 |
| CN106292337A (zh) * | 2016-10-17 | 2017-01-04 | 安徽大学 | 基于正弦加速度函数的永磁球形电机点到点运动轨迹规划方法及其应用 |
-
2017
- 2017-06-20 CN CN201710470172.9A patent/CN107370426B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101527491A (zh) * | 2009-04-21 | 2009-09-09 | 安徽大学 | 三自由度运动永磁球形步进电动机 |
| US20130307450A1 (en) * | 2012-05-15 | 2013-11-21 | Robert Lyle Fuller | Electrical Control System |
| CN104724301A (zh) * | 2015-03-19 | 2015-06-24 | 天津大学 | 一种基于球形电机三自由度动量交换的卫星姿态控制方法 |
| CN105629732A (zh) * | 2016-01-29 | 2016-06-01 | 北京航空航天大学 | 一种考虑控制受限的航天器姿态输出反馈跟踪控制方法 |
| CN106292337A (zh) * | 2016-10-17 | 2017-01-04 | 安徽大学 | 基于正弦加速度函数的永磁球形电机点到点运动轨迹规划方法及其应用 |
Non-Patent Citations (1)
| Title |
|---|
| 张小红等: "《磁悬浮球形主动关节转子位置与姿态检测研究》", 《徐州工程学院学报（自然科学版）》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108061855A (zh) * | 2017-11-30 | 2018-05-22 | 天津大学 | 一种基于mems传感器的球形电机转子位置检测方法 |
| CN107911060A (zh) * | 2017-12-19 | 2018-04-13 | 天津大学 | 一种永磁球形电机定子线圈激励方法 |
| CN107911060B (zh) * | 2017-12-19 | 2020-04-21 | 天津大学 | 一种永磁球形电机定子线圈激励方法 |
| CN109284569A (zh) * | 2018-10-15 | 2019-01-29 | 天津大学 | 一种永磁电机周期定子的组合不稳定分析及预测方法 |
| CN109284569B (zh) * | 2018-10-15 | 2023-04-07 | 天津大学 | 一种永磁电机周期定子的组合不稳定分析及预测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107370426B (zh) | 2020-04-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103853891B (zh) | 一种基于有限元分析的变参数永磁同步电动机建模方法 | |
| CN103823379B (zh) | 基于迭代学习的高频率角振动转台滑模控制方法 | |
| CN107370426A (zh) | 一种永磁球形电机四元数反馈线性化的运动控制方法 | |
| CN109542087B (zh) | 一种永磁同步电机失磁故障模拟方法、系统及介质 | |
| CN104724301B (zh) | 一种基于球形电机三自由度动量交换的卫星姿态控制方法 | |
| Liu et al. | Adaptive backstepping sliding mode control for 3-DOF permanent magnet spherical actuator | |
| CN107179682A (zh) | 一种被动式负载模拟器及多余力矩抑制方法 | |
| CN107070102A (zh) | 电机、微型云台以及微型云台的电机控制方法 | |
| TW202107807A (zh) | 馬達及其驅動方法 | |
| Dušek et al. | Mathematical model of differentially steered mobile robot | |
| CN106292337A (zh) | 基于正弦加速度函数的永磁球形电机点到点运动轨迹规划方法及其应用 | |
| Guo et al. | Dynamic analysis and current calculation of a permanent magnet spherical motor for point‐to‐point motion | |
| Yue et al. | Modeling and robust adaptive control of a 3-axis motion simulator | |
| Shi et al. | A new adaptive analytical model for the spherical reluctance motor based on hybrid trigonometric function–power function | |
| CN105619394A (zh) | 一种基于误差四元数反馈的rov姿态控制方法 | |
| Li et al. | Identification of Euler angles of permanent magnet spherical motor rotor based on Hall sensors array | |
| Li et al. | Magnetic field model for permanent magnet spherical motor with double polyhedron structure | |
| Li et al. | Friction torque field distribution of a permanent‐magnet spherical motor based on multi‐physical field coupling analysis | |
| CN106992732B (zh) | 一种电机磁共能模型建立系统 | |
| CN206640471U (zh) | 电机及微型云台 | |
| CN108073087B (zh) | 一种多学科综合仿真平台环境下的机电伺服系统控制方法 | |
| Lee et al. | Hybrid multi-dof motor for multi-copter unmanned aerial vehicle | |
| Wang et al. | Simulation analysis of a new electric dynamic load simulator based on double-stator permanent-magnet synchronous motor | |
| Zhao et al. | Time-optimal trajectory planning of permanent magnet spherical motor based on genetic algorithm | |
| Li et al. | Magnetic field calculation and dynamics simulation of a permanent magnetic hybrid driven 3-dof motor |
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 |
Granted publication date: 20200417 Termination date: 20210620 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |