高效精确的辐射屏蔽优化设计方法研究报告

一、 研究作者、机构及发表信息

本研究由哈尔滨工业大学的韩宇、应涛、朱赫、杨剑群（IEEE会员）和李兴冀（IEEE高级会员）共同完成。论文标题为“Efficient and Accurate Optimal Design Method for Radiation Shielding”，发表于 IEEE Transactions on Nuclear Science 期刊，卷71，第11期，出版时间为2024年11月。

二、 学术背景与研究目标

本研究属于航天器辐射防护工程与计算物理交叉领域，具体聚焦于航天器电子元器件在轨任务期间所遭受的总电离剂量（Total Ionizing Dose， TID）的精确计算与屏蔽结构的优化设计。

研究背景与动因：随着商业航天的快速发展，现代航天工业对航天器设计提出了低成本、高效率、短周期的迫切要求。空间辐射环境（如地球、火星或行星际空间中的带电粒子）会对航天器上的电气、电子、机电和电光（EEEE）设备造成干扰、损伤甚至毁坏，其中TID作为一种累积效应，直接影响航天器的质量预算（基于屏蔽结构设计）和任务寿命（基于TID暴露限值）的估算。因此，为敏感元器件设计有效的辐射屏蔽至关重要。传统的主流三维辐射传输仿真工具主要采用两种方法：光线追踪（Ray-Tracing， RT）方法和蒙特卡洛（Monte-Carlo， MC）方法（包括正向蒙特卡洛（FMC）和逆向蒙特卡洛（Reverse Monte-Carlo， RMC）方法）。RT方法计算速度极快（秒级），但其基于“直线穿透”假设，对于电子等散射明显的粒子，计算结果与更精确的MC方法相比可能存在较大误差（文献指出可达-70%至300%）。而MC方法（尤其是FMC）虽然精度高，但计算资源消耗巨大，对于复杂几何结构耗时过长，难以满足工程快速迭代需求。RMC方法在精度上与FMC接近且速度更快，但若将其反复用于屏蔽设计的迭代优化，其计算成本和时间依然可观，难以满足现代航天工业对设计效率的要求。

研究目标：针对上述挑战，本研究旨在开发一种新的辐射屏蔽优化设计方法，以有机融合RT方法和RMC方法的优势，弥补各自单独使用时的不足。核心目标是：在确保敏感元器件TID仿真结果准确性的前提下，显著提升辐射屏蔽增强设计的效率，并最终实现屏蔽结构质量的优化，即用更少的附加屏蔽质量达到相同的防护效果，从而节省宝贵的航天器有效载荷资源。

三、 详细研究流程与方法

本研究提出了一种创新的辐射屏蔽优化设计流程，该方法并非单一实验，而是一个结合了仿真计算与迭代优化的系统性设计框架。其核心是分阶段、协同使用RT和RMC两种仿真方法，具体流程如下：

流程概述：整个优化设计流程包含四个主要步骤，并通过迭代循环实现最终优化。研究使用的核心软件是哈尔滨工业大学开发的极端环境辐射效应技术计算机辅助设计（eretcad）软件包，该软件的RMC和FMC模块基于Geant4开发，RT模块包含固体斜程算法和法向壳层算法。

详细步骤：

步骤一：基准TID计算与需求判定 * 研究对象：一个包含敏感元器件的简化航天器几何模型。模型中，航天器外壳被简化为一个铝制立方体（1000×1000×1000 mm，壁厚0.8 mm），内部有一个带上下孔洞的铝制立方体壳（100×100×100 mm，壁厚10 mm），敏感元器件（硅球体，半径0.1 mm）置于航天器中心。 * 研究方法与操作：首先，使用RMC方法计算敏感元器件在原始航天器几何结构下的沉积剂量，记为 rmc-tid。同时，设定一个任务要求的TID设计阈值，记为 rmc-tid_threshold。这是后续所有优化工作的基准和目标。若初始rmc-tid已满足要求，则无需优化；否则，进入优化设计流程。在本研究的案例中，初始rmc-tid为9.4520×10³ rad(Si)，假设设计阈值设定为后续传统方法优化后的结果（6.5410×10³ rad(Si)），则需要进行屏蔽增强设计。

