机械系统动力学报告 题目：电梯机械系统的动态特性分析 姓名： 专业： 学号：
电梯机械系统的动态特性分析 一、课题背景介绍 随着社会的快速发展，城市人口密度越来越大，高层建筑不断涌现，因此，现在对电梯的提出了更高的要求，随着科技的进步，在满足客观需求的基础上，电梯向着舒适性，高速，高效的方向发展。在电梯的发展过程中，安全性和功能性一直是电梯公司首要考虑的因素，其中舒适性也要包含在电梯的设计中，避免出现速度或者加速度出现突变，或者电梯运行过程中的振动引起人们的不适。因此，在电梯的设计过程中，对电梯进行动态特性分析是十分必要的。 二、在MATLAB中编程、绘图。 通过同组小伙伴的努力，已经得到了该系统的简化模型与运动方程。因此进行编程： 该系统的微分方程：[][][]{}[]Q x k x c x M= + ? ? ? ? ? ? + ? ? ? ? ? ?? ? ? ，其中矩阵[M]、 [C]、[K]、[Q]都已知。 该系统的微分方程是一个二阶一元微分方程，在MATLAB中，提供有求解常微分方程数值解的函数，其中在MATLAB中常用的求微分方程数值解的有7个：ode45，ode23，ode113，ode15s，ode23s，ode23t，ode23tb 。 ode是MATLAB专门用于解微分方程的功能函数。该求解器有变步长（variable-step）和定步长（fixed-step）两种类型。不同类型有着不同的求解器，其中ode45求解器属于变步长的一种，采用Runge-Kutta
算法；和他采用相同算法的变步长求解器还有ode23。 ode45表示采用四阶，五阶Runge-Kutta单步算法,截断误差为(Δx)^3。解决的是Nonstiff(非刚性)常微分方程。 ode45是解决数值解问题的首选方法，若长时间没结果，应该就是刚性的，可换用ode23试试。 Ode45函数调用形式如下：[T,Y]=ode45(odefun,tspan,y0) 相关参数介绍如下： 通过以上的了解，并对该微分方程进行变换与降阶，得出程序。MATLAB程序： （1）建立M函数文件来定义方程组如下： function dy=func(t,y) dy=zeros(10,1); dy(1)=y(2); dy(2)=1/1660*(-0.006*y(2)+0.003*y(4)-0.0006*y(10)-1.27*10^7*y(1)+1.27*10^7*y (3)+2.54*10^6*y(9)); dy(3)=y(4); dy(4)=1/1600*(+0.03*y(2)-0.007*y(4)+0.003*y(6)+1.27*10^7*y(1)-7.274*10^8*y(3 )+1.27*10^7*y(5)); dy(5)=y(6);
第2章多体系统动力学基本理论
本章主要介绍多体系统动力学的基本理论，包括多刚体系统动力学建模、多柔体系统动力学建模、多体系统动力学方程求解及多体系统动力学中的刚性（Stiff）问题。通过本章的学习可以对多体系统动力学的基本理论有较深入的了解，为具体软件的学习打下良好的理论基础。 2.1 多体系统动力学研究状况 多体系统动力学的核心问题是建模和求解问题，其系统研究开始于20世纪60年代。从60年代到80年代，侧重于多刚体系统的研究，主要是研究多刚体系统的自动建模和数值求解；到了80年代中期，多刚体系统动力学的研究已经取得一系列成果，尤其是建模理论趋于成熟，但更稳定、更有效的数值求解方法仍然是研究的热点；80年代之后，多体系统动力学的研究更偏重于多柔体系统动力学，这个领域也正式被称为计算多体系统动力学，它至今仍然是力学研究中最有活力的分支之一，但已经远远地超过一般力学的涵义。 本节将叙述多体系统动力学发展的历史和目前国内外研究的现状。 2.1.1 多体系统动力学研究的发展 机械系统动力学分析与仿真是随着计算机技术的发展而不断成熟的，多体系统动力学是其理论基础。计算机技术自其诞生以来，渗透到了科学计算和工程应用的几乎每一个领域。数值分析技术与传统力学的结合曾在结构力学领域取得了辉煌的成就，出现了以ANSYS、NASTRAN等为代表的应用极为广泛的结构有限元分析软件。计算机技术在机构的静力学分析、运动学分析、动力学分析以及控制系统分析上的应用，则在二十世纪八十年代形成了计算多体系统动力学，并产生了以ADAMS和DADS为代表的动力学分析软件。两者共同构成计算机辅助工程（CAE）技术的重要内容。 多体系统是指由多个物体通过运动副连接的复杂机械系统。多体系统动力学的根本目的是应用计算机技术进行复杂机械系统的动力学分析与仿真。它是在经典力学基础上产生的新学科分支，在经典刚体系统动力学上的基础上，经历了多刚体系统动力学和计算多体系统动力学两个发展阶段，目前已趋于成熟。 多刚体系统动力学是基于经典力学理论的，多体系统中最简单的情况——自由质点和一般简单的情况——少数多个刚体，是经典力学的研究内容。多刚体系统动力学就是为多个刚体组成的复杂系统的运动学和动力学分析建立适宜于计算机程序求解的数学模型，并寻求高效、稳定的数值求解方法。由经典力学逐步发展形成了多刚体系统动力学，在发展过程中形成了各具特色的多个流派。 早在1687年，牛顿就建立起牛顿方程解决了质点的运动学和动力学问题；刚体的概念最早由欧拉于1775年提出，他采用反作用力的概念隔离刚体以描述铰链等约束，并建立了
“机械动力学”课程教学大纲 英文名称：Mechanical Dynamics 课程编号：MACH3441 学时：32 （理论学时：32 实验学时：课外学时：2实验） 学分：2 适用对象：机械设计、机械制造及自动化、机械电子工程、流体机械、电机、电器、材料工程等本科生高年级。 先修课程：高等数学、普通物理学、理论力学、材料力学、线性代数使用教材及参考书： [1] 石端伟主编. 机械动力学. 北京：中国电力出版社，2007. [2] 张策主编. 机械动力学.北京：高等教育出版社, 2008. [3] 倪振华主编. 振动力学. 西安交通大学出版社，1988. 一、课程性质和目的 性质：专业课 目的： 1.了解机械动力学的研究内容、发展历史以及最新研究进展。 2.培养机械系统动力学分析的基本能力。 3.了解机械系统动力学分析相关的CAE软件。 4.了解机械系统动态测试有关技术。 5.培养查阅和运用相关科技文献进行动力学分析的初步能力。 6.培养创新思维以及解决工程实际问题的能力。 7.培养科学、严谨的工作作风。
二、课程内容简介 随着现代机械装备朝着高精度、高效、大功率的方向发展，其动态性能指标的优劣越来越受到广泛关注和高度重视。机械动力学已日益成为现代机械设计与制造工程领域不可或缺的基础知识。本课程主要介绍机械系统动力分析的基本理论、分析方法、测试与控制技术以及典型机械系统动力学分析方法。通过课程的学习，培养学生能够在机械系统动力分析方面具有明确的基本概念、必要的专业基础知识、一定的机械系统动力分析能力与计算能力。 三、教学基本要求 1．了解相关机械系统动力学分析的新理论、新方法及发展趋向。 2. 掌握有关机械系统动力学分析的基本概念、基本理论与方法。 3. 了解典型机械系统动力学分析流程，具有进行工程实际问题分析的初步能力。 4. 建立正确的机械系统动力分析的思维方式，理论联系实际，具备一定的科研创新精神； 5. 课后需要查阅文献，并开展讨论，完成作业。 四、教学内容及安排 第一章：绪论 1．熟悉研究机械动力学的意义。 2．熟悉机械动力学的主要研究内容。 教学安排及教学方式
多体动力学优化方法 3 李庆国1 ,曾庆良1 ,范文慧 2 (1.山东科技大学机电学院,青岛山东266510;2.清华大学国家C I M S 工程技术研究中心,北京100084) 摘 要: 介绍一种多体动力学优化设计方法,基于I SI GHT 软件集成Pr o /E 和Ada m s,建立优化设 计平台。夹紧装置优化设计实例,验证了该平台的有效性和合理性。关键词: 多体动力学优化;多学科设计优化(MDO );I SI GHT 中图分类号:O313.3 文献标识码:B 文章编号:1001-0874(2007)03-0089-02 A Me thod ofMulti 2body Dynam i c Op ti m i za ti o n L I Q ing 2guo 1 , ZEN G Q ing 2liang 1 , FAN W en 2hui 2 (1.College of Mechanical &Electric Engineering,,Shandong University of Science and Technol ogy,Q ingdao 266510,China; 2.Nati onal C I M S Engineering Research Center of Tsinghua University,Beijing 100084,China ) Ab s trac t: This paper intr oduces a method of multi 2body dyna m ic op ti m izati on design and builds an op ti m izati on design p latfor m based on I SI GHT s oft w are integrati on Pr o /E and Ada m s .The effectiveness and reliability of the p latfor m is validated by taking the op ti m izati on design of chucking fixture as exa mp le .Keywo rd s: multi 2body dyna m ic op ti m izati on;multidisci p linary design op ti m izati on;I SI GHT 3国家自然科学基金资助项目(编号:60474059) 1 多体动力学和MDO 多体系统是多个相互运动的物体通过运动副相联的多刚体系统和多柔体系统。上世纪80年代初,多刚体系统动力学计算机仿真已广泛应用于工程领域,通常用来研究系统的位移、速度、加速度与其受力之间的关系。随着计算机技术的飞速发展,仿真、优化技术已在多体系统设计中得到大量应用。为了解决不同学科间的协同设计问题,人们提出了多学科设计优化的思想。 多学科设计优化(Multidisci p linary Design Op ti 2m izati on 简称MDO )是一种设计复杂系统和子系统的方法论。通过充分利用各个学科(子系统)之间相互作用所产生的协同效应,获得系统的整体最优解[1] 。 多学科设计优化问题,在数学形式上可表达为:寻找:x 最小化:f =f (x,y ) 约束:h i (x,y )=0 (i =1,2,3,…,n ) g i (x,y )≤0 (j =1,2,3,…,m ) 其中f 为目标函数,x 为设计变量,y 是状态变量,h i (x,y )是等式约束,g i (x,y )是不等式约束。 MDO 的研究主要分为三个方面:面向设计的多 学科分析设计软件的集成;有效的MDO 算法,实现多学科并行设计,获得系统最优解;MDO 分布式计算的环境支持。目前,已经出现了较成熟的商业软件,I SI GHT 就是典型代表。2 多学科优化软件I SI GHT I SI GHT 是一个通过软件协同驱动产品设计优 化的软件。特色是融合了优化设计中需要的三大主要功能:自动化功能、集成化功能和最优化功能。 (1)自动化功能 I SI GHT 的过程集成(Pr ocess I ntegrati on )功能可 以对各种CAD 或CAE 软件进行自动化启动、监视和控制,文件解析(File Parser )功能可以自动地编辑、生成输入文件和自动处理输出文件及读取计算结果。 (2)集成化功能 ? 98?2007年第3期 煤 矿 机 电
有限元软件与多体动力学软件 数值分析技术与传统力学的结合在结构力学领域取得了辉煌的成就，出现了以ANSYS 、NASTRAN 等为代表的应用极为广泛的结构有限元分析软件。计算机技术在机构的静力学分析、运动学分析、动力学分析以及控制系统分析上的应用，则在二十世纪八十年代形成了计算多体系统动力学，并产生了以ADAMS 和DADS 为代表的动力学分析软件。两者共同构成计算机辅助工程（CAE ）技术的重要内容。 商业通用软件的广泛应用给我们工程师带来了极大的便利，很多时候我们不需要精通工程问题中的力学原理，依然可以通过商业软件来解决问题，不过理论基础的缺失还是会给我们带来不少的困扰。随着动力有限元与柔性多体系统分析方法的成熟，有时候正确区分两者并不是很容易。 机械领域应用比较广泛的有两类软件，一类是有限元软件，代表的有：ANSYS, NASTRAN, ABAQUS, LS-DYNA, Dytran 等；另一类是多体动力学软件，代表的有ADAMS, Recurdyn ， Simpack 等。在使用时，如何选用这两类软件并不难，但是如果深究这两类软件根本区别并不容易。例如，有限元软件可以分析静力学问题，也可以分析“动力学”问题，这里的“动力学”与多体动力学软件里面的动力学一样吗？有限元软件在分析动力学问题时，可以模拟物体的运动，它与多体动力学软件中模拟物体运动相同吗？多体动力学软件也可以分析柔性体的应力、应变等，这与有限元软件分析等价吗？ 1 有限元软件 有限单元法是一种数学方法，不仅可以计算力学问题，还可以计算声学，热，磁等多种问题，我们这里只探讨有限元法在机械领域的应用。 计算结构应力、应变等的力学基础是弹性力学，弹性力学亦称为弹性理论，主要研究弹性体在外力作用或温度变化等外界因素下所产生的应力、应变和位移，从而为工程结构或构件的强度、刚度设计提供理论依据和计算方法。也就是说用有限元软件分析力学问题时，是用有限元法计算依据弹性力学列出的方程。 考虑下面这个问题，在()0t ， 时间内给一个结构施加一个随时间变化的载荷()P t ，我们希望得到结构的应力分布，在刚刚施加载荷的时候，结构中的应力会有波动，应力场是变化的，但很久以后，应力场趋于稳定。 如果我们想得到载荷施加很久以后，稳定的应力场分布，那么应该用静力学分析方法分析
一曲柄滑块机构运动学仿真 1、设计任务描述 通过分析求解曲柄滑块机构动力学方程，编写matlab程序并建立Simulink 模型，由已知的连杆长度和曲柄输入角速度或角加速度求解滑块位移与时间的关系，滑块速度和时间的关系，连杆转角和时间的关系以及滑块位移和滑块速度与加速度之间的关系，从而实现运动学仿真目的。 2、系统结构简图与矢量模型 下图所示是只有一个自由度的曲柄滑块机构，连杆与长度已知。 图2-1 曲柄滑块机构简图 设每一连杆（包括固定杆件）均由一位移矢量表示，下图给出了该机构各个杆件之间的矢量关系 图2-2 曲柄滑块机构的矢量环
3．匀角速度输入时系统仿真 3.1 系统动力学方程 系统为匀角速度输入的时候，其输入为输出为；。 (1) 曲柄滑块机构闭环位移矢量方程为: (2)曲柄滑块机构的位置方程 (3)曲柄滑块机构的运动学方程 通过对位置方程进行求导，可得 由于系统的输出是与，为了便于建立A*x=B形式的矩阵，使x=[], 将运动学方程两边进行整理，得到 将上述方程的v1与w3提取出来，即可建立运动学方程的矩阵形式 3.2 M函数编写与Simulink仿真模型建立 3.2.1 滑块速度与时间的变化情况以及滑块位移与时间的变化情况 仿真的基本思路：已知输入w2与，由运动学方程求出w3和v1，再通过积分，即可求出与r1。 （1）编写Matlab函数求解运动学方程 将该机构的运动学方程的矩阵形式用M函数compv（u）来表示。 设r2=15mm，r3=55mm，r1（0）=70mm，。 其中各个零时刻的初始值可以在Simulink模型的积分器初始值里设置
M函数如下： function[x]=compv(u) %u(1)=w2 %u(2)=sita2 %u(3)=sita3 r2=15; r3=55; a=[r3*sin(u(3)) 1;-r3*cos(u(3)) 0]; b=[-r2*u(1)*sin(u(2));r2*u(1)*cos(u(2))]; x=inv(a)*b; （2）建立Simulink模型 M函数创建完毕后，根据之前的运动学方程建立Simulink模型，如下图： 图3-1 Simulink模型 同时不要忘记设置r1初始值70，如下图： 图3-2 r1初始值设置
1. 绝对节点坐标法 传统有限元方法建立的单元为非等参数单元，其使用节点处的位移梯度来描述物体的无限小的转动，但在物体发生大变形时，节点处的位移梯度已不能准确描述物体的转动变形，从而极大影响到计算的精度。 Shabana [1]提出了绝对节点坐标法（Absolute nodal coordinate formulation, ANCF ），其理论基础主要是有限元和连续介质力学理论。该方法将物体的单元节点坐标定义在全局坐标系下，使用节点处的斜率(slope)矢量作为节点坐标而不是节点处的无限小转动[2]，不需要另外计算刚体位移与柔性变形之间的耦合，能较精确地计算大变形的多体系统动力学问题。其最终推导出的多体系统的微分代数方程组（DAEs ）中，质量矩阵是一个常数矩阵，但刚度矩阵将是一个非线性的时间函数。 1.1梁单元的绝对节点坐标法 Shabana 首先推导出一维梁单元的绝对节点坐标法模型[1][3]。在这种模型中，梁单元用中性轴来简化，如图1所示，其上面任意一点P 在全局坐标系下的坐标表达为： 23101232320123r =Se r a a x a x a x r b b x b x b x ??+++??==????+++???? 图1 其中，x 为沿轴线的单元局部坐标，[]0,x l ∈，l 为梁单元初始长度；S 为单元形函数；e 为含有8个单元节点坐标的广义坐标矢量。 123456781102205162e []|,|,|,|, T x x x l x l e e e e e e e e e r e r e r e r ========= 1 2 1 2 304078,,,x x x l x l r r r r e e e e x x x x ====????====????
单自由度杆机构的Adams动力学仿真 摘要：文章分析了单自由度的铰链机构的动力学问题，已知原动件曲柄的转矩，绘制输出件摆杆的运动曲线。首先在Adams软件中构造连杆，添加三个连杆，使其成一定角度，相互连接。再在两杆之间添加转动副，并且头尾连杆与地相连。并在曲柄处加转矩，最后进行仿真，并绘出相应图表。 关键词：铰链机构；Adams仿真 1、机构模型的建立 根据题目要求，选择一个铰链四杆机构——曲柄摇杆机构为模型，其结构简图如图1所示。其中，曲柄1为原动件。 图1曲柄摇杆机构简图 在Adams软件中，建立该曲柄摇杆机构的模型如图2所示。 图2 Adams中的曲柄摇杆机构模型
曲柄摇杆机构各连杆的惯性参数参考表1。杆件的材料均选择钢材（密度ρ=7.801×10-6 kg?mm-3，杨氏模量E=2.07×105 N?mm-2，泊松比μ=0.29）。 表1 传动导杆机构各部件惯性参数 2、利用Adams软件添加约束和力矩 杆1和地之间有转动副，杆1和杆2、杆2和杆3之间有转动副，杆3和地之间有转动副。杆1为原动件，在杆1上添加转矩。转矩大小为30。 图3约束与转矩 3、进行仿真 点击仿真按钮，开始仿真，选择仿真时间为2s，可以观察到该机构各个时间的运动状态如图4和图5所示。
（a）T=0时刻（b）T=1时刻 图4仿真过程中机构模型的运动状态 （a）T=1.2时刻（b）T=2时刻 图5仿真过程中机构模型的运动状态 结论 当原动件曲柄的转矩取为30时，点击“后处理”，可以绘制出输出件摆杆的位移曲线、角速度曲线、加速度曲线分别如图10、图11和图12所示。 图10输出件摆杆的位移曲线
硕士学位课程考试试卷 考试科目：电液伺服控制 考生姓名：刘双龙 考生学号：20140713189 学院：机械工程学院专业：机械工程 考生成绩： 任课老师(签名) 考试日期：2014年1月20日午时至时
考试主题：电液伺服（比例）系统 考试题目： 1、为什么把液压控制阀称为液压放大元件？ 2、什么叫阀的工作点？零位工作点的条件是什么? 3、电液伺服阀由哪几部分组成?各部分的作用是什么？ 4、什么是液压固有频率?在阀控缸系统中液压固有频率与活塞位 置有关吗?为什么? 5、为什么电液伺服系统一般都要加校正装置? 6、结合自己研究领域，写一篇液压伺服系统建模、分析的论文， 字数不少于2000字。 注：要求独立完成，不允许抄袭。 交作业时间： 最迟2015年第一个学期的第一周交到7教136，交纸质档。
三自由度平台液压伺服系统建模 摘要： 我的专业是机械工程，主要方向是机械设计，所以本文选择了与我专业方向有关的一个机构进行建模。本文开始对机构进行了说明（采用已有的机构，并非自己设计），然后对其进行运动学分析，从而的到上平台和下平台的速度及加速度，和雅可比矩阵及液压缸速度。然后对驱动机构进行电液伺服系统建模。其中 一：自由度运动平台系统简介 本文所研究的三自由度运动平台类似与六自由度平台是由一个上平台(动平台)、地基(下平台)、三个支杆、三个线性作动器以及若干关节连接而成的。上平台装有负载，完成既定的位置、速度、加速度运动要求，进而实现刑于道路状况的复现。其结构示意图如图1．1所示。 图 1三自由度运动平台的结构图 该平台的结构如下：上平台与地面之间以三个支杆(strut)来约束并起支撑作用，并以三个液压缸作为驱动部件进行驱动。每个液压缸两端为关节轴承，中间为一个移动副和一个转动副连接；每根支杆两端也是采用关节轴承分别与地面和上平台相连中间一个转动副。通过计算可知每个支杆所在的支路都具有5个自由度，每个支路对上平台提供一个约束；每个液压作动器所在的支路都具有6个自由度，对于上平台没有约束。通过每个分支对上平台的约束很容易计算得出其自由度为3。因此，通过三套液压作动器的驱动，上平台能够实现对于给定运动的跟踪复现。 简单直观的对运动进行分析可得到：由于三根支杆的限制作用，上平台平动受到限制：而转动自由度相对更为自由，运动范围更大。当两竖直作动器差动动
第十四章 机械系统动力学 14-11、在图14-19中，行星轮系各轮齿数为123z z z 、、，其质心与轮心重合，又齿轮1、2对质心12O O 、的转动惯量为12J J 、，系杆H 对的转动惯量为H J ，齿轮2的质量为2m ，现以齿轮1为等效构件，求该轮系的等效转动惯量J ν。 2222 2121221 12323121 13212 1 13222 12311212213121313 ( )()()()1()()()( )()()()o H H H o H J J J J m z z z z z z z z z O O z z z z z z z O O J J J J m z z z z z z z z νννωωω ωωωω ωω ωωωωνω=+++=-= += +=+-=++++++解： 14-12、机器主轴的角速度值1()rad ?从降到时2()rad ?，飞轮放出的功 (m)W N ，求飞轮的转动惯量。 max min 122 2 121 ()2 2F F Wy M d J W J ?ν??ωωωω==-=-? 解： 14-15、机器的一个稳定运动循环与主轴两转相对应，以曲柄和连杆所组成的转动副A 的中心为等效力的作用点，等效阻力变化曲线c A F S ν-如图14-22所示。等效驱动力a F ν为常数，等效构件（曲柄）的平均角速度值25/m rad s ?=， 3 H 1 2 3 2 1 H O 1 O 2
不均匀系数0.02δ=，曲柄长度0.5OA l m =，求装在主轴（曲柄轴）上的飞轮的转动惯量。 (a) W v 与时间关系图 （b ）、能量指示图 a 2 24()2 3015m Wy=25N m 25 6.28250.02 c va OA vc OA OA va F W W F l F l l F N Mva N J kg m νν=∏?∏=∏+==∏= =?解：稳定运动循环过程 14-17、图14-24中各轮齿数为12213z z z z =、，，轮1为主动轮，在轮1上加力矩1M =常数。作用在轮 2 上的阻力距地变化为： 2r 22r 020M M M ??≤≤∏==∏≤≤∏=当时，常数；当时，，两轮对各自中心的转动惯量为12J J 、。轮的平均角速度值为m ω。若不均匀系数为δ，则：（1）画出以轮1为等效构件的等效力矩曲线M ν?-；（2）求出最大盈亏功；（3）求飞轮的转动惯量F J 。 图14-24 习题14-17图 40Nm 15∏ 12.5∏ 22.5∏ 15Nm ∏ 2∏ 2.5∏ 4∏ 25∏ 1 1 z 2 z 2 r M 2 M ∏ 2∏ 2?
太原理工大学研究生试题 姓名： 学号： 专业班级： 机械工程2014级 课程名称: 《机械系统动力学》 考试时间: 120分钟 考试日期: 题号 一 二 三 四 五 六 七 八 总分 分数 1 圆柱型仪表悬浮在液体中，如图1所示。仪表质量为m ，液体的比重为ρ，液体的粘性阻尼系数为r ，试导出仪表在液体中竖直方向自由振动方程式，并求固有频率。（10分） 2 系统如图2所示，试计算系统微幅摆动的固有频率，假定OA 是均质刚性杆，质量为m 。（10分） 3 图3所示的悬臂梁，单位长度质量为ρ，试用雷利法计算横向振动的周期。假定梁的 变形曲线为?? ? ?? -=x L y y M 2cos 1π（y M 为自由端的挠度）。（10分） 4 如图4所示的系统，试推导质量m 微幅振动的方程式并求解θ(t)。（10分） 5 一简支梁如图5所示，在跨中央有重量W 为4900N 电机，在W 的作用下，梁的静挠度δst=，粘性阻尼使自由振动10周后振幅减小为初始值的一半，电机n=600rpm 时，转子不平衡质量产生的离心惯性力Q=1960N ，梁的分布质量略去不计，试求系统稳态受迫振动的振幅。（15分） 6 如图6所示的扭转摆，弹簧杆的刚度系数为K ，圆盘的转动惯量为J ，试求系统的固有频率。(15分) 7如图7一提升机，通过刚度系数m N K /1057823?=的钢丝绳和天轮（定滑轮）提升货载。货载重量N W 147000=，以s m v /025.0=的速度等速下降。求提升机突然制动时的钢丝绳最大张力。(15分) 8某振动系统如图8所示，试用拉个朗日法写出动能、势能和能量散失函数。(15分) 太原理工大学研究生试题纸
Assignment Vehicle system dynamics simulation 学院：机电学院 专业：机械工程及自动化 姓名： 指导教师：
The model we are going to analys: The FBD of the suspension system is shown as follow:
According to the New's second Law, we can get the equation: 2 )()(221211mg z z c z z k z m --+-=???? 221212)()(z k mg z z c z z k z m w +-----=? ??? 0)()()()(222111222111=-++--+-++--+? ? ? ? ? ? ? ?w w w w z L z k z L z k z L z c z L z c z m χχχχ 0)()()()(2222111122221111=-++----++---? ? ? ? ? ? ? ?w w w w z L z L k z L z L k z L z L c z L z L c J χχχχχ d w w w w Q z L z k z L z c z m ,111111111)()(-=------? ? ? ? ?χχ d w w w w Q z L z k z L z c z m ,222222222)()(-=-+--+-? ????χχ When there is no excitation we can get the equation: 2)()(221211mg z z c z z k z m --+-=???? 2 21212)()(z k mg z z c z z k z m w +-----=? ??? Then we substitude the data into the equation, we write a procedure to simulate the system: Date: ???? ?? ??? ??==?==?===MN/m 0.10k m 25.1s/m kN 0.20MN/m 0.1m kg 3020kg 2100kg 3250w 2l c k I m m by w b
图4 机器安装示意图 88、一个质量20Kg 的机器，按图4所示方式安装。若弹簧的总刚度 为17KN/m ，总阻尼为300m s N ?。试求初始条mm x 250=，s mm x 3000= 时的振动响应。 88、解：由0=++kx x c x m 代入数据后得 08501501017300203=++=?++x x x x x x （8分） 其中，152=a ，8502=n ω，计算阻尼比和固有圆频率 17.2826.012.291126.02 .295.722=-?=-=<===ζωωωζn d n a （4分） 将初始条件代入 00020020arctan )(ax x x ax x x A d d +=++= ω?ω （4分） 得： o d d ax x x mm ax x x A 3.555.25.730017.2825arctan arctan )(4.30)17.2825.7300(25)(0002220020?+?=+==?++=++= ω?ω（2分）
则系统的振动响应为 4. 305.7+ =-t x t（2分）e sin( 28 ) 96 .0 . 17
1. “机械动力学”主要研究哪些内容，请以任一机器为对象举例说明研究内容及其相互关系。 答：机械动力学是研究机械在力的作用下的运动和机械在运动中产生的力，并从力与运动的相互作用的角度进行机械设计和改进的科学。动力学主要研究内容概括起来有：1，共振分析；2，振动分析与动载荷计算；3，计算机与现代测试技术的运用；4，减震与隔振。柴油机上的发动机，发动机不平衡时会产生很强的地面波，从而产生噪声，而承受震动的结构，发动机底座，会由于振动引起的交变应力而导致材料的疲劳失效，而且振动会加剧机械零部件的磨损，如轴承和齿轮的磨损等，并使机械中的紧固件如螺母等变松。在加工时还会导致零件加工质量变差。通过对共振的研究和分析，使机械的运转频率避免共振区，避免机械共振事故的发生，通过振动分析与动载荷计算可以求出在外力作用下机械的真实运动，运用计算机和现代测试技术对机械的运行状态进行检测，以及故障诊断，模态分析以及动态分析，现实中机器运转时由于各种激励因素的存在，不可避免发生振动，为了减小振动，通常在机器底部加装弹簧，橡胶等隔振材料。 2.简述在刚性运动前提下，如何进行运动构件的真实运动分析求解（请列出步骤）？ 答：首先建立等效力学模型，将复杂的机械系统简化为一个构件，即等效构件，根据质点系动能定理，将作用于机械系统上的所有外力和外力矩、所有构件的质量和转动惯量，都向等效构件转化；其次计算等效构件上的等效量（包括等效力矩，等效力，等效质量，等效转动惯量）；再次建立等效构件的运动方程式，有两种形式，能量形式和力矩形式；最后通过方程式求出等效构件的角速度函数和角加速度函数，这样便可以求出机械系统的真实运动规律。 3.在弹性运动假设下，有哪些弹性动力学建模方法，各有什么特点？请解释“瞬时刚化” 的概念。） 答：弹性动力学模型有集中参数模型和有限元模型。集中参数模型建立起的运动方程为常微分方程，但是由于质量简化过多，模型粗糙，精度比较差；有限元建立的运动方程也为常微分方程，但相较集中参数模型精确，适应性广，可以模拟复杂形状的构件，运算模型统一。瞬时刚化：机构在运动到循环中的某一位置时，可将机构的形状和作用在其上的载荷瞬时冻结起来，从而可瞬时的将机构看做一个刚体结构。
机械系统动力学试题 一、 简答题： 1．机械振动系统的固有频率与哪些因素有关？关系如何？ 2．简述机械振动系统的实际阻尼、临界阻尼、阻尼比的联系与区别。 3．简述无阻尼单自由度系统共振的能量集聚过程。 4. 简述线性多自由度系统动力响应分析方法。 5. 如何设计参数，使减振器效果最佳？ 二、 计算题： 1、 单自由度系统质量Kg m 10=, m s N c /20?=, m N k /4000=, m x 01.00=, 00=? x ,根据下列条件求系统的总响应。 （a ） 作用在系统的外激励为t F t F ωcos )(0=,其中N F 1000=， s rad /10=ω。 （b ） 0)(=t F 时的自由振动。 2、 质量为m 的发电转子，它的转动惯量J 0的确定采用试验方法：在转子径向R 1的地方附加一小质量m 1。试验装置如图2所示，记录其振动周期。 a ）求发电机转子J 0。 b ）并证明R 的微小变化在R 1=（m/m 1+1）·R 时有最小影响。 3、 如图3所示扭转振动系统，忽略阻尼的影响 J J J J ===321，K K K ==21 （1）写出其刚度矩阵； （2）写出系统自由振动运动微分方程； （2）求出系统的固有频率； （3）在图示运动平面上，绘出与固有频率对应的振型图。 1 θ（图2）
（图3） 4、求汽车俯仰振动（角运动）和跳振（上下垂直振动）的频率以及振 动中心（节点）的位置（如图4）。参数如下：质量m=1000kg，回转半径r=0.9m，前轴距重心的距离l1=0.1m,后轴距重心的距离l2=1.5m,前弹簧刚度k1=18kN/m,后弹簧刚度k2=22kN/m （图4） 5、如5图所示锻锤作用在工件上的冲击力可以近似为矩形脉冲。已知 工件，铁锤与框架的质量为m1=200 Mg，基础质量为m2=250Mg，弹簧垫的刚度为k1=150MN/m，土壤的刚度为k2=75MN/m.假定各质量的初始位移与速度均为零，求系统的振动规律。
多体系统动力学简介
多体系统动力学研究对象——机构 工程中的对象是由大量零部件构成的系统。在对它们进行设计优化与性态分析时可以分成两大类 一类为结构 ——正常工况下构件间没有相对运动（房屋建筑，桥梁等） ——关心的是这些结构在受到载荷时的强度、刚度与稳定 一类为机构 ——系统在运动过程中这些部件间存在相对运动（汽车，飞机起落架。机器人等）——力学模型为多个物体通过运动副连接的系统，称为多体系统 多体系统动力学俄研究的对象——机构（复杂机械系统）
不考虑系统运动起因的情况下研究各部件的位置与姿态及其变化速度和加速度的关系 典型案例：平面和空间机构的运动分析 系统各部件间通过运动副与驱动装置连接在一起 数学模型：各部件的位置与姿态坐标的非线性代数方程，以及速度与加速度的线性代数方程
当系统受到静载荷时，确定在运动副制约下的系统平衡位置以及运动副静反力 典型案例：机车或汽车中安装有大量的弹簧阻尼器，整车设计中必须考虑系统在静止状态下车身的位置与姿态，为平稳性与操纵稳定性的研究打下基础 数学模型：非线性微分代数方程组
讨论载荷和系统运动的关系 研究复杂机械系统在载荷作用下各部件的动力学响应是工程设计中的重要问题 动力学正问题——已知外力求系统运动的问题 动力学逆问题——已知系统运动确定运动副的动反力，是系统各部件强度分析的基础 动力学正逆混合问题——系统的某部分构件受控，当它们按照某已知规律运动时，讨论在外载荷作用下系统其他构件如何运动 数学模型：非线性微分代数方程组
机械系统的多体系统力学模型 在对复杂机械系统进行运动学与动力学分析前需要建立它的多体系统力学模型。对系统如下四要素进行定义： ?物体 ?铰链 ?外力（偶） ?力元 实际工程中的机械系统多体系统力学模型的定义取决于研究的目的 模型定义的要点是以能揭示系统运动学与动力学性态的最简模型为优 性态分析的求解规模与力学模型的物体与铰的个数有关
车辆动力学主要仿真软件 I960年，美国通用汽车公司研制了动力学软件DYNA主要解决多自由度无约束的机械系统的动力学问题，进行车辆的“质量－弹簧－阻尼”模型分析。作为第一代计算机辅助设计系统的代表，对于解决具有约束的机械系统的动力学问题，工作量依然巨大，而且没有提供求解静力学和运动学问题的简便形式。 随着多体动力学的诞生和发展，机械系统运动学和动力学软件同时得到了迅速的发展。1973 年,美国密西根大学的N.Orlandeo 和，研制的ADAM 软件，能够简单分析二维和三维、开环或闭环机构的运动学、动力学问题，侧重于解决复杂系统的动力学问题，并应用GEAR 刚性积分算法，采用稀疏矩阵技术提高计算效率。1977年,美国Iowa 大学在，研究了广义坐标分类、奇异值分解等算法并编制了DADS软件，能够顺利解决柔性体、反馈元件的空间机构运动学和动力学问题。随后，人们在机械系统动力学、运动学的分析软件中加入了一些功能模块，使其可以包含柔性体、控制器等特殊元件的机械系统。 德国航天局DLF早在20世纪70年代，Willi Kort tm教授领导的团队就开始从事MBS软件的开发，先后使用的MBS软件有Fadyna (1977)、MEDYNA1984),以及最终享誉业界的SIMPAC(1990). 随着计算机硬件和数值积分技术的迅速发展，以及欧洲航空航天事业需求的增长，DLR决定停止开发基于频域求解技术的MED YN软件，并致力于基于时域数值积分技术的发展。1985年由DLR开发的相对坐标系递归算法的SIMPACI软件问世，并很快应用到欧洲航空航天工业，掀起了多体动力学领域的一次算法革命。 同时，DLR首次在SIMPAC嗽件中将多刚体动力学和有限元分析技术结合起来，开创了多体系统动力学由多刚体向刚柔混合系统的发展。另外，由于SIMPACI算法技术的优势，成功地将控制系统和多体 计算技术结合起来，发展了实时仿真技术。
IRB 7600大功率机器人运动仿真
摘要 (2) 1引言 (3) 2机器人发展概述 (3) 2.1机器人的三大定律产生 (3) 2.2工业机器人的发展和特点 (3) 2.3工业机器人现状与前景 (5) 3 ABB机器人和大功率机器人的发展概述 (5) 3.1 ABB公司的发展 (6) 3.2 ABB工业机器人的现状 (6) 3.3简述IRB 7600机器人特点 (6) 3.４IRB 7600机器人的主要参数和应用 (7) 4. 基于ADAMS的IRB 7600大功率机器人运动学仿真 (8) 4.1 IRB 7600大功率机器人的运动学分析 (8) 4.1.1分析IRB 7600大功率机器人得到简图，建立方程 (9) 4.1.2 IRB 7600大功率机器人正向运动学解 (11) 4.2ADAMS中的的运动仿真 (12) 4.2.1在ADAMS中建立IRB 7600机器人的模型 (12) 4.2.2运动的施加 (14) 4.2.4运动结果分析 (16) 总结 (19) 参考文献 (20)
摘要 现代机器人技术飞速发展，其中工业机器人的应用也越来越广泛，成为高科技中极为重要的组成部分。本文主要针对ABB机器中的IRB 7600大功率机器人，对其运动进行仿真探究，学习机器人的一般运动方法。 ABB大功率机器人系列开辟了全新的应用领域，该机器人有多种版本，最大承重能力高达650kg。IRB 7600适合用于各行业重载场合，大转矩、大惯性、刚性结构以及卓越的加速性能等优良特性使这款市场主导产品声誉日隆。用于装配、清洁/喷涂、切割/去毛刺、研磨/抛光、机械管理、物料搬运、货盘堆跺、扳弯机管理、点焊，应用前景广。通过对IRB 7600的模型建立，基于ADAMS的点焊机器人运动学仿真，得到了机器人的仿真运动曲线和模型图。对模型和曲线分析，初步的了解到大功率机器人的运动和工作方式。 关键字：IRB 7600、ABB、ADAMS、仿真
《机械动力学》期末复习题及答案1、判断 1.机构平衡问题在本质上是一种以动态静力分析为基础的动力学综合，或动力学设计。 答案:正确 2.优化平衡就是采用优化的方法获得一个绝对最佳解。 答案:错误 3.惯性力的计算是建立在主动构件作理想运动的假定的基础上的。 答案:正确 4.等效质量和等效转动惯量与机械驱动构件的真实速度无关。 答案:正确 5.作用于等效构件上的等效力（或等效力矩）所作的功等于作用于系统上的外力所作的功。答案: 错误 6.两点动代换后的系统与原有系统在静力学上是完全等效的。 答案:错误 7.对于不存在多余约束和多个自由度的机构，动态静力分析是一个静定问题。 答案:错误 8.摆动力的完全平衡常常会导致机械结构的简单化。 答案:错误 9.机构摆动力完全平衡的条件是：机构运动时，其总质心作变速直线运动。 答案:错误 10.等效质量和等效转动惯量与质量有关。 答案:错误 11.平衡是在运动设计完成之前的一种动力学设计。 答案:错误 12.在动力分析中主要涉及的力是驱动力和生产阻力。 答案:正确 13.当取直线运动的构件作为等效构件时，作用于系统上的全部外力折算到该构件上得到等效力。答案:正确 14.摆动力的平衡一定会导致机械结构的复杂化。 答案:错误 15.机器人操作机是一个多自由度的闭环的空间机构。 答案:错误 16.质量代换是将构件的质量用若干集中质量来代换，使这些代换质量与原有质量在运动学上等效答案:正确 17.弹性动力分析考虑构件的弹性变形。 答案:正确 18.机构摆动力矩完全平衡的条件为机构的质量矩为常数。 答案:错误
19.拉格朗日方程是研究约束系统静力动力学问题的一个普遍的方法。 答案:正确 20.在不含有变速比传动而仅含定速比传动的系统中，传动比为常数。 答案:正确 21.平衡分析着眼于全部消除或部分消除引起震动的激振力。 答案:正确 22.通路定理是用来判断能否实现摆动力完全平衡的理论。 答案:错误 23.无论如何，等效力与机械驱动构件的真实速度无关。 答案:正确 24.综合平衡不仅考虑机构在机座上的平衡，同时也考虑运动副动压力的平衡和输入转矩的平衡。答案:正确 25.速度越快，系统的固有频率越大。 答案:错误 26.平衡的实质就是采用构件质量再分配等手段完全地或部分地消除惯性载荷。 答案:正确 27.优化综合平衡是一个多目标的优化问题，是一种部分平衡。 答案:正确 28.机构摆动力完全平衡的条件为机构的质量矩为常数。 答案:正确 29.当以电动机为原动机时，驱动力矩是速度的函数。 答案:错误 30.为了使得等效构件的运动与机构中该构件的运动一致，要将全部外力等效地折算到该机构上这 一折算是依据功能原理进行的。 答案:正确 2、单选 1.动力学反问题是已知机构的（），求解输入转矩和各运动副反力及其变化规律。 A.运动状态 B.运动状态和工作阻力 C.工作阻力 D.运动状态或工作阻力 答案:B 2.平衡的实质就是采用构件质量再分配等手段完全地或部分地消除（）。 A.加速度 B.角加速度 C.惯性载荷 D.重力 答案: C 3.摆动力的完全平衡常常会导致机械结构的（）。 A.简单化