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

24小时热门版块排行榜    

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

超级小蜗牛

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

已搞定  谢了

发自小木虫IOS客户端
回复此楼
21楼2016-01-15 23:11:36
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