2024年5月4日发(作者：root一键权限获取)

第４２卷增刊　

２００８年９月　

原子能科学技术　

Ｖｏ１．４２，Ｓｕｐｐ１．　

Ｓｅｐ．２００８　

Ａｔｏｍｉｃ　Ｅｎｅｒｇｙ　Ｓｃｉｅｎｃｅ　ａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ　

ＦＬＵＫＡ在永久磁铁辐射退磁实验模拟中的应用　

邱　睿　，李君利　，毕　垒　，李文茜　。　，周婧劫　

（１．清华大学工程物理系，北京　１０００８４；２．粒子技术与辐射成像教育部重点实验室（清华大学），北京　１０００８４；　

３．高能辐射成像国防重点学科实验室，北京　１０００８４）　

摘要：利用ＦＬＵＫＡ程序，对韩国浦项加速器实验室的２．５　ＧｅＶ电子直线加速器辐射退磁实验系统进行　

蒙特卡罗模拟，得到了永久磁铁处各种粒子（包括光子、电子、中子及质子）的能谱以及注量的空间分布。　

根据实验中的束流强度，估算出永久磁铁所受的吸收剂量率和ｌ　ＭｅＶ等效中子注量率分别为　

１８．７　Ｇｙ／ｓ和６．３７×１０　ｃｍ　・Ｓ～。　

关键词：蒙特卡罗方法；ＦＬＵＫＡ程序；辐射退磁；电子直线加速器　

中图分类号：ＴＬ９９　文献标志码：Ａ　文章编号：１０００—６９３１（２００８）ｓｏ一０３７５—０５　

Ａｐｐｌｉｃａｔｉｏｎ　ｏｆ　ＦＬＵＫＡ　ｉｎ　Ｓｉｍｕｌａｔｉｏｎ　ｏｆ　Ｒａｄｉａｔｉｏｎ—Ｉｎｄｕｃｅｄ　

Ｄｅｍａｇｎｅｔｉｚａｔｉｏｎ　Ｅｘｐｅｒｉｍｅｎｔ　ｏｆ　Ｐｅｒｍａｎｅｎｔ　Ｍａｇｎｅｔ　

ＱＩＵ　Ｒｕｉ　＇　，ＬＩ　Ｊｕｎ—ｌｉ　，　，ＢＩ　Ｌｅｉ　＇。　，ＬＩ　Ｗｅｎ－ｑｉａｎ　，。　，ＺＨＯＵ　Ｊｉｎｇ—ｊｉｅ　，　，。　

（１．Ｄｅｐａｒｔｍｅｎｔ　ｏｆ　Ｐｈｙｓｉｃｓ　Ｅｎｇｉｎｅｅｒｉｎｇ，Ｔｓｉｎｇｈｕａ　Ｕｎｉｖｅｒｓｉｔｙ，Ｂｅｉｊｉｎｇ　１　０００８４，Ｃｈｉｎａ；　

２．Ｋｅｙ　Ｌａｂｏｒａｔｏｒｙ　ｏｆ　Ｐａｒｔｉｃｌｅ＆Ｒａｄｉａｔｉｏｎ　Ｉｍａｇｉｎｇ（Ｔｓｉｎｇｈｕａ　Ｕｎｉｖｅｒｓｉｔｙ），Ｂｅｉｊｉｎｇ　１０００８４，Ｃｈｉｎａ；　

３．Ｍｉｎｉｓｔｒｙ　ｏｆ　Ｅｄｕｃａｔｉｏｎ　Ｋｅｙ　Ｌａｂｏｒａｔｏｒｙ　ｏｆ　Ｈｉｇｈ　Ｅｎｅｒｇｙ　Ｒａｄｉａｔｉｏｎ　Ｉｍａｇｉｎｇ　Ｆｕｎｄａｍｅｎｔａｌ　Ｓｃｉｅｎｃｅ　

ｆｏｒ　Ｎａｔｉｏｎａｌ　Ｄｅｆｅｒｉｓｅ，Ｂｅｉｊｉｎｇ　１０００８４，Ｃｈｉｎａ）　

Ａｂｓｔｒａｃｔ：Ｍｏｎｔｅ—Ｃａｒｌｏ　ｓｉｍｕｌａｔｉｏｎ　ｏｆ　ｔｈｅ　ｒａｄｉａｔｉｏｎ—ｉｎｄｕｃｅｄ　ｄｅｍａｇｎｅｔｉｚａｔｉｏｎ　ｅｘｐｅｒｉｍｅｎｔ　

