CN106475908A - Follow grinding process lathe operation energy consumption Forecasting Methodology based on standard G code - Google Patents

Abstract

The invention discloses a kind of follow grinding process lathe operation energy consumption Forecasting Methodology based on standard G code.This method comprises the following steps：First, the movement velocity of the kinematic axiss according to involved by G code calculates each row machining code and the working time, then draw out using the time as abscissa, each motion axle speed is schemed as the V t of vertical coordinate.Then, the gear ratio according to each shaft transmission, is calculated the rotating speed of each axle motor, and the motion of each for lathe axle is rotated by motor to characterize.And draw out using the time as abscissa, motor speed is the n t figure of vertical coordinate.Then, using parameters such as known motor speed, Rated motor torques, calculate and draw out power of motor time dependent P t figure.Finally, lathe operation energy consumption can just be obtained through mathematic integral with summation.The present invention can before reality processing quantitative analyses lathe operation energy consumption, for processing scheme optimized choice provide foundation, to reach the effect of energy-saving and emission-reduction.

Description

Follow grinding process lathe operation energy consumption Forecasting Methodology based on standard G code

Technical field

The invention belongs to energy-conserving and environment-protective field is and in particular to a kind of follow grinding process lathe based on standard G code is transported Row energy consumption Forecasting Methodology.

Background technology

At present, industrial department carbon emission accounts for the 70% about of China's total emission volumn, and machine cut is the important of industry manufacture Ingredient, has tremendous influence to carbon emission amount.Lathe, as the of paramount importance master tool of industrial circle, is used for process In energy consumption evaluated, significance will be produced to energy-saving and emission-reduction.

Part mainly includes lathe standby energy consumption, workpiece cutting energy consumption and lathe operation in the energy consumption of course of processing lathe Energy consumption.Wherein lathe standby energy consumption is held essentially constant in process.Workpiece cutting energy consumption and lathe specification and workpiece The close relation such as material, but this part energy consumption is relatively low with respect to other two groups of magnitudes, for the impact very little of overall energy consumption.And Lathe operation energy consumption with work pieces process code relation closely, and affects on total energy consumption very greatly it is therefore necessary to find effective ways Lathe operation energy consumption is estimated.For lathe energy consumption, current research is concentrated mainly on lathe power consumption state and monitors in machine Field, is realized the collection of data, obtains lathe energy consumption after treated calculating by install sensor.However, due to sensor Costly, and more sensitive to the noise in environment, vibration, so all having with demographic data's disposal ability to Financial cost Very high requirement, just leads to lathe energy consumption cannot popularize in common process workshop in machine monitoring system.

Mainly for assessment of lathe operation energy consumption, this part is lathe operation energy consumption Forecasting Methodology based on standard G code In total energy consumption, proportion is larger and changes obvious part.This method can depart from data collecting system, by standard G generation Code is analyzed calculating, you can obtain part course of processing lathe operation energy consumption.Realize before reality processing, quantitative understanding is whole The required energy consumption of each axle motion in the individual course of processing, the optimized choice for processing scheme provides strong foundation.

To sum up shown, the prediction of lathe operation energy consumption is realized energy-saving and emission-reduction to reduction carbon emission and can be played positive role, and With respect to energy-consumption monitoring system, there is the higher suitability.

Content of the invention

Present invention aims to the deficiency of prior art, a kind of follow grinding mistake based on standard G code is proposed The Forecasting Methodology of journey lathe operation energy consumption, can be before reality processing, and quantitative analyses lathe operation energy consumption, is the excellent of processing scheme Change and select to provide foundation, to reach the purpose of technical ability reduction of discharging.

For achieving the above object, the present invention adopts the following technical scheme that：Process G code by analyzing part, to part plus During work, lathe operation energy consumption is predicted.Comprise the following steps that：

Step one：Moved according to each axle of G Command Resolution, analyze and be calculated each axle of the course of processing and do exercises information, Including forms of motion, speed and time.During follow grinding, the resolution of velocity method of common G instruction is as follows：

1) position G00

Typical numerical control code：G00Ua Wc；

Numerical control code implication：Reached by specified location with the mobile cutter of quick translational speed.Cutter moves along the x-axis a, along Z Axle moves c.

Following two situations are divided into according to cutter path：

