直动式固定凸轮与连杆机构的设计

1、湖南工业大学直动式固定凸轮与连杆机构的设计设计者所在院（系）：湖南工业大学专业：机械设计制造及其自动化班级学号：指导老师：时间：2015年12月27日湖南工业大学耳录一、课程设计的目的3二、设计内容与步骤41、设计内容42.设计步骤4三、设计要求6四、设计指导71、概述72、基本参数93、设计步聚111）确定驱动方案112）确定e113）确定h124）确定：125）确定'126）求算b1、b2127）设计凸轮廊线148）检验压力角16五、参数优化18六、结论19七、参考文献20八、附图21摘要包装设计课程设计是在完成机械设计课程学习后，一次重要的实践性教学环节。是高等工科院校大多数专业

2、学生第一次较全面的设计能力训练，也是对机械设计课程的全面复习和实践。其目的是培养理论联系实际的设计思想，训练综合运用机械设计和有关选修课程的理论，结合生产实际分析和解决工程实际问题的能力，巩固、加深和扩展有关机械设计方面的知识。本次设计的题目是直动式固定凸轮与连杆机构的设计。根据题目要求和机械设计的特点作者做了以下几个方面的工作：根据有关参数进行计算或编写有关设计计算程序；利用程序设计的方法输生结果并自动生成图形；画由装配图及其主要零件图；完成设计计算说明书。一、课程设计的目的包装机械设计课程设计是本课程各教学环节中重要的一环，它让学习者联系实际进一步深入理解、掌握所学的理论知识。其基本目的是

3、：培养理论联系实际的设计思想，训练综合运用包装机械和有关先修课程的理论，结合生产实际分析和解决工程实际问题的能力，巩固、加深和扩展有关包装机械设计方面的知识。通过制订设计方案，合理选择裹包机中块状物品推送机构和零件类型，正确计算零件工作能力、确定尺寸和选择材料，以及较全面地考虑制造工艺、使用和维护等要求，之后进行结构设计，达到了解和掌握机械零件、包装机械经常采用的机构的设计过程和方法。进行设计基本技能的训练。例如计算、绘图、熟悉和运用设计资料（手册、图册、标准和规范等）以及使用经验数据、进行经验估算和处理数据的能力。10二、设计内容与步骤1、设计内容以裹包机中块状物品推送机构的典型机构一一固定

4、凸轮与连杆组合机构为题。课程设计通常包括如下内容：读懂块状物品推送机构典型机构一一固定凸轮与连杆组合机构，了解设计题目要求；分析该块状物品推送机构设计的可能方案；具体计算和设计该方案中机构的基本参数；进行机体结构及其附件的设计；绘制装配图及零件工作图；编写计算说明书以及进行设计答辩。2.设计步骤(1)设计准备认真研究设计任务书，明确设计要求、条件、内容和步骤；通过阅读有关资料、图纸、参观实物或模型、观看电视教学片、挂图以及推送机构进行拆装实验等，了解设计对象；复习有关课程内容，熟悉零部件的设计方法和步骤；准备好设计需要的图书、资料和用具；拟定设计计划等。(2)推送机构装置的总体设计决定推送机构

5、装置的方案；选择机构的类型，计算机构装置的运动参数。(3)装配图设计计算和选择机构的参数；确定机体结构和有关尺寸；绘制装配图草图；选择计算轴承和进行支承结构设计；进行机体结构及其附件的设计；完成装配图的其他要求；审核图纸。(4)零件工作图设计(5)整理和编写计算说明书(6)设计总结和答辩三、设计要求在课程设计之前，准备好必要的设计手册或参考资料，以便在设计过程中逐步去学习查阅资料。确定设计题目后，至少应复习在课程中学过的相关内容。完成本课程设计的具体要求如下：1、设计说明书要全面反映设计思想、设计过程和结论性认识。其工艺设计要有文字、计算、公式来源、参数选取的资料名称或代号、图表（草图）。说明

6、书用A4纸打印，约20页左右，并装订成册。2、设计图样按“机械制图”、“公差与配合”等国家标准完成。3、零件图按生产图样要求完成，零件的有关精度和技术要求要有合理的标注或说明。设计过程中，提倡独立思考、深入钻研，主动地、创造性地进行设计，反对不求甚解、照抄照搬或依赖老师。要求设计态度严肃认真、有错必改，反对敷衍塞责，容忍错误的存在。只有这样，才能保证课程设计达到教学基本要求，在设计思想、设计方法和设计技能等方面得到良好的训练。四、设计指导裹包机所包装的产品，绝大多数是单件或多件集合而成的块状物品。包装作业线中前后机之间物品的输送、换向、排列组合，及单机内部的物品移动等，需要用各种各样的机构或装

7、置完成。以下是几种典型的推送块状物品的组合机构一一固定凸轮与连杆组合机构。1、概述图1所示，是该机构的结构简图，用于香皂、糖果等裹包机中，将物品向上推送较大距离。原动杆件AB按逆时针方向转动，驱动较销C上的滚动轴承6在固定梢凸轮4的梢内运动，再通过连杆CD使推送杆（即滑块）2按预定规律作上下往复移动。这种直动从杆类型的固定凸轮与连杆组合机构相当于连杆长度可变的曲柄滑块机构，曲柄为AR滑块为推送杆，连杆为BD,在运动过程中连杆BD的长度是变化的。图1直动从动杆类型的固定凸轮和连杆组合的推送机构结构简图1推送板2推送杆3导轨4固定槽凸轮5支座6滚动轴承7导轨图2固定凸轮与连杆组合机构示意图1-推料

8、板2-推料杆3-固定凸轮4-滚子以所示推送机构，还有2种运动方式：一是是曲柄AB推着杆件BC运动，杆件BC承受压力；二是曲柄AB拉着BC杆运动，杆件BC承受拉力。这是两种不同的驱动方案。当然，无论是前者还是后者，都可以在两种驱动方案中任意选择，本题我们选择直动从动杆类型的固定凸轮和连杆组合机构。2、基本参数为研究方便，特规定：以曲柄回转中心A为坐标的原点，并作x、y轴。对于直动从动类型（见图1所示），y轴与从动杆的运动方向平行；对于摆动从动杆类型（见图3所示）y轴与钱销D的两个运动极限位置之连线D0D平行。考虑到曲柄有两种转向，又规定y轴的正轴逆着曲柄转向旋转900后所得轴为厘轴的正轴，于是，

9、前者1.轴的正轴向右，而手者则向左。基本参数有：ey轴与DD线的间距，简称偏心距；h一钱销D至7轴的最小距离；a一曲柄AB长；bi、b2一杆件BGCD<6从动杆升程运动起始时刻的曲柄位置AB和y轴负轴的夹角，6=180°-ZBAY钱销B和D的距离用b表示，b=BD它的最大值和最小值分别用bmax、bmin表不。已知参数:从动杆的行程%升程和降程的运动规律100130二0©130：100二升程和降程按五次多项式运动表1初始参数表固定凸轮与连杆组合机构的特点是，从动杆的运动可以象凸轮机构的从动杆那样实现停留和按照定规律(如五次多项式)运动。从动杆的行程、动停时间、运动速度

10、由工艺要求预先给定。这样，当参数e、h、口、8,确定后，每一运动时刻的b值及bmax、bmin值也随之确定。显然，bi、b2应满足下式h+b2=bmax：*(1)b2一bi=bmin湖南工业大学因此，应根据从动杆的运动规律和确定的c、h、口、s、l值,先计算出bmax、bm%然后用下式求算bl、b?值:bl = 2 (bmax1.，b2 = (bmax2- bmin )bmin )(2)3、设计步聚1）确定驱动方案它对凸轮的压力角机构的传动效率影响较大。应根据运动要求确定之。用下列符号表示运动要求：S一分别为直动总行程；邛i一升程运动对应的曲柄转角；%最高位置停留对应的曲柄转角；我一降程运动对

11、应的曲柄转角；Q一最低位置停留对应的曲柄转角；Q+中2+%+Q=3600,当QA我时，先用曲柄AB拉着杆件BC运动的方案；当Q时，应选用曲柄AB推着BC杆运动的方案。巴=九时,可任选其中一种方案。2）确定e直动从动杆，取e=0ii湖南工业大学3)确定h从结构紧凑和减小凸轮压力角考虑，应将h值取小些。但h值愈小，对从动杆驱动力的压力角也愈大。通常取hnSi,h=3Smo4)确定3若a值过小，会使凸轮压力角明显增大，甚至不能实现预期动动。可取a=0.8Sm。5)确定5其值对凸轮的压力角影响极大，6过小，尤其是过大，会使压力角急剧增加。在前述参数确定后，最好将6优化，目标函数为aim(')、

12、aimmin式中aim为凸轮的最大压力角。8=-46)求算bi、b2须先求算bmax、bmin。(1)直动从动杆类型参阅图4,依据较销B、D的坐标，可建立它们之间距离的公式。B的坐标为XB=asin()(3)yB:-acos()D的坐标为:Xd=e（4）yD=hS式中中一曲柄转角，取升程起始时的中=0;S为与中相对应的从动杆位移，即较销D至其最低位置的距离。S值分为升程（邛二0%）、最高位置停留（中二中1邛1+中2）、降程（邛二中什干2邛1+中2+中3）、最低位置停留（中二中什中2+93360°）四个阶段求算。b值为：b=（xB-Xd）（yB-yD）（5）将式（3）、（4）代入式（5

13、）,求算bnax、bmin,然后用式（2）算得b1、b2O因此,应根据从动杆的运动规律和确定的e,h,a（1,先计算出bmax,bmin,利用MATLAB强大的可视化功能，绘出b随。的曲线图，见图4,程序如附录程序一.可以很方便地得到：bmax=397.8748mm,bmin=256.9556mm。图4 B D长度曲线图然后再利用（2）式求得4、4分别为：bl=70.4613mm,b2=327.4171mm7）设计凸轮廊线固定凸轮的理论廊线就是滚子中心C的运动轨迹线，根据较销RD的位置及bi、b2值可确定C的位置。令钱销B,D的连线BD与D0D1线（或y轴）的夹角为3BD与CD的夹角为（3,则

14、：r.XB-XDXB-XD(6)口-arcsinarctanbyarcosbb2bb2显然，xB>xD时0为正值，反之则为负值，而（3始终为正值钱销C的坐标为:xC=xDb2sinG-）yc=yDb2cosG一：该式运算符号“+”和“一”的确定原则是：令b=bmax时的力为想，b=bmin时的4为。，在e=%区间，取“-”号；在4=0*和4=em360:区间，取“+”。由图4可得到©mC一,一,、，'C一一.、，=82.8496（1.4460弧度），0m=216.7442（3.7829弧度）。在MATLAB,编程计算得到凸轮理论轮廓即C点的轨迹上各点的坐标，并绘出C点的

15、轨迹线，见图5,程序如附录程序二。：806040200*20-6040-200204060 S0图5凸轮轮廓线图8)检验压力角凸轮机构设计中，初始参数选择不当会导致凸轮轮廓设计不合理从而导致从动件运动不能按预定规律进行，或者运行中动力学性能欠佳。(1)凸轮的压力角a参阅图2,a1为Pc和Vc的夹角。Pc为驱动校销运动的力(不考虑摩擦力)，BC重合，vc为钱链C的运动方向一与C点的凸轮廊线切线重合。用K、K2分别表示Pcvc的斜率，则(9)(10)Nc-VbKi=Xc-Xbdyc_dyc-dK2一-dxcdxc''da1 = arctgK2 - K11 K1K2(11)37计算式

16、K2时需求函数yc(xc)或ycT),XcS)的一阶导数，这些函数很难用显fi方程表示，求解一阶导数更加繁琐，因此本文利用MATLAB方便的向量运算功能来求解K2。将凸轮细分为若干个点，第i点的斜率可采用式(12)来求解(12)")二'。Xc(i-1)-Xc(i1)在MATLABtt式(12)进行求解非常方便，其程序代码简短易懂如下:Fori=2:length(fa)-1%fa凸轮某一段的转角K2(i)=(yc(i-1)-yc(i+1)/(Xc(i-1)-xji+1)end应保证ai的最大值不超过许用值，即aimw(Q)。可取)=450利用附录程序三，通过matlab得到该机

17、构的alm=44.1427<45,设计合理。(2)从动杆的压力角a2a?为CD杆对较销D的驱动力(不考虑摩擦力)与D的运动方向的夹角。对于直动从动杆，显然8=日土P(13)对于AB推动BC的驱动方案，在中=外邛m区间取“一”号，其它区间取“+”号。应使a2的最大值a2maxE瓦。直动从动杆以=30°。利用附录程序四，通过matlab计算得到该机构的a2max=27.3036<30通过对压力角的检验可知，以上参数的设计合理。五、参数优化5值是从动杆升程运动起始时刻的曲柄位置AB0和y轴负轴的夹角，通过改变5值，观察凸轮压力角和从动杆压力角的变化，见表2。6/()543210

18、-1-2-3almax61.3058.6756.3154.1852.2250.4148.7247.1245.60a2max25.3525.725.7926.0226.2426.4526.6726.8827.096/()-4-5-6-7-8-9-10-11-12almax44.1444.2245.0245.8246.6247.6148.4349.4450.45a2max27.3027.5127.7527.9228.1228.2328.4228.6128.83表26与压力角的关系S值对从动杆压力角影响很小，但对凸轮压力角影响极大。当5为5°时，凸轮压力角&超出了许用范围；随着5的

19、减小，a变小；当"为-4°时，为达到最小值44.14。；随着S继续减小，为增大。可见当以上其它参数不变的情况下，S为-4°,该机构的动力性能达到最优。六、结论1）运用MATLAB强大的编程功能实现了凸轮连杆机构的计算机辅助设计计算，得到了连杆的长度和凸轮曲线的轨迹，方便地检验了凸轮压力角和从动杆压力角，并分析了关键参数S值与压力角的关系，最后找到最优方案，从而提高了机构的动力性能。2）借用MATLAB简单的数组计算功能，提出了凸轮压力角有效的计算方法。七、参考文献1许林成彭国勋.包装机械M.长沙：湖南大学出版社，1989.2许林成.包装机械设计与原理M.上海：上海

20、科学技术出版社，1988.3孙智慧晏祖根.包装机械概论M.北京：印刷工业出版社，2012,(7).4肖乾周新建.凸轮机构的计算机辅助设计与运动仿真分析J华东交通大学学报,2006,23(4):103-105.5周开利邓春晖.MATLA廊用教程M.北京：北京大学出版社，2007,36张义智郭连考.给纸机递纸吸嘴机构共腕凸轮设计J.包装工程,2007,28：104-105.7王家文王皓刘海.MATLAB7.0编程基础M.北京：机械工业出版社,2005.8孙智慧高德.包装机械M.北京：中国轻工业出版社,2010,(1).9刘筱霞.包装机械M.北京：化学工业出版社,2007,(1).10高德.包装机械

21、设计M.北京：化学工业出版社,2005,(8).11刘月花郭谆钦.机械设计基础M.长春：东北师范大学出版社，2010,(3).八、附录附录程序一：clc;clear;globalblb2FAMAXFAMINsitaEfa1=130*pi/180;fa2=0*pi/180;fa3=130*pi/180;fa4=100*pi/180;sm=120;E=0;h=1.5*sm;a=0.8*sm;sita=-8*pi/180;faO1=0:0.001:fa1;faa=fa01/fa1;s=sm.*(6.*(faa)八5-15.*(faa)八4+10.*(faa)八3);xb=a.*sin(sita+fa

22、01);yb=-a.*cos(sita+fa01);xd=E;yd=h+s;b01=sqrt(xb-xd).八2+(yb-yd).八2);figure;plot(fa01,b01);holdon;fa02=fa2+fa1:0.001:fa3+fa1+fa2;s=sm.*(1-(6.*(fa02-fa1-fa2)/fa3)八5-15.*(fa02-fa1-fa2)/fa3)八4+10.*(fa02-fa1-fa2)/fa3)八3);xb=a*sin(sita+fa02);yb=-a*cos(sita+fa02);xd=E;yd=h+s;b02=sqrt(xb-xd).八2+(yb-yd).八2)

23、;plot(fa02,b02);holdon;bmin=min(b01,b02)bmax=max(b01,b02)b1=(bmax-bmin)/2b2=(bmax+bmin)/2famax=b01(1);fo门=2:2095iffamax<b01(i)famax=b01(i);FAMAX=fa01(i);endendfamax,FAMAXfamin=b02(1);fori=2:2095iffamin>b02(i)famin=b02(i);FAMIN=fa02(i);endendfamin,FAMIN附录程序二:famax=1.3990;famin=3.8449;b1=70.4613

24、;b2=327.4171;fa1=130*pi/180;fa2=0*pi/180;fa3=130*pi/180;fa4=100*pi/180;sm=100;e=0;h=3*sm;a=0.8*sm;sita=-4*pi/180;fa01=0:0.002:famax-0.002;faa=fa01/fa1;s=sm.*(6.*(faa)八5-15.*(faa)八4+10.*(faa)八3);xb=a.*sin(sita+fa01);yb=-a.*cos(sita+fa01);xd=e;yd=h+s;b=sqrt(xb-xd).八2+(yb-yd).八2);theta1=asin(xb-xd)./b)

25、;beta1=acos(b.A2+b2.A2-b1.A2)./(2.*b.*b2);xc=xd+b2*sin(theta1+beta1);yc=yd-b2*cos(theta1+beta1);plot(xc,yc,'r');holdon;fa02=famax+0.002:0.002:fa1;faa=fa02/fa1;s=sm.*(6.*(faa)八5-15.*(faa)八4+10.*(faa)八3);%plot(fa02,s)xb=a.*sin(sita+fa02);yb=-a.*cos(sita+fa02);xd=e;yd=h+s;b=sqrt(xb-xd).八2+(yb-y

26、d).八2);theta2=asin(xb-xd)./b);beta2=acos(b.A2+b2.A2-b1.A2)./(2.*b.*b2);xc=xd+b2*sin(theta2-beta2);yc=yd-b2*cos(theta2-beta2);plot(xc,yc,'r');holdon;fa03=fa2+fa1:0.002:famin;s=sm.*(1-(6.*(fa03-fa1-fa2)/fa3)A5-15.*(fa03-fa1-fa2)/fa3)A4+10.*(fa03-fa1-fa2)/fa3)A3)；%piot(fa03,s)xb=a*sin(sita+fa03

27、);yb=-a*cos(sita+fa03);xd=e;yd=h+s;b=sqrt(xb-xd).A2+(yb-yd).A2);theta=asin(xb-xd)./b);beta3=acos(b.A2+b2A2-b1.A2)./(2.*b.*b2);xc=xd+b2*sin(theta-beta3);yc=yd-b2*cos(theta-beta3);plot(xc,yc,'r');holdon;fa04=famin+0.002:0.002:fa1+fa2+fa3;s=sm.*(1-(6.*(fa04-fa1-fa2)/fa3)八5-15.*(fa04-fa1-fa2)/fa

28、3)八4+10.*(fa04-fa1-fa2)/fa3)八3);%plot(fa04,s)xb=a*sin(sita+fa04);yb=-a*cos(sita+fa04);xd=e;yd=h+s;b=sqrt(xb-xd).八2+(yb-yd).八2);theta=asin(xb-xd)./b);beta4=acos(b.A2+b2.A2-b1.A2)./(2.*b.*b2);xc=xd+b2*sin(theta+beta4);yc=yd-b2*cos(theta+beta4);plot(xc,yc,'r');holdon;fa05=fa1+fa2+fa3:0.002:2*pi

29、;s=0;%plot(fa05,s);xb=a*sin(sita+fa05);yb=-a*cos(sita+fa05);xd=e;yd=h+s;b=sqrt(xb-xd).八2+(yb-yd).八2);theta=asin(xb-xd)./b);beta5=acos(b.A2+b2.A2-b1.A2)./(2.*b.*b2);xc=xd+b2*sin(theta+beta5);yc=yd-b2*cos(theta+beta5);plot(xc,yc,'r');holdon;附录程序三：clc;clear;globalsitaEhasmfamax=1.4090;famin=3.7

30、329;b1=67.8307;b2=182.2473;fa1=130*pi/180;fa2=0*pi/180;fa3=130*pi/180;fa4=100*pi/180;sm=100;e=0;h=1.5*sm;a=0.8*sm;sita=-2*pi/180;faO1=0:0.001:famax;faa=fa01/fa1;s=sm.*(6.*(faa)八5-15.*(faa)八4+10.*(faa)八3);xb=a.*sin(sita+fa01);yb=-a.*cos(sita+fa01);xd=e;yd=h+s;b=sqrt(xb-xd).八2+(yb-yd).八2);theta1=asin(

31、xb-xd)./b);beta1=acos(b.A2+b2.A2-b1.A2)./(2.*b.*b2);xc=xd+b2*sin(theta1+beta1);yc=yd-b2*cos(theta1+beta1);holdon;k11=(yc-yb)./(xc-xb);fori=2:length(fa01)-1k21(i)=(yc(i-1)-yc(i+1)./(xc(i-1)-xc(i+1);endk21(1)=k21(2);fori=1:length(fa01)-1a1(i尸abs(atan(k21(i)-k11(i)./(1+k11(i)*k21(i)*180/pi;endfa01(leng

32、th(fa01)二口;a1(ength(fa01)=a1(ength(fa01)-1);al；plot(fa01,a1);holdon;fa02=famax:0.001:fa1;faa=fa02/fa1;s=sm.*(6.*(faa)八5-15.*(faa)八4+10.*(faa)八3);%plot(fa02,s)xb=a.*sin(sita+fa02);yb=-a.*cos(sita+fa02);xd=e;yd=h+s;b=sqrt(xb-xd).八2+(yb-yd).八2);theta2=asin(xb-xd)./b);beta2=acos(b.A2+b2.A2-b1.A2)./(2.*b

33、.*b2);xc=xd+b2*sin(theta2-beta2);yc=yd-b2*cos(theta2-beta2);holdon;k12=(yc-yb)./(xc-xb);fori=2:length(fa02)-1k22(i)=(yc(i-1)-yc(i+1)./(xc(i-1)-xc(i+1);endk22(1)=k22(2);fori=1:length(fa02)-1a2(i尸abs(atan(k22(i)-k12(i)./(1+k12(i)*k22(i)*180/pi;endfa02(length(fa02)二口；a2(length(fa02)=a2(length(fa02)-1);

34、a2;plot(fa02,a2);holdon;fa03=fa2+fa1:0.001:famin;s=sm.*(1-(6.*(fa03-fa1-fa2)/fa3)八5-15.*(fa03-fa1-fa2)/fa3)八4+10.*(fa03-fa1-fa2)/fa3)八3);%plot(fa03,s)xb=a*sin(sita+fa03);yb=-a*cos(sita+fa03);xd=e;yd=h+s;b=sqrt(xb-xd).八2+(yb-yd).八2);theta=asin(xb-xd)./b);beta3=acos(b.A2+b2.A2-b1.A2)./(2.*b.*b2);xc=xd

35、+b2*sin(theta-beta3);yc=yd-b2*cos(theta-beta3);holdon;k13=(yc-yb)./(xc-xb);fori=2:length(fa03)-1k23(i)=(yc(i-1)-yc(i+1)./(xc(i-1)-xc(i+1);endk23(1)=k23(2);fori=1:length(fa03)-1a3(i尸abs(atan(k23(i)-k13(i)./(1+k13(i)*k23(i)*180/pi;endfa03(length(fa03)=;a3(length(fa03)=a3(length(fa03)-1);a3;plot(fa03,a

36、3);holdon;fa04=famin:0.001:fa1+fa2+fa3;s=sm.*(1-(6.*(fa04-fa1-fa2)/fa3)A5-15.*(fa04-fa1-fa2)/fa3)A4+10.*(fa04-fa1-fa2)/fa3)A3);%plot(fa04,s)xb=a*sin(sita+fa04);yb=-a*cos(sita+fa04);xd=e;yd=h+s;b=sqrt(xb-xd).八2+(yb-yd).八2);theta=asin(xb-xd)./b);beta4=acos(b.A2+b2.A2-b1.A2)./(2.*b.*b2);xc=xd+b2*sin(th

37、eta+beta4);yc=yd-b2*cos(theta+beta4);holdon;k14=(yc-yb)./(xc-xb);fori=2:length(fa04)-1k24(i)=(yc(i-1)-yc(i+1)./(xc(i-1)-xc(i+1);endk24(1)=k24(2);fori=1:length(fa04)-1a4(i尸abs(atan(k24(i)-k14(i)./(1+k14(i)*k24(i)*180/pi;endfa04(length(fa04)=;a4(length(fa04)=a4(length(fa04)-1);a4;plot(fa04,a4);holdon;

38、fa05=fa1+fa2+fa3:0.001:2*pi;s=0;%plot(fa05,s);xb=a*sin(sita+fa05);yb=-a*cos(sita+fa05);xd=e;yd=h+s;b=sqrt(xb-xd).八2+(yb-yd).八2);theta=asin(xb-xd)./b);beta5=acos(b.A2+b2.A2-b1.A2)./(2.*b.*b2);xc=xd+b2*sin(theta+beta5);yc=yd-b2*cos(theta+beta5);holdon;k15=(yc-yb)./(xc-xb);fori=2:length(fa05)-1k25(i)=(

39、yc(i-1)-yc(i+1)./(xc(i-1)-xc(i+1);endk25(1)=k25(2);fori=1:length(fa05)-1a5(i尸abs(atan(k25(i)-k15(i)./(1+k15(i)*k25(i)*180/pi;endfa05(length(fa05)二口;a5(length(fa05)=a5(length(fa05)-1);a5;plot(fa05,a5);holdon;max(a1,a2,a3,a4,a5)min(a1,a2,a3,a4,a5)附录程序四：clc;clear;globalsitaEhasmfamax=1.3990;famin=3.844

40、9;b1=70.4613;b2=327.4171;fa1=130*pi/180;fa2=0*pi/180;fa3=130*pi/180;fa4=100*pi/180;sm=100;e=0;h=3*sm;a=0.8*sm;sita=-4*pi/180;fa01=0:0.002:famax-0.002;faa=fa01/fa1;s=sm.*(6.*(faa)八5-15.*(faa)八4+10.*(faa)八3);xb=a.*sin(sita+fa01);yb=-a.*cos(sita+fa01);xd=e;yd=h+s;b=sqrt(xb-xd).八2+(yb-yd).八2);theta1=asi

41、n(xb-xd)./b);beta1=acos(b.A2+b2.A2-b1.A2)./(2.*b.*b2);xc=xd+b2*sin(theta1+beta1);yc=yd-b2*cos(theta1+beta1);holdon;%+从动件压力角红色c1=abs(theta1-beta1)*180/pi;plot(fa01,c1,'r-*');holdon;fa02=famax+0,002:0.002:fa1;faa=fa02/fa1;s=sm.*(6.*(faa)八5-15.*(faa)八4+10.*(faa)八3);%plot(fa02,s)xb=a.*sin(sita+fa02);yb="a.*cos(sita+fa0