ｓｙｓｔｅｍ　ｗｉｔｈ　２．５　ＧｅＶ　ｅｌｅｃｔｒｏｎ　ｌｉｎｅａｒ　ａｃｃｅｌｅｒａｔｏｒ　ａｔ　Ｋｏｒｅａ　Ｐｏｈａｎｇ　Ａｃｃｅｌｅｒａｔｏｒ　Ｌａｂｏｒａｔｏ—　

ｒｙ（ＰＡＬ）ｗａｓ　ｃａｒｒｉｅｄ　ｏｕｔ　ｗｉｔｈ　ＦＬＵＫＡ　ｃｏｄｅ．Ｔｈｅ　ｓｐｅｃｔｒａ　ｏｆ　ｖａｒｉｏｕｓ　ｐａｒｔｉｃｌｅｓ，ｉｎｃｌｕ—　

ｄｉｎｇ　ｐｈｏｔｏｎ，ｅｌｅｃｔｒｏｎ，ｎｅｕｔｒｏｎ　ａｎｄ　ｐｒｏｔｏｎ，ａｎｄ　ｓｐａｔｉａｌ　ｄｉｓｔｒｉｂｕｔｉｏｎ　ｏｆ　ｆｌｕｅｎｃｅ　ｗｅｒｅ　

ｄｅｒｉｖｅｄ．Ｆｕｒｔｈｅｒｍｏｒｅ　ｔｈｅ　ａｂｓｏｒｂｅｄ　ｄｏｓｅ　ｒａｔｅ　ａｎｄ　１　ＭｅＶ　ｅｑｕｉｖａｌｅｎｔ　ｎｅｕｔｒｏｎ　ｆｌｕｅｎｃｅ　ｏｆ　

ｐｅｒｍａｎｅｎｔ　ｍａｇｎｅｔ　ｗｅｒｅ　ｅｓｔｉｍａｔｅｄ　ａｓ　１８．７　Ｇｙ／ｓ　ａｎｄ　６．３７×１０　ｃｍ一　・ｓ一　ａｃｃｏｒｄｉｎｇ　ｔｏ　

ｔｈｅ　ｂｅａｍ　ｉｎｔｅｎｓｉｔｙ　ｉｎ　ｔｈｅ　ｅｘｐｅｒｉｍｅｎｔ．　

Ｋｅｙ　ｗｏｒｄｓ：Ｍｏｎｔｅ—Ｃａｒｌｏ　ｍｅｔｈｏｄ；ＦＬＵＫＡ　ｃｏｄｅ；ｒａｄｉａｔｉｏｎ　ｄｅｍａｇｎｅｔｉｚａｔｉｏｎ；ｅｌｅｃｔｒｏｎ　

ａｃｃｅｌｅｒａｔｏｒ　

Ｎｄ　Ｆｅ　Ｂ型永久磁铁由于其优异的磁性　

插入件中　引。插入件是第三代同步辐射装置　

和第四代光源最重要的组成部分之一，由一组　

能被广泛应用于加速器领域，其中，非常重要的　

一

个应用是用于同步辐射装置和第四代光源的　

极性交替按周期排列的永久磁铁组成嘲。　

收稿日期：２００８－０６—３０；修回日期：２００８—０８－２０　

基金项目：国家自然科学基金资助项目（１０７７５０８４）；高能质子加速器辐射防护关键问题研究资助项目　

作者简介：邱睿（１９８１一），女，湖北枣阳人，博士后，辐射防护专业　

３７６　

Ｎｄ　Ｆｅ　Ｂ型永久磁铁具有较强的辐射敏感性，　

而且长期处于高能电子加速器所产生的混合强　

辐射环境中，极有可能遭受辐射损伤而导致退　

磁，从而对其应用的系统造成影响甚至导致严　

重的后果　］。永久磁铁的辐射退磁效应是同　

步辐射装置以及第四代光源的一个关键问题。　

此外，永久磁铁的辐射退磁效应研究在空间科　

学研究、核能利用开发、磁约束受控热核聚变研　

究以及未来新概念武器和空间对抗等众多领域　

