CN113688499B - 一种基于动静态测试的三相12/8极开关磁阻电机建模方法 - Google Patents

一种基于动静态测试的三相12/8极开关磁阻电机建模方法 Download PDF

Info

Publication number
CN113688499B
CN113688499B CN202110777173.4A CN202110777173A CN113688499B CN 113688499 B CN113688499 B CN 113688499B CN 202110777173 A CN202110777173 A CN 202110777173A CN 113688499 B CN113688499 B CN 113688499B
Authority
CN
China
Prior art keywords
phase
inductance
reluctance motor
unsaturated
theta
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.)
Active
Application number
CN202110777173.4A
Other languages
English (en)
Other versions
CN113688499A (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.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
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 Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN202110777173.4A priority Critical patent/CN113688499B/zh
Publication of CN113688499A publication Critical patent/CN113688499A/zh
Application granted granted Critical
Publication of CN113688499B publication Critical patent/CN113688499B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/02Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

本发明公开了一种基于动静态测试的三相12/8极开关磁阻电机建模方法。由于开关磁阻电机的非线性与饱和性,开关磁阻电机建模需要使用其磁链和转矩特性数据。使用转矩平衡法静态测试得到四条磁链曲线,在电机旋转时进行动态测试得到开关磁阻电机不饱和电感。基于以上测试得到的曲线进行拟合计算得到开关磁阻电机所有位置与电流下的磁链和转矩特性,完成建模。所述建模方法能够有效提升转矩平衡法的精确度,降低转子位置固定法的复杂度。

Description

一种基于动静态测试的三相12/8极开关磁阻电机建模方法
技术领域
本发明涉及一种基于动静态测试的三相12/8极开关磁阻电机建模方法,属于电机控制领域。
背景技术
开关磁阻电机结构简单,转子无永磁体及绕组,转矩由定、转子间气隙磁阻变化产生。调速范围宽、可靠性高、适用于恶劣环境等优点,使其具有广阔的应用前景。然而,双凸极结构及磁路饱和特性,导致很难通过常规电磁和物理特性推导获得其精确的非线性数学模型,因此需要使用开关磁阻电机磁链特性建立其数学模型。目前开关磁阻电机磁链特性测试方法有两种:直接法与间接法。直接法是使用固定在电机上的磁链传感器直接测量磁链特性,该方法成本高误差大。间接法是通过测量电压与电流间接计算磁链,测试方法可分为两种:转子位置固定法与转矩平衡法。其中,转子位置固定法需要转子固定装置,测量电流与电压可得到每个位置的磁链特性。转矩平衡法是先通过给单相通电将转子旋转至对其位置,然后利用开关磁阻电机的对称结构,保持合转矩相同,给各相通入不同脉冲电压,可测量特殊位置的磁链特性。综上所述,转子位置固定法与转矩平衡法均有明显的缺点:转子位置固定法需要设计测量平台,测量成本高,耗时长,难以操作;转矩平衡法可得到磁链特性过少,拟合结果误差大。
发明内容
针对转子位置固定法实现成本高、操作繁琐以及转矩平衡法建模误差较大的缺陷,本发明提出了一种基于动静态测试的开关磁阻电机建模方法,通过转矩平衡法静态测试获得的特殊位置磁链特性,以及动态测量获得的开关磁阻电机不饱和电感,根据相电感与不饱和电感的关系,计算得到开关磁阻电机所有位置的磁链特性。技术方案如下:
步骤一:使用转矩平衡法,利用开关磁阻电机在对齐位置与非对齐位置转矩为零的特点,通过给各相通入不同的脉冲电压,获得0°,7.5°,15°和22.5°四个特殊位置电感特性。
步骤二:开关磁阻电机低速运行时,给非导通相注入脉冲电压,记录相绕组的电流波形。根据母线电压、开关管压降、二极管压降以及各相电流斜率计算非饱和电感,非饱和电感的计算由式(1)求解,结合所记录的相应位置信号,再根据电感的对称性得到开关磁阻电机不饱和电感特性。
式中Vdc、VD、VT、iph分别表示直流母线电压、二极管压降、开关管压降、相电流。
步骤三:通过式(2)计算开关磁阻电机特殊角度,并将电感特性分为三个阶段。其中θ1表示开关磁阻电机的定子前沿与转子后沿重合时的转子位置,θ2表示开关磁阻电机的定子前沿与转子前沿重合时的转子位置,θa表示开关磁阻电机的定子与转子中心对齐时的转子位置,θhr为转子前沿与定子中心线相重合时的转子位置,βr为转子极弧,βs为定子极弧,Nr为转子极对数。
步骤四:在0°<θ<θ1位置,认为相电感等于不饱和电感。
步骤五:在θ1<θ<θhr位置,认为不饱和电感与相电感为线性关系,可使用一次函数对不饱和电感与相电感进行拟合,如式(3)。基于不饱和电感得到θ1到θhr位置相电感。
L(θ)=k(i)Lun(θ,i)+L0(i) (3)
步骤六:在θhr<θ<θa位置,认为不饱和电感与相电感为二次关系,可使用二次函数对不饱和电感与相电感进行拟合,如式(4)。完成不饱和电感与相电感拟合后,基于不饱和电感数据与二次函数求根公式反解相电感数据,得到θhr到θ2位置相电感。
Lun(θ)=a(i)L(θ,i)2+b(i)L(θ,i)+c(i) (4)
步骤七:整合上述步骤得到[0,θa]的电感特性,根据开关磁阻电机电感特性的对称性可以得到完整的电感特性,由式(5)计算得到磁链特性,由式(6)计算得到转矩特性。
ψk(θ,ik)=Lk(θ,ik)ik (5)
其中,ik,Tk,ψk分别表示第k相绕组电流、转矩和磁链。
步骤八:基于步骤七得到的开关磁阻电机磁链和转矩特性,完成开关磁阻电机数学模型建立。
本发明的有益效果:①弥补了转矩平衡法精度误差过大的缺陷,能够得到较为精确的开关磁阻电机磁链和转矩特性;②测量过程操作简单,耗时短,成本低,能够实现开关磁阻电机建模。
附图说明
图1为转矩平衡法得到的四条磁链特性图。
图2为不饱和电感图。
图3为不饱和电感与相电感的关系曲线图。
图4为计算磁链与实验测试对比图。
图5为计算转矩与实验测试对比图。
图6为转速1560rpm,开关角5°,关断角15°电动状态下,相电流对比图。
图7为转速1000rpm,开关角30°,关断角40°发电状态下,相电流对比图。
具体实施方式
以下结合附图和具体实例,对本发明的技术方案进行详细说明。实例所用电机为一个1kW三相12/8极开关磁阻电机。
步骤一:使用转矩平衡法,利用开关磁阻电机在对齐位置与非对齐位置转矩为零的特点,通过给各相通入不同的脉冲电压,获得0°,7.5°,15°,22.5°位置电感特性,如图1所示;相电感计算公式如式(7)所示。
其中U为相电压,i为相电阻,R为绕组电阻。
步骤二:开关磁阻电机低速运行时,给非导通相注入脉冲电压,记录相绕组的电流波形和对应的位置信号。根据母线电压、开关管压降、二极管压降以及各相电流斜率计算非饱和电感,非饱和电感的计算由式(1)求解,结合所记录的相应位置信号,再根据电感的对称性得到开关磁阻电机不饱和电感特性。如图2所示。
步骤三:计算开关磁阻电机特殊角度,并将电感特性分为三个阶段,如图3所示。特殊角度的计算公式如式(2)所示。
步骤四:在0°<θ<θ1位置,相电感等于不饱和电感。
步骤五:在θ1<θ<θhr位置,建立相电感与不饱和电感的线性函数关系如式(3)所示。
将L(θhr,i)、L(θ1,i)、Lunhr,i)、Lun1,i)带入线性函数解算得到k(i)和L0(i),其中Lun1,i)和L(θ1,i)可通过不饱和电感线性区计算得到,L(θhr,i)近似为15°位置相电感。具体解算公式如式(8)所示。
步骤六:在θhr<θ<θa位置,建立相电感与不饱和电感的二次函数关系如式(4)所示。将L(θhr,i)、L(θa,i)、Lunhr)、Luna)和k(i)带入二次函数如式(9)所示,根据式(10)解算得到a(i)、b(i)、c(i)。其中L(θa,i)为22.5°位置相电感。得到不饱和电感与相电感二次函数关系后,根据不饱和电感的值使用二次函数求根公式反解相电感的值,从两根中取较大正值。具体解算公式如式(11)所示。
步骤七:利用上述步骤得到的完整电感特性,计算得到磁链特性与转矩特性曲线如图4、图5所示。
步骤八:基于步骤七得到的开关磁阻电机磁链和转矩特性,完成开关磁阻电机数学模型建立。
将使用本发明得到的开关磁阻电机与仿真和实验结果进行对比验证。图6所示为转速1560rpm,开关角5°,关断角15°电动状态下,开关磁阻电机相电流对比,实线为实验结果,点线为基于本发明计算得到特性仿真结果。图7所示为转速1000rpm,开关角30°,关断角40°发电状态下,开关磁阻电机相电流对比,实线为实验结果,点线为基于本发明计算得到特性仿真结果。以上对比验证证明本发明的准确性。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (4)

