永久电流模式解析
永久电流模式的理论基础
概要
老师！今天是永久电流模式解析的话题吧？那是什么东西呢？
超导环路的电阻接近零时电流持续很长时间的工作模式。接合部电阻的电流衰减评估。MRI主磁铁的稳定运行。
支配方程
离散化手法
这个方程在计算机上到底要怎么解呢？
使用有限元法（FEM）进行空间离散化。组装单元刚度矩阵，构造整体刚度方程。
矩阵求解算法
矩阵求解算法具体是怎么回事呢？
通过直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解联立方程组。对于大规模问题，带预处理的迭代法是有效的。
| 求解方法 | 分类 | 内存使用量 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小～中规模 |
| Cholesky分解 | 直接法（对称正定） | O(n²) | 小～中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模·非对称 |
| AMG预处理 | 预处理 | O(n) | 超大规模 |
也就是说，有限元法这个地方要是做得不到位，后面就要吃亏！我铭记在心！
商用工具中的实现
那么做永久电流模式解析用什么软件呢？
| 工具名称 | 开发者/现在 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| ANSYS Maxwell | ANSYS Inc. | .aedt, .maxwell |
| JMAG-Designer | JSOL Corporation | .jmag, .jproj |
供应商系统及产品整合的历程
各个软件的成长历程，有没有什么戏剧性的故事呢？
COMSOL Multiphysics
给我讲讲"COMSOL Multiphysics"吧！
1986年在瑞典成立。以MATLAB联动的FEMLAB开始，后来改名为COMSOL。擅长多物理场。
现在的所属：COMSOL AB
ANSYS Maxwell
给我讲讲"ANSYS Maxwell"吧！
Ansoft Maxwell。低频电磁场解析。2008年整合入ANSYS。
现在的所属：ANSYS Inc.
JMAG-Designer
JMAG具体是怎么回事呢？
由日本JSOL公司开发。电磁场分析工具，专注于电气设备设计。
现在的所属：JSOL Corporation
啊，我明白了！年在瑞典成立这样的结构我终于理解了。
文件格式与互操作性
不同软件间数据交换的时候有什么要注意的吗？
| 格式 | 扩展名 | 种类 | 概述 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | 符合ISO 10303的3D CAD数据交换格式。形状+PMI对应。 |
| IGES | .igs/.iges | 中立CAD | 早期的CAD数据交换规范。曲面数据兼容性存在问题。向STEP过渡。 |
| VTK | .vtk/.vtu | 可视化 | 可视化工具包格式。在ParaView等中使用。 |
在不同求解器间转换模型时，要素类型的对应关系、材料模型的兼容性、荷载和边界条件的表示差异需要注意。特别是高阶要素或特殊要素（内聚要素、用户定义要素等）多数情况下不能在求解器间直接转换。
看起来格式看似简单，其实里面藏着很深的东西啊。
实务上的注意事项
教科书里没有的"现场的智慧"有吗？
网格收敛性的确认、边界条件的妥当性验证、材料参数的灵敏度分析非常重要。
嗯，进展不错！实际上动手是最好的学习方法。有不懂的地方随时问我。
永久电流模式——"关掉开关磁场仍然持续"的超导特性
"永久电流模式（PCM）"是超导最引人瞩目的特性之一。从外部电源断开超导磁铁后，磁场仍持续——通过"持久电流开关（PCS）"对线圈两端短接，关闭外部电源时，超导环路内电流永不停止。实际上"磁通蠕变"会导致电流每小时衰减不到几ppm，但MRI磁铁可以维持多年以上的稳定磁场。PCM在MRI·NMR分光仪中采用，超过外部电源稳定度的磁场稳定性有助于医疗诊断·科学计测的精度。
永久电流模式的数值计算方法
数值手法的详细说明
具体用什么算法解永久电流模式解析呢？
哇～永久电流模式解析的话题，超级有意思！再给我讲讲。
离散化的制定
用形状函数 $N_i$ 近似未知量：
用这个式子来表示。
基础方程的离散形式
用这个式子来表示。
嗯～只看公式有点看不懂呢… 是什么意思啊？
连续体的支配方程离散化后，得到如下代数方程系统：
这里$[K]$是整体刚度矩阵（或等效系统矩阵），$\{u\}$是未知节点变量向量，$\{F\}$是外力向量。
啊，我明白了！连续体的支配方程这样的机制我终于理解了。
要素技术
"要素技术"听说过，但可能没理解透…
| 要素类型 | 阶数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低（剪力锁定） | 低 |
| 四面体2阶 | 二次 | 10 | 高 | 中 |
| 六面体1阶 | 线性 | 8 | 中 | 中 |
| 六面体2阶 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中～高 | 中 |
积分方案
积分方案具体是怎么回事呢？
听到这里，我终于理解为什么要素类型那么重要了！
收敛性和稳定性
收敛不了的时候，首先要检查什么？
收敛速度：二次要素以$O(h^2)$的阶数减小误差（光滑解的情况）
看起来网格细分很简单，其实深藏着很深的东西呢。
求解器设置的建议
具体用什么算法解永久电流模式解析呢？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数基准 |
| 预处理方法 | ILU(0) or AMG | 按问题规模而定 |
| 最大迭代次数 | 1000 | 不收敛时修改设置 |
| 内存模式 | In-core | 尽可能采用 |
边要素（Nedelec要素）
专为电磁场解析的要素。自动确保切线分量连续性，排除虚假模。3D高频解析的标准。
节点要素
用于标量势定式。在静磁场的标量势法或静电场解析中有效。
FEM vs BEM（边界要素法）
FEM：对非线性材料、非均质介质对应。BEM：能自然处理无限领域（开领域问题）。混合FEM-BEM也有效。
非线性收敛（磁气饱和）
用牛顿-拉弗逊法处理B-H曲线非线性。残差基准：$||R||/||R_0|| < 10^{-4}$是一般水平。
频域解析
通过时间谐波假设降为稳态问题。需要复数运算，但宽带特性用时域解析获取。
时间域的时间步长
需要最高频率分量的1/20以下的时间步长。隐式时间积分可以取较大步长，但要注意精度。
永久电流模式的实务应用
实践指南
老师，"实践指南"给我讲讲！
阐述永久电流模式解析的实务解析流程及注意事项。
解析流程
从最初的一步开始教教我！首先要做什么呢？
1. 预处理 (Pre-processing)
- 导入CAD数据并简化形状
- 定义材料特性
- 网格生成（要素类型、尺寸决定）
- 设置边界条件和荷载条件
2. 求解 (Solving)
- 求解器设置（解法、收敛基准、输出控制）
- 提交作业并执行计算
- 收敛监视
3. 后处理 (Post-processing)
- 结果可视化（位移、应力等物理量）
- 结果验证和妥当性确认
- 报告编写
网格生成的最佳实践
网格的好坏怎么判断呢？
要素品质指标
"要素品质指标"给我讲讲！
| 指标 | 理想值 | 允许范围 | 影响 |
|---|---|---|---|
| 纵横比 | 1.0 | < 5.0 | 精度下降 |
| 雅可比行列式比 | 1.0 | > 0.3 | 要素退化 |
| 翘曲 | 0° | < 15° | 精度下降 |
| 斜度 | 0° | < 45° | 收敛性恶化 |
| 锥形比 | 0 | < 0.5 | 精度下降 |
网格密度的决定
网格密度的决定具体是怎么回事呢？
边界条件的设置指南
边界条件，听说这个地方错了全部就完蛋…
啊，我明白了！过度约束警戒这样的机制我终于理解了。
商用工具的逐一实现步骤
各种各样的软件吧？各自的特点给我讲讲！
| 工具名称 | 开发者/现在 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| ANSYS Maxwell | ANSYS Inc. | .aedt, .maxwell |
| JMAG-Designer | JSOL Corporation | .jmag, .jproj |
COMSOL Multiphysics
"COMSOL Multiphysics"给我讲讲！
1986年在瑞典成立。以MATLAB联动的FEMLAB开始，后来改名为COMSOL。擅长多物理场。
现在的所属：COMSOL AB
ANSYS Maxwell
"ANSYS Maxwell"给我讲讲！
Ansoft Maxwell。低频电磁场解析。2008年整合入ANSYS。
现在的所属：ANSYS Inc.
老师的说明好容易理解！工具名字的疑惑解开了。
常见故障及对策
新手容易犯什么样的失误呢？想预先知道！
| 现象 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不良、不恰当的边界条件 | 改善网格、重新检查约束条件 |
| 应力异常大 | 应力奇点、网格依赖 | 避免奇点、局部网格细分 |
| 位移不现实 | 材料常数误差、单位系不统一 | 确认输入数据 |
| 计算时间过长 | 不必要的细分、低效解法 | 网格最优化、并行计算 |
品质保证检查清单
教科书里没有的"现场的智慧"有吗？
嗯，进展不错！实际上动手是最好的学习方法。有不懂的地方随时问我。
"磁场掉了"——MRI永久电流模式中的场漂移对策
MRI永久电流模式的"磁场漂移（Field Drift）"与图像品质的长期稳定性直接相关。漂移主因是①磁通蠕变（Anderson-Kim型，通常<1 ppm/hr）②机械振动引起的磁通移动③温度变化导致的Jc变化。临床MRI（1.5T/3T）要求磁场稳定性在0.1 ppm/hr以内。CAE中模拟永久电流模式的长期蠕变，通过线圈设计（电流余量·接头最优化）满足要求值。此外"有源屏蔽磁铁"设计中，屏蔽线圈的蠕变特性与主线圈的差异导致磁场分布随时间变化，需要进行对称蠕变设计。
永久电流模式的软件比较
商用工具比较
各种各样的软件吧？各自的特点给我讲讲！
对应永久电流模式解析的主要商用CAE工具功能比较及各产品的历史背景进行详述。
哇～永久电流模式解析的话题，超级有意思！再给我讲讲。
支持工具一览
那么做永久电流模式解析用什么软件呢？
| 工具名称 | 开发者/现在 | 主要文件格式 |
|---|---|---|
| COMSOL Multiphysics | COMSOL AB | .mph |
| ANSYS Maxwell | ANSYS Inc. | .aedt, .maxwell |
| JMAG-Designer | JSOL Corporation | .jmag, .jproj |
COMSOL Multiphysics
"COMSOL Multiphysics"给我讲讲！
1986年在瑞典成立。以MATLAB联动的FEMLAB开始，后来改名为COMSOL。擅长多物理场。
现在的所属：COMSOL AB
ANSYS Maxwell
"ANSYS Maxwell"给我讲讲！
Ansoft Maxwell。低频电磁场解析。2008年整合入ANSYS。
现在的所属：ANSYS Inc.
JMAG-Designer
JMAG具体是怎么回事呢？
由日本JSOL公司开发。电磁场分析工具，专注于电气设备设计。
现在的所属：JSOL Corporation
所以说，年在瑞典成立这样的结构才是这样啊。
功能比较矩阵
预算和时间都有限，综合最强是哪个？
| 功能 | COMSOL | Maxwell | JMAG |
|---|---|---|---|
| 基本功能 | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ |
| GPU对应 | △ | △ | ○ |
转换时的风险
转换时的风险具体是怎么回事呢？
啊，我明白了！不同工具间的模型这样的机制我终于理解了。
许可证形式
"许可证形式"听说过，但可能没理解透…
| 工具 | 许可证 | 特点 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 高价但有官方支持 |
| OpenFOAM | GPL | 免费但支持需付费 |
| COMSOL | 节点锁定/浮动 | 以模块单位购买 |
| Code_Aster | GPL | EDF开发的开源求解器 |
选择指南
最后到底该选哪个，判断基准给我教教？
在永久电流模式解析的工具选定中需要考虑以下几点：
嗯，进展不错！实际上动手是最好的学习方法。有不懂的地方随时问我。
永久电流模式解析工具——ANSYS Maxwell Transient vs COMSOL
永久电流模式的解析需要"超导线圈+电路（PCS）+热"的多物理场耦合，ANSYS Maxwell（Transient Solver+Circuit）和COMSOL（AC/DC+Heat Transfer）是代表性工具。Maxwell通过外部Simplorer（ANSYS Twin Builder）可以进行耦合，在持久电流开关的电气行为进行线路级精密解析，能够仿真整个充放电过程。COMSOL能够用自定义PDE实现超导的E-J特性，适合包括蠕变的长时间行为计算。MRI磁铁制造商（日立、东芝、Bruker）采用的是在MAXWELL/COMSOL用自主设计代码验证的二阶段方法。
永久电流模式的先端研究
先端话题与研究动向
永久电流模式解析领域，今后会怎么发展呢？
来看永久电流模式解析的最新研究动向与先进手法。
哇～永久电流模式解析的话题，超级有意思！再给我讲讲。
最新的数值手法
接下来是最新数值手法的话题吧。什么内容呢？
嗯～只看公式有点看不懂呢…是什么意思啊？
高性能计算 (HPC) 的应对
| 并行化方法 | 概述 | 支持求解器 |
|---|---|---|
| MPI (领域分割) | 分布式内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 多数求解器 |
| GPU (CUDA/OpenCL) | GPGPU活用。特别是陈述法有效 | LS-DYNA, Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
永久电流模式的故障应对
故障排除
哇～永久电流模式解析的话题，超级有意思！再给我讲讲。
常见错误及对策
老师也因为永久电流模式解析通宵调试过吗？（笑）
1. 收敛失败
收敛失败具体是怎么回事呢？
现象：求解器未在指定迭代次数内收敛而异常终止
可能的原因：
- 网格品质不足（过度扭曲的要素）
- 材料参数设置不当
- 不恰当的初值条件
- 非线性过强（荷载步不足）
对策：
- 实施网格品质检查（纵横比、雅可比行列式）
- 确认材料参数的单位系
- 荷载分成多个步长（增加子步数）
- 放宽收敛判定基准（但要注意精度）
也就是说，收敛失败的地方要是做得不到位，后面就要吃亏吧。铭记在心！
2. 非物理的结果
接下来是非物理的结果的话题吧。什么内容呢？
现象：应力/位移/温度等显示非现实值
可能的原因：
- 边界条件设置误差
- 单位系混用（SI单位与工程单位混淆）
- 不适当的要素类型选择
- 存在应力奇点
对策：
- 确认反力合计（力的平衡）
- 确认单位系一致性
- 重新检查要素类型的适切性
- 奇点排除或子建模
前辈说"至少要把收敛失败做好"的意思我终于明白了。
3. 计算时间超过
计算时间超过具体是怎么回事呢？
现象：计算耗时超过预计的好几倍
对策：
- 优化网格粗密分布
- 活用对称性（1/2、1/4模型）
- 优化求解器设置（迭代法、预处理选择）
- 活用并行计算
4. 内存不足
"内存不足"给我讲讲！
现象：Out of Memory 错误
前辈说"至少要把收敛失败做好"的意思我终于明白了。
对策：
- 使用核外解法
- 减少网格规模
- 确认64位版求解器的使用
- 增加内存分配
哇～收敛失败的话题，超级有意思！再给我讲讲。
Nastran代表性错误
代表性错误具体是怎么回事呢？
Abaqus代表性错误
"