CN106777812A - 一种变压器绕组短路振动特性的计算方法 - Google Patents
一种变压器绕组短路振动特性的计算方法 Download PDFInfo
- Publication number
- CN106777812A CN106777812A CN201710030273.4A CN201710030273A CN106777812A CN 106777812 A CN106777812 A CN 106777812A CN 201710030273 A CN201710030273 A CN 201710030273A CN 106777812 A CN106777812 A CN 106777812A
- Authority
- CN
- China
- Prior art keywords
- equivalent
- transformer
- model
- ansys
- short circuit
- 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
Classifications
-
- 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
- G06F30/23—Design 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)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种变压器绕组短路振动特性的计算方法,首先建立等效模型:用ANSYS建立变压器二维有限元模型,由于二维模型结构表示的不完整性,参数要做一定等效处理。根据材料特性的等效性原理,计算并赋予变压器各部分结构的等效参数。进行振动特性分析:采用ANSYS磁‑结构直接耦合技术,使用变压器二维等效模型,设置系统重力加速度和阻尼参数,并对模型施加边界条件;分别以各线饼的等效电流密度为激励,采用瞬态分析方法,多载荷步加载并求解;利用ANSYS后处理器,输出绕组振动过程中的位移、速度和加速度。本方法实现了对变压器绕组振动问题的降围,以及短路振动特性的准确动态求解。
Description
技术领域
本发明属于电力变压器短路特性计算技术领域,涉及一种变压器绕组短路振动特性的计算方法。
背景技术
电力变压器是电力系统中的重要设备之一,其运行状况直接影响着电网的稳定运行。近年来变压器事故统计显示,其抗短路能力不足,是威胁变压器安全运行的主要因素。其中,变压器突发短路故障时,绕组在短路力作用下,而产生剧烈振动,出现结构松动或变形,是造成电力变压器损坏和事故的重要原因。
绕组轴向振动特性的研究始于70年代初期,而后一直为各国学者所关注。但是由于变压器绕组结构复杂,以及当时研究工具的限制,对于绕组轴向振动特性的解析大多都作了很大程度的简化。
目前,对于变压器短路故障的轴向振动模型,大多分为以下两种:一种是采用经典的“质量-弹簧-阻尼”系统模型,该方法能够快速计算振动特性;另一种则是采用有限元模型。对于现有数值解法,有限元法已成为分析这类问题最有效的方法。其中,较多人采用变压器有限元模型,加载绕组短路力求解动态响应;也有少数人开始使用三维磁-结构耦合分析方法。
由于变压器短路引起的绕组振动是一个复杂的机电过程,对于经典的弹簧系统模型和实际仍存在一定差距;采用有限元模型加载短路力计算方法,仍忽略了电磁场和结构场间的动态相互影响;而三维磁-结构耦合分析方法,虽然较准确地解析了绕组短路振动过程,但由于计算占用资源巨大,工程应用过于复杂。
因此,寻找一种准确快速地计算变压器绕组短路振动特性的方法,为变压器设计和快速故障分析,提供了技术支持。
发明内容
本发明的目的在于提供一种变压器绕组短路振动特性的计算方法,解决了目前变压器绕组短路振动特性计算方法过于复杂,计算结果不准确的问题。
本发明所采用的技术方案是按照以下步骤进行:
步骤1:根据变压器尺寸和部件,在ANSYS中建立变压器二维有限元模型;
步骤2:由于二维模型的等效性,分别计算变压器各部分结构的等效属性;
(1)计算绕组线饼、垫块层和端圈等结构的等效密度,分别均匀等效处理;
(2)计算不同位置绝缘纸板在轴向预紧力作用下的弹性模量;
(3)计算各绝缘层的等效弹性模量;
步骤3:在ANSYS中定义并分配各结构的材料属性:电磁特性和结构特性;等效结构要赋予步骤2中计算的等效参数:等效密度、等效弹性模量;
步骤4:根据模型不同部位精细化程度的不同,分区域合理划分网格结构,并控制网格形状,避免出现畸变;
步骤5:定义轴向坐标系下的重力加速度g=9.8m/s2;考虑系统阻尼的影响,根据具体变压器阻尼大小,定义系统瑞利阻尼常数α和β;
步骤6:设置模型边界条件,轴对称边界施加磁通平行条件,固定件施加位移约束条件,分别对压板上边界和托板下边界施加轴向位移约束;
步骤7:分别以各线饼的等效电流密度作为激励,以时间历程定义载荷步;设置故障持续0.2s的一段动态载荷,按照1ms为步长,分别对高、低压绕组各线饼加载;
步骤8:采用瞬态分析方法,打开大形变效应和瞬态效应,划分载荷子步,设置写入结果文件的频率;
步骤9:利用ANSYS后处理器,输出故障过程中绕组振动的位移、速度和加速度随时间变化曲线。
进一步,步骤1中ANSYS中建立变压器二维有限元模型是指在ANSYS MechanicalAPDL软件环境下,选择PLANE13磁-结构耦合有限元单元,设置单元具有UX、UY、AZ自由度,直接耦合磁场和结构场,选取二维轴对称建模方式,根据变压器实际尺寸,建立二维轴对称模型。
进一步,步骤2中由于垫块、端圈、压板和压木等结构均采用绝缘纸板制成,在一定范围内,其应力-应变关系可表示为:
σ=aε+bε3 (1)
式中:σ为应力;ε为应变;a为线性常数;b为硬化系数。
通过实验可测得:a=105MPa;b=1750MPa。
在轴向预紧力作用下,其弹性模量可表示为:
压板和托板的弹性模量采用此参数直接定义。
考虑各相同性的两种材料,在载荷作用下产生相同应变,建立了Voigt
模型:
Eφ=φ1E1+φ2E2 (3)
垫块层和端部压木层采用Voigt模型等效。
考虑两相承受相同的应力,建立了Reuss模型:
端圈为多层层压式结构,具有两种结构属性,将其按结构特点纵向划分为几部分,采用Voigt模型和Reuss模型组合计算。
式中:Eφ为材料的等效弹性模量;E1为绝缘纸板的弹性模量;E2为空隙的弹性模量;φ1和φ2分别为绝缘纸板和空隙的体积分数,且φ1+φ2=1。
进一步,步骤4中模型整体网格尺寸设置为0.02,绕组及周围增大剖分网格数量,尺寸设置为0.015,整体采用四边形剖分。
进一步,步骤8中以0.1ms划分载荷子步,每隔5个子步保存一次计算结果。
进一步,利用ANSYS后处理器,输出故障过程中绕组振动的位移、速度和加速度随时间变化曲线。
本发明的有益效果是采用ANSYS二维有限元等效模型,实现了变压器绕组振动问题的降围;采用ANSYS磁-结构直接耦合分析技术,实现了短路振动特性的准确动态求解。
附图说明
图1为一种变压器模型结构图。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
本发明的计算方法,包括以下步骤:
1.根据变压器尺寸,建立变压器二维模型;
在ANSYS Mechanical APDL软件环境下,选择PLANE13磁-结构耦合有限元单元,设置单元具有UX、UY、AZ自由度,直接耦合磁场和结构场,选取二维轴对称建模方式。根据变压器结构尺寸,建立二维轴对称模型,包括铁心、高/低压绕组(线饼、垫块层)、端圈、压板等。如图1所示为一种变压器模型结构图。
2.由于二维模型的等效性,分别计算变压器各部分结构的等效属性。
1)计算变压器各部分结构的等效密度,包括线饼、垫块层和端圈等结构,对其分别均匀等效处理。
2)由于垫块、端圈、压板和压木等结构均采用绝缘纸板制成,属于非线性材料,计算不同位置绝缘纸板在轴向预紧力作用下的弹性模量。绝缘纸板在一定范围内,其应力-应变关系可表示为:
σ=aε+bε3 (1)
式中:σ为应力;ε为应变;a为线性常数;b为硬化系数。
通过实验可测得:a=105MPa;b=1750MPa。
在轴向预紧力作用下,其弹性模量可表示为:
压板和托板的弹性模量采用此参数直接定义。
3)计算二维等效模型其余绝缘层的等效弹性模量。
根据垫块层和端部绝缘的结构特点,可将其分别视作一种复合材料,即绝缘纸板和空隙两种材料的复合,对其采用Voigt模型和Reuss模型等效计算。
考虑各相同性的两种材料,在载荷作用下产生相同应变,建立了Voigt模型,可表示为:
Eφ=φ1E1+φ2E2 (3)
垫块层和端部压木层采用Voigt模型等效。
考虑两相承受相同的应力,建立了Reuss模型,可表示为:
端圈为多层层压式结构,具有两种结构属性,将其按结构特点纵向划分为几部分,采用Voigt模型和Reuss模型组合计算;类垫块层采用Voigt模型等效,而各层之间采用Reuss模型等效。
式中:Eφ为材料的等效弹性模量;E1为两相材料中相1即绝缘纸板的弹性模量;E2为两相材料中相2即空隙的弹性模量;φ1和φ2分别为复合材料中相1和相2的体积分数,且φ1+φ2=1。
3.在ANSYS中定义并分配各结构的材料属性:电磁特性和结构特性,包括如相对磁导率、密度、弹性模量、泊松比;等效结构部分要赋予之前计算的等效参数:等效密度、各绝缘层的等效弹性模量。
4.根据模型不同部位精细化程度的不同,分区域合理划分网格结构。模型整体网格尺寸设置为0.02,绕组及周围增大剖分网格数量,尺寸设置为0.015,整体采用四边形剖分,并控制网格形状,避免出现畸变。
5.考虑重力的影响,定义轴向坐标系下的重力加速度g=9.8m/s2;考虑系统阻尼的影响,根据具体变压器阻尼大小,定义系统瑞利阻尼常数α和β。
6.设置模型边界条件。轴对称边界施加磁通平行条件;固定件施加位移约束条件,分别对压板上边界和托板下边界施加轴向位移约束。
7.分别以各线饼的等效电流密度作为激励,以时间历程定义载荷步。设置故障持续0.2s的一段动态载荷,按照1ms为步长,分别对高、低压绕组各线饼加载。
8.采用瞬态分析方法,打开大形变效应和瞬态效应,以0.1ms划分载荷子步,每隔5个子步保存一次计算结果。
9.利用ANSYS后处理器,输出故障过程中绕组振动的位移、速度和加速度随时间变化曲线。
以上所述仅是对本发明的较佳实施方式而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施方式所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。
Claims (6)
1.一种变压器绕组短路振动特性的计算方法,其特征在于按照以下步骤进行:
步骤1:根据变压器尺寸和部件,在ANSYS中建立变压器二维有限元模型;
步骤2:由于二维模型的等效性,分别计算变压器各部分结构的等效属性;
(1)计算绕组线饼、垫块层和端圈等结构的等效密度,分别均匀等效处理;
(2)计算不同位置绝缘纸板在轴向预紧力作用下的弹性模量;
(3)计算各绝缘层的等效弹性模量;
步骤3:在ANSYS中定义并分配各结构的材料属性:电磁特性和结构特性;等效结构要赋予步骤2中计算的等效参数:等效密度、等效弹性模量;
步骤4:根据模型不同部位精细化程度的不同,分区域合理划分网格结构,并控制网格形状,避免出现畸变;
步骤5:定义轴向坐标系下的重力加速度g=9.8m/s2;考虑系统阻尼的影响,根据具体变压器阻尼大小,定义系统瑞利阻尼常数α和β;
步骤6:设置模型边界条件,轴对称边界施加磁通平行条件,固定件施加位移约束条件,分别对压板上边界和托板下边界施加轴向位移约束;
步骤7:分别以各线饼的等效电流密度作为激励,以时间历程定义载荷步;设置故障持续0.2s的一段动态载荷,按照1ms为步长,分别对高、低压绕组各线饼加载;
步骤8:采用瞬态分析方法,打开大形变效应和瞬态效应,划分载荷子步,设置写入结果文件的频率;
步骤9:利用ANSYS后处理器,输出故障过程中绕组振动的位移、速度和加速度随时间变化曲线。
2.按照权利要求1所述一种变压器绕组短路振动特性的计算方法,其特征在于:所述步骤1中ANSYS中建立变压器二维有限元模型是指在ANSYSMechanical APDL软件环境下,选择PLANE13磁-结构耦合有限元单元,设置单元具有UX、UY、AZ自由度,直接耦合磁场和结构场,选取二维轴对称建模方式,根据变压器实际尺寸,建立二维轴对称模型。
3.按照权利要求1所述一种变压器绕组短路振动特性的计算方法,其特征在于:所述步骤2中由于垫块、端圈、压板和压木等结构均采用绝缘纸板制成,在一定范围内,其应力-应变关系可表示为:
σ=aε+bε3 (1)
式中:σ为应力;ε为应变;a为线性常数;b为硬化系数。
通过实验可测得:a=105MPa;b=1750MPa。
在轴向预紧力作用下,其弹性模量可表示为:
压板和托板的弹性模量采用此参数直接定义。
考虑各相同性的两种材料,在载荷作用下产生相同应变,建立了Voigt模型:
Eφ=φ1E1+φ2E2 (3)
垫块层和端部压木层采用Voigt模型等效。
考虑两相承受相同的应力,建立了Reuss模型:
端圈为多层层压式结构,具有两种结构属性,将其按结构特点纵向划分为几部分,采用Voigt模型和Reuss模型组合计算。
式中:Eφ为材料的等效弹性模量;E1为绝缘纸板的弹性模量;E2为空隙的弹性模量;φ1和φ2分别为绝缘纸板和空隙的体积分数,且φ1+φ2=1。
4.按照权利要求1所述一种变压器绕组短路振动特性的计算方法,其特征在于:所述步骤4中模型整体网格尺寸设置为0.02,绕组及周围增大剖分网格数量,尺寸设置为0.015,整体采用四边形剖分。
5.按照权利要求1所述一种变压器绕组短路振动特性的计算方法,其特征在于:所述步骤8中以0.1ms划分载荷子步,每隔5个子步保存一次计算结果。
6.按照权利要求1所述一种变压器绕组短路振动特性的计算方法,其特征在于:利用ANSYS后处理器,输出故障过程中绕组振动的位移、速度和加速度随时间变化曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710030273.4A CN106777812B (zh) | 2017-01-17 | 2017-01-17 | 一种变压器绕组短路振动特性的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710030273.4A CN106777812B (zh) | 2017-01-17 | 2017-01-17 | 一种变压器绕组短路振动特性的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106777812A true CN106777812A (zh) | 2017-05-31 |
CN106777812B CN106777812B (zh) | 2020-06-09 |
Family
ID=58946274
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710030273.4A Expired - Fee Related CN106777812B (zh) | 2017-01-17 | 2017-01-17 | 一种变压器绕组短路振动特性的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106777812B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112364443A (zh) * | 2020-09-03 | 2021-02-12 | 国网天津市电力公司电力科学研究院 | 一种变压器绕组轴向固有振动特性定量分析方法 |
CN112541289A (zh) * | 2020-12-15 | 2021-03-23 | 西南交通大学 | 一种轨道车辆车体铝型材结构的声振等效建模方法 |
CN113806974A (zh) * | 2021-08-06 | 2021-12-17 | 南方电网科学研究院有限责任公司 | 一种变压器套管电连接的稳定性评估方法、装置及系统 |
CN116933565A (zh) * | 2023-09-12 | 2023-10-24 | 合肥工业大学 | 考虑磁-结构耦合影响的绕组轴向短路电磁力分析方法 |
CN117540539A (zh) * | 2023-10-25 | 2024-02-09 | 湖南大学 | 一种基于振动加速度的变压器优化设计方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011090456A (ja) * | 2009-10-21 | 2011-05-06 | Jfe Steel Corp | トランス用方向性電磁鋼板の受注・出荷方法 |
CN105095609A (zh) * | 2015-09-21 | 2015-11-25 | 武汉大学 | 基于有限元法的变压器电磁振动噪声计算方法 |
CN105224724A (zh) * | 2015-09-10 | 2016-01-06 | 云南电网有限责任公司电力科学研究院 | 一种变压器绕组短路冲击位移与振动加速度的计算方法 |
-
2017
- 2017-01-17 CN CN201710030273.4A patent/CN106777812B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011090456A (ja) * | 2009-10-21 | 2011-05-06 | Jfe Steel Corp | トランス用方向性電磁鋼板の受注・出荷方法 |
CN105224724A (zh) * | 2015-09-10 | 2016-01-06 | 云南电网有限责任公司电力科学研究院 | 一种变压器绕组短路冲击位移与振动加速度的计算方法 |
CN105095609A (zh) * | 2015-09-21 | 2015-11-25 | 武汉大学 | 基于有限元法的变压器电磁振动噪声计算方法 |
Non-Patent Citations (2)
Title |
---|
王勇: "场路结合并考虑耦合的磁力机械分析与设计方法研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
田云德: "复合材料等效弹性模量的改进混合律方法", 《西南交通大学学报》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112364443A (zh) * | 2020-09-03 | 2021-02-12 | 国网天津市电力公司电力科学研究院 | 一种变压器绕组轴向固有振动特性定量分析方法 |
CN112364443B (zh) * | 2020-09-03 | 2022-05-24 | 国网天津市电力公司电力科学研究院 | 一种变压器绕组轴向固有振动特性定量分析方法 |
CN112541289A (zh) * | 2020-12-15 | 2021-03-23 | 西南交通大学 | 一种轨道车辆车体铝型材结构的声振等效建模方法 |
CN112541289B (zh) * | 2020-12-15 | 2023-04-07 | 西南交通大学 | 一种轨道车辆车体铝型材结构的声振等效建模方法 |
CN113806974A (zh) * | 2021-08-06 | 2021-12-17 | 南方电网科学研究院有限责任公司 | 一种变压器套管电连接的稳定性评估方法、装置及系统 |
CN113806974B (zh) * | 2021-08-06 | 2023-10-17 | 南方电网科学研究院有限责任公司 | 一种变压器套管电连接的稳定性评估方法、装置及系统 |
CN116933565A (zh) * | 2023-09-12 | 2023-10-24 | 合肥工业大学 | 考虑磁-结构耦合影响的绕组轴向短路电磁力分析方法 |
CN116933565B (zh) * | 2023-09-12 | 2023-12-01 | 合肥工业大学 | 考虑磁-结构耦合影响的绕组轴向短路电磁力分析方法 |
CN117540539A (zh) * | 2023-10-25 | 2024-02-09 | 湖南大学 | 一种基于振动加速度的变压器优化设计方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106777812B (zh) | 2020-06-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106777812A (zh) | 一种变压器绕组短路振动特性的计算方法 | |
Hirane et al. | On the layerwise finite element formulation for static and free vibration analysis of functionally graded sandwich plates | |
Li et al. | Vibration and damping study of multifunctional grille composite sandwich plates with an IMAS design approach | |
CN103559418B (zh) | 一种基于故障物理的电子产品任务可靠度计算方法 | |
Dung et al. | The Third‐Order Shear Deformation Theory for Modeling the Static Bending and Dynamic Responses of Piezoelectric Bidirectional Functionally Graded Plates | |
Carrera et al. | Component-wise method applied to vibration of wing structures | |
Noori et al. | Dynamic analysis of functionally graded porous beams using complementary functions method in the Laplace domain | |
Chen et al. | Vibration analysis of functionally graded Timoshenko beams | |
CN106768559B (zh) | 一种短路故障下变压器绕组的动态受力分析方法 | |
Chen et al. | Isogeometric three-dimensional vibration of variable thickness parallelogram plates with in-plane functionally graded porous materials | |
Gao et al. | Robustly tuning bandgaps in two-dimensional soft phononic crystals with criss-crossed elliptical holes | |
CN106650007A (zh) | 变压器稳态时漏磁场及绕组振动位移计算方法 | |
Carrera et al. | On the role of large cross-sectional deformations in the nonlinear analysis of composite thin-walled structures | |
Tran et al. | A nonlocal IGA numerical solution for free vibration and buckling analysis of porous sigmoid functionally graded (P-SFGM) nanoplate | |
Wei et al. | Research on dynamic model of printed circuit board based on finite element method | |
CN107862130B (zh) | 一种导致继电器贮存退化的内部关键因素的分析方法 | |
Zhu et al. | Fast electrothermal coupling calculation method for supporting digital twin construction of electrical equipment | |
Shahmohammadi et al. | Buckling of multilayered CNT/GPL/fibre/polymer hybrid composite plates resting on elastic support using modified nonlocal first-order plate theory | |
Nguyen et al. | A new finite element procedure for the dynamic analysis of BDFGS plates located on pasternak foundation subjected to the moving oscillator load | |
CN106682284A (zh) | 复合材料构件热校形工艺的仿真模拟方法 | |
Sun et al. | Warpage simulation and DOE analysis with application in package-on-package development | |
Le et al. | Vibration Analysis of Functionally Graded Beams Using a Novel Two-Variable Model with Various Boundary Conditions | |
Pei et al. | Prediction of wire sweep during the encapsulation of IC packaging with wire density effect | |
Ranjan et al. | Thermo-elastic free vibration analysis of functionally graded flat panel with temperature gradient along thickness | |
CN105138746A (zh) | 通过扬声器振动部件的劲度系数反推其杨氏模量的方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200609 Termination date: 20210117 |