《曲柄滑块机构的运动规律.doc》由会员分享，可在线阅读，更多相关《曲柄滑块机构的运动规律.doc（12页珍藏版）》请在金锄头文库上搜索。
1、黔南民族师范学院数学实验报告题目曲柄滑块机构的运动规律制作成员班级学号任务吴佳丽B11数应（2）1105015241设计实验目的，实验问题张玉萍1105015211上机实验结果的分析与结论姜光云1105015242建立模型，设计程序2013年4月22日题目数学实验报告题目一、 实验目的 1. 本实验主要涉及微积分中对函数特性的研究，通过实验复习函数的求导，Taylor公式和其他有关知识，着重介绍运用建立近似模型并进行数值计算来研究函数的方法。 二、 实验问题 1. 给定一机构如图1.2所示，设连杆QP长度L=300mm ，曲柄OQ的长为r=100mm，距离e=20mm，曲柄的角速度w=240转
2、/min。 对在一个周期（即【0.2】）中 1、计算滑块的位移；2、计算滑块的行程；3、计算滑块的速度；4、计算滑块的加速度；5、计算滑块的摆角及其最值.三、 建立数学模型1. 取O点为坐标原点，沿x轴向右方向为正，P在x轴上的坐标为x,用x表示滑块的位移，利用三角关系有： 由于，故有 而于是滑块的速度 (1.4)从而，得到滑块的加速度为 由关系式得摆角的表达式为四、 问题求解和程序设计流程利用滑块位移的解析式（1.1），可以应用数学软件MATLAB进行计算，编制名为m1_1.m的m文件：function m1_1(t)r=100;l=300;e=20;x=r*cos(t)+sqrt(l2-(
3、r*sin(t)-e).2);end然后再命令窗口键入m1_1(0:pi/12:2*pi) 滑块的加速度在0,内变化时的滑块速度 编制MATLAB的m文件m1_2.m:function m1_2(t)r=100;l=300;w=240/60*2*pi;e=20;v=-w*r*sin(t)+(w*r*cos(t)*(r*sin(t)-e)/(sqrt(l2-(r*sin(t)-e).2) end end然后在命令窗口键入m1_2(0:pi/12:pi)滑块的加速度：编制MATLAB的m文件m1_3.m:function m1_3(t)r=100;l=300;w=240/60*2*pi;e=20;
4、 a=w2*(-r*cos(t)-(r*r*cos(2*t)+r*e*sin(t)*(l2-(r*sin(t)-e).2)+r*r*cos(t).2*(r*sin(t)-e).2)/sprt(l2-(r*sin(t)-e).2).3) end然后在命令窗口键入format long g%数字显示方式为长格式（15位有效数字）m1_3(0:pi/12:pi)摆角的加速度及其最值：编制MATLAB的m文件m1_4.m:function m1_4(t)r=100;l=300;w=240/60*2*pi;e=20;b=arcsin(r*sin(t)-e)./l)后在命令窗口键入format long
5、g%数字显示方式为长格式（15位有效数字） m1_1(0:pi/6:pi) 五、 上机实验结果的分析与结论m1_1(0:pi/12:2*pi)x = Columns 1 through 10 399.3326 396.5349 385.0988 366.3937 342.5134 315.9398 289.1366 264.1760 242.5134 224.9723 Columns 11 through 20 211.8937 203.3498 199.3326 199.8781 205.1165 215.2466 230.4209 250.5348 274.9545 302.2986 Co
6、lumns 21 through 25 330.4209 356.6680 378.3216 393.0632 399.3326方便观察 ， 我们列出了从0变化到时位移的一些相应数据：/radx/mm0399.3326/12396.5342/12385.09883/12366.39374/12342.51345/12315.9398/2 289.13667/12264.17608/12242.51349/12224.972310/12211.893711/12203.3498199.332613/12199.878114/12205.116515/12215.246616/12230.4209
7、17/12250.534818/12274.954519/12302.298620/12 330.420921/12356.668022/12 378.321623/12393.06322399.3326速度：v = Columns 1 through 6 167.925190836271 -698.08944422326 -1475.38949547774 -2081.9412336388 -2462.68361477916 -2599.40275042012 Columns 7 through 12 -2513.27412287183 -2255.8700172317 -1890.4348
8、5938297 -1472.36511688789 -1037.8846273941 -602.876972901868 Columns 13 through 18 -167.925190836272 274.78133763591 734.356633416719 1213.40895844524 1698.84631194499 2153.26248260396 Columns 19 through 24 2513.27412287183 2702.01028504786 2654.27216221713 2340.89739208146 1778.91748945512 1026.185
9、07948922 Column 25 167.925190836272加速度：a = Columns 1 through 6 -84361.7758140926 -80348.6233660269 -67559.9492334591 -48000.0272284981 -24805.7860760063 -1797.80226341147 Columns 7 through 12 17476.9871791866 30899.0500454472 38359.6820909656 41329.4345268796 41845.8509156108 41677.6906982094 Column
10、s 13 through 18 41969.1605198511 43188.5486127271 45106.6956437737 46677.2239721538 45933.0244155622 40280.6743815471 Columns 19 through 24 27567.6703983307 7583.82207268844 -17232.4437514097 -42652.2377832238 -64299.1045052961 -78837.7654515092 Column 25 -84361.7758140926摆角的结果:b= Columns 1 through
11、6 -0.0667161484102253 0.0196076047273053 0.10016742116156 0.169851100110178 0.223873861713631 0.258166876496835 Columns 7 through 12 0.269932795833403 0.258166876496835 0.223873861713631 0.169851100110178 0.10016742116156 0.0196076047273054 Columns 13 through 18 -0.0667161484102252 -0.15354226997510
12、7 -0.2355042367208 -0.307176937338054 -0.36327970779016 -0.399157208905875 Columns 19 through 24 -0.411516846067488 -0.399157208905875 -0.36327970779016 -0.307176937338054 -0.2355042367208 -0.153542269975107 Column 25 -0.0667161484102253摆角的最大值为：max(b)=0.269932795833403; 摆角的最小值为：min(b)=-0.399157208905875; 六、 实验总结与体会通过本次实验，了解了关于曲柄滑块的运动规律，并熟悉了matlab的上机操作，让自己学会了处理模型的方法及步骤。 说明：
