曲柄滑块机构的运动律.docVIP

黔南民族师范学院 数学实验报告 题目 曲柄滑块机构的运动规律 制作成员班级学号任务吴佳丽B11数应（2）1105015241设计实验目的，实验问题张玉萍1105015211上机实验结果的分析与结论姜光云1105015242建立模型，设计程序 2013年4月22日 黔南民族师范学院数学实验报告 题目 PAGE VIII PAGE VII 数学实验报告题目 实验目的 1. 本实验主要涉及微积分中对函数特性的研究，通过实验复习函数的求导，Taylor公式和其他有关知识，着重介绍运用建立近似模型并进行数值计算来研究函数的方法。 实验问题 1. 给定一机构如图1.2所示，设连杆QP长度L=300mm ，曲柄OQ的长为r=100mm，距离e=20mm，曲柄的角速度w=240转/min。 对在一个周期（即【0.2】）中 1、计算滑块的位移； 2、计算滑块的行程； 3、计算滑块的速度； 4、计算滑块的加速度； 5、计算滑块的摆角及其最值. 建立数学模型 取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(l^2-(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(l^2-(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; a=w^2*((-r*cos(t)-(((r*r*cos(2*t)+r*e*sin(t))*(l^2-(r*sin(t)-e).^2)+r*r*cos(t).^2*(r*sin(t)-e).^2)/sprt((l^2-(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 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 Columns 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/