一种磁共振Halbach型匀场线圈的设计方法发明专利.docxVIP

一种磁共振Halbach型匀场线圈的设计方法 技术领域 本发明涉及属于核磁共振匀场线圈技术领域，具体涉及一种磁共振Halbach型匀场线圈的设计方法。 背景技术 Halbach型（又称为海尔贝克，哈尔巴赫等）磁体是由20世纪80 年代美国物理学家K.Halbach提出了一种新型磁体设计。Halbach型磁体由多个具有不同的充磁方向的磁块按照一定的规律组合而成，能利用较少的磁材产生较均匀的单边磁场，如图1所示。与常规磁体设计相比，在磁体的有效空间和场强相当情况下，Halbach型磁体具有体积更小所消耗的材料更少等优点，在小型化便携化核磁共振（NMR）、磁共振成像（MRI）以及电子顺磁共振（EPR）等磁共振仪器领域有较大的应用前景。但Halbach型磁体也存在较大的主磁场不均匀的问题，而NMR仪器、MRI仪器以及EPR仪器对磁体主磁场均匀度有较高要求，因此应用于磁共振仪器领域的Halbach型磁体，需要通过匀场技术提高其主磁场均匀度。提高主磁场均匀度的匀场技术可分为有源匀场和无缘匀场两种，其中有源匀场是通过特定电流分布的匀场线圈来实现的。 本发明专利为了解决Halbach型磁体系统主磁场不均匀性的问题，提出了一种磁共振Halbach型匀场线圈的设计方法，使用该方法设计的匀场线圈可以消除Halbach型磁体系统主磁场各阶不均匀磁场分量，提高主磁场的均匀性（如图1所示）。 发明内容 本发明的目的是针对现有技术存在的上述问题，提供一种磁共振Halbach型匀场线圈的设计方法，能消除Halbach型磁体系统主磁场各阶不均匀磁场分量，提高主磁场的均匀性。 为了解决上述技术问题，本发明所采用的技术方案是： 一种磁共振Halbach型匀场线圈的设计方法，包括以下步骤： 步骤1、首先对Halbach匀场线圈所在的线圈圆柱面进行网格划分获得各个网格节点； 步骤2、在目标圆球面上选取目标点； 步骤3、建立网格节点电流密度流函数?和网格节点电流密度的有限差分关系； 步骤4、设置目标点磁场强度； 步骤5、建立网格节点电流密度流函数?和目标点磁场强度z轴方向分量之间的方程组； 步骤6、构造功耗和局部均匀分布组合约束函数H，根据功耗和局部均匀分布组合约束函数获得约束矩阵 L，利用Tikhonov正则化方法求解步骤5中的网格节点电流密度流函数?和目标点磁场强度z轴方向分量之间的方程组，得到网格节点电流密度流函数?的具体数值； 其中，R1为线圈圆柱面的底面半径，?为线圈圆柱面的高度, R1为线圈圆柱面的底面半径，?为绕圆周的角度方向网格电流密度，?为z轴方向网格电流密度，zs为网格节点的z轴高度变量，?为轴向划分线的方位角变量； 步骤7、根据网格节点电流密度流函数?在对应的线圈圆柱面上的分布的等势线做等差值排布即得到匀场线圈的绕线模式。 如上所述的步骤1中： 线圈圆柱面由多个轴向划分线和多个周向划分圆划分，各个轴向划分线位于线圈圆柱面上且与线圈圆柱面的中心轴平行，各个周向划分圆位于线圈圆柱面上且垂直于线圈圆柱面的中心轴，轴向划分线的编号用i表示；周向划分圆的编号用j表示。 如上所述的所述的步骤3中，网格节点电流密度包括绕圆周的角度方向网格电流密度?和z轴方向网格电流密度?，网格节点电流密度流函数?和网格节点电流密度的有限差分关系为： 其中，网格节点?表示网格节点（i,j）和网格节点（i,j+1）之间的中间节点；网格节点?表示网格节点（i,j）和网格节点（i+1,j）之间的中间节点；?表示网格节点（i,j）处电流密度流函数，?表示网格节点（i+1，j）处电流密度流函数，?表示网格节点（i，j+1）处电流密度流函数，?为网格节点?所在的线圈圆柱面的半径，?为相邻轴向划分线的方位角之间的角度差，?为相邻周向划分圆之间的间距。 如上所述的步骤5包括以下步骤： 步骤5.1、建立目标点磁场强度z轴方向分量?的控制方程的离散方程； 其中，zs为网格节点的z轴高度变量，?f= 1，2…N，N为目标点的总数目，?为目标点在直角坐标系下的坐标值，?为真空磁导率， 步骤5.2、根据目标点磁场强度z轴方向分量?的控制方程的离散方程，建立网格节点电流密度流函数?和目标点磁场强度z轴方向分量之间的方程组： 为第一个网格节点到第M个网格节点电流密度流函数，?，……，?为网格节点电流密度流函数的系数，M为线圈圆柱面网格划分得到的所有网格节点的总数目，N为目标点的总数目，?为第一个目标点磁场强度z轴方向分量到第N个目标点磁场强度z轴方向分量。 本发明相对于现有技术，具有以下有益效果：本发明简单有效，为设计Halbach型匀场线圈提供了新的思路，设计的Halbach型匀场线圈可有效消除Halbach型磁体系统由于生产制造误差而引入的主磁场各阶不均匀磁场分量，进而提高Halbach型磁体系统主磁场的