牛头刨床运动分析 机械原理剖析

合集下载

机械原理大作业——牛头刨床

机械原理大作业——牛头刨床

机械原理大作业——牛头刨床大作业,一,平面连杆机构的运动分析题号: 6班级 : 姓名 : 学号 : 同组者 :成绩 :完成时间 :目录题目、原始数据及要求 ..................................................................... .......................1 一平面连杆机构运动分析方程 ..................................................................... . (1)1.1速度计算公式 ..................................................................... .. (2)1.2加速度计算公式 ..................................................................... ..............2 二程序 ..................................................................... (3)2.1计算程序框图 ..................................................................... (3)2.2计算源程序 ..................................................................... .........................4 三 3.1 (一组数据 Lab =200mm)计算结果 (9)3.2运动线图 ..................................................................... . (10)3.3 体会 ..................................................................... .................................... 12 四 4.1(第二组数据 Lab =150mm)计算结果 . (12)4.2 运动线图 ..................................................................... .. (13)4.3 体会 ..................................................................... .................................... 15 五 5.1(第三组数据 Lab =220mm)计算结果 . (16)5.2 运动线图 ..................................................................... (17)5.3 体会 ..................................................................... ...................................... 21 六参考资料 ..................................................................... (21)题目、原始数据及要求:图所示为一牛头刨床(?级机构)。

牛头刨床运动分析

牛头刨床运动分析

连杆机构运动分析机构运动分析的解析法有多种,其中比较常用的有矢量方程解析法、矩阵法和复数矢量法等。

用解析法作机构运动分析时,首先建立机构的位置方程,然后将其对时间求一次、二次导数,可以得到机构的速度方程和加速度方程,完成运动分析。

例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其方向与相反。

牛头刨床刨刀往复运动机构的分析与设计

牛头刨床刨刀往复运动机构的分析与设计

机械工程学院机械原理课程设计说明书设计题目:牛头刨床刨刀往复运动机构的分析与设计专业:机械设计制造及其自动化班级:13级姓名:学号指导教师:侍红岩2016年 1 月 4 日目录1 设计任务 (1)1.1 设计题目 (1)1.2 工作原理及工艺动作过程 (1)1.3 原始数据及设计要求 (1)1.4 设计任务 (2)2 系统传动方案设计 (3)2.1 曲柄滑块机构与摆动导杆机构 (3)2.2 齿轮和摆动导杆机构 (4)2.3 执行机构方案的比较 (5)2.4 执行机构方案的确定 (5)3 机构运动简图及数据分析 (7)3.1 机械结构简图 (7)3.2 牛头刨床数据分析 (8)4 机构运动分解 (10)5 主机构受力分析 (11)5.1 各运动副反力 (11)5.2 曲柄机构平衡力矩 (14)参考文献 (16)1 设计任务1.1 设计题目牛头刨床刨刀往复运动机构的分析与设计。

1.2 工作原理及工艺动作过程牛头刨床是一种用于平面切削加工的机床,整个机构的运转是由原动件1带动杆2的,通过连杆3推动滑块4运动;从而实现刨刀的往复运动。

刨头右行时,刨刀进行切削,称工作行程,此时要求速度较低并且均匀;刨头左行时,刨刀不切削,称空回行程,此时要求速度较高,以提高生产率。

1.3 原始数据及设计要求图1-1已知行程比系数K=1.4原动件曲柄1转速n1=60r/min,刨刀5行程H=400mm,其它参数为,L4=220mm,L5=180mm,L6=350mm,L |3=L3/2,H1=100mm,H2尺寸应满足传动角尽可能大;故刨刀5移动导路位于D点圆弧轨迹弦高一半处;构件重量分别为G2=200N,G5=700N,质心位于S3、S5处;构件3绕质心转动惯量J S3=1.1kg.m2,回程阻力为零,其它忽略不计。

刨刀工作阻力如图1-1所示,回程阻力为零,其它条件忽略不计。

表 11.4 设计任务(1)绘制机构运动简图。

牛头刨床机械系统及工作原理

牛头刨床机械系统及工作原理
牛头刨床机械系统及工作原理
牛头刨床解析
§ 牛头刨床工作时,装有刀架的滑枕3由床身内部的摆杆 带动,沿床身顶部的导轨作直线往复运动,使刀具实现 切削过程的主运动,通过调整变速手柄5可以改变滑枕 的运动速度,行程长度则可通过滑枕行程调节柄6调节。
§ 刀具安装在刀架2前端的抬刀板上,转动刀架上方的手 轮,可使刀架沿滑枕前端的垂直导轨上下移动。刀架还 可沿水平轴偏转,用以刨削侧面和斜面。滑枕回程时, 抬刀板可将刨刀朝前上方抬起,以免刀具擦伤已加工表 面。夹具或工件则安装在工作台1上,并可沿横梁8上的 导轨作间歇的横向移动,实现切削过程的进给运动。横 梁8还可沿床身的竖直导轨上、下移动,以调整工件与 刨刀的相对位置。
V带传动特点分析
§ 当量摩擦因数较大,V带传动产生的最大摩擦力 大约是平带的3倍,允许包角小,传动比较大,中 心距较小,预紧力较小;传动效率85%~95%, 带速范围20 ~ 30m/s,传递功率可达700KW, 只适合开口传动。
PPT文档演模板
牛头刨床机械系统及工作原理
牛头刨床V带传动特点
§ 1、优点:(1)具有良好的弹性,能起吸振缓 冲作用,因而传动平稳,噪音小;(2)过载时, 带与带轮会出现打滑,防止其它零件损坏;(3) 结构简单,成本低,加工和维护方便;(4)、 适用于两轴中心距较大的传动。
PPT文档演模板
牛头刨床机械系统及工作原理
双曲柄机构实例
PPT文档演模板
牛头刨床机械系统及工作原理

(其余两杆长度之和)
时————双摇杆机构
PPT文档演模板
牛头刨床机械系统及工作原理
连杆机构的急回特性
PPT文档演模板
牛头刨床机械系统及工作原理
1.凸轮机构特点及运用:

机械原理课程设计——牛头刨床.

机械原理课程设计——牛头刨床.

一:课程设计题目、内容及其目的课题:牛头刨床内容1.对机构进行运动分析已知:曲柄每分钟转数错误!未找到引用源。

,各构件尺寸及质心位置。

作机构1~2个位置的速度多边形和加速度多边形,作滑块的运动线图,以上内容与后面动态静力分析一起画在1号图纸上。

2.对机构进行动态静力分析已知:各构件的重量G(曲柄1、滑块2、和连杆5的重量都可以忽略不计),导杆3的转动惯量错误!未找到引用源。

及切削力错误!未找到引用源。

变化规律如下图。

确定构件一个位置的各运动副反力及应加于曲柄上的平衡力矩。

3、用UG进行模拟运动仿真校核机构运动分析和动态静力分析的结果4、电动机功率的确定与型号的选择5、齿轮减速机构设计目的:1:学会机械运动见图设计的步骤和方法;2:巩固所学的理论知识,掌握机构分析与综合的基本方法;3:培养学生使用技术资料,计算作图及分析与综合能力;4:培养学生进行机械创新的能力。

二:牛头刨床简介和机构的要求1:牛头刨床简介牛头刨床是一种用于平面切削加工的机床,如图1。

电动机经皮带和齿轮传动,经过减速机构减速从而带动曲柄1。

刨床工作时,由导杆3 经过连杆4 带动刨刀5 作往复运动。

刨头左行时,刨刀进行切削,称工作行程,此时要求速度较低并且均匀,以减少电动机容量和提高切削质量,刨头右行时,刨刀不切削,称空行程,此时要求速度较高,以提高生产率。

为此刨床采用有急回作用的导杆机构。

刨刀每切削完一次,利用空回行程的时间,通过棘轮带动螺旋机构(图中未画),使工作台连同工件作一次进给运动,以便刨刀继续切削。

刨头在工作行程中,受到很大的切削阻力(在切削的前后各有一段约0.05H 的空刀距离),而空回行程中只有摩擦阻力。

