matlab中PETlb的应用.doc

MATLAB R2009a\help\product heIp\Partial Differential Equation Toolbox\Graphical User Interface this chapter discusses the graphical user interface (GUI) pdetooL The main components of the GUI are the menus，the dialog boxes, and the toolbar. PDE Toolbox 求解椭 PDE Toolbox 求解椭 、抛物、双曲方程的基本步骤 第一步： 在MATLAB命令窗口，输入命令 ? pedtool 进入 PDE Toolbox 窗口 第二步：建立几何模型 第二步：建立几何模型 丨□田CD①》 在toolbar中选择几何图形，如建立单位圆，则点击第3或第4个按钮，用鼠标的右键，click-and-drag， 创立一个圆cl 双点击圆cl，弹出dialog box，输入几何参数 第三步：输入边界条件 点击 输入边界条件。 若所有边界条件都是齐次的第一类边界条件，此步可省略。 第四步：输入偏微分方程参数 点击 PDE 输入偏微分方程参数。 93 liWBlJH 蠹ou 讀 C i W 1 — 第五步：划分网格 点击 ，?-\ 一次，初步划分网格。 点击 多次，细划网格。 单元圆细化两次的结果 第六步：输入初始条件 若方程是椭圆的，此步可省略。 若方程是抛物的，则需输入初始位移条件及计算参数若方程是双曲的，则需输入初始位移条件、初始速度条件及计算参数 若方程是抛物的，则需输入初始位移条件及计算参数 若方程是双曲的，则需输入初始位移条件、初始速度条件及计算参数 第七步：求解有限元方程 点击 = 求解有限元方程。 第八步：绘图 点击 绘图。 利用绘图选项，将计算结果进行可视化输出。 下面是两个编程算例。 泊松方程边值问题 单位圆域内泊松方程齐次边值问题 )=1， X2 + y2 1 w = 0, x2 + y2 = \ 其精确解为 M = i(l_x2_y2) 试求泊松方程问题的数值解并与精确解比较。 解(1)建立有限元模型 在MATLAB命令窗口，输入命令 ? pedtool 进入 PDE Toolbox 銜口 11 o m石运 在toolbar中点击第3或第4个按钮，用鼠标的右键，click-and-drag,创立一个岡cl 双点击圆cl,输入几何参数 Lol 回 Lol 回 0-1 Q Object Dialog CreteObject type Crete X-ccntcr Y-c?nter Radus Nome OK Cancel 输入边界条件（略）。 因为所有边界条件都是齐次的第一类边界条件。 输入偏微分方程参数。 点击 PDE 输入偏微分方程参数。 划分网格。点击 划分网格。 点击 I 一次，初步划分网格。 点击 1 2次，细划网格。 结点数 结点数2097 （2）求解有限元模型并输出计算结果点击点击然后,求解有限元方程。,选择绘图参数 （2）求解有限元模型并输出计算结果 点击 点击 然后, 求解有限元方程。 ,选择绘图参数 绘图（plot） 数值解的最大值为0.25,与精确解相同。 注：程序 unitcircle Lm 波动方程初边值问题 方形域波动方程初边值问题 utt — Aw = 0,71Lo 二 arctan(cos —X) utt — Aw = 0, 71 Lo 二 arctan(cos —X)， 71 3 sin 7Cx exp(sin —)，)， 0, -1 x 1，-1 )， 1,0 z 5 -lxl-lyl 一1x1，一 l_yl 0z5 du dn 0, 0z5 试求波动方程问题的数值解。 解（1）建立有限元模型 在MATLAB命令窗口，输入命令 ? pedtool 进入 PDE Toolbox 窗 U 在toolbar中点击第1或第2个按钮，用鼠标的右键，click-and-drag,创立一，个A形cl 双点击圆SQ1,输入几何参数 双点击圆SQ1,输入几何参数 输入边界条件。 点击 _bn 输入边界条件。 左右两条边的边界条件都是齐次的第一类边界条件（红色线条），可省略设置。 分别双点击上边和下边，选择齐次的第二类边界条件。 进行网 进行网 点击输入偏微分方程参数， 点击 输入偏微分方程参数， 网格划分（结点数10145） 输入初始条件。 点击菜单 Solve 中的^arameters选项。 方程是双曲的，则需输入初始位移条件、初始速度条件及计算参数 H Solve Farame*efs I― （2）求解有限元模型并输出计算结果 点击 求解有限