一种基于变压器三维电磁耦合模型的绕组应力评估方法
一种基于变压器三维电磁耦合模型的绕组应力评估方法
【技术领域】
[0001] 本发明是一种基于变压器Ξ维电磁禪合模型的绕组应力评估方法，应用于电力变 压器绕组受力分析与结构稳定性评估。
【背景技术】
[0002] 电力变压器是电网和电力传输中最重要的设备之一，其绕组承受短路的能力直接 影响电网能否安全可靠的运行。在单次短路冲击作用下，绕组的稳定性取决于整定值与实 际受力值，当整定值大于实际受力值时，绕组是稳定的。在实际运行中变压器可能受到多次 短路冲击，绕组的抗短路可能会随着短路冲击次数的增加而降低，同时导线的应力可能会 增大，表现出一定的累积效应。由于在每次短路冲击作用下，绕组的结构形状都会发生形 变，绕组的不同形状又会使得整个磁场分布发生变化，进而导致绕组所受到的力发生变化。 根据结构力学可知，不同形状下绕组发生失稳的条件也将变化。因此研究不同累积情况下 绕组所受到的应力情况，运对电力变压器绕组结构稳定性评估有一定的指导意义和实际价 值。
【发明内容】
[0003] 本发明的目的是，提供一种科学合理，真实有效，实用价值高的基于变压器Ξ维电 磁禪合模型的绕组应力评估方法，能够根据实际变压器参数，建立变压器Ξ维电磁禪合模 型并求解，从与Ξ维模型对应的二维模型中计算得到绕组受力值，然后将受力值加载到对 应的Ξ维模型中计算得到绕组受到的应力。
[0004] 本发明的目的是由W下技术方案来实现的：一种基于变压器Ξ维电磁禪合模型的 绕组应力评估方法，其特征是，它包括W下步骤：
[0005] 1)变压器漏磁场的有限元计算
[0006] 基于矢量磁位节点有限元法对变压器漏磁场进行计算，矢量磁位节点有限元法采 用矢量磁位A，忽略磁滞效应和满流效应，根据Maxwel 1方程得到变压器非线性磁场方程：
[0007]
(1)
[000引公式（1)中，μ为媒质的磁导率;J为电流密度，需要通过电路等效参数计算获得； "▽X"代表旋度运算，
[0009] 矢量磁位节点有限元法中节点单元的自由度为节点i的矢量磁位Ai，采用标量形 状函数Ni，A为未知变量，则单元插值函数：
[0010]
(2)
[0011]公式（2)中，nnode为单元节点数，
[0012] 求解场域的整体插值函数：
[0013]
(3)
[0014] 公式(3)中，Μη(η=1，2,…，nn)为基函数序列，由相关单元形状函数Ni对应叠加形 成，η为基函数序列通项编号，nn为总项数，即总节点数，An为单元标量磁位；
[0015] 对公式(2)应用格林定理，得伽迂金加权余量方程：
[0016]
(4)
[0017] 公式(4)中，Mm(m=l，2,…，Πη)为权函数序列，m为权函数序列通项编号，V表示边界 体积分量，S表示边界面积分量，边界面积分项为边界面磁场强度的切向分量的贡献，en为 边界面单位法向分量，在变压器计算模型中只设及该项为零的情况，所W公式(4)整理得： [001 引
[0019] 将公式(3)代入公式(5)，针对全部权函数，可W将加权余量方程离散形成代数方 程组，通过求解可得所有节点上的矢量磁位Α，利用磁场的微分公式B = VXA，则区域内任 意位置的磁通密度都可W表示：
[0020] 基于矢量磁位节点有限元方法计算漏磁场，可得变压器中漏磁场的最大漏磁场强 度公式，
[0021]
(6)
[0022] 公式(6)中，I表示电流有效值，Bm为漏磁场强度，Hm为漏磁感应强度，Z为线圈应数， P为罗氏系数，μ〇为真空磁导率，Η为线圈高度；
[0023] 2)变压器绕组应力的计算
[0024] 变压器绕组应力计算主要包含，绕组的平均环行拉应力、压缩应力和撑条或垫块 之间的导线福向弯曲应力，
[0025] 绕组受到福向电磁力时，在直径为Dn的圆周上的压力：
[0026]
饼
[0027] 公式(7)中，Fx为两垫块间导线上的福向力，Dn为圆周的平均直径，出为线饼高度，P 为单位面积上的压力；
[0028] 绕组导线材料中产生的应力，
[0029]
(s)
[0030] 公式(8)中，β为线饼的福向宽度;οχ为导线福向应力，
[0031] 导线上的最大弯矩：
[0032]
側
[0033] 公式(9)中，Fy为两垫块间导线上的轴向力，S为两垫块间导线上的长度，
[0034] 导线为矩形截面的抗弯截面系数：
[0035]
(10)
[0036] 公式(10)中，d为单根导线福向厚度;b为单根导线轴向高度，
[0037] 幅向垫块之间跨度内导线轴向弯曲应力：
[00；3 引
(Π )
[0039] 公式(11)中，OY为导线轴向应力。
[0040] 本发明的基于变压器Ξ维电磁禪合模型的绕组应力评估方法，能够根据实际变压 器参数，建立变压器Ξ维电磁禪合模型并求解，W矢量磁位节点有限元法为基础计算变压 器漏磁场，通过模拟变压器内部磁场分布，对绕组所受到的应力进行计算分析。由于绕组发 生形变后，使得绕组的几何结构不再规则，导致利用加载电流原理而编制的磁结构禪合模 型不再适用。所W从与Ξ维模型对应的二维模型中计算绕组受力值，然后将受力值通过相 关的程序加载到对应的Ξ维模型中来运算绕组所受应力。为了更好地模拟绕组实际情况， 在绕组模型中设定了撑条，并且通过两撑条间绕组凹陷程度的不同来体现绕组的形变，进 而观察不同累积情况下，模型所受到的应力情况。能够为变压器绕组结构稳定性的判定提 供依据，具有科学合理，真实有效，实用价值高等优点。
【附图说明】
[0041] 图1是未发生形变的变压器绕组模型立体结构示意图；
[0042] 图2是图1的左视示意图；
[0043] 图3是发生形变的变压器绕组模型立体结构示意图；
[0044] 图4是图3的左视示意图；
[0045] 图5是变压器内绕组累积效应最大应力示意图；
[0046] 图6是变压器内绕组累积效应最大应力细化示意图；
[0047] 图7是不同型号变压器内绕组累积效应最大应力示意图；
[004引图8是变压器外绕组累积效应最大应力示意图。
【具体实施方式】
[0049] 下面结合附图和具体实施例对本发明作进一步描述：
[0050] 本发明的一种基于变压器Ξ维电磁禪合模型的绕组应力评估方法，包括W下步 骤：
[0051 ] 1)变压器漏磁场的有限元计算
[0052]基于矢量磁位节点有限元法对变压器漏磁场进行计算，矢量磁位节点有限元法采 用矢量磁位A，忽略磁滞效应和满流效应，根据Maxwel 1方程得到变压器非线性磁场方程： [0化3]
(1)
[0054] 公式（1)中，μ为媒质的磁导率;J为电流密度，需要通过电路等效参数计算获得； "▽X"代表旋度运算，
[0055] 矢量磁位节点有限元法中节点单元的自由度为节点i的矢量磁位Ai，采用标量形 状函数Ni，A为未知变量，则单元插值函数：
[0化6]
樹
[0057]公式（2)中，rinode为单元节点数，
[005引求解场域的整体插值函数：
[0化9]
(3)
[0060」公式(3)中，Μη(η=1，2,…，nn)为基函数序列，由相关单元形状函数Ni对应叠加形 成，η为基函数序列通项编号，Πη为总项数，即总节点数，An为单元标量磁位；
[0061] 对公式(2)应用格林定理，得伽迂金加权余量方程：
[0062]
(4)
[0063] 公式(4)中，Mm(m= 1,2,…，Πη)为权函数序列，m为权函数序列通项编号，V表示边界 体积分量，S表示边界面积分量，边界面积分项为边界面磁场强度的切向分量的贡献，en为 边界面单位法向分量，在变压器计算模型中只设及该项为零的情况，所W公式(4)整理得：
[0064]
(勾
[0065] 将公式(3)代入公式(5)，针对全部权函数，可W将加权余量方程离散形成代数方 程组，通过求解可得所有节点上的矢量磁位A，利用磁场的微 分公式B = VXA，则区域内任 意位置的磁通密度都可W表示：
[0066] 基于矢量磁位节点有限元方法计算漏磁场，可得变压器中漏磁场的最大漏磁场强 度公式，
[0067]
批）
[0068] 公式(6)中，I表示电流有效值，Bm为漏磁场强度，Hm为漏磁感应强度，Z为线圈应数， P为罗氏系数，μ〇为真空磁导率，Η为线圈高度；
[0069] 2)变压器绕组应力的计算
[0070] 变压器绕组应力计算主要包含，绕组的平均环行拉应力、压缩应力和撑条或垫块 之间的导线福向弯曲应力，
[0071] 绕组受到福向电磁力时，在直径为Dn的圆周上的压力：
[0072]
(7)
[0073] 公式(7)中，Fx为两垫块间导线上的福向力，Dn为圆周的平均直径，出为线饼高度，P 为单位面积上的压力；
[0074] 绕组导线材料中产生的应力，
[0075]
做
[0076] 公式(8)中，β为线饼的福向宽度;οχ为导线福向应力，
[0077] 导线上的最大弯矩：
[007引
(9)
[0079] 公式(9)中，Fy为两垫块间导线上的轴向力，S为两垫块间导线上的长度，
[0080] 导线为矩形截面的抗弯截面系数： 「00811
(10)
[0082] 公式(10)中，d为单根导线福向厚度;b为单根导线轴向高度，
[0083] 幅向垫块之间跨度内导线轴向弯曲应力：
[0084]
(11)
[0085] 公式（11)中，OY为导线轴向应力。
[0086] 变压器内绕组累积效应作用下的最大应力分析：
[0087] 1)变压器内绕组累积效应模型的建立
[0088] 假设档距之间凹陷的最大变形为2mm，即每经过一次短路冲击作用后，绕组凹陷形 变的最大处在原来的基础上形变2mm，W此为基础得到绕组的二维模型，然后将计算得到的 受力值施加在相对应的Ξ维模型上，进而得出所需的结论。未发生形变的变压器绕组模型 立体结构图和左视图，如图1、图2所示；发生形变的变压器绕组模型立体结构图和左视图， 如图3、图4所示。
[0089] 2)变压器内绕组累积效应
[0090] 本发明对变压器内绕组的累积情况均是W2mm的形变量做研究的，并且假设变压 器内绕组经受过五次短路冲击，则变压器内绕组在不同形变量下的最大应力计算结果，如 表1所示。
[0091] 表1变压器内绕组累积效应计算结果
[0092]
[0093] 由表1得到变压器内绕组累积效应最大应力图，如图5和图6所示。
[0094] 由变压器内绕组累积效应最大应力图可知绕组所受应力在形变量为8mm附近处发 生突变。因此为了使研究结果更准确，本发明在上述形变的基础上，对拐点附近(8-lOmm)处 进行0.5mm细化;对突变点附近(9-9.6mm)处进行0.1mm细化。则变压器内绕组在不同形变量 下的最大应力计算结果，如表2、表3所示。
[0095] 表2拐点附近(8-lOmm)处进行0.5mm细化结果
[0096]
[0099] 根据表2和表3得到变压器内绕组累积效应最大应力细化图，如图6所示。
[0100] 为了更好地验证变压器内绕组在多次短路冲击下累积效应的普遍性，本发明对不 同型号的电力变压器进行建模验证，则该变压器内绕组在不同形变量下的最大应力计算结 果，如表4所不。
[0101] 表4不同型号变压器内绕组累积效应计算结果
[0102]
[0103] 根据表4得到变压器内绕组累积效应最大应力图，如图5所示。
[0104] 从图7中可W看出，不同型号变压器内绕组在多次短路冲击下所受的最大应力同 样具有突变现象。运说明不同型号变压器在多次短路冲击下内绕组都具有相似的累积效 应，即变压器内绕组在多次短短冲击下所受最大应力有突变现象。
[0105] 变压器外绕组累积效应作用下的最大应力分析：
[0106] 变压器外绕组累积效应也是W2mm的形变量做研究，用相同的方法对变压器外绕 组在多次短路冲击下的累积效应进行研究，则不同形变量下的最大应力计算结果，如表5所 /J、- 〇
[0107] 表5变压器外绕组累积效应计算结果 [010 引
[0109] 根据表5得到变压器外绕组累积效应最大应力图，如图8所示。
[0110] 从图8可W看出，变压器外绕组在多次短路冲击下所受最大应力没有突变现象。
[0111] 本发明的基于变压器Ξ维电磁禪合模型的绕组应力评估方法，经过仿真计算和分 析的结果表明，该方法能够有效模拟变压器内部绕组磁场分布和受力情况，实现了本发明 目的和达到了所述的效果。
【主权项】
1. 一种基于变压器三维电磁耦合模型的绕组应力评估方法，其特征是，它包括以下步 骤： 1)变压器漏磁场的有限元计算 基于矢量磁位节点有限元法对变压器漏磁场进行计算，矢量磁位节点有限元法采用矢 量磁位A，忽略磁滞效应和涡流效应，根据Maxwel 1方程得到变压器非线性磁场方程：公式（1)中，μ为媒质的磁导率;J为电流密度，需要通过电路等效参数计算获得；X"， 代表旋度运算， 矢量磁位节点有限元法中节点单元的自由度为节点i的矢量磁位Ai，采用标量形状函数 N1，A为未知变量，则单元插值函数：公式(2)中，rinocb为单元节点数， 求解场域的整体插值函数：公式(3)中，Mn(n = l，2,…，nn)为基函数序列，由相关单元形状函数^对应叠加形成，η 为基函数序列通项编号，~为总项数，即总节点数，Αη为单元标量磁位； 对公式(2)应用格林定理，得伽辽金加权余量方程：公式(4)中，Mm(m= 1，2，…，ηη)为权函数序列，m为权函数序列通项编号，V表示边界体积 分量，S表示边界面积分量，边界面积分项为边界面磁场强度的切向分量的贡献，en为边界 面单位法向分量，在变压器计算模型中只涉及该项为零的情况，所以公式(4)整理得：将公式(3)代入公式(5)，针对全部权函数，可以将加权余量方程离散形成代数方程组， 通过求解可得所有节点上的矢量磁位A，利用磁场的微分公式β=ν^4,则区域内任意位置 的磁通密度都可以表示： 基于矢量磁位节点有限元方法计算漏磁场，可得变压器中漏磁场的最大漏磁场强度公 式，公式(6)中，I表示电流有效值，Bm为漏磁场强度，Hm为漏磁感应强度，Z为线圈匝数，P为 罗氏系数，μ〇为真空磁导率，H为线圈高度； 2)变压器绕组应力的计算 变压器绕组应力计算主要包含，绕组的平均环行拉应力、压缩应力和撑条或垫块之间 的导线辐向弯曲应力， 绕组受到辐向电磁力时，在直径为〇"的圆周上的压力：公式(7)中，Fx为两垫块间导线上的辐向力，Dn为圆周的平均直径，Hd为线饼高度，p为单 位面积上的压力； 绕组导线材料中产生的应力，公式(8)中，β为线饼的辐向宽度;σχ为导线辐向应力， 导线上的最大弯矩：公式(9)中，Fy为两垫块间导线上的轴向力，S为两垫块间导线上的长度， 导线为矩形截面的抗弯截面系数：公式(10)中，d为单根导线辐向厚度;b为单根导线轴向高度， 幅向垫块之间跨度内导线轴向弯曲应力：公式(11)中，σγ为导线轴向应力。
【专利摘要】一种基于变压器三维电磁耦合模型的绕组应力评估方法，其特点是，建立变压器三维电磁耦合模型并求解，以矢量磁位节点有限元方法为基础计算变压器漏磁场，通过模拟变压器内部磁场分布，对绕组所受到的应力进行计算分析。由于绕组发生形变后，使得绕组的几何结构不再规则，导致利用加载电流原理而编制的磁结构耦合模型不再适用。所以从与三维模型对应的二维模型中计算绕组受力值，然后将受力值通过相关的程序加载到对应的三维模型中来运算绕组所受应力。通过两撑条间绕组凹陷程度的不同来体现绕组的形变，进而观察不同累积情况下，模型所受到的应力情况。为变压器绕组结构稳定性的判定提供依据，具有科学合理，真实有效，实用价值高等优点。
【IPC分类】G06F17/50
【公开号】CN105488290
【申请号】CN201510916798
【发明人】潘超, 金明权, 吕嘉慧, 宋云东
【申请人】东北电力大学, 国网辽宁省电力有限公司电力科学研究院
【公开日】2016年4月13日
【申请日】2015年12月13日
【技术领域】
[0001] 本发明是一种基于变压器Ξ维电磁禪合模型的绕组应力评估方法，应用于电力变 压器绕组受力分析与结构稳定性评估。
【背景技术】
[0002] 电力变压器是电网和电力传输中最重要的设备之一，其绕组承受短路的能力直接 影响电网能否安全可靠的运行。在单次短路冲击作用下，绕组的稳定性取决于整定值与实 际受力值，当整定值大于实际受力值时，绕组是稳定的。在实际运行中变压器可能受到多次 短路冲击，绕组的抗短路可能会随着短路冲击次数的增加而降低，同时导线的应力可能会 增大，表现出一定的累积效应。由于在每次短路冲击作用下，绕组的结构形状都会发生形 变，绕组的不同形状又会使得整个磁场分布发生变化，进而导致绕组所受到的力发生变化。 根据结构力学可知，不同形状下绕组发生失稳的条件也将变化。因此研究不同累积情况下 绕组所受到的应力情况，运对电力变压器绕组结构稳定性评估有一定的指导意义和实际价 值。
【发明内容】
[0003] 本发明的目的是，提供一种科学合理，真实有效，实用价值高的基于变压器Ξ维电 磁禪合模型的绕组应力评估方法，能够根据实际变压器参数，建立变压器Ξ维电磁禪合模 型并求解，从与Ξ维模型对应的二维模型中计算得到绕组受力值，然后将受力值加载到对 应的Ξ维模型中计算得到绕组受到的应力。
[0004] 本发明的目的是由W下技术方案来实现的：一种基于变压器Ξ维电磁禪合模型的 绕组应力评估方法，其特征是，它包括W下步骤：
[0005] 1)变压器漏磁场的有限元计算
[0006] 基于矢量磁位节点有限元法对变压器漏磁场进行计算，矢量磁位节点有限元法采 用矢量磁位A，忽略磁滞效应和满流效应，根据Maxwel 1方程得到变压器非线性磁场方程：
[0007]
(1)
[000引公式（1)中，μ为媒质的磁导率;J为电流密度，需要通过电路等效参数计算获得； "▽X"代表旋度运算，
[0009] 矢量磁位节点有限元法中节点单元的自由度为节点i的矢量磁位Ai，采用标量形 状函数Ni，A为未知变量，则单元插值函数：
[0010]
(2)
[0011]公式（2)中，nnode为单元节点数，
[0012] 求解场域的整体插值函数：
[0013]
(3)
[0014] 公式(3)中，Μη(η=1，2,…，nn)为基函数序列，由相关单元形状函数Ni对应叠加形 成，η为基函数序列通项编号，nn为总项数，即总节点数，An为单元标量磁位；
[0015] 对公式(2)应用格林定理，得伽迂金加权余量方程：
[0016]
(4)
[0017] 公式(4)中，Mm(m=l，2,…，Πη)为权函数序列，m为权函数序列通项编号，V表示边界 体积分量，S表示边界面积分量，边界面积分项为边界面磁场强度的切向分量的贡献，en为 边界面单位法向分量，在变压器计算模型中只设及该项为零的情况，所W公式(4)整理得： [001 引
[0019] 将公式(3)代入公式(5)，针对全部权函数，可W将加权余量方程离散形成代数方 程组，通过求解可得所有节点上的矢量磁位Α，利用磁场的微分公式B = VXA，则区域内任 意位置的磁通密度都可W表示：
[0020] 基于矢量磁位节点有限元方法计算漏磁场，可得变压器中漏磁场的最大漏磁场强 度公式，
[0021]
(6)
[0022] 公式(6)中，I表示电流有效值，Bm为漏磁场强度，Hm为漏磁感应强度，Z为线圈应数， P为罗氏系数，μ〇为真空磁导率，Η为线圈高度；
[0023] 2)变压器绕组应力的计算
[0024] 变压器绕组应力计算主要包含，绕组的平均环行拉应力、压缩应力和撑条或垫块 之间的导线福向弯曲应力，
[0025] 绕组受到福向电磁力时，在直径为Dn的圆周上的压力：
[0026]
饼
[0027] 公式(7)中，Fx为两垫块间导线上的福向力，Dn为圆周的平均直径，出为线饼高度，P 为单位面积上的压力；
[0028] 绕组导线材料中产生的应力，
[0029]
(s)
[0030] 公式(8)中，β为线饼的福向宽度;οχ为导线福向应力，
[0031] 导线上的最大弯矩：
[0032]
側
[0033] 公式(9)中，Fy为两垫块间导线上的轴向力，S为两垫块间导线上的长度，
[0034] 导线为矩形截面的抗弯截面系数：
[0035]
(10)
[0036] 公式(10)中，d为单根导线福向厚度;b为单根导线轴向高度，
[0037] 幅向垫块之间跨度内导线轴向弯曲应力：
[00；3 引
(Π )
[0039] 公式(11)中，OY为导线轴向应力。
[0040] 本发明的基于变压器Ξ维电磁禪合模型的绕组应力评估方法，能够根据实际变压 器参数，建立变压器Ξ维电磁禪合模型并求解，W矢量磁位节点有限元法为基础计算变压 器漏磁场，通过模拟变压器内部磁场分布，对绕组所受到的应力进行计算分析。由于绕组发 生形变后，使得绕组的几何结构不再规则，导致利用加载电流原理而编制的磁结构禪合模 型不再适用。所W从与Ξ维模型对应的二维模型中计算绕组受力值，然后将受力值通过相 关的程序加载到对应的Ξ维模型中来运算绕组所受应力。为了更好地模拟绕组实际情况， 在绕组模型中设定了撑条，并且通过两撑条间绕组凹陷程度的不同来体现绕组的形变，进 而观察不同累积情况下，模型所受到的应力情况。能够为变压器绕组结构稳定性的判定提 供依据，具有科学合理，真实有效，实用价值高等优点。
【附图说明】
[0041] 图1是未发生形变的变压器绕组模型立体结构示意图；
[0042] 图2是图1的左视示意图；
[0043] 图3是发生形变的变压器绕组模型立体结构示意图；
[0044] 图4是图3的左视示意图；
[0045] 图5是变压器内绕组累积效应最大应力示意图；
[0046] 图6是变压器内绕组累积效应最大应力细化示意图；
[0047] 图7是不同型号变压器内绕组累积效应最大应力示意图；
[004引图8是变压器外绕组累积效应最大应力示意图。
【具体实施方式】
[0049] 下面结合附图和具体实施例对本发明作进一步描述：
[0050] 本发明的一种基于变压器Ξ维电磁禪合模型的绕组应力评估方法，包括W下步 骤：
[0051 ] 1)变压器漏磁场的有限元计算
[0052]基于矢量磁位节点有限元法对变压器漏磁场进行计算，矢量磁位节点有限元法采 用矢量磁位A，忽略磁滞效应和满流效应，根据Maxwel 1方程得到变压器非线性磁场方程： [0化3]
(1)
[0054] 公式（1)中，μ为媒质的磁导率;J为电流密度，需要通过电路等效参数计算获得； "▽X"代表旋度运算，
[0055] 矢量磁位节点有限元法中节点单元的自由度为节点i的矢量磁位Ai，采用标量形 状函数Ni，A为未知变量，则单元插值函数：
[0化6]
樹
[0057]公式（2)中，rinode为单元节点数，
[005引求解场域的整体插值函数：
[0化9]
(3)
[0060」公式(3)中，Μη(η=1，2,…，nn)为基函数序列，由相关单元形状函数Ni对应叠加形 成，η为基函数序列通项编号，Πη为总项数，即总节点数，An为单元标量磁位；
[0061] 对公式(2)应用格林定理，得伽迂金加权余量方程：
[0062]
(4)
[0063] 公式(4)中，Mm(m= 1,2,…，Πη)为权函数序列，m为权函数序列通项编号，V表示边界 体积分量，S表示边界面积分量，边界面积分项为边界面磁场强度的切向分量的贡献，en为 边界面单位法向分量，在变压器计算模型中只设及该项为零的情况，所W公式(4)整理得：
[0064]
(勾
[0065] 将公式(3)代入公式(5)，针对全部权函数，可W将加权余量方程离散形成代数方 程组，通过求解可得所有节点上的矢量磁位A，利用磁场的微 分公式B = VXA，则区域内任 意位置的磁通密度都可W表示：
[0066] 基于矢量磁位节点有限元方法计算漏磁场，可得变压器中漏磁场的最大漏磁场强 度公式，
[0067]
批）
[0068] 公式(6)中，I表示电流有效值，Bm为漏磁场强度，Hm为漏磁感应强度，Z为线圈应数， P为罗氏系数，μ〇为真空磁导率，Η为线圈高度；
[0069] 2)变压器绕组应力的计算
[0070] 变压器绕组应力计算主要包含，绕组的平均环行拉应力、压缩应力和撑条或垫块 之间的导线福向弯曲应力，
[0071] 绕组受到福向电磁力时，在直径为Dn的圆周上的压力：
[0072]
(7)
[0073] 公式(7)中，Fx为两垫块间导线上的福向力，Dn为圆周的平均直径，出为线饼高度，P 为单位面积上的压力；
[0074] 绕组导线材料中产生的应力，
[0075]
做
[0076] 公式(8)中，β为线饼的福向宽度;οχ为导线福向应力，
[0077] 导线上的最大弯矩：
[007引
(9)
[0079] 公式(9)中，Fy为两垫块间导线上的轴向力，S为两垫块间导线上的长度，
[0080] 导线为矩形截面的抗弯截面系数： 「00811
(10)
[0082] 公式(10)中，d为单根导线福向厚度;b为单根导线轴向高度，
[0083] 幅向垫块之间跨度内导线轴向弯曲应力：
[0084]
(11)
[0085] 公式（11)中，OY为导线轴向应力。
[0086] 变压器内绕组累积效应作用下的最大应力分析：
[0087] 1)变压器内绕组累积效应模型的建立
[0088] 假设档距之间凹陷的最大变形为2mm，即每经过一次短路冲击作用后，绕组凹陷形 变的最大处在原来的基础上形变2mm，W此为基础得到绕组的二维模型，然后将计算得到的 受力值施加在相对应的Ξ维模型上，进而得出所需的结论。未发生形变的变压器绕组模型 立体结构图和左视图，如图1、图2所示；发生形变的变压器绕组模型立体结构图和左视图， 如图3、图4所示。
[0089] 2)变压器内绕组累积效应
[0090] 本发明对变压器内绕组的累积情况均是W2mm的形变量做研究的，并且假设变压 器内绕组经受过五次短路冲击，则变压器内绕组在不同形变量下的最大应力计算结果，如 表1所示。
[0091] 表1变压器内绕组累积效应计算结果
[0092]
[0093] 由表1得到变压器内绕组累积效应最大应力图，如图5和图6所示。
[0094] 由变压器内绕组累积效应最大应力图可知绕组所受应力在形变量为8mm附近处发 生突变。因此为了使研究结果更准确，本发明在上述形变的基础上，对拐点附近(8-lOmm)处 进行0.5mm细化;对突变点附近(9-9.6mm)处进行0.1mm细化。则变压器内绕组在不同形变量 下的最大应力计算结果，如表2、表3所示。
[0095] 表2拐点附近(8-lOmm)处进行0.5mm细化结果
[0096]
[0099] 根据表2和表3得到变压器内绕组累积效应最大应力细化图，如图6所示。
[0100] 为了更好地验证变压器内绕组在多次短路冲击下累积效应的普遍性，本发明对不 同型号的电力变压器进行建模验证，则该变压器内绕组在不同形变量下的最大应力计算结 果，如表4所不。
[0101] 表4不同型号变压器内绕组累积效应计算结果
[0102]
[0103] 根据表4得到变压器内绕组累积效应最大应力图，如图5所示。
[0104] 从图7中可W看出，不同型号变压器内绕组在多次短路冲击下所受的最大应力同 样具有突变现象。运说明不同型号变压器在多次短路冲击下内绕组都具有相似的累积效 应，即变压器内绕组在多次短短冲击下所受最大应力有突变现象。
[0105] 变压器外绕组累积效应作用下的最大应力分析：
[0106] 变压器外绕组累积效应也是W2mm的形变量做研究，用相同的方法对变压器外绕 组在多次短路冲击下的累积效应进行研究，则不同形变量下的最大应力计算结果，如表5所 /J、- 〇
[0107] 表5变压器外绕组累积效应计算结果 [010 引
[0109] 根据表5得到变压器外绕组累积效应最大应力图，如图8所示。
[0110] 从图8可W看出，变压器外绕组在多次短路冲击下所受最大应力没有突变现象。
[0111] 本发明的基于变压器Ξ维电磁禪合模型的绕组应力评估方法，经过仿真计算和分 析的结果表明，该方法能够有效模拟变压器内部绕组磁场分布和受力情况，实现了本发明 目的和达到了所述的效果。
【主权项】
1. 一种基于变压器三维电磁耦合模型的绕组应力评估方法，其特征是，它包括以下步 骤： 1)变压器漏磁场的有限元计算 基于矢量磁位节点有限元法对变压器漏磁场进行计算，矢量磁位节点有限元法采用矢 量磁位A，忽略磁滞效应和涡流效应，根据Maxwel 1方程得到变压器非线性磁场方程：公式（1)中，μ为媒质的磁导率;J为电流密度，需要通过电路等效参数计算获得；X"， 代表旋度运算， 矢量磁位节点有限元法中节点单元的自由度为节点i的矢量磁位Ai，采用标量形状函数 N1，A为未知变量，则单元插值函数：公式(2)中，rinocb为单元节点数， 求解场域的整体插值函数：公式(3)中，Mn(n = l，2,…，nn)为基函数序列，由相关单元形状函数^对应叠加形成，η 为基函数序列通项编号，~为总项数，即总节点数，Αη为单元标量磁位； 对公式(2)应用格林定理，得伽辽金加权余量方程：公式(4)中，Mm(m= 1，2，…，ηη)为权函数序列，m为权函数序列通项编号，V表示边界体积 分量，S表示边界面积分量，边界面积分项为边界面磁场强度的切向分量的贡献，en为边界 面单位法向分量，在变压器计算模型中只涉及该项为零的情况，所以公式(4)整理得：将公式(3)代入公式(5)，针对全部权函数，可以将加权余量方程离散形成代数方程组， 通过求解可得所有节点上的矢量磁位A，利用磁场的微分公式β=ν^4,则区域内任意位置 的磁通密度都可以表示： 基于矢量磁位节点有限元方法计算漏磁场，可得变压器中漏磁场的最大漏磁场强度公 式，公式(6)中，I表示电流有效值，Bm为漏磁场强度，Hm为漏磁感应强度，Z为线圈匝数，P为 罗氏系数，μ〇为真空磁导率，H为线圈高度； 2)变压器绕组应力的计算 变压器绕组应力计算主要包含，绕组的平均环行拉应力、压缩应力和撑条或垫块之间 的导线辐向弯曲应力， 绕组受到辐向电磁力时，在直径为〇"的圆周上的压力：公式(7)中，Fx为两垫块间导线上的辐向力，Dn为圆周的平均直径，Hd为线饼高度，p为单 位面积上的压力； 绕组导线材料中产生的应力，公式(8)中，β为线饼的辐向宽度;σχ为导线辐向应力， 导线上的最大弯矩：公式(9)中，Fy为两垫块间导线上的轴向力，S为两垫块间导线上的长度， 导线为矩形截面的抗弯截面系数：公式(10)中，d为单根导线辐向厚度;b为单根导线轴向高度， 幅向垫块之间跨度内导线轴向弯曲应力：公式(11)中，σγ为导线轴向应力。
【专利摘要】一种基于变压器三维电磁耦合模型的绕组应力评估方法，其特点是，建立变压器三维电磁耦合模型并求解，以矢量磁位节点有限元方法为基础计算变压器漏磁场，通过模拟变压器内部磁场分布，对绕组所受到的应力进行计算分析。由于绕组发生形变后，使得绕组的几何结构不再规则，导致利用加载电流原理而编制的磁结构耦合模型不再适用。所以从与三维模型对应的二维模型中计算绕组受力值，然后将受力值通过相关的程序加载到对应的三维模型中来运算绕组所受应力。通过两撑条间绕组凹陷程度的不同来体现绕组的形变，进而观察不同累积情况下，模型所受到的应力情况。为变压器绕组结构稳定性的判定提供依据，具有科学合理，真实有效，实用价值高等优点。
【IPC分类】G06F17/50
【公开号】CN105488290
【申请号】CN201510916798
【发明人】潘超, 金明权, 吕嘉慧, 宋云东
【申请人】东北电力大学, 国网辽宁省电力有限公司电力科学研究院
【公开日】2016年4月13日
【申请日】2015年12月13日
最新回复(0)