CN109738709A - 一种大型汽轮发电机端部电磁场、温度场计算方法 - Google Patents
一种大型汽轮发电机端部电磁场、温度场计算方法 Download PDFInfo
- Publication number
- CN109738709A CN109738709A CN201811416189.7A CN201811416189A CN109738709A CN 109738709 A CN109738709 A CN 109738709A CN 201811416189 A CN201811416189 A CN 201811416189A CN 109738709 A CN109738709 A CN 109738709A
- Authority
- CN
- China
- Prior art keywords
- model
- temperature
- field
- turbogenerator
- calculated
- 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.)
- Granted
Links
Landscapes
- Control Of Eletrric Generators (AREA)
- Manufacture Of Motors, Generators (AREA)
Abstract
本发明公开了一种大型汽轮发电机端部电磁场、温度场计算方法,其特征在于,包括建立计算模型、建立端部电磁场分布计算模块、建立温度场分布计算模块、结果分析等步骤。本发明通过电‑热顺序耦合分析技术计算,使设计人员能够预测发电机端部磁场及端部结构件的温升是否保持在绝缘材料和金属材料所允许的限度之内,从而验证设计结构以及材料选取的合理性;另外,本发明计算过程模块化,易于使用便于修改,可重复调用,顺序耦合计算模式有效降低了求解方程的规模,保证计算精度的情况下可使计算速度大为提高。
Description
技术领域
本发明涉及一种大型汽轮发电机端部电磁场温度场的数值计算方法,可用于大型汽轮发电机定子端部结构件的温升计算及结构设计研究中,属于发电机部件电磁场及温升计算方法技术领域。
背景技术
大型汽轮发电机中各部件的温升对机组的运行性能有很大影响,同时也是限制发电机容量增加的一个重要原因。随着发电机单机容量的不断增加,定子线负荷已超过每厘米2000安以上,这将导致发电机端部区域的磁密增加;另外,由于电网需求,发电机可能运行于深度进相工况,这也使得发电机端部磁场会更为集中,端部结构件中的感应涡流有所增大,损耗可能会有所增加,这不仅影响发电机的效率,还可能造成发电机端部结构件的局部过热,影响发电机运行的可靠性。因此,在开发之初需对大型汽轮发电机端部磁场、结构件的温度场进行分析研究,确保其温升能够保持在所选绝缘材料和金属材料允许的限度之内。
发明内容
本发明所要解决的问题是:提供一种通用的在稳态运行工况下汽轮发电机端部电磁场、端部结构件温度场的计算方法。
为了解决上述问题,本发明的技术方案是提供了一种大型汽轮发电机端部电磁场、温度场计算方法,其特征在于,包括以下步骤:
步骤1:建立计算模型:
按照发电机端部线圈及结构件的图纸建立模型;
步骤2:建立端部电磁场分布计算模块:
计算发电机端部电磁场及涡流在发电机端部结构中的分布情况,计算焦耳热分布,取焦耳热分布中设定时长的瞬时损耗进行平均,保存计算得到的平均损耗结果;
步骤3:建立温度场分布计算模块:
在电流场分布计算的基础上计算发电机端部结构件温度场的分布情况,将计算得到的温度分布结果进行保存;
步骤4:结果分析:
通过结果处理获取压圈结构上具体的温度分布状况,进而得到各部分的温升情况。
优选地,所述步骤1中,按照发电机端部线圈及结构件的图纸采用CAD计算机辅助设计软件NX建立三维模型。
更优选地,所述三维模型包括部分铁心、阶梯段铁心、压圈、磁屏蔽、端部线圈、转子结构的三维模型,将三维模型导入CAE计算机辅助分析软件Infolytica或Magnet中,并运行于计算机上。
优选地,所述步骤2包括以下步骤:
步骤2.1:给定端部线棒、压圈、边段铁心、转子线圈、转轴结构的电阻率ρ及相对磁导率μr,包括设定自由空间的相对磁导率为1;
步骤2.2:根据模型结构中的不同部件,用四面体电磁单元对步骤1建立的模型进行网格划分,将发电机端部结构模型离散为数值模型;
步骤2.3:在铁心端面、自由空间的外边界施加磁力线平行条件,即设定相应方向的矢量磁位A=0以模拟磁场有效界限;
步骤2.4:按三相划分定子绕组,建立外部电路,采用外部电路法给模型中A、B、C三相的每根定子线圈导体施加具体的电流;
步骤2.5:采用外部电路法给转子绕组施加直流励磁电流;
步骤2.6:设定求解方式为瞬态旋转法,转子按步进方法进行旋转,设定计算时长及步长,以此得到瞬态的端部磁场,并计算得到压圈焦耳热qv,取最后2个周期的瞬时损耗进行平均,将平均损耗结果保存。
优选地,所述步骤3包括:
步骤3.1:根据需要计算的结构件给定材料属性及冷却参数,给定压圈的导热系数;
步骤3.2:设定压圈的热边界面,根据各边界面的风速、风温计算出各表面传热系数为h,在热边界面设定传热系数h;
步骤3.3:根据热传导微分方程进行求解,式中,T为待求的温度,▽为向量微分算子,λcu为导热系数,将计算得到的温度分布结果进行保存。
本发明通过电-热顺序耦合分析技术计算,使设计人员能够预测发电机端部磁场及端部结构件的温升是否保持在绝缘材料和金属材料所允许的限度之内,从而验证设计结构以及材料选取的合理性;另外,本发明计算过程模块化,易于使用便于修改,可重复调用,顺序耦合计算模式有效降低了求解方程的规模,保证计算精度的情况下可使计算速度大为提高。
附图说明
图1为本发明提供的计算方法的流程图。
具体实施方式
为使本发明更明显易懂,兹以优选实施例,并配合附图作详细说明如下。
实施例
一种大型汽轮发电机端部电磁场、温度场计算方法。本发明的计算过程采用模块化,将发电机端部电磁场、损耗场的分布计算与温度场的分布计算分模块进行,采用顺序耦合的方式,电磁场及损耗计算结果保存后供温度场计算调用(计算流程如图1所示),避免了电磁场与温度场同时耦合计算时由于自由度多而导致的求解方程庞大、计算缓慢、计算程序不易重复调用、程序不易阅读等缺点,具体步骤为:
步骤1、建立计算模型
按照发电机端部线圈及结构件的图纸采用CAD计算机辅助设计软件NX建立模型,包括部分铁心、阶梯段铁心、压圈、磁屏蔽、端部线圈、转子结构的三维模型,并将CAD模型导入CAE计算机辅助分析软件Infolytica或Magnet内,运行于HP工作站,即计算机上。
步骤2、建立端部电磁场分布计算模块
计算端部电磁场及涡流在端部结构中的分布情况,计算焦耳热分布,采用Infolytica/Magnet软件的交互式界面进行设置,并采用VB脚本语言对损耗结果的处理进行模块化编程,具体步骤为:
步骤2.1、给定端部线棒、压圈、边段铁心、转子线圈、转轴结构的电阻率ρ及相对磁导率μr,包括设定自由空间的相对磁导率为1;
步骤2.2、根据模型结构中的不同部件,用四面体电磁单元对步骤1建立的计算模型进行网格划分,将发电机端部结构模型离散为数值模型;
步骤2.3、在铁心端面、自由空间的外边界施加磁力线平行条件,即设定相应方向的矢量磁位A=0以模拟磁场有效界限;
步骤2.4、按三相划分定子绕组,建立外部电路,采用外部电路法给模型中A、B、C三相的每根定子线圈导体施加具体的电流;
步骤2.5、采用外部电路法给转子绕组施加直流励磁电流;
步骤2.6、由于端部结构为非对称结构,设定求解方式为瞬态旋转法,转子按步进方法进行旋转;
步骤2.7、设定计算时长及步长,由于发电机是稳态运行,计算时长设定为6~10周期即可达到稳定,步长按2°电角度设定;
步骤2.8、以此得到瞬态的端部磁场,并计算得到压圈焦耳热qv,取最后2个周期的瞬时损耗进行平均,将平均损耗结果保存,为后续的温度场计算提供数据。
步骤3、建立温度场分布计算模块
在电流场分布计算的基础上计算端部结构件温度场的分布情况,具体步骤为:
步骤3.1、根据需要计算的结构件(压圈、齿压板等)给定材料属性及冷却参数:给定压圈的导热系数;
步骤3.2、设定压圈的热边界面,根据各边界面的风速、风温计算出各表面传热系数为h,在热边界面设定传热系数h;
步骤3.3、根据热传导微分方程进行求解,式中,T即为待求的温度,▽为向量微分算子,λcu为导热系数,将计算得到的温度分布结果进行保存;
步骤4、结果分析
计算完成后,通过结果处理获取压圈结构上具体的温度分布状况,进而得到各部分的温升情况。
Claims (5)
1.一种大型汽轮发电机端部电磁场、温度场计算方法,其特征在于,包括以下步骤:
步骤1:建立计算模型:
按照发电机端部线圈及结构件的图纸建立模型;
步骤2:建立端部电磁场分布计算模块:
计算发电机端部电磁场及涡流在发电机端部结构中的分布情况,计算焦耳热分布,取焦耳热分布中设定时长的瞬时损耗进行平均,保存计算得到的平均损耗结果;
步骤3:建立温度场分布计算模块:
在电流场分布计算的基础上计算发电机端部结构件温度场的分布情况,将计算得到的温度分布结果进行保存;
步骤4:结果分析:
通过结果处理获取压圈结构上具体的温度分布状况,进而得到各部分的温升情况。
2.如权利要求1所述的大型汽轮发电机端部电磁场、温度场计算方法,其特征在于,所述步骤1中,按照发电机端部线圈及结构件的图纸采用CAD计算机辅助设计软件NX建立三维模型。
3.如权利要求2所述的大型汽轮发电机端部电磁场、温度场计算方法,其特征在于,所述三维模型包括部分铁心、阶梯段铁心、压圈、磁屏蔽、端部线圈、转子结构的三维模型,将三维模型导入CAE计算机辅助分析软件Infolytica或Magnet中,并运行于计算机上。
4.如权利要求1所述的大型汽轮发电机端部电磁场、温度场计算方法,其特征在于,所述步骤2包括以下步骤:
步骤2.1:给定端部线棒、压圈、边段铁心、转子线圈、转轴结构的电阻率ρ及相对磁导率μr,包括设定自由空间的相对磁导率为1;
步骤2.2:根据模型结构中的不同部件,用四面体电磁单元对步骤1建立的模型进行网格划分,将发电机端部结构模型离散为数值模型;
步骤2.3:在铁心端面、自由空间的外边界施加磁力线平行条件,即设定相应方向的矢量磁位A=0以模拟磁场有效界限;
步骤2.4:按三相划分定子绕组,建立外部电路,采用外部电路法给模型中A、B、C三相的每根定子线圈导体施加具体的电流;
步骤2.5:采用外部电路法给转子绕组施加直流励磁电流;
步骤2.6:设定求解方式为瞬态旋转法,转子按步进方法进行旋转,设定计算时长及步长,以此得到瞬态的端部磁场,并计算得到压圈焦耳热qv,取最后2个周期的瞬时损耗进行平均,将平均损耗结果保存。
5.如权利要求4所述的大型汽轮发电机端部电磁场、温度场计算方法,其特征在于,所述步骤3包括:
步骤3.1:根据需要计算的结构件给定材料属性及冷却参数,给定压圈的导热系数;
步骤3.2:设定压圈的热边界面,根据各边界面的风速、风温计算出各表面传热系数为h,在热边界面设定传热系数h;
步骤3.3:根据热传导微分方程进行求解,式中,T为待求的温度,▽为向量微分算子,λcu为导热系数,将计算得到的温度分布结果进行保存。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811416189.7A CN109738709B (zh) | 2018-11-26 | 2018-11-26 | 一种大型汽轮发电机端部电磁场、温度场计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811416189.7A CN109738709B (zh) | 2018-11-26 | 2018-11-26 | 一种大型汽轮发电机端部电磁场、温度场计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109738709A true CN109738709A (zh) | 2019-05-10 |
CN109738709B CN109738709B (zh) | 2021-01-05 |
Family
ID=66358229
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811416189.7A Active CN109738709B (zh) | 2018-11-26 | 2018-11-26 | 一种大型汽轮发电机端部电磁场、温度场计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109738709B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111062155A (zh) * | 2019-11-22 | 2020-04-24 | 上海电气电站设备有限公司 | 一种大型调相机端部钢制结构件的涡流损耗计算方法 |
CN113267689A (zh) * | 2021-04-21 | 2021-08-17 | 中国人民解放军海军工程大学 | 无线电能传输电、磁、热及温升测试系统、测试方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005347686A (ja) * | 2004-06-07 | 2005-12-15 | National Institute For Materials Science | 高臨界電流特性を有する超伝導材料及びその製造方法 |
CN103678835A (zh) * | 2014-01-15 | 2014-03-26 | 三峡大学 | 一种电机在电磁场-流场-温度场耦合计算中的建模方法 |
CN104008228A (zh) * | 2014-04-28 | 2014-08-27 | 北京交通大学 | 一种求解发电机定子端部绕组表面散热系数的新方法 |
CN104182580A (zh) * | 2014-08-21 | 2014-12-03 | 内蒙古工业大学 | 一种永磁风力发电机中涡流损耗分析方法 |
CN105245041A (zh) * | 2015-10-15 | 2016-01-13 | 北京交通大学 | 一种求解大型电机定转子通风沟表面散热系数的方法 |
CN106354971A (zh) * | 2016-09-14 | 2017-01-25 | 中国石油大学(华东) | 基于场路耦合分析的变压器状态预估新方法 |
CN107391822A (zh) * | 2016-08-20 | 2017-11-24 | 三峡大学 | 一种基于自适应时间步长的瞬变电磁‑温度场耦合计算方法 |
-
2018
- 2018-11-26 CN CN201811416189.7A patent/CN109738709B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005347686A (ja) * | 2004-06-07 | 2005-12-15 | National Institute For Materials Science | 高臨界電流特性を有する超伝導材料及びその製造方法 |
CN103678835A (zh) * | 2014-01-15 | 2014-03-26 | 三峡大学 | 一种电机在电磁场-流场-温度场耦合计算中的建模方法 |
CN104008228A (zh) * | 2014-04-28 | 2014-08-27 | 北京交通大学 | 一种求解发电机定子端部绕组表面散热系数的新方法 |
CN104182580A (zh) * | 2014-08-21 | 2014-12-03 | 内蒙古工业大学 | 一种永磁风力发电机中涡流损耗分析方法 |
CN105245041A (zh) * | 2015-10-15 | 2016-01-13 | 北京交通大学 | 一种求解大型电机定转子通风沟表面散热系数的方法 |
CN107391822A (zh) * | 2016-08-20 | 2017-11-24 | 三峡大学 | 一种基于自适应时间步长的瞬变电磁‑温度场耦合计算方法 |
CN106354971A (zh) * | 2016-09-14 | 2017-01-25 | 中国石油大学(华东) | 基于场路耦合分析的变压器状态预估新方法 |
Non-Patent Citations (1)
Title |
---|
夏海霞: "汽轮发电机端部电磁-温度耦合场的研究及其数值分析", 《浙江大学博士学位论文》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111062155A (zh) * | 2019-11-22 | 2020-04-24 | 上海电气电站设备有限公司 | 一种大型调相机端部钢制结构件的涡流损耗计算方法 |
CN113267689A (zh) * | 2021-04-21 | 2021-08-17 | 中国人民解放军海军工程大学 | 无线电能传输电、磁、热及温升测试系统、测试方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109738709B (zh) | 2021-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Simpson et al. | Additive manufacturing of shaped profile windings for minimal AC loss in electrical machines | |
Bardalai et al. | Reduction of winding AC losses by accurate conductor placement in high frequency electrical machines | |
Ge et al. | Fast iron loss and thermal prediction method for power density and efficiency improvement in switched reluctance machines | |
Zhang et al. | Thermal analysis and management for doubly salient brushless DC generator with flat wire winding | |
Kaehler et al. | Transient 3-D FEM computation of eddy-current losses in the rotor of a claw-pole alternator | |
Hsieh et al. | Design of brushless doubly-fed machines based on magnetic circuit modeling | |
Zhu et al. | Coupled electromagnetic-thermal-fluidic analysis of permanent magnet synchronous machines with a modified model | |
Di Leonardo et al. | Finite elements model co-simulation of an induction motor drive for traction application | |
CN109738709A (zh) | 一种大型汽轮发电机端部电磁场、温度场计算方法 | |
Zhao et al. | Exact modeling and multiobjective optimization of vernier machines | |
Liang et al. | Analysis of copper loss of permanent magnet synchronous motor with formed transposition winding | |
Zhiyong et al. | Thermal analysis of PMSM based on lumped parameter thermal network method | |
Chung et al. | Modeling of the armature slotting effect in the magnetic field distribution of a linear permanent magnet motor | |
Liang et al. | Analytical algorithm for strand end leakage reactance of transposition bar in AC machine | |
Zhao et al. | A generalized mesh-based thermal network model for SPM machines combining coupled winding solution | |
Carbonieri et al. | A fast and direct analysis of three-phase induction motors using finite element | |
Khan et al. | Improved performance design realization of a fractional kilowatt induction motor with predictive current control for water pumping | |
Golovanov et al. | Analytical methodology for modelling of circulating current loss in synchronous electrical machines with permanent magnets | |
Kefalas et al. | Finite element transient thermal analysis of PMSM for aerospace applications | |
Ho et al. | Reduction of numerical errors of time-stepping finite element analysis for dynamic simulation of electric machines | |
Qiu et al. | Analytical analysis of sleeve permeability for output performance of high speed permanent magnet generators driven by micro gas turbines | |
Handgruber et al. | Three-dimensional eddy current loss modeling in steel laminations of skewed induction machines | |
Mabhula et al. | Computational efficient parameter and performance prediction of wound-rotor induction motor | |
Stermecki et al. | Calculation of load-dependent equivalent circuit parameters of squirrel cage induction motors using time-harmonic FEM | |
Mabhula et al. | Gradient-based multi-objective design optimisation formulation of grid-connected wound-rotor induction motors |
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 |