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

1.3 MATLAB编程

运动方程如下：

由初始条件确定X11 编程如下：

m=1.6632; k=15; q10=0; q20=10;

p1=sqrt((3-sqrt(5))*k/(2*m)); p2=sqrt((3+sqrt(5))*k/(2*m)); u1=(p1^2*m-2*k)/(-k);

第 25 页 共 29 页

X12 ?1 ?2

u2=(p2^2*m-2*k)/(-k);

X12=(1/(u2-u1))*sqrt((u2*q10-q20)^2); X11=(1/(u1-u2))*sqrt((u1*q10-q20)^2); a1=-pi/2; a2=-pi/2; t=0:0.01:5;

x=-(u1*X11*sin(p1*t+a1)+u2*X12*sin(p2*t+a2)); Plot(t,x)

2、受迫振动

2.1 Solidworks运动仿真

第 26 页 共 29 页

2.2 simulink仿真

2.3 MATLAB编程

运动方程如下：

第 27 页 共 29 页

由初始条件确定X11 编程如下：

X12 ?1 ?2

m=1.6632; k=15; F0=0.5; w=2*pi; q10=0; q20=0;

X1=k*F0/((k-m*w^2)*(k-m*w^2)-k^2);

X2=(2*k-m*w^2)*F0/((2*k-m*w^2)*(k-m*w^2)-k^2); p1=sqrt((3-sqrt(5))*k/(2*m)); p2=sqrt((3+sqrt(5))*k/(2*m)); u1=(p1^2*m-2*k)/(-k); u2=(p2^2*m-2*k)/(-k);

X11=(1/(u1-u2))*sqrt((u2*X1*w-w*X2)^2/p1^2); X12=(1/(u2-u1))*sqrt((u1*X1*w-w*X2)^2/p2^2); a1=0; a2=0; t=0:0.01:5;

x=u1*X11*sin(p1*t+a1)+u2*X12*sin(p2*t+a2)+X2*sin(w*t); Plot(t,x);

出图如下：

第 28 页 共 29 页