Gear carburizing and quenching deformation calculation method, system, equipment and storage medium

Abstract

The present application discloses a method, system, device and storage medium for calculating gear carburizing and quenching deformation. The method for calculating gear carburizing and quenching deformation includes obtaining a first node coordinate set of the gear, material thermal conductivity and material specific heat capacity, calculating the carbon content of a target time period based on the first node coordinate set by using Fick's second law, and obtaining a TTT curve by testing according to the carbon content; calculating the temperature of the target time period according to the material thermal conductivity and the material specific heat capacity; calculating the first volume fraction and the second volume fraction of the target time period according to the carbon content, the TTT curve and the temperature; calculating the plastic strain increment of the target time period according to the first volume fraction and the second volume fraction; calculating the second node coordinate set of the gear through an Abaqus solver based on the plastic strain increment; and calculating the deformation value of the gear according to the first node coordinate set and the second node coordinate set, thereby improving the accuracy of the calculation of gear carburizing and quenching deformation.

Description

Technical Field

[0001] The present application relates to a technical field related to gear carburizing and quenching deformation calculation, and in particular to a gear carburizing and quenching deformation calculation method, system, device and storage medium. Background Art

[0002] Carburizing and quenching process is an important part in the heat treatment process of low carbon steel. After carburizing, the surface carbon concentration increases, the hardenability is enhanced, and a high carbon solid is formed on the surface, which can improve the surface hardness and wear resistance of parts and extend the service life of parts. However, the phase change and temperature non-uniformity during carburizing and quenching cause parts to deform, further hindering the improvement of part precision. When the deformation of parts is too large, it will not only increase the trial production cycle of parts, but also reduce the fatigue life of parts during service.

[0003] At present, for the carburizing process, only the effect of carbon content on the material's thermophysical parameters, phase transition kinetics and mechanical properties is considered. In the actual process, carbon atoms dissolve into the gaps in the metal crystal lattice to form a solid solution. Due to the different atomic sizes of the solvent and the solute, the lattice is locally distorted. For interstitial solid solutions, as solute atoms are added, the solid solution lattice constant always increases with the addition of solute atoms. In order to accurately calculate the additional deformation caused by carburizing and quenching, the interstitial solid solution deformation caused by carburizing cannot be ignored. Therefore, the existing numerical models cannot accurately predict carburizing and quenching deformation. Summary of the invention

[0004] The present application aims to at least solve the technical problems existing in the prior art. To this end, the present application proposes a gear carburizing and quenching deformation calculation method, system, device and storage medium, which can consider the phase change thermoelasticity of carburizing expansion strain and improve the accuracy of gear carburizing and quenching deformation calculation.

[0005] In a first aspect of the present application, a method for calculating gear carburizing and quenching deformation is provided, comprising the following steps:

[0006] Obtaining a first node coordinate set, material thermal conductivity, and material specific heat capacity of the gear, wherein the first node coordinate set is an initial coordinate set of all nodes after meshing the gear before deformation;

[0007] Calculating the carbon content of the target period by using Fick's second law based on the first node coordinate set, and performing a test according to the carbon content to obtain a TTT curve;

[0008] Calculating the temperature of the target time period according to the thermal conductivity of the material and the specific heat capacity of the material;

[0009] Calculating a first volume fraction and a second volume fraction for a target period according to the carbon content, the TTT curve, and the temperature;

[0010] Calculating a plastic strain increment in a target period according to the first volume fraction and the second volume fraction;

[0011] Based on the plastic strain increment, calculating a second node coordinate set of the gear by an Abaqus solver;

[0012] The deformation value of the gear is calculated according to the first node coordinate set and the second node coordinate set.

[0013] The control method according to the embodiment of the present application has at least the following beneficial effects:

[0014] The method obtains the first node coordinate set of the gear, the thermal conductivity of the material and the specific heat capacity of the material, wherein the first node coordinate set is the initial coordinate set of all nodes after the gear before deformation is meshed; the carbon content of the target time period is calculated based on the first node coordinate set by Fick's second law, and the TTT curve is obtained by testing according to the carbon content; the temperature of the target time period is calculated according to the thermal conductivity of the material and the specific heat capacity of the material; the first volume fraction and the second volume fraction of the target time period are calculated according to the carbon content, the TTT curve and the temperature; the plastic strain increment of the target time period is calculated according to the first volume fraction and the second volume fraction; based on the plastic strain increment, the second node coordinate set of the gear is calculated by the Abaqus solver; the deformation value of the gear is calculated according to the first node coordinate set and the second node coordinate set. The present application improves the accuracy of the calculation of the carburizing and quenching deformation of the gear by considering the phase change thermoelastic-plasticity of the carburizing expansion strain.

[0015] According to some embodiments of the present application, the calculating the carbon content of the target period by using Fick's second law based on the first node coordinate set includes:

[0016] Calculating the carbon diffusion coefficient according to the thermal conductivity of the material and the specific heat capacity of the material;

[0017] The carbon content of the target period is calculated according to the carbon diffusion coefficient and the first node coordinate set by Fick's second law using the following formula:

[0018]

[0019] Where C is the carbon content, t s is the carburizing time corresponding to the target period, x is the X coordinate of the first node coordinate set, y is the Y coordinate of the first node coordinate set, and z is the Z coordinate of the first node coordinate set.

[0020] According to some embodiments of the present application, calculating the temperature of the target time period according to the thermal conductivity of the material and the specific heat capacity of the material includes:

[0021] Obtaining the released latent heat value during the target period;

[0022] The temperature of the target period is calculated according to the released latent heat value, the thermal conductivity of the material and the specific heat capacity of the material by the following formula:

[0023]

[0024] Where T is the temperature, t q is the quenching time corresponding to the target period, ρ is the gear density, λ is the thermal conductivity of the gear material, c is the specific heat capacity of the gear material, is the Laplace operator, and W is the released latent heat.

[0025] According to some embodiments of the present application, calculating the first volume fraction and the second volume fraction of the target time period according to the carbon content, the TTT curve and the temperature includes:

[0026] Calculate the first moment, the second moment, the first transition amount and the second transition amount corresponding to the temperature T;

[0027] The first volume fraction of the target time period is calculated according to the temperature by the following formula:

[0028]

[0029] Where b is the first parameter, n is the second parameter, V is the first volume fraction of the diffusion-type phase transition, and t q is the quenching time corresponding to the target period, t1 is the first moment corresponding to the temperature T in the TTT curve, t2 is the second moment corresponding to the temperature T in the TTT curve, V1 is the first transition value corresponding to t1, and V2 is the second transition value corresponding to t2;

[0030] calculating martensite mobility based on the temperature and the carbon content;

[0031] Calculate phase transformation kinetic parameters and martensite start transformation temperature value according to the carbon content and the TTT curve;

[0032] The second volume fraction of the target time period is calculated according to the martensite mobility, the phase transformation kinetic parameter and the martensite start transformation temperature value by the following formula:

[0033]

[0034] Among them, V M is the second volume fraction of the non-diffusion phase transition, v m (T,C) is the martensite mobility, α is the phase transition kinetic parameter 1, β is the phase transition kinetic parameter 2, M s is the temperature at which martensite transformation begins.

[0035] According to some embodiments of the present application, calculating the plastic strain increment of the target time period according to the first volume fraction and the second volume fraction includes:

