基于Comsol的单相变压器绕组与铁芯振动形变仿真全流程解析

写变压器振动仿真这个方向，我算是跟它磨合了挺长时间。玩Comsol的人很多，但大多数是在做电磁场或者散热，真正把“绕组+铁芯”的振动形变完整跑通的人不算多。今天就把我基于Comsol做的单相变压器绕组及铁芯振动形变仿真完整梳理一遍。围绕Comsol、单相变压器、绕组、铁芯、振动形变仿真这几个关键词，我会把从几何建模、物理场设置、网格剖分到结果解读的完整链路，包括那些在文档里根本查不到的坑，一块儿交代清楚。这篇东西适合正在用Comsol做多物理场耦合的工程师，也适合刚入门的硕士生照着落地一个能跑通的案例。

1. 从物理机理出发：绕组和铁芯为什么会振动

1.1 两个振动源，两条完全不同的物理路径

单相变压器在运行时的振动，主要来源并不是很多人以为的“电”本身，而是两个电磁-结构耦合效应：

第一是铁芯的磁致伸缩。硅钢片在外加交变磁场下，沿磁化方向会发生微小的尺寸变化，这个效应叫磁致伸缩。虽然相对形变通常只有几个ppm（百万分之一），但铁芯柱和铁轭的整体尺寸有几十厘米量级，叠加出来的微米级位移在100Hz的激励下持续作用，足以让变压器发出明显的“嗡嗡”声，并且长期下去会导致铁芯叠片松动、硅钢片间磨损。

第二是绕组的洛伦兹力振动。绕组导线处在漏磁场中，电流和磁场相互作用产生体积力。单相变压器负载运行时，低压侧大电流绕组在漏磁场下受到的电磁力尤其明显。特别是当发生外部短路时，绕组受到的电磁力可以比正常运行大几十倍到上百倍，这也是绕组变形、绝缘磨损甚至匝间短路的主要诱因。

这两个源的物理路径不一样：磁致伸缩力是铁芯内部的体积应变力，和磁场强度的平方近似成正比；绕组电磁力是外部磁场与载流导体的相互作用力，是典型的洛伦兹体积力。所以在Comsol里做仿真，必须把这两种力的加载方式分开处理，不能一股脑堆在一起，否则后处理时根本说不清哪个位移是哪个源贡献的。

1.2 为什么拿单相变压器做起点最合适

三相变压器当然更贴近工业现场，但做振动仿真探索时，我强烈建议先从单相变压器入手。原因有几点：

单相变压器的磁场分布和绕组排布在空间上更规整，激励源简单，方便你验证模型是否自洽。三相变压器还要额外处理三相磁路耦合、各相绕组之间的互感、磁路不平衡等问题，这些因素会直接干扰你对“振动形变仿真本身”的理解。单相变压器可以只用一相磁路+两个绕组（原边/副边）就能把核心问题跑通，计算规模也小很多。

另外，很多基础研究文献里的实验数据都是单相变压器上测的，便于你拿自己的仿真结果去对标。我的第一个完整模型就是一台1kVA、220V/110V的单相心式变压器，跑通了之后再迁移到三相模型就顺手多了。

1.3 多物理场耦合的整体思路

在Comsol里做这个仿真，本质上是电磁场 → 力 → 结构形变的单向或双向耦合。振动形变仿真最常用的做法是：

• 用**磁场（mf）**接口计算稳态或频域下的磁场分布、电流分布；
• 利用Comsol的多物理场耦合节点，把Maxwell应力张量和洛伦兹力计算出来；
• 把力作为载荷传递给**固体力学（solid）**接口，求解绕组和铁芯的形变与应力；
• 如果要看共振风险，再做预应力特征频率分析，把电磁力产生的预应力叠加到刚度矩阵上，计算结构在通电状态下的固有频率和振型。

多数情况下我先跑频域分析，因为工频50Hz下的电磁力是简谐的，频域直接算稳态更省时间。如果要做短路冲击或非线性饱和分析，再换成瞬态研究。

2. 几何建模与材料参数：这一步省了，后面全崩

2.1 几何细节做到什么程度合适

