Finite element numerical simulation method and device, computer equipment and storage medium

Abstract

The invention relates to a finite element numerical simulation method and device, computer equipment and a storage medium, and belongs to the technical field of material mechanics analys.The method comprises the steps that a geometric model corresponding to a target standard part is constructed based on a main program of finite element analysis software; constructing a finite element model corresponding to the geometric model; solving the second finite element model, determining a state variable corresponding to the current increment step, and sending the state variable to a subprogram of finite element analysis software; calculating a heuristic stress based on the subprogram, determining a correction result of the state variable according to the heuristic stress, and returning the correction result to the main program when the correction result converges; and in response to the detected correction result, utilizing the main program to solve the state variable of the next increment step until the state variables of the preset number of increment steps are solved. According to the invention, the test time and economic cost can be effectively reduced.

Description

Technical Field

[0001] The present application relates to the technical field of material mechanics analysis, and in particular to a finite element numerical simulation method, device, computer equipment and storage medium. Background Art

[0002] When analyzing engineering materials, numerical simulation, experimental methods, theoretical methods, etc. are mainly used. The experimental method is only for simple composite materials, and experimental operations such as stretching, bending, and impact are performed on them. Through a large number of experiments, the experimental data are summarized and relatively reliable empirical formulas are fitted. Theoretical methods are currently mostly for objects with simple structures and regular shapes, and theoretical derivation and calculation of structural responses are performed. Numerical simulation methods are the most practical and common. Finite element software is used to establish material models and input known material parameters to perform discrete analysis on the structure. This method is simple, fast, and cost-saving. ABAQUS (finite element analysis software) is one of the most advanced finite element software in the world. It has a wealth of constitutive models suitable for material analysis, up to hundreds of various unit types, and extraordinary nonlinear analysis and coupled field analysis capabilities, making it a powerful tool in the field of material analysis. However, some commonly used constitutive models in engineering materials are not included in ABAQUS, which limits the application of ABAUQS in material engineering and makes it impossible to quickly conduct further analysis of specific materials.

[0003] Therefore, how to conduct specific analysis on specific materials to effectively reduce time and economic costs is an urgent problem to be solved. Summary of the invention

[0004] Based on this, it is necessary to provide a finite element numerical simulation method, device, computer equipment and storage medium that can effectively reduce time cost and economic cost in response to the above technical problems.

[0005] In one aspect, a finite element numerical simulation method is provided, the method comprising:

[0006] Based on the main program of the finite element analysis software, the geometric model corresponding to the target standard part is constructed;

[0007] Meshing the geometric model according to the target unit type to obtain a first finite element model, setting boundary constraints and load conditions for the first finite element model to obtain a second finite element model;

[0008] Solving the second finite element model, determining the state variables corresponding to the current increment step, and sending the state variables to a subroutine of the finite element analysis software;

[0009] Calculating a trial stress based on the subroutine, determining a correction result of the state variable according to the trial stress, and returning the correction result to the main program when the correction result converges;

[0010] In response to detecting the correction result, the main program is used to solve the state variables of the next incremental step until the state variables of a preset number of incremental steps are solved.

[0011] Optionally, the method for generating the first finite element model includes:

[0012] Acquiring geometric parameters and material parameters of the target standard part;

[0013] Using the main program of the finite element analysis software, an initial geometric model corresponding to the target standard part is constructed according to the geometric parameters and the preset proportions;

[0014] Inputting the material parameters into the initial geometric model to obtain a target geometric model corresponding to the target standard part;

[0015] The target geometric model is meshed according to the target unit type to generate the first finite element model.

[0016] Optionally, the method for generating the second finite element model includes:

[0017] Based on the first preset mapping table, obtaining a target unit type for meshing the target standard part;

[0018] Meshing the target geometric model according to the target unit type to obtain the first finite element model;

[0019] Based on the second preset mapping table, the boundary constraint conditions and the load conditions corresponding to the target standard part are determined, and the boundary constraint conditions and the load conditions are set for the first finite element model to obtain the second finite element model.

[0020] Optionally, before determining the correction result of the state variable according to the test stress, the method further includes:

[0021] Based on the state variables corresponding to the current incremental step, calculate the current moment's trial stress;

[0022] comparing the test stress with the stress at a previous moment;

[0023] In response to detecting that the probe stress is less than the stress, storing the probe stress;

[0024] In response to detecting that the tentative stress is greater than or equal to the stress, the state variable corresponding to the current incremental step is corrected based on the stress.

[0025] Optionally, in response to detecting that the probe stress is greater than or equal to the stress, correcting the state variable based on the stress includes:

[0026] The deviator tensor of the stress is determined, and the equivalent stress model is determined based on the deviator tensor:

[0027] According to the correlation between the equivalent stress value and the equivalent trial stress value, a target scaling factor model is constructed;

[0028] Iteratively solving the target scaling factor model based on the Newton-Raphson algorithm, and updating the state variable before each iterative solution step;

[0029] In response to detecting that the value obtained by iterative solution is less than a preset threshold, solving the equivalent stress model according to the updated state variable to obtain the equivalent stress;

[0030] Based on the equivalent stress, a stress increment is determined, and the state variable corresponding to the current increment step is corrected according to the stress increment.

[0031] Optionally, according to the correlation between the equivalent stress value and the equivalent trial stress value, constructing the target scaling factor model includes:

[0032] In response to detecting that the correlation between the equivalent stress value and the equivalent test stress value is When , constructing the target scaling factor model includes:

[0033]

[0034] Where δ represents the scaling factor, represents the equivalent stress value, represents the equivalent test stress value, G represents the shear modulus, D0 represents the limiting shear strain rate, Z represents the internal state variable of the material, n represents the material strain rate sensitivity constant, and Δt represents the incremental step size;

[0035] Before each iterative solution, the state variable is updated by an update function, and the update function includes:

[0036]

[0037]

[0038]

[0039]

[0040]

[0041]

[0042] Among them, Z I represents the isotropic hardening component, Δt represents the incremental step size, t represents the time, and Z D represents the kinematic hardening component, m1 represents Z I The evolution rate from the initial value Z0 to the final value Z1, represents plastic power, A1 and A2 represent thermal recovery coefficients, r1 and r2 represent thermal recovery indexes, Z2 represents an extreme value, m2 represents the evolution rate of β from the initial value Z0 to the final value Z3, u and V represent the directional tensor of stress σ and the second-order tensor β, respectively. Represents Z I The time derivative, represents the derivative with respect to β time, Z represents the internal state variable of the material, represents the time derivative of T, and T represents temperature.

