6. Post-traitement#

\[\newcommand{\vector}[1]{\underline{#1}} \newcommand{\vectorZero}{\vector{0}} \newcommand{\tensTwo}[1]{\boldsymbol{#1}} \newcommand{\tensTwoZero}{\tensTwo{0}} \newcommand{\tensTwoUnit}{\tensTwo{I}} \newcommand{\tensFour}[1]{\mathbb{#1}} \newcommand{\inverse}[1]{{#1}^{-1}} \newcommand{\transpose}[1]{{#1}^{T}} \newcommand{\inverseTranspose}[1]{ {#1}^{-T}} \newcommand{\vectorCmp}[2]{#1_{\left(#2\right)}} \newcommand{\vectorCmpCO}[2]{#1_{\left(#2\right)}} \newcommand{\vectorCmpCT}[2]{#1^{\left(#2\right)}} \newcommand{\tensTwoCmpCO}[3]{#1_{\left(#2 #3\right)}} \newcommand{\tensTwoCmpCT}[3]{#1^{\left(#2 #3\right)}} \newcommand{\tensTwoInva}[2]{I^{#1}_{#2}} \newcommand{\tensTwoDevia}[1]{\boldsymbol{\tilde{#1}}} \newcommand{\tensFourCmpCO}[5]{#1_{\left(#2 #3 #4 #5\right)}} \newcommand{\tensFourCmpCT}[5]{#1^{\left(#2 #3 #4 #5\right)}} \newcommand{\domain}{\Omega} \newcommand{\domainRefe}{\domain^{0}} \newcommand{\domainCurr}{\domain^{t}} \newcommand{\bound}{\partial\domain} \newcommand{\boundRefe}{\bound^{0}} \newcommand{\boundCurr}{\bound^{t}} \newcommand{\boundN}{\bound_{\textrm{N}}} \newcommand{\boundNRefe}{\boundN^{0}} \newcommand{\boundNCurr}{\boundN^{t}} \newcommand{\boundD}{\bound_{\textrm{D}}} \newcommand{\boundDRefe}{\boundD^{0}} \newcommand{\boundDCurr}{\boundD^{t}} \newcommand{\normal}{\vector{n}} \newcommand{\normalRefe}{\normal^{0}} \newcommand{\normalCurr}{\normal^{t}} \newcommand{\posi}{\vector{x}} \newcommand{\posiRefe}{\posi^{0}} \newcommand{\posiCurr}{\posi^{t}} \newcommand{\disp}{\vector{u}} \newcommand{\dispVirt}{\delta\vector{v}} \newcommand{\funcTransfor}{\vector{\varphi}^{t}} \newcommand{\gradTransfor}{\tensTwo{F}} \newcommand{\jacobTransfor}{J} \newcommand{\posiIncr}{\vector{dx}} \newcommand{\posiIncrRefe}{\posiIncr^{0}} \newcommand{\posiIncrCurr}{\posiIncr^{t}} \newcommand{\strainCmp}{\varepsilon} \newcommand{\strain}{\tensTwo{\strainCmp}} \newcommand{\EGLCmp}{E} \newcommand{\ECGDroiteCmp}{C} \newcommand{\ECGGaucheCmp}{B} \newcommand{\ELOGCmp}{{E}_{ln}} \newcommand{\EGL}{\tensTwo{\EGLCmp}} \newcommand{\ECGDroite}{\tensTwo{\ECGDroiteCmp}} \newcommand{\ECGGauche}{\tensTwo{\ECGGaucheCmp}} \newcommand{\ELOG}{\tensTwo{\ELOGCmp}} \newcommand{\EDilCmp}{C^{\star}} \newcommand{\EDil}{\tensTwo{\EDilCmp}} \newcommand{\divTensTwo}[1]{\vector{\nabla} {\cdot} #1 } \newcommand{\divVector}[1]{\vector{\nabla} {\cdot} #1 } \newcommand{\gradScal}[1]{\vector{\nabla} #1 } \newcommand{\gradVector}[1]{\vector{\nabla} \times #1 } \DeclareMathOperator{\rand}{rand} \DeclareMathOperator{\round}{round} \DeclareMathOperator{\trace}{trace} \DeclareMathOperator{\rank}{rang} \DeclareMathOperator{\dev}{dev} \DeclareMathOperator{\cot}{cot} \newcommand{\metric} {\tensTwo{g}} \newcommand{\vectorBaseCV}[1]{\vector{g_{#1}}} \newcommand{\vectorBaseCT}[1]{\vector{g^{#1}}} \newcommand{\metricRefe} {\metric^{t}} \newcommand{\metricCurr} {\metric^{0}} \newcommand{\discVect}[1]{\lbrace #1 \rbrace} \newcommand{\discVectLigne}[1]{\langle #1 \rangle} \newcommand{\discVectZero}{\discVect{0}} \newcommand{\discVectCmp}[2]{\lbrace #1 \rbrace_{(#2)}} \newcommand{\discVectLigneCmp}[2]{\langle #1 \rangle_{(#2)}} \newcommand{\discMatr}[1]{\left[ #1 \right]} \newcommand{\discMatrZero}{\discMatr{0}} \newcommand{\discMatrCmp}[3]{\left[ #1 \right]_{(#2#3)}} \newcommand{\onQuadPoint}[2]{{#1}_{#2}} \newcommand{\onNode}[2]{{#1}^{#2}} \newcommand{\quadOrder}{k_{Q}} \newcommand{\nbQuadPoint}{n_Q} \newcommand{\quadWeight}[1]{\omega_{#1}} \newcommand{\quadPointIndex}{i_{\textrm{pg}}} \newcommand{\nodeIndex}{i_{\textrm{no}}} \newcommand{\pres}{p} \newcommand{\presRefe}{\pres^{0}} \newcommand{\presCurr}{\pres^{t}} \newcommand{\stressCmp}{\sigma} \newcommand{\stress}{\tensTwo{\stressCmp}} \newcommand{\stressPKTwoCmp}{S} \newcommand{\stressPKTwo}{\tensTwo{\stressPKTwoCmp}} \newcommand{\work}[1]{W^{#1}} \newcommand{\enerInterne}{\Psi_i} \newcommand{\yieldStress}{\sigmaCmp_{Y}} \newcommand{\youngModulus}{E} \newcommand{\youngModulusCplx}{E^{\star}} \newcommand{\poissonCoef}{\nu} \newcommand{\poissonCoefCplx}{\nu^{\star}} \newcommand{\shearModulus}{G} \newcommand{\bulkModulus}{K} \newcommand{\modulusTangent}{\tensFour{K}} \newcommand{\anglDila}{\psi} \newcommand{\anglFric}{\varphi} \newcommand{\cohesion}{c} \newcommand{\presCapi}{\pres_{c}} \newcommand{\presEau}{\pres_{w}} \newcommand{\presAir}{\pres_{a}} \newcommand{\shapeFunc}{\Phi} \newcommand{\shapeDFunc}{B} \newcommand{\measLine}{\, \mathrm{d} l} \newcommand{\measDomain}{\, \mathrm{d} \domain} \newcommand{\measBound}{\, \mathrm{d} \Gamma} \newcommand{\normEucl}[1]{\Vert{#1} \Vert } \newcommand{\valeAbs}[1]{\vert{#1} \vert } \newcommand{\derivee}[2]{#1_{, #2}} \newcommand{\soundSpeed}{c} \newcommand{\soundSpeedComp}{c_{P}} \newcommand{\soundSpeedCisa}{c_{S}} \newcommand{\posiTang}{\posi_{\tau}} \newcommand{\posiNorm}{{x}_{n}} \newcommand{\dispTang}{\disp_{\tau}} \newcommand{\dispNorm}{{u}_{n}} \newcommand{\basisVector}[1]{\vector{e}_{#1}} \newcommand{\sigmVector}{\vector{t}} \newcommand{\sigmVectorTang}{\sigmVector_{\tau}} \newcommand{\sigmVectorNorm}{t_{n}} \newcommand{\lameLambda}{\lambda} \newcommand{\lameMu}{\mu} \newcommand{\temp}{T} \newcommand{\therCp}{C_{p}} \newcommand{\therFlux}{\vector{q}} \newcommand{\therCondTher}{\lambda} \newcommand{\therSource}{s} \newcommand{\tempVirt}{\delta T} \newcommand{\beamArea}{S} \newcommand{\beamLength}{L} \newcommand{\setC}{\mathbb{C}} \newcommand{\setR}{\mathbb{R}}\]

Dans le cas où TYPE_RESU est différent de “FICHIER”, “FICHIER_TEMPS” et “TABLE_CONTROL”, les fichiers résultats de Miss3D sont post-traités par CALC_MISS afin de fournir la réponse harmonique ou temporelle de la structure, ou des évolutions des grandeurs caractéristiques (déplacement, vitesse, accélération, spectre d’oscillateur) en certains points de post-traitement.

6.1. La définition du chargement#

Il y a trois façon différentes d’imposer le chargement à l’aide de CALC_MISS:

  • On peut appliquer la force sismique équivalente évaluée au préalable avec Miss3D en multipliant cette force par une fonction multiplicatrice. Il s’agit de l’option EXCIT_FORC

  • On peut appliquer un chargement générique sur le modèle qui a été produit préalablement avec AFFE_CHAR_MECA à l’aide de l’option EXCIT_GENE

  • Enfin, on peut imposer un chargement mono-appui à la base du modèle via l’option EXCIT_MONO

Avertissement

Compte-tenu de la différente nature des chargements, on ne peut imposer qu’un seul type pour calcul:

  • En effet, l’option EXCIT_FORC produit un résultat en repère absolu, tandis que EXCIT_MONO produit un résultat en repère relatif. Le mélange de chargement n’est donc pas raisonnable

  • De l’autre coté, l’option EXCIT_GENE permet de produit un chargement générique (que ça soit un déplacement imposé dans le repère absolu ou un chargement type mono-appui)

6.1.1. Opérande TYPE_EXCIT#

Si on est en présence d’un chargement EXCIT_MONO ou EXCIT_FORC, il faut declarer la nature du chargement fourni par l’utilisateur. Deux choix sont possibles:

  • “INST” pour un chargement en temporel

  • “FREQ” pour un chargement en fréquentiel

6.1.2. Opérandes INTERPOL, PAS_INST, INST_FIN#

Dans le cas où l’utilisateur fournit un chargement temporel (TYPE_EXCIT=”INST”), on a la possibilité de ré-interpoler ce chargement (INTERPOL = “OUI”, par défaut “NON”). Il faudra alors fournir deux paramètres supplémentaires:

  • PAS_INST, c’est-à-dire le pas de temps d’interpolation

  • INST_FIN, c’est-à-dire l’instant final du chargement

Avertissement

Dans tous les cas, INST_FIN ne peut pas être plus grand que l’instant final de la fonction initiale. En ne connaissant pas la nature du chargement, on ne choisit pas la nature d’un éventuel prolongement. INST_FIN a donc plutôt un fonction de coupure du signal d’entrée .

Opérande MODELEe

Il s’agit du modèle de la structure (transmis à DYNA_VIBRA).

6.2. Calcul de la réponse harmonique ou temporelle de la structure#

On se trouve dans le cas TYPE_RESU=”HARM_GENE” (réponse harmonique) ou “TRAN_GENE” (réponse temporelle). On calcule alors la réponse harmonique de la structure au chargement fourni (EXCIT_GENE, EXCIT_MONO ou EXCIT_FORC).

Avertissement

Dans le cas “TRAN_GENE”, on effectue la restitution temporelle en utilisant l’opérateur REST_SPEC_TEMP (option PROL_ZERO).

6.2.1. Opérandes EXCIT_FORC, ACCE_X, ACCE_Y, ACCE_Z, DEPL_X, DEPL_Y, DEPL_Z#

L’option EXCIT_FORC est utilisée pour imposer le chargement via les forces sismiques équivalentes calculées au préalable avec Miss3D. Dans ce cas on fournira:

  • Soit un accélérogramme dans une ou plusieurs directions (ACCE_X, ACCE_Y, ACCE_Z)

  • Soit des déplacements imposés dans une ou plusieurs directions (DEPL_X, DEPL_Y, DEPL_Z)

Avertissement

En fonction de la nature du chargement il faut fournir soit des fonctions réelles (TYPE_EXCIT=”INST”) soit des fonctions complexes (TYPE_EXCIT=”FREQ”).

6.2.2. Opérandes EXCIT_HARMO, VECT_ASSE, FONC_MULT, COEF_MULT, FONC_MULT_C, COEF_MULT_C, PHAS_DEG=SIMP, PUIS_PULS#

EXCIT_HARMO permet de définir une charge générique dans CALC_MISS, charge qui sera fourni en entrée de DYNA_VIBRA. Pour plus de détails sur la définition d’une charge générique dans DYNA_VIBRA on peut se référer à l’option EXCIT dans la documentation u4.53.03.

6.2.3. Opérandes EXCIT_MONO, DIRECTION, ACCE#

EXCIT_MONO permet de définir un chargement mono-appui à la base du modèle. Pour cela, il faut specifier:

  • La direction (mot-clé DIRECTION) du mouvement d’entraînement (“X”, “Y” ou “Z”)

  • La fonction multiplicatrice (mot-clé ACCE). Cette fonction multiplicatrice peut être réelle (TYPE_EXCIT=”INST”) ou complexe (TYPE_EXCIT=”FREQ”)

6.3. Calcul des évolutions en certains points#

On est ainsi dans le cas TYPE_RESU=”TABLE”. Dans ce cas, on calcule la réponse harmonique de la structure à une accélération unitaire (dans la ou les directions demandées). Puis, pour chaque chargement, on recombine en chaque lieu de post-traitement \(M\) les contributions fréquentielles unitaires (équation (6.1)):

(6.1)#\[{u}_{M}(f)={u}_{x}\times FFT(a_x)+{u}_{y}\times FFT(a_y)+{u}_{z}\times FFT(a_z)\]

On calcule également la \(FFT\) de cette réponse et le spectre d’oscillateur fourni par CALC_FONCTION/SPEC_OSCI. On fait de même pour \({\dot{u}}_{M}\) et \({\ddot{u}}_{M}\). Toutes ces fonctions sont stockées dans une table. On retrouve ainsi pour chaque cas de chargement (pour le premier NUME_CAS = 0):

  • Sur la première ligne, les «fonctions chargement», c’est-à-dire les accélérogrammes de l’excitation (temporelle, NOM_PARA=”INST”) dans les 3 directions: FONC_X, FONC_Y, FONC_Z

  • Sur la deuxième ligne, les \(FFT\) de ces signaux (NOM_PARA=”FREQ”)

  • Puis pour chaque point, l’évolution du déplacement, vitesse et accélération. Avec pour chacun, le signal, sa \(FFT\) et le spectre d’oscillateur

6.3.1. Opérandes EXCIT_FORC, ACCE_X, ACCE_Y, ACCE_Z#

L’option EXCIT_FORC est utilisée pour imposer le chargement via les forces sismiques équivalentes calculées au préalable avec Miss3D. Dans ce cas on fournira un accélérogramme dans une ou plusieurs directions (ACCE_X, ACCE_Y, ACCE_Z).

6.3.2. Opérande NORME, AMOR_SPEC_OSCI, LIST_FREQ_SPEC_OSCI#

Ces paramètres sont transmis à CALC_FONCTION pour l’option SPEC_OSCI (u4.32.04) où AMOR_REDUIT a été renommé en AMOR_SPEC_OSCI pour ne pas confondre avec la liste d’amortissements utilisés pour le calcul harmonique. De même LIST_FREQ a aussi été renommé ici en LIST_FREQ_SPEC_OSCI pour éviter les confusions avec le mot-clé LIST_FREQ .

6.4. Post-traitement des résultats aux points de contrôle#

On est ainsi dans le cas TYPE_RESU=”TABLE_CONTROL”.

6.4.1. Opérande GROUP_MA_CONTROL#

Il s’agit du groupe des mailles ponctuelles localisant les points de contrôle (transmis à IMPR_MACR_ELEM). Lors du post-traitement, des fonctions réponses sont créées pour chacun des points qui sont pris dans l’ordre de définition de ce groupe de mailles. Ainsi, dans la table, le point désigné PC1 ne correspond pas de manière générale à un nœud ou groupe de nœud nommé PC1. Il s’agit de la première maille ponctuelle de GROUP_MA_CONTROL.

6.4.2. Opérande TOUT_CHAM#

Si cet opérande est absent, on ne post-traite en temps que les accélérations aux points de contrôle. S’il est présent avec la valeur TOUT_CHAM=”OUI”, on va également post-traiter en temps les champs de vitesse et déplacement.

6.4.3. Opérandes ACCE_X, ACCE_Y, ACCE_Z, NORME, AMOR_SPEC_OSCI, LIST_FREQ_SPEC_OSCI#

Identiques aux paragraphes de l’option TYPE_RESU=”TABLE”.

6.4.4. Table produite#

Le chargement appliqué dans le calcul Miss3D est une accélération harmonique unitaire. Les deux premières lignes correspondent aux accélérations ACCE_X/Y/Z fournies par l’utilisateur, interpolées avec le pas de temps fourni, et sa \(FFT\):

  • En chaque point de contrôle, on récupère la fonction de transfert dans les trois directions à cette sollicitation. Il s’agit des lignes avec TRANSFERT/FREQ

  • Ensuite, on a la combinaison \({a}_{Mx}(f)={ft}_{x}(f)\times FFT({a}_{x})\) et même chose en y et z en fonction du chargement appliqué

  • On calcule également la \(FFT\) de cette réponse et le spectre d’oscillateur fourni par CALC_FONCTION/SPEC_OSCI

  • Toutes ces fonctions sont stockées dans la table produite

Avertissement

Le paramètre de la table désignant le point de contrôle est nommé GROUP_NO pour être homogène au cas TABLE. Comme on l’a vu plus haut, il s’agit simplement d’un numéro de point dans le groupe de mailles des points de contrôle.

6.5. La définition de la liste des fréquences de calcul - l’opérande LIST_FREQ_CALC#

La définition de la liste des fréquences pour le calcul fréquentiel est automatiquement gérée dans le cas où TYPE_RESU vaut “TRAN_GENE”. Dans le cas où la charge est fourni en temporel (TYPE_EXCIT=”INST”) on aurait:

  • Un pas de fréquence \(\Delta f=\frac{1}{n_{pas}\cdot\Delta t}\)

  • Une fréquence de calcul maximale \(f_{max}=\frac{1}{2\Delta t}\)

où:

  • \(n_{pas}\) est le nombre de pas de temps de la liste des instants

  • \(\Delta t\) est le pas de temps

Dans le cas où la charge est fourni en fréquentiel (TYPE_EXCIT=”FREQ”), la plage de fréquences utilisées est celle du signal harmonique. Dans le cas où TYPE_RESU vaut “HARM_GENE”,”TABLE” ou “TABLE_CONTROL” on peut définir la liste des fréquences via le mot-clé LIST_FREQ_CALC.

Avertissement

Le mot-clé LIST_FREQ_CALC n’est pas à renseigner si TYPE_RESU vaut “TRAN_GENE”.

6.5.1. Les paramètres AUTO, FREQ_MIN, FREQ_MAX, FREQ_PAS, FREQ et LIST_FREQ#

Si AUTO vaut “OUI”, la liste de fréquences est automatiquement gérée comme dans le cadre où TYPE_RESU vaut “TRAAN_GENE”. Si AUTO vaut “NON”, l’utilisateur peut définir la liste de fréquences:

  • Soit il utilise FREQ_MIN, FREQ_MAX et FREQ_PAS pour renseigner respectivement la fréquence minimale, la fréquence maximale et le pas de fréquence

  • Soit il utilise FREQ pour renseigner une liste de fréquences à la main

  • Soit il utilise LIST_FREQ pour renseigner une liste de fréquence construite avec DEFI_LIST_REEL