很多新手拿到变压器就恨不得把螺栓、垫片、绝缘纸全建出来，结果网格剖分失败或者计算量爆炸。我自己的经验是：第一版模型做“干净的简化几何”，只保留对电磁和结构响应影响最大的部件。

以我用的1kVA单相心式变压器为例，几何清单如下：

• 铁芯：心柱+上下铁轭，矩形截面，尺寸按实际硅钢片叠片后的外轮廓来，不要建每片硅钢片的间隙，那会让网格数量失控；
• 原边绕组：低压侧靠近铁芯柱，线圈区域做成长方体或圆柱体，材料设置为铜，等效电导率按实际导线填充率折算；
• 副边绕组：高压侧套在原边外面，同样做等效块体；
• 油道/气道或绝缘层：在绕组与铁芯之间留出空气域，这个空气域对漏磁场分布影响很大，不能省；
• 外壳和夹件：第一版完全可以不建，等你把核心模型跑通了再慢慢加。

提示：几何建模一定要开启“自动去除薄小区域”的思路。绕组绝缘纸、漆膜这类尺寸极小的薄层，在电磁分析中容易引起网格畸形，但在振动形变里贡献不大，第一版直接去掉。

2.2 材料参数：B-H曲线比你想的更重要

Comsol材料库里有内置铜和结构钢，但铁芯的硅钢片B-H曲线必须自己手动输入。默认的“Iron”材料B-H曲线往往和实际取向硅钢片差别很大，特别是饱和区，直接决定铁芯磁阻和绕组电感，进而影响电磁力的计算精度。

我用的硅钢片参数大概是这样：

B-H曲线建议至少给10~15个点，从B=0一直到B=1.9T以上。数据可以从硅钢片手册里摘，也可以直接引用同类文献。Comsol里用的是插值函数，注意横纵坐标别填反了。

绕组铜的区域我通常按“等效铜块”处理：把铜导线的截面积、匝数和总几何体积折算成等效电导率。公式很简单：

[
\sigma_{eq} = \sigma_{cu} \times \frac{A_{cu}}{A_{window}}
]

其中 (A_{cu}) 是导线实际铜截面积总和，(A_{window}) 是绕组窗口面积。这样做的原因是让你不需要逐匝建模，大幅降低网格量，但电磁力积分仍然比较准。

2.3 空气域和边界条件的设置细节

单相变压器仿真的空气域不需要太大，我通常取变压器外轮廓向外扩展30%~50%即可。磁场接口的默认边界条件是磁绝缘，这已经近似认为空气域外边界磁通不穿出，作为变压器这类闭合磁路设备是合理的。

关键是绕组区域的激励设置。单相变压器原边接220V交流电压源，副边接负载。在Comsol的磁场接口里，可以用“线圈”特征给绕组加上电压激励，副边可以用“线圈”设置成电流源或接外部电路。我的做法是：

• 原边线圈（220V侧）：使用线圈类型“电压”，设定匝数、导线截面和电压220V；
• 副边线圈（110V侧）：使用线圈类型“电流”，设定电流大小按负载折算；
• 两个线圈要勾选“线圈方向”由几何体现，Comsol会自动根据几何方向计算电流方向。

如果不接入外部电路，直接在频域激励下跑，也能得到一个近似的额定工作点。但如果要做负载特性分析，建议用“电路”接口把两个线圈和负载电阻连起来，这样电流才能随负载自适应。

2.4 绕组和铁芯之间的机械连接关系

振动形变仿真的结构边界条件要和真实变压器尽量贴近。我第一版仿真时偷懒，把绕组和铁芯接触面全设成“固定约束”，算出来的位移几乎全被铁芯憋死了，和实测数据差距非常大。

实际上，绕组是通过绝缘筒、撑条和压板固定在铁芯上的，并不是刚性地“焊死”在铁芯表面。绕组和铁芯之间存在接触和预紧力。简化的做法是：铁芯底部和夹件接触的位置设置固定约束，铁芯与绕组之间设置成“接触”或加一个虚拟的弹性边界条件，用弹簧基础模拟绝缘支撑的等效刚度。