1.一种基于动静态测试的三相12/8极开关磁阻电机建模方法,其特征在于:使用转矩平衡法静态测试得到四条磁链曲线,在电机旋转时进行动态测试得到开关磁阻电机不饱和电感,基于以上测试得到的曲线进行拟合计算得到开关磁阻电机所有位置与电流下的磁链和转矩特性,完成建模,所述建模方法能够有效提升转矩平衡法的精确度,降低转子位置固定法的复杂度,该方法实现步骤如下:
步骤一:使用转矩平衡法,利用开关磁阻电机在对齐位置与非对齐位置转矩为零的特点,通过给各相通入不同的脉冲电压,获得0°、7.5°、15°和22.5°四个特殊位置电感特性,
步骤二:开关磁阻电机低速运行时,给非导通相注入脉冲电压,记录相绕组的电流波形,根据母线电压、开关管压降、二极管压降以及各相电流斜率计算非饱和电感,非饱和电感的计算由式(1)求解,结合所记录的相应位置信号,再根据电感的对称性得到开关磁阻电机不饱和电感特性,
式中Vdc、VD、VT、iph分别表示直流母线电压、二极管压降、开关管压降、相电流,
步骤三:通过式(2)计算开关磁阻电机特殊角度,并将电感特性分为三个阶段,其中θ1表示开关磁阻电机的定子前沿与转子后沿重合时的转子位置,θ2表示开关磁阻电机的定子前沿与转子前沿重合时的转子位置,θa表示开关磁阻电机的定子与转子中心对齐时的转子位置,θhr为转子前沿与定子中心线相重合时的转子位置,βr为转子极弧,βs为定子极弧,Nr为转子极对数,
步骤四:在0°<θ<θ1位置,认为相电感等于不饱和电感,
步骤五:在θ1<θ<θhr位置,认为不饱和电感与相电感为线性关系,可使用一次函数对不饱和电感与相电感进行拟合,如式(3),基于不饱和电感得到θ1到θhr位置相电感,
L(θ)=k(i)Lun(θ,i)+L0(i) (3)
步骤六:在θhr<θ<θa位置,认为不饱和电感与相电感为二次关系,可使用二次函数对不饱和电感与相电感进行拟合,如式(4),完成不饱和电感与相电感拟合后,基于不饱和电感数据与二次函数求根公式反解相电感数据,得到θhr到θ2位置相电感,
Lun(θ)=a(i)L(θ,i)2+b(i)L(θ,i)+c(i) (4)
步骤七:整合上述步骤得到[0,θa]的电感特性,根据开关磁阻电机电感特性的对称性可以得到完整的电感特性,由式(5)计算得到磁链特性,由式(6)计算得到转矩特性,
ψk(θ,ik)=Lk(θ,ik)ik (5)
其中,ik ,Tk,ψk分别表示第k相绕组电流、转矩和磁链,
步骤八:基于步骤七得到的开关磁阻电机的磁链和转矩特性,完成开关磁阻电机数学模型建立。
2.根据权利要求1所述一种基于动静态测试的三相12/8极开关磁阻电机建模方法,其特征在于:步骤一中通过转矩平衡法获取开关磁阻电机四个特殊位置的相电感。
3.根据权利要求1所述一种基于动静态测试的三相12/8极开关磁阻电机建模方法,其特征在于:步骤五中在θ1<θ<θhr位置,相电感与不饱和电感为线性关系,使用一次函数L(θ)=k(i)Lun(θ,i)+L0(i)对不饱和电感与相电感进行拟合,求解得到相电感。
4.根据权利要求1所述一种基于动静态测试的三相12/8极开关磁阻电机建模方法,其特征在于:步骤六中在θhr<θ<θa位置,不饱和电感与相电感为二次函数关系,使用二次函数Lun(θ)=a(i)L(θ,i)2+b(i)L(θ,i)+c(i)对不饱和电感与相电感进行拟合,利用二次函数求根公式反解相电感。
CN202110777173.4A 2021-07-09 2021-07-09 一种基于动静态测试的三相12/8极开关磁阻电机建模方法 Active CN113688499B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110777173.4A CN113688499B (zh) 2021-07-09 2021-07-09 一种基于动静态测试的三相12/8极开关磁阻电机建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110777173.4A CN113688499B (zh) 2021-07-09 2021-07-09 一种基于动静态测试的三相12/8极开关磁阻电机建模方法

