METHOD AND SYSTEM FOR CALCULATING STORED ENERGY FIELD OF PRIMARY SHEAR ZONE DURING STEADY-STATE CUTTING

Abstract

A method and system for calculating a primary shear zone stored energy field during steady-state cutting, the method including: fitting parameters of a workpiece material stored energy evolution model; discretizing the primary shear zone into infinitesimals on a main shear plane. The infinitesimals are small enough, a strain, strain rate, and temperature are assumed constant; introducing an equivalent cutting edge model simplifying three-dimensional cutting into two-dimensional cutting, calculating element strain and strain rate using a shear plane model, and analyzing element temperature using a heat conduction equation; deriving a differential equation of stored energy versus location in the primary shear zone using stored energy evolution, strain rate distribution, strain distribution, and temperature distribution models; and solving the differential equation for each infinitesimal using an initial shear plane as a model boundary, obtaining stored energy at each location to obtain a stored energy field distribution of the primary shear zone.

Claims

1. A method for calculating a stored energy field of a primary shear zone during steady-state cutting, the method comprising steps of: fitting parameters of a stored energy evolution model of a workpiece material; performing infinitesimal division on the primary shear zone; simplifying actual three-dimensional cutting into two-dimensional cutting, performing analysis to obtain a shear plane model, calculating a strain and a strain rate of each infinitesimal, and analyzing a temperature of the each infinitesimal; deriving a differential equation of stored energy versus location by using the stored energy evolution model, a strain rate distribution model, a strain distribution model, and a temperature distribution model; and solving the differential equation of stored energy versus location for the each infinitesimal by using an initial shear plane of the primary shear zone as a model boundary, to obtain stored energy at each location, so as to obtain a stored energy field distribution of the primary shear zone.

2. The method for calculating a stored energy field of the primary shear zone during steady-state cutting according to claim 1, wherein during the fitting of the parameters of the stored energy evolution model of the workpiece material, fitting parameters of a stored energy evolution model of the workpiece material about to the temperature, the strain, and the strain rate based on stress-strain curves of the workpiece material in different deformation conditions.

3. The method for calculating a stored energy field of the primary shear zone during steady-state cutting according to claim 1, wherein after the infinitesimal division on the primary shear zone, the strain, the strain rate, and the temperature in the each infinitesimal are set as constants.

4. The method for calculating a stored energy field of the primary shear zone during steady-state cutting according to claim 1, wherein before calculation of the strain and the strain rate of the each infinitesimal and analysis of the temperature of the each infinitesimal, an equivalent cutting edge model is introduced to simplify the actual three-dimensional cutting into the two-dimensional cutting, the strain and the strain rate of the each element are calculated according to the shear plane model, and the temperature of the each element is analyzed according to a heat conduction equation.

5. The method for calculating a stored energy field of the primary shear zone during steady-state cutting according to claim 1, wherein during the analysis of the temperature of the each infinitesimal, according to the heat conduction equation, a temperature value of a K.sup.th plane is represented by a temperature value of a (K−1).sup.th plane.

6. The method for calculating a stored energy field of the primary shear zone during steady-state cutting according to claim 1, wherein a line connecting two end points of an actual cutting edge projected on a base plane is defined as an equivalent cutting edge, an equivalent angle of the equivalent cutting edge is calculated by using a cutting tool angle of the equivalent cutting edge as an actual cutting tool angle, and then a strain distribution and a strain rate distribution of the primary shear zone are calculated by using a normal rake angle of the equivalent cutting edge as an input parameter of the shear plane model.

7. The method for calculating a stored energy field of the primary shear zone during steady-state cutting according to claim 1, wherein during calculation of stored energy of each discrete plane of the primary shear zone, stored energy of a K.sup.th plane is obtained by integration of stored energy of a (K−1).sup.th plane.

8. The method for calculating a stored energy field of the primary shear zone during steady-state cutting according to claim 1, wherein stored energy-based shear stress field prediction is based on a mapping relationship between stored energy and a dislocation density and a mapping relationship between a shear flow stress and a dislocation density.

9. A system for calculating a stored energy field of a primary shear zone during steady-state cutting, the system comprising: a fitting unit, configured to fit parameters of a stored energy evolution model of a workpiece material; an infinitesimal generation unit, configured to perform infinitesimal division on the primary shear zone; a conversion unit, configured to simplify actual three-dimensional cutting into two-dimensional cutting; and a solving module, configured to receive data outputted by the fitting unit, the infinitesimal generation unit, and the conversion unit, calculate a strain and a strain rate of each infinitesimal and analyze a temperature of the each infinitesimal according to the data outputted by the conversion unit, derive a differential equation of stored energy versus location of the primary shear zone by using the stored energy evolution model of the fitting unit, and solve the differential equation of stored energy versus location for the each infinitesimal divided by the infinitesimal generation unit, to obtain a stored energy field distribution of the primary shear zone.