如果只做定性分析，用“弹簧基础”是最省事的，弹簧刚度按绝缘撑条的压缩刚度估算。如果做定量研究，就得建接触对，缺点是非线性求解时间明显增加。我在第一个探索模型里用的就是弹簧基础，结果和实验位移幅值能对到同一个数量级，已经足够指导后续设计。

3. 物理场与耦合设置：把电磁力正确地嫁接到结构上

3.1 磁场接口的核心设置

磁场接口我一般选择“稳态”或“频域”。“频域”可以直接研究50Hz激励下的正弦稳态响应，默认求解复数形式，后处理时显示的位移幅值就是绕组铁芯的振动幅值。

绕组区域要设置安培定律，铜的电导率按上文说的方法折算。铁芯区域如果是叠片结构，在磁场里可以勾选“等效叠片”，让Comsol自动处理叠片方向的涡流抑制效应，避免三维模型计算出虚假的大涡流。这个小细节影响非常大，我一开始没勾选，铁芯涡流损耗计算值比实际大了好几倍，电磁力也跟着失真。

求解时，非线性B-H曲线要求使用“非线性”求解器。Comsol默认会自动选择牛顿法，但对B-H曲线这种强非线性材料，最好在求解器设置里打开“辅助扫描”，把电流从0逐步加载到额定值，这样收敛要稳得多。

3.2 电磁力的三种加载方式，我用下来最顺手的方式

Comsol里把电磁场计算结果转化为结构力，有三种常见做法：

1. 多物理场耦合节点：直接在“多物理场”节点添加“电磁力”耦合，选择磁场和固体力学接口，Comsol会自动计算边界上的Maxwell应力张量和导体区域内的洛伦兹力。这是我最推荐的方式，因为作用域会正确耦合，不需要手动干预。

2. 手动添加边界载荷：在固体力学接口里，给铁芯表面添加边界载荷，载荷表达式用磁场模块计算出的Maxwell表面应力。好处是你能清楚看到力的大小和方向，坏处是容易漏掉体积力部分，绕组区域的洛伦兹力不能靠这种方式加。

3. 体积力手动计算：在绕组区域添加体积力，表达式用 ( \mathbf{J} \times \mathbf{B} )。这个方法最直观，也很适合教学演示，但需要自己对坐标分量写表达式，方向很容易搞错。

我个人的建议是：体积力用多物理场耦合节点，如果你要研究磁致伸缩，再单独叠加一个各向异性应变源。

3.3 磁致伸缩的等效处理：从没有数据到跑起来

磁致伸缩数据很难找，因为硅钢片的磁致伸缩系数受材料牌号、热处理、轧制方向影响很大，而且是非线性的。文献里常见值在 (5\times 10^{-6} \sim 30\times 10^{-6}) 左右，也就是微应变级别。

Comsol其实没有内置磁致伸缩材料模型，但可以用热膨胀类比来做：磁致伸缩本质上是磁场诱导的应变，而热膨胀是温度诱导的应变，数学形式完全一样。具体做法是：

• 添加一个“热膨胀”节点，但温度场不实际求解；
• 把“温度变化”替换成磁场强度的函数，比如用 ( \Delta T = \lambda_s \cdot (B/B_{sat})^2 )；
• 各向异性的体现：在铁芯柱方向（轧制方向）设置膨胀系数为磁致伸缩系数，其他方向设置的小一些或接近0。

虽然这只是个工程近似，但对振动形变仿真来说完全够用。跑出来的铁芯振动位移量级大概是几微米，频率是100Hz，这个和实测的变压器铁芯振动特征非常吻合。

3.4 频域计算后再做预应力模态分析

要判断变压器会不会发生共振，必须在通电状态下做模态分析。不能直接用无应力状态下的固有频率，因为电磁力会在结构里产生预应力，而预应力会改变结构刚度，进而改变固有频率。

操作流程是：先做一次频域稳态计算，得到电磁力分布；然后在“研究”里添加一个“特征频率”研究，在“固体力学”接口的“特征频率”设置里勾选“包含预应力”，并选择刚才频域解的载荷作为预应力来源。这样得到的固有频率就是工频工作点附近的真实固有频率。

