Heavy water reactor pressure tube irradiation deformation modeling method, system, apparatus, and storage medium

Abstract

The application provides a heavy water reactor pressure tube irradiation deformation modeling method, system, device and storage medium, the method comprises the following steps: obtaining and saving input data through a script, the input data comprising the geometric parameters of the pressure tube to be calculated, irradiation parameters and multiple sets of working condition parameters, the multiple sets of working condition parameters comprising the fluence distribution, temperature distribution and pressure distribution of multiple measuring points of the pressure tube to be calculated; performing a viscoplastic self-consistent algorithm through a script, and saving the output constitutive model parameters, the constitutive model parameters comprising thermal creep anisotropy coefficients, irradiation creep anisotropy coefficients and irradiation growth strain coefficients; arranging and saving the input data and the constitutive model parameters into an input file through a script; calling a first interface through a script, obtaining and saving a finite element geometric model generated according to the input file and a preconfigured file; and calling a second interface through a script, obtaining and saving multiple sets of deformation data corresponding to the multiple measuring points respectively calculated according to the finite element geometric model.

Description

Technical Field

[0001] This application mainly relates to the field of pressure tube irradiation deformation analysis technology, and in particular to a method, system, equipment and storage medium for modeling pressure tube irradiation deformation of heavy water reactors. Background Technology

[0002] The pressure tube is a core component of a heavy water reactor, serving as a fuel container, cooling channel, pressure boundary, and safety barrier. Specifically, the pressure tube provides structural support for the fuel rod bundles and precisely positions them within the reactor core, while simultaneously securing them to the piping via positioning components. High-temperature, high-pressure heavy water coolant flows through the pressure tube, carrying the heat generated by nuclear fission out of the core and to the steam generator for power generation. As a critical component of the reactor's primary loop pressure boundary, the pressure tube maintains its structural integrity under high temperature, high pressure, and high neutron irradiation conditions, acting as the first physical barrier to prevent radioactive material leakage. Furthermore, the pressure tube is independent of the fuel rod bundles and can slide axially to enable refueling operations without reactor shutdown. The pressure tube also isolates the high-temperature, high-pressure coolant inside from the low-temperature, low-pressure moderator outside, preventing mutual interference.

[0003] Every component of the pressure tube in a heavy water reactor is subject to the high temperature, high pressure, and high-neutron radiation environment within the reactor. Deformation of the pressure tube can lead to various problems, such as fuel rod bypass, pressure tube thinning, and contact between the pressure tube and the runoff pipes, thus limiting the normal operation of the heavy water reactor. Therefore, the irradiation deformation of the pressure tube must be carefully considered at every stage of heavy water reactor design, manufacturing, deployment, and operation to ensure its continued effectiveness and extend the lifespan of the heavy water reactor unit. However, current methods typically employ manual modeling and calculation of deformation results. This approach still requires significant manual intervention, resulting in a large workload, low reuse rate, and high cost. Summary of the Invention

[0004] The technical problem to be solved by this application is to provide a method, system, device and storage medium for modeling the irradiation deformation of heavy water reactor pressure tubes, which can realize automated modeling of the irradiation deformation of heavy water reactor pressure tubes, as well as batch calculation and visualization output of deformation data.

[0005] To address the aforementioned technical problems, this application provides a method for modeling the irradiation deformation of a pressure tube in a heavy water reactor, comprising: acquiring and saving input data via a script, the input data including geometric parameters, irradiation parameters, and multiple sets of operating condition parameters of the pressure tube to be calculated, wherein the multiple sets of operating condition parameters include flux distribution, temperature distribution, and pressure distribution at multiple measurement points of the pressure tube to be calculated; executing a viscoplastic self-consistent algorithm via the script and saving the constitutive model parameters output by the viscoplastic self-consistent algorithm, the constitutive model parameters including thermal creep anisotropy coefficient, irradiation creep anisotropy coefficient, and irradiation growth strain coefficient; organizing and saving the input data and constitutive model parameters into an input file via the script; calling a first interface via the script to acquire and save a finite element geometric model generated based on the input file and a pre-configuration file; and calling a second interface via the script to acquire and save multiple sets of deformation data calculated based on the finite element geometric model corresponding to the multiple measurement points respectively.

[0006] In one embodiment of this application, the heavy water reactor pressure tube irradiation deformation modeling method further includes constructing a visualization interface through the script. The step of constructing the visualization interface through the script includes: establishing a first coordinate system and a second coordinate system through the script, wherein the horizontal axis of the first coordinate system and the second coordinate system represents the plurality of measurement points, the vertical axis of the first coordinate system represents the tube wall thinning value, and the vertical axis of the second coordinate system represents the inner diameter increase value; plotting the coordinate points of the tube wall thinning value in the plurality of deformation data corresponding to each measurement point in the first coordinate system through the script; and plotting the coordinate points of the inner diameter increase value in the plurality of deformation data corresponding to each measurement point in the second coordinate system through the script.

[0007] In one embodiment of this application, the step of constructing a visualization interface through the script further includes: obtaining user interaction actions through the script; when the user interaction action includes selecting a coordinate point, displaying the pipe wall thinning value and inner diameter growth value of the measurement point corresponding to the selected coordinate point in the first coordinate system and the second coordinate system, respectively, through the script.

[0008] In one embodiment of this application, the geometric parameters of the pressure tube to be calculated include the tube length, inner diameter, and outer diameter, and the irradiation parameters of the pressure tube to be calculated include the continuous irradiation time and the reconstructed irradiation time.

[0009] In one embodiment of this application, the first interface includes a command-line interface of finite element analysis software. The step of generating and saving a finite element geometric model based on the input file and the pre-configuration file by calling the first interface through the script includes: submitting the input file and the pre-configuration file to the command-line interface through the script; and obtaining and saving the finite element geometric model returned by the command-line interface through the script.

[0010] In one embodiment of this application, the second interface includes a cross-sectional strain calculation interface, a cross-sectional stress calculation interface, an axial growth calculation interface, and an inner and outer diameter dimension change calculation interface of the finite element analysis software. The step of calling the second interface through the script to calculate and save multiple sets of deformation data corresponding to the multiple measurement points based on the finite element geometric model includes: submitting the finite element geometric model to the cross-sectional strain calculation interface through the script and obtaining the output cross-sectional strain corresponding to the multiple measurement points; submitting the finite element geometric model to the cross-sectional stress calculation interface through the script and obtaining the output cross-sectional stress corresponding to the multiple measurement points; submitting the finite element geometric model to the axial growth calculation interface through the script and obtaining the output axial growth value corresponding to the multiple measurement points; submitting the finite element geometric model to the inner and outer diameter dimension change calculation interface through the script and obtaining the output inner and outer diameter dimension changes corresponding to the multiple measurement points; and saving the cross-sectional strain, cross-sectional stress, axial growth value, and inner and outer diameter dimension changes corresponding to the multiple measurement points through the script.

