Code_Aster
®
Version
8.2
Titrate:
Macro-control
MACR_INFO_MAIL
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.03.02-D1
Page:
1/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
Organization (S):
EDF-R & D/SINETICS
Instruction manual
U7.0- booklet: Data exchanges
Document: U7.03.02
Macro-control
MACR_INFO_MAIL
1 Goal
To give information on the quality of a mesh.
This macro-control makes it possible to obtain information on a mesh, of degree 1 or 2, container
mesh-points, segments, triangles, quadrangles or tetrahedrons. The impression is
directed on the file of messages.
One finds initially an assessment total of the mesh: dimension, degree, a many nodes, number
extreme elements, dimensions.
Then four information is accessible in an optional way:
·
a summary on the quality of the elements,
·
a control of nonthe interpenetration of the elements,
·
an assessment on the connexity of the mesh,
·
a summary of the size of different the under-field.
This analysis is made by the software of adaptation of mesh LOBSTER. For more details, to see:
http://www.code-aster.org/outils/homard or to refer to the documents quoted in
bibliography.
Code_Aster
®
Version
8.2
Titrate:
Macro-control
MACR_INFO_MAIL
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.03.02-D1
Page:
2/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
2 Syntax
MACR_INFO_MAIL (
MESH
=
my
[mesh]
NON_SIMPLEXE =/
0
/1
/
2
[DEFECT]
NUMBERS
=/
“YES”
[DEFECT]
/
“NOT”
QUALITY =/
“YES”
[DEFECT]
/
“NOT”
CONNEXITY =/
“YES”
[DEFECT]
/
“NOT”
CUT
=/
“YES”
[DEFECT]
/
“NOT”
INTERPENETRATION
=
/
“YES”
[DEFECT]
/
“NOT”
MAILLAGE_FRONTIERE = maf
[mesh]
GROUP_MA
=
l_grma [l_gr_maille]
LANGUAGE
=/
“FRENCH” [DEFECT]
“FRENCH”
“ENGLISH”
“ENGLISH”
VERSION_HOMARD
=
/
“V7_7”
[DEFECT]
“V7_N'
“V7_N_PERSO”
INFORMATION =
/
1
[DEFECT]
/2
),
Code_Aster
®
Version
8.2
Titrate:
Macro-control
MACR_INFO_MAIL
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.03.02-D1
Page:
3/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
3 Operands
3.1 Operand
MESH
MESH
= my
Mesh of the type [mesh] to analyze. Attention, the analysis will relate only to the nodes, them
mesh-points, segments, triangles, quadrangles or tetrahedrons, in degree 1 or 2. If
one provides a comprising mesh of other elements, two cases of figure are possible: that is to say one
stop in error, is an adaptation on the area in simplex and a restitution with identical of the remainder
mesh. The choice between these two operating modes is made by the key word
NON_SIMPLEXE.
3.2 Operand
NON_SIMPLEXE
NON_SIMPLEXE
=/0
/1
/
2
[DEFECT]
In its current version, LOBSTER can read all the types of elements but does not make carry the analysis
that on some: mesh-points, segments, triangles, quadrangles and tetrahedrons.
By retaining option 0, the transmission of a mesh containing another thing that simplexes
a stop in error will involve.
By choosing option 1, one will be able to transmit a mesh comprising any type
of element. The analysis will relate only to the area in simplexes or quadrangles.
With option 2, one authorizes the analysis only for one mesh which contains exclusively
elements which LOBSTER can treat: segments, triangles, quadrangles and tetrahedrons. The presence
hexahedrons or pentahedrons will cause a stop in error. It is the default option.
3.3 Operand
NUMBERS
NUMBERS
=
/
“YES”
[DEFECT]
/
“NOT”
If the choice is “NOT”, nothing occurs.
If the choice is “YES”, an assessment of the numbers of nodes and elements are printed on the file of
messages.
3.4 Operand
QUALITY
QUALITY
=
/
“YES”
[DEFECT]
/
“NOT”
If the choice is “NOT”, nothing occurs.
If the choice is “YES”, an assessment of the quality of the elements is printed on the file of message.
The quality of a triangle is defined as being the relationship between the length on the largest side and
the radius of the inscribed circle. The quality of a quadrangle is defined like the quotient of the product of
the biggest length and of the averages on the sides and the diagonals by smallest of
surfaces of the triangles intern with the quadrangles. In the same way, the quality of a tetrahedron is defined
as being the relationship between the length on the largest side and the radius of the registered sphere. These
reports/ratios are standardized to be worth 1 in the case of an equilateral triangle, of a square, or one
equilateral tetrahedron. For any nonequilateral element, quality is higher than 1. See
reference [bib1] for detailed explanations.
The result is presented in the form of tables, with the extreme values.
Code_Aster
®
Version
8.2
Titrate:
Macro-control
MACR_INFO_MAIL
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.03.02-D1
Page:
4/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
3.5 Operand
INTERPENETRATION
INTERPENETRATION
=
/
“YES”
[DEFECT]
/
“NOT”
If the choice is “NOT”, nothing occurs.
If the choice is “YES”, it is checked that the mesh is correct from the point of view of the recovery:
no triangle is with horse on another triangle, no quadrangle is not with horse on another
quadrangle, no tetrahedron enters another tetrahedron.
3.6 Operand
CONNEXITY
CONNEXITY
=
/
“YES”
[DEFECT]
/
“NOT”
If the choice is “NOT”, nothing occurs.
If the choice is “YES”, an assessment of the connexities is printed on the file of messages. One specifies
then if the segments, the elements 2D (joined together triangles and quadrangles) or the tetrahedrons are
of only one holding or divided into several blocks.
3.7 Operand
CUT
CUT
=
/
“YES”
[DEFECT]
/
“NOT”
If the choice is “NOT”, nothing occurs.
If the choice is “YES”, an assessment of the sizes of the under-fields is printed on the file of
messages. A under-field is defined like a whole of of the same meshs dimension and
belonging to the same groups.
3.8 Operand
MAILLAGE_FRONTIERE
MAILLAGE_FRONTIERE = maf
In dimension 2, the choice of this option means that the segments forming the edge of the mesh of
calculation are attached to a fine description of this edge. This fine mesh of the edge is transmitted here.
link is done by membership of the segments to the same groups.
3.9 Operand
GROUP_MA
GROUP_MA = l_grma
If this option misses, the link between the segments is established for all the groups present in
mesh of the border. If the link is established only for some groups, they are indicated here.
3.10 Operand
LANGUAGE
LANGUAGE
=
/
“FRENCH” [DEFECT]
“FRENCH”
“ENGLISH”
“ENGLISH”
This operand specifies the language in which the messages resulting from LOBSTER are printed.
Code_Aster
®
Version
8.2
Titrate:
Macro-control
MACR_INFO_MAIL
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.03.02-D1
Page:
5/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
3.11 Operand
VERSION_HOMARD
VERSION_HOMARD
=
“V7_7”
[DEFECT]
“V7_N'
“V7_N_PERSO”
This operand makes it possible to select the LOBSTER version which is used for this analysis.
By defect, LOBSTER 7.7 is launched. It is the version of reference. Choice “V7_N' activates
version 7.n of LOBSTER which is the version of development. Active choice “V7_N_PERSO”
a version of development specific to the user. This option in fact is reserved for the team of
development of LOBSTER to develop new functionalities.
3.12 Operand
INFORMATION
INFORMATION
=
/1
/2
If INFORMATION is worth 2, the entirety of the LOBSTER exit is included in the file of message.
If not, nothing in particular takes place.
Code_Aster
®
Version
8.2
Titrate:
Macro-control
MACR_INFO_MAIL
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.03.02-D1
Page:
6/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
4 Example
MACR_INFO_MAIL (MESH = HAMMER,
CUT
=
“YES”,
CONNEXITY
=
“YES”,
INTERPENETRATION
=
“YES”,
QUALITY =
“YES”
)
This sequence will write information general (dimension, degree, a number of nodes and
elements,…) then a diagnosis on the quality and the interpenetration of the elements used in
mesh
HAMMER
as well as the size and the connexity of the various under-fields. The version of
LOBSTER by defect is used.
ANALYZE MESH
===================
Mesh has to analyze
Creation date: Thursday September 2, 2004 has 9. 57 mn 11 S
Dimension: 2
Degree: 1
It is a starting mesh.
Direction | Unit | Minimum | Maximum
-------------------------------------------------------------
X | Unknown factor | 0.0000 | 20.000
y | Unknown factor | 0.0000 | 11.000
The diagnoses on the interpenetrations of elements points out the precision retained for the control of
coincidence of the nodes.
INTERPENETRATION OF THE ELEMENTS
=============================
… Coincidence of the nodes
--------------------------
Two nodes are known as coincidents if the absolute deviation
between their co-ordinates is lower a: 0.10E-39
… Interpenetration of the edges
-------------------------------
… Interpenetration of the triangles
----------------------------------
**********************************************************
**
* No problem was meets. *
**
**********************************************************
Code_Aster
®
Version
8.2
Titrate:
Macro-control
MACR_INFO_MAIL
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.03.02-D1
Page:
7/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
The diagnosis on the quality of the elements points out initially the extreme values recorded in
analyzed mesh. Here the minimum is 1,0044, very near to the absolute minimum theoretical which is 1,
and the maximum met is 1,5788. Then one presents the distribution by equidistant section to
to leave the optimum value, 1. It is seen that 16 triangles have a quality ranging between 1 and 1,05, that is to say
13,56% of the total number of triangles. 26 triangles have a quality ranging between 1,05 and 1,10, are
22,03% of the total number of triangles. In cumulated, one thus notes that 42 (=16+26) triangles have one
quality better than 1,10, is 35,59% of the total. And so on. For example, 98,31% of the elements
have a quality better than 1,45.
QUALITY OF THE ELEMENTS
====================
**********************************************************
* Quality of the triangles of the mesh of calculation *
* Recall: quality is equal to the report/ratio of the diameter *
* of the triangle on the radius of the inscribed circle, *
* standardizes has 1 for a regular triangle. *
**********************************************************
* Minimum: 1.0044 Maximum: 1.5788 *
**********************************************************
**********************************************************
* Function of distribution *
**
* Values * a Number of elements *
* Minicomputer < < Maximum * by class * office plurality *
** in %. numbers * in %. numbers *
**********************************************************
* 1.00 < 1.05 * 13.56. 16 * 13.56. 16 *
* 1.05 < 1.10 * 22.03. 26 * 35.59. 42 *
* 1.10 < 1.15 * 22.03. 26 * 57.63. 68 *
* 1.15 < 1.20 * 18.64. 22 * 76.27. 90 *
* 1.20 < 1.25 * 9.32. 11 * 85.59. 101 *
* 1.25 < 1.30 * 6.78. 8 * 92.37. 109 *
* 1.30 < 1.35 * 3.39. 4 * 95.76. 113 *
* 1.35 < 1.40 * 1.69. 2 * 97.46. 115 *
* 1.40 < 1.45 * 0.85. 1 * 98.31. 116 *
* 1.45 < 1.50 * 0.00. 0 * 98.31. 116 *
* 1.50 < 1.55 * 0.00. 0 * 98.31. 116 *
* 1.55 < 1.60 * 1.69. 2 * 100.00. 118 *
* 1.60 < 1.65 * 0.00. 0 * 100.00. 118 *
* 1.65 < 1.70 * 0.00. 0 * 100.00. 118 *
* 1.70 < 1.75 * 0.00. 0 * 100.00. 118 *
* 1.75 < 1.80 * 0.00. 0 * 100.00. 118 *
* 1.80 < 1.85 * 0.00. 0 * 100.00. 118 *
* 1.85 < 1.90 * 0.00. 0 * 100.00. 118 *
* 1.90 < 1.95 * 0.00. 0 * 100.00. 118 *
* 1.95 < 2.00 * 0.00. 0 * 100.00. 118 *
* 2.00 < 2.05 * 0.00. 0 * 100.00. 118 *
* 2.05 < 2.10 * 0.00. 0 * 100.00. 118 *
* 2.10 < 2.15 * 0.00. 0 * 100.00. 118 *
* 2.15 < 2.20 * 0.00. 0 * 100.00. 118 *
* 2.20 < 2.25 * 0.00. 0 * 100.00. 118 *
* 2.25 < 2.30 * 0.00. 0 * 100.00. 118 *
* 2.30 < 2.35 * 0.00. 0 * 100.00. 118 *
* 2.35 < 2.40 * 0.00. 0 * 100.00. 118 *
* 2.40 < 2.45 * 0.00. 0 * 100.00. 118 *
* 2.45 < 2.50 * 0.00. 0 * 100.00. 118 *
* 2.50 < inf. * 0.00. 0 * 100.00. 118 *
**********************************************************
Code_Aster
®
Version
8.2
Titrate:
Macro-control
MACR_INFO_MAIL
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.03.02-D1
Page:
8/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
On the following figure, one can see the representation in the form of histogram of the percentages
elements in each range of quality concerned.
0
5
10
15
20
25
1,00/1,05
1,05/1,
1
1,15/1,
2
1,20/1,25
1,30/1,35
1,35/1,
4
1,45/1,
5
1,50/1,55
1,60/1,
6
Range of quality
Percentage of élém
ents
Information on the number of entities of calculation is classified by type: nodes, mesh-points,
edges, etc. One makes the sorting between the elements of edges and the others.
A NUMBER Of ENTITIES OF CALCULATION
==========================
**********************************************************
* Nodes *
**********************************************************
* Numbers total * 76 *
**********************************************************
**********************************************************
* Edges *
**********************************************************
* Numbers total * 32 *
*. of which edges isolees * 0 *
*. of which edges of edge of areas 2D * 32 *
*. of which edges intern with the faces/volumes * 0 *
**********************************************************
**********************************************************
* Triangles *
**********************************************************
* Numbers total * 118 *
**********************************************************
Code_Aster
®
Version
8.2
Titrate:
Macro-control
MACR_INFO_MAIL
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.03.02-D1
Page:
9/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
The control of the connexity makes it possible to make sure that the field of calculation is only one piece. It
can arrive that this field is in several disjoined parts, for example for problems of
contact. But the diagnosis of several pieces often reveals a problem during creation of
mesh: under-parts were badly restuck. And as the visual monitoring always does not allow
to realize it, calculation is false… while seeming to be right. The worst situation which can
to arrive. It is thus strongly recommended to use this option.
CONNEXITY OF THE ENTITIES OF CALCULATION
===============================
**********************************************************
* The triangles are in only one block. *
**********************************************************
The display of the sizes of the under-parts of the mesh is done by dimensions. Tri A place roughly speaking
according to the groups of meshs declared; there can be a finer subdivision when the groups
recut. This presentation makes it possible to control that one has quite with a grid what one wanted. In
private individual, that allows to flush out errors of sticking together or creation of groups at the time of
creation of the mesh. Here still, it is strongly recommended to use this option.
SIZES OF THE SOUS-DOMAINES OF CALCULATION
===================================
Direction | Unit
-----------------------------------
X | Unknown factor
y | Unknown factor
**********************************************************
* Under-fields 2D *
**********************************************************
* Number * Name * Surface *
**********************************************************
* - 12 * MANCHE__________________________ * 20.000 *
* - 11 * MASSE___________________________ * 128.00 *
**********************************************************
* Total: * 148.00 *
**********************************************************
**********************************************************
* 1D Under-fields *
**********************************************************
* Number * Name * Length *
**********************************************************
* - 2 * BORD_2__________________________ * 10.000 *
* - 1 * BORD_1__________________________ * 4.0000 *
* - 3 * BORD_3__________________________ * 16.000 *
* - 5 * BORD_5__________________________ * 8.0000 *
* - 4 * BORD_4__________________________ * 26.246 *
**********************************************************
* Total: * 64.246 *
**********************************************************
Code_Aster
®
Version
8.2
Titrate:
Macro-control
MACR_INFO_MAIL
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.03.02-D1
Page:
10/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
5 Bibliography
[1]
G. NICOLAS: “Software LOBSTER - Volume 1 - general Presentation”, report/ratio EDF
HI-23/04/005, February 2005.
[2]
G. NICOLAS: “Software LOBSTER - Volume 2 Algorithmes of refinement and déraffinement
mesh ", report/ratio EDF HI-23/04/006, February 2005.
[3]
G. NICOLAS: “Software LOBSTER - Volume 3 Interfaces with the computer codes”, report/ratio
EDF HI-23/04/007, February 2005.