CN114355964A - 多自由度单腿运动学求解方法、装置和机器人 - Google Patents

多自由度单腿运动学求解方法、装置和机器人 Download PDF

Info

Publication number
CN114355964A
CN114355964A CN202111642269.6A CN202111642269A CN114355964A CN 114355964 A CN114355964 A CN 114355964A CN 202111642269 A CN202111642269 A CN 202111642269A CN 114355964 A CN114355964 A CN 114355964A
Authority
CN
China
Prior art keywords
joint
hip
ankle
anterior
angle
Prior art date
Legal status (The legal status 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 status listed.)
Granted
Application number
CN202111642269.6A
Other languages
English (en)
Other versions
CN114355964B (zh
Inventor
周江琛
陈春玉
刘益彰
葛利刚
闫萌
熊友军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Youbixuan Intelligent Robot Co ltd
Original Assignee
Ubtech Robotics Corp
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 Ubtech Robotics Corp filed Critical Ubtech Robotics Corp
Priority to CN202111642269.6A priority Critical patent/CN114355964B/zh
Publication of CN114355964A publication Critical patent/CN114355964A/zh
Application granted granted Critical
Publication of CN114355964B publication Critical patent/CN114355964B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Manipulator (AREA)

Abstract

本申请实施例提供一种多自由度单腿运动学求解方法、装置和机器人,所述单腿包括依次设置的三个旋转轴不相交的髋侧关节、髋转关节和髋前关节,及膝前关节,以及两个旋转轴相交的踝前关节和踝侧关节,该方法包括:从所有关节的旋转轴中选取若干节点,将所有节点依次连接,得到基于脚踝局部坐标系建立的单腿的简化运动学模型,根据获取的起始节点与末端节点在世界坐标系中的位姿,基于简化运动学模型中各个节点的位姿与关节角度之间的几何关系,计算各个关节的关节角度。该方法可以有效解决具有三旋转轴不相交的X‑Z‑Y构型的多连杆机构的逆运动学求解问题。

Description

多自由度单腿运动学求解方法、装置和机器人
技术领域
本申请涉及机器人控制技术领域,尤其涉及一种多自由度单腿运动学求解方法、装置和机器人。
背景技术
双足仿人机器人是一种能够模仿人类站立和行走跳跃等运动特点的机器人。在逆运动学求解过程中,通常采用给定机器人的位置和姿态,利用动力学模型来进行步态规划设计,得到各个关节的期望旋转角度(亦称关节角度),然后通过电机驱动相应关节的旋转轴进行旋转。
双足机器人的左右腿通常为设为各有六自由度(亦称关节)的多连杆结构,其中,髋关节具有3个自由度,膝关节具有1个自由度,踝关节有2个自由度,传统的逆运动学求解方法一般要求腿部髋关节的各旋转轴相交于一点。然而,对于腿部髋关节采用了先翻滚(roll)后偏摆(yaw)再俯仰(pitch)的自由度结构设计,导致髋关节的各旋转轴不相交一点,对于这类特殊构型的机器人,利用传统的DH参数方法无法进行有效求解。
发明内容
本申请实施例提供一种多自由度单腿运动学求解方法、装置和机器人,该方法可以有效解决具有三个旋转轴不相交的X-Z-Y构型的机器人的逆运动学求解问题。
第一方面,本申请实施例提供一种多自由度单腿运动学求解方法,所述单腿包括依次设置的三个旋转轴不相交的髋侧关节、髋转关节和髋前关节,及膝前关节,以及两个旋转轴相交的踝前关节和踝侧关节,所述方法包括:
从所有关节的旋转轴中选取若干节点,将所有节点依次连接,得到基于脚踝局部坐标系建立的所述单腿的简化运动学模型,其中,所述若干节点包括作为起始节点的所述踝前关节和所述踝侧关节的旋转轴交点、作为第一中间节点的所述膝前关节的旋转轴中点、作为第二中间节点的所述髋前关节与所述髋转关节的旋转轴交点、以及作为末端节点的所述髋侧关节和所述髋转关节的旋转轴交点;
根据获取的所述起始节点与所述末端节点在世界坐标系中的位姿,基于所述简化运动学模型中各个节点的位姿与关节角度之间的几何关系,计算所述髋侧关节、所述髋转关节、所述髋前关节、所述膝前关节、所述踝前关节和所述踝侧关节各自的关节角度。
第二方面,本申请实施例还提供一种机器人运动控制方法,所述机器人包括至少一条多自由度单腿结构,所述方法包括:
获取机器人的躯干与至少一足底末端在世界坐标系下的位姿;
根据所述躯干的所述位姿确定机器人各腿的髋侧关节和髋转关节的旋转轴交点在世界坐标系下的位姿,以及根据对应足底末端的所述位姿确定对应腿的踝前关节和踝侧关节的旋转轴交点在世界坐标系下的位姿;
通过上述的多自由度单腿运动学求解方法计算所述机器人每条腿的各个关节角度;
根据计算的所有关节角度对机器人进行相应控制。
第三方面,本申请实施例还提供一种多自由度单腿运动学求解装置,所述单腿包括依次设置的三个旋转轴不相交的髋侧关节、髋转关节和髋前关节,及膝前关节,以及两个旋转轴相交的踝前关节和踝侧关节,所述装置包括:
简化模型构建模块,用于从所有关节的旋转轴中选取若干节点,将所有节点依次连接,得到基于脚踝局部坐标系建立的所述单腿的简化运动学模型,其中,所述若干节点包括作为起始节点的所述踝前关节和所述踝侧关节的旋转轴交点、作为第一中间节点的所述膝前关节的旋转轴中点、作为第二中间节点的所述髋前关节与所述髋转关节的旋转轴交点、以及作为末端节点的所述髋侧关节和所述髋转关节的旋转轴交点;
关节角度求解模块,用于根据获取的所述起始节点与所述末端节点在世界坐标系中的位姿,基于所述简化运动学模型中各个节点的位姿与关节角度之间的几何关系,计算所述髋侧关节、所述髋转关节、所述髋前关节、所述膝前关节、所述踝前关节和所述踝侧关节各自的关节角度。
第四方面,本申请实施例还提供一种机器人,所述机器人包括处理器和存储器,所述存储器存储有计算机程序,所述处理器用于执行所述计算机程序以实施上述的多自由度单腿运动学求解方法或机器人运动控制方法。
第五方面,本申请实施例还提供一种可读存储介质,其存储有计算机程序,所述计算机程序在处理器上执行时,实施上述的多自由度单腿运动学求解方法或机器人运动控制方法。
本申请的实施例具有如下有益效果:
本申请实施例的多自由度单腿运动学求解方法针对多自由度的单腿结构,该腿的髋侧关节、髋转关节和髋前关节三个旋转轴不相交,而踝前关节和踝侧关节的两个旋转轴相交,该运动学求解方法通过从所有关节的旋转轴中选取若干节点,以构建单腿的简化运动学模型,进而根据获取的起始节点与末端节点在世界坐标系中的位姿,基于简化运动学模型中各个节点的位姿与关节角度之间的几何关系,计算各个关节的关节角度。该方法可以有效解决具有三旋转轴不相交的X-Z-Y构型的多连杆机构的逆运动学求解问题,直观简洁,且计算效率更高。
附图说明
为了更清楚地说明本申请实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本申请的某些实施例,因此不应被看作是对范围的限定,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他相关的附图。
图1示出了本申请实施例的机器人的结构示意图;
图2示出了本申请实施例的髋关节的三个自由度不交于一点的双足机器人模型;
图3示出了本申请实施例的多自由度单腿运动学求解方法的流程图;
图4示出了本申请实施例的六自由度单腿的简化运动学模型;
图5示出了本申请实施例的计算所有关节角度的流程图;
图6示出了本申请实施例的包含等效角度在内的三个关节角度求解的流程图;
图7示出了本申请实施例的求解三个俯仰角度的流程图;
图8示出了本申请实施例的机器人运动控制方法的流程图;
图9示出了本申请实施例多自由度单腿运动学求解装置的结构示意图;
图10示出了本申请实施例的机器人运动学求解装置的结构示意图。
具体实施方式
下面将结合本申请实施例中附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。
通常在此处附图中描述和示出的本申请实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本申请的实施例的详细描述并非旨在限制要求保护的本申请的范围,而是仅仅表示本申请的选定实施例。基于本申请的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本申请保护的范围。
在下文中,可在本申请的各种实施例中使用的术语“包括”、“具有”及其同源词仅意在表示特定特征、数字、步骤、操作、元件、组件或前述项的组合,并且不应被理解为首先排除一个或更多个其它特征、数字、步骤、操作、元件、组件或前述项的组合的存在或增加一个或更多个特征、数字、步骤、操作、元件、组件或前述项的组合的可能性。
此外,术语“第一”、“第二”、“第三”等仅用于区分描述,而不能理解为指示或暗示相对重要性。
除非另有限定,否则在这里使用的所有术语(包括技术术语和科学术语)具有与本申请的各种实施例所属领域普通技术人员通常理解的含义相同的含义。所述术语(诸如在一般使用的词典中限定的术语)将被解释为具有与在相关技术领域中的语境含义相同的含义并且将不被解释为具有理想化的含义或过于正式的含义,除非在本申请的各种实施例中被清楚地限定。
下面结合附图,对本申请的一些实施方式作详细说明。在不冲突的情况下,下述的实施例及实施例中的特征可以相互结合。
请参照图1,为本申请实施例提供的机器人的一种结构示意图。在本实施例中,该机器人10包括处理器11、存储器12和感知单元13等,其中,存储器12和感知单元13均与处理器11连接,存储器12存储有相应的计算机程序,处理器11用于执行该计算机程序以实施本申请实施例中的多自由度单腿运动学求解方法或机器人运动控制方法,从而有效解决具有三个旋转轴不相交的X-Z-Y构型的机器人的逆运动学求解问题。
其中,处理器11可以是一种具有信号的处理能力的集成电路芯片。处理器可以是通用处理器,包括中央处理器(Central Processing Unit,CPU)、图形处理器(GraphicsProcessing Unit,GPU)及网络处理器(Network Processor,NP)、数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件中的至少一种。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等,可以实现或者执行本申请实施例中的公开的各方法、步骤及逻辑框图。
其中,存储器12可以但不限于是,随机存取存储器(Random Access Memory,RAM),只读存储器(Read Only Memory,ROM),可编程只读存储器(Programmable Read-OnlyMemory,PROM),可擦除只读存储器(Erasable Programmable Read-Only Memory,EPROM),电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory,EEPROM)等。其中,存储器用于存储计算机程序,处理器在接收到执行指令后,可相应地执行计算机程序。
其中,感知单元13主要指在机器人身上设置的传感器,这些传感器用于将必要的外部信息以及自身状态信息传递给机器人的控制系统,从而为机器人的决策提供必要的条件。在一些情况下,这些传感器可分为内部传感器和外部传感器,其中,内部传感器用于测量机器人自身状态的功能元件,其功能是测量运动学量和力学量,用于机器人感知自身的运动状态,使得机器人可以按照规定的位置、轨迹和速度等参数运动。例如,内部传感器包括但不限于,设于相应关节或躯干处的位置或角度传感器、六维力传感器、压力传感器等。而外部传感器主要是感知机器人自身和环境之间的相互信息,例如包括视觉等信息,具体可包括但不限于接近传感器、视觉传感器等。应当明白的是,上述列举的这些传感器可根据机器人的实际应用场景而选择性设置。
在本申请实施例中,该机器人10可以指具有一个多自由度连杆结构的机械臂结构,如工业用机械臂、医学用机械臂等,也可以指具有两个及两个以上的多自由度连杆结构的机器人,如双足机器人(亦称仿人机器人)、三足机器人、四足机器人、六足机器人等,这里不作限定。值得注意的是,本申请实施例中的机器人应当包含存在连续的三个自由度的旋转轴不相交于一点的X-Z-Y构型,以便能采用本申请实施例的多自由度单腿运动学求解方法来有效求解出各个关节角度。
为了更直观地理解上述的三旋转轴不相交于一点的X-Z-Y构型的结构,这里以一双足机器人为例,该双足机器人包括左腿和右腿,每条腿均具有6个自由度,其中,髋关节具有3个自由度,膝关节具有1个自由度,踝关节有2个自由度。为了让机器人的整体结构紧凑,如图2所示,该双足机器人每条腿的髋关节的三个自由度采用了先翻滚(roll)后偏摆(yaw)再俯仰(pitch)的结构设计,对应为三个旋转轴不相交的髋侧关节、髋转关节和髋前关节;膝关节具有俯仰自由度,又称膝前关节;而踝关节具有俯仰和翻滚两个自由度,分别为踝前关节和踝侧关节,其中,踝前关节和踝侧关节的旋转轴相交。
为了方便理解,对于机器人的运动,笛卡尔坐标系建立如下:X轴方向指向机器人的正前方,Y轴方向指向机器人的左侧,Z轴方向指向机器人的躯干上方。对于上述的髋关节,其三个自由度的旋转角度依次可称为翻滚角、偏摆角和俯仰角。可以理解,翻滚角是指在笛卡尔空间坐标系中绕X轴旋转的角度;偏摆角是指在笛卡尔空间坐标系中绕Z轴旋转的角度;俯仰角是指在笛卡尔空间坐标系中绕Y轴旋转的角度。
应当明白的是,本申请实施例的多自由度单腿运动学求解方法除了可以应用于上述具有三旋转轴不相交的X-Z-Y构型的双足机器人腿部结构外,该方法也可以扩展到具有相似构型的其他机械臂,这里不作限定。
基于上述的机器人结构,为了方便理解,这里以具有多自由度单腿结构为例进行运动学求解过程进行说明。下面对本申请实施例的多自由度单腿运动学求解方法进行具体阐述。该方法可以有效完成此类特殊构型的多自由度连杆机构的运动学逆解,相比于传统的DH方法更为直观简洁,相比雅可比迭代的方法计算效率更高,具备良好的普适性和实用性。
请参照图3,本实施例提出一种多自由度单腿运动学求解方法的流程图。示范性地,该多自由度单腿运动学求解方法包括步骤S110~S120:
S110,从所有关节的旋转轴中选取若干节点,将所有节点依次连接,得到基于脚踝局部坐标系建立的该单腿的简化运动学模型。
其中,节点是指在连杆结构中选取的一些关键点,通过这些关键点的依次连接,可以将该多自由度的腿部进行结构简化,以便于进行运动学分析。
在本实施例中,该节点的数量不作限定,可以根据实际需求来选取,通常地,选取的节点越多,描述的结构信息也越多,计算量也往往越大,因此可以根据实际需求来选取合适数量的节点。例如,对于具有六自由度的单腿结构,可选取3~5个节点不等,这里不作具体限定。
其中,脚踝局部坐标系是指对腿部的脚踝这一位置建立的局部坐标系。可以理解,这里由脚踝局部坐标系来建立简化模型,当机器人在运动时,局部坐标系随其依附的机器人连杆的位置而变化,以便对每个时刻的运动状态进行相应分析。
当机器人开始运动时,以具有六自由度的单腿结构为例,腿部的各个关节角度将变化,对于任意时刻的腿部姿态,均可通过本实施例提出的简化运动学模型来描述。在一种实施例中,示范性地,可从单腿结构中选取四个节点并将这些节点依次连接,得到该单腿的简化运动学模型。具体地,包括踝前关节和踝侧关节的旋转轴交点、膝前关节的旋转轴中点、髋前关节与髋转关节的旋转轴交点、髋侧关节和髋转关节的旋转轴交点,如图4所示,从脚踝关节到髋关节,这四个节点依次定义为起始节点O、第一中间节点A、第二中间节点B和末端节点C。进而,由线段O、A、B和C依次连接所呈现的构型即为六自由度单腿的简化运动学模型。
应当明白的是,将脚踝位置处选取的节点定义为起始节点,这样可以与脚踝局部坐标系的原点选取为同一点,以方便后续分析。可以理解,这些节点的顺序并不仅限于上述示例,也可以从髋关节到脚踝关节来依次定义,例如,起始节点也可以是髋侧关节和髋转关节的旋转轴交点,而末端节点为踝前关节和踝侧关节的旋转轴交点等。
可以理解,该简化运动学模型中各个节点的位姿与关节角度之间存在相应映射关系。如图4所示,在构建的简化运动学模型中,可以描述出六个关节的旋转角度(即关节角度)q1~q6。如图4所示,其中,q1、q2和q3分别表示髋侧关节、髋转关节和髋前关节的关节角度,q4表示膝前关节的关节角度,q5和q6分别表示踝前关节和踝侧关节的关节角度。可以明白的是,由于髋前关节、膝前关节和踝前关节均为俯仰角,均绕Y轴旋转而不会超过XOZ平面,因此,这三个关节的旋转角度可以投影到同一平面以进行合并处理。在简化运动学模型中,如图4所示,表现为节点A、B和C投影到同一平面XOZ’上,合并后的旋转角度记为等效角度,该等效角度即为由节点O指向节点B的连线与Z’轴的夹角α。
于是,根据齐次变换原理,对于合并得到的等效角度,与髋前关节、膝前关节和踝前关节的关节角度,将满足如下等式关系:
Rα=R5*R4*R3。
其中,Rα为等效角度α对应的旋转矩阵,R3为髋前关节的关节角度q3对应的旋转矩阵,R4为膝前关节的关节角度q4对应的旋转矩阵,R5为踝前关节的关节角度q5对应的旋转矩阵。可以理解,一个向量绕旋转轴旋转给定的角度后,该旋转角度可通过罗德里格斯(Rodrigues)转换公式转换为对应的旋转矩阵。
此外,由于踝侧关节与髋侧关节均为绕X轴旋转的翻滚角,根据上述模型中末端节点和起始节点的几何位置,此时踝侧关节的关节角度与节点A相对于节点O的位置存在映射关系,具体地,该踝侧关节的关节角度q6的正切值等于末端节点与起始节点之间的相对位置中的Z方向分量与Y方向分量的比值。若末端节点A与起始节点O之间的相对位置为(xa,ya,za),其中,xa,ya,za分别为在X、Y和Z方向上的分量,通过表达式描述上述映射关系,即有:
q6=arctan(ya/za)。
于是,基于上述相应的几何等式关系来计算各个关节角度。
S120,根据获取的起始节点与末端节点在世界坐标系中的位姿,基于该简化运动学模型中各个节点的位姿与关节角度之间的几何关系,计算各个关节角度。
其中,起始节点与末端节点在世界坐标系中的位姿可通过预先给定或规划的机器人相应部位的位姿直接得到或间接计算得到,具体地,节点O在世界坐标系中的位姿与机器人脚底的位姿有关,例如,在某些情况下,可将脚底的位姿直接选取为踝前关节和踝侧关节的旋转轴交点(即起始节点O)的位姿;而节点A在世界坐标系中的位姿与机器人躯干的位姿有关,躯干通常与髋侧关节直接相连,故可利用躯干的位姿来得到髋侧关节的位姿。
在一种实施方式中,对于上述步骤S120,如图5所示,可包括以下子步骤S210~250:
S210,根据起始节点与末端节点在世界坐标系中的位姿,计算末端节点与起始节点之间的相对位置和相对姿态、以及末端节点相对于起始节点对应的第一旋转矩阵。
其中,第一旋转矩阵也即起始节点旋转到末端节点时对应的旋转矩阵。示范性地,根据给定的节点A与节点O各自的位姿,可计算得到节点A相对于节点O的位置(xa,ya,za)与姿态
Figure BDA0003444105130000121
其中,
Figure BDA0003444105130000122
分别为绕X轴、Z轴和Y轴旋转的角度。进而,根据相对姿态通过Rodrigues公式可计算出末端节点相对于起始节点对应的第一旋转矩阵,这里记为ROA
S220,根据起始节点与末端节点之间的相对位置基于第一等式方程,计算踝侧关节的关节角度,进而计算踝侧关节的关节角度对应的第二旋转矩阵。
其中,第一等式方程是指上述的踝侧关节的关节角度与起始节点与末端节点之间的相对位置之间的映射关系。示范性地,在计算出节点A相对于节点O的位置后,可通过上述映射关系计算出关节角度q6。进而,在得到关节角度q6后,同样可利用Rodrigues公式可计算出对应的旋转矩阵R6,即上述的第二旋转矩阵。
S230,根据第一旋转矩阵和第二旋转矩阵,基于髋侧关节到踝侧关节的姿态转换对应的旋转矩阵关系,计算第三旋转矩阵。
考虑到髋侧关节到踝侧关节的姿态转换,根据机器人运动学原理,相应地,各个关节角度的旋转矩阵将满足如下关系:
ROA=R6R5R4R3R2R1
其中,R1~R6依次为关节角度q1~q6对应的旋转矩阵。
由于关节角度q3~q5可以在同一平面合并得到一等效角度α,故上述表达式也可描述为:
ROA=R6RαR2R1
进一步地,设中间变量Rtemp=RαR2R1,即有:Rtemp=R6 -1ROA
示范性地,在得到旋转矩阵ROA和R6后,可计算得到中间变量Rtemp,即上述的第三旋转矩阵。
S240,根据第三旋转矩阵,计算髋侧关节的关节角度、髋转关节的关节角度、以及由髋前关节、膝前关节和踝前关节的关节角度在同一平面合并后的等效角度。
在一种实施方式中,对于步骤S240,如图6所示,可包括如下子步骤S241~S242:
S241,基于第三旋转矩阵等于髋侧关节的关节角度、膝前关节的关节角度和等效角度各自对应的旋转矩阵的乘积关系,从第三旋转矩阵中分别确定与髋侧关节的关节角度、髋转关节的关节角度和等效角度相关的向量。
由上述等式Rtemp=RαR2R1可知,故这三个角度对应的旋转矩阵经矩阵运算后相关参数将体现在Rtemp中相应向量中,由于关节角度q1、q2和等效角度α的旋转轴不相交,故可利用相应向量中的相应元素来分别求解各个关节角度。示范性地,可分别利用第三旋转矩阵的第二列向量中的相应元素来计算髋侧关节和髋转关节的关节角度,以及利用第三旋转矩阵中的第一行向量的相应元素来计算等效角度。
可以理解,本实施例中的旋转矩阵中的元素位置与从机器人上选取的起始节点和末端终点的位置有关,当起始节点和末端终点的位置对换时,旋转矩阵中的与对应关节角度或等效角度相关的向量位置也将相应变化。对应地,用于计算相应角度的元素位置也将适应性调整。
S242,根据髋侧关节、髋转关节和等效角度各自对应的向量中的相应元素,分别计算髋侧关节的关节角度、髋转关节的关节角度和该等效角度。
以得到的第三旋转矩阵Rtemp为例,在一种实施方式中,髋侧关节的关节角度q1可通过如下公式计算:
q1=arctan(Rtemp(3,2)/Rtemp(2,2));
其中,Rtemp(i,j)表示取矩阵Rtemp中的第i行第j列的元素。
在一种实施方式中,髋转关节的关节角度q2可通过如下公式计算:
Figure BDA0003444105130000141
可以理解,在计算关节角度q2时,增加的负号是根据关节角度的方向来确定的。
在一种实施方式中,等效角度α可通过如下公式计算:
α=arctan(Rtemp(1,3)/Rtemp(1,1))。
S250,计算等效角度对应的等效旋转矩阵,并根据等效旋转矩阵基于第二等式方程,计算髋前关节、膝前关节、踝前关节各自的关节角度。
同理,利用Rodrigues公式可计算出等效角度α对应的等效旋转矩阵Rα,进而,利用等效旋转矩阵Rα中的相应元素来计算得到上述三个俯仰角。
在一种实施方式中,对于步骤S250,如图7所示,可包括如下子步骤S251~S252:
S251,根据上述第二等式方程,从等效旋转矩阵中确定与髋前关节、膝前关节、踝前关节的关节角度相关的向量。
由上述等式Rα=R5*R4*R3可知,故这三个角度对应的旋转矩阵经矩阵运算后相关参数将体现在Rα的相应向量中,示范性地,可分别利用等效旋转矩阵的第一行向量中的相应元素来计算髋前关节的关节角度,利用等效旋转矩阵的第二列向量的相应元素来计算膝前关节的关节角度,以及等效旋转矩阵的第二列向量的相应元素来计算踝前关节的关节角度。
S252,根据髋前关节、膝前关节和踝前关节各自对应的向量中的相应元素,分别计算髋前关节、膝前关节、踝前关节各自的关节角度。
以上述的等效旋转矩阵Rα为例,在一种实施方式中,髋前关节的关节角度q3可通过如下公式计算:
q3=arctan(Rα(1,3)/Rα(1,1));
其中,Rα(i,j)表示取矩阵Rα中的第i行第j列的元素。
在一种实施方式中,膝前关节的关节角度q4可通过如下公式计算:
Figure BDA0003444105130000151
在一种实施方式中,踝前关节的关节角度q5可通过如下公式计算:
q5=arctan(Rα(3,2)/Rα(2,2))。
至此,计算得到六自由度单腿结构的6个关节角度。可以理解,上述的各个关节角度的计算方式是通过相应旋转矩阵中的相应元素来直接计算得到的,而旋转矩阵的获取与相应的关节角度或等效角度的设置顺序有关。
本申请实施例的多自由度单腿运动学求解方法对于具有三个旋转轴不相交于一点的X-Z-Y构型的多自由度连杆结构,通过根据当前的腿部姿态来选取一些关键节点以构建相应的简化运动学模型,进而利用模型满足的相应几何关系来计算各个关节角度,相比于传统的DH方法更为直观简洁,相比利用雅可比矩阵进行迭代求解的方法计算效率更高等。
请参照图8,基于上述实施例的多自由度单腿运动学求解方法,本申请实施例还提出一种机器人运动控制方法,该机器人包括至少一条多自由度单腿结构,例如,可以是如图2所示的双足机器人等。在该机器人运动控制方法中,这里对于机器人的每条腿均采用上述的运动学求解方法进行关节角度的求解。
示范性地,该机器人运动控制方法包括步骤S310~330:
S310,获取机器人的躯干与至少一足底末端在世界坐标系下的位姿。
在机器人运动控制中,通常会预先规划机器人的躯干和/或足底末端在世界坐标系下的运动轨迹。示范性地,从规划的运动轨迹中可以获得机器人的躯干以及足底末端在世界坐标系下的位姿。可以理解,该机器人可以是仅包含一条具有上述实施例中的三旋转轴不相交的X-Z-Y构型的多自由度腿,也可以包含多条具有上述特殊构型的多自由度腿,这里不作限定。
以双足机器人为例,双足机器人的双腿均具有上述特殊构型,由于双腿为对称结构,对此,本实施例可以针对每条腿采用上述实施例中的多自由度单腿运动学求解方法来计算出当前腿中的各个关节角度,从而得到机器人的所有关节角度。
S320,根据躯干的位姿确定机器人各腿的髋侧关节和髋转关节的旋转轴交点在世界坐标系下的位姿,以及根据足底末端的位姿确定对应腿的踝前关节和踝侧关节的旋转轴交点在世界坐标系下的位姿。
由于机器人的躯干与每条腿的髋侧关节是直接相连的,在一种实施方式中,可直接将躯干的位置以及躯干质心到髋侧关节的旋转轴中点的距离,计算出髋侧关节的旋转轴中心的位置,同时将躯干的姿态作为髋侧关节的旋转轴中心的姿态,由此得到髋侧关节的旋转轴中心的位置与姿态,也即髋侧关节和髋转关节的旋转轴交点在世界坐标系下的位姿。
同理,由于机器人的足底末端与踝侧关节是直接相连的,在一种实施方式中,若将踝关节的局部坐标系作为上述的足底末端,而踝关节的局部坐标系又选取为踝前关节和踝侧关节的旋转轴交点,则可直接将足底末端的位姿作为该踝侧关节的旋转轴中点的位姿,也为踝前关节和踝侧关节的旋转轴交点的位姿。
S330,通过多自由度单腿运动学求解方法计算每条腿的各个关节角度。
S340,根据计算的所有关节角度对机器人进行相应控制。
于是,根据上述步骤确定的简化运动学模型中的起始节点和末端节点在世界坐标系下的位姿,可采用上述实施例的多自由度单腿运动学求解方法进行逆运动学求解,得到当前腿中各个关节的关节角度。进而,根据该关节角度来控制对应关节的旋转轴旋转,以使机器人运动到期望位置,从而实现相应的运动控制。
请参照图9,基于上述实施例的多自由度单腿运动学求解方法,本申请实施例还提出一种多自由度单腿运动学求解装置100,其中,该单腿包括依次设置的三个旋转轴不相交的髋侧关节、髋转关节和髋前关节,及膝前关节,以及两个旋转轴相交的踝前关节和踝侧关节。下面对本申请实施例提供的多自由度单腿运动学求解装置100的具体组成进行相应描述。
简化模型构建模块110,用于从所有关节的旋转轴中选取若干节点,将所有节点依次连接,得到基于脚踝局部坐标系建立的所述单腿的简化运动学模型,其中,所述若干节点包括作为起始节点的所述踝前关节和所述踝侧关节的旋转轴交点、作为第一中间节点的所述膝前关节的旋转轴中点、作为第二中间节点的所述髋前关节与所述髋转关节的旋转轴交点、以及作为末端节点的所述髋侧关节和所述髋转关节的旋转轴交点;
关节角度求解模块120,用于根据获取的所述起始节点与所述末端节点在世界坐标系中的位姿,基于所述简化运动学模型中各个节点的位姿与关节角度之间的几何关系,计算所述髋侧关节、所述髋转关节、所述髋前关节、所述膝前关节、所述踝前关节和所述踝侧关节各自的关节角度。
可以理解,本实施例的装置对应于上述实施例的多自由度单腿运动学求解方法,上述实施例中的可选项同样适用于本实施例,在此不再重复描述。
请参照图10,基于上述实施例的机器人运动控制方法,本申请实施例还提出一种机器人运动控制装置200,下面对本申请实施例提供的机器人运动控制装置200的具体组成进行相应描述。
获取模块210,用于获取机器人的躯干与至少一足底末端在世界坐标系下的位姿。
确定模块220,用于根据所述躯干的所述位姿确定机器人各腿的髋侧关节和髋转关节的旋转轴交点在世界坐标系下的位姿,以及根据对应足底末端的所述位姿确定对应腿的踝前关节和踝侧关节的旋转轴交点在世界坐标系下的位姿。
求解模块230,用于通过上述的多自由度单腿运动学求解方法计算所述机器人每条腿的各个关节角度。
控制模块240,用于根据计算的所有关节角度对机器人进行相应控制。
可以理解,本实施例的装置对应于上述实施例的机器人运动控制方法,上述实施例中的可选项同样适用于本实施例,在此不再重复描述。
本申请还提供了一种可读存储介质,用于储存上述机器人中使用的所述计算机程序。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的,例如,附图中的流程图和结构图显示了根据本申请的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分,所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在作为替换的实现方式中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个连续的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,结构图和/或流程图中的每个方框、以及结构图和/或流程图中的方框的组合,可以用执行规定的功能或动作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
另外,在本申请各个实施例中的各功能模块或单元可以集成在一起形成一个独立的部分,也可以是各个模块单独存在,也可以两个或更多个模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是智能手机、个人计算机、服务器、或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本申请的保护范围之内。