重点检查100Hz及其倍频（200Hz、300Hz、400Hz）附近有没有固有频率。如果某一阶固有频率正好落在这些激励频率附近，意味着共振风险很大，后续就要考虑加筋、改变支撑刚度或者调整绕组压紧力来错开频率。

4. 网格剖分与求解器调优：耐心是唯一的捷径

4.1 网格剖分的三个“必须加密”区域

网格剖分直接影响电磁力和结构位移的求解精度。我在调试过程中总结出三个必须重点加密的区域：

第一是气隙和绕组间绝缘区域。漏磁场主要集中在这里，绕组电磁力的源头就是漏磁场，如果网格太粗，电磁力会在空间上被平滑掉，绕组受力算小甚至算成负的。这个区域我通常用边界层网格，至少要有三层以上单元落在气隙里。

第二是铁芯和绕组表面。Maxwell应力张量是通过表面场值积分得到的，表面网格越细，应力积分越准。我习惯在铁芯表面加一层扫掠网格，法向尺寸控制在0.5mm以内。

第三是绕组厚度方向。洛伦兹体积力随径向位置变化很大，如果绕组区域只用一层单元，算出来的体积力空间分布是“平均场”，形变结果会偏刚性。给绕组厚度方向上切两到三层网格，位移曲线就平滑多了。

4.2 网格质量检查与单元类型选择

电磁-结构耦合模型，我统一用二阶单元。一阶单元在弯曲变形问题上太硬，算出来的位移偏小；三阶单元成本翻好几倍，收益有限。二阶四面体就能满足大部分需求。

网格做好后一定跑一遍“统计”，重点看最差单元质量，不要低于0.1，低于这个值基本就是几何里有小尖角、小薄片导致的。出现这种问题优先修几何，不要硬画网格。

4.3 求解器设置的实战经验

非线性B-H曲线带来的最大问题是迭代不收敛。我的经验是打开“辅助扫描”，让电流从20%、40%、60%、80%、100%分步加载。每步都用上一步的解做初值，收敛概率会大幅提升。

频域求解时，如果网格量在20万自由度以上，默认的直接求解器MUMPS可能内存吃紧。可以换PARDISO，它对内存管理更好，速度也更快。如果内存还是不够，再考虑迭代求解器，但需要仔细调预条件子，否则容易发散。

结构部分的形变计算，由于位移量级是微米级，结构刚度矩阵数值相对电磁场来说大得离谱，我建议把两个物理场分开求解（单向耦合），先算磁场，再算结构。用“分离式”求解器，两个物理场交替迭代，比全耦合稳定得多，速度快好几倍。

4.4 常见问题排查速查表

5. 从仿真结果到设计决策：振动形变不能只停在云图

5.1 绕组电磁力分布的三个典型特征

跑通模型之后，你会看到绕组电磁力分布非常不均匀。以我那个单相心式变压器模型为例，典型的特征有三个：

绕组端部受力明显大于中部。这是端部漏磁场弯曲导致的，端部磁场径向分量大，和电流作用产生的轴向力就更强。这就是为什么实际变压器绕组端部要额外加绝缘支撑和端圈压紧，仿真结果可以清楚告诉你多大的力应该由端部绝缘承担。

低压绕组和高压绕组受力方向相反。因为两个绕组电流方向相反，洛伦兹力方向也相反。这意味着低压绕组被向内压，高压绕组被向外推。在短路工况下这对力可能达到数千牛顿级别，直接导致绕组径向变形。

铁芯柱高度方向的轴向力累积明显。从铁芯底部到顶部，轴向力是逐匝累积的，越靠近端部累计力越大。这个轴向力是导致铁芯上下振动和绕组轴向位移的主要来源，在结构边界条件里要重点考虑压板的压紧力是否能平衡这部分力。

5.2 铁芯形变特征：100Hz是主频，但别忽略倍频

在频域结果中，铁芯的振动形变主频率是100Hz，这是磁致伸缩的典型特征：无论磁场方向怎么变化，磁致伸缩的长度变化都是两个周期内完成一次伸缩，所以激励频率翻倍。

