版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

sas4444

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 3791.2
红花: 1
帖子: 281
在线: 133.5小时
虫号: 1171982
注册: 2010-12-17
性别: GG
专业: 微/纳机械系统

[交流] rayleigh-plesset仿真已有1人参与

function Rx=RP(t,R)
P0=100000;    %初始压强P0=0.1MPa
T0=298.16;    %初始温度T0=273.16+25
Hv=2257.2;    %焓Hv=2257.2KJ/Kg
R0=0.01;       %初始半径R0=0.02m
rho=1000;    %液体密度ρ=1000Kg/cm2
nu=1*10^-6;    %运动粘度ν=1.0*10^-6m2/s
gama=0.0725; %表面张力γ=0.0725N/m
m=0.018;       %摩尔质量m=0.018Kg/mol
Rg=8.3144;    %通用气体常数Rg=8.3144KJ/(Kg·mol·K)
alpha=R0/8.86;  %常数α=R0/8.86
r=0.008;       %距离r=8*10^-3m
k=1.66;       %绝热系数κ=1.66
Rx=zeros(2,1);
Rx(1)=R(2);
Rx(2)=(P0*exp((Hv*m/Rg)*(1/T0-(R(1)^3-alpha^3)^(k-1)/(T0*R0^(3*(k-1)))) ...
      +2*gama/R0*(R0^(3*(k-1))/(R(1)^3-alpha^3)^(k-1))*(R0/R(1))^(3*k)-(P0-rho/2*(R(1)^2/r^2*R(2))^2))/rho ...
      -(1.5*R(2)^2+4*nu*R(2)/R(1)+2*gama/(rho*R(1))))/R(1);

end
然后
[t,R]=ode45('RP',[0,100],[0.000007,2000]);
如何避免的出来的是复数或者是无穷大

回复此楼

» 猜你喜欢

售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有7人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有3人回复
明天应该可查了！？已经有5人回复
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》已经有7人回复
没有任何消息-是不是就凉了已经有8人回复
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急已经有3人回复
人气不行了已经有10人回复
我面上完蛋了已经有10人回复
能否退出参与的面上项目解除限项已经有24人回复
filecode，4个jtjc了已经有17人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2012-08-28 20:48:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

迷失黑夜里、

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 5.3
帖子: 41
在线: 13.3小时
虫号: 5851414
注册: 2017-03-06
性别: MM
专业: 机械工程

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

您好，您这个计算出来了吗，我近期也在用matlab解RP方程求处R和R的一阶导数

赞一下