1、 pedtool进入PDE Toolbox 窗口第二步:建立几何模型在toolbar中选择几何图形,如建立单位圆,则点击第3或第4个按钮,用鼠标的右键,click-and-drag,创立一个圆c1双点击圆c1,弹出dialog box,输入几何参数第三步:输入边界条件点击 输入边界条件。若所有边界条件都是齐次的第一类边界条件,此步可省略。第四步:输入偏微分方程参数点击 输入偏微分方程参数。第五步:划分网格点击 一次,初步划分网格。点击 多次,细划网格。 单元圆细化两次的结果第六步:输入初始条件若方程是椭圆的,此步可省略。点击菜单 中的 Parameters 选项。若方程是抛物的,则需输入初始位
2、移条件及计算参数若方程是双曲的,则需输入初始位移条件、初始速度条件及计算参数第七步:求解有限元方程点击 求解有限元方程。第八步:绘图点击 绘图。利用绘图选项,将计算结果进行可视化输出。下面是两个编程算例。泊松方程边值问题单位圆域内泊松方程齐次边值问题其精确解为试求泊松方程问题的数值解并与精确解比较。解 (1) 建立有限元模型在toolbar中点击第3或第4个按钮,用鼠标的右键,click-and-drag,创立一个圆c1双点击圆c1,输入几何参数输入边界条件(略)。因为所有边界条件都是齐次的第一类边界条件。输入偏微分方程参数。划分网格。点击 2次,细划网格。 结点数2097(2)求解有限元模型
3、并输出计算结果点击 ,选择绘图参数然后,绘图(plot)数值解的最大值为0.25,与精确解相同。注:程序 unitcircle_1.m波动方程初边值问题方形域波动方程初边值问题 试求波动方程问题的数值解。在toolbar中点击第1或第2个按钮,用鼠标的右键,click-and-drag,创立一个方形c1双点击圆SQ1,输入几何参数输入边界条件。左右两条边的边界条件都是齐次的第一类边界条件(红色线条),可省略设置。分别双点击上边和下边,选择齐次的第二类边界条件。进行网格划分。 网格划分(结点数10145)输入初始条件。方程是双曲的,则需输入初始位移条件、初始速度条件及计算参数若欲进行动画演示,则在绘图选项中,选择 Animation程序 wavesquare.m