[0043] Optionally, using the main program to detect whether the calibration result meets a preset standard includes:

[0044] Performing a convergence test on the correction result based on the subroutine;

[0045] In response to detecting that the test result meets a preset standard, determining that the correction result converges;

[0046] In response to detecting that the test result does not meet the preset standard, it is determined that the correction result does not converge.

[0047] On the other hand, a finite element numerical simulation device is provided, the device comprising:

[0048] A geometric model building module is used to build a geometric model corresponding to the target standard part based on the main program of the finite element analysis software;

[0049] A finite element model building module, used for meshing the geometric model according to a target unit type to obtain a first finite element model, setting boundary constraints and load conditions for the first finite element model to obtain a second finite element model;

[0050] A solution module, used for solving the second finite element model, determining the state variables corresponding to the current increment step, and sending the state variables to a subroutine of the finite element analysis software;

[0051] A correction module, used for calculating the trial stress based on the subroutine, determining the correction result of the state variable according to the trial stress, and returning the correction result to the main program when the correction result converges;

[0052] The loop solving module is used to solve the state variables of the next incremental step by using the main program when the correction result is detected, until the state variables of a preset number of incremental steps are solved.

[0053] In another aspect, a computer device is provided, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the following steps are implemented:

[0054] Based on the main program of the finite element analysis software, the geometric model corresponding to the target standard part is constructed;

[0055] Meshing the geometric model according to the target unit type to obtain a first finite element model, setting boundary constraints and load conditions for the first finite element model to obtain a second finite element model;

[0056] Solving the second finite element model, determining the state variables corresponding to the current increment step, and sending the state variables to a subroutine of the finite element analysis software;

[0057] Calculating the trial stress based on the subroutine, determining the correction result of the state variable according to the trial stress, and returning the correction result to the main program when the correction result converges;

[0058] In response to detecting the correction result, the main program is used to solve the state variables of the next incremental step until the state variables of a preset number of incremental steps are solved.

[0059] In another aspect, a computer-readable storage medium is provided, on which a computer program is stored, and when the computer program is executed by a processor, the following steps are implemented:

[0060] Based on the main program of the finite element analysis software, the geometric model corresponding to the target standard part is constructed;

[0061] Meshing the geometric model according to the target unit type to obtain a first finite element model, setting boundary constraints and load conditions for the first finite element model to obtain a second finite element model;

[0062] Solving the second finite element model, determining the state variables corresponding to the current increment step, and sending the state variables to a subroutine of the finite element analysis software;

[0063] Calculating the trial stress based on the subroutine, determining the correction result of the state variable according to the trial stress, and returning the correction result to the main program when the correction result converges;

[0064] In response to detecting the correction result, the main program is used to solve the state variables of the next incremental step until the state variables of a preset number of incremental steps are solved.

[0065] The above-mentioned finite element numerical simulation method, device, computer equipment and storage medium, the method includes: based on the main program of the finite element analysis software, constructing a geometric model corresponding to the target standard part; meshing the geometric model according to the target unit type to obtain a first finite element model, setting boundary constraints and load conditions for the first finite element model to obtain a second finite element model; solving the second finite element model to determine the state variables corresponding to the current incremental step, and sending the state variables to the subroutine of the finite element analysis software; calculating the trial stress based on the subroutine, determining the correction result of the state variable according to the trial stress, and returning the correction result to the main program when the correction result converges; in response to detecting the correction result, using the main program to solve the state variables of the next incremental step until the state variables of a preset number of incremental steps are solved. This application is based on finite element software, and subroutines with different material properties are developed for different materials, which can effectively solve and analyze specific problems, thereby effectively reducing the time and economic cost of the experiment. BRIEF DESCRIPTION OF THE DRAWINGS

[0066] Figure 1 is an application environment diagram of a finite element numerical simulation method in an embodiment;

[0067] Figure 2 A schematic diagram of a finite element numerical simulation method in one embodiment;

[0068] Figure 3 A schematic diagram of the ABAQUS main program flow of a finite element numerical simulation method in an embodiment;

[0069] Figure 4 A schematic diagram of a UMAT subroutine flow chart of a finite element numerical simulation method in one embodiment;

[0070] Figure 5 It is a schematic diagram of convergence calculation results in a constant stiffness method of a finite element numerical simulation method in one embodiment;

[0071] Figure 6 It is a schematic diagram of divergent calculation results in the constant stiffness method of the finite element numerical simulation method in one embodiment;

[0072] Figure 7 is a structural block diagram of a finite element numerical simulation device in one embodiment;

[0073] Figure 8 FIG. 4 is a diagram showing the internal structure of a computer device in one embodiment. DETAILED DESCRIPTION

[0074] In order to make the purpose, technical solutions and advantages of this application clearer, the technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the drawings in the embodiments of this application. Obviously, the described embodiments are only part of the embodiments of this application, not all of them. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0075] It should be understood that in the description of the present application, unless the context clearly requires otherwise, words such as "include", "comprises", and the like throughout the specification should be interpreted as including rather than being exclusive or exhaustive; that is, as including but not limited to.

[0076] It should also be understood that the terms "first", "second", etc. are only used for descriptive purposes and cannot be understood as indicating or implying relative importance. In addition, in the description of this application, unless otherwise specified, "plurality" means two or more.

[0077] It should be noted that the terms "S1", "S2", etc. are only used for the purpose of describing the steps, and do not specifically refer to the order or sequence, nor are they used to limit the present application. They are only for the convenience of describing the method of the present application, and cannot be understood as indicating the order of the steps. In addition, the technical solutions between the various embodiments can be combined with each other, but they must be based on the ability of ordinary technicians in this field to implement them. When the combination of technical solutions is contradictory or cannot be implemented, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection required by the present application.

[0078] The finite element numerical simulation method provided in this application can be applied to Figure 1 In the application environment shown, the terminal 102 communicates with the data processing platform set on the server 104 through the network, wherein the terminal 102 can be but not limited to various personal computers, laptops, smart phones, tablet computers and portable wearable devices, and the server 104 can be implemented by an independent server or a server cluster composed of multiple servers.

[0079] In one embodiment, Figure 2 As shown, a finite element numerical simulation method is provided, which is applied to Figure 1 The terminal in is used as an example to illustrate, including the following steps:

[0080] S1: Based on the main program of the finite element analysis software, the geometric model corresponding to the target standard part is constructed.

