CN106844950A - 一种220kV三相变压器A相短路电磁场的计算方法及装置 - Google Patents

一种220kV三相变压器A相短路电磁场的计算方法及装置 Download PDF

Info

Publication number
CN106844950A
CN106844950A CN201710039957.0A CN201710039957A CN106844950A CN 106844950 A CN106844950 A CN 106844950A CN 201710039957 A CN201710039957 A CN 201710039957A CN 106844950 A CN106844950 A CN 106844950A
Authority
CN
China
Prior art keywords
phase
short circuit
subdivision
transformers
phase transformers
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.)
Pending
Application number
CN201710039957.0A
Other languages
English (en)
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.)
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Original Assignee
Electric Power Research Institute of Guangdong Power Grid 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 Electric Power Research Institute of Guangdong Power Grid Co Ltd filed Critical Electric Power Research Institute of Guangdong Power Grid Co Ltd
Priority to CN201710039957.0A priority Critical patent/CN106844950A/zh
Publication of CN106844950A publication Critical patent/CN106844950A/zh
Pending legal-status Critical Current

Links

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)
  • Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)

Abstract

本发明实施例公开了一种220kV三相变压器A相短路电磁场的计算方法及装置,用于解决由于目前的大型220kV变压器的模型比较复杂,材料的非线性特性难以获得,求解的计算量比较大,而导致的大型220kV变压器三维电磁场求解效率低的技术问题。本发明实施例方法包括:建立220kV三相变压器简化模型;将所述220kV三相变压器简化模型剖分为若干个三维网格剖分单元;对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析。

Description

