资源描述:
《课程设计--牛头刨床的运动分析—matlab程序》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、牛头刨床的运动分析—Matlab程序摘要:文中建立了用解析法作牛头刨床的运动分析,基于机构的位置方程式,将位置方程对时间求一次和二次倒数,得到机构的速度和加速的方程,然后就得到位移、速度、加速度,并通过复数矢量法和矩阵法,并且借助Matlab等软件对牛头刨床进行运动分析,并对牛头刨床的运动学规律进行了仿真,为该机床的优化设计提供了理论支撑。Abstract:Thispaperestablishestheanalyticalmethodusedforshapingmotionanalysis,basedonthepositionofthemech
2、anismequations,thepositionequationontimeandoneofthefirstandtwocountdown,bethevelocityandaccelerationequationsofdisplacement,velocity,andthenget,acceleration,andbycomplexnumbervectormethodandmatrixmethod,andwiththeaidofMatlabsoftwareshapermotionanalysis,andtheshaperofthekinem
3、aticsofthemachinetoolforsimulation,optimizationdesignhasprovidedthetheorysupport.1.引言:对牛头刨床进行运动分析的任务是在已知机构尺寸及原动件运动规律的情况下,确定机构中其他构件上某些点的轨迹、位移、速度、加速度和构件的角位移、角速度及角加速度。而这些内容,不论是设计新的机器,还是为了了解现有机械的运动性能,都是十分必要的,而且它还是研究机械动力性能的必要基础和前提。对于牛头刨床来说,其刨刀在有效工作行程中应该接近等速运动,而回程速度却应高于前者速度,从而提高效率
4、。为了确定牛头刨床的设计是否满足要求,就必须对其进行运动分析。本文建立了牛头刨床的力学模型,并推导了其运动微分方程,通过数值方法求得了牛头刨床的运动学规律并开发了一款用于牛头刨床的通用运动学分析软件。2.机构结构分析:1.自由度计算:F=3n-2-F=5*3-7*2=12.结构分析:5个可动构件,两个二级杆组组成,5个转动副2个移动副,为二级机构。3.机构运动分及程序设计1.如图先建立坐标系,已知=125mm,=600mm,=150mm,原动件1的方位角~和等角速度=1rad/s如图所示:解:1.位置分析:由于有四个位置量,为了求解需要建立两个
5、封闭矢量方程。有图形ABCA可写出矢量方程:可变为:…………(1)由封闭图形CDEGC可写出矢量方程:…………(2)得到方程:由上式可得:2.速度和加速度分析式(1)和(2)对时间取一次二次导数写成矩阵形式,即得一下速度方程:=加速度方程=+二、matlab源程序:程序1clearall;clc;w1=1;l1=0.125;l3=0.600;l4=0.150;l6=0.275;l61=0.575;form=1:3601o1(m)=pi*(m-1)/1800;o31(m)=atan((l6+l1*sin(o1(m)))/(l1*cos(o1(m)
6、)));ifo31(m)>=0o3(m)=o31(m);elseo3(m)=pi+o31(m);end;s3(m)=(l1*cos(o1(m)))/cos(o3(m));o4(m)=pi-asin((l61-l3*sin(o3(m)))/l4);se(m)=l3*cos(o3(m))+l4*cos(o4(m));ifo1(m)==pi/2o3(m)=pi/2;s3(m)=l1+l6;endifo1(m)==3*pi/2o3(m)=pi/2;s3(m)=l6-l1;endA1=[cos(o3(m)),-s3(m)*sin(o3(m)),0,0;s
7、in(o3(m)),s3(m)*cos(o3(m)),0,0;0,-l3*sin(o3(m)),-l4*sin(o4(m)),-1;0,l3*cos(o3(m)),l4*cos(o4(m)),0];B1=w1*[-l1*sin(o1(m));l1*cos(o1(m));0;0];D1=A1B1;E1(:,m)=D1;ds(m)=D1(1);w3(m)=D1(2);w4(m)=D1(3);ve(m)=D1(4);A2=[cos(o3(m)),-s3(m)*sin(o3(m)),0,0;sin(o3(m)),s3(m)*cos(o3(m)),0,
8、0;0,-l3*sin(o3(m)),-l4*sin(o4(m)),-1;0,l3*cos(o3(m)),l4*cos(o4(m)),0];B2=-[