[0081] It should be noted that the finite element analysis software used in this application is ABAQUS, which is a powerful finite element software for engineering simulation. The range of problems it can solve ranges from relatively simple linear analysis to many complex nonlinear problems. ABAQUS includes a rich unit library that can simulate any geometric shape, and has various types of material model libraries that can simulate the performance of typical engineering materials, including metals, rubbers, polymer materials, composite materials, reinforced concrete, compressible hyperelastic foam materials, and geological materials such as soil and rock. As a general simulation tool, ABAQUS can not only solve a large number of structural (stress / displacement) problems, but also simulate many problems in other engineering fields, such as heat conduction, mass diffusion, thermoelectric coupling analysis, acoustic analysis, geotechnical mechanics analysis (fluid penetration / stress coupling analysis) and piezoelectric medium analysis. In the finite element analysis software ABAQUS, although there are already a large number of unit material property libraries and solution models, ABAQUS lacks corresponding data and technical support for some special and advanced materials. Therefore, the finite element analysis software in this application includes not only the main program, but also the subprograms (User Defined Material The UMAT subroutine, UMAT, can help users to carry out numerical simulation and analysis for specific materials and specific problems to a great extent. The UMAT subroutine written in FORTRAN (Formula Translation) language can be well connected with the ABAQUS software main program interface to realize data exchange.

[0082] Furthermore, ABAQUS software requires the user to use the Cauchy stress tensor (a mathematical representation of the stress state) in the constitutive model adopted by the user, so as to facilitate the subsequent integral algorithm solution and numerical simulation analysis, and use the Newton-Raphson method to solve the nonlinear finite element equilibrium equation:

[0083] ∑(∫B T JB k =P(t+Δt)-∑∫B T [σ(t)+Δσ(Δu k )]dV (4.1)

[0084] Δu k+1 =Δu k +C k (4.2)

[0085] Among them, B is the unit geometry matrix, and the node displacement increment can be converted into strain increment through the unit displacement matrix:

[0086] Δε=B·Δu (4.3)

[0087] Among them, the stress increment is related to the strain increment. From formula (4.3), it can be seen that the stress increment is also related to the displacement increment. When using the Newton-Raphson iterative algorithm to solve the equilibrium equation of each load step, the Jacobian matrix needs to be given:

[0088]

[0089] It is represented as the DDSDDE matrix in the UMAT subroutine. Usually, the DDSDDE matrix only affects the calculation convergence speed but not the calculation results. Therefore, considering the calculation cost and program simplification, the elastic stiffness matrix is ​​used as the Jacobian matrix in solving the nonlinear finite element equilibrium equation, that is, the "constant stiffness method".

[0090] Furthermore, when the constant stiffness method is used to solve nonlinear finite element problems, its essence is to solve a nonlinear system of equations. The overall degrees of freedom of the finite element problem can be expressed as:

[0091]

[0092] In the formula, is the basic unknown vector, for linear problems, is a first-order term, and H is It is irrelevant, equation (4.5) is a linear problem, which can be easily solved to obtain For nonlinear problems, it is no longer a simple linear term, but a multi-term related to . Equation (4.5) is a nonlinear problem. In the constant stiffness method, the nth step is usually used. Solve the H of step n+1, and then use equation (4.6) to solve the H of step n+1

[0093]

[0094] From equation (4.6), we can see that an initial value needs to be given when starting to solve, and for nonlinear problems, the stiffness matrix needs to be inverted. There will be "convergence" and "divergence" in the solution process. Incorrect solution methods can easily lead to "divergence" of calculation results, such as Figure 5-Figure 6 As shown in the figure, the criterion for determining the completion of the solution is that the solution process must be "converged", which is reflected in the iterative process as the difference between the unknown quantity results obtained in the two iterations before and after meets the given tolerance, indicating that the iteration is completed and the solution is finished.

[0095] In some specific embodiments, Figure 3 As shown, the method for generating the first finite element model includes:

[0096] Obtaining geometric parameters and material parameters of the target standard parts, wherein the target standard parts refer to mechanical parts that have been fully standardized, such as fasteners, connectors, transmission parts, seals, hydraulic components, etc. The geometric parameters may be parameters such as the circumference, area, and diameter of the target standard parts, and the material parameters refer to key indicators that describe the physical, chemical, and mechanical properties of the materials, such as density, hardness, ductility, strength, thermal stability, heat resistance, corrosion resistance, thermal conductivity, etc.;

[0097] Using the main program of the finite element analysis software, an initial geometric model corresponding to the target standard part is constructed according to the geometric parameters and the preset proportion, wherein the preset proportion can be set according to actual needs, such as constructing the initial geometric model according to the actual proportion;

[0098] Inputting the material parameters into the initial geometric model, that is, adding material properties to the initial geometric model, to obtain a target geometric model corresponding to the target standard part;

[0099] The target geometric model is meshed according to the target unit type to generate the first finite element model, wherein the finite element model is a solid model that divides a continuous geometric model into a finite number of units through a finite element method, and the division method is a meshing method.

[0100] S2: Meshing the geometric model according to the target unit type to obtain a first finite element model, setting boundary constraints and load conditions for the first finite element model to obtain a second finite element model.

[0101] It should be noted that the unit types may include bar unit, beam unit, shell unit, plate unit, plane stress unit, plane strain unit, etc.

[0102] In some specific embodiments, the method for generating the second finite element model includes:

[0103] Based on the first preset mapping table, obtaining a target unit type for meshing the target standard part;

[0104] Meshing the target geometric model according to the target unit type to obtain the first finite element model;

[0105] Based on a second preset mapping table, boundary constraints and load conditions corresponding to the target standard part are determined, and the boundary constraints and load conditions are set for the first finite element model to obtain the second finite element model, wherein the boundary constraints may include inertial loads, structural loads, structural constraints, and thermal loads, and the load conditions are to apply an initial load to the finite element model.

[0106] Among them, the first preset mapping table at least includes a mapping relationship formed by a standard part and a unit type, and the second preset mapping table at least includes a mapping relationship formed by a standard part and a boundary constraint condition and a load condition. The values ​​in the first preset mapping table and the second preset mapping table are obtained based on expert experience. For example, if the standard part is a connecting part, if its commonly used unit type is a rod unit, then the standard part forms a mapping relationship with the rod unit. If its boundary constraint condition is a thermal load and the initial load is X, then the standard part forms a mapping relationship with the thermal load and X. After determining the second finite element model through the preset mapping table, the staff can also adjust it according to actual conditions to ensure the accuracy of the finite element model construction.

[0107] S3: Solving the second finite element model, determining the state variables corresponding to the current incremental step, and sending the state variables to a subroutine of the finite element analysis software.

[0108] It should be noted that the solution result of the second finite element model can be obtained through the specific steps of step S1. In each incremental step solution process, the finite element model is solved mainly by the ABAQUS main program according to boundary constraints, load application conditions, meshing strategies, and analysis step settings, and the strain increment, time step, and various state variable group values ​​of each incremental step are calculated. When calling the UMAT subroutine, the ABAQUS main program passes the above variables to it.

[0109] S4: Calculating the trial stress based on the subroutine, determining the correction result of the state variable according to the trial stress, and returning the correction result to the main program when the correction result converges.

[0110] It should be noted that the trial stress can be obtained by calculating the input variables. In this application, the stress is mainly corrected. The subroutine mainly includes the constructed viscoplastic constitutive equation and solving it by adapting the Euler implicit integration algorithm. The stress increment can be obtained to correct the variable value obtained by the main program.

[0111] Specifically, starting from the generalized form, the model of the viscoplastic unified constitutive theory can be expressed as:

[0112]

[0113] Among them, Y is the solution vector, which contains stress tensor, strain tensor or internal variable tensor, and is a function of time. The initial condition is Y(t0) = Y0. The rate form constitutive equation is rewritten as an incremental form formula, and it can be easily solved by using Euler implicit integration, Euler explicit integration or mixed integration. Equation (4.7) is rewritten as an incremental form:

[0114] Yt+Δt =Y t +ΔY(Δt) (4.8)

[0115] Among them, Y t represents the solution vector at time t (previous time), ΔY represents the change in the solution vector caused by the time increment, and Y t+Δt is the solution vector at time t+Δt (the next moment).

[0116] The main task of the constitutive equation integration algorithm is not only to calculate the results, but also to ensure the stability and correctness of the calculation process. The viscoplastic unified constitutive theory has a "stiffness" effect in mathematical form, that is, the function form used in the constitutive theory causes the relevant variables in the constitutive differential equation to be sensitive to changes in independent variables or time increments. The intuitive feeling is that the significant increase in the inelastic strain rate during the numerical simulation loading process makes the equation integration of the viscoplastic unified constitutive theory unstable or even incorrect. Therefore, the adaptive Euler implicit integration algorithm is selected to ensure the premise of accurate calculation results. Specifically:

[0117] The method of numerical solution according to implicit formula is the adaptive backward Euler implicit integration algorithm. Usually, the implicit formula cannot be used directly. It is necessary to obtain the initial value through the explicit formula first, and then use the implicit formula to solve. When dealing with linear problems, the implicit integration algorithm has absolute stability. When dealing with nonlinear problems, it is necessary to solve through linear approximation. Among them, Newton-Raphson is the most widely used iterative algorithm.

[0118] Considering the equivalent inelastic strain rate:

[0119]

[0120] Under uniaxial stress state, in, Denotes the equivalent stress, rewrite equation (4.9) to obtain the inelastic strain rate expression and equivalent inelastic strain rate expression under uniaxial stress state:

[0121]

[0122]

[0123] Assume that the variable at time t is σ t , ε t , and Z t It is known that the strain increment Δε within Δt is calculated according to the load conditions t+Δt Based on elastic prediction (stress trial) and plastic correction combined with Newton-Raphson iterative algorithm, the stress increment Δσ t+Δt Solving can update the stresses.

[0124] First, assuming that all strain increments are elastic, the stress update formula is:

[0125]

[0126] Compare the test stress at the next moment The stress at the previous moment like This means that the material undergoes elastic deformation and can be directly used Updated as final stress if This indicates that the material has plastic flow and needs plastic correction:

[0127]

[0128] Under high pressure, the average normal stress only causes elastic volume deformation without plastic volume deformation, while the stress deviator will produce distortion without volume deformation. For isotropic materials, they have the same elastic properties in all directions. Considering the relationship between the stress deviator increment and the strain deviator increment in the elastic stage, the deviator tensor form of equation (4.13) is:

[0129]

[0130] Based on the associated plastic flow law (Prandtl-Reuss flow law), (4.14) is rewritten as:

[0131]

[0132] During the plastic correction process, the radial return direction always remains in the same direction, so there is Combining (4.15) we have:

[0133]

[0134] When solving the equivalent plastic strain, the convergence will be affected by the numerical value. Therefore, a scaling factor δ (0 < δ < 1) is considered. Assume that the equivalent stress at time t + Δt has the following relationship with the equivalent test: Construct the δ equation, that is, the target scaling factor model:

[0135]

[0136] Where δ represents the scaling factor, represents the equivalent stress value, represents the equivalent test stress value, G represents the shear modulus, D0 represents the limiting shear strain rate, Z represents the internal state variable of the material, n represents the material strain rate sensitivity constant, and Δt represents the incremental step size;

[0137] According to the Newton-Raphson algorithm, equation (4.17) is iteratively solved starting from δ = 1 until F(δ) is less than a certain value, which means that the solution δ meets the requirements. The Newton-Raphson iterative algorithm process is:

[0138]

[0139] Before each iterative solution, the internal variables are updated according to the update function, which is:

[0140]

[0141]

[0142]

[0143]

[0144]

[0145]

[0146] Among them, Z I represents the isotropic hardening component, Δt represents the incremental step size, t represents the time, and Z D represents the kinematic hardening component, m1 represents Z I The evolution rate from the initial value Z0 to the final value Z1, represents plastic power, A1 and A2 represent thermal recovery coefficients, r1 and r2 represent thermal recovery indexes, Z2 represents an extreme value, m2 represents the evolution rate of β from the initial value Z0 to the final value Z3, u and V represent the directional tensor of stress σ and the second-order tensor β, respectively. Represents Z I The time derivative, represents the derivative with respect to β time, Z represents the internal state variable of the material, represents the derivative of T with respect to time, T represents temperature,

[0147] After determining that the required δ is obtained, the equivalent stress at time t+Δt is obtained. Solving equations (4.9), (4.10) and (4.11) together gives the plastic strain increment: and stress increment Δσ t+Δt , and update the stresses while updating other variables:

[0148] σ t+Δt =σ t +Δσ t+Δt (4.22)

[0149] Among them, the thermomechanical fatigue of metal materials usually includes two load modes, one is the same direction positive phase (same phase) in which the load and temperature maintain the same change trend, and the other is the reverse phase (anti-phase) in which the load and temperature maintain opposite change trends. Therefore, in the thermomechanical fatigue cyclic plastic constitutive model of metal materials considering the temperature history effect, the sensitivity of the internal variables to temperature changes is as well as To characterize the changes in the internal state variables of the material when it is subjected to mechanical loads and temperature loads, and to affect the final results;

[0150] Furthermore, the intrinsic variables consist of an isotropic hardening component and a kinematic hardening component:

[0151] Z=Z I +Z D (4.23)

[0152] Z I The evolution equation under proportional loading conditions is:

[0153]

[0154] The first term on the right side of the equation represents the work hardening term, and the second term represents the thermal recovery term;

[0155] Z D The evolution equation under proportional loading conditions is expressed by the inner product of the second-order tensor β and the stress direction tensor u:

[0156] Z D =β:u (4.25)

[0157]

[0158] The first term on the right side of the equation represents the work hardening term, and the second term represents the thermal recovery term;

[0159]

[0160] Based on the temperature history effect, that is, the influence of temperature changes on the thermomechanical fatigue of metal materials, terms related to temperature changes are added to the isotropic hardening component and the kinematic hardening component. Therefore, Equations (4.24) and (4.26) are rewritten in the form of Equations (4.20) and (4.21), and the third term on the right side of the equation represents the temperature change term.

[0161] In some specific embodiments, Figure 4 As shown, before determining the correction result of the state variable according to the test stress, the method further includes:

[0162] Based on the state variables corresponding to the current increment, the trial stress at the current moment is calculated. The state variables are the internal variables passed into the main program, such as strain ε, inelastic strain ε in , strain increment Δε, stress σ and other stored internal variables, calculate the trial stress according to the above steps

[0163] The test stress is compared with the stress at the previous moment, that is, With σ t Make comparisons;

[0164] In response to detecting that the test stress is less than the stress, the test stress is saved, that is, When , it means that the material does not undergo plastic flow and directly Perform stress updates;

[0165] In response to detecting that the test stress is greater than or equal to the stress, the state variable corresponding to the current incremental step is corrected based on the stress, that is, It indicates that the material undergoes plastic flow and needs plastic correction.

[0166] In some specific embodiments, in response to detecting that the probe stress is greater than or equal to the stress, correcting the state variable based on the stress includes:

[0167] The deviator tensor of the stress is determined, and the equivalent stress model is determined based on the deviator tensor:

[0168] According to the correlation between the equivalent stress value and the equivalent trial stress value, a target scaling factor model is constructed;

[0169] Iteratively solving the target scaling factor model based on the Newton-Raphson algorithm, and updating the state variable before each iterative solution step;

[0170] In response to detecting that the value obtained by iterative solution is less than a preset threshold, the equivalent stress model is solved according to the updated state variable to obtain the equivalent stress, wherein the preset threshold can be set according to actual needs, such as 0.01, that is, when F(δ) < 0.001, it means that the solved δ meets the requirements;

[0171] Based on the equivalent stress, a stress increment is determined, and the state variable corresponding to the current increment step is corrected according to the stress increment.

[0172] In some specific implementations, according to the correlation between the equivalent stress value and the equivalent trial stress value, constructing the target scaling factor model includes:

[0173] In response to detecting that the correlation between the equivalent stress value and the equivalent test stress value is When , construct a model about the target scaling factor;

[0174] Before each iterative solution step, the state variable is updated by an update function;

[0175] Specifically, based on the above-mentioned metal plastic flow law and introducing the reduction factor δ, a nonlinear scalar equation about (i.e., the target scaling factor model) is constructed, and δ is solved by the Newton-Raphson iterative algorithm. After the error criterion is judged (i.e., the value obtained by the iterative solution is less than the preset threshold), the internal variables are updated with those that meet the error requirements, and the predicted response to the temperature change is achieved by updating the internal variables considering the temperature effect.

[0176] In some specific implementations, using the subroutine to detect whether the correction result converges includes:

[0177] Performing a convergence test on the correction result based on the subroutine;

[0178] In response to detecting that the test result meets a preset standard, determining that the correction result converges;

[0179] In response to detecting that the test result does not meet the preset standard, it is determined that the correction result does not converge.

[0180] The correction result detection is to detect whether the correction result converges in the solution process through a convergence test. The preset standard can be as follows: Figure 5-Figure 6 to define.

[0181] S5: In response to detecting the correction result, the main program is used to solve the state variables of the next incremental step until the state variables of a preset number of incremental steps are solved.

[0182] It should be noted that after the plastic correction is completed, the stored stress and other state variables as well as the DDSDDE matrix are updated, and the corrected state variables are sent to the main program. When the main program receives the state variables, it determines that the calculation of the current incremental step is completed and enters the calculation of the next incremental step until the calculation of all incremental steps is completed.

[0183] In the above-mentioned finite element numerical simulation method, the method includes: constructing a geometric model corresponding to the target standard part based on the main program of the finite element analysis software; meshing the geometric model according to the target unit type to obtain a first finite element model, setting boundary constraints and load conditions for the first finite element model to obtain a second finite element model; solving the second finite element model to determine the state variables corresponding to the current incremental step, and sending the state variables to a subroutine of the finite element analysis software; calculating the trial stress based on the subroutine, determining the correction result of the state variable according to the trial stress, and returning the correction result to the main program ; When the main program detects that the correction result meets the preset standard, the state variables of the next incremental step are solved until the state variables of the preset number of incremental steps are solved. This application can solve practical engineering problems well by effectively combining the constitutive model with the numerical analysis method and based on the finite element analysis technology. Through the finite element software, subroutines with different material properties are developed for different materials, which can effectively solve and analyze specific problems. According to this numerical simulation method, the isothermal and thermomechanical fatigue mechanical behavior of the material can be described on the basis of the cyclic constitutive model, which can effectively reduce the time and economic cost of the test.

[0184] It should be understood that although Figure 2-Figure 4 The steps in the flowchart are shown in sequence as indicated by the arrows, but these steps are not necessarily executed in the order indicated by the arrows. Unless otherwise specified in this document, there is no strict order restriction for the execution of these steps, and these steps can be executed in other orders. Moreover, Figure 2-Figure 4 At least part of the steps may include multiple sub-steps or multiple stages. These sub-steps or stages are not necessarily executed at the same time, but can be executed at different times. The execution order of these sub-steps or stages is not necessarily sequential, but can be executed in turn or alternately with other steps or at least part of the sub-steps or stages of other steps.

[0185] In one embodiment, Figure 7 As shown, a finite element numerical simulation device is provided, comprising: a geometric model building module, a finite element model building module, a solution module, a correction module and a loop solution module, wherein:

[0186] A geometric model building module is used to build a geometric model corresponding to the target standard part based on the main program of the finite element analysis software;

[0187] A finite element model building module, used for meshing the geometric model according to a target unit type to obtain a first finite element model, setting boundary constraints and load conditions for the first finite element model to obtain a second finite element model;

[0188] A solution module, used for solving the second finite element model, determining the state variables corresponding to the current increment step, and sending the state variables to a subroutine of the finite element analysis software;

[0189] A correction module, used for calculating the trial stress based on the subroutine, determining the correction result of the state variable according to the trial stress, and returning the correction result to the main program when the correction result converges;

[0190] The loop solving module is used to solve the state variables of the next incremental step by using the main program when the correction result is detected, until the state variables of a preset number of incremental steps are solved.

[0191] As a preferred implementation, in the embodiment of the present invention, the finite element model building module is specifically used for:

[0192] Acquiring geometric parameters and material parameters of the target standard part;

[0193] Using the main program of the finite element analysis software, an initial geometric model corresponding to the target standard part is constructed according to the geometric parameters and the preset proportions;

[0194] Inputting the material parameters into the initial geometric model to obtain a target geometric model corresponding to the target standard part;

[0195] The target geometric model is meshed according to the target unit type to generate the first finite element model.

[0196] As a preferred implementation, in the embodiment of the present invention, the finite element model building module is further used for:

[0197] Based on the first preset mapping table, obtaining a target unit type for meshing the target standard part;

[0198] Meshing the target geometric model according to the target unit type to obtain the first finite element model;

[0199] Based on the second preset mapping table, the boundary constraint conditions and the load conditions corresponding to the target standard part are determined, and the boundary constraint conditions and the load conditions are set for the first finite element model to obtain the second finite element model.

[0200] As a preferred implementation, in the embodiment of the present invention, the correction module is specifically used for:

[0201] Based on the state variables corresponding to the current incremental step, calculate the current moment's trial stress;

[0202] comparing the test stress with the stress at a previous moment;

[0203] In response to detecting that the probe stress is less than the stress, storing the probe stress;

[0204] In response to detecting that the tentative stress is greater than or equal to the stress, the state variable corresponding to the current incremental step is corrected based on the stress.

[0205] As a preferred implementation, in the embodiment of the present invention, the correction module is further used for:

[0206] The deviator tensor of the stress is determined, and the equivalent stress model is determined based on the deviator tensor:

[0207] According to the correlation between the equivalent stress value and the equivalent trial stress value, a target scaling factor model is constructed;

[0208] Iteratively solving the target scaling factor model based on the Newton-Raphson algorithm, and updating the state variable before each iterative solution step;

[0209] In response to detecting that the value obtained by iterative solution is less than a preset threshold, solving the equivalent stress model according to the updated state variable to obtain the equivalent stress;

[0210] Based on the equivalent stress, a stress increment is determined, and the state variable corresponding to the current increment step is corrected according to the stress increment.

[0211] As a preferred implementation, in the embodiment of the present invention, the correction module is further used for:

[0212] In response to detecting that the correlation between the equivalent stress value and the equivalent test stress value is When , constructing the target scaling factor model includes:

[0213]

[0214] Where δ represents the scaling factor, represents the equivalent stress value, represents the equivalent test stress value, G represents the shear modulus, D0 represents the limiting shear strain rate, Z represents the internal state variable of the material, n represents the material strain rate sensitivity constant, and Δt represents the incremental step size;

[0215] Before each iterative solution, the state variable is updated by an update function, and the update function includes:

[0216]

[0217]

[0218]

[0219]

[0220]

[0221]

[0222] Among them, Z I represents the isotropic hardening component, Δt represents the incremental step size, t represents the time, and Z D represents the kinematic hardening component, m1 represents Z I The evolution rate from the initial value Z0 to the final value Z1, represents plastic power, A1 and A2 represent thermal recovery coefficients, r1 and r2 represent thermal recovery indexes, Z2 represents an extreme value, m2 represents the evolution rate of β from the initial value Z0 to the final value Z3, u and V represent the directional tensor of stress σ and the second-order tensor β, respectively. Represents Z I The time derivative, represents the derivative with respect to β time, Z represents the internal state variable of the material, represents the time derivative of T, and T represents temperature.

[0223] As a preferred implementation, in the embodiment of the present invention, the correction module is specifically used for:

[0224] Performing a convergence test on the correction result based on the subroutine;

[0225] In response to detecting that the test result meets a preset standard, determining that the correction result converges;

[0226] In response to detecting that the test result does not meet the preset standard, it is determined that the correction result does not converge.

[0227] For the specific definition of the finite element numerical simulation device, please refer to the definition of the finite element numerical simulation method above, which will not be repeated here. Each module in the above-mentioned finite element numerical simulation device can be implemented in whole or in part by software, hardware and a combination thereof. The above-mentioned modules can be embedded in or independent of the processor in the computer device in the form of hardware, or can be stored in the memory of the computer device in the form of software, so that the processor can call and execute the operations corresponding to the above modules.

[0228] In one embodiment, a computer device is provided. The computer device may be a terminal, and its internal structure diagram may be as follows: Figure 8As shown. The computer device includes a processor, a memory, a network interface, a display screen and an input device connected through a system bus. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The network interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, a finite element numerical simulation method is implemented. The display screen of the computer device can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer device can be a touch layer covering the display screen, or a button, trackball or touchpad set on the computer device housing, or an external keyboard, touchpad or mouse, etc.

[0229] Those skilled in the art will understand that Figure 8 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0230] In one embodiment, a computer device is provided, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the following steps when executing the computer program:

[0231] S1: Based on the main program of the finite element analysis software, the geometric model corresponding to the target standard part is constructed;

[0232] S2: Meshing the geometric model according to the target unit type to obtain a first finite element model, setting boundary constraints and load conditions for the first finite element model to obtain a second finite element model;

[0233] S3: solving the second finite element model, determining the state variables corresponding to the current incremental step, and sending the state variables to a subroutine of the finite element analysis software;

[0234] S4: calculating the trial stress based on the subroutine, determining the correction result of the state variable according to the trial stress, and returning the correction result to the main program when the correction result converges;

[0235] S5: In response to detecting the correction result, the state variables of the next incremental step are solved using the main program until the state variables of a preset number of incremental steps are solved.

[0236] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0237] Acquiring geometric parameters and material parameters of the target standard part;

