CN106685294A - 异步电机转子电阻及励磁电感解耦矫正方法 - Google Patents
异步电机转子电阻及励磁电感解耦矫正方法 Download PDFInfo
- Publication number
- CN106685294A CN106685294A CN201611201979.4A CN201611201979A CN106685294A CN 106685294 A CN106685294 A CN 106685294A CN 201611201979 A CN201611201979 A CN 201611201979A CN 106685294 A CN106685294 A CN 106685294A
- Authority
- CN
- China
- Prior art keywords
- rotor
- flux linkage
- vector
- formula
- asynchronous motor
- 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
- 230000005284 excitation Effects 0.000 title claims abstract description 58
- 238000012937 correction Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000004907 flux Effects 0.000 claims abstract description 139
- 239000013598 vector Substances 0.000 claims abstract description 128
- 230000003068 static effect Effects 0.000 claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 238000005314 correlation function Methods 0.000 claims abstract description 4
- 238000005070 sampling Methods 0.000 claims description 23
- 230000004069 differentiation Effects 0.000 claims description 12
- 238000006467 substitution reaction Methods 0.000 claims description 12
- 238000013178 mathematical model Methods 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 7
- 238000013461 design Methods 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 239000011541 reaction mixture Substances 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims description 2
- 230000001052 transient effect Effects 0.000 abstract description 3
- 230000010354 integration Effects 0.000 abstract 1
- 230000008859 change Effects 0.000 description 5
- 238000002474 experimental method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003750 conditioning effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
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/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
- Control Of Electric Motors In General (AREA)
Abstract
本发明公开了一种异步电机转子电阻及励磁电感解耦矫正方法,本方法根据采集并运算得到的两相静止坐标系αβ下电机定子电压向量V、定子电流向量i和转子电角速度ωr,通过高阶滑模观测器获得准确的反电动势e,依据反电动势和转子磁链之间的关系建立异步电机状态空间表达式,实现了转子磁链向量的准确观测,根据电流磁链模型获得转子磁链矢量计算值进而获得相位差Δθ、幅值差分析推导出参数误差和磁链误差的关联函数,获得加权系数,并利用耦矫正函数一矫正系统中给定转子电阻解耦矫正函数二矫正系统中给定励磁电感本发明克服了磁链获得的积分问题的同时,不存在磁链暂态问题,且能够获得准确的转子时间常数,实现了励磁电感的在线辨识。
Description
技术领域
本发明涉及一种异步电机参数在线矫正方法,具体地说是异步电机转子电阻及励磁电感解耦矫正方法。
背景技术
异步电机是一类高阶多变量非线性强耦合系统,矢量控制的应用使其获得了直流电机调速性能。间接矢量控制通过转速和转差频率的加和获得同步速度,进而获得磁链定向角。该矢量控制算法简单,易于工程实现,且具有较好的稳定型,因此是应用较为广泛的一种矢量控制方案。转差频率计算的准确性,完全取决于转子时间常数的准确性,而该参数受温度、磁饱和、涡流等因素的影响通常偏离其标称值,制约着矢量定向的准确和间接矢量控制系统的性能。另外,在电动汽车驱动系统中,通常需要对整车控制系统下发的转矩指令进行跟踪,使得系统控制精度对参数变化更为敏感。当系统控制用转子电阻、励磁电感等参数偏离其真实值时,转子磁场定向将会失准,系统的稳态控制精度和动态响应能力都将下降。
为提升矢量控制性能,关键在于转子电阻和励磁电感的准确辨识。为此,人们做出了各种努力,如题为“一种基于无功功率的异步电机矢量控制转子磁场准确定向方法”(陆海峰等,中国电机工程学报,2005年第25卷第6期116-120页)的文章和题为“基于转子磁链q轴分量的异步电机间接矢量控制转差频率校正”(樊扬等,中国电机工程学报,2009年第29卷第9期62-66页)的文章。这两篇文章分别提出了基于无功功率和转子磁链交轴分量的磁链定向在线校正方案,前者利用无功功率的偏差对转子时间常数进行在线校正,后者则利用转子磁链交轴分量实现校正,但这两种矫正算法都受励磁电感影响。
题为“一种利用转矩观测矫正异步电机转子磁场定向的方法”(燕俊峰等,中国电机工程学报,2015年第35卷第7期4517-4523页)的文章。该文利用转矩模型参考自适应对转差频率进行校正,但该方案受负载的影响较大。
题为“基于参数在线校正的电动汽车异步电机间接矢量控制”(张杰等,电工技术学报,2014第29卷第7期90-96页)的文章。该文通过电压的稳态模型获得磁链,并基于此设计了励磁电感和转子时间常数在线校正方案,虽然克服了磁链获得的积分问题,但其校正方案的暂态性能受到磁链计算的影响。
综上所述,现有技术均未能较好的实现参数及定向的在线矫正。
发明内容
本发明要解决的技术问题是克服上述各种参数矫正的局限性,针对电机运行过程中转子电阻及励磁电感的偏差,提供一种较为准确且快速的参数同时矫正方法。
为了实现本发明的目的,本发明提供了一种异步电机转子电阻及励磁电感解耦矫正方法,包括如下步骤:
步骤1,采集异步电机在静止坐标系αβ下的定子电压向量V、定子电流向量i和转子电角速度ωr;
步骤2,建立在静止坐标系αβ下异步电机状态空间表达式为:
在公式(1)中,e为反电动势向量,为e的微分,为转子磁链向量,为的微分,ωs为转子磁链电角速度,Lm为励磁电感,Lr为转子电感,Ls为定子电感,αr=Rr/Lr为转子时间常数倒数,ε=δLsLr/Lm,其中Rr为转子电阻,Rs为定子电阻, 为反对称矩阵;
步骤3,根据公式(2)建立转子磁链滑模观测器,对转子磁链向量进行观测:
在公式(2)中,为反电动势向量e的观测值,为的微分,为转子磁链向量观测值,为的微分, 为反电动势误差,为系统中给定励磁电感,为系统中给定转子电感,为系统中给定定子电感,为系统中给定转子时间常数倒数,其中为系统中给定转子电阻,为系统中给定定子电阻,k为设定的观测器增益一,g为设定的观测器增益二,k和g的数值均为负数;
步骤4,根据步骤1中得到定子电压向量V、定子电流向量i和转子电角速度ωr,经过参数偏差时电流磁链模型计算获得转子磁链矢量计算值
步骤5,根据步骤3中获得的转子磁链向量观测值和步骤4中获得的转子磁链矢量计算值经过叉乘获得相位差正弦量sin(Δθ),并设相位差Δθ≈sin(Δθ),其中所述相位差 为转子磁链矢量计算值的相位,θ为转子磁链向量观测值的相位;
步骤6,根据步骤3中获得的转子磁链向量观测值和步骤4中获得的转子磁链矢量计算值两者幅值作差获得幅值差
步骤7,利用电流磁链模型推导获得加权系数一f1、加权系数二f2、加权系数三f3、加权系数四f4;
步骤8,根据步骤5中获得的相位差Δθ、步骤6中获得的幅值差和步骤7中获得的四个加权系数f1、f2、f3和f4,构建解耦矫正函数一用于矫正系统中给定转子电阻构建解耦矫正函数二用于矫正系统中给定励磁电感
具体的,当系统中给定转子电阻和系统中给定励磁电感分别相对于转子电阻Rr和励磁电感Lm产生偏差时,进行以下两种操作:
1)在系统中给定转子电阻中减去经过PI调节器的解耦矫正函数一将其矫正回转子电阻Rr;
2)在系统中给定励磁电感中减去经过PI调节器的解耦矫正函数二将其矫正回励磁电感Lm。
优选地,步骤1中所述的异步电机在静止坐标系αβ下的定子电压向量V的采集方式包括以下两种:
第一种,采样得到实时异步电机线电压Uab、Ucb,经过公式(3)的坐标变换获得定子电压向量V;
第二种,直接采用电机控制器运算单元计算出逆变器调制信号V′代替定子电压向量V。
优选地,步骤1中所述的异步电机在静止坐标系αβ下的定子电流向量i的采集步骤如下:
1)采样得到实时异步电机三相定子电流iA、iB、iC;
2)利用公式(4)的坐标变换获得异步电机在静止坐标系αβ下的定子电流向量i。
优选地,步骤1中所述的异步电机在静止坐标系αβ下的转子电角速度ωr的采集步骤如下:
1)在一个预设采样周期T内采样安装在电机轴上的光电式旋转编码器发出的脉冲数N;
2)根据转子电角速度ωr和光电式旋转编码器发出的脉冲数N以及预设采样周期T之间的关系计算出转子电角速度ωr,其计算公式为:
在公式(5)中,M为光电式旋转编码器旋转一周所产生的脉冲数,P为异步电机极对数,T为预设采样周期。
优选地,步骤2中所述的反电动势e按如下步骤获得:
1)将异步电机的数学模型表示成如下形式:
在公式(6)中Γ=Lr/Lm,λ=LrRs/Lm,为定子电流向量i的微分;
将公式(6)中的反电动势e表示成如下形式:
根据公式(6),将高阶电流滑模观测器设计成如下形式:
2)将公式(8)中的控制项U设计成如下形式:
U=Ueq+Un (9)
在公式(9)中,其中, 为电流观测误差;kp为高阶滑模观测器的比例系数,且满足kp>0;ki为高阶滑模观测器的积分系数,且满足ki>0;
在公式(9)中,Un设计成如下形式:
在公式(10)中,k′为高阶滑模观测器增益,为Un的微分,s为非奇异终端滑模面,
γ=diag(γαγβ),γα>0,γβ>0;
p,q为奇数,且满足1<p/q<2,η10>0,η11>0;
3)当满足时,非奇异终端滑模观测器收敛,反电动势向量e由(11)获得,即:
e=U (11)。
优选地,步骤2中所述的异步电机状态空间表达式(1)按如下步骤获得:
1)将异步电机的数学模型表示成如下形式:
在公式(12)中Γ=Lr/Lm,λ=LrRs/Lm,为i的微分;
在公式(12)中,反电动势向量e按照下式计算:
令转子电角速度ωr的微分则反电动势向量e的微分表示成如下形式:
在公式(14)中,转子磁链电角速度ωs由反电动势e经过锁相环获得;
2)将异步电机数学模型式(12)中的第1行代入公式(14)以消除电流微分项则反电动势向量e的微分进一步表述成如下形式:
3)将公式(15)和公式(12)中第2行联合,构成异步电机状态空间表达式(1)。
优选地,步骤4中所述转子磁链矢量计算值的获得步骤如下:
1)将电流磁链模型表示成如下形式:
其中为转子磁链标称值,为转子磁链标称值的微分;
2)为建立参数偏差时电流磁链模型,将公式(16)中转子时间常数倒数αr替换为系统中给定转子时间常数倒数励磁电感Lm替换成系统中给定励磁电感转子磁链标称值及转子磁链标称值的微分分别替换成转子磁链矢量计算值和转子磁链矢量计算值的微分则转子磁链矢量计算值按如下公式获得:
优选地,步骤5中所述相位差正弦量sin(Δθ)按如下公式获得:
公式(18)中,为转子磁链向量观测值的幅值,为转子磁链矢量计算值的幅值。
优选地,步骤6中幅值差按如下公式获得:
公式(19)中,为转子磁链向量观测值的幅值,为转子磁链矢量计算值的幅值。
优选地,步骤7中加权系数一f1、加权系数二f2、加权系数三f3、加权系数四f4按如下步骤获得:
1)电流磁链模型表示为如下形式:
其中为转子磁链标称值,为转子磁链标称值的微分,定义电感比例系数设电感比例系数Cmr不受磁饱和影响且为常值,电流磁链模型公式(20)简化为:
将公式(21)中转子电阻Rr替换为系统中给定转子电阻励磁电感Lm替换成系统中给定励磁电感转子磁链标称值及转子磁链标称值的微分分别替换成转子磁链矢量计算值和转子磁链矢量计算值的微分获得参数偏差时简化电流磁链模型表示为如下公式:
2)将公式(22)减去公式(21),得到:
其中为转子电阻误差,为励磁电感误差;
3)将公式(23)的相关量表示成如下向量形式:
其中θ1为转子磁链标称值相角,θ2为定子电流相角,j为虚部单位;
将公式(24)带入(23)并设简化获得如下关联函数表达式:
其中ωs1=ωs-ωr为转差频率,为励磁电流给定值,为转矩电流给定值,id为励磁电流实际值,iq为转矩电流实际值;
4)当仅有励磁电感变化时,即ΔRr=0时,公式(25)化简为:
设取加权系数一加权系数二
5)当仅有转子电阻变化时,即ΔLm=0时,公式(25)化简为:
设
取加权系数三加权系数四
本发明相对于现有技术的有益效果是:
1、与基于无功功率和转子磁链q轴分量矫正转子时间常数方案相比,本发明不仅获得了准确的转子时间常数,且实现了励磁电感的在线辨识。
2、与利用转矩矫正方案相比,本发明在全转矩范围内均能实现较好的参数辨识,不存在轻载问题。
3、与通过稳态下磁链模型矫正方案相比,本发明通过转子磁链观测器设计,克服了磁链获得的积分问题的同时,不存在磁链暂态问题。
4、该方案将转子电阻和励磁电感矫正近似解耦,减小了参数矫正之间的相互影响,加快了参数矫正过程。
附图说明
图1为本发明中选取的静止坐标系。
图2为本发明中矫正原理结构框图。
图3为本发明中仅励磁电感变化的矫正效果实验波形。
图4为本发明中仅转子电阻变化的矫正效果实验波形。
图5为本发明中参数同时变化的矫正效果实验波形。
具体实施方式
下面结合附图和实施例对本发明做进一步阐述。
信号采集部分,定子电压向量V和定子电流向量i是通过采样定子线电压Uab、定子线电压Ucb、定子A相电流iA、定子B相电流iB、定子C相电流iC,并经过三相静止坐标系到两相静止坐标系变换获得,实际转子电角速度ωr是利用光电式旋转编码器获得。
图1为本发明中选取的静止坐标系,图2为本发明中矫正原理结构框图。参见图1和图2,本实施例按如下步骤进行。
步骤1,采集异步电机在静止坐标系αβ下的定子电压向量V、定子电流向量i和转子电角速度ωr。
(一)定子电压向量V
定子电压向量V的采集方式包括以下两种:
第一种,采样得到实时异步电机线电压Uab、Ucb,经过公式(3)的坐标变换获得定子电压向量V;
在过程中,首先利用霍尔电压传感器采样定子A、B两相之间的线电压和定子C、B两相之间的线电压,然后输入带有低通滤波器功能的采样调理电路,实现电压信号的采样,定子电压信号Uab和Ucb;在数字信号处理芯片中通过编程对获得的定子电压信号Uab和Ucb进行如式(3)所示的坐标变换,获得定子电压信号的所需形式定子电压向量V。
第二种,直接采用电机控制器运算单元计算出逆变器调制信号V′代替定子电压向量V。
(二)定子电流向量i
1)利用霍尔电流传感器采集定子A相电流iA、定子B相电流iB、定子C相电流iC;
2)然后将其采集的数据输入电流采样通道;在数字信号处理芯片中通过编程对获得的电流采样值iA、iB、iC利用公式(4)的坐标变换获得异步电机在静止坐标系αβ下的定子电流向量i。
(三)转子电角速度ωr
1)在一个预设采样周期T内采样安装在电机轴上的光电式旋转编码器发出的脉冲数N;
2)根据转子电角速度ωr和光电式旋转编码器发出的脉冲数N以及预设采样周期T之间的关系计算出转子电角速度ωr,其计算公式为:
在公式(5)中,M为光电式旋转编码器旋转一周所产生的脉冲数,P为异步电机极对数,T为预设采样周期。
步骤2,建立在静止坐标系αβ下异步电机状态空间表达式为:
在公式(1)中,e为反电动势向量,为e的微分,为转子磁链向量,为的微分,ωs为转子磁链电角速度,Lm为励磁电感,Lr为转子电感,Ls为定子电感,αr=Rr/Lr为转子时间常数倒数,ε=δLsLr/Lm,其中Rr为转子电阻,Rs为定子电阻, 为反对称矩阵。
(一)所述的反电动势e按如下步骤获得:
1)将异步电机的数学模型表示成如下形式:
在公式(6)中Γ=Lr/Lm,λ=LrRs/Lm,为定子电流向量i的微分;
将公式(6)中的反电动势e表示成如下形式:
根据公式(6),将高阶电流滑模观测器设计成如下形式:
2)将公式(8)中的控制项U设计成如下形式:
U=Ueq+Un (9)
在公式(9)中,其中, 为电流观测误差;kp为高阶滑模观测器的比例系数,且满足kp>0;ki为高阶滑模观测器的积分系数,且满足ki>0;
在公式(9)中,Un设计成如下形式:
在公式(10)中,k′为高阶滑模观测器增益,为Un的微分,s为非奇异终端滑模面,
γ=diag(γαγβ),γα>0,γβ>0;
p,q为奇数,且满足1<p/q<2,η10>0,η11>0;
3)当满足时,非奇异终端滑模观测器收敛,反电动势向量e由(11)获得,即:
e=U (11)。
(二)所述的异步电机状态空间表达式(1)按如下步骤获得:
1)将异步电机的数学模型表示成如下形式:
在公式(12)中Γ=Lr/Lm,λ=LrRs/Lm,为i的微分;
在公式(12)中,反电动势向量e按照下式计算:
虑到转速等机械量的动态变化相比于定子电流和转子磁链等电气量的变化较慢,在电气时间常数范围内,令转子电角速度ωr的微分则反电动势向量e的微分表示成如下形式:
在公式(14)中,转子磁链电角速度ωs由反电动势e经过锁相环获得;
2)将异步电机数学模型式(12)中的第1行代入公式(14)以消除电流微分项则反电动势向量e的微分进一步表述成如下形式:
3)将公式(15)和公式(12)中第2行联合,构成异步电机状态空间表达式(1)。
步骤3,根据公式(2)建立转子磁链滑模观测器,对转子磁链向量进行观测:
在公式(2)中,为反电动势向量e的观测值,为的微分,为转子磁链向量观测值,为的微分, 为反电动势误差,为系统中给定励磁电感,为系统中给定转子电感,为系统中给定定子电感,为系统中给定转子时间常数倒数,其中为系统中给定转子电阻,为系统中给定定子电阻,k为设定的观测器增益一,g为设定的观测器增益二,k和g的数值均为负数。
步骤4,根据步骤1中得到定子电压向量V、定子电流向量i和转子电角速度ωr,经过参数偏差时电流磁链模型计算获得转子磁链矢量计算值具体步骤如下;
1)将电流磁链模型表示成如下形式:
其中为转子磁链标称值,为转子磁链标称值的微分;
2)为建立参数偏差时电流磁链模型,将公式(16)中转子时间常数倒数αr替换为系统中给定转子时间常数倒数励磁电感Lm替换成系统中给定励磁电感转子磁链标称值及转子磁链标称值的微分分别替换成转子磁链矢量计算值和转子磁链矢量计算值的微分则转子磁链矢量计算值按如下公式获得:
步骤5,根据步骤3中获得的转子磁链向量观测值和步骤4中获得的转子磁链矢量计算值经过叉乘获得相位差正弦量sin(Δθ),并设相位差Δθ≈sin(Δθ),其中所述相位差 为转子磁链矢量计算值的相位,θ为转子磁链向量观测值的相位。
所述相位差正弦量sin(Δθ)按如下公式获得:
公式(18)中,为转子磁链向量观测值的幅值,为转子磁链矢量计算值的幅值,且sin(Δθ)≈Δθ。
步骤6,根据步骤3中获得的转子磁链向量观测值和步骤4中获得的转子磁链矢量计算值两者幅值作差获得幅值差其计算公式为;
公式(19)中,为转子磁链向量观测值的幅值,为转子磁链矢量计算值的幅值。
步骤7,利用电流磁链模型推导获得加权系数一f1、加权系数二f2、加权系数三f3、加权系数四f4。具体的,按照如下步骤获得:
1)电流磁链模型表示为如下形式:
其中为转子磁链标称值,为转子磁链标称值的微分,定义电感比例系数设电感比例系数Cmr不受磁饱和影响且为常值,电流磁链模型公式(20)简化为:
将公式(21)中转子电阻Rr替换为系统中给定转子电阻励磁电感Lm替换成系统中给定励磁电感转子磁链标称值及转子磁链标称值的微分分别替换成转子磁链矢量计算值和转子磁链矢量计算值的微分获得参数偏差时简化电流磁链模型表示为如下公式:
2)将公式(22)减去公式(21),得到:
其中为转子电阻误差,为励磁电感误差;
3)将(23)相关量表示成如下向量形式:
其中θ1为转子磁链标称值相角,θ2为定子电流相角,j为虚部单位;
将公式(24)带入(23)并设简化获得如下关联函数表达式:
其中ωs1=ωs-ωr为转差频率,为励磁电流给定值,为转矩电流给定值,id为励磁电流实际值,iq为转矩电流实际值;
4)当仅有励磁电感变化时,即ΔRr=0时,公式(25)化简为:
设取加权系数一加权系数二
5)当仅有转子电阻变化时,即ΔLm=0时,公式(25)化简为:
设
取加权系数三加权系数四
步骤8,根据步骤5中获得的相位差Δθ、步骤6中获得的幅值差和步骤7中获得的四个加权系数f1、f2、f3和f4,构建解耦矫正函数一用于矫正系统中给定转子电阻构建解耦矫正函数二用于矫正系统中给定励磁电感
具体的,当系统中给定转子电阻和系统中给定励磁电感分别相对于转子电阻Rr和励磁电感Lm产生偏差时,进行以下两种操作:
1)在系统中给定转子电阻中减去经过PI调节器的解耦矫正函数一将其矫正回转子电阻Rr;
2)在系统中给定励磁电感中减去经过PI调节器的解耦矫正函数二将其矫正回励磁电感Lm。
图3-5给出了本实施例在转子电阻和励磁电感辨识过程中的性能表现。
实验电机参数:额定功率pN=15kW,额定电压UN=180V,转子电阻Rr=0.0122Ω,定子电阻Rs=0.027Ω,励磁电感Lm=2.18mΗ,转子电感Lr=2.271mΗ,定子电感Ls=2.271mΗ,极对数P=2,额定频率fN=120HZ。实验中观测器参数的给定值为:kp=0.1,ki=100,γα=0.001,γβ=0.001,p=13,q=11,η10=100,η11=100;k′=3448,k=-1000,g=-0.1。
在试验中为了更好地研究参数校正方案的性能,避免电机磁饱和因素的结果分析的影响,实验时对电机进行了弱磁运行,让电机工作在磁路的线性区域。给定电角速度10Hz,给定励磁电流35A,转矩电流35A。图3中仅励磁电感变化30%,可以发现启用矫正后,励磁电感能较快地矫正回真实值,转子电阻受影响很小;图4中仅转子电阻变化50%,同样,启用矫正后转子电阻能较快回到真实在,励磁电感受影响较小;图5转子电阻和励磁电感同时变化,该算法能较快地将两个参数矫正回真实值。
Claims (10)
1.一种异步电机转子电阻及励磁电感解耦矫正方法,其特征在于,包括如下步骤:
步骤1,采集异步电机在静止坐标系αβ下的定子电压向量V、定子电流向量i和转子电角速度ωr;
步骤2,建立在静止坐标系αβ下异步电机状态空间表达式为:
在公式(1)中,e为反电动势向量,为e的微分,为转子磁链向量,为的微分,ωs为转子磁链电角速度,Lm为励磁电感,Lr为转子电感,Ls为定子电感,αr=Rr/Lr为转子时间常数倒数,ε=δLsLr/Lm,其中Rr为转子电阻,Rs为定子电阻, 为反对称矩阵;
步骤3,根据公式(2)建立转子磁链滑模观测器,对转子磁链向量进行观测:
在公式(2)中,为反电动势向量e的观测值,为的微分,为转子磁链向量观测值,为的微分, 为反电动势误差,为系统中给定励磁电感,为系统中给定转子电感,为系统中给定定子电感,为系统中给定转子时间常数倒数,其中为系统中给定转子电阻,为系统中给定定子电阻,k为设定的观测器增益一,g为设定的观测器增益二,k和g的数值均为负数;
步骤4,根据步骤1中得到定子电压向量V、定子电流向量i和转子电角速度ωr,经过参数偏差时电流磁链模型计算获得转子磁链矢量计算值
步骤5,根据步骤3中获得的转子磁链向量观测值和步骤4中获得的转子磁链矢量计算值经过叉乘获得相位差正弦量sin(Δθ),并设相位差Δθ≈sin(Δθ),其中所述相位差 为转子磁链矢量计算值的相位,θ为转子磁链向量观测值的相位;
步骤6,根据步骤3中获得的转子磁链向量观测值和步骤4中获得的转子磁链矢量计算值两者幅值作差获得幅值差
步骤7,利用电流磁链模型推导获得加权系数一f1、加权系数二f2、加权系数三f3、加权系数四f4;
步骤8,根据步骤5中获得的相位差Δθ、步骤6中获得的幅值差和步骤7中获得的四个加权系数f1、f2、f3和f4,构建解耦矫正函数一用于矫正系统中给定转子电阻构建解耦矫正函数二用于矫正系统中给定励磁电感
具体的,当系统中给定转子电阻和系统中给定励磁电感分别相对于转子电阻Rr和励磁电感Lm产生偏差时,进行以下两种操作:
1)在系统中给定转子电阻中减去经过PI调节器的解耦矫正函数一将其矫正回转子电阻Rr;
2)在系统中给定励磁电感中减去经过PI调节器的解耦矫正函数二将其矫正回励磁电感Lm。
2.根据权利要求1所述的一种异步电机转子电阻及励磁电感解耦矫正方法,其特征在于,步骤1中所述的异步电机在静止坐标系αβ下的定子电压向量V的采集方式包括以下两种:
第一种,采样得到实时异步电机线电压Uab、Ucb,经过公式(3)的坐标变换获得定子电压向量V;
第二种,直接采用电机控制器运算单元计算出逆变器调制信号V′代替定子电压向量V。
3.根据权利要求1所述的一种异步电机转子电阻及励磁电感解耦矫正方法,其特征在于,步骤1中所述的异步电机在静止坐标系αβ下的定子电流向量i的采集步骤如下:
1)采样得到实时异步电机三相定子电流iA、iB、iC;
2)利用公式(4)的坐标变换获得异步电机在静止坐标系αβ下的定子电流向量i。
4.根据权利要求1所述的一种异步电机转子电阻及励磁电感解耦矫正方法,其特征在于,步骤1中所述的异步电机在静止坐标系αβ下的转子电角速度ωr的采集步骤如下:
1)在一个预设采样周期T内采样安装在电机轴上的光电式旋转编码器发出的脉冲数N;
2)根据转子电角速度ωr和光电式旋转编码器发出的脉冲数N以及预设采样周期T之间的关系计算出转子电角速度ωr,其计算公式为:
在公式(5)中,M为光电式旋转编码器旋转一周所产生的脉冲数,P为异步电机极对数,T为预设采样周期。
5.根据权利要求1所述的一种异步电机转子电阻及励磁电感解耦矫正方法,其特征在于,步骤2中所述的反电动势e按如下步骤获得:
1)将异步电机的数学模型表示成如下形式:
在公式(6)中Γ=Lr/Lm,λ=LrRs/Lm,为定子电流向量i的微分;
将公式(6)中的反电动势e表示成如下形式:
根据公式(6),将高阶电流滑模观测器设计成如下形式:
2)将公式(8)中的控制项U设计成如下形式:
U=Ueq+Un (9)
在公式(9)中,其中, 为电流观测误差;kp为高阶滑模观测器的比例系数,且满足kp>0;ki为高阶滑模观测器的积分系数,且满足ki>0;
在公式(9)中,Un设计成如下形式:
在公式(10)中,k′为高阶滑模观测器增益,为Un的微分,s为非奇异终端滑模面,
γ=diag(γα γβ),γα>0,γβ>0;
p,q为奇数,且满足1<p/q<2,η10>0,η11>0;
3)当满足时,非奇异终端滑模观测器收敛,反电动势向量e由(11)获得,即:
e=U (11)。
6.根据权利要求1所述的一种异步电机转子电阻及励磁电感解耦矫正方法,其特征在于,步骤2中所述的异步电机状态空间表达式(1)按如下步骤获得:
1)将异步电机的数学模型表示成如下形式:
在公式(12)中Γ=Lr/Lm,λ=LrRs/Lm,为i的微分;
在公式(12)中,反电动势向量e按照下式计算:
令转子电角速度ωr的微分则反电动势向量e的微分表示成如下形式:
在公式(14)中,转子磁链电角速度ωs由反电动势e经过锁相环获得;
2)将异步电机数学模型式(12)中的第1行代入公式(14)以消除电流微分项则反电动势向量e的微分进一步表述成如下形式:
3)将公式(15)和公式(12)中第2行联合,构成异步电机状态空间表达式(1)。
7.根据权利要求1所述的一种异步电机转子电阻及励磁电感解耦矫正方法,其特征在于,步骤4中所述转子磁链矢量计算值的获得步骤如下:
1)将电流磁链模型表示成如下形式:
其中为转子磁链标称值,为转子磁链标称值的微分;
2)为建立参数偏差时电流磁链模型,将公式(16)中转子时间常数倒数αr替换为系统中给定转子时间常数倒数励磁电感Lm替换成系统中给定励磁电感转子磁链标称值及转子磁链标称值的微分分别替换成转子磁链矢量计算值和转子磁链矢量计算值的微分则转子磁链矢量计算值按如下公式获得:
8.根据权利要求1所述的一种异步电机转子电阻及励磁电感解耦矫正方法,其特征在于,步骤5中所述相位差正弦量sin(Δθ)按如下公式获得:
公式(18)中,为转子磁链向量观测值的幅值,为转子磁链矢量计算值的幅值。
9.根据权利要求1所述的一种异步电机转子电阻及励磁电感解耦矫正方法,其特征在于,步骤6中幅值差按如下公式获得:
公式(19)中,为转子磁链向量观测值的幅值,为转子磁链矢量计算值的幅值。
10.根据权利要求1所述的一种异步电机转子电阻及励磁电感解耦矫正方法,其特征在于,步骤7中加权系数一f1、加权系数二f2、加权系数三f3、加权系数四f4按如下步骤获得:
1)电流磁链模型表示为如下形式:
其中为转子磁链标称值,为转子磁链标称值的微分,定义电感比例系数设电感比例系数Cmr不受磁饱和影响且为常值,电流磁链模型公式(20)简化为:
将公式(21)中转子电阻Rr替换为系统中给定转子电阻励磁电感Lm替换成系统中给定励磁电感转子磁链标称值及转子磁链标称值的微分分别替换成转子磁链矢量计算值和转子磁链矢量计算值的微分获得参数偏差时简化电流磁链模型表示为如下公式:
2)将公式(22)减去公式(21),得到:
其中为转子电阻误差,为励磁电感误差;
3)将公式(23)的相关量表示成如下向量形式:
其中θ1为转子磁链标称值相角,θ2为定子电流相角,j为虚部单位;
将公式(24)带入(23)并设简化获得如下关联函数表达式:
其中ωs1=ωs-ωr为转差频率,为励磁电流给定值,为转矩电流给定值,id为励磁电流实际值,iq为转矩电流实际值;
4)当仅有励磁电感变化时,即ΔRr=0时,公式(25)化简为:
设取加权系数一加权系数二
5)当仅有转子电阻变化时,即ΔLm=0时,公式(25)化简为:
设
取加权系数三加权系数四
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611201979.4A CN106685294B (zh) | 2016-12-23 | 2016-12-23 | 异步电机转子电阻及励磁电感解耦矫正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611201979.4A CN106685294B (zh) | 2016-12-23 | 2016-12-23 | 异步电机转子电阻及励磁电感解耦矫正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106685294A true CN106685294A (zh) | 2017-05-17 |
CN106685294B CN106685294B (zh) | 2019-05-07 |
Family
ID=58870158
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611201979.4A Active CN106685294B (zh) | 2016-12-23 | 2016-12-23 | 异步电机转子电阻及励磁电感解耦矫正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106685294B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113696736A (zh) * | 2021-08-06 | 2021-11-26 | 上汽通用五菱汽车股份有限公司 | 车辆旋变失效保护方法、装置及计算机可读存储介质 |
CN114389497A (zh) * | 2022-03-24 | 2022-04-22 | 希望森兰科技股份有限公司 | 一种异步电机电压电流混合磁链观测器定向误差补偿方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102340278A (zh) * | 2011-09-30 | 2012-02-01 | 哈尔滨工业大学(威海) | 矢量变频器中电机定子磁链的估算方法 |
CN103326659A (zh) * | 2013-07-05 | 2013-09-25 | 合肥工业大学 | 一种异步电机转子磁场定向校正的控制方法 |
US9024569B2 (en) * | 2013-06-28 | 2015-05-05 | Eaton Corporation | System and method of rotor time constant online identification in an AC induction machine |
CN105281630A (zh) * | 2015-11-08 | 2016-01-27 | 浙江大学 | 异步电机无速度传感器系统中在线辨识定转子电阻的方法 |
CN105610369A (zh) * | 2016-03-04 | 2016-05-25 | 合肥工业大学 | 一种基于滑模观测器的异步电机磁链观测方法 |
-
2016
- 2016-12-23 CN CN201611201979.4A patent/CN106685294B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102340278A (zh) * | 2011-09-30 | 2012-02-01 | 哈尔滨工业大学(威海) | 矢量变频器中电机定子磁链的估算方法 |
US9024569B2 (en) * | 2013-06-28 | 2015-05-05 | Eaton Corporation | System and method of rotor time constant online identification in an AC induction machine |
CN103326659A (zh) * | 2013-07-05 | 2013-09-25 | 合肥工业大学 | 一种异步电机转子磁场定向校正的控制方法 |
CN105281630A (zh) * | 2015-11-08 | 2016-01-27 | 浙江大学 | 异步电机无速度传感器系统中在线辨识定转子电阻的方法 |
CN105610369A (zh) * | 2016-03-04 | 2016-05-25 | 合肥工业大学 | 一种基于滑模观测器的异步电机磁链观测方法 |
Non-Patent Citations (1)
Title |
---|
杨淑英等: ""基于磁链观测器的异步电机转矩闭环控制"", 《太阳能学报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113696736A (zh) * | 2021-08-06 | 2021-11-26 | 上汽通用五菱汽车股份有限公司 | 车辆旋变失效保护方法、装置及计算机可读存储介质 |
CN114389497A (zh) * | 2022-03-24 | 2022-04-22 | 希望森兰科技股份有限公司 | 一种异步电机电压电流混合磁链观测器定向误差补偿方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106685294B (zh) | 2019-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110198150B (zh) | 一种永磁同步电机多参数在线辨识方法 | |
CN108551287B (zh) | 车用内置式永磁同步电机驱动系统转矩闭环控制方法 | |
CN110350835B (zh) | 一种永磁同步电机无位置传感器控制方法 | |
CN105763119B (zh) | 一种cmg框架永磁同步电机的控制系统及控制方法 | |
CN103178769B (zh) | 永磁同步电机静止状态下参数离线辨识方法 | |
CN107482982B (zh) | 一种基于铁损模型的异步电机矢量控制方法 | |
CN105610369B (zh) | 一种基于滑模观测器的异步电机磁链观测方法 | |
CN103326659B (zh) | 一种异步电机转子磁场定向校正的控制方法 | |
CN108900129A (zh) | 一种无位置传感器永磁同步电机控制方法、装置和系统 | |
CN101286724A (zh) | 一种永磁同步电机直接转矩控制方法 | |
JP2003513594A (ja) | 不均衡状態に起因するモータトルクリップルの最小化 | |
CN103401503A (zh) | 一种在谐波平面在线辨识双三相电机参数的方法及装置 | |
CN106452256B (zh) | 基于转子磁链观测器的异步电机参数在线矫正方法 | |
CN113691186A (zh) | 一种永磁同步电机无位置传感器控制转子位置角补偿方法 | |
CN105680752B (zh) | 一种永磁同步电机饱和电感参数的辨识方法及系统 | |
CN113131817A (zh) | 一种永磁同步电机在线参数辨识系统及方法 | |
CN101557193B (zh) | 非正弦反电动势表面式交流永磁电机矢量控制方法 | |
CN111769779A (zh) | 基于改进型Luenberger观测器的PMSM直接转矩控制方法 | |
CN115864928A (zh) | 一种基于校正电流预测的pmsm模型参考自适应转速估算方法 | |
CN106685294B (zh) | 异步电机转子电阻及励磁电感解耦矫正方法 | |
CN117081448B (zh) | 一种内置式永磁同步电机高精度动态建模方法 | |
Yang et al. | Stability analysis of q-axis rotor flux based model reference adaptive system updating rotor time constant in induction motor drives | |
Horch et al. | Nonlinear integral backstepping control for induction motor drive with adaptive speed observer using super twisting strategy | |
Guziński et al. | Sensorless low speed PMSM motor control with cogging torque compensation | |
CN114567227A (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 |