CN106021179A - NC code-based linear cutting energy consumption prediction method - Google Patents
NC code-based linear cutting energy consumption prediction method Download PDFInfo
- Publication number
- CN106021179A CN106021179A CN201610328991.5A CN201610328991A CN106021179A CN 106021179 A CN106021179 A CN 106021179A CN 201610328991 A CN201610328991 A CN 201610328991A CN 106021179 A CN106021179 A CN 106021179A
- Authority
- CN
- China
- Prior art keywords
- energy consumption
- cutting
- thickness
- unit interval
- wire
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
Abstract
The invention discloses an NC code-based linear cutting energy consumption prediction method. The method comprises the following steps of 1, inputting a preprocessed part program code, a part material and the thickness; 2, analyzing the code to obtain parameters of empty wire cutting, linear cutting and arc cutting parts; 3, calling a basic database to obtain basic parameters; and 4, calling an energy consumption prediction model to obtain a predicted energy consumption value. According to the method, the required total energy consumption can be predicted through an NC code, the part material and the thickness. The basic database in the method is obtained through an experimental method.
Description
[technical field]
The present invention relates to the technical field of energy consumption prediction, a kind of line based on NC code cutting energy consumption Forecasting Methodology
Technical field.
[background technology]
Line cutting is that a kind of fine wire (referred to as wire electrode) utilizing continuous moving is made electrode and workpiece is carried out pulse
Spark discharge ablation metal, the processing method of excision forming.During Wire EDM work, carry out between wire electrode and workpiece
Pulsed discharge.The negative pole of the electrode wire connection pulse power, workpiece connects the positive pole of the pulse power.When carrying out an electric pulse, at electrode
Producing a spark discharge between silk and workpiece, central temperature at discharge channel is instantaneous may be up to more than 10000 DEG C, and high temperature makes
Workpiece metal melts, and even has a small amount of gasification, and high temperature also makes the working solution part between wire electrode and workpiece produce gasification, these
Working solution and metal vapors moment after gasification expand rapidly, and have the characteristic of blast.This thermal expansion and local microburst,
Dish out fusing and the metal material that gasified and realize workpiece material is carried out galvanic corrosion cutting processing.
The energy consumption main source of line cutting has host computer energy consumption ECO(Ecomputer), wire transporting energy consumption EW(Ewire-
Feed), water pump bath energy consumption EP(Epump), illuminating lamp energy consumption EL(Elamp) energy consumption of constant power and cutting energy consumption ECU
(Ecutting) energy consumption of Variable power.Energy consumption is cut for line, if can directly be predicted according to NC code, then in processing
Before, directly can be obtained energy consumption by program.
[summary of the invention]
The purpose of the present invention solves the problems of the prior art exactly, proposes a kind of line based on NC code cutting energy consumption pre-
Survey method, it is possible to by NC code, part material and thickness, it was predicted that go out required total energy consumption.Basic database in the present invention leads to
Cross laboratory method to obtain.
For achieving the above object, the present invention propose a kind of line based on NC code cutting energy consumption Forecasting Methodology, including with
Lower step:
The first step: input preprocessing parts program code, part material and thickness；
Second step: code analysis, it is thus achieved that empty wire, straight cuts, circular arc cutting each several part parameter；
3rd step: call basic database, it is thus achieved that underlying parameter；
4th step: call energy consumption forecast model, it is thus achieved that prediction of energy consumption value: empty wire energy consumption EKWPass throughSolve, wherein EDKWFor the empty wire energy consumption under the unit time, XKWFor empty wire length, vwWalk for sky
Silk speed；Straight cuts energy consumption EZCUPass throughSolve, wherein EDZCUThick for the unit time
Straight cuts energy consumption under Du, XZCUFor straight cuts total length, vCUFor cutting wire travelling speed, ηSFor velocity coeffficient, ηHFor thickness
Coefficient, ηCFor material coefficient；Circular arc cutting energy consumption EYCUPass throughSolve, wherein
EDYCUEnergy consumption, X is cut for the circular arc under unit time thicknessYCUTotal length, η is cut for circular arcQFor coefficient of curvature；Host computer energy
Consumption ECOPass throughSolve, wherein EDCOFor unit time host computer energy consumption；
Wire transporting energy consumption EWPass throughSolve, wherein EDWFor unit time wire transporting energy consumption；
Water pump bath energy consumption EPPass throughSolve, wherein EDPRush for unit time water pump
Water consumption；Illuminating lamp energy consumption ELPass throughSolving, wherein a is 0 or 1, and a is 0 table
Showing and do not turn on light, a is that 1 expression is turned on light, EDLFor unit time illuminating lamp energy consumption；Total energy consumption EZPass through EZ=EKW+EZCU+EYCU+ECO+EW
+EP+ELSolve.
As preferably, the code analysis method in described second step is: judge whether have in preprocessing parts program code
B, distinguishes and uses 3B format program processing method and ISO format program processing method, it is thus achieved that empty wire, straight cuts, circular arc are cut
Cut each several part parameter: empty wire length, straight cuts total length, circular arc cut total length, circular arc cut radius.
As preferably, the basic database in described 3rd step includes material coefficient ηC, thickness coefficient ηH, velocity coeffficient ηS、
Empty wire travelling speed vw, cutting wire travelling speed vCU, empty wire energy consumption E under the unit intervalDKW, straight line under unit interval thickness cuts
Cut energy consumption EDZCU, circular arc cutting energy consumption E under unit interval thicknessDYCU, coefficient of curvature ηQ, unit interval host computer energy consumption
EDCO, unit interval wire transporting energy consumption EDW, unit interval water pump bath energy consumption EDP, unit interval illuminating lamp energy consumption EDL, empty wire speed
Degree vw, cutting wire travelling speed vCUAcquisition according to lathe property, remaining parameter is all obtained by laboratory method, by data acquisition packaging
Putting the power curve obtaining wire cutting machine tool, the integration of power over time is energy consumption, and experimental procedure is:
The first step: open host computer, now measured energy consumption is E divided by the timeDCO；
Second step: opening water pump bath, the energy consumption increased in the unit interval is EDP；
3rd step: wire transporting, the energy consumption increased in the unit interval is EDW；
4th step: open illuminating lamp, the energy consumption increased in the unit interval is EDL；
5th step: carry out wire, the energy consumption increased in cutting to the unit interval before workpiece is EDKW；
6th step: straight cuts thickness is that the 4th step that compares in the unit interval during stainless steel work-piece of 1mm increases
Energy consumption is EDZCU；In the case of straight cuts, change cutting thickness is nmm, obtains the unit interval straight cuts energy under thickness nmm
Consumption EDnZCU, thickness coefficient ηHMeet EDnZCU=n*EDZCU*ηH, thus obtain the thickness coefficient η under nmm thicknessH；
7th step: straight cuts thickness be the Cutting Length in the unit interval during stainless steel work-piece of 1mm be X1D, straight line
Cutting thickness be the Cutting Length in the unit interval during stainless steel work-piece of nmm be XnD, thus obtain the speed under nmm thickness
Coefficient
8th step: circular arc cutting thickness is the stainless steel work-piece of 1mm, cut radius is the circle of 1mm, unit interval internal phase
The energy consumption that relatively the 4th step increases is EDYCU；In the case of circular arc cutting, the arc radius changing cutting is rmm, obtains radius rmm
Under unit interval circular arc cutting energy consumption EDrYCU, coefficient of curvature meets EDrYCU=r*EDYCU*ηQ, thus obtain under rmm radius
Coefficient of curvature ηQ；
9th step: straight cuts thickness is that the 4th step that compares in the unit interval during other materials workpiece of 1mm increases
Energy consumption be EDQZCU, the material coefficient of other materials is
As preferably, described velocity coeffficient ηS, material coefficient ηC, thickness coefficient ηH, coefficient of curvature ηQ, a be immeasurable
Guiding principle parameter, empty wire length XKW, straight cuts total length XZCU, circular arc cutting total length XYCU, straight cuts thickness be 1mm not
The Cutting Length in unit interval during rust steel workpiece is X1D, straight cuts thickness unit when being the stainless steel work-piece of nmm time
Interior Cutting Length XnDUnit be mm, empty wire travelling speed vw, cutting wire travelling speed vCUUnit be mm/s, empty wire
Energy consumption EKW, empty wire energy consumption E under unit length thicknessDKW, straight cuts energy consumption EZCU, straight cuts under unit length thickness
Energy consumption EDZCU, circular arc cutting energy consumption EYCU, circular arc cutting energy consumption E under unit length thicknessDYCU, host computer energy consumption ECO, unit
Time host computer energy consumption EDCO, wire transporting energy consumption EW, unit interval wire transporting energy consumption EDW, water pump bath energy consumption EP, unit interval water pump
Bath energy consumption EDP, illuminating lamp energy consumption EL, unit interval illuminating lamp energy consumption EDL, total energy consumption EZUnit be J.
Beneficial effects of the present invention: the present invention can pass through NC code, part material and thickness, it was predicted that goes out required total energy
Consumption.Basic database in the present invention is obtained by laboratory method.
Inventive feature and advantage will combine accompanying drawing by embodiment and be described in detail.
[accompanying drawing explanation]
Fig. 1 is the method figure of a kind of line based on NC code of present invention cutting energy consumption Forecasting Methodology；
Fig. 2 is the code analysis methods figure of a kind of line based on NC code of present invention cutting energy consumption Forecasting Methodology；
Fig. 3 is the 3B format program processing method figure of a kind of line based on NC code of present invention cutting energy consumption Forecasting Methodology；
Fig. 4 is the 3B format program circular arc processing method of a kind of line based on NC code of present invention cutting energy consumption Forecasting Methodology
Figure；
Fig. 5 is the ISO format program processing method figure of a kind of line based on NC code of present invention cutting energy consumption Forecasting Methodology；
Fig. 6 is the ISO format program circular arc process side of a kind of line based on NC code of present invention cutting energy consumption Forecasting Methodology
Method figure；
Fig. 7 is the basic database method for building up figure of a kind of line based on NC code of present invention cutting energy consumption Forecasting Methodology；
Fig. 8 is the embodiment cutting image shape graph of a kind of line based on NC code of present invention cutting energy consumption Forecasting Methodology.
[detailed description of the invention]
Refering to Fig. 1～Fig. 8, the present invention, comprise the following steps:
The first step: input preprocessing parts program code, part material and thickness；
Second step: code analysis, it is thus achieved that empty wire, straight cuts, circular arc cutting each several part parameter；
3rd step: call basic database, it is thus achieved that underlying parameter；
4th step: call energy consumption forecast model, it is thus achieved that prediction of energy consumption value: empty wire energy consumption EKWPass throughSolve, wherein EDKWFor the empty wire energy consumption under the unit time, XKWFor empty wire length, vwWalk for sky
Silk speed；Straight cuts energy consumption EZCUPass throughSolve, wherein EDZCUThick for the unit time
Straight cuts energy consumption under Du, XZCUFor straight cuts total length, vCUFor cutting wire travelling speed, ηSFor velocity coeffficient, ηHFor thickness
Coefficient, ηCFor material coefficient；Circular arc cutting energy consumption EYCUPass throughSolve, wherein
EDYCUEnergy consumption, X is cut for the circular arc under unit time thicknessYCUTotal length, η is cut for circular arcQFor coefficient of curvature；Host computer energy
Consumption ECOPass throughSolve, wherein EDCOFor unit time host computer energy
Consumption；Wire transporting energy consumption EWPass throughSolve, wherein EDWFor unit time wire transporting energy
Consumption；Water pump bath energy consumption EPPass throughSolve, wherein EDPFor unit time water pump
Bath energy consumption；Illuminating lamp energy consumption ELPass throughSolving, wherein a is 0 or 1, and a is 0
Representing and do not turn on light, a is that 1 expression is turned on light, EDLFor unit time illuminating lamp energy consumption；Total energy consumption EZPass through EZ=EKW+EZCU+EYCU+ECO+
EW+EP+ELSolve.
Concrete, the code analysis method in described second step is: judge whether have B in preprocessing parts program code,
Distinguish and use 3B format program processing method and ISO format program processing method, it is thus achieved that empty wire, straight cuts, circular arc cutting
Each several part parameter: empty wire length, straight cuts total length, circular arc cutting total length, circular arc cut radius.
Concrete, the basic database in described 3rd step includes material coefficient ηC, thickness coefficient ηH, velocity coeffficient ηS, empty
Wire travelling speed vw, cutting wire travelling speed vCU, empty wire energy consumption E under the unit intervalDKW, straight cuts under unit interval thickness
Energy consumption EDZCU, circular arc cutting energy consumption E under unit interval thicknessDYCU, coefficient of curvature ηQ, unit interval host computer energy consumption EDCO、
Unit interval wire transporting energy consumption EDW, unit interval water pump bath energy consumption EDP, unit interval illuminating lamp energy consumption EDL, empty wire travelling speed vw、
Cutting wire travelling speed vCUAcquisition according to lathe property, remaining parameter is all obtained by laboratory method, is obtained by data acquisition unit
Obtaining the power curve of wire cutting machine tool, the integration of power over time is energy consumption, and experimental procedure is:
The first step: open host computer, now measured energy consumption is E divided by the timeDCO；
Second step: opening water pump bath, the energy consumption increased in the unit interval is EDP；
3rd step: wire transporting, the energy consumption increased in the unit interval is EDW；
4th step: open illuminating lamp, the energy consumption increased in the unit interval is EDL；
5th step: carry out wire, the energy consumption increased in cutting to the unit interval before workpiece is EDKW；
6th step: straight cuts thickness is that the 4th step that compares in the unit interval during stainless steel work-piece of 1mm increases
Energy consumption is EDZCU；In the case of straight cuts, change cutting thickness is nmm, obtains the unit interval straight cuts energy under thickness nmm
Consumption EDnZCU, thickness coefficient ηHMeet EDnZCU=n*EDZCU*ηH, thus obtain the thickness coefficient η under nmm thicknessH；
7th step: straight cuts thickness be the Cutting Length in the unit interval during stainless steel work-piece of 1mm be X1D, straight line
Cutting thickness be the Cutting Length in the unit interval during stainless steel work-piece of nmm be XnD, thus obtain the speed under nmm thickness
Coefficient
8th step: circular arc cutting thickness is the stainless steel work-piece of 1mm, cut radius is the circle of 1mm, unit interval internal phase
The energy consumption that relatively the 4th step increases is EDYCU；In the case of circular arc cutting, the arc radius changing cutting is rmm, obtains radius rmm
Under unit interval circular arc cutting energy consumption EDrYCU, coefficient of curvature meets EDrYCU=r*EDYCU*ηQ, thus obtain under rmm radius
Coefficient of curvature ηQ；
9th step: straight cuts thickness is that the 4th step that compares in the unit interval during other materials workpiece of 1mm increases
Energy consumption be EDQZCU, the material coefficient of other materials is
Concrete, described velocity coeffficient ηS, material coefficient ηC, thickness coefficient ηH, coefficient of curvature ηQ, a be dimensionless
Parameter, empty wire length XKW, straight cuts total length XZCU, circular arc cutting total length XYCU, straight cuts thickness be the stainless of 1mm
The Cutting Length in unit interval during steel workpiece is X1D, straight cuts thickness unit interval when being the stainless steel work-piece of nmm
Interior Cutting Length XnDUnit be mm, empty wire travelling speed vw, cutting wire travelling speed vCUUnit be mm/s, empty wire energy
Consumption EKW, empty wire energy consumption E under unit length thicknessDKW, straight cuts energy consumption EZCU, straight cuts energy under unit length thickness
Consumption EDZCU, circular arc cutting energy consumption EYCU, circular arc cutting energy consumption E under unit length thicknessDYCU, host computer energy consumption ECO, unit time
Between host computer energy consumption EDCO, wire transporting energy consumption EW, unit interval wire transporting energy consumption EDW, water pump bath energy consumption EP, unit interval water pump punching
Water consumption EDP, illuminating lamp energy consumption EL, unit interval illuminating lamp energy consumption EDL, total energy consumption EZUnit be J.
Work process of the present invention:
A kind of line based on NC code of present invention cutting energy consumption Forecasting Methodology in the course of the work, is processed so that accompanying drawing 8 is required
Diagram shape thickness 3mm cast iron part as a example by illustrate, only know that need to process thickness is 3mm, material is cast iron, NC program
As follows, obtain empty wire, straight cuts, circular arc cutting each several part parameter in two kinds of situation, it is known that 3B format program and
Know ISO format program.One unit length, for 1mm；One unit thickness is 1mm；One unit interval is 1s.
Situation one:
3B format program is:
N001 B B10000 B10000 GY L4；
N002 B B25000 B50000 GX SR1；
N003 B3 B3 B30000 GX L3；
N004 B3 B3 B30000 GX L2；
N005 B10000 B24000 B42000 GX NR4；
N006 B70000 B B70000 GX L1；
N007 B B10000 B10000 GY L2；
N008 D。
The first step: identification N is alphabetic word, the number of N is n=8；
Second step: completion program, for each line program, if without space after identifying B, then supplements numeral 0；Identify number after B
Value is units, then identify the number of the 3rd B numerical value 0, add this number 0, become after above defect numerical value:
N001 B0 B10000 B10000 GY L4；
N002 B0 B25000 B50000 GX SR1；
N003 B30000 B30000 B30000 GX L3；
N004 B30000 B30000 B30000 GX L2；
N005 B10000 B24000 B42000 GX NR4；
N006 B70000 B0 B70000 GX L1；
N007 B0 B10000 B10000 GY L2；
N008 D；
3rd step: contrast the first row and the 7th line program, obtains empty wire length XKWFor 15mm, straight cuts a length of two
After the 3rd B of row, numerical value deducts 15mm after removing 1000 and is 5mm；
4th step: the program line of band L is carried out straight cuts parameter extraction operation, extracting method is: the number after first B
According to note B1, the data note B2 after second B, the radical sign value of a length of B1, B2 quadratic sum, unit um；3rd line program, the 4th stroke
Sequence, the 6th line program carry out length extractionExtract
To parameter add up as straight cuts total length
5th step: the program line of band R is carried out the extraction of circular arc cutting parameter and operates, the data accumulation that will extract successively,
The program of final stage band R until having added up, thus obtain circular arc cutting total length XYCU: the first step: obtain starting point place as
Limit, cut direction: the 2nd line program is band SR1, and starting point is at first quartile, and cut direction is clockwise；5th line program is band NR4,
Starting point is in fourth quadrant, and cut direction is counterclockwise；Second step: solve this quadrant long at this counting direction length value, circular arc, geometry
Method solves: the first quartile inside counting direction length value method of geometry of the 2nd line program solves and obtains 25000um, and circular arc is long
It is 12500 π um, r=25000um；The fourth quadrant inside counting direction length value method of geometry of the 5th line program solves and obtains
16000um, circular arc is a length ofR=26000um；3rd step: judge counting step and second step institute
Seek length relation, and solve: the counting step of the 2nd line program is 50000um, 50000um 25000um, then length in this quadrant
Length 12500 π um required by second step, remaining quadrant is obtained by following operation；The counting step of the 5th line program is
42000um, 42000um > 16000um, then length required by a length of second step in this quadrantIts
Afterimage limit is obtained by following operation；4th step: judge the quadrant district of process: the 2nd line program, (counting step 50000um-second
Required length 25000um of step)/radius 25000um obtains more than 10, then judges to cut 14 liang clockwise according to 14321432 sequence numbers
Individual quadrant；5th line program, (length 16000um required by counting step 42000um-second step)/radius 26000um obtains more than 10, then
Judge to cut 41 two quadrants counterclockwise according to 12341234 sequence numbers；5th step: a point quadrant solves, and the whole quadrant that comprises is
/ 4th of full circle, the most whole quadrant comprised utilizes geometrical relationship to solve: the fourth quadrant Inner arc of the 2nd line program is a length of
12500 π um, a length of 13000 π um of first quartile Inner arc of the 5th line program；6th step: accumulative all quadrants arc length obtains
Counting step is more than circular arc Cutting Length during length required by second step, and the 2nd line program is 25000 π um=25 π mm, the 5th stroke
Sequence is
Add up the parameter extracted to cut total length for circular arc
Situation two:
ISO format program is:
N001 G00 X0 Y0；
N002 G01 Y-10000；
N003 G02 Y-60000 J-35000；
N004 G01 X-30000 Y-90000；
N005 G01 X-60000 Y-60000；
N006 G03 X-70000 Y-10000 I-70000 J-36000；
N007 G01 Y-10000；
N008 G00 X0 Y0；
N009 M02。
The first step: identification N is alphabetic word, the number of N be n, n be 9；
Second step: completion program, for each line program, identifies that one only occur in X, Y, then supplement another letter and at this
Numeral 0 is supplemented after letter；Identify that one only occur in I, J, then supplement another letter and after this letter, supplement numeral 0, becoming:
N001 G00 X0 Y0；
N002 G01 X0 Y-10000；
N003 G02 X0 Y-60000 I0 J-35000；
N004 G01 X-30000 Y-90000；
N005 G01 X-60000 Y-60000；
N006 G03 X-70000 Y-10000 I-70000 J-36000；
N007 G01 X0 Y-10000；
N008 G00 X0 Y0；
N009 M02；
3rd step: the 3rd step: relatively front two row and last two line programs, obtains empty wire length XKWFor 15mm, straight cuts
A length of: the radical sign value of Xn-2, Yn-2 quadratic sum of the radical sign value of second segment running length X2, Y2 quadratic sum+last three row-
15000, unit is um,
4th step: will carry out straight cuts parameter extraction operation except the program line of front two row and the band G01 of last three row, extracts
Method is: assume to extract the i-th step data, and (Xi-Xi-1), the radical sign value of (Yi-Yi-1) quadratic sum are the length of this step operation, unit
For um；The data accumulation extracted is straight cuts total length XZCU；4th line program, the 5th line program, the 7th line program are except front two
The program line of the band G01 of row and last three row, the straight cuts obtaining three row according to said method is a length of:
70mm, thus add up to obtain straight cuts total length
5th step: the program line recognizing G02, G03 carrying out circular arc cutting parameter and extracts operation, the data extracted are tired out
It is calculated as circular arc cutting total length XYCU: the first step: obtain starting point coordinate, central coordinate of circle, terminal point coordinate, radius；Starting point coordinate extracts
Numerical value composition after X, Y of lastrow program, central coordinate of circle is numerical value composition after I, J of this line program, and terminal point coordinate is this
Numerical value composition after X, Y of line program, radius is the distance of starting point and the center of circle, the 3rd line program, starting point coordinate (0 ,-10000),
Central coordinate of circle (0 ,-35000), terminal point coordinate (0 ,-60000), radius 25000；6th line program, starting point coordinate (-60000 ,-
60000), central coordinate of circle (-70000 ,-36000), terminal point coordinate (-70000 ,-10000), radius 26000；Second step: Gui Yuan
Change: being converted into zero point coordinate, starting point coordinate, central coordinate of circle, the transverse and longitudinal coordinate of terminal point coordinate all deduct the numerical value of central coordinate of circle, the
3 line programs, starting point coordinate (0,25000), central coordinate of circle (0,0), terminal point coordinate (0 ,-25000)；6th line program, starting point coordinate
(10000 ,-24000), central coordinate of circle (0,0), terminal point coordinate (0,26000)；3rd step: judge starting point coordinate, terminal point coordinate institute
At quadrant, it is judged that method is: abscissa, vertical coordinate are just for first quartile, abscissa for negative, vertical coordinate for be just second as
Limit, abscissa, vertical coordinate are to be born as third quadrant, abscissa for just, vertical coordinate for bearing as fourth quadrant, the 3rd line program,
Point is at first quartile, and terminal is in fourth quadrant；6th line program, starting point is in fourth quadrant, and terminal is at first quartile；4th step: sentence
The quadrant district that cracked ends is crossed, the program segment of band G02 judges according to 14321432 sequence numbers, and terminus place quadrant numeral goes out in order
Existing word string, it is desirable to less than four numerals；Program segment with G03 judges according to 12341234 sequence numbers, terminus place as
The word string that limit numeral occurs in order, it is desirable to less than four numerals；3rd line program, crosses first quartile, fourth quadrant；6th row
Program, crosses fourth quadrant, first quartile；5th step: a point quadrant solves, whole comprise that quadrant is full circle 1/4th, the most whole
The individual quadrant comprised utilizes geometrical relationship to solve, the 3rd line program, and first quartile is 12500 π um, and fourth quadrant is 12500 π
um；6th line program, fourth quadrant isFirst quartile is 13000 π um；6th step: accumulative all quadrants circle
Arc length degree, the 3rd line program is 25000 π um=25 π mm；6th line program is
Add up the parameter extracted to cut total length for circular arc
After obtaining empty wire, straight cuts, circular arc cutting each several part parameter according to 3B program or ISO program, can call
Basic database, it is thus achieved that underlying parameter, underlying parameter includes sky wire travelling speed vw, cutting wire travelling speed vCU, velocity coeffficient ηS, thick
Degree coefficient ηH, material coefficient ηC, coefficient of curvature ηQ.Above-mentioned basic database is except empty wire travelling speed vwWith cutting wire travelling speed vCUBy
Outside lathe property obtains, remaining is built by laboratory method.Call energy consumption forecast model, prediction of energy consumption value can be obtained.
The present invention, it is possible to by NC code, part material and thickness, it was predicted that go out required total energy consumption.Basis in the present invention
Data base is obtained by laboratory method.
Above-described embodiment is the description of the invention, is not limitation of the invention, any to simple transformation of the present invention after
Scheme belong to protection scope of the present invention.
Claims (4)
1. line based on a NC code cutting energy consumption Forecasting Methodology, it is characterised in that: comprise the following steps:
The first step: input preprocessing parts program code, part material and thickness；
Second step: code analysis, it is thus achieved that empty wire, straight cuts, circular arc cutting each several part parameter；
3rd step: call basic database, it is thus achieved that underlying parameter；
4th step: call energy consumption forecast model, it is thus achieved that prediction of energy consumption value: empty wire energy consumption EKWPass throughAsk
Solve, wherein EDKWFor the empty wire energy consumption under the unit time, XKWFor empty wire length, vwFor empty wire travelling speed；Straight cuts energy consumption
EZCUPass throughSolve, wherein EDZCUFor the straight cuts energy consumption under unit time thickness,
XZCUFor straight cuts total length, vCUFor cutting wire travelling speed, ηSFor velocity coeffficient, ηHFor thickness coefficient, ηCFor material coefficient；Circle
Arc cutting energy consumption EYCUPass throughSolve, wherein EDYCUFor under unit time thickness
Circular arc cutting energy consumption, XYCUTotal length, η is cut for circular arcQFor coefficient of curvature；Host computer energy consumption ECOPass throughSolve, wherein EDCOFor unit time host computer energy consumption；Wire transporting energy consumption
EWPass throughSolve, wherein EDWFor unit time wire transporting energy consumption；Water pump is washed by water
Energy consumption EPPass throughSolve, wherein EDPFor unit time water pump bath energy consumption；According to
Bright lamp energy consumption ELPass throughSolving, wherein a is 0 or 1, and a is that 0 expression is not turned on light,
A is that 1 expression is turned on light, EDLFor unit time illuminating lamp energy consumption；Total energy consumption EZPass through EZ=EKW+EZCU+EYCU+ECO+EW+EP+ELAsk
Solve.
A kind of line based on NC code cutting energy consumption Forecasting Methodology, it is characterised in that: described second
Code analysis method in step is: judge whether to have in preprocessing parts program code B, distinguishes and uses 3B format program process side
Method and ISO format program processing method, it is thus achieved that empty wire, straight cuts, circular arc cutting each several part parameter: empty wire length, straight
Line cutting total length, circular arc cutting total length, circular arc cut radius.
A kind of line based on NC code cutting energy consumption Forecasting Methodology, it is characterised in that: the described 3rd
Basic database in step includes material coefficient ηC, thickness coefficient ηH, velocity coeffficient ηS, empty wire travelling speed vw, cutting wire travelling speed
vCU, empty wire energy consumption E under the unit intervalDKW, straight cuts energy consumption E under unit interval thicknessDZCU, under unit interval thickness
Circular arc cutting energy consumption EDYCU, coefficient of curvature ηQ, unit interval host computer energy consumption EDCO, unit interval wire transporting energy consumption EDW, unit
Time water pump bath energy consumption EDP, unit interval illuminating lamp energy consumption EDL, empty wire travelling speed vw, cutting wire travelling speed vCUAcquisition root
According to lathe property, remaining parameter is all obtained by laboratory method, is obtained the power curve of wire cutting machine tool by data acquisition unit,
The integration of power over time is energy consumption, and experimental procedure is:
The first step: open host computer, now measured energy consumption is E divided by the timeDCO；
Second step: opening water pump bath, the energy consumption increased in the unit interval is EDP；
3rd step: wire transporting, the energy consumption increased in the unit interval is EDW；
4th step: open illuminating lamp, the energy consumption increased in the unit interval is EDL；
5th step: carry out wire, the energy consumption increased in cutting to the unit interval before workpiece is EDKW；
6th step: straight cuts thickness is the energy consumption that the 4th step increases that compares in the unit interval during stainless steel work-piece of 1mm
For EDZCU；In the case of straight cuts, change cutting thickness is nmm, obtains the unit interval straight cuts energy consumption under thickness nmm
EDnZCU, thickness coefficient ηHMeet EDnZCU=n*EDZCU*ηH, thus obtain the thickness coefficient η under nmm thicknessH；
7th step: straight cuts thickness be the Cutting Length in the unit interval during stainless steel work-piece of 1mm be X1D, straight cuts
Thickness be the Cutting Length in the unit interval during stainless steel work-piece of nmm be XnD, thus obtain the velocity coeffficient under nmm thickness
8th step: circular arc cutting thickness is the stainless steel work-piece of 1mm, cut radius is the circle of 1mm, compares in the unit interval
The energy consumption that 4th step increases is EDYCU；In the case of circular arc cutting, the arc radius changing cutting is rmm, obtains under radius rmm
Unit interval circular arc cutting energy consumption EDrYCU, coefficient of curvature meets EDrYCU=r*EDYCU*ηQ, thus obtain the curvature under rmm radius
Coefficient ηQ；
9th step: straight cuts thickness is the energy that the 4th step increases that compares in the unit interval during other materials workpiece of 1mm
Consumption is EDQZCU, the material coefficient of other materials is
A kind of line based on NC code cutting energy consumption Forecasting Methodology, its feature exists
In: described velocity coeffficient ηS, material coefficient ηC, thickness coefficient ηH, coefficient of curvature ηQ, a be dimensionless group, empty wire
Length XKW, straight cuts total length XZCU, circular arc cutting total length XYCU, straight cuts thickness is when being the stainless steel work-piece of 1mm
Cutting Length in unit interval is X1D, cutting in straight cuts thickness unit interval when being the stainless steel work-piece of nmm long
Degree XnDUnit be mm, empty wire travelling speed vw, cutting wire travelling speed vCUUnit be mm/s, empty wire energy consumption EKW, unit
Empty wire energy consumption E under length thicknessDKW, straight cuts energy consumption EZCU, straight cuts energy consumption E under unit length thicknessDZCU, circle
Arc cutting energy consumption EYCU, circular arc cutting energy consumption E under unit length thicknessDYCU, host computer energy consumption ECO, unit interval computer master
Function consumption EDCO, wire transporting energy consumption EW, unit interval wire transporting energy consumption EDW, water pump bath energy consumption EP, unit interval water pump bath energy consumption
EDP, illuminating lamp energy consumption EL, unit interval illuminating lamp energy consumption EDL, total energy consumption EZUnit be J.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610328991.5A CN106021179B (en) | 2016-05-18 | 2016-05-18 | A kind of wire cutting energy consumption prediction technique based on NC code |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610328991.5A CN106021179B (en) | 2016-05-18 | 2016-05-18 | A kind of wire cutting energy consumption prediction technique based on NC code |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106021179A true CN106021179A (en) | 2016-10-12 |
| CN106021179B CN106021179B (en) | 2019-01-22 |
Family
ID=57098588
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610328991.5A Expired - Fee Related CN106021179B (en) | 2016-05-18 | 2016-05-18 | A kind of wire cutting energy consumption prediction technique based on NC code |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106021179B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106475908A (en) * | 2016-11-08 | 2017-03-08 | 上海大学 | Follow grinding process lathe operation energy consumption Forecasting Methodology based on standard G code |
| CN108255132A (en) * | 2016-12-28 | 2018-07-06 | 浙江科技学院 | A kind of method for waveform identification based on wire cutting power curve mutation waveform database |
| CN109613890A (en) * | 2019-01-10 | 2019-04-12 | 湖南科技大学 | A kind of machine tooling efficiency prediction technique |
| CN110560921A (en) * | 2019-08-22 | 2019-12-13 | 浙江科技学院 | total energy consumption prediction method for laser cutting based on shortest distance |
| CN110560922A (en) * | 2019-08-22 | 2019-12-13 | 浙江科技学院 | Machining energy consumption prediction method for laser cutting based on geometric features |
| CN112276263A (en) * | 2020-10-14 | 2021-01-29 | 宁波市博虹机械制造开发有限公司 | G code-based special motion control method for electric spark forming machine |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103390195A (en) * | 2013-05-28 | 2013-11-13 | 重庆大学 | Machine workshop task scheduling energy-saving optimization system based on reinforcement learning |
| CN104834787A (en) * | 2015-05-18 | 2015-08-12 | 江南大学 | Perception and analysis method for energy consumption of bearing grinding machine production workshop |
| CN105159237A (en) * | 2015-06-29 | 2015-12-16 | 中国西电电气股份有限公司 | Energy consumption prediction method for digitalized workshop numerical control machine tool |
| CN105204433A (en) * | 2015-09-30 | 2015-12-30 | 重庆大学 | Machine tool energy-saving running method for switching machine tool states within random machining interval |
| CN105259791A (en) * | 2015-11-16 | 2016-01-20 | 哈尔滨工业大学 | Machining parameter optimization method based on general cutting energy consumption model |
-
2016
- 2016-05-18 CN CN201610328991.5A patent/CN106021179B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103390195A (en) * | 2013-05-28 | 2013-11-13 | 重庆大学 | Machine workshop task scheduling energy-saving optimization system based on reinforcement learning |
| CN104834787A (en) * | 2015-05-18 | 2015-08-12 | 江南大学 | Perception and analysis method for energy consumption of bearing grinding machine production workshop |
| CN105159237A (en) * | 2015-06-29 | 2015-12-16 | 中国西电电气股份有限公司 | Energy consumption prediction method for digitalized workshop numerical control machine tool |
| CN105204433A (en) * | 2015-09-30 | 2015-12-30 | 重庆大学 | Machine tool energy-saving running method for switching machine tool states within random machining interval |
| CN105259791A (en) * | 2015-11-16 | 2016-01-20 | 哈尔滨工业大学 | Machining parameter optimization method based on general cutting energy consumption model |
Non-Patent Citations (4)
| Title |
|---|
| 丁春伟等: "智能数控提高机床能效的研究", 《机电产品开发与创新》 * |
| 何彦等: "数控机床多能量源的动态能耗建模与仿真方法", 《机械工程学报》 * |
| 徐立云等: "机加工过程一类能耗估算模型", 《同济大学学报（自然科学版）》 * |
| 郑军等: "基于 Petri 网的砂型铸造工序能耗计算模型", 《中国机械工程》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106475908A (en) * | 2016-11-08 | 2017-03-08 | 上海大学 | Follow grinding process lathe operation energy consumption Forecasting Methodology based on standard G code |
| CN108255132A (en) * | 2016-12-28 | 2018-07-06 | 浙江科技学院 | A kind of method for waveform identification based on wire cutting power curve mutation waveform database |
| CN108255132B (en) * | 2016-12-28 | 2021-01-08 | 浙江科技学院 | Waveform identification method based on linear cutting power curve mutation waveform database |
| CN109613890A (en) * | 2019-01-10 | 2019-04-12 | 湖南科技大学 | A kind of machine tooling efficiency prediction technique |
| CN110560921A (en) * | 2019-08-22 | 2019-12-13 | 浙江科技学院 | total energy consumption prediction method for laser cutting based on shortest distance |
| CN110560922A (en) * | 2019-08-22 | 2019-12-13 | 浙江科技学院 | Machining energy consumption prediction method for laser cutting based on geometric features |
| CN112276263A (en) * | 2020-10-14 | 2021-01-29 | 宁波市博虹机械制造开发有限公司 | G code-based special motion control method for electric spark forming machine |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106021179B (en) | 2019-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106021179A (en) | NC code-based linear cutting energy consumption prediction method | |
| CN106001811B (en) | A kind of wire cutting energy consumption Forecasting Methodology based on geometric properties | |
| CN109657039B (en) | Work history information extraction method based on double-layer BilSTM-CRF | |
| CN105045927B (en) | Construction project labor and materials machine data automatic coding and system | |
| CN103592214A (en) | Method for evaluating corrosion state of carbon steel material in atmospheric environment | |
| CN112547807B (en) | Decision tree algorithm-based hot rolled strip quality accurate judgment method | |
| CN107807910A (en) | A kind of part-of-speech tagging method based on HMM | |
| CN107544450B (en) | Process industry network model construction method and system based on data | |
| CN103971019A (en) | Method for predicting workpiece machining energy consumption based on geometrical characteristics | |
| CN110287456A (en) | Bulk lots volume rolled surface defect analysis method based on data mining | |
| CN106020430B (en) | A kind of basic database method for building up of wire cutting energy consumption prediction | |
| CN109102538A (en) | Method, the synoptic analysis method of piston ring land characteristic point and land features line are extracted using isopleth data | |
| CN108021878B (en) | Contour line negative-direction landform intelligent identification method | |
| CN106529741A (en) | Space relevant characteristic-based ultra-short-period wind power prediction method | |
| Shi et al. | An improved faster R-CNN for steel surface defect detection | |
| CN109239809A (en) | A kind of weather piston ring land line drawing method based on chord contact method | |
| CN109472038A (en) | A kind of matched method and system of smelting iron and steel technological parameter tracing pattern | |
| CN106127340A (en) | A kind of computational methods being applicable to ultra-short term | |
| CN112906215A (en) | Pipe tool path generation method based on SolidWorks secondary development | |
| Hanjie et al. | EDM machining quality control and parameters optimization | |
| CN110717044A (en) | Text classification method for research and report text | |
| CN103451707A (en) | Method for coating nano-ceramic coating on surface of precision welded pipe | |
| Gilbert et al. | Visualisation of probabilistic access forecasts for offshore operations | |
| CN103901819B (en) | The process operation fast searching method of feature based mark | |
| CN104021312A (en) | Flood process similarity analytic method and system based on process linear centroid |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190122 Termination date: 20200518 |