[0011] In one embodiment of this application, the viscoplastic self-consistent algorithm includes a first algorithm and a second algorithm. The first algorithm is used to calculate the thermal creep anisotropy coefficient, and the second algorithm is used to calculate the irradiation creep anisotropy coefficient and the irradiation growth strain coefficient. The step of executing the viscoplastic self-consistent algorithm through the script and saving the constitutive model parameters output by the viscoplastic self-consistent algorithm includes: obtaining the input of the first algorithm through the script, the input of the first algorithm including a first parameter configuration file, a pressure tube texture file, a first single crystal constitutive file, and a first loading parameter file; obtaining the input of the second algorithm through the script, the input of the second algorithm including a second parameter configuration file, the pressure tube texture file, a second single crystal constitutive file, and a second loading parameter file; executing the first algorithm and the second algorithm through the script respectively, and saving the thermal creep anisotropy coefficient output by the first algorithm, and the irradiation creep anisotropy coefficient and the irradiation growth strain coefficient output by the second algorithm.

[0012] In one embodiment of this application, the step of organizing and saving the input data and constitutive model parameters into an input file using the script includes: obtaining, through the script, the thermal creep anisotropy coefficient, irradiation creep anisotropy coefficient, and irradiation growth strain coefficient selected by the user from the constitutive model parameters, respectively corresponding to the tube head and tube tail of the pressure tube to be calculated; and generating and saving the input file based on the input data, the thermal creep anisotropy coefficient, irradiation creep anisotropy coefficient, and irradiation growth strain coefficient corresponding to the tube head and tube tail of the pressure tube to be calculated, respectively.

[0013] This application also proposes a heavy water reactor pressure tube irradiation deformation modeling system, comprising: a data input module configured to acquire and save input data, the input data including geometric parameters, irradiation parameters, and multiple sets of operating condition parameters of the pressure tube to be calculated, wherein the multiple sets of operating condition parameters include flux distribution, temperature distribution, and pressure distribution at multiple measurement points; a parameter calculation module configured to execute a viscoplastic self-consistent algorithm and save constitutive model parameters output by the viscoplastic self-consistent algorithm, the constitutive model parameters including thermal creep anisotropy coefficient, irradiation creep anisotropy coefficient, and irradiation growth strain coefficient; a file generation module configured to organize and save an input file based on user-selected saved constitutive model parameters and input data; and a deformation calculation module configured to call a first interface to generate and save a finite element geometric model based on the input file and pre-configuration file, and to call a second interface to calculate and save multiple sets of deformation data corresponding to the multiple measurement points based on the finite element geometric model.

[0014] In one embodiment of this application, the heavy water reactor pressure tube irradiation deformation modeling system further includes a visualization module. The visualization module is configured to: establish a first coordinate system and a second coordinate system using a visualization drawing library, wherein the horizontal axis of the first coordinate system and the second coordinate system represents the plurality of measurement points, the vertical axis of the first coordinate system represents the tube wall thinning value, and the vertical axis of the second coordinate system represents the inner diameter increase value; plot the coordinate points of the tube wall thinning value in the plurality of deformation data corresponding to each measurement point in the first coordinate system; and plot the coordinate points of the inner diameter increase value in the plurality of deformation data corresponding to each measurement point in the second coordinate system.

[0015] In one embodiment of this application, the visualization module is further configured to: acquire user interaction actions; when the user interaction action includes selecting a coordinate point, display the pipe wall thinning value and inner diameter growth value of the measurement point corresponding to the selected coordinate point in the first coordinate system and the second coordinate system, respectively.

[0016] This application also proposes an electronic device comprising: a memory for storing instructions executable by a processor; and a processor for executing the instructions to implement the method described above.

[0017] This application also proposes a computer storage medium storing computer program code that, when executed by a processor, implements the method described above.

[0018] Compared with the prior art, this application has the following advantages: (1) The entire process of obtaining input data through script, executing viscoplastic self-consistent algorithm and automatically organizing it into input file, as well as automatically calling the first and second interfaces to calculate multiple sets of deformation data, realizes automated modeling and calculation from input data to deformation results, reduces the workload of manual modeling, manual modification of input files, repeated submission of finite element analysis jobs and manual extraction of deformation results, thereby improving the automation level of the heavy water reactor pressure tube irradiation deformation analysis process; (2) The modeling method of calling the first interface to generate the finite element geometric model and calling the second interface to calculate the deformation data is adopted, which facilitates the replacement of the finite element analysis program used for modeling, and improves the reusability of script code and the portability of script. Attached Figure Description

[0019] The accompanying drawings are included to provide a further understanding of this application; they are incorporated into and constitute a part of this application. The drawings illustrate embodiments of this application and, together with this specification, serve to explain the principles of this application. In the drawings: Figure 1 This is a flowchart illustrating a heavy water reactor pressure tube irradiation deformation modeling method according to an embodiment of this application; Figure 2 This is a flowchart illustrating the first algorithm of an embodiment of this application; Figure 3 This is a flowchart illustrating the second algorithm according to an embodiment of this application; Figure 4 This is a schematic flowchart illustrating the execution of a viscoplastic self-consistent algorithm according to an embodiment of this application; Figure 5 yes Figure 1 A flowchart illustrating the sub-steps included in step S130; Figure 6 yes Figure 1 A flowchart illustrating the sub-steps included in step S140; Figure 7 yes Figure 1 A flowchart illustrating the sub-steps included in step S150; Figure 8 This is a flowchart illustrating a method for modeling the deformation of a pressure tube in a heavy water reactor under irradiation, according to another embodiment of this application. Figure 9 This is a schematic diagram illustrating the process of constructing a visual interface according to an embodiment of this application; Figure 10 This is a flowchart illustrating the construction process of a visual interface according to another embodiment of this application; Figure 11 This is a schematic diagram of the structure of a heavy water reactor pressure tube irradiation deformation modeling system according to an embodiment of this application; Figure 12 This is a schematic diagram of the first page corresponding to the data input module of one embodiment of this application; Figure 13 This is a schematic diagram of the second page corresponding to the parameter calculation module of one embodiment of this application; Figure 14 This is a schematic diagram of the third page corresponding to the file generation module of one embodiment of this application; Figure 15 This is a schematic diagram of the fourth page corresponding to the deformation calculation module of one embodiment of this application; Figure 16 This is a schematic diagram of the structure of an electronic device according to an embodiment of this application.

[0020] Reference numerals: Modeling system 400, Data input module 411, Parameter calculation module 412, File generation module 413, Deformation calculation module 414, Visualization module 415, First page 421, Second page 422, Third page 423, Fourth page 424, Electronic device 500, Internal communication bus 501, Processor 502, Read-only memory 503, Random access memory 504, Communication port 505, Hard disk 506. Detailed Implementation

[0021] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are merely some examples or embodiments of this application. For those skilled in the art, these drawings can be applied to other similar scenarios without creative effort. Unless obvious from the context or otherwise specified, the same reference numerals in the drawings represent the same structures or operations.

[0022] As indicated in this application and claims, unless the context clearly indicates otherwise, the words "a," "an," "an," and / or "the" are not specifically singular and may include plural forms. Generally speaking, the terms "comprising" and "including" only indicate the inclusion of explicitly identified steps and elements, which do not constitute an exclusive list, and the method or apparatus may also include other steps or elements.

[0023] Unless otherwise specifically stated, the relative arrangement, numerical expressions, and values ​​of the components and steps described in these embodiments do not limit the scope of this application. It should also be understood that, for ease of description, the dimensions of the various parts shown in the drawings are not drawn to actual scale. Techniques, methods, and devices known to those skilled in the art may not be discussed in detail, but where appropriate, such techniques, methods, and devices should be considered part of the specification. In all examples shown and discussed herein, any specific values ​​should be interpreted as merely exemplary and not as limitations. Therefore, other examples of exemplary embodiments may have different values. It should be noted that similar reference numerals and letters in the following drawings denote similar items; therefore, once an item is defined in one drawing, it need not be further discussed in subsequent drawings.

