求教阻尼比如何计算
我用MATLAB对悬臂梁建模，该如何求各阶阻尼比呢。 有关动力学＼模态理论方面书籍中有这方面的公式，不妨查查！好运！ 在动力学方程中,阻尼有特定的意义！但如要仿真测试阻尼比，又有不同方法所对应的阻尼比公式，如伴功率法、指数衰减法、频率增量法等！ 我也想知道这个问题。 能用仿真求出阻尼比？ 第一次听说！ 可以计算的，我原来下载了一个matlab程序，粘贴如下，由于忘记从那下载的了，不能注明出处了
clear;clc;
%--------------------
%约束阻尼结构单层
f=;
%-------------------------------------------------------------------------
%钢板基层
H1=6*10^(-3);
ro1=7700;
E11=19.5*10^(10);
eta1=[.1 .095 .09 .08 .065 .05 .03 .01 .004 .0024 .0019 .0015 .0012 .001 .0009 .00085 .00083 .0008 .00076 .00073 .0007 .00068 .00067 .00066 .00065];
%--------------------------------------------------------------------------
%阻尼层--------------------------------------------------------------------------
H2=4.5*10^(-3);
% %氯丁橡胶
% ro2=1300;E22=1*4*10^8;eta2=1*0.2;mu2=0.45;
% %QZD型阻尼橡胶
% ro2=1400;mu2=0.45;
% E22=[.205 .21 .215 .22 .225 .235 .256 .27 .295 .315 .337 .36 .40 .44 .48 .52 .58 .65 .72 .8 .9 1.0 1.1 1.2 1.3]*10^8;
% eta2=[.33 .35 .376 .4 .42 .45 .478 .5 .53 .565 .6 .63 .66 .68 .694 .7 .705 .706 .7 .697 .688 .67 .65 .62 .59];
%SA-3型阻尼橡胶
ro2=1370;
mu2=0.45;
eta2=1*[.8 .95 1.09 1.2 1.18 1.14 1.11 1.05 1.0 0.9 .83 .78 .7 0.64 0.58 .52.470.4.37 .31 .29 .25 .22 .21 .2];
E22=1*[.9 .95 1.1 1.3 1.6 2.0 2.5 3.0 4.0 5.0 6.4 8.0 9.4 11.5 14.0 16.5 19.5 22.0 26.0 30. 31. 32. 33. 33. 33.]*10^8;
%--------------------------------------------------------------------------
% %约束层铝板
% H3=1.0*10^(-3);
% ro3=2700;E33=7.1*10^(10);
% eta3=1*10^(-3);
%约束层钢板
H3=2*10^(-3);
ro3=7700;
E33=19.5*10^(10);
eta3=[.1 .095 .09 .08 .065 .05 .03 .01 .004 .0024 .0019 .0015 .0012 .001 .0009 .00085 .00083 .0008 .00076 .00073 .0007 .00068 .00067 .00066 .00065];
%--------------------------------------------------------------------------
H21=(H1+H2)/2;
H31=(H1+H3)/2+H2;
m1=ro1*H1;
m=m1+ro2*H2+ro3*H3;
m2m1=m/m1;
for n=1:25
E1(n)=E11*(1+j*eta1(n));
E2(n)=E22(n)*(1+j*eta2(n));
E3(n)=E33*(1+j*eta3(n));
K3(n)=E3(n)*H3;
K2(n)=E2(n)*H2;
K1(n)=E1(n)*H1;
BB(n)=real((K1(n)*H1^2+K3(n)*H3^2)/12);
w=2*pi*f(n);
G2=E2(n)/(2*(1+mu2));
p=(m/BB(n))^0.5*w;
g=G2/H2/K3(n)/p;
X(n)=(K1(n)*H31+K2(n)*(H31-H21))/(K1(n)+K2(n)/2+g*(K1(n)+K2(n)+K3(n)));%阻尼层的剪切参数
D(n)=real((K1(n)*(H21-H31/2)+(K2(n)*H21+K3(n)*H31)*g)/(K1(n)+K2(n)/2+g*(K1(n)+K2(n)+K3(n))));%第一层离中心面的距离% 常数
B1(n)=K1(n)*H1^2/12+K2(n)*H2^2/12+K3(n)*H3^2/12;
B2(n)=K1(n)*D(n)^2+K2(n)*(H21-D(n))^2+K3(n)*(H31-D(n))^2;
B3(n)=(K2(n)*H2/12+K2(n)*(H21-D(n))/2+K3(n)*(H31-D(n)))*X(n);
B(n)=B1(n)+B2(n)-B3(n);
BB1(n)=real(B(n));
BB2(n)=imag(B(n));
eta(n)=BB2(n)/BB1(n);
end
%---------------------------------------------------------------
loglog(f,eta,'k-'),grid on;hold on;xlabel('频率/Hz');ylabel('抑振量/dB');
hold on; loglog(f,eta1,'r-')
回复 6楼 msnoy 的帖子
楼上的程序好像是给定阻尼比？并非仿真求出阻尼比？
页:
[1]