CN102809364B - Method for determining complex curved surface profile error - Google Patents

Method for determining complex curved surface profile error Download PDF

Info

Publication number
CN102809364B
CN102809364B CN201210236183.8A CN201210236183A CN102809364B CN 102809364 B CN102809364 B CN 102809364B CN 201210236183 A CN201210236183 A CN 201210236183A CN 102809364 B CN102809364 B CN 102809364B
Authority
CN
China
Prior art keywords
error
complex
described
curved profile
amp
Prior art date
Application number
CN201210236183.8A
Other languages
Chinese (zh)
Other versions
CN102809364A (en
Inventor
何改云
张美�
刘欣
贾红洋
Original Assignee
天津大学
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 天津大学 filed Critical 天津大学
Priority to CN201210236183.8A priority Critical patent/CN102809364B/en
Publication of CN102809364A publication Critical patent/CN102809364A/en
Application granted granted Critical
Publication of CN102809364B publication Critical patent/CN102809364B/en

Links

Abstract

The invention discloses a method for determining a complex curved surface profile error, which relates to the field of complex curved surface evaluation. The method comprises the following steps of transmitting a current coordinate position to a data processing terminal by a numerical control system; building a complex curved surface profile error model by the data processing terminal through the current coordinate position; obtaining a position pose initial parameter of an ideal curved surface through the rough matching of a measuring coordinate system and a design coordinate system; obtaining a complex curved surface profile initial error through the position pose initial parameter and the complex curved surface profile error model by the data processing terminal; and optimizing the position pose initial parameter by the data processing terminal, obtaining a complex curved surface profile current error again, and using the complex curved surface profile current error as a complex curved surface profile final error and outputting the complex curved surface profile final error by the data processing terminal until a difference value between the complex curved surface profile current error and the complex curved surface profile initial error is less than a threshold value. The complex curved surface profile final error obtained through the method is close to a theoretical value which is accordant with a minimum condition, the workpiece delay waste rate and the production cost are reduced, and the detection accuracy is increased.

Description

A kind of method of determining complex-curved profile error

Technical field

The present invention relates to complex-curved error evaluation field, relate in particular to a kind of method of determining complex-curved profile error.

Background technology

In fields such as Aero-Space, shipbuilding and auto industrys, many important parts are often made up of the free profile of high precision, life-span and the security of these parts to whole equipment plays an important role, so free form surface is being played the part of more and more important role in engineering reality, due to the process complexity of free form surface, processing cost is too high, so its error evaluation process is had to higher requirement, require to evaluate exactly wrap-around error as far as possible, make it approach true error, give up to reduce production costs by avoiding by mistake.

In NC Machining Process, the classic method that workpiece is measured will be used various conventional measurers (clock gauge and slide calliper rule etc.) or equipment (three coordinate measuring machine etc.).For simple geometric measurement, operator can be directly with gage measuring on numerically-controlled machine; For the measurement problem of the more complicated such as geometric properties or locus of workpiece, because conventional measurer cannot solve, operator must carry out workpiece transfer to three coordinate measuring machine.But while using this measuring method, will carry workpiece, the sequence of operations such as clamping and location again, tend in the meantime produce error, this can exert an influence to real error, reduces the efficiency of evaluation.

On the other hand along with CAD/CAM/CAI(computer aided design/computer aided machine/computer aided detection) development and machining precision from micron to nano level raising, GPS(product geometric sense technical manual) also there is huge variation in system.New generation of GPS is by the comprehensive call that provides emphatically to be suitable for every technology such as CAX(CAD/CAM/CAI) the integration environment, the geometric tolerances definition of clear and definite, system specifications more.Its most outstanding feature be exactly emphasize with mathematical method to the geometric error of product be described, definition, modeling and information transmission.The research of the mathematical definition theory to tolerance is carried out early in the U.S., and has promulgated first tolerance mathematical definition standard in the world.

How accurately to obtain profile error value and depend on concisely reasonably error model and fast and effeciently optimized algorithm, in recent years, scholar both domestic and external for provide various geometric configuratioies tolerance mathematical definition, set up profile error mathematical model, find efficiently fast optimized algorithm all made further investigation, for example:

(1) U.Roy and B.Li studied form and position tolerance mathematical modeling (referring to Roy U, Li B, for annotation and the expression of polytopic plants geometric tolerances (1. form tolerance), Computer Aided Design, 1998,30 (2): 151~161.)

(2) University Of Tianjin studied the mathematical definition of sphericity, complex position degree tolerance and error model (referring to He Gaiyun. the approximation theory of morpheme error and algorithm research [D]. Tianjin: University Of Tianjin, 2006.)

(3) Cai Min, poplar by the mathematical model of having studied CYLINDRICAL FEATURE such as new (referring to Cai Min, poplar by new etc., the research [J] of the CYLINDRICAL FEATURE mathematical model based on mathematical definition. mechanical engineering journal, 2003,39 (12): 86~90)

(4) Jiang Guohua etc. has studied cylindricity and has wanted geometric error evaluation algorithm, by conic projection on plane, so the evaluation of deviation from cylindrical form is converted into the evaluation of deviation from circular from, rotational plane again, obtain minimum error amount (referring to S.Hossein, Guohua Jiang, etal, evaluation [J] the .Precision Engineering of cylinder key element geometric properties, 2003,27 (2): 195~204)

(5) Zhu Limin etc. has been set up the mathematical model of point to curved surface directed distance, and has obtained function the second Taylor series formula by the method for difference, and is optimized taking expansion as objective function, thereby obtains error amount to be asked.(referring to Limin Zhu, Zhenhua Xiong, etal, complex-curved location based on distance function and profile error assessment method [J] .Journal of Manufacturing Science and Engineering, 2004,126 (3): 542~554), (Li Min Zhu, Xiao Ming Zhang, etal, the geometric meaning of signed distance function with and application [J] in curved surface approaches, Journal of Computing and Information Science in Engineering, 2010,10 (4): 1-10)

(6) Liao's equality has been set up point based on cutting apart approximation theory to complex-curved distance model, and has evaluated complex-curved profile error by genetic algorithm.(referring to Liao Ping, based on genetic algorithm with cut apart the complex-curved profile error of approximatioss accurate Calculation [J]. mechanical engineering journal, 2010,46 (10): 1~7)

(7) set up complex-curved thick coupling location model based on rigid motion theory under Euclidean Norm (referring to Xu Yi, Li Zexiang, the matching detection new method [J] of free form surface. Harbin Institute of Technology's journal, 2010,42 (1): 106~108), (Liu Yuanpeng, Liu Jing etc., complex-curved measurement data optimum matching problem research [J], China Mechanical Engineering, 16 (12): 1080~1082)

(8) how Wen Xiulan has studied and to have solved flatness and Linear Error (referring to Xiulan Wen by genetic algorithm, Aiguo Song, for the improved genetic algorithm of one [J] the .International Journal of Machine Tools & Manufacture of Evaluation plane degree and straightness error, 2003,43 (11): 1157~1162)

Inventor is realizing in process of the present invention, finds at least to exist in prior art following shortcoming and defect:

Above-mentioned achievement in research major part has only been set up the error evaluation of simple profile, though some has completed complex-curved error evaluation, but the error amount assessing is far longer than given tolerance range, causes workpiece mistake rate higher, has improved production cost.

Summary of the invention

The invention provides a kind of method of determining complex-curved profile error, the present invention has reduced workpiece mistake rate and production cost, described below:

A method of determining complex-curved profile error, said method comprising the steps of:

(1) the survey termination of contact three-dimensional gauge head contacts the digital control system of sending trigger pip after surface of the work and send to by signal receiving device machining center, the current coordinate position latch that described contact three-dimensional gauge head is surveyed end by described digital control system is got off, and sends described current coordinate position to data processing terminal;

(2) described data processing terminal is set up complex-curved profile error model by described current coordinate position;

(3) described data processing terminal, by measuring the thick coupling of coordinate system and design coordinate system, obtains the pose initial parameter of ideal surface;

(4) described data processing terminal obtains complex-curved profile tolerance initial error by described pose initial parameter and described complex-curved profile error model;

(5) described data processing terminal is optimized described pose initial parameter, and again execution step (4), obtain complex-curved profile tolerance error current, until the difference of described complex-curved profile tolerance error current and described complex-curved profile tolerance initial error is while being less than threshold value, using described complex-curved profile tolerance error current as complex-curved profile tolerance final error and output.

Described data processing terminal is set up complex-curved profile error model by described current coordinate position and is specifically comprised:

1) described data processing terminal obtains complex-curved profile tolerance tolerance range;

| ( p - q ) · n q | ≤ t 2

Wherein, p point is measurement point, and q point is mathematical point, n qthe normal vector that represents q on ideal surface, t is tolerance value;

2) described data processing terminal is set up described complex-curved profile error model by described complex-curved profile tolerance tolerance range;

