牛头刨床设计基础实验设计说明书.doc

辽 宁 工 业 大 学 课 程 设 计 说 明 书（论 文）一、设计题目与原始数据1题目：牛头刨床的综合设计与分析2原始数据：刨头的行程 h=550mm行程速比系数 k=1.6机架长 lo2o3=400mm质心与导杆的比值 lo3s4/lo3b=0.5连杆与导杆的比值 lbf/lo3b=0.3刨头重心至f点距离 xs6=160mm导杆的质量 m4=15 kg刨头的质量 m6=58 kg导杆的转动惯量 js4=0.7 kgm切割阻力 fc=1300n切割阻力至o2的距离 yp=175mm构件2的转速 n2=80 rpm许用速度不均匀系数 =1/40齿轮z1、z2的模数 m12=15 mm小齿轮齿数 z1=18大齿轮齿数 z2=46凸轮机构的最大摆角 max=16凸轮的摆杆长 lo4c=140mm凸轮的推程运动角 0=60凸轮的远休止角 01=10凸轮的回程运动角 0=60凸轮机构的机架长 lo2o4=150mm凸轮的基圆半径 ro=55mm凸轮的滚子半径 rr=15mm二、牛头刨床示意图 图1三、导杆机构设计1、已知：行程速比系数 k=1.6刨头的行程 h=550mm机架长度 lo2o3=400mm连杆与导杆的比 lbf/lo3b=0.52、各杆尺寸设计如下a、求导杆的摆角：max =180（k-1）/（k+1）=180（1.6-1）/（1.6+1）=41.54 b、求导杆长：lo3b1=h/2sin（max/2）=550/2sin（41.54/2）=776mmc、求曲柄长：lo2a =lo2o3sin（max/2）=400sin21=142mm d、求连杆长：lbf=lo3blbf/lo3b=7760.3=233mme、求导路中心到o3的距离：lo3m =lo3b-lde/2=lo3b1-1-cos(max/2)/2=750mm f、取比例尺：l=0.005m/mm 在1#图纸中央画机构位置图，机构位置图见1#图纸。大致图形如下： 图2四、机构的运动分析已知：曲柄转速：n2=80rpm各构件的重心： 构件6的重心：xs6=160mm第5点：a、速度分析求va3va3 =va2=lo2a2= lo2an/30=(142/1000)803.14/30=116.58.37=1.19m/s 求va4 = + 大小： ？ 1.19 ？方向：o3a o2a o3a取v= 0.025 ，在1#图纸的左下方画速度多边形求vb用速度影像求vb=560.025=1.4m/s求vf = + 大小： ？ 0.98 ？方向： 水平 导路 bf接着画速度多边形，由速度多边形求得：vf=v=550.025=1.38m/s 方向：水平向右求44=3=va4/lo3a=390.025/（10800.005）=1.81 rad/s 方向：顺时针求va4a3va4a3= v=500.025=0.125m/s 方向：如速度图所示 b、加速度分析求aka4a3aka4a3=24va4a3=21.810.125=0.453m/s2 方向：如速度图所示 求aa3aa3=aa2=22lo2a=8.382142/1000=9.97m/s2 方向：ao2 求ana4ana4=23lo3a=1.81210800.005=1.77m/s2 方向：ao3 求aa4 + = + + 大小：1.77 ？ 9.97 0.453 ?方向：ao3 o3a o3a o3a取a= 0.10在1#图的左下方画加速度多边形，由加速度多边形求得：aa4= a=190.10=1.9m/s2 方向：如图所示求ab用加速度影像求ab=270.10=2.7m/s2 方向：如图所示求af = + + 大小： ？ 2.7 0.45 ?方向：水平 fb bfanfb= v2bf/ lbf=(12.90.025) 2 /(233/1000)=0.45 m/s2接着画加速度多边形，由加速度多边形求得：af=a =90.10=0.9 m/s2 方向：水平向右 第8点：a、速度分析求va3va3 =va2=lo2a2= lo2an/30=(142/1000)803.14/30=116.58.37=1.19m/s求va4 = + 大小： ？ 1.19 ？方向：o3a o2a o3a取v=0.025 ，在1#图纸的左下方画速度多边形求vb用速度影像求vb=0求vf = + 大小： ？ 0 ？方向：水平 导路 bf接着画速度多边形，由速度多边形求得：vf = v=0求4 4=3=va4/lo3a=0求va4a3v=400.025=1m/s 方向：如速度图所示b、加速度分析求aka4a3aka4a3=24va4a3=0 方向：如速度图所示求aa3aa3=aa2=22lo2a=8.382142/1000=9.97m/s2 方向：ao2 求ana4 ana4=23lo3a=0求aa4 ana4 + ata4 = aa3 + aka4a3 + ara4a3大小： 0 ？ 9.97 0 ?方向：ao3 o3a 如图 o3a取a=0.2 在1#图纸的左下方画加速度多边形大致图形如该图所示aa4= a=490.2=9.8m/s2求ab 用加速度影像求ab=1040.2=20.8m/s2 方向：如图所示 求af = + + 大小： ？ 20.8 0 ?方向：水平 fb bfanfb= v2bf/ lbf=0接着画加速度多边形,由加速度多边形得：af=a=1000.2=20m/s2 方向：水平向左 第12点：a、速度分析求va3va3 =va2=lo2a2= lo2an/30=(142/1000)803.14/30=0.1428.37=1.19m/s 求va4 = + 大小 ? 1.19 0方向：o3a o2a o3a取v=0.025 ,在1#图纸的左下方画速度多边形求vb用速度影像求vb=520.025=1.3m/s求vf= + 大小： ? 1.3 ?方向：水平 导路 bfvf=v=530.025=1.33m/s 方向：水平向左求44=3=va4/lo3a=210.025/（600.005）=1.75 rad/s 方向：逆时针求va4a3va4a3=v=330.025=0.825m/s 方向：如速度图所示b、加速度分析求aka4a3aka4a3 =24va4a3=21.750.825=2.89m/s2 方向：如速度图所示 求aa3aa3=aa2=22lo2a=8.372142/1000=9.97m/s2 方向：ao2求ana4ana4=23lo3a=1.752600.005=0.918m/s2 方向：ao3求aa4 + = + + 大小：0.918 ？ 9.97 2.89 ?方向：ao3 o3a 如图 o3a取a=0.20在1#图纸的左下方画加速度多边形，由加速度多边形求得：aa4= a=540.20=10.8m/s2求ab用加速度影像求ab=1420.2=28.4m/s2 方向：如图所示求af = + + 大小： ？ 28.4 0.775 ?方向：水平 fb bfanfb= v2bf/ lbf=(170.025) 2 /(233/1000)=0.775 m/s2接着画加速度多边形，由加速度多边形得：af= a=1360.20=27.2m/s2 方向：水平向左收集同组同学的位移、速度、加速度的数据并汇编在如下表：在1#图纸左上角绘制刨头的运动线图。刨头的运动线图见1#图纸。 曲柄位置名称 结果1234456sf00.040.110.200.280.33 0.41 vf01.171.171.561.71.61.43af14.58.906.22.25-0.09-1.16-4.20 曲柄位置名称 结果788910101112sf0.490.550.530.420.2750.210.05vf1.100-1.21-2.66-3.6-3.16-1.3af-8.2-16.3-26-28.2-1417.1327五、机构的动态静力分析已知:导杆的质量 m4=15kg 刨头的质量 m6=58kg (其余质量忽略不计) 导杆绕重心的转动惯量 js4=0.7kgm 切削阻力为常数大小为 fc=1300n1.确定惯性力、惯性力矩第5点：p16=-m6af=-581.16=-67.28npi4=-m4as=-151.35=-20.25nm14=-js44=-0.71.28=-0.896n mlh4 =m14/ pi4=0.044 m 第12点：p16=-m6af=-5827.2=-1577.6npi4=-m4as=-1514.1=-211.5nm14=-js44 =-0.736.54=-25.58nmlh4=m14/ pi4=0.121m将计算结果汇总在如下表中：曲柄位置导杆4刨头pi4m14lh4p165点20.250.90.04467.2812点211.525.58/0.1211577.62.确定齿轮2的重量查指导书得齿轮2的重量g2=500n3.确定各运动副反力 第5点:a、取构件5、6为示力体在机构位置图上方绘制示力体图比例尺为：l=0.005m/mm大致图形如图： + + + + =0 上式中只有fr45、fr76的大小未知取力比例尺：p=20n/mm 在机构位置图下面画力多边形 大致图形如图求得：fr45=p=6920= 1380n 方向与力多边形中的方向一致 fr76=p=3.4520=387.8n方向：垂直导路向上mf=0：fc（lo2m-yp）+g6xs6=fr76h76h76=fc(lo2m-yp)+g6xs6/r76 =0.611mb、取构件3、4为示力体在机构位置图右侧绘制示力体图比例尺为：l=0.005m/mm 大致图形如图其平衡方程为： + + + + = 0mo3=0 (确定fr23的大小)：fr23h23+f14hp+g4h4=fr54h54量得:hp =0.16m ；h4=0.015m； h54=0.745mfr23=(fr54h54+f14hp+g4h4)/h23=1682n 矢量式中fr74的大小和方向未知 仍取力比例尺p=20n/mm接着画力多边形图求得:fr74=p=2920=580n方向与力多边形中的方向一致 b、取构件2为示力体在机构位置图右下方绘示力体图比例尺为:l=0.005m/mm其平衡方程为： + + + = 0 mo2=0 （确定pb的大小）：fr32h32=pbrbpb=fr32h32/rb=640n式中的r72大小和方向未知仍然取力比例尺p=20n/mm接着画力多边形图，求得： fr72=p=10120=2020n 方向与为多边形中的方向一致 第12点:a、取构件5、6为示力体在机构位置图上方绘制示力体图 比例尺为：l=0.005m/mm + + + =0上式中只有r45、r76的大小未知取力比例尺：p=fc/=20n/mm在机构位置图下面画力多边形图求得：fr45=p=8920 =1780n 方向与力多边形中的方向一致 fr76=p=3220 =740n 方向：垂直导路向上 mf=0：g6xs6+=fr76h76h76=g6xs6+/fr76=0.433m b、取构件3、4为示力体在机构位置图右侧绘制示力体图比例尺为：l=0.005m/mm 大致图形如图其平衡方程为： + + + + = 0mo3=0 (确定r23的大小)：fr54h54+fi4hi4+g4h4=fr23l23量得:hp =0.505m； h4=0.11m； h54=0.72m fr23=(fr54h54+ p14hp+g4h4)/lo3a=4789n 矢量式中fr74的大小和方向未知 仍取力比例尺p=20n/mm接着画力多边形图,求得:fr74=p=6420=1280n 方向与力多边形中的方向一致b、取构件2为示力体在机构位置图右下方绘示力体图 比例尺为:l=0.005m/mm其平衡方程为： + + + = 0mo2=0 （确定pb的大小）：fr32h32=pbrb量得:h32=0.07m算得:rb =0.305mpb=fr32h32/rb=1099n仍然取力比例尺p=20n/mm接着画力多边形图,求得：fr72=p=26320=5260n方向与力多边形中的方向一致4、将各运动副反力汇总如下： 位置反力指定的两个位置第5点第12点fr7220205260fr745801280fr76387.8740fr4513801780fr3416824789 fr23168247895、计算平衡力偶矩并汇总如下：曲柄位置123456mb0172318342.3220.2195曲柄位置789101112mb11841.2158515-254-3136、绘制平衡力偶矩曲线mb-2该曲线在1#图纸的右上角 纵坐标比例尺：mb=10nm/mm 横坐标比例尺：2=2度/毫米平衡力偶矩曲线mb-2见1#图纸。六、飞轮设计已知：许用速度不均匀系数 =1/40 平衡力矩曲线 mb-2 驱动力矩为常数 曲柄的转数 n2=80rpm 飞轮装在齿轮z1的o1轴上1、作等效阻力矩曲线mr1-1由于飞轮准备装在z1的o1轴上，因此|mr|=|mb/i12|可由mb-2曲线直接画出mr11曲线（见1#图）。为了使图形一样，其比例尺选为：mr=mb/i12=10/2.7=3.8nm/mmi12=z2/ z1=46/18=2.62、求功曲线wr1-1取极距h=30mm图解积分mr1-1得wr1-1曲线。纵坐标比例尺为：w =mr1h/180=3.80.0930=10.26j /mm 3、求功曲线wd1-1根据一个稳定运转循环中能量变化为零，以及md=常数的条件可作出wd1-1曲线。比例尺仍为：w=10.26j/mm4、求驱动力矩曲线md1-1仍取极距h=30mm图解微分wd1-1得md1-1曲线。纵坐标比例尺为：mr=3.8nm/mm 得驱动力矩：md1=hmr=153.8=57nm 5、确定最大盈亏功将功曲线变成动能曲线。量取：w，=45 mm最大盈亏功为：w=4510.26=461.7j 6、求飞轮的转动惯量n1= n2i12=802.6=208 rpmjf=900w/2n12=900461.740/(22082)=38.9 kgm 7、确定飞轮尺寸b=4gjf/d3h材料用灰铸铁=7104n/m3取飞轮直径d=0.5m取轮缘的高宽比为h/b=1.5b2=4gjf/1.5d3=49.838.9/（3.141.50.537104)b=192mmh=1.5b=288mm图3七、设计凸轮轮廓曲线已知：推杆的运动规律为等加速等减速上升和等加速等减速下降,凸轮与曲柄共轴,顺时回转；凸轮机构的最大摆角 max=16凸轮的摆杆长 lo4c=140mm凸轮的推程运动角 0=60凸轮的远休止角 01=10凸轮的回程运动角 0=60凸轮机构的机架长 lo2o4=150mm凸轮的基圆半径 ro=55mm凸轮的滚子半径 rr=15mm在2#图纸上绘制摆杆的角位移曲线和凸轮轮廓曲线图。摆杆的角位移曲线和凸轮轮廓曲线图形见2#图纸。八、齿轮设计及绘制啮合图已知：齿轮1的尺数 z1=18齿轮2的尺数 z2=46模数 m12=15压力角 =20齿顶高系数 h*a=1径向间隙系数 c*=0.251、 列表计算几何尺寸等（表格在17页） 2、绘制齿廓啮合图在2#图纸上绘制齿廓啮合图。取比例尺为：l=0.001m/mm齿廓啮合图见2#图纸。九、解析法1导杆机构设计已知：（1）行程速比系数k；（2）刨头和行程h；（3）机架长lo2o3（4）连杆与导杆的比lbf/lo3b求解: (1)求导杆的摆角：max=180（k-1）/（k+1）(2)求导杆长：lo3b1=h/2sin（max/2）(3)求曲柄长：lo2a=lo2o3sin（max/2）(4)求连杆长：lbf=lo3blbf/lo3b(5)求导路中心到o3的垂直距离lo3m：从受力情况（有较大的传动角）出发，刨头导路o3b线常取为通过b1b2挠度de的中点m.即： lo3m=lo3b-lde/2将上述已知条件和公式编入程序。（源程序和运行结果见附录a）结果分析：与图解法比较，误差在毫米以下，不用修改。2机构的运动分析已知： (1)曲柄转速2； (2)各构件的长度。名称符号计算公式计算结果小齿轮分度圆直径d1d1=mz1270大齿轮分度圆直径d2d2=mz2690小齿轮齿顶圆直径da1da1=d1+2ha300大齿轮齿顶圆直径da2da2=d2+2ha720小齿轮齿根圆直径df1df1=d1-2hf232.6大齿轮齿根圆直径df2df2=d2-2hf652.6小齿轮基圆直径db1db1=d1cos253.8大齿轮基圆直径db2db2=d2cos648.4分度圆齿距pp=m47.1基圆齿距pbpb=pcos44.3分度圆齿厚ss=p/223.6分度圆齿槽宽ee=p/223.6径向间隙cc=c*m3.25标准中心距aa=m(z1+z2)/2480实际中心距aa=a480传动比ii=z2/z12.56重合度=b1b2/pb1.56求解：、建立机构的运动方程式如图所示：选定直角坐标系xoy。标出各杆的矢量和转角。各构件矢量所组成的封闭矢量方程式为: + = a b其中令：ll=lo2o3；y=l03m；s=l03a；将a式分别投影在x和y轴上得l2cosf2=s cos f4 cll+l2 sin f2=s sin f4 d两式相除则得tgf4=(ll+l2sinf2)l2cosf2 (1)在三角形a0203中s2=llll+l2l22l1l2cos(90+f2) (2)将c d两式对时间求导一次得l2w2sinf2=sw4sinf4+vrcosf4 el2w2cosf2=sw4cosf4+vrsinf4 f将坐标xoy绕o点转f4角(也就是将e f两式中的f2角f4角分别减去f4)，经整理后可分别得到vr=l2 w2sin(f2f4) (3)w4=l2 w2 cos(f2-f4)s (4) 再将e f二式方别对时同求导一次后，同样将坐标xoy绕0点转f4角(也就是将式中的f2角f4角分别成去f4)，经整理后可分别得到ar=sw4w4l2w2w2cos(f2f4) (5)ak=2 vr w4 (6)e4=2 vr w 4+ l2w2w2sin(f2一f4) (7)将b式分别投|影在x和y轴上得x：l4 cos f4十l5 cos f5 (8)y：l4 sin f4十l5 sin f5 (9)由(9)式可直接得sin f5=（yl4sinf4）l5 (10)对(9)式求导，一次可得l4w4cosf4=l5w5cosf5 于是由g式可得w5=(l4w4cosf4)l5cosf5 (11)对g式求导一次经整理可得e5=（l4e4cosf4+l4w4 w4sinf4+l5w5w5sinf5)l5cosf5 (12)(8)式中的x是坐标值，要想得到f点的位移xf应该是xf=xx0xf=l4 cos f4+l5 cos f5一(l4 cos f40+l5 cos f50) (13)式中f40 f50是导杆4处在左极限位置l时，导杆4和连杆5与坐标的正向夹角对（13）式求导一次可得：vf=l4w4sinf4l5 w5sinf5 (14)对（14）式求导一次可得：af=l4cosf4w4w4l4sinf4e4l5cosf5 w5w5l5sinf5e5 (15)角度的分析关于f4和f5两个角度的分析 当曲柄2运动到第一象限和第四象限时，导杆4在第一象限。此时得出的f4就是方位角。当曲柄2运动到第二象限和第三象限时导杆4是在第二象限，得出的f4是负值，所以方位角应该是f4=180+f4由于计算机中只有反正切，由（10）式是不能直接求出f5因此要将其再转换成反正切的形式f5=atn(gsqr(1g*g) (16) 式中g=sin f5=(yl4*sin f4) l5无论曲柄2运动到第几象限。连杆5都是在第二第三象限，由于在第二象限时f5是负值，在第三象限时f5是正值，因此在转换方位角时可以用一个公式来表示即：f5=180+f5 (17)开始计算是在左极限l的位置。因此f2的初值应该是： f2=fq=195(fq为起始角)运行到8时导杆处在右极限终止位置，因此f2的数值应该是：f2=fz=345 (fz为终止角)编写程序。（源程序和运行结果见附录b）结果分析： 上述结果与图解法比较，除加速度略有点误差外其余各结果均无误差。因此验证了图解法和解析法的运算结果都是正确的。加速度的误差尽管很小但也进行了查找修正。3机构动态静力分析已知: (1)各构件的质量； (2)导杆绕自身重心的曲转动惯量为js4； (3)切削阻力fc的变化规律； (4)齿轮2的重量g2=500 n。求解：一、建立直角坐标系(与运动分析时的坐标相一致如图所示) 二、求出刨头6和导杆4质心点的加速度和导杆4的角加速度。 三、确定刨头6和导杆4的惯性力和导杆4的惯性力矩。 四、远离原动件拆基本杆组建立机构的力方程式根据已知条件求出各运动副反力及加在原动件上的平衡力和平衡力矩。 （1）取构件56为示力体可得到方程式： (在回程和工作行程的两端0.05h处fc=0) 向x轴投影 r45cos(f45)+p16fc=0 向y轴投影r45sin（f45）+r76g6=0向f点取距r76h76=fc(l02my p)+g6*xs6（2）取构件34为示力体可得到方程式： 向o3点取矩求r23r23lo3+p14hp+g4h4=r54h54 （3）取构件2为示力体可得到方程式： 编写程序。（源程序及运行结果见附录c）结果分析： 上述结果与图解法比较，除某个运动副反力略有点误差外其余各结果均无误差。因此验证了图解法和解析法的运算结果都是正确的。力的误差尽管很小但也进行了查找修正。4凸轮机构的轮廓曲线设计已知(1)从动件8的运动规律及0、01、；(2)从动件8的长度lo4c；(3)从动件的最大摆角max=18；(4)从动件与凸轮的中心距lo2o4；(5)凸轮理论轮廓的基圆半径r0；(6)滚子半径rr；(7)凸轮与曲柄共轴,顺时针回转。1.建立数学模型选取xoy坐标系，b0点为凸轮起始点。开始时推杆滚子中心处于b0点处，当凸轮转过角度时，推杆相应地产生位移角。根据反转法原理,此时滚子中心应处于b点，其直角坐标为：y = a sinlsin(0) ax = a coslcos(0) b式中a为凸轮轴心o与摆动推杆轴心a之间的距离,l为摆动推杆的长度。在oa0b0中0=arc cos (a2l2r20)/2al ca式和b式即为凸轮理论轮廓线的方程式。凸轮的实际廓线与理论廓线的距离处处相等，为其理论廓线的等距曲线，且等于滚子半径rr，故当已知理论廓线上任意一点b(x，y)时，只要沿理论廓线在该点的法线方向取距离为rr，即得实际廓线上得相应点b(x,y)。由高等数学知，理论廓线b点处法线nn得斜率(与切线斜率互为负倒数)，应为：tgdx/dy(dx/d)/(dy/d) d式中dx/d、dy/d可根据a式和b式求得：dy/da coslcos(0)(1-d/d) edx/da sinlcos(0)(1-d/d) f代入d式可求出。此处应注意：角在0至360之间变化，当式中分子分母均大于0时，角在0至90之间；分子分母均小于0时，角在180至270之间；如果dy/d0则角在90至180之间；又如dy/d0，dx/d0，则角在270至360之间。当求出角后，实际廓线上对应b(x,y)的坐标可由下式求出：x=xrrcos gy=yrrcos h式中“”号为内等距曲线，“”号为外等距曲线。g式和h式即为凸轮的实际轮廓线方程式。在e式和f式中的即为给定的运动规律中的摆角推杆的角位移，d/d为对凸轮转角的导数，根据给定的等加速等减速上升和等加速等减速返回的运动规律其公式如下：1推程推杆在等加速运动阶段的方程式(在0至0/2之间变化)=2max(/0)2d/d=4max/02推杆在等减速运动阶段的方程式(在0/2至0之间变化)max2max(0)2/02d/d=4max(0)/022回程推杆在等加速运动阶段的方程式(在0至0/2之间变化)=2max(/0)2d/d=4max/02推杆在等减速运动阶段的方程式(在0/2至0之间变化)max2max(0)2/02d/d=4max(0)/02编写理论廓线的极坐标(向径)和理论廓线及实际廓线的直角坐标计算程序。（理论廓线的极坐标(向径)和理论廓线及实际廓线的直角坐标计算程序和运行结果见附录d）十、本设计的思想体会时光飞逝，转瞬间，劳累繁忙但是充实的课程设计就要结束了，回顾这十几天的过程，感到受益匪浅。通过这次设计,在同学的互相帮助，积极交流以及老师的悉心教导下，我对机械基础的基本原理以及广泛使用有了深刻理解，并在一定程度上可以洞悉大部分机械的工作原理 形成。真实的提高了在课堂上所学只是理论体系的内涵。并大大提高了我对这些知识的综合运用能力.如不仅用上了连杆机构,齿轮机构,凸轮机构,结构分析,运动分析,受力分析以及飞轮设计等,且知道如何把他们结合起来综合应用,深刻理解到只有全面的掌握了这些知识,将来才能真正的搞好运动设计。图解法直观明了，运动分析准确，计算机准确计算，只有有把两种方法都学好,都掌握了结合起来用才是最完善的。3张大图纸，一本厚厚的书，再加上老师的悉心教导，同学的热心帮助积极交流，我受益匪浅，这次的课程设计让我难以忘怀。参考文献1张国荣编著、机械原理课程设计、西北工业大学、2012.32东北工业大学、机械学教研室 编著、机械原理课程设计指导书、2008.33西北工业大学、孙恒，陈作模 主编、机械原理（第七版）、高等教育出版社、2006.54蒋加伏，沈岳 主编、计算机文化基础、北京邮电大学出版社、2003.55杨路明 主编、c语言程序设计教程（第二版）、北京邮电大学出版社、2005.12附录附录a导杆机构设计源程序及运行结果：#includemain()float pi,p,k,h,l1,bfo3b,qmax,l4,l2,l5,y;pi=3.1415926;p=pi/180;k=1.6;h=0.55;l1=0.4;bfo3b=0.3;qmax=180*(k-1)/(k+1)*p;l4=h/(2*sin(qmax/2);l2=l1*sin(qmax/2);l5=l4*bfo3b;y=l4-l4*(1-cos(qmax/2)/2;printf(l1=%.4fnl2=%.4fn,l1,l2);printf(l4=%.4fnl5=%.4fny =%.4fn,l4,l5,y);getch ();运行结果l1=0.4000l2=0.1418l4=0.7755l5=0.2327y =0.7503附录b机构运动分析源程序及运行结果：#includefloat pi,p,k,h,l1,bfo3b,qmax,l4,l2,l5,fq,y;float f2,f40,g,f50,n,f4,i,fz,w2,s,w4,v,ak,ar,e4;float l,f5,w5,e5,xf,vf,af,ls,at,an,as,fs;gan()pi=3.1415926;p=pi/180;k=1.6;h=0.55;l1=0.4;bfo3b=0.3;qmax=180*(k-1)/(k+1)*p;l4=h/(2*sin(qmax/2);l2=l1*sin(qmax/2);l5=l4*bfo3b;y=l4-l4*(1-cos(qmax/2)/2;fq=180+qmax/2/p;fz=360-qmax/2/p;main()gan();f2=fq*p;f40=atan(l1+l2*sin(f2)/l2/cos(f2);f40=180*p+f40;g=y/l5-l4/l5*sin(f40);f50=atan(-g/sqrt(1-g*g);f50=180*p+f50;printf( n f4 ak);printf( w4 sf vf af);printf( w4 as fsn);for(n=1;n=13;n=+n-0.5)if(n=4.5)elseif(n=10.5)f2=270*p;f4=90*p;elseif(n=8.5)f2=fz*p;f4=atan(l1+l2*sin(f2)/l2/cos(f2);elseif(n!=4.5&n!=8.5&n!=10.5)i=30*(n-1);f2=(fq-i)*p;if(f20)f2=360*p+f2;f4=atan(l1+l2*sin(f2)/l2/cos(f2);if(f2/p90)f4=180*p+f4;w2=-2*pi*80/60;s=sqrt(l2*l2+l1*l1+2*l1*l2*sin(f2);w4=l2*w2*cos(f2-f4)/s;v=-l2*w2*sin(f2-f4);ak=2*v*w4;ar=s*w4*w4-l2*w2*w2*cos(f2-f4);e4=-(ak+l2*w2*w2*sin(f2-f4)/s;l=l4/l5;g=y/l5-l*sin(f4);f5=atan(-g/sqrt(1