专利名称：双凸极电机无位置传感器控制方法
技术领域：
双凸极电机无位置传感器控制方法，属特种电机数字控制方法。
背景技术：
对于双凸极电机驱动系统(DSM)而言，实时而准确的转子位置信息是其可靠运行的必要前提。目前双凸极电机系统中，一般采用轴位置传感器(见图1)或者其它类型的探测式位置检测器来获得位置信息(见图2)，这不仅会增大体积，提高系统成本和复杂程度，更重要的是会降低DSM系统结构的坚固性，并且维修困难，影响整个系统的可靠运行，尤其限制了其在环境较恶劣场合的应用。因此如何让它去掉位置检测器，直接利用电机的电压和电流信息间接确定转子位置，从而使系统结构更加坚固，运行更加可靠、高效，成本更加低廉，无疑是一个很有潜力的研究方向。
迄今为止，国内外对这一领域的研究尚属起步阶段，研究对象多集中于电机结构相类似的开关磁阻电机无位置传感器。其研究的思路有两种(1)脉冲注入定位法利用空闲相，人为地注入低幅高频的模拟测试信号从而产生需要的电流等信息以得到位置信息。脉冲注入法的算法尽管相对比较简单，但基于高频脉冲的输入使其不免有着内在固有的速度限制，而且测试电流可能带来负转矩，其对整个系统出力和效率的影响也是很大的不足；(2)反电势法通过检测端电压的过零点，得到电机的位置信号。但是这种方法需要辅助起动设施，需要进行相位补偿，还存在计数值转速电流限幅的匹配问题。因为这些缺点的存在，最后以上的方法都没有能够很好地实现。
如何实时精确获得转子位置角度信息一直是双凸极无位置传感器技术的难题。国内外尚没有有效的双凸极电机无位置传感器技术。
发明内容
本发明的目的是提供一种直接利用电机的电枢电流和电压信息间接地确定电机转子位置，实现电机无位置传感器(见图3)，达到电机系统结构简单、运行可靠、体积小、坚固、效率高、成本低的新型双凸极电机无位置传感器的控制方法。
本发明的双凸极电机无位置传感器控制方法，直接利用电机的电枢电流和电压信息间接确定电机转子位置，实现双凸极电机无位置传感器的运行。具体方法包括1.通过实测电机的电枢电流和电压，求得在不同电流下随电机转子位置变化的电感值；2.将同一电流下电机转子位置角作为自变量，电感值作为因变量，通过MATLAB软件进行曲线拟合，绘制出同一电流下电感值随电机转子位置角度变化的曲线；3.利用MATLAB软件中的polyfit函数实现电感曲线的多项式拟合，以得到电感值随电机转子位置角度变化的多项式函数表达式
Lp(θ)＝a0*θ5+a1*θ4+a2*θ3+a3*θ2+a4*θ+a5；4.拟合出不同电流下电感值随电机转子位置角度变化的曲线族；5.利用电感曲线最小值在不同电流下几乎不变的特点，用电流为0时的电感拟合多项式减去电感的最小值后乘以与电流有关的比值函数(f(i))，再加上这个电感最小值，构造电感值随电流和电机转子位置角度变化的三维函数Lp(θ，i)＝(a0*θ5+a1*θ4+a2*θ3+a3*θ2+a4*θ+a5-L0min)*f(i)+L0min；6.利用数字信号处理器(DSP)或单片机芯片将采样到的电流电压信息进行处理得到参考电感值的大小，再把所得参考电感值与第5条所述的函数值作比较，确定不同时刻的转子位置信息；7.根据双凸极电机的运行机理，不同的电机转子位置信息对应于不同的电感区间，在电感上升区相绕组通入正电流，在电感下降区相绕组通入负电流，实现双凸极电机的无位置传感器的运行。
本发明的双凸极电机无位置传感器控制方法，不需要人为地注入低幅高频的模拟测试信号，也不存在相位补偿，以及电流转速匹配的问题，而是直接以电机运行时的电流电压信息为基础，根据电机的三维电感曲线反演求出位置信息，利用DSP或单片机芯片实现间接位置检测，可使双凸极电机系统更加简单体积更小结构更坚固运行可靠效率高等优点，此方法对无位置传感器双凸极电机的应用有较大的实用价值。
图1是采用轴位置传感器的双凸极电机截面图。
图2是采用轴位置传感器获得的三相位置信号示意图。
图3是本发明的无位置传感器双凸极电机的结构示意图。
图4是计算机采样电流电压示意图。图4中符号名称U——电压，i——电流，T——周期，k——周期数图5是电感曲线族(电流范围±50A)。
图6是多个电感周期下电感曲线族。
图7是近似电感曲线族。
图8是电感最大值比拟合曲线。
图9是相绕组自感的三维曲面图。
图10是电感曲面中取出的电感曲线的误差示意图。
图11是电流范围为±200A时由构造的三维函数得到的电感曲面图。
图12是电机转子位置获取的计算机流程图。
图13是电励磁双凸极电机无位置传感器运行机理图。
具体实施例方式
1实际电感值估算1.1磁链估算以一台12/8极电磁式双凸极电机为例，本专利技术首先利用电机运行时实测电压和电流来估算磁链值。
双凸极电机一相电压方程为u=ir+dψdt---(1)]]>由式(1)解出一相绕组磁链表达式ψ=∫0t(u-ir)dt+ψ(0)---(2)]]>如果已知从时刻0到时刻t间每一时刻的电压和电流值i以及时刻0时的初始磁链ψ(0)，就可以积分计算出绕组当前时刻的实际磁链。对(2)式进行离散得到ψj=ψ(0)+Σk=0N[u(k)-ij(k)r(k)]·T]]>⇒ψj(N)=ψj(N-1)+[u(N)-ij(N)r(N)]·T---(3)]]>式中，T为离散积分磁链的计算周期，u(k)、i(k)、r(k)均为第k个离散周期内绕组的电压、电流和电阻值。由式(3)可知1个周期T中磁链的变化量Δψ为Δψ＝[u(N)-i(N)r(N)]T(4)对电流采用PWM控制，且PWM斩波周期采样周期和磁链估算的离散周期T相同且保持同步，因而电源电压和绕组电流在第k个周期内的离散值就是该采用周期内的采样值Ud(k)、i(k)。若当前斩波周期的占空比为α，则如图4所示在第k个周期T里绕组电流不变为i(k)，绕组电压在开通时段αT内为电压电压Ud(k)，在零电压续流斩波时段(1-α)T内绕组电压为零，假设绕组电阻为常数r，且不随时间变化，则在第k个周期内的磁链增量为Δψ(k)＝[Ud(k)-i(k)r]α(k)T+
[1-α(k)]T(5)＝[Ud(k)α(k)-i(k)r]T令Δψ1(k)＝Ud(k)α(k)T，Δψ2(k)＝i(k)rT得Δψ(k)＝Δψ1(k)-Δψ2(k) (6)分别算出Δψ1(k)和Δψ2(k)两项，再代入式(6)得到第k个离散周期T内的磁链增量，最后得到t＝kT时电机当前导通相的磁链估算值，如式(7)所示ψ(k)＝ψ(k-1)+Δψ(k) (7)当到达最后一个采样周期(k＝N)，可利用单片机或DSP记录下磁链随角度变化的离散数据值。不断提高测试电压，得到不断提高的稳态测试电流，进而可以得到磁链ψ随角度θ和电流i变化而变化的离散数据值。
1.2电感估算由磁链—电感公式ψ＝Li (8)变形可得L＝ψ/i (9)因此，在获得磁链的离散数据值之后，由(9)可以得到电感随位置角度和电流变化的离散数据值表。
2单条非线性电感曲线拟合的实现方法由估算出的电感值数据表，选择样条插值的方法对电感数据进行处理。先将同一电流下转子位置角作为自变量，电感值作为应变量，经过曲线拟合得到同一电流下电感值随转子位置角度变化的曲线。利用样条插值的方法得到不同电流值下电感值的拟合曲线族(电流分别取-50、-40、-30、-20、-10、0、10、20、30、40、50A)，见图5所示。最后在样条插值的基础上，利用MATLAB软件的polyfit函数可以实现电感曲线的多项式拟合。选择5次拟合后的结果Lp(θ)＝a0*θ5+a1*θ4+a2*θ3+a3*θ2+a4*θ+a5(10)其中，Lp(θ)表示相绕组自感，θ表示转子位置，a0……a5为拟合多项式的系数。(10)式即为电感随转子位置角度变化的函数。
3构造函数原理介绍本专利在对电感拟合曲线进行研究的基础上，利用电感曲线最小值几乎不变的特点，把电流与位置角度以及电感的关系用数学的方法联系起来，构造出电感随角度位置以及电流变化的高次函数。
由式(10)可以得到电感随角度变化的函数关系式，同时得到电感曲线(电流从-50A到+50A)，从时间轴上选取多个电感值周期，见图6。这一系列的电感曲线可以近似抽象地看做是若干个等腰梯形连接而成的，见图7(d)。选取图7(d)曲线族中任意一条曲线(如图7(a)所示)，把它向下平移，使曲线的最小值与X轴重合，可以得到如图7(b)所示的曲线。这条曲线的表达式乘以一系列的常数就可以得到形如图7(c)所示的曲线族。图7(c)再向上平移就可以得到形如图7(d)的曲线族。
因此只要知道电感曲线族中任意一条曲线的函数表达式，把它减去电感最小值，再乘以一系列系数，最后在加上电感最小值就可以近似得到形如图6的曲线族。并且由图7(c)可以看出当电流均匀变化时，每条电感曲线之间近似存在一定的比值关系，且这个比值关系主要由电感的最大值决定。这就是利用图形特点构造三维函数的原理基础。
4电感最大值比曲线拟合由于双凸极电机相绕组通入正电流和负电流对电感的影响不一样，一个起增磁作用，另一个起去磁作用，因此正负电流对于电感的影响应该分别考虑。本文选择电流为0时的电感最大值减去电感最小值后的值作为基准，用其它电流情况下电感曲线的最大值减去最小值后的值与之比较，得到一系列的比值pn，见式(11)pn＝(Lnmax-Lmin)/(L0max-Lmin) (11)其中，pn表示不同电流下电感最大值减去电感最小值后的比值，Lnmax表示不同电流下电感的最大值，Lmin表示电感最小值(电感最小值在不同电流下几乎不变)，L0max表示电流为0时电感的最大值，n＝1，2，3……。这一系列比值可以拟合成一条曲线(见图8)，这条连续的曲线反应了不同电流下电感曲线之间近似的比值变化关系。由MATLAB的曲线拟合函数可以得到上图中曲线的表达式f(i)＝-0.0035i+1.0005(12)5.构造三维电感函数用电流为0时的电感拟合多项式(10)减去电感的最小值后乘以式(12)，再加上这个电流最小值，可以近似构造出的电感随电流和位置角度变化的连续函数Lp(θ，i)＝(a0*θ5+a1*θ4+a2*θ3+a3*θ2+a4*θ+a5-L0min)*f(i))+L0min(13)其中，L0min为电枢电流为0时的电感曲线最小值。由式(13)可以得到电感随电流和角度变化的三维曲面，如图9所示(电流范围-50A～50A)。
把估算电感值与构造的三维函数所得电感曲线进行比较，见图10。平滑曲线表示的是由构造的函数所得的电感曲线族(电流范围±50A)，“*”表示的是估算电感值。由图10可以看出构造函数法所得电感值与实际估算电感值非常接近，曲线甚至几乎重合，误差与实际估算值相比很小，说明该模型有较高的精确性，能够较好地反映实际电感的变化。本专利技术方案中构造的三维函数在电感最大值与电感最小值附近的精度最高，并且结构简单(只包含加法和乘法运算)，编程简单，因此非常合适用来确定转子位置角度，实现无位置检测。当电流范围再扩大，如扩大到±200A，构造函数得到的电感三维曲面图见图11。
6转子位置角度估算利用DSP或单片机将采样到的电流电压信息进行处理得到参考电感大小，再把所得参考电感值与之前构造的三维函数值作比较，当两者差的绝对值在一定的误差范围之内，记录下此时的转子位置位置θ，就可以确定不同时刻转子的位置角度。位置角度求取的流程框图见图12。其中θ表示转子位置角度(本例子中θ范围0～45°)，L1表示由实测电压和电流估算出的电感参考值，L2表示由三维函数计算的电感值，W表示L1与L2之差，e表示一无穷小量。
7无位置传感器的实现由双凸极电机运行机理，不同的转子位置θ分别对应于电感曲线的上升区(图13中I区域)和下降区(图13中III区域)，电机的运行需要在电感上升区相绕组通入正电流，在电感下降区相绕组通入负电流。本专利技术实现对位置信息的获取后，对软件记录下的θ值进行分类，不同的θ值分别对应于不同的电感上升区和电感下降区，如图13所示，在θ1～θ2时绕组通入正点流ia，在θ3～θ4时相绕组通入负电流-ia(由于电感的阻碍作用，提前换向可以增加电流的出力，即可以在θ2～θ4时相绕组通入负电流)，从而实现双凸极电机无位置传感器运行。此无位置传感器方案不需要人为地注入低幅高频的模拟测试信号，而是直接以电机运行时的电流电压信息为基础，根据电机的三维电感曲线反演求出位置信息，利用DSP或单片机芯片可以实现间接位置检测。此方案对无位置传感器双凸极电机的应用具有较大的实用价值。
权利要求
1.一种双凸极电机无位置传感器控制方法，其特征在于，直接利用电机的电枢电流和电压信息间接确定电机转子位置，实现双凸极电机无位置传感器的运行，具体方法包括(1)通过实测电机的电枢电流和电压，求得在不同电流下随电机转子位置变化的电感值；(2)将同一电流下电机转子位置角作为自变量，电感值作为因变量，通过MATLAB软件进行曲线拟合，绘制出同一电流下电感值随电机转子位置角度变化的曲线；(3)利用MATLAB软件中的polyfit函数实现电感曲线的多项式拟合，以得到电感值随电机转子位置角度变化的多项式函数表达式Lp(θ)＝a0*θ5+a1*θ4+a2*θ3+a3*θ2+a4*θ+a5；(4)拟合出不同电流下电感值随电机转子位置角度变化的曲线族；(5)利用电感曲线最小值在不同电流下几乎不变的特点，用电流为0时的电感拟合多项式减去电感的最小值后乘以与电流有关的比值函数(f(i))，再加上这个电感最小值，构造电感值随电流和电机转子位置角度变化的三维函数Lp(θ，i)＝(a0*θ5+a1*θ4+a2*θ3+a3*θ2+a4*θ+a5-L0min)*f(i)+L0min；(6)利用数字信号处理器(DSP)或单片机芯片将采样到的电流电压信息进行处理得到参考电感值的大小，再把所得参考电感值与第5条所述的函数值作比较，确定不同时刻的转子位置信息；(7)根据双凸极电机的运行机理，不同的电机转子位置信息对应于不同的电感区间，在电感上升区相绕组通入正电流，在电感下降区相绕组通入负电流，实现双凸极电机的无位置传感器的运行。
全文摘要
一种双凸极电机无位置传感器控制方法，属特种电机数字控制方法。本方法是直接利用电机的电枢电流和电压信息间接地确定电机转子位置，实现电机无位置传感器的可靠运行。不需要人为地注入低幅高频的模拟测试信号，也不存在相位补偿及电流转速匹配问题，而是直接以电机运行时的电流电压信息，根据电机的三维电感曲线反演求出转子位置信息，利用数字信号处理器或单片机芯片实现间接位置检测，使双凸极电机系统更简单、体积更小、结构更坚固、运行更可靠、效率更高等优点。此方法对无位置传感器双凸极电机的应用有较大的实用价值。
文档编号H02P6/14GK1725625SQ200510040720
公开日2006年1月25日 申请日期2005年6月23日 优先权日2005年6月23日
发明者张乐, 周波 申请人:南京航空航天大学