《机械设计仿真软件.ppt》由会员分享,可在线阅读,更多相关《机械设计仿真软件.ppt(180页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、机械系统设计建模与仿真主讲人:张玉华 机械工程学院 车辆工程系第一页,编辑于星期日:十四点 五十四分。教学安排教学安排计划学时计划学时 48 其中其中 20学时实验学时实验教材:教材:机械系统设计建模与仿真(兼上机实验指导书)机械系统设计建模与仿真(兼上机实验指导书)张玉华张玉华 主编主编 上课时间:上课时间:119周周 20周考试周考试实验安排:实验实验安排:实验1 ADAMS基本操作基本操作 第第14周周 实验实验2 几何建模与参数化几何建模与参数化 第第15周周 实验实验3 机构约束与施加载荷机构约束与施加载荷 第第16周周 实验实验4 编辑样机模型编辑样机模型 第第17周周 实验实验5
2、 样机仿真分析样机仿真分析 第第18周周 地点:地点:机械楼机械楼 2层层 CAD中心中心 胡老师指导胡老师指导第二页,编辑于星期日:十四点 五十四分。上课要点上课要点1.上课精力集中,上课精力集中,认真思考认真思考2.认真做好笔记,按时完成作业认真做好笔记,按时完成作业3.遵守课堂纪律遵守课堂纪律 4.(不迟到,不早退,不开手机)(不迟到,不早退,不开手机)第三页,编辑于星期日:十四点 五十四分。第一章 绪论11 机械系统的设计机械系统的设计12 多刚体系统动力学多刚体系统动力学13 牛顿牛顿-欧拉方法欧拉方法14 虚拟样机技术虚拟样机技术第四页,编辑于星期日:十四点 五十四分。1 11 1
3、 机械系统的设计机械系统的设计1.机器机器2.传动传动3.机构机构4.零件零件单缸内燃机单缸内燃机第五页,编辑于星期日:十四点 五十四分。牛头刨床1.机器机器2.传动传动3.机构机构4.零件零件第六页,编辑于星期日:十四点 五十四分。机器、机构、机械系统机器、机构、机械系统机构:机构:是由两个以上具有相对运动的构件组成是由两个以上具有相对运动的构件组成的系统,机构的作用在于传递运动或改变运的系统,机构的作用在于传递运动或改变运动的形式。动的形式。机机器器:是是由由若若干干机机构构组组成成的的系系统统。例例如如,内内燃燃机机包包含含曲曲柄柄滑滑块块机机构构、齿齿轮轮机机构构和和控控制制进气与排气
4、的凸轮机构。进气与排气的凸轮机构。机械系统:是机构与机器的总称。它由许多构件机械系统:是机构与机器的总称。它由许多构件和零件组成。和零件组成。第七页,编辑于星期日:十四点 五十四分。构件与构件与零件的区别零件的区别构件是运动的单元;构件是运动的单元;零件是制造的单元。零件是制造的单元。构件构件:组成机构的各个组成机构的各个相对运动部分相对运动部分称为称为构件构件。构件构件 可以是单一的整体,也可以是几个元件的可以是单一的整体,也可以是几个元件的刚性组合刚性组合。零件零件:组成构件的元件则称为:组成构件的元件则称为零件零件。第八页,编辑于星期日:十四点 五十四分。机械系统设计的基本问题 机械系统
5、设计的基本问题是机构机械系统设计的基本问题是机构的的综合综合、运动学运动学和和动力学动力学分析与设计。分析与设计。机构机构综综合合着重研究着重研究创创造性构思、造性构思、发发明、明、创创新新设计设计新机构的理新机构的理论论和方法。和方法。而而机构的运机构的运动动学学和和动动力学分析力学分析,一方面是用于,一方面是用于现现有机械系有机械系统统的性的性能分析与改能分析与改进进,另一方面是,另一方面是为为机构的机构的综综合提供理合提供理论论依据。因而它依据。因而它们们是机械系是机械系统设计统设计中重点研究的内容,也是本中重点研究的内容,也是本书书要重点介要重点介绍绍的内的内容。容。第九页,编辑于星期
6、日:十四点 五十四分。本课程的任务本课程的任务熟悉熟悉多刚体系统动力学多刚体系统动力学的基本概念、基本理论,掌握建的基本概念、基本理论,掌握建立机构的运动分析和动力分析数学模型的方法。立机构的运动分析和动力分析数学模型的方法。以多刚体系统动力学为理论指导,以多刚体系统动力学为理论指导,虚拟样机技术虚拟样机技术为设计为设计手段,研究机械系统的建模方法。手段,研究机械系统的建模方法。熟悉熟悉ADAMS软件的基本操作,掌握机械系统虚拟样软件的基本操作,掌握机械系统虚拟样机的建模和仿真分析方法,提高机械系统的设计质量,机的建模和仿真分析方法,提高机械系统的设计质量,提高机械产品的性能,提高自主知识产权
7、产品的核心提高机械产品的性能,提高自主知识产权产品的核心竞争力。竞争力。第十页,编辑于星期日:十四点 五十四分。12 多刚体系统动力学多刚体系统动力学 关于关于刚刚体的假体的假设设是不考是不考虑虑物体的物体的变变形。但是物体形。但是物体总总是有是有变变形的,形的,而物体的而物体的变变形形对对系系统统的运的运动动也是有影响的,有也是有影响的,有时则时则有决定性的影响,有决定性的影响,因此,因此,严严格地格地讲讲多多刚刚体系体系统应为统应为多体系多体系统统即柔性体系即柔性体系统统。目前,。目前,国内外已从多国内外已从多刚刚体系体系统统的研究的研究扩扩展到展到多体系多体系统统(包括柔体系包括柔体系统
8、统)的研究的研究。但是,在某些情况,比如构件的但是,在某些情况,比如构件的变变形很小,且构件的形很小,且构件的变变形形对对系系统统的的动动力学特性影响不大,仍然可以将力学特性影响不大,仍然可以将这类这类系系统视为统视为多多刚刚体系体系统统。我。我们仅们仅研研究多究多刚刚体系体系统统并以此作并以此作为为研究多体系研究多体系统动统动力学的基力学的基础础。工程中的机械系统大多由许多构件组成,研究这些复杂系统时,工程中的机械系统大多由许多构件组成,研究这些复杂系统时,往往可以将构成系统的各构件简化为往往可以将构成系统的各构件简化为刚体刚体,而刚体之间靠,而刚体之间靠运动副运动副连接,连接,从而得到从而
9、得到“多刚体系统多刚体系统”。例如自行车、曲柄滑块机构、汽车中的转向机构、飞机的起例如自行车、曲柄滑块机构、汽车中的转向机构、飞机的起落架、工业机器人等落架、工业机器人等说说明明:第十一页,编辑于星期日:十四点 五十四分。运动副运动副 连连接构件的运接构件的运动动副,可以是副,可以是圆圆柱柱铰链铰链(两两刚刚体之体之间间有一个相有一个相对转动对转动的自由度的自由度),万向万向联轴节联轴节(两个相两个相对转动对转动自由度自由度),球球铰铰(三个相三个相对转动对转动的自由度的自由度),也可以是其它形式的也可以是其它形式的运运动动学学约约束束(如棱柱形如棱柱形约约束允束允许许一个相一个相对对滑滑动动
10、的自由度的自由度),甚至没有物理意,甚至没有物理意义义上的运上的运动动学学约约束,而只有力的束,而只有力的作用作用(如如弹弹簧簧连连接接),即所,即所谓谓的的广广义铰义铰。第十二页,编辑于星期日:十四点 五十四分。多刚体系统类型 多多刚刚体系体系统统从从结结构上可以分构上可以分为为两两类类:树树状状结结构构和和非非树树状状结结构构。两两类结类结构的区分取决于构的区分取决于“通路通路”的概念。的概念。如果系如果系统统中任意两中任意两刚刚体之体之间间都只有一个通路存在,都只有一个通路存在,则则称系称系统为统为树树状状结结构构,图图中的中的(a)(a)、(c)(c)。如果系。如果系统统中至少有两个中
11、至少有两个刚刚体之体之间间存在两存在两个个(或更多的或更多的)通路,通路,则则称系称系统为统为非非树树状状结结构构,图图中的中的(b)(b),这时这时,从,从BiBi到到BjBj的两个通路构成一个的两个通路构成一个闭闭合合链链。第十三页,编辑于星期日:十四点 五十四分。多刚体系统结构示例 机械系机械系统统中,机械手,空中,机械手,空间飞间飞行器以及人体步行行器以及人体步行时时的的摆动摆动相相都可以都可以视为树视为树状状结结构系构系统统。自行自行车车、曲柄滑、曲柄滑块块机构以及人体站立机构以及人体站立时时的支撑相的支撑相则则可可视为视为非非树树状状结结构系构系统统。第十四页,编辑于星期日:十四点
12、 五十四分。树状结构的分类树状结构的分类 树状结构树状结构是研究多体系统动力学的基础,因为任何非树状结构是研究多体系统动力学的基础,因为任何非树状结构均可将其闭合链打开加上某些附加约束而视为树状结构。树状结构均可将其闭合链打开加上某些附加约束而视为树状结构。树状结构又可以分为两类:又可以分为两类:系系统统中中某某刚刚体体(编编号号为为B1)B1)与与一一运运动动已已知知的的刚刚体体(通通常常称称之之为为基基座座,编号为编号为B B0 0)相铰接,此类称为相铰接,此类称为有根树有根树。典型的如工业机械手。典型的如工业机械手。系系统统中任一中任一刚刚体都不与基座相体都不与基座相连连此此类类称称为为
13、悬悬空空树树。如。如卫卫星、星、腾腾空的运空的运动员动员等。等。第十五页,编辑于星期日:十四点 五十四分。多刚体动力学的特点多刚体动力学的特点 多多刚刚体体动动力学的研究内容同力学的研究内容同样样也分也分为为运运动动学和学和动动力学两部力学两部分,与分,与经经典力学的区典力学的区别别之之处处在于多在于多刚刚体系体系统统是十分复是十分复杂杂的系的系统统,其,其自由度数大自由度数大,且各,且各构构件的运件的运动动一般都有一般都有大位移大位移变变化化,因此,因此,不但运不但运动动微分方程数多微分方程数多,且有大量的,且有大量的非非线线性性项项,一般很,一般很难难求得求得解解析解析解,而必,而必须须借
14、助借助计计算机作数算机作数值计值计算。算。第十六页,编辑于星期日:十四点 五十四分。多刚体动力学的主要研究多刚体动力学的主要研究 寻求建立多刚体系统运动微分方程的解析方法。这种方法应是寻求建立多刚体系统运动微分方程的解析方法。这种方法应是一种规格化的方法,能方便、快捷地统一处理各类问题、面向计算机一种规格化的方法,能方便、快捷地统一处理各类问题、面向计算机的分析方法。的分析方法。发发展展与与各各种种分分析析方方法法配配套套的的算算法法,以以实实现现复复杂杂非非线线性性常微分方程常微分方程(ODE)(ODE)或微分或微分代数方程代数方程(DAE)(DAE)的数值积分。的数值积分。根根据据计计算算
15、结结果果提提供供易易于于分分析析的的各各种种输输出出形形式式,如如曲曲线线、图象、动画等。图象、动画等。应应用用以以上上方方法法对对具具体体系系统统进进行行分分析析,并并解解决决力力学学性性能能分分析析、参数优化、寻求最优控制规律等力学问题。参数优化、寻求最优控制规律等力学问题。第十七页,编辑于星期日:十四点 五十四分。13 牛顿-欧拉方法 牛顿牛顿欧拉法是一种规格化的方法,能方便、欧拉法是一种规格化的方法,能方便、快捷地统一处理各类问题、面向计算机的分析方法。快捷地统一处理各类问题、面向计算机的分析方法。虽然方程数较多,但建立方程的过程却十分简单,虽然方程数较多,但建立方程的过程却十分简单,
16、而且易于编程上机计算。而且易于编程上机计算。下面讨论如何采用牛顿下面讨论如何采用牛顿欧拉方法对曲柄滑块机欧拉方法对曲柄滑块机构进行动力学建模方法。构进行动力学建模方法。第十八页,编辑于星期日:十四点 五十四分。曲柄滑块机构动力学建模曲柄滑块机构动力学建模将将曲曲柄柄滑滑块块机机构构看看作作由由B B1 1和和B B2 2组组成成的的系系统统,解解除除约约束束,如如图图所所示示,X X1 1,Y Y1 1,-X-X1 1,-Y-Y1 1与与Y Y2 2均均为为约束反力。约束反力。列出列出B B1 1、B B2 2的运动微分方程。的运动微分方程。B B1 1只有转动只有转动 B B2 2既有移动又
17、有转动既有移动又有转动 (1.1)(1.1)(1.2)(1.2)第十九页,编辑于星期日:十四点 五十四分。曲柄滑块机构动力学建模曲柄滑块机构动力学建模 从从(1.2)(1.2)式至式至(1.3)(1.3)式,前四个式,前四个为为微分方程,后三个微分方程,后三个为为代数方程,代数方程,共七个方程构成一封共七个方程构成一封闭闭方程,可求得七个未知量方程,可求得七个未知量 约约束条件是束条件是A A1 1与与A A2 2点重合以及点重合以及D D点在点在x x轴轴上,上,由此得到由此得到约约束方程束方程为为 (1.3)所以,采用这种方法将使微分代数方程组中的方程数目增多,所以,采用这种方法将使微分代
18、数方程组中的方程数目增多,但每个方程的建立则要简单得多。但每个方程的建立则要简单得多。第二十页,编辑于星期日:十四点 五十四分。14 虚拟样机技术虚拟样机技术 虚拟样机技术又称为虚拟样机技术又称为机械系统动态仿真技术机械系统动态仿真技术,是国际上,是国际上2020世纪世纪8080年代随着计算机技术的发展而迅速发展起来的年代随着计算机技术的发展而迅速发展起来的一项计算机辅助工程(一项计算机辅助工程(CAECAE)技术。工程师在计算机上)技术。工程师在计算机上建建立样机模型立样机模型,对模型进行各种,对模型进行各种动态性能分析动态性能分析,然后改进,然后改进样机样机设计方案设计方案,用数字化形式代
19、替传统的实物,用数字化形式代替传统的实物样机实验。样机实验。运用虚拟样机技术,可以运用虚拟样机技术,可以大大简化大大简化机械产品的设计开发过机械产品的设计开发过程,程,大幅度缩短大幅度缩短产品开发周期,产品开发周期,大量减少大量减少产品开发费用和成产品开发费用和成本,明显本,明显提高提高产品质量,产品质量,提高提高产品的系统级性能,获得最优产品的系统级性能,获得最优化和创新的设计产品。化和创新的设计产品。第二十一页,编辑于星期日:十四点 五十四分。机械系统动力学自动分析软件机械系统动力学自动分析软件ADAMSADAMS(Automatic Dynamic Analysis of Mechani
20、cal Systems)(Automatic Dynamic Analysis of Mechanical Systems)是美国是美国MDIMDI公司公司(Mechanical Dynamics lnc.)(Mechanical Dynamics lnc.)开发的著名的虚拟样机分开发的著名的虚拟样机分析软件。析软件。ADAMS ADAMS一方面是虚拟样机分析的应用软件,用户可以运用该软一方面是虚拟样机分析的应用软件,用户可以运用该软件非常方便地对虚拟机械系统进行件非常方便地对虚拟机械系统进行静力学静力学、运动学运动学和和动力学动力学分分析析;另一方面,它又是虚拟样机分析另一方面,它又是虚拟样
21、机分析开发工具开发工具,其开放性的程,其开放性的程序结构和多种接口,可以成为特殊行业用户进行特殊类型虚序结构和多种接口,可以成为特殊行业用户进行特殊类型虚拟样机分析的拟样机分析的二次开发工具平台二次开发工具平台。虚拟样机虚拟样机仿真分析基本步骤仿真分析基本步骤如图如图1-31-3所示。所示。第二十二页,编辑于星期日:十四点 五十四分。仿真分析基本步骤仿真分析基本步骤虚拟样机仿真分析虚拟样机仿真分析机械系统机械系统建模建模l几何建模几何建模l施加约束施加约束运动副和运动参数运动副和运动参数l施加载荷施加载荷仿真仿真分析分析l设置测量和仿真输出设置测量和仿真输出l进行仿真分析进行仿真分析仿真结果仿
22、真结果分析分析l回放仿真结果回放仿真结果l绘制仿真结果曲线绘制仿真结果曲线第二十三页,编辑于星期日:十四点 五十四分。仿真分析基本步骤仿真分析基本步骤否否是是改进机械改进机械系统模型系统模型l增加摩擦力,改进载荷函数增加摩擦力,改进载荷函数l定义柔性物体和连接定义柔性物体和连接l定义控制定义控制验证仿真验证仿真分析结果分析结果分析分析l输入实验数据输入实验数据l添加实验数据曲线添加实验数据曲线与实验结果一致与实验结果一致重复仿真重复仿真分析分析l设置可变参数点设置可变参数点l定义设计变量定义设计变量机械系统机械系统优化分析优化分析l进行主要设计影响因素研究进行主要设计影响因素研究l进行试验设计
23、研究进行试验设计研究l进行优化研究进行优化研究第二十四页,编辑于星期日:十四点 五十四分。第二章第二章 动力学动力学基本概念基本概念 2 21 1 非自由系统的约束非自由系统的约束 2 21 11 1 完整约束与非完整约束完整约束与非完整约束 2 21 12 2 定常约束与非定常约束定常约束与非定常约束 2 22 2 广义坐标和自由度广义坐标和自由度 2 22 21 1 广义坐标广义坐标 2 22 22 2 用广义坐标表示的非完整约束方程用广义坐标表示的非完整约束方程 2 22 23 3 坐标变分和自由度坐标变分和自由度 第二十五页,编辑于星期日:十四点 五十四分。21 非自由系统的约束 多多
24、个个质质点点的的集集合合可可以以组组成成一一个个质质点点系系统统,根根据据系系统统的的运运动动是是否否受受到到预预先先规规定定的的几几何何及及运运动动条条件件的的制制约约,可可以以分分为为自自由由系系统统和和非非自自由由系统系统。对对于于非非自自由由系系统统,那那些些预预先先规规定定的的、与与初初始始条条件件及及受受力力条条件件无无关关的的、限限制制系系统统的的几几何何位位置置或或(和和)速速度度的的运运动动学学条条件件称称为为约约束束。约约束束有有多多种种形形式式,这这里里只只介介绍绍其其中中两两类。类。第二十六页,编辑于星期日:十四点 五十四分。2 21 11 1 完整约束与非完整约束完整
25、约束与非完整约束 仅仅仅仅限限制制系系统统的的几几何何位位置置(也也称称位位形形)的的约约束束称称为为完完整约束。整约束。完整约束又称为完整约束又称为几何约束几何约束。若不仅限制系统的若不仅限制系统的位形位形而且还限制系统的而且还限制系统的运动速度运动速度,这样的约束称为这样的约束称为非完整约束非完整约束。完整约束与非完整约束的表达完整约束与非完整约束的表达 第二十七页,编辑于星期日:十四点 五十四分。约束方程的一般表达约束方程的一般表达式式 若若用用x xi i、y yi i、z zi i表表示示系系统统中中某某质质点点的的笛笛卡卡尔尔直直角角坐坐标标,那那么么N N个个质质点点组组成成的的
26、质质点点系系统统的的完完整整约约束束的的约约束束方程方程可写作可写作 非完整约束的约束非完整约束的约束方程取微分的形式。一个由方程取微分的形式。一个由N N个质个质点组成的系统的非完整约束方程可写作点组成的系统的非完整约束方程可写作(2.1.2)(2.1.2)f fk k(x(x1 1,y y1 1,z,z1 1,x,x2 2,y y2 2,z,z2 2,x xN N,y yN N,z,zN N,),)0 0 (k k1 1,2 2,3,3,r r r r3N3N)(2.1.1)(2.1.1)第二十八页,编辑于星期日:十四点 五十四分。图2-1 轮子的约束例例2.1 一个半径为一个半径为r r
27、的轮子沿斜面向下作纯滚动,分析轮的轮子沿斜面向下作纯滚动,分析轮子所受的约束。子所受的约束。第二十九页,编辑于星期日:十四点 五十四分。解:轮子所受的解:轮子所受的几何约束几何约束为为 (2.1.3)(2.1.3)又又运动条件的限制运动条件的限制是轮子作纯滚动时是轮子作纯滚动时P P点的速度为零点的速度为零,即即 (2.1.4)(2.1.4)或或 (2.1.5)(2.1.5)这一约束方程显然是可积分的,即这一约束方程显然是可积分的,即 (2.1.6)(2.1.6)故而轮子仍受故而轮子仍受完整约束完整约束,其约束方程为,其约束方程为(2.1.3)(2.1.3)式和式和(2(21 16)6)式。式
28、。纯滚动时轮子的约束纯滚动时轮子的约束第三十页,编辑于星期日:十四点 五十四分。例例2.2 质点质点m1和和m2由一长由一长为为l的刚性杆相连,设该的刚性杆相连,设该系统在图系统在图2-22-2所示所示xoy平面平面内运动。若要求杆中点内运动。若要求杆中点C的速度保持沿杆轴方向,的速度保持沿杆轴方向,分析该系统的约束情况。分析该系统的约束情况。图图2-2 平面运动杆的约束平面运动杆的约束第三十一页,编辑于星期日:十四点 五十四分。解:由于杆是刚性的,所以解:由于杆是刚性的,所以m1与与m2必须满足的必须满足的几何约束几何约束是是 (x1x2)2十十(y1y2)2l2 (217)而而运动约束运动
29、约束是是C点的速度必须沿杆轴方向,即点的速度必须沿杆轴方向,即平面运动杆的约束平面运动杆的约束(2(21 18)8)(2(21 18)8)式说明系统受到一个式说明系统受到一个非完整约束非完整约束。代入代入Ml,M2 2的坐标即为的坐标即为第三十二页,编辑于星期日:十四点 五十四分。我们经常遇到的系统一般我们经常遇到的系统一般是非完整系统是非完整系统。非完整约束又分。非完整约束又分为为一阶线性非完整约束一阶线性非完整约束、一阶非线性非完整约束一阶非线性非完整约束、二阶非完整二阶非完整约束约束等。等。N N个质点的系统受到个质点的系统受到k k个一阶线性非完整约束时,其约束个一阶线性非完整约束时,
30、其约束方程可以写作方程可以写作 非完整约束的类型非完整约束的类型或写成或写成 (2.1.10)(2.1.10)(2.1(2.19)9)第三十三页,编辑于星期日:十四点 五十四分。2 21 12 2 定常约束与非定常约束定常约束与非定常约束 约束方程中约束方程中不显含不显含时间时间t t的约束称为的约束称为定常约束定常约束。约。约束方程中束方程中显含显含时间时间t t的约束称为的约束称为非定常约束非定常约束。例如由方程例如由方程 所确定的约束为所确定的约束为非定常约束。非定常约束。(2(21 112)12)(2(21 111)11)例如由方程例如由方程 所确定的约束为所确定的约束为定常约束定常约
31、束。例例2.3第三十四页,编辑于星期日:十四点 五十四分。例2.3 设质点设质点M所系绳子穿过所系绳子穿过o点,如图点,如图2-3所示,绳子另一端以一匀速所示,绳子另一端以一匀速v拉动使拉动使M在在xy平面内运动。试讨论平面内运动。试讨论M的约束。的约束。图图2-3 质点质点M的非定常约束的非定常约束解解:设设M的起始位置为的起始位置为l0 0,则它到,则它到o点点的距离的距离l将随时间变化。其约束方程将随时间变化。其约束方程为为 x2 2+y2 2(l0 0-vt)2 2 (2 (21 113)13)显然,显然,M M所受的约束所受的约束是非定常约束是非定常约束。第三十五页,编辑于星期日:十
32、四点 五十四分。2 22 2 广义坐标和自由度广义坐标和自由度 第三十六页,编辑于星期日:十四点 五十四分。图图2-4 动点动点M的位置的位置2 22 21 1 广义坐标广义坐标 我们习惯于用笛卡尔直角坐标系来描述系统的几何位置即位形。我们习惯于用笛卡尔直角坐标系来描述系统的几何位置即位形。然而,根据问题的不同,然而,根据问题的不同,不一定非得不一定非得采用长度坐标参数来描述系统采用长度坐标参数来描述系统的几何位置。的几何位置。例例如如,描描述述作作平平面面运运动动的的动动点点M的的几何位置几何位置的参数可以用:的参数可以用:直角坐标直角坐标(x,y),极坐标极坐标(,r),参数参数(A,),
33、等等。,等等。这这就就是是说说,动动点点M的的几几何何位位置置可可以以用用不不同同的的参参数数组组来来描描述述,即即有有了了选选择择参参数数的的余余地地。为为此此,引引入入广广义义坐坐标标的概念。的概念。第三十七页,编辑于星期日:十四点 五十四分。广义坐标的概念广义坐标的概念 所所谓谓广广义义坐坐标标,就就是是选选择择一一组组互互相相独独立立的的参参数数q1 1,q2 2,,qn只只要要它它们们能能够够确确定定系系统统的的位位形形,而而不不管管这这些些参参数数的的几几何何意意义义如如何何。这这样样的的一一组组参参数数就就称称为为广广义义坐坐标标。因因此此,上上述述中中的的(x,y),(,r),
34、(A,)等等都都可可以以作作为为描描述述M点点的的位位形形的的广广义义坐坐标标。可可见见,广广义义坐坐标标对对于于某某一一系系统统来来讲讲不不是是唯唯一一的的,或或者者说说,可可以以任任意意选选取取。广广义义坐坐标标可可以用下面的以用下面的通式通式表示表示 ri iri i(q1 1,q2 2,qn,t)(2)(22 21)1)式式中中,ri表表示示系系统统中中第第i个个质质点点的的位位形形;qj(j1 1,2,2,n)和和t是是广广义义坐坐标标。第三十八页,编辑于星期日:十四点 五十四分。2 22 22 2 用广义坐标表示的非完整约束方程用广义坐标表示的非完整约束方程 一个由一个由N N个质
35、点组成的系统的非完整约束方程可写作微分形式。个质点组成的系统的非完整约束方程可写作微分形式。xi ixi i(q1 1,q2 2,qn,t)yi iyi i(q1 1,q2 2,qn,t)zi izi i(q1 1,q2 2,qn,t)速度的广义坐标表示速度的广义坐标表示 第三十九页,编辑于星期日:十四点 五十四分。(1)速度的广义坐标表示 设设N个个质质点点组组成成的的系系统统有有n个个广广义义坐坐标标qj(j1,1,n),),且且qjqj(t),则系统中第,则系统中第i个质点的速度是个质点的速度是式中,相应地式中,相应地 称为称为广义速度广义速度。v可以写作如下投影形式可以写作如下投影形式
36、(2(22 22)2)(2(22 23)3)第四十页,编辑于星期日:十四点 五十四分。定常系统定常系统对于定常系统,因对于定常系统,因 (2 (22 24)4)所以,所以,(2 (22 25)5)第四十一页,编辑于星期日:十四点 五十四分。图图2-5 点点M的速度的速度例2.4 空空间间中中的的一一动动点点M,若若选选取取极极坐坐标标r、为为广广义义坐坐标标,如图如图2-52-5所示,求所示,求M点在笛卡尔直角坐标系中的点在笛卡尔直角坐标系中的位置和速度。位置和速度。(2(22 27)7)(2(22 26)6)M点的位置是点的位置是(2(22 28)8)M点的速度为点的速度为于是于是M点的速度
37、为点的速度为第四十二页,编辑于星期日:十四点 五十四分。(2)(2)用广义坐标表示的非完整约束方程用广义坐标表示的非完整约束方程 一阶线性非完整约束方程已由一阶线性非完整约束方程已由(2(21 19)9)式给出:式给出:把第把第i个质点的速度的广义坐标分量代入该式得到个质点的速度的广义坐标分量代入该式得到第四十三页,编辑于星期日:十四点 五十四分。图图2-6 微分和变分微分和变分2 22 23 3 坐标变分和自由度坐标变分和自由度 坐标的变分与坐标的微分是两个不同的概念。坐标的变分与坐标的微分是两个不同的概念。设某系统运动的微分方程的解是设某系统运动的微分方程的解是坐标的变分坐标的变分则是指在
38、某一时刻则是指在某一时刻t,qj本身在约束许可条件下的任意的本身在约束许可条件下的任意的无限小增量。也就是系统的可能运无限小增量。也就是系统的可能运动动(图中的虚线所示图中的虚线所示)与真实运动在与真实运动在某时刻的差,记作某时刻的差,记作qj 既有既有不同点不同点,也有,也有共同点共同点。所谓坐标的微分是指在上式所描所谓坐标的微分是指在上式所描述的真实运动中坐标的无限小变述的真实运动中坐标的无限小变化,即经过化,即经过dtdt时间之后发生的坐时间之后发生的坐标变化标变化dqdqj j (图中实线部分图中实线部分)由于都是坐标的无限小变由于都是坐标的无限小变化,故变分也表现出微分化,故变分也表
39、现出微分的形式,并且和微分的形式,并且和微分具有具有相同的运算规则相同的运算规则。第四十四页,编辑于星期日:十四点 五十四分。自由度计算 我我们们把系把系统统独立的坐独立的坐标变标变分数称分数称为为系系统统的自由度的自由度。如果系如果系统统是是自由的自由的,则则其位形的确定要其位形的确定要3N3N个坐个坐标标。这这些坐些坐标标自然自然相互独立,其相互独立,其变变分也相互独立,故分也相互独立,故自由度自由度为为3N3N。对对于于N N个个质质点点组组成的力学系成的力学系统统,如何计算自由度呢?,如何计算自由度呢?如果系如果系统统受到受到k k个完整个完整约约束束,那么在,那么在3N3N个坐个坐标
40、标中,只有中,只有3N-k3N-k个相互个相互独立,并且它独立,并且它们们的的变变分也相互独立,故其分也相互独立,故其自由度自由度为为3N-k3N-k个。个。如果系如果系统为统为非完整系非完整系统统、假、假设该设该系系统统除了除了k个完整个完整约约束束之外,之外,还还受受到到l个非完整个非完整约约束束,该该系系统统独立的坐独立的坐标标数数为为3N-k个,但其独立的个,但其独立的坐坐标变标变分数只有分数只有3N-k-l个个 (由于由于l个微分形式个微分形式约约束的存在束的存在),故,故系系统统的的自由度自由度为为3N-k-l个个。广广义义坐坐标标数数为为n,独立的坐,独立的坐标标数数,独立的独立
41、的坐坐标变标变分数分数,系系统统的自由度的自由度之间的关系之间的关系。第四十五页,编辑于星期日:十四点 五十四分。自由度计算综综上所述,若一个系上所述,若一个系统统的的广广义义坐坐标标数数为为n,则则:完整系完整系统统:n=独立的坐独立的坐标标数数 独立的坐独立的坐标变标变分数分数 系系统统的自由度。的自由度。非完整系非完整系统统:n独立的坐独立的坐标标数数 独立的坐独立的坐标变标变分数系分数系统统的自由度。的自由度。n 系系统统的自由度的自由度第四十六页,编辑于星期日:十四点 五十四分。例例25 一一平平面面曲曲柄柄滑滑块块机机构构,A、B两两点点的的位位置置可可确确定定系系统统的的位位形形
42、,分析其自由度。分析其自由度。图图2-7第四十七页,编辑于星期日:十四点 五十四分。例例25 解解:这是一个平面机构,这是一个平面机构,A、B共有共有2 2N4个坐标,系统要满足个坐标,系统要满足3 3个完个完整约束整约束 该系统没有非完整约束,因此是一个完整系统,其自由度数为该系统没有非完整约束,因此是一个完整系统,其自由度数为4 43 31 1,独立的坐标数也是,独立的坐标数也是1 1。若选取。若选取为广义坐标,当为广义坐标,当给定时,整个给定时,整个系统的位形也就确定了。系统的位形也就确定了。(2.2.152.2.15)(2.2.162.2.16)第四十八页,编辑于星期日:十四点 五十四
43、分。作业P12 习题习题2-1,2-4第四十九页,编辑于星期日:十四点 五十四分。第五十页,编辑于星期日:十四点 五十四分。3 31 1 刚体绕定点转动的欧拉定理刚体绕定点转动的欧拉定理 3 32 2 描述刚体定点转动的解析法描述刚体定点转动的解析法第三章 刚体定点转动运动学 刚刚体定点体定点转动转动的方向余弦描述的方向余弦描述 刚刚体定点体定点转动转动的欧拉角描述的欧拉角描述 刚刚体定点体定点转动转动的广的广义义欧拉角描述欧拉角描述 第五十一页,编辑于星期日:十四点 五十四分。3 31 1 刚体绕定点转动的欧拉定理刚体绕定点转动的欧拉定理 方位和方位变化方位和方位变化设设o为刚体的固定点,刚
44、体上某为刚体的固定点,刚体上某ABo可完可完全确定刚体的方位。全确定刚体的方位。今今ABo转到转到ABO,存在一通过固定点,存在一通过固定点的轴的轴OC,当,当0A绕绕OC转过一转过一角到达角到达0A时,时,ABO与与ABO一定完全重合,一定完全重合,这种转动这种转动通常称为通常称为刚体的一次转动刚体的一次转动或或欧拉转动欧拉转动,OC即为即为一次转轴一次转轴或或欧拉转轴欧拉转轴。具有具有固定点的刚体固定点的刚体由某一由某一方位方位到另一方位的到另一方位的方位变化方位变化永永远等价于远等价于绕通过固定点的某轴绕通过固定点的某轴的一个的一个有限有限(转角转角)的转动的转动,这就是这就是刚体绕定点
45、转动的欧拉定理刚体绕定点转动的欧拉定理。第五十二页,编辑于星期日:十四点 五十四分。静锥和动锥如果将刚体的转动过程分为如果将刚体的转动过程分为若干时间间隔若干时间间隔,每一时刻欧拉,每一时刻欧拉转轴的位置显然是不同的转轴的位置显然是不同的 。在某一在某一时时刻刻t ti i,当,当时间间时间间隔隔t00时时,oci称称为刚为刚体在体在ti时时刻刻的的瞬瞬时转动轴时转动轴,平均角速度向量的极,平均角速度向量的极值值i称称为为瞬瞬时时角速度角速度向量向量。瞬瞬时转轴时转轴位置的不断位置的不断变变化化在空在空间间形成了以定点形成了以定点O为顶为顶点点的的锥锥面,称之面,称之为为静瞬静瞬时锥时锥面面,
46、简简称称静静锥锥。同同时时它在它在刚刚体内部留下了体内部留下了轨轨迹,构成了迹,构成了动动瞬瞬时时锥锥面面,它也是以,它也是以O为顶为顶点的点的锥锥面,面,简简称称动锥动锥。第五十三页,编辑于星期日:十四点 五十四分。刚刚体体绕绕定点定点转动转动的的过过程程 刚刚体体绕绕定点定点转动转动的的过过程可以程可以看成是看成是一系列以角速度一系列以角速度i绕绕瞬瞬时转动轴转动时转动轴转动的合成的合成。也也可以可以说说,刚刚体做定点体做定点转动时转动时,动动瞬瞬时锥时锥面在静瞬面在静瞬时锥时锥面上以面上以角速度角速度(t)作无滑作无滑动动的的滚动滚动,见图见图3 32 2。第五十四页,编辑于星期日:十四
47、点 五十四分。定点转动刚体上点的速度和加速度当刚体相对某动参考系以当刚体相对某动参考系以1转动而此动参考系又以转动而此动参考系又以2相对相对定参考系转动,则刚体的运动可以看成绕某个定参考系转动,则刚体的运动可以看成绕某个OC轴以角速度轴以角速度1 1十十2 2作转动,作转动,OC即为即为的方向。这就是说,刚体绕相交轴的方向。这就是说,刚体绕相交轴转动合成时,角速度的合成服从向量加法。转动合成时,角速度的合成服从向量加法。设设刚刚体体的的瞬瞬时时角角速速度度为为,则则刚刚体体上上相相对对定定点点的的向向径径为为r的的点点的速度为的速度为(3.1.1)(3.1.1)(3.1.2)(3.1.2)其中
48、,其中,为刚体的角加速度;为刚体的角加速度;称为转动加速称为转动加速度;度;称为向心加速度。称为向心加速度。第五十五页,编辑于星期日:十四点 五十四分。32 描述刚体定点转动的解析法 上一节的讨论实际上是刚体定点转动的一种简单上一节的讨论实际上是刚体定点转动的一种简单的、几何的、定性的描述,本节详细介绍刚体定点的、几何的、定性的描述,本节详细介绍刚体定点转动的定量的描述。转动的定量的描述。刚刚体定点转动的方向余弦描述刚刚体定点转动的方向余弦描述 刚刚体定点转动的欧拉角描述刚刚体定点转动的欧拉角描述 刚刚体定点转动的广义欧拉角描述刚刚体定点转动的广义欧拉角描述 第五十六页,编辑于星期日:十四点
49、五十四分。图图 3-3 i和和j 坐标系坐标系(1)方向余弦矩阵 假设以参考空间某一点假设以参考空间某一点O为原点,有两个笛卡尔直角坐标系为原点,有两个笛卡尔直角坐标系 o (简称简称i系系)和和oxyz(简称简称j系系),各坐标轴之间夹角的余弦值构成了一个方向余弦各坐标轴之间夹角的余弦值构成了一个方向余弦矩阵矩阵A,它可以表示两坐标系之间的,它可以表示两坐标系之间的空间关系空间关系。第五十七页,编辑于星期日:十四点 五十四分。图图 3-3 i和和j 坐标系坐标系两坐标系之间的两坐标系之间的空间关系空间关系 如果以如果以j j系为参考系,系为参考系,i i系是由系是由j j系绕系绕O O点点转
50、动后的结果;同理,如果以转动后的结果;同理,如果以i i系为参考系,系为参考系,j j系是由系是由i i系绕系绕O O点转动后的结果点转动后的结果.i相对相对j系的方向系的方向余弦矩阵余弦矩阵j相对相对i系的方系的方向余弦矩阵向余弦矩阵 x y z 展开:展开:第五十八页,编辑于星期日:十四点 五十四分。两矩阵之间的两矩阵之间的关系关系它们是两个正交矩阵,即它们是两个正交矩阵,即第五十九页,编辑于星期日:十四点 五十四分。矢量Q在不同空间中的表达和转换 假设在假设在j系和系和i系的原点有一系的原点有一空间向量空间向量Q(见图见图3-3-3)3)。用。用Qi(Q,Q,Q)表示表示Q在在i系中的位