CN102179725B - Arrangement method of heat characteristic monitoring measurement points of numerical control machine - Google Patents
Arrangement method of heat characteristic monitoring measurement points of numerical control machine Download PDFInfo
- Publication number
- CN102179725B CN102179725B CN2011100653659A CN201110065365A CN102179725B CN 102179725 B CN102179725 B CN 102179725B CN 2011100653659 A CN2011100653659 A CN 2011100653659A CN 201110065365 A CN201110065365 A CN 201110065365A CN 102179725 B CN102179725 B CN 102179725B
- Authority
- CN
- China
- Prior art keywords
- measuring point
- temperature
- control machine
- numerical control
- sequence
- 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.)
- Expired - Fee Related
Links
- 238000012544 monitoring process Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000005259 measurement Methods 0.000 title abstract 13
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 230000035945 sensitivity Effects 0.000 claims description 27
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 238000012216 screening Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 abstract description 5
- 230000008901 benefit Effects 0.000 abstract description 2
- 206010057040 Temperature intolerance Diseases 0.000 abstract 2
- 230000008543 heat sensitivity Effects 0.000 abstract 2
- 238000003801 milling Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 238000003754 machining Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000005144 thermotropism Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000020169 heat generation Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 230000007935 neutral effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
Images
Landscapes
- Automatic Control Of Machine Tools (AREA)
Abstract
The invention discloses an arrangement method of heat characteristic monitoring measurement points of a numerical control machine. The method comprises the following steps of: (1) selecting measurement points on a numerical control machine body, and acquiring temperature values of the measurement points under different time and deformation of a tool nose position relative to a workpiece under different time according to the frequency to form a temperature change sequence X of the measurement points and a deformation sequence delta of the tool nose position relative to the workpiece; (2) calculating temperature correlation coefficients according to the temperature change sequence X of the measurement points, and dividing the measurement points into groups; (3) performing heat sensitivity calculation on the measurement points in each group according to the data of the temperature change sequence X of the measurement points and the deformation sequence delta of the tool nose position relative to the workpiece, and selecting the measurement points with highest heat sensitivity in each group to form a selection group; and (4) selecting the measurement points serving as the heat characteristic monitoring measurement points of the numerical control machine from the selection group according to different arrangement number of the heat characteristic monitoring measurement points of the numerical control machine. The method disclosed by the invention has the advantages of simplicity, effectiveness, clear principle, simplicity in operation, good controllability and the like.
Description
Technical field
The present invention relates to Digit Control Machine Tool thermal characteristic monitoring field, especially relate to the method for arranging of measuring point in a kind of Digit Control Machine Tool thermal characteristic monitoring.
Background technology
Machining accuracy of NC machine tool is to weigh the important indicator of Digit Control Machine Tool service behaviour; A large amount of researchs show that hot error is the worst error source of precision processing machine such as Digit Control Machine Tool, accounts for about 70% of overall error; And for ultra-precision machine tool, even reach nearly 90%.Because Digit Control Machine Tool exists inside and outside thermal source, particularly internal heat resource inevitably in process, must cause the generation of thermal deformation errors.Thereby the Digit Control Machine Tool thermal characteristic monitored, control that the thermal deformation influence just seems extremely important in its process.
Measure inconvenience in real time owing to carry out thermal deformation during to machine tooling; In reality,, combine the model of temperature and thermal deformation to extrapolate lathe in Current Temperatures distortion situation after the match through the variations in temperature of analyzing lathe usually through some some temperature sensors of thermo-responsive location arrangements at lathe.At present be applied to set up the model that concerns between temperature and the thermal deformation more relevant report is arranged as various modeling algorithms such as linear regression, SVMs, neutral nets.
Consider complexity, the diversity of Digit Control Machine Tool thermal source, in the actual processing of lathe, start from cost consideration simultaneously, the arrangement requirement of temperature sensor measuring point was during the machine tool thermal step response was monitored: (1) sensor will lack as much as possible; (2) to confirm that the sensor of being arranged can reflect as far as possible accurately that lathe overall thermal step response changes; (3) in order to give full play to the performance of sensor, the temperature data that each sensor records is will the temperature data degree of coupling independent as far as possible and that other sensors record little; (4) in order to improve the accuracy of subsequent thermal The deformation calculation to greatest extent, it is relatively more responsive with respect to thermal deformation to require sensor to record the position temperature, to reduce the influence of measure error to thermal deformation.
Consider this these requirements, reasonably select the measuring point of the hot monitoring of lathe just to seem extremely important.At present, mainly still depend on experience, the selection of measuring point is also lacked guide of theory and efficient ways, exist technical problem for the selection of Digit Control Machine Tool thermal characteristic monitoring stations.
Summary of the invention
The invention provides a kind of method for arranging of numerical control machine heat characteristic monitoring stations; Behind the temperature field and thermal deformation data of measuring Digit Control Machine Tool; Calculating and the temperature-sensitive sensitivity of using data processing model to carry out the temperature correlation property coefficient of measuring point are calculated; Filter out the measuring point of optimum number controlled machine thermal characteristics monitoring, this method have simple effectively, clear principle, simple to operate, characteristics such as controllability is good.
A kind of method for arranging of numerical control machine heat characteristic monitoring stations may further comprise the steps:
1) on the Digit Control Machine Tool body, chooses the measuring point that is used for the detected temperatures variation; Press the temperature on the frequency collection measuring point; Obtain the temperature value of measuring point under the different time, described temperature value in chronological sequence is combined to form the variations in temperature sequence X of measuring point, in the temperature on gathering measuring point; With respect to the deflection of workpiece, described deflection in chronological sequence is combined to form the deflection sequence Δ of position of tool tip with respect to workpiece according to position of tool tip under the described frequency collection different time;
2) carry out after the temperature dependency coefficient calculations measuring point being divided into each group according to the variations in temperature sequence X of measuring point in the step 1);
3) according to the variations in temperature sequence X of measuring point in the step 1) and position of tool tip with respect to the deflection sequence Δ of workpiece to step 2) in each measuring point in each group carry out the temperature-sensitive sensitivity and calculate; Through relatively after the screening measuring point of the maximum heat susceptibility in every group being selected, form the selection group;
4) if numerical control machine heat characteristic monitoring stations arranges that number equates with measuring point number in the selection group; All measuring points in the then selection group are promptly as numerical control machine heat characteristic monitoring stations; If numerical control machine heat characteristic monitoring stations arranges that number is less than the measuring point number in the selection group; From the selection group, press the descending select progressively and the measuring point of arranging the number equal amount of temperature-sensitive sensitivity of measuring point, the measuring point of selection is promptly as numerical control machine heat characteristic monitoring stations.
In order to obtain better invention effect, the present invention is carried out preferably:
The described measuring point that is used for the detected temperatures variation is distributed near lathe bed, the main shaft and the leading screw vicinity.Consider that main shaft and leading screw are the thermal source of lathe; Distortion is the main cause of thermal deformation of machine tool and lathe bed is heated afterwards; For the measuring point that makes selection can better react the whole thermal deformation situation of lathe, all arranged a series of temperature point at lathe bed, spindle jacket and lead screw shaft bearing sleeve place, preferred; The described number that is used for the measuring point that detected temperatures changes is more than or equal to 10, and is uniformly distributed in lathe bed, spindle jacket and lead screw shaft bearing sleeve place.The numerical control machine heat characteristic monitoring stations of after temperature dependency coefficient calculations and maximum heat susceptibility calculate, selecting like this can better be monitored the thermal characteristics of Digit Control Machine Tool, helps follow-up thermal deformation to Digit Control Machine Tool and carries out modeling and error compensation.
Temperature on the described measuring point adopts temperature sensor to detect, and the TEMP implement body is optional to be selected often, and preferred temperature sensor can be selected thermocouple and thermal resistance for use.
Described position of tool tip adopts displacement transducer to detect with respect to the deflection of workpiece; Displacement transducer can be selected laser displacement sensor for use, and position of tool tip is with respect to the deflection of workpiece, the hot error of promptly being brought by thermal deformation; Hot source of error in Z axle, X axle and three directions of Y axle (be Z to, X to Y to) thermal deformation that brings; Z accounts for the largest percentage in hot error to thermal deformation (thermal stretching of main shaft), and generally for easy to detect and from cost consideration, described position of tool tip is preferably the deflection of position of tool tip with respect to the workpiece Z-direction with respect to the deflection of workpiece; Can use a laser displacement sensor to the Digit Control Machine Tool position of tool tip with respect to workpiece Z to deflection detect; The situation that not only can reflect numerical control machining tool heat error has also been practiced thrift use cost to a certain extent, can certainly use three laser displacement sensors that the deflection of Z axle, X axle, three directions of Y axle is detected; Average then; Can better detect the hot error of position of tool tip, can guarantee that the result that the inventive method obtains has better accuracy, use two laser displacement sensors to detect the deflection of both direction with respect to workpiece; And then obtain the deflection of position of tool tip with respect to workpiece, also can meet the demands.
For convenience that detects and the processing of being convenient to data, described frequency is preferably fixed frequency, under fixed frequency; Be the same the blanking time of the collection of the temperature on the measuring point, can guarantee the stability of Data Detection on the one hand like this, on the other hand; The data that detect are the time periods that are dispersed in detection; In the calculating of carrying out temperature correlation property coefficient and temperature-sensitive sensitivity, better accuracy rate can be arranged, thereby guarantee the accuracy that numerical control machine heat characteristic monitoring stations is chosen.Frequency has been represented the blanking time of data acquisition; In the inventive method; So long as the blanking time of temperature acquisition and position of tool tip are equal just with respect to the blanking time of the deflection of workpiece; Just can reach the corresponding relation that has that makes corresponding temperature and thermal deformation, the temperature that promptly some time inscribes on the measuring point is to the deflection of the position of tool tip under should be constantly with respect to workpiece.Therefore, frequency be not necessarily need fixing, can frequency conversion, but the frequency that the frequency that will guarantee the measuring point temperature acquisition and position of tool tip are gathered with respect to the deflection of workpiece is equal, so just can be to satisfy the object of the invention.
Described step 2) may further comprise the steps:
(a) in measuring point, appoint and get a measuring point, note is made i measuring point, is benchmark with i measuring point, gathers N time altogether by frequency, and then the variations in temperature sequence X of i measuring point is expressed as X
i={ x
I1, x
I2... X
IN, will carry out the temperature dependency coefficient calculations except that i measuring point and i measuring point the measuring point, as if being k measuring point except that i the measuring point the measuring point, the variations in temperature sequence X of k measuring point is expressed as X
k={ x
K1, x
K2... X
KN, calculate the temperature correlation property coefficient between i measuring point and k the measuring point, the variations in temperature sequence X
i={ x
I1, x
I2... X
INIn maximum be x
Imax, the variations in temperature sequence X
k={ x
K1, x
K2... X
KNIn maximum be x
Kmax
(b) the temperature dependency coefficients by using between i measuring point and k the measuring point
Formula calculates, wherein,
Obtain C
IkValue;
(c) C
Ik>=0.7, represent that these two measuring points have obvious correlation, are placed on one group, C with k measuring point and i measuring point
Ik＜0.7 these two measuring points of expression do not have obvious correlation; K measuring point and i measuring point are not placed on one group; Then to carrying out the temperature dependency coefficient calculations except that k measuring point and i other measuring points and i measuring point the measuring point; Obtain by with measuring point and i the group that measuring point is formed of i measuring point temperature correlation property coefficient>=0.7, the measuring point that divides into groups that fails is formed remaining set;
(d) appoint from the remaining set of step (c) that to get a point be benchmark, the repeating step (a) and (b) with (c), obtain a new group again, finish grouping up to all measuring points.
In the step 3), described temperature-sensitive sensitivity is calculated, and may further comprise the steps:
(e) in measuring point, appoint and get a measuring point, note is made h measuring point, gathers N time altogether by frequency, and then the variations in temperature sequence X of h measuring point is expressed as X
h={ x
H1, x
H2... X
HN, position of tool tip is represented Δ={ δ with respect to the deflection sequence Δ of workpiece
1, δ
2... δ
N, the variations in temperature sequence X
h={ x
H1, x
H2... X
HNIn maximum be x
Hmax, position of tool tip is with respect to the deflection sequence Δ={ δ of workpiece
1, δ
2... δ
NIn maximum be δ
Max
(f) the temperature-sensitive sensitivity of h measuring point is defined as
obtains the temperature-sensitive sensitivity of h measuring point after calculating;
(g) repeating step (e) and (f) calculates carry out the temperature-sensitive sensitivity except that h other measuring points the measuring point.
In the step 4); Numerical control machine heat characteristic monitoring stations layout number can be taken all factors into consideration according to the factors such as accuracy of cost, detection; Generally speaking, adjust according to actual conditions, numerical control machine heat characteristic monitoring stations arranges that number equates with measuring point number in the selection group; Then such numerical control machine heat characteristic monitoring stations arranges that carrying out the follow-up hot error-detecting that causes for thermal deformation has accuracy preferably; Consider the cost reason, can select numerical control machine heat characteristic monitoring stations layout number, from the selection group, press the descending select progressively measuring point of temperature-sensitive sensitivity of measuring point less than the measuring point number in the selection group; The accuracy of the hot error that can guarantee to greatest extent like this to detect, thus help follow-up thermal deformation of machine tool being carried out modeling and error compensation.
Compared with prior art, the useful effect that has of the present invention is:
Numerical control machine heat characteristic monitoring stations method for arranging of the present invention has solved present Digit Control Machine Tool thermal characteristic monitoring stations and has arranged main dependence experience, lacks problems such as rationale and effective method.
Numerical control machine heat characteristic monitoring stations method for arranging of the present invention have simple effectively, clear principle, easy to operate, advantage such as controllability is good, can reduce Digit Control Machine Tool thermal characteristic monitoring stations number effectively, avoid temperature strong correlation between measuring point.
The inventive method obtains numerical control machine heat characteristic monitoring stations, can reflect the thermal deformation that Digit Control Machine Tool is whole more comprehensively, exactly, helps follow-up thermal deformation to Digit Control Machine Tool and carries out modeling and error compensation, thereby guaranteed the accuracy of error compensation.
Description of drawings
Fig. 1 is the flow chart of Digit Control Machine Tool thermal characteristics monitoring stations method for arranging of the present invention;
Fig. 2 is the sketch map of point position on the Digit Control Machine Tool in the embodiment of the invention 1, and wherein, T1～T15 is 15 measuring points that embodiment 1 is arranged in Digit Control Machine Tool;
Fig. 3 is the temperature variation curve that T2, T7, T9, T11 and T12 are ordered in the embodiment of the invention 1;
Fig. 4 be the embodiment of the invention 1 position of tool tip with respect to workpiece Z to the change curve of deflection.
The specific embodiment
Embodiment 1
Below in conjunction with accompanying drawing and embodiment the present invention is described further.
As shown in Figure 2; Machine adopted CNC milling machine XK713, the three-dimensional stroke of its X-direction, Y direction and Z-direction is 600 * 400 * 500mm, near lathe bed, main shaft, the leading screw vicinity all arranged a series of temperature point; On the CNC milling machine machine body, select and be used for different measuring point T1～T15 that detected temperatures changes; T1～T5 is on the column of CNC milling machine lathe bed, and T6～T8 is on the base of CNC milling machine lathe bed, and T9～T10 is on the lead screw shaft bearing of CNC milling machine; T11 is at the main shaft bottom end cover of CNC milling machine; T12～T15 and then installs 15 thermal resistance PT100 temperature sensors at measuring point T1～T15 on the rack sleeve of CNC milling machine, be used for collecting temperature and change.
With a laser displacement sensor (LK-G5000) measure position of tool tip with respect to workpiece Z to deflection, consider that the thermal stretching (Z thermotropism error) of main shaft is in thermal deformation, to account for the largest percentage, present embodiment is the example explanation with Z thermotropism error.
Start CNC milling machine, state is unloaded, and the testing time is 2.5 hours, and fixed frequency is 1min
-1The temperature that is measuring point adopts collection in 1 minute 1 time; Position of tool tip with respect to workpiece Z to deflection also adopt and gathered in 1 minute 1 time, when temperature sensor begins the variation to measuring point T1～T15 collecting temperature, laser displacement sensor also with same frequency measure position of tool tip with respect to workpiece Z to deflection; The time is 2.5 hours during test, and fixed frequency is 1min
-1Gathered altogether 150 times; Each measuring point among measuring point T1～T15 all has the concrete data of 150 temperature; The concrete data of measuring point T1～T15 are in chronological sequence formed the variations in temperature sequence X that together forms measuring point T1～T15; Measure position of tool tip with respect to workpiece Z to deflection 150 data are also arranged, with these deflection data in chronological sequence be combined to form position of tool tip with respect to workpiece Z to deflection sequence Δ, the concrete data of variations in temperature sequence X of measuring point T1～T15 are (unit ℃):
Measuring point T1:
X
1＝{18.3，18.7，18.9，19.2，19.4，19.7，20.0，20.1，20.2，20.4，...27.8，27.8，27.9，27.9，27.9}；
Measuring point T2:
X
2＝{18.4，18.8，19.1，19.4，19.7，19.9，20.1，20.2，20.4，20.5，...27.9，28.0，28.0，28.0，28.0}；
Measuring point T3:
X
3＝{18.3，18.7，19.0，19.3，19.4，19.8，20.0，20.1，20.2，20.4，...27.8，27.8，27.9，27.9，27.9}；
Measuring point T14:
X
14＝{19.1，20.1，21.0，21.75，22.5，23.1，23.5，23.9，24.4，24.7，...40.8，40.8，40.9，41，41，41}；
Measuring point T15:
X
15＝{18.9，19.7，20.5，21.1，21.8，22.2，22.6，22.9，23.3，23.5，...39.9，39.9，40，40，40，40}；
Position of tool tip with respect to workpiece Z to deflection sequence Δ (the μ m of unit):
Δ＝{1.9，2.2，4.1，3.1，4.7，5.6，6.3，6.8，7.5，8.3，...43.8，43.9，43.8，43.7，43.7}；
T11, T2, T9 and T8 variations in temperature sequence X are as shown in Figure 3, position of tool tip with respect to workpiece Z to deflection sequence Δ as shown in Figure 4.
According to above-mentioned variations in temperature sequence X
1～X
15Data carry out after the temperature dependency coefficient calculations measuring point T1～T15 being divided into groups to measuring point T1～T15, are benchmark with measuring point T1, at first calculate the temperature correlation property coefficient between measuring point T1 and the measuring point T2, adopt
Formula calculates, wherein,
The variations in temperature sequence X
1=18.3,18.7,18.9,19.2,19.4,19.7,20.0,20.1,20.2,20.4 ... and 27.8,27.8,27.9,27.9, the maximum among the 27.9} is 27.9, the variations in temperature sequence X
2=18.4,18.8,19.1,19.4,19.7,19.9,20.1,20.2,20.4,20.5 ... and 27.9,28.0,28.0,28.0, the maximum among the 28.0} is 28.0.Data are brought into calculate, obtain C
12Value, C
12>=0.7, therefore, measuring point T2 and measuring point T1 are assigned in same group, promptly the 1st group, and then successively measuring point T3～T15 and measuring point T1 are carried out temperature dependency coefficient calculations, result of calculation C
13, C
14And C
15>=0.7, C
16, C
17, C
18, C
19, C
110, C
111, C
112, C
113, C
114And C
115＜0.7, then T3, T4 and T5 are increased in the 1st group of measuring point T1 place, measuring point T6～T15 of dividing into groups of success then, promptly measuring point T6～T15 carries out the temperature dependency coefficient calculations once more.
Measuring point T6 to fail in dividing into groups is a benchmark, and measuring point T7～T15 and measuring point T6 are carried out the calculating of temperature correlation property coefficient, after calculating, obtains C
67And C
68>=0.7, C
69, C
610, C
611, C
612, C
613, C
614And C
615＜0.7, measuring point T6, T7 and T8 are assigned in the 2nd group, again T9～the T15 of not success grouping is carried out the coefficient calculations of temperature dependency.
Measuring point T9 to fail in dividing into groups is a benchmark, and measuring point T10～T15 and measuring point T9 are carried out the calculating of temperature correlation property coefficient, after calculating, obtains C
910>=0.7, C
911, C
912, C
913, C
914And C
915＜0.7, measuring point T9 and T10 are assigned in the 3rd group, again T11～the T15 of not success grouping is carried out the coefficient calculations of temperature dependency.
Measuring point T11 to fail in dividing into groups is a benchmark, and measuring point T12～T15 and measuring point T11 are carried out the calculating of temperature correlation property coefficient, after calculating, obtains C
1112, C
1113, C
1114And C
1115>=0.7, measuring point T11, T12, T13, T14 and T15 are assigned in the 4th group, at this moment measuring point T1～T15 has divided group, and above-mentioned grouping situation is listed in table 1.
Table 1
| Group number | Form |
| The 1st group | ?T1，T2，T3，T4，T5 |
| The 2nd group | ?T6，T7，T8 |
| The 3rd group | ?T9，T10 |
| The 4th group | ?T11，T12，T13，T14，T15 |
T1～T15 variations in temperature sequence X according to measuring point
1～X
15With position of tool tip with respect to workpiece Z to deflection sequence Δ, the measuring point T1～T15 after dividing into groups is carried out the temperature-sensitive sensitivity calculates, at first the temperature-sensitive sensitivity of the measuring point T1 in the 1st group is calculated the temperature-sensitive sensitivity of T1
Obtain η
1, then other measuring point T2～T5 in the 1st group is carried out the calculating of temperature-sensitive sensitivity, obtain η
2, η
3, η
4And η
5, be η to the maximum through the temperature-sensitive sensitivity that relatively filters out in the 1st group
2, the temperature-sensitive sensitivity of measuring point T2 is maximum in the 1st group.In like manner, again the measuring point in the 2nd group, the 3rd group and the 4th group is carried out the calculating of temperature-sensitive sensitivity, the temperature-sensitive sensitivity of measuring point T8 is maximum in the 2nd group; The temperature-sensitive sensitivity of measuring point T9 is maximum in the 3rd group, and the temperature-sensitive sensitivity of measuring point T11 is maximum in the 4th group, and T2, T8, T9 and T11 are formed the selection group; Wherein, the temperature-sensitive sensitivity of these 4 measuring points is followed successively by T11, T2, T9, T8 from big to small, when numerical control machine heat characteristic monitoring stations arranges that number is 4; T2, T8, T9 and T11 are as numerical control machine heat characteristic monitoring stations, if numerical control machine heat characteristic monitoring stations layout number is pressed the descending select progressively measuring point of temperature-sensitive sensitivity of measuring point less than 4 from the selection group; If the number of selecting is 3; Then selecting T11, T2 and T9 is numerical control machine heat characteristic monitoring stations, if the number of selecting is 2, then selecting T11 and T2 is numerical control machine heat characteristic monitoring stations; If the number of selecting is 1, then selecting T11 is numerical control machine heat characteristic monitoring stations.
And this result also relatively coincide with the actual conditions of lathe, promptly the spindle heat generation amount of conventional CNC milling machine is maximum, and especially the lower ball cover place of main shaft (being the T11 position) heating is the most obvious.
Claims (6)
1. the method for arranging of a numerical control machine heat characteristic monitoring stations may further comprise the steps:
1) on the Digit Control Machine Tool body, chooses the measuring point that is used for the detected temperatures variation; Press the temperature on the frequency collection measuring point; Obtain the temperature value of measuring point under the different time, described temperature value in chronological sequence is combined to form the variations in temperature sequence X of measuring point, in the temperature on gathering measuring point; With respect to the deflection of workpiece, described deflection in chronological sequence is combined to form the deflection sequence Δ of position of tool tip with respect to workpiece according to position of tool tip under the described frequency collection different time;
2) carry out after the temperature dependency coefficient calculations measuring point being divided into each group according to the variations in temperature sequence X of measuring point in the step 1);
Described step 2) may further comprise the steps:
(a) in measuring point, appoint and get a measuring point, note is made i measuring point, is benchmark with i measuring point, gathers N time altogether by frequency, and then the variations in temperature sequence X of i measuring point is expressed as X
i={ x
I1, x
I2... X
IN, will carry out the temperature dependency coefficient calculations except that i measuring point and i measuring point the measuring point, as if being k measuring point except that i the measuring point the measuring point, the variations in temperature sequence X of k measuring point is expressed as X
k={ x
K1, x
K2... X
KN, calculate the temperature correlation property coefficient between i measuring point and k the measuring point, the variations in temperature sequence X
i={ x
I1, x
I2... X
INIn maximum be x
Imax, the variations in temperature sequence X
k={ x
K1, x
K2... X
KNIn maximum be x
Kmax
(b) the temperature dependency coefficients by using between i measuring point and k the measuring point
Formula calculates, wherein,
Obtain C
IkValue;
(c) C
Ik>=0.7, represent that these two measuring points have obvious correlation, are placed on one group, C with k measuring point and i measuring point
Ik＜0.7 these two measuring points of expression do not have obvious correlation; K measuring point and i measuring point are not placed on one group; Then to carrying out the temperature dependency coefficient calculations except that k measuring point and i other measuring points and i measuring point the measuring point; Obtain by with measuring point and i the group that measuring point is formed of i measuring point temperature correlation property coefficient>=0.7, the measuring point that divides into groups that fails is formed remaining set;
(d) appoint from the remaining set of step (c) that to get a point be benchmark, the repeating step (a) and (b) with (c), obtain a new group again, finish grouping up to all measuring points;
3) according to the variations in temperature sequence X of measuring point in the step 1) and position of tool tip with respect to the deflection sequence Δ of workpiece to step 2) in each measuring point in each group carry out the temperature-sensitive sensitivity and calculate; Through relatively after the screening measuring point of the maximum heat susceptibility in every group being selected, form the selection group;
In the step 3), described temperature-sensitive sensitivity is calculated, and may further comprise the steps:
(e) in measuring point, appoint and get a measuring point, note is made h measuring point, gathers N time altogether by frequency, and then the variations in temperature sequence X of h measuring point is expressed as X
h={ x
H1, x
H2... X
HN, position of tool tip is represented Δ={ δ with respect to the deflection sequence Δ of workpiece
1, δ
2... δ
N, the variations in temperature sequence X
h={ x
H1, x
H2... X
HNIn maximum be x
Hmax, position of tool tip is with respect to the deflection sequence Δ={ δ of workpiece
1, δ
2... δ
NIn maximum be δ
Max
(f) the temperature-sensitive sensitivity of h measuring point is defined as
obtains the temperature-sensitive sensitivity of h measuring point after calculating;
(g) repeating step (e) and (f) calculates carry out the temperature-sensitive sensitivity except that h other measuring points the measuring point;
4) if numerical control machine heat characteristic monitoring stations arranges that number equates with measuring point number in the selection group; All measuring points in the then selection group are promptly as numerical control machine heat characteristic monitoring stations; If numerical control machine heat characteristic monitoring stations arranges that number is less than the measuring point number in the selection group; From the selection group, press the descending select progressively and the measuring point of arranging the number equal amount of temperature-sensitive sensitivity of measuring point, the measuring point of selection is promptly as numerical control machine heat characteristic monitoring stations.
2. according to the method for arranging of the numerical control machine heat characteristic monitoring stations described in the claim 1, it is characterized in that the described measuring point that is used for the detected temperatures variation is distributed in lathe bed, spindle jacket and lead screw shaft bearing sleeve place.
3. according to the method for arranging of the numerical control machine heat characteristic monitoring stations described in claim 1 or 2, it is characterized in that the described number that is used for the measuring point that detected temperatures changes is more than or equal to 10.
4. according to the method for arranging of the numerical control machine heat characteristic monitoring stations described in the claim 1, it is characterized in that the temperature on the described measuring point adopts temperature sensor to detect;
Described position of tool tip adopts displacement transducer to detect with respect to the deflection of workpiece.
5. according to the method for arranging of the numerical control machine heat characteristic monitoring stations described in claim 1 or 4, it is characterized in that described position of tool tip is the deflection of position of tool tip with respect to the workpiece Z-direction with respect to the deflection of workpiece.
6. according to the method for arranging of the numerical control machine heat characteristic monitoring stations described in the claim 1, it is characterized in that described frequency is a fixed frequency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2011100653659A CN102179725B (en) | 2011-03-17 | 2011-03-17 | Arrangement method of heat characteristic monitoring measurement points of numerical control machine |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN2011100653659A CN102179725B (en) | 2011-03-17 | 2011-03-17 | Arrangement method of heat characteristic monitoring measurement points of numerical control machine |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102179725A CN102179725A (en) | 2011-09-14 |
| CN102179725B true CN102179725B (en) | 2012-11-07 |
Family
ID=44566116
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN2011100653659A Expired - Fee Related CN102179725B (en) | 2011-03-17 | 2011-03-17 | Arrangement method of heat characteristic monitoring measurement points of numerical control machine |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102179725B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102689234B (en) * | 2012-06-08 | 2014-04-30 | 合肥工业大学 | Device and method for detecting thermal errors of hobbing machine tool |
| CN103837336A (en) * | 2012-11-23 | 2014-06-04 | 江西昌河航空工业有限公司 | Method for detecting heat stability of machine tool main shaft component |
| CN103567815B (en) * | 2013-11-12 | 2016-04-27 | 沈阳机床（集团）设计研究院有限公司 | Based on the test of numerical control machine tool cutting Thermal Error and the evaluation method of milling aperture |
| CN105666244B (en) * | 2016-01-06 | 2018-07-13 | 北京工业大学 | The method of boring bar thermal stretching error temperature point yojan under numerical control borer fuel factor |
| EP3444686B1 (en) * | 2017-08-15 | 2021-12-22 | GF Machining Solutions AG | Method for using a geometrical probe with a spindle of a machine tool, and machine tool configured to carry out such a method |
| CN108608016A (en) * | 2018-04-27 | 2018-10-02 | 北京科技大学 | A kind of discrimination method and its system of electro spindle rapid warm raising |
| CN108445839B (en) * | 2018-05-06 | 2020-08-21 | 北京工业大学 | Machine tool precision sensitivity analysis method based on error increment |
| CN110737959B (en) * | 2019-10-17 | 2021-04-30 | 山东大学 | Synchronous design method for multi-machine tool selection and structural topological configuration in additive manufacturing |
| CN113977353A (en) * | 2021-11-30 | 2022-01-28 | 江苏大卫精工科技有限公司 | Method for testing thermal elongation of boring machine spindle |
| CN116304939A (en) * | 2023-02-22 | 2023-06-23 | 中广核智能科技(深圳)有限责任公司 | Multi-monitoring data evolution rule identification method and system based on waveform alignment |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003019640A (en) * | 2001-07-03 | 2003-01-21 | Okuma Corp | Thermal displacement monitoring device for machine tool |
| CN101794138A (en) * | 2010-04-14 | 2010-08-04 | 华中科技大学 | Dynamic characteristic test and analysis system for numerical control machine tool |
| CN101804581A (en) * | 2010-03-23 | 2010-08-18 | 四川普什宁江机床有限公司 | Implementation method of automatic compensation for thermal deformation of machine tool |
-
2011
- 2011-03-17 CN CN2011100653659A patent/CN102179725B/en not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003019640A (en) * | 2001-07-03 | 2003-01-21 | Okuma Corp | Thermal displacement monitoring device for machine tool |
| CN101804581A (en) * | 2010-03-23 | 2010-08-18 | 四川普什宁江机床有限公司 | Implementation method of automatic compensation for thermal deformation of machine tool |
| CN101794138A (en) * | 2010-04-14 | 2010-08-04 | 华中科技大学 | Dynamic characteristic test and analysis system for numerical control machine tool |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102179725A (en) | 2011-09-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102179725B (en) | Arrangement method of heat characteristic monitoring measurement points of numerical control machine | |
| CN101797704B (en) | Method for thermal deformation error compensation of digital control gear hobbing machine | |
| CN107553220B (en) | A kind of numerically-controlled machine tool composition error real-time compensation method | |
| CN105700473B (en) | A kind of full workbench curved surface thermal error compensation method of precise numerical control machine | |
| TWI448353B (en) | Method and apparatus of machine tools for intelligently compensating thermal error | |
| CN101804581A (en) | Implementation method of automatic compensation for thermal deformation of machine tool | |
| CN105397560B (en) | One kind is dry to cut chain digital control gear hobbing machine bed and workpiece method for thermal deformation error compensation | |
| CN105759719A (en) | Numerical-control machine-tool thermal error prediction method based on unbiased estimation splitting model and system thereof | |
| CN104216334A (en) | Selection optimization method of temperature measurement point combination for positioning errors of numerically-controlled machine tool under thermal effect | |
| CN103034169A (en) | Modeling and compensation method of heat error of numerical control machine tool | |
| CN108188821B (en) | A kind of Ball-screw in NC Machine Tools feed system Thermal Error prediction technique | |
| CN104950808B (en) | Machine tool thermal error compensation method based on augmented naive Bayes network | |
| CN105109094A (en) | High-precision servo press fitting method based on computer | |
| CN103885386A (en) | Gray model thermal error data processing method based on Kalman filtering | |
| CN101629801B (en) | Method for confirming guide rail thermal error of numerical control grinder | |
| CN112475410A (en) | Correlation analysis system and method for milling temperature and multivariate influence factors | |
| CN103273380B (en) | Based on the Thermal Deformation of NC Machine Tool error on-line monitoring method of simulation model | |
| CN105975673B (en) | The Nonlinear dynamics parameter identification method of vibrating type energy collector | |
| CN105706010A (en) | Production line simulation device | |
| CN111708323B (en) | Five-axis small gantry numerical control machining center with thermal deformation error compensation function | |
| CN108646670B (en) | Real-time prediction method for temperature of parts of numerical control machine tool | |
| CN107219819A (en) | A kind of machine tool motion compensation method analyzed based on mean temperature difference | |
| CN103886191A (en) | Straightness compensation method for machine tool body | |
| Polyakov et al. | Static and thermal characteristics of numerically controlled machine tools. | |
| CN109884983B (en) | Method and system for separating thermal error from machine tool positioning error |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| 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: 20121107 Termination date: 20170317 |