[0024] Furthermore, it should be noted that the use of terms such as "first" and "second" to define components is merely for the purpose of distinguishing the corresponding components. Unless otherwise stated, these terms have no special meaning and therefore should not be construed as limiting the scope of protection of this application. In addition, although the terminology used in this application is selected from commonly known and used terms, some terms mentioned in this application's specification may have been chosen by the applicant according to his or her judgment, and their detailed meanings are explained in the relevant sections of this description. Moreover, this application should be understood not only through the actual terms used, but also through the meaning implied by each term.

[0025] It should be understood that when a component is referred to as "on another component," "connected to another component," "coupled to another component," or "in contact with another component," it can be directly on, connected to, coupled to, or in contact with that other component, or there may be an intervening component. In contrast, when a component is referred to as "directly on another component," "directly connected to," "directly coupled to," or "directly in contact with" another component, there is no intervening component. Similarly, when a first component is referred to as "electrically contacting" or "electrically coupled to" a second component, there is an electrical path between the first and second components that allows current to flow. This electrical path may include capacitors, coupled inductors, and / or other components that allow current to flow, even if there is no direct contact between the conductive components.

[0026] Flowcharts are used in this application to illustrate the operations performed by the system according to embodiments of this application. It should be understood that the preceding or following operations are not necessarily performed in exact order. Instead, various steps can be processed in reverse order or simultaneously. Furthermore, other operations may be added to these processes, or one or more steps may be removed from these processes.

[0027] The following specific embodiments illustrate the irradiation deformation modeling method for the pressure tube of a heavy water reactor in this application.

[0028] refer to Figure 1 The flowchart shown is a schematic diagram of a heavy water reactor pressure tube irradiation deformation modeling method according to an embodiment. The heavy water reactor pressure tube irradiation deformation modeling method includes... Figure 1 Steps S110 to S150 will be explained in detail below.

[0029] In step S110, input data is acquired and saved through a script. The input data includes the geometric parameters, irradiation parameters, and multiple sets of operating condition parameters of the pressure tube to be calculated. The multiple sets of operating condition parameters include the injection distribution, temperature distribution, and pressure distribution of each of the multiple measurement points of the pressure tube to be calculated.

[0030] In this application, a script can be understood as a collection of programs consisting of one or more script files, executed by a runtime environment. The methods or steps of the embodiments of this application can be executed independently based on the collection of programs stored in a single script file, or the methods or steps of the embodiments of this application can be executed in combination based on the collection of programs stored in multiple script files. The multiple script files can call each other and transfer their internally stored data and functions. The script files can be stored in a computer-readable storage medium, and are invoked and run by the processor during execution. For example, script files include Python code files, C language files, and Java language files.

[0031] For example, in step S110, the script can be deployed on a terminal device to directly acquire user input data collected by the terminal device and save the input data in a readable storage medium on the terminal device. The script can also be deployed in the cloud to acquire user input data uploaded to the cloud and save the input data in a readable storage medium on the cloud. For example, the script can save the input data in an input.txt file and store the file in a readable storage medium on both the terminal device and the cloud.

[0032] The pressure tube of a heavy water reactor contains multiple independent fuel rod bundles arranged in series along the axial direction. Due to differences in operating parameters (including flux, temperature, and pressure) at different axial positions, the degree of irradiation deformation may vary across different parts of the pressure tube. Therefore, multiple measurement points are set along the axial direction of the pressure tube, and operating parameters are collected at each point to obtain local operating conditions at different axial positions. These results are then used to calculate and analyze deformation data at different axial positions in subsequent steps. Taking the CANDU-6 heavy water reactor as an example, a single pressure tube contains 12 fuel rod bundles, arranged in series along the axial direction. For instance, for a pressure tube like the CANDU-6 heavy water reactor to be calculated, 12 measurement points can be set at equal intervals along the axial direction of the pressure tube, and the flux distribution, temperature distribution, and pressure distribution at these 12 measurement points can be collected.

[0033] The geometric parameters of the pressure tube to be calculated represent its structural dimensions, providing a geometric basis for the subsequent generation of the finite element geometric model. The irradiation parameters of the pressure tube to be calculated represent its irradiation history (including irradiation time). These irradiation parameters can be set according to actual analysis needs, providing irradiation parameter input for the subsequent generation of the finite element geometric model.

[0034] In some embodiments, the geometric parameters of the pressure tube to be calculated include pipe length, inner diameter, and outer diameter, and the irradiation parameters include continuous irradiation time and reconfiguration irradiation time. The pipe length refers to the axial extension length of the pressure tube to be calculated; the inner diameter is the inner wall diameter; and the outer diameter is the outer wall diameter. The inner and outer diameters are used to determine the wall thickness of the pressure tube to be calculated. Continuous irradiation time refers to the cumulative irradiation duration of the pressure tube to be calculated since it was put into operation. Reconfiguration irradiation time refers to the cumulative irradiation duration corresponding to the change in irradiation growth direction from towards the tail to towards the head when the fixed end of the pressure tube to be calculated is switched from the head to the tail. For example, at the start of irradiation, the head of the pressure tube to be calculated is fixed, and the pressure tube grows towards the tail under irradiation; when the cumulative irradiation duration reaches the reconfiguration irradiation time, the head of the pressure tube to be calculated is loosened and the tail is fixed, and the pressure tube grows towards the head under irradiation.

[0035] In step S120, the viscoplastic self-consistent algorithm is executed by script, and the constitutive model parameters output by the viscoplastic self-consistent algorithm are saved. The constitutive model parameters include thermal creep anisotropy coefficient, irradiation creep anisotropy coefficient, and irradiation growth strain coefficient.

[0036] The inputs required by the Visco-Plastic Self-Consistent (VPSC) algorithm can be stored in files, including parameter configuration files, pressure tube texture files, single-crystal constitutive files, and loading parameter files. The parameter configuration file stores the single-crystal viscoplastic constitutive parameters, including slip system type, strain rate sensitivity index, and hardening parameters. The pressure tube texture file stores the texture information of the pressure tube material. The single-crystal constitutive file stores the relationship between deformation strain and stress (corresponding to creep behavior) or the relationship between deformation strain and irradiation flux (corresponding to irradiation growth behavior). The loading parameter file stores the macroscopic loading boundaries, including strain boundary conditions or stress boundary conditions.

[0037] In the process of executing the viscoplastic self-consistent algorithm via script, firstly, the script drives the reading of parameter configuration files, pressure tube texture files, single crystal constitutive files, and loading parameter files to obtain the input data required for the viscoplastic self-consistent algorithm; then, the script drives the execution of the viscoplastic self-consistent algorithm based on the input data to calculate the constitutive model parameters that reflect the anisotropic deformation behavior of the material.

[0038] The constitutive model parameters output by the viscoplastic self-consistent algorithm include thermal creep anisotropy coefficients, irradiation creep anisotropy coefficients, and irradiation growth strain coefficients. These parameters are saved to a readable storage medium via a script for later use in generating finite element geometric models. For example, the script saves the thermal creep anisotropy coefficients to the ANISOTH.OUT file, the irradiation creep anisotropy coefficients to the ANISOIR.OUT file, and the irradiation growth strain coefficients to the GROWIR.OUT file, and then stores these files to a readable storage medium.

