Opérateur PROJ_MESU_MODAL
Contents
1. Ce mot cl renseigne le nom du n ud du mod le num rique que l on souhaite associer au n ud d observation no1 Manuel d utilisation Fascicule u4 73 Traitement de mesure Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster dl Titre Op rateur PROJ_MESU_MODAL Date 01 03 2013 Page 6 10 Responsable Harinaivo ANDRIAMBOLOLONA Cl UA 73 01 R vision 10613 3 5 Mot cl facteur RESOLUTION On d finit ici la m thode de r solution utiliser et les param tres associ s cette m thode 3 5 1 Op rande METHODE METHODE LU SVD On propose la m thode LU d composition en LU Lower Upper et la m thode SVD d composition en valeurs singuli res pour le calcul de la matrice inverse Pour la m thode SVD le nombre de valeurs singuli res prendre en compte d pend de la valeur de eps que l utilisateur renseigne sous l op rande EPS Par d faut on adopte la m thode LU 3 5 2 Op rande EPS Ce mot cl est utilis si on choisit la m thode SVD 0 EPS 0 eps Ce mot cl donne la valeur partir de laquelle une valeur singuli re est consid r e comme nulle Il d termine ainsi le nombre de valeurs singuli res exploiter lors de la r solution Un eps gal z ro signifie que toutes les valeurs singuli res sont prendre en compte eps gal 1 signifie que l on ne consid re2. Fascicule u4 73 Traitement de mesure Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl htmi
3. RESOLUTION ei ETHODE 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 Manuel d utilisation Fascicule u4 73 Traitement de mesure Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster Fo Titre Op rateur PROJ MESU MODAL Date 01 03 2013 Page 10 10 Responsable Harinaivo ANDRIAMBOLOLONA Cl U4 73 01 R vision 10613 6 Bibliographie _ C VARE Extrapolation de r sultats de mesures exp rimentales sur un mod le num rique en dynamique Sp cification des d veloppements dans le Code Aster Note EDF DER HP 54 98 063 B S AUDEBERT Evaluation comparative de diff rentes m thodes d inversion Note EDF DER HP 62 93 036 A TIKHONOV V ARSENINE M thodes de r solution de probl mes mal pos s Ed Mir 1976 M BONNET Traitement num rique de probl mes inverses de source en acoustique lin aire Contrat EDF Convention P55L08 1E5240 A TARANTOLA Inverse problem theory Methods for data fitting and model parameter estimation Elsevier 1987 Manuel d utilisation
4. rande REGU coito a nas 6 3 5 40p randes COEF_PONDER et COEF PONDER Fi 7 4Phase de v rification et d ex cution 7 4 1Calcul de la base d expansion restreinte aux degr s de libert mesur s 7 4 2Calcul des coordonn es g n ralis es E EEEE E AEAEE EEEE EEE eEnnt 7 5Exemple d utilisation de PROJ_MESU_MODAL ooocccccccnncnccnnnonnnonnnnconcnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnannno 8 GBIBlIOg aphie iii ia 10 Manuel d utilisation Fascicule u4 73 Traitement de mesure Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl htmi Code Aster Version default Titre Op rateur PROJ_MES U MODAL Date 01 03 2013 Page 3 10 Responsable Harinaivo ANDRIAMBOLOLONA Cl U4 73 01 R vision 10613 2 Syntaxe repgen gene PROJ MESU MODAL MODELE CALCUL _F BASE base mode meca MODELE mocalc modele MODELE MESURE _F MESURE mesure dyna trans dyna harmo A mode meca mode meca cl Y MODELE mostru modele 0 NOM CHAM DEPL DEFAUT VITE ACCE SIEF_NOEU EPSI_NOEU 0 NOM PARA lpara 1 Kn CORR MANU _F NOEU MESURE nol noeud NOEU CALCUL no2 noeud
5. rateur PROJ MESU MODAL Date 01 03 2013 Page 7 10 Responsable Harinaivo ANDRIAMBOLOLONA Cl U4 73 01 R vision 10613 Il est d conseill d utiliser ce mot cl lorsque le champ mesur mesure est de type mode meca 3 5 4 Op randes COEF _PONDER et COEF PONDER F Ce mot cl correspond au poids affect l information a priori Il est utilis si on applique une r gularisation la solution n 0 COEF PONDER coef Liste des coefficients de pond ration sur la solution a priori m thode de r gularisation de Tikhonov bib3 COEF PONDER F coef f Liste des fonctions de pond ration sur la solution a priori m thode de r gularisation de Tikhonov Les variables de ces fonctions sont les m mes que celles du champ mesur mesure Si le nombre de coefficients ou de fonctions de pond ration donn est inf rieur au nombre de vecteurs de base utilis s dans la base d expansion les coefficients ou fonctions de pond ration des vecteurs suppl mentaires sont pris gaux au dernier coefficient ou la derni re fonction de la liste 4 Phase de v rification et d ex cution 4 1 Calcul de la base d expansion restreinte aux degr s de libert mesur s Dans un premier temps le maillage de la mesure est projet sur le maillage du mod le num rique On d termine ensuite la participation des n uds du mod le num rique pour chaque n ud de mesure via la fonction de forme de l l ment qui c
6. Code Aster Pa Titre Op rateur PROJ MESU MODAL Date 01 03 2013 Page 1 10 Responsable Harinaivo ANDRIAMBOLOLONA Cl U4 73 01 R vision 10613 Op rateur PROJ MESU MODAL 1 But Extrapoler des mesures exp rimentales sur un mod le num rique en dynamique Les donn es exp rimentales peuvent tre des d placements des vitesses des acc l rations des d formations ou des contraintes Elles sont d finies en tant que fonction du temps ou de la fr quence ou sous forme de liste Il s agit de calculer les coordonn es g n ralis es de la mesure relatives une base d expansion d finie sur le mod le num rique Cette base d expansion d form es contraintes ou d formations est calcul e au pr alable suivant un concept de type mode meca Les vecteurs de base sont ensuite restreints aux degr s de libert mesur s L association spatiale entre les points de mesure et les n uds du maillage num rique peut se r aliser manuellement ou et automatiquement Le calcul des coordonn es g n ralis es est r alis e par r solution d un probl me de minimisation de type moindres carr s ventuellement r gularis suivant la m thode de Tikhonov Est applicable tout type de mod le 1D 2D et 3D Produit une structure de donn es de type tran gene harm gene ou mode gene Manuel d utilisation Fascicule u4 73 Traitement de mesure Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL ht
7. RESOLUTION E METHODE LU DEFAUT SVD Si METHODE SVD alors 0 EPS 0 DEFAUT eps R Q REGUL f NON DEFAUT NORM_MIN TIK_RELA Si REGUL NON alors Q COEF PONDER O DEFAUT o w IR COEF PONDER F w f 1 fonction Si mesure dyna trans alors repgene tran gene Si mesure dyna harmo alors repgene harm gene Si mesure mode meca alors repgene mode gene Si mesure mode meca c alors repgene mode gene Manuel d utilisation Fascicule u4 73 Traitement de mesure Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl html Code Aster a Titre Op rateur PROJ MESU MODAL Date 01 03 2013 Page 4 10 Responsable Harinaivo ANDRIAMBOLOLONA Cl U4 73 01 R vision 10613 3 Op randes 3 1 Mot cl facteur MODELE CALCUL Ce mot cl facteur regroupe les caract ristiques du mod le num rique sur lequel on veut extrapoler la mesure Il ne doit appara tre qu une seule fois 3 1 1 Op rande MODELE MODELE mocalc Nom du mod le num rique sur lequel est construite la base d expansion 3 1 2 Op rande BASE BASE base Nom de la base d expansion Cette base est de type mode meca Ce concept a ventuellement t enrichi via la commande CALC CHAMP par les champs de d formations et ou de contraintes modales calcul s aux n uds 3 2 Mot cl fac
8. ante a le m me poids lors de l inversion Op rande NOM_PARA 0 NOM PARA lpara Liste des noms symboliques des param tres des donn es mesur es que l on souhaite transmettre au mod le g n ralis Mot cl facteur CORR MANU Ce mot cl facteur permet l utilisateur de d finir surcharger manuellement la correspondance entre le noeud d observation et le noeud analogue du mod le num rique Ce mot cl facteur est facultatif mais il peut appara tre autant de fois que n cessaire Par contre les op randes sous ce mot cl facteur vont par paire Un NOEU MESURE doit avoir son NOEU CALCUL associ Si ce mot cl facteur est absent l association spatiale entre les points de mesure et les noeuds du maillage num rique s effectue automatiquement en utilisant la fonction de forme de l l ment du mod le num rique pour d terminer la valeur du champ sur le point de mesure Op rande NOEU MESURE NOEU MESURE nol Ce mot cl renseigne le nom du n ud d observation que l on souhaite associer au noeud du mod le num rique no2 Dans certains cas le fichier de maillage associ la mesure est au format universel format Ideas on ne peut donc pas savoir a priori le nom Aster associ au n ud ll est donc n cessaire dans ce cas de lire le maillage issu de PRE IDEAS par LIRE MAILLAGE afin de pouvoir affecter le nom du n ud Op rande NOEU CALCUL NOEU CALCUL no2
9. ontient le n ud de mesure La correspondance obtenue entre les n uds est fournie dans le fichier MESSAGE de l tude Aster Le deuxi me traitement consiste calculer la composante du champ base d expansion au n ud de mesure suivant les degr s de libert mesur s 4 2 Calcul des coordonn es g n ralis es La solution de l quation de minimisation est donn e par 1_ T T n 0 D rum B num Di exp 0 1 AT 7 HT n i PB Oi de FAN priori Avec n i coordonn es g n ralis es pour le num ro d ordre i ti ou fi Gel mesure au num ro d ordre i base d expansion restreinte aux degr s de libert de mesure gt ali coefficients permettant de sp cifier le poids affect l information a priori au num ro d ordre i Ces variables ou fonctions sont d finies par l utilisateur dans les op randes COEF PONDER OU COEF PONDER F du mot cl facteur RESOLUTION Elles sont introduites sous la forme d une liste de r els ou de fonctions et correspondent terme terme chaque vecteur de la base d expansion retenue Manuel d utilisation Fascicule u4 73 Traitement de mesure Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copylett fdl html Code Aster dl Titre Op rateur PROJ MESU MODAL Date 01 03 2013 Page 8 10 Responsable Harinaivo ANDRIAMBOLOLONA Cl U4 73 01 R vision 10613 Sui
10. que la plus grande valeur singuli re Par d faut on choisit EPS 0 3 5 3 Op rande REGUL REGUL NON NORME_MIN TIK RELA REGUL permet de pr ciser la m thode de r gularisation que l on veut utiliser Par d faut on ne rajoute pas de r gularisation pas de contrainte suppl mentaire sur la solution REGUL NON Actuellement deux types de r gularisation sont disponibles norme minimale REGUL NORM MIN ou Tikhonov d ordre O et Tikhonov relatif REGUL TIK_RELA On cherche minimiser pour chaque num ro d ordre du champ mesur la fonctionnelle suivante par rapport n 7 2 2 ldap Diin n T a N N prioril avec d n coordonn es g n ralis es relatives la base d expansion e dexp Champ mesur suivant les degr s de libert d observation D m base d expansion restreinte aux degr s de libert d observation a coefficients de pond ration permettant de sp cifier le poids affect l information a priori sur la solution Suivant la m thode utilis e les param tres de la fonctionnelle pr c dente se d clinent comme suit Sans r gularisation x 0 Norme minimale NORM_MIN N priori Tikhonov relatif TIK_RELA N priori Solution trouv e au num ro d ordre pr c dent Manuel d utilisation Fascicule u4 73 Traitement de mesure Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster Mo Titre Op
11. t l affecte au mod le correspondant Manuel d utilisation Fascicule u4 73 Traitement de mesure Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fal html Code Aster default Titre Op rateur PROJ MESU MODAL Date 01 03 2013 Page 9 10 Responsable Harinaivo ANDRIAMBOLOLONA Cl U4 73 01 R vision 10613 mesure LIRE RESU FORMAT IDEAS DS58 UNITE 33 MAILLAGE mailmesu TYPE RE DYNA TRANS NOM CHAM SIEF NOEU Y em D finition de la base d expansion La base d expansion doit tre du type mode meca Ce concept peut tre issu de MODE ITER SIMULT 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 ESURE mesure NOM CHAM SIEF NOEU CORR MANU _F NOEU MESURE nol NOEU CALCUL no2
12. teur MODELE MESURE Ce mot cl facteur regroupe les informations sur le champ mesur observ que l on souhaite extrapoler sur le mod le num rique Il ne doit appara tre qu une seule fois 3 2 1 Op rande MODELE MODELE mostru Nom du mod le associ l observation 3 2 2 Op rande MESURE MESURE mesure Nom du champ mesur Ce mot cl d termine le type de concept produit par l op rateur PROJ MESU MODAL Si mesure est de type dyna trans le concept produit est de type tran gene Si mesure est de type dyna harmo le concept produit est de type harm gene Si mesure est de type mode meca ou mode meca c le concept produit est de type mode gene Manuel d utilisation Fascicule u4 73 Traitement de mesure Copyright 2015 EDF R amp D Document diffus sous licence GNU FDL http www gnu org copyleft fdl htmi Code Aster Mo Titre Op rateur PROJ_MESU_MODAL Date 01 03 2013 Page 5 10 Responsable Harinaivo ANDRIAMBOLOLONA Cl U4 73 01 R vision 10613 3 2 3 Op rande NOM CHAM 3 3 3 4 3 4 1 3 4 2 DEPL VITI ACC SIEF NOEU EPSI NOEU 0 NOM CHAM DEFAUT 0 El Ce mot cl permet de choisir le ou les champs lire et extrapoler Les composantes du champ consid r es sont celles qui ont t mesur es observ es et lues dans mesure Pour l instant on ne affecte pas de coefficient de pond ration sur les diff rentes composantes du champ chaque compos
13. tp www gnu org copylett fdl html Code Aster dl Titre Op rateur PROJ MESU MODAL Date 01 03 2013 Page 2 10 Responsable Harinaivo ANDRIAMBOLOLONA Cl U4 73 01 R vision 10613 Table des mati res ITP TER e a a A a A 1 PAS E E T E A E E E E T EE ATNEN 3 S181 S E RA 4 3 1Mot cl facteur MODELE_CALCUb o ooooocccicicccocooonononononcnncnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnicinnncnns 4 3 1 10p rand MODELE ciiai a a td anna 4 3 1 20p rand BASE ut lee tn et Moine Ar A LR M ne cn a a Er 4 3 2Mot cl facteur MODELE_MESURE ooooococccococcccccccccnnnncnccccccncnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnrnnnnnnnnaneniincnns 4 3 2 10p rande MODELE iaiia AA te entente rte 4 3 2 20p rande MESURE octal lt 4 3 2 30p rande NOM_CHAM nn nn a aai a a a a a a aaa eta 5 3 30p rande NOM PARA iii it a 5 3 4Mot cl facteur CORR_MAN U ooooooocccccccccccccococnnnnccccncnnnnnnnnnnnnn nro nnnnnnnnnnnnnnnnnnnnnnnnnnnrinnnnnnnnnnnnncanicines 5 3 4 10p rande NOEU_MESURE ooooooooccccccccconcncnncccccnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnniniss 5 3 4 20p rande NOEU_CALCUbocoooonccccocoooconccnccccnnnncnnnnnncnnncnncnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnnncnnnnnnnnins 5 3 5Mot cl facteur R SOLUTION nn nn nro nnnnnnnnnnnnnnnnnnnrnnnnnnnnnnnnnnnnnennnnnes 6 3 5 10p rande METHODE iiiiiiieeeeeeeeeeeeereeeeeeeeeeeeeneeneneeeeneeee 6 3 5 20p rande EPS inner bn nent cd dd 6 3 5 30p
14. uds de mesure peuvent avoir t d finies Elles nont 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 l deas Pour le transformer en format Aster on utilise l op rateur PRE _IDEAS PRE IDEAS UNITE IDEAS 19 UNITE MAILLAGE 2173 mailmesu LIRE MAILLAGE UNITE 21 e 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 gt 3 degr s de libert par n ud DX DY et DZ ou affectation d une mod lisation DIS TR discret en translation rotation gt 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 e 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 e
15. vant la m thode utilis e les param tres pr c dents se d clinent comme suit Sans r gularisation amp 0 Norme minimale NORM MIN N priori Tikhonov relatif TIK RELA N priori Ni 1 Remarque 1 Si un coefficient de pond ration est n gatif le traitement s arr te en erreur fatale Remarque 2 Si tous les coefficients de pond ration sont nuls pour un num ro d ordre donn et que le nombre de mesures est strictement inf rieur au nombre de vecteurs de base un message d alarme est mis pour pr venir du risque de matrice singuli re en effet dans ce cas il n y a pas unicit de la solution A l issue du calcul les coordonn es g n ralis es identifi es sont d riv es afin de calculer les vitesses et les acc l rations correspondantes Le r sultat de l inversion est un concept de type tran gene harm gene OU mode gene 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 e 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 l deas ou prochainement MED Des mailles reliant les n
Download Pdf Manuals
Related Search