均有强烈的需求背景和重要意义　。］。　

从２Ｏ世纪８０年代开始，国内外的研究者们　

开始用不同的放射源对永久磁铁的辐射退磁效　

应进行研究¨８　¨。第三代同步辐射装置和第四　

代光源的电子加速器均处于ＧｅＶ量级的高能　

段，在这一能量段的电子加速器领域，永久磁铁　

的辐照效应研究却少有报道。因此，清华大学工　

程物理系和韩国浦项加速器实验室一起合作，共　

同开展高能电子加速器中永久磁铁辐射退磁效　

应的研究，在２．５　ＧｅＶ电子直线加速器上进行了　

永久磁铁辐射退磁实验。由于永久磁铁所处的　

辐射场非常强，永久磁铁所在位置的能谱、吸收　

剂量率以及中子注量率等辐射物理量无法直接　

测量得到，因此，本工作利用ＦＬＵＫＡ程序，对该　

辐射退磁实验系统进行蒙特卡罗模拟。　

ｌ　永久磁铁高能电子加速器辐射退磁　

实验　

１．１　实验装置　

永久磁铁辐射损伤实验在韩国浦项加速器　

实验室２．５　ＧｅＶ电子直线加速器的隧道末端　

实验区进行。２．５　ＧｅＶ电子通过束流位置监测　

器ＢＰＭ（ｂｅａｍ　ｐｏｓｉｔｉｏｎ　ｍｏｎｉｔｏｒ）入射在靶上，　

Ｎｄ　Ｆｅ—Ｂ型永久磁铁置于靶后８　ｍｍ处铝制的　

支架上，电子束的几何中一１，正对磁铁的几何中　

心（图１）。在电子入射前、后，分别用数字特斯　

拉计（ｄｉｇｉｔａｌ　ｔｅｓｌａｍｅｔｅｒ）沿ｙ轴测量磁铁上方　

的磁感应强度。　

１．２　辐射物理分析　

电子损失在铜块上会产生一系列次级粒　

子。高能电子引发电磁簇射，从而产生大量且　

能谱相当宽的韧致辐射光子及次级电子。韧致　

辐射光子除了可以通过光电效应、康普顿效应　

和电子对效应产生次级电子，还会因轫致辐射　

原子能科学技术　第４２卷　

Ｂ　

图１永久磁铁辐射退磁实验辐照装置示意图　

Ｆｉｇ．１　Ｓｃｈｅｍｅ　ｏｆ　ｉｒｒａｄｉａｔｉｏｎ　ｄｅｖｉｃｅｓ　

ｉｎ　ｒａｄｉａｔｉｏｎ　ｄｅｍａｇｎｅｔｉｚａｔｉｏｎ　ｅｘｐｅｒｉｍｅｎｔ　

光子能量较高，达到或超过了光核反应的阈　

能ｌ＿ｌ　，产生中子、质子、ａ粒子等。　

因此，高能电子加速器插入件中永久磁铁　

所处的辐射环境是一个非常复杂的混合场，包　

括大量能谱分布非常宽的光子及电子，还有一　

定数量且能谱也相当宽的中子及质子。由于永　

久磁铁所处的辐射场非常强，永久磁铁所在位　

置的能谱、吸收剂量率以及中子注量率无法直　

接测量得到，本工作利用ＦＬＵＫＡ程序口　对该　

辐射退磁实验系统进行了蒙特卡罗模拟。　

２　ＦＬＵＫＡ模拟计算　

高能电子加速器辐射场模拟中所采用的几　

何模型及材料设定基于永久磁铁辐射退磁实验　

装置的真实结构和材料。靶的几何尺寸为　

２５　ｍｍ×１００　ｍｍ×４０　ｍｍ，材料为Ｃｕ，密度为　

８．９６　ｇ／ｃｍ。。永久磁铁的材料组成设为　

Ｎｄ　Ｆｅ　Ｂ，也就是只考虑Ｎｄ—Ｆｅ－Ｂ型永磁材料　

的基本相，几何尺寸为１２　ｍｍ×４６　ｍｍ×　

８　ｍｍ，密度为７．６　ｇ／ｃｍ。。　

入射电子束的相关参数采用韩国浦项加速器　

实验室２．５　ＧｅＶ电子直线加速器的参数。电子束　

