当前位置:首页>> >>


平面连杆机构的动态仿真.rar

收藏

资源目录
    文档预览:
    编号:20180802233205124    类型:共享资源    大小:1.19MB    格式:RAR    上传时间:2018-08-02
      
    29
    金币
    关 键 词:
    平面 连杆机构 动态 仿真
    资源描述:
    南昌航空大学科技学院学士学位论文1目 录1 绪论 ...............................................................................................................................(2)1.1 引言 ........................................................................................................................(2)1.2 平面连杆机构及杆组概述 ....................................................................................(2)1.3 进行杆组系统仿真的意义 ...................................................................................(3)1.4 仿真软件的发展状况与应用 ...............................................................................(3)1.5 MATLAB 概述 ......................................................................................................(3)2 7R 六杆 Ⅲ级机构运动学仿真 .........................................................................(5)2.1 曲柄原动件运动学分析 .....................................................................................(5)2.2 6RⅢ级杆组运动学分析 ...................................................................................(6)2.3 7R 六杆Ⅲ级机构 MATLAB 仿真积分模块初值的确定 ...............................(11)2.4 7R 六杆Ⅲ级机构运动学仿真模型及结果 ......................................................(16)3 7R 六杆 Ⅲ级机构动力学仿真 .......................................................................(23)3.1 曲柄原动件动力学数学模型的建立 ...............................................................(23)3.2 6RⅢ级杆组动力学数学模型的建立 ...............................................................(25)3.3 需要引用的函数 .............................................................................................(30)3.4 7R 六杆Ⅲ级机构运动学仿真模型及结果 ......................................................(32)4 结论 ...........................................................................................................................(39)参考文献 .......................................................................................................................(40)致谢 .................................................................................................................................(41)南昌航空大学科技学院学士学位论文21 绪论1.1 引言大学的四年生活,通过老师的讲解和我自己的学习,我收获了很多,我也深深的喜欢上了机械这个行业,对机械加工和制造方面尤为感兴趣,我觉得通过自己的努力和思考来改变工艺规程来提高生产效率,提高经济效益很有成就感。我所研究的课题就是给了这样的机会我可以通过我的努力来优化工艺规程,提高经济效益。此次毕业设计,是在我们学完了机械制造工艺学、工艺装备设计等课程,进行了生产实习之后,进行的一个重要的实践性环节。这要求我们把所学的工艺理论和实践知识,在实际的工艺、夹具设计中综合地加以运用,这有助与提高了我们分析和解决生产实际问题的能力,为以后从事相关的技术工作奠定的基础。1.2 平面连杆机构及杆组概述平面连杆机构是将各构件用转动副或移动副联接而成的平面机构。最简单的平面连杆机构是由四个构件组成的,简称平面四杆机构。它的应用非常广泛,而且是组成多杆机构的基础。全部用回转副组成的平面四杆机构称为铰链四杆机构。仅能在某一角度摆动的连架杆,称为摇杆。对于铰链四杆机构来说,机架和连杆总是存在的,因此可按照连架杆是曲柄还是摇杆,将铰链四杆机构分为三种基本型式:曲柄摇杆机构、双曲柄机构和双摇杆机构。在实际机械中,平面连杆机构的型式是多种多样的,但其中绝大多数是在铰链四杆机构的基础上发展和演化而成。如曲柄滑块机构、导杆机构等。任何机构都是由原动件、机架和从动件构成的系统。由于机架的自由度为零,一般每个原动件的自由度为 1,且根据运动链成为机构的条件可知,机构的自由度与原动件为应相等,所以,从动件系统的自由度数必为零。机构的从动件系统还可以进一步分解成若干个不可再分的自由度为零的构件组合,这种组合称为杆组。设n 表示活动构件数,PL 表示低副个数,根据 n 的取值不同,村级可分为Ⅱ级杆组和Ⅲ级杆组。其中Ⅱ级杆组分为 5 种:RRRⅡ级杆组、RRPⅡ级杆组、RPRⅡ级杆组、PRPⅡ 级杆组以及 RPPⅡ 级杆组。任何机构都可以看作是由若干个基本村级依次联接于原动件和机架而构成的,这就是所谓机构的组成原理。通常,把由最高级别为Ⅱ级杆组的基本杆组构成的机构称为Ⅱ级机构;把最高级为Ⅲ级杆组的基本杆组构成的机构称为Ⅲ级机构。南昌航空大学科技学院学士学位论文31.3 进行杆组系统仿真的意义系统仿真是建立在控制理论、相似理论、信息处理技术和计算机初等理论基础之上的,以计算机和其他专用物理效应设备为工具,利用系统模型对真实或假设的系统进行试验,并借助于专家的经验知识、统计数据和信息资料对实验结果进行分析研究,进而做出决策的一门综合的实验性学科。仿真技术是分析、研究各种系统,尤其是复杂系统的重要工具。随着机械行业的迅速发展,对研究、设计的机械设备越来越复杂,用于制造各种零件的材料价格越来越昂贵,不可能每一步都采取试制再修改的方法进行设计,采用仿真的方法可以在一定程度上克服这种不足的不足,降低研究成本,提高效率。而连杆机构作为常见的传动机构,对其进行运动学和动力学仿真,建立起基本杆组模块的仿真模型,无疑对日后的设计大有裨益。一般机构的运动分析,使用Quik BASIC语言或Fortran语言编写程序进行计算,其缺点“透明性”差,修改麻烦等.而用MATLAB 对机构进行运动仿真,利用MATLAB的simulink仿真模型的数据可视化的特点,就可以很容易观察到运动参数是如何变化的,极其简便.同时,用MATLAB建立和修改仿真模型具有方便、快捷、很容易扩展等优点.MATLAB仿真求解器提供很多解不同微分方程的方法,可以根据不同的微分方程类型选择相应的求解方法.机构的动力学分析,由已知工作阻力,求出运动副的约束反力和驱动力(或力矩) ,为选择和设计轴承和零部件强度的计算及选择原动机提供理论依据。1.4 仿真软件的发展状况与应用早期的计算机仿真技术大致经历了几个阶段:20 世纪 40 年代模拟计算机仿真;50 年代初数字仿真;60 年代早期仿真语言的出现等。80 年代出现的面向对象仿真技术为系统仿真方法注入了活力。我国早在 50 年代就开始研究仿真技术了,当时主要用于国防领域,以模拟计算机的仿真为主。70 年代初开始应用数字计算机进行仿真[4] 。随着数字计算机的普及,近 20 年以来,国际、国内出现了许多专门用于计算机数字仿真的仿真语言与工具,如 CSMP,ACSL , SIMNOM, MATLAB/Simulink, Matrix/System Build, CSMP-C 等。1.5 MATLAB 概述MATLAB 是国际上仿真领域最权威、最实用的计算机工具。它是 MathWork 公司于 1982 年推出的一套高性能的数值计算和可视化数学软件,被誉为“巨人肩上的南昌航空大学科技学院学士学位论文4工具” 。MATLAB 是一种应用于计算技术的高性能语言。它将计算,可视化和编程结合在一个易于使用的环境中,此而将问题解决方案表示成我们所熟悉的数学符号,其典型的使用包括:.数学计算.运算法则的推导.模型仿真和还原.数据分析,采集及可视化.科技和工程制图.开发软件,包括图形用户界面的建立MATLAB 是一个交互式系统,它的基本数据元素是矩阵,且不需要指定大小。通过它可以解决很多技术计算问题,尤其是带有矩阵和矢量公式推导的问题,有时还能写入非交互式语言如 C 和 Fortran 等。MATLAB 的名字象征着矩阵库。它最初被开发出来是为了方便访问由LINPACK 和 EISPAK 开发的矩阵软件,其代表着艺术级的矩阵计算软件。MATLAB 在拥有很多用户的同时经历了许多年的发展时期。在大学环境中,它作为介绍性的教育工具,以及在进阶课程中应用于数学,工程和科学。在工业上它是用于高生产力研究,开发,分析的工具之一。Simulink 概述 Simulink 是用于仿真建模及分析动态系统的一组程序包,它支持线形和非线性系统,能在连续时间,离散时间或两者的复合情况下建模。系统也能采用复合速率,也就是用不同的部分用不同的速率来采样和更新。Simulink 提供一个图形化用户界面用于建模,用鼠标拖拉块状图表即可完成建模。在此界面下能像用铅笔在纸上一样画模型。相对于以前的仿真需要用语言和程序来表明不同的方程式而言有了极大的进步。Simulink 拥有全面的库,如接收器,信号源,线形及非线形组块和连接器。同时也能自己定义和建立自己的块。模块有等级之分,因此可以由顶层往下的步骤也可以选择从底层往上建模。可以在高层上统观系统,然后双击模块来观看下一层的模型细节。这种途径可以深入了解模型的组织和模块之间的相互作用。南昌航空大学科技学院学士学位论文52 7R六杆Ⅲ级机构运动学仿真2.1 曲柄原动件运动学分析2.1.1 曲柄原动件运动学数学模型的建立如图1所示,在复数坐标系中,曲柄AB复向量的模 rj为常数、幅角θj为变量,通过转动副A与机架连接,转动副A的复向量的模r i为常量、幅角θi为常量,曲柄AB端点 B的位移、速度和加速度的推导如下:(2.1)jijijjji errABe,将方程 2.1 两边对时间 t 求两次导数得:(2.2))2/()2/( jj jje由式 2.2 写成矩阵形式有(2.3) )sin()2/sin(cocoImRe jjjj jjjj rrB2.1.2 曲柄 MATLAB 运动学仿真模块 M 函数根据式(2.3)编写曲柄原动件 MATLAB 的 M 函数如下:function y=crank(x)%x(1)=rj 曲柄杆长%x(2)=thetaj 曲柄与水平方向夹角%x(3)=dthetaj 曲柄角速度%x(4)=ddthetaj 曲柄角加速度%y(1)=Re[ddB] 转动副 B 加速度实轴分量图 1 曲柄的复数坐标系南昌航空大学科技学院学士学位论文6%y(2)=Im[ddB] 转动副 B 加速度虚轴分量ddB=[x(1)*x(4)*cos(x(2)+pi/2)+x(1)*x(3)^2*cos(x(2)+pi);x(1)*x(4)*sin(x(2)+pi/2)+x(1)*x(3)^2*sin(x(2)+pi)];y=ddB;此函数模块用于计算转动副 B 的加速度的水平分量和垂直分量输入参数为曲柄的长度、角位移、角速度和角加速度;输出参数为曲柄端部(转动副 B)的加速度的水平分量和垂直分量。2.2 6RⅢ级杆组运动学分析2.2.1 6RⅢ级杆组运动学数学模型的建立如图2所示,在复数坐标系中,由3个外转动副(B,C,D) 和3个内转动副(E,F,G),4 个构件(BE,CF,DG和EFG)构成1个6RⅢ级杆组,构件BE,CF,DG的长度分别为r i, rj, rk构件EFG的3个边为e,f ,g,方向如图所示,规定所有复向量与实轴正方向逆时针夹角为θ,并用相应的下标来区别,用 B,C,D,E,F和G 分别表示该转动副的复数坐标,则各个构件的运动参数推导如下:(2.4)kjijjierDGCFerBE(2.5)gfejjjeEFgfe图 2 6RⅢ级杆组的位置参数南昌航空大学科技学院学士学位论文7将式(2.4)代入式(2.5)并整理得:(2.6)BDfereCgrfki ekjji jjj jjjj 将式(2.4)至式 (2.6)三式合并成矩阵得:         BDCfeerer eer gefeerer eer gegfekji jjkjijj jjigfekji jjkji jj jji fi kji fi kji222 )()()()()( )()( )2/()2/()2/( )2/()2/( )2/()/( 000 000 000     (2.7)将式(2.7)展开整理得:    0)2sin(0)2sin(0)2sin( cococo 0)2sin()i()2cos(0 0co2s )2sin(00)()2sin( co02cosco    fk fekiiii gjjjj rrrrrr南昌航空大学科技学院学士学位论文8  BDCB frr er rrgfekji fkii jj k iijjiigfekjiImRIe 0)cos(0)cos(0)sn(co )sin(i)si(cc )cs()s( o00222 (2.8)点 E,F,G 的加速度分别为      222)sin(co)2sin(coImReIe )sin(co)2sin(coIeIe )sin(co)2sin(coImReIe kkkkk jjjjjj iiiiii rrDG rrCF rrB (2.9)2.2.2 6RⅢ级杆组 MATLAB 运动学仿真模块 M 函数根据式(2.9)编写 6RⅢ级杆组 MATLAB 的 M 函数如下:function y=R6ki(x)%x(1)=ri BE 杆长%x(2)=rj CF 杆长%x(3)=rk DG 杆长南昌航空大学科技学院学士学位论文9%x(4)=e FG 杆长%x(5)=f GE 杆长%x(6)=g EF 杆长%x(7)=theta-i BE 杆与水平方向夹角%x(8)=theta-j DF 杆与水平方向夹角%x(9)=theta-k DG 杆与水平方向夹角%x(10)=theta-e FG 杆与水平方向夹角%x(11)=theta-f GE 杆与水平方向夹角%x(12)=theta-g EF 杆与水平方向夹角%x(13)=dtheta-i BE 杆角速度%x(14)=dtheta-j EF 杆角速度%x(15)=dtheta-k DG 杆角速度%x(16)=dtheta-e GE 杆角速度%x(17)=dtheta-f GE 杆角速度%x(18)=dtheta-g EF 杆角速度%x(19)=Re[ddB] 转动副 B 加速度实轴分量%x(20)=Im[ddB] 转动副 B 加速度虚轴分量%x(21)=Re[ddC] 转动副 C 加速度实轴分量%x(22)=Im[ddC] 转动副 C 加速度虚轴分量%x(23)=Re[ddD] 转动副 D 加速度实轴分量%x(24)=Im[ddD] 转动副 D 加速度虚轴分量%y(1)=ddtheta-i BE 杆角加速度%y(2)=ddtheta-j CF 杆角加速度%y(3)=ddtheta-k DG 杆角加速度%y(4)=ddtheta-e FG 杆角加速度%y(5)=ddtheta-f GE 杆角加速度%y(6)=ddtheta-f EF 杆角加速度%y(7)=Re[ddE] 转动副 E 加速度实轴分量%y(8)=Im[ddE] 转动副 E 加速度虚轴分量%y(9)=Re[ddF] 转动副 F 加速度实轴分量%y(10)=Im[ddF] 转动副 F 加速度虚轴分量南昌航空大学科技学院学士学位论文10%y(11)=Re[ddG] 转动副 G 加速度实轴分量%y(12)=Im[ddG] 转动副 G 加速度虚轴分量a=[-x(1)*cos(x(7)+pi/2) x(2)*cos(x(8)+pi/2) 0 0 0 -x(6)*cos(x(12)+pi/2);-x(1)*sin(x(7)+pi/2) x(2)*sin(x(8)+pi/2) 0 0 0 -x(6)*sin(x(12)+pi/2);0 -x(2)*cos(x(8)+pi/2) x(3)*cos(x(9)+pi/2) -x(4)*cos(x(10)+pi/2) 0 0;0 -x(2)*sin(x(8)+pi/2) x(3)*sin(x(9)+pi/2) -x(4)*sin(x(10)+pi/2) 0 0;x(1)*cos(x(7)+pi/2) 0 -x(3)*cos(x(9)+pi/2) 0 -x(5)*cos(x(11)+pi/2) 0;x(1)*sin(x(7)+pi/2) 0 -x(3)*sin(x(9)+pi/2) 0 -x(5)*sin(x(11)+pi/2) 0];c=-[-x(1)*cos(x(7)+pi) x(2)*cos(x(8)+pi) 0 0 0 -x(6)*cos(x(12)+pi);-x(2)*sin(x(7)+pi) x(2)*sin(x(8)+pi) 0 0 0 -x(6)*sin(x(12)+pi);0 -x(2)*cos(x(8)+pi) x(3)*cos(x(9)+pi) -x(4)*cos(x(10)+pi) 0 0;0 -x(2)*sin(x(8)+pi) x(3)*sin(x(9)+pi) -x(4)*sin(x(10)+pi) 0 0;x(1)*cos(x(7)+pi) 0 -x(3)*cos(x(9)+pi) 0 -x(5)*cos(x(11)+pi) 0;x(1)*sin(x(7)+pi) 0 -x(3)*sin(x(9)+pi) 0 -x(5)*sin(x(11)+pi) 0];b1=c*[x(13)^2;x(14)^2;x(15)^2;x(16)^2;x(17)^2;x(18)^2];b2=[x(19)-x(21);x(20)-x(22);x(21)-x(23);x(22)-x(24);x(23)-x(19);x(24)-x(20)];b=b1+b2;ddtheta=inv(a)*b;y(1)=ddtheta(1);y(2)=ddtheta(2);y(3)=ddtheta(3);y(4)=ddtheta(4);y(5)=ddtheta(5);y(6)=ddtheta(6);y(7)=x(19)+x(1)*ddtheta(1)*cos(x(7)+pi/2)+x(1)*x(13)^2*cos(x(7)+pi);y(8)=x(20)+x(1)*ddtheta(1)*sin(x(7)+pi/2)+x(1)*x(13)^2*sin(x(7)+pi);y(9)=x(21)+x(2)*ddtheta(2)*cos(x(8)+pi/2)+x(2)*x(14)^2*cos(x(8)+pi);y(10)=x(22)+x(2)*ddtheta(2)*sin(x(8)+pi/2)+x(2)*x(14)^2*sin(x(8)+pi);y(11)=x(23)+x(3)*ddtheta(3)*cos(x(9)+pi/2)+x(3)*x(15)^2*cos(x(9)+pi);y(12)=x(24)+x(3)*ddtheta(3)*sin(x(9)+pi/2)+x(3)*x(15)^2*sin(x(9)+pi);这个模块用于求Ⅲ级杆组中各杆的加速度的水平及垂直分量。输入参数为构件
    展开阅读全文
    1
      金牌文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    0条评论

    还可以输入200字符

    暂无评论,赶快抢占沙发吧。

    关于本文
    本文标题:平面连杆机构的动态仿真.rar
    链接地址:http://www.gold-doc.com/p-37396.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们
    copyright@ 2014-2018 金牌文库网站版权所有
    经营许可证编号:浙ICP备15046084号-3
    收起
    展开