d ( p , S ) = x p x n q - x q x n q + y p y n q - y q y n q + z p z n q - z q z n q

+ y p z n q - z p y n q z p x n q - x p z n q x p y n q - y p x n q - x n q - y n q - z n q T α β γ δ x δ y δ z

+ α β γ δ x δ y δ z T - 2 ( z q z n q + y q y n q ) x q y n q + y q x n q x q z n q + z q x n q 0 - z n q y n q x q y n q + y q x n q - 2 ( x q x n q + z q z n q ) y q z n q + z q y n q z n q 0 - x n q x q z n q + z q x n q y q z n q + z q y n q - 2 ( y q y n q + x q x n q ) - y n q x n q 0 0 z n q - y n q 0 0 0 - z n q 0 x n q 0 0 0 y n q - x n q 0 0 0 0 α β γ δ x δ y δ z

Wherein, any one complex-curved around x, y, the angle of z axle rotation is respectively α, beta, gamma; Along x, y, the translational movement of z axle is δ x, δ y, δ z; The coordinate of measurement point p is (x p, y p, z p), the coordinate of mathematical point q is (x q, y q, z q), normal vector n qcoordinate be

Described data processing terminal obtains complex-curved profile tolerance initial error by described pose initial parameter and described complex-curved profile error model and is specially:

Described data processing terminal is brought into described pose initial parameter in described complex-curved profile error model and gets the ultimate range of described measurement point p to described ideal surface, using the ultimate range of 2 times as described complex-curved profile tolerance initial error.

The beneficial effect of technical scheme provided by the invention is: this method, by pose initial parameter is optimized, is obtained complex-curved profile tolerance error current; In the time that difference is less than threshold value, the also output using complex-curved profile tolerance error current as complex-curved profile tolerance final error; The final error getting by this method approaches given tolerance range, has reduced workpiece mistake rate and production cost, has improved accuracy of detection.

Brief description of the drawings

Fig. 1 is the structural representation of machining center provided by the invention;

Fig. 2 is the process flow diagram of a kind of method of determining complex-curved profile error provided by the invention.

Embodiment

For making the object, technical solutions and advantages of the present invention clearer, below in conjunction with accompanying drawing, embodiment of the present invention is described further in detail.

In order to reduce workpiece mistake rate and production cost, the embodiment of the present invention provides a kind of method of determining complex-curved profile error, referring to Fig. 1 and Fig. 2, described below:

101: the survey termination of contact three-dimensional gauge head contacts the digital control system of sending trigger pip after surface of the work and send to by signal receiving device machining center, the current coordinate position latch that contact three-dimensional gauge head is surveyed end by digital control system is got off, and sends current coordinate position to data processing terminal;