一种220kV三相变压器A相短路电磁场的计算方法及装置
技术领域
本发明涉及电力技术领域,尤其涉及一种220kV三相变压器A相短路电磁场的计算方法及装置。
背景技术
随着电力建设的不断发展,电力设备朝着大型化方向发展。随着变压器电压等级的大幅提升,作为电力系统主要载体之一的变压器的安全稳定运行,关系着电网的整体抗风险能力。如何保证变压器长期处于良好的工作状态一直是电力企业设备管理的重中之重。近些年,电网系统屡次发生因外部短路冲击造成主变重瓦斯保护跳闸、压力释放阀动作喷油等。当因外部电路短路而造成电流冲击时绕组导线中所通过的短路电流数值可达到而定数值的15~20倍。变压器在极短的时间之内将产生较大的电磁力和电磁损耗。在电磁力的作用下,绕组可能会发生振动,从而影响油箱中的油流运动,产生油流涌动。要准确的确定油箱中的油流运动,需要对变压器在短路状态产生的电磁场强度分布、绕组受到的电磁力以及油流的温升有精确的了解。
要对变压器中绕组的受力进行准确的计算,首先应对变压器短路时的电磁场进行精确的计算。因此,有限元法在很早就被引入到了变压器的磁场分析中。对于变压器的损耗计算,常有的计算手段是运用解析公式和半经验公式的时域法,正交分解合成法,但局限性太大。在采用计算机辅助分析方面,一般也都是采用2D有限元法。
这些方法都能在一定程度上解决工程问题。但是这些方法本身含有的简化太多,计算的结果适用性比较小,提供的数据难以支持进一步的变压器振动或油流涌动的分析。因此,针对此种大型220kV变压器的三维电磁场求解,进行分析的很少见。
究其原因,其大型220kV变压器的模型比较复杂,材料的非线性特性难以获得,求解的计算量比较大,从而导致了大型220kV变压器三维电磁场求解效率低的技术问题。
发明内容
本发明实施例提供了一种220kV三相变压器A相短路电磁场的计算方法及装置,解决了由于目前的大型220kV变压器的模型比较复杂,材料的非线性特性难以获得,求解的计算量比较大,而导致的大型220kV变压器三维电磁场求解效率低的技术问题。
本发明实施例提供的一种220kV三相变压器A相短路电磁场的计算方法,包括:
建立220kV三相变压器简化模型;
将所述220kV三相变压器简化模型剖分为若干个三维网格剖分单元;
对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析。
可选地,将所述220kV三相变压器简化模型剖分为若干个三维网格剖分单元具体包括:
将所述220kV三相变压器简化模型剖分为若干个四面体二阶单元。
可选地,对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析具体包括:
对若干个所述四面体二阶单元进行与A相对应的绕组短路处理,并设置外电路短路瞬态激励源;
对所述四面体二阶单元的棱边的矢量位自由度进行一阶元计算,对所述四面体二阶单元的节点的标量位自由度进行二阶元计算,求解出三维瞬态磁场;
对所述三维瞬态磁场进行计算结果分析。
可选地,对所述三维瞬态磁场进行计算结果分析具体包括:
根据求解出的所述三维瞬态磁场进行分析确定所述220kV三相变压器简化模型对应的内部磁场分布特性、各相绕组和铁芯欧姆损耗,以及各相绕组受力情况。
可选地,根据求解出的所述三维瞬态磁场进行分析还包括:
获取到在A相短路状态下所述220kV三相变压器简化模型对应的三相电流瞬态波形、预置时刻下所述220kV三相变压器简化模型对应的磁感应强度分布云图和在指定路径位置下考察得到的磁感应强度分布曲线图。
本发明实施例提供的一种220kV三相变压器A相短路电磁场的计算装置,其特征在于,包括:
建立单元,用于建立220kV三相变压器简化模型;
剖分单元,用于将所述220kV三相变压器简化模型剖分为若干个三维网格剖分单元;
分析单元,用于对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析。
可选地,剖分单元,具体用于将所述220kV三相变压器简化模型剖分为若干个四面体二阶单元。
可选地,分析单元包括:
短路处理子单元,用于对若干个所述四面体二阶单元进行与A相对应的绕组短路处理,并设置外电路短路瞬态激励源;
计算子单元,用于对所述四面体二阶单元的棱边的矢量位自由度进行一阶元计算,对所述四面体二阶单元的节点的标量位自由度进行二阶元计算,求解出三维瞬态磁场;
分析子单元,用于对所述三维瞬态磁场进行计算结果分析。
可选地,分析子单元包括:
第一分析模块,用于根据求解出的所述三维瞬态磁场进行分析确定所述变压器简化模型对应的内部磁场分布特性、各相绕组和铁芯欧姆损耗,以及各相绕组受力情况。
可选地,分析子单元还包括:
第二分析模块,用于根据求解出的所述三维瞬态磁场进行分析,获取到在A相短路状态下所述220kV三相变压器简化模型对应的三相电流瞬态波形、预置时刻下所述220kV三相变压器简化模型对应的磁感应强度分布云图和在指定路径位置下考察得到的磁感应强度分布曲线图。
从以上技术方案可以看出,本发明实施例具有以下优点:
本发明实施例提供的一种220kV三相变压器A相短路电磁场的计算方法及装置,包括:建立220kV三相变压器简化模型;将所述220kV三相变压器简化模型剖分为若干个三维网格剖分单元;对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析,解决了由于目前的大型220kV变压器的模型比较复杂,材料的非线性特性难以获得,求解的计算量比较大,而导致的大型220kV变压器三维电磁场求解效率低的技术问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。
图1为本发明实施例提供的一种220kV三相变压器A相短路电磁场的计算方法的一个实施例的流程示意图;
图2为本发明实施例提供的一种220kV三相变压器A相短路电磁场的计算方法的另一个实施例的流程示意图;
图3为一个铁芯模型三维图;
图4为另一个铁芯模型三维图;
图5为变压器模型图;
图6为铁芯BH曲线图;
图7为铁芯BP曲线图;
图8为A相短路瞬态电流波形图;
图9为铁芯网格剖分示意图;
图10为铁芯、线圈、拉板及夹件网格剖分示意图;
图11为整体网格剖分正视图;
图12为网格统计量示意图;
图13为二阶四面体网格示意图;
图14为A相短路时高压侧电压及电流测试波形示意图;
图15为A相高中压短路时外电路原理图;
图16为A相短路时高压侧电流波形示意图;
图17为A相短路时中压电流波形示意图;
图18为A相短路某时刻磁感应强度分布云图;
图19为A相短路某时刻磁感应强度矢量图;
图20为A相短路时变压器铁芯场强分布云图;
图21为A相短路时变压器铁芯磁矢量分布图;
图22为本发明实施例提供的一种220kV三相变压器A相短路电磁场的计算装置的一个实施例的结构示意图;
图23为本发明实施例提供的一种220kV三相变压器A相短路电磁场的计算装置的另一个实施例的结构示意图。
具体实施方式
本发明实施例提供了一种220kV三相变压器A相短路电磁场的计算方法及装置,解决了由于目前的大型220kV变压器的模型比较复杂,材料的非线性特性难以获得,求解的计算量比较大,而导致的大型220kV变压器三维电磁场求解效率低的技术问题。
为使得本发明的发明目的、特征、优点能够更加的明显和易懂,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,下面所描述的实施例仅仅是本发明一部分实施例,而非全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
本次电磁计算采用ANSYS Maxwell R17.0低频电磁场有限元分析软件进行仿真,ANSYS Maxwell R17.0集成于ANSYS Electronic Desktop R17.0电子桌面。ANSYS Maxwell是全球最著名的电磁场有限元分析软件,广泛应用于电机、作动器、高压电器设备、高压开关等机电产品和电力产品的电磁场计算和方案优化。
请参阅图1,本发明实施例提供的一种220kV三相变压器A相短路电磁场的计算方法的一个实施例包括:
101、建立220kV三相变压器简化模型;
本实施例中,当需要对三相变压器三维电磁场进行分析时,首先获取到通过Maxwell工具中的变压器基元模型建立的220kV三相变压器简化模型。
102、将220kV三相变压器简化模型剖分为若干个三维网格剖分单元;
当获取到通过Maxwell工具中的变压器基元模型建立的220kV三相变压器简化模型,需要通过Maxwell工具的静态自适应网格将220kV三相变压器简化模型剖分为若干个三维网格剖分单元。
103、对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析。
当通过Maxwell工具的静态自适应网格将变压器简化模型剖分为若干个三维网格剖分单元之后,需要通过Maxwell工具对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析。
请参阅图2,本发明实施例提供的一种三相变压器三维电磁场的计算方法的另一个实施例包括:
201、建立220kV三相变压器简化模型;
本实施例中,当需要对三相变压器三维电磁场进行分析时,首先获取到通过Maxwell工具中的变压器基元模型建立的220kV三相变压器简化模型。
202、将220kV三相变压器简化模型剖分为若干个四面体二阶单元;
当获取到通过Maxwell工具中的变压器基元模型建立的220kV三相变压器简化模型,需要通过Maxwell工具的静态自适应网格将220kV三相变压器简化模型剖分为若干个四面体二阶单元。
203、对若干个四面体二阶单元进行与A相对应的绕组短路处理,并设置外电路短路瞬态激励源;
当通过Maxwell工具的静态自适应网格将220kV三相变压器简化模型剖分为若干个四面体二阶单元之后,需要通过Maxwell工具对若干个四面体二阶单元进行与A相对应的绕组短路处理,并设置外电路短路瞬态激励源。
204、对四面体二阶单元的棱边的矢量位自由度进行一阶元计算,对四面体二阶单元的节点的标量位自由度进行二阶元计算求解出三维瞬态磁场;
当通过Maxwell工具对若干个四面体二阶单元进行对应绕组短路处理之后,需要对四面体二阶单元的棱边的矢量位自由度进行一阶元计算,对四面体二阶单元的节点的标量位自由度进行二阶元计算求解出三维瞬态磁场。
205、根据求解出的三维瞬态磁场进行分析确定220kV三相变压器简化模型对应的内部磁场分布特性、各相绕组和铁芯欧姆损耗,以及各相绕组受力情况;
求解出三维瞬态磁场之后,需要根据求解出的三维瞬态磁场进行分析确定220kV三相变压器简化模型对应的内部磁场分布特性、各相绕组和铁芯欧姆损耗,以及各相绕组受力情况。
206、获取到在A相短路状态下220kV三相变压器简化模型对应的三相电流瞬态波形、预置时刻下220kV三相变压器简化模型对应的磁感应强度分布云图和在指定路径位置下考察得到的磁感应强度分布曲线图。
求解出三维瞬态磁场之后,需要获取到短路状态下变压器简化模型对应的三相电流瞬态波形、预置时刻下变压器简化模型对应的磁感应强度分布云图和在指定路径位置下考察得到的磁感应强度分布曲线图。
下面以一具体应用场景进行详细的描述,应用例包括:
变压器高、中、低压非短路相绕组均为开路,调压线圈未通电流,单相(取A相)高中压绕组短路,进行瞬态磁场分析。
电力变压器等级110KV,额定容量150MVA,铭牌和结构参数见表1。
表1
变压器内部铁芯采用冷轧硅钢片叠压制成,计算采用Maxwell的建模功能对仿真对象(即220kV三相变压器)完成模型的建立。线圈采用饼式结构,结构参数如表1,计算采用Maxwell自带变压器线圈基元模型,模型如图3和图4所示,首先对变压器铁芯及线圈进行建模。
其他结构,如拉板、夹件、油箱结构参数如表2所示。
表2
基于Maxwell建立的变压器整体模型如图5所示。
材料参数见表3:
表3
变压器发生A相高中压绕组短路时,其激励源的瞬态电流波形如图8所示。
在三维瞬态场中仍旧采用的是T-Ω算法,对于低频瞬态磁场,麦克斯韦方程组可以写为式(1)—式(3)的形式。
ANSYS Maxwell软件磁场有限元计算分析满足下列Maxwell方程:
B=μH (3)
在式(1)—式(3)的基础上,可以构造出两个恒等式,如式(4),(5)所示。
Maxwell在求解三维瞬态磁场时,其棱边上的矢量位自由度采用了一阶元计算,而节点上的标量位自由度采用二阶元计算。
在三维瞬态磁场中,可以调用电压源或电流源作为模型激励源,而绕组又分为两种,一种是绞线型绕组,一种是实体绕组。其中绞线绕组不考虑涡流分布,认为绕组内的电密是完全均匀的,而实体绕组则需要计算其趋肤效应。在施加电压源时,绕组上的电流并不知道,所以需要进行绕组回路上的电压计算。
对于绞线型绕组,因其电阻值是一个集中参数,所以在外加电压源的时候,可以直接由电压源和直流电阻计算得到电流的数值,而对于实体线圈,因为其电阻值与频率、材料等有关,所以该电阻也称作交流电阻。在对实体线圈施加电压源的时候其交流电阻值并不知道,Maxwell是先按照式(6)进行计算;
Vi=∫∫∫RiJ0i(E+v×B)dR (6)
其中,J0i是第i个回路上的电密,同样包含X、Y、Z三个方向上的分量。
如果是变化的磁场,那么还应该计算绕组反电势,反电势可以按照式(7)计算。
Ei=∫∫∫RiHi·BidR (7)
在三维瞬态磁场中,瞬态包含两个含义,一个是瞬态的电磁过程,一个是瞬态的机械过程。由于此次计算未涉及在机械瞬态过程,故在此不再对处理方式进行说明。
电磁场的分析和计算通常归结为求微分方程的解。对于常微分方程,只要由辅助条件决定任意常数之后,其解就成为唯一的。对于偏微分方程,使其能成为唯一的辅助条件可以分为两种,一种是表达场的边界所处的物理情况,称为边界条件;一种种是确定场的初始状态,称为初始条件。边界条件和创始条件合称为定解条件。泛定方程和定解条件作为一个整体,称为定解问题。电磁场求解过程中有各种各样的边界条件本次求解过程中,具体应用了以下两类:
(1)诺依曼边界条件
其中,τ为诺依曼边界;n为边界Γ的外法线矢量。
f(τ)和h(τ)为一般函数,可以为常数和零,当为零时称齐次诺依曼边界条件。一般电磁场问题中将诺依曼边界条件称为第二类边界条件,它规定了边界处势的法向导数分布,在Maxwell计算中提到的齐次诺依曼边界条件,即法向导数为零,为默认边界条件,不需要用户指定。
(2)自然边界条件
媒质分解面上的边界条件,即不同媒质交界面场量的切向和法向边界条件属于自然边界条件,在Maxwell计算中是系统的默认边界条件,不需要用户指定。
整个模型统一将工作线圈均考虑在内,出于工程分析的角度出发,本次分析对模型进行了简化,简化的假设条件有:
(1)整个模型建模时对局部细节进行简化;
(2)模型结构参数不明确处由经验建模;
(3)铁芯材料导磁特性BH曲线与损耗特性BP曲线均由提供资料图片读图描点得到;
(4)变压器单相高中压短路瞬态激励由外电路提供,某些参数需根据计算进行调整;
(5)除短路相绕组激励外,其余激励源均设为开路;
(6)由于线圈饼数非常多,大大增加了后续流体计算的复杂程度,所以在结果提取时对每相各等级绕组做了简化求和处理。
这些假设或简化对分析的精确度是有所影响,由于这些因素带来的误差不予计及。
分析过程的重难点解决方法如下:
(1)电磁流体一体化耦合
此项目目的是在变压器发生单相高中压短路时,对其内部油流影响的分析,需要将发生短路时各相绕组各饼所受电磁力作为油流分析的源条件。各饼受力可以在流体仿真软件Fluent中直接赋值,然后进行计算得到相应的油流分析结果。但是由于变压器的结构比较复杂,运行时的实际工况与直接添加的电磁力可能存在比较大的偏差,导致油流分析的结果不准确。
此时需要采用多物理场耦合功能,利用电磁仿真软件Maxwell计算变压器线圈上的电动力,并将其直接映射到Fluent软件,计算对变压器内部油流分布的影响。需要特别指出的是,基于ANSYS Workbench进行多物理场耦合时,所有的模型、仿真结果等信息在不同软件之前进行传递时,均通过Workbench的数据映射功能自动实现,而无需再通过其他第三方软件,从而保证了工作的效率,以及数据传递的准确性。
(2)模型处理
由于变压器的外形尺寸较大,内部结构复杂,且未提供仿真模型。变压器全部模型均从零开始建立,过程中采用了Maxwell自带基元对铁芯和高、中、低及调压线圈建模,大大提高了工作效率。
(3)材料属性、工况条件等前处理工作
计算所涉及材料属性,如铁芯导磁特性、损耗特性均只提供了图片,读图工作不仅繁琐,并且对精度要求很高。计算采用Maxwell自带的图片捕捉工具,直接读取电子图片中曲线,省去纸质手动描点,大大简化了工作繁琐程度。但需说明,此工具仍与工程师操作、图片质量,有较大依赖性,读图工作得到数据造成的误差,仿真计算不予考虑。
短路瞬态激励由外电路提供,某些参数视需根据计算结果进行调整所得瞬态电流施加于绕组,以更加真实的仿真实际情况。
(4)网格剖分
Maxwell是基于有限元算法,有限元法是将整个求解区域离散化,分割成许多小的区域,称之为“单元”或“有限元”。由此,网格剖分的质量将直接影响计算结果的精度。而变压器外形尺寸较大,内部结构复杂,给网格剖分带来极大挑战。且网格剖分工作较难把握,尤其对于经验不是很丰富的工程师。
本次计算采用Maxwell静态场具备的自动自适应网格剖分功能与动态网格导入功能,将静态自适应网格应用到瞬态磁场中进行计算。这种方法不仅提高工作效率,且大大降低了仿真精度对工程师与软件操作熟练程度的依赖性,提高了可操作性与计算结果的可信度。图9-图12为网剖结果展示。
(5)计算量
Maxwell是基于有限元算法,四分之一模型剖分网格量在93万左右,而且Maxwell所采用的三维网格剖分单元是有限元理论中最为稳定的四面体二阶单元,即每个网格除了将其4个顶点作为节点外,同时还将每条棱中心也作为其计算节点,每个网格共有10个节点,如图13所示,所以计算精度很高。
另外,由于本次电磁计算的主要目的是对变压器各相绕组各饼受力的分析,所以在基本计算量的基础上,添加了15多个受力计算矩阵,显然,计算量非常庞大。
采用eHPC并行计算,进行求解,计算使用了10个核同时运行。计算采用一台12核,64G内存配置的服务器进行求解。
(6)结果提取
由于变压器模型复杂,各相各等级绕组饼数较多,每饼受力数据为140多组,计算量非常多,且大大增加了后续流体计算中的复杂程度。
最终决定结果提取时对每相各等级绕组做了简化求和处理。
变压器单相高中压短路瞬态激励由外电路提供,某些参数需根据计算进行调整,其中,中压绕组工作时起调压作用,调压线圈不在工作状态。工况为单相高、中压绕组短路,A相高、中压绕组短路时高压线圈电压、电流如图14所示,B、C相绕组设为开路即电流为0,外电路搭建如图15所示。
对变压器单相高中压绕组短路瞬态过程进行电磁分析,计算变压器单相高、中压绕组短路时,内部磁场分布特性、各相绕组和铁芯欧姆损耗,以及各相绕组受力情况。
(1)激励波形
变压器在发生单相短路时,高压侧A、B、C三相电流瞬态波形如图16所示,其中波形为短路相瞬态电流。短路相高、中压侧电流波形如图17所示。
(2)场量结果
整个结果可以显示任何时刻变压器全部结构的磁场分布,图18-图21分别是不同相高中压短路时,同一时刻变压器内部磁感应强度分布云图和矢量图。
请参阅图22,本发明实施例中提供的一种220kV三相变压器A相短路电磁场的计算装置的一个实施例包括:
建立单元301,用于建立220kV三相变压器简化模型;
剖分单元302,用于将220kV三相变压器简化模型剖分为若干个三维网格剖分单元;
分析单元303,用于对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析。
本实施例中,建立单元301通过获取到通过Maxwell工具中的变压器基元模型建立的220kV三相变压器简化模型,剖分单元302通过Maxwell工具的静态自适应网格将220kV三相变压器简化模型剖分为若干个三维网格剖分单元,分析单元303通过Maxwell工具对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析,解决了由于目前的大型220kV变压器的模型比较复杂,材料的非线性特性难以获得,求解的计算量比较大,而导致的大型220kV变压器三维电磁场求解效率低的技术问题。
上面是对三相变压器A相短路状态下电磁场的计算装置的各单元进行详细的描述,下面将对子单元进行详细的描述,请参阅图23,本发明实施例中提供的一种220kV三相变压器A相短路电磁场的计算装置的另一个实施例包括:
建立单元401,用于建立220kV三相变压器简化模型;
剖分单元402,用于将220kV三相变压器简化模型剖分为若干个三维网格剖分单元,剖分单元402,具体用于将220kV三相变压器简化模型剖分为若干个四面体二阶单元;
分析单元403,用于对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析。
分析单元403具体包括:
短路处理子单元4031,用于对若干个所述四面体二阶单元进行与A相对应的绕组短路处理,并设置外电路短路瞬态激励源;
计算子单元4032,用于对四面体二阶单元的棱边的矢量位自由度进行一阶元计算,对四面体二阶单元的节点的标量位自由度进行二阶元计算,求解出三维瞬态磁场;
分析子单元4033,用于对三维瞬态磁场进行计算结果分析。
分析子单元4033具体包括:
第一分析模块4033a,用于根据求解出的三维瞬态磁场进行分析确定变压器简化模型对应的内部磁场分布特性、各相绕组和铁芯欧姆损耗,以及各相绕组受力情况;
第二分析模块4033b,用于根据求解出的三维瞬态磁场进行分析,获取到在A相短路状态下220kV三相变压器简化模型对应的三相电流瞬态波形、预置时刻下所述220kV三相变压器简化模型对应的磁感应强度分布云图和在指定路径位置下考察得到的磁感应强度分布曲线图。
本实施例中,建立单元401通过获取到通过Maxwell工具中的变压器基元模型建立的220kV三相变压器简化模型,剖分单元402通过Maxwell工具的静态自适应网格将220kV三相变压器简化模型剖分为若干个三维网格剖分单元,分析单元403通过Maxwell工具对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析,解决了由于目前的大型220kV变压器的模型比较复杂,材料的非线性特性难以获得,求解的计算量比较大,而导致的大型220kV变压器三维电磁场求解效率低的技术问题。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统,装置和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在本申请所提供的几个实施例中,应该理解到,所揭露的系统,装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-OnlyMemory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种220kV三相变压器A相短路电磁场的计算方法,其特征在于,包括:
建立220kV三相变压器简化模型;
将所述220kV三相变压器简化模型剖分为若干个三维网格剖分单元;
对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析。
2.根据权利要求1所述的220kV三相变压器A相短路电磁场的计算方法,其特征在于,将所述220kV三相变压器简化模型剖分为若干个三维网格剖分单元具体包括:
将所述220kV三相变压器简化模型剖分为若干个四面体二阶单元。
3.根据权利要求2所述的220kV三相变压器A相短路电磁场的计算方法,其特征在于,对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析具体包括:
对若干个所述四面体二阶单元进行与A相对应的绕组短路处理,并设置外电路短路瞬态激励源;
对所述四面体二阶单元的棱边的矢量位自由度进行一阶元计算,对所述四面体二阶单元的节点的标量位自由度进行二阶元计算,求解出三维瞬态磁场;
对所述三维瞬态磁场进行计算结果分析。
4.根据权利要求3所述的220kV三相变压器A相短路电磁场的计算方法,其特征在于,对所述三维瞬态磁场进行计算结果分析具体包括:
根据求解出的所述三维瞬态磁场进行分析确定所述220kV三相变压器简化模型对应的内部磁场分布特性、各相绕组和铁芯欧姆损耗,以及各相绕组受力情况。
5.根据权利要求4所述的220kV三相变压器A相短路电磁场的计算方法,其特征在于,根据求解出的所述三维瞬态磁场进行分析还包括:
获取到在A相短路状态下所述220kV三相变压器简化模型对应的三相电流瞬态波形、预置时刻下所述220kV三相变压器简化模型对应的磁感应强度分布云图和在指定路径位置下考察得到的磁感应强度分布曲线图。
6.一种220kV三相变压器A相短路电磁场的计算装置,其特征在于,包括:
建立单元,用于建立220kV三相变压器简化模型;
剖分单元,用于将所述220kV三相变压器简化模型剖分为若干个三维网格剖分单元;
分析单元,用于对若干个所述三维网格剖分单元进行与A相对应的绕组短路处理,并设置激励源,进行电磁瞬态分析。
7.根据权利要求6所述的220kV三相变压器A相短路电磁场的计算装置,其特征在于,剖分单元,具体用于将所述220kV三相变压器简化模型剖分为若干个四面体二阶单元。
8.根据权利要求7所述的220kV三相变压器A相短路电磁场的计算装置,其特征在于,分析单元包括:
短路处理子单元,用于对若干个所述四面体二阶单元进行与A相对应的绕组短路处理,并设置外电路短路瞬态激励源;
计算子单元,用于对所述四面体二阶单元的棱边的矢量位自由度进行一阶元计算,对所述四面体二阶单元的节点的标量位自由度进行二阶元计算,求解出三维瞬态磁场;
分析子单元,用于对所述三维瞬态磁场进行计算结果分析。
9.根据权利要求8所述的220kV三相变压器A相短路电磁场的计算装置,其特征在于,分析子单元包括:
第一分析模块,用于根据求解出的所述三维瞬态磁场进行分析确定所述变压器简化模型对应的内部磁场分布特性、各相绕组和铁芯欧姆损耗,以及各相绕组受力情况。
10.据权利要求9所述的220kV三相变压器A相短路电磁场的计算装置,其特征在于,分析子单元还包括:
第二分析模块,用于根据求解出的所述三维瞬态磁场进行分析,获取到在A相短路状态下所述220kV三相变压器简化模型对应的三相电流瞬态波形、预置时刻下所述220kV三相变压器简化模型对应的磁感应强度分布云图和在指定路径位置下考察得到的磁感应强度分布曲线图。
CN201710039957.0A 2017-01-18 2017-01-18 一种220kV三相变压器A相短路电磁场的计算方法及装置 Pending CN106844950A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710039957.0A CN106844950A (zh) 2017-01-18 2017-01-18 一种220kV三相变压器A相短路电磁场的计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710039957.0A CN106844950A (zh) 2017-01-18 2017-01-18 一种220kV三相变压器A相短路电磁场的计算方法及装置

