机器人动力学 - 机器人学中的惯性矩阵坐标转换及在SolidWorks中的测量
发布日期:2021-07-01 04:00:05 浏览次数:2 分类:技术文章

本文共 4513 字,大约阅读时间需要 15 分钟。

文章目录

题记

很早就想谈谈这个话题,奈何时间不允许。最近相对比较闲,所以来侃一侃机器人学中的惯性矩阵这点事儿。

机器人学中的惯性张量

惯性张量是表述机器人本体在转动过程中状态改变的难易程度的一个量,与之对应的是平动中的质量。

对于惯性张量的表示,在不同的坐标系下其值是不同的,我们可以通过坐标旋转和平行移轴定理来将不同坐标系下的惯性张量矩阵联系在一起。

假设现有一坐标系 A {A} A,其惯性张量矩阵如下:

A I = [ I x x − I x y − I x z − I x y I y y − I y z − I x z − I y z I z z ] ^{A} \mathrm{I}=\left[\begin{array}{ccc}{I_{x x}} & {-I_{x y}} & {-I_{x z}} \\ {-I_{x y}} & {I_{y y}} & {-I_{y z}} \\ {-I_{x z}} & {-I_{y z}} & {I_{z z}}\end{array}\right] AI=IxxIxyIxzIxyIyyIyzIxzIyzIzz
w h e r e where where I x x = ∭ V ( y 2 + z 2 ) ρ d v I y y = ∭ V ( x 2 + z 2 ) ρ d v I z z = ∭ V ( x 2 + y 2 ) ρ d v I x y = ∭ V ( x y ) ρ d v I x z = ∭ V ( x z ) ρ d v I y z = ∭ V ( y z ) ρ d v \begin{aligned} I_{x x} &=\iiint_{V}\left(y^{2}+z^{2}\right) \rho d v \\ I_{y y} &=\iiint_{V}\left(x^{2}+z^{2}\right) \rho d v \\ I_{z z} &=\iiint_{V}\left(x^{2}+y^{2}\right) \rho d v \\ I_{x y} &=\iiint_{V}(x y) \rho d v \\ I_{x z} &=\iiint_{V}(x z) \rho d v \\ I_{y z} &=\iiint_{V}(y z) \rho d v \end{aligned} IxxIyyIzzIxyIxzIyz=V(y2+z2)ρdv=V(x2+z2)ρdv=V(x2+y2)ρdv=V(xy)ρdv=V(xz)ρdv=V(yz)ρdv
刚体有单元 d v dv dv组成,其密度为 ρ ρ ρ,每个单元位置由矢量 [ x , y , z ] T [x ,y, z]^T [xyz]T指定。
在这里插入图片描述

惯性张量在不同坐标系下的转换

  1. 假设坐标系 1 {1} 1的惯性张量矩阵为 I 1 I_1 I1,坐标系 2 {2} 2的惯性张量矩阵为 I 2 I_2 I2,两个坐标系原点 o {o} o重合,由坐标系 1 {1} 1变换到坐标系 2 {2} 2的矩阵为 R 12 R_{12} R12,则 I 1 I_1 I1 I 2 I_2 I2的关系如下:

    I 1 = R 12 I 2 ( R 12 ) T I_1 = R_{12} I_2 (R_{12})^T I1=R12I2(R12)T

    下面来举个栗子!

    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

  2. 假设有一输出坐标系为坐标系 1 {1} 1,质点为 c {c} c,其在坐标系 1 {1} 1中的坐标为 P c = [ x c , y c , z c ] T P_c = [x_c ,y_c, z_c]^T Pc=[xcyczc]T,对齐坐标系 1 {1} 1的质心坐标系为 C {C} C,则由平行移轴定理可得:

    I 1 = I c + m ( P c T P c I 3 × 3 − P c P c T ) I_1 = I_c + m (P_{c}^{T}P_cI_{3×3} - P_cP_{c}^{T}) I1=Ic+m(PcTPcI3×3PcPcT)

下面来举个栗子!

  • 例1 方块实体
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
  • 例2 机械臂连杆
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述
    在这里插入图片描述

注意:

  • SolidWorks测量的惯性张量值,都是绝对的,因此在写惯性张量矩阵时要在除对角线位置外添加负号!!!

几个概念

