简化管束设备室整体临界安全分析方法研究.pdf

简化管束设备室整体临界安全分析方法研究陈 添,张毅诚,胡小利,费钧天,邵 增,易 璇,杨海峰(中国核电工程有限公司,北京1 0 0 8 4 0)摘 要:对处理溶液类易裂变材料设备室的整体临界安全分析,提出了简化管束设备室整体临界安全分析方法,保守地将设备室内盛装易裂变材料的管道的料液集中布置成为简化管束,并紧贴设备室内反应性最大的设备上,以全面保守地考虑管道对临界安全的影响。基于简化管束模型,针对地震工况下的管道位移,采用保守且合理的模型进行地震工况下设备室整体临界安全分析。计算结果表明,简化管束整体临界安全分析方法比按实际管道布置精细化建模的方法的计算分析过程更简便,是一种合理且保守的考虑管道的设备室整体临界安全分析方法。关键词:简化管束;临界安全分析;设备室;管道中图分类号:T L 3 4 9;T L 7 3 1文献标志码:A D O I:1 0.1 2 0 6 1/j.i s s n.2 0 9 5 6 2 2 3.2 0 2 4.0 1 0 4 0 4A n a l y s i s M e t h o d o f C r i t i c a l i t y S a f e t y o f E n t i r e E q u i p m e n t R o o m U s i n g S i m p l i f i e d P i p e B u n d l e M o d e lCHE N T i a n Z HANG Y i c h e n g HU X i a o l i F E I J u n t i a n S HAO Z e n g Y I X u a n YANG H a i f e n g收稿日期:2 0 2 3 0 6 0 1;修回日期:2 0 2 4 0 1 1 6基金项目:国防科工局乏燃料后处理科研专项基金资助项目(KY 2 0 0 0 8)作者简介:陈添(1 9 9 2-),男,江苏连云港人,工程师,硕士,主要从事反应堆物理研究。E-m a i l:c h e n t i a n 3 2 11 2 6.c o m C h i n a N u c l e a r P o w e r E n g i n e e r i n g C o L t d B e i j i n g 1 0 0 8 4 0 C h i n a A b s t r a c t T o i n c l u s i v e l y c o n s i d e r t h e i m p a c t o f p i p e s t r a n s p o r t i n g o n c r i t i c a l i t y s a f e t y t h e c r i t i c a l i t y s a f e t y a n a l y s i s m e t h o d o f t h e e n t i r e e q u i p m e n t r o o m u s i n g s i m p l i f i e d p i p e b u n d l e m o d e l i s p r o p o s e d i n t h i s p a p e r T h i s m e t h o d s i m p l i f i e s t h e p i p e s t h a t t r a n s p o r t d i f f e r e n t s o l u t i o n i n t o a h y p o t h e t i c a l p i p e b u n d l e m o d e l a d j a c e n t t o t h e e q u i p m e n t w i t h t h e l a r g e s t r e a c t i v i t y B a s e d o n s i m p l i f i e d p i p e b u n d l e m o d e l a n d t h e e a r t h q u a k e c o n d i t i o n d i s p l a c e m e n t d a t a t h e c r i t i c a l i t y s a f e t y a n a l y s i s o f t h e e n t i r e e q u i p m e n t r o o m i s c o n d u c t e d T h e c a l c u l a t i o n r e s u l t s s h o w t h a t c r i t i c a l i t y s a f e t y a n a l y s i s m e t h o d o f t h e e n t i r e e q u i p m e n t r o o m u s i n g s i m p l i f i e d p i p e b u n d l e m o d e l i s m o r e c o n s e r v a t i v e t h a n t h a t u s i n g a c t u a l p i p e l i n e l a y o u t m o d e l W i t h a s i m p l e a n a l y s i s p r o c e s s t h i s m e t h o d i s a r e a s o n a b l e c o n s e r v a t i v e p i p e-c o n s i d e r e d c r i t i c a l i t y s a f e t y a n a l y s i s m e t h o d o f t h e e n t i r e e q u i p m e n t r o o m K e y w o r d s s i m p l i f i e d p i p e b u n d l e c r i t i c a l i t y s a f e t y a n a l y s i s e q u i p m e n t r o o m p i p e s 在处理溶液类易裂变材料的设备室中,除了与临界安全相关的主要设备外,还存在着大量传输、转运易裂变材料铀、钚料液的管道,设备室整体临界安全分 析必须考虑 管道 的 影 响。国 内 已1-404010第1 5卷 第1期2 0 2 4年2月现 代 应 用 物 理MO D E R N A P P L I E D P HY S I C SV o l.1 5,N o.1F e b.2 0 2 4有对管道本 身临界安全 分析评价 的 相 关 标 准1和分析研究23,但是还未 有考虑管 道的设 备 室整体临界安 全分析方法 的相关 研 究。设 备 室 中的管道布置复杂、数量众多、粗细不一,转运的料液种类繁多,有的管道只在部分时间传输料液。在进行设备室整体临界安全分析时,常用精细化建模的方法进行计算,根据每根管道的起始位置、直径、壁厚、长度、管壁材料及 料液成分等参数进行精细化逐一建模,并使用3维蒙特卡罗程序计算。精 细化建模计 算的优 点 是 考 虑 了 管道的实际 布 置,缺 点 是 建 模 工 作 十 分 烦 琐。此外,使用蒙特卡罗程序临界计算时,源中子的抽样需要合理地 分布到含 易裂变材料 的每个管 道或设备,无 法 保 证 所 有 的 区 域 或 管 道 都 被 抽 样到。对于松散耦合系统,蒙特卡罗抽样过程容易收敛到局部,抽样的随机 性无法保 证4,还 需 要消耗大量的 计算时 间5。精 细 化 分 析 的 结 果 不具有包络性,当管道布置设计或现场实际施工布置出现变更时,需要在庞大的计算模型中准确定位变更之处并更新,这是一个几乎不可实现的任务。另一种常 用的方法是 将含易 裂 变 材 料 的 料液简单聚拢为简单几何体形状,例如溶液墙。这种保守方法的优点是建模简便,缺点是这种过度保守的方法在管道料液总量较大时,计算结果很难满足临界安全的要求。因此,在进行设备室整体临界安全分析中,考虑管道对临界安全的影响是需要解决的难点问题。为此,本文提出了一种适度保守并且 建模较为 简便的考虑 管道布置 的简化管束设备室整体临界安全分析方法。1简化管束设备室整体临界安全分析流程简化管束设备室整体临界安全分析方法,需集中布置设备室内盛装易裂变材料的管道,并紧邻反应性最大的设备进行临界安全分析。设备室整体临界安全分析,除了考虑主要的工艺设备外,还需要考虑大量管道本身的临界安全以及管道与设备之间的相互作用。在处理溶液类易裂变材料的设备室中,管道的布置相对分散,但管道间要求有最小布置间距。与实际分散的管道相比,集中布置的管道使中子相互作用显著增强;同时,集中布置的管束紧邻反应性最大的设备,使管道和设备间的中子相互作用也显著增强,这二者让临界安全分析具有保守性。因为管道壁厚和管径的数量关系,所以需保守考虑管壁的腐蚀量。设备室中很多低浓度的料液或酸溶液、水溶液通过管道,这些管道主要起屏蔽作用,对设备室整体反应性的影响是负效果。基于上述的一些分析原则,简化管束设备室整体临界安全分析方法的流程主要分为管道筛选处理、建立简化管束模型和基于简化管束模型的设备室整体临界安全分析3部分。本文采用3维蒙特卡罗中子输运计算程序MONK-1 0 A进行临界安全计算分析。1.1管道筛选处理首先根据管道参数对管道进行筛选处理,管道参数包括起始点、长度、半径、介质及旋转角度等。梳理设备室内管道介质的种类(水、酸、有机相及易裂变料液等),将料液为气体或不含易裂变材料溶液的管道筛除。经分析,纯水管道对设备室整体临界分析中所起的屏蔽作用大于慢化反射作用,会降低整体设备室反应性,因此设备室整体临界分析中可以保守地不考虑纯水等不含有易裂变材料管道。管道内铀浓度和或钚浓度小于0.1 gL-1的料液,反应性很小,对设备室整体的反应性贡献与纯水基本相同,可以不考虑这些料液传输的管道。对有可能传输多种介质的管道,选取其中最大反应性的料液状态。仅部分时段传输料液的管道,也按照同时充满进行考虑,并且选取其中最大反应性的料液状态。在有效管道中,有一些直径较大、长度较小的管道接口和小设备,将这些管道接口、小设备单独进行建模,可以不用作为简化管束的最大直径。按照上述条件筛选后的有效管道,包含了反应性较大的铀钚浓度料液的管道,统计得到这些管道料液总体积。在计算体积时,使用有效管道的外径计算每个管道的料液体积,全面保守地考虑管壁被腐蚀的情况。有效管道内料液总体积、最粗管道的直径及其对应的最小布置间距等数据都被用于简化管束的建模。1.2简化管束模型简化管束模型,即将设备室内的所有内部存有料液的有效管道集中布置成为一个假想的管束模型。简化管束管径取有效管道中最大管径,管道的间距为有效管道中最大管径管道的最小布置间距,管道料液中铀和钚的浓度取有效管道料液中反应性最大的铀和钚的浓度组合。简化管束模型中所有管道的长度均相同,管道壁厚需考虑一定的腐蚀裕量。简化管束模型中所有的管道按照三角形阵列紧密排列,组成轮廓为六边形的管束形状,图1为简化管束模型示意图。2-404010 陈 添 等:简化管束设备室整体临界安全分析方法研究第1期图1简化管束模型示意图F i g.1 S c h e m a t i c d i a g r a m o f s i m p l i f i e d p i p e b u n d l e m o d e l为了得到反应性最大的简化管束模型,保持管束模型的料液总体积,即有效管道料液总体积不变,逐步增大六边形管束的圈数,相应地逐步减小管束的长度,通过搜索计算最终得到反应性最大的六边形简化管束模型。对某个料液总体积下的简化管束进行了不同管径下的全面搜索计算,简化管束不同管径下最大反应性尺寸搜索结果如图2所示。其中,l1为管束的长度,l2为六边形管束的对边距。由图2可见,l1/l2在1附近时简化管束的有效增殖因子ke f f最大,另外简化管束的ke f f随着管径的增大而增大,因此选取有效管道中的最大管径进行简化管束的建模是保守的。图2简化管束不同管径下最大反应性尺寸搜索结果F i g.2 M a x i m u m r e a c t i v i t y o f s i m p l i f i e d p i p e b u n d l e m o d e l w i t h d i f f e r e n t d i a m e t e r s o f p i p e s 按照上述流程建立简化管束模型,并假设设备室所有有效管道内的料液聚在一起,按照有效管道的最高料液浓度、最大管径、最小间距来建模,并且考虑了最大反应性的简化管束尺寸,因此简化管束模型是较为保守的,包络设备室中管道之间的相互作用。在设备室整体临界安全分析中,将得到的最大反应性尺寸的简化管束放入设备室中,进行设备室整体临界安全分析。1.3设备室整体临界安全分析以相互作用最大的方式,将搜索计算得到的具有最大反应性的简化管束模型紧贴布置在设备室内反应性最大的设备旁。图3为简化管束设备室整体临界计算模型横截面图。计算设备室内各设备单体的ke f f,选择ke f f最大的设备,将简化管束模型紧邻放置于该设备附近。若该设备附近还存在其他ke f f较大的设备,考虑简化管束模型同时紧邻或朝向其他ke f f较大的设备放置模型;若该设备附近没有其他ke f f较大的设备,考虑墙体反射作用的影响,简化管束模型同时紧邻或朝向混凝土墙体放置模型。若设备室内存在多个ke f f较大值设备,以简化管束模型紧 邻 不 同ke f f较 大 设 备 的 情 况 为 参 考 来 放 置模型。图3简化管束设备室整体临界计算模型横截面图F i g.3 C r o s s-s e c t i o n d i a g r a m o f c a l c u l a t i o n m o d e l o f e q u i p m e n t r o o m w i t h s i m p l i f i e d p i p e b u n d l e 简化管束设备室整体临界安全分析方法还可以进行地震工况下的临界安全分析。当采用精细化建模方法时,由于管道布置错综复杂,无法确定每根管道地震位移的方向,难以得到一种保守情形;保守考虑地震工况下有效管道最大变形量对设备室整体临界安全的影响,采用简化管束设备室整体临界安全分析方法。在建模过程中,将管道的位移量应用在简化管束模型中。图4为地震工况下简化管束设备室整体临界计算模型横截面图。除最靠近设备的一层管道外,考虑所有管道同时朝最具反应性的设备发生最大位移,更加靠近该设备,如图4(a)所示。另外,对地震工况下的一个极限情形补充分析:保守考虑管束环绕紧贴反应性最大的设备,由于紧贴设备,简化管束中最内一圈管道产生形变。如图4(b)3-404010第1 5卷现 代 应 用 物 理所示。分析中,保持总料液体积不变,对简化管束的高径比进行进一步反应性搜索,包括管道在地震工况下的可能变形后与设备相互作用的变化情况。(a)P i p e b u n d l e d e f o r m a t i o n(b)P i p e b u n d l e s u r r o u n d i n g e q u i p m e n t图4地震工况下简化管束设备室整体临界计算模型横截面图F i g.4 C r o s s-s e c t i o n d i a g r a m o f c a l c u l a t i o n m o d e l o f e q u i p m e n t r o o m w i t h s i m p l i f i e d p i p e b u n d l e u n d e r e a r t h q u a k e c o n d i t i o n简化管束设备室整体临界安全分析方法采取了多种保守假设,当计算结果不能满足临界安全限值要求时,可考虑释放分析中的部分保守度。比如,将设备室所有易裂变料液按照简化管束集中布置时,将反应性最大设备所在的一面墙体附近的料液集中布置,其余三面墙体的料液集中布置在其余三面墙体反应性最大的设备附近;再如,由于简化管束模型采用了有效管道中的最高料液浓度和最大管径,若最大管径或最大料液浓度的有效管道数目较少,可以将这些有效管道精细化建模,再以剩余的有效管道进行简化管束建模,这么做可降低简化管束的料液浓度或减小管径,释放一定的保守度。在释放保守度后,仍不能满足临界安全限值的情况下,考虑使用精细化建模的方法分析设备室的整体临界安全。确认不能满足临界安全要求的情况下,应考虑变更设备室的管道布置,使其最终满足临界安全要求。设备室整体临界安全分析的总体思路和分析流程如图5所示。图5简化管束设备室整体临界安全分析的总体思路和分析流程F i g.5 A n a l y s i s p r o c e s s o f t h e e q u i p m e n t r o o m c r i t i c a l i t y s a f e t y a n a l y s i s m e t h o d u s i n g s i m p l i f i e d p i p e b u n d l e m o d e l2简化管束设备室整体临界安全分析方法 的应用以某设备室为例,应用简化管束设备室整体临界安全分析方法。该设备室的管道共包含8 0 0余根管道和2 3个小设备,经过有效管道筛选,共有4 5 0根有效管道,对其中半径较大的管道接口和小设备进行精细化建模处理。统计该设备室所有有效管道含料液管道的总体积,根据最粗管径按照最小布置间距建立六边形的简化管束模型,管道料液为反应性最大的料液,该料液钚浓度较高。逐步增加管道的层数,当 管 束 有8层1 6 9根 管 道 时ke f f最 大,为0.4 7 2 30.0 0 0 5。在设备室中选取ke f f最大的设备,将得到的ke f f最大时的管束(8层1 6 9根D N 4 0管道)紧贴放置于4-404010 陈 添 等:简化管束设备室整体临界安全分析方法研究第1期该设备与临近的另一个设备之间。此外,还对管束在设备室中的其他布置形式进行了分析,得到六边形,方形,矩形管束的ke f f分别为0.7 5 5 00.0 0 0 5,0.7 4 0 60.0 0 0 5,0.7 4 1 30.0 0 0 5。由计算结果可知,方形管束的结果略小于六边形的管束结果。对该设备室进行管道精细化建模计算,对该设备室内的设备和所有管道根据实际布置逐一建模。将简化管道设备室整体临界安全分析的结果与精细化建模计算的结果及不考虑管道的设备室临界计算结果进行对比,不同模型计算得出的ke f f如表1所列。由表1可以看出,简化管束设备室整体临界安全分析方法的结果大于不考虑管道的设备室结果和管道精细化建模结果,是具有保守性的。由表1还可知,管道精细化建模结果与不考虑管道的设备室结果基本相等,这是因为与主要的临界安全相关设备相比,管道对设备室整体的临界安全计算结果的影响较小。表1不同模型计算得出的ke f fT a b.1 C a l c u l a t i o n r e s u l t s o f ke f f b y u s i n g d i f f e r e n t p i p e l i n e l a y o u t m o d e l sC a l c u l a t i o n c o n d i t i o nke f fE q u i p m e n t r o o m w i t h o u t p i p e s0.7 4 0 70.0 0 0 5E n t i r e e q u i p m e n t r o o m w i t h p i p e s u s i n g a c t u a l p i p e l i n e l a y o u t m o d e l0.7 4 0 70.0 0 0 5E n t i r e e q u i p m e n t r o o m w i t h p i p e s u s i n g s i m p l i f i e d p i p e b u n d l e m o d e l0.7 5 5 00.0 0 0 5采用简化 管 束 设 备 室 整 体 临 界 安 全 分 析 方法进行地震工况下的分析。经过力学计算,地震工况下的有效管道最大位移为3 8.1 mm,考虑简化管束,除最靠近设备的一层管道外,所有管道同时朝最具反应性的设备发生最大位移,同时设备室的设备 都考虑地震 变形,ke f f的 计 算 结 果 为0.7 7 8 60.0 0 0 5。作为地震工况下一个极限情形的补充分析,设备室各设备都考虑设备形变,保守考虑管束环绕紧贴反应性最大的设备,改变环绕管束的高 度和半径 找到环绕管 束与设备 最大相互作用的情形,管束环绕设备时管束高度变化的临界计算结果如表2所列。由表2可见:ke f f最大为0.7 9 5 00.0 0 0 5,该结果比简化管束变形的结果更保守。表2管束环绕设备时管束高度变化的临界计算结果T a b.2 C i r t i c a l c a c u l a t i o n r e s u l t s o f e q u i p m e n t s u r r o u n d e d b y p i p e b u n d l e w h e n c h a n g i n g p i p e b u n d l e h e i g h tS i m p l i f i e d p i p e b u n d l e h e i g h t/c mO u t e r r a d i u s o f s i m p l i f i e d p i p e b u n d l e/c mke f f1 2 01 1 2.8 90.7 9 3 30.0 0 0 51 3 01 0 9.1 80.7 9 3 50.0 0 0 51 4 01 0 5.8 90.7 9 4 80.0 0 0 51 5 01 0 2.9 60.7 9 3 60.0 0 0 51 6 01 0 0.3 20.7 9 5 00.0 0 0 51 7 09 7.9 40.7 9 3 20.0 0 0 51 8 09 5.7 70.7 9 2 40.0 0 0 51 9 09 3.7 80.7 9 0 90.0 0 0 53结论本文提出 了 一 种 简 化 管 束 设 备 室 整 体 临 界安全分析方法,可以对设备室中的管道对设备室的影响进行临 界安全分析。该方 法 假 设 设 备 室内的所有盛装易 裂变材料的 管道集中布 置成 简化管束并紧邻反应性最大的设备,以此为基础还可以对地震工 况下的情形 进行 分 析。简 化 管 束设备室整体临界 安全分析方 法采用了很 多保 守假设,根据需要可以释放一定的保守度,在应用该方法时具有 一定的灵活 性。通 过 对 简 化 管 束设备室整体临界 安全分析方 法在某设备 室的 实际应用,计算结果证明该方法比所有管道根据实际布置逐一建模的精细化建模方法更保守,并且可以避免精细化建模方法的建模烦琐、临界安全分析程序的建模 能力不足及 蒙特卡罗抽 样的 局限性等问题,与精细化建模方法相比,还能够对地震工况情形 进行合理保 守的 模 拟。另 外 该 方法也可以避免将 管道料液简 单聚拢导致 无法 满足临界安全要求的问题。因此,简化管束设备室整体临界安全分析方法是一种合理保守、简便且灵活的 考 虑 管 道 的 设 备 室 整 体 临 界 安 全 分 析方法。参考文献 1 反应堆外易裂变材料的核临界安全 第4部分 含易裂变物质水溶液的钢质管道交接的核临界安全准则 G B 1 5 1 4 6 4-1 9 9 4 S 北京 中国标准出版社 1 9 9 5 N u c l e a r c r i t i c a l i t y s a f e t y f o r f i s s i l e m a t e r i a l s o u t s i d e r e a c t o r-N u c l e a r c r i t i c a l i t y s a f e t y c r i t e r i a f o r s t e e l-p i p e i n t e r s e c t i o n s c o n t a i n i n g a q u e o u s s o l u t i o n s o f f i s s i l e m a t e r i a l G B 1 5 1 4 6 4-1 9 9 4 S B e i j i n g 5-404010第1 5卷现 代 应 用 物 理S t a n d a r d s P r e s s o f C h i n a 1 9 9 5 2 于淼 易璇 管道临界安全分析初步研究 C 第十七届反应堆数值计算和粒子输运学术会议暨2 0 1 8年反应堆物理会议 深 圳 2 0 1 8 YU M i a o Y I X u a n P r e l i m i n a r y s t u d y o n c r i t i c a l i t y s a f e t y a n a l y s i s o f p i p e s C C O R P HY 2 0 1 8 S h e n z h e n 2 0 1 8 3 B A R D E L AY A C A P L I N G S A C KO M N u c l e a r c r i t i c a l i t y s a f e t y o f p i p e s n e t w o r k s c o n t a i n i n g u r a n i u m C I n t e r n a t i o n a l C o n f e r e n c e o n N u c l e a r C r i t i c a l i t y S a f e t y C h a r l o t t e U S A 2 0 1 6 4 吴晓飞 吴海成 刘振华 等 蒙特卡罗方法计算特征值问题收敛性研究 J 中国原子能科学研究院年报 2 0 0 6 1 2 1 4 WU X i a o-f e i WU H a i-c h e n g L I U Z h e n-h u a e t a l S t u d y o n c o n v e r g e n c e o f e i g e n v a l u e p r o b l e m u s i n g M o n t e C a r l o m e t h o d J A n n u a l R e p o r t o f C h i n a I n s t i t u t e o f A t o m i c E n e r g y 2 0 0 6 1 2 1 4 5 潘流俊 蒙特卡罗临 界 计 算 中 的 裂 变 矩 阵 加 速 方 法 和p-CMF D方法研 究 D 绵 阳 中 国 工 程 物 理 研 究 院 2 0 1 6 P AN L i u-j u n S t u d y o n t h e f i s s i o n m a t r i x a c c e l e r a t i o n m e t h o d a n d t h e CMF D a c c e l e r a t i o n m e t h o d i n c r i t i c a l i t y c a l c u l a t i o n u s i n g M o n t e C a r l o m e t h o d D M i a n y a n g C h i n a A c a d e m y o f E n g i n e e r i n g P h y s i c s 2 0 1 6 上接第0 1 0 4 0 3-5页 8 张玉芳 裴香涛 尉伟 等 合肥光源贮存环四六极磁铁组合支撑设计及样机振动测试 J 中国科学技术大学学报 2 0 1 3 4 3 6 4 8 6 4 9 0 Z HAN G Y u-f a n g P E I X i a n g-t a o WE I W e i e t a l D e s i g n o f q u a d r u p o l e a n d s e x t u p o l e m a g n e t s s u p p o r t i n g m e c h a n i s m f o r HL S s t o r a g e r i n g a n d v i b r a t i o n t e s t s o f p r o t o t y p e J J o u r n a l o f U n i v e r s i t y o f S c i e n c e a n d T e c h n o l o g y o f C h i n a 2 0 1 3 4 3 6 4 8 6 4 9 0 9 王梓豪 李春华 王之琢 等 基于模态匹配法的接触刚度 J 强激光与粒子束 2 0 1 9 3 1 8 0 8 5 1 0 1 WAN G Z i-h a o L I C h u n-h u a WANG Z h i-z h u o e t a l C o n t a c t s t i f f n e s s m e a s u r e m e n t b a s e d o n m o d a l m a t c h i n g m e t h o d J H i g h P o w e r L a s e r a n d P a r t i c l e B e a m s 2 0 1 9 3 1 8 0 8 5 1 0 1 1 0 陈超祥 胡其登 S o l i d W o r k s S i m u l a t i o n 高级教程 M 北京 机械工业出版 社 2 0 1 8 C HE N C h a o-x i a n g HU Q i-d e n g S o l i d W o r k s S i m u l a t i o n M B e i j i n g M e c h a n i c a l E n g i n e e r i n g P r e s s 2 0 1 8 1 1 黄晴晴 杨振 罗涛 等 四代光源贮存环磁铁支撑设计及仿真优 化 J 强 激 光 与 粒 子 束 2 0 2 4 3 6 1 0 1 4 0 0 1 HUAN G Q i n g-q i n g YANG Z h e n L UO T a o e t a l D e s i g n s i m u l a t i o n a n d o p t i m i z a t i o n o f m a g n e t s u p p o r t s i n 4 t h g e n e r a t i o n s y n c h r o t r o n l i g h t s o u r c e s J H i g h P o w e r L a s e r a n d P a r t i c l e B e a m s 2 0 2 4 3 6 1 0 1 4 0 0 1 6-404010 陈 添 等:简化管束设备室整体临界安全分析方法研究第1期