v5.02.105 SDNL105 - Sous-structuration transitoire non linéaire : choc de 3 poutres entre elles#
Résumé:
Le domaine d’application de ce test concerne la dynamique des structures, et plus particulièrement le calcul de réponse transitoire non linéaire par sous-structuration dynamique.
Il s’agit de calculer la réponse transitoire non linéaire de 3 poutres en flexion avec chocs au centre des poutres. Les poutres sont modélisées par des éléments de type POU_D_E (modèle d’Euler).
Les résultats de référence sont issus d’un calcul transitoire direct par recombinaison modale. Ce test permet donc de valider les outils de calcul de réponse transitoire par sous‑structuration, dans le cas de la prise en compte de non-linéarités de type choc entre structures mobiles.
Solution de référence#
Méthode de calcul utilisée pour la solution de référence#
La solution de référence est donnée par un calcul transitoire direct par superposition modale (modélisation A).
Résultats de référence#
Valeur des déplacements et vitesse des nœuds des 3 poutres selon la direction \(Y\) et à l’instant \(t=1s\) .
Déplacement (\(m\) ) |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
|
Schéma d’intégration d’Euler |
||
Poutre de gauche |
1.64 10–2 |
2.54 10–2 |
Poutre milieu |
1.12 10–2 |
4.49 10–2 |
Poutre de droite |
5.90 10–3 |
1.05 10–1 |
Schéma d’intégration de Devogelaere |
||
Poutre de gauche |
1.64 10–2 |
2.54 10–2 |
Poutre milieu |
1.12 10–2 |
4.41 10–2 |
Poutre de droite |
5.89 10–3 |
1.05 10–1 |
Schéma d’intégration à pas de temps adaptatif d’ordre 2 |
||
Poutre de gauche |
1.64 10–2 |
2.55 10–2 |
Poutre milieu |
1.12 10–2 |
4.41 10–2 |
Poutre de droite |
5.91 10–3 |
1.05 10–1 |
Incertitude sur la solution#
Solution numérique.
Modélisation A#
Caractéristiques de la modélisation#
La poutre est maillée en segments auxquels sont affectés des éléments de type POU_D_E.
Le problème transitoire traité, projeté sur la base des 15 premiers modes propres de la structure, est résolu directement par l’opérateur de calcul transitoire par recombinaison modale.
Caractéristiques du maillage#
Nombre de nœuds: 41
Nombre de mailles et types : 42 SEG2
Valeurs observées : référence pour modélisation B#
Identification |
Aster |
Schéma d’intégration d’Euler |
|
Poutre de gauche : Déplacement (\(m\) ) |
1.64 10–2 |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
2.54 10–2 |
Poutre de milieu : Déplacement (\(m\) ) |
1.12 10–2 |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
4.43 10–2 |
Poutre de droite : Déplacement (\(m\) ) |
5.90 10–3 |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
1.05 10–1 |
Schéma d’intégration de Devogelaere |
|
Poutre de gauche : Déplacement (\(m\) ) |
1.64 10–2 |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
2.54 10–2 |
Poutre de milieu : Déplacement (\(m\) ) |
1.12 10–2 |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
4.41 10–2 |
Poutre de droite : Déplacement (\(m\) ) |
5.89 10–3 |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
1.05 10–1 |
Schéma d’intégration à pas de temps adaptatif d’ordre 2 |
|
Poutre de gauche : Déplacement (\(m\) ) |
1.64 10–2 |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
2.55 10–2 |
Poutre de milieu : Déplacement (\(m\) ) |
1.12 10–2 |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
4.41 10–2 |
Poutre de droite : Déplacement (\(m\) ) |
5.91 10–3 |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
1.05 10–1 |
Modélisation B#
Caractéristiques de la modélisation#
La sous-structuration dynamique permet de calculer le comportement vibratoire des 3 poutres à partir des caractéristiques dynamiques d’une seule poutre. Celle-ci est maillée en segments auxquels sont affectés des éléments de type POU_D_E.
La structure est étudiée à l’aide de la méthode de sous-structuration avec interfaces de type « Craig‑Bampton » (interfaces bloquées).
La base des 15 premiers modes propres de la structure complète est calculée par sous-structuration. Puis le problème transitoire, projeté sur cette base, est résolu par l’opérateur de calcul transitoire par recombinaison modale.
Caractéristiques du maillage#
Nombre de nœuds: 15
Nombre de mailles et types : 14 SEG2
Grandeurs testées et résultats#
Identification |
Référence |
Aster |
% différence |
Schéma d’intégration d’Euler |
|||
Poutre de gauche : Déplacement (\(m\) ) |
1.64 10–2 |
1.64 10–2 |
|
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
2.54 10–2 |
2.54 10–2 |
|
Poutre de milieu : Déplacement (\(m\) ) |
1.12 10–2 |
1.12 10–2 |
< 0.01 % |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
4.43 10–2 |
4.43 10–2 |
|
Poutre de droite : Déplacement (\(m\) ) |
5.90 10–3 |
5.90 10–3 |
|
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
1.05 10–1 |
1.05 10–1 |
|
Schéma d’intégration de Devogelaere |
|||
Poutre de gauche : Déplacement (\(m\) ) |
1.64 10–2 |
1.64 10–2 |
|
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
2.54 10–2 |
2.54 10–2 |
|
Poutre de milieu : Déplacement (\(m\) ) |
1.12 10–2 |
1.12 10–2 |
< 0.01 % |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
4.41 10–2 |
4.41 10–2 |
|
Poutre de droite : Déplacement (\(m\) ) |
5.89 10–3 |
5.89 10–3 |
|
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
1.05 10–1 |
1.05 10–1 |
|
Schéma d’intégration à pas de temps adaptatif d’ordre 2 |
|||
Poutre de gauche : Déplacement (\(m\) ) |
1.64 10–2 |
1.64 10–2 |
|
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
2.55 10–2 |
2.55 10–2 |
|
Poutre de milieu : Déplacement (\(m\) ) |
1.12 10–2 |
1.12 10–2 |
< 0.01 % |
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
4.41 10–2 |
4.41 10–2 |
|
Poutre de droite : Déplacement (\(m\) ) |
5.91 10–3 |
5.91 10–3 |
|
Vitesse (\({\mathrm{m.s}}^{-1}\) ) |
1.05 10–1 |
1.05 10–1 |
Synthèse des résultats#
La précision sur les déplacements et les vitesses des nœuds milieux des 3 poutres à l’instant \(t=1s\) est excellente (erreur relative \(<\mathrm{0.01s}\) ).
Ce test valide donc les opérateurs de calcul transitoire non-linéaire par sous-structuration dynamique.