基于Pareto熵粒子群算法的机器人性能优化方法技术

本发明专利技术基于Pareto熵粒子群算法的机器人性能优化方法，通过对一种串联机器人进行运动学分析，得到了机器人的多目标优化模型。将机器人运动灵活性指标条件数的全域平均性能指标、全域波动性能指标、全域最差灵活性指标作为三个优化目标。利用多目标粒子群算法，引入Pareto效率对它们进行优化。将优化后的Pareto第一前端最优解集进行数据分析，最终经优化计算得到了可行的杆长参数，实现了机器人的机构优化，使机器人性能有了较大提高。

【技术实现步骤摘要】
基于Pareto熵粒子群算法的机器人性能优化方法
本专利技术涉及机器人性能优化方法领域，特别是涉及基于Pareto熵粒子群算法的机器人性能优化方法。
技术介绍
机器人的机构设计是一个困难而复杂的问题，而性能评价指标是设计机器人的关键问题之一。在设计过程中不可避免地要以一些性能指标作为基础，包括机器人的条件数、可操作度、最小奇异值等一些指标。在机器人性能指标的研究上，国内外许多学者做了大量的研究工作，并取得了许多成果。YOSHIKAWA定义了机器人的可操作度，并利用可操作度椭球对其几何意义进行了阐述。SALISBURY等提出了条件数，并用其对机器手手指工作空间进行了评价。但是雅可比矩阵条件数依赖于机器人的位姿，只是一个局部指标，指标的优化往往具有局限性。因此在对机器人运动学多个性能评价指标之间的耦合和制约问题的研究基础上，将运动学参数进行多目标性能优化，对机器人性能的提高具有重要意义。一些学者已经在多指标优化问题上做了相关研究工作。ZhangLufan采用NSGA-II算法优化了宏微观运动平台连接机构的总质量、关键点最小变形以及第1阶固有频率。ZhangDan利用粒子群优化算法对土拨鼠机器人的刚度和工作空间进行了优化。粒子群优化算法因形式简洁、收敛快速和参数调节机制灵活等优点，同时一次运行可得到多个解，且能逼近非凸或不连续的Pareto最优前端，因而被认为是求解多目标优化问题最具潜力的方法之一。
技术实现思路
为了解决上述存在的问题，本专利技术提供一种基于Pareto熵粒子群算法的机器人性能优化方法，采用改进的多目标粒子群的优化方法进行一种焊接机器人的运动学性能多目标优化。通过对机器人的运动学分析，建立其运动学性能的多目标优化模型。以机器人运动灵活性能指标条件数的全域运动平均性能、全域运动波动性能、全域运动最差性能3个指标为优化目标，以机器人的杆长参数为优化变量进行多目标优化，比较优化结果与初始杆长时机器人的性能，进而对机器人性能优化进行改进，为达此目的，本专利技术提供基于Pareto熵粒子群算法的机器人性能优化方法,其特征在于：1)建立机器人运动模型；机器人为六自由度机器人采用，由6个转动副构成，按照D-H标准建立坐标系；根据机构影响系数法建立一阶影响系数矩阵，Si表示连杆i轴线方向的单位矢量，Ri表示连杆i坐标系与固定坐标系两原点之间的距离，(P-Ri)表示连杆i对研究点P的矢径，则机器人一阶影响系数雅可比矩阵J为：2)设定优化目标；雅可比矩阵的条件数代表了雅可比转换矩阵向各个方向的变化均一性，在设计机器人机构和对机器人进行控制时，要使雅可比矩阵尽量各向同性，即尽量“均匀”，就要控制雅可比矩阵的条件数尽量接近1，为此需设定了一个指标η：η表示全域性能指标，它等于机器人工作空间w内所有观测点的局部指标条件数k倒数的平均值，在η的基础上又设定了全域性能波动指标σ：以上两个指标可以提供机器人运动灵活性的；增加一个描述最差灵活性指标的δ：δ＝η-kmin(4)；它由表示整体灵活性的指标η与表示最差灵活性指标kmin的差来评价，δ越小，说明在任务空间上灵活性变化幅度越小；3)问题描述；采用目标函数的最小方向寻优，对三个指标η、σ与δ进行变化得到目标函数f1、f2、f3进行多目标优化，最后找到杆长l1和l2一些最优的组合；可用数学表示如下：4)优化方法；其算法流程如下：输入：具有M个目标的待优化问题MOP，具有D个决策变量的搜索空间SD；初始化参数：最大外部档案容量F，种群粒子个数PN，最大迭代次数T_max；输出：存储在外部档案中的近似Pareto最优解集Gpa；步骤1：初始化：产生初始化种群P，令迭代计数器n＝0，在搜索空间SD中按均匀分布随机生成PN个粒子，计算每个粒子的M个目标函数值，初始化每个粒子的外部档案Pa，更新Gpa；步骤2：更新迭代计数器：t＝t+1；步骤3：评估进化环境：计算G_pa的Pareto熵和差熵，评估种群进化状态Status，计算G_pa中每个Pareto最优解的个体格距离密度和个体格占优强度；步骤4：自适应调节粒子运动参数,计算当前粒子运动参数；步骤5：更新种群：选择粒子i的全局最优解；gBest(i),从Pa(i)中选择一个与gBest(i)空间距；离最近的成员作为粒子i的个体最优解pBest(i)；更新粒子i的速度和位置，评估粒子i的M个目标函；数值，更新粒子i的个体外部档案Pa(i)，更新全局；外部档案Gpa：步骤6：检测算法终止条件：如果n>T_max，则输出Gpa；否则，跳转步骤(2)。作为本专利技术进一步运用，本专利技术所述六自由度机器人为弧焊机器人或涂胶机器人或切割机器人或喷漆机器人。本专利技术基于Pareto熵粒子群算法的机器人性能优化方法，通过对一种串联机器人进行运动学分析，得到了机器人的多目标优化模型。将机器人运动灵活性指标条件数的全域平均性能指标、全域波动性能指标、全域最差灵活性指标作为三个优化目标。利用多目标粒子群算法，引入Pareto效率对它们进行优化。将优化后的Pareto第一前端最优解集进行数据分析，最终经优化计算得到了可行的杆长参数，实现了机器人的机构优化，使机器人性能有了较大提高。附图说明图1是本专利技术机器人连杆参数示意图；图2是本专利技术算法流程示意图；图3是本专利技术优化指标1的等高线分布图；图4是本专利技术优化指标2的等高线分布图；图5是本专利技术优化指标3的等高线分布图；图6是本专利技术优化结果示意图；图7是本专利技术优化结果全域性能均值比较图；图8是本专利技术优化结果全域性能波动值比较图；图9是本专利技术优化结果全域性能最差值比较图；具体实施方式下面结合附图与具体实施方式对本专利技术作进一步详细描述：本专利技术提供一种基于Pareto熵粒子群算法的机器人性能优化方法，采用改进的多目标粒子群的优化方法进行一种焊接机器人的运动学性能多目标优化。通过对机器人的运动学分析，建立其运动学性能的多目标优化模型。以机器人运动灵活性能指标条件数的全域运动平均性能、全域运动波动性能、全域运动最差性能3个指标为优化目标，以机器人的杆长参数为优化变量进行多目标优化，比较优化结果与初始杆长时机器人的性能，进而对机器人性能优化进行改进。1)建立机器人运动模型；本专利技术机器人为六自由度机器人，由6个转动副构成，该机器人是一种典型的串联机器人。由于它的灵活性，可以广泛应用于弧焊、涂胶、切割、喷漆等多种领域。按照D-H标准建立坐标系，可得出该机器人的D-H参数如图1。为了分析机器人的运动学，根据机构影响系数法建立一阶影响系数矩阵，Si表示连杆i轴线方向的单位矢量，Ri表示连杆i坐标系与固定坐标系两原点之间的距离，(P-Ri)表示连杆i对研究点P的矢径，则机器人一阶影响系数雅可比矩阵J为：2)设定优化目标；雅可比矩阵的条件数代表了雅可比转换矩阵向各个方向的变化均一性，在设计机器人机构和对机器人进行控制时，要使雅可比矩阵尽量各向同性，即尽量“均匀”，就要控制雅可比矩阵的条件数尽量接近1，由于该指标只是用于表示机器人的局部的一个性能值，并不能反映工作空间的整体性能指标，为此需定义了一个新指标η：η表示全域性能指标，它等于机器人工作空间w内所有观测点的局部指标条件数k倒数的平均值，在η的基础上又提出了全域性能波动指标σ本文档来自技高网...

【技术保护点】
基于Pareto熵粒子群算法的机器人性能优化方法,其特征在于：1)建立机器人运动模型；机器人为六自由度机器人采用，由6个转动副构成，按照D‑H标准建立坐标系；根据机构影响系数法建立一阶影响系数矩阵，Si表示连杆i轴线方向的单位矢量，Ri表示坐标系i与固定坐标系两原点之间的距离，(P‑Ri)表示连杆i对研究点P的矢径，则机器人一阶影响系数雅可比矩阵J为：2)设定优化目标；雅可比矩阵的条件数代表了雅可比转换矩阵向各个方向的变化均一性，在设计机器人机构和对机器人进行控制时，要使雅可比矩阵尽量各向同性，即尽量“均匀”，就要控制雅可比矩阵的条件数尽量接近1，为此需设定了一个指标η：η表示全域性能指标，它等于机器人工作空间w内所有观测点的局部指标条件数k倒数的平均值，在η的基础上设定全域性能波动指标σ：以上两个指标可以提供机器人运动灵活性的；增加一个描述最差灵活性指标的δ：δ＝η‑kmin      (4)；它由表示整体灵活性的指标η与表示最差灵活性指标kmin的差来评价，δ越小，说明在任务空间上灵活性变化幅度越小；3)问题描述；采用目标函数的最小方向寻优，对三个指标η、σ与δ进行变化得到目标函数f1、f2、f3进行多目标优化，最后找到杆长l1和l2一些最优的组合；可用数学表示如下：4)优化方法；其算法流程如下：输入：具有M个目标的待优化问题MOP，具有D个决策变量的搜索空间SD；初始化参数：最大外部档案容量F，种群粒子个数PN，最大迭代次数T_max；输出：存储在外部档案中的近似Pareto最优解集G_pa；步骤1：初始化：产生初始化种群P，令迭代计数器n＝0，在搜索空间SD中按均匀分布随机生成PN个粒子，计算每个粒子的M个目标函数值，初始化每个粒子的外部档案Pa，更新Gpa；步骤2：更新迭代计数器：t＝t+1；步骤3：评估进化环境：计算G_pa的Pareto熵和差熵，评估种群进化状态Status，计算G_pa中每个Pareto最优解的个体格距离密度和个体格占优强度；步骤4：自适应调节粒子运动参数,计算当前粒子运动参数；步骤5：更新种群：选择粒子i的全局最优解；gBest(i),从Pa(i)中选择一个与gBest(i)空间距；离最近的成员作为粒子i的个体最优解pBest(i)；更新粒子i的速度和位置，评估粒子i的M个目标函；数值，更新粒子i的个体外部档案Pa(i)，更新全局；外部档案G_pa：步骤6：检测算法终止条件：如果n>T_max，则输出G_pa；否则，跳转步骤(2)。...

【技术特征摘要】
1.基于Pareto熵粒子群算法的机器人性能优化方法,其特征在于：1)建立机器人运动模型；机器人为六自由度机器人采用，由6个转动副构成，按照D-H标准建立坐标系；根据机构影响系数法建立一阶影响系数矩阵，Si表示连杆i轴线方向的单位矢量，Ri表示连杆i坐标系与固定坐标系两原点之间的距离，(P-Ri)表示连杆i对研究点P的矢径，则机器人一阶影响系数雅可比矩阵J为：2)设定优化目标；雅可比矩阵的条件数代表了雅可比转换矩阵向各个方向的变化均一性，在设计机器人机构和对机器人进行控制时，要使雅可比矩阵尽量各向同性，即尽量“均匀”，就要控制雅可比矩阵的条件数尽量接近1，为此需设定了一个指标η：η表示全域性能指标，它等于机器人工作空间w内所有观测点的局部指标条件数k倒数的平均值，在η的基础上设定全域性能波动指标σ：以上两个指标可以提供机器人运动灵活性的；增加一个描述最差灵活性指标的δ：δ＝η-kmin(4)；它由表示整体灵活性的指标η与表示最差灵活性指标kmin的差来评价，δ越小，说明在任务空间上灵活性变化幅度越小；3)问题描述；采用目标函数的最小方向寻优，对三个指标η...

【专利技术属性】
技术研发人员：应明峰，莫晓晖，姜劲，张静，吴敏，
申请(专利权)人：金陵科技学院，
类型：发明
国别省市：江苏;32