CN109492332B - 针对钠冷快堆单根组件热变形的模拟方法 - Google Patents
针对钠冷快堆单根组件热变形的模拟方法 Download PDFInfo
- Publication number
- CN109492332B CN109492332B CN201811513849.3A CN201811513849A CN109492332B CN 109492332 B CN109492332 B CN 109492332B CN 201811513849 A CN201811513849 A CN 201811513849A CN 109492332 B CN109492332 B CN 109492332B
- Authority
- CN
- China
- Prior art keywords
- discrete unit
- axial
- section
- component
- discrete
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
一种预测钠冷快堆单根组件热变形的模拟方法,1、确定组件的三维温度场数据、几何参数以及物性参数,2按二维数组排列三维温度场数据并进行合理差值,3、对组件的横截面进行离散化处理,将其离散为12个计算单元,4、求解每一个离散单元的静矩,5、求解每一个离散单元的热应变,推算每个离散单元的热弯矩以及该轴向节点的总热弯矩,6、推算每个轴向节点的横向热变形位移,7、推算该轴向节点的热轴力,8、推算每个轴向节点的轴向热变形位移,9、获得组件在三维空间内的热变形情况。
Description
技术领域
本发明涉及钠冷快堆单根组件热变形技术领域,具体涉及一种针对钠冷快堆单根组件热变形的模拟方法。
背景技术
钠冷快堆中的组件在整个钠冷快堆的运行过程中占有非常重要的地位。在钠冷快堆的运行过程中,燃料组件的活性区中会发生裂变反应释放大量的热量,而燃料组件周围充当冷却剂的钠会通过流动的方式将燃料组件内裂变反应释放的热量带出堆芯外。在我国自主设计研发的钠冷快堆中,由于中子通量分布不均、燃料富集度差异以及冷井等多重因素的影响,快堆中组件释热功率分布是不均匀的,并且在堆芯内部不同组件区域的温度差异非常大。对于堆芯中的某一根组件而言,由于其与毗邻组件的间隙非常小,因此该组件的温度场分布会受到毗邻组件温度场分布的影响,导致在该组件的不同高度以及组件的不同周向方向上温度分布也会有较大的区别。在不均匀温度场的作用下,组件局部的热膨胀特性会出现较大的差异,从而使得组件发生热弯曲变形。组件发生热弯曲变形后就会通过垫块与周围的毗邻组件发生接触。如果弯曲变形量较大的话随之产生的接触力也会变大,这样会给钠冷快堆的安全运行带来不利影响,如增加组件插拔力、组件倒换困难、加速组件破损、堆芯象限功率倾斜等。
由于钠冷快堆中的组件变形现象在快堆安全分析中占有非常重要的地位,并且组件的弯曲变形通常在毫米的量级,弯曲变形量通过实验的方式难以进行测量。基于上述原因,需要研发钠冷快堆组件的热变形模拟技术来预测组件在快堆堆芯中的热变形行为,从而为快堆组件的安全分析提供充分的理论依据。
由于我国钠冷快堆组件在结构设计上具有特殊性,因此必须开发出适用于我国钠冷快堆的组件热变形分析的模拟技术来满足相应的快堆组件安全分析。
发明内容
为了克服上述现有技术存在的问题,本发明的目的在于提供针对钠冷快堆单根组件热变形的模拟方法,通过组件的热力耦合本构关系式来预测组件在不均匀三维温度场分布中的热变形情况。
为了实现上述目的,本发明采取了以下技术方案:
一种针对钠冷快堆单根组件热变形的模拟方法,该方法包括以下步骤:
步骤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轴两侧互为相反数。综上所述,同一横截面上各离散单元之间存在如上式所示的关系,因此只需求解相邻的三个离散单元的静矩值即可通过几何关系推得所有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为参考温度;
Tij为每一个离散单元的热态温度;
Ti为轴向位置为xi的横截面的热态温度;
A为离散单元的横截面积;
为离散单元的静矩即一阶矩;
从上式可以看出,组件的总热弯矩主要与组件的杨氏模量、离散单元的静矩还有离散单元的热应变相关。由于离散单元的静矩已经求得,因此求解总热弯矩的核心在于求解离散单元的热应变。而离散单元的热应变是单元温度场以及热膨胀系数的函数,因此在可以调用组件的三维温度场数据,进而通过参考温度推算每个离散单元的温变,最终获得每个离散单元的热应变。
获得某一轴向节点的总热弯矩以后可以遍历所有的轴向节点,求得所有轴向节点所在截面的总热弯矩。
步骤5:通过下式获得每一个轴向节点的横向热变形位移(挠度);
其中:
MT为离散单元截面的总热弯矩;
w为横向热变形位移即挠度;
xi为离散单元的轴向坐标;
EI为抗弯刚度;
h为步长即轴向节点间距。
该式说明任意内节点的横向位移与相邻接点的横向位移以及该节点的总热弯矩之间存在强烈相关性,因此可以根据上式遍历所有的轴向节点获得所有轴向节点的横向热变形位移。
步骤6:通过下式求解截面的热轴力;
其中:
NT为截面的热轴力;
xi为离散单元的轴向坐标;
E为杨氏模量;
εT为轴向位置为xi的横截面的热应变;
αij为每一个离散单元的热膨胀系数;
Tf为参考温度;
Tij为每一个离散单元的热态温度;
A为套管的横截面积;
为离散单元的横截面积;
从上式可以看出,与热弯矩类似,组件的热轴力同样主要与组件的杨氏模量、还有离散单元的热应变相关,但是由于热轴力的积分中不含有z项,因此不需要求解离散单元的静矩,而是直接获得离散单元的面积即可。热应变在热弯矩求解中已求得,在此不再进行赘述。
获得某一轴向节点的热轴力以后可以遍历所有的轴向节点,求得所有轴向节点所在截面的热轴力。
步骤7:通过下式获得每一个轴向节点的轴向热变形位移:
其中:
EA为抗拉刚度;
u为轴向热变形位移;
xi为离散单元的轴向坐标;
NT为热轴力;
f为重力项;
h为步长即轴向节点间距。
该式说明任意内节点的轴向位移与相邻接点的轴向位移以及该节点的热轴力之间存在强烈相关性,因此可以根据上式遍历所有的轴向节点获得所有轴向节点的轴向位移。
步骤8:整理不同方向的热变形位移信息。通过整合每一个轴向节点在不同方向的位移矢量来获得组件在三维空间内的热变形情况。
Claims (2)
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轴的距离;
εT为轴向位置为xi的横截面的热应变;
α i 为轴向位置为xi的横截面的热膨胀系数;
i为轴向坐标的索引号;
j为同一截面上离散单元的索引号;
α ij 为每一个离散单元的热膨胀系数;
Tf为参考温度;
Ti为轴向位置为xi的横截面的热态温度;
Tij为每一个离散单元的热态温度;
A为套管的横截面积;
为离散单元的静矩即一阶矩;
然后遍历所有的轴向节点,求得所有轴向节点所在截面的总热弯矩;
步骤5:通过下式获得每一个轴向节点的横向热变形位移即挠度;
其中:
MT为离散单元截面的总热弯矩;
w为横向热变形位移即挠度;
xi为离散单元的轴向坐标;
EI为抗弯刚度;
h为步长即轴向节点间距;
步骤6:通过下式求解截面的热轴力;
其中:
NT为截面的热轴力;
xi为离散单元的轴向坐标;
E为杨氏模量;
εT为轴向位置为xi的横截面的热应变;
α ij为每一个离散单元的热膨胀系数;
Tij为每一个离散单元的热态温度;Tf为参考温度;
A为套管的横截面积;
为离散单元的横截面积;
然后遍历所有的轴向节点,求得所有轴向节点所在截面的热轴力;
步骤7:通过下式获得每一个轴向节点的轴向热变形位移;
其中:
EA为抗拉刚度;
u为轴向热变形位移;
xi为离散单元的轴向坐标;
NT为热轴力;
f为重力项;
h为步长即轴向节点间距;
步骤8:整理不同方向的热变形位移信息,获得组件在三维空间内的热变形情况。
2.根据权利要求1所述的一种针对钠冷快堆单根组件热变形的模拟方法,其特征在于:步骤2所述的具体位置标记为每一层标高共有12个采样点,其中6个采样点在组件六条棱的位置,剩余6个采样点在组件六个面的中心位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811513849.3A CN109492332B (zh) | 2018-12-11 | 2018-12-11 | 针对钠冷快堆单根组件热变形的模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811513849.3A CN109492332B (zh) | 2018-12-11 | 2018-12-11 | 针对钠冷快堆单根组件热变形的模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109492332A CN109492332A (zh) | 2019-03-19 |
CN109492332B true CN109492332B (zh) | 2019-08-27 |
Family
ID=65698408
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811513849.3A Active CN109492332B (zh) | 2018-12-11 | 2018-12-11 | 针对钠冷快堆单根组件热变形的模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109492332B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108345741A (zh) * | 2017-12-13 | 2018-07-31 | 湘潭大学 | 基于无网格rkpm的各向异性材料结构二维热变形和热应力分析方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10163534B2 (en) * | 2014-01-27 | 2018-12-25 | Terrapower, Llc | Modeling for fuel element deformation |
CN106484984B (zh) * | 2016-09-29 | 2019-11-12 | 北京大学 | 星载柔性附件热致微振动响应仿真分析平台 |
CN107301265B (zh) * | 2017-05-11 | 2021-06-11 | 北京空间飞行器总体设计部 | 基于关节热应变释放冲击的航天器动力学建模方法 |
CN108334661B (zh) * | 2017-12-29 | 2021-08-06 | 武汉华中数控股份有限公司 | 一种进给轴热变形预测方法 |
CN108595827A (zh) * | 2018-04-20 | 2018-09-28 | 安徽工业大学 | 一种C-Mn-Al高强度钢热变形组织演变机制及热加工性能的确定方法 |
-
2018
- 2018-12-11 CN CN201811513849.3A patent/CN109492332B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108345741A (zh) * | 2017-12-13 | 2018-07-31 | 湘潭大学 | 基于无网格rkpm的各向异性材料结构二维热变形和热应力分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109492332A (zh) | 2019-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111291494B (zh) | 用于核反应堆triso燃料颗粒的多尺度多物理场耦合模拟方法 | |
Saxena | Synthesis of compliant mechanisms for path generation using genetic algorithm | |
Wang et al. | Two-stage size-layout optimization of axially compressed stiffened panels | |
CN109063235A (zh) | 一种用于反应堆模拟的多物理耦合系统及方法 | |
Song et al. | Tokamak engineering mechanics | |
CN115270660B (zh) | 空间热离子反应堆瞬态行为多尺度多物理场耦合分析方法 | |
Sahu et al. | Analytical and semi-analytical models of conduction controlled rewetting: a state of the art review | |
CN114913936B (zh) | 一种针对铀钚混合氧化物燃料的多物理燃料性能分析方法 | |
Yu et al. | Preliminary coupling of the Thermal/Hydraulic solvers in the Monte Carlo code MCS for practical LWR analysis | |
CN109492332B (zh) | 针对钠冷快堆单根组件热变形的模拟方法 | |
He et al. | Development of a dynamics model for graphite-moderated channel-type molten salt reactor | |
Lu et al. | Coupling methods for HTR-PM primary circuit | |
CN109753691B (zh) | 针对钠冷快堆单排组件耦合热变形的模拟方法 | |
Luo et al. | Development and application of a multi-physics and multi-scale coupling program for lead-cooled fast reactor | |
Liangqi et al. | Response surface methodology-based hybrid robust design optimization for complex product under mixed uncertainties | |
Fanning et al. | Multi-resolution modeling of subassembly pin bundles for advanced fast reactor safety simulations | |
Esmaili et al. | Prediction of temperature distribution in annular fuels using orthogonal collocation method | |
Zhou et al. | A coupling analysis method of the thermal hydraulics and neutronics based on inverse distance weighted method | |
CN110781621B (zh) | 一种含几何形状缺陷薄壁结构承载响应的快速重分析方法 | |
Zhou | Operator split, Picard iteration and JFNK methods based on nonlinear CMFD for transient full core models in the coupling multiphysics environment | |
Fan et al. | An optimized CFD method for conceptual flow design of water cooled ceramic blanket | |
Han et al. | Development of a thermal hydraulic analysis code for gas-cooled reactors with annular fuels | |
Robert et al. | Impact of thermal coupling on a pebble bed reactor equilibrium from hyper-fidelity depletion | |
Gamble et al. | A feasibility study on the use of the MOOSE computational framework to simulate three-dimensional deformation of CANDU reactor fuel elements | |
Bindra et al. | A Computational-experimental Study to Simulate Mixing and Thermal Stratification in SFRs |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |