[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN109815562B - Assembly pose optimization method based on tolerance space - Google Patents

Assembly pose optimization method based on tolerance space Download PDF

Info

Publication number
CN109815562B
CN109815562B CN201910012090.9A CN201910012090A CN109815562B CN 109815562 B CN109815562 B CN 109815562B CN 201910012090 A CN201910012090 A CN 201910012090A CN 109815562 B CN109815562 B CN 109815562B
Authority
CN
China
Prior art keywords
assembly
pose
tolerance
space
uncertainty
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910012090.9A
Other languages
Chinese (zh)
Other versions
CN109815562A (en
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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201910012090.9A priority Critical patent/CN109815562B/en
Publication of CN109815562A publication Critical patent/CN109815562A/en
Application granted granted Critical
Publication of CN109815562B publication Critical patent/CN109815562B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/30Computing systems specially adapted for manufacturing

Landscapes

  • Automatic Assembly (AREA)
  • Numerical Control (AREA)

Abstract

Firstly, establishing an assembly constraint feature model, wherein the assembly constraint feature model is a mathematical model of assembly features formed by geometric features of an assembly object, and can calculate the numerical value of corresponding features under the condition of given assembly pose parameters; secondly, determining the uncertainty of assembly positioning of an assembly system; thirdly, an assembly pose evaluation model based on a tolerance space is established, wherein the tolerance space is a feasible domain space formed by the tolerance of assembly constraint features and the assembly pose constraint. For a given assembly pose parameter, evaluating the given assembly pose according to the evaluation basis by taking the distance from the pose to the space boundary of the feasible region; fourth, the optimal assembly pose is solved according to the assembly pose evaluation model.

Description

Assembly pose optimization method based on tolerance space
Technical Field
The invention relates to an assembly technology, in particular to an assembly pose optimization technology, and specifically relates to an assembly pose optimization method based on a tolerance space.
Background
In the digital assembly of the aircraft parts, the state of an assembly object is acquired through sensors such as laser and vision and is transmitted in the form of digital quantity in the process that the parts are positioned to the target pose through an assembly fixture. For parts that are not so rigid or structurally deformed, the pose of the assembled object is usually expressed by geometric kinematic parameters (rigid body translation amount and rotation amount). The actual state of the component deviates from the theoretical digital-analog due to component machining errors and component-to-component assembly cumulative errors. In the assembly of the aircraft parts, the deviation cannot be completely eliminated, the pose can be optimized only through an assembly coordination technology, and the deviation distribution is adjusted to meet the assembly process requirement.
Typically, the assembly object is located by multiple constraint features, so pose optimization of the assembly object is essentially a multi-objective constraint solving problem. Because of the difficulty in solving the multi-objective problem in engineering, the traditional method is to convert the multi-objective problem into a single-objective problem by means of weighted summation. The weight value is usually determined by experience or repeated experiments, and for the assembly of some complex structures, the weight value is often difficult to accurately set, and the unmanageable weight value can influence the optimization result of the assembly pose.
Disclosure of Invention
The invention aims at solving the problem that the weight value is difficult to accurately set in the assembly pose optimization and the assembly pose optimization result is affected, and provides an assembly pose optimization method based on a tolerance space.
The technical scheme of the invention is as follows:
the assembly pose optimization method based on the tolerance space is characterized by comprising the following steps of:
firstly, establishing an assembly constraint feature model;
secondly, determining the uncertainty of assembly positioning of an assembly system;
thirdly, an assembly pose evaluation model based on a tolerance space is established;
fourth, the optimal assembly pose is solved according to the assembly pose evaluation model.
The assembly constraint feature model is a mathematical model of assembly features formed by geometric features of an assembly object, and can calculate the numerical value of the corresponding features under the condition of given assembly pose parameters. In the product design stage, for product quality control, tolerance is defined for each assembly constraint feature, and the assembly constraint feature value after the assembly pose optimization is within the tolerance range.
The assembly positioning uncertainty of the assembly system is expressed by pose parameters of an assembly object.
The tolerance space is a feasible domain space formed by the constraint of tolerance of assembly constraint features on the pose of the assembly. And for the given assembly pose parameters, evaluating the given assembly pose according to the evaluation basis by taking the distance from the pose to the space boundary of the feasible region.
The invention has the beneficial effects that:
(1) The invention realizes the optimization of the assembly pose. The multi-constraint assembly pose solving is converted from the multi-objective optimizing problem to the single-objective optimizing problem, so that the solving difficulty is effectively reduced.
(2) The assembly pose evaluation model based on the tolerance space is used for evaluating the assembly pose quality without depending on experience or repeated experiments, so that the influence of incorrect weight value on the pose optimization result is avoided.
(3) The invention not only can acquire the optimized assembly pose, but also can judge whether the assembly result under the optimized pose can meet the assembly constraint condition according to the assembly pose evaluation parameters.
Drawings
FIG. 1 is a two-dimensional tolerance space diagram of the present invention.
FIG. 2 is a three-dimensional multi-constraint part tolerance space schematic of the present invention.
Detailed Description
The invention is further described below with reference to the drawings and examples.
As shown in fig. 1-2.
An assembly pose optimization method based on tolerance space comprises the following steps:
firstly, establishing an assembly constraint feature model;
secondly, determining the uncertainty of assembly positioning of an assembly system;
thirdly, an assembly pose evaluation model based on a tolerance space is established;
fourth, the optimal assembly pose is solved according to the assembly pose evaluation model.
Wherein: the assembly constraint feature model is a mathematical model of assembly features formed by geometric features of an assembly object, and can calculate the numerical value of the corresponding features under the condition of given assembly pose parameters. The specific modeling method of the assembly constraint characteristic model needs to be determined according to the specific geometric characteristics of the assembly object, and the established characteristic model is expressed as C ij (t i )∈[C ijmin ,C ijmax ]Wherein C ij Is the characteristic value, t i =[x y z α β γ] T Representing the assembly pose parameters, C ijmin ,C ijmax Upper and lower tolerance limits for assembly constraint features.
The positioning uncertainty space is represented by six-dimensional pose parameters as shown in a formula (6) when the positioning uncertainty δt of the assembly system is described by the position uncertainty δd and the pose uncertainty δθ, the uncertainty parameters can be determined according to the performance evaluation of the actual positioning system, and the change of the parameters can influence the assembly feasibility judgment but cannot influence the numerical value of the optimal pose.
Figure GDA0004066509290000031
Wherein δx δyδzδαδβδγ represents six-dimensional pose parameters of uncertainty space, and R represents real number domain.
The assembly pose evaluation model is established through the following processes:
first, a definition of the tolerance space is given, and as shown in fig. 1 (a) and 1 (b), a point P (x, y) ∈r is given on a two-dimensional plane 2 Is of the feasible region of (1)
Figure GDA0004066509290000032
For any one feasible solution P in the feasible domain i E.OMEGA, set P i To locate the theoretical target position, the actual locating position P is due to uncertainty in locating i ' deviations may exist. Definition: all P's meeting constraints under a certain positioning uncertainty i The' set of formations S is the tolerance space. Point-to-point P i In other words, if the uncertainty of positioning is consistent in all directions, the inscribed circle of the feasible region Ω is the point P i Is not limited by the tolerance space of (a). To ensure the integrity of the solution, for points P outside the feasible region i A similar definition is also given, as in fig. 1 (c), represented as an circumscribed circle of the feasible region Ω.
To quantify the tolerance space, the tolerance space is assessed in units of the set positioning uncertainty δt, defining this parameter as ease:
E=F(S,δt) (7)
the greater the ease E indicates the easier the positioning of the pose with the current solution as the target, the lower the positioning accuracy requirement. At point P i Taking the example as the case, the circle radius of the area corresponding to the positioning uncertainty δt is r i The circle radius of the corresponding area of the tolerance space S is R i (negative value of corresponding radius if S is corresponding to the circumscribed circle), the easiness of tolerance space correspondence is expressed as:
Figure GDA0004066509290000033
the visibility is a dimensionless constant for quantifying the tolerance space, from which definition the following law can be derived:
1) If E is more than or equal to 1, the current solution is in a feasible domain, and the positioning uncertainty requirement is met;
2) If E <1 is more than or equal to 0, the current solution is in a feasible region, and the actual position possibly exceeds the feasible region in consideration of positioning uncertainty;
3) If E <0, it indicates that the current solution is outside the feasible region, or that the feasible region does not exist.
Expanding the tolerance space assessment model into the part assembly pose assessment of the three-dimensional space, as shown in fig. 2, the feasible region omega is formed by various constraint features and the tolerance C thereof ij (t i )∈[C ijmin ,C ijmax ]Determining, wherein t i =[x y z α β γ] T Representing the fitting pose of the parts. To quantitatively describe the part positioning tolerance space, a scale variable ζ is introduced:
Δt i =(δt i ,ζ)=[ζ·δx ζ·δy ζ·δz δα δβ δγ] T (9)
then the part is in position t according to the definition of the tolerance space i Ease of (E) i Can be described as to any δt i E V, satisfy t i ′=t i +Δt i The maximum ζ still belonging to the tolerance space Ω:
Figure GDA0004066509290000041
if ζ i Ease of absence E i Described as the presence of δt i E V satisfies t i ′=t i +Δt i A negative value of the minimum ζ belonging to the tolerance space Ω.
The optimal assembly pose can be solved by adopting a particle swarm algorithm as a standard constraint optimization problem.
The invention is not related in part to the same as or can be practiced with the prior art.

