机械振动大作业-求初始激的自由振动响应.docxVIP

图示系统中, m1＝m2＝m3＝m, k1＝k2＝k3＝k, 设初始位移为1, 初始速度为0, 求初始激励的自由振动响应。要求：利用影响系数法求解刚度阵K和质量阵M，建立控制方程；(15分)求解系统固有频率和基准化振型；(13分)求解对初始激励的响应（运动方程）；(12分)利用软件仿真对初始激励响应曲线(Matlab，simulink，excel均可)，给出仿真程序(或框图)、分析结果；尝试对m、k赋值，分析曲线变化；(10分)浅谈对本课程的理解、体会，对授课的意见、建议；(10分)字迹清晰，书写规整。(10分)（1）利用影响系数法求解刚度阵K和质量阵M，建立控制方程；①求解刚度矩阵K令，则弹簧变形量δ=[1 1 0]T，在此条件下系统保持平衡，按定义需加于三物块的力如图所示根据平衡条件可得同理，令得令得故刚度矩阵为②求解质量矩阵M令得，，令得，，令得，，故质量矩阵为③建立控制方程应用叠加原理可得：（2）求解系统固有频率和基准化振型；①求解固有频率令则将代入解得②求解基准化振型将ω1、ω2、ω3分别代入([K]-ω2[M]){X}=0得：则（3）求解对初始激励的响应（运动方程）；对初始条件标准化：标准坐标下的初始激励响应：广义坐标下的初始激励响应（4）利用软件仿真对初始激励响应曲线(Matlab，simulink，excel均可)，给出仿真程序(或框图)、分析结果；尝试对m、k赋值，分析曲线变化；利用MATLAB对初始激励响应曲线在m=1, k=1; m=1, k=10; k=10, m=1; m=1, k=10 四种情况下进行仿真，仿真源程序与仿真结果见附录。当m=1，k=1时，仿真结果如图1所示，可以看出x1，x2，x3均呈现周期性变化；当m=1，k=10时，仿真结果如图2所示,1、2两图相比较可知，变化周期与k负相关；当m=10，k=1时，仿真结果如图3所示,1、3两图相比较可知，变化周期与m正相关；当m=10，k=10时，仿真结果如图4所示,1、4两图相比较可知，m、k增大倍数相同时，变化周期不变。附录：MATLAB仿真程序k=1 m=1 时的初始激励响应曲线仿真程序如下：clear allclose allclck=1;m=1;w1=sqrt(0.198*k/m)w2=sqrt(1.555*k/m)w3=sqrt(3.247*k/m)t=linspace(0,100,10000);x1=0.543.*cos(w1.*t)+0.349.*cos(w2.*t)+0.107.*cos(w3.*t)x2=0.979.*cos(w1.*t)+0.155.*cos(w2.*t)-0.134.*cos(w3.*t)x3=0.122.*cos(w1.*t)-0.28.*cos(w2.*t)+0.059.*cos(w3.*t)plot(t,x1,m-,t,x2,r-,t,x3,k-);title(k=1 m=1时的初始激励响应曲线);xlabel(t/s);ylabel(x(t)/m);legend(x1,x2,x3);②k=1 m=10时的初始激励响应曲线clear allclose allclck=1;m=10;w1=sqrt(0.198*k/m)w2=sqrt(1.555*k/m)w3=sqrt(3.247*k/m)t=linspace(0,100,10000);x1=0.543.*cos(w1.*t)+0.349.*cos(w2.*t)+0.107.*cos(w3.*t)x2=0.979.*cos(w1.*t)+0.155.*cos(w2.*t)-0.134.*cos(w3.*t)x3=0.122.*cos(w1.*t)-0.28.*cos(w2.*t)+0.059.*cos(w3.*t)plot(t,x1,m-,t,x2,r-,t,x3,k-);title(k=1 m=10时的初始激励响应曲线);xlabel(t/s);ylabel(x(t)/m);legend(x1,x2,x3);③k=10 m=1时的初始激励响应曲线clear allclose allclck=10;m=1;w1=sqrt(0.198*k/m)w2=sqrt(1.555*k/m)w3=sqrt(3.247*k/m)t=linspace(0,100,10000);x1=0.543.*cos(w1.*t)+0.349.*cos(w2.*t)+0.107.*cos(w3.*t)x2=0.979.*cos(w1.*t)+0.155.*cos(w2.*t)-0.134.*cos(w3.*t)x3=0.122.*cos(w1.*t)-0.28.*cos(w2.*t)+0.059.*cos(w3.*t)plot(t,x1,m-,t,x2,r-,t,x3,k-);