As shown in Figure 1, contact three-dimensional gauge head is arranged on machine tool chief axis by the taper shank of standard.The connection of the two is very easy, but it should be noted the centering of gauge head and main-shaft axis.The major function of contact three-dimensional gauge head is to adopt a little at surface of the work, after the survey termination of contact three-dimensional gauge head contacts surface of the work, sends trigger pip and sends signal receiving device to.Signal receiving device is arranged on the appropriate location of machining center, and the trigger pip of sending for receiving contact three-dimensional gauge head carries out trigger pip being sent to after necessary signal condition the digital control system of machining center.Owing to carrying out signal transmission with infrared mode between contact three-dimensional gauge head and signal receiving device, so need to consider not stop in the course of work of machining center the signal transmission between Contacting three coordinates measurement gauge head and receiving trap when signal receiving device is installed.Signal receiving device send to machining center digital control system be only contact three-dimensional gauge head touch surface of the work produce trigger pip, digital control system receives the servo-drive system stop motion at controlled working center immediately after trigger pip, and the current coordinate position latch that the contact three-dimensional gauge head that is arranged on machining center main shaft place is surveyed end is got off.Digital control system receives the information of the point of the needs measurement transmitting by communication line on the one hand, sends the geological information of the surface of the work point of having measured to data processing terminal on the one hand.

102: data processing terminal is set up complex-curved profile error model by current coordinate position;

Wherein, this step specifically comprises:

1) data processing terminal obtains complex-curved profile tolerance tolerance range;

Refer to that for complex-curved key element profile tolerance tolerance range diameter equals tolerance value t, the centre of sphere is positioned at the region that two enveloping surfaces of a series of balls on ideal surface S limit, so find, in space, some the distance function to tested curved surface key element is most important, Euclidean distance d (the p of nearest any of what point can be similar to the distance of ideal surface S be expressed as this point of distance on point to curved surface, S)

d(p,S)=‖p-q‖ 2

Wherein, the point that q point is upper for ideal surface S and p point is nearest, p point is measurement point, q point is mathematical point.

Further, by d (p, S) linear parsing at q point place

d(p,S)≈d(q)+▽d(q) T·(p-q)

Wherein d (q) represents the distance of some q to ideal surface S, ▽ d (q) trepresent the Grad of d (p, S) function at q.The eikonal equation(eikonal equation of introducing) be the nonlinear partial differential equation running into during ripple is propagated, there is following form | ▽ u (x) |=F (x), x ∈ Ω, Ω is R nin (n ties up real number space), be communicated with open region, F (x) is the Positive Function of input, and ▽ refers to gradient, || represent European norm, represent the border of Ω.Physically, what the solution u of equation (x) represented is the bee-line from the border of Ω to Ω internal point, and in the time of F=1, what this separated representative is the directed distance from Ω border to Ω inside.If d is (p, S) be u (x), can be using the outside ideal surface border as being communicated with open region Ω, ▽ d (q, S)=1, d (q, S)=0, the solution of d (p, S) is the outer bee-line that a bit arrives ideal surface S of ideal surface so, note ▽ d (q, S)=1=n q, n qthe normal vector that represents q on ideal surface, can be expressed as d (p, S)=(p-q) n to the signed distance function of ideal surface S so q.If measurement point p is in tolerance range, this measurement point p should meet t is as implied above, representative be tolerance value, this formula can be defined as to complex-curved profile tolerance tolerance range.

2) data processing terminal is set up complex-curved profile error model by complex-curved profile tolerance tolerance range.

In the time that the pose of ideal surface S changes, while being rotated translation motion, the mathematical point q corresponding with ideal surface S external pelivimetry point p also can change.

If one is complex-curved around x, y, the angle of z axle rotation is respectively α, beta, gamma, along x, y, the translational movement of z axle is δ x, δ y, δ z, ω=[α, beta, gamma], v=[δ x, δ y, δ z], volume coordinate rotation matrix is;

R = cos α cos β cos β sin γ - sin β sin α sin β cos γ - cos α sin γ sin α sin β sin γ + cos α cos γ sin α cos β cos α sin β cos γ + sin α sin γ cos α sin β sin γ - sin α cos γ cos α cos β

Work as α, beta, gamma is enough little, utilizes Taylor expansion formula to omit higher order term, can obtain the approximate matrix R of R *

R * = 1 γ - β - γ 1 α β - α 1

For convenience of calculating, introduce matrix, that is:

ω ^ = R * - I = 0 γ - β γ 0 α β - α 0 , Wherein I is unit matrix

R *can be reduced to R * = 1 γ - β - γ 1 α β - α 1 = I + ω ^