步骤二：RT方法分析与初步薄弱环节识别 * 研究方法与操作：使用RT方法对同一原始几何结构进行仿真。一方面，计算得到RT方法下的TID值（rt-tid）。另一方面，也是更关键的一步，RT方法被用来分析敏感点周围空间的等效铝屏蔽厚度分布。具体操作是：以敏感点为中心，将周围4π立体角空间均匀划分为大量角扇区（本研究案例中极角120分，方位角240分，共28,800个基本立体角）。从中心点向每个扇区发射射线，计算每条射线穿透几何结构所经历的等效铝厚度。这样，就得到了一个围绕敏感点的、基于方向的厚度分布图。这个分布图直观地揭示了航天器结构在哪些角度方向上对辐射的防护较为薄弱（等效厚度小），这些薄弱方向是后续进行针对性屏蔽补强的关键依据。

步骤三：基于RT分析的屏蔽增强初步设计 * 研究方法与操作：此步骤是优化设计的核心创新环节。目标是设计出附加的屏蔽补丁块，以最小的质量代价将TID降低到目标值以下。 1. 设定RT阈值：将RMC的设计阈值 rmc-tid_threshold 暂时作为RT方法的优化目标阈值，记为 rt-tid_threshold。 2. 剂量分析与补强厚度计算：基于步骤二得到的每个角扇区的等效厚度和相应的剂量深度曲线（Dose Depth Curve， DDC），分析各方向的剂量贡献。由于DDC是单调递减函数，等效厚度越小的方向贡献的剂量越高。算法从等效厚度最小的射线开始，逐步选择更厚的等效厚度值，将所有厚度低于此值的射线“分配”到该厚度（即认为在这些薄弱方向上增加了相应厚度的屏蔽材料），并重新计算整体TID。通过二分查找法，迭代调整这个“分配厚度”，直到计算出的整体RT-TID接近 rt-tid_threshold。此时，每条射线都被赋予了一个所需的“补丁厚度”值。 3. 空间映射与补丁生成：为了将基于射线角度的补强方案转化为实际可制造的几何结构，需要在敏感元器件周围定义一个待补强的薄壳表面（本研究案例中是一个50×50×50 mm的薄铝壳A）。将该表面划分为网格（如6×6或10×10）。对于每个网格面元，其所需的补丁厚度计算为：所有与该面元相交的射线的补丁厚度值，按其空间立体角比例进行加权平均。由此，生成一个非均匀的、与辐射薄弱方向精准对应的三维屏蔽补丁结构。计算此迭代后几何结构下的RT-TID（iterated-rt-tid）。

步骤四：RMC精度验证与迭代优化 * 研究方法与操作：使用RMC方法对步骤三设计出的带补丁结构的航天器几何模型进行高精度TID计算，得到 iterated-rmc-tid。 * 迭代判断：比较 iterated-rmc-tid 与最初的 rmc-tid_threshold。 * 如果 iterated-rmc-tid > rmc-tid_threshold，说明基于RT初步设计的屏蔽仍不足。此时，需要降低 rt-tid_threshold（即设定一个更严格的RT优化目标），然后返回步骤三重新进行补强设计。 * 如果 iterated-rmc-tid 足够接近且不高于 rmc-tid_threshold，则优化完成。 这个迭代过程通常只需很少的几次（1-2次），因为RT分析已经为补强提供了精准的方向性指导，避免了RMC方法的盲目试错，从而极大提高了整体优化效率。

四、 主要研究结果

本研究通过上述流程，对一个简化的航天器模型进行了优化设计，并与传统方法（在敏感元件外围立方体六个面上简单附加等厚平板）进行了对比。

1. 传统方法结果：在薄壳A的六个外表面各附加3 mm厚的铝板后，RMC计算得到的TID为6.5410×10³ rad(Si)，附加屏蔽体积为一定值（作为比较基准）。
2. 优化方法结果（6×6网格）：采用本文提出的优化方法，以传统方法的结果（6.5410×10³ rad(Si)）作为 rmc-tid_threshold。经过优化设计后，生成的补丁结构不再是均匀平板，而是在辐射薄弱区域（特别是对应立方体壳上下孔洞的方向）集中增厚，其他区域较薄甚至无需增厚。优化后的结构经RMC验证，其TID为6.5406×10³ rad(Si)，完全满足防护要求。关键成果：在达到相同防护水平的前提下，优化方法所需的附加屏蔽体积比传统方法减少了 20.37%，质量相应同比减少。
3. 优化方法结果（10×10网格）：为进一步优化，将网格加密至10×10。优化设计出的补丁结构更加精细。结果显示，其RMC-TID进一步降至6.5403×10³ rad(Si)，同时附加屏蔽体积相比传统方法减少了 24.75%。这证明更精细的网格划分能实现更优的轻量化效果。
4. 非对称结构验证：为证明方法在复杂场景下的适用性，研究还测试了一个非对称几何结构（内部立方体壳在x轴方向两侧厚度不同）。优化结果显示，即使在非对称条件下，该方法依然能有效工作，相比传统均匀附加平板的方法，实现了 18.59% 的屏蔽质量减少。
5. 角分布可视化分析：研究提供了屏蔽厚度角分布图进行对比。结果显示，传统方法只是均匀增加了整体厚度，并未改变原始几何的厚度分布模式；而优化方法则精准地填补了沿y轴方向（特别是孔洞区域）的薄弱区域，使角向厚度分布更加均衡；网格加密后，厚度分布则更为均匀。这直观地证明了优化方法的“精准补强”特性。

