Small deformation mechanical analysis system and method based on reactor numerical calculation framework and related products

Abstract

The invention discloses a small deformation mechanical analysis system and method based on a reactor numerical calculation framework and related products, and relates to the technical field of nuclear reactor fuel performance analysis. On the basis of a reactor numerical development framework, a unified tangent modulus and stress data interface is constructed according to a finite element theory and mechanical constitutive models such as elasticity, plasticity and creep; a radial regression algorithm is used for calculating the plasticity / creep strain of a nonlinear material, consistent tangent modulus and stress data are updated and transmitted into a nonlinear solver, a small deformation mechanical analysis system of a reactor numerical development framework is formed, and theoretical model support is provided for a fuel performance analysis program based on small deformation.

Description

Technical Field

[0001] This invention relates to the field of nuclear reactor fuel performance analysis technology, and in particular to small deformation mechanics analysis systems, methods and related products based on reactor numerical calculation framework. Background Technology

[0002] During the operation of a nuclear reactor, structural materials are subjected to extreme conditions such as high temperature, high pressure, and radiation, which causes the material behavior to exhibit significant nonlinear characteristics, such as plastic deformation, creep, and fatigue. This makes it difficult for traditional material models to accurately describe the material behavior, thereby affecting the overall performance and safety of the reactor.

[0003] Due to the complex geometry of reactors, multiphysics coupling effects, and nonlinear variations in material properties, traditional methods often struggle to simultaneously meet the requirements of high reusability, high accuracy, and high efficiency. For example, traditional finite element analysis methods typically employ a single material model, making it difficult to handle constitutive relations of multiple complex materials simultaneously, and resulting in low computational efficiency in multi-condition coupled analysis. While some commercial software supports extending to more complex constitutive relations and material models, the inconsistent calculation formats for stress and tangent modulus across different material models necessitate frequent data conversions during the solver's iteration process, significantly increasing computational overhead. Furthermore, interface conversions may introduce numerical errors, leading to a decrease in computational accuracy.

[0004] Therefore, existing numerical calculation frameworks are insufficient to meet the rapid iteration requirements of multi-material, multi-condition nuclear reactor design. Summary of the Invention

[0005] This invention provides a small deformation mechanics analysis system, method, and related products based on a reactor numerical calculation framework. Building upon this framework, it achieves unified management of consistent tangent modulus and stress data interfaces for mechanical materials. Furthermore, it transmits consistent tangent modulus and stress residuals to the finite element numerical discretization interface according to different mechanical materials. This avoids reduced computational efficiency under multiple materials and operating conditions, accelerates the iteration of the numerical calculation process, and improves reactor design efficiency. This addresses the problem that existing numerical calculation frameworks struggle to adapt to the rapid iteration requirements of multi-material, multi-condition nuclear reactor design.

[0006] This invention is achieved through the following technical solution: A first aspect of the present invention provides a small deformation mechanics analysis system based on a reactor numerical computation framework, comprising: The finite element model module is used to build a finite element calculation model of the reactor. The solver module is used to iteratively solve the finite element calculation model based on Newton's iterative formula to obtain the displacement field of the reactor. The material model module is used to calculate stress and uniform tangent modulus using a radial regression algorithm based on the strain increment in each iteration of the solver module, and the calculation follows a unified data interface specification; the unified data interface specification defines the transmission protocol for stress and uniform tangent modulus. The finite element numerical discretization interface is used to input the stress and uniform tangent modulus calculated by the material model module into the solution engine module, so that the solution engine module can perform the next iteration calculation based on the stress and uniform tangent modulus.

[0007] In some implementations, the finite element model module includes: The geometric modeling module is used to create the geometric model of the reactor; A mesh generation module is used to divide the geometric model into finite element mesh elements; The attribute setting module is used to define the material properties, boundary conditions, and load conditions of the finite element mesh element. The shape function setting module is used to select the element type and shape function of the finite element mesh element to approximately represent the displacement field of the finite element mesh element.

[0008] In some implementations, the solver engine module employs a nonlinear solver to perform nonlinear iterative solving, and the solution steps are as follows: S1, Initialize displacement increment; S2, calculate the strain increment based on the displacement increment; S3, The strain increment is input into the material model module for strain correction to obtain updated stress and consistent tangent modulus; S4, Receive the stress and uniform tangent modulus passed in by the material model module through the finite element numerical discretization interface, and assemble the residual vector and Jacobian matrix based on the stress and uniform tangent modulus, and update the Newton iteration formula; S5, call the nonlinear solver to solve the Newton iteration formula to obtain the updated displacement increment; repeat steps S2~S5 until convergence.

[0009] In some embodiments, the material model module includes a radial regression calculation module for plastic materials and a radial regression calculation module for creep materials; The radial regression calculation module for plastic materials is used to update the stress and uniform tangent modulus of plastic materials; The radial regression calculation module for creep materials is used to update the stress and uniform tangent modulus of creep materials.

[0010] In some embodiments, the material model module further includes a material type determination module, which is used to determine the material type based on material properties; The material model module is also used to update the stress and uniform tangent modulus by calling the radial regression calculation module for plastic materials when the material is a plastic material; and to update the stress and uniform tangent modulus by calling the radial regression calculation module for creep materials when the material is a creeping material.

[0011] A second aspect of the present invention provides a small deformation mechanical analysis method based on a reactor numerical computation framework, applicable to the small deformation mechanical analysis system based on a reactor numerical computation framework as described in any one of the first aspects, wherein the small deformation mechanical analysis method includes: Establish a finite element calculation model of the reactor within the numerical calculation framework of the reactor; The solution engine module is invoked to iteratively solve the finite element calculation model to obtain the displacement field of the reactor. In each iteration of the calculation: The material model module is invoked to calculate the stress and uniform tangent modulus based on the strain increment of the current iteration; The stress and uniform tangent modulus calculated by the material model module are input into the solution engine module through the finite element numerical discretization interface, so that the solution engine module can perform the next iteration calculation based on the stress and uniform tangent modulus.

[0012] In some implementations, establishing the finite element calculation model of the reactor includes: Establish a geometric model of the reactor within the framework of reactor numerical computation; The geometric model is divided into finite element meshes; Define the material properties, boundary conditions, and load conditions of the finite element mesh element; The element type and shape function of the finite element mesh are selected to approximate the displacement field of the finite element mesh.

[0013] In some implementations, the material model module calculates the stress and uniform tangent modulus based on the strain increment calculated in the current iteration as follows: Determine the material type based on its properties; When the material is a plastic material, the radial regression algorithm for plastic materials is used to update the stress and the uniform tangent modulus; When the material is a creeping material, the radial regression algorithm for creeping materials is used to update the stress and the uniform tangent modulus.

[0014] In some embodiments, the radial regression algorithm for plastic materials includes the following steps: S1-1, assuming that the strain increments in the current iteration are all elastic strains, calculate the test stress; S1-2, Calculate the Von-Mise equivalent test stress based on the test stress; S1-3, Based on the equivalent test stress, determine whether the yield function is greater than zero; if it is greater than zero, then correct the plastic strain increment according to S1-4; otherwise, take the assumed strain increment and the test stress as the strain and stress under the real conditions. S1-4, Based on the equivalent test stress, the equivalent plastic strain increment is solved using the Newton method; S1-5, update the plastic consistency tangent modulus.

[0015] In some implementations, the radial regression algorithm for creeping materials includes the following steps: S2-1, assuming that the strain increments in the current step are all elastic strains, calculate the test stress; S2-2, Calculate the Von-Mise equivalent test stress based on the test stress; S2-3, Based on the equivalent test stress, the equivalent creep strain increment is solved using the Newton method; S2-4, Update creep-consistent tangent modulus.

[0016] A third aspect of the present invention provides an electronic device including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor, when executing the computer program, implements the small deformation mechanical analysis method based on a reactor numerical calculation framework as described in any of the second aspects of the present invention.

