CN110162892B - 一种大型同步调相机定子端部漏感计算方法 - Google Patents

一种大型同步调相机定子端部漏感计算方法 Download PDF

Info

Publication number
CN110162892B
CN110162892B CN201910441097.2A CN201910441097A CN110162892B CN 110162892 B CN110162892 B CN 110162892B CN 201910441097 A CN201910441097 A CN 201910441097A CN 110162892 B CN110162892 B CN 110162892B
Authority
CN
China
Prior art keywords
stator
finite element
element model
magnetic
solving
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
CN201910441097.2A
Other languages
English (en)
Other versions
CN110162892A (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.)
Chongqing University
State Grid Hubei Electric Power Co Ltd
Original Assignee
Chongqing University
State Grid Hubei Electric Power Co Ltd
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 Chongqing University, State Grid Hubei Electric Power Co Ltd filed Critical Chongqing University
Priority to CN201910441097.2A priority Critical patent/CN110162892B/zh
Publication of CN110162892A publication Critical patent/CN110162892A/zh
Application granted granted Critical
Publication of CN110162892B publication Critical patent/CN110162892B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

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 Ac Motors In General (AREA)

Abstract

本发明涉及一种大型同步调相机定子端部漏感计算方法,属于同步电机参数计算领域。该方法包括:S1:同步电机电磁有限元建模:分别建立2D有效段有限元模型和3D端部有限元模型;S2:求解饱和状态下的同步电机磁密及磁导率分布;S3:锁定对应饱和状态的磁导率,线性求解磁场储能;S4:基于磁场储能计算2D/3D模型电感,并对3D端部模型电感扣除有效段含量,从而分离出端部漏感。本发明能够考虑实际运行中磁路饱和对端部磁场分布的影响,得到更加准确的大型同步调相机定子端部漏感参数,从而为更好地预判电机运行性能、分析电机动态行为提供。

Description

