CN102183205A - Method for matching optimal assembly poses of large-sized parts - Google Patents
Method for matching optimal assembly poses of large-sized parts Download PDFInfo
- 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
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
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
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
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:
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:
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
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:
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
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
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:
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:
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:
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
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:
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:
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
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:
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.
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)
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)
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 |
-
2011
- 2011-03-07 CN CN 201110053699 patent/CN102183205A/en active Pending
Patent Citations (5)
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)
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 |