CN100424989C - 用于异步旋转电气机械的一种直进式控制方法和控制单元和一种用于其的存储介质 - Google Patents

用于异步旋转电气机械的一种直进式控制方法和控制单元和一种用于其的存储介质 Download PDF

Info

Publication number
CN100424989C
CN100424989C CNB2006100848514A CN200610084851A CN100424989C CN 100424989 C CN100424989 C CN 100424989C CN B2006100848514 A CNB2006100848514 A CN B2006100848514A CN 200610084851 A CN200610084851 A CN 200610084851A CN 100424989 C CN100424989 C CN 100424989C
Authority
CN
China
Prior art keywords
centerdot
rightarrow
point
matrix
vector
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.)
Expired - Fee Related
Application number
CNB2006100848514A
Other languages
English (en)
Other versions
CN1870419A (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.)
Alstom Transport Technologies SAS
Original Assignee
Alstom Transport SA
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 Alstom Transport SA filed Critical Alstom Transport SA
Publication of CN1870419A publication Critical patent/CN1870419A/zh
Application granted granted Critical
Publication of CN100424989C publication Critical patent/CN100424989C/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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 Ac Motors In General (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

控制异步旋转机械的这种直进式控制方法包括:a)建立设置点电流的步骤(60),该设置点电流适于在间隔T的末端处实现机械力设置点和设置点磁通,并且基于在所述设置点电流、设置点磁通及在间隔T的开始处的初始电流和磁通值之间的关系而被建立;和b)计算在间隔T期间连续施加的平均控制向量的步骤(64),以在间隔T的末端处得到与在步骤a)期间建立的设置点电流相等的电流。

Description