Provable: x is vectorial arbitrarily

If curved surface just carries out small variation, translation and rotation can be expressed as signed distance function when all very little d ( p , S ) = [ p - ( I + ω ^ ) · q - v ] · ( I + ω ^ ) · n q .

Its representative be point to the signed distance function of movement curved surface, utilize calculus of differences can obtain the single order differential term of signed distance function, the paper of mentioning in itself and technical background (5) is made comparisons:

Δd ( p , S ) = ( p - ( I + ω ^ ) · q - v ) · ( I + ω ^ ) · n q - ( p - q ) · n q

= ( p - q - ω ^ q - v ) · ( I + ω ^ ) · n q - ( p - q ) · n q

= ( p - q ) · n q + ( p - q ) · ω ^ n q - ( ω ^ q + v ) · n q - ( ω ^ q + v ) · ω ^ n q - ( p - q ) · n q

In the model of paper, the line of p and q is and n qthe line segment paralleling, so Section 2 is zero in above formula, and because rotation translational movement is very little, item more than second order can be ignored, so above formula becomes

Δd ( p , S ) = - ( ω ^ q + v ) · n q

= - n q · v - ( q × n q ) · ω

This has just obtained and result same in paper, has proved the correctness of model, using the signed distance function with rotation translational movement as objective function.After formula is arranged, obtaining signed distance function differential expressions is:

d ( p , S ) = x p x n q - x q x n q + y p y n q - y q y n q + z p z n q - z q z n q

+ y p z n q - z p y n q z p x n q - x p z n q x p y n q - y p x n q - x n q - y n q - z n q T α β γ δ x δ y δ z

+ α β γ δ x δ y δ z T - 2 ( z q z n q + y q y n q ) x q y n q + y q x n q x q z n q + z q x n q 0 - z n q y n q x q y n q + y q x n q - 2 ( x q x n q + z q z n q ) y q z n q + z q y n q z n q 0 - x n q x q z n q + z q x n q y q z n q + z q y n q - 2 ( y q y n q + x q x n q ) - y n q x n q 0 0 z n q - y n q 0 0 0 - z n q 0 x n q 0 0 0 y n q - x n q 0 0 0 0 α β γ δ x δ y δ z

Mathematical model using above formula as Evaluating Profile Error of Complex Surfaces.

Wherein, the coordinate of measurement point p is (x p, y p, z p), the coordinate of mathematical point q is (x q, y q, z q), normal vector n qcoordinate be

103: data processing terminal, by measuring the thick coupling of coordinate system and design coordinate system, obtains the pose initial parameter of ideal surface;

Because measurement point in model and mathematical point are under the same coordinate system, thus need two processes for the evaluation of complex-curved profile error, the one, need to align with design coordinate system measuring coordinate system, to eliminate the errors of bringing due to coordinate system difference; The 2nd, thus utilize optimized algorithm error model to be optimized to the result that obtains evaluation.

Align for coordinate system, i.e. thick coupling, take the coordinate matching method based on unique point, because matching process does not change the inherent geometric properties of measuring curved surface or ideal surface, so can regard rigid body as by measuring curved surface or ideal surface, rotary flat shifting movement is as rigid motion, rigid motion does not change the result of inner product and apposition computing, so the distance of measuring on curved surface or ideal surface between points can not change because the generation of rotation translation, so just can extract unique point, P, Q represents respectively to measure point set and theoretical point set, O (x lz, y lz, z lz) mid point of representation theory point set Q, O ' (x cz, y cz, z cz) represent that the mid point idiographic flow of measurement point set P is as follows:

1) get theoretical point set Q and the coordinate of measuring point set P mid point;

O ( x lz , y lz , z lz ) = Σ i = 1 n p i ( x , y , z ) n , O ′ = Σ i = 1 n q i ( x , y , z ) n

Wherein, p irepresent to measure the arbitrfary point in point set P, q iarbitrfary point in representation theory point set Q.

2) in whole theoretical point set Q and measurement point set P, the nearest point of detection range mid point is put p as First Characteristic respectively tz1, q tz1;