1. non-linkage interpolation positioning：Cutter is positioned to each axle respectively with each axle motion maximum speed.Tool motion path one As be not straight line.

In code implementation, the speed of each axle motion and time such as formula (1) are to formula (4) Suo Shi：

V_x=V_x-max(2)

V_z=V_z-max(4)

Wherein t_x、t_zRepresent that X-axis and Z axis move to the time required for target location, V by current location respectively_x、V_zPoint Do not represent the translational speed of X-axis and Z axis in the implementation procedure of this section of code, V_x-max、V_z-maxRepresent that Machine Manufacture business sets respectively The highest translational speed that fixed X-axis and Z axis can reach.

2. linkage interpolation positioning：Cutter reaches the point specified, positioning speed along straight line movement within the time the shortest Degree is less than the quickest translational speed of each axle.

IfThe then speed of each axle motion and the time such as formula (5) in code implementation To formula (8) Suo Shi：

V_x=V_x-max(6)

t_z=t_x(7)

Otherwise, in code implementation the speed of each axle motion and time such as formula (9) to formula (12) Suo Shi：

t_x=t_z(9)

V_z=V_z-max(12)

Wherein t_x、t_zRepresent that X-axis and Z axis move to the time required for target location, V by current location respectively_x、V_zPoint Do not represent the translational speed of X-axis and Z axis in the implementation procedure of this section of code, V_x-max、V_z-maxRepresent that Machine Manufacture business sets respectively The highest translational speed that fixed X-axis and Z axis can reach.

2) linear interpolation G01

For entering the single shaft linear interpolation of withdrawing, incision mill；For the special-shaped outline grinding such as the connecting rod neck of bent axle, cam C-X two-axle interlocking interpolation and the X-Z linkage for English class teaching.

1. single axial movement

Typical numerical control code：G01Ua Ff；

Numerical control code implication：Cutter moves a in X direction with feed speed f.

As shown in formula (13) along the feed speed of X-direction movement and time in code implementation：

V_x=f (14)

Wherein, t_xRepresent the movement time of X-axis, V_xRepresent the movement velocity of X-axis.This form equally can be applied mechanically to Z axis, C Axial motion is controlled, and feed speed asks method similar with X-direction with the time.

2. C-X linkage

Typical numerical control code：G91G01Xa CγFf；

Numerical control code implication：Cutter moves a in X direction with feed speed f, and C axle rotates γ.

Along the feed speed of each coordinate axess movement and time such as formula (15) to formula (18) institute in code implementation Show：

Wherein, t_x、t_cIt is respectively the movement time of X-axis and C axle, V_x、V_cSpeed for X-axis and C axle.

3. X-Z linkage

Typical numerical control code：G91 G01 Xa Zc Ff；

Numerical control code implication：Cutter moves a in X direction with feed speed f, moves c along Z-direction.

Along the feed speed of each coordinate axess movement and time such as formula (19) to formula (22) institute in code implementation Show：

Wherein, t_x、t_zIt is respectively the movement time of X-axis and Z axis, V_x、V_zSpeed for X-axis and C axle.

Step 2：Calculate and draw each axle movement velocity and change over V-t figure.By the calculating of step one, can obtain To every line code execution required time and each motion axle speed, (if in this time, certain root axle does not move then it is assumed that speed Degree is zero).Draw out further in process using the time as abscissa, each axle speed is as the V-t image of vertical coordinate (wherein, the speed unit of linear axis is mm/min, and the speed unit of rotary shaft is r/min).

Step 3：Calculate and draw out each work drive motor rotating speed and change over n-t figure.By drive mechanism gear ratio and fortune Moving axis speed substitutes into formula (23) or (24) and can calculate and be plotted in whole process, using the time as abscissa, respectively Rotating speed is as the n-t image of vertical coordinate.The calculating reference formula (23) of linear axis, the calculating reference formula (24) of rotary shaft.

Wherein, V is the translational speed of linear axis (X-axis or Z axis), P_hFor the helical pitch of ball screw, N is ball screw line Number, n is the rotating speed of rotary shaft (C axle or grinding wheel spindle), and i is drive mechanism gear ratio, n_mFor motor speed.

