CN107368637A - Precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory - Google Patents

Precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory Download PDF

Info

Publication number
CN107368637A
CN107368637A CN201710544527.4A CN201710544527A CN107368637A CN 107368637 A CN107368637 A CN 107368637A CN 201710544527 A CN201710544527 A CN 201710544527A CN 107368637 A CN107368637 A CN 107368637A
Authority
CN
China
Prior art keywords
error
interval
geometric
theory
lathe
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.)
Pending
Application number
CN201710544527.4A
Other languages
Chinese (zh)
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.)
Tianjin University
Original Assignee
Tianjin University
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 Tianjin University filed Critical Tianjin University
Priority to CN201710544527.4A priority Critical patent/CN107368637A/en
Publication of CN107368637A publication Critical patent/CN107368637A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Abstract

The invention discloses a kind of precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory, it comprises the following steps:(1) it is theoretical based on Multibody Kinematics, the mapping model established between each geometric error of lathe and lathe space error;(2) interval theory is based on, carries out carrying out sensitivity analysis to each error source;(3) interval theory is based on, establishes range optimization distribution model;(4) genetic algorithm is based on, is solved with MATLAB softwares so as to realize geometric accuracy optimization distribution.The error source that the present invention makes to be difficult to realize is relaxed to the feasible section of maximum, greatly reduces production cost, improves lathe production efficiency of assembling.

Description

Precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory
Technical field
The present invention relates to Digit Control Machine Tool geometric accuracy optimization distribution field, especially a kind of precision based on interval theory crouches Formula machining center geometric accuracy optimizing distribution method.
Background technology
Precise numerical control machine is the machine-tool of basic manufacturing equipment, and its design directly affects manufacturing with working ability Development.The quality of high-grade, digitally controlled machine tools quality, the accuracy of manufacture and assembly precision of lathe are directly shown, and machine tool accuracy designs With the quality of precision distribution, there is significant impact to its quality.
For machine tool accuracy assignment problem, domestic and foreign scholars carried out it is relatively broad and deep probe into, achieve certain Progress.Mainly have component tolerances weight European norm maximum method, Monte Carlo Method, dimension chain method, based on tolerance-cost The methods of precision distribution method of model, these methods are directed to parallel institution mostly, and can not directly instruct high-grade numerical control machine The manufacture assembling process of bed.At present, domestic Machine Tool Enterprises lead to when carrying out precise numerical control machine accuracy Design and precision distribution Frequently with the method for traditional analogy, the inquiry of lathe handbook and designer's estimation, do not tackle the problem at its root, design Lathe can too rely on the experience of designer unavoidably, can not be formed it is scientific, theorize, the design method of systematization.
The space error of precise numerical control machine point to be cut to point of a knife point directly reflects the machining accuracy of lathe, and space Error is by caused by each kinematic pair geometric error coupling.Therefore, how on the premise of complete machine space error is met, Rationally, it is economical and the problem of rapidly optimization distribution geometric error is current precise numerical control machine assembling manufacturing urgent need to resolve.
The content of the invention
A kind of precision horizontal processing based on interval theory is provided it is an object of the invention to solve above-mentioned technical problem Center geometric accuracy optimizing distribution method, its error source that can make to be difficult to realize are relaxed to the feasible section of maximum.
In order to solve the above-mentioned technical problem, the present invention adopts the following technical scheme that:
A kind of precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory, comprises the following steps:
Step (1):Digit Control Machine Tool is regarded as to the multi-body system being in series by multiple rigid bodies first;Secondly many body system is used System kinematical theory describes the topological relation on Digit Control Machine Tool between each moving component;Reuse homogeneous coordinate transformation matrix phase Multiply the pose transformational relation between adjacent component disjunctor coordinate system on expression Digit Control Machine Tool;Finally set up each geometry of Digit Control Machine Tool Mapping model between error and lathe space error;
Step (2):Based on interval theory, sensitivity is defined according to the interval extension factor, line sensitivity of going forward side by side is analyzed with anti- Reflect influence degree of each error source to complete machine space error;
Step (3):Using each error source interval width as design variable, complete machine space body diagonal error be constraints, Each error source sensitivity coefficient and the algebraical sum of interval width inverse square product are object function, establish range optimization distribution mould Type;
Step (4):Based on genetic algorithm, solved with MATLAB softwares so as to realize geometric accuracy optimization distribution.
The beneficial effects of the invention are as follows:The present invention proposes the geometry based on interval theory for precise horizontal machining center Precision optimizing distribution method, the complete machine space error that lathe geometric error is established based on Multibody Kinematics theory map mould Type, sensitivity analysis is carried out based on interval theory and establishes range optimization distribution model, be finally based on genetic algorithm utilization MATLAB softwares distribute to solve so as to realize that geometric accuracy optimizes, and the error source for making to be difficult to realize is relaxed to the feasible region of maximum Between, production cost is greatly reduced, improves lathe production efficiency of assembling.
Brief description of the drawings
Fig. 1 is the flow of the precise horizontal machining center geometric accuracy optimizing distribution method of the invention based on interval theory Figure;
Fig. 2 is u2000/800H precise horizontal machining center three-dimensional model diagrams;
Fig. 3 is principle of genetic algorithm figure;
Fig. 4 is each geometric error source precision optimizing allocation result.
Wherein, 1- lathe beds 2-X is to slide unit 3-B to turntable 4- workpiece 5- column 6- main spindle box 7- main shaft 8- cutters.
Embodiment
The present invention is further detailed explanation with reference to the accompanying drawings and examples:
As shown in figure 1, the precise horizontal machining center geometric accuracy optimization distribution provided by the invention based on interval theory Method, comprise the following steps:
(1) mapping model established between each geometric error of lathe and lathe space error;
The geometric error of Digit Control Machine Tool, according to whether it is related to machine tool motion position, position correlation geometric error can be divided into And position independence geometric error (PIGEs) (PDGEs).The present invention is by taking u2000/800H precise horizontal machining centers as an example, model As shown in Fig. 2 being illustrated to the method for the present invention, above-mentioned machining center is common in the case where only considering translation shaft geometric error There are 21 geometric errors:9 site errors related to three axis coordinate positions and 9 angular errors, and position are independent Three axis between 3 error of perpendicularitys, 18 position correlation geometric errors, 3 position independence geometric errors altogether, such as Shown in table 1.
21 geometric errors of the Digit Control Machine Tool of table 1
Reference axis PDGEs PIGEs
X-axis δx(x),δy(x),δz(x),εx(x),εy(x),εz(x) --
Y-axis δx(y),δy(y),δz(y),εx(y),εy(y),εz(y) εxy
Z axis δx(z),δy(z),δz(z),εx(z),εy(z),εz(z) εxzyz
δ in tablex、δy、δzSite error is represented, subscript represents the direction of error;εx、εy、εzRepresent angular errors, inferior horn Mark represents the direction of angular errors rotation axis;X, y, z in bracket represents the direction of motion of X, Y, Z axis.εxy、εyz、εzxRespectively The error of perpendicularity between denotation coordination axle X and Y, Y and Z, X and Z." -- " represents no error.
Digit Control Machine Tool is a kind of multi-body system being in series by multiple rigid bodies, and theory of multi body system can be used to describe number The topological relation between each moving component on lathe is controlled, adjacent component on Digit Control Machine Tool is represented with homogeneous coordinate transformation matrix multiple Pose transformational relation between disjunctor coordinate system, so as to set up 21 geometric errors of Digit Control Machine Tool and lathe end position and attitude error Between mapping model (hereinafter referred to as machine tool error model).
Preferable static, the motion spy of adjacent body is obtained according to the topological structure of Fig. 2 u2000/800H horizontal Machining centers Matrix and preferable static, kinematic error eigenmatrix are levied, as shown in table 2.
The preferable static, kinematic matrix of the adjacent body of table 2 and preferable static, kinematic error matrix
T in tableijp、ΔTijp、Tijs、ΔTijsPreferable static feature matrix, the ideal movements feature of adjacent body are represented respectively Matrix, preferable Quiet Error eigenmatrix, ideal movements error character matrix;X, y, z represents machine tool motion coordinate respectively;xijp、 yijp、zijpMachine tool structure parameter is represented respectively;I4×4Represent quadravalence unit matrix;I, j=1,2,3 ..., 8;1 represents lathe bed, 2 tables Show X to slide unit, 3 represent B to turntable, and 4 represent workpiece, and 5 represent column, and 6 represent main spindle box, and 7 represent main shaft, and 8 represent cutter.
If homogeneous coordinates value P of the point of a knife point under main shaft coordinate systemtFor:
Pt=[0 0 t 1]T
T represents tool length in formula.
If homogeneous coordinates value P of the point to be processed under workpiece coordinate system on workpiecewFor:
Pw=[xw yw zw 1]T
xw,yw,zwRepresent coordinate value of the point to be processed under workpiece coordinate system on workpiece.
Homogeneous coordinates value P of the point of a knife point under preferable lathe bed coordinate systemt-idealFor:
Pt-ideal=T04pT04sT45pT45sT56pT67pPt
Homogeneous coordinates value P of the point to be processed under preferable lathe bed coordinate system on workpiecew-idealFor:
Pw-ideal=T01pT01sT12pT12sT23pPw
Under ideal conditions, point to be processed is to overlap on point of a knife point and workpiece, then has:
Pt-ideal=Pw-ideal
Homogeneous coordinates value P of the point of a knife point under actual lathe bed coordinate systemt-actualFor:
Pt-actual=T04pΔT04pT04sΔT04sT45pΔT45pT45sΔT45sT56pΔT56pT67pΔT67pPt
Homogeneous coordinates value P of the point to be processed under preferable lathe bed coordinate system on workpiecew-actualFor:
Pw-actual=T01pΔT01pT01sΔT01sT12pΔT12pT12sΔT12sT23pΔT23pPw
Space error is:
In formula:A、Respectively geometric error mapping matrix and geometric error vector.
By the geometric error of the related geometric error of above-mentioned lathe geometric error model opsition dependent and position independence write as Lower form:
In formulaRepresentation space error;AD、AIThe respectively related geometric error mapping matrix in position and position independence Geometric error mapping matrix;εD、εIThe respectively geometric error vector of the related geometric error vector sum position independence in position.
(2) interval theory is based on, carries out carrying out sensitivity analysis to each error source;
Interval number is defined as:
Superscript I expressions interval number, the superscript I implications occurred below are identical in formula;The upper and lower of section is represented respectively Limit.
By geometric error spatial model section:
Superscript T represents transposition in formula.
Interval extension factor ηi,jIt is defined as:
λ is used to describe the positive and negative of the interval extension factor in formula;The section of i-th of space error component is represented respectively Bound;Represent the section bound of j-th of geometric error component of a vector;Ai,jRepresent the of geometric error mapping matrix I rows, the element of jth row.
Sensitivity definition is:
For site error source, sensitivity definition:
For angular error source, due to the influence of Abbe error, sensitivity definition:
In formulaRepresent sensitivity when lathe is located at coordinate and tool length is;Represent lathe positioned at seat The interval extension factor of mark and tool length when being;Represent that lathe is located at the area when origin of coordinates and tool length are zero Between broadening factor.
(3) interval theory is based on, establishes range optimization distribution model:
It is about by design variable, complete machine space body diagonal error and each error source border of each error source interval width Beam condition, each error source sensitivity coefficient and the algebraical sum of interval width inverse square product are object function, and it is excellent to establish section Change distribution model.
Min represents to minimize in formula;F represents object function;XIRepresent the set of interval variable;N represents interval variable Sum;μiRepresent the sensitivity of i-th of interval variable;W represents interval width;Represent i-th of interval variable; Δ x, Δ y, Δ z represent X, Y, the space error of Z-direction respectively;[xI] represent interval variable constraint section.
(4) genetic algorithm (GA) is based on, is solved with the GAs Toolbox of MATLAB softwares so as to realize geometry Precision optimizing distributes;
Based on genetic algorithm, schematic diagram is as shown in figure 3, the genetic algorithm instrument of the MATLAB softwares with R2014a versions Case, genetic parameter use default value, carry out emulation solution, and operation result is as shown in figure 4, the optimization distribution of each error source geometric accuracy As a result it is as shown in table 3.
Each error source geometric accuracy optimization allocation result of table 3
Error source δx(x) δx(x) δx(x) δx(x) δx(x) δx(x) δx(x) δx(x) δx(x)
Empirical value/um 6 5 5 6 5 5 6 5 5
Optimal value/urad 8.9 7.4 7.4 8.9 7.4 7.4 8.9 7.4 7.4
Relax/% 48.3 48 48 48.3 48 48 48.3 48 48
Error source εx(x) εx(x) εx(x) εx(x) εx(x) εx(x) εx(x) εx(x) εx(x)
Empirical value/um 8 8 8 8 8 8 8 8 8
Optimal value/urad 11.5 12 11.7 11.9 10.7 9.5 11.2 11 11.9
Relax/% 43.8 50 46.3 48.8 33.8 18.8 40 37.5 48.8
Error source εxy εxy εxy
Empirical value/um 10 10 10
Optimal value/urad 14.9 11.1 12.8
Relax/% 49 11 28
Although the preferred embodiments of the present invention are described above in conjunction with accompanying drawing, the invention is not limited in upper The embodiment stated, above-mentioned embodiment is only schematical, be not it is restricted, this area it is common Technical staff in the case of present inventive concept and scope of the claimed protection is not departed from, may be used also under the enlightenment of the present invention By make it is many in the form of, these are belonged within protection scope of the present invention.