[0039] This application does not limit the specific implementation of the viscoplastic self-consistent algorithm. Any viscoplastic self-consistent algorithm that can output constitutive model parameters based on single-crystal constitutive relations is within the protection scope of this application. In specific implementation, a matching viscoplastic self-consistent algorithm can be selected according to the analysis requirements.

[0040] In some embodiments, the viscoplastic self-consistent algorithm includes a first algorithm and a second algorithm. The first algorithm is used to calculate the thermal creep anisotropy coefficient, and the second algorithm is used to calculate the irradiation creep anisotropy coefficient and the irradiation growth strain coefficient.

[0041] refer to Figure 2 As shown, the calculation process of the first algorithm includes: Step S201: For a single grain, under known macroscopic stress... In the case of assuming the stress of the grain Step S202: Calculate the strain rate of the grain. viscoplastic creep compliance tensor Post-strain rate term Step S203: Calculate the macroscopic strain rate and macroscopic creep compliance tensor Step S204: Based on the macroscopic creep compliance tensor And grain shape calculation of Eshelby tensor; Step S205: Calculate interaction tensor and localized tensors , Step S206: Recalculate the macroscopic creep compliance tensor based on the self-consistency condition. Step S207: Determine the macroscopic creep compliance tensor Whether it converges; if it does not converge, jump to step S204; otherwise, proceed to the next step; Step S208: Recalculate the macroscopic strain rate using the macroscopic constitutive model. Step S209: Recalculate the stress and strain rate of the grain based on the self-consistent condition to obtain the updated stress. and strain rate Step S210: Determine stress and strain rate If convergence fails, proceed to step S202; otherwise, proceed to the next step. Step S211: Update grain orientation and shape, and update rotation tensor. Step S212: Update critical resolved shear stress (CRSS). Step S213: Switch to the next grain and repeat steps S201 to S213.

[0042] refer to Figure 3 As shown, the calculation process of the second algorithm includes: Step S301: For a single grain, under known macroscopic stress... In the case of assuming the stress of the grain Step S302: Calculate the irradiation growth rate of the grain. strain rate viscoplastic creep compliance tensor Post-strain rate term Step S303: Calculate the macroscopic strain rate and macroscopic creep compliance tensor Step S304: Based on the macroscopic creep compliance tensor And grain shape calculation of Eshelby tensor; Step S305: Calculate interaction tensor and localized tensors , Step S306: Recalculate the macroscopic creep compliance tensor based on the self-consistency condition. Step S307: Determine the macroscopic creep compliance tensor Whether it converges; if it does not converge, jump to step S304; otherwise, proceed to the next step; Step S308: Recalculate the macroscopic strain rate using the macroscopic constitutive model. Step S309: Recalculate the stress and strain rate of the grain based on the self-consistent condition to obtain the updated stress. and strain rate Step S310: Determine stress and strain rate If convergence is not achieved, proceed to step S302; otherwise, proceed to the next step. Step S311: Update the grain orientation and shape, and update the rotation tensor. Step S312: Switch to the next grain and re-execute steps S301 to S312.

[0043] In one embodiment, step S120 further includes Figure 4 The steps S121 to S123 are shown below. In step S121, the input to the first algorithm is obtained through a script. The input to the first algorithm includes a first parameter configuration file, a pressure tube texture file, a first single crystal constitutive file, and a first loading parameter file. In step S122, the input to the second algorithm is obtained through a script. The input to the second algorithm includes a second parameter configuration file, a pressure tube texture file, a second single crystal constitutive file, and a second loading parameter file. In step S123, the first algorithm and the second algorithm are executed through a script, and the thermal creep anisotropy coefficient output by the first algorithm, as well as the irradiation creep anisotropy coefficient and irradiation growth strain coefficient output by the second algorithm, are saved.

[0044] In step S130, the input data and constitutive model parameters are organized and saved as an input file using a script.

[0045] In this step, the input data and constitutive model parameters are organized into input files with a preset format so that the input files can be submitted to the first interface and used to generate the finite element geometric model. The preset format refers to the input format supported by the first interface. The preset format can be designed according to the type of the first interface. For example, when the first interface is provided by the finite element analysis software ABAQUS, the input files include files with the .f and .inp extensions. In these file formats, the preset format includes key-value pairs, along with the rules for writing these files.

[0046] In practice, firstly, the script reads the file saved in the readable storage medium to obtain the input data obtained in step S110 and the constitutive model parameters calculated in step S120; then, the input data and constitutive model parameters are organized according to a preset format, including embedding the values ​​of each parameter into the corresponding keys; then, the organized input file is saved to the readable storage medium as a prerequisite for the finite element calculation of the pressure tube irradiation deformation, for use in subsequent steps.

[0047] As an example of step S130, the input data and constitutive model parameters are processed by a script to generate the files para_zirconium.f and Pressure_tube.inp. The para_zirconium.f file contains records of the constitutive model parameters, which are used by the first interface to model the anisotropic deformation behavior of the pressure tube to be calculated. The Pressure_tube.inp file contains records of the input data, which are used by the first interface to model the geometry, irradiation history, and operating conditions of the pressure tube to be calculated.

[0048] In some embodiments, such as Figure 5 As shown, step S130 further includes steps S131 and S132. In step S131, the thermal creep anisotropy coefficient, irradiation creep anisotropy coefficient, and irradiation growth strain coefficient, respectively corresponding to the tube head and tube tail of the pressure tube to be calculated, are obtained by the user from the constitutive model parameters via a script. In step S132, the input file is generated and saved by the script based on the input data, the thermal creep anisotropy coefficient, irradiation creep anisotropy coefficient, and irradiation growth strain coefficient corresponding to the tube head and tube tail of the pressure tube to be calculated, respectively. For example, the three types of coefficients corresponding to the tube head and tube tail are written into the corresponding tube head and tube tail record areas of the input file according to a preset format via a script.

[0049] For example, the script obtains the user's selected operation through an interactive interface, which includes checkboxes and / or text input boxes. The user can select the file path of the saved constitutive model parameters through the checkboxes and / or enter the file path through the text input boxes. The script reads the corresponding three types of parameters according to the file path and writes the three types of coefficients corresponding to the tube head and tube tail respectively into the corresponding tube head and tube tail record areas of the input file according to a preset format.

[0050] In step S140, the first interface is called through the script to obtain and save the finite element geometric model generated based on the input file and the pre-configuration file.

[0051] In practice, a script calls the first interface, submitting the input file and pre-configuration file to it. Upon receiving the input file and pre-configuration file, the first interface executes the irradiation deformation macroscopic semi-empirical C6 constitutive equation to generate a finite element geometric model. The expression for the irradiation deformation macroscopic semi-empirical C6 constitutive equation is as follows:

