METHOD FOR STEADY-STATE CONTROL OF CUTTING STATE

Abstract

A method for steady-state control of cutting state, implemented based on computer numerical control (CNC) machine tools or cutting lathes. The method comprises: fitting parameters of a stored energy evolution model of a workpiece material to be machined; discretizing primary shear zone into multiple infinitesimals along normal direction of main shear plane; introducing equivalent cutting edge model, inputting pre-used cutting parameters, calculating strain and strain rate of each infinitesimal and analyzing temperature of each infinitesimal; deducing and solving differential equation of each infinitesimal of stored energy to position of the primary shear zone by taking initial shear plane of the primary shear zone as model boundary; determining application values of cutting parameters according to solved results; and, controlling and adjusting the cutting tool in the actual cutting process to cut the workpiece material to be machined with the application values.

Claims

1. A method for steady-state control of cutting state, being implemented based on a computer numerical control (CNC) machine tool or a cutting lathe, and comprising steps of: fitting parameters of a stored energy evolution model of a workpiece material to be machined based on stress-strain curves of the workpiece material to be machined in different deformation conditions; performing infinitesimal division on the primary shear zone; simplifying an actual three-dimensional (3D) cutting into a two-dimensional (2D) cutting, obtaining a shear plane model by analyzing the 2D cutting; inputting multiple groups of data of cutting speed and cutting depth pre-used in an actual cutting process into the shear plane model, calculating a strain and a strain rate of each of infinitesimals, and analyzing temperature of the each of the infinitesimals; 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; solving the differential equation of stored energy versus location for the each of the infinitesimals 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; analyzing the stored energy field distribution of the primary shear zone, predicting a cutting force, cutting temperature, and a chip morphology during the actual cutting process; when the predicted cutting force, cutting temperature, and/or chip morphology obtained from a certain prediction calculation meet the predefined steady-state requirements of the cutting state, the pre-used data of the cutting speed and cutting depth input in the certain prediction calculation are used as application values of steady-state control of a cutting tool in the actual cutting process; and controlling and adjusting a cutting state of the cutting tool in the actual cutting process to process the workpiece material to be machined by using the application values of the cutting speed and the cutting depth, so as to maintain the actual cutting process in a steady-state and obtain expected material properties of a machined workpiece material.

2. The method for steady-state control of cutting state according to claim 1, wherein fitting of the parameters of the stored energy evolution model of the workpiece material to be machined comprises: fitting the parameters of the stored energy evolution model of the workpiece material to be machined about to a temperature, a strain, and a strain rate based on the stress-strain curves of the workpiece material to be machined in the different deformation conditions.

3. The method for steady-state control of cutting state 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 steady-state control of cutting state according to claim 1, wherein before calculation of the strain and the strain rate of the each of the infinitesimals and analysis of the temperature of the each of the infinitesimals, an equivalent cutting edge model is introduced to simplify the actual three-dimensional (3D) cutting into the two-dimensional (2D) cutting, the strain and the strain rate of the each of the infinitesimals are calculated according to the shear plane model, and the temperature of the each of the infinitesimals is analyzed according to a heat conduction equation.

5. The method for steady-state control of cutting state according to claim 1, wherein during the analysis of the temperature of the each of the infinitesimals, according to the heat conduction equation, a temperature value of a K.sup.th plane is represented by a temperature value of a (K1).sup.th plane.

6. The method for steady-state control of cutting state according to claim 1, wherein a line connecting two end points of an cutting edge of the cutting tool in the actual cutting process 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 a cutting angle of the cutting tool in the actual cutting process, 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 steady-state control of cutting state 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 (K1).sup.th plane.

