- 积分
- 0
- 注册时间
- 2011-8-25
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2011-9-20 14:45:11
|
显示全部楼层
来自 韩国
本帖最后由 南北左右 于 2011-9-20 15:37 编辑
TBE_Legend 发表于 2011-9-17 18:58
你把word文档压缩成 rar 文件就能上传了。
C无穷大 是多少??
你好,根据你的程序,我稍微调整了下参数
g = 98/10;
\[Rho]L = 800;
mu = 25/100;
h = 20;
d = 12 10^-10;
\[Sigma] = 3/1000;
P0 = 300000;
Tb = 293;
T0 = 298;
\[Rho]g0 = 668/1000;
Cs = 296/1000;
Cinf = 125/10000;
re = 2 \[Rho]L R[t] u[t]/mu;
Cd = 16/re (1 + (8/re + 1/2 (1 + (3315/1000)/re^(1/2)))^-1);
Pb = (\[Rho]g[t] Tb)/(\[Rho]g0 T0) 101300;
Sh = 2/Sqrt[Pi] (1 - 2/(3 (1 + 9/100 re^(2/3))^(3/4)))^(1/2) (25/10 +
Sqrt[(2 R[t] u[t])/d]);
ode = {u'[t] == 2 g - 3/(4 R[t]) Cd u[t]^2,
R'[t] ==
R[t]/(4 mu) (Pb -
P0 - \[Rho]L g (h - s[t]) - (2 \[Sigma])/R[t]), \[Rho]g'[
t] == -((3 \[Rho]L Sh d (Cs - Cinf))/(2 R[t]^2)) - (3 \[Rho]g[
t])/R[t] R'[t], s'[t] == u[t], u[0] == 2/1000,
R[0] == 2/1000, \[Rho]g[0] == 32/10, s[0] == 1/10};
sol = First[NDSolve[ode, {u, \[Rho]g, s, R}, {t, 0, 1}]]
Table[Plot[f /. sol, {t, 0, 1}, PlotRange -> All, Frame -> True,
PlotLabel -> f], {f, {u[t], R[t], \[Rho]g[t], s[t]}}]
出现了这样的错误NDSolve::ndsz: At t == 0.5757468138991848`, step size is effectively zero; singularity or stiff system suspected. >>
图片趋势是预想中的,但是如何能够画出完整的图呢,是不是由于参数设置的问题 导致算到t == 0.5757468138991848就进行不下去了,改了一天了 没啥进展,实在没办法,只好还向你请教下,实在不好意思~~·
还是或者因为我的方程是刚性的,这种情况该怎么办呢~
|
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|