[0052] In the formula, express strain rate in the direction ( (representing the radial, lateral, and axial dimensions of the pressure tube to be calculated, respectively), units are... , , and They represent The strain rates of thermal creep, irradiation creep, and irradiation growth are considered. The expression for the strain rate of thermal creep is as follows:

[0053] The expression for the strain rate of irradiation creep is as follows:

[0054] The expression for the strain rate during irradiation growth is as follows:

[0055] in, It represents temperature, and the unit is Kelvin (K). Represents fast neutron flux, in units of , This indicates the irradiation time, expressed in seconds (s). , , , , and It is a constant. and Represents the anisotropy parameter of thermal creep. Indicates the anisotropy parameters of irradiated growth, where It can be calculated based on the texture and microstructure of the pressure tube to be calculated. and This represents the equivalent stress due to thermal creep, expressed in MPa. Indicates the irradiation creep rate as a function of axial position A changing function, Indicates the anisotropy parameter of irradiation creep as a function of axial position. A changing function, This represents the equivalent stress of irradiation creep, which is determined by radial stress, transverse stress, and axial stress, respectively. Indicates the irradiation growth rate as a function of axial position. and irradiation time A changing function, , , and This indicates the activation temperature, measured in Kelvin (K).

[0056] The pre-configuration file includes preset values ​​for some parameters of the semi-empirical C6 constitutive equation for macroscopic irradiation deformation. It can be understood that by solving the aforementioned semi-empirical C6 constitutive equation for macroscopic irradiation deformation, the radial, transverse, and axial strains of the pressure tube to be calculated are obtained to generate a finite element geometric model. Furthermore, based on the finite element geometric model and combined with the original structural dimensions of the pressure tube to be calculated, the deformation results of the pressure tube (such as wall thinning, inner diameter increase, and axial growth) can be obtained.

[0057] In some embodiments, the first interface includes a command-line interface for finite element analysis software. The command-line interface refers to the interface provided by the finite element analysis software that initiates the finite element analysis process by running commands. Through the command-line interface, input files and pre-configuration files can be batch-processed into the finite element analysis software without opening a graphical interface, and the software can then be started to perform modeling and calculation tasks.

[0058] In this embodiment, step S140 further includes Figure 6 The steps S141 and S142 are shown. In step S141, the input file and pre-configuration file are submitted to the command-line interface via a script; in step S142, the finite element geometric model returned by the command-line interface is obtained and saved via a script. For example, when the first interface is the command-line interface provided by the finite element analysis software ABAQUS, the first interface takes the following form: @echo off abaqus job=ModelName user=UmatSubFile cpus=cpu_count interactive pause” The keyword `job` specifies the name of the analysis task, the keyword `user` specifies the filename of the subroutine (i.e., the pre-configuration file), and the keyword `cpus` specifies the number of CPU cores used for parallel solving.

[0059] As an example, a script assembles commands, including input files and pre-configuration files, according to the first interface, and executes these commands in the runtime environment to start the finite element analysis software. Using the finite element autonomous analysis tool UMAT (User Material Subroutine), the win32api.ShellExecute command is executed to call the calculation program of the ABAQUS finite element analysis software. Based on the read input file and pre-configuration file, the finite element solution for thermal creep, irradiation creep, and irradiation growth constitutive relations in the macroscopic semi-empirical C6 constitutive equation of irradiation deformation is achieved, generating a finite element geometric model (odb geometric model), which is returned as an output file (.odb database file). During the finite element solution process, the execution status of the command-line interface is monitored by a script, and the output file of the finite element geometric model is obtained and saved after execution.

[0060] It is understood that embodiments of this application may use the modeling interfaces provided by different finite element analysis software as the first interface. When changing finite element analysis software, only the modeling interface provided by the new finite element analysis software needs to be used as the first interface, without modifying the script, thereby improving the efficiency of replacement and code reusability.

[0061] In step S150, the second interface is called through the script to obtain and save multiple sets of deformation data corresponding to multiple measurement points, calculated based on the finite element geometric model.

[0062] Specifically, the script calls the second interface to submit the finite element geometric model generated and saved in step S140. After receiving the finite element geometric model, the second interface performs post-processing calculations to obtain deformation data corresponding to multiple measurement points. The script then saves the multiple sets of deformation data returned by the second interface as deformation result files in a readable storage medium. Each set of deformation data corresponds to a measurement point and includes at least one of the following at the measurement point: pipe wall thinning value, inner diameter increase value, and axial growth value. The second interface may include one or more sub-interfaces, each used to calculate different deformation data based on the finite element geometric model. The deformation result files can be in .xlsx or .txt format. Furthermore, the script can also save the output multiple sets of deformation results to an SQLite3 database for easy subsequent reading and querying.

[0063] In some embodiments, the second interface includes a cross-sectional strain calculation interface, a cross-sectional stress calculation interface, an axial growth calculation interface, and an inner and outer diameter dimension change calculation interface of the finite element analysis software. In this embodiment, step S150 further includes... Figure 7The steps S151 to S155 are shown below. In step S151, the finite element geometric model is submitted to the cross-sectional strain calculation interface via a script, and the output cross-sectional strain corresponding to multiple measurement points is obtained. In step S152, the finite element geometric model is submitted to the cross-sectional stress calculation interface via a script, and the output cross-sectional stress corresponding to multiple measurement points is obtained. In step S153, the finite element geometric model is submitted to the axial growth calculation interface via a script, and the output axial growth value corresponding to multiple measurement points is obtained. In step S154, the finite element geometric model is submitted to the inner and outer diameter dimension change calculation interface via a script, and the output inner and outer diameter dimension changes corresponding to multiple measurement points are obtained. In step S155, the cross-sectional strain, cross-sectional stress, axial growth value, and inner and outer diameter dimension changes corresponding to multiple measurement points are saved via a script.

[0064] For example, the finite element geometric model generated in step S140 is an odb geometric model generated by the finite element analysis software ABAQUS and stored in a .odb database file. The second interface is an interface implemented in Python that can call the odbAccess module provided by ABAQUS. The second interface opens the .odb database file storing the odb geometric model through the odbAccess module, extracts the strain components, stress components, axial displacement, and radial displacement corresponding to each measurement point, and calculates the deformation data corresponding to each measurement point based on the original structural dimensions of the pressure tube to be calculated.

[0065] As another embodiment of the heavy water reactor pressure tube irradiation deformation modeling method of this application, such as Figure 8 As shown, in addition to steps S110 to S150 described above, it also includes step S160: constructing a visual interface using a script. (See reference...) Figure 9 The illustrated flowchart shows a method for constructing a visual interface using scripts, which further includes... Figure 9 Steps S161 to S163 in the process.

[0066] In step S161, a first coordinate system and a second coordinate system are established by a script. The horizontal axis of the first coordinate system and the second coordinate system represent multiple measurement points, the vertical axis of the first coordinate system represents the pipe wall thinning value, and the vertical axis of the second coordinate system represents the inner diameter increase value.

[0067] Specifically, the script calls the plotting functions of the visualization plotting library to establish a first coordinate system and a second coordinate system on two separate canvases. Both the first and second coordinate systems use the identifiers of multiple measurement points (such as measurement point numbers) as the scales on the horizontal axis to distinguish measurement points at different axial positions. The vertical axes of both the first and second coordinate systems are numerical axes. The vertical axis of the first coordinate system represents the pipe wall thinning value, reflecting the degree of wall thinning at each measurement point. The vertical axis of the second coordinate system represents the inner diameter increase value, reflecting the degree of inner diameter increase at each measurement point. By establishing two independent coordinate systems, the pipe wall thinning value and the inner diameter increase value are displayed in their respective coordinate systems, avoiding the mixing of different types of deformation data in the same coordinate system.

[0068] In step S162, the script plots the coordinates of the pipe wall thinning value from multiple sets of deformation data corresponding to each measurement point in the first coordinate system. Specifically, the script reads the saved multiple sets of deformation data and plots the coordinates in the first coordinate system with the measurement point identifier as the abscissa and the corresponding pipe wall thinning value as the ordinate. The coordinates in the first coordinate system reflect the changing trend of the pipe wall thinning value along the axial direction of the pressure pipe.

[0069] In step S163, a script is used to plot the coordinates of the inner diameter growth value from multiple sets of deformation data corresponding to each measurement point in the second coordinate system. Specifically, the script reads the saved multiple sets of deformation data and plots the coordinates in the second coordinate system with the measurement point identifier as the abscissa and the corresponding inner diameter growth value as the ordinate. The coordinates in the second coordinate system can reflect the changing trend of the inner diameter growth value of the pressure tube along the axial direction.

[0070] As an example, this code uses a Python script to call the Matplotlib open-source plotting library to create scatter plots of pipe wall thinning and inner diameter growth. Specifically, it uses the pyplot module to create two subplots, corresponding to the first and second coordinate systems respectively. In the subplot for the first coordinate system, the `Scatter` function is used to plot points with the measurement point's identifier as the x-axis and the corresponding pipe wall thinning value as the y-axis. In the subplot for the second coordinate system, the `Scatter` function is used to plot points with the measurement point's identifier as the x-axis and the corresponding inner diameter growth value as the y-axis. Furthermore, the color, marker style, and labels of the coordinate points can be set using keywords in the `Scatter` function.

[0071] It should be noted that, in addition to Figure 9In addition to the embodiments shown, any method that can construct a visual interface through scripts to enable users to intuitively obtain deformation data is within the scope of protection of this application. For example, multiple tables can be created on the interface through scripts, and the deformation data corresponding to each measurement point can be filled into the corresponding cells. Another example is to create a three-dimensional diagram of the pressure tube to be calculated on the interface through scripts, and display the deformation data corresponding to each measurement point at the corresponding position in the three-dimensional diagram.

[0072] In some embodiments, the method of building a visual interface via scripts includes, in addition to, Figure 9 The steps S161 to S163 shown also include Figure 10 Steps S164 and S165 are executed after step S163.

[0073] In step S164, user interaction actions are acquired via a script. Specifically, after drawing is completed, the script continuously monitors the user's interaction actions in the first and second coordinate systems. User interaction actions include, but are not limited to, clicking coordinate points and hovering the mouse over coordinate points.

[0074] In step S165, when the user interaction includes selecting a coordinate point, the script displays the wall thinning value and inner diameter increase value of the measurement point corresponding to the selected coordinate point in the first coordinate system and the second coordinate system, respectively. Specifically, the script identifies the user interaction; when the user interaction includes selecting a coordinate point (e.g., clicking the coordinate point or hovering the mouse over the coordinate point), the script obtains the measurement point corresponding to that coordinate point and reads multiple sets of saved deformation data to obtain the wall thinning value and inner diameter increase value corresponding to that measurement point. In the first coordinate system, the corresponding wall thinning value is displayed next to the coordinate point with the identifier of the measurement point as the abscissa. In the second coordinate system, the corresponding inner diameter increase value is displayed next to the coordinate point with the identifier of the measurement point as the abscissa.

[0075] It is understood that the heavy water reactor pressure tube irradiation deformation modeling method of this application embodiment can model and analyze a single pressure tube to be calculated. Furthermore, the heavy water reactor pressure tube irradiation deformation modeling method of this application embodiment can also model and analyze multiple pressure tubes to be calculated in batches. For example, taking the number of pressure tubes to be calculated as a unit, the script uses the QTime component in PySide2.QtCore to implement independent loop calculations for a single pressure tube to be calculated, using the calculation status of multiple sets of deformation data as the judgment condition.

[0076] Some embodiments of this application also propose a heavy water reactor pressure tube irradiation deformation modeling system (hereinafter referred to as the "modeling system"). For example... Figure 11As shown, the modeling system 400 includes a data input module 411, a parameter calculation module 412, a file generation module 413, and a deformation calculation module 414. The data input module 411 is configured to acquire and save input data. The parameter calculation module 412 is configured to execute a viscoplastic self-consistent algorithm and save the constitutive model parameters output by the algorithm. The file generation module 413 is configured to organize and save the input file based on the user-selected, saved constitutive model parameters and input data. The deformation calculation module 414 is configured to call a first interface to generate and save a finite element geometric model based on the input file and a pre-configuration file, and to call a second interface to calculate and save multiple sets of deformation data corresponding to multiple measurement points based on the finite element geometric model.

[0077] In some embodiments, such as Figure 11 As shown, the modeling system 400 also includes a front-end page for user interaction. The front-end page includes a first page 421, a second page 422, a third page 423, and a fourth page 424. The first page 421 corresponds to the data input module 411, used to receive input data and transfer it to the data input module 411 for saving. The second page 422 corresponds to the parameter calculation module 412, used to receive the save address of the file selected by the user and transfer this save address to the parameter calculation module 412 to execute the viscoplastic self-consistent algorithm. Furthermore, the second page 422 also displays the data saved in the user-selected file. The third page 423 corresponds to the file generation module 413, used to receive the constitutive model parameters and the save address of the input data selected by the user, and transfer this save address to the file generation module 413 to generate an input file. Furthermore, the third page 423 also displays the data saved in the input file. Page 424 corresponds to the deformation calculation module 414. It is used to receive input files and call the first and second interfaces in sequence to perform irradiation deformation modeling and analysis, and output multiple sets of deformation data corresponding to multiple measurement points.

[0078] In other embodiments, the modeling system 400 further includes a visualization module 415. The visualization module 415 is configured to perform steps S161 to S163 as described above to establish a first coordinate system and a second coordinate system in the fourth page 424, and to display the pipe wall thinning value and inner diameter growth value corresponding to each measurement point in the first coordinate system and the second coordinate system, respectively.

[0079] refer to Figure 12 The diagram shown is a schematic of the first page corresponding to the data input module of one embodiment. Figure 13 The diagram shown is a schematic of the second page corresponding to the parameter calculation module in one embodiment. Figure 14 The diagram shown illustrates the third page corresponding to the file generation module in one embodiment. Figure 15This diagram illustrates the fourth page corresponding to the deformation calculation module and visualization module in one embodiment. The first to fourth pages are interactive visualization pages developed based on the PySide2 open-source library.

[0080] like Figure 12 As shown, the title of the first page is "Data Input," corresponding to the data input module. On the left side of the first page are text input boxes for pipe length, inner wall diameter (mm), outer wall diameter (mm), continuous irradiation time (h), reconstructed irradiation time (h), flux distribution (n / m²s), temperature distribution (K), and pressure distribution (MPa) corresponding to the 12 measurement points, as well as buttons for "Load Data," "Preview," "Save," "Save As," and "Cancel." On the right side of the first page is a preview box for the input data. Users can input the above eight types of data sequentially through the first page and save the input data to a readable storage medium using the "Save" button; alternatively, they can load the saved eight types of data from a readable storage medium using the "Load Data" button.

[0081] like Figure 13 As shown, the title of the second page is "Parameter Calculation," corresponding to the parameter calculation module. On the left side of the second page are text input boxes for the save paths of the first parameter configuration file, pressure tube texture file, first single-crystal constitutive file, first loading parameter file, second parameter configuration file, second single-crystal constitutive file, and second loading parameter file, respectively, as well as "Load" and "Cancel" buttons. On the right side of the second page is a preview box displaying the data saved in the file where the mouse is hovering. If the data in the selected file is correct, the user can click the "Load" button to calculate and save the thermal creep anisotropy coefficient, irradiation creep anisotropy coefficient, and irradiation growth strain coefficient; otherwise, clicking the "Cancel" button will delete the selected file and allow for a new selection.

[0082] like Figure 14 As shown, the title of the third page is "File Generation," corresponding to the file generation module. On the left side of the third page, text input boxes are provided for the save paths of the input data and constitutive model parameters, respectively, along with "Query," "OK," "f-inp Generation," and "Cancel" buttons. In this embodiment, the user can select the constitutive model parameters for the "pipe head" and "pipe tail" according to their modeling needs. After selection, clicking the "OK" button will display the selected input data and constitutive model parameters on the right side of the third page. "f-inp Generation" indicates that the input files (such as the para_zirconium.f file and the Pressure_tube.inp file) are generated and saved. Clicking "f-inp Generation" generates and saves the input files, providing input for the deformation calculation module.

[0083] like Figure 15 As shown, the title of the fourth page is "Deformation Calculation," corresponding to the deformation calculation module and the visualization module. On the left side of the fourth page, there is a text input box for the save path of the input files, as well as "Calculate," "Cancel," "Report," and "Backup" buttons. The input files include para_zirconium.f and Pressure_tube.inp. The para_zirconium.f file contains records of constitutive model parameters, which are used by the first interface to model the anisotropic deformation behavior of the pressure tube to be calculated. The Pressure_tube.inp file contains records of input data, which are used by the first interface to model the geometry, irradiation history, and operating conditions of the pressure tube to be calculated. After loading the save path of the input files, the user can click the "Calculate" button to call the first interface, which generates and saves the finite element geometric model based on the input files and pre-configuration files, and then call the second interface to calculate and save multiple sets of deformation data corresponding to 12 measurement points based on the finite element geometric model. A visualization area is provided on the right side of the fourth page. The visualization module establishes a first coordinate system with the horizontal axis representing the measurement point marker and the vertical axis representing the pipe wall thinning value, and a second coordinate system with the horizontal axis representing the measurement point marker and the vertical axis representing the inner diameter increase value. Multiple sets of deformation data are then displayed in the first and second coordinate systems.

[0084] Next, refer to Figure 15 When the user hovers the mouse over a coordinate point, the visualization module captures the user's interaction and displays the pipe wall thinning value and inner diameter increase value corresponding to that coordinate point in both the first and second coordinate systems. Furthermore, the lower part of the fourth page displays data such as transverse cross-sectional strain, axial cross-sectional strain, transverse cross-sectional stress, axial cross-sectional stress, and axial growth (mm) for the corresponding measurement point.

[0085] Based on the first to fourth pages described above, users can conveniently, quickly, and accurately model and analyze the irradiation deformation of heavy water reactor pressure tubes. It is understood that the above-described implementation is merely a typical example of this application, and not the only feasible implementation. For example, implementing the front-end page using Java language files and the back-end module using Python script files is also a feasible implementation of this application.

[0086] An embodiment of this application also proposes a method such as Figure 16 The electronic device 500 is shown. According to... Figure 16 The electronic device 500 may include an internal communication bus 501, a processor 502, a read-only memory 503, a random access memory 504, and a communication port 505. When used in a personal computer, the electronic device may also include a hard disk 506.

[0087] The internal communication bus 501 enables data communication between components of the electronic device 500. The processor 502 can perform judgments and issue prompts. In some embodiments, the processor 502 may consist of one or more processors. The communication port 505 enables data communication between the electronic device 500 and external devices. In some embodiments, the electronic device 500 can send and receive information and data from a network through the communication port 505.

[0088] Electronic device 500 may also include different forms of program storage units and data storage units, such as hard disk 506, read-only memory 503, and random access memory 504, capable of storing various data files used for computer processing and / or communication, as well as possible program instructions executed by processor 502. The processor executes these instructions to implement the main parts of the methods described above. The results of processor processing are transmitted to the user equipment via a communication port and displayed on the user interface.

[0089] In another aspect, this application also proposes a computer storage medium storing computer program code that, when executed by a processor, implements the method described above.

[0090] The basic concepts have been described above. Obviously, for those skilled in the art, the above disclosure is merely illustrative and does not constitute a limitation of this application. Although not explicitly stated herein, those skilled in the art may make various modifications, improvements, and corrections to this application. Such modifications, improvements, and corrections are suggested in this application, and therefore remain within the spirit and scope of the exemplary embodiments of this application.

[0091] Furthermore, this application uses specific terms to describe embodiments of the application. For example, "an embodiment," "one embodiment," and / or "some embodiments" refer to a particular feature, structure, or characteristic related to at least one embodiment of the application. Therefore, it should be emphasized and noted that "an embodiment," "one embodiment," or "an alternative embodiment" mentioned twice or more in different locations in this specification do not necessarily refer to the same embodiment. In addition, certain features, structures, or characteristics in one or more embodiments of the application can be appropriately combined.

[0092] Similarly, it should be noted that, in order to simplify the description of the present application and thus aid in the understanding of one or more embodiments, the foregoing description of the embodiments of the present application sometimes combines multiple features into a single embodiment, drawing, or description thereof. However, this disclosure method does not imply that the subject matter of the present application requires more features than those mentioned in the claims. In fact, the embodiments contain fewer features than all the features of the single embodiments disclosed above.

[0093] In some embodiments, numbers describing the quantity of components and attributes are used. It should be understood that such numbers used in the description of embodiments are modified in some examples with the terms "approximately," "approximately," or "generally." Unless otherwise stated, "approximately," "approximately," or "generally" indicates that the numbers are allowed to vary by ±20%. Accordingly, in some embodiments, the numerical parameters used in the specification and claims are approximate values, which may be changed depending on the characteristics required by individual embodiments. In some embodiments, numerical parameters should take into account specified significant digits and employ a general method of digit reservation. Although the numerical ranges and parameters used to confirm their breadth of scope in some embodiments of this application are approximate values, in specific embodiments, such values ​​are set as precisely as feasible.

[0094] Although this application has been described with reference to specific embodiments, those skilled in the art should recognize that the above embodiments are only used to illustrate this application, and various equivalent changes or substitutions can be made without departing from the spirit of this application. Therefore, any changes or modifications to the above embodiments within the essential spirit of this application will fall within the scope of the claims of this application.

Claims

1. A method for modeling the deformation of a pressure tube in a heavy water reactor under irradiation, characterized in that, include: The script acquires and saves input data, which includes the geometric parameters, irradiation parameters, and multiple sets of operating condition parameters of the pressure tube to be calculated. The multiple sets of operating condition parameters include the flux distribution, temperature distribution, and pressure distribution at multiple measurement points of the pressure tube to be calculated. The script executes the viscoplastic self-consistent algorithm and saves the constitutive model parameters output by the viscoplastic self-consistent algorithm. The constitutive model parameters include the thermal creep anisotropy coefficient, the irradiation creep anisotropy coefficient, and the irradiation growth strain coefficient. The script organizes and saves the input data and constitutive model parameters into an input file. The script calls the first interface to obtain and save the finite element geometric model generated based on the input file and the pre-configuration file; and The script calls the second interface to obtain and save multiple sets of deformation data corresponding to the multiple measurement points, calculated based on the finite element geometric model.

2. The heavy water reactor pressure tube irradiation deformation modeling method as described in claim 1, characterized in that, It also includes building a visualization interface through the script, the step of building the visualization interface through the script includes: The script establishes a first coordinate system and a second coordinate system. The horizontal axis of the first coordinate system and the second coordinate system represents the plurality of measurement points. The vertical axis of the first coordinate system represents the pipe wall thinning value, and the vertical axis of the second coordinate system represents the inner diameter increase value. The script is used to plot the coordinates of the pipe wall thinning value in the multiple sets of deformation data corresponding to each measurement point in the first coordinate system. The script is used to plot the coordinates of the inner diameter growth value in the multiple sets of deformation data corresponding to each measurement point in the second coordinate system.

3. The method for modeling the deformation of a pressure tube in a heavy water reactor under irradiation as described in claim 2, characterized in that, The steps for constructing a visual interface using the script also include: The script is used to obtain user interaction actions; When the user interaction action includes selecting a coordinate point, the script displays the pipe wall thinning value and inner diameter increase value of the measurement point corresponding to the selected coordinate point in the first coordinate system and the second coordinate system, respectively.

4. The method for modeling the deformation of a pressure tube in a heavy water reactor under irradiation as described in claim 1, characterized in that, The geometric parameters of the pressure tube to be calculated include the tube length, inner diameter, and outer diameter, and the irradiation parameters of the pressure tube to be calculated include the continuous irradiation time and the reconstruction irradiation time.

5. The method for modeling the deformation of a pressure tube in a heavy water reactor under irradiation as described in claim 1, characterized in that, The first interface includes a command-line interface for finite element analysis software. The steps of generating and saving a finite element geometric model based on the input file and pre-configuration file by calling the first interface through the script include: The script submits the input file and the pre-configuration file to the command-line interface. The script retrieves and saves the finite element geometric model returned by the command-line interface.

6. The method for modeling the deformation of a pressure tube in a heavy water reactor under irradiation as described in claim 1, characterized in that, The second interface includes interfaces for calculating cross-sectional strain, cross-sectional stress, axial growth, and inner / outer diameter changes in finite element analysis software. The steps of calling the second interface through the script to calculate and save multiple sets of deformation data corresponding to the multiple measurement points based on the finite element geometric model include: The script submits the finite element geometric model to the cross-sectional strain calculation interface and obtains the cross-sectional strain corresponding to the multiple measurement points. The script submits the finite element geometric model to the section stress calculation interface and obtains the output section stress corresponding to the multiple measurement points. The script submits the finite element geometric model to the axial growth calculation interface and obtains the axial growth values ​​corresponding to the multiple measurement points. The script submits the finite element geometric model to the interface for calculating the inner and outer diameter changes, and obtains the output inner and outer diameter changes corresponding to the multiple measurement points. The script saves the cross-sectional strain, cross-sectional stress, axial growth value, and inner and outer diameter changes corresponding to the multiple measurement points.

7. The method for modeling the deformation of a pressure tube in a heavy water reactor under irradiation as described in claim 1, characterized in that, The viscoplastic self-consistent algorithm includes a first algorithm and a second algorithm. The first algorithm is used to calculate the thermal creep anisotropy coefficient, and the second algorithm is used to calculate the irradiation creep anisotropy coefficient and the irradiation growth strain coefficient. The steps of executing the viscoplastic self-consistent algorithm through the script and saving the constitutive model parameters output by the viscoplastic self-consistent algorithm include: The input to the first algorithm is obtained through the script. The input to the first algorithm includes a first parameter configuration file, a pressure tube texture file, a first single crystal constitutive file, and a first loading parameter file. The input to the second algorithm is obtained through the script. The input to the second algorithm includes a second parameter configuration file, the pressure tube texture file, a second single crystal constitutive file, and a second loading parameter file. The script executes the first algorithm and the second algorithm respectively, and saves the thermal creep anisotropy coefficient output by the first algorithm, as well as the irradiation creep anisotropy coefficient and irradiation growth strain coefficient output by the second algorithm.

8. The method for modeling the deformation of a pressure tube in a heavy water reactor as described in claim 1, characterized in that, The steps of organizing and saving the input data and constitutive model parameters into an input file using the script include: The script obtains the thermal creep anisotropy coefficient, irradiation creep anisotropy coefficient, and irradiation growth strain coefficient, respectively, selected by the user from the constitutive model parameters, corresponding to the tube head and tube tail of the pressure tube to be calculated. The script generates and saves the input file based on the input data, the thermal creep anisotropy coefficients, the irradiation creep anisotropy coefficients, and the irradiation growth strain coefficients corresponding to the pipe head and tail of the pressure tube to be calculated, respectively.

9. A heavy water reactor pressure tube irradiation deformation modeling system, characterized in that, include: The data input module is configured to acquire and save input data, which includes the geometric parameters, irradiation parameters and multiple sets of operating condition parameters of the pressure tube to be calculated, wherein the multiple sets of operating condition parameters include the injection distribution, temperature distribution and pressure distribution of multiple measurement points; The parameter calculation module is configured to execute the viscoplastic self-consistent algorithm and save the constitutive model parameters output by the viscoplastic self-consistent algorithm, including the thermal creep anisotropy coefficient, the irradiation creep anisotropy coefficient, and the irradiation growth strain coefficient. The file generation module is configured to organize and save the input file based on the user-selected, saved constitutive model parameters and input data; and The deformation calculation module is configured to call a first interface to generate and save a finite element geometric model based on the input file and pre-configuration file, and to call a second interface to calculate and save multiple sets of deformation data corresponding to the multiple measurement points based on the finite element geometric model.

10. The heavy water reactor pressure tube irradiation deformation modeling system as described in claim 9, characterized in that, It also includes a visualization module, which is configured as follows: A first coordinate system and a second coordinate system are established using a visualization drawing library. The horizontal axis of the first coordinate system and the second coordinate system represents the plurality of measurement points. The vertical axis of the first coordinate system represents the pipe wall thinning value, and the vertical axis of the second coordinate system represents the inner diameter increase value. Plot the coordinates of the pipe wall thinning value in the multiple sets of deformation data corresponding to each measurement point in the first coordinate system; Plot the coordinates of the inner diameter growth value in the multiple sets of deformation data corresponding to each measurement point in the second coordinate system.

11. The heavy water reactor pressure tube irradiation deformation modeling system as described in claim 10, characterized in that, The visualization module is also configured to: Get user interaction actions; When the user interaction action includes selecting a coordinate point, the wall thinning value and inner diameter increase value of the measurement point corresponding to the selected coordinate point are displayed in the first coordinate system and the second coordinate system, respectively.

12. An electronic device, characterized in that, include: Memory is used to store instructions that can be executed by the processor; as well as A processor for executing the instructions to implement the method as described in any one of claims 1-8.

13. A computer storage medium storing computer program code, characterized in that, The computer program code, when executed by a processor, implements the method as described in any one of claims 1-8.

Images