8. The method for steady-state control of cutting state 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. The method for steady-state control of cutting state according to claim 1, wherein controlling and adjusting the cutting state of the cutting tool in the actual cutting process to process the workpiece material to be machined by using the application values of the cutting speed and the cutting depth comprises: inputting the application values of the cutting speed and cutting depth into the CNC machine tool, or operating corresponding control levers on the cutting lathe by an operator according to the application values of the cutting speed and cutting depth, to control and adjust the cutting tool to perform a steady-state machining on the workpiece material to be machined according to the input application values during the actual cutting process, and obtain the expected material properties of the machined workpiece material.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0031] 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.

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

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

[0034] 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.

[0035] 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.

[0036] 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.

[0037] FIG. 6 is a schematic diagram of a method for steady-state control of cutting state according to one or more implementations of the present invention.

[0038] 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

[0039] 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.

[0040] 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 include and/or including used in this specification indicate that there are features, steps, operations, devices, components, and/or combinations thereof.

[0041] 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.

[0042] 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.

[0043] 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.

Example 1

[0044] In a typical implementation of the present invention, the present example discloses a method for calculating a stored energy field of the primary shear zone during steady-state cutting. The method includes the following steps:

[0045] (1) Fitting 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:

[00001]$E_s=b^2{}_{\mathrm{total}}/,$ [0046] wherein, a is a dislocation interaction parameter having a value of 0.5, u 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:

[00002]${}_{\mathrm{total}}=\left(1-f\right){}_s+f{}_{L\mathrm{ABs}},$ [0047] wherein, 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:

[00003]$f={\left[\left(D_{cell}-\right)/D_{cell}\right]}^3.$

[0048] The dislocation cell wall thickness is 1.2810.sup.9 m, and the dislocation cell structure diameter is expressed as:

[00004]$D_{\mathrm{cell}}=k_{\mathrm{cell}}/\sqrt{{}_{\mathrm{total}}},$ [0049] wherein, 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:

[00005]$d{}_s/d={\left(b{D}_{\mathrm{cell}}\right)}^{-1}\sqrt{{}_s}-r{}_s,$$d{}_{\mathrm{LABs}}/d=\sqrt{2{}_m/3}{D}_{\mathrm{cell}}/b,$ [0050] wherein, 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.

