机械原理课程设计-牛头刨床的综合设计与分析

1、全套图纸加V信 sheji1120或扣 3346389411目 录一、设计题目与原始数据 1二、牛头刨床示意图 2三、导杆机构设计 2四、机构的运动分析 4五、机构的动态静力分析 9六、飞轮设计14七、凸轮设计15八、齿轮设计16九、解析法171.导杆机构设计172.机构的运动分析183.机构动态静力分析184.凸轮机构的轮廓曲线设计23十、本设计的思想体会 24十一、附录24十二、参考文献321、 设计题目与数据1题目牛头刨床的综合设计与分析2原始数据刨头的行程 H=600mm行程速比系数 K=1.8机架长 LO2O3=370mm质心与导杆的比值 LO3S4/LO3B=0.5连杆与导杆的比值

2、 LBF/LO3B=0.3刨头重心至F点距离 XS6=210mm导杆的质量 m4=20kg刨头的质量 m6=52kg导杆的转动惯量 JS4=0.9kgm切割阻力 FC=1400N切割阻力至O2的距离 YP=165mm构件2的转速 n2=80rpm许用速度不均匀系数 =1/30齿轮Z1、Z2的模数 m12=16小齿轮齿数 Z1=18大齿轮齿数 Z2=42凸轮机构的最大摆角 max=18º凸轮的摆杆长 LO4C=130mm凸轮的推程运动角 0=60º凸轮的远休止角 01=10º凸轮的回程运动角 0'=60º凸轮机构的机架长 Lo2o4=140mm凸轮

3、的基圆半径 ro=50mm凸轮的滚子半径 rr=15mm二、牛头刨床示意图图1三、导杆机构设计1、已知：行程速比系数 K=1.8刨头的行程 H=600mm机架长度 LO2O3=370mm连杆与导杆的比 LBF/LO3B=0.32、各杆尺寸设计如下A、求导杆的摆角：max =180°×（K-1）/（K+1）=180°×（1.8-1）/（1.8+1）=51.4° B、求导杆长：LO3B1=H/2sin（max/2）=550/2sin（51.4°/2）=691.8mm C、求曲柄长：LO2A =LO2O3×sin（max/2）=4

4、00×sin25.7°=160.5mm D、求连杆长LBF=LO3B×LBF/LO3B=691.8×0.3=160.5mm E、求导路中心到O3的距离LO3M =LO3B-LDE/2=LO3B1-1-cos(max/2)/2=655mm F、取比例尺 L=0.005m/mm 在A1图纸中央画机构位置图，大致图形如下： 图2四、机构的运动分析已知：曲柄转速n2=80rpm各构件的重心： 构件6的重心：XS6=210mm第1点：A、速度分析求VA3VA3 =VA2=LO2An/30=0.1605×80/30=1.34m/S = + 大小： ？ 1.

5、34 ？方向：O3A O2A O3A取v=VA3/Pa3=0.02 在A1图的左下方画速度多边形大致图形如A1图A 7所示。求VB用速度影像求VB=67*0.02=1.34m/s;求VF = + 大小： ？ 1.34 ？方向：水平 导路 BF VF =pf×V=0m/s 求4 4=3=VA4/LO3A=0rad/s 求VA4A3VA4A3 = ×V=1.34m/s B、加速度分析求aKA4A3aKA4A3=24VA4A3=0m/s2 求aA3A3=aA2=22×LO2A=11.2m/s2 方向：AO2 求anA4anA4=23×LO3A=0m/s2 方向

6、：AO3 求aA4 + = + + 大小 0 ？ 11.2 0 ?方向：AO3 O3A 如图6 O3A取a=aA3/=0.38 在A1图的左下方画加速度多边形，大致图形如B 7图所示，则：aA4=pa4×a=30*0.38=11.4m/ s2 ；求aB 用加速度影像求aB=61×0.38=23.18m/s2 方向如A1图 所示 ；求aF = + + 大小： ？ 23.18 4.14 ?方向：水平 FB BF接着画加速度多边形。由加速度多边形得：aF=×a=51*0.38=19.38m/s2 方向：水平向右 。 第4点：A、速度分析求VA3VA3 =VA2=LO2A

7、n/30=0.1605×80/30=1.34m/s 求VA4 = + 大小： ？ 1.34 ？方向：O3A O2A O3A取V=VA3/Pa3=0.02 在A1图的左下方画速度多边形求VB用速度影像求VB=165×0.02=1.73m/s求VF = + 大小： ？ 1.77 ？方向： 水平 导路 BF接着画速度多边形。由速度多边形求得：VF=V=2.1m/s 方向水平向右；求44=3=VA4/LO3A=2.53rad/S 求VA4A3VA4A3= ×V=0.42m/s 方向如速度图A1 左下A 4所示 B、加速度分析求aKA4A3aKA4A3=24VA4A3=1.

8、8m/s2 求aA3aA3=aA2=22×LO2A=11.2m/s2 方向：AO2 求anA4anA4=23×LO3A=3.28m/s2 方向：AO3 求aA4 + = + + 大小：3.28 ？ 11.2 1.8 ?方向：AO3 O3A O3A O3A取a=aA3/pa3=0.1 在A1图的左下方画加速度多边形。 大致图形如A左下B 4aA4=pa4×a=11.2m/s2 求aB方向如A1图C 4所示用加速度影像求aB=pb×a=54*0.1=5.4m/s2 = + + 大小： ？ 5.4 0 ?方向：水平 FB BF接着画加速度多边形。由加速度多边形

9、求得：aF=pf×a=22*0.1=2.2m/s2 水平向右。 第6点：A、速度分析求VA3VA3 =VA2=LO2An/30=0.1605×80/30=1.34m/S 求VA4 = + 大小 ? 1.34 0方向：O3A O2A O3A取V=VA3/Pa3=0.02 在A1图的左下方画速度多边形，大致图形如A 12图 所示。求VB用速度影像求VB=85*0.02=1.7m/s。求VF= + 大小： ? 1.7 ?方向：水平 导路 BFVF=V=82*0.02=1.64m/s 方向水平向左 ； 求44=3=VA4/LO3A=62*0.02/105*0.005=2.36rad

10、/S 方向：逆时针；求VA4A3VA4A3=×V=0.54m/s 方向如速度A1图11所示。 B、加速度分析求aKA4A3aKA4A3 =24VA4A3=2*2.36*28*0.02=2.64m/s2求aA3aA3=aA2=22×LO2A=11.25m/s2 方向：AO2求anA4anA4 =23×LO3A=2.92m/s2 方向：AO3求aA4 + = + + 大小：2.92 ？ 11.2 2.64 ?方向：AO3 O3A 如图 O3A取a=aA3/pa3 =0.15 在A1图的左下方画加速度多边形 大致图形如该B 12图 所示aA4= ×a=24*0

11、.15=3.6m/s2 求aB用加速度影像求aB=32*0.15=4.8m/s2 方向如A1图11所示求aF = + + 大小： ？ 4.8 0.5 ?方向：水平 FB BF接着画加速度多边形。由加速度多边形得：aF= ×a =22*0.15=3.3m/s2 方向：水平向右 在A1图纸左上角绘制刨头的运动线图。大致图形如A1图左上所示。收集同组同学的位移、速度、加速度得数据并混编在下表中：曲柄位置1234456SF00.030.110.210.2980.320. 43VF00.981.471.7251.771.761.56aF19.3810.45.9252.0-0.5-1.05-3.

12、6曲柄位置788910101112SF0.5170.620.60.5850.420.2950.220.045VF1.060.80-0.26-1.5-4.37-1.96-1.56aF-7.2-13.0-22.8-16.8-25.2-4.6-30.037.0五、机构的动态静力分析已知：导杆的质量： m4=20Kg； 刨头的质量： m6=52Kg； (其余质量忽略不计) 导杆绕重心的转动惯量： JS4=0.9gm 切削阻力为常数大小为：；FC=1400N1.确定惯性力、惯性力矩第4点：P16=-m6×aF=-52×(-22*0.1)=114.4NP14=-m4×as=-

13、20×54*0.1/2=-54NM14=-2.5Nmh =M14/F14=-2.5/-45=0.005m第6点：P16=-m6×aF=-52×22*0.15=-171.6N P14=-m4×aS=-20×32/2*0.15=-48N M14=-JS4×c =-3.5Nmh =M14/F14=-3.5/-48=0.073m将计算结果汇总在如下表中：曲柄位置导杆4刨头P14M14Lh4P164点-54-2.50.005114.46点-48-3.50.073-171.62.确定齿轮2的重量查指导书得齿轮2的重量G2=500N 3.确定各运动

14、副反力 第4点:A、取构件5、6为示力体在机构位置图上方绘制示力体图，如C 7图：比例尺为：L=0.005m/mm + + + + =0上式中只有FR45、FR76的大小未知。取力比例尺：P=Fc/=15N/mm； 在机构位置图右下方面画力多边形，大致图形如A1图F 7。如图求得：FR45=×P=101*15= 1515N方向与力多边形中的方向一致 ；FR76=×P=21*15=315N； 方向垂直导路向上。因为：MF=0；FC（LO2M-YP）+G6×XS6=FR76h76h76=Fc×(LO2M-YP)+G6×XS6/R76 =0.84m

15、B、取构件3、4为示力体在机构位置图右侧绘制示力体图，比例尺为：L=0.005m/mm，大致图形如A1图D7所示。其平衡方程为： + + + + = 0 MO3=0 (确定FR23的大小)：FR23h23+F14hp+G4h4=FR54h54量得：hI4 =38.5*0.005=0.1925m ；h4=6*0.005=0.03m； h54=135.5*0.005=0.6775m；FR23=(FR54h54+F'14hP+G4h4)/h23=1991.70N 矢量式中FR74的大小和方向未知，仍取力比例尺P=15N/mm。 接着画力多边形图求得：FR74=he×P=33.9&#

16、215;20=678N ； 方向与力多边形中he的方向一致。C、取构件2为示力体在机构位置图右方绘示力体图，比例尺为:L=0.005m/mm 大致图形如A1图E 7 其平衡方程为： + + + = 0MO2=0 （确定Pb的大小）：FR32h32=Pbrb量得:h32=0.005*30=0.15m 算得:rb =0.3158m Pb=FR32h32/rb=984.42N 式中的R72大小和方向未知仍然取力比例尺P=15N/mm 接着画力多边形图，求得： FR72=lm×P=160×15=2400N 方向与为多边形中ij的方向一致 第6点:A、取构件5、6为示力体 在机构位置

17、图上方绘制示力体图, 如A1图 C 12所示： 比例尺为：L=0.005m/mm + + + =0上式中只有R45、R76的大小未知取力比例尺：P=FI6/=25N/mm在机构位置图下面画力多边形图，大致图形如A1图 F 12 求得：FR45=×P=49×25 =1225N 方向与力多边形中的方向一致 FR76=ea×P=19×25 =475N 方向：垂直导路向上 MF=0：G6×XS6=FR76h76h76=G6×XS6/FR76=0.60m B、取构件3、4为示力体在机构位置图右侧绘制示力体图，比例尺为：L=0.005m/mm，大

18、致图形如A1 图D12其平衡方程为： + + + + = 0MO3=0 (确定FR23的大小)：FR54h54+FI4hI4+G4h4=FR23L23量得：hI4 =52x0.005=0.26； h4=12*0.005=0.06m； h54=135*0.005=0.68m FR23=(FR54h54+ P14hP+G4h4)/LO3A=1593.04N 矢量式中FR74的大小和方向未知，仍取力比例尺P=25N/mm；接着画力多边形图,求得:FR74=fg×P=8*25=200N 方向与力多边形中fg的方向一致C、取构件2为示力体在机构位置图右下方绘示力体图 比例尺为:L=0.005m

19、/mm大致图形如A1图E12：其平衡方程为： + + + = 0MO2=0 （确定Pb的大小）：FR32h32=Pbrb量得:h32=30*0.005=0.15m 算得:rb =63*0.005=0.315m Pb=FR32h32/rb=758.59N 仍然取力比例尺P=25N/mm 接着画力多边形图,求得：FR72=lm×P=62×25=1550N方向与为多边形中lm的方向一致。 4、将各运动副反力汇总如下： 位置反力指定的两个位置第4点第6点FR7224001550.00FR74678200.00FR76315475.00FR4515151225.00FR231991.

20、701593.04 FR321991.701593.045、计算平衡力偶矩并汇总如下：曲柄位置123456Mb0227.3320.1314.1271.7226.3曲柄位置789101112Mb132.6-16.868.75388.3-147.2-478.1 6、绘制平衡力偶矩曲线Mb-2该曲线在A1图的右上角 纵坐标比例尺：Mb=10Nm/mm 横坐标比例尺：2=2度/毫米 见A1图 I六、飞轮设计已知：许用速度不均匀系数 =1/40 平衡力矩曲线 Mb-2 驱动力矩为常数 曲柄的转数 n2=80rpm 飞轮装在齿轮Z1的O1轴上1、作等效阻力矩曲线Mr-r由于飞轮准备装在Z1的O1轴上，因此

21、|Mr|=|Mb/i12|可由Mb-2曲线直接画出Mr-1曲线（见A1图I）。为了使图形一样，其比例尺选为：Mr=Mb/i12=4.3Nm/mm 。2、求功曲线Wr-1取极距H=40mm； 图解积分Mr-1得Wr-1曲线。纵坐标比例尺为：W =M××H×/180°=3.91×5.1×40×/180°=14.0J/mm 3、求功曲线Wd-1根据一个稳定运转循环中能量变化为零，以及Md=常数的条件，可作出Wd-1曲线。比例尺仍为：W=14.0J/mm。4、求驱动力矩曲线Md-1仍取极距H=40mm 图解微分Wd-1得M

22、d-1曲线。纵坐标比例尺为：M=3.91Nm/mm 得驱动力矩：Md =h×M=10×3.91=39.1Nm 5、确定最大盈亏功为：W=32×14.0=522J 6、求飞轮的转动惯量JF=900W/2n2=(900×522)/(3.142 ×802)=22.3Kgm 7、确定飞轮尺寸b=4gJF/D3H，材料用灰铸铁=7×104N/m3，取飞轮直径D=0.5m， 取轮缘的高宽比为H/b=1.5； b2=4gJF/D3H=(4×9.8×22.3)/(3.14×1.5×0.53×7

23、5;104)b=0.146m=146mm H=1.5b=1.5×146=219mm 图3七、设计凸轮轮廓曲线已知：推杆的运动规律为等加速等减速上升和等加速等减速下降，凸轮与曲柄共轴，顺时凸轮机构的最大摆角 max=18°；凸轮的摆杆长 LO4C=130mm；凸轮的推程运动角 0=60°；凸轮的远休止角 01=10°；凸轮的回程运动角 0'=60°；凸轮机构的机架长 Lo2o4=140mm；凸轮的基圆半径 ro=50mm；凸轮的滚子半径 rr=15mm；摆杆的角位移曲线以及凸轮轮廓曲线的设计已绘制在A2图纸上。八、齿轮设计及绘制啮合图已知

24、：齿轮1的尺数Z1=18；齿轮2的尺数Z2=42；模数m12=16；压力角=20°；齿顶高系数h*a=1；径向间隙系数C*=0.25。1、 列表计算齿轮尺寸如下：名称符号计算公式计算结果小齿轮分度圆直径d1d1=mz1288大齿轮分度圆直径d2d2=mz2672小齿轮齿顶圆直径da1da1=d1+2ha320大齿轮齿顶圆直径da2da2=d2+2ha704小齿轮齿根圆直径df1df1=d1-2hf248大齿轮齿根圆直径df2df2=d2-2hf632小齿轮基圆直径db1db1=d1cos270.63大齿轮基圆直径db2db2=d2cos631.48分度圆齿距PP=m50.26基圆齿距

25、pbpb=pcos47.23分度圆齿厚ss=p/225.13分度圆齿槽宽ee=p/225.13径向间隙cc=c*m4标准中心距aa=m(z1+z2)/2480实际中心距aa=a480传动比ii=z2/z12.33重合度=B1B2/Pb1.33九、解析法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=LO2O3×sin（max/2）；(4)求连

26、杆长：LBF=LO3B×LBF/LO3B； (5)求导路中心到O3的垂直距离LO3M：从受力情况（有较大的传动角）出发，刨头导路O3B线常取为通过B1B2 挠度DE的中点M。即： LO3M=LO3B-LDE/2将上述已知条件和公式编入程序（见附录）。与图解法比较，误差在毫米以下，不用修改。2机构运动分析已知： (1)曲柄转速2； (2)各构件的长度。求解：、建立机构的运动方程式如图所示：选定直角坐标系XOY，标出各杆的矢量和转角。各构件矢量所组成的封闭矢量方程式为： 图22+ = a b其中令：Ll=LO2O3；Y=L03M；S=L03A；将a式分别投影在x和y轴上得：L2

27、cosF2=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二式方别对时间求导一次后

28、，同样将坐标XOY绕O点转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式求导一次经整理可

29、得：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运动到

30、第一象限和第四象限时，导杆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=1

31、95°(Fq为起始角)，运行到8时导杆处在右极限终止位置，因此F2的数值应该是：F2=FZ=345°(FZ为终止角)。编写程序及运行结果见附录。结果分析： 上述结果与图解法比较，除加速度略有点误差外其余各结果均无误差。因此验证了图解法和解析法的运算结果都是正确的。加速度的误差尽管很小但也进行了查找修正3机构动态静力分析已知：(1)各构件的质量； (2)导杆绕自身t心曲转动惯量为J。； (3)切削阻力F的变化规律； (4)齿轮2的重量G=500 N。求解： 一、建立直角坐标系(与运动分析时的坐标相一致如图23所示) 二、求出刨头6和导杆4质心点的加速度和导杆4的角加速度。 三

32、、确定刨头6和导杆4的惯性力和导杆4的惯性力矩。 四、远离原动件拆基本杆组建立机构的力方程式根据已知条件求出各运动副反力及加在原动件上的平衡力和平衡力矩。 1取构件56为示力体可得到方程式： (在回程和工作行程的两端0.05H处FC=0) 向X轴投影 R45×cos(f45)+P16FC=0 向Y轴投影R45×sin（f45）+R76G6=0向F点取距R76×h76=FC·(L02mY p)+G6*Xs62取构件34为示力体可得到方程式： 向O3点取矩求R23R23LO3+P14hP+G4h4=R54h543取构件2为示力体可得到方程式： 源程序及运行结

33、果见附录结果分析： 上述结果与图解法比较，除某个运动副反力略有点误差外其余各结果均无误差。因此验证了图解法和解析法的运算结果都是正确的。力的误差尽管很小但也进行了查找修正.4凸轮机构的轮廓曲线设计已知(1)从动件8的运动规律及0、01、；(2)从动件8的长度LO4C；(3)从动件的最大摆角max=16º；(4)从动件与凸轮的中心距LO2O4；(5)凸轮理论轮廓的基圆半径r0；(6)滚子半径rr；(7)凸轮与曲柄共轴,顺时针回转。1.建立数学模型选取XOY坐标系，B0点为凸轮起始点。开始时推杆滚子中心处于B0点处，当凸轮转过角度时，推杆相应地产生位移角。根据反转法原理,此时滚子中心应处

34、于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、

35、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/d<0，dx/d>0则角在90°至180°之间；又如dy/d>0，dx/d<0，则角在270°至360°之间。当求出角后，实际廓线上对应B(x,y)的坐标可由下式求出：x=x±

36、;rrcos gy=y±rrcos 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推杆在等减速运

37、动阶段的方程式(在0/2至0之间变化)max2max(0)2/02d/d=4max(0)/02理论廓线的极坐标(向径)和理论廓线及实际廓线的直角坐标计算程序和运行结果见附录。十、设计体会附录1导杆机构设计程序及运行结果Private Sub Command1_Click()Dim Qmax%, K!, pi!, rad!, H!, L4!, L2!, L1!, bfo3b!, Y!pi = 3.1415926K = InputBox("请输入行程速比系数K")H = InputBox("请输入行程H")L1 = InputBox("请输入机架长

38、L1")bfo3b = InputBox("请输入连杆与导杆之比LBF/Lo3B")Qmax = 180 * (K - 1) / (K + 1)Label1.Caption = "导杆摆角：" & Qmaxrad = (pi / 180) * QmaxL4 = H / (2 * Sin(rad / 2)Label2.Caption = "导杆Lo3b1=" & L4L2 = L1 * in(rad / 2)Label3.Caption = "曲柄长Lo2a=" & L2L5 = L

39、4 * bfo3bLabel4.Caption = "连杆Lbf=" & L5Y = L4 - L4 * (1 - Cos(rad / 2) / 2Label5.Caption = "导路中心线XX到O3点的垂直距离Lo3m=" & YEnd SubPrivate Sub Form_Load()End Sub运行结果：导杆摆角：51°导杆LO3b1=0.6968462曲柄长：Lo2A=0.15928连杆Lbf=0.2090导路中心线xx到O3点的垂直距离LO3m=0.66292机构运动分析源程序及运行结果Private Sub

40、Form_Click()Dim Qmax!, K!, pi!, p!, H!, L4!, L2!, L1!, bfo3b!, Y!, Fq!, Fz!, F4!, F40!, F50!, L!Dim F2!, g!, L5!, ak!, w4!, xf!, vf!, af!, e4!, fs!, n!, a3!, i%, w2!, s!, ar!, F5!, w5!, e5!, at!, an!pi = 3.1415926: p = pi / 180K = InputBox("请输入行程速比系数K")H = InputBox("请输入行程H")L1 =

41、InputBox("请输入机架长L1")bfo3b = InputBox("请输入连杆与导杆之比LBF/Lo3B")Qmax = 180 * (K - 1) / (K + 1) * pL4 = H / (2 * Sin(Qmax / 2)L2 = L1 * Sin(Qmax / 2)L5 = L4 * bfo3bY = L4 - L4 * (1 - Cos(Qmax / 2) / 2Fq = 180 + Qmax / 2 / pFz = 360 - Qmax / 2 / pF40 = Atn(L1 + L2 * Sin(195 * p) / L2 / Cos(195 * p)F40 = 180 * p + F40F50 = 177 * pPrint Spc(0.5); "