- 积分
- 7
- 注册时间
- 2005-5-28
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2005-8-12 02:08:30
|
显示全部楼层
来自 浙江杭州
Re:推荐一款偏微分方程专解软件--Flexpde。
>方程为五维方程组。这个例子的源代码可以在 ...FlexPDE4\Samples\steady_state\fluids\buoyant中找到。对于一般偏微分方程,Flexpde提供一个有效的描述办法:
· TITLE – 用于在解图像中显示标题
· VARIABLES –描述变量
· DEFINITIONS –定义,可以在这里定义常量或者代换式,楼上附件中所给方程后两项就可以在这里定义,具体语句大家可参考源代码。代码中 “u = dy(psi)” 表示“定义u为psi对y的偏导数”,至于psi是几维到几维的映射Flexpde则会自己判断。
· EQUATIONS – 每一个方程都包括两部分,变量声明(例如为Temp)和含此变量的方程,含变量的方程一般来说应当将被声明的变量单独写在等号的左边,例如写做
“div(k*grad(temp)) = rho0*cp*[u*dx(temp) + v*dy(temp)]”
而不写做“div(k*grad(temp)) - rho0*cp*[u*dx(temp) - v*dy(temp)]=0”
以便快速构造刚度矩阵。但这不是必须的,只是影响速度和精度。
· BOUNDARIES – 边界条件,这里只需要学会定义解区间和定义边界条件。解区间定义是通过Star, Line to,Arc to等函数实现的,例如
"star(0,0)
...
Line to(1,0)
...
Line to(1,1)
...
Line to(0,1)
...
Line to finish
"即可通过定义四边形的四个顶点定义一个正方行。
详细定义方法可见Help
· MONITORS and PLOTS – 定义解的输出方式,例如:
grid(x,y)--画出解区间的网格
contour(temp) as "Temperature" painted--画出temp的彩色等势图,并以"Temperature" 为标题。
vector(u,v) as "Flow Velocity" norm--画出向量场(u,v)的箭头图,并以"Flow Velocity" 为标题。
等等
· END – 写完.
>接下来只要点击菜单栏中Run即可了,结果呢,大家自己看好了,我权限只能上传200kb图片,就没法为人民服务了。
>Flexpde和Matlab是有接口的,接口程序叫“FPMat”是国外一个牛人写的,有400kb我就不传了,可以在Math网站上搜到,懒得搜又想要的的给我发邮件好了(xueqilsj@163.com)。
光盘制作积分奖励 |
评分
-
1
查看全部评分
-
|