24小时热门版块排行榜

>论坛更新日志 (2439)
>虫友互识 (212)
>导师招生 (102)
>考博 (88)
>文献求助 (76)
>休闲灌水 (68)
>仿真模拟 (45)
>基金申请 (26)
>文学芳草园 (18)
>论文投稿 (17)
>博后之家 (15)
>催化 (15)
>考研 (14)
>找工作 (11)
>公派出国 (11)
>硕博家园 (10)

返回列表

sas4444

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 3795.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]);
如何避免的出来的是复数或者是无穷大

回复此楼