CN107562968B - 一种足式机器人动力学建模的混合计算方法 - Google Patents
一种足式机器人动力学建模的混合计算方法 Download PDFInfo
- Publication number
- CN107562968B CN107562968B CN201610503559.5A CN201610503559A CN107562968B CN 107562968 B CN107562968 B CN 107562968B CN 201610503559 A CN201610503559 A CN 201610503559A CN 107562968 B CN107562968 B CN 107562968B
- Authority
- CN
- China
- Prior art keywords
- matrix
- robot
- physical parameters
- mass
- center
- 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.)
- Active
Links
Images
Abstract
本发明提供一种足式机器人动力学建模的混合计算方法,通过预先存储的足式机器人的基本物理参数,对足式机器人各个关节的角度和角速度进行采样,获得采样反馈的采样数值;根据采样数值计算出坐标旋转的转移矩阵;根据基本物理参数,根据拉格朗日方法中对惯性矩阵的计算方法,对每个连杆的惯性矩阵进行推导,通过矩阵求和的方式获得足式机器人的惯性矩阵;根据采样数值、基本物理参数,在牛顿欧拉方法的基础上进行改进,通过迭代计算,获得足式机器人的离心力、科里奥利力和重力的相关向量,以完成足式机器人的动力学建模,综合拉格朗日方法和牛顿欧拉方法各自的计算优点并在不同阶段加以修改,提高动力学建模计算效率。
Description
技术领域
本发明涉及足式机器人动力学建模计算领域,具体涉及一种足式机器人动力学建模的混合计算方法。
背景技术
自从机器人诞生以来,人类的生活、工作方式发生了巨大的变化,机器人在各种场合辅助人类完成复杂的、重复的工作。近年来,国内外掀起了研究仿人机器人的热潮,其研究内容主要包括稳定行走、灵巧作业、人机交互等方面。其中机器人在运动过程中能否稳定的行走取决于其机器人动力学建模方法的是否高效和便于分析。
目前,常用的机器人动力学建模方法主要分为两种,即拉格朗日法和牛顿欧拉法。其中,拉格朗日法是应用系统的动能和势能,通过构建拉格朗日函数并结合虚位移的方法建立动力学模型,由于所构建的动力学模型具有清晰的数学表达式,所以非常便于对整个系统的动力学分析,然而当系统的结构变得十分复杂时,其构建的动力学模型表达式中的各个矩阵变得非常复杂,从而计算效率很低;而牛顿欧拉法则是以牛顿第二定律和转矩对于角动量的作用为基础,通过迭代计算的方法来构建动力学模型,对于复杂的机器人系统而言,虽然采用了迭代的计算方法,使得计算效率大大提高,但是由于得到的模型形式不够清晰,致使此方法不适合于进行动力学的相关分析。
随着近些年人们对机器人控制理论的不断创新和探索,寻求一种既适合于动力学分析又拥有高效率计算的动力学建模方法已经迫在眉睫并具有较高的应用价值与推广意义。
发明内容
本发明针对现有足式机器人的动力学建模方法的处理效率低,计算复杂等问题,提供了一种足式机器人动力学建模的混合计算方法。
一方面,本发明提供一种足式机器人动力学建模的混合计算方法,
足式机器人的动力学模型的数学形式如下:
混合计算方法包括:
预先存储所述足式机器人的基本物理参数;
对所述足式机器人各个关节的角度和角速度进行采样,获得采样反馈的采样数值;
根据所述采样数值计算出坐标旋转的转移矩阵;
根据所述基本物理参数,结合拉格朗日方法中对惯性矩阵的计算方法,对所述足式机器人各个连杆的惯性矩阵进行计算,通过矩阵求和运算获得所述足式机器人的惯性矩阵D(q);
作为一种可选的实施方式,结合拉格朗日方法中对惯性矩阵的计算方法,对所述足式机器人各个连杆的惯性矩阵进行计算,通过矩阵求和运算获得所述足式机器人的惯性矩阵D(q),包括:
根据所述基本物理参数和所述坐标旋转的转移矩阵,计算出所述各个连杆的质心线速度的雅克比矩阵,和所述各个连杆的质心角速度的雅克比矩阵;
根据所述基本物理参数,所述各个连杆的质心的线速度的雅克比矩阵,以及所述各个连杆的质心角速度的雅克比矩阵,经过矩阵运算得到所述足式机器人各个连杆的惯性矩阵;再通过矩阵求和的方式计算出所述机器人的惯性矩阵,所述拉格朗日方法为应用所述惯性矩阵运算。
作为一种可选的实施方式,根据所述采样数值,所述基本物理参数,对牛顿欧拉方法进行改进,通过迭代计算方式,获得所述足式机器人的所述离心力、科里奥利力和重力的相关向量,以完成所述足式机器人的动力学建模,包括:
根据所述采样数值,所述基本物理参数,以及所述坐标旋转的转移矩阵,计算出所述各个连杆的质心的线速度和所述各个关节的连杆围绕质心的角速度;
根据所述基本物理参数,所述坐标旋转的转移矩阵,所述各个连杆的所述质心的线速度,以及所述各个关节的连杆围绕质心的角速度,计算出所述各个关节的质心的线加速度中的陀螺效应部分,和所述各个关节的所述质心的角加速度中的陀螺效应部分;
根据所述基本物理参数,所述坐标旋转的转移矩阵,所述各个连杆的所述质心的线加速度,以及所述各个连杆的所述质心的角加速度,计算出所述各个关节的所受的作用力,和所述各个关节的所受的作用转矩;
根据所述坐标旋转的转移矩阵,所述各个关节所受的作用力,以及所述各个关节所受的作用转矩的计算,去映射到所述足式机器人的所述各个关节上,计算出所述足式机器人的所述各个关节的转矩,最终等效出所述足式机器人的离心力、科里奥利力和重力相关向量,以完成所述足式机器人的动力学建模;所述对牛顿欧拉方法进行改进为应用所述离心力、科里奥利力和重力相关向量的运算。
本发明提供一种足式机器人动力学建模的混合计算方法,通过预先存储的所述足式机器人的基本物理参数;对所述足式机器人各个关节的角度和角速度进行采样,获得采样反馈的采样数值;根据所述采样数值计算出坐标旋转的转移矩阵;根据所述基本物理参数,结合拉格朗日方法中惯性矩阵的计算方法,对机器人各个连杆的惯性矩阵进行计算,并根据矩阵求和的方式获得所述足式机器人的惯性矩阵;根据所述采样数值、所述基本物理参数,对牛顿欧拉方法进行改进,通过迭代计算,获得所述足式机器人的离心力、科里奥利力和重力的相关向量,以完成所述足式机器人的动力学建模,结合拉格朗日方法和牛顿欧拉方法的计算优点,统一计算获得足式机器人的动力学建模,综合了拉格朗日方法和牛顿欧拉方法的计算优点,使得动力学建模过程的计算效率大大提高,计算出的模型更加稳定,计算步骤更加简化,更加便于建模分析。
附图说明
图1为本发明实施例提供的一种足式机器人动力学建模的混合计算方法流程图;
图2为本发明实施例提供的另一种足式机器人动力学建模的混合计算方法流程图;
图3为本发明实施例提供的另一种足式机器人动力学建模的混合计算方法流程图。
具体实施方式
下面阐述的实施例代表允许本领域技术人员实践本发明的必要信息,并且示出实践本发明的最佳方式。一旦根据附图阅读了以下的描述,本领域技术人员就将理解本发明的构思并且将认识到此处未特别阐明的这些构思的应用。应当理解,这些构思和应用落入本公开和所附权利要求书的范围。下面结合实施例对本发明进一步说明。
请参见图1,图1为本发明实施例提供的一种足式机器人动力学建模的混合计算方法流程图。如图1所示,本实施例提供的足式机器人动力学建模的混合计算方法包括:
110、预先存储足式机器人的基本物理参数;
120、对足式机器人各个关节的角度和角速度进行采样,获得采样反馈的采样数值;
130、根据采样数值计算出坐标旋转的转移矩阵;
140、根据基本物理参数,结合拉格朗日方法对惯性矩阵的计算,推导出各个连杆的惯性矩阵,通过矩阵求和的方式获得足式机器人的惯性矩阵;
其中,足式机器人的动力学模型的数学形式如下:
作为一种可选的实施方式,请参见图2,图2为本发明实施例提供的另一种足式机器人动力学建模的混合计算方法流程图,如图2所示,基于图1所示的步骤140进一步包括以下步骤,即根据基本物理参数,以及结合拉格朗日方法中对惯性矩阵的计算方法,对足式机器人各个连杆的惯性矩阵进行计算,并通过矩阵求和运算,获得足式机器人的惯性矩阵,包括:
141、根据基本物理参数和坐标旋转的转移矩阵,计算出各个连杆质心的线速度的雅克比矩阵,以及各个连杆质心角速度的雅克比矩阵;
142、根据基本物理参数,各个连杆的质心的线速度的雅克比矩阵,以及各个连杆的质心角速度的雅克比矩阵,经过矩阵运算得到足式机器人的各个连杆的惯性矩阵,通过矩阵求和运算得到整个机器人系统的惯性矩阵D(q)。在本实施例中,拉格朗日方法为应用计算机器人的惯性矩阵。
作为一种可选的实施方式,请参见图3,图3为本发明实施例提供的另一种足式机器人动力学建模的混合计算方法流程图,如图3所示,基于图2所示的步骤150进一步包括以下步骤,即根据采样数值、基本物理参数以及改进牛顿欧拉方法,通过迭代计算,获得足式机器人的离心力、科里奥利力和重力的相关向量,以完成足式机器人的动力学建模,包括:
151、根据采样数值、基本物理参数和坐标旋转的转移矩阵,计算出各个连杆的质心的线速度和各个关节的连杆围绕质心的角速度;
152、根据基本物理参数、坐标旋转的转移矩阵、各个连杆的质心的线速度和各个关节的连杆围绕质心的角速度,计算出各个连杆的质心的线加速度的陀螺效应部分和各个连杆的质心的角加速度的陀螺效应部分;
153、根据基本物理参数、坐标旋转的转移矩阵、各个关节的质心的线加速度和各个关节的质心的角加速度,计算出各个关节的质心所受的作用力和各个关节的质心所受的作用转矩;
154、根据坐标旋转的转移矩阵、各个关节所受的作用力和各个关节所受的作用转矩,去映射到足式机器人的各个关节上,从而计算出足式机器人的各个关节的转矩,最终由此等效出足式机器人的离心力、科里奥利力和重力相关向量,以完成足式机器人的动力学建模。
在本实施例中,对牛顿欧拉方法进行改进为应用计算离心力、科里奥利力和重力相关向量。
综上描述,本发明实施例提供的足式机器人动力学建模的混合计算方法,通过预先存储的足式机器人的基本物理参数;对足式机器人各个关节的角度和角速度进行采样,获得采样反馈的采样数值;根据采样数值计算出坐标旋转的转移矩阵;根据基本物理参数,以及结合拉格朗日方法中惯性矩阵的计算方法,对机器人各个连杆的惯性矩阵进行计算,通过矩阵求和的方式获得足式机器人的惯性矩阵;根据采样数值、基本物理参数以及改进牛顿欧拉方法,通过迭代计算,获得足式机器人的离心力、科里奥利力和重力的相关向量,以完成足式机器人的动力学建模,结合拉格朗日方法和牛顿欧拉方法的计算优点,统一计算获得足式机器人的动力学建模,综合了拉格朗日方法和牛顿欧拉方法的计算优点,使得动力学建模过程的计算效率大大提高,计算出的模型更加稳定,计算步骤更加简化,更加便于建模分析。
显然,上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而这些属于本发明的精神所引伸出的显而易见的变化或变动仍处于本发明的保护范围之中。
Claims (1)
1.一种足式机器人动力学建模的混合计算方法,其特征在于,
所述足式机器人的动力学模型的数学形式如下:
所述混合计算方法包括:
预先存储所述足式机器人的基本物理参数,对所述足式机器人各个关节的角度和角速度进行采样,获得采样反馈的采样数值;
根据所述采样数值计算出坐标旋转的转移矩阵;
根据所述基本物理参数,结合拉格朗日方法中对惯性矩阵的计算方法,对所述足式机器人各个连杆的惯性矩阵进行计算,通过矩阵求和运算获得所述足式机器人的惯性矩阵D(q),具体步骤包括:
根据所述基本物理参数和所述坐标旋转的转移矩阵,计算出所述各个连杆的质心线速度的雅克比矩阵,和所述各个连杆的质心角速度的雅克比矩阵;
根据所述基本物理参数,所述各个连杆的质心的线速度的雅克比矩阵,以及所述各个连杆的质心角速度的雅克比矩阵,经过矩阵运算得到所述足式机器人各个连杆的惯性矩阵;通过矩阵求和的方式得到机器人的惯性矩阵D(q),所述拉格朗日方法为应用所述转移矩阵进行矩阵运算;
根据所述采样数值,所述基本物理参数,以及所述坐标旋转的转移矩阵,计算出所述各个连杆的质心的线速度以及所述各个关节的连杆围绕所述质心的角速度;
根据所述基本物理参数,所述坐标旋转的转移矩阵,所述各个连杆的所述质心的线速度,以及所述各个关节的连杆围绕质心的角速度,计算出所述各个连杆的质心的线加速度中的陀螺效应部分,和所述各个连杆的质心的角加速度中的陀螺效应部分;
根据所述基本物理参数,所述坐标旋转的转移矩阵,所述各个连杆的质心的加速度,以及所述各个关节的质心的角加速度,计算出所述各个关节所受的作用力,和所述各个关节所受的作用转矩;
根据所述坐标旋转的转移矩阵,所述各个关节所受的作用力,以及所述各个关节所受的作用转矩的计算逻辑,去映射到所述足式机器人的各个关节上,计算出所述足式机器人的所述各个关节的转矩,最终等效出所述足式机器人的离心力、科里奥利力和重力的相关向量以完成所述足式机器人的动力学建模;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610503559.5A CN107562968B (zh) | 2016-06-30 | 2016-06-30 | 一种足式机器人动力学建模的混合计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610503559.5A CN107562968B (zh) | 2016-06-30 | 2016-06-30 | 一种足式机器人动力学建模的混合计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107562968A CN107562968A (zh) | 2018-01-09 |
CN107562968B true CN107562968B (zh) | 2021-01-22 |
Family
ID=60969040
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610503559.5A Active CN107562968B (zh) | 2016-06-30 | 2016-06-30 | 一种足式机器人动力学建模的混合计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107562968B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110850834B (zh) * | 2019-12-02 | 2021-08-03 | 合肥工业大学 | 并联机器人的建模方法、建模系统、控制方法及控制系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103426196A (zh) * | 2013-08-30 | 2013-12-04 | 电子科技大学 | 一种流体环境下的关节动画建模技术 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101439929B1 (ko) * | 2012-04-30 | 2014-09-17 | 장성훈 | 로봇의 삼차원 시뮬레이션 방법 및 시스템 |
-
2016
- 2016-06-30 CN CN201610503559.5A patent/CN107562968B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103426196A (zh) * | 2013-08-30 | 2013-12-04 | 电子科技大学 | 一种流体环境下的关节动画建模技术 |
Non-Patent Citations (2)
Title |
---|
Dynamic Modeling of Biped Robot using Lagrangian and Recursive Newton-Euler Formulations;Hayder F. N. Al-Shuka 等;《International Journal of Computer Applications》;20140930;全文 * |
Gough-Stewart;吕帮俊 等;《华中科 大学学报(自然科学版)》;全文;20111031;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107562968A (zh) | 2018-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110315543B (zh) | 一种双足机器人步态生成与优化方法 | |
US11579591B2 (en) | Feedforward control method for floating base dynamics, computer-readable storage medium and robot | |
JP7324932B2 (ja) | 動的計画コントローラ | |
Wang et al. | A UKF-based predictable SVR learning controller for biped walking | |
Irawan et al. | Force threshold-based omni-directional movement for hexapod robot walking on uneven terrain | |
CN107562968B (zh) | 一种足式机器人动力学建模的混合计算方法 | |
Wang et al. | Realization of a real-time optimal control strategy to stabilize a falling humanoid robot with hand contact | |
US20230191604A1 (en) | Robot stability control method, robot and computer-readable storage medium | |
Mahapatra et al. | Energy-efficient inverse dynamic model of a Hexapod robot | |
Ding et al. | Dynamic modeling and locomotion control for quadruped robots based on center of inertia on SE (3) | |
Wang et al. | The simulation of nonlinear model predictive control for a human-following mobile robot | |
Xie et al. | Gait optimization and energy-based stability for biped locomotion using large-scale programming | |
Lee et al. | Development and motion control of a biped walking robot based on passive walking theory | |
Zhang et al. | Optimal energy gait planning for humanoid robot using geodesics | |
Chen et al. | Optimized 3D stable walking of a bipedal robot with line-shaped massless feet and sagittal underactuation | |
Carnier et al. | Pseudospectral Optimization of Planar Robotic Walking using Whole-body Dynamics | |
Liu et al. | Deep deterministic policy gradient with constraints for gait optimisation of biped robots | |
Yilmaz et al. | Circular arc-shaped walking trajectory generation for bipedal humanoid robots | |
JP4457617B2 (ja) | モーション生成システム及びロボットのモーション評価方法 | |
WO2024021767A1 (zh) | 控制足式机器人的方法、装置、设备、足式机器人、计算机可读存储介质和计算机程序产品 | |
WO2024021744A1 (zh) | 一种控制足式机器人的方法、装置、电子设备、计算机可读存储介质、计算机程序产品及足式机器人 | |
Liu et al. | Dynamic Modeling of Quadrupedal Robot Based on the Screw Theory | |
Tan et al. | Dynamics Modeling and Gait Analysis of Transport-Oriented Hexapod Robot | |
US20230030054A1 (en) | Method and apparatus for controlling motion of legged robot, device, and medium | |
Roy et al. | Study on energy consumption in turning motion of hexapod walking robots |
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 |