变压器铁心励磁特性的有限元计算方法、装置及终端与流程

1.本发明涉及变压器技术领域，尤其涉及一种变压器铁心励磁特性的有限元计算方法、装置及终端。


背景技术：

2.随着我国电力网络的不断成长和发展，风能和太阳能等新能源发电形式的分布式电源越来越多的并网运行，大量电力电子器件投入应用，例如，大功率整流器件、变频器、工业电弧炉负载等，这些电力电子器件的大量应用会带来严重的间谐波电能污染，导致电力系统波形越来越复杂。
3.含铁心绕组类变压器作为电力系统中的重要设备，在含间谐波电压的影响下会发生绕组电流畸变，导致铁心达到过饱和状态从而产生励磁电流畸变，引起损耗增加，还可能会产生噪声问题和振动问题，影响变压器的正常工作。因此，需要进行含间谐波电压激励的频率、含量变化对铁心励磁特性影响分析及计算；
4.现有技术中，求取励磁特性往往通过多物理场耦合的方式，计算过程复杂、耗时长且准确度不高，尤其是对于不同激励下的励磁情况，往往需要对不同的激励条件分别重新计算。


技术实现要素：

5.本发明提供一种变压器铁心励磁特性的有限元计算方法、装置及终端，以省时、高效的实现含铁心绕组类变压器励磁特性的计算。
6.第一方面，本发明提供了一种变压器铁心励磁特性的有限元计算方法，包括：
7.建立变压器的结构模型；
8.对变压器的结构模型进行网格剖分，得到三维有限元磁场模型；
9.基于三维有限元磁场模型，通过电路计算获得变压器绕组的瞬态电流值和动态电感值；
10.基于三维有限元磁场模型，通过磁场计算获得变压器铁心的励磁电流值和磁通密度分布。
11.在一种可能的实现方式中，在获得变压器绕组的瞬态电流值和动态电感值，以及，变压器铁心的励磁电流值和磁通密度分布之后，方法还包括；
12.基于变压器绕组的瞬态电流值和变压器铁心的磁通密度分布，利用变压器电磁耦合特性建立电磁耦合模型；
13.利用变压器绕组的动态电感值、变压器铁心的励磁电流值为参数进行循环迭代，实现变压器磁场和电路的间接耦合。
14.在一种可能的实现方式中，建立变压器的结构模型，包括；
15.获取窗口参数，窗口参数包括铁心窗口高度和铁心窗口宽度；
16.获取铁心参数，铁心参数包括铁心叠厚、铁心宽度和铁心高度；
17.获取绕组参数，绕组参数包括绕组匝数和绕组尺寸；
18.获取电场参数，电场参数包括铁心密度、铁心相对介电常数、铁心电导率、绕组相对介电常数、绕组电导率、空气相对介电常数和空气电导率；
19.获取磁场参数，磁场参数包括铁心沿轧制方向的b-h磁化曲线、铁心沿垂直轧制方向的b-h磁化曲线、绕组相对磁导率和空气相对磁导率；
20.基于窗口参数进行结构模型的窗口设置，基于铁心参数进行结构模型的铁心参数设置，基于绕组参数进行结构模型的绕组参数设置，基于电场参数进行结构模型的电场参数设置，基于磁场参数进行结构模型的磁场参数设置，得到变压器的结构模型。
21.在一种可能的实现方式中，对变压器的结构模型进行网格剖分，包括；
22.对铁心部分采用四面体网格进行剖分，对绕组线圈部分采用扫掠式网格进行剖分，对空气域部分采用四面体网格进行剖分；
23.其中，对铁心部分进行剖分的四面体网格细于对空气域部分进行剖分的四面体网格。
24.在一种可能的实现方式中，基于三维有限元磁场模型，通过电路计算获得变压器绕组的瞬态电流值和动态电感值，包括：
25.基于三维有限元磁场模型和预设的变压器电路微分方程求解获得绕组的瞬态电流值和动态电感值；
26.变压器电路微分方程为：
[0027][0028]
r为电路的电阻矩阵，v(t)为电源电压瞬时值，ld为电路的动态电感矩阵，i为绕组的瞬态电流值，t为时间。
[0029]
在一种可能的实现方式中，基于三维有限元磁场模型，通过磁场计算获得变压器铁心的励磁电流值和磁通密度分布，包括：
[0030]
将瞬态电流值代入麦克斯韦方程计算出变压器铁心的磁通密度分布及励磁电流值；
[0031]
麦克斯韦方程为：
[0032][0033]
式中：b为磁通密度，j为电流密度，h为磁场强度，e为电场强度，t为时间。
[0034]
第二方面，本发明提供了一种变压器铁心励磁特性的有限元计算装置，包括：
[0035]
模型建立模块，用于建立变压器的结构模型；
[0036]
网格剖分模块，用于对变压器的结构模型进行网格剖分，得到三维有限元磁场模型；
[0037]
数据获取模块，用于基于三维有限元磁场模型，通过电路计算获得变压器绕组的瞬态电流值和动态电感值；
[0038]
数据运算模块，用于基于三维有限元磁场模型，通过磁场计算获得变压器铁心的励磁电流值和磁通密度分布。
[0039]
在一种可能的实现方式中，所述模型建立模块，还用于：
[0040]
获取窗口参数，窗口参数包括铁心窗口高度和铁心窗口宽度；
[0041]
获取铁心参数，铁心参数包括铁心叠厚、铁心宽度和铁心高度；
[0042]
获取绕组参数，绕组参数包括绕组匝数和绕组尺寸；
[0043]
获取电场参数，电场参数包括铁心密度、铁心相对介电常数、铁心电导率、绕组相对介电常数、绕组电导率、空气相对介电常数和空气电导率；
[0044]
获取磁场参数，磁场参数包括铁心沿轧制方向的b-h磁化曲线、铁心沿垂直轧制方向的b-h磁化曲线、绕组相对磁导率和空气相对磁导率；
[0045]
基于窗口参数进行结构模型的窗口设置，基于铁心参数进行结构模型的铁心参数设置，基于绕组参数进行结构模型的绕组参数设置，基于电场参数进行结构模型的电场参数设置，基于磁场参数进行结构模型的磁场参数设置，得到变压器的结构模型。
[0046]
第三方面，本发明提供了一种终端，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上第一方面或第一方面的任一种可能的实现方式所述方法的步骤。
[0047]
第四方面，本发明提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如上第一方面或第一方面的任一种可能的实现方式所述方法的步骤。
[0048]
本发明提供一种变压器铁心励磁特性的有限元计算方法、装置及终端，首先，建立变压器的结构模型，然后，对变压器的结构模型进行网格剖分，得到三维有限元磁场模型；之后，基于三维有限元磁场模型，通过电路计算获得变压器绕组的瞬态电流值和动态电感值；最后，基于三维有限元磁场模型，通过磁场计算获得变压器铁心的励磁电流值和磁通密度分布。本发明通过建立变压器的结构模型并进行网格剖分得到三维有限元磁场模型，可以更准确的反映含铁心绕组类变压器的特性，通过电路计算和磁场计算，可以更快速高效的反映变压器铁心在不同激励条件下的励磁特性变化。
附图说明
[0049]
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
[0050]
图1是本发明实施例提供的变压器铁心励磁特性的有限元计算方法的实现流程图；
[0051]
图2是本发明实施例提供的变压器的结构模型的示意图；
[0052]
图3是本发明实施例提供的b-h磁化曲线的示意图；
[0053]
图4是本发明实施例提供的三维有限元磁场模型的示意图；
[0054]
图5是本发明实施例提供的三维有限元磁场模型的磁通密度分布的示意图；
[0055]
图6是本发明实施例提供的附加在1.575v/5hz含间谐波电压时励磁电流的计算值和理论值的示意图；
[0056]
图7是本发明实施例提供的附加在1.575v/10hz含间谐波电压时励磁电流的计算
值和理论值的示意图；
[0057]
图8是本发明实施例提供的变压器铁心励磁特性的有限元计算装置的结构示意图；
[0058]
图9是本发明实施例提供的终端的示意图。
具体实施方式
[0059]
以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本发明实施例。然而，本领域的技术人员应当清楚，在没有这些具体细节的其它实施例中也可以实现本发明。在其它情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本发明的描述。
[0060]
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图通过具体实施例来进行说明。
[0061]
参见图1，其示出了本发明实施例提供的变压器铁心励磁特性的有限元计算方法的实现流程图，详述如下：
[0062]
步骤s110、建立变压器的结构模型。
[0063]
本实施例中，使用comsol仿真软件根据变压器铁心的窗口参数、变压器的铁心参数、变压器铁心的绕组参数、变压器的电场参数和变压器的磁场参数，建立变压器的结构模型，如图2所示。
[0064]
具体为，获取窗口参数，窗口参数包括铁心窗口高度和铁心窗口宽度；示例性的，铁心窗口高度是130mm，铁心窗口宽度是70mm。
[0065]
获取铁心参数，铁心参数包括铁心叠厚、铁心宽度和铁心高度；示例性的，铁心叠厚是66mm，铁心宽度是200mm，铁心高度是260mm。
[0066]
获取绕组参数，绕组参数包括绕组匝数和绕组尺寸；示例性的，绕组匝数是20匝，绕组尺寸是高60mm厚4mm的方形线圈结构。
[0067]
获取电场参数，电场参数包括铁心密度、铁心相对介电常数、铁心电导率、绕组相对介电常数、绕组电导率、空气相对介电常数和空气电导率；示例性的，铁心密度为7650kg/m3,铁心相对介电常数是1，铁心电导率是5s/m，绕组相对介电常数是1，绕组电导率5s/m，空气相对介电常数为1，空气电导率是5s/m。
[0068]
获取磁场参数，磁场参数包括铁心沿轧制方向的b-h磁化曲线、铁心沿垂直轧制方向的b-h磁化曲线、绕组相对磁导率和空气相对磁导率；示例性的，铁心沿轧制方向的b-h磁化曲线、铁心沿垂直轧制方向的b-h磁化曲线如图3所示，绕组相对磁导率是5s/m，空气相对磁导率是5s/m。
[0069]
基于窗口参数进行结构模型的窗口设置，基于铁心参数进行结构模型的铁心参数设置，基于绕组参数进行结构模型的绕组参数设置，基于电场参数进行结构模型的电场参数设置，基于磁场参数进行结构模型的磁场参数设置，得到变压器的结构模型。
[0070]
步骤s120、对变压器的结构模型进行网格剖分，得到三维有限元磁场模型。
[0071]
本实施例中，使用comsol仿真软件对变压器的结构模型进行网格剖分，得到三维有限元磁场模型，如图4所示。
[0072]
具体的，结构模型包括铁心部分、绕组线圈部分和空气域部分。
[0073]
具体的，对铁心部分采用四面体网格进行剖分，对绕组线圈部分采用扫掠式网格进行剖分，对空气域部分采用四面体网格进行剖分，得到三维有限元磁场模型。
[0074]
具体的，对铁心部分进行剖分的四面体网格细于对空气域部分进行剖分的四面体网格。
[0075]
示例性的，对变压器的结构模型进行网格剖分，完整网格包含69432个域单元、8732个边界元和887个边单元。
[0076]
步骤s130、基于三维有限元磁场模型，通过电路计算获得变压器绕组的瞬态电流值和动态电感值。
[0077]
本实施例中，使用comsol仿真软件，根据三维有限元磁场模型和预设的变压器电路微分方程，通过电路仿真计算获得变压器绕组的瞬态电流值和动态电感值。
[0078]
具体的，基于三维有限元磁场模型和预设的变压器电路微分方程求解获得绕组的瞬态电流值和动态电感值；
[0079]
变压器电路微分方程为：
[0080][0081]
r为电路的电阻矩阵，v(t)为电源电压瞬时值，ld为电路的动态电感矩阵，i为绕组的瞬态电流值，t为时间。
[0082]
示例性的，为了提高计算精确度，设置pardiso作为求解器，选择直接全耦合的计算模式，设置求解步长0.2ms和求解时间2000ms。其中，求解步长越小，求解精度越高。
[0083]
步骤s140、基于三维有限元磁场模型，通过磁场计算获得变压器铁心的励磁电流值和磁通密度分布。
[0084]
本实施例中，使用comsol仿真软件，根据瞬态电流值和麦克斯韦方程通过磁场仿真的方式计算变压器铁心的励磁电流值和磁通密度分布。
[0085]
具体为，将瞬态电流值代入麦克斯韦方程计算出变压器铁心的磁通密度分布及励磁电流值；
[0086]
麦克斯韦方程为：
[0087][0088]
式中：b为磁通密度，j为电流密度，h为磁场强度，e为电场强度，t为时间，
[0089]
将定义矢量磁位a，令得到变压器铁心电磁场微分方程，计算出变压器铁心的磁通密度分布及励磁电流值；
[0090]
变压器铁心电磁场微分方程为：
[0091][0092]
式中：jc线圈电流密度；μ为磁阻率；η为电导率。
[0093]
可选的，在获得变压器绕组的瞬态电流值和动态电感值，以及，变压器铁心的励磁电流值和磁通密度分布之后，
[0094]
基于变压器绕组的瞬态电流值和变压器铁心的磁通密度分布，利用变压器电磁耦合特性建立电磁耦合模型；
[0095]
利用变压器绕组的动态电感值、变压器铁心的励磁电流值为参数进行循环迭代，实现变压器磁场和电路的间接耦合。
[0096]
在一种可能的实现方式中，在绕组施加频率为50hz，幅值为52.5v的正弦电压，可以得到三维有限元磁场模型的磁通密度分布如图5所示，铁心的磁通密度分布为1.3t，和铁心的理论值基本相等，可证明三维有限元模型可以准确计算出铁心磁通密度分布。
[0097]
在一种可能的实现方式中，在绕组施加频率为50hz，幅值为52.5v分别附加在：0.525v(基频1％)/5hz含间谐波电压、1.575v(基频3％)/5hz含间谐波电压、1.35v(基频3％)/10hz含间谐波电压、1.35v(基频3％)/7hz含间谐波电压、2.625v(基频5％)/5hz含间谐波电压的激励。基于三维有限元磁场模型计算磁通密度分布，对计算值与理论值进行比对，得到误差情况如表1，可证明本模型可以准确反映铁心在含间谐波电压的磁通密度变化。
[0098][0099]
表1
[0100]
在一种可能的实现方式中,在绕组施加频率为50hz，幅值为52.5v，附加在1.575v/5hz含间谐波电压时，基于三维有限元磁密模型计算励磁电流，得到励磁电流的计算值和理论值如图6所示；在绕组施加频率为50hz，幅值为52.5v，附加在1.575v/10hz含间谐波电压时，基于三维有限元磁密模型计算励磁电流，得到励磁电流的计算值和理论值如图7所示，可证明三维有限元磁场模型可以准确反映铁心在含间谐波电压下励磁电流的变化。
[0101]
在本发明实施例中，首先，建立变压器的结构模型。之后，对变压器的结构模型进行网格剖分，得到三维有限元磁场模型。然后，基于三维有限元磁场模型，通过电路计算获得变压器绕组的瞬态电流值和动态电感值。最后，基于三维有限元磁场模型，通过磁场计算获得变压器铁心的励磁电流值和磁通密度分布。如此，通过磁通密度分布、动态电感值、瞬态电流值和励磁电流值可以反映出含铁心绕组类变压器在含间谐波电压影响下发生的变化。
[0102]
应理解，上述实施例中各步骤的序号的大小并不意味着执行顺序的先后，各过程的执行顺序应以其功能和内在逻辑确定，而不应对本发明实施例的实施过程构成任何限定。
[0103]
基于上述实施例提供的变压器铁心励磁特性的有限元计算方法，相应地，本发明还提供了应用于变压器铁心励磁特性的有限元计算方法的变压器铁心励磁特性的有限元计算方法装置的具体实现方式。请参见以下实施例。
[0104]
图8示出了本发明实施例提供的变压器铁心励磁特性的有限元计算装置的结构示意图，为了便于说明，仅示出了与本发明实施例相关的部分，详述如下：
[0105]
如图8所示，变压器铁心励磁特性的有限元计算装置800包括：模型建立模块810、网格剖分模块820、数据获取模块830、数据运算模块840。
[0106]
模型建立模块810，用于建立变压器的结构模型。
[0107]
网格剖分模块820，用于对变压器的结构模型进行网格剖分，得到三维有限元磁场模型。
[0108]
数据获取模块830，用于基于三维有限元磁场模型，通过电路计算获得变压器绕组的瞬态电流值和动态电感值。
[0109]
数据运算模块840，用于基于三维有限元磁场模型，通过磁场计算获得变压器铁心的励磁电流值和磁通密度分布。
[0110]
在一种可能的实现方式中，所述模型建立模块810，具体用于：
[0111]
获取窗口参数，窗口参数包括铁心窗口高度和铁心窗口宽度；
[0112]
获取铁心参数，铁心参数包括铁心叠厚、铁心宽度和铁心高度；
[0113]
获取绕组参数，绕组参数包括绕组匝数和绕组尺寸；
[0114]
获取电场参数，电场参数包括铁心密度、铁心相对介电常数、铁心电导率、绕组相对介电常数、绕组电导率、空气相对介电常数和空气电导率；
[0115]
获取磁场参数，磁场参数包括铁心沿轧制方向的b-h磁化曲线、铁心沿垂直轧制方向的b-h磁化曲线、绕组相对磁导率和空气相对磁导率；
[0116]
基于窗口参数进行结构模型的窗口设置，基于铁心参数进行结构模型的铁心参数设置，基于绕组参数进行结构模型的绕组参数设置，基于电场参数进行结构模型的电场参数设置，基于磁场参数进行结构模型的磁场参数设置，得到变压器的结构模型。
[0117]
在一种可能的实现方式中，所述网格剖分模块820，具体用于：
[0118]
结构模型包括铁心部分、绕组线圈部分和空气域部分；
[0119]
对变压器的结构模型进行网格剖分包括：
[0120]
对铁心部分采用四面体网格进行剖分，对绕组线圈部分采用扫掠式网格进行剖分，对空气域部分采用四面体网格进行剖分；
[0121]
其中，对铁心部分进行剖分的四面体网格细于对空气域部分进行剖分的四面体网格。
[0122]
在一种可能的实现方式中，所述数据获取模块830，具体用于：
[0123]
基于三维有限元磁场模型和预设的变压器电路微分方程求解获得瞬态电流值和动态电感值；
[0124]
变压器电路微分方程为：
[0125][0126]
r为电路的电阻矩阵，v(t)为电源电压瞬时值，ld为电路的动态电感矩阵，i为绕组的瞬态电流值，t为时间。
[0127]
在一种可能的实现方式中，所述数据运算模块840，具体用于：
[0128]
将瞬态电流值代入麦克斯韦方程计算出变压器铁心的磁通密度分布及励磁电流值；
[0129]
麦克斯韦方程为：
[0130][0131]
式中：b为磁通密度，j为电流密度，h为磁场强度，e为电场强度，t为时间。
[0132]
在一种可能的实现方式中，所述数据运算模块840，还具体用于：
[0133]
在获取变压器绕组的瞬态电流值和动态电感值，以及，变压器铁心的励磁电流值和磁通密度分布之后；
[0134]
基于变压器绕组的瞬态电流值和变压器铁心的磁通密度分布，利用变压器电磁耦合特性建立电磁耦合模型；
[0135]
利用变压器绕组的动态电感值、变压器铁心的励磁电流值为参数进行循环迭代，实现变压器磁场和电路的间接耦合。
[0136]
图9是本发明实施例提供的终端的示意图。如图9所示，该实施例的终端9包括：处理器90、存储器91以及存储在所述存储器91中并可在所述处理器90上运行的计算机程序92。所述处理器90执行所述计算机程序92时实现上述各个变压器铁心励磁特性的有限元计算方法实施例中的步骤，例如图1所示的步骤s110至步骤s140。或者，所述处理器90执行所述计算机程序92时实现上述各装置实施例中各模块/单元的功能，例如图8所示模块810至840的功能。
[0137]
示例性的，所述计算机程序92可以被分割成一个或多个模块/单元，所述一个或者多个模块/单元被存储在所述存储器91中，并由所述处理器90执行，以完成本发明。所述一个或多个模块/单元可以是能够完成特定功能的一系列计算机程序指令段，该指令段用于描述所述计算机程序92在所述终端9中的执行过程。例如，所述计算机程序92可以被分割成图8所示的模块810至840。
[0138]
所述终端9可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。所述终端9可包括，但不仅限于，处理器90、存储器91。本领域技术人员可以理解，图9仅仅是终端9的示例，并不构成对终端9的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如所述终端还可以包括输入输出设备、网络接入设备、总线等。
[0139]
所称处理器90可以是中央处理单元(central processing unit，cpu)，还可以是其他通用处理器、数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit，asic)、现场可编程门阵列(field-programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
[0140]
所述存储器91可以是所述终端9的内部存储单元，例如终端9的硬盘或内存。所述存储器91也可以是所述终端9的外部存储设备，例如所述终端9上配备的插接式硬盘，智能存储卡(smart media card，smc)，安全数字(secure digital，sd)卡，闪存卡(flash card)等。进一步地，所述存储器91还可以既包括所述终端9的内部存储单元也包括外部存储设备。所述存储器91用于存储所述计算机程序以及所述终端所需的其他程序和数据。所述存储器91还可以用于暂时地存储已经输出或者将要输出的数据。
[0141]
所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，仅以上述各功
能单元、模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能单元、模块完成，即将所述装置的内部结构划分成不同的功能单元或模块，以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中，上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。另外，各功能单元、模块的具体名称也只是为了便于相互区分，并不用于限制本技术的保护范围。上述系统中单元、模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
[0142]
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。
[0143]
本领域普通技术人员可以意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
[0144]
在本发明所提供的实施例中，应该理解到，所揭露的装置/终端和方法，可以通过其它的方式实现。例如，以上所描述的装置/终端实施例仅仅是示意性的，例如，所述模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口，装置或单元的间接耦合或通讯连接，可以是电性，机械或其它的形式。
[0145]
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
[0146]
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
[0147]
所述集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个变压器铁心励磁特性的有限元计算方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、u盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(read-only memory，rom)、随机存取存储器(random access memory，ram)、电载波信号、电信信号以及软件分发介质等。需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括是电载波信号和电信信号。
[0148]
以上所述实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围，均应包含在本发明的保护范围之内。