CN107194064A - 一种变压器绕组固有频率计算方法 - Google Patents
一种变压器绕组固有频率计算方法 Download PDFInfo
- Publication number
- CN107194064A CN107194064A CN201710358122.1A CN201710358122A CN107194064A CN 107194064 A CN107194064 A CN 107194064A CN 201710358122 A CN201710358122 A CN 201710358122A CN 107194064 A CN107194064 A CN 107194064A
- Authority
- CN
- China
- Prior art keywords
- transformer winding
- model
- equivalent
- winding
- natural frequency
- 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
Links
- 238000004804 winding Methods 0.000 title claims abstract description 50
- 238000004364 calculation method Methods 0.000 title claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 22
- 238000009413 insulation Methods 0.000 claims abstract description 11
- 239000000463 material Substances 0.000 claims abstract description 9
- 125000006850 spacer group Chemical group 0.000 claims abstract description 9
- 238000006073 displacement reaction Methods 0.000 claims abstract description 8
- 230000000694 effects Effects 0.000 claims abstract description 6
- 239000011111 cardboard Substances 0.000 claims description 9
- 230000036316 preload Effects 0.000 claims description 8
- 241001421775 Thereus Species 0.000 claims description 7
- 239000011800 void material Substances 0.000 claims description 6
- 239000002131 composite material Substances 0.000 claims description 2
- 239000011087 paperboard Substances 0.000 claims description 2
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 1
- 239000000123 paper Substances 0.000 description 1
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)
- Housings And Mounting Of Transformers (AREA)
Abstract
本发明公开了一种变压器绕组固有频率计算方法,在ANSYS中建立变压器绕组的二维轴对称模型;计算变压器绕组结构的等效密度、等效弹性模量,在ANSYS中定义并分配材料属性参数。划分网格结构,控制网格形状,并避免出现畸变。设置绕组位移约束条件,分别对压板的上边界和托板的下边界施加全位移约束。对模型进行模态分析,输出计算结果。修改垫块层、端圈和端绝缘的等效弹性模量,等效不同预紧力的作用效果,计算绕组各阶固有频率。本发明的有益效果是可实现快速等效计算不同预紧力下的各阶固有频率。
Description
技术领域
本发明属于变压器技术领域,涉及一种变压器绕组固有频率计算方法。
背景技术
电力变压器是电力系统中的重要设备之一。电力变压器因外部短路而损坏原因众多,情况复杂。从近年来统计显示,变压器本身抗短路能力不足,是造成其损坏主要原因。其中,由于绕组轴向预紧力不足,绕组在轴向力作用下出现结构松动,绕组变形的事故也呈上升趋势。
电力变压器绕组固有频率是变压器设计的一个重要技术指标。如果变压器绕组固有频率和短路力频率相近,则易发生谐振。振动可能造成变压器绕组松动,甚至严重破坏。所以,变压器设计时绕组固有频率应避开短路力频率,避免发生短路谐振。变压器绕组固有频率目前采用三维有限元法计算,但是三维有限元法计算变压器绕组固有频率过程复杂,运算量大。
发明内容
本发明的目的在于提供一种变压器绕组固有频率计算方法,解决了目前三维有限元法计算变压器绕组固有频率过程复杂的问题。
本发明所采用的技术方案是按照以下步骤进行:
步骤1:在ANSYS中建立变压器绕组的二维轴对称模型;
步骤2:计算变压器绕组结构的等效密度、等效弹性模量,在ANSYS中定义并分配材料属性参数。
步骤3:划分网格结构,控制网格形状,并避免出现畸变。
步骤4:设置绕组位移约束条件,分别对压板的上边界和托板的下边界施加全位移约束。
步骤5:对模型进行模态分析,输出计算结果。
步骤6:修改垫块层、端圈和端绝缘的等效弹性模量,等效不同预紧力的作用效果,计算绕组各阶固有频率。
进一步,步骤1中ANSYS中建立变压器二维有限元模型是指在ANSYS MechanicalAPDL软件环境下,选择PLANE13有限元单元,设置单元具有UX、UY自由度,选取二维轴对称建模方式,根据变压器绕组结构尺寸,建立二维轴对称模型。
进一步,步骤2中的变压器绕组各结构包括绕组线饼、垫块层、端圈、端绝缘。
进一步,步骤2中的垫块、端圈和端绝缘假设均采用绝缘纸板制成,应力-应变关系表示为:
σ=aε+bε3 (1)
式中:σ为应力;ε为应变;a为线性常数;b为硬化系数;
其中a=105MPa;b=1750MPa;
弹性模量表示为:
考虑轴向预紧力作用影响,根据应力-应变关系,得不同预紧力下绝缘纸
板的弹性模量;
考虑两相在载荷作用下产生相同的应变,建立了Voigt模型:
Eφ=φ1E1+φ2E2 (3)
考虑两相承受相同的应力,建立了Reuss模型:
式中:Eφ为材料的等效弹性模量;E1为绝缘纸板的弹性模量;E2为空隙的弹性模量;φ1和φ2分别为绝缘纸板和空隙的体积分数,且φ1+φ2=1;垫块层采用Voigt模型等效;端圈采用多层层压式结构,具有两种结构属性,将其按结构特点纵向划分为几部分,采用Voigt模型和Reuss模型的组合计算。
进一步,步骤3中模型整体网格尺寸设置为0.02,采用四边形剖分。
进一步,步骤5中模态分析采用BlockLanczos法,计算前20阶固有频率。
本发明的有益效果是可实现快速等效计算不同预紧力下的各阶固有频率。
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
本发明变压器绕组固有频率计算方法的计算步骤如下:
步骤1:在ANSYS中建立变压器绕组的二维轴对称模型。在ANSYS Mechanical APDL软件环境下,选择PLANE13有限元单元,设置单元具有UX、UY自由度,选取二维轴对称建模方式,根据变压器绕组结构尺寸,建立二维轴对称模型。
步骤2:计算变压器绕组结构的等效材料属性。
(1)计算绕组线饼、垫块层、端圈和端绝缘的等效密度;
(2)计算在轴向预紧力作用下垫块层、端圈和端绝缘的等效弹性模量。
在ANSYS中定义并分配各结构的材料属性:等效密度、等效弹性模量、泊松比。
由于垫块、端圈和端绝缘等结构,均采用绝缘纸板制成,在一定范围内,其应力-应变关系可表示为:
σ=aε+bε3 (1)
式中:σ为应力;ε为应变;a为线性常数;b为硬化系数。
通过实验可测得:a=105MPa;b=1750MPa。
其弹性模量可表示为:
考虑轴向预紧力作用影响,根据应力-应变关系,可得不同预紧力下绝缘纸板的弹性模量。
考虑各相同性的两种材料的结合方式不同,分别建立了Voigt模型和Reuss模型。
考虑两相在载荷作用下产生相同的应变,建立了Voigt模型:
Eφ=φ1E1+φ2E2 (3)
考虑两相承受相同的应力,建立了Reuss模型:
式中:Eφ为材料的等效弹性模量;E1为绝缘纸板的弹性模量;E2为空隙的弹性模量;φ1和φ2分别为绝缘纸板和空隙的体积分数,且φ1+φ2=1。垫块层可以采用Voigt模型等效。端圈采用多层层压式结构,具有两种结构属性,将其按结构特点纵向划分为几部分,采用Voigt模型和Reuss模型的组合计算。
步骤3:划分网格结构,控制网格形状,并避免出现畸变。本发明模型整体网格尺寸设置为0.02,采用四边形剖分。
步骤4:设置绕组位移约束条件,分别对压板的上边界和托板的下边界施加全位移约束。
步骤5:对模型进行模态分析,输出计算结果。模态分析采用BlockLanczos法,计算前20阶固有频率。
步骤6:修改垫块层、端圈和端绝缘的等效弹性模量,等效不同预紧力的作用效果,计算绕组各阶固有频率。
本发明公开了一种变压器绕组固有频率计算方法,主要包括以下步骤:用ANSYS建立变压器绕组二维有限元模型;根据绕组结构和轴向预紧力大小,赋予其等效参数;采用模态分析法,计算其各阶固有频率;通过修改绕组等效弹性模量参数,考虑其不同预紧力作用下的各阶固有频率。本方法实现了变压器绕组的二维动力学等效,及快速计算不同预紧力下的各阶固有频率。
以上所述仅是对本发明的较佳实施方式而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施方式所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。
Claims (6)
1.一种变压器绕组固有频率计算方法,其特征在于按照以下步骤进行:
步骤1:在ANSYS中建立变压器绕组的二维轴对称模型;
步骤2:计算变压器绕组结构的等效密度、等效弹性模量,在ANSYS中定义并分配材料属性参数。
步骤3:划分网格结构,控制网格形状,并避免出现畸变。
步骤4:设置绕组位移约束条件,分别对压板的上边界和托板的下边界施加全位移约束。
步骤5:对模型进行模态分析,输出计算结果。
步骤6:修改垫块层、端圈和端绝缘的等效弹性模量,等效不同预紧力的作用效果,计算绕组各阶固有频率。
2.按照权利要求1所述一种变压器绕组固有频率计算方法,其特征在于:所述步骤1中ANSYS中建立变压器二维有限元模型是指在ANSYS Mechanical APDL软件环境下,选择PLANE13有限元单元,设置单元具有UX、UY自由度,选取二维轴对称建模方式,根据变压器绕组结构尺寸,建立二维轴对称模型。
3.按照权利要求1所述一种变压器绕组固有频率计算方法,其特征在于:所述步骤2中的变压器绕组各结构包括绕组线饼、垫块层、端圈、端绝缘。
4.按照权利要求1所述一种变压器绕组固有频率计算方法,其特征在于:所述步骤2中的垫块、端圈和端绝缘假设均采用绝缘纸板制成,应力-应变关系表示为:
σ=aε+bε3 (1)
式中:σ为应力;ε为应变;a为线性常数;b为硬化系数;
其中a=105MPa;b=1750MPa;
弹性模量表示为:
考虑轴向预紧力作用影响,根据应力-应变关系,得不同预紧力下绝缘纸板的弹性模量;
考虑两相在载荷作用下产生相同的应变,建立了Voigt模型:
Eφ=φ1E1+φ2E2 (3)
考虑两相承受相同的应力,建立了Reuss模型:
式中:Eφ为材料的等效弹性模量;E1为绝缘纸板的弹性模量;E2为空隙的弹性模量;φ1和φ2分别为绝缘纸板和空隙的体积分数,且φ1+φ2=1;垫块层采用Voigt模型等效;端圈采用多层层压式结构,具有两种结构属性,将其按结构特点纵向划分为几部分,采用Voigt模型和Reuss模型的组合计算。
5.按照权利要求1所述一种变压器绕组固有频率计算方法,其特征在于:所述步骤3中模型整体网格尺寸设置为0.02,采用四边形剖分。
6.按照权利要求1所述一种变压器绕组固有频率计算方法,其特征在于:所述步骤5中模态分析采用Block Lanczos法,计算前20阶固有频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710358122.1A CN107194064A (zh) | 2017-05-19 | 2017-05-19 | 一种变压器绕组固有频率计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710358122.1A CN107194064A (zh) | 2017-05-19 | 2017-05-19 | 一种变压器绕组固有频率计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107194064A true CN107194064A (zh) | 2017-09-22 |
Family
ID=59874172
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710358122.1A Pending CN107194064A (zh) | 2017-05-19 | 2017-05-19 | 一种变压器绕组固有频率计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107194064A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110554278A (zh) * | 2019-07-30 | 2019-12-10 | 中国电力科学研究院有限公司 | 一种基于模态分析对变压器绕组进行检测的方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120063039A1 (en) * | 2010-09-15 | 2012-03-15 | Shah Saurabh B | Systems and methods for protection of components in electrical power delivery systems |
CN103267907A (zh) * | 2013-04-19 | 2013-08-28 | 上海交通大学 | 一种变压器绕组模态参数识别方法 |
CN105608260A (zh) * | 2015-12-17 | 2016-05-25 | 上海电器科学研究院 | 一种交流电动机有绕组定子模态精确计算方法 |
-
2017
- 2017-05-19 CN CN201710358122.1A patent/CN107194064A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20120063039A1 (en) * | 2010-09-15 | 2012-03-15 | Shah Saurabh B | Systems and methods for protection of components in electrical power delivery systems |
CN103267907A (zh) * | 2013-04-19 | 2013-08-28 | 上海交通大学 | 一种变压器绕组模态参数识别方法 |
CN105608260A (zh) * | 2015-12-17 | 2016-05-25 | 上海电器科学研究院 | 一种交流电动机有绕组定子模态精确计算方法 |
Non-Patent Citations (5)
Title |
---|
YU YANG 等: "Broadband frequency response analysis of transformer windings", 《IEEE TRANSACTIONS ON DIELECTRICS AND ELECTRICAL INSULATION》 * |
于怀金: "大型变压器漏磁场及短路电动力的计算与分析", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
李凯琪: "电力变压器抗短路能力核算模型与评估技术研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
田云德 等: "复合材料等效弹性模量的改进混合律方法", 《西南交通大学学报》 * |
罗惠引: "基于光纤法珀传感器的变压器铁心振动信号测试研究", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110554278A (zh) * | 2019-07-30 | 2019-12-10 | 中国电力科学研究院有限公司 | 一种基于模态分析对变压器绕组进行检测的方法及系统 |
CN110554278B (zh) * | 2019-07-30 | 2023-05-19 | 中国电力科学研究院有限公司 | 一种基于模态分析对变压器绕组进行检测的方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106777812B (zh) | 一种变压器绕组短路振动特性的计算方法 | |
CN109413974B (zh) | 一种多层结构吸波材料及其制备方法 | |
CN109033626A (zh) | 一种特高压并联电抗器振动噪声计算方法 | |
CN203397878U (zh) | 一种变压器式可控电抗器磁集成装置 | |
CN203968362U (zh) | 扬声器和音频播放设备 | |
CN112287583A (zh) | 一种变压器的噪声仿真预测方法 | |
CN105118647A (zh) | 大容量高频变压器最优工作频率的确定方法 | |
CN107194064A (zh) | 一种变压器绕组固有频率计算方法 | |
CN203277043U (zh) | 一种降低油浸式非晶合金变压器噪音结构 | |
CN206780796U (zh) | 一种可调节电磁碳纤维复合材料成型装置 | |
Cai et al. | Thermal-mechanical-electrical coupled design of multilayer energy storage ceramic capacitors | |
CN203433957U (zh) | 一种低噪声sf6气体绝缘变压器 | |
CN206163266U (zh) | 一种新型高频三相电感 | |
CN105335575B (zh) | 一种电子天线结构一体化层压板的有限元建模方法 | |
CN115562526A (zh) | 一种变电站的参数化建模预处理方法、系统及电子设备 | |
CN204289037U (zh) | 一种铁心电抗器的铁心 | |
CN114779030A (zh) | 一种确定变压器出线装置的绝缘强度的方法及系统 | |
Li et al. | Axial stability analysis and simulation of power transformer windings | |
CN204257369U (zh) | 一种干式电力变压器 | |
Tong et al. | Electromagnetic-thermal coupled analysis of a 10 kV silicone rubber cast dry-type transformer | |
CN203871281U (zh) | 一种高压断路器绝缘支架 | |
Chen et al. | Research on Electrical Characteristics of Metal Rubber | |
CN204117804U (zh) | 一种电感 | |
CN103440673A (zh) | 基于元胞自动机的刻蚀表面演化模型压缩方法 | |
CN203673844U (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170922 |
|
RJ01 | Rejection of invention patent application after publication |