Publications (2)

Publication Number Publication Date
CN113688499A CN113688499A (zh) 2021-11-23
CN113688499B true CN113688499B (zh) 2024-02-27

Family

ID=78576999

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110777173.4A Active CN113688499B (zh) 2021-07-09 2021-07-09 一种基于动静态测试的三相12/8极开关磁阻电机建模方法

Country Status (1)

Country Link
CN (1) CN113688499B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5107195A (en) * 1991-02-11 1992-04-21 General Electric Company Rotor position estimator for a switched reluctance machine using a lumped parameter flux/current model
CN102082544A (zh) * 2010-11-26 2011-06-01 江苏大学 无轴承同步磁阻电机转矩和悬浮力直接控制器及构造方法
CN104316876A (zh) * 2014-10-08 2015-01-28 西北工业大学 一种考虑互感耦合的获取三相12/8极srm磁链特性的快速方法
CN106017758A (zh) * 2016-05-19 2016-10-12 哈尔滨工程大学 一种电机动静态转矩同步在线测试装置与测试方法
CN108429504A (zh) * 2018-03-01 2018-08-21 西北工业大学 一种基于低成本位置传感器的开关磁阻电机转矩控制方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10637386B2 (en) * 2018-07-26 2020-04-28 Enedym Inc. Torque ripple reduction in switched reluctance machine

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5107195A (en) * 1991-02-11 1992-04-21 General Electric Company Rotor position estimator for a switched reluctance machine using a lumped parameter flux/current model
CN102082544A (zh) * 2010-11-26 2011-06-01 江苏大学 无轴承同步磁阻电机转矩和悬浮力直接控制器及构造方法
CN104316876A (zh) * 2014-10-08 2015-01-28 西北工业大学 一种考虑互感耦合的获取三相12/8极srm磁链特性的快速方法
CN106017758A (zh) * 2016-05-19 2016-10-12 哈尔滨工程大学 一种电机动静态转矩同步在线测试装置与测试方法
CN108429504A (zh) * 2018-03-01 2018-08-21 西北工业大学 一种基于低成本位置传感器的开关磁阻电机转矩控制方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于自适应模糊神经网络的开关磁阻电动机建模与仿真;丁文, 周会军, 鱼振民;中小型电机;20051010(05);全文 *
开关磁阻电机静态特性测试研究;曲兵妮;宋建成;郑丽君;郑建斌;;微电机;20090828(08);全文 *

