WO2024011779A1 - Swing angle planning method for group-hole measurement, and readable medium and device - Google Patents

Swing angle planning method for group-hole measurement, and readable medium and device Download PDF

Info

Publication number
WO2024011779A1
WO2024011779A1 PCT/CN2022/125794 CN2022125794W WO2024011779A1 WO 2024011779 A1 WO2024011779 A1 WO 2024011779A1 CN 2022125794 W CN2022125794 W CN 2022125794W WO 2024011779 A1 WO2024011779 A1 WO 2024011779A1
Authority
WO
WIPO (PCT)
Prior art keywords
seed
measurement
cluster
normal vector
clustering
Prior art date
Application number
PCT/CN2022/125794
Other languages
French (fr)
Chinese (zh)
Inventor
张桂
高鑫
沈昕
姜振喜
毛一砚
王斌利
胡立
游莉萍
Original Assignee
成都飞机工业(集团)有限责任公司
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 成都飞机工业(集团)有限责任公司 filed Critical 成都飞机工业(集团)有限责任公司
Publication of WO2024011779A1 publication Critical patent/WO2024011779A1/en

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00Pattern recognition
    • G06F18/20Analysing
    • G06F18/23Clustering techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • 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

Definitions

  • the present invention relates to the technical field of on-machine detection. Specifically, it relates to a group hole measurement swing angle planning method, readable media and equipment, which are used to plan the measurement swing angle reasonably and efficiently during on-machine detection.
  • On-machine inspection can complete the measurement of part dimensions without disassembling the workpiece, thereby avoiding problems such as installation errors caused by secondary clamping, covering part deformation, loss of measurement efficiency, etc., to improve measurement quality and measurement efficiency, and in aviation manufacturing It has been increasingly widely used in the industry.
  • the stylus needs to enter the inside of the hole. Due to the poor openness in the hole, the angle between the stylus swing angle and the hole axis needs to be strictly controlled. If the angle is too large, measurement interference is likely to occur; if the angle is too small, , the measured swing angle will increase sharply and the measurement efficiency will decrease. Therefore, how to determine the measurement swing angle for on-machine inspection of hole features is a key technology for on-machine inspection.
  • the pendulum angle planning method is a brute force solution, that is, all possible measurement pendulum angles are listed in a row. out, and then optimize the swing angle according to the principle of "the most hole features that can be measured” until the swing angle planning of all hole features is completed.
  • this method can realize swing angle planning, its sequential calculation in an enumerated manner is still inefficient, and there is a possibility of omission or repeated calculation during the actual operation.
  • the object of the present invention is to solve the deficiencies in the existing technology and provide a method, readable medium and equipment for group hole measurement swing angle planning of hole characteristics, so as to solve the inconvenience of on-machine hole feature measurement swing angle calculation and group hole measurement swing angle calculation.
  • the problem of low on-machine detection efficiency caused by too many numbers by obtaining the hole feature normal vector set, stylus parameters and other parameters, calculate the angle threshold between the stylus axis and the hole axis when the hole and stylus do not interfere , and set the initial conditions and termination conditions of the clustering iteration; then use the hole feature normal vector set as the data set to perform clustering iteration, and obtain the measurement normal vector set. Stop the clustering iteration until the clustering converges, and convert the measurement normal vector set.
  • the planning process of measuring the swing angle is completed at this point; through the above method, the automatic planning of the swing angle measurement is realized, simplifying the operation steps and avoiding omissions or duplications in the planning process.
  • the present invention provides a group hole measurement swing angle planning method for on-machine hole feature detection, which includes the following steps:
  • Step S1 Obtain the hole feature normal vector set, stylus parameters, measurement depth and tolerance information
  • Step S2 Based on the stylus parameters, measurement depth and tolerance information, calculate the angle threshold between the stylus axis and the hole axis when the hole and stylus do not interfere; set the iterative clustering initial conditions, initialize the seed set, and set Iterative clustering termination conditions;
  • Step S3 Perform clustering: use the hole feature normal vector set as the data set, calculate the clustering set and the measurement normal vector set, and calculate candidate seeds;
  • Step S4 Determine the clustering convergence: if the clustering converges, execute step S6; if the clustering does not converge, execute step S5;
  • Step S5 Update the seed set first; then, if the current iteration number is less than or equal to the maximum iteration number, increase the current iteration number by one, and re-execute steps S3 and S4 in sequence; if the current iteration number is greater than the maximum iteration number, the current clustering Increase the scale by one, that is, add candidate seeds to the current seed set to form a new seed set, and re-execute steps S3 and S4;
  • Step S6 Convert the set of measured normal vectors into a set of measured pendulum angles, output the set of measured pendulum angles and merge to end the program.
  • the hole feature normal vector set is expressed as
  • the tolerance information is expressed as tol
  • the measurement depth is expressed as h
  • the angle threshold is expressed as T Agl
  • the cluster set is expressed as M c
  • the cluster scale is expressed as T
  • the number of cluster subsets in the cluster set M c is expressed by the cluster Determined by the scale T
  • the cluster subset is represented as M ci
  • the index i of the cluster subset M ci represents the sequence number in its cluster set M c
  • the seed set is represented as M s
  • the elements in the seed set M s are seeds.
  • the subscript j represents the seed
  • the sequence number in the seed set M s the seed in the seed set M s
  • the number is equal to the clustering scale T.
  • the seeds to be added to the seed set M s are candidate seeds and are represented by S new ; the current number of iterations is represented by n, and the maximum number of iterations is represented by T n ; the measurement normal vector set is represented by M n , and the measurement pendulum
  • the set of angles is denoted M A .
  • the stylus parameters include the diameter of the probe and the diameter of the probe.
  • the diameter of the side rod is represented by d and the diameter of the probe is represented by D.
  • the value of the tolerance information is greater than or equal to 0.5°.
  • the tolerance information tol side rod diameter d, probe diameter D. Measure the depth h and calculate the angle threshold T Agl through the following formula:
  • the initial clustering scale is 1, the current iteration number is 1; the initialization seed set is: the hole feature normal vector is centered on the first element As the only seed of the seed set; set the iterative clustering termination condition: determine the maximum number of iterations.
  • step S3 the cluster set and the measurement normal vector set are calculated through the following specific steps.
  • the specific steps are as follows:
  • Step S31 Add the angles formed by the i-th element in the hole feature normal vector set and various sub-elements in the seed set as elements into the i-th included angle set;
  • Step S32 If the j-th seed in the seed set forms an angle with the i-th element in the hole feature normal vector set, and is the element with the smallest value in the i-th angle set; then add the i-th element in the hole feature normal vector set to the cluster.
  • the jth cluster subset in the class set corresponds to adding the jth seed in the seed set to the measurement normal vector set; where, ⁇ i ⁇ Z
  • Step S33 Traverse the cluster set, and add the angle formed by each element of the j-th cluster subset in the cluster set and the j-th seed in the seed set to the candidate set as an element.
  • the corresponding cluster when the element in the candidate set reaches the maximum value is The specific elements of the jth cluster subset in the class set are used as candidate seeds; where ⁇ j ⁇ Z
  • the cluster subsets in the cluster set correspond to the seeds in the seed set in sequence.
  • the method for updating the seed set in step S5 is: converting each cluster subset in the cluster set into various sub-sets in the new seed set one by one; specifically, using each cluster subset
  • the average value of the elements in the class subset is sequentially used as various subsets in the new seed set, which can be expressed by the following formula:
  • the seed S j is expressed in the form of S j (x j , y j , z j ).
  • the subscript j represents the sequence number of the seed in the seed set.
  • index i (x j , y j , z j ) defines the seed S
  • the index j in the above formula is the same as The index i is equal.
  • the measurement normal vector set is M n
  • the elements of the measurement normal vector set M n are
  • Each element (A v , C v ) in the measured swing angle set M A is obtained by the following formula:
  • i represents the coordinate value of the measurement normal vector on the X-axis
  • j represents the coordinate value of the measurement normal vector on the Y-axis
  • k represents the coordinate value of the measurement normal vector on the Z-axis
  • (i, j, k) defines the measurement normal vector Proportional values in the X, Y, and Z directions
  • (A v , C v ) represents the measurement swing angle
  • a v represents the working angle of the A-axis during measurement
  • C v represents the working angle of the C-axis during measurement
  • V is the labeled measurement method vector
  • the present invention provides a computer-readable medium on which a computer program is stored.
  • the computer program is executed, the group hole measurement swing angle planning method described in any one of the above is implemented.
  • An electronic device including:
  • the processor is configured to execute any one of the above group hole measurement swing angle planning methods by executing the executable instructions.
  • the present invention optimizes the planning of hole feature measurement swing angles and clusters similar swing angles, thereby reducing measurement angles, saving probe calibration time, and improving on-machine detection efficiency.
  • Figure 1 is a schematic flow chart of the group hole measurement swing angle planning method provided by the present invention.
  • Figure 2 is a schematic diagram of the state when the probe is in the hole to be measured in the present invention.
  • This embodiment provides a group hole measurement swing angle planning method for on-machine hole feature detection, which specifically includes the following steps:
  • Step S1 Obtain the hole feature normal vector set, stylus parameters, measurement depth and tolerance information
  • Step S2 Based on the stylus parameters, measurement depth and tolerance information, calculate the angle threshold between the stylus axis and the hole axis when the hole and stylus do not interfere; set the iterative clustering initial conditions, initialize the seed set, and set Iterative clustering termination conditions;
  • Step S3 Perform clustering: use the hole feature normal vector set as the data set, calculate the clustering set and the measurement normal vector set, and calculate candidate seeds;
  • Step S4 Determine the clustering convergence: if the clustering converges, execute step S6; if the clustering does not converge, execute step S5;
  • Step S5 Update the seed set first; then, if the current iteration number is less than or equal to the maximum iteration number, increase the current iteration number by one, and re-execute steps S3 and S4 in sequence; if the current iteration number is greater than the maximum iteration number, the current clustering The scale is increased by one, that is, candidate seeds are added to the current seed set to form a new seed set, and then the current number of iterations returns to the initial value, and steps S3 and S4 are re-executed;
  • Step S6 Convert the set of measured normal vectors into a set of measured pendulum angles, output the set of measured pendulum angles and merge to end the program.
  • the hole feature normal vector set is expressed as Referring to Figure 2, the tolerance information is expressed as tol and the tolerance information is the allowed angular deviation range of the stylus axis, the measurement depth is expressed as h, and the angle threshold is expressed as T Agl ; the clustering set is expressed as M c, and the clustering set is expressed as M c .
  • the class scale is expressed as T.
  • the number of cluster subsets in the cluster set M c is determined by the cluster scale T.
  • the cluster subset is expressed as M ci .
  • the index i of the cluster subset M ci is expressed in its cluster set M
  • the sequence number in c the seed set is expressed as M s
  • the elements in the seed set M s are seeds
  • the seeds are expressed as
  • the subscript j represents the seed
  • the sequence number in the seed set M s the seed in the seed set M s
  • the number is equal to the clustering scale T.
  • the seeds to be added to the seed set M s are candidate seeds and are represented by S new ; the current number of iterations is represented by n, and the maximum number of iterations is represented by T n ; the measurement normal vector set is represented by M n , and the measurement pendulum
  • the set of angles is denoted M A .
  • the stylus parameters include the diameter of the probe and the diameter of the probe rod.
  • the diameter of the side rod is represented by d and the diameter of the probe is represented by D.
  • the value of the tolerance information is greater than or equal to 0.5°, that is, tol ⁇ 05°.
  • the initial value of the clustering scale is 1, and the initial value of the number of iterations is 1;
  • the initialization seed set is: use the first element in the hole feature normal vector as the seed The unique seed of the set; determines the maximum number of iterations followed by each cluster size, where the maximum number of iterations is 5.
  • step S3 the cluster set and the measurement normal vector set are calculated through the following specific steps.
  • the specific steps are as follows:
  • Step S31 Add the angles formed by the i-th element in the hole feature normal vector set and various sub-elements in the seed set as elements into the i-th included angle set;
  • Step S32 If the j-th seed in the seed set forms an angle with the i-th element in the hole feature normal vector set, and is the element with the smallest value in the i-th angle set; then add the i-th element in the hole feature normal vector set to the cluster.
  • the jth cluster subset in the class set corresponds to adding the jth seed in the seed set to the measurement normal vector set; where, ⁇ i ⁇ Z
  • Step S33 Traverse the cluster set, and add the angle formed by each element of the j-th cluster subset in the cluster set and the j-th seed in the seed set to the candidate set as an element.
  • the corresponding cluster when the element in the candidate set reaches the maximum value is The specific elements of the jth cluster subset in the class set are used as candidate seeds; where ⁇ j ⁇ Z
  • the cluster subsets in the cluster set and the seeds in the seed set correspond one to one in sequence.
  • the method for updating the seed set in step S5 is: converting each cluster subset in the cluster set into various subsets in the new seed set one by one; specifically, using the average value of the elements in each cluster subset As various children in the new seed set in turn, they can be expressed by the following formula:
  • the seed S j is expressed in the form of S j (x j , y j , z j ).
  • the subscript j represents the serial number of the seed in the seed set. (x j , y j , z j ) defines the seed.
  • i represents the coordinate value of the measurement normal vector on the X-axis
  • j represents the coordinate value of the measurement normal vector on the Y-axis
  • k represents the coordinate value of the measurement normal vector on the Z-axis
  • (i, j, k) defines the measurement normal vector Proportional values in the X, Y, and Z directions
  • (A v , C v ) represents the measurement swing angle
  • a v represents the working angle of the A-axis during measurement
  • C v represents the working angle of the C-axis during measurement
  • V is the labeled measurement method vector
  • This embodiment uses the group hole measurement swing angle planning method provided in Embodiment 1, and implements it under the following implementation conditions:
  • the machine tool structure is an AC swing angle five-axis machine tool, and the forming range is: A: -90° to 90°, C: -180° to 180°.
  • This embodiment provides a group hole measurement swing angle planning method, as shown in the flow chart of Figure 1, in which the specific implementation steps based on the implementation conditions provided by this embodiment are as follows:
  • Step S1 Obtain the hole feature normal vector set
  • Step S3 Execute clustering: use the hole feature normal vector set as the data set and the measured swing angle set as the target set, calculate the clustering set M c and the measured normal vector set M n , and calculate the candidate seed S new .
  • Step S31 Calculate the hole feature normal vector set elements in With various seeds in the seed set M s angle between is the first included angle;
  • Step S4 Determine the clustering convergence: If the maximum value T max of the elements in the candidate set is greater than the angle threshold Then the clustering does not converge, otherwise it converges; at this time, T max >, the clustering does not converge;
  • Step S10 Repeat steps S3 to S6 again,
  • the clustering converges, the loop exits, and M n is output.
  • Step S11 Convert the measurement normal vector set M n into the measurement pendulum angle set M A .
  • This set is the hole feature normal vector set The corresponding set of pendulum angles.
  • Step S1 Obtain the hole feature normal vector set
  • Step S3 Discretize the machine tool swing angle range according to 2T Agl to generate machine tool node A1:
  • Step S4 For every element in Solve the problem that the angle between A1 and A1 is not greater than Elements:
  • the angle between (0.0000, 0.0000, 1.0000) is not greater than 3.9272° (0.32, 1.24)
  • the angle between (0.0000, 0.0712, 0.9503) is not greater than 3.9272° (0.32, 1.24)
  • the angle between (0.0000, 0.0508, 0.9802) is not greater than 3.9272° (0.32, 1.24)
  • the middle element is the node with the most corresponding elements among all nodes in A1.
  • Step S6 Y1 is not empty, from Delete Y1;
  • Step S7 From Delete Y1 in If not empty, execute S4;
  • Step S9 Any element of A1 corresponds to 0 Medium element, Y1 is the empty set;
  • Step S10 The remaining nodes cannot measure the remaining holes, so the intervals need to be reduced and divided again to generate denser nodes;
  • Step S11 Execute step S3 to re-discretize the machine tool swing angle range according to T Agl to generate machine tool node A2:
  • Step S12 For every element in Solve for the elements in A2 whose angle is no greater than T Agl :
  • angles with (0.0000, 0.0000, 1.0000) are not greater than 3.9272°: (-0.25, 1.24) (3.33, 1.24), (-0.25, 5.34), (3.33, 5.34);
  • angles with (0.0000, 0.0712, 0.9503) are not greater than 3.9272°: (-0.25, 1.24), (3.33, 1.24), (-0.25, 5.34), (-3.43, 5.34);
  • angles with (0.0000, 0.0508, 0.9802) are not greater than 3.9272°: (-0.25, 1.24), (3.33, 1.24), (-0.25, 5.34), (-3.43, 5.34);
  • angles with (0.3536, -0.6124, 0.7071) are not greater than 3.9272°: (45.23, 28.22), (45.23, 31.25);
  • A21 (-0.25, 1.24) corresponds to 3
  • the middle element is the node with the most corresponding elements among all nodes in A2.
  • Step S14 Y2 is not empty, from Delete Y2;
  • Step S15 From Delete Y2 in If not empty, execute S12
  • Step S16 In A2, the angles with (0.3536, -0.6124, 0.7071) are not greater than 3.9272° including (45.23, 28.22), (45.23, 31.25);
  • the middle element is the node with the most corresponding elements among all nodes in A2.
  • Step S18 Y2 is not empty, from Delete Y2;
  • Step S19 From After deleting Y2 in If it is empty, the measurement and swing angle planning is completed.
  • Step S20 The final planned measured swing angle is
  • this embodiment greatly improves the accuracy of measurement angles through swing angle planning, reduces the number of required measured swing angles, and reduces the number of holes.
  • the planning time of measuring the swing angle is greatly reduced, and the hole characteristic measurement swing angle can be quickly obtained, which saves the measurement preparation time and improves the processing and measurement efficiency.