结果逻辑与结论贡献：上述结果链条清晰地展示了研究流程的有效性。步骤一和步骤二建立了基准并识别了问题（薄弱方向）。步骤三利用RT的快速分析能力，生成了一个理论上（基于RT模型）能达标且质量更优的初步设计。步骤四通过RMC这一高精度“裁判”进行验证和微调。最终数据证明，通过RT与RMC的协同，确实能够在保证最终TID计算结果（由RMC确认）准确性的前提下，设计出比传统经验方法质量显著降低的屏蔽结构。这些具体的质量减少百分比数据，直接支撑了论文关于该方法能“显著减少辐射屏蔽质量、节约有效载荷资源”的核心结论。

五、 研究结论与价值

本研究成功提出并验证了一种高效、精确的辐射屏蔽优化设计方法。其核心结论是：通过有机融合RT方法（高效、能提供方向性薄弱信息）和RMC方法（高精度、可作为验证基准）的优势，可以克服各自单独使用的缺点，实现航天器辐射屏蔽设计的快速优化与精确保证。

科学价值与应用价值： * 方法学价值：为解决工程中“精度”与“效率”的矛盾提供了一个创新性的框架。它改变了传统要么牺牲精度求速度（只用RT）、要么牺牲速度求精度（只用MC迭代）的困境，通过分工协作（RT指导方向，RMC把关精度），实现了两者兼顾。 * 工程应用价值：直接服务于现代航天工业，特别是商业航天对低成本、轻量化、短周期的迫切需求。该方法能帮助辐射工程师在更短的设计周期内，设计出质量更优的屏蔽方案，直接降低航天器发射质量，节约发射成本，提升任务效益。 * 软件工具价值：该方法在eretcad软件包中实现，为航天器设计单位提供了一个可操作、一体化的专业工具，有助于提升我国在航天器抗辐射设计领域的自动化与智能化水平。

六、 研究亮点

1. 创新性的方法融合：本研究最大的亮点在于提出了一个系统性的流程，将RT与RMC两种主流仿真方法不是作为替代选项，而是作为设计流程中前后衔接、功能互补的环节，充分发挥各自长处。
2. 显著的工程效益：研究通过具体案例量化证明了该方法的优势：在确保防护性能的前提下，可实现高达24.75% 的附加屏蔽质量减少。这对于“克克计较”的航天器设计具有重大现实意义。
3. 流程的实用性与高效性：优化流程清晰、可操作。利用RT快速定位薄弱环节，极大减少了需要高精度RMC仿真的迭代次数，从而在整体上大幅缩短了优化设计周期，满足了工程快速迭代的需求。
4. 解决了传统方法的两个关键问题：论文明确指出传统方法的两大问题（RT精度不足、RMC迭代优化效率低下），并针对性地提出了解决方案，体现了明确的问题导向。

七、 其他有价值内容与未来展望

论文在讨论部分坦诚指出了当前方法的局限性并对未来工作进行了展望： * 局限性：目前使用的RMC方法尚未包含质子核相互作用的物理过程，因此TID结果未计入核反应次级粒子的剂量贡献。虽然这部分贡献通常较小，但完善物理过程将进一步提升RMC仿真精度。此外，当前方法主要针对TID进行优化设计，对于位移损伤剂量（DDD）或单粒子效应（SEE）等其他辐射效应，该方法能提供屏蔽补强方向，但尚不能给出定量的优化设计。 * 未来工作：作者计划从两方面拓展：一是开发适用于多个敏感点同时优化的屏蔽设计方法，以应对航天器分系统级的需求；二是将方法应用于更复杂的真实几何结构进行验证和优化，以增强方法的普适性和工程实用性。

本研究是一项紧密结合工程实际需求、方法具有显著创新性、且结果证明有效的优秀工作，为航天器辐射防护设计领域贡献了一种兼具高效性与精确性的新思路和新工具。