Material multi-section linear plastic hardening model construction and implementation method based on near-field dynamics

Abstract

The invention belongs to the related technical field of engineering mechanics and materials, and discloses a material multi-section linear plastic hardening model construction and implementation method based on near-field dynamics, which is used for effectively solving a nonlinear hardening behavior of a material. The method comprises the following steps: constructing a near-field dynamics initial numerical discrete model; setting boundary conditions needing to be applied to the initial numerical discrete model; calculating a current configuration parameter of the Nth time step; constructing a multi-section linear plastic hardening model parameter matrix; constructing a yield judgment criterion; the yield state is judged on the basis of the multi-section linear plastic hardening model parameter matrix and the yield judgment criterion; calculating equivalent plastic strain increment according to the yield state; and updating a near-field dynamic force density vector and a plastic correlation quantity according to the equivalent plastic strain increment. According to the method, on the premise that the precision is guaranteed, the calculation efficiency of the near-field dynamic elastoplasticity problem is remarkably improved, and calculation obstacles are cleared away for applying high-precision fracture mechanical analysis to a large-scale engineering structure.

Description

Technical Field

[0001] This invention belongs to the technical field of engineering mechanics and materials, and more specifically, relates to a method for constructing and implementing a multi-segment linear plastic hardening model of materials based on peri-field dynamics. Background Technology Structural fracture is one of the leading causes of shipwrecks. Accurate fracture mechanics analysis and fatigue life assessment can predict the strength and lifespan of a structure under various load conditions, thus preventing catastrophic accidents. During the fracture process, a ship undergoes plastic strengthening and crack propagation, ultimately leading to overall fracture. Therefore, numerical simulation of the plastic strengthening and crack propagation processes in metals is crucial. Simulating crack initiation and propagation using numerical methods has always been a highly challenging scientific subject.

[0002] In continuum mechanics, the constitutive relations of materials are typically described using differential equations. However, when a crack tip or crack initiation occurs, the governing equations, expressed in differential form, lose their physical meaning at these locations. This presents numerous challenges for traditional continuum mechanics when solving problems involving discontinuous fields. In recent years, a novel meshless method based on nonlocality, namely peri-field dynamics (PD), has been widely used to simulate crack propagation problems. In this method, matter is discretized into material points, each interacting with material points in its neighboring region. By using integral governing equations, the simulation of discontinuous regions is effectively addressed.

[0003] Traditional nonlinear plastic hardening methods require numerical methods such as Newton's iteration method to solve for the equivalent plastic strain increment in the peri-field dynamics. However, this process involves complex numerical iterations, which significantly reduces the computational efficiency of dynamics (PD) and affects fracture and fatigue calculations at the scale of actual ships. The enormous computational cost of this iterative process severely restricts the application of peri-field dynamics methods in large-scale engineering problems at the scale of actual ships.

[0004] To address the aforementioned issues and achieve rapid plasticity calculation of PD at the scale of actual ships, this invention proposes a method for constructing and implementing a multi-segment linear plasticity hardening model for materials based on near-field dynamics. This method effectively solves the nonlinear hardening behavior of materials by discretizing the nonlinear hardening curve into multiple linear segments, avoiding iterative calculations and significantly improving computational efficiency. Summary of the Invention

[0005] In view of the above-mentioned defects or improvement needs of the existing technology, the present invention provides a method for constructing and implementing a multi-segment linear plastic hardening model of materials based on peri-field dynamics. Its purpose is to quickly solve the nonlinear hardening behavior of materials, thereby solving the technical problem that the existing plastic hardening model requires complex iterations when calculating the equivalent plastic strain increment, resulting in low computational efficiency.