[0238] Using the main program of the finite element analysis software, an initial geometric model corresponding to the target standard part is constructed according to the geometric parameters and the preset proportions;

[0239] Inputting the material parameters into the initial geometric model to obtain a target geometric model corresponding to the target standard part;

[0240] The target geometric model is meshed according to the target unit type to generate the first finite element model.

[0241] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0242] Based on the first preset mapping table, obtaining a target unit type for meshing the target standard part;

[0243] Meshing the target geometric model according to the target unit type to obtain the first finite element model;

[0244] Based on the second preset mapping table, the boundary constraint conditions and the load conditions corresponding to the target standard part are determined, and the boundary constraint conditions and the load conditions are set for the first finite element model to obtain the second finite element model.

[0245] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0246] Based on the state variables corresponding to the current incremental step, calculate the current moment's trial stress;

[0247] comparing the test stress with the stress at a previous moment;

[0248] In response to detecting that the probe stress is less than the stress, storing the probe stress;

[0249] In response to detecting that the tentative stress is greater than or equal to the stress, the state variable corresponding to the current incremental step is corrected based on the stress.

[0250] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0251] The deviator tensor of the stress is determined, and the equivalent stress model is determined based on the deviator tensor:

[0252] According to the correlation between the equivalent stress value and the equivalent trial stress value, a target scaling factor model is constructed;

