# Generate Mesh

The `generateMesh`

function creates a triangular mesh for a 2-D geometry and a tetrahedral mesh for a 3-D geometry. By default, the mesh generator uses internal algorithms to choose suitable sizing parameters for a particular geometry. You also can use additional arguments to specify the following parameters explicitly:

Target maximum mesh edge length, which is an approximate upper bound on the mesh edge lengths. Note that occasionally, some elements can have edges longer than this parameter.

Target minimum mesh edge length, which is an approximate lower bound on the mesh edge lengths. Note that occasionally, some elements can have edges shorter than this parameter.

Mesh growth rate, which is the rate at which the mesh size increases away from the small parts of the geometry. The value must be between 1 and 2. This ratio corresponds to the edge length of two successive elements. The default value is 1.5, that is, the mesh size increases by 50%.

Quadratic or linear geometric order. A quadratic element has nodes at its corners and edge centers, while a linear element has nodes only at its corners.

Create a model and include the geometry of a plate with a square hole in its center.

`model = femodel(Geometry="PlateSquareHolePlanar.stl");`

Plot the geometry.

pdegplot(model)

Generate a default mesh. For this geometry, the default target maximum and minimum mesh edge lengths are 8.9443 and 4.4721, respectively.

model = generateMesh(model); model.Geometry.Mesh

ans = FEMesh with properties: Nodes: [2x1210 double] Elements: [6x570 double] MaxElementSize: 8.9443 MinElementSize: 4.4721 MeshGradation: 1.5000 GeometricOrder: 'quadratic'

View the mesh.

figure pdemesh(model)

For comparison, create a mesh with the target maximum element edge length of 20.

model = generateMesh(model,Hmax=20); model.Geometry.Mesh

ans = FEMesh with properties: Nodes: [2x298 double] Elements: [6x132 double] MaxElementSize: 20 MinElementSize: 10 MeshGradation: 1.5000 GeometricOrder: 'quadratic'

figure pdemesh(model)

Now create a mesh with the target minimum element edge length of 0.1.

model = generateMesh(model,Hmin=0.1); model.Geometry.Mesh

ans = FEMesh with properties: Nodes: [2x1390 double] Elements: [6x660 double] MaxElementSize: 8.9443 MinElementSize: 0.1000 MeshGradation: 1.5000 GeometricOrder: 'quadratic'

figure pdemesh(model)

Create a mesh, specifying both the maximum and minimum element edge lengths instead of using the default values.

model = generateMesh(model,Hmax=20,Hmin=0.1); model.Geometry.Mesh

ans = FEMesh with properties: Nodes: [2x450 double] Elements: [6x208 double] MaxElementSize: 20 MinElementSize: 0.1000 MeshGradation: 1.5000 GeometricOrder: 'quadratic'

View the mesh.

figure pdemesh(model)

Create a mesh with the same maximum and minimum element edge lengths, but with the growth rate of 1.9 instead of the default value of 1.5.

model = generateMesh(model,Hmax=20, ... Hmin=0.1, ... Hgrad=1.9); model.Geometry.Mesh

ans = FEMesh with properties: Nodes: [2x374 double] Elements: [6x170 double] MaxElementSize: 20 MinElementSize: 0.1000 MeshGradation: 1.9000 GeometricOrder: 'quadratic'

figure pdemesh(model)

You also can choose the geometric order of the mesh. The toolbox can generate meshes made up of quadratic or linear elements. By default, it uses quadratic meshes, which have nodes at both the edge centers and corner nodes.

model = generateMesh(model,Hmax=50); figure pdemesh(model,NodeLabels="on") hold on plot(model.Geometry.Mesh.Nodes(1,:), ... model.Geometry.Mesh.Nodes(2,:), ... "ok","MarkerFaceColor","g")

To save memory, override the default quadratic geometric order.

model = generateMesh(model,Hmax=50, ... GeometricOrder="linear"); figure pdemesh(model,NodeLabels="on") hold on plot(model.Geometry.Mesh.Nodes(1,:), ... model.Geometry.Mesh.Nodes(2,:), ... "ok","MarkerFaceColor","g")