因此刨头在整个运动循环中,受力变化是很大的,这就影响了主轴的匀速运转,故需安装飞轮来减小主轴的速度波动,以提高切削质量和减小电动机容量。

2:机构的要求牛头刨床的主传动的从动机构是刨头,在设计主传动机构时,要满足所设计的机构要能使牛头刨床正常的运转,同时设计的主传动机构的行程要有急回运动的特性,刨削速度尽可能为匀速运动,以及很好的动力特性。

牛头刨床设计及分析best

牛头刨床设计及分析best

机械原理课程设计实习报告专业:机械设计制造及其自动化姓名:xxx学号:2010100xxxx班级:07210x指导老师:曾小慧实习时间:2012.7.5-7.13目录一、设计任务 (3)二、牛头刨床工作原理 (3)三、原始参数 (3)四、设计步骤及要求 (4)1.导杆机构的运动综合 (4)2.用解析法作导杆机构的运动分析 (6)3.导杆机构的动态静力分析 (9)4.行星轮系设计 (15)5.变位齿轮设计 (17)五、心得体会 (19)六、附录 (20)七、参考文献 (26)一、设计任务1牛头刨床刀杆机构的运动综合、运动分析和动态静力分析;2对牛头刨床传动装置中行星轮机构、齿轮机构进行综合。

二、机构工作原理牛头刨床是一种用于平面切削加工的机床,电动机经行星轮系和齿轮Z4、Z5减速带动曲柄2转动。

刨床工作时,由导杆机构2-3-4-5-6带动刨头和刨刀作往复运动。

刨头向左时,刨刀进行切削,这个行程称工作行程,刨头受到较大的切削力。

刨头右行时,刨刀不切削,称空回行程,此时要求速度较高,以提高生产力。

三、原始参数H:刨头行程;K:行程速比系数;Fc切削阻力;m4 m5 m6分别为导杆、连杆及刨头的质量;J4、J5分别分别为导杆4及导杆5绕各自质心的转动惯量;m1、m H分别为行星减速器中心轮及齿轮4、5的模数;Z4,Z5为齿轮4及5的齿数;n1:电机转速;n2:曲柄2及齿轮5的转速;k:行星轮个数。

选取第五组数据。

导杆机构的运动分析和运动综合导杆机构的动力分析H K l O2O3l O3O4/l O3B l BF/l O3B l BS5/l BF m4m5m6Js4Js5F C单位mm mm kg kg m2kg 5600 1.8 370 0.5 0.3 0.5 22 3 52 0.9 0.015 1400行星轮设计变位齿轮n 1 n 2 K 类型 m 1 Z 4Z 5m Hα单位rpmmm mm 5 1000 8032K-H514491620四、设计步骤及要求1.导杆机构的运动综合设L O3B =L 3 L BF =L 4 L O3D =L '6 L O2A =L 1 L O3O2=L 6 L O3A =S 3 L DE =S E 1、导杆的摆角ψ K=1.8180k 51.43180-︒+ψ=⇒ψ=︒︒ψ2、导杆的长度L 33H/2H 600mm L 691.4mm sin /2=⇒==ψ 3、连杆的长度L 443L 0.3L 207.4mm =⨯=4、刨头导路中心线xx 至O3点的垂直距离L '6O3E 3L L cos 2622.9mm =⨯ψ=根据已知xx 被认为通过圆弧BB ’的绕度ME 的中点D 知O E'33O3M DM 63L L L L L L 657.2mm 2-=-=-= 5、曲柄的长度L 1616L 370mm L L sin /2160.5mm =⇒=⨯ψ=6、切削越程长度0.05H ,如图所示则切削越程长度为0.05H=0.05×600=30mm 7、机构运动简图8、计算机构的自由度 F=3×5-2×7=12.用解析法作导杆机构的运动分析建立导杆机构的运动分析数学模型:如图所示,先建立一直角坐标系,并标出各杆矢量及其方位角。

牛头刨床机械原理课程设计报告3点和6点

牛头刨床机械原理课程设计报告3点和6点