[0253] Iteratively solving the target scaling factor model based on the Newton-Raphson algorithm, and updating the state variable before each iterative solution step;

[0254] In response to detecting that the value obtained by iterative solution is less than a preset threshold, solving the equivalent stress model according to the updated state variable to obtain the equivalent stress;

[0255] Based on the equivalent stress, a stress increment is determined, and the state variable corresponding to the current increment step is corrected according to the stress increment.

[0256] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0257] In response to detecting that the correlation between the equivalent stress value and the equivalent test stress value is When , constructing the target scaling factor model includes:

[0258]

[0259] Where δ represents the scaling factor, represents the equivalent stress value, represents the equivalent test stress value, G represents the shear modulus, D0 represents the limiting shear strain rate, Z represents the internal state variable of the material, n represents the material strain rate sensitivity constant, and Δt represents the incremental step size;

[0260] Before each iterative solution, the state variable is updated by an update function, and the update function includes:

[0261]

[0262]

[0263]

[0264]

[0265]

[0266]

[0267] Among them, Z I represents the isotropic hardening component, Δt represents the incremental step size, t represents the time, and Z D represents the kinematic hardening component, m1 represents Z I The evolution rate from the initial value Z0 to the final value Z1, represents plastic power, A1 and A2 represent thermal recovery coefficients, r1 and r2 represent thermal recovery indexes, Z2 represents an extreme value, m2 represents the evolution rate of β from the initial value Z0 to the final value Z3, u and V represent the directional tensor of stress σ and the second-order tensor β, respectively. Represents Z I The time derivative, represents the derivative with respect to β time, Z represents the internal state variable of the material, represents the time derivative of T, and T represents temperature.