3) ergodic theory point set Q finds respectively with First Characteristic point distance point farthest and puts p as Second Characteristic with measurement point set P tz2, q tz2;

4) ergodic theory point set Q finds respectively with Second Characteristic point distance point farthest respectively as the 3rd unique point p with measurement point set P tz3, q tz3;

5) taking three unique points as Foundation local coordinate system:

The X-axis of two local coordinate systems is respectively: y-axis is e 2 = e 1 × ( p tz 1 - p tz 2 ) | e 1 × ( p tz 1 - p tz 2 ) | , e 2 ′ = e 2 × ( q tz 1 - q tz 3 ) | e 2 × ( q tz 1 - q tz 3 ) | , Z axis is e 3=e 1× e 2, e 3 ′ = e 1 ′ × e 2 ′ .

6) pass through can determine the rotation matrix that is tied to design coordinate system from measuring coordinate, then pass through K=p tz1-q tz1r, can determine translation vector.

Pass through above-mentioned steps, can complete and measure coordinate system and design aliging of coordinate system, after thick coupling, measuring point set P and theoretical point set Q is in little deviation range, meet the requirement to parameter area in model completely, in meticulous coupling below, just carry out error evaluation using the complex-curved profile error model of setting up as objective function above.The process of above-mentioned thick coupling, for conventionally known to one of skill in the art, when specific implementation, can also adopt other thick matching process, and the embodiment of the present invention does not limit this.

104: data processing terminal obtains complex-curved profile tolerance initial error by pose initial parameter and complex-curved profile error model;

This step is specially: pose initial parameter is brought in complex-curved profile error model and gets the ultimate range of measurement point p to ideal surface S, using the ultimate range of 2 times as complex-curved profile tolerance initial error.

105: data processing terminal is optimized pose initial parameter, and again perform step 104, obtain complex-curved profile tolerance error current, until the difference of complex-curved profile tolerance error current and complex-curved profile tolerance initial error is while being less than threshold value, using complex-curved profile tolerance error current as complex-curved profile tolerance final error and output.

Wherein, the embodiment of the present invention gets complex-curved profile tolerance final error by continuous iteration, and threshold value is set according to the needs in practical application, and when specific implementation, the embodiment of the present invention does not limit this.

Wherein, data processing terminal is optimized as conventionally known to one of skill in the art pose initial parameter, can adopt multiple optimization method to realize, and the embodiment of the present invention does not repeat at this.

Verify the feasibility of a kind of method of determining complex-curved profile error that the embodiment of the present invention provides below with concrete test, described below:

In order to verify the validity of above-mentioned model, design is experiment as follows, taking the worm tooth-surface often used in producing and ellipsoid as example, describes.

Two kinds of curved surfaces are carried out respectively to 5 groups, amounted to 10 groups of experiments, respectively two measurement points have been added to (0 of Normal Distribution, 0.004), (0,0.005), (0,0.006), (0,0.007), (0,0.008) stochastic error is similar to mismachining tolerance, evaluates respectively by this method and least square method, as shown in the table:

The error amount that this method assesses is significantly less than the error amount that least square method assesses, and so just can approach actual tolerance range, obtains real error amount, obtains high accuracy of detection, and the manufacturer of making a living provides real examining report, reduces the useless rate of mistake and production cost.

In sum, the embodiment of the present invention provides a kind of method of determining complex-curved profile error, and this method, by pose initial parameter is optimized, is obtained complex-curved profile tolerance error current; In the time that difference is less than threshold value, the also output using complex-curved profile tolerance error current as complex-curved profile tolerance final error; The final error getting by this method approaches given tolerance range, has reduced workpiece mistake rate and production cost, has improved accuracy of detection.

It will be appreciated by those skilled in the art that accompanying drawing is the schematic diagram of a preferred embodiment, the invention described above embodiment sequence number, just to describing, does not represent the quality of embodiment.

The foregoing is only preferred embodiment of the present invention, in order to limit the present invention, within the spirit and principles in the present invention not all, any amendment of doing, be equal to replacement, improvement etc., within all should being included in protection scope of the present invention.

