基于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单相心式变压器为例，几何清单如下：

• 铁芯：心柱+上下铁轭，矩形截面，尺寸按实际硅钢片叠片后的外轮廓来，不要建每