圆柱壳固有频率计算
不知道为什么我的程序无法运行，不提示错误也不出结果，请高手帮忙看看，万分感谢%受点力激振后的响应
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=;
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 "程序无法运行"? 可以啊!
的确没报错, 当然不提示错误
"不出结果"? LZ没要输出, 程序怎会自动输出结果!?:@) 回复 2 # Happy99 的帖子
这个要是要输出应该下面加什么啊？？？求的是固有频率lwn(r) 回复 2 # Happy99 的帖子
这个是输出结果的前面一部分，要输入什么才能得到这个结果呢
回复 3 # zhufeng9237 的帖子
很难想像写出这程序了, 但不知道如何输出!:@)
在command window中直接打LZ要的参数! 回复 2 # Happy99 的帖子
这个我知道怎么回事了，谢谢了。但是知道lwn怎么找出对应的a(r),b(r)呢？？ 回复 6 # zhufeng9237 的帖子
不是一样吗!? (same as 5F)
matlab应该属於直译交谈式的软件, 只要有基本的矩阵数学概念, 真的非常容易上手, 但总是要"亲自试一试"!
页:
[1]