[0017] A third aspect of the present invention provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the small deformation mechanical analysis method based on a reactor numerical calculation framework as described in any of the second aspects of the present invention.

[0018] A fourth aspect of the present invention provides a computer program product comprising a computer program / instruction that, when executed by a processor, implements the small deformation mechanics analysis method based on a reactor numerical computation framework as described in any of the second aspects of the present invention.

[0019] Compared with the prior art, the present invention has the following advantages and beneficial effects: By unifying the management of consistent tangent modulus and stress data interfaces for mechanical materials, and by transmitting consistent tangent modulus and stress residuals to the finite element numerical discretization interface according to different mechanical materials, the computational efficiency of multi-material and multi-condition calculations is avoided, the iteration of the numerical calculation process is accelerated, and the efficiency of reactor design is improved. By constructing a finite element discrete formula applicable to a variety of mechanical materials, the influence of plastic and creep nonlinear materials on the calculation of consistent tangent modulus and stress is accurately reflected, the mechanism reflected is more comprehensive, and the efficiency and accuracy of stress and strain calculation are improved. Uniform tangent modulus can more accurately reflect the stiffness characteristics of a material under the current stress state, improve the convergence speed of nonlinear finite element analysis, enhance the stability and robustness of numerical models, and provide important support for the efficient solution of complex engineering problems. Attached Figure Description

[0020] To more clearly illustrate the technical solutions of the exemplary embodiments of the present invention, the accompanying drawings used in the embodiments will be briefly described below. It should be understood that the following drawings only show some embodiments of the present invention and should not be considered as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort. In the drawings: Figure 1 This is a schematic diagram of the structure of a small deformation mechanical analysis system proposed in an embodiment of the present invention; Figure 2 This is a flowchart of a nonlinear solution method for a solution engine module proposed in an embodiment of the present invention; Figure 3 This is a flowchart of a deformation mechanics analysis method proposed in an embodiment of the present invention; Figure 4 This is a flowchart of a displacement field comprehensive solution method proposed in an embodiment of the present invention. Detailed Implementation

[0021] To make the objectives, technical solutions, and advantages of the present invention clearer, the present invention will be further described in detail below with reference to the embodiments and accompanying drawings. The illustrative embodiments and descriptions of the present invention are only used to explain the present invention and are not intended to limit the present invention.

[0022] It should be noted that the terms "comprising" and "having" and any variations thereof in the specification, claims, and accompanying drawings of this invention are intended to cover non-exclusive inclusion, for example, a process, method, system, product, or device that includes a series of steps or units is not necessarily limited to other steps or units inherent in the device.

[0023] The terminology used in the various embodiments of the invention is for the purpose of describing particular embodiments only and is not intended to limit the various embodiments of the invention. As used herein, the singular form is intended to include the plural form as well, unless the context clearly indicates otherwise. Unless otherwise defined, all terms used herein (including technical and scientific terms) have the same meaning as commonly understood by one of ordinary skill in the art to which the various embodiments of the invention pertain. The terms (such as those defined in commonly used dictionaries) are to be interpreted as having the same meaning as in the context of the relevant technical field and are not to be interpreted as having an idealized or overly formal meaning, unless clearly defined in the various embodiments of the invention.

[0024] Embodiments of the present invention provide a small deformation mechanics analysis system, method, and related products based on a reactor numerical calculation framework, used for strain / stress calculation of small deformation elastic, plastic, and creep mechanical materials. By constructing a unified and consistent tangent modulus and stress data interface, it can be applied to the finite element discrete formulas of various mechanical materials. This not only avoids the reduction in calculation efficiency for multiple materials and multiple working conditions, but also significantly accelerates the iterative process of numerical calculation, thereby improving the efficiency of reactor design.

