具体实施方式
在参照图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)写成:
通过把两个左手项乘以P-1,在几次变换后得到:
机械的这种新分立状态表示因为其简单性是有益的:eD·T和D-1·(eD·T-I)是用于给定机械的对角矩阵,并且非常容易由特征值μi形式计算。事实上,它们的相应对角分别由如下元素组成:
这个公式表示状态向量到特征向量空间,即状态“特征向量” 以及控制“特征矩阵”:[P-1·B]的空间中的投影。作为特征化系统的特性的特征值,这些状态特征向量具有特定的物理意义,因为它们的坐标,即“特征值”,与控制定律分离,因为它们表示在特征向量基中。这同样反映诸个矩阵的两个是对角的事实。
变量的如下变化简化状态方程(7):
控制向量的计算然后包括:对于不可得到的参数填入来自测量和/或观察器的初始特征状态
填入设置点特征状态
作为要实现的目标的函数;及最后计算能够在选择时间间隔T的末端处实现设置点特征状态的控制向量。
对于旋转电气机械,由于目标是得到机械力,所以实际上是,选择由机械的激励电流的状态和机械的磁通的状态组成的状态向量
事实上,通过简单的向量乘积,在电流向量和磁通向量之间的幅值和角度把力目标变换成关于在选择的轴线系统中的那些向量的坐标,以简化其表达式。这里轴线系统是其横坐标轴线与磁通向量对准的正交轴线系统。在这种情况下,Laplace定律向量乘积由磁通的模量和在相同时刻电流向量到纵坐标轴线上的投影的模量的乘积表示,忽略取决于机械的常数。况且,磁通的模量和在相同时刻到横坐标轴线上的电流向量的投影的模量的乘积定义电流和磁通向量的标量乘积,并因此定义在Park表示中的机械的磁能,忽略取决于机械的乘法常数。
矩阵Γ0定义为是在初始轴线系统中电流向量的投影矩阵,该系统的横坐标轴与初始磁通向量对准。这里,“初始”是指标记间隔T开始的时刻tn。这个矩阵能以两个矩阵Γc和ΓT乘积的形式写出,因为在第一个情况下,把电流向量投影到其横坐标轴与设置点磁通向量对准的设置点轴线系统中,并且在第二个情况下,把该轴线系统的电流向量向初始轴线系统投影。这种变换ΓT当然取决于间隔T的值,由于它从设置点轴线系统到初始轴线系统(假定系统对于这种形式的表达是因果的)。因此满足如下关系:
ΓT·Γc=Γ0 (12)
把关系(11)和(12)结合产生:
这种形式表示设置点状态特征向量
和初始状态特征向量
分别表示在设置点轴线系统中和在初始轴线系统中,这意味着可非常简单地表达向量和标量乘积。
矩阵(D、P、P
-1、B)取决于所述机械,并且可以用分析方法或数值方法计算。控制间隔T基于其它标准选择。向量
取决于在间隔T的开始处机械的初始状态,并且向量
是设置点向量。对于每个控制步骤这两个向量
与三个投影矩阵:(Γ
c、Γ
0、Γ
T)一起必须计算。Γ
T依据应用可能在未知变量中确定为测量和/或设置点值的函数。
I.2-设置点电流的表达
这里设置点向量被定义,以能够同时实现例如以主转矩设置点C
c的形式表达的机械力设置点F
mc、和在间隔T的末端处要实现的设置点磁通
在由电流向量
和磁场向量
定义的平面中,具有至少四维的状态向量选择成定义机械的状态:
应用于设置点的这个状态向量基于磁性状态和机械力被定义。
为了简单化在各向同性异步机械的情况下的计算形式,电流和磁通向量可以由复数代表,每个在轴线的直接正交系统中具有两个坐标,其中实数横坐标轴在磁场向量的方向上,并因此在磁通向量的方向上。
系统(13)保持有四阶,但现在在复数平面中由二阶系统表示。
μ1和μ2定义状态方程的离散系统的特性的复数矩阵的两个复数特征值;复数控制矩阵和复数逆通过矩阵用这些维可写成如下:
使Pjk为μ1、μ2的函数(16)
控制向量是一阶的和复数的,并且由关系 表达。便利的是,选择相对于控制单元固定的轴线系统的坐标(α,β)。
借助于以上定义,C=D-1·(eD·T-I)[P-1·B]成为:
让我们使用这些定义细化状态特征向量的方程系统:
其中:
让:
方程系统成为:
它然后划分成两个复数方程:
这可等效地写成:
然后在方程(22)和(23)之间消去复数控制向量:
这种关系使得显然的是:设置点电流向量被链接到设置点磁通向量上,并且都取决于这个方程的第二项的已知初始状态。
取决于所述应用的系数(a2·p11-a1·p21)是非零的,因为旋转电气机械的磁通取决于机械的激励电流。因此有可能写出如下关系:
(27)
注意,如果系数(a2·p11-a1·p21)是零,则跨过初始状态的设置点电流和设置点磁通的相互依赖性完全消失,这在实际中意味着不可能由电流向量控制磁通向量。磁通因此或者是恒定的或者以某种其它方式控制,并且这里描述的直进式控制方法不适用。
让:
得到关系(27)的更简明形式:
其中电流和磁通向量在设置点轴线系统中由复数表达式表示。
其中Icd是磁化电流,并且Icq是与电机的磁通相结合产生机械力的电流。它然后有可能写成:
其中
·代表向量的标量乘积的实数标量;
·代表向量乘积的模量,它也是实数。我们因此可以写出:
我们现在计算:电流和磁通向量的标量乘积,它等于磁能Q,忽略常数kq;和相同向量的向量乘积,它等于机械力Fm,忽略常数km:
这种复数关系划分成两个实数方程:
(37)
并且,因为选择的轴线系统的意义:
(38)
其中Re{...}和Im{...}分别指示复数的实部和虚部。
然后有可能以两种不同的形式通过消去ε0从这两个方程提取Icd:
第一种形式:
ε0由第二方程计算,并且它由其在第一方程中的表达式替换:
然后:
(42)
应用的上下文消去Arcsin()的两个主要确定的一个,具有负余弦的那个保留。
第二种形式:
两个方程的两项被平方和求和,以由其数值替换cos2(ε0)+sin2(ε0)=1:
这种新关系被扩展以计算Icd作为二次方程的解:
(44)
平方引入了必须借助于所述应用的上下文消去的增根,只保留在根式前面具有负号的那个根。
进一步注意:
这个关系因此用来把机械力设置点F
mc和磁通模量设置点
变换成用于把电流向量投影到与磁通向量轴正交的轴上的设置点:
对于Icd的一个或另一个那么取决于:
·应用的参数:取决于特征值的a1、a2、p11、p12、p21、p21;
把磁化电流链接到主要设置点、参数及初始状态上的关系的所述两种形式因此能够实现设置点电流的第二坐标的计算、设置点电流向量到磁通向量上的投影。
或
(48)
总之,我们使用以前关系从(Fmc、)得到(Icd、Icq)和(Φcd,0),使我们能够建造投影到设置点轴线系统中的设置状态特征向量:
其来自关系:
以这种方式建立的设置点电流向量使我们在间隔T的末端处同时得到机械力设置点F
mc和磁通设置点
我们类似地可写出:
初始状态由测量得知。
I.3-建立用于异步电动机的矩阵Γ
T
的例子
因为控制向量这里在相对于定子固定的轴线系统中施加到三相旋转电动机上,所以选择由定子电流向量
和转子磁通向量
组成的状态向量,如在轴线正交系统中从定子看到的那样,该轴线正交系统通常叫做(α,β),在减小电机的磁极对数目之后相对于电机的定子固定。在这种基准正交轴线系统(α,β)中,Concordia变换用于从三相系统到两相系统的这种变化。α轴线可以在例如相位R的电极性方向上选择。
定子电流向量可借助于布置在电机外部的、并且与电机的相或一部分相串联的电流传感器而测量。转子磁通向量可从电机电压和电流得知,或者可被观察。这两个向量用来估计或计算初始状态。如有必要,它们也联合用来通过它们的向量乘积并由此通过它们的坐标以及磁能设置点来移动机械转矩设置点。
选择的状态向量因此在这种情况下是4维的。有四个电气自由度。它可以用选择成简化符号的两个向量或者用在轴线复数系统(α,β)中代表它们的复数互换地指示,β是纯虚轴,以通过把矩阵的秩从四减小到二(在复数平面中)来简化特征值的计算。
使用复数的第二个优点在于,它通过简单的指数复数标量而不是使用2维平方复数矩阵变换坐标轴的转动。
矩阵ΓT然后能以如下复数形式写出:
其中Δθp代表在间隔T期间磁通的预测角转动。
借助于状态变化变量的这种选择,可规定初始状态特征向量:
变量
代表投影到由相对于固定基准轴线系统(α,β)的变换
所限定的初始轴线系统中的状态特征向量的初始状态。初始值
和
由测量的电流和观察的磁通得知。
投影到基准的初始和设置点框架中的状态特征向量的一般方程成为:
(55)
这个系统能以两个复数方程的形式再次细化:
(56)
(57)
在这两个复数方程之间再次消去控制向量:
(58)
或者:
现在可计算在初始轴线系统与设置点轴线系统之间的变换。
是4阶(实数)或两阶(复数),并且已经只使用了两次:在旋转电机的情况下一次一般计算Icd和一次计算Δθp。
I.4-控制向量的计算
状态方程系统的唯一未知量现在是控制向量的两个实数坐标,它们是如下两个复数方程之一的解:
它们可通过现在消去ΓT简单地形成:
(63)
这是复数关系,该复数关系供给最后两个未知量:控制向量的坐标。
现在可得出结论:在各向同性旋转电机的情况下已经转化4维状态方程系统,并且我们已经从主设置点(F
mc,
)求出用于四个未知量:I
cd、Δθ
p及
的准确分析表达式。
部分II-对于异步旋转电机的应用
异步旋转电机是各向同性旋转电机。
在异步电机的情况下,在Concordia三相到两相变换之后,使用图1的等效电路,它把转子和定子的漏电感与电机的定子归到一起。这个等效电路具有简单表示的优点,把定子电流分离成其两个分量:产生转子磁通的电流和产生转矩的电流。
II.1-异步电机状态表示
使用对应电气微分方程,容易表明,连续状态方程的线性系统可以以如下形式写出:
并因此:
有:
·定子电阻: Rs
·转子电阻: Rr
·定子漏电感: ls
·转子漏电感: lr
·互感: Lm
·定子电感: Ls=Lm+ls
·转子电感: Lr=Lm+lr
·定子时间常数:
·转子时间常数:
·极对数目: NP
·扩散系数:
机械和电气角频率的定义和符号:
·转子的机械角频率:Ω
·极性机械角频率:ω=Np·Ω
·定子电气角频率:ωs
·转子电气角频率:ωr
·转差:
及:
λ=σ·Ls (70)
II.2-在特征向量基中状态方程系统的表达
给出矩阵A的特征值μ1的特征方程是:
μ2+(α+γ-i·ω)·μ+(γ-α·β·Lm)·(α-i·ω)=0(72)
让Δ指示方程的判别式:
Δ=(α-γ-i·ω)2+4·α·β·Lm·(α-i·ω) (73)
两个特征值因此由下式表示:
这里注意,特征值取决于极性机械角频率。它们因此必须由电机的模型、测量速度及极对数目计算。
ω=Np·Ω (76)
对角化矩阵A是:
与特征值相对应的特征向量 是方程的一个解:
对于第一特征向量,定义到在乘法常数内:
矩阵方程的第二行产生如下关系:
α·Lm·p11=(μ1+α-i·ω)·p21 (80)
它使我们能够选择:
p21=1 (82)
以类似方式,对于第二特征向量,定义到在乘法常数内:
p22=1
通过矩阵 从而:因此可以写A·P=P·D
或者,使用特征值:
并且因此可以写出其逆矩阵:
或者,使用特征值:
我们因此可写出:
A=P·D·P-1 (89)
对于使用特征值的符号,A成为:
使用由结合中间变量生成的如下关系:
(μ1+γ)·(μ2+γ)=α·β·Lm·(μ1+μ2+γ) (91)
使用如下减小变量:
概括成:
用F=P·eD·T·P-1替换F=eA·T,我们可现在写出:
其中:
得到特别简单的形式:
我们识别方法的一般形式:
其中:
状态“特征向量”的坐标可以通过用它们的向量坐标替换状态向量而使显得非常简单:
状态的坐标、初始和设置点“特征向量”现在是:
它通过设置产生如下关系:
当电机的参数已知时,对于给定机械速度ω和对于给定预测水平T计算在方程中涉及的各种减小参数是简单问题。初始状态特征向量
在轴线系统(α,β)中也由观察器所测量或估计的测量电流
和磁通
得知。因而,为了能够计算未知量:在轴线系统(α,β)中的控制向量
只要保持把设置点状态特征向量固定在预测水平T处。
II.3-建立设置点电流:
为了固定对于异步电机特定的转子磁通和电磁转矩设置点,现在必须由Laplace定律,固定产生磁能的电流和磁通向量的标量乘积、和它们产生转矩的向量乘积。
现在注意,这里提到的电流和磁通向量在同一时刻-在间隔T的末端考虑。此外我们准备计算待施加的控制向量,从而通过模型考虑到电气时间常数在间隔T的末端处同时得到这些设置点。计算因此是必须没有两个轴线的随动或去耦的准确计算。
对于这种计算,初始轴线系统(d
0,q
0)和设置点轴线系统
如在图2A和2B中指示的那样分别定义。为了更精确,轴线系统(d
0,q
0)是其横坐标轴与初始磁通向量
对准的正交轴线系统,并且轴线系统
是其横坐标轴与设置点磁通向量
对准的正交轴线系统。
因为两个设置点向量定位在水平T处,即在间隔T的末端处,所以在间隔T的末端处预测轴线系统
由此以新符号预测。其位置是瞬时位置。
在初始时刻的转子磁通向量的幅角ρ
0是已知的,由于在此时刻的磁通是已知的。在水平T处的转子磁通向量的幅角
是未知的。然而,磁通向量转动期间的间隔T是已知的,因为它按照其它考虑而选择的。让我们把在这个间隔期间的磁通向量的瞬时角速度ω
s当作新的末知量。由此我们推出:
我们把在状态特征向量之间的关系的两项乘以:
ΓT·Γc=Γ0 (109)
在这种情况下,从关系(108)选择:
它等效于反向转动,或等效于从轴线系统(α,β)到对于设置点的设置点轴线系统
和到对于初始状态的初始轴线系统(d
0,q
0)的坐标轴变化。
其中:
设置点向量现在作为要求转矩的函数和作为必需的磁通的模量、两个主设置点和原始状态向量的函数而容易计算,并且原始状态向量是已知的。
如下关系通过消去控制向量
从关系(111)至(114)导出:
(115)
第二项完全是已知的;我们按如下定义模量和幅角:
它使我们能够简化以上表达式:
这种表达式通过取实部和虚部解出:
并且消去ωs:
(121)
这是对于瞬态和稳定状态条件用来把设置点瞬态激励电流链接到设置点磁通和链接到设置点转矩上的所需的关系。
这种关系的另一种形式可以通过平方以及求和实数和虚数部分得到,这消去角度:
(122)
这个方程解出非常简单:
我们在这种特定情形下识别已经以一般方式建立的解的两种形式。
作为设置点转矩和设置点磁通模量的函数的异步电机的电流设置点
由如下关系给出:
有
II.4建立矩阵Γ
T
:
矩阵ΓT由如下关系给出:
II.5计算控制向量:
控制向量由如下关系给出:
我们识别与电机对应的控制向量的一般解。
现在可在其中初始向量被测量的轴线系统(α,β)中写出初始向量,留下在其中计算设置点向量的设置点系统中的设置点向量。
部分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期间建立的初始磁通向量
的估计计算幅角ρ
0.
在步骤56中计算初始特征向量
和
为了更准确,在步骤56中借助于如下关系计算系数
和
在与步骤50至56并行执行的步骤60期间,借助于关系(123)和(124)建立电流设置点。
在步骤62中建立特征向量
和
为了更准确,在步骤62中借助于如下关系计算因数
和
一旦在步骤64中已经建立设置点和初始特征向量,就借助于关系127计算控制向量。
在步骤66中,转换器14通过脉冲宽度调制在间隔T中连续地施加这个计算出的平均控制向量。
然后对于新间隔T重复步骤42至66。
这种控制系统和方法的多种其它实施例是可行的。例如,初始磁通向量可被测量,而不是借助于观察器进行估计。
这里在其中实现产生电流设置点和控制向量的准确分析解的特定情形下,描述了该方法和系统。可选择地,作为预计的特定用途的函数,可省去这些准确分析解的可忽略项。例如,如果T与电机的时间常数相比很小,则可以把关系(123)的平方根作为|η0|处理。
在其中电机由转换器控制的特定情形下,已经描述了系统10。然而,该描述同样适用于由能够产生所计算出的控制电压的开关模式控制器、放大器或整流器所控制的电机。
最后,如在引言中指示的那样,如果描述电机操作的状态方程系统不是线性的,则有可能把它在多个操作点周围线性化,并且把以上描述的方法和系统应用于每个线性化系统。