一种饱和条件下无取向硅钢片铁心磁导率张量表示方法
【专利摘要】本发明提供一种饱和条件下无取向硅钢片铁心磁导率张量表示方法，首先根据计算模型、通过迭代计算确定等效平均磁导率具体数值；然后计及叠片结构所引入的磁特性各向异性，建立磁导率张量表达式。由于计及了叠片铁心各向异性磁特性并在张量磁导率中引入计及铁磁材料饱和效应的等效平均磁导率分量，所建立的磁导率张量模型符合无取向硅钢片铁心实际磁化过程，并将实际非线性电磁场方程简化为线性方程，从而一定程度上削减计算所需时间和人力成本，为电机和变压器等电气设备的性能预测和优化设计高效实施提供了一种有效方案。
【专利说明】
-种饱和条件下无取向括钢片铁心磁导率张量表示方法
技术领域
[0001] 本发明设及一种娃钢片铁屯、磁导率表示方法，尤其设及一种饱和条件下无取向娃 钢片铁屯、磁导率张量表示方法。
【背景技术】
[0002] 无取向娃钢片铁屯、是各种类型电机、小型变压器等电气设备的基本部件，并通过 叠装或卷装工艺形成叠层铁屯、，主要用W构建磁路和安装骨架，并在各种电气设备中所占 体积和重量比例最大。对电机和变压器等电气设备进行性能预测和优化设计时，需要借助 有效的叠片铁屯、电磁特性数学模型尤其是磁导率模型，在此基础上进行电磁场及禪合场计 算。在对上述电气设备进行物理场分析过程中，由于计算资源的限制，无法将所有的娃钢片 独立建模和分析，因此面对的直接分析对象只能是铁屯、整体，并通常将其等效为电磁特性 各向异性的连续体。
[0003] 另一方面，在叠片铁屯、中很难避免出现局部磁饱和现象，尤其是在过载等条件下。 铁磁材料中的饱和效应通常需要借助相关材料的磁化曲线进行模拟，由此导致所需求解的 电磁场控制方程为非线性方程。电磁场方程的求解方案主要分为有限元法为代表的数值方 法和W分离变量法为代表的解析方法。在实际应用过程中，非线性数值分析的受制约因素 很多，特别是在铁屯、深度饱和情况下，非线性迭代计算的收敛性、计算时间W及计算过程稳 定性等问题尤为突出。而若选择使用解析法进行场分析，除非特殊情况，一般很难实现对非 线性控制方程的有效求解。
[0004] 综上所述，有必要通过合理表述铁磁材料中的饱和效应，实现对实际非线性电磁 场问题进行线性化处理，从而达到节省数值分析计算代价的基本目标，同时使得采用解析 法分析电磁场饱和问题成为可能。
【发明内容】
[0005] 本发明的目的是提供一种无取向娃钢片叠片铁屯、磁导率张量表示方法，计及了由 于叠片结构所引入的铁屯、磁属性各向异性特性，并通过在张量磁导率中引入计及铁磁材料 饱和效应的等效平均磁导率分量，将非线性控制方程简化为线性方程，从而一定程度上削 减电磁场计算所需时间和人力成本，为电机和变压器等电气设备的性能预测和优化设计高 效实施提供了一种有效方案。
[0006] 本发明的技术解决方案是：
[0007] -种饱和条件下无取向娃钢片叠片铁屯、磁导率张量表示方法，包括如下步骤：
[000引S1、借助相关计算模型、通过迭代计算确定等效平均磁导率Weq具体数值；
[0009] S2、计及叠片结构所引入的磁特性各向异性，若假设娃钢片厚度方向即叠片方向 沿直角坐标系的y轴方向，则所述磁导率张量基本形式为：
[0010]
[001。 上式中，F为铁屯、叠片系数，μ日为真空磁导率。
[0012]进一步地，步骤S1中，通过相关计算模型及迭代计算求取等效平均磁导率Weq，具体 步骤为：
[OOU] 811、设定饱和磁感应强度Bs的第k次尝试值皆)；
[0014] S12、根据娃钢片材料直流磁化曲线及关系式5=哉，^=。"&确定对应娃钢片表面 磁场强度幅值巧。^;
[0015] S13、利用相关计算模型确定乂
[0016] S14、计算磁场分布情况并确定Bs;
[0017] S15、判断是否满足|武"|-巧|《6(其中误差容限6取值为〇.〇11');若满足，则本次等 效平均磁导率尝试值即为有效的磁导率值，即乂，若不满足，改变Bs尝试值数值大小， 并重复W上过程，直至满足判别式条件为止。
[0018] 进一步地，步骤S13中，等效平均磁导率Weq通过求解W下非线性方程得到：
[0019]
[0020] 上式中，d为单片娃钢片厚度，δΑ和δ分别为基于极限非线性磁化模型和线性磁化 模型的集肤深度表述，且有
[0023] 其中σ为娃钢片材料电导率，ω为角频率。
[0024] 本发明的有益效果是：该饱和条件下无取向娃钢片叠片铁屯、磁导率张量表示方 法，充分考虑了叠片结构所引入的铁屯、磁属性各向异性特性，同时计入了叠片铁屯、中的磁 饱和效应，所建立的叠片铁屯、磁导率张量模型符合无取向娃钢片铁屯、实际磁化过程;基于 该铁屯、磁导率张量表示方法，能够实现对实际非线性电磁场问题的线性化简化处理，从而 一定程度上削减电磁场计算所需时间和人力成本。
【附图说明】
[0025] 图1是本发明的计算过程流程图；
[0026] 图2是本发明中等效平均磁导率Weq的计算流程图；
[0027] 图3是娃钢片叠片铁屯、及坐标系示意图；
[002引图4是算例中分析的叠片铁屯、模型实物图。
【具体实施方式】
[0029] 下面结合附图对本发明的实施例进行详细阐述，W使本专利的优点和特征能更易 于被本领域的技术人员理解，从而对本发明的保护范围做出更为清楚明确的界定。
[0030] 本发明一种饱和条件下无取向娃钢片叠片铁屯、磁导率张量表示方法，考虑了叠片 结构所引入的铁屯、磁属性各向异性特性，同时计入了叠片铁屯、中的磁饱和效应，所建立的 叠片铁屯、磁导率张量模型贴近叠片铁屯、实际磁化过程。
[0031] 该种饱和条件下无取向娃钢片叠片铁屯、磁导率张量表示方法，如图1，包括：
[0032] S1、借助相关计算模型、通过迭代计算确定等效平均磁导率μβ。具体数值。
[0033] 如图2，磁饱和情况下，求取等效平均磁导率Weq的具体步骤如下：
[0034] S11、设定饱和磁感应强度Bs的第k次尝试值公,W;
[0035] S12、根据娃钢片材料直流磁化曲线及关系式确定对应娃钢片表面 磁场强度幅值好
[0036] S13、利用相关计算模型确定乂
[0037] 等效平均磁导率Weq通过求解W下非线性方程得到：
[00；3 引
[0039]上式中，d为单片娃钢片厚度，δΑ和δ分别为基于极限非线性磁化模型和线性磁化 模型的集肤深度表述，且有
[0042] 其中σ为娃钢片材料电导率，ω为角频率；
[0043] S14、计算磁场分布情况并确定Bs;
[0044] S15、判断是否满足|啤'-5|^6(其中误差容限6取值为〇.〇11');若满足，则本次等 效平均磁导率尝试值即为有效的磁导率值，即;若不满足，改变Bs尝试值数值大小， 并重复W上过程，直至满足判别式条件为止。由于单值磁化曲线的连续性，总能找到合适的 化和BsW满足判别准则，因此整个计算过程不存在收敛性问题。
[0045] S2、计及叠片结构所引入的磁特性各向异性，若假设娃钢片厚度方向即叠片方向 沿直角坐标系的y轴方向，如图3，则所述磁导率张量基本形式为：
[0046]
[0047] 上式中，F为铁屯、叠片系数，μ〇为真空磁导率。
[004引有效性分析
[0049] 通过一个算例验证本发明提出的磁导率张量表示方法的有效性。
[0050] 分析对象为一个尺寸为200mm X 240mm X 72mm(外框宽度X高度X厚度）的框型铁 屯、模型，由140片厚度为0.5mm、宽度为40mm娃钢片（型号：50W600;叠片系数:0.97;电导率： 2.13MS/m)叠装而成。铁屯、励磁线圈和测量线圈应数均为220应，分别绕制于铁屯、上臂和下 臂。叠片铁屯、模型实物图如图4所示。测量过程使用可编程电源进行励磁，使用功率分析仪 进行电压、电流、功率等电参数测量。
[0051] 基于本发明提出的饱和条件下无取向娃钢片铁屯、磁导率张量表示方法，对图4所 示叠片铁屯、模型实施了线性时谐场Ξ维有限元分析。由于饱和情况下，励磁电流严重偏离 正弦波形，算例中Ξ维有限元分析过程选择W方均根(rms)电流测量值作为励磁载荷。感应 电动势和磁密幅值的测量值与计算结果比较见表1，表中选取了初始饱和及深度饱和两种 典型状态。如表1所示，基于本发明提出的张量磁导率所得到的有限元计算结果与测量数据 非常接近，验证了本发明的有效性和精确性。
[0052] 表1感应电动势和磁密幅值的测量值与计算值比较
[0化3]
[0054] W上对本发明的【具体实施方式】进行了描述，但本发明并不限于W上描述。对于本 领域的技术人员而言，任何对本技术方案的同等修改和替代都是在本发明的范围之中。因 此，在不脱离本发明的精神和范围下所作的均等变换和修改，都应涵盖在本发明的范围内。
【主权项】
1. 一种饱和条件下无取向硅钢片铁心磁导率张量表示方法，其特征在于，包括如下步 骤： 51、 借助相关计算模型、通过迭代计算确定等效平均磁导率具体数值； 52、 计及叠片结构所引入的磁特性各向异性，若假设硅钢片厚度方向即叠片方向沿直 角坐标系的y轴方向，则所述磁导率张量基本形式为：上式中，F为铁心叠片系数，μ〇为真空磁导率。2. 如权利要求1所述的饱和条件下无取向硅钢片铁心磁导率张量表示方法，其特征在 于，所述步骤S1中，通过相关计算模型及迭代计算求取等效平均磁导率μ%，具体步骤为： SI 1、设定饱和磁感应强度Bs的第k次尝试值尽1< >: 512、 根据硅钢片材料直流磁化曲线及关系式确定对应硅钢片表面磁场 强度幅值 513、 利用相关计算模型确定 514、 计算磁场分布情况并确定&尝试值； 515、 判断是否满足其中误差容限ε取值为0.01T);若满足，则本次等效平 均磁导率尝试值即为有效的磁导率值，即若不满足，改变队尝试值数值大小，并重 复以上过程，直至满足判别式条件为止。3. 如权利要求1所述的饱和条件下无取向硅钢片铁心磁导率张量表示方法，其特征在 于，所述步骤S13中，等效平均磁导率μ%通过求解以下非线性方程得到：〇上式中，d为单片硅钢 片厚度，SA和δ分别为基于极限非线性磁化模型和线性磁化模型的集肤深度表述，且有其中σ为硅钢片材料电导率，ω为角频率。 (3) (4)
【文档编号】G06F17/11GK106066846SQ201610400204
【公开日】2016年11月2日
【申请日】2016年6月8日 公开号201610400204.3, CN 106066846 A, CN 106066846A, CN 201610400204, CN-A-106066846, CN106066846 A, CN106066846A, CN201610400204, CN201610400204.3
【发明人】王坚, 李宏胜, 蒋春容
【申请人】南京工程学院