10. The system for calculating the stored energy field of the primary shear zone during steady-state cutting according to claim 9, wherein the infinitesimals are a series of discrete infinitesimals generated in a normal direction of the main shear plane on the primary shear zone.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0020] The accompanying drawings constituting a part of the present invention are used to provide further understanding of the present invention. The exemplary examples of the present invention and descriptions thereof are used to explain the present invention, and do not constitute an improper limitation to the present invention.

[0021] FIG. 1 is a schematic diagram of a shear zone model according to one or more implementations of the present invention.

[0022] FIG. 2 is a block diagram of a stored energy field calculation program according to one or more implementations of the present invention.

[0023] FIG. 3 is a schematic diagram of stored energy field distributions of the primary shear zone in the conditions of different cutting parameters according to one or more implementations of the present invention.

[0024] FIG. 4 is a schematic diagram of a stress field prediction result of the primary shear zone according to one or more implementations of the present invention.

[0025] FIG. 5 is a schematic diagram of a temperature field prediction result of the primary shear zone according to one or more implementations of the present invention.

[0026] The spacing or dimensions between each part are exaggerated to show the position of each part, and the schematic diagrams are used only for illustrative purposes.

DETAILED DESCRIPTION

[0027] It should be pointed out that the following detailed descriptions are all illustrative and are intended to provide further descriptions of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meanings as those usually understood by a person of ordinary skill in the art to which the present invention belongs.

[0028] It should be noted that the terms used herein are merely used for describing specific implementations, and are not intended to limit exemplary implementations of the present invention. As used herein, the singular form is also intended to include the plural form unless the present invention clearly dictates otherwise. In addition, it should be further understood that, terms “comprise” and/or “include” used in this specification indicate that there are features, steps, operations, devices, components, and/or combinations thereof.

[0029] For convenience of description, the terms “above”, “below”, “left”, and “right” only indicate directions consistent with those of the accompanying drawings, are not intended to limit the structure, and are used only for ease and brevity of illustration and description, rather than indicating or implying that the mentioned device or element needs to have a particular orientation or needs to be constructed and operated in a particular orientation. Therefore, such terms should not be construed as a limitation on the present invention.

[0030] For the part of term explanation, terms in the present invention such as “mount”, “connect”, “connection”, and “fix” should be understood in a broad sense. For example, the connection may be a fixed connection, a detachable connection, or an integral connection, a mechanical connection, an electrical connection, a direct connection, an indirect connection by using an intermediate medium, an interior connection between two components, or interaction between two components. A person of ordinary skill in the art may understand specific meanings of the foregoing terms in the present invention according to a specific situation.

[0031] As described in the background, in the traditional research of a cutting mechanism, the stress field has directionality, excessive characteristic parameters are present, and the research process is complex. In view of the shortcomings of the prior art, the present invention provides a method and system for calculating a stored energy field of the primary shear zone during steady-state cutting, to predict a cutting force, a cutting temperature, a chip morphology, and material properties by using a stored energy distribution of the primary shear zone.

Embodiment 1

[0032] In a typical implementation of the present invention, this embodiment discloses a method for calculating a stored energy field of the primary shear zone during steady-state cutting.

[0033] The method includes the following steps:

[0034] (1) Fit parameters of a stored energy evolution model of a workpiece material based on stress-strain curves of the workpiece material in different deformation conditions, where the model is related to the temperature, strain, and strain rate. Stored energy E.sub.s may be represented by a dislocation density ρ.sub.total. That is to say,

E.sub.s=αμb.sup.2ρ.sub.total/χ.

[0035] α is a dislocation interaction parameter having a value of 0.5, μ is a shear modulus of the material, b is a value of a Burgers vector, and χ is a proportionality coefficient allowing for alloying elements and having a value of 0.6. Therefore, the establishment of the stored energy evolution model depends on a dislocation density evolution model. According to a dislocation density type, the dislocation density includes a statistical stored dislocation density ρ.sub.s and a geometrically necessary dislocation density ρ.sub.LABs. Therefore, the dislocation density is:

ρ.sub.total=(1−ƒ)ρ.sub.s+ƒρ.sub.LABs.

