面向航天器大型舱体构件加工的多机器人协同规划方法

未命名 08-12 阅读:192 评论:0


1.本发明涉及多机器人协同规划技术领域,具体涉及面向航天器大型舱体构件加工的多机器人协同规划方法。


背景技术:

2.大型航天器舱体结构件具有尺寸大、刚性弱等特点,且待加工支架数量多,不均匀地分布在舱体外表面。目前,此类构件主要采用龙门式多轴数控机床进行加工,但数控机床存在尺寸大、制造成本高、加工对象单一、加工行程有限等问题,且在该种加工模式下支架需要反复拆装,影响了航天器大型结构件的加工精度和加工效率。针对上述制造问题,一种基于多机器人系统的原位制造模式应运而生,即舱体固定、机器人移动至不同位置对舱体进行加工,该种加工模式适应性强,可以加工尺寸跨度较大的零件,同时也避免了舱体支架的反复拆装,采用多机器人协同作业可以大幅提高大型航天器舱体构件产品的生产效率。
3.但由于舱体构件尺寸大,受机器人加工范围限制,在不改变舱体旋转角度的情况下机器人无法加工到舱体周向的所有支架,往往需要借助变位器对舱体旋转使得机器人能够覆盖到所有支架任务,而舱体旋转变位规划归结于为分区加工制造问题。文献“ma j w,lu x,li g l,et al.toolpath topology design based on vector field of tool feeding direction in sub-regional processing for complex curved surface[j].journal of manufacturing processes,2020,52:44-57.”将复杂曲面划分为子区域,根据曲面的几何特征或加工技术确定细分准则,从而实现不同子区域内可变工艺参数的加工。该方法主要面向于机床加工模式,且主要集中在对于形状确定、加工特征单一的大型部件局部区域的分区研究,对于具有待加工支架朝向与分布随舱体的旋转而发生改变特点的航天器舱体结构难以适用。而且对于航天器舱体结构件,其尺寸大、加工范围广,在加工时机器人往往需要改变站位,一是为了能够加工到所有任务;二是为了调整加工任务时的姿态,使得机器人团队加工任务的效率和质量得到提升。文献“wei t,dai j,zhou w,et al.process planning and control technology on multi-station working mode of robot drilling and riveting system with auxiliary axis[j].china mechanical engineering,2014,25(01):23-27.”提出了一种面向工业机器人的分站式规划策略,对机器人的转站机制进行了建模与优化。该方法充分考虑了单机器人的移动换站问题,但忽略了多机器人协同作业过程中的空间争用问题,无法直接应用于大型舱体构件多机器人协同加工这种存在作业空间争用的情形。任务分配的均衡性对多机器人团队的作业效率起到了至关重要的作用,文献“wang z,gombolay m.learning scheduling policies for multi-robot coordination with graph attention networks[j].ieee robotics and automation letters,2020,5(3):4509-4516”针对该问题提出了一种基于网络的多机器人协调学习调度策略,该方法可对各种规模的机器人团队进行快速调度优化,并从多任务中找到高质量的解决方案,但此方法只适用于任务固定、加工信息不随时间改变的作业,无法处理舱外支架的朝向随舱体旋转而发生改变的复杂任务情况。
[0004]
现有发明存在的不足之处:
[0005]
1.针对大型舱体构件的加工制造需求,现有发明多采用数控机床加工模式,该模式存在制造成本高、加工行程受限等问题,且制造周期长,无法满足大型舱体构件的高效高精制造需求。
[0006]
2.现有发明适用于具有单一特征且加工状态信息固定的加工对象,而对于加工设备为多台机器人、加工对象为大型弱刚性薄壁构件的双弱结构现有发明鲜有提及。


技术实现要素:

[0007]
发明目的:本发明所要解决的技术问题是针对现有技术的不足,提供面向航天器大型舱体构件加工的多机器人协同规划方法,包括:
[0008]
步骤1,导入航天器大型舱体三维数模,获取航天器舱体数据和舱体表面待加工支架数据,所述航天器舱体数据包含舱体长度l和舱体半径r;所述舱体表面待加工支架数据包含待加工支架的个数n和第i个待加工支架的位姿ti={xi,yi,zi,ai,bi,ci},i=1,2,

,n,其中xi、yi,、zi为第i个待加工支架的空间位置坐标,ai、bi、ci分别为旋转角、俯仰角和偏转角,形成待加工支架任务矩阵t={t1,t2,

,tn};
[0009]
步骤2,采用蒙特卡洛法求解机器人工作空间包络,以能够覆盖舱体长度l为约束确定多机器人作业系统单元组成,使多机器人系统能够加工舱体轴向方向的所有待加工支架;
[0010]
步骤3,基于步骤1获取的航天器舱体数据和舱体表面待加工支架数据,以机器人运动灵巧性指标可操作度ω=(det j(θ)j(θ)
t
)
1/2
为评价指标,确定多机器人系统的空间布局方案,使多机器人系统的运动灵活性最优,其中j(θ)为机器人的雅可比矩阵,det是计算矩阵行列式的函数,t是矩阵转置;
[0011]
步骤4,设置舱体旋转角度步长α(单位:
°
),确定舱体旋转位置数量k以及舱体旋转位置集合{r1,r2,

,rk},判断机器人对k个舱体旋转位置下支架位姿的可达性,得到在k个舱体旋转位置下的可加工任务集合{x1,x2,

,xk},rk表示第k个舱体旋转位置,xk表示第k个舱体旋转位置下的可加工任务集合,采用机器人笛卡尔空间刚度矩阵k=j(θ)-tkθ
j(θ)-1
对机器人加工k个舱体旋转位置下的可加工任务集合{x1,x2,

,xk}的加工性能进行评价,其中k
θ
为机器人的关节刚度矩阵,得到k个舱体旋转位置下的机器人加工任务的刚度值{k1,k2,

,kk},kk表示第k个舱体旋转位置下的机器人加工任务的刚度值;
[0012]
步骤5,获得在所有舱体旋转位置下的总加工任务集合xr={x1∪x2∪

∪xk},采用遗传算法从舱体旋转位置集合{r1,r2,

,rk}中筛选能够加工到总加工任务集合xr的所有舱体旋转变位组合方案{rot1,rot2,

,rotn},rotn表示第n个舱体旋转变位组合方案,其中任意一个可加工到总加工任务集合xr的舱体旋转变位组合方案roti={ra,rb,

,rc}包括ki个舱体旋转角度,ra、rb、rc分别代表舱体第a、b、c个旋转位置,其中a,b,c≤k且a,b,c∈n*,k
min
≤ki≤k,k
min
为舱体的最少旋转变位次数,采用舱体旋转变位性能评价函数rot_fun=(ka+kb+

+kc)/ki,其中,ka、kb、kc分别表示第a、b、c个舱体旋转位置下的机器人加工任务的刚度值,选取rot_fun结果最优的方案rot={ra,rb,

,rc}为最终舱体旋转变位方案;
[0013]
步骤6,在选定最优的方案rot={ra,rb,

,rc}的基础上,设定机器人沿舱体轴向方向的移动步长s,设定机器人可选择站位的数量j及机器人站位集合{p1,p2,

,pj},pj表
示第j个机器人站位,遍历所有机器人站位,判断机器人对j个机器人站位下支架位姿的可达性,得到在舱体旋转位置ra下机器人站位{p1,p2,

,pj}对应的可加工任务集合{x
1a
,x
2a
,

,x
ja
},并采用机器人笛卡尔空间刚度矩阵k=j(θ)-tkθ
j(θ)-1
对机器人加工j个机器人站位下的可加工任务集合{x
1a
,x
2a
,

,x
ja
}的加工性能进行评价,得到j个机器人站位下的机器人加工任务的刚度值{k
1a
,k
2a
,

,k
ja
},k
ja
表示第j个机器人站位下的机器人加工任务的刚度值;
[0014]
获得舱体旋转位置rb、rc下机器人站位{p1,p2,

,pj}对应的可加工任务集合{x
1b
,x
2b
,

,x
jb
}、{x
1c
,x
2c
,

,x
jc
},以及得到对应在j个机器人站位下的机器人加工任务的刚度值{k
1b
,k
2b
,

,k
jb
}、{k
1c
,k
2c
,

,k
jc
};
[0015]
步骤7,根据最优的方案rot={ra,rb,

,rc},以及当前舱体旋转位置ra下所有机器人站位{p1,p2,

,pj}对应的总加工任务集合xa={x
1a
∪x
2a


∪x
ja
},筛选能够加工到当前舱体旋转位置ra下的所有任务xa的机器人站位组合方案{pos
1a
,pos
2a
,

,pos
naa
},并得到对应机器人站位组合方案的刚度值{k
1a
,k
2a
,

,k
naa
},pos
naa
表示在舱体旋转位置ra下的第na个机器人站位组合方案,k
naa
表示在机器人站位组合方案pos
naa
下的机器人加工任务的刚度值,其中在舱体旋转位置ra下任意一个机器人站位组合方案pos
ia
={pd,pe,

,pf}包括m个机器人站位,pd、pe、pf分别表示第d、e、f个机器人站位,d,e,f≤j且d,e,f∈n*;
[0016]
获得舱体旋转位置rb、rc下机器人站位{p1,p2,

,pj}对应的总加工任务集合xb={x
1b
∪x
2b


∪x
jb
}、xc={x
1c
∪x
2c


∪x
jc
},筛选能够加工到当前舱体旋转位置rb、rc下的所有任务xb、xc的机器人站位组合方案{pos
1b
,pos
2b
,

,pos
nbb
}、{pos
1c
,pos
2c
,

,pos
ncc
},并得到对应机器人站位组合方案的刚度值{k
1b
,k
2b
,

,k
nbb
}、{k
1c
,k
2c
,

,k
ncc
},pos
nbb
表示在舱体旋转位置rb下的第nb个机器人站位组合方案,k
ncc
表示在机器人站位组合方案pos
ncc
下的机器人加工任务的刚度值;
[0017]
根据最优的方案rot={ra,rb,

,rc},其中任意一个机器人站位方案posi={pos
ga
,pos
hb
,

,pos
lc
},g,h,l≤j且g,h,l∈n*,pos
ga
代表在舱体旋转位置ra时采用机器人站位方案pos
ga
,pos
hb
代表在舱体旋转位置rb时采用机器人站位方案pos
hb
,pos
lc
代表在舱体旋转位置rc时采用机器人站位方案pos
lc
,计算机器人站位方案{pos
ga
,pos
hb
,

,pos
lc
}下的刚度值{k
ga
,k
hb
,

,k
lc
},采用机器人站位性能评价函数pos_fun=k
ga
+k
hb
+

+k
lc
,选取pos_fun结果最优的方案pos={pos
ga
,pos
hb
,

,pos
lc
}为最终机器人站位方案;
[0018]
步骤8,在最终机器人站位方案pos={pos
ga
,pos
hb
,

,pos
lc
}的基础上,得到机器人ri的可加工任务集合x
ri
={x
riag
,x
ribh
,

,x
ricl
},i=1,2,

,m,x
riag
为机器人ri在舱体旋转位置ra时采用机器人站位方案pos
ga
时的可加工任务集合,x
ribh
为机器人ri在舱体旋转位置rb时采用机器人站位方案pos
hb
时的可加工任务集合,x
ricl
为机器人ri在舱体旋转位置rc时采用机器人站位方案pos
lc
时的可加工任务集合,进而对各个机器人的加工任务进行任务分配;
[0019]
步骤9,基于市场原理的任务拍卖算法对加工任务进行分配,得到机器人ri的任务分配结果a
ri
={a
riag
,a
ribh
,

,a
ricl
},a
riag
为机器人ri在舱体旋转位置ra、机器人站位方案pos
ga
下的加工任务分配结果,a
ribh
为机器人ri在舱体旋转位置rb、机器人站位方案pos
ha
下的加工任务分配结果,a
ricl
为机器人ri在舱体旋转位置rc、机器人站位方案pos
la
下的加工任务分配结果,优化机器人加工任务的顺序,形成各个机器人的加工任务时序调度表,最终形
成多机器人团队的加工调度方案。
[0020]
进一步的,步骤1中,待加工支架的位姿会随舱体的旋转而发生改变,进而影响机器人对待加工支架的可达性。
[0021]
进一步的,步骤2中,所述多机器人作业系统单元组成包括机器人型号、agv型号、机器人数量m、末端执行器构造、刀具型号和规格,其中,航天器舱体的尺寸决定了机器人型号、agv型号和机器人数量m,舱体表面待加工支架的工艺要求决定了末端执行器构造、刀具型号和规格。
[0022]
进一步的,步骤3中,所述多机器人系统的空间布局方案包括机器人在舱体两侧的分布情况、机器人与舱体的相对位置关系;所述机器人在舱体两侧的分布情况影响机器人之间的相对位置关系,进而影响任务分配与任务时序调度环节;所述机器人与舱体的相对位置关系指机器人距舱体的轴线的水平距离以及舱体轴线距离地面的高度,机器人距离舱体的远近影响机器人的加工姿态,从而影响机器人的加工质量。
[0023]
进一步的,步骤4中,舱体旋转位置数量k计算公式如下:
[0024]
k=360/α。
[0025]
进一步的,步骤5中,舱体的最少旋转变位次数k
min
是指在能够覆盖总加工任务集合xr中所有任务前提下的最少的舱体旋转变位次数。
[0026]
进一步的,步骤6中,机器人可选择站位的数量j计算公式如下:
[0027]
j=2l/s
[0028]
所述机器人可选择站位的数量j与舱体的长度l(单位:m)与机器人沿舱体轴向方向的移动步长s(单位:m)有关,存在j=2l/s的等式关系,可选择站位的设置要保证各个站位之间的间隔距离,保障两个以上机器人同时加工过程的安全性。
[0029]
进一步的,步骤8中,所述对各个机器人的加工任务进行任务分配是指去除“重叠任务”并将“唯一任务”分配给对应机器人;“重叠任务”指可被同侧机器人加工以及能够在不同舱体旋转角度下被加工的任务,“唯一任务”指只能在某一个舱体旋转角度下的某一个机器人站位下被加工的任务。将“唯一任务”分配给对应机器人,计算当前机器人工作量,然后依次将“重叠任务”分配给能加工此重叠任务且工作量最少的机器人,完成多机器人团队的任务分配过程。
[0030]
进一步的,步骤9中,所述各个机器人的加工任务时序调度表是指各个机器人ri加工对应任务集合a
ri
={a
riag
,a
ribh
,

,a
ricl
}中任务的顺序;所述多机器人团队的加工调度方案是指形成以时间为轴,确定在不同时间节点下的舱体旋转角度、机器人加工站位以及在不同站位下的机器人加工任务列表和顺序。
[0031]
本发明还提供了一种存储介质,存储有计算机程序或指令,当所述计算机程序或指令被运行时,实现所述面向航天器大型舱体构件加工的多机器人协同规划方法。
[0032]
本发明还提供了面向航天器大型舱体构件加工的多机器人协同规划系统,包括航天器舱体旋转变位规划系统、多机器人移动换站规划系统、多机器人作业时序规划系统,航天器舱体旋转变位规划系统用于确定舱体的旋转次数以及每次旋转的角度,多机器人移动换站规划系统用于在确定舱体旋转方案下规划各个机器人的换站次数以及空间位置,多机器人作业时序规划系统用于在确定舱体变位方案和机器人换站方案下规划机器人加工任务的数量及顺序。
[0033]
本发明有益效果在于:
[0034]
1、本发明可快速系统的对面向大型构件的多机器人协同加工过程进行规划,包含舱体旋转变位规划、机器人移动换站规划以及加工时序规划,特别适用于大型部件多机器人协同作业领域,具有覆盖范围广、通用性强的特点。
[0035]
2、本发明形成了面向加工过程的多机器人协同作业加工时序调度表,明确指出了机器人在不同时间点的加工任务,对多机器人系统的加工过程具有指导性意义,本发明能够显著提高多机器人系统的加工效率,特别适用于航空航天、船舶等重大工程领域的核心装备制造过程。
附图说明
[0036]
下面结合附图和具体实施方式对本发明做更进一步的具体说明,本发明的上述和/或其他方面的优点将会变得更加清楚。
[0037]
图1为本发明方法流程图。
[0038]
图2为本发明的多机器人加工单元布局图。
[0039]
图3为本发明舱体旋转变位规划系统的流程图。
[0040]
图4为本发明机器人站位规划系统的流程图。
[0041]
图5为本发明的机器人可选择站位设置示例图。
[0042]
图6为本发明加工任务时序规划系统的流程图。
[0043]
图7为本发明多机器人系统加工任务时序规划甘特图。
具体实施方式
[0044]
如图1所示,本发明提供了面向航天器大型舱体构件加工的多机器人协同规划方法,包括:
[0045]
步骤1:设置舱体尺寸为直径2.6m,长度6m,50个支架分布在舱体表面,将支架任务分别编号为1到50,所有的支架都为平面铣削任务,设定每个任务的加工时间均为1min。
[0046]
步骤2:根据舱体尺寸选取4个型号为kuka kr500-3mt的机器人作为加工主体,选取4个agv移动小车作为搭载机器人的移动设备,根据支架的任务类型,选择末端执行器为铣削末端执行器,设定机器人末端的移动速度为1m/s,以上设备构成面向航天器大型舱体的多机器人作业系统。
[0047]
步骤3:根据航天器舱体尺寸和支架任务分布情况,为了覆盖舱体轴线方向的所有任务,选择再舱体轴线方向两侧各布局两个机器人,规定各个机器人只能在舱体一侧进行转站移动,不允许机器人从舱体一侧转站至另一侧,多机器人加工单元布局以及舱体的旋转方向如图2所示。
[0048]
步骤4:舱体旋转变位规划系统的流程如图3所示,首先设置舱体的旋转角度步长,并计算在此步长下的舱体旋转总次数,得到在所有舱体旋转角度下的可加工任务集合。以舱体初始位置为0
°
位置,选取舱体旋转变位的步长为α=30
°
,在此步长下舱体的旋转变位次数n=360/α=12,得到在12个舱体旋转变位角度下的可加工任务集合如表1所示。
[0049]
表1
[0050][0051][0052]
步骤5:然后,基于获得以上12个舱体旋转位置下的总加工任务编号集合结果,筛选可以加工到所有任务的舱体旋转变位组合方式,得到可加工到所有任务的舱体旋转角度组合方式共有12组,以序号组合方式代表舱体的旋转角度组合分别为(1,3,11),(1,9,11),(2,4,6),(2,4,12),(3,5,7),(3,6,12),(3,7,11),(4,6,8),(4,8,12),(5,7,9),(7,9,11),(8,10,12)。
[0053]
进一步,设定舱体旋转方案评价指标,以舱体的旋转变位次数和可加工任务的数量为评价指标,舱体旋转变位次数越少、可加工任务的数量越多,则此种舱体旋转变位方案最优,利用遗传算法优化舱体旋转变位方案的选择过程,得到所有方案中的最优方案。
[0054]
根据上述舱体旋转评价指标选择最优的舱体旋转变位组合序号为(4,6,8),即舱体旋转角度组合方式为(90
°
,150
°
,210
°
),在此种方案下舱体每次的旋转角度及在各个角度可以加工到的任务如表2所示。
[0055]
表2
[0056][0057]
步骤6:机器人站位规划系统的流程图如图4所示,首先,在选定舱体旋转变位方案的基础上,设定机器人的可选择站位数量,确定机器人的站位布局,舱体长度为6m,以舱体的两个端面为起始点与终点,沿舱体轴线方向,步长为1.5m设置机器人的可选择加工站位,舱体每侧各设置5个站位,两侧共计10个可选择站位,分别标号为1到10,站位设置如图5所示,由于机器人布局为舱体轴线两侧各两个机器人,所以机器人选择站位时需要在站位1~5和站位6~10中分别选择两个不同的站位。在确定舱体旋转变位方案的情况下,得到机器人在不同机器人站位下的可加工任务集合如表3所示。
[0058]
表3
[0059][0060]
步骤7:然后,基于预先选定的舱体旋转变位方案,得到可以加工到所有任务的机器人站位组合方案,可加工到所有任务的机器人站位组合方案为{(1,3,6,7)
90
,(2,4,8,9)
150
,(3,4,7,8)
210
},{(1,3,6,8)
90
,(2,4,7,9)
150
,(2,4,7,8)
210
},

,{(2,5,8,10)
90
,(3,5,7,8)
150
,(2,5,8,9)
210
}等,由于满足要求的机器人站位组合方案数量众多,此处只列举了个别机器人站位组合方案,以机器人站位组合方案{(1,3,6,7)
90
,(2,4,8,9)
150
,(3,4,7,8)
210
}为例,此种表述代表在舱体旋转90
°
的位置下,四个机器人的站位分别为站位1、站位3、站位6和站位7,在舱体旋转150
°
的位置下,四个机器人的站位分别为站位2、站位4、站位8和站位9,在舱体旋转210
°
的位置下,四个机器人的站位分别为站位3、站位4、站位7和站位8。
[0061]
进一步,设定机器人站位方案评价指标,以机器人的换站次数、可加工任务的数量和机器人加工任务时的作业刚度为评价指标,机器人的换站次数越少、可加工任务的数量
越多、机器人加工任务时的作业刚度值越大,则此种机器人站位选择方案最优,利用遗传算法优化机器人站位方案的选择过程,得到所有方案中的最优方案。
[0062]
根据上述机器人站位方案评价指标选择机器人的站位组合方式为{(1,3,6,10)
90
,(1,4,6,9)
150
,(2,3,8,9)
210
},即在舱体旋转90
°
的位置下,四个机器人的站位分别为站位1、站位3、站位6和站位10,在舱体旋转150
°
的位置下,四个机器人的站位分别为站位1、站位4、站位6和站位9,在舱体旋转210
°
的位置下,四个机器人的站位分别为站位2、站位3、站位8和站位9。在此种方案下舱体每次的旋转角度及在各个角度可以加工到的任务如表4所示。
[0063]
表4
[0064][0065]
步骤8:加工任务时序规划系统的流程图如图6所示,首先,根据预先选定的舱体旋转变位方案和机器人站位方案得到各个机器人在对应站位下的可加工任务集合。然后,对各个机器人的加工任务进行任务分配,去除可被同侧机器人加工以及可以在不同舱体旋转角度下被加工的“重叠任务”,得到初始任务分配结果。初始任务分配涉及的任务类型为在固定舱体旋转方案和机器人站位方案下,只能在某一个舱体旋转角度下的某一个机器人站位下被加工的任务,该任务类型被定义为“唯一任务”,唯一任务无法再在其他舱体旋转位置下被加工,初始任务分配结果如表5所示,表中包含的任务均为“唯一任务”。
[0066]
表5
[0067][0068]
本加工方案下的“重叠任务”分布情况如表6所示,表示为在该舱体旋转角度的该机器人站位下此任务可以被加工。
[0069]
表6
[0070][0071]
进一步,计算获得初始任务后各个机器人的工作量,任务再分配过程依次将“重叠任务”分配给能加工此任务且工作量最少的机器人,得到最终的任务分配结果如表7所示。
[0072]
表7
[0073][0074]
步骤9:最后,根据各个机器人的任务分配结果,优化机器人在各个舱体旋转角度和站位下的加工任务顺序,找到一个最短的加工任务排序使得完成所有任务时机器人末端移动的距离最短,使得多机器人系统的加工时间最短,进而形成各个机器人的加工任务时序调度表,最终形成多机器人团队的加工调度方案。
[0075]
多机器人团队的加工调度方案是指形成以时间为轴,确定在不同时间节点下的舱体旋转角度、机器人加工站位以及在不同站位下的机器人加工任务列表和顺序,多机器人系统的加工任务时序规划甘特图如图7所示,其中表格中填写的数字为任务的标号,以序号1为例,代表任务1在舱体旋转90
°
情况下,由机器人1位于加工位置1下在0~1分钟进行加工。根据加工任务时序规划甘特图可清晰地观察到多机器人系统不同时间点下的加工状态和多机器人系统的任务完成时间,此种面向大型舱体构件加工的多机器人协同任务规划系统与方法可以为多机器人的加工过程提供科学指导。
[0076]
具体实现中,本技术提供计算机存储介质以及对应的数据处理单元,其中,该计算机存储介质能够存储计算机程序,所述计算机程序通过数据处理单元执行时可运行本发明提供的面向航天器大型舱体构件加工的多机器人协同规划方法的发明内容以及各实施例中的部分或全部步骤。所述的存储介质可为磁碟、光盘、只读存储记忆体(read-only memory,rom)或随机存储记忆体(random access memory,ram)等。
[0077]
本领域的技术人员可以清楚地了解到本发明实施例中的技术方案可借助计算机程序以及其对应的通用硬件平台的方式来实现。基于这样的理解,本发明实施例中的技术方案本质上或者说对现有技术做出贡献的部分可以以计算机程序即软件产品的形式体现出来,该计算机程序软件产品可以存储在存储介质中,包括若干指令用以使得一台包含数据处理单元的设备(可以是个人计算机,服务器,单片机。muu或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
[0078]
本发明提供了面向航天器大型舱体构件加工的多机器人协同规划方法,具体实现该技术方案的方法和途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。本实施例中未明确的各组成部分均可用现有
技术加以实现。

技术特征:
1.面向航天器大型舱体构件加工的多机器人协同规划方法,其特征在于,包括以下步骤:步骤1,导入航天器大型舱体三维数模,获取航天器舱体数据和舱体表面待加工支架数据,所述航天器舱体数据包含舱体长度l和舱体半径r;所述舱体表面待加工支架数据包含待加工支架的个数n和第i个待加工支架的位姿t
i
={x
i
,y
i
,z
i
,a
i
,b
i
,c
i
},i=1,2,

,n,其中x
i
、y
i
,、z
i
为第i个待加工支架的空间位置坐标,a
i
、b
i
、c
i
分别为旋转角、俯仰角和偏转角,形成待加工支架任务矩阵t={t1,t2,

,t
n
};步骤2,采用蒙特卡洛法求解机器人工作空间包络,以能够覆盖舱体长度l为约束确定多机器人作业系统单元组成,使多机器人系统能够加工舱体轴向方向的所有待加工支架;步骤3,基于步骤1获取的航天器舱体数据和舱体表面待加工支架数据,以机器人运动灵巧性指标可操作度ω=(det j(θ)j(θ)
t
)
1/2
为评价指标,确定多机器人系统的空间布局方案,使多机器人系统的运动灵活性最优,其中j(θ)为机器人的雅可比矩阵,det是计算矩阵行列式的函数,t是矩阵转置;步骤4,设置舱体旋转角度步长α,确定舱体旋转位置数量k以及舱体旋转位置集合{r1,r2,

,r
k
},判断机器人对k个舱体旋转位置下支架位姿的可达性,得到在k个舱体旋转位置下的可加工任务集合{x1,x2,

,x
k
},r
k
表示第k个舱体旋转位置,x
k
表示第k个舱体旋转位置下的可加工任务集合,采用机器人笛卡尔空间刚度矩阵k=j(θ)-t
k
θ
j(θ)-1
对机器人加工k个舱体旋转位置下的可加工任务集合{x1,x2,

,x
k
}的加工性能进行评价,其中k
θ
为机器人的关节刚度矩阵,得到k个舱体旋转位置下的机器人加工任务的刚度值{k1,k2,

,k
k
},k
k
表示第k个舱体旋转位置下的机器人加工任务的刚度值;步骤5,获得在所有舱体旋转位置下的总加工任务集合x
r
={x1∪x2∪

∪x
k
},采用遗传算法从舱体旋转位置集合{r1,r2,

,r
k
}中筛选能够加工到总加工任务集合x
r
的所有舱体旋转变位组合方案{rot1,rot2,

,rot
n
},rot
n
表示第n个舱体旋转变位组合方案,其中任意一个可加工到总加工任务集合x
r
的舱体旋转变位组合方案rot
i
={r
a
,r
b
,

,r
c
}包括k
i
个舱体旋转角度,r
a
、r
b
、r
c
分别代表舱体第a、b、c个旋转位置,其中a,b,c≤k且a,b,c∈n*,k
min
≤k
i
≤k,k
min
为舱体的最少旋转变位次数,采用舱体旋转变位性能评价函数rot_fun=(k
a
+k
b
+

+k
c
)/k
i
,其中,k
a
、k
b
、k
c
分别表示第a、b、c个舱体旋转位置下的机器人加工任务的刚度值,选取rot_fun结果最优的方案rot={r
a
,r
b
,

,r
c
}为最终舱体旋转变位方案;步骤6,在选定最优的方案rot={r
a
,r
b
,

,r
c
}的基础上,设定机器人沿舱体轴向方向的移动步长s,设定机器人可选择站位的数量j及机器人站位集合{p1,p2,

,p
j
},p
j
表示第j个机器人站位,遍历所有机器人站位,判断机器人对j个机器人站位下支架位姿的可达性,得到在舱体旋转位置r
a
下机器人站位{p1,p2,

,p
j
}对应的可加工任务集合{x
1a
,x
2a
,

,x
ja
},并采用机器人笛卡尔空间刚度矩阵k=j(θ)-t
k
θ
j(θ)-1
对机器人加工j个机器人站位下的可加工任务集合{x
1a
,x
2a
,

,x
ja
}的加工性能进行评价,得到j个机器人站位下的机器人加工任务的刚度值{k
1a
,k
2a
,

,k
ja
},k
ja
表示第j个机器人站位下的机器人加工任务的刚度值;获得舱体旋转位置r
b
、r
c
下机器人站位{p1,p2,

,p
j
}对应的可加工任务集合{x
1b
,x
2b
,

,x
jb
}、{x
1c
,x
2c
,

,x
jc
},以及得到对应在j个机器人站位下的机器人加工任务的刚度值{k
1b
,k
2b
,

,k
jb
}、{k
1c
,k
2c
,

,k
jc
};
步骤7,根据最优的方案rot={r
a
,r
b
,

,r
c
},以及当前舱体旋转位置r
a
下所有机器人站位{p1,p2,

,p
j
}对应的总加工任务集合x
a
={x
1a
∪x
2a


∪x
ja
},筛选能够加工到当前舱体旋转位置r
a
下的所有任务x
a
的机器人站位组合方案{pos
1a
,pos
2a
,

,pos
naa
},并得到对应机器人站位组合方案的刚度值{k
1a
,k
2a
,

,k
naa
},pos
naa
表示在舱体旋转位置r
a
下的第n
a
个机器人站位组合方案,k
naa
表示在机器人站位组合方案pos
naa
下的机器人加工任务的刚度值,其中在舱体旋转位置r
a
下任意一个机器人站位组合方案pos
ia
={p
d
,p
e
,

,p
f
}包括m个机器人站位,p
d
、p
e
、p
f
分别表示第d、e、f个机器人站位,d,e,f≤j且d,e,f∈n*;获得舱体旋转位置r
b
、r
c
下机器人站位{p1,p2,

,p
j
}对应的总加工任务集合x
b
={x
1b
∪x
2b


∪x
jb
}、x
c
={x
1c
∪x
2c


∪x
jc
},筛选能够加工到当前舱体旋转位置r
b
、r
c
下的所有任务x
b
、x
c
的机器人站位组合方案{pos
1b
,pos
2b
,

,pos
nbb
}、{pos
1c
,pos
2c
,

,pos
ncc
},并得到对应机器人站位组合方案的刚度值{k
1b
,k
2b
,

,k
nbb
}、{k
1c
,k
2c
,

,k
ncc
},pos
nbb
表示在舱体旋转位置r
b
下的第n
b
个机器人站位组合方案,k
ncc
表示在机器人站位组合方案pos
ncc
下的机器人加工任务的刚度值;根据最优的方案rot={r
a
,r
b
,

,r
c
},其中任意一个机器人站位方案pos
i
={pos
ga
,pos
hb
,

,pos
lc
},g,h,l≤j且g,h,l∈n*,pos
ga
代表在舱体旋转位置r
a
时采用机器人站位方案pos
ga
,pos
hb
代表在舱体旋转位置r
b
时采用机器人站位方案pos
hb
,pos
lc
代表在舱体旋转位置r
c
时采用机器人站位方案pos
lc
,计算机器人站位方案{pos
ga
,pos
hb
,

,pos
lc
}下的刚度值{k
ga
,k
hb
,

,k
lc
},采用机器人站位性能评价函数pos_fun=k
ga
+k
hb
+

+k
lc
,选取pos_fun结果最优的方案pos={pos
ga
,pos
hb
,

,pos
lc
}为最终机器人站位方案;步骤8,在最终机器人站位方案pos={pos
ga
,pos
hb
,

,pos
lc
}的基础上,得到机器人r
i
的可加工任务集合x
ri
={x
riag
,x
ribh
,

,x
ricl
},i=1,2,

,m,x
riag
为机器人r
i
在舱体旋转位置r
a
时采用机器人站位方案pos
ga
时的可加工任务集合,x
ribh
为机器人r
i
在舱体旋转位置r
b
时采用机器人站位方案pos
hb
时的可加工任务集合,x
ricl
为机器人r
i
在舱体旋转位置r
c
时采用机器人站位方案pos
lc
时的可加工任务集合,进而对各个机器人的加工任务进行任务分配;步骤9,对加工任务进行分配,得到机器人r
i
的任务分配结果a
ri
={a
riag
,a
ribh
,

,a
ricl
},a
riag
为机器人r
i
在舱体旋转位置r
a
、机器人站位方案pos
ga
下的加工任务分配结果,a
ribh
为机器人r
i
在舱体旋转位置r
b
、机器人站位方案pos
ha
下的加工任务分配结果,a
ricl
为机器人r
i
在舱体旋转位置r
c
、机器人站位方案pos
la
下的加工任务分配结果,优化机器人加工任务的顺序,形成各个机器人的加工任务时序调度表,最终形成多机器人团队的加工调度方案。2.根据权利要求1所述的方法,其特征在于,步骤1中,待加工支架的位姿会随舱体的旋转而发生改变。3.根据权利要求2所述的方法,其特征在于,步骤2中,所述多机器人作业系统单元组成包括机器人型号、agv型号、机器人数量m、末端执行器构造、刀具型号和规格,其中,航天器舱体的尺寸决定了机器人型号、agv型号和机器人数量m,舱体表面待加工支架的工艺要求决定了末端执行器构造、刀具型号和规格。4.根据权利要求3所述的方法,其特征在于,步骤3中,所述多机器人系统的空间布局方案包括机器人在舱体两侧的分布情况、机器人与舱体的相对位置关系;所述机器人在舱体两侧的分布情况影响机器人之间的相对位置关系,进而影响任务分配与任务时序调度环
节;所述机器人与舱体的相对位置关系指机器人距舱体的轴线的水平距离以及舱体轴线距离地面的高度,机器人距离舱体的远近影响机器人的加工姿态,从而影响机器人的加工质量。5.根据权利要求4所述的方法,其特征在于,步骤4中,舱体旋转位置数量k计算公式如下:k=360/α。6.根据权利要求5所述的方法,其特征在于,步骤5中,舱体的最少旋转变位次数k
min
是指在能够覆盖总加工任务集合x
r
中所有任务前提下的最少的舱体旋转变位次数。7.根据权利要求6所述的方法,其特征在于,步骤6中,机器人可选择站位的数量j计算公式如下:j=2l/s。8.根据权利要求7所述的方法,其特征在于,步骤8中,所述对各个机器人的加工任务进行任务分配是指去除重叠任务并将唯一任务分配给对应机器人;重叠任务指可被同侧机器人加工以及能够在不同舱体旋转角度下被加工的任务,唯一任务指只能在一个舱体旋转角度下的一个机器人站位下被加工的任务;将唯一任务分配给对应机器人,计算当前机器人工作量,然后依次将重叠任务分配给能加工重叠任务且工作量最少的机器人,完成多机器人团队的任务分配过程。9.根据权利要求8所述的方法,其特征在于,步骤9中,所述各个机器人的加工任务时序调度表是指各个机器人r
i
加工对应任务集合a
ri
={a
riag
,a
ribh
,

,a
ricl
}中任务的顺序;所述多机器人团队的加工调度方案是指形成以时间为轴,确定在不同时间节点下的舱体旋转角度、机器人加工站位以及在不同站位下的机器人加工任务列表和顺序。10.一种存储介质,其特征在于,存储有计算机程序或指令,当所述计算机程序或指令被运行时,实现如权利要求1至9中任一项所述的方法。

技术总结
本发明提供了面向航天器大型舱体构件加工的多机器人协同规划方法,包括:针对航天器大型舱体构件多机器人协同加工过程,首先规划舱体旋转变位次数与每次旋转的角度,其次规划多机器人站位布局与换站策略,最后在舱体位置和机器人站位固定不变情况下规划多机器人加工任务时序。本发明可以快速系统地对大型舱体构件加工过程进行规划,合理选择舱体旋转变位和机器人站位方案,提高机器人协同作业刚度;合理规划加工时序使得多机器人加工任务过程更加紧凑,可缩短机器人空闲时间,提高多机器人协同加工效率。人协同加工效率。人协同加工效率。


技术研发人员:李波 林佳美 田威 廖文和 刘少睿
受保护的技术使用者:南京航空航天大学
技术研发日:2023.05.17
技术公布日:2023/8/9
版权声明

本文仅代表作者观点,不代表航家之家立场。
本文系作者授权航家号发表,未经原创作者书面授权,任何单位或个人不得引用、复制、转载、摘编、链接或以其他任何方式复制发表。任何单位或个人在获得书面授权使用航空之家内容时,须注明作者及来源 “航空之家”。如非法使用航空之家的部分或全部内容的,航空之家将依法追究其法律责任。(航空之家官方QQ:2926969996)

航空之家 https://www.aerohome.com.cn/

飞机超市 https://mall.aerohome.com.cn/

航空资讯 https://news.aerohome.com.cn/

分享:

扫一扫在手机阅读、分享本文

相关推荐