Claims (11)

1.一种多自由度单腿运动学求解方法,其特征在于,所述单腿包括依次设置的三个旋转轴不相交的髋侧关节、髋转关节和髋前关节,及膝前关节,以及两个旋转轴相交的踝前关节和踝侧关节,所述方法包括:
从所有关节的旋转轴中选取若干节点,将所有节点依次连接,得到基于脚踝局部坐标系建立的所述单腿的简化运动学模型,其中,所述若干节点包括作为起始节点的所述踝前关节和所述踝侧关节的旋转轴交点、作为第一中间节点的所述膝前关节的旋转轴中点、作为第二中间节点的所述髋前关节与所述髋转关节的旋转轴交点、以及作为末端节点的所述髋侧关节和所述髋转关节的旋转轴交点;
根据获取的所述起始节点与所述末端节点在世界坐标系中的位姿,基于所述简化运动学模型中各个节点的位姿与关节角度之间的几何关系,计算所述髋侧关节、所述髋转关节、所述髋前关节、所述膝前关节、所述踝前关节和所述踝侧关节各自的关节角度。
2.根据权利要求1所述的多自由度单腿运动学求解方法,其特征在于,所述髋侧关节、所述髋转关节、所述髋前关节、所述膝前关节、所述踝前关节和所述踝侧关节各自的关节角度的计算,包括:
根据所述起始节点与所述末端节点在世界坐标系中的位姿,计算所述末端节点与所述起始节点之间的相对位置和相对姿态、以及所述末端节点相对于所述起始节点对应的第一旋转矩阵;
根据所述起始节点与所述末端节点之间的相对位置基于第一等式方程,计算所述踝侧关节的关节角度,进而计算所述踝侧关节的关节角度对应的第二旋转矩阵;
根据所述第一旋转矩阵和所述第二旋转矩阵,基于所述髋侧关节到所述踝侧关节的姿态转换对应的旋转矩阵关系,计算第三旋转矩阵;
根据所述第三旋转矩阵,计算所述髋侧关节的关节角度、所述髋转关节的关节角度、以及由所述髋前关节、所述膝前关节和所述踝前关节的关节角度在同一平面合并后的等效角度;
计算所述等效角度对应的等效旋转矩阵,并根据所述等效旋转矩阵基于第二等式方程,计算所述髋前关节、所述膝前关节、所述踝前关节各自的关节角度。
3.根据权利要求2所述的多自由度单腿运动学求解方法,其特征在于,所述第一等式方程通过预先构建得到,包括:
根据所述末端节点和所述起始节点的几何位置,构建所述踝侧关节的关节角度的正切值等于所述末端节点与所述起始节点之间的相对位置中的Z方向分量与Y方向分量的比值的方程以作为第一等式方程。
4.根据权利要求2所述的多自由度单腿运动学求解方法,其特征在于,所述第二等式方程通过预先构建得到,包括:
将所述末端节点、所述第一中间节点和所述第二中间节点投影到同一平面,构建所述髋前关节、所述膝前关节和所述踝前关节各自的关节角度在同一平面合并后的等效角度对应的等效旋转矩阵等于所述髋前关节、所述膝前关节、所述踝前关节各自的关节角度对应的旋转矩阵的乘积的方程以作为第二等式方程。
5.根据权利要求2所述的多自由度单腿运动学求解方法,其特征在于,所述根据所述第三旋转矩阵,计算所述髋侧关节的关节角度、所述髋转关节的关节角度、以及所述等效角度,包括:
基于所述第三旋转矩阵等于所述髋侧关节的关节角度、所述髋转关节的关节角度和所述等效角度各自对应的旋转矩阵的乘积关系,从所述第三旋转矩阵中分别确定与所述髋侧关节的关节角度、所述髋转关节的关节角度和所述等效角度相关的向量;
根据所述髋侧关节、所述髋转关节和所述等效角度各自对应的所述向量中的相应元素,分别计算所述髋侧关节的关节角度、所述髋转关节的关节角度和所述等效角度。
6.根据权利要求2或4所述的多自由度单腿运动学求解方法,其特征在于,所述根据所述等效旋转矩阵基于第二等式方程,计算所述髋前关节、所述膝前关节、所述踝前关节各自的关节角度,包括:
根据所述第二等式方程,从所述等效旋转矩阵中确定与所述髋前关节、所述膝前关节、所述踝前关节的关节角度相关的向量;
根据所述髋前关节、所述膝前关节和所述踝前关节各自对应的所述向量中的相应元素,分别计算所述髋前关节、所述膝前关节、所述踝前关节的关节角度。
7.根据权利要求2所述的多自由度单腿运动学求解方法,其特征在于,所述踝侧关节的关节角度对应的第二旋转矩阵、以及所述等效角度对应的等效旋转矩阵均通过罗德里格斯公式计算得到。
8.一种机器人运动控制方法,其特征在于,所述机器人包括至少一条多自由度单腿结构,所述方法包括:
获取机器人的躯干与至少一足底末端在世界坐标系下的位姿;
根据所述躯干的所述位姿确定机器人各腿的髋侧关节和髋转关节的旋转轴交点在世界坐标系下的位姿,以及根据对应足底末端的所述位姿确定对应腿的踝前关节和踝侧关节的旋转轴交点在世界坐标系下的位姿;
通过如权利要求1至7中任一项所述的多自由度单腿运动学求解方法计算所述机器人每条腿的各个关节角度;
根据计算的所有关节角度对机器人进行相应控制。
9.一种多自由度单腿运动学求解装置,其特征在于,所述单腿包括依次设置的三个旋转轴不相交的髋侧关节、髋转关节和髋前关节,及膝前关节,以及两个旋转轴相交的踝前关节和踝侧关节,所述装置包括:
简化模型构建模块,用于从所有关节的旋转轴中选取若干节点,将所有节点依次连接,得到基于脚踝局部坐标系建立的所述单腿的简化运动学模型,其中,所述若干节点包括作为起始节点的所述踝前关节和所述踝侧关节的旋转轴交点、作为第一中间节点的所述膝前关节的旋转轴中点、作为第二中间节点的所述髋前关节与所述髋转关节的旋转轴交点、以及作为末端节点的所述髋侧关节和所述髋转关节的旋转轴交点;
关节角度求解模块,用于根据获取的所述起始节点与所述末端节点在世界坐标系中的位姿,基于所述简化运动学模型中各个节点的位姿与关节角度之间的几何关系,计算所述髋侧关节、所述髋转关节、所述髋前关节、所述膝前关节、所述踝前关节和所述踝侧关节各自的关节角度。
10.一种机器人,其特征在于,所述机器人包括处理器和存储器,所述存储器存储有计算机程序,所述处理器用于执行所述计算机程序以实施权利要求1-7中任一项所述的多自由度单腿运动学求解方法或权利要求8中所述的机器人运动控制方法。
11.一种可读存储介质,其特征在于,其存储有计算机程序,所述计算机程序在处理器上执行时,实施根据权利要求1-7中任一项所述的多自由度单腿运动学求解方法或权利要求8中所述的机器人运动控制方法。
CN202111642269.6A 2021-12-29 2021-12-29 多自由度单腿运动学求解方法、装置和机器人 Active CN114355964B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111642269.6A CN114355964B (zh) 2021-12-29 2021-12-29 多自由度单腿运动学求解方法、装置和机器人

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111642269.6A CN114355964B (zh) 2021-12-29 2021-12-29 多自由度单腿运动学求解方法、装置和机器人

