CN114238850A - 换流变有载分接开关切换芯子磁-热耦合计算方法及系统 - Google Patents

换流变有载分接开关切换芯子磁-热耦合计算方法及系统 Download PDF

Info

Publication number
CN114238850A
CN114238850A CN202111236750.5A CN202111236750A CN114238850A CN 114238850 A CN114238850 A CN 114238850A CN 202111236750 A CN202111236750 A CN 202111236750A CN 114238850 A CN114238850 A CN 114238850A
Authority
CN
China
Prior art keywords
changer
load tap
magnetic
thermal coupling
analysis
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
CN202111236750.5A
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.)
Shanghai Jiaotong University
Shanghai Huaming Power Equipment Co Ltd
Maintenance and Test Center of Extra High Voltage Power Transmission Co
Original Assignee
Shanghai Jiaotong University
Shanghai Huaming Power Equipment Co Ltd
Maintenance and Test Center of Extra High Voltage Power Transmission Co
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 Shanghai Jiaotong University, Shanghai Huaming Power Equipment Co Ltd, Maintenance and Test Center of Extra High Voltage Power Transmission Co filed Critical Shanghai Jiaotong University
Priority to CN202111236750.5A priority Critical patent/CN114238850A/zh
Publication of CN114238850A publication Critical patent/CN114238850A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations

Abstract

本发明公开了一种换流变有载分接开关切换芯子磁‑热耦合计算方法及系统，该方法包括步骤：获取有载分接开关各个部件的尺寸信息，并获取有载分接开关各个部件对应的磁导率；基于所获取的尺寸信息和磁导率建立有载分接开关的三维模型；基于所述三维模型确定有载分接开关磁‑热耦合分析的线性方程组，并从线性方程组确定有载分接开关各个部件的磁‑热耦合分布。通过有载分接开关的磁‑热耦合分析，根据线性方程组计算有载分接开关，可确定其磁热分布用于产品设计与运行中异常过热分析，从而判断其过热集中区域和热老化薄弱点，便于对有载分接开关各零部件的设计做出优化，以提高有载分接开关的质量；同时，提高对各部件的磁热耦合检测精度。

Description