[0036] ƒ is a volume fraction of a geometrically necessary dislocation, which may be expressed by a dislocation cell structure diameter D.sub.cell and the dislocation cell wall thickness δ as:

ƒ=[(D.sub.cell−δ)/D.sub.cell].sup.3.

[0037] The dislocation cell wall thickness δ is 1.28×10.sup.−9 m, and the dislocation cell structure diameter is expressed as:

D.sub.cell=k.sub.cell/√{square root over (ρ.sub.total)}.

[0038] k.sub.cell is a material constant. Evolution equations of the statistical stored dislocation density ρ.sub.s and the geometrically necessary dislocation density ρ.sub.LABs with a strain γ are:

dρ.sub.s/dγ=(bD.sub.cell).sup.−1√{square root over (ρs)}−rρ.sub.s

dρ.sub.LABs/dγ=√{square root over (2ρm/3)}D.sub.cell/δb.

[0039] r is a recovery coefficient, which is obtained by fitting results of stress-strain curves at different strain rates and different temperatures, and ρ.sub.m is a mobile dislocation density.

[0040] (2) As shown in FIG. 1, a zone CDFE is the primary shear zone, CD is an initial shear plane, EF is a final shear plane, and AB is a main shear plane. In order to reduce the complexity of the analysis, the primary shear zone is discretized into N infinitesimals from CD to EF in a normal direction of AB, that is, N+1 planes. When N is large enough, a strain rate and a temperature in each infinitesimal may be assumed as constants.

[0041] (3) Introduce an equivalent cutting edge model, as shown in FIG. 2, define two ends of an actual cutting tool as an equivalent cutting edge, simplify actual three-dimensional cutting into two-dimensional cutting, and calculate distributions of a strain rate {dot over (γ)} and a strain γ of the primary shear zone according to a shear plane model:

[00001]$\begin{array}{c}\stackrel{.}{\gamma }=\{\begin{array}{cc}\frac{{\stackrel{.}{\gamma }}_{\max }}{\mathrm{ks}}y& 0<y<\mathrm{ks}\\ \frac{{\stackrel{.}{\gamma }}_{\max }}{s-\mathrm{ks}}\left(s-y\right)& \mathrm{ks}<y<s\end{array}\\ \gamma =\{\begin{array}{cc}\frac{{\stackrel{.}{\gamma }}_{\max }}{2\mathrm{ksV}\sin {\varphi }_e}{y}^2& 0<y<\mathrm{ks}\\ \frac{{\stackrel{.}{\gamma }}_{\max }}{2\left(1-k\right)\mathrm{sV}\sin {\varphi }_e}\left(2\mathrm{sy}-y^2-2\mathrm{ks}\right)& \mathrm{ks}\le y<s\end{array}.\end{array}$

[0042] s is a thickness of the primary shear zone, k is a ratio of a distance from CD to AB to a distance from AB to EF, y is a distance from a point in a cutting zone to CD, {dot over (γ)}.sub.max is a maximum strain rate in the primary shear zone, ϕ.sub.e is a shear angle, and V is a cutting speed. {dot over (γ)}.sub.max and k may be expressed as:

[00002]$\begin{array}{c}{\stackrel{.}{\gamma }}_{\max }=\frac{2V\cos {\gamma }_e}{\sqrt{3}s\cos \left({\varphi }_e-{\gamma }_e\right)}\\ k=\frac{\cos {\varphi }_e\cos \left({\varphi }_e-{\gamma }_e\right)}{\cos {\gamma }_e}.\end{array}$

[0043] A thickness of the shear zone is obtained according to the empirical formula of oxley:

[00003]$s=\frac{a_c}{5.9\sin {\varphi }_e}.$

[0044] In order to obtain the shear angle, experiments are usually required to obtain a deformation coefficient. For simplicity, in this embodiment, an approximate shear angle formula of Merchant is used:

[00004]${\varphi }_e=\frac{\pi }{4}-\frac{\beta }{2}+\frac{{\gamma }_e}{2}.$

[0045] β is a friction angle. The model does not allow for the influence of the cutting speed on a friction coefficient. Therefore, β is a constant.

[0046] According to the distribution laws of the strain and the strain rate, the formula is substituted into a center position of the each infinitesimal, to obtain an average strain and an average strain rate of the each infinitesimal. The average strain and the average strain rate are respectively used as a strain feature value and a strain rate feature value of the each infinitesimal.

[0047] According to the heat conduction equation, the temperature value of the K.sup.th plane is represented by a temperature value of a (K−1).sup.th plane:

[00005]$T_K=T_{K-1}+\frac{R_1}{\rho \mathrm{cV}\sin {\varphi }_e}{\int }_{y_{K-1}}^{y_K}\mathrm{Qdy}.$

[0048] R.sub.1 is a ratio of mass transfer to heat transfer during the cutting, ρ and c are respectively a density and a specific heat capacity of the workpiece material, y.sub.K and y.sub.K-1 are respectively coordinates of the K.sup.th plane and the (K−1).sup.th plane, and Q is heat generated per unit time and per unit volume in the each infinitesimal.

[0049] (4) Derive a differential equation of stored energy versus location y by using the stored energy evolution model:

[00006]$\frac{{\mathrm{dE}}_s}{\mathrm{dy}}=\frac{{\mathrm{dE}}_{s}d\gamma }{d\gamma \mathrm{dy}}+\frac{{\mathrm{dE}}_{s}d\stackrel{.}{\gamma }}{d\stackrel{.}{\gamma }\mathrm{dy}}+\frac{{\mathrm{dE}}_s\mathrm{dT}}{\mathrm{dTdy}}.$

[0050] Since the strain rate and the temperature in the each infinitesimal are both regarded as constants, the differential equation in the each infinitesimal may be simplified. Finally, stored energy E.sub.s|K of the K.sup.th plane is calculated by stored energy E.sub.s|K-1 of the (K−1).sup.th plane:

[00007]$E_{s.Math.K}=E_{s.Math.K-1}+\frac{d\psi \left({\stackrel{.}{\gamma }}_{\mathrm{AVE}},\stackrel{\_}{T}\right)d\gamma }{d\gamma \mathrm{dy}}\mathrm{dy}.$

[0051] The initial shear plane of the primary shear zone is used as a model boundary. Stored energy of a next plane is calculated in the each infinitesimal according to the foregoing formula, to obtain stored energy of all (N+1).sup.th planes. The stored energy of all planes is used as the stored energy of the location in the primary shear zone. That is to say, the stored energy field distribution of the primary shear zone is obtained. The specific stored energy calculation process is shown in the process block diagram in FIG. 2.

[0052] (5) Predict a stress field and a temperature field of the primary shear zone based on the stored energy field, and then analyze the cutting force, the cutting temperature, a chip forming law, and the material modification.

[0053] According to the foregoing technical solutions, the coefficient k.sub.cell in step (1) may be obtained in two ways. In the first way, the coefficient may be estimated as μ/200. In the second way, the dislocation density ρ.sub.total and the dislocation cell diameter D.sub.cell after the machining are measured by experiments, and then k.sub.cell=D.sub.cell√{square root over (ρ.sub.total)} is calculated.

[0054] The recovery factor r is a function of the strain rate and the temperature, which is expressed as:

[00008]$r=\left[a\exp \left(-\frac{\stackrel{.}{\gamma }}{{\stackrel{.}{\gamma }}_0}\right)+b\right]\exp \left(\mathrm{mT}\right).$

[0055] A, B, and m are all parameters obtained by fitting the stress-strain curve.

[0056] According to the foregoing technical solutions, the primary shear zone in the steady-cutting process is discretized into N infinitesimals in a normal direction of the main shear plane, and a strain rate and a temperature in each infinitesimal are regarded as constants to simplify a solving process of a differential equation.

[0057] According to the foregoing technical solutions, a line connecting two end points of an actual cutting edge projected on a base plane is defined as an equivalent cutting edge, a cutting tool angle of the equivalent cutting edge is calculated according to a geometric relationship and is used as an actual cutting tool angle, and then a strain distribution and a strain rate distribution of the primary shear zone are calculated by using a normal rake angle of the equivalent cutting edge as an input parameter of the shear plane model.

[0058] According to the foregoing technical solutions, stored energy of the discrete planes in the primary shear zone in step (4) may be obtained by solving the following mathematical physical problem:

[0059] Boundary conditions are as follows: the initial shear plane is used as a model boundary, that is, an 0.sup.th plane in the (N+1).sup.th planes obtained by division, and it is assumed according to the experimental results and the models that a temperature of the initial shear plane is a room temperature (25° C.), a strain and a strain rate are both 0, and stored energy is stored energy of an initial material.

[0060] A mathematical physical equation is:

[00009]$E_{s.Math.K}=E_{s.Math.K-1}+\frac{d\psi \left({\stackrel{.}{\gamma }}_{\mathrm{AVE}},\stackrel{\_}{T}\right)d\gamma }{d\gamma \mathrm{dy}}\mathrm{dy}.$

[0061] The equation may be explained as follows: The stored energy of the K.sup.th plane may be obtained by integration of the stored energy of the (K−1).sup.th plane, where {dot over (γ)}.sub.AVE and T are respectively a feature strain rate and a feature temperature of a K.sup.th element.

[0062] According to the foregoing technical solutions, the stored energy-based shear stress field prediction in step (5) is based on a mapping relationship between stored energy and a dislocation density and a mapping relationship between a shear flow stress and a dislocation density. The mapping relationship between a shear stress and a dislocation density may be expressed as:

[00010]$\tau =\frac{1}{\sqrt{3}}\left({\tau }_0+M{\alpha }_s\mu b\sqrt{{\rho }_f}\right)+{\tau }_1\mu {\left\{1-{\left[\frac{k_{b}T}{\Delta F}\ln \left(\frac{{\stackrel{.}{\epsilon }}_p}{{\stackrel{.}{\epsilon }}_0}\right)\right]}^{\frac{1}{p}}\right\}}^{\frac{1}{q}}).$

[0063] In the formula, the first term τ.sub.0 is a non-thermal stress independent of the strain rate during movement of a dislocation, the second term is a long-range stress of the movement of the dislocation, and the third term is the short-range stress during the movement of the dislocation. M is a Taylor coefficient, α.sub.s is a constant generally having a value of 0.3-0.5, τ.sub.1μ is a stress value required for the dislocation to cross an obstacle without the assistance of thermal activation, where a value thereof mainly depends on a strength of the obstacle during the movement of the dislocation, the constants p and q satisfy 0<p≤1 and 1≤q≤2, and values thereof depend on a shape of the energy barrier, ΔF=ƒ.sub.2μb.sup.3 is Helmholtz free energy required by the dislocation to overcome the short-range barrier without the assistance of an external force, k.sub.b is a Boltzmann constant, T is an absolute temperature, and {dot over (ε)}.sub.p and {dot over (ε)}.sub.0 are respectively a plastic strain rate and a reference strain rate. According to the mapping relationship between stored energy and a dislocation density, a shear flow stress prediction formula based on the stored energy is obtained as:

[00011]$\tau ={\tau }_0+M\left(\mathrm{\alpha \mu }b\sqrt{E_s\chi /{\alpha }_s\mu b^2}\right)+{\tau }_1\mu {\left\{1-{\left[\frac{k_{b}T}{\Delta F}\ln \left(\frac{{\stackrel{.}{\epsilon }}_p}{{\stackrel{.}{\epsilon }}_0}\right)\right]}^{\frac{1}{p}}\right\}}^{\frac{1}{q}}.$

[0064] B is a lattice damping coefficient, and ρ.sub.m is a mobile dislocation density.

[0065] According to the foregoing technical solutions, a normal stress c perpendicular to the main shear plane is calculated based on a Hencky equation:

[00012]$\frac{\partial \kappa }{\partial y_e}+2\tau \frac{\partial \psi }{\partial y_e}-\frac{\partial \tau }{\partial x_e}=0.$

[0066] ψ is a turn of a shear line. Therefore, by means of integration, a distribution of κ can be obtained.

[0067] A normal stress κ.sub.A at a free surface may be estimated by a stress at the free surface. Therefore,

[00013]${\kappa }_A=\tau \left[1+2\left(\frac{\pi }{4}-{\varphi }_e\right)\right].$

[0068] A chip morphology prediction method is as follows: It is assumed the loading stage exists before the main shear plane and an unloading stage exists after the main shear plane AB. If a stored energy peak occurs at a location before AB, it may be considered that adiabatic shear occurs in the cutting condition and sawtooth chips are formed.

[0069] According to the foregoing technical solutions, a microhardness prediction formula based on the stored energy is:

[00014]$\mathrm{Hv}=\frac{{\sigma }_0+M\mathrm{\alpha \mu }b\sqrt{\chi E_s/\mathrm{\alpha \mu }{b}^2}}{C_H}.$

[0070] In the formula, σ.sub.0=Mτ.sub.0, M is a Taylor coefficient having a value of 3.06, and C.sub.H is a material constant.

[0071] A microscopic residual stress may be directly represented by the stored energy.

[0072] This embodiment fills up the gap in the existing cutting mechanism research technology. Based on the shear plane model in the traditional research of the cutting mechanism, the equivalent cutting edge model was first introduced to simplify three-dimensional cutting into two-dimensional cutting. Then, the established stored energy evolution model of the machined material and the strain rate distribution model, the strain distribution model, and the temperature distribution model of the primary shear zone during the cutting are substituted into the shear plane model, and the differential equation of stored energy versus location in the primary shear zone is obtained to obtain the stored energy distribution of the primary shear zone. Finally, the stress field and the temperature field of the primary shear zone are analyzed based on the obtained stored energy distribution, and the cutting force, the cutting temperature, and the material modification are further predicted. The stored energy field runs through the entire cutting process, and the cutting mechanism is explained more deeply, completely, and clearly in a simpler way, greatly promoting the research of the cutting mechanism.

[0073] For example, the workpiece material is nickel-based alloy inconel718, the cutting tool is a cemented carbide cutting tool, and parameters of the cutting tool are as follows: a rake angle of 0°, a tool cutting edge angle of 90°, a tool cutting edge inclination of 0°, and a tool tip arc radius of 0.8 mm. The following cutting parameters are used:

TABLE-US-00001 Group No.: Cutting speed (m/min) Feed (mm/r) Depth of cut (mm) 1 18 0.1 1.0 2 24 0.1 1.0 3 30 0.1 1.0 4 36 0.1 1.0 5 42 0.1 1.0

[0074] A method for calculating a stored energy field of the primary shear zone during steady-state cutting of the nickel-based alloy inconel718 includes the following steps:

[0075] 1) Fit parameters of a stored energy evolution model of the material based on stress-strain curves of the workpiece material in different deformation conditions, where the model is related to a temperature, a strain, and a strain rate;