Abstract

A swing angle planning method for group-hole measurement, and a readable medium and a device. The method comprises: acquiring feature parameters of a hole to be subjected to measurement and feature parameters of a measurement probe; calculating a measurement included-angle threshold, and determining an initial condition for iteration and a termination condition for same; executing clustering, updating a cluster set and a measurement normal vector set, and calculating candidate seeds; updating a seed set, and performing iterative optimization; determining whether the number of iterations exceeds a limit; adding the candidate seeds to the seed set; performing iterative optimization until a cluster is converged; and converting the measurement normal vector set into a measurement swing angle set. The time for calibrating a measurement probe can be reduced, thereby improving the on-machine inspection efficiency.

Description

群孔测量摆角规划方法、可读介质和设备Group hole measurement swing angle planning method, readable media and equipment 技术领域Technical field
本发明涉及在机检测技术领域,具体的说,是群孔测量摆角规划方法、可读介质和设备,用于在机检测时对测量摆角进行合理高效地规划。The present invention relates to the technical field of on-machine detection. Specifically, it relates to a group hole measurement swing angle planning method, readable media and equipment, which are used to plan the measurement swing angle reasonably and efficiently during on-machine detection.
背景技术Background technique
在机检测能够在不拆卸工件的条件下完成对零件尺寸的测量,从而避免二次装夹带来的安装误差、包庇零件变形、测量效率损失等问题,以提高测量质量和测量效率,在航空制造业中得到了越来越广泛的应用。但在测量孔特征时,测针需要进入孔内部,由于孔内开敞性差,需严格控制测针摆角和孔轴线夹角,若夹角太大,容易发生测量干涉;若夹角太小,测量摆角将急剧增加,测量效率降低。因此,如何确定孔特征的在机检测的测量摆角,是在机检测的关键技术。On-machine inspection can complete the measurement of part dimensions without disassembling the workpiece, thereby avoiding problems such as installation errors caused by secondary clamping, covering part deformation, loss of measurement efficiency, etc., to improve measurement quality and measurement efficiency, and in aviation manufacturing It has been increasingly widely used in the industry. However, when measuring hole characteristics, the stylus needs to enter the inside of the hole. Due to the poor openness in the hole, the angle between the stylus swing angle and the hole axis needs to be strictly controlled. If the angle is too large, measurement interference is likely to occur; if the angle is too small, , the measured swing angle will increase sharply and the measurement efficiency will decrease. Therefore, how to determine the measurement swing angle for on-machine inspection of hole features is a key technology for on-machine inspection.
目前,在确定孔位测量摆角时,主要有如下两种方式:一是直接沿法式方向测量孔位,此方法确定摆角容易,但测量角度多,测针标定效率低,对结构复杂、孔数量庞大的零件,标定时间甚至远超加工时间与测量时间;二是对轴线相近的孔,大致判断可行的测量摆角,使用仿真软件进行干涉仿真,最终确定可行的摆角;此方法可有效减少测量角度,但需反复尝试与修正,效率较低,对于复杂零件,人工不可能完成。At present, when determining the hole position and measuring the swing angle, there are two main ways: one is to directly measure the hole position in the French direction. This method is easy to determine the swing angle, but there are many measurement angles and the stylus calibration efficiency is low. It is difficult for complex structures and For parts with a large number of holes, the calibration time even far exceeds the processing time and measurement time; second, for holes with similar axes, roughly determine the feasible measurement swing angle, use simulation software to perform interference simulation, and finally determine the feasible swing angle; this method can It can effectively reduce the measurement angle, but it requires repeated attempts and corrections, and the efficiency is low. For complex parts, it is impossible to complete manually.
如现有技术如公开号为CN112033331A的中国专利《一种基于三坐标测针的群孔测量摆角规划方法》所公开的摆角规划方法是一种暴力解法,即将所有可能的测量摆角列出,然后按照“能够测量的孔特征最多”为原则进行摆角优选,直到完成所有孔特征的摆角规划为止。该方法虽能实现摆角规划,但其以枚举的方式进行依次运算仍旧效率较低、且实际操作过程中存在遗漏或者重复计算的可能。For example, in the existing technology, such as the Chinese patent "A Group Hole Measurement Pendulum Angle Planning Method Based on Three-Coordinate Stylus" published by the Chinese patent No. CN112033331A, the pendulum angle planning method is a brute force solution, that is, all possible measurement pendulum angles are listed in a row. out, and then optimize the swing angle according to the principle of "the most hole features that can be measured" until the swing angle planning of all hole features is completed. Although this method can realize swing angle planning, its sequential calculation in an enumerated manner is still inefficient, and there is a possibility of omission or repeated calculation during the actual operation.
发明内容Contents of the invention
本发明的目的在于现有技术所存在的缺陷,提供一种孔特征的群孔测量摆角规划方法、可读介质和设备,解决在机检测孔特征测量摆角计算不便、群孔测量摆角数量过多而引起的在机检测效率低下的问题;通过获取孔特征法矢集合、测针参数等参数,计算出孔与测针不发生干涉时测针轴线与孔轴线之间的夹角阈值,并设置聚类迭代的初始条件和终止条件;继而以孔特征法矢集合为数据集进行聚类迭代,获得测量法矢集合,直至聚类收敛时停止聚类迭代,将测量发矢集合转换为测量摆角集合,至此,完成测量摆角的规划过程;通过上述方法,实现了测量摆角的自动规划,简化操作步骤同时避免了规划过程出现遗漏或重复。The object of the present invention is to solve the deficiencies in the existing technology and provide a method, readable medium and equipment for group hole measurement swing angle planning of hole characteristics, so as to solve the inconvenience of on-machine hole feature measurement swing angle calculation and group hole measurement swing angle calculation. The problem of low on-machine detection efficiency caused by too many numbers; by obtaining the hole feature normal vector set, stylus parameters and other parameters, calculate the angle threshold between the stylus axis and the hole axis when the hole and stylus do not interfere , and set the initial conditions and termination conditions of the clustering iteration; then use the hole feature normal vector set as the data set to perform clustering iteration, and obtain the measurement normal vector set. Stop the clustering iteration until the clustering converges, and convert the measurement normal vector set. In order to measure the swing angle set, the planning process of measuring the swing angle is completed at this point; through the above method, the automatic planning of the swing angle measurement is realized, simplifying the operation steps and avoiding omissions or duplications in the planning process.
本发明通过下述技术方案实现:The present invention is realized through the following technical solutions:
首先,本发明提供了一种群孔测量摆角规划方法,用于在机检测孔特征,包括以下步骤:First, the present invention provides a group hole measurement swing angle planning method for on-machine hole feature detection, which includes the following steps:
步骤S1:获取孔特征法矢集合、测针参数、测量深度和容差信息;Step S1: Obtain the hole feature normal vector set, stylus parameters, measurement depth and tolerance information;
步骤S2:根据测针参数、测量深度和容差信息,计算出孔与测针不发生干涉时测针轴线与孔轴线之间的夹角阈值;设置迭代聚类初始条件、初始化种子集合、设置迭代聚类终止条件;Step S2: Based on the stylus parameters, measurement depth and tolerance information, calculate the angle threshold between the stylus axis and the hole axis when the hole and stylus do not interfere; set the iterative clustering initial conditions, initialize the seed set, and set Iterative clustering termination conditions;
步骤S3:执行聚类:以孔特征法矢集合为数据集,计算聚类集合和测量法矢集合,并计算候选种子;Step S3: Perform clustering: use the hole feature normal vector set as the data set, calculate the clustering set and the measurement normal vector set, and calculate candidate seeds;
步骤S4:判断聚类收敛性:若聚类收敛,则执行步骤S6;若聚类不收敛,则执行步骤S5;Step S4: Determine the clustering convergence: if the clustering converges, execute step S6; if the clustering does not converge, execute step S5;
步骤S5:先更新种子集合;然后,若当前迭代次数小于等于最大迭代次数,则当前迭代次数增加一次,并重新依次执行步骤S3、步骤S4;若当前迭代次数大于最大迭代次数,则当前聚类规模加一,即增加候选种子至当前种子集合形成新的种子集合,并重新执行步骤S3、步骤S4;Step S5: Update the seed set first; then, if the current iteration number is less than or equal to the maximum iteration number, increase the current iteration number by one, and re-execute steps S3 and S4 in sequence; if the current iteration number is greater than the maximum iteration number, the current clustering Increase the scale by one, that is, add candidate seeds to the current seed set to form a new seed set, and re-execute steps S3 and S4;
步骤S6:将测量法矢集合转化为测量摆角集合,输出测量摆角集合并结束程序。Step S6: Convert the set of measured normal vectors into a set of measured pendulum angles, output the set of measured pendulum angles and merge to end the program.
本发明中,孔特征法矢集合表示为
Figure PCTCN2022125794-appb-000001
容差信息表示为tol、测量深度表示为h、夹角阈值表示为T Agl;聚类集合表示为M c,聚类规模表示为T,聚类集合M c中聚类子集数量由聚类规模T决定,聚类子集表示为M ci,聚类子集M ci的角标i表示在其聚类集合M c中的序号;种子集合表示为M s,种子集合M s中元素为种子,并将种子表示为
Figure PCTCN2022125794-appb-000002
其中角标j表示种子
Figure PCTCN2022125794-appb-000003
在种子集合M s中的序号,种子集合M s中种子
Figure PCTCN2022125794-appb-000004
数量与聚类规模T相等,待加入种子集合M s中种子为候选种子并用S new表示;当前迭代次数表示为n,最大迭代次数表示为T n;测量法矢集合表示为M n,测量摆角集合表示为M A
In the present invention, the hole feature normal vector set is expressed as
Figure PCTCN2022125794-appb-000001
The tolerance information is expressed as tol, the measurement depth is expressed as h, and the angle threshold is expressed as T Agl ; the cluster set is expressed as M c , the cluster scale is expressed as T, and the number of cluster subsets in the cluster set M c is expressed by the cluster Determined by the scale T, the cluster subset is represented as M ci , and the index i of the cluster subset M ci represents the sequence number in its cluster set M c ; the seed set is represented as M s , and the elements in the seed set M s are seeds. , and express the seed as
Figure PCTCN2022125794-appb-000002
The subscript j represents the seed
Figure PCTCN2022125794-appb-000003
The sequence number in the seed set M s , the seed in the seed set M s
Figure PCTCN2022125794-appb-000004
The number is equal to the clustering scale T. The seeds to be added to the seed set M s are candidate seeds and are represented by S new ; the current number of iterations is represented by n, and the maximum number of iterations is represented by T n ; the measurement normal vector set is represented by M n , and the measurement pendulum The set of angles is denoted M A .
为了更好地完善本发明,进一步优化地,所述测针参数包括测头直径和测杆直径,侧杆直径表示为d、测头直径表示为D。In order to better improve the present invention, further optimally, the stylus parameters include the diameter of the probe and the diameter of the probe. The diameter of the side rod is represented by d and the diameter of the probe is represented by D.
为了更好地完善本发明,进一步优化地,所述容差信息的取值大于等于0.5°为了更好地完善本发明,进一步优化地,根据容差信息tol、侧杆直径d、测头直径D、测量深度h,通过以下公式计算夹角阈值T AglIn order to better perfect the present invention, further optimally, the value of the tolerance information is greater than or equal to 0.5°. In order to better perfect the present invention, further optimally, according to the tolerance information tol, side rod diameter d, probe diameter D. Measure the depth h and calculate the angle threshold T Agl through the following formula:
Figure PCTCN2022125794-appb-000005
Figure PCTCN2022125794-appb-000005
为了更好地完善本发明,进一步优化地,设置所述迭代聚类初始条件为:初始聚类 规模为1、当前迭代次数为1;初始化种子集合为:将孔特征法矢中第一个元素作为种子集合的唯一种子;设置所述迭代聚类终止条件:确定最大迭代次数。In order to better improve the present invention, further optimally set the initial conditions of the iterative clustering as follows: the initial clustering scale is 1, the current iteration number is 1; the initialization seed set is: the hole feature normal vector is centered on the first element As the only seed of the seed set; set the iterative clustering termination condition: determine the maximum number of iterations.
为了更好地完善本发明,进一步优化地,在步骤S3中,通过以下具体步骤计算得到聚类集合和测量法矢集合,具体步骤如下:In order to better improve the present invention, in further optimization, in step S3, the cluster set and the measurement normal vector set are calculated through the following specific steps. The specific steps are as follows:
步骤S31:以孔特征法矢集合中第i元素与种子集合中各种子分别形成的夹角,加入第i夹角集合中作为元素;Step S31: Add the angles formed by the i-th element in the hole feature normal vector set and various sub-elements in the seed set as elements into the i-th included angle set;
步骤S32:若种子集合中第j种子,与孔特征法矢集合中第i元素形成夹角,为第i夹角集合取值最小的元素;则将孔特征法矢集合中第i元素加入聚类集合中第j聚类子集,对应将种子集合中第j种子加入测量法矢集合;其中,{i∈Z|1≤i≤N},N为孔特征法矢集合中元素的数量;Step S32: If the j-th seed in the seed set forms an angle with the i-th element in the hole feature normal vector set, and is the element with the smallest value in the i-th angle set; then add the i-th element in the hole feature normal vector set to the cluster. The jth cluster subset in the class set corresponds to adding the jth seed in the seed set to the measurement normal vector set; where, {i∈Z|1≤i≤N}, N is the number of elements in the hole feature normal vector set;
步骤S33:遍历聚类集合,将聚类集合中第j聚类子集的各元素与种子集合中第j种子形成的夹角加入候选集合作为元素,候选集合中元素取得最大值时对应的聚类集合中第j个聚类子集的特定元素作为候选种子;其中,{j∈Z|1≤j≤M},M为种子集合中种子的数量。Step S33: Traverse the cluster set, and add the angle formed by each element of the j-th cluster subset in the cluster set and the j-th seed in the seed set to the candidate set as an element. The corresponding cluster when the element in the candidate set reaches the maximum value is The specific elements of the jth cluster subset in the class set are used as candidate seeds; where {j∈Z|1≤j≤M}, M is the number of seeds in the seed set.
为了更好地完善本发明,进一步优化地,所述聚类集合中聚类子集与种子集合中种子顺次一一对应。In order to better improve the present invention, further optimally, the cluster subsets in the cluster set correspond to the seeds in the seed set in sequence.
为了更好地完善本发明,进一步优化地,当所述候选集合中元素取得的最大值大于夹角阈值,则判定聚类不收敛,否则判定聚类收敛。In order to better improve the present invention, it is further optimized that when the maximum value obtained by the elements in the candidate set is greater than the angle threshold, it is determined that the clustering does not converge, otherwise it is determined that the clustering has converged.
为了更好地完善本发明,进一步优化地,步骤S5中更新种子集合的方法为:将聚类集合中各聚类子集一一转换为新的种子集合中各种子;具体为使用各聚类子集中元素平均值依次作为新的种子集合中的各种子,可用如下公式表示:
Figure PCTCN2022125794-appb-000006
此式中种子S j以S j(x j,y j,z j)的形式表示,角标j表示该种子在种子集合中的序号,(x j,y j,z j)定义了种子S j在空间直角坐标系X、Y、Z方向上的比例值;聚类子集M ci中,角标i表示该聚类子集在聚类集合中的序号,上式中的角标j与角标i相等。
In order to better improve the present invention, further optimally, the method for updating the seed set in step S5 is: converting each cluster subset in the cluster set into various sub-sets in the new seed set one by one; specifically, using each cluster subset The average value of the elements in the class subset is sequentially used as various subsets in the new seed set, which can be expressed by the following formula:
Figure PCTCN2022125794-appb-000006
In this formula, the seed S j is expressed in the form of S j (x j , y j , z j ). The subscript j represents the sequence number of the seed in the seed set. (x j , y j , z j ) defines the seed S The proportional value of j in the X, Y, and Z directions of the space rectangular coordinate system; in the cluster subset M ci , the index i represents the serial number of the cluster subset in the cluster set. The index j in the above formula is the same as The index i is equal.
为了更好地完善本发明,进一步优化地,其中测量法矢集合为M n,测量法矢集合M n的元素为
Figure PCTCN2022125794-appb-000007
所述测量摆角集合M A中各元素(A v,C v)通过以下公式求得:
In order to better improve the present invention, in a further optimized manner, the measurement normal vector set is M n , and the elements of the measurement normal vector set M n are
Figure PCTCN2022125794-appb-000007
Each element (A v , C v ) in the measured swing angle set M A is obtained by the following formula:
Figure PCTCN2022125794-appb-000008
Figure PCTCN2022125794-appb-000008
Figure PCTCN2022125794-appb-000009
其中,式中
Figure PCTCN2022125794-appb-000010
表示测量法矢,i表示测量法矢在X轴坐标值,j表示测量法矢在Y轴坐标值,k表示测量法矢在Z轴坐标值,(i,j,k)定义了测量法矢
Figure PCTCN2022125794-appb-000011
在X、Y、Z方向上的比例值;(A v,C v)表示测量摆角,A v表示测量时A轴工作角度,C v表示测量时C轴工作角度;V为标注测量法矢和测量摆角序号的角标,测量法矢和测量摆角依顺序一一对应。
Figure PCTCN2022125794-appb-000009
Among them, in the formula
Figure PCTCN2022125794-appb-000010
represents the measurement normal vector, i represents the coordinate value of the measurement normal vector on the X-axis, j represents the coordinate value of the measurement normal vector on the Y-axis, k represents the coordinate value of the measurement normal vector on the Z-axis, (i, j, k) defines the measurement normal vector
Figure PCTCN2022125794-appb-000011
Proportional values in the X, Y, and Z directions; (A v , C v ) represents the measurement swing angle, A v represents the working angle of the A-axis during measurement, C v represents the working angle of the C-axis during measurement; V is the labeled measurement method vector There is a one-to-one correspondence with the corner mark of the measured pendulum angle serial number, the measured normal vector and the measured pendulum angle in sequence.
其次,本发明提供了一种计算机可读介质,其上存储有计算机程序,所述计算机程序被执行时实现以上任一项所述的群孔测量摆角规划方法。Secondly, the present invention provides a computer-readable medium on which a computer program is stored. When the computer program is executed, the group hole measurement swing angle planning method described in any one of the above is implemented.
一种电子设备,包括:An electronic device including:
处理器;processor;
存储器,用于存储所述处理器的可执行指令;memory for storing executable instructions for the processor;
其中,所述处理器配置为经由执行所述可执行指令来执行以上任一项所述的群孔测量摆角规划方法。Wherein, the processor is configured to execute any one of the above group hole measurement swing angle planning methods by executing the executable instructions.
有益效果:本发明优化规划孔特征测量摆角,并将相近摆角进行聚类,减少测量角度,节省测头标定时间,提高在机检测效率。Beneficial effects: The present invention optimizes the planning of hole feature measurement swing angles and clusters similar swing angles, thereby reducing measurement angles, saving probe calibration time, and improving on-machine detection efficiency.
附图说明Description of drawings
下面将结合附图对技术方案进行清楚、完整地描述,显然所描述的实施例是本发明一部分实施例,而不是全部的实施例。The technical solution will be clearly and completely described below in conjunction with the accompanying drawings. It is obvious that the described embodiments are some, not all, of the embodiments of the present invention.
图1为本发明提供的群孔测量摆角规划方法的流程示意图;Figure 1 is a schematic flow chart of the group hole measurement swing angle planning method provided by the present invention;
图2为本发明中测针在被测孔内时的状态示意图。Figure 2 is a schematic diagram of the state when the probe is in the hole to be measured in the present invention.
具体实施方式Detailed ways
以下结合实施例的具体实施方式,对本发明创造的上述内容再做进一步的详细说明。但不应将此理解为本发明上述主题的范围仅限于以下的实例。在不脱离本发明上述技术思想情况下,根据本领域普通技术知识和惯用手段作出的各种替换或变更,均应包括在本发明的范围内。The above content created by the present invention will be further described in detail below in conjunction with the specific implementation of the embodiments. However, this should not be understood to mean that the scope of the above subject matter of the present invention is limited to the following examples. Without departing from the above technical ideas of the present invention, various substitutions or changes made based on common technical knowledge and conventional means in the art should be included in the scope of the present invention.
以下实施例通过将本发明的方法与现有技术所公开的方法进行对比进行演示说明。The following examples demonstrate the method of the present invention by comparing it with methods disclosed in the prior art.
实施例1:Example 1:
本实施例提供一种群孔测量摆角规划方法,用于在机检测孔特征,具体包括以下步骤:This embodiment provides a group hole measurement swing angle planning method for on-machine hole feature detection, which specifically includes the following steps:
步骤S1:获取孔特征法矢集合、测针参数、测量深度和容差信息;Step S1: Obtain the hole feature normal vector set, stylus parameters, measurement depth and tolerance information;
步骤S2:根据测针参数、测量深度和容差信息,计算出孔与测针不发生干涉时测针轴线与孔轴线之间的夹角阈值;设置迭代聚类初始条件、初始化种子集合、设置迭代聚类终止条件;Step S2: Based on the stylus parameters, measurement depth and tolerance information, calculate the angle threshold between the stylus axis and the hole axis when the hole and stylus do not interfere; set the iterative clustering initial conditions, initialize the seed set, and set Iterative clustering termination conditions;
步骤S3:执行聚类:以孔特征法矢集合为数据集,计算聚类集合和测量法矢集合,并计算候选种子;Step S3: Perform clustering: use the hole feature normal vector set as the data set, calculate the clustering set and the measurement normal vector set, and calculate candidate seeds;
步骤S4:判断聚类收敛性:若聚类收敛,则执行步骤S6;若聚类不收敛,则执行步骤S5;Step S4: Determine the clustering convergence: if the clustering converges, execute step S6; if the clustering does not converge, execute step S5;
步骤S5:先更新种子集合;然后,若当前迭代次数小于等于最大迭代次数,则当前迭代次数增加一次,并重新依次执行步骤S3、步骤S4;若当前迭代次数大于最大迭代次数,则当前聚类规模加一,即增加候选种子至当前种子集合形成新的种子集合,然后当前迭代次数回归初始值,并重新执行步骤S3、步骤S4;Step S5: Update the seed set first; then, if the current iteration number is less than or equal to the maximum iteration number, increase the current iteration number by one, and re-execute steps S3 and S4 in sequence; if the current iteration number is greater than the maximum iteration number, the current clustering The scale is increased by one, that is, candidate seeds are added to the current seed set to form a new seed set, and then the current number of iterations returns to the initial value, and steps S3 and S4 are re-executed;
步骤S6:将测量法矢集合转化为测量摆角集合,输出测量摆角集合并结束程序。Step S6: Convert the set of measured normal vectors into a set of measured pendulum angles, output the set of measured pendulum angles and merge to end the program.
本发明中,孔特征法矢集合表示为
Figure PCTCN2022125794-appb-000012
参阅图2容差信息表示为tol且该容差信息为示测针轴线所被允许的角度偏离范围、测量深度表示为h、夹角阈值表示为T Agl;聚类集合表示为M c,聚类规模表示为T,聚类集合M c中聚类子集数量由聚类规模T决定,聚类子集表示为M ci,聚类子集M ci的角标i表示在其聚类集合M c中的序号;种子集合表示为M s,种子集合M s中元素为种子,并将种子表示为
Figure PCTCN2022125794-appb-000013
其中角标j表示种子
Figure PCTCN2022125794-appb-000014
在种子集合M s中的序号,种子集合M s中种子
Figure PCTCN2022125794-appb-000015
数量与聚类规模T相等,待加入种子集合M s中种子为候选种子并用S new表示;当前迭代次数表示为n,最大迭代次数表示为T n;测量法矢集合表示为M n,测量摆角集合表示为M A
In the present invention, the hole feature normal vector set is expressed as
Figure PCTCN2022125794-appb-000012
Referring to Figure 2, the tolerance information is expressed as tol and the tolerance information is the allowed angular deviation range of the stylus axis, the measurement depth is expressed as h, and the angle threshold is expressed as T Agl ; the clustering set is expressed as M c, and the clustering set is expressed as M c . The class scale is expressed as T. The number of cluster subsets in the cluster set M c is determined by the cluster scale T. The cluster subset is expressed as M ci . The index i of the cluster subset M ci is expressed in its cluster set M The sequence number in c ; the seed set is expressed as M s , the elements in the seed set M s are seeds, and the seeds are expressed as
Figure PCTCN2022125794-appb-000013
The subscript j represents the seed
Figure PCTCN2022125794-appb-000014
The sequence number in the seed set M s , the seed in the seed set M s
Figure PCTCN2022125794-appb-000015
The number is equal to the clustering scale T. The seeds to be added to the seed set M s are candidate seeds and are represented by S new ; the current number of iterations is represented by n, and the maximum number of iterations is represented by T n ; the measurement normal vector set is represented by M n , and the measurement pendulum The set of angles is denoted M A .
进一步完善本实施例:所述测针参数包括测头直径和测杆直径,侧杆直径表示为d、测头直径表示为D。To further improve this embodiment: the stylus parameters include the diameter of the probe and the diameter of the probe rod. The diameter of the side rod is represented by d and the diameter of the probe is represented by D.
进一步完善本实施例:所述容差信息的取值大于等于0.5°,即tol≥05°。To further improve this embodiment: the value of the tolerance information is greater than or equal to 0.5°, that is, tol≥05°.
进一步完善本实施例:请参阅图2,根据容差信息tol、侧杆直径d、测头直径D、测量深度h,通过以下公式计算夹角阈值T AglFurther improve this embodiment: Please refer to Figure 2. According to the tolerance information tol, side rod diameter d, probe diameter D, and measurement depth h, calculate the angle threshold T Agl through the following formula:
Figure PCTCN2022125794-appb-000016
Figure PCTCN2022125794-appb-000016
进一步完善本实施例:设置所述迭代聚类初始条件为:聚类规模的初始值为1、迭代次数的初始值为1;初始化种子集合为:将孔特征法矢中第一个元素作为种子集合的唯一种子;确定各聚类规模共同遵循的最大迭代次数,其中最大迭代次数为5。Further improve this embodiment: set the initial conditions of the iterative clustering as follows: the initial value of the clustering scale is 1, and the initial value of the number of iterations is 1; the initialization seed set is: use the first element in the hole feature normal vector as the seed The unique seed of the set; determines the maximum number of iterations followed by each cluster size, where the maximum number of iterations is 5.
进一步完善本实施例:在步骤S3中,通过以下具体步骤计算得到聚类集合和测量法矢集合,具体步骤如下:To further improve this embodiment: in step S3, the cluster set and the measurement normal vector set are calculated through the following specific steps. The specific steps are as follows:
步骤S31:以孔特征法矢集合中第i元素与种子集合中各种子分别形成的夹角,加入第i夹角集合中作为元素;Step S31: Add the angles formed by the i-th element in the hole feature normal vector set and various sub-elements in the seed set as elements into the i-th included angle set;
步骤S32:若种子集合中第j种子,与孔特征法矢集合中第i元素形成夹角,为第i夹角集合取值最小的元素;则将孔特征法矢集合中第i元素加入聚类集合中第j聚类子集,对应将种子集合中第j种子加入测量法矢集合;其中,{i∈Z|1≤i≤N},N为孔特征法矢集合中元素的数量;Step S32: If the j-th seed in the seed set forms an angle with the i-th element in the hole feature normal vector set, and is the element with the smallest value in the i-th angle set; then add the i-th element in the hole feature normal vector set to the cluster. The jth cluster subset in the class set corresponds to adding the jth seed in the seed set to the measurement normal vector set; where, {i∈Z|1≤i≤N}, N is the number of elements in the hole feature normal vector set;
步骤S33:遍历聚类集合,将聚类集合中第j聚类子集的各元素与种子集合中第j种子形成的夹角加入候选集合作为元素,候选集合中元素取得最大值时对应的聚类集合中第j个聚类子集的特定元素作为候选种子;其中,{j∈Z|1≤j≤M},M为种子集合中种子的数量。Step S33: Traverse the cluster set, and add the angle formed by each element of the j-th cluster subset in the cluster set and the j-th seed in the seed set to the candidate set as an element. The corresponding cluster when the element in the candidate set reaches the maximum value is The specific elements of the jth cluster subset in the class set are used as candidate seeds; where {j∈Z|1≤j≤M}, M is the number of seeds in the seed set.
进一步完善本实施例:所述聚类集合中聚类子集与种子集合中种子顺次一一对应。To further improve this embodiment: the cluster subsets in the cluster set and the seeds in the seed set correspond one to one in sequence.
进一步完善本实施例:当所述候选集合中元素取得的最大值大于夹角阈值,则判定聚类不收敛,否则判定聚类收敛。To further improve this embodiment: when the maximum value obtained by the elements in the candidate set is greater than the angle threshold, it is determined that the clustering does not converge, otherwise it is determined that the clustering has converged.
进一步完善本实施例:步骤S5中更新种子集合的方法为:将聚类集合中各聚类子集一一转换为新的种子集合中各种子;具体为使用各聚类子集中元素平均值依次作为新的种子集合中的各种子,可用如下公式表示:To further improve this embodiment: the method for updating the seed set in step S5 is: converting each cluster subset in the cluster set into various subsets in the new seed set one by one; specifically, using the average value of the elements in each cluster subset As various children in the new seed set in turn, they can be expressed by the following formula:
Figure PCTCN2022125794-appb-000017
此上式中种子S j以S j(x j,y j,z j)的形式表示,角标j表示该种子在种子集合中的序号,(x j,y j,z j)定义了种子S j在空间直角坐标系X、Y、Z方向上的比例值;聚类子集M ci中,角标i表示该聚类子集在聚类集合中的序号,上式中的角标j与角标i相等。
Figure PCTCN2022125794-appb-000017
In the above formula, the seed S j is expressed in the form of S j (x j , y j , z j ). The subscript j represents the serial number of the seed in the seed set. (x j , y j , z j ) defines the seed. The proportional value of S j in the X, Y, and Z directions of the space rectangular coordinate system; in the cluster subset M ci , the index i represents the serial number of the cluster subset in the cluster set, and the index j in the above formula It is equal to the index i.
进一步完善本实施例:其中测量法矢集合为M n,测量法矢集合M n的元素为
Figure PCTCN2022125794-appb-000018
所述测量摆角集合M A中各元素(A v,C v)通过以下公式求得:
Further improve this embodiment: where the measurement normal vector set is M n , and the elements of the measurement normal vector set M n are
Figure PCTCN2022125794-appb-000018
Each element (A v , C v ) in the measured swing angle set M A is obtained by the following formula:
Figure PCTCN2022125794-appb-000019
Figure PCTCN2022125794-appb-000019
Figure PCTCN2022125794-appb-000020
其中,上式中
Figure PCTCN2022125794-appb-000021
表示测量法矢,i表示测量法矢在X轴坐标值,j表示测量法矢在Y轴坐标值,k表示测量法矢在Z轴坐标值,(i,j,k)定义了测量法矢
Figure PCTCN2022125794-appb-000022
在X、Y、Z方向上的比例值;(A v,C v)表示测量 摆角,A v表示测量时A轴工作角度,C v表示测量时C轴工作角度;V为标注测量法矢和测量摆角序号的角标,测量法矢和测量摆角依顺序一一对应。
Figure PCTCN2022125794-appb-000020
Among them, in the above formula
Figure PCTCN2022125794-appb-000021
represents the measurement normal vector, i represents the coordinate value of the measurement normal vector on the X-axis, j represents the coordinate value of the measurement normal vector on the Y-axis, k represents the coordinate value of the measurement normal vector on the Z-axis, (i, j, k) defines the measurement normal vector
Figure PCTCN2022125794-appb-000022
Proportional values in the X, Y, and Z directions; (A v , C v ) represents the measurement swing angle, A v represents the working angle of the A-axis during measurement, C v represents the working angle of the C-axis during measurement; V is the labeled measurement method vector There is a one-to-one correspondence with the corner mark of the measured pendulum angle serial number, the measured normal vector and the measured pendulum angle in sequence.
实施例2:Example 2:
本实施例利用实施例1所提供群孔测量摆角规划方法,并在如下实施条件下进行具体实施:This embodiment uses the group hole measurement swing angle planning method provided in Embodiment 1, and implements it under the following implementation conditions:
①孔特征法矢集合
Figure PCTCN2022125794-appb-000023
① Hole feature normal vector set
Figure PCTCN2022125794-appb-000023
Figure PCTCN2022125794-appb-000024
Figure PCTCN2022125794-appb-000024
②测量深度h=3mm、容差信息tol=0.5°、测杆直径d=1.5mm、测头直径D=2mm;②Measurement depth h=3mm, tolerance information tol=0.5°, rod diameter d=1.5mm, probe diameter D=2mm;
③机床结构为AC摆角五轴机床,形成范围为:A:-90°至90°,C:-180°至180°。③The machine tool structure is an AC swing angle five-axis machine tool, and the forming range is: A: -90° to 90°, C: -180° to 180°.
本实施例提供了群孔测量摆角规划方法,如图1所示的流程图,其中以本实施例提供的实施条件进行具体实施操作步骤如下:This embodiment provides a group hole measurement swing angle planning method, as shown in the flow chart of Figure 1, in which the specific implementation steps based on the implementation conditions provided by this embodiment are as follows:
步骤S1:获取孔特征法矢集合
Figure PCTCN2022125794-appb-000025
Step S1: Obtain the hole feature normal vector set
Figure PCTCN2022125794-appb-000025
Figure PCTCN2022125794-appb-000026
Figure PCTCN2022125794-appb-000026
获取测量深度h=3mm、容差信息tol=0.5°,获取测针参数,测针参数包括测杆直径d=1.5mm、测头直径D=2mm;Obtain the measurement depth h=3mm, the tolerance information tol=0.5°, and obtain the stylus parameters. The stylus parameters include the diameter of the measuring rod d=1.5mm and the diameter of the measuring head D=2mm;
步骤S2:设置迭代聚类初始条件:初始聚类规模T=1,当前迭代次数n=1;初始化种子集合M s:将孔特征法矢
Figure PCTCN2022125794-appb-000027
中第一个元素作为种子集合M s的唯一种子
Figure PCTCN2022125794-appb-000028
可以表示为
Figure PCTCN2022125794-appb-000029
Figure PCTCN2022125794-appb-000030
设置迭代聚类终止条件:各聚类规模下的共同遵循的最大迭代次数T n=5;计算测量不干涉条件下测针轴线与孔轴线的夹角阈值T Agl,得到T Agl=3.9272°。
Step S2: Set the initial conditions for iterative clustering: initial clustering scale T = 1, current iteration number n = 1; initialize seed set M s : change the hole feature normal vector
Figure PCTCN2022125794-appb-000027
The first element in is used as the only seed of the seed set M s
Figure PCTCN2022125794-appb-000028
It can be expressed as
Figure PCTCN2022125794-appb-000029
Figure PCTCN2022125794-appb-000030
Set the iterative clustering termination condition: the maximum number of iterations T n =5 that is common to each cluster scale; calculate the angle threshold T Agl between the stylus axis and the hole axis under the condition of non-interference measurement, and obtain T Agl =3.9272°.
步骤S3:执行聚类:以孔特征法矢集合为数据集,以测量摆角集合为目标集,计算聚类集合M c与测量法矢集合M n,并计算候选种子S newStep S3: Execute clustering: use the hole feature normal vector set as the data set and the measured swing angle set as the target set, calculate the clustering set M c and the measured normal vector set M n , and calculate the candidate seed S new .
步骤S31:计算孔特征法矢集合
Figure PCTCN2022125794-appb-000031
中各元素
Figure PCTCN2022125794-appb-000032
与种子集合M s中各种子
Figure PCTCN2022125794-appb-000033
之间的夹角
Figure PCTCN2022125794-appb-000034
为第一夹角;
Step S31: Calculate the hole feature normal vector set
Figure PCTCN2022125794-appb-000031
elements in
Figure PCTCN2022125794-appb-000032
With various seeds in the seed set M s
Figure PCTCN2022125794-appb-000033
angle between
Figure PCTCN2022125794-appb-000034
is the first included angle;
Figure PCTCN2022125794-appb-000035
Figure PCTCN2022125794-appb-000035
此时存在种子结集合M s中唯一种子
Figure PCTCN2022125794-appb-000036
与孔特征法矢集合
Figure PCTCN2022125794-appb-000037
中各元素的夹角取得第i夹角集合中取值最小的元素,表示为:
Figure PCTCN2022125794-appb-000038
则将
Figure PCTCN2022125794-appb-000039
加入到聚类集合M c的第一个聚类子集M c1,将
Figure PCTCN2022125794-appb-000040
加入到测量法矢集合M n
At this time, there is a unique seed in the seed set M s
Figure PCTCN2022125794-appb-000036
Set with hole feature normal vectors
Figure PCTCN2022125794-appb-000037
The angle between each element in is obtained as the element with the smallest value in the i-th angle set, which is expressed as:
Figure PCTCN2022125794-appb-000038
then will
Figure PCTCN2022125794-appb-000039
The first cluster subset M c1 added to the cluster set M c will
Figure PCTCN2022125794-appb-000040
Add to the measurement method vector set M n ;
Figure PCTCN2022125794-appb-000041
Figure PCTCN2022125794-appb-000041
步骤S33:遍历聚类集M c,聚类子集中
Figure PCTCN2022125794-appb-000042
各元素到对应种子
Figure PCTCN2022125794-appb-000043
的夹角加入候选集合作为元素,候选集合的元素
Figure PCTCN2022125794-appb-000044
取得最大值为T max=45.0021°,并从取得最大值T max的元素
Figure PCTCN2022125794-appb-000045
Figure PCTCN2022125794-appb-000046
取出作为候选种子S new
Step S33: Traverse the cluster set M c and select the cluster subset
Figure PCTCN2022125794-appb-000042
Each element to the corresponding seed
Figure PCTCN2022125794-appb-000043
The angle between is added to the candidate set as an element, and the elements of the candidate set
Figure PCTCN2022125794-appb-000044
Obtain the maximum value T max =45.0021°, and obtain the element with the maximum value T max from
Figure PCTCN2022125794-appb-000045
middle
Figure PCTCN2022125794-appb-000046
Take out S new as the candidate seed,
Figure PCTCN2022125794-appb-000047
Figure PCTCN2022125794-appb-000047
步骤S4:判断聚类收敛性:若候选集合中的元素最大值T max大于夹角阈值
Figure PCTCN2022125794-appb-000048
则聚类不收敛,反之则收敛;此时T max>,聚类不收敛;
Step S4: Determine the clustering convergence: If the maximum value T max of the elements in the candidate set is greater than the angle threshold
Figure PCTCN2022125794-appb-000048
Then the clustering does not converge, otherwise it converges; at this time, T max >, the clustering does not converge;
步骤S5:聚类不收敛则更新种子集合:更新后的种子集合M s中唯一种子
Figure PCTCN2022125794-appb-000049
Figure PCTCN2022125794-appb-000050
当前迭代次数n=1≤5=T n,则迭代次数加1,记迭代次数n=n+1=2,开始在聚类规模T=1不变的情况下的重复执行步骤S3至步骤S4,直到n=6时:
Step S5: If the clustering does not converge, update the seed set: the only seed in the updated seed set M s
Figure PCTCN2022125794-appb-000049
Figure PCTCN2022125794-appb-000050
If the current number of iterations n=1≤5= Tn , then add 1 to the number of iterations, record the number of iterations n=n+1=2, and start to repeat steps S3 to S4 with the clustering scale T=1 unchanged. , until n=6:
Figure PCTCN2022125794-appb-000051
Figure PCTCN2022125794-appb-000051
Figure PCTCN2022125794-appb-000052
Figure PCTCN2022125794-appb-000052
T max=34.269°;S new=(0.3536,-0.6124,0.7071) T max =34.269°; S new =(0.3536, -0.6124, 0.7071)
聚类仍不收敛,且此时n>T n。此时,令聚类规模加一,记新的聚类规模T=T+1=2,增加候选种子S new至原种子集合M s,则新的种子集合M s=M s+S newThe clustering still does not converge, and n>T n at this time. At this time, increase the clustering scale by one, record the new clustering scale T=T+1=2, and add candidate seeds S new to the original seed set M s , then the new seed set M s =M s +S new :
Figure PCTCN2022125794-appb-000053
Figure PCTCN2022125794-appb-000053
重新执行步骤S3至步骤S4,Re-execute step S3 to step S4,
Figure PCTCN2022125794-appb-000054
Figure PCTCN2022125794-appb-000054
Figure PCTCN2022125794-appb-000055
Figure PCTCN2022125794-appb-000055
Figure PCTCN2022125794-appb-000056
Figure PCTCN2022125794-appb-000056
T max=14.577°>T Agl T max =14.577°>T Agl
聚类不收敛,n=n+1=2,更新种子集合:Clustering does not converge, n=n+1=2, update the seed set:
Figure PCTCN2022125794-appb-000057
Figure PCTCN2022125794-appb-000057
步骤S10:再次重复步骤S3至步骤S6,Step S10: Repeat steps S3 to S6 again,
Figure PCTCN2022125794-appb-000058
Figure PCTCN2022125794-appb-000058
T max=3.885°<T Agl T max =3.885°<T Agl
聚类收敛,退出循环,输出M nThe clustering converges, the loop exits, and M n is output.
步骤S11:将测量法矢集合M n转化为测量摆角集合M AStep S11: Convert the measurement normal vector set M n into the measurement pendulum angle set M A .
Figure PCTCN2022125794-appb-000059
Figure PCTCN2022125794-appb-000059
取整后:After rounding:
Figure PCTCN2022125794-appb-000060
Figure PCTCN2022125794-appb-000060
该集合即为孔特征法矢集合
Figure PCTCN2022125794-appb-000061
对应的摆角集合。
This set is the hole feature normal vector set
Figure PCTCN2022125794-appb-000061
The corresponding set of pendulum angles.
对比例:Comparative ratio:
本对比例提供的是,根据公开号为CN112033331A的中国专利《一种基于三坐标测针的群孔测量摆角规划方法》所公开的技术方案,所作出的群孔测量摆角规划方法,具体步骤如下:What this comparative example provides is a group hole measurement swing angle planning method based on the technical solution disclosed in the Chinese patent "A Group Hole Measurement Swing Angle Planning Method Based on a Three-Coordinate Stylus" with the publication number CN112033331A. Specifically, Proceed as follows:
步骤S1:获取孔特征法矢集合
Figure PCTCN2022125794-appb-000062
Step S1: Obtain the hole feature normal vector set
Figure PCTCN2022125794-appb-000062
Figure PCTCN2022125794-appb-000063
Figure PCTCN2022125794-appb-000063
测量深度h=3mm、容差信息tol=0.5°、测杆直径d=1.5mm、测头直径D=2mm;Measuring depth h=3mm, tolerance information tol=0.5°, rod diameter d=1.5mm, probe diameter D=2mm;
步骤S12:计算测量不干涉条件下,设置初始聚类规模T=1,初始化种子集合
Figure PCTCN2022125794-appb-000064
Figure PCTCN2022125794-appb-000065
设置各聚类规模下的最大聚类迭代次数T n=5,当前迭代次数n=1,测针轴线与孔轴线的夹角阈值T Agl
Step S12: Under the condition of non-interference in calculation and measurement, set the initial clustering scale T = 1 and initialize the seed set.
Figure PCTCN2022125794-appb-000064
Figure PCTCN2022125794-appb-000065
Set the maximum number of clustering iterations under each clustering scale T n =5, the current number of iterations n =1, and the angle threshold T Agl between the stylus axis and the hole axis:
T Agl=3.9272°。 T Agl =3.9272°.
步骤S3:将机床摆角范围按2T Agl离散化,产生机床节点A1: Step S3: Discretize the machine tool swing angle range according to 2T Agl to generate machine tool node A1:
Figure PCTCN2022125794-appb-000066
Figure PCTCN2022125794-appb-000066
步骤S4:对于
Figure PCTCN2022125794-appb-000067
中每一个元素
Figure PCTCN2022125794-appb-000068
求解A1中与其夹角不大于
Figure PCTCN2022125794-appb-000069
的元素:
Step S4: For
Figure PCTCN2022125794-appb-000067
every element in
Figure PCTCN2022125794-appb-000068
Solve the problem that the angle between A1 and A1 is not greater than
Figure PCTCN2022125794-appb-000069
Elements:
在A1中,与(0.0000,0.0000,1.0000)夹角不大于3.9272°的有(0.32,1.24)In A1, the angle between (0.0000, 0.0000, 1.0000) is not greater than 3.9272° (0.32, 1.24)
与(0.0000,0.0712,0.9503)夹角不大于3.9272°的有(0.32,1.24)The angle between (0.0000, 0.0712, 0.9503) is not greater than 3.9272° (0.32, 1.24)
与(0.0000,0.0508,0.9802)夹角不大于3.9272°的有(0.32,1.24)The angle between (0.0000, 0.0508, 0.9802) is not greater than 3.9272° (0.32, 1.24)
没有与(0.3536,-0.6124,0.7071)夹角不大于3.9272°的元素。There is no element whose angle with (0.3536, -0.6124, 0.7071) is not greater than 3.9272°.
步骤S5:A11=(0.32,1.24)对应3个
Figure PCTCN2022125794-appb-000070
中元素,是A1所有节点对应元素最多的节点,
Figure PCTCN2022125794-appb-000071
Step S5: A11=(0.32, 1.24) corresponds to 3
Figure PCTCN2022125794-appb-000070
The middle element is the node with the most corresponding elements among all nodes in A1.
Figure PCTCN2022125794-appb-000071
步骤S6:Y1不为空,从
Figure PCTCN2022125794-appb-000072
中删除Y1;
Step S6: Y1 is not empty, from
Figure PCTCN2022125794-appb-000072
Delete Y1;
步骤S7:从
Figure PCTCN2022125794-appb-000073
中删除Y1,
Figure PCTCN2022125794-appb-000074
Figure PCTCN2022125794-appb-000075
不为空,执行S4;
Step S7: From
Figure PCTCN2022125794-appb-000073
Delete Y1 in
Figure PCTCN2022125794-appb-000074
Figure PCTCN2022125794-appb-000075
If not empty, execute S4;
步骤S8;在A1中,没有与(0.3536,-0.6124,0.7071)夹角不大于3.9272°的元素。Step S8; In A1, there is no element whose angle with (0.3536, -0.6124, 0.7071) is not greater than 3.9272°.
步骤S9:A1任意元素均对应0个
Figure PCTCN2022125794-appb-000076
中元素,Y1为空集;
Step S9: Any element of A1 corresponds to 0
Figure PCTCN2022125794-appb-000076
Medium element, Y1 is the empty set;
步骤S10:剩余节点已经不能测量剩余的孔,则需要减小间隔再次划分产生更密集的节点;Step S10: The remaining nodes cannot measure the remaining holes, so the intervals need to be reduced and divided again to generate denser nodes;
步骤S11:执行步骤S3,将机床摆角范围按T Agl重新离散化,产生机床节点A2: Step S11: Execute step S3 to re-discretize the machine tool swing angle range according to T Agl to generate machine tool node A2:
Figure PCTCN2022125794-appb-000077
Figure PCTCN2022125794-appb-000077
步骤S12:对于
Figure PCTCN2022125794-appb-000078
中每一个元素
Figure PCTCN2022125794-appb-000079
求解A2中与其夹角不大于T Agl的元素:
Step S12: For
Figure PCTCN2022125794-appb-000078
every element in
Figure PCTCN2022125794-appb-000079
Solve for the elements in A2 whose angle is no greater than T Agl :
在A1中,与(0.0000,0.0000,1.0000)夹角不大于3.9272°的有(-0.25,1.24)(3.33,1.24),(-0.25,5.34),(3.33,5.34);In A1, the angles with (0.0000, 0.0000, 1.0000) are not greater than 3.9272°: (-0.25, 1.24) (3.33, 1.24), (-0.25, 5.34), (3.33, 5.34);
与(0.0000,0.0712,0.9503)夹角不大于3.9272°的有(-0.25,1.24),(3.33,1.24),(-0.25,5.34),(-3.43,5.34);The angles with (0.0000, 0.0712, 0.9503) are not greater than 3.9272°: (-0.25, 1.24), (3.33, 1.24), (-0.25, 5.34), (-3.43, 5.34);
与(0.0000,0.0508,0.9802)夹角不大于3.9272°的有(-0.25,1.24),(3.33,1.24),(-0.25,5.34),(-3.43,5.34);The angles with (0.0000, 0.0508, 0.9802) are not greater than 3.9272°: (-0.25, 1.24), (3.33, 1.24), (-0.25, 5.34), (-3.43, 5.34);
与(0.3536,-0.6124,0.7071)夹角不大于3.9272°的有(45.23,28.22),(45.23,31.25);The angles with (0.3536, -0.6124, 0.7071) are not greater than 3.9272°: (45.23, 28.22), (45.23, 31.25);
S13:A21=(-0.25,1.24)对应3个
Figure PCTCN2022125794-appb-000080
中元素,是A2所有节点对应元素最多的节点,
Figure PCTCN2022125794-appb-000081
S13: A21=(-0.25, 1.24) corresponds to 3
Figure PCTCN2022125794-appb-000080
The middle element is the node with the most corresponding elements among all nodes in A2.
Figure PCTCN2022125794-appb-000081
步骤S14:Y2不为空,从
Figure PCTCN2022125794-appb-000082
中删除Y2;
Step S14: Y2 is not empty, from
Figure PCTCN2022125794-appb-000082
Delete Y2;
步骤S15:从
Figure PCTCN2022125794-appb-000083
中删除Y2,
Figure PCTCN2022125794-appb-000084
Figure PCTCN2022125794-appb-000085
不为空,执行S12
Step S15: From
Figure PCTCN2022125794-appb-000083
Delete Y2 in
Figure PCTCN2022125794-appb-000084
Figure PCTCN2022125794-appb-000085
If not empty, execute S12
步骤S16:在A2中,与(0.3536,-0.6124,0.7071)夹角不大于3.9272°的有(45.23,28.22),(45.23,31.25);Step S16: In A2, the angles with (0.3536, -0.6124, 0.7071) are not greater than 3.9272° including (45.23, 28.22), (45.23, 31.25);
步骤S17:A22=(45.23,28.22)对应1个
Figure PCTCN2022125794-appb-000086
中元素,是A2所有节点对应元素最多的节点,
Figure PCTCN2022125794-appb-000087
Step S17: A22=(45.23, 28.22) corresponds to 1
Figure PCTCN2022125794-appb-000086
The middle element is the node with the most corresponding elements among all nodes in A2.
Figure PCTCN2022125794-appb-000087
步骤S18:Y2不为空,从
Figure PCTCN2022125794-appb-000088
中删除Y2;
Step S18: Y2 is not empty, from
Figure PCTCN2022125794-appb-000088
Delete Y2;
步骤S19:从
Figure PCTCN2022125794-appb-000089
中删除Y2后,
Figure PCTCN2022125794-appb-000090
为空,测量摆角规划完成。
Step S19: From
Figure PCTCN2022125794-appb-000089
After deleting Y2 in
Figure PCTCN2022125794-appb-000090
If it is empty, the measurement and swing angle planning is completed.
步骤S20:最终规划的测量摆角为Step S20: The final planned measured swing angle is
Figure PCTCN2022125794-appb-000091
Figure PCTCN2022125794-appb-000091
本实施例中提供的群孔测量摆角规划方法与对比例中提供的摆角规划方法相比,本实施例通过摆角规划,测量角度准确性大幅提高,所需测量摆角数量降低,孔测量摆角规划时间大幅度降低,能够快速获得孔特征测量摆角,节约了测量准备时间,提高了加工及测量效率。Compared with the swing angle planning method for group hole measurement provided in this embodiment and the swing angle planning method provided in the comparative example, this embodiment greatly improves the accuracy of measurement angles through swing angle planning, reduces the number of required measured swing angles, and reduces the number of holes. The planning time of measuring the swing angle is greatly reduced, and the hole characteristic measurement swing angle can be quickly obtained, which saves the measurement preparation time and improves the processing and measurement efficiency.
以上所述,仅是本发明的较佳实施例,并非对本发明做任何形式上的限制,凡是依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化,均落入本发明的保护范围之内。The above are only preferred embodiments of the present invention and do not limit the present invention in any form. Any simple modifications or equivalent changes made to the above embodiments based on the technical essence of the present invention fall within the scope of the present invention. within the scope of protection.

