磁気浮上
理论与物理
磁悬浮原理
老师，磁悬浮列车是怎么浮起来的？
主要有3种方式。
1. EDS（电磁感应方式） — 超导磁铁通过地面线圈时，涡流产生排斥力。JR式磁悬浮
2. EMS（电磁吸引方式） — 电磁铁的吸引力将列车拉向轨道。Transrapid
3. 涡流排斥方式 — 永磁体/电磁铁在导体板上移动产生悬浮
涡流排斥力：
速度$v$越大，悬浮力越大。但同时也会产生阻力（制动力）。
所以需要达到一定速度以上才能悬浮起来啊。
对于EDS方式，悬浮起始速度约为150 km/h。低速时依靠车轮行驶。
总结
- EDS — 涡流排斥。超导磁铁+地面线圈
- EMS — 电磁铁吸引。必须进行反馈控制
- 悬浮力$\propto v B^2$ — 速度依赖
磁悬浮力学——恩绍定理与“静磁场无法悬浮”的原因
斯科特·恩绍于1842年证明的“恩绍定理”表明，仅靠静态电磁力无法使带电体/磁性体稳定悬浮。这意味着磁悬浮本质上是一项动态的（依赖于控制的）技术，而超导体实现的完全抗磁性（迈斯纳效应）或抗磁性体的悬浮是作为该定理的例外而成立的。磁悬浮列车的超导磁悬浮也可以从这个角度理解，这是“无需控制即可实现稳定悬浮”这一超导优势的理论依据。
各项的物理意义
- 电场项 $\nabla \times \mathbf{E} = -\partial \mathbf{B}/\partial t$：法拉第电磁感应定律。随时间变化的磁通密度产生电动势。【日常例子】自行车发电机（发电机）通过旋转磁铁使附近的线圈产生电压——这是磁场随时间变化会感应出电场这一定律的直接应用。IH电磁炉也基于相同原理，高频磁场的变化在锅底感应出涡流，通过焦耳热加热。
- 磁场项 $\nabla \times \mathbf{H} = \mathbf{J} + \partial \mathbf{D}/\partial t$：安培-麦克斯韦定律。电流和位移电流产生磁场。【日常例子】电线通电时周围会产生磁场——这就是安培定律。电磁铁基于此原理工作，通过给线圈通电产生强磁场。智能手机的扬声器也应用了此定律：电流→磁场→振膜的力。在高频（GHz频段天线等）情况下，位移电流 $\partial D/\partial t$ 不可忽略，用于描述电磁波辐射。
- 高斯定律 $\nabla \cdot \mathbf{D} = \rho_v$：表明电荷是电通量的发散源。【日常例子】用垫板摩擦头发会产生静电使头发竖起——带电的垫板（电荷）向外辐射电力线，对轻质的头发施加力。电容器设计中，电极间的电场分布用此定律计算。ESD（静电放电）防护也基于高斯定律的电场分析。
- 磁通守恒 $\nabla \cdot \mathbf{B} = 0$：表示不存在磁单极子。【日常例子】将条形磁铁切成两半也无法得到只有N极或只有S极的磁铁——N极和S极总是成对出现。这意味着磁力线描绘的是“无始无终的闭合回路”。在数值分析中，为了满足此条件，采用矢量势 $\mathbf{B} = \nabla \times \mathbf{A}$ 的公式化方法，自动保证磁通守恒。
假设条件与适用范围
- 线性材料假设：磁导率·介电常数不依赖于磁场·电场强度（饱和区域需要非线性B-H曲线）
- 准静态近似（低频）：位移电流项可忽略（$\omega \varepsilon \ll \sigma$）。涡流分析中常用
- 2D假设（截面分析）：电流方向均匀，可忽略端部效应时有效
- 各向同性假设：各向异性材料（如硅钢板的轧制方向等）需要定义方向特性
- 不适用的案例：等离子体（电离气体）、超导体、非线性光学材料需要额外的本构关系
数值解法与实现
FEM中的悬浮力分析
磁悬浮的FEM分析怎么做？
求解运动磁铁/线圈与导体的相互作用。方法：
- 滑移网格 — 分离运动体与固定体的网格，在界面处耦合
- 对流项 — 固定导体，用$\sigma(\mathbf{v}\times\mathbf{B})$表示运动
- 时域瞬态分析 — 包含速度变化的动态分析
悬浮力通过麦克斯韦应力张量或虚功法计算。
EMS的控制仿真也用FEM吗？
EMS中电磁铁的电流通过PID控制反馈。将FEM与控制电路（MATLAB/Simulink等）耦合。JMAG具备Simulink联动功能。
总结
磁悬浮的FEM——非线性磁力的查找表生成
在磁悬浮系统的FEM分析中，通常会预先计算电磁铁电流值与悬浮间隙组合对应的电磁力，生成查找表（LUT）。LUT的分辨率（电流·间隙的网格数）和插值精度决定了控制仿真的精度。为了准确纳入非线性磁饱和，需要在每个点收敛牛顿法，这需要进行数百至数千个FEM计算并行执行的大规模计算。
边单元（Nedelec单元）
专用于电磁场分析的单元。自动保证切向分量的连续性，排除伪模式。是3D高频分析的标准。
节点单元
用于标量势公式化。在静磁场的标量势法或静电场分析中有效。
FEM vs BEM（边界元法）
FEM: 对应非线性材料·非均匀介质。BEM: 自然处理无限域（开域问题）。混合FEM-BEM也有效。
非线性收敛（磁饱和）
用牛顿-拉夫森法处理B-H曲线的非线性。残差标准: $||R||/||R_0|| < 10^{-4}$ 为通用标准。
频域分析
通过时间谐波假设归结为稳态问题。需要进行复数运算，但宽带特性需通过时域分析获取。
时域的时间步长
需要最高频率成分的1/20以下的时间步长。隐式时间积分中也可使用更大的步长，但需注意精度。
频域与时域的使用区分
频域分析类似于“将收音机调到特定频率”——可以高效计算单一频率下的响应。时域分析类似于“同时录制所有频道”——可以再现包含所有频率成分的瞬态现象，但计算成本高。
实践指南
实际业务中的设计
磁悬浮列车、磁轴承、输送装置、磁悬浮熔炼是典型应用。
实际业务检查清单
- [ ] 悬浮力是否充分超过自重（安全系数1.5以上）
- [ ] 横向恢复力（导向力）是否得到保证
- [ ] 推进系统能否补偿制动力导致的速度损失
- [ ] EMS的控制响应速度能否跟上间隙变化
- [ ] 超导磁铁是否采取了失超防护措施
中央新干线的超导悬浮——将磁悬浮“实用化”的真实数据
中央新干线超导磁悬浮系统以约10cm的间隙悬浮车体，以时速500km以上行驶。通过在超导线圈中通入永久电流在车上产生强磁场，利用与地面推进·悬浮·导向线圈的电磁力一体化控制悬浮·推进·导向。要使用FEM分析这种复杂的电磁力相互作用，需要数十万至数百万单元的模型，JRTT（原交通·观光·气象）的分析成果是日本电磁CAE技术最大规模的应用案例。
分析流程的比喻
电机的电磁场分析感觉上类似于“给吉他调音”。调整琴弦粗细（线圈匝数）和琴桥位置（磁铁配置），以引出最美妙的音色（高效的扭矩特性）。改变一个参数，整体平衡就会改变——所以参数化研究很重要。
初学者容易陷入的误区
“空气区域？为什么要用网格划分空气？”——这是几乎所有初次接触电磁场分析的人都会产生的疑问。答案是“因为磁力线也会延伸到铁芯之外”。如果将分析区域刚好设在铁芯边缘，无处可去的磁通会“撞上”边界并反射，产生实际中不可能出现的磁通集中。想象一下房间太小，球不断弹到墙壁上的状态。
边界条件的思考方式
远场边界条件看似不起眼但极其重要。需要在数值上表达“从这里开始是无限延伸的空间”。如果设置错误，磁通就会像撞上“看不见的墙”一样被反射。
软件比较
工具
| 工具 | 特点 |
|---|---|
| JMAG | Motion耦合。悬浮力的速度依赖计算 |
| Ansys Maxwell | 3D涡流+Motion。瞬态分析 |
| COMSOL AC/DC | 移动网格。多物理场耦合 |
| Opera (Dassault) | 在超导磁铁+悬浮分析方面有实际业绩 |
磁悬浮设计工具——JMAG-RT与MATLAB/Simulink的联动
在磁悬浮系统设计中，使用FEM（JMAG）预先计算电磁力的非线性映射，并将其结果嵌入MATLAB/Simulink的控制系统仿真中，这种工作流程已成为标准。使用JMAG-RT（实时模型），可以在保持FEM计算精度的同时进行实时控制仿真。磁悬浮式洗衣机电机（夏普等）的设计中也使用了此方法，是CAE在民用产品设计中应用扩展的案例。
选型时最重要的三个问题
- “要解决什么问题”：所需的物理模型·单元类型是否支持磁悬浮。例如，流体方面是否有LES支持，结构方面接触·大变形的处理能力会成为差异点。
- “谁来使用”：新手团队适合GUI完善的工具，有经验者适合脚本驱动的灵活工具。类似于汽车的AT车（GUI）和MT车（脚本）的区别。
- “未来扩展到什么程度”：着眼于未来的分析规模扩大（HPC支持）、向其他部门扩展、与其他工具的联动进行选择，有助于长期的成本削减。
尖端技术
なった
詳しく
報告