用于异步旋转电气机械的一种直进式控制方法和控制单元和一种用于其的存储介质
技术领域
本发明涉及用于异步旋转电气机械的一种直进式(deadbeat)控制方法和一种直进式控制单元,并且涉及一种用在该直进式控制方法中的存储介质。
背景技术
在说明书的其余部分中,术语“电机”指异步旋转电气机械。通过下面给出的对于说明书的前序部分,这里应该指出,异步旋转电气机械是各向同性电机,即其中直接电感等于正交电感的电机。在各向同性电机的情况下,通过把矩阵和向量写在复数平面中使下面所定义的维数mn较小。
直进式控制方法可使用在现有技术中已知的状态方程形式,该状态方程形式基于原理:电机的状态完全由代表其自由度的参数所取的已知值描述。电机的状态可由n维状态向量
Figure C20061008485100051
特征化,其中n等于自由度的数量。由适于控制机械的m维控制向量
Figure C20061008485100052
代表的控制输入生成的在机械中的变化那么由如下线性状态方程系统描述,这对于在自动控制领域的技术人员是熟悉的:
X → · = A · X → + B · V → - - - ( 1 )
其中:
·
Figure C20061008485100054
是状态向量
Figure C20061008485100055
相对于时间的导数;
·
Figure C20061008485100056
是瞬时控制向量;
·A是在没有任何控制输入的情况下机械的自由行为的n.n维矩阵;及
·B是n.m维控制矩阵。
在各向同性电机的情况下,通过把矩阵和向量写在复数平面中使下面所定义的尺寸mn较小。
两个矩阵A和B代表电机的线性模型,并且从对于给定动态的电机的电气微分方程得到。不是线性的模型必须在操作点周围线性化,并且在这种情况下,因此必须具有适用的多个模型。
为了确定电机在平均控制向量
Figure C20061008485100061
的连续应用间隔T的末端处的状态,必须通过间隔T在时间中间隔开的两个时刻tn与tn+1之间积分方程(1)。结果可以放入到在现有持术中已知的状态方程的离散系统的形式中:
X → ( t n + 1 ) = F ( T ) · X → ( t n ) + G ( T ) · V → ‾ ( t n → n + 1 ) - - - ( 2 )
其中:
·F(T)是由F=eA·T定义的电机的n.n维转换矩阵;
·G(T)是由G=A-1·(eA·T-Inn)·B定义的n.m维控制矩阵,其中Inn是n.n维单位矩阵;
·
Figure C20061008485100064
是分别在时刻tn+1和tn处的状态向量;及
·
Figure C20061008485100065
是在间隔T期间,即从时刻tn至时刻tn+1,施加的平均控制向量。
平均控制向量
Figure C20061008485100066
或者是在间隔T期间连续施加的瞬时向量,或者是在间隔T期间直接施加到电机上的瞬时控制向量的时间延续的平均值。在瞬时向量的时间延续的情况下,每个瞬时向量施加与电机的时间常数相比非常短的时间段,结果是这种瞬时向量的连续施加产生与在相同时间间隔T期间连续施加一个选定相位和幅值的单个瞬时向量的相同的效果。
常常使用连续的瞬时向量,因为大多数电机驱动器只能产生有限数量个振幅和相位的瞬时控制向量。例如,三相转换器只能产生六个不同的非零瞬时控制向量。这是因为,为了能够从三相转换器得到任意振幅和相位的控制向量,标准的实践是,把时间连续的瞬时控制向量直接施加到电机上,该瞬时控制向量在时刻tn与tn+1之间的平均值等于一个其相位和振幅随意选择的平均控制向量。例如,脉冲宽度调制方法产生一个其相位和振幅可以随意从电源装置选择的平均控制向量,只能够产生有限数量个振幅和相位的瞬时控制向量。
在描述的其余部分中,除非另有说明,术语“控制向量”指平均控制向量。
控制向量典型地是把特定电压同时施加到电机的所有相上的电压向量。
如果例如基于诸如Kalman滤波器之类的观察器可以测量或估计的电机的模型(A,B)和初始状态 X → 0 = X → ( t n ) 是已知的,则变得有可能预测在间隔T期间控制向量
Figure C20061008485100072
的连续施加之后的新状态 X → P = X → ( t n + T ) .
相反,有可能计算谋求预测状态与设置点状态一致的控制向量
Figure C20061008485100074
这相当于用设置点状态
Figure C20061008485100075
替换
Figure C20061008485100076
状态方程的离散系统然后成为:
X → c = F ( T ) · X → 0 + G ( T ) · V → ‾ - - - ( 3 )
然而,矩阵G(T)不是方矩阵,并且因此不能求逆,结果是我们不知道如何用分析方法计算是方程系统(3)的解的控制向量
Figure C20061008485100078
情况是,方程系统按常规求逆,并且借助于近似方法计算控制向量。在其中磁通变化缓慢的旋转电机或旋转电气机械的情况下,近似包括假定电机的磁通在稳定状态条件下建立。
借助于这个种类的近似,然后有可能计算能够在间隔T的末端处实现设置点转矩的控制向量。欧洲专利申请EP-A-1 045 514公开了以上种类的直进式控制方法的一个例子。
然而,当在实际中使用这种直进式控制方法时,电机被逐渐退磁。这些控制方法因此具有如下缺陷:必须添加环路使电机的磁通跟随设置点磁通,以防止电机通过直进式控制方法退磁。
发明内容
本发明目的在于,通过提出一种不必使用环路使电机的磁通随动的直进式控制方法来解决这个问题。
因而本发明包括一种控制电机的直进式控制方法,该电机的磁通取决于其激励电流,该方法包括:
a)建立设置点电流的步骤,该设置点电流适于在间隔T的末端处实现机械力设置点和设置点磁通,并且基于在设置点电流、设置点磁通及在间隔T的开始处的初始电流和磁通值之间的关系而被建立;和
b)计算在间隔T期间连续施加的平均控制向量的步骤,以在间隔T的末端处得到与在步骤a)期间建立的设置点电流相等的电流。
与现有技术直进式控制方法相反,以上方法在计算用来实现设置点电流的控制向量之前,首先建立用来实现机械力设置点和设置点磁通的电流设置点。电机的磁化因此被保留,而不必添加用于这个目的的磁通控制环路,因为磁通保持等于设置点磁通。
这种控制方法的实施例可以具有如下特征的一个或多个:
·设置点电流是在复数平面中如下关系的解:
( a 2 · p 11 - a 1 · p 21 ) · I → c + ( a 2 · p 12 - a 1 · p 22 ) · Φ → c =
( a 2 · e 1 · p 11 - a 1 · e 2 · p 21 ) · I → 0 + ( a 2 · e 1 · p 12 - a 1 · e 2 · p 22 ) · Φ → 0
其中:
·pij是由如下关系A=P·D·P-1定义的特征值的矩阵P的逆矩阵的系数,其中D是包括矩阵A的特征值μi的对角矩阵,P-1是矩阵P的逆矩阵,及矩阵A是在模型化电机操作的状态方程的线性系统中电机的自由特性的矩阵;
·系数ai是由如下关系C=D-1·(eD·T-I)·(P-1·B)定义的矩阵C的系数,其中矩阵I是单位矩阵,并且矩阵B是在状态方程的线性系统中的电机的控制矩阵;
·向量
Figure C20061008485100083
Figure C20061008485100084
分别代表在间隔T的开始处电机的电流的和磁通的初始状态;
·向量
Figure C20061008485100085
Figure C20061008485100086
分别是设置点电流向量和设置点磁通向量;及
·系数ei由如下关系 e i = e μ i · T 定义;
·控制向量是电机的状态方程的如下离散系统的准确解:
X → c = F ( T ) · X → 0 + G ( T ) · V → ‾
其中:
·向量
Figure C20061008485100091
是定义在间隔T的开始处的电机的电流和磁通的初始状态的状态向量;
·向量
Figure C20061008485100092
是定义电机的设置点电流和设置点磁通的设置点状态向量;
·向量
Figure C20061008485100093
是电机的控  制向量;
·矩阵F(T)是由关系F=eA·T定义的电机的转换矩阵,其中A是在模型化电机操作的状态方程的线性系统中的电机的自由特性的矩阵;及
·矩阵G(T)是电机的控制矩阵,该控制矩阵的值是间隔T的值的函数;
·控制向量是如下关系的结果:
Figure C20061008485100094
其中:
·θ0是初始轴线系统的角位置,该系统的横坐标轴在间隔T的开始处与电机的初始磁通向量对准;
·系数pij是由关系A=P·D·P-1定义的特征值的矩阵P的逆矩阵的系数,其中D是包括矩阵A的特征值μi的对角矩阵,P-1是矩阵P的逆矩阵,及矩阵A是在模型化电机操作的状态方程的线性系统中电机的自由特性的矩阵;
Figure C20061008485100095
是在初始轴线系统中在间隔T的开始处初始电流向量的坐标,该系统的横坐标轴与初始磁通向量的方向对准;
Figure C20061008485100097
是在初始轴线系统中初始磁通向量的横坐标;
·Icd、Icq是在设置点轴线系统中在间隔T的末端处设置点电流向量的坐标,该系统的横坐标轴与设置点磁通向量的方向对准;及
·Φcd是在设置点轴线系统中设置点磁通向量的横坐标;
·模型化电机操作的状态方程的线性系统的状态向量由定子电流向量和转子磁通向量形成。
这种方法的实施例也具有如下优点:
·设置点电流是以上方程的解的事实,保证它将在间隔T的末端处实现机械力设置点和要求的设置点磁通;
·控制向量也是以上状态方程系统的准确解的事实,也保证控制方法将在间隔T的末端处实现已经固定的设置点电流;
·容易测量定子电流向量,这便于控制方法的实施例;及
·转子磁通向量变化缓慢,这便于借助于观察器估计它。
本发明也提供一种用来控制在其激励电流的基础上受控的旋转电机的直进式控制单元,所述单元包括:
·用来建立设置点电流的模块,该设置点电流适于在间隔T的末端处实现机械力设置点和设置点磁通,并且基于在设置点电流、设置点磁通及在间隔T的开始处的初始电流和磁通值之间的关系而被建立;和
·用来计算在间隔T期间连续施加的控制向量的模块,以在间隔T的末端处得到与用来建立的设置点电流的模块建立的设置点电流相等的电流。
本发明还提供一种包含指令的信息存储介质,用来当由电子计算机执行所述指令时执行以上控制方法。
附图说明
在阅读仅通过例子给出的如下描述和参考附图时,可更清楚地理解本发明,在附图中:
图1是异步电动机的等效电路;
图2A和2B分别限定在直进式控制方法中使用的初始轴线系统和设置点轴线系统;
图3是三相异步电动机控制系统的构造图;及
图4是在图3系统中使用的直进式控制方法的流程图。
具体实施方式
在参照图3和4详细描述直进式控制系统和方法的具体实施例(见部分III)之前,在各向同性旋转电气机械的情形中描述在系统和方法的开发背后的数学理论(见部分I),并且然后在异步电动机的特定情况下更详细地描述在该说明书中开发的一般方程(见部分II)。
在该说明书的引言部分中引入的符号也用在说明书的其余部分中。
部分I):直进式控制方程-一般情形
I.1-在对于机械专用的向量基础中的状态方程系统的表达
为了计算各向同性旋转电气机械的控制向量,假定系统(3)适于自动控制。
作为一般规则,两个矩阵F(T)和G(T)形式上是未知的。有可能以各种方式确定它们。下面基于系统的计算特征值:μj,j∈{1,n}描述的计算方法表明简单和准确的分析解,因为系统的模式是足够熟知的。
让D是A的对角矩阵,A的对角线包括n个特征值μi;并且P是特征值的矩阵(通过矩阵),从而:
A=P·D·P-1    (4)
转换矩阵那么是简单的以计算:
F=eA·T=P·eD·T·P-1    (5)
其中eD·T是包括特征值μi乘以间隔T的指数的对角矩阵。使用这种方法,离散状态方程系统(3)写成:
X → c = P · e D · T · P - 1 · X → 0 + A - 1 · ( P · e D · T · P - 1 - I ) · B · V → ‾ - - - ( 6 )
通过把两个左手项乘以P-1,在几次变换后得到:
[ P - 1 · X → c ] = e D · T · [ P - 1 · X → 0 ] + D - 1 · ( e D · T - I ) · [ P - 1 · B ] · V → ‾ - - - ( 7 )
机械的这种新分立状态表示因为其简单性是有益的:eD·T和D-1·(eD·T-I)是用于给定机械的对角矩阵,并且非常容易由特征值μi形式计算。事实上,它们的相应对角分别由如下元素组成:
Figure C20061008485100113
Figure C20061008485100114
对于j∈{1,...,n}    (8)
这个公式表示状态向量到特征向量空间,即状态“特征向量” Ψ → = P - 1 · X → , 以及控制“特征矩阵”:[P-1·B]的空间中的投影。作为特征化系统的特性的特征值,这些状态特征向量具有特定的物理意义,因为它们的坐标,即“特征值”,与控制定律分离,因为它们表示在特征向量基中。这同样反映诸个矩阵的两个是对角的事实。
变量的如下变化简化状态方程(7):
Ψ → c = P - 1 · X → c - - - ( 8 )
Ψ → 0 = P - 1 · X → 0 - - - ( 9 )
Figure C20061008485100123
Ψ → c = e D · T · Ψ → 0 + C · V → ‾ - - - ( 11 )
控制向量的计算然后包括:对于不可得到的参数填入来自测量和/或观察器的初始特征状态
Figure C20061008485100125
填入设置点特征状态
Figure C20061008485100126
作为要实现的目标的函数;及最后计算能够在选择时间间隔T的末端处实现设置点特征状态的控制向量。
对于旋转电气机械,由于目标是得到机械力,所以实际上是,选择由机械的激励电流的状态和机械的磁通的状态组成的状态向量
Figure C20061008485100127
事实上,通过简单的向量乘积,在电流向量和磁通向量之间的幅值和角度把力目标变换成关于在选择的轴线系统中的那些向量的坐标,以简化其表达式。这里轴线系统是其横坐标轴线与磁通向量对准的正交轴线系统。在这种情况下,Laplace定律向量乘积由磁通的模量和在相同时刻电流向量到纵坐标轴线上的投影的模量的乘积表示,忽略取决于机械的常数。况且,磁通的模量和在相同时刻到横坐标轴线上的电流向量的投影的模量的乘积定义电流和磁通向量的标量乘积,并因此定义在Park表示中的机械的磁能,忽略取决于机械的乘法常数。
矩阵Γ0定义为是在初始轴线系统中电流向量的投影矩阵,该系统的横坐标轴与初始磁通向量对准。这里,“初始”是指标记间隔T开始的时刻tn。这个矩阵能以两个矩阵Γc和ΓT乘积的形式写出,因为在第一个情况下,把电流向量投影到其横坐标轴与设置点磁通向量对准的设置点轴线系统中,并且在第二个情况下,把该轴线系统的电流向量向初始轴线系统投影。这种变换ΓT当然取决于间隔T的值,由于它从设置点轴线系统到初始轴线系统(假定系统对于这种形式的表达是因果的)。因此满足如下关系:
ΓT·Γc=Γ0    (12)
把关系(11)和(12)结合产生:
Γ T · [ Γ c · Ψ → c ] = e D · T · [ Γ 0 · Ψ → ] + [ Γ 0 · C ] · V → ‾ - - - ( 13 )
这种形式表示设置点状态特征向量
Figure C20061008485100132
和初始状态特征向量
Figure C20061008485100133
分别表示在设置点轴线系统中和在初始轴线系统中,这意味着可非常简单地表达向量和标量乘积。
矩阵(D、P、P-1、B)取决于所述机械,并且可以用分析方法或数值方法计算。控制间隔T基于其它标准选择。向量
Figure C20061008485100134
取决于在间隔T的开始处机械的初始状态,并且向量是设置点向量。对于每个控制步骤这两个向量
Figure C20061008485100136
与三个投影矩阵:(Γc、Γ0、ΓT)一起必须计算。ΓT依据应用可能在未知变量中确定为测量和/或设置点值的函数。
I.2-设置点电流的表达
这里设置点向量被定义,以能够同时实现例如以主转矩设置点Cc的形式表达的机械力设置点Fmc、和在间隔T的末端处要实现的设置点磁通
Figure C20061008485100137
这里只有设置点磁通
Figure C20061008485100138
的模量||Φc||是必需的。
在由电流向量
Figure C20061008485100139
和磁场向量
Figure C200610084851001310
定义的平面中,具有至少四维的状态向量选择成定义机械的状态:
X → = I → Φ → - - - ( 14 )
应用于设置点的这个状态向量基于磁性状态和机械力被定义。
为了简单化在各向同性异步机械的情况下的计算形式,电流和磁通向量可以由复数代表,每个在轴线的直接正交系统中具有两个坐标,其中实数横坐标轴在磁场向量的方向上,并因此在磁通向量的方向上。
系统(13)保持有四阶,但现在在复数平面中由二阶系统表示。
μ1和μ2定义状态方程的离散系统的特性的复数矩阵的两个复数特征值;复数控制矩阵和复数逆通过矩阵用这些维可写成如下:
B = b 1 b 2 - - - ( 15 )
P - 1 = p 11 p 22 p 21 p 22 使Pjk为μ1、μ2的函数(16)
控制向量是一阶的和复数的,并且由关系 V → ‾ = V → α + i · V → β · I 1 表达。便利的是,选择相对于控制单元固定的轴线系统的坐标(α,β)。
借助于以上定义,C=D-1·(eD·T-I)[P-1·B]成为:
C = e μ 1 · T - 1 μ 1 · ( p 11 · b 1 + p 12 · b 2 ) e μ 2 · T - 1 μ 2 ( p 21 · b 1 + p 22 · b 2 ) = a 1 a 2 - - - ( 17 )
让我们使用这些定义细化状态特征向量的方程系统:
Ψ → c = p 11 p 12 p 21 p 22 · I → c Φ → c = Ψ → 1 c Ψ → 2 c - - - ( 18 )
Ψ → 0 = p 11 p 12 p 21 p 22 · I → 0 Φ → 0 = Ψ → 10 Ψ → 20 - - - ( 19 )
Ψ → c = e μ 1 · T 0 0 e μ 2 · T · Ψ → 0 + a 1 a 2 · V → ‾ - - - ( 20 )
其中:
·
Figure C20061008485100149
分别是电流和设置点磁通向量;
·
Figure C200610084851001410
Figure C200610084851001411
分别是电流和磁通初始向量。
让:
e 1 = e μ 1 · T e 2 = e μ 2 · T
方程系统成为:
p 11 p 12 p 21 p 22 · I → c Φ → c = e 1 · p 11 e 1 · p 12 e 2 · p 21 e 2 · p 22 · I → 0 Φ → 0 + a 1 a 2 · V → ‾ - - - ( 21 )
它然后划分成两个复数方程:
p 11 · I → c + p 12 · Φ → c = e 1 · ( p 11 · I → 0 + p 12 · Φ → 0 ) + a 1 · V → ‾ - - - ( 22 )
p 21 · I → c + p 22 · Φ → c = e 2 · ( p 21 · I → 0 + p 22 · Φ → 0 ) + a 2 · V → ‾ - - - ( 23 )
这可等效地写成:
Ψ → 1 c = e 1 · Ψ → 10 + a 1 · V → ‾ - - - ( 24 )
Ψ → 2 c = e 2 · Ψ → 20 + a 2 · V → ‾ - - - ( 25 )
然后在方程(22)和(23)之间消去复数控制向量:
( a 2 · p 11 - a 1 · p 21 ) · I → c + ( a 2 · p 12 - a 1 · p 22 ) · Φ → c =
( a 2 · e 1 · p 11 - a 1 · e 2 · p 21 ) · I → 0 + ( a 2 · e 1 · p 12 - a 1 · e 2 · p 22 ) · Φ → 0 - - - ( 26 )
这种关系使得显然的是:设置点电流向量被链接到设置点磁通向量上,并且都取决于这个方程的第二项的已知初始状态。
取决于所述应用的系数(a2·p11-a1·p21)是非零的,因为旋转电气机械的磁通取决于机械的激励电流。因此有可能写出如下关系:
I → c + a 2 · p 12 - a 1 · p 22 a 2 · p 11 - a 1 · p 21 · Φ → c = ( a 2 · e 1 · p 11 - a 1 · e 2 · p 21 ) · I → 0 + ( a 2 · e 1 · p 12 - a 1 · e 2 · p 22 ) · Φ → 0 a 2 · p 11 - a 1 · p 21
(27)
注意,如果系数(a2·p11-a1·p21)是零,则跨过初始状态的设置点电流和设置点磁通的相互依赖性完全消失,这在实际中意味着不可能由电流向量控制磁通向量。磁通因此或者是恒定的或者以某种其它方式控制,并且这里描述的直进式控制方法不适用。
让:
E → 0 = ( a 2 · e 1 · p 11 - a 1 · e 2 · p 21 ) · I → 0 + ( a 2 · e 1 · p 12 - a 1 · e 2 · p 22 ) · Φ → 0 a 2 · p 11 - a 1 · p 21 = E 0 · e i · ϵ 0 - - - ( 28 )
得到关系(27)的更简明形式:
I → c + a 2 · p 12 - a 1 · p 22 a 2 · p 11 - a 1 · p 21 · Φ → c = E 0 · e i · ϵ 0 - - - ( 29 )
其中电流和磁通向量在设置点轴线系统中由复数表达式表示。
I → c = I cd + i · I cq - - - ( 30 )
Φ → c = Φ cd = | | Φ → c | | - - - ( 31 )
其中Icd是磁化电流,并且Icq是与电机的磁通相结合产生机械力的电流。它然后有可能写成:
I cd + i · I cq + a 2 · p 12 - a 1 · p 22 a 2 · p 11 - a 1 · p 21 · Φ c = E 0 · e i · ϵ 0 - - - ( 32 )
考虑两个复数向量:
Figure C20061008485100157
Figure C20061008485100158
Figure C20061008485100159
Figure C200610084851001510
共轭。已知的是:
其中
·代表向量的标量乘积的实数标量;
·代表向量乘积的模量,它也是实数。我们因此可以写出:
Re { M &RightArrow; * &CenterDot; N &RightArrow; } = < M &RightArrow; , N &RightArrow; > - - - ( 34 )
Figure C200610084851001515
我们现在计算:电流和磁通向量的标量乘积,它等于磁能Q,忽略常数kq;和相同向量的向量乘积,它等于机械力Fm,忽略常数km
I &RightArrow; c * &CenterDot; &Phi; &RightArrow; c = E 0 &CenterDot; e - i &CenterDot; &epsiv; 0 &CenterDot; &Phi; &RightArrow; c - a 2 * &CenterDot; p 12 * - a 1 * &CenterDot; p 22 * a 2 * &CenterDot; p 11 * - a 1 * &CenterDot; p 21 * &CenterDot; | | &Phi; &RightArrow; c | | 2 = Q k q + i &CenterDot; F m k m - - - ( 36 )
这种复数关系划分成两个实数方程:
Re { I &RightArrow; c * &CenterDot; &Phi; &RightArrow; c } = < I &RightArrow; c , &Phi; &RightArrow; c > = Q k q = I cd &CenterDot; | | &Phi; &RightArrow; c | | = E 0 &CenterDot; cos ( &epsiv; 0 ) &CenterDot; | | &Phi; &RightArrow; c | | - Re { a 2 &CenterDot; p 12 - a 1 &CenterDot; p 22 a 2 &CenterDot; p 11 - a 1 &CenterDot; p 21 } &CenterDot; | | &Phi; &RightArrow; c | | 2
(37)
并且,因为选择的轴线系统的意义:
Figure C20061008485100163
(38)
其中Re{...}和Im{...}分别指示复数的实部和虚部。
I cd + Re { a 2 &CenterDot; p 12 - a 1 &CenterDot; p 22 a 2 &CenterDot; p 11 - a 1 &CenterDot; p 21 } &CenterDot; | | &Phi; &RightArrow; c | | = E 0 &CenterDot; cos ( &epsiv; 0 ) - - - ( 39 )
I cq + I m { a 2 &CenterDot; p 12 - a 1 &CenterDot; p 22 a 2 &CenterDot; p 11 - a 1 &CenterDot; p 21 } &CenterDot; | | &Phi; &RightArrow; c | | = E 0 &CenterDot; sin ( &epsiv; 0 ) - - - ( 40 )
然后有可能以两种不同的形式通过消去ε0从这两个方程提取Icd
第一种形式:
ε0由第二方程计算,并且它由其在第一方程中的表达式替换:
&epsiv; 0 = arcsin ( I cq E 0 + Im { a 2 &CenterDot; p 12 - a 1 &CenterDot; p 22 a 2 &CenterDot; p 22 - a 1 &CenterDot; p 21 } &CenterDot; | | &Phi; &RightArrow; c | | E 0 ) - - - ( 41 )
然后:
I cd = - Re { a 2 &CenterDot; p 12 - a 1 &CenterDot; p 22 a 2 &CenterDot; p 11 - a 1 &CenterDot; p 21 } &CenterDot; | | &Phi; &RightArrow; c | | + E 0 &CenterDot; cos [ arcsin ( I cq E 0 + Im { a 2 &CenterDot; p 12 - a 1 &CenterDot; p 22 a 2 &CenterDot; p 11 - a 1 &CenterDot; p 21 } | | &Phi; &RightArrow; c | | E 0 ) ]
(42)
应用的上下文消去Arcsin()的两个主要确定的一个,具有负余弦的那个保留。
第二种形式:
两个方程的两项被平方和求和,以由其数值替换cos20)+sin20)=1:
[ I cd + Re { a 2 &CenterDot; p 12 - a 1 &CenterDot; p 22 a 2 &CenterDot; p 11 - a 1 &CenterDot; p 21 } &CenterDot; | | &Phi; &RightArrow; c | | ] 2 + [ I cq + Im { a 2 &CenterDot; p 12 - a 1 &CenterDot; p 22 a 2 &CenterDot; p 11 - a 1 &CenterDot; p 21 } &CenterDot; | | &Phi; &RightArrow; c | | ] 2 = E 0 2 - - - ( 43 )
这种新关系被扩展以计算Icd作为二次方程的解:
I cd = - Re { a 2 &CenterDot; p 12 - a 1 &CenterDot; p 22 a 2 &CenterDot; p 11 - a 1 &CenterDot; p 21 } &CenterDot; | | &Phi; &RightArrow; c | | &PlusMinus; E 0 2 - [ I cq + Im { a 2 &CenterDot; p 12 - a 1 &CenterDot; p 22 a 2 &CenterDot; p 11 - a 1 &CenterDot; p 21 } &CenterDot; | | &Phi; &RightArrow; c | | ] 2
(44)
平方引入了必须借助于所述应用的上下文消去的增根,只保留在根式前面具有负号的那个根。
进一步注意:
I cq &CenterDot; | | &Phi; &RightArrow; c | | = F m k m - - - ( 45 )
这个关系因此用来把机械力设置点Fmc和磁通模量设置点
Figure C20061008485100173
变换成用于把电流向量投影到与磁通向量轴正交的轴上的设置点:
I cq = F mc k m &CenterDot; | | &Phi; &RightArrow; c | | - - - ( 46 )
对于Icd的一个或另一个那么取决于:
·设置点Fmc
Figure C20061008485100175
·应用的参数:取决于特征值的a1、a2、p11、p12、p21、p21
·测量的初始条件:通过
Figure C20061008485100176
的中段的
Figure C20061008485100177
Figure C20061008485100178
把磁化电流链接到主要设置点、参数及初始状态上的关系的所述两种形式因此能够实现设置点电流的第二坐标的计算、设置点电流向量到磁通向量上的投影。
I cd = - Re { a 2 &CenterDot; P 12 - a 1 &CenterDot; P 22 a 2 &CenterDot; P 11 - a 1 P 21 } &CenterDot; | | &Phi; &RightArrow; c | |
+ E 0 &CenterDot; cos [ arcsin ( F mc k m &CenterDot; | | &Phi; &RightArrow; c | | E 0 &CenterDot; + Im { a 2 &CenterDot; P 12 - a 1 &CenterDot; P 22 a 2 &CenterDot; P 11 - a 1 &CenterDot; P 21 } &CenterDot; | | &Phi; &RightArrow; c | | E 0 ) ] - - - ( 47 )
I cd = Re { a 2 &CenterDot; P 12 - a 1 &CenterDot; P 22 a 2 &CenterDot; P 11 - a 1 &CenterDot; P 21 } &CenterDot; | | &Phi; &RightArrow; c | | &PlusMinus; E 0 2 - [ F mc k m &CenterDot; | | &Phi; &RightArrow; c | | + Im { a 2 &CenterDot; P 12 - a 1 &CenterDot; P 22 a 2 &CenterDot; P 11 - a 1 &CenterDot; P 21 } &CenterDot; | | &Phi; &RightArrow; c | | ] 2
(48)
总之,我们使用以前关系从(Fmc)得到(Icd、Icq)和(Φcd,0),使我们能够建造投影到设置点轴线系统中的设置状态特征向量:
&Gamma; c &CenterDot; &Psi; &RightArrow; c = P 11 &CenterDot; ( I cd + i &CenterDot; I cq ) + P 12 &CenterDot; &Phi; cd P 21 &CenterDot; ( I cd + i &CenterDot; I cq ) + P 22 &CenterDot; &Phi; cd = &Psi; &RightArrow; 1 c dq &Psi; &RightArrow; 2 c dq - - - ( 49 )
其来自关系:
&Gamma; T &CenterDot; [ &Gamma; c &CenterDot; &Psi; &RightArrow; c ] = e D &CenterDot; T &CenterDot; [ &Gamma; 0 &CenterDot; &Psi; &RightArrow; 0 ] + [ &Gamma; 0 &CenterDot; C ] &CenterDot; V &RightArrow; &OverBar; - - - ( 50 )
以这种方式建立的设置点电流向量使我们在间隔T的末端处同时得到机械力设置点Fmc和磁通设置点
Figure C20061008485100183
我们类似地可写出:
&Gamma; 0 &CenterDot; &Psi; &RightArrow; 0 = P 11 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + P 12 &CenterDot; &Phi; 0 d 0 P 21 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + P 22 &CenterDot; &Phi; 0 d 0 = &Psi; &RightArrow; 10 d 0 q 0 &Psi; &RightArrow; 20 d 0 q 0 - - - ( 51 )
初始状态由测量得知。
I.3-建立用于异步电动机的矩阵Γ T 的例子
因为控制向量这里在相对于定子固定的轴线系统中施加到三相旋转电动机上,所以选择由定子电流向量和转子磁通向量
Figure C20061008485100186
组成的状态向量,如在轴线正交系统中从定子看到的那样,该轴线正交系统通常叫做(α,β),在减小电机的磁极对数目之后相对于电机的定子固定。在这种基准正交轴线系统(α,β)中,Concordia变换用于从三相系统到两相系统的这种变化。α轴线可以在例如相位R的电极性方向上选择。
定子电流向量可借助于布置在电机外部的、并且与电机的相或一部分相串联的电流传感器而测量。转子磁通向量可从电机电压和电流得知,或者可被观察。这两个向量用来估计或计算初始状态。如有必要,它们也联合用来通过它们的向量乘积并由此通过它们的坐标以及磁能设置点来移动机械转矩设置点。
选择的状态向量因此在这种情况下是4维的。有四个电气自由度。它可以用选择成简化符号的两个向量或者用在轴线复数系统(α,β)中代表它们的复数互换地指示,β是纯虚轴,以通过把矩阵的秩从四减小到二(在复数平面中)来简化特征值的计算。
X &RightArrow; = I &RightArrow; s &Phi; &RightArrow; r = I s&alpha; + i &CenterDot; I s&beta; &Phi; r&alpha; + i &CenterDot; &Phi; r&beta; - - - ( 52 )
使用复数的第二个优点在于,它通过简单的指数复数标量而不是使用2维平方复数矩阵变换坐标轴的转动。
矩阵ΓT然后能以如下复数形式写出:
&Gamma; T = e i &CenterDot; &Delta; &theta; p - - - ( 53 )
其中Δθp代表在间隔T期间磁通的预测角转动。
借助于状态变化变量的这种选择,可规定初始状态特征向量:
&Gamma; 0 &CenterDot; &Psi; &RightArrow; 0 = &Psi; &RightArrow; 10 d 0 q 0 &Psi; &RightArrow; 20 d 0 q 0 = p 11 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + p 12 &CenterDot; &Phi; 0 q 0 p 21 &CenterDot; ( I 0 d 0 + I 0 q 0 ) + p 22 &CenterDot; &Phi; 0 d 0 - - - ( 54 )
变量
Figure C20061008485100193
代表投影到由相对于固定基准轴线系统(α,β)的变换 &Gamma; 0 = e - i &CenterDot; &theta; 0 所限定的初始轴线系统中的状态特征向量的初始状态。初始值
Figure C20061008485100195
&Gamma; 0 = e - i &CenterDot; &theta; 0 由测量的电流和观察的磁通得知。
投影到基准的初始和设置点框架中的状态特征向量的一般方程成为:
e i &CenterDot; &Delta;&theta; p &CenterDot; P 11 &CenterDot; ( I cd + i &CenterDot; I cq ) + P 12 &CenterDot; &Phi; cd P 21 &CenterDot; ( I cd + i &CenterDot; I cq ) + P 22 &CenterDot; &Phi; cd = e &mu; 1 &CenterDot; T 0 0 e &mu; 2 &CenterDot; T &CenterDot; P 11 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + P 12 &CenterDot; &Phi; 0 d 0 P 21 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + P 22 &CenterDot; &Phi; 0 d 0 + e - i &CenterDot; &theta; 0 &CenterDot; a 1 a 2 &CenterDot; V &RightArrow; &OverBar;
(55)
这个系统能以两个复数方程的形式再次细化:
e i &CenterDot; &Delta; &theta; p &CenterDot; [ p 11 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 12 &CenterDot; &Phi; cd ] = e &mu; 1 &CenterDot; T &CenterDot; [ p 11 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + p 12 &CenterDot; &Phi; 0 d 0 ] + e - i &CenterDot; &theta; 0 &CenterDot; a 1 &CenterDot; V &RightArrow; &OverBar;
(56)
e i &CenterDot; &Delta; &theta; p &CenterDot; [ p 21 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 22 &CenterDot; &Phi; cd ] = e &mu; 2 &CenterDot; T &CenterDot; [ p 21 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + p 22 &CenterDot; &Phi; 0 d 0 ] + e - i &CenterDot; &theta; 0 &CenterDot; a 2 &CenterDot; V &RightArrow; &OverBar;
(57)
在这两个复数方程之间再次消去控制向量:
&Gamma; T = e i &CenterDot; &Delta; &theta; p = e &mu; 1 &CenterDot; T &CenterDot; a 2 &CenterDot; [ p 11 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + p 12 &CenterDot; &Phi; 0 d 0 ] - e &mu; 2 &CenterDot; T &CenterDot; a 1 &CenterDot; [ p 21 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + p 22 &CenterDot; &Phi; 0 d 0 ] a 2 &CenterDot; [ p 11 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 12 &CenterDot; &Phi; cd ] - a 1 &CenterDot; [ p 21 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 22 &CenterDot; &Phi; cd ]
(58)
或者:
&Gamma; T = e &mu; 1 &CenterDot; T &CenterDot; a 2 &CenterDot; &Psi; &RightArrow; 10 d 0 q 0 - e &mu; 2 &CenterDot; T &CenterDot; a 1 &CenterDot; &Psi; &RightArrow; 20 d 0 q 0 a 2 &CenterDot; &Psi; &RightArrow; 1 c dq - a 1 &CenterDot; &Psi; &RightArrow; 2 c dq - - - ( 59 )
现在可计算在初始轴线系统与设置点轴线系统之间的变换。
总之,
Figure C200610084851001912
现在完全由主设置点(Fmc
Figure C200610084851001913
)定义。系统
&Gamma; T &CenterDot; [ &Gamma; c &CenterDot; &Psi; &RightArrow; c ] = e D &CenterDot; T &CenterDot; [ &Gamma; 0 &CenterDot; &Psi; &RightArrow; 0 ] + [ &Gamma; 0 &CenterDot; C ] &CenterDot; V &RightArrow; &OverBar; - - - ( 60 )
是4阶(实数)或两阶(复数),并且已经只使用了两次:在旋转电机的情况下一次一般计算Icd和一次计算Δθp
I.4-控制向量的计算
状态方程系统的唯一未知量现在是控制向量的两个实数坐标,它们是如下两个复数方程之一的解:
e i &CenterDot; &Delta; &theta; p &CenterDot; [ p 11 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 12 &CenterDot; &Phi; cd ] = e &mu; 1 &CenterDot; T &CenterDot; [ p 11 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + p 12 &CenterDot; &Phi; 0 d 0 ] + e - i &CenterDot; &theta; 0 &CenterDot; a 1 &CenterDot; V &RightArrow; &OverBar; - - - ( 61 )
e i &CenterDot; &Delta; &theta; p &CenterDot; [ p 21 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 22 &CenterDot; &Phi; cd ] = e &mu; 2 &CenterDot; T &CenterDot; [ p 21 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + p 22 &CenterDot; &Phi; 0 d 0 ] + e - i &CenterDot; &theta; 0 &CenterDot; a 2 &CenterDot; V &RightArrow; &OverBar; - - - ( 62 )
它们可通过现在消去ΓT简单地形成:
V &RightArrow; &OverBar; = e i &CenterDot; &theta; 0 &CenterDot; e &mu; 1 &CenterDot; T &CenterDot; [ p 11 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + p 12 &CenterDot; &Phi; 0 d 0 ] &CenterDot; [ p 21 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 22 &CenterDot; &Phi; cd ] - e &mu; 2 &CenterDot; T &CenterDot; [ p 21 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + p 22 &CenterDot; &Phi; 0 d 0 ] &CenterDot; [ p 11 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 12 &CenterDot; &Phi; cd ] a 2 [ p 11 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 12 &CenterDot; &Phi; cd ] - a 1 &CenterDot; [ p 21 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 22 &CenterDot; &Phi; cd ]
(63)
V &RightArrow; &OverBar; = e i &CenterDot; &theta; 0 &CenterDot; e &mu; 1 &CenterDot; T &CenterDot; &Psi; &RightArrow; 10 d 0 q 0 &CenterDot; &Psi; &RightArrow; 2 c dq - e &mu; 2 &CenterDot; T &CenterDot; &Psi; &RightArrow; 20 d 0 q 0 &CenterDot; &Psi; &RightArrow; 1 c dq a 2 &CenterDot; &Psi; &RightArrow; 1 c dq - a 1 &CenterDot; &Psi; &RightArrow; 2 c dq - - - ( 64 )
这是复数关系,该复数关系供给最后两个未知量:控制向量的坐标。
现在可得出结论:在各向同性旋转电机的情况下已经转化4维状态方程系统,并且我们已经从主设置点(Fmc
Figure C20061008485100205
)求出用于四个未知量:Icd、Δθp
Figure C20061008485100206
的准确分析表达式。
部分II-对于异步旋转电机的应用
异步旋转电机是各向同性旋转电机。
在异步电机的情况下,在Concordia三相到两相变换之后,使用图1的等效电路,它把转子和定子的漏电感与电机的定子归到一起。这个等效电路具有简单表示的优点,把定子电流分离成其两个分量:产生转子磁通的电流和产生转矩的电流。
II.1-异步电机状态表示
使用对应电气微分方程,容易表明,连续状态方程的线性系统可以以如下形式写出:
I &CenterDot; s&alpha; + i &CenterDot; I &CenterDot; s&beta; &Phi; &CenterDot; r&alpha; + i &CenterDot; &Phi; &CenterDot; r&beta; = - &gamma; &beta; &CenterDot; ( &alpha; - i &CenterDot; &omega; ) &alpha; &CenterDot; L m - ( &alpha; - i &CenterDot; &omega; ) &CenterDot; I s&alpha; + i &CenterDot; I s&beta; &Phi; r&alpha; + i &CenterDot; &Phi; r&beta; + 1 &lambda; 0 &CenterDot; ( V s&alpha; + i &CenterDot; V s&beta; ) - - - ( 65 )
并因此:
A = - &gamma; &beta; &CenterDot; ( &alpha; - i &CenterDot; &omega; ) &alpha; &CenterDot; L m - ( &alpha; - i &CenterDot; &omega; ) B = 1 &lambda; 0 - - - ( 66 )
有:
·定子电阻:      Rs
·转子电阻:      Rr
·定子漏电感:    ls
·转子漏电感:    lr
·互感:          Lm
·定子电感:      Ls=Lm+ls
·转子电感:      Lr=Lm+lr
·定子时间常数: &tau; s = L s R s
·转子时间常数: &tau; r = L r R r
·极对数目:      NP
·扩散系数: &sigma; = 1 - L 2 m L &Gamma; &CenterDot; L s
机械和电气角频率的定义和符号:
·转子的机械角频率:Ω
·极性机械角频率:ω=Np·Ω
·定子电气角频率:ωs
·转子电气角频率:ωr
·转差: g = ( &omega; s - &omega; ) &omega; s = &omega; r &omega; s
及:
R sr = R s + R r &CenterDot; L m 2 L r 2 - - - ( 67 )
&alpha; = 1 &tau; r - - - ( 68 )
&beta; = L m &sigma; &CenterDot; L s &CenterDot; L r = 1 - &sigma; &sigma; &CenterDot; 1 L m - - - ( 69 )
λ=σ·Ls    (70)
&gamma; = 1 &sigma; &CenterDot; [ 1 &tau; s + ( 1 - &sigma; ) &tau; r ] = R sr &lambda; - - - ( 71 )
II.2-在特征向量基中状态方程系统的表达
给出矩阵A的特征值μ1的特征方程是:
μ2+(α+γ-i·ω)·μ+(γ-α·β·Lm)·(α-i·ω)=0(72)
让Δ指示方程的判别式:
Δ=(α-γ-i·ω)2+4·α·β·Lm·(α-i·ω)    (73)
两个特征值因此由下式表示:
&mu; 1 = - 1 2 &CenterDot; ( &alpha; + &gamma; - i &CenterDot; &omega; - &Delta; ) - - - ( 74 )
&mu; 2 = - 1 2 &CenterDot; ( &alpha; + &gamma; - i &CenterDot; &omega; + &Delta; ) - - - ( 75 )
这里注意,特征值取决于极性机械角频率。它们因此必须由电机的模型、测量速度及极对数目计算。
ω=Np·Ω    (76)
对角化矩阵A是:
D = &mu; 1 0 0 &mu; 2 - - - ( 77 )
与特征值相对应的特征向量 ( &Pi; &RightArrow; 1 = p 11 p 21 , &Pi; &RightArrow; 2 = p 12 p 22 ) 是方程的一个解:
( A - &mu; 1 &CenterDot; I ) &CenterDot; &Pi; &RightArrow; i = 0 - - - ( 78 )
对于第一特征向量,定义到在乘法常数内:
- &gamma; &beta; &OverBar; &CenterDot; ( &alpha; - i &CenterDot; &omega; ) &alpha; &CenterDot; L m - ( &alpha; - i &CenterDot; &omega; ) &CenterDot; p 11 p 21 = &mu; 1 &CenterDot; p 11 p 21 - - - ( 79 )
矩阵方程的第二行产生如下关系:
α·Lm·p11=(μ1+α-i·ω)·p21    (80)
它使我们能够选择:
p 11 = 2 &CenterDot; ( &mu; 1 + &alpha; - i &CenterDot; &omega; ) 2 &CenterDot; &alpha; &CenterDot; L m = &alpha; - &gamma; - i &CenterDot; &omega; + &Delta; 2 &CenterDot; &alpha; &CenterDot; L m - - - ( 81 )
p21=1    (82)
以类似方式,对于第二特征向量,定义到在乘法常数内:
p 12 = 2 &CenterDot; ( &mu; 2 + &alpha; - i &CenterDot; &omega; ) 2 &CenterDot; &alpha; &CenterDot; L m = &alpha; - &gamma; - i &CenterDot; &omega; - &Delta; 2 &CenterDot; &alpha; &CenterDot; L m - - - ( 83 )
p22=1
通过矩阵 P = &Pi; &RightArrow; 1 &Pi; &RightArrow; 2 = p 11 p 12 p 21 p 22 从而:因此可以写A·P=P·D
P = &alpha; - &gamma; - i &CenterDot; &omega; + &Delta; 2 &CenterDot; &alpha; &CenterDot; L m &alpha; - &gamma; - i &CenterDot; &omega; - &Delta; 2 &CenterDot; &alpha; &CenterDot; L m 1 1 - - - ( 85 )
或者,使用特征值:
P = 1 &alpha; &CenterDot; L m &CenterDot; - ( &mu; 2 + &gamma; ) - ( &mu; 1 + &gamma; ) &alpha; &CenterDot; L m &alpha; &CenterDot; L m - - - ( 86 )
并且因此可以写出其逆矩阵:
P - 1 = 1 2 &Delta; &CenterDot; 2 &CenterDot; &alpha; &CenterDot; L m - ( &alpha; - &gamma; - i &CenterDot; &omega; - &Delta; ) - 2 &CenterDot; &alpha; &CenterDot; L m ( &alpha; - &gamma; - i &CenterDot; &omega; + &Delta; ) - - - ( 87 )
或者,使用特征值:
P - 1 = 1 &mu; 1 - &mu; 2 &CenterDot; &alpha; &CenterDot; L m ( &mu; 1 + &gamma; ) - &alpha; &CenterDot; L m - ( &mu; 2 + &gamma; ) - - - ( 88 )
我们因此可写出:
A=P·D·P-1    (89)
对于使用特征值的符号,A成为:
A = - &gamma; - ( &mu; 1 + &gamma; ) &CenterDot; ( &mu; 2 + &gamma; ) &alpha; &CenterDot; L m &alpha; &CenterDot; L m ( &mu; 1 + &mu; 2 + &gamma; ) - - - ( 90 )
使用由结合中间变量生成的如下关系:
1+γ)·(μ2+γ)=α·β·Lm·(μ12+γ)    (91)
A = - &gamma; - ( &mu; 1 + &gamma; ) &CenterDot; ( &mu; 2 + &gamma; ) &alpha; &CenterDot; L m &alpha; &CenterDot; L m ( &mu; 1 + &gamma; ) &CenterDot; ( &mu; 2 + &gamma; ) &alpha; &CenterDot; L m &CenterDot; &beta; - - - ( 92 )
使用如下减小变量:
&xi; 0 = &mu; 1 - &mu; 2 &alpha; &CenterDot; L m , &xi; 1 = &mu; 1 + &gamma; &alpha; &CenterDot; L m , &xi; 2 = &mu; 2 + &gamma; &alpha; &CenterDot; L m - - - ( 93 )
概括成:
A = &alpha; &CenterDot; L m &CenterDot; - &xi; - &xi; 1 &CenterDot; &xi; 2 1 &xi; 1 &CenterDot; &xi; 2 &beta; , D = &mu; 1 0 0 &mu; 2 - - - ( 94 )
P = - &xi; 2 - &xi; 1 1 1 , P - 1 = 1 &xi; 0 &CenterDot; 1 &xi; 1 - 1 - &xi; 2 - - - ( 95 )
用F=P·eD·T·P-1替换F=eA·T,我们可现在写出:
[ P - 1 &CenterDot; X &RightArrow; c ] = e D &CenterDot; T &CenterDot; [ P - 1 &CenterDot; X &RightArrow; 0 ] + D - 1 &CenterDot; ( e D &CenterDot; T - 1 ) &CenterDot; [ P - 1 &CenterDot; B ] &CenterDot; V &RightArrow; &OverBar; - - - ( 96 )
其中:
e D &CenterDot; T = e &mu; 1 &CenterDot; T 0 0 e &mu; 2 &CenterDot; T D - 1 = 1 &mu; 1 0 0 1 &mu; 2 P - 1 &CenterDot; B = 1 &lambda; &CenterDot; &xi; 0 &CenterDot; 1 - 1 - - - ( 97 )
得到特别简单的形式:
1 &xi; 1 1 &xi; 2 &CenterDot; X &RightArrow; c = e &mu; 1 &CenterDot; T 0 0 e &mu; 2 &CenterDot; T &CenterDot; 1 &xi; 1 1 &xi; 2 &CenterDot; X &RightArrow; 0 + V &RightArrow; &OverBar; &lambda; &CenterDot; e &mu; 1 &CenterDot; T - 1 &mu; 1 e &mu; 2 &CenterDot; T - 1 &mu; 2 - - - ( 98 )
我们识别方法的一般形式:
&Psi; &RightArrow; c = e D &CenterDot; T &CenterDot; &Psi; &RightArrow; 0 + C &CenterDot; V &RightArrow; &OverBar; - - - ( 99 )
其中:
&Psi; &RightArrow; c = 1 &xi; 1 1 &xi; 2 &CenterDot; X &RightArrow; c = &Psi; &RightArrow; 1 c &Psi; &RightArrow; 2 c &Psi; &RightArrow; 0 = 1 &xi; 1 1 &xi; 2 &CenterDot; X &RightArrow; 0 = &Psi; &RightArrow; 10 &Psi; &RightArrow; 20 - - - ( 100 )
e D &CenterDot; T = e &mu; 1 &CenterDot; T 0 0 e &mu; 2 &CenterDot; T , C = 1 &lambda; &CenterDot; e &mu; 1 &CenterDot; T - 1 &mu; 1 e &mu; 2 &CenterDot; T - 1 &mu; 2 - - - ( 101 )
状态“特征向量”的坐标可以通过用它们的向量坐标替换状态向量而使显得非常简单:
I &RightArrow; sc + &xi; 1 &CenterDot; &Phi; &RightArrow; rc I &RightArrow; sc + &xi; 2 &CenterDot; &Phi; &RightArrow; rc = e &mu; 1 &CenterDot; T &CenterDot; ( I &RightArrow; s 0 + &xi; 1 &CenterDot; &Phi; &RightArrow; r 0 ) e &mu; 2 &CenterDot; T &CenterDot; ( I &RightArrow; s 0 + &xi; 2 &CenterDot; &Phi; &RightArrow; r 0 ) + V &RightArrow; &OverBar; &lambda; &CenterDot; e &mu; 1 &CenterDot; T - 1 &mu; 1 e &mu; 2 &CenterDot; T - 1 &mu; 2 - - - ( 102 )
状态的坐标、初始和设置点“特征向量”现在是:
&Psi; &RightArrow; 10 = I &RightArrow; s 0 + &xi; 1 &CenterDot; &Phi; &RightArrow; r 0 , &Psi; &RightArrow; 20 = I &RightArrow; s 0 + &xi; 2 &CenterDot; &Phi; &RightArrow; r 0 - - - ( 103 )
&Psi; &RightArrow; 1 c = I &RightArrow; sc + &xi; 1 &CenterDot; &Phi; &RightArrow; rc , &Psi; &RightArrow; 2 c = I &RightArrow; sc + &xi; 2 &CenterDot; &Phi; &RightArrow; rc - - - ( 104 )
它通过设置产生如下关系:
a 1 = e &mu; 1 &CenterDot; T - 1 &lambda; &CenterDot; &mu; 1 , a 2 = e &mu; 2 &CenterDot; T - 1 &lambda; &CenterDot; &mu; 2 - - - ( 105 )
&Psi; &RightArrow; 1 c = e &mu; 1 &CenterDot; T &CenterDot; &Psi; &RightArrow; 10 + a 1 &CenterDot; V &RightArrow; &OverBar; - - - ( 106 )
&Psi; &RightArrow; 2 c = e &mu; 2 &CenterDot; T &CenterDot; &Psi; &RightArrow; 20 + a 2 &CenterDot; V &RightArrow; &OverBar; - - - ( 107 )
当电机的参数已知时,对于给定机械速度ω和对于给定预测水平T计算在方程中涉及的各种减小参数是简单问题。初始状态特征向量在轴线系统(α,β)中也由观察器所测量或估计的测量电流
Figure C200610084851002413
和磁通
Figure C20061008485100251
得知。因而,为了能够计算未知量:在轴线系统(α,β)中的控制向量
Figure C20061008485100252
只要保持把设置点状态特征向量固定在预测水平T处。
II.3-建立设置点电流:
为了固定对于异步电机特定的转子磁通和电磁转矩设置点,现在必须由Laplace定律,固定产生磁能的电流和磁通向量的标量乘积、和它们产生转矩的向量乘积。
现在注意,这里提到的电流和磁通向量在同一时刻-在间隔T的末端考虑。此外我们准备计算待施加的控制向量,从而通过模型考虑到电气时间常数在间隔T的末端处同时得到这些设置点。计算因此是必须没有两个轴线的随动或去耦的准确计算。
对于这种计算,初始轴线系统(d0,q0)和设置点轴线系统如在图2A和2B中指示的那样分别定义。为了更精确,轴线系统(d0,q0)是其横坐标轴与初始磁通向量
Figure C20061008485100254
对准的正交轴线系统,并且轴线系统
Figure C20061008485100255
是其横坐标轴与设置点磁通向量
Figure C20061008485100256
对准的正交轴线系统。
因为两个设置点向量定位在水平T处,即在间隔T的末端处,所以在间隔T的末端处预测轴线系统
Figure C20061008485100257
由此以新符号预测。其位置是瞬时位置。
在初始时刻的转子磁通向量的幅角ρ0是已知的,由于在此时刻的磁通是已知的。在水平T处的转子磁通向量的幅角
Figure C20061008485100258
是未知的。然而,磁通向量转动期间的间隔T是已知的,因为它按照其它考虑而选择的。让我们把在这个间隔期间的磁通向量的瞬时角速度ωs当作新的末知量。由此我们推出:
&rho; ~ = &omega; s &CenterDot; T + &rho; 0 - - - ( 108 )
我们把在状态特征向量之间的关系的两项乘以:
ΓT·Γc=Γ0    (109)
在这种情况下,从关系(108)选择:
e i &CenterDot; &omega; s &CenterDot; T &CenterDot; e - i &CenterDot; &rho; &RightArrow; = e - i &CenterDot; &rho; 0 - - - ( 110 )
它等效于反向转动,或等效于从轴线系统(α,β)到对于设置点的设置点轴线系统
Figure C200610084851002511
和到对于初始状态的初始轴线系统(d0,q0)的坐标轴变化。
e i &CenterDot; &omega; s &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 1 c &CenterDot; e - i &CenterDot; &rho; ~ ] = e &mu; 1 &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 10 &CenterDot; e - i &rho; ~ 0 ] + a 1 &CenterDot; e - i &CenterDot; &rho; ~ 0 &CenterDot; V &RightArrow; &OverBar; - - - ( 111 )
e i &CenterDot; &omega; s &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 2 c &CenterDot; e - i &CenterDot; &rho; ~ ] = e &mu; 2 &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 20 &CenterDot; e - i &rho; ~ 0 ] + a 2 &CenterDot; e - i &CenterDot; &rho; ~ 0 &CenterDot; V &RightArrow; &OverBar; - - - ( 112 )
其中:
&Psi; &RightArrow; jc &CenterDot; e - i &CenterDot; &rho; ~ = I c d ~ + i &CenterDot; I c q ~ + &xi; j &CenterDot; &Phi; c d ~ j &Element; { 1 ; 2 } - - - ( 113 )
&Psi; &RightArrow; j 0 &CenterDot; e - i &CenterDot; &rho; 0 = I 0 d 0 + i &CenterDot; I 0 q 0 + &xi; j &CenterDot; &Phi; 0 d 0 - - - ( 114 )
设置点向量现在作为要求转矩的函数和作为必需的磁通的模量、两个主设置点和原始状态向量的函数而容易计算,并且原始状态向量是已知的。
如下关系通过消去控制向量
Figure C20061008485100265
从关系(111)至(114)导出:
e i &CenterDot; &omega;s &CenterDot; T &CenterDot; [ ( I c d &OverBar; + i &CenterDot; I c q &OverBar; ) + &Phi; c d ~ &CenterDot; &xi; 1 &CenterDot; a 2 - &xi; 2 &CenterDot; a 1 a 2 - a 1 ] = e &mu; 1 &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 10 &CenterDot; e - i &CenterDot; &rho; 0 ] &CenterDot; a 2 - e &mu; 2 &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 20 &CenterDot; e - i &CenterDot; &rho; 0 ] &CenterDot; a 1 a 2 - a 1
(115)
第二项完全是已知的;我们按如下定义模量和幅角:
&eta; 0 &CenterDot; e i &CenterDot; &zeta; = e &mu; 1 &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 10 &CenterDot; e - i &CenterDot; &rho; 0 ] &CenterDot; a 2 - e &mu; 2 &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 20 &CenterDot; e - i &CenterDot; &rho; 0 ] a 1 a 2 - a 1 = E 0 &CenterDot; e i &CenterDot; &epsiv; 0 &CenterDot; e - i &CenterDot; &rho; 0 - - - ( 116 )
E 0 &CenterDot; e i &CenterDot; &epsiv; 0 = e &mu; 1 &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 10 ] &CenterDot; a 2 - e &mu; 2 &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 20 ] &CenterDot; a 1 a 2 - a 1 - - - ( 117 )
它使我们能够简化以上表达式:
( I c d &OverBar; + i &CenterDot; I c q &OverBar; ) &Phi; c d &OverBar; &CenterDot; &xi; 1 &CenterDot; a 2 - &xi; 2 &CenterDot; a 1 a 2 - a 1 = &eta; 0 &CenterDot; e i &CenterDot; ( &zeta; - &omega;s &CenterDot; T ) - - - ( 118 )
这种表达式通过取实部和虚部解出:
I c d ~ = - &Phi; c d ~ &CenterDot; Re { &xi; 1 &CenterDot; a 2 - &xi; 2 &CenterDot; a 1 a 2 - a 1 } + &eta; 0 &CenterDot; cos ( &zeta; - &omega;s &CenterDot; T ) - - - ( 119 )
I c q ~ = - &Phi; c d ~ &CenterDot; Im { &xi; 1 &CenterDot; a 2 - &xi; 2 &CenterDot; a 1 a 2 - a 1 } + &eta; 0 &CenterDot; sin ( &zeta; - &omega;s &CenterDot; T ) - - - ( 120 )
并且消去ωs
I c d ~ = - &Phi; c d &OverBar; &CenterDot; Re { &xi; 1 &CenterDot; a 2 - &xi; 2 &CenterDot; a 1 a 2 - a 1 } + &eta; 0 &CenterDot; cos ( arcsin [ I c q &OverBar; &eta; 0 + &Phi; c q ~ &eta; 0 &CenterDot; Im { &xi; 1 &CenterDot; a 2 - &xi; 2 &CenterDot; a 1 a 2 - a 1 } ] )
(121)
这是对于瞬态和稳定状态条件用来把设置点瞬态激励电流链接到设置点磁通和链接到设置点转矩上的所需的关系。
这种关系的另一种形式可以通过平方以及求和实数和虚数部分得到,这消去角度:
( I c d ~ + &Phi; rc &CenterDot; Re { &xi; 1 &CenterDot; a 2 - &xi; 2 &CenterDot; a 1 a 2 - a 1 } ) 2 + ( I c q ~ + &Phi; rc &CenterDot; Im { &xi; 1 &CenterDot; a 2 - &xi; 2 &CenterDot; a 1 a 2 - a 1 } ) 2 - &eta; 0 2 = 0
(122)
这个方程解出非常简单:
I c d ~ = - Re { &xi; 1 &CenterDot; a 2 - &xi; 2 &CenterDot; a 1 a 2 - a 1 } &CenterDot; &Phi; rc - &eta; 0 2 - ( I c q ~ + Im { &xi; 1 &CenterDot; a 2 - &xi; 2 &CenterDot; a 1 a 2 - a 1 } &CenterDot; &Phi; rc ) 2 - - - ( 123 )
我们在这种特定情形下识别已经以一般方式建立的解的两种形式。
作为设置点转矩和设置点磁通模量的函数的异步电机的电流设置点
Figure C20061008485100273
由如下关系给出:
I c q ~ = C c N p &CenterDot; L m L r &CenterDot; | | &Phi; &RightArrow; r | | | | &Phi; &RightArrow; r | | = &Phi; rc - - - ( 124 )
只剩下计算两个未知量
Figure C20061008485100276
和ωs
II.4建立矩阵Γ T
矩阵ΓT由如下关系给出:
&Gamma; T = e i &CenterDot; &omega; s &CenterDot; T = e &mu; 1 &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 10 &CenterDot; e - i &CenterDot; &rho; 0 ] &CenterDot; a 2 - e &mu; 2 &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 20 &CenterDot; e - i &CenterDot; &rho; 0 ] &CenterDot; a 1 [ &Psi; &RightArrow; 1 c &CenterDot; e - i &CenterDot; &rho; &OverBar; ] &CenterDot; a 2 - [ &Psi; &RightArrow; 2 c &CenterDot; e - i &CenterDot; &rho; &OverBar; ] &CenterDot; a 1 - - - ( 125 )
II.5计算控制向量:
控制向量由如下关系给出:
V &RightArrow; &OverBar; = e i &CenterDot; &rho; 0 &CenterDot; e &mu; 1 &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 10 &CenterDot; e - i &CenterDot; &rho; 0 ] &CenterDot; [ &Psi; &RightArrow; 2 c &CenterDot; e - i &CenterDot; &rho; ~ ] - e &mu; 2 &CenterDot; T &CenterDot; [ &Psi; &RightArrow; 20 &CenterDot; e - i &CenterDot; &rho; 0 ] &CenterDot; [ &Psi; &RightArrow; 1 c &CenterDot; e - i &CenterDot; &rho; ~ ] [ &Psi; &RightArrow; 1 c &CenterDot; e - i &CenterDot; &rho; &OverBar; ] &CenterDot; a 2 - [ &Psi; &RightArrow; 2 c &CenterDot; e - i &CenterDot; &rho; &OverBar; ] &CenterDot; a 1 - - - ( 126 )
我们识别与电机对应的控制向量的一般解。
现在可在其中初始向量被测量的轴线系统(α,β)中写出初始向量,留下在其中计算设置点向量的设置点系统中的设置点向量。
V &RightArrow; &OverBar; = e &mu; 1 &CenterDot; T &CenterDot; &Psi; &RightArrow; 10 &CenterDot; [ &Psi; &RightArrow; 2 c &CenterDot; e - i &CenterDot; &rho; &OverBar; ] - e &mu; 2 &CenterDot; T &CenterDot; &Psi; &RightArrow; 20 &CenterDot; [ &Psi; &RightArrow; 1 c &CenterDot; e - i &CenterDot; &rho; &OverBar; ] [ &Psi; &RightArrow; 1 c &CenterDot; e - i &CenterDot; &rho; ~ ] &CenterDot; a 2 - [ &Psi; &RightArrow; 2 c &CenterDot; e - i &CenterDot; &rho; ~ ] &CenterDot; a 1 - - - ( 127 )
部分III-具体实施例:
图3表示这里是带有定子和转子的三相异步旋转电动机的旋转电气机械12的直进式控制系统10。
控制系统10包括由脉冲宽度调制控制的电压转换器14,以产生供给电机12的每一相的电流和电压。更准确地说,在控制单元18的控制下,转换器14把来自直流电压源16的直流电压转换成用于电机12的三相供给电压。
控制单元18适于接收设置点Φc、转子磁通、转矩设置点Cc及间隔T的值,并且也适于接收由最多三个电流传感器20至22和最多三个电压传感器24至26实现的测量。传感器20至22的每一个适于测量在电机12的一相中的定子电流。传感器24至26的每一个适于测量在电机12的相应相中的定子电压。
单元18也适于接收由传感器28实现的电机12的转子的角速度的测量。然而,转子的机械速度的测量可以由其估计代替,例如由扩展KALMAN估计代替。
由这些测量、间隔T的值及设置点Φc和Cc,单元18能够计算待施加到用于电机12的转换器14上的控制向量,以在间隔T的末端处到达转矩设置点Cc和磁通设置点Φc
在电压转换器的情况下,控制向量是其坐标在转换器的固定轴线系统中定义的电压向量。
为了计算这个控制向量,单元18特别包括:模块30,用来建立能够在间隔T的末端处实现设置点Φc和Cc的设置点电流;和模块32,用来计算为在该间隔T的末端处得到与由模块30建立的设置点电流相等的电流而在间隔T期间连续施加的控制向量。
单元18也包括适于从电流和/或电压测量估计电机12的转子磁通的观察器33。
单元18典型地基于适于执行在信息存储介质34上存储的指令的常规可编程计算。为此,介质34存储用来当由单元18执行这些指令时用来执行图4的指令。
其次参照图4的方法描述系统10如何工作。
当设计用于电机12的直进式控制方法时,在步骤40中,由连接到转换器14上的电机12的微分电气方程,确定由转换器14控制的电机12的连续状态方程的线性系统,在该步骤40期间,常数σ、τr、α、β、λ及γ被计算,并且存储在例如介质34上。
当执行控制方法时,并且在间隔T期间计算控制向量的连续施加之前,在步骤42中选择设置点Φc、Cc及间隔T的值。
然后在步骤44期间,传感器28测量转子的机械角频率Ω。
然后在步骤46中,计算其值是机械角频率Ω的函数的电机参数。为了更准确,在步骤46期间,计算极性机械角频率ω、判别式Δ(由关系73定义)、特征值μ1和μ2(分别为关系74和75)、系数ξ1和ξ2(关系93)及系数a1和a2(关系105,也取决于T)。
在步骤50期间测量初始定子电流,并且在步骤52期间由观察器33估计初始转子磁通。
在步骤54期间,由在步骤52期间建立的初始磁通向量
Figure C20061008485100291
的估计计算幅角ρ0.
在步骤56中计算初始特征向量
Figure C20061008485100292
Figure C20061008485100293
为了更准确,在步骤56中借助于如下关系计算系数
Figure C20061008485100294
Figure C20061008485100295
&Psi; &RightArrow; 10 &CenterDot; e - i &CenterDot; &rho; 0 = I &RightArrow; s 0 &CenterDot; e - i &CenterDot; &rho; 0 + &xi; 1 &CenterDot; | | &Phi; &RightArrow; r 0 | | - - - ( 128 )
&Psi; &RightArrow; 20 &CenterDot; e - i &CenterDot; &rho; 0 = I &RightArrow; s 0 &CenterDot; e - i &CenterDot; &rho; 0 + &xi; 2 &CenterDot; | | &Phi; &RightArrow; r 0 | | - - - ( 129 )
在与步骤50至56并行执行的步骤60期间,借助于关系(123)和(124)建立电流设置点。
在步骤62中建立特征向量
Figure C20061008485100298
Figure C20061008485100299
为了更准确,在步骤62中借助于如下关系计算因数
Figure C200610084851002910
Figure C200610084851002911
&Psi; &RightArrow; 1 c &CenterDot; e - i &CenterDot; &rho; &OverBar; = I d &OverBar; c + i &CenterDot; I q &RightArrow; c + &xi; 1 &CenterDot; | | &Phi; &RightArrow; rc | | - - - ( 130 )
&Psi; &RightArrow; 2 c &CenterDot; e - i &CenterDot; &rho; &OverBar; = I d &OverBar; c + i &CenterDot; I q &RightArrow; c + &xi; 2 &CenterDot; | | &Phi; &RightArrow; rc | | - - - ( 140 )
一旦在步骤64中已经建立设置点和初始特征向量,就借助于关系127计算控制向量。
在步骤66中,转换器14通过脉冲宽度调制在间隔T中连续地施加这个计算出的平均控制向量。
然后对于新间隔T重复步骤42至66。
这种控制系统和方法的多种其它实施例是可行的。例如,初始磁通向量可被测量,而不是借助于观察器进行估计。
这里在其中实现产生电流设置点和控制向量的准确分析解的特定情形下,描述了该方法和系统。可选择地,作为预计的特定用途的函数,可省去这些准确分析解的可忽略项。例如,如果T与电机的时间常数相比很小,则可以把关系(123)的平方根作为|η0|处理。
在其中电机由转换器控制的特定情形下,已经描述了系统10。然而,该描述同样适用于由能够产生所计算出的控制电压的开关模式控制器、放大器或整流器所控制的电机。
最后,如在引言中指示的那样,如果描述电机操作的状态方程系统不是线性的,则有可能把它在多个操作点周围线性化,并且把以上描述的方法和系统应用于每个线性化系统。