Claims (1)

  1. A kind of 1. precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory, it is characterised in that including Following steps:
    Step (1):Digit Control Machine Tool is regarded as to the multi-body system being in series by multiple rigid bodies first;Secondly transported using multi-body system Dynamic theory describes the topological relation between each moving component on Digit Control Machine Tool;Reuse homogeneous coordinate transformation matrix multiple table Show the pose transformational relation between adjacent component disjunctor coordinate system on Digit Control Machine Tool;Finally set up each geometric error of Digit Control Machine Tool With the mapping model between lathe space error;
    Step (2):Based on interval theory, sensitivity is defined according to the interval extension factor, line sensitivity of going forward side by side analysis is each to reflect Influence degree of the error source to complete machine space error;
    Step (3):It is constraints, respectively misses by design variable, complete machine space body diagonal error of each error source interval width Poor source sensitivity coefficient and the algebraical sum of interval width inverse square product are object function, establish range optimization distribution model;
    Step (4):Based on genetic algorithm, solved with MATLAB softwares so as to realize geometric accuracy optimization distribution.
CN201710544527.4A 2017-07-06 2017-07-06 Precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory Pending CN107368637A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710544527.4A CN107368637A (en) 2017-07-06 2017-07-06 Precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710544527.4A CN107368637A (en) 2017-07-06 2017-07-06 Precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory

Publications (1)

Publication Number Publication Date
CN107368637A true CN107368637A (en) 2017-11-21

Family

ID=60306070

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710544527.4A Pending CN107368637A (en) 2017-07-06 2017-07-06 Precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory

Country Status (1)

Country Link
CN (1) CN107368637A (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108052747A (en) * 2017-12-17 2018-05-18 北京工业大学 A kind of geometric precision of machine tool optimization method based on Method of valuo analysis
CN108873810A (en) * 2018-07-12 2018-11-23 沈阳机床股份有限公司 A kind of critical error source discrimination influencing the decay of three axis machining center precision
CN108920773A (en) * 2018-06-08 2018-11-30 华中科技大学 A kind of the ultraprecise kinematic system detailed protocol design method and system of kinetics-driven
CN110666590A (en) * 2019-09-12 2020-01-10 天津大学 Machine tool body diagonal error measuring method based on multi-beam laser interferometer
CN112558547A (en) * 2021-02-19 2021-03-26 成都飞机工业(集团)有限责任公司 Quick optimization method for geometric error compensation data of translational shaft of five-axis numerical control machine tool
CN115638754A (en) * 2022-10-03 2023-01-24 北京工业大学 Three-coordinate measuring machine precision distribution method based on inter-zone analytic method
CN117032072A (en) * 2023-08-28 2023-11-10 哈尔滨工业大学 Geometric error distribution method for five-axis ultra-precision machine tool

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6167634B1 (en) * 1998-03-28 2001-01-02 Snu Precision Co., Ltd. Measurement and compensation system for thermal errors in machine tools
CN103034166A (en) * 2012-11-26 2013-04-10 北京工业大学 Recognition method of critical geometrical error source of machine tool
CN104156519A (en) * 2014-07-30 2014-11-19 北京工业大学 Method for designing geometric accuracy of multi-axis numerical control machine tool to improve processing accuracy and reliability
CN105094047A (en) * 2015-07-31 2015-11-25 北京工业大学 Extended Fourier amplitude based extraction method of machine tool important geometric error source
CN105760662A (en) * 2016-02-03 2016-07-13 北京工业大学 Machine tool machining precision reliability and sensitivity analyzing method based on quick Markov chain

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6167634B1 (en) * 1998-03-28 2001-01-02 Snu Precision Co., Ltd. Measurement and compensation system for thermal errors in machine tools
CN103034166A (en) * 2012-11-26 2013-04-10 北京工业大学 Recognition method of critical geometrical error source of machine tool
CN104156519A (en) * 2014-07-30 2014-11-19 北京工业大学 Method for designing geometric accuracy of multi-axis numerical control machine tool to improve processing accuracy and reliability
CN105094047A (en) * 2015-07-31 2015-11-25 北京工业大学 Extended Fourier amplitude based extraction method of machine tool important geometric error source
CN105760662A (en) * 2016-02-03 2016-07-13 北京工业大学 Machine tool machining precision reliability and sensitivity analyzing method based on quick Markov chain

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
田文杰: "精密卧式加工中心几何精度设计及误差补偿方法研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108052747A (en) * 2017-12-17 2018-05-18 北京工业大学 A kind of geometric precision of machine tool optimization method based on Method of valuo analysis
CN108052747B (en) * 2017-12-17 2021-11-05 北京工业大学 Machine tool geometric precision optimization method based on value analysis method
CN108920773A (en) * 2018-06-08 2018-11-30 华中科技大学 A kind of the ultraprecise kinematic system detailed protocol design method and system of kinetics-driven
CN108920773B (en) * 2018-06-08 2020-06-02 华中科技大学 Design method and system for detailed scheme of dynamics-driven ultra-precise motion system
CN108873810A (en) * 2018-07-12 2018-11-23 沈阳机床股份有限公司 A kind of critical error source discrimination influencing the decay of three axis machining center precision
CN108873810B (en) * 2018-07-12 2021-03-12 沈阳机床股份有限公司 Key error source identification method influencing precision decay of triaxial machining center
CN110666590A (en) * 2019-09-12 2020-01-10 天津大学 Machine tool body diagonal error measuring method based on multi-beam laser interferometer
CN112558547A (en) * 2021-02-19 2021-03-26 成都飞机工业(集团)有限责任公司 Quick optimization method for geometric error compensation data of translational shaft of five-axis numerical control machine tool
CN112558547B (en) * 2021-02-19 2021-06-08 成都飞机工业(集团)有限责任公司 Quick optimization method for geometric error compensation data of translational shaft of five-axis numerical control machine tool
CN115638754A (en) * 2022-10-03 2023-01-24 北京工业大学 Three-coordinate measuring machine precision distribution method based on inter-zone analytic method
CN115638754B (en) * 2022-10-03 2024-03-15 北京工业大学 Three-coordinate measuring machine precision distribution method based on interval analytic hierarchy process
CN117032072A (en) * 2023-08-28 2023-11-10 哈尔滨工业大学 Geometric error distribution method for five-axis ultra-precision machine tool

Similar Documents

Publication Publication Date Title
CN107368637A (en) Precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory
CN110287553B (en) Machining error model global sensitivity analysis method based on Monte Carlo simulation
CN107966957B (en) A kind of numerical control worm gear-grinding machine space error decoupling compensation method
CN109765848B (en) Method for calculating geometric error contribution value of motion axis of five-axis machine tool and evaluating influence
CN105425727B (en) Five axis Flank machining cutter path method for fairing
CN108445839B (en) Machine tool precision sensitivity analysis method based on error increment
CN104865897B (en) The processing method of curved surface part and the process equipment of curved surface part
CN104460515B (en) One kind is based on the axle tool length compensation method of postpositive disposal five
CN108052747B (en) Machine tool geometric precision optimization method based on value analysis method
CN106502203A (en) A kind of Geometric Error for Computerized Numerical Control Milling Machine modeling method
CN103034166A (en) Recognition method of critical geometrical error source of machine tool
CN102621929B (en) Method for optimizing fixture height and machining path of double-rotary-table five-axis linkage numerical control machine tool
CN204584869U (en) The process equipment of curved surface part
CN106292545A (en) The method utilizing macroprogram digital control processing Cylinder Surface
CN107081638A (en) A kind of lathe work Space Thermal error identification method for establishing model
CN105243218A (en) Thermal error precision conversion and model establishing method of machine tool
CN106873522A (en) A kind of numerical control turning cutter path planing method of non-axis symmetry sweeping surface
CN109471408A (en) The three axis lathe geometric error compensation methods based on NC code refactoring
CN105458372A (en) Side milling error compensation device based on non-extended straight-line surface and cutter spacing planning method of side milling error compensation device
CN104536385A (en) Method for correcting machining program of numerical control machine tool
CN106378478A (en) Micro-milling machining method of regular pyramid micro structure
CN110989490A (en) Method for acquiring optimal installation position of workpiece based on contour error
CN106959665B (en) The geometric accuracy maintaining method of machine tool product multiple weighing value fusion based on big data
CN109839920A (en) A kind of five-axis machine tool kinematic axis Sensitivity Analysis Method
CN116909209A (en) Numerical control machine tool error modeling and predicting method considering dynamic thermal error

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20171121

WD01 Invention patent application deemed withdrawn after publication