Difference between revisions of "FEA Mesh Generation Commands"

From OSUPDOCS
Jump to navigation Jump to search
 
(62 intermediate revisions by the same user not shown)
Line 1: Line 1:
The best way to generate an FEA mesh is to use [[Keypoint]], [[Path]], and [[Area]] commands, with an optional [[BMP, Region, and Hole]] commands to assign material properties from a digitized image of the object and from shape commands.
The best way to generate an FEA mesh is to use [[Keypoint Command|Keypoint]], [[Path Command|Path]], and [[Area Command|Area]] commands.
__TOC__
== Scripted Input Files ==


== Scripted Input Files ==
In scripted input files, you build the mesh by using multiple [[#mesh_cmds|mesh generation commands]]. These commands are usually grouped together, but may be at various locations in the input file. For example, a simple 2D, rectangular mesh of width <tt>#length</tt>, height <tt>#height</tt>, and thickness <tt>#thick</tt> with <tt>#horiz</tt> elements in the horizontal direction and <tt>#vert</tt> elements in the vertical direction (where terms beginning in number signs are [[Variable Names|script variables]]) can be created with:
 
Area #matname$,#thick
  Path "Bottom",#horiz
    Keypoint "BotLeft",0,0
    Keypoint "BotRight",#length,0
  EndPath
  Path "Right",#vert
    Keypoints "BotRight"
    Keypoint "TopRight",#length,#height
  EndPath
  Path "Top",#horiz
    Keypoints "TopRight"
    Keypoint "TopLeft",0,#height
  EndPath
  Path "Left",#vert
    Keypoints "TopLeft","BotRight"
  EndPath
EndArea
Resequence "BotRight"
 
<span id="mesh_cmds"></span>The following commands are used for generating the mesh in scripted files.


In scripted input files, you build the mesh by using multiple [[Mesh Generation Commands]]. These commands are usually grouped together, but may be at various locations in the input file.
* [[Area Command|Area]]: Define an area for mesh generation.
* [[Path Command|Path]]: Define a path used in defining [[Area Command|areas]].
* [[Path Command|Paths]]: Add multiple defined paths to an [[Area Command|area]].
* [[Keypoint Command|Keypoint]]: Define a keypoint used in defining [[Path Command|paths]]
* [[Keypoint Command|Keypoints]]: Add multiple defined keypoints to a [[Path Command|path]].
* [[Resequence Command]]: To turn on the option to resequence the nodes, which will minimize the bandwidth of the problem and speed the calculations.


== XML Input Files ==
== XML Input Files ==


In <tt>XML</tt> input files, all [[Mesh Generation Commands]] must be within a single Mesh block:
In <tt>XML</tt> input files, you build the mesh with [[#mesh_cmds_xml|mesh generation commands]] that must be within a single <tt><Mesh></tt> block. Here is a <tt><Mesh></tt> block that creates at 50X10 rectangular mesh with 6 elements in the horizontal direction and 3 in the vertical direction:


<pre>
<pre>
&lt;Mesh&gt;
&lt;Mesh&gt;
   &lt;<a href="#Keypoints"><code>Keypoints</code></a>&gt;
   &lt;Keypoints&gt;
     &lt;pt x='0' y='0' id='lowerleft'/&gt;
     &lt;pt x='0' y='0' id='lowerleft'/&gt;
     &lt;pt x='0' y='10' id='upperleft'/&gt;
     &lt;pt x='0' y='10' id='upperleft'/&gt;
Line 18: Line 46:
   &lt;/Keypoints&gt;
   &lt;/Keypoints&gt;


   &lt;<a href="#Paths"><code>Path</code></a> id="bottom" intervals='6' ratio='1'&gt;
   &lt;Path id="bottom" intervals='6' ratio='1'&gt;
     &lt;keypt id='lowerleft'/&gt;
     &lt;keypt id='lowerleft'/&gt;
     &lt;keypt id='lowerright'/&gt;
     &lt;keypt id='lowerright'/&gt;
Line 35: Line 63:
   &lt;/Path&gt;
   &lt;/Path&gt;


   &lt;<a href="#Areas"><code>Area</code></a> mat="1" thick="10" type="4" flip="0" angle="0"&gt;
   &lt;Area mat="1" thick="10" type="4" flip="0" angle="0"&gt;
     &lt;path id="bottom"/&gt;
     &lt;path id="bottom"/&gt;
     &lt;path id="right"/&gt;
     &lt;path id="right"/&gt;
Line 42: Line 70:
   &lt;/Area&gt;
   &lt;/Area&gt;
&lt;/Mesh&gt;
&lt;/Mesh&gt;
&lt;!-- image and shape methods to set element properties --&gt;
&lt;BMP name="MatImage.bmp" width="50" angles="MatAngles.bmp"&gt;
  &lt;Origin x="0" y="0"/&gt;
  &lt;Intensity mat="2" imin="76" imax="255"&gt;
    &lt;Thickness units="mm"&gt;1&lt;/Thickness&gt;
  &lt;/Intensity&gt;
  &lt;Intensity mat="3" imin="1" imax="75"&gt;
    &lt;Thickness units="mm"&gt;1&lt;/Thickness&gt;
  &lt;/Intensity&gt;
  &lt;Intensity imin="0" imax="255" minAngle="0.000000" maxAngle="90.000000"/&gt;
&lt;/BMP&gt;
&lt;Body mat='1' angle='0' thick='1'&gt;
  &lt;Oval units='mm' xmin='10' xmax='40' ymin='10' ymax='40'/&gt;
  &lt;Rect units='mm' xmin='10' xmax='40' ymin='10' ymax='40'/&gt;   
&lt;/Body&gt;
&lt;Hole&gt;
  &lt;Rect xmin='20' xmax='30' ymin='20' ymax='30'/&gt;   
&lt;/Hole&gt;
</pre>
</pre>


==== Keypoints ====
<span id="mesh_cmds_xml"></span>The fundamental <tt>XML</tt> commands for generating a mesh are:
 
The <tt>&lt;Keypoints&gt;</tt> section must come first. It can define any number of key points. Each subordinate <tt>&lt;pt&gt;</tt> command defines a key point, provides its <tt>x</tt> and <tt>y</tt> coordinates and gives it an <tt>id</tt>. The <tt>id</tt>'s can be text or numbers and must all be unique among key points.
 
==== Paths ====
 
Following the <tt>&lt;Keypoints&gt;</tt>, there are a series of <tt>&lt;Path&gt;</tt> commands that link the key points into lines or arcs and define mesh information along those paths. Subordinate to each <tt>&lt;Path&gt;</tt> command are a series of <tt>&lt;keypt&gt;</tt> commands that list the key points that define the path. Use two key points to define a line or three key points to define an arc. Currently the path can only have two or three key points. The key points are selected by their <tt>id</tt>. The attributes for a <tt>&lt;Path&gt;</tt> command are:
 
*<tt>id</tt>: An <tt>id</tt> for the path. The <tt>id</tt>'s can be text or numbers and must be unique among all paths.
*<tt>intervals</tt>: Number of elements along this path when it is meshed into an area.
* <tt>ratio</tt>: Ratio of the size of the first element along the path to the last one. If <tt>ratio<0</tt>, then the absolute size of the first element on the path with be <tt>|ratio|</tt>. The default is 1.0 which gives equally sized elements.
 
==== Areas ====
 
Finally, the FEA mesh is created by a series of <tt>&lt;Area&gt;</tt> commands. An <tt>&lt;Area&gt;</tt> command links a series of paths into a region to be meshed. Subordinate to each <tt>&lt;Area&gt;</tt> command are a series of <tt>&lt;path&gt;</tt> commands that list the paths that define the area. Currently [[NairnFEA]] only supports areas with four paths, to mesh quadrilateral areas, or areas with two paths, to mesh interfaces.
 
For quadrilateral areas, the end of the last path must connect to the start of the first path. The direction of each path does not matter as they will be rearranged if needed to define the area; but the paths must circumnavigate the area in the counter-clockwise direction. If <tt>n<sub>1</sub></tt>, <tt>n<sub>2</sub></tt>, <tt>n<sub>3</sub></tt>, and <tt>n<sub>4</sub></tt> are the number of nodes along each of the four paths, the meshing algorithm requires that <tt>n<sub>1</sub>+n<sub>2</sub>=n<sub>3</sub>+n<sub>4</sub></tt>.
 
For interface meshing, the two paths must be different, but indentical paths. In other words, they must connect different keypoints, but the coordinates of the keypoints, the number of intervals along the path, and the path ratios must be identical. The order and orientation of the two paths in the area command does not matter. Finally, interfaces only function when they connect two areas. To insure this connection, the paths for an interface must have each been previously used in exactly one quadrilateral area.
 
The attributes for a <tt>&lt;Area&gt;</tt> command are:


* <tt>mat</tt>: Material number in the list of defined materials. Quadrilateral areas cannot use an <a href="definematl.html#interface">interface material</a> while interface areas must use an <a href="definematl.html#interface">interface material</a>. When materials and angle will be set using a <code>BMP</code> file, the <code>mat</code> should be set to 0 to indicate no material assigned until the <a href="#bmpcmd"><code>&lt;BMP&gt;</code></a> command.
* [[Keypoint Command#XML Keypoints Block|<tt><Keypoints></tt>]]: Define all key point to be used in defining [[Path Command|paths]]
* <tt>matname</tt>: Alternatively, materials can be specied by name. See <a href="definematl.html#usemats">Defining Materials</a> for details on this method.
* [[Path Command#XML Path Command|<tt><Path></tt>]]: Define a path to be used in defining [[Area Command#XML Area Block|areas]].
* <tt>angle</tt>: Optional material angle relevant to <a href="definematl.html">anisotropic materials</a>. It can be entered as a number or as a <a href="function.html">user-defined function</a>. A function is evaluated at the midpoint of an element. Triangular elements are created in pairs and the function is evaluated at the midpoint of each pair.
* [[Area Command#XML Area Block|<tt><Area></tt>]]: Define an area for mesh generation.
* <tt>thick</tt>: Thickness in mm for plane stress or plain strain analysis; not needed for axisymmetric analysis.
* [[Resequence Command]]: To turn on the option to resequence the nodes, which will minimize the bandwidth of the problem and speed the calculations.
* <tt>type</tt>: The type of element to create. The options are:
<ul class="lev2args">
<li><code>1</code>: 3-node constant-stain triangular elements (linear)</li>
<li><code>2</code>: 4-node isoparametric rectangular elements(linear)</li>
<li><code>3</code>: 8-node isoparametric rectangular elements (quadratic)</li>
<li><code>4</code>: 6-node isoparametric triangular elements (quadratic)</li>
<li><code>5</code>: 4-node imperfect interface elements (linear)</li>
<li><code>6</code>: 6-node imperfect interface elements (quadratic)</li>
<li><code>8</code>: 9-node Langrangian rectangular element (quadratic) (9 point Gaussian quadrature)


The first <code>&lt;Area&gt;</code> command must define an element type. Subsequent <code>&lt;Area&gt;</code> commands can omit the element type unless the element type is changed to have an analysis with mixed elements. When mixed elements are used, they must be compatible elements (<i>i.e.</i>, you can mix linear elements (types 1, 2, and 5) with each other or quadratic elements (types 3, 4, 6, and 8) with each other, but cannot mix linear elements with quadratic elements). Whenever a <code>type</code> attribute is used, it will apply to that area and all subsequent areas until it is changed. Finally, quadrilateral areas can only use solid elements (types 1, 2, 3, 4, and 8) and interface areas can use only interface elements (type 5 and 6).</li>
== Using Images and Shapes ==
<li><code>flip</code>: For triangular elements, the initially created quadrilaterials are split into two elements by drawing a diagonal across the quadrilaterial. The <code>flip</code> option can change which diagonal is used. It should be 0 or 1. Once it is set, it will apply to that area and all subsequent areas until it is changed.</li>
</ul>


== Mesh Generation Commands ==
If some elements in the generated mesh are set the material id _NONE_ (see [[Area Command]]) or mat = 0 (see [[Area Command#XML Area Block|<tt><Area></tt> Command]]), you can follow mesh creation with [[Create a Mesh with Images and Shapes|image and shape commands]] to fill in material types, angles, and define the final mesh geometry. This option is described [[Create a Mesh with Images and Shapes|here]].

Latest revision as of 13:17, 2 June 2015

The best way to generate an FEA mesh is to use Keypoint, Path, and Area commands.

Scripted Input Files

In scripted input files, you build the mesh by using multiple mesh generation commands. These commands are usually grouped together, but may be at various locations in the input file. For example, a simple 2D, rectangular mesh of width #length, height #height, and thickness #thick with #horiz elements in the horizontal direction and #vert elements in the vertical direction (where terms beginning in number signs are script variables) can be created with: 

Area #matname$,#thick
  Path "Bottom",#horiz
    Keypoint "BotLeft",0,0
    Keypoint "BotRight",#length,0
  EndPath
  Path "Right",#vert
    Keypoints "BotRight"
    Keypoint "TopRight",#length,#height
  EndPath
  Path "Top",#horiz
    Keypoints "TopRight"
    Keypoint "TopLeft",0,#height
  EndPath
  Path "Left",#vert
    Keypoints "TopLeft","BotRight"
  EndPath
EndArea
Resequence "BotRight"

The following commands are used for generating the mesh in scripted files.

  • Area: Define an area for mesh generation.
  • Path: Define a path used in defining areas.
  • Paths: Add multiple defined paths to an area.
  • Keypoint: Define a keypoint used in defining paths
  • Keypoints: Add multiple defined keypoints to a path.
  • Resequence Command: To turn on the option to resequence the nodes, which will minimize the bandwidth of the problem and speed the calculations.

XML Input Files

In XML input files, you build the mesh with mesh generation commands that must be within a single <Mesh> block. Here is a <Mesh> block that creates at 50X10 rectangular mesh with 6 elements in the horizontal direction and 3 in the vertical direction: 

<Mesh>
  <Keypoints>
    <pt x='0' y='0' id='lowerleft'/>
    <pt x='0' y='10' id='upperleft'/>
    <pt x='50' y='0' id='lowerright'/>
    <pt x='50' y='10' id='upperright'/>
  </Keypoints>

  <Path id="bottom" intervals='6' ratio='1'>
    <keypt id='lowerleft'/>
    <keypt id='lowerright'/>
  </Path>
  <Path id="right" intervals='3'>
    <keypt id='lowerright'/>
    <keypt id='upperright'/>
  </Path>
  <Path id="top" intervals='6'>
    <keypt id='upperright'/>
    <keypt id='upperleft'/>
  </Path>
  <Path id="left" intervals='3'>
    <keypt id='upperleft'/>
    <keypt id='lowerleft'/>
  </Path>

  <Area mat="1" thick="10" type="4" flip="0" angle="0">
    <path id="bottom"/>
    <path id="right"/>
    <path id="top"/>
    <path id="left"/>
  </Area>
</Mesh>

The fundamental XML commands for generating a mesh are:

  • <Keypoints>: Define all key point to be used in defining paths
  • <Path>: Define a path to be used in defining areas.
  • <Area>: Define an area for mesh generation.
  • Resequence Command: To turn on the option to resequence the nodes, which will minimize the bandwidth of the problem and speed the calculations.

Using Images and Shapes

If some elements in the generated mesh are set the material id _NONE_ (see Area Command) or mat = 0 (see <Area> Command), you can follow mesh creation with image and shape commands to fill in material types, angles, and define the final mesh geometry. This option is described here.

Retrieved from "http://osupdocs.forestry.oregonstate.edu/index.php?title=FEA_Mesh_Generation_Commands&oldid=4102"