[发明专利]针对钠冷快堆单排组件耦合热变形的模拟方法

说明书

一种针对钠冷快堆单排组件耦合热变形的模拟方法，1、确定单排组件的组件个数、目标组件的物性参数以及几何参数，2、将单排组件按照由堆芯朝堆外的方向用数字编号进行标记，3、获取每一根组件在所有轴向位置的自由热变形位移，4、对每一根组件在上、下垫块处的热变形位移进行取值，5、推导出四种不同的接触力与垫块之间的位移模式，6、列出每一根组件在上、下垫块处的平衡方程，7、列出与所有组件相关的变形协调方程，8、将组件的位移参数进行合并处理，9、获得上、下垫块位置的总位移；10、求解每相邻两根组件之间上、下垫块的接触力，11、通过组件的位移叠加原则求得每一根组件的最终位移；本发明方法具备计算量小，收敛性好、计算结果准确的优点。

技术领域

本发明涉及钠冷快堆单排组件耦合热变形技术领域，具体涉及一种针对钠冷快堆单排组件耦合热变形的模拟方法。

背景技术

钠冷快堆中的组件在整个钠冷快堆的运行过程中占有非常重要的地位。在钠冷快堆的运行过程中，燃料组件的活性区中会发生裂变反应释放大量的热量，而燃料组件周围充当冷却剂的钠会通过流动的方式将燃料组件内裂变反应释放的热量带出堆芯外。在我国自主设计研发的钠冷快堆中，由于中子通量分布不均、燃料富集度差异以及堆芯冷井等多重因素的影响，快堆中组件释热功率分布是不均匀的，并且在堆芯内部不同组件区域的温度差异非常大。对于堆芯中的某一根组件而言，由于其与毗邻组件的间隙非常小，因此该组件的温度场分布会受到毗邻组件温度场分布的影响，导致在该组件的不同高度以及组件的不同周向方向上温度分布也会有较大的区别。在不均匀温度场的作用下，组件局部的热膨胀特性会出现较大的差异，从而使得组件发生热弯曲变形。组件发生热弯曲变形后就会通过垫块与周围的毗邻组件发生接触。如果弯曲变形量较大的话随之产生的接触力也会变大，这样会给钠冷快堆的安全运行带来不利影响，如增加组件插拔力、组件倒换困难、加速组件破损、堆芯象限功率倾斜等。

在钠冷快堆中，中子源组件位于堆芯的中部，因此靠近堆芯的环形区域会有较大的中子通量，因而会释放更多的裂变能，从而该区域的组件会具有较高的温度。越靠近堆芯的外围，中子通量降低，裂变反应释放的裂变能降低，因而外围区域的组件温度也会相应降低。并且这种从堆芯中部到堆外的温度分布趋势沿着堆芯的周向方向上存在着一定的对称性。因此可以独立取出从堆芯中部到堆芯外围的一排组件进行研究，并且会发现这排组件的温度呈现出由堆芯向堆外逐渐降低的趋势。由于温度的不均匀分布，这排组件的会朝着温度降低的方向发生热弯曲变形，并且由于组件之间的距离非常近，组件与组件之间会发生接触，并产生很大的接触力。此外组件之间的相互接触也会阻碍组件原本的热弯曲变形，使得组件最终的变形状态更加难以预测。

由于单排组件的受限热变形过程是一个热力强耦合的过程，因此采用传统的有限元分析软件来分析该过程计算量巨大，且极易发散，而单根组件自由热变形分析程序又无法模拟及预测单排组件热变形的动态耦合过程。基于上述原因，有必要研发钠冷快堆单排组件的耦合热变形模拟技术来预测组件在快堆堆芯中的热变形行为，从而为快堆组件的安全分析提供充分的理论依据。

发明内容

为了克服上述现有技术存在的问题，本发明的目的在于提供针对钠冷快堆单排组件的耦合热变形的模拟方法，通过求解组件的变形协调方程组来预测单排组件在热梯度作用下的耦合变形行为。

为了实现上述目的，本发明采取了以下技术方案：

一种针对的钠冷快堆单排组件耦合热变形的模拟方法，该方法包括以下步骤：

步骤1：根据用户需求，输入单排组件的组件个数、目标组件的物性参数以及几何参数，物性参数指组件的杨氏模量和热膨胀系数，几何参数指组件的长度、套管厚度、外套管的对边距、组件的上下垫块到组件固支端的距离、相邻组件之间的轴心距以及相邻组件上下垫块之间的间隙；

步骤2：将单排组件按照由堆芯朝堆外的方向用数字编号进行标记，并以最靠近堆芯的组件位置为零点，根据相邻组件的轴心距确定单排组件中每一根组件在弯曲的一维方向的初始位置坐标；