[0036] Get the total strain increment, stress increment and stress generated by UMATHT;

[0037] Calculating a third volume fraction of austenite transformation according to the first volume fraction and the second volume fraction;

[0038] The elastic strain increment is calculated according to the first volume fraction, the second volume fraction, the third volume fraction, the stress increment and the stress by the following formula:

[0039]

[0040] in, is the elastic strain increment, is the inverse matrix of the elasticity matrix, dσ ij is the stress increment, σ is the stress, when k is 1, V k is the third volume fraction. When k is 2, V k is the first volume fraction. When k is 3, V k is the second volume fraction;

[0041] The temperature strain increment is calculated according to the first volume fraction, the second volume fraction and the third volume fraction by the following formula:

[0042]

[0043] in, is the temperature strain increment, T0 is the temperature corresponding to the initial time of the target period, T1 is the temperature corresponding to the end time of the target period, when k is 1, α k is the thermal expansion coefficient of austenite phase transformation. When k is 2, α k is the thermal expansion coefficient of the diffusion phase transition. When k is 3, α k is the thermal expansion coefficient of the non-diffusion phase transition, δ ij is the Kronecker function;

[0044] The phase change strain increment is calculated according to the first volume fraction and the second volume fraction by the following formula:

[0045]

[0046] in, is the phase change strain increment. When k is 2, β k is the structural expansion of the diffusion phase transition. When k is 3, β kStructural expansion for non-diffusion phase transition;

[0047] The phase transformation plastic strain increment is calculated according to the first volume fraction and the second volume fraction by the following formula:

[0048]

[0049] in, is the phase transformation plastic strain increment, R is the phase transformation plasticity coefficient, S ij is the deviatoric stress tensor;

[0050] The carburizing strain increment is calculated according to the carbon content by the following formula:

[0051]

[0052] in, is the carburizing strain increment, α c is the carburizing parameter value corresponding to the carbon content;

[0053] The plastic strain increment of the target time period is calculated according to the carburizing strain increment, the phase transformation plastic strain increment, the phase transformation strain increment, the temperature strain increment, the elastic strain increment and the total strain increment by the following formula:

[0054]

[0055] in, is the plastic strain increment, dε ij is the total strain increment.

[0056] According to some embodiments of the present application, the calculating the second node coordinate set of the gear by using an Abaqus solver based on the plastic strain increment includes:

[0057] Based on the plastic strain increment, the carburizing and quenching deformation values ​​of all nodes of the gear are calculated by the Abaqus solver;

[0058] The second node coordinate set is calculated according to the carburizing and quenching deformation values ​​of all the nodes.

[0059] According to some embodiments of the present application, the deformation value of the gear is calculated by the following formula:

[0060]

[0061] Among them, d is the deformation value of the gear node, x0 is the X coordinate of the first node coordinate, y0 is the Y coordinate of the first node coordinate, z0 is the Z coordinate of the first node coordinate, x1 is the X coordinate of the second node coordinate, y1 is the Y coordinate of the second node coordinate, and z1 is the Z coordinate of the second node coordinate.

[0062] In a second aspect of the present application, a gear carburizing and quenching deformation calculation system is provided, and the gear carburizing and quenching deformation calculation system comprises:

[0063] A data acquisition module, used to acquire a first node coordinate set, material thermal conductivity and material specific heat capacity of the gear, wherein the first node coordinate set is an initial coordinate set of all nodes after meshing the gear before deformation;

[0064] A carbon content calculation module, used to calculate the carbon content of a target period by using Fick's second law based on the first node coordinate set, and to obtain a TTT curve by testing according to the carbon content;

[0065] A temperature calculation module, used for calculating the temperature of a target time period according to the thermal conductivity of the material and the specific heat capacity of the material;

[0066] a volume fraction calculation module, configured to calculate a first volume fraction and a second volume fraction in a target time period according to the carbon content, the TTT curve and the temperature;

[0067] A plastic strain increment calculation module, used for calculating the plastic strain increment of a target time period according to the first volume fraction and the second volume fraction;

[0068] A second node coordinate set calculation module, used for calculating a second node coordinate set of the gear through an Abaqus solver based on the plastic strain increment;

[0069] A deformation value calculation module is used to calculate the deformation value of the gear according to the first node coordinate set and the second node coordinate set.

[0070] The system obtains the first node coordinate set of the gear, the thermal conductivity of the material and the specific heat capacity of the material, wherein the first node coordinate set is the initial coordinate set of all nodes after the gear before deformation is meshed; the carbon content of the target time period is calculated based on the first node coordinate set by Fick's second law, and the TTT curve is obtained by testing according to the carbon content; the temperature of the target time period is calculated according to the thermal conductivity of the material and the specific heat capacity of the material; the first volume fraction and the second volume fraction of the target time period are calculated according to the carbon content, the TTT curve and the temperature; the plastic strain increment of the target time period is calculated according to the first volume fraction and the second volume fraction; based on the plastic strain increment, the second node coordinate set of the gear is calculated by the Abaqus solver; the deformation value of the gear is calculated according to the first node coordinate set and the second node coordinate set. The application improves the accuracy of the calculation of the carburizing and quenching deformation of the gear by considering the phase change thermal elastic-plasticity of the carburizing expansion strain.

[0071] According to a third aspect of the present application, there is provided an electronic device for calculating gear carburizing and quenching deformation, comprising at least one control processor and a memory for communicating with the at least one control processor; the memory stores instructions executable by the at least one control processor, and the instructions are executed by the at least one control processor so that the at least one control processor can execute the above-mentioned gear carburizing and quenching deformation calculation method.

[0072] In a fourth aspect of the present application, a computer-readable storage medium is provided, wherein the computer-readable storage medium stores computer-executable instructions, and the computer-executable instructions are used to enable a computer to execute the above-mentioned gear carburizing and quenching deformation calculation method.

[0073] It should be noted that the beneficial effects between the second to fourth aspects of the present application and the prior art are the same as the beneficial effects between the above-mentioned gear carburizing and quenching deformation calculation system and the prior art, and will not be described in detail here.

[0074] Additional aspects and advantages of the present application will be given in part in the description below, and in part will become apparent from the description below, or will be learned through the practice of the present application. BRIEF DESCRIPTION OF THE DRAWINGS

[0075] The above and / or additional aspects and advantages of the present application will become apparent and easily understood from the description of the embodiments in conjunction with the following drawings, in which:

[0076] Figure 1 It is a flow chart of an embodiment of a method for calculating gear carburizing and quenching deformation provided by the present application;

[0077] Figure 2 It is a schematic diagram of gear mesh division in an embodiment of a gear carburizing and quenching deformation calculation method provided by the present application;

[0078] Figure 3 It is a schematic diagram of gear carburizing field distribution in an embodiment of the gear carburizing quenching deformation calculation method provided by the present application;

[0079] Figure 4 It is a schematic diagram of the gear temperature field distribution of an embodiment of the gear carburizing and quenching deformation calculation method provided by the present application;

[0080] Figure 5 It is a schematic diagram of the gear organization field distribution of an embodiment of the gear carburizing and quenching deformation calculation method provided by the present application;

[0081] Figure 6 It is a schematic diagram of gear bainite and martensite in an embodiment of the gear carburizing and quenching deformation calculation method provided by the present application;

