Code_Aster
®
Version
8.2
Titrate:
Procedure
IMPR_FICO_HOMA
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.04.01-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.04.01
Procedure
IMPR_FICO_HOMA
1 Goal
To print the file of LOBSTER configuration.
This procedure makes it possible to write the whole of the data necessary to the LOBSTER execution:
·
type of use,
·
thresholds of refinement and déraffinement of mesh,
·
analyze mesh.
Note:
A priori, this procedure should never be employed directly by an end-user.
It is activated through the macro-controls of adaptation of mesh.
Code_Aster
®
Version
8.2
Titrate:
Procedure
IMPR_FICO_HOMA
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.04.01-D1
Page:
2/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
2 Syntax
IMPR_FICO_HOMA
(
PROCESSING =
_F (
# choice of a type of processing
/
ADAPTATION =/“RAFF_DERA”
,
/
“REFINEMENT”
,
/
“DERAFFINEMENT”
,
/
UNIFORM
=
/
“REFINEMENT”
,
/
“DERAFFINEMENT”
,
/
“NOTHING”
,
/
INFORMATION
=
“YES”
,
# Name MED of the mesh in input
NOM_MED_MAILLAGE_N = nmmaillagen,
[TXM]
# If ADAPTATION or UNIFORM
# name MED of the mesh at exit
NOM_MED_MAILLAGE_NP1 = nmmaillagenp1
,
[TXM]
# number of the current iteration
NITER = niter,
[I]
MAJ_CHAM =
/
“YES”,
/
“NOT”
,
[DEFECT]
#
Finsi
# If ADAPTATION =
'
REFINEMENT
'
or ADAPTATION =
'
RAFF_DERA
'
/CRIT_RAFF_PE
= crp
,
[R]
/
CRIT_RAFF_ABS
=
CRA
,
[R]
/
CRIT_RAFF_REL
=
crr
,
[R]
#
Finsi
# If ADAPTATION =
'
DERAFFINEMENT
'
or ADAPTATION =
'
RAFF_DERA
'
/CRIT_DERA_PE
= cdp
,
[R]
/
CRIT_DERA_ABS
=
half-value layer
,
[R]
/
CRIT_DERA_REL
=
cdr
,
[R]
#
Finsi
# If refinement, adaptive or uniform
NIVE_MAX
=
nivmax
,
[I]
#
Finsi
# So déraffinement, adaptive or uniform
NIVE_MIN
=
nivmin
,
[I]
#
Finsi
# If ADAPTATION, names of fields MED related to the indicator
NOM_MED =
nmindica
,
[TXM]
COMPONENT =
mcindica
,
[TXM]
RESULT
=
nrindica
,
[TXM]
NOM_CHAM
=
nhindica
,
[TXM]
Selection of the temporal parameter
/
NUME_ORDRE
=
noindica
[I]
/
INST
=
inst
[R]
I
PRECISION
=
/
prec [R]
/
1.0E-3
[DEFECT]
I CRITERION =/“RELATIVE”
[DEFECT]
/
“ABSOLUTE”
#
Finsi
)
#
Follow-up of a border
NOM_MED_MAILLAGE_FRONTIERE
=
maf
[TXM]
GROUP_MA
=
l_grma [l_gr_maille]
Code_Aster
®
Version
8.2
Titrate:
Procedure
IMPR_FICO_HOMA
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.04.01-D1
Page:
3/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
# Updated of fields on the new mesh
MAJ_CHAM = _F (
NOM_MED =
nmchamp
[TXM]
Selection of the temporal parameter
/
RESULT
=
resu [evol_elas]
[evol_noli]
[evol_ther]
[evol_thme]
NOM_CHAM
=
nomsymb [K16]
/
NUME_ORDRE
=
command
[I]
/
INST
=
moment
[R]
I
PRECISION
=
/
prec [R]
/
1.0E-3
[DEFECT]
I CRITERION =/“RELATIVE”
[DEFECT]
/
“ABSOLUTE”
)
# Management of the types of elements
NON_SIMPLEXE
=
/
0
[DEFECT]
/1
/2
ANALYZE =
_F (
/A NUMBER
=
/
“YES”
,
[DEFECT]
“NOT”,
/
QUALITY =/
“YES”
,
[DEFECT]
/
“NOT”
,
/
INTERPENETRATION
=
/
“YES”
,
[DEFECT]
/
“NOT”
,
/
CONNEXITY =/
“YES”
,
[DEFECT]
/
“NOT”,
/
CUT
=/
“YES”
,
[DEFECT]
“NOT”,
)
LANGUAGE
=
/
“French”
,
[DEFECT]
/
“FRENCH”
,
/
“ENGLISH”
,
/
“ENGLISH”
,
INFORMATION =
/
1
,
[DEFECT]
/
2
,
# Reference symbol of the files necessary to LOBSTER
FICHIER_CONF
=
f_conf
[TXM]
FICHIER_DONN
=
f_donn
[TXM]
)
Code_Aster
®
Version
8.2
Titrate:
Procedure
IMPR_FICO_HOMA
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.04.01-D1
Page:
4/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
3 Operands
3.1 Key word
PROCESSING
PROCESSING = _F (
This operand makes it possible to choose the mode of use of LOBSTER. The selection is done initially
among three headings: free adaptation of a mesh, uniform adaptation of a mesh,
information on a mesh.
3.1.1 Operand
ADAPTATION
/ADAPTATION =/“RAFF_DERA”,
/
“REFINEMENT”
,
/
“DERAFFINEMENT”
,
This operand is employed to make free adaptation of a mesh. In other words,
the decision of (die) refining an element is caught according to the value of an indicator
of error calculated before.
The choice can be done between three alternatives:
·
“RAFF_DERA”: the mesh is refined and déraffiné according to the indicator
of error,
·
“REFINEMENT”: only the function of refinement is activated. Elements with
low level of error are not déraffinés,
·
“DERAFFINEMENT”: it is the reverse, only the function of déraffinement is activated.
The elements on high level of error are not refined.
3.1.2 Operand
UNIFORM
/UNIFORM =
/“REFINEMENT”,
/
“DERAFFINEMENT”
,
/
“NOTHING”
,
This operand is employed to make a uniform adaptation of a mesh. In others
terms, all the elements of the mesh are treated same manner, without holding account
of an indicator of error. The choice can be done between three alternatives:
·
“REFINEMENT”: all the elements are refined,
·
“DERAFFINEMENT”: all the elements are déraffinés,
·
“NOTHING”: all the elements are preserved; the mesh is the same one at the exit as with
the input.
3.1.3 Operand
INFORMATION
/
INFORMATION
=
“YES”
In this case, LOBSTER does not adapt the mesh but delivers information on sound
contents.
3.1.4 Operand
NOM_MED_MAILLAGE_N
NOM_MED_MAILLAGE_N = nmmaillagen
It is the name according to conventions' MED under which the mesh of input in the file is known
transferred between Code_Aster and LOBSTER.
Code_Aster
®
Version
8.2
Titrate:
Procedure
IMPR_FICO_HOMA
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.04.01-D1
Page:
5/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
3.1.5 Operand
NOM_MED_MAILLAGE_NP1
NOM_MED_MAILLAGE_NP1 = nmmaillagenp1
It is the name according to conventions' MED under which will be recorded the mesh of exit in
file transferred between LOBSTER and Code_Aster. This data is obligatory when one is done
adaptation, free or uniform.
3.1.6 Operand
NITER
NITER = niter
Number of the iteration of the mesh of input in LOBSTER.
This number is worth 0 for a mesh not having undergone any adaptation. It is worth N for a mesh
obtained after N adaptations. This data is obligatory when an adaptation, free is made or
uniform.
3.1.7 Operand
MAJ_CHAM
MAJ_CHAM =
/“YES”,
/
“NOT”
,
So “YES”, the function of update of fields on the new mesh will be activated in
LOBSTER.
So “NOT”, nothing is done.
3.1.8 Operands of the indicator
It is necessary to transmit to LOBSTER the characteristics of the indicator of error. It is obligatory to specify
under which field name and which name of component the indicator is stored in file MED.
Then the alternative is as follows:
·
that is to say the sequence number was provided in the options of MACR_ADAP_MAIL; it is necessary then
to inform NUME_ORDRE,
·
that is to say one moment was provided in the options of MACR_ADAP_MAIL; it is then necessary to inform it
key word INST,
·
either nothing was provided; Aster will choose that which exists for the field specified by RESULT and
NOM_CHAM.
3.1.8.1 Operand
NOM_MED
NOM_MED = nmindica
Name according to conventions' MED of the field of the indicator of error such as it was written in
file transferred between Code_Aster and LOBSTER. It is a channel of 32 characters.
3.1.8.2 Operand
COMPONENT
COMPONENT = ncindica
Name of the component of the field of the indicator of error which will be actually used for
to adapt the mesh. It is a channel of 8 characters.
3.1.8.3 Operand
NUME_ORDRE
NUME_ORDRE
=
noindica
It is the sequence number of the field of indicator of error in its file MED.
3.1.8.4 Operand
INST
INST = inst
It is the moment when the indicator of error is retained.
Code_Aster
®
Version
8.2
Titrate:
Procedure
IMPR_FICO_HOMA
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.04.01-D1
Page:
6/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
3.1.8.5 Operand
PRECISION
PRECISION = prec
Allows to test the equality between the actual value of the moment requested and that presents in
file.
3.1.8.6 Operand
CRITERION
CRITERION = crit
Specify if the test of equality takes place in absolute or relative.
3.1.8.7 Operand
RESULT
RESULT
= nrindica
It is the name of the concept result associated with the indicator of error.
3.1.8.8 Operand
NOM_CHAM
NOM_CHAM
= nhindica
It is the field name symbolic system associated with the indicator with error.
3.1.9 Operand
CRIT_RAFF_xxxx
In the case of free adaptation implying of the refinement of mesh, it is necessary to define a high criterion
of error. All the elements whose error is higher than this criterion will be refined. Three alternatives are
possible:
3.1.9.1 Operand
CRIT_RAFF_PE
/CRIT_RAFF_PE = crp
The criterion is defined by a proportion of elements to refine. It is a real number included/understood
between 0 and 1. The process is as follows:
·
calculation of the number of elements N corresponding to the proportion defined by crp is
N = crp X numbers total elements,
·
N elements with the strongest error will be refined.
3.1.9.2 Operand
CRIT_RAFF_ABS
/CRIT_RAFF_ABS = CRA
The criterion is defined by an absolute value of the error. All elements with an error
higher than this value will be refined.
3.1.9.3 Operand
CRIT_RAFF_REL
/CRIT_RAFF_REL = crr
The criterion is defined by a relative value of the error. It is a number ranging between 0 and 1.
The process is as follows:
·
calculation of the minimal and maximum values of the indicator of error,
·
calculation of the value corresponding to the proportion of error: v = vmin + crr (vmax
vmin),
·
all the elements whose error is higher than this value will be refined.
3.1.10 CRIT_DERA_xxxx operand
In the case of free adaptation implying of déraffinement, it is necessary to define a low criterion of error. All
the elements whose error is lower than this criterion will be déraffinés. Three alternatives are possible.
Code_Aster
®
Version
8.2
Titrate:
Procedure
IMPR_FICO_HOMA
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.04.01-D1
Page:
7/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
3.1.10.1 Operand CRIT_DERA_PE
/CRIT_DERA_PE = cdp
The criterion is defined by a proportion of elements in déraffiner. It is a number included/understood enters
0 and 1. The process is as follows:
·
calculation of the number of elements N corresponding to the proportion defined by cdp is N =
cdp X numbers total elements
·
N elements with the weakest error will be déraffinés.
3.1.10.2 Operand CRIT_DERA_ABS
/CRIT_DERA_ABS = half-value layer
The criterion is defined by an absolute value of the error. All elements with an error
lower than this value will be déraffinés.
3.1.10.3 Operand CRIT_DERA_REL
/CRIT_DERA_REL = cdr
The criterion is defined by a relative value of the error. It is a number ranging between 0 and 1.
The process is as follows:
·
calculation of the minimal and maximum values of the indicator of error,
·
calculation of the error value V corresponding to the proportion of cdr error such as: v
= vmin + cdr (vmax vmin),
·
all the elements whose error is lower than this value will be déraffinés.
3.1.11 Operand NIVE_MAX
NIVE_MAX = nivmax
It is the maximum level of refinement of the mesh. In other words an element of the initial mesh
could not be divided more nivmax time as a whole of the process.
3.1.12 Operand NIVE_MIN
NIVE_MIN = nivmin
It is the minimal level of déraffinement mesh. I.e. only elements resulting
of at least nivmin cuttings of mesh can be déraffinés.
3.1.13 Operand NOM_MED_MAILLAGE_FRONTIERE
NOM_MED_MAILLAGE_FRONTIERE = nmmaf
It is the name under which the mesh of the border in file MED is known.
3.1.14 Operand GROUP_MA
GROUP_MA = l_grma
Give the list of the groups of meshs defining the border.
3.2 Word
key
MAJ_CHAM
MAJ_CHAM = _F
(
This key word is to be employed once as many as one has fields to update between the old one
mesh and adapted mesh.
Code_Aster
®
Version
8.2
Titrate:
Procedure
IMPR_FICO_HOMA
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.04.01-D1
Page:
8/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
3.2.1 Operand
NOM_MED
NOM_MED
=
nmchamp
Name MED of the field to be updated.
3.2.2 Selection of the temporal parameter
The selection of the sequence number associated with the field to be interpolated is done by the designation of a number
of command or a value of moment or by identifying the field in the structure of result.
3.2.2.1 Operand
RESULT
/
RESULT
= resu
Name of the concept [result] containing the field to be updated.
3.2.2.2 Operand
NOM_CHAM
NOM_CHAM
= nomch
Reference symbol of the field to be updated.
3.2.2.3 Operand
NUME_ORDRE
/
NUME_ORDRE
=
nochamp
It is the sequence number of the field to be updated in its file MED.
3.2.2.4 Operand
INST
/INST = inst
It is the moment when the field is retained.
3.3 Word
key
NON_SIMPLEXE
NON_SIMPLEXE
=
/
0
[DEFECT]
1
2
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.4 Word
key
ANALYZE
ANALYZE = _F (
Allows to specify the type of analysis wished for the mesh.
3.4.1 Operand
QUALITY
/
QUALITY
=
/
“YES”
,
[DEFECT]
/
“NOT”
,
So “YES”, a summary of the quality of the elements of the mesh will be published.
So “NOT”, nothing is done.
Code_Aster
®
Version
8.2
Titrate:
Procedure
IMPR_FICO_HOMA
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.04.01-D1
Page:
9/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
3.4.2 Operand
INTERPENETRATION
/
INTERPENETRATION
=
/
“YES”
,
[DEFECT]
/
“NOT”
,
So “YES”, LOBSTER tests if the elements of the mesh do not overlap.
So “NOT”, nothing is done.
3.4.3 Operand
NUMBERS
/
NUMBERS
=
/
“YES”
,
[DEFECT]
“NOT”
,
So “YES”, a summary of the numbers of nodes and elements are printed.
So “NOT”, nothing is done.
3.4.4 Operand
CONNEXITY
/
CONNEXITY
=
/
“YES”
,
[DEFECT]
“NOT”
,
So “YES”, information on the connexity of the elements of the mesh is printed.
So “NOT”, nothing is done.
3.4.5 Operand
CUT
/
CUT
=
/
“YES”
,
[DEFECT]
“NOT”
,
So “YES”, the sizes of the various under-fields are printed.
So “NOT”, nothing is done.
3.5 Operand
LANGUAGE
LANGUAGE
=
“French”
,
[DEFECT]
“FRENCH”
,
“ENGLISH”
,
“ENGLISH”
,
The operand defines the language used for the messages resulting from LOBSTER.
3.6 Operand
FICHIER_CONF
FICHIER_CONF = f_conf
Reference symbol of the file of configuration for LOBSTER which will be written by this procedure.
3.7 Operand
FICHIER_DONN
FICHIER_DONN = f_donn
Reference symbol of the data file for LOBSTER which will be written by this procedure.
3.8 Operand
INFORMATION
INFORMATION =
/1,
/
2
,
If INFORMATION is worth 2, some impressions of déboggage take place.
If not, nothing takes place.
Code_Aster
®
Version
8.2
Titrate:
Procedure
IMPR_FICO_HOMA
Date:
31/01/06
Author (S):
G. NICOLAS
Key
:
U7.04.01-D1
Page:
10/10
Instruction manual
U7.0- booklet: Data exchanges
HT-62/06/004/A
4 Example
IMPR_FICO_HOMA (PROCESSING = _F (NOM_MED_MAILLAGE_N=' MZERO', INFORMATION=' OUI'),
ANALYZE = _F (QUALITY = “YES”,
INTERPENETRATION
=
“YES”),
ITEM
=
“MZERO_INFO”)
This control will involve the writing of a file of configuration requiring of LOBSTER to analyze
a mesh. This file will reside in sub-directory MZERO_INFO of the current directory.