1、基于 Mathematica 的机构运动分析研究研究小组成员：龚加庆 付忠贵金晓东 卢航赵航基于 Mathematica 的机构运动分析研究摘要：机构运动分析的任务是在一直尺寸及原动件运动规律的情况下，确定机构中其他构件上某些点的轨迹、位移、速度及加速度和构件角位移、角速度、角加速度。这些分析无论是设计新机械，还是为了了解现有机构的性能，都是十分重要的，而且是究机械动力性能的重要前提。现有比较重要的两种结构运动分析方法：图解法和解析法。基于解析法的精确性同时利用计算机的标准运算等优点，我们选择利用 Mathematica 对机构进行运动的分析。我们利用解析法对要研究的四连杆机构建立方程，给定原
2、动件及构件的尺寸参数，在通过 Mathematica 编程来求解相关构件上一点的运动规律，并对其进行画图。然后再利用 admas 软件进行相应的仿真并画出对应点的运动参量的图线来验证我们利用Mathematica 求解出的点的运动规律。结果两款软件给出的对应点的运动参量（位移，速度，加速度）图线完全一致。因此我们得出结论基于 Mathematica 的机构运动分析研究是十分科学可靠的，并且有很强的实用性，有很大的利用价值。Abstract: The organization is the task of motion analysis has been moving pieces of the
3、 original size and movement of thecircumstances, determine the organization of certain other components in the locus ofpoints, displacement, velocity and acceleration and angular displacement components, angular velocity, angular acceleration. The analysis of both the design of new machinery, or to
4、understand the performance of existing institutions are very important to study mechanical power and is an important prerequisite for the performance. Two of the more important structure of the existing motion analysis: graphical and analytical method. Analytical method based on the accuracy of the
5、standard operations while leveraging the advantages of the computer, we chose to use Mathematica on the bodies of sport. We use analytical method to study the establishment of four-bar linkage equation, given the original size of moving parts and component parameters, related by Mathematica program
6、to solve a point on the movement of components, and its drawing. Then use admas software simulation and draw the corresponding point of the motion parameters corresponding to the graph to verify that we use Mathematica to solve the movement of the point. The results of two corresponding points given
7、 software movement parameters (displacement, velocity, acceleration) of lines exactly the same. Therefore, we conclude based on kinematic analysis of Mathematicas study is very scientific and reliable, and has a strong practical, there is great value in use.引言：机构运动分析的任务是在一直尺寸及原动件运动规律的情况下，确定机构中其他构件上某
8、些点的轨迹、位移、速度及加速度和构件角位移、角速度、角加速度。这些分析无论是设计新机械，还是为了了解现有机构的性能，都是十分重要的，而且是究机械动力性能的重要前提。现在我们我们研究的内容是基于 Mathematica 的机构运动分析研究。因为 Mathematica 矢量分析方法利用了复数的运算十分简便的优点，不仅对任何机构包括复杂的连杆机构运动分析和动力分析，而且用来进行机构的运动的综合。Mathematica 的用户群中最主要的是科技工作者和其它专业人士。但是，Mathematica 还被广泛地用于教学中。而且近些年很多机械公司再设计机构时也有很多利用 Mathematica 进行仿真或者
9、利用其对运动轨迹的设计等。我们课题内容是利用 Mathematica 对四连杆机构的中的一根杆件上一点进行分析，并得出其运动位移，速度，加速度，然后利用 admas 或 working model 软件进行仿真来验证我们的研究结果。Mathematica 简介Mathematica 是一款科学计算软件，很好地结合了数值和符号计算引擎、图形系统、编程语言、文本系统、和与其他应用程序的高级连接。很多功能在相应领域内处于世界领先地位，截至 2009 年，它也是为止使用最广泛的数学软件之一。Mathematica 的发布标志着现代科技计算的开始。自从 20 世纪 60 年代以来，在数值、代数、图形、和
10、其它方面应用广泛，Mathematica 是世界上通用计算系统中最强大的系统。自从 1988 发布以来，它已经对如何在科技和其它领域运用计算机产生了深刻的影响。1、 Mathematica 软件的主要功能（1） 符号运算Mathematica 是美国 Wolfram Research 公司开发的著名数学软件，以符号运算为主，能像人一样进行带字母的运算，得到的是精确的结果。符号运算功能可以分成 4 大类：（A） 初等数学可以进行各种数和初等函数式的计算与化简。（B） 微积分可以求极限、导数（包括高阶导数和偏导数等） 、不定积分和定积分（包括多重积分） ，将函数展成幂级数，进行无穷级数求和及积分变
11、换。（C） 线性代数可进行行列式的计算、矩阵的各种运算（加法、乘法、求逆矩阵等） 、解线性方程组、求特征值和特征向量、进行矩阵分解。（D ） 解方程组能解各类方程组（包括微分方程组） 。（2） 数值计算Mathematica 的数值计算更具科学性。它与通常的数值计算程序有所不同，允许用户指定任意精度，例如能轻而易举地求出 的 300 位近似值；也可以求出 1000！和 210000 等整数的准确值，令人惊叹。Mathematica 具有众多的数值计算函数，能满足线性代数、插值与拟合、数值积分、微分方程数值解、求极值、线性规划及概率统计等方面的常用计算需求（3） 绘图它的绘图功能也很出色，能绘制
12、各种二维和三维彩色图形，自动化程度很高。（4） 编程在 Mathematica 中，用户可以自己编写各种程序（文本文件） ，开发新的功能。用户开发的功能可以在软件启动时被调入，与软件本身的功能一样使用。Mathematica 4.0 版自带一百多个专门的程序包，都是另外编写的程序文件，补充并完善了 Mathematica 的功能。Mathematica 良好的底层环境极具吸引力，为数学编程提供了广阔天地。2、应用 mathematica 的优点（1）具有简单易学的交互式方法（2）具有非常强大的数值计算能力（3）具有非常强大的符号计算能力，能够进行各种复杂的数学计算，而且速度很快，准确，节约我们
13、很多计算时间（4）具有人工智能列表处理能力（5）结构化设计语言（6）绘图功能也很出色，能绘制各种二维和三维彩色图形，自动化程度很高。3、mathematica 在工程和科学研究中的应用情况主要使用者是从事理论研究的数学工作者和其它科学工作者、以及从事实际工作的工程技术人员。但是，Mathematica 还被广泛地用于教学中。从高中到研究生院的数以百计的课程都使用它。此外，随着学生版的出现，Mathematica 已经在全世界的学生中流行起来，成为了一个著名的工具。广泛的应用-一些典型例子1 处理含有上百万语汇的复杂符号计算 2 上载、分析、和视觉化数据 3 解数值和符号方程、微分方程、和最小化
14、问题 4 建立各类数值模型和仿真，包括从简单的控制系统到星系碰撞、金 5 融衍生物、复杂生物系统、化学反应、环境影响研究、和粒子加速器中的磁场 6 促进工程公司和金融企业迅速发展的应用程序 7 制作高质量的互动技术报告或论文，用于电子发行和印刷发行 8 向各层次的学生阐述数学或科学概念，可应用于从小学生至研究生Mathematica 可以用于解决各种领域的涉及复杂的符号计算和数值计算的问题。它代替了许多以前仅仅只能靠纸和笔解决的工作，这种思维和解题工具的革新可能对各种研究领域和工程领域产生深远的影响。机构的结构分析下图是我们研究的机构图机构结构分析根据结构图，我们画的其结构简图如下图1、机构结
15、构分析从简图中可以知道，此机构共有 5 个构件，由原动件 a，构件 b 、c、e ，以及机架 d 组成，其中机架 d 分别与原动件 a，构件 c 分别构成一个转动副，原动件 a 和构件 b 构成一个转动副，构件 b 和构件 c 和构件 e 构成二个转动副。2、 自由度计算（1）计算公式F3 n(2 )lph式中：n 为机构的活动构件数目；为机构的低副数目；lp为机构的高副数目。hF3X3(2X4 0)=1故其自由度为 1。4 机构运动分析及程序设计1、位置分析由结构简图，用分析法有由位置分析列机构矢量封闭方程 dcba， 消去求 解 adcb cos2)cos(2)cos(222 adaa 0
16、sin 22 bdac令 A= sin2acB= dcosC= cos222 acb则公式简化为 0cssinCBA解得 CBtg2222同理可以求的 程序见程序2 机构的速度分析对 求导dcba得ttt bac213 eee对 点积得2e 2321 ett c对 点积得3 03231 eett ba化解得 )sin()sin( ac)sin()sin( ba推出 )sin(cac)si(bab3 机构的加速度分析对 求导得ttt bac213 eee tnntn bbac 221332 ee 程序见程序对 点积得2e 222122323 eeee nbnatcnc 化简 bacc 22)os
17、( )sin( 得 )sin()cos()c( 222 cbac同理对 点积得3e)sin( )cos()cos( 222 b cbab l程序见程序=xr_:=(2+222+2+ 442428+422+422 4222428+422+422+ 4222428+422+422+ 4222428+422+422 832428+422+422+(22(4242+42+428)2242(4222+4+222222+4+222222222+443+42+42434222+4222)(428+422+422)/(428+422+422)/(2+2)+(2+222+2+ 442428+422+422 42
18、22428+422+422+ 4222428+422+422+ 4222428+422+422 832428+422+422+(22(4242+42+428)2242(4222+4+222222+4+222222222+443+42+42434222+4222)(428+422+422)/(428+422+422)/(2+2)/(2+222+2+ 442428+422+422 4222428+422+422+ 4222428+422+422+ 4222428+422+422 832428+422+422+(22(4242+42+428)2242(4222+4+222222+4+22222222
19、2+443+42+42434222+4222)(428+422+422)/(428+422+422)/(2+2)+)+=zr_:=(2+222+2+ 442428+422+422 4222428+422+422+ 4222428+422+422+ 4222428+422+422 832428+422+422+(22(4242+42+428)2242(4222+4+222222+4+222222222+443+42+42434222+4222)(428+422+422)/(428+422+422)/(2+2)+(2+222+2+ 442428+422+422 4222428+422+422+ 4
20、222428+422+422+ 4222428+422+422 832428+422+422+(22(4242+42+428)2242(4222+4+222222+4+222222222+443+42+42434222+4222)(428+422+422)/(428+422+422)/(2+2)/(2+222+2+ 442428+422+422 4222428+422+422+ 4222428+422+422+ 4222428+422+422 832428+422+422+(22(4242+42+428)2242(4222+4+222222+4+222222222+443+42+4243422
21、2+4222)(428+422+422)/(428+422+422)/(2+2)+)+上面为程序中 E 点的位置在 x-z 的坐标，即位置矢量。根据上面的坐标分别对 角求导便可以得到对应方向的速度，在对速度求导就可以得到加速度。 （我们球的的速度和加速度是对 求导的而不是对t 求导的因为原动件是匀角速度运动的，所以我们只需进行对 的求导即可，只不多时相应的速度差了角速度 这个系数，而加速度差了 2）=m()=: =n()=:x, z,=Dm(), =Dn( ), xa za球出了相应变量的公式后，我们可以利用画图函数对其进行画图分析。机构分析的结果：为了配合 admas 软件的方针验证，我们挑
22、选了一组参数 h=1.5 mod=0 a=1 b=2.6 c=2.5d=2.7 来进行计算 E 点的 x，z 向的分速度和分加速度图像为:xzAxAz为了对分析结果进行验证我们还利用 admas 对同一组参数进行仿真得到的 E 点的 x，z 向的分速度和分加速度图像为:xzAxAz对这两种软件分析的结果进行比较，他们的曲线走势完全一样，由此得出我们基于 mathmatica 软件的机构运动分析结果是正确可靠的。参考文献：【1】： 美蒂芬 沃尔夫雷姆 . 特别版 MATHEMATICA 全书 M.西安： 西安交通大学出版社，2002：96-130【2】：孙恒，陈作模，葛文杰 .机械原理M.高等教
23、育出版社，2005（7）： 30-60【3】：薛定宇，陈阳泉 .基于 mathmatica 的系统仿真技术与应用M.北京：清华大学出版社，2002.【4】李军，邢俊文，谭恩浩 .ADAMS 实力教程M. 北京：北京理工出版社，2002.研究感想本次课题研究确实让我们的能力得到了锻炼，mathmatica 这款软件我们前都没有接触过。正因为这样大家对他的兴趣就更大了，通过老师在课堂上我们做了精炼的介绍，我们对其实用方法有了初步的了解，在课下我们也到图书馆还有网站上查找相关资料，对mathmatica 的基本功能及使用方法有了初步认识之后，我们又通过对老师给的程序的分析，逐句逐函数的解释，对每个函
24、数的功能及用途有了很好的掌握。准备工作做好了，我们就开始了我们机构的运动研究，通过某一个函数就可以得到我们想要的结果，在收获中体验了研究的乐趣。于此同时我也对研究工作有了初步的认识，它应该是真正有我们小组成员亲自动手，自己设计得出的实验结果，与以往做课题报告时对已有的文献复制粘贴是有着本质的区别的。做研究一定要认真，要有一定的理论基础，研究过程可以学到很多的方法，比如比较法，验证法等。总结来说本次研究对我们的能力的提高有了一定的促进作用，我们通过研究掌握了 mathmatica 和 admas 软件的使用方法。这是一个宝贵的财富，我们以后也会更好的利用它，在小学期我们会利用mathmatica 进行更加复杂的设计和仿真，也希望在这个过程中有更多的收获。结论： 总结用 mathematica 进行机构运动分析的优越性结果比较情况