能量为２．５　ＧｅＶ。电子束的脉冲形状在ｚ、　方向　

均呈高斯分布，半高宽分别为２．５和７．５　ｒｎｉｎ。　

３　计算结果　

３．１　磁铁处粒子能谱　

定义永久磁铁所处的辐射场，最重要的信　

息就是永久磁铁所在区域的粒子类型及其能　

量、注量以及空间分布等。因此，利用ＦＩ　ＵＫＡ　

程序计算了Ｓ１表面上各种粒子的能谱，包括光　

子、电子、质子和中子（图２）。图２中，计算结　

果归一到１个初始入射电子的贡献。　

增刊　邱睿等：ＦＬＵＫＡ在永久磁铁辐射退磁实验模拟中的应用　３７７　

不同的发展趋势：光子和电子数目不断增加，在　

永久磁铁的右边界达到永久磁铁内部的最大　

吕　

ｏ　

值；中子和质子的注量均是先增大随后降低，在　

永久磁铁的中心区域达到最大值。这清楚地体　

现了电磁簇射在永久磁铁内部发展的过程。　

３．３　吸收剂量　

：　

＞　

０　

一　

＼　

勰　

利用ＦＬＵＫＡ程序对永久磁铁内的能量　

沉积及其空间分布进行了模拟计算。一个初始　

电子在永久磁铁内的总能量沉积为２．６９×　

图２磁铁表面Ｓ１上各种粒子的能谱　

Ｆｉｇ．２　Ｓｐｅｃｔｒａ　ｏｆ　ｖａｒｉｏｕｓ　ｐａｒｔｉｃｌｅｓ　

ａｔ　ｍａｇｎｅｔ　ｓｕｒｆａｃｅ　Ｓ１　

１０　ＧｅＶ／ｃｍ。，转换为吸收剂量即为５．６７×　

１０　。Ｇｙ。实验过程中，束流强度约为３．３×　

１０　Ｓ一，因此，可以估算出永久磁铁所受的剂　

量率为１８．７　Ｇｙ／ｓ。　

图４示出３个不同平面内能量沉积的二维　

Ｍ　

■——光子；●——质子；▲——中子；ｖ——电子　

由图２可见，在各个能段，光子的注量最　

空间分布。这里，沉积能量归一到１个初始入　

射电子的贡献。　

在ｚ轴方向上，能量沉积量逐渐增大，在　

一

高，然后依次为电子、中子和质子。１００　ＭｅＶ　

以上中子和质子能谱非常接近，应为准氘核过　

程在光中子的产生机制中占主导地位的结果。　

３．２　磁铁处粒子注量空间分布　

８　ｍｍ处达到最大值，这与电磁簇射的发展　

过程相一致。在ｘ轴和Ｙ轴方向上，能量沉积　

通过模拟各种粒子在永久磁铁内部的输　

运过程，得到各种粒子的注量在磁铁内部的　

空间分布。图３为永久磁铁内部各种粒子注　

量在ｙ－ｚ平面（ｚ一０平面）内的空间分布。粒　

子注量单位为ｃｍ

电子的贡献。　

由图３可见，各种粒子的注量变化呈现出　

注量／ｃｍ　

４５　

１２Ｏ　

４０　４０　

分别沿ｚ一０和Ｙ一０平面呈对称分布。与粒　

子的注量分布类似，在ｃｒ－ｙ平面上，能量沉积主　

要集中在磁铁中心的一个椭圆形区域。　

３．４　１　ＭｅＶ等效中子注量　

ＡＳＴＭ（Ａｍｅｒｉｅａｎ　Ｓｏｃｉｅｔｙ　ｆｏｒ　Ｔｅｓｔｉｎｇ　ａｎｄ　

ｌ　

８　

，即归一到１个初始入射　

Ｍａｔｅｒｉａｌｓ）制定了标准ＡＳＴＭ　Ｅ７７２—９４，可将　

具有不同能谱的中子注量折算到１　ＭｅＶ等效　

注量／ｃｍ～　

０．Ｏ３５　

注量／ｃｍ　

４５　

１Ｏ４｝主量／ｃｍ一　

０．０３０　

３５　

３Ｏ　

８Ｏ　

１ＯＯ　３５　

１０　３０　

０．０２５　

ｌ０　

１５　

３０　

ｌ　２５　

ｉ　

２０　

＼

