Code_Aster ®
Version
8.2
Titrate:
Procedure IMPR_MATRICE
Date:
31/01/06
Author (S):
J. Key PELLET
:
U7.04.32-F1 Page:
1/10
Organization (S): EDF-R & D/AMA
Handbook of Utilization
U7.0- booklet: Data exchanges
Document: U7.04.32
Procedure IMPR_MATRICE
1 Goal
To print on a file with format IDEAS or RESULTAT of the elementary matrices and assemblies.
The impressions with format RESULTAT are especially useful for the developers of the code at ends of
“débuggage”.
Handbook of Utilization
U7.0- booklet: Data exchanges
HT-62/06/004/A
Code_Aster ®
Version
8.2
Titrate:
Procedure IMPR_MATRICE
Date:
31/01/06
Author (S):
J. Key PELLET
:
U7.04.32-F1 Page:
2/10
2 Syntax
IMPR_MATRICE
(
|
MATR_ELEM =_F (
#
impressions in the form of a universal file IDEAS.
/
FORMAT
=
“IDEAS”
,
[DEFAUT]
MATRICE = mel
,
[matr_elem_DEPL_R]
UNITE
=
/
30,
[DEFAUT]
/
iunit,
[I]
VERSION =
5,
[DEFAUT]
# impressions in a listing (débuggage):
/
FORMAT
=
“RESULTAT”,
MATRICE = mel
,
/
[matr_elem_ *]
/
[vect_elem_ *]
UNITE
=
/
8, [DEFAUT]
/
iunit,
[I]
#
choice
meshs
:
/ALL = “YES”,
/ |
NOEUD =
l_no,
[l_noeud]
| GROUP_NO
= l_gno,
[l_group_no]
|
MAILLE =
l_ma,
[l_maille]
| GROUP_MA
= l_gma,
[l_group_ma]
NOM_CMP
=
l_cmp,
[l_Kn]
GRAIN
=/
“VALEUR”,
[DEFAUT]
/
“NOEUD”
,
/
“MAILLE”,
NB_CHIFFRE =/4,
[DEFAUT]
/
nchif,
[I]
/
999
,
#
precision
max.
)
Handbook of Utilization
U7.0- booklet: Data exchanges
HT-62/06/004/A
Code_Aster ®
Version
8.2
Titrate:
Procedure IMPR_MATRICE
Date:
31/01/06
Author (S):
J. Key PELLET
:
U7.04.32-F1 Page:
3/10
|
MATR_ASSE =_F (
# impressions in the form of a universal file IDEAS.
/
FORMAT
=
“IDEAS”
,
[DEFAUT]
MATRICE = chechmate
,
[matr_asse_DEPL_R]
UNITE
=
/
30,
[DEFAUT]
/
iunit,
[I]
VERSION =
5,
[DEFAUT]
# impressions in a listing (débuggage):
/
FORMAT
=
“RESULTAT”
,
MATRICE = mel
,
[matr_asse_ *]
UNITE
=
/
8, [DEFAUT]
/
iunit,
[I]
#
choice
nodes
:
/ALL = “YES”,
/ |
NOEUD =
l_no,
[l_noeud]
| GROUP_NO
= l_gno,
[l_group_no]
|
MAILLE =
l_ma,
[l_maille]
| GROUP_MA
= l_gma,
[l_group_ma]
OPTION=/
“SOUS_MATRICE”
, [DEFAUT]
/
“LIGNE”
,
/
“COLONNE”
,
NOM_CMP=
l_cmp,
[l_Kn]
GRAIN =
/“VALEUR”,
[DEFAUT]
/
“NOEUD”
,
NB_CHIFFRE =/4,
[DEFAUT]
/
nchif,
[I]
/
999
,
#
precision
max.
VALE_ZERO
=
/
0.
,
[DEFAUT]
/
eps,
[R]
)
)
Handbook of Utilization
U7.0- booklet: Data exchanges
HT-62/06/004/A
Code_Aster ®
Version
8.2
Titrate:
Procedure IMPR_MATRICE
Date:
31/01/06
Author (S):
J. Key PELLET
:
U7.04.32-F1 Page:
4/10
3 Operands
3.1
Impressions in the form of a universal file IDEAS
3.1.1 Key word factor MATR_ELEM
Allows to specify to it (S) concept (S) of matr_elem type to be printed as well as the format of impression.
3.1.1.1 Operand
FICHIER
FICHIER =
nom_fich
[KN]
Name of the file in which one wishes to see printing the data to be transferred. This name of
file must be defined beforehand by procedures DEBUT [U4.11.01], POURSUITE
[U4.11.03] or DEFI_FICHIER [U4.12.03].
3.1.1.2 Operand
FORMAT
Format under which one wishes to print the data.
3.1.1.3 Operand
VERSION
VERSION =
5
[DEFAUT]
Number of the version of IDEAS corresponding to the format of the file used to print them
data. Only version 5 (default value) is supported for the moment.
3.1.1.4 Operand
MATRICE
MATRICE = mel
Name of the elementary matrices (standard matr_elem_DEPL_R) that one wishes to print on one
file.
3.1.2 Key word factor MATR_ASSE
Allows to specify to it (S) concept (S) of type matr_asse to print as well as the format of impression.
3.1.2.1 Operand
FICHIER
FICHIER =
nom_fich
[KN]
Name of the file in which one wishes to see printing the data to be transferred. This name of
file must be defined beforehand by procedures DEBUT, POURSUITE or
DEFI_FICHIER.
3.1.2.2 Operand
FORMAT
FORMAT
=
“IDEAS” [DEFECT]
Format under which one wishes to print the data. Only the format by defect “IDEAS” is
supported for the moment.
3.1.2.3 Operand
VERSION
VERSION =
5
[DEFAUT]
Number of the version of IDEAS corresponding to the format of the file used to print them
data. Only version 5 (default value) is supported for the moment.
Handbook of Utilization
U7.0- booklet: Data exchanges
HT-62/06/004/A
Code_Aster ®
Version
8.2
Titrate:
Procedure IMPR_MATRICE
Date:
31/01/06
Author (S):
J. Key PELLET
:
U7.04.32-F1 Page:
5/10
3.1.2.4 Operand
MATRICE
MATRICE = mel
Name of the assembled matrix (standard matr_asse_DEPL_R) that one wishes to print on one
file.
3.2
Impressions on listing
3.2.1 Operand
FORMAT = “RESULT”
It is the choice of the format RESULTAT which makes it possible to print in “light” (on listing) the contents of
matrices (or elementary vectors).
3.2.2 What can one print?
A priori all concepts of the type MATR_ASSE_ * except the MATR_ASSE_GENE or MATR_ELEM_ * or
VECT_ELEM_ *
For the MATR_ASSE_ *, one will use the key word factor MATR_ASSE,
for the MATR_ELEM_ * (or VECT_ELEM_ *), one will use the key word factor MATR_ELEM.
3.2.3 Operand
FICHIER
FICHIER =
nom_fich
[KN]
Name of the file in which one wishes to see printing the treated data. By defect, it acts
file “RESULTAT”.
3.2.4 Operand
MATRICE
MATRICE = chechmate
Name of the assembled or elementary matrix which one wishes to print on a file.
Note:
One cannot print the matr_asse_gene but they are simple objects
(full matrices) for which command IMPR_CO is sufficient.
3.2.5 To select the meshs or the nodes which one wants to print
3.2.5.1 Case
RESUELEM (matrix or vector)
Choice of the elements which one wants to print:
has)
/list of meshs: lma
only meshs of lma
b)/
TOUT=' OUI'
all meshs
c)
/list of nodes: lno
all meshs resting on at least
a node of the list of nodes lno
Note:
The possibility c) makes it possible to print the matrices (or the vectors) elementary of “Lagrange” which
force a physical node which one knows the name.
Handbook of Utilization
U7.0- booklet: Data exchanges
HT-62/06/004/A
Code_Aster ®
Version
8.2
Titrate:
Procedure IMPR_MATRICE
Date:
31/01/06
Author (S):
J. Key PELLET
:
U7.04.32-F1 Page:
6/10
3.2.5.2 Case
MATR_ASSE
Choice of the nodes which one wants to print:
has)
/list of meshs: lma
nodes of the meshs of lma
b)/
TOUT=' OUI'
all “physical " nodes carrying ddls
c)
/list of nodes: lno
only nodes of lno.
3.2.5.3 Word
key
OPTION
By defect (OPTION = “MATRICE”), one prints only the terms of the MATR_ASSE defined by
“under” matrix corresponding to the choice of the nodes which one has just made. If the user asks
(OPTION = “LIGNE”), one will print all the contents of the lines corresponding to the nodes chosen,
(in the same way for OPTION = “COLONNE”).
3.2.6 To select
CMPS to be printed
NOM_CMP = l_kn list of the components which one wants to print (by defect: all).
3.2.7 Format
of impression
3.2.7.1 General
The matrices, even if they are symmetrical, will be entirely printed: they are “rectangles”
full.
NB_CHIFFRE = nchif
nchif is the number of decimals to be written for each floating number.
VALE_ZERO = eps
One will print only the values whose absolute value is higher than eps. This key word is not
valid that for the MATR_ASSE, to limit the impressions.
3.2.7.2 Case
RESUELEM (matr_elem or vect_elem)
GRAIN = “VALUE”
one prints a value by line of listing.
Example: RESUELEM (matrix) = mel MAILLE = (“M1”, “m2”), NOM_CMP= “TEMP”
M1/NO1/“TEMP”
M1/NO1/“TEMP”
12.
M1/NO1/“TEMP”
M1/NO4/“TEMP”
12.
M1/NO1/“TEMP”
M1/NO7/“TEMP”
12.
M1/NO4/“TEMP”
M1/NO4/“TEMP”
12.
M1/NO4/“TEMP”
M1/NO7/“TEMP”
12.
…
M2/NO4/“TEMP”
M2/NO4/“TEMP”
7.
…
Handbook of Utilization
U7.0- booklet: Data exchanges
HT-62/06/004/A
Code_Aster ®
Version
8.2
Titrate:
Procedure IMPR_MATRICE
Date:
31/01/06
Author (S):
J. Key PELLET
:
U7.04.32-F1 Page:
7/10
Example: RESUELEM (vector) = mel, MAILLE = “M1”, NOM_CMP= “TEMP”
M1/NO1/“TEMP”
12.
M1/NO4/“TEMP”
12.
M1/NO7/“TEMP”
12.
GRAIN = “NODE”
In this case, the “grain” of impression is the node.
Example: RESUELEM (matrix) = mel, MAILLE = “M1”, NOM_CMP= (“DX”, “DZ”)
M1/NO1 X NO1
DX
DZ
DX 12.
13.
DX 13.
15.
M1/NO1 X NO4
DX
DZ
DX 12.
13.
DX 13.
15.
…
GRAIN = “MESH”
In this case, the “grain” of impression is the mesh.
Example: RESUELEM = mel, MAILLE = “M1”, NOM_CMP= (“DX”, “DZ”)
M1
N1/DX N1/DZ N3/DX N3/DZ…
N1/DX 12. 13. 14. 15. …
N1/DZ 13. 13. 14. 15. …
N3/DX 14. 14. 14. 15. …
…
3.2.7.3 MATR_ASSE
For a MATR_ASSE, one will print only the “grains” for which at least a value is nonnull
(see key word VALE_ZERO in paragraph 3.2.7.1).
GRAIN = “VALUE”
one prints a value by line of listing.
Example: MATR_ASSE=matass, NOEUD= (“NO1”, “NO3”, “NO7”), NOM_CMP= (“DX”, “DZ”),
OPTION
=
“SOUS_MATRICE”
NO1/“DX”
NO1/“DX”
12.
NO1/“DX”
NO1/“DZ”
12.
NO1/“DX”
NO3/“DX”
12.
NO1/“DX”
NO3/“DZ”
12.
NO1/“DX”
NO7/“DX”
12.
NO1/“DX”
NO7/“DZ”
12.
NO3/“DX”
NO3/“DX”
12.
NO3/“DX”
NO3/“DZ”
12.
…
NO7/“DZ”
NO7/“DZ”
12.
Handbook of Utilization
U7.0- booklet: Data exchanges
HT-62/06/004/A
Code_Aster ®
Version
8.2
Titrate:
Procedure IMPR_MATRICE
Date:
31/01/06
Author (S):
J. Key PELLET
:
U7.04.32-F1 Page:
8/10
Example: MATR_ASSE = matass, NOEUD = (“NO1”), NOM_CMP= (“DX”, “DZ”),
OPTION
=
“COLONNE”
NO1/“DX”
NO1/“DX”
12.
NO1/“DX”
NO1/“DZ”
12.
NO1/“DZ”
NO1/“DX”
12.
NO1/“DZ”
NO1/“DZ”
12.
NO2/“DX”
NO1/“DX”
12.
NO2/“DX”
NO1/“DZ”
12.
NO2/“DZ”
NO1/“DX”
12.
NO2/“DZ”
NO1/“DZ”
12.
NO3/“DX”
NO1/“DX”
12.
NO3/“DX”
NO1/“DZ”
12.
…
NO999
/“DZ”
NO1/“DX”
12.
NO999
/“DZ”
NO1/“DZ”
12.
Example: MATR_ASSE = matass, NOEUD = (“NO1”), NOM_CMP= (“DX”, “DZ”),
OPTION
=
“LIGNE”
NO1/“DX” NO1
/“DX”
12.
NO1/“DZ” NO1
/“DX”
12.
NO1/“DX” NO1
/“DZ”
12.
…
NO1/“DX”
NO999/“DZ”
12.
NO1/“DZ”
NO999/“DZ”
12.
GRAIN = “NODE”
In this case, the “grain” of impression is the node.
Example: MATR_ASSE= matass, NOEUD = (“NO1”, “NO3”), NOM_CMP= (“DX”, “DZ”)
NO1 X NO1
DX
DZ
DX 12.
13.
DX 13.
15.
NO1 X NO3
DX
DZ
DX 12.
13.
DX 13.
15.
NO3 X NO3
DX
DZ
DX 12.
13.
DX 13.
15.
…
Handbook of Utilization
U7.0- booklet: Data exchanges
HT-62/06/004/A
Code_Aster ®
Version
8.2
Titrate:
Procedure IMPR_MATRICE
Date:
31/01/06
Author (S):
J. Key PELLET
:
U7.04.32-F1 Page:
9/10
4 Examples
One wants:
1) elementary matrices corresponding to the meshs S1 and S2,
2) 1 table by mesh,
3) 2 significant digits for the floating numbers,
4) one is interested only in the components DY and DRZ.
matel = CALC_MATR_ELEM (…)
IMPR_MATRICE (
MATR_ELEM = _F (
MATRICE
=
matel,
FORMAT = “RESULT”,
FILE = “RESULT”,
MAILLE = (“S1”, “S2”),
GRAIN = “MESH”,
NOM_CMP = (“DY”, “DRX”),
NB_CHIFFRE = 2))
One obtains in file RESULTAT:
_____________________________________________________
! LIGREL = MO .MODELE OPTION = RIGI_MECA_GE!
_____________________________________________________
_______________________________________________________________
!S1! N1/DY! N1/DRX! N2/DY! N2/DRX!
_______________________________________________________________
!N1/DY! - 5.3E-03! 5.8E-02! 5.3E-03! 5.8E-02!
_______________________________________________________________
!N1/DRX! 5.8E-02! - 2.1E+01! - 5.8E-02! 1.1E+01!
_______________________________________________________________
!N2/DY! 5.3E-03! - 5.8E-02! - 5.3E-03! - 5.8E-02!
_______________________________________________________________
!N2/DRX! 5.8E-02! 1.1E+01! - 5.8E-02! - 2.1E+01!
_______________________________________________________________
_______________________________________________________________
!S2! N2/DY! N2/DRX! N3/DY! N3/DRX!
_______________________________________________________________
!N2/DY! - 5.3E-03! 5.8E-02! 5.3E-03! 5.8E-02!
_______________________________________________________________
!N2/DRX! 5.8E-02! - 2.1E+01! - 5.8E-02! 1.1E+01!
_______________________________________________________________
!N3/DY! 5.3E-03! - 5.8E-02! - 5.3E-03! - 5.8E-02!
_______________________________________________________________
!N3/DRX! 5.8E-02! 1.1E+01! - 5.8E-02! - 2.1E+01!
_______________________________________________________________
Handbook of Utilization
U7.0- booklet: Data exchanges
HT-62/06/004/A
Code_Aster ®
Version
8.2
Titrate:
Procedure IMPR_MATRICE
Date:
31/01/06
Author (S):
J. Key PELLET
:
U7.04.32-F1 Page:
10/10
Intentionally white left page.
Handbook of Utilization
U7.0- booklet: Data exchanges
HT-62/06/004/A
Outline document