惯性张量惯性积惯性矩转动惯量,几个概念傻傻分不清……

  1. 首先我们来看下转动惯量,转动惯量(此处为质量转动惯量)只决定于刚体的形状、质量分布和转轴的位置,其表达式为:
    I = ∑ i m i r i 2 I=\sum_{i} m_{i} r_{i}^{2} I=imiri2
    假设刚体质量连续分布,那么表达式可以写成: I = ∭ V r 2 d m = ∭ V r 2 ρ d V I=\iiint_{V} r^{2} \mathrm{d} m=\iiint_{V} r^{2} \rho \mathrm{d} V I=Vr2dm=Vr2ρdV
    w h e r e where where m i m_i mi表示质元的质量, r r r表示质元到转轴的垂直距离, ρ ρ ρ为密度。
    在SI单位制中,它的单位是 k g ⋅ m 2 kg·m^2 kgm2
    在这里插入图片描述2. 对于面积转动惯量,假设转轴为 z z z,那么平面积A对z轴的转动惯量(又称极转动惯量)为:
    I z = ∫ ( x 2 + y 2 ) d A I_{z}=\int\left(x^{2}+y^{2}\right) d A Iz=(x2+y2)dA
    在SI单位制中,它的单位是 m 4 m^4 m4
    在这里插入图片描述
    3.对于惯性积,质量惯性积是刚体动力学中一个重要的质量几何性质。刚体中的质量微元 Δ m i Δm_i Δmi与这微元的两个直角坐标的乘积对刚体的总和。其表达式为:
    I x y = ∫ x y d m I_{x y}=\int x y d m Ixy=xydm
    w h e r e where where x i x_i xi y i y_i yi为刚体微元 d m dm dm x x x y y y轴上的坐标。
    在SI单位制中,它的单位是 k g ⋅ m 2 kg·m^2 kgm2
    4.面积惯性积是截面的一个重要几何性质。平面积A对评估面内互相垂直的 x x x y y y轴的惯性积为:
    I x y = ∫ x y d A I_{x y}=\int x y d A Ixy=xydA
    w h e r e where where x x x y y y为面元 d A dA dA的坐标。
    在SI单位制中,它的单位是 m 4 m^4 m4
    5.对于惯性矩,惯性矩(moment of inertia of an area)是一个几何量,通常被用作描述截面抵抗弯曲的性质。面积元素 d A dA dA与其至 z z z轴或 y y y轴距离平方的乘积 y 2 d A y^2dA y2dA z 2 d A z^2dA z2dA,分别称为该面积元素对于 z z z轴或 y y y轴的惯性矩或截面二次轴矩。
    z z z轴惯性矩: I Z = ∫ A y 2 d A I_{Z}=\int_{A} y^{2} d A IZ=Ay2dA
    y y y轴惯性矩: I y = ∫ A z 2 d A I_{y}=\int_{A} z^{2} d A Iy=Az2dA
    对任意一对互相垂直轴的惯性矩之和,等于截面对该二轴交点的极惯性矩: I p = ∫ A ρ 2 d A I_{p}=\int_{A} \rho^{2} d A Ip=Aρ2dA
    在这里插入图片描述
  • 跟质量相关的都是空间上的量,而跟面积相关的只是截面上的量。前者描述了物体转动的难易程度,或者说是抵抗物体发生位移(对于旋转来说)的能力(刚体力学量),对于平动中该量对应的是质量,后者表述的是物体截面抵抗外力变形的能力(材料力学中的刚度,是弹性力学量)。

从以上搬的砖的内容来看,惯性矩与惯性积是一个东西,极惯性矩与面积转动惯量相同。

  • 二者之间的联系在于:面积惯性矩跟质量无关, I = ∫ r i 2 d A I=\int r_i^2 dA I=ri2dA,如果引入材料的密度 ρ ρ ρ,则该式变为: m I = ∫ r i 2 d m mI=\int r_i^2 dm mI=ri2dm,即为转动惯量的积分表达式。这说明:惯性矩可以表征材料转动的难以程度,这种局部转动的宏观表现即为弯曲变形和扭转变形;惯性矩仅体现出几何形状对物体抵抗自身变形的影响,而转动惯量同时体现质量分布和几何形状的影响。

博主瞎扯淡(猜错请批评指正)

各种文献中的叫法都不一样,从看过的文献中,博主猜测:

  1. 惯性矩和惯性积是描述截面抗形变能力的,惯性矩是类似于 I y = ∫ x 2 d A I_{y}=\int x^2dA Iy=x2dA I x = ∫ y 2 d A I_{x}=\int y^2dA Ix=y2dA这种的, 这是极惯性矩 I z = ∫ ( x 2 + y 2 ) d A I_{z}=\int\left(x^{2}+y^{2}\right) d A Iz=(x2+y2)dA,惯性积的形式为 I x y = ∫ x y d A I_{xy}=\int xydA Ixy=xydA

在这里插入图片描述

2. 惯量矩和惯量积是描述转动难易程度,惯量矩是类似于 I x x = ∫ ( x 2 + y 2 ) d m I_{xx}=\int (x^2+y^2)dm Ixx=(x2+y2)dm I y y = ∫ ( x 2 + z 2 ) d m I_{yy}=\int (x^2+z^2)dm Iyy=(x2+z2)dm这种的,惯量积的形式为 I x y = ∫ x y d m I_{xy}=\int xydm Ixy=xydm
在这里插入图片描述

总结:

  • 叫“性”的跟截面积有关,叫“量”的跟质量有关
  • 叫“矩”的跟单轴或者两轴的平方和相关,叫“积”的跟两轴乘积有关

参考文献:

https://baike.baidu.com/item/%E6%83%AF%E6%80%A7%E7%A7%AF/2294932?fr=aladdin

https://baike.baidu.com/item/%E6%83%AF%E6%80%A7%E7%9F%A9/8155407?fr=aladdin

https://baike.baidu.com/item/%E6%83%AF%E6%80%A7%E7%A7%AF/2294932?fr=aladdin

https://baike.baidu.com/item/%E6%83%AF%E6%80%A7%E5%BC%A0%E9%87%8F/5322910?fr=aladdin

http://muchong.com/html/201311/6637340.html

转载地址:https://miracle.blog.csdn.net/article/details/103579030 如侵犯您的版权,请留言回复原文章的地址,我们会给您删除此文章,给您带来不便请您谅解!

上一篇:SolidWorks转Simscape Multibody(两种方案)
下一篇:MATLAB - 符号运算(持续更新)

发表评论

最新留言

哈哈,博客排版真的漂亮呢~
[***.90.31.176]2024年05月05日 19时42分29秒