一种基于COMSOL的铁心振动特性计算方法技术

本发明专利技术提供了一种基于COMSOL的铁心振动特性计算方法，该方法利用COMSOL Multiphysics仿真软件，在间谐波激励下，仿真计算含铁心绕组类设备铁心的磁通密度以及振动位移，包括S1：获取铁心模型的结构参数、电气参数；S2：几何建模与材料参数设置；S3：物理场设置、网格剖分及求解计算共三个步骤。本发明专利技术可以更加省时、高效对比多种不同激励下振动情况，并且在结构力场分析中考虑了各向异性，计算磁致伸缩时更加准确。

A calculation method of core vibration characteristics based on COMSOL

【技术实现步骤摘要】
一种基于COMSOL的铁心振动特性计算方法
本专利技术属于含铁心绕组类设备铁心振动特性计算领域，尤其涉及间谐波激励下，铁心的磁通密度以及振动位移的计算方法。
技术介绍
随着我国电力系统的不断发展和成长，电力电子元件、非线性元件及负荷接入电网，带来严重的间谐波污染。作为电力系统中的重要电气设备，含铁心绕组类设备在间谐波的影响下发生电流畸变，并在很大程度上产生更严重的噪声问题，影响附近居民生产生活。利用电路-磁场-固体力学耦合仿真计算方式能够较准确反映铁心的振动位移。因此，可以通过使用计算机软件对铁心进行模型的建立，并计算间谐波激励下铁心振动特性。本专利技术的目的在于填补现有对于间谐波激励下变压器铁心振动仿真的空白，完成COMSOLMultiphysics仿真软件下铁心在磁场与物理场下的模型建立，并保证其耦合的准确性，并对间谐波激励下，铁心磁通密度分布及振动位移进行仿真计算。
技术实现思路
本专利技术提供了一种基于COMSOL的铁心振动特性计算方法，利用COMSOLMultiphysics仿真软件，仿真计算铁心在间谐波激励下的磁通密度分布以及振动位移，包括以下步骤：步骤1：获取铁心模型的结构参数、电气参数；根据实际应用场景，通过现场测量获取含铁心绕组类设备铁心的窗口高度及宽度，铁心模型厚度、宽度和高度，绕组匝数、高度，利用单片测量仪器获取硅钢片H-B和λ-B曲线，根据生产商提供参数获取硅钢片的杨氏模量、泊松比、密度、相对介电常数、电导率以及绕组的相对介电常数、电导率、相对磁导率，根据资料获取空气相对介电常数、电导率、相对磁导率；步骤2：基于COMSOL进行铁心模型几何建模与材料参数设置；根据步骤1中根据获取铁心的窗口高度及宽度，设置COMSOLMultiphysics软件中的窗口参数；根据所获取的铁心模型叠厚、宽度、高度，设置COMSOLMultiphysics软件中的铁心模型参数；根据获取的绕组匝数、尺寸，设置COMSOLMultiphysics软件中的绕组参数；根据获取的硅钢片型号设置COMSOLMultiphysics软件中的电场、磁场、固体力学场参数；步骤3：在COMSOL几何模型中进行物理场设置、网格剖分及求解计算；根据步骤2所建立的铁心的几何模型，在COMSOLMultiphysics软件中通过电路-磁场-固体力学-声学耦合实现计算；将绕组域、空气域、铁心部分进行网格剖分，设置求解步长0.5ms和求解时间3s；仿真计算不同间谐波激励条件下，铁心的磁通密度分布和振动位移。其中，步骤3中所述的在COMSOLMultiphysics软件中通过电路-磁场-固体力学-声学耦合实现计算的过程包括：在计算电场和磁场耦合过程时，引入矢量磁位A，使A满足：另外，电流密度还满足方程：式中，Je是源电流密度，N是绕组线圈匝数，Vcoil是在绕组上施加的电压激励，Rcoil是绕组电阻，ecoil是载流线圈各点的电流密度的方向向量，S是绕组线圈截面积；将变压器铁心看作稳态场，忽略位移电流，得到计算方程：式中，μ0、μr分别为真空磁导率和相对磁导率，t为时间(s)，σ为介质电导率；联立方程式(1)、(2)、(3)，完成电路和磁场的耦合；根据步骤1中对硅钢片单片测量得到的磁致伸缩特性，将轧制方向和垂直于轧制方向的磁通密度B和磁致伸缩峰峰值λ的函数作为初始预应变ε0加入模型；再通过设置铁心材料杨氏模量和泊松比，利用胡克定律，得到线弹性材料应力与应变关系，在坐标轴x轴方向上满足方程(4)：式中，εx表示x轴方向应变，E、ν分别表示变压器铁心杨氏模量和泊松比，σx，σy，σz分别表示x轴，y轴，z轴所受应力；物理场中铁心任一点满足的微分方程(5)：式中，u表示物体内任一点位移矢量，ρ表示物体密度，表示物体任一点各方向应力分量，Fv表示磁场物理场计算得到的洛仑兹力；在物理场求解时，将铁心设为自由振动，同时为模拟重力影响，需要将铁心底部设置为固定约束，即满足方程(6)：u＝0(6)从而计算铁心振动完成磁场与固体力学的耦合；假设空气为理想气体，即空气域为静止、均匀、连续、绝热的理想流体介质，且假设声波是小振幅波，即振幅相对于波长为无穷小，利用波动方程(7)表示声音在空气中传播状态：式中，ρc表示流体密度，pt表示声压，qd表示偶极源，Qm表示单极源；将声源简化为单极源和偶极源两种声源，单极源是进行均匀涨缩振动的球面声源，球源表面上各点沿径向做同振幅，同相位振动，且半径远小于声波波长，偶极源是两个振幅相位相反，而大小相等的相互靠近的简单源，令k＝ω/c，k表示波数，ω表示角频率，c表示声速，(7)可以简化为亥姆霍兹方程(8)：为求解铁心在周围空间中产生的声压，将固体力学中求解得到的铁心法向振动加速度作为初始值进行计算，如方程(9)所示：假设模型中所有边界均为完美刚体，相当于法向速率为0，添加硬声场边界作为边界条件，满足方程(10)：通过振动加速度和声压级完成了固体力学和声场的耦合。优选地，步骤1中的铁心的窗口尺寸为宽70mm、高130mm。优选地，步骤1中的铁心模型为高度宽200mm、长260mm、厚66mm的单相两柱机构；两绕组为内部直径106mm、外部直径120mm、高100mm的管体结构。更优选地，步骤2中所设置的电场、磁场、固体力学场参数为：铁心、绕组和空气域相对介电常数都为1，电导率为10S/m，绕组和空气域相对磁导率为1，铁心材料杨氏模量为2×1011Pa，泊松比为0.28，密度为7650kg/m3。进一步优选的，步骤3中将激励条件设置为45V/50Hz，分别在附加1％含量的0.45V/5Hz的间谐波、3％含量的1.35V/5Hz的间谐波、3％含量的1.35/20Hz的间谐波三种情形下，仿真计算铁心的磁通密度分布和振动位移。COMSOLMultiphysics软件作为电磁仿真软件，相比其他有限元仿真优于多物理场耦合计算，可以更加省时、高效对比多种不同激励下振动情况，并且在结构力场分析中考虑了各向异性，计算磁致伸缩时更加准确。附图说明图1为本专利技术一种实施例的仿真计算流程图。图2为基于COMSOLMultiphysics建立的铁心模型。图3为铁心设置的B-H磁化曲线。图4为铁心设置的λ-B曲线。图5a为45V/50Hz激励下铁心模型磁通密度分布(单位/T)。图5b为45V/50Hz激励下铁心振动位移云图(单位/微米)。图6a为45V/50Hz附加1.35V/5Hz间谐波激励下铁心模型磁通密度分布(单位/T)。图6b为45V/50Hz附加1.35V/5Hz间谐波激励下铁心振动位移云图(单位/微米)。具体实施方式下面结合附图和实施例本文档来自技高网...

