机械系统动力学三级项目报告(4)

2.3.2 MATLAB编程

编程如下：

m=1.6632;k=15; wn=sqrt(k/m); n=wn;

x0=10;v0=0; c1=x0;

c2=v0+n*c1; t=0:0.02:5;

x=exp(-n*t).*(c1+c2*t);

Plot(t,x);

所得曲线如下：

第 15 页 共 29 页

2.3.3 simulink仿真：

第 16 页 共 29 页

2.4大阻尼

2.4.1 SolidWorks运动仿真

2.4.2 MATLAB编程 编程如下：

m=1.6632;k=15; wn=sqrt(k/m); n=5;

[c1,c2]=solve('-c1-c2=10','1.0023*c1+8.9977*c2=0','c1','c2'); t=0:0.01:5;

x1=(-c1*exp(-n*t+t*sqrt(n^2-wn^2))); x2=(-c2*exp(-n*t+t*sqrt(n^2-wn^2))); x=x1+x2; Plot(t,x)；

所得曲线如下：

第 17 页 共 29 页

2.4.3 simulink仿真：

第 18 页 共 29 页

3、受迫振动 3.1无阻尼

3.1.1Soloidworks运动仿真

3.1.2 MATLAB编程

运动方程如下：

由初始条件可求得A和θ。

程序如下：

m=1.6632; k=15; H=1;

wn=sqrt(k/m);

第 19 页 共 29 页