Step 4：Calculate and draw out each work drive motor power and change over P-t figure.Wherein, grinding wheel spindle motor Shown in power calculation algorithms such as formula (26), it is permanent torque output when less than flex point rotary speed movement, higher than flex point rotary speed movement When be constant power output.The calculating reference formula (25) of flex point rotating speed.Headstock revolution (C axle) AC servo motor, workbench are indulged Power calculation algorithms such as formula to movement (Z axis) AC servo motor, grinding carriage transverse shifting (X-axis) AC servo motor (27) shown in.

Wherein, f is ac frequency, N_poleFor motor pole number, n_inflectFor spindle motor flex point rotating speed, T is motor volume Determine moment of torsion, n_mFor motor speed, P₀For motor rated power, P is power of motor.

Step 5：Computer bed operating energy consumption.By each motor, power over time is integrated obtaining electricity in process Function consumes, and each spindle motor energy consumption is added the lathe operation energy consumption obtaining in whole process.

Wherein, P_spindleFor grinding wheel spindle motor power, P_{x_axis}For grinding carriage transverse shifting AC servo motor power, P_{z_axis}Vertically move AC servo motor power, P for workbench_{c_axis}Turn round AC servo motor power, t for the headstock₁、t₂Point Not Wei program beginning and end time.

The present invention compared with prior art, has and projects substantive distinguishing features and notable technological progress as follows：

The present invention departs from data collecting system, by being analyzed to standard G code calculating, you can obtain part processed Journey lathe operation energy consumption.Before reality processing, the required energy consumption of each axle motion in quantitative understanding whole process, for processing The optimized choice of scheme provides strong foundation.Lathe operation energy consumption prediction to reduce carbon emission, realize energy-saving and emission-reduction can play long-pending Pole acts on.

Brief description

Fig. 1 is the FB(flow block) based on the follow grinding process lathe operation energy consumption Forecasting Methodology of standard G code for the present invention.

Fig. 2 is embodiment code application workpiece figure.

Fig. 3 is embodiment code application workpiece pictorial diagram.

Fig. 4 is interpolation axle V-T figure during workpiece tarry matter.

Fig. 5 is interpolation axle n-t figure during workpiece tarry matter.

Fig. 6 is interpolation axle P-t figure during workpiece tarry matter.

Specific embodiment

The invention will be further described with reference to the accompanying drawings and examples, but the scope of application of the present invention is not limited to down State embodiment.In the case of without departing from the above-mentioned thought of the present invention, this technological know-how and customary means are departed from according to this area, does Go out various replacements and change, all should include within the scope of the present invention.

Embodiment one：The present invention is a kind of follow grinding process lathe operation energy consumption Forecasting Methodology based on standard G code, Referring to Fig. 1.Step is as follows：

(1) moved according to each axle of G Command Resolution, analyze and be calculated each axle of the course of processing and do exercises information, including Forms of motion, speed and time；

(2) draw each axle movement velocity and change over V-t figure：Entered according to each axle speed of G Command Resolution using step (1) Line code parses；It is calculated every line code execution required time and each motion axle speed, draw out using the time as horizontal seat Mark, each axle speed is as the V-t image of vertical coordinate；

(3) calculate and draw out each work drive motor rotating speed change over n-t figure；

(4) calculate and draw out each work drive motor power change over P-t figure；

(5) computer bed operating energy consumption.

Embodiment two：The present embodiment is essentially identical with embodiment one, and special feature is as follows：

(1) in described step (1), code analysis are carried out according to each axle speed of G Command Resolution：Common G instruction in follow grinding Resolution of velocity method as follows：

In described step (1), code analysis are carried out according to each axle speed of G Command Resolution；Common G instruction in follow grinding Resolution of velocity method is as follows：

1) position G00

Typical numerical control code：G00 Ua Wc；

Numerical control code implication：Reached by specified location with the mobile cutter of quick translational speed, cutter moves along the x-axis a, along Z Axle moves c；

Following two situations are divided into according to cutter path：

1. non-linkage interpolation positioning：Cutter is positioned to each axle respectively with each axle motion maximum speed, tool motion path one As be not straight line；

In code implementation, the speed of each axle motion and time such as following equation (1) are to formula (4) Suo Shi：

V_x=V_x-max(2)

V_z=V_z-max(4)

Wherein t_x、t_zRepresent that X-axis and Z axis move to the time required for target location, V by current location respectively_x、V_zPoint Do not represent the translational speed of X-axis and Z axis in the implementation procedure of this section of code, V_x-max、V_z-maxRepresent that Machine Manufacture business sets respectively The highest translational speed that fixed X-axis and Z axis can reach；

2. linkage interpolation positioning：Cutter reaches the point specified, positioning speed along straight line movement within the time the shortest Degree is less than the quick translational speed of each axle；

IfThe then speed of each axle motion and the time such as following public affairs in code implementation Formula (5) is to formula (8) Suo Shi：

V_x=V_x-max(6)

t_z=t_x(7)

Otherwise, in code implementation the speed of each axle motion and time such as following equation (9) to formula (12) Suo Shi：

t_x=t_z(9)

V_z=V_z-max(12)

Wherein t_x、t_zRepresent that X-axis and Z axis move to the time required for target location, V by current location respectively_x、V_zPoint Do not represent the translational speed of X-axis and Z axis in the implementation procedure of this section of code, V_x-max、V_z-maxRepresent that Machine Manufacture business sets respectively The highest translational speed that fixed X-axis and Z axis can reach；

2) linear interpolation G01

For entering the single shaft linear interpolation of withdrawing, incision mill；For the special-shaped outline grinding such as the connecting rod neck of bent axle, cam C-X two-axle interlocking interpolation and the X-Z linkage for English class teaching；

1. single axial movement

Typical numerical control code：G01Ua Ff；

Numerical control code implication：Cutter moves a in X direction with feed speed f；

Along the feed speed of X-direction movement and time such as following equation (13) to formula (14) institute in code implementation Show：

V_x=f (14)

Wherein, t_xRepresent the movement time of X-axis, V_xRepresent the movement velocity of X-axis, equally can apply mechanically this form to Z axis, C Axial motion is controlled, and feed speed asks method similar with X-direction with the time；

2. C-X linkage

Typical numerical control code：G91G01Xa CγFf；

Numerical control code implication：Cutter moves a in X direction with feed speed f, and C axle rotates γ；

Along the feed speed of each coordinate axess movement and time such as following equation (15) to formula in code implementation (18) shown in：

Wherein, t_x、t_cIt is respectively the movement time of X-axis and C axle, V_x、V_cSpeed for X-axis and C axle；

3. X-Z linkage

Typical numerical control code：G91 G01 Xa Zc Ff；

Numerical control code implication：Cutter moves a in X direction with feed speed f, moves c along Z-direction；

Along the feed speed of each coordinate axess movement and time such as following equation (19) to formula in code implementation (22) shown in：

Wherein, t_x、t_zIt is respectively the movement time of X-axis and Z axis, V_x、V_zSpeed for X-axis and C axle.

(2) in described step (3), motor working speed is obtained by kinematic axiss speed calculation, formula is as follows：

Wherein, V is the translational speed of linear axis (X-axis or Z axis), P_hFor the helical pitch of ball screw, N is ball screw line Number, n is the rotating speed of rotary shaft (C axle or grinding wheel spindle), and i is drive mechanism gear ratio, n_mFor motor speed.

