作业:如图1所示是一个单自由系统附件一个减振器形成的的两自由振动系统,已知m1=105kg,m2=7kg,k1=10000n/m,k2=410n/m,c2=1.15n·m-1·s,f1(t)=f1ejωt。
求:(简化为粘性比例阻尼进行实模态分析)
1. 物理坐标下的振动微分方程;
2. 频响函数矩阵;
3. 频响函数的模态展式矩阵;
4. 脉冲相应函数;
5. 画出h11(ω)的幅频特性曲线,相频特性曲线,实频特性曲线,虚频特性曲线,nyquist图,bode图;
6. 固有频率,阻尼固有频率;
7. 画出振型图;
8. 模态坐标系下的振动微分方程;
9. 模态参数:复模态质量,复模态刚度,复模态阻尼。
10.按实模态系统,给出灵敏度分析。
11.集全班同学的数据(必要的话再补做不同m2,k2,c2参数下的数据,画出x1的最大振幅与m2,k2,c2,的变化曲线,从而分析出减振器的最佳参数。
解:1.振动微分方程。
对质量m1、m2绘分离体图(如图1-1),用牛二定律列分离体在铅垂方向的力平衡方程得。
将(1.1)整理可得:
且m1=105、m2=7、k1=10000、k2=410、c2=1.15,代入(1.2)得:
可以得出此二自由度系统振动微分方程为:
其中m=;c=;k=;f(t)= 图1-1、系统的分离体图。
2.频响函数矩阵。
由书p25(1.4-58)公式可知,此二自由度系统频响函数矩阵为一2×2方阵,其表达式为:,其中m=;c=;k2.1)
写成矩阵形式:
3.频响函数的模态展式矩阵。
1)求解瑞利阻尼矩阵。
由于粘性阻尼矩阵c无法进行正交性对角化,故不能直接应用坐标变换将(1.3)解耦。由于在该题中,粘性阻尼相对很小,对于小阻尼振动系统,可以利用瑞利比例阻尼来代替粘性阻尼,以获得可对角化的阻尼矩阵。
1)瑞利比例阻尼系数的确定。
瑞利比例阻尼:,其中;、为瑞利比例阻尼系数。
瑞利比例阻尼系数存在以下关系:
其中为圆频率(为系统固有频率,书中表示为);为阻尼比。
将上式写为矩阵形式:
可得:,其中3.1)
由此可知,只要我们确定了一个系统任意两阶的固有频率及其阻尼比,就可以确定出瑞利比例阻尼系数,从而得到瑞利比例阻尼矩阵。
2)求该二阶系统的。
一、二阶固有频率及其阻尼比。
利用求解该系统振动微分方程的特征值来确定固有频率及其阻尼比。由书p23(1.4-43)-(1.4-46)公式为求解步骤,下面利用matlab来计算固有频率和阻尼比:
编写matlab程序求特征方程,程序如下:
syms x;
m1=105; m2=7; k1=10000; k2=410; c2=1.15;
m=[m1 0;
0 m2];
c=[c2 -c2;
c2 c2];
k=[k1+k2 -k2;
k2 k2];
y=det(m*x^2+c*x+k)
解以上求得的多项式:
> p=[735 644 115920 11500 4100000];
> x0=roots(p)
由特征值可得:、
3)求瑞利比例阻尼系数及瑞利比例阻尼矩阵。
根据公式(3.1)编写matlab程序求解特征方程,程序如下:
function cr=rayleigh()
--计算瑞利阻尼系数alpha和beta--
xi1=0.0396; xi2=0.0707; f1=7.2650; f2=10.2882;
omega1=2*pi*f1;
omega2=2*pi*f2;
a=[1/(2*omega1) omega1/2;
1/(2*omega2) omega2/2];
xi=[xi1;
xi2];x=inv(a)*xi;
alpha=x(1,1)
beta=x(2,1)
--计算瑞利阻尼矩阵cr(2*2)--
alpha=-1.8801;
beta=0.0026;
m1=105; m2=7; k1=10000; k2=410;
m=[m1 0;
0 m2];
k=[k1+k2 -k2;
k2 k2];
cr=alpha*m+beta*k;
可知:瑞利比例阻尼系数、
瑞利比例阻尼矩阵。
2)求解模态矩阵(及特征矢量矩阵)
因此,我们可以利用求此二阶系统在无阻尼振动下的微分方程的特征矢量更简单的得出模态矩阵。
改写matlab程序求解此二阶系统在无阻尼振动下的微分方程的特征方程,程序如下:
syms x;
m1=100; m2=5; k1=10000; k2=500;
m=[m1 0;
0 m2];
k=[k1+k2 -k2;
k2 k2];
y=det(k-x^2*m)
解以上求得的多项式:
> p=[735 -115920 4100000];
> x0=roots(p)
可知:、。将其分别代入回,可得:
求得模态矩阵。
3)求解频响函数的模态展式矩阵。
1)求模态质量矩阵、模态刚度矩阵和模态阻尼矩阵。
2)由此可得频响函数的模态展式为:
写成矩阵形式为:
将所求、、、代入:
4.脉冲响应函数。
对(3.3)作傅立叶逆变换,得到脉冲响应函数矩阵:
5.的幅频、相频、实频、虚频特性曲线以及导纳图和博德图。
1)的幅频特性曲线:与的关系。
其中,为阻尼比。代入可得:
编写matlab程序画图,程序如下:
omega=5:.01:15;
o1=omega/sqrt(52.7802); o2=omega/sqrt(105.8471);
k1=56721.8062; k2=121340; xi1=0.0396; xi2=0.0707;
y11=1./(k1.*sqrt((1-o1.
^2).^2+(2.*xi1.
*o1).^2))+1./(k2.
*sqrt((1-o2.^2).^2+(2.
*xi2.*o2).^2));
plot(omega,y11,'linewidth',2);
grid on
xlabel('频率 hz')
ylabel('幅值 mm')
title('m1的一阶幅频特性')
输出图形:2)的相频特性曲线:与的关系。代入可得:
编写matlab程序画图,程序如下:
omega=5:.01:15;
o1=omega/sqrt(52.7802); o2=omega/sqrt(105.8471);
k1=56721.8063; k2=121340; xi1=0.0396; xi2=0.0707;
y11=atan((-2.*xi1.*o1)./1-o1.^2))+atan((-2.*xi2.*o2)./1-o2.^2));
plot(omega,y11,'linewidth',2);
axis([5 15 -2 2])
grid on
xlabel('频率 hz')
ylabel('相位角')
title('m1的一阶相频特性')
输出图形:3)的实频特性曲线:与的关系。代入可得:
编写matlab程序画图,程序如下:
omega=5:.01:15;
o1=omega/sqrt(52.7802); o2=omega/sqrt(105.8471);
k1=56721.8062; k2=121340; xi1=0.0396; xi2=0.0707;
y11=(1-o1.^2)./k1.
*(1-o1.^2).^2+(2.
*xi1.*o1).^2))+1-o2.
^2)./k2.*(1-o2.
^2).^2+(2.*xi2.
*o2).^2));
plot(omega,y11,'linewidth',2);
grid on
xlabel('频率 hz')
ylabel('幅值 mm')
title('m1的一阶实频特性')
输出图形:4)的虚频特性曲线:与的关系。代入可得:
编写matlab程序画图,程序如下:
omega=5:.01:15;
o1=omega/sqrt(52.7802); o2=omega/sqrt(105.8471);
k1=56721.8062; k2=121340; xi1=0.0396; xi2=0.0707;
y11=(1-o1.^2)./k1.
*(1-o1.^2).^2+(2.
*xi1.*o1).^2))+1-o2.
^2)./k2.*(1-o2.
^2).^2+(2.*xi2.
*o2).^2));
plot(omega,y11,'linewidth',2);
axis ([5 15 -0.0003 0.0003])
grid on
xlabel('频率 hz')
ylabel('幅值 mm')
title('m1的一阶虚频特性')
输出图形:5)的导纳图:与的关系。
编写matlab程序画图,程序如下:
omega=0:.01:15;
机械结构分析设计作业
柳州职业技术学院课程设计。任务书。系 部专业班学生学号。一 课程设计题目 带式运输机上的单级圆柱齿轮减速器设计。二 课程设计进行的日期 年月日起至年月日止 三 课程设计指导教师。四 任务书的内容 1 目的 1 通过课程设计培养正确的设计思想,掌握机械设计的基本方法,巩固 提高及综合运用本课程及先修课...
场上作业机械项目运营分析报告
场上作业机械项目。运营分析报告。规划设计 投资分析 第一章项目总体情况说明。一 经营环境分析。1 工业4.0是一场全球范围的技术革命,也是中国制造业转型升级的必经之路。作为中国在新一轮科技革命和产业变革背景下针对制造业发展提出的重要战略举措,中国制造2025 不只是2015年短期的主题性投资,更有其...
农业机械跨区作业距离的分析
吴利华 刘占良王。辉。内蒙古农机推广站 河北农业大学。摘要 农机跨区作业是在市场经济体制下农机资源优化配置的模式。文中从宏观和微观经。济的角度分析了影响农机跨区作业范围因素 并建立了分析模型。分析表明 农业机械跨区作业最远距离与机械的年折旧费用成正比,而与单位转移路途的耗费成反比 在市场经济条件下,...