一种面向大部件群孔加工设备精度指标优化方法与流程

1.本发明涉及机械加工领域，具体地说，涉及一种面向大部件群孔加工设备精度指标优化方法。


背景技术：

2.在飞机装配过程中，结构件之间的装配主要以机械连接为主，连接孔的加工质量对于提高连接精度、提高飞机性能及使用寿命具有非常重要的影响。以人工为主的装配模式存在效率低、稳定性差、工艺知识继承难等问题，难以在有限资源支撑下保证优质、高效、低成本批量生产。结合现代飞机数字化制造的发展趋势，飞机装配技术从手工劳动、半机械化、机械化逐步向数字化、柔性化、自动化的方向发展。因此多轴联动数控机床越来越广泛的应用于飞机装配制造领域。
3.在实际制孔应用中，由于对加工质量需求，对多轴联动数控机床精度的要求往往也比较高。而在机床设计、招标过程中，机床的精度关系到机床的制造难度以及制造成本，盲目制定过高的技术指标会带来技术成本、经济成本的增加。因此如何确定机床各轴的精度指标使得精度既能满足加工要求又能相对宽松是一项重要的待解决的问题。
4.多轴联动数控机床的精度指标包含各轴定位精度、重复定位精度、直线度、垂直度、角度偏差等精度指标，总计多达几十余种。在传统工作模式下，这些精度指标确定方式依靠人工经验、加工反馈进行确定，存在较大的不确定性、制定难度以及较长的制定周期，制定过程缺乏量化的、科学化的依据。在实际应用中，由于技术条件、实施条件等实际情况限制，必须对某一些精度指标进行固定，对剩余的精度指标进行优化。在该种情况下，依靠人工经验来确定精度指标实现难度高，且缺乏科学依据。因此亟需一种科学化的机床精度指标优化方法，用以确定一组定量的精度指标，使其满足既能加工精度要求，又能相对宽松，降低机床的制造难度、制造成本。


技术实现要素：

5.本发明针对现有技术中依靠人工经验来确定精度指标实现难度高，且缺乏科学依据的问题，提出一种面向大部件群孔加工设备精度指标优化方法，通过从数控机床群孔加工孔位误差模型中提取了不同种类的精度指标，并根据这些指标，建立了数控机床精度指标评价模型、精度指标适应度模型，并根据建立的模型，以精度指标为粒子，通过计算出精度指标粒子的全局最优点gbest
t
和精度指标粒子的局部最优点pbest
kt
，结合机床加工过程中群孔孔位位置误差，计算出一组定量的精度指标，使得在降低机床制造难度、制造成本的情况下，又能保证机床满足精度需求。
6.本发明具体实现内容如下：一种面向大部件群孔加工设备精度指标优化方法，包括以下步骤：步骤1：根据数控机床类型，建立数控机床的拓扑结构；步骤2：根据步骤1建立的数控机床的拓扑结构，建立数控机床刀尖点运动的空间
定位误差模型p
error
；步骤3：根据步骤2建立的空间定位误差模型p
error
，以及制孔孔位误差与理想轴线垂直平面的投影关系，建立数控机床群孔加工孔位误差模型；步骤4：根据步骤3建立的数控机床群孔加工孔位误差模型，提取数控机床群孔加工孔位误差模型中的机床的精度指标；步骤5：计算出步骤4提取出的精度指标的合集，建立数控机床精度指标评价模型；步骤6：根据步骤5建立的精度指标评价模型，判断是否满足设备制孔精度要求p，并根据判断结果建立精度指标适应度模型；步骤7：根据步骤6建立的精度指标适应度模型，计算出优化后的精度指标粒子速度v
kt+1
和优化后的精度指标粒子位置g
kt+1
；步骤8：以精度指标为粒子，设置最大适应度和最大迭代次数k，初始化粒子群粒子，计算出精度指标粒子的初始位置g
k1
和精度指标粒子的初始速度v
k1
；步骤9：根据步骤7计算出的优化后的精度指标粒子速度v
kt+1
和优化后的精度指标粒子位置g
kt+1
和步骤8计算出的精度指标粒子的初始位置g
k1
和精度指标粒子的初始速度v
k1
，计算出精度指标粒子的局部最优点pbest
kt
、精度指标粒子的全局最优点gbest
t
；步骤10：根据步骤9计算出的全局最优点gbest
t
和局部最优点pbest
kt
以及上一轮迭代的精度指标粒子速度，更新精度指标粒子速度和精度指标粒子位置，得到更新后的精度指标粒子速度v’kt+1
、更新后的精度指标粒子位置x
kt+1
；步骤11：判断是否满足结束条件，满足则当前精度指标为最优，不满足则重复步骤9，直至满足条件。
7.为了更好地实现本发明，进一步地，所述步骤2具体包括以下步骤：步骤21：根据步骤1建立的机床拓扑结构，建立相邻运动单元之间的运动传递矩阵m
ij
；步骤22：根据步骤21建立的相邻运动单元之间的运动传递矩阵m
ij
，结合刀尖点的初始定位p
initial
，计算出机床在进给轴输入参数作用下机床刀尖的空间定位向量p
ideal
；步骤23：根据步骤21建立的相邻运动单元之间的运动传递矩阵m
ij
，结合刀具在a摆动轴坐标系下的初始姿态v
initial
，计算出机床在进给轴输入参数作用下机床刀尖的空间姿态向量v
ideal
；步骤24：根据机床进给轴的误差矩阵，计算出线性进给轴u的误差传递矩阵e
uij
和转动进给轴v的误差传递矩阵e
vij
；步骤25：根据步骤24计算出的线性进给轴u的误差传递矩阵e
uij
、刀尖点的初始定位p
initial
、运动传递矩阵m
ij
，计算出实际空间定位向量p
actual
；步骤26：根据步骤22计算出的空间定位向量p
ideal
、步骤25计算出的实际空间定位向量p
actual
，建立数控机床刀尖点运动的空间定位误差模型p
error
。
8.为了更好地实现本发明，进一步地，所述步骤3具体包括以下步骤：步骤31：根据制孔孔位误差在与理想轴线垂直平面的投影关系，步骤23计算出的空间姿态向量v
ideal
、步骤26建立的空间定位误差模型p
error
，计算出制孔位置处位置偏差p
error_v
；步骤32：根据步骤31计算出的制孔位置处位置偏差p
error_v
，计算出当前机床制孔
时孔位误差re，建立数控机床群孔加工孔位误差模型。
9.为了更好地实现本发明，进一步地，所述步骤4的具体操作为：将线性进给轴u的误差传递矩阵e
uij
和转动进给轴v的误差传递矩阵e
vij
中包含的定位精度指标组成集合g1、重复定位精度指标组成集合g2、直线度定位精度指标组成集合g3、角度偏差指标组成集合g4、垂直度指标组成集合g5。
10.为了更好地实现本发明，进一步地，所述步骤5具体包括以下步骤：步骤51：计算出定位精度指标组成集合g1、重复定位精度指标组成集合g2、直线度定位精度指标组成集合g3、角度偏差指标组成集合g4、垂直度指标组成集合g5的合集，得到当前设备的精度指标集合q；步骤52：根据构成当前设备的精度指标集合q的精度指标向量g、制孔过程中所制每个孔各轴进给量pi，计算出每个制孔点位的孔位误差re(g，pi)；步骤53：根据步骤52计算出的每个制孔点位的孔位误差re(g，pi)，计算出所有加工孔位的平均误差r
ea
；步骤54：判断平均误差r
ea
是否满足设备制孔精度要求p，若满足，计算出单组精度指标评价系数ji(gi)；步骤55：根据步骤54计算出的单组精度指标评价系数ji(gi)，计算出总体精度指标评价系数j(g)，建立数控机床精度指标评价模型。
11.为了更好地实现本发明，进一步地，所述步骤6的具体操作为：根据步骤5建立的精度指标评价模型，判断所有加工孔位的平均误差r
ea
是否满足设备制孔精度要求p，若所有加工孔位的平均误差r
ea
满足精度要求p，则适应度为总体精度指标评价系数j(g)，若所有加工孔位的平均误差r
ea
不满足精度要求p，则适应度为负的所有加工孔位的平均误差r
ea
，根据判断结果建立精度指标适应度模型。
12.为了更好地实现本发明，进一步地，所述步骤7具体包括以下步骤：步骤71：根据步骤6建立的建立精度指标适应度模型、第t次迭代时设备的精度指标g
kt
、每个精度指标粒子的速度v
kt
，计算出精度指标粒子的局部最优点pbest
kt
、精度指标粒子的全局最优点gbest
t
；步骤72：根据步骤71计算出的精度指标粒子的局部最优点pbest
kt
、精度指标粒子的全局最优点gbest
t
、惯性因子、学习因子c1、学习因子c2，计算出优化后的精度指标粒子的速度v
kt+1
；步骤73：根据优化后的精度指标粒子的速度v
kt+1
，计算出优化后的精度指标粒子的位置g
kt+1
。
13.为了更好地实现本发明，进一步地，所述步骤8具体包括以下步骤：步骤81：以精度指标为粒子，初始化粒子群粒子，根据设定的初始精度指标、最大适应度和最大迭代次数k、最大速度系数，计算出相对最大速度v
max
；步骤82：根据设定的初始精度指标，计算出初始局部最优点pbest
k0
、初始全局最优点gbest0、精度指标粒子的初始位置g
k1
、精度指标粒子的初始速度v
k1
。
14.为了更好地实现本发明，进一步地，所述步骤11的具体操作为：判断是否满足结束条件，若全局最优点gbest
t
的适应度大于最大适应度或迭代次数t大于最大迭代次数k，则满足结束条件，将当前精度指标作为最优精度指标，结束迭代，若不满足结束条件，则将
迭代次数加一，重复步骤9，直至满足结束条件。
15.为了更好地实现本发明，进一步地，所述步骤1的具体操作为：根据低序体阵列法对机床的拓扑结构进行简化，并对相对应运动单元进行低序体编号。
16.本发明具有以下有益效果：（1）本发明通过建立数控机床精度指标优化评价模型，为后续利用粒子群算法得到满足要求的设备精度指标提供保障，使得各组精度指标之间具有差异性又能保证每组精度指标内部相对接近；（2）本发明通过建立精度指标的适应度模型，实现了在精度限制条件下评价精度指标好坏；（3）本发明针对机床不同精度指标数量级存在差异的特点，使用相对最大速度概念很好的使各组指标速度分离，互不影响，大大提高了算法优化成功率；（4）本发明提出的面对大部件群孔加工设备精度指标优化方法，可在工艺设计阶段对多轴数控机床的极限精度指标进行评估，依据该指标进行机床的设计、制造及采购；用以确定一组定量的精度指标，使其满足既能加工精度要求，又能相对宽松，降低机床的制造难度、制造成本。
附图说明
17.图1为本发明的面向大部件群孔加工设备精度指标优化方法技术路线图。
18.图2为本发明的制孔孔位误差示意图。
19.图3为本发明的精度指标优化粒子速度更新示意图。
20.图4为本发明的精度指标优化粒子位置更新示意图。
21.图5为本发明的精度指标优化粒子群优化过程示意图。
22.图6为本发明的精度指标优化粒子群优化结果示意图。
具体实施方式
23.为了更清楚地说明本发明实施例的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，应当理解，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例，因此不应被看作是对保护范围的限定。基于本发明中的实施例，本领域普通技术工作人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
24.在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“设置”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；也可以是直接相连，也可以是通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
25.实施例1：本实施例提出一种面向大部件群孔加工设备精度指标优化方法，包括以下步骤：步骤1：根据数控机床类型，建立数控机床的拓扑结构；步骤2：根据步骤1建立的数控机床的拓扑结构，建立数控机床刀尖点运动的空间
定位误差模型p
error
；步骤3：根据步骤2建立的空间定位误差模型p
error
，以及制孔孔位误差与理想轴线垂直平面的投影关系，建立数控机床群孔加工孔位误差模型；步骤4：根据步骤3建立的数控机床群孔加工孔位误差模型，提取数控机床群孔加工孔位误差模型中的机床的精度指标；步骤5：计算出步骤4提取出的精度指标的合集，建立数控机床精度指标评价模型；步骤6：根据步骤5建立的精度指标评价模型，判断是否满足设备制孔精度要求p，并根据判断结果建立精度指标适应度模型；步骤7：根据步骤6建立的精度指标适应度模型，计算出优化后的精度指标粒子速度v
kt+1
和优化后的精度指标粒子位置g
kt+1
；步骤8：以精度指标为粒子，设置最大适应度和最大迭代次数k，初始化粒子群粒子，计算出精度指标粒子的初始位置g
k1
和精度指标粒子的初始速度v
k1
；步骤9：根据步骤7计算出的优化后的精度指标粒子速度v
kt+1
和优化后的精度指标粒子位置g
kt+1
和步骤8计算出的精度指标粒子的初始位置g
k1
和精度指标粒子的初始速度v
k1
，计算出精度指标粒子的局部最优点pbest
kt
、精度指标粒子的全局最优点gbest
t
；步骤10：根据步骤9计算出的全局最优点gbest
t
和局部最优点pbest
kt
以及上一轮迭代的精度指标粒子速度，更新精度指标粒子速度和精度指标粒子位置，得到更新后的精度指标粒子速度v’kt+1
、更新后的精度指标粒子位置x
kt+1
；步骤11：判断是否满足结束条件，满足则当前精度指标为最优，不满足则重复步骤9，直至满足条件。
26.工作原理：本实施例通过从数控机床群孔加工孔位误差模型中提取了不同种类的精度指标，并根据这些指标，建立了数控机床精度指标评价模型、精度指标适应度模型，并根据建立的模型，以精度指标为粒子，通过计算出精度指标粒子的全局最优点gbest
t
和精度指标粒子的局部最优点pbest
kt
，结合机床加工过程中群孔孔位位置误差，计算出一组定量的精度指标，使得在降低机床制造难度、制造成本的情况下，又能保证机床满足精度需求。
27.实施例2：本实施例在上述实施例1的基础上，如图2、图3、图4、图5、图6所示，对具体步骤进行说明。
28.所述步骤1的具体操作为：根据低序体阵列法对机床的拓扑结构进行简化，并对相对应运动单元进行低序体编号。
29.所述步骤2具体包括以下步骤：步骤21：根据步骤1建立的机床拓扑结构，建立相邻运动单元之间的运动传递矩阵m
ij
；步骤22：根据步骤21建立的相邻运动单元之间的运动传递矩阵m
ij
，结合刀尖点的初始定位p
initial
，计算出机床在进给轴输入参数作用下机床刀尖的空间定位向量p
ideal
；步骤23：根据步骤21建立的相邻运动单元之间的运动传递矩阵m
ij
，结合刀具在a摆动轴坐标系下的初始姿态v
initial
，计算出机床在进给轴输入参数作用下机床刀尖的空间姿态向量v
ideal
；
步骤24：根据机床进给轴的误差矩阵，计算出线性进给轴u的误差传递矩阵e
uij
和转动进给轴v的误差传递矩阵e
vij
；步骤25：根据步骤24计算出的线性进给轴u的误差传递矩阵e
uij
、刀尖点的初始定位p
initial
、运动传递矩阵m
ij
，计算出实际空间定位向量p
actual
；步骤26：根据步骤22计算出的空间定位向量p
ideal
、步骤25计算出的实际空间定位向量p
actual
，建立数控机床刀尖点运动的空间定位误差模型p
error
。
30.所述步骤3具体包括以下步骤：步骤31：根据制孔孔位误差在与理想轴线垂直平面的投影关系，步骤23计算出的空间姿态向量v
ideal
、步骤26建立的空间定位误差模型p
error
，计算出制孔位置处位置偏差p
error_v
；步骤32：根据步骤31计算出的制孔位置处位置偏差p
error_v
，计算出当前机床制孔时孔位误差re，建立数控机床群孔加工孔位误差模型。
31.所述步骤4的具体操作为：将线性进给轴u的误差传递矩阵e
uij
和转动进给轴v的误差传递矩阵e
vij
中包含的定位精度指标组成集合g1、重复定位精度指标组成集合g2、直线度定位精度指标组成集合g3、角度偏差指标组成集合g4、垂直度指标组成集合g5。
32.所述步骤5具体包括以下步骤：步骤51：计算出定位精度指标组成集合g1、重复定位精度指标组成集合g2、直线度定位精度指标组成集合g3、角度偏差指标组成集合g4、垂直度指标组成集合g5的合集，得到当前设备的精度指标集合q；步骤52：根据构成当前设备的精度指标集合q的精度指标向量g、制孔过程中所制每个孔各轴进给量pi，计算出每个制孔点位的孔位误差re(g，pi)；步骤53：根据步骤52计算出的每个制孔点位的孔位误差re(g，pi)，计算出所有加工孔位的平均误差r
ea
；步骤54：判断平均误差r
ea
是否满足设备制孔精度要求p，若满足，计算出单组精度指标评价系数ji(gi)；步骤55：根据步骤54计算出的单组精度指标评价系数ji(gi)，计算出总体精度指标评价系数j(g)，建立数控机床精度指标评价模型。
33.所述步骤6的具体操作为：根据步骤5建立的精度指标评价模型，判断所有加工孔位的平均误差r
ea
是否满足设备制孔精度要求p，若所有加工孔位的平均误差r
ea
满足精度要求p，则适应度为总体精度指标评价系数j(g)，若所有加工孔位的平均误差r
ea
不满足精度要求p，则适应度为负的所有加工孔位的平均误差r
ea
，根据判断结果建立精度指标适应度模型。
34.所述步骤7具体包括以下步骤：步骤71：根据步骤6建立的建立精度指标适应度模型、第t次迭代时设备的精度指标g
kt
、每个精度指标粒子的速度v
kt
，计算出精度指标粒子的局部最优点pbest
kt
、精度指标粒子的全局最优点gbest
t
；步骤72：根据步骤71计算出的精度指标粒子的局部最优点pbest
kt
、精度指标粒子的全局最优点gbest
t
、惯性因子、学习因子c1、学习因子c2，计算出优化后的精度指标粒子的速度v
kt+1
；
步骤73：根据优化后的精度指标粒子的速度v
kt+1
，计算出优化后的精度指标粒子的位置g
kt+1
。
35.所述步骤8具体包括以下步骤：步骤81：以精度指标为粒子，初始化粒子群粒子，根据设定的初始精度指标、最大适应度和最大迭代次数k、最大速度系数，计算出相对最大速度v
max
；步骤82：根据设定的初始精度指标，计算出初始局部最优点pbest
k0
、初始全局最优点gbest0、精度指标粒子的初始位置g
k1
、精度指标粒子的初始速度v
k1
。
36.所述步骤11的具体操作为：判断是否满足结束条件，若全局最优点gbest
t
的适应度大于最大适应度或迭代次数t大于最大迭代次数k，则满足结束条件，将当前精度指标作为最优精度指标，结束迭代，若不满足结束条件，则将迭代次数加一，重复步骤9，直至满足结束条件。
37.本实施例的其他部分与上述实施例1相同，故不再赘述。
38.实施例3：本实施例在上述实施例1-2任一项的基础上，如图1、图2、图3、图4、图5、图6所示，以龙门式ca摆角五轴数控机床为例进行详细说明。
39.本实施例提出一种面向大部件群孔加工设备精度指标优化方法，该方法将大部件群孔加工设备精度指标计算进行了定量化，包括以下步骤：步骤s1：确定数控机床类型，建立数控机床的拓扑结构；步骤s2：建立数控机床刀尖点运动的空间定位误差模型；步骤s3：建立数控机床群孔加工孔位误差模型；步骤s4：提取误差模型中的机床的精度指标；步骤s5：建立精度指标评价模型；步骤s6：建立适应度模型；步骤s7：设备精度指标优化算法粒子参数与模型定义；步骤s8：以精度指标为粒子，初始化粒子群粒子及相关参数；步骤s9：根据粒子群算法，计算粒子适应度，局部最优点及全局最优点；步骤s10：更新粒子速度和粒子位置；步骤s11：判断是否满足算法结束条件，满足则该组指标为最优，不满足则重复步骤（8），直至满足条件。
40.步骤s1中所述确定数控机床类型，建立数控机床的拓扑结构，根据低序体阵列法对机床的拓扑结构进行简化，并对相对应运动单元进行低序体编号，本实施例以龙门式ca摆角五轴数控机床为例，机床床身编号0、x轴运动单元编号1、y轴运动单元编号2、z轴运动单元编号3、c轴运动单元编号4、a轴运动单元编号5及刀具单元编号6，其中z轴运动单元的负方向作为加工产品方向，运动范围为：
41.所述步骤s2中建立数控机床刀尖点运动的空间定位模型，通过机床拓扑结构建立相邻运动单元之间的运动传递矩阵：
42.结合刀尖点的初始定位：
43.其中，l为刀长；计算得出机床在进给轴输入参数作用下机床刀尖的空间定位向量：
44.刀具在a摆动轴坐标系下初始姿态，计算得出机床在进给轴输入参数作用下机床刀尖的空间姿态向量：
45.其中和分别为机床为进给轴的运动矩阵，x、y、z、c、a分别为机床各进给轴的运动参数。
46.考虑在实际工作情况，刀尖点实际运动轨迹是机床理想运动与误差因素综合作用的结果。那么根据多体运动学理论，任意刚体在空间的运动状态可以通过对应方向坐标矩阵的变换关系进行表达，线性进给轴u和转动进给轴v对应的误差传递矩阵可以分别表示为：
47.其中分别为机床各进给轴的误差矩阵，转动进给轴误差情况同线性进给轴误差类似，同时进给轴传递过程中需考虑垂直度误差的影响，那么实际空间定位向量：
48.其中根据龙门式ca摆角五轴数控机床拓扑结构，涉及的垂直度误差包含x轴移动对y轴的垂直度为、z轴移动对x轴的垂直度为、z轴移动对y轴的垂直度为、c轴移动对x轴的垂直度为、c轴移动对y轴的垂直度为、c轴移动对a轴的垂直度为， 通过刀尖点实际位置与理论位置之间距离，建立孔位误差模型：
49.所述步骤s3中建立数控机床群孔加工孔位误差模型，依据制孔孔位误差在与理想轴线垂直平面的投影关系，可以计算在制孔位置处位置偏差的计算表达式为：
50.并求出该机床制孔时孔位误差为：
51.所述步骤s4中提取误差模型中的机床的精度指标，将所有的线性进给轴u和转动进给轴v对应的误差传递矩阵中包含的定位精度指标组成集合g1、重复定位精度指标组成集合g2、直线度定位精度指标组成集合g3、角度偏差指标组成集合g4、垂直度指标组成集合g5。
52.所述步骤s5中建立精度指标评价公式，将这些指标集合求并集，组成为该设备的精度指标集合,，m为指标集合个数：
53.将这些精度指标构成向量表示为：
54.将设备制孔过程中所制孔总数设为n，每个孔各轴进给量设为：
55.每个点位的孔位误差为：
56.所有加工孔位的平均误差设为：
57.本实施例提出了基于粒子群算法的机床精度指标优化方法，利用改进的粒子群算法可找出一组指标使得精度指标尽可能宽松且满足条件，p为设备制孔精度要求。
58.构造单组精度指标评价系数计算公式：
59.其中为均值，s为标准差，为控制系数，用于减小标准差太小时带来的影响。该公式用于评价每组精度指标的好坏，指标均值越大、标准差越小指标的评分越高，反之评分越低。
60.总体精度指标评价系数计算公式，以每组精度指标评价系数的均值作为总体精度指标评价系数：
61.所述步骤s6中建立适应度公式，构造适应度计算公式：
62.该公式用于粒子群算法中适应度的计算，当综合误差满足精度要求p时，适应度即为总体精度指标评价系数，当综合误差不满足要求时适应度为负的综合误差。构建完适应度公式后，即可使用粒子群算法对设备精度指标进行优化。
63.所述步骤s7中设备精度指标优化算法粒子参数与公式定义，设粒子群粒子个数为m，适应度计算函数为，在第t次迭代时每个粒子的位置为设备的精度指标，每个精度指标粒子的速度为
64.定义局部最优点：
65.定义全局最优点：
66.定义精度指标粒子的相对最大速度，设粒子预设精度指标为，相对最大速度定义为，c为相对最大速度系数。
67.定义精度指标粒子速度更新公式：
[0068][0069]
那么每个精度指标粒子更新后的速度为:
[0070]
其中为学习因子，为惯性因子，为随机数。
[0071]
定义位置更新公式：
[0072]
所述步骤s8中以精度指标为粒子，初始化粒子群粒子及相关参数，给定设备初始精度指标，并设置最大适应度和最大迭代次数k，设置相对在最大速度系数，并计算相对最大速度：
[0073]
设置初始局部最优点：
[0074]
全局最优点：
[0075]
设置初始化m个精度指标粒子，迭代次数t=1，位置为：
[0076]
初始速度为：
[0077]
所述步骤s9中根据粒子群算法，计算粒子适应度，局部最优点及全局最优点。
[0078]
计算局部最优点：
[0079]
计算全局最优点：
[0080]
所述步骤s10）中更新精度指标粒子速度和粒子位置；利用计算出的全局最优点和、局部最优点、上一轮迭代的例子速度，依据速度更新公式对精度指标粒子速度进行更新：
[0081][0082]
每个精度指标粒子更新后的速度为:
[0083]
依据精度指标粒子更新后的速度更新精度指标粒子的位置：
[0084]
所述步骤s11中判断是否满足算法结束条件：或迭代次数t》k，若满足条件，找到一组设备精度指标，算法结束，否则迭代次数t=t+1，重复步骤s9。
[0085]
按照本实施例的上述方法，可找到一组设备精度指标使得平均误差，既能满足加工需求，又能使尽可能接近p，利用到设备全部的加工能力，且找到的定位精度指标集合g1、重复定位精度指标集合g2、直线度定位精度指标集合g3、角度偏差指标集合g4、垂直度指标集合g5，在每组指标集合内指标的分布都相对均匀有着较小的标准差。
[0086]
本实施例的其他部分与上述实施例1-2任一项相同，故不再赘述。
[0087]
实施例4：本实施例在上述实施例1-3任一项的基础上，提出一种面向大部件群孔加工设备精度指标优化方法，基于改进型粒子群算法对机床精度指标进行优化，如图1所示，具体包括以下步骤：步骤（1），根据低序体阵列法对机床的拓扑结构进行简化，并对相对应运动单元进行低序体编号，本实施例以龙门式ca摆角五轴数控机床为例，机床床身编号0、x轴运动单元编号1、y轴运动单元编号2、z轴运动单元编号3、c轴运动单元编号4、a轴运动单元编号5及刀具单元编号6，其中z轴运动单元的负方向作为加工产品方向，运动范围为。
[0088]
步骤（2）和步骤（3），如图2所示，通过机床拓扑结构建立相邻运动单元之间的运动传递矩阵，结合刀尖点的初始定位，l为刀长，计算得出机床在进给轴输入参数作用下机床刀尖的空间定位向量，如式(1)所示：
ꢀꢀꢀ
（1）刀具在a摆动轴坐标系下初始姿态，计算得出机床在进给轴输入参数作用下机床刀尖的空间姿态向量，如式(2)所示：
ꢀꢀꢀ
（2）
其中式(1)和式(2)中：
[0089][0090][0091][0092][0093]
分别为机床为进给轴的运动矩阵，x、y、z、c、a分别为机床各进给轴的运动参数。
[0094]
考虑在实际工作情况，刀尖点实际运动轨迹是机床理想运动与误差因素综合作用的结果。那么根据多体运动学理论，任意刚体在空间的运动状态可以通过对应方向其次坐标矩阵的变换关系进行表达，线性进给轴u和转动进给轴v对应的误差传递矩阵，如式(3)所示：
ꢀꢀꢀ
（3）其中式(3)中：
[0095]
其中分别为机床各进给轴的误差矩阵，转动进给轴误差情况同线性进给轴误差类似，同时进给轴传递过程中需考虑垂直度误差的影响，那么实际空间定位向量，如式(4)所示：（4）其中根据龙门式ca摆角五轴数控机床拓扑结构，涉及的垂直度误差如下：x轴移动对y轴的垂直度为：
[0096]
z轴移动对x轴的垂直度为：
[0097]
z轴移动对y轴的垂直度为：
[0098]
c轴移动对x轴的垂直度为：
[0099]
c轴移动对y轴的垂直度为：
[0100]
c轴移动对a轴的垂直度为：
[0101]
结合式(1)和式(4)，通过刀尖点实际位置与理论位置之间距离，建立孔位误差模型，如式(5)所示：
ꢀꢀ
（5）依据制孔孔位误差在与理想轴线垂直平面的投影关系，结合式(5)，计算在制孔位置处位置偏差，如式(6)所示：
ꢀꢀꢀ
（6）结合式(6)，求出该机床制孔时孔位误差为：
ꢀꢀꢀ
（7）步骤（4）、（5）及（6），将所有的线性进给轴u和转动进给轴v对应的误差传递矩阵中包含的定位精度指标组成集合g1、重复定位精度指标组成集合g2、直线度定位精度指标组成集合g3、角度偏差指标组成集合g4、垂直度指标组成集合g5。
[0102]
将这些指标集合求并集，组成为该设备的精度指标集合，如式(8)所示，m为指标集合个数：
ꢀꢀꢀ
（8）将这些精度指标构成向量如式(9)所示：
ꢀꢀꢀ
（9）将设备制孔过程中所制孔总数设为n，每个孔各轴进给量如式(10)所示：
ꢀꢀꢀ
（10）结合式(7)式(9)及式(10)，每个点位的孔位误差如式(11)所示：
ꢀꢀꢀ
（11）结合式(11)，所有加工孔位的平均误差如式(12) 所示：
ꢀꢀꢀ
（12）结合式(12)，利用改进的粒子群算法可找出一组指标使得精度指标尽可能宽松且满足条件，p为设备制孔精度要求，造单组精度指标评价系数计算公式，如式3 (1)所示：
ꢀꢀꢀ
（13）其中为均值，s为标准差，为控制系数，用于减小标准差太小时带来的影响。该公式用于评价每组精度指标的好坏，指标均值越大、标准差越小指标的评分越高，反之评分越低。
[0103]
结合式(13)，总体精度指标评价系数计算公式，以每组精度指标评价系数的均值作为总体精度指标评价系数：
ꢀꢀꢀ
（14）结合式式(12)、式(14)，构造适应度计算公式，如式(15)所示：
ꢀꢀꢀ
（15）式(15)用于粒子群算法中适应度的计算，当综合误差满足精度要求p时，适应度即为总体精度指标评价系数，当综合误差不满足要求时适应度为负的综合误差。构建完适应度公式后，即可使用粒子群算法对设备精度指标进行优化。
[0104]
步骤（7），结合式(15)设粒子群精度指标粒子个数为m，适应度计算函数为，在第t次迭代时每个精度指标粒子的位置为设备的精度指标，如式(16)所示：
ꢀꢀꢀ
（16）每个精度指标粒子的速度，如式(17)所示：
ꢀꢀꢀ
（17）结合式(16)，定义精度指标局部最优点：
ꢀꢀꢀ
（18）结合式(18)定义精度指标全局最优点： （19）结合式(17)，定义精度指标粒子的相对最大速度，设粒子预设精度指标如式(20)所示：
ꢀꢀꢀ
（20）相对最大速度定义如式(21)所示，其中c为相对最大速度系数：
ꢀꢀꢀ
（21）定义速度更新公式如式(22)所示：
[0105]
ꢀꢀꢀ
（22）那么结合(22)式每个精度指标粒子更新后的速度如式(23)所示:
ꢀꢀꢀ
（23）其中为学习因子，为惯性因子，为随机数。
[0106]
结合式(20)和式(23)定义指标粒子位置更新公式如式(24)：（24）步骤（8），以精度指标为粒子，初始化粒子群粒子及相关参数，给定设备初始精度指标，并设置最大适应度和最大迭代次数k，设置相对最大速度系数，并计算相对最大速度如式(25)所示：
ꢀꢀꢀ
（25）结合式(25)，设置设备初始精度指标局部最优点如式(26)所示：
ꢀꢀꢀ
（26）设备初始精度指标全局最优点如式(27)所示：
ꢀꢀꢀ
（27）设置初始化m个精度指标粒子，迭代次数t=1，初始位置如式(28)所示：
ꢀꢀꢀ
（28）初始速度如式(29)所示：
ꢀꢀꢀ
（29）随后进行步骤（9），计算精度指标粒子适应度，精度指标局部最优点及精度指标全局最优点。
[0107]
结合式(13)(14)(15)，计算设备精度指标局部最优点，如式(30)所示：
ꢀꢀꢀ
（30）结合式(13)(14)(15)，计算设备精度指标全局最优点，如式(31)所示：
ꢀꢀꢀ
（31）结合式(30)、式(31)，进行步骤（10），利用计算出的精度指标全局最优点、精度指标局部最优点和上一轮迭代的设备精度指标粒子速度，依据速度更新公式对精度指标粒子速度进行更新，如式(32)所示：
[0108]
ꢀꢀꢀ
（32）每个精度指标粒子速度更新示意图如图3所示，式(32)表示计算出粒子速度后，通过每个精度指标粒子的相对最大速度对粒子速度进行限制，当计算出的粒子速度且时，将粒子速度限制为，当例子速度且时，将粒子速度限制为，更新后的速度如式(33)所示:
ꢀꢀꢀ
（33）依据精度指标粒子速度更新粒子位置，示意图如图4所示，更新后位置如式(34)所示：
ꢀꢀꢀ
（34）对所有精度指标粒子进行速度更新和位置更新后，示意图如图5所示，最后进行步骤（11）判断是否满足算法结束条件：或迭代次数t》k，若满足条件，找到一组设备精度指标，算法结束，否则迭代次数t=t+1，重复步骤（9）。按照上述方法，经过多次迭代后全局最优点逐渐逼近实际最优点，如图6所示，此时可找到一组设备精度指标使得平均误差，既能满足加工需求，又能使尽可能接近p，利用到设备全部的加工能力，且找到的定位精度指标集合g1、重复定位精度指标集合g2、直线度定位精度指标集合g3、角度偏差指标集合g4、垂直度指标集合g5，在每组指标集合内指标的分布都相对均匀有着较小的标准差。
[0109]
本实施例的其他部分与上述实施例1-3任一项相同，故不再赘述。
[0110]
以上所述，仅是本发明的较佳实施例，并非对本发明做任何形式上的限制，凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化，均落入本发明的保护范围之内。

技术特征：
1.一种面向大部件群孔加工设备精度指标优化方法，其特征在于，包括以下步骤：步骤1：根据数控机床类型，建立数控机床的拓扑结构；步骤2：根据步骤1建立的数控机床的拓扑结构，建立数控机床刀尖点运动的空间定位误差模型p
error
；步骤3：根据步骤2建立的空间定位误差模型p
error
，以及制孔孔位误差与理想轴线垂直平面的投影关系，建立数控机床群孔加工孔位误差模型；步骤4：根据步骤3建立的数控机床群孔加工孔位误差模型，提取数控机床群孔加工孔位误差模型中的机床的精度指标；步骤5：计算出步骤4提取出的精度指标的合集，建立数控机床精度指标评价模型；步骤6：根据步骤5建立的精度指标评价模型，判断是否满足设备制孔精度要求p，并根据判断结果建立精度指标适应度模型；步骤7：根据步骤6建立的精度指标适应度模型，计算出优化后的精度指标粒子速度v
kt+1
和优化后的精度指标粒子位置g
kt+1
；步骤8：以精度指标为粒子，设置最大适应度α和最大迭代次数k，初始化粒子群粒子，计算出精度指标粒子的初始位置g
k1
和精度指标粒子的初始速度v
k1
；步骤9：根据步骤7计算出的优化后的精度指标粒子速度v
kt+1
和优化后的精度指标粒子位置g
kt+1
和步骤8计算出的精度指标粒子的初始位置g
k1
和精度指标粒子的初始速度v
k1
，计算出精度指标粒子的局部最优点pbest
kt
、精度指标粒子的全局最优点gbest
t
；步骤10：根据步骤9计算出的全局最优点gbest
t
和局部最优点pbest
kt
以及上一轮迭代的精度指标粒子速度，更新精度指标粒子速度和精度指标粒子位置，得到更新后的精度指标粒子速度v’kt+1
、更新后的精度指标粒子位置x
kt+1
；步骤11：判断是否满足结束条件，若满足结束条件，则将当前精度指标作为为最优的精度指标，不满足则重复步骤9，直至满足条件。2.如权利要求1所述的一种面向大部件群孔加工设备精度指标优化方法，其特征在于，所述步骤2具体包括以下步骤：步骤21：根据步骤1建立的机床拓扑结构，建立相邻运动单元之间的运动传递矩阵m
ij
；步骤22：根据步骤21建立的相邻运动单元之间的运动传递矩阵m
ij
，结合刀尖点的初始定位p
initial
，计算出机床在进给轴输入参数作用下机床刀尖的空间定位向量p
ideal
；步骤23：根据步骤21建立的相邻运动单元之间的运动传递矩阵m
ij
，结合刀具在a摆动轴坐标系下的初始姿态v
initial
，计算出机床在进给轴输入参数作用下机床刀尖的空间姿态向量v
ideal
；步骤24：根据机床进给轴的误差矩阵，计算出线性进给轴u的误差传递矩阵e
uij
和转动进给轴v的误差传递矩阵e
vij
；步骤25：根据步骤24计算出的线性进给轴u的误差传递矩阵e
uij
、刀尖点的初始定位p
initial
、运动传递矩阵m
ij
，计算出实际空间定位向量p
actual
；步骤26：根据步骤22计算出的空间定位向量p
ideal
、步骤25计算出的实际空间定位向量p
actual
，建立数控机床刀尖点运动的空间定位误差模型p
error
。3.如权利要求2所述的一种面向大部件群孔加工设备精度指标优化方法，其特征在于，所述步骤3具体包括以下步骤：
步骤31：根据制孔孔位误差在与理想轴线垂直平面的投影关系，步骤23计算出的空间姿态向量v
ideal
、步骤26建立的空间定位误差模型p
error
，计算出制孔位置处位置偏差p
error_v
；步骤32：根据步骤31计算出的制孔位置处位置偏差p
error_v
，计算出当前机床制孔时孔位误差r
e
，建立数控机床群孔加工孔位误差模型。4.如权利要求3所述的一种面向大部件群孔加工设备精度指标优化方法，其特征在于，所述步骤4的具体操作为：将线性进给轴u的误差传递矩阵e
uij
和转动进给轴v的误差传递矩阵e
vij
中包含的定位精度指标组成集合g1、重复定位精度指标组成集合g2、直线度定位精度指标组成集合g3、角度偏差指标组成集合g4、垂直度指标组成集合g5。5.如权利要求4所述的一种面向大部件群孔加工设备精度指标优化方法，其特征在于，所述步骤5具体包括以下步骤：步骤51：计算出定位精度指标组成集合g1、重复定位精度指标组成集合g2、直线度定位精度指标组成集合g3、角度偏差指标组成集合g4、垂直度指标组成集合的合集g5，得到当前设备的精度指标集合q；步骤52：根据构成当前设备的精度指标集合q的精度指标向量g、制孔过程中所制每个孔各轴进给量p
i
，计算出每个制孔点位的孔位误差r
e
(g，p
i
)；步骤53：根据步骤52计算出的每个制孔点位的孔位误差r
e
(g，p
i
)，计算出所有加工孔位的平均误差r
ea
；步骤54：判断平均误差r
ea
是否满足设备制孔精度要求p，若满足，计算出单组精度指标评价系数j
i
(g
i
)；步骤55：根据步骤54计算出的单组精度指标评价系数j
i
(g
i
)，计算出总体精度指标评价系数j(g)，建立数控机床精度指标评价模型。6.如权利要求5所述的一种面向大部件群孔加工设备精度指标优化方法，其特征在于，所述步骤6的具体操作为：根据步骤5建立的精度指标评价模型，判断所有加工孔位的平均误差r
ea
是否满足设备制孔精度要求p，若所有加工孔位的平均误差r
ea
满足精度要求p，则适应度为总体精度指标评价系数j(g)，若所有加工孔位的平均误差r
ea
不满足精度要求p，则适应度为负的所有加工孔位的平均误差r
ea
，根据判断结果建立精度指标适应度模型。7.如权利要求6所述的一种面向大部件群孔加工设备精度指标优化方法，其特征在于，所述步骤7具体包括以下步骤：步骤71：根据步骤6建立的建立精度指标适应度模型、第t次迭代时设备的精度指标g
kt
、每个精度指标粒子的速度v
kt
，计算出精度指标粒子的局部最优点pbest
kt
、精度指标粒子的全局最优点gbest
t
；步骤72：根据步骤71计算出的精度指标粒子的局部最优点pbest
kt
、精度指标粒子的全局最优点gbest
t
、惯性因子、学习因子c1、学习因子c2，计算出优化后的精度指标粒子的速度v
kt+1
；步骤73：根据优化后的精度指标粒子的速度v
kt+1
，计算出优化后的精度指标粒子的位置g
kt+1
。8.如权利要求7所述的一种面向大部件群孔加工设备精度指标优化方法，其特征在于，所述步骤8具体包括以下步骤：步骤81：以精度指标为粒子，初始化粒子群粒子，根据设定的初始精度指标、最大适应
度α和最大迭代次数k、最大速度系数，计算出相对最大速度v
max
；步骤82：根据设定的初始精度指标，计算出初始局部最优点pbest
k0
、初始全局最优点gbest0、精度指标粒子的初始位置g
k1
、精度指标粒子的初始速度v
k1
。9.如权利要求8所述的一种面向大部件群孔加工设备精度指标优化方法，其特征在于，所述步骤11的具体操作为：判断是否满足结束条件，若全局最优点gbest
t
的适应度α大于最大适应度或迭代次数t大于最大迭代次数k，则满足结束条件，将当前精度指标作为最优的精度指标，结束迭代，若不满足结束条件，则将迭代次数加一，重复步骤9，直至满足结束条件。10.如权利要求1所述的一种面向大部件群孔加工设备精度指标优化方法，其特征在于，所述步骤1的具体操作为：根据低序体阵列法对机床的拓扑结构进行简化，并对相对应运动单元进行低序体编号。

技术总结
本发明涉及机械加工领域，具体地说，涉及一种面向大部件群孔加工设备精度指标优化方法，通过从数控机床群孔加工孔位误差模型中提取了不同种类的精度指标，并根据这些指标，建立了数控机床精度指标评价模型、精度指标适应度模型，并根据建立的模型，以精度指标为粒子，通过计算出精度指标粒子的全局最优点gbest


技术研发人员：李杰 刘东君 马振博 刘大炜 蒋超 樊西锋 王浩 冯如一 张松 李春明 王斌
受保护的技术使用者：成都飞机工业（集团）有限责任公司
技术研发日：2023.08.08
技术公布日：2023/9/9