24小时热门版块排行榜    

查看: 851  |  回复: 6
本帖产生 1 个 机械EPI ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

Remantic

木虫 (正式写手)

元始天尊

★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
茉莉新西兰: 金币+5, 机械EPI+1, 非常专业的应助,辛苦了,奖励EPI一枚 2013-07-18 22:43:20
实际上应是a1与kb的关系曲线,因为kb只是a1的函数,属于变下限积分问题,偶编了个MATLAB程序,你看看对不对:
% 小木虫
% http://muchong.com/bbs/viewthread.php?tid=6130700&pid=1#pid1
% date 2013-07-18
clc;clear;
kb = zeros(size(0:0.05:0.5));
a2 = 0.0826;
E =2.1e11;L=90;
i=0;
for a1 =0:0.05:0.5;
    i=i+1;
    FH = @(a)3*(1+cos(a1).*((a2-a)).*sin(a)-cos(a)).^2.*(a2-a).*cos(a)./(2*E*L*(sin(a)+(a2-a).*cos(a)).^3);
    kb(i) = quadl(FH,-a1,a2,1e-6);
end
plot(0:0.05:0.5,1./kb,'r','linewidth',3);
xlabel('\alpha ');ylabel('Kb');
Matlab编程
untitled.jpg

» 本帖已获得的红花(最新10朵)

最艰难的道路往往通往最高处。
4楼2013-07-17 19:29:39
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 285178048 的主题更新
信息提示
请填处理意见