[0051] (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.

[0052] (3) Introducing 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 3D cutting into 2D cutting, and calculate distributions of a strain rate {dot over ()} and a strain of the primary shear zone according to a shear plane model:

[00006]$\stackrel{}{}=\{\begin{array}{cc}\frac{{\stackrel{}{}}_{\max }}{\mathrm{ks}}y& 0<y<\mathrm{ks}\\ \frac{{\stackrel{}{}}_{\max }}{s-\mathrm{ks}}\left(s-y\right)& \mathrm{ks}<y<s\end{array},$$=\{\begin{array}{cc}\frac{{\stackrel{}{}}_{\max }}{2\mathrm{ksV}\sin {}_e}{y}^2& 0<y<\mathrm{ks}\\ \frac{{\stackrel{}{}}_{\max }}{2\left(1-k\right)\mathrm{sV}\sin {}_e}\left(2\mathrm{sy}-y^2-2\mathrm{ks}\right)& \mathrm{ks}y<s\end{array},$ [0053] wherein, 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:

[00007]${\stackrel{}{}}_{\max }=\frac{2V\cos {}_e}{\sqrt{3}s\cos \left({}_e-{}_e\right)},$$k=\frac{\cos {}_e\cos \left({}_e-{}_e\right)}{\cos {}_e}.$

[0054] A thickness of the shear zone is obtained according to the empirical formula of Oxley:

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

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

[00009]${}_e=\frac{}{4}-\frac{}{2}+\frac{{}_e}{2},$ [0056] wherein, is a friction angle. The model does not allow for the influence of the cutting speed on a friction coefficient. Therefore, is a constant.

[0057] 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.

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

[00010]$T_K=T_{K-1}+\frac{R_1}{\mathrm{cV}\sin {}_e}{}_{y_{K-1}}^{y_K}\mathrm{Qdy},$ [0059] wherein, 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.K1 are respectively coordinates of the K.sup.th plane and the (K1).sup.th plane, and Q is heat generated per unit time and per unit volume in the each infinitesimal.

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

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

[0061] 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|K1 of the (K1).sup.th plane:

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

[0062] 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.

[0063] (5) Predicting 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.

[0064] 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.

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

[00013]$r=\left[\mathrm{Aexp}\left(-\frac{\stackrel{}{}}{{\stackrel{}{}}_0}\right)+B\right]\exp \left(\mathrm{mT}\right),$ [0066] wherein, A, B, and m are all parameters obtained by fitting the stress-strain curve.

[0067] 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.

[0068] 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.

[0069] 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:

[0070] 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.

[0071] A mathematical physical equation is:

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

[0072] 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 (K1).sup.th plane, wherein {dot over ()}.sub.AVE and T are respectively a feature strain rate and a feature temperature of a K.sup.th element.

[0073] 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:

[00015]$=\frac{1}{\sqrt{3}}\left({}_0+M\left({}_{s}b\sqrt{{}_f}\right)+{}_1{\left\{1-{\left[\frac{k_{b}T}{F}\ln \left(\frac{{\stackrel{.}{}}_p}{{\stackrel{}{}}_0}\right)\right]}^{\frac{1}{p}}\right\}}^{\frac{1}{q}}\right).$

[0074] In the formula, the first term to 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, as is a constant generally having a value of 0.3-0.5, Tiu 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<p1 and 1q2, and values thereof depend on a shape of the energy barrier, F=f.sub.2b.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:

[00016]$={}_0+M\left(b\sqrt{E_s/{}_s{b}^2}\right)+{}_1{\left\{1-{\left[\frac{k_{b}T}{F}\ln \left(\frac{{\stackrel{.}{}}_p}{{\stackrel{}{}}_0}\right)\right]}^{\frac{1}{p}}\right\}}^{\frac{1}{q}},$ [0075] wherein, B is a lattice damping coefficient, and .sub.m is a mobile dislocation density.

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

[00017]$\frac{}{y_e}+2\frac{}{y_e}-\frac{}{x_e}=0,$ [0077] wherein, is a turn of a shear line. Therefore, by means of integration, a distribution of can be obtained.

[0078] A normal stress K.sub.A at a free surface may be estimated by a stress at the free surface. Therefore:

[00018]${}_A=\left[1+2\left(\frac{}{4}-{}_e\right)\right].$

[0079] 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.

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

[00019]$H_v=\frac{{}_0+Mb\sqrt{E_s/b^2}}{C_H},$ [0081] wherein, .sub.0=M.sub.0, M is a Taylor coefficient having a value of 3.06, and C.sub.H is a material constant.

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

[0083] The present example 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.

[0084] 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 Cutting Depth Group speed Feed of cut No.: (m/min) (mm/r) (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

[0085] 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:

[0086] 1) Fitting 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;

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

[00020]$E_s=b^2{}_{\mathrm{total}}/.$

[0088] For the inconel718 workpiece material, =0.5, =80000 MPa, b=2.56e7 mm, and =0.6. A total dislocation density is

[00021]${}_{\mathrm{total}}=\left(1-f\right){}_s+f{}_{\mathrm{LABs}},$ [0089] wherein, 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 8 as:

[00022]$f={\left[\left(D_{\mathrm{cell}}-\right)/D_{\mathrm{cell}}\right]}^3.$

[0090] The dislocation cell wall thickness is 1.2810.sup.9 m, and the dislocation cell structure diameter is expressed as:

[00023]$D_{\mathrm{cell}}=k_{\mathrm{cell}}/\sqrt{{}_{\mathrm{total}}},$ [0091] wherein, 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:

[00024]$d{}_s/d={\left({\mathrm{bD}}_{\mathrm{cell}}\right)}^{-1}\sqrt{{}_s}-r{}_s,$$d{}_{\mathrm{LABs}}/d=\sqrt{2{}_m/3}{D}_{\mathrm{cell}}/b.$

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

[00025]$r=\left[87.84\exp \left(-\frac{\stackrel{.}{}}{2040}\right)+16.65\right]\exp \left(0.00045T\right),$${}_m=2.4{10}^6{\mathrm{mm}}^{-2}.$

[0093] 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.

[0094] 3) Introducing an equivalent cutting edge model, simplify actual 3D cutting into 2D cutting, and calculate distributions of a strain rate j and a strain y of the primary shear zone according to a shear plane model:

[00026]$\stackrel{.}{}=\{\begin{array}{cc}\frac{{\stackrel{.}{}}_{\max }}{\mathrm{ks}}y& 0<y<\mathrm{ks}\\ \frac{{\stackrel{.}{}}_{\max }}{s-\mathrm{ks}}\left(s-y\right)& \mathrm{ks}<y<s\end{array},$$=\{\begin{array}{cc}\frac{{\stackrel{.}{}}_{\max }}{2\mathrm{ksV}\sin {}_e}{y}^2& 0<y<\mathrm{ks}\\ \frac{{\stackrel{.}{}}_{\max }}{2\left(1-k\right)\mathrm{sV}\sin {}_e}\left(2\mathrm{sy}-y^2-2\mathrm{ks}\right)& \mathrm{ks}y<s\end{array},$ [0095] wherein, {acute over ()}.sub.max and k may be expressed as:

[00027]${\stackrel{.}{}}_{\max }=\frac{2V\cos {}_e}{\sqrt{3}s\cos \left({}_e-{}_e\right)},$$k=\frac{\cos {}_e\cos \left({}_e-{}_e\right)}{\cos {}_e}.$

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

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

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

[00029]${}_e=\frac{}{4}-\frac{}{2}+\frac{{}_e}{2},$ [0098] wherein, 32 arc tan(f) 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.

[0099] 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.

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

[00030]$T_K=T_{K-1}+\frac{R_1}{\mathrm{cV}\sin {}_e}{}_{y_{K-1}}^{y_K}\mathrm{Qdy},$ [0101] wherein, 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.k1 are respectively coordinates of the K.sup.th plane and the (K1).sup.th plane, and Q is heat generated per unit time and per unit volume in the each infinitesimal.

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

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

[0103] 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|K1 of the (K1) th plane:

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

[0104] 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.

[0105] 5) Predicting 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.

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

[00033]$={}_0+M\left(b\sqrt{E_s/{}_s{b}^2}\right)+{}_1{\left\{1-{\left[\frac{k_{b}T}{F}\ln \left(\frac{{\stackrel{.}{}}_p}{{\stackrel{.}{}}_0}\right)\right]}^{\frac{1}{p}}\right\}}^{\frac{1}{q}},$

[0107] 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.

[0108] 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) [00034]$\frac{.Math.\mathrm{Predicted}\mathrm{value}-\mathrm{Experimental}\mathrm{value}.Math.}{\mathrm{Experimental}\mathrm{value}}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%

[0109] 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 micro-hardness and stored energy, and a microscopic stored energy distribution of the cut shear zone may be further obtained according to the stored energy distribution.

Example 2

[0110] In a typical implementation of the present invention, the present example discloses a system for calculating a stored energy field of the primary shear zone during steady-state cutting. The system includes: [0111] 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 [0112] 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.

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

[0114] 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 the present example 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 the present example 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.

[0115] 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.

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

[0117] It may be understood that, in the present example, 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).

Example 3

[0118] In a typical implementation of the present invention, the present example discloses a method for steady-state control of cutting state, including the following steps: [0119] fitting parameters of a stored energy evolution model of a workpiece material to be machined based on stress-strain curves of the workpiece material to be machined in different deformation conditions; [0120] performing infinitesimal division on the primary shear zone; [0121] simplifying an actual three-dimensional (3D) cutting into a two-dimensional (2D) cutting, obtaining a shear plane model by analyzing the 2D cutting; [0122] inputting multiple groups of data of cutting speed and cutting depth pre-used in an actual cutting process into the shear plane model, calculating a strain and a strain rate of each of infinitesimals, and analyzing temperature of the each of the infinitesimals; [0123] 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; [0124] solving the differential equation of stored energy versus location for the each of the infinitesimals 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; [0125] analyzing the stored energy field distribution of the primary shear zone, predicting a cutting force, cutting temperature, and a chip morphology during the actual cutting process; [0126] when the predicted cutting force, cutting temperature, and/or chip morphology obtained from a certain prediction calculation meet the predefined steady-state requirements of the cutting state, the pre-used data of the cutting speed and cutting depth input in the certain prediction calculation are used as application values of steady-state control of a cutting tool in the actual cutting process; and [0127] controlling and adjusting a cutting state of the cutting tool in the actual cutting process to process the workpiece material to be machined by using the application values of the cutting speed and the cutting depth, so as to maintain the actual cutting process in a steady-state and obtain expected material properties of a machined workpiece material.

[0128] The method for the steady-state control of the cutting state in an actual cutting process, specially comprises the following steps:

[0129] 1) Fitting 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;