ｌ　２５　

　

８差２５　

６　

２０　

０．０２ｏ　

ｌ　２０　

＾２５　

３０　

６０　２Ｏ　

０．Ｏ１　５　

１５　１　５　

４Ｏ　１　５　

４　

１０　

０．０１０　

１０　

３５　

４０　

１Ｏ　

２０　

５　５　

Ｏ．００５　

４５　

２　４　６　８　２　４　６　８　２　４　６　８　２　４　６　８　

ｚ｝ｍｍ　ｚ｜ｍｍ　２｛ｆｉｌｍ　ｚｆｍｍ　

图３各种粒子在永久磁铁ｙ－ｚ平面（　—Ｏ）内的空间分布　

Ｆｉｇ．３　Ｓｐａｔｉａｌ　ｄｉｓｔｒｉｂｕｔｉｏｎｓ　ｏｆ　ｐａｒｔｉｃｌｅｓ　ｉｎ　ｙ－ｚ（ｚ—Ｏ）ｐｌａｎｅ　ｏｆ　ｐｅｒｍａｎｅｎｔ　ｍａｇｎｅｔ　

ａ——光子．ｂ——电子；ｃ——中子；ｄ——质子　

３７８　原子能科学技术　第４２卷　

中子注量ｌ＿】　。由于高能电子加速器所产生的　

中子能谱非常宽（图２），相对于总的中子注量，　

对质子和中子分别造成的１　ＭｅＶ等效中　

子注量进行求和，得到总的１　ＭｅＶ等效中子注　

量为１．９３×１０　ｃｍ＿。。根据实验过程中束流　

１　ＭｅＶ等效中子注量能够更好地用于衡量辐　

射效应，其定义公式为：　

ｒ∞　

强度约为３．３×ｌＯ加Ｓ＿。，可估算出永久磁铁所　

受的１　ＭｅＶ等效中子注量率为６．３７×　

（１）　——一

１　ＭｅＶ　

　ｌ（Ｅ）Ｄ（Ｅ）ｄＥ　

　ｑ，１　ＭｅＶ一　—　

１０　ｃｍ　・ｓ。同时可看到，中子对总的　

其中：　（Ｅ）为中子的能谱；Ｄ（Ｅ）是能量为Ｅ的　

中子对硅材料位移损伤的比释动能因子；Ｄ　、，　

为１　ＭｅＶ中子对硅材料位移损伤的比释动能　

因子。ＡＳＴＭ　Ｅ７７２—９４标准给出了１Ｏ　。～　

２０　ＭｅＶ中子的Ｄ（Ｅ），来源于Ｐ．Ｊ．Ｇｒｉｆｆｉｎ　

１　ＭｅＶ等效中子注量的贡献占大多数，达到　

８５．４　。　

另外，根据１　ＭｅＶ等效中子注量的定义　

（式（１））可以推出，能量为Ｅ的粒子所造成的　

１　ＭｅＶ等效中子注量为：　

ｑｌｌ　

ｑ）（Ｅ）Ｄ（Ｅ）ｄＥ　

Ｅ）一生　（２）　

等　获得的数据。对于质子，利用文献［１６］，　

也可将质子的注量等效为１　ＭｅＶ中子注量。　

利用ＦＬＵＫＡ计算得到的质子和中子的　

能谱以及质子和中子的位移损伤比释动能因　

根据中子和质子的能谱以及它们相应的位　

移损伤比释动能因子，计算得到了具有不同能　

量Ｅ的粒子所造成的１　ＭｅＶ等效中子注量　

（图５）。　

子，根据式（１）可以计算得到它们所对应的　

１　ＭｅＶ等效中子注量（表１）。计算结果归一到　

１个初始入射电子的贡献。　

表１　不同粒子造成的１　ＭｅＶ等效中子注量　

Ｔａｂｌｅ　１　１　ＭｅＶ　ｅｑｕｉｖａｌｅｎｔ　ｎｅｕｔｒｏｎ　ｆｌｕｅｎｃｅ　

由图５可见，在２．５　ＧｅＶ电子加速器的辐　

射场中，由不同能量的中子所造成的１　ＭｅＶ等　

效中子注量差别较大，相差４个数量级甚至更　

多。其中，０．２～１０　ＭｅＶ能量段的中子造成的　