Claims (10)

  1. 群孔测量摆角规划方法,用于在机检测孔特征,其特征在于,包括以下步骤:The group hole measurement swing angle planning method is used to detect hole characteristics on the machine, which is characterized by including the following steps:
    步骤S1:获取孔特征法矢集合、测针参数、测量深度和容差信息;Step S1: Obtain the hole feature normal vector set, stylus parameters, measurement depth and tolerance information;
    步骤S2:根据测针参数、测量深度和容差信息,计算出孔与测针不发生干涉时测针轴线与孔轴线之间的夹角阈值;设置聚类规模的初始值、迭代次数的初始值、初始化种子集合、并确定最大迭代次数;Step S2: Based on the stylus parameters, measurement depth and tolerance information, calculate the angle threshold between the stylus axis and the hole axis when the hole and stylus do not interfere; set the initial value of the clustering scale and the initial number of iterations. value, initialize the seed collection, and determine the maximum number of iterations;
    步骤S3:执行聚类:以孔特征法矢集合为数据集,计算聚类集合和测量法矢集合,并计算候选种子;Step S3: Perform clustering: use the hole feature normal vector set as the data set, calculate the clustering set and the measurement normal vector set, and calculate candidate seeds;
    步骤S4:判断聚类收敛性:若聚类收敛,则执行步骤S6;若聚类不收敛,则执行步骤S5;Step S4: Determine the clustering convergence: if the clustering converges, execute step S6; if the clustering does not converge, execute step S5;
    步骤S5:先更新种子集合;然后,若迭代次数小于等于最大迭代次数,则迭代次数增加一次,并重新依次执行步骤S3、步骤S4;若迭代次数大于最大迭代次数,则聚类规模加一,增加候选种子至当前种子集合形成新的种子集合,当前迭代次数回归初始值,并重新执行步骤S3、步骤S4;Step S5: Update the seed set first; then, if the number of iterations is less than or equal to the maximum number of iterations, increase the number of iterations by one, and re-execute steps S3 and S4 in sequence; if the number of iterations is greater than the maximum number of iterations, increase the clustering size by one, Add candidate seeds to the current seed set to form a new seed set, return the current iteration number to the initial value, and re-execute steps S3 and S4;
    步骤S6:将测量法矢集合转化为测量摆角集合,输出测量摆角集合并结束程序。Step S6: Convert the set of measured normal vectors into a set of measured pendulum angles, output the set of measured pendulum angles and merge to end the program.
  2. 根据权利要求1所述的群孔测量摆角规划方法,其特征在于:所述测针参数包括测头直径和测杆直径。The method for planning the swing angle of group hole measurement according to claim 1, wherein the stylus parameters include a probe diameter and a probe rod diameter.
  3. 根据权利要求1所述的群孔测量摆角规划方法,其特征在于:所述容差信息的取值大于等于0.5°。The group hole measurement swing angle planning method according to claim 1, characterized in that: the value of the tolerance information is greater than or equal to 0.5°.
  4. 根据权利要求1所述的群孔测量摆角规划方法,其特征在于:聚类规模的初始值为1、迭代次数的初始值为1;初始化种子集合为:将孔特征法矢中第一个元素作为种子集合的唯一种子。The group hole measurement swing angle planning method according to claim 1, characterized in that: the initial value of the clustering scale is 1, the initial value of the number of iterations is 1; the initialization seed set is: the hole feature normal vector is the first Element serves as the unique seed of the seed collection.
  5. 根据权利要求1所述的群孔测量摆角规划方法,其特征在于:在步骤S3中,通过以下具体步骤计算得到聚类集合和测量法矢集合,具体步骤如下:The group hole measurement swing angle planning method according to claim 1, characterized in that: in step S3, the cluster set and the measurement normal vector set are calculated through the following specific steps, and the specific steps are as follows:
    步骤S31:以孔特征法矢集合中第i元素与种子集合中各种子分别形成的夹角,加入第i夹角集合中作为元素;Step S31: Add the angles formed by the i-th element in the hole feature normal vector set and various sub-elements in the seed set as elements into the i-th included angle set;
    步骤S32:若种子集合中第j种子,与孔特征法矢集合中第i元素形成夹角,为第i夹角集合取值最小的元素;则将孔特征法矢集合中第i元素加入聚类集合中第j聚类子集,对应将种子集合中第j种子加入测量法矢集合;其中,{i∈Z|1≤i≤N},N为孔特征法矢集合中元素的数量;Step S32: If the j-th seed in the seed set forms an angle with the i-th element in the hole feature normal vector set, and is the element with the smallest value in the i-th angle set; then add the i-th element in the hole feature normal vector set to the cluster. The jth cluster subset in the class set corresponds to adding the jth seed in the seed set to the measurement normal vector set; where, {i∈Z|1≤i≤N}, N is the number of elements in the hole feature normal vector set;
    步骤S33:遍历聚类集合,将聚类集合中第j聚类子集的各元素与种子集合中第j种子形成的夹角加入候选集合作为元素,候选集合中元素取得最大值时对应的聚类集合中第j个聚类 子集的特定元素作为候选种子;其中,{j∈Z|1≤j≤M},M为种子集合中种子的数量。Step S33: Traverse the cluster set, and add the angle formed by each element of the j-th cluster subset in the cluster set and the j-th seed in the seed set to the candidate set as an element. The corresponding cluster when the element in the candidate set reaches the maximum value is The specific elements of the jth cluster subset in the class set are used as candidate seeds; where {j∈Z|1≤j≤M}, M is the number of seeds in the seed set.
  6. 根据权利要求5所述的群孔测量摆角规划方法,其特征在于:所述聚类集合中聚类子集与种子集合中种子顺次一一对应。The group hole measurement swing angle planning method according to claim 5, characterized in that: the cluster subsets in the cluster set correspond to the seeds in the seed set in sequence.
  7. 根据权利要求5所述的群孔测量摆角规划方法,其特征在于,当所述候选集合中元素取得的最大值大于夹角阈值,则判定聚类不收敛,否则判定聚类收敛。The group hole measurement swing angle planning method according to claim 5, characterized in that when the maximum value obtained by the elements in the candidate set is greater than the angle threshold, it is determined that the clustering does not converge, otherwise it is determined that the clustering has converged.
  8. 根据权利要求1所述的群孔测量摆角规划方法,其特征在于:步骤S5中更新种子集合的方法为:将聚类集合中各聚类子集一一转换为新的种子集合中各种子,具体为使用各聚类子集中元素平均值依次作为新的种子集合中的各种子。The group hole measurement swing angle planning method according to claim 1, characterized in that: the method for updating the seed set in step S5 is to convert each cluster subset in the cluster set one by one into various types in the new seed set. Specifically, the average value of the elements in each cluster subset is used as the various children in the new seed set.
  9. 一种计算机可读介质,其上存储有计算机程序,其特征在于:所述计算机程序被执行时实现上述权利要求1至8中任一项所述的群孔测量摆角规划方法。A computer-readable medium with a computer program stored thereon, characterized in that when the computer program is executed, the method for planning the group hole measurement swing angle described in any one of claims 1 to 8 is implemented.
  10. 一种电子设备,其特征在于,包括:An electronic device, characterized by including:
    处理器;processor;
    存储器,用于存储所述处理器的可执行指令;memory for storing executable instructions for the processor;
    其中,所述处理器配置为经由执行所述可执行指令来执行权利要求1至8中任意一项所述的群孔测量摆角规划方法。Wherein, the processor is configured to execute the group hole measurement swing angle planning method according to any one of claims 1 to 8 by executing the executable instructions.
