牛头刨床的运动和动力分析
- 格式:doc
- 大小:277.50 KB
- 文档页数:15
二、牛头刨床的运动、动力分析[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.送料系统。
刨削后的木材会被送料链条送出牛头刨床,完成整个刨削过程。
在刨削过程中,送料系统需要确保木材能够顺利送出,同时也需要对木材进行支撑,防止出现倾斜和晃动的情况,保证刨削后的木材表面平整。
总结。
牛头刨床通过送料系统将木材固定和送入刨削系统,刨削系统对木材进行刨削
操作,最后通过送料系统将刨削后的木材送出。
这就是牛头刨床的工作原理。
牛头刨床在木材加工中起着重要作用,能够有效提高木材的加工效率和加工质量,广泛应用于家具制造、建筑装饰等行业。
希望本文对牛头刨床的工作原理有所帮助,感谢阅读!。
一、概述1.课程设计的题目此次课程设计的题目是:牛头刨床的主传动结构的设计. 2.课程设计的任务和目的1)任务:1 牛头刨床的机构选型、运动方案的确定;2 导杆机构进行运动分析;3 导杆机构进行动态静力分析;4.飞轮设计;5.凸轮机构设计。
2)目的:机械原理课程设计是培养学生掌握机械系统运动方案设计能力的技术基础课程,它是机械原理课程学习过程中的一个重要实践环节。
其目的是以机械原理课程的学习为基础,进一步巩固和加深所学的基本理论、基本概念和基本知识,培养学生分析和解决与本课程有关的具体机械所涉及的实际问题的能力,使学生熟悉机械系统设计的步骤及方法,其中包括选型、运动方案的确定、运动学和动力学的分析和整体设计等,并进一步提高计算、分析,计算机辅助设计、绘图以及查阅和使用文献的综合能力。
.3.课程设计的要求牛头刨床的主传动的从动机构是刨头,在设计主传动机构时,要满足所设计的机构要能使牛头刨床正常的运转,同时设计的主传动机构的行程要有急回运动的特性,以及很好的动力特性。
尽量是设计的结构简单,实用,能很好的实现传动功能。
二.机构简介与设计数据1,机构简介牛头刨床是一种用于平面切削加工的机床,如图4-1。
电动机经皮带和齿轮传动,带动曲柄2和固结在其上的凸轮8。
刨床工作时,由导杆机构2-3-4-5-6带动刨头6和刨刀7作往复运动。
刨头右行时,刨刀进行切削,称工作行程,此时要求速度较低并且均匀,以减少电动机容量和提高切削质量;刨头左行时,刨刀不切削,称空回行程,此时要求速度较高,以提高生产效率。
为此刨床采用有急回作用的导杆机构。
刨刀每切削完一次,利用空回行程的时间,凸轮8通过四杆机构1-9-10-11与棘轮带动螺旋机构(图中没有画出),使工作台连同工件一次进级运动,以便刨刀继续切削。
刨头在工作行程过程中,受到很大的切削阻力(在切削的前后各有一段约0.05H的空刀距离,见图4-1,b)而空回行程中则没有切削阻力。
目录一、实验目的 (1)二、实验任务 (1)三、实验原理 (1)四、实验设备 (2)五、实验内容 (3)1、牛头刨床机构 (3)2、内燃机机构 (4)3、精压机机构 (4)4、两齿轮-曲柄摇杆机构 (5)5、两齿轮-曲柄摆块机构 (6)6、喷气织机开口机构 (6)7、双滑块机构 (7)8、冲压机构 (8)9、插床机构 (8)10、筛料机构 (9)11、凸轮-连杆组合机构 (9)12、凸轮-五连杆机构 (10)13、行程放大机构 (11)14、冲压机构 (11)15、双摆杆摆角放大机构 (12)六、实验方法与步骤 (12)一、实验目的1、机构创新设计与运动分析实验是综合性实验,在掌握机构组成原理、基本机构的类型、结构、设计知识的基础上,以ZBS-C机构创意设计实验台作为操作平台,进行机构创新设计实验;2、培养学生运用创新思维方法,遵循机械设计的基本法则,对机构运动系统方案进行设计与研究。
以期通过实验使学生创新意识、综合设计能力得到加强,实验技能得到提高。
二、实验任务1、选用工程机械中各种平面机构运动简图,在ZBS-C机构创新设计实验台搭接、运行,满足机构运动要求。
2、根据设计机构创新方案、画出机构运动简图,在ZBS-C机构创新设计实验台搭接、运行,满足机构设计要求。
三、实验原理1、杆组的概念由于平面机构具有确定运动的条件是机构的原动件数目与机构的自由度数相等,因此机构由机架、原动件和自由度为零的从动件系统通过运动副联接而成。
将从动件系统拆成若干个不可再分的自由度为零的运动链,称为基本杆组,简称杆组。
根据杆组的定义,组成平面机构杆组的条件是:F=3n-2pL-pH=0其中构件数n,高副数PH和低副数PL都必需是整数。
由此可以获得各种类型的杆组。
当n=1,PL=1,PH=1 时即可获得单构件高副杆组,常见的有如下几种:图 3-1 单构件高副杆组因此满足上式的构件数和运动副数的组合为:n=2,4,6……,PL=3,6,9……。
已知:电机的功率W=5.5KW,电机的高速转动(960 转/分)牛头刨床曲柄转速50r/min,工作行程310mm,行程速比系数1.46.连杆与导杆之比LDE/LCD=0.25 曲柄与水平线的夹角120°设计要求:电动机轴与曲柄轴2平行,刨刀刀刃D点与铰链点C的垂直距离为50mm,使用寿命10年,每日一班制工作,载荷有轻微冲击。
允许曲柄2转速偏差为±5%。
要求导杆机构的最大压力角应为最小值;凸轮机构的最大压力角应在许用值[α]之内,摆动从动件9的升、回程运动规律均为等加速等减速运动。
执行构件的传动效率按0.95计算,系统有过载保护。
按小批量生产规模设计。
1、确定传动装置的类型,画出机械系统传动简图。
2、选择电动机,进行传动装置的运动和动力参数计算。
3、传动装置中的传动零件设计计算。
4、绘制传动装置中减速器装配图一张(A0)。
5、绘制减速器箱体、齿轮及轴的零件图各一张(A1)。
6、编写设计计算说明书一份。
完成以上工作后准备机械设计部分的答辩。
牛头刨床机构的分析与综合一、课程设计的目的和任务1、目的机械原理课程设计是培养学生掌握机械系统运动方案设计能力的技术基础课程,它是机械原理课程学习过程中的一个重要实践环节。
其目的是以机械原理课程的学习为基础,进一步巩固和加深所学的基本理论、基本概念和基本知识,培养学生分析和解决与本课程有关的具体机械所涉及的实际问题的能力,使学生熟悉机械系统设计的步骤及方法,其中包括选型、运动方案的确定、运动学和动力学的分析和整体设计等,并进一步提高计算、分析,计算机辅助设计、绘图以及查阅和使用文献的综合能力。
2、任务本课程设计的任务是对牛头刨床的机构选型、运动方案的确定;对导杆机构进行运动分析和动态静力分析。
并在此基础上确定飞轮转惯量,设计牛头刨床上的凸轮机构和齿轮机构。
二、工作原理与结构组成牛头刨床的简介牛头刨床是用于加工中小尺寸的平面或直槽的金属切削机床,多用于单件或小批量生产。
刨床进给机构设计学生姓名学号所属学院机械电气化工程学院专业农业机械化及其自动化班级指导老师日期前言牛头刨床是刨削类机床中应用较广的一种,主要用于单件小批量生产中刨削中小型工件上的平面、成形面和沟槽等。
一般牛头刨床工作原理是当曲柄匀速转动时,摇杆左右摆动,使牛头刨床的刨刀沿着固定的轨迹运动,通过曲柄转动驱动刨头作往复移动。
牛头刨床机构具有急回的特性,即牛头刨床工作时刨刀慢慢向前移动,完成一次工作后刨刀急速返回原来的位置。
如何实现刨头在切削行程中速度平稳,难以凭经验确定。
而且机构的几何参数对切削行程刨头速度的平稳性的影响,也难以直接判断。
为了确定牛头刨床进给机构的设计是否满足要求,就必须对其进行分析。
牛头刨床是用棘轮机构完成进给运动。
工作台进给量的大小是用棘轮外圈加一档环,在摇杆固定摆角的范围内盖住棘轮上一部分牙齿来进行调整。
由于正反向进给要来回调整棘轮外圈挡环,操作比较麻烦。
为此,可以对进给机构进行改进,不用挡环来回调整,而只要将棘爪调整180度,就能实现正反进给。
牛头刨床主要用于单件小批量生产中刨削中小型工件上的平面、成形面和沟槽及一些成形面。
滑枕带着刨刀作直线住复运动的刨床,因滑枕前端的刀架形似牛头而得名。
中小型牛头刨床的主运动大多采用曲柄摇杆机构(见曲柄滑块机构)传动,故滑枕的移动速度是不均匀的。
大型牛头刨床多采用液压传动,滑枕基本上是匀速运动。
滑枕的返回行程速度大于工作行程速度。
由于采用单刃刨刀加工,且在滑枕回程时不切削,牛头刨床的生产率较低。
刨床进给机构实现的功能是除水平进给可自动外,垂直进给、非工作时工作台的水平移动和垂直升降均需手动。
关键词:刨床;进给;工作台;水平移动;垂直升降目录1 牛头刨床简要介绍 (1)1.1 牛头刨床外形图 (1)1.2 牛头刨床解析 (1)1.3 牛头刨床工作原理 (1)1.4 牛头刨床的组成 (2)2 牛头刨床进给机构总体设计 (2)2.1 牛头刨床进给机构方案的确定 (2)2.2 牛头刨床进给机构的组成 (3)2.3 牛头刨床进给机构工作原理 (3)2.4 牛头刨床垂直进给机构方案的确定 (4)3 连杆机构部分设计 (4)3.1 曲柄摇杆机构解析 (4)3.2 曲柄摇杆机构运动分析 (5)3.3 建立数学模型 (5)4 棘轮机构部分的设计 (6)4.1 棘轮机构工作原理 (6)4.2 棘轮机构的特点及应用 (7)4.3 棘轮机构的设计要点 (7)4.4 棘轮机构方案的确定 (8)4.5 棘轮和棘爪的设计 (8)4.6 棘爪架 (10)4.7 棘轮罩 (10)5 螺旋机构(丝杠副)部分的设计 (10)5.1 螺旋机构的选型 (10)5.2 螺旋机构基本参数的确定 (11)6 工作台垂直进给机构部分的设计 (11)6.1 工作台垂直进给机构分析 (11)6.2 工作台垂直进给机构原理 (11)总结 (13)致谢 (14)参考文献 (15)1 牛头刨床简要介绍1.1 牛头刨床外形图牛头刨床主要由床身、滑枕、刀架、工作台、横梁等组成如图1-1所示。
单自由度机械系统动力学作业题目:图1所示为一牛头刨床。
各构件长度为:1110L mm =,3540L mm =,4135L mm =;尺寸580H mm =,1380H mm =。
导杆3重量3200G N =,质心3S 位于导杆中心,导杆绕3S 的转动惯量23 1.1J kg m =⋅。
滑枕5的重量5700G N =。
其余构件重量均可不计。
电动机型号为Y100L2-4,电动机轴至曲柄1的传动比23.833i =,电动机转子及传动齿轮等折算到曲柄上的转动惯量21133.3J kg m =⋅。
刨床的平均传动效率0.85η=。
空行程时作用在滑枕上的摩擦阻力50f F N =,切削某工件时的切削力和摩擦阻力如图2所示。
1)求空载启动后曲柄的稳态运动规律; 2)求开始刨削工件的加载过程,直至稳态。
图1 牛头刨床 图2 牛头刨床加工某工件时的负载图 解:(1)运动分析可以用解析法列出各杆角速度、各杆质心速度的表达式。
但为简便起见,现调用改自课本附录Ⅰ中的Matlab 子程序来进行计算。
图1中给出了构件和运动副的编号。
先调用子程序crank 分析点②的运动学参数,再调用子程序vosc 进行滑块2—导杆3这一杆组的运动学分析,然后再调用子程序vguide 进行小连杆4—滑枕5这一杆组的运动学分析。
这一段的Matlab 程序如下:crank(1,2,L(1),TH(1),W(1)); vosc(2,3,4,L(3)); vguide(4,5,L(4)); 其中:L(i)、TH(i)、W(i)分别表示第i 个杆的长度、位置角、角速度。
(2)等效转动惯量和等效力矩取曲柄1为等效构件,等效转动惯量为2223335513111()()()S e J J J G v G v g g ωωωω=+++ (a) 式中:g 为重力加速度,3S v 为导杆3质心的速度,5v 为滑枕的速度。
等效驱动力矩可由电动机机械特性导出,设m M 、de M 分别为电动机输出力矩和等效驱动力矩,两者有如下关系:de m M iM = (b)式中i 为电动机轴和曲轴间的传动比。
牛头刨床的连杆机构运动分析0 前言机构运动分析的任务是对于结构型式及尺寸参数已定的具体机构,按主动件的位置、速度和加速度来确定从动件或从动件上指定点的位置、速度和加速度。
许多机械的运动学特性和运动参数直接关系到机械工艺动作的质量,运动参数又是机械动力学分析的依据,所以机构的运动分析是机械设计过程中必不可少的重要环节。
以计算机为手段的解析方法,由于解算速度快,精确度高,程序有一定的通用性,已成为机构运动分析的主要方法。
连杆机构作为在机械制造特别是在加工机械制造中主要用作传动的机构型式,同其他型式机构特别是凸轮机构相比具有很多优点。
连杆机构采用低副连接,结构简单,易于加工、安装并能保证精度要求。
连杆机构可以将主动件的运动通过连杆传递到与执行机构或辅助机构直接或间接相连的从动件,实现间歇运动,满足给定的运动要求,完成机器的工艺操作。
牛头刨床是一种利用工作台的横向运动和纵向往复运动来去除材料的一种切削加工机床。
工作台的纵向往复运动是机床的主运动,实现工件的切削。
工作台的横向运动即是进给运动,实现对切削精度的控制。
本文中只分析纵向运动的运动特性。
牛头刨床有很多机构组成,其中实现刨头切削运动的六连杆机构是一个关键机构。
刨床工作时,通过六杆机构驱动刨刀作往复移动。
刨刀右行时,当刨刀处于工作行程时;要求刨刀的速度较低且平稳,以减小原动机的容量和提高切削质量。
当刨刀处于返回行程时,刨刀不工作,称为空行程,此时要求刨刀的速度较高以提高生产率。
由此可见,牛头刨床的纵向运动特性对机床的性能有决定性的影响。
1 牛头刨床的六连杆机构牛头刨床有很多机构组成,其中实现刨头切削运动的六杆机构是一个关键机构。
图1所示的为一牛头刨床的六连杆机构。
杆1为原动件,刨刀装在C点上。
假设已知各构件的尺寸如表1所示,原动件1以等角速度ω1=1rad/s沿着逆时针方向回转,要求分析各从动件的角位移、角速度和角加速度以及刨刀C点的位移、速度和加速度的变化情况。
二、牛头刨床的运动、动力分析[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)作为计算式。
(2) 速度分析将式(5)对时间t 求导可得:)2(44)2(3343πϕπϕϕϕ+++=i i E e l e l x(10)应用欧拉公式将(7)的虚部和实部分离得:444333c o s c o s 0ϕϕϕϕ l l --= (11) 444333s i n s i n ϕϕϕϕ l l x E --= (12)由式(11)可得:444333cos cos ϕϕϕϕl l =- 对设计方案进行分析可知0cos 4≠ϕ,否则无法完成行程所以 443334cos cos ϕϕϕϕl l -= (13)又由式(4)可得:()Sl AB 3113cos ϕϕωϕ-= (14)其中的S 的表达式可由(2)求得:31sin sin ϕϕAB l a S +=(由设计方案显然3sin ϕ不为零) (15)由式(12)、(1)可得到由1ϕ确定的E v 的公式:()44333cos sin ϕϕϕϕ--== l xv E E (16)(3) 加速度分析将式(9)对时间t 求导可得:()()⎪⎭⎫ ⎝⎛++⎪⎭⎫ ⎝⎛+++++=2442442332334433πϕπϕπϕπϕϕϕϕϕi i i i E e l e l e l e l x(17)应用欧拉公式将(15)的虚部和实部分离得:34432443333233sin cos sin cos ϕϕϕϕϕϕϕϕ l l l l x E ----= (18) 44442443333233cos sin cos sin 0ϕϕϕϕϕϕϕϕ l l l l +-+-= (19) 由式(16)可得:44424433332334cos sin cos sin ϕϕϕϕϕϕϕϕl l l l +-= (20) 由式(15)、(17)可得到由1ϕ确定的E a 的公式:()()4424343233433cos cos sin ϕϕϕϕϕϕϕϕl l l xa E E +---== (21) 其中由式(14)可求得()SS l 31312132sin ϕϕϕωϕ--=由式(15)可求得()3111sin ϕϕω--=l S 2、推导出曲柄所加平衡力矩111()M M ϕ=的数学式。
对机构进行受力分析可知,假设在整个力的传递过程中损耗为零,则驱动力矩做的功应当与阻力所做的功相等,若曲柄沿逆时针旋转,令当刨头到达极限位移时杆3与y 轴夹角为θ 则有:3322sin l Hl H==θ即:32arcsinl H=θ 所以当[)[]πθπθπϕ2,2,01-⋃+∈时,阻抗力1000-=r F ;当[)θπθπϕ-+∈2,1时,阻抗力0=r F ; 由受力分析可得:E r vF M =11ϕ所以可得:()41433311cos sin ϕωϕϕϕω--==l F v F M r Er四、打印源程序代码 运用的语言种类:C 语言程序名称:牛头刨床的运动、动力分析#include<stdlib.h> #include<stdio.h> #include<math.h> #define NUM 37main() {float XE[NUM]; float VE[NUM]; float AE[NUM]; float M1[NUM];float a=270, b=520,L3=550,L4=100,H=420,N1=27,PI=3.14159;double PQ0,Q1,PQ1,PQ3,PQ4,S,L1,D1PQ3,D2PQ3,D1PQ4,D1S,W1,FR; int i,j;/* 计算曲柄长度L1和角速度W1*/ L1=(a*H)/(2*L3); W1=(N1*2*PI)/60;printf("\n Q1 XE VE AE M1 ");for(i=0;i<=36;i++){ Q1=10*i;PQ1=(Q1*PI)/180;/*分情况讨论PQ3的取值*/if(Q1==90||Q1==270) PQ3=PI/2; else{ if(atan((a+L1*sin(PQ1))/(L1*cos(PQ1)))<0)PQ3=atan((a+L1*sin(PQ1))/(L1*cos(PQ1)))+PI;elsePQ3=atan((a+L1*sin(PQ1))/(L1*cos(PQ1)));}PQ4=PI-asin((b-L3*sin(PQ3))/ L4);/*计算XE*/XE[i]=L3*cos(PQ3)+ L4*cos(PQ4);/*分情况讨论S的取值*/if(Q1==90)S=a+L1;else{if(Q1==270)S=a-L1;elseS=(L1*cos(PQ1))/cos(PQ3);}D1PQ3=(L1*W1*cos(PQ1-PQ3))/S;/*计算VE的值并将其值存放在VE[i]中*/VE[i]=(-D1PQ3*L3*sin(PQ3-PQ4))/cos(PQ4);/*计算AE的值并将其值存放在AE[i]中*/D1PQ4=(-L3*D1PQ3*cos(PQ3))/(L4*cos(PQ4));D1S=-W1*L1*sin(PQ1-PQ3);D2PQ3=(W1*W1*L1*sin(PQ3-PQ1)-2*D1PQ3*D1S)/S;AE[i]=(L3*D2PQ3*sin(PQ4-PQ3)-L3*D1PQ3*D1PQ3*cos(PQ4-PQ3)+L4*D1PQ4*D1PQ4)/cos(PQ4);/*计算并分区间讨论FR的值*/PQ0=asin(H/(2*L3));if(PQ1>=(PI+PQ0)&& PQ1<=(2*PI-PQ0)) FR=0;else FR=-1000;/*计算M的值并将其值存放在M[i]中*/M1[i]=(-FR*D1PQ3*L3*sin(PQ3-PQ4))/(W1*cos(PQ4));/*输出计算结果*/printf("\n%8.0f",Q1);printf("\t%8.2f",XE[i]);printf("\t%8.2f",VE[i]);printf("\t%8.2f",AE[i]);printf("\t%10.2f",M1[i]);}printf("\n Please enter the option to exit(1):");scanf("%d",&j);if(j==1) exit;return 0;} 程序运行结果:Q1(°) XE(mm)VE(mm/s) AE(mm/s*s) M(N*mm)0 96.38 -189.22 -1116 66922.05 10 82.92 -244.43 -931.12 86448.69 20 66.46 -287.24 -783.92 101592.03 30 47.64 -321.09 -661.65 113561.47 40 26.94 -348.6 -555.14 123291.81 50 4.7 -371.61 -458.05 131430.81 60 -18.86 -391.18 -366.28 138351.81 70 -43.53 -407.65 -277.56 144178.36 80 -69.12 -420.75 -190.86 148810.81 90 -95.39 -429.69 -105.57 151973.62 100 -122.06 -433.37 -20.74 153273.72 110 -148.77 -430.67 65.36 152317.5 120 -175.08 -420.5 155.22 148723.31 130 -200.51 -402.13 251.79 142226.23140 -224.55 -375.14 358.26 132679.95 150 -246.65 -339.43 478.26 120047.27 160 -266.27 -295.08 616.45 104361.7 170 -282.9 -242.18 779.33 85655.47 180 -295.99 -180.48 975.78 63830.48 190 -304.98 -108.67 1216.98 38433.07 200 -309.14 -23.5 1514.77 8312.74 210 -307.48 81.22 1876.72 0 220 -298.52 214.45 2291.31 0 230 -280.25 383.83 2687.4 0 240 -250.46 585.24 2865.3 0 250 -207.95 787.77 2483.9 0 260 -154.42 931.03 1282.27 0 270 -95.39 960.49 -527.49 0 280 -38.2 877.34 -2229.77 0 290 11.57 728.24 -3192.72 0 300 51.25 556.14 -3361.88 0 310 80.25 384.84 -3048.03 0 320 99.04 226.76 -2558.31 0 330 108.66 88.83 -2074.39 0 340 110.49 -25.88 -1669.96 9152.98 350 105.94 -117.66 -1355.59 41612.64 360 96.38 -189.22 -1116 66921.41刨头位置曲线图-350-300-250-200-150-100-5005010015015913172125293337Q1(10度)X E (m m )刨头速度变化曲线图-600-400-20002004006008001000120014710131619222528313437Q1(10度)V E (m m /s )刨头加速度变化曲线图-4000-3000-2000-100001000200030004000135791113151719212325272931333537Q1(10度)A E (m m /(s *s ))曲柄所加平衡力矩变化曲线图020000400006000080000100000120000140000160000180000135791113151719212325272931333537Q1(10度)M 1(N *m m )。