第5章--典型机械系统的建模PPT课件.ppt

控制相关专业研究生选修课程 系统建模方法 马宏军 东北大学信息学院控制理论与导航技术研究所 2013年3月 逸夫楼203 1 第5章典型机械系统的建模 机械系统遍及工程技术和社会各个领域 除机械设备与装置外 还是构成其他复杂系统的基础和基本环节 如控制系统地执行机构 飞机舵面传动装置 导弹发射架 飞行模拟器的运动平台等 这些系统建模目标多是建立选定参考坐标系下的系统运动方程和动力学方程 属于 白箱 问题 因此 采用的建模方法不外乎是机理分析法或图解法 对复杂的机械系统还可能应用辨识方法 在建模中 主要将利用牛顿力学定律 拉格朗日函数 并结合能量守恒原理及有关近似理论等 针对特殊的机械系统 机器人 其运动学及动力学分析的数学建模和仿真与传统的机械动态特性研究因其多运动自由度特点 多体动力学理论基础在机器人运动动力学分析中特别适用 2 5 1基于力学理论的机械系统建模 一 空间任意力系的平衡方程 由理论力学可知 空间任意力系平衡的必要和充分条件是 力系中所有各力在三坐标轴中每一轴上的投影和分别等于零 又这些力对于这些轴的力矩的代数和也分别等于零 其数学表达式为 3 二 牛顿第二定律数学表达式 牛顿第二定律告诉我们 物体受外力作用时 所获得的加速度大小与合力大小成正比 与物体的质量成反比 加速度的方向与合外力的方向相同 其数学表达式为 4 例5 1测量转动惯量实验装置如右图一个转动物体 它的质量为m 由两根垂直的绳索 无弹性 挂起 每根绳索的长度为h 绳索相距为2a 重心位于通过连接绳索两点的中点的垂线上 假设物体绕通过重心的垂直轴转一个小的角度 然后释放 求摆动周期T 物体通过重心的垂直轴转的转动惯量J 假设物体绕通过重心的垂直轴转一个小的角度时 夹角和夹角间存在下列关系 因此 注意 每根绳索的受力F的垂直分量等于mg 2 F的水平分量为mg 2 两根绳索的F的水平分量产生扭矩mga使物体转动 因此 摆动的运动方程为 5 或写成 由此求得摆动周期为 得到转动惯量J 1 隐含的假定 2 系统的阶次 3 模型是否合理 利用常识判断 4 非线性的情况 利用计算机辅助判断 5 系统参数的获取 6 例5 2单摆系统下图所示的单摆系统为输入力矩 为输出摆角 m为小球质量 L为摆长 根据力系平衡建立系统方程 这是一个非线性方程 根据Taylor级数展开得 当很小时 高阶小数可以忽略 则 非线性系统方程可简化成线性系统方程 第一次作业 7 例5 3设一个弹簧 质量 阻尼系统安装在一个不计质量的小车上 如下图所示 推导系统数学模型 假设t 0时小车静止不动 并且安装在小车上的系统也处于静止状态 在这个系统中 u t 是小车的位移 并且是系统的输入量 不计小车的质量 得到 1 验证 阶次 稳定性等 2 拓展 非线性 忽略因素 3 参数确定 白箱实验 时域频域指标反解 8 例5 4有一质量 弹簧 阻尼系统如图所示 运用力学方法建立该系统的数学模型 系统图 力分解图 根据力平衡原理 建立系统方程 9 10 例5 5机械式加速度计 下图给出机械式加速度计测量悬浮试验橇加速度的示意图 试验橇采取磁悬浮方式以较小的高度e悬浮于导轨上方 由于质量M相对于及速度计箱体的位移y与箱体的 即试验橇的 加速度成正比 因而加速度计能测得试验橇的加速度 我们的目的是设计一个具有合理动态响应的加速度计 它能在可以接受的时间内测得所需要的特征量 y t qa t q为常数 11 分析质量M的受力情况 我们有 或 或 12 例5 6倒立摆系统 左下图为人手保持倒摆平衡的问题 相应的平衡条件为 右下图表示的是小车上的倒摆控制问题 小车必须处于运动状态才能保持质量m始终处于小车上方 系统状态变量应当与旋转角以及小车的位移有关 人手到立摆的平衡小车和倒摆 13 设M m 旋转角 足够小 于是可以对运动方程做线性近似处理 这样 系统水平方向受力之和将为 其中 u t 等于施加在小车上的外力 l是质量到铰接点的距离 铰接点处的转矩之和为 选定两个2阶系统的状态变量为 将a b两式写成状态变量的形式 可得 a b c d 14 为得到1阶微分方程组 解出式 d 中的 代入式 c 并注意到M m 则有 e 再解出式 c 中的 并代入式 d 可得 于是 4个1阶微分方程为 15 系统状态方程则为 1 验证 阶次 稳定性等 2 拓展 非线性 忽略因素 3 参数确定 白箱实验 时域频域指标反解 16 5 2能量法推导运动方程 一 功 能 功率如果力被认为是努力的度量 那么功就是成就的度量 而能量就是做功的能力 功的概念没有考虑时间的因素 就要引入功率的概念 功机械系统中的功等于力与力作用的距离的乘积 或力矩与角位移的乘积 力与距离要在同一方向上度量 设力F作用于a至b连接路径中运动的质点m上 那么F所作的功可一般描述为 能量一般情况下 能量可以定义为做功的能力 机械系统中能有势能和动能两种形式 功率是做功的速率 即 dW表示在dt时间间隔内所作的功 17 二 能量法推导运动方程 能量法推导运动方程的根本就是能量守恒定律 如果系统没有能量输入和输出 我们从系统总能量保持相等这一事实出发来推导运动方程 例5 7如右图表示一个半径为R 质量为m的均质圆柱体 它可以绕其转轴自由转动并通过一个弹簧与墙壁连接 假设圆柱体纯滚动而无滑动 求系统的动能和势能并导出系统运动方程 圆柱体的动能等于质心移动动能和绕质心转动的动能之和 系统由于弹簧变形所产生的势能为 系统总能量为 18 考虑到圆柱体做无滑动的滚动 因此 并且注意到转动惯量J等于 我们得到 考虑到能量守恒定律 总能量为常数 即总能量导数为零 得到 注意到 并不总为0 因此必须恒等于0 即 如果将以上方程转为转动运动 只要把代入得到 能量不守恒怎么办 19 经典动力学的两个发展方面 拓宽研究领域 矢量动力学又称为牛顿 欧拉动力学 牛顿运动定律由单个自由质点 受约束质点和质点系 以达朗贝尔原理为基础 欧拉将牛顿运动定律 刚体和理想流体 寻求新的表达形式 将虚位移原理和达朗贝尔原理综合应用于动力学 建立分析力学的新体系 拉格朗日力学 20 考察由N个质点的 具有理想约束的系统 根据达朗贝尔原理 有 令系统有任意一组虚位移 系统的总虚功为 动力学普遍方程 21 系统的总虚功为 利用理想约束条件 得到 动力学普遍方程 任意瞬时作用于具有理想的系统上的主动力与惯性力在系统的任意虚位移上的元功之和等于零 22 动力学普遍方程的直角坐标形式 动力学普遍方程适用于具有理想约束的系统 动力学普遍方程既适用于具有定常约束的系统 也适用于具有非定常约束的系统 动力学普遍方程既适用于具有完整约束的系统 也适用于具有非完整约束的系统 动力学普遍方程既适用于具有有势力的系统 也适用于具有无势力的系统 23 动力学普遍方程主要应用于求解动力学第二类问题 即 已知主动力求系统的运动规律 应用动力学普遍方程求解系统运动规律时 重要的是正确分析运动 并在系统上施加惯性力 由于动力学普遍方程中不包含约束力 因此 不需要解除约束 也不需要将系统拆开 应用动力学普遍方程 需要正确分析主动力和惯性力作用点的虚位移 并正确计算相应的虚功 动力学普遍方程的应用 24 拉格朗日 Lagrange 方程 主动力 虚位移 广义坐标 第i个质点的位矢 由动力学普遍方程 得 广义力 25 26 第一个Lagrange经典关系 消点 27 对任意一个广义坐标qj求偏导数 如果将位矢对任意一个广义坐标qj求偏导数 再对时间求导数 则得到 第二个拉格朗日关系式 28 29 此即拉格朗日方程 或称为第二类拉格朗日方程 如果作用在系统上的主动力都是有势力 根据有势力的广义主动力 30 引入拉格朗日函数 L T V 得到主动力为有势力的拉格朗日方程 31 对于只具有完整约束 自由度为N的系统 可以得到由N个拉格朗日方程组成的方程组 应用拉格朗日方程 一般应遵循以下步骤 首先 要判断约束性质是否完整 主动力是否有势 决定采用哪一种形式的拉格朗日方程 其次 要确定系统的自由度 选择合适的广义坐标 按照所选择的广义坐标 写出系统的动能 势能或广义力 将动能或拉格朗日函数 广义力代入拉格朗日方程 拉格朗日方程的应用 32 5 3拉格朗日方程 多自由度系统 将作为n个自由度系统的一套广义坐标 系统的运动由n个微分方程表示 其中广义坐标是因变量 时间为自变量 令作为系统在任意瞬时的势能 令作为系统在同瞬时的动能 拉格朗日函数定义为 设广义坐标是独立的 令是广义坐标的变分 非保守力 外力和摩擦力等 在广义坐标上的虚功可以写成 拉格朗日方程为 33 例5 8系统如图所示 运用拉格朗日方程建立该系统的数学模型 解 选择y1 y2为广义坐标系 其系统动能和势能分别为 34 35 例5 9某行星滚动机构中有一质量为m 半径为r的实心圆柱在半径为R 质量为M的圆筒内无滑动地滚动 已知圆柱和圆筒对轴心O的转动惯量分别为 圆柱对轴心O 的转动惯量为 建立圆筒绕其轴心转动时 该系统运动数学模型 分析 该系统为两自由度系统 取广义坐标分别为圆筒转角 和圆柱轴心偏离角 由于圆柱与圆筒间的运动是无滑动纯滚动 故在接触点A处它们具有相同的线速度 系统动能T为圆柱滚动和圆筒转动所具有的动能 36 系统的动力为重力 圆筒的势能等于零 则系统的势能为 于是有拉格朗日函数 代入拉格朗日方程有 即为该行星滚动机构的运动数学模型 37 例5 10用拉格朗日方程建立图示系统运动的微分方程 用 1 2和x作为广义坐标 以矩阵的形式写出微分方程 解 系统在任意时刻的动能为 系统在同一时刻的势能为 拉格朗日函数为 38 利用拉格朗日方程可得 39 5 4机器人静力分析与动力学 计算机技术的不断进步和发展使机器人技术的发展一次次达到一个新水平 上至太空舱 宇宙飞船 下至微机器人 深海开发 机器人技术已拓展到全球经济发展的诸多领域 成为高科技中极为重要的组成部分 人类文明的发展 科技的进步已和机器人的研究 应用产生了密不可分的关系 人类社会的发展已离不开机器人技术 而机器人技术的进步又对推动科技发展起着不可替代的作用 40 18世纪瑞士的写字偶人 哈工大爬壁机器人 爬缆索机器人 仿人机器人 北航仿生鱼 管道机器人 41 排雷机器人 索杰纳 火星车引导机器人 工业机器人 42 机器人 特别是其中最有代表性的关节型机器人 实质上是由一系列关节连接而成的空间连杆开式链机构 要研究机器人 就必须对其运动学和动力学有一个基本的了解 稳态下研究的机器人运动学分析只限于静态位置问题的讨论 未涉及机器人运动的力 速度 加速度等动态过程 实际上 机器人是一个复杂的动力学系统 机器人系统在外载荷和关节驱动力矩 驱动力 的作用下将取得静力平衡 在关节驱动力矩 驱动力 的作用下将发生运动变化 机器人的动态性能不仅与运动学因素有关 还与机器人的结构形式 质量分布 执行机构的位置 传动装置等对动力学产生重要影响的因素有关 43 机器人动力学主要研究机器人运动和受力之间的关系 目的是对机器人进行控制 优化设计和仿真 机器人动力学主要解决动力学正问题和逆问题两类问题 动力学正问题是根据各关节的驱动力 或力矩 求解机器人的运动 关节位移 速度和加速度 主要用于机器人的仿真 动力学逆问题是已知机器人关节的位移 速度和加速度 求解所需要的关节力 或力矩 是实时控制的需要 本节首先通过实例介绍与机器人速度和静力有关的雅可比矩阵 在机器人雅可比矩阵分析的基础上进行机器人的静力分析 讨论动力学的基本问题 对机器人的动态特性作简要论述 以便为机器人编程 控制等打下基础 44 一 机器人雅可比矩阵 机器人雅可比矩阵 简称雅可比 揭示了操作空间与关节空间的映射关系 雅可比不仅表示操作空间与关节空间的速度映射关系 也表示二者之间力的传递关系 为确定机器人的静态关节力矩以及不同坐标系间速度 加速度和静力的变换提供了便捷的方法 1 机器人雅可比的定义 在机器人学中 雅可比是一个把关节速度向量变换为手爪相对基坐标的广义速度向量v的变换矩阵 在机器人速度分析和静力分析中都将用到雅可比 现通过一个例子来说明 45 下图为二自由度平面关节型机器人 2R机器人 端点位置X Y与关节 1 2的关系为 即 图5 1二自由度平面关节型机器人简图 将其微分得 即 46 令 可将上式简写为 J称为图示2R机器人的速度雅可比 它反映了关节空间微小运动d 与手部作业空间微小位移dX的关系 若对式J进行运算 则图示2R机器人的雅可比可写为 从J中元素的组成可见 J阵的值是关于 1及 2的函数 5 2 5 1 47 推而广之 对于n自由度机器人 关节变量可用广义关节变量q表示 当关节为转动关节时 当关节为移动关节时 反映了关节空间的微小运动 机器人末端在操作空间的位置和方位可用末端手爪的位姿X表示 它是关节变量的函数 X X q 并且是一个6维列矢量 反映了操作空间的微小运动 它由机器人末端微小线位移和微小角位移 微小转动 组成 因此 式5 1可写为 式中 J q 是6 n维偏导数矩阵 称为n自由度机器人速度雅可比 可表示为 5 3 48 5 4 49 2 机器人速度分析 利用机器人速度雅可比可对机器人进行速度分析 对式 5 3 左 右两边各除以dt得 5 5 5 6 式中 v为机器人末端在操作空间中的广义速度 为机器人关节在关节空间中的关节速度 J q 为确定关节空间速度与操作空间速度v之间关系的雅可比矩阵 50 式中 右边第一项表示仅由第一个关节运动引起的端点速度 右边第二项表示仅由第二个关节运动引起的端点速度 总的端点速度为这两个速度矢量的合成 因此 机器人速度雅可比的每一列表示其他关节不动而某一关节运动产生的端点速度 对于图5 1所示2R机器人而言 J q 是式 5 2 所示的2 2矩阵 若令J1 J2分别为式 5 2 所示雅可比的第1列矢量和第2列矢量 则式 5 6 可写为 51 式中 J 1称为机器人逆速度雅可比 图5 1所示二自由度机器人手部的速度为 假如已知的及是时间的函数 即 则可求出该机器人手部在某一时刻的速度v f t 即手部瞬时速度 反之 假如给定机器人手部速度 可由式 5 6 解出相应的关节速度为 5 7 52 例5 11如图5 2所示的二自由度机械手 手部沿固定坐标系X0轴正向以1 0m s的速度移动 杆长l1 l2 0 5m 设在某瞬时 1 30 2 60 求相应瞬时的关节速度 图5 2二自由度机械手手爪沿X0方向运动示意图 解由式 5 2 知 二自由度机械手速度雅可比为 53 因此 逆雅可比为 由式 5 7 可知 且 即vX 1m s vY 0 因此 54 二 机器人动力学方程 机器人动力学的研究有牛顿 欧拉法 拉格朗日法 高斯法 凯恩法及罗伯逊 魏登堡法等 本节介绍动力学研究常用的牛顿 欧拉方程和拉格朗日方程 1 欧拉方程 欧拉方程又称为牛顿 欧拉方程 应用欧拉方程建立机器人机构的动力学方程是指 研究构件质心的运动使用牛顿方程 研究相对于构件质心的转动使用欧拉方程 欧拉方程表征了力 力矩 惯性张量和加速度之间的关系 质量为m 质心在C点的刚体 作用在其质心的力F的大小与质心加速度aC的关系为 式中 F aC为三维矢量 式 2 21 称为牛顿方程 5 8 55 欲使刚体得到角速度为 角加速度为 的转动 则作用在刚体上力矩M的大小为 5 9 式中 M 均为三维矢量 为刚体相对于原点通过质心C并与刚体固结的刚体坐标系的惯性张量 式 5 9 即为欧拉方程 在三维空间运动的任一刚体 其惯性张量可用质量惯性矩IXX IYY IZZ和惯性积IXY IYZ IZX为元素的3 3阶矩阵或4 4阶齐次坐标矩阵来表示 通常将描述惯性张量的参考坐标系固定在刚体上 以方便刚体运动的分析 这种坐标系称为刚体坐标系 简称体坐标系 56 2 拉格朗日方程在机器人的动力学研究中 主要应用拉格朗日方程建立起机器人的动力学方程 这类方程可直接表示为系统控制输入的函数 若采用齐次坐标 递推的拉格朗日方程也可建立比较方便而有效的动力学方程 对于任何机械系统 拉格朗日函数L定义为系统总动能Ek与总势能Ep之差 即 L Ek Ep 5 10 由拉格朗日函数L所描述的系统动力学状态的拉格朗日方程为 5 11 式中 n为连杆数目 qi为系统选定的广义坐标 Fi为作用在第i个坐标上的广义力或力矩 57 3 平面关节机器人动力学分析机器人是一个非线性的复杂动力学系统 动力学问题的求解比