(3) in described step (4), power of motor is calculated by motor working speed.The work(of grinding wheel spindle motor Shown in rate computational methods such as formula (26), it is permanent torque output when less than flex point rotary speed movement, during higher than flex point rotary speed movement For constant power output；The calculating reference formula (25) of flex point rotating speed；Headstock revolution (C axle) AC servo motor, workbench are longitudinally (Z axis) move AC servo motor, the power calculation algorithms such as formula (27) of the mobile AC servo motor of grinding carriage horizontal (X-axis) Shown：

Wherein, f is ac frequency, N_poleFor motor pole number, n_inflectFor spindle motor flex point rotating speed, T is motor volume Determine moment of torsion, n_mFor motor speed, P₀For motor rated power, P is power of motor.

(4) described step (5) computer bed operating energy consuming process, will each motor power over time enters in process Row integration obtains energy consumption of electrical machinery, and each spindle motor energy consumption is added the lathe operation energy consumption obtaining in whole process, such as formula (28) shown in：

Wherein, P_spindleFor grinding wheel spindle motor power, P_{x_axis}For grinding carriage transverse shifting AC servo motor power, P_{z_axis}Vertically move AC servo motor power, P for workbench_{c_axis}Turn round AC servo motor power, t for the headstock₁、t₂Point Not Wei program beginning and end time.

Embodiment three：The present embodiment is eccentric shaft-like work tarry matter program, is flat with the H367 grinding machine that Shanghai Machine Tool Plant produces Platform, using corundum wheel, diameter of work is 20.05mm, and eccentric throw is 2.15mm, referring to figure (2) and figure (3)；Grinding speed 40m/s, grinding wheel radius 489m.Table 2 show embodiment code, and code execution step is：First, C axle, X-axis and Z axis move respectively Move to initial Working position；Then carry out the tarry matter of workpiece；Finally, emery wheel moves along X-axis and returns to security bit.In whole process The motor being related to includes：Longitudinally (Z axis) are mobile for grinding wheel spindle motor, headstock gyration (C axle) AC servo motor, workbench Laterally (X-axis) moves AC servo motor for AC servo motor, grinding carriage.The parameter of electric machine is as follows：

Table 1 parameter of electric machine

Embodiment code is more tediously long, so intercepting part illustrates, as shown in table 2 herein：

Table 2 embodiment code (part)

Energy consumption prediction comprises the following steps：

1), judge the interpolation algorithm that every line code is used, select suitable formula generation by formula (1) to formula (22) Enter to calculate, calculate each motion axle speed and the working time involved by each row machining code.Result of calculation is as shown in the table, Lathe initial position is C180Z30.0X50.0；

Table 3 step one content

2) obtain using the time as abscissa, each axle speed is as the V-t image of vertical coordinate.Grinding wheel spindle is processed whole During speed be held constant at 1560r/min.Headstock revolution (C axle) speed is 30r/min in angular positioning, in grinding During be held constant at 10r/min.Longitudinal (Z axis) translational speed of workbench is only when Z-direction moves to Working position Time rises to 1000mm/min, remains stationary in other times.Laterally (X-axis) translational speed is complex for workbench.Workpiece light During mill shown in interpolation axle V-t Fig. 4；

3) understand with reference to lathe handbook, grinding wheel spindle motor drives emery wheel rotation (grinding wheel spindle rotation) by belt pulley, Gear ratio is 1.30.Headstock revolution AC servo motor reaches fixture by belt pulley, and being started building by clamp strap, (C axle turns for part rotation Dynamic), gear ratio is 0.56.Grinding carriage transverse shifting AC servo motor drives ball wire bar pair by yielding coupling, because of rolling Lunella coronata mother be bound up with grinding carriage, so that grinding carriage realizes feed motion (X-axis movement), ball screw is single thread, Helical pitch is 5mm.Workbench vertically moves AC servo motor and drives ball screw to select by shaft coupling, with driving and workbench The ball nut being bound up moves, and realizes workbench and vertically moves (Z axis movement), and ball screw is single thread, and helical pitch is 5mm.Above parameter is substituted in formula (23) or (24) and obtains the change in time of rotating speed during workpiece tarry matter for each motor Change.Wherein grinding wheel spindle motor keeps invariablenes turning speed in the course of the work is 1200r/min.Interpolation axle n- during workpiece tarry matter Shown in t Fig. 5；

4) understand and each parameter of electric machine with reference to lathe handbook, it is as follows that substitution formula (26) can be calculated each Rated motor torque Shown in table：

Table 4 Rated motor torque

The flex point rotating speed of grinding wheel spindle motor (4 poles, power frequency 50Hz) is 1500r/min, and grinding wheel spindle motor is in tarry matter During rotating speed remain 1200r/min, so understand grinding wheel spindle motor during workpiece tarry matter for permanent torque work.By Formula (26) can be calculated grinding wheel spindle, and power remains 6Kw in the process.Headstock revolution (C axle) AC servo motor, work Station vertically moves (Z axis) AC servo motor, grinding carriage transverse shifting (X-axis) AC servo motor also keeps in the process Permanent torque exports, and substitutes into each interpolation shaft power during formula (26) can be calculated workpiece tarry matter and changes over P-t Fig. 6 institute Show；

5) each power of motor is changed over function and substitute in formula (28), if program starts as timeorigin, then t₁ =0, t₂=15.958.It is calculated the tarry matter program of this eccentric shaft, in operation, the predictive value of energy consumption is 100.307KJ.

