Improved Johnson-Cook constitutive model embedding development implementation method for hardware wear simulation
By improving the Johnson-Cook constitutive model in the implicit algorithm material library of the finite element software and writing the UMAT subroutine, the problems of large errors and low computational efficiency in the wear simulation of transmission line hardware were solved, and a more accurate and stable wear simulation effect was achieved.
Patent Information
- Application Number
- CN202510901747.2
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-07-01
- Publication Date
- 2025-10-17
AI Technical Summary
The existing Johnson-Cook constitutive model has large errors and low computational efficiency when describing the wear behavior of transmission line hardware under impact loads such as instantaneous strong winds. The explicit algorithm may lead to unstable calculation results and poor convergence when dealing with complex problems.
The Johnson-Cook constitutive model is improved in the implicit algorithm material library of the finite element software, and the UMAT subroutine is written. The improved Johnson-Cook constitutive model and the finite element software main program are combined with the implicit algorithm to perform hardware wear simulation, thereby improving the accuracy and convergence of wear behavior.
It achieves a more accurate description of the wear characteristics of transmission line hardware under impact loads, improves the accuracy and computational efficiency of wear simulation, and solves the error and convergence problems of traditional models.
Smart Images
Figure CN120805567A_ABST
Abstract
Description
Technical Field
[0001] The present invention belongs to the technical field of power transmission and transformation equipment operation and maintenance, and in particular relates to a method for embedding and developing an improved Johnson-Cook constitutive model for hardware wear simulation. Background Art
[0002] When conducting finite element simulation of the wear process of transmission line hardware, the explicit algorithm material library of the finite element software provides users with the traditional Johnson-Cook constitutive model. There are two disadvantages to using the explicit algorithm and the traditional Johnson-Cook constitutive model: (1) The parameter identification process of the traditional Johnson-Cook constitutive model is complex and requires optimization through a combination of experimental data and numerical simulation. The process is tedious and time-consuming, and it cannot accurately describe the wear behavior of transmission line hardware under impact loads such as instantaneous strong winds, and there are large errors; (2) The explicit algorithm may encounter convergence problems when dealing with complex problems, especially when the time step is not properly selected, which can easily lead to unstable calculation results and a significant increase in the amount of calculation, resulting in low calculation efficiency. Therefore, this field urgently needs to improve the traditional Johnson-Cook constitutive model based on the above shortcomings. Summary of the Invention
[0003] In view of this, the present invention aims to overcome the shortcomings of the above-mentioned problems in the prior art, and proposes a method for embedding and developing an improved Johnson-Cook constitutive model for hardware wear simulation, and implements the improved Johnson-Cook constitutive model subroutine in the implicit algorithm material library of the finite element software. This method can more accurately predict the wear behavior of transmission line hardware under impact loads such as instantaneous strong winds, and improve the accuracy of transmission line hardware wear simulation.
[0004] To achieve the above object, the technical solution of the present invention is achieved as follows:
[0005] The first aspect of the present invention provides a method for developing and implementing an improved Johnson-Cook constitutive model for hardware wear simulation, comprising:
[0006] The improved Johnson-Cook constitutive model is
[0007]
[0008] Among them, σ eq is the equivalent stress, ε eq is the equivalent strain, ε eq * =ε eq / ε0 is the dimensionless equivalent plastic strain rate, ε0 is the reference strain rate, T * =(TT r) / (T m -T r ) is a dimensionless temperature, T is the current test temperature, T r is the reference temperature of the fitting material, T m is the melting point temperature of the fitting material, A, B1, B2, C, k, m are material parameters, wherein A, B1 and B2 are plastic deformation term parameters, and C, k and m are strain rate term parameters;
[0009] The improved Johnson-Cook constitutive model is realized by using the implicit algorithm material library of the finite element software, and the main program and the subroutine of the finite element software are reprogrammed, and the analysis step is set in the finite element software.
[0010] Further, at the beginning of the analysis step, the main program calls the UMAT subroutine, and the main program passes the initial values of the necessary variables to the corresponding variables of the UMAT subroutine, and the updated strain and Jacobian matrix are passed to the main program when the UMAT subroutine ends.
[0011] Further, the UMAT subroutine step is:
[0012] S1: reading the state data at time t;
[0013] S2: calculating the elastic matrix;
[0014] S3: calculating the trial stress at time t+Δt;
[0015] S4: calculating the stress value of the kth iteration;
[0016] S5: judging whether the yield condition is met;
[0017] S6: updating the plastic parameter;
[0018] S7: updating the strain and stress;
[0019] S8: updating the Jacobian matrix and updating the state variable PROPS.
[0020] Further, in the step S1, the elastic modulus, the Poisson's ratio, the improved Johnson-Cook constitutive model parameters and the Johnson-Cook damage model parameters at time t are read.
[0021] Further, in the step 3, the hardness and the trial stress are calculated according to the plastic deformation term of the improved Johnson-Cook constitutive model, and the trial stress calculation formula is as follows:
[0022]
[0023] Wherein ε p(0) is the strain plastic part at the 0th iteration, εp n is the plasticity part of strain at the nth analysis step, Δλ (0) is the plasticity parameter change at the 0th iteration, δ (0) is the stress value at the 0th iteration, C is the stiffness tensor of the fitting material, ε n+1 is the strain at the nth+1 analysis step.
[0024] Further, in the step S4, the stress value is implicitly updated according to the J2 flow criterion, and the stress value formula of the kth iteration is calculated as follows:
[0025] σ (k) = σ (0) - Δλ (k) C: r (k)
[0026] wherein δ (k) is the stress value at the kth iteration, Δλ (k) is the plasticity parameter change at the kth iteration, r (k) is the plasticity flow direction at the kth iteration.
[0027] Further, in the step S5, the Mises equivalent stress value is taken as the standard value to judge whether to yield or not, and the formula is as follows:
[0028] f (k) = σ (k) - σ Y (ε (k) )
[0029] wherein δ Y (ε (k) ) is the Mises equivalent stress value, if the convergence condition is met, the iteration is completed, and the strain is updated; if the convergence condition is not met, the step S4 is performed.
[0030] Further, in the step S6, the plasticity parameter is updated by using the Newton-Raphson method, and the formula is as follows:
[0031] Δλ (k+1) = Δλ (k) + δλ (k) .
[0032] Further, in the step S7, the strain and stress are updated, and the formula is as follows:
[0033] ε p(k+1) = ε p(k) - δλ (k) r (k)
[0034] σ (k+1) = C: (ε n+1 - ε p(k+1) )
[0035] where ε p(k+1) is the plastic part of strain at the k+1 iteration, ε p(k) is the plastic part of strain at the k iteration, δλ (k) is the plastic parameter increment at the k iteration, the strain, stress and plastic parameter after the k+1 iteration are substituted into step S5 to judge the yield condition until the convergence condition is met.
[0036] Further, in step S8, the Jacobian matrix and the state variable PROPS are updated as the initial values for the next time step.
[0037] The second aspect of the present application provides an electronic device, comprising a processor and a memory connected with the processor and used for storing executable instructions of the processor, and the processor is used for executing the above-mentioned method for developing and implementing an improved Johnson-Cook constitutive model embedded in a hardware wear simulation.
[0038] The third aspect of the present application provides a computer readable storage medium, which stores a computer program, and the computer program is executed by a processor to realize the above-mentioned method for developing and implementing an improved Johnson-Cook constitutive model embedded in a hardware wear simulation.
[0039] Compared with the prior art, the method for developing and implementing an improved Johnson-Cook constitutive model embedded in a hardware wear simulation has the following advantages:
[0040] (1) When simulating the wear of hardware of a power transmission line, the existing constitutive model in the finite element software is the traditional Johnson-Cook constitutive model. The traditional Johnson-Cook constitutive model has a large error in describing the wear behavior of the hardware of the power transmission line under the action of instantaneous strong wind and other impact loads. The present application proposes an improved Johnson-Cook constitutive model which can be embedded in the finite element software, and the model can more accurately describe the wear characteristics of the hardware of the power transmission line under the action of impact loads.
[0041] (2) The traditional Johnson-Cook constitutive model only exists in the material library of the explicit algorithm of the finite element software. Compared with the implicit algorithm, the explicit algorithm has lower accuracy in describing the wear behavior of the hardware of the power transmission line. The present application realizes the UMAT subprogram of the improved Johnson-Cook constitutive model in the material library of the implicit algorithm, and further improves the accuracy and convergence of the wear simulation. BRIEF DESCRIPTION OF DRAWINGS
[0042] The accompanying drawings, which form a part of this application, are included to provide a further understanding of the application and are incorporated in and constitute a part of this application. The embodiments of the application illustrate the application and together with the description serve to explain the principles of the application. In the drawings:
[0043] Figure 1 The flow chart of the improved Johnson-Cook constitutive model subprogram realized by the finite element software implicit algorithm material library in the embodiment of the application. DETAILED DESCRIPTION
[0044] It should be noted that the embodiments in the application and the features in the embodiments can be combined with each other without conflict.
[0045] The embodiments of the present disclosure are described in detail below, and examples of the embodiments are shown in the drawings, wherein the same or similar notations represent the same or similar elements or elements having the same or similar functions throughout. The embodiments described below by referring to the drawings are exemplary and are intended to explain the present disclosure, and cannot be understood as a limitation of the present disclosure.
[0046] Embodiment 1:
[0047] The application provides a method for embedding and developing an improved Johnson-Cook constitutive model for fitting wear of fittings, and the improved Johnson-Cook constitutive model is
[0048]
[0049] wherein σ eq is an equivalent stress, ε eq is an equivalent strain, ε eq * = ε eq / ε0 is a dimensionless equivalent plastic strain rate, ε0 is a reference strain rate, T * = (T-T r ) / (T m -T r ) is a dimensionless temperature, T is a current test temperature, T r is a reference temperature of the fitting material, T m is a melting point temperature of the fitting material, A, B1, B2, C, k, and m are material parameters, wherein A, B1, and B2 are all plastic deformation term parameters, and C, k, and m are all strain rate term parameters.
[0050] The improved Johnson-Cook constitutive model is realized by the finite element software implicit algorithm material library, and then the finite element software main program and subprogram are re-written, and the finite element software analysis step is set.
[0051] In the embodiment of the present application, the improved Johnson-Cook constitutive model subroutine block diagram realized by the finite element software implicit algorithm material library is shown in Fig. 1. When the analysis step starts, the main program calls the UMAT subroutine, the main program passes the initial values of necessary variables to the corresponding variables of the UMAT subroutine, and the updated strain and Jacobian matrix of the UMAT subroutine at the end are passed to the main program. Figure 1
[0052] The corresponding UMAT subroutine variable declaration part is as follows:
[0053] SUBROUTINE UMAT (STRESS, STATEV, DDSDDE, SSE, SPD, SCD, 1RPL, DDSDDT, DRPLDE, DRPLDT, STRAN, DSTRAN, TIME, 2DTIME, TEMP, DTEMP, PREDEF, DPRED, MATERL, NDI, NSHR, 3NTENS, NSTATV, PROPS, NPROPS, COORDS, DROT, PNEWDT, 4CELENT, DFGRD0, DFGRD1, NOEL, NPT, KSLAY, KSPT, KSTEP, 5KINC) INCLUDE 'ABA_PARAM.INC'
[0054] …
[0055] In the embodiment of the present application, the basic steps of the UMAT subroutine are as follows:
[0056] (1) In step S1, the state data at time t is read, and the reading program is as follows:
[0057] E = PROPS (1)! Elastic modulus
[0058] V = PROPS (2)! Poisson's ratio
[0059] A = PROPS (3)! Improved Johnson-Cook constitutive model parameter A
[0060] B1 = PROPS (4)! Improved Johnson-Cook constitutive model parameter B1
[0061] B2 = PROPS (5)! Improved Johnson-Cook constitutive model parameter B2
[0062] C = PROPS (6)! Improved Johnson-Cook constitutive model parameter C
[0063] K = PROPS (7)! Improved Johnson-Cook constitutive model parameter k
[0064] M = PROPS(8)! Improved Johnson-Cook constitutive model parameter m
[0065] D1 = PROPS(9)! Johnson-Cook damage model parameter D1
[0066] D2 = PROPS(10)! Johnson-Cook damage model parameter D2
[0067] D3 = PROPS(11)! Johnson-Cook damage model parameter D3
[0068] D4 = PROPS(12)! Johnson-Cook damage model parameter D4
[0069] D5 = PROPS(13)! Johnson-Cook damage model parameter D5
[0070] EDOT = PROPS(14)! Reference strain rate
[0071] TR = PROPS(15)! Reference temperature
[0072] TM = PROPS(16)! Melting temperature
[0073] After reading the data, the Lame constant λ and G are calculated.
[0074] (2) In step S2, the elastic matrix C is constructed according to the isotropic criterion;
[0075] (3) In step S3, the trial stress at time t+Δt is calculated.
[0076]
[0077] where ε p(0) is the strain plastic part at the 0th iteration, ε p n is the strain plastic part at the nth analysis step, Δλ (0) is the plastic parameter change amount at the 0th iteration, δ (0) is the stress value at the 0th iteration, C is the stiffness tensor of the hardware material, ε n+1 is the strain at the n+1th analysis step.
[0078] The elastic part calculation subroutine is:
[0079] SUBROUTINE USERHARD(SYIELD, HARD, EQPLAS, TABLE)
[0080] INCLUDE 'ABA_PARAM.INC'
[0081] DIMENSION TABLE(3)
[0082] A = TABLE(1)
[0083] B1 = TABLE(2)
[0084] B2 = TABLE(3)
[0085] HARD = 0.0
[0086] IF (EQPLAS.EQ. 0.0) THEN
[0087] SYIELD = A
[0088] ELSE HARD = B1 + 2*B2*EQPLAS SYIELD = A + B1*EQPLAS + B2*EQPLAS**2
[0089] END IF
[0090] RETURN
[0091] END
[0092] (4) In step S4, the stress value of the kth iteration is calculated.
[0093] σ (k) = σ (0) - Δλ (k) C:r (k) (3)
[0094] wherein δ (k) is the stress value of the kth iteration, Δλ (k) is the plastic parameter variation amount of the kth iteration, and r (k) is the plastic flow direction of the kth iteration.
[0095] (5) In step S5, it is judged whether the yield condition is satisfied.
[0096] f (k) = σ (k) - σ Y (ε (k) ) (4)
[0097] wherein δ Y (ε (k) ) is the Mises equivalent stress value. If the convergence condition is satisfied, the iteration is completed and the strain is updated; if the convergence condition is not satisfied, the fourth step is performed.
[0098] (6) In step S6, the plastic parameter
[0099] Δλ(k+1) = Δλ (k) + δλ (k) (5)
[0100] where the numerical iteration part subroutine is SYIELD = SYIELD0
[0101] DEQPL = (SMISES - SYIELD) / EG3
[0102] DEQMIN = HALF * DTIME * EXP(1.0D-3 / D)
[0103] DO 87 KEWTON = 1, NEWTON
[0104] DEQPL = MAX(DEQPL, DEQMIN)
[0105] CALL USERHARD(SYIELD, HARD, EQPLAS + DEQPL, PROPS(4))
[0106] TVP = ONE + (DEQPL / DTIME) ** D
[0107] HARD1 = HARD * TVP + SYIELD * D * (DEQPL / DTIME) ** (D-1) / DEQPL SYIELD = SYIELD * TVP
[0108] RHS = SMISES - EG3 * DEQPL - SYIELD * DEQPL DEQPL = DEQPL + RHS / (EG3 + HARD1)
[0109] (7) In step S7, the strain and stress are updated.
[0110] ε p(k+1) = ε p(k) - δλ (k) r (k) (6)
[0111] σ (k+1) = C: (ε n+1 - ε p(k+1) ) (7)
[0112] where ε p(k+1) is the strain plasticity part at the k+1th iteration, ε p(k) is the strain plasticity part at the kth iteration, and δλ (k) is the plasticity parameter increment at the kth iteration. The strain, stress and plasticity parameter after the k+1th iteration are substituted into equation (4) to judge the yield condition until the convergence condition is met.
[0113] (8) In step S8, the Jacobian matrix is updated, and the state variable PROPS is updated.
[0114] Embodiment 2
[0115] An electronic device includes a processor and a memory connected to the processor in communication, and used for storing instructions executable by the processor, the processor being configured to execute the above-mentioned method for developing and implementing a Johnson-Cook constitutive model embedded in a simulation of wear of a fitting.
[0116] Embodiment 3
[0117] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the above-mentioned method for developing and implementing a Johnson-Cook constitutive model embedded in a simulation of wear of a fitting.
[0118] The collection, storage, use, processing, transmission, provision and disclosure of user personal information involved in the present disclosure comply with relevant laws and regulations and do not violate public order and good customs.
[0119] It should be noted that personal information from users should be collected for legitimate and reasonable purposes and not shared or sold outside these legitimate uses. In addition, such collection / sharing should be carried out after the user's informed consent is received, including but not limited to informing the user to read the user agreement / user notice before the user uses the function, and signing an agreement / authorization including authorization of relevant user information. In addition, any necessary steps should be taken to protect and secure access to such personal information data and ensure that other people with access to personal information data comply with their privacy policies and processes.
[0120] The present disclosure contemplates that users may opt in to the collection of personal information data across multiple devices, applications, and websites owned or operated by the present disclosure. The present disclosure contemplates that users may opt in to the sharing of personal information data, including but not limited to, across multiple devices, applications, and websites owned or operated by the present disclosure. For example, a user may opt in to the sharing of personal information data from a mobile device to a website owned or operated by the present disclosure. The present disclosure contemplates that users may opt in to the sharing of personal information data from a website to a mobile device. The present disclosure contemplates that users may opt in to the sharing of personal information data from a website to an application. The present disclosure contemplates that users may opt in to the sharing of personal information data from an application to a website. The present disclosure contemplates that users may opt in to the sharing of personal information data from one mobile device to another mobile device. The present disclosure contemplates that users may opt in to the sharing of personal information data from one mobile application to another mobile application. The present disclosure contemplates that users may opt in to the sharing of personal information data from a mobile device to a website or application associated with a different user.
[0121] The acquisition, transmission, storage, use, processing, etc. of data in the technical solutions of the present disclosure comply with relevant provisions of national laws and regulations.
[0122] It should be noted that in the embodiments of the present disclosure, some industry existing solutions, components, models, etc. may be mentioned, which should be considered as exemplary, and the purpose is only to illustrate the feasibility in the implementation of the technical solutions of the present disclosure, but it does not mean that the applicant has or will necessarily use the solutions.
[0123] In the descriptions of the aforementioned embodiments, the reference terms "one embodiment", "some embodiments", "example", "specific example", or "some examples" mean that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present disclosure. In this specification, the schematic expressions of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described may be combined in any one or more embodiments or examples in a suitable manner. In addition, those skilled in the art may combine and combine the different embodiments or examples described in this specification and the features of the different embodiments or examples, unless they are mutually inconsistent.
[0124] Furthermore, the terms "first" and "second" are used for descriptive purposes only and should not be construed as indicating or implying relative importance or implicitly specifying the number of technical features being referred to. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one such feature. Throughout the present disclosure, "plurality" means at least two, such as two, three, etc., unless otherwise specifically defined.
[0125] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, segment or portion of code comprising one or more executable instructions for implementing the steps of a custom logical function or process, and the scope of the preferred embodiments of the present disclosure includes additional implementations in which functions may be performed out of the order shown or discussed, including performing functions in a substantially simultaneous manner or in the reverse order depending on the functions involved, which should be understood by those skilled in the art to which the embodiments of the present disclosure belong.
[0126] The logic and / or steps represented in the flowcharts and / or described herein, for example, can be considered as a sequence of instructions to implement logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, apparatus, or device, such as a computer-based system, processor- based system, or other system that can fetch the instructions from the instruction execution system, apparatus, or device and execute the instructions. For purposes of this specification, a "computer-readable medium" can be any apparatus that can contain, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. The computer-readable medium can be a computer- readable storage medium or a computer-readable signal medium. The computer-readable storage medium can be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any suitable combination of the foregoing. More specific examples (a non-exhaustive list) of the computer-readable storage medium include the following: an electrical connection having one or more wires (electrical connections), a portable computer diskette (a magnetic device), a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), an optical fiber, and a portable compact disc read-only memory (CDROM). In addition, the computer-readable medium can even be paper or another suitable medium upon which the program is printed, as the program can be electronically captured, for example, via optical scanning of the paper or other medium, then compiled, interpreted, or otherwise processed in a suitable manner, if necessary, and then stored in a computer memory.
[0127] It should be understood that aspects of the present disclosure can be implemented in hardware, software, firmware or combinations thereof. In the above embodiments, the various steps or methods can be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. As such, if implemented in hardware, and in another embodiment, any of the following technologies, known in the art, or their combinations can be used: discrete logic circuitry having logic gates for implementing logic functions on data signals, application specific integrated circuits having appropriate combinational logic gates, programmable gate arrays (PGA), field programmable gate arrays (FPGA), and the like.
[0128] Those skilled in the art can understand that all or part of the steps carried out by the above-mentioned embodiment method can be completed by programs instructing related hardware, and the programs can be stored in a computer-readable storage medium. When the programs are executed, one or a combination of the steps of the method embodiment is included.
[0129] In addition, each of the functional units in the various embodiments of the present disclosure can be integrated in one processing module, or each unit can exist physically separately, or two or more units can be integrated in one module. The integrated module can be implemented in the form of hardware or in the form of a software functional module. When the integrated module is implemented in the form of a software functional module and sold or used as an independent product, it can also be stored in a computer readable storage medium.
[0130] The storage medium mentioned above can be a read-only memory, a magnetic disk or an optical disk, etc. Although the embodiments of the present disclosure have been shown and described above, it should be understood that the above embodiments are exemplary and should not be construed as limiting the present disclosure, and those skilled in the art can make changes, modifications, replacements and variations to the above embodiments within the scope of the present disclosure.
Claims
1. A method for developing and implementing an improved Johnson-Cook constitutive model embedded in hardware wear simulation, characterized by: include The improved Johnson-Cook constitutive model is Among them, σ eq is the equivalent stress, ε eq is the equivalent strain, ε eq * =ε eq / ε0 is the dimensionless equivalent plastic strain rate, ε0 is the reference strain rate, T * =(TT r ) / (T m -T r ) is the dimensionless temperature, T is the current test temperature, T r is the reference temperature of the hardware material, T m is the melting point temperature of the hardware material, A, B1, B2, C, k, and m are material parameters, among which A, B1, and B2 are plastic deformation parameters, and C, k, and m are strain rate parameters; The improved Johnson-Cook constitutive model is implemented using the implicit algorithm material library of the finite element software, and then the main program and subroutines of the finite element software are written, and the analysis steps are set in the finite element software.
2. The method for embedding and developing an improved Johnson-Cook constitutive model for hardware wear simulation according to claim 1, characterized in that: At the beginning of the analysis step, the main program calls the UMAT subroutine, and the main program passes the necessary initial values of variables to the corresponding variables of the UMAT subroutine. At the end of the UMAT subroutine, the updated strain and Jacobian matrix are passed to the main program.
3. The method for embedding and developing an improved Johnson-Cook constitutive model for hardware wear simulation according to claim 2, characterized in that: The UMAT subroutine steps are: S1: Read the status data at time t; S2: Calculate and construct the elasticity matrix; S3: Calculate the trial stress at time t+Δt; S4: Calculate the stress value of the kth iteration; S5: Determine whether the yield condition is met; S6: Update plasticity parameters; S7: Update strain and stress; S8: Update the Jacobian matrix and the state variable PROPS.
4. The method for embedding and developing an improved Johnson-Cook constitutive model for hardware wear simulation according to claim 3, characterized in that: In the step S1, the elastic modulus, Poisson's ratio, improved Johnson-Cook constitutive model parameters, and Johnson-Cook damage model parameters at time t are read.
5. The method for embedding and developing an improved Johnson-Cook constitutive model for hardware wear simulation according to claim 3, characterized in that: In step 3, the hardness and trial stress are calculated according to the plastic deformation term of the improved Johnson-Cook constitutive model. The trial stress calculation is shown as follows: Among them, ε p(0) is the strain plastic part at the 0th iteration, ε p n is the strain plastic part at the nth analysis step, Δλ (0) is the change of plasticity parameters at the 0th iteration, δ (0) is the stress value at the 0th iteration, C is the stiffness tensor of the hardware material, ε n+1 is the strain at the n+1th analysis step.
6. The method for embedding and developing an improved Johnson-Cook constitutive model for hardware wear simulation according to claim 5, characterized in that: In step S4, the stress value is implicitly updated according to the J2 flow criterion, and the stress value of the kth iteration is calculated as follows: s (k) =s (0) -Dl (k) C:r (k) where δ (k) is the stress value at the kth iteration, Δλ (k) is the change of plasticity parameter at the kth iteration, r (k) is the plastic flow direction at the kth iteration.
7. The method for embedding and developing an improved Johnson-Cook constitutive model for hardware wear simulation according to claim 3, characterized in that: In step S5, the Mises equivalent stress value is used as the standard value to determine whether yield occurs, and the formula is as follows: f (k) =s (k) -s Y (e (k) ) Among them, δ Y (ε (k) ) is the Mises equivalent stress value. If the convergence condition is met, the iteration is completed and the strain is updated; if the convergence condition is not met, step S4 is performed.
8. The method for embedding and developing an improved Johnson-Cook constitutive model for hardware wear simulation according to claim 6, characterized in that: In step S6, the Newton-Raphson method is used to update the plasticity parameters, and the formula is as follows: Dl (k+1) =Dl (k) +sl (k) 。 9. The method for embedding and developing an improved Johnson-Cook constitutive model for hardware wear simulation according to claim 6, characterized in that: In step S7, the strain and stress formulas are updated as follows: e p(k+1) =e p(k) -dl (k) r (k) Among them, ε p(k+1) is the strain plastic part at the k+1th iteration, ε p(k) is the strain plasticity part at the kth iteration, δλ (k) is the plastic parameter increment at the kth iteration. The strain, stress and plastic parameters after the k+1th iteration are substituted into step S5 to judge the yield condition until the convergence condition is met.
10. The method for embedding and developing an improved Johnson-Cook constitutive model for hardware wear simulation according to claim 3, characterized in that: In step S8, the Jacobian matrix and state variables PROPS are updated as initial values for the next time step.