CN115638754B - Three-coordinate measuring machine precision distribution method based on interval analytic hierarchy process - Google Patents
Three-coordinate measuring machine precision distribution method based on interval analytic hierarchy process Download PDFInfo
- Publication number
- CN115638754B CN115638754B CN202211217412.1A CN202211217412A CN115638754B CN 115638754 B CN115638754 B CN 115638754B CN 202211217412 A CN202211217412 A CN 202211217412A CN 115638754 B CN115638754 B CN 115638754B
- Authority
- CN
- China
- Prior art keywords
- cmm
- interval
- precision
- weight
- axis
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000008569 process Effects 0.000 title claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims abstract description 27
- 239000013598 vector Substances 0.000 claims abstract description 27
- 238000013461 design Methods 0.000 claims abstract description 17
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 5
- 238000012360 testing method Methods 0.000 claims abstract description 4
- 238000005259 measurement Methods 0.000 claims description 11
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 229910000831 Steel Inorganic materials 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 3
- 239000010959 steel Substances 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000008901 benefit Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
The invention discloses a precision distribution method of a three-coordinate measuring machine based on interval analytic hierarchy process, which comprises the steps of firstly establishing a CMM quasi-rigid body model; then solving 21 geometric errors under the precision requirement by using a LASO algorithm and a QR decomposition method in combination with the precision design requirement of the CMM whole machine; based on the mechanical structure of the CMM, establishing a precision hierarchy of the whole CMM; then constructing interval judgment matrixes of different parts of the CMM by using an interval analytic hierarchy process; consistency test is carried out on the constructed interval judgment matrix; solving interval weight vectors of 21 items of geometric errors relative to main parts of the CMM by using a computing characteristic root vector method; then, the interval weight is defuzzified, and the weight vector is normalized; and (3) carrying out top-down distribution calculation on the weight of the precision of the whole machine of the CMM through a vector calculation method, and finally obtaining the comprehensive weight of key parts of the CMM relative to the precision of the whole machine. The invention can effectively improve the design efficiency of the whole CMM and realize the aim of lowest cost.
Description
Technical Field
The invention relates to a precision distribution method of a three-coordinate measuring machine (Coordinate Measuring Machine, CMM for short), in particular to a precision distribution method of a three-coordinate measuring machine based on interval analytic hierarchy process.
Background
In the scheme design stage of the CMM product, the precision of main parts of the CMM product is designed according to the precision design requirement of the whole machine. And the precision of the key parts is reasonably allocated by comprehensively considering factors such as processing and assembling manufacturability, efficiency and the like by taking the precision of the whole machine as a target, and the precision of the parts is optimized by comprehensively considering factors such as cost on the basis. At present, the precision distribution of the CMM at home and abroad mainly comprises the following methods: one method is to perform precision allocation by establishing a geometric error model; another approach is to analogize to similar designs through existing CMM designs; the precision distribution design is also carried out by adopting the traditional modes of analogy, inquiry, experience estimation and the like. Because the method is difficult to efficiently finish the product precision distribution design and depends on experience and subjective preference of individuals, the balance between practical benefit and economic benefit of the CMM product is difficult to ensure.
Therefore, the precision distribution method of the three-coordinate measuring machine based on the interval analytic hierarchy process is necessary to be invented. The precision distribution of the main parts of the CMM is completed in the scheme design stage, so that the cost performance of the CMM product is improved.
Disclosure of Invention
The invention provides a method, which refers to a method for solving the geometric errors of a CMM (application number/patent number: CN 202010896543.1 'method for solving 21 items of geometric errors of a three-coordinate measuring machine based on LASSO algorithm'). The method has the characteristics of simplicity in operation, high efficiency, convenience and the like.
In order to achieve the purpose, the invention is realized by adopting the following technical scheme:
a three-coordinate measuring machine precision distribution method based on interval analytic hierarchy process includes the following steps:
step one: a quasi-rigid body model of the CMM is constructed.
Without complex deformations, the parts of the CMM can be regarded approximately as a rigid body. The CMM has 4 parts that move relative to each other: the device comprises a workbench, a bridge frame, a sliding frame and a main shaft. A rectangular coordinate system is established on the 4 parts which do relative motion: OXYZ, O 1 X 1 Y 1 Z 1 ，O 2 X 2 Y 2 Z 2 ，O 3 X 3 Y 3 Z 3 . The coordinate system is respectively connected with the movable bridge and the slide through a certain connecting rodThe frames are connected. The measuring head is fixedly arranged on the main shaft, and the measuring end has a coordinate (x) p0 ,y p0 ,z p0 ). When measuring the coordinates of a point on the workpiece on the CMM, the measurement is determined by the relative position of the measuring end and the workpiece. If the measurement end P0 is defined by the point (x p0 ，y p0 ，z p0 ) Moving X, Y, Z along the X, Y, Z axes, respectively, the position of the measuring tip in the ozz coordinate system becomes (X ', Y', Z '), and (X', Y ', Z') and (x+x), taking into account the existence of CMM geometrical errors p0 ,y+y p0 ,z+z p0 ) The difference is the measurement error. And according to the geometric errors and the motion relations of different motion axes, combining space coordinate transformation to obtain a relation model of the volume errors and 21 geometric errors.
Step two: and solving 21-term geometric errors of the CMM based on the quasi-rigid body model and the design accuracy requirement of the whole CMM.
The maximum allowable indication error of the CMM is MPE E = (a+b/1000) μm. The volume errors of the X axis, Y axis and Z axis of the CMM are all equal to the MPE E And consistent. Point a to be measured in CMM measurement space p (x p ,y p ,z p ) To the point A to be measured 1 (x 1 ,y 1 ,z 1 ) The displacement of the three coordinate axes is x respectively p1 ＝x p -x 1 ,y p1 ＝y p -y 1 ,z p1 ＝z p -z 1 Bringing the volume error and displacement of the to-be-measured point into a quasi-rigid body model of the CMM to obtain 21 geometric error solving equation sets:
A p x p ＝b p (1)
wherein:
b p ＝[Δx p ,Δy p ,Δz p ] T
the CMM21 geometry error is then solved using LASSO algorithm and QR decomposition.
Step three: and constructing a hierarchical structure of CMM precision by combining an interval hierarchical analysis method.
Step four: and constructing an interval judgment matrix of different parts of the CMM.
Number of scale intervals according to (1-9)Assigning values to the relative importance degrees among the elements, and constructing an interval judgment matrix of each element under the same criterion as follows:
wherein:i.e. < ->n represents n elements in total under the parent criterion k; />An arbitrary number of intervals in the range of 1-9 or the reciprocal thereof may be selected to represent the importance ratio between the ith element and the jth element under the parent criterion k; i is more than or equal to 1 and less than or equal to n, j is more than or equal to 1 and less than or equal to n; />
Wherein: the (1-9) scale concept is as follows:
a ij =1, meaning that element i and element j are of the same importance to the upper layer factor.
a ij =3, meaning that element i is slightly more important than element j.
a ij =5, representing an elementi is more important than element j.
a ij =7, meaning that element i is more important than element j.
a ij =9, meaning that element i is more important than element j.
a ij =2n, n=1, 2,3,4, meaning that element i and element j are of importance between a ij =2n_1 and a ij Between =2n+1.
Step five: and carrying out consistency test on the interval judgment matrix to judge the rationality of the interval judgment matrix.
When the following expression (3) is satisfied, it means that the consistency of the section judgment matrix is good, and it is possible to employ:
0≤α≤1≤β (3)
wherein:
step six: the interval weight vector of 21 items of geometric errors relative to the main parts of the CMM is calculated.
Adopting a characteristic root vector method to respectively calculate A k- 、A k+ And record the weight vectors of (a) as x respectively k- And x k+ 。
Then calculate A k The interval weight of (2) is:
W k ＝(w) 1×n ＝[αx k- ,βx k+ ]＝[W k- ,W k+ ] (4)
respectively taking weight vectors W k- ,W k+ As a local interval weight vector element.
Step seven: the interval weight defuzzifies.
And (3) taking the difference value of the upper limit and the lower limit of the weight interval into consideration, and performing defuzzification on the interval weight, wherein the formula 5 is a weight interval defuzzification formula.
Order theThen
ω i ＝M i +(2γ i -1)·D i (5)
Wherein: m is M i Represents the median value of the weight interval, D i Representing half of the difference between the left and right intervals of the weight, wherein gamma is more than or equal to 0 i ≤1,1≤i≤n，γ i Representing the subjective preferences of the designer for the accuracy design requirements. And carrying out normalization processing on the weight vector to obtain the comprehensive weight of the CMM key parts relative to the precision of the whole machine. The weight vector is:
W i ＝[ω 1 ,ω 2 ,...,ω n ] (6)
wherein: w (W) i Represents a weight vector, ω 1 -ω n Representing the weights of the different components.
Step eight: and calculating the distribution precision of the key parts of the CMM.
In summary, the method provided by the invention can rapidly obtain the distribution precision of the key parts of the CMM, effectively improve the design efficiency of the whole machine and realize the aim of lowest cost.
Drawings
FIG. 1 is a diagram of a CMM model of the structure type YFXZ.
Fig. 2a is a graph of X-axis positioning error and straightness motion error.
Fig. 2b is an X-axis angular motion error plot.
Fig. 2c is a graph of Y-axis positioning error and straightness motion error.
Fig. 2d is a Y-axis angular motion error plot.
Fig. 2e is a graph of Z-axis positioning error and straightness motion error.
Fig. 2f is a Z-axis angular motion error plot.
Figure 3 is a hierarchical diagram of the accuracy of a CMM.
In the figure: 1. the device comprises a base, 2, a workbench, 3, a bridge frame, 4, a sliding frame, 5 and a Z axis.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be further described in detail below by referring to the accompanying drawings and examples, so that those skilled in the art can implement the present invention by referring to the text of the specification.
The experiment adopts the structure type shown in figure 1 asThe accuracy requirement of the YFXZ-type CMM is as follows: MPE (MPE) E =2 μm+l/180 μm. A three-coordinate measuring machine precision distribution method based on interval analytic hierarchy process comprises the following steps:
step one: a quasi-rigid body model of the CMM is constructed.
The YFXZ-type CMM quasi-rigid body model is as follows:
wherein delta x (x) Positioning error for X axis; delta y (x) Is the straightness error of the X axis in the Y direction; delta z (x) Is the straightness error of the X axis in the Z direction; epsilon x (x) Is X-axis rolling error; epsilon y (x) Pitch error for the X-axis; epsilon z (x) Is X-axis deflection error; delta y (Y) is Y-axis positioning error; delta x (Y) is the Y-axis X-direction straightness error; delta z (Y) is the Y-axis Z-direction straightness error; epsilon x (Y) is a Y-axis roll pitch error; epsilon y (Y) is a Y-axis roll error; epsilon z (Y) is the Y-axis yaw error; delta z (Z) is a Z-axis positioning error; delta x (Z) is Z-axis X-direction straightness error; delta y (Z) is the Z-axis Y-direction straightness error; epsilon x (Z) is a Z-axis pitch error; epsilon y (Z) is a Z-axis pitch error; epsilon z (Z) is a Z-axis roll error; alpha xy Is the perpendicularity error of the X axis and the Y axis; alpha xz Is the perpendicularity error of the X axis and the Z axis; alpha yz Is the perpendicularity error of the Y axis and the Z axis.
Step two: and solving 21-term geometric errors of the CMM based on the quasi-rigid body model and the design accuracy requirement of the whole CMM.
Maximum allowable error of CMMThe difference is MPE E Volume errors of X-axis, Y-axis and Z-axis of CMM are all equal to MPE =2 μm+l/180 μm E And consistent. Point a to be measured in CMM measurement space i (x i ,y i ,z i ) To the point A to be measured 1 (x 1 ,y 1 ,z 1 ) The displacement of the three coordinate axes is x respectively i1 ＝x i -x 1 ,y i1 ＝y i -y 1 ,z i1 ＝z i -z 1 And (3) bringing the volume error and displacement of the to-be-measured point into quasi-rigid body models (1) - (3) of the CMM to obtain the following steps:
A i x i ＝b i (10)
wherein:
b i ＝[Δx i ,Δy i ,Δz i ] T (13)
let the number of geometric errors contained in 21 geometric errors be f, n measurement points in total, and arrange the equation set to obtain:
because the coefficient matrix a of the equation set (14) is singular, the equation set cannot be solved by the conventional least square method. And solving the quasi-rigid body model formulas (1) to (3) by using an LASSO algorithm in machine learning. Because epsilon z (x),ε x (z),ε y (z),ε z The coefficients of the four angular motion errors (z) are composed of the coordinates of the initial measurement points, the initial points being the coordinatesThe origin, the coefficient of the four angular motion errors is zero. So 17 geometric errors can be found using LASSO.
Geometric model based on volume error and uniaxial geometric error, and angular motion error epsilon can be obtained by utilizing QR decomposition z (x),ε x (z),ε y (z),ε z (z)。
The 21-term geometry error results are shown in fig. 2.
Step three: and constructing a hierarchical structure of CMM precision by combining an interval hierarchical analysis method.
Layering the whole machine precision of the CMM by using a layer-by-layer distribution idea: the uppermost layer is the overall design constraint of the accuracy of the CMM machine, i.e. MPE of the CMM E The method comprises the steps of carrying out a first treatment on the surface of the The second layer divides the CMM into X, Y, Z three motion modules according to the motion axis; the third layer is a main component layer of the CMM and comprises a CMM base, a cross beam, a Z axis, a precision ball screw, a grating ruler, a steel wire toothed belt and the like; the fourth layer is the 21 term geometry error of the CMM. The constructed precision hierarchy is shown in fig. 3.
Step four: and constructing an interval judgment matrix of different parts of the CMM.
With ball screw C GX As a father criterion, the importance weight interval judgment matrix between the corresponding CMM geometric errors is as follows:
since the mapping relation of the corresponding elements (geometric errors of the fourth layer) of the guide rails of the three motion axes of the CMM is the same, the guide rail C of the X, Y, Z axis is D ,C H ,C L The importance interval judgment matrix of the corresponding geometric errors is the same. Taking the guide rail as a father criterion, judging a matrix C in an importance weight interval between corresponding CMM geometric errors D The method comprises the following steps:
step five: and carrying out consistency test on the interval judgment matrix to judge the rationality of the interval judgment matrix.
And (3) calculating a section judgment matrix consistency parameter by using the formula, wherein under the ball screw rule, alpha= 0.9537 and beta= 1.0382, and the constructed section judgment matrix meets the consistency requirement.
Step six: the interval weight vector of 21 items of geometric errors relative to the main parts of the CMM is calculated.
1) Matrix C GX Is decomposed into two matrices of which the matrix is composed,respectively called a left judgment matrix and a right judgment matrix;
2) Respectively writing out judgment matrixAnd then calculating the maximum eigenvalue and the corresponding maximum eigenvector of the left and right judgment matrixes respectively by utilizing a method for solving matrix eigenvectors, and carrying out normalization processing to obtain the corresponding weight vector.
3) Combining the interval number weight calculation formula: w (W) k ＝(w) 1×n ＝[αx k- ,βx k+ ]＝[W k- ,W k+ ]Calculating to obtain C GX A corresponding weight vector. The weight vectors of the ball screw of the Y motion axis and the Z axis steel wire toothed belt can be obtained by the same way as the three motion axis transmission mechanisms have the same effect.
Since the three motion axis guide rails have the same function, the importance interval judgment matrix corresponding to the geometric errors is the same.
Step seven: the interval weight is defuzzified, and the obtained weight values of the precision distribution of the ball screw and the guide rail are shown in tables 1 and 2.
Table 1 weights of ball screw corresponding to 6 geometric errors
Table 2 weights of the corresponding geometric errors of the shaft guide
Step eight: the distribution accuracy of the key parts of the CMM is calculated, and the results are shown in table 3.
Table 3 CMM key component precision assignment results
The above examples are merely illustrative of preferred embodiments of the present invention and are not intended to limit the scope of the present invention, as any changes, modifications and improvements made within the principles of the present invention shall fall within the scope of the invention.
Claims (5)
1. The three-coordinate measuring machine precision distribution method based on the interval analytic hierarchy process is characterized by comprising the following steps of: the method comprises the following steps:
step one: constructing a quasi-rigid body model of the CMM;
the components of the CMM are considered rigid bodies; the CMM has four parts that move relative to each other: the device comprises a workbench, a bridge frame, a sliding frame and a main shaft; a rectangular coordinate system is established on the four parts which do relative motion: OXYZ, O 1 X 1 Y 1 Z 1 ，O 2 X 2 Y 2 Z 2 ，O 3 X 3 Y 3 Z 3 The method comprises the steps of carrying out a first treatment on the surface of the The measuring head is fixedly arranged on the main shaft, and the measuring end has coordinates (x) p0 ,y p0 ,z p0 ) The method comprises the steps of carrying out a first treatment on the surface of the When the coordinate of a certain point of the workpiece is measured on the CMM, the measurement result is determined by the relative position of the measuring end and the workpiece; if the measuring end P 0 From the point (x p0 ，y p0 ，z p0 ) Moving x, y, z along X, Y, Z axes respectively, the position of the measuring tip in the xyz coordinate system becomes (x ', y', z '), (x', y ', z') and (x+x) taking into account the existence of CMM geometrical errors p0 ,y+y p0 ,z+z p0 ) The difference is the measurement error; according to the geometric errors and the motion relations of different motion axes, combining space coordinate conversion to obtain a relation model of the volume errors and 21 geometric errors;
step two: solving a relation model of 21 items of geometric errors of the CMM based on a quasi-rigid body model and the design precision requirement of the whole CMM;
the maximum allowable indication error of the CMM is MPE E Volume errors of X-axis, Y-axis and Z-axis of CMM are all equal to MPE for = (a+b/1000) μm E Consistent; point a to be measured in CMM measurement space p (x p ,y p ,z p ) To the point A to be measured 1 (x 1 ,y 1 ,z 1 ) The displacement of the three coordinate axes is x respectively p1 ＝x p -x 1 ,y p1 ＝y p -y 1 ,z p1 ＝z p -z 1 Bringing the volume error and displacement of the to-be-measured point into a quasi-rigid body model of the CMM to obtain 21 geometric error solving equation sets:
A p x p ＝b p (1)
wherein:
b p ＝[Δx p ,Δy p ,Δz p ] T
then solving the geometric error of the CMM21 by using an LASSO algorithm and QR decomposition;
step three: constructing a hierarchical structure of CMM precision by combining an interval hierarchical analysis method; layering the whole machine precision of the CMM by using a layer-by-layer distribution idea: the uppermost layer is the overall design constraint of the accuracy of the CMM machine, i.e. MPE of the CMM E The method comprises the steps of carrying out a first treatment on the surface of the The second layer divides the CMM into X, Y, Z three motion modules according to the motion axis; the third layer is a CMM component layer and comprises a CMM base, a cross beam, a Z axis, a precision ball screw, a grating ruler and a steel wire toothed belt; the fourth layer is the 21 term geometry error of the CMM;
step four: constructing an interval judgment matrix of different parts of the CMM;
according to the number of 1-9 scale intervalsAssigning values to the relative importance degrees among the elements, and constructing an interval judgment matrix of each element under the same criterion as follows:
wherein:i.e. < ->n represents n elements in total under the parent criterion k; />Representing the importance ratio between the ith element and the jth element under the father criterion k, and selecting any interval number in the range of 1-9 or the reciprocal thereof; i is more than or equal to 1 and less than or equal to n, j is more than or equal to 1 and less than or equal to n; />
Wherein: the 1-9 scale concept is as follows:
a ij =1, meaning that element i and element j have the same importance to the upper layer factor;
a ij =3, meaning that element i is slightly more important than element j;
a ij =5, indicating that element i is more important than element j;
a ij =7, meaning that element i is more important than element j;
a ij =9, meaning that element i is more important than element j;
a ij =2n, n=1, 2,3,4, meaning that element i and element j are of importance between a ij =2n_1 and a ij Between =2n+1;
step five: consistency test is carried out on the interval judgment matrix, and the rationality is judged;
when the following formula (3) is satisfied, it means that the consistency of the section judgment matrix is good, and the following is adopted:
0≤α≤1≤β (3)
wherein:
step six: calculating interval weight vectors of 21 geometric errors relative to the CMM parts;
adopting a characteristic root vector method to respectively calculate A k- 、A k+ And record the weight vectors of (a) as x respectively k- And x k+ ；
Then calculate A k The interval weight of (2) is:
W k ＝(w) 1×n ＝[αx k- ,βx k+ ]＝[W k- ,W k+ ] (4)
respectively taking weight vectors W k- ,W k+ As a local interval weight vector element;
step seven: disambiguation of interval weights;
the upper limit and the lower limit of the weight interval are considered to carry out defuzzification on the interval weight, and the formula (5) is a weight interval defuzzification formula;
order theThen
ω i ＝M i +(2γ i -1)·D i (5)
Wherein: m is M i Represents the median value of the weight interval, D i Representing half of the difference between the left and right intervals of the weight, wherein gamma is more than or equal to 0 i ≤1,1≤i≤n，γ i Representing subjective preferences of designers for accuracy design requirements; normalizing the weight vector to obtain the comprehensive weight of the key parts of the CMM relative to the precision of the whole machine; the weight vector is:
W i ＝[ω 1 ,ω 2 ,...,ω n ] (6)
wherein: w (W) i Represents a weight vector, ω 1 -ω n The weights of different parts are represented;
step eight: calculating the distribution precision of key parts of the CMM; and (3) providing the precision requirements of different parts by using the 21 geometric errors obtained in the second step and the weight vector obtained in the seventh step in a vector multiplication summation mode.
2. The method for distributing precision of three-coordinate measuring machine based on interval analytic hierarchy process as claimed in claim 1, wherein: and (3) constructing a quasi-rigid body model of the YFXZ type three-coordinate measuring machine, namely, the movement sequence is Y axis-stylus-X axis-Z axis.
3. The method for distributing precision of three-coordinate measuring machine based on interval analytic hierarchy process as claimed in claim 1, wherein: in the second step, on the basis of the constructed quasi-rigid body model, the maximum allowable indication error of the CMM is added to an X axis, a Y axis and a Z axis and used as a volume error, volume error data of a planning point in a measurement space is generated, and a relation model of 21-term geometric errors of the CMM is solved by using an LASO algorithm and QR decomposition.
4. The method for distributing precision of three-coordinate measuring machine based on interval analytic hierarchy process as claimed in claim 1, wherein:
in the fourth step, according to the size and experience of the 21 geometric errors obtained by solving in the second step, an interval judgment matrix of different parts of the CMM is constructed.
5. The method for distributing precision of three-coordinate measuring machine based on interval analytic hierarchy process as claimed in claim 1, wherein: in the seventh step of the method, after a weight interval allocated by the CMM precision is obtained, the interval weight is defuzzified, and the weight vector is normalized to obtain reasonable weights of corresponding elements of different parts.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211217412.1A CN115638754B (en) | 2022-10-03 | 2022-10-03 | Three-coordinate measuring machine precision distribution method based on interval analytic hierarchy process |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211217412.1A CN115638754B (en) | 2022-10-03 | 2022-10-03 | Three-coordinate measuring machine precision distribution method based on interval analytic hierarchy process |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115638754A CN115638754A (en) | 2023-01-24 |
| CN115638754B true CN115638754B (en) | 2024-03-15 |
Family
ID=84942526
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211217412.1A Active CN115638754B (en) | 2022-10-03 | 2022-10-03 | Three-coordinate measuring machine precision distribution method based on interval analytic hierarchy process |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115638754B (en) |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107368637A (en) * | 2017-07-06 | 2017-11-21 | 天津大学 | Precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory |
| CN109522643A (en) * | 2018-11-14 | 2019-03-26 | 北京工业大学 | Precision of NC Machine Tool distributes Multipurpose Optimal Method |
| CN110209119A (en) * | 2019-06-26 | 2019-09-06 | 重庆大学 | A kind of Precision of NC Machine Tool appraisal procedure and life-span prediction method based on metaaction unit and the subjective and objective weight of integration |
| CN110345867A (en) * | 2019-07-15 | 2019-10-18 | 北京工业大学 | A method of the space CMM Point Coordinates error is corrected based on inverse distance weighting |
| CN110597177A (en) * | 2019-09-18 | 2019-12-20 | 重庆大学 | Numerical control machine tool precision control method based on precision mapping |
| CN111144023A (en) * | 2019-12-31 | 2020-05-12 | 华能如东八仙角海上风力发电有限责任公司 | Reliability optimization distribution method based on improved Pareto artificial bee colony algorithm |
| CN112229321A (en) * | 2020-08-31 | 2021-01-15 | 北京工业大学 | Method for solving 21-item geometric errors of three-coordinate measuring machine based on LASSO algorithm |
| CN113932707A (en) * | 2021-08-29 | 2022-01-14 | 北京工业大学 | Method for compensating geometric error of three-coordinate measuring machine based on elastic network algorithm |
| CN114034244A (en) * | 2021-11-08 | 2022-02-11 | 北京工业大学 | Method for obtaining geometric error of any point in CMM measuring range based on Makima algorithm |
| CN114757572A (en) * | 2022-05-10 | 2022-07-15 | 重庆大学 | Method for determining key quality characteristics of electromechanical product by considering uncertainty |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003333355A (en) * | 2002-05-10 | 2003-11-21 | Canon Inc | Color evaluation apparatus and method |
| JP6295070B2 (en) * | 2013-12-05 | 2018-03-14 | オークマ株式会社 | Geometric error identification method for multi-axis machine tools and multi-axis machine tools |
| CN116091734A (en) * | 2016-09-29 | 2023-05-09 | 北京理工大学 | Geometric error characterization method based on NURBS curved surface reconstruction |
-
2022
- 2022-10-03 CN CN202211217412.1A patent/CN115638754B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107368637A (en) * | 2017-07-06 | 2017-11-21 | 天津大学 | Precise horizontal machining center geometric accuracy optimizing distribution method based on interval theory |
| CN109522643A (en) * | 2018-11-14 | 2019-03-26 | 北京工业大学 | Precision of NC Machine Tool distributes Multipurpose Optimal Method |
| CN110209119A (en) * | 2019-06-26 | 2019-09-06 | 重庆大学 | A kind of Precision of NC Machine Tool appraisal procedure and life-span prediction method based on metaaction unit and the subjective and objective weight of integration |
| CN110345867A (en) * | 2019-07-15 | 2019-10-18 | 北京工业大学 | A method of the space CMM Point Coordinates error is corrected based on inverse distance weighting |
| CN110597177A (en) * | 2019-09-18 | 2019-12-20 | 重庆大学 | Numerical control machine tool precision control method based on precision mapping |
| CN111144023A (en) * | 2019-12-31 | 2020-05-12 | 华能如东八仙角海上风力发电有限责任公司 | Reliability optimization distribution method based on improved Pareto artificial bee colony algorithm |
| CN112229321A (en) * | 2020-08-31 | 2021-01-15 | 北京工业大学 | Method for solving 21-item geometric errors of three-coordinate measuring machine based on LASSO algorithm |
| CN113932707A (en) * | 2021-08-29 | 2022-01-14 | 北京工业大学 | Method for compensating geometric error of three-coordinate measuring machine based on elastic network algorithm |
| CN114034244A (en) * | 2021-11-08 | 2022-02-11 | 北京工业大学 | Method for obtaining geometric error of any point in CMM measuring range based on Makima algorithm |
| CN114757572A (en) * | 2022-05-10 | 2022-07-15 | 重庆大学 | Method for determining key quality characteristics of electromechanical product by considering uncertainty |
Non-Patent Citations (2)
| Title |
|---|
| Calibration Method for Angular Positioning Deviation of a High-Precision Rotary Table Based on the Laser Tracer Multi-Station Measurement System;Hongfang Chen等;《applied sciences》;全文 * |
| 基于区间层次分析法的高速送料机可靠性综合分配方法;朱庆鹏等;《锻压技术》;第40卷(第4期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115638754A (en) | 2023-01-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103390082B (en) | The sane excellent method of completing the square of a kind of gang tool geometric accuracy | |
| CN104375460B (en) | A kind of Digit Control Machine Tool machining accuracy reliability sensitivity analysis method | |
| Zhang et al. | Measurement and compensation for volumetric positioning errors of CNC machine tools considering thermal effect | |
| Wu et al. | A geometric accuracy analysis and tolerance robust design approach for a vertical machining center based on the reliability theory | |
| CN112697041B (en) | Monte carlo method-based pre-evaluation method for measurement accuracy of assembly pose | |
| Ma et al. | Assembly errors analysis of linear axis of CNC machine tool considering component deformation | |
| Li et al. | A sensitivity method to analyze the volumetric error of five-axis machine tool | |
| CN104200063B (en) | The uncertainty description of lathe Space processing error and Forecasting Methodology | |
| Baum et al. | Thermally induced volumetric error compensation by means of integral deformation sensors | |
| Majda | Relation between kinematic straightness errors and angular errors of machine tool | |
| CN113910001B (en) | Numerical control machine tool space error identification method | |
| CN115638754B (en) | Three-coordinate measuring machine precision distribution method based on interval analytic hierarchy process | |
| CN113656902A (en) | Error sensitivity analysis method for working end pose of multi-axis motion platform | |
| Sun et al. | Analysis and optimization of assembly precision-cost model based on 3D tolerance expression | |
| CN109883381A (en) | A kind of three-dimensional space large-scale metrology method of articulated coordinate machine | |
| Fan et al. | A novel global sensitivity analysis method for vital geometric errors of five-axis machine tools | |
| Li et al. | Sculptured surface-oriented machining error synthesis modeling for five-axis machine tool accuracy design optimization | |
| CN113359609A (en) | Key geometric error optimization proportioning compensation method for five-axis numerical control machine tool | |
| Zou et al. | Error distribution of a 5-axis measuring machine based on sensitivity analysis of geometric errors | |
| Sun et al. | The grouping and matching of three dimensional deviation based on homogeneous transformation and Taguchi theory | |
| CN113932707B (en) | Method for compensating geometric error of three-coordinate measuring machine based on elastic network algorithm | |
| Li et al. | A novel error equivalence model on the kinematic error of the linear axis of high-end machine tool | |
| Brecher et al. | Hybrid modeling of thermo-elastic behavior of a three-axis machining center using integral deformation sensors | |
| CN117921685B (en) | Calibration method for co-positioning precision of double robots | |
| CN109933950A (en) | A kind of guideway abrasion analysis and prediction technique based on multiple dimensioned bridge domain method |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |