Code_Aster ®
Version
7.4
Titre :
Opérateur COMB_CHAM_ELEM
Date :
04/01/05
Auteur(s) :
X. DESROCHES, L. VIVAN Clé
:
U4.72.03-G Page
: 1/6
Organisme(s) : EDF-R&D/AMA, CS SI
Manuel d'Utilisation
Fascicule U4.7- : Opérations sur les résultats et les champs
Document : U4.72.03
Opérateur COMB_CHAM_ELEM
1 But
Effectuer la combinaison linéaire de champs aux éléments.
Tous les concepts de type cham_elem_* à combiner doivent être de même structure (mêmes listes
de groupes d'éléments, provenir de la même option de calcul).
Cet opérateur sert également à effectuer la recombinaison de modes de Fourier dans une direction
donnée (pour les déformations et contraintes en mécanique, les flux en thermique).
Le concept produit (cham_elem_*) peut être réentrant.
Manuel d'Utilisation
Fascicule U4.7- : Opérations sur les résultats et les champs
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur COMB_CHAM_ELEM
Date :
04/01/05
Auteur(s) :
X. DESROCHES, L. VIVAN Clé
:
U4.72.03-G Page
: 2/6
2 Syntaxe
chml
[cham_elem_*] = COMB_CHAM_ELEM
(
reuse = chml,
/ COMB_R=_F( CHAM_ELEM= ch , / [cham_elem_SIEF_R]
/
[cham_elem_FLUX_R]
/
[cham_elem_EPSI_R]
/
[cham_elem_ENER_R]
/
[cham_elem_CRIT_R]
/
[cham_elem_DBEL_R]
/
[cham_elem_PRES_R]
/
[cham_elem_SIEF_C]
COEF_R=
r
,
[R]
PARTIE= / 'REEL'
, [K]
/
'IMAG'
,
),
/
COMB_C
=_F(
CHAM_ELEM
=
ch,
/
[cham_elem_SIEF_R]
/
[cham_elem_SIEF_C]
/
COEF_R=
r
,
[R]
/
COEF_C=
c
,
[C]
),
/
COMB_FOURIER
=_F( CHAM_ELEM = ch,
/ [cham_elem_SIEF_R]
/
[cham_elem_FLUX_R]
/
[cham_elem_EPSI_R]
NUME_MODE=
nume,
[I]
TYPE_MODE= / 'SYME' , [K]
/ 'ANTI' ,
COEF_R= /
r
,
[R]
/ 1.,
[DEFAUT]
),
ANGL
=
angl
[R]
)
si COMB_R et CHAM_ELEM
:
[cham_elem_SIEF_R]
alors [*]
-> SIEF_R
[cham_elem_FLUX_R]
[*]
-> FLUX_R
[cham_elem_EPSI_R]
[*]
-> EPSI_R
[cham_elem_ENER_R]
[*]
-> ENER_R
[cham_elem_CRIT_R]
[*]
-> CRIT_R
[cham_elem_DBEL_R]
[*]
-> DBEL_R
[cham_elem_PRES_R]
[*]
-> PRES_R
[cham_elem_SIEF_C]
[*]
-> SIEF_R
si COMB_C et CHAM_ELEM
:
[cham_elem_SIEF_R]
alors [*]
-> SIEF_C
[cham_elem_SIEF_C]
[*]
-> SIEF_C
si COMB_FOURIER et CHAM_ELEM :
[cham_elem_SIEF_R] alors [*]
-> SIEF_R
[cham_elem_FLUX_R]
[*]
-> FLUX_R
[cham_elem_EPSI_R]
[*]
-> EPSI_R
Manuel d'Utilisation
Fascicule U4.7- : Opérations sur les résultats et les champs
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur COMB_CHAM_ELEM
Date :
04/01/05
Auteur(s) :
X. DESROCHES, L. VIVAN Clé
:
U4.72.03-G Page
: 3/6
3 Opérandes
3.1
Rappel sur la syntaxe des valeurs complexes
Les valeurs complexes peuvent être déclarées de deux façons différentes :
·
sous la forme
a + ib avec la syntaxe RI a, b
où
a et b sont des nombres réels.
·
sous la forme
(module, phase) avec MP module, phase
où
module et phase sont des nombres réels (phase en degré).
3.2 Mot
clé
COMB_R
/ COMB_R = _F
Description des termes de la combinaison linéaire produisant un cham_elem_* à coefficients
réels.
CHAM_ELEM = ch,
Nom du concept de type cham_elem_* à combiner.
COEF_R = r,
Coefficient réel à appliquer au concept ch.
PARTIE
=
Dans le cas de CHAM_ELEM complexe, on prend la partie réelle ou imaginaire.
3.3 Mot
clé
COMB_C
/ COMB_C = _F
Description des termes de la combinaison linéaire produisant un cham_elem_* à valeurs
complexes.
CHAM_ELEM = ch,
Nom du concept de type cham_elem_* à combiner.
/ COEF_R = r,
Coefficient réel à appliquer au concept ch.
/
COEF_C
=
c,
Coefficient complexe à appliquer au concept argument de CHAM_ELEM.
Manuel d'Utilisation
Fascicule U4.7- : Opérations sur les résultats et les champs
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur COMB_CHAM_ELEM
Date :
04/01/05
Auteur(s) :
X. DESROCHES, L. VIVAN Clé
:
U4.72.03-G Page
: 4/6
3.4 Mot
clé
COMB_FOURIER
/
COMB_FOURIER = _F
Recombinaison de FOURIER avec description harmonique par harmonique des termes de la
série de FOURIER.
En élasticité, la recombinaison de FOURIER sur les déformations (resp. les contraintes )
s'écrit :
l I4
04,2
sinl I4
04,2
() =
cos
s +
a
l
0
sin l I
0
cosl I
l
l L
-
2,4
2
2,4
2
où L est l'ensemble des modes (symétriques et/ou antisymétriques) qui ont été calculés.
In est la matrice identité d'ordre n,
sl : l est le numéro de l'harmonique (s pour symétrique),
al : l est le numéro de l'harmonique (a pour antisymétrique),
: angle de la section où sont calculées les déformations (resp. contraintes).
Dans la formule ci-dessus, tous les coefficients r valent 1. (cas standard).
En thermique, la recombinaison de FOURIER sur les flux s'écrit :
l I2
02,1
sinl I2
02,1
() =
cos
s +
a
l
0
sin l
0
cosl
l
l L
-
,
1 2
,
1 2
avec les mêmes notations que ci-dessus.
3.4.1 Opérande
CHAM_ELEM
CHAM_ELEM = ch,
Nom du concept de type cham_elem_* à recombiner.
3.4.2 Opérande
COEF_R
COEF_R = r,
Coefficient réel à appliquer à l'harmonique considérée (par défaut 1.)
3.4.3 Opérande
NUME_MODE
NUME_MODE = nume,
Numéro de l'harmonique.
3.4.4 Opérande
TYPE_MODE
TYPE_MODE = / 'SYME' ,
/
'ANTI' ,
Harmonique SYMEtrique ou ANTIsymétrique.
3.5 Opérande
ANGL
ANGL = angl,
Angle en degré où a lieu la recombinaison de FOURIER.
Manuel d'Utilisation
Fascicule U4.7- : Opérations sur les résultats et les champs
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur COMB_CHAM_ELEM
Date :
04/01/05
Auteur(s) :
X. DESROCHES, L. VIVAN Clé
:
U4.72.03-G Page
: 5/6
4
Phases de vérification / exécution
Phase de vérification
On vérifie que le "concept produit" n'est pas un opérande d'entrée de l'opérateur.
Phase d'exécution
·
A l'exécution l'opérateur vérifie avant d'effectuer la combinaison linéaire, que le type des
valeurs du résultat est compatible avec le type des valeurs des concepts cham_elem_* à
combiner.
·
On effectue également les vérifications de cohérence : les arguments de type cham_elem_*
doivent avoir les mêmes ligrels, c'est-à-dire avoir été calculés sur les mêmes listes de
groupes d'éléments.
5 Exemples
5.1
Combinaison linéaire classique
Utot = COMB_CHAM_ELEM ( COMB_C = (_F(CHAM_ELEM= U1, COEF_R =1.),
_F(CHAM_ELEM
=U2,
COEF_C
=RI
0.,
1.),),
)
Le concept produit Utot est de type cham_elem_*_C (complexe) :
Utot = U1 + i U2
Le concept Utot peut être utilisé à nouveau comme concept produit pour une nouvelle combinaison
avec U3 cham_elem_*_R (réels) :
Utot = COMB_CHAM_ELEM ( reuse = Utot,
COMB_C = (_F
(CHAM_ELEM = Utot, COEF_R = 1.),
_F (CHAM_ELEM
=
U3,
COEF_R
=
2.),),
)
Ce qui produit :
Utot = Utot + 2. U3
5.2
Recopie de concept de type cham_elem
CH_sauv = COMB_CHAM_ELEM (COMB_R = _F (CHAM_ELEM = ch, COEF_R = 1.)
)
5.3
Recombinaison de FOURIER en élasticité
sig45 = COMB_CHAM_ELEM
( COMB_FOURIER=
(_F ( CHAM_ELEM = sig1,
NUME_MODE
=
1,
TYPE_MODE
=
'SYME'
,),
_F ( CHAM_ELEM = sig2,
NUME_MODE
=
0,
TYPE_MODE
=
'ANTI'
,),
ANGLE
=
45.,
)
On calcule la recombinaison de FOURIER d'un champ de contraintes pour la section = 45. d'une
harmonique résultant d'un chargement mode 1 symétrique et d'une harmonique résultant d'un
chargement en torsion pure (mode 0 antisymétrique).
Manuel d'Utilisation
Fascicule U4.7- : Opérations sur les résultats et les champs
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur COMB_CHAM_ELEM
Date :
04/01/05
Auteur(s) :
X. DESROCHES, L. VIVAN Clé
:
U4.72.03-G Page
: 6/6
5.4
Recombinaison de FOURIER en thermique
flux0 = COMB_CHAM_ELEM
( COMB_FOURIER=
(_F ( CHAM_ELEM = flux1,
NUME_MODE
=
1,
TYPE_MODE
=
'SYME'
,),
_F ( CHAM_ELEM = flux2,
NUME_MODE
=
2,
TYPE_MODE
=
'SYME'
,),
ANGLE
=
0.,
)
On calcule la recombinaison de FOURIER pour la section = 0. de 2 modes symétriques : mode 1 et
mode 2.
Chacun des champs de flux (flux1 et flux2) provient d'une résolution sur le mode considéré.
Manuel d'Utilisation
Fascicule U4.7- : Opérations sur les résultats et les champs
HT-66/05/004/A
Document Outline