Case02 Mesh Regions to define element size

 

In the present case the logic behind mesh regions will be demonstrated. Mesh_region data structure may be used to define regions in the space with a target element size. Then all geometry points of the domain located within such mesh region will be assigned such element size. Two cases with identical mesh generation data including:

 

1. Default element size of 80 m

2.A "rectangular" Mesh_region (a rectangular prism in 3D) located in one of the corners of the cube geometry with a target element size of 20 m

 

Will be used to demonstrate Mesh_region usage and limitations (note that the refereed method for defining mesh with regions of different element sizes is via adaptivity) . The two cases provided in Mesh_002\Case02\Data differ in the geometry definition data as follows:

 

Mesh_002_Case02a.dat : Same cube as in previous cases, defined by 8 nodes, 12 lines, 6 surfaces and 1 volume

Mesh_002_Case02b.dat : Similar cube where 3 extra nodes and 3 extra lines have been defined which are located within the defined mesh region in the corner of the cube (see figure below).

 

Mesh_002_Fig05

Geometry for the two cases with the nodes defining the geometry

 

 

Below a description of the data files is provided. It should be noted that only the data relative to mesh generation will be described in detail while for the remaining data a brief description is provided.

 

Datafile Description

 

The basic data includes:

 

1.Geometry definition data (Nodal_data , Geometry_line , Geometry_surface and Geometry_volume ) at the end of the file after the END DATA command to define the cube geometry.

2.Group_data and Group_control_data for a single group with TET4V elements (note that the Element_type defined must be valid for the type of mesh being generated)

3.Material_data and Fluid_properties defining an elastic material

4.Support_data defining displacement fixities in the 3 directions for the whole cube volume (by assigning the corresponding Displacement_codes set to the internal Geometry_set "All_volumes" )

5.Mesh generation data (Mesh_control_data, Unstructured_mesh_data and Mesh_region) with identical data for both cases.

6.Control_data defining the Duration of the simulation stage, the Target_number_time_steps for the stage as well as several output options.

 

 

Mesh Data

 

Data File


 

* Mesh_control_data

! ---------------------------------

 Generation_algorithm    2          

 Mesh_generation_flag    0

 

 

* Unstructured_mesh_data

! ---------------------------------

 Default_element_size         80.0

 Element_size_bounds IDM=2

  /Max. Element Size/        80.0

  /Min. Element Size/       20.0

 Mesh_regions  IDM=1

   1

 

* Mesh_region  NUM=1

! ---------------------------------

 Region_type     "Rectangular"

 Element_size    20.0

 Region_size  IDM=6

  /X-min/    -100.0

  /Y-min/    -100.0

  /Z-min/    600.0

  /X-max/    400.0

  /Y-max/    400.0

  /Z-max/    1100.0

 

1.Mesh_control_data defines the options for mesh generation. In this case this has been set to:

a.Generation_algorithm is set to 2 (unstructured mesh)

b.Mesh_generation_flag is set to 0 in order to perform the analysis. Note that in this case it would not be required as we are only interested in generating the mesh.

 

2.Unstructured_mesh_data with:

a.A default element size of 80 m.

b.Maximum and minimum mesh size bounds set to 80.0 m and 20.0 m respectively

c.Assignment of the Mesh_region with ID number NUM=1 via Mesh_regions keyword

 

3.The Mesh_region data is defined as:

a.The Region_type is set as "Rectangular" which in 3D corresponds to a Box

b.The target Element_size within the region has been set to 20 m

c.Region_size is used to define the topology / coordinates of the mesh region. The data required is dependent on the Region_type. In the present case the region is rectangular so the data required comprises definition of the two corner points diagonal in the box with minimum and maximum X, Y and Z coordinates respectively (see picture below).

 

Mesh_002_Fig07

Mesh region considered in this case shown with the cube geometry

 

 

 

 

 

Results

 

The results for the two cases are provided in: Mesh_002\Case02\Results. In the figure below it can be seen that in case02a the mesh region only has an effect in prescribing an element size of 20 m in the geometry point at (x, y, z) = (0, 0, 1000). In case02b the three extra points / nodes used to define the geometry are also located within the mesh region and therefore will have a prescribed element size of 20 m. Consequently the lines joining the points within the mesh region have a consistent element size of 20 m as indicated by the red arrows.

 

As may be anticipated, defining a mesh with different target element sizes at different parts of the domain may be cumbersome via Mesh_region. Thus a more flexible and easy way to achieve that is by remeshing the initial meshed geometry via adaptivity as demonstrated in Case03.

 

Mesh_002_Fig06

Results for the two cases with the mesh region in orange and the geometry points affected by the mesh region in green. The unaffected points are shown in grey