换流变有载分接开关切换芯子磁-热耦合计算方法及系统
技术领域
本发明涉及换流变压器有载分接开关领域，尤其是涉及一种换流变有载分接开关切换芯子磁-热耦合计算方法及系统。
背景技术
换流变是高电压直流输电技术中最主要的一次电力设备，主要作用是提供特殊要求的电源，其主要参数由直流系统的特殊要求以及所联结的交流系统参数而确定。在整流站，用换流变压器将交流系统和直流系统隔离，通过换流装置将交流网路的电能转换为高压直流电能，送到高压直流输电线路；在逆变站，通过换流装置将直流电能转换为交流电能，再通过换流变压器转换为正常交流正弦电压，送到交流电网；从而实现交流输电网路与高压直流输电线路的联络。
换流变压器有载分接开关动作频率高、动作时序复杂与快速，导致有载分接开关设计与制造加工精度高、装配精度远高于一般的电气设备。
现有中一般采用电流产生的焦耳热直接施加于分接开关的油室进行温升计算，因分接开关采用大量的铸铝机构件，实际运行中存在触头电流导致的涡流附加损耗，当前的单一场的计算方法无法准确计算该涡流附加损耗引起的附加温升，限制了分接开关的可靠设计，因此需要对此进行改进。
发明内容
为了提高对有载分接开关各部件的磁热耦合检测精度，本申请提供一种换流变有载分接开关切换芯子磁-热耦合计算方法及系统。
第一方面，本申请提供一种换流变有载分接开关切换芯子磁-热耦合计算方法，采用如下的技术方案：
一种换流变有载分接开关切换芯子磁-热耦合计算方法，包括以下步骤：
获取有载分接开关各个部件的尺寸信息，并获取有载分接开关各个部件对应的磁导率；
基于所获取的尺寸信息和磁导率建立有载分接开关的三维模型；
基于所述三维模型确定有载分接开关磁-热耦合分析的线性方程组，并从线性方程组确定有载分接开关各个部件的磁-热耦合分布。
通过采用上述技术方案，有载分接开关由各部分零部件组成，获取组成有载分接开关各部分零件尺寸和磁导率等数据，根据所获取的数据建立关于有载分接开关的三维模型，采用三维模型中的相关数据得出有载分接开关的磁-热耦合分析的线性方程组，可确定有载分接开关磁热分布用于换流变有载分接开关的产品设计与运行中异常过热分析，从而判断换流变有载分接开关过热集中区域和热老化薄弱点，便于对有载分接开关各零部件的设计做出优化，以提高有载分接开关的质量；同时，该方法计算精度高，耦合性能好，以及求解效率高，提高对有载分接开关各部件的磁热耦合检测精度。
可选的，所述线性方程组通过以下方法计算：
对所述三维模型进行自适应网格剖分，以获取有载分接开关磁-热耦合分析的节点数N和单元数M；
根据所述节点数N和单元数M确定所述线性方程组：
Figure DEST_PATH_IMAGE002
；其中，
Figure DEST_PATH_IMAGE004
表示磁分析的系数矩阵，包含N*N个元素，
Figure 293749DEST_PATH_IMAGE004
中每个元素用
Figure DEST_PATH_IMAGE006
表示；H表示温度分布的系数矩阵，包含N*N个元素，H为中每个元素用
Figure DEST_PATH_IMAGE008
表示；
Figure DEST_PATH_IMAGE010
表示磁分析的已知向量，大小为N，
Figure 643697DEST_PATH_IMAGE010
中每个元素用
Figure DEST_PATH_IMAGE012
表示；
Figure DEST_PATH_IMAGE014
表示温度分布的已知向量，大小为N，
Figure 528607DEST_PATH_IMAGE014
中每个元素用
Figure DEST_PATH_IMAGE016
表示；
Figure DEST_PATH_IMAGE018
表示磁分析的未知向量，大小为N，
Figure 14821DEST_PATH_IMAGE018
中每个元素用
Figure DEST_PATH_IMAGE020
表示；
Figure DEST_PATH_IMAGE022
表示温度分布的未知向量，大小为N，K中每个元素用
Figure DEST_PATH_IMAGE024
表示；
Figure 225354DEST_PATH_IMAGE006
、
Figure 485434DEST_PATH_IMAGE008
、
Figure 95407DEST_PATH_IMAGE012
和
Figure 373810DEST_PATH_IMAGE016
通过以下方法计算：
Figure DEST_PATH_IMAGE026
；
Figure DEST_PATH_IMAGE028
；
Figure DEST_PATH_IMAGE030
；
Figure DEST_PATH_IMAGE032
；
其中，
Figure DEST_PATH_IMAGE034
表示数组B中
Figure DEST_PATH_IMAGE036
行描述的空间区域，
Figure DEST_PATH_IMAGE038
表示数组A中
Figure 348719DEST_PATH_IMAGE036
行数据建立的形状函数，
Figure DEST_PATH_IMAGE040
表示数组A中
Figure DEST_PATH_IMAGE042
行数据建立的形状函数，
Figure DEST_PATH_IMAGE044
表示有载分接开关的磁边界条件，
Figure DEST_PATH_IMAGE046
表示有载分接开关热边界条件
Figure DEST_PATH_IMAGE048
通过采用上述技术方案，对有载分接开关的三维模型进行自适应网格剖分，以得到有载分接开关磁-热耦合分析的节点数N和单元数M，节点数N表示该模型运用了N个三维空间点描述有载分接开关磁-热耦合分析的三维图形模型，单元数M表示该模型运用了M个三维四面体单元描述有载分接开关磁-热耦合分析的三维图形模型，分别用数组A和数组B存储节点和单元基本信息，数组A的每一行包含一个节点编号和该点相对有载分接开关主轴底部的坐标值，数组B的每一行包含一个三维四面体单元编号和4个不同节点编号来描述部分有载分接开关磁-热耦合分析的三维图形模型，数组B中所有单元构成整个有载分接开关磁-热耦合分析的三维图形模型；将N和M代入到公式中进行求解，便可以得到线性方程组。
可选的，所述线性方程组中的未知向量通过以下方法计算：
根据所述节点初始向量C0和H0，
将磁分析的系数矩阵D分解成M、B和N三个矩阵，通过以下方法分解：
如果i=j，则
Figure DEST_PATH_IMAGE050
=
Figure DEST_PATH_IMAGE052
，否则
Figure 814205DEST_PATH_IMAGE050
=0；
如果i>j，则
Figure DEST_PATH_IMAGE054
=
Figure DEST_PATH_IMAGE056
，否则
Figure 551172DEST_PATH_IMAGE054
=0；
如果i<j，则
Figure DEST_PATH_IMAGE058
=
Figure 716705DEST_PATH_IMAGE056
，否则
Figure 190412DEST_PATH_IMAGE058
=0；
其中，i为自然数，且0<i<N；j为自然数，且0<j<N；
将温度分布的系数矩阵H分解成S、V和Q三个矩阵，通过以下方法分解：
如果i=j，则
Figure DEST_PATH_IMAGE060
=
Figure DEST_PATH_IMAGE062
，否则
Figure 471089DEST_PATH_IMAGE060
=0；
如果i>j，则
Figure DEST_PATH_IMAGE064
=
Figure DEST_PATH_IMAGE066
，否则
Figure 485182DEST_PATH_IMAGE064
=0；
如果i<j，则
Figure DEST_PATH_IMAGE068
=
Figure 239642DEST_PATH_IMAGE066
，否则
Figure 618671DEST_PATH_IMAGE068
=0；
其中，i为自然数，且0<i<N；j为自然数，且0<j<N；
得到从
Figure 75060DEST_PATH_IMAGE050
、
Figure 830526DEST_PATH_IMAGE054
、
Figure 938029DEST_PATH_IMAGE058
、
Figure 753538DEST_PATH_IMAGE060
、
Figure 697223DEST_PATH_IMAGE064
和
Figure 990801DEST_PATH_IMAGE068
，计算
Figure DEST_PATH_IMAGE070
和
Figure DEST_PATH_IMAGE072
的公式：
Figure DEST_PATH_IMAGE074
（1-1）
Figure DEST_PATH_IMAGE076
（1-2）
Figure 47750DEST_PATH_IMAGE070
=
Figure DEST_PATH_IMAGE078
（1-3）
Figure 814587DEST_PATH_IMAGE072
=
Figure DEST_PATH_IMAGE080
--
Figure DEST_PATH_IMAGE082
（1-4）
其中，w为自然数，且0<12000，重复迭代计算上述（1-1）-（1-4）的公式，直到
Figure DEST_PATH_IMAGE084
且
Figure DEST_PATH_IMAGE086
。
通过采用上述技术方案，将磁分析的系数矩阵和温度分布的系数矩阵分别分解成三个矩阵，将分解得到的
Figure 120934DEST_PATH_IMAGE050
、
Figure 218203DEST_PATH_IMAGE054
、
Figure 34718DEST_PATH_IMAGE058
、
Figure 660872DEST_PATH_IMAGE060
、
Figure 579149DEST_PATH_IMAGE064
和
Figure 214530DEST_PATH_IMAGE068
代入到公式中，通过数字迭代方法计算出
Figure 901863DEST_PATH_IMAGE070
和
Figure 980809DEST_PATH_IMAGE072
，从而得出有载分接开关磁热分析的未知向量。
可选的，所述方法还包括：
通过以下公式计算有载分接开关三维模型中任一点的电场信息：
Figure DEST_PATH_IMAGE088
；
其中，
Figure DEST_PATH_IMAGE090
表示有载分接开关磁-热耦合分析中的某一点电位，
Figure DEST_PATH_IMAGE092
表示空间中j单元的体积,M为三维图形模型中的剖分的单元数，
Figure 963546DEST_PATH_IMAGE038
和
Figure 668197DEST_PATH_IMAGE040
是根据所述数组A中i行和j行数据建立的形状函数；
根据所述电场信息生成有载分接开关的电场云图。
通过采用上述技术方案，根据有载分接开关磁-热耦合分析的未知向量的解，对应于有载分接开关磁-热耦合分析的三维图形模型中的剖分的节点N和单元数M,有载分接开关磁-热耦合分析的三维图形模型中任意一点电场分布可通过上述计算公式获得，并将值绘制成空间等位云图。
第二方面，本申请提供一种换流变有载分接开关切换芯子磁-热耦合计算系统，采用如下的技术方案：
所述系统包括参数获取模块、三维建模模块和磁-热耦合分布计算模块：
所述参数获取模块用于获取有载分接开关各个部件的尺寸信息，并获取有载分接开关各个部件对应的磁导率；
所述三维建模模块与参数获取模块连接，用于接收所述尺寸信息和磁导率，基于所述尺寸信息和磁导率建立有载分接开关的三维模型；
所述磁-热耦合分布计算模块与三维建模模块连接，用于接收所述三维模型，基于所述三维模型确定有载分接开关磁-热耦合分析的线性方程组，并从所述线性方程组确定有载分接开关各个部件的磁-热耦合分布。
通过采用上述技术方案，有载分接开关由各部分零部件组成，获取组成有载分接开关各部分零件尺寸和磁导率等数据，根据所获取的数据建立关于有载分接开关的三维模型，采用三维模型中的相关数据得出有载分接开关的磁-热耦合分析的线性方程组，可确定有载分接开关磁热分布用于换流变有载分接开关的产品设计与运行中异常过热分析，从而判断换流变有载分接开关过热集中区域和热老化薄弱点，便于对有载分接开关各零部件的设计做出优化，以提高有载分接开关的质量；同时，该方法计算精度高，耦合性能好，以及求解效率高，提高对有载分接开关各部件的磁热耦合检测精度。
第三方面，本申请提供一种换流变有载分接开关切换芯子磁-热耦合计算系统，采用如下的技术方案：
一种应用上述任一所述的一种换流变有载分接开关切换芯子磁-热耦合计算方法，该方法用于±500kV直流输电工程整流侧或者逆变侧换流变有载分接开关磁热分析确定。
可选的，该方法用于±800kV直流输电工程整流侧或者逆变侧换流变有载分接开关磁热分析确定。
第四方面，本申请提供一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述任一项所述一种换流变有载分接开关切换芯子磁-热耦合计算方法的步骤。
第五方面，本申请提供一种计算机可读存储介质，存储有能够被处理器加载并执行上述任一种所述方法的计算机程序。
综上所述，本申请具有以下有益效果：
通过有载分接开关的磁-热耦合分析，有载分接开关由各部分零部件组成，获取组成有载分接开关各部分零件尺寸和磁导率等数据，根据所获取的数据建立关于有载分接开关的三维模型，采用三维模型中的相关数据得出有载分接开关的磁-热耦合分析的线性方程组，可确定有载分接开关磁热分布用于换流变有载分接开关的产品设计与运行中异常过热分析，从而判断换流变有载分接开关过热集中区域和热老化薄弱点，便于对有载分接开关各零部件的设计做出优化，以提高有载分接开关的质量；同时，该方法计算精度高，耦合性能好，以及求解效率高，提高对有载分接开关各部件的磁热耦合检测精度。
附图说明
图1是本申请其中一实施例中换流变压器有载分接开关切换芯子磁-热耦合计算方法的原理图；
图2是本申请其中一实施例中换流变压器有载分接开关的整体结构示意图；
图3是本申请其中一实施例中换流变压器有载分接开关的三维模型自适应网格剖分示意图；
图4是本申请其中一实施例中换流变压器有载分接开关的温度场分布示意图；
图5是本申请其中一实施例中单侧、双侧和中性点触头的温度场分布示意图；
图6是本申请其中一实施例中换流变压器有载分接开关切换芯子磁-热耦合计算系统的框图。
具体实施方式
以下结合附图1-6对本申请作进一步详细说明。
本申请实施例公开一种换流变有载分接开关切换芯子磁-热耦合计算方法，参见图1，该方法包括以下步骤：
S1、获取有载分接开关各个部件的尺寸信息，并获取有载分接开关各个部件对应的磁导率。
在本实施例中，各个部件是指构成有载分接开关的各个零部件；磁导率是指各个零部件磁介质磁性的物理量。
具体地，有载分接开关模型参数的获取工具，包括测量有载分接开关几何模型尺寸的长卷尺和各部件相对介电常数的测量仪。有载分接开关几何模型尺寸包含以下方面：测量有载分接开关绝缘筒的长、宽和高；测量单侧触头的长、宽和高；测量双侧触头的长、宽和高；测量中性点触头的长、宽和高；测量真空灭弧室的直径和高度；测量绝缘支架的长、宽和高；测量隔离刀闸的内径、外径和高度。
换流变有载分接开关各部件的磁导率，包括有载分接开关绝缘筒的磁导率，单侧触头、双侧触头、中性点触头和绝缘支架的磁导率，以及真空灭弧室的磁导率。
其中测量换流变压器有载分接开关几何模型尺寸的长卷尺、各部件相对介电常数的测量仪为市售产品。方法步骤如下：
换流变压器有载分接开关模型参数的获取：通过长卷尺测量换流变有载分接开关尺寸如表1:
表1：换流变几何模型尺寸
Figure DEST_PATH_IMAGE094
通过磁导率测量仪测量有载分接开关各部件磁导率如表2：
表2：换流变各部件磁导率
测量部位 相对介电常数
绝缘筒 /
绝缘支架 /
单侧触头、双侧触头、中性点触头 2.5
真空灭弧室 2
S2、基于所获取的尺寸信息和磁导率建立有载分接开关的三维模型。
在本实施例中，有载分接开关的三维模型包括各零部件组成的单元模块以及有载分接开关整体的三维模型。
参见图2，具体地，使用现有的三维建模软件，在本实施例中，通过三维建模软件绘制换流变压器的有载分接开关磁-热耦合分析的三维图形模型，输入所获取的有载分接开关的尺寸信息和磁导率，建立关于有载分接开关的三维模型。
换流变铁芯的三维图形绘制包括根据有载分接开关绝缘筒的长、宽和高绘制绝缘筒的三维图形。换流变绕组的三维图形绘制包括根据测量单侧触头的长、宽和高绘制单侧触头的三维图形；根据测量双侧触头的长、宽和高绘制双侧触头的三维图形；根据测量中性点触头的长、宽和高绘制中性点触头的三维图形。换流变铁芯的三维图形绘制包括根据有载分接开关真空灭弧室的直径和高度绘制真空灭弧室的三维图形。换流变铁芯的三维图形绘制包括根据有载分接开关测量隔离刀闸的内径、外径和高度绘制隔离刀闸的三维图形。
S3、基于三维模型确定有载分接开关磁-热耦合分析的线性方程组，并从线性方程组确定有载分接开关各个部件的磁-热耦合分布。
在本实施中，磁-热耦合分析是指有载分接开关的各零部件在组装状态下以及运动状态下的具体磁热分布情况。
具体地，参照图3，使用三维建模软件的分解功能对有载分接开关的三维模型进行自适应网格剖分，得到有载分接开关磁-热耦合分析的节点数N和单元数M，节点数N表示该模型运用了N个三维空间点描述有载分接开关磁-热耦合分析的三维图形模型，单元数M表示该模型运用了M个三维四面体单元描述有载分接开关磁-热耦合分析的三维图形模型，分别用数组A和数组B存储节点和单元基本信息，数组A的每一行包含一个节点编号和该点相对有载分接开关主轴底部的坐标值，数组B的每一行包含一个三维四面体单元编号和4个不同节点编号来描述部分有载分接开关磁-热耦合分析的三维图形模型，数组B中所有单元构成整个有载分接开关磁-热耦合分析的三维图形模型。
根据节点数N和单元数M确定线性方程组：
Figure 85403DEST_PATH_IMAGE002
；其中，
Figure 584517DEST_PATH_IMAGE004
表示磁分析的系数矩阵，包含N*N个元素，
Figure 461075DEST_PATH_IMAGE004
中每个元素用
Figure 438258DEST_PATH_IMAGE006
表示；H表示温度分布的系数矩阵，包含N*N个元素，H为中每个元素用
Figure 100184DEST_PATH_IMAGE008
表示；
Figure 239041DEST_PATH_IMAGE010
表示磁分析的已知向量，大小为N，
Figure 353628DEST_PATH_IMAGE010
中每个元素用
Figure 885234DEST_PATH_IMAGE012
表示；
Figure 401666DEST_PATH_IMAGE014
表示温度分布的已知向量，大小为N，
Figure 242583DEST_PATH_IMAGE014
中每个元素用
Figure 844466DEST_PATH_IMAGE016
表示；
Figure 678299DEST_PATH_IMAGE018
表示磁分析的未知向量，大小为N，
Figure 49237DEST_PATH_IMAGE018
中每个元素用
Figure 61056DEST_PATH_IMAGE020
表示；
Figure 884655DEST_PATH_IMAGE022
表示温度分布的未知向量，大小为N，K中每个元素用
Figure 23643DEST_PATH_IMAGE024
表示；
Figure 717930DEST_PATH_IMAGE006
、
Figure 431808DEST_PATH_IMAGE008
、
Figure 8283DEST_PATH_IMAGE012
和
Figure 195637DEST_PATH_IMAGE016
通过以下方法计算：
Figure 275588DEST_PATH_IMAGE026
；
Figure 894789DEST_PATH_IMAGE028
；
Figure 161822DEST_PATH_IMAGE030
；
Figure 626301DEST_PATH_IMAGE032
；
其中，
Figure 577071DEST_PATH_IMAGE034
表示数组B中
Figure 367172DEST_PATH_IMAGE036
行描述的空间区域，
Figure 387081DEST_PATH_IMAGE038
表示数组A中
Figure 920830DEST_PATH_IMAGE036
行数据建立的形状函数，
Figure 224642DEST_PATH_IMAGE040
表示数组A中
Figure 920065DEST_PATH_IMAGE042
行数据建立的形状函数，
Figure 692849DEST_PATH_IMAGE044
表示有载分接开关的磁边界条件，
Figure 764710DEST_PATH_IMAGE046
表示有载分接开关热边界条件。
进一步地，线性方程组中的未知向量通过以下方法计算：
根据节点初始向量C0和H0，初始向量C0和H0如下：
Figure DEST_PATH_IMAGE096
将磁分析的系数矩阵D分解成M、B和N三个矩阵，通过以下方法分解：
如果i=j，则
Figure 955651DEST_PATH_IMAGE050
=
Figure 87556DEST_PATH_IMAGE052
，否则
Figure 82056DEST_PATH_IMAGE050
=0；
如果i>j，则
Figure 206876DEST_PATH_IMAGE054
=
Figure 236012DEST_PATH_IMAGE056
，否则
Figure 273238DEST_PATH_IMAGE054
=0；
如果i<j，则
Figure 20614DEST_PATH_IMAGE058
=
Figure 185010DEST_PATH_IMAGE056
，否则
Figure 599811DEST_PATH_IMAGE058
=0；
其中，i为自然数，且0<i<N；j为自然数，且0<j<N；
将温度分布的系数矩阵H分解成S、V和Q三个矩阵，通过以下方法分解：
如果i=j，则
Figure 276780DEST_PATH_IMAGE060
=
Figure 245873DEST_PATH_IMAGE062
，否则
Figure 712496DEST_PATH_IMAGE060
=0；
如果i>j，则
Figure 450644DEST_PATH_IMAGE064
=
Figure 95252DEST_PATH_IMAGE066
，否则
Figure 551642DEST_PATH_IMAGE064
=0；
如果i<j，则
Figure 57840DEST_PATH_IMAGE068
=
Figure 916075DEST_PATH_IMAGE066
，否则
Figure 731584DEST_PATH_IMAGE068
=0；
其中，i为自然数，且0<i<N；j为自然数，且0<j<N；
得到从
Figure 409690DEST_PATH_IMAGE050
、
Figure 749274DEST_PATH_IMAGE054
、
Figure 275064DEST_PATH_IMAGE058
、
Figure 261474DEST_PATH_IMAGE060
、
Figure 692456DEST_PATH_IMAGE064
和
Figure 789725DEST_PATH_IMAGE068
，计算
Figure 606240DEST_PATH_IMAGE070
和
Figure 294710DEST_PATH_IMAGE072
的公式：
Figure 947408DEST_PATH_IMAGE074
（1-1）
Figure 599101DEST_PATH_IMAGE076
（1-2）
Figure 286434DEST_PATH_IMAGE070
=
Figure 614647DEST_PATH_IMAGE078
（1-3）
Figure 754642DEST_PATH_IMAGE072
=
Figure 442981DEST_PATH_IMAGE080
--
Figure 984821DEST_PATH_IMAGE082
（1-4）
其中，w为自然数，且0<12000，重复迭代计算上述（1-1）-（1-4）的公式，直到
Figure DEST_PATH_IMAGE097
且
Figure 15093DEST_PATH_IMAGE086
。
进一步地，参照图4，根据有载分接开关磁-热耦合分析的未知向量的解对应于有载分接开关磁-热耦合分析的三维模型中的剖分的节点N和单元数M。通过以下公式计算有载分接开关三维模型中任一点的电场信息：
Figure DEST_PATH_IMAGE098
；其中，
Figure 924275DEST_PATH_IMAGE090
表示有载分接开关磁-热耦合分析中的某一点电位，
Figure 167037DEST_PATH_IMAGE092
表示空间中j单元的体积,M为三维图形模型中的剖分的单元数，
Figure 828963DEST_PATH_IMAGE038
和
Figure 748246DEST_PATH_IMAGE040
是根据数组A中i行和j行数据建立的形状函数；根据电场信息生成有载分接开关的电场云图。并绘制成空间等位云图，如图5，具体包括单侧触头、双侧触头和中性点触头。
举例说明，在本实施例中，对有载分接开关进行自适应网格剖分之后，得到换流变有载分接开关电场分析的节点数246782和单元数1345231，并形成节点和单元基本信息数值A和B（其中矩阵中的M表示的因数值太多，无法一一列举的省略符号）；
Figure DEST_PATH_IMAGE100
Figure DEST_PATH_IMAGE102
对应的，便可以得到线性方程组为：
Figure DEST_PATH_IMAGE104
Figure DEST_PATH_IMAGE106
根据有载分接开关磁-热耦合分析的线性方程组，依据如下流程与步骤求解未知向量：
Figure DEST_PATH_IMAGE108
Figure DEST_PATH_IMAGE110
本申请实施例还公开一种换流变有载分接开关切换芯子磁-热耦合计算系统，参照图6，该系统包括参数获取模块、三维建模模块和磁-热耦合分布计算模块：
参数获取模块用于获取有载分接开关各个部件的尺寸信息，并获取有载分接开关各个部件对应的磁导率；
三维建模模块与参数获取模块连接，用于接收尺寸信息和磁导率，基于尺寸信息和磁导率建立有载分接开关的三维模型；
磁-热耦合分布计算模块与三维建模模块连接，用于接收三维模型，基于三维模型确定有载分接开关磁-热耦合分析的线性方程组，并从线性方程组确定有载分接开关各个部件的磁-热耦合分布。
进一步地，另一实施例中，磁-热耦合分布计算模块包括磁温向量计算单元，对三维模型进行自适应网格剖分以获取有载分接开关磁-热耦合分析的节点数N和单元数M；根据节点数N和单元数M确定线性方程组：
Figure 733606DEST_PATH_IMAGE002
；其中，
Figure DEST_PATH_IMAGE111
表示磁分析的系数矩阵，包含N*N个元素，
Figure 311218DEST_PATH_IMAGE111
中每个元素用
Figure 76917DEST_PATH_IMAGE006
表示；H表示温度分布的系数矩阵，包含N*N个元素，H为中每个元素用
Figure 917834DEST_PATH_IMAGE008
表示；
Figure 254138DEST_PATH_IMAGE010
表示磁分析的已知向量，大小为N，
Figure 838703DEST_PATH_IMAGE010
中每个元素用
Figure 960374DEST_PATH_IMAGE012
表示；
Figure 441034DEST_PATH_IMAGE014
表示温度分布的已知向量，大小为N，
Figure 530212DEST_PATH_IMAGE014
中每个元素用
Figure 652889DEST_PATH_IMAGE016
表示；
Figure 143913DEST_PATH_IMAGE018
表示磁分析的未知向量，大小为N，
Figure 575901DEST_PATH_IMAGE018
中每个元素用
Figure 886796DEST_PATH_IMAGE020
表示；
Figure DEST_PATH_IMAGE112
表示温度分布的未知向量，大小为N，K中每个元素用
Figure 344322DEST_PATH_IMAGE024
表示；
Figure 440586DEST_PATH_IMAGE006
、
Figure 263048DEST_PATH_IMAGE008
、
Figure 61240DEST_PATH_IMAGE012
和
Figure 791298DEST_PATH_IMAGE016
通过以下方法计算：
Figure 991336DEST_PATH_IMAGE026
；
Figure 765125DEST_PATH_IMAGE028
；
Figure DEST_PATH_IMAGE113
；
Figure 581772DEST_PATH_IMAGE032
；
其中，
Figure 115521DEST_PATH_IMAGE034
表示数组B中
Figure 920797DEST_PATH_IMAGE036
行描述的空间区域，
Figure 616221DEST_PATH_IMAGE038
表示数组A中
Figure 389005DEST_PATH_IMAGE036
行数据建立的形状函数，
Figure 460866DEST_PATH_IMAGE040
表示数组A中
Figure 619184DEST_PATH_IMAGE042
行数据建立的形状函数，
Figure 751088DEST_PATH_IMAGE044
表示有载分接开关的磁边界条件，
Figure 745589DEST_PATH_IMAGE046
表示有载分接开关热边界条件。
进一步地，另一实施例中，磁-热耦合分布计算模块包括未知向量计算单元，通过以下方法计算：
根据节点初始向量C0和H0，初始向量C0和H0如下：
Figure DEST_PATH_IMAGE096A
将磁分析的系数矩阵D分解成M、B和N三个矩阵，通过以下方法分解：
如果i=j，则
Figure 434190DEST_PATH_IMAGE050
=
Figure 197747DEST_PATH_IMAGE052
，否则
Figure 500552DEST_PATH_IMAGE050
=0；
如果i>j，则
Figure 231617DEST_PATH_IMAGE054
=
Figure 645280DEST_PATH_IMAGE056
，否则
Figure 528923DEST_PATH_IMAGE054
=0；
如果i<j，则
Figure 2629DEST_PATH_IMAGE058
=
Figure 253613DEST_PATH_IMAGE056
，否则
Figure 470968DEST_PATH_IMAGE058
=0；
其中，i为自然数，且0<i<N；j为自然数，且0<j<N；
将温度分布的系数矩阵H分解成S、V和Q三个矩阵，通过以下方法分解：
如果i=j，则
Figure DEST_PATH_IMAGE114
=
Figure DEST_PATH_IMAGE115
，否则
Figure 989543DEST_PATH_IMAGE114
=0；
如果i>j，则
Figure DEST_PATH_IMAGE116
=
Figure DEST_PATH_IMAGE117
，否则
Figure 696468DEST_PATH_IMAGE116
=0；
如果i<j，则
Figure DEST_PATH_IMAGE118
=
Figure 434748DEST_PATH_IMAGE117
，否则
Figure 190214DEST_PATH_IMAGE118
=0；
其中，i为自然数，且0<i<N；j为自然数，且0<j<N；
得到从
Figure DEST_PATH_IMAGE119
、
Figure DEST_PATH_IMAGE120
、
Figure DEST_PATH_IMAGE121
、
Figure 891192DEST_PATH_IMAGE114
、
Figure 441122DEST_PATH_IMAGE116
和
Figure 135539DEST_PATH_IMAGE118
，计算
Figure DEST_PATH_IMAGE122
和
Figure DEST_PATH_IMAGE123
的公式：
Figure DEST_PATH_IMAGE124
（1-1）
Figure DEST_PATH_IMAGE125
（1-2）
Figure 549158DEST_PATH_IMAGE122
=
Figure DEST_PATH_IMAGE126
（1-3）
Figure 793058DEST_PATH_IMAGE123
=
Figure DEST_PATH_IMAGE127
--
Figure DEST_PATH_IMAGE128
（1-4）
其中，w为自然数，且0<12000，重复迭代计算上述（1-1）-（1-4）的公式，直到
Figure 592518DEST_PATH_IMAGE097
且
Figure DEST_PATH_IMAGE129
。
进一步地，另一实施例中，该系统还包括：等位云图生成模块，
通过以下公式计算有载分接开关三维模型中任一点的电场信息：
通过以下公式计算有载分接开关三维模型中任一点的电场信息：
Figure 757920DEST_PATH_IMAGE098
；
其中，
Figure DEST_PATH_IMAGE130
表示有载分接开关磁-热耦合分析中的某一点电位，
Figure DEST_PATH_IMAGE131
表示空间中j单元的体积,M为三维图形模型中的剖分的单元数，
Figure DEST_PATH_IMAGE132
和
Figure DEST_PATH_IMAGE133
是根据数组A中i行和j行数据建立的形状函数；
根据电场信息生成有载分接开关的电场云图。
本申请实施例还公开了一种计算机设备，该计算机设备可以是服务器。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。该计算机程序被处理器执行时以实现一种换流变有载分接开关切换芯子磁-热耦合计算方法，该方法包括以下步骤：
S1、获取有载分接开关各个部件的尺寸信息，并获取有载分接开关各个部件对应的磁导率；
S2、基于所获取的尺寸信息和磁导率建立有载分接开关的三维模型；
S3、基于三维模型确定有载分接开关磁-热耦合分析的线性方程组，并从线性方程组确定有载分接开关各个部件的磁-热耦合分布。
本申请实施例还公开一种计算机可读存储介质。在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
S1、获取有载分接开关各个部件的尺寸信息，并获取有载分接开关各个部件对应的磁导率；
S2、基于所获取的尺寸信息和磁导率建立有载分接开关的三维模型；
S3、基于三维模型确定有载分接开关磁-热耦合分析的线性方程组，并从线性方程组确定有载分接开关各个部件的磁-热耦合分布。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围，均应包含在本发明的保护范围之内。