牛头刨床中导杆机构的运动分析及动态静力分析第一章机械原理课程设计的目的和任务1课程设计的目的:机械原理课程设计是高等工业学校机械类学生第一次全面的机械运动学和动力学分析与设计的训练,是本课程的一个重要教学环节。

起目的在于进一步加深学生所学的理论知识,培养学生的独立解决有关课程实际问题的能力,使学生对于机械运动学和动力学的分析和设计有一个比较完整的概念,具备计算,和使用科技资料的能力。

在次基础上,初步掌握电算程序的编制,并能使用电子计算机来解决工程技术问题。

2课程设计的任务:机械原理课程设计的任务是对机器的主题机构进行运动分析。

动态静力分析,并根据给定的机器的工作要求,在次基础上设计;或对各个机构进行运动设计。

要求根据设计任务,绘制必要的图纸,编制计算程序和编写说明书等。

第二章、机械原理课程设计的方法机械原理课程设计的方法大致可分为图解法和解析法两种。

图解法几何概念比较清晰、直观;解析法精度较高。

第三章、机械原理课程设计的基本要求1.作机构的运动简图,再作机构两个位置的速度,加速度图,列矢量运动方程;2.作机构两位置之一的动态静力分析,列力矢量方程,再作力的矢量图;3.用描点法作机构的位移,速度,加速度与时间的曲线。

第四章机械原理课程设计的已知条件1、机构简介图1表1 设计数据牛头刨床是一种用于平面切削加工的机床,如图1所示。

电动机经过皮带和齿轮传动,带动曲柄2和固结在其上的凸轮8。

刨床工作时,由导杆机构2-3-4-5-6带动刨头6和刨刀7作往复运动。

刨头右行时,刨刀不切削,称为空回行程,此时要求速度较高,以提高生产率。

为此刨床采用有急回运动的导杆机构。

刨刀每切削完一次,利用空回行程的时间,凸轮8通过四杆机构1-9-10-11与棘轮带动螺旋机构(图1中未画),使工作台连同工件作一次进给运动,以便刨刀继续切削。

刨头在工作过程中,受到很大的切削阻力(在切削的前后各有一段0.05H的空刀距离,简图1,b),而空回行程中则没有切削阻力。

牛头刨床机械原理课程设计5、12点

牛头刨床机械原理课程设计5、12点

课程设计说明书—牛头刨床1. 机构简介牛头刨床是一种用于平面切削加工的机床。

电动机经皮带和齿轮传动,带动曲柄2和固结在其上的凸轮8。

刨床工作时,由导杆机构2-3-4-5-6带动刨头6和刨刀7作往复运动。

刨头右行时,刨刀进行切削,称工作行程,此时要求速度较低并且均匀,以减少电动机容量和提高切削质量;刨头左行时,刨刀不切削,称空回行程,此时要求速度较高,以提高生产率。

为此刨床采用有急回作用的导杆机构。

刨刀每次削完一次,利用空回行程的时间,凸轮8通过四杆机构1-9-10-11与棘轮带动螺旋机构,使工作台连同工件作一次进给运动,以便刨刀继续切削。

刨头在工作行程中,受到很大的切削阻力,而空回行程中则没有切削阻力。

因此刨头在整个运动循环中,受力变化是很大的,这就影响了主轴的匀速运转,故需安装飞轮来减少主轴的速度波动,以提高切削质量和减少电动机容量。

图1-11.导杆机构的运动分析已知曲柄每分钟转数n2,各构件尺寸及重心位置,且刨头导路x-x位于导杆端点B所作圆弧高的平分线上。

要求作机构的运动简图,并作机构两个位置的速度、加速度多边形以及刨头的运动线图。

以上内容与后面动态静力分析一起画在1号图纸上。

1.1设计数据牛头刨床是一种用于平面切削加工的机床。

电动机经皮带和齿轮传动,带动曲柄2和固结在其上的凸轮8。

刨床工作时,由导杆机构2-3-4-5-6带动刨头6和刨刀7作往复运动。

刨头右行时,刨刀进行切削,称工作切削。