Publications (1)

Publication Number Publication Date
CN106844950A true CN106844950A (zh) 2017-06-13

Family

ID=59124195

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710039957.0A Pending CN106844950A (zh) 2017-01-18 2017-01-18 一种220kV三相变压器A相短路电磁场的计算方法及装置

Country Status (1)

Country Link
CN (1) CN106844950A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107273614A (zh) * 2017-06-15 2017-10-20 广东电网有限责任公司电力科学研究院 220kV变压器绕组受到短路电力时受力计算方法
CN107423528A (zh) * 2017-08-29 2017-12-01 广东电网有限责任公司电力科学研究院 一种电力变压器绕组短路瞬态受力分析方法及装置

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105550435A (zh) * 2015-12-13 2016-05-04 东北电力大学 一种基于变压器分段模型的绕组漏磁与电动力计算方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105550435A (zh) * 2015-12-13 2016-05-04 东北电力大学 一种基于变压器分段模型的绕组漏磁与电动力计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
谢道伟: ""油浸式三相电力变压器电磁场及温度场的瞬态仿真分析"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *
陈孝荣: ""自供电磁流变减振系统研究"", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107273614A (zh) * 2017-06-15 2017-10-20 广东电网有限责任公司电力科学研究院 220kV变压器绕组受到短路电力时受力计算方法
CN107273614B (zh) * 2017-06-15 2020-07-07 广东电科院能源技术有限责任公司 220kV变压器绕组受到短路电力时受力计算方法
CN107423528A (zh) * 2017-08-29 2017-12-01 广东电网有限责任公司电力科学研究院 一种电力变压器绕组短路瞬态受力分析方法及装置