PCT/CN2022/125794 2022-07-15 2022-10-18 Swing angle planning method for group-hole measurement, and readable medium and device WO2024011779A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202210830966.2A CN114897110B (en) 2022-07-15 2022-07-15 Group hole measurement swing angle planning method, readable medium and equipment
CN202210830966.2 2022-07-15

Publications (1)

Publication Number Publication Date
WO2024011779A1 true WO2024011779A1 (en) 2024-01-18

Family

ID=82729327

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2022/125794 WO2024011779A1 (en) 2022-07-15 2022-10-18 Swing angle planning method for group-hole measurement, and readable medium and device

Country Status (2)

Country Link
CN (1) CN114897110B (en)
WO (1) WO2024011779A1 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114897110B (en) * 2022-07-15 2022-11-18 成都飞机工业(集团)有限责任公司 Group hole measurement swing angle planning method, readable medium and equipment
CN115577279B (en) * 2022-11-11 2023-03-10 季华实验室 Punching clustering method and device, and punching process design and processing method of punching part

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102176211A (en) * 2010-12-26 2011-09-07 北京航空航天大学 Path planning method of measuring system based on projected grating phase method
EP3435162A1 (en) * 2017-07-28 2019-01-30 ASML Netherlands B.V. Metrology method and apparatus and computer program
CN111401413A (en) * 2020-02-29 2020-07-10 同济大学 Optimization theory-based parallel clustering method with scale constraint
CN112033331A (en) * 2020-07-28 2020-12-04 成都飞机工业(集团)有限责任公司 Group hole measurement swing angle planning method based on three-coordinate measuring needle
CN114193231A (en) * 2022-02-16 2022-03-18 成都飞机工业(集团)有限责任公司 Bottom hole orifice measuring method for numerical control countersink
CN114227026A (en) * 2021-12-10 2022-03-25 西安交通大学 Ultrafast laser hole-controllable group hole precision machining device and method
CN114897110A (en) * 2022-07-15 2022-08-12 成都飞机工业(集团)有限责任公司 Group hole measurement swing angle planning method, readable medium and equipment

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103900510A (en) * 2014-03-04 2014-07-02 成都飞机工业(集团)有限责任公司 Detection path modeling method of coordinate measuring machine
CN107554813B (en) * 2017-09-05 2019-06-14 西北工业大学 A kind of drilling interference plug pin integrated apparatus and method
CN111664813B (en) * 2020-05-15 2021-08-03 成都飞机工业(集团)有限责任公司 Free surface arbitrary hole normal vector measuring device and method and compensation method
CN111859026B (en) * 2020-06-29 2022-08-09 北京航空航天大学 XML-based three-coordinate measurement and evaluation information transfer method and system
CN111815611B (en) * 2020-07-14 2024-03-22 南京航空航天大学 Round hole feature extraction method for rivet hole measurement point cloud data
CN112496863B (en) * 2020-11-23 2022-10-28 中国航发沈阳黎明航空发动机有限责任公司 Automatic measuring method for numerical control machining angle
CN113239580A (en) * 2020-12-21 2021-08-10 南京航空航天大学 Laser radar measuring station position planning method for large structural member profile detection
CN114186645B (en) * 2022-02-15 2022-04-26 中国长江三峡集团有限公司 Blasting half-porosity detection method and device, electronic equipment and storage medium
CN114723160B (en) * 2022-04-22 2023-08-04 成都飞机工业(集团)有限责任公司 Measurement path planning method for on-machine detection

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102176211A (en) * 2010-12-26 2011-09-07 北京航空航天大学 Path planning method of measuring system based on projected grating phase method
EP3435162A1 (en) * 2017-07-28 2019-01-30 ASML Netherlands B.V. Metrology method and apparatus and computer program
CN111401413A (en) * 2020-02-29 2020-07-10 同济大学 Optimization theory-based parallel clustering method with scale constraint
CN112033331A (en) * 2020-07-28 2020-12-04 成都飞机工业(集团)有限责任公司 Group hole measurement swing angle planning method based on three-coordinate measuring needle
CN114227026A (en) * 2021-12-10 2022-03-25 西安交通大学 Ultrafast laser hole-controllable group hole precision machining device and method
CN114193231A (en) * 2022-02-16 2022-03-18 成都飞机工业(集团)有限责任公司 Bottom hole orifice measuring method for numerical control countersink
CN114897110A (en) * 2022-07-15 2022-08-12 成都飞机工业(集团)有限责任公司 Group hole measurement swing angle planning method, readable medium and equipment

