5. Exemple d’utilisation de PROJ_MESU_MODAL#
Pour les exemples d’utilisation, il est vivement conseillé de se rapporter aux cas tests SDLD104 et SDLV122.
On présente dans ce paragraphe les différentes étapes pour l’expansion de la mesure sur le modèle numérique.
Lecture du maillage constitué des points de mesure:
Cette opération a pour objectif de lire la position des points de mesure dans un fichier de type maillage. Le format de ce fichier doit être lisible par le Code_Aster (format GIBI, universel (I-deas) ou prochainement MED). Des mailles reliant les nœuds de mesure peuvent avoir été définies. Elles n’ont évidemment aucune signification physique mais permettront éventuellement de visualiser les résultats lors de la phase de post-traitement.
Dans la majorité des cas, le maillage est issu d’un code de mesure expérimentale qui fournit un fichier au format universel (format I-deas). Pour le transformer en format Aster, on utilise l’opérateur PRE_IDEAS.
PRE_IDEAS(UNITE_IDEAS = 19, UNITE_MAILLAGE = 21,)
mailmesu = LIRE_MAILLAGE ( UNITE = 21,)
Affectation d’un modèle mécanique au maillage:
Cette opération a pour objectif de définir le modèle du support des nœuds du maillage constitué des points de mesure. Deux cas peuvent être envisagés: affectation d’une modélisation DIS_T (discret en translation => 3 degrés de liberté par nœud: DX, DY et DZ) ou affectation d’une modélisation DIS_TR (discret en translation - rotation => 6 degrés de liberté par nœud: DX, DY, DZ, DRX, DRY et DRZ) dans le cas où des mesures de rotation sont réalisées.
modlmesu = AFFE_MODELE ( MAILLAGE = mailmesu,
AFFE = _F( GROUP_NO = “noeumesu”, MODELISATION = ‘DIS_T’, PHENOMENE = ‘MECANIQUE’,),
)
Lecture de la mesure:
La mesure peut être lue via l’opérateur LIRE_RESU. Cet opérateur permet de lire un fichier au format universel (dataset 58). Il récupère la composante du champ observé et l’affecte au modèle correspondant.
mesure = LIRE_RESU ( FORMAT = 'IDEAS_DS58',
UNITE = 33,
MAILLAGE = mailmesu,
TYPE_RESU = 'DYNA_TRANS',
NOM_CHAM = 'SIEF_NOEU',)
Définition de la base d’expansion:
La base d’expansion doit être du type mode_meca. Ce concept peut être issu de CALC_MODES ou de DEFI_BASE_MODALE.
Si on veut extrapoler un champ de déformation ou de contrainte, la base doit être enrichie, via la commande CALC_CHAMP, par les champs de déformation ou de contrainte calculés aux nœuds.
Calcul des coordonnées généralisées:
Le calcul des coordonnées généralisées relatives à la base d’expansion est assuré par l’opérateur PROJ_MESU_MODAL.
repgene = PROJ_MESU_MODAL (
MODELE_CALCUL = _F( MODELE = modlcalc,
BASE = base,),
MODELE_MESURE = _F( MODELE = modlmesu,
MESURE = mesure,
NOM_CHAM = 'SIEF_NOEU',),
CORR_MANU = _F( NOEU_MESURE = 'no1',
NOEU_CALCUL = 'no2',),
RESOLUTION = _F( METHODE = 'SVD',
EPS = 1.E-4,),
)
Expansion sur le modèle numérique:
Cette expansion consiste à calculer sur tous les nœuds du modèle numérique, le champ compatible avec le champ observé sur les degrés de liberté mesurés. Cette expansion est réalisée par la commande REST_GENE_PHYS.
reponse = REST_GENE_PHYS ( RESU_GENE = repgene,
TOUT_CHAM = 'OUI',)