Claims (5)

1. the follow grinding process lathe operation energy consumption Forecasting Methodology based on standard G code it is characterised in that：By known number Control machining code is predicted to the lathe operation energy consumption of follow grinding process, and step is as follows：

(1) moved according to each axle of G Command Resolution, analyze and be calculated each axle of the course of processing and do exercises information, include moving Form, speed and time；

(2) draw each axle movement velocity and change over V-t figure：Using step (1), in generation, is carried out according to each axle speed of G Command Resolution Code parsing；It is calculated every line code execution required time and each motion axle speed, draw out using the time as abscissa, respectively Axle speed is as the V-t image of vertical coordinate；

(3) calculate and draw out each work drive motor rotating speed change over n-t figure；

(4) calculate and draw out each work drive motor power change over P-t figure；

(5) computer bed operating energy consumption.

2. the follow grinding process lathe operation energy consumption Forecasting Methodology based on standard G code according to claim 1, it is special Levy and be：In described step (1), code analysis are carried out according to each axle speed of G Command Resolution；The speed of common G instruction in follow grinding Degree decomposition method is as follows：

1) position G00

Typical numerical control code：G00 Ua Wc；

Numerical control code implication：Reached by specified location with the mobile cutter of quick translational speed, cutter moves along the x-axis a, moves along Z axis Dynamic c；

Following two situations are divided into according to cutter path：

1. non-linkage interpolation positioning：Cutter is positioned to each axle respectively with each axle motion maximum speed, and tool motion path is general not It is straight line；

In code implementation, the speed of each axle motion and time such as following equation (1) are to formula (4) Suo Shi：

$t_x=\frac{a}{V_{x-max}}---\left(1\right)$

V_x=V_x-max(2)

$t_z=\frac{c}{V_{z-max}}---\left(3\right)$

V_z=V_z-max(4)

Wherein t_x、t_zRepresent that X-axis and Z axis move to the time required for target location, V by current location respectively_x、V_zRepresent respectively The translational speed of X-axis and Z axis, V in the implementation procedure of this section of code_x-max、V_z-maxRepresent the X that Machine Manufacture business sets respectively The highest translational speed that axle and Z axis can reach；

2. linkage interpolation positioning：Cutter reaches, along straight line movement, the point specified within the time the shortest, and locating speed is not Exceed the quick translational speed of each axle；

IfThe then speed of each axle motion and the time such as following equation (5) in code implementation To formula (8) Suo Shi：

$t_x=\frac{a}{V_{x-max}}---\left(5\right)$

V_x=V_x-max(6)

t_z=t_x(7)

$V_z=V_x\×{\tan }^{-1}\frac{c}{a}---\left(8\right)$

Otherwise, in code implementation the speed of each axle motion and time such as following equation (9) to formula (12) Suo Shi：

t_x=t_z(9)

$V_x=\frac{V_z}{{\tan }^{-1}\frac{c}{a}}---\left(10\right)$

$t_z=\frac{c}{V_{z-max}}---\left(11\right)$

V_z=V_z-max(12)

Wherein t_x、t_zRepresent that X-axis and Z axis move to the time required for target location, V by current location respectively_x、V_zRepresent respectively The translational speed of X-axis and Z axis, V in the implementation procedure of this section of code_x-max、V_z-maxRepresent the X that Machine Manufacture business sets respectively The highest translational speed that axle and Z axis can reach；

2) linear interpolation G01

For entering the single shaft linear interpolation of withdrawing, incision mill；C-X for the special-shaped outline grinding such as the connecting rod neck of bent axle, cam Two-axle interlocking interpolation and the X-Z linkage for English class teaching；

1. single axial movement

Typical numerical control code：G01 Ua Ff；

Numerical control code implication：Cutter moves a in X direction with feed speed f；

Along the feed speed of X-direction movement and time such as following equation (13) to formula (14) Suo Shi in code implementation：

$t_x=\frac{a}{f}---\left(13\right)$

V_x=f (14)

Wherein, t_xRepresent the movement time of X-axis, V_xRepresent the movement velocity of X-axis, equally can apply mechanically this form to Z axis, C axle side To motion be controlled, feed speed asks method similar with X-direction with the time；

2. C-X linkage

Typical numerical control code：G91 G01 Xa C γ Ff；

