System and Method for Acquiring Inherent Energy Efficiency Factor Function of Computerized Numerically Controlled Machine Tool

Abstract

The disclosure discloses a system for acquiring an inherent energy efficiency factor function of a CNC machine tool, the system comprises an equipment information management module, a test parameter setting module, a NC code generation module, a field test management module, a data analysis module, a validity verification module and power sensors. The NC code generation module is used to generate the NC code to control the operation of the CNC machine tool to be tested according to the test parameters and information of the CNC machine tool to be tested. The field test management module is used to analyze the power data information of power sensors and generate the input power sets of CNC machine tool in each test operation process; The data analysis module is used to generate the inherent energy efficiency factor function according to the input power set.

Claims

1. A system for acquiring an inherent energy efficiency factor function of a computerized numerically controlled (CNC) machine tool, wherein the CNC machine tool comprises the following energy consumption subsystems: a CNC system, a spindle system, a feed shaft system, an automatic tool changing system and an auxiliary system, the system comprises: an equipment information management module, a test parameter setting module, a numerically controlled (NC) code generation module, a field test management module, a data analysis module and power sensors configured to install on the CNC machine tool to be tested and obtain the power data of the CNC machine tool to be tested; the equipment information management module is configured to input the basic information, spindle system information, feed shaft system information, automatic tool automatic tool changing system information and auxiliary system information of the CNC machine tool to be tested; the test parameter setting module is configured to set test parameters comprising a state-running duration t.sub.o, a state-switching-mark duration t.sub.d, a state-switching-delay duration t.sub.dc, number of training samples N.sub.tr and an initial distance of a feed shaft d.sub.I-o; the NC code generation module is configured to generate NC codes to control a operation according to the test parameters and the basic information, the spindle system information, the feed shaft system information, the tool automatic tool changing system information and the auxiliary system information of the NC machine tool to be tested, the NC codes comprises NC codes configured to control the spindle system running in a test process, NC codes configured to control the feed shaft system running in the test process, the NC codes configured to control the automatic tool changing system running in the test process and NC codes configured to control the auxiliary system running in the test process; the field test management module is configured to set parameters of the power sensors, read power data information of the power sensors, parse power data information of the power sensors and generate power set; the data analysis module is configured to generate an inherent energy efficiency factor function of each energy consumption subsystem according to the power set of the CNC machine tool to be tested during the test operation, and generate an inherent energy efficiency factor function of each operation stage of the CNC machine tool according to the inherent energy efficiency factor function of each energy consumption subsystem.

2. The system as claimed in claim 1, wherein the basic information comprises a machine type, a machine model, and a CNC system type, the spindle system information comprises a maximum spindle speed n.sub.S-max, a minimum spindle speed n.sub.S-min, a rated spindle speed n.sub.S-r, a speed interval ΔS-n-u above a rated speed and a speed interval ΔS-n-l below the rated speed, the feed axis system information comprises a fast feed speed fv.sub.I-max in a direction of I axis, a axial stroke d.sub.I, and a maximum cutting feed speed cv.sub.I-max, the subscript I∈[X, Y, Z], X, Y, Z represent a X-axis, a Y-axis, and a Z-axis of the machine tool respectively, the automatic tool automatic tool changing system information comprises the number of tool positions in a tool library N.sub.t, the auxiliary system information comprises a total number of auxiliary systems controlled by the CNC system N.sub.Au, a name of each auxiliary system, and a start-stop control code of each auxiliary system.

