CN113886987B - 一种基于改进绕组函数法的笼型转子感应电动机建模方法 - Google Patents
一种基于改进绕组函数法的笼型转子感应电动机建模方法 Download PDFInfo
- Publication number
- CN113886987B CN113886987B CN202111158461.8A CN202111158461A CN113886987B CN 113886987 B CN113886987 B CN 113886987B CN 202111158461 A CN202111158461 A CN 202111158461A CN 113886987 B CN113886987 B CN 113886987B
- Authority
- CN
- China
- Prior art keywords
- rotor
- winding
- function
- motor
- stator
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/60—Other road transportation technologies with climate change mitigation effect
- Y02T10/64—Electric machine technologies in electromobility
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Windings For Motors And Generators (AREA)
Abstract
本专利提出了一种基于改进绕组函数法的笼型转子感应电动机建模方法,它涉及电机动态分析领域,其特征在于适用于笼型转子感应电动机,基于改进绕组函数法和转子多回路法,所使用的转子多回路法做了改进,通过对称性进行简化减少了计算量。首先,建立适用于绕组函数法笼型转子绕组转子多回路模型,并充分考虑电机定转子开槽造成的影响,利用改进绕组函数法求取电机的定子绕组和转子绕组的电感参数,并将求得的电感参数带入感应电动机状态方程中进行求解。本发明的优势是仅需要知道电机的几何尺寸和绕组参数即可对电机进行动态分析,适用于电机设计初期及需要做大量参数调整时的分析,同时编程简单,计算速度快。
Description
技术领域
本发明涉及电机动态分析领域,对传统的感应电机转子多回路法做了改进,使其更适用于绕组函数法的电感求取,同时借助感应电动机的状态方程确定电机的性能。
背景技术
感应电动机因为结构简单,成本低廉,运行可靠被广泛应用在工业领域中,目前基于等效电路和等效磁路结合的感应电动机稳态性能的计算已经十分成熟,但在故障分析,启动性能等方面仍然需要感应电动机的动态分析。电感参数是沟通电与磁,反应机电能量转换的重要参数,在感应电动机的动态分析中具有十分重要的地位。
目前,电感参数计算的方法包括解析法和有限元法,尽管有限元的计算更加准确,但在研究某些问题时需要重复修改物理模型,仿真往往需要耗费大量的时间。解析法则更加简单,易于程序化,求解速度快,满足一定准确性的条件下,更加适合研究这类问题。
绕组函数法是基于安培环路定律,以磁链为中间变量的一种电感计算方法,其计算简单,时间快,准确性能够满足工程上的需要。对于感应电动机来说,怎样建立适用于绕组函数法的笼型绕组电路模型是研究这类问题的关键,同时为了电感计算的准确性,怎样充分考虑齿槽效应和饱和效应也是该领域的难点之一。
发明内容
针对上述技术背景,本文提出了一种基于改进绕组函数法的笼型转子感应电动机建模方法,其特征在于:所述方法基于改进绕组函数法和转子多回路法,所使用的转子多回路法做了改进,通过对称性进行简化减少了计算量。本方法适用于笼型转子感应电动机,首先建立了笼型转子的简化转子多回路模型,然后建立了考虑了定转子开槽影响的改进绕组函数求取定转子各相绕组电感,将所求电感带入到状态方程中求解,得到电机性能参数。
本发明的具体实现过程如下:
步骤一:
建立笼型转子多回路模型:首先利用感应电动机的笼型转子的对称性,按照电机的极数2p将转子划分为p对极,每对极占据360°电角度,其中相差180°电角度的两根导条为一匝线圈的两个线圈边,每两对极下的相差360°电角度的线圈互相串联,将拥有Q2根导条的笼型绕组划分为Q2/2p相,每相有p个线圈的绕组;
步骤二:
根据电机的几何尺寸,建立考虑定转子开槽影响的气隙函数:首先选择参考坐标系,将电机气隙划分为定子气隙gs、转子气隙gr和均匀气隙长度g0,假设电机的定子槽口高为h11,槽口宽为b11,转子槽口高为h21,槽口宽为b21;在定子槽口对应的区域,设定子气隙长度gs=h11,在转子槽口对应的区域,设转子气隙长度gr=h21,其他区域,gs和gr均为0,最终考虑定转子开槽影响的总气隙长度为:
g=gs+gr+g0 (1)
步骤三:
根据电机的绕组参数,建立定转子绕组的匝数函数:定子绕组可以分为m相,对于单层绕组,每相绕组可以分为Q1/m个线圈,双层绕组每相绕组可以分为2Q1/m个线圈,每个线圈由n1匝构成;转子绕组由Q2/2p相绕组构成,每相绕组由p个线圈构成,每个线圈由两根导条构成,其匝数为一匝;首先选定参考方向,假设线圈中电流的流动方向,每个线圈边中的电流与参考方向相同则匝数为正,与参考方向相反匝数为负,每个线圈边的匝数函数为阶跃函数,其幅值为该线圈边的匝数,将每相各个线圈边的匝数函数进行线性叠加即为该相的匝数函数,依次求出定转子绕组各相的匝数函数n。
步骤四:
根据电机绕组的各相匝数函数和气隙函数,建立改进绕组函数模型:电机的改进组函数与气隙函数和匝数函数有关,能够充分考虑到气隙的不均匀,改进绕组函数为
M(φ,θ)=n(φ,θ)-<M(φ,θ)〉 (2)
其中:M(φ,θ)为改进绕组函数;n(φ,θ)为匝数函数;<g-1(φ,θ)>为气隙逆函数的平均值。
步骤五:
根据电机绕组的绕组函数,求取电机各相绕组的自感与互感:电感由绕组函数、匝数函数和气隙函数计算得到,其计算公式为:
其中μ0为真空磁导率,r为转子外径,l为电机轴向长度,依次求取定转子各相绕组自感和各相绕组之间的互感。
步骤六
根据电机的电感参数,建立感应电动机的电压和转矩方程:首先根据步骤五求出来的电感参数,建立电感矩阵,根据电感矩阵和电阻矩阵,建立感应电动机的电压平衡和转矩方程;感应电动机的电压和转矩方程为:
其中u为电压矩阵,i为电流矩阵,R为电阻矩阵,L为电感矩阵,ψ为磁链矩阵,Te为电磁转矩,TL为负载转矩,RΩ为旋转阻力系数,J为转动惯量,Ω为机械角速度,θ为转子角位移;
步骤七:
根据电压和转矩方程,建立以磁链为状态变量的感应电动机状态方程:
步骤八:
用四阶龙格—库塔方法求解状态方程,得到感应电动机的电磁转矩Te,定子电流i和转子转速Ω。
与现有技术相比,本发明的有益效果是:
1、本发明所提出的方法仅需要电机的几何尺寸和绕组参数,即可进行电机的动态分析,使用简单。
2、本发明所提出的方法适用于电机设计初期及需要做大量参数调整时的动态分析,计算速度快。
3、本发明适合计算机辅助设计,易于编程。
附图说明
图1为一种基于改进绕组函数法的笼型转子感应电动机建模方法流程图。
图2为考虑定转子开槽影响的气隙函数。
图3为定子A相绕组的匝数函数。
图4为定子C相绕组的绕组函数。
图5为定子A相绕组的自感与位置角之间的关系。
图6为定子A相绕组与转子第一相绕组之间的互感。
图7为电机的定子A相电流随时间变化的曲线。
图8为电机的电磁转矩随时间变化的曲线。
图9为电机的转子转速随时间变化的曲线。
具体实施方式
下面以一台三相感应电动机为对象,结合具体实例来详细描述本发明专利所叙述的方法,该方法的流程图如图1所示。
表1电机基本参数
定子内径 | 136.0mm | 定子槽数 | 36 |
转子外径 | 135.2mm | 转子槽数 | 32 |
气隙长度 | 0.4mm | 定子槽口宽 | 3.5mm |
轴向长度 | 160mm | 转子槽口宽 | 1mm |
极数 | 4 | 定子槽口高 | 0.8mm |
每槽导体数 | 35 | 转子槽口高 | 0.5mm |
步骤一:
建立笼型转子多回路模型:首先利用感应电动机的笼型转子的对称性,将4极电机划分为2对极,每对极占据360°电角度,其中相差180°电角度的两根导条为一匝线圈的两个线圈边,每两对极下的相差360°的线圈互相串联,将拥有32根导条的笼型绕组划分为有8相,每相绕组有2个线圈的转子等效绕组。
步骤二:
根据电机的几何尺寸,建立考虑开槽影响的气隙函数:首先选择定子参考坐标系,将电机气隙划分为定子气隙gs、转子气隙gr和均匀气隙长度g0。
根据公式(1),总气隙长度为:
总气隙长度如图2所示,本例中的气隙函数与转子位置有关,每一圆周位置都对应一个气隙值。
步骤三:
根据电机的绕组参数,建立定转子绕组的匝数函数:定子绕组可以分为3相,每相绕组可以分为12个线圈,每个线圈由35匝构成;转子绕组由8相绕组构成,每相绕组由2个线圈构成,每个线圈由两根导条构成,其匝数为一匝;首先选定参考方向,假设线圈中电流的流动方向,每个线圈边中的电流与参考方向相同则匝数为正,与参考方向相反匝数为负,每个线圈边的匝数函数为阶跃函数,其幅值为该线圈边的匝数,将每相各个线圈边的匝数函数进行线性叠加即为该相的匝数函数,依次求出定转子各相的匝数函数A相绕组的匝数函数如图3所示,B相和C相匝数函数分别滞后A相120°和240°电角度。
步骤四:
根据电机绕组的各相匝数函数和气隙函数,建立改进绕组函数模型:根据公式(2)和公式(3),求取各相绕组的改进绕组函数,其中C相绕组的改进绕组函数如图4所示,B相和A相的改进绕组函数与C相分别相差120°和240°电角度。
步骤五:
根据电机绕组的绕组函数,求取电机各相绕组的自感与互感:根据公式(4)计算各相绕组的自感和互感。
以A相和B相绕组为例,A相自感与A相和B相之间的互感计算公式如公式(8)、(9)所示:
其中μ0为真空磁导率,r为转子外径,l为电机轴向长度,可按照电机参数表查取对应值;依次求取定转子各相绕组自感和各相绕组之间的互感。A相自感与A相和B相之间的互感如图5和图6所示,绕组自感在一固定值波动,这是由于定转子开槽的影响,而定转子互感则随转子位置的改变而改变。
步骤六:
根据电机的电感,建立感应电动机的电压和转矩方程:首先根据步骤五求出来的电感系数,建立电感矩阵,根据电感矩阵和电阻矩阵,依据公式(5)建立感应电动机的电压平衡和转矩方程,其中
其中us为定子电压矩阵;ur为转子电压矩阵;is为定子电流矩阵;ir为转子电流矩阵;Rs为定子电阻矩阵;Rr为转子电阻矩阵;Lss为定子自感矩阵;Lrr为转子自感矩阵;Lsr和Lrs为定转子互感矩阵;ψs和ψr分别为定转子磁链矩阵。
步骤七:
根据电压和转矩方程,按照公式(6)建立以磁链为状态变量的感应电动机状态方程:
步骤八:
用四阶龙格—库塔方法求解状态方程,得到感应电动机的定子A相电流iA,电磁转矩Te和转子转速Ω,分别如图7、图8和图9所示。该结果是在电机空载状态下的得出的,其稳定后定子电流幅值为31.56A,转速为1498r/min,接近空载转速1500r/min,电磁转矩为3.43N·m,接近0。
以上实例描述了本发明专利的基本逻辑与运行思路,但本行业的技术人员应该了解,本发明专利不受上述实例的限制,上述实例和说明书中描述的只是说明本专利的逻辑与思路,在不脱离本发明专利精神和范围的前提下,本发明专利还会有各种变化和改进,这些变化和改进都将落入要求保护的本发明专利范围内。
Claims (1)
1.一种基于改进绕组函数法的笼型转子感应电动机建模方法,其特征在于:所述方法基于改进绕组函数法和转子多回路法,本方法适用于笼型转子感应电动机,所述方法的具体实现过程为:
步骤一:
建立笼型转子多回路模型:首先利用感应电动机的笼型转子的对称性,按照电机的极数2p将转子划分为p对极,每对极占据360°电角度,其中相差180°电角度的两根导条为一匝线圈的两个线圈边,每两对极下的相差360°电角度的线圈互相串联,将拥有Q2根导条的笼型绕组划分为Q2/2p相,每相有p个线圈的绕组;
步骤二:
根据电机的几何尺寸,建立考虑定转子开槽影响的气隙函数:首先选择参考坐标系,将电机气隙划分为定子气隙gs、转子气隙gr和均匀气隙长度g0,假设电机的定子槽口高为h11,槽口宽为b11,转子槽口高为h21,槽口宽为b21;在定子槽口对应的区域,设定子气隙长度gs=h11,在转子槽口对应的区域,设转子气隙长度gr=h21,其他区域,gs和gr均为0,最终考虑定转子开槽影响的总气隙长度为:
g=gs+gr+g0
步骤三:
根据电机的绕组参数,建立定转子绕组的匝数函数:定子绕组可以分为m相,对于单层绕组,每相绕组可以分为Q1/m个线圈,双层绕组每相绕组可以分为2Q1/m个线圈,每个线圈由n1匝构成;转子绕组由Q2/2p相绕组构成,每相绕组由p个线圈构成,每个线圈由两根导条构成,其匝数为一匝;首先选定参考方向,假设线圈中电流的流动方向,每个线圈边中的电流与参考方向相同则匝数为正,与参考方向相反匝数为负,每个线圈边的匝数函数为阶跃函数,其幅值为该线圈边的匝数,将每相各个线圈边的匝数函数进行线性叠加即为该相的匝数函数,依次求出定转子绕组各相的匝数函数n;
步骤四:
根据电机绕组的各相匝数函数和气隙函数,建立改进绕组函数模型:电机的改进绕组函数与气隙函数和匝数函数有关,能够充分考虑到气隙的不均匀,改进绕组函数为
M(φ,θ)=n(φ,θ)-<M(φ,θ)>
其中:M(φ,θ)为改进绕组函数;n(φ,θ)为匝数函数;<g-1(φ,θ)>为气隙逆函数的平均值;
步骤五:
根据电机绕组的绕组函数,求取电机各相绕组的自感与互感:电感由绕组函数、匝数函数和气隙函数计算得到,其计算公式为:
其中μ0为真空磁导率,r为转子外径,l为电机轴向长度,依次求取定转子各相绕组自感和各相绕组之间的互感;
步骤六:
根据电机的电感参数,建立感应电动机的电压和转矩方程:首先根据步骤五求出来的电感参数,建立电感矩阵,根据电感矩阵和电阻矩阵,建立感应电动机的电压平衡和转矩方程;感应电动机的电压和转矩方程为:
其中u为电压矩阵,i为电流矩阵,R为电阻矩阵,L为电感矩阵,ψ为磁链矩阵,Te为电磁转矩,TL为负载转矩,RΩ为旋转阻力系数,J为转动惯量,Ω为机械角速度,θ为转子角位移;
步骤七:
根据电压和转矩方程,建立以磁链为状态变量的感应电动机状态方程:
步骤八:
用四阶龙格—库塔方法求解状态方程,得到感应电动机的电磁转矩Te,定子电流i和转子转速Ω。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111158461.8A CN113886987B (zh) | 2021-09-29 | 2021-09-29 | 一种基于改进绕组函数法的笼型转子感应电动机建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111158461.8A CN113886987B (zh) | 2021-09-29 | 2021-09-29 | 一种基于改进绕组函数法的笼型转子感应电动机建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113886987A CN113886987A (zh) | 2022-01-04 |
CN113886987B true CN113886987B (zh) | 2023-06-30 |
Family
ID=79004581
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111158461.8A Active CN113886987B (zh) | 2021-09-29 | 2021-09-29 | 一种基于改进绕组函数法的笼型转子感应电动机建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113886987B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2687868A1 (es) * | 2018-03-08 | 2018-10-29 | Universidade De Santiago De Compostela | Procedimiento y producto de programa informático para acelerar el cálculo del estado estacionario de un motor de inducción de jaula de ardilla |
CN110737990A (zh) * | 2019-10-21 | 2020-01-31 | 河南理工大学 | 直驱永磁同步电机定子绕组故障状态数学模型建立方法 |
CN110765581A (zh) * | 2019-09-23 | 2020-02-07 | 武汉船用电力推进装置研究所(中国船舶重工集团公司第七一二研究所) | 一种十二相永磁同步电机的建模方法 |
KR20200078861A (ko) * | 2018-12-24 | 2020-07-02 | 한국산업기술대학교산학협력단 | 유도 전동기의 회로 파라미터 추출 장치 및 그 방법 |
CN113094911A (zh) * | 2021-04-16 | 2021-07-09 | 江苏大学 | 一种磁场调制永磁容错电机高功率因数设计方法 |
CN113452209A (zh) * | 2021-03-08 | 2021-09-28 | 中国人民解放军海军工程大学 | 一种多相笼型感应电机磁场计算方法及系统 |
-
2021
- 2021-09-29 CN CN202111158461.8A patent/CN113886987B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2687868A1 (es) * | 2018-03-08 | 2018-10-29 | Universidade De Santiago De Compostela | Procedimiento y producto de programa informático para acelerar el cálculo del estado estacionario de un motor de inducción de jaula de ardilla |
KR20200078861A (ko) * | 2018-12-24 | 2020-07-02 | 한국산업기술대학교산학협력단 | 유도 전동기의 회로 파라미터 추출 장치 및 그 방법 |
CN110765581A (zh) * | 2019-09-23 | 2020-02-07 | 武汉船用电力推进装置研究所(中国船舶重工集团公司第七一二研究所) | 一种十二相永磁同步电机的建模方法 |
CN110737990A (zh) * | 2019-10-21 | 2020-01-31 | 河南理工大学 | 直驱永磁同步电机定子绕组故障状态数学模型建立方法 |
CN113452209A (zh) * | 2021-03-08 | 2021-09-28 | 中国人民解放军海军工程大学 | 一种多相笼型感应电机磁场计算方法及系统 |
CN113094911A (zh) * | 2021-04-16 | 2021-07-09 | 江苏大学 | 一种磁场调制永磁容错电机高功率因数设计方法 |
Non-Patent Citations (2)
Title |
---|
张经纬 ; 祝后权 ; 黄振华 ; 杨高 ; 戴勇 ; .多相整距绕组感应电机的机电动力系统模型研究.舰船科学技术.2010,(第08期),全文. * |
鲍晓华 ; 程志恒 ; 王汉丰 ; 狄冲 ; 朱庆龙 ; .基于改进绕组函数法的感应电机偏心电感参数计算.电工技术学报.2016,(第06期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN113886987A (zh) | 2022-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10848025B2 (en) | Permanent magnet brushless motor having high winding factor, and design and fault-tolerant control methods thereof | |
Zarko et al. | Analytical calculation of magnetic field distribution in the slotted air gap of a surface permanent-magnet motor using complex relative air-gap permeance | |
WO2022110274A1 (zh) | 一种磁场调制永磁电机损耗分析与抑制方法 | |
CN110851949B (zh) | 分析多层磁障永磁辅助式同步磁阻电机电磁性能的方法 | |
CN113094911B (zh) | 一种磁场调制永磁容错电机高功率因数设计方法 | |
More et al. | Power density improvement of three phase flux reversal machine with distributed winding | |
Zhao et al. | Performance analysis of a coreless axial-flux PMSM by an improved magnetic equivalent circuit model | |
CN112737160B (zh) | 一种集中绕组外转子游标电机提升功率因数的方法 | |
CN104917348A (zh) | 电动车用大功率奇数分数槽电机 | |
Mallampalli et al. | Six-phase pole-changing winding induction machines with improved performance | |
CN108206660A (zh) | 一种新型五相表嵌式永磁同步电机的转矩解析方法 | |
Shen et al. | Design optimization and performance investigation of linear doubly salient slot permanent magnet machines | |
Awah et al. | Impact of structural dimensions and poles on the torque performance of dual-stator permanent magnet machines | |
CN113886987B (zh) | 一种基于改进绕组函数法的笼型转子感应电动机建模方法 | |
CN110765581B (zh) | 一种十二相永磁同步电机的建模方法 | |
Røkke et al. | Analytical calculation of yoke flux patterns in fractional-slot permanent magnet machines | |
Shima et al. | Steady-state magnetic circuit analysis of salient-pole synchronous machines considering cross-magnetization | |
Siregar et al. | A New Technique to Reduce the Cogging Torque of Integral Slot Number in Permanent Magnet Synchronous Machine | |
CN108984844B (zh) | 一种无刷双馈电机性能的分析方法及装置 | |
CN107729628B (zh) | 一种三相电励磁双凸极电机非线性电感建模方法 | |
CN115395854A (zh) | 一种基于主动式无位置传感器策略的永磁容错轮毂电机及其驱动、设计方法 | |
CN110401273B (zh) | 低谐波的分数槽集中绕组设计方法 | |
Cao et al. | Cogging torque reduction for outer rotor interior permanent magnet synchronous motor | |
Jeong et al. | Performance comparison of conventional and segmental rotor type switched reluctance motor | |
Mirnikjoo et al. | Winding function model of a 6/7 variable flux reluctance machine |
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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Xia Yunyan Inventor after: Wang Yanbo Inventor after: Xu Yongming Inventor after: Ai Mengmeng Inventor before: Wang Yanbo Inventor before: Xia Yunyan Inventor before: Xu Yongming Inventor before: Ai Mengmeng |
|
GR01 | Patent grant | ||
GR01 | Patent grant |