Publications (2)

Publication Number Publication Date
CN114355964A true CN114355964A (zh) 2022-04-15
CN114355964B CN114355964B (zh) 2023-08-18

Family

ID=81103039

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111642269.6A Active CN114355964B (zh) 2021-12-29 2021-12-29 多自由度单腿运动学求解方法、装置和机器人

Country Status (1)

Country Link
CN (1) CN114355964B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110039544A (zh) * 2019-04-28 2019-07-23 南京邮电大学 基于三次样条插值的仿人足球机器人步态规划
CN110989585A (zh) * 2019-11-28 2020-04-10 深圳市优必选科技股份有限公司 避奇异的步态规划方法、装置、可读存储介质及机器人
CN111538949A (zh) * 2020-07-10 2020-08-14 深圳市优必选科技股份有限公司 冗余机器人逆运动学求解方法、装置和冗余机器人
CN111914416A (zh) * 2020-07-24 2020-11-10 之江实验室 一种高能效轻量化结构双足机器人的逆运动学求解方法
CN112051741A (zh) * 2020-09-16 2020-12-08 北京理工大学 一种双足机器人动态运动生成与控制方法
CN112882485A (zh) * 2021-01-25 2021-06-01 河南大学 一种双足机器人逆运动学的几何计算方法
CN113835429A (zh) * 2021-09-08 2021-12-24 南方科技大学 一种仿生双足机器人的控制方法及仿生双足机器人
US20220203526A1 (en) * 2020-12-30 2022-06-30 Ubtech Robotics Corp Ltd Humanoid robot control method, humanoid robot using the same, and computer readable storage medium
WO2022227457A1 (zh) * 2021-04-30 2022-11-03 深圳市优必选科技股份有限公司 质心位姿估计方法, 装置, 计算机可读存储介质及机器人

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110039544A (zh) * 2019-04-28 2019-07-23 南京邮电大学 基于三次样条插值的仿人足球机器人步态规划
CN110989585A (zh) * 2019-11-28 2020-04-10 深圳市优必选科技股份有限公司 避奇异的步态规划方法、装置、可读存储介质及机器人
CN111538949A (zh) * 2020-07-10 2020-08-14 深圳市优必选科技股份有限公司 冗余机器人逆运动学求解方法、装置和冗余机器人
CN111914416A (zh) * 2020-07-24 2020-11-10 之江实验室 一种高能效轻量化结构双足机器人的逆运动学求解方法
CN112051741A (zh) * 2020-09-16 2020-12-08 北京理工大学 一种双足机器人动态运动生成与控制方法
US20220203526A1 (en) * 2020-12-30 2022-06-30 Ubtech Robotics Corp Ltd Humanoid robot control method, humanoid robot using the same, and computer readable storage medium
CN112882485A (zh) * 2021-01-25 2021-06-01 河南大学 一种双足机器人逆运动学的几何计算方法
WO2022227457A1 (zh) * 2021-04-30 2022-11-03 深圳市优必选科技股份有限公司 质心位姿估计方法, 装置, 计算机可读存储介质及机器人
CN113835429A (zh) * 2021-09-08 2021-12-24 南方科技大学 一种仿生双足机器人的控制方法及仿生双足机器人

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LI, SHIQI 等: "Human-like redundancy resolution: An integrated inverse kinematics scheme for anthropomorphic manipulators with radial elbow offset", 《ADVANCED ENGINEERING INFORMATICS》 *
周江琛 等: "基于变质心高度策略的双足机器人变步长步态规划", 《中南大学学报(自然科学版)》 *
朱雅光: "基于阻抗控制的多足步行机器人腿部柔顺控制研究", 《CNKI》 *