Claims (3)

1. The assembly pose optimization method based on the tolerance space is characterized by comprising the following steps of:
firstly, establishing an assembly constraint feature model; the assembly constraint feature model is a mathematical model of assembly features formed by geometric features of an assembly object, and can calculate the numerical value of the corresponding features under the condition of given assembly pose parameters; in the product design stage, tolerance is defined for each assembly constraint feature for product quality control, and the assembly constraint feature value after the optimization of the assembly pose is within the tolerance range;
secondly, determining the uncertainty of assembly positioning of an assembly system;
thirdly, an assembly pose evaluation model based on a tolerance space is established;
fourthly, solving the optimal assembly pose according to the assembly pose evaluation model;
the assembly positioning uncertainty δt of the assembly system is described by the position uncertainty δd and the attitude uncertainty δθ, a positioning uncertainty space is represented by six-dimensional pose parameters as shown in a formula (6), the uncertainty can be determined according to the performance evaluation of an actual positioning system, and the change of the uncertainty can influence the assembly feasibility judgment but cannot influence the numerical value of the optimal pose;
Figure FDA0004066509280000011
wherein δx δy δz δα δβδγ respectively represent six-dimensional pose parameters of an uncertainty space, and R represents a real number domain;
the tolerance space is a feasible region space formed by the constraint of tolerance of assembly constraint features on the pose of the assembly; for a given assembly pose parameter, evaluating the given assembly pose according to the evaluation basis by taking the distance from the pose to the space boundary of the feasible region;
the method for establishing the assembly pose evaluation model comprises the following steps:
first, a definition of the tolerance space is given, and the point P (x, y) εR is given on the two-dimensional plane 2 Is of the feasible region of (1)
Figure FDA0004066509280000012
For any one feasible solution P in the feasible domain i E.OMEGA, set P i To locate the theoretical target position, the actual locating position P is due to uncertainty in locating i ' there may be a bias; definition: all P's meeting constraints under a certain positioning uncertainty i ' the set S of formations is a tolerance space; point-to-point P i In other words, if the uncertainty of positioning is consistent in all directions, the inscribed circle of the feasible region Ω is the point P i Is a tolerance space of (2); to ensure the integrity of the solution, for points P outside the feasible region i Similar definitions are also made;
to quantify the tolerance space, the tolerance space is assessed in units of the set positioning uncertainty δt, defining this parameter as ease:
E=F(S,δt) (2)
the greater the easiness E is, the easier the positioning of the pose with the current solution as the target pose is, and the lower the positioning precision requirement is; at point P i Taking the example as the case, the circle radius of the area corresponding to the positioning uncertainty δt is r i The circle radius of the corresponding area of the tolerance space S is R i If S is a negative value of the corresponding radius for the corresponding circumscribed circle, the easiness of the tolerance space correspondence is expressed as:
Figure FDA0004066509280000021
ease is a dimensionless constant used to quantify the tolerance space, from its definition, the following law can be derived:
1) If E is more than or equal to 1, the current solution is in a feasible domain, and the positioning uncertainty requirement is met;
2) If E <1 is more than or equal to 0, the current solution is in a feasible region, and the actual position possibly exceeds the feasible region in consideration of positioning uncertainty;
3) If E <0, the current solution is outside the feasible domain, or the feasible domain does not exist;
expanding the tolerance space evaluation model to part assembly pose evaluation in three-dimensional space, wherein the feasible domain omega is formed by various constraint features and tolerance C thereof ij (t i )∈[C ijmin ,C ijmax ]Determining, wherein t i =[x y z α β γ] T Representing the fitting pose of the parts; to quantitatively describe the part positioning tolerance space, a scale variable ζ is introduced:
Δt i =(δt i ,ζ)=[ζ·δx ζ·δy ζ·δz δα δβ δγ] T (4)
then the part is in position t according to the definition of the tolerance space i Ease of (E) i Can be described as to any δt i E V, satisfy t i ′=t i +Δt i The maximum ζ still belonging to the tolerance space Ω:
E i =F(S i ,δt i )=maxζ i
Figure FDA0004066509280000022
if ζ i Ease of absence E i Described as the presence of δt i E V satisfies t i ′=t i +Δt i A negative value of the minimum ζ belonging to the tolerance space Ω.
2. The method for optimizing assembly pose as claimed in claim 1, wherein said mathematical model of assembly features is C ij (t i )∈[C ijmin ,C ijmax ]Wherein C ij Is the characteristic value, t i =[x y z α β γ] T Representing the assembly pose parameters, C ijmin ,C ijmax Upper and lower tolerance limits for assembly constraint features.
3. The assembly pose optimization method according to claim 1, wherein the assembly positioning uncertainty of the assembly system is represented by pose parameters of the assembly object.
CN201910012090.9A 2019-01-07 2019-01-07 Assembly pose optimization method based on tolerance space Active CN109815562B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910012090.9A CN109815562B (en) 2019-01-07 2019-01-07 Assembly pose optimization method based on tolerance space

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910012090.9A CN109815562B (en) 2019-01-07 2019-01-07 Assembly pose optimization method based on tolerance space