[0006] According to one aspect of the present invention, a method for constructing and implementing a multi-segment linear plastic hardening model of a material based on near-field dynamics is provided, comprising the following steps: S1. Model Construction and Initialization: Based on the material entity structure to be analyzed, an initial numerical discretization model of near-field dynamics is constructed, the structure is discretized into multiple material points, and a near-field domain is defined for each material point; S2. Applying boundary conditions: Based on the external loads and constraints, set corresponding boundary conditions on the initial numerical discrete model. S3. Time stepping and deformation calculation: Under the explicit time integration framework, based on the boundary conditions and the state of the previous time step, calculate the position, velocity and current configuration parameters of all material points in the current N+1 time step, including the near-field bond elongation, volume expansion and cumulative equivalent plastic strain. S4. Construction of a multi-segment linear hardening model: Based on material properties and nonlinear hardening rate, a parameter matrix of a multi-segment linear plastic hardening model is constructed. The construction process includes: obtaining the true stress-plastic strain curve of the material, discretizing the curve from the yield point into several continuous linear segments, determining the tangent modulus of each segment and its intercept on the plastic strain axis, and forming a parameter lookup table. S5. Explicit judgment of yield state: Construct a yield judgment criterion based on the current force vector state bias and cumulative equivalent plastic strain, and determine the current state of each material point as elastic loading / unloading, plastic unloading or plastic loading based on the current configuration parameters obtained in step S3. S6. Plastic Response Branching: Based on the judgment result of step S5, perform branch calculations: S6.1 If it is an elastic loading / unloading or plastic unloading state, then set the equivalent plastic strain increment of the current time step to zero, and calculate the force density vector based on the pure elastic constitutive relation; S6.2 If it is a plastic loading state, calculate the equivalent plastic strain increment of the current time step based on the current configuration parameters, yield judgment criteria and multi-segment linear plastic hardening model parameter matrix. S7. Force state update and system response solution: Based on the equivalent plastic strain increment determined in step S6, update the force density vector and plasticity-related quantities of all near-field bonds; solve the acceleration of all material points by integrating the updated plasticity-related quantities, thereby updating their velocity and position, and completing the calculation of the current time step. S8. Iterative Convergence and Result Output: Repeat steps S3 to S7, advancing the time step until the preset convergence condition or simulation termination condition is met, and output the simulation results including the material's equivalent plastic strain increment, plastic deformation, near-field dynamic elongation, force density vector response, and crack propagation path.

[0007] Preferably, the near-field bond elongation in step S3 is... The calculation formula is as follows:

[0008] Where Y represents the particle coordinates of the current configuration, and X represents the particle coordinates of the reference configuration; Volume expansion θ The calculation formula is as follows:

[0009] in, ω It is a weight function. x The L2 norm represents the reference configuration coordinates. V j Represents the volume of a point mass. j Representing matter points of the same family, N Indicates the number of matter points of the same group; The cumulative equivalent plastic strain is inherited from the previous step.

[0010] Preferably, the expression for the multi-segment linear plastic hardening model in step S4 is as follows:

[0011] in, This indicates the elasticity phase of the model. E For elastic modulus, For elastic strain, This indicates the plastic stage of the model. This represents the intersection point of the nth line segment and the plastic stress axis. The superscript represents the tangent modulus of the nth line segment. p Indicates plasticity, It represents plastic strain.

[0012] Preferably, the yield judgment criterion in step S5 is defined by the following expression:

[0013] in, Represents the set of all possible plastic states. Represents the yield function. Indicates plastic component t d Related association functions, This represents a material function related to cumulative plastic strain. Represents plastic component t d Related European-style spaces; Preferably, in step S6.1, the force density vector The calculation formula is as follows:

[0014] in, These are the near-field dynamic coefficients. For the weight function, This is the judgment value of the elongation elasticity deviation at the (N+1)th time step.

[0015] Preferably, in step S6.2, the calculation formula for the equivalent plastic strain increment is as follows:

[0016] in, For the equivalent plastic strain increment, This is the judgment value for the force density vector. Let be the yield value of the Nth plastic segment. Let be the tangent modulus of the Nth plastic segment. Let be the yield value of the (i+1)th plastic segment. Let be the yield value of the i-th plastic segment. Let be the tangent modulus of the i-th plastic segment. Let be the tangent modulus of the i-th plastic segment. For equivalent plastic strain, This is the tangent modulus of the current plastic segment. , These are the near-field dynamic coefficients.

[0017] Preferably, the specific steps for updating the force density vector in step S7 include: S7.1, Update the cumulative equivalent plastic strain; S7.2. Update the plastic bias portion of the near-field bond elongation according to the correlation flow rule; S7.3 Calculate the force vector state bias component and plasticity-related quantities at the (N+1)th time step. The specific calculation formula is as follows:

[0018] in, The equivalent plastic stress at the current time step, For the equivalent plastic strain increment, This is the judgment value for the force density vector. This represents the increment of the elongation plasticity at step N+1. , These are the near-field dynamic coefficients.

[0019] According to another aspect of the present invention, a material plasticity simulation system based on peri-field dynamics is provided for implementing the above-mentioned method, comprising: a model initialization module for constructing a discrete model and applying boundary conditions; a material constitutive module for providing a multi-segment linear plastic hardening model parameter matrix; an explicit solver module for performing time stepping, yield determination, plasticity increment calculation, and force state update; and a data output module for outputting and visualizing simulation results.

[0020] According to another aspect of the present invention, an engineering structure fracture and fatigue assessment device is provided, comprising the above-mentioned material plasticity simulation system based on peri-field dynamics, for performing elastoplastic fracture simulation on ship and aerospace structures to predict their crack initiation and propagation behavior and fatigue life.

[0021] In summary, compared with the prior art, the multi-segment linear plasticity strengthening constitutive method based on conventional ground-state near-field dynamics provided by this invention has the following beneficial effects: 1. High computational efficiency: By replacing the nonlinear hardening model with a multi-segment linear model, the solution of the equivalent plastic strain increment is transformed from a complex nonlinear iteration into a simple linear calculation, replacing the Newton-Raphson iterative loop. While ensuring accuracy, it significantly improves the computational efficiency of near-field dynamic elastoplastic problems (especially those involving hardening behavior), clearing the computational obstacles for applying high-precision fracture mechanics analysis to large-scale engineering structures.

[0022] 2. Stable algorithm: Explicitly solving the formula avoids convergence problems that may be caused by iteration, thus improving the robustness of numerical computation.

[0023] 3. High versatility: The model construction and implementation method is applicable to various J2 plasticity models based on isotropic strengthening, kinematic strengthening or hybrid strengthening criteria, and is easy to integrate into existing peri-field dynamics solvers.

[0024] 4. Controllable accuracy: By increasing the number of linear segments, the original nonlinear hardening curve can be approximated infinitely, providing a flexible means of balancing computational efficiency and simulation accuracy. Attached Figure Description

[0025] Figure 1 This is the overall algorithm flowchart in an embodiment of the present invention.

[0026] Figure 2 This is a discrete schematic diagram of the discretization of the multi-segment linear plastic hardening model in an embodiment of the present invention.

[0027] Figure 3 This is a comparison chart of the multi-segment linear calculation results and the benchmark data in an embodiment of the present invention. Detailed Implementation

[0028] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention. Furthermore, the technical features involved in the various embodiments of this invention described below can be combined with each other as long as they do not conflict with each other.

[0029] Please see Figure 1 The overall algorithm flowchart is shown in the figure. Specifically, this invention provides a method for constructing and implementing a multi-segment linear plastic hardening model of materials based on near-field dynamics, including the following steps: S1. Model Construction and Initialization: Based on the material entity structure to be analyzed, an initial numerical discretization model of near-field dynamics is constructed, the structure is discretized into multiple material points, and a near-field domain is defined for each material point; Specifically, a geometric model is performed on the engineering structure to be analyzed (such as a notched steel plate). Then, spatial discretization is performed using a perifield dynamics method: the structure is divided into a large number of material points with volume (V) and mass (m). Each material point i is identified by its initial coordinates. A perifield region with radius δ is defined for each material point, and all other material points j within this region are called its "family points". Material points only interact mechanically with their family points within their perifield region through "bonds". The perifield radius δ is closely related to the material properties and scale characteristics of the problem under study, and is typically taken as 3-4 times the average spacing between material points. The displacement, velocity, and acceleration of all material points are initialized to zero, and historical variables (such as equivalent plastic strain) are also zero.

[0030] Each matter point interacts with the matter points in its near-field region. The governing equations for these matter points can be expressed as:

[0031] in, Indicates the density of a substance. Representing a point of matter x At any moment t acceleration, H Represents the near-field region of a matter point. Indicates the material point in the reference configuration x With matter point x′ The relative positions between them Indicates the material points in the current configuration. x With matter point x′ The relative displacement between them This represents the volume of a near-field matter point. This represents the near-field force of a point mass on its own family of particles. This represents the near-field force of a point of the same family on a near-field matter point. Represents the physical force of a material point.

[0032] S2. Applying boundary conditions: Based on the external loads and constraints, set corresponding boundary conditions on the initial numerical discrete model. S3. Time stepping and deformation calculation: Under the explicit time integration framework, based on the boundary conditions and the state of the previous time step, calculate the position, velocity and current configuration parameters of all material points in the current N+1 time step, including the near-field bond elongation, volume expansion and cumulative equivalent plastic strain. Near-field bond elongation in step S3 The calculation formula is as follows:

[0033] Where Y represents the particle coordinates of the current configuration, and X represents the particle coordinates of the reference configuration; Volume expansion θ The calculation formula is as follows:

[0034] in, ω It is a weight function. x The L2 norm represents the reference configuration coordinates. V j Represents the volume of a point mass. j Representing matter points of the same family, N Indicates the number of matter points of the same group; The cumulative equivalent plastic strain is inherited from the previous step.

[0035] S4. Construction of a multi-segment linear hardening model: Based on material properties and nonlinear hardening rate, a parameter matrix of a multi-segment linear plastic hardening model is constructed. The construction process includes: obtaining the true stress-plastic strain curve of the material, discretizing the curve from the yield point into several continuous linear segments, determining the tangent modulus of each segment and its intercept on the plastic strain axis, and forming a parameter lookup table. Specifically, it includes the following steps: S4.1 Obtain the true stress-plastic strain data of the material through uniaxial tensile testing; S4.2, such as Figure 2 As shown, starting from the yield point, its hardening curve is approximated by four straight line segments (K=4); S4.3 Determine the parameters for each line segment; S4.4 Establish a parameter lookup table, as shown in Table 1, to clarify the hardening parameters corresponding to each plastic strain range.

[0036] Table 1 Parameters of the discretized linear segment

[0037] The expression for the multi-segment linear plastic hardening model in S4 is as follows:

[0038] in, This indicates the elasticity phase of the model. E For elastic modulus, For elastic strain, This indicates the plastic stage of the model. This represents the intersection point of the nth line segment and the plastic stress axis. The superscript represents the tangent modulus of the nth line segment. p Indicates plasticity, It represents plastic strain.

[0039] S5. Explicit judgment of yield state: Construct a yield judgment criterion based on the current force vector state bias and cumulative equivalent plastic strain, and determine the current state of each material point as elastic loading / unloading, plastic unloading or plastic loading based on the current configuration parameters obtained in step S3. The yield judgment criterion is defined by the following expression:

[0040] in, Represents the set of all possible plastic states. Represents the yield function. Indicates plastic component t d Related association functions, This represents a material function related to cumulative plastic strain. Represents plastic component t d Related European-style spaces; like ≤0, the material has not yielded and is in an elastic loading and unloading state; like If the value is greater than 0 and the material has already yielded in the previous step, then it is in a state of plastic loading.

[0041] S6. Plastic Response Branching: Based on the judgment result of step S5, perform branch calculations: S6.1 If it is an elastic loading / unloading or plastic unloading state, then set the equivalent plastic strain increment of the current time step to zero, and calculate the force density vector based on the pure elastic constitutive relation; Force density vector The calculation formula is as follows:

[0042] in, Here is the near-field dynamic constant. For the weight function, The elastic deviation value for the elongation at the (N+1)th time step; S6.2 If it is a plastic loading state, calculate the equivalent plastic strain increment of the current time step based on the current configuration parameters, yield judgment criteria and multi-segment linear plastic hardening model parameter matrix. The formula for calculating the equivalent plastic strain increment is as follows:

[0043] in, For the equivalent plastic strain increment, This is the judgment value for the force density vector. Let be the yield value of the Nth plastic segment. Let be the tangent modulus of the Nth plastic segment. Let be the yield value of the (i+1)th plastic segment. Let be the yield value of the i-th plastic segment. Let be the tangent modulus of the i-th plastic segment. Let be the tangent modulus of the i-th plastic segment. For equivalent plastic strain, This is the tangent modulus of the current plastic segment. , These are the near-field dynamic coefficients.

[0044] S7. Force state update and system response solution: Based on the equivalent plastic strain increment determined in step S6, update the force density vector and plasticity-related quantities of all near-field bonds; solve the acceleration of all material points by integrating the updated plasticity-related quantities, thereby updating their velocity and position, and completing the calculation of the current time step. S7.1, Update the cumulative equivalent plastic strain; S7.2. Update the plastic bias portion of the near-field bond elongation according to the correlation flow rule; S7.3 Calculate the force vector state bias component and plasticity-related quantities at the (N+1)th time step; the specific calculation formula is as follows:

[0045] in, The equivalent plastic stress at the current time step, For the equivalent plastic strain increment, This is the judgment value for the force density vector. This represents the increment of the elongation plasticity at step N+1. , These are the near-field dynamic coefficients.

[0046] S8. Iterative Convergence and Result Output: Repeat steps S3 to S7, advancing the time step until the preset convergence condition or simulation termination condition is met, and output the simulation results including the material's equivalent plastic strain increment, plastic deformation, near-field dynamic elongation, force density vector response, and crack propagation path.

[0047] This embodiment also provides a material plasticity simulation system for implementing the above method, including: a model initialization module for constructing a discrete model and applying boundary conditions; a material constitutive module for providing a multi-segment linear plastic hardening model parameter matrix; an explicit solver module for performing time stepping, yield judgment, plasticity increment calculation, and force state update; and a data output module for outputting and visualizing simulation results. Furthermore, it can be encapsulated as an engineering structure fracture and fatigue assessment device, incorporating the aforementioned near-field dynamics-based material plasticity simulation system, for performing elastoplastic fracture simulation of ship and aerospace structures to predict their crack initiation and propagation behavior and fatigue life, for structural safety design and life assessment in fields such as shipbuilding and aviation.

[0048] In this embodiment Figure 1 The differences in computational flow between this method and existing methods are also clarified. The red arrows and solid red boxes represent the computational flow unique to this method; the blue arrows and dashed blue boxes represent the computational flow of existing methods. In this method, only the step of constructing the parameter matrix of the multi-segment linear plastic hardening model is needed when calculating the current configuration parameters; this step is no longer required in subsequent calculations, and the equivalent plastic strain increment for the current step can be directly obtained when calculating the equivalent plastic strain increment. Existing methods, on the other hand, require numerical iterations using methods such as Newton's method or the bisection method when calculating the equivalent plastic strain increment.

[0049] The system in this embodiment can output real-time full-field stress, strain, plastic strain distribution, and damage cloud map. Figure 3 The materials used are two aluminum alloys with different properties; specific parameters are shown in Table 2. Figure 3 The results of the plastic hardening process and nonlinear benchmark solution for two materials calculated using this method are presented. The stress-strain response calculated by the method of this invention agrees well with the high-precision benchmark solution without losing too much numerical accuracy, and can significantly shorten the calculation time. The accurate agreement of the results verifies the accuracy of the method.

[0050] Table 2. Two aluminum alloys with different properties

[0051] Those skilled in the art will readily understand that the above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, and improvements 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 method for constructing and implementing a multi-segment linear plastic hardening model for materials based on peri-field dynamics, characterized in that: Includes the following steps: S1. Model Construction and Initialization: Based on the material entity structure to be analyzed, an initial numerical discretization model of near-field dynamics is constructed, the structure is discretized into multiple material points, and a near-field domain is defined for each material point; S2. Applying boundary conditions: Based on the external loads and constraints, set corresponding boundary conditions on the initial numerical discrete model. S3. Time stepping and deformation calculation: Under the explicit time integration framework, based on the boundary conditions and the state of the previous time step, calculate the position, velocity and current configuration parameters of all material points in the current N+1 time step, including the near-field bond elongation, volume expansion and cumulative equivalent plastic strain. S4. Construction of a multi-segment linear hardening model: Based on material properties and nonlinear hardening rate, a parameter matrix of a multi-segment linear plastic hardening model is constructed. The construction process includes: obtaining the true stress-plastic strain curve of the material, discretizing the curve from the yield point into several continuous linear segments, determining the tangent modulus of each segment and its intercept on the plastic strain axis, and forming a parameter lookup table. S5. Explicit judgment of yield state: Construct a yield judgment criterion based on the current force vector state bias and cumulative equivalent plastic strain, and determine the current state of each material point as elastic loading / unloading, plastic unloading or plastic loading based on the current configuration parameters obtained in step S3. S6. Plastic Response Branching: Based on the judgment result of step S5, perform branch calculations: S6.1 If it is an elastic loading / unloading or plastic unloading state, then set the equivalent plastic strain increment of the current time step to zero, and calculate the force density vector based on the pure elastic constitutive relation; S6.2 If it is a plastic loading state, calculate the equivalent plastic strain increment of the current time step based on the current configuration parameters, yield judgment criteria and multi-segment linear plastic hardening model parameter matrix. S7. Force state update and system response solution: Based on the equivalent plastic strain increment determined in step S6, update the force density vector and plasticity-related quantities of all near-field bonds; solve the acceleration of all material points by integrating the updated plasticity-related quantities, thereby updating their velocity and position, and completing the calculation of the current time step. S8. Iterative Convergence and Result Output: Repeat steps S3 to S7, advancing the time step until the preset convergence condition or simulation termination condition is met, and output the simulation results including material plastic deformation, stress-strain response and crack propagation path.

2. The method for constructing and implementing a multi-segment linear plastic hardening model of materials based on peri-field dynamics as described in claim 1, characterized in that: Near-field bond elongation in step S3 The calculation formula is as follows: Where Y represents the particle coordinates of the current configuration, and X represents the particle coordinates of the reference configuration; Volume expansion θ The calculation formula is as follows: in, ω It is a weight function. x The L2 norm represents the reference configuration coordinates. V j Represents the volume of a point mass. j Representing matter points of the same family, N This indicates the number of matter points of the same family.

3. The method for constructing and implementing a multi-segment linear plastic hardening model of materials based on peri-field dynamics as described in claim 1, characterized in that: The expression for the multi-segment linear plastic hardening model in step S4 is as follows: in, This indicates the elasticity phase of the model. E For elastic modulus, For elastic strain, This indicates the plastic stage of the model. This represents the intersection point of the nth line segment and the plastic stress axis. The superscript represents the tangent modulus of the nth line segment. p Indicates plasticity, It represents plastic strain.

4. The method for constructing and implementing a multi-segment linear plastic hardening model of materials based on peri-field dynamics as described in claim 1, characterized in that: The yield judgment criterion described in step S5 is defined by the following expression: in, Represents the set of all possible plastic states. Represents the yield function. Indicates plastic component t d Related association functions, This represents a material function related to cumulative plastic strain. Represents plastic component t d Related European-style spaces.

5. The method for constructing and implementing a multi-segment linear plastic hardening model of materials based on peri-field dynamics as described in claim 1, characterized in that: In step S6.1, the force density vector The calculation formula is as follows: in, These are the near-field dynamic coefficients. For the weight function, This is the judgment value of the elongation elasticity deviation at the (N+1)th time step.

6. The method for constructing and implementing a multi-segment linear plastic hardening model of materials based on peri-field dynamics as described in claim 1, characterized in that: In step S6.2, the formula for calculating the equivalent plastic strain increment is as follows: in, For the equivalent plastic strain increment, This is the judgment value for the force density vector. Let be the yield value of the Nth plastic segment. Let be the tangent modulus of the Nth plastic segment. Let be the yield value of the (i+1)th plastic segment. Let be the yield value of the i-th plastic segment. Let be the tangent modulus of the i-th plastic segment. Let be the tangent modulus of the i-th plastic segment. For equivalent plastic strain, This is the tangent modulus of the current plastic segment. , These are the near-field dynamic coefficients.

7. The method for constructing and implementing a multi-segment linear plastic hardening model of materials based on peri-field dynamics as described in claim 1, characterized in that: The specific steps for updating the force density vector in step S7 include: S7.1, Update the cumulative equivalent plastic strain; S7.

2. Update the plastic bias portion of the near-field bond elongation according to the correlation flow rule; S7.3 Calculate the force vector state bias component and plasticity-related quantities at the (N+1)th time step; the specific calculation formula is as follows: in, The equivalent plastic stress at the current time step, For the equivalent plastic strain increment, This is the judgment value for the force density vector. This represents the increment of the elongation plasticity at step N+1. , These are the near-field dynamic coefficients.

8. A material plasticity simulation system based on peri-field dynamics, used to implement the method as described in any one of claims 1-7, characterized in that, include: The model initialization module is used to construct the discrete model and apply boundary conditions; the material constitutive module is used to provide the parameter matrix of the multi-segment linear plastic hardening model. The explicit solver module is used to perform time stepping, yield determination, plastic increment calculation, and force state update. The data output module is used to output and visualize simulation results.

9. A device for assessing fracture and fatigue in engineering structures, characterized in that, The material plasticity simulation system based on peri-field dynamics as described in claim 8 is used to perform elastoplastic fracture simulation on ship and aerospace structures to predict their crack initiation and propagation behavior and fatigue life.

Images