6. Modélisation A

6.1. 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.

6.2. 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).

6.3. Paramètres du calcul

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

6.4. Grandeurs testées et résultats

Identification	Transitoire	Harmonique	Différence
\(\mathrm{P3}\) – \(\mathrm{DX}\) (2.33 s) iter=1	-4.58502E-2	-4.58588E-2	0.019 %
\(\mathrm{P3}\) – \(\mathrm{DY}\) (2.33 s) iter=1	5.67299E-3	5.67541E-3	0.043 %
\(\mathrm{P3}\) – \(\mathrm{DX}\) (2.34 s) iter=2	-4.82202E-2	-4.82436E-2	0.048 %
\(\mathrm{P3}\) – \(\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) :

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