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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 材料区 » 材料综合 » 材料其他 » 流体力学H-B模型 距离和流速之间的关系
51/1返回列表
查看: 3603  |  回复: 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的回帖

yy1257035418

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

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

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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +26 医学老男孩 2026-08-13 56/2800 2026-08-15 13:47 by 医学老男孩
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +7 majunge000 2026-08-11 9/450 2026-08-15 13:34 by kiz2011
[基金申请] 时间戳又变了8-15 +5 archvillain 2026-08-15 9/450 2026-08-15 12:06 by archvillain
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +4 7lpolszZVXgi 2026-08-14 8/400 2026-08-15 11:53 by KxMI1BYBxWX1
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +5 HFw0lei2R37i 2026-08-14 9/450 2026-08-15 11:33 by KxMI1BYBxWX1
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +6 Ldrop2023 2026-08-13 6/300 2026-08-15 10:09 by tuanggou
[公派出国] 售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
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
关闭
请填处理意见
关闭