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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 材料区 » 材料综合 » 材料其他 » 流体力学H-B模型 距离和流速之间的关系
51/1返回列表
查看: 3595  |  回复: 22
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

超级小蜗牛

新虫 (小有名气)
已领完
流体力学H-B模型 距离和流速之间的关系
领取红包 (小木虫手机app专属红包)

扫一扫，下载小木虫客户端

流体力学H-B模型 距离和流速之间的关系 Matlab总是算不出结果% H-B模型平板流动速度
clc;
clear;
global n P l R
t0=3e3;
R=0.05;
n=0.6;
m=1.26;
l=0.1;
k=3.89;
P=5.5e4;
H=(1-n)*R;
w=(1+n)*R*pi;
r=0.03:0.0004:0.05;
c=length(r);
u1=zeros(1,c);
center=(1+n)*R/2;
for i=1:c
    h=t0*l/P              %   此处是hc=b-a的大小的 一半
    tw=H*P/(2*l)
    if r(i)<=center-h
        u1(i)=-l*(-t0-P*(r(i)-center)/l).^(m+1)/((m+1)*P*k.^m)+(-t0+tw).^(m+1)*l/((m+1)*P*k.^m);  %u1
    else
      if  r(i)>=center+h
        u1(i)=l*(-t0+tw).^(m+1)/((m+1)*P*k.^m)-l*((r(i)-center)*P/l-t0).^(m+1)/((m+1)*P*k.^m);  %u2
      else u1(i)=l*(tw-t0).^(m+1)/((m+1)*P*k.^m);          %uc
      end
    end
end
figure;
plot(r,u1,'--r');
grid on;
axis([0.03 0.05 0 15]);
xlabel('Distance from center linr(m)');
ylabel('Velocity(m/s)')

hold on
n=0.6;
R=0.05;
m=1.26;
x0=[0.037,0.047];
x=fsolve(@(x)myfun2(x,t0),x0);
r=0.03:0.0004;0.05;
c=length(r);
u2=zeros(1,c);

A=P*x(2)*x(1)/(2*l);

for  i=1:c
            
if r(i)<=x(1)
     u2(i)=quad(@(r)(abs(1/k*(-P*r/(2*l)+A./r-t0))).^m,n*R,r(i));            %u1
    else
        if r(i)>=x(2)
        u2(i)=quad(@(r)(abs((-1/k.^m)*(P*r/(2*l)-A./r-t0))).^m,R,r(i));        %u2
        else
            u2(i)=quad(@(r)(abs(1/k*(-P*r/(2*l)+A./r-t0))).^m,(n*R),x(1));   %uc
        end
    end
end
figure;
plot(r,u2,'-b');
grid on;
axis([0.03 0.05 0 15]);
xlabel('Distance from center linr(m)');
ylabel('Velocity(m/s)')
                  
function y=myfun2(x,t0)
global R n m l P
y=[x(2)-x(1)-2*l*t0/P;
    integral(@(r)(abs(-P*r./(2*l)+P*x(1)*x(2)./(2*l*r)-t0)).^m,(n*R),x(1))+integral(@(r)(abs(P*r./(2*l)-P*x(1)*x(2)./(2*r.*l)-t0)).^m,R,x(2))];
end
求大神帮忙

发自小木虫IOS客户端
回复此楼

» 猜你喜欢

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2016-01-15 22:10:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangliang521

木虫 (著名写手)
看不懂……

发自小木虫IOS客户端
回复此楼
16楼2016-01-15 22:43:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 23 个回答

yy1257035418

金虫 (正式写手)
看不懂啊

发自小木虫Android客户端
回复此楼
12楼2016-01-15 22:13:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woshilaoaa

至尊木虫 (文坛精英)
教头

发自小木虫Android客户端
回复此楼
13楼2016-01-15 22:20:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

woodtower

铁杆木虫 (著名写手)
帮顶

发自小木虫Android客户端
回复此楼
14楼2016-01-15 22:24:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
v颂华v2楼
2016-01-15 22:11   回复  
超级小蜗牛(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
牛少好轻功3楼
2016-01-15 22:11   回复  
超级小蜗牛(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
qiaoyiyan4楼
2016-01-15 22:11   回复  
超级小蜗牛(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
阿狸买卖猫5楼
2016-01-15 22:11   回复  
超级小蜗牛(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
wahaha91156楼
2016-01-15 22:11   回复  
超级小蜗牛(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
sioc0077楼
2016-01-15 22:11   回复  
超级小蜗牛(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
zhke5238楼
2016-01-15 22:11   回复  
超级小蜗牛(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
兽医做科研9楼
2016-01-15 22:11   回复  
超级小蜗牛(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
小小房屋910楼
2016-01-15 22:11   回复  
超级小蜗牛(金币+1): 谢谢参与
已获得1个金币 发自小木虫IOS客户端
xizheng11楼
2016-01-15 22:11   回复  
超级小蜗牛(金币+1): 谢谢参与
已获得1个金币 发自小木虫Android客户端
查看全部 23 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 07:09 by 4wMiSEwB6436
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 k0dTPqJtl0jt 2026-08-14 5/250 2026-08-15 04:45 by 4wMiSEwB6436
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 k0dTPqJtl0jt 2026-08-14 4/200 2026-08-15 01:28 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 HFw0lei2R37i 2026-08-14 7/350 2026-08-15 00:40 by 4wMiSEwB6436
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 HFw0lei2R37i 2026-08-14 8/400 2026-08-15 00:16 by 4wMiSEwB6436
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 HFw0lei2R37i 2026-08-14 5/250 2026-08-14 22:32 by 4wMiSEwB6436
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分？ +11 yufeiwaner 2026-08-09 13/650 2026-08-12 08:17 by yufeiwaner
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
关闭
请填处理意见
关闭