此时要求速度较低且均匀,以减少电动机容量和提高切削质量;刨头左行时,刨刀不切削,称空回行程,此时要求速度较高,以提高生产效率。

为此刨床采用急回作用得导杆机构。

刨刀每切削完一次,利用空回行程的时间,凸轮8通过四杆机构1-9-10-11与棘轮机构带动螺旋机构,使工作台连同工件作一次进给运动,以便刨刀继续切削。

刨头在工作行程中,受到很大的切削阻力,而空回行程中则没有切削阻力。

因此刨头在整个运动循环中,受力变化是很大的,这就影响了主轴的匀速运转,故需装飞轮来减小株洲的速度波动,以减少切削质量和电动机容量。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

机械原理大作业——10A班级:机械113姓名:姚小龙学号:201106263位置方程利用两个封闭图形ABDEA 和EDCGE ,建立两个封闭矢量方程,由此可得:⎭⎬⎫+=++=+' s l l s l l l l56431643 (1) 把(1)式分别向x 轴、y 轴投影得:⎪⎪⎭⎪⎪⎬⎫=+=++=++=+ h l l s l l l h s l l h s l 33445334411133441123344sin sin cos cos sin sin sin cos cos cos θθθθθθθθθθ(2) 在(2)式中包含3s 、5s 、3θ、4θ四个未知数,消去其中三个可得到只含4θ一个未知数 方程:[][]{}[][]sin sin sin 2sin cos cos sin sin 244111234424224244112244111=-+--+-++-+θθθθθθθθl l h l hl h l l l h l l h(3)当1θ取不同值时,用牛顿迭代法解(3)式,可以求出每个4θ的值,再根据方程组(2)可以求出其他杆件的位置参数3s 、5s 、3θ的值:⎪⎪⎭⎪⎪⎬⎫-+=+=-= 3441113334453443sin sin sin cos cos )sin arcsin(θθθθθθθl l h s l l s l l h (4)速度方程对(2)式对时间求一次导数并把结果写成矩阵的形式得:⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-----00cos sin 0cos cos 01sin sin 00cos cos sin 0sin sin cos 1111143443344334433344333θθωωωθθθθθθθθθθl l v v l l l l l s l s C e B (5)其中C v 为刨刀的水平速度,v eB 为滑块2相对于杆3的速度。

由于每个1θ对应的3s 、3θ、4θ已求出,方程组式(5)的系数矩阵均为常数,采用按列选主元的高斯消去法可求解(式5)可解得角速度ω3、ω4、eB v 、C v加速度方程把(5)对时间求导得矩阵式:⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--+⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡----------=⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-----00sin cos 0sin sin 00cos cos 00sin sin cos cos 0cos cos sin sin 0cos cos 01sin sin 00cos cos sin 0sin sin cos 11111114344433344433344433333344433333343443344334433344333θωθωωωωθωθωθωθωθωθωθθωθωθωθθωααl l v v l l l l l s v l s v a a θl θl θl θl θl θs θθl θs θC eBe B eBC e B(6)同样采用按列选主元的高斯消去法可求解(6)可得角加速度3α、4α、eB a 、C a %主程序开始clear;clc;l1=180; %L1=lab l3=960; %l3=lCD l4=160; %l4=lED h=900; h1=460; h2=110; du=180/pi; omega1=1; alpha1=0;theta1=linspace(0,35*pi/18,36);%定义常量和已知参数,l1代表杆1的长度,l3代表杆3的长度,l4代表杆4的长度,h 表示EG 的长度,h1表示AE 的竖直距离,h2表示AE 的水平距离,theta1表示角θ1的不同值。

theta3=zeros(1,36);theta4=zeros(1,36);s3=zeros(1,36);s5=zeros(1,36);test=zeros(1,36);vBe=zeros(1,36);vc=zeros(1,36);omega1=ones(1,36);omega3=zeros(1,36);omega4=zeros(1,36);aBe=zeros(1,36);ac=zeros(1,36);alpha1=zeros(1,36);alpha3=zeros(1,36);alpha4=zeros(1,36);A=zeros(4,4); dA=zeros(4,1);%定义最终的结果数据,当θ1取不同值时,theta3表示θ3的值,theta4表示θ4的值,s3表示BD 的长度,s5表示GC 的长度,vBe 表示B 点在杆3上运动的速度,vc 表示杆5的运动速度,即牛头刨刀的速度,omega3表示杆3的转动角速度,omega4表示杆4的转动角速度,aBe 表示B 点在杆3上运动的角加速度,ac 表示杆5的加速度,即牛头刨刀的加速度,alpha3表示杆3的角加速度,alpha4表示杆4的角加速度,矩阵A,dA 表示线性方程组的系数矩阵i=0; %i 为循环变量,在循环结构中使用 syms THETA1 THETA4 %定义符号变量,为以下计算做准备 fun1=((h1+l1*sin(THETA1)-l4*sin(THETA4))^2+(h2+l1*cos(THETA1)-l4*cos(THETA4))^2)*(l4^2*sin(THETA4)^2+h^2-2*h*l4*sin(THETA4))-l3^2*(h1+l1*sin(THETA1)-l4*sin(THETA4))^2; %定义迭代法中要求解的关于THETA4的方程。

x0=0; %定义在牛顿迭代法中的变量THEA4的初值。

for i=1:36 %用循环结构求当theta1取不同值时,theta3值。

fun2=subs(fun1,THETA1,theta1(i));%把不同的THETA1的值代入要求解的方程[theta4(i),EA,it]=NEWTON(fun2,'THETA4',x0,0.0001,1000);%用牛顿迭代法求得THEATA4,并赋值到theta4的数组中x0=theta4(i); %把这次计算的解作为下一次计算的初值。

endfor i=1:36%用循环结构求当theta1的值取不同值时,theta3、s3、s5的取值。

因为theta3的值可能的取值范围为[0,π],对theta3求解时应分以下两种情况讨论if sign(h2+l1*cos(theta1(i))-l4*cos(theta4(i)))>0 %theta3<π/2theta3(i)=asin((h-l4*sin(theta4(i)))/l3);else theta3(i)=pi-asin((h-l4*sin(theta4(i)))/l3); %theta3>π/2endtest(i)=h1+l1*sin(theta1(i))-l4*sin(theta4(i));s5(i)=l4*cos(theta4(i))+l3*cos(theta3(i));s3(i)=(h1+l1*sin(theta1(i))-l4*sin(theta4(i)))/sin(theta3(i));endfor i=1:36%用循环结构求当theta1的值取不同值时vBe、omega3、omega4、vc的值。

A(1,1)=cos(theta3(i));A(1,2)=-s3(i)*sin(theta3(i));A(1,3)=-l4*sin(theta4(i));A(2,1)=sin(theta3(i));A(2,2)=s3(i)*cos(theta3(i));A(2,3)=l4*cos(theta4(i));A(3,2)=-l3*sin(theta3(i));A(3,3)=-l4*sin(theta4(i));A(3,4)=-1;A(4,2)=l3*cos(theta3(i));A(4,3)=l4*cos(theta4(i));dA(1,1)=-omega1(1,1)*l1*sin(theta1(i));dA(2,1)=omega1(1,2)*l1*cos(theta1(i));x=gauss(A,dA); %用按列选主元的高斯消去法求解vBe(i)=x(1);omega3(i)=x(2);omega4(i)=x(3);vc(i)=x(4);%把求得的结构赋值给各物理量endfor i=1:36 %用循环结构求当theta1的值取不同值时aBe、alpha3、alpha4、vc的值。