[0082] Figure 7is a schematic diagram of gear deformation distribution in an embodiment of a method for calculating gear carburizing and quenching deformation provided by the present application;

[0083] Figure 8 It is a schematic diagram of mesh division of an embodiment of a method for calculating gear carburizing and quenching deformation provided by the present application;

[0084] Fig. 9 It is a schematic diagram of the process route of an embodiment of the gear carburizing and quenching deformation calculation method provided by the present application;

[0085] Fig.10 is a schematic diagram of carburizing field distribution in an embodiment of a method for calculating gear carburizing and quenching deformation provided by the present application;

[0086] Fig.11 It is a schematic diagram of temperature field distribution of an embodiment of a method for calculating gear carburizing and quenching deformation provided by the present application;

[0087] Fig.12 It is a schematic diagram of the tissue field distribution of an embodiment of the gear carburizing and quenching deformation calculation method provided by the present application;

[0088] Fig.13 It is a structural schematic diagram of an embodiment of a gear carburizing and quenching deformation calculation system provided by the present application;

[0089] Fig.14 It is a schematic diagram of the structure of an embodiment of the electronic device provided by the present application. DETAILED DESCRIPTION

[0090] The embodiments of the present application are described in detail below, and examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are only used to explain the present application, and cannot be understood as limiting the present application.

[0091] In the description of this application, if there is a description of first, second, etc., it is only for the purpose of distinguishing technical features, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features or implicitly indicating the order of the indicated technical features.

[0092] In the description of the present application, it should be understood that the descriptions involving orientation, such as the orientation or positional relationship indicated as up, down, etc., are based on the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the present application and simplifying the description, and do not indicate or imply that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore should not be understood as a limitation on the present application.

[0093] In the description of this application, it should be noted that, unless otherwise clearly defined, terms such as setting, installing, connecting, etc. should be understood in a broad sense, and technicians in the relevant technical field can reasonably determine the specific meaning of the above terms in this application based on the specific content of the technical solution.

[0094] Carburizing and quenching process is an important part in the heat treatment process of low carbon steel. After carburizing, the surface carbon concentration increases, the hardenability is enhanced, and a high carbon solid is formed on the surface, which can improve the surface hardness and wear resistance of parts and extend the service life of parts. However, the phase change and temperature non-uniformity during carburizing and quenching cause parts to deform, further hindering the improvement of part precision. When the deformation of parts is too large, it will not only increase the trial production cycle of parts, but also reduce the fatigue life of parts during service.

[0095] At present, for the carburizing process, only the effect of carbon content on the material's thermophysical parameters, phase transition kinetics and mechanical properties is considered. In the actual process, carbon atoms dissolve into the gaps in the metal crystal lattice to form a solid solution. Due to the different atomic sizes of the solvent and the solute, the lattice is locally distorted. For interstitial solid solutions, as solute atoms are added, the solid solution lattice constant always increases with the addition of solute atoms. In order to accurately calculate the additional deformation caused by carburizing and quenching, the interstitial solid solution deformation caused by carburizing cannot be ignored. Therefore, the existing numerical models cannot accurately predict carburizing and quenching deformation.

[0096] In order to solve the above technical defects, the embodiments of the present application provide a method, system, device and storage medium for calculating gear carburizing and quenching deformation.

[0097] See also Figure 1 , is a flow chart of a method for calculating gear carburizing and quenching deformation provided by an embodiment of the present application, the method is applied to electronic equipment, and the electronic equipment may be a server, etc. Figure 1 As shown, the gear carburizing and quenching deformation calculation method includes:

[0098] Step S101, obtaining a first node coordinate set, material thermal conductivity and material specific heat capacity of the gear, wherein the first node coordinate set is an initial coordinate set of all nodes after meshing the gear before deformation;

[0099] Step S102, calculating the carbon content of the target period by using Fick's second law based on the first node coordinate set, and performing a test according to the carbon content to obtain a TTT curve;

[0100] Step S103, calculating the temperature of the target period according to the thermal conductivity and specific heat capacity of the material;

[0101] Step S104, calculating the first volume fraction and the second volume fraction of the target time period according to the carbon content, the TTT curve and the temperature;

[0102] Step S105, calculating the plastic strain increment of the target time period according to the first volume fraction and the second volume fraction;

[0103] Step S106, calculating a second node coordinate set of the gear through an Abaqus solver based on the plastic strain increment;

[0104] Step S107: Calculate the deformation value of the gear according to the first node coordinate set and the second node coordinate set.

[0105] The method obtains the first node coordinate set of the gear, the thermal conductivity of the material and the specific heat capacity of the material, wherein the first node coordinate set is the initial coordinate set of all nodes after the gear before deformation is meshed; the carbon content of the target time period is calculated based on the first node coordinate set by Fick's second law, and the TTT curve is obtained by testing according to the carbon content; the temperature of the target time period is calculated according to the thermal conductivity of the material and the specific heat capacity of the material; the first volume fraction and the second volume fraction of the target time period are calculated according to the carbon content, the TTT curve and the temperature; the plastic strain increment of the target time period is calculated according to the first volume fraction and the second volume fraction; based on the plastic strain increment, the second node coordinate set of the gear is calculated by the Abaqus solver; the deformation value of the gear is calculated according to the first node coordinate set and the second node coordinate set. The present application improves the accuracy of the calculation of the carburizing and quenching deformation of the gear by considering the phase change thermoelastic-plasticity of the carburizing expansion strain.

[0106] In some embodiments, in step S102, the carbon content of the target period is calculated based on the first node coordinate set by using Fick's second law, including:

[0107] Step S201, calculating the carbon diffusion coefficient according to the thermal conductivity and specific heat capacity of the material;

[0108] Step S202: Calculate the carbon content of the target period according to the carbon diffusion coefficient and the first node coordinate set using Fick's second law using the following formula:

[0109]

[0110] Where C is the carbon content, t s is the carburizing time corresponding to the target period, x is the X coordinate of the first node coordinate set, y is the Y coordinate of the first node coordinate set, and z is the Z coordinate of the first node coordinate set.

[0111] Specifically, refer to Figure 2 and Figure 3 , taking gears as the research object, the gear mesh division structure is as follows Figure 1As shown, the initial coordinates of each node of the gear after meshing are obtained. Then, the carburizing temperature, carburizing time, specific heat capacity and thermal conductivity of the gear material, and the initial carbon content of carburizing are used as known conditions to input Abaqus' UMATHT for use in a custom subroutine. Then, the Abaqus processor solves the equation of carbon diffusion to obtain the calculated carbon content at each position of the gear.

[0112] The carbon diffusion equation is as follows, satisfying Fick's second law:

[0113]

[0114] Where C is the carbon content, t s is the carburizing time corresponding to the target period, x is the X coordinate of the first node coordinate set, y is the Y coordinate of the first node coordinate set, and z is the Z coordinate of the first node coordinate set.

[0115] In some embodiments, in step S103, calculating the temperature of the target time period according to the thermal conductivity and specific heat capacity of the material includes:

[0116] Step S301, obtaining the released latent heat value in the target period;

[0117] Step S302: Calculate the temperature of the target period according to the released latent heat value, the thermal conductivity of the material and the specific heat capacity of the material by the following formula:

[0118]

[0119] Where T is the temperature, t q is the quenching time corresponding to the target period, ρ is the gear density, λ is the thermal conductivity of the gear material, c is the specific heat capacity of the gear material, is the Laplace operator, and W is the released latent heat.

[0120] Specifically, refer to Figures 4 to 6 The carbon content of each node of the above gear, carburizing time, carburizing temperature, quenching temperature, quenching time, heat transfer coefficient of quenching oil, heat transfer coefficient of air, TTT curve of gear material (gear material performance parameters, which can be obtained by looking up the table), thermal conductivity, specific heat capacity and density of gear material are used as known conditions. The UMATHT of Abaqus is input into the custom subroutine, and then the Abaqus processor solves the temperature field and organization field of the whole process of carburizing and quenching. The final calculation results are as follows: Figures 4 to 6 shown.

[0121] The temperature equation for heat conduction is:

[0122]

[0123] Where T is the temperature, tq is the quenching time corresponding to the target period, ρ is the gear density, λ is the thermal conductivity of the gear material, c is the specific heat capacity of the gear material, is the Laplace operator, and W is the released latent heat.

[0124] In some embodiments, in step S104, calculating the first volume fraction and the second volume fraction of the target time period according to the carbon content, the TTT curve and the temperature includes:

[0125] Step S401, calculating the first moment, the second moment, the first transition amount, and the second transition amount corresponding to the temperature T;

[0126] Step S402: Calculate the first volume fraction of the target period according to the temperature using the following formula:

[0127]

[0128] Where b is the first parameter, n is the second parameter, V is the first volume fraction of the diffusion-type phase transition, and t q is the quenching time corresponding to the target period, t1 is the first moment corresponding to the temperature T in the TTT curve, t2 is the second moment corresponding to the temperature T in the TTT curve, V1 is the first transition value corresponding to t1, and V2 is the second transition value corresponding to t2;

[0129] Step S403, calculating martensite mobility according to temperature and carbon content;

[0130] Step S404, calculating the phase transformation kinetic parameters and the martensite starting transformation temperature value according to the carbon content and the TTT curve;

[0131] Step S405: Calculate the second volume fraction of the target time period according to the martensite mobility, phase transformation kinetic parameters and martensite start transformation temperature using the following formula:

[0132]

[0133] Among them, V M is the second volume fraction of the non-diffusion phase transition, v m (T,C) is the martensite mobility, α is the phase transition kinetic parameter 1, β is the phase transition kinetic parameter 2, M s is the temperature at which martensite transformation begins.

[0134] Specifically, for the calculation of the organizational field, the main focus is on the diffusion-type phase transformation and non-diffusion-type phase transformation of austenite. Taking temperature as input, if the temperature reaches the temperature at which the diffusion-type phase transformation begins, the volume fraction of the diffusion-type phase transformation is calculated:

[0135]

[0136] Among them, V is the first volume fraction of the diffusion-type phase transition, b is the first parameter, and n is the second parameter, which can be obtained through the TTT curve. The expression is as follows:

[0137]

[0138] Among them, t1 is the first moment corresponding to the temperature T, t2 is the second moment corresponding to the temperature T, V1 is the first transition amount, and V2 is the second transition amount.

[0139] If the temperature reaches the temperature at which the non-diffusion phase transformation (i.e., martensitic phase transformation) begins, calculate the volume fraction of the non-diffusion phase transformation:

[0140]

[0141] Among them, V M is the second volume fraction of the non-diffusion phase transition, v m (T,C) is the martensite mobility, α is the phase transition kinetic parameter 1, β is the phase transition kinetic parameter 2, M s is the temperature at which martensite transformation begins.

[0142] In some embodiments, in step S105, calculating the plastic strain increment of the target time period according to the first volume fraction and the second volume fraction includes:

[0143] Step S501, obtaining the total strain increment, stress increment and stress generated by UMATHT;

[0144] Step S502, calculating a third volume fraction of austenite transformation according to the first volume fraction and the second volume fraction;

[0145] Step S503, calculating the elastic strain increment according to the first volume fraction, the second volume fraction, the third volume fraction, the stress increment and the stress by the following formula:

[0146]

[0147] in, is the elastic strain increment, is the inverse matrix of the elasticity matrix, dσ ij is the stress increment, σ is the stress, when k is 1, V k is the third volume fraction. When k is 2, V k is the first volume fraction. When k is 3, V k is the second volume fraction;

[0148] Step S504: Calculate the temperature strain increment according to the first volume fraction, the second volume fraction and the third volume fraction by using the following formula:

[0149]

[0150] in, is the temperature strain increment, T0 is the temperature corresponding to the initial time of the target period, T1 is the temperature corresponding to the end time of the target period, when k is 1, α k is the thermal expansion coefficient of austenite phase transformation. When k is 2, α k is the thermal expansion coefficient of the diffusion phase transition. When k is 3, α k is the thermal expansion coefficient of the non-diffusion phase transition, δ ij is the Kronecker function;

[0151] Step S505: Calculate the phase change strain increment according to the first volume fraction and the second volume fraction using the following formula:

[0152]

[0153] in, is the phase change strain increment. When k is 2, β k is the structural expansion of the diffusion phase transition. When k is 3, β k Structural expansion for non-diffusion phase transition;

[0154] Step S506: Calculate the phase transformation plastic strain increment according to the first volume fraction and the second volume fraction using the following formula:

[0155]

[0156] in, is the phase transformation plastic strain increment, R is the phase transformation plasticity coefficient, S ij is the deviatoric stress tensor;

[0157] Step S507: Calculate the carburizing strain increment according to the carbon content using the following formula:

[0158]

[0159] in, is the carburizing strain increment, α c is the carburizing parameter value corresponding to the carbon content;

[0160] Step S508: Calculate the plastic strain increment of the target period according to the carburizing strain increment, phase transformation plastic strain increment, phase transformation strain increment, temperature strain increment, elastic strain increment and total strain increment by the following formula:

[0161]

[0162] in, is the plastic strain increment, dεij is the total strain increment.

[0163] In some embodiments, in step S106, based on the plastic strain increment, a second node coordinate set of the gear is calculated by an Abaqus solver, including:

[0164] Step S601, based on the plastic strain increment, the carburizing and quenching deformation values ​​of all nodes of the gear are calculated by the Abaqus solver;

[0165] Step S602: Calculate a second node coordinate set according to the carburizing and quenching deformation values ​​of all nodes.

[0166] In some embodiments, in step S107, the deformation value of the gear is calculated by the following formula:

[0167]

[0168] Among them, d is the deformation value of the gear node, x0 is the X coordinate of the first node coordinate, y0 is the Y coordinate of the first node coordinate, z0 is the Z coordinate of the first node coordinate, x1 is the X coordinate of the second node coordinate, y1 is the Y coordinate of the second node coordinate, and z1 is the Z coordinate of the second node coordinate.

[0169] Specifically, refer to Figure 7 , taking the above-obtained organizational field and carburizing field as input, and taking the elastic modulus, stress-strain curve, elasticity, temperature, phase change, phase change plasticity and strain increment caused by carburizing as known conditions, input Abaqus' UMAT for use in the custom subroutine, and then the Abaqus processor solves the strain field of the entire carburizing and quenching process, that is, the carburizing and quenching deformation of the gear. The final deformation result is as follows Figure 7 shown.

