Post-traitement =============== Analyse du comportement global de la structure ---------------------------------------------- Déplacements et efforts maximaux ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ L’analyse de l’évolution des déplacements et des réactions aux appuis (efforts tranchants, efforts normaux et moments fléchissants globaux) fournit des éléments importants pour appréhender le comportement global d’une structure. **Remarques** *Il est pratique d’utiliser le mot-clé* *OBSERVATIONafin de sauvegarder l’ensemble des résultats pour un point donné et non pour toute la structure.* .. _Ref181537223: Détermination du spectre d'oscillateur ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ La détermination du spectre d’oscillateur d'un accélérogramme :math:`a(t)` est disponible par l’opérateur CALC_FONCTION[U6.62.04] avec le mot clé SPEC_OSCI: il est obtenu par intégration numérique de l’équation de Duhamel par la méthode de Nigam [:ref:`R5.05.01 `]. Cette commande fournit le spectre de pseudo accélération absolue et, sur demande, le spectre de pseudo vitesse ou le spectre de déplacement relatif. D'une manière générale, l'utilisateur doit toujours vérifier que la plage fréquentielle employé pour le calcul du spectre d'oscillateur est cohérente avec le contenu fréquentiel du signal d'entrée. Cela est d'autant plus crucial si l'on cherche à obtenir, par exemple, la valeur asymptotique de la pseudo‑accélération. Dans certains cas, la liste de fréquence par défaut du mot clé SPEC_OSCI peut être trop restrictive. .. _Ref181537306: Fréquences propres ~~~~~~~~~~~~~~~~~~ L’analyse de l’évolution des fréquences propres d’une structure en cours de chargement permet d’avoir une bonne indication sur l’état d’endommagement de celle-ci. Les phénomènes de fissuration font diminuer la raideur globale de la structure et on observe une chute de fréquence. Si on détermine les fréquences propres et l’amortissement de la structure endommagée, on peut ensuite réaliser un calcul de réponse sismique spectral sur celle-ci. La détermination des fréquences propres de la structure endommagée est possible à l’aide du mot-clé MODE_VIBR de l’opérateur DYNA_NON_LINE [:external:ref:`U4.53.01 `]. Il s’agit d’une méthode permettant de réaliser une analyse modale vibratoire (fréquences propres et modes associés) sur les matrices de raideur et de masse globales de DYNA_NON_LINE. Comme on doit calculer les fréquences propres de la structure pour une matrice de raideur élastique endommagée, il est nécessaire d’avoir remis la structure au repos avant de réaliser le calcul des fréquences propres par l’opérateur MODE_VIBR. On se reportera au cas-test SDNV106 [:external:ref:`V5.03.106 `] pour plus de détails. Une autre possibilité consiste à réaliser, après le chargement sismique et après la mise au repos de la structure, un chargement de type «bruit blanc» à faible niveau de sollicitations à l’aide de *Code_Aster* . On réalise alors une analyse de Fourier de la réponse obtenue afin de déterminer les fréquences propres de la structure endommagée. Bilan d’énergies dissipées ~~~~~~~~~~~~~~~~~~~~~~~~~~ Le calcul des énergies sur tout ou partie de la structure se fait par l’opérateur POST_ELEM :ref:`u4.81.22`: l’énergie cinétique, l’énergie de déformation élastique, l’énergie de déformation totale, le travail des efforts extérieurs, l’énergie de dissipation pour les éléments de plaque de la famille DKTG avec les lois de comportement global de béton armé GLRC_DM et DHRC. Cependant, l’énergie d’amortissement visqueux n’est actuellement pas calculée, ni l’énergie de déformation pour les lois adoucissantes telles que ENDO_ISOT_BETON. De plus, le calcul du travail des efforts extérieurs s’effectue à partir d’une partie des forces nodales et ne tient pas compte des forces d’inertie ou d’amortissement. Le bilan énergétique est donc à réaliser avec précaution. On peut aussi avec le mot-clé ENERGIEde l’opérateur DYNA_NON_LINEfaire le calcul du bilan d'énergie sur toute la structure, y compris l’énergie d’amortissement visqueux. Il est stockédans la table de nom PARA_CALC. Il est extrait avec la commande RECU_TABLE [:external:ref:`U4.71.02 `]. Analyse du comportement local de la structure --------------------------------------------- Introduction ~~~~~~~~~~~~ Les options de calcul de champs citées dans ce chapitre sont définies dans CALC_CHAMP [:external:ref:`U4.81.04 `]. **Béton** L’analyse des cartes d’isovaleurs d’endommagement permet d’appréhender le comportement local d’une structure. On peut déterminer les niveaux maximaux de contraintes et de fissuration atteints en certains points de la structure. **Remarque** *Si on observe que seuls quelques éléments finis sont très fortement endommagés, il peut s’agir d’un phénomène de localisation de l’endommagement dû au modèle de comportement adoucissant utilisé. Ce phénomène entraîne une dépendance de la réponse au maillage. Afin de déterminer s’il y localisation ou non, on préconise de reprendre le calcul en remaillant plus finement la zone suspectée. Si on observe que l’endommagement a de nouveau lieu sur une seule bande de maille (réduction de la zone endommagée), c’est qu’on a localisation.* *Il est à noter que les modèles de délocalisation des lois de comportement par régularisation de la déformation ne sont pas encore actuellement disponibles en dynamique non linéaire dans Code_Aster (cf. §* :ref:`4.4.4 ` *).* **Armatures** L’analyse des cartes d’isovaleurs des déformations plastiques cumulées dans les armatures permet de déterminer le niveau de plastification atteint dans les aciers. Modélisation locale (éléments massifs) ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Il n’y a pas de spécificité liée au post-traitement dans ce cas. Il n’est donc pas détaillé ici. **Béton** On calcule directement les contraintes aux nœuds par élément SIGM_ELNO. **Armatures** De même, les éléments de grilles membranes (GRILLE_MEMBRANE) ou les éléments de barre (BARRE) se post-traitent en calculant les contraintes associées SIGM_ELNO. **Remarque** *En HPP, on calcule directement les déformations aux nœuds par éléments* *EPSI_ELNOà partir des déplacements.* .. _Ref181537832: Éléments de plaques et coques multicouches ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ **Béton** Pour les coques multicouches, on post-traite les contraintes par éléments aux nœuds (SIGM_ELNO) à partir des contraintes aux points d’intégration de chaque couche (SIEF_ELGA) calculées lors d’un calcul non linéaire. Si l'on souhaite obtenir ces contraintes dans une couche et une position particulière, il faut utiliser l'opérateur POST_CHAMP. Ces contraintes sont calculées dans le repère local de la coque défini par l’utilisateur dans la commande AFFE_CARA_ELEM. On peut post-traiter de la même manière les variables internes (VARI_ELNO). **Armatures** Les éléments de grilles excentrées (GRILLE_EXCENTRE) ou les éléments de barre (BARRE), se post-traitent simplement en calculant les contraintes associées SIGM_ELNO. .. _Ref181537718: Éléments de poutres multifibres ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Pour les poutres multifibres, on ne dispose pas actuellement d’option de calcul des contraintes dans les fibres directement dans CALC_CHAMP. On peut toutefois récupérer les contraintes dans les fibres en utilisant les opérandes de localisation d’un champ dans l’opérateur RECU_FONCTION: .. code-block:: text / NOEUD = no, [noeud] / GROUP_NO = grno, [gr_noeud] / MAILLE = ma, [maille] / GROUP_MA = grma, [gr_maille] / POINT = nupoint,[I] / SOUS_POINT = nusp, [I] L'entier nupointprécise le numéro du point de Gauss dont on veut récupérer la valeur (cas des cham_elem "aux points de GAUSS"). L'entier nuspprécise le numéro du sous point, c’est-à-dire de la fibre, sur laquelle on vient relever les contraintes (cf. le maillage de la section § :ref:`3.4 ` ). On peut ainsi analyser les contraintes et les déformations dans les fibres béton et dans les fibres acier. Éléments de plaques et coques « globales » ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Dans cette modélisation, on ne dispose pas directement des contraintes et des déformations dans l’épaisseur de la coque. Il est toutefois possible de reconstruire le champ de déformation (mais pas le champ de contrainte) pour une cote donnée à partir des champs d’efforts et de déformations généralisés. Comme le champ de déformation est considéré linéaire dans la section, on a: :math:`\lbrace \begin{array}{c}{\varepsilon}_{x}(z)={\varepsilon}_{x}^{\mathrm{tot}}+z{\kappa}_{x}^{\mathrm{tot}}\\ {\varepsilon}_{y}(z)={\varepsilon}_{y}^{\mathrm{tot}}+z{\kappa}_{y}^{\mathrm{tot}}\end{array}` avec :math:`{\varepsilon}_{x}^{\mathrm{tot}}` , :math:`{\varepsilon}_{y}^{\mathrm{tot}}` , :math:`{\kappa}_{x}^{\mathrm{tot}}` et :math:`{\kappa}_{y}^{\mathrm{tot}}` les déformations généralisées et :math:`{\varepsilon}_{x}(z)` et :math:`{\varepsilon}_{y}(z)` les déformations dans l’épaisseur de la coque globale. Éléments discrets porteurs de lois de comportement linéaire ou non ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Dans cette modélisation, on dispose directement des déplacements (voire rotations) des nœuds extrémité des éléments discrets, mais aussi des efforts qui y transitent. De plus, en cas de loi de comportement non linéaire, *cf.* [:ref:`R5.03.17 `], on a accès aux valeurs de diverses variables internes utiles pour analyser le fonctionnement détaillé (contribution anélastique, énergie dissipée, etc).