应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求大神帮我看看以下matlab程序
11/1返回列表
查看: 337  |  回复: 0
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

xiaowei619

新虫 (初入文坛)

[求助] 求大神帮我看看以下matlab程序

我这有一个求timoshenko梁固有频率的程序，如下：
clc
clear all;
syms eigen;
E = 2.1*10^5;
u = 0.3;
d = 7.8*10^-6;
G = E/2*(1+u);
rout = 66;
rin = 54;
m = rin/rout;
A = pi*(rout^2-rin^2);
L = 26;
I = pi*(rout^4-rin^4)/4;
k = 6*(1+u)*(1+m^2)^2/((7+6*u)*(1+m^2)^2+(20+12*u)*m^2);
ssqu = E*I/(k*A*G*L^2);
rsqu = I/(A*L^2);
bsqu = d*A*L^4/(E*I)*eigen;
Q = bsqu*(ssqu+rsqu)/2;
yib = sqrt(0.25*bsqu^2*(ssqu+rsqu)^2-bsqu*(bsqu*ssqu*rsqu-1));
afa = sqrt(Q+yib);
bda = sqrt(-Q+yib);
lam = (afa-bsqu*ssqu/afa);
det = (bda+bsqu*ssqu/bda);
D11 = (afa-lam)*(cos(afa)-cosh(bda));
D12 = (lam-afa)*sin(afa)+lam*afa*(bda-det)*sinh(bda)/(bda*det);
D21 = -lam*afa*sin(afa)+(afa-lam)*det*bda*sinh(bda)/(det-bda) ；
D22 = lam*afa*(cosh(bda)-cos(afa));
solve(D11*D22-D12*D21);
未知变量eigen，方程D11*D22-D12*D21=0;
solve里面用eigen的式子带进去貌似也得不到结果？
求大神给我指点，不甚感激！

[ Last edited by xiaowei619 on 2014-3-4 at 11:24 ]
回复此楼

» 猜你喜欢
1楼 2014-03-03 16:36:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

科研从小木虫开始，人人为我，我为人人
相关版块跳转 我要订阅楼主 xiaowei619 的主题更新
11/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