Claims (7)

1. 一种控制异步旋转电机的直进式控制方法,该电机的磁通取决于其激励电流,该方法包括在时间间隔T上连续地施加一个平均控制向量或施加其平均值,以在间隔T的末端处得到与预定机械力设置点相等的、由所述电机产生的机械力,该方法的特征在于它包括:
a)建立设置点电流的步骤(60),该设置点电流适于在间隔T的末端处实现机械力设置点和设置点磁通,并且基于在所述设置点电流、设置点磁通及在间隔T的开始处的初始电流和磁通值之间的关系而被建立;和
b)计算要在间隔T期间连续施加的所述平均控制向量的步骤(64),以在间隔T的末端处得到与在步骤a)期间建立的设置点电流相等的电流。
2. 根据权利要求1所述的方法,其特征在于,设置点电流是如下复数关系的解:
( a 2 &CenterDot; p 11 - a 1 &CenterDot; p 21 ) &CenterDot; I &RightArrow; c + ( a 2 &CenterDot; p 12 - a 1 &CenterDot; p 22 ) &CenterDot; &Phi; &RightArrow; c =
( a 2 &CenterDot; e 1 &CenterDot; p 11 - a 1 &CenterDot; e 2 &CenterDot; p 21 ) &CenterDot; I &RightArrow; 0 + ( a 2 &CenterDot; e 1 &CenterDot; p 12 - a 1 &CenterDot; e 2 &CenterDot; p 22 ) &CenterDot; &Phi; &RightArrow; 0
其中:
·pij是由如下关系A=P·D·P-1定义的特征值的矩阵P的逆矩阵的系数,其中D是包括矩阵A的特征值μi的对角矩阵,P-1是矩阵P的逆矩阵,及矩阵A是在模拟所述电机操作的状态方程的线性系统中所述电机的自由特性的矩阵;
·系数ai是由如下关系C=D-1·(eD·T-I)·(P-1·B)定义的矩阵C的系数,其中矩阵I是单位矩阵,并且矩阵B是在状态方程的线性系统中的所述电机的控制矩阵;
·向量
Figure C2006100848510002C4
分别代表在间隔T的开始处所述电机的电流的和磁通的初始状态;
·向量
Figure C2006100848510002C5
Figure C2006100848510002C6
分别是设置点电流向量和设置点磁通向量;及
·系数ei由如下关系 e i = e &mu; i &CenterDot; T 定义。
3. 根据权利要求1或权利要求2所述的方法,其特征在于,所述平均控制向量是所述电机的状态方程的如下离散系统的准确解:
X &RightArrow; c = F ( T ) &CenterDot; X &RightArrow; 0 + G ( T ) &CenterDot; V &RightArrow; &OverBar;
其中:
·向量是定义在间隔T的开始处的所述电机的电流和磁通的初始状态的状态向量;
·向量
Figure C2006100848510003C3
是定义所述电机的设置点电流和设置点磁通的设置点状态向量;
·向量
Figure C2006100848510003C4
是所述电机的平均控制向量;
·矩阵F(T)是由关系F=eA·T定义的电机的转换矩阵,其中A是在模拟所述电机操作的状态方程的线性系统中的所述电机的自由特性的矩阵;及
·矩阵G(T)是所述电机的控制矩阵,该控制矩阵的值是间隔T的值的函数。
4. 根据权利要求2所述的方法,其特征在于,所述平均控制向量是如下关系的结果:
V &RightArrow; &OverBar; = e i &CenterDot; &theta; 0 &CenterDot; e &mu; i &CenterDot; T [ p 11 ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + p 12 &CenterDot; &Phi; 0 d 0 ] &CenterDot; [ p 21 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 22 &CenterDot; &Phi; cd ] - e &mu; i &CenterDot; T &CenterDot; [ p 21 &CenterDot; ( I 0 d 0 + i &CenterDot; I 0 q 0 ) + p 22 &CenterDot; &Phi; 0 d 0 ] &CenterDot; [ p 11 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 12 &CenterDot; &Phi; cd ] a 2 &CenterDot; [ p 11 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 12 &CenterDot; &Phi; cd ] - a 1 &CenterDot; [ p 21 &CenterDot; ( I cd + i &CenterDot; I cq ) + p 22 &CenterDot; &Phi; cd ]
其中:
·θ0是初始轴线系统的角位置,该系统的横坐标轴在间隔T的开始处与所述电机的初始磁通向量对准;
·系数pij是由关系A=P·D·P-1定义的特征值的矩阵P的逆矩阵的系数,其中D是包括矩阵A的特征值μi的对角矩阵,P-1是矩阵P的逆矩阵,及矩阵A是在模拟所述电机操作的状态方程的线性系统中所述电机的自由特性的矩阵;
Figure C2006100848510003C6
是在初始轴线系统中在间隔T的开始处初始电流向量的坐标,该系统的横坐标轴与初始磁通向量的方向对准;
Figure C2006100848510003C7
是在初始轴线系统中初始磁通向量的横坐标;
·Icd、Icq是在设置点轴线系统中在间隔T的末端处设置点电流向量的坐标,该系统的横坐标轴与设置点磁通向量的方向对准;及
·Φcd是在设置点轴线系统中设置点磁通向量的横坐标。
5. 根据权利要求2所述的方法,用于具有定子和转子的异步旋转电机,其特征在于,模拟所述电机操作的状态方程的线性系统的状态向量由定子电流向量和转子磁通向量形成。
6. 根据权利要求4所述的方法,其中,在其横坐标轴与设置点磁通向量的方向对准的设置点轴线系统中,设置点电流在间隔T的末端处具有坐标Icd和Icq;并且由坐标Icq和设置点磁通来建立坐标Icd
7. 一种用来控制异步旋转电机的直进式控制单元,该电机的磁通基于其激励电流被控制,其特征在于,所述单元包括:
·计算模块(30),用来建立设置点电流,该设置点电流适于在间隔T的末端处实现机械力设置点和设置点磁通,并且基于在所述设置点电流、设置点磁通及在间隔T的开始处的初始电流和磁通值之间的关系而被建立;和
·用来计算要在间隔T期间连续施加,以在间隔T的末端处得到与由所述计算模块(30)建立的设置点电流相等的电流的平均控制向量的模块(32)。
CNB2006100848514A 2005-05-23 2006-05-23 用于异步旋转电气机械的一种直进式控制方法和控制单元和一种用于其的存储介质 Expired - Fee Related CN100424989C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR0505137 2005-05-23
FR0505137A FR2886075B1 (fr) 2005-05-23 2005-05-23 Procede et unite de commande a reponse pile d'une machine asynchrone, support d'enregistrement pour ce procede