[0130] Stored energy E.sub.s may be represented by a dislocation density .sub.total, which is:

[00035]$E_s=b^2{}_{\mathrm{total}}/.$

[0131] Use a workpiece material of inconel718 as an example.

[0132] For the inconel718 workpiece material, =0.5, =80000 MPa, b=2.56e7 mm, and =0.6. A total dislocation density is

[00036]${}_{\mathrm{total}}=\left(1-f\right){}_s+f{}_{\mathrm{LABs}},$ [0133] wherein, 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:

[00037]$f={\left[\left(D_{\mathrm{cell}}-\right)/D_{\mathrm{cell}}\right]}^3.$

[0134] The dislocation cell wall thickness is 1.2810.sup.9 m, and the dislocation cell structure diameter is expressed as:

[00038]$D_{\mathrm{cell}}=k_{\mathrm{cell}}/\sqrt{{}_{\mathrm{total}}},$ [0135] wherein, k.sub.cell=6.4. Evolution equations of the statistical stored dislocation density Ps and the geometrically necessary dislocation density .sub.LABs with a strain y are:

[00039]$d{}_s/d={\left(b{D}_{\mathrm{cell}}\right)}^{-1}\sqrt{{}_s}-r{}_s,d{}_{\mathrm{LABs}}/d=\sqrt{2{}_m/3}{D}_{\mathrm{cell}}/b.$

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

