Halbach阵列磁体设计与匀场的历史方法

《Journal of Magnetic Resonance Open》：Historical Methods for Designing and Shimming Halbach Array Magnets

　　

本文以课堂教学式教程呈现，提供推导细节与算例；研究人员对Klaus Halbach的二维束线磁体设计方法进行了解释与解读；研究人员将Robert Gluckstern关于有限长度磁体的分析以公共领域形式呈现；并对实际磁体材料与被动匀场进行了定量讨论。

Halbach阵列磁体设计与匀场的历史方法——论文解读

**研究背景与意义**
Halbach阵列（Halbach array）或Halbach环磁体（Halbach ring magnet）能够产生垂直于环轴的横向磁场，因其对永磁材料的高效利用、无铁结构带来的低重量以及内部场强高且外部杂散场小等优点，广泛用于需要相对均匀磁场的装置中。然而，该方案的主要实际弱点在于永磁块无法高精度制造：受粉末冶金工艺限制，即使两块磁体外部尺寸相同，其磁化强度也可能存在1-3%的差异，磁化方向存在±3°的偏差。这种缺陷会严重劣化磁场质量，必须通过测量分选、径向位置调整、附加小磁片等匀场（shimming）手段进行补偿。Klaus Halbach和Robert Gluckstern在数十年前发展出精密的解析方法，这些方法基于复变函数分析，能够给出精确的磁场表达式。但当代设计者往往直接依赖数值模拟，对这些解析结果及其物理内涵逐渐陌生。为了传承这些经典方法并指导高均匀度磁体设计，研究人员在《Journal of Magnetic Resonance Open》上发表了这篇教学式论文，系统梳理并扩展了Halbach的二维（2D）设计方法和Gluckstern的有限长度磁体分析，并通过算例展示了其在真实磁体误差评估与被动匀场中的应用。

**关键技术方法**
研究人员主要采用了以下技术方法：1) 复变函数方法：以复数表示二维空间位置与磁场分量，将磁场表示为复量B*=B_x?iB_y，利用积分变换导出永磁块的磁场公式；2) 多极场展开：将磁场积分中的被积函数展开为幂级数，得到多极场系数b_n，从而分析分块结构允许的谐波阶次；3) Gluckstern轴向间隙法：将有限长环磁体看作连续磁化分布，推导轴上磁场表达式，并通过将磁体分段、调整环间轴向间隙来消除低阶场误差；4) 随机误差模拟：对磁块的磁化强度与取向加入高斯分布随机误差，模拟真实材料的不一致性；5) 最小二乘拟合匀场：基于各磁块径向移动对磁场的影响图（field map），用线性最小二乘求解最优径向移动量，以最小化残余不均匀度。上述方法均为解析计算与数值拟合结合，不涉及具体实验样本队列。

**研究结果**
文章主体分为Halbach二维设计方法论、Gluckstern有限长度效应和真实磁体与匀场三大部分，下面分小节概述。

**2. Halbach的2D设计方法论**：该部分从单根线电流的磁场复表达式出发，推导了永磁块磁场的积分公式。利用稀土钴（SmCo）和钕铁硼（NdFeB）材料的线性B?H特性、易轴取向和相对磁导率接近1的性质，将磁化强度用复数B_r表示。通过对磁场积分进行分部积分，得到简洁的磁场公式B*(z₀)=1/(2π)∫∫B_r/(z₀?z)²dxdy。在此基础上，将积分核展开为z₀的幂级数，定义多极场系数b_n。对于环形连续磁化结构，若磁化方向随极角线性变化，则只有当旋转指数满足特定条件时对应的多极场才不为零；对于偶极子设计，磁化矢量绕环旋转两圈时，理论上只能在孔内产生均匀场。实际中磁体由有限个梯形块组成，研究者推导了每块磁体的贡献，并证明了分块结构只允许产生满足(1?n)/N_s为整数的谐波。例如8段偶极磁体，最低阶允许误差项为b₉，其相对强度与分段数及工作半径密切相关。当计划使用50%孔径时，N_s=4的磁体会产生约10%的四阶误差；而N_s=8可将误差降至1%以下。对更高均匀度（如ppm级）应用，需限制孔径利用率或增加段数。

**3. Gluckstern的有限长度方法**：Halbach未讨论有限长度效应。Gluckstern针对MRI用环偶极磁体，推导了有限长连续磁化环的轴上磁场表达式，并开发了轴向间隙方法。研究人员给出了轴上磁场的解析形式，并分析显示，场不均匀度对磁体长度更敏感，而非径向厚度。通过将磁体轴向分割为两个环并调整间隙，可以消除偶次场误差，例如将两个内径40 mm、外径60 mm、长度50 mm的环分开约6.5 mm即可使二阶场误差归零。增加环的数量可提供更多自由度以消除更高阶误差。

**4. 真实磁体与匀场**：研究人员以12块矩形磁块构成的偶极环为模型，计算了完美情况下的场图，并加入1%的强度误差和1°的方向误差进行模拟，发现峰峰值不均匀度约为13000 ppm，远超分段本身导致的误差。为了校正这种误差，他们研究了径向移动匀场。计算表明，单块磁体径向移动0.01个单位（内半径的1%）在r=0.5圆上产生的场变化约5000 ppm，足以抵消上述误差。通过最小二乘拟合各磁块移动量，残余峰峰值不均匀度可降至815 ppm，接近完美磁体的第一允许误差水平。进一步分析指出，若要将最终均匀度控制在50 ppm，径向位置控制精度需达到内半径的0.01%，例如对于内径40 mm的磁体，需4 μm的定位精度。该研究还评估了轴向间隙对机械误差的敏感性，并指出其他被动匀场方法（如按钮式匀场）和主动匀场虽可进一步改善均匀度，但已超出本文的历史方法范围。

**总结与讨论**
作者结合Halbach和Gluckstern的方法，提出了一套完整的设计工作流程：首先根据场强和尺寸要求确定环的内外半径，然后利用Gluckstern公式优化轴向分段与间隙，接着选择方位角分段数（如N_s=12）并调整段几何以补偿场强变化，最后预留径向移动空间以进行实测匀场。作者强调，在现代高精度磁体设计中，解析方法仍具有重要价值，因为永磁体本身的性质误差往往会抵消过度精细的优化，因此实际工程中，有效的匀场设计往往比磁体几何优化更为关键。

**研究结论**
（翻译自原文Conclusion）这些现在在NMR和MRI界流行的环形磁体的创始人，尤其是Klaus Halbach和Robert Gluckstern，对这些磁体有着深刻的理解，这些理解体现在他们的论文、专利和私人笔记中。这些方法发展于大规模计算机计算普及之前，依赖于已经不再常用的纸笔技术，因此他们的结果显得晦涩难懂。研究人员以教学形式呈现了他们的解析结果，希望现代设计者能够从他们的见解中受益。文中以足够的细节给出了应用实例，以指导那些希望将这些不熟悉的方法应用于其他感兴趣几何结构的人。这里开发的用于计算磁体内部任意位置磁场的方程，对于探索和理解永磁材料中不可避免的磁化误差的影响将特别有用。这些方程还将有助于制定针对这些误差的有效缓解策略。对于那些寻求改善设计均匀度的人来说，一个实际的教训是，永磁块的性质很可能会抵消广泛的优化努力，并破坏依赖微妙几何效应的设计。