ｃｏｎｔｒｉｂｕｔｅｄ　ｂｙ　ｄｉｆｆｅｒｅｎｔ　ｐａｒｔｉｃｌｅｓ　

１　ＭｅＶ等效中子注量较大。相对于中子，不同　

能量的质子造成的１　ＭｅＶ等效中子注量差别　

较小，尤其是１０　ＭｅＶ以下的质子；２　ＭｅＶ以　

下，质子造成的１　ＭｅＶ等效中子注量变化比较　

小，２　ＭｅＶ以上不断增加在１００　ＭｅＶ左右达　

到了最大值，并形成一个峰，之后逐渐降低。　

沉积能量／（ＧｅＶ・ｃｍ。）　

１２　

１　０　

１０　

ｇ　８　

６　

４　

沉积能量／（ＧｃＶ・ｃｍ　）　

２　

２　４　６　８　

５　ｌＯ　１５　２０　２５　

’｜ｍｍ　

３０　３５　４０　４５　

ｚ／ｍｍ　

沉积能量／（ＧｅＶ・ｃｍ　）　

２　

ｌ　４　

６　

８　

５　１０　１５　２０　２５　．３０　３５　４０　４５　

ｙ／ｍｍ　

图４永久磁铁内部能量沉积分布　

Ｆｉｇ．４　Ｄｉｓｔｒｉｂｕｔｉｏｎ　ｏｆ　ｅｎｅｒｇｙ　ｄｅｐｏｓｉｔｉｏｎ　ｉｎ　ｐｅｒｍａｎｅｎｔ　ｍａｇｎｅｔ　

ａ——　ｚ平面．ｂ——　平面；ｃ—　ｚ平面　

增刊　邱睿等：ＦＬＵＫＡ在永久磁铁辐射退磁实验模拟中的应用　３７９　

目　

ＥｆＧｅＶ　

图５不同能量粒子造成的１　ＭｅＶ等效中子注量　

Ｆｉｇ．５　１　ＭｅＶ　ｅｑｕｉｖａｌｅｎｔ　ｎｅｕｔｒｏｎ　ｆｌｕｅｎｃｅ　

ｃｏｎｔｒｉｂｕｔｅｄ　ｂｙ　ｐａｒｔｉｃｌｅｓ　ｗｉｔｈ　ｄｉｆｆｅｒｅｎｔ　ｅｎｅｒｇｙ　

■——质子；▲——中子　

４　结论　

本工作利用ＦＬＵＫＡ程序，对该辐射退磁　

实验的系统进行了蒙特卡罗模拟，对永久磁铁　

所在的辐射场及相关辐射物理量进行了模拟计　

算和分析。通过模拟高能电子在靶以及永久磁　

铁中的输运过程，计算得到各种粒子包括光子、　

电子、中子及质子的能谱以及通量的空间分布，　

根据所得结果分析了电磁簇射在永久磁铁内的　

发展情况。计算了１个初始电子所造成的永久　

磁铁的吸收剂量和１　ＭｅＶ等效中子注量，根据　

实验中的束流强度，估算出吸收剂量率和　

１　ＭｅＶ等效中子注量率分别为１８．７　Ｇｙ／ｓ和　

６．３７×１０　ｃｒｎ　・ｓ－。。另外，对不同能量的粒　

子所造成的１　ＭｅＶ等效中子注量进行了计算　

和分析，发现中子对总ｌ　ＭｅＶ等效中子注量的　

贡献达到８５．４　，其中，０．２～１０　ＭｅＶ能量段　

的中子造成的１　ＭｅＶ等效中子注量较大。　

参考文献：　

［１］周寿增，董清飞．超强永磁体——稀土铁系永磁　

材料［Ｍ］．２版．北京：冶金工业出版社，２００４：　

３４—５８．　

［２３　

ＢＲ０ＷＮ　Ｇ，ＨＡＬＢＡＣＨ　Ｋ，ＨＡＲＲＩＳ　Ｊ，ｅｔ　ａ１．　

Ｗｉｇｇｅｒ　ａｎｄ　ｕｎｕｌａｔｏｒ　ｍａｇｎｅｔｓ：Ａ　ｒｅｖｉｅｗ［Ｊ］．Ｎｕ—　

ｃｌ　Ｉｎｓｔｒｕｍ　Ｍｅｔｈｏｄｓ，１９８３，２０８：６５－７７．　

