CN104216334B - Selection optimization method of temperature measurement point combination for positioning errors of numerically-controlled machine tool under thermal effect - Google Patents

Selection optimization method of temperature measurement point combination for positioning errors of numerically-controlled machine tool under thermal effect Download PDF

Info

Publication number
CN104216334B
CN104216334B CN201410471517.9A CN201410471517A CN104216334B CN 104216334 B CN104216334 B CN 104216334B CN 201410471517 A CN201410471517 A CN 201410471517A CN 104216334 B CN104216334 B CN 104216334B
Authority
CN
China
Prior art keywords
temperature
point
lathe
machine tool
combination
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.)
Active
Application number
CN201410471517.9A
Other languages
Chinese (zh)
Other versions
CN104216334A (en
Inventor
程强
祁卓
李凯
刘志峰
蔡力钢
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201410471517.9A priority Critical patent/CN104216334B/en
Publication of CN104216334A publication Critical patent/CN104216334A/en
Application granted granted Critical
Publication of CN104216334B publication Critical patent/CN104216334B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The invention provides a selection optimization method of a temperature measurement point combination for positioning errors of a numerically-controlled machine tool under thermal effect. The selection optimization method is capable of identifying the influence of the temperature measurement point in each position on the positioning errors of the machine tool based on a grey correlation policy and a rough set theory. The selection optimization method comprises the following steps: k temperature sensors are mounted in special positions of the machine tool to measure the real-time temperature values, changing over time, of the machine tool during operation, and meanwhile, a laser interferometer is used for measuring positioning error values affected by temperatures; n sensitive temperature measurement point positions are screened out by use of the grey correlation policy; the positioning errors and the temperature data of the machine tool are preprocessed according to the principle of the rough set theory and a policy table is established; m feasible temperature point combinations are obtained by use of rough set reduction software; the optimal temperature measurement point combination of the machine tool is identified by virtue of comprehensive analysis. The selection optimization method of the temperature measurement point combination for the positioning errors of the numerically-controlled machine tool under the thermal effect is capable of solving the problem of excessive temperature measurement points or poor compensation model robustness in the positioning error compensation modeling process of the numerically-controlled machine tool.

Description

Under a kind of heat effect with regard to Digit Control Machine Tool, the selection of position error temperature point combination is excellent Change method
Technical field
The present invention relates to the measurement of lathe position error and error compensation modeling in a kind of Computerized Numerical Control Cutting Processes The optimization method of temperature variable combination used, belongs to NC Machine Error analysis technical field.
Background technology
Position accuracy for CNC machine tools be lathe moving component move under digital control system control can reach position essence Degree, in general it is simply that referring to the order of accuarcy that lathe navigates to the point of a knife of cutter impact point in program.Machine finish is It is to be determined by the relative displacement that lathe is fixed a cutting tool and workpiece between eventually.Therefore reduce position error to the processing essence improving lathe Degree is most important.
Machine tool error penalty method due to its good economy performance, feasibility is high and become topmost raising precision, reduction at present The means of error.And set up the key link that position error forecast model is error compensation, the form of model and accuracy are direct The speed of impact error compensation and effect.Set up accurate position error forecast model must obtain related to position error Lathe Temperature Distribution, this is accomplished by each position on lathe and arranges substantial amounts of temperature sensor, for measuring machine bed operating mistake Real time temperature distribution in journey.In general, in heat error compensation the selection of temperature point from hundreds of to several.
However, temperature point excessively not only make arrangement measuring point workload increase, and temperature point be arranged to too close The output signal of adjacent measuring point also can be made to have larger dependency, affect computational accuracy on the contrary.So, select several key temperatures Measuring point is realized being accurately positioned error modeling and is just seemed particular importance, but how to select temperature point to be the modeling of lathe position error One of and the key issue in compensation technique.
Content of the invention
Present invention aims to existing issue, based on rough set theory, on the basis of grey correlation strategy, point Each temperature point significance level that machine tooling position error is affected in analysis machine tool temperature field distribution is it is proposed that according to coarse Set analysis software (ROSETTA) carries out yojan to lathe temperature, error information, and comprehensive analysis are found out and affected spy to position error The sensor combinations of insensitive several measurement points optimize the purpose of lathe position error temperature point to reach, that is, find out optimum Temperature point combines.
For achieving the above object, the technical solution used in the present invention is a kind of with regard to position error under Digit Control Machine Tool heat effect The selection optimization method of temperature point combination, for solving how to optimize temperature point combination during digital control machine tool positioning error compensates Technical problem.
The comprising the following steps that of the method,
Step 1, time dependent temperature variable and position error amount in collection Digit Control Machine Tool running;
First, the critical positions k temperature sensor of installation in Digit Control Machine Tool carries out temperature survey, described Digit Control Machine Tool It is electronic that critical positions mainly include the exemplary position of front-end of spindle and rear end, spindle box body front and back end and upper end, three axial filament thick sticks Machine, bearing, guide rail, operating position.
Then, first measurement and positioning error under lathe cold conditions (just starting shooting), after measurement, quick moving movement axle makes lathe Temperature raise, then measure again, again temperature rise so repeat to tend towards stability to each temperature change of lathe, that is, lathe reaches thermal balance shape State terminates to measure.Can be obtained by running lathe:The change of temperature T that 1. temperature sensor of k position records t in time Amount T { T1(t), T2(t) ..., TK(t)};2. lathe position error amount Y (t) that laser interferometer records;
Step 2, application grey correlation analysis filters out m sensitive temperature point position:
Set up reference sequence (position error data) using grey correlation analysis and compare ordered series of numbers (k temperature point data) Between coefficient of association ξ0kWith degree of association γ0k, and these degrees of association are arranged in order from big to small, represent these measuring points respectively Position temperature change produces the impact size of position error to lathe;Set a threshold gamma ', in general, defining threshold value is:When degree of association γ0kDuring more than γ ', temperature point position is retained;And The temperature change of remaining position is very small on position error impact, is all cast out, will successfully be reduced to m by k temperature point Individual sensitive measuring point T ' { T '1(t), T '2(t) ..., T 'm(t)}.
Step 3, according to the principle of rough set theory, carries out pre- place to position error under lathe heat effect and temperature data Reason, constitutes a decision table;
Using the temperature of the k position surveyed as conditional attribute C, i.e. C={ T1(t), T2(t) ..., Tk(t) }, surveyed Position error displacement attribute D as a result, i.e. D={ Y (t) }, thus establishing a system decision-making table K=(U, C ∪ D), and This decision table is established as an Excel table.
Step 4, draws n feasible temperature point combination using Rough Set Analysis software (ROSETTA):
The Excel table of the system decision-making being established above table K=(U, C ∪ D) is input to Rough Set Analysis software (ROSETTA) in, carry out Data Reduction process by after Data-parallel language, Data Discretization, obtain the feasible temperature of n kind and survey Point combination, machine tool temperature field distribution situation can be intactly expressed in the combination of these temperature points.
Step 5, comprehensive analysis identification lathe Optimal Temperature measuring point combination, complete to select optimization method:
M sensitive temperature measuring point of Integrated comparative yojan and n temperature point combination, filter out and comprise sensitive temperature measuring point Most and degree of association highest temperature point combines, and as required Optimal Temperature measuring point combines.
Compared with prior art, the invention has the beneficial effects as follows:The present invention misses in measurement Digit Control Machine Tool temperature field and positioning On the basis of difference, go out m and the high transducer arrangements point of lathe position error dependency using grey correlation Policy Filtering, then will Temperature data and position error data are established as decision table according to rough set theory, pass through Rough Set Analysis software afterwards (ROSETTA) yojan is optimized to decision table, draws n feasibility temperature point combination;Find out finally by comprehensive analysis Optimal Temperature measuring point combines, and determines the installation site that several transducer arrangements points compensate as the modeling of lathe position error.Compare Traditional method finding lathe key temperatures location point by many experiments number of times based on engineering judgement, the present invention has and saves time Efficiently, save temperature sensor, simplify modeling process, the robustness of lathe Model of locating error and accuracy high the advantages of.
Brief description
Fig. 1 is workflow diagram of the present invention;
Fig. 2 is numerically controlled lathe schematic diagram and temperature sensor thermometric arrangement schematic diagram;
Fig. 3 is numerically controlled lathe position error measurement arrangement schematic diagram;
In figure:1st, the generating laser of laser interferometer, 2, the main shaft of lathe, 3, workbench, 4, the connecing of laser interferometer Receive device.
Specific embodiment
The present invention is described further with implementation process below in conjunction with the accompanying drawings.
As Figure 1-3, lathe location error compensation of the present invention models the selection optimization method of temperature point, and it is A kind of combined selection method based on grey correlation analysis and rough set theory, realizes according to following steps:
Consider first to produce the correlative factor of position error under heat effect, move back and forth including machine tool motion part and produce heat Amount, the impact of motor running heating, lathe each part heat-generation and heat and ambient temperature, analysis according to this determines harvester bed temperature The position of sensor in degrees of data experiment.As shown in table 1, table 1 is referred in 29 transducer arrangements positions:Numbering is 1,2,3, 4th, 6,8,14 and No. 18 sensors are arranged on front-end of spindle and the exemplary position of rear end, each four of front end rear end, and each biography Sensor spacing is equal, is looped around on axle head;Sensor is equidistantly installed, it is to avoid distance closely interferes very much, away from detection too far away Not comprehensive, installed identical with lower sensor.5th, 7,11 and No. 15 are arranged on spindle box body front and back end and upper end, 12,19, 21st, 22,25,28 and No. 29 are arranged on X, Y, Z tri- axial filament thick stick motor, bearing, shaft coupling and nut, 20,23,24,26 and No. 27 are arranged on X, Y, Z axis direction guiding rail, and 9,10,13,16 and No. 17 are arranged in vertical slide plate, workbench and ambient temperature.
1 29 transducer arrangements location tables of table
Laser interferometer is fixed on lathe, concrete installation is as shown in Figure 3:The generating laser of laser interferometer 1st, the receptor 4 of the main shaft 2 of lathe, workbench 3, laser interferometer;First generating laser 1 is arranged on workbench 3 front end Smooth place, afterwards the receptor 4 of laser interferometer is arranged on the main shaft 2 of lathe.Laser interference is adjusted after installation The laser head of the generating laser 1 of instrument, makes measurement axis point-blank or parallel with the axis of machine tool movement, will light path Harmonize straight;During standby bed operating, on request the relevant parameter of lathe is measured.
Then run lathe and carry out data acquisition.The data recycling grey correlation analysis to collect, filters out m correlation The big transducer arrangements point of coefficient.Then with Rough Set Analysis software (ROSETTA), temperature data and position error data are entered Row is processed, and carries out yojan to temperature point, obtains n feasible temperature point combination.Finally, comprehensive analysis m temperature is sensitive The temperature point combination of point and n, filters out and comprises sensitive temperature measuring point at most and the combination of degree of association highest temperature point, that is, For required Optimal Temperature measuring point combination.
The step that implements of the present embodiment is:
Step 1, time dependent temperature variable and position error amount in collection Digit Control Machine Tool running:
First, the critical positions in Digit Control Machine Tool install k (can select from above-mentioned 29 temperature sensor location) temperature Degree sensor carries out temperature survey, laser interferometer is fixed on lathe and carries out position error measurement (generating laser It is fixed on guide rail, receptor is fixed on lathe cutter saddle);
Then, first measurement and positioning error under lathe cold conditions (just starting shooting), after measurement, quick moving movement axle makes lathe Temperature raise, then measure again, again temperature rise so repeat to tend towards stability to each temperature change of lathe, that is, lathe reaches thermal balance shape State terminates to measure.Can be obtained by running lathe:1. the temperature that the temperature sensor of k position records is measured over time T{T1(t), T2(t) ..., Tk(t)};2. lathe position error amount Y (t) that laser interferometer records;
Step 2, application grey correlation analysis filters out m sensitive temperature point position:
Grey correlation analysis is the similar close degree mathematical theory institute according to characteristic parameter each in system between serial The systematic analysiss carrying out.
In processing cutting process, lathe position error changes obvious, choosing at the discontinuous point that lathe runs The detected value at material time node take lathe in start, shutting down carrys out modeling analysis, also contemplates Grey Relational Model simultaneously Even time interval to data acquisition.Specifically, l group can be taken (typically to take with certain access time unit interval node Twenty or thirty group is advisable) data.
Because the implication of each influence factor and purpose are different, thus desired value generally has different dimension data levels, If the data between two sequences differs greatly in size, the effect of decimal value sequence is easily covered by big sequence of values, For the ease of comparing, ensure, between each factor, there is equivalence and same sequence, initial data need to be processed, be allowed to nondimensionalization And normalization.In general, generally using first value conversion, equalization conversion and extreme differenceization, three kinds can be converted to initial data Processing method.Choose alternative approach here according to corresponding requirements.
If reference sequence is position error data Y (t)={ Y (j) | j=1,2 ..., l }, comparing ordered series of numbers is k temperature survey Point data Ti={ Ti(j) | i=1,2 ..., k;J=1,2 ..., l }.Then Y (t) is to TiCoefficient of association in jth point:
In formula, Δ0iJ () is jth point Y (t) and TiThe absolute value of the difference of (t), Δ0i(j)=| Y (t)-Ti(t)|; miniminjΔ0iJ () is the two poles of the earth lowest difference;maximaxjΔ0iJ () is that the two poles of the earth are maximum poor;ρ is resolution ratio, ρ ∈ [0,1], one As take ρ=0.5, in concrete operation, ρ value can be adjusted, to increase contrast according to the degree of association between each data sequence The resolution capability of analysis.
Degree of association γ between two sequences0iThe coefficient of association ξ in available two each moment of sequence0iJ the meansigma methodss of () are counting Calculate, i.e.
In formula, γ0iThe degree of association for subsequence i and auxiliary sequence;L is the data amount check of two comparative sequences.
Finally degree of association γ to same auxiliary sequence by each subsequence0iSequentially line up by size, i.e. composition association Sequence, it directly reflects each subsequence " primary and secondary " relation to same auxiliary sequence, that is, represents that these point position temperature become Change the impact size that lathe is produced with position error.Now, set a threshold value When degree of association γ0kDuring more than γ ', temperature point position is retained, and remaining temperature point position is then cast out, will k Individual temperature point is successfully reduced to the sensitive measuring point T ' { T ' of m1(t), T '2(t) ..., T 'm(t)}.
One sequence of grey correlation is established by above-mentioned analytical calculation, γ is taken to inteerelated order0kMore than the factor of γ ', then have
T1> T6> T8> T10> T13> T5> T14> T18> T11> T19.
This sequence of grey correlation represents the arranging situation of the temperature point maximum on position error impact.
Step 3, according to the principle of rough set theory, carries out pre- place to position error under lathe heat effect and temperature data Reason, constitutes a decision table:
Using the temperature of the k position surveyed as conditional attribute C, i.e. C={ T1(t), T2(t) ..., Tk(t) }, surveyed Position error displacement attribute D as a result, i.e. D={ Y (t) }, thus establishing a system decision-making table K=(U, C ∪ D), and This decision table is established as an Excel table.
Step 4, draws n feasible temperature point combination using Rough Set Analysis software (ROSETTA):
The Excel table of the system decision-making being established above table K=(U, C ∪ D) is input to Rough Set Analysis software (ROSETTA) in, carry out Data Reduction process by after Data-parallel language, Data Discretization, obtain the feasible temperature of n kind and survey Point combination, the combination of these temperature points can be than more complete expression machine tool temperature field distribution situation.
Following feasible combination can be obtained through Rough Set Analysis software (ROSETTA) analysis:
{T1T5T6T8T10T13T14}、{T1T6T8T10T13T14T18}、{T5T6T8T10T13T14T18}、 {T1T5T6T8T10T11T19}、{T1T5T6T8T10T13T18}、{T1T5T6T8T10T13T14T18}.This six kinds Combination can react whole Digit Control Machine Tool thermo parameters method situation and on position error impact situation than more complete.
Step 5, comprehensive analysis identification lathe Optimal Temperature measuring point combination, complete to select optimization method:
M sensitive temperature measuring point of Integrated comparative yojan and n temperature point combination, filter out and comprise sensitive temperature measuring point Most and degree of association highest temperature point combines, and as required Optimal Temperature measuring point combines.
By comparing sequence of grey correlation and temperature point combination it can be deduced that an Optimal Temperature measuring point combines:{T1, T5, T6, T8, T10, T13, T14}.
After above-mentioned 5 steps complete, an optimum point position combination can be obtained, reached cost-effective, simple Change the purpose that location error compensation measures experimental implementation and improves the robustness of Model of locating error.

Claims (3)

1. under a kind of heat effect with regard to Digit Control Machine Tool position error temperature point combination selection optimization method it is characterised in that:
The comprising the following steps that of the method,
Step 1, time dependent temperature variable and position error amount in collection Digit Control Machine Tool running;
First, the critical positions in Digit Control Machine Tool are installed k temperature sensor and are carried out temperature survey, described Digit Control Machine Tool important Position mainly includes exemplary position, spindle box body front and back end and upper end, three axial filament thick stick motor, the axle of front-end of spindle and rear end Hold, guide rail, operating position;
Then, first measurement and positioning error under lathe cold conditions, after measurement, quick moving movement axle makes lathe temperature raise, then Measure again, again temperature rise so repeat to tend towards stability to each temperature change of lathe, that is, lathe reach thermal equilibrium state terminate measure;Logical Cross operation lathe can obtain:Variable quantity T { the T of temperature T that 1. temperature sensor of k position records t in time1(t), T2 (t) ..., Tk(t)};2. lathe position error amount Y (t) that laser interferometer records;
Step 2, application grey correlation analysis filters out m sensitive temperature point position:
Set up reference sequence using grey correlation analysis and compare the coefficient of association ξ between ordered series of numbers0kWith degree of association γ0k, and by this A little degrees of association are arranged in order from big to small, represent that these point position temperature changes produce the impact of position error to lathe respectively Size;Set a threshold gamma ', defining threshold value is:When degree of association γ0kIt is more than During γ ', temperature point position is retained;And the temperature change of remaining position is very small on position error impact, all given up Go, will successfully be reduced to the individual sensitive measuring point T ' { T ' of m by k temperature point1(t), T '2(t) ..., T 'm(t)};
Step 3, according to the principle of rough set theory, carries out pretreatment, structure to position error under lathe heat effect and temperature data Become a decision table;
Using the temperature of the k position surveyed as conditional attribute C, i.e. C={ T1(t), T2(t) ..., Tk(t) }, the positioning surveyed Error displacement attribute D as a result, i.e. D={ Y (t) }, thus establishing a system decision-making table K=(U, C ∪ D), and by this Decision table is established as an Excel table;
Step 4, draws n feasible temperature point combination using Rough Set Analysis software:
The Excel table of the system decision-making being established above table K=(U, C ∪ D) is input in Rough Set Analysis software, by number Carry out Data Reduction process according to after polishing, Data Discretization, obtain the feasible temperature point combination of n kind, these temperature points Machine tool temperature field distribution situation can be intactly expressed in combination;
Step 5, comprehensive analysis identification lathe Optimal Temperature measuring point combination, complete to select optimization method:
M sensitive temperature measuring point of Integrated comparative yojan and the combination of n temperature point, filter out that to comprise sensitive temperature measuring point most And degree of association highest temperature point combines, as required Optimal Temperature measuring point combines;
Described lathe location error compensation models the selection optimization method of temperature point, it be a kind of based on grey correlation analysis and The combined selection method of rough set theory, step 1 is further comprising the steps of,
Consider first to produce the correlative factor of position error under heat effect, move back and forth including machine tool motion part and produce heat, electricity Motivation runs the impact of heating, lathe each part heat-generation and heat and ambient temperature, and analysis according to this determines the harvester bed temperature number of degrees Factually test the position of middle sensor;As shown in table 1, table 1 is referred in 29 transducer arrangements positions, numbering is 1,2,3,4,6, 8th, 14 and No. 18 sensors are arranged on front-end of spindle and the exemplary position of rear end, each four of front end rear end, and each sensor Spacing is equal, is looped around on axle head;Sensor is equidistantly installed, it is to avoid distance closely interferes very much, incomplete away from detection too far away Face, is installed identical with lower sensor;5th, 7,11 and No. 15 are arranged on spindle box body front and back end and upper end;12、19、21、22、 25th, 28 and No. 29 are arranged on X, Y, Z tri- axial filament thick stick motor, bearing, shaft coupling and nut;20th, 23,24,26 and No. 27 peaces It is contained on X, Y, Z axis direction guiding rail;9th, 10,13,16 and No. 17 are arranged in vertical slide plate, workbench and ambient temperature;
1 29 transducer arrangements location tables of table
Laser interferometer is fixed on lathe:The generating laser (1) of laser interferometer, the main shaft (2) of lathe, work Platform (3), the receptor (4) of laser interferometer;First generating laser (1) is arranged on the smooth place of workbench (3) front end, it Afterwards the receptor (4) of laser interferometer is arranged on the main shaft (2) of lathe;The laser of laser interferometer is adjusted after installation The laser head of emitter (1), makes measurement axis point-blank or parallel with the axis of machine tool movement, will light path harmonize Directly;During standby bed operating, on request the relevant parameter of lathe is measured;
Then run lathe and carry out data acquisition;The data recycling grey correlation analysis to collect, filters out m correlation coefficient Big transducer arrangements point;Then with Rough Set Analysis software, temperature data and position error data are processed, to temperature Measuring point carries out yojan, obtains n feasible temperature point combination;Finally, the temperature of comprehensive analysis m temperature sensitive point and n Measuring point combines, and filters out and comprises sensitive temperature measuring point at most and the combination of degree of association highest temperature point, as required optimum Temperature point combines.
2. under a kind of heat effect with regard to Digit Control Machine Tool according to claim 1, the selection of position error temperature point combination is excellent Change method it is characterised in that:
Step 2 also includes applying grey correlation analysis to filter out m sensitive temperature point position:
Grey correlation analysis is that the similar close degree between serial is carried out with mathematical theory according to characteristic parameter each in system Systematic analysiss;
In processing cutting process, lathe position error changes obvious, selection machine at the discontinuous point that lathe runs Detected value at start, the material time node of shutdown for the bed carrys out modeling analysis, also contemplates Grey Relational Model logarithm simultaneously According to the even time interval of collection, with certain access time unit interval node, take l group data;
Because the implication of each influence factor and purpose are different, thus desired value generally has different dimension data levels, if Data between two sequences differs greatly in size, then the effect of decimal value sequence is easily covered by big sequence of values, in order to It is easy to compare, ensures, between each factor, there is equivalence and same sequence, initial data need to be processed, be allowed to nondimensionalization and return One change;
If reference sequence is position error data Y (t)={ Y (j) | j=1,2 ..., l }, comparing ordered series of numbers is k temperature point number According to Ti={ Ti(j) | i=1,2 ..., k;J=1,2 ..., l };Then Y (t) is to TiCoefficient of association in jth point:
ξ 0 i ( j ) = min i min j Δ 0 i ( j ) + ρmax i max j Δ 0 i ( j ) Δ 0 i ( j ) + ρmax i max j Δ 0 i ( j ) - - - ( 1 )
In formula, Δ0iJ () is jth point Y (t) and TiThe absolute value of the difference of (t), Δ0i(j)=| Y (t)-Ti(t)|;miniminjΔ0i J () is the two poles of the earth lowest difference;maximaxjΔ0iJ () is that the two poles of the earth are maximum poor;ρ is resolution ratio, ρ=0.5, in concrete operation, root According to the degree of association between each data sequence, ρ value is adjusted, to increase the resolution capability of relative analyses;
Degree of association γ between two sequences0iThe coefficient of association ξ in available two each moment of sequence0iJ the meansigma methodss of () are calculating, i.e.
γ 0 i = 1 l Σ j = 1 l ξ 0 i ( j ) , j = 1 , 2 , ... , l - - - ( 2 )
In formula, γ0iThe degree of association for subsequence i and auxiliary sequence;L is the data amount check of two comparative sequences;
Finally degree of association γ to same auxiliary sequence by each subsequence0iSequentially line up by size, that is, form inteerelated order, it Directly reflect the primary-slave relation to same auxiliary sequence for each subsequence, that is, represent these point position temperature changes to machine Bed produces the impact size of position error;Now, set a threshold valueWork as association Degree γ0kDuring more than γ ', temperature point position is retained, and remaining temperature point position is then cast out, will k temperature Measuring point is successfully reduced to the sensitive measuring point T ' { T ' of m1(t), T '2(t) ..., T 'm(t)};
One sequence of grey correlation is established by above-mentioned analytical calculation, γ is taken to inteerelated order0kMore than the factor of γ ', then have
T1> T6> T8> T10> T13> T5> T14> T18> T11> T19
This sequence of grey correlation represents the arranging situation of the temperature point maximum on position error impact;
Step 3 also includes the principle according to rough set theory, carries out pre- place to position error under lathe heat effect and temperature data Reason, constitutes a decision table:
Using the temperature of the k position surveyed as conditional attribute C, i.e. C={ T1(t), T2(t) ..., Tk(t) }, the positioning surveyed Error displacement attribute D as a result, i.e. D={ Y (t) }, thus establishing a system decision-making table K=(U, C ∪ D), and by this Decision table is established as an Excel table;
Step 4 also includes drawing n feasible temperature point combination using Rough Set Analysis software:
The Excel table of the system decision-making being established above table K=(U, C ∪ D) is input in Rough Set Analysis software, by number Carry out Data Reduction process according to after polishing, Data Discretization, obtain the feasible temperature point combination of n kind, these temperature points Combination can be than more complete expression machine tool temperature field distribution situation;
Obtain following feasible combination through Rough Set Analysis software analysis:
{T1T5T6T8T10T13T14}、{T1T6T8T10T13T14T18}、{T5T6T8T10T13T14T18}、{T1 T5T6T8T10T11T19}、{T1T5T6T8T10T13T18}、{T1T5T6T8T10T13T14T18};
This six kinds of combinations can react whole Digit Control Machine Tool thermo parameters method situation and on position error impact feelings than more complete Condition;
Step 5 also includes comprehensive analysis identification lathe Optimal Temperature measuring point combination, completes to select optimization method:
M sensitive temperature measuring point of Integrated comparative yojan and the combination of n temperature point, filter out that to comprise sensitive temperature measuring point most And degree of association highest temperature point combines, as required Optimal Temperature measuring point combines;
By comparing sequence of grey correlation and temperature point combination, draw an Optimal Temperature measuring point combination:{T1, T5, T6, T8, T10, T13, T14}.
3. under a kind of heat effect with regard to Digit Control Machine Tool according to claim 1, the selection of position error temperature point combination is excellent Change method it is characterised in that:In described step 2, usual to initial data using first value conversion, equalization conversion and extreme difference Three kinds of processing methods of conversion, choose alternative approach according to corresponding requirements.
CN201410471517.9A 2014-09-16 2014-09-16 Selection optimization method of temperature measurement point combination for positioning errors of numerically-controlled machine tool under thermal effect Active CN104216334B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410471517.9A CN104216334B (en) 2014-09-16 2014-09-16 Selection optimization method of temperature measurement point combination for positioning errors of numerically-controlled machine tool under thermal effect

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410471517.9A CN104216334B (en) 2014-09-16 2014-09-16 Selection optimization method of temperature measurement point combination for positioning errors of numerically-controlled machine tool under thermal effect

Publications (2)

Publication Number Publication Date
CN104216334A CN104216334A (en) 2014-12-17
CN104216334B true CN104216334B (en) 2017-02-22

Family

ID=52097950

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410471517.9A Active CN104216334B (en) 2014-09-16 2014-09-16 Selection optimization method of temperature measurement point combination for positioning errors of numerically-controlled machine tool under thermal effect

Country Status (1)

Country Link
CN (1) CN104216334B (en)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN107219566B (en) * 2017-05-02 2019-05-07 西北工业大学 Cloud prediction and forecasting procedure based on GM (1,1) gray model
CN108133110B (en) * 2017-12-28 2021-07-02 重庆大学 Method for measuring and calculating temperature field of tool rest unit in dry cutting and hobbing process
CN108188821B (en) * 2018-01-08 2019-04-26 东北大学 A kind of Ball-screw in NC Machine Tools feed system Thermal Error prediction technique
CN108803485B (en) * 2018-07-05 2020-08-14 大连理工大学 I5iport protocol-based external thermal error compensation method for numerical control machine tool
CN109765846B (en) * 2019-03-11 2022-03-18 哈尔滨理工大学 Method for selecting thermal error modeling temperature sensitive points of double-turntable five-axis numerical control machine tool
CN110276400B (en) * 2019-06-24 2021-08-17 重庆大学 Tool and fixture optimization method based on AHP-grey correlation analysis algorithm
CN110579999A (en) * 2019-08-27 2019-12-17 东莞市巨冈机械工业有限公司 z-direction zero drift error compensation method based on triaxial drilling and tapping numerical control machine tool, electronic equipment and computer readable storage medium
CN111338291B (en) * 2020-04-07 2023-07-14 湖北文理学院 Abbe error compensation method and system based on machine tool positioning precision measurement
CN112926200B (en) * 2021-02-23 2023-03-14 内蒙古工业大学 Numerical control machine tool spindle thermal state information mining and analyzing method and device and electronic equipment
CN113297723B (en) * 2021-04-22 2022-11-08 哈尔滨理工大学 Mean shift-grey correlation analysis-based optimization method for electric spindle temperature measurement point
CN113218536B (en) * 2021-05-14 2024-04-02 湘潭大学 Method for selecting temperature measuring points of battery pack of electric automobile
CN113238480B (en) * 2021-05-17 2022-04-26 合肥工业大学 Parameterized regulating and controlling system and method for metal cutting machining cooling gas jet
CN116437649B (en) * 2023-06-13 2023-09-22 浙江德塔森特数据技术有限公司 Machine room safety operation and maintenance method and device based on blockchain and readable storage medium
CN117170308B (en) * 2023-11-02 2024-02-02 中国机械总院集团云南分院有限公司 Machine tool dynamic error compensation method and system based on instruction sequence analysis

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0735443A3 (en) * 1995-03-28 1997-12-29 Siemens Aktiengesellschaft Neural network for compensating a position deviation in machining equipments for work pieces
JP4082598B2 (en) * 2003-08-01 2008-04-30 キタムラ機械株式会社 Method and apparatus for correcting thermal displacement of numerically controlled machine tool
CN101620431A (en) * 2009-06-22 2010-01-06 江俊逢 Digital control method and system for carrying out real-time compensation for thermal deformation errors
CN102096672A (en) * 2009-12-09 2011-06-15 西安邮电学院 Method for extracting classification rule based on fuzzy-rough model
CN103576604A (en) * 2012-07-25 2014-02-12 上海睿涛信息科技有限公司 Dynamic real-time compensation system for positioning errors of numerical-control machine tool
CN103926874A (en) * 2013-12-03 2014-07-16 江苏齐航数控机床有限责任公司 Selection optimization method of numerically-controlled machine tool thermal error compensation modeling temperature measuring point combination
CN103984287A (en) * 2014-03-12 2014-08-13 江苏齐航数控机床有限责任公司 Numerically-controlled machine tool thermal error compensation grey neural network modeling method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0735443A3 (en) * 1995-03-28 1997-12-29 Siemens Aktiengesellschaft Neural network for compensating a position deviation in machining equipments for work pieces
JP4082598B2 (en) * 2003-08-01 2008-04-30 キタムラ機械株式会社 Method and apparatus for correcting thermal displacement of numerically controlled machine tool
CN101620431A (en) * 2009-06-22 2010-01-06 江俊逢 Digital control method and system for carrying out real-time compensation for thermal deformation errors
CN102096672A (en) * 2009-12-09 2011-06-15 西安邮电学院 Method for extracting classification rule based on fuzzy-rough model
CN103576604A (en) * 2012-07-25 2014-02-12 上海睿涛信息科技有限公司 Dynamic real-time compensation system for positioning errors of numerical-control machine tool
CN103926874A (en) * 2013-12-03 2014-07-16 江苏齐航数控机床有限责任公司 Selection optimization method of numerically-controlled machine tool thermal error compensation modeling temperature measuring point combination
CN103984287A (en) * 2014-03-12 2014-08-13 江苏齐航数控机床有限责任公司 Numerically-controlled machine tool thermal error compensation grey neural network modeling method

Also Published As

Publication number Publication date
CN104216334A (en) 2014-12-17

Similar Documents

Publication Publication Date Title
CN104216334B (en) Selection optimization method of temperature measurement point combination for positioning errors of numerically-controlled machine tool under thermal effect
CN105759719B (en) A kind of numerical control machining tool heat error prediction technique and system splitting model based on unbiased esti-mator
CN103926874B (en) The selection optimization method of numerical control machine heat error compensation modeling temperature point combination
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
CN102854841B (en) Shape and position error in-situ compensating and processing method for curved surface parts
CN102672527B (en) Full working stroke thermal error compensation method of numerically-controlled machine tool feeding system and implementation system thereof
Wang et al. Thermal sensor selection for the thermal error modeling of machine tool based on the fuzzy clustering method
CN102736558A (en) Numerical control machine tool thermal error real-time compensation modeling method based on time series algorithm
CN106444628A (en) Numerically-controlled machine tool spindle thermal extension error real-time compensation method
CN106842922B (en) Numerical control machining error optimization method
CN105269404B (en) Numerically-controlled machine tool point of a knife dynamic characteristic accuracy detecting device and its method
CN102658499A (en) Spindle thermal distortion compensation method for precision horizontal machining center
CN110039373B (en) Method and system for predicting thermal deformation of spindle of numerical control machine tool
CN101804581A (en) Implementation method of automatic compensation for thermal deformation of machine tool
CN105666244B (en) The method of boring bar thermal stretching error temperature point yojan under numerical control borer fuel factor
CN103984287A (en) Numerically-controlled machine tool thermal error compensation grey neural network modeling method
CN104950808B (en) Machine tool thermal error compensation method based on augmented naive Bayes network
Zhang et al. Temperature variable optimization for precision machine tool thermal error compensation on optimal threshold
CN109765846A (en) The selection method of double turntable five-axle number control machine tool thermal error modeling temperature sensitive points
CN106483928A (en) A kind of dry chain digital control gear hobbing machine bed radial heat distortion error of cutting is in machine real-time compensation method
CN105397560A (en) Thermal deformation error compensation method for dry-cutting numerically-controlled gear hobbing machine tool and workpieces
CN107066775A (en) A kind of Forecasting Methodology of cutter turning temperature rise average
CN104166373A (en) Online detecting method and system for numerical control machine tool machining error
CN102193045A (en) Method for checking consistency of analog quantity signal channels

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