【技术保护点】
1.一种基于COMSOL的铁心振动特性计算方法，利用COMSOL Multiphysics仿真软件，仿真计算含铁心绕组类设备铁心在间谐波激励下的磁通密度分布以及振动位移，包括以下步骤：/n步骤1：获取铁心模型的结构参数、电气参数；/n根据实际应用场景，通过现场测量获取含铁心绕组类设备铁心的窗口高度及宽度，铁心模型厚度、宽度和高度，绕组匝数、高度,利用单片测量仪器获取硅钢片H-B和λ-B曲线，根据生产商提供参数获取硅钢片的杨氏模量、泊松比、密度、相对介电常数、电导率以及绕组的相对介电常数、电导率、相对磁导率，根据资料获取空气相对介电常数、电导率、相对磁导率；/n步骤2：基于COMSOL进行铁心模型几何建模与材料参数设置；/n根据步骤1中根据获取铁心的窗口高度及宽度，设置COMSOL Multiphysics软件中的窗口参数；根据所获取的铁心模型叠厚、宽度、高度，设置COMSOL Multiphysics软件中的铁心模型参数；根据获取的绕组匝数、尺寸，设置COMSOL Multiphysics软件中的绕组参数；根据获取的硅钢片型号设置COMSOL Multiphysics软件中的电场、磁场、固体力学场参数；/n步骤3：在COMSOL几何模型中进行物理场设置、网格剖分及求解计算；/n根据步骤2所建立的铁心的几何模型，在COMSOL Multiphysics软件中通过电路-磁场-固体力学-声学耦合实现计算；/n将绕组域、空气域、铁心部分进行网格剖分，设置求解步长0.5ms和求解时间3s；仿真计算不同间谐波激励条件下，铁心的磁通密度分布和振动位移。/n...

