颗粒团聚行为对弥散型核燃料芯体失效的影响分析.pdf

第 卷第期原子能科学技术V o l ,N o 年月A t o m i cE n e r g yS c i e n c ea n dT e c h n o l o g yA p r 颗粒团聚行为对弥散型核燃料芯体失效的影响分析董颖璇,吕俊男,李群,(西安交通大学 航天航空学院 机械结构强度与振动国家重点实验室,陕西 西安 ;中国核动力研究设计院 反应堆燃料及材料重点实验室,四川 成都 )摘要:针对燃料颗粒团聚行为对弥散型核燃料芯体失效的影响,开发了弥散型核燃料元件代表性体元模型的参数化建模及数值计算脚本,综合考虑燃料基体界面损伤层、燃料颗粒辐照肿胀以及环境压力等关键特征进行了数值建模,系统分析了燃料颗粒团聚体积分数、颗粒团聚位置、运行环境静水压力以及基体材料性质对弥散型核燃料芯体失效的影响规律,发现燃料颗粒团聚区域是弥散型核燃料芯体失效的起点,而运行环境静水压力对芯体应力集中有抑制作用,揭示了燃料颗粒团聚现象对弥散型核燃料元件力学性能的“短板效应”.团聚区域的燃料颗粒通过密集堆叠的方法进行几何建模,局部燃料颗粒团聚体积分数由最小颗粒间距定义.有限元计算结果表明,燃料颗粒团聚体积分数越大,弥散型核燃料芯体越容易失效.而包壳外的环境压力会降低芯体内燃料颗粒团聚区域的应力水平,但随着局部团聚程度的增加,环境压力对芯体中M i s e s应力的影响逐渐减小.此外,燃料颗粒团聚区域相对弥散型核燃料芯体厚度方向位置改变对芯体的最大M i s e s应力几乎没有影响;燃料颗粒团聚区域在弥散型核燃料芯体的面内分布位置和环境压力共同决定弥散型核燃料芯体的失效行为.本研究可为弥散型核燃料元件的失效条件预测、可靠性评估和结构优化设计提供分析方法和数值依据.关键词:弥散型核燃料元件;颗粒团聚;强度分析;短板效应;基体失效中图分类号:T L 文献标志码:A文章编号:()收稿日期:;修回日期:基金项目:国家自然科学基金(,);中核集团青年科技创新团队项目(J T )通信作者:李群;吕俊男d o i:/y z k y o u x i a n A n a l y s i so fP a r t i c l eA g g l o m e r a t i o nE f f e c t o nF a i l u r eo fD i s p e r s i o nN u c l e a rF u e lM e a tD ONGY i n g x u a n,L YUJ u n n a n,L IQ u n,(S t a t eK e yL a b o r a t o r yf o rS t r e n g t ha n dV i b r a t i o no fM e c h a n i c a lS t r u c t u r e s,S c h o o l o fA e r o s p a c e,X ia nJ i a o t o n gU n i v e r s i t y,X ia n ,C h i n a;S c i e n c e a n dT e c h n o l o g yo nR e a c t o rF u e la n dM a t e r i a l sL a b o r a t o r y,N u c l e a rP o w e rI n s t i t u t eo fC h i n a,C h e n g d u ,C h i n a)A b s t r a c t:B yd i s p e r s i n g t h e f u e l p a r t i c l e s i na n i n e r tm e t a lm a t r i x,t h ed i s p e r s i o nn u c l e a r f u e le l e m e n th a st h ea d v a n t a g e so fh i g hb u r n u p,h i g ht h e r m a lc o n d u c t i v i t y,g r e a ts a f e t ya n ds t r o n gd e s i g n a b i l i t y,e t c I t h a sb e e ns u c c e s s f u l l ya p p l i e d i n t h e t e s t r e a c t o r s W i t h i nad i s p e r s i o nn u c l e a rf u e le l e m e n t,t h ef a i l u r eo f f u e lm e a t i ss i g n i f i c a n t l yi n f l u e n c e db yt h ep a r t i c l ea g g l o m e r a t i o nb e h a v i o r,i n c l u d i n gt h ea g g l o m e r a t i o ne x t e n t,t h el o c a t i o no f a g g l o m e r a t i o nr e g i o n,f u e lp a r t i c l es p a c i n g,e t c T h ep r e s e n ts t u d yi n v e s t i g a t e s t h ec r u c i a lf a c t o r si n f l u e n c i n gt h ef a i l u r eb e h a v i o ri nt h ed i s p e r s i o nn u c l e a rf u e lm e a t,f o c u s i n go n t h ep a r t i c l e a g g l o m e r a t i o ne f f e c t o n t h e f a i l u r eo f f u e lm e a t B yd e v e l o p i n gas c r i p tf o rp a r a m e t r i cm o d e l i n ga n dn u m e r i c a lc a l c u l a t i o n so ft h er e p r e s e n t a t i v ev o l u m ee l e m e n t i nt h ed i s p e r s i o nn u c l e a r f u e l e l e m e n t,t h eM i s e ss t r e s s f i e l do fd i s p e r s i o nn u c l e a r f u e lm e a t i n c o r p o r a t i n gw i t hd i f f e r e n tp a r t i c l ea g g l o m e r a t i o nc h a r a c t e r i s t i c sw a ss i m u l a t e di nt h i sp a p e r N u m e r i c a lm o d e l i n g w a sp e r f o r m e db yc o n s i d e r i n gt h ed a m a g e l a y e ra t t h e f u e l m a t r i x i n t e r f a c e,i r r a d i a t i o n i n d u c e df u e l p a r t i c l e ss w e l l i n ga n dt h ee n v i r o n m e n t a lp r e s s u r eo u t s i d et h ec l a d d i n g T h ee f f e c t so fp a r t i c l ea g g l o m e r a t i o nv o l u m e f r a c t i o n,l o c a t i o no ft h ea g g l o m e r a t ef u e lp a r t i c l e s,a m b i e n th y d r o s t a t i cp r e s s u r ea n dm a t r i xm a t e r i a lp r o p e r t i e so nt h ef a i l u r eo fd i s p e r s i o nn u c l e a rf u e lm e a tw e r ea n a l y z e ds y s t e m a t i c a l l y T h ef u e lp a r t i c l e si na g g l o m e r a t i o nr e g i o nw e r eg e o m e t r i c a l l ym o d e l e dt h r o u g hd e n s es t a c k i n g,w h i l e t h e l o c a lp a r t i c l ea g g l o m e r a t i o nv o l u m ef r a c t i o nw a sd e f i n e db y t h em i n i m u mp a r t i c l e s p a c i n g T h e r e s u l t s s h o wt h a t t h e f a i l u r eb e h a v i o ri nd i s p e r s i o nn u c l e a r f u e lm e a t i se n h a n c e dw i t ht h e i n c r e a s i n go fp a r t i c l ea g g l o m e r a t i o nv o l u m e f r a c t i o n T h ee n v i r o n m e n t a lp r e s s u r eo u t s i d e t h ec l a d d i n gc a n i n h i b i t t h es t r e s sc o n c e n t r a t i o n i nt h ed i s p e r s i o nf u e lm e a t,b u t t h i s i n f l u e n c ed e c l i n e sw i t ht h ee x t e n to fl o c a l f u e la g g l o m e r a t i o ni n c r e a s i n g F u r t h e r m o r e,t h ep a r t i c l ea g g l o m e r a t i o nl o c a t i o nb e t w e e nt h eu p p e r a n d l o w e r c l a d d i n g l a y e r s o f t h e t h i np l a t e i n d i c a t e s l i t t l e e f f e c t o n t h em a x i m u m M i s e ss t r e s s i nt h ed i s p e r s i o nn u c l e a rf u e lm e a t H o w e v e r,c o m p a r e dt ot h ed i s p e r s i o nm e a ta d j a c e n tt ot h es u r r o u n d i n gc l a d d i n g,f u e lp a r t i c l e sa g g l o m e r a t i n gi nc e n t r a lm e a t(i e o n l yc o n s t r a i n e db yt h eu p p e ra n d l o w e r c l a d d i n g)a r em o r es e n s i t i v et oe n v i r o n m e n t a lp r e s s u r e T h el o c a lf u e lp a r t i c l ea g g l o m e r a t i o nh a sas h o r t e s ts t a v ee f f e c to nt h e f a i l u r eo fd i s p e r s i o nn u c l e a r f u e lm e a t T h ev o l u m e f r a c t i o no fa g g l o m e r a t i o nw i l l s e r i o u s l ye x a c e r b a t e t h e s t r e s s c o n c e n t r a t i o n,a n d t h e a g g l o m e r a t i o n l o c a t i o n i nt h e i n t e r i o ro fd i s p e r s i o nf u e lm e a tw i l lf a i lf i r s t T h e r e f o r e,i ti sc o n c l u d e dt h a tt h ee x t e n to ff u e lp a r t i c l ea g g l o m e r a t i o n,t h el o c a t i o no fa g g l o m e r a t i o nr e g i o nw i t h i nt h ed i s p e r s i o nn u c l e a r f u e lm e a t,a n dt h ee n v i r o n m e n t a l p r e s s u r ea r ec r u c i a l f a c t o r s i n f l u e n c i n gt h e f a i l u r eb e h a v i o r o f f u e lm e a t T h e f u e l p a r t i c l e a g g l o m e r a t i o n r e g i o n i s t h eo n s e tl o c a t i o nf o rt h el o c a lm a c r o s c o p i cf a i l u r e,a n dt h ee n v i r o n m e n t a lp r e s s u r eo u t s i d et h ec l a d d i n gc a nd e c r e a s et h ep l a s t i cs t r e s sc o n c e n t r a t i o ni nt h ed i s p e r s i o nf u e lm e a tt oi m p r o v et h es a f e t yo fr e a c t o rc o r e T h i s w o r kp r o v i d e sa na n a l y t i c a la p p r o a c ha n dn u m e r i c a l r e f e r e n c ef o rt h ef a i l u r ec o n d i t i o np r e d i c t i o n,t h er e l i a b i l i t ye v a l u a t i o n,a n dt h es t r u c t u r a l o p t i m i z a t i o nd e s i g no f t h ep l a t ed i s p e r s i o nn u c l e a r f u e l e l e m e n t K e yw o r d s:d i s p e r s i o nn u c l e a r f u e l e l e m e n t;p a r t i c l ea g g l o m e r a t i o n;s t r e n g t ha n a l y s i s;s h o r t e s t s t a v ee f f e c t;m a t r i xf a i l u r e弥散型核燃料元件是将易裂变燃料相以小颗粒形式弥散分布在相对惰性的非裂变基体材料中的一种混合型燃料元件,具有导热性能好、抑制燃料肿胀、包容裂变碎片、阻止裂变物质泄漏等优点,可提高反应堆堆芯的燃耗和安全可靠性 .然而,现行混料工艺条件无法保证燃料颗粒在芯体中的均匀分布,这必然导致核燃料颗粒在芯体中的局部团聚现象.燃料颗第期董颖璇等:颗粒团聚行为对弥散型核燃料芯体失效的影响分析粒团聚将造成弥散型核燃料芯体中局部传热性能变差,颗粒之间的肿胀干涉效应变强,致使芯体更容易在颗粒团聚区域出现裂纹萌生和断裂等失效现象,威胁反应堆堆芯的安全运行.因此,核燃料颗粒在基体中的局部团聚行为(即团聚程度、分布位置等)是影响弥散型核燃料元件性能的关键因素,分析颗粒团聚行为对弥散型核燃料芯体失效的影响可为优化燃料元件设计奠定基础.颗粒团聚是颗粒增强复合结构中非均匀分布颗粒的局部聚集现象,故燃料颗粒团聚代表弥散型核燃料元件中燃料相的局部非均匀性增加.弥散型核燃料芯体在结构上属于复合结构.颗粒团聚对复合结构的等效力学性能、失效强度、断裂韧性以及损伤行为等 都有影响,复合结构的断口处裂纹源通常还伴有增强相颗粒的团聚,这种现象可统称为颗粒团聚对复合结构等效力学性能的“短板效应”.颗粒分布均匀性好的复合结构强度高,并且颗粒团聚会降低材料的韧性.弥散型核燃料元件的强度和失效同样取决于芯体中的颗粒团聚区域.同时,燃料元件中的应力场受反应堆一回路中冷却剂产生的环境静水压力作用,而环境压力对燃料颗粒团聚“短板效应”的影响尚无定论.为精确评估弥散型核燃料芯体的失效行为,需要分析燃料颗粒团聚行为对芯体应力场的影响,并讨论环境压力对燃料颗粒团聚“短板效应”的作用.本研究通过构建弥散型核燃料元件的代表性体元(R V E)模型,分析陶瓷燃料颗粒辐照肿胀作用下燃料颗粒团聚对芯体失效的影响.综合考虑数值模型的计算性和收敛性,对三维团聚模型进行合理的等效简化,将非团聚区域进行等效处理,聚焦于燃料颗粒团聚特征(包括团聚体积分数和团聚位置)的影响.其中,非团聚区域的等效处理考虑燃料颗粒和基体材料的热膨胀系数、体积模量和剪切模量.辐照效应考虑基体所受燃料颗粒的辐照肿胀、燃料基体界面的损伤层(I L)以及随辐照剂量和温度变化的材料性质等因素.弥散型核燃料芯体中的应力场及失效载荷条件主要来源于燃料颗粒的辐照肿胀变形及包壳外的环境压力作用.本研究针对燃料颗粒团聚行为对弥散型核燃料元件芯体失效的影响展开分析,通过对弥散型核燃料元件团聚R V E模型的数值模拟,揭示燃料颗粒团聚体积分数、颗粒团聚位置、运行环境静水压力等关键因素对芯体失效的影响机理.数值方法 有限元模型弥散型核燃料元件的结构如图 a所示,由金属包壳、金属基体和非均匀分布的核燃料颗粒组成.局部燃料颗粒分布间距过近导致的颗粒团聚现象会改变芯体中的应力场分布,进而影响弥散型核燃料元件的失效行为.因此,本研究从弥散型核燃料元件中抽象出R V E模型,通过力学等效过程建立含颗粒团聚行为的R V E模型,等效过程如图 b d所示.a 弥散型核燃料元件;b 局部弥散型核燃料元件代表性体元模型;c 非团聚区域等效简化的弥散型核燃料元件代表性体元模型;d 网格划分结果图弥散型核燃料元件的代表性体元模型等效过程F i g E q u i v a l e n tp r o c e s s f o r r e p r e s e n t a t i v ev o l u m ee l e m e n tm o d e l o fd i s p e r s i o nn u c l e a r f u e l e l e m e n t原子能科学技术第 卷从图 a所示的弥散型核燃料元件中抽象出图 b所示的R V E模型,模型包含上下包壳及弥散芯体两部分.R V E模型在整个燃料元件中所处位置不同,其边界条件也不同.为了聚焦燃料颗粒团聚行为对整个芯体强度的影响,假设每个R V E模型中只存在一个团聚区域,芯体其余部分燃料相等效均匀化.根据M o r i T a n a k a等效理论,含非均匀相的复合材料性质可通过等效处理求得其等效材料常数.将图 b中的非团聚区域等效为图 c所 示 的 均 匀 材 料,则 材 料 性 质 可 通 过M o r i T a n a k a等效理论获得 .此外,由于复合结构的破坏应力与颗粒间的距离相关,图 c所示的团聚区域中采用团聚颗粒的最小间距定义团聚程度.基体所受燃料颗粒的辐照肿胀压力作用可通过解析方法获得,取代燃料颗粒实体作用于基体内部.由于燃料颗粒裂变产生的固态裂变产物的扩散作用,金属基体和燃料颗粒间存在约 m的第三相物质损伤层 .本文在颗粒和基体中间构建了 m的损伤层模型.同时,上下包壳表面受反应堆中冷却剂流动导致的环境压力作用.有限元模型的网格划分结果如图 d所示,对几何模型中的损伤层进行网格细化,以提高计算收敛性.a 立方排布;b 密集堆叠图团聚区域的颗粒分布形式F i g P a r t i c l ed i s t r i b u t i o nf o r mi na g g l o m e r a t i o nr e g i o n 燃料颗粒团聚建模方法燃料颗粒团聚建模形式有多种,如立方排布和密集堆叠等形式,如图所示.利用图 a所示的立方排布形式计算得到的颗粒团聚体积分数在 左右,而根据图 b所示密集堆叠排列构建颗粒分布模型的颗粒团聚体积分数相对最高,可达到,故本研究采用密集堆叠排列方式进行弥散型核燃料芯体含颗粒团聚特征的R V E建模.颗粒团聚程度由团聚区域的最小颗粒间距(dt)定义,根据图 b所示的几何构型与团聚区域的局部燃料颗粒体积分数(Vt)的关系,团聚区域的最小颗粒间距可表达为:dtRf(Vt)/()其中,Rf为燃料颗粒半径.非团聚区域等效处理方法非团聚区域采用等效处理的方法进行均一化(图),以聚焦燃料颗粒团聚行为对弥散型核燃料芯体失效的“短板效应”影响进行研究.均一化等效处理时,非团聚区域的材料属性根据M o r i T a n a k a等效理论 进行计算.图非团聚区域中材料等效转化F i g E q u i v a l e n t t r a n s f o r m a t i o no fm a t e r i a lp r o p e r t i e s i nn o n a g g l o m e r a t i o nr e g i o n根据M o r i T a n a k a方法,含异质夹杂的复合结构的等效体积模量(Ke f f)表达式为:Ke f fKmVf(Kp/Km)(Vf)(Kp/Km)()其中:Km和Kp分别为金属基体和燃料颗粒的体积模量,MP a;Vf为非团聚区域的燃料颗粒团聚体积分数;为与夹杂形状和位置相关的参数,根据式()进行计算.KmGmKpGm()等效剪切模量(Ge f f)采用式()表征.Ge f fGmVf(Gp/Gm)(Vf)(Gp/Gm)()其中:Gm和Gp分别为金属基体和燃料颗粒的剪切模量,MP a;为与夹杂形状和位置相关的参数,根据式()进行计算.Gm(KmGm)Gm(KmGm)Gp(KmGm)Gm(KmGm)()而等效热膨胀系数(e f f)采用式()进行计算.e f fp(mp)/Ke f f/Kp/Km/Kp()第期董颖璇等:颗粒团聚行为对弥散型核燃料芯体失效的影响分析其中,m和p分别为金属基体和燃料颗粒的热膨胀系数.燃料颗粒辐照肿胀模拟在弥散型核燃料元件中,芯体的应力场因受燃料颗粒辐照肿胀作用和包壳外的环境水压作用而产生.金属基体在燃料颗粒辐照肿胀压力pf作用下变形失效,pf可通过解析模型 求解获得.单颗粒辐照肿胀压力解析模型考虑了辐照条件、材料性质以及燃料颗粒高燃耗结构的非均匀分布气孔等因素的影响.为提高含颗粒团聚行为的R V E模型计算性能,采用陶瓷燃料颗粒辐照肿胀的力学解析模型 计算基体所受肿胀压力pf,用肿胀压力替代燃料颗粒的作用,使金属基体在肿胀压力pf载荷下变形失效.基体所受陶瓷核燃料球体的辐照肿胀压力pf按式()进行计算.pfEfvf(vf)fvf(pTpvfEf)(kBU)/()其中:vf和Ef分别为陶瓷燃料的泊松比和弹性模量(MP a);k为 常 系 数;BU为 燃 耗 深 度,F I MA;f为燃料 颗粒气 孔 率;T为 温 度,;p为气孔内压,MP a,可根据式()求解.pngR TVg(ngaVgngb)()其中:ng、Vg和R分别为裂变气体浓度、裂变气体体积和气体常数;a和b为实际气体状态方程参数.单颗粒陶瓷燃料中裂变气体的总浓度ng通常满足:ng RfDf BU/Mf()其中:Df为燃料颗粒密度,g/c m;为裂变气体产额;Mf为陶瓷燃料的摩尔质量,g/m o l.UO陶瓷具有耐高温、耐辐照、肿胀率低、包容裂变产物性能好等优点,广泛用于设计和生产弥散型燃料元件.另一方面,具有良好力学性能的核级 不锈钢(S S)是反应堆的关键结构材料,可用作弥散型燃料元件的金属基体和包壳材料.针对UO燃料 S S基体这种弥散型核燃料芯体结构,基于式()所示的肿胀压力pf表达 式,采用UO陶瓷 燃料和 S S的材料参数对pf值进行解析计算和影响因素分析.温度和燃耗对UO陶瓷和 S S的弹性模量和泊松比的影响关系 如下:vm(T)T()Em(T)T T()Ef(T,BU)/(BU)(T )()其中,vm和Em分别为 S S的泊松比和弹性模量(MP a).S S基体和燃料颗粒间的损伤层材料参数取值参考文献 中的研究结果,不考虑韧性,弹性模量约为基体材料 S S的倍,泊松比约为基体材料的 倍.在确定了非团聚区域等效材料常数和团聚区域的燃料颗粒建模方法后,为快速生成具有不同团聚程度的弥散型燃料元件代表性体元模型,开发了含颗粒团聚行为的R V E模型的参数化建模脚本,程序流程如图所示.在图所示的建模流程中,通过输入陶瓷燃料颗粒燃耗、辐照温度、燃料颗粒团聚体积分数和环境压力等参数,可根据式()()和式()()直接计算不同辐照条件下的材料属性和非团聚区域的等效材料性能,根据式()中团聚颗粒间距dt计算结果生成团聚区域的燃料颗粒几何构型,设置损伤层和金属基体间的接触,并将根据式()得到的燃料颗粒辐照肿胀压力解析值pf和环境压力q作为载荷条件施加在模型上.约束与边界条件随R V E模型在弥散型核燃料元件中代表位置的不同而改变.最终通过调控关键参数实现具有不同燃料颗粒团聚程度R V E模型的程序化数值计算.数值结果与讨论 团聚体积分数的影响燃料颗粒团聚行为主要涉及颗粒的团聚程度和团聚区在弥散型核燃料芯体的分布位置,其中颗粒的团聚程度通过团聚颗粒间距确定.固定陶瓷燃料颗粒的燃耗和温度分别为 F I MA和 ,非团聚区域的燃料颗粒体积分数为,弥散型燃料芯体强度分析结果如图所示.根据计算结果,分别探讨燃料颗粒团聚体积分数、环境压力以及团聚区域位置对芯体失效的影响.原子能科学技术第 卷图含颗粒团聚行为的R V E模型参数化建模流程F i g S c h e m a t i cd i a g r a mf o rp a r a m e t r i cm o d e l i n go fR V Em o d e lw i t hp a r t i c l ea g g l o m e r a t i o nb e h a v i o r图弥散型核燃料芯体强度分析结果F i g S t r e n g t ha n a l y s i sr e s u l to fd i s p e r s i o nn u c l e a r f u e lm e a t首先,对于代表性体元位于弥散型燃料元件中心的情况,研究燃料颗粒团聚体积分数对芯体中最大M i s e s应力 的影响,计算 结果如图 a所示.在不同燃料颗粒团聚体积分数Vf下,最大M i s e s应力随环境压力q的增加而减小,但不同曲线的斜率不同,即不同团聚体积分数下最大M i s e s应力的变化率不同.燃料颗粒团聚体积分数的变化范围在(,)区间,环境压力q从 MP a增加至 MP a.根据计算结果可知,弥散型核燃料芯体中的最大M i s e s应力随燃料颗粒团聚体积分数的增加而增大.但当燃料颗粒团聚体积分数大于 时,M i s e s应力基本不变.故燃料颗粒团聚是弥散型核燃料芯体强度分析的短板.而环境压力会抑制颗粒团聚区域的应力集中,使得芯体的最大M i s e s应力减小,弥散型核燃料芯体的失效概率降低.但燃料颗粒团聚体积分数越大的模型,最大M i s e应力改变量越小,即随着燃料颗粒团聚程度的增加,芯体中应力状态对环境压力的敏感性降低.由于反应堆中冷却剂循环流动,燃料元件的包壳外不可避免地存在冷却剂诱导的环境压力.环境压力对芯体中最大M i s e s应力的影响如图 b所示.随着团聚区颗粒间距的增大,燃料颗粒团聚体积分数变大,基体中最大M i s e s应力随环境压力的改变量非线性减小.当环境压力由 MP a增加到 MP a时,燃料颗粒团聚R V E模型中基体最大M i s e s应力第期董颖璇等:颗粒团聚行为对弥散型核燃料芯体失效的影响分析分别降低约 MP a和 MP a.这是由于颗粒间距更小的模型应力集中程度更强,对环境压力更为敏感.环境压力实际上是对弥散型核燃料芯体的力学约束.但由于环境压力的约束作用与颗粒间距减小的促进作用间相互竞争的影响,不同燃料颗粒团聚程度下芯体最大M i s e s应力随环境压力的变化率不同.故冷却剂产生的环境压力将提高芯体的安全性,最大M i s e s应力随环境压力的增大而减小.因此,冷却剂对芯体应力状态的抑制作用对降低弥散型核燃料元件的团聚短板效应具有重要意义.团聚区域位置的影响给定燃料颗粒团聚体积分数为 ,非团聚区域的燃料颗粒体积分数为 ,分析燃料颗粒团聚区域位置变化对芯体最大M i s e s应力的影响规律,结果如图所示.a 两种极限位置示意图;b 极限位置下最大M i s e s应力随环境压力的变化图团聚区域距离包壳位置对芯体中最大M i s e s应力的影响F i g E f f e c to fd i s t a n c e f r o ma g g l o m e r a t i o n l o c a t i o nt oc l a d d i n go nm a x i m u m M i s e ss t r e s s i nm a t r i x该 问 题 将 分 两 种 情 况 进 行 分 析:)当R V E模型位于弥散型核燃料芯体中心位置(即只具有上下包壳的结构)时,燃料颗粒团聚区在上下包壳面间沿z轴的位置变化;)当燃料颗粒团聚区位于代表性体元中心位置时,考虑R V E模型在 整个弥散型 燃料元件不 同位置变化.首先,针对R V E模型位于弥散型核燃料芯体中心的情况,研究燃料颗粒团聚区沿z轴位置变化对芯体强度的影响.如图 a所示,在燃料元件的代表性体元中,团聚区域存在两种极限位置:)位置,R V E模型中间,最远离包壳的位置;)位置,最接近包壳的位置.通过计算不同燃料颗粒团聚位置下最大M i s e s应力 随 环 境 压 力 的 变 化(图 b)可 知,最 大M i s e s应力随环境压力的增大而减小,不同团聚位置下最大M i s e s应力的变化趋势相同.故燃料颗粒团聚区域在弥散型核燃料元件上下包壳中的位置对芯体中的应力状态影响不大.其次,针对燃料颗粒团聚区位于代表性体元中心位置不变的情况,研究R V E模型在整个弥散型核燃料元件不同位置(即团聚燃料颗粒在x y面中位置变化)时的芯体失效行为.R V E模型在弥散型核燃料元件中的位置有种类型,如图 a所示.种典型位置下的几何模型不变,通过不同边界条件定义种位置.位置处于芯体中间位置,R V E除包壳上下表面受环境压力外,芯体四周施加对称性边界条件;位置为单边毗邻包壳,施加单向位移约束边界条件,另外个毗邻芯体的表面施加对称性边界条件;位置为弥散型核燃料元件的边角位置,两个相邻表面施加单向位移约束边界条件,另外两个毗邻芯体的表面施加对称性边界条件约束.个位置R V E模型的最大M i s e s应力计算结果如图 b所示.根据计算结果,位置和位置中毗邻包壳位置的弥散型核燃料芯体的最大M i s e s应力基本相等.无环境压力作用时,位置的M i s e s应力最大.随着环境压力的增加,所有位置的最大M i s e s应力均减 小,但不同位 置的最大M i s e s应力变化率不同.弥散型核燃料芯体中心位置(位置)的颗粒团聚对环境压力更敏原子能科学技术第 卷感,最大M i s e s应力随环境压力的改变量最大.位置和位置的最大M i s e s应力在环境压力作用下的改变量较小,环境压力对毗邻包壳位置的R V E模型的影响相似.根据强度理论,弥散型核燃料芯体中心位置的颗粒团聚位置更容易发生失效.图R V E模型在弥散型核燃料元件中不同位置对芯体最大M i s e s应力的影响F i g E f f e c to f a g g l o m e r a t i o nl o c a t i o ni nd i s p e r s i o nn u c l e a r f u e l e l e m e n to nm a x i m u m M i s e ss t r e s s i nm a t r i x当环境压力为 MP a时,弥散型核燃料芯体x y面的中心位置(只有上下包壳约束)最大M i s e s应力最大;当环境压力为 MP a时,弥散型 核 燃 料 芯 体x y面 的 角 落 位 置 最 大M i s e s应力最大.故弥散型核燃料芯体的失效行为是环境压力和燃料颗粒团聚分布位置共同影响的结果.当环境压力小于 MP a时,位置的最大M i s e s应力最大;但当环境压力大于 MP a时,位置和位置的最大M i s e s应力均大于位置.这说明在弥散型核燃料元件中,芯体中心位置的应力场对环境压力更敏感,在冷却剂诱导环境压力作用下其安全性提升最多.当环境压力超过 MP a时,整个燃料元件中不同位置的局部M i s e s应力集中将被均匀化,抑制到相似的水平.故反应堆一回路中冷却剂产生的环境静水压力可以降低颗粒团聚区域的M i s e s应力场,提高弥散型核燃料芯体的失效条件,增加反应堆安全性.而燃料颗粒团聚行为是导致弥散型核燃料芯体失效的关键因素,综合分析图可知,当颗粒团聚程度大于 时,环境静水压力对芯体局部失效将不再有抑制作用,尤其是当团聚位置在图 a所示的毗邻包壳位置时.燃料/基体力学性质的影响将燃料颗粒弹性模量Ef固定,金属基体弹性模量Em作为参考变量,讨论金属基体材料性质对弥散型核燃料芯体失效的影响.通常,陶瓷燃料颗粒的弹性模量大于金属基体,令基体弹性模量Em的取值范围在陶瓷燃料弹性模量Ef的 倍之间变化.图为弥散型燃料芯体中最大M i s e s应力随基体弹性模量Em的变化.由图可见,芯体中