但铁芯位移的频谱里也会有50Hz成分和150Hz、200Hz成分。50Hz成分主要来自铁芯磁路的非线性——当B-H曲线进入饱和区，磁通波形发生畸变，磁致伸缩就不再是标准的“平方关系”，而是会产生基频分量。这在设计铁芯时很有参考价值：如果铁芯设计得过于饱和，除了损耗上升，铁芯振动的50Hz成分也会变大，人耳听起来会有一种额外的“低沉振动”感。

5.3 结构形变的工程解读：位移几微米也要管吗

很多同学第一次跑出位移结果后，看到铁芯形变量才几个微米，就觉得“这么小，应该没事吧”。这个判断需要修正。微米级的位移本身不会直接损坏结构，但以下几点值得关注：

绝缘磨损是长期疲劳问题。绕组在100Hz振动下，即使振幅只有几十微米，经过几万小时运行，也会和绝缘撑条、垫块产生微动磨损。磨损产生的粉末会降低绝缘性能，还可能堵塞油道。所以变压器的振动标准关注的是振动速度和加速度，而不是绝对位移。

固有频率错开比振幅大小更重要。只要绕组的固有频率和100Hz及其倍频错开一定距离（一般要求错开10%以上），即使电磁力很大，结构也不会持续积累能量。反过来，如果恰好共振，几十微米的幅值会放大到数百微米，那就真的会出问题。

正因如此，预应力模态分析的结论往往比振动位移云图更能指导设计。我做完仿真后，给变压器厂家提的建议通常是调整绝缘垫块的刚度和位置，把绕组系统的固有频率从95Hz附近挪到110Hz以上，效果非常明显。

5.4 从单相到三相：这套方法可以直接迁移

单相变压器模型跑通了，往三相扩展的思路是清晰的：

• 几何上增加三柱铁芯（或壳式铁芯）和三个绕组；
• 电激励上用三相电压源，相位分别差120度；
• 磁场和结构耦合方式不变；
• 分析重点要增加各相之间磁路耦合带来的振动相位差。

三相变压器的一个有趣现象是：如果三相磁路完全对称，铁芯的某些振动分量会相互抵消；但实际制造中叠片工艺必然存在差异，对称性并不完美，所以实测铁芯振动会比理想模型更复杂。要研究这种非对称效应，三维模型就绕不过去了。

另外还可以继续扩展的方向包括：把结构响应进一步和声学接口耦合，用频域声压计算变压器远场噪声；或者做短路瞬态分析，把瞬态电磁力导入结构瞬态求解器，观察绕组在短路冲击下的动态变形过程。这些方向一旦跑通，基本就是一篇学位论文的核心内容了。

6. 实操心得：那些文档里不写的经验

最后再说几个我实际踩过、后来彻底改过来的习惯。

先跑二维轴对称，再上三维。 单相变压器完全可以用二维轴对称简化分析，计算时间从几十分钟降到几十秒。先用二维摸清物理规律，确定力的大小和方向都合理了再上三维做精细结构响应。我见过太多人一上来就建全三维模型，结果边界条件设错，白跑了好几天。

材料参数一定要写注释和来源。 Comsol模型文件里可以给材料节点写描述，我强烈建议把B-H曲线的数据来源、硅钢片牌号和折算电导率的计算过程都记进去。否则三个月后回来看模型，你根本记不清那组数据是怎么来的，更不要说当审稿人或同事质问你数据来源时你会非常被动。

用实验数据“反向标定”模型。 如果你手头有加速度传感器或者激光测振仪，哪怕是简化的实验结果，也要拿来做对标。我通常把仿真位移幅值和实验振动速度做换算对比。对不上不要急着改网格，先检查边界条件和激励幅值。很多时候误差来自材料参数和实际设备有出入，这时候要把“材料等效”这一步做仔细，而不是盲目加密网格。

这套模型跑通之后，你会明显感觉到自己对变压器运行时“看不见的力”有了更直观的认识。绕组和铁芯的振动不是玄学，而是可以从电磁场强度、材料磁致伸缩系数和结构刚度一步步算出来的工程问题。把这套链路固化下来，后面无论做电机、电抗器还是变压器，都只是换一个几何和激励源的事。