CN111754567B - Comprehensive compensation method for static and dynamic errors in grinding and polishing processing of aircraft composite member robot - Google Patents
Comprehensive compensation method for static and dynamic errors in grinding and polishing processing of aircraft composite member robot Download PDFInfo
- Publication number
- CN111754567B CN111754567B CN202010504997.XA CN202010504997A CN111754567B CN 111754567 B CN111754567 B CN 111754567B CN 202010504997 A CN202010504997 A CN 202010504997A CN 111754567 B CN111754567 B CN 111754567B
- Authority
- CN
- China
- Prior art keywords
- robot
- coordinate system
- polishing
- processing
- grinding
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B24—GRINDING; POLISHING
- B24B—MACHINES, DEVICES, OR PROCESSES FOR GRINDING OR POLISHING; DRESSING OR CONDITIONING OF ABRADING SURFACES; FEEDING OF GRINDING, POLISHING, OR LAPPING AGENTS
- B24B1/00—Processes of grinding or polishing; Use of auxiliary equipment in connection with such processes
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J11/00—Manipulators not otherwise provided for
- B25J11/005—Manipulators for mechanical processing tasks
- B25J11/0065—Polishing or grinding
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B25—HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
- B25J—MANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
- B25J9/00—Programme-controlled manipulators
- B25J9/16—Programme controls
- B25J9/1656—Programme controls characterised by programming, planning systems for manipulators
- B25J9/1661—Programme controls characterised by programming, planning systems for manipulators characterised by task planning, object-oriented languages
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Robotics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Manipulator (AREA)
Abstract
本发明属于复材构件机器人磨抛加工领域,并具体公开了一种飞机复材构件机器人磨抛加工静动态误差综合补偿方法。该方法包括:搭建飞机复材构件机器人磨抛加工坐标系测量系统;建立机器人加工系统全局坐标系;构建机器人动力学递推模型;辨识机器人加工系统静态参数,建立机器人距离误差标定模型,以补偿机器人本体几何误差;构建机器人加工动态误差补偿模型;根据机器人加工动态误差理论数据集和先验知识数据集对机器人加工动态误差补偿模型进行优化训练;根据优化后的机器人加工动态误差补偿模型对机器人末端位姿进行误差补偿。本发明在动力学模型约束下精准构建坐标系全局控制网,便于误差链的全面分析,保证机器人磨抛加工型面精度。
The invention belongs to the field of robotic grinding and polishing processing of composite material components, and specifically discloses a comprehensive compensation method for static and dynamic errors in the robotic grinding and polishing processing of aircraft composite material components. The method includes: building a coordinate system measurement system for robot grinding and polishing of aircraft composite components; establishing a global coordinate system for the robot processing system; constructing a robot dynamics recursion model; identifying static parameters of the robot processing system, and establishing a robot distance error calibration model to compensate Robot body geometric error; construct a dynamic error compensation model for robot processing; optimize and train the dynamic error compensation model for robot processing based on the theoretical data set and prior knowledge data set for robot processing; train the robot based on the optimized dynamic error compensation model for robot processing Error compensation is performed on the terminal pose. This invention accurately constructs a coordinate system global control network under the constraints of the dynamic model, which facilitates comprehensive analysis of the error chain and ensures the accuracy of the robot grinding and polishing surface.
Description
技术领域Technical field
本发明属于复材构件机器人磨抛加工领域,更具体地,涉及一种飞机复材构件机器人磨抛加工静动态误差综合补偿方法。The invention belongs to the field of robotic grinding and polishing of composite components, and more specifically, relates to a comprehensive compensation method for static and dynamic errors in the robotic grinding and polishing of aircraft composite components.
背景技术Background technique
碳纤维复合材料因其耐高温、耐摩擦、高强度、密度小、耐腐蚀、轻质高强、性能可设计等优良性能,已广泛应用于航空航天结构材料,是实现航空航天装备轻量化和高效能化的关键,其用量已成为航空航天飞行器先进性和国际竞争力的标志之一。据统计,Z10直升机复合材料(以下简称“复材”)占比达90%,J20战机复材占比约40%,CR929客机复材占比约51%,C919客机复材占比约30%。Carbon fiber composite materials have been widely used in aerospace structural materials due to their excellent properties such as high temperature resistance, friction resistance, high strength, low density, corrosion resistance, lightweight and high strength, and designable performance. They are an important tool to achieve lightweight and high-performance aerospace equipment. The key to globalization, and its usage has become one of the symbols of the advancement and international competitiveness of aerospace vehicles. According to statistics, Z10 helicopter composite materials (hereinafter referred to as "composites") account for 90%, J20 fighter composite materials account for approximately 40%, CR929 passenger aircraft composite materials account for approximately 51%, and C919 passenger aircraft composite materials account for approximately 30% .
由于碳纤维复合材料易老化、屏敝率差,在飞机复材构件成型后其表面需打磨以便喷涂防护层来提升构件防腐蚀性、耐候性、耐磨性等防护性能。目前,已有的飞机复材构件表面打磨技术主要是人工打磨和人工助力打磨。需要工人手持或辅助打磨器械加工大型飞机复材构件,其存在加工精度低、产品一致性差、加工效率不高等诸多弊端,严重制约航空复材构件自动化、智能化生产的发展,成为航空制造技术领域的难点。Since carbon fiber composite materials are prone to aging and have poor barrier properties, after the aircraft composite components are formed, their surfaces need to be polished so that a protective layer can be sprayed on them to improve the component's anti-corrosion, weather resistance, wear resistance and other protective properties. At present, the existing surface grinding technologies for aircraft composite components are mainly manual grinding and manual-assisted grinding. Workers are required to process large-scale aircraft composite components with hand-held or auxiliary grinding instruments. This has many disadvantages such as low processing accuracy, poor product consistency, and low processing efficiency, which seriously restricts the development of automated and intelligent production of aviation composite components and has become a field of aviation manufacturing technology. difficulties.
相较于前两种飞机复材构件的表面打磨方式,机器人磨抛加工可发挥其柔性好、操作空间大、可扩展多种外部传感单元等优势,有望实现飞机复材构件的智能化加工。但目前机器人磨抛加工过程全局误差链建立不全面,难以有效捕捉姿态、标定、测量及加工等误差,严重影响工件的加工型面精度,限制了机器人磨抛加工在飞机复材构件智能化生产中的应用。Compared with the first two surface grinding methods of aircraft composite components, robot grinding and polishing can take advantage of its advantages such as good flexibility, large operating space, and the ability to expand a variety of external sensing units, and is expected to realize intelligent processing of aircraft composite components. . However, the current global error chain in the robot grinding and polishing process is not fully established, and it is difficult to effectively capture errors in attitude, calibration, measurement and processing, which seriously affects the processing surface accuracy of the workpiece and limits the use of robot grinding and polishing in the intelligent production of aircraft composite components. applications in.
基于上述缺陷和不足,本领域亟需提出一种飞机复材构件机器人磨抛加工静动态误差综合补偿方法,以解决实际飞机复材构件机器人磨抛加工中精度低、产品一致性差、加工效率不高等问题。Based on the above defects and deficiencies, there is an urgent need in this field to propose a comprehensive compensation method for static and dynamic errors in the robotic grinding and polishing of aircraft composite components to solve the problem of low precision, poor product consistency, and low processing efficiency in the actual robotic grinding and polishing of aircraft composite components. Advanced questions.
发明内容Contents of the invention
针对现有技术的以上缺陷或改进需求,本发明提供了一种飞机复材构件机器人磨抛加工静动态误差综合补偿方法,其中结合机器人磨抛加工静动态误差自身的特征及其机器人加工动态误差补偿模型优化的特点,相应的基于机器人加工系统距离误差标定模型补偿机器人本体几何参数误差,并在此基础上采用融合先验知识的新型支持向量机算法建立误差补偿模型,实时矫正机器人位姿及补偿末端加载变形误差,减少机器人加工动态误差,使得所加工的飞机复材构件精度高、产品一致性好以及加工效率高等特点。In view of the above defects or improvement needs of the existing technology, the present invention provides a comprehensive compensation method for static and dynamic errors in robot grinding and polishing of aircraft composite components, which combines the characteristics of static and dynamic errors in robot grinding and polishing with its dynamic errors in robot processing. According to the characteristics of compensation model optimization, the robot body geometric parameter error is compensated based on the distance error calibration model of the robot processing system. On this basis, a new support vector machine algorithm that integrates prior knowledge is used to establish an error compensation model to correct the robot pose and posture in real time. It compensates the end loading deformation error and reduces the dynamic error of robot processing, so that the processed aircraft composite components have high precision, good product consistency and high processing efficiency.
为实现上述目的,本发明提出了一种飞机复材构件机器人磨抛加工静动态误差综合补偿方法,包括以下步骤:In order to achieve the above objectives, the present invention proposes a comprehensive static and dynamic error compensation method for robot grinding and polishing of aircraft composite components, which includes the following steps:
S1搭建飞机复材构件机器人磨抛加工坐标系测量系统;S1 builds a coordinate system measurement system for robot grinding and polishing of aircraft composite components;
S2建立机器人加工系统全局坐标系;S2 establishes the global coordinate system of the robot processing system;
S3基于所述机器人加工系统全局坐标系构建机器人动力学递推模型;S3 builds a robot dynamics recursive model based on the global coordinate system of the robot processing system;
S4根据所述机器人动力学递推模型辨识机器人加工系统静态参数,并根据该机器人加工系统静态参数建立机器人距离误差标定模型,以补偿机器人本体几何误差;S4 identifies the static parameters of the robot processing system based on the robot dynamics recursion model, and establishes a robot distance error calibration model based on the static parameters of the robot processing system to compensate for the robot body geometric error;
S5根据机器人本体几何误差构建机器人加工系统全局坐标系下机器人加工动态误差补偿模型;S5 builds a dynamic error compensation model for robot processing in the global coordinate system of the robot processing system based on the geometric error of the robot body;
S6根据机器人加工动态误差理论数据集和先验知识数据集对机器人加工动态误差补偿模型进行优化训练;S6 optimizes and trains the robot processing dynamic error compensation model based on the robot processing dynamic error theory data set and the prior knowledge data set;
S7根据优化后的机器人加工动态误差补偿模型对机器人末端位姿进行误差补偿。S7 performs error compensation on the robot end pose based on the optimized robot processing dynamic error compensation model.
作为进一步优选的,步骤S1中,所述机器人磨抛加工坐标系测量系统包括激光跟踪仪和双目测量设备。As a further preference, in step S1, the robot grinding and polishing processing coordinate system measurement system includes a laser tracker and binocular measurement equipment.
作为进一步优选的,步骤S2中,所述机器人加工系统全局坐标系包括工件坐标系、工具坐标系、测量坐标系以及机器人全局基准坐标系。As a further preference, in step S2, the global coordinate system of the robot processing system includes a workpiece coordinate system, a tool coordinate system, a measurement coordinate system and a robot global reference coordinate system.
作为进一步优选的,步骤S2中,所述机器人加工系统全局坐标系采用基于离散型差分进化的非线性优化算法建立。As a further preference, in step S2, the global coordinate system of the robot processing system is established using a nonlinear optimization algorithm based on discrete differential evolution.
作为进一步优选的,步骤S5中,采用融合先验知识的新型支持向量机回归算法构建机器人加工动态误差补偿模型。As a further preference, in step S5, a new support vector machine regression algorithm that integrates prior knowledge is used to construct a dynamic error compensation model for robot processing.
作为进一步优选的,步骤S6中,采用混沌粒子群优化算法对机器人加工动态误差补偿模型进行优化训练。As a further preference, in step S6, a chaotic particle swarm optimization algorithm is used to optimize and train the robot processing dynamic error compensation model.
作为进一步优选的,步骤S6具体包括以下步骤:As further preferred, step S6 specifically includes the following steps:
S61选取一组机器人磨抛加工空间内分散的目标点集合A,并通过所述机器人磨抛加工坐标系测量系统测量得到所述目标点集合A在测量坐标系W1下的位姿PAM,获取所述目标点集合A在机器人全局基准坐标系W2下的理论位姿PA,然后求取所述测量坐标系W1与所述机器人全局基准坐标系W2间的坐标系转换矩阵T;S61 selects a set of target points A scattered in the robot grinding and polishing processing space, and measures the pose P AM of the target point set A in the measurement coordinate system W 1 through the robot grinding and polishing processing coordinate system measurement system. Obtain the theoretical pose P A of the target point set A in the robot's global reference coordinate system W 2 , and then obtain the coordinate system transformation matrix T between the measurement coordinate system W 1 and the robot's global reference coordinate system W 2 ;
S62确定一组机器人末端待补偿目标点集合B,并将该待补偿目标点集合B按预设比例分成机器人加工动态误差先验知识数据集Bp以及剩余数据集CBBp;S62 determines a set B of target points to be compensated at the end of the robot, and divides the set B of target points to be compensated into a robot processing dynamic error prior knowledge data set B p and a remaining data set C B B p according to a preset ratio;
S63在所述测量坐标系W1下测量得到所述先验知识数据集Bp中各点的目标位姿PBpm,然后通过所述坐标系转换矩阵T得到所述先验知识数据集Bp中各点在所述机器人全局基准坐标系W2下的实际位姿PBpr;S63 Measure the target pose P Bpm of each point in the prior knowledge data set B p under the measurement coordinate system W 1 , and then obtain the prior knowledge data set B p through the coordinate system transformation matrix T. The actual pose P Bpr of each point in the robot's global reference coordinate system W 2 ;
S64根据所述先验知识数据集Bp中各点的实际位姿PBpr和理论位姿PBpt得到所述先验知识数据集Bp中各点对应的矫正位姿PBpr’;S64 obtains the corrected pose P Bpr ' corresponding to each point in the prior knowledge data set B p according to the actual pose P Bpr and the theoretical pose P Bpt of each point in the prior knowledge data set B p ;
S65将所述先验知识数据集Bp中各点的实际位姿PBpr作为机器人加工动态误差补偿模型的输入,所述矫正位姿PBpr’作为机器人加工动态误差补偿模型的输出,对机器人加工动态误差补偿模型进行优化训练;S65 uses the actual pose P Bpr of each point in the prior knowledge data set B p as the input of the robot processing dynamic error compensation model, and the corrected pose P Bpr ' as the output of the robot processing dynamic error compensation model. The machining dynamic error compensation model is optimized and trained;
S66将所述剩余数据集CBBp的目标位姿作为优化训练后的机器人加工动态误差补偿模型的输入,并由该模型输出矫正后目标点集合B’,并将该矫正后目标点集合B’作为机器人磨抛加工静动态误差补偿的依据。S66 uses the target pose of the remaining data set C B B p as the input of the optimized and trained robot processing dynamic error compensation model, and the model outputs the corrected target point set B', and the corrected target point set B' serves as the basis for static and dynamic error compensation in robot grinding and polishing.
作为进一步优选的,步骤S61中,目标点集合A的获取方法如下:As a further preference, in step S61, the target point set A is obtained as follows:
S611设定构成集合的样本点的个数;S611 sets the number of sample points that constitute the set;
S612对机器人磨抛加工空间进行区域划分;S612 divides the robot grinding and polishing processing space into regions;
S613根据机器人磨抛加工路径以及构成集合的样本点的个数,在机器人磨抛加工空间划分后的各区间随机选取样本点,从而得到目标点集合A。S613 randomly selects sample points in each interval after the robot grinding and polishing processing space is divided according to the robot grinding and polishing processing path and the number of sample points that constitute the set, thereby obtaining the target point set A.
作为进一步优选的,步骤S7中,对机器人末端位姿进行误差补偿包括对机器人末端位置以及姿态的补偿。As a further preference, in step S7, error compensation for the robot end pose includes compensation for the robot end position and attitude.
作为进一步优选的,机器人末端位置的补偿通过位置补偿器实现,机器人末端姿态的补偿通过姿态补偿器实现。As a further preference, the robot end position is compensated by a position compensator, and the robot end attitude is compensated by an attitude compensator.
总体而言,通过本发明所构思的以上技术方案与现有技术相比,主要具备以下的技术优点:Generally speaking, compared with the existing technology, the above technical solution conceived by the present invention mainly has the following technical advantages:
1.本发明在所述全局坐标系基础上建立动力学递推模型,辨识机器人加工系统的静态参数,依据所述静态参数得到机器人加工系统距离误差标定模型,补偿机器人本体几何参数误差;最终采用融合先验知识的新型支持向量机算法建立误差补偿模型,实时矫正机器人位姿及补偿末端加载变形误差,减少机器人加工动态误差。本发明在动力学模型约束下精准构建坐标系全局控制网,便于误差链的全面分析,保证机器人磨抛加工型面精度。1. The present invention establishes a dynamic recursion model based on the global coordinate system, identifies the static parameters of the robot processing system, obtains the distance error calibration model of the robot processing system based on the static parameters, and compensates the robot body geometric parameter error; finally adopts The new support vector machine algorithm that integrates prior knowledge establishes an error compensation model, corrects the robot's posture and compensates for end loading deformation errors in real time, and reduces dynamic errors in robot processing. This invention accurately constructs a coordinate system global control network under the constraints of the dynamic model, which facilitates comprehensive analysis of the error chain and ensures the accuracy of the robot grinding and polishing surface.
2.本发明建立的距离误差标定模型对机器人加工系统动态补偿的初始值进行了优化,显著的提高了迭代效率和迭代精度。2. The distance error calibration model established by the present invention optimizes the initial value of dynamic compensation of the robot processing system, significantly improving iteration efficiency and iteration accuracy.
3.本发明采用混沌粒子群优化算法具有结构简单、学习速度快和泛化性能突出等特点,实现了机器人位姿和加工轨迹的实时动态补偿,提高了机器人磨抛加工系统的准确性和鲁棒性。3. The present invention adopts the chaotic particle swarm optimization algorithm, which has the characteristics of simple structure, fast learning speed and outstanding generalization performance. It realizes real-time dynamic compensation of robot posture and processing trajectory, and improves the accuracy and robustness of the robot grinding and polishing processing system. Great sex.
4.本发明对机器人磨抛加工空间进行区域划分,并根据机器人磨抛加工路径以及构成集合的样本点的个数,提高了样本空间的遍历性,使得优化后的机器人加工动态误差补偿模型适应性更强,获取的机器人末端位姿补偿更加精确。4. The present invention divides the robot grinding and polishing processing space into regions, and improves the traversability of the sample space according to the robot grinding and polishing processing path and the number of sample points constituting the set, so that the optimized robot processing dynamic error compensation model can adapt to The accuracy is stronger, and the obtained robot end pose compensation is more accurate.
附图说明Description of drawings
图1为本发明实施例涉及的一种飞机复材构件机器人磨抛加工静动态误差综合补偿方法的流程示意图;Figure 1 is a schematic flow chart of a comprehensive compensation method for static and dynamic errors in robotic grinding and polishing of aircraft composite components according to an embodiment of the present invention;
图2是本发明实施例中涉及的机器人位姿及末端加载变形误差补偿流程图。Figure 2 is a flow chart of robot posture and terminal loading deformation error compensation involved in the embodiment of the present invention.
具体实施方式Detailed ways
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。In order to make the purpose, technical solutions and advantages of the present invention more clear, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described here are only used to explain the present invention and are not intended to limit the present invention. In addition, the technical features involved in the various embodiments of the present invention described below can be combined with each other as long as they do not conflict with each other.
如图1和图2所示,本发明实施例提供的一种飞机复材构件机器人磨抛加工静动态误差综合补偿方法,采用融合先验知识的新型支持向量机算法实现机器人末端目标位姿与修正位姿的映射,使得补偿后的实际位姿与目标位姿保持一致,以提高机器人磨抛加工的精度。在本发明中,先验知识指机器人末端目标点的理论位姿。具体实现方法如下:As shown in Figures 1 and 2, an embodiment of the present invention provides a comprehensive compensation method for static and dynamic errors in robot grinding and polishing of aircraft composite components, using a new support vector machine algorithm that integrates prior knowledge to achieve the robot end target pose and posture. Correct the mapping of pose so that the actual pose after compensation is consistent with the target pose to improve the accuracy of robot grinding and polishing. In the present invention, prior knowledge refers to the theoretical pose of the robot's end target point. The specific implementation method is as follows:
步骤一,搭建飞机复材构件机器人磨抛加工坐标系测量系统。Step 1: Build a coordinate system measurement system for robot grinding and polishing of aircraft composite components.
本发明中,坐标系测量系统包括激光跟踪仪和双目测量设备。或者,为了进一步提高测量的精度、速度和分辨率,坐标系测量系统为由多个相机、多个标定球以及标定板构成的Optitrack系统。In the present invention, the coordinate system measurement system includes a laser tracker and binocular measurement equipment. Or, in order to further improve the accuracy, speed and resolution of measurement, the coordinate system measurement system is an Optitrack system composed of multiple cameras, multiple calibration balls and calibration plates.
步骤二,基于离散型差分进化的非线性优化算法建立机器人加工系统全局坐标系。Step 2: Establish the global coordinate system of the robot processing system based on the nonlinear optimization algorithm of discrete differential evolution.
根据步骤一搭建的坐标系测量系统,采用离散型差分进化的非线性优化算法建立机器人加工系统全局坐标系。其中,机器人加工系统全局坐标系至少包括工件坐标系、工具坐标系、测量坐标系以及机器人全局基准坐标系。工件坐标系为待加工的飞机复杂构件所在的坐标系,工具坐标系为机器人末端夹持的工具所在的坐标系,测量坐标系为测量系统所在的坐标系,机器人全局基准坐标系为机器人的基坐标系。通过相对位置关系以及机器人的运动学关系,获取各个坐标系之间的转化关系,进而可测量获取机器人末端位姿信息,即为机器人末端的实际位姿。在本发明中,首先构建离散型差分进化的非线性优化算法模型,然后将坐标系测量系统测量得到的工件、工具、机器人以及参照物之间的位姿关系等作为可变参数输入离散型差分进化的非线性优化算法模型,得到工件坐标系、工具坐标系、测量坐标系以及机器人全局基准坐标系及各坐标系之间的转换关系。坐标测量系统测量得到的机器人末端的实际位姿与控制系统输入的机器人末端理论位姿之间的差值即为机器人末端位姿误差。其中,位姿误差包括位置误差和姿态误差,位置误差包括机器人末端沿X轴、Y轴以及Z轴的误差,姿态误差包括机器人末端α角、β角以及γ角的角度误差。According to the coordinate system measurement system built in step 1, the nonlinear optimization algorithm of discrete differential evolution is used to establish the global coordinate system of the robot processing system. Among them, the global coordinate system of the robot processing system at least includes the workpiece coordinate system, the tool coordinate system, the measurement coordinate system and the robot global reference coordinate system. The workpiece coordinate system is the coordinate system where the complex aircraft components to be processed are located, the tool coordinate system is the coordinate system where the tool held by the end of the robot is located, the measurement coordinate system is the coordinate system where the measurement system is located, and the robot global reference coordinate system is the base of the robot. Coordinate System. Through the relative position relationship and the kinematic relationship of the robot, the transformation relationship between each coordinate system is obtained, and then the robot end pose information can be measured and obtained, which is the actual pose of the robot end. In the present invention, a nonlinear optimization algorithm model of discrete differential evolution is first constructed, and then the pose relationships between workpieces, tools, robots and reference objects measured by the coordinate system measurement system are input into the discrete differential algorithm as variable parameters. The evolved nonlinear optimization algorithm model obtains the workpiece coordinate system, tool coordinate system, measurement coordinate system, robot global reference coordinate system and the conversion relationship between each coordinate system. The difference between the actual pose of the robot end measured by the coordinate measurement system and the theoretical pose of the robot end input by the control system is the robot end pose error. Among them, the pose error includes position error and attitude error. The position error includes the error of the robot end along the X-axis, Y-axis and Z-axis. The attitude error includes the angle error of the α angle, β angle and γ angle of the robot end.
步骤三,基于所述机器人加工系统全局坐标系构建机器人动力学递推模型,以获取机器人关节力矩关系式。该机器人关节力矩关系式包含了机器人关节力矩、角度、角速度以及角的加速度。Step 3: Construct a robot dynamics recursion model based on the global coordinate system of the robot processing system to obtain the robot joint torque relationship. The robot joint torque relationship includes the robot joint torque, angle, angular velocity and angular acceleration.
步骤四,根据所述机器人动力学递推模型辨识机器人加工系统静态参数,并根据该机器人加工系统静态参数建立机器人距离误差标定模型,以补偿机器人本体几何误差,并根据该距离误差标定模型对机器人本体几何误差进行补偿。一般而言,可采用理论辨识法辨识机器人加工系统静态参数。其中机器人本体几何误差主要由机器人安装误差、工具安装误差等因几何构建安装而产生的误差。通过基于所述机器人加工系统全局坐标系构建机器人动力学递推模型,并以此来辨别机器人加工系统静态参数,可实现对机器人静态误差的控制和补偿。Step 4: Identify the static parameters of the robot processing system based on the robot dynamics recursive model, and establish a robot distance error calibration model based on the static parameters of the robot processing system to compensate for the geometric error of the robot body, and calibrate the robot based on the distance error calibration model. Body geometry errors are compensated. Generally speaking, the theoretical identification method can be used to identify the static parameters of the robot processing system. Among them, the geometric error of the robot body is mainly caused by the robot installation error, tool installation error and other errors caused by geometric construction and installation. By constructing a robot dynamics recursion model based on the global coordinate system of the robot processing system, and using this to identify the static parameters of the robot processing system, the control and compensation of the static errors of the robot can be achieved.
步骤五,根据机器人本体几何误差构建机器人加工系统全局坐标系下机器人加工动态误差补偿模型。其中,机器人加工动态误差补偿模型采用融合先验知识的新型支持向量机回归算法构建。即将先验知识与新型支持向量机回归算法融合,构建机器人加工系统全局坐标系下机器人加工动态误差补偿模型。在本发明中,先验知识指机器人末端目标点的理论位姿,一组机器人末端目标点的理论位姿构成先验知识数据集。在本发明中,首先基于先验知识的新型支持向量机回归算法构建机器人加工动态误差补偿模型,并将机器人加工动态参数、机器人末端加工变形量、传动误差、控制误差、阈值参数、权值系数等作为基于先验知识的新型支持向量机回归算法的参数,构建机器人加工动态误差补偿模型。Step 5: Construct a dynamic error compensation model for robot processing in the global coordinate system of the robot processing system based on the geometric error of the robot body. Among them, the robot machining dynamic error compensation model is constructed using a new support vector machine regression algorithm that integrates prior knowledge. That is, prior knowledge is integrated with a new support vector machine regression algorithm to construct a dynamic error compensation model for robot processing in the global coordinate system of the robot processing system. In the present invention, prior knowledge refers to the theoretical pose of the robot's end target point, and a group of theoretical poses of the robot's end target point constitutes a priori knowledge data set. In the present invention, a new support vector machine regression algorithm based on prior knowledge first constructs a dynamic error compensation model for robot processing, and combines the robot processing dynamic parameters, robot end processing deformation, transmission error, control error, threshold parameters, and weight coefficients. etc. are used as parameters of a new support vector machine regression algorithm based on prior knowledge to construct a dynamic error compensation model for robot machining.
步骤六,根据机器人加工动态误差理论数据集和先验知识数据集对机器人加工动态误差补偿模型进行优化训练。整个迭代优化训练的过程采用混沌粒子群优化算法进行。其具体包括以下步骤:Step 6: Optimize and train the robot processing dynamic error compensation model based on the robot processing dynamic error theory data set and the prior knowledge data set. The entire iterative optimization training process is carried out using the chaotic particle swarm optimization algorithm. It specifically includes the following steps:
(1)选取一组机器人磨抛加工空间内分散的目标点集合A,并通过所述机器人磨抛加工坐标系测量系统测量得到所述目标点集合A在测量坐标系W1下的位姿PAM,获取所述目标点集合A在机器人全局基准坐标系W2下的理论位姿PA,然后求取所述测量坐标系W1与所述机器人全局基准坐标系W2间的坐标系转换矩阵T;(1) Select a set of target points A scattered in the robot grinding and polishing processing space, and measure the pose P of the target point set A in the measurement coordinate system W 1 through the robot grinding and polishing processing coordinate system measurement system. AM , obtain the theoretical pose P A of the target point set A in the robot's global reference coordinate system W 2 , and then obtain the coordinate system transformation between the measurement coordinate system W 1 and the robot's global reference coordinate system W 2 matrixT;
(2)确定一组机器人末端待补偿目标点集合B,并将该待补偿目标点集合B按预设比例分成机器人加工动态误差先验知识数据集Bp以及剩余数据集CBBp;(2) Determine a set B of target points to be compensated at the end of the robot, and divide the set B of target points to be compensated into the robot machining dynamic error prior knowledge data set B p and the remaining data set C B B p according to a preset ratio;
(3)在所述测量坐标系W1下测量得到所述先验知识数据集Bp中各点的目标位姿PBpm,然后通过所述坐标系转换矩阵T得到所述先验知识数据集Bp中各点在所述机器人全局基准坐标系W2下的实际位姿PBpr;(3) Measure the target pose P Bpm of each point in the prior knowledge data set B p under the measurement coordinate system W 1 , and then obtain the prior knowledge data set through the coordinate system transformation matrix T The actual pose P Bpr of each point in B p under the robot's global reference coordinate system W 2 ;
(4)根据所述先验知识数据集Bp中各点的实际位姿PBpr和理论位姿PBpt得到所述先验知识数据集Bp中各点对应的矫正位姿PBpr’;(4) According to the actual pose P Bpr and theoretical pose P Bpt of each point in the prior knowledge data set B p , obtain the corrected pose P Bpr ' corresponding to each point in the prior knowledge data set B p ;
(5)将所述先验知识数据集Bp中各点的实际位姿PBpr作为机器人加工动态误差补偿模型的输入,所述矫正位姿PBpr’作为机器人加工动态误差补偿模型的输出,对机器人加工动态误差补偿模型进行优化训练;(5) Use the actual pose P Bpr of each point in the prior knowledge data set B p as the input of the robot processing dynamic error compensation model, and the corrected pose P Bpr ' as the output of the robot processing dynamic error compensation model, Optimize and train the dynamic error compensation model for robot processing;
(6)将所述剩余数据集CBBp的目标位姿作为优化训练后的机器人加工动态误差补偿模型的输入,并由该模型输出矫正后目标点集合B’,并将该矫正后目标点集合B’作为机器人磨抛加工静动态误差补偿的依据。(6) Use the target pose of the remaining data set C B B p as the input of the optimized and trained robot processing dynamic error compensation model, and output the corrected target point set B' from the model, and use the corrected target Point set B' serves as the basis for static and dynamic error compensation in robot grinding and polishing processing.
作为本发明的优选方案,步骤(6)中,将所述剩余数据集CBBp中各点的目标位姿作为步骤(5)优化训练后的机器人加工动态误差补偿模型的输入,训练初步优化训练后的机器人加工动态误差补偿模型输出的补偿值向矫正位姿PCBpr’转化的能力,使得补偿值与矫正位姿PCBpr’的差值在阈值范围内,完成机器人加工动态误差补偿模型的优化。其中,矫正位姿PCBpr’为剩余数据集CBBp中各点的实际位姿PCBpm和理论位姿PCBpt得到所述剩余数据集CBBp中各点对应的矫正位姿PCBpr’。实际位姿PCBpr为在所述测量坐标系W1下测量得到所述剩余数据集CBBp中各点的位姿PCBpm,然后通过所述坐标系转换矩阵T得到所述剩余数据集CBBp中各点在所述机器人全局基准坐标系W2下的实际位姿PCBpr。As a preferred solution of the present invention, in step (6), the target pose of each point in the remaining data set C B B p is used as the input of the robot processing dynamic error compensation model optimized and trained in step (5), and the preliminary training Optimize the ability of the trained robot processing dynamic error compensation model to transform the compensation value output into the corrected pose P CBpr ', so that the difference between the compensation value and the corrected pose P CBpr ' is within the threshold range, completing the robot processing dynamic error compensation model. Optimization. Among them, the corrected pose P CBpr ' is the actual pose P CBpm and the theoretical pose P CBpt of each point in the remaining data set C B B p to obtain the corrected pose P corresponding to each point in the remaining data set C B B p CBpr '. The actual pose P CBpr is the pose P CBpm of each point in the remaining data set C B B p measured in the measurement coordinate system W 1 , and then the remaining data set is obtained through the coordinate system transformation matrix T The actual pose P CBpr of each point in C B B p under the robot's global reference coordinate system W 2 .
更具体的,在本发明中,为了提高训练样本的适用性,需要设定构成目标点集合A的样本个数,原则上样本个数越多,结果越精确。然后对机器人磨抛加工空间进行区域划分。根据机器人磨抛加工路径以及构成集合的样本点的个数,在机器人磨抛加工空间划分后的各区间随机选取样本点,从而得到目标点集合A。优选的,根据机器人磨抛加工路径覆盖量,选取样本点的个数,如即,覆盖量越大,选取的样本点越多。More specifically, in the present invention, in order to improve the applicability of training samples, it is necessary to set the number of samples that constitute the target point set A. In principle, the more the number of samples, the more accurate the result. Then the robot grinding and polishing processing space is divided into regions. According to the robot grinding and polishing processing path and the number of sample points that constitute the set, sample points are randomly selected from each interval after the robot grinding and polishing processing space is divided, thereby obtaining the target point set A. Preferably, the number of sample points is selected according to the coverage amount of the robot grinding and polishing processing path. For example, the greater the coverage amount, the more sample points are selected.
步骤七,根据优化后的机器人加工动态误差补偿模型对机器人末端位姿进行误差补偿。即使用所述机器人加工动态误差补偿模型矫正机器人位姿及补偿末端加载变形误差。Step 7: Perform error compensation on the end position of the robot based on the optimized robot processing dynamic error compensation model. That is, the robot machining dynamic error compensation model is used to correct the robot posture and compensate for the terminal loading deformation error.
图2所示为本发明实施例中的机器人位姿及末端加载变形误差补偿流程图,通过所述机器人磨抛加工坐标系测量系统测量得到一组测量坐标系W1下目标点集合A的机器人末端位姿PAM,同时通过机器人软件端获得与之对应的一组机器人全局基准坐标系W2下的理论位姿PA,利用LM算法求得所述测量坐标系W1与所述机器人全局基准坐标系W2间的坐标系转换矩阵T。其次,对于待补偿目标点集合B,从中随机选取小部分目标点构成先验知识数据集Bp,测得其在所述测量坐标系W1下的位姿PBpm,利用所述坐标系转换矩阵T求得所述先验知识数据集Bp在所述机器人全局基准坐标系W2下的实际位姿PBpr,将所述实际位姿PBpr与其相对应的理论位姿PBpt相比较得到所述先验知识数据集Bp中对应目标点的矫正位姿PBpr’,以所述先验知识数据集Bpre中目标点的实际坐标PBpr作为输入,以与其相对应的矫正位姿PBpr’作为输出,训练融合先验知识的新型支持向量机回归算法建立的机器人加工动态误差补偿模型。最终,将所述待补偿目标点集合B的剩余目标点所构成的集合,即剩余数据集CBBp作为输入,经由融合先验知识的新型支持向量机回归算法建立的模型输出得到矫正后的目标点集合B’,将所述矫正后的目标点集合B’中点的位姿输入机器人控制器,补偿机器人对末端误差。Figure 2 shows a flow chart of robot posture and end loading deformation error compensation in the embodiment of the present invention. Through the measurement of the robot grinding and polishing processing coordinate system measurement system, a set of robots with target point set A under the measurement coordinate system W1 is obtained. At the same time, the end pose P AM is obtained through the robot software end, and the corresponding theoretical pose P A in a set of robot global reference coordinate systems W 2 is obtained. The LM algorithm is used to obtain the measurement coordinate system W 1 and the robot global The coordinate system transformation matrix T between the reference coordinate system W 2 . Secondly, for the target point set B to be compensated, a small number of target points are randomly selected to form a priori knowledge data set B p , and its pose P Bpm in the measurement coordinate system W 1 is measured, and the coordinate system is used to transform The matrix T obtains the actual pose P Bpr of the prior knowledge data set B p in the robot global reference coordinate system W 2 , and compares the actual pose P Bpr with its corresponding theoretical pose P Bpt . Obtain the corrected pose P Bpr ' of the corresponding target point in the prior knowledge data set B p , use the actual coordinates P Bpr of the target point in the prior knowledge data set B pre as input, and use its corresponding corrected position Pose P Bpr ' is used as the output to train a robot machining dynamic error compensation model established by a new support vector machine regression algorithm that integrates prior knowledge. Finally, the set composed of the remaining target points of the target point set B to be compensated, that is, the remaining data set C B B p , is used as input, and the model output established by the new support vector machine regression algorithm integrating prior knowledge is corrected. The target point set B' is input into the robot controller to compensate for the end error of the robot.
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。It is easy for those skilled in the art to understand that the above descriptions are only preferred embodiments of the present invention and are not intended to limit the present invention. Any modifications, equivalent substitutions and improvements, etc., made within the spirit and principles of the present invention, All should be included in the protection scope of the present invention.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010504997.XA CN111754567B (en) | 2020-06-05 | 2020-06-05 | Comprehensive compensation method for static and dynamic errors in grinding and polishing processing of aircraft composite member robot |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010504997.XA CN111754567B (en) | 2020-06-05 | 2020-06-05 | Comprehensive compensation method for static and dynamic errors in grinding and polishing processing of aircraft composite member robot |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111754567A CN111754567A (en) | 2020-10-09 |
| CN111754567B true CN111754567B (en) | 2024-02-27 |
Family
ID=72674799
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010504997.XA Active CN111754567B (en) | 2020-06-05 | 2020-06-05 | Comprehensive compensation method for static and dynamic errors in grinding and polishing processing of aircraft composite member robot |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111754567B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112643681A (en) * | 2020-12-16 | 2021-04-13 | 湖南涉外经济学院 | Intelligent path planning device and method for industrial mechanical arm |
| CN115493635B (en) * | 2022-10-12 | 2025-09-12 | 哈尔滨理工大学 | An experimental platform and adjustment method for hole position compensation and posture detection |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110385720A (en) * | 2019-07-26 | 2019-10-29 | 南京航空航天大学 | A kind of robot localization error compensating method based on deep neural network |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6812665B2 (en) * | 2002-04-19 | 2004-11-02 | Abb Ab | In-process relative robot workcell calibration |
-
2020
- 2020-06-05 CN CN202010504997.XA patent/CN111754567B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110385720A (en) * | 2019-07-26 | 2019-10-29 | 南京航空航天大学 | A kind of robot localization error compensating method based on deep neural network |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111754567A (en) | 2020-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109773786B (en) | Industrial robot plane precision calibration method | |
| CN110193829B (en) | Robot precision control method for coupling kinematics and rigidity parameter identification | |
| CN112873199B (en) | Robot absolute positioning precision calibration method based on kinematics and spatial interpolation | |
| CN106737855B (en) | A Robot Accuracy Compensation Method Combining Pose Error Model and Stiffness Compensation | |
| CN109822577B (en) | Mobile robot high-precision machining method based on visual servo | |
| CN108015808B (en) | A Kinematic Calibration Method for Hybrid Robots | |
| CN104608129B (en) | Robot Calibration Method Based on Plane Constraints | |
| CN104215206B (en) | Base coordinate calibration method of two-robot collaboration system | |
| CN112558547A (en) | Quick optimization method for geometric error compensation data of translational shaft of five-axis numerical control machine tool | |
| CN113580148B (en) | Kinematics Calibration Method for Parallel Robot Based on Equivalent Kinematic Chain | |
| CN110228067B (en) | Double-arm robot combined operation testing method based on laser tracker | |
| CN102183205A (en) | Method for matching optimal assembly poses of large-sized parts | |
| CN114952858B (en) | Industrial robot trajectory tracking method and system based on friction compensation control | |
| CN106777656B (en) | A PMPSD-based Absolute Precision Calibration Method for Industrial Robots | |
| CN109591019B (en) | An accurate spatial positioning method for objects without deterministic positioning features | |
| CN113359609A (en) | Key geometric error optimization proportioning compensation method for five-axis numerical control machine tool | |
| CN104890013A (en) | Pull-cord encoder based calibration method of industrial robot | |
| CN111754567B (en) | Comprehensive compensation method for static and dynamic errors in grinding and polishing processing of aircraft composite member robot | |
| CN115229796A (en) | A robot end pose tracking compensation method for the milling of large aerospace components | |
| CN105643620B (en) | A kind of industrial robot simple calibrating method based on cross rod piece | |
| CN115091455A (en) | Industrial robot positioning error compensation method | |
| CN111390914A (en) | Robot zero position and tool coordinate calibration method | |
| CN114754672A (en) | A method for determining the coordinate transformation relationship between a robot and a measuring instrument | |
| CN110470223A (en) | A kind of three dimensional space coordinate measurement method based on 1D displacement sensor | |
| CN112157654B (en) | Optimization method for positioning error of robot machining 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 |