CN102183205A - Method for matching optimal assembly poses of large-sized parts - Google Patents

Method for matching optimal assembly poses of large-sized parts Download PDF

Info

Publication number
CN102183205A
CN102183205A CN 201110053699 CN201110053699A CN102183205A CN 102183205 A CN102183205 A CN 102183205A CN 201110053699 CN201110053699 CN 201110053699 CN 201110053699 A CN201110053699 A CN 201110053699A CN 102183205 A CN102183205 A CN 102183205A
Authority
CN
China
Prior art keywords
pose
assembling
parts
assembly
detection
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
CN 201110053699
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.)
Beihang University
Original Assignee
Beihang University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beihang University filed Critical Beihang University
Priority to CN 201110053699 priority Critical patent/CN102183205A/en
Publication of CN102183205A publication Critical patent/CN102183205A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Automatic Assembly (AREA)

Abstract

The invention discloses a measurement data processing method for pose measurement and adjustment in digital assembly of large-sized parts, comprising the following steps: comparing the current pose of a digital theory model matching assembly part according to actually measured data; and calculating actual optimal target assembly poses according to the current pose. The position of each critical characteristic on the part and corresponding detection points of each assembly detection item are measured by using a laser tracking instrument. Measurement data after processed by the method can provide accurate assembly position and pose data and positioning adjustment parameters for a large-sized part digital flexible assemble system, and the optimal assembly effect of parts is achieved by matching real-time measurement.

Description