【技术特征摘要】
1.一种基于COMSOL的铁心振动特性计算方法，利用COMSOLMultiphysics仿真软件，仿真计算含铁心绕组类设备铁心在间谐波激励下的磁通密度分布以及振动位移，包括以下步骤：
步骤1：获取铁心模型的结构参数、电气参数；
根据实际应用场景，通过现场测量获取含铁心绕组类设备铁心的窗口高度及宽度，铁心模型厚度、宽度和高度，绕组匝数、高度,利用单片测量仪器获取硅钢片H-B和λ-B曲线，根据生产商提供参数获取硅钢片的杨氏模量、泊松比、密度、相对介电常数、电导率以及绕组的相对介电常数、电导率、相对磁导率，根据资料获取空气相对介电常数、电导率、相对磁导率；
步骤2：基于COMSOL进行铁心模型几何建模与材料参数设置；
根据步骤1中根据获取铁心的窗口高度及宽度，设置COMSOLMultiphysics软件中的窗口参数；根据所获取的铁心模型叠厚、宽度、高度，设置COMSOLMultiphysics软件中的铁心模型参数；根据获取的绕组匝数、尺寸，设置COMSOLMultiphysics软件中的绕组参数；根据获取的硅钢片型号设置COMSOLMultiphysics软件中的电场、磁场、固体力学场参数；
步骤3：在COMSOL几何模型中进行物理场设置、网格剖分及求解计算；
根据步骤2所建立的铁心的几何模型，在COMSOLMultiphysics软件中通过电路-磁场-固体力学-声学耦合实现计算；
将绕组域、空气域、铁心部分进行网格剖分，设置求解步长0.5ms和求解时间3s；仿真计算不同间谐波激励条件下，铁心的磁通密度分布和振动位移。


2.如权利要求1所述的一种基于COMSOL的铁心振动特性计算方法，其中，步骤3中所述的在COMSOLMultiphysics软件中通过电路-磁场-固体力学-声学耦合实现计算的过程包括：
在计算电场和磁场耦合过程时，引入矢量磁位A，使A满足：



另外，电流密度还满足方程：



式中，Je是源电流密度，N是绕组线圈匝数，Vcoil是在绕组上施加的电压激励，Rcoil是绕组电阻，ecoil是载流线圈各点的电流密度的方向向量，S是绕组线圈截面积；
将变压器铁心看作稳态场，忽略位移电流，得到计算方程：



式中，μ0、μr分别为真空磁导率和相对磁导率，t为时间(s)，σ为介质电导率；
联立方程式(1)、(2)、(3)，完成电路和磁场的耦合；
根据步骤1中对硅钢片单片测量得到的磁致伸缩特性，将轧制方向和垂直于轧制方向的磁通密度B和磁致伸缩峰峰值λ的函数作为初始预应变ε0加入模型；
再通过设置铁心材料杨氏模量和泊松比，利用胡克定律，得到线弹性材料应力与应变关系，在坐标轴x轴方向上满足方程(4)：


<...

【专利技术属性】
技术研发人员：韩学，宋文乐，李慧奇，杨光，蒋曦，刘翔宇，刘东，代会荣，王思涵，王笑凯，
申请(专利权)人：国网河北省电力有限公司沧州供电分公司，华北电力大学保定，
类型：发明
国别省市：河北;13