Similar Documents

Publication Publication Date Title
US5692158A (en) Methods for generating models of non-linear systems and components and for evaluating parameters in relation to such non-linear models
Annakkage et al. Dynamic system equivalents: A survey of available techniques
CN106407499A (zh) 一种三相变压器三维电磁场的计算方法及装置
Liu et al. A real-time nonlinear hysteretic power transformer transient model on FPGA
CN106055812A (zh) 一种三维瞬态电磁场数值模拟分析方法及装置
CN107423528B (zh) 一种电力变压器绕组短路瞬态受力分析方法及装置
CN104200055B (zh) 特高压变压器空载合闸情况下的励磁涌流仿真方法及装置
CN101794321A (zh) 考虑励磁阻抗非线性影响的单相三绕组自耦变压器模型
CN102609575A (zh) 一种基于隐式数值积分的电力系统暂态稳定仿真方法
CN106407576A (zh) 一种设备电磁辐射特性参数的提取方法和系统
CN104698235A (zh) 继电保护测试中暂态故障波形的生成方法
CN106650007A (zh) 变压器稳态时漏磁场及绕组振动位移计算方法
CN103018534A (zh) 确定谐波电压的方法及系统
CN108291939A (zh) 低磁芯损耗内置磁体电机设计
CN103500245A (zh) 一种基于多回路法的场路瞬态-瞬态耦合仿真方法
CN110135050A (zh) 一种气体绝缘变电站空间电磁场计算方法
CN106844950A (zh) 一种220kV三相变压器A相短路电磁场的计算方法及装置
CN106682362A (zh) 一种220kV三相变压器C相短路电磁场的计算方法及装置
CN109145447A (zh) 一种可控磁饱和电抗器建模方法
CN106446372A (zh) 一种三相变压器a相短路状态下电磁场的计算方法及装置
CN104679936B (zh) 永磁饱和型故障限流器磁流热耦合建模方法
CN103258124A (zh) 一种高压变电站电磁设备工频磁场评估方法
CN105243272A (zh) 一种变压器绕组结构电感计算方法
CN106815436A (zh) 一种220kV三相变压器B相短路电磁场的计算方法及装置
CN106372351A (zh) 一种三相变压器c相短路状态下电磁场的计算方法及装置

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20170613