[0268] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0269] Performing a convergence test on the correction result based on the subroutine;

[0270] In response to detecting that the test result meets a preset standard, determining that the correction result converges;

[0271] In response to detecting that the test result does not meet the preset standard, it is determined that the correction result does not converge.

[0272] In one embodiment, a computer readable storage medium is provided, on which a computer program is stored, and when the computer program is executed by a processor, the following steps are implemented:

[0273] S1: Based on the main program of the finite element analysis software, the geometric model corresponding to the target standard part is constructed;

[0274] S2: Meshing the geometric model according to the target unit type to obtain a first finite element model, setting boundary constraints and load conditions for the first finite element model to obtain a second finite element model;

[0275] S3: solving the second finite element model, determining the state variables corresponding to the current incremental step, and sending the state variables to a subroutine of the finite element analysis software;

[0276] S4: calculating the trial stress based on the subroutine, determining the correction result of the state variable according to the trial stress, and returning the correction result to the main program when the correction result converges;

[0277] S5: In response to detecting the correction result, the state variables of the next incremental step are solved using the main program until the state variables of a preset number of incremental steps are solved.

[0278] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0279] Acquiring geometric parameters and material parameters of the target standard part;

[0280] Using the main program of the finite element analysis software, an initial geometric model corresponding to the target standard part is constructed according to the geometric parameters and the preset proportions;

