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
- equivalent
- model
- winding
- calculation method
- 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 20
- 238000009413 insulation Methods 0.000 claims abstract description 11
- 230000000694 effects Effects 0.000 claims abstract description 9
- 238000006073 displacement reaction Methods 0.000 claims abstract description 8
- 238000004458 analytical method Methods 0.000 claims abstract description 7
- 238000010276 construction Methods 0.000 claims abstract description 5
- 238000000034 method Methods 0.000 claims description 10
- 238000013461 design Methods 0.000 claims description 5
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000005516 engineering process Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000000007 visual effect Effects 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 |
|---|---|---|
| CN112287491B (zh) | 一种复合型点阵材料及其设计方法 | |
| CN102096748B (zh) | 一种轻量化大刚度精冲压力机机身的优化设计方法 | |
| CN106768559B (zh) | 一种短路故障下变压器绕组的动态受力分析方法 | |
| CN110555263B (zh) | 一种用于曲壳结构优化设计的水平集拓扑优化方法 | |
| CN108629140B (zh) | 一种基于测地线距离的带孔复合材料结构设计优化方法 | |
| CN106777812A (zh) | 一种变压器绕组短路振动特性的计算方法 | |
| CN101567021A (zh) | 用于振动发电的矩形悬臂梁压电振子优化设计有限元方法 | |
| Xing et al. | An eigenelement method of periodical composite structures | |
| CN105183988B (zh) | 一种接地螺栓通电后温度及应力应变有限元计算分析方法 | |
| CN112446163B (zh) | 基于参数化水平集的能量有限元拓扑优化方法 | |
| CN107194064A (zh) | 一种变压器绕组固有频率计算方法 | |
| CN111079327A (zh) | 一种基于各向异性过滤技术的异形曲面加筋拓扑优化方法 | |
| CN115577603A (zh) | 降低单元矩阵维度的仿真方法、系统及相关设备 | |
| CN108170944A (zh) | 一种半导体器件的压力均衡制作参数优化方法及制作方法 | |
| Wei et al. | Research on dynamic model of printed circuit board based on finite element method | |
| CN105956347B (zh) | 纤维织物增强陶瓷基复合材料力学行为的仿真预测方法 | |
| Zhou et al. | A ready-to-manufacture optimization design of 3d chiral auxetics for additive manufacturing | |
| CN113033033A (zh) | 一种波纹管结构优化方法及波纹管 | |
| Vukovic et al. | Simulation platform for flexible electronics | |
| CN117807833B (zh) | 一种压电俘能器的拓扑结构优化方法 | |
| CN117787208B (zh) | 印刷电路板变形仿真方法、装置、电子设备和存储介质 | |
| CN108491595A (zh) | 一种声-固耦合结构的高频局部响应预示方法 | |
| Sladek et al. | Plate bending problems using the nonsingular boundary element formulation and C1-continuous elements | |
| CN111883216A (zh) | 一种基于材料属性优化的仿真网格划分方法及系统 | |
| CN117669314A (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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170922 |