[0025] The Workbench for Integrated Nuclear Gengeral-Pdesimulation (WINGS) is a solution for the development of numerical computation software. It constructs a complete numerical computation system process through framework design principles. The small deformation mechanics analysis system of this invention is developed based on the Workbench for Integrated Nuclear Gengeral-Pdesimulation (WINGS).

[0026] Please see Figure 1 , Figure 1 The diagram shows the structure of the small deformation mechanics analysis system proposed in this invention, including a finite element model module, a solver engine module, a material model module, and a finite element numerical discretization interface. The finite element model module can call the finite element analysis module of the existing computing framework to construct the finite element calculation model of the reactor.

[0027] The solver module constructs finite element discretization formulas applicable to various mechanical materials, and performs iterative solutions on the finite element calculation model based on Newton's iteration formula to calculate the displacement field of the reactor. The solver module can call existing nonlinear solvers to execute the iterative solution steps of the displacement field.

[0028] In each iteration of the calculation, the material model module obtains the strain increment of the current iteration step and uses the radial regression algorithm to calculate the stress and uniform tangent modulus.

[0029] A unified data interface for consistent tangent modulus and stress is constructed in the material model module. During each iteration of the solver module, this interface acquires the strain increment of the current iteration and updates the material state variables based on constitutive models of elasticity, plasticity, and creep, using a radial regression algorithm to calculate stress and consistent tangent modulus. By defining a unified data interface in the material model module, the state variable update calculations follow a unified data interface specification. The unified data interface defines the transmission protocol for stress and consistent tangent modulus, such as the tensor format for stress and tangent modulus.

[0030] The finite element numerical discretization interface receives output data from the material model module, transmitting the stress and uniform tangent modulus calculated by the material model module to the solver module. This allows the solver module to perform the next iteration based on these stresses and uniform tangent moduli. By transmitting the uniform tangent modulus and stress residuals to the finite element numerical discretization interface according to different mechanical materials, the computational efficiency is avoided for multiple materials and operating conditions, thus accelerating the iterative numerical calculation process and improving reactor design efficiency.

[0031] Since the behavior of nuclear reactor materials under extreme environments such as high temperature, high pressure and radiation is usually nonlinear, including plastic deformation and creep, this invention uses a nonlinear iterative algorithm to discretize the mechanical equilibrium equations and construct finite element discretization formulas applicable to a variety of mechanical materials.

[0032] First, the weak form of the solid mechanics equilibrium equations is obtained, which are expressed as follows:

[0033] in, Indicates stress, This represents the gradient operator.

[0034] The weak form of the equilibrium equations is expressed as:

[0035] in, n It is the normal vector. This is a test function.

[0036] The mechanical equilibrium equations are discretized using a nonlinear iterative algorithm, and the residual vector in the nonlinear iteration is represented as follows:

[0037] The Jacobian matrix in nonlinear iteration is:

[0038] This method enables a unified approach to handling the nonlinear behavior of different materials under extreme environments such as high temperature, high pressure, and radiation, including plastic deformation and creep, thereby improving the accuracy and reliability of calculations. By using consistent tangent modulus and stress formulas, the complex mechanical properties of materials can be described more accurately, avoiding calculation errors caused by nonlinearity in material models.

[0039] In some implementations, the aforementioned finite element model module includes a geometric modeling module, a mesh generation module, a property setting module, and a shape function setting module. The geometric modeling module is used to establish the geometric model of the reactor; the mesh generation module is used to divide the geometric model into finite element mesh elements; the property setting module is used to define the material properties, boundary conditions, and load conditions of the finite element mesh elements; and the shape function setting module is used to select the element type and shape function of the finite element mesh elements to approximately represent the displacement field of the finite element mesh elements.

[0040] After preprocessing the numerical calculation model of the reactor using the finite element model module, the solution engine module is invoked, and the system will automatically perform nonlinear solutions and material state variable update calculations.

[0041] In some implementations, see Figure 2 As shown, the solver module calls the nonlinear solver to perform nonlinear iterative solving, including the following steps.

