有限元的MATLAB解法.pdfVIP

有限元的MATLAB解法

1。打开MATLAB。

2。输入“pdetool”再回车，会跳出PDEToolbox的窗口(PDE意为

偏微分方程,是partialdifferentialequations的缩写），需要的话可点击

Options菜单下Grid命令，打开栅格.

3。完成平面几何模型：在PDEToolbox的窗口中，点击工具栏下的

矩形几何模型进行制作模型，可画矩形R，椭圆E，圆C，然后在

Setformula栏进行编辑并（如双脊波导R1+R2+R3改为RI—R2—R3，

设定a、b、s/a、d/b的值从而方便下步设定坐标）

用算术运算符将图形对象名称连接起来，若还需要，可进行储存，

形成M文件。

4.用左键双击矩形进行坐标设置：将大的矩形left和bottom都设为

0，width是矩形波导的X轴的长度，height是矩形波导的y轴的长

度,以大的矩形左下角点为原点坐标为参考设置其他矩形坐标。

5。进行边界设置：点击“Boundary”中的“BoundaryMode”,再点

5

击“Boundary”中的“SpecifyBoundaryConditions”，选择符合的边

界条件,Neumann为诺曼条件,Dirichlet为狄利克雷条件,边界颜色显

示为红色.

6。进入PDE模式：点击点击菜单下“PDEMode”命令，进入PDE

模式，单击“PDESpecification”,设置方程类型,“Elliptic”为椭圆

型，“Parabolic”为抛物型，““为双曲型，““

为特征值问题.

7。对模型进行剖分:点击“Mesh”中““进行初次剖分,

若要剖的更细,再点击“RefineMesh”进行网格加密.

8。进行计算：点击“Solve”中“SolvePDE”,解偏微分方程并显示

图形解,u值即为Hz或者Ez。

9。单击“Plot”菜单下“Parameters”选项，打开“PlotPlot

对话框.选中Color,Height(3-Dplot)和Showmesh三项，然后单击

“Plot”按钮，显示三维图形解。

10.如果要画等值线图和矢量场图,单击“Plot”菜单下“Parameters”

选项,打开“PlotSelection”对话框。选中Contour和Arrows两项，

然后单击Plot按钮，可显示解的等值线图和矢量场图.

11。将计算结果条件和边界导入MATLAB中：点击“ExportSolution”,

再点击“Mesh”中“ExportMesh”。

6

12。在MATLAB中将编好的计算程序导入，按F5运行.

备注：

Property（属性）用于画图时选用相应的绘图类型

u方程的解

abs（grad（u)）每个三角形的中心的▽u的绝对值

abs(c*grad（u)）每个三角形的中心的c·▽u的绝对值

—grad（u）u的负梯度-▽u

我们也可以用MATLAB程序求解PDE问题,同时显示解的图形；

7

一个长直接接地金属矩形槽,其侧壁与底面电位均为0,顶盖电

位为100V，求槽内的电位分布：

100V

0V0V