[0281] Inputting the material parameters into the initial geometric model to obtain a target geometric model corresponding to the target standard part;

[0282] The target geometric model is meshed according to the target unit type to generate the first finite element model.

[0283] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0284] Based on the first preset mapping table, obtaining a target unit type for meshing the target standard part;

[0285] Meshing the target geometric model according to the target unit type to obtain the first finite element model;

[0286] Based on the second preset mapping table, the boundary constraint conditions and the load conditions corresponding to the target standard part are determined, and the boundary constraint conditions and the load conditions are set for the first finite element model to obtain the second finite element model.

[0287] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0288] Based on the state variables corresponding to the current incremental step, calculate the current moment's trial stress;

[0289] comparing the test stress with the stress at a previous moment;

[0290] In response to detecting that the probe stress is less than the stress, storing the probe stress;

[0291] In response to detecting that the tentative stress is greater than or equal to the stress, the state variable corresponding to the current incremental step is corrected based on the stress.

[0292] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0293] The deviator tensor of the stress is determined, and the equivalent stress model is determined based on the deviator tensor:

[0294] According to the correlation between the equivalent stress value and the equivalent trial stress value, a target scaling factor model is constructed;

[0295] Iteratively solving the target scaling factor model based on the Newton-Raphson algorithm, and updating the state variable before each iterative solution step;

[0296] In response to detecting that the value obtained by iterative solution is less than a preset threshold, solving the equivalent stress model according to the updated state variable to obtain the equivalent stress;

[0297] Based on the equivalent stress, a stress increment is determined, and the state variable corresponding to the current increment step is corrected according to the stress increment.

[0298] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0299] In response to detecting that the correlation between the equivalent stress value and the equivalent test stress value is When , constructing the target scaling factor model includes:

[0300]

[0301] Where δ represents the scaling factor, represents the equivalent stress value, represents the equivalent test stress value, G represents the shear modulus, D0 represents the limiting shear strain rate, Z represents the internal state variable of the material, n represents the material strain rate sensitivity constant, and Δt represents the incremental step size;

[0302] Before each iterative solution, the state variable is updated by an update function, and the update function includes:

[0303]

[0304]

[0305]

[0306]

[0307]

[0308]

