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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 材料区 » 材料综合 » 材料其他 » 流体力学H-B模型 距离和流速之间的关系
231/1返回列表
查看: 3544  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 职称评审，求友友推荐见刊最快的期刊 +3 工厂打螺丝 2026-08-08 3/150 2026-08-11 13:42 by 会议编辑
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] 帮忙看看fileCode +6 wwncly 2026-08-10 12/600 2026-08-11 13:01 by wwncly
[基金申请] filecode +9 documentary 2026-08-10 9/450 2026-08-11 10:08 by wk7465
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +15 Lanmanbaby 2026-08-09 23/1150 2026-08-11 09:51 by demonstreets
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题，想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode，很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 好消息？这个有何含义？？？ +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
关闭
请填处理意见
关闭