3. The system as claimed in claim 2, wherein the NC codes configured to control the operation of the spindle system in the test process comprises the spindle starting speed setting code, and a total number of spindle speed test is 2 N.sub.tr, and the i-th spindle speed ns.sub.S-tr[i] is set according to a following formula: ${\mathrm{ns}}_{S-\mathrm{tr}}\left[i\right]=\{\begin{array}{c}{n}_{S-r}-{\Delta }_{S-n-i}.Math.\frac{n_{S-\min }-n_{s-r}}{N_{tr}{\Delta }_{S-n-l}}.Math.\left(i-N_{tr}\right),0\le i<N_{tr}\\ n_{S-r}+{\Delta }_{S-n-u}.Math.\frac{n_{S-\max }-n_{S-r}}{N_{tr}{\Delta }_{S-n-u}}.Math.\left(i-N_{tr}\right),N_{tr}\le i\le 2N_{tr}\end{array},$ and [ ] represents the integer operation; the NC codes configured to control the running of the feed shaft system in the test process comprises a rapid feed&return distance setting code, a total number of rapid feed&return distance tests is N.sub.tr, and the i-th rapid feed&return distance d.sub.I-tr[i] is set according to a following formula: $d_{I-\mathrm{tr}}\left[i\right]=d_{I-o}+\frac{d_I-2d_{I-o}}{N_{tr}}i,$ and 1≤i≤N.sub.tr, d.sub.I is a axial travel of the feed shaft in a direction of I axis, d.sub.I-o represents an initial distance of the feed shaft in the direction of I axis; the NC codes configured to control the feed shaft system in the test process also comprises a cutting feed speed setting code, the total number of cutting feed rate tests is N.sub.tr, and the i-th cutting feed rate cv.sub.I-tr[i] is set according to a following formula: $c{v}_{I-\mathrm{tr}}\left[i\right]=\frac{{\mathrm{cv}}_{I=\max }}{N_{tr}}i,$ and 1≤i≤N.sub.tr, cv.sub.I-max represents a maximum cutting feed rate in the I-axis direction.

4. The system as claimed in claim 2, wherein the power set generated by the field test management module according to the power data information collected by the power sensors in the test process of CNC machine tool comprises: an input power set of the CNC machine tool in a power supply opening process D.sub.SP; an input power set of the CNC machine tool in a CNC system operation process D.sub.SC; the input power set D.sub.Au-i of the CNC machine tool during an operation of an i-th auxiliary system, 1≤i≤N.sub.Au; the input power set D.sub.PS-i of the CNC machine tool, in a start-up process of an i-th spindle speed, 1≤i≤2N.sub.tr; an input power set D.sub.US-i of the CNC machine tool in a process of an i-th spindle speed empty operation; an input power set D.sub.I-PF-i of the CNC machine tool in a process of an i-th rapid feed&return distance in the feed shaft I, 1≤i≤N.sub.tr; an input power set D.sub.I-UF-i of the CNC machine tool in a feed axis I cutting a feed process at an i-th cutting speed, 1≤i≤N.sub.tr; an input power set D.sub.PT-n.sub.t of the CNC machine tool in a process of changing an n.sub.t-th tool position by the automatic tool changing system.

5. The system as claimed in claim 1, wherein the data analysis module is configured to establish an inherent energy efficiency factor discrete function of each energy consumption subsystem of the CNC machine tool by a discrete modeling method, comprising following steps: step 2.1: calculating a power mean P.sub.D, a set time t.sub.D and a set energy consumption E.sub.D for each power set, according to following general formula: $\stackrel{\_}{P_D}=\frac{1}{N_D}\underset{k=1}{\overset{N_D}{.Math.}}{P}_k,P_k\in D;t_D=\frac{N_D}{fs};E_D=\frac{1}{fs}\underset{k=1}{\overset{N_D}{.Math.}}{P}_k,P_k\in D,$ wherein D represents the power set, P.sub.k represents a k-th element in the power set D, N.sub.D represents the number of elements in the power set D, and fs represents a sampling frequency of the power sensors; step 2.2: setting up the inherent energy efficiency factor discrete function of each energy consumption subsystem of the CNC machine tool: wherein a power consumption PSP of the CNC machine tool in a power supply opening process:
P.sub.SP=P.sub.in-SP, and P.sub.in-SP represents a power average of the power set D.sub.SP of the CNC machine tool in the power supply opening process, a power consumption of a control system of the CNC machine tool P.sub.Cs:
P.sub.Cs=P.sub.in-SC−P.sub.SP, and P.sub.in-SC represents a power mean of the power set D.sub.SC of the CNC machine tool in the operation process of the CNC system, a power consumption of the auxiliary system P.sub.Au-i:
P.sub.Au-i=P.sub.in-Au-i−P.sub.in-SC, and P.sub.in-Au-t represents a power average of the power set D.sub.Au-i of the CNC machine tool during the operation of the i-th auxiliary system, 1≤i≤N.sub.Au, a discrete function of startup time and energy consumption of the spindle system is C.sub.S-PS: $C_{S-PS}=\left[\begin{array}{ccc}n& t_{S-\mathrm{PS}}& E_{S-PS}\end{array}\right]=\left[\begin{array}{ccc}n{s}_{S-tr}\left[1\right]& t_{S-PS}\left[1\right]& E_{S-PS}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ n{s}_{S-tr}\left[i\right]& t_{S-PS}\left[i\right]& E_{S-PS}\left[i\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{ns}}_{S-\mathrm{tr}}\left[2N_{\mathrm{tr}}\right]& t_{S-\mathrm{PS}}\left[2N_{\mathrm{tr}}\right]& E_{S-\mathrm{PS}}\left[2N_{\mathrm{tr}}\right]\end{array}\right]$ and n represents a spindle speed and is an independent variable; both t.sub.S-PS and E.sub.S-PS are dependent variables, which are a starting time of the spindle system and a energy consumption of the spindle system. the t.sub.S-PS[i] is a running time of the power set D.sub.PS-i in the starting process of the CNC machine tool in the start-up process of the i-th spindle speed, 1≤i≤2 N.sub.tr, a calculation formula of E.sub.S-PS[i] is as follows:
E.sub.S-PS[i]=E.sub.in-PS[i]−P.sub.in-SCt.sub.S-PS[i], and E.sub.S-PS[i] is the set energy consumption of the power set D.sub.PS-i of the CNC machine tool in a start-up process of an i-th spindle speed, 1≤i≤2 N.sub.tr, an idling power discrete function of the spindle system in a process of air operation C.sub.S-US: $C_{S-\mathrm{US}}=\left[\begin{array}{cc}n& P_{S-\mathrm{US}}\end{array}\right]=\left[\begin{array}{cc}n{s}_{S-tr}\left[1\right]& P_{S-\mathrm{US}}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}\\ n{s}_{S-tr}\left[i\right]& P_{S-\mathrm{US}}\left[i\right]\\ \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{ns}}_{S-\mathrm{tr}}\left[2N_{\mathrm{tr}}\right]& P_{S-\mathrm{US}}\left[2N_{\mathrm{tr}}\right]\end{array}\right],$ and n represents the spindle speed and is the independent variable, P.sub.S-US represents a power function of a spindle system in an air operation process as a dependent variable, a calculation power formula of P.sub.S-US[i] is as follows:
P.sub.S-US[i]=P.sub.in-US[i]−P.sub.in-SC, and P.sub.in-US[i] is a mean power of the power set D.sub.US-i of the CNC machine tool in the process of the i-th spindle speed empty operation, 1≤i≤2 N.sub.tr, a discrete function C.sub.I-PF of a rapid feed&return time and an energy consumption of the rapid feed&return of the feed shaft system of the CNC machine tool is: $C_{I-PF}=\left[\begin{array}{ccc}{d}_{I-\mathrm{tr}}& t_{I-PF}& E_{I-PF}\end{array}\right]=\left[\begin{array}{ccc}{d}_{I­tr}\left[1\right]& t_{I-PF}\left[1\right]& E_{I­PF}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ d_{I­tr}\left[i\right]& t_{I-PF}\left[i\right]& E_{I-\mathrm{PF}}\left[i\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ d_{I­tr}\left[N_{\mathrm{tr}}\right]& t_{I-PF}\left[N_{tr}\right]& E_{I-\mathrm{PF}}\left[N_{\mathrm{tr}}\right]\end{array}\right],$ and the subscripts I∈[X, Y, Z], X, Y, Z represent the X axis, the Y axis and the Z axis of the machine tool respectively, d.sub.I-tr represents a rapid feed&return distance and is an independent variable, t.sub.I-PF and E.sub.I-PF are dependent variables, t.sub.I-PF means a rapid feed&return time, E.sub.I-PF means a rapid feed&return energy consumption, t.sub.I-PF[i] is a running time of a feed axis I in the power set D.sub.I-PF-i of the CNC machine tool in a process of an i-th rapid feed&return distance in the feed shaft I, 1≤i≤N.sub.tr, a calculation formula of E.sub.I-FF[i] is as follows:
E.sub.I-FF[i]=E.sub.in-I-PF[i]−P.sub.in-SCt.sub.I-PF[i], and E.sub.in-I-PF[i] is a set energy consumption of the power set D.sub.I-PF-i of the CNC machine tool in a process of an i-th rapid feed&return distance in the feed shaft I, 1≤i≤N.sub.tr, a discrete function C.sub.I-UF of feed power in a cutting feed operation process of the feed shaft system is: $C_{I-\mathrm{UF}}=\left[\begin{array}{cc}{\mathrm{cv}}_{I-\mathrm{tr}}& P_{I-\mathrm{UF}}\end{array}\right]=\left[\begin{array}{cc}{\mathrm{cv}}_{I-\mathrm{tr}}\left[1\right]& P_{I-\mathrm{UF}}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{cv}}_{I-\mathrm{tr}}\left[i\right]& P_{I-\mathrm{UF}}\left[i\right]\\ \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{cv}}_{I-\mathrm{tr}}\left[N_{\mathrm{tr}}\right]& P_{I-\mathrm{UF}}\left[N_{\mathrm{tr}}\right]\end{array}\right],$ and cv.sub.I-tr represents a cutting feed rate of I axis, which is an independent variable, P.sub.I-UF represents a feed power of I-axis cutting and is a dependent variable, a calculation formula of P.sub.I-UF[i] is as follows:
P.sub.I-UF[i]=P.sub.in-I-UF[i]−P.sub.in-SC, and P.sub.in-I-UF[i] is a mean power of the power set D.sub.I-UF-i of the CNC machine tool in the feed axis I cutting the feed process at the i-th cutting speed, 1≤i≤N.sub.tr, a discrete function C.sub.PT of a tool changing time and energy consumption of the automatic tool changing system during a tool changing operation is: $C_{\mathrm{PT}}=\left[\begin{array}{ccc}{n}_t& t_{\mathrm{PT}}& E_{\mathrm{PT}}\end{array}\right]=\left[\begin{array}{ccc}1& t_{\mathrm{PT}}\left[1\right]& E_{\mathrm{PT}}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ n_t& t_{\mathrm{PT}}\left[n_t\right]& E_{\mathrm{PT}}\left[n_t\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ N_t& t_{\mathrm{PT}}\left[N_t\right]& E_{\mathrm{PT}}\left[N_t\right]\end{array}\right],$ and n.sub.i represents a tool change position and is an independent variable, t.sub.PT and E.sub.PT are dependent variables, t.sub.PT means a tool change time, E.sub.PT means a tool change energy consumption of the tool changing system, t.sub.PT[n.sub.t] is a running time for the power set D.sub.PT-nt of the CNC machine tool in a process of changing an n.sub.t-th tool position by the automatic tool changing system, a calculation formula of E.sub.PT[n.sub.t] is as follows:
E.sub.PT[n.sub.t]=E.sub.in-PT[n.sub.t]−P.sub.in-SCt.sub.PT[n.sub.t], and E.sub.in-PT [n.sub.t] is a set energy consumption of the power set D.sub.PT-n.sub.t of the CNC machine tool in a process of changing an n.sub.t-th tool position by the automatic tool changing system, 1≤n.sub.t≤N.sub.t.

6. The system as claimed in claim 5, wherein the data analysis module establishes an inherent energy efficiency factor fitting function of each energy consumption subsystem of CNC machine tool by adaptive fitting modeling method based on discrete modeling, comprising following steps: step 2.3: regarding the inherent energy efficiency factor discrete function of each energy consumption subsystem of the CNC machine tool established by the discrete modeling method as a data set for fitting, and dividing the data set into a training set and a test set by a cross-validation method in machine learning; step 2.4: fitting a first fitting function and a second fitting function by a least square method according to the training set; step 2.5: by the test set, calculating errors of the first fitting function and the second fitting function according to an error function; step 2.6: If the error of the first fitting function is less than the error of the second fitting function, selecting the first fitting function as the inherent energy efficiency factor function, otherwise, selecting the second fitting function as the inherent energy efficiency factor function; therefore, establishing following fitting functions of each energy consumption subsystem: the spindle system of CNC machine tool: a fitting function of a starting time length t.sub.S-PS(n), a fitting function of a starting energy consumption E.sub.S-PS(n) and a fitting function of an air operation power P.sub.S-US(n) which are with the spindle speed n as the independent variable respectively; the feed shaft system of CNC machine tool: a rapid feed&return time fitting function t.sub.I-PF(d.sub.I) which is with a rapid feed&return distance d.sub.I as an independent variable, a rapid feed&return energy consumption fitting function E.sub.I-PF(d.sub.I) which is with the rapid feed&return distance d.sub.I as the independent variable, a cutting feed power fitting function P.sub.I-UF(cv.sub.I) with a cutting feed rate cv.sub.I as a independent variable; the automatic tool changing system: a tool change time fitting function t.sub.PT(n.sub.t) and a tool change energy consumption fitting function E.sub.PT(n.sub.t) which are with the tool change position n.sub.t as an independent variable.

7. The system as claimed in claim 6, wherein the inherent energy efficiency factor functions of the CNC machine tool comprise a power function P.sub.in-S(s.sub.Cs) of the CNC machine tool in a standby stage, a power function P.sub.in-PA(s.sub.Au-i) of the CNC machine tool in an opening stage of the auxiliary system, an energy consumption function E.sub.in-PS(n,s.sub.Au-i) of the CNC machine tool in a starting stage of the spindle system, an energy consumption function E.sub.in-PF(d.sub.I,s.sub.Au-i) of the CNC machine tool in a rapid feed&return stage, an energy consumption function E.sub.in-PT(n.sub.t,s.sub.Au-i) of the CNC machine tool in an automatic tool change stage and a power function P.sub.in-U(n,cv.sub.I,s.sub.Au-i) of the CNC machine tool in an idling stage, specific expressions are as follows: the power function P.sub.in-S(s.sub.Cs) of the CNC machine tool in the standby stage:
P.sub.in-S(s.sub.Cs)=P.sub.SP+s.sub.CsP.sub.Cs, and P.sub.SP represents a power consumption when a CNC machine tool power is turned on, and P.sub.Cs represents a new power consumption when the CNC machine tool power is turned on, s.sub.Cs means a state of the CNC system, s.sub.Cs=0 means the state of the CNC system is closed, s.sub.Cs=1 means the state of the CNC system is open; the power function of the CNC machine tool in an opening stage of the auxiliary system P.sub.in-PA(s.sub.Au-i):
P.sub.in-PA(s.sub.Au-i)P.sub.S+Σs.sub.Au-iP.sub.Au-i, P.sub.S represents a standby power of CNC machine tool, P.sub.S=P.sub.in-S(s.sub.Cs=1); P.sub.Au-i means a new machine tool power consumption when a first auxiliary system is opened, 1≤i≤N.sub.Au; s.sub.Au-i represents a state of the i-th auxiliary system, 1≤i≤N.sub.Au, s.sub.Au-i=0 represents a closed state of the i-th auxiliary system, and s.sub.Au-i=1 represents an open state of the i-th auxiliary system; the energy consumption function E.sub.in-PS(n,s.sub.Au-i) of CNC machine tool in the starting stage of the spindle system:
E.sub.in-PS(n,s.sub.Au-t)=(P.sub.S+Σs.sub.Au-tP.sub.Au-t)t.sub.S-PS(n)+E.sub.S-PS(n), and n denotes the spindle speed, t.sub.S-PS(n) and E.sub.S-PS(n) are a starting time function and an energy consumption function of the spindle system, respectively; the energy consumption function E.sub.in-PF(d.sub.I,s.sub.Au-i) of CNC machine tool in the rapid feed&return stage:
E.sub.in-PF(d.sub.I,s.sub.Au-i)=(P.sub.S+Σs.sub.Au-iP.sub.Au-i)Σt.sub.I-PF(d.sub.I)+ΣE.sub.I-PF(d.sub.I), and t.sub.PT(n.sub.t) and E.sub.PT(n.sub.t) represent the rapid feed&return time fitting function, the rapid feed&return energy consumption fitting function of I axis respectively; the energy consumption function E.sub.in-PT(n.sub.t,s.sub.Au-i) of CNC machine tool in the automatic tool change stage:
E.sub.in-PT(n.sub.t,s.sub.Au-i)=(P.sub.S+Σs.sub.Au-iP.sub.Au-i)t.sub.PT(n.sub.t)+E.sub.PT(n.sub.t), and n.sub.t represents the number of cutter positions that need to be rotated during automatic tool change; t.sub.PT(n.sub.t) and E.sub.PT(n.sub.t) represents the tool change time fitting function and the tool change energy consumption fitting function of the automatic tool changing system respectively; the power function P.sub.in-U(n,cv.sub.I,s.sub.Au-i) in the idling stage of the CNC machine tool:
P.sub.in-U(n,cv.sub.I,s.sub.Au-t)=P.sub.S+P.sub.S-US(n)+Σs.sub.Au-iP.sub.Au-i+ΣP.sub.I-UF(cv.sub.I), Where P.sub.S-US(n) represents the a fitting function of an air operation power of the spindle system; cv.sub.I represents the feed rate and P.sub.I-UF(cv.sub.I) represents the a cutting feed power fitting function of the I axis.

8. The system as claimed in claim 1, wherein the system comprises: a validity verification module, configured to verify a effectiveness of the inherent energy efficiency factor function, an judge whether the inherent energy efficiency factor function is effective according to an error between a verification value and a predicted value, wherein the verification value is obtained by a verification experiment, and the predicted value is calculated by the data analysis module according to operation conditions of a verification experiment and the inherent energy efficiency factor function of the CNC machine tool in each operation stage, to carry out the verification experiment, the NC code generation module and the field test management module are improved as follows: the NC code generation module is also configured to generate a NC code for a verification experiment test according to the test parameters and the basic information, the spindle system information, the feed shaft system information, the tool automatic tool changing system information, and the auxiliary system information, the NC code for the verification experiment test comprises a NC code for controlling the operation of the auxiliary system in the verification experiment, a NC code for controlling a starting and an empty operation of the spindle system according to a specified speed n.sub.S-V in the verification experiment, a NC code for controlling the feed shaft system according to a specified rapid feed&return distance d.sub.I-V and a specified cutting feed speed cv.sub.I-V in the verification experiment, and a NC code for controlling the automatic tool changing system according to a specified tool position n.sub.t-V in the verification experiment; the power set generated by the field test management module according to the power data information collected by the power sensors in the verification experiment of the CNC machine tool comprises: an input power set D.sub.SC-V of the CNC machine tool in an operation process of the CNC system, an input power set D.sub.Au-V of the CNC machine tool in a whole operation process of the auxiliary system, an input power set D.sub.PS-V of the CNC machine tool in the spindle system according to the specified speed n.sub.S-V, an input power set D.sub.US-V of the CNC machine tool in an empty operation process of the spindle system according to the specified speed n.sub.S-V, an input power set D.sub.I-PF-V of the CNC machine tool in the feed shaft I (I∈[X, Y, Z]) in the feed shaft d.sub.I-V.

9. A method for obtaining an inherent energy efficiency factors of a CNC machine tool, comprising, establishing an inherent energy efficiency factor function of each energy consumption subsystem and an inherent energy efficiency factor function of each operation stage of a CNC machine tool by a system for acquiring an inherent energy efficiency factor function of a CNC machine tool according to claim 1, and calculating inherent energy efficiency factors of the CNC machine tool in different operation processes according to the inherent energy efficiency factor function of the CNC machine tool, specifically comprising the following steps: step A1: equipment information management module inputting basic information, active system information, and auxiliary system information of the CNC machine tool to be tested; step A2: setting test parameters in a test parameter setting module, wherein the test parameters comprise: a state-running duration t.sub.o, a state switching flag time t.sub.d, a state-switching-delay duration t.sub.dc, the number of test samples N.sub.tr and an initial distance of a feed shaft d.sub.I-o; step A3: generating NC codes of a test experiment in a NC code generation module according to the test parameters in step A2, the basic information, the spindle system information, the feed shaft system information, the tool automatic tool changing system information, and the auxiliary system information of the CNC machine tool in step A1; step A4: setting power sensors parameters in a field test module, and inputting the NC codes of the test experiment into a CNC system of the CNC machine tool to be tested; step A5: the CNC machine tool running the NC codes of the test experiment and controlling the CNC machine tool to be tested to run according to the NC codes of the test experiment, at the same time, collecting power data information of each energy consumption subsystem in the test experiment of the CNC machine tool to be tested by the power sensors installed on the CNC machine tool to be tested, generating a power set of each energy consumption subsystem in the test experiment by the field test management module according to the power data information collected by the power sensors; step A6: a data analysis module establishing the inherent energy efficiency factor function of each energy consumption subsystem according to the power set in step A5; step A7: the data analysis module establishing the inherent energy efficiency factor function of CNC machine tool in each operation stage according to the inherent energy efficiency factor function of each energy consumption subsystem in step A6; step A8: according to the inherent energy efficiency factor function of CNC machine tool in each operation stage and operating conditions of CNC machine tool, calculating inherent energy efficiency factors in each operation stage of CNC machine tool.

10. A method for obtaining inherent energy efficiency factors of a CNC machine tool, comprising, establishing an inherent energy efficiency factor function of each energy consumption subsystem and an inherent energy efficiency factor function of each operation stage of a CNC machine tool by a system for acquiring an inherent energy efficiency factor function of a CNC machine tool according to claim 8, and calculating inherent energy efficiency factors of the CNC machine tool in different operation processes according to the inherent energy efficiency factor function of the CNC machine tool, specifically comprising the following steps: step B1: equipment information management module inputting basic information, active system information, and auxiliary system information of the CNC machine tool to be tested; step B2: setting test parameters in a test parameter setting module, wherein the test parameters comprise: a state-running duration t.sub.o, a state switching flag time t.sub.d, a state-switching-delay duration t.sub.dc, the number of test samples N.sub.tr, and an initial distance of a feed shaft d.sub.I-o; step B3: generating NC codes of a test experiment and NC codes of a verification experiment in a NC code generation module according to the test parameters in step B2, the basic information, the spindle system information, the feed shaft system information, the tool automatic tool changing system information, and the auxiliary system information of the CNC machine tool in step B1; step B4: setting power sensors parameters in a field test module, and inputting the NC codes of the test experiment and the NC codes of the verification experiment into a CNC system of the CNC machine tool to be tested; step B5: the CNC machine tool running the NC codes of the test experiment and controlling the CNC machine tool to be tested to run according to the NC codes of the test experiment, at the same time, collecting power data information of each energy consumption subsystem in the test experiment of the CNC machine tool to be tested by the power sensors installed on the CNC machine tool to be tested, generating a power set of each energy consumption subsystem in the test experiment by the field test management module according to the power data information collected by the power sensors; step B6: the CNC machine tool running the NC codes of the verification experiment and controlling the CNC machine tool to be tested to run according to the NC codes of the verification experiment, at the same time, collecting power data information of each energy consumption subsystem in the test experiment of the CNC machine tool to be tested by the power sensors installed on the CNC machine tool to be tested, generating a power set of each energy consumption subsystem in the verification experiment by the field test management module according to the power data information collected by the power sensors; step B7: a data analysis module establishing the inherent energy efficiency factor function of each energy consumption subsystem according to the power set of the test experiment in step B5; step B8: the data analysis module establishing the inherent energy efficiency factor function of CNC machine tool in each operation stage according to the inherent energy efficiency factor function of each energy consumption subsystem in step A6; step B9: a validity verification module carrying out a validity verification of the inherent energy efficiency factor function of the CNC machine tool generated by the data analysis module, if the inherent energy efficiency factor function passes the validity verification, Entering step B10. if the inherent energy efficiency factor function fails to pass the validity verification, repeating steps B1 to B9, if still fail to pass the validity verification, contact technical personnel to eliminate the problem; step B10: according to the inherent energy efficiency factor function of CNC machine tool in each operation stage and operating conditions of CNC machine tool, calculating inherent energy efficiency factors in each operation stage of CNC machine tool.

11. The method as claimed in claim 9, wherein the basic information comprises a machine type, a machine model, and a CNC system type, the spindle system information comprises a maximum spindle speed n.sub.S-max, a minimum spindle speed n.sub.S-min, a rated spindle speed n.sub.S-r, a speed interval ΔS-n-u above a rated speed and a speed interval ΔS-n-l below the rated speed, the feed axis system information comprises a fast feed speed fv.sub.I-max in a direction of I axis, a axial stroke d.sub.I, and a maximum cutting feed speed cv.sub.I-max, the subscript I∈[X, Y, Z], X, Y, Z represent a X-axis, a Y-axis, and a Z-axis of the machine tool respectively, the automatic tool automatic tool changing system information comprises the number of tool positions in a tool library N.sub.t, the auxiliary system information comprises a total number of auxiliary systems controlled by the CNC system N.sub.Au, a name of each auxiliary system, and a start-stop control code of each auxiliary system.

12. The method as claimed in claim 11, wherein the NC codes configured to control the operation of the spindle system in the test process comprises the spindle starting speed setting code, and a total number of spindle speed test is 2 N.sub.tr, and the i-th spindle speed ns.sub.S-tr[i] is set according to a following formula: $n{s}_{S-tr}\left[i\right]=\{\begin{array}{c}{n}_{S-r}-{\Delta }_{S-n-l}.Math.\frac{n_{S-\min }-n_{S-r}}{N_{tr}{\Delta }_{S-n-l}}.Math.\left(i-N_{tr}\right),0\le i<N_{tr}\\ n_{S-r}+{\Delta }_{S-n-u}.Math.\frac{n_{S-\max }-n_{S-r}}{N_{tr}{\Delta }_{S-n-u}}.Math.\left(i-N_{tr}\right),N_{tr}\le i\le 2N_{tr}\end{array},$ and [ ] represents the integer operation; the NC codes configured to control the running of the feed shaft system in the test process comprises a rapid feed&return distance setting code, a total number of rapid feed&return distance tests is N.sub.tr, and the i-th rapid feed&return distance d.sub.I-tr[i] is set according to a following formula: $d_{I-\mathrm{tr}}\left[i\right]=d_{I-o}+\frac{d_I-2d_{I-o}}{N_{tr}}i,$ and 1≤i≤N.sub.tr, d.sub.I is a axial travel of the feed shaft in a direction of I axis, d.sub.I-o represents an initial distance of the feed shaft in the direction of I axis; the NC codes configured to control the feed shaft system in the test process also comprises a cutting feed speed setting code, the total number of cutting feed rate tests is N.sub.tr, and the i-th cutting feed rate Cv.sub.I-tr[i] is set according to a following formula: ${\mathrm{cv}}_{I-\mathrm{tr}}\left[i\right]=\frac{c{v}_{I-\max }}{N_{tr}}i,$ and 1≤i≤N.sub.tr, Cv.sub.I-max represents a maximum cutting feed rate in the I-axis direction.

13. The method as claimed in claim 11, wherein the power set generated by the field test management module according to the power data information collected by the power sensors in the test process of CNC machine tool comprises: an input power set of the CNC machine tool in a power supply opening process D.sub.SP; an input power set of the CNC machine tool in a CNC system operation process D.sub.SC; the input power set D.sub.Au-i of the CNC machine tool during an operation of an i-th auxiliary system, 1≤i≤N.sub.Au; the input power set D.sub.PS-i of the CNC machine tool, in a start-up process of an i-th spindle speed, 1≤i≤2 N.sub.tr; an input power set D.sub.US-i of the CNC machine tool in a process of an i-th spindle speed empty operation; an input power set D.sub.I-PF-i of the CNC machine tool in a process of an i-th rapid feed&return distance in the feed shaft I, 1≤i≤N.sub.tr; an input power set D.sub.I-UF-i of the CNC machine tool in a feed axis I cutting a feed process at an i-th cutting speed, 1≤i≤N.sub.tr; an input power set D.sub.PT-n.sub.t of the CNC machine tool in a process of changing an n.sub.t-th tool position by the automatic tool changing system.

14. The method as claimed in claim 9, wherein the data analysis module is configured to establish an inherent energy efficiency factor discrete function of each energy consumption subsystem of the CNC machine tool by a discrete modeling method, comprising following steps: step 2.1: calculating a power mean P.sub.D, a set time t.sub.D and a set energy consumption E.sub.D for each power set, according to following general formula: $\stackrel{\_}{P_D}=\frac{1}{N_D}\underset{k=1}{\overset{N_D}{.Math.}}{P}_k,P_k\in D;t_D=\frac{N_D}{fs};E_D=\frac{1}{fs}\underset{k=1}{\overset{N_D}{.Math.}}{P}_k,P_k\in D,$ wherein D represents the power set, P.sub.k represents a k-th element in the power set D, N.sub.D represents the number of elements in the power set D, and fs represents a sampling frequency of the power sensors; step 2.2: setting up the inherent energy efficiency factor discrete function of each energy consumption subsystem of the CNC machine tool: wherein a power consumption PSP of the CNC machine tool in a power supply opening process:
P.sub.SP=P.sub.in-SP, and P.sub.in-SP represents a power average of the power set D.sub.SP of the CNC machine tool in the power supply opening process, a power consumption of a control system of the CNC machine tool P.sub.Cs:
P.sub.Cs=P.sub.in-SC−P.sub.SP, and P.sub.in-SC represents a power mean of the power set D.sub.SC of the CNC machine tool in the operation process of the CNC system, a power consumption of the auxiliary system P.sub.Au-i:
P.sub.Au-i=P.sub.in-Au-i−P.sub.in-SC, and P.sub.in-Au-i represents a power average of the power set D.sub.Au-I of the CNC machine tool during the operation of the i-th auxiliary system, 1≤i≤N.sub.Au, a discrete function of startup time and energy consumption of the spindle system is C.sub.S-PS: $C_{S-\mathrm{PS}}=\left[\begin{array}{ccc}n& t_{S-\mathrm{PS}}& E_{S-\mathrm{PS}}\end{array}\right]=\left[\begin{array}{ccc}n{s}_{S-\mathrm{tr}}\left[1\right]& t_{S-\mathrm{PS}}\left[1\right]& E_{S-\mathrm{PS}}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ n{s}_{S-\mathrm{tr}}\left[i\right]& t_{S-\mathrm{PS}}\left[i\right]& E_{S-\mathrm{PS}}\left[i\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ n{s}_{S-\mathrm{tr}}\left[2N_{\mathrm{tr}}\right]& t_{S-\mathrm{PS}}\left[2N_{\mathrm{tr}}\right]& E_{S-\mathrm{PS}}\left[2N_{\mathrm{tr}}\right]\end{array}\right],$ and n represents a spindle speed and is an independent variable; both t.sub.S-PS and E.sub.S-PS are dependent variables, which are a starting time of the spindle system and a energy consumption of the spindle system. the t.sub.S-PS[i] is a running time of the power set D.sub.PS-i in the starting process of the CNC machine tool in the start-up process of the i-th spindle speed, 1≤i≤2 N.sub.tr, a calculation formula of E.sub.S-PS[i] is as follows:
E.sub.S-PS[i]=E.sub.in-PS[i]−P.sub.in-SCt.sub.S-PS[i], and E.sub.S-PS[i] is the set energy consumption of the power set D.sub.PS-i of the CNC machine tool in a start-up process of an i-th spindle speed, 1≤i≤2 N.sub.tr, an idling power discrete function of the spindle system in a process of air operation C.sub.S-US: $C_{S-\mathrm{US}}=\left[\begin{array}{cc}n& P_{S-\mathrm{US}}\end{array}\right]=\left[\begin{array}{cc}n{s}_{S-\mathrm{tr}}\left[1\right]& P_{S-\mathrm{US}}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}\\ n{s}_{S-\mathrm{tr}}\left[i\right]& P_{S-\mathrm{US}}\left[i\right]\\ \mathrm{.Math.}& \mathrm{.Math.}\\ n{s}_{S-\mathrm{tr}}\left[2N_{\mathrm{tr}}\right]& P_{S-\mathrm{US}}\left[2N_{\mathrm{tr}}\right]\end{array}\right],$ and n represents the spindle speed and is the independent variable, P.sub.S-US represents a power function of a spindle system in an air operation process as a dependent variable, a calculation power formula of P.sub.S-US[i] is as follows:
P.sub.S-US[i]=P.sub.in-US[i]−P.sub.in-SC, and P.sub.in-US[i] is a mean power of the power set D.sub.US-i of the CNC machine tool in the process of the i-th spindle speed empty operation, 1≤i≤2 N.sub.tr, a discrete function C.sub.I-PF of a rapid feed&return time and an energy consumption of the rapid feed&return of the feed shaft system of the CNC machine tool is: $C_{I-\mathrm{PF}}=\left[\begin{array}{ccc}{d}_{I-\mathrm{tr}}& t_{I-PF}& E_{I-PF}\end{array}\right]=\left[\begin{array}{ccc}{d}_{I-\mathrm{tr}}\left[1\right]& t_{I-PF}\left[1\right]& E_{I-PF}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ d_{I-\mathrm{tr}}\left[i\right]& t_{I-PF}\left[i\right]& E_{I-PF}\left[i\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ d_{I-\mathrm{tr}}\left[N_{\mathrm{tr}}\right]& t_{I-PF}\left[N_{\mathrm{tr}}\right]& E_{I-PF}\left[N_{\mathrm{tr}}\right]\end{array}\right]$ and the subscripts I∈[X, Y, Z], X, Y, Z represent the X axis, the Y axis and the Z axis of the machine tool respectively, d.sub.I-tr represents a rapid feed&return distance and is an independent variable, t.sub.I-PF and E.sub.I-PF are dependent variables, t.sub.I-PF means a rapid feed&return time, E.sub.I-PF means a rapid feed&return energy consumption, t.sub.I-PF is a running time of a feed axis I in the power set D.sub.I-PF-i of the CNC machine tool in a process of an i-th rapid feed&return distance in the feed shaft I, 1≤i≤N.sub.tr, a calculation formula of E.sub.I-FF[i] is as follows:
E.sub.I-FF[i]=E.sub.in-I-PF[i]−P.sub.in-SCt.sub.I-PF[i], and E.sub.in-I-PF[i] is a set energy consumption of the power set D.sub.I-PF-i of the CNC machine tool in a process of an i-th rapid feed&return distance in the feed shaft I, 1≤i≤N.sub.tr, a discrete function C.sub.I-UF of feed power in a cutting feed operation process of the feed shaft system is: $C_{I-\mathrm{UF}}=\left[\begin{array}{cc}{\mathrm{cv}}_{I-\mathrm{tr}}& P_{I-\mathrm{UF}}\end{array}\right]=\left[\begin{array}{cc}{\mathrm{cv}}_{I-\mathrm{tr}}\left[1\right]& P_{I-\mathrm{UF}}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{cv}}_{I-\mathrm{tr}}\left[i\right]& P_{I-\mathrm{UF}}\left[i\right]\\ \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{cv}}_{I-\mathrm{tr}}\left[N_{\mathrm{tr}}\right]& P_{I-\mathrm{UF}}\left[N_{\mathrm{tr}}\right]\end{array}\right],$ and cv.sub.I-tr represents a cutting feed rate of I axis, which is an independent variable, P.sub.I-UF represents a feed power of I-axis cutting and is a dependent variable, a calculation formula of P.sub.I-UF[i] is as follows:
P.sub.I-UF[i]=P.sub.in-I-UF[i]−P.sub.in-SC, and P.sub.in-I-UF[i] is a mean power of the power set D.sub.I-UF-i of the CNC machine tool in the feed axis I cutting the feed process at the i-th cutting speed, 1≤i≤N.sub.tr, a discrete function C.sub.PT of a tool changing time and energy consumption of the automatic tool changing system during a tool changing operation is: $C_{\mathrm{PT}}=\left[\begin{array}{ccc}{n}_t& t_{\mathrm{PT}}& E_{\mathrm{PT}}\end{array}\right]=\left[\begin{array}{ccc}1& t_{\mathrm{PT}}\left[1\right]& E_{\mathrm{PT}}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ n_t& t_{\mathrm{PT}}\left[n_t\right]& E_{\mathrm{PT}}\left[n_t\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ N_t& t_{\mathrm{PT}}\left[N_t\right]& E_{\mathrm{PT}}\left[N_t\right]\end{array}\right],$ and n.sub.t represents a tool change position and is an independent variable, t.sub.PT and E.sub.PT are dependent variables, t.sub.PT means a tool change time, E.sub.PT means a tool change energy consumption of the tool changing system, t.sub.PT[n.sub.t] is a running time for the power set D.sub.PT-nt of the CNC machine tool in a process of changing an n.sub.t-th tool position by the automatic tool changing system, a calculation formula of E.sub.PT[n.sub.t] is as follows:
E.sub.PT[n.sub.t]=E.sub.in-PT[n.sub.t]−P.sub.in-SCt.sub.PT[n.sub.t], and E.sub.in-PT [n,] is a set energy consumption of the power set D.sub.PT-n.sub.t of the CNC machine tool in a process of changing an n.sub.t-th tool position by the automatic tool changing system, 1≤n.sub.t≤N.sub.t.

15. The method as claimed in claim 14, wherein the data analysis module establishes an inherent energy efficiency factor fitting function of each energy consumption subsystem of CNC machine tool by adaptive fitting modeling method based on discrete modeling, comprising following steps: step 2.3: regarding the inherent energy efficiency factor discrete function of each energy consumption subsystem of the CNC machine tool established by the discrete modeling method as a data set for fitting, and dividing the data set into a training set and a test set by a cross-validation method in machine learning; step 2.4: fitting a first fitting function and a second fitting function by a least square method according to the training set; step 2.5: by the test set, calculating errors of the first fitting function and the second fitting function according to an error function; step 2.6: If the error of the first fitting function is less than the error of the second fitting function, selecting the first fitting function as the inherent energy efficiency factor function, otherwise, selecting the second fitting function as the inherent energy efficiency factor function; therefore, establishing following fitting functions of each energy consumption subsystem: the spindle system of CNC machine tool: a fitting function of a starting time length t.sub.S-PS(n), a fitting function of a starting energy consumption E.sub.S-PS(n) and a fitting function of an air operation power P.sub.S-US(n) which are with the spindle speed n as the independent variable respectively; the feed shaft system of CNC machine tool: a rapid feed&return time fitting function t.sub.I-PF(d.sub.I) which is with a rapid feed&return distance d.sub.I as an independent variable, a rapid feed&return energy consumption fitting function E.sub.I-PF(d.sub.I) which is with the rapid feed&return distance d.sub.I as the independent variable, a cutting feed power fitting function P.sub.I-UF(cv.sub.I) with a cutting feed rate cv.sub.I as a independent variable; the automatic tool changing system: a tool change time fitting function t.sub.PT(n.sub.t) and a tool change energy consumption fitting function E.sub.PT(n.sub.t) which are with the tool change position n.sub.t as an independent variable.

16. The method as claimed in claim 15, wherein the inherent energy efficiency factor functions of the CNC machine tool comprise a power function P.sub.in-S(s.sub.Cs) of the CNC machine tool in a standby stage, a power function P.sub.in-PA(s.sub.Au-i) of the CNC machine tool in an opening stage of the auxiliary system, an energy consumption function E.sub.in-PS(n,s.sub.Au-i) of the CNC machine tool in a starting stage of the spindle system, an energy consumption function E.sub.in-PF(d.sub.I,s.sub.Au-i) of the CNC machine tool in a rapid feed&return stage, an energy consumption function E.sub.in-PT(n.sub.t,s.sub.Au-i) of the CNC machine tool in an automatic tool change stage and a power function P.sub.in-U(n,cv.sub.I,s.sub.Au-i) of the CNC machine tool in an idling stage, specific expressions are as follows: the power function P.sub.in-S(s.sub.Cs) of the CNC machine tool in the standby stage:
P.sub.in-S(s.sub.Cs)=P.sub.SP+s.sub.CsP.sub.Cs, and P.sub.SP represents a power consumption when a CNC machine tool power is turned on, and P.sub.Cs represents a new power consumption when the CNC machine tool power is turned on, s.sub.Cs means a state of the CNC system, s.sub.Cs=0 means the state of the CNC system is closed, s.sub.Cs=1 means the state of the CNC system is open; the power function of the CNC machine tool in an opening stage of the auxiliary system P.sub.in-PA(s.sub.Au-i):
P.sub.in-PA(s.sub.Au-i)=P.sub.S+Σs.sub.Au-iP.sub.Ai-i, P.sub.S represents a standby power of CNC machine tool, P.sub.S=P.sub.in-S(s.sub.Cs=1); P.sub.Au-i means a new machine tool power consumption when a first auxiliary system is opened, 1≤i.Math.N.sub.Au; s.sub.Au-i represents a state of the i-th auxiliary system, 1≤i≤N.sub.Au, S.sub.Au-i=0 represents a closed state of the i-th auxiliary system, and s.sub.Au-i=1 represents an open state of the i-th auxiliary system; the energy consumption function E.sub.in-PS(n,s.sub.Au-i) of CNC machine tool in the starting stage of the spindle system:
E.sub.in-PS(n.sub.t,S.sub.Au-i)=(P.sub.S+Σs.sub.Au-iP.sub.Au-i)t.sub.S-PS(n)+E.sub.S-PS(n), and n denotes the spindle speed, t.sub.S-PS(n) and E.sub.S-PS(n) are a starting time function and an energy consumption function of the spindle system, respectively; the energy consumption function E.sub.in-PF(d.sub.I,s.sub.Au-i) of CNC machine tool in the rapid feed&return stage:
E.sub.in-PF(d.sub.I,s.sub.Au-i)=(P.sub.S+Σs.sub.Au-iP.sub.Au-i)Σt.sub.I-PF(d.sub.I)+Σ.sub.I-PF(d.sub.I), and t.sub.PT(n.sub.t) and E.sub.PT(n.sub.t) represent the rapid feed&return time fitting function, the rapid feed&return energy consumption fitting function of I axis respectively; the energy consumption function E.sub.in-PT(n.sub.t,s.sub.Au-i) of CNC machine tool in the automatic tool change stage:
E.sub.in-PT(n.sub.t,s.sub.Au-i)=(P.sub.S+Σs.sub.Au-iP.sub.Au-i)t.sub.PT(n.sub.t)+E.sub.PT(n.sub.t), and n.sub.t represents the number of cutter positions that need to be rotated during automatic tool change; t.sub.PT(n.sub.t) and E.sub.PT(n.sub.t) represents the tool change time fitting function and the tool change energy consumption fitting function of the automatic tool changing system respectively; the power function P.sub.in-U(n,cv.sub.I,s.sub.Au-i) in the idling stage of the CNC machine tool:
P.sub.in-U(n,cv.sub.I,s.sub.Au-i)=P.sub.S+P.sub.S-US(n)+Σs.sub.Au-iP.sub.Au-i+ΣP.sub.I-UF(cv.sub.I), Where P.sub.S-US(n) represents the a fitting function of an air operation power of the spindle system; cv.sub.I represents the feed rate and P.sub.I-UF(cv.sub.I) represents the a cutting feed power fitting function of the I axis.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0092] FIG. 1 shows the inherent energy efficiency factor function acquisition system framework of the CNC machine tool;

[0093] FIG. 2 shows the installation of the power sensors.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0094] FIG. 1 is the inherent energy efficiency factor function acquisition system framework of the CNC machine tool; The vertical CNC milling center XK714D is tested by using the inherent energy efficiency factor function acquisition system and acquisition method of the CNC machine tool. The process is as follows:

[0095] Step 1: The equipment information management module inputs the basic information, the spindle system information, the feed shaft system information, the automatic tool automatic tool changing system information, and the auxiliary system information of the CNC machine tool to be tested, as follows:

TABLE-US-00001 TABLE 1 Basic information on CNC machine tool Basic information of CNC machine tool Machine type Vertical CNC milling center Machine model XK714D CNC system FANU

TABLE-US-00002 TABLE 2 Spindle system information Spindle system information Rated speed r/min 750 Minimum rotation speed r/min 60 Maximum rotation speed r/min 80000 Speed interval above rated speed r/min 250 Speed interval below rated speed r/min 50

TABLE-US-00003 TABLE 3 Feed shaft system information Feed shaft system information X Y Z Rapid feed&return speed mm/min 20000 20000 15000 Cutting feed rate mm/min 5000 5000 5000 Range of motion mm 400 300 300

TABLE-US-00004 TABLE 4 Automatic Tool automatic tool changing system information The maximum number of cutting tools in 12 Automatic Tool Changing System

TABLE-US-00005 TABLE 5 Auxiliary System Information Auxiliary System Information 1 2 3 Name compressed air cutting fluid Chip system system transportation system Start code M07 M08 M45 Stop code M09 M09 M46

[0096] Step 2: Set the test parameters in the test parameter setting module, comprising the following test parameters: state-running duration t.sub.o, state switching flag time t.sub.d, state-switching-delay duration t.sub.dc, the number of test samples N.sub.tr, and the initial distance of feed shaft d.sub.I-o; As shown in Table 6:

TABLE-US-00006 TABLE 6 Parameter information Test parameter State-running duration t.sub.o 10 s State-switching-mark duration t.sub.d 5 s State-switching-delay duration t.sub.dc 20 s Number of training samples N.sub.tr 10 s Initial distance d.sub.I-o of feed shaft i 50 mm

[0097] Step 3: According to the test parameters in Step 2 and the basic information, active power system information, and auxiliary system information of CNC machine tool in Step 1, the NC code of the test experiment and the NC code of the verification experiment are generated in the NC code generation module; NC code generation rules for each axis are consistent, and here, for example only X-axis:

[0098] The NC code generation rule of the spindle system is:

[0099] i values from 1 to 2 N.sub.tr, generating the following code in turn: [0100] G97 S ns.sub.S-tr[i]; [0101] M03; [0102] G04 X t.sub.o; [0103] M05; [0104] G04 X t.sub.o

[0105] Among them, the code ‘G97 S ns.sub.S-tr[i]’ indicates that the spindle starting speed is set to ns.sub.S-tr[i] The code ‘M03’ represents the spindle forward start: The code ‘M05’ indicates that the spindle stops rotating; The code ‘G04 X t.sub.o’ means running the next code after the t.sub.o is suspended; At the same time, the calculation formula of ns.sub.S-tr[i] is as follows:

[00010]$\begin{array}{cc}n{s}_{S-\mathrm{tr}}\left[i\right]=\{\begin{array}{c}{n}_{S-r}-{\Delta }_{S-n-l}.Math.\frac{n_{S-\min }-n_{S-r}}{N_{tr}{\Delta }_{S-n-l}}.Math.\left(i-N_{tr}\right),0\le i<N_{tr}\\ n_{S-r}+{\Delta }_{S-n-u}.Math.\frac{n_{S-\max }-n_{S-r}}{N_{tr}{\Delta }_{S-n-u}}.Math.\left(i-N_{tr}\right),N_{tr}\le i\le 2N_{tr}\end{array};& \end{array}$

[0106] The total number of speed tests is 2 N.sub.tr; 0<i≤2N.sub.tr; Where └ ┘ denotes the integer operation;

[0107] NC code generation rules of the feed shaft system include rapid feed&return test code generation rules and cutting feed generation rules, rapid feed&return test code generation rules are as follows:

[0108] i values from 1 to 2 N.sub.tr, generating the following code in turn: [0109] G00 X d.sub.X-tr[i]; [0110] G04 X t.sub.o; [0111] G04X d.sub.X-o; [0112] G04 X t.sub.d

[0113] Where ‘G00X d.sub.X-tr[i]’ denotes the specified distance d.sub.X-tr[i] for fast forward and backward in the X direction; At the same time, the calculation formula of d.sub.X-tr[i] is as follows:

[00011]$\begin{array}{cc}{d}_{X-\mathrm{tr}}\left[i\right]=d_{X-o}+\frac{d_X-2d_{X-o}}{N_{tr}}i;& \end{array}$

[0114] Among them, the total number of tests of rapid feed&return specified distance is N.sub.tr; 1≤i≤N.sub.tr.

[0115] The generation rules for cutting feed test code are as follows:

[0116] i values from 1 to 2 N.sub.tr, generating the following code in turn: [0117] G01 X d.sub.X-tr[j]F cv.sub.X-tr[i]; [0118] G04 X t.sub.o; [0119] G00 X d.sub.X-oF cv.sub.X-tr[i]; [0120] G04 X t.sub.d

[0121] The code ‘G01X d.sub.X-tr[j] F cv.sub.X-tr[i]’ represents the specified feeding speed cvX-tr[i] to the direction of the X-axis dX-tr[i] The calculation formula CvX-tr[i] is as follows:

[00012]${\mathrm{cv}}_{X-\mathrm{tr}}\left[i\right]=\frac{{\mathrm{cv}}_{X-\max }}{N_{\mathrm{tr}}}i;$

[0122] The total number of cutting feed rate tests is N.sub.tr, 1≤i≤N.sub.tr; cv.sub.X-max represents the maximum cutting feed rate in the X-axis direction.

[0123] The NC code generation rules for testing the automatic tool changing system are:

[0124] i values from 2 to N.sub.t, generating the following code in turn: [0125] Ti; M06; [0126] G04 X t.sub.o; [0127] T1; M06 [0128] G04 X t.sub.d

[0129] Among them, code ‘Ti; M06’ means replacing the i-th knife.

[0130] NC code generation rules for other auxiliary systems of machine tool control start and stop are:

[0131] i values from 1 to N.sub.Au, generating the following code in turn: [0132] M NC.sub.i-s; [0133] G04 X t.sub.0; [0134] M NC.sub.t-p; [0135] G04 X t.sub.d

[0136] Among them, M NC.sub.i-s represents the startup NC code of the other auxiliary system i (1≤i≤N.sub.Au), M NC.sub.t-p represents the shutdown NC code of the other auxiliary system i.

[0137] The NC code generation rules for validation experiments are:

[0138] i values from 1 to N.sub.Au, generating the following code in turn: [0139] M NC.sub.i-s

[0140] After the above code, add the following code:

[00013]$G04;X{t}_o;$$T\frac{N_t}{2};M06;$$G04X{t}_d$$G00{\mathrm{Xd}}_{X-o}{\mathrm{Yd}}_{Y-o}{\mathrm{Zd}}_{Z-o}$$G04;X{t}_d;$$G97S{n}_{s-r};$$M03;$$G04;X{t}_d;$$G00X2d_{X-o}Y2d_{Y-o}Z2d_{Z-o}F{\mathrm{cv}}_{\max }$$G04;X{t}_d;$$M05;$

[0141] Then, add the code generated by the following rules:

[0142] i values from 1 to N.sub.Au, generating the following code in turn: [0143] M NC.sub.i-p

[0144] Specific NC code generated according to the above NC code generation rules can be referred to in Table 7-8 below:

TABLE-US-00007 TABLE 7 NC code preview 1 Auxiliary system Spindle test Rapid feed&return detection (partial code) (partial code) M07; . . . . . . G04 X10; G97 S700; G00 X80; M09; M03; G04 X10; G04 X5; G04 X10; G00 X50; M08; M05; G04 X5; G04 X10; G04 X5; G00 X110; M09; G97 S750; G04 X10; G04 X5; M03; G00 X50; M45; G04 X10; G04 X5; G04 X10; M05; G00 X140; M46; G04 X5; G04 X10; G04 X5; . . . . . .

TABLE-US-00008 TABLE 8 NC code preview 2 Automatic Tool Cutting feed Change (Part of Confirmatory (partial code) Code) experiment . . . . . . M07; M08; M45 G01 X80 F2000; T3; M06; G04; X5; G04 X10; G04 X10; T6; M06; G01 X50 F2000; T1 M06; G04; X5; G04 X5; G04 X5; G00 X100 Y100 Z100 G01 X110 F2500; T4 M06; G04; X5; G04 X10; G04 X10; G00 X250 Y200 Z200 G01 X50 F2500; T1; M06; G04; X5; G04 X5; G04 X5; G97 S750; G01 X140 F3000; T5; M06; M03; G04 X10; G04 X10; G04; X5; G01 X100 Y100 Z100 F5000 . . . . . . G04; X5; M05; M09; M46

[0145] Step 4: Set the power sensors parameters in the field test module, and input the NC code of the test experiment and the NC code of the verification experiment into the CNC system of the machine tool to be tested; Power sensors parameters are shown in Table 9:

TABLE-US-00009 TABLE 9 Power sensors information Power sensors information Model HC-33C3 Sampling frequency fs 20 HZ Precision ±0.5%

[0146] Field test management module is also used to generate field test steps comprising the following steps according to test parameters:

[0147] Step X. 1: The power sensors is installed at the power input end of the CNC machine tool, and the installation position of the power sensors is shown in FIG. 2. The data output end of the power sensors is connected to the data input end of the field test management module through the conversion interface.

[0148] Step X. 2: Start the total power supply of the CNC machine tool and run continuously for t.sub.o time;

[0149] Step X. 3: Turn off the total power supply of CNC machine tool and wait for t.sub.d;

[0150] Step X. 4: Start the total power supply of the NC machine tool and run continuously for time;

[0151] Step X. 5: Start CNC system and wait for initialization;

[0152] Step X. 6: After the NC system is initialized, t.sub.o the duration shall be continuously run;

[0153] Step X. 7: Close the CNC system and wait for t.sub.d;

[0154] Step X. 8: Start the CNC system and wait for initialization;

[0155] Step X. 9: Input NC code generated by the NC code generation module into the CNC system;

[0156] Step X. 10: CNC system runs the NC code used to control the spindle system in the test experiment;

[0157] Step X. 11: Adjust the rate knob of the operating panel of the CNC machine tool, set the rate to 50%, and the CNC system runs the NC code used to control the operation of the feed shaft system in the test experiment.

[0158] Step X. 12: Adjust the rate knob of the operating panel of the CNC machine tool, set the rate to 100%, and the CNC system runs the NC code used to control the operation of the feed shaft system in the test experiment.

[0159] Step X. 13: CNC system runs the NC code used to control the operation of an automatic tool changing system in the test experiment;

[0160] Step X. 14: NC system runs the NC code used to control other auxiliary systems controlled by machine tools in the test experiment;

[0161] Step X. 15: CNC system runs the NC code for verification experiment;

[0162] Step X. 16: Turn off the CNC system of the machine tool and turn off the power supply of the machine tool.

[0163] Step 5: CNC machine tool run the NC code of the test experiment and control the test machine tools to run according to the NC code of the test experiment. At the same time, the power data information of each energy consumption subsystem in the test experiment is collected by the power sensors installed on the test machine tools. The power set of each energy consumption subsystem in the test experiment is generated by the field test management module according to the power data information collected by the power sensors.

[0164] The test power set includes: the input power set of the CNC machine tool in the power supply opening process D.sub.SP, the input power set of the CNC machine tool in the numerical control system running process D.sub.SC, the input power set of the CNC machine tool in the i-th (1≤i≤N.sub.Au) auxiliary system running process D.sub.Au-i, the input power set of the CNC machine tool in the i-th (1≤i≤2 N.sub.tr) spindle speed starting process D.sub.PS-i, the input power set of the CNC machine tool in the i-th (1≤i≤2N.sub.tr) spindle speed running process D.sub.US-i, the input power set of the CNC machine tool in the feed shaft I (I∈[X, Y, Z]) in the rapid feed&return process with i-th (1≤i≤N.sub.tr) rapid feed&return distance D.sub.I-PF-i, the input power set of the CNC machine tool in the feed shaft I (I∈[X, Y, Z]) in the i-th (1≤i≤N.sub.tr) cutting feed speed cutting feed process D.sub.I-UF-i and the input power set of the CNC machine tool in the tool changing system for the n.sub.t (1≤n.sub.t≤N.sub.t) tool position tool process D.sub.PT-n.sub.t.

[0165] Step 6: CNC machine tool run the NC code of the verification experiment and control the CNC machine tool to be tested to run according to the NC code of the verification experiment. At the same time, the power data information of each energy consumption subsystem in the verification experiment is collected by the power sensors installed on the CNC machine tool to be tested. The power set of each energy consumption subsystem in the verification experiment is generated by the field test management module according to the power data information collected by the power sensors.

[0166] The verification experiment power set includes: the input power set D.sub.SC_V of NC machine tool in the operation process of NC system, the input power set D.sub.Au-V of NC machine tool in the whole operation process of the auxiliary system, the input power set D.sub.PS-V of NC machine tool in the spindle system starting at the specified speed n.sub.s-V, the input power set D.sub.US-V of CNC machine tool in the spindle system running at the specified speed n.sub.S-V, the input power set D.sub.I-PF-V of CNC machine tool in the feed shaft I (I∈[X, Y, Z]) in the process of rapid feed&return at the specified speed d.sub.I-V, the input power set D.sub.I-PF-V of CNC machine tool in the feed shaft I (I∈[X, Y, Z]) in the process of cutting at the specified feed speed cv.sub.I-V, and the input power set D.sub.PT-V of CNC machine tool in the automatic tool changing system.

[0167] Step 7: Data analysis module establishes the inherent energy efficiency factor function of each energy consumption subsystem according to the power set of the test experiment in step 5; The disclosure provides two function modeling methods: 1) discrete modeling method; 2) Discrete modeling method combined with an adaptive fitting modeling method; The following two function modeling methods are described respectively.

[0168] 1) Discrete Modeling Method

[0169] In the data analysis module, the discrete modeling method is used to establish the discrete function of the inherent energy efficiency factors of each energy consumption subsystem of CNC machine tool, and the discrete function of the inherent energy efficiency factors is used as the inherent energy efficiency factor function, comprising the following steps:

[0170] Step 2.1: Calculate the power average P.sub.D, run time t.sub.D, and energy consumption E.sub.D for each power set according to the following general formula:

[00014]$\stackrel{\_}{P_D}=\frac{1}{N_D}\underset{k=1}{\overset{N_D}{.Math.}}{P}_k,P_k\in D;$$t_D=\frac{N_D}{fs};$$E_D=\frac{1}{fs}\underset{k=1}{\overset{N_D}{.Math.}}{P}_k,P_k\in D$

[0171] Among them, D represents the power set, P.sub.k represents the k.sup.th element in the power set D, N.sub.D represents the number of elements in the power set D, and fs represents the sampling frequency of the power sensors;

[0172] Step 2.2: Set up the discrete function of inherent energy efficiency factors of energy consumption subsystem below CNC machine tool:

[0173] The power consumption of the energy dissipation subsystem of the CNC machine tool power supply opening process PSP:

P.sub.SP=P.sub.in-SP;

[0174] Among them, P.sub.in-SP represent the average power of the power set D.sub.SP of CNC machine tool in the power supply opening process.

[0175] Power Consumption of CNC Machine Tool Control System P.sub.Cs;

P.sub.Cs=P.sub.in-SC−P.sub.SP;

[0176] Among them, P.sub.in-SC represent the average power of the power set D.sub.SC of CNC machine tool in the CNC system operation process.

[0177] Power consumption of the CNC machine tool auxiliary system P.sub.Au-i:

P.sub.Au-i=P.sub.in-Au-t−P.sub.in-SC

[0178] Among them, P.sub.in-Au-i represent the power mean value of the power set D.sub.Au-i of CNC machine tool during the operation of the i (1≤i≤N.sub.Au) auxiliary system.

[0179] The discrete function of startup time and startup energy consumption of CNC machine tool spindle system C.sub.S-PS;

[00015]$C_{S-PS}=\left[\begin{array}{ccc}n& t_{S-\mathrm{PS}}& E_{S-PS}\end{array}\right]=\left[\begin{array}{ccc}{\mathrm{ns}}_{S-\mathrm{tr}}\left[1\right]& t_{S-\mathrm{PS}}\left[1\right]& E_{S-PS}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{ns}}_{S-\mathrm{tr}}\left[i\right]& t_{S-\mathrm{PS}}\left[i\right]& E_{S-PS}\left[i\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{ns}}_{S-\mathrm{tr}}\left[2N_{\mathrm{tr}}\right]& t_{S-\mathrm{PS}}\left[2N_{\mathrm{tr}}\right]& E_{S-PS}\left[2N_{\mathrm{tr}}\right]\end{array}\right]$

[0180] Where n represents the spindle speed, which is the independent variable; Both t.sub.S-PS and E.sub.S-PS are dependent variables, which are the starting time of the spindle system and the energy consumption of the spindle system. The t.sub.S-PS [i] is the running time of the power set D.sub.PS-i in the starting process of the i (1≤i≤2N.sub.tr) spindle speed of CNC machine tool. The calculation formula E.sub.S-PS[i] is as follows:

E.sub.S-PS[i]=E.sub.in-PS[i]−P.sub.in-SCt.sub.S-PS[i]

[0181] Among them, E.sub.in-PS[i] is the set energy consumption of the power set D.sub.PS-i of the CNC machine tool during the starting process of the i (1≤i≤2N.sub.tr) spindle speed.

[0182] The idling power discrete function of the spindle system of CNC machine tool in the idling period is C.sub.S-US:

[00016]$C_{S-\mathrm{US}}=\left[\begin{array}{cc}n& P_{S-\mathrm{US}}\end{array}\right]=\left[\begin{array}{cc}{\mathrm{ns}}_{S-\mathrm{tr}}\left[1\right]& P_{S-\mathrm{US}}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{ns}}_{S-\mathrm{tr}}\left[i\right]& P_{S-\mathrm{US}}\left[i\right]\\ \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{ns}}_{S-\mathrm{tr}}\left[2N_{\mathrm{tr}}\right]& P_{S-\mathrm{US}}\left[2N_{\mathrm{tr}}\right]\end{array}\right]$

[0183] Where n represents the spindle speed, which is the independent variable; P.sub.S-the US represents the power function of the spindle system in the air operation process as the dependent variable; The calculation power formula of P.sub.S-US[i] is as follows:

P.sub.S-US[i]=P.sub.in-US−P.sub.in-SC

[0184] Where P.sub.in-US[i] is the mean power of the power set D.sub.US-i of the CNC machine tool in the i (1≤i≤2N.sub.tr) spindle speed empty operation process.

[0185] The discrete function C.sub.I-PF(I∈[X, Y, Z], X, Y, Z denotes the X-axis, Y-axis, and Z-axis of the machine tool respectively) of the feed axis system of the CNC machine tool in the process of rapid feed&return:

[00017]$C_{I-\mathrm{PF}}=\left[\begin{array}{ccc}{d}_{I-\mathrm{tr}}& t_{I-\mathrm{PF}}& E_{I-\mathrm{PS}}\end{array}\right]=\left[\begin{array}{ccc}{d}_{I-\mathrm{tr}}\left[1\right]& t_{I-\mathrm{PF}}\left[1\right]& E_{I-\mathrm{PS}}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ d_{I-\mathrm{tr}}\left[i\right]& t_{I-\mathrm{PF}}\left[i\right]& E_{I-\mathrm{PS}}\left[i\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ d_{I-\mathrm{tr}}\left[N_{\mathrm{tr}}\right]& t_{I-\mathrm{PF}}\left[N_{\mathrm{tr}}\right]& E_{I-\mathrm{PS}}\left[N_{\mathrm{tr}}\right]\end{array}\right]$

[0186] Among them, d.sub.I-tr represents rapid feed&return distance as independent variables; t.sub.I-PF and E.sub.I-PF are dependent variables, t.sub.I-PF means rapid feed&return time, E.sub.I-PF means rapid feed&return energy consumption; t.sub.I-PF[i] is the running time of feed axis I (I∈[X, Y, Z]) in the power set D.sub.I-PF-i with i.sup.th (1≤i≤N.sub.tr) rapid feed&return distance in rapid feed&return process; The calculation formula of E.sub.I-FF[i] is as follows:

E.sub.I-FF[i]=E.sub.in-I-PF[i]−P.sub.in-SCt.sub.I-PF[i];

[0187] In the formula, E.sub.in-I-PF[i] is the set energy consumption of the power set D.sub.I-PF-i of the feed axis I (I∈[X, Y, Z]) with the i.sup.th(1≤i≤N.sub.tr) fast forward fast-backward distance in the fast forward fast-backward process.

[0188] The discrete function C.sub.I-UF of feed power in the cutting feed operation process of the CNC machine tool feed shaft system is:

[00018]$C_{I-\mathrm{UF}}=\left[\begin{array}{cc}{\mathrm{cv}}_{I-\mathrm{TR}}& P_{I-\mathrm{UF}}\end{array}\right]=\left[\begin{array}{cc}{\mathrm{cv}}_{I-\mathrm{TR}}\left[1\right]& P_{I-\mathrm{UF}}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{cv}}_{I-\mathrm{TR}}\left[i\right]& P_{I-\mathrm{UF}}\left[i\right]\\ \mathrm{.Math.}& \mathrm{.Math.}\\ {\mathrm{cv}}_{I-\mathrm{TR}}\left[N_{\mathrm{tr}}\right]& P_{I-\mathrm{UF}}\left[N_{\mathrm{tr}}\right]\end{array}\right]$

[0189] Where cv.sub.I-tr represents the cutting feed rate of I axis, which is the independent variable; P.sub.I-UF represents the feed power of I-axis cutting as the dependent variable;

[0190] The calculation formula of P.sub.I-UF[i] is as follows:

P.sub.I-UF[i]=P.sub.in-I-UF[i]−P.sub.in-SC

[0191] In the formula, P.sub.in-I-UF[i] is the mean power of D.sub.I-UF-i is the power set of the CNC machine tool in the feed axis I (I∈[X, Y, Z]) cutting feed process at the i (1≤i≤N.sub.tr) cutting speed.

[0192] The discrete function C.sub.PT of the tool changing time and energy consumption of the CNC tool changing system during the tool changing operation is:

[00019]$C_{PT}=\left[\begin{array}{ccc}{n}_t& t_{\mathrm{PT}}& E_{PT}\end{array}\right]=\left[\begin{array}{ccc}1& t_{\mathrm{PT}}\left[1\right]& E_{\mathrm{PT}}\left[1\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ n_t& t_{\mathrm{PT}}\left[n_t\right]& E_{\mathrm{PT}}\left[n_t\right]\\ \mathrm{.Math.}& \mathrm{.Math.}& \mathrm{.Math.}\\ N_t& t_{\mathrm{PT}}\left[N_{\mathrm{tr}}\right]& E_{\mathrm{PT}}\left[N_{\mathrm{tr}}\right]\end{array}\right]$

[0193] Where n.sub.t represents the tool change position, which is the independent variable; t.sub.PT and E.sub.PT are dependent variables, t.sub.PT means tool change time, E.sub.PT means tool change energy consumption of tool changing system; t.sub.PT[n.sub.t] is the running time for the power set D.sub.PT-nt of the tool process of the automatic tool changing system for the n.sub.t.sup.th (1≤n.sub.t≤N.sub.t) tool position; The calculation formula of E.sub.PT[n.sub.t] is as follows:

E.sub.PT[n.sub.t]=E.sub.in-PT[n.sub.t]−P.sub.in-SCt.sub.PT[n.sub.t]

[0194] In the expression, E.sub.in-PT[n.sub.t] is the set energy consumption of the power set D.sub.PT-n.sub.t of the tool process of the automatic tool changing system for the n.sub.t.sup.th (1≤n.sub.t≤N.sub.t) tool position.

[0195] The discrete function of inherent energy efficiency factors can be transformed into tabular form to facilitate the calculation and query of inherent energy efficiency factors. For example, the discrete function of inherent energy efficiency factors C.sub.S-PS for the spindle system startup process and air operation process of the vertical CNC milling center XK714D in this specific implementation method is transformed into the following table 10:

TABLE-US-00010 TABLE 10 The discrete function C.sub.S − PS of inherent energy efficiency factors in the starting process of CNC machine tool spindle system Rotational speed(r/min) 250 300 350 400 450 500 550 600 650 700 Duration(s) 0.75 0.75 0.75 0.75 0.8 0.8 0.8 0.8 0.8 0.8 Energy 307 390 412 468 603 692 690 801 817 831 consumption(J) Rotational speed(r/min) 750 1500 2250 3000 3750 4500 5250 6000 6750 7500 Duration(s) 0.8 0.85 0.90 0.95 1.15 1.35 1.75 2.15 2.6 3.4 Energy 967 1757 2615 3365 4703 6257 8249 10393 13089 16464 consumption(J)

[0196] Through the table, you can directly query the operating conditions for speed 250r/min, 300r/min, 350r/min . . . 7500r/min, the starting time t.sub.S-PS and the starting energy consumption E.sub.S-PS of the inherent energy efficiency factors in the spindle starting the process; For example, under the operating condition of 250 r/min rotational speed, the starting time of the spindle starting process is 0.75 s, and the starting energy consumption of the spindle starting process is 307 J. If it is necessary to calculate the inherent energy efficiency factors of the spindle system in the starting process of untested speed n.sub.un, the two speeds n.sub.+1 and n.sub.−1 closest to nun are first found in the tested speed nun; Then, 0.5(E.sub.S-PS (n.sub.+1)+E.sub.S-PS (n.sub.+1)) is approximated to the start-up energy consumption corresponding to the unknown speed nun, and 0.5(t.sub.S-PS (n.sub.+1)+t.sub.S-PS (n.sub.+l)) is approximated to the start-up time corresponding to the unknown speed nun.

[0197] 2) Discrete Modeling Method Combined with an Adaptive Fitting Modeling Method

[0198] Based on discrete modeling, the data analysis module uses the adaptive fitting modeling method to establish the inherent energy efficiency factor function of each energy consumption subsystem of CNC machine tool, comprising the following steps:

[0199] Step 2.3: The discrete function of each energy consumption subsystem of CNC machine tool established by discrete modeling method is regarded as a data set for fitting; Split the dataset into the training set and test set; 70% of the discrete data in the discrete function of the inherent energy efficiency factor can be used as the training set, and the rest as the test set;

[0200] Step 2.4: First fitting function and the second fitting function are fitted by the least square method according to the training set. The general formula of the first fitting function and the second fitting function is as follows:

[00020]$y\left(x,\omega \right)=\underset{j=0}{\overset{q}{.Math.}}{w}_j{x}^j,q=1,2;$

[0201] Step 2.5: Using the test set, the errors of the first fitting function, and the second fitting function are calculated according to the error function; The error function is as follows:

[00021]$E\left(\omega \right)=\frac{1}{2}\underset{n=1}{\overset{N}{.Math.}}{\left\{y\left(x_n,\omega \right)-y_n\right\}}^2$

[0202] Step 2.6: If the error E(ω.sub.q=1) of the first fitting function is smaller than that of the second fitting function E(ω.sub.q=2), the first fitting function is selected as the function of the inherent energy efficiency factor; Otherwise, the quadratic fitting function is selected as the inherent energy efficiency factor function.

[0203] For the vertical CNC milling center XK714D, the inherent energy efficiency factor function of the CNC machine tool energy consumption subsystem obtained by the adaptive fitting modeling method is as follows:

[0204] The spindle system of CNC machine tool: starting time fitting function t.sub.S-PS (n), starting energy consumption fitting function E.sub.S-PS (n), and air power fitting function P.sub.S-US (n) with spindle speed n as independent variables;

