针对钠冷快堆单根组件热变形的模拟方法技术

技术编号:20623013 阅读:48 留言:0更新日期:2019-03-20 14:37
一种预测钠冷快堆单根组件热变形的模拟方法,1、确定组件的三维温度场数据、几何参数以及物性参数,2按二维数组排列三维温度场数据并进行合理差值,3、对组件的横截面进行离散化处理,将其离散为12个计算单元,4、求解每一个离散单元的静矩,5、求解每一个离散单元的热应变,推算每个离散单元的热弯矩以及该轴向节点的总热弯矩,6、推算每个轴向节点的横向热变形位移,7、推算该轴向节点的热轴力,8、推算每个轴向节点的轴向热变形位移,9、获得组件在三维空间内的热变形情况。

Simulation Method for Thermal Deformation of Single Module of Sodium Cooled Fast Reactor

A simulation method for predicting thermal deformation of single component of sodium cooled fast reactor is presented. 1. Determining three-dimensional temperature field data, geometric parameters and physical parameters of the component, 2 arranging three-dimensional temperature field data according to two-dimensional array and making reasonable difference, 3. Discrete the cross section of the component into 12 calculation units, 4. Solving the static moment of each discrete element, 5. Solving each discrete element. The thermal strain of each discrete element is calculated. The thermal bending moment of each discrete element and the total thermal bending moment of the axial node are calculated. 6. The transverse thermal deformation displacement of each axial node is calculated. 7. The thermal axial force of the axial node is calculated. 8. The axial thermal deformation displacement of each axial node is calculated. 9. The thermal deformation of the component in three-dimensional space is obtained.