[0170] This application can also be used to calculate the carburizing and quenching deformation of various samples.

[0171] Specifically, in some embodiments, the carburizing and quenching process is calculated based on the established carburizing and quenching multi-field quantity calculation module.

[0172] Reference Figure 8 and Fig. 9 , Fig. 9 The schematic diagram of the sample size process route is as follows: A cylindrical sample of φ12.7×38mm is established. The sample material is 8620RH. The sample is oil quenched to 110℃ after strong infiltration, diffusion, and heat preservation, and then air-cooled to room temperature. The cylindrical sample size and heat treatment process route are as follows Fig. 9 The mesh of the surface carburized area of ​​the sample is refined, the mesh thickness is 200mm, the total depth is 3mm, and the rest is a loose mesh area with a mesh thickness of 335mm. The mesh division is as shown in Figure 8 shown.

[0173] Reference Fig.10 The final calculated carburizing result cloud diagram is as follows Fig.10 As shown, the carbon content distribution gradient in different carburizing steps is related to the environmental carbon content. The final carbon content is distributed in a gradient along the cross section. The surface carbon content is 0.75%, and the 0.35% carbon content is an indicator of the effective hardened layer depth. The effective hardened layer depth is about 0.9 mm.

[0174] Reference Fig.11 , Fig.11 The temperature change history of different positions of the carburized quenched cylindrical specimen is shown. The surface cooling rate is the fastest at the beginning of quenching, and the cooling rate is the slowest at the core. The temperature at each position of the cylinder is basically the same around 20s, dropping to 200°C. The temperature rise caused by the latent heat of phase change can be observed in the temperature change over time curve.

[0175] Reference Fig.12 The final calculated tissue distribution cloud diagram is as follows Fig.12 As shown in the figure, after carburizing and quenching, the cylindrical specimen is mainly martensite, with the surface martensite content of about 90.6%, and the subsurface martensite content of up to 98.3%. There is a small amount of residual austenite in the surface layer, which gradually decreases with the increase of layer depth, and the core bainite content is up to about 10.1%.

[0176] Specifically, in order to facilitate understanding by those skilled in the art, a set of optimal embodiments is provided below:

[0177] 1. Data acquisition:

[0178] The first node coordinate set, material thermal conductivity and material specific heat capacity of the gear are obtained, wherein the first node coordinate set is an initial coordinate set of all nodes after meshing the gear before deformation.

[0179] 2. Calculation of carbon content:

[0180] Based on the first node coordinate set, the carbon content of the target period is calculated by Fick's second law, and the TTT curve is obtained by testing according to the carbon content, specifically:

[0181] Calculate the carbon diffusion coefficient based on the thermal conductivity and specific heat capacity of the material;

[0182] Step S202: Calculate the carbon content of the target period according to the carbon diffusion coefficient and the first node coordinate set using Fick's second law using the following formula:

[0183]

[0184] Where C is the carbon content, t sis the carburizing time corresponding to the target period, x is the X coordinate of the first node coordinate set, y is the Y coordinate of the first node coordinate set, and z is the Z coordinate of the first node coordinate set.

[0185] 3. Temperature calculation:

[0186] The temperature of the target period is calculated based on the thermal conductivity and specific heat capacity of the material, specifically:

[0187] Obtaining the released latent heat value during the target period;

[0188] The temperature of the target period is calculated by the following formula based on the released latent heat value, material thermal conductivity and material specific heat capacity:

[0189]

[0190] Where T is the temperature, t q is the quenching time corresponding to the target period, ρ is the gear density, λ is the thermal conductivity of the gear material, c is the specific heat capacity of the gear material, is the Laplace operator, and W is the released latent heat.

[0191] 4. Volume fraction calculation:

[0192] The first volume fraction and the second volume fraction of the target period are calculated according to the carbon content, TTT curve and temperature, specifically:

[0193] Calculate the first moment, the second moment, the first transition amount and the second transition amount corresponding to the temperature T;

[0194] The first volume fraction of the target period is calculated according to the temperature using the following formula:

[0195]

[0196]

[0197] Where b is the first parameter, n is the second parameter, V is the first volume fraction of the diffusion-type phase transition, and t q is the quenching time corresponding to the target period, t1 is the first moment corresponding to the temperature T in the TTT curve, t2 is the second moment corresponding to the temperature T in the TTT curve, V1 is the first transition value corresponding to t1, and V2 is the second transition value corresponding to t2;

[0198] Calculate martensite mobility based on temperature and carbon content;

[0199] Calculate the phase transformation kinetic parameters and martensite start transformation temperature according to carbon content and TTT curve;

[0200] The second volume fraction of the target period is calculated according to the martensite mobility, phase transformation kinetic parameters and martensite start transformation temperature using the following formula:

[0201]

[0202] Among them, V M is the second volume fraction of the non-diffusion phase transition, v m (T,C) is the martensite mobility, α is the phase transition kinetic parameter 1, β is the phase transition kinetic parameter 2, M s is the temperature at which martensite transformation begins.

[0203] 5. Calculation of plastic strain increment:

[0204] The plastic strain increment of the target period is calculated according to the first volume fraction and the second volume fraction, specifically:

[0205] Get the total strain increment, stress increment and stress generated by UMATHT;

[0206] Calculating a third volume fraction of austenite transformation based on the first volume fraction and the second volume fraction;

[0207] The elastic strain increment is calculated according to the first volume fraction, the second volume fraction, the third volume fraction, the stress increment and the stress by the following formula:

[0208]

[0209] in, is the elastic strain increment, is the inverse matrix of the elasticity matrix, dσ ij is the stress increment, σ is the stress, when k is 1, V k is the third volume fraction. When k is 2, V k is the first volume fraction. When k is 3, V k is the second volume fraction;

[0210] The temperature strain increment is calculated according to the first volume fraction, the second volume fraction and the third volume fraction using the following formula:

[0211]

[0212] in, is the temperature strain increment, T0 is the temperature corresponding to the initial time of the target period, T1 is the temperature corresponding to the end time of the target period, when k is 1, α k is the thermal expansion coefficient of austenite phase transformation. When k is 2, α k is the thermal expansion coefficient of the diffusion phase transition. When k is 3, α kis the thermal expansion coefficient of the non-diffusion phase transition, δ ij is the Kronecker function;

[0213] The phase transformation strain increment is calculated according to the first volume fraction and the second volume fraction using the following formula:

[0214]

[0215] in, is the phase change strain increment. When k is 2, β k is the structural expansion of the diffusion phase transition. When k is 3, β k Structural expansion for non-diffusion phase transition;

[0216] The phase transformation plastic strain increment is calculated according to the first volume fraction and the second volume fraction using the following formula:

[0217]

[0218] in, is the phase transformation plastic strain increment, R is the phase transformation plasticity coefficient, S ij is the deviatoric stress tensor;

[0219] The carburizing strain increment is calculated according to the carbon content using the following formula:

[0220]

[0221] in, is the carburizing strain increment, α c is the carburizing parameter value corresponding to the carbon content;

[0222] The plastic strain increment of the target period is calculated by the following formula based on the carburizing strain increment, phase transformation plastic strain increment, phase transformation strain increment, temperature strain increment, elastic strain increment and total strain increment:

