Modélisation A
==============


Caractéristiques de la modélisation
-----------------------------------

Dans cette modélisation, la résolution du problème dynamique s’opère dans une boucle de calculs linéaires où on recalcule à chaque fois sur toute la plage temporelle le complément de forces nodales dû à la non linéarité du décollement. On utilise, soit un calcul harmonique avec une évolution fréquentielle et retour en temps par transformée de Fourier dans la méthode temps-fréquence classique, soit un calcul transitoire sur base physique avec une méthode strictement temporelle alternative.


Caractéristiques du maillage
----------------------------

Le modèle est composé de 55 nœuds (285 ddls), 42 éléments (32 éléments plaques DKT et 10éléments discrets DIS_T).


Paramètres du calcul
--------------------

Chaque calcul dynamique transitoire est mené sur un intervalle de :math:`5s` par pas de temps de :math:`0.005s` archivé tous les 2 pas. Chaque calcul harmonique est mené avec un pas de :math:`1/20.48\mathrm{Hz}` qui permet de restituer une fenêtre temporelle de :math:`20.48s` suffisante pour bien calculer la FFT de la force de pesanteur constante dans le temps; la fréquence maximale de calcul vaut :math:`25\mathrm{Hz}` et celle prolongée est de :math:`50\mathrm{Hz}` afin d’obtenir un pas de temps de :math:`0.01s` dans la fenêtre temporelle restituée par FFT.


Grandeurs testées et résultats
------------------------------


.. csv-table::

    "**Identification**", "**Transitoire**", "**Harmonique**", "**Différence**"
    ":math:`\mathrm{P3}` – :math:`\mathrm{DX}` (2.33 s) iter=1", "-4.58502E-2", "-4.58588E-2", "0.019 %"
    ":math:`\mathrm{P3}` – :math:`\mathrm{DY}` (2.33 s) iter=1", "5.67299E-3", "5.67541E-3", "0.043 %"
    ":math:`\mathrm{P3}` – :math:`\mathrm{DX}` (2.34 s) iter=2", "-4.82202E-2", "-4.82436E-2", "0.048 %"
    ":math:`\mathrm{P3}` – :math:`\mathrm{DY}` (2.34 s) iter=2", "6.55566E-3", "6.54736E-3", "0.127 %"


Comparaison de la FFT fortran (REST_SPEC_TEMP) et de la FFT python (CALC_FONCTION) :


.. csv-table::

    "**Identification**", "**REST_SPEC_TEMP**", "**CALC_FONCTION**", "**Écart absolu**"
    ":math:`\mathrm{P3}` – :math:`\mathrm{DX}` (2.34 s) iter=2", "-4.82202E-2", "-4.82202E-2", "2.5673907444E-16"
    ":math:`\mathrm{P3}` – :math:`\mathrm{DY}` (2.34 s) iter=2", "6.55566E-3", "6.55566E-3", "-7.4593109467E-17"