上海应用技术学院
数学实验报告
题目： 曲柄滑块机构的运动规律
姓名： 周 玲
院系： 理学院数学与应用数学系
学号： 1112211115
指导老师： 许建强
2015年3月30日
目录
一、 实验目的........................................................................................... 3 二、 实际问题........................................................................................... 3 三、 数学模型........................................................................................... 3 四、 数值积分方法 .................................................................................. 2 五、 实验任务........................................................................................... 4
任务一 . ................................................................................................ 4 任务二 . ................................................................................................ 5 任务三 . ................................................................................................ 7 任务四 . ............................................................. 错误！未定义书签。
一、 实验目的
本实验主要涉及微积分中对函数特性的研究。通过实验复习函数求导法, Taylor 公式和其他有关知识。着重介绍运用建立近视似模型并进行数值计算来研究讨论函数的方法。
二、 实际问题
曲柄滑块机构是一种常用的机械结构, 它将曲柄的转动转化为滑块在直线上的往复运动, 是气压机、冲床、活塞式水泵等机械的主机构。右图为其示意图。
记曲柄OQ 的长为r , 连杆QP 的长为l , 当曲柄绕固定点O 以角速度w 旋转时, 由连杆带动滑块P 在水平槽内做往复直线运动。假设初始时刻曲柄的端点Q 位于水平线段OP 上, 曲柄从初始位置起转动的角度为θ, 而连杆QP 与OP 的锐夹角为β(称为摆角) 。在机械设计中要研究滑块的运动规律和摆角的变化规律, 确切的说, 要研究滑块的位移, 速度和加速度关于θ的函数关系, 摆角β及其角速度和角加速度关于θ的函数关系, 进而
（1）求出滑块的行程s (即滑块往复运动时左、右极限位置间的距离);
（2）求出滑块的最大和最小加速度(绝对值), 以了解滑块在水平方向上的作用力;
（3）求出β的最大和最小角加速度(绝对值), 以了解连杆的转动惯量对滑块的影响;
在求解上述问题时, 我们假定:
r =100(mm ), l =3r =300(mm ), ω=240(转/min)
符号说明:r －曲柄OQ 的长; l －连杆PQ 的长度; β－摆角(连杆PQ 与OP 的锐夹角); ω－角速度; P －滑块; x －滑块的位移; a －滑块的加速度。
三、 数学模型
取O 点为坐标原点，OP 方向为x 轴正方向，P 在x 轴上的坐标为x ，那么可用x 表示滑块的位移。利用三角关系，立即得到
x =r cos θ+l 2-r 2sin 2θ （1.1）
由于θ=ωt ，故有
dx dx d θdx ==ω （1.2） dt d θdt d θ
而
dx r 2sin θcos θ （1.3） =-r sin θ-
222d θl -r sin θ
于是滑块的速度
⎛r cos θ v =-ωr sin 1+
l 2-r 2sin 2θ⎝
进而，可以得到滑块的加速度为
a =
⎫
⎪ （1.4） ⎪⎭
dv dv =ωdt d θ
⎡⎤ （1.5） 224
r (l cos(2θ) +r sin θ) ⎥
=-ω2r ⎢cos θ+3⎢⎥
2222⎢⎥(l -r sin θ) ⎣⎦
同样，基于关系式
l sin β=r sin θ （1.6）
我们有摆角的表达式
β=arcsin sin θ⎪ （1.7）
式（1.6）对t 求导， l cos β可得
由此再得
⎛r
⎝l ⎫⎭
d βd θ=r cos θ=r ωcos θ dt dt
d βr ωcos θ
= （1.8） dt l cos β
d β⎫⎛ωsin θcos β-cos θsin β ⎪d 2βr ω dt ⎪=- （1.9） 22
l ⎪dt cos β ⎪⎝⎭
利用（1.6），不难由上两式导出
d βr ωcos θ= （1.10）
222dt l -r sin θ和
d 2βr ω2sin θ(l 2-r 2) =-3 （1.11） dt 2
(l 2-r 2sin 2θ) 2
至此，我们得到了滑块位移x 和连杆摆角β运动规律中有关变量依赖θ的表达式。
四、 数值积分方法
将位移的表达式（1.1）改写为
r 2
x =r cos θ+l (1-2sin 2θ) 2
l
r 2
一般而言，2是远比1小的数，于是利用
l
(1+ε) α=1+αε+ , 得到滑块位移的近似模型为
1
ε
r 2
x 1=r cos θ+l -sin 2θ （1.13）
2l
从而有相应的近似速度
⎛⎫dx 1dx 1d θr 2
⎪v 1===ω -r sin θ-sin(2θ) ⎪dt d θdt 2l ⎝⎭
（1.14）
r ⎛⎫
=-ωr sin θ-sin(2θ) ⎪
2l ⎝⎭
和近似加速度
dv 1r ⎛⎫
=-ω2r cos θ+cos(2θ) ⎪ （1.15） dt l ⎝⎭
这里的速度v 1和加速度a 1是直接对近似位移模型x 1求导得来的，而不是对v 和a 的精确表达式（1.4）和（1.5）的近似。当然我们也可以直接从滑块速度的解析式（1.4）进行近似。仍利用公式（1.12）有
a 1=
⎫1⎛r 2
= 1-2sin θ⎪⎪222
l -r sin θl ⎝l ⎭
1
2
-
1
2
⎫1⎛r 22
≈ 1+2sin θ⎪⎪ l ⎝2l ⎭
把上式代入（1.4），就得到滑块速度的近似模型
⎡r cos θ⎛⎫⎤r 22
v 2=-ωr sin θ⎢1+1+2sin θ⎪ ⎪⎥l 2l ⎝⎭⎦⎣
（1.16）
32
⎛r sin(2θ) r sin θsin(2θ) ⎫
⎪=-ωr sin θ++3 ⎪2l 4l ⎝⎭
从（1.16）出发，又可得近似加速度
⎡r cos(2θ) r 3(sin2(2θ) +2sin 2θcos(2θ) ⎤
a 2=-ωr ⎢cos θ++⎥ （1.17） 3
l 4l ⎣⎦
对摆角β可以利用幂级数展开的Maclaurin 公式
2
arcsin ε=ε+
ε3
6
+ , ε
得到摆角的近似模型。粗略一些，可以取
r
β1=sin θ （1.19）
l
r
（当较小时可用此式）。而必要时可以取
l
r r 3
β2=sin θ+3sin 3θ （1.20）
l 6l
相应的近似角速度为
d β1r
=ωcos θ （1.21） dt l
或
⎛r ⎫d β2r 32 =ω cos θ+3sin θcos θ⎪ （1.22） ⎪dt l 2l ⎝⎭
近似角加速度为
d 2β12r =-ωsin θ （1.23） 2
l dt
或
⎤d 2β2r 32⎡r 3
=-ω⎢sin θ+3(sinθ-sin (2θ) cos θ) ⎥ （1.24） dt 22l ⎣l ⎦
五、 实验任务
任务一
试用摆角的角加速度的三种表达式, 即式（1.11）、（1.23）和（1.24），取步长为
π
，r ，l ，ω的值如前，计算当θ∈[0,π]变化时角加速度的值, 并列表12
加以比较。
实验程序：
function m1_2(t)
r=100;l=300;w=240/60*2*pi;
a0=-r*w^2*sin(t)*(l^2-r^2)./(l^2-r^2*sin(t).^2).^(3/2) a1=-w^2*r*sin(t)/l
a2=-w^2*(r*sin(t)/l+r^3*(sin(t).^3-sin(2*t).*cos(t))/(2*l^3))
>> m1_1([0:pi/12:pi])
从结果中可以看出误差的大小, 取决于近似表达式的精度, 在利用泰勒公式求
d 2β2
近似模型时, 如果展开的精度越高, 则误差就越小, 在数据表中也可以看出,
dt
d 2β1
取得精度比高, 所以结果与真实值相差的更小。
dt
任务二
利用（1.12）式，对角摆角的角速度（1.10）式和角速度（1.11）式进行简化，将结果与（1.21）~(1.24)式进行比较，并与上题的计算结果相比较。 解析： 由式 (1+ε) α=1+αε+ , ε
d βr ωcos θ
(1.10) =
222dt l -r sin θ
⎫d β3r ω⎛r 2
⎪可以化简为 (1. 10) ' =cos θ+sin(2θ) sin θ2 ⎪dt l ⎝4l ⎭
d 2β3r ω2sin θ(l 2-r 2)
(1.11) =-3
dt 2
(l 2-r 2sin 2θ) 2
可以化简为
d 2βr ω2(l 2-r 2) ⎛3r 23⎫ (1. 11) ' =-sin θ+sin θ⎪232 ⎪dt l 2l ⎝⎭
将化简结果与（1.21）~(1.24)式进行比较，可以发现有类似的项。
实验程序：
function m2_1(t)
r=100;l=300;w=240/60*2*pi;
b0=r*w*cos(t)./sqrt(l^2-r^2*sin(t).^2) b1=w*r*cos(t)/l
b2=w*(r*cos(t)/l+r^3*sin(t).^2.*cos(t)/(2*l^3)) b3=r*w/l*(cos(t)+r^2*sin(2*t).*sin(t)/(4*l^2))
a0=-r*w^2*sin(t)*(l^2-r^2)./(l^2-r^2*sin(t).^2).^(3/2) a1=-w^2*r*sin(t)/l
a2=-w^2*(r*sin(t)/l+r^3*(sin(t).^3-sin(2*t).*cos(t))/(2*l^3))
a3=-r*w^2*(l^2-r^2)/l^3*(sin(t)+3*r^2*sin(t).^3/(2*l^2)) >> m2_1([0:pi/12:pi])
运行结果如下：
d 2β3d β3d βd β2d β
从上表中可看出与最接近真值，最接近真值3。由此看来，
dt dt dt dt dt
方案三最优。
任务三
给定一机构如右图所示。设连杆QP 长度l=300mm，曲柄OQ 的长为r=100mm，距离e=20mm，曲柄的角速度w=240转/min。对θ在一个周期（即[0, 2π]）中计算滑块的位移、行程、速度、加速度和摆角及其最值。
解析：这个机构的特点是:滑块的运动轨迹仍然在原来的平面上, 且与轴线Ox 平行, 但运动轨迹与Ox 有距离e (称为偏心距) 。这样进程时间将与退程时间不同。
由于P 点始终在直线y =e 上, 所以我们只需要考虑滑块在x 方向上的位移, 不需要再考虑在y 轴上的位移。取O 点为坐标原点，沿x 轴向右方向为正，P 在x 轴上的坐标为x, 用x 表示滑块的位移，利用三角关系有：
x =r cos θ+2-(r sin θ-e ) 2 （2.1）
由于θ=ωt ，故有
dx dx d θdx ==ω （2.2） dt d θdt d θ
dx
求导程序： d θ
>> syms r l e t
>> x=r*cos(t)+sqrt(l^2-(r*sin(t)-e).^2); >> diff(x,t)
dx r cos θ(e -r sin θ)
得 =-r sin θ+
22d θl -(e -r sin θ)
即
dx r cos θ(r sin θ-e )
（2.3） =-r sin θ-
22d θl -(r sin θ-e )
于是滑块的速度 v =-ωr sin θ-
从而，得到滑块的加速度为
a =
ωr cos θ(r sin θ-e )
l -(r sin θ-e )
2
2
(2.4)
dv dv =ωdt d θ
(2.5)
222222⎫⎛(r cos 2θ+re sin θ)(l -(r sin θ-e ) ) +r cos θ(r sin θ-e ) ⎪
=ω2 -r cos θ-
⎪(l 2-(r sin θ-e ) 2) 3⎝⎭
由关系式
l sin β+e =r sin θ (2.6）
得摆角的表达式为
β=arcsin(
r sin θ-e
) (2.7) l
滑块的行程：s =x max -x min =(l +r ) 2-e 2-(l -r ) 2-e 2=200.[1**********]
实验程序： （1）滑块的位移： function m3_1(t) r=100; l=300; e=20;
x=r*cos(t)+sqrt(l^2-(r*sin(t)-e).^2) >>m3_1([0:pi/12:2*pi])
（2）滑块的速度： function m3_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) >>m3_2([0:pi/12:2*pi])
（3）滑块的加速度： function m3_3(t) r=100; l=300;
w=240/60*2*pi; e=20;
a=w^2*(-r*cos(t)-((r^2*cos(2*t)+r*e*sin(t)).*(l^2-(r*sin(t)-e).^2)+r^2*cos(t).^2.*(r*sin(t)-e).^2)./(l^2-(r*sin(t)-e).^2).^(3/2)) >> m3_3([0:pi/12:2*pi])
（4）摆角及其最值： function m3_4(t) r=100; l=300;
w=240/60*2*pi; e=20;
b=asin((r*sin(t)-e)/l) >> m3_4([0:pi/12:2*pi])
运行结果如下：
3
由上表可看出，摆角β，在π处取得最小值0.0667；在π处取得最大值
2
0.4115。
任务四
设T1和T2分别为滑块进程和退程所需时间，试根据任务二中的数据求出T1和T2，两者是否相等？在l 和r 不变的前提下，令k=T1/T2。如果要求k=1.2，试求偏心距e 的值。
答：T 1>T 2。
由 0. 2π+2. 2arcsin
e e -2. 2arcsin =0， 200400
实验程序及结果：
>> x=fzero(@(x) 0.2*pi+2.2*asin(x./200),2) x =
-56.3465
上海应用技术学院
数学实验报告
题目： 曲柄滑块机构的运动规律
姓名： 周 玲
院系： 理学院数学与应用数学系
学号： 1112211115
指导老师： 许建强
2015年3月30日
目录
一、 实验目的........................................................................................... 3 二、 实际问题........................................................................................... 3 三、 数学模型........................................................................................... 3 四、 数值积分方法 .................................................................................. 2 五、 实验任务........................................................................................... 4
任务一 . ................................................................................................ 4 任务二 . ................................................................................................ 5 任务三 . ................................................................................................ 7 任务四 . ............................................................. 错误！未定义书签。
一、 实验目的
本实验主要涉及微积分中对函数特性的研究。通过实验复习函数求导法, Taylor 公式和其他有关知识。着重介绍运用建立近视似模型并进行数值计算来研究讨论函数的方法。
二、 实际问题
曲柄滑块机构是一种常用的机械结构, 它将曲柄的转动转化为滑块在直线上的往复运动, 是气压机、冲床、活塞式水泵等机械的主机构。右图为其示意图。
记曲柄OQ 的长为r , 连杆QP 的长为l , 当曲柄绕固定点O 以角速度w 旋转时, 由连杆带动滑块P 在水平槽内做往复直线运动。假设初始时刻曲柄的端点Q 位于水平线段OP 上, 曲柄从初始位置起转动的角度为θ, 而连杆QP 与OP 的锐夹角为β(称为摆角) 。在机械设计中要研究滑块的运动规律和摆角的变化规律, 确切的说, 要研究滑块的位移, 速度和加速度关于θ的函数关系, 摆角β及其角速度和角加速度关于θ的函数关系, 进而
（1）求出滑块的行程s (即滑块往复运动时左、右极限位置间的距离);
（2）求出滑块的最大和最小加速度(绝对值), 以了解滑块在水平方向上的作用力;
（3）求出β的最大和最小角加速度(绝对值), 以了解连杆的转动惯量对滑块的影响;
在求解上述问题时, 我们假定:
r =100(mm ), l =3r =300(mm ), ω=240(转/min)
符号说明:r －曲柄OQ 的长; l －连杆PQ 的长度; β－摆角(连杆PQ 与OP 的锐夹角); ω－角速度; P －滑块; x －滑块的位移; a －滑块的加速度。
三、 数学模型
取O 点为坐标原点，OP 方向为x 轴正方向，P 在x 轴上的坐标为x ，那么可用x 表示滑块的位移。利用三角关系，立即得到
x =r cos θ+l 2-r 2sin 2θ （1.1）
由于θ=ωt ，故有
dx dx d θdx ==ω （1.2） dt d θdt d θ
而
dx r 2sin θcos θ （1.3） =-r sin θ-
222d θl -r sin θ
于是滑块的速度
⎛r cos θ v =-ωr sin 1+
l 2-r 2sin 2θ⎝
进而，可以得到滑块的加速度为
a =
⎫
⎪ （1.4） ⎪⎭
dv dv =ωdt d θ
⎡⎤ （1.5） 224
r (l cos(2θ) +r sin θ) ⎥
=-ω2r ⎢cos θ+3⎢⎥
2222⎢⎥(l -r sin θ) ⎣⎦
同样，基于关系式
l sin β=r sin θ （1.6）
我们有摆角的表达式
β=arcsin sin θ⎪ （1.7）
式（1.6）对t 求导， l cos β可得
由此再得
⎛r
⎝l ⎫⎭
d βd θ=r cos θ=r ωcos θ dt dt
d βr ωcos θ
= （1.8） dt l cos β
d β⎫⎛ωsin θcos β-cos θsin β ⎪d 2βr ω dt ⎪=- （1.9） 22
l ⎪dt cos β ⎪⎝⎭
利用（1.6），不难由上两式导出
d βr ωcos θ= （1.10）
222dt l -r sin θ和
d 2βr ω2sin θ(l 2-r 2) =-3 （1.11） dt 2
(l 2-r 2sin 2θ) 2
至此，我们得到了滑块位移x 和连杆摆角β运动规律中有关变量依赖θ的表达式。
四、 数值积分方法
将位移的表达式（1.1）改写为
r 2
x =r cos θ+l (1-2sin 2θ) 2
l
r 2
一般而言，2是远比1小的数，于是利用
l
(1+ε) α=1+αε+ , 得到滑块位移的近似模型为
1
ε
r 2
x 1=r cos θ+l -sin 2θ （1.13）
2l
从而有相应的近似速度
⎛⎫dx 1dx 1d θr 2
⎪v 1===ω -r sin θ-sin(2θ) ⎪dt d θdt 2l ⎝⎭
（1.14）
r ⎛⎫
=-ωr sin θ-sin(2θ) ⎪
2l ⎝⎭
和近似加速度
dv 1r ⎛⎫
=-ω2r cos θ+cos(2θ) ⎪ （1.15） dt l ⎝⎭
这里的速度v 1和加速度a 1是直接对近似位移模型x 1求导得来的，而不是对v 和a 的精确表达式（1.4）和（1.5）的近似。当然我们也可以直接从滑块速度的解析式（1.4）进行近似。仍利用公式（1.12）有
a 1=
⎫1⎛r 2
= 1-2sin θ⎪⎪222
l -r sin θl ⎝l ⎭
1
2
-
1
2
⎫1⎛r 22
≈ 1+2sin θ⎪⎪ l ⎝2l ⎭
把上式代入（1.4），就得到滑块速度的近似模型
⎡r cos θ⎛⎫⎤r 22
v 2=-ωr sin θ⎢1+1+2sin θ⎪ ⎪⎥l 2l ⎝⎭⎦⎣
（1.16）
32
⎛r sin(2θ) r sin θsin(2θ) ⎫
⎪=-ωr sin θ++3 ⎪2l 4l ⎝⎭
从（1.16）出发，又可得近似加速度
⎡r cos(2θ) r 3(sin2(2θ) +2sin 2θcos(2θ) ⎤
a 2=-ωr ⎢cos θ++⎥ （1.17） 3
l 4l ⎣⎦
对摆角β可以利用幂级数展开的Maclaurin 公式
2
arcsin ε=ε+
ε3
6
+ , ε
得到摆角的近似模型。粗略一些，可以取
r
β1=sin θ （1.19）
l
r
（当较小时可用此式）。而必要时可以取
l
r r 3
β2=sin θ+3sin 3θ （1.20）
l 6l
相应的近似角速度为
d β1r
=ωcos θ （1.21） dt l
或
⎛r ⎫d β2r 32 =ω cos θ+3sin θcos θ⎪ （1.22） ⎪dt l 2l ⎝⎭
近似角加速度为
d 2β12r =-ωsin θ （1.23） 2
l dt
或
⎤d 2β2r 32⎡r 3
=-ω⎢sin θ+3(sinθ-sin (2θ) cos θ) ⎥ （1.24） dt 22l ⎣l ⎦
五、 实验任务
任务一
试用摆角的角加速度的三种表达式, 即式（1.11）、（1.23）和（1.24），取步长为
π
，r ，l ，ω的值如前，计算当θ∈[0,π]变化时角加速度的值, 并列表12
加以比较。
实验程序：
function m1_2(t)
r=100;l=300;w=240/60*2*pi;
a0=-r*w^2*sin(t)*(l^2-r^2)./(l^2-r^2*sin(t).^2).^(3/2) a1=-w^2*r*sin(t)/l
a2=-w^2*(r*sin(t)/l+r^3*(sin(t).^3-sin(2*t).*cos(t))/(2*l^3))
>> m1_1([0:pi/12:pi])
从结果中可以看出误差的大小, 取决于近似表达式的精度, 在利用泰勒公式求
d 2β2
近似模型时, 如果展开的精度越高, 则误差就越小, 在数据表中也可以看出,
dt
d 2β1
取得精度比高, 所以结果与真实值相差的更小。
dt
任务二
利用（1.12）式，对角摆角的角速度（1.10）式和角速度（1.11）式进行简化，将结果与（1.21）~(1.24)式进行比较，并与上题的计算结果相比较。 解析： 由式 (1+ε) α=1+αε+ , ε
d βr ωcos θ
(1.10) =
222dt l -r sin θ
⎫d β3r ω⎛r 2
⎪可以化简为 (1. 10) ' =cos θ+sin(2θ) sin θ2 ⎪dt l ⎝4l ⎭
d 2β3r ω2sin θ(l 2-r 2)
(1.11) =-3
dt 2
(l 2-r 2sin 2θ) 2
可以化简为
d 2βr ω2(l 2-r 2) ⎛3r 23⎫ (1. 11) ' =-sin θ+sin θ⎪232 ⎪dt l 2l ⎝⎭
将化简结果与（1.21）~(1.24)式进行比较，可以发现有类似的项。
实验程序：
function m2_1(t)
r=100;l=300;w=240/60*2*pi;
b0=r*w*cos(t)./sqrt(l^2-r^2*sin(t).^2) b1=w*r*cos(t)/l
b2=w*(r*cos(t)/l+r^3*sin(t).^2.*cos(t)/(2*l^3)) b3=r*w/l*(cos(t)+r^2*sin(2*t).*sin(t)/(4*l^2))
a0=-r*w^2*sin(t)*(l^2-r^2)./(l^2-r^2*sin(t).^2).^(3/2) a1=-w^2*r*sin(t)/l
a2=-w^2*(r*sin(t)/l+r^3*(sin(t).^3-sin(2*t).*cos(t))/(2*l^3))
a3=-r*w^2*(l^2-r^2)/l^3*(sin(t)+3*r^2*sin(t).^3/(2*l^2)) >> m2_1([0:pi/12:pi])
运行结果如下：
d 2β3d β3d βd β2d β
从上表中可看出与最接近真值，最接近真值3。由此看来，
dt dt dt dt dt
方案三最优。
任务三
给定一机构如右图所示。设连杆QP 长度l=300mm，曲柄OQ 的长为r=100mm，距离e=20mm，曲柄的角速度w=240转/min。对θ在一个周期（即[0, 2π]）中计算滑块的位移、行程、速度、加速度和摆角及其最值。
解析：这个机构的特点是:滑块的运动轨迹仍然在原来的平面上, 且与轴线Ox 平行, 但运动轨迹与Ox 有距离e (称为偏心距) 。这样进程时间将与退程时间不同。
由于P 点始终在直线y =e 上, 所以我们只需要考虑滑块在x 方向上的位移, 不需要再考虑在y 轴上的位移。取O 点为坐标原点，沿x 轴向右方向为正，P 在x 轴上的坐标为x, 用x 表示滑块的位移，利用三角关系有：
x =r cos θ+2-(r sin θ-e ) 2 （2.1）
由于θ=ωt ，故有
dx dx d θdx ==ω （2.2） dt d θdt d θ
dx
求导程序： d θ
>> syms r l e t
>> x=r*cos(t)+sqrt(l^2-(r*sin(t)-e).^2); >> diff(x,t)
dx r cos θ(e -r sin θ)
得 =-r sin θ+
22d θl -(e -r sin θ)
即
dx r cos θ(r sin θ-e )
（2.3） =-r sin θ-
22d θl -(r sin θ-e )
于是滑块的速度 v =-ωr sin θ-
从而，得到滑块的加速度为
a =
ωr cos θ(r sin θ-e )
l -(r sin θ-e )
2
2
(2.4)
dv dv =ωdt d θ
(2.5)
222222⎫⎛(r cos 2θ+re sin θ)(l -(r sin θ-e ) ) +r cos θ(r sin θ-e ) ⎪
=ω2 -r cos θ-
⎪(l 2-(r sin θ-e ) 2) 3⎝⎭
由关系式
l sin β+e =r sin θ (2.6）
得摆角的表达式为
β=arcsin(
r sin θ-e
) (2.7) l
滑块的行程：s =x max -x min =(l +r ) 2-e 2-(l -r ) 2-e 2=200.[1**********]
实验程序： （1）滑块的位移： function m3_1(t) r=100; l=300; e=20;
x=r*cos(t)+sqrt(l^2-(r*sin(t)-e).^2) >>m3_1([0:pi/12:2*pi])
（2）滑块的速度： function m3_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) >>m3_2([0:pi/12:2*pi])
（3）滑块的加速度： function m3_3(t) r=100; l=300;
w=240/60*2*pi; e=20;
a=w^2*(-r*cos(t)-((r^2*cos(2*t)+r*e*sin(t)).*(l^2-(r*sin(t)-e).^2)+r^2*cos(t).^2.*(r*sin(t)-e).^2)./(l^2-(r*sin(t)-e).^2).^(3/2)) >> m3_3([0:pi/12:2*pi])
（4）摆角及其最值： function m3_4(t) r=100; l=300;
w=240/60*2*pi; e=20;
b=asin((r*sin(t)-e)/l) >> m3_4([0:pi/12:2*pi])
运行结果如下：
3
由上表可看出，摆角β，在π处取得最小值0.0667；在π处取得最大值
2
0.4115。
任务四
设T1和T2分别为滑块进程和退程所需时间，试根据任务二中的数据求出T1和T2，两者是否相等？在l 和r 不变的前提下，令k=T1/T2。如果要求k=1.2，试求偏心距e 的值。
答：T 1>T 2。
由 0. 2π+2. 2arcsin
e e -2. 2arcsin =0， 200400
实验程序及结果：
>> x=fzero(@(x) 0.2*pi+2.2*asin(x./200),2) x =
-56.3465