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

24小时热门版块排行榜    

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

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的回帖

爱草莓

至尊木虫 (著名写手)
看着都眼晕

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

wangliang521

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

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

超级小蜗牛

新虫 (小有名气)
引用回帖:
12楼: Originally posted by yy1257035418 at 2016-01-15 22:13:01
看不懂啊

已搞定  谢了

发自小木虫IOS客户端
回复此楼
17楼2016-01-15 23:10:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

超级小蜗牛

新虫 (小有名气)
引用回帖:
16楼: Originally posted by wangliang521 at 2016-01-15 22:43:20
看不懂……

已搞定  谢了

发自小木虫IOS客户端
回复此楼
18楼2016-01-15 23:11:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

超级小蜗牛

新虫 (小有名气)
引用回帖:
15楼: Originally posted by 爱草莓 at 2016-01-15 22:36:49
看着都眼晕

已搞定  谢了

发自小木虫IOS客户端
回复此楼
19楼2016-01-15 23:11:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

超级小蜗牛

新虫 (小有名气)
引用回帖:
14楼: Originally posted by woodtower at 2016-01-15 22:24:06
帮顶

已搞定  谢了

发自小木虫IOS客户端
回复此楼
20楼2016-01-15 23:11:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

超级小蜗牛

新虫 (小有名气)
引用回帖:
13楼: Originally posted by woshilaoaa at 2016-01-15 22:20:54
教头

已搞定  谢了

发自小木虫IOS客户端
回复此楼
21楼2016-01-15 23:11:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangliang521

木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
老兄，可否把这个案例分享一？

发自小木虫IOS客户端
赞一下 回复此楼
22楼2016-01-16 08:16:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

超级小蜗牛

新虫 (小有名气)
引用回帖:
22楼: Originally posted by wangliang521 at 2016-01-16 08:16:20
老兄，可否把这个案例分享一？

不是什么大问题  只是有个标点符号错了

发自小木虫IOS客户端
赞一下 回复此楼
23楼2016-01-16 13:22:32
已阅   回复此楼   关注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客户端
相关版块跳转 我要订阅楼主 超级小蜗牛 的主题更新
231/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:39 by 9i7PCSZyMzot
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:15 by 9i7PCSZyMzot
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 01:03 by 9i7PCSZyMzot
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 E4V9bALS54sj 2026-09-19 3/150 2026-09-20 00:51 by 9i7PCSZyMzot
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 3/150 2026-09-20 00:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:51 by 9i7PCSZyMzot
[基金申请] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:39 by 9i7PCSZyMzot
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:27 by 9i7PCSZyMzot
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-19 4/200 2026-09-19 23:03 by 9i7PCSZyMzot
[公派出国] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 4pRyVGCVerAA 2026-09-18 7/350 2026-09-19 21:03 by 9i7PCSZyMzot
[基金申请] 国社科系统bug了，是不是要放榜了？ +7 kynobel 2026-09-16 9/450 2026-09-19 14:52 by 接好运婷
[考博] 申博发邮件 +8 Lmengk 2026-09-15 9/450 2026-09-19 09:27 by 中间体超市
[有机交流] 各位大神，目前国内有哪些比较好用的逆合成软件？ 20+3 zapen 2026-09-15 9/450 2026-09-18 18:17 by zapen
[考博] 上海工程技术大学激光智能制造课题组｜2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +8 两三岁ss 2026-09-14 8/400 2026-09-16 09:46 by 水士口
信息提示
关闭
请填处理意见
关闭