Numerical control code implication：Cutter moves a in X direction with feed speed f, and C axle rotates γ；

Along the feed speed of each coordinate axess movement and time such as following equation (15) to formula (18) institute in code implementation Show：

$t_x=t_c=\frac{L}{f}---\left(15\right)$

$V_x=\frac{a}{t_x}---\left(16\right)$

$V_c=\frac{\mathrm{\γ}}{t_c}---\left(17\right)$

$L=\sqrt{a^2+{\mathrm{\γ}}^2}---\left(18\right)$

Wherein, t_x、t_cIt is respectively the movement time of X-axis and C axle, V_x、V_cSpeed for X-axis and C axle；

3. X-Z linkage

Typical numerical control code：G91 G01 Xa Zc Ff；

Numerical control code implication：Cutter moves a in X direction with feed speed f, moves c along Z-direction；

Along the feed speed of each coordinate axess movement and time such as following equation (19) to formula (22) institute in code implementation Show：

$t_x=t_z=\frac{L}{f}---\left(19\right)$

$V_x=\frac{a}{t_x}---\left(20\right)$

$V_z=\frac{c}{t_z}---\left(21\right)$

$L=\sqrt{a^2+c^2}---\left(22\right)$

Wherein, t_x、t_zIt is respectively the movement time of X-axis and Z axis, V_x、V_zSpeed for X-axis and C axle.

3. the follow grinding process lathe operation energy consumption Forecasting Methodology based on standard G code according to claim 1, it is special Levy and be：In described step (3), motor working speed is obtained by kinematic axiss speed calculation, formula is as follows：

$n_m=\frac{V}{P_h\×N}---\left(23\right)$

$n_m=\frac{n}{i}---\left(24\right)$

Wherein, V is the translational speed of linear axis (X-axis or Z axis), P_hFor the helical pitch of ball screw, N is ball screw line number, and n is The rotating speed of rotary shaft (C axle or grinding wheel spindle), i is drive mechanism gear ratio, n_mFor motor speed.

4. the follow grinding process lathe operation energy consumption Forecasting Methodology based on standard G code according to claim 1, it is special Levy and be：In described step (4), power of motor, the energy meter of grinding wheel spindle motor are calculated by motor working speed Shown in calculation method such as formula (26), it is permanent torque output when less than flex point rotary speed movement, higher than being perseverance during flex point rotary speed movement Power output；The calculating reference formula (25) of flex point rotating speed；Headstock revolution (C axle) AC servo motor, workbench are longitudinally (Z axis) Mobile AC servo motor, grinding carriage are laterally shown in the power calculation algorithms such as formula (27) of (X-axis) mobile AC servo motor：

$n_{\inf lect}=\frac{f\×120}{N_{pole}}---\left(25\right)$

$P=\left\{\begin{array}{cc}\frac{T\×n_m}{9550},& n_m\≤n_{\inf lect}\\ P_0,& n_m>n_{\inf lect}\end{array}\right.---\left(26\right)$

$P=\frac{T\×n_m}{9550}---\left(27\right)$

Wherein, f is ac frequency, N_poleFor motor pole number, n_inflectFor spindle motor flex point rotating speed, T is Rated motor torsion Square, n_mFor motor speed, P₀For motor rated power, P is power of motor.

5. the follow grinding process lathe operation energy consumption Forecasting Methodology based on standard G code according to claim 1, it is special Levy and be：Described step (5) computer bed operating energy consuming process, will each motor power over time is amassed in process Get energy consumption of electrical machinery, each spindle motor energy consumption is added the lathe operation energy consumption obtaining in whole process, such as formula (28) Shown：

$E_{run-time}={\∫}_{t_1}^{t_2}{P}_{spindle}dt+{\∫}_{t_1}^{t_2}{P}_{x\_axis}dt+{\∫}_{t_1}^{t_2}{P}_{z\_axis}dt+{\∫}_{t_1}^{t_2}{P}_{c\_axis}dt---\left(28\right)$

Wherein, P_spindleFor grinding wheel spindle motor power, P_{x_axis}For grinding carriage transverse shifting AC servo motor power, P_{z_axis} Vertically move AC servo motor power, P for workbench_{c_axis}Turn round AC servo motor power, t for the headstock₁、t₂It is respectively journey The time of sequence beginning and end.