The method of the best assembling of a kind of large parts pose coupling
Technical field:
The present invention relates to the method for the best assembling pose coupling in a kind of large parts digitizing alignment measurement process.
Background technology:
Measuring auxiliary mounting technology is to realize large product digitizing assembling, one of effective way of control assembly quality and key characteristic, utilize portable large-scale metrology equipment in the product assembling, each part pose to be carried out real-time follow-up, check the position of key feature, to reach the technical requirement of product assembly quality.The accurate adjustment of pose is to guarantee one of focus of product assembly precision with control in the large parts assembling, and especially the effect when parts butt joint, product final assemble is more outstanding.
Assembling pose coupling is a key link of Measurement and Data Processing in the large parts digitizing assembling, being intended to provides the accurate pose data of assembly so that its accurate location for the enforcement of measuring auxiliary mounting technology, and provide optimum target assembling pose to guarantee the assembling accuracy, reach best practical set effect.
On the current pose of assembly and optimum target pose definite, general processing mode is simply to utilize each critical size and form and position tolerance whether to be in the tolerances in design scope to judge, is in fact rough a kind of finding the solution in large component assembling.Therefore, to realize under this mode that high-precision assembly quality still needs to adopt the raising Allowance Design to require to control, this means that the manufacturing accuracy of parts has requirements at the higher level, all not ideal enough on effect and the economy.The present invention proposes finding the solution of can be in the tolerances in design allowed band accurate more object pose of realizing best assembling effect of method, adjusts parameter for the digitizing locating device provides better location.
Summary of the invention:
The objective of the invention is needs, a kind of comparatively desirable pose matching process of proposition at the deficiency and the auxiliary assembling and positioning of large parts digitizing of the general processing mode existence of the data in the assembling pose measurement.
The method of best assembling pose coupling in a kind of large parts assembling, it comprises that current pose coupling of parts and parts optimum target assembling pose mated for 2 stages:
Current pose matching stage comprises following basic step:
1) locus of fixed reference reference point in the mensuration erecting yard;
2) the reference point is fitted to benchmark assembling coordinate system, unify with global coordinate system in the digitizing assembling model;
3) locus of each position reference point on the mensuration removable assembly to be positioned;
4) nominal position of each position reference point in the contrast digitizing assembling model uses current pose matching algorithm to ask for six independent parameters of the current spatial pose of assembly;
Optimum target assembling pose matching stage comprises following basic step:
5) measure on the oriented benchmark assembly position with the corresponding measuring point of each assembling accuracy detection, fit to corresponding geometric element according to the type of reference data geometric element;
6) measure the position of assembling the corresponding measuring point of accuracy detection on the removable assembly to be positioned with each, and measuring point data is transformed under the local parts coordinate system of assembly;
7) each point position of connecting firmly with assembly to be positioned of supposition is along with assembly moves by rule in the certain space scope, the accuracy error of each detection when the calculating assembly is in the series of discrete position and attitude, and calculate the comprehensive overall accuracy Z-factor of multiple goal to the single goal mapping;
8) calculate each individual event accuracy error and overall accuracy Z-factor and whether satisfy the setting requirement, otherwise repeating step 7, until the spatial pose of the removable assembly to be positioned that obtains to meet the demands, this pose is the actual optimum target assembling pose of assembly to be positioned.
The measurement of pose reference point and each detection measuring point all adopts laser tracker to carry out the data acquisition operations on described measurement erecting yard reference point, the parts.Laser tracker is positioned over the field position that can observe each common point and corresponding measuring point.
The step that the reference measurement point is fitted to the frame of reference is: earlier each benchmark measuring point is fitted to reference plane XOY, the projection of selected wherein two somes straight line on XOY plane is as X-direction, selected a bit in the projection of XOY plane as true origin O, determine the basic parameter of match coordinate system thus, and this original match coordinate system is carried out the certain space conversion as overall reference for assembling coordinate system.
The pose of parts can be represented by the locus of this parts coordinate system under overall situation assembling coordinate system, its hexa-atomic group of being made of 6 independent parameters
Figure BDA0000049056230000021
The expression parts,, rotate around the z axle around y axle rotation θ degree successively around the x of reference frame axle rotation ψ degree from initial pose to current pose
Figure BDA0000049056230000022
Degree, at last at x, y, each translation p of z direction x, p y, p zObtain, carry out spatial pose and these 6 independent parameter argument tables are shown when calculating 4 * 4 spatial alternation matrix T and participate in computing again, wherein R is the rotation matrix component, and P is a translational component:
T = R P 0 1
Figure BDA0000049056230000024
P = p x p y p z
Described current pose matching algorithm, the coupling of measuring the locus of parts B under global coordinate system and the nominal position of digitizing master pattern mid point with retrain the form that all is expressed as least square, adopt the model of multiple-objection optimization, constraint distributes weights to a coupling and each, the optimization aim function as:
min{sum(ε i)}
ε wherein iThe matching error of representing each pose reference point, represent with the range deviation that physical location of putting and model are in this following nominal position of spatial pose state:
Constraint condition: E I_low<ε i<E I_up, E I_low, E I_upBe each reference point of regulation and the position deviation bound of theoretical value.{ P I-0Be nominal position point group on the entire assembly model, coordinate figure is with respect to local parts coordinate system, { P I-1Being location point group on the actual measurement parts, coordinate figure is with respect to assembling coordinate system, T 1For waiting to ask the parts spatial pose matrix of representing with 6 independent parameters, on this basis, the optimization aim function further is refined as:
min{sum(ε i)}=min{sum(P i-1-T 1P i-0)}
Each measuring point name and measured data all are expressed as the form (x, y, z, 1) of homogeneous coordinates TParticipate in calculating, the computing formula of each point coupling deviation is:
ϵ i = ϵ ix ϵ iy ϵ iz 0 = x i - 1 y i - 1 z i - 1 1 - T 1 x i - 0 y i - 0 z i - 0 1
Described current pose matching algorithm uses heuristic particle cluster algorithm to realize finding the solution, and on behalf of a possible spatial pose, each particle separate, 6 dimension x of particle i=(x I1, x I2..., x ID) respectively expression determine 6 independent parameter amounts of pose
Figure BDA0000049056230000032
Each particle i is included as the position vector x of a D dimension i=(x I1, x I2..., x ID) and velocity vector v i=(v I1, v I2..., v ID), during particle i search solution space, preserve its optimum that searches experience position p i=(p I1, p I2..., p ID), at the beginning of each iteration, according to self inertia and the optimum experience of experience and colony position p g=(p G1, p G2..., p GD) adjust oneself velocity vector to adjust self-position, until the optimum solution that obtains to satisfy condition.R is to position renewal the time, and a coefficient in that speed adds previously is constraint factor.c 1, c 2For positive constant, be referred to as speedup factor.ξ, η is for all there being the random number of distribution in [0 1].W is the inertia weight factor.
Each particle position changes and is undertaken by following formula:
x id=x id+rv id
v id=wv id+c 1ξ(p id-x id)+c 2η(p gd-x id)
D is the dimension in the D dimension in the formula, and this value is 6 in the pose coupling.
Described optimum target assembling pose matching stage, each assembles the accuracy detection and comprises general size, form and position tolerance and the detection content that is provided with by special matching requirements.Wherein must earlier the reference measurement data be fitted to the standard basis geometric element by the least square rule, comprise three types of point, line, surface the detection that has reference to require.
Described optimum target assembling pose matching stage, basic procedure is will respectively to detect corner group coordinate figure { P on the removable parts to be positioned earlier kBe transformed into { P under the local coordinate system of parts from global coordinate system KL, these measuring points that convert then move to global position coordinate figure { P different pose state T under by rule with these parts in the space KS, then in conjunction with benchmark measuring point group numerical value { P KDThese parts be in the space this detection assessment result that may assemble under the pose state predict that its computing formula is:
ε I-k=f I-k(P kD,P kS)
f I-kFor calculating the function universal expression form of detection numerical result, determine by type of detection.The accuracy error of each detection when calculating parts to be positioned and being in the series of discrete position, and calculate the comprehensive overall accuracy accuracy error of multiple goal to the single goal mapping, find the spatial pose Ta of overall accuracy accuracy error minimum to assemble pose as optimum target.
Described best assembling pose matching process, algorithm model is for all being expressed as the matching constraint of positioning reference parts and parts assembling accuracy detection of each key characteristic under global coordinate system to be positioned the form of detection deviate, adopt the Optimization Model of multiple goal, each detection coupling and each constraint are distributed weights W to the unified mapping of single goal kAfter ask the weighted comprehensive deviation,, the optimization aim function as:
min { E } = min { Σ k = 1 m W k ϵ k / Σ k = 1 m W k }
= min { Σ k = 1 m W k f k ( P kD , P kS ) / Σ k = 1 m W k }
= min { Σ k = 1 m W k f k ( P kD , TT 1 - 1 P k ) / Σ k = 1 m W k }
Constraint condition: E K_low<ε I-k<E K_upBe each matching requirements detection normal value bound,
ε wherein kThe matching error of representing k item assembling detection, T 1Result of calculation for the rapid current pose matching stage of previous step.The assign weight method of suing for peace behind the coefficient of employing is transformed into the value that belongs under the same dimension with all kinds of deviations and compares differentiating the size of whole synthesis assembling accuracy precision, and synthesis precision is high more, and then to be considered as the assembly quality effect good more.
Described optimum target assembling pose matching algorithm equally with current pose matching algorithm uses heuristic particle cluster algorithm to realize finding the solution, and difference is that the parameters of bounding algorithm rule is provided with different.
The invention has the advantages that: the flow process and the data processing method that 1) provide detailed large parts to assemble pose measurement; 2) finding the solution of can be in the tolerances in design allowed band accurate more object pose of realizing best assembling effect adjusted parameter for the digitizing locating device provides better location; 3) by once or the data of minority pose measurement just can get access to comparatively ideal pose and adjust data, make assembling reach the effect of actual the best;
Description of drawings
Fig. 1 is the current pose matching algorithm flow process according to embodiment of the present invention;
Fig. 2 is the optimum target assembling pose matching algorithm flow process according to embodiment of the present invention;
Fig. 3 is the scene according to the invention process large parts alignment measurement example;
Fig. 4 is that the corresponding with it measurement of correlation point of cylindrical compartment section parts 1 to be positioned among the embodiment distributes;
Fig. 5 is the corresponding with it relevant measuring point distribution plan of oriented middle part interface frame, mating frame parts 3 among the embodiment;
Fig. 6 is the corresponding with it relevant measuring point distribution plan of oriented conical section parts 2 among the embodiment;
Among the figure: cylindrical compartment section parts 1, circular cone cabin section parts 2, interface frame, mating frame 3, supporting tool 4, laser tracker 5, cylindrical compartment fragment position hole measuring point group 6, cylindrical compartment section and circular cone cabin section interface depth of parallelism measuring point group 7, cylindrical compartment section and circular cone cabin section outside surface right alignment measuring point group 8, assembling coordinate system benchmark measuring point group 9.
Embodiment:
The method of best assembling pose coupling in the large parts assembling, it comprises that current pose coupling of parts and parts optimum target assembling pose mated for 2 stages:
Current pose matching stage comprises following basic step:
1) locus of fixed reference reference point in the mensuration erecting yard;
2) the reference point is fitted to benchmark assembling coordinate system, unify with global coordinate system in the digitizing assembling model;
3) locus of each position reference point on the mensuration removable assembly to be positioned;
4) nominal position of each position reference point in the contrast digitizing assembling model uses current pose matching algorithm to ask for six independent parameters of the current spatial pose of assembly;
Optimum target assembling pose matching stage comprises following basic step:
5) measure on the oriented benchmark assembly position with the corresponding measuring point of each assembling accuracy detection, fit to corresponding geometric element according to the type of reference data geometric element;
6) measure the position of assembling the corresponding measuring point of accuracy detection on the removable assembly to be positioned with each, and measuring point data is transformed under the local parts coordinate system of assembly;
7) each point position of connecting firmly with assembly to be positioned of supposition is along with assembly moves by rule in the certain space scope, the accuracy error of each detection when the calculating assembly is in the series of discrete position and attitude, and calculate the comprehensive overall accuracy Z-factor of multiple goal to the single goal mapping;
8) calculate each individual event accuracy error and overall accuracy Z-factor and whether satisfy the setting requirement, otherwise repeating step 7, until the spatial pose of the removable assembly to be positioned that obtains to meet the demands, this pose is the actual optimum target assembling pose of assembly to be positioned.
At first with reference to figure 3 Fig. 5, disclosure example relates to and is used for method and the process that large parts digitizing assembling establishing criteria model and measured data are carried out the pose coupling.As used herein, term " part ", " parts ", " assembly " are intended to comprise the assembly of independent part, the assembling of a plurality of part.In an illustrated embodiment, parts comprise the big assembly that is used to build aircraft, and this method and process can be used to and be applicable to the part assemblies of various other types of widespread adoption.
As shown in Figure 3, the Measurement and Data Processing system at the large parts digitizing flexible assembly scene that one embodiment of the present of invention are used for the best assembling pose matching process that proposes according to the present invention, the robotization that this conjunction measuring system of system, posture adjustment control system and digitizing locating device are finished a plurality of parts is located and is docked.The situation of the butt joint of 2 large-scale cabin sections in certain the satellite agent structure that relates among the embodiment, cabin section 1 is cylindrical, and cabin section 2 is conical, and the cabin is intersegmental realizes being located by connecting by interface frame, mating frame 3, before being connected and fixed, need each cabin section is navigated under the desirable assembling pose state.Its middle deck section 2 has supported and has navigated on the frock, and cabin section 1 need be carried out pose measurement and is located by connecting on the interface frame, mating frame 3 that fixedly mounts cabin section 1 by optimum target assembling pose.
The measurement of pose reference point and each detection measuring point all adopts laser tracker to measure on described measurement erecting yard reference point, the parts.For guaranteeing measuring accuracy, adopt the scheme of many laser tracker combined measurements, each instrument is installed on the field position that should be able to observe each common point and corresponding measuring point mutually by measuring programme.Laser tracker adopts Switzerland Leica LTD640 laser tracker in the example, and the measurement radius of this instrument can reach 40m, and resolution is 0.001mm, and measuring accuracy can reach 15 μ+5 μ/m in the total travel scope.
The step that in the example reference measurement point is fitted to the frame of reference is: 4 fixed reference measuring points (some group 9) on the frock are fitted to reference plane XOY, the projection of selected wherein adjacent 2 somes straight line on XOY plane is as X-direction, selected wherein 1 projection in XOY plane is as true origin O, determine the basic parameter of match coordinate system thus, and this original match coordinate system is carried out spatial alternation arrive this four point geometry center as overall reference for assembling coordinate system O-XYZ, all raw measurement datas are all unified to be transformed into based under this overall situation assembling coordinate system, and assemble aliging of coordinate system in realization and the digitizing master pattern.
Set up global coordinate system, will make measuring system, pose adjust software systems, servo-control system and digital detent mechanism and under unified assembling coordinate system, work exactly, make the data of transmitting each other under same coordinate system.Each steady arm of control system cooperative motion under unified coordinate system is realized the adjustment of stressless large parts attitude.
The pose of parts 1 can be represented by the locus of this parts rigid body coordinate system under overall situation assembling coordinate system, its hexa-atomic group of being made of 6 independent parameters
Figure BDA0000049056230000061
Form, the expression parts,, rotate around the z axle around y axle rotation θ degree successively around the x of reference frame axle rotation ψ degree from initial pose to current pose
Figure BDA0000049056230000062
Degree, at last at x, y, each translation p of z direction x, p y, p zObtain.
Carry out spatial pose coupling and these 6 independent parameter argument tables are shown when calculating 4 * 4 spatial alternation matrix T and participate in computing again, wherein R is the rotation matrix component, and P is a translational component, and expression formula is as follows:
T = R P 0 1 ,
Figure BDA0000049056230000064
P = p x p y p z ,
Described current pose matching algorithm, the coupling of measuring the nominal position of each reference point (some group 6) in the locus of parts 2 under global coordinate system and the digitizing master pattern with retrain the form that all is expressed as least square, adopt the model of multiple-objection optimization, constraint distributes weights to a coupling and each, the optimization aim function as:
min{sum(ε i)}
ε wherein iThe matching error of representing each pose reference point, represent with the range deviation that physical location of putting and model are in this following nominal position of spatial pose state:
Constraint condition: E I_low<ε i<E I_up, E I_low, E I_upBe each reference point of regulation and the position deviation bound of theoretical value.
Among the embodiment, { P I-0Be nominal position point group on the model, coordinate figure is with respect to local parts coordinate system, { P I-1Being location point group on the actual measurement parts 1, coordinate figure is with respect to assembling coordinate system, T 1For waiting to ask the parts spatial pose of representing with 6 independent parameters.On this basis, the optimization aim function further is refined as:
min{sum(ε i)}=min{sum(P i-1-T 1P i-0)}
Each measuring point name and measured data all are expressed as the form (x, y, z, 1) of homogeneous coordinates TParticipate in calculating, the computing formula of each pose reference point coupling deviation is:
ϵ i = ϵ ix ϵ iy ϵ iz 0 = x i - 1 y i - 1 z i - 1 1 - T 1 x i - 0 y i - 0 z i - 0 1
Described current pose matching algorithm uses heuristic particle cluster algorithm to realize finding the solution, and on behalf of a possible spatial pose, each particle separate, 6 dimension x of particle i=(x I1, x I2..., x ID) respectively expression determine 6 independent parameter amounts of pose Each particle i is included as the position vector x of a D dimension i=(x I1, x I2..., x ID) and velocity vector v i=(v I1, v I2..., v ID), during particle i search solution space, preserve its optimum that searches experience position p i=(p I1, p I2..., p ID), at the beginning of each iteration, according to self inertia and the optimum experience of experience and colony position p g=(p G1, p G2..., p GD) adjust oneself velocity vector to adjust self-position, until the optimum solution that obtains to satisfy condition.R is to position renewal the time, and a coefficient in that speed adds previously is constraint factor.c 1, c 2For positive constant, be referred to as speedup factor.ξ, η is for all there being the random number of distribution in [0 1].W is the inertia weight factor.
Each particle position changes and is undertaken by following formula:
x id=x id+rv id
x id=wv id+c 1ξ(p id-x id)+c 2η(p gd-x id),
D is the dimension in the D dimension in the formula, and this value is 6 in the pose coupling.
Among the described embodiment, by using optimum target assembling pose matching process of the present invention, the utilization particle cluster algorithm carries out current pose matching algorithm and finds the solution, through after about 450 iteration, fitness function is restrained fully, the average coupling deviation for each position reference point of result of calculation is 0.108mm, and maximum deviation is 0.159mm; Adhere to specification.
Described optimum target assembling pose matching stage, each assembles the accuracy detection and comprises: size, Geometrical Tolerance Principle reach by the accurate requirement of special assembling.Form and position tolerance comprises form tolerance and position of related features two big classes, position of related features need be considered the reference element requirement, wherein must earlier the reference measurement data be fitted to the standard basis geometric element by the least square rule, comprise three types of point, line, surface the detection reference.
Described optimum target assembling pose matching stage, basic procedure is will respectively to detect corner group coordinate figure { P on the parts 1 to be positioned earlier kBe transformed into { P under the local coordinate system of parts from global coordinate system KL, these measuring points that convert then move to global position coordinate figure { P different pose state T under by rule with parts 1 in the space KS, then in conjunction with benchmark measuring point group numerical value { P KDCalculate parts 1 and be in the space this detection assessment result that may assemble under the pose state and predict that its computing formula is:
ε I-k=f I-k(P kD,P kS)
f I-kFor calculating the function universal expression form of detection numerical result, determine by type of detection.The accuracy error of each detection when calculating assembling parts 1 are in the series of discrete position, and calculate the comprehensive overall accuracy accuracy error coefficient of multiple goal to the single goal mapping, find the spatial pose Ta of overall accuracy accuracy error coefficient minimum to assemble pose by rule as optimum target.
Shown in Figure 6 as Fig. 4, assembling accuracy detection comprises among the described embodiment: the depth of parallelism of cylindrical compartment section 1 and butt end frame 3 interfaces, circular cone cabin section 2 and cylindrical compartment section 1 outside surface right alignment, the position degree of a plurality of location holes on the cylindrical compartment section 1.Wherein the corresponding measuring point of depth of parallelism detection is the some group 7 that is distributed on the interface; The corresponding measuring point of right alignment detection is some group 8; Each location hole measuring point is the some group 6 of each center, hole on the cylindrical compartment section 1.
Described best assembling pose matching process, algorithm model is for all being expressed as the matching constraint of positioning reference parts and parts assembling accuracy detection of each key characteristic under global coordinate system to be positioned the form of detection deviate, adopt the Optimization Model of multiple goal, each detection coupling and each constraint are distributed weights W to the unified mapping of single goal kAfter ask the weighted comprehensive deviation, the optimization aim function as:
min { E } = min { Σ k = 1 m W k ϵ k / Σ k = 1 m W k }
= min { Σ k = 1 m W k f k ( P kD , P kS ) / Σ k = 1 m W k }
= min { Σ k = 1 m W k f k ( P kD , TT 1 - 1 P k ) / Σ k = 1 m W k } ,
Constraint condition: E K_low<ε I-k<E K_upBe each matching requirements detection normal value bound,
ε wherein kThe matching error of representing k item assembling detection, T 1Result of calculation for current pose matching stage.The assign weight method of suing for peace behind the coefficient of employing is transformed into the value that belongs under the same dimension with all kinds of deviations and compares differentiating the size of whole synthesis assembling accuracy precision, and synthesis precision is high more, and then to be considered as the assembly quality effect good more.
Among the described embodiment, by using optimum target assembling pose matching process of the present invention, the utilization particle cluster algorithm carries out optimum target assembling pose matching algorithm and finds the solution, through after about 500 iteration, fitness function is restrained fully, matching result and digitizing master pattern name pose side-play amount are [0.009mm, 0.03mm, 0.0mm ,-0.002rad ,-0.004rad,-0.002rad], this shows that detect to require the actual optimum target assembling pose asked for the measured data of assembling parts very approaching with the nominal value of standard digital model according to assembling, this illustrates the validity of matching algorithm.After 1 realization is located to the cabin section according to this optimum target assembling pose, each accuracy error that assembles the accuracy matching result is respectively: parallel misalignment is 0.156mm between interface, the mean place degree deviation of each location hole is 0.184mm, two cabin section appearance coaxiality deviations are 0.129mm, adhere to specification the validity of this explanation matching algorithm.

Claims (8)

1. the method for best assembling pose coupling during a large parts assembles, it comprises that current pose coupling of parts and parts optimum target assembling pose mated for 2 stages:
Current pose matching stage comprises following basic step:
1) locus of fixed reference reference point in the mensuration erecting yard;
2) the reference point is fitted to benchmark assembling coordinate system, unify with global coordinate system in the digitizing assembling model;
3) locus of each position reference point on the mensuration removable assembly to be positioned;
4) nominal position of each position reference point in the contrast digitizing assembling model uses current pose matching algorithm to ask for six independent parameters of the current spatial pose of assembly;
Optimum target assembling pose matching stage comprises following basic step:
5) measure on the oriented benchmark assembly position with the corresponding measuring point of each assembling accuracy detection, fit to corresponding geometric element according to the type of reference data geometric element;
6) measure the position of assembling the corresponding measuring point of accuracy detection on the removable assembly to be positioned with each, and measuring point data is transformed under the local parts coordinate system of assembly;
7) each point position of connecting firmly with assembly to be positioned of supposition is along with assembly moves by rule in the certain space scope, the accuracy error of each detection when the calculating assembly is in the series of discrete position and attitude, and calculate the comprehensive overall accuracy Z-factor of multiple goal to the single goal mapping;
8) calculate each individual event accuracy error and overall accuracy Z-factor and whether satisfy the setting requirement, otherwise repeating step 7, until the spatial pose of the removable assembly to be positioned that obtains to meet the demands, this pose is the actual optimum target assembling pose of assembly to be positioned.
2. best assembling pose matching process according to claim 1 is characterized in that actual measurement data and digitizing standard assembling model are mated calculating, obtains the pose process of parts under global coordinate system:
The parts pose can be represented with the hexa-atomic group of v that rotation matrix and translation vector constitute The expression parts,, rotate around the z axle around y axle rotation θ degree successively around the x of reference frame axle rotation ψ degree from initial pose to current pose
Figure FDA0000049056220000012
Degree, at last at x, y, each translation p of z direction x, p y, p zObtain, spatial pose 6 independent parameters thus is that basic variable is expressed as 4 * 4 spatial alternation matrix T:
Figure FDA0000049056220000014
Figure FDA0000049056220000015
Wherein, R represents the rotation matrix component, and P represents translation vector, and with T presentation-entity spatial pose the time, the each point coordinate figure all is transformed into and participates in data operation behind the homogeneous form and handle, coordinate promptly by originally (x, y, z) TForm increases one dimension and becomes (x, y, z, 1) TFormal representation.
3. best assembling pose matching process according to claim 1, current pose matching stage, it is characterized in that the spatial pose of measurement parts to be positioned under global coordinate system all is expressed as the form of least square with reference to the coupling of the nominal position of eyeball and digitizing master pattern mid point and constraint, the optimization aim function as:
min{sum(ε i)},
ε wherein iThe matching error of expression each point, represent with the range deviation that physical location of putting and model are in the following nominal position of spatial pose state undetermined:
Figure FDA0000049056220000021
Constraint condition: E I_low<ε i<E I_up, E I_low, E I_upBe each reference point of regulation and the position deviation bound of theoretical value.
4. best assembling pose matching process according to claim 1, current pose matching stage is characterized in that eyeball and name point two groups of data computing processing procedures, wherein { P I-0Be nominal position point group on the entire assembly model to be measured, coordinate figure is with respect to local parts coordinate system, { P I-1Be eyeball group in position on the UUT, coordinate figure is with respect to assembling coordinate system, T 1For waiting to ask the parts spatial pose of representing with 6 independent parameters to be measured.On this basis, the optimization aim function further is refined as:
min{sum(ε i)}=min{sum(P i-1-T 1P i-0)}。
5. best assembling pose matching process according to claim 1, it is characterized in that using particle cluster algorithm to realize that current pose mates the calculating of mating with optimum target assembling pose and finds the solution, wherein each particle represents one of the parts spatial pose to separate, and it has 6 dimension x i=(x I1, x I2..., x ID) respectively representative show 6 independent variables of spatial pose
Figure FDA0000049056220000022
6. best assembling pose matching process according to claim 1, optimum target assembling pose matching stage is characterized in that the assembling accuracy detection constraint that relates to optimum target assembling pose comprises general dimensional tolerence, form and position tolerance and the special detection item type of formulating according to the product matching requirements.
7. best assembling pose matching process according to claim 1, optimum target assembling pose matching stage is characterized in that earlier each detection point coordinate value { P on the removable parts to be positioned kBe transformed into { P under the local coordinate system of parts from global coordinate system KL, calculate these measuring points then and in the space, move to global position coordinate figure { P under the different pose state T by rule with these parts KS, then in conjunction with benchmark measuring point group numerical value { P KDCalculate parts and be in that each may assemble the detection assessment result under the pose state in the space, its computing formula is:
ε I-k=f I-k(P kD,P kS),
f I-kFor calculating the function of detection numerical result, determine by type of detection.
8. best assembling pose matching process according to claim 1, optimum target assembling pose matching stage, it is characterized in that the matching constraint of positioning reference parts and the removable parts to be positioned assembling accuracy detection of each key characteristic under global coordinate system all is expressed as the form of detection deviate, adopt the Optimization Model of multiple goal, each detection coupling and each constraint are distributed weights W to the unified mapping of single goal kAfter ask the weighted comprehensive deviation, the optimization aim function as:
Figure FDA0000049056220000031
Figure FDA0000049056220000032
Figure FDA0000049056220000033
Constraint condition: E K_low<ε I-k<E K_up isEach matching requirements detection normal value bound,
ε wherein kThe matching error of representing k item assembling detection, T 1Result of calculation for current pose matching stage.
CN 201110053699 2011-01-19 2011-03-07 Method for matching optimal assembly poses of large-sized parts Pending CN102183205A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110053699 CN102183205A (en) 2011-01-19 2011-03-07 Method for matching optimal assembly poses of large-sized parts

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN201110021132.9 2011-01-19
CN201110021132 2011-01-19
CN 201110053699 CN102183205A (en) 2011-01-19 2011-03-07 Method for matching optimal assembly poses of large-sized parts

Publications (1)

Publication Number Publication Date
CN102183205A true CN102183205A (en) 2011-09-14

Family

ID=44569454

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110053699 Pending CN102183205A (en) 2011-01-19 2011-03-07 Method for matching optimal assembly poses of large-sized parts

Country Status (1)

Country Link
CN (1) CN102183205A (en)

Cited By (38)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102528398A (en) * 2012-02-14 2012-07-04 四川大学 Segmentation device for super-huge workpiece and machining method for segmentation device
CN103116668A (en) * 2013-01-25 2013-05-22 北京航空航天大学 Communication method of iGPS measurement system and computer-graphics aided three-dimensional interactive application (CATIA) software data
CN103164584A (en) * 2013-03-29 2013-06-19 江西洪都航空工业集团有限责任公司 Calculation method of coordination accuracy based on key characteristics
CN103246785A (en) * 2013-05-21 2013-08-14 西安电子科技大学 Force-directed virtual assembly method
CN104197835A (en) * 2014-09-11 2014-12-10 北京航空航天大学 Spatial position simulation and calibration method
CN104567681A (en) * 2015-01-08 2015-04-29 航天东方红卫星有限公司 Precise measurement method for satellite precise benchmark truss structure device
CN104625720A (en) * 2015-01-07 2015-05-20 北京卫星环境工程研究所 Automatic solar wing butt joint method based on laser tracker and adjusting parallel platform
CN104729455A (en) * 2014-12-04 2015-06-24 北京航空航天大学 Pipeline flexible welding and assembling robot posture calculation method based on measurement data
CN105157634A (en) * 2015-05-28 2015-12-16 无锡塔尔基热交换器科技有限公司 Metal product reference correction measurement method
CN105215669A (en) * 2015-09-17 2016-01-06 中国商用飞机有限责任公司 Attitude-adjusting method is measured in docking based on iGPS
CN105468834A (en) * 2015-11-20 2016-04-06 南昌航空大学 Assembly quality control method based on measurement point error ellipsoids
CN106200554A (en) * 2016-08-29 2016-12-07 广东工业大学 A kind of part processing pose optimization method and system
CN106247988A (en) * 2015-06-11 2016-12-21 北京卫星环境工程研究所 Celestial body attitude based on laser tracker and solar wing spreading frame accuracy measurement method
CN106897473A (en) * 2015-12-21 2017-06-27 北京航空航天大学 A kind of pose uncertainty evaluation method
CN107016215A (en) * 2017-04-27 2017-08-04 北京理工大学 A kind of assembling pose computational methods based on potential energy minimum
CN107045578A (en) * 2017-04-25 2017-08-15 江苏科技大学 A kind of chip pocket processing grinding wheel pose method for solving based on NPSO algorithms
CN105910535B (en) * 2016-05-10 2018-06-08 天津大学 A kind of single binocular pose bias measurement method for bay section automatic butt
CN108592852A (en) * 2018-05-18 2018-09-28 桂林电子科技大学 A kind of assessment method of the fast steady letter of the concentricity based on dual Tolerance Principle
CN108680129A (en) * 2018-05-18 2018-10-19 桂林电子科技大学 A kind of assessment method of the concentricity based on dual Tolerance Principle
CN108801135A (en) * 2016-05-31 2018-11-13 哈尔滨工业大学 Nuclear fuel rod pose automatic identification equipment
CN109191527A (en) * 2018-11-15 2019-01-11 凌云光技术集团有限责任公司 A kind of alignment method and device based on minimum range deviation
CN109766597A (en) * 2018-12-26 2019-05-17 北京理工大学 A kind of assembly pose high accuracy linearization method for solving considering geometric error
CN109773593A (en) * 2019-01-25 2019-05-21 重庆理工大学 A kind of method for grinding under the constraint condition based on surplus
CN109945824A (en) * 2017-12-21 2019-06-28 中国科学院长春光学精密机械与物理研究所 A kind of method and system of carrier aircraft heading reference measurement and transmitting
CN110296680A (en) * 2019-06-11 2019-10-01 清华大学 The attitude measurement method and system of robot crawl kinds of goods
CN110394686A (en) * 2019-06-26 2019-11-01 东华大学 It is a kind of actively to seek position and correction installation system
CN111145247A (en) * 2019-12-18 2020-05-12 配天机器人技术有限公司 Vision-based position detection method, robot and computer storage medium
CN111595238A (en) * 2020-06-03 2020-08-28 中国工程物理研究院机械制造工艺研究所 Laser tracker precision field evaluation system based on multi-station method
CN111907729A (en) * 2020-06-30 2020-11-10 成都飞机工业(集团)有限责任公司 Optimal pose fitting method for large part pose adjustment of airplane
CN112066877A (en) * 2020-08-28 2020-12-11 浙江西子势必锐航空工业有限公司 Method for assembling and adjusting airplane type frame
CN112197696A (en) * 2020-08-24 2021-01-08 北京卫星制造厂有限公司 Load interface hole making template adjusting and testing method based on integral optimization algorithm
CN112307561A (en) * 2020-10-30 2021-02-02 泉州装备制造研究所 Pose coordination method and system for assembling multi-body complex structure on large equipment
CN112348878A (en) * 2020-10-23 2021-02-09 歌尔光学科技有限公司 Positioning test method and device and electronic equipment
CN112504197A (en) * 2020-11-09 2021-03-16 航天材料及工艺研究所 Multi-station efficient assembly method based on unified standard
CN113405496A (en) * 2021-05-10 2021-09-17 中航西安飞机工业集团股份有限公司 Method for improving overall pose accuracy of component by optimizing and constructing measuring points
CN113405497A (en) * 2021-05-10 2021-09-17 中航西安飞机工业集团股份有限公司 Method for improving local pose accuracy of component by optimizing and constructing measuring points
CN113536459A (en) * 2021-07-15 2021-10-22 一汽奔腾轿车有限公司 Power assembly parameterization setting method in automobile digital design process
CN117984074A (en) * 2024-03-29 2024-05-07 成都飞机工业(集团)有限责任公司 Gesture-adjusting positioning method for large three-dimensional assembly without fixed measuring points

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20030062721A (en) * 2002-01-18 2003-07-28 한국항공우주산업 주식회사 Apparatus and method for fuselage automated splicing
WO2007125097A1 (en) * 2006-04-28 2007-11-08 Airbus Deutschland Gmbh Method and apparatus for ensuring the dimensional constancy of multisegment physical structures during assembly
CN101363715A (en) * 2008-09-26 2009-02-11 浙江大学 Computation method for attitude of aircraft fuselage based on laser tracking instrument
CN101456452B (en) * 2008-12-25 2010-10-13 浙江大学 Aircraft fuselage flexiblization and automation attitude-adjusting method
CN101362511B (en) * 2008-09-19 2010-11-10 浙江大学 Synergetic control method of aircraft part pose alignment based on four locater

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20030062721A (en) * 2002-01-18 2003-07-28 한국항공우주산업 주식회사 Apparatus and method for fuselage automated splicing
WO2007125097A1 (en) * 2006-04-28 2007-11-08 Airbus Deutschland Gmbh Method and apparatus for ensuring the dimensional constancy of multisegment physical structures during assembly
CN101362511B (en) * 2008-09-19 2010-11-10 浙江大学 Synergetic control method of aircraft part pose alignment based on four locater
CN101363715A (en) * 2008-09-26 2009-02-11 浙江大学 Computation method for attitude of aircraft fuselage based on laser tracking instrument
CN101456452B (en) * 2008-12-25 2010-10-13 浙江大学 Aircraft fuselage flexiblization and automation attitude-adjusting method

Cited By (57)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102528398B (en) * 2012-02-14 2014-02-19 四川大学 Segmentation device for super-huge workpiece and machining method for segmentation device
CN102528398A (en) * 2012-02-14 2012-07-04 四川大学 Segmentation device for super-huge workpiece and machining method for segmentation device
CN103116668A (en) * 2013-01-25 2013-05-22 北京航空航天大学 Communication method of iGPS measurement system and computer-graphics aided three-dimensional interactive application (CATIA) software data
CN103164584B (en) * 2013-03-29 2016-03-30 江西洪都航空工业集团有限责任公司 Based on the coordinating order of accuracy computing method of key characteristic
CN103164584A (en) * 2013-03-29 2013-06-19 江西洪都航空工业集团有限责任公司 Calculation method of coordination accuracy based on key characteristics
CN103246785A (en) * 2013-05-21 2013-08-14 西安电子科技大学 Force-directed virtual assembly method
CN104197835A (en) * 2014-09-11 2014-12-10 北京航空航天大学 Spatial position simulation and calibration method
CN104197835B (en) * 2014-09-11 2017-01-18 北京航空航天大学 Spatial position simulation and calibration method
CN104729455A (en) * 2014-12-04 2015-06-24 北京航空航天大学 Pipeline flexible welding and assembling robot posture calculation method based on measurement data
CN104625720B (en) * 2015-01-07 2017-01-18 北京卫星环境工程研究所 Automatic solar wing butt joint method based on laser tracker and adjusting parallel platform
CN104625720A (en) * 2015-01-07 2015-05-20 北京卫星环境工程研究所 Automatic solar wing butt joint method based on laser tracker and adjusting parallel platform
CN104567681A (en) * 2015-01-08 2015-04-29 航天东方红卫星有限公司 Precise measurement method for satellite precise benchmark truss structure device
CN105157634A (en) * 2015-05-28 2015-12-16 无锡塔尔基热交换器科技有限公司 Metal product reference correction measurement method
CN106247988A (en) * 2015-06-11 2016-12-21 北京卫星环境工程研究所 Celestial body attitude based on laser tracker and solar wing spreading frame accuracy measurement method
CN106247988B (en) * 2015-06-11 2018-07-20 北京卫星环境工程研究所 Celestial body posture based on laser tracker and solar wing spreading frame accuracy measurement method
CN105215669A (en) * 2015-09-17 2016-01-06 中国商用飞机有限责任公司 Attitude-adjusting method is measured in docking based on iGPS
CN105468834A (en) * 2015-11-20 2016-04-06 南昌航空大学 Assembly quality control method based on measurement point error ellipsoids
CN106897473A (en) * 2015-12-21 2017-06-27 北京航空航天大学 A kind of pose uncertainty evaluation method
CN105910535B (en) * 2016-05-10 2018-06-08 天津大学 A kind of single binocular pose bias measurement method for bay section automatic butt
CN108801135A (en) * 2016-05-31 2018-11-13 哈尔滨工业大学 Nuclear fuel rod pose automatic identification equipment
CN109000557A (en) * 2016-05-31 2018-12-14 哈尔滨工业大学 A kind of nuclear fuel rod pose automatic identifying method
CN109000557B (en) * 2016-05-31 2019-07-05 哈尔滨工业大学 A kind of nuclear fuel rod pose automatic identifying method
CN106200554B (en) * 2016-08-29 2018-09-04 广东工业大学 A kind of part processing pose optimization method and system
CN106200554A (en) * 2016-08-29 2016-12-07 广东工业大学 A kind of part processing pose optimization method and system
CN107045578A (en) * 2017-04-25 2017-08-15 江苏科技大学 A kind of chip pocket processing grinding wheel pose method for solving based on NPSO algorithms
CN107045578B (en) * 2017-04-25 2018-12-21 江苏科技大学 A kind of chip pocket processing grinding wheel pose method for solving based on NPSO algorithm
CN107016215A (en) * 2017-04-27 2017-08-04 北京理工大学 A kind of assembling pose computational methods based on potential energy minimum
CN107016215B (en) * 2017-04-27 2019-10-01 北京理工大学 A kind of assembly pose calculation method based on potential energy minimum
CN109945824B (en) * 2017-12-21 2022-03-22 中国科学院长春光学精密机械与物理研究所 Method and system for measuring and transmitting course reference of carrier
CN109945824A (en) * 2017-12-21 2019-06-28 中国科学院长春光学精密机械与物理研究所 A kind of method and system of carrier aircraft heading reference measurement and transmitting
CN108592852A (en) * 2018-05-18 2018-09-28 桂林电子科技大学 A kind of assessment method of the fast steady letter of the concentricity based on dual Tolerance Principle
CN108680129A (en) * 2018-05-18 2018-10-19 桂林电子科技大学 A kind of assessment method of the concentricity based on dual Tolerance Principle
CN109191527A (en) * 2018-11-15 2019-01-11 凌云光技术集团有限责任公司 A kind of alignment method and device based on minimum range deviation
CN109191527B (en) * 2018-11-15 2021-06-11 凌云光技术股份有限公司 Alignment method and device based on minimum distance deviation
CN109766597A (en) * 2018-12-26 2019-05-17 北京理工大学 A kind of assembly pose high accuracy linearization method for solving considering geometric error
CN109773593A (en) * 2019-01-25 2019-05-21 重庆理工大学 A kind of method for grinding under the constraint condition based on surplus
CN110296680B (en) * 2019-06-11 2020-11-24 清华大学 Method and system for measuring attitude of robot for grabbing goods
CN110296680A (en) * 2019-06-11 2019-10-01 清华大学 The attitude measurement method and system of robot crawl kinds of goods
CN110394686A (en) * 2019-06-26 2019-11-01 东华大学 It is a kind of actively to seek position and correction installation system
CN111145247B (en) * 2019-12-18 2023-07-07 配天机器人技术有限公司 Position degree detection method based on vision, robot and computer storage medium
CN111145247A (en) * 2019-12-18 2020-05-12 配天机器人技术有限公司 Vision-based position detection method, robot and computer storage medium
CN111595238A (en) * 2020-06-03 2020-08-28 中国工程物理研究院机械制造工艺研究所 Laser tracker precision field evaluation system based on multi-station method
CN111595238B (en) * 2020-06-03 2021-08-31 中国工程物理研究院机械制造工艺研究所 Laser tracker precision field evaluation system based on multi-station method
CN111907729A (en) * 2020-06-30 2020-11-10 成都飞机工业(集团)有限责任公司 Optimal pose fitting method for large part pose adjustment of airplane
CN112197696A (en) * 2020-08-24 2021-01-08 北京卫星制造厂有限公司 Load interface hole making template adjusting and testing method based on integral optimization algorithm
CN112066877B (en) * 2020-08-28 2022-02-18 浙江西子势必锐航空工业有限公司 Method for assembling and adjusting airplane type frame
CN112066877A (en) * 2020-08-28 2020-12-11 浙江西子势必锐航空工业有限公司 Method for assembling and adjusting airplane type frame
CN112348878A (en) * 2020-10-23 2021-02-09 歌尔光学科技有限公司 Positioning test method and device and electronic equipment
CN112348878B (en) * 2020-10-23 2023-03-21 歌尔科技有限公司 Positioning test method and device and electronic equipment
CN112307561A (en) * 2020-10-30 2021-02-02 泉州装备制造研究所 Pose coordination method and system for assembling multi-body complex structure on large equipment
CN112504197B (en) * 2020-11-09 2022-09-02 航天材料及工艺研究所 Multi-station efficient assembly method based on unified standard
CN112504197A (en) * 2020-11-09 2021-03-16 航天材料及工艺研究所 Multi-station efficient assembly method based on unified standard
CN113405496A (en) * 2021-05-10 2021-09-17 中航西安飞机工业集团股份有限公司 Method for improving overall pose accuracy of component by optimizing and constructing measuring points
CN113405497A (en) * 2021-05-10 2021-09-17 中航西安飞机工业集团股份有限公司 Method for improving local pose accuracy of component by optimizing and constructing measuring points
CN113405496B (en) * 2021-05-10 2022-07-12 中航西安飞机工业集团股份有限公司 Method for improving overall pose accuracy of component by optimizing and constructing measuring points
CN113536459A (en) * 2021-07-15 2021-10-22 一汽奔腾轿车有限公司 Power assembly parameterization setting method in automobile digital design process
CN117984074A (en) * 2024-03-29 2024-05-07 成都飞机工业(集团)有限责任公司 Gesture-adjusting positioning method for large three-dimensional assembly without fixed measuring points

Similar Documents

Publication Publication Date Title
CN102183205A (en) Method for matching optimal assembly poses of large-sized parts
CN111964575B (en) Digital twin modeling method for milling of mobile robot
CN111590581B (en) Positioning compensation method and device for robot
CN109848983A (en) A kind of method of highly conforming properties people guided robot work compound
CN109304730B (en) Robot kinematic parameter calibration method based on laser range finder
Sun et al. A review of robot control with visual servoing
CN109366503B (en) Large-scale component-oriented processing method based on mobile series-parallel robot
CN111238375B (en) Laser tracker-based appearance reconstruction method for large-scale component of mobile detection robot
CN109163675B (en) Method for detecting angular pendulum shaft position accuracy based on laser tracker
CN111426270B (en) Industrial robot pose measurement target device and joint position sensitive error calibration method
CN107457785B (en) Robot position compensation method based on joint feedback
CN114523475B (en) Automatic calibration and compensation device and method for errors of robot assembly system
CN109911241B (en) Seven-degree polynomial-based multi-section automatic posture adjusting method
CN109176487A (en) A kind of cooperating joint section scaling method, system, equipment, storage medium
Hu et al. Robot positioning error compensation method based on deep neural network
CN109591019A (en) A kind of Space Precision Orientation Method of no certainty location feature object
Peng et al. Development of an integrated laser sensors based measurement system for large-scale components automated assembly application
CN109883381A (en) A kind of three-dimensional space large-scale metrology method of articulated coordinate machine
Zheng et al. A novel algorithm of posture best fit based on key characteristics for large components assembly
CN112157654A (en) Optimization method for positioning error of robot machining system
Ding et al. Development of a novel mobile robotic system for large-scale manufacturing
Qi et al. Accuracy Improvement Calibrations for the Double‐Position 4‐PPPS Aircraft Docking System
Kong et al. Online measurement method for assembly pose of gear structure based on monocular vision
CN113932707B (en) Method for compensating geometric error of three-coordinate measuring machine based on elastic network algorithm
CN116245944A (en) Cabin automatic docking method and system based on measured data

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20110914