圆柱壳固有频率计算

zhufeng9237

不知道为什么我的程序无法运行，不提示错误也不出结果，请高手帮忙看看，万分感谢
%受点力激振后的响应
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

Happy99

"程序无法运行"? 可以啊!
的确没报错, 当然不提示错误
"不出结果"? LZ没要输出, 程序怎会自动输出结果!?:@)

zhufeng9237

回复 2 # Happy99 的帖子

这个要是要输出应该下面加什么啊？？？求的是固有频率lwn(r)

zhufeng9237

回复 2 # Happy99 的帖子

这个是输出结果的前面一部分，要输入什么才能得到这个结果呢

Happy99

回复 3 # zhufeng9237 的帖子

很难想像写出这程序了, 但不知道如何输出!:@)
在command window中直接打LZ要的参数!

zhufeng9237

回复 2 # Happy99 的帖子

这个我知道怎么回事了，谢谢了。但是知道lwn怎么找出对应的a(r),b(r)呢？？

Happy99

回复 6 # zhufeng9237 的帖子

不是一样吗!? (same as 5F)
matlab应该属於直译交谈式的软件, 只要有基本的矩阵数学概念, 真的非常容易上手, 但总是要"亲自试一试"!