Master Mécanique Spécialité Innovations
en Conception et Matériaux (ICM) à
Polytech Orléans
TD 9
Module : Simulation numérique avancée 2 Partran/Marc
Encadré par : Pr. Benoit MAGNAIN
Réaliser par :
EL-HOR Nadia
Année Universitaire 2024/2
Table de matière
I . Travail réalisé : ............................................................................................................................ 3
II.Méthode de Newmark :............................................................................................................. 3
I. Le Schéma différences finies : .............................................................................................. 5
Conclusion .............................................................................................................................................. 5
Introduction :
L'objectif principal de ce TP est d’analyser l’impact des choix numériques sur les résultats
obtenus dans le cadre de la résolution du problème de flexion. En particulier, nous mettrons en
évidence les effets des différentes méthodes d’intégration numérique (explicite ou implicite),
de la taille du pas de temps, ainsi que du type de matrice de masse (concentrée ou répartie) sur
la précision et la stabilité de la simulation.
La résolution numérique de ce type de problème nécessite l’utilisation de méthodes de calcul
adaptées à la discrétisation du système, tant dans le domaine temporel que spatial. Ces choix
influencent directement la qualité des résultats obtenus, ainsi que la complexité des calculs. Par
ailleurs, une gestion appropriée du pas de temps et de la matrice de masse est essentielle pour
assurer à la fois la précision des résultats et la stabilité de la simulation.
Ce travail permettra de mieux comprendre l'importance de ces paramètres dans la modélisation
numérique des structures et de tirer des conclusions sur les meilleures pratiques en fonction des
différents scénarios rencontrés.
I . Travail réalisé :
La poutre est modélisée en 2D avec des dimensions de L = 1 m et h = 0.03 m. Elle est encastrée à son
extrémité gauche et soumise à une excitation verticale de F = 100 N à son extrémité droite. Le
matériau utilisé est un acier standard, caractérisé par un module de Young E = 210 GPa et un
coefficient de Poisson ν = 0.33.
Figure 1 : la poutre en flexion console
II.Méthode de Newmark :
La méthode de Newmark est une technique d'intégration numérique largement utilisée pour la
résolution des équations différentielles en dynamique des structures. Elle permet de résoudre
l'équation du mouvement sous forme discrétisée dans le temps et peut être formulée sous
différentes variantes selon les paramètres choisis.
Avant d'aborder les aspects dynamiques et la résolution numérique du problème, il est
essentiel de comprendre le comportement de la poutre dans le cas statique. Dans cette
première étape, nous avons visualisé la poutre encastrée à une extrémité, soumise à une
charge verticale en son extrémité libre.
Figure 2: Visualisation de cas élastique
Dans notre cas, la méthode de Newmark est utilisée pour simuler la réponse dynamique de la
poutre sous excitation verticale. Selon le choix des paramètres β et γ, nous pouvons comparer
les résultats obtenus entre les schémas implicites et explicites en termes de stabilité et de
précision.
Figure 3 : Configuration du job
Le schéma de Newmark est une méthode implicite utilisée pour résoudre les équations
différentielles du déplacement. Il est inconditionnellement stable lorsque 0.52β≥ γ ≥0.5. En
faisant varier le pas de temps, il est possible d'observer l'évolution de la flèche à l'extrémité
droite en fonction du temps pour différentes valeurs de pas de temps.
Figure 4 : la flèche en fonction des différents pas de temps
L'analyse des courbes montre que plus le pas de temps Δt est petit, plus les résultats gagnent en
précision, permettant ainsi de mieux capturer les vibrations de la poutre. Toutefois, une
réduction du pas de temps s'accompagne d'une augmentation du temps de calcul, ce qui peut
devenir contraignant pour des simulations longues.
Dans notre cas, un pas de temps de Δt=0.001 s constitue un bon compromis entre précision et
performance computationnelle.
Par ailleurs, la diagonalisation de la matrice de masse n’a pas d’impact significatif sur le temps
de calcul. Les résultats obtenus avec une matrice de masse diagonalisée sont similaires à ceux
obtenus avec la matrice complète, ce qui suggère que cette simplification peut être utilisée sans
perte de précision notable.
[Link] Schéma différences finies :
Ce schéma, étant conditionnellement stable, impose l'utilisation d'un pas de temps très réduit,
de l'ordre de 10−7s, en raison de la taille des éléments adoptée. Cette contrainte entraîne un
temps de calcul extrêmement long, rendant la simulation difficilement exploitable.
En conséquence, il n’a pas été possible d’obtenir des résultats comparables avec d’autres
méthodes, car la stabilité du schéma explicite dépend fortement du choix du pas de temps, qui
doit rester suffisamment petit pour éviter toute divergence des calculs.
Conclusion
Ce rapport compare les schémas implicite (Newmark) et explicite (différences finies) pour
l’analyse dynamique d’une poutre. Le schéma de Newmark, grâce à sa stabilité inconditionnelle
sous la condition 2β≥γ≥0.52, offre un bon compromis entre précision et temps de calcul avec
un pas de temps optimal de Δt=0.001s.
À l’inverse, le schéma explicite, bien que théoriquement utilisable, nécessite un pas de temps
extrêmement faible (10−7 s), ce qui entraîne un coût de calcul prohibitif et empêche son
application pratique dans notre étude.