Also Published As

Publication number Publication date
CN114897110B (en) 2022-11-18
CN114897110A (en) 2022-08-12

Similar Documents

Publication Publication Date Title
WO2024011779A1 (en) Swing angle planning method for group-hole measurement, and readable medium and device
CN107332240B (en) Method for searching static voltage stability domain boundary of power system based on optimization model
CN111784778B (en) Binocular camera external parameter calibration method and system based on linear solving and nonlinear optimization
CN109493389B (en) Camera calibration method and system based on deep learning
CN110398219B (en) Joint arm type coordinate measuring machine parameter calibration method based on hybrid optimization algorithm
CN113211445B (en) Robot parameter calibration method, device, equipment and storage medium
WO2017107865A1 (en) Image retrieval system, server, database, and related method
CN111783343A (en) Method for selecting direction template of unstructured grid and flow field solving method
WO2022193640A1 (en) Robot calibration method and apparatus, and robot and storage medium
WO2019114631A1 (en) Method and device for acquiring tcp coordinates of robot
CN111723441B (en) Plastic part pre-deformation part modeling method, system, device and storage medium
CN111126268B (en) Key point detection model training method and device, electronic equipment and storage medium
CN111487440A (en) Calibration method of five-hole probe
CN115512343A (en) Method for correcting and recognizing reading of circular pointer instrument
CN114374931B (en) Fingerprint positioning method based on metric learning of neighbor component analysis
CN110222382B (en) Aircraft coordinate system optimization fitting method
CN116402144A (en) Determination method and determination device for multi-quantum bit measurement result and quantum computer
CN109631813B (en) Calibration method of large-size articulated arm type coordinate measuring machine
CN116136388A (en) Calibration method, device, equipment and storage medium for robot tool coordinate system
CN112184561A (en) Rocket stage pose deviation rectifying method and device based on multivariate search
CN112182739A (en) Aircraft structure non-probability credible reliability topological optimization design method
CN115847392B (en) Force sensor no-load reading compensation method, device and storage medium
CN111637879B (en) Double-star-sensitive weighted attitude determination method based on multi-dimensional differential error characteristic distribution
CN112596382B (en) Geometric parameter optimization calibration method and system for series servo mechanism
TWI793044B (en) Eye-hand calibration method and eye-hand calibration device for robot arm

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 22950889

Country of ref document: EP

Kind code of ref document: A1