[0309] Among them, Z I represents the isotropic hardening component, Δt represents the incremental step size, t represents the time, and Z D represents the kinematic hardening component, m1 represents Z I The evolution rate from the initial value Z0 to the final value Z1, represents plastic power, A1 and A2 represent thermal recovery coefficients, r1 and r2 represent thermal recovery indexes, Z2 represents an extreme value, m2 represents the evolution rate of β from the initial value Z0 to the final value Z3, u and V represent the directional tensor of stress σ and the second-order tensor β, respectively. Represents ZI The time derivative, represents the derivative with respect to β time, Z represents the internal state variable of the material, represents the time derivative of T, and T represents temperature.

[0310] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0311] Performing a convergence test on the correction result based on the subroutine;

[0312] In response to detecting that the test result meets a preset standard, determining that the correction result converges;

[0313] In response to detecting that the test result does not meet the preset standard, it is determined that the correction result does not converge.

[0314] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to memory, storage, database or other media used in the embodiments provided in this application can include non-volatile and / or volatile memory. Non-volatile memory can include read-only memory (ROM), programmable ROM (PROM), electrically programmable ROM (EPROM), electrically erasable programmable ROM (EEPROM) or flash memory. Volatile memory can include random access memory (RAM) or external cache memory. As an illustration and not limitation, RAM is available in many forms, such as static RAM (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double data rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), synchronous link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM (RDRAM), direct memory bus dynamic RAM (DRDRAM), and memory bus dynamic RAM (RDRAM).

[0315] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0316] The above-mentioned embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the invention patent. It should be pointed out that, for a person of ordinary skill in the art, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the patent of the present application shall be subject to the attached claims.

Claims

1. A finite element numerical simulation method, characterized in that: The method comprises: Based on the main program of the finite element analysis software, the geometric model corresponding to the target standard part is constructed; Meshing the geometric model according to the target unit type to obtain a first finite element model, setting boundary constraints and load conditions for the first finite element model to obtain a second finite element model; Solving the second finite element model, determining the state variables corresponding to the current increment step, and sending the state variables to a subroutine of the finite element analysis software; Calculating a trial stress based on the subroutine, determining a correction result of the state variable according to the trial stress, and returning the correction result to the main program when the correction result converges; In response to detecting the correction result, the main program is used to solve the state variables of the next incremental step until the state variables of a preset number of incremental steps are solved.

2. The finite element numerical simulation method according to claim 1, characterized in that: The method for generating the first finite element model comprises: Acquiring geometric parameters and material parameters of the target standard part; Using the main program of the finite element analysis software, an initial geometric model corresponding to the target standard part is constructed according to the geometric parameters and the preset proportions; Inputting the material parameters into the initial geometric model to obtain a target geometric model corresponding to the target standard part; The target geometric model is meshed according to the target unit type to generate the first finite element model.

3. The finite element numerical simulation method according to claim 1, characterized in that: The method for generating the second finite element model comprises: Based on the first preset mapping table, obtaining a target unit type for meshing the target standard part; Meshing the target geometric model according to the target unit type to obtain the first finite element model; Based on the second preset mapping table, the boundary constraint conditions and the load conditions corresponding to the target standard part are determined, and the boundary constraint conditions and the load conditions are set for the first finite element model to obtain the second finite element model.

4. The finite element numerical simulation method according to claim 1, characterized in that: Before determining the correction result of the state variable according to the trial stress, the method further includes: Based on the state variables corresponding to the current incremental step, calculate the current moment's trial stress; comparing the test stress with the stress at a previous moment; In response to detecting that the probe stress is less than the stress, storing the probe stress; In response to detecting that the tentative stress is greater than or equal to the stress, the state variable corresponding to the current incremental step is corrected based on the stress.

5. The finite element numerical simulation method according to claim 4, characterized in that: In response to detecting that the probe stress is greater than or equal to the stress, correcting the state variable based on the stress includes: The deviator tensor of the stress is determined, and the equivalent stress model is determined based on the deviator tensor: According to the correlation between the equivalent stress value and the equivalent trial stress value, a target scaling factor model is constructed; Iteratively solving the target scaling factor model based on the Newton-Raphson algorithm, and updating the state variable before each iterative solution step; In response to detecting that the value obtained by iterative solution is less than a preset threshold, solving the equivalent stress model according to the updated state variable to obtain the equivalent stress; Based on the equivalent stress, a stress increment is determined, and the state variable corresponding to the current increment step is corrected according to the stress increment.

6. The finite element numerical simulation method according to claim 5, characterized in that: According to the correlation between the equivalent stress value and the equivalent trial stress value, the target scaling factor model is constructed including: In response to detecting that the correlation between the equivalent stress value and the equivalent test stress value is When , constructing the target scaling factor model includes: Where δ represents the scaling factor, represents the equivalent stress value, represents the equivalent test stress value, G represents the shear modulus, D0 represents the limiting shear strain rate, Z represents the internal state variable of the material, n represents the material strain rate sensitivity constant, and Δt represents the incremental step size; Before each iterative solution, the state variable is updated by an update function, and the update function includes: Among them, Z I represents the isotropic hardening component, Δt represents the incremental step size, t represents the time, and Z D represents the kinematic hardening component, m1 represents Z I The evolution rate from the initial value Z0 to the final value Z1, represents plastic power, A1 and A2 represent thermal recovery coefficients, r1 and r2 represent thermal recovery indexes, Z2 represents an extreme value, m2 represents the evolution rate of β from the initial value Z0 to the final value Z3, u and V represent the directional tensor of stress σ and the second-order tensor β, respectively. Represents Z I The time derivative, represents the derivative with respect to β time, Z represents the internal state variable of the material, represents the time derivative of T, and T represents temperature.

7. The finite element numerical simulation method according to claim 1, characterized in that: Using the subroutine to detect whether the correction result converges includes: Performing a convergence test on the correction result based on the subroutine; In response to detecting that the test result meets a preset standard, determining that the correction result converges; In response to detecting that the test result does not meet the preset standard, it is determined that the correction result does not converge.

8. A finite element numerical simulation device, characterized in that: The device comprises: A geometric model building module is used to build a geometric model corresponding to the target standard part based on the main program of the finite element analysis software; A finite element model building module, used for meshing the geometric model according to a target unit type to obtain a first finite element model, setting boundary constraints and load conditions for the first finite element model to obtain a second finite element model; A solution module, used for solving the second finite element model, determining the state variables corresponding to the current increment step, and sending the state variables to a subroutine of the finite element analysis software; A correction module, used for calculating the trial stress based on the subroutine, determining the correction result of the state variable according to the trial stress, and returning the correction result to the main program when the correction result converges; The loop solving module is used to solve the state variables of the next incremental step by using the main program when the correction result is detected, until the state variables of a preset number of incremental steps are solved.

9. A computer 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, the method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images