CN113932754A - Method for quickly evaluating verticality of error measurement workpiece based on single reference plane - Google Patents

Method for quickly evaluating verticality of error measurement workpiece based on single reference plane Download PDF

Info

Publication number
CN113932754A
CN113932754A CN202010602658.5A CN202010602658A CN113932754A CN 113932754 A CN113932754 A CN 113932754A CN 202010602658 A CN202010602658 A CN 202010602658A CN 113932754 A CN113932754 A CN 113932754A
Authority
CN
China
Prior art keywords
plane
point
measuring point
key
measured
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
CN202010602658.5A
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.)
Guilin University of Electronic Technology
Original Assignee
Guilin University of Electronic 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 Guilin University of Electronic Technology filed Critical Guilin University of Electronic Technology
Priority to CN202010602658.5A priority Critical patent/CN113932754A/en
Publication of CN113932754A publication Critical patent/CN113932754A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/22Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring angles or tapers; for testing the alignment of axes

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

The invention belongs to the field of precision measurement and computer application, and relates to a method for evaluating hole shaft parts by combining a virtual gauge with measurement data in a mathematical model mode, in particular to a method for evaluating the perpendicularity of an L-shaped error measurement workpiece based on a single reference plane, which comprises the following steps: 1: acquiring a reference plane measuring point set and a measured element point set, and establishing a characteristic row vector set and a state element set; 2: obtaining an initial key sequence number set; 3: establishing an analysis matrix and an analysis column vector; 4: performing rank analysis; 5: determining an optimizing direction; 6: solving new key points and updating coordinates of the measuring point set; 7: judging whether the evaluation of the reference elements meets the qualification conditions; 8: updating the point set of the measured elements; 9: updating the real-time state element set; 10: updating the key sequence number set of the measured elements; 11: establishing an analysis matrix and an analysis column vector; 12: performing rank analysis; 13: determining an optimizing direction; 14: solving new key points and updating coordinates of the measuring point set; 15: and (5) judging the qualification.

Description