Also Published As

Publication number Publication date
CN114355964B (zh) 2023-08-18

Similar Documents

Publication Publication Date Title
CN111558941B (zh) 浮动基动力学前馈控制方法、装置和多足机器人
Kofinas et al. Complete analytical inverse kinematics for NAO
JP6137155B2 (ja) 干渉回避方法、制御装置及びプログラム
CN113618742B (zh) 一种机器人避障方法、装置和机器人
JP4693643B2 (ja) ロボットの教示支援装置及びそのためのプログラム
KR20020086468A (ko) 각식 이동 로봇을 위한 동작 제어 장치 및 동작 제어방법, 및 로봇 장치
JP2003039353A (ja) マニピュレータの制御方法
CN113835429A (zh) 一种仿生双足机器人的控制方法及仿生双足机器人
CN113146638B (zh) 质心位姿估计方法、装置、计算机可读存储介质及机器人
US20230117928A1 (en) Nonlinear trajectory optimization for robotic devices
Serrezuela et al. A consistent methodology for the development of inverse and direct kinematics of robust industrial robots
JP2014014920A (ja) 機構制御装置、プログラム及び機構制御方法
CN116038647A (zh) 代步车、机械臂的控制方法、电子设备及存储介质
CN112757301B (zh) 机器人抗扰动控制方法、装置、电子设备及存储介质
CN112720479B (zh) 一种机器人姿态控制方法、装置及机器人
CN114355964B (zh) 多自由度单腿运动学求解方法、装置和机器人
JP6515828B2 (ja) 干渉回避方法
CN115994288B (zh) 一种双足机器人逆运动学求解方法及装置
CN113119102A (zh) 基于浮动基飞轮倒立摆的仿人机器人建模方法和装置
CN113858207B (zh) 人形机器人的倒地控制方法、装置、电子设备及存储介质
CN114326769B (zh) 机器人运动矫正方法及装置、机器人控制设备和存储介质
CN114454180A (zh) 移动机器人的运动控制方法、装置和移动机器人
CN114237233B (zh) 机器人下棋方法、装置、计算机可读存储介质及机器人
Fadli et al. Design and implementation of walking pattern and trajectory compensator of NAO humanoid robot
US11787045B2 (en) Robot with an inverse kinematics (IK)-based controller for retargeting input motions

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20231204

Address after: Room 601, 6th Floor, Building 13, No. 3 Jinghai Fifth Road, Beijing Economic and Technological Development Zone (Tongzhou), Daxing District, Beijing, 100176

Patentee after: Beijing Youbixuan Intelligent Robot Co.,Ltd.

Address before: 518000 16th and 22nd Floors, C1 Building, Nanshan Zhiyuan, 1001 Xueyuan Avenue, Nanshan District, Shenzhen City, Guangdong Province

Patentee before: Shenzhen Youbixuan Technology Co.,Ltd.

TR01 Transfer of patent right