[0223]

[0224] in, is the plastic strain increment, dε ij is the total strain increment.

[0225] 6. Calculation of the second node coordinate set:

[0226] Based on the plastic strain increment, the Abaqus solver is used to calculate the second node coordinate set of the gear, which is:

[0227] Based on the plastic strain increment, the carburizing and quenching deformation values ​​of all nodes of the gear are calculated by the Abaqus solver;

[0228] The second node coordinate set is calculated according to the carburizing and quenching deformation values ​​of all nodes.

[0229] 7. Calculation of deformation value:

[0230] The deformation value of the gear is calculated according to the first node coordinate set and the second node coordinate set, wherein the calculation formula for calculating the deformation value of the gear according to the first node coordinate set and the second node coordinate set is:

[0231]

[0232] Among them, d is the deformation value of the gear node, x0 is the X coordinate of the first node coordinate, y0 is the Y coordinate of the first node coordinate, z0 is the Z coordinate of the first node coordinate, x1 is the X coordinate of the second node coordinate, y1 is the Y coordinate of the second node coordinate, and z1 is the Z coordinate of the second node coordinate.

[0233] This application uses carburizing temperature and carburizing time as parameter inputs, based on the material performance parameters (thermal conductivity and specific heat capacity) that vary with temperature and carbon content, and through the Abaqus subroutine interface UMATHT to discretize Fick's second law, calculates the carburizing field distribution under different carburizing process parameters.

[0234] This application is based on the material TTT curve and the calculated carburizing field distribution, and uses the Scheil superposition law to realize the calculation of the quenching process structure field. Based on the current temperature and structure fraction, the linear superposition law is used to calculate the material thermophysical property parameters, and the heat transfer equation is brought into the solution of the temperature field. The carburizing quenching temperature-structure field calculation is realized by combining Fortran with the UMATHT subroutine.

[0235] This application establishes a phase change thermal elastic-plastic constitutive model considering carburizing expansion strain. The carburizing field and temperature field calculated in the above steps are used as input, and the tangent stiffness method is used in the UMATHT subroutine to update the elastic-plastic matrix to realize the carburizing quenching deformation calculation.

[0236] In addition, refer to Fig.13 One embodiment of the present application provides a gear carburizing and quenching deformation calculation system, including a data acquisition module 1100, a carbon content calculation module 1200, a temperature calculation module 1300, a volume fraction calculation module 1400, a plastic strain increment calculation module 1500, a second node coordinate set calculation module 1600 and a deformation value calculation module 1700, wherein:

[0237] The data acquisition module 1100 is used to acquire a first node coordinate set, material thermal conductivity, and material specific heat capacity of the gear, wherein the first node coordinate set is an initial coordinate set of all nodes after meshing the gear before deformation;

[0238] The carbon content calculation module 1200 is used to calculate the carbon content of the target period based on the first node coordinate set by using Fick's second law, and to obtain a TTT curve by testing according to the carbon content;

[0239] The temperature calculation module 1300 is used to calculate the temperature of the target period according to the thermal conductivity and specific heat capacity of the material;

[0240] The volume fraction calculation module 1400 is used to calculate the first volume fraction and the second volume fraction of the target period according to the carbon content, the TTT curve and the temperature;

[0241] The plastic strain increment calculation module 1500 is used to calculate the plastic strain increment of the target time period according to the first volume fraction and the second volume fraction;

[0242] The second node coordinate set calculation module 1600 is used to calculate the second node coordinate set of the gear through the Abaqus solver based on the plastic strain increment;

[0243] The deformation value calculation module 1700 is used to calculate the deformation value of the gear according to the first node coordinate set and the second node coordinate set.

[0244] The system obtains the first node coordinate set of the gear, the thermal conductivity of the material and the specific heat capacity of the material, wherein the first node coordinate set is the initial coordinate set of all nodes after the gear before deformation is meshed; the carbon content of the target time period is calculated based on the first node coordinate set by Fick's second law, and the TTT curve is obtained by testing according to the carbon content; the temperature of the target time period is calculated according to the thermal conductivity of the material and the specific heat capacity of the material; the first volume fraction and the second volume fraction of the target time period are calculated according to the carbon content, the TTT curve and the temperature; the plastic strain increment of the target time period is calculated according to the first volume fraction and the second volume fraction; based on the plastic strain increment, the second node coordinate set of the gear is calculated by the Abaqus solver; the deformation value of the gear is calculated according to the first node coordinate set and the second node coordinate set. The application improves the accuracy of the calculation of the carburizing and quenching deformation of the gear by considering the phase change thermal elastic-plasticity of the carburizing expansion strain.

[0245] It should be noted that the present system embodiment and the above-mentioned method embodiment are based on the same inventive concept, so the relevant contents of the above-mentioned method embodiment are also applicable to the present system embodiment and will not be repeated here.

[0246] Fig.14 A schematic diagram of the hardware structure for calculating gear carburizing and quenching deformation provided in an embodiment of the present application is shown.

[0247] The gear carburizing and quenching deformation calculation device may include a processor 901 and a memory 902 storing computer program instructions.

[0248] Specifically, the processor 901 may include a central processing unit (CPU), or an application specific integrated circuit (ASIC), or may be configured to implement one or more integrated circuits of the embodiments of the present application.

[0249] The memory 902 may include a large capacity memory for data or instructions. By way of example and not limitation, the memory 902 may include a hard disk drive (HDD), a floppy disk drive, a flash memory, an optical disk, a magneto-optical disk, a magnetic tape, or a universal serial bus (USB) drive or a combination of two or more of these. In appropriate cases, the memory 902 may include a removable or non-removable (or fixed) medium. In appropriate cases, the memory 902 may be inside or outside the integrated gateway disaster recovery device. In a specific embodiment, the memory 902 is a non-volatile solid-state memory.

[0250] In some embodiments, the memory 902 may include a read-only memory (ROM), a random access memory (RAM), a magnetic disk storage media device, an optical storage media device, a flash memory device, an electrical, optical or other physical / tangible memory storage device. Thus, in general, the memory includes one or more tangible (non-transitory) computer-readable storage media (e.g., a memory device) encoded with software including computer-executable instructions, and when the software is executed (e.g., by one or more processors), it is operable to perform the operations described with reference to the method according to an aspect of the present disclosure.

[0251] The processor 901 reads and executes the computer program instructions stored in the memory 902 to implement any one of the gear carburizing and quenching deformation calculation methods in the above embodiments.

[0252] In one example, the gear carburizing and quenching deformation calculation device may further include a communication interface 903 and a bus 910. Fig.14 As shown, the processor 901, the memory 902, and the communication interface 903 are connected via a bus 910 and communicate with each other.

[0253] The communication interface 903 is mainly used to implement communication between various modules, devices, units and / or equipment in the embodiments of the present application.

[0254] Bus 910 includes hardware, software or both, and the parts of gear carburizing quenching deformation computing equipment are coupled to each other.For example, but not limitation, bus may include accelerated graphics port (AGP) or other graphics bus, enhanced industrial standard architecture (EISA) bus, front side bus (FSB), hypertransmission (HT) interconnection, industrial standard architecture (ISA) bus, infinite bandwidth interconnection, low pin count (LPC) bus, memory bus, micro channel architecture (MCA) bus, peripheral component interconnection (PCI) bus, PCI-Express (PCI-X) bus, serial advanced technology attachment (SATA) bus, video electronics standard association local (VLB) bus or other suitable bus or two or more of these combinations. In appropriate cases, bus 910 may include one or more buses. Although the present application embodiment describes and shows a specific bus, the application considers any suitable bus or interconnection.