[0042] S1, Initialize displacement increment After invoking the solver module, incremental step calculations begin, initializing the incremental data of the state variables to 0. S2, Calculation of strain increment Δ based on displacement increment ε The strain increment can be calculated based on the relationship between strain and stress.

[0043] S3 inputs the strain increment into the material model module for strain correction, resulting in updated stress and consistent tangent modulus.

[0044] S4 receives the stress and uniform tangent modulus from the material model module via the finite element numerical discretization interface, and assembles the residual vector based on the stress and uniform tangent modulus. And Jacobi matrix Update Newton's iterative formula .

[0045] S5, call the nonlinear solver to solve the Newton iteration formula from step S4, and obtain the updated displacement increment. Repeat steps S2 to S5 until the residuals converge.

[0046] In some implementations, the material model module includes a radial regression calculation module for ductile materials and a radial regression calculation module for creep materials. The radial regression calculation module for ductile materials is used to update the stress and uniform tangent modulus of ductile materials; the radial regression calculation module for creep materials is used to update the stress and uniform tangent modulus of creep materials.

[0047] Furthermore, the material model module also includes a material type determination module, which is used to determine the material type based on material properties; when the material is a plastic material, the plastic material radial regression calculation module is called to update the stress and uniform tangent modulus; when the material is a creep material, the creep material radial regression calculation module is called to update the stress and uniform tangent modulus.

[0048] Specifically, for ductile materials, stress and uniform tangent modulus are updated using a radial regression algorithm for ductile materials.

[0049] When a material is subjected to an external force exceeding its elastic limit, the permanent deformation that occurs is plastic deformation. Plastic deformation does not completely return to its original state after the external force is removed. The stress in the residual vector and Jacobian matrix is ​​calculated using a radial regression algorithm. and consistent tangent modulus The radial return algorithm for plastic materials specifically includes the following calculation steps.

[0050] 1. Elastic Prediction: Assuming that the strain increments in the current increment step are all elastic strains, calculate the test stress:

[0051] 2. Plastic state assessment: Calculate the Von-Mise equivalent test stress:

[0052] in, To test the stress, I It is the identity matrix. This is the equivalent test stress.

[0053] 3. Calculate the yield function:

[0054] in, p For equivalent plastic strain, Let be the hardening function with respect to equivalent plastic strain. This is the yield stress. If the yield function is greater than 0, permanent deformation occurs, and the plastic strain increment needs to be calculated. Otherwise, the assumed strain increment and trial stress are used as the strain and stress under actual conditions.

[0055] 4. Calculation of plastic correction; The equivalent plastic strain increment is solved using the Newton method. The m-th step of the Newton iteration is:

[0056] in, G Shear modulus This represents the equivalent plastic strain increment.

[0057] 5. Update the plastic consistency tangent modulus:

[0058]

[0059] In the formula, It is a second-order unit tensor. It is a 4th-order symmetric unit tensor.

[0060] For creeping materials, stress and uniform tangent modulus are updated using a creeping material radial regression algorithm.

[0061] Creep refers to the phenomenon where the strain of a solid material increases over time while the stress remains constant. It differs from plastic deformation, which typically occurs only after the stress exceeds the elastic limit, while creep can occur even when the stress is below the elastic limit, provided the stress has been applied for a sufficiently long time. The stress in the residual vector and Jacobian matrix is ​​calculated using a radial regression algorithm. and consistent tangent modulus The radial return algorithm for creep materials specifically includes the following calculation steps.

[0062] 1. Elastic prediction; Assuming that the strain increments in the current increment step are all elastic strains, calculate the test stress. .

[0063] 2. Calculate the Von-Mise equivalent test stress. .

[0064] 3. Creep correction calculation; The equivalent creep strain increment is solved using the Newton method. The m-th step of the Newton iteration is:

[0065] in, For the equivalent plastic strain increment, For time step.

[0066] 4. Update creep-consistent tangent modulus:

[0067]

[0068] The basic steps of the radial regression algorithm for creeping materials are the same as those for ductile materials. For creeping materials, the stress state is updated through radial regression to satisfy the creep condition. For ductile materials, a ductile state determination step is added. This involves using the yield function to determine if the material is in a ductile state. If it is not, the assumed strain increment and trial stress are used as the strain and stress under the actual conditions for the next increment step of iterative calculation. If it is in a ductile state, the stress state is updated through the radial regression algorithm to satisfy the yield condition. This process ensures numerical stability and computational efficiency.

[0069] Embodiments of the present invention also provide a small deformation mechanics analysis method applying the above-described small deformation mechanics analysis system based on a reactor numerical computation framework, such as... Figure 3 As shown, the steps include the following.

[0070] S10, Establish the finite element calculation model of the reactor within the reactor numerical calculation framework.

[0071] S20, the solver module is invoked to iteratively solve the finite element model to obtain the displacement field of the reactor; in each iteration: S201, calls the material model module to calculate stress and consistent tangent modulus based on the strain increment calculated in the current iteration; S202, through the finite element numerical discretization interface, the stress and uniform tangent modulus calculated by the material model module are passed into the solution engine module, so that the solution engine module can perform the next iteration calculation based on the stress and uniform tangent modulus.

[0072] The establishment of the finite element calculation model of the reactor includes: Establish a geometric model of the reactor within the framework of reactor numerical computation; Divide the geometric model into finite element meshes; Define the material properties, boundary conditions, and load conditions of the finite element mesh elements; Choose the element type and shape function of the finite element mesh to approximate the displacement field of the finite element mesh.

[0073] In some implementations, the material model module calculates stress and uniform tangent modulus based on the strain increment calculated in the current iteration as follows: the material type is determined based on the material properties; when the material is a plastic material, the stress and uniform tangent modulus are updated using the radial regression algorithm for plastic materials; when the material is a creeping material, the stress and uniform tangent modulus are updated using the radial regression algorithm for creeping materials.

[0074] The radial regression algorithm for plastic materials includes the following steps: S1-1, assuming that the strain increments in the current iteration are all elastic strains, calculate the test stress; S1-2, Calculate the Von-Mise equivalent test stress based on the test stress; S1-3: Based on the equivalent test stress, determine whether the yield function is greater than zero; if it is greater than zero, then correct the plastic strain increment according to S1-4; otherwise, take the assumed strain increment and test stress as the strain and stress under the real conditions. S1-4, based on the equivalent test stress, the equivalent plastic strain increment is solved using the Newton method; S1-5, update the plastic consistency tangent modulus.

[0075] The radial regression algorithm for creep materials includes the following steps: S2-1, Assuming that the strain increments in the current iteration are all elastic strains, calculate the trial stress; S2-2, Calculate the Von-Mise equivalent test stress based on the test stress; S2-3, Based on the equivalent test stress, the equivalent creep strain increment is solved using the Newton method; S2-4, Update creep-consistent tangent modulus.

[0076] This invention designs radial regression algorithms for plastic and creep materials based on a reactor numerical development framework to accurately calculate plastic and creep strains and accelerate computation.

[0077] See Figure 4 As shown, a complete solution process combining the finite element calculation model, the solution engine module, and the material model module is described below.

[0078] (1) Mesh generation: Use mesh generation software to divide the geometric model into finite element elements; (2) The mechanical equations are determined, including material properties, boundary conditions and load conditions; (3) Select the element type and interpolation function, and choose a suitable element type shape function to approximate the displacement field within the element; (4) Start incremental step calculation, and initialize the incremental data of the state variables to 0; (5) Calculate the elastic strain increment and elastic stress during the test; (6) If it is a plastic material, determine whether to perform plastic adjustment by using the yield function. If it is in a non-plastic state, proceed to the next incremental step; if it is in a plastic state, perform plastic adjustment. (7) Calculate plastic strain, uniform tangent modulus and stress using radial regression algorithm; (8) If it is a creeping material, use the radial regression algorithm to calculate the creep strain, and calculate the uniform tangent modulus and stress; (9) Consistent tangent modulus and stress data are fed into the nonlinear solver for solving.

[0079] Embodiments of the present invention also provide an electronic device comprising a processor and a memory, wherein the number of processors may be one or more. The memory, as a computer-readable storage medium, can be used to store software programs, computer-executable programs, and modules. The processor executes various functional applications and data processing of the electronic device by running the software programs, instructions, and modules stored in the memory, thereby realizing the small deformation mechanical analysis method based on a reactor numerical computation framework according to any of the above embodiments of the present invention.

[0080] The memory may primarily comprise a program storage area and a data storage area. The program storage area may store the operating system and at least one application program required for a given function; the data storage area may store data created based on terminal usage. Furthermore, the memory may include high-speed random access memory (RAM) and non-volatile memory, such as at least one disk storage device, flash memory, or other non-volatile solid-state storage device. In some instances, the memory may further include memory remotely located relative to the processor, which can be connected to the electronic device via a network. Examples of such networks include, but are not limited to, the Internet, intranets, local area networks (LANs), mobile communication networks, and combinations thereof.

[0081] Embodiments of the present invention also provide a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the small deformation mechanical analysis method based on a reactor numerical computation framework according to any embodiment of the present invention.

[0082] The computer storage medium of this invention can be any combination of one or more computer-readable media. A computer-readable medium can be a computer-readable signal medium or a computer-readable storage medium. A computer-readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of computer-readable storage media (a non-exhaustive list) include: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, portable compact disk read-only memory (CD-ROM), optical storage device, magnetic storage device, or any suitable combination thereof. In this invention, a computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, apparatus, or device.

[0083] Computer-readable signal media may include data signals propagated in baseband or as part of a carrier wave, carrying computer-readable program code. Such propagated data signals may take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. Computer-readable signal media may also be any computer-readable medium other than computer-readable storage media, capable of sending, propagating, or transmitting programs for use by or in connection with an instruction execution system, apparatus, or device.

[0084] Embodiments of the present invention also provide a computer program product that, when run on a computer, causes the computer to execute the small deformation mechanics analysis method based on the reactor numerical calculation framework of any of the above embodiments of the present invention.

[0085] The specific embodiments described above further illustrate the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above description is only a specific embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A small deformation mechanics analysis system based on a reactor numerical calculation framework, characterized in that, include: The finite element model module is used to build a finite element calculation model of the reactor. The solver module is used to iteratively solve the finite element calculation model based on Newton's iterative formula to obtain the displacement field of the reactor. The material model module is used to calculate stress and uniform tangent modulus using a radial regression algorithm based on the strain increment in each iteration of the solver module, and the calculation follows a unified data interface specification; the unified data interface specification defines the transmission protocol for stress and uniform tangent modulus. The finite element numerical discretization interface is used to input the stress and uniform tangent modulus calculated by the material model module into the solution engine module, so that the solution engine module can perform the next iteration calculation based on the stress and uniform tangent modulus.

2. The small deformation mechanical analysis system based on the reactor numerical calculation framework according to claim 1, characterized in that, The finite element model module includes: The geometric modeling module is used to create the geometric model of the reactor; A mesh generation module is used to divide the geometric model into finite element mesh elements; The attribute setting module is used to define the material properties, boundary conditions, and load conditions of the finite element mesh element. The shape function setting module is used to select the element type and shape function of the finite element mesh element to approximately represent the displacement field of the finite element mesh element.

3. The small deformation mechanics analysis system based on the reactor numerical calculation framework according to claim 1, characterized in that, The solution engine module uses a nonlinear solver to perform nonlinear iterative solving. The solution steps are as follows: S1, Initialize displacement increment; S2, calculate the strain increment based on the displacement increment; S3, The strain increment is input into the material model module for strain correction to obtain updated stress and consistent tangent modulus; S4, Receive the stress and uniform tangent modulus passed in by the material model module through the finite element numerical discretization interface, and assemble the residual vector and Jacobian matrix based on the stress and uniform tangent modulus, and update the Newton iteration formula; S5, call the nonlinear solver to solve the Newton iteration formula to obtain the updated displacement increment; repeat steps S2~S5 until convergence.

