资源描述：
《湖大机械院牛头刨床matlab运动分析》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、HUNANUNIVERSITY标题:牛头刨床运动分析-matlab姓名:杨国周学号:201204010620专业班级：12级机自六班学院：机械与运载工程学院二零一四年六月一、工作原理及工艺动作过程牛头刨床是一种用于切削平面的加工机床，它是依靠刨刀的往复运动和支撑并固定工件的工作台的单向间歇移动来实现对平面的切削加工。刨刀向左运动时为切削工件，向右运动时为空回。且切削时刨刀的移动速度低于空行程速度，即刨刀具有急回现象。刨刀可随小刀架作不同进给量的垂直进给；安装工件的工作台应具有不同进给量的横向进给，以完成平面的加工，工作台还应具有升降功能，以适应不同高度的工件加工。二、matlab简介ATL
2、AB是一个包含大量计算算法的集合。其拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能。在通常情况下，可以用它来代替底层编程语言，在计算要求相同的情况下，使用MATLAB的编程工作量会大大减少。函数所能解决的问题其大致包括矩阵运算和线性方程组的求解、微分方程及偏微分方程的组的求解、符号运算、傅立叶变换和数据的统计分析、工程中的优化问题、稀疏矩阵运算、复数的各种运算、三角函数和其他初等数学运算、多维数组操作以及建模动态仿真等。三、CAD建立主运动模型各构建的尺寸为：l1=125mm, l3=600mm，l4=150mm,原动件1的方位角==-?和等角速度w1=1r
3、ad/s。用矩阵法求该机构中各从动的方位角、角速度、和角加速度以及E点的位移、速度和加速度的运动曲线图。四、题意分析建立以C为原点，以CG为y轴，以垂直于CG并过C的线为x轴的平面坐标系。记AB为向量l1,CB为向量S3，CD为向量l3，DE为向量l4，GE为向量SE,CA为向量l6，CG为向量l7。写出投影方程：S3*cosθ3=l1*cosθ1；S3*sinθ3=l6+l1*sinθ1；l3*cosθ3+l4cosθ4-Se；l3sinθ3+l4sinθ4=l7由以上各式即可求得S3、θ3、θ4、Se四个运动变量，而滑块2的方位角θ2=θ3，然后，分别将上列各式对时间取一次二次导数，并
4、写成矩阵形式，即得一下速度和加速度方程：五、将矩阵运算转化为Matlab源程序clearall;clcl1=125;l3=600;l4=150;l6=275;l7=575;w1=1;y1=0:2*pi/3600:2*pi;s3=sqrt(l6^2+l1^2-2*l1*l6*cos((pi/2)+y1));y3=acos(l1*cos(y1)./s3);y4=pi-asin((l7-l3*sin(y3))/l4);se=l3*cos(y3)+l4*cos(y4);fori=1:3601a11(i)=cos(y3(i));a12(i)=-s3(i)*sin(y3(i));a13(i)=0;a1
5、4(i)=0;a21(i)=sin(y3(i));a22(i)=s3(i)*cos(y3(i));a23(i)=0;a24(i)=0;a31(i)=0;a32(i)=-l3*sin(y3(i));a33(i)=-l4*sin(y4(i));a34(i)=-1;a41(i)=0;a42(i)=l3*cos(y3(i));a43(i)=l4*cos(y4(i));a44(i)=0;b11=-l1*sin(y1(i));b21=l1*cos(y1(i));b23=0;b24=0;b=[b11;b21;b23;b24];a1=[a11(i),a12(i),a13(i),a14(i)];a2=[a2
6、1(i),a22(i),a23(i),a24(i)];a3=[a31(i),a32(i),a33(i),a34(i)];a4=[a41(i),a42(i),a43(i),a44(i)];a=[a1;a2;a3;a4];v=ab;vv1(i)=v(1,1);vv2(i)=v(2,1);vv3(i)=v(3,1);vv4(i)=v(4,1);aa11(i)=-vv2(i)*sin(y3(i));aa12(i)=-vv1(i)*sin(y3(i))-vv2(i)*s3(i)*cos(y3(i));aa13(i)=0;aa14(i)=0;aa21(i)=vv2(i)*cos(y3(i));aa2
7、2(i)=vv1(i)*cos(y3(i))-s3(i)*vv2(i)*sin(y3(i));aa23(i)=0;aa24(i)=0;aa31(i)=0;aa32(i)=-l3*vv2(i)*cos(y3(i));aa33(i)=-l4*vv3(i)*cos(y4(i));aa34(i)=0;aa41(i)=0;aa42(i)=-l3*vv2(i)*sin(y3(i));aa43(i)=-l4*vv3(i)*sin(y4
