哈工大机械振动大作业.docVIP

Harbin Institute of Technology 机械振动大作业 课程名称： 机械振动 班 级： 1208105 班 作 者： 付帅磊 学 号： 1120810523 指导教师： 陈照波 于东 设计时间： 2015.11.12 哈尔滨工业大学 第一次大作业 1.The amplitude of vibration of an undamped system is measured to be 1 mm. the phase shift is measured to be 2 rad and the frequency 5 rad/sec. Calculate the initial conditions. 解：由题意知，无阻尼振动位移方程为： 对无阻尼单自由度振动，，。 其中，。 解得，。 所以 图1. 2.Using the equation: “”to evaluate the constant and in terms of the initial conditions. 解： 所以，。 An automobile is modeled as 1000 kg mass supported by a stiffness k=400000N/m. When it oscillates, the maximum deflection is 10 cm. when loaded with the passengers, the mass becomes 1300 kg. Calculate the change in the frequency, velocity, amplitude, and acceleration if the maximum deflection remain 10 cm. 解：设时，汽车处于平衡位置且有向上运动趋势，即，。 则， 其中cm，， 则，，。 当车上载有乘客时，，cm。 设，。 ，， 。 图2. Use the given data to plot the response of the SDOF system. ，，。 解：对于有阻尼的单自由度振动系统 用matlab进行计算并绘图如图3. 图3. Solve the equation and plot the response. 特征方程，解得 方程的通解为，由得 所以。 图4. 第二次大作业 1.二自由度受迫振动微分方程为： 其中，。 设方程的一组稳态解为，代入以上方程得 解得： 则稳态解为： 下面求全解。 对于二自由度系统的受迫振动，其解包括齐次方程通解（自由振动的解）和特解（稳态解）。 对于方程 设，代入上式得： 考虑有非零解，则，解得： 振幅比，，取1，2 则方程的全解可写成： 其中 式中4个未知数，，，，同时有四个初始条件，即为 所以可解得： 所以振动的图像如图5所示 图5. 2. 先求稳态解 ，， 令，解得， 由得。 则 所以 ，，， m1=9;m2=1; K = [27 -3;-3 3]; C = [2.7 -0.3;-0.3 0.3]; dx = @(t,x)[x(3:4); (-K*x(1:2)-C*x(3:4)+[0;sin(3*t)])./[m1;m2]]; % 初始条件（x3~x4 为 x1~x2） x0 = [1 0 0 0]; % 求解，直接绘图 [t,x] = ode45(dx,[0 150],x0); plot(t,x(:,1:2)); lstr = arrayfun(@(i){sprintf(x%i,i)},1:2); legend(lstr{:}); legend(lstr{:},Location,S,Orientation,horizontal); grid 图7. 图8 图9. 对其局部放大，可见三种方法较为接近。 图10 对进行分析 图11. [1]任兴民,工程振动基础,北京,机械工业出版社,2006.7 [2]赵玫,机械振动与噪声学,北京,科学出版社,2004.9 [3]张义民,机械振动学基础,北京,高等教育出版社,2010.5