Claims (9)

1.一种换流变有载分接开关切换芯子磁-热耦合计算方法，其特征在于，包括以下步骤：
获取有载分接开关各个部件的尺寸信息，并获取有载分接开关各个部件对应的磁导率；
基于所获取的尺寸信息和磁导率建立有载分接开关的三维模型；
基于所述三维模型确定有载分接开关磁-热耦合分析的线性方程组，并从线性方程组确定有载分接开关各个部件的磁-热耦合分布。
2.根据权利要求1所述的换流变有载分接开关切换芯子磁-热耦合计算方法，其特征在于，所述线性方程组通过以下方法计算：
对所述三维模型进行自适应网格剖分，以获取有载分接开关磁-热耦合分析的节点数N和单元数M；
根据所述节点数N和单元数M确定所述线性方程组：
Figure DEST_PATH_IMAGE001
；其中，
Figure 798181DEST_PATH_IMAGE002
表示磁分析的系数矩阵，包含N*N个元素，
Figure 552642DEST_PATH_IMAGE002
中每个元素用
Figure DEST_PATH_IMAGE003
表示；H表示温度分布的系数矩阵，包含N*N个元素，H为中每个元素用
Figure 728408DEST_PATH_IMAGE004
表示；
Figure DEST_PATH_IMAGE005
表示磁分析的已知向量，大小为N，
Figure 965223DEST_PATH_IMAGE005
中每个元素用
Figure 720690DEST_PATH_IMAGE006
表示；
Figure DEST_PATH_IMAGE007
表示温度分布的已知向量，大小为N，
Figure 844504DEST_PATH_IMAGE007
中每个元素用
Figure 410745DEST_PATH_IMAGE008
表示；
Figure DEST_PATH_IMAGE009
表示磁分析的未知向量，大小为N，
Figure 885589DEST_PATH_IMAGE009
中每个元素用
Figure 444746DEST_PATH_IMAGE010
表示；
Figure DEST_PATH_IMAGE011
表示温度分布的未知向量，大小为N，K中每个元素用
Figure 469072DEST_PATH_IMAGE012
表示；
Figure 658745DEST_PATH_IMAGE003
、
Figure 371617DEST_PATH_IMAGE004
、
Figure 468886DEST_PATH_IMAGE006
和
Figure 301713DEST_PATH_IMAGE008
通过以下方法计算：
Figure DEST_PATH_IMAGE013
；
Figure 973871DEST_PATH_IMAGE014
；
Figure DEST_PATH_IMAGE015
；
Figure 423307DEST_PATH_IMAGE016
；
其中，
Figure DEST_PATH_IMAGE017
表示数组B中
Figure 606158DEST_PATH_IMAGE018
行描述的空间区域，
Figure DEST_PATH_IMAGE019
表示数组A中
Figure 559071DEST_PATH_IMAGE018
行数据建立的形状函数，
Figure 887284DEST_PATH_IMAGE020
表示数组A中
Figure 276546DEST_PATH_IMAGE022
行数据建立的形状函数，
Figure DEST_PATH_IMAGE023
表示有载分接开关的磁边界条件，
Figure 512355DEST_PATH_IMAGE024
表示有载分接开关热边界条件。
3.根据权利要求2所述的换流变有载分接开关切换芯子磁-热耦合计算方法，其特征在于，所述线性方程组中的未知向量通过以下方法计算：
根据所述节点初始向量C0和H0，
将磁分析的系数矩阵D分解成M、B和N三个矩阵，通过以下方法分解：
如果i=j，则
Figure DEST_PATH_IMAGE025
=
Figure 336085DEST_PATH_IMAGE026
，否则
Figure 835200DEST_PATH_IMAGE025
=0；
如果i>j，则
Figure DEST_PATH_IMAGE027
=
Figure 993649DEST_PATH_IMAGE028
，否则
Figure 497398DEST_PATH_IMAGE027
=0；
如果i<j，则
Figure DEST_PATH_IMAGE029
=
Figure 424902DEST_PATH_IMAGE028
，否则
Figure 94918DEST_PATH_IMAGE029
=0；
其中，i为自然数，且0<i<N；j为自然数，且0<j<N；
将温度分布的系数矩阵H分解成S、V和Q三个矩阵，通过以下方法分解：
如果i=j，则
Figure 960237DEST_PATH_IMAGE030
=
Figure DEST_PATH_IMAGE031
，否则
Figure 272270DEST_PATH_IMAGE030
=0；
如果i>j，则
Figure 788702DEST_PATH_IMAGE032
=
Figure DEST_PATH_IMAGE033
，否则
Figure 941203DEST_PATH_IMAGE032
=0；
如果i<j，则
Figure 543086DEST_PATH_IMAGE034
=
Figure 878383DEST_PATH_IMAGE033
，否则
Figure 249322DEST_PATH_IMAGE034
=0；
其中，i为自然数，且0<i<N；j为自然数，且0<j<N；
得到从
Figure 261140DEST_PATH_IMAGE025
、
Figure 84740DEST_PATH_IMAGE027
、
Figure 722263DEST_PATH_IMAGE029
、
Figure 947708DEST_PATH_IMAGE030
、
Figure 130428DEST_PATH_IMAGE032
和
Figure 706903DEST_PATH_IMAGE034
，计算
Figure DEST_PATH_IMAGE035
和
Figure 446320DEST_PATH_IMAGE036
的公式：
Figure DEST_PATH_IMAGE037
（1-1）
Figure 306697DEST_PATH_IMAGE038
（1-2）
Figure 925897DEST_PATH_IMAGE035
=
Figure DEST_PATH_IMAGE039
（1-3）
Figure 537139DEST_PATH_IMAGE036
=
Figure 267197DEST_PATH_IMAGE040
--
Figure DEST_PATH_IMAGE041
（1-4）
其中，w为自然数，且0<12000，重复迭代计算上述（1-1）-（1-4）的公式，直到
Figure 732814DEST_PATH_IMAGE042
且
Figure DEST_PATH_IMAGE043
。
4.根据权利要求3所述的换流变有载分接开关切换芯子磁-热耦合计算方法，其特征在于，所述方法还包括：
通过以下公式计算有载分接开关三维模型中任一点的电场信息：
Figure 303341DEST_PATH_IMAGE044
；
其中，
Figure DEST_PATH_IMAGE045
表示有载分接开关磁-热耦合分析中的某一点电位，
Figure 119987DEST_PATH_IMAGE046
表示空间中j单元的体积,M为三维图形模型中的剖分的单元数，
Figure 138890DEST_PATH_IMAGE019
和
Figure 193434DEST_PATH_IMAGE020
是根据所述数组A中i行和j行数据建立的形状函数；
根据所述电场信息生成有载分接开关的电场云图。
5.一种换流变有载分接开关切换芯子磁-热耦合计算系统，其特征在于，所述系统包括参数获取模块、三维建模模块和磁-热耦合分布计算模块：
所述参数获取模块用于获取有载分接开关各个部件的尺寸信息，并获取有载分接开关各个部件对应的磁导率；
所述三维建模模块与参数获取模块连接，用于接收所述尺寸信息和磁导率，基于所述尺寸信息和磁导率建立有载分接开关的三维模型；
所述磁-热耦合分布计算模块与三维建模模块连接，用于接收所述三维模型，基于所述三维模型确定有载分接开关磁-热耦合分析的线性方程组，并从所述线性方程组确定有载分接开关各个部件的磁-热耦合分布。
6.一种应用如权利要求1-4任一所述的一种换流变有载分接开关切换芯子磁-热耦合计算方法，其特征在于，所述方法用于±500kV直流输电工程整流侧或者逆变侧换流变有载分接开关磁热分析确定。
7.根据权利要求6所述的换流变有载分接开关切换芯子磁-热耦合计算方法，其特征在于，所述方法用于±800kV直流输电工程整流侧或者逆变侧换流变有载分接开关磁热分析确定。
8.一种电子设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，其特征在于：所述处理器执行所述计算机程序时实现如权利要求1-4任一项所述的一种换流变有载分接开关切换芯子磁-热耦合计算系统方法的步骤。
9.一种计算机可读存储介质，其特征在于，存储有能够被处理器加载并执行如权利要求1-4中任一种方法的计算机程序。
CN202111236750.5A 2021-10-23 2021-10-23 换流变有载分接开关切换芯子磁-热耦合计算方法及系统 Pending CN114238850A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111236750.5A CN114238850A (zh) 2021-10-23 2021-10-23 换流变有载分接开关切换芯子磁-热耦合计算方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111236750.5A CN114238850A (zh) 2021-10-23 2021-10-23 换流变有载分接开关切换芯子磁-热耦合计算方法及系统

