WO2020088375A1 - 燃料棒包壳蠕变有限长管修正方法的建立方法 - Google Patents
燃料棒包壳蠕变有限长管修正方法的建立方法 Download PDFInfo
- Publication number
- WO2020088375A1 WO2020088375A1 PCT/CN2019/113442 CN2019113442W WO2020088375A1 WO 2020088375 A1 WO2020088375 A1 WO 2020088375A1 CN 2019113442 W CN2019113442 W CN 2019113442W WO 2020088375 A1 WO2020088375 A1 WO 2020088375A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- cladding
- cladding tube
- tube
- length
- fuel rod
- Prior art date
Links
Images
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
Definitions
- the invention relates to the technical field of creep analysis of fuel rod cladding in a reactor, and in particular to a method for establishing a correction method of a fuel rod cladding creep finite length tube.
- the ceramic pellets will appear densely at the beginning of irradiation, and the axial height of the pellet column will become shorter, and the fuel rod cladding will become longer due to irradiation growth. Since the coolant pressure is greater than the internal pressure of the fuel rod and the creep behavior of the cladding, after the pellet contacts the cladding, a cladding pipe section without supporting the pellet may appear in the axial direction of the pellet column, that is, an unsupported cladding pipe section. Under the effect of external pressure and creep of the cladding, the finite length cladding tube that loses the core support creeps inward, and the cladding tube deformation continues to increase. The cladding deformation reaches a certain level or the cladding is flattened. The cladding creep collapsed.
- the corrected length of the unsupported cladding tube section in the fuel rod is limited, the two ends are supported by the core block, and the cladding creep model is obtained from the experiment on the empty cladding tube.
- the empty cladding tube used in the experiment is relatively
- the supported cladding pipe section is infinitely long, so the creep behavior of the unsupported pipe section cannot be calculated directly from the cladding creep model, but the cladding creep model needs to be corrected.
- the correction factor is called the finite length tube correction factor.
- the correction factor of the cladding creep finite length tube directly affects the cladding creep behavior, and the cladding creep is the core of the fuel rod cladding creep collapse analysis, so the finite length tube correction factor is the cladding creep analysis Important parameters directly affect the rationality and reliability of cladding creep analysis.
- the technical problem to be solved by the present invention is to provide a method for establishing a correction method for a fuel rod cladding creep finite length tube in view of the above-mentioned defects.
- the technical solution adopted by the present invention to solve its technical problems is to provide a method for establishing a fuel rod cladding creep finite length tube correction method, which includes the following steps:
- the parameters include the outer diameter, thickness, tube length, initial ovality, material type of the cladding tube, and the pressure difference P between the cladding tube and the cladding tube.
- stiffness K include limited stiffness when finite length of the cladding tube and the cladding tube is infinite stiffness K infinite;
- the initial ellipticity and the pressure difference P between the cladding tube and the cladding tube on the correction factor of the cladding tube determine the curve of the correction factor of the finite tube and establish the correction method of the finite tube.
- the physical property parameters of the cladding tube are determined according to the type of cladding tube material; the physical property parameters of the cladding tube include Young's modulus and Poisson's ratio.
- step S1 the pressure difference P inside and outside the cladding tube is equal everywhere.
- step S2 includes:
- step S2.1 includes:
- the cladding tube is also set for contact, and a non-slip, hard contact method is selected.
- the physical property parameters of the cladding tube include Young's modulus and Poisson's ratio; the interface property parameters of the cladding tube include the cladding tube interface model and thickness.
- the stiffness K is limited to the ratio of the internal and external pressure difference P of the cladding tube and the ovality of the cladding tube ty for a finite length;
- step S4 the cladding tube correction factor F can be obtained by the following formula (2):
- the relationship between the initial ovality and the correction factor is determined according to the relationship curve of the cladding tube correction factors of the initial ovality of different cladding tubes under the same cladding tube length;
- the relationship between the cladding tube internal and external pressure differential P and the correction factor is determined.
- the beneficial effects of the present invention By establishing a fuel rod cladding creep finite length tube correction method, the problem of the influence of the end boundary conditions on the cladding creep behavior in the fuel rod cladding creep collapse analysis is solved, and the fuel rod cladding is improved.
- the rationality and reliability of shell creep analysis ensure the reliability of fuel rod cladding creep collapse analysis in engineering analysis and design, which indirectly improves the safety and economy of nuclear power plants.
- FIG. 1 is a flowchart of a method for establishing a fuel rod cladding creep finite length tube correction method according to an embodiment of the present invention
- FIG. 2 is a schematic diagram of the setting process of the cladding tube geometric model in the present invention.
- a method for establishing a fuel rod cladding creep finite length tube correction method may include the following steps:
- the parameters of the cladding tube include the outer diameter, thickness, tube length (length of the cladding tube), initial ovality, material type of the cladding tube, pressure difference P between the cladding tube and the like.
- the cladding tubes of different lengths are determined according to the needs, and the destabilizing force of the cladding tubes of different lengths is analyzed by finite element software to determine the equivalent length of the cladding tube of infinite length, ie the maximum The length of the cladding tube.
- the length of the cladding tube is 10 mm, 20 mm, 40 mm, and 60 mm, and it is assumed that the tube length of 60 mm is infinite.
- the cladding tubes with different initial ovality are determined according to the needs, and the maximum initial ovality of the cladding tube is determined according to the processing and manufacturing of the fuel rod cladding tube.
- the initial ovality of the cladding tube is 0.02mm, 0.04mm, and 0.06mm, assuming that 0.06mm is the maximum initial ovality.
- the material type of the cladding tube needs to be determined, for example, Zr-4.
- the internal and external pressure difference P is used as the magnitude of the load on the cladding tube, for example, 9 MPa.
- the pressure difference P inside and outside the cladding tube is equal everywhere.
- the finite element software can be implemented using finite element software such as ABAQUS.
- This step S2 may further include:
- this step S2.1 may include S2.1.1, calculating the long axis and short axis of the cladding tube according to the outer diameter, thickness and initial ovality of the cladding tube, and according to the tube length of the cladding tube and the obtained length Axis and short axis, establish the geometric model of the cladding tube; S2.1.2. Divide the cladding tube into four parts, and set the midpoints of the four dividing lines as nodes, to obtain four nodes on the ellipse of the middle section of the cladding tube.
- the physical parameters of the cladding tube include Young's modulus and Poisson's ratio.
- the interface property parameters of the cladding tube include the cladding tube interface model and thickness.
- the cladding tube is also set to contact after the analysis step is set according to the need, and the non-slip, hard contact method is selected, as shown in Figure 2 (b) As shown.
- the setting of the boundary conditions includes constraining the two ends of the cladding tube and setting the outer surface of the cladding tube to be under pressure, as shown in (c) in FIG. 2.
- one end of the cladding tube constrains 6 degrees of freedom, and the other end constrains 5 degrees of freedom.
- the axial (z-direction) displacement of the cladding tube is not constrained.
- the pressure on the outer surface of the cladding tube is consistent with the pressure difference between the inside and outside of the cladding tube.
- the cladding tube ovality ty under different cladding tube lengths and different initial ovality can be calculated.
- the initial ovality of the cladding tube is 0.02mm, 0.04mm and 0.06mm respectively;
- the tube length is infinite tube length (ie 60mm)
- the initial ovality is 0.02mm and 0.04mm respectively
- the ovality of the cladding tube of 0.06mm is 0.02mm and 0.04mm respectively.
- stiffness of the cladding tube under different tube length and different initial ovality is obtained ; stiffness of the cladding tube comprises a finite time for the finite length of the cladding tube and the cladding tube stiffness K is infinite stiffness K infinity.
- the stiffness of the cladding tube is defined as the ratio of the pressure difference between the inside and outside of the cladding tube to the ellipticity, where the stiffness K is limited to the pressure difference P and the cladding of the cladding tube when the length is limited.
- the ratio of tube ovality ty, unit is MPa / mm; the stiffness K is infinite is the ratio of internal and external pressure difference P of cladding tube and cladding tube ovality ty at infinite length, unit is MPa / mm.
- the stiffness of the cladding tube is an inherent characteristic of the material or structure.
- the cladding tube rigidity is only related to the length of the cladding tube. Because the stiffness of the cladding tube reflects the ability of the cladding tube to resist deformation, the correction factor of the finite-length tube can be analyzed and calculated through the difference between the stiffness of the finite-length and infinite-length cladding tubes, and the creep model of the cladding can be corrected.
- the cladding tube correction factor F can be obtained by the following formula (2):
- the correction factor of the cladding tube under different tube lengths and initial ovality can be calculated.
- the initial ellipticity is 0.02mm, 0.04mm and 0.06mm cladding tube correction factor
- the tube length is infinite tube length (that is, 60mm)
- the initial ellipticity is 0.02mm, 0.04mm And 0.06mm cladding tube correction factor.
- the relationship between the correction factors of the cladding tubes with different initial ovality of the cladding tube under the same cladding tube length the relationship between the initial ovality and the correction factor can be determined: the initial ovality has very little effect on the correction factor.
- the relationship between the cladding tube correction factors of different cladding tube lengths under the same initial ovality of the cladding tube can be determined: as the tube length increases, the correction factor continues to increase, exceeding a certain No further changes after the value.
- the cladding tube correction factor is only related to the tube length, so the finite length tube can be determined
- the curve of the correction factor is used to establish a finite-length tube correction method for the correction of the cladding creep behavior in the analysis of fuel rod cladding creep collapse.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
Abstract
一种燃料棒包壳蠕变有限长管修正方法的建立方法，包括：S1、确定包壳管的参数；S2、分析处理获得在不同管长和不同初始椭圆度下包壳管的椭圆度；S3、根据获得的不同管长和不同初始椭圆度下包壳管的椭圆度、结合包壳管内外压差，获得不同管长和不同初始椭圆度下包壳管的刚度；S4、获得包壳管修正因子F；S5、根据包壳管管长、初始椭圆度和包壳管内外压差对包壳管修正因子的影响性分析，确定有限长管修正因子的曲线，建立有限长管修正方法。该方法解决了燃料棒包壳蠕变坍塌分析中端部边界条件对包壳蠕变行为影响的问题，提高了燃料棒包壳蠕变分析的合理性和可靠性，保证了燃料棒在工程分析和设计中包壳蠕变坍塌分析可靠性。
Description
本发明涉及一种反应堆中燃料棒包壳蠕变分析技术领域，尤其涉及一种燃料棒包壳蠕变有限长管修正方法的建立方法。
在上世纪70～80年代，燃料棒包壳蠕变坍塌是燃料棒失效的主要原因之一，因此燃料棒蠕变坍塌成为是燃料棒设计准则之一。
在压水堆运行过程中，陶瓷芯块在辐照初会发生密实现象，并导致芯块柱轴向高度变短，而燃料棒包壳由于辐照生长将会变长。由于冷却剂压力大于燃料棒内压以及包壳的蠕变行为，芯块与包壳接触后可能在芯块柱轴向出现没有芯块支撑的包壳管段，即未支撑包壳管段。在包壳外压和蠕变作用下，失去芯块支撑的有限长包壳管向内蠕变，并导致包壳管变形不断增加，包壳变形到达一定程度或包壳被压扁，即发生包壳蠕变坍塌。
由于燃料棒中未支撑包壳管段长度是有限的，两端受芯块的支撑作用，而包壳蠕变模型是在空包壳管上进行实验获得的，实验所用空包壳管相对于未支撑包壳管段是无限长的，因此不能直接通过包壳蠕变模型计算未支撑管段的蠕变行为，而需对包壳蠕变模型进行修正，其修正因子称为有限长管修正因子。
由于包壳蠕变有限长管修正因子直接影响包壳蠕变行为，而包壳蠕变是 燃料棒包壳蠕变坍塌分析中的核心，因此有限长管修正因子是包壳蠕变分析中的重要参数，直接影响包壳蠕变分析的合理性和可靠性。
发明内容
本发明要解决的技术问题在于，针对现有的上述缺陷，提供一种燃料棒包壳蠕变有限长管修正方法的建立方法。
本发明解决其技术问题所采用的技术方案是：提供一种燃料棒包壳蠕变有限长管修正方法的建立方法，包括以下步骤：
S1、确定包壳管的参数，所述参数包括包壳管的外直径、厚度、管长、初始椭圆度、包壳管材料类型以及包壳管内外压差P；
S2、通过有限元软件分析处理获得在不同管长和不同初始椭圆度下包壳管的椭圆度ty；
S3、根据获得的不同管长和不同初始椭圆度下包壳管的椭圆度ty、结合包壳管内外压差P，获得不同管长和不同初始椭圆度下包壳管的刚度；包壳管的刚度包括包壳管为有限长时的刚度K
有限以及包壳管为无限长时的刚度K
无限；
S4、根据下式(1)获得包壳管修正因子F：
S5、根据包壳管管长、初始椭圆度和包壳内外压差P对包壳管修正因子的影响性分析，确定有限长管修正因子的曲线，建立有限长管修正方法。
优选地，步骤S1中，根据包壳管材料类型确定包壳管的物性参数；所述包壳管的物性参数包括杨氏模量和泊松比。
优选地，步骤S1中，所述包壳管内外压差P各处相等。
优选地，步骤S2包括：
S2.1、建立包壳管几何模型，根据包壳管两端椭圆环形的长轴和短轴，将包壳管分为四部分，并将四条分割线的中点设置为节点，得到包壳管中间截面椭圆上的四个节点；
S2.2、设置包壳管的物性参数、包壳管界面属性参数；
S2.3、分析步设置：设置包壳管的分析模型为静力分析，设置输出四个节点在x、y和z三个方向的位移；
S2.4、边界条件设置：对包壳管的两端分别进行约束设置，设置包壳管外表面为均匀受压，压力与包壳管内外压差一致；
S2.5、网格划分：将包壳管几何模型进行网格划分；
S2.6、计算分析和后处理：通过求解器对包壳管几何模型进行分析计算，计算完成后，得到包壳管变形云图；将包壳管上四个节点位移导出，计算包壳管椭圆度ty。
优选地，步骤S2.1包括：
S2.1.1、根据包壳管的外直径、厚度和初始椭圆度计算获得包壳管的长轴和短轴，根据包壳管的管长以及获得的长轴和短轴，建立包壳管几何模型；
S2.1.2、将包壳管分为四部分，并将四条分割线的中点设置为节点，得到包壳管中间截面椭圆上的四个节点。
优选地，在分析步设置后还对包壳管进行接触设置，选用无滑移、硬接触的方式。
优选地，步骤S2.2中，所述包壳管的物性参数包括杨氏模量和泊松比；包壳管的界面属性参数包括包壳管界面模型、厚度。
优选地，步骤S2.6中，根据导出的包壳管上四个节点的位移，结合包壳 管变形前的长轴和短轴，获得包壳管变形后的长轴D
L和短轴D
S，根据ty＝D
L-D
S获得包壳管椭圆度ty。
优选地，步骤S3中，刚度K
有限为有限长时包壳管内外压差P和包壳管椭圆度ty的比值；
刚度K
无限为无限长时包壳管内外压差P和包壳管椭圆度ty的比值；
步骤S4中，包壳管修正因子F可由下式(2)获得：
优选地，结合步骤S3、S4，根据相同包壳管管长下不同包壳管初始椭圆度的包壳管修正因子关系曲线，确定初始椭圆度和修正因子的关系；
根据相同包壳管初始椭圆度下不同包壳管管长的包壳管修正因子关系曲线，确定包壳管管长和修正因子的关系；
根据相同包壳管管长和初始椭圆度下不同包壳内外压差P的包壳管修正因子关系曲线，确定包壳管内外压差P和修正因子的关系。
本发明的有益效果：通过建立燃料棒包壳蠕变有限长管修正方法，解决了燃料棒包壳蠕变坍塌分析中端部边界条件对包壳蠕变行为影响的问题，提高了燃料棒包壳蠕变分析的合理性和可靠性，保证了燃料棒在工程分析和设计中包壳蠕变坍塌分析可靠性，间接提高核电站的安全性和经济性。
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1是本发明一实施例的燃料棒包壳蠕变有限长管修正方法的建立方法流程图；
图2是本发明中包壳管几何模型的设置过程示意图。
结合图1、图2，本发明一实施例的燃料棒包壳蠕变有限长管修正方法的建立方法，可包括以下步骤：
S1、确定包壳管的参数。
包壳管的参数包括包壳管的外直径、厚度、管长(包壳管长度)、初始椭圆度、包壳管材料类型、包壳管内外压差P等。
其中，为了分析管长对修正因子的影响，根据需要确定不同长度的包壳管，并通过有限元软件分析不同长度包壳管的失稳压力，确定无限长包壳管等效长度，即最大的包壳管长度。例如包壳管长度为10mm、20mm、40mm、60mm，假设60mm管长为无限管长。
为了分析初始椭圆度对修正因子的影响，根据需要确定不同初始椭圆度的包壳管，根据燃料棒包壳管加工制造等情况，确定包壳管最大初始椭圆度。例如包壳管初始椭圆度为0.02mm、0.04mm和0.06mm，假设0.06mm为最大初始椭圆度。
由于不同的包壳管材料的物性参数(杨氏模量、泊松比)不同，因此需要确定包壳管的材料类型，例如Zr-4。
确定包壳管的内外压差P，该内外压差P作为包壳管所受载荷大小，例如9MPa。包壳管内外压差P各处相等。
S2、通过有限元软件分析处理获得在不同管长和不同初始椭圆度下包壳管的椭圆度ty。有限元软件可以采用如ABAQUS等有限元软件实现。
该步骤S2进一步可包括：
S2.1、建立包壳管几何模型，根据包壳管两端椭圆环形的长轴和短轴，将包壳管分为四部分，并将四条分割线的中点设置为节点，得到包壳管中间截面椭圆上的四个节点，如图2中(a)所示。
具体地，该步骤S2.1可包括S2.1.1、根据包壳管的外直径、厚度和初始椭圆度计算获得包壳管的长轴和短轴，根据包壳管的管长以及获得的长轴和短轴，建立包壳管几何模型；S2.1.2、将包壳管分为四部分，并将四条分割线的中点设置为节点，得到包壳管中间截面椭圆上的四个节点。
S2.2、设置包壳管的物性参数、包壳管界面属性参数，该设置的各参数赋予上述建立的包壳管几何模型中。
包壳管的物性参数包括杨氏模量和泊松比等。包壳管的界面属性参数包括包壳管界面模型、厚度等。
S2.3、分析步设置：设置包壳管的分析模型为静力分析，根据需要设置分析时间和时间步长，将包壳管中间截面椭圆上的四个节点创建为集合，设置输出四个节点在x、y和z三个方向的位移。
考虑到包壳管在受压过程中，其内表面可能接触，根据需要在分析步设置后还对包壳管进行接触设置，选用无滑移、硬接触的方式，如图2中(b)所示。
S2.4、边界条件设置。
边界条件设置包括对包壳管的两端分别进行约束设置、设置包壳管外表面为均有受压，如图2中(c)所示。
其中，包壳管一端约束6个方向的自由度，另一端约束5个方向自由度，包壳管的轴向(z方向)位移不约束。包壳管外表面受压的压力与包壳管内外压差一致。
S2.5、网格划分：将包壳管几何模型进行网格划分。
将网格单元设置为双曲薄/厚壳、减缩积分、沙漏控制及有限膜应变单元，并赋予给包壳管几何模型，对包壳管几何模型进行网格划分，划分为四边形网格，如图2中(d)所示。
S2.6、计算分析和后处理：通过求解器对包壳管几何模型进行分析计算，计算完成后，得到包壳管变形云图，从包壳管变形云图中可观察包壳管的变形情况，包括包壳管上网格的位移。将包壳管上四个节点位移导出，计算包壳管椭圆度ty。
从包壳管变形云图显示可知，包壳管自中间位置到两端位置，变形(包壳管上各个节点、网格的位移)逐渐减小。
根据导出的包壳管上四个节点的位移，结合包壳管变形前的长轴和短轴，获得包壳管变形后的长轴D
L和短轴D
S，根据ty＝D
L-D
S获得包壳管椭圆度ty。
根据上述包壳管椭圆度ty的计算方法，可以计算不同包壳管管长和不同初始椭圆度下包壳管的椭圆度ty。例如，管长为10mm时，初始椭圆度分别为0.02mm、0.04mm和0.06mm的包壳管椭圆度；管长为无限管长(即60mm)时，初始椭圆度分别为0.02mm、0.04mm和0.06mm的包壳管椭圆度。
S3、根据获得的不同管长和不同初始椭圆度下包壳管的椭圆度ty、结合包壳管内外压差P(单位MPa)，获得不同管长和不同初始椭圆度下包壳管的刚度；包壳管的刚度包括包壳管为有限长时的刚度K
有限以及包壳管为无限长时的刚度K
无限。
为了准确反映包壳管变形与刚度的关系，将包壳管刚度定义为包壳管内外压差与椭圆度的比值，其中，刚度K
有限为有限长时包壳管内外压差P和包 壳管椭圆度ty的比值，单位为MPa/mm；刚度K
无限为无限长时包壳管内外压差P和包壳管椭圆度ty的比值，单位为MPa/mm。
包壳管刚度是材料或结构的固有特性。对于受压包壳管，当外部载荷、材料性质、几何形状和边界条件等工况条件相同，包壳管刚度只与包壳管长有关。由于包壳管刚度反应了包壳管抵抗变形的能力，因此可以通过有限长和无限长包壳管刚度的不同，分析计算有限长管修正因子，进而实现对包壳蠕变模型的修正。
S4、根据下式(1)获得包壳管修正因子F：
进一步地，包壳管修正因子F可由下式(2)获得：
根据上式(1)、(2)即可计算获得不同管长和初始椭圆度下包壳管修正因子。例如，管长为10mm时，初始椭圆度分别为0.02mm、0.04mm和0.06mm的包壳管修正因子；管长为无限管长(即60mm)时，初始椭圆度分别为0.02mm、0.04mm和0.06mm的包壳管修正因子。
结合上述步骤S3、S4，即可：
根据相同包壳管管长下不同包壳管初始椭圆度的包壳管修正因子关系曲线，可以确定初始椭圆度和修正因子的关系：初始椭圆度对修正因子影响非常小。
根据相同包壳管初始椭圆度下不同包壳管管长的包壳管修正因子关系曲线，可以确定包壳管管长和修正因子的关系：随着管长增加，修正因子不断 增加，超过一定值后不再变化。
根据相同包壳管管长和初始椭圆度下不同包壳内外压差P的包壳管修正因子关系曲线，确定包壳管内外压差P和修正因子的关系：包壳管内外压差对修正因子影响非常小。
S5、根据包壳管管长、初始椭圆度和包壳内外压差P对包壳管修正因子的影响性分析，可获得包壳管修正因子只与其中管长相关，因此可确定有限长管修正因子的曲线，建立有限长管修正方法，以用于燃料棒包壳蠕变坍塌分析中包壳蠕变行为的修正。
以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。
Claims (10)
- 一种燃料棒包壳蠕变有限长管修正方法的建立方法，其特征在于，包括以下步骤：S1、确定包壳管的参数，所述参数包括包壳管的外直径、厚度、管长、初始椭圆度、包壳管材料类型以及包壳管内外压差P；S2、有限元分析处理：通过有限元软件分析处理获得在不同管长和不同初始椭圆度下包壳管的椭圆度ty；S3、获得包壳管刚度：根据获得的不同管长和不同初始椭圆度下包壳管的椭圆度ty、结合包壳管内外压差P，获得不同管长和不同初始椭圆度下包壳管的刚度；包壳管的刚度包括包壳管为有限长时的刚度K 有限以及包壳管为无限长时的刚度K 无限；S4、获得包壳管修正因子：根据下式(1)获得包壳管修正因子F；
S5、建立有限长管修正方法：根据包壳管管长、初始椭圆度和包壳内外压差P对包壳管修正因子的影响性分析，确定有限长管修正因子的曲线，建立有限长管修正方法。 - 根据权利要求1所述的燃料棒包壳蠕变有限长管修正方法，其特征在于，步骤S1中，根据包壳管材料类型确定包壳管的物性参数；所述包壳管的物性参数包括杨氏模量和泊松比。
- 根据权利要求1所述的燃料棒包壳蠕变有限长管修正方法，其特征在于，步骤S1中，所述包壳管内外压差P各处相等。
- 根据权利要求1所述的燃料棒包壳蠕变有限长管修正方法，其特征在于，步骤S2包括：S2.1、建立包壳管几何模型，根据包壳管两端椭圆环形的长轴和短轴，将包壳管分为四部分，并将四条分割线的中点设置为节点，得到包壳管中间截面椭圆上的四个节点；S2.2、设置包壳管的物性参数、包壳管界面属性参数；S2.3、分析步设置：设置包壳管的分析模型为静力分析，设置输出四个节点在x、y和z三个方向的位移；S2.4、边界条件设置：对包壳管的两端分别进行约束设置，设置包壳管外表面为均匀受压，压力与包壳管内外压差一致；S2.5、网格划分：将包壳管几何模型进行网格划分；S2.6、计算分析和后处理：通过求解器对包壳管几何模型进行分析计算，计算完成后，得到包壳管变形云图；将包壳管上四个节点位移导出，计算包壳管椭圆度ty。
- 根据权利要求4所述的燃料棒包壳蠕变有限长管修正方法，其特征在于，步骤S2.1包括：S2.1.1、根据包壳管的外直径、厚度和初始椭圆度计算获得包壳管的长轴和短轴，根据包壳管的管长以及获得的长轴和短轴，建立包壳管几何模型；S2.1.2、将包壳管分为四部分，并将四条分割线的中点设置为节点，得到包壳管中间截面椭圆上的四个节点。
- 根据权利要求4所述的燃料棒包壳蠕变有限长管修正方法，其特征在于，在分析步设置后还对包壳管进行接触设置，选用无滑移、硬接触的方式。
- 根据权利要求4所述的燃料棒包壳蠕变有限长管修正方法，其特征在 于，步骤S2.2中，所述包壳管的物性参数包括杨氏模量和泊松比；包壳管的界面属性参数包括包壳管界面模型、厚度。
- 根据权利要求4所述的燃料棒包壳蠕变有限长管修正方法，其特征在于，步骤S2.6中，根据导出的包壳管上四个节点的位移，结合包壳管变形前的长轴和短轴，获得包壳管变形后的长轴D L和短轴D S，根据ty＝D L-D S获得包壳管椭圆度ty。
- 根据权利要求1所述的燃料棒包壳蠕变有限长管修正方法，其特征在于，步骤S3中，刚度K 有限为有限长时包壳管内外压差P和包壳管椭圆度ty的比值；刚度K 无限为无限长时包壳管内外压差P和包壳管椭圆度ty的比值；步骤S4中，包壳管修正因子F可由下式(2)获得：
- 根据权利要求1-9任一项所述的燃料棒包壳蠕变有限长管修正方法，其特征在于，结合步骤S3、S4，根据相同包壳管管长下不同包壳管初始椭圆度的包壳管修正因子关系曲线，确定初始椭圆度和修正因子的关系；根据相同包壳管初始椭圆度下不同包壳管管长的包壳管修正因子关系曲线，确定包壳管管长和修正因子的关系；根据相同包壳管管长和初始椭圆度下不同包壳内外压差P的包壳管修正因子关系曲线，确定包壳管内外压差P和修正因子的关系。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811269670.8 | 2018-10-29 | ||
| CN201811269670.8A CN109583022B (zh) | 2018-10-29 | 2018-10-29 | 燃料棒包壳蠕变有限长管修正方法的建立方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2020088375A1 true WO2020088375A1 (zh) | 2020-05-07 |
Family
ID=65920700
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/CN2019/113442 WO2020088375A1 (zh) | 2018-10-29 | 2019-10-25 | 燃料棒包壳蠕变有限长管修正方法的建立方法 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN109583022B (zh) |
| WO (1) | WO2020088375A1 (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109583022B (zh) * | 2018-10-29 | 2020-06-23 | 中广核研究院有限公司 | 燃料棒包壳蠕变有限长管修正方法的建立方法 |
| CN114676523A (zh) * | 2022-03-31 | 2022-06-28 | 中广核研究院有限公司 | 有限长包壳管蠕变变形方法、装置和计算机设备 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101243521A (zh) * | 2005-08-10 | 2008-08-13 | 法国原子能委员会 | 大结构板式燃料元件 |
| US20120321031A1 (en) * | 2009-12-04 | 2012-12-20 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Nuclear fuel rod and method of manufacturing pellets for such a rod |
| CN104931348A (zh) * | 2015-06-08 | 2015-09-23 | 西南交通大学 | 圆环径向压缩能量预测材料单轴本构关系测定方法 |
| CN105934797A (zh) * | 2014-01-27 | 2016-09-07 | 泰拉能源公司 | 用于燃料元件变形的建模 |
| CN109583022A (zh) * | 2018-10-29 | 2019-04-05 | 中广核研究院有限公司 | 燃料棒包壳蠕变有限长管修正方法的建立方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108446505B (zh) * | 2018-03-23 | 2021-06-15 | 东北大学 | 一种漏斗结晶器内铸坯凝固传热计算方法 |
-
2018
- 2018-10-29 CN CN201811269670.8A patent/CN109583022B/zh active Active
-
2019
- 2019-10-25 WO PCT/CN2019/113442 patent/WO2020088375A1/zh active Application Filing
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101243521A (zh) * | 2005-08-10 | 2008-08-13 | 法国原子能委员会 | 大结构板式燃料元件 |
| US20120321031A1 (en) * | 2009-12-04 | 2012-12-20 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Nuclear fuel rod and method of manufacturing pellets for such a rod |
| CN105934797A (zh) * | 2014-01-27 | 2016-09-07 | 泰拉能源公司 | 用于燃料元件变形的建模 |
| CN104931348A (zh) * | 2015-06-08 | 2015-09-23 | 西南交通大学 | 圆环径向压缩能量预测材料单轴本构关系测定方法 |
| CN109583022A (zh) * | 2018-10-29 | 2019-04-05 | 中广核研究院有限公司 | 燃料棒包壳蠕变有限长管修正方法的建立方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109583022A (zh) | 2019-04-05 |
| CN109583022B (zh) | 2020-06-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2020088375A1 (zh) | 燃料棒包壳蠕变有限长管修正方法的建立方法 | |
| WO2021139287A1 (zh) | 负泊松比结构的非线性本构关系分析方法、系统及装置 | |
| Murtaza et al. | The effects of thermal stresses on the elliptical surface cracks in PWR reactor pressure vessel | |
| Li et al. | Buckling behavior of large-scale thin-walled ellipsoidal head under internal pressure | |
| WO2023184883A1 (zh) | 有限长包壳管蠕变变形计算方法、装置和计算机设备 | |
| CN105260560A (zh) | 一种带沟槽的复合材料管的优化设计方法 | |
| CN106777780B (zh) | 管子-管板接头柔性静压胀接机理分析方法 | |
| Cheney | Pressure buckling of ring encased in cavity | |
| US20220245313A1 (en) | Method for calculating fluid-structure interaction response of ceramic matrix composites | |
| Rosenzweig et al. | Capillary-elastic instabilities of liquid-lined lung airways | |
| Wu et al. | Progressive Collapse Resistance of Formwork Support System with Couplers. | |
| Valdeolivas et al. | A full 3D finite element model for buckling analysis of stiffened steel liners in hydroelectric pressure tunnels | |
| Baas et al. | Influence of Nozzle Transition on Leak Before Break Parameters | |
| Valdeolivas et al. | Consideration of geometric imperfections in three-dimensional finite element model analysis of stiffened steel liners subjected to external pressure | |
| CN111883216B (zh) | 一种基于材料属性优化的仿真网格划分方法及系统 | |
| CN203867066U (zh) | 在钢管和混凝土中间加设记忆合金层的自填充钢管混凝土 | |
| CN115061234B (zh) | 一种声压高灵敏度实芯光子晶体光纤、制备方法及水声器 | |
| CN113887014A (zh) | 临界热流密度关系式的确定方法、可读存储介质及开发装置 | |
| del Coz Díaz et al. | Non-linear Analysis of Sheet Cover (Umbrella) of Reinforced Concrete of 40 m diameter | |
| Gu et al. | Stability analysis and test research of single-layer latticed dome considering bar buckling | |
| Sumesh et al. | Structural Deformations on Critical Cracked Pressurised Pipe Bends: Implication on the Von Mises Stresses | |
| Nomura et al. | Design condition of a sustainable two-layer circular tube with energy absorbing capacity | |
| Tallavo et al. | A Comparative Evaluation of Finite Element Modeling of Creep Deformation of Fuel Channels in CANDU® Nuclear Reactors | |
| CN103117521B (zh) | 电缆附件扩张形变与结构恢复方法 | |
| Yu | Three-dimensional modeling of Candu nuclear fuel elements in bowing analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 19878730 Country of ref document: EP Kind code of ref document: A1 |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 19878730 Country of ref document: EP Kind code of ref document: A1 |
|
| 32PN | Ep: public notification in the ep bulletin as address of the adressee cannot be established |
Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 09.09.2021) |
|
| 122 | Ep: pct application non-entry in european phase |
Ref document number: 19878730 Country of ref document: EP Kind code of ref document: A1 |