声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1606|回复: 6

[混合编程] 圆柱壳固有频率计算

[复制链接]
发表于 2010-12-9 20:17 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

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

查看全部评分

回复
分享到:

使用道具 举报

发表于 2010-12-9 21:04 | 显示全部楼层
"程序无法运行"? 可以啊!
的确没报错, 当然不提示错误
"不出结果"? LZ没要输出, 程序怎会自动输出结果!?:@)
 楼主| 发表于 2010-12-10 09:32 | 显示全部楼层
回复 2 # Happy99 的帖子

这个要是要输出应该下面加什么啊???求的是固有频率lwn(r)
 楼主| 发表于 2010-12-10 10:03 | 显示全部楼层
回复 2 # Happy99 的帖子

这个是输出结果的前面一部分,要输入什么才能得到这个结果呢
23.jpg
发表于 2010-12-10 11:49 | 显示全部楼层
回复 3 # zhufeng9237 的帖子

很难想像写出这程序了, 但不知道如何输出!:@)
在command window中直接打LZ要的参数!
 楼主| 发表于 2010-12-10 15:49 | 显示全部楼层
回复 2 # Happy99 的帖子

这个我知道怎么回事了,谢谢了。但是知道lwn怎么找出对应的a(r),b(r)呢??
发表于 2010-12-10 22:33 | 显示全部楼层
回复 6 # zhufeng9237 的帖子

不是一样吗!? (same as 5F)
matlab应该属於直译交谈式的软件, 只要有基本的矩阵数学概念, 真的非常容易上手, 但总是要"亲自试一试"!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 08:17 , Processed in 0.064134 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表