［３］　

马礼敦，杨福家．同步辐射应用概论ｌ－Ｍ］．上海：　

复旦大学出版社，２００１．　

［４］　

ＨＡＲＤＥＫＯＰＦ　Ｒ．Ｐｅｒｍａｎｅｎｔ　ｍａｇｎｅｔ　ｒａｄｉａｔｉｏｎ　

ｒｅｓｉｓｔａｎｃｅ　ｉｎ　ｔｈｅ　ＳＮＳ　Ｌｉｎａｃ［Ｒ］．Ｔｅｎｎｅｓｓｅｅ，ＵＳ：　

Ｓｐａｌｌａｔｉｏｎ　Ｎｅｕｔｒｏｎ　Ｓｏｕｒｃｅ，２００１　

［５］ＣＨＡＶＡＮＮＥ　Ｊ，ＥＬＬＥＡＵＭＥ　Ｐ，ｖａｎ　ＶＡＥＲＥ—　

ＮＢＥＲＧＨ　Ｐ．Ｐａｒｔｉａｌ　ｄｅｍａｇｎｅｔｉｚａｔｉｏｎ　ｏｆ　ＩＤ６　ａｎｄ　

ｄｏｓｅ　ｍｅａｓｕｒｅｍｅｎｔｓ　ｏｎ　ｃｅｒｔａｉｎ　ＩＤ，Ｍａｃｈｉｎｅ　Ｔｅｃｈ—　

ｎｉｃａｌ　Ｎｏｔｅ　１－１９９６／ＩＤ［Ｒ］．Ｆｒａｎｃｅ：Ｅｕｒｏｐｅａｎ　

Ｓｙｎｃｈｒｏｔｒｏｎ　Ｒａｄｉａｔｉｏｎ　Ｆａｃｉｌｉｔｙ，１９９６．　

［６］黄本诚，马有礼．航天器空间环境试验技术　

［Ｍ］．北京：国防工业出版社，２００１：２４—２８．　

［７３　ＣＨＥＮ　Ｃ　Ｈ，ＴＡＬＮＡＧＩ　Ｊ，ＬＩＵ　Ｊ，ｅｔ　ａ１．Ｔｈｅ　

ｅｆｆｅｃｔ　ｏｆ　ｎｅｕｔｒｏｎ　ｉｒｒａｄｉａｔｉｏｎ　ｏｎ　ＮｄＦｅＢ　ａｎｄ　

Ｓｍ２　Ｃｏ１７　ｂａｓｅｄ　ｈｉｇｈ—ｔｅｍｐｅｒａｔｕｒｅ　ｍａｇｎｅｔｓ［Ｊ］．　

ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｍａｇ，２００５，４１（１０）：３　８３２—３　８３４．　

［８］ＡＬＤＥＲＭＡＮ　Ｊ，Ｊ０Ｂ　Ｐ　Ｋ，ＰＵＨＬ　Ｊ．Ｍｅａｓｕｒｅ—　

ｍｅｎｔ　ｏｆ　ｒａｄｉａｔｉｏｎ－ｉｎｄｕｃｅｄ　ｄｅｍａｇｎｅｔｉｚａｔｉｏｎ　ｏｆ　

ＮｄＦｅＢ　ｐｅｒｍａｎｅｎｔ　ｍａｇｎｅｔｓ［Ｊ］．Ｎｕｃｌ　Ｉｎｓｔｒｕｍ　

Ｍｅｔｈｏｄｓ，２００２，Ａ４８１：９－２８．　

［９］ＯＫＵＤＡ　Ｓ，ＯＨＡＳＨＩ　Ｋ，ＫＯＢＡＹＡＳＨＩ　Ｎ．　

Ｎｕｃｌ　Ｉｎｓｔｒｕｍ　Ｍｅｔｈｏｄｓ，１９９４，Ｂ９４：２２７—２３２．　

Ｅｌ０］ＣＯＳＴ　Ｊ　Ｒ，ＢＲＯＷＮ　Ｒ　Ｄ，ＧＩＯＲＧＩ　Ａ　Ｌ，ｅｔ　ａ１．　

Ｅｆｆｅｃｔｓ　ｏｆ　ｎｅｕｔｒｏｎ　ｉｒｒａｄｉａｔｉｏｎ　ｏｎ　Ｎｄ－Ｆｅ—Ｂ　ｍａｇ—　