Claims (2)

1. a method of determining complex-curved profile error, is characterized in that, said method comprising the steps of:
(1) the survey termination of contact three-dimensional gauge head contacts the digital control system of sending trigger pip after surface of the work and send to by signal receiving device machining center, the current coordinate position latch that described contact three-dimensional gauge head is surveyed end by described digital control system is got off, and sends described current coordinate position to data processing terminal;
(2) described data processing terminal is set up complex-curved profile error model by described current coordinate position;
(3) described data processing terminal, by measuring the thick coupling of coordinate system and design coordinate system, obtains the pose initial parameter of ideal surface;
(4) described data processing terminal obtains complex-curved profile tolerance initial error by described pose initial parameter and described complex-curved profile error model;
(5) described data processing terminal is optimized described pose initial parameter, and again execution step (4), obtain complex-curved profile tolerance error current, until the difference of described complex-curved profile tolerance error current and described complex-curved profile tolerance initial error is while being less than threshold value, using described complex-curved profile tolerance error current as complex-curved profile tolerance final error and output;
Described data processing terminal is set up complex-curved profile error model by described current coordinate position and is specifically comprised:
1) described data processing terminal obtains complex-curved profile tolerance tolerance range;
| ( p - q ) n q | ≤ t 2
Wherein, p point is measurement point, and q point is mathematical point, n qthe normal vector that represents q on ideal surface, t is tolerance value;
2) described data processing terminal is set up described complex-curved profile error model by described complex-curved profile tolerance tolerance range;
Wherein, any one complex-curved around x, y, the angle of z axle rotation is respectively α, beta, gamma; Along x, y, the translational movement of z axle is δ x, δ y, δ z; The coordinate of measurement point p is (x p, y p, z p), the coordinate of mathematical point q is (x q, y q, z q), normal vector n qcoordinate be .
2. a kind of method of determining complex-curved profile error according to claim 1, it is characterized in that, described data processing terminal obtains complex-curved profile tolerance initial error by described pose initial parameter and described complex-curved profile error model and is specially:
Described data processing terminal is brought into described pose initial parameter in described complex-curved profile error model and gets the ultimate range of described measurement point p to described ideal surface, using the ultimate range of 2 times as described complex-curved profile tolerance initial error.
CN201210236183.8A 2012-07-09 2012-07-09 Method for determining complex curved surface profile error CN102809364B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210236183.8A CN102809364B (en) 2012-07-09 2012-07-09 Method for determining complex curved surface profile error

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210236183.8A CN102809364B (en) 2012-07-09 2012-07-09 Method for determining complex curved surface profile error

Publications (2)

Publication Number Publication Date
CN102809364A CN102809364A (en) 2012-12-05
CN102809364B true CN102809364B (en) 2014-10-22

Family

ID=47233144

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210236183.8A CN102809364B (en) 2012-07-09 2012-07-09 Method for determining complex curved surface profile error

Country Status (1)

