Code_Aster ®
Version
7.4
Titre :
Opérateur POST_DYNA_MODA_T
Date :
27/01/05
Auteur(s) :
E. BOYERE, L. VIVAN Clé
:
U4.84.02-F Page
: 1/8
Organisme(s) : EDF-R&D/AMA, CS SI
Manuel d'Utilisation
Fascicule U4.8- : Post traitement et analyses dédiées
Document : U4.84.02
Opérateur POST_DYNA_MODA_T
1 But
Post-traiter les résultats en coordonnées généralisées produits par DYNA_TRAN_MODAL. Deux options
sont disponibles : post traitement des non linéarités de choc ou des relations effort-déplacement. Dans
le premier cas, on choisit un diagnostic d'usure ou une meilleure connaissance des chocs survenus
pendant l'analyse transitoire. Le contenu de la table produite est imprimable sur le fichier RESULTAT
par la commande IMPR_TABLE [U4.91.03].
Manuel d'Utilisation
Fascicule U4.8- : Post traitement et analyses dédiées
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur POST_DYNA_MODA_T
Date :
27/01/05
Auteur(s) :
E. BOYERE, L. VIVAN Clé
:
U4.84.02-F Page
: 2/8
2 Syntaxe
table [table] = POST_DYNA_MODA_T (
RESU_GENE
=
tg,
[tran_gene]
/ CHOC
=_F(
INST_INIT =/ -1.0, [DEFAUT]
/
t0,
[R]
INST_FIN = / 999., [DEFAUT]
/
t1,
[R]
NB_BLOC = / 1,
[DEFAUT]
/
nb,
[I]
SEUIL_FORCE = / 0.0 [DEFAUT]
/ s, [R]
DUREE_REPOS = / 0.0,[DEFAUT]
/ d, [R]
OPTION = /'USURE', [DEFAUT]
/'IMPACT',
NB_CLASSE
= / 10,
[DEFAUT]
/
nc,
[I]
),
/ RELA_EFFO_DEPL
=_F(
NOEUD = noeu,
[noeud]
NOM_CMP = noncmp, [K8]
),
INFO
=
/
1,
[DEFAUT]
/ 2,
TITRE
= titre,
[l_Kn]
)
Manuel d'Utilisation
Fascicule U4.8- : Post traitement et analyses dédiées
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur POST_DYNA_MODA_T
Date :
27/01/05
Auteur(s) :
E. BOYERE, L. VIVAN Clé
:
U4.84.02-F Page
: 3/8
3 Opérandes
3.1 Opérande
RESU_GENE
RESU_GENE = tg
Résultat d'un calcul transitoire par recombinaison modale, produit par l'opérateur
DYNA_TRAN_MODAL [U4.53.21].
3.2 Mot
clé
CHOC
CHOC
Mot clé facteur permettant de spécifier que l'on souhaite une analyse des non-linéarités de choc
prises en compte dans le calcul transitoire modal.
3.2.1 Opérande
INST_INIT
INST_INIT = t
0
Instant de début du moyennage des signaux et d'analyse des chocs. Par défaut t correspond au
0
premier instant du calcul transitoire par recombinaison modale.
3.2.2 Opérande
INST_FIN
INST_FIN = t
1
Instant de fin du moyennage des signaux et d'analyse des chocs.
(t1 = 999. valeur par défaut),
3.2.3 Opérande
NB_BLOC
NB_BLOC = nb
Nombre de blocs temporels de découpage de l'intervalle [t , ] pour le moyennage des signaux
0 t1
(1 par défaut).
Le mot clé n'est pas utilisé pour l'option 'IMPACT'.
3.2.4 Opérande
SEUIL_FORCE
SEUIL_FORCE = s
Seuil caractérisant une phase de contact, ( fn > s contact établi)
(s = 0. valeur par défaut).
3.2.5 Opérande
DUREE_REPOS
DUREE_REPOS = d
d
durée minimale de repos caractérisant la fin d 'un choc,
t
fin de choc
t
,
+
f ,
c
[tc tc d]
s
n
d = 0.
valeur par défaut.
Manuel d'Utilisation
Fascicule U4.8- : Post traitement et analyses dédiées
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur POST_DYNA_MODA_T
Date :
27/01/05
Auteur(s) :
E. BOYERE, L. VIVAN Clé
:
U4.84.02-F Page
: 4/8
3.2.6 Opérande
OPTION
OPTION = / 'USURE'
/
'IMPACT'
Ce mot clé permet de choisir entre un post-traitement en vue d'un diagnostic d'usure (mot clé
'USURE') ou d'une meilleure connaissance des éventuels chocs survenus lors de l'analyse
transitoire (mot clé 'IMPACT'). Ce dernier post traitement est adapté au calcul des internes des
centrales REP voir [§4].
3.2.7 Opérande
NB_CLASSE
NB_CLASSE = nc
Nombre de classes que l'utilisateur veut distinguer lors de l'élaboration de l'histogramme décrivant
les maxima de forces d'impact. Par défaut, nc = 10.
3.3 Mot
clé
RELA_EFFO_DEPL
RELA_EFFO_DEPL
Mot clé facteur permettant une analyse des relations de non-linéarité effort-déplacement.
On archive dans la table résultat le nom de la relation au noeud observé, les instants d'analyse et
le maximum atteint par la composante observée au cours de l'analyse.
On détermine toutes les phases à comportement non linéaire et on archive dans la table résultat
pour chacune :
·
l'instant initial et l'instant final de l'intervalle de la phase non linéaire,
·
le maximum atteint et l'instant associé dans cet intervalle.
3.3.1 Opérande
NOEUD
NOEUD = noeu
Nom du noeud de la structure sur lequel la relation non linéaire à traiter a été définie.
3.3.2 Opérande
NOM_CMP
NOM_CMP = noncmp
Nom de la composante traitée au noeud noeu.
3.4 Opérande
INFO
INFO = imp, niveau des impressions
/ 1 aucune impression sur le fichier message (tout est stocké dans la table produite)
/ 2 impression du nombre de pas de temps de calcul
3.5 Opérande
TITRE
TITRE = titre
titre est le titre du calcul. Il sera imprimé en tête des résultats. Il est stocké dans le concept
table.
Manuel d'Utilisation
Fascicule U4.8- : Post traitement et analyses dédiées
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur POST_DYNA_MODA_T
Date :
27/01/05
Auteur(s) :
E. BOYERE, L. VIVAN Clé
:
U4.84.02-F Page
: 5/8
4
Vérification - Exécution
4.1
Pour le mot clé facteur CHOC
La valeur de INST_FIN est comparée à l'instant final t f du résultat
tran_gene. La valeur de
INST_FIN retenue est
(
min t , t
f
1) .
Si la valeur de INST_INIT t0 est supérieure à la valeur de INST_FIN, on s'arrête en erreur.
F choc
Bloc 1
Bloc 2
SEUIL_FORCE
t
t
0
1
d
INST_INIT
INST_FIN
DUREE_REPOS
Bloc 1 : 1 choc
Bloc 2 : 2 chocs
: signal pris en compte dans les statistiques
4.1.1 Exécution avec option 'IMPACT'
Pour chaque non linéarité de choc, on calcule et on archive dans la table résultat :
·
les chocs et grandeurs associées :
pour chaque choc, on a la valeur de l'instant où la force est maximale, la valeur de la force
maximale et de l'impulsion, la durée du choc, la vitesse d'impact et le nombre de rebonds,
·
les données globales du choc :
sur l'ensemble des chocs constatés, on précise : le maximum absolu de force de choc, la
valeur moyenne des maxima de forces de choc et l'écart type des extrema de force de chocs,
·
l'histogramme décrivant les max des forces d'impact :
on a le nombre de classes de l'histogramme, les valeurs de ses abscisses (force min et max
de chaque classe) et la densité de probabilité de la force maximale de chaque classe.
Manuel d'Utilisation
Fascicule U4.8- : Post traitement et analyses dédiées
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur POST_DYNA_MODA_T
Date :
27/01/05
Auteur(s) :
E. BOYERE, L. VIVAN Clé
:
U4.84.02-F Page
: 6/8
4.1.2 Exécution avec l'option 'USURE'
Pour chaque non linéarité de choc, on calcule et on archive dans la table résultat :
·
les valeurs moyennes min, max, écart-type, RMS des déplacements relatifs des noeuds de
chocs dans leur repère global,
·
les valeurs moyennes et RMS (sur le temps de choc et temps total) ainsi que min et max des
forces normales et tangentielles de choc,
·
le nombre de chocs moyens sur chaque liaison de choc, le temps de choc moyen, le temps
de rebond moyen,
·
la puissance d'usure moyenne calculée au sens d'ARCHARD [bib1] :
1 T
P
=
F (t).V (t) dt
usure
T
N
T
0
4.2
Pour le mot clé facteur RELA_EFFO_DEPL
On vérifie que le noeud noeu correspond à une relation non linéaire.
On imprime à l'aide de la commande IMPR_TABLE :
---------------------------------------------------------------------------
ASTER 4.03.15 CONCEPT TABLE CALCULE LE 03/02/98 A 12:08:19 DE TYPE
TABL_POST_DYNA
RELATION
NOEUD
CMP
PHASE
INST_INIT
INST_FIN
MAXI
INST_MAXI
NON_LIN
N100
DRY
1
1.24000E+00
1.26600E+00
1.80716E-05 1.25200E+00
NON_LIN
N100
DRY
2
2.02600E+00
2.09000E+00
-2.00433E-05 2.05800E+00
NON_LIN
N100
DRY
3
3.00000E+00
3.04000E+00
-1.89110E-05 3.02200E+00
NON_LIN
N100
DRY
4
3.10400E+00
3.20400E+00
3.50715E-05 3.15400E+00
NON_LIN
N100
DRY
5
3.26000E+00
3.33600E+00
-2.91359E-05 3.30000E+00
NON_LIN
N100
DRY
6
3.41400E+00
3.43200E+00
1.73099E-05 3.42400E+00
.
Manuel d'Utilisation
Fascicule U4.8- : Post traitement et analyses dédiées
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur POST_DYNA_MODA_T
Date :
27/01/05
Auteur(s) :
E. BOYERE, L. VIVAN Clé
:
U4.84.02-F Page
: 7/8
5 Exemples
d'utilisation
5.1
Exemple avec l'option IMPACT : lancer d'un oscillateur avec choc
On présente ici le test SDND101. Il s'agit de calculer la réponse d'un système masse-ressort de
vitesse initiale V0 non nulle, pouvant impacter sur une butée à choc. Le jeu initial est nul.
K
m
Kchoc
V0
On compare les valeurs des instants de force maximale, valeur de force maximale, durée du temps de
choc, valeur de l'impulsion et de la vitesse d'impact ainsi que le nombre d'impacts élémentaires pour
les deux premières oscillations du système aux valeurs analytiques.
Fichier de commande
POURSUITE
()
# Le systeme masse-ressort est lache avec une vitesse initiale V0
#
VITEPHYS = AFFE_CHAM_NO (
MAILLAGE= poutre,
GRANDEUR=
'DEPL_R',
CHAM_NO
=
vectass,
AFFE=
(
TOUT=
'OUI',
NOM_CMP=
(
'DX'
),
VALE_R=
V0
)
)
#
VITINI
=
PROJ_VECT_BASE
( BASE= MODES, VECT_ASSE= VITEPHYS,
NUME_DDL_GENE=
NUMEGE,
TYPE_VECT=
'VITE')
#
MASSEGEN = PROJ_MATR_BASE ( BASE= MODES,
NUME_DDL_GENE=
NUMEGE,
MATR_ASSE=
MATRMASS,
)
#
RIGIDGEN = PROJ_MATR_BASE ( BASE= MODES,
NUME_DDL_GENE=
NUMEGE,
MATR_ASSE=
MATRRIGI,
)
#
PLANZ = DEFI_OBSTACLE ( TYPE= 'PLAN_Z' )
PLANY = DEFI_OBSTACLE ( TYPE= 'PLAN_Y' )
#
DYNAMODA
=
DYNA_TRAN_MODAL
( METHODE = 'EULER',
MASSE_GENE
=
MASSEGEN,
RIGI_GENE
=
RIGIDGEN,
AMOR_REDUIT=
0.,
ETAT_INIT=_F(VITE_INIT_GENE=VITINI),
INCREMENT=_F(
INST_INIT= 0.,
INST_FIN=
0.5,
PAS=
0.0005
),
Manuel d'Utilisation
Fascicule U4.8- : Post traitement et analyses dédiées
HT-66/05/004/A
Code_Aster ®
Version
7.4
Titre :
Opérateur POST_DYNA_MODA_T
Date :
27/01/05
Auteur(s) :
E. BOYERE, L. VIVAN Clé
:
U4.84.02-F Page
: 8/8
CHOC=
_F(
NOEUD_1=
'NO1',
OBSTACLE=
PLANZ,
ORIG_OBST=
(-1.,
0.,
0.),
NORM_OBST=
(0.
0.
1.),
JEU=
1.00,
RIGI_NOR=
1000000.,
RIGI_TAN=
0.,
COULOMB=
0.
)
)
# Post-traitement des efforts de choc
TT
= POST_DYNA_MODA_T
( RESU_GENE = DYNAMODA,
CHOC=
_F(
INST_INIT
= 0.,
INST_FIN
=
0.495,
SEUIL_FORCE
=
0.,
DUREE_REPOS
=
0.,
OPTION
=
'IMPACT',
NB_CLASSE
=
8
),
INFO
=
1
)
IMPR_TABLE ( TABLE = TT,
FILTRE
=
(
NOM_PARA=
'CALCUL',
VALE_K
=
'GLOBAL'
),
NOM_PARA=
(
'NOEUD',
'F_MAX_ABS',
'F_MAX_MOY',
'F_MAX_ETYPE'
))
6 Bibliographie
[1]
ARCHARD : The wear of metals under unlubrificated conditions Proc-Roy-Soc. (1956).
Manuel d'Utilisation
Fascicule U4.8- : Post traitement et analyses dédiées
HT-66/05/004/A
Document Outline