4. The small deformation mechanical analysis system based on the reactor numerical calculation framework according to claim 1, characterized in that, The material model module includes a radial regression calculation module for plastic materials and a radial regression calculation module for creep materials; The radial regression calculation module for plastic materials is used to update the stress and uniform tangent modulus of plastic materials; The radial regression calculation module for creep materials is used to update the stress and uniform tangent modulus of creep materials.

5. The small deformation mechanical analysis system based on the reactor numerical calculation framework according to claim 4, characterized in that, The material model module also includes a material type determination module, which is used to determine the material type based on material properties; The material model module is also used to update the stress and uniform tangent modulus by calling the radial regression calculation module for plastic materials when the material is a plastic material; and to update the stress and uniform tangent modulus by calling the radial regression calculation module for creep materials when the material is a creeping material.

6. A small deformation mechanical analysis method based on a reactor numerical calculation framework, characterized in that, The small deformation mechanical analysis system based on the reactor numerical calculation framework as described in any one of claims 1-5, wherein the small deformation mechanical analysis method comprises: Establish a finite element calculation model of the reactor within the numerical calculation framework of the reactor; The solution engine module is invoked to iteratively solve the finite element calculation model to obtain the displacement field of the reactor. In each iteration of the calculation: The material model module is invoked to calculate the stress and uniform tangent modulus based on the strain increment calculated in the current iteration; The stress and uniform tangent modulus calculated by the material model module are input into the solution engine module through the finite element numerical discretization interface, so that the solution engine module can perform the next iteration calculation based on the stress and uniform tangent modulus.

7. The small deformation mechanical analysis method based on the reactor numerical calculation framework according to claim 6, characterized in that, The material model module calculates stress and uniform tangent modulus based on the strain increment calculated in the current iteration as follows: Determine the material type based on its properties; When the material is a plastic material, the radial regression algorithm for plastic materials is used to update the stress and the uniform tangent modulus; When the material is a creeping material, the radial regression algorithm for creeping materials is used to update the stress and the uniform tangent modulus.

8. The small deformation mechanical analysis method based on the reactor numerical calculation framework according to claim 7, characterized in that, The radial regression algorithm for plastic materials includes the following steps: S1-1, assuming that the strain increments in the current iteration are all elastic strains, calculate the test stress; S1-2, Calculate the Von-Mise equivalent test stress based on the test stress; S1-3, Based on the equivalent test stress, determine whether the yield function is greater than zero; if it is greater than zero, then correct the plastic strain increment according to S1-4; otherwise, take the assumed strain increment and the test stress as the strain and stress under the real conditions. S1-4, Based on the equivalent test stress, the equivalent plastic strain increment is solved using the Newton method; S1-5, update the plastic consistency tangent modulus.

9. The small deformation mechanical analysis method based on the reactor numerical calculation framework according to claim 7, characterized in that, The radial regression algorithm for creep materials includes the following steps: S2-1, Assuming that the strain increments in the current iteration are all elastic strains, calculate the trial stress; S2-2, Calculate the Von-Mise equivalent test stress based on the test stress; S2-3, Based on the equivalent test stress, the equivalent creep strain increment is solved using the Newton method; S2-4, Update creep-consistent tangent modulus.

10. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the small deformation mechanical analysis method based on the reactor numerical calculation framework as described in any one of claims 6-9.

11. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the small deformation mechanical analysis method based on the reactor numerical calculation framework as described in any one of claims 6-9.

12. A computer program product comprising a computer program / instructions, characterized in that, When the computer program / instruction is executed by the processor, it implements the small deformation mechanical analysis method based on the reactor numerical calculation framework as described in any one of claims 6-9.

Images