Publications (2)

Publication Number Publication Date
CN109815562A CN109815562A (en) 2019-05-28
CN109815562B true CN109815562B (en) 2023-04-28

Family

ID=66603904

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910012090.9A Active CN109815562B (en) 2019-01-07 2019-01-07 Assembly pose optimization method based on tolerance space

Country Status (1)

Country Link
CN (1) CN109815562B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111044289B (en) * 2019-12-26 2021-09-03 哈尔滨工业大学 Large-scale high-speed rotation equipment alignment error measuring method based on closed-loop dynamic measurement
CN117648730B (en) * 2024-01-26 2024-04-19 北京适创科技有限公司 Method, device, equipment and storage medium for determining state of assembly

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106406277B (en) * 2016-09-23 2019-01-25 贵州珞石三盛科技有限公司 Robot kinematics' parameter error Optimization Compensation method and device

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于容差约束的机翼最优位姿评价算法;王青等;《机械工程学报》;20151031;第124-127页 *

Also Published As

Publication number Publication date
CN109815562A (en) 2019-05-28

Similar Documents

Publication Publication Date Title
CN111274671B (en) Precise repair method for complex product assembly process based on digital twin and operation system thereof
Jiang et al. A practical sampling method for profile measurement of complex blades
CN110695991B (en) Automatic hand-eye calibration method based on variation particle swarm optimization
CN104867136B (en) A kind of workpiece point cloud matching algorithm minimum based on distance variance
RU2729914C2 (en) Application of plant automation to aircraft assembly and method of its assembly
CN113741426B (en) Robot machining path planning method based on local point cloud curve fitting
CN102183205A (en) Method for matching optimal assembly poses of large-sized parts
CN109163675B (en) Method for detecting angular pendulum shaft position accuracy based on laser tracker
CN109597864A (en) Instant positioning and map constructing method and the system of ellipsoid boundary Kalman filtering
CN109815562B (en) Assembly pose optimization method based on tolerance space
CN106152935A (en) It is applicable to the laser measurement system without feature space curved surface blank and processing positioning method thereof
JP2018017826A (en) Autonomous moving body and environment map update device
Zhou et al. Off-line programming system of industrial robot for spraying manufacturing optimization
Wang et al. A pose estimation system based on deep neural network and ICP registration for robotic spray painting application
CN103106632A (en) Fusion method of different-accuracy three-dimension point cloud data based on mean shift
Wu et al. Adaptive location of repaired blade for multi-axis milling
CN112904358A (en) Laser positioning method based on geometric information
CN110900610A (en) Industrial robot calibration method based on LM algorithm and particle filter algorithm optimization
Khisamutdinov et al. Automation system goals for the creation and operation of the tool
Hartmann et al. Geometrical compensation of deterministic deviations for part finishing in bulk forming
Liu et al. Comparison of initial registration algorithms suitable for ICP algorithm
CN110412982B (en) Robot path robustness planning method based on monocular camera ranging uncertainty
Zhao et al. Profile and thickness constrained adaptive localization for manufacturing curved thin-walled parts based on on-machine measurement
CN108595373B (en) Uncontrolled DEM registration method
Berry et al. Closed-loop coordinate metrology for hybrid manufacturing system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant