CN105479459A - 机器人零力控制方法以及系统 - Google Patents

Abstract

一种机器人零力控制方法以及系统，方法包括：S1、基于惯性力、离心力、科里奥利力、粘摩擦力、静摩擦力、重力，建立机器人动力学模型；S2、基于所述机器人动力学模型以及反馈数据，计算机器人各个关节的外力力矩；S3、根据步骤S2中计算的外力力矩计算速度指令；S4、根据速度指令以及位置反馈计算位置指令。本发明通过建立动力学模型，可以直接计算各关节的外力力矩，无需借助力/转矩传感器，降低了系统成本和复杂性；采用位置指令控制的位置控制方式，而非采用直接转矩控制的方式，降低了系统安全性与稳定性的设计难度，另外，由于在计算外力力矩时考虑了惯性力，该发明还可适用于自重较大的机器人。

Description

机器人零力控制方法以及系统

技术领域

本发明涉及零力控制领域，尤其涉及一种机器人零力控制方法以及系统。

背景技术

目前工业机器人的示教过程主要依赖于示教器，这样的示教方式具有工作效率低、过程繁琐不直观、对操作人员知识水平要求高的特点。因此，有必要设计一种直观高效的机器人直接示教系统。直接示教系统的开发依赖于零力控制技术的实现。当前零力控制技术主要有两种实现方式：基于位置控制的零力控制技术和基于直接转矩控制的零力控制技术。

基于位置控制的零力控制技术需借助于力/转矩传感器，测量出施加在机器人各关节的外力，然后将该外力与经动力学模型计算的力相加并转换为相应的位置指令。由于该方案需要力/转矩传感器，所以会极大提高系统的成本和复杂性。

基于直接转矩控制的零力控制技术，首先通过动力学模型计算出各关节对应的重力及摩擦力，然后控制各关节输出对应大小的转矩，机器人在外力的作用下克服惯性力而运动。该方案不需要知道精确的加速度信息，但是由于其直接对力矩进行控制，对系统的稳定性设计提出了较大的挑战，同时，由于这样的实现方式不对惯性力做补偿，无法适用于自重较大的机器人。

发明内容

本发明要解决的技术问题在于，针对现有技术的上述缺陷，提供一种机器人零力控制方法以及系统。

本发明解决其技术问题所采用的技术方案是：构造一种机器人零力控制方法，包括：

S1、基于惯性力、离心力、科里奥利力、粘摩擦力、静摩擦力、重力，建立机器人动力学模型；

S2、基于所述机器人动力学模型以及反馈数据，计算机器人各个关节的外力力矩；

S3、根据步骤S2中计算的外力力矩计算速度指令；

S4、根据速度指令以及位置反馈计算位置指令。

在本发明所述的机器人零力控制方法中，所述动力学模型为：

$B\left(q\right)\stackrel{\·\·}{q}+C_1\left(q\right)\stackrel{\·}{q}+C_2(q,\stackrel{\·}{q})\stackrel{\·}{q}+F_v\stackrel{\·}{q}+F_{s}sign\left(\stackrel{\·}{q}\right)+G\left(q\right)=\mathrm{\τ}---\left(1\right)$

其中，B(q)为N×N惯量矩阵，C₁(q)为N×N离心力矩阵，为N×N科里奥利力矩阵，F_v为N×N粘摩擦力矩阵，F_s为N×N静摩擦力矩阵，G(q)为N×1重力向量，τ为N×1关节力矩向量，表示以为自变量的符号函数，为加速度信息，为速度信息，q为位置信息，N为机器人关节数。

在本发明所述的机器人零力控制方法中，所述步骤S2包括：

S21、将公式(1)中静摩擦力设为0，并将动力学模型中的设为期望加速度，设为速度反馈q设为位置反馈q_f，得到如公式(2)所示的关节作用力力矩τ_nsf：

$B\left(q_f\right){\stackrel{\·\·}{q}}_r+C_1\left(q_f\right){\stackrel{\·}{q}}_f+C_2(q_f,\stackrel{\·}{q_f}){\stackrel{\·}{q}}_f+F_v{\stackrel{\·}{q}}_f+G\left(q_f\right)={\mathrm{\τ}}_{nsf}---\left(1\right);$