一种大型同步调相机定子端部漏感计算方法
技术领域
本发明属于同步电机参数计算领域,涉及一种大型同步调相机定子端部漏感的计算方法。
背景技术
相对于SVC/SVG等无功补偿装置,大容量调相机具备更高的可靠性和更强暂态能力,能有效改善直流换相失败、抑制系统暂态过电压,因此在直流输电系统中得到重要应用。调相机的暂态无功能力和对系统电压稳定的支持能力直接取决于其暂态参数,而定子绕组的端部漏感是超瞬变电感的重要组成部分。能否准确地计算端部漏感,是准确评估调相机暂态无功能力和预测系统电压安全的关键。然而,大型同步调相机定子绕组端部结构及电磁关系十分复杂,绕组型式及布置、端部结构件、定转子边段铁心等都会对定子绕组端部漏磁产生影响,准确计算电机端部漏感较为困难。目前在电机设计中,一般采用传统经验公式计算端部漏电感,近年来更多地采用三维数值分析方法,然而这些方法大都未计及端部磁路饱和的影响,认为端部漏感是一个定值,实际上端部结构的磁饱和会影响端部漏电感的大小,忽略饱和影响必然会造成端部漏感计算的不准确。因此,为了获取准确的调相机定子端部漏感参数,以评估调相机暂态无功能力、预测系统电压安全性能,亟需发明一种考虑饱和效应影响的大型调相机定子端部漏感计算方法。
发明内容
有鉴于此,本发明的目的在于提供一种大型同步调相机定子端部漏感的计算方法,通过考虑实际运行中磁路饱和对端部磁场分布的影响,获取更加准确的同步调相机定子端部漏感参数,为更好地分析调相机动态行为、评估调相机暂态无功能力、预测系统电压安全性能提供支撑。
为达到上述目的,本发明提供如下技术方案:
一种大型同步调相机定子端部漏感计算方法,具体包括以下步骤:
S1:同步电机电磁有限元建模:分别建立2D有效段有限元模型和3D端部有限元模型;
S2:分别基于步骤S1建立的有限元模型,进行所需工况的非线性时域电磁场仿真,此时模型中应设置合理的定、转子电流及其相位关系;求解饱和状态下的同步电机磁密及磁导率分布,并将其锁定导出;
S3:将导出的磁导率分布再导入至新的求解模型,仅在定子侧加载三相对称电流,进行冻结磁导率状态下的线性电磁场仿真,线性求解定子磁链或电机储能;
S4:基于磁链或储能法分别求解2D有效段有限元模型和3D端部有限元模型的电感,并将3D端部有限元模型的电感扣除包含2D有效段有限元模型的部分,从而分离出定子端部饱和漏感。
进一步,所述步骤S1中,2D有效段有限元模型包括有2D求解域,3D端部有限元模型包括有3D求解域;
2D求解域对应定子有效段长度lactive=ltotal-2lend,其中ltotal为电机定子铁心总长,lend为两侧定子端部铁心轴向长度;
3D求解域包括能够考虑端部效应的足够长的边段铁心、压指、压板等结构件,以及能反映其空间分布的定子端部绕组。
进一步,所述步骤S3中具体包括:
S31:将导出的磁导率分布再导入至新的求解模型,激励加载方式为:
Ia=Im cos(ωt+θ)
Figure BDA0002072041810000021
Figure BDA0002072041810000022
If=0
其中,Ia、Ib、Ic是定子三相电流,Im为电流峰值,ω为电角频率,θ为定子A相初始相位角,If为励磁电流;
S32:进行冻结磁导率状态下的线性电磁场仿真,求解定子某相最大磁链或磁场储能;
磁链计算方法:
Figure BDA0002072041810000023
磁场储能计算方法:
Figure BDA0002072041810000024
其中,J是一匝内的电流密度,A是矢量位,v是区域的体积,n是匝数,B是磁通密度,H是磁场强度。
进一步,所述步骤S4中具体包括:分离出端部漏感的方法:
Lσ_end=Ltotal-Lactivelend/lactive
其中,Lσ_end是单边端部漏感,Ltotal是单边端部总电感,Lactive是2D有效段有限元模型计算出的铁心有效段对应的电感,lend是3D端部有限元模型的铁心长度。
本发明的有益效果在于:本发明能够准确考虑实际运行中磁路饱和对端部磁场分布的影响,获取更加准确的同步调相机定子端部漏感参数,为更好地分析调相机动态行为、评估调相机暂态无功能力、预测系统电压安全性能提供支撑。
本发明的其他优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其他优点可以通过下面的说明书来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作优选的详细描述,其中:
图1为本发明所述方法的流程图;
图2为本发明的同步电机定子端部有限元模型的2D/3D求解域;
图3为2D有效段有限元模型图;
图4为3D端部有限元模型图;
图5为2D模型额定工况饱和磁导率分布图;
图6为3D模型额定工况饱和磁导率分布图;
图7为本发明所述方法计算出的端部漏感标幺值结果。
具体实施方式
以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需要说明的是,以下实施例中所提供的图示仅以示意方式说明本发明的基本构想,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。
如图1所示,本发明所述的一种大型同步调相机定子端部漏感计算方法,具体包括以下步骤:
S1:同步电机电磁有限元建模:分别建立2D有效段有限元模型和3D端部有限元模型;
S2:分别基于步骤S1建立的有限元模型,进行所需工况的非线性时域电磁场仿真,此时模型中应设置合理的定、转子电流及其相位关系;求解饱和状态下的同步电机磁密及磁导率分布,并将其锁定导出;
S3:将导出的磁导率分布再导入至新的求解模型,仅在定子侧加载三相对称电流,进行冻结磁导率状态下的线性电磁场仿真,线性求解定子磁链或电机储能;
S4:基于磁链或储能法分别求解2D有效段有限元模型和3D端部有限元模型的电感,并将3D端部有限元模型的电感扣除包含2D有效段有限元模型的部分,从而分离出定子端部饱和漏感。
图2为本发明的同步电机定子端部有限模型的2D/3D求解域,具体为:
电机定子铁心总长为ltotal,所求解的端部铁心轴向长度为lend,则2D模型的有效段长度应为:
lactive=ltotal-2lend
3D有限元求解域包括能够考虑端部效应的足够长的边段铁心及压指、压板等结构件,以及能反映其空间分布的定子端部绕组。
图3为2D有效段有限元模型图。图4为3D端部有限元模型图。图5为2D额定工况饱和磁导率分布图。图6为3D额定工况模型饱和磁导率分布图,左侧为定子铁心,右侧为转子铁心。
图7为本发明的改进方法计算出的端部漏感标幺值结果。图7中对比了传统方法计算出的不饱和值和本发明方法计算出的三个对比工况的饱和值。可以看出,本发明结果能够反映不同工况饱和程度对端部漏感的影响,且d轴和q轴饱和程度有所不同。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (1)

1.一种大型同步调相机定子端部漏感计算方法,其特征在于,该方法具体包括以下步骤:
S1:同步电机电磁有限元建模:分别建立2D有效段有限元模型和3D端部有限元模型;
S2:分别基于步骤S1建立的有限元模型,进行所需工况的非线性时域电磁场仿真,设置合理的定、转子电流及其相位关系;求解饱和状态下的同步电机磁密及磁导率分布,并将其锁定导出;
S3:将导出的磁导率分布再导入至新的求解模型,仅在定子侧加载三相对称电流,进行冻结磁导率状态下的线性电磁场仿真,线性求解定子磁链或电机储能;
S4:基于磁链或储能法分别求解2D有效段有限元模型和3D端部有限元模型的电感,并将3D端部有限元模型的电感扣除包含2D有效段有限元模型的部分,从而分离出定子端部饱和漏感;
所述步骤S1中,2D有效段有限元模型包括有2D求解域,3D端部有限元模型包括有3D求解域;
2D求解域对应定子有效段长度lactive=ltotal-2lend,其中ltotal为电机定子铁心总长,lend为两侧定子端部铁心轴向长度;
3D求解域包括能够考虑端部效应的足够长的边段铁心、压指、压板,以及能反映其空间分布的定子端部绕组;
所述步骤S3中具体包括:
S31:将导出的磁导率分布再导入至新的求解模型,激励加载方式为:
Ia=Im cos(ωt+θ)
Figure FDA0004054502600000011
Figure FDA0004054502600000012
If=0
其中,Ia、Ib、Ic是定子三相电流,Im为电流峰值,ω为电角频率,θ为定子A相初始相位角,If为励磁电流;
S32:进行冻结磁导率状态下的线性电磁场仿真,求解定子某相最大磁链或磁场储能;
磁链计算方法:
Figure FDA0004054502600000013
磁场储能计算方法:
Figure FDA0004054502600000021
其中,J是一匝内的电流密度,A是矢量磁位,v是区域的体积,n是匝数,B是磁通密度,H是磁场强度;
所述步骤S4中具体包括:分离出端部漏感的方法:
Lσ_end=Ltotal-Lactivelend/lactive
其中,Lσ_end是单边端部漏感,Ltotal是单边端部总电感,Lactive是2D有效段有限元模型计算出的铁心有效段对应的电感,lend是3D端部有限元模型的铁心长度,lactive为2D求解域对应定子有效段长度。
CN201910441097.2A 2019-05-24 2019-05-24 一种大型同步调相机定子端部漏感计算方法 Active CN110162892B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910441097.2A CN110162892B (zh) 2019-05-24 2019-05-24 一种大型同步调相机定子端部漏感计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910441097.2A CN110162892B (zh) 2019-05-24 2019-05-24 一种大型同步调相机定子端部漏感计算方法

Publications (2)

Publication Number Publication Date
CN110162892A CN110162892A (zh) 2019-08-23
CN110162892B true CN110162892B (zh) 2023-04-07

Family

ID=67632741

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910441097.2A Active CN110162892B (zh) 2019-05-24 2019-05-24 一种大型同步调相机定子端部漏感计算方法

Country Status (1)

Country Link
CN (1) CN110162892B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113239601B (zh) * 2021-07-13 2021-12-10 国网江西省电力有限公司电力科学研究院 一种调相机暂稳态参数实时提取方法及装置

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108563912A (zh) * 2018-05-11 2018-09-21 北京理工大学 一种表贴式永磁同步电机气隙磁场的解析方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CH244943A (de) * 1945-10-04 1946-10-15 Oerlikon Maschf Mechanischer Stromrichter.
CN107368653B (zh) * 2017-07-20 2020-07-28 云南电网有限责任公司电力科学研究院 一种变压器漏磁场的二维有限元仿真方法和系统
CN107317331A (zh) * 2017-07-25 2017-11-03 华北电力大学 一种特高压变压器轴对称直流偏磁仿真模型

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108563912A (zh) * 2018-05-11 2018-09-21 北京理工大学 一种表贴式永磁同步电机气隙磁场的解析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
林楠 ; 王东 ; 魏锟 ; 张庆湖 ; 易新强 ; .新型混合励磁同步电机的数学模型与等效分析.电工技术学报.(第03期),全文. *

Also Published As

Publication number Publication date
CN110162892A (zh) 2019-08-23

Similar Documents

Publication Publication Date Title
Liu et al. Nonlinear equivalent magnetic network of a linear permanent magnet vernier machine with end effect consideration
Liang et al. Analytical model and design of spoke-type permanent-magnet machines accounting for saturation and nonlinearity of magnetic bridges
Rahim et al. Performance analysis of salient-pole self-excited reluctance generators using a simplified model
Wang et al. A phase-domain synchronous machine model with constant equivalent conductance matrix for EMTP-type solution
Wu et al. Axial‐flux permanent‐magnet synchronous generator with coreless armature and non‐integral coil–pole ratio
Pechlivanidou et al. A comparative study on 2D and 3D magnetic field analysis of permanent magnet synchronous motor using FEM simulations
CN110162892B (zh) 一种大型同步调相机定子端部漏感计算方法
Liang et al. Numerical calculation of circulating current losses in stator transposition bar of large hydro‐generator
Zhang et al. Improved calculation method for inductance value of the air-gap inductor
Rasilo et al. Experimental determination and numerical evaluation of core losses in a 150‐kVA wound‐field synchronous machine
Zhou et al. Electromagnetic design of high-speed and high-thrust cross-shaped linear induction motor
von Pfingsten et al. Highly efficient approach to the simulation of variable‐speed induction motor drives
Jun et al. Voltage regulation rate and THD optimization analysis of coreless axial flux PM synchronous generator for wind power generation
Li et al. Parametric study for the design of the end region of large synchronous generators based on three-dimensional transient finite element analysis
Taran et al. A hybrid analytical and fe-based method for calculating ac eddy current winding losses taking 3d effects into account
Yu et al. Analytical model for magnetic field calculation of SPMSM with chamfered pole considering iron core saturation
Kazemi Sangdehi et al. Analysis and optimization of a 12/14 double‐stator flux‐switching machine using low cost magnet
Lee et al. Influence of electrical steel characteristics on efficiency of industrial induction motors
Xuan et al. Modeling for the design of fractional slot PM machines with concentrated windings protected from demagnetization during three-phase short circuit
Gholami et al. An extended winding function model for induction machine modelling considering saturation effect
Nie et al. Compensation analysis of longitudinal end effect in three‐phase switched reluctance linear machines
Chen et al. Hybrid analytical model for air‐gap magnetic field prediction of surface‐mounted permanent magnet motors with a quasi‐regular polygon rotor
Escobar et al. Continuous-domain semi-analytical method for tolerance analysis of axial flux permanent magnet machines
Wang et al. Design optimization and analysis of a segmented-rotor switched reluctance machine for BSG application in HEVs
Aguiar et al. A new approach for computation of magnetic core losses in large hydro electrical generator

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
TA01 Transfer of patent application right

Effective date of registration: 20201019

Address after: 400044 Shapingba street, Shapingba District, Chongqing, No. 174

Applicant after: Chongqing University

Applicant after: STATE GRID HUBEI ELECTRIC POWER Co.,Ltd.

Address before: 400044 Shapingba street, Shapingba District, Chongqing, No. 174

Applicant before: Chongqing University

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant