Modélisation D
==============

Géométrie
---------

On peut créer la géométrie en définissant les points :math:`A,C,D,B` , puis les deux droites :math:`\mathrm{AC}` et :math:`\mathrm{DB}` , et l’arc de cercle :math:`\mathrm{CD}` . Ceci peut se faire de plusieurs façons:

•    En définissant les coordonnées de chaque point, y compris le point milieu de l’arc :math:`\mathit{CD}` (Menu New Entity → Basic → …);

•    De façon plus simple en utilisant l'outil «Sketcher» (Menu New Entity → Basic→2D Sketch).

Après avoir construit l’ensemble de la ligne, il faut définir les groupessous ce «Wire»: il suffira pour ces travaux pratiques de créer les groupes :math:`A` , :math:`B` , :math:`\mathrm{AC}` , :math:`\mathrm{CD}` , et :math:`\mathrm{DB}` ( Menu New Entity → Group → Create Group ).


Maillage
--------

Lancer le module Mesh de la plate-forme Salome-Meca.

Le maillage est défini par le menu Mesh → Create Mesh. Sélectionner la géométrie à mailler, puis l’algorithme et l'hypothèse de discrétisation par dimension: 1D Wire Discretisation avec l’hypothèse de 15 segments par edge.


Création et lancement du cas de calcul (via AsterStudy)
-------------------------------------------------------

Lancer le module AsterStudyde la plate-forme Salome-Meca.

Puis en colonne gauche, cliquer surl'onglet Case View.

On définit le fichier de commandes du cas de calcul.

Nota:ajouter des commandes par le menu Commands→ Show All.

L’étude nécessite un premier calcul thermique transitoire suivi d’un calcul mécanique.

Les principales étapes pour la création et le lancement du cas de calcul sont les suivantes:


*   Lire le maillage au format MED: Commande LIRE_MAILLAGE.


*   Définir les éléments finis utilisés: Commande AFFE_MODELE.


*   Définir le matériau: Commande DEFI_MATERIAU.


*   Affecter le matériau: Commande AFFE_MATERIAU.


*   Affecter les caractéristiques des éléments poutres: Commande AFFE_CARA_ELEM.


*   Affecter les conditions aux limites mécaniques et le chargement : Commande AFFE_CHAR_MECA/ Enforce DOF (DDL_IMPO). Le point :math:`A` est encastré.


*   Définirles matrices du problème élastique : CommandeASSEMBLAGE / MASS_MECA et RIGI_MECA.


*   Calculer les 5 premiers modes propres : Commande CALC_MODES.


*   Imprimer les modes propres au format MED: Commande IMPR_RESU. On imprimera le maillage et les modes.

Pour lancerle cas de calcul, en colonne gauche, cliquer sur l'onglet History View.

*Analyse transitoire:*


*   Construction de la force ponctuelle


1.
    *
        *
            *   Affecter la charge «force au point :math:`B` » : Commande AFFE_CHAR_MECA/ FORCE_NODALE.


            *   Calculer les vecteurs forces élémentaires : Commande CALC_VECT_ELEM.


            *   Construire un champ aux nœuds par assemblage du vecteur force : Commande ASSE_VECTEUR.


            *   Définir la fonction évolution du temps : CommandeFORMULE.


*   Transitoire sur base modale


1.  Projection du problème assemblé sur la base des modes propres : Commande PROJ_BASE.


2.  Calculer le transitoire par recombinaison modale : Commande DYNA_VIBRA.


3.  Extraire les déplacements en :math:`Y` en :math:`B` : Commande RECU_FONCTION.


4.  Imprimer ces fonctions au format TABLEAU: Commande IMPR_FONCTION.