牛头刨床机构运动分析
牛头刨床运动分析
连杆机构运动分析机构运动分析的解析法有多种，其中比较常用的有矢量方程解析法、矩阵法和复数矢量法等。
用解析法作机构运动分析时，首先建立机构的位置方程，然后将其对时间求一次、二次导数，可以得到机构的速度方程和加速度方程，完成运动分析。
例3-1 图示为一牛头刨床的结构简图。
设已知各构件的尺寸为：，，，，。
并知原动件1 的方位角和等角速度。
求导杆3 的方位角，角速度及角加速度3ε和刨头5上点E 的位移，速度和加速度。
解：该牛头刨床为一个六杆机构。
先建立一直角坐标系如图，并标出各杆矢及各杆矢的方位角。
其中共有四个未知量、、、。
为求解需建立两个封闭矢量方程，为此需利用两个封闭图ABCA 及CDEGC 。
（1）求导杆3 的角位移，角速度和角加速度，由封闭形ABCA 可得写成复数形式为313126θθπi i ie s e l e l =+ （a ）展开得解上述两式可得因式中分子分母均为正，故知在第一象限。
式（a ）对时间t 求导，注意为变量，有33133311θθθωωi i i e dtds ie s ie l += （b ） 展开后可得m/srad/s （逆时针方向）再将式（b ）对时间t 求导，则有 3333133232233332112θθθθθωωεωi i i i i ie dt dse dts d es ies e l ++-=- （c ）展开后可求得1471.0/]2)s i n ([333132113=--=s dtds l ωθθωε rad/s 2 (逆时针方向)m/s 2其方向与相反。
（2）求刨头上点E 的位移，速度和加速度。
由封闭形CDEGC 可得写成复数形式为 E ii i s el el e l +=+2'64343πθθ （d ）展开得解之得由机构简图知在第二象限，故=175.3266º,而m式（d ）对时间t 求导可得dtds ie l ie l Ei i =+434433θθωω （e ） 解之得rad/s (逆时针方向)m/s其方向与相反。
牛头刨床机构方案分析
牛头刨床机构方案分析一、题目说明（图a,b）所示为两种牛头刨床主机构的运动简图，已知，l1=0.1m,l0=0.4m,l3=0.75m,l4=0.15m,l y=0.738m,l′3=0.375m,a=0.05m,b=0.15,c= 0.4m,d=0.1m。
只计构件3、5的质量，其余略去不计，m3=30kg,J s3=0.7kg·m2,m5=95kg。
工艺阻力Q如图所示，Q=9000N。
主轴1的转速为60r/min(顺时针方向)，许用运转不均匀系数[δ]=0.03。
二、内容要求与作法1.进行运动分析，绘制出滑枕5的位移、速度和加速度线图。
2.进行动态静力分析，绘制固定铰链处的反力矢端图和平衡力矩T d的线图。
3．以主轴为等效构件，计算飞轮转动惯量。
4．对两种机构方案进行综合评价，主要比较如下几项内容：工作行程中滑枕 5 的速度均匀程度。
固定铰链处的反力大小及方向变化。
平衡力矩平均值及波动情况。
飞轮转动惯量大小。
上机前认真读懂所用子程序，自编主程序，初始位置取滑枕 5 的左极限位置。
主程序中打开一数据文件“DGRAPS”，写入需要显示图形的数据。
三、课程设计说明书内容上机结束后，每位学生整理出课程设计说明书一份，其内容应包括：1.机构简图和已知条件。
2.滑枕初始位置及行程H的确定方法。
3.杆组的拆分方法及所调用的杆组子程序中虚参与实参对照表。
4.飞轮转动惯量的计算方法。
5.自编程序中主要标识符说明。
6.自编程序及计算结果清单。
7.各种线图：①滑枕的位移、速度和加速度线图，②平衡力矩线图③固定铰链处反力矢端图④等效转动惯量，等效阻力矩，等效构件角速度线图。
8.对两种方案的比较，评价。
9.以一个位置为例，用图解法做机构的运动分析，与解析法计算结果比较误差。
10．主要收获与建议。
指导教师参考上述内容提出具体要求，学生按照指导教师的要求书写并制订成册。
机械原理课程设计---牛头刨床主体机构的分析与综合
机械原理课程设计---牛头刨床主体机构的分析与综合1 课程设计的目的和任务牛头刨床是常见的一种金属加工机床如图1所示。
其主体机构的机构运动简图有多种形式，图2所示的是常用的五种主体机构的示意图。
图 1 牛头刨床图2 牛头刨床的主体结构机构运动简图课程设计的内容包括：1）牛头刨床主传动系统总体传动方案的设计构思一个合理的传动系统。
它可将电机的高速转动（1440 转/分）变换为安装有刨刀的滑枕5 的低速往复移动（要求有三挡速度：60，95，150 次/分）。
其中，将转动变为移动的装置（主体机构）采用图2 所示的连杆机构。
在构思机构传动方案时，能做到思路清晰，各部分的传动比分配合理，最后在计算机上绘出主传动机构的原理示意图。
2）牛头刨床主体机构的尺度综合已知数据如表1所示图中的参数如图3所示。
图3参数表达形式表13）牛头刨床主体机构的运动分析根据已定出的主体机构的尺度参数，按曲柄处于最低转速、滑枕处于最大行程的工况对主体机构进行运动分析。
设各具有旋转运动的构件对x 轴的转角分别为i i θ , ( 为旋转构件的标号），相应的角速度和角加速度分别为ωi ,εi ；用解析法求出当曲柄转角θ1 从刨刀处于最右侧时起，沿逆时针方向转动每隔100 计算一组运动参数，其中包括：各杆的角位置、角速度、角加速度及刨刀的位置刀s （以最右点为零点）、速度刀v 和加速度刀a ，应用计算机在同一幅图中绘出刨刀的位移曲线、速度曲线和加速度曲线，并分析计算结果的合理性。
4）牛头刨床主体机构的受力分析杆的受力以及质量如表2所示。
已知数据其余构件的质量和转动惯量以及运动的摩擦忽略不计。
假定刨刀在空回行程不受力，在工作行程中所受的阻力为水平力，其大小见图3。
用解析法求出机构处于不同位置时应加在曲柄上的驱动力矩TN 以及各运动副的约束总反力的大小和方向。
图3 刨刀的有效阻力课程设计的主要内容包括：设计任务（包括设计条件和要求）；②传动方案的确定；③机构综合的方法和结果；④运动分析的方法和结果；⑤受力分析的方法和结果；⑥结束语；⑦主要参考文献；⑧附件（计算机程序等）。
牛头刨床的连杆机构运动分析
牛头刨床的连杆机构运动分析0 前言机构运动分析的任务是对于结构型式及尺寸参数已定的具体机构，按主动件的位置、速度和加速度来确定从动件或从动件上指定点的位置、速度和加速度。
许多机械的运动学特性和运动参数直接关系到机械工艺动作的质量，运动参数又是机械动力学分析的依据，所以机构的运动分析是机械设计过程中必不可少的重要环节。
以计算机为手段的解析方法，由于解算速度快，精确度高，程序有一定的通用性，已成为机构运动分析的主要方法。
连杆机构作为在机械制造特别是在加工机械制造中主要用作传动的机构型式，同其他型式机构特别是凸轮机构相比具有很多优点。
连杆机构采用低副连接，结构简单，易于加工、安装并能保证精度要求。
连杆机构可以将主动件的运动通过连杆传递到与执行机构或辅助机构直接或间接相连的从动件，实现间歇运动，满足给定的运动要求，完成机器的工艺操作。
牛头刨床是一种利用工作台的横向运动和纵向往复运动来去除材料的一种切削加工机床。
工作台的纵向往复运动是机床的主运动，实现工件的切削。
工作台的横向运动即是进给运动，实现对切削精度的控制。
本文中只分析纵向运动的运动特性。
牛头刨床有很多机构组成，其中实现刨头切削运动的六连杆机构是一个关键机构。
刨床工作时，通过六杆机构驱动刨刀作往复移动。
刨刀右行时，当刨刀处于工作行程时;要求刨刀的速度较低且平稳，以减小原动机的容量和提高切削质量。
当刨刀处于返回行程时，刨刀不工作，称为空行程，此时要求刨刀的速度较高以提高生产率。
由此可见，牛头刨床的纵向运动特性对机床的性能有决定性的影响。
1 牛头刨床的六连杆机构牛头刨床有很多机构组成，其中实现刨头切削运动的六杆机构是一个关键机构。
图1所示的为一牛头刨床的六连杆机构。
杆1为原动件，刨刀装在C点上。
假设已知各构件的尺寸如表1所示，原动件1以等角速度ω1=1rad/s沿着逆时针方向回转，要求分析各从动件的角位移、角速度和角加速度以及刨刀C点的位移、速度和加速度的变化情况。
牛头刨床分析报告
2、各齿轮参数
齿顶高系数ha* 模数 压力角 齿顶高系数 2.5 20 Z3 2.5 Z4 Z5 Z6 36 Z7 42 顶隙系数c* 顶隙系数 0.625 Z8 48 Z9 71
齿轮 Z1 Z2
齿数 24 36 30 48 25
4、变速箱各档齿轮啮合
档位 齿轮传动示意图 啮合关系线图 传动比
nⅢ
四、牛头刨床连杆机构分析
1、牛头刨床连杆机构运动简图
自由度计算
F = 3 n − 2 PL − PH
F=自由度 F= n=活动构件数 (n=N-1，N=构件数） PL=低副数 PH=高副数
2牛头刨床主机构工作原理
• 牛头刨床是刨削类机床中应用较广的一种，工作时滑枕带 着刨刀，作直线住复运动，因滑枕前端的刀架形似牛头而 得名。它适合刨削中、小型零件。滑枕的移动速度是不均 匀的。滑枕的返回行程速度大于工作行程速度。由于采用 单刃刨刀加工，且在滑枕回程时不切削，牛头刨床的生产 率较低。普通牛头刨床由滑枕带着刨刀作水平直线住复运 动，刀架可在垂直面内回转一个角度，并可手动进给，工 作台带着工件作间歇的横向或垂直进给运动，常用于加工 平面、沟槽和燕尾面等。 • 牛头刨床的主运动 主运动为：电动机→变速机构（带传动机构、 主运动 齿轮传动机构）→摇杆机构→滑枕往复运动→刨刀工作； • 牛头刨床的进给运动 进给运动为：电动机→变速机构（带传动机构、 进给运动 齿轮机构）→棘轮进给机构→工作台横向进给运动（螺旋 机构）。
三、牛头刨床工作台快速进给机 构传动简图
工作台快速进给机构传动分析
• 啮合关系线图
1—2=3—4=5—6=8
传动比计算
i=z2z4z6/z1z3z5=35×17× × × 23/2×34×23=8.75 × × 螺杆转向分析
毕业设计(论文)-牛头刨床六杆机构运动分析
毕业设计（论文）-牛头刨床六杆机构运动分析河南理工大学本科毕业设计,论文,摘要在工程技术领域，经常会遇到一些需要反复操作，重复性很高的工作，如果能有一个供反复操作且操作简单的专用工具，图形用户界面就是最好的选择。
如在本设计中对于牛头刨床平面六杆机构来说，为了保证结构参数与运动参数不同的牛头刨床的运动特性，即刨刀在切削过程中接近于等速运动从而保证加工质量和延长刀具寿命，以及刀具的急回性能从而提高生产率，这样的问题如果能够通过设计一个模型平台，之后只需改变参量就可以解决预期的问题，这将大大的提高设计效率。
设计本设计中正是通过建立牛头刨床六杆机构的数学模型，然后用MATLAB程序出一个友好的人机交互的图形界面，并将数学模型参数化，使用户只需改变牛头刨床的参数就可以方便的实现运动分析和运动仿真，用户可以形象直观地观察到牛头刨床的运动轨迹、速度变化及加速度变化规律。
关键词:牛头刨床六杆机构 MATLAB 运动仿真程序开发1河南理工大学本科毕业设计,论文,AbstractIn the engineering area, often repeatedly encountered some operational needs, repetitive highly, and if the operation can be repeated for a simple operation and dedicated tool graphical user interface is the best choice. As in the planer graphic design for six pole bodies, and campaigns to ensure the structural parameters of different parameters planer movement characteristics, planning tool inthe process of cutting close to equal campaign to ensure processing quality and extended life cutlery and cutlery rush back to the performance enhancing productivity, If such issues can be adopted to design a model platform parameter can be changed only after the expected settlement, which will greatly enhance the efficiency of the design. It is through the establishment of this design planer six pole bodies mathematical model, and then use MATLAB to devise procedures of a friendly aircraft in the world graphics interface, and mathematical models of the parameters, so that users only need to change the parameters planer can facilitate the realization of movement analysis and sports simulation, Users can visual image observed in planer movement trajectories, speed changes and acceleration changes.Keywords:Planer 6 pole bodies MATLAB Campaign simulation Procedure development.2河南理工大学本科毕业设计,论文,目录1 绪论 (4)2牛头刨床六杆机构运动分析程序设计2.1 MATLAB介绍 (5)2.2 MATLAB的特点 (6)2.3 用MATLAB处理工程问题优缺点................................7 3牛头刨床运动分析的模型3.1 基本概念与原理 (9)3.2 牛头刨床六杆机构的数学模型 .................................9 4图形用户界面GUI4.1界面设计的原则 (13)4.2 功能要求 (16)4.3界面结构设计 (17)4.4 程序框图的设计 .............................................19 5运动仿真程序界面设计与编程实现5.1 句柄图形体系 (21)5.1.1 图形对象、对象句柄和句柄图形树结构 (22)5.1.2 对象属性 (23)5.1.3 对象句柄的获取方法 (23)5.1.4 对象句柄的获取和设置 (25)5.2 主界面参数含义 (27)5.3 界面制作步骤 (27)6总结 (49)7致谢 (50)8参考书目 (51)9附录程序源代码 (52)3河南理工大学本科毕业设计,论文,1 绪论1.1本课题的意义机构运动分析是不考虑引起机构运动的外力的影响，而仅从几何角度出发，根据已知的原动件的运动规律(通常假设为匀速运动)，确定机构其它构件上各点的位移、速度、加速度，或构件的角位移、角速度、角加速度等运动参数。
牛头刨床机构运动分析
目录一、概述1.1、课程设计的目的——————————————— 21.2、工作原理—————————————————— 21.3、设计要求—————————————————— 31.4、设计数据—————————————————— 41.5、创新设计内容及工作量———————————— 4二、牛头刨床主传动机构的结构设计与分析2.1、方案分析—————————————————— 52.2、主传动机构尺寸的综合与确定————————— 52.2、杆组拆分—————————————————— 62.4、绘制刀头位移曲线图————————————— 7三、牛头刨床主传动机构的运动分析及程序3.1、解析法进行运动分析————————————— 83.2、程序编写过程（计算机C语言程序）—————— 103.3、计算数据结果——————————————— 123.4、位移、速度和加速度运动曲线图与分析————— 13四、小结心得体会——————————————————— 18五、参考文献参考文献——————————————————— 19一、概述1.1、课程设计的目的目的：机械课程创新设计是培养学生机械系统方案设计能力的技术基础课程，他是机制专业课程学习过程中的一个重要实践环节。
其目的是以机制专业课程的学习为基础，进一步巩固和加深所学的基本理论、基本概念和基本知识，培养学生分析和解决与本专业课程有关的具体机械所涉及的实际问题的能力，使学生熟悉机械系统设计的步骤及方法，其中包括选型、运动方案的确定、运动学和动力学的分析和整体设计等，并进一步提高计算、分析、计算机辅助设计、绘图以及查阅和使用文献的综合能力。
1.2、工作原理牛头刨床是一种靠刀具的往复直线运动及工作台的间歇运动来完成工件的平面切削加工的机床。
图1为其参考示意图。
电动机经过减速传动装置（皮带和齿轮传动）带动执行机构（导杆机构和凸轮机构）完成刨刀的往复运动和间歇移动。
(完整版)牛头刨床运动分析实例
例： 如图所示为一牛头刨床的机构运动简图。
设已知各构件尺寸为：1125mm l =，3600mm l =，4150mm l =，原动件1的方位角1=0~360θ︒︒和等角速度1=1rad/s w 。
试用矩阵法求该机构中各从动件的方位角、角速度和角加速度以及E 点的位移、速度和家速度的运动线图。
解：先建立一直角坐标系，并标出各杆矢量及方位角。
其中共有四个未知量3θ、4θ、3s 及E s 。
为求解需建立两个封闭矢量方程，为此需利用两个封闭图形ABCA 及CDEGC ，由此可得，613346,'E l l s l l l s +=+=+(1-1)写成投影方程为： 331133611334433446cos cos sin sin cos cos 0sin sin 'E s l s l l l l s l l l θθθθθθθθ==++-=+= (1-2)解上面方程组，即可求得3θ、4θ、3s 及E s 四个位置参数，其中23θθ=。
将上列各式对时间取一次、二次导数，并写成矩阵形式，即可得以下速度和加速度方程式。
速度方程式：3331133331131334443344cos sin 00sin sin cos 00cos 0sin sin 1000cos cos 0E s l s s l w w l l w l l v θθθθθθθθθθ⎡⎤--⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥---⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦ (1-3)机构从动件的位置参数矩阵：33333333443344cos sin 00sin cos 000sin sin 10cos cos 0s s l l l l θθθθθθθθ-⎡⎤⎢⎥⎢⎥⎢⎥---⎢⎥⎣⎦机构从动件的的速度列阵：334E s w w v ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦机构原动件的位置参数矩阵：1111sin cos 00l l θθ-⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦1w ：机构原动件的角速度加速度方程式：333333333344433443333333333333333344433344cos sin 00sin cos 000sin sin 10cos cos 0sin sin cos 00cos cos sin 000cos cos 00sin sin E s s s l l l l w s s w w s s w l w l w l w l w θθθθαθθαθθαθθθθθθθθθθ⎡⎤-⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥---⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦----=-----11131113144cos sin 000E l w s l w w w w v θθ⎡⎤⎡⎤-⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥+⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎣⎦(1-4)机构从动件的位置参数矩阵求导：33333333333333333444333444sin sin cos 00cos cos sin 000cos cos 00sin sin 0w s s w w s s w l w l w l w l w θθθθθθθθθθ⎡⎤---⎢⎥⎢⎥-⎢⎥⎢⎥--⎢⎥--⎢⎥⎣⎦机构从动件的的加速度列阵：334E s ααα⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦机构原动件的位置参数矩阵求导：111111cossinl wl wθθ-⎡⎤⎢⎥-⎢⎥⎢⎥⎢⎥⎣⎦主程序（matlab）：%牛头刨床运动分析主程序s;%x(1)——代表3θ；%x(2)——代表构件3的转角3θ；%x(3）——代表构件4的转角4s;%x(4)——代表E点的线位移El；%x(5)——代表1l;%x(6)——代表3l;%x(7)——代表4l;%x(8)——代表6l;%x(9）——代表'6w。
牛头刨床主传动机构运动方案分析
三、机构选型、方案分析及方案的确定方案一的运动分析及评价（1）运动是否具有确定的运动该机构中构件n=5。
在各个构件构成的的运动副中Pl=6,Ph=1.凸轮和转子、2杆组成运动副中有一个局部自由度，即F'=1。
机构中不存在虚约束。
.由以上条件可知：机构的自由度 F=3n-(2Pl+Ph-p')-F'=1机构的原动件是凸轮机构，原动件的个数等于机构的自由度，所以机构具有确定的运动。
（2）机构传动功能的实现在原动件凸轮1带动杆2会在一定的角度范围内摇动。
通过连杆3推动滑块4运动，从而实现滑块(刨刀)的往复运动。
（3）主传动机构的工作性能凸轮1的角速度恒定，推动2杆摇摆，在凸轮1随着角速度转动时，连杆3也随着杆2的摇动不断的改变角度，使滑块4的速度变化减缓，即滑块4的速度变化在切削时不是很快，速度趋于匀速；在凸轮的回程时，只有惯性力和摩擦力，两者的作用都比较小，因此，机构在传动时可以实现刨头的工作行程速度较低，而返程的速度较高的急回运动。
传动过程中会出现最小传动角的位置，设计过程中应注意增大基圆半径，以增大最小传动角。
机构中存在高副的传动，降低了传动的稳定性。
（4）机构的传力性能要实现机构的往返运动，必须在凸轮1和转子间增加一个力，使其在回转时能够顺利的返回，方法可以是几何封闭或者是力封闭。
几何封闭为在凸轮和转子设计成齿轮形状，如共扼齿轮，这样就可以实现其自由的返回。
机构在连杆的作用下可以有效的将凸轮1的作用力作用于滑块4。
但是在切削过程中连杆3和杆2也受到滑块4的作用反力。
杆2回受到弯力，因此对于杆2的弯曲强度有较高的要求。
同时，转子与凸轮1的运动副为高副，受到的压强较大。
所以该机构不适于承受较大的载荷，只使用于切削一些硬度不高的高的小型工件。
该机构在设计上不存在影响机构运转的死角，机构在运转过程中不会因为机构本身的问题而突然停下。
（5）机构的动力性能分析。
由于凸轮的不平衡，在运转过程中，会引起整个机构的震动，会影响整个机构的寿命。
牛头刨床导杆机构的运动分析
牛头刨床导杆机构的运动分析1. 引言牛头刨床导杆机构是一种常用于刨削加工的机床设备，具有高精度和高效率的特点。
在刨床导杆机构中，导杆起着重要的作用，它能够提供刨床刀架的平行导向和支撑作用。
导杆机构的运动特性对于刨削加工的精度和效率有着重要的影响。
因此，对于牛头刨床导杆机构的运动分析具有重要意义。
2. 导杆机构的构成牛头刨床导杆机构一般由导杆、导轨、导轨座以及导杆座等部件组成。
导杆是导轨上移动的零件，通过导杆座和刨床刀架相连。
导轨是固定在刨床床身上的导向轨道，导轨座则用于支撑导杆并保证其在导轨上的平稳移动。
3. 运动分析导杆机构的运动分析主要包括导杆的运动规律和导轨座的运动约束等方面。
3.1 导杆的运动规律导杆在刨床加工过程中的运动规律可由导杆机构的几何关系来描述。
假设导杆长度为L，导轨座的前进长度为x，导杆与导轨的接触点为A，刨床刀架的前进位置为y。
根据几何关系可知，导杆的运动规律满足以下方程：x = L - y3.2 导轨座的运动约束为了保证导杆在导轨上的平稳运动，导轨座需要满足一定的运动约束条件。
一般情况下，导轨座的运动约束可由导轨座与导轨的接触点A、导杆与导轨的接触点B以及导杆长度L所构成的闭合平面链进行分析。
根据平面链的运动学原理，可以得到以下方程：L = AB + BC其中，AB代表导轨座的前进长度，BC代表导杆的伸出长度。
根据实际情况，可以确定导轨座的运动范围，并进行相应的设计。
4. 结论通过对牛头刨床导杆机构的运动分析，可以得出以下结论：1.导杆的运动规律可由导杆机构的几何关系来描述，其中导杆的前进位置与导轨座的前进长度满足特定的关系。
2.导轨座的运动约束需要满足导轨座与导轨接触点、导杆与导轨接触点以及导杆长度所构成的闭合平面链的特定关系。
综上所述，牛头刨床导杆机构的运动分析对于提高刨削加工的精度和效率具有重要意义，可以为刨床设计和制造提供理论依据和技术支持。
参考文献：[1] 张三. 牛头刨床导杆机构运动分析研究[J]. 机床与液压, 2010, 12(4): 45-50.。
牛头刨床导杆机构的设计及运动分析
机械原理课程设计编程说明书设计题目：牛头刨床导杆机构的设计及运动分析（4）目录一、计任务及要求……………………………………………………二、数学模型的建立………………………………………………三、程序框……………………………………………………四、程序中符号说明…………………………………………五、程序清单及运行结果………………………………………六、课程设计总结………………………………………………七、参考文献………………………………………………牛头刨床导杆机构运动分析一、设计任务及要求:已知：牛头刨床的导杆机构的曲柄每分钟转速n2，各构件尺寸及重心位置，且报头导路x-x位于导杆端点B所作圆弧的平分线上。
数据如下表所示：设计内容导杆机构运动分析单位r/min mm符号n1 LO2O4LO2A LO4B LBE LO4S4数据70 400 95 800 ２５６４００要求：1、用c语言编写计算程序,对机构进行动态分析和动态显示。
2、上机调试程序并打印结果。
3、画出导杆4的角位移，角速度，角加速度的曲线。
4、编写设计计算说明书。
二、 数学模型的建立：ABCXa b d Z 2Z 1Z 3βαY该牛头刨床导杆机构为六杆机构，拆分成两个四杆机构：1）摆动导杆机构；2） 曲柄滑块机构。
求导杆4的角位移，角速度，角加速度，分析摆动导杆机构。
如图所示建立坐标。
三个向量构成封闭图形，所以可得：O Z Z Z =−→−+−→−-−→−342 （1）按复数形式可以写成：上式可以简化为,900)sin (cos )sin (cos )sin (cos 333︒==+++-+θθθϕϕααi d i b i a，求解得）式对时间求二阶导数将（得）对时间求一阶导数，将（联立解得）（）（相等得）式中实部、虚部分别根据（）（。
2)8)(sin()7)(cos(.2)6(sin 2)5(cos sin )4(),3(40sin sin 30cos cos 220)sin (cos )sin (cos )sin (cos 2242233ϕαωϕαωωαααϕϕαϕαθθϕϕααϕ--=-=++=+==+-=-=+++-+=a b v ad d a b a a d arctgd b a b a i d i b i a A ba)9](2)sin()cos([142224ωϕαωϕαξϕξA v a a b----==。
牛头刨床横向进给运动分析.
横向进给机构及进给量的调整
1、2—齿轮 3—连杆 4—棘爪 5—棘轮 6—丝杆 7—棘轮护盖
B6065牛头刨床横向进给运动及调整 由变速机构将运动传递给齿轮2，齿轮2带动齿轮1（曲柄）转动， 通过连杆3驱使棘爪4（摇杆）往复摆动。正向进给时，棘爪4拨动 棘轮5使丝杆6转一个角度，实现横向进给。反向时，由于棘爪后面 是斜的，爪内弹簧被压缩，棘爪从棘轮顶滑过。因此工作台横向自 动进给是间歇的。
进给量分析例题：
已知牛头刨床工作台横向进给单线螺纹的螺距 P为6mm，最小进给量Lmin为0.2mm，最大横向 进给量Lmax为2.4mm，求该棘轮齿数为Z多少？ 当进给量为1mm时，棘轮应转过多少个齿？进给 量最大时转过多少齿，转过多少度？
解： 由棘轮和单线螺纹的相关知识可知：
P=Lmin×Z
牛头刨床横向进给量分析
牛头刨床横向进给量分析
导入
通过学习牛头刨床的知识 我们知道牛头刨床的切削运动， 普通牛头刨床由滑枕带着刨刀 作水平直线往复运动，刀架可 在垂直面内回转一个角度，并 可手动进给，工作台带着工件 作间歇的横向或者垂直运动， 常用于加工平面、沟槽和燕尾 面等。
牛头刨床的进给运动为： 电动机——变速机构——齿轮 进给机构——工作台横向进给 机构（棘轮机构和螺旋机构）
则Z=P/Lmin=6÷0.2=30
当进给量为1mm时其转过的齿当最大进给量时转过的齿数为 Lmax÷Lmin=2.4÷0.2=12；
当最大进给量时转过的角度为 360°×12/30=144°
横向进给机构及进给量的调整
棘轮机构：将棘爪的往复摆动转化为棘轮的间歇转动 螺旋机构：将棘轮（丝杆）的间歇转动转化为工作台（螺母） 的间歇直线移动 结论：工作台的每次进给量取决于棘轮每次转过的角度，即棘 爪每次拔动棘轮转过的齿数。
牛头刨床机构的运动特点
牛头刨床机构的运动特点
牛头刨床机构的运动特点是指它在工作时产生的机构运动特征。
牛头刨床机构是一种将直线运动转化为往复运动的机械结构，其主要特点包括：
1. 往复运动：牛头刨床机构通过曲轴和连杆组成，使得工作台面能够以往复运动的方式进行刨削操作。
这种往复运动是机构的基本特点，能够有效地实现工作过程中的刨削加工。
2. 快速往复：牛头刨床机构的刨床行程通常较长，且具有快速往复的特点。
它能够在短时间内完成一次完整的刨削过程，提高生产效率。
3. 平稳可靠：牛头刨床机构的运动稳定性较好，能够在工作过程中保持平稳的运动状态，避免产生过大的振动和冲击力，保证加工质量和切削工具的寿命。
4. 节能高效：牛头刨床机构在运动时能够利用较少的能量来完成工作，实现能耗的节约。
同时，其运动速度较快，能够提高生产效率，减少加工时间。
5. 方向可逆：牛头刨床机构的运动方向是可逆的，可以根据需要进行前后或上下刮削，具有较好的灵活性和适应性。
总的来说，牛头刨床机构的运动特点是往复、快速、平稳可靠、节能高效和方向可逆。
这些特点使得牛头刨床机构成为一种常用的刨削设备，广泛应用于金属加工和木工加工等领域。
牛头刨床机构的运动特点
牛头刨床机构的运动特点
牛头刨床是一种常用的木工机床，用于加工木料的平整表面。
其机构的运动特点主要包括刨刀的直线来回运动、工作台的往复运动以及刨床本身的旋转运动。
牛头刨床机构的运动特点之一是刨刀的直线来回运动。
刨刀是刨床上最重要的部件之一，它负责将木料上的凸起部分切削掉，使其表面变得平整。
刨刀的运动是通过传动装置实现的，通常是由电机驱动传动带或传动链带动刨刀轴旋转，使刨刀进行直线来回运动。
这种运动方式可以保证刨刀在切削过程中保持稳定的切削速度和切削力，从而得到平整的切削表面。
牛头刨床机构的运动特点之二是工作台的往复运动。
工作台是刨床上承载木料的部件，它的运动方式通常是沿着刨刀运动方向往复移动。
工作台的往复运动可以通过传动装置实现，通常是由电机驱动传动带或传动链带动工作台进行往复运动。
工作台的往复运动可以使木料在刨刀的作用下平稳地进行切削，保证切削表面的平整度和精度。
牛头刨床机构的运动特点之三是刨床本身的旋转运动。
刨床通常由机身和床身组成，机身上安装有刨刀和传动装置，床身上安装有工作台和支撑装置。
刨床本身可以通过传动装置实现旋转运动，通常是由电机驱动传动带或传动链带动刨床旋转。
刨床的旋转运动可以
使木料在刨刀的作用下均匀地进行切削，保证切削表面的平整度和精度。
牛头刨床机构的运动特点主要包括刨刀的直线来回运动、工作台的往复运动以及刨床本身的旋转运动。
这些特点使得牛头刨床能够有效地进行木料的切削加工，保证切削表面的平整度和精度。
在实际应用中，根据不同的加工要求和木料特性，可以调整刨刀的运动速度、工作台的运动频率和刨床的旋转速度，以获得最佳的切削效果。
机械原理牛头刨床课程设计--牛头刨床导杆机构的运动分析、动态静力分析
青岛理工大学琴岛学院课程设计说明书课题名称：机械原理课程设计学院：机电工程系专业班级：机械113学号：20110201083学生：张三指导老师：李燕青岛理工大学教务处2013 年 12月 27日《机械原理课程设计》评阅书摘要选取方案三，利用图解法对1点和6电状态时牛头刨床导杆机构进行运动分析、动态静力分析，并汇总本方案所得各位置点的速度、加速度、机构受力数据绘制曲线图。
进行方案比较，确定最佳方案。
将一个班级分为 3 组，每组11人左右，一组选择一个备选方案进行如下分析工作：课程设计内容：牛头刨床导杆机构的运动分析、动态静力分析；（1）绘制机构运动简图（两个位置）；（2）速度分析、加速度分析；（3）机构受力分析（求平衡力矩）；（4）绘制运动线图。
（上述三项作在一张A1号图纸上）精选文档目录摘要 (I)1设计任务 (1)2导杆机构的基本尺寸确定 (2)3 导杆机构的运动分析 (4)3.1 速度分析 (4)3.2 加速度分析 (5)4导杆机构的动态静力分析 (8)4.1 运动副反作用力分析 (8)4.2 曲柄平衡力矩分析 (10)总结 (11)致谢 (12)参考文献 (13)1设计任务一、课程设计的性质、目的和任务机械原理课程设计是高等工业学校机械类专业学生第一次较全面的机械运动学和动力学分析与设计的训练，是本课程的一个重要教学环节。
其意义和目的在于：以机械系统运动方案设计为结合点，把机械原理课程设计的各章理论和方法融会贯通起来，进一步巩固和加深学生所学的理论知识；培养学生独立解决有关本课程实际问题的能力，使学生对于机械运动学和动力学的分析和设计有一个较完整的概念，具备计算、制图和使用技术资料的能力。
二、课程设计教学的内容和要求将一个班级进行分组，每组10人左右，一组选择一个备选方案进行如下分析工作：课程设计内容：牛头刨床导杆机构的运动分析、动态静力分析；（1）绘制机构运动简图；（2）速度分析、加速度分析；（1张1号图纸）（3）机构动态静力分析；（4）绘制运动线图。
牛头刨床运动分析实例
牛头刨床运动分析实例
牛头刨床是一种钻床，其切削过程是将刀具在工件上沿着轴向切削，切削时刀具的进给和转速都可以调整。
牛头刨床的切削过程主要涉及到以下几个运动：
1. 主轴转动运动：主轴是牛头刨床的重要组成部分，它承担着刀具的转动任务。
主轴转动时，刀具也会随之转动，将工件上的材料切削掉。
2. 工件进给运动：工件的进给运动与主轴转动相反，它将工件沿着切削方向向刀具推进。
工件进给速度可以控制，不同材料需要不同的进给速度。
3. 横向移动运动：牛头刨床的刀具可以在横向方向上移动，这个运动通常用于切削不规则形状的工件。
切削刀具可以通过手动操作或数控系统控制横向移动。
以上三个运动相互协调，在不同的切削任务中产生不同的切削效果。
牛头刨床是一种精密的切削工具，通常用于制造高精度零件，例如轴承、齿轮等。
机械原理课程设计说明书-牛头刨床的运动分析与设计
机械原理课程设计说明书-牛头刨床的运动分析与设计一、设计目标本机械原理课程设计的目标是对牛头刨床进行运动分析与设计，通过分析刨床的运动原理和结构特点，设计出合理的刨床结构，确保刨床的运动稳定性和工作效率。
二、刨床的运动分析1. 刨床的基本运动牛头刨床的基本运动包括主轴转动、工作台进给运动和刀架进给运动。
主轴转动通过电动机驱动刨刀进行旋转，实现刨削工作。
工作台进给运动使工件在水平平面上进行进给运动，供刀架进行刨削。
刀架进给运动使刀架在垂直于工作台的方向上进行进给，并在工件刨削时左右平移，调整刨削的位置。
2. 刨床的运动传动刨床的运动传动主要通过齿轮传动和导轨传动实现。
主轴转动通过电动机通过齿轮传动带动主轴实现。
工作台进给运动通过齿轮和导轨的组合实现，工作台在导轨上进行水平移动。
刀架进给运动通过螺杆和导轨的组合实现，螺杆带动刀架进行垂直平移，并在导轨上进行水平移动。
三、刨床结构设计基于上述运动分析，对牛头刨床进行结构设计如下：1. 主轴结构：主轴采用直径大、刚度高的优质轴承，保证刨床的稳定性和工作效率。
主轴和电动机通过齿轮传动连接，确保刨床主轴的转动平稳。
2. 工作台结构：工作台采用结实的铸铁材料，设计为可拆卸结构，方便工件的放置和取出。
工作台通过导轨和齿轮传动实现水平进给运动，导轨和齿轮选用耐磨材料，减小运动阻力。
3. 刀架结构：刀架采用铸铁材料，设计为可调节结构，方便调整刨削位置。
刀架通过螺杆和导轨的组合实现垂直进给运动和水平进给运动，确保刀具与工件的接触面平整。
四、设计流程1. 进行刨床的运动分析，确定刨床的基本运动和运动传动方式。
2. 根据运动分析结果，进行刨床的结构设计，包括主轴结构、工作台结构和刀架结构。
3. 设计刨床各部件的尺寸和连接方式，确保结构的牢固性和可拆卸性。
4. 进行刨床的总体装配和调试，确保刨床的运动平稳和工作效率。
5. 测试刨床的性能和稳定性，进行必要的调整和改进。
五、安全注意事项1. 在使用刨床时，应仔细阅读操作指南，并按照操作规程进行操作。
牛头刨床导杆机构的运动分析、动态静力分析
摘要——牛头刨床运动和动力分析一、机构简介与设计数据1、机构简介牛头刨床是一种用于平面切削加工的机床，如图1-1a。
电动机经皮带和齿轮传动，带动曲柄2和固结在其上的凸轮8。
刨床工作时，由导杆机构2 –3 –4 –5 –6 带动刨头6和刨刀7作往复运动。
刨头右行时，刨刀进行切削，称工作行程，此时要求速度较低并且均匀，以减少电动机容量和提高切削质量；刨头左行时，刨刀不切削，称空回行程，此时要求速度较高，以提高生常率。
为此刨床采用有急回作用的导杆机构。
刨刀每切削完一次，利用空回行程的时间，凸轮8通过四杆机构1 – 9 – 10 – 11 与棘轮带动螺旋机构（图中未画），使工作台连同工件做一次进给运动，以便刨刀继续切削。
刨头在工作行程中，受到很大的切削阻力（在切削的前后各有一段约0.05H的空刀距离，图1-1b），而空回行程中则没有切削阻力。
因此刨头在整个运动循环中，受力变化是很大的，这就影响了主轴的匀速运转，故需安装飞轮来减小主轴的速度波动，以提高切削质量和减少电动机容量a b图目录摘要 (III)1设计任务 (1)2 导杆机构的运动分析 (2)导杆机构的动态静力分析 (4)3.1运动副反作用力分析 (4)3.2力矩分析 (6)4方案比较 (7)5总结 (10)6参考文献 (10)《机械原理课程设计》说明书1设计任务机械原理课程设计的任务是对机器的主题机构进行运动分析。
动态静力分析，确定曲柄平衡力矩，并对不同法案进行比较，以确定最优方案。
要求根据设计任务，绘制必要的图纸和编写说明书等。
2 导杆机构的运动分析2.1 速度分析取曲柄位置1’对其进行速度分析，因为2和3在以转动副相连，所以V A2=V A3，其大小等于ω2l02A，指向于ω2相同。
取构件3和4的重合点A进行速度分析。
列速度矢量方程，得υA4 = υA3 + υA4A3大小 ? √ ?方向⊥O4A ⊥O2A ∥O4B选比例尺μv=0.004（m/s）/mm，做出速度矢量图（见图a）νA4=0.088m/sνA3=0.816m/s取5构件作为研究对象，列速度矢量方程，得υC5 = υB5 + υC5B5大小 ? √ ?方向∥XX ⊥O4B ⊥BC取速度极点p，选比例尺μv=0.004（m/s）/mm，做出速度矢量图（见图a）νC5=0.16m/sνC5B5=0.044m/s2.2 加速度分析取曲柄位置“1”进行加速度分析。
牛头刨床的运动和动力分析
二、牛头刨床的运动、动力分析[H,L1]=solve('H=300','L1=(270*H)/(2*550)',' H','L1')H =300L1 =810/11>> [N1,W1]=solve('N1=23',' W1=(N1*2*PI)/60','N1','W1')N1 =23W1 =23/30*PI一、任务根据牛头刨床的机构简图及必要的数据，进行机构的运动学和动力学分析，并给出刨头的位移、速度、加速度和曲柄平衡力矩的曲线。
二、已知条件1、机构运动简图2、机构尺寸mm a 270=， mm b 520=， mm l 5503=， mm l 1004=3、刨头行程和曲柄转速行程mm H 420=， 转速min /891r n =4、刨头的切削阻力工作行程始终为1000N ，空程为0N 。
三、 表达式推导如图所示以A 点为坐标原点，平行刨头运动方向为x 轴，建立直角坐标系，标出各杆矢量及方位角。
由机构的结构分析有：过D 和D '做刨头所在导轨的垂线DG 和D 'G '，从图形中的角度关系易证明GE= G 'E '，所以有EE '=DD '及EE '=θsin 23l =H al l =132，因而我们可以得到312l aHl =1、推导出刨头()1ϕE E x x =，()1ϕE E v v =，()1ϕE E a a =的数学表达式。
（1） 位置分析由矢量封闭三角形ABC 可得封闭矢量方程为CB AB l l a=+即： 312ii i AB ae l e Se πϕϕ+= （1）应用欧拉公式θθθsin cos i e i +=，将（1）的虚部和实部分离得： 31sin sin ϕϕS l a AB =+ （2）31cos cos ϕϕS l AB = （3） 由上面两式求解可得： 当0cos 1=ϕ即21πϕ=或231πϕ=时，由（3）得 0cos 3=ϕ 及 23πϕ=当0cos 1≠ϕ时：113cos sin tan ϕϕϕAB AB l l a +=（4）此时按照机构结构简图及反正切的定义范围易得： 当0tan 3<ϕ时： 113c o s s i n a r c t a n ϕϕπϕAB AB l l a ++=当0tan 3>ϕ时： 113c o s s i n a r c t a nϕϕϕAB AB l l a +=由矢量封闭图形CFED 可得封闭矢量方程为D E CD FE CF l l l l+=+即： 43432ϕϕπi i iE e l e l be x +=+ （5）应用欧拉公式将（4）的虚部和实部分离得：4433cos cos ϕϕl l x E += （6） 4433sin sin ϕϕl l b += （7）由式（7）可得 4334sin sin l l b ϕϕ-=此时按照机构结构简图及反正弦函数的定义范围易得： 4334sin arcsinl l b ϕπϕ--= （8） 由这两个式子可以消去4ϕ，得到由1ϕ确定的E x 的公式：()()23323324sin cos ϕϕl b l x l E -+-= （9） 显然式（9）作为计算式时难以确定E x 的符号，因此在编程时我选择式（6）作为计算式。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
高等机构学题目: 牛头刨床机构运动分析院系名称：机械与动力学院专业班级：机械工程学生姓名：学号：学生姓名：学号：学生姓名：学号：指导教师：2015年12 月17日目录一问题描述................................................................................................................................ - 1 -二运动分析................................................................................................................................ - 1 -2.1矢量法构建机构独立位置方程 ............................................................................. - 1 -2.2机构速度分析 ............................................................................................................. - 2 -2.3机构加速度分析......................................................................................................... - 2 -2.4机构运动线图绘制.................................................................................................... - 3 -三总结......................................................................................................................................... - 4 -附录一：Matlab程序............................................................................................................... - 4 -牛头刨床机构运动分析一 问题描述如图1-1所示的牛头刨床机构中，800h mm =，1360h mm =，2120h mm =，200AB l mm =，960CD l mm =，160DE l mm =。
设曲柄以等角速度15/rad s ω=逆时针方向回转，试对其进行运动分析，求出该机构中各从动件的方位角、角速度和角加速度以及各机构的运动线图。
图1-1 牛头刨床机构二 运动分析2.1矢量法构建机构独立位置方程如图2-1所示，以E 为坐标原点建立直角坐标系，并标出各杆矢量及其方位角。
其中共有四个未知量334,,,c S S θθ。
图2-1 坐标系建立以两个封闭图形ABDEA 和EDCFE 为基准构建两个封闭矢量位置方程，即：134AE l S l +=+34c h S l l +=+将上述矢量方程分别沿X 轴和Y 轴进行投影，得牛头刨床机构的独立位置方程如下：3344211cos cos cos s l h l θθθ+=+ 3344111sin sin sin s l h l θθθ+=+3344cos cos c l l s θθ+= 3344sin sin l l h θθ+=利用Matlab 进行编程求解，可求得各机构的位置，程序见附录一。
2.2机构速度分析将机构的位置方程对时间求一次导数，并写成矩阵的形式，得机构的速度方程如下：33344311333443111334443344cos sin sin 0sin sin cos cos 0cos 0sin sin 100cos cos 00c s l s l s l l l l l l s θθθθθθθωθωθθωθθ---⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦ 利用Matlab 进行编程求解，可求得各机构的角速度或速度，程序见附录一。
2.3机构加速度分析将机构的速度方程对时间求一次导数，并写成矩阵的形式，得机构的加速度方程如下：33344333333344433344333333344433444333443344cos sin sin 0sin sin cos cos 0sin cos cos 0cos cos sin sin 00sin sin 10cos 0cos cos 0c s l s s s l s l s s l l l l l l l s θθθθθωθωθθθθαθθωθωθθθαωθωθθ------⎡⎤⎡⎤⎢⎥⎢⎥--⎢⎥⎢⎥=-⎢⎥⎢⎥----⎢⎥⎢⎥⎣⎦⎣⎦33443334441111111cos 00sin sin 0cos sin 00c s l l s l l ωθωωθωθωθωθω⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦-⎡⎤⎢⎥-⎢⎥+⎢⎥⎢⎥⎣⎦利用Matlab 进行编程求解，可求得各机构的角加速度或加速度，程序见附录一。
2.4机构运动线图绘制通过Matlab进行计算求解，得到各构件的位置、速度和加速度，如表2-1所示。
根据所求得的各构件的位置、速度及加速度，进行机构运动线图的绘制，如图2-2所示。
程序见附录一。
表2-1 各构件的位置、速度和加速度图2-2 机构的运动线图三 总结通过对牛头刨床机构的运动分析，让我们学会了如何使用矩阵法建立平面机构的运动方程。
对机构进行运动分析的关键是独立位置方程的建立和求解，由于独立位置方程是一个非线性方程组，计算难度较大。
本文借用了Matlab 软件进行编程求解独立位置方程，同时对牛头刨床机构进行了运动仿真，并绘制了牛头刨床机构的运动线图，完成了从理论分析到编程求解的运动分析过程。
附录一：Matlab 程序(1)子函数PosionFun.mfunction f=Position_Fun(x,theta1,h,h1,h2,l1,l3,l4) f=[ x(1)*cos(x(2))+l4*cos(x(3))-h2-l1*cos(theta1); x(1)*sin(x(2))+l4*sin(x(3))-h1-l1*sin(theta1); l3*cos(x(2))+l4*cos(x(3))-x(4);位移线图时间/s 角位移 / ︒位移/m角速度线图时间/s角速度 / r a d ⋅s -1速度/m ⋅s -1角加速度线图时间/s角加速度 / r a d ⋅s -2加速度/m ⋅s -2牛头刨床运动仿真mmm ml3*sin(x(2))+l4*sin(x(3))-h];end(2)子函数Six_Bar.mfunction [theta,omega,alpha]=Six_Bar(theta0,theta1,omega1,alpha1,h,h1,h2,l1,l3,l4) theta=fsolve(@(x)Position_Fun(x,theta1,h,h1,h2,l1,l3,l4),theta0);S3=theta(1);theta3=theta(2);theta4=theta(3);Sc=theta(4);%计算连杆3、连杆4、滑块2和C点的速度A=[ cos(theta3) -S3*sin(theta3) -l4*sin(theta4) 0;sin(theta3) S3*cos(theta3) l4*cos(theta4) 0;0 -l3*sin(theta3) -l4*sin(theta4) 1;0 l3*cos(theta3) l4*cos(theta4) 0];B=[-l1*sin(theta1);l1*cos(theta1);0;0];omega=A\(omega1*B);v3=omega(1);omega3=omega(2);omega4=omega(3);vc=omega(4);%计算连杆3、连杆4的角加速度，滑块2及C点的加速度A=[ cos(theta3) -S3*sin(theta3) -l4*sin(theta4) 0;sin(theta3) S3*cos(theta3) l4*cos(theta4) 0;0 -l3*sin(theta3) -l4*sin(theta4) 1;0 l3*cos(theta3) l4*cos(theta4) 0];At=[-sin(theta3) -v3*sin(theta3)-S3*omega3*cos(theta3) -l4*omega4*cos(theta4) 0;cos(theta3) v3*cos(theta3)-S3*omega3*sin(theta3) -l4*omega4*sin(theta4) 0;0 -l3*omega3*cos(theta3) -l4*omega4*cos(theta4) 0;0 -l3*omega3*sin(theta3) -l4*omega4*sin(theta4) 0];B=[-l1*sin(theta1);l1*cos(theta1);0;0];Bt=[-l1*omega1*cos(theta1);-l1*omega1*sin(theta1);0;0];alpha=A\(-At*omega+alpha1*B+omega1*Bt);a3=alpha(1);alpha3=alpha(2);alpha4=alpha(3);ac=alpha(4);end(3)主程序SixBar_main.m %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 牛头刨床机构运动分析%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%输入已知数据clear;l1=0.2;l3=0.96;l4=0.16;h=0.8;h1=0.36;h2=0.12;omega1=5;alpha1=0;hd=pi/180;du=180/pi;theta0=[0.3;60*hd;270*hd;0.45];%%%调用子函数Six_Bar 计算牛头刨床机构位移，角速度，角加速度for n1=1:459theta1(n1)=-2*pi+5.8119+(n1-1)*hd;[theta,omega,alpha]=Six_Bar(theta0,theta1(n1),omega1,alpha1,h,h1,h2,l1,l3,l4);S3(n1)=theta(1); %滑块2相对于CD杆的位移theta3(n1)=theta(2); %杆3转过的角度theta4(n1)=theta(3); %杆4转过的角度Sc(n1)=theta(4); %杆5的位移v3(n1)=omega(1); %滑块2相对于CD杆的速度omega3(n1)=omega(2); %杆3转过的角速度omega4(n1)=omega(3); %杆4转过的角速度vc(n1)=omega(4); %杆5的速度a3(n1)=alpha(1); %滑块2相对于CD杆的加速度alpha3(n1)=alpha(2); %杆3转过的角加速度alpha4(n1)=alpha(3); %杆4转过的角加速度ac(n1)=alpha(4); %杆5的加速度theta0=theta;endthetaOmegaAlpha=[theta3'*du,theta4'*du,Sc',omega3',omega4',vc',alpha3',alpha4', ac'];xlswrite('Positon_Speed_Acceleration.xls',thetaOmegaAlpha,'sheet1','b1:j459');%%% 位移，角速度，角加速度和四杆机构图形输出figure(1);n1=1:459;t=(n1-1)*2*pi/360;% 绘角位移和位移线图subplot(2,2,1);plot(t,theta3*du,'r-.','LineWidth',1.5);hold on;grid on;axis auto;[haxes,hline1,hline2]=plotyy(t,theta4*du,t,Sc);set(hline1,'LineWidth',1.5);set(hline2,'LineWidth',1.5);grid on;hold on;title('位移线图');xlabel('时间/s');axes(haxes(1));ylabel('角位移/ \circ');axes(haxes(2));ylabel('位移/m');hold on;grid on;text(2.75,-0.4,'\theta_3');text(3,0.65,'\theta_4');text(5,-0.25,'S_c');% 绘角速度及速度线图subplot(2,2,2);plot(t,omega3,'r-.','LineWidth',1.5);grid on;hold on;axis auto;[haxes,hline1,hline2]=plotyy(t,omega4,t,vc); set(hline1,'LineWidth',1.5);set(hline2,'LineWidth',1.5);grid on;hold on;title('角速度线图');xlabel('时间/s');axes(haxes(1));ylabel('角速度/ rad\cdots^{-1}');axes(haxes(2));ylabel('速度/m\cdots^{-1}');grid on;hold on;text(1.25,0.55,'\omega_3');text(4.65,2.25,'\omega_4');text(5,-2.85,'v_c');% 绘角加速度和加速度线图subplot(2,2,3);plot(t,alpha3,'r-.','LineWidth',1.5);grid on;hold on;[haxes,hline1,hline2]=plotyy(t,alpha4,t,ac);set(hline1,'LineWidth',1.5);set(hline2,'LineWidth',1.5);grid on;hold on;title('角加速度线图');xlabel('时间/s');axes(haxes(1));ylabel('角加速度/ rad\cdots^{-2}');axes(haxes(2));ylabel('加速度/m\cdots^{-2}');grid on;hold on;text(3,6.5,'\alpha_3');text(4.25,17.5,'\alpha_4');text(1.25,-4.5,'a_c');%绘制牛头刨床机构subplot(2,2,4);n1=20;x(1)=0;y(1)=0;x(2)=l4*1000*cos(theta4(n1));y(2)=l4*1000*sin(theta4(n1));x(3)=l4*1000*cos(theta4(n1))+(S3(n1)*1000-50)*cos(theta3(n1)); y(3)=l4*1000*sin(theta4(n1))+(S3(n1)*1000-50)*sin(theta3(n1)); x(4)=h2*1000;y(4)=h1*1000;x(5)=x(4)+l1*1000*cos(theta1(n1));y(5)=y(4)+l1*1000*sin(theta1(n1));x(6)=x(3)+100*cos(theta3(n1));y(6)=y(3)+100*sin(theta3(n1));x(7)=l4*1000*cos(theta4(n1))+l3*1000*cos(theta3(n1));y(7)=l4*1000*sin(theta4(n1))+l3*1000*sin(theta3(n1));x(8)=x(7)-900;y(8)=h*1000;x(9)=x(7)+600;y(9)=h*1000;x(10)=l4*1000*cos(theta4(n1))+(S3(n1)*1000-50)*cos(theta3(n1)); y(10)=l4*1000*sin(theta4(n1))+(S3(n1)*1000-50)*sin(theta3(n1)); x(11)=x(10)+25*cos(pi/2-theta3(n1));y(11)=y(10)-25*sin(pi/2-theta3(n1));x(12)=x(11)+100*cos(theta3(n1));y(12)=y(11)+100*sin(theta3(n1));x(13)=x(12)-50*cos(pi/2-theta3(n1));y(13)=y(12)+50*sin(pi/2-theta3(n1));x(14)=x(10)-25*cos(pi/2-theta3(n1));y(14)=y(10)+25*sin(pi/2-theta3(n1));x(15)=x(10);y(15)=y(10);x(16)=0;y(16)=0;x(17)=x(4);y(17)=y(4);k=1:3;plot(x(k),y(k));hold on;k=4:5;plot(x(k),y(k));hold on;k=6:9;plot(x(k),y(k));hold on;k=10:15;plot(x(k),y(k));hold on;k=16:17;plot(x(k),y(k),'-.');hold on;grid on;axis([-350 800 -250 950]);title('牛头刨床运动仿真');grid on;xlabel('mm');ylabel('mm');plot(x(1),y(1),'o');plot(x(2),y(2),'o');plot(x(4),y(4),'o');plot(x(5),y(5),'o');plot(x(7),y(7),'o');%%%牛头刨床机构运动仿真figure(2)m=moviein(20);j=0;for n1=1:5:360j=j+1;clf;x(1)=0;y(1)=0;x(2)=l4*1000*cos(theta4(n1));y(2)=l4*1000*sin(theta4(n1));x(3)=l4*1000*cos(theta4(n1))+(S3(n1)*1000-50)*cos(theta3(n1));y(3)=l4*1000*sin(theta4(n1))+(S3(n1)*1000-50)*sin(theta3(n1));x(4)=h2*1000;y(4)=h1*1000;x(5)=x(4)+l1*1000*cos(theta1(n1));y(5)=y(4)+l1*1000*sin(theta1(n1));x(6)=x(3)+100*cos(theta3(n1));y(6)=y(3)+100*sin(theta3(n1));x(7)=l4*1000*cos(theta4(n1))+l3*1000*cos(theta3(n1));y(7)=l4*1000*sin(theta4(n1))+l3*1000*sin(theta3(n1));x(8)=x(7)-900;y(8)=h*1000;x(9)=x(7)+600;y(9)=h*1000;x(10)=l4*1000*cos(theta4(n1))+(S3(n1)*1000-50)*cos(theta3(n1)); y(10)=l4*1000*sin(theta4(n1))+(S3(n1)*1000-50)*sin(theta3(n1)); x(11)=x(10)+25*cos(pi/2-theta3(n1));y(11)=y(10)-25*sin(pi/2-theta3(n1));x(12)=x(11)+100*cos(theta3(n1));y(12)=y(11)+100*sin(theta3(n1));x(13)=x(12)-50*cos(pi/2-theta3(n1));y(13)=y(12)+50*sin(pi/2-theta3(n1));x(14)=x(10)-25*cos(pi/2-theta3(n1));y(14)=y(10)+25*sin(pi/2-theta3(n1));x(15)=x(10);y(15)=y(10);x(16)=0;y(16)=0;x(17)=x(4);y(17)=y(4);k=1:3;plot(x(k),y(k));hold on;k=4:5;plot(x(k),y(k));hold on;k=6:9;plot(x(k),y(k));。