Country Link
CN (1) CN102809364B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104345687B (en) * 2013-08-09 2017-10-20 中国科学院沈阳计算技术研究所有限公司 A kind of tool sweep volume modeling method
CN103591912B (en) * 2013-11-11 2016-03-16 沈阳黎明航空发动机(集团)有限责任公司 The measuring method of the special-shaped feature locations degree of a kind of ring-shaped work pieces
CN103617369B (en) * 2013-12-09 2016-03-30 国核(北京)科学技术研究院有限公司 Ultra-large type curved surface parts machining precision assessment method and equipment
CN104089599B (en) * 2014-07-04 2017-01-25 北京工业大学 Quasi morphological filtering method for extracting two-dimensional contour in contact measuring head measurement
CN105651216B (en) * 2014-11-10 2018-08-03 沈阳黎明航空发动机(集团)有限责任公司 A kind of no positioning reference plane blade parts coordinate system measurement method
CN104515478B (en) * 2014-12-11 2017-08-01 华中科技大学 A kind of automatic method for three-dimensional measurement of high-precision blade of aviation engine and system
CN104634294A (en) * 2015-02-04 2015-05-20 天津大学 Method for detecting and evaluating geometric error of grooved pulley of curved groove
CN104750977B (en) * 2015-03-10 2018-03-23 北京理工大学 A kind of method and device of complex position degree error evaluation
CN105468834A (en) * 2015-11-20 2016-04-06 南昌航空大学 Assembly quality control method based on measurement point error ellipsoids
CN105910570A (en) * 2015-12-31 2016-08-31 广东冠邦科技有限公司 Irregular variable cross section roller hole type curved surface error measurement method
CN107084698B (en) * 2017-04-26 2020-03-31 中车唐山机车车辆有限公司 Rail vehicle wheel detection system and detection method

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101413793A (en) * 2008-11-27 2009-04-22 天津大学 Method for assessing composite space type face geometric error
CN102168965A (en) * 2010-12-21 2011-08-31 北京信息科技大学 A method and a device for obtaining a profile tolerance error of a helical surface

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10277888A (en) * 1997-04-01 1998-10-20 Toshiba Corp Processing data modifying method and device

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101413793A (en) * 2008-11-27 2009-04-22 天津大学 Method for assessing composite space type face geometric error
CN102168965A (en) * 2010-12-21 2011-08-31 北京信息科技大学 A method and a device for obtaining a profile tolerance error of a helical surface

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
JP特开平10-277888A 1998.10.20
何改云 等.复合位置度误差的最小条件解析算法.《计量技术》.2004,6-9.
何改云.形位误差的逼近原理及算法研究.《 中国博士学位论文全文数据库》.2007,全文.
复合位置度误差的最小条件解析算法;何改云 等;《计量技术》;20041231;6-9 *
形位误差的逼近原理及算法研究;何改云;《 中国博士学位论文全文数据库》;20070515;全文 *

Also Published As

Publication number Publication date
CN102809364A (en) 2012-12-05

Similar Documents

Publication Publication Date Title
Sładek Coordinate Metrology
CN103363949B (en) Mixed measurement analysis method for satellite antenna
US6922599B2 (en) System and method for producing an assembly by directly implementing three-dimensional computer-aided design component definitions
Kunzmann et al. A uniform concept for calibration, acceptance test, and periodic inspection of coordinate measuring machines using reference objects
Yu et al. Adaptive sampling method for inspection planning on CMM for free-form surfaces
CN103256916B (en) Evaluation method of part flatness error based on minimum area
CN103878641B (en) The rotating shaft geometric error discrimination method that a kind of five-axle number control machine tool is general
CN101957175B (en) Three-point micro-plane-based normal detection method
Mutilba et al. Traceability of on-machine tool measurement: a review
CN101342664B (en) On-machine quality detecting system for complex space type surface
CN103328154B (en) Error measure device and error measure method
CN100462677C (en) Three-coordinate measuring machine-tool error compensation system and method
CN106354094B (en) Lathe slave laser scanning coordinate scaling method based on space criteria ball
CN102944197B (en) A kind of method for detecting precision of five-spindle machining center of double-rotary table structure
CN102785129A (en) On-line detection method of curved surface machining precision of complex part
CN101913103A (en) Method for measuring angular errors of rotating table of numerical control machine
CN101758422B (en) Detection analysis device for technical indexes of numerical control device
Fu et al. Product-of-exponential formulas for precision enhancement of five-axis machine tools via geometric error modeling and compensation
Xiang et al. Using a double ball bar to identify position-independent geometric errors on the rotary axes of five-axis machine tools
CN103389038A (en) Targeting multi-station measuring method for detecting geometric accuracy of numerical control machine tool through laser tracker
CN104308657B (en) A kind of rotary axis of machine tool geometric error six measured based on ball bar encloses discrimination method
CN107186548B (en) A kind of five-axle number control machine tool rotating shaft geometric error detection method
CN102721393B (en) On-site self-calibration method for measurement system error of precise rotary table
CN101907455B (en) Three-coordinate measuring method for positioning pin not vertical to reference surface
CN102785128B (en) The part processing precision on-line detecting system of NC Machine lathe and detection method

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model