偏微分方程数值解的matlab 实现,偏微分方程数值解的Matlab 实现

【实例简介】

工程中有许多问题可以归结为偏微分方程问题，如弹塑性力学中研究对象(结构、边坡等)内部的应力应变问题、地下水渗流问题等。这些由偏微分方程及边界条件、初始条件等组合成的数学模型，只有在十分特殊的条件下才能求得解析解。因此，在很长一段时间内，人们对于这一类问题是无能为力的。随着计算机技术的发展，各种数值方法应运而生，如有限元法、有限差分法、离散元法、拉格朗日元法等等。利用数值法，可以求得这些问题的数值解。它不是问题的精确解，但可以无限接近精确解。Matlab采用有限元法求解偏微分方程的数值解

偏微吩方程数宽解的 Matlab罢砚

foxdog制作

中输入公式:C1-C2。

3.3定义边界条件

在 Boundary菜单中选择 Boundary mode选项或直接在工具条中单击a按钮,则显示几何模型的外

边界和内边界。如图24所示。

回区

-1;

11l=3,re,命,=|÷正m=:2-1

〔!1〔2〕001〔E

图24定义边界条件

在图中单击边界,则边界线变成黑色,表示它被选中。按住 Shift键,可以连续选择多条边界线段。

选择要定义的边界以后,双击边界或在 Specify Boundary conditions…近项,打开 Boundary Condition

对话框,如图25所示。

吧ndar丁 Condition

瓦[

F IrI dNrs I: ' ll 1(u-I

cue

11

IuI-H

图25 Boundary Condition对话框

对话框中各选项的意义分别为:

Boundary cond. cquation标签显示边界条件方程。

Condition type在下面的两个单选钮中进行选择,确定边界条件类型。

Neumann单选钮选择此项,将边界条件类型确定为 Neuman条仵。

Dirichlet单选钮缺省时选此项。选择此项时,将边界条件类型确定为

Dirichlet边界条件。

● Coelficient value栏在该栏中的对应文丕框中输入边界条件公式卬的系数值。

本例中选择缺省设置。

在 Boundary菜单中近择 Show Edge labels选项和 Show Subdomain labels选项,可以在图中显示

边缘标签和子域标签。

3.4定义PDE类型和PDE系数

偏微分方覆憝宽解的 Matlab寥砚

foxdog制作

在PDE菜单中选择 PdE mode选项,图形将显示为PE模式。选择 Show subdomain labels选项,将

在图中显示子域标签。

在工具条中单击r按钮或在PDE菜单中单击 PDE Specification…选项,可以打开PDE

Specification对话框,如图26所丌。

E