A(1,1)=cos(theta3(i));A(1,2)=-s3(i)*sin(theta3(i));A(1,3)=-l4*sin(theta4(i));A(2,1)=sin(theta3(i));A(2,2)=s3(i)*cos(theta3(i));A(2,3)=l4*cos(theta4(i));A(3,2)=-l3*sin(theta3(i));A(3,3)=-l4*sin(theta4(i));A(3,4)=-1;A(4,2)=l3*cos(theta3(i));A(4,3)=l4*cos(theta4(i));dA(1,1)=-omega3(i)*sin(theta3(i))*vBe(i)*2-s3(i)*omega3(i)^2*cos(theta3(i))-l4*omega4(i)^ 2*cos(theta4(i))-l1*cos(theta1(i));dA(2,1)=omega3(i)*cos(theta3(i))*vBe(i)*2-s3(i)*omega3(i)^2*sin(theta3(i))-l4*omega4(i)^2*sin( theta4(i))-l1*sin(theta1(i));dA(3,1)=-l3*omega3(i)^2*cos(theta3(i))-l4*omega4(i)^2*cos(theta4(i));dA(4,1)=-l3*omega3(i)^2*sin(theta3(i))-l4*omega4(i)^2*sin(theta4(i));%构造速度方程的系数矩阵x=gauss(A,dA); %用按列选主元的高斯消去法求解aBe(i)=x(1);alpha3(i)=x(2);alpha4(i)=x(3);ac(i)=x(4);%把求得的结构赋值给各物理量end%主程序结束%出图程序figure(1);i=1:10:360;%l3角位移图subplot(2,2,1);plot(i,theta3*du,'r');title('角位移图');xlabel('曲柄转角\theta_1/\circ');ylabel('角位移/\circ');grid on;hold on;text(200,110,'\theta3');%l4角位移图subplot(2,2,2);plot(i,theta4*du,'r');title('角位移图');xlabel('曲柄转角\theta_1/\circ');ylabel('角位移/\circ');grid on;hold on;text(150,10,'\theta4');%滑块2位移subplot(2,2,3);plot(i,s3,'r');title('位置');xlabel('滑块位置\s3/\circ');ylabel('毫米/\circ');grid on;hold on;text(150,500,'s3');%c点位移subplot(2,2,4);plot(i,s5,'r');title('位置');xlabel('滑块位置\s3/\circ');ylabel('毫米/\circ');grid on;hold on;text(150,0,'s5');figure(2);%l3角速度subplot(2,2,1);plot(i,omega3,'r');title('角速度图');xlabel('曲柄转角\theta_1/\circ') ylabel('角速度/rad\cdots^{-1}') grid on;hold on;text(150,0,'\omega_3');%l4角速度subplot(2,2,2);plot(i,omega4,'r');title('角速度图');xlabel('曲柄转角\theta_1/\circ') ylabel('角速度/rad\cdots^{-1}') grid on;hold on;text(150,0.2,'\omega_4');%c点速度subplot(2,2,3);plot(i,vc,'r');title('速度图');xlabel('曲柄转角\theta_1/\circ'); ylabel('速度mm/s');grid on;hold on;text(200,0,'Vc');%l4的速度subplot(2,2,4);plot(i,vBe,'r');title('速度图');xlabel('曲柄转角\theta_1/\circ'); ylabel('速度mm/s');grid on;hold on;text(200,-100,'Vbe');figure(3);%l3 角加速度图subplot(2,2,1);plot(i,alpha3,'r');title('角加速度图');xlabel('曲柄转角\theta_1/\circ')ylabel('角加速度/rad\cdots^{-2}')grid on;hold on;text(200,-0.1,'\alpha_3');%l4 角加速度图subplot(2,2,2);plot(i,alpha4,'r');title('角加速度图');xlabel('曲柄转角\theta_1/\circ')ylabel('角加速度/rad\cdots^{-2}')grid on;hold on;text(200,-0.5,'\alpha_4');%c点加速度图subplot(2,2,3);plot(i,ac,'r');title('加速度图');xlabel('曲柄转角\theta_1/\circ')ylabel('加速度/m\cdots^{-2}')grid on;hold on;text(200,150,'ac');%c点加速度图subplot(2,2,4);plot(i,aBe,'r');title('加速度图');xlabel('曲柄转角\theta_1/\circ')ylabel('加速度/m\cdots^{-2}')grid on;hold on;text(200,100,'Abe');%牛顿迭代法的函数定义function [r,ea,iter]=NEWTON(fun,x,x0,es,maxit)%定义函数名和输入输出的参数。

相关文档
最新文档