Method for quickly evaluating verticality of error measurement workpiece based on single reference plane
Technical Field
The invention belongs to the field of precision metering and computer application, and particularly relates to a method for evaluating the verticality of an error measurement workpiece based on a single reference plane, which is stable, rapid and simple in form, can be used for the qualification detection and error evaluation of the verticality of an L-shaped measurement workpiece of the single reference plane, and the evaluation result can provide useful guidance for the improvement of a machining process.
Background
The invention takes an L-shaped error measurement workpiece as an example, and verifies a single-datum-plane-based quick evaluation method for perpendicularity of the L-shaped error measurement workpiece. The specific placing position of the workpiece is that the bottom surface of the L-shaped workpiece is attached to a horizontal plane, and the side surface of the L-shaped workpiece is vertical to the bottom surface. The bottom surface of the L-shaped workpiece is used as a reference surface, and the vertical surface is a measured plane.
The size error and the shape and position error (short for shape error and position error) directly influence the product quality, the assembly and the service life of the product, and the method has important significance for quickly and accurately calculating the part error. The perpendicularity error of the angle iron part based on the single reference plane refers to the fact that a theoretical correct angle is formed between the angle iron part and the reference plane, and an extraction (actual) surface is limited between two parallel planes which are perpendicular to the reference plane and have a tolerance value h. However, the four types of currently popular evaluation methods are difficult to directly apply to the quick evaluation of the perpendicularity of the angle iron parts with a single reference plane.
The first, the optical gap method. By comparing the right-angle measuring tool with the measured element, whether the part meets the verticality error requirement can be directly measured through the optical gap (or the feeler gauge). During detection, one surface of the square is placed on a reference surface, so that the two surfaces are tightly attached, and the square can be adjusted to a position meeting the minimum condition when necessary. And then, the vertical edge of the square is contacted with the measured actual surface, and the maximum gap f is formed between the vertical edge of the square and the measured actual surface, namely the perpendicularity error of the part. The gap value can be measured directly by the optical gap method or with a feeler gauge. The variation between the measured actual element and the square can be measured by using the measuring blocks, namely two groups of measuring blocks with different sizes are respectively padded at the maximum and minimum gaps between the measured actual surface and the square. And adjusting the thicknesses of the two groups of gauge blocks until the gaps at the gauge blocks are eliminated, wherein the size difference of the two groups of gauge blocks is the perpendicularity error of the part. The measurement method mainly depends on manual measurement, has low measurement efficiency and is not suitable for measurement in mass part production.
The second category, indicator methods. And (3) taking two adjacent vertical surfaces of the square box as ideal simulation elements, and comparing the ideal simulation elements with the measured elements through an indicator to measure the verticality error. According to the perpendicularity tolerance of the upper surface to the left side surface of the given part. During detection, the reference surface of the detected part is fixed on the right-angle seat, and the difference of the indicating and measuring values close to the reference detected surface is adjusted to be minimum. And then sliding the whole surface for measurement, and taking the difference between the maximum and minimum indication values in each measurement point as the verticality error of the part. If necessary, the verticality error can be evaluated by the minimum directional area. The indicator used in the method has errors and has no obvious advantage compared with the first method.
And in the third category, the horizontal plane is used as a reference, two adjacent vertical surfaces of the frame type level meter are used as ideal simulation elements, and the perpendicularity error is measured. According to the perpendicularity tolerance of the vertical guide rail to the plane of the workbench. During detection, a frame type level meter is used for roughly adjusting the reference surface to a horizontal position, then the frame type level meter is used for sequentially measuring the reference actual surface according to a certain point distribution distance, the vertical side surface of the frame type level meter is used for measuring the measured actual surface, and indication values of all measuring points are recorded. And converting the measured angle value into a line value, thereby solving the perpendicularity error of the measured part. The method can not solve the problem of low measurement efficiency, is not simple and convenient to operate, has no obvious advantages compared with the first two methods, and is mainly suitable for measuring large parts.
The fourth category, artificial intelligence/biological intelligence algorithms. The advantage of this type of method over the third type of method is to analyze the "objective function with complex gradient or no apparent analytic expression" and to find the "global optimum". The method also realizes standard solutions in a plurality of software at present, so the method is easy to popularize. Although these methods are relatively hot at present, they are not suitable for use in the diameter evaluation of the smallest circumscribed cylinder. This is because the gradient of the objective function for the diameter assessment of the minimum circumscribing cylinder is the sum of a large number of simple analytical expressions, and part errors are typically small. Therefore, the "local optimal value" of the objective function can be considered as the "global optimal value", and the fourth method has no obvious advantages over the third method.
And a fifth category, specialized geometric assessment methods. And gradually searching for perpendicularity errors meeting the definition according to the geometric properties between the reference elements and the measured elements and the national standard and ISO standard of the angle iron parts. In the method, the form of a mathematical model is complex, the speed change is uneven, and the method is not easy to be popularized to various tolerance evaluations, and certainly, the method is also difficult to be popularized to the verticality evaluation of the angle iron parts based on the single reference plane.
In conclusion, when the existing geometric evaluation method is applied to the perpendicularity of the angle iron parts with a single reference plane, stability, quickness and simple form cannot be simultaneously considered.
Disclosure of Invention
The purpose of the invention is:
aiming at the problems in the prior art, the invention provides a stable, quick and simple-form quick evaluation method for L-shaped measurement workpiece perpendicularity errors based on a single reference plane, and the evaluation result can also provide useful guidance for subsequent improved machining processes.
The scheme adopted by the invention is as follows:
according to the method, on the basis that the evaluation of the reference plane is completed and qualified, the verticality error of the angle iron part based on the single reference plane is evaluated.
A method for quickly evaluating the verticality of a measured workpiece based on a single reference plane is mainly realized by the following steps:
step 1: firstly, measuring points of a reference surface A are obtained and are used to form a reference measuring point setu j And according to au j Establishing a characteristic line vector setW j Greatb j }; obtaining initial measuring point of measured element and using it to form measuring point setu i }; wherein:
j=1, 2, 3, …, Mjthe serial number of the reference measuring point is,Mthe total number of the measuring points of the reference plane is;
i=1, 2, 3, …, Nithe serial number of the measuring point of the measured segment,Nthe total number of the measuring points of the measured plane is;
u j ={x j , y j , z j is the measurement pointjAnd the Z-axis of the coordinate system is perpendicular to the reference plane, the XOY plane of the coordinate systemClose to the ideal reference plane;
t j =z j all state elementst j Is a set of state elementst j };
When z is j >At the time of 0, the number of the first,W j =([-1, -y j , x j ]) (ii) a When z is j <At the time of 0, the number of the first,W j =([1, y j ,- x j ]);W j is a feature row vector, all feature row vectorsW j Is a set of characteristic line vectorsW j };
b j =bIs a real number greater than 0, all boundary elementsb j Is a set of boundary elementsb j };
After the step 1 is finished, performing a step 2;
step 2: gett j Maximum valuet max Corresponding serial numberr 2 Is numbered for the key point and willr 2 Last page added to key serial number setrIn (c) } the reaction solution is,t max the distance from the reference measuring point to the upper boundary; gett j Minimum valuet min Corresponding serial numberr 1 Is numbered for the key point and willr 1 Last page added to key serial number setrIn (c) } the reaction solution is,t min is the distance from the reference point to the lower boundary.
Step 3 is carried out after step 2 is finished;
and step 3: according to the key sequence numberrEstablishment of an analysis matrixWAnd analyzing the column vectorsbWherein:
W=[…, A p T, …, A q T, …]Tis aFA matrix of rows and 4 columns,Fis a set of key pointsrThe number of the elements in the (C),p, qis a set of key pointsrThe elements in (1);
b=[…, b p , …, b q , …]Tis aFA column vector of rows.
Step 4 is carried out after step 3 is finished;
and 4, step 4: for analysis matrixWAnd an augmented analysis matrixW, b]Performing rank analysis;
computing an analysis matrixWRank ofR W =rank(W) Extended analysis matrix [ alpha ], [ beta ] anW, b]Rank ofR Wb =rank([W, b]) And compareR W AndR Wb there are only two cases:
the first condition is as follows: if it is notR W =R Wb Then, the optimization should be continued, jumping to step 5;
case two: if it is notR W <R Wb Then, the counter C =1 is initialized and step 4.1 is performed;
step 4.1: attempting to slave an analysis matrixWAnd analyzing the column vectorsbMiddle removed key point setrGet the line corresponding to one element in the matrix to get the reduced matrixW r - And reducing the column vectorb r - (ii) a WhereinW r - Is oneJA matrix of rows and 4 columns,b r - is oneJA column vector of rows.
Step 4.2: solving linear equationsW r - v r - = b r - Solution of (2)v r - =v r-0 Then calculateb r =W r v r-0 (ii) a If it is notb r - > b r Then, the matrix will be reducedW r - And reducing the column vectorb r - Respectively asWMatrix and scoreDisjunctive column vectorbWill elementrMoving out key point setrAnd jumping to the step 5; if it is notb r - <= b r Then step 4.3 is performed;
step 4.3: judging counterCWhether or not equal toJIf, ifCJJumping to step 4.1; if it is notC=JIf the optimization should be ended, jump to step 7
And 5: calculating motion vector of measuring pointv 0I.e. linear equationsW v= bA solution ofv=v 0Wherein, in the step (A),v=[ v 1, v 2, v 3]T v 0=[ v 0,1, v 0,2, v 0,3,]T
step 6 is carried out after step 5 is finished;
step 6: solving a new key point by the tracing problem, and updating the state of the measuring point of the reference segment;
firstly, calculating the linear velocity of the measuring pointv j =W j v 0
Then, calculating the dynamic tracking time of each measuring point of the reference segmentτ j τ j,down =(t j -t min )÷(b j - v j ) | j 1,2, … M=τ j,up =(t max-t j )÷(b j - v j ) | j 1,2, … M=(ii) a Whereinτ j ={τ j,down ,τ j,up }。τ j,down Is the dynamic tracking time from the reference measuring point to the lower boundary,τ j,up is the dynamic time of tracing from the reference point to the upper boundary.
Then, the decision key point: dynamic tracking time of each measuring point of reference segmentτ j Minimum value ofτ j,min Corresponding measurementThe point is the next key pointr 3 And adding the corresponding serial number to the key serial number setrIn (c) }.
Then, updating a reference segment measuring point state set, wherein the measuring point state element set below the XOY plane is updated as follows:
t j =τ j,min v j -t min (ii) a The measuring point state element set of the measuring point positioned above the XOY surface of the coordinate is updated intot j =t max -τ j,min v i
Finally, according tou j = u j +τ min v j Updating reference segment measuring point setu j }; and simultaneously, updating the coordinates of the measurement point sets of all the measured segments in the same way.
Finishing one-time optimization after the step 6 is finished, and performing the step 2;
and 7: determine whether the reference section meets the flatness error requirement, i.e.t max -t min TAnd T is the tolerance value of the reference plane. If yes, jumping to step 8; if not, the evaluation is finished, and the part is not qualified.
And 8: obtainingu i ={x i , y i , z i I.e. measured plane fitting pointiAnd fitting the measured surface to the XOZ plane of the rectangular spatial coordinate system, i.e. the Z-axis of the coordinate system is perpendicular to the reference surface, in order to make the center plane of the measured surface close to the XOZ plane of the coordinate system. At the same time, the XOY plane of the coordinate system is close to the ideal reference plane. First, for measurement convenience, all the measurement points of the measured surface are projected to the XOY plane. Secondly, two contraction planes of the virtual gauge are established by points of the measuring points which are farthest away from the Y axis and the negative Y axis of the central plane of the virtual gauge respectively. Then, obtaining the initial key of the tested sectionSerial numberl 1 Andl 2 will bel 1 Andl 2 adding key serial number setlIn (1) };
l 1 is defined asu i Iny i Point of maximum valueu l, 1The serial number of (a) is included,l 2 is defined asu i Iny i Point of minimum valueu l, 2The serial number of (2);
then, setting boundary positions H and G;
H=y l,1y l, 1is composed ofu l, 1Is/are as followsyCoordinate values; h is the left boundary position, the left boundary is the crossing point [0, H, 0 ]]And parallel to the XOZ plane;
G=y l,2y l, 2is composed ofu l, 2Is/are as followsyCoordinate values; g is the right boundary position, the right boundary is the crossing point [0, G, 0 ]]And parallel to the XOZ plane;
obtaining the middle position midle of the left and right boundaries according to the above; midle = (H + G)/2;
finally, according to the measured point setu i Establishing state element set of measuring pointsc i Great Chinese character and feature line vector setW i }; wherein:
c i= y i status elements of all stationsc i Is a set of state elementsc i };
When y is i >At the time of the middle of the message,W i =([1, x i ,-1, -x i ])/y i ;y i <at the time of the middle of the message,W i =([-1, -x i ,1, x i ,])/y i W i is a special oneSyndrome row vector, characteristic row vector of all measurement pointsW i Is a set of characteristic line vectorsW i };
Step 9 is performed after step 8 is completed;
and step 9: according tou i Establishing a state element set (R) of a measured plane measuring point i And { E } and i };
R i = H - y i ;R i the distance from a measuring point of the measured plane to the left boundary;
E i =y i - G;E i the distance from a measuring point of the measured plane to the right boundary is obtained;
after step 9, performing step 10;
step 10: adding the measuring point serial number of a key point into a key serial number setlIn (1) }; if step 14 has not been performed, then a set of real-time state elements is takenc i Maximum value ofc max Corresponding measuring pointu l, 3Is a key point and the serial number of the measuring point isl 3Last page added to key serial number setlIn (1) }; at this time D/2 =c max D is the distance size from the boundary to the center; if step 14 generates a keypointu l,4Then, the key pointu l,4Will replace the measuring pointu l, 3Number of its measuring pointl 4Last page added to key serial number setlIn (1) };
step 11: according to the above-mentioned key point setlEstablishment of an analysis matrixWAnd analyzing the column vectorsbWherein:bis a number greater than zero, and:
W=[…, W u T, …, W v T, …]Tis a matrix with L rows and 4 columns, and L is a key point setlThe number of the elements in the (C),u, vis a set of key pointslThe elements in (1);
b=[…, b u , …, b v , …]Tis a column vector of L rows;
step 12 is performed after step 11 is completed;
step 12: for the above-mentioned established analysis matrixWAnd an augmented analysis matrixW, b]Analyzing the rank;
first computing analysis matrixWRank ofr W =rank(W) Extended analysis matrix [ alpha ], [ beta ] anW, b]Rank ofr Wb =rank([W, b]) And further compare the rankr W Sum rankr Wb By verification, there are only two cases:
the first situation is as follows: if it is notr W = r Wb If yes, the optimization should be continued, and the step 13 is skipped;
case two: if it is notr W < r Wb Then go to step 12.1;
step 12.1: let the key point collectlGet rid of the firstkAn elementl k Then obtainl q According tol q Establishment of a reduction matrixW q And reducing the column vectorb q Wherein:kis 1;
W q =[…, W u T, …, W v T, …]Tis aQA matrix of rows and 4 columns,Qis a setl q The number of the elements in the (C),uv isl q The elements in (1);
b q =[…, b u , …, b v , …]Tis aQA column vector of rows;
step 12.2: solving linear equationsW q f q = b q Solution of (2)f q =f o Then calculateb k =W q f o (ii) a If it is notb k > b q Then the matrix will be reducedW q And reducing the column vectorb q Respectively asWMatrix and analysis column vectorbWill elementl k Moving out key point setlAnd jumping to the step 13; if it is notb k <=b q Then step 12.3 is performed;
step 12.3: judgment ofkWhether or not equal toQIf, ifkIs not equal toQThen jump to step 11.1; if it is notkIs equal toQIf the optimization should be finished, jump to step 15; wherein the content of the first and second substances,f q =[f q,1 , f q,2 , f q,3 , f q,4 ]Tf o =[f q0,1, f q0,2, f q0,3, f q0,4]T
step 13: motion vector of measuring point and virtual gauge is solvedf o I.e. linear equationsWf = bA solution off=f o Wherein:
f=[f 1, f 2, f 3, f 4]Tf o =[f o,1, f o,2, f o,3, f o,4]T
step 14 is performed after step 13 is completed;
step 14: finding new key point by tracing problemu i Updating is carried out;
first, calculating the relative speed between each measuring point and the virtual gaugef i }:f i =W i f o
Then, the time of pursuit of each measurement point is calculatedt i Andt i ,,
t i =R i / f i | i=1, 2, ...Mt i trace the left boundary with the measuring pointu i The time of (d);
t i ,,=E i / f i | i=1, 2, ...Mt i ,,trace the right border with the measuring pointu i The time of (d);
t i gett i Andt i ,,the smaller of these;
then, the decision key point: to catch up timet i Minimum value in the part of greater than zerot minThe corresponding measuring points are key pointsu l ,3And corresponding the serial number to the pointl 3Adding key serial number setl};
Then, according tou i =u i + t min f i Updating measuring point setu i };
Finally, according toH= H - t min bUpdating the position of the left boundaryH(ii) a According toG= G + t min bUpdating right boundary positionG
Finishing one-time optimization after the step 14 is finished, and performing the step 9;
step 15: the actual size of the measured plane isH - GIf (a), (b) isH - G)< dIf the measured plane meets the verticality requirement, the measured plane is qualified; otherwise, the product is not qualified. Wherein the content of the first and second substances,dis the tolerance value of verticality.
To facilitate numerical calculations, the following optimizations may be performed:
order tobTaking a value greater than 0, may be, but is not limited to, 1.
The measured elements are side surface measuring point sets of the measured workpiece.
The invention has the beneficial effects that:
1. the geometrical characteristics of the L-shaped measurement workpiece verticality based on the single reference plane are fully considered, the computer is used for replacing the considered manual evaluation, and the evaluation mode is simple and convenient, so that the method has higher efficiency and is easier to popularize compared with the first type of evaluation method. 2. The geometrical characteristics of the perpendicularity of the L-shaped measuring workpiece based on the single datum plane are fully considered, a better value is obtained through mature linear operation in each iteration, and the minimum limit equivalent size can be finally obtained, so that the algorithm is stable, and the problem that the evaluation indicator of the second method has errors is solved. 3. The fact that most of measuring points are invalid measuring points in the evaluation of the perpendicularity of the shaft hole type part based on the single datum plane is implied, and the invalid measuring points are not added with iteration, so that the iteration times are less, and the method is more suitable for evaluating large-sized workpieces compared with a third evaluation method. 4. Calculating the optimizing direction by considering only the keypoint setkAnd (4) corresponding measuring points, so that the operation amount of each iteration is small, the mathematical model is established simply, and compared with the fourth type of evaluation method, the method has the advantages of equivalent total operation speed, lower mathematical modeling complexity and worth popularizing.
The invention provides a method for quickly evaluating the perpendicularity of an L-shaped measuring workpiece based on a single reference plane, which is stable, quick and simple in form, can be used for evaluating the error of the L-shaped measuring workpiece based on the perpendicularity of the single reference plane, and provides guidance for the improvement of the machining process, so that the method has industrial possibility.
Drawings
FIG. 1 is a baseline assessment summary flow chart A of the present invention.
FIG. 2 is a flowchart B of the measured element evaluation method of the present invention.
FIG. 3 is a tolerance layout for parts in an exemplary embodiment.
Detailed Description
The following are specific examples of the present invention, and the aspects of the present invention will be further described with reference to the drawings, but the present invention is not limited to these examples.
And evaluating the verticality error of the L-shaped measuring workpiece based on the single datum plane, wherein the tolerance design specification is shown in figure 3.
Step 1: measuring point set for obtaining reference segmentu j Great moment and measuring point set of the segment to be measuredu i The method comprises the following steps:
Figure DEST_PATH_IMAGE001
the initial key serial numbers at this time are 6 and 8;
at this time, the upper boundary position of the reference plane is z'6=0.0514 and lower boundary is z'8=0.0014;
Obtaining a middle position of the boundary as midle = 0.0264;
establishing a characteristic line vector set of reference plane measuring pointsW j The method comprises the following steps:
Figure 534985DEST_PATH_IMAGE002
after step 1, step 2 is performed.
Step 2: according tou j Establishing a reference state element sett j Ant' j }:
Figure DEST_PATH_IMAGE003
After step 2, step 3 is performed.
And step 3: establishing an analysis matrix according to the key sequence numberWAnd analyzing the column vectorsb
W= [ W 6 T , W 8 T ]T =
Figure 452125DEST_PATH_IMAGE004
b=[1, 1]T
After step 3, step 4 is performed.
And 4, step 4:R W =rank(W) =2,R W b =rank([W, b])=2,R W =R W b and jumping to the step 5.
After step 4, step 5 is performed.
And 5: linear equation of equationsW v = bIs thatv = v 0 =[ 0, 0.1000, 0.2000]T
Step 5 ends and proceeds to step 6.
Step 6: calculating the normal linear velocity of the measuring pointv j v j =W j v 0The results are as follows:
Figure DEST_PATH_IMAGE005
then, the minimum tracking time of each measuring point is calculatedτ j And deciding a key point: to catch up timeτ j Of greater than zero
Minimum value ofτ j,min The measured point corresponding to the value of =0.00033 is a key pointu 13
Finally, updating the reference measuring point setu j Greatu j The method comprises the following steps:
Figure 44912DEST_PATH_IMAGE006
updating the upper and lower boundary positions of the reference, wherein the upper boundary is z'6=0.0511, lower boundary z'8=0.0018。
And 6, finishing the first optimization of the benchmark evaluation after the step 6 is finished, and returning to the step 2.
After the benchmark evaluation optimization is completed according to the invention, all the key sequence number sets are obtained as a great facer}={6,8,13,15};
Step 6 obtains a new key point ofu 15
Step 6, after finishing step 2;
step 2: according tou j Establishing a reference state element sett j Ant' j };
And step 2 is finished and step 3 is carried out.
And step 3:W = [ W 6 T , W 8 T , W 13 T ,W 15 T ]T=
Figure DEST_PATH_IMAGE007
b= [1,1,1,1]T
after step 3, step 4 is performed.
And 4, step 4:R W = rank(W) = 3,R Wb =rank([W, b])= 4,R W R Wb step 4.1 is performed.
Step 4.1: making key serial number collect lastrRemoving the 1 st element 6 to obtain {8, 13, 15} = {6, 8, 13, 15}, and establishing a reduction matrix according to {8, 13, 15}W r- And reducing the column vectorb r- Respectively is as follows:
W r- =[ W 8 T , W 13 T ,W 15 T ]T=
Figure 502438DEST_PATH_IMAGE008
b= [1,1,1]T
step 4.2: linear equation of equationsW r- v r- = b r- Solution of (2)v r- =[5.5000,-0.1000,0.1500]Calculatingbr=W r- v r- = -3<1。
Step 4.3 is performed.
Step 4.3: counter C =1<4, jump to step 4.1.
Step 4.1: making key serial number collect lastrRemoval of the 2 nd element 8 yields a retaining openingr} = {6, 13, 15}, and a reduction matrix is established according to {6, 13, 15}W r- And reducing the column vectorb r- Respectively is as follows:
W r- =[ W 6 T , W 13 T ,W 15 T ]T=
Figure DEST_PATH_IMAGE009
b= [1,1,1]T
step 4.2: linear equation of equationsW r- v r- = b r- Solution of (2)v r- =[-0.5000,-0.1000,-0.0500]Calculatingbr=W r- v r- = -3.1862<1。
Step 4.3 is performed.
Step 4.3: counter C =2<4, jump to step 4.1.
Step 4.1: making key serial number collect lastrRemoving the 3 rd element 13 to obtain a retaining openingr} = {6, 8, 15}, and a reduction matrix is established according to {6, 8, 15}W r- And reducing the column vectorb r- Respectively is as follows:
W r- =[ W 6 T , W 8 T ,W 15 T ]T=
Figure 316810DEST_PATH_IMAGE010
b= [1,1,1]T
step 4.2: linear equation of equationsW r- v r- = b r- Solution of (2)v r- =[-1.5000,0.1000,0.0500]Calculatingbr=W r- v r- = -3.2954<1。
Step 4.3 is performed.
Step 4.3: counter C =3<4, jump to step 4.1.
Step 4.1: making key serial number collect lastrRemoval of the 4 th element 15 yields a retaining openingr} = {6, 8, 13}, and a reduction matrix is established according to {6, 8, 13}W r- And reducing the column vectorb r- Respectively is as follows:
W r- =[ W 6 T , W 8 T ,W 13 T ]T=
Figure DEST_PATH_IMAGE011
b= [1,1,1]T
step 4.2: linear equation of equationsW r- v r- = b r- Solution of (2)v r- =[-1.5000,0.1000,0.0500]Calculatingbr=W r- v r- = -1.3495<1. Step 4.3 is performed.
Step 4.3: counter C =4, jump to step 7.
And 7: judging whether the reference section meets the requirement of flatness error:t max -t min =0.0501- 0.0028=0.0473≦0.05=Tthe reference plane meets the flatness requirement, wherein T is the tolerance value of the reference plane,t max the distance from the measuring point to the upper boundary z,t min is the distance from the measured point to the lower boundary z. Step 8 is performed.
And 8: measured plane measuring point set after benchmark evaluation is completed{u i The method comprises the following steps:
Figure 946463DEST_PATH_IMAGE012
the initial key serial numbers at this time are 7 and 15;
the left boundary of the reference plane is H =0.0449 and the right boundary is G = -0.0036 at this time;
obtaining a boundary middle position as midle = 0.0207;
establishing a set of state elementsc i Anc' i The method comprises the following steps:
Figure DEST_PATH_IMAGE013
establishing a feature line vector setWiThe method comprises the following steps:
Figure 275813DEST_PATH_IMAGE014
step 8 ends and step 9 is performed.
And step 9: updating a real-time state element set { R } i And { E } and i therein ofy i<At the time of the middle of the message,c i = R i y i>at the time of the middle of the message,c i = E i
step 9 ends and proceeds to step 10.
Step 10: if step 14 has not been performed, a set of real-time state elements is takenc i Maximum value y of }max =0.0485 corresponding measuring pointu 15The measuring point serial number 15 of the key point is listed as a set of key points during the initialization processlIn (1) }; at this time D =0.0485, D being the virtual gauge size.
Step 10 ends and proceeds to step 11.
Step 11:W q = [ W 7 T , W 15 T ]T =
Figure DEST_PATH_IMAGE015
b q =[1, 1]T
step 11 ends and proceeds to step 12.
Step 12:r W =rank(W) =2,r Wb =rank([W, b])=2,r W = r Wb and a jump is made to step 13.
Step 13: linear equation of equationsWq f q = b q Is thatf q = f o =[0.0656 , -0.0021 , 0 , 0]T
Step 13 ends and step 14.
Step 14: first, according tof i =W i f o Calculating the relative speed between each measuring point and the virtual gaugef i The method comprises the following steps:
Figure 287762DEST_PATH_IMAGE016
secondly, a new key point is found by tracing the problem,t i =R i / f i | i=1, 2, ...M
then, the decision key point: to catch up timet i Minimum value in the part of greater than zerot i minThe measured point corresponding to the value of =0.00083 is a key pointu 5And adding the serial number of the key measuring point into the key point setlIn (1) };
then, according tou i =u i + t min f i Updating measuring point setu i As follows:
Figure DEST_PATH_IMAGE017
and finishing one-time optimization after the step is finished, and performing the step 9.
And step 9: updating a set of real-time state elementsc i }
Figure 956641DEST_PATH_IMAGE018
After step 9, step 10 is performed.
Step 10: according to the invention, new key points are searched in the tracking time, and after all the optimizations are completed, a key sequence number set is obtainedl} = {7, 15, 5}, when D =c max =0.0485 distance size of left and right boundaries of gauge;
after step 10, step 11 is performed.
Step 11:W=[W 7 TW 15 TW 5 T]T=
Figure DEST_PATH_IMAGE019
b=[1,1,1]T
step 11 ends and proceeds to step 12.
Step 12:r W = rank(W) =2,r Wb = rank([W , b])=3,r W r Wb step 12.1 is performed.
Step 12.1: the key sequence number set {7, 15, 5} gets {15, 5} after the first element 7 is removed, and a reduced matrix is built according to {15, 5}W q And reducing the column vectorb q Respectively is as follows:
W q =[W 15 TW 5 T]T=
Figure 12322DEST_PATH_IMAGE020
b q =[1,1]T
after step 12.1 is finished, step 12.2 is performed.
Step 12.2: solving linear equationsW q f q = b q Is solved byf q = f o =[-0.0046,0.00027,0,0]Then calculateb k =W q f o =-2.5<1; namely, it isb k <=b q Then step 12.3 is performed.
Step 12.3: counter with a memoryk =1<3, jumping to step 12.1.
Step 12.1: the second element 15 of the key sequence number set {7, 15, 5} is removed to obtain {7, 5}, and a reduced matrix is established according to {7, 5}W q And reducing the column vectorb q Respectively is as follows:
W q =[W 7 TW 5 T]T=
Figure DEST_PATH_IMAGE021
b q =[1,1]T
after step 12.1 is finished, step 12.2 is performed.
Step 12.2: solving linear equationsW q f q = b q Is solved byf q = f o =[114.9624,-11.4964,0,0]Then calculateb k =W q f o =-6.8074<1; namely, it isb k <=b q Then step 12.3 is performed.
Step 12.3: counter with a memoryk =2<3, jumping to step 12.1.
Step 12.1: order keyThe sequence number set {7, 15, 5} gets {7, 15} after the third element 5 is removed, and a reduced matrix is built according to {7, 15}W q And reducing the column vectorb q Respectively is as follows:
W q =[W 7 TW 15 T]T=
Figure 579700DEST_PATH_IMAGE022
b q =[1,1]T
after step 12.1 is finished, step 12.2 is performed.
Step 12.2: solving linear equationsW q f q = b q Is solved byf q = f o =[0.0656,-0.0021,0,0]Then calculateb k =W q f o =-24.1818<1; namely, it isb k <=b q Then step 12.3 is performed.
Step 12.3: counter with a memoryk =3, jump to step 15
Step 15: at the moment, the effective size D =0.0485 between the left plane and the right plane of the measured plane is less than 0.05, the measured plane meets the verticality requirement, and the measured plane is qualified.
To conclude: the tested L-shaped error measurement workpiece based on single reference meets the verticality requirement given by the figure 3
In the above description, the present invention has been described by way of specific embodiments, but those skilled in the art will appreciate that various modifications and variations can be made within the spirit and scope of the invention as hereinafter claimed.
The verticality measurement result of the L-shaped error measurement workpiece meets the definition and requirement of the verticality in GB/T1182-2018.
The measurement result of the invention shows that the L-shaped error measurement workpiece can be used for measuring the verticality error of the part with the precision requirement.