Publications (2)

Publication Number Publication Date
CN1870419A CN1870419A (zh) 2006-11-29
CN100424989C true CN100424989C (zh) 2008-10-08

Family

ID=35583393

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2006100848514A Expired - Fee Related CN100424989C (zh) 2005-05-23 2006-05-23 用于异步旋转电气机械的一种直进式控制方法和控制单元和一种用于其的存储介质

Country Status (6)

Country Link
US (1) US7298112B2 (zh)
EP (1) EP1727272B1 (zh)
CN (1) CN100424989C (zh)
CA (1) CA2547023C (zh)
ES (1) ES2446110T3 (zh)
FR (1) FR2886075B1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2916585B1 (fr) * 2007-05-25 2009-08-21 Alstom Transport Sa Procede de regulation d'une tension ou d'un courant d'un filtre rlc, support d'enregistrement et vehicules pour ce procede.
FR3054754B1 (fr) * 2016-08-01 2020-02-21 Supergrid Institute Procede de controle d'un convertisseur modulaire multi-niveaux
CN113359454B (zh) * 2021-06-15 2022-11-22 中国科学院长春光学精密机械与物理研究所 一种提高多自由度运动系统控制精度的方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1045514A1 (fr) * 1999-04-13 2000-10-18 Alstom Procédé de régulation d'une machine tournante et dispositif d'alimentation d'une telle machine

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2614481B1 (fr) * 1987-02-13 1990-08-31 Pk I Procede de commande d'un moteur asynchrone et entrainement electrique mettant ce procede en application
US4894599A (en) * 1988-03-31 1990-01-16 International Business Machines Corporation deadbeat control of disk drive actuator
US5334923A (en) * 1990-10-01 1994-08-02 Wisconsin Alumni Research Foundation Motor torque control method and apparatus
ATE181472T1 (de) * 1994-03-23 1999-07-15 Abb Daimler Benz Transp Verfahren zur drehmomentregelung einer asynchronmaschine
JPH11235099A (ja) * 1997-12-08 1999-08-27 Meidensha Corp 誘導機のベクトル制御方式
JP4475368B2 (ja) * 2000-03-10 2010-06-09 富士電機システムズ株式会社 速度センサレスベクトル制御装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1045514A1 (fr) * 1999-04-13 2000-10-18 Alstom Procédé de régulation d'une machine tournante et dispositif d'alimentation d'une telle machine

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A current control technique for voltage-fed induction motordrives. Sheng-Ming Yang, Chen-Haur Lee.Industrial Electronics Society, 1999. IECON '99 Proceedings. The 25th Annual Conference of the IEEE,Vol.3 . 1999
A current control technique for voltage-fed induction motordrives. Sheng-Ming Yang, Chen-Haur Lee.Industrial Electronics Society, 1999. IECON '99 Proceedings. The 25th Annual Conference of the IEEE,Vol.3 . 1999 *

Also Published As

Publication number Publication date
US20060273754A1 (en) 2006-12-07
EP1727272A3 (fr) 2008-03-05
FR2886075B1 (fr) 2007-09-14
CA2547023A1 (fr) 2006-11-23
CA2547023C (fr) 2009-08-18
ES2446110T3 (es) 2014-03-06
US7298112B2 (en) 2007-11-20
EP1727272B1 (fr) 2014-01-15
EP1727272A2 (fr) 2006-11-29
FR2886075A1 (fr) 2006-11-24
CN1870419A (zh) 2006-11-29

Similar Documents

Publication Publication Date Title
US11146197B2 (en) Method of controlling a permanent magnent synchronous motor and a motor circuit
CN105075104A (zh) 用于确定感应电机的励磁曲线和转子电阻的系统以及制造其的方法
CN110112974B (zh) 电机控制方法、控制器、存储介质及电机驱动系统
CN109946978B (zh) 一种考虑延时环节的伺服系统分数阶模型辨识方法
CN100424989C (zh) 用于异步旋转电气机械的一种直进式控制方法和控制单元和一种用于其的存储介质
JP2012222890A (ja) モータ制御装置
Wallscheid et al. Design and empirical identification of a lumped parameter thermal network for permanent magnet synchronous motors with physically motivated constraints
Bramerdorfer et al. Measurement-based nonlinear modeling of PMSMs
Liu et al. Kriging assisted on-line torque calculation for brushless DC motors used in electric vehicles
Wang et al. Sliding mode variable structure I/O feedback linearization design for the speed control of PMSM with load torque observer
Ounnas et al. A combined methodology of H∞ fuzzy tracking control and virtual reference model for a PMSM
Gevorkov et al. Hardware-in-the-loop simulator of a flow control system for centrifugal pumps
Gaona et al. Sensitivity analysis of a permanent magnet temperature observer for PM synchronous machines using the Monte Carlo method
Salah-eddine et al. Microcontroller based data acquisition and system identification of a dc servo motor using arx, armax, oe, and bj models
Lešić et al. State and parameter estimation for field-oriented control of induction machine based on unscented Kalman filter
Vu et al. Observer-based LQR for wind energy conversion systems with single measurement
Susdorf et al. Synthesis of power-efficient series motor control
Wójcik et al. Application of iterative learning control for ripple torque compensation in PMSM drive
Ab Rahman et al. System identification for a mathematical model of DC motor system
Okada et al. Performance Improvement of Magnet Temperature Estimation using Kernel Method based Non-Linear Parameter Estimator for Variable leakage flux IPMSMs
Lodha et al. Robust control of DC drive using sliding mode control
Yoo Actuator performance envelope through nonlinear test bed
Lachheb et al. Nonlinear Analytic Modeling for Novel Linear Variable Reluctance Motors
Gutiérrez-Villalobos et al. Hybrid artificial neural network for induction motor parameter estimation
Harkat et al. Improved control of permanent magnet synchronous motors using adaptive nonlinear control with Deadbeat Observer

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: ALSTOM TRANSPORT TECHNOLOGY CO., LTD.

Free format text: FORMER OWNER: ALSTOM TRANSP SA

Effective date: 20150527

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20150527

Address after: French Levallois Perret

Patentee after: Alstom transport scientific & technical corporation

Address before: French Levallois Perret

Patentee before: Alstom Transport S. A.

CP02 Change in the address of a patent holder

Address after: France 93400 Darren Albert Saint Ouen Street No. 48

Patentee after: Alstom transport scientific & technical corporation

Address before: French Levallois Perret

Patentee before: Alstom transport scientific & technical corporation

CP02 Change in the address of a patent holder
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20081008

Termination date: 20210523

CF01 Termination of patent right due to non-payment of annual fee