【技术实现步骤摘要】
针对钠冷快堆单根组件热变形的模拟方法
本专利技术涉及钠冷快堆单根组件热变形
,具体涉及一种针对钠冷快堆单根组件热变形的模拟方法。
技术介绍
钠冷快堆中的组件在整个钠冷快堆的运行过程中占有非常重要的地位。在钠冷快堆的运行过程中,燃料组件的活性区中会发生裂变反应释放大量的热量,而燃料组件周围充当冷却剂的钠会通过流动的方式将燃料组件内裂变反应释放的热量带出堆芯外。在我国自主设计研发的钠冷快堆中,由于中子通量分布不均、燃料富集度差异以及冷井等多重因素的影响,快堆中组件释热功率分布是不均匀的,并且在堆芯内部不同组件区域的温度差异非常大。对于堆芯中的某一根组件而言,由于其与毗邻组件的间隙非常小,因此该组件的温度场分布会受到毗邻组件温度场分布的影响,导致在该组件的不同高度以及组件的不同周向方向上温度分布也会有较大的区别。在不均匀温度场的作用下,组件局部的热膨胀特性会出现较大的差异,从而使得组件发生热弯曲变形。组件发生热弯曲变形后就会通过垫块与周围的毗邻组件发生接触。如果弯曲变形量较大的话随之产生的接触力也会变大,这样会给钠冷快堆的安全运行带来不利影响,如增加组件插拔力、组件倒换困难、加速组件破损、堆芯象限功率倾斜等。由于钠冷快堆中的组件变形现象在快堆安全分析中占有非常重要的地位,并且组件的弯曲变形通常在毫米的量级,弯曲变形量通过实验的方式难以进行测量。基于上述原因,需要研发钠冷快堆组件的热变形模拟技术来预测组件在快堆堆芯中的热变形行为,从而为快堆组件的安全分析提供充分的理论依据。由于我国钠冷快堆组件在结构设计上具有特殊性,因此必须开发出适用于我国钠冷快堆的组件热变形分析的模拟技术来满足相应的快堆组件安全分析。
技术实现思路
为了克服上述现有技术存在的问题,本专利技术的目的在于提供针对钠冷快堆单根组件热变形的模拟方法,通过组件的热力耦合本构关系式来预测组件在不均匀三维温度场分布中的热变形情况。为了实现上述目的,本专利技术采取了以下技术方案:一种针对钠冷快堆单根组件热变形的模拟方法,该方法包括以下步骤:步骤1:根据用户需求,输入目标组件的物性参数以及几何参数,物性参数包括组件的杨氏模量以及热膨胀系数,几何参数包括组件的长度,套管厚度以及外套管的对边距;步骤2:输入目标组件的三维温度场数据:三维温度场数据以矩阵的形式给出,其中列标为组件轴向不同位置温度采样点的标高,行标为在组件某一固定标高处进行温度采点的具体位置标记;获得三维温度场数据以后采用插值方法获得轴向均匀分布的温度场数据从而满足计算的输入条件;步骤3:对于目标组件的截面进行离散化处理,将其分为12个关于z轴对称的单元,并通过下式获得12个离散的截面单元的静矩即一阶矩;其中:A1到A12为每一个离散单元的横截面积;z为离散单元的积分域到z轴的距离;由于组件的截面为中心对称图形,因此在同一横截面上各离散单元之间存在如上所示的关系,因此只需求解相邻的三个离散单元的静矩值即可通过几何关系推得所有12个离散单元的静矩值;其中相邻的三个离散单元静矩的计算方法如下所示:其中:a1为套管外对边距,a2为套管内对边距;y为离散单元的积分域到y轴的距离;y1和y2分别为离散单元7的左、右积分边界;y3和y4分别为离散单元8的右侧子单元的左、右积分边界;y5为离散单元9的右积分边界;步骤4:获得所有离散单元的静矩以后通过下式求解离散单元截面的总热弯矩;其中:MT为离散单元所在截面的总热弯矩;xi为离散单元的轴向坐标;E为杨氏模量;z为离散单元的积分域到z轴的距离;εT为轴向位置为xi的横截面的热应变;αi为轴向位置为xi的横截面的热膨胀系数;i为轴向坐标的索引号;j为同一截面上离散单元的索引号Tf为参考温度;Ti为轴向位置为xi的横截面的热态温度;Tij为每一个离散单元的热态温度;A为套管的横截面积;为离散单元的静矩即一阶矩;然后遍历所有的轴向节点,求得所有轴向节点所在截面的总热弯矩;步骤5:通过下式获得每一个轴向节点的横向热变形位移即挠度;其中:MT为离散单元截面的总热弯矩;w为横向热变形位移即挠度;xi为离散单元的轴向坐标;EI为抗弯刚度;h为步长即轴向节点间距;步骤6:通过下式求解截面的热轴力;其中:NT为截面的热轴力;xi为离散单元的轴向坐标;E为杨氏模量;εT为轴向位置为xi的横截面的热应变;αij为每一个离散单元的热膨胀系数;Tf为参考温度;Tij为每一个离散单元的热态温度;A为套管的横截面积;为离散单元的横截面积;然后遍历所有的轴向节点,求得所有轴向节点所在截面的热轴力;步骤7:通过下式获得每一个轴向节点的轴向热变形位移;其中:EA为抗拉刚度;u为轴向热变形位移;xi为离散单元的轴向坐标;NT为热轴力;f为重力项;h为步长即轴向节点间距;步骤8:整理不同方向的热变形位移信息,获得组件在三维空间内的热变形情况。与现有技术相比,本专利技术有如下突出特点:1.本专利技术将组件的横截面离散为12个单元进行求解,因此相比于传统方法精度更高;2.通过对组件的热轴力以及热弯矩的求解可以准确获得组件在三维空间内的热变形情况;3.本专利技术模拟方法基于Euler-Bernoulli梁理论进行构建,因此相比于传统的有限元方法具有计算效率更高,计算结果更加准确的特点;4.本专利技术模拟方法不仅仅适用钠冷快堆组件,针对于不同的计算对象以及物理问题,本专利技术模拟方法可以在该基础之上进行扩展。附图说明图1为组件横截面离散单元示意图。具体实施方式下面结合附图和具体实施方式对本专利技术作进一步详细说明:本专利技术一种针对钠冷快堆单根组件热变形的模拟方法,通过组件的热力耦合本构关系式来预测组件在不均匀三维温度场分布中的热变形情况。该方法具体流程包括以下方面:步骤1:根据用户需求,输入目标组件的物性参数以及几何参数,物性参数包括组件的杨氏模量以及热膨胀系数,几何参数包括组件的长度,套管厚度以及外套管的对边距等。组件的温度场数据、几何参数以及物性参数为该模拟方法的必要输入条件,有了这些参数以后便可以进行热变形的模拟计算工作。这些参数可以通过实验数据获得,也可以通过程序、软件的计算以及文献调研获得。步骤2:输入目标组件的三维温度场数据:三维温度场数据以矩阵的形式给出,其中列标为组件轴向不同位置温度采样点的标高,行标为在组件某一固定标高处进行温度采点的具体位置标记(每一层标高共有12个采样点,包括组件的六条棱以及六个面的中心位置)。获得三维温度场数据以后采用插值方法获得轴向均匀分布的温度场数据从而满足程序模拟的输入条件。所有的温度点数据均会保存在预先设置好的空文件中(一般为txt格式),等所有温度点数据记录完整后,会自动调用差值子程序,从而通过差值方法补全所缺的温度点数据,并自动生成新的文件将所有温度点数据记录在其中。步骤3:对于目标组件的截面进行离散化处理,将其分为12个单元,并通过下式获得12个离散的截面单元的静矩(一阶矩);其中:A1到A12为每一个离散单元的横截面积;z为离散单元的积分域到z轴的距离;参考图1中的组件横截面离散单元示意图可知,由于组件的截面为中心对称图形,因此关于z轴对称的面积区域的静矩积分在z轴两侧相等,而关于y轴对称的面积区域的静矩积分在y轴两侧互为相反数。综上所述,同一横截面上各离散单元之间存在如上式所本文档来自技高网
...

【技术保护点】
1.一种针对钠冷快堆单根组件热变形的模拟方法,其特征在于:包括以下步骤:步骤1:根据用户需求,输入目标组件的物性参数以及几何参数,物性参数包括组件的杨氏模量以及热膨胀系数,几何参数包括组件的长度,套管厚度以及外套管的对边距;步骤2:输入目标组件的三维温度场数据:三维温度场数据以矩阵的形式给出,其中列标为组件轴向不同位置温度采样点的标高,行标为在组件某一固定标高处进行温度采点的具体位置标记;获得三维温度场数据以后采用插值方法获得轴向均匀分布的温度场数据从而满足计算的输入条件;步骤3:对于目标组件的截面进行离散化处理,将其分为12个关于z轴对称的单元,并通过下式获得12个离散的截面单元的静矩即一阶矩;

【技术特征摘要】
1.一种针对钠冷快堆单根组件热变形的模拟方法,其特征在于:包括以下步骤:步骤1:根据用户需求,输入目标组件的物性参数以及几何参数,物性参数包括组件的杨氏模量以及热膨胀系数,几何参数包括组件的长度,套管厚度以及外套管的对边距;步骤2:输入目标组件的三维温度场数据:三维温度场数据以矩阵的形式给出,其中列标为组件轴向不同位置温度采样点的标高,行标为在组件某一固定标高处进行温度采点的具体位置标记;获得三维温度场数据以后采用插值方法获得轴向均匀分布的温度场数据从而满足计算的输入条件;步骤3:对于目标组件的截面进行离散化处理,将其分为12个关于z轴对称的单元,并通过下式获得12个离散的截面单元的静矩即一阶矩;其中:A1到A12为每一个离散单元的横截面积;z为离散单元的积分域到z轴的距离;由于组件的截面为中心对称图形,因此在同一横截面上各离散单元之间存在如上所示的关系,因此只需求解相邻的三个离散单元的静矩值即可通过几何关系推得所有12个离散单元的静矩值;其中相邻的三个离散单元静矩的计算方法如下所示:其中:y为离散单元的积分域到y轴的距离;a1为套管外对边距,a2为套管内对边距;y1和y2分别为离散单元7的左、右积分边界;y3和y4分别为离散单元8的右侧子单元的左、右积分边界;y5为离散单元9的右积分边界;步骤4:获得所有离散单元的静矩以后通过下式求解离散单元截面的总热弯矩;其中:MT为离散单元所在截面的总热弯矩;xi为离散单元的轴向坐标;E为杨氏模量;z为离散单元的积分域到z轴的距离...

【专利技术属性】
技术研发人员:苏光辉马振辉马泽华巫英伟田文喜秋穗正
申请(专利权)人:西安交通大学
类型:发明
国别省市:陕西,61

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1