Claims (6)

1. The invention takes an L-shaped error measurement workpiece as an example, and verifies a single-datum-plane-based quick evaluation method for perpendicularity of the L-shaped error measurement workpiece; the specific placing position of the workpiece is that the bottom surface of the L-shaped workpiece is attached to a horizontal plane, and the side surface of the L-shaped workpiece is vertical to the bottom surface;
the bottom surface of the L-shaped workpiece is used as a reference surface, and the vertical surface is a measured plane.
2. A single datum plane based method for quickly evaluating the verticality of an error measurement workpiece is characterized by comprising the following steps of:
step 1: firstly, measuring points of a reference surface A are obtained and are used to form a reference measuring point setu j And according to au j Establishing a characteristic line vector setW j Greatb j }; obtaining initial measuring point of measured element and using it to form measuring point setu i }; wherein:
j=1, 2, 3, …, Mjthe serial number of the reference measuring point is,Mthe total number of the measuring points of the reference plane is;
i=1, 2, 3, …, Nithe serial number of the measuring point of the measured segment,Nthe total number of the measuring points of the measured plane is;
u j ={x j , y j , z j is the measurement pointjAnd the Z-axis of the coordinate system is perpendicular to the reference plane, the XOY plane of the coordinate system is close to the ideal reference plane;
t j =z j all state elementst j Is a set of state elementst j };
When z is j >At the time of 0, the number of the first,W j =([-1, -y j , x j ]) (ii) a When z is j <At the time of 0, the number of the first,W j =([1, y j ,- x j ]);W j is a feature row vector, all feature row vectorsW j Is a set of characteristic line vectorsW j };
b j =bIs a real number greater than 0, all boundary elementsb j Is a set of boundary elementsb j };
After the step 1 is finished, performing a step 2;
step 2: gett j Maximum valuet max Corresponding serial numberr 2 Is numbered for the key point and willr 2 Last page added to key serial number setrIn (c) } the reaction solution is,t max the distance from the reference measuring point to the upper boundary; gett j Minimum valuet min Corresponding serial numberr 1 Is numbered for the key point and willr 1 Last page added to key serial number setrIn (c) } the reaction solution is,t min the distance from the reference measuring point to the lower boundary;
step 3 is carried out after step 2 is finished;
and step 3: according to the key sequence numberrEstablishment of an analysis matrixWAnd analyzing the column vectorsbWherein:
W=[…, A p T, …, A q T, …]Tis aFA matrix of rows and 4 columns,Fis a set of key pointsrThe number of the elements in the (C),p, qis a set of key pointsrThe elements in (1);
b=[…, b p , …, b q , …]Tis aFA column vector of rows;
step 4 is carried out after step 3 is finished;
and 4, step 4: for analysis matrixWAnd an augmented analysis matrixW, b]Performing rank analysis;
computing an analysis matrixWRank ofR W =rank(W) Extended analysis matrix [ alpha ], [ beta ] anW, b]Rank ofR Wb =rank([W, b]) And compareR W AndR Wb there are only two cases:
the first condition is as follows: if it is notR W =R Wb Then, the optimization should be continued, jumping to step 5;
case two: if it is notR W <R Wb Then, the calculator C =1 is initialized and step 4.1 is performed;
step 4.1: attempting to slave an analysis matrixWAnd analyzing the column vectorsbMiddle removed key point setrGet the line corresponding to one element in the matrix to get the reduced matrixW r- And reducing the column vectorb r- (ii) a WhereinW r- Is oneJA matrix of rows and 4 columns,br-is oneJA column vector of rows;
step 4.2: solving linear equationsW r- v r- = b r- Solution of (2)v r- =v r-0 Then calculateb r =W r v r-0 (ii) a If it is notb r- > b r Then, the matrix will be reducedW r- And reducing the column vectorb r- Respectively asWMatrix and analysis column vectorbWill elementrMoving out key point setrAnd jumping to the step 5; if it is notb r- <= b r Then step 4.3 is performed;
step 4.3: judging counterCWhether or not equal toJIf, ifCJJumping to step 4.1; if it is notC=JShould end the optimization and jump to stepStep 7
And 5: calculating motion vector of measuring pointv 0I.e. linear equationsW v= bA solution ofv=v 0Wherein, in the step (A),v=[ v 1, v 2, v 3]T v 0=[ v 0,1, v 0,2, v 0,3,]T
step 6 is carried out after step 5 is finished;
step 6: solving a new key point by the tracing problem, and updating the state of the measuring point of the reference segment;
firstly, calculating the linear velocity of the measuring pointv j =W j v 0(ii) a Then, calculating the dynamic tracking time of each measuring point of the reference segmentτ j
τ j,down =(t j -t min )÷(b j - v j ) | j 1,2, … M=τ j,up =(t max-t j )÷(b j - v j ) | j 1,2,… M=(ii) a Whereinτ j ={τ j,down ,τ j,up };τ j,down Is the dynamic tracking time from the reference measuring point to the lower boundary,τ j,up is the dynamic time of tracing from the reference measuring point to the upper boundary;
secondly, the decision key points are: dynamic tracking time of each measuring point of reference segmentτ j Minimum value ofτ j,min The corresponding measuring point is the next key pointr 3 And adding the corresponding serial number to the key serial number setrIn (1) }; then, updating a reference segment measuring point state set, wherein the measuring point state element set below the XOY plane is updated as follows:t j =τ j,min v j -t min (ii) a The measuring point state element set of the measuring point positioned above the XOY surface of the coordinate is updated intot j =t max -τ j,min v i
Finally, according tou j = u j +τ minv j Updating reference segment measuring point setu j }; simultaneously, updating the coordinates of the measured point sets of all the measured sections in the same way;
finishing one-time optimization after the step 6 is finished, and performing the step 2;
and 7: determine whether the reference section meets the flatness error requirement, i.e.t max -t min TWherein T is a tolerance value of the reference plane; if yes, jumping to step 8; if not, finishing the evaluation, and enabling the part not to meet the requirements;
and 8: obtainingu i ={x i , y i , z i I.e. measured plane fitting pointiFitting the measured surface on an XOZ plane of a space rectangular coordinate system, namely, the Z axis of the coordinate system is vertical to a reference surface, so that the central plane of the measured surface is close to the XOZ plane of the coordinate system; meanwhile, the XOY plane of the coordinate system is close to an ideal reference plane; firstly, projecting all measuring points of a measured surface to an XOY plane for convenient measurement; secondly, establishing two contraction planes of the virtual gauge by using points at which the measuring points are farthest away from the Y axis and the negative Y axis of the central plane of the virtual gauge respectively; then, obtaining the initial key sequence number of the segment to be testedl 1 Andl 2 will bel 1 Andl 2 adding key serial number setlIn (1) };
l 1 is defined asu i Iny i Point of maximum valueu l, 1The serial number of (a) is included,l 2 is defined asu i Iny i Point of minimum valueu l,2The serial number of (2);
then, setting boundary positions H and G;
H=y l,1y l, 1is composed ofu l,1Is/are as followsyCoordinate values; h is the left boundary position, the left boundary is the crossing point [0, H, 0 ]]And parallel to the XOZ plane;
G=y l,2y l, 2is composed ofu l,2Is/are as followsyCoordinate values; g is the right boundary position, the right boundary is the crossing point [0, G, 0 ]]And parallel to the XOZ plane;
obtaining the middle position midle of the left and right boundaries according to the above; midle = (H + G)/2;
finally, according to the measured point setu i Establishing state element set of measuring pointsc i Great Chinese character and feature line vector setW i }; wherein:
c i= y i status elements of all stationsc i Is a set of state elementsc i };
When y is i >At the time of the middle of the message,W i =([1, x i ,-1, -x i ])/y i ;y i <at the time of the middle of the message,W i =([-1, -x i ,1, x i ,])/y i W i is a characteristic line vector, the characteristic line vectors of all measuring pointsW i Is a set of characteristic line vectorsW i };
Step 9 is performed after step 8 is completed;
and step 9: according tou i Establishing a state element set (R) of a measured plane measuring point i And { E } and i };
R i = H - y i ;R i the distance from a measuring point of the measured plane to the left boundary;
E i =y i - G;E i the distance from a measuring point of the measured plane to the right boundary is obtained;
after step 9, performing step 10;
step 10: adding the measuring point serial number of a key point into a key serial number setlIn (1) }; if step 14 has not been performed, then a set of real-time state elements is takenc i Maximum value ofc max Corresponding measuring pointu l,3Is a key point and the serial number of the measuring point isl 3Last page added to key serial number setlIn (1) }; at this time D/2 =c max D is the distance size from the boundary to the center; if step 14 generates a keypointu l,4Then, the key pointu l,4Will replace the measuring pointu l,3Number of its measuring pointl 4Last page added to key serial number setlIn (1) };
step 11: according to the above-mentioned key point setlEstablishment of an analysis matrixWAnd analyzing the column vectorsbWherein:bis a number greater than zero, and:
W=[…, W u T, …, W v T, …]Tis a matrix with L rows and 4 columns, and L is a key point setlThe number of the elements in the (C),u, vis a set of key pointslThe elements in (1);
b=[…, b u , …, b v , …]Tis a column vector of L rows;
step 12 is performed after step 11 is completed;
step 12: for the above-mentioned established analysis matrixWAnd an augmented analysis matrixW, b]Analyzing the rank;
first computing analysis matrixWRank ofr W =rank(W) Enhanced analysis matrix[W, b]Rank ofr Wb =rank([W, b]) And further compare the rankr W Sum rankr Wb By verification, there are only two cases:
the first situation is as follows: if it is notr W = r Wb If yes, the optimization should be continued, and the step 13 is skipped;
case two: if it is notr W < r Wb Then go to step 12.1;
step 12.1: let the key point collectlGet rid of the firstkAn elementl k Then obtainl q According tol q Establishment of a reduction matrixW q And reducing the column vectorb q Wherein:kis 1;
W q =[…, W u T, …, W v T, …]Tis aQA matrix of rows and 4 columns,Qis a setl q The number of the elements in the (C),uv isl q The elements in (1);
b q =[…, b u , …, b v , …]Tis aQA column vector of rows;
step 12.2: solving linear equationsW q f q = b q Solution of (2)f q =f o Then calculateb k =W q f o (ii) a If it is notb k > b q Then the matrix will be reducedW q And reducing the column vectorb q Respectively asWMatrix and analysis column vectorbWill elementl k Moving out key point setlAnd jumping to the step 13; if it is notb k <=b q Then step 12.3 is performed;
step 12.3: judgment ofkWhether or not equal toQIf, ifkIs not equal toQThen jump to step 11.1; if it is notkIs equal toQIf the optimization should be finished, jump to step 15; wherein the content of the first and second substances,f q =[f q,1 , f q,2 , f q,3 , f q,4 ]Tf o =[f q0,1, f q0,2, f q0,3, f q0,4]T
step 13: motion vector of measuring point and virtual gauge is solvedf o I.e. linear equationsWf= bA solution off=f o Wherein:
f=[f 1, f 2, f 3, f 4]Tf o =[f o,1, f o,2, f o,3, f o,4]T
step 14 is performed after step 13 is completed;
step 14: finding new key point by tracing problemu i Updating is carried out;
first, calculating the relative speed between each measuring point and the virtual gaugef i }:f i =W i f o
Then, the time of pursuit of each measurement point is calculatedt i Andt i ,,
t i =R i / f i | i=1, 2, ...Mt i trace the left boundary with the measuring pointu i The time of (d);
t i ,,=E i / f i | i=1, 2, ...Mt i ,,trace the right border with the measuring pointu i The time of (d);
t i gett i Andt i ,,the smaller of these;
then, the decision key point: to catch up timet i Minimum value in the part of greater than zerot minThe corresponding measuring points are key pointsu l, 3And corresponding the serial number to the pointl 3Adding key serial number setl};
Then, according tou i =u i + t minf i Updating measuring point setu i };
Finally, according toH= H - t min bUpdating the position of the left boundaryH(ii) a According toG= G + t min bUpdating right boundary positionG
Finishing one-time optimization after the step 14 is finished, and performing the step 9;
step 15: the actual size of the measured plane isH - GIf (a), (b) isH - G)< dIf the measured plane meets the verticality requirement, the measured plane is qualified; otherwise, the product is unqualified; wherein the content of the first and second substances,dis the tolerance value of verticality.
3. The single reference plane-based error measurement workpiece perpendicularity rapid assessment method according to claim 1, wherein the measurement data of the general reference plane and the measured plane fitting data are subjected to coordinate transformation to obtain the XOZ plane of the measured plane three-dimensional center plane close to the coordinate system, the reference plane close to the XOY plane of the coordinate systemMeasuring point set of faceu i }。
4. The method for rapidly evaluating the verticality of a single-datum-plane-based error measurement workpiece according to claim 2, wherein the coordinate transformation is as follows: moving according to the average value of the coordinates; or two: moving according to the extreme value of the coordinate; or three: and moving according to the root-mean-square minimum principle of coordinates.
5. The method for rapidly evaluating the verticality of a workpiece based on single reference plane error measurement according to claim 1, wherein the method comprises the steps ofb =1。
6. The single-datum-plane-based error measurement workpiece perpendicularity rapid assessment method according to claim 1, wherein the measured plane is a vertical plane perpendicular to the measured plane, and the datum plane is a vertical plane;
the verticality measurement result of the L-shaped error measurement workpiece meets the definition and requirement of verticality in GB/T1182-2018; the measurement result of the invention shows that the L-shaped error measurement workpiece can be used for measuring the verticality error of the part with the precision requirement.
CN202010602658.5A 2020-06-29 2020-06-29 Method for quickly evaluating verticality of error measurement workpiece based on single reference plane Pending CN113932754A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010602658.5A CN113932754A (en) 2020-06-29 2020-06-29 Method for quickly evaluating verticality of error measurement workpiece based on single reference plane

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010602658.5A CN113932754A (en) 2020-06-29 2020-06-29 Method for quickly evaluating verticality of error measurement workpiece based on single reference plane