Publications (1)

Publication Number Publication Date
CN114238850A true CN114238850A (zh) 2022-03-25

Family

ID=80743224

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111236750.5A Pending CN114238850A (zh) 2021-10-23 2021-10-23 换流变有载分接开关切换芯子磁-热耦合计算方法及系统

Country Status (1)

Country Link
CN (1) CN114238850A (zh)

Similar Documents

Publication Publication Date Title
CN107958125B (zh) 一种三相变压器电磁暂态中低频模型建模方法及系统
Kumbhar et al. Reduction of loss and local overheating in the tank of a current transformer
CN115561564B (zh) 一种电缆接头动态载流量的arima序列预测方法
CN107103161B (zh) 一种累积效应对变压器抗短路能力影响的仿真计算方法
Acero et al. Magnetic vector potential based model for eddy-current loss calculation in round-wire planar windings
CN109543246B (zh) 一种确定变压器直流偏磁温度场分布的方法和系统
Girgis et al. Calculation of core hot-spot temperature in power and distribution transformers
Dey et al. Improving thermal performance of high frequency power transformers using bobbinless transformer design
CN114238850A (zh) 换流变有载分接开关切换芯子磁-热耦合计算方法及系统
Urooj et al. Optimal design of power transformer with advance core material using ANSYS technique
CN107274088B (zh) 一种地下配电室的风险场计算方法及系统
Valdivieso et al. Time-domain homogenization of foil windings in 2-D axisymmetric finite-element models
Demirdelen A heuristic approach for volume calculation of oil-type power transformers: firefly algorithm
Tahir et al. Performance evaluation of a tape-wound core transformer using meta-model based scaling laws
Zhang et al. 2-D coupled fluid-thermal analysis of oil-immersed power transformers based on finite element method
Ahmad et al. Short circuit stress calculation in power transformer using finite element method on high voltage winding displaced vertically
Cano-Plata et al. Hybrid simulation of power quality assessment: an application for power ground grid in arc furnace systems
Naderian et al. A fast method for calculation of transformers leakage reactance using energy technique
Naderian-Jahromi et al. Calculation of distribution transformer leakage reactance using energy technique
Díaz Flórez et al. Novel technique for the calculation of eddy current losses and Lorentz forces in foil winding transformers
CN110910968A (zh) 一种高压电气设备放电分解产物空间分布情况的计算方法
Dursun et al. Harmonic load losses in power transformer windings using Finite Element methods
Shazly et al. Extensions to the Finite Element Technique for the magneto-thermal analysis of aged oil cooled-insulated power transformers
CN110850172B (zh) 一种调匝式消弧线圈直流电阻合格程度的评价方法
CN109241586A (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