|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
不知道为什么我的程序无法运行,不提示错误也不出结果,请高手帮忙看看,万分感谢
%受点力激振后的响应
clear
%m %轴向=s 行(1,2,3....)
%n %周向 列(0,1,2,3...)
%参数输入,
sigma=0.24; %泊松比
L=1.009; %长度
h=0.0059; %厚度
R=0.228; %半径
E=19.6e+10; %杨氏模量
q=7850; %密度
%============================================
beta=(h^2)/(12*R^2); %刚度因子
cl=(E/(q*(1-sigma^2)))^0.5; %平面振动波相速度
wr=cl/R; %环频率
%计算固有频率
r=1;
m0=1:12;
n0=0:12;
for m=1:length(m0)
for n=1:length(n0)
km=m0(m)*pi/L;
A0=((1-sigma)/2)*((1-sigma^2)*(km*R)^4+beta*((km*R)^2+n0(n)^2)^4);
A1=((1-sigma)/2)*(((km*R)^2+n0(n)^2)^2+n0(n)^2+(3+2*sigma)*(km*R)^2)+beta*((3-sigma)/2)*((km*R)^2+n0(n)^2)^3;
A2=1+((3-sigma)/2)*((km*R)^2+n0(n)^2)+beta*((km*R)^2+n0(n)^2)^2;
%固有频率方程x^3-A2*x^2+A1*x-A0=0
P=[1,-A2,A1,-A0];
X=roots(P);
Omiga(r,3)=X(3)^0.5; %无量纲频率
wn=wr*Omiga(r,3);
lwn(r)=wn/(2*pi);%固有频率
a(r)=m0(m);
b(r)=n0(n);
r=r+1;
end
end
%排序开始
M1=length(m0)*length(n0);
for i1 = 1:M1
for j1=i1+1:M1
if lwn(i1)>=lwn(j1)
d=lwn(i1);
lwn(i1)=lwn(j1);
lwn(j1)=d;
e=a(i1);
a(i1)=a(j1);
a(j1)=e;
f=b(i1);
b(i1)=b(j1);
b(j1)=f;
end
end
end
for k1=1:M1
m1(k1)=a(k1);n1(k1)=b(k1);%排序后的模态序数
fm(k1)=lwn(k1);wm(k1)=2*pi*lwn(k1);%排序后频率
end |
评分
-
1
查看全部评分
-
|