Publications (1)

Publication Number Publication Date
CN113932754A true CN113932754A (en) 2022-01-14

Family

ID=79272650

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010602658.5A Pending CN113932754A (en) 2020-06-29 2020-06-29 Method for quickly evaluating verticality of error measurement workpiece based on single reference plane

Country Status (1)

Country Link
CN (1) CN113932754A (en)

Similar Documents

Publication Publication Date Title
CN110285781B (en) Rapid assessment method for plane parallelism relative to reference plane
CN105783786B (en) Part chamfering measuring method and device based on structured light vision
CN106354942B (en) Assembly contact stress distribution evaluation method based on entropy theory
CN107492120B (en) Point cloud registration method
CN104613897B (en) The adaptively sampled device of vehicle exterior cladding element mold freedom curved surface and measuring method
CN107063060A (en) A kind of method and device for determining surface planarity
CN106406237B (en) A kind of processing method with free form surface metal parts
CN103256916A (en) Evaluation method of part flatness error based on minimum area
CN113516695B (en) Point cloud registration strategy in laser profiler flatness measurement
CN103791878A (en) Numerically-controlled machine tool geometric accuracy identification method
CN102980532B (en) Method for measuring large-diameter aspheric surface shapes in splicing manner by adopting three-coordinate measuring machine
CN103868455A (en) Method for achieving visual reconstruction of space coordinates of target point in water tank
CN113932754A (en) Method for quickly evaluating verticality of error measurement workpiece based on single reference plane
CN105389476A (en) Interpolation algorithm for intensity-modulated radiation therapy plan dose data based on gradient features
CN103411562A (en) Structured light laser strip center extraction method based on dynamic programming and mean-shift
CN108267106B (en) Quick, stable and simple cylindricity error evaluation method
CN113267122B (en) Industrial part size measurement method based on 3D vision sensor
CN109508482A (en) A kind of calculation method for complex-curved surface profile degree error uncertainty
TWI444586B (en) System and method for detecting form-position tolerances of an object
CN109446539A (en) In conjunction with surface error body modeling method and device
CN112419394A (en) Turnout contour loss detection method and device
CN113483721A (en) Method for rapidly evaluating straightness in any direction by applying maximum entity requirement
CN110619139A (en) Rapid position degree evaluation method based on single-datum hole shaft part
CN113483718A (en) Method for rapidly evaluating straightness in any direction based on minimum entity requirement
CN113884059A (en) Single-datum plane shaft hole part gradient evaluation method based on rank analysis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20220114

WD01 Invention patent application deemed withdrawn after publication