[0255] The gear carburizing and quenching deformation calculation device can execute the gear carburizing and quenching deformation calculation method in the embodiment of the present application based on the three-dimensional design model, thereby realizing the combination Figure 1 and Fig.13 The invention describes a method and system for calculating gear carburizing and quenching deformation.

[0256] In addition, in combination with the gear carburizing and quenching deformation calculation method in the above embodiment, the embodiment of the present application can provide a computer storage medium for implementation. The computer storage medium stores computer program instructions; when the computer program instructions are executed by the processor, any one of the gear carburizing and quenching deformation calculation methods in the above embodiment is implemented.

[0257] It should be clear that the present application is not limited to the specific configuration and processing described above and shown in the figures. For the sake of simplicity, a detailed description of the known method is omitted here. In the above embodiments, several specific steps are described and shown as examples. However, the method process of the present application is not limited to the specific steps described and shown, and those skilled in the art can make various changes, modifications and additions, or change the order between the steps after understanding the spirit of the present application.

[0258] The functional blocks shown in the above-described block diagram can be implemented as hardware, software, firmware or a combination thereof. When implemented in hardware, it can be, for example, an electronic circuit, an application specific integrated circuit (ASIC), appropriate firmware, a plug-in, a function card, etc. When implemented in software, the elements of the present application are programs or code segments that are used to perform the required tasks. The program or code segment can be stored in a machine-readable medium, or transmitted on a transmission medium or a communication link by a data signal carried in a carrier wave. "Machine-readable medium" can include any medium capable of storing or transmitting information. Examples of machine-readable media include electronic circuits, semiconductor memory devices, ROM, flash memory, erasable ROM (EROM), floppy disks, CD-ROMs, optical disks, hard disks, optical fiber media, radio frequency (RF) links, etc. The code segment can be downloaded via a computer network such as the Internet, an intranet, etc.

[0259] It should also be noted that the exemplary embodiments mentioned in this application describe some methods or systems based on a series of steps or devices. However, this application is not limited to the order of the above steps, that is, the steps can be performed in the order mentioned in the embodiment, or in a different order from the embodiment, or several steps can be performed simultaneously.

[0260] Aspects of the present disclosure are described above with reference to the flowchart and / or block diagram of the method, device (system) and computer program product according to the embodiment of the present disclosure. It should be understood that each box in the flowchart and / or block diagram and the combination of each box in the flowchart and / or block diagram can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device to produce a machine so that these instructions executed by the processor of the computer or other programmable data processing device enable the implementation of the function / action specified in one or more boxes of the flowchart and / or block diagram. Such a processor can be, but is not limited to, a general-purpose processor, a special-purpose processor, a special application processor, or a field programmable logic circuit. It can also be understood that each box in the block diagram and / or flowchart and the combination of boxes in the block diagram and / or flowchart can also be implemented by dedicated hardware that performs a specified function or action, or can be implemented by a combination of dedicated hardware and computer instructions.

[0261] The above is only a specific implementation of the present application. Those skilled in the art can clearly understand that for the convenience and simplicity of description, the specific working processes of the systems, modules and units described above can refer to the corresponding processes in the aforementioned method embodiments, and will not be repeated here. It should be understood that the protection scope of the present application is not limited to this. Any technician familiar with the technical field can easily think of various equivalent modifications or replacements within the technical scope disclosed in this application, and these modifications or replacements should be included in the protection scope of this application.

Claims

1. A method for calculating gear carburizing and quenching deformation, characterized in that: The gear carburizing and quenching deformation calculation method includes: Obtaining a first node coordinate set, material thermal conductivity, and material specific heat capacity of the gear, wherein the first node coordinate set is an initial coordinate set of all nodes after meshing the gear before deformation; Calculating the carbon content of the target period by using Fick's second law based on the first node coordinate set, and performing a test according to the carbon content to obtain a TTT curve; The temperature of the target time period is calculated according to the thermal conductivity of the material and the specific heat capacity of the material, specifically: Obtaining the released latent heat value during the target period; The temperature of the target period is calculated according to the released latent heat value, the thermal conductivity of the material and the specific heat capacity of the material by the following formula: Where T is the temperature, t q is the quenching time corresponding to the target period, ρ is the gear density, λ is the thermal conductivity of the gear material, c is the specific heat capacity of the gear material, is the Laplace operator, W is the released latent heat value; The first volume fraction and the second volume fraction of the target time period are calculated according to the carbon content, the TTT curve and the temperature, specifically: Calculate the first moment, the second moment, the first transition amount and the second transition amount corresponding to the temperature T; The first volume fraction of the target time period is calculated according to the temperature by the following formula: Where b is the first parameter, n is the second parameter, V is the first volume fraction of the diffusion-type phase transition, and t q is the quenching time corresponding to the target period, t1 is the first moment corresponding to the temperature T in the TTT curve, t2 is the second moment corresponding to the temperature T in the TTT curve, V1 is the first transition value corresponding to t1, and V2 is the second transition value corresponding to t2; calculating martensite mobility based on the temperature and the carbon content; Calculate phase transformation kinetic parameters and martensite start transformation temperature value according to the carbon content and the TTT curve; The second volume fraction of the target time period is calculated according to the martensite mobility, the phase transformation kinetic parameter and the martensite start transformation temperature value by the following formula: Among them, V M is the second volume fraction of the non-diffusion phase transition, v m (T,C) is the martensite mobility, α is the phase transition kinetic parameter 1, β is the phase transition kinetic parameter 2, M s is the martensite transformation starting temperature; The plastic strain increment of the target time period is calculated according to the first volume fraction and the second volume fraction, specifically: Get the total strain increment, stress increment and stress generated by UMATHT; Calculating a third volume fraction of austenite transformation according to the first volume fraction and the second volume fraction; The elastic strain increment is calculated according to the first volume fraction, the second volume fraction, the third volume fraction, the stress increment and the stress by the following formula: in, is the elastic strain increment, is the inverse matrix of the elasticity matrix, dσ ij is the stress increment, σ is the stress, when k is 1, V k is the third volume fraction. When k is 2, V k is the first volume fraction. When k is 3, V k is the second volume fraction; The temperature strain increment is calculated according to the first volume fraction, the second volume fraction and the third volume fraction by the following formula: in, is the temperature strain increment, T0 is the temperature corresponding to the initial time of the target period, T1 is the temperature corresponding to the end time of the target period, when k is 1, α k is the thermal expansion coefficient of austenite phase transformation. When k is 2, α k is the thermal expansion coefficient of the diffusion phase transition. When k is 3, α k is the thermal expansion coefficient of the non-diffusion phase transition, δ ij is the Kronecker function; The phase change strain increment is calculated according to the first volume fraction and the second volume fraction by the following formula: in, is the phase change strain increment. When k is 2, β k is the structural expansion of the diffusion phase transition. When k is 3, β k Structural expansion for non-diffusion phase transition; The phase transformation plastic strain increment is calculated according to the first volume fraction and the second volume fraction by the following formula: in, is the phase transformation plastic strain increment, R is the phase transformation plasticity coefficient, S ij is the deviatoric stress tensor; The carburizing strain increment is calculated according to the carbon content by the following formula: in, is the carburizing strain increment, α c is the carburizing parameter value corresponding to the carbon content; The plastic strain increment of the target time period is calculated according to the carburizing strain increment, the phase transformation plastic strain increment, the phase transformation strain increment, the temperature strain increment, the elastic strain increment and the total strain increment by the following formula: in, is the plastic strain increment, dε ij is the total strain increment; Based on the plastic strain increment, calculating a second node coordinate set of the gear by an Abaqus solver; The deformation value of the gear is calculated according to the first node coordinate set and the second node coordinate set.

2. The method for calculating gear carburizing and quenching deformation according to claim 1, characterized in that: The calculating the carbon content of the target period by using Fick's second law based on the first node coordinate set includes: Calculating the carbon diffusion coefficient according to the thermal conductivity of the material and the specific heat capacity of the material; The carbon content of the target period is calculated according to the carbon diffusion coefficient and the first node coordinate set by Fick's second law using the following formula: Where C is the carbon content, t s is the carburizing time corresponding to the target period, x is the X coordinate of the first node coordinate set, y is the Y coordinate of the first node coordinate set, and z is the Z coordinate of the first node coordinate set.

3. The method for calculating gear carburizing and quenching deformation according to claim 2, characterized in that: The step of calculating the second node coordinate set of the gear by using an Abaqus solver based on the plastic strain increment includes: Based on the plastic strain increment, the carburizing and quenching deformation values ​​of all nodes of the gear are calculated by the Abaqus solver; The second node coordinate set is calculated according to the carburizing and quenching deformation values ​​of all the nodes.

4. The method for calculating gear carburizing and quenching deformation according to claim 1, characterized in that: The deformation value of the gear is calculated by the following formula: Among them, d is the deformation value of the gear node, x0 is the X coordinate of the first node coordinate, y0 is the Y coordinate of the first node coordinate, z0 is the Z coordinate of the first node coordinate, x1 is the X coordinate of the second node coordinate, y1 is the Y coordinate of the second node coordinate, and z1 is the Z coordinate of the second node coordinate.

5. A gear carburizing and quenching deformation calculation system, characterized in that: The gear carburizing and quenching deformation calculation system comprises: A data acquisition module, used to acquire a first node coordinate set, material thermal conductivity and material specific heat capacity of the gear, wherein the first node coordinate set is an initial coordinate set of all nodes after meshing the gear before deformation; A carbon content calculation module, used to calculate the carbon content of a target period by using Fick's second law based on the first node coordinate set, and to obtain a TTT curve by testing according to the carbon content; The temperature calculation module is used to calculate the temperature of the target time period according to the thermal conductivity of the material and the specific heat capacity of the material, specifically: Obtaining the released latent heat value during the target period; The temperature of the target period is calculated according to the released latent heat value, the thermal conductivity of the material and the specific heat capacity of the material by the following formula: Where T is the temperature, t q is the quenching time corresponding to the target period, ρ is the gear density, λ is the thermal conductivity of the gear material, c is the specific heat capacity of the gear material, is the Laplace operator, W is the released latent heat value; A volume fraction calculation module is used to calculate the first volume fraction and the second volume fraction of the target time period according to the carbon content, the TTT curve and the temperature, specifically: Calculate the first moment, the second moment, the first transition amount and the second transition amount corresponding to the temperature T; The first volume fraction of the target time period is calculated according to the temperature by the following formula: Where b is the first parameter, n is the second parameter, V is the first volume fraction of the diffusion-type phase transition, and t q is the quenching time corresponding to the target period, t1 is the first moment corresponding to the temperature T in the TTT curve, t2 is the second moment corresponding to the temperature T in the TTT curve, V1 is the first transition value corresponding to t1, and V2 is the second transition value corresponding to t2; calculating martensite mobility based on the temperature and the carbon content; Calculate phase transformation kinetic parameters and martensite start transformation temperature value according to the carbon content and the TTT curve; The second volume fraction of the target time period is calculated according to the martensite mobility, the phase transformation kinetic parameter and the martensite start transformation temperature value by the following formula: Among them, V M is the second volume fraction of the non-diffusion phase transition, v m (T,C) is the martensite mobility, α is the phase transition kinetic parameter 1, β is the phase transition kinetic parameter 2, M s is the martensite transformation starting temperature; A plastic strain increment calculation module is used to calculate the plastic strain increment of a target time period according to the first volume fraction and the second volume fraction, specifically: Get the total strain increment, stress increment and stress generated by UMATHT; Calculating a third volume fraction of austenite transformation according to the first volume fraction and the second volume fraction; The elastic strain increment is calculated according to the first volume fraction, the second volume fraction, the third volume fraction, the stress increment and the stress by the following formula: in, is the elastic strain increment, is the inverse matrix of the elasticity matrix, dσ ij is the stress increment, σ is the stress, when k is 1, V k is the third volume fraction. When k is 2, V k is the first volume fraction. When k is 3, V k is the second volume fraction; The temperature strain increment is calculated according to the first volume fraction, the second volume fraction and the third volume fraction by the following formula: in, is the temperature strain increment, T0 is the temperature corresponding to the initial time of the target period, T1 is the temperature corresponding to the end time of the target period, when k is 1, α k is the thermal expansion coefficient of austenite phase transformation. When k is 2, α k is the thermal expansion coefficient of the diffusion phase transition. When k is 3, α k is the thermal expansion coefficient of the non-diffusion phase transition, δ ij is the Kronecker function; The phase change strain increment is calculated according to the first volume fraction and the second volume fraction by the following formula: in, is the phase change strain increment. When k is 2, β k is the structural expansion of the diffusion phase transition. When k is 3, β k Structural expansion for non-diffusion phase transition; The phase transformation plastic strain increment is calculated according to the first volume fraction and the second volume fraction by the following formula: in, is the phase transformation plastic strain increment, R is the phase transformation plasticity coefficient, S ij is the deviatoric stress tensor; The carburizing strain increment is calculated according to the carbon content by the following formula: in, is the carburizing strain increment, α c is the carburizing parameter value corresponding to the carbon content; The plastic strain increment of the target time period is calculated according to the carburizing strain increment, the phase transformation plastic strain increment, the phase transformation strain increment, the temperature strain increment, the elastic strain increment and the total strain increment by the following formula: in, is the plastic strain increment, dε ij is the total strain increment; A second node coordinate set calculation module, used for calculating a second node coordinate set of the gear through an Abaqus solver based on the plastic strain increment; A deformation value calculation module is used to calculate the deformation value of the gear according to the first node coordinate set and the second node coordinate set.

6. A gear carburizing and quenching deformation calculation device, characterized in that: comprising at least one control processor and a memory for communicatively coupling with the at least one control processor; The memory stores instructions that can be executed by the at least one control processor, and the instructions are executed by the at least one control processor so that the at least one control processor can execute the gear carburizing and quenching deformation calculation method as described in any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer-executable instructions, and the computer-executable instructions are used to enable a computer to execute a gear carburizing and quenching deformation calculation method as described in any one of claims 1 to 4.

Images