[0076] Stored energy E.sub.s may be represented by a dislocation density ρ.sub.total. That is to say,

E.sub.s=αμb.sup.2ρ.sub.total/χ.

[0077] For the inconel718 workpiece material, α=0.5, μ=80000 MPa, b=2.56e-7 mm, and χ=0.6. A total dislocation density is

ρ.sub.total=(1-f)ρ.sub.s+fρ.sub.LABs.

[0078] f is a volume fraction of a geometrically necessary dislocation, which may be expressed by a dislocation cell structure diameter D.sub.cell and the dislocation cell wall thickness δ as:

f=[(D.sub.cell−δ)/D.sub.cell].sup.3.

[0079] The dislocation cell wall thickness δ is 1.28×10.sup.−9 m, and the dislocation cell structure diameter is expressed as:

D.sub.cell=k.sub.cell/√{square root over (ρ.sub.total)}.

[0080] k.sub.cell=6.4. Evolution equations of the statistical stored dislocation density ρ.sub.s and the geometrically necessary dislocation density ρ.sub.LABs with a strain γ are:

dρ.sub.s/dγ=(bD.sub.cell).sup.−1√{square root over (ρs)}−rρ.sub.s

dρ.sub.LABs/dγ=√{square root over (2ρm/3)}D.sub.cell/δb.

[0081] By fitting results of stress-strain curves at different strain rates and different temperatures, the following is obtained:

[00015]$\begin{array}{c}r=\left[87.84\exp \left(-\frac{\stackrel{.}{\gamma }}{2040}\right)+16.65\right]\exp \left(0.00045T\right)\\ {\rho }_m=2.4\times {10}^6{\mathrm{mm}}^{-2}.\end{array}$

[0082] 2) In order to reduce the complexity of the analysis, the primary shear zone is discretized into N infinitesimals from an initial shear plane CD to a final shear plane EF in a normal direction of a main shear plane AB, that is, N+1 planes, where when N is large enough, a strain rate and a temperature in each infinitesimal may be assumed as constants, and in consideration of calculation efficiency and calculation accuracy, N=200.

[0083] 3) Introduce an equivalent cutting edge model, simplify actual three-dimensional cutting into two-dimensional cutting, and calculate distributions of a strain rate {dot over (γ)} and a strain γ of the primary shear zone according to a shear plane model:

[00016]$\begin{array}{c}\stackrel{.}{\gamma }=\{\begin{array}{cc}\frac{{\stackrel{.}{\gamma }}_{\max }}{\mathrm{ks}}y& 0<y<\mathrm{ks}\\ \frac{{\stackrel{.}{\gamma }}_{\max }}{s-\mathrm{ks}}\left(s-y\right)& \mathrm{ks}<y<s\end{array}\\ \gamma =\{\begin{array}{cc}\frac{{\stackrel{.}{\gamma }}_{\max }}{2\mathrm{ksV}\sin {\varphi }_e}{y}^2& 0<y<\mathrm{ks}\\ \frac{{\stackrel{.}{\gamma }}_{\max }}{2\left(1-k\right)\mathrm{sV}\sin {\varphi }_e}\left(2\mathrm{sy}-y^2-2\mathrm{ks}\right)& \mathrm{ks}\le y<s\end{array}.\end{array}$

{dot over (γ)}max and k may be expressed as:

[00017]$\begin{array}{c}{\stackrel{.}{\gamma }}_{\max }=\frac{2V\cos {\gamma }_e}{\sqrt{3}s\cos \left({\varphi }_e-{\gamma }_e\right)}\\ k=\frac{\cos {\varphi }_e\cos \left({\varphi }_e-{\gamma }_e\right)}{\cos {\gamma }_e}.\end{array}$

[0084] A thickness of the shear zone is obtained according to the empirical formula of oxley:

[00018]$s=\frac{a_c}{5.9\sin {\varphi }_e}.$

[0085] In order to obtain the shear angle, experiments are usually required to obtain a deformation coefficient. For simplicity, in this embodiment, an approximate shear angle formula of Merchant is used:

[00019]${\varphi }_e=\frac{\pi }{4}-\frac{\beta }{2}+\frac{{\gamma }_e}{2}.$

[0086] In the formula, β=arctan (f.sub.μ) is a friction angle, and f.sub.μ=0.35. A cutting speed V, an equivalent rake angle γ.sub.e, a cutting thickness a.sub.c, and β are substituted into the foregoing calculation formula, to obtain the strain distribution and the strain rate distribution of the primary shear zone.

[0087] According to the distribution laws of the strain and the strain rate, the formula is substituted into a center position of the each infinitesimal, to obtain an average strain and an average strain rate of the each infinitesimal. The average strain and the average strain rate are respectively used as a strain feature value and a strain rate feature value of the each infinitesimal.

[0088] According to the heat conduction equation, the temperature value of a K.sup.th plane is represented by a temperature value of a (K−1).sup.th plane:

[00020]$T_K=T_{K-1}+\frac{R_1}{\rho cV\sin {\varphi }_e}{\int }_{y_{K-1}}^{y_K}\mathrm{Qdy}.$

[0089] R.sub.1 is a ratio of mass transfer to heat transfer during the cutting, p and c are respectively a density and a specific heat capacity of the workpiece material, γ.sub.K and γ.sub.K-1 are respectively coordinates of the K.sup.th plane and the (K−1).sup.th plane, and Q is heat generated per unit time and per unit volume in the each infinitesimal.

[0090] 4) Derive a differential equation of stored energy versus location y by using the stored energy evolution model:

[00021]$\frac{{\mathrm{dE}}_s}{\mathrm{dy}}=\frac{{\mathrm{dE}}_{s}d\gamma }{d\gamma \mathrm{dy}}+\frac{{\mathrm{dE}}_{s}d\stackrel{.}{\gamma }}{d\stackrel{.}{\gamma }\mathrm{dy}}+\frac{{\mathrm{dE}}_s\mathrm{dT}}{\mathrm{dTdy}}.$

[0091] Since the strain rate and the temperature in the each infinitesimal are both regarded as fixed values, the differential equation in the each infinitesimal may be simplified. Finally, stored energy E.sub.s|K of the K.sup.th plane is calculated by stored energy E.sub.s|K of the (K−1).sup.th plane:

[00022]$E_{s.Math.K}=E_{s.Math.K-1}+\frac{d\psi \left({\stackrel{.}{\gamma }}_{\mathrm{AVE}},\stackrel{\_}{T}\right)d\gamma }{d\gamma \mathrm{dy}}\mathrm{dy}.$

[0092] The initial shear plane of the primary shear zone is used as a model boundary. Stored energy of a next plane is calculated in the each infinitesimal according to the foregoing formula, to obtain stored energy of all (N+1).sup.th planes. The stored energy of all planes is used as the stored energy of the location in the primary shear zone. That is to say, the stored energy field distribution of the primary shear zone is obtained. FIG. 3 shows a stored energy field result obtained according to the cutting parameters in the table.

[0093] 5) Predict a stress field and a temperature field of the primary shear zone based on the stored energy field, and then analyze the cutting force, the cutting temperature, a chip forming law, and the material modification.

[0094] A shear flow stress prediction formula based on the stored energy is obtained as:

[00023]$\tau ={\tau }_0+M\left(\mathrm{\alpha \mu }b\sqrt{E_s\chi /{\alpha }_s\mu b^2}\right)+{\tau }_1\mu {\left\{1-{\left[\frac{k_{b}T}{\Delta F}\ln \left(\frac{{\stackrel{.}{\epsilon }}_p}{{\stackrel{.}{\epsilon }}_0}\right)\right]}^{\frac{1}{p}}\right\}}^{\frac{1}{q}}.$

[0095] A stress field distribution of the primary shear zone is predicted according to the stored energy field result, as shown in FIG. 4. A temperature field distribution of the primary shear zone is predicted according to the stored energy field result, as shown in FIG. 5.

[0096] A predicted stress is used as an input and is substituted into a main cutting force prediction formula. A comparison between a predicted result and an experimental result of the cutting force is shown in the table. It may be learned that the predicted result is quite similar with the experimental result.

TABLE-US-00002 Group No.: Predicted value (N) Experimental value (N) [00024]$\frac{.Math.\mathrm{Predicted}\mathrm{value}-\mathrm{Experimental}\mathrm{value}.Math.}{\mathrm{Experimental}\mathrm{value}}\times 100\%$ 1 482.343581 579.1666667 16.72% 2 443.7292291 557.5980392 20.42% 3 441.7968373 546.8137255 19.21% 4 427.8047855 536.0294118 20.19% 5 419.3042459 533.3333333 21.38%

[0097] A chip morphology prediction method is as follows: It is assumed the loading stage exists before the main shear plane and an unloading stage exists after the main shear plane AB. If a stored energy peak occurs at a location before AB, it may be considered that adiabatic shear occurs in the cutting condition and sawtooth chips are formed. In consideration of a hysteresis effect of the adiabatic shear, a critical speed at which the sawtooth chips are generated is about 30 m/min, which is quite in line with the experimental result. Hardness variations of the cut shear zone may be obtained based on a mapping relationship between microhardness and stored energy, and a microscopic stored energy distribution of the cut shear zone may be further obtained according to the stored energy distribution

Embodiment 2

[0098] In a typical implementation of the present invention, this embodiment discloses a system for calculating a stored energy field of the primary shear zone during steady-state cutting.

[0099] The system includes:

[0100] a fitting unit, an infinitesimal generation unit, and a conversion unit, where the fitting unit is configured to fit a stored energy evolution model of a material, the infinitesimal generation unit is configured to divide the primary shear zone into infinitesimals, and the conversion unit is configured to convert a three-dimensional cutting model to a two-dimensional cutting model; and

[0101] a solving module, configured to receive data outputted by the fitting unit, the infinitesimal generation unit, and the conversion unit, calculate a strain and a strain rate of each infinitesimal and analyze the temperature of the each infinitesimal according to the data outputted by the conversion unit, derive a differential equation of stored energy versus location of the primary shear zone by using the stored energy evolution model of the fitting unit, and solve the differential equation of stored energy versus location for the each infinitesimal divided by the infinitesimal generation unit, to obtain a stored energy field distribution of the primary shear zone.

[0102] The infinitesimals are a series of discrete infinitesimals generated in a normal direction of the main shear plane on the primary shear zone.

[0103] It may be understood that, a pre-processing module may be an existing processor. The processor is connected to a memory in which the program codes of the fitting unit, the infinitesimal generation unit, and the conversion unit in this embodiment are burned. Alternatively, the pre-processing module includes three processors. The three processors each are connected to a memory in which the program codes of the fitting unit, the infinitesimal generation unit, and the conversion unit in this embodiment are burned. An output conversion module is further connected to the pre-processing module for reading, analyzing, organizing, assembling, converting, and drawing results outputted by the pre-processing module, and mapping various attributes to inherent names and attributes.

[0104] The solution unit is configured to establish a differential equation and a definite solution condition for a mathematical engineering model according to the mathematical engineering model, a numerical discrete algorithm, and a numerical solution method that are preselected, perform calculation of discretized regions for continuous time physical quantities and continuous space physical quantities, and establish an algebraic equation by the numerical solution method to form a solution result.

[0105] It may be understood that the solution unit may be an existing processor.

[0106] It may be understood that, in this embodiment, the system described in Embodiment 1 may be used for calculation. Specifically, in the pre-processing module, the fitting unit is configured to perform step (1), the infinitesimal generation unit is configured to perform step (2), the conversion unit is configured to perform step (3), and the solving module is configured to perform step (4) and step (5).

[0107] The above descriptions are merely preferred embodiments of the present invention and are not intended to limit the present invention. A person skilled in the art may make various alterations and variations to the present invention. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention shall fall within the protection scope of the present invention.