S22、根据力矩常数以及反馈电流计算反馈力矩τ_f：K_tI_f＝τ_f；

S23、根据步骤S21中计算的关节作用力力矩τ_nsf以及步骤S22中计算的反馈力矩τ_f计算得到外力力矩τ_e：τ_f-τ_nsf＝τ_e。

在本发明所述的机器人零力控制方法中，所述步骤S3包括：

S31、将公式(1)中静摩擦力以外的作用力设为0，并将动力学模型中的设为速度反馈得到静摩擦力力矩τ_sf： $F_{s}sign\left(\stackrel{\·}{q_f}\right)={\mathrm{\τ}}_{sf};$

S32、根据所述步骤S41中的静摩擦力力矩τ_sf以及步骤S21中的关节作用力力矩τ_nsf计算外力力矩死区作用域[τ_dn,τ_dp]，其中，τ_dn,τ_dp基于以下公式确定：

$\left\{\begin{array}{c}{k}_1{\mathrm{\τ}}_{nsf}-\left|{\mathrm{\τ}}_{sf}\right|={\mathrm{\τ}}_{dn}\\ k_2{\mathrm{\τ}}_{nsf}+\left|{\mathrm{\τ}}_{sf}\right|={\mathrm{\τ}}_{dp}\end{array}\right.$

其中，k₁、k₂表示力矩误差校正系数；

S33、根据步骤S23中计算得到的外力力矩τ_e，基于以下公式计算得到速度指令

$\left\{\begin{array}{cc}{k}_{3}sign\left({\mathrm{\τ}}_e\right){\left|{\mathrm{\τ}}_e\right|}^p={\stackrel{\·}{q}}_r,& {\mathrm{\τ}}_e\∉\[{\mathrm{\τ}}_{dn},{\mathrm{\τ}}_{dp}\]\\ 0={\stackrel{\·}{q}}_r,& {\mathrm{\τ}}_e\∈\[{\mathrm{\τ}}_{dn},{\mathrm{\τ}}_{dp}\]\end{array}\right.$

其中，k₃为力矩速度映射系数，p∈(0,1)。

在本发明所述的机器人零力控制方法中，所述步骤S4基于以下公式计算得到位置指令q_r；

$q_f+{\stackrel{\·}{\hat{q}}}_r\mathrm{\Δ}t=q_r$

其中，Δt为位置插补周期，q_f为位置反馈，为速度指令。

本发明还公开了一种机器人零力控制系统，包括：

动力学模型模块，用于基于惯性力、离心力、科里奥利力、粘摩擦力、静摩擦力、重力，建立机器人动力学模型；

外力力矩计算模块，用于基于所述动力学模型模块中的机器人动力学模型以及反馈数据，计算机器人各个关节的外力力矩；

力矩速度映射模块，用于根据外力力矩计算模块中计算的外力力矩计算速度指令；

速度位置映射模块，用于根据速度指令以及位置反馈计算位置指令。

在本发明所述的机器人零力控制系统中，所述动力学模型为：

$B\left(q\right)\stackrel{\·\·}{q}+C_1\left(q\right)\stackrel{\·}{q}+C_2(q,\stackrel{\·}{q})\stackrel{\·}{q}+F_v\stackrel{\·}{q}+F_{s}sign\left(\stackrel{\·}{q}\right)+G\left(q\right)=\mathrm{\τ}---\left(1\right)$

其中，B(q)为N×N惯量矩阵，C₁(q)为N×N离心力矩阵，为N×N科里奥利力矩阵，F_v为N×N粘摩擦力矩阵，F_s为N×N静摩擦力矩阵，G(q)为N×1重力向量，τ为N×1关节力矩向量，表示以为自变量的符号函数，为加速度信息，为速度信息，q为位置信息，N为机器人关节数。

在本发明所述的机器人零力控制系统中，所述外力力矩计算模块包括：

关节作用力力矩计算单元，用于将公式(1)中静摩擦力设为0，并将动力学模型中的设为期望加速度，设为速度反馈q设为位置反馈q_f，得到如公式(2)所示的关节作用力力矩τ_nsf：

$B\left(q_f\right){\stackrel{\·\·}{q}}_r+C_1\left(q_f\right){\stackrel{\·}{q}}_f+C_2(q_f,\stackrel{\·}{q_f}){\stackrel{\·}{q}}_f+F_v{\stackrel{\·}{q}}_f+G\left(q_f\right)={\mathrm{\τ}}_{nsf}---\left(1\right);$

反馈力矩计算单元，用于根据力矩常数以及反馈电流计算反馈力矩τ_f：K_tI_f＝τ_f；

外力力矩计算单元，用于根据关节作用力力矩计算单元计算的关节作用力力矩τ_nsf以及反馈力矩计算单元计算的反馈力矩τ_f计算得到外力力矩τ_e：τ_f-τ_nsf＝τ_e。

在本发明所述的机器人零力控制系统中，所述力矩速度映射模块包括：

静摩擦力力矩计算单元，用于将公式(1)中静摩擦力以外的作用力设为0，并将动力学模型中的设为速度反馈得到静摩擦力力矩τ_sf：

外力力矩死区作用域计算单元，用于根据所述静摩擦力力矩计算单元计算的静摩擦力力矩τ_sf以及关节作用力力矩计算单元计算的关节作用力力矩τ_nsf计算外力力矩死区作用域[τ_dn,τ_dp]，其中，τ_dn,τ_dp基于以下公式确定：

$\left\{\begin{array}{c}{k}_1{\mathrm{\τ}}_{nsf}-\left|{\mathrm{\τ}}_{sf}\right|={\mathrm{\τ}}_{dn}\\ k_2{\mathrm{\τ}}_{nsf}+\left|{\mathrm{\τ}}_{sf}\right|={\mathrm{\τ}}_{dp}\end{array}\right.$

其中，k₁、k₂表示力矩误差校正系数；

速度指令计算单元，用于根据外力力矩计算单元中计算得到的外力力矩τ_e，基于以下公式计算得到速度指令

$\left\{\begin{array}{cc}{k}_{3}sign\left({\mathrm{\τ}}_e\right){\left|{\mathrm{\τ}}_e\right|}^p={\stackrel{\·}{q}}_r,& {\mathrm{\τ}}_e\∉\[{\mathrm{\τ}}_{dn},{\mathrm{\τ}}_{dp}\]\\ 0={\stackrel{\·}{q}}_r,& {\mathrm{\τ}}_e\∈\[{\mathrm{\τ}}_{dn},{\mathrm{\τ}}_{dp}\]\end{array}\right.$

其中，k₃为力矩速度映射系数，p∈(0,1)。

在本发明所述的机器人零力控制系统中，所述的根据速度指令以及位置反馈计算位置指令包括：基于以下公式计算得到位置指令q_r；

$q_f+{\stackrel{\·}{\hat{q}}}_r\mathrm{\Δ}t=q_r$

其中，Δt为位置插补周期，q_f为位置反馈，为速度指令。

实施本发明的机器人零力控制方法以及系统，具有以下有益效果：本发明通过建立动力学模型，可以直接计算各关节的外力力矩，无需借助力/转矩传感器，降低了系统成本和复杂性；采用位置指令控制的位置控制方式，而非采用直接转矩控制的方式，降低了系统安全性与稳定性的设计难度，另外，由于在计算外力力矩时考虑了惯性力，该发明还可适用于自重较大的机器人。

附图说明

下面将结合附图及实施例对本发明作进一步说明，附图中：

图1是本发明的机器人零力控制方法的流程图；

图2是本发明的机器人零力控制系统的结构示意图；

图3是速度反馈与计算得到的外力力矩之间关系示意图。

具体实施方式

为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。

如图1所示，是本发明的机器人零力控制方法的流程图；

本发明的机器人零力控制方法，包括：

S1、基于惯性力、离心力、科里奥利力、粘摩擦力、静摩擦力、重力，建立机器人动力学模型；

S2、基于所述机器人动力学模型以及反馈数据，计算机器人各个关节的外力力矩；

S3、根据步骤S2中计算的外力力矩计算速度指令；

S4、根据速度指令以及位置反馈计算位置指令。

其中，步骤S1中的所述动力学模型为：

$B\left(q\right)\stackrel{\·\·}{q}+C_1\left(q\right)\stackrel{\·}{q}+C_2(q,\stackrel{\·}{q})\stackrel{\·}{q}+F_v\stackrel{\·}{q}+F_{s}sign\left(\stackrel{\·}{q}\right)+G\left(q\right)=\mathrm{\τ}---\left(1\right)$

其中，B(q)为N×N惯量矩阵，C₁(q)为N×N离心力矩阵，为N×N科里奥利力矩阵，F_v为N×N粘摩擦力矩阵，F_s为N×N静摩擦力矩阵，G(q)为N×1重力向量，τ为N×1关节力矩向量，表示以为自变量的符号函数，为加速度信息，为速度信息，q为位置信息，N为机器人关节数。

在公式(1)的基础上，所述步骤S2包括：

S21、将公式(1)中静摩擦力设为0，并将动力学模型中的设为期望加速度，设为速度反馈q设为位置反馈q_f，得到如公式(2)所示的关节作用力力矩τ_nsf：

$B\left(q_f\right){\stackrel{\·\·}{q}}_r+C_1\left(q_f\right){\stackrel{\·}{q}}_f+C_2(q_f,\stackrel{\·}{q_f}){\stackrel{\·}{q}}_f+F_v{\stackrel{\·}{q}}_f+G\left(q_f\right)={\mathrm{\τ}}_{nsf}---\left(1\right);$

S22、根据力矩常数以及反馈电流计算反馈力矩τ_f：K_tI_f＝τ_f；

S23、根据步骤S21中计算的关节作用力力矩τ_nsf以及步骤S22中计算的反馈力矩τ_f计算得到外力力矩τ_e：τ_f-τ_nsf＝τ_e。

因此，通过建立动力学模型，可以直接计算各关节的外力力矩，无需借助力/转矩传感器，降低了系统成本和复杂性。

在公式(1)的基础上，所述步骤S3包括：

S31、将公式(1)中静摩擦力以外的作用力设为0，并将动力学模型中的设为速度反馈得到静摩擦力力矩τ_sf：因为真正使得机器人运动的外力力矩不包括静摩擦力，所以在估算外力力矩时需要将静摩擦力去掉。

S32、根据所述步骤S41中的静摩擦力力矩τ_sf以及步骤S21中的关节作用力力矩τ_nsf计算外力力矩死区作用域[τ_dn,τ_dp]，其中，τ_dn,τ_dp基于以下公式确定：

$\left\{\begin{array}{c}{k}_1{\mathrm{\τ}}_{nsf}-\left|{\mathrm{\τ}}_{sf}\right|={\mathrm{\τ}}_{dn}\\ k_2{\mathrm{\τ}}_{nsf}+\left|{\mathrm{\τ}}_{sf}\right|={\mathrm{\τ}}_{dp}\end{array}\right.$

其中，k₁、k₂表示力矩误差校正系数；

S33、根据步骤S23中计算得到的外力力矩τ_e，基于以下公式计算得到速度指令

$\left\{\begin{array}{cc}{k}_{3}sign\left({\mathrm{\τ}}_e\right){\left|{\mathrm{\τ}}_e\right|}^p={\stackrel{\·}{q}}_r,& {\mathrm{\τ}}_e\∉\[{\mathrm{\τ}}_{dn},{\mathrm{\τ}}_{dp}\]\\ 0={\stackrel{\·}{q}}_r,& {\mathrm{\τ}}_e\∈\[{\mathrm{\τ}}_{dn},{\mathrm{\τ}}_{dp}\]\end{array}\right.$

其中，k₃为力矩速度映射系数，p∈(0,1)，因此当外力力矩在0附近时，速度变化较快，当外力力矩较大时，嘟嘟变化较慢，也就说机器人反应变迟钝，更安全。此公式将力与速度在其他领域中比较常见的关系引入本发明。

具体的，所述步骤S4包括：基于以下公式计算得到位置指令q_r；

$q_f+{\stackrel{\·}{\hat{q}}}_r\mathrm{\Δ}t=q_r$

其中，Δt为位置插补周期，q_f为位置反馈，为速度指令。

因此本发明采用位置指令的位置控制方式，而非采用直接转矩控制的方式，降低了系统安全性与稳定性的设计难度，同时，由于公式(1)中的机器人动力学模型之中考虑了惯性力，可以适用于自重较大的机器人。

相应的本发明公开了一种机器人零力控制系统，包括：

动力学模型模块，用于基于惯性力、离心力、科里奥利力、粘摩擦力、静摩擦力、重力，建立机器人动力学模型；

外力力矩计算模块，用于基于所述动力学模型模块中的机器人动力学模型以及反馈数据，计算机器人各个关节的外力力矩；

力矩速度映射模块，用于根据外力力矩计算模块中计算的外力力矩计算速度指令；

速度位置映射模块，用于根据速度指令以及位置反馈计算位置指令。

其中，所述动力学模型为：

$B\left(q\right)\stackrel{\·\·}{q}+C_1\left(q\right)\stackrel{\·}{q}+C_2(q,\stackrel{\·}{q})\stackrel{\·}{q}+F_v\stackrel{\·}{q}+F_{s}sign\left(\stackrel{\·}{q}\right)+G\left(q\right)=\mathrm{\τ}---\left(1\right)$

其中，B(q)为N×N惯量矩阵，C₁(q)为N×N离心力矩阵，为N×N科里奥利力矩阵，F_v为N×N粘摩擦力矩阵，F_s为N×N静摩擦力矩阵，G(q)为N×1重力向量，τ为N×1关节力矩向量，表示以为自变量的符号函数，为加速度信息，为速度信息，q为位置信息，N为机器人关节数。

其中，所述外力力矩计算模块包括：

关节作用力力矩计算单元，用于将公式(1)中静摩擦力设为0，并将动力学模型中的设为期望加速度，设为速度反馈q设为位置反馈q_f，得到如公式(2)所示的关节作用力力矩τ_nsf：

$B\left(q_f\right){\stackrel{\·\·}{q}}_r+C_1\left(q_f\right){\stackrel{\·}{q}}_f+C_2(q_f,\stackrel{\·}{q_f}){\stackrel{\·}{q}}_f+F_v{\stackrel{\·}{q}}_f+G\left(q_f\right)={\mathrm{\τ}}_{nsf}---\left(1\right);$

反馈力矩计算单元，用于根据力矩常数以及反馈电流计算反馈力矩τ_f：K_tI_f＝τ_f；

外力力矩计算单元，用于根据关节作用力力矩计算单元计算的关节作用力力矩τ_nsf以及反馈力矩计算单元计算的反馈力矩τ_f计算得到外力力矩τ_e：τ_f-τ_nsf＝τ_e。

其中，所述力矩速度映射模块包括：

静摩擦力力矩计算单元，用于将公式(1)中静摩擦力以外的作用力设为0，并将动力学模型中的设为速度反馈得到静摩擦力力矩τ_sf：

外力力矩死区作用域计算单元，用于根据所述静摩擦力力矩计算单元计算的静摩擦力力矩τ_sf以及关节作用力力矩计算单元计算的关节作用力力矩τ_nsf计算外力力矩死区作用域[τ_dn,τ_dp]，其中，τ_dn,τ_dp基于以下公式确定：

$\left\{\begin{array}{c}{k}_1{\mathrm{\τ}}_{nsf}-\left|{\mathrm{\τ}}_{sf}\right|={\mathrm{\τ}}_{dn}\\ k_2{\mathrm{\τ}}_{nsf}+\left|{\mathrm{\τ}}_{sf}\right|={\mathrm{\τ}}_{dp}\end{array}\right.$

其中，k₁、k₂表示力矩误差校正系数；

速度指令计算单元，用于根据外力力矩计算单元中计算得到的外力力矩τ_e，基于以下公式计算得到速度指令

$\left\{\begin{array}{cc}{k}_{3}sign\left({\mathrm{\τ}}_e\right){\left|{\mathrm{\τ}}_e\right|}^p={\stackrel{\·}{q}}_r,& {\mathrm{\τ}}_e\∉\[{\mathrm{\τ}}_{dn},{\mathrm{\τ}}_{dp}\]\\ 0={\stackrel{\·}{q}}_r,& {\mathrm{\τ}}_e\∈\[{\mathrm{\τ}}_{dn},{\mathrm{\τ}}_{dp}\]\end{array}\right.$

其中，k₃为力矩速度映射系数，p∈(0,1)。

其中，所述的根据速度指令以及位置反馈计算位置指令包括：基于以下公式计算得到位置指令q_r；

$q_f+{\stackrel{\·}{\hat{q}}}_r\mathrm{\Δ}t=q_r$

其中，Δt为位置插补周期，q_f为位置反馈，为速度指令。

下面结合一个具体的实施例详细说明本发明。

现以某SCARA机器人第一、二关节为零力控制实施对象，首先根据步骤S1建立动力学模型如下

$\begin{array}{c}(I_{l_1}+m_{l_1}{l}_1^2+I_{l_2}+m_{l_2}(a_1^2+l_2^2+2a_1{l}_2{c}_2\left)\right)+\stackrel{\·\·}{q_1}+(I_{l_2}+m_{l_2}(l_2^2+a_1{l}_2{c}_2\left)\right)\stackrel{\·\·}{q_2}\\ -2m_{l_2}{a}_1{l}_2{s}_2\stackrel{\·}{q_1}\stackrel{\·}{q_2}-m_{l_2}{a}_1{l}_2{s}_2{\stackrel{\·}{q_2}}^2+f_{v1}\stackrel{\·}{q_1}+f_{s1}sign\left(\stackrel{\·}{q_1}\right)={\mathrm{\τ}}_1,\end{array}$

$\begin{array}{c}(I_{l_2}+m_{l_2}(l_2^2+a_1{l}_2{c}_2\left)\right)+\stackrel{\·\·}{q_1}+(I_{l_2}+m_{l_2}{l}_2^2)\stackrel{\·\·}{q_2}\\ +m_{l_2}{a}_1{l}_2{s}_2{\stackrel{\·}{q_1}}^2+f_{v2}\stackrel{\·}{q_2}+f_{s2}sign\left(\stackrel{\·}{q_2}\right)={\mathrm{\τ}}_2\end{array}$

其中，为第1，2连杆绕质心坐标系Z轴的转动惯量，为第1，2连杆质量(含电机等转配部件)，a₁，a₂为第1，2连杆长度，l₁，l₂为第1，2连杆质心距离连杆起始处的长度，f_v1、f_v2、f_s1、f_s2分别为对应第一个连杆的粘摩擦和静摩擦系数，q₁，q₂，为第1，2连杆旋转的角度，速度及加速度，s_i，c_i为第i个连杆位置的正余弦运算，τ₁，τ₂为第1，2连杆的力矩。

然后按照步骤S21-S23计算得出外力力矩，如图3中实线曲线所示；根据上述步骤S31-S32并取k₁＝k₂＝0.05，得出外力力矩作用于死区，根据步骤S33，分别针对第一二关节取k₃＝0.4106，p＝0.5及k₃＝0.6584，p＝0.5，根据步骤S3，取Δt＝0.0005s，映射外力力矩为速度指令，在通过步骤S4转换为位置指令发送至伺服，驱动机器人运动。

为验证本发明的正确性，通过采集机器人速度反馈，如图3虚线曲线所示。通过图3可以看出，当外力力矩一直明显偏离0值时，关节速度呈明显增加或明显减小的趋势，一旦外力力矩在0值附近波动时，对应的反馈速度也有相应的波动，这表明本发明的正确性在实际机器人之上得到了验证。

综上所述，实施本发明的机器人零力控制方法以及系统，具有以下有益效果：本发明通过建立动力学模型，可以直接计算各关节的外力力矩，无需借助力/转矩传感器，降低了系统成本和复杂性；采用位置控制的方式，而非采用直接转矩控制的方式，降低了系统安全性与稳定性的设计难度，另外，由于在计算外力力矩时考虑了惯性力，该发明还可适用于自重较大的机器人。

上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。

Claims (10)

1.一种机器人零力控制方法，其特征在于，包括：

S1、基于惯性力、离心力、科里奥利力、粘摩擦力、静摩擦力、重力，建立机器人动力学模型；

S2、基于所述机器人动力学模型以及反馈数据，计算机器人各个关节的外力力矩；

S3、根据步骤S2中计算的外力力矩计算速度指令；

S4、根据速度指令以及位置反馈计算位置指令。

2.根据权利要求1所述的机器人零力控制方法，其特征在于，所述动力学模型为：

$B\left(q\right)\stackrel{\·\·}{q}+C_1\left(q\right)\stackrel{\·}{q}+C_2(q,\stackrel{\·}{q})\stackrel{\·}{q}+F_v\stackrel{\·}{q}+F_{s}sign\left(\stackrel{\·}{q}\right)+G\left(q\right)=\mathrm{\τ}---\left(1\right)$

其中，B(q)为N×N惯量矩阵，C₁(q)为N×N离心力矩阵，为N×N科里奥利力矩阵，F_v为N×N粘摩擦力矩阵，F_s为N×N静摩擦力矩阵，G(q)为N×1重力向量，τ为N×1关节力矩向量，表示以为自变量的符号函数，为加速度信息，为速度信息，q为位置信息，N为机器人关节数。

3.根据权利要求2所述的机器人零力控制方法，其特征在于，所述步骤S2包括：

S21、将公式(1)中静摩擦力设为0，并将动力学模型中的设为期望加速度，设为速度反馈q设为位置反馈q_f，得到如公式(2)所示的关节作用力力矩τ_nsf：

$B\left(q_f\right){\stackrel{\·\·}{q}}_r+C_1\left(q_f\right){\stackrel{\·}{q}}_f+C_2(q_f,\stackrel{\·}{q}){\stackrel{\·}{q}}_f+F_v{\stackrel{\·}{q}}_f+G\left(q_f\right)={\mathrm{\τ}}_{nsf}---\left(1\right);$

S22、根据力矩常数以及反馈电流计算反馈力矩τ_f：K_tI_f＝τ_f；

S23、根据步骤S21中计算的关节作用力力矩τ_nsf以及步骤S22中计算的反馈力矩τ_f计算得到外力力矩τ_e：τ_f-τ_nsf＝τ_e。

4.根据权利要求3所述的机器人零力控制方法，其特征在于，所述步骤S3包括：

S31、将公式(1)中静摩擦力以外的作用力设为0，并将动力学模型中的设为速度反馈得到静摩擦力力矩τ_sf：

S32、根据所述步骤S41中的静摩擦力力矩τ_sf以及步骤S21中的关节作用力力矩τ_nsf计算外力力矩死区作用域[τ_dn,τ_dp]，其中，τ_dn,τ_dp基于以下公式确定：

$\left\{\begin{array}{c}{k}_1{\mathrm{\τ}}_{nsf}-\left|{\mathrm{\τ}}_{sf}\right|={\mathrm{\τ}}_{dn}\\ k_2{\mathrm{\τ}}_{nsf}+\left|{\mathrm{\τ}}_{sf}\right|={\mathrm{\τ}}_{dp}\end{array}\right.$

其中，k₁、k₂表示力矩误差校正系数；

S33、根据步骤S23中计算得到的外力力矩τ_e，基于以下公式计算得到速度指令

$\left\{\begin{array}{cc}{k}_{3}sign\left({\mathrm{\τ}}_e\right)|{\mathrm{\τ}}_e{|}^p={\stackrel{\·}{q}}_r,& {\mathrm{\τ}}_e\∉\[{\mathrm{\τ}}_{dn},{\mathrm{\τ}}_{dp}\]\\ 0={\stackrel{\·}{q}}_r,& {\mathrm{\τ}}_e\∈\[{\mathrm{\τ}}_{dn},{\mathrm{\τ}}_{dp}\]\end{array}\right.$

其中，k₃为力矩速度映射系数，p∈(0,1)。

5.根据权利要求1所述的机器人零力控制方法，其特征在于，所述步骤S4基于以下公式计算得到位置指令q_r；

$q_f+{\stackrel{\·}{\hat{q}}}_r\mathrm{\Δ}t=q_r$

其中，Δt为位置插补周期，q_f为位置反馈，为速度指令。

6.一种机器人零力控制系统，其特征在于，包括：

动力学模型模块，用于基于惯性力、离心力、科里奥利力、粘摩擦力、静摩擦力、重力，建立机器人动力学模型；

外力力矩计算模块，用于基于所述动力学模型模块中的机器人动力学模型以及反馈数据，计算机器人各个关节的外力力矩；

力矩速度映射模块，用于根据外力力矩计算模块中计算的外力力矩计算速度指令；

速度位置映射模块，用于根据速度指令以及位置反馈计算位置指令。

7.根据权利要求6所述的机器人零力控制系统，其特征在于，所述动力学模型为：

$B\left(q\right)\stackrel{\·\·}{q}+C_1\left(q\right)\stackrel{\·}{q}+C_2(q,\stackrel{\·}{q})\stackrel{\·}{q}+F_v\stackrel{\·}{q}+F_{s}sign\left(\stackrel{\·}{q}\right)+G\left(q\right)=\mathrm{\τ}---\left(1\right)$

其中，B(q)为N×N惯量矩阵，C₁(q)为N×N离心力矩阵，为N×N科里奥利力矩阵，F_v为N×N粘摩擦力矩阵，F_s为N×N静摩擦力矩阵，G(q)为N×1重力向量，τ为N×1关节力矩向量，表示以为自变量的符号函数，为加速度信息，为速度信息，q为位置信息，N为机器人关节数。

8.根据权利要求7所述的机器人零力控制系统，其特征在于，所述外力力矩计算模块包括：

关节作用力力矩计算单元，用于将公式(1)中静摩擦力设为0，并将动力学模型中的设为期望加速度，设为速度反馈q设为位置反馈q_f，得到如公式(2)所示的关节作用力力矩τ_nsf：

$B\left(q_f\right){\stackrel{\·\·}{q}}_r+C_1\left(q_f\right){\stackrel{\·}{q}}_f+C_2(q_f,{\stackrel{\·}{q}}_f){\stackrel{\·}{q}}_f+F_v{\stackrel{\·}{q}}_f+G\left(q_f\right)={\mathrm{\τ}}_{nsf}---\left(1\right);$

反馈力矩计算单元，用于根据力矩常数以及反馈电流计算反馈力矩τ_f：K_tI_f＝τ_f；

外力力矩计算单元，用于根据关节作用力力矩计算单元计算的关节作用力力矩τ_nsf以及反馈力矩计算单元计算的反馈力矩τ_f计算得到外力力矩τ_e：τ_f-τ_nsf＝τ_e。

9.根据权利要求8所述的机器人零力控制系统，其特征在于，所述力矩速度映射模块包括：

静摩擦力力矩计算单元，用于将公式(1)中静摩擦力以外的作用力设为0，并将动力学模型中的设为速度反馈得到静摩擦力力矩τ_sf：

外力力矩死区作用域计算单元，用于根据所述静摩擦力力矩计算单元计算的静摩擦力力矩τ_sf以及关节作用力力矩计算单元计算的关节作用力力矩τ_nsf计算外力力矩死区作用域[τ_dn,τ_dp]，其中，τ_dn,τ_dp基于以下公式确定：

$\left\{\begin{array}{c}{k}_1{\mathrm{\τ}}_{nsf}-\left|{\mathrm{\τ}}_{sf}\right|={\mathrm{\τ}}_{dn}\\ k_2{\mathrm{\τ}}_{nsf}+\left|{\mathrm{\τ}}_{sf}\right|={\mathrm{\τ}}_{dp}\end{array}\right.$

其中，k₁、k₂表示力矩误差校正系数；

速度指令计算单元，用于根据外力力矩计算单元中计算得到的外力力矩τ_e，基于以下公式计算得到速度指令

$\left\{\begin{array}{cc}{k}_{3}sign\left({\mathrm{\τ}}_e\right)|{\mathrm{\τ}}_e{|}^p={\stackrel{\·}{q}}_r,& {\mathrm{\τ}}_e\∉\[{\mathrm{\τ}}_{dn},{\mathrm{\τ}}_{dp}\]\\ 0={\stackrel{\·}{q}}_r,& {\mathrm{\τ}}_e\∈\[{\mathrm{\τ}}_{dn},{\mathrm{\τ}}_{dp}\]\end{array}\right.$

其中，k₃为力矩速度映射系数，p∈(0,1)。

10.根据权利要求6所述的机器人零力控制系统，其特征在于，所述的根据速度指令以及位置反馈计算位置指令包括：基于以下公式计算得到位置指令q_r；

$q_f+{\stackrel{\·}{\hat{q}}}_r\mathrm{\Δ}t=q_r$

其中，Δt为位置插补周期，q_f为位置反馈，为速度指令。