[00040]$r=\left[87\text{.84}\exp \left(-\frac{\stackrel{}{}}{2040}\right)+16.65\right]\exp \left(0.00045T\right),{}_m=2.410^6{\mathrm{mm}}^{-2}.$

[0137] 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.

[0138] 3) Introducing an equivalent cutting edge model, simplify actual 3D cutting into 2D cutting, and calculate distributions of a strain rate j and a strain y of the primary shear zone according to a shear plane model:

[00041]$\stackrel{}{}=\{\begin{array}{cc}\frac{{\stackrel{.}{}}_{\mathrm{ma}x}}{\mathrm{ks}}y& 0<y<\mathrm{ks}\\ \frac{{\stackrel{.}{}}_{\mathrm{ma}x}}{s-\mathrm{ks}}\left(s-y\right)& \mathrm{ks}<y<s\end{array},=\{\begin{array}{cc}\frac{{\stackrel{.}{}}_{\mathrm{ma}x}}{2\mathrm{ksV}\sin {}_e}{y}^2& 0<y<\mathrm{ks}\\ \frac{{\stackrel{.}{}}_{\mathrm{ma}x}}{2\left(1-k\right)\mathrm{sV}\sin {}_e}\left(2\mathrm{sy}-y^2-2\mathrm{ks}\right)& \mathrm{ks}y<s\end{array},$ [0139] wherein, {dot over ()}.sub.max and k may be expressed as:

[00042]${\stackrel{}{}}_{\mathrm{ma}x}=\frac{2V\cos {}_e}{\sqrt{3}s\cos \left({}_e-{}_e\right)},k=\frac{\cos {}_e\cos \left({}_e-{}_e\right)}{\cos {}_e}.$

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

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

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

[00044]${}_e=\frac{}{4}-\frac{}{2}+\frac{}{2},$ [0142] wherein, B=arc tan(f.sub.) is a friction angle, and f.sub.=0.35. A cutting speed V, an equivalent rake angle .sub.e, a cutting thickness .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.

[0143] 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.

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

[00045]$T_K=T_{K-1}+\frac{R_1}{\mathrm{cV}\sin {}_e}{}_{y_{K-1}}^{y_K}\mathrm{Qdy},$ [0145] wherein, 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.K1 are respectively coordinates of the K.sup.th plane and the (K1).sup.th plane, and Q is heat generated per unit time and per unit volume in the each infinitesimal.

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

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

[0147] 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|K1 of the (K1).sup.th plane:

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

[0148] 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 stored energy fields result obtained according to the cutting parameters in the table.

[0149] 5) Predicting 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.

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

[00048]$={}_0+M\left(b\sqrt{E_s/{}_s{b}^2}\right)+{}_1{\left\{1-{\left[\frac{k_{b}T}{F}\ln \left(\frac{{\stackrel{.}{}}_p}{{\stackrel{}{}}_0}\right)\right]}^{\frac{1}{p}}\right\}}^{\frac{1}{q}},$

[0151] 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.

[0152] (6) Using a predicted stress as an input and is substituted into a main cutting force prediction formula, to analyze and judge the cutting force prediction results.

[0153] Wherein, the main cutting force F.sub.c prediction formula is as follows:

[00049]$F_c=\frac{a_w{a}_c\cos \left(-{}_e\right)}{\cos \left({}_e+-{}_e\right)\sin {}_e},$ [0154] wherein, .sub.w is a cutting width and .sub.c is a cutting depth.

[0155] When the predicted cutting force, cutting temperature, and/or chip morphology obtained from a certain prediction calculation meet the predefined steady-state requirements of the cutting state, the pre-used data of the cutting speed and cutting depth input in the certain prediction calculation are used as application values of steady-state control of a cutting tool in the actual cutting process.

[0156] 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. Hardness variations of the cut shear zone may be obtained based on a mapping relationship between micro-hardness and stored energy, and a microscopic stored energy distribution of the cut shear zone may be further obtained according to the stored energy distribution.

[0157] (7) In the actual cutting process, the application values of the cutting speed and cutting depth can be input into a computer numerical control (CNC) machine tool, or operating corresponding control levers on a regular cutting lathe by an operator according to the application values of the cutting speed and cutting depth, the cutting tool can be controlled and adjusted to perform a steady-state machining on the workpiece material to be machined according to the input application values during the actual cutting process, and ultimately obtain the expected material properties of the machined workpiece material.

[0158] It can be understood that the predefined steady-state requirements for cutting states described in the present example may vary depending on the workpiece material to be machined, as well as the type and model of the cutting tool. Those skills in the art can set and adjust them based on experience, actual needs, and/or special requirements. The present invention is not particularly limited.

[0159] It can be understood that the example can adopt the example described in Example 1, but is not limited to the example described. Specifically, the workpiece material to be machined is the 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-00003 Cutting Depth Group speed Feed of cut No.: (m/min) (mm/r) (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

[0160] The cutting force, the cutting temperature, and the forming law of the cutting chip, etc., are predicted. Based on the prediction results, and taking the adiabatic shear hysteresis effect into consideration, the critical velocity for generating serrated chips is approximately 30 m/min.

[0161] 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.