[00022]$t_{S-\mathrm{PS}}\left(n\right)=\{\begin{array}{c}0\mathrm{.8},60\le n\le 8000\\ 7\times 10^{-8}{n}^2-2\times 10^{-4}n+1.02,750\le n\le 8000,\left(R^2=0.995\right)\end{array}{E}_{S-\mathrm{PS}}\left(n\right)=\{\begin{array}{c}9\times 10^{-4}{n}^2+2.16n-192,60<n<750,\left(R^2=0.973\right)\\ 3\times 10^{-4}{n}^2-0.15n+1205,750\le n\le 8000,\left(R^2=0.999\right)\end{array}{E}_{S-\mathrm{US}}\left(n\right)=\{\begin{array}{c}2\times 10^{-1}{n}^2-0.063n+113,60<n<750,\left(R^2=0.973\right)\\ 4\times 10^{-6}{n}^2+0.0176n+222,750\le n\le 8000,\left(R^2=0.962\right)\end{array}$

[0205] CNC machine tool feed axis system: The fitting function t.sub.I-PF(d.sub.I) of rapid feed&return time length of X-axis, Y-axis, and Z-axis of CNC machine tool with the rapid feed&return distance d.sub.I as the independent variable, rapid feed&return energy consumption fitting function E.sub.I-PF(d.sub.I); The following I∈[X, Y, Z], X, Y, Z represent the X-axis, Y-axis, and Z-axis of the machine tool respectively.

[0206] When the ratio is 100%, the fitting results of the test experiment are as follows: X-axis rapid feed&return energy consumption fitting function E.sub.X-PF(d.sub.X) and rapid feed&return time fitting function t.sub.X-PF(d.sub.X), Y-axis forward and fast backward energy consumption fitting function E.sub.Y-PF(d.sub.Y) and rapid feed&return time fitting function t.sub.Y-PF (d.sub.Y), Z-axis rapid feed&return energy consumption fitting function E.sub.Z-PF(d.sub.Z) and rapid feed&return time fitting function t.sub.Z-PF(d.sub.Z).

[00023]$E_{X-PF}\left(d_X\right)=\{\begin{array}{c}0.8345d_X+29,\mathrm{positive}\mathrm{direction}\left(R^2=0.989\right)\\ 0.7299d_X+28,\mathrm{negative}\mathrm{direction}\left(R^2=0.990\right)\end{array};t_{X-FF}\left(d_X\right)=\{\begin{array}{c}3\times 10^{-3}{d}_X+0.06,\mathrm{positive}\mathrm{direction}\left(R^2=0.993\right)\\ 2.8\times 10^{-3}{d}_X+0.05,\mathrm{negative}\mathrm{direction}\left(R^2=0.979\right)\end{array};E_{Y-FF}\left(d_Y\right)=\{\begin{array}{c}-2.5\times 10^{-3}{d}_Y^2+1.3195d_Y+3,\mathrm{positive}\mathrm{direction}\left(R^2=0.989\right)\\ -3.6\times 10^{-3}{d}_Y^2+1.6404d_Y-8,\mathrm{positive}\mathrm{direction}\left(R^2=0.979\right)\end{array};t_{Y-FF}\left(d_Y\right)=\{\begin{array}{c}2.6\times 10^{-3}{d}_Y+0.09,\mathrm{positive}\mathrm{direction}\left(R^2=0.993\right)\\ 2.7\times 10^{-3}{d}_Y+0.07,\mathrm{negative}\mathrm{direction}\left(R^2=0.988\right)\end{array};E_{Z-FF}\left(d_Z\right)=\{\begin{array}{c}1.2015d_Z+11,\mathrm{positive}\mathrm{direction}\left(R^2=0.993\right)\\ 0.4376d_Z+15,\mathrm{negative}\mathrm{direction}\left(R^2=0.990\right)\end{array};t_{Z-PF}\left(d_Y\right)=\{\begin{array}{c}3.5\times 10^{-3}{d}_Z+0.11,\mathrm{positive}\mathrm{direction}\left(R^2=0.995\right)\\ 7\times 10^{-6}{d}_Z^2+2.4\times 10^{-3}{d}_Z+0.07,\mathrm{negative}\mathrm{direction}\left(R^2=0.985\right)\end{array};$

[0207] When the ratio is 50%, the fitting results of the test experiment are as follows: X-axis fast forward fast-backward energy consumption fitting function E.sub.X-PF(d.sub.X) and fast forward fast-backward time fitting function t.sub.X-PF(d.sub.X), Y-axis fast forward fast-backward energy consumption fitting function E.sub.Y-PF(d.sub.Y) and fast forward fast-backward time fitting function t.sub.Y-PF(d.sub.Y), Z-axis fast forward fast-backward energy consumption fitting function E.sub.Z-PF(d.sub.Z) and fast forward fast-backward time fitting function t.sub.Z-PF(d.sub.Z).

[00024]$E_{X-PF}\left(d_X\right)=\{\begin{array}{c}0.8263d_X+8,\mathrm{positive}\mathrm{direction}\left(R^2=0.998\right)\\ 0.7469d_X+7,\mathrm{negative}\mathrm{direction}\left(R^2=0.998\right)\end{array}{t}_{X-\mathrm{PF}}\left(d_X\right)=\{\begin{array}{c}5.8\times 10^{-3}{d}_X+0.07,\mathrm{positive}\mathrm{direction}\left(R^2=0.998\right)\\ 6.1\times 10^{-3}{d}_X+0.01,\mathrm{negative}\mathrm{direction}\left(R^2=0.996\right)\end{array}{E}_{Y-PF}\left(d_Y\right)=\{\begin{array}{c}0.6961d_Y+4,\mathrm{positive}\mathrm{direction}\left(R^2=0.994\right)\\ 0.7378d_Y+6,\mathrm{positive}\mathrm{direction}\left(R^2=0.993\right)\end{array}{t}_{Y-\mathrm{PF}}\left(d_Y\right)=\{\begin{array}{c}6.3\times 10^{-3}{d}_Y+0.01,\mathrm{positive}\mathrm{direction}\left(R^2=0.992\right)\\ 6.1\times 10^{-3}{d}_Y+0.03,\mathrm{positive}\mathrm{direction}\left(R^2=0.992\right)\end{array}{E}_{Z-\mathrm{PF}}\left(d_Z\right)=\{\begin{array}{c}1.4358d_Z+3,\mathrm{positive}\mathrm{direction}\left(R^2=0.998\right)\\ 0.2869d_Z+7,\mathrm{negative}\mathrm{direction}\left(R^2=0.987\right)\end{array}{t}_{Z-PF}\left(d_Z\right)=\{\begin{array}{c}8.1\times 10^{-3}{d}_Z+0.08,\mathrm{positive}\mathrm{direction}\left(R^2=0.997\right)\\ 8.2\times 10^{-3}{d}_Z+0.01,\mathrm{negative}\mathrm{direction}\left(R^2=0.998\right)\end{array}$

[0208] CNC machine tool feed axis system: cutting feed power fitting function P.sub.I-UF(cv.sub.I) with cutting feed speed cv.sub.I as the independent variable; where the subscript I∈[X, Y, Z], and X, Y, Z represent the X-axis, Y-axis, and Z-axis of the machine tool, respectively.

[00025]$P_{X-UF}\left(c{v}_X\right)=\{\begin{array}{c}0.0109c{v}_X+20,\mathrm{positive}\mathrm{direction}\left(R^2=0.997\right)\\ 0.0098c{v}_X,\mathrm{negative}\mathrm{direction}\left(R^2=0.994\right)\end{array}{P}_{Y-UF}\left(c{v}_Y\right)=\{\begin{array}{c}0.0107c{v}_Y,\mathrm{positive}\mathrm{direction}\left(R^2=0.977\right)\\ 0.0115c{v}_Y+4,\mathrm{negative}\mathrm{direction}\left(R^2=0.979\right)\end{array}{P}_{Z-UF}\left(c{v}_Z\right)=\{\begin{array}{c}0.0212c{v}_Z+22,\mathrm{positive}\mathrm{direction}\left(R^2=0.996\right)\\ 0.0048c{v}_Z+3,\mathrm{negative}\mathrm{direction}\left(R^2=0.965\right)\end{array}$

[0209] Automatic tool changing system for CNC machine tool: tool change time fitting function t.sub.PT(n.sub.t) and tool change energy consumption fitting function E.sub.PT(n.sub.t) with tool change data n t as independent variables.

[00026]$t_{PT}\left(n_t\right)=\{\begin{array}{c}1.257n_t+8.9,\mathrm{positive}\mathrm{rotations}\left(R^2=0.999\right)\\ 1.2n_t+9.2,\mathrm{negative}\mathrm{rotations}\left(R^2=0.999\right)\end{array}{E}_{\mathrm{PT}}\left(n_t\right)=\{\begin{array}{c}130n_t+2424,\mathrm{positive}\mathrm{rotations}\left(R^2=0.984\right)\\ 229n_t+2179,\mathrm{negative}\mathrm{rotations}\left(R^2=0.987\right)\end{array}$

[0210] Step 8: Data analysis module establishes the inherent energy efficiency factor function of CNC machine tool in each operation stage according to the inherent energy efficiency factor function of each energy consumption subsystem in step 7;

[0211] Power function of CNC machine tool in the standby stage P.sub.in-S(s.sub.Cs):

P.sub.in-PA(s.sub.cf,s.sub.ct)=250s.sub.Cs270;

[0212] The power function P.sub.in-PA(s.sub.cf,s.sub.ct) of CNC machine tool in the opening stage of the auxiliary system:

P.sub.in-PA(s.sub.cf,s.sub.ct)=530+220s.sub.cf+170s.sub.ct;

[0213] Among them, s.sub.cf indicates the state of the cutting fluid system, s.sub.cf=0 indicates that the cutting fluid system is closed, and s.sub.cf=1 indicates that the cutting fluid system is open. s.sub.ct means the state of the chip transport system, s.sub.ct=0 means the chip transport system is closed, s.sub.ct=1 means the chip transport system is open; Besides, the opening and closing of compressed air systems do not result in significant power consumption, which is omitted.

[0214] Energy Consumption Function of CNC Machine Tool in Starting Stage of Spindle System E.sub.in-PS(n, s.sub.cf, s.sub.ct):

E.sub.in-PS(n,s.sub.cf,s.sub.ct):=(530+220s.sub.cf+170s.sub.ct)t.sub.S-PS(n)+E.sub.S-PS(n);

[0215] Energy Consumption Function of CNC machine tool in the Stage of Fast Forward and Backward E.sub.in-PF(d.sub.I,s.sub.cf,s.sub.ct):

E.sub.in-PF(d.sub.I,s.sub.cf,s.sub.ct)=(530+220s.sub.cf+170s.sub.ct)Σt.sub.I-PF(d.sub.I)+ΣE.sub.I-PF(d.sub.I);

[0216] Energy Consumption Function of CNC Machine Tool in Tool Changing Stage E.sub.in-PT(n.sub.t,s.sub.cf,s.sub.ct):

E.sub.in-PT(n.sub.t,s.sub.cf,s.sub.ct)=(530+220s.sub.cf+170s.sub.ct)t.sub.PT(n.sub.t)+E.sub.PT(n.sub.t);

[0217] Power function of CNC machine tool in the idling stage P.sub.in-U(n,cv.sub.I,s.sub.cf,s.sub.ct):

P.sub.in-U(n,cv.sub.I,s.sub.cf,s.sub.ct)=530+220s.sub.cf+170s.sub.ct+P.sub.S-US(n)+ΣP.sub.I-UF(cv.sub.I),

[0218] Step 9: Validity verification of the inherent energy efficiency factor function of CNC machine tool generated by the data analysis module is carried out through the validity verification module. If the inherent energy efficiency factor function passes the validity verification, step B10 is entered. If the inherent energy efficiency factor function fails to pass the validity verification, then repeat step B1 to step B9. If it still fails to pass the validity verification of the function, contact technical personnel to eliminate the problem.

[0219] In the process of validity verification, the data analysis module calculates the verification value of the inherent energy efficiency factor according to the general formula in Step 2.1, comprising the verification power P.sub.in-SC-V in the operation stage of the numerical control system, the verification power P.sub.in-Au-V in the whole operation stage of the auxiliary system, the verification energy consumption E.sub.in-PS-V in the startup stage of the spindle system, the verification power P.sub.in-US-V in the empty operation stage of the spindle system, the verification energy consumption E.sub.in-I-PF-V in the rapid feed&return stage of the feed shaft I (I∈[X, Y, Z]), the verification power P.sub.in-US-V the cutting feed stage of the feed shaft I (I∈[X, Y, Z]), and the verification energy consumption E.sub.in-PT-V in the tool changing stage of the automatic tool changing system.

[0220] Among them, P.sub.in-SC-V is to verify the average power of the power set D.sub.SC-V in the experiment; P.sub.in-Au-V is to verify the average power of the power set D.sub.Au-V in the experiment; E.sub.in-PS-V is the set energy consumption of the power set D.sub.PS-V in the verification experiment; P.sub.in-US-V is to verify the average power of the power set D.sub.US-V in the experiment; E.sub.in-I-PF-V is the set energy consumption of the power set D.sub.I-PF-V in the verification experiment (I∈[X, Y, Z]); P.sub.in-US-V is the mean power of the power set D.sub.I-UF-V in the verification experiment (I∈[X, Y, Z]); E.sub.in_PT-V is the set energy consumption of the power set D.sub.PT-V in the verification experiment.

[0221] The data analysis module calculates the predicted values of the inherent energy efficiency factors under the verification experimental conditions according to the operation conditions of the verification experiment and the inherent energy efficiency factor function of the CNC machine tool, comprising the predicted power P.sub.in-S(s.sub.Cs=1) in the operation stage of the CNC system, the predicted power P.sub.in-PA(s.sub.Au-i=1) in the whole operation stage of the auxiliary system, the predicted energy consumption E.sub.in-PS(n=n.sub.V,s.sub.Au-i=1) in the start stage of the spindle system, the predicted power P.sub.in-U(n=n.sub.V,cv.sub.I=0,s.sub.Au-i=1) in the rapid feed&return stage of the feed shaft I (I∈[X, Y, Z]), and the predicted power P.sub.in-U(n=0,cv.sub.I=cv.sub.I-V,s.sub.Au-i=1) in the cutting feed stage of the feed shaft I (I∈[X, Y, Z]), and the predicted energy consumption E.sub.in-PT(n.sub.t=n.sub.t-V,s.sub.Au-i=1) in the tool changing stage.

[0222] The relative error value is used to verify the effectiveness of the detected inherent energy efficiency elements of CNC machine tool in the standby operation process, the operation process of the auxiliary system, the startup process of the spindle system under the specified speed, the empty operation process of the spindle system under the specified speed, the feed shaft system in the specified fast feed and fast retreat process, the feed shaft system in the specified cutting feed process and the automatic tool changing system in the specified tool change process. If the relative error is within 5%, it can be regarded as a pass. The verification results of this specific implementation method are as follows:

TABLE-US-00011 TABLE 10 Validation Table Appointed Verification Acquiring parameters of value of system Awaiting proved validation verification predictive Relative Effective or stage experimental experiment value error not Standby stage s.sub.Cs = 1 528 W 520 W −1.52% effective Auxiliary system s.sub.cf = 1, s.sub.ct = 1 918 W 910 W −0.87% effective operation phase Starting stage of the .sup.S.sub.cf = 1, s.sub.ct = 1 925 J 928 J 0.32% effective spindle system n.sub.v = 750 r/min Running open stage of s.sub.cf = 1, s.sub.ct = 1 1114 W 1110 W −0.36% effective spindle system n.sub.v = 750 r/min cv.sub.I = 0 Fast forward and X-axis s.sub.cf = 1, s.sub.ct = 1 611 J 609 J −0.33% effective backward d.sub.X − v = 150 mm stage of feed Y-axis s.sub.cf = 1, s.sub.ct = 1 426 J 421 J −0.17% effective shaft system d.sub.Y − v = 100 mm Z-axis s.sub.cf = 1, s.sub.ct = 1 532 J 539 J 1.32% effective d.sub.Z − v = 100 mm The cutting feed X-axis s.sub.cf = 1, s.sub.ct = 1 997 W 986 W −1.10% effective stage of feed cv.sub.X − v = 5000 mm/min shaft system Y-axis s.sub.cf = 1, s.sub.ct = 1 997 W 986 W −1.10% effective cv.sub.Y − v = 5000 mm/min Z-axis s.sub.cf = 1, s.sub.ct = 1 956 W 962 W 0.63% effective cv.sub.Z − v = 5000 mm/min Tool changing stage of s.sub.cf = 1, s.sub.ct = 1 18228 J 18275 J 0.26% effective tool changing system n.sub.t − v = 6

[0223] It can be seen from Table 10 that the relative errors are all less than 5%. The inherent energy efficiency factor function (discrete modeling method combined with an adaptive fitting modeling method) of the specific implementation method has passed the validity verification. Therefore, the inherent energy efficiency factor function that has passed the validity verification can be used to calculate the inherent energy efficiency factor under the specified operating conditions, and step 10 can be entered.

[0224] Step 10: According to the inherent energy efficiency factor function of CNC machine tool in each operation process and the operating conditions of CNC machine tool, the inherent energy efficiency factors in each operation process of CNC machine tool are calculated.

[0225] According to the inherent energy efficiency factor function of each energy consumption subsystem in the test experiment and the operating conditions of each energy consumption subsystem in the test experiment, the inherent energy efficiency factors of each energy consumption subsystem in different operating conditions and different operating processes are calculated, to obtain the inherent energy efficiency factors of CNC machine tool to be tested. Operating conditions are composed of operating parameters, comprising the following operating parameters: CNC system operating state s.sub.Cs, spindle speed n, feed shaft fast feed fast retreat distance d.sub.I, feed shaft cutting speed cv.sub.I, tool change position n.sub.t, the auxiliary system operating state s.sub.Au-i.