Also Published As

Publication number Publication date
CN113688499A (zh) 2021-11-23

Similar Documents

Publication Publication Date Title
Rahman et al. Problems associated with the direct torque control of an interior permanent-magnet synchronous motor drive and their remedies
EP1535389B1 (en) Motor drive control
CN103618485B (zh) 一种无位置传感器无刷直流电机初始位置检测方法
CN103856139B (zh) 无速度传感器永磁同步电机转子磁极初始位置识别方法
CN110441643B (zh) 永磁同步电机控制***中逆变器功率管断路故障诊断方法
CN110212819B (zh) 一种用于高速无刷直流电机的换相误差补偿方法
CN107834935A (zh) 一种永磁同步电机转子的初始绝对位置检测方法
CN113300647B (zh) 一种永磁同步电机静止型交直轴电感辨识方法
CN110071674B (zh) 一种无位置传感器永磁同步电机最大转矩电流比控制方法
CN115001183B (zh) 磁悬浮开关磁阻电机及其悬浮力控制装置、方法和***
CN104316876A (zh) 一种考虑互感耦合的获取三相12/8极srm磁链特性的快速方法
CN110112965B (zh) 一种永磁同步电机反电动势观测方法
CN114337434A (zh) 一种考虑电感饱和效应的永磁电机参数离线辨识方法
CN113688499B (zh) 一种基于动静态测试的三相12/8极开关磁阻电机建模方法
Pengcheng et al. Offline parameter identification strategy of permanent magnet synchronous motor considering the inverter nonlinearities
CN108847802A (zh) 一种场定向无感启动的转子位置估算方法
TWI720657B (zh) 估測永磁同步馬達d-q軸電感的方法
CN114337450B (zh) 一种电流滞环环宽与电压自适应调节器的交流电机参数辨识方法
Morimoto et al. Realization of high torque density encoderless servo drive system
Szalay et al. Saliency model extension for sensorless initial position and polarity detection of permanent magnet synchronous motors
CN110112976B (zh) 一种利用牛顿拉普逊法的永磁同步电机参数辨识方法
CN113179066A (zh) 一种无传感器永磁同步电机逆变器死区的时间补偿方法
CN114696676A (zh) 基于脉冲电压法的离线参数辨识方法
CN113300655B (zh) 一种基于误差补偿的开关磁阻电机无位置传感器控制方法
Zhang et al. High Performance Sensorless Control of Brushless DC Motor Based on A Novel Sliding Mode Observer

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