ｎｅｔｉｃ　ｐｒｏｐｅｒｔｉｅｓ［Ｊ］．ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｍａｇ，１　９８８，　

２４（３）：２　０１６－２　０１９．　

［１１］ＢＬＡＣＫＭＯＲＥ　Ｅ　Ｗ．Ｒａｄｉａｔｉｏｎ　ｅｆｆｅｃｔｓ　ｏｆ　ｐｒｏｔｏｎｓ　

ｏｎ　ｓａｍａｒｉｕｍ－ｃｏｂａｌｔ　ｐｅｒｍａｎｅｎｔ　ｍａｇｎｅｔｓ［Ｊ］．　

ＩＥＥＥ　Ｔｒａｎｓ　ｏｎ　Ｍａｇ，ｌ９８５，３２（５）：３　６６９—３　６７１．　

［１２］阿卡别尔金，张礼．低能光核过程［Ｍ］．北京：　

人民教育出版社，１９６０．　

［１３］ＦＥＲＲＡＲＩ　Ａ，ＳＡＬＡ　Ｐ　Ｒ，ＦＡＳＳＯＡ，ｅｔ　ａ１．　

ＦＬＵＫＡ：Ａ　ｍｕｌｔｉ－ｐａｒｔｉｃｌｅ　ｔｒａｎｓｐｏｒｔ　ｃｏｄｅ［ＣＰ］．　

Ｓｗｉｔｚｅｒｌａｎｄ：ＣＥＲＮ，２００５．　

［１４］ＡＳＴＭ．ＡＳＴＭ　Ｅ７２２—９４　Ｓｔａｎｄａｒｄ　ｐｒａｃｔｉｃｅ　ｆｏｒ　

ｃｈａｒａｃｔｅｒｉｚｉｎｇ　ｎｅｕｔｒｏｎ　ｅｎｅｒｇｙ　ｆｌｕｅｎｃｅ　ｓｐｅｃｔｒａ　ｉｎ　

ｔｅｒｍｓ　ｏｆ　ａｎ　ｅｑｕｉｖａｌｅｎｔ　ｍｏｎｏｅｎｅｒｇｅｔｉｃ　ｎｅｕｔｒｏｎ　ｆｌｕ—　

ｅｎｃｅ　ｆｏｒ　ｒａｄｉａｔｉｏｎ：Ｈａｒｄｎｅｓｓ　ｔｅｓｔｉｎｇ　ｏｆ　ｅｌｅｃｔｒｏｎ—　

ｉｃｓ［Ｓ］．ＵＳ：Ａｍｅｒｉｃａｎ　Ｓｏｃｉｅｔｙ　ｆｏｒ　Ｔｅｓｔｉｎｇ　ａｎｄ　

Ｍａｔｅｒｉａｌｓ，２００２．　

［１５－１　ＧＲＩＦＦＩＮ　Ｐ　Ｊ，ＫＥＬＬＹ　Ｊ　Ｇ，ＬＵＥＲＡ　Ｔ　Ｆ，ｅｔ　ａ１．　

ＳＮＬＲＭＬ　ｒｅｃｏｍｍｅｎｄｅｄ　ｄｏｓｉｍｅｔｒｙ　ｃｒｏｓｓ　ｓｅｃｔｉｏｎ　

ｃｏｍｐｅｎｄｉｕｍ，ＳＡＮＤ９２—００９４　ＵＣ－７１３［Ｒ］．ＵＳ：　

Ｓａｎｄｉａ　Ｎａｔｉｏｎａｌ　Ｌａｂｏｒａｔｏｒｉｅｓ，１９９３．　

［１６］ＨＵＨＴＩＮＥＮ　Ｍ，ＡＡＲＮ１０　Ｐ　Ａ．Ｐｉｏｎ　ｉｎｄｕｃｅｄ　

ｄｉｓｐｌａｃｅｍｅｎｔ　ｄａｍａｇｅ　ｉｎ　ｓｉｌｉｃｏｎ　ｄｅｖｉｃｅｓ［Ｊ￣．Ｎｕｃｌ　

Ｉｎｓｔｒｕｍ　Ｍｅｔｈｏｄｓ，１９９３，Ａ３３５：５８０—５８２．