Préparation à l’Agrégation de Mécanique.
Mécanique des Solides Indéformables Préparation à l’Agrégation de Mécanique
S. Pommier
Contacts
Sylvie Pommier Zone 1A2002 Tél. : 01 47 40 28 69 Mèl :
[email protected] :
[email protected] han.fr
Vincent Wendling Zone 1A2004 Tél. : 01 47 40 21 87 Mèl :
[email protected]
École Normale Supérieure de Cachan, 61, Avenue du président Wilson, Wilson, 94235 CACHAN
Préparation à l’Agrégation de Mécanique.
INTRODUCTION
1
Introduction_____________________________________________________________________2 ◊
Plan du cours___________________ cours__________________________ _______________ _______________ ______________ _______________ _____________ _____ 2
◊
Exemple Exemple 1 : Vibration Vibration d’une aile d’avion d’avion ________ _______________ ______________ ______________ ______________ __________ ___ 3
◊
Exemple Exemple 2 : Ecoulemen Ecoulementt d’un milieu milieu granulaire_____________ granulaire____________________ _______________ ______________ ______ 4
◊
Sources Sources bibliograp bibliographiques hiques ________ _______________ ______________ _______________ _______________ _______________ ______________ ______ 4
RAPPELS DE MATHEMATIQUES
1
Calcul Calcul Vectoriel Vectoriel _______ ______________ _______________ _______________ _______________ _______________ ______________ _______________ ______________6 ______6 1.1
Opérations Opérations sur les vecteurs vecteurs _______ _______________ _______________ ______________ _______________ _______________ ______________ _________ 6 ◊
Produit Produit scalaire scalaire _______ _______________ _______________ ______________ _______________ _______________ ______________ ______________ _________ 6
◊
Produit Produit vectoriel vectoriel ________ _______________ ______________ _______________ _______________ ______________ _______________ ______________ ______ 6
◊
Produit Produit mixte mixte _______ _______________ ________________ _______________ _______________ _______________ _______________ _______________ _________ 7
◊
Division Division vectoriell vectoriellee _______ ______________ ______________ _______________ _______________ ______________ _______________ _____________ _____ 7
1.2
Champs Champs de vecteurs______ vecteurs_____________ ______________ ______________ _______________ _______________ ______________ _______________ __________ 8 ◊
Glisseur Glisseur _______ ______________ _______________ _______________ _______________ _______________ _______________ _______________ _____________ ______ 8
◊
Moment Moment en un point d’un glisseur________ glisseur________________ _______________ ______________ _______________ _______________ _______ 8
◊
Moment Moment d’un glisseur glisseur par rapport rapport à un axe_________________ axe_________________________ _______________ _____________ ______ 9
◊
Ensembles Ensembles de glisseurs glisseurs _______ _______________ _______________ _______________ _______________ ______________ _______________ __________ 9
1.3
Torseurs_____ Torseurs____________ _______________ _______________ ______________ _______________ _______________ _______________ _______________ _________ __ 10 ◊
Définition Définition ________ _______________ _______________ _______________ _______________ ________________ _______________ _______________ __________ 10
◊
Torseur Torseur associé associé à un ensemble ensemble de glisseur glisseur _______ ______________ ______________ _______________ _______________ _________ 10
◊
Invariants Invariants du torseur._________ torseur._________________ _______________ ______________ _______________ _______________ _______________ ________ 10
◊
Point central, central, axe central, central, moment moment central central d’un torseur torseur ______ ______________ _______________ _____________ ______ 11
◊
Symétrie Symétrie du champ champ des moments moments d’un torseur. torseur. Origine Origine du mot « torseur torseur ». ______ ___________ _____ 11
1.4
Opérations Opérations sur les torseurs torseurs ________ _______________ ______________ _______________ _______________ ______________ ______________ _______ 12 ◊
Addition_____ Addition_____________ _______________ _______________ _______________ _______________ _______________ _______________ ______________ ______ 12
◊
Multiplicati Multiplication on par un réel ________ _______________ _______________ _______________ ______________ _______________ _____________ _____ 12
◊
Décompositi Décomposition on ________ _______________ _______________ _______________ _______________ _______________ _______________ _____________ _____ 12
◊
Produit Produit ou co-moment co-moment de deux torseurs torseurs _______ _______________ _______________ ______________ ______________ __________ ___ 13
◊
Torseur Torseur à structure_____ structure____________ _______________ _______________ _______________ _______________ _______________ ______________ ______ 13
1.5
Champ Champ equiprojecti equiprojectiff de vecteurs___________ vecteurs__________________ ______________ ______________ _______________ _______________ _______ 14 ◊
Définition Définition ________ _______________ _______________ _______________ _______________ ________________ _______________ _______________ __________ 14
École Normale Supérieure de Cachan, 61, Avenue du président Wilson, Wilson, 94235 CACHAN
Préparation à l’Agrégation de Mécanique. ◊
Propriétés Propriétés________ _______________ _______________ _______________ _______________ _______________ _______________ _______________ __________ ___ 14
◊
Champ Champ des moments moments d’un torseur. torseur. _______ _______________ _______________ ______________ ______________ ______________ _______ 14
2
Dérivation vectorielle ______________________________ ____________________________________________________________15 ______________________________15 2.1
Dérivée Dérivée d’un vecteur vecteur ________ _______________ ______________ _______________ _______________ ______________ _______________ ____________ ____ 15 ◊
Définition Définition ________ _______________ _______________ _______________ _______________ ________________ _______________ _______________ __________ 15
◊
Propriétés Propriétés________ _______________ _______________ _______________ _______________ _______________ _______________ _______________ __________ ___ 15
2.2
Changemen Changementt de base de dérivation dérivation _______ ______________ ______________ _______________ _______________ ______________ _________ __ 15 ◊
Vocabulaire Vocabulaire _______ _______________ _______________ _______________ _______________ _______________ _______________ ______________ _________ __ 15
◊
Dérivée Dérivée d’un vecteur vecteur exprimé exprimé dans la base de dérivation dérivation_______ ______________ ______________ ____________ _____ 16
◊
Dérivée Dérivée d’un vecteur vecteur exprimé exprimé dans une base distincte distincte de la base de dérivation. dérivation. _______ _________ 16
◊
Propriété Propriétéss du vecteur vecteur vitesse vitesse de rotation. rotation. _______ _______________ _______________ ______________ _______________ __________ 17
3
A retenir retenir _______ _______________ ________________ _______________ _______________ _______________ _______________ _______________ _______________ __________18 __18 ◊
Champ Champ équiprojectif équiprojectif : _______ _______________ _______________ ______________ _______________ _______________ ______________ _________ __ 18
◊
Torseur : {T }=
◊
Produit Produit de deux torseurs torseurs définis définis au même même point A : ________ _______________ ______________ ______________ _________ 18
◊
Invariants Invariants du torseur torseur _______ ______________ _______________ _______________ ______________ _______________ _______________ __________ ___ 18
◊
Axe et point centraux centraux _______ ______________ _______________ _______________ ______________ _______________ _______________ _________ __ 18
◊
Changement Changement de base de dérivation dérivation _______ _______________ _______________ ______________ ______________ ______________ _______ 18
R __________________________________________________ _____________________ 18 _____________________________ M A A
CINEMATIQUE CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
1
Cinematique du point _______________________________ ____________________________________________________________2 _____________________________200 ◊
Vecteur Vecteur position position d’un point_________________ point________________________ ______________ _______________ _______________ __________ ___ 20
◊
Vecteur Vecteur vitesse vitesse d’un point_________________ point_________________________ _______________ _______________ _______________ __________ ___ 20
◊
Vecteur Vecteur accélérat accélération ion d’un point ________ _______________ _______________ _______________ ______________ _______________ ________ 20
2
Le solide indeformable _____________________________ ___________________________________________________________20 ______________________________20 ◊
Définition Définition ________ _______________ _______________ _______________ _______________ ________________ _______________ _______________ __________ 20
◊
Equivalence Equivalence Repère-So Repère-Solide lide _______ _______________ _______________ ______________ _______________ _______________ ____________ _____ 21
3
Paramétrage Paramétrage de la position relative de deux solides __________________________________ _____________________________________21 ___21 3.1
Définition Définition des coordonnée coordonnéess de l’origine l’origine d’un repère. repère. ________ _______________ ______________ ______________ _________ __ 21 ◊
Coordonnée Coordonnéess cartésienn cartésiennes es _______ ______________ _______________ _______________ ______________ _______________ ______________ ______ 21
◊
Coordonnée Coordonnéess cylindriqu cylindriques es ________ _______________ ______________ _______________ _______________ _______________ _____________ _____ 22
◊
Coordonnée Coordonnéess sphériques sphériques _______ _______________ _______________ _______________ _______________ ______________ ______________ _______ 22
3.2
4
Définition Définition de l’orientat l’orientation ion relative relative de deux bases__________________ bases_________________________ ______________ _________ __ 23 ◊
Angles d’Euler d’Euler _______ _______________ _______________ ______________ _______________ _______________ ______________ ______________ _______ 23
◊
Angles de Cardan Cardan _______ ______________ ______________ _______________ _______________ ______________ _______________ _____________ _____ 24
◊
Quaternions Quaternions _______ _______________ _______________ _______________ _______________ ______________ _______________ _______________ _________ __ 25 Cinematique du solide______________________________________________ solide_______________ _____________________________________________28 ______________28 École Normale Supérieure de Cachan, 61, Avenue du président Wilson, Wilson, 94235 CACHAN
Préparation à l’Agrégation de Mécanique.
4.1
Introduction, notations ______________________________________________________ 28
4.2
Champ des vecteurs vitesse des points d’un solide : torseur cinématique________________ 28 ◊
Propriété : Equiprojectivité. ________________________________________________ 29
◊
Calcul du vecteur rotation instantanée. ________________________________________ 29
◊
Exemple : Mouvement de translation _________________________________________ 30
◊
Exemple : Mouvement de rotation instantanée __________________________________ 30
4.3
Champ des vecteurs accélération des points d’un solide. ____________________________ 30
5
Composition des mouvements _____________________________________________________31 5.1
Introduction_______________________________________________________________ 31
5.2
Composition des vecteurs vitesse ______________________________________________ 32 ◊
Relation de composition des vecteurs vitesses __________________________________ 32
◊
Définitions______________________________________________________________ 32
◊
Généralisation ___________________________________________________________ 33
◊
Composition des torseurs cinématiques _______________________________________ 33
◊
Exemple________________________________________________________________ 33
5.3
Composition des vecteurs accélération __________________________________________ 34 ◊
Relation de composition des vecteurs accélération _______________________________ 34
◊
Définitions______________________________________________________________ 35
6
A retenir ______________________________________________________________________35 ◊
Cinématique du point : ____________________________________________________ 35
◊
Cinématique du solide :____________________________________________________ 35
◊
Formules de changement de point____________________________________________ 36
◊
Formules de composition des mouvements. ____________________________________ 36
CINEMATIQUE DES SYSTEMES DE SOLIDES
1
Système de solides ______________________________________________________________38 ◊
Définitions______________________________________________________________ 38
◊
Types de liaisons_________________________________________________________ 38
◊
Représentation d’une liaison. _______________________________________________ 39
◊
Degrés de liberté d’une liaison.______________________________________________ 39
2
Tableau des liaisons normalisées ___________________________________________________39
3
Cinématique du contact entre deux solides ____________________________________________41
4
Modélisation cinématique _________________________________________________________42 4.1
Introduction_______________________________________________________________ 42
4.2
Graphe cinématique_________________________________________________________ 42
4.3
Graphe de structure ou graphe des liaisons _______________________________________ 42
4.4
Mobilité d’un système_______________________________________________________ 42 ◊
Fermeture géométrique ____________________________________________________ 43
École Normale Supérieure de Cachan, 61, Avenue du président Wilson, 94235 CACHAN
Préparation à l’Agrégation de Mécanique. ◊
Fermeture cinématique ____________________________________________________ 43
◊
Calcul de la mobilité ______________________________________________________ 43
4.5
Exemple : Presse de modélisme _______________________________________________ 44 ◊
Plan du mécanisme _______________________________________________________ 44
◊
Construction du graphe de structure __________________________________________ 45
◊
Réduction du graphe de structure ____________________________________________ 45
◊
Construction du graphe cinématique__________________________________________ 46
◊
Mobilité du système ______________________________________________________ 46
5
A retenir ______________________________________________________________________47
CONSERVATION DE LA MASSE ET CINETIQUE
1
Système matériel à masse conservative_______________________________________________49 ◊
Définitions______________________________________________________________ 49
◊
Conséquences ___________________________________________________________ 49
2
Torseur cinétique, torseur dynamique et energie cinétique________________________________50 ◊
Torseur cinétique_________________________________________________________ 50
◊
Torseur dynamique _______________________________________________________ 50
◊
Energie cinétique_________________________________________________________ 50
◊
Autres cas. Exemple : Action de la pesanteur ___________________________________ 51
3
Centre d’inertie, opérateur d’inertie _________________________________________________51 3.1
Centre d’inertie G __________________________________________________________ 51
3.2
Opérateur d’inertie JA(Σ)_____________________________________________________ 51 ◊
Définition ______________________________________________________________ 51
◊
Relation entre JA(Σ) et JG(Σ) ou théorème de Huyghens généralisé . _________________ 51
◊
Expression dans la base (O,x,y,z) ____________________________________________ 52
3.3
Exemples_________________________________________________________________ 53 ◊
Opérateur d’inertie en O d’un disque D de rayon R, de centre O et de masse m. ________ 53
◊
Opérateur d’inertie en O d’un cylindre C de rayon R, de hauteur h et de centre O. ______ 53
◊
Opérateur d’inertie en O d’un cône de révolution C de rayon R et de hauteur h. ________ 54
◊
Opérateur d’inertie en O d’une sphère creuse S de centre O, de rayon R et de masse m.__ 54
◊
Opérateur d’inertie en O d’une sphère pleine S de centre O, de rayon R et de masse m. __ 54
4
Consequences du principe de conservation de la masse __________________________________54 4.1
Torseur cinétique___________________________________________________________ 54 ◊
Expression de la résultante cinétique. _________________________________________ 54
◊
Expression du moment cinétique. ____________________________________________ 55
4.2
Torseur dynamique _________________________________________________________ 55 ◊
Expression de la résultante dynamique. _______________________________________ 55
◊
Relation entre les moments cinétique et dynamique. _____________________________ 56
École Normale Supérieure de Cachan, 61, Avenue du président Wilson, 94235 CACHAN
Préparation à l’Agrégation de Mécanique.
4.3
Énergie cinétique___________________________________________________________ 56
4.4
Autres cas. Exemple : Action de la pesanteur _____________________________________ 57
5
A retenir ______________________________________________________________________58 ◊
Conséquence du principe de conservation de la masse ____________________________ 58
◊
Torseur cinétique_________________________________________________________ 58
◊
Torseur dynamique _______________________________________________________ 58
◊
Energie cinétique_________________________________________________________ 58
CONSERVATION DE L’ENERGIE
1
Introduction____________________________________________________________________61
2
Energétique ____________________________________________________________________61 2.1
Torseur des actions mécaniques extérieures à un solide _____________________________ 61
2.2
Puissance_________________________________________________________________ 61 ◊
Puissance associée à des actions extérieures. ___________________________________ 61
◊
Puissance associée à des actions réciproques.___________________________________ 62
2.3
Travail ___________________________________________________________________ 63
2.4
Energie Potentielle _________________________________________________________ 63 ◊
Energie potentielle associée à des efforts extérieurs ______________________________ 63
◊
Energie potentielle associée à des actions mutuelles______________________________ 64
◊
Exemple : énergie potentielle associée aux inter-efforts gravitationnels. ______________ 64
◊
Quelques actions mutuelles avec énergie potentielle associée ______________________ 65
2.5
Energie cinétique___________________________________________________________ 66
3
Conservation de l’énergie : Théorème de l’énergie cinétique______________________________66 ◊
Théorème de l’énergie pour un solide S _______________________________________ 66
◊
Théorème de l’énergie cinétique pour deux solides S1 et S2 _______________________ 66
◊
Théorème de l’énergie cinétique pour un système Σ de n solides____________________ 67
◊
Intégrale première de l’énergie cinétique : Système conservatif_____________________ 67
4
A retenir ______________________________________________________________________68
PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
1
Introduction____________________________________________________________________70
2
Actions mécaniques ou efforts _____________________________________________________70 2.1
Torseur des actions mécaniques extérieures à un système de solides Σ _________________ 70
2.2
Exemple d’action à distance : la pesanteur _______________________________________ 71
2.3
Actions de contact : Loi de Coulomb ___________________________________________ 71 ◊
Lois de Coulomb relatives à la résultante.______________________________________ 72
École Normale Supérieure de Cachan, 61, Avenue du président Wilson, 94235 CACHAN
Préparation à l’Agrégation de Mécanique. ◊
Lois de Coulomb relatives au moment en I. ____________________________________ 73
◊
Contact non ponctuel, exemple d’une liaison pivot sans frottement __________________ 73
◊
Liaisons parfaites : Actions de contact sans frottement. ___________________________ 74
◊
Lois de comportement de liaisons. ___________________________________________ 75
2.4
Graphe d’analyse___________________________________________________________ 76 ◊
3
Exemple : principe d’un étouffeur de vibrations_________________________________ 76 Principe fondamental de la dynamique _______________________________________________77
3.1
Introduction : un peu d’histoire________________________________________________ 77
3.2
Enoncé du principe fondamental de la dynamique _________________________________ 78
3.3
Conséquences _____________________________________________________________ 78 ◊
Théorème de la résultante dynamique_________________________________________ 78
◊
Théorème du moment dynamique____________________________________________ 79
◊
Théorème des actions mutuelles _____________________________________________ 79
◊
Cas particulier de la statique ________________________________________________ 79
◊
Equations du mouvement __________________________________________________ 79
3.4
Référentiels Galiléens/non Galiléens ___________________________________________ 79 ◊
Référentiels Galiléens _____________________________________________________ 80
◊
Référentiel non Galiléen ___________________________________________________ 80
◊
Exemple : Accélération de la pesanteur. _______________________________________ 81
4
Principe des puissances virtuelles ___________________________________________________83 4.1
Introduction : un peu d’histoire. _______________________________________________ 83
4.2
Enoncé du principe des puissances virtuelles ou PPV. ______________________________ 83
4.3
Choix de torseurs virtuels particuliers et théorèmes de la dynamique __________________ 84
4.3.1
Torseur global quelconque : Equivalence du principe des puissances virtuelles et du principe
fondamental de la dynamique ___________________________________________________________ 84
5
4.3.2
Torseur des vitesses galiléennes : Théorème de l’énergie cinétique __________________ 84
◊
Théorème de l’énergie cinétique pour un solide unique.___________________________ 84
◊
Théorème de l’énergie cinétique pour un système de solides. ______________________ 84
4.3.3
Torseurs de Lagrange : équations de Lagrange__________________________________ 85
◊
Définition des torseurs de Lagrange : _________________________________________ 85
◊
Conséquences ___________________________________________________________ 85
◊
Application du PPV à un unique solide S.______________________________________ 86
◊
Application du PPV à un système de solides Σ : _________________________________ 86
◊
Fonction de force. ________________________________________________________ 87
4.3.4
Combinaison des torseurs de Lagrange : Equation de Painlevé _____________________ 89
◊
Définition du torseur de Painlevé:____________________________________________ 89
◊
Conséquences ___________________________________________________________ 90
◊
Intégrale première de Painlevé ______________________________________________ 90 A retenir ______________________________________________________________________92
École Normale Supérieure de Cachan, 61, Avenue du président Wilson, 94235 CACHAN
Préparation à l’Agrégation de Mécanique. APPLICATION, THEORIE DES PETITS MOUVEMENTS, THEORIE DES CHOCS
1
Introduction____________________________________________________________________96
2
Méthodologie:__________________________________________________________________96
3
Mise en forme du problème. _______________________________________________________97 3.1
Linéarisation ______________________________________________________________ 97
3.1.1
Linéarisation des équations de Lagrange ______________________________________ 97
◊
Premier membre des équations de Lagrange____________________________________ 97
◊
Second membre des équations de Lagrange : Fonction de force_____________________ 98
◊
Second membre des équations de Lagrange : Loi visqueuse________________________ 99
◊
Système d’équations du mouvement après linéarisation :__________________________ 99
3.1.2
Exemple : L’étouffeur de vibrations __________________________________________ 99
◊
Calcul des Qqi : _________________________________________________________ 100
◊
Calcul des Pqi : __________________________________________________________ 101
◊
Equations de Lagrange : __________________________________________________ 101
◊
Equations linéarisées _____________________________________________________ 101
3.2
Rappels : Résolution de systèmes linéaires d’équations différentielles ________________ 102
3.2.1
Résolution numérique. ___________________________________________________ 102
3.2.2
Résolution analytique ____________________________________________________ 103
◊
Cas N°1 : Le second membre F (t) est nul. ____________________________________ 104
◊
Cas N°2 : Le second membre F (t) est nul et C est nulle. _________________________ 105
◊
Cas N°3 : Le second membre F (t) est non-nul et C est nulle. _____________________ 106
◊
Cas N°4 : Le second membre F (t) est nul et C est une combinaison linéaire de M et K. 106
◊
Critère de Routh ________________________________________________________ 107
◊
Rappel : Calcul du déterminant d’une matrice _________________________________ 108
◊
Rappel : Matrice 2x2 valeurs propres, vecteurs propres __________________________ 109
4
Equilibre et stabilité ____________________________________________________________109 4.1
Introduction______________________________________________________________ 109
4.2
Systèmes conservatifs ______________________________________________________ 109 ◊
Equilibre ______________________________________________________________ 110
◊
Stabilité : Définition _____________________________________________________ 110
◊
Stabilité : Théorème de Lejeune Dirichlet_____________________________________ 110
◊
Extension aux systèmes visqueux ___________________________________________ 111
◊
Equilibre et stabilité pour un système à un seul paramètre ________________________ 112
◊
Equilibre et stabilité pour un système à n paramètres ____________________________ 112
◊
Exemple, étouffeur de vibration ____________________________________________ 113
4.3
Cas général, méthode directe ou méthode de Liapounov ___________________________ 114 ◊
Etat de mouvement ______________________________________________________ 114
◊
Détermination de l’équilibre par la méthode directe _____________________________ 114
◊
Stabilité au sens de Liapounov _____________________________________________ 114 École Normale Supérieure de Cachan, 61, Avenue du président Wilson, 94235 CACHAN
Préparation à l’Agrégation de Mécanique. ◊
Théorème de Liapounov __________________________________________________ 114
◊
Stabilité asymptotique____________________________________________________ 115
◊
Stabilité orbitale ________________________________________________________ 115
5
Vibrations autour d’une position d’équilibre stable ____________________________________115 ◊
Introduction ____________________________________________________________ 115
◊
Vibrations libres ________________________________________________________ 116
◊
Vibrations forcées _______________________________________________________ 116
◊
Vibrations amorties ______________________________________________________ 116
6
Mécanique des Chocs-Percussions _________________________________________________116 6.1
Introduction______________________________________________________________ 116
6.2
Cas d’un point matériel _____________________________________________________ 117
6.3
Cas d’un solide ou d’un système de solides _____________________________________ 117 ◊
Remarque 1 : ___________________________________________________________ 118
◊
Remarque 2 : ___________________________________________________________ 118
◊
Remarque 3 : ___________________________________________________________ 118
6.3.2
Percussion de liaison_____________________________________________________ 118
6.3.3
Choc sans frottement entre deux solides ______________________________________ 119
◊
Définition de e__________________________________________________________ 119
◊
Propriété de e et P : 0 ≤ e ≤ 1 et 0 ≤ P ________________________________________ 119
◊
Cas particuliers importants ________________________________________________ 121
École Normale Supérieure de Cachan, 61, Avenue du président Wilson, 94235 CACHAN
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables INTRODUCTION
1
Introduction _______________________________________________________ 2 ◊
Plan du cours_____________________________________________________________ 2
◊
Exemple 1 : Vibration d’une aile d’avion _______________________________________ 3
◊
Exemple 2 : Ecoulement d’un milieu granulaire __________________________________ 4
◊
Sources bibliographiques ___________________________________________________ 4
1
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables INTRODUCTION
1
INTRODUCTION
La mécanique du solide indéformable, comme les autres branches de la mécanique, procède d’une schématisation des mouvements réels ou potentiels à l’intérieur du système étudié. Le choix d’un schéma cinématique plutôt qu’un autre dépend du niveau de simplification recherché, des matériaux et de l’échelle à laquelle le problème est traité (voir les deux exemples ci-dessous). Ainsi si le champ des vitesses eulériennes à l’intérieur du système étudié est : Champ équiprojectif
Solide indéformable
Champ équiprojectif par morceaux tridimensionnels
Système de solides indéformables
Champ équiprojectif par morceaux bidimensionnels
Poutres de la mécanique des structures
Champ équiprojectif par morceaux unidimensionnels
Plaques et coques
Des objets mathématiques bien adaptés à chacun de ces schémas cinématiques (mécanique des milieux continus, coques, plaque, poutres, solides) ont été développés afin de pouvoir exprimer les principes fondamentaux de la mécanique sous forme d’équations. L’objet mathématique privilégié de la mécanique des solides indéformables est le torseur. Les principes fondamentaux sont les suivants : • • • •
Conservation de la masse. Conservation de l’énergie (premier principe de la thermodynamique) Conservation de la quantité de mouvement (d’Alembert). Second principe de la thermodynamique.
L’écriture des trois premiers principes conduit systématiquement à un système d’équations pour lequel le nombre d’équations est inférieur au nombre d’inconnues. Les équations complémentaires sont données par les lois de comportement, dont on s’assure qu’elles permettent de vérifier le second principe de la thermodynamique. Ces lois de comportement seront par exemple dans le cadre de la mécanique du solide indéformable : • • • •
Comportement rigide indéformable pour les solides Lois de contact entre solides (lois de Coulomb) Comportement de liaisons entre solides (liaison parfaites ou liaisons élastiques). Lois d’action à distance (attraction gravitationnelle, par exemple)
L’objet de ce cours est d’apporter les outils et les méthodes de travail permettant la résolution de problèmes mécaniques dans le cadre de la mécanique du solide indéformable. Pour cela, la première partie sera consacrée à la description de la cinématique dans le cadre de la mécanique des solides indéformables. Puis les principes fondamentaux seront exprimés en utilisant le formalisme associé à ce schéma cinématique. Il n’y aura pas de chapitre spécifique consacré aux lois de comportement, étant donné que les lois usuellement employées pour décrire les interactions entre solides sont peu nombreuses et bien connues. Une fois que les outils permettant de mettre les problèmes en équations auront été présentés, des méthodes de résolution des systèmes d’équations obtenus seront présentées, dans le cadre des petits mouvements autour d’une position connue.
◊ Plan du cours • • •
•
A – Introduction. B – Rappels de mathématiques. C – Schématisation de la cinématique. C1 – Cinématique du solide ind éformable. C2 – Cinématique des systèmes de solides indéformables. D – Expression des principes fondamentaux. D1 – Conservation de la masse.
2
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables INTRODUCTION
•
D2 – Conservation de l’énergie. D3 – Conservation de la quantité de mouvement. E – Méthodes de résolution. E1 – Equilibre et Stabilité E2 – Vibrations libres ou forcées E3 – Chocs
◊ Exemple 1 : Vibration d’une aile d’avion
Le dimensionnement mécanique d’une aile d’avion se fait dans le cadre de la mécanique des structures. L’aile peut être schématisée comme une poutre à section complexe et variable, et la portance, proportionnelle au carré de la vitesse d’une section, peut être assimilée à une charge linéique. L’aile se fléchit significativement en vol, l’amplitude de battement en bout d’aile est d’environ un mètre en fonctionnement normal mais peut être bien plus élevée, après un trou d’air, par exemple. L’approche « poutre » permet de calculer le moment fléchissant à l’attache de l’aile sur la cellule en fonction de la déflection β en bout d’aile, c'est-à-dire de calculer la raideur de la structure. Si l’on connaît la déflection maximale β, on peut en déduire les contraintes au niveau de l’attache de l’aile sur la cellule et dimensionner cette attache. Cependant, pour estimer cette déflection maximale β, il est nécessaire de connaître le comportement dynamique de l’avion complet. En effet, les moteurs, par exemple, ont une masse très importante par rapport à celle de l’avion (masse d’un moteur CFM56-3 = 2 tonnes, masse d’un A320 hors moteur : 37 tonnes). Au cours de certaines manœuvres, des oscillations des moteurs, couplés au battement des ailes peuvent apparaître. Pour connaître, par exemple, la déflection maximale β en bout d’aile et dimensionner l’attache de l’aile sur la cellule, il faut connaître l’amplitude de ces mouvements. Ceci nécessite une étude du comportement dynamique de l’avion complet, qui inclue la cellule, les ailes, l’empennage, les moteurs et même le chargement de l’avion (mobile ou non). Pour ce type d’étude, le détail des déformations internes à chacun des éléments de l’avion est négligé, pour se limiter à l’étude des mouvements relatifs entre ces éléments. L’aile par exemple pourra être assimilée à un ou plusieurs éléments rigides attachés à l’avion et liés entre eux par des liaisons pivot élastiques. Le moteur pourra être assimilé à une masse M, attachée à l’aile par une liaison pivot élastique. Dans cet exemple particulier, les ailes sont déformables, mais sont assimilées à un système de solides rigides afin de pouvoir traiter le comportement dynamique de l’avion complet. Le choix des « solides » du système résulte donc d’une schématisation du champ de déplacement des points de l’aile.
3
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables INTRODUCTION
◊ Exemple 2 : Ecoulement d’un milieu granulaire
L’écoulement d’un milieu granulaire (sable, poudres …), est traité de manière différente selon le point de vue de l’observateur. A l’échelle de l’écoulement complet, le milieu granulaire peut être traité en première approximation comme un matériau continu déformable. On utilisera donc soit la mécanique des milieux continus, soit la mécanique des fluides pour traiter le problème, en utilisant des lois de comportement appropriées. A l’échelle des grains, le milieu granulaire est un empilement irrégulier de grains. Individuellement, les grains peuvent être considérés comme des solides indéformables. Ces solides sont en contact et glissent ou roulent les uns par rapport aux autres. Lors de l’écoulement et après l’écoulement, les grains s’arrangent en voûtes entre lesquelles il reste du vide. Cet édifice peut être déstabilisé. C’est ce qui se produit par exemple sur une pente enneigée lors d’une avalanche. L’étude des conditions de stabilité des édifices de grains se fait à l’aide de la mécanique des solides indéformables, avec des lois de contact entre solides appropriées. Dans cet exemple particulier, le champ est « solidifiant » par morceaux, c'est-à-dire sur chacun des grains. Le choix de la mécanique des milieux continus ou de la mécanique des solides se fait en fonction de la « dimension des morceaux » vis-à-vis de l’échelle du problème à traiter.
◊ Sources bibliographiques •
• •
• •
•
•
Mécanique 1,2,3, Cours et exercices (1995), Yves Brémont, Paul Réocreux, collection Sciences Industrielles, Ed. Ellipses, Paris. Dynamique, Cours et exercices, (2002) Robert Lassia et Christophe Bard, Ed. Ellipses, Paris. Cours de Physique, tome I : Mécanique, (1965), C. Kittel, W. Knight, M. Ruderman, traduit par P. Lallemand, Ed Dunod, Paris. Mécanique Classique, (1971), J.J. Moreau, Ed. Masson et Cie, Paris. Mécanique générale, (2001), Jean Claude Bône, Michel Boucher, Polycopiés de tronc commun, ECP. Quelques Compléments de Mécanique générale, (1997), Jean Pierre Pelle, polycopié ENS pour la préparation à l’agrégation. Equations différentielles et systèmes dynamiques, (1999), J. Hubbard, B. West, traduit par V. Gautheron, Ed. Cassini, Paris.
4
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES
1
Calcul Vectoriel ____________________________________________________ 6 1.1
Opérations sur les vecteurs ____________________________________________ 6
◊ ◊ ◊ ◊ 1.2
1.3
Glisseur _________________________________________________________________ Moment en un point d’un glisseur_____________________________________________ Moment d’un glisseur par rapport à un axe ______________________________________ Ensembles de glisseurs _____________________________________________________
8 8 9 9
Champ equiprojectif de vecteurs ______________________________________ 10
◊ ◊ 1.4
Définition ______________________________________________________________ 10 Propriétés_______________________________________________________________ 10
Torseurs __________________________________________________________ 10
◊ ◊ ◊ ◊ ◊ ◊ 1.5
Définition ______________________________________________________________ Champ des moments d’un torseur. ___________________________________________ Torseur associé à un ensemble de glisseur _____________________________________ Invariants du torseur.______________________________________________________ Point central, axe central, moment central d’un torseur ___________________________ Symétrie du champ des moments d’un torseur. Origine du mot « torseur ». ___________
11 11 11 11 12 12
Opérations sur les torseurs ___________________________________________ 13
◊ ◊ ◊ ◊ ◊
Addition________________________________________________________________ Multiplication par un réel __________________________________________________ Décomposition __________________________________________________________ Produit ou co-moment de deux torseurs _______________________________________ Torseur à structure________________________________________________________
13 13 13 14 14
Dérivation vectorielle_______________________________________________ 15 2.1
Dérivée d’un vecteur ________________________________________________ 15
◊ ◊ 2.2
Définition ______________________________________________________________ 15 Propriétés_______________________________________________________________ 15
Changement de base de dérivation_____________________________________ 15
◊ ◊ ◊ ◊
3
6 6 7 7
Champs de vecteurs __________________________________________________ 8
◊ ◊ ◊ ◊
2
Produit scalaire ___________________________________________________________ Produit vectoriel __________________________________________________________ Produit mixte_____________________________________________________________ Division vectorielle ________________________________________________________
Vocabulaire _____________________________________________________________ Dérivée d’un vecteur exprimé dans la base de dérivation __________________________ Dérivée d’un vecteur exprimé dans une base distincte de la base de dérivation. ________ Propriétés du vecteur vitesse de rotation. ______________________________________
15 16 16 17
A retenir _________________________________________________________ 18 ◊ Champ équiprojectif : _____________________________________________________ 18 R Torseur : {T }= ◊ _____________________________ _____________________ 18 M A A ◊ ◊ ◊ ◊
Produit de deux torseurs définis au même point A : ______________________________ Invariants du torseur ______________________________________________________ Axe et point centraux _____________________________________________________ Changement de base de dérivation ___________________________________________
5
18 18 18 18
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES
1
1.1
◊ -
CALCUL VECTORIEL
Opérations sur les vecteurs
Produit scalaire Définition Le produit scalaire des deux vecteurs U et V est le nombre réel suivant noté U.V :
U ⋅ V = U V Cos (U ,V )
•
-
-
Propriétés : •
Symétrie : U ⋅ V = V ⋅ U
•
Distributivité : U ⋅ (V + W ) = U ⋅ V + U ⋅ W
•
Multiplication par un réel : λ U ⋅ α V = λα U ⋅ V
Expression analytique Dans une base orthonormée (x,y,z) le produit scalaire des deux vecteurs V1(x1,y1,z1) et V2(x2,y2,z2) s’écrit :
V 1 ⋅ V 2 = x1 . x 2 + y1 . y 2 + z 1 . z 2
◊ -
Produit vectoriel Définition Le produit vectoriel des deux vecteurs U et V est le vecteur noté U ∧ V tel que, U ∧ V soit
perpendiculaire au plan (U,V), le trièdre (U,V, U ∧ V ) soit direct, et la norme de U ∧ V soit égale à : •
-
U ∧ V = U V Sin(U , V )
Interprétation géométrique La norme du produit vectoriel U ∧ V , représente la surface du parallélogramme défini par les deux vecteurs U et V :
-
Propriétés •
Antisymétrie : U ∧ V = −V ∧ U
•
Distributivité par rapport à l’addition : U ∧ (V + W ) = U ∧ V + U ∧ W
•
Multiplication par un réel : λ U ∧ α V = λα (U ∧ V )
6
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES •
•
-
Application à une base orthonormée directe (x,y,z) :
x ∧ x = 0
y ∧ x = − z z ∧ x = y
x ∧ y = z
y ∧ y = 0
z ∧ y = − x
x ∧ z = − y
y ∧ z = x
z ∧ z = 0
Double produit vectoriel (Formule de Gibbs) : U ∧ (V ∧ W ) = (U ⋅ W ) V − (U ⋅ V )W
Expression analytique Dans une base orthonormée (x,y,z) le produit vectoriel des deux vecteurs V1(x1,y1,z1) et V2(x2,y2,z2) s’écrit :
V 1 ∧ V 2 = ( y1 . z 2 − z 1 . y 2 ). x + ( z 1 . x 2 − x1 . z 2 ). y + ( x1 . y 2 − y1 . x 2 ). z
◊ -
Produit mixte Définition Le produit mixte des trois vecteurs U, V et W est le nombre réel suivant, et noté (U,V,W) : •
-
(U ,V ,W ) = U ⋅ (V ∧ W )
Interprétation géométrique La valeur absolue du produit mixte (U,V,W) représente le volume du parallélépipède défini par les trois vecteurs.
-
-
Propriétés •
Permutation des opérateurs : (U , V , W ) = U ⋅ (V ∧ W ) = (U ∧ V ) ⋅ W
•
Distributivité par rapport à l’addition : (U + X , V , W ) = (U , V , W ) + ( X , V , W )
•
Multiplication par un réel : (λ U , α V , γ W ) = λαγ (U , V , W )
•
Permutation des vecteurs : (U , V , W ) =
•
Permutation circulaire :
−(V , U , W ) (U ,V ,W ) = (V ,W ,U ) = (W ,U ,V )
Expression analytique. Dans une base orthonormée (x,y,z) le produit mixte des trois vecteurs V1(x1,y1,z1), V 2(x2,y2,z2) et V3(x2,y2,z2) se calcule comme le déterminant suivant :
x1 x 2
(V 1 ,V 2 ,V 3 ) = y1
y 2 y 3 = − x3 . y 2 . z 1 + x 2 . y 3 z . 1 + x3 . y1 . z 2 − x1 . y 3 . z 2 − x 2 . y1 z . 3 + x1 . y 2 . z 3
z 1 z 2
◊ -
x3 z 3
Division vectorielle Définition
7
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES
Soient deux vecteurs A et B non nuls et orthogonaux, le résultat de la division vectorielle est l’ensemble des vecteurs X tels que : •
-
A ∧ X = B
Solution générale L’ensemble X est défini de la manière suivante, α étant un réel : •
1.2
◊ -
X =
B ∧ A A ⋅ A
+ α . A
Champs de vecteurs Glisseur Définition Un glisseur est défini par un vecteur V et un point P quelconque de son support et noté (P,V). Un représentant de ce glisseur est un bipoint, appartenant à la droite (D) définie par (P,V). Ici, par exemple, le bipoint AB ou le bipoint CD.
◊ -
Moment en un point d’un glisseur Définition On appelle moment au point A du glisseur (P,V), noté MA(P,V) ou MA(V) le vecteur suivant : •
-
M A (V ) = AP ∧ V = V ∧ PA
Interprétation géométrique Soit H le pied de la perpendiculaire abaissée de A sur la droite (D) définie par le glisseur. La norme du moment du glisseur (P,V) au point A, est égale à :
8
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES
M A ( P ,V ) = AH V
•
-
Propriétés Le moment au point A du glisseur (P,V) est indépendant du choix du point P sur le support (D) du glisseur. Champ des moments :
•
•
M B ( P , V ) = BP ∧ V = ( BA + AP ) ∧ V = AP ∧ V + BA ∧ V = M A ( P ,V ) + BA ∧ V M B ( P ,V ) = M A ( P ,V ) + BA ∧ V
◊ -
Moment d’un glisseur par rapport à un axe Définition On appelle moment par rapport à l’axe ∆(A,x) du glisseur (P,V) le nombre réel suivant : •
-
M ∆ ( P ,V ) = x ⋅ M A ( P ,V ) = ( x, AP , V )
Interprétation géométrique Le moment du glisseur (P,V) par rapport à l’axe ∆(A,x) est égal au produit du bras de levier OH, par la composante W du vecteur V, perpendiculaire à la fois au bras de levier et à l’axe. Sur la figure ci-dessous :
M ∆ ( P , V ) = HO W
◊ -
et
M A ( P ,V ) = AN V
Ensembles de glisseurs Ensemble fini de glisseurs. Si l’on considère un ensemble fini de n glisseurs (Pi,Vi), deux grandeurs peuvent être définies, la résultante, R, et le moment au point A, MA, de l’ensemble fini de glisseurs. •
R =
n
∑V i
i =1 •
M A =
n
∑ AP ∧ V i =1
i
i
Alors le champ des moments de l’ensemble fini de glisseurs vérifie aussi : •
M B = M A + R ∧ AB
9
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES
-
Ensemble infini de glisseurs. Si l’on considère un ensemble infini de glisseurs (P,F(P)), où F(P) est une densité de champ de vecteurs définie en tout point P d’un domaine E. Deux grandeurs peuvent être définies, la résultante, R, et le moment au point A, MA, de l’ensemble infini de glisseurs.
R =
•
∫ F ( P )dv P ∈ E
M A =
•
∫ AP ∧ F ( P )dv P ∈ E
Alors le champ des moments de l’ensemble fini de glisseurs vérifie encore :
M B = M A + R ∧ AB
•
1.3
◊
Champ equiprojectif de vecteurs Définition Un champ de vecteur V est equiprojectif si : ∀ A, ∀ B
◊
Propriétés •
•
1.4
AB ⋅ V A = AB ⋅ V B
Si un champ de vecteur equiprojectif est connu en trois points non alignés de l’espace, alors il est connu en tout point P (voir figure ci-dessous) Par ailleurs si deux champs de vecteur V1 et V2 sont equiprojectif alors aV1+bV2 est equiprojectif aussi quel que soient les deux réels a et b choisis.
Torseurs Le torseur est l’outil privilégié de la mécanique du solide. Il est utilisé pour représenter le mouvement d’un solide, à caractériser une action mécanique, à formuler le principe fondamental de la dynamique de manière générale, etc…
10
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES
◊
Définition Un torseur est un ensemble défini par ses d eux éléments dits « éléments de réduction » : •
Un vecteur noté R appelé la résultante du torseur
•
Un champ de vecteur M vérifiant la relation : ∀ A, ∀ B M B
= M A + R ∧ AB
MA est appelé le moment au point A du torseur T Le torseur T se note de la façon suivante au point A :
• •
R
{T }= M A
◊
A
Champ des moments d’un torseur. Le champ des moments d’un torseur est equiprojectif et réciproquement, tout champ de vecteur equiprojectif est le champ des moments d’un torseur.
-
Démonstration :
R R = M M A A B B
Si l’on prend le Torseur T tel que : {T }=
Le champ des moments de ce torseur vérifie par définition : M B
= M A + BA ∧ R
En appliquant un produit scalaire par AB a cette relation, on retrouve bien que ce champ des moments est equiprojectif.
M B ⋅ AB = M A ⋅ AB + ( BA ∧ R ) ⋅ AB 14 4 244 3
=0
-
Remarque : Si un solide est indéformable, le champ des vecteurs vitesse des points de ce solide est nécessairement equiprojectif. Par conséquent il est représentable par un torseur, dont le vecteur moment est le vecteur vitesse du point considéré. On verra plus loin que le résultante du torseur est le vecteur rotation de ce solide.
◊
Torseur associé à un ensemble de glisseur Un ensemble de glisseur fini ou infini répond à la définition du torseur, par conséquent :
•
•
◊
n ∑ V i Ensemble fini de glisseurs : {T }= n i =1 ∑ AP i ∧ V i A i =1 F ( P ) dv ∫ Ensemble infini de glisseurs : {T }= P ∈ E ( ) AP F P dv ∧ ∫ A P ∈ E
Invariants du torseur. Entre deux points quelconques A et B de l’espace, deux composantes du torseur sont conservées, qui constitue les deux invariants du torseur :
11
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES • •
Premier invariant : La résultante R Second invariant : La projection du moment du torseur sur sa résultante :
∀ A, ∀ B M B = M A + R ∧ AB ⇒ R ⋅ M A = R ⋅ M B
◊ -
Point central, axe central, moment central d’un torseur Point central : Point où le moment du torseur à la même direction que la résultante.
-
Axe central : Ensemble des points centraux On se propose de montrer que les points centraux sont alignés sur un même axe, pour un torseur T, qui se note au point A :
R
{T }= M A
A
Le moment au point A du torseur peut se décomposer en deux termes, U et W, où U est la composante de MA selon R et W est orthogonal a R alors :
M A = U + W avec U ∧ R = 0 et W ⋅ R = 0 Si B est un point central, du fait du second invariant, ∀ A, M B Comme M B
= U .
= M A + R ∧ AB ⇒ U = U + V + R ∧ AB ⇒ R ∧ AB = −V . Par division
vectorielle, on en déduit :
AB =
R ∧ M A R 2
+ α R
ainsi si l’on pose : AH =
R ∧ M A R 2
alors les points centraux
s’alignent sur une droite de même direction que la résultante du torseur R et passant par le point H.
-
Moment central : Le moment central est le moment du torseur en un point quelconque de son axe central. La norme du moment d’un torseur est minimale pour les points centraux. Par conséquent si le moment d’un torseur est nul en un point, ce point appartient à l’axe central. L’axe central se définit alors à l’aide de ce point et de la résultante.
◊
Symétrie du champ des moments d’un torseur. Origine du mot « torseur ». Soit R(B,x,y,z) un repère orthonormé direct, dont l’axe (B,z) est confondu avec l’axe central d’un torseur T. Posons alors : r
R = R. z {T }= r M M . z = B B B Si l’on choisit un axe (H,u) quelconque dans un plan orthogonal à z et qui rencontre l’axe (B,z) au point H, et un nouveau repère associé à cet axe R’(H,u,v,z), alors pour un point A quelconque appartenant à cet axe : r
r
r
r
BH = h. z et HA = r .u ⇒ M A = M B + A B ∧ R = Mb. z + R.r .v D’où l’expression du moment du torseur T au point A : r
M A = Mb. z + R.r .v
12
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES
• •
•
1.5
◊
Lorsque la distance r de A à l’axe central (B,z) est nulle MA=MB. Lorsque la distance r augmente le moment du torseur au point A tourne progressivement dans le plan (v,z) jusqu’à s’aligner avec la direction v. Ainsi, observe t’on une « torsion » du moment du torseur au point A, lorsque le point A s’éloigne de l’axe central du torseur, d’où l’origine du mot torseur.
Opérations sur les torseurs Addition La somme de deux torseurs {T} et {T’} est le torseur {T+T’}. Pour faire la somme de deux torseurs, il faut au préalable les écrire au même point :
R
R' et {T '}= ' M A A A A
{T }= M
{T } + {T '} = {T + T '}
On vérifie ensuite que {T+T’} est bien un torseur. C'est-à-dire que son champ des moments vérifie bien la relation suivante :
M B {T + T '} = M A {T + T '} + BA ∧ R{T + T '} ∀ A, ∀ B Démonstration :
M B {T + T '} = M B {T } + M B {T '} = ( M A {T } + BA ∧ R{T }) + ( M A {T '} + BA ∧ R{T '}) M B {T + T '} = M A {T } + M A {T '} + BA ∧ ( R{T } + R{T '}) = M A {T + T '} + BA ∧ R{T + T '}
◊
Multiplication par un réel Soit {T} un torseur et λ un réel :
λ . R λ .{T }= = {λ T } λ M A A
◊ -
Décomposition Torseur Couple Un torseur couple est un torseur dont la résultante est nulle.
0
{T }= M A
A
13
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES
Le moment d’un torseur couple est le même en tout point de l’espace. Un torseur couple peut être représenté par un ensemble de glisseur de direction parallèle, de même norme et de sens contraire. En effet, si l’on considère deux glisseurs (P,V) et (A,-V), alors le torseur associé à cet ensemble de glisseur vaut :
0 − V V + = = ∧ M AP V ∧ 0 AP V A A A A
{T }= -
Torseur à résultante Un torseur à résultante est un torseur dont le moment central est nul.
R
{T }= 0 B
-
avec R ≠ 0
et B ∈ ∆ , où ∆ est l’axe central du torseur.
Décomposition d’un torseur Tout torseur est donc en général la somme d’un torseur couple et d’un torseur résultante. En effet, on peut écrire:
R
0
R
{T }= = + M M 0 A
◊
A
A
A
A
Produit ou co-moment de deux torseurs Le produit de deux torseurs {T} et {T’} définis au même point A, est le nombre réel suivant :
R
R'
{T }= et {T '}= ' M M A
A
A
A
{T }⋅ {T '} = R ⋅ M A' + R '⋅ M A
En outre, le produit de deux torseurs est commutatif : {T } ⋅ {T '} = {T '} ⋅ {T }
◊
Torseur à structure Un torseur à structure est un torseur défini à partir d’une densité de champ infini de vecteurs F(P), c'est-à-dire de la forme :
F ( P ) dv P ∈∫ E {T [ F ( P )]}= ∫ AP ∧ F ( P )dv A P ∈ E Le produit d’un torseur à structure par un torseur quelconque se met sous la forme :
F ( P ) dv ∫ R {T [ F ( P )]}⋅ {T }= P ∈ E ⋅ = R ⋅ ∫ AP ∧ F ( P )dv + M A ⋅ ∫ F ( P ) dv M ( ) ∧ AP F P dv ∫ A A P ∈ E P ∈ E ∈ P E A
{T [ F ( P )]}⋅ {T } =
∫ ( R, AP , F ( P ))dv + ∫ ( M
P
P ∈ E
{T [ F ( P )]}⋅ {T } =
P ∈ E
∫ ( R, AP , F ( P ))dv + ∫ M
P
P ∈ E
+ R ∧ PA) ⋅ F ( P ) dv
P ∈ E
14
⋅ F ( P )dv +
∫ ( R, PA, F ( P ))dv P ∈ E
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES
∫ M
Soit finalement : {T F ( P ) }⋅ {T } =
[
]
P
⋅ F ( P )dv
P ∈ E
2
2.1
◊
DERIVATION VECTORIELLE
Dérivée d’un vecteur
Définition Par définition la dérivée d’un vecteur V(t) par rapport à la variable t, dans l’espace vectoriel E est le vecteur suivant :
V (t + h ) − V (t ) d = V lim dt h h→ 0 E {
Par conséquent, la dérivée d’un vecteur V, dépend du choix de l’espace vectoriel de référence dans lequel est exprimé le vecteur. En pratique, il est donc nécessaire de toujours préciser par rapport à quel référentiel du mouvement est effectuée la dérivée.
◊
Propriétés •
•
•
•
•
d d d ( ) + = + V V V 1 2 1 dt dt dt V 2 R R R df d d ( ) ( ) = + f t V f t V V 1 Produit par une fonction scalaire f : 1 dt dt dt R R d d d Dérivée du produit scalaire : (V 1 ⋅ V 2 ) = V 1 ⋅ V 2 + V 1 ⋅ V 2 dt dt R dt R d d d ( ) ∧ = ∧ + ∧ V V V V V Dérivée d’un produit vectoriel : 2 1 2 1 dt 1 dt V 2 R dt R R Somme :
Dérivée d’un produit mixte :
d dt •
2.2
◊
d d d V 1 , V 2 , V 3 + V 1 , V 2 , V 3 + V 1 , V 2 , V 3 dt R dt R dt R
(V 2 , V 2 , V 3 ) =
Dérivée d’une fonction de fonction :
d d d θ ( ) = V θ t [ ] dt d θ V dt R R
Changement de base de dérivation
Vocabulaire La base dans laquelle on exprime les vecteurs sera indifféremment appelée, base de calcul ou base de projection.
15
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES
La base dans laquelle est effectuée la dérivation, sera indifféremment appelée base de dérivation, référentiel du mouvement ou repère.
◊
Dérivée d’un vecteur exprimé dans la base de dérivation Dans ce cas particulier, la base de projection est confondue avec le référentiel du mouvement choisi. Alors, si un vecteur V s’exprime dans cette base R(0,x,y,z) à l’aide de trois composantes a,b,c, comme les trois vecteurs unitaires de cette base sont constants :
da db dc d = + V x y z dt dt dt dt R
◊
Dérivée d’un vecteur exprimé dans une base distincte de la base de dérivation. Supposons une base de projection R 1(x1,x2,x3) dans laquelle le vecteur V est exprimé à l’aide de trois composantes (a1,a2,a3). Supposons aussi une base R(e1,e2,e3) attachée au référentiel du mouvement et distincte de la première. Lors de la dérivation du vecteur V par rapport au référentiel du mouvement R, il faut tenir compte du fait que les vecteurs unitaires de la base R 1 dans laquelle est exprimé le vecteur V ne sont pas constants dans la base de dérivation R. Ainsi : 3 dai dxi d = + V x a ∑ dt i i dt dt i =1 R R
Soit en rassemblant les termes : 3 d d dt V = dt V + ∑ ai R R1 i =1
dxi dt R
A ce stade nous avons besoin de l’expression des dérivées des vecteurs unitaires de la base R1 par rapport au référentiel du mouvement R. L’orientation d’un base par rapport à une autre se défini à l’aide trois angles de rotation (αk , k=1,3). Alors :
∂ x [dxi ] R = ∑ i d α k ∀i k =1 ∂α k R 3
et
3 ∂ xi d α k dxi = ∀i ∑ ∂ dt α dt R k =1 k R
Par ailleurs, les paramètres αk étant des angles de rotation, on a :
∂ xi = ek ∧ xi . α ∂ k R Par suite :
3 d α k dxi ek ∧ xi ∀i = ∑ dt dt k =1 R
Si l’on définit un vecteur Ω de la façon suivante : 3
Ω( R1 R ) = ∑ r
k =1
r dx ek il vient i = Ω( R1 R ) ∧ xi dt dt R
d α k
On en déduit alors la formule de changement de base de dérivation : r r d r d r = + Ω ∧ ( R1 R ) V dt V dt V R R1
16
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES
Où Ω est le vecteur vitesse de rotation de la base R1 par rapport à la base R.
◊
Propriétés du vecteur vitesse de rotation. r
-
Composition des rotations :
r
r
Ω( R3 R1) = Ω( R3 R 2) + Ω( R 2 R1)
Etant donné un vecteur V, on peut écrire successivement : r r d r d r = + Ω ∧ V R 2 R 1 V ( ) V dt dt R1 R 2 r r d r d r ( ) = + Ω ∧ V V R 3 R 2 V dt dt R 2 R 3
Alors : r r r d r d r ( ) ( ) = + Ω + Ω ∧ V V R 3 R 2 R 2 R 1 V ( ) dt dt R1 R 3
Soit : r
r
r
Ω( R3 R1) = Ω( R3 R 2 ) + Ω( R 2 R1) r
-
Inversion des bases de dérivations :
r
Ω( R1 R 2) = −Ω( R 2 R1)
r r d r d r ( ) = + Ω ∧ V V R 2 R 1 V dt dt R1 R 2 r r d r d r = − Ω ∧ V V R 2 R 1 V ( ) dt dt R 2 R1
Soit : r
r
Ω( R1 R 2) = −Ω( R 2 R1)
17
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables RAPPELS DE MATHEMATIQUES
3
◊
A RETENIR
Champ équiprojectif : •
◊
◊
Un champ de vecteur V est equiprojectif si : ∀ A, ∀ B
•
R (résultante), M (champ des moments) tel que: ∀ A, ∀ B M B
•
Le champ des moments d’un torseur est equiprojectif
r
r
= M A + R ∧ AB
Produit de deux torseurs définis au même point A :
R
A
R'
A
A
A
{T }⋅ {T '} = R ⋅ M A' + R '⋅ M A
Invariants du torseur • •
◊
r
R Torseur : {T }= M A A
{T }= et {T '}= ' M M ◊
r
A B ⋅ V A = A B ⋅ V B
Premier invariant : La résultante R Second invariant : La projection du moment du torseur sur sa résultante.
Axe et point centraux Un point est dit central pour le torseur T, si en ce point son moment et sa résultante ont même direction. Les points centraux s’alignent sur un axe, dit axe central. La norme du moment est minimale sur l’axe central. La direction de l’axe central est celle de la résultante, et l’axe passe par le point H, défini à partir d’un point A quelconque par :
R
{T }= M A
◊
A
AH =
R ∧ M A R 2
Changement de base de dérivation Si le mouvement d’une base R 1(x1,x2,x3) par rapport à un référentiel R(e1,e2,e3) est défini par trois angles (αk , k=1,3). •
∂ xi = ek ∧ xi α ∂ k R
Dérivée des vecteurs de la base R1 par rapport au référentiel R : 3
Ω( R1 R ) = ∑ r
•
Vecteur vitesse de rotation de R1/R :
k =1 r
•
•
r
d α k dt r
ek
Ω( R3 R1) = Ω( R3 R 2) + Ω( R 2 R1) r r d r d r = + Ω ∧ V V R 1 R V ( ) Changement de base de dérivation : dt dt R R1 Composition des rotations :
18
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
1
Cinematique du point ______________________________________________ 20
◊ ◊ ◊ 2
Le solide indeformable _____________________________________________ 20
◊ ◊ 3
Définition des coordonnées de l’origine d’un repère. ______________________ 21
◊ ◊ ◊ 3.2
Coordonnées cartésiennes __________________________________________________ 21 Coordonnées cylindriques __________________________________________________ 22 Coordonnées sphériques ___________________________________________________ 22
Définition de l’orientation relative de deux bases _________________________ 23
◊ ◊ ◊
Angles d’Euler __________________________________________________________ 23 Angles de Cardan ________________________________________________________ 24 Quaternions _____________________________________________________________ 25
Cinematique du solide ______________________________________________ 28 4.1
Introduction, notations ______________________________________________ 28
4.2
Champ des vecteurs vitesse des points d’un solide : torseur cinématique _____ 28
◊ ◊ ◊ ◊ 4.3
5
Définition ______________________________________________________________ 20 Equivalence Repère-Solide _________________________________________________ 21
Paramétrage de la position relative de deux solides_______________________ 21 3.1
4
Vecteur position d’un point _________________________________________________ 20 Vecteur vitesse d’un point__________________________________________________ 20 Vecteur accélération d’un point _____________________________________________ 20
Propriété : Equiprojectivité. ________________________________________________ Calcul du vecteur rotation instantanée. ________________________________________ Exemple : Mouvement de translation _________________________________________ Exemple : Mouvement de rotation instantanée __________________________________
29 29 30 30
Champ des vecteurs accélération des points d’un solide.___________________ 30
Composition des mouvements ________________________________________ 31 5.1
Introduction _______________________________________________________ 31
5.2
Composition des vecteurs vitesse ______________________________________ 32
◊ ◊ ◊ ◊ ◊ 5.3
Relation de composition des vecteurs vitesses __________________________________ Définitions______________________________________________________________ Généralisation ___________________________________________________________ Composition des torseurs cinématiques _______________________________________ Exemple________________________________________________________________
32 32 33 33 33
Composition des vecteurs accélération _________________________________ 34
◊ ◊
Relation de composition des vecteurs accélération _______________________________ 34 Définitions______________________________________________________________ 35
6 A retenir _________________________________________________________ 35
◊ ◊ ◊ ◊
Cinématique du point : ____________________________________________________ Cinématique du solide : ____________________________________________________ Formules de changement de point ____________________________________________ Formules de composition des mouvements. ____________________________________
19
35 35 36 36
Préparation à l’agrégation l’agrégation de mécanique mécanique : Polycopié Polycopié de Mécanique Mécanique des Solides Indéformab Indéformables les CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
1
CINEMATIQUE DU POINT
On rappelle ici les définitions de la position, de la vitesse et de l’accélération d’un point P point P par rapport à un repère R(O,x,y,z).
◊
Vecteur position d’un point Le vecteur position du point P dans le repère R(O,x,y,z), à l’instant t, est le vecteur OP(t). La trajectoire (T) du point P est l’ensemble des points P(t) obtenu lorsque t varie.
◊
Vecteur vitesse d’un point Le vecteur vitesse du point P par rapport au repère R(O,x,y,z), à l’instant t, est la dérivée du vecteur position OP(t) par rapport à t, dans R. r
d −−> OP (t ) dt R
V ( P / R ) =
Le vecteur vitesse du point P à l’instant t est tangent à la trajectoire en P(t) −− > r Si l’on considère un point A fixe dans R, et le point M(t) tel que : AM (t ) = V ( P / R )t , alors la trajectoire (H) du point M est appelée hodographe relatif au point A du vecteur vitesse de P par rapport au repère R.
◊
Vecteur accélération d’un point Le vecteur accélération du point P par rapport au repère R(O,x,y,z), à l’instant t, est la dérivée du vecteur vitesse V(P/R) par rapport à t, dans R. r d r Γ( P / R ) = V ( P / R ) dt R
Le vecteur accélération du point P par rapport au repère R à l’instant t est tangent à l’hodographe (H) du vecteur vitesse du point P au point M(t).
2
◊
LE SOLIDE INDEFORMABLE
Définition Un solide est dit indéformable lorsque quels que soient les points A et B de ce solide, la distance AB reste constante au cours du mouvement. On se limitera par la suite à appeler « solide » un solide indéformable.
20
Préparation à l’agrégation l’agrégation de mécanique mécanique : Polycopié Polycopié de Mécanique Mécanique des Solides Indéformab Indéformables les CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
◊
Equivalence Repère-Solide Dans un repère, la position relative des axes est invariante au cours du temps. C’est pourquoi un repère est équivalent à un solide. L’étude du mouvement du solide S2 par rapport au solide S1 est identique à l’étude du mouvement du repère R2 attaché au solide S2 par rapport au repère R1 attaché au solide S1.
3
PARAMETRAGE DE LA POSITION RELATIVE DE DEUX SOLIDES
Positionner le solide 2 par rapport au solide 1 revient donc à positionner le repère R2(O2,x2,y2,z2) attaché au solide 2, par rapport au repère R1(O1,x1,y1,z1) par rapport au solide 1. La position du repère R2(O2,x2,y2,z2) par rapport à R1(O1,x1,y1,z1) est complètement déterminée si l’on se fixe : •
•
3.1
les coordonnées de l’origine O2 du repère R2(O2,x2,y2,z2) dans le repère R1(O1,x1,y1,z1). Il existe plusieurs façons de définir ces coordonnées (cartésienne, cylindrique et sphérique). Dans tous les cas 3 paramètres indépendants sont nécessaires pour positionner O2 dans R1(O1,x1,y1,z1). l’orientation de la base R2(O2,x2,y2,z2) par rapport à R1(O1,x1,y1,z1). La base R2(O2,x2,y2,z2) est définie par 2 vecteurs (6 paramètres) unitaires (2 équations) et orthogonaux (1 équation), le troisième vecteur se déduisant des deux autres par un produit vectoriel. Trois paramètres indépendants sont également nécessaires pour positionner l’orientation de la base (x2,y2,z2) par rapport à la base (x1,y1,z1).
Définition des coordonnées de l’origine d’un repère. On cherche en premier lieu à positionner l’origine O2 du repère R2(O2,x2,y2,z2), dans un repère R1(O1,x,y,z). Il existe trois systèmes de coordonnées classiques, les coordonnées, cartésiennes, cylindriques et sphériques. n.b. L’origine O2 du repère R2, est choisie de façon arbitraire.
◊
Coordonnées cartésiennes Les coordonnées cartésiennes, notées ( x, x y, , y z ,z ) du point O2 dans le repère R1(O1,x,y,z) sont les projections du vecteur O1O2 sur chacun des axes (x,y,z).
→ → r r r r O1O2 = x x + y y + z z x = O1O2 ⋅ x
→ r y = O1O2 ⋅ y
→ r z = O1O2 ⋅ z
En coordonnées cartésiennes, la vitesse du point O2, par rapport au repère R1(O1,x,y,z) s’écrit:
→
dO1O2 dt
r
r
r
= dx x + dy y + dz z dt
dt
dt
21
Préparation à l’agrégation l’agrégation de mécanique mécanique : Polycopié Polycopié de Mécanique Mécanique des Solides Indéformab Indéformables les CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
◊
Coordonnées cylindriques Pour définir les coordonnées cylindriques, il faut d’abord définir la projection H du point O2 dans le plan (O1,x,y), puis un vecteur unitaire u de direction O1H. Les coordonnées cylindriques du point O2 sont alors ______
r = = O1 H , la mesure algébrique de la distance de O1 à H r r
θ = ( x , u ) , angle orienté dans le plan de normale z z, projection de O1O2 sur l’axe z.
Relation entre les coordonnées cartésiennes et cylindriques : x = r . cos θ et y = r . sin θ En coordonnées cylindriques, la vitesse du point O2, par rapport au repère R1(O1,x,y,z) s’écrit: → dO1O 2 dt
◊
=
r r r dr r d θ r dz r u + r . v + z , où v = z ∧ u dt dt dt
Coordonnées sphériques Pour définir les coordonnées sphériques, il faut d’abord définir la projection H du point O2 dans le plan (O1,x,y), puis un vecteur unitaire u de direction O1H, et enfin un vecteur unitaire w de direction O1O2. Les coordonnées sphériques du point O2 sont alors ______
ρ = O1O2 , la mesure algébrique de la distance de O1 à O2 r r
θ = ( x , u ) , angle orienté dans le plan de normale z r r
r
φ = ( z , w) , angle orienté dans le plan de normale v, où v
22
r
r
= z ∧ u
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
Relation entre les coordonnées cartésiennes et sphériques :
x = ρ . sin φ . cosθ , y
= ρ . sin φ . sin θ et z = ρ . cos φ
En coordonnées sphériques, la vitesse du point O2, par rapport au repère R1(01,x,y,z) s’écrit :
→
dO1O2 dt
3.2
=
d ρ r d φ r d θ r w + ρ .Sinφ . v − ρ .Sinφ . z dt dt dt
Définition de l’orientation relative de deux bases Dans un deuxième temps on cherche à définir l’orientation de la base (X,Y,Z) du repère R2(O2,X,Y,Y) par rapport à la base (x1,y1,z1) du repère R1(O1,x1,y1,z1). Pour cela plusieurs méthodes existent. La plus communément utilisée est le paramétrage par les angles d’Euler. Mais on trouvera aussi les angles de Cardan dit aussi de « roulis, tangage, lacet ». Ces méthodes s’appliquent bien à de petites variations d’angles autour d’une position connue. Lorsque de grandes rotations apparaissent on utilise alors un paramétrage par la méthode des Quaternions.
◊
Angles d’Euler Les trois angles d’Euler correspondent à la composition de trois rotations planes successives qui permettent de faire coïncider la base (x1,y1,z1) avec la base (X,Y,Z). La première rotation d’angle ψ, autour de l’axe z1 permet de passer à une première base intermédiaire (x2,y2,z2=z1). L’angle ψ est appelé « angle de précession ». Une seconde rotation d’angle θ, est alors appliquée autour de l’axe x2, de la première base intermédiaire, ce qui permet de définir une seconde base intermédiaire (x3=x2,y3,z3). L’angle θ est appelé « angle de nutation ». Une dernière rotation d’angle β est appliquée autour de l’axe z3 de la seconde base intermédiaire, ce qui permet de positionner la base (X,Y,Z=z3). L’angle β est appelé « angle de rotation propre ».
23
Préparation à l’agrégation l’agrégation de mécanique mécanique : Polycopié Polycopié de Mécanique Mécanique des Solides Indéformab Indéformables les CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
La composition de rotations planes successives permet de dessiner des figures de projection qui sont souvent très utiles pour la résolution des problèmes.
Figures de projection correspondant correspondant aux trois angles d’Euler d’Euler
◊
Angles de Cardan Les trois angles de cardan, ou « roulis, tangage, lacet » correspondent à la composition de trois rotations planes successives qui permettent de faire coïncider la base (x1,y1,z1) avec la base (x4,y4,z4). La première rotation d’angle θ, autour de l’axe x1 permet de passer à une première base intermédiaire (x2=x1,y2,z2=z1). L’angle θ est appelé « angle de roulis ». Une seconde rotation d’angle ψ, est alors appliquée autour de l’axe y2, de la première base intermédiaire, ce qui permet de définir une seconde base intermédiaire (x3,y3=y2,z3). L’angle ψ est appelé « angle de tangage ». 24
Préparation à l’agrégation l’agrégation de mécanique mécanique : Polycopié Polycopié de Mécanique Mécanique des Solides Indéformab Indéformables les CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
Une dernière rotation d’angle µ est appliquée autour de l’axe z3 de la seconde base intermédiaire, ce qui permet de positionner la base (x4,y4,z4=z3). L’angle µ est appelé « angle de lacet ».
Figures de projection correspondant correspondant aux trois angles de Cardan Ce paramétrage est habituellement employé pour paramétrer de petits mouvements du solide autour d’une base (x1,y1,z1) définie à l’aide de la trajectoire du centre de gravité du solide dans le référentiel du mouvement R. La direction x1, axe de roulis, est confondue avec la direction du vecteur vitesse du point O1 par rapport au référentiel R. La direction y1, axe de tangage, est orthogonale à x1, dans le plan local défini par la trajectoire du point O1. La direction z1, axe de lacet, est orthogonale au plan local défini par la trajectoire trajectoire du point O1.
◊ -
Quaternions Position du Problème Les angles d’Euler et de Cardan conduisent à des difficultés lors de la résolution numérique des problèmes de mécanique dans lesquels apparaissent apparaissent de grandes rotations. En effet, la composition de trois rotations planes successives permet de construire la matrice de passage R(t) entre le repère attaché au 25
Préparation à l’agrégation l’agrégation de mécanique mécanique : Polycopié Polycopié de Mécanique Mécanique des Solides Indéformab Indéformables les CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
solide et le repère attaché au référentiel du mouvement. Les éléments de cette matrice sont des fonctions non linéaires des trois angles choisis. Ceci complique l’écriture d’algorithmes de résolution numérique des problèmes. En effet, très schématiquement une résolution numérique est conduite de la façon suivante. A un instant t donné, donné, on suppose connus les trois angles d’Euler, l’orientation du solide est alors définie par une matrice de passage R(t ) fonction des trois angles d’Euler ou de Cardan. La résolution des équations du problème à l’instant t , permet de déterminer la vitesse de rotation instantanée du solide dR dt par rapport au référentiel référentiel du mouvement.
Aussi, l’orientation du solide à t+∆ t+∆t se se calcule t’elle comme suit: R(t + ∆t ) = I + ∆t .
dR . R dt
Pour pouvoir poursuivre le calcul numérique sur l’intervalle de temps suivant, la vitesse de rotation dans la nouvelle orientation doit alors être calculée. Si les équations du problème sont écrites en fonction des trois angles d’Euler, il faut donc déterminer les trois angles d’Euler paramétrant la nouvelle orientation. Attention ! Comme la matrice de passage contient des fonctions non linéaires des angles d’Euler, si ∆t est grand : θ (t + ∆t ) ≠ θ (t ) + ∆t .
d θ dt t
., quel que soit l’angle d’Euler θ considéré
En conséquence, soit le calcul est effectué à partir d’une succession de rotations infinitésimales (long !), soit les paramètres définissant l’orientation du solide à (t+∆t) doivent être déterminé en inverse à chaque itération. La matrice de passage qui appliquée à un vecteur défini dans le repère attaché au solide donne ses composantes dans le référentiel du mouvement pour un paramétrage en angles d’Euler (ψ,θ,β), s’écrit de la façon suivante :
Connaissant P(t), le calcul inverse des angles d’Euler est difficile, en particulier du fait d’indéterminations sur les angles (une même orientation peut être paramétrée par au minimum huit jeux d’angles d’Euler différents). Une nouvelle méthode de paramétrage du mouvement a donc été établie, dite méthode des quaternions, qui est bien adaptée à la résolution numérique des problèmes de mécanique des solides en grandes rotations.
-
Principe Une rotation Ω peut être simplement décrite à l’aide d’un axe de rotation (vecteur unitaire u) et d’un angle de rotation θ . On peut donc imaginer d’utiliser simplement comme paramètres du mouvement les coordonnées du vecteur unitaire u et l’angle θ. Cependant, avec ce choix de paramètres une indétermination demeure car une même rotation peut être obtenue avec un vecteur unitaireu unitaire u et un angle θ ou avec un vecteur unitaire – u et un angle 2π-θ. Cette indétermination peut être levée en utilisant comme paramètres le produit des coordonnées (a,b,c) du vecteur unitaire u avec le sinus de l’angle de rotation θ divisé par deux. Le paramétrage est alors mis sous la forme d’une matrice carrée Q d’ordre 4 appelée quaternion:
Si : {q 0 , q1 , q 2 , q 3 } = Cos Alors : Q
θ 2
, a.Sin
θ 2
, b.Sin
θ 2
, c.Sin
θ
2
= qo .e + q1 .i + q 2 . j + q3 .k
1 0 Avec : e= 0 0
0
0
1
0
0 1 0
0
0 − 1 0 , 1 0 i= 0 0 0 1 0 0 0
0 0 0
−1
0
0 0 , 0 j = 1 1 0 0
26
0
−1
0
0
0
0
1
0
0 − 1 et 0 k = 0 0 0 − 1 0
0
0
1
0
−1
0
1
0
0
0
0
0
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
Ces quatre matrices permettent de définir une base (e,i,j,k) d’un sous espace vectoriel dont les r éléments sont les quaternions Q, qui se décomposent en une partie réelle Re(Q) et une partie pure Q assimilée à un vecteur sur la base (i,j,k).
r
Re(Q ) = qo .e et
Si Q = qo .e + q1.i + q2 . j + q3 .k alors
Q
Des opérations de base ont été définies dans cette base. Avec Q
Q
= qo .e + q1 .i + q 2 . j + q3 .k
= qo .e − q1 .i − q2 . j − q3 .k
•
Conjugaison :
•
Produit : vu les matrices (e,i,j,k) on peut construire la table des produits de base :
i×i
= j × j = k × k = −e i × j = − j × i = k j × k = −k × j = i k × i = −i × k = j
e×e = e e×i
= i×e = i e × j = j × e = j e × k = k × e = k •
Produit vectoriel de parties pures : Compte tenu de l’analogie des trois derniers produits de base avec le produit vectoriel, on construit un produit vectoriel de parties pures par analogie avec le 3 produit vectoriel dans l’espace R : Soit
Q = qo .e + q1.i + q2 . j + q3 .k et Q * r
alors •
r
Q ∧ Q*
Produit scalaire de parties pures : même remarque
r r
= (q1q1* ).i + (q2 q2* ). j + (q3 q3* ).k r
On peut alors en déduire que : •
= qo * .e + q1* .i + q 2 * . j + q3* .k
= (q 2 q3* − q3 q2 * ).i + (q3 q1* − q1q3* ). j + (q1q2 * − q2 q1* ).k Q ⋅ Q*
r
Q ⋅ Q*
= Re Q.Q * )
Avec ces notations le produit de deux quaternions peut aussi s’écrire :
Q × Q* -
= q1 .i + q 2 . j + q3 .k
r
r
r
r
= qo .Q * + qo * .Q − (Q ⋅ Q * )e + (Q ∧ Q * )
Propriétés •
Le quaternion associé à la composition de rotations successives est le produit de leurs quaternions, dans l’ordre de leur application.
Q = Q1 × Q2 × Q3 •
r
=
1 q12
+ q22 + q32
q1 q2 q 3
ϕ Cos = q o 2 ϕ Sin ϕ = q12 + q 22 + q 32 2
et
La dérivée temporelle d’un quaternion s’effectue comme celle d’un vecteur
dQ dt •
Pour trois rotations, dans l’ordre 1,2,3
L’expression du vecteur rotation Ω à l’aide du vecteur unitaire u définissant l’axe de la rotation complète et l’angle de rotation j se calcule comme suit à partir du quaternion:
u
•
= {q0 , q1 , q2 , q3 }
=
dq o dt
.e +
dq1 dt
.i +
dq 2 dt
. j +
dq3 dt
.k
Le vecteur vitesse de rotation Ω d’une base par rapport au référentiel du mouvement se calcule de la façon suivante si Q est le quaternion associé à l’orientation de la base à l’instantt :
27
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
Ω=2
dQ dt
×Q
Exemple : Rotation d’angle θ, et d’axe z, dans un repère R(0,x,y,z). A l’instant t , l’orientation de la base est définie par : •
Q
θ θ θ θ = Cos , 0, 0, Sin = Cos .e + Sin .k 2 2 2 2 Alors :
4
4.1
dQ dt
θ &
θ θ = − Sin .e + Cos .k 2 2 2
et donc
et
2
Q
dQ dt
Q
= Cos
θ 2
.e − Sin
θ 2
.k
= θ &.k
CINEMATIQUE DU SOLIDE
Introduction, notations Soit un point P1 d’un solide S1 en mouvement par rapport au repère R(O,x,y,z). Les vecteurs vitesse et accélération du point P1 par rapport au repère R sont alors notés : r
V ( P 1 ∈ S 1 / R )
et
r
Γ( P 1 ∈ S 1 / R )
Cette notation permet de distinguer la vitesse d’un point appartenant à un solide de la vitesse d’un point de l’espace n’appartenant à aucun solide, comme par exemple le point de contact P entre les solides S1 et S2. La vitesse du point P sera alors notée : r
V ( P / R )
4.2
Champ des vecteurs vitesse des points d’un solide : torseur cinématique La formule de changement de base de dérivation, permet de définir la vitesse d’un point d’un solide, par rapport au référentiel du mouvement. Supposons un référentiel du mouvement R1(O1,x1,y1,z1) et un solide S2 en mouvement par rapport à ce référentiel, auquel est attaché un repère R2(O2,x2,y2,z2). La base attachée à R2 a une vitesse de rotation Ω par rapport à la base attachée à R1. Supposons deux points quelconques A et B du solide S2, alors :
28
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
→ → → → → d O A d O O d O A d O O d O A → r 1 = 1 2 + 2 = 1 2 + 2 + Ω ( R2 / R1) ∧ O2 A dt dt dt dt dt R1 R1 R1 R1 R 2 De même :
→ → → → → d O B d O O d O B d O O d O B → r 1 = 1 2 + 2 = 1 2 + 2 + Ω ( R2 / R1) ∧ O2 B dt dt dt dt dt R1 R1 R1 R1 R 2 Comme le solide S2 est indéformable :
→ → → → d O B d O A d O B d O A → r 2 = 2 = 0 ⇒ 1 = 1 + BA∧ Ω ( S 2 / R1) dt dt dt dt R 2 R 2 R1 R1
∀ A, ∀ B
Soit :
∀ A, ∀ B ∈ S 2
r
→
r
r
V ( B ∈ S 2 / R1) = V ( A ∈ S 2 / R1) + BA∧ Ω ( S 2 / R1 )
Le champ des vecteurs vitesses des points du solide S2 par rapport à R1(O1,x1,y1,z1), peut donc être représenté par un torseur, dit torseur cinématique, dont la résultante est la vitesse de rotation de la base (x2,y2,z2) par rapport à la base (x1,y1,z1), et le moment en un point A, la vitesse du point A appartenant au solide 2, par rapport au repère R1 : r
Ω( S 2 / R1) A {V ( S 2 / R1)}= r V ( A ∈ S 2 / R1) A
◊
Propriété : Equiprojectivité. Comme les champs des vecteurs vitesses d’un solide se représente par un torseur, c’est également un champ équiprojectif. Ceci signifie que quels que soient deux points A et B d’un solide S2 :
∀ A, ∀ B ∈ S 2
−−> r −−> r AB⋅ V ( A ∈ S 2 / R1 ) = AB⋅ V ( B ∈ S 2 / R1 )
On peut aussi le montrer directement, à partir de la propriété d’indéformabilité du solide : Si S2 est un solide indéformable, quels que soient A et B deux points appartenant à S2 alors la distance AB reste constante au cours du temps. Ceci s’écrit :
d AB ⋅ AB ) = 0 = AB ( dt dt d
⋅ AB
R1
Ceci s’écrit aussi : r r d d ( ) ( A ∈ S 2 / R1 )]⋅ AB = 0 − ⋅ = ∈ − O B O A AB [ V B S / R V 1 1 2 1 dt dt R1
Dont on déduit :
∀ A, ∀ B ∈ S 2
−− > r −−> r AB⋅ V ( A ∈ S 2 / R1 ) = AB⋅ V ( B ∈ S 2 / R1 )
Le champ des vitesses des points d’un solide indéformable est donc bien équiprojectif.
◊
Calcul du vecteur rotation instantanée. D’après ce qui a été dit plus haut, la dérivée des vecteurs unitaires (x2,y2,z2) de la base du repère R2(O2,x2,y2,z2) par rapport au référentiel R1, se calcule comme suit : 29
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
→ → d x2 r = Ω ( S 2 / R1) ∧ x2 dt R1 On en déduit alors, si l’on connaît les dérivées des vecteurs unitaires :
→ → → → d x → d y → d z 2 2Ω ( S 2 / R1) = x2 ∧ + y2 ∧ 2 + z 2 ∧ 2 dt dt dt R1 R1 R1 r
◊
Exemple : Mouvement de translation Si le mouvement du solide S2 par rapport à R1 se représente par un torseur couple : r
0 { V ( S 2 / R1)}= r A V ( A ∈ S 2 / R1) A Alors : ∀ B ∈ S 2
r
r
V ( B ∈ S 2 / R1) = V ( A ∈ S 2 / R1)
Le solide S2 est donc en translation par rapport à R1.
◊
Exemple : Mouvement de rotation instantanée Si le mouvement du solide S2 par rapport à R1 se représente au point A, par un torseur résultante : r
Ω(S 2 / R1) r A {V ( S 2 / R1)}= 0 A r
→
r
Alors : ∀ B ∈ S 2 V ( B ∈ S 2 / R1) = BA∧ Ω ( S 2 / R1 ) Le solide S2 est donc en rotation par rapport à R1 autour de l’axe central, (∆) du torseur cinématique. Cet axe central passe par le point A et sa direction est alignée avec Ω.
4.3
Champ des vecteurs accélération des points d’un solide. La formule de changement de base de dérivation, permet également de définir le champ des vecteurs accélération des points d’un solide, par rapport au référentiel du mouvement. Supposons un référentiel du mouvement R1(O1,x1,y1,z1) et un solide S2 en mouvement par rapport à ce référentiel, auquel est attaché un repère R2(O2,x2,y2,z2). La base attachée à R2 a une vitesse de rotation Ω par rapport à la base attachée à R1. Supposons deux points quelconques A et B du solide S2, alors :
→ r r r ∀ A, ∀ B ∈ S 2 V ( B ∈ S 2 / R1) = V ( A ∈ S 2 / R1) + BA∧ Ω ( S 2 / R1 ) La relation entre les vecteurs accélération des points A et B du solide S2 dans son mouvement par rapport au repère R1, s’obtient en dérivant les deux membres de cette égalité par rapport à t dans le repère R1.
d → r d V r ( B S 2 / R1) d V r ( A S 2 / R1) BA ( S / R ) ∈ = ∈ + ∧ Ω dt dt 2 1 dt R1 R1 R1 Soit :
30
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
d → r r r Γ( B ∈ S 2 / R1) = Γ( A ∈ S 2 / R1) + BA∧ Ω ( S 2 / R1 ) dt R1 En développant :
→ d r d → r r r Γ( B ∈ S 2 / R1) = Γ( A ∈ S 2 / R1) + BA ∧ Ω ( S 2 / R1 ) + BA∧ Ω ( S 2 / R1 ) dt R1 dt R1 Pour calculer la dérivée de BA par rapport à t, on utilise la formule de changement de base de dérivation, soit :
→ r → d → d → BA = BA + Ω(S 2 / R1) ∧ BA = 0 + Ω(S 2 / R1) ∧ BA dt R1 dt S 2 D’où la relation cherchée :
→ r → d r r Γ( B ∈ S 2 / R1) = Γ( A ∈ S 2 / R1) + Ω(S 2 / R1) ∧ BA ∧ Ω (S 2 / R1 ) + BA∧ Ω (S 2 / R1 ) dt R1 r
r
Où encore :
→ r → r r r d r Γ( B ∈ S 2 / R1) = Γ( A ∈ S 2 / R1) + Ω (S 2 / R1 ) ∧ AB+ Ω ( S 2 / R1 ) ∧ Ω(S 2 / R1) ∧ AB dt R1 Ainsi, le champ des vecteurs accélérations des points d’un solide ne peut pas être représenté par un torseur, du fait de l’existence du dernier terme.
5
5.1
COMPOSITION DES MOUVEMENTS
Introduction Soit un point P, appartenant à un solide S2, en mouvement à la fois par rapport à un repère R1(O1,x1,y1,z1) et par rapport à un repère R(O,x,y,z). On va chercher la relation entre les vecteurs vitesses V(P/R1) et V(P/R), ainsi que la relation entre les vecteurs accélération Γ(P/R1) et Γ(P/R). Ces relations, dites « de composition du mouvement », sont particulièrement utiles lorsqu’on étudie des mécanismes dans lesquels les mouvements relatifs mutuels des pièces sont connus, mais pas la cinématique d’ensemble du mécanisme.
31
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
5.2
Composition des vecteurs vitesse On cherche à définir en premier lieu la relation entre V(P/R) et V(P/R1) r
d → OP dt R
V ( P / R ) =
◊
d → O1 P dt R1
r
V ( P / R1) =
et
Relation de composition des vecteurs vitesses Donc : r
d → d → OP = OO1 dt R dt R
d → + O1 P dt R
V ( P / R ) =
d → OO1 = V (O1 ∈ R1 / R ) dt R
or :
Par ailleurs :
Par suite :
→ → r d → d → O1 P = O1 P + Ω( R1 / R ) ∧ O1 P = V ( P / R1) + Ω( R1 / R) ∧ O1 P dt R dt R1
→ r r V ( P / R ) = V ( P / R1) + V (O1 ∈ R1 / R ) + Ω( R1 / R ) ∧ O1 P
Si maintenant on considère le point de R1 qui coïncide avec P à l’instant t : → r r V ( P ∈ R1 / R ) = V (O1 ∈ R1 / R ) + Ω( R1 / R ) ∧ O1 P Et donc : r
r
r
V ( P / R ) = V ( P / R1) + V ( P ∈ R1 / R )
◊
Définitions Dans le mouvement du point P par rapport aux deux repères R et R1, on appelle : r
•
Vecteur vitesse absolue : V ( P / R )
32
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
r
•
Vecteur vitesse relative : V ( P / R1)
•
Vecteur vitesse d’entraînement : V ( P ∈ R1 / R )
r
◊
Généralisation Soit un point P mobile par rapport à n repères Ri(i=1,n). On peut écrire successivement : r
r
r
V ( P / R n −1 ) = V ( P / R n ) + V ( P ∈ Rn / R n −1 )
Ainsi : r
n
r
V ( P / R1 ) = V ( P / Rn ) +
∑=
r
V ( P ∈ Ri / Ri −1 )
i 2
◊
Composition des torseurs cinématiques Il a déjà été montré en utilisant la formule de changement de base de dérivation que lors de la composition des mouvements par rapport à n repères, n
r
Ω( Rn R1) =
∑=
r
Ω( Ri Ri −1 )
i 2
Comme par ailleurs : r
r
V ( P / R1 ) = V ( P / Rn ) +
n
∑=
r
V ( P ∈ Ri / Ri −1 )
i 2
On peut donc écrire la relation de composition des torseurs cinématiques :
{
r
r
Avec, V ( Ri
Ω( Ri / Ri −1 ) / Ri −1 )} = r V ( P ∈ Ri / Ri −1 )
La relation de composition des torseurs cinématiques s’écrit : n
r
r
{V ( Rn / R1 )} = ∑ {V ( Ri / Ri −1 )} i =2
◊
Exemple r
L’objectif de cet exemple est d’illustrer la distinction entre le vecteur vitesse absolue d’un point I, r
V ( I / R ) et son vecteur vitesse d’entraînement V ( I ∈ S 1 / R ) par un solide S1.
Supposons deux roues de friction S1 et S2. S1 est en rotation autour de l’axe (O,z) et S2 autour de l’axe (A,z). On pose :
r
Ω(S 1 / R ) = ω 1 z
et
r
Ω(S 2 / R ) = ω 2 z
→ r Les deux roues de friction sont en contact au point I : OI = r 1 y
et
→ r AI = −r 2 y
Le vecteur vitesse absolue du point de contact I par rapport au repère R s’obtient en dérivant le vecteur position du point I : r
d → d r OI = r 1 y dt R dt R
V ( I / R ) =
Soit :
r
r
V ( I / R ) = 0
Le vecteur vitesse du point du solide S1 qui à l’instant t coïncide avec I, noté IœS1, est le vecteur vitesse d’un point qui décrit un cercle de centre O et de rayon r1 à la vitesse ω1.
33
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
Soit :
r
r
V ( I ∈ S 1 / R ) = −ω 1 r 1 x
Le vecteur vitesse du point du solide S2 qui à l’instant t coïncide avec I, noté IœS2, est le vecteur vitesse d’un point qui décrit un cercle de centre A et de rayon r2 à la vitesse ω2. Soit :
5.3
r
r
V ( I ∈ S 1 / R ) = ω 2 r 2 x
Composition des vecteurs accélération Soit un point P, appartenant à un solide S2, en mouvement à la fois par rapport à un repère R1(O1,x1,y1,z1) et par rapport à un repère R(O,x,y,z). On va chercher maintenant la relation entre les vecteurs accélération Γ(P/R1) et Γ(P/R).
◊
Relation de composition des vecteurs accélération Il a été montré aux paragraphes précédents que : r
r
r
V ( P / R ) = V ( P / R1) + V ( P ∈ R1 / R )
Ce qui s’écrit aussi : → r r V ( P / R ) = V ( P / R1) + V (O1 ∈ R1 / R ) + Ω( R1 / R ) ∧ O1 P Dérivons chaque terme par rapport au temps dans le repère R : → d d r d r d = + ∈ + Ω ∧ V P / R V P / R 1 V O R 1 / R R 1 / R O ( ) ( ) ( ) ( ) 1 P dt dt dt 1 dt R R R R Soit, en appliquant la formule de changement de base de dérivation au premier terme : r r r d r d r ( ) ( ) ( ) ( ) ( ) ( ) ( P / R1) = + Ω ∧ = Γ + Ω ∧ V P / R 1 V P / R 1 R 1 / R V P / R 1 P / R 1 R 1 / R V dt dt R R1
Comme O1 est l’origine du repère R1 :
d dt V (O1 ∈ R1 / R ) = Γ(O1 ∈ R1 / R ) R Enfin, en appliquant la formule de changement de base de dérivation au dernier terme :
34
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
d → → → Ω( R1 / R ) ∧ O1 P = d Ω( R1 / R ) ∧ O1 P + Ω( R1 / R ) ∧ d O1 P dt R dt R dt R d → → → Ω( R1 / R ) ∧ O1 P = d Ω( R1 / R ) ∧ O1 P + Ω( R1 / R ) ∧ V ( P / R1) + Ω( R1 / R) ∧ O1 P dt R R dt Comme, par ailleurs, d’après la formule de changement de point du champ de vecteurs accélérations:
→ r → r d r Γ( P ∈ R1 / R) = Γ(O1 ∈ R1 / R) + Ω ( R1 / R) ∧ O1 P + Ω ( R1 / R) ∧ Ω ( R1 / R) ∧ O1 P dt R r
r
On en déduit en regroupant tous les termes : r
r
r
r
Γ( P / R ) = Γ( P / R1) + Γ( P ∈ R1 / R) + 2.Ω ( R1 / R) ∧ V ( P / R1)
◊
Définitions Dans le mouvement du point P par rapport aux repères R et R1, on appelle : r
-
Vecteur accélération absolue :
Γ( P / R )
-
Vecteur accélération relative :
Γ( P / R1)
-
Vecteur accélération d’entraînement :
-
Vecteur accélération de Coriolis : 2.Ω ( R1 / R ) ∧ V ( P / R1)
r
r
Γ( P ∈ R1 / R)
r
6
A RETENIR
-
Indéformabilité équivalence repère/solide
-
Paramétrage complet de la position d’un solide par rapport à un référentiel Coordonnées de l’origine (3 paramètres) Orientation de la base (3 paramètres)
◊
Cinématique du point : r
d −−> OP (t ) dt R
-
Vecteur vitesse du point P, par rapport à R, à l’instant t : V ( P / R ) =
-
Vecteur accélération du point P, par rapport à R, à l’instant t :
◊
r d r Γ( P / R ) = V ( P / R ) dt R
Cinématique du solide : r
-
Vecteur vitesse du point P, appartenant au solide S1 par rapport à R, à l’instant t : V ( P ∈ S 1 / R )
35
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DU POINT, CINEMATIQUE DU SOLIDE
-
◊ -
Vecteur accélération du point P, appartenant au solide S1,par rapport à R, à l’instant t : Γ( P ∈ S 1 / R)
Formules de changement de point Formules de changement de point du champ des vecteurs vitesse d’un solide
∀ A, ∀ B ∈ S 1
r
→
r
r
V ( B ∈ S 1 / R) = V ( A ∈ S 1 / R ) + BA∧ Ω ( S 1 / R )
Le champ des vecteurs vitesses des points du solide S1 en mouvement par rapport à R, se représente par un torseur, dit torseur cinématique : r
Ω( S 1 / R ) A {V ( S 1/ R )}= r V ( A ∈ S 1 / R ) A -
Formules de changement de point du champ des vecteurs accélération d’un solide
→ r → r r r d r Γ( B ∈ S 1 / R ) = Γ( A ∈ S 1 / R ) + Ω ( S 1 / R ) ∧ AB + Ω ( S 1 / R ) ∧ Ω(S 1 / R ) ∧ AB dt R Le champ des vecteurs accélération des points du solide S1 en mouvement par rapport à R, ne se représente pas par un torseur.
◊
Formules de composition des mouvements. On suppose un point P en mouvement par rapport à un solide S1 auquel est attaché un repère R1, lui-même en mouvement par rapport au référentiel du mouvement R. On note (PœR1) le point du solide S1 qui à l’instant t, coïncide avec P alors : r
r
r
1 424 3
1 424 3
14 4 244 3
absolue
relative
V ( P / R ) = V ( P / R1) + V ( P ∈ R1 / R ) r
r
entraînement
r
r
Γ ( P / R ) = Γ ( P / R1) + Γ ( P ∈ R1 / R) + 2.Ω ( R1 / R) ∧ V ( P / R1) 1 424 3 1 424 3 14 4 244 3 144 4 4 244 4 4 3 absolue
relative
entraînement
Coriolis
36
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DES SYSTEMES DE SOLIDES
1
Système de solides _________________________________________________ 38
◊ ◊ ◊ ◊
Définitions______________________________________________________________ Types de liaisons _________________________________________________________ Représentation d’une liaison. _______________________________________________ Degrés de liberté d’une liaison. ______________________________________________
38 38 39 39
2
Tableau des liaisons normalisées _____________________________________ 39
3
Cinématique du contact entre deux solides _____________________________ 41
4
Modélisation cinématique ___________________________________________ 42 4.1
Introduction _______________________________________________________ 42
4.2
Graphe cinématique_________________________________________________ 42
4.3
Graphe de structure ou graphe des liaisons _____________________________ 42
4.4
Mobilité d’un système _______________________________________________ 42
◊ ◊ ◊ 4.5
Exemple : Presse de modélisme _______________________________________ 44
◊ ◊ ◊ ◊ ◊ 5
Fermeture géométrique ____________________________________________________ 43 Fermeture cinématique ____________________________________________________ 43 Calcul de la mobilité ______________________________________________________ 43 Plan du mécanisme _______________________________________________________ Construction du graphe de structure __________________________________________ Réduction du graphe de structure ____________________________________________ Construction du graphe cinématique__________________________________________ Mobilité du système ______________________________________________________
44 45 45 46 46
A retenir _________________________________________________________ 47
37
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DES SYSTEMES DE SOLIDES
1
◊
SYSTEME DE SOLIDES
Définitions Lorsque la mécanique du solide est appliquée à des mécanismes, les mouvements relatifs entre solides sont limités par l’existence de liaisons entre les différentes pièces du mécanisme. Ainsi, un système de solide est il constitué de deux sous-ensembles, l’ensemble des solides indéformables et l’ensemble des liaisons entre solides. Par la suite, le système de solides sera enrichi d’un troisième sous-ensemble, l’ensemble des actions mécaniques. Le système de solides pourra donc être représenté par des graphes, dont l’analyse permet de définir le nombre d’inconnues cinématiques du système. Par ailleurs, cette représentation permet d’aider au choix des sous-systèmes à isoler, des théorèmes généraux à appliquer et des projections pertinentes à effectuer. Graphe cinématique : les liaisons constituent les sommets, les solides constituent les arcs. Graphe de structure : les solides constituent les sommets, les liaisons constituent les arcs.
◊ -
Types de liaisons Liaison unilatérale/bilatérale Lorsqu’une liaison du fait même de sa réalisation technologique ne peut pas être rompue (sauf par destruction du système), elle est dite bilatérale. Dans le cas contraire, la liaison est dite unilatérale.
-
Liaison holonome/non holonome Les liaisons pour lesquelles l’équation de liaison est uniquement fonction des paramètres de position (équation holonome), est dite liaison cinématique. Sinon, la l’équation de liaison est dite nonholonome. Par exemple, une liaison pivot, autorise la rotation autour de l’axe du pivot, mais interdit les autres mouvements, translations, ou rotations autour des deux autres axes. Si les mouvements relatifs entre les deux solides S1 et S2 sont paramétrés par trois paramètres de translation et trois paramètres de rotation : Translations : X . x + Y . y + Z . z r
r
r
Rotations :
α . x + β . y + γ . z r
r
r
Alors l’existence d’une liaison pivot impose les cinq équations ho lonomes suivantes : X=Xo, Y=Yo, Z=Zo, β = β o, γ =γ o
Si cette liaison est motorisée et que le moteur impose une vitesse de rotation w(t), cette d α = ω (t ) . motorisation impose une dernière équation de liaison qui, cette fois, est non-holonome: dt En règle générale, un actionneur impose l’évolution temporelle d’un paramètre de position et conduit donc, quelles que soient les conditions de fonctionnement à des équations non-holonomes.. 38
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DES SYSTEMES DE SOLIDES
◊
Représentation d’une liaison. Par rapport au repère local attaché à une liaison entre deux pièces, le champ des vitesses relatives entre les deux pièces (S1 et S2) reliées par la liaison en question peut être représenté par un torseur. Le repère de la liaison est en général choisi de telle sorte que l’axe x, soit un axe central pour la liaison, aligné avec Ω. Alors en tout point A de cet axe :
Ω(S 2 / S 1) A {V ( S 2/S 1)}= V A r
r
◊
∀ A ∈ (O, x ) r
Degrés de liberté d’une liaison. Dans le repère local associé à la liaison entre deux solides, les mouvements relatifs des deux solides sont limités à trois translations et trois rotations au maximum. Parmi ces 6 mouvements élémentaires, le nombre de mouvements élémentaires indépendants autorisés par la liaison définit le degré de liberté de cette liaison. La liaison pivot ci-dessus est à un seul degré de liberté. Une liaison hélicoïdale (vis-écrou) permet deux mouvements, rotation autour de l’axe de la liaison et translation le long de ce même axe. Mais dans ce cas, rotation et translation sont proportionnelles. Par conséquent la liaison hélicoïdale est également à un seul degré de liberté.
2
TABLEAU DES LIAISONS NORMALISEES
Les liaisons normalisées présentées dans le tableau qui suit sont classées par degré de liberté croissant. La représentation graphique de la liaison ne présume pas de la réalisation technologique de ces liaisons, mais est une représentation schématique normalisée des mouvements autorisés par cette liaison entre les deux solides Sk et Si..
39
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DES SYSTEMES DE SOLIDES
40
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DES SYSTEMES DE SOLIDES
3
CINEMATIQUE DU CONTACT ENTRE DEUX SOLIDES
Dans le cas du contact ponctuel, une terminologie particulière est employée, qui est décrite cidessous. Pour un contact surfacique ou linéique, tout ce qui est écrit pour le contact ponctuel reste valable en chaque point de la surface ou de la ligne de contact.
Supposons deux solides S1 et S2 en contact en un point I. Il existe un unique plan tangent Π entre les deux solides, défini par la normale n21, à S2 ou à S1, au point de contact I, dirigée de S2 vers S1. Alors les éléments de réduction du torseur cinématique en I du mouvement de S1 par rapport à S2 sont :
Ω( S 1 / S 2) I {V ( S 1/ R 2 )}= V ( I ∈ S 1 / S 2) I r
r
On note alors :
Ω( S 1 / S 2) = (Ω( S 1 / S 2) ⋅ n 21 )n 21 + n 21 ∧ (Ω( S 1 / S 2) ∧ n 21 ) r
r
r
r
r
r
r
144 244 3
144 4 4 244 4 4 3
Pivotement
Roulement
La projection du vecteur rotation sur la normale au plan de contact, est le terme depivotement. La composante du vecteur rotation appartement au plan tangent est le terme de roulement . r
Enfin, le moment en I du torseur V ( I ∈ S 1 / S 2) est le glissement en I de S1 par rapport à S2. On r
le note souvent aussi : G ( I ∈ S 1 / S 2) . D’après la formule de composition des vitesses : r
r
r
V ( I / S 2 ) = V ( I / S 1) + V ( I ∈ S 1 / S 2) 1 4 24 3
1 4 24 3
14 4 244 3
absolue
relative
entraînement
Ce qui donne donc ici : r
r
r
G ( I ∈ S 1 / S 2) = V ( I / S 2 ) − V ( I / S 1)
Ainsi le glissement en I de S1 par rapport à S2 est un vecteur parallèle au plan de contact.
41
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DES SYSTEMES DE SOLIDES
4
4.1
MODELISATION CINEMATIQUE
Introduction Un système de solides indéformables est composé de deux sous-ensembles : les solides et les liaisons entre les solides. Une représentation peut donc en être faite au moyen d’un graphe. Un graphe est aussi composé de deux sous-ensembles : des points appelés sommets du graphe et des lignes appelés arcs qui relient certains sommets entre eux. Il y a donc deux bijections possibles entre un système de solides indéformables et un graphe.
4.2
Graphe cinématique Lorsqu’on représente un système de solides par un graphe cinématique, les sommets du graphe représentent les liaisons et les arcs, les solides. Les sommets représentatifs des liaisons sont dessinés en respectant la normalisation et les positions spatiales relatives des entités géométriques caractéristiques. Un graphe cinématique est donc en général (le cas particulier est relatif à un problème plan) un graphe tridimensionnel et se dessine en perspective. La fonction principale du graphe cinématique est d’aider à la compréhension du fonctionnement du système, à la visualisation du paramétrage et au calcul.
4.3
Graphe de structure ou graphe des liaisons A l’inverse, lorsqu’on représente un système de solides par un graphe de structure, les sommets du graphe représentent les solides et les arcs, les liaisons. Ce graphe pourra être complété par la suite par des arcs parallèles figurant les actions mécaniques. Sur chaque arc, il y a le nom de la liaison qu’il représente ainsi que les caractéristiques géométriques. Aux sommets sont placés les symboles alphanumériques désignant les solides. Le graphe de structure a deux fonctions principales : - aider à la détermination de la mobilité du système c’est à dire du nombre minimal de paramètres permettant de décrire complètement la cinématique du système, - aider au choix des sous-systèmes à isoler, des théorèmes généraux de la dynamique à utiliser, des projections à effectuer pour répondre à un problème posé.
4.4
Mobilité d’un système La mobilité d’un système correspond au nombre minimal de paramètres indépendants nécessaires pour décrire totalement la cinématique du système. Dans un mécanisme, chaque liaison présente un certains nombre de degrés de liberté. Mais la mobilité du système complet n’est pas égale à la somme des degrés de liberté de chacune des liaisons. Le graphe de structure sera généralement employé pour déterminer la mobilité du système et choisir les paramètres indépendants du problème. En effet, lorsque le graphe présente des fermetures, des équations supplémentaires entre les paramètres apparaissent, ce qui diminue d’autant la mobilité du système.
42
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DES SYSTEMES DE SOLIDES
◊
Fermeture géométrique Lorsque dans le graphe de structure apparaît un chemin fermé, (S1, S2, … ; Sn-1, Sn, S1) alors, la fermeture géométrique de ce chemin s’écrit : O1O 2
+
....
r
+ On − 1On + On O1 = 0 et
P (b1 / b n ) ⋅ P (b1 / bn −1 ) ⋅ ... ⋅ P (b2 / b1 ) = Identité
où (Oi,bi) est le repère, d’origine Oi et de base bi, attaché à chaque solide Si, et où P(bi+1,bi) est une matrice de changement de base Les équations scalaires obtenues sont des équations holonomes.
◊
Fermeture cinématique Si le chemin fermé possède des liaisons cinématiques, il faut alors écrire une équation de fermeture cinématique portant sur le torseur cinématique du chemin fermé : {V(Sn/S1)}+{V(S1/S2)}+ … + {V(Sn-1/Sn)}=0
Les équations scalaires obtenues sont des équations non-holonomes. Il est toujours possible d’écrire une fermeture cinématique à la place d’une fermeture géométrique. Les équations non holonomes de la fermeture cinématique forment un système équivalent à celui obtenu par dérivation temporelle des équations holonomes de la fermeture géométrique. Le choix d’utiliser une fermeture géométrique ou une fermeture cinématique sera guidé par des conditions de simplicité de mise en œuvre et conduira souvent à une procédure mixte.
◊
Calcul de la mobilité C’est le nombre minimal de paramètres indépendants nécessaires pour décrire totalement la cinématique du système. Pour l’obtenir, il faut suivre la procédure suivante : En premier lieu, déterminer le nombre maximal de chemins fermés indépendants. Ce nombre s’appelle le nombre cyclomatique et vaut : µ = nl - ns + 1, avec nl=nombre de liaisons et ns=nombre de sommets du graphe. En second lieu, pour chacun de ces µ chemins fermés, il faut écrire les équations de fermeture et déterminer le rang r du système d’équations obtenu. Si on note n p, nombre total de paramètres de position, qui est égal à la somme des degrés de liberté de toutes les liaisons du système, la mobilité m du système est alors p ar définition : m = n p - r Cette procédure est systématique, mais généralement fastidieuse. Dans la pratique, il n’est pas nécessaire de paramétrer explicitement tous les degrés de liberté de toutes les liaisons et d’expliciter ensuite toutes les équations de fermeture. On peut souvent remplacer les chemins fermés du graphe de structure par une liaison équivalente (voir l’exemple ci-dessous). Ceci permet de réduire le graphe de structure et d’en déduire un graphe de structure minimal et un paramétrage minimal. On peut alors calculer la mobilité en appliquant la procédure décrite ici au graphe de structure minimal ou bien la déterminer en imaginant le blocage d’un degré de liberté d’une liaison. On regarde si le système reste mobile ou non. S’il reste mobile, on ajoute un deuxième blocage d’un nouveau degré de liberté et ainsi de suite jusqu’à immobilité complète du système. La mobilité est alors le nombre de blocages effectués.
-
Mobilité utile et mobilité interne On peut classer les n p paramètres de position en deux catégories suivant qu’ils sont associés à des liaisons avec l’extérieur du système ou à des liaisons internes au système. Par commodité, on parlera de paramètres utiles et de paramètres internes. La mobilité utile peut être trouvée en utilisant la « procédure » du blocage. On observe le système sous la forme d’une boîte noire dont les seuls degrés de liberté observables et accessibles (donc blocables) sont ceux des liaisons externes.
43
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DES SYSTEMES DE SOLIDES
-
Bilan : Choix d’un paramétrage • • • • • • • •
4.5
Identifier les solides, identifier les liaisons. Tracer le graphe de structure complet. Calculer le nombre maximal de chemins fermés indépendants. Réduire le graphe de structure en remplaçant autant de chemins fermés indépendants que possible par une liaison équivalente. Tracer le graphe de structure minimal. Choisir un paramétrage minimal associé au graphe de structure minimal. Tracer les figures de projection associées au paramétrage choisi. Expliciter les équations de fermeture restantes.
Exemple : Presse de modélisme Le cas d’une presse de modélisme est présenté ici pour illustrer les principes de la modélisation cinématique qui ont été évoqués plus haut. Le plan du mécanisme est présenté ci-dessous. Cet exemple est issu de Mecanique 1, Yves Brémont/Paul Réocreux.
◊
Plan du mécanisme
La presse est constituée d’un bâti, constitué d’une embase 00, d’un plan d’appui 05, de deux colonnes 01 et 02 et d’une bague supérieure 04. Ces pièces n’ont aucun mouvement relatif. L’ensemble de ces pièces sera donc noté par (0).
44
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DES SYSTEMES DE SOLIDES
Par ailleurs, la traverse 10, les deux bagues 11, et les deux rondelles 12 et 13 n’ont également aucun mouvement relatif. L’ensemble de ces pièces sera noté (1).
◊
Construction du graphe de structure Nous avons donc 6 solides principaux et 8 liaisons. Ce qui permet de dessiner le graphe de structure et de choisir les paramètres du mouvement pour chacune des liaisons. On tient compte du fait que le problème est plan.
Dans ce graphe apparaissent des chemins fermés. On peut calculer l e nombre maximal de chemins fermés indépendant comme suit : Nombre de liaisons nl=8 Nombre de sommets ns=6 Nombre de chemins fermés indépendants ou nombre cyclomatique µ = nl -ns +1=3. Ces chemins fermés permettent d’écrire des équations de fermeture et donc de réduire le nombre de paramètres nécessaire à la modélisation complète de la cinématique du système. On peut paramétrer chacune des liaisons puis poser les équations et réduire le nombre de paramètres, ou analyser le problème et remplacer les chemins fermés par des liaisons équivalentes
◊
Réduction du graphe de structure L’existence de deux liaisons pivots parallèles entre le bâti (0) et la traverse (10) interdit la rotation autour de ces axes. Ainsi ces deux liaisons parallèles peuvent elles être remplacées par une liaison glissière. On élimine ainsi un premier chemin fermé. Ensuite, la tige 20 est liée à la traverse 10 par deux branches parallèles. Dans chaque branche on trouve une liaison rotule de centre O (40/10) ou (30/10) puis une liaison appui plan (20/30) ou (20/40). La mise en série d’une rotule et d’un appui plan est équivalente à une liaison ponctuelle. Deux liaisons ponctuelles au même point, équivalent à une seule. Ainsi le second chemin fermé est-il ramené à une unique liaison ponctuelle. Enfin, le dernier chemin fermé est naturellement réduit en considérant que la traverse sommet (03), encastrée au bâti, fait partie du bâti. On peut alors dessiner un g raphe de structure simplifié, pour lequel est aussi choisi un paramétrage. Il reste encore un chemin fermé donc des équations de fermeture à poser.
45
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DES SYSTEMES DE SOLIDES
◊
Construction du graphe cinématique La modélisation cinématique retenue peut être également représentée à l’aide du schéma cinématique, pour lequel on place les liaisons aux sommets et les solides sur les arcs du graphe. Ce graphe permet une meilleure compréhension du fonctionnement du système.
◊
Mobilité du système • • • • • • • •
Nombre total de paramètres n p=7 Nombre de liaisons nl=3 Nombre de sommets ns=3 Nombre de chemins fermés indépendants ou nombre cyclomatique µ = nl -ns +1=1 Nombre d’équations scalaires de fermeture à écrire n=6 : vecteur translation et vecteur rotation projetés sur les axes x,y,z. Mobilité m=np-n=1 Mobilité utile : Le paramètre d’entrée est ωz, paramètre de sortie az. Les deux paramètres sont liés, la mobilité utile est égale à un. Mobilité interne : La mobilité interne est alors égale à zéro. Les deux déplacements bx et by et les trois rotations rx, ry et rz sont bloquées.
46
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CINEMATIQUE DES SYSTEMES DE SOLIDES
5
-
A RETENIR
Système de solides ensemble {{solides}, {liaisons}}. Se modélise à l’aide d’un graphe. Graphe de structure : Sommets=solides, arcs=liaisons. Graphe cinématique : Sommets=liaisons, arcs=solides. La modélisation permet le calcul de la mobilité et le choix d’un paramétrage. Le choix d’un paramétrage permet la mise en place des figures de projections.
-
Liaison Se représente par un torseur cinématique Bilatérale/Unilatérale Holonome/Non Holonome
-
Notations spécifiques au contact Torseur cinématique au point de contact I entre les solides S1 et S2, où S1 est en mouvement par rapport à S2 et où n21 est la normale au plan de contact au point I entre S1 et S2 dirigée de 2 vers 1: •
Ω( S 1 / S 2) 1 2 { } ( / ) = S R V I V ( I ∈ S 1 / S 2) I
•
Pivotement :
r
•
r
r
Ω( S 1 / S 2) ⋅ n 21 Roulement : n 21 ∧ (Ω( S 1 / S 2) ∧ n 21 ) Glissement : G ( I ∈ S 1 / S 2) = V ( I ∈ S 1 / S 2) = V ( I / S 2 ) − V ( I / S 1) r
r
•
-
r
r
r
r
r
r
Bilan : Choix d’un paramétrage • • • • • • • •
Identifier les solides, identifier les liaisons. Tracer le graphe de structure complet. Calculer le nombre maximal de chemins fermés indépendants. Réduire le graphe de structure en remplaçant autant de chemins fermés indépendants que possible par une liaison équivalente. Tracer le graphe de structure minimal. Choisir un paramétrage minimal associé au graphe de structure minimal. Tracer les figures de projection associées au paramétrage choisi. Expliciter les équations de fermeture restantes.
47
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE LA MASSE ET CINETIQUE
1
Système matériel à masse conservative_________________________________ 49 ◊ Définitions______________________________________________________________ 49 ◊ Conséquences ___________________________________________________________ 49
2
Torseur cinétique, torseur dynamique et energie cinétique_________________ 50 ◊ Torseur cinétique_________________________________________________________ 50 ◊ Torseur dynamique _______________________________________________________ 50 ◊ Energie cinétique_________________________________________________________ 50 ◊ Autres cas. Exemple : Action de la pesanteur ___________________________________ 51
3
Centre d’inertie, opérateur d’inertie ___________________________________ 51
3.1
Centre d’inertie G __________________________________________________ 51
3.2
Opérateur d’inertie JA( ) ____________________________________________ 51
◊ ◊ ◊
3.3
Exemples__________________________________________________________ 53 ◊ ◊ ◊ ◊ ◊
4
Opérateur d’inertie en O d’un disque D de rayon R, de centre O et de masse m. ________ Opérateur d’inertie en O d’un cylindre C de rayon R, de hauteur h et de centre O. ______ Opérateur d’inertie en O d’un cône de révolution C de rayon R et de hauteur h. ________ Opérateur d’inertie en O d’une sphère creuse S de centre O, de rayon R et de masse m.__ Opérateur d’inertie en O d’une sphère pleine S de centre O, de rayon R et de masse m. __
53 53 54 54 54
Consequences du principe de conservation de la masse ___________________ 54
4.1
Torseur cinétique ___________________________________________________ 54 ◊ ◊
4.2
Expression de la résultante cinétique. _________________________________________ 54 Expression du moment cinétique. ____________________________________________ 55
Torseur dynamique _________________________________________________ 55 ◊ ◊
5
Définition ______________________________________________________________ 51 Relation entre JA(Σ) et JG(Σ) ou théorème de Huyghens généralisé . _________________ 51 Expression dans la base (O,x,y,z) ____________________________________________ 52
Expression de la résultante dynamique. _______________________________________ 55 Relation entre les moments cinétique et dynamique. _____________________________ 56
4.3
Énergie cinétique ___________________________________________________ 56
4.4
Autres cas. Exemple : Action de la pesanteur ____________________________ 57
A retenir _________________________________________________________ 58 ◊ Conséquence du principe de conservation de la masse ____________________________ 58 ◊ Torseur cinétique_________________________________________________________ 58 ◊ Torseur dynamique _______________________________________________________ 58 ◊ Energie cinétique_________________________________________________________ 58
48
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE LA MASSE ET CINETIQUE
1
◊
SYSTEME MATERIEL A MASSE CONSERVATIVE
Définitions •
•
Système matériel : Un système matériel est un système sur lequel est définie une mesure de la masse. La masse est habituellement définie à l’aide d’une densité volumique, surfacique ou linéique. Système matériel à masse conservative : Un système matériel D est dit à masse conservative si toute partie d de D, qu’on suit au cours du temps, a une masse constante :
∀d ∈ D ∀t
◊
m(d , t ) = constante = m(d )
Conséquences
Si l’on exprime la masse à l’aide d’une densité volumique :
∀d ∈ D
∫( )ρ ( P )dv( P ) = constante
m(d , t ) =
d t
Alors :
∀d ∈ D
d dt
m(d , t ) =
d dt ρ ( P ) + ρ ( P ).divV t ( P ) dv( P ) = 0 d (t )
∫
r
Ceci implique que : d
r
ρ ( P ) + ρ ( P ).divV t ( P ) = 0 dt
Ce qui est l’équation de continuité de la mécanique des milieux continus tridimensionnels. Plus généralement, si l’on intègre des fonctions régulières relativement à la densité volumique : H (t ) =
∫ h( P ). ρ ( P ).dv( P ) Dt
Alors : d
H (t ) = dt
d dt [h( P ). ρ ( P )] + h( P ). ρ ( P ).divV t ( P ).dv( P ) D
∫
r
t
Si l’on développe :
d d d d . . . . ( ) ( ) ( ) ( ) ( ) ( ) ( ) H t = ρ P h P + h P ρ P + ρ P divV t P dv( P ) = ρ ( P ). h( P ) .dv( P ) dt dt dt dt Dt Dt =0
∫
∫
r
144 4 4 244 4 4 3
Ainsi, pour un système matériel à masse conservative et pour toute fonction h régulière : H (t ) =
∫ Dt
h( P ). ρ ( P ).dv ( P ) ⇒
d
H (t ) = dt
49
∫ Dt
d h( P ) dt
.ρ ( P ).dv( P )
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE LA MASSE ET CINETIQUE
2
TORSEUR CINETIQUE, TORSEUR DYNAMIQUE ET ENERGIE CINETIQUE
Soit un système matériel Σ en mouvement par rapport à un repère R, différentes intégrales faisant intervenir la mesure de la masse peuvent être définies sur ce système matériel. On définira en particulier, les quantités de mouvements, les quantités d’accélérations et l’énergie cinétique. Le principe de conservation de la masse permettra ensuite de simplifier ces expressions et d’écrire des relations entre elles.
◊
Torseur cinétique
Le torseur cinétique, ou torseur des quantités de mouvement, du système matériel Σ dans son mouvement par rapport au repère R est, en un point A quelconque, le torseur à structure suivant :
V ( P R ).dm( P ) P ∈Σ {C (Σ R )}= → AP ∧ V ( P R ).dm( P ) A P ∈Σ
∫
r
∫
r
-
La résultante du torseur cinétique est appelée résultante cinétique ou encore quantité de mouvement deΣ dans son mouvement par rapport à R.
-
Le moment du torseur cinétique est appelé moment cinétique, au point A de Σ dans son mouvement par rapport à R. On le note habituellement : σ A (Σ R ) r
◊
Torseur dynamique
Le torseur dynamique, ou torseur des quantités d’accélération, du système matériel Σ dans son mouvement par rapport au repère R est, en un point A quelconque, le torseur à structure suivant :
Γ( P R ).dm( P ) P ∈Σ {D (Σ R )}= → AP ∧ Γ( P R ).dm( P ) A P ∈Σ
∫
r
∫
r
-
La résultante du torseur dynamique est appelée résultante dynamique ou encore quantité d’accélération de Σ dans son mouvement par rapport à R.
-
Le moment du torseur dynamique est appelé moment dynamique, au point A de Σ dans son mouvement par rapport à R. On le note habituellement : δ A (Σ R )
-
La formule de changement de point du moment dynamique est, par construction :
r
r
r
r
δ A (Σ / R) = δ B (Σ / R) + AB ∧ m(Σ).Γ( B / R)
◊
∀ A, B
Energie cinétique
L’énergie cinétique du système matériel Σ dans son mouvement par rapport à un repère R est le scalaire suivant : T (Σ R ) =
1 2
2 [ ] ( ) V P R dm( P ) ∫ r
P ∈Σ
50
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE LA MASSE ET CINETIQUE
◊
Autres cas. Exemple : Action de la pesanteur
Il existe d’autres cas où une quantité faisant intervenir la mesure de la masse est intégrée sur un solide. L’action de la pesanteur sur un solide S, par exemple, est le torseur à structure suivant :
g .dm( P ) P ∈Σ {F ( pes → S )}= → AP ∧ g .dm( P ) A P ∈Σ
∫
r
∫
3
r
CENTRE D’INERTIE, OPERATEUR D’INERTIE
Avant de pouvoir employer le principe de conservation de la masse pour simplifier les expressions du torseur cinétique, du torseur dynamique et de l’énergie cinétique, les notions de centre d’inertie et d’opérateur d’inertie vont être rappelées.
3.1
Centre d’inertie G Pour un système matériel Σ de masse m(Σ), en mouvement par rapport à un repère R(O,x,y,z), la position du centre d’inertie G est définie de la façon suivante : → m(Σ ). OG =
∫
→ OP .dm( P )
P ∈Σ
3.2
◊
Opérateur d’inertie J A(Σ) Définition
Pour un système matériel Σ de masse m(Σ), en mouvement par rapport à un repère R(O,x,y,z), l’opérateur d’inertie du système Σ au point A est l’opérateur linéaire JA(Σ) qui a tout vecteur u fait correspondre le vecteur suivant :
∀ u → J A (Σ ).u = r
r
∫
→ r → AP ∧ u ∧ AP .dm( P )
P ∈Σ
On peut également le noter de la façon suivante :
∀ u → J A (Σ ).u = − r
r
∫
P ∈Σ
◊
→ → r AP ∧ AP ∧ u .dm( P ) = −
2 → AP ∧ .dm( P ).u P ∈Σ
∫
r
Relation entre J A( Σ ) et J G( Σ ) ou théorème de Huyghens généralisé .
L’opérateur d’inertie du système Σ en un point A quelconque, JA(Σ) est tel que :
51
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE LA MASSE ET CINETIQUE
∀ u → J A (Σ).u = r
r
→ r → AP ∧ u ∧ AP .dm( P ) =
∫
P ∈Σ
→ → → AG + GP ∧ u ∧ AP .dm( P ) P ∈Σ
∫
r
Soit : J A (Σ ).u
r
∫ ∈Σ
=
→ r → AG ∧ u ∧ AP .dm( P ) +
P
→ r J A (Σ ).u = AG ∧ u ∧ r
∫ ∈Σ
→ r → GP ∧ u ∧ AP .dm( P )
P
→ → → → → AO + OP .dm( P ) + GP ∧ u ∧ AG + GP .dm( P ) P ∈Σ P ∈Σ
∫
∫
r
Alors, G étant le centre d’inertie de Σ : J A (Σ ).u
r
→ → → → → → = m(Σ) AG ∧ u ∧ AG + GP .dm( P ) ∧ u ∧ AG + GP ∧ u ∧ GP .dm( P ) P ∈Σ P ∈Σ
∫
r
∫
r
r
Dont on déduit le théorème de Huyghens généralisé : → r → r r J A (Σ ).u = I G (Σ ).u + m(Σ). AG ∧ u ∧ AG
◊
Expression dans la base (O,x,y,z) →
Si l’on note : EP = x. x + y. y + z . z r
r
r
Alors l’expression de l’opérateur linéaire d’inertie JE(Σ) au point E du système matériel Σ, dans la base (O,x,y,z) est la suivante :
( y 2 + z 2 )dm − ( xy )dm − ( zx )dm Σ Σ Σ 2 2 ( z + x )dm − ( yz )dm J E (Σ ) = − ( xy )dm Σ Σ Σ 2 2 − ( zx )dm − ( yz )dm ( x + y )dm Σ Σ Σ ( x, y, z )
∫
∫
∫ ∫
∫ ∫
∫
∫
∫
r
r r
On note alors :
A − F − E − D J E (Σ ) = − F B − E − D C ( x , y, z ) r
r r
Où A,B et C sont les moments d’inertie respectivement par rapport aux axes (E,x), (E,y), (E,z) et où D, E et F sont appelée les produits d’inertie. -
Base principale d’inertie : Il existe toujours une base orthonormée dans laquelle l’opérateur d’inertie JE(Σ) est diagonal.
-
Influence des symétries : •
A − F 0 Symétrie par rapport au plan z : P(x,y,z)=P’(x,y,-z) : J E (Σ ) = − F B 0 0 0 C ( x , y, z ) r
52
r r
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE LA MASSE ET CINETIQUE
•
A 0 0 Symétrie de révolution autour de l’axe (E,z) : J E (Σ ) = 0 A 0 0 0 C (−,−, z ) r
3.3
Exemples
◊
Opérateur d’inertie en O d’un disque D de rayon R, de centre O et de masse m.
◊
Opérateur d’inertie en O d’un cylindre C de rayon R, de hauteur h et de centre O.
53
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE LA MASSE ET CINETIQUE
◊
Opérateur d’inertie en O d’un cône de révolution C de rayon R et de hauteur h.
◊
Opérateur d’inertie en O d’une sphère creuse S de centre O, de rayon R et de masse m.
1 0 0 2 I O (S ) = mR 2 0 1 0 3 0 0 1 ( x , y, z ) r
◊
r r
Opérateur d’inertie en O d’une sphère pleine S de centre O, de rayon R et de masse m.
1 0 0 2 I O (S ) = mR 2 0 1 0 5 0 0 1 ( x , y, z ) r
4
4.1
r r
CONSEQUENCES DU PRINCIPE DE CONSERVATION DE LA MASSE
Torseur cinétique Le torseur cinétique, du système matériel Σ dans son mouvement par rapport au repère R est, en un point A quelconque, le torseur à structure suivant :
V ( P R ).dm( P ) P ∈Σ {C (Σ R )}= → σ A = AP ∧ V ( P R ).dm( P ) P ∈Σ A
∫
∫
r
◊
r
r
Expression de la résultante cinétique.
O étant l’origine du repère R(O,x,y,z), le centre d’inertie du système S est défini de la façon suivante : → m(Σ ). OG =
∫
→ OP .dm( P )
P ∈Σ
alors
d dt
→ m(Σ ). OG
=
d dt
∫
→ OP .dm( P )
P ∈Σ
54
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE LA MASSE ET CINETIQUE
D’après le principe de conservation de la masse :
→ d → OP .dm( P ) = OP .dm( P ) dt dt R P ∈Σ R P ∈Σ d
∫
∫
Donc la résultante cinétique s’exprime de la façon suivante :
∫
r
r
V ( P R ).dm( P ) = m.V (G R)
P ∈Σ
◊
Expression du moment cinétique.
De même, on cherche l’expression du moment cinétique en A attaché à S : σ A (Σ / R ) = r
→ r AP ∧ V ( P ∈ Σ R ).dm( P )
∫∈Σ
P r
r
→
r
Or : V ( P ∈ Σ R ) = V ( A ∈ Σ R ) + Ω(Σ / R) ∧ AP Donc : σ A (Σ / R ) = r
∫
→ r AP ∧ V ( A ∈ Σ R ).dm( P ) +
P ∈Σ
Donc :
∫
→ r → AP ∧ Ω(Σ / R ) ∧ AP .dm( P )
P ∈Σ
→ r r r σ A (Σ / R ) = m(Σ ). AG ∧ V ( A ∈ Σ R ) + J A (Σ).Ω(Σ / R) r
-
Remarque 1 : si A est fixe dans R alors σ A (Σ / R ) = J A (Σ).Ω(Σ / R )
-
Remarque 2 : si A est le centre d’inertie G alors σ G (Σ / R ) = J G (Σ).Ω(Σ / R )
r
r
r
4.2
Torseur dynamique Le torseur dynamique système matériel Σ dans son mouvement par rapport au repère R s’écrit en un point A quelconque :
Γ( P R ).dm( P ) P ∈Σ {D (Σ R )}= → δ A = AP ∧ Γ( P R ).dm( P ) P ∈Σ A
∫
∫
r
◊
r
r
Expression de la résultante dynamique.
La résultante cinétique est donnée par la relation suivante :
∫
r
P ∈Σ
Si l’on dérive de part et d’autre par rapport au temps :
V ( P R ).dm( P ) = m(Σ) d [V (G R)] R dt dt P ∈Σ R d
∫
r
r
En appliquant maintenant le principe de conservation de la masse :
55
r
V ( P R ).dm( P ) = m(Σ).V (G R)
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE LA MASSE ET CINETIQUE
V ( P R ).dm( P ) = d [V ( P R )] R .dm( P ) = Γ( P R).dm( P ) dt dt P ∈Σ P ∈Σ R P ∈Σ d
∫
∫
r
∫
r
r
On en déduit alors l’expression de la résultante dynamique :
∫
r
r
Γ( P R ).dm( P ) = m(Σ).Γ(G R)
P ∈Σ
◊
Relation entre les moments cinétique et dynamique.
Le moment cinétique en A s’écrit de la façon suivante : σ A (Σ / R ) = r
→ r AP ∧ V ( P R ).dm( P )
∫∈Σ
P
Si l’on dérive de part et d’autre par rapport au temps :
→ AP ∧ V ( P R ).dm( P ) [σ A (Σ / R )] R = dt dt P ∈Σ R d
d
r
∫
r
En appliquant maintenant le principe de conservation de la masse : d dt
[σ A (Σ / R )] R = r
∫ ∈Σ
d → r AP ∧ V ( P R ) dt
∫ ∈Σ
d → AP dt
P
.dm( P )
R
Et en développant : d dt
[σ A (Σ / R )] R = r
P
r
∧ V ( P R ).dm( P ) +
R
→ d r AP ∧ V ( P R ) R .dm( P ) dt P ∈Σ
[
∫
]
Soit : d dt
[σ A (Σ / R )] R = r
∫ (V ( P R) − V ( A R))∧ V ( P R).dm( P ) + δ A (Σ / R) r
r
r
r
P ∈Σ
On en déduit la relation entre le moment cinétique et le moment dynamique : r
δ A (Σ / R ) =
d dt
r
r
[σ A (Σ / R)] R + m(Σ)V ( A R) ∧ V (G R) r
r
Remarque : si A est le centre d’inertie G alorsσ G (Σ / R ) = J G (Σ).Ω(Σ / R) r
-
r
Donc δ G (Σ / R ) =
4.3
d
[ J G (Σ).Ω(Σ / R)] R dt r
Énergie cinétique L’énergie cinétique du système matériel Σ en mouvement par rapport à un repère R s’écrit : T (Σ R ) =
1 2
2 ∫ [V ( P R)] dm( P ) r
P ∈Σ
Si l’on connaît les éléments de réduction du torseur cinématique en un point Q de Σ, on peut écrire :
56
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE LA MASSE ET CINETIQUE
→ V Q R R QP ( ) ( ) 2 T (Σ R ) = + Ω Σ / ∧ P ∈Σ
∫
2
r
dm( P )
On obtient alors en développant l’expression de l’énergie cinétique en un point Q deΣ : → 2 2 T (Σ R ) = m(Σ )[V (Q R )] + Ω(Σ / R ) ⋅ J G (Σ )Ω(Σ / R ) + 2 m(Σ)V (Q R) ⋅ Ω(Σ / R) ∧ QG r
r
Si maintenant le point Q choisi est le centre d’inertie G : 2
2 T (Σ R ) = m(Σ )[V (G R )] + Ω(Σ / R ) ⋅ J G (Σ )Ω(Σ / R ) r
4.4
Autres cas. Exemple : Action de la pesanteur Les simplifications appliquées ci-dessus aux quantités de mouvements, d’accélération ou à l’énergie cinétique s’appliquent également dans les autres cas, comme par exemple dans le cas du torseur de l’action de la pesanteur :
g .dm( P ) m. g P ∈Σ = → A {F ( pes → S )}= AP ∧ g .dm( P ) A AG ∧ m. g A P ∈Σ
∫
∫
r
r
r
r
57
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE LA MASSE ET CINETIQUE
5
◊
A RETENIR
Conséquence du principe de conservation de la masse
Pour un système matériel à masse conservative, définie à l’aide d’une densité volumique ρ, et pour toute fonction régulière h : H (t ) =
∫
h( P ). ρ ( P ).dv( P ) ⇒
Dt
◊
d
H (t ) = dt
∫
d h( P )
Dt
dt
.ρ ( P ).dv( P )
Torseur cinétique
V ( P R ).dm( P ) P ∈Σ {C (Σ R )}= → σ A (Σ / R ) = AP ∧ V ( P R ).dm( P ) P ∈Σ A
∫
r
∫
r
r
∫∈Σ
r
r
V ( P R ).dm( P ) = m.V (G R )
→
P
→ r r r σ A (Σ / R ) = m(Σ). AG ∧ V ( A ∈ Σ R ) + J A (Σ).Ω(Σ / R)
→
r
σ G (Σ / R ) = J G (Σ ).Ω(Σ / R ) r
→
◊
Torseur dynamique
Γ( P R ).dm( P ) P ∈Σ {D (Σ R )}= → δ A (Σ / R ) = AP ∧ Γ( P R ).dm( P ) P ∈Σ A
∫
r
∫
r
r
∫∈Σ
r
r
Γ( P R ).dm( P ) = m(Σ).Γ(G R)
→
P r
→
δ A (Σ / R ) =
→
δ G (Σ / R ) =
r
◊
d dt
r
r
[σ A (Σ / R )] R + m(Σ )V ( A R ) ∧ V (G R) r
d
[ J G (Σ).Ω(Σ / R)] R dt r
Energie cinétique T (Σ R ) =
1 2
2 [ V ( P R )] dm( P ) ∫ r
P ∈Σ
[
r
→
]2 + Ω(Σ / R )⋅ J G (Σ)Ω(Σ / R) + 2 m(Σ)V (Q R) ⋅ Ω(Σ / R) ∧ QG r
Q ∈ Σ → 2 T (Σ R ) = m(Σ ) V (Q R )
58
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE LA MASSE ET CINETIQUE
2
2 T (Σ R ) = m(Σ )[V (G R )] + Ω(Σ / R ) ⋅ J G (Σ )Ω(Σ / R ) r
59
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE L’ENERGIE
1
Introduction ______________________________________________________ 61
2
Energétique ______________________________________________________ 61 2.1
Torseur des actions mécaniques extérieures à un solide ___________________ 61
2.2
Puissance__________________________________________________________ 61
◊ ◊
Puissance associée à des actions extérieures. ___________________________________ 61 Puissance associée à des actions réciproques. ___________________________________ 62
2.3
Travail____________________________________________________________ 63
2.4
Energie Potentielle __________________________________________________ 63
◊ ◊ ◊ ◊ 2.5
Energie potentielle associée à des efforts extérieurs ______________________________ Energie potentielle associée à des actions mutuelles____________________________ __ Exemple : énergie potentielle associée aux inter-efforts gravitationnels. ______________ Quelques actions mutuelles avec énergie potentielle associée ______________________
63 64 64 65
Energie cinétique ___________________________________________________ 66
3
Conservation de l’énergie : Théorème de l’énergie cinétique _______________ 66 ◊ Théorème de l’énergie pour un solide S _______________________________________ 66 ◊ Théorème de l’énergie cinétique pour deux solides S1 et S2 _______________________ 66 Théorème de l’énergie cinétique pour un système Σ de n solides ____________________ 67 ◊ Intégrale première de l’énergie cinétique : Système conservatif_____________________ 67 ◊
4
A retenir _________________________________________________________ 68
60
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE L’ENERGIE
1
INTRODUCTION
Dans cette partie nous allons rappeler rapidement les principales grandeurs employées en énergétique (puissance, travail, énergie potentielle et énergie cinétique) et les exprimer en employant le formalisme défini pour la mécanique des milieux indéformables. Dans un second temps nous écrirons le principe de conservation de l’énergie, en employant le formalisme associé à la cinématique des systèmes de solides indéformables.
2
2.1
ENERGETIQUE
Torseur des actions mécaniques extérieures à un solide Les efforts extérieurs à un solide, sont représentés par des champs vectoriels définis par rapport à une mesure, volumique, surfacique ou linéique. On peut donc leur associer des torseurs. Supposons un solide S. Les actions mécaniques extérieures à S sur S sont représentées par un champ vectoriel Ψ.dP. Où dP est une mesure volumique, surfacique ou linéique. S’il s’agit d’actions à distance (attraction gravitationnelle, par exemple), l’intégrale est une intégrale volumique. S’il s’agit d’actions de contact, l’intégrale est alors une intégrale de surface. Alors le torseur des actions de S sur S s’écrit au point A :
ΨS → S .dP S F {S → S }= AP ∧ ΨS → S .dP A S
∫
∫
2.2
◊
Puissance
Puissance associée à des actions extérieures. Se étant un système matériel appartenant à Si . Les actions mécaniques de Se sur Si sont représentées par un champ vectoriel Ψ.dP. Où dP est une mesure volumique, surfacique ou linéique. La puissance développée par les actions mécaniques de Se sur Si, pour la loi k (contact, attraction gravitationnelle …), dans le mouvement de Si par rapport au repère R, à la date t, est par définition :
k
P Se → Si / R =
∫
r
ΨSe→Si ⋅ V ( P i / R )dP
Si
Ceci s’écrit aussi :
k
k
P Se → Si / R = F Se → Si ×V ( Ri / R )
61
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE L’ENERGIE
Alors :
k
k
P Se → Si / R1 =
P Se → Si / R 2 =
∫
r
ΨSe→Si ⋅V ( P i / R1)dP
Si
∫
r
ΨSe→ Si ⋅V ( P i / R 2)dP
Si
Par conséquent :
k
k
k
k
∫
P Se → Si / R1 − P Se → Si / R 2 =
ΨSe→ Si ⋅ [V ( P i / R1) − V ( P i / R 2)]dP r
r
Si
Soit :
k
∫
P Se → Si / R1 − P Se → Si / R 2 =
ΨSe→ Si ⋅ [V ( P i ∈ R1 / R 2)]dP r
Si
Soit encore :
k
k
P Se → Si / R1 − P Se → Si / R 2 = F Se → Si × V ( R1 / R 2)
Σ est un système de solides Si, et Se étant un système matériel appartenant à Σ . La puissance développée par les actions mécaniques de Se sur Σ, pour la loi k (contact, attraction gravitationnelle …), dans le mouvement de Σ par rapport au repère Galiléen R, à la date t, est par définition :
k
P Se → Σ / R =
n
∑∫Ψ
⋅ V ( P j / Rg )dP r
Se→ Sj
j =1 Si
Ce qui s’écrit aussi :
k n k P Se → Σ / R = F Se → Sj ×V ( Rj / R ) j =1
∑
◊
Puissance associée à des actions réciproques. Soit deux systèmes disjoints S1 et S2 en mouvement par rapport à R. La puissance des actions réciproques entre S1 et S2, pour la lo i d’interaction k, s’écrit :
k
k
k
P S 1 ↔ S 2 / R = P S 1 → S 2 / R + P S 2 → S 1 / R
Par conséquent :
k
k
k
k
k
k
P S 1 ↔ S 2 / R1 = P S 1 → S 2 / R1 + P S 2 → S 1 / R1
De même :
k
k
P S 1 ↔ S 2 / R 2 = P S 1 → S 2 / R 2 + P S 2 → S 1 / R 2
Donc :
k
k
P S 1 ↔ S 2 / R1 − P S 1 ↔ S 2 / R 2 = F S 1 → S 2 ×V ( R1 / R 2) + F S 2 → S 1 ×V ( R1 / R 2)
62
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE L’ENERGIE
k
k
Or : F S 1 → S 2 = −F S 2 → S 1
k
k
Donc P S 1 ↔ S 2 / R1 − P S 1 ↔ S 2 / R 2 = 0
La puissance des actions réciproques entre deux solides pou toute loi d’interaction k est
k
indépendante du référentiel du mouvement choisi. On la note donc : P S 1 ↔ S 2
2.3
Travail
Se étant un système matériel appartenant à Si . Le travail entre les dates t1 et t2 des actions mécaniques de Se sur Si, pour la loi k (contact, attraction gravitationnelle …), dans le mouvement de Si par rapport au repère R, est par définition : t =t 2
t 2
k k W Se → Si / R = P Se → Si / R dt t 1 t =t 1
2.4
∫
Energie Potentielle On suppose deux solides S1 et S2 en mouvement par rapport à un référentiel R.
◊
Energie potentielle associée à des efforts extérieurs La puissance développée par les efforts exercés par S2 sur S1, en accord avec la loi d’interaction k, s’écrit :
k
P S 2 → S 1 / R =
∫
r
r
ΨS k 2→S 1 ( P 2, P 1).V ( P 1 / R )dP 1.dP 2
S 1×S 2
On procède à une double intégration. On intègre l’effet de chacun des points P2 de S2 sur un même point P1 de S1. Puis on intègre sur l’ensemble des points P1 de S1. Alors, si quel que soit le mouvement de S1 par rapport à R on peut écrire :
k
P S 2 → S 1 / R = −
d dt
∫ Θ( P 1, P 2).dP 1.dP 2 S 1×S 2
On peut associer une énergie potentielle aux actions mécaniques de S2 sur S1 dans son mouvement par rapport à R, avec :
k
∫ Θ( P 1, P 2).dP 1.dP 2
E p S 2 → S 1 / R =
S 1×S 2
Et dans ce cas :
k
P S 2 → S 1 / R = −
d
k
E p S 2 → S 1 / R dt
63
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE L’ENERGIE
◊
Energie potentielle associée à des actions mutuelles De même, si quels que soient les mouvements de S1 et S2 par rapport à R quelconque, la puissance des actions mutuelles entre S1 et S2 peut s’écrire explicitement :
k
P S 2 ↔ S 1 = −
d dt
∫ Θ( P 1, P 2).dP 1.dP 2 S 1×S 2
On pourra alors associer une énergie potentielle aux actions mutuelles entre S1 et S2 :
k
∫ Θ( P 1, P 2).dP 1.dP 2
E p S 2 ↔ S 1 =
S 1×S 2
Et comme précédemment :
k
P S 2 ↔ S 1 = −
d
k
E p S 2 ↔ S 1 dt
Dans tout les cas l’énergie potentielle n’est définie qu’à une constante additive près qu’on omet d’écrire en général.
◊
Exemple : énergie potentielle associée aux inter-efforts gravitationnels. L’attraction gravitationnelle d’un point P1 du solide S1 sur un point P2 du solide S2 en mouvement par rapport à un référentiel Rg, s’écrit comme suit : r
g Ψ P 1→ P 2
P 1 P 2
( P 1 , P 2 ) = − G o .dm( P 1 )dm( P 2 )
P 1 P 2
3
Où G o est la constante de gravitation universelle.
Alors l’action du solide S1 sur un point P2 du solide S2 s’écrit donc :
∫
r
( P ) = − G o . dm( P 1 )dm( P 2 ) Ψ g S 1→ P 2 2
P 1 P 2 P 1 P 2
S 1
3
Et la puissance développée par l’attraction gravitationnelle de S1 sur S2, où S2 est en mouvement par rapport au repère Rg, s’écrit : r g g P S 1 → S 2 / Rg = Ψ S 1→ P 2 ( P 2 ) ⋅ V ( P 2 / Rg ) S 2 r
∫
r
Par ailleurs, la puissance développée par l’attraction gravitationnelle de S2 sur S1, où S1 est en mouvement par rapport au repère Rg, s’écrit : r g g P S 2 → S 1 / Rg = ΨS 2→ P 1 ( P 1 ) ⋅ V ( P 1 / Rg ) S 2 r
∫
r
64
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE L’ENERGIE
Par conséquent la puissance développée par les actions mutuelles entre S1 et S2 est :
g g g P S 2 ↔ S 1 = P S 1 → S 2 / Rg + P S 2 → S 1 / Rg r
r
r
En développant : r r g P 1 P 2 P S 2 ↔ S 1 = − G o . ( V ( P 2 / Rg ) − V ( P 1 / Rg ))dm ( P 1 )dm( P 2 ) 3 S 1×S 2 P P r
∫
1 2
g P 1 P 2 P S 2 ↔ S 1 = G o . 3 P S 1× S 2 P 1 2 r
∫
d ⋅ P 1 P 2 dm( P 1 )dm( P 2 ) dt Rg
Soit finalement :
g d 1 P S 2 ↔ S 1 = G o . dt P 1 P 2 S 1× S 2 r
∫
dm( P 1 )dm( P 2 ) Rg
Et comme il s’agit de systèmes à masse conservative :
g d dm( P 1 )dm( P 2 ) P S 2 ↔ S 1 = G o . P dt S 1× S 2 1 P 2 r
∫
Donc il existe une énergie potentielle associée aux actions mutuelles de l’attraction gravitationnelle entre deux solides S1 et S2 qui s’écrit :
g dm( P 1 )dm( P 2) E p S 2 ↔ S 1 = − G o . P 1 P 2 S 1× S 2 r
∫
◊
Quelques actions mutuelles avec énergie potentielle associée
-
g dm( P 1 )dm( P 2) Attraction gravitationnelle entre deux solides S1 et S2 : E p S 2 ↔ S 1 = − G o . P 1 P 2 S 1× S 2
-
Attraction gravitationnelle à la surface de la terre entre la terre et un solide S : approximation !
r
∫
g g r E p Terre ↔ S = E p Pes → S / Terre = −m(S ). g ⋅ OG r
r
O étant le point de référence à la surface de la terre, g le vecteur accélération de la pesanteur en O et G le centre d’inertie du solide S.
-
r
Ressort de traction-compression : E p S 2 ↔ S 1 =
k 2
(l − l o )2
Où r est un ressort de traction compression, de raideur k, et de longueur libre lo, et où l paramètre la distance entre les points d’attache du ressort sur les deux solides S1 et S2.
-
r
Ressort de torsion : E p S 2 ↔ S 1 =
C 2
(θ − θ o )2
Où r est un ressort de torsion, de raideur C, et de rotation libre θo, et où θ paramètre la torsion du ressort.
65
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE L’ENERGIE
2.5
Energie cinétique L’énergie cinétique du système matériel Σ de centre d’inertie G dans son mouvement par rapport à un repère R est le scalaire suivant : T (Σ R ) =
1 2
2 ∫ [V ( P R)] dm( P ) r
P ∈Σ
Ceci s’écrit aussi : T (Σ R ) =
3
◊
1 2
[
]2 + 12 Ω(Σ / R ) ⋅ J G (Σ)Ω(Σ / R)
r
m(Σ ) V (G R )
CONSERVATION DE L’ENERGIE : THEOREME DE L’ENERGIE CINETIQUE
Théorème de l’énergie pour un solide S Le principe de conservation de l’énergie conduit à écrire, pour un solide S, que la puissance développée par les actions mécaniques extérieures au solide S, est égale à la variation d’énergie cinétique du solide soit :
(
)
P S → S / Rg =
◊
dT (S / Rg ) dt
Théorème de l’énergie cinétique pour deux solides S1 et S2 Si l’on considère maintenant un système Σ constitué de deux solides S1 et S2, on peut appliquer le principe de conservation de l’énergie à chacun des deux solides :
(
)
P S 1 → S 1 / Rg =
Donc :
Soit :
dT (Σ / Rg ) dt
dT (Σ / Rg ) dt
dT (Σ / Rg ) dt
dT (S 1 / Rg ) dt
=
dT (S 1 / Rg ) dt
(
+
)
P S 2 → S 2 / Rg =
et
dT (S 2 / Rg ) dt
dT (S 2 / Rg ) dt
= P (S 1 → S 1 / Rg ) + P (S 2 → S 2 / Rg )
= F (S 1 → S 1)×V (S 1 / Rg ) +F (S 2 → S 2)×V (S 2 / Rg )
[
]
[
Donc :
]
= F (Σ → S 1)+F (S 2 → S 1) ×V (S 1 / Rg ) + F (Σ → S 2)+F (S 1 → S 2) ×V (S 2 / Rg )
Comme F (S 2 → S 1) = −F (S 1 → S 2) On en déduit :
Ou encore
dT (Σ / Rg ) dt
=
∑F (Σ → Si)×V (Si / Rg ) +F (S 1 → S 2)× [V (S 2 / Rg ) − V (S 1 / Rg )] i =1,2
dT (Σ / Rg ) dt
= P (Σ → Σ / Rg )+ P (S 1 ↔ S 2)
La puissance développée par les efforts extérieurs et intérieurs au système Σ est égale à la variation de l’énergie cinétique du système. Ceci se généralise naturellement à un système de n solides.
66
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE L’ENERGIE
◊
Théorème de l’énergie cinétique pour un système Σ de n solides La généralisation à un système Σ de n solides donne : dT (Σ / Rg ) dt
n
=
∑ F (Σ → S )×V (S i
n
i
/ Rg ) +
∑∑ F (S
i =1
j
→ S i )×V (S i / S j )
i =1 j >i
Le théorème de l’énergie cinétique pour un système Σ de n solides s’écrit donc : dT (Σ / Rg ) dt
◊
=
(
)
P Σ → Σ / Rg
(
+
P S i ↔ S j
14 4 244 3
Puissance des efforts extérieurs
)
14243
Puissance des efforts intérieurs
Intégrale première de l’énergie cinétique : Système conservatif Si les puissances des efforts intérieurs et extérieurs sont nulles ou bien se calculent par dérivation d’une énergie potentielle, dont on notera la somme E p(Σ/Rg), alors le théorème de l’énergie cinétique s’écrit : dT (Σ / Rg ) dt
=−
dE p (Σ / Rg ) dt
Et par conséquent, il existe une intégrale première du mouvement, appelée intégrale première de l’énergie cinétique : T (Σ / Rg ) + E p (Σ / Rg ) = Constante
Cette intégrale première traduit directement la conservation de l’énergie du système. Il y a transformation de l’énergie potentielle en énergie cinétique et réciproquement. Le système est alors dit « conservatif ».
67
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables CONSERVATION DE L’ENERGIE
4
A RETENIR
•
Torseur des actions mécanique de S sur S, au point A :
ΨS → S .dP F {S → S }= S AP ∧ ΨS → S .dP A S
∫
∫
•
Energie potentielle pour quelques lois d’interaction :
•
g dm( P 1 )dm( P 2) Attraction gravitationnelle entre deux solides S1 et S2 : E p S 2 ↔ S 1 = − G o . P 1 P 2 S 1× S 2 r
•
∫
Attraction gravitationnelle à la surface de la terre, entre la terre et un solide S en mouvement par rapport à la terre : g g r E p Terre ↔ S = E p Pes → S / Terre = −m(S ). g ⋅ OG r
•
r
r 1 1 2 2 Ressort de traction-torsion : E p S 1 ↔ S 2 = K ( L − Lo ) + C (θ − θ o ) 2 2
Où r est un ressort de traction torsion, de raideur k en traction-compression, de raideur C en torsion, et de longueur libre lo, et où l paramètre la distance entre les points d’attache du ressort sur les deux solides S1 et S2 et de torsion libre θo où θ paramètre la torsion du ressort. •
Théorème de l’énergie cinétique pour un système Σ de n solides dT (Σ / Rg ) dt
•
(
=
)
P Σ → Σ / Rg
(
+
Puissance des efforts extérieurs
)
P S i ↔ S j 14243
14 4 244 3
Puissance des efforts intérieurs
Intégrale première de l’énergie cinétique : Conservation de l’énergie
(
)
(
)
Si P Σ → Σ / Rg = 0 et P S i ↔ S j = 0 Si P Σ → Σ / Rg = −
d
(
)
E p Σ → Σ / Rg dt
Alors T (Σ / Rg ) = Constante et-
Alors T (Σ / Rg ) + E p (Σ / Rg ) = Constante
68
(
)
P S i ↔ S j = −
d
(
E p S i ↔ S j dt
)
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
1
Introduction ______________________________________________________ 70
2
Actions mécaniques ou efforts _______________________________________ 70 2.1
Torseur des actions mécaniques extérieures à un système de solides
2.2
Exemple d’action à distance : la pesanteur ______________________________ 71
2.3
Actions de contact : Loi de Coulomb ___________________________________ 71
◊ ◊ ◊ ◊ ◊ 2.4
72 73 73 74 75
Exemple : principe d’un étouffeur de vibrations_________________________________ 76
Principe fondamental de la dynamique ________________________________ 77 3.1
Introduction : un peu d’histoire _______________________________________ 77
3.2
Enoncé du principe fondamental de la dynamique________________________ 78
3.3
Conséquences ______________________________________________________ 78
◊ ◊ ◊ ◊ ◊ 3.4
Théorème de la résultante dynamique_________________________________________ Théorème du moment dynamique____________________________________________ Théorème des actions mutuelles _____________________________________________ Cas particulier de la statique ________________________________________________ Equations du mouvement __________________________________________________
78 79 79 79 79
Référentiels Galiléens/non Galiléens ___________________________________ 79
◊ ◊ ◊ 4
Lois de Coulomb relatives à la résultante.________________________________ ______ Lois de Coulomb relatives au moment en I. ____________________________________ Contact non ponctuel, exemple d’une liaison pivot sans frottement __________________ Liaisons parfaites : Actions de contact sans frottement. ___________________________ Lois de comportement de liaisons. ___________________________________________
Graphe d’analyse ___________________________________________________ 76
◊ 3
_______ 70
Référentiels Galiléens _____________________________________________________ 80 Référentiel non Galiléen ___________________________________________________ 80 Exemple : Accélération de la pesanteur. _______________________________________ 81
Principe des puissances virtuelles _____________________________________ 83 4.1
Introduction : un peu d’histoire. ______________________________________ 83
4.2
Enoncé du principe des puissances virtuelles ou PPV. _____________________ 83
4.3
Choix de torseurs virtuels particuliers et théorèmes de la dynamique ________ 84
4.3.1 Torseur global quelconque : Equivalence du principe des puissances virtuelles et du principe fondamental de la dynamique ___________________________________________________________ 84 4.3.2 Torseur des vitesses galiléennes : Théorème de l’énergie cinétique __________________ 84 ◊ Théorème de l’énergie cinétique pour un solide unique.___________________________ 84 ◊ Théorème de l’énergie cinétique pour un système de solides. ______________________ 84 4.3.3 Torseurs de Lagrange : équations de Lagrange__________________________________ 85 ◊ Définition des torseurs de Lagrange : _________________________________________ 85 ◊ Conséquences ___________________________________________________________ 85 ◊ Application du PPV à un unique solide S.______________________________________ 86 ◊ Application du PPV à un système de solides Σ : _________________________________ 86 ◊ Fonction de force. ________________________________________________________ 87 4.3.4 Combinaison des torseurs de Lagrange : Equation de Painlevé _____________________ 89 ◊ Définition du torseur de Painlevé:____________________________________________ 89 ◊ Conséquences ___________________________________________________________ 90 ◊ Intégrale première de Painlevé ______________________________________________ 90
5
A retenir _________________________________________________________ 92
69
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
1
INTRODUCTION
Le principe fondamental de la dynamique a été énoncé pour la première fois vers la fin du XVII siècle par Newton. Il permet de prévoir avec une très bonne précision les phénomènes mécaniques classiques. Cette présentation traditionnelle est actuellement remplacée par une présentation équivalente, le principe des puissances virtuelles qui permet de faire le lien entre les différentes branches de la mécanique, mécanique des structures, des plaques et coques et mécanique des milieux continus ... Dans ce chapitre, nous allons d’abord présenter le principe fondamental de la dynamique dans sa forme classique puis le principe des puissances virtuelles et les différents théorèmes qui en découlent. Le choix de l’une ou l’autre des expressions de ce principe sera généralement guidé par la simplicité de mise en œuvre pour un problème précis. Le principe des puissances virtuelles et les équations de Lagrange seront généralement employés pour déterminer les mouvements de systèmes de solides.
2
2.1
ACTIONS MECANIQUES OU EFFORTS
Torseur des actions mécaniques extérieures à un système de solides Σ Les efforts extérieurs à un solide ou à un système de solides, sont représentés par des champs vectoriels Ψ.dP définis par rapport à une mesure dP, volumique, surfacique ou linéique. S’il s’agit d’action à distance (attraction gravitationnelle, par exemple), l’intégrale est une intégrale volumique. S’il s’agit d’actions de contact, l’intégrale est alors une intégrale de surface. On peut leur associer des torseurs. Supposons un solide S et un solide Σe. Le torseur des actions mécanique de point A :
Σe sur S s’écrit au
ΨΣe→S .dP Si F {Σ e → S }= AP ∧ ΨΣe→S .dP A Σ
∫
∫
Lorsque plusieurs solides interviennent, les contributions de chacun des solides sont distinguées. Dans le cas particulier ci-dessous, si l’on note Σ l’extérieur de Σ, Σ est ici constitué des deux solides Sa et Sb. Alors :
F {Σ → Σ} = F {Sa → Σ} + F {Sb → Σ}
Par ailleurs, sur chacun des solides, le contact se produit en deux zones, Ca1 et Ca2 pour Sa et Cb1 et Cb2 pour Sb :
70
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
ΨSa→Σ .dP + ΨSa→Σ .dP Ca1 Ca 2 Ca 1 Ca 2 F {Sa → Σ}= = F Sa → Σ + F Sa → Σ AP ∧ ΨSa→Σ .dP + AP ∧ ΨSa→Σ .dP Ca1 Ca 2 A
∫
∫ ∫
∫
On écrira ainsi pour chacune des liaisons entre solides, et pour chacune des actions à distances, un torseur des actions mécaniques associé.
2.2
Exemple d’action à distance : la pesanteur Le torseur, en un point A, des actions de la pesanteur sur un solide S de masse m et de densité (homogène) ρ d’écrira comme suit : r ρ . g .dv g S F pes → S = r AP ∧ ρ . g .dv A S
∫
r
r
r
g m. g Ou encore : F pes → S = r A AG ∧ m g
∫
n.b. On reviendra sur ce cas particulier
2.3
Actions de contact : Loi de Coulomb
Supposons deux solides S1 et S2 en contact en un point I. Il existe un unique plan tangent Π entre les deux solides, défini par la normale n21, à S2 ou à S1, au point de contact I, dirigée de S2 vers S1. Les éléments de réduction du torseur cinématique en I du mouvement de S1 par rapport à S2 sont : r
Ω( S 1 / S 2) { V ( S 1/ R 2 )}= r I G ( I ∈ S 1 / S 2) I Par ailleurs le torseur des actions de contact au point I se note :
c Rr [F ( S 2 → S 1)] F ( S 2 → S 1)= r I M [I,F ( S 2 → S 1)] I
71
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
On pose alors : •
r
r
= n12 ⋅ R[F ( S 2 → S 1)] r r r r T 21 = n12 ∧ ( R[F ( S 2 → S 1)] ∧ n12 ) r r M 21 ( I ) = M [I,F ( S 2 → S 1)] N 21
N21 : Composante normale
•
T21 : Composante tangentielle :
•
M21(I) : Moment en I :
•
M21P(I) : Moment de résistance au pivotement en I :
r
(r
)r
M 21 P ( I ) = n12 ⋅ M [ I ,F (S 2 → S 1)] n12 •
M21R (I) : Moment de résistance au roulement en I :
r
M 21 R ( I ) = n12
◊
r
r
∧ ( M [ I ,F (S 2 → S 1)] ∧ n12 )
Lois de Coulomb relatives à la résultante. La nature du contact entre les deux solides est caractérisée par un coefficient de frottement sans dimension f (coefficient de frottement de glissement). Naturellement ces lois n’ont de sens que si le contact est maintenu, soit N21>0. r
On admet dans ce cas
M [I,F ( S 2 → S 1)] = 0
Si f =0 (liaison parfaite)
T 21
r
r
=0
Si f ∫0 r
Condition de non-glissement en I :
r G ( I ∈ S 1 / S 2 )
Condition de glissement en I :
r r r G ( I ∈ S 1 / S 2) ≠ 0 ⇒ T 21
=
0
⇒
r
T 21
≤
f . N 21
= f . N 21
Et dans le cas où on a glissement : r
r
r
T 21 ∧ G ( I ∈ S 1 / S 2) = 0 r
r
T 21 ⋅ G ( I ∈ S 1 / S 2) < 0
Soit :
r
T 21
=
r G ( I ∈ S 1 / S 2) − f . N 21 . r G ( I ∈ S 1 / S 2)
En pratique on en déduit deux équations scalaires dans le plan de contact. Ainsi, la résultante du torseur des actions de contact appartient t’elle à un cône, dit cône de frottement (voir la figure ci-dessous). Le coefficient de frottement f est la tangente de l’ouverture de ce cône ou angle de frottement φ.
72
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
◊
Lois de Coulomb relatives au moment en I. Comme dans le cas des lois de Coulomb relative à la résultante, dans le cas des lois de Coulomb relative au moment en I, on décompose le moment en I du torseur des actions mécaniques de S2 sur S1 en un terme normal au plan de contact, appelé résistance au pivotement et en un terme appartenant au plan de contact, appelé résistance au roulement. Par ailleurs la résultante du torseur cinématique se décompose en un terme de roulement et un terme de pivotement. Si le torseur cinématique en I du mouvement de S1 par rapport à S2 s’écrit : r
Ω(S 1 / S 2) r I {V ( S 1/ R 2 )}= G ( I ∈ S 1 / S 2) I On note alors : r
r
r
r
r
r
r
Ω(S 1 / S 2) = (1 Ω4 ( S 1 / S 2) ⋅ n 21 )n 21 + n 21 ∧ (Ω( S 1 / S 2) ∧ n 21 ) 4 244 3 144 4 4 244 4 4 3 r
P 1 / 2 : Pivotement
R1 / 2 : Roulement
R1 / 2
Condition de roulement en I : Alors
M 21 R ( I ) ⋅ R1 / 2
<0
≠0
et
Condition de non-roulement en I :
R1 / 2
M 21 R ( I )
= h. N 21
M 21 R ( I )
< h. N 21
=0
Alors
Des lois équivalentes sont écrites pour le pivotement : P 1 / 2
Condition de pivotement en I : Alors
r
M 21 P ( I ) ⋅ P 1 / 2 n 21
<0
Condition de non-pivotement en I : Alors
≠0
et P 1 / 2
M 21 P ( I )
= k . N 21
M 21 P ( I )
< k . N 21
=0
Le coefficient de résistance au roulement h et le coefficient de résistance au pivotement k sont homogènes à une longueur. En outre les moments résistants au roulement et au pivotement sont proportionnels à l’effort normal. En effet, les solides sont considérés ici comme indéformables, mais ils sont en réalité toujours déformables. De ce fait, plus l’effort normal est grand, plus l’aire de contact entre les deux solides est importante et plus la conformation des deux surfaces est marquée. Le roulement nécessite de quitter cette configuration adaptée et de conformer de nouvelles surfaces ce qui nécessite d’autant plus d’énergie que la surface de contact est grande, c'est-à-dire que l’effort normal est grand. Naturellement, la dimension de l’aire de contact dépend également de la déformabilité des solides au niveau du contact. Vous pouvez, par exemple, expérimenter facilement à l’aide d’un vélo et d’une côte que le coefficient de résistance au roulement h pour le contact du pneu avec la route diminue avec la pression des pneus.
◊
Contact non ponctuel, exemple d’une liaison pivot sans frottement Dans cet exemple « théorique » la liaison pivot est réalisée à l’aide d’une surface de contact de révolution s, d’axe (O,δ), entre deux solides S1 et S2, qui est non cylindrique.
73
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
Comme la liaison est sans frottement, f =0, alors en tout point P de la surface de contact s , la composante tangentielle des actions de contact de S1 sur S2 est nulle donc :
r
ψ S 1→ S 2 ( P ) =
r
Ψ S 1→ S 2 ( P ).n12 ( P ).d
Alors le torseur en O des actions de contact de S1 sur S2 s’écrit : c r r ( → F R S 1 S 2) = ΨS 1→S 2 ( P ).n12 ( P ).ds( P ) c s ( → ) F S 1 S 2 = r c r M O,F ( S 1 → O S 2) = OP ∧ ΨS 1→S 2 ( P ).n12 ( P ).ds( P ) s
∫
∫
O
On ne peut rien dire de particulier de la résultante, en revanche comme il s’agit d’une surface de révolution le moment dynamique en O est aligné avec l’axe de la liaison : r Γ( P S 2).dm ( P ) P ∈S 1 → r O {D (S 1 S 2)}= r δ O = OP ∧ Γ(S 1 S 2).dm( P ) P ∈S 1 O
∫
∫
Donc on peut écrire pour le moment du torseur des actions de contact : r
r
c
δ O ⋅ M O,F ( S 1 → S 2)
r
r
c
=
r
r ( ∫ 144 4 244 4 3)
δ O ⋅ OP ∧ .n12 ( P ) ΨS 1→ S 2 ( P ).ds( P )
=0
s
Soit : δ O ⋅ M O,F ( S 1 → S 2) = 0
◊
Liaisons parfaites : Actions de contact sans frottement. On procède de même pour chacune des liaisons normalisées par contact sans frottement entre S1 et S2, on pourra noter qu’il existe autant d’équations que de degrés de liberté pour la liaison. Les différentes équations de liaison pour les liaisons parfaites (sans frottement) sont regroupées dans le tableau cidessous. A part la liaison hélicoïdale, les particularités des torseurs d’actions de contact se trouvent facilement en mettant en évidence des propriétés géométriques liant les caractéristiques géométriques des liaisons et l’ensemble des normales aux surfaces de contact. Noms
Caractéristiques
Equations de liaisons r r
r
glissière
x
r
pivot
( A, x )
hélicoïdale
( A, x ) , pas réduit p
r
74
r
R ⋅ x
=0
r
M A ⋅ x
r
=0
r
r
( M A
+ p R) ⋅ x = 0
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
A,
r r
r r
x , y avec x ⋅ y = 0
rotule
r
M A ⋅ x
( A, x )
rotule à doigt
r
r
M A ⋅ x
r
( A, x ) ,
r
r
M A r
r
linéique rectiligne
r r
= 0 , R ⋅ x = 0 r
= 0 , M A ⋅ y = 0
A
appui plan
◊
r
r
pivot glissant
r
=0 r
r
r
=0
r
r
M A ⋅ z = 0 , R ∧ x
r
r r
z avec x ⋅ z = 0
r
M A
r
linéique annulaire
A,
z
ponctuelle
P,
z
r
r
M A
r
r
r
∧ ( x ∧ z ) = 0 , R ∧ x = 0 r
r
r r
r
r
= 0 , R ⋅ z = 0 r
M P = 0 , R = R z
Lois de comportement de liaisons. En dehors des liaisons parfaites, il existe aussi des liaisons élastiques, assimilables à des ressorts de traction compression pour les degrés de liberté en translation et à des ressorts de torsion pour les degrés de liberté en rotation. Il existe également des liaisons dissipatives, telles que des amortisseurs par exemple. Les équations de liaison associées à ces lois de comportement de liaisons sont présentées cidessous.
-
Ressort de traction compression Supposons deux solides S1 et S2 liés par une liaison glissière élastique de direction u. La projection de la résultante des actions mécaniques de S2 sur S1 selon la direction de la liaisonu s’écrit alors, où λ=l-lo est l’allongement du ressort :
r
r
R (F (S 2 → S 1)) ⋅ u
-
= −k λ
Ressort de torsion Supposons deux solides S1 et S2 liés par une liaison pivot élastique de direction u. La projection du moment en un point de l’axe A des actions mécaniques de S2 sur S1 selon la direction de la liaison u s’écrit alors, où δθ=θ-θo est l’angle de torsion du ressort : r
r
M (A ,F (S 2 → S 1)) ⋅ u
-
= −k δθ
Amortisseur de traction-compression Supposons deux solides S1 et S2 liés par une liaison glissière avec amortissement de direction u. La projection de la résultante des actions mécaniques de S2 sur S1 selon la direction de la liaisonu s’écrit, où A est le point d ‘attache de l’amortisseur sur le solide S1 : r
r
R (F (S 2 → S 1)) ⋅ u
-
r
r
= −d (V ( A ∈ S 1 / S 2) ⋅ u )
Amortisseur de torsion De même, supposons deux solides S1 et S2 liés par une liaison pivot avec amortissement de direction u. La projection du moment en A des actions mécaniques de S2 sur S1 selon la direction de la liaison u s’écrit, où A est le point d ‘attache de l’amortisseur sur le solide S1 : r
r
M (A,F (S 2 → S 1)) ⋅ u
r
r
= −d (Ω (S 1 / S 2) ⋅ u )
75
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
2.4
Graphe d’analyse Dès que les actions mécaniques appliquées à un système de solides sont bien identifiées, le graphe de structure peut être complété afin de construire un graphe d’analyse. Dans un graphe de structure, les sommets sont les solides et les arcs les liaisons. A chacun des arcs on pourra donc associer autant d’équations de liaison qu’il y a de degré de liberté dans la liaison. Dans un graphe d’analyse on fait également figurer les actions mécaniques, et les lois de comportement. Ce graphe permet de bien identifier les équations à mettre en place, les axes de projection à choisir et les sous-systèmes de solides à identifier. L’intérêt de ce graphe est illustré ci-dessous à l’aide d’un exemple, qui sera repris par la suite, lorsque les théorèmes fondamentaux auront été tous exposés.
◊
Exemple : principe d’un étouffeur de vibrations Le système étudié est constitué d’un bâti So. Sur ce bâti est fixé un support de machine tournante S1 par une liaison glissière d’axe (O,xo). Un ressort de rappel est monté entre S1 et S0. Par ailleurs, sur le support machine est monté sur une liaison pivot d’axe (G1,y0), un moteur dont la masse est incluse dans la masse M1 du support. La partie tournante du moteur présente un bâlourd S3, de masse M3, dont le centre d’inertie G3 est à la distance R du centre d’inertie G1 du support de la machine tournante. Enfin un pendule simple S2, l’étouffeur, est également attaché par une liaison pivot d’axe (G1,y0) au support S1. La masse de l’étouffeur est notée M2 et la distance G1G2 est notée d .
-
Objectifs : Déterminer les caractéristiques à donner au pendule S2 (masse M2, distance d) pour qu’il jour effectivement un rôle d’étouffeur de vibrations, le reste du système étant supposé donné.
-
Graphe de structure : On a quatre solides en présence, S0, S1, S2 et S3. On construit le graphe de structure en précisant les équations de liaison associées et les paramètres du mouvement. On peut alors compléter le graphe de manière à faire apparaître les actions mécaniques. Ici : l’action de la pesanteur et l’action du ressort monté entre S0 et S1. Dans cet exemple, le degré de mobilité est de trois. Nous avons trois inconnues cinématiques, X le déplacement horizontal de S1/S0, θ la rotation de S2/S1 et φ la rotation de S3/S1. Nous devons donc établir trois équations scalaires.
-
Analyse En utilisant le graphe d’analyse, on identifie les équations nécessaires et les sous-systèmes à isoler. Il nous faut trois équations, puisque le degré de mobilité est de trois. Il suffit de regarder le graphe d’analyse que nous venons de construire pour déterminer de quel type d’équation nous aurons besoin : 1 : une équation portant sur la résultante des actions mécaniques de S0 sur (S1 U S2 U S3). Cette équation sera projetée sur l’axe xo pour lequel n’apparaît pas d’inconnue de liaison.
76
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
2 : une équation portant sur le moment en G1 des actions mécaniques de S3 sur S1. Cette équation sera projetée sur l’axe yo pour lequel n’apparaît pas d’inconnue d e liaison. 3 : une équation portant sur le moment en G1 des actions mécaniques de S2 sur S1. Cette équation sera projetée sur l’axe yo pour lequel n’apparaît pas d’inconnue d e liaison. Nous avons donc identifié en utilisant le graphe d’analyse, les sous-systèmes à isoler, les équations à écrire et les projections à effectuer. Ce travail est à faire préalablement à tout calcul.
3
3.1
PRINCIPE FONDAMENTAL DE LA DYNAMIQUE
Introduction : un peu d’histoire ème
Au début du 16 siècle, la notion de force restait obscure. Une force, à la différence d’un mouvement, est difficile à mettre en évidence. Dans le système de Ptolémée, par exemple, les trajectoires et les mouvements relatifs des planètes sont relativement bien décrits, quoique de façon complexe, par les systèmes d’épicycles, mais les lois qui régissent ces mouvements restent encore inconnues. Le principe d’inertie, énoncé par Galilée, puis le principe fondamental de la dynamique, énoncé par Newton, sont à la base de la mécanique classique. La dynamique étant le régime où dominent les effets d’inertie, par opposition à la statique. Le principe d'inertie: Ce principe est peut-être le plus grand apport de Galilée (1564-1642) à la physique. En faisant des expériences avec des billes qui roulent sur des plans de différentes natures, il observe que si le plan est très rugueux, la bille s'arrête rapidement, par contre, si le plan est très lisse ou recouvert d'huile par exemple, la bille parcourt une distance beaucoup plus grande avant de s'arrêter. Galilée eut alors l'idée de forces de frottement : le plan rugueux frotte très fortement sur la bille et l 'oblige à s'arrêter rapidement, en revanche, sur le plan lisse les forces de frottement sont très faibles et n'empêchent pas la bille de rouler. Dans la vie de tous les jours, les forces de frottement sont partout présentes et obligent les corps à stopper leur mouvement, c'est pour cela que pour entretenir ce mouvement on doit constamment appliquer une force extérieure à un corps pour contrebalancer ces forces de frottement. 77
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
Mais si on pouvait réduire ces forces de frottement à zéro, alors, le corps conserverait son mouvement indéfiniment. C'est en faisant une telle extrapolation que, Galilée donne une première formulation du principe d'inertie : « Tout corps possède une certaine "inertie" qui l'oblige à conserver sa vitesse, à moins qu'une force extérieure, une force de frottement par exemple, ne l'oblige à arrêter ce mouvement, i.e. à modifier cette vitesse ». C'est finalement à partir de ce principe que Newton (1642-1727) introduit en 1687 la description des causes du mouvement et introduira la notion fondamentale de force: le mouvement d'un système isolé étant rectiligne et uniforme, toute déviation par rapport à cette loi du mouvement est caractérisée par une force. En l'absence de force, le corps poursuit sa trajectoire et conserve sa vitesse. Il en découle d’ailleurs que l’étude de petites déviations à une loi du mouvement d’un corps permet l’étude des forces auxquelles ce corps est soumis, ce qui sera à la base des principes variationnels. Le principe fondamental de la dynamique: La formalisation des principes de la mécanique est publiée dans l’ouvrage principal de I. Newton: les « Philosophiae naturalis principia mathematica » (1687), où sont définies les notions de masse, de quantité de mouvement, de forces (motrice, centripète). Newton introduit un espace et un temps absolus, et énonce les principes fondamentaux de l'inertie, de l'égalité des actions réciproques ainsi que la «relation fondamentale de la dynamique» exprimant chez lui l'égalité entre force et dérivée de la quantité de mouvement
3.2
Enoncé du principe fondamental de la dynamique Il existe au moins un repère Rg appelé repère Galiléen, tel que pour tout système matériel Σ en mouvement par rapport à Rg, le torseur dynamique de Σ dans son mouvement par rapport à Rg soit égal au torseur des actions mécaniques extérieures à Σ. Si l’on note Σ l’extérieur de Σ le principe fondamental de la dynamique s’écrit alors :
{D (Σ / R g )}= {F (Σ → Σ)}
3.3
Conséquences En écrivant qu’en tout point de l’espace, les deux torseurs intervenant dans le principe fondamental de la dynamique sont égaux, on obtient deux équations vectorielles appelées théorèmes généraux de la dynamique. Soit m la masse et G le centre d’inertie du système matériel S en mouvement par rapport au référentiel galiléen Rg. On pose en un point A quelconque : r r
m Γ(G / R g ) r δ A (Σ / R g ) A
R (Σ → Σ ) r M (Σ → Σ ) A A
{D (Σ / R g )}=
◊
{F (Σ → Σ )}=
Théorème de la résultante dynamique Alors, la résultante du torseur dynamique de Σ dans son mouvement par rapport à Rg est égale à la résultante du torseur des actions mécaniques extérieures à Σ. r r
(
m Γ G / R g
) = R(Σ → Σ)
78
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
◊
Théorème du moment dynamique Par ailleurs, le moment en tout point A du torseur dynamique de Σ dans son mouvement par rapport à Rg est égal au moment du torseur des actions mécaniques extérieures à Σ au même point. r r
δ A
◊
(Σ / R g ) = M A (Σ → Σ )
Théorème des actions mutuelles Soit un système matériel Σ constitué de deux sous partie Σ1 et Σ2. On peut appliquer successivement le principe fondamental de la dynamique aux deux sous parties Σ1 et Σ2 de Σ. Alors :
D Σ1 / R g + D Σ2 / R g = D Σ / R g = {F (Σ → Σ )} = {F (Σ → Σ1)} + {F (Σ → Σ2 )} Par ailleurs :
D Σ1 / R g = {F (Σ1 → Σ1)} = {F (Σ → Σ1)} + {F (Σ2 → Σ1)}
Et :
D Σ2 / R g = {F (Σ 2 → Σ 2)} = {F (Σ → Σ2 )} + {F (Σ1 → Σ2 )}
En recombinant ces équations, on en déduit :
{F (Σ1 → Σ2)} = −{F (Σ2 → Σ1)}
◊
Cas particulier de la statique Le principe fondamental de la dynamique se ramène à celui de la statique lorsque le torseur dynamique du système matériel Σ par rapport à Rg est nul. Il est nul en particulier dans les trois cas suivants : • • •
◊
Les effets d’inertie sur le système matériel Σ sont négligés La masse du système matériel Σ est supposée nulle Le système matériel Σ est animé d’un mouvement de translation rectiligne uniforme par rapport au repère galiléen Rg.
Equations du mouvement L’application du principe fondamental de la dynamique conduit à écrire deux équations vectorielles dont les projections sur une base orthonormée donnent au maximum six équations scalaires indépendantes. Ces équations scalaires sont des équations différentielles du second ordre, en général non linéaires de la forme : &&i (t )) = 0 , où les paramètres qi sont les i paramètres du mouvement. f (q i (t ), q& i (t ), q
Outre les qi, ces équations peuvent également contenir des composantes inconnues d’actions mécaniques. Par définition, une équation du mouvement est une équation différentielle du second ordre des paramètres du mouvement ne contenant aucune composante inconnue d’action mécanique. Par définition une intégrale première du mouvement est une équation différentielle du premier ordre de la forme f (q i (t ), q& i (t )) = Constante , obtenue par intégration d’un équation du mouvement.
3.4
Référentiels Galiléens/non Galiléens En mécanique classique le temps est le même pour tous les observateurs, il est donné par des mouvements particuliers ou des processus naturels appelés horloges (oscillation d’un quartz, mouvement de certains astres, désintégration de la matière …). On parle alors de Chronologie Galiléenne.
79
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
◊
Référentiels Galiléens On montre que tout repère R en translation rectiligne uniforme par rapport à un repère Galiléen Rg est également un repère Galiléen. Le choix d’un repère Galiléen est fonction du problème posé. Un repère Galiléen est un repère dans lequel le principe fondamental de la dynamique est vérifié avec une bonne approximation, pour une étude donnée.
-
Repère de Copernic Ce repère est défini par le centre d’inertie du système solaire (sensiblement le centre du soleil) et par les directions stellaires. Il constitue une excellente approximation d’un repère galiléen pour l’étude des fusées interplanétaires par exemple.
-
Repère défini par le centre de la terre et les directions stellaires Ce repère constitue une excellente approximation d’un repère galiléen pour l’étude du mouvement de systèmes matériels restant dans le domaine terrestre, à l’intérieur d’une sphère centrée sur le centre de la terre et de rayon maximal égal à quatre ou cinq fois celui de la terre.
-
Repère terrestre Pour la plupart des problèmes terrestres, un repère lié à la terre constitue une très bonne approximation d’un repère galiléen.
◊
Référentiel non Galiléen Supposons un repère R ayant un mouvement quelconque mais connu par rapport à un repère galiléen Rg. Le principe fondamental de la dynamique, appliqué au système matériel Σ dans son mouvement par rapport au repère galiléen Rg s’écrit :
D Σ / R g = {F (Σ → Σ )} Le torseur dynamique a pour éléments de réduction en un point A quelconque : r Γ( P / Rg ) dm P ∈Σ = r ( ) AP P / Rg dm ∧ Γ P ∈Σ
{D (Σ / R g )}
∫
∫
On utilise alors la relation de composition des vecteurs accélération : r
r
r
r
r
Γ( P / Rg ) = Γ( P / R) + Γ( P ∈ R / Rg ) + 2.Ω ( R / Rg ) ∧ V ( P / R ) Ceci permet de construire le torseur « des forces d’inertie d’entraînement » de mouvement par rapport à R et Rg.
son
r Γ( P ∈ R / Rg ) dm P ∈Σ = − r ( ) AP P R / Rg dm ∧ Γ ∈ P ∈Σ
∫
{D (Σ ∈ R / R g )} ie
Σ dans
∫
En utilisant les formules issues du principe de conservation de la masse : r
m.Γ(G / Rg ) r r = − d [ ] ( ) ( ) ( ) AG m . G / Rg J . R / Rg ∧ Γ + Σ Ω G dt
{D (Σ ∈ R / R g )} ie
On construit de même le torseur « des forces d’inertie de Coriolis » de Σ dans son mouvement par rapport à R et Rg.
80
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
r r 2.Ω ( R / Rg ) ∧ V ( P / R ) dm P ∈Σ = − r r AP ∧ [2.Ω ( R / Rg ) ∧ V ( P / R )]dm P ∈Σ
{D (Σ ∈ R / R g )} ic
∫
∫
Avec ces conventions le principe fondamental de la dynamique s’écrit alors dans le repère nongaliléen R :
{D (Σ / R )} = {F (Σ → Σ)} + D ie Σ ∈ R / R g + D ic Σ ∈ R / R g
◊
Exemple : Accélération de la pesanteur. Considérons un fil à plomb en équilibre à la surface de la terre, matérialisant la verticale en un lieu donné. On choisi un repère galiléen Rg (O,xg,yg,zg), centré sur le centre de la terre et dont les axes pointent vers des directions stellaires. Un second repère R(O,x,y,zg) est attaché à la terre et en rotation uniforme par rapport à Rg autour de zg à la vitesse ω (un tour en 24 heures). On pose θ=(xg,x) avec θ=ω.t. Le fil à plomb, constitué d’un fil de masse négligeable et d’une sphère S, de masse m, de centre r
r r
d’inertie G est situé dans le plan méridien (O,z g,x). On pose : OG = r n , α = (n, x ) . n est un vecteur unitaire du plan méridien (O,zg,x) normal à la terre et dirigé vers le ciel, t est un vecteur unitaire du plan méridien (O,zg,x) normal à n, α est la latitude du point G. Le fil à plomb est en équilibre dans le repère R soit :
r
r
Ω(S / Rg ) = ω . z g . On peut donc considérer S et R comme équivalents pour la dérivation.
La sphère S est soumise à deux actions mécaniques, représentables au point G par deux forces, T : la tension du fil, de même direction que le fil, Fa : la force d’attraction Newtonienne de la terre, telle que : r mM r F a = − K n , Où m est la masse de la sphère, M, la masse de la terre (M= 5,9742 10 24 kg), K la 2 r constante de gravitation universelle (K=6.67 10-11) et r la distance entre les centres d’inertie de la terre et de la sphère. Le torseur dynamique de la sphère S au point G vaut : r
m.Γ (G ∈ S / Rg ) r r r = d [ G (S ).Ω(S / Rg )]= 0 δ G (S / Rg ) = dt J
{D (S / R g )}
81
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
Par ailleurs les actions mécaniques appliquées au point G se réduisent, à la tension du fil T et à la force d’attraction gravitationnelle Fa. Soit : r r
{F (S → S )} = T +r F a
0
En appliquant maintenant le principe de la résultante dynamique : r r r
m.Γ(G ∈ S / Rg ) = T + F a
Par conséquent : r
r
r
T = m.Γ(G ∈ S / Rg ) − F a
r
= m.Γ(G ∈ S / Rg ) + K
mM r n r 2
Dans le repère terrestre, le fil à plomb est en équilibre sous l’action de deux forces, la tension du fil et le poids de la sphère P=mg. Ceci permet d’en déduire l’expression de l’accélération de la pesanteur dans le repère terrestre : r
r
r
r
r
T + P = T + m g . =0.
fl
r M r r g = − K n − Γ(G ∈ S / Rg ) r 2
En utilisant maintenant la formule de changement de point du vecteur accélération :
→ r → r r r d r Γ(G ∈ S / Rg ) = Γ(O ∈ S / Rg ) + Ω ( S / Rg ) ∧ OG + Ω ( S / Rg ) ∧ Ω (S / Rg ) ∧ OG dt R On en déduit : r
r
r
r
r
Γ(G ∈ S / Rg ) = ω . z g ∧ [ω . z g ∧ r .n ] = −r .ω 2 . cos α . x Soit :
r r r r r M r M g = − K n + r .ω 2 . cos α . x ⇒ g = − K r .ω 2 . cos 2 α .n − r .ω 2 . cos α . sin α .t + r 2 r 2
Ainsi, le vecteur accélération de la pesanteur, relatif au repère terrestre n’est-il pas exactement dirigé vers le centre de la terre sauf à l’équateur et aux pôl es. La composante tangentielle de l’accélération de la pesanteur, gt, reste inférieure à 1% de sa composante normale en deçà de 5000 km d’altitude.
82
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
4
4.1
PRINCIPE DES PUISSANCES VIRTUELLES
Introduction : un peu d’histoire. En 1743, d’Alembert (1717-1783) publie le Traité de dynamique, dans lequel il expose le principe « de la conservation des forces vives » aujourd’hui appelé principe de d’Alembert : L’énoncé du principe de d’Alembert est le suivant : « Si l'on imagine un système de corps en mouvement, liés entre eux d'une manière quelconque, et réagissant les uns sur les autres au moyen de ces liaisons, de manière à modifier les mouvements que chaque corps isolé prendrait en vertu des seules forces qui l'animent, on pourra considérer ces mouvements comme composés : 1° des mouvements que les corps prennent effectivement, en vertu des forces qui les animent séparément, combinées avec les réactions du système 2° d'autres mouvements qui sont détruits par suite des liaisons du système. d'où il résulte que les mouvements ainsi détruits doivent être tels, que les corps animés de ces seuls mouvements se feraient équilibre au moyen des liaisons du système ». On pourra noter que par rapport au principe fondamental de la dynamique de Newton, l’objet du traité de dynamique de d’Alembert est le comportement de systèmes de solides présentant des liaisons ou en contact mutuel. En 1788, Lagrange publie la « Mécanique analytique », qui est fondé sur un nouveau principe variationnel, issu de l'association du principe des vitesses virtuelles et de celui de d'Alembert. Lagrange écrit que le travail total des forces appliquées et des forces d'inertie est nul pour tout déplacement virtuel compatible avec les liaisons. Lagrange aboutit ainsi, pour les systèmes conservatifs, à un système d'équations différentielles universel liant les dérivées de la «force vive» (assimilable à notre énergie cinétique) et celles de la fonction de force (qui joue le même rôle que l'énergie potentielle).
4.2
Enoncé du principe des puissances virtuelles ou PPV. A tout instant t, pour tout système matériel S, et pour tout champ de vitesse virtuelle, la puissance virtuelle des quantités d’accélération galiléenne est égale à la somme de la puissance virtuelle des actions extérieures et de la puissance virtuelle des actions intérieures à S. En mécanique des solides indéformables, on utilise des champs de vitesses virtuelles dont la restriction sur chacun des solides est le champ des moments d’un torseur, afin de respecter la rigidité de chacun de ces solides. En conséquence, la puissance virtuelle des actions intérieure à chacun des solides d’un système de solides est nulle, la puissance virtuelle des efforts intérieurs se réduit donc à la puissance virtuelle des interactions entre les solides du système. Ce qui se note pour un système de solides Σ, et un champ de vitesses virtuelles quelconque ∗ V dont la restriction sur chacun des solides Si constituant Σ est le champ des moments d’un torseur noté r
∗ V i : n
∑= i 1
* D (S i / Rg )×V i
n
=
∑=
*
P
n
(Σ → S i , Rg ) + ∑ P * (S i ↔ S j ) i =1
i 1
83
r
∀V ∗
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
4.3 4.3.1
Choix de torseurs virtuels particuliers et théorèmes de la dynamique Torseur global quelconque : Equivalence du principe des puissances virtuelles et du principe fondamental de la dynamique r Si le champ des vitesses virtuelles V ∗ est en totalité le champ des moments d’un torseur V ∗ quelconque, le principe des puissances virtuelles conduit à :
D (Σ / Rg )×V ∗ = F (Σ → Σ )×V ∗
∀V ∗
On retrouve ainsi le principe fondamental de la dynamique sous la forme :
D (Σ / Rg ) = F (Σ → Σ )
4.3.2
◊
Torseur des vitesses galiléennes : Théorème de l’énergie cinétique Théorème de l’énergie cinétique pour un solide unique. r Si le champ des vitesses virtuelles V ∗ est le champ des moments du torseur des vitesses galiléennes V S j /Rg , le principe des puissances virtuelles pour un solide Si, conduit à :
D (S i / Rg )× V (S i / Rg ) = F S i → S i × V (S i / Rg ) r r d 1 r Γ( P / Rg )dm( P ) ⋅ V ( P / Rg ) = V ( P / Rg ) 2 dm( P ) dt 2 Rg S S
Or : D (S i / Rg ) ×V (S i / Rg ) =
∫
∫
i
i
En appliquant le principe de conservation de la masse :
D (Si / Rg )×V (S i / Rg ) =
d 1 r 2 V ( P / Rg ) dm( P ) dt 2
∫
S i
=
dT (S i / Rg ) dt
Rg
On retrouve ainsi le théorème de l’énergie cinétique pour un solide Si : dT (S i / Rg ) dt
◊
= P(S i → S i / Rg )
Théorème de l’énergie cinétique pour un système de solides. Si maintenant on considère un système de solide Σ, constitué de n solides Si, il suffit de faire la somme des équations issues du théorème de l’énergie pour chacun des n solides, soit : dT (Σ / Rg ) dt dT (Σ / Rg ) dt
n
=
∑= i 1
dT (S i / Rg ) dt
n
=
i
i
i
i 1
F (Σ → S i ) + F (S j → S i ) ×V (S i / Rg ) i =1 j ≠i n
=
∑= F (S → S )×V (S / Rg )
∑
∑
En rassemblant terme à terme : dT (Σ / Rg ) dt
n
=
∑= F (Σ → S )×V (S i
i 1
n
i
/ Rg ) +
∑= ∑> F (S
j
→ S i )×V (S i
i 1 j i
/ Rg ) +
∑> F (S → S )×V (S i
j i
Soit :
84
j
j
)
/ Rg
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
dT (Σ / Rg ) dt
n
=
n
F (S ∑= F (Σ → S )×V (S / Rg ) + ∑∑ = > i
i
j
i 1
→ S i )× [V (S i / Rg ) − V (S j / Rg )]
i 1 j i
Ce qui donne finalement : dT (Σ / Rg ) dt
n
=
∑=
n
(
F Σ → S i
)×V (Si / Rg ) + ∑∑ F (S j → S i )×V (S i / S j ) i =1 j >i
i 1
Le théorème de l’énergie cinétique pour un système de solides s’écrit donc : dT (Σ / Rg ) dt
4.3.3
◊
(
=
)
P Σ → Σ / Rg 14 4 244 3
Puissance des efforts extérieurs
(
+
P S i
↔ S j )
14243
Puissance des efforts intérieurs
Torseurs de Lagrange : équations de Lagrange Définition des torseurs de Lagrange : Supposons que la position d’un solide S dans Rg soit paramétrée par n paramètres qi et la variable temps t . Puisque le champ des vitesses de S est le champ des moments d’un torseur il vérifie :
∀( P , Q ) ∈ S
r
r
r
V (Q / Rg ) = V ( P / Rg ) + Ω(S / Rg ) ∧ PQ
Alors quel que soit le paramètre du mouvement qi :
∀( P , Q ) ∈ S
r ∂ r ∂ r ∂ r ∂ Ω(S / Rg ) ∧ PQ + Ω(S / Rg ) ∧ V (Q / Rg ) = V ( P / Rg ) + PQ ∂q& i ∂q& i ∂q& i ∂q& i
Le dernier terme est nul donc :
∂ r ∂ r ∂ r Ω(S / Rg ) ∧ PQ V (Q / Rg ) = V ( P / Rg ) + ∂q& i ∂q& i ∂q& i
∀( P , Q ) ∈ S
r
∂V Ainsi, les champs sont également des torseurs, appelés torseurs de Lagrange et notés ∂q& i V qi(S/Rg).
◊ -
Conséquences Conséquence N°1 : Comme les paramètres qi sont les paramètres du mouvement : r
∀ P ∈ S V ( P / Rg ) =
∂OP dqi ∂ qi Rg dt i =1 n
∑
∂OP + ∂t Rg
En dérivant une seconde fois, les paramètres ( q i , q& i ,t) étant indépendants, on trouve :
∀ P ∈ S
∂OP ∂ r V ( P / Rg ) = ∂q& i ∂qi
Rg
-
Conséquence N°2 : Le produit du torseur dynamique du solide S par un torseur de Lagrange est noté Pqi(S/Rg) et vaut: P q (S / Rg ) = D S / R g i
×V qi (S / Rg )
Compte tenu de ce qui vient d’être montré :
85
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
r ∂OP Γ( P / Rg ) ⋅ P qi (S / Rg ) = .dm( P ) ∂qi P ∈S
∫
Or :
∂ ∂q& i
∂OP 1 r r V ( P / Rg )2 = V ( P / Rg ) ⋅ ∂qi 2
Et par suite : d r d ∂OP ∂ 1 r ∂OP r + V ( P / Rg ) ⋅ V ( P / Rg )2 = [V ( P / Rg )] Rg ⋅ dt ∂q& i 2 dt ∂q i ∂qi Rg dt Rg d
Soit : r ∂ 1 r ∂OP r d ∂OP 2 = Γ ⋅ + ⋅ ( ) ( ) ( ) V P / Rg P / Rg V P / Rg ∂qi dt ∂q&i 2 dt ∂qi Rg
d
Rg
Par ailleurs : r
∂ 1 r ∂V ( P / Rg ) r d ∂OP r = V ( P / Rg ) ⋅ V ( P / Rg )2 = V ( P / Rg )⋅ ∂qi 2 ∂qi dt ∂qi
Rg
On en déduit donc :
r d ∂ ∂ ∂OP Γ( P / Rg ) ⋅ P qi (S / Rg ) = .dm( P ) = .T (S / Rg ) − ∂qi dt ∂q& i ∂qi P ∈S
∫
◊
Application du PPV à un unique solide S. Si, pour écrire le principe des puissances virtuelles, on choisit comme champ de vitesses virtuelles ∗ V le champ des moments du torseur de Lagrange V qi(S/Rg), on obtient alors pour chacun des paramètres qi, l’équation suivante : r
D (S / R g )×V qi (S / Rg ) = F (S → S )×V qi (S / Rg )
(
)
Le second membre est alors noté Qqi S → S / Rg et est appelée la ième composante des efforts extérieurs généralisés sur S. Le principe des puissances virtuelles s’écrit alors sous la forme d’équations de Lagrange :
(
)
P qi (S / Rg ) = Qqi S → S / Rg
d ∂ ∂ .T (S / Rg ) et Qqi (S → S / Rg ) = F (S → S )×V qi (S / Rg ) − dt ∂q& i ∂qi
Avec P qi (S / Rg ) =
On note qu’on obtient autant d’équations de Lagrange qu’il y a de paramètres du mouvement, ce qui permet de résoudre le problème.
◊
Application du PPV à un système de solides Σ : On écrit alors que la puissance virtuelle des quantité d’accélération est égale à la puissance virtuelle des efforts extérieurs plus la puissance virtuelle des efforts intérieurs. Soit : P qi (S / Rg ) = Qqi
(Σ → Σ / Rg ) + Qq (int Σ ) i
Comme nous l’avons fait pour l’énergie cinétique, on trouve la puissance virtuelle des efforts intérieurs en écrivant le PPV pour chacun des solides puis en regroupant les termes deux à deux, alors : 86
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES n
Qqi (int Σ ) =
F (S ∑∑ = >
j
→ S i )× [V qi (Si / Rg ) − V qi (S j / Rg )] =
i 1 j i
◊ -
n
F (S ∑∑ = >
j
→ S i )×V qi (S i / S j )
i 1 j i
Fonction de force. Définition On dit que les actions mécaniques (*) d’une partie Se extérieure à Σ sur une partie Si intérieure à Σ dérivent d’une fonction de force si il existe une fonction numérique U des variables (q ,t i ) telle que :
∀i = 1, n Qqi (S e → S i / Rg ) =
∂U (S e → S i / Rg ) ∂qi
L’existence d’une fonction de force simplifie considérablement le calcul de Qqi.
-
Relation entre fonction de force et énergie potentielle On rappelle peut associer une énergie potentielle aux actions extérieures (exercées par un solide Se) sur un solide Si en mouvement par rapport à R, si il existe une fonction scalaire Ep telle que : P (S e Où P (S e
→ S i / R ) = −
d
E p (S e dt
→ S i / R )
→ S i / R ) est la puissance des actions extérieures à Si (exercées par Se) sur Si.
En développant par rapport aux paramètres du mouvement :
n ∂ E p (S e → S i / R ) ∂ E p & P (S e → S i / R ) = − E p (S e → S i / R ) = − qi + ∂qi ∂t dt i =1
∑
d
L’existence d’une énergie potentielle E p (S e
→ S i / R ) entraîne l’existence d’une fonction de force U (S e → S i / R ) , telle que U (S e → S i / R ) = − E p (S e → S i / R ) . Attention la réciproque n’est pas vraie !
-
Illustration On étudie un ressort de traction compression R (raideur k, longueur libre l o, masse négligeable) monté entre deux solides, T le cylindre, et M la masse. La liaison T-M est une glissière. On suppose que T est en mouvement par rapport à Ro(O,xo,yo,zo) d’un mouvement rectiligne vibratoire connu : r
r
OH = h. sinω t . y o . La position de M dans Ro est repérée par le paramètre y tel que : OP = y y o
r
Alors :
0 {V ( M / Ro)} = r , y& y o
r
Et :
0 {V ( M / T )} = r ( y& − hω cos ω t ) y o r
Par ailleurs :
− k ( y − h sinr ω t − l o ) y o {F (r → M )} = 0
87
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
Donc: P (T ← r → M ) = {F (r → M )}× ({V ( M / Ro)} − {V (T / Ro)}) = {F (r → M )}× {V ( M / T )} Soit : P (T ← r → M ) = −k ( y − h sin ω t − l o )( y& − hω cos ω t ) =
d k
− ( y − h sin ω t − l o ) 2 dt 2 Ro
Il existe donc une énergie potentielle E p pour les inter-efforts entre T et M par l’intermédiaire du ressort r. On a donc : E p (T ← r → M ) =
k 2
( y − h sin ω t − l o )2
et
U (T ← r → M ) = −
k 2
( y − h sin ω t − l o )2
Alors : Q y (T ← r → M ) =
∂U (T ← r → M ) = −k ( y − h sin ω t − l o ) ∂ y
On cherche maintenant à déterminer les Qqi associés aux efforts extérieurs exercés sur M par l’intermédiaire du ressort r : P (r → M / Ro ) = {F (r → M )}× {V ( M / Ro )} Soit
P (r → M / Ro ) = − k ( y − h sin ω t − l o ) y&
Il n’existe pas d’énergie potentielle ici telle q ue P (r → M / Ro ) = −
d
E p (r → M / Ro) dt
On calcule donc directement Qy, pour cela il faut calculer le torseur de Lagrange : r
0 {V ( M / Ro )} = r y& y o
r
Donc
{V y ( M / Ro )} = ∂∂& {V ( M / Ro)} = r0 y y o
Alors : Q y (r → M / Ro ) = F (r → M ) ×V y ( M / Ro ) Soit : Q y (r → M / Ro ) = − k ( y − h sin ω t − l o ) =
1 − k ( y − h sin ω t − l o ) dy 2 d
Donc, bien qu’il n’existe pas d’énergie potentielle E p (r → M / Ro ) dans ce cas, il existe néanmoins une fonction de force U (r → M / Ro ) . Si l’on suppose maintenant que le mouvement de T par rapport à M est libre et paramétré par λ tel r
r
que : OH = λ . y o , tandis que le mouvement de M par rapport à T est paramétré par γ tel que : HP = l . y o r
0 {V ( M / Ro)} = & & r (λ + l ) y o
Alors :
r
Donc
0 {V λ ( M / Ro )} = {V l ( M / Ro )} = r y o
r
Par ailleurs :
− k (l −rl o ) y o {F (r → M )} = 0
Alors :
Qλ (r → M / Ro ) = Ql (r → M / Ro ) = − k (l − l o )
Donc :
∂ Q (r → M / Ro ) = − k ∂l λ
Donc :
et
∂ Q (r → M / Ro) = 0 ∂λ l
∂ ∂ Qλ (r → M / Ro ) ≠ Q (r → M / Ro ) ∂l ∂λ l
Donc la fonction de force U (r → M / Ro ) n’existe pas. Cet exemple montre qu’il faut être très rigoureux dans la définition des grandeurs à déterminer et que cela implique des notations complètes.
-
Action de la pesanteur
88
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
Dans un repère R(O,x,y,z) lié à la terre. On considère que l’accélération de la pesanteur g est r r dirigée selon la verticale à la surface (O,z) (ce qui est une approximation) et dirigée vers le bas. g = − g . z . Alors, il existe une énergie potentielle associée à l’action de la pesanteur sur le solide S de masse m et de centre d’inertie G, dans son mouvement par rapport à la terre : r
g r E p pes → S / R = − m g ⋅ OG Il existe donc également une fonction de force U associée : r
g r U pes → S / R = m g ⋅ OG Si le mouvement de G par rapport à R est paramétré par i paramètres q i alors le calcul des Qqi est immédiat : r
∀i = 1, n
-
g pes → S / R ∂ U r g r ∂ OG Q q pes → S / R = m g ⋅ = i ∂q i ∂q i
Ressort de traction-torsion Deux solides S1 et S2 sont liés par une liaison pivot glissant dans laquelle est intercalée un ressort de traction-torsion de raideur K en traction-compression, de raideur C en torsion, de longueur à vide Lo et de torsion à vide θo. Le paramètre L mesure la distance entre les deux solides le long de l’axe de la liaison et le paramètre θ mesure la rotation du solide S2 par rapport au solide S1. L’ énergie potentielle pour les actions mutuelles de S1 et S2 s’écrit : 1 E p (S 1 ← r → S 2) = K ( L − Lo )2 2
+
1 2
C (θ − θ o ) 2
On en déduit la fonction de force pour les inter-efforts entre S1 et S2 1 U (S 1 ← r → S 2) = − K ( L − Lo )2 2
−
1 2
C (θ − θ o ) 2
Les mouvements relatifs de S1 et S2 sont paramétrés ici par les deux paramètres L et θ, donc :
4.3.4
◊
Q L (S 1 ← r → S 2) =
∂U (S 1 ← r → S 2) = − K ( L − Lo ) ∂ L
Qθ (S 1 ← r → S 2 ) =
∂U (S 1 ← r → S 2) = −C (θ − θ o ) ∂θ
Combinaison des torseurs de Lagrange : Equation de Painlevé Définition du torseur de Painlevé:
∂OP dqi comme champ de vitesses virtuelles pour écrire le ∂ q dt i Rg i =1 n
On choisit le champ
∑
principe des puissances virtuelles. Ce champ est le champ des moments d’un torseur, appelé torseur de Painlevé, étant donné qu’il est une combinaison linéaire de torseur de Lagrange. r
On rappelle que V ( P / Rg ) =
∂OP dqi ∂ qi Rg dt i =1 n
∑
∂OP + . ∂t Rg
89
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
Donc le torseur de Painlevé est distinct du torseur des vitesses galiléennes. L’équation de Painlevé est donc aussi distincte du théorème de l’énergie cinétique, à condition que le mouvement du système dépende explicitement du temps. Autrement dit, si le système contient des actionneurs qui imposent partiellement la cinématique du système.
◊ -
Conséquences Premier membre de l’équation de Painlevé. Pour chacun des torseurs de Lagrange on a :
d ∂ ∂ .T (Σ / Rg ) − dt ∂q& i ∂q i
P qi (Σ / Rg ) =
Donc pour le torseur de Painlevé, le premier membre du principe des puissances virtuelles s’écrit :
d ∂ ∂ dq i & − .T (Σ / Rg ) dt dt ∂q i ∂q i i =1 n
P P (Σ / Rg ) =
∑
Si l’on développe : n
P P (Σ / Rg ) =
∑= dt d ∂∂q& i 1
dq i
T (Σ / Rg )
i
dt
−
dq i dt
∂ T (Σ / Rg ) ∂q i
L’énergie cinétique se développe comme suit vis-à-vis des paramètres qi :
n ∂OP dq i ∂OP + T (Σ / Rg ) = V 2 ( P / Rg ) dm = 2 2 q i dt t ∂ ∂ Σ Σ i =1 1
∫
r
1
∫∑
2
dm
Si l’on développe : 2 n ∂OP ∂OP dq dq j n ∂OP ∂OP dq i ∂OP 1 i dm + dm + ⋅ ⋅ T (Σ / Rg ) = dm ∂ ∂ ∂ 2 q i ∂q j dt dt qi t dt 2 ∂t 4 3 i , j =1 i =1 3 1Σ44 24 144 4 4 244 4 4 1Σ44 4 4 4 4 244 4 4 4 4 3 Σ
1
∫∑
∫ ∑
∫
T 1(Σ / Rg )
T 2(Σ / Rg )
T 0(Σ / Rg )
Avec ces notations on peut alors écrire : P P (Σ / Rg ) =
-
d dt
[T 2(Σ / Rg ) − T 0(Σ / Rg )] +
∂ T (Σ / Rg ) ∂t
Second membre de l’équation de Painlevé. Le second membre ne présente aucune particularité. Au final, l’équation de Painlevé s’écrit : n dq ∂ [T 2(Σ / Rg ) − T 0(Σ / Rg )] + T (Σ / Rg ) = (Qqi (Σ → Σ / Rg ) + Qqi (int Σ )) i ∂t dt dt i =1
d
◊
∑
Intégrale première de Painlevé L’intérêt principal de l’équation de Painlevé est de fournir dans certaines conditions une intégrale première du mouvement : • • • •
toutes les liaisons géométriques sont parfaites toutes les liaisons géométriques extérieures le sont avec des solides ayant un mouvement imposé par rapport à Rg. Toutes les actions mécaniques dérivent d’une fonction de force U T(S/Rg) et U ne dépendent pas explicitement du temps. 90
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
Alors : T 2(Σ / Rg ) − T 0(Σ / Rg ) = U + constante
91
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
5
A RETENIR
•
Torseur des actions des Σi sur Σ au point A :
ΨΣi →Σ .dP Σ F {Σ i → Σ}= AP ∧ ΨΣi →Σ .dP A Σ
∫
∫
•
Principe fondamental de la dynamique : il existe au moins un repère Rg appelé repère Galiléen, tel que pour tout système matériel Σ en mouvement par rapport à Rg :
D Σ / R g = {F (Σ → Σ )} •
Dans un référentiel non Galiléen R en mouvement (connu) par rapport à Rg, le principe fondamental de la dynamique s’écrit :
{D (Σ / R )} = {F (Σ → Σ )} + D ie Σ ∈ R / R g + D ic Σ ∈ R / R g Avec : r Γ( P ∈ R / Rg ) dm P ∈ Σ = − r AP ∧ Γ( P ∈ R / Rg ) dm P ∈Σ
∫
{D (Σ ∈ R / R g )} ie
∫
r r 2.Ω ( R / Rg ) ∧ V ( P / R ) dm P ∈Σ = − r r AP ∧ [2.Ω ( R / Rg ) ∧ V ( P / R )]dm P ∈Σ
∫
{D (Σ ∈ R / R g )} ic
•
∫
Principe des puissances virtuelles (PPV) : pour un système de solides Σ, et un champ de vitesses r virtuelles V ∗ dont la restriction sur chacun des solides Si constituant Σ est le champ des ∗ moments d’un torseur noté V i quelconque s’écrit : n
∑=
* D (S i / Rg )×V i
i 1
•
=
∑=
*
P
i 1
n
(Σ → S i , Rg ) + ∑ P * (S i ↔ S j )
r
∀V ∗
i =1
r PPV avec V ∗ égal au champ des vitesses Galiléennes : théorème de l’énergie cinétique
dT (Σ / Rg ) dt
•
n
=
(
)
P Σ → Σ / Rg
(
+
14 4 244 3
Puissance des efforts extérieurs
P S i
↔ S j )
14243
Puissance des efforts intérieurs
r PPV avec V ∗ égal au champ des moments de l’un quelconque des torseurs de Lagrange : équation de Lagrange
(
)
P qi (S / Rg ) = Qqi S → S / Rg
d ∂ ∂ .T (S / Rg ) et − dt ∂q& i ∂q i
Avec P qi (S / Rg ) =
(
) = F (S → S )×V q (S / Rg )
Q qi S → S / Rg
92
i
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables PRINCIPE FONDAMENTAL DE LA DYNAMIQUE, PRINCIPE DES PUISSANCES VIRTUELLES
Avec V qi (S / Rg ) → •
∀( P , Q ) ∈ S
∂ r ∂ r ∂ r V (Q / Rg ) = V ( P / Rg ) + Ω(S / Rg ) ∧ PQ ∂q& i ∂q& i ∂q& i
Fonction de force : les actions mécaniques (*) d’une partie Se extérieure à Σ sur une partie Si intérieure à Σ dérivent d’une fonction de force, s’il existe une fonction numérique U des variables (q ,t i ) telle que :
∀i = 1, n
Q q (S e i
→ S i
/ Rg ) =
∂U (S e → S i ∂q i
/ Rg )
•
Quelques fonctions de force bien utiles :
•
g dm( P 1 )dm( P 2) Attraction gravitationnelle entre deux solides S1 et S2 : U S 2 ↔ S 1 = G o . P 1 P 2 S 1× S 2
•
r
∫
Attraction gravitationnelle à la surface de la terre, entre la terre et un solide S en mouvement par rapport à la terre : r
r
g g r U Terre ↔ S = U Pes → S / Terre = m(S ). g ⋅ OG
•
r 1 Ressort de traction-torsion : U S 1 ↔ S 2 = − K ( L − Lo )2 2
−
1 2
C (θ − θ o )2
Où r est un ressort de traction torsion, de raideur k en traction-compression, de raideur C en torsion, et de longueur libre lo, et où l paramètre la distance entre les points d’attache du ressort sur les deux solides S1 et S2 et de torsion libre θo où θ paramètre la torsion du ressort.
93
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
1
Introduction ______________________________________________________ 96
2
Méthodologie: ____________________________________________________ 96
3
Mise en forme du problème. _________________________________________ 97 3.1
Linéarisation_______________________________________________________ 97
3.1.1
◊ ◊ ◊ ◊ 3.1.2
◊ ◊ ◊ ◊ 3.2
Rappels : Résolution de systèmes linéaires d’équations différentielles_______ 102
3.2.1 3.2.2
◊ ◊ ◊ ◊ ◊ ◊ ◊ 4
Résolution numérique. ___________________________________________________ 102 Résolution analytique ____________________________________________________ 103 Cas N°1 : Le second membre F (t) est nul. ____________________________________ 104 Cas N°2 : Le second membre F (t) est nul et C est nulle. _________________________ 105 Cas N°3 : Le second membre F (t) est non-nul et C est nulle. _____________________ 106 Cas N°4 : Le second membre F (t) est nul et C est une combinaison linéaire de M et K. 106 Critère de Routh ________________________________________________________ 107 Rappel : Calcul du déterminant d’une matrice _________________________________ 108 Rappel : Matrice 2x2 valeurs propres, vecteurs propres __________________________ 109
Equilibre et stabilité_______________________________________________ 109 4.1
Introduction ______________________________________________________ 109
4.2
Systèmes conservatifs_______________________________________________ 109
◊ ◊ ◊ ◊ ◊ ◊ ◊ 4.3
Equilibre ______________________________________________________________ Stabilité : Définition _____________________________________________________ Stabilité : Théorème de Lejeune Dirichlet_______________________________ ______ Extension aux systèmes visqueux ___________________________________________ Equilibre et stabilité pour un système à un seul paramètre ________________________ Equilibre et stabilité pour un système à n paramètres ____________________________ Exemple, étouffeur de vibration ____________________________________________
110 110 110 111 112 112 113
Cas général, méthode directe ou méthode de Liapounov__________________ 114
◊ ◊ ◊ ◊ ◊ ◊ 5
Linéarisation des équations de Lagrange ______________________________________ 97 Premier membre des équations de Lagrange____________________________________ 97 Second membre des équations de Lagrange : Fonction de force_____________________ 98 Second membre des équations de Lagrange : Loi visqueuse________________________ 99 Système d’équations du mouvement après linéarisation : __________________________ 99 Exemple : L’étouffeur de vibrations __________________________________________ 99 Calcul des Qqi : _________________________________________________________ 100 Calcul des Pqi : __________________________________________________________ 101 Equations de Lagrange : __________________________________________________ 101 Equations linéarisées _____________________________________________________ 101
Etat de mouvement ______________________________________________________ Détermination de l’équilibre par la méthode directe _____________________________ Stabilité au sens de Liapounov _____________________________________________ Théorème de Liapounov __________________________________________________ Stabilité asymptotique ____________________________________________________ Stabilité orbitale ________________________________________________________
114 114 114 114 115 115
Vibrations autour d’une position d’équilibre stable______________________ 115
◊ ◊ ◊ ◊
Introduction ____________________________________________________________ Vibrations libres ________________________________________________________ Vibrations forcées _______________________________________________________ Vibrations amorties ______________________________________________________
115 116 116 116
6 Mécanique des Chocs-Percussions ___________________________________ 116 6.1
Introduction ______________________________________________________ 116
6.2
Cas d’un point matériel _____________________________________________ 117
94
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
6.3
Cas d’un solide ou d’un système de solides _____________________________ 117
◊ ◊ ◊ 6.3.2 6.3.3
◊ ◊ ◊
Remarque 1 : ___________________________________________________________ Remarque 2 : ___________________________________________________________ Remarque 3 : ___________________________________________________________ Percussion de liaison _____________________________________________________ Choc sans frottement entre deux solides______________________________________ Définition de e__________________________________________________________ Propriété de e et P : 0 ≤ e ≤ 1 et 0 ≤ P ______________________________________ Cas particuliers importants ________________________________________________
95
118 118 118 118 119 119 119 121
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
1
INTRODUCTION
Dans les chapitres précédents nous nous sommes attachés à décrire la cinématique retenue, et à exprimer les principes fondamentaux dans le formalisme de la mécanique des solides indéformables. L’écriture de ces principes permet de mettre en place un système d’équations différentielles contenant autant d’équations que d’inconnues. L’objectif de cette partie est d’exposer différentes méthodes permettant de résoudre les systèmes d’équations différentielles obtenus pour des applications particulières. On s’intéressera aux problèmes d’équilibre, de stabilité des équilibres et de vibrations autour des positions d’équilibre. Généralement, on procèdera d’abord à une linéarisation des équations puis on appliquera une procédure de résolution numérique. Cependant, on verra que des solutions analytiques peuvent être obtenues dans de nombreux cas. Dans un premier temps, nous allons rapidement résumer les chapitres précédents afin de dégager une méthodologie à suivre pour la mise en place des équations du système. Puis nous verrons ensuite comment résoudre les systèmes d’équations différentielles obtenus.
2
-
METHODOLOGIE:
A : Objectifs. •
-
Bien préciser les objectifs de l’étude afin de choisir une schématisation cinématique et un paramétrage adaptés au problème.
B : Cinématique. •
Tracer un graphe de structure
•
Déterminer la mobilité du système Réduire éventuellement le graphe de structure. Mettre en place le paramétrage et les figures de projections associées.
• •
Nombre d’équations du mouvement à trouver.
fl
-
C :Analyse. • •
Compléter le graphe de structure avec les actions mécaniques, tracer le graphe d’analyse. Compte tenu de l’objectif visé et en s’appuyant sur le graphe d’analyse choisir :
Choisir :
Le ou les sous systèmes à isoler
La ou les projections à effectuer.
Le ou les théorèmes fondamentaux à utiliser.
96
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL • •
Préciser clairement quelles équations doivent être écrites et pourquoi avant de les expliciter. N’effectuer les projections qu’au dernier moment et se contenter du minimum nécessaire, pour éviter des calculs fastidieux.
fl Système d’équations différentielles non-linéaires.
3
MISE EN FORME DU PROBLEME.
3.1
Linéarisation Dans le cas général, on ne sait pas trouver une solution analytique au système d’équations différentielles obtenu. Aussi utilise t’on une méthode approchée. On procède à une linéarisation au premier ordre du système d’équations différentielles obtenu afin de pouvoir le résoudre. Supposons qu’à l’instant t o on connaisse les paramètres du mouvement, on note q(t o ) le vecteur dont les composantes sont les n paramètres du mouvement à l’instant to et q(t) le vecteur dont les composantes sont les n paramètres du mouvement à l’instant t . On note ε la variation infinitésimale des paramètres du mouvements correspondant à une variation infinitésimale du temps dt , autour du point (t o, q(t o )) : q (t ) = q(t o
+ dt ) = q(t o ) + ε
alors
dq dt
=
d ε
d 2 q
et
dt
dt 2
=
d 2 ε dt 2
On suppose pour la linéarisation que ε est très petit et on suppose également que les dérivées successives par rapport au temps de ε restent aussi très petites. Cette dernière hypothèse n’est pas toujours vérifiée. Lors d’un choc en particulier, de très grandes variations de vitesses sont observées pour de très petites variations de position. D’autres techniques de résolution devront alors être employées.
3.1.1
Linéarisation des équations de Lagrange Nous allons voir ce qu’implique cette opération de linéarisation, lorsque les paramètres du mouvement d’un système de solides sont n paramètres du mouvement indépendants, et que les équations du mouvement sont données par les équations de Lagrange :
(
)
P qi (S / Rg ) = Qqi S → S / Rg
d ∂ ∂ .T (S / Rg ) et − dt ∂q& i ∂q i
Avec P qi (S / Rg ) =
◊
(
) = F (S → S )×V q (S / Rg )
Q qi S → S / Rg
i
Premier membre des équations de Lagrange On rappelle que, l’énergie cinétique d’un solide S, se développe comme suit en fonction des dérivées partielle vis-à-vis des paramètres du mouvement :
n ∂OP dq k ∂OP 2 T (S / Rg ) = V ( P / Rg ) dm = + q k dt t ∂ ∂ 2 2 S S k =1 1
∫
r
1
∫ ∑
2
dm
Si l’on développe et que l’on simplifie en considérant que les produits deux à deux des vitesses de variations des paramètres du mouvement sont nulles, puisqu’on linéarise, on obtient alors :
97
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL n
P qi (S / Rg ) =
∑=
j 1
d 2 q j dt 2
∫ S
n dq ∂OP ∂OP ∂OP ∂ 2 OP ∂OP ∂ 2 OP j dm + dm + ⋅ ⋅ ⋅ 2 dm 2 ∂q i ∂q j ∂ ∂ dt q i ∂q j ∂t qi ∂t j =1
∑ ∫
∫
S
S
Si le système ne contient aucun actionneur, alors la position du point P ne dépend pas explicitement du temps. Le premier membre des équations de Lagrange s’écrit alors comme suit : 2
n d q ∂OP r j = 0 ⇒ P qi (S / Rg ) = 2 ∂t j =1 dt
∑
∫ S
n ∂OP ∂OP dm = M ij ⋅ ∂q i ∂q j j =1
∑
d 2 q j dt 2
n
=
∑= M
ij
j 1
d 2 ε j dt 2
Compte tenu de la commutativité du p roduit scalaire : Mij=M ji Le premier membre du système formé par les n équations de Lagrange linéarisées, peut donc se noter dans ce cas particulier : M (q o , t ).
d 2 ε dt 2
Où M est une matrice carrée symétrique définie positive de dimension nxn, où n est la mobilité du système et ε un vecteur de dimension n, contenant les n variations des paramètres du mouvement autour d’une position qo. Dans le cas général, le premier membre du système formé par les équations de Lagrange linéarisées, se note de la façon suivante : M (q o , t )
d 2 ε 2
dt
+ C (q o , t )
d ε dt
+ F (q o , t )
Où M et C sont des matrices carrées de dimension nxn, où n est la mobilité du système et ε un vecteur de dimension n, contenant les n variations des paramètres du mouvement autour d’une position qo. M est symétrique, définie positive.
◊
Second membre des équations de Lagrange : Fonction de force Lorsque les actions mécaniques extérieure à S appliquées sur S du fait d’une loi (*), dérivent d’une fonction de force, il existe une fonction numérique U des variablesqi telle que :
(*) ∂U S → S / Rg (*) Q q S → S / Rg = i ∂q i La linéarisation par rapport aux n variations infinitésimales par rapport à la position qo s’écrit donc :
(*) ∂U S → S / Rg (*) Q qi S → S / Rg = ∂q i
(*) 2 U S → S / Rg ∂ m + .ε j ∂ q i .∂q j j =1
∑
q = qo
La partie du second membre du système formé par les équations de Lagrange linéarisées, issue d’une action mécanique dérivant d’une fonction de force U, se note donc : K (q o , t ).ε + Q(q o , t )
Où K est une matrice carrée symétrique de dimension nxn, où n est la mobilité du système, Q un vecteur de dimension n et ε un vecteur de dimension n, contenant les n variations des paramètres du mouvement autour de la position de départ qo.
98
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
◊
Second membre des équations de Lagrange : Loi visqueuse Lorsque les actions mécaniques extérieure à S appliquées sur S sont décrite par une loi de comportement de type visqueux (efforts aérodynamiques par exemple), le second membre des équations de Lagrange se calcule alors comme suit :
(
) = {F (S → S )}×V q (S / Rg )
Q qi S → S / Rg
i
(
Soit aussi : Q qi S → S / R g
r
) = R{F (S → S )}⋅
∂ r V ( P / Rg ) ∂q& i
Si l’on admet que les efforts découlent d’une loi visqueuse telle que, au point P : r
{ (
r
)} = −C v .V ( P / Rg )
R F S → S
d OP dt R
r
V ( P / Rg ) =
Et comme :
Alors
(
Qqi S → S / R g
et
) = −C v . d OP
dt Rg
⋅
∂ r ∂OP V ( P / Rg ) = ∂q& i ∂q i
∂OP ∂qi
Soit encore :
(
n
) = −C v .∑ ∂OP ⋅ ∂OP
Qqi S → S / R g
∂q j
j =1
dq j
∂qi dt
− C v
∂OP ∂OP ⋅ ∂t ∂qi
La partie du second membre, issue d’effets visqueux, du système formé par les équations de Lagrange linéarisées, peut donc se noter C (q o , t ).
d ε dt
+F (q o , t )
Où C est une matrice carrée symétrique définie positive de dimension nxn et F un vecteur de dimension n, où n est la mobilité du système et ε un vecteur de dimension n, contenant les n variations des paramètres du mouvement autour d’une position qo. On peut remarquer que dans ce cas particulier la matrice C est proportionnelle à la matrice de M.
◊
Système d’équations du mouvement après linéarisation : Dans le cas général, le système formé par les n équations de Lagrange s’écrit donc après linéarisation : M (t , q o ).
d 2 ε 2
dt
+ C (t , q o ).
dw dt
ε + K (t , q o ).ε = F (t , q o )
Où ε est un vecteur de dimension égale à la mobilité du système m, dont les composantes sont les m variations infinitésimales autour d’une position qo repérées par m paramètres du mouvement qi supposés connus à l’instant t . M est appelée la matrice de masse, C la matrice d’amortissement et K la matrice de raideur. La matrice de masse est toujours symétrique, définie positive, de par sa définition.
3.1.2
Exemple : L’étouffeur de vibrations Afin d’illustrer cette procédure de linéarisation, on revient au problème de l’étouffeur de vibrations déjà exposé au chapitre précédent. Le système est constitué d’un moteur rotatif, dont la partie fixe est montée sur un support machine S1, de masse M1 et de centre d’inertie G1, lui-même attaché au bâti S0 par une liaison glissière élastique (paramètre du mouvement X). Par ailleurs, la partie tournante du moteur présente un balourd S3, de masse M3, dont le centre d’inertie G3 est situé à une distance R de G1,
99
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
et dont la position est repérée par l’angle φ. Si le moteur tourne à vitesse constante, le mouvement de S2/S1 est connu : d φ dt
= ω = constante .
On ajoute au système un étouffeur de vibration S2, de masse M2, dont le centre d’inertie G2 est situé à une distance l de G1 et dont la position par rapport à S1 est repérée par l’angle θ.
On va utiliser la méthode de Lagrange pour écrire les équations du mouvement. Le principe des puissances virtuelles s’écrit en utilisant les torseurs de Lagrange, où Σ désigne (S0 U S1 U S2 U S3) : P qi (S / Rg ) = Qqi
(Σ → Σ / Rg ) + Qq (int Σ) i
d ∂ ∂ .T (S / Rg ) et Qqi (S → S / Rg ) = F (S → S )×V qi (S / Rg ) − dt ∂q& i ∂qi
Avec P qi (S / Rg ) =
◊
Calcul des Qqi : Les liaisons pivots entre S3 et S1 et entre S2 et S1 ainsi que la liaison glissière entre S1 et S0 sont supposées parfaites, elles ne développent donc pas de puissances virtuelles. Seuls le poids et l’action du ressort développent de la puissance. Dans les deux cas, on dispose d’une fonction de force : 1 Ressort : U (So ← r → S 1) = − K ( X − X o ) 2 2 Soit : Q X (So ← r → S 1) = − K ( X − X o ) , Qθ (So ← r → S 1) = Qφ (So ← r → S 1) = 0 Pesanteur : L’axe (O,G1) reste confondu avec x, donc : r
r
g r ∂OG1 Q X pes → S 1 / Rg = M 1. g ⋅ = 0, ∂ X Par ailleurs : r
g r ∂OG 2 = − M 2. g .l .Sinθ Qθ pes → S 2 / Rg = M 2. g ⋅ θ ∂ r
r
g g Qφ pes → S 1 / Rg = Qθ pes → S 1 / Rg = 0
r
g g Q X pes → S 2 / Rg = Qφ pes → S 2 / Rg = 0 100
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
◊
Calcul des P qi : L’énergie cinétique du système se calcule comme suit : T (Σ / Rg ) = T (S 1 / Rg ) + T (S 2 / Rg ) + T (S 3 / Rg )
Avec :
dX T (S 1 / Rg ) = M 1. 2 dt 1
2
Et :
dX 1 T (S 2 / Rg ) = M 2. dt 2
2
2 dX d θ d θ + l . + 2 lCosθ dt dt dt
Enfin :
dX 2 dX 2 [ ] + + T (S 3 / Rg ) = M 3. R . ω 2 ω R . . Cos ω t 2 dt dt 1
d ∂ ∂ .T (S / Rg ) − dt ∂q& i ∂qi
Par conséquent, comme : P qi (S / Rg ) =
d 2θ d 2 X + P θ (Σ / Rg ) = M 2.l l . Cosθ dt 2 dt 2 Et P X (Σ / Rg ) = M 1.
◊
2 d 2θ d θ − M 3. R.ω 2 . sin ω t + M 2.l Cosθ − Sin θ 2 2 dt dt dt
d 2 X
Equations de Lagrange : On obtient donc deux équations de Lagrange :
( ) → M 1.
L X
2 d 2θ d θ − M 3. R.ω 2 . sin ω t = − K .( X − X ) + M 2.l Cosθ − Sin θ o 2 2 dt dt dt
d 2 X
d 2θ d 2 X + Cosθ = − M 2. g .l .Sinθ 2 dt 2 dt
( ) → M 2.l l .
L θ
◊
Equations linéarisées Pour pouvoir résoudre ce système, on se place en un point (Xi,θi) et on fait un développement limité au premier ordre autour de ce point des équations du mouvement obtenues. On rappelle que les développements limités de Cosθ et Sinθ autour d’un point θi s’écrivent : Cos(θ i
+ ε ) = Cosθ i − ε .Sinθ i + O(ε 2
Sin(θ i
+ ε ) = Sinθ i + ε .Cosθ i + O ε 2
Pour linéariser, on suppose que les paramètres X et θ sont peu différents de la position de référence (Xi,θi), et par ailleurs, on fait également l’hypothèse que les dérivées temporelle des paramètres
101
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
du mouvement sont petites devant celles au point de référence. Après linéarisation, les deux équations du mouvement s’écrivent : En posant : X = X i + χ
( ) → M 1.
L X
d 2 χ 2
dt
θ = θ i
et
+ M 2.l .Cosθ i
( ) → M 2.l .Cosθ i
L θ
d 2 χ dt 2
d 2 ε 2
dt 2
+ M 2.l .
+ ε
+ K . χ = − K .( X i − X o ) + M 3. R.ω 2 . sin ω t
d 2 ε dt 2
+ M 2.l . g .ε .Cosθ i = − M 2. g .lSinθ i
Le système obtenu peut alors se mettre sous la forme matricielle suivante : M .
d 2 w dt 2
+ K .w +F (t ) = 0
Avec :
χ w = , ε
M 2.l .Cosθ i M 1 , M = 2 2 . . 2 . M l Cos θ M l i
et
F
2
. .ω (t ) = − K ( X i − X o ) + M 3 R
− M 2. g .l .Sinθ i
3.2
0 K K = 0 2 . . . M g l Cos θ i
. sin ω t
Rappels : Résolution de systèmes linéaires d’équations différentielles Dans le cas général, nous aurons donc à résoudre des systèmes d’équations différentielles linéaires de la forme suivante : M (t , q o ).
d 2 w 2
dt
+ C (t , q o ).
dw dt
+ K (t , q o ).w = F (t , q o )
Où w est un vecteur de dimension égale à la mobilité du système m, dont les composantes sont les m variations infinitésimales autour d’une position qo repérées par m paramètres du mouvement qi supposés connus à l’instant t . La matrice de masse M est toujours symétrique, définie positive. Nous distinguerons deux cas. Si les matrices M, C et K sont à coefficients constants le système pourra être résolu analytiquement. Sinon, le problème sera résolu numériquement. La résolution de systèmes d’équations différentielles linéaires est un problème commun à de nombreuses disciplines, cette partie vise donc essentiellement à donner quelques rappels rapides sur les différentes techniques à employer selon les cas rencontrés. n.b. En mécanique du solide indéformable, dans le cas où l’on applique une procédure de résolution numérique, il vaut mieux paramétrer les rotations à l’aide d es coordonnées des quaternions plutôt qu’à l’aide d’angles de rotations.
3.2.1
Résolution numérique. On considère un système d’équations différentielles linéaires de dimension n. Dans ce cas, nous supposons connues à un instant donné t , toutes les caractéristiques du mouvement, et nous cherchons à les déterminer à l’instant t+h. h étant le pas de temps choisi. Divers algorithmes de résolution numériques existent, qui ne sont pas l’objet de ce cours. On donne juste ici un aperçu des méthodes les plus courantes. Les méthodes les plus employées sont les méthodes du point milieu et de Runge-Kutta. Pour appliquer ces méthodes, il faut transformer le système d’équations différentielles linéaires du second ordre de dimension n en un système d’équations différentielles linéaires du premier ordre de dimension 2n. On
102
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
pose un vecteur y dont les composantes à l’instant t sont, les n paramètres du mouvement et les n vitesses de variation de ces paramètres :
w dw dt
y =
Alors, M étant inversible, le système : M .
d 2 w dt 2
+ C
1 0 = −1 −1 y dt − M K − M C
dy
dw dt
+ K .w = F (t ) s’écrit maintenant:
0 + −1 M F (t )
Ayant fixé un pas de temps h, on construit alors une suite :
t i +1 = t i + h y i +1 = y i + m.h Où m est une pente moyenne sur l’intervalle (ti, ti+1), si l’on pose : m1
=
dy dt (t , x ) i
i
m2
=
dy
m3
dt (t + h , x + h m ) i
2
i
1
2
=
dy dt (t + h , x + h m ) i i 2 2
m4
2
=
dy dt (t + h, x + hm ) i i 3
Illustration graphique de la méthode du point milieu Pour la méthode du point milieu, on choisit m=m2 comme pente moyenne sur l’intervalle (ti, t i+1). Cette méthode est illustrée graphiquement ci-dessus. Pour la méthode de Runge Kutta (plus précise), la pente moyenne sur l’intervalle (ti, ti+1) se calcule comme suit: m=
3.2.2
1 6
(m1 + 2.m2 + 2.m3 + m 4 )
Résolution analytique Un système d’équations différentielles du second ordre à coefficients constants, et dont le second membre est périodique et ne dépend que du temps, peut, dans certains cas, être résolu de manière analytique : M .
d 2 w 2
dt
+ C
dw dt
+ K .w = F (t )
103
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
M, C et K sont des matrices carrées de dimension mxm où m est la mobilité du système, w un vecteur inconnu de dimension m et F (t) un vecteur connu de dimension m. La matrice M est inversible.
◊
Cas N°1 : Le second membre F (t) est nul. Le système s’écrit alors : M .
d 2 w dt 2
+ C
dw dt
+ K .w = 0
On cherche des solutions de la forme : w(t ) = r .e ω .t , où r est un vecteur de dimension m, et où ω est un nombre complexe. Etant donné que e0=1, le vecteur r contient les valeurs initiales des paramètres du mouvement. Il reste à déterminer ω. En remplaçant dans l’équation différentielle ci-dessus il vient : M .
+ K .w = 0 ( .ω 2 + C .ω + K ).r .eω .t = 0 dt ⇒ M dt 2 w(t ) = r .e ω .t
d 2 w
+ C
dw
En dehors de la solution triviale r =0, on a des solutions à cette équation à condition que :
∆ = D(ω ) = Det M ( .ω 2 + C .ω + K = 0 Cette équation, dite équation caractéristique, est un polynôme de degré 2.m, où m est la mobilité du système. Elle admet donc 2.m racines (distinctes ou non) qui sont soit réelles, soit complexes conjuguées. On obtient l’ensemble des fonctions à valeurs réelles solution du système en prenant les parties réelles des solutions complexes conjuguées. Trois cas particuliers se présentent :
-
Solution imaginaire pure : ω = i b : b réel. La solution est alors une solution périodique de type Cos(b t ) . Cette solution est compatible avec l’hypothèse de petits mouvements qui est sous-jacente à l’opération de linéarisation.
-
Solution à partie réelle négative : ω =a +i b, avec a et b réels et a<0. La solution est alors une fonction composée de type : e a t .Cos(b t ) avec a < 0. Cette fonction est décroissante cycliquement et reste également compatible avec l’hypothèse des petits mouvements.
104
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
-
Solution à partie réelle positive : ω =a +i b, avec a et b réels, a > 0. La solution est alors une fonction composée de type : e a t .Cos(b t ) avec a > 0. Cette fonction est croît cycliquement de manière exponentielle et conduit donc à des mouvements qui deviennent rapidement incompatibles avec l’hypothèse des petits mouvements, qui a permit d’établir le système d’équations différentielles linéaires étudié.
L’étude du signe de la partie réelle des solutions de l’équation caractéristique est donc primordiale pour juger de la stabilité de la validité de la solution obtenue. Le critère de Routh qui est donné plus loin permet de discuter du signe des parties réelles des solutions de l’équation caractéristique sans résoudre cette équation.
◊
Cas N°2 : Le second membre F (t) est nul et C est nulle. Le système s’écrit alors : M .
d 2 w 2
+ K .w = 0
dt
Comme précédemment le système admet des solutions de la forme w(t ) = r .e ω .t , où ω est cette fois ci imaginaire pur , ω =i.b, avec b réel, et où r est un vecteur de dimension m contenant les valeurs
(
initiales des paramètres du mouvement. Les solutions l'équation caractéristique Det M .b 2 − K
= 0 sont
appelées les pulsations propres du système. On obtient m solutions imaginaires pures conjuguées. Comme les matrices M et K sont symétriques définie positives, et compte tenu des solutions trouvées, le système linéaire peut être réécrit de la façon suivante : M −1 K .w(t ) = b 2 .w(t ) Les solutions cherchées sont donc des vecteurs propres de la matrice D = −1 K , c'est-à-dire que les solutions w sont telles que : D.w=λ.w, avec ici λ=b2 réel et positif. Les vecteurs propres sont aussi appelés « modes propres ». On peut donc définir une base propre à partir des vecteurs propres de la matrice D, et une matrice de rotation Q constituée de ces vecteurs propres : Q = (w1 ,..., wm ) . Si x est un vecteur exprimé dans la base propre, x est solution si : λ x(t ) = Q T . D.Q x(t )
où
Q T . D.Q est diagonale
105
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
Cette propriété implique que, écrites dans la base propre, les équations du mouvement sont complètement découplées. Ce qui se dit aussi plus succinctement : « les modes propres sont découplés ». Toutes les combinaisons linéaires de modes propres sont également solutions. Une combinaison linéaire de modes propres se définit à l’aide de ses « coordonnées modales » dans la base propre.
◊
Cas N°3 : Le second membre F (t) est non-nul et C est nulle. Le système s’écrit alors : M .
d 2 w dt 2
+ K .w = F (t )
Un moyen simple de résoudre un tel système est de le projeter dans la base des modes propres du système sans second membre, définit par la matrice Q (voir Cas N°2). Alors on obtient un système de m équations du mouvement découplées, de la forme suivante, où ε est un vecteur de la base modale : Q T M . .Q
d 2 ε 2
dt
+ Q T . K .Q.ε = Q T .F (t )
& & On note alors la ième équation à résoudre : mi .ε i
+ k i .ε i = f i (t )
Chaque équation se résout indépendamment en cherchant des solutions du type : ε i
= ai .Cos(ω i .t + ϕ i ) + ε ip (t )
avec
=
ω i
k i mi
Où εi p est une solution particulière et ωi la pulsation propre associée au i
ème
mode propre.
Si le second membre f i(t) est périodique de pulsation ω, on peut le décomposer en séries de Fourier, et le système s’écrit alors :
(
&&i mi . ε
+ ω i2 .ε i
+∞
) = ∑ C n .e i.n.ω .t i
n = −∞
Dans le cas le plus simple, où la fonction est sinusoïdale :
(
&&i mi . ε
+ ω i2 .ε i ) = C f .Cos (ω f .t + ϕ f )
La solution de cette équation est alors la suivante, ai et ϕ i se déterminant à l’aide des valeurs initiales : ε i (i ) = a i .Cos(ω i .t + ϕ i ) +
C f .Cos ω f .t + ϕ f
(
mi ω i2
− ω f 2
)
Cette solution présente une singularité pour ωi=ωf , qu’on appelle une résonance
◊
Cas N°4 : Le second membre F (t) est nul et C est une combinaison linéaire de M et K. Le système s’écrit alors, avec a et b deux scalaires réels : M .
d 2 w dt 2
. + b K . ) + (a M
dw dt
+ K .w = 0
Ce système s’écrit aussi : d 2 w 2
dt
+ (aI + b M . −1 K )
dw dt
+ M −1 K .w = 0
106
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
Les vecteurs propres de D = −1 K sont encore vecteurs propres de combinaisons linéaires quelconques de D et de l’identité, donc le système est toujours découplé une fois projeté dans la base modale du système : M .
d 2 w dt 2
+ K .w = 0
Après projection le système s’écrit comme suit : T
Q M . .Q
d 2 ε 2
dt
. + b K . )Q + Q T (a M
d ε dt
+ Q T K .Q.ε = 0
Il reste alors à résoudre m équations différentielles découplées, du type : && ε i
+
ci mi
.ε &i
+ ω i2 .ε i = 0
La solution de chacune de ces équations s’écrit comme suit : ε i (t ) = e
− d i .t
Re a e
−t d i2 −ω i2
+ b e t
d i −ω i 2
2
avec
d i
=
ci 2.mi
Ceci correspond en général à une solution sinusoïdale amortie.
◊
Critère de Routh Le critère de Routh permet de discuter le signe de la partie réelle des solutions d’une équation algébrique à coefficients réels, sans résoudre cette équation. Supposons une équation, telle que ao>0 : a o . y n
+ a1 . y n −1 + a 2 . y n−2 + ... + a n = 0
On construit alors le tableau suivant, par « déterminant » successifs, jusqu’à n’avoir que des zéros en bas de chacune des colonnes. Equation de degré pair : construction de la table : a0
a2
a4
a6
…
a1
a3
a5
a7
…
a1.a2-a0.a3
a1.a4-a0.a5
a1.a6-a0.a7
…
(a1.a2-a0.a3).a3-( a1.a4-a0.a5).a1
(a1.a2-a0.a3).a5-( a1.a6-a0.a7).a1
…
…
…
…
Les deux première lignes du tableau contiennent les coefficients des termes d’exposant pair ou impair. Le degré de l’équation donne le type des coefficients de la première ligne (pair ou impair). Si le degré est impair, la première ligne contient les coefficients des termes d’exposants impairs.. Pour que les racines de l’équation aient des parties réelles négatives, il faut et il suffit que les éléments de la première colonne soient tous positifs Si les éléments de la première colonne changent de signe, le nombre de changement de signe donne le nombre de racines dont la partie réelle est positive. On en déduit, que les parties réelles sont toutes négatives si : • • •
Pour n=1 : a0, a1>0 Pour n=2, a0, a1, a2>0 Pour n=3, a0, a1, a2, a3>0 et a0.a3-a1.a2>0
107
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL •
◊
Pour n=4, a0, a1, a2, a3, a4 >0 et a3.(a1a2-a0a3)-a1.(a1.a4) >0
Rappel : Calcul du déterminant d’une matrice Pour une matrice 2 × 2, on montre que la matrice inverse est donnée par :
Le nombre ad - bc est appelé déterminant de la matrice A, noté :
La matrice inverse A-1 n'existe donc que si det A est différent de zéro. La matrice A est singulière si det A = 0, régulière dans le cas contraire. Ce résultat se généralise à une matrice de dimension quelconque. Le déterminant peut se calculer de manière récursive. Par exemple, pour n = 3, on a, en développant par rapport à la première ligne :
Dans ce développement, chaque déterminant d'ordre 2 est appelé mineur du terme qui le précède. Par exemple, le mineur de a est :
On peut développer le déterminant par rapport à n'importe quelle ligne ou colonne. Pour chaque élément aij de la ligne ou colonne choisie : L mineur est le déterminant de la sous-matrice obtenue en supprimant la ligne i et la colonne j Le signe du produit est donné par le tableau ci-dessous : +
-
+
-
+
-
+
-
+
108
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
◊
4
4.1
Rappel : Matrice 2x2 valeurs propres, vecteurs propres
EQUILIBRE ET STABILITE
Introduction L’objectif de cette partie est de présenter diverses méthodes permettant de définir et d’étudier la stabilité d’une position d’équilibre. Deux méthodes peuvent être envisagées. La méthode directe, qui s’appuie sur le théorème de Lejeune-Dirichlet, permet de définir les positions d’équilibre indépendamment des équations du mouvement. Cette méthode est bien adaptée aux système dont les forces dérivent d’un potentiel indépendant du temps. La méthode générale, consiste à écrire les équations du mouvement linéarisées au voisinage des positions d’équilibre ce qui permet de discuter de leur stabilité.
4.2
Systèmes conservatifs Si les puissances des efforts intérieurs et extérieurs sont nulles ou bien se calculent par dérivation d’une énergie potentielle qui ne dépend pas du temps, dont on notera la somme E p(Σ/Rg), alors l’intégrale première de l’énergie cinétique traduit directement la conservation de l’énergie du système : T (Σ / Rg ) + E p (Σ / Rg ) = Constante
Le système est alors dit « conservatif ». L’énergie potentielle se transforme en énergie cinétique et réciproquement au cours du mouvement. Les équations de Lagrange pour un système conservatif, donnent :
109
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
(
)
P qi (S / Rg ) = Qqi S → S / Rg
(*) ∂U S → S / Rg d ∂ ∂ (*) − . T (S / Rg ) et Qqi S → S / Rg = Avec P qi (S / Rg ) = ∂q i dt ∂q& i ∂qi Soit : d
∂ ∂ (T + U ) = 0 T (S / Rg ) − dt ∂q& i ∂qi Et comme, pour un système conservatif, la fonction de force U ne dépend que des paramètres de position qi :
d ∂ ∂ & − (T + U ) = 0 dt ∂qi ∂qi On définit alors la fonction L, dite Lagrangien du système, telle que : L = T + U Alors les équations de Lagrange, pour un système conservatif s’écrivent aussi :
d ∂ ∂ & − L = 0 dt ∂q i ∂qi
◊
Equilibre Les positions d’équilibre sont telles que les paramètres du mouvement restent constants. Dans cette configuration, l’énergie cinétique est nulle et les dérivées par rapport au temps sont nulles. On déduit alors des équations de Lagrange, pour toute position d’équilibre :
d ∂ ∂ ∂ L ∀i ∈ (1, n ) → − L = 0 ⇒ =0 ∂qi q = q dt ∂q& i ∂qi e Et comme, à l’équilibre, l’énergie cinétique est nulle, ceci équivaut à écrire :
∀i ∈ (1, n ) →
∂ E p ∂U =0⇔ =0 ∂qi q =q ∂qi q =q e e
Les positions d’équilibre sont donc obtenues pour les extremums (maximums, minimums ou points d’inflexion locaux) de la fonction de force U ou de l’énergie potentielle Ep.
◊
Stabilité : Définition Une position d’équilibre d’un système matériel est dite stable si le mouvement induit par une perturbation (conditions initiales en vitesse et en position très petites) reste au voisinage de cette position.
◊
Stabilité : Théorème de Lejeune Dirichlet D’après l’intégrale première du théorème de l’énergie cinétique, pour un système conservatif Σ, nous pouvons écrire : T (Σ / Rg )to + E p (Σ / Rg )to
= T (Σ / Rg )t + E p (Σ / Rg )t = Constante
Supposons que l’énergie potentielle présente un minimum local strict dans la position d’équilibre étudiée. Comme l’énergie potentielle est définie à une constante près, on peut choisir une expression de cette énergie potentielle telle que le minimum local de l’énergie potentielle soit nul dans la position d’équilibre étudiée.
110
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
Si l’on choisit une petite perturbation, telle que à to : T (Σ / Rg )to
et E p (Σ / Rg )to
≤ α
T (Σ / Rg )to + E p (Σ / Rg )to
≤ α soit
≤ 2α
α étant très petit. Comme à tout instant : 0 ≤ T (Σ / Rg )t et 0 ≤ E p (Σ / Rg )t On en déduit que par la suite, à tout instant : 0 ≤ T (Σ / Rg )t ≤ 2α et 0 ≤ E p (Σ / Rg )t ≤ 2α •
Comme l’énergie potentielle ne dépend que de la position, 0 ≤ E p (Σ / Rg )t ≤ 2α , signifie que tous les q i
− q ie restent très petits.
•
Comme dans un système conservatif, l’énergie cinétique est une forme quadratique des vitesses de variation des paramètres du mouvement, (les liaisons sont indépendantes du temps), dq i alors : 0 ≤ T (Σ / Rg )t ≤ 2α signifie que les restent également très petits. dt
•
La position d’équilibre est donc stable.
Théorème de Lejeune Dirichlet : Etant donné un système conservatif, tout maximum local strict de la fonction de force (tout minimum local strict de l’énergie potentielle), définit une position d’équilibre stable. Pour un système conservatif, la recherche des positions d’équilibre et l’étude de leur stabilité peut donc se faire directement à l’aide de l’énergie potentielle du système, sans qu’il soit besoin d’écrire complètement les équations du mouvement.
◊
Extension aux systèmes visqueux On considère le cas particulier ou le système développe une résistance passive opposée à la vitesse (force visqueuse, traînée aérodynamique). Les équations de Lagrange pour un tel système, donnent toujours :
(
)
P qi (S / Rg ) = Qqi S → S / Rg
d ∂ ∂ ∂U r visq ∂ r ∂ & − ∂ .T (S / Rg ) et: Qqi (S → S / R g ) = ∂ + R F S → S ⋅ ∂ & V ( P / Rg ) . qi dt qi qi qi
Avec P qi (S / Rg ) =
En tenant compte de la nature visqueuse de la résistance passive :
(
Qqi S → S / R g
) = ∂∂U − C v V ( P / Rg ) ⋅ ∂∂& r
qi
qi
r
V ( P / Rg )
111
Préparation à l’agrégation l’agrégation de mécanique mécanique : Polycopié Polycopié de Mécanique Mécanique des Solides Indéformab Indéformables les APPLICATIONS, TECHNIQUES TECHNIQUES DE CALCUL CALCUL
En définissant alors la fonction L, dite Lagrangien du système, telle que : L = T + U ; les équations de Lagrange, pour un système visqueux s’écrivent alorsi : r d ∂ ∂ ∂ r − = − ⋅ ( ) L C V P / Rg V ( P / Rg ) & v ∂q& i dt ∂qi ∂qi
Dans toute position d’équilibre la vitesse est nulle et le second membre de cette équation est donc nul également, ce qui nous ramène au cas précédent pour la détermination des positions d’équilibres. Pour un système visqueux, les positions d’équilibre correspondent encore aux extremums de la fonction de force du système. Pour l’étude de la stabilité, on peut encore appliquer le théorème de Lejeune Dirichlet. En effet, d’après le principe de conservation de l’énergie mécanique : dT dt
dU
=
dt
+ P visq
Entre deux instants t et to on peut toujours écrire : t
dT
∫ dt
t
dt =
to
dU
∫ dt
t
dt +
to
∫
P visq
to
dt
dt
soit
T t − T to
= U t − U to + W visq
Ce qui s’écrit aussi : T (Σ / Rg )to + E p (Σ / Rg )to
+ W visq = T (Σ / Rg )t + E p (Σ / Rg )t
avec dans tous les cas W visq
<0
Si, comme précédemment, on choisit une petite perturbation α, telle que à to : T (Σ / Rg )to
≤ α
et E p (Σ / Rg )to
T (Σ / Rg )to + E p (Σ / Rg )to
≤ α soit
≤ 2α
Comme à tout instant : 0 ≤ T (Σ / Rg )t et 0 ≤ E p (Σ / Rg )t On en déduit que par la suite, à tout instant : 0 ≤ T (Σ / Rg )t ≤ 2α + W visq et 0 ≤ E p (Σ / Rg )t ≤ 2α + W visq Donc Wvisq étant négative, la position d’équilibre est encore stable.
◊
Equilibre et stabilité pour un système à un seul paramètre La position d’équilibre est un aussi un minimum local si :
∂U =0 ∂q q =q e
◊
et
∂ 2U <0 2 ∂q q =q e
Equilibre et stabilité pour un système à n paramètres La position qe est une position d’équilibre si :
∂U = 0 ∀i ∂qi q =qe Cette position d’équilibre est stable si la matrice négatives.
112
∂ 2U ∂qi ∂q j
= aij est à valeurs propres toutes q = qe
Préparation à l’agrégation l’agrégation de mécanique mécanique : Polycopié Polycopié de Mécanique Mécanique des Solides Indéformab Indéformables les APPLICATIONS, TECHNIQUES TECHNIQUES DE CALCUL CALCUL
Il est généralement plus simple de procéder à un développement limité à l’ordre 2 de la fonction de force U au voisinage de la position d’équilibre. Soit : n
n n ∂U ∂ 2U 1 ( (q i − qi e ) + U = U q = qe + q j − q j )(q i − q i e ) e ∂qi q =qe 2 j =1 i =1 ∂qi ∂q j i =1
∑
∑∑
Or la condition d’équilibre impose :
q =qe
∂U = 0 ∀i ∂qi q =qe
Donc le développement limité s’écrit : U = U q = qe
+
n
n
2
∂ U ( q j − q j )(q i − q i e ) ∑ ∑ e ∂q ∂q 2
1
i
j =1 i =1
j q = qe
Pour que la position d’équilibre soit stable, il suffira alors que dans tous les cas : n
n
∂ 2U <0 (q j − q j e )(qi − qi e ) ∂q ∂q i j j =1 i =1 q = qe
∑∑ ◊
Exemple, étouffeur de vibration Si l’on revient au cas de l’étouffeur de vibration déjà vu plus haut. Dans le cas où, la vitesse de rotation ω du balourd est nulle, on peut chercher des positions d’équilibre à l’aide du théorème de Lejeune Dirichlet étant donné que toutes les forces dérivent d’un potentiel. On peut alors écrire pour le système complet : Q X (int ) = − K ( X − Xo) =
∂U ressort ∂ X
Qθ (int ) = − M 2 . g .l .Sinθ =
∂U pesanteur ∂θ
Les positions d’équilibre sont telles que Qqi=0 soit : X = Xo et Sinθ = 0 Il existe donc deux positions d’équilibre, l’une pour laquelle la masse M2 est en haut (θ=π) et l’autre pour laquelle elle est en bas (θ=0)
∂ 2U = − K Alors : a11 = ∂ X 2
a 22
∂ 2U = 2 = − M 2 . g .l .Cosθ ∂θ
a12
∂ 2U = a 21 = =0 ∂θ .∂ X
En appliquant le critère précédent, la position d’équilibre est donc stable si : est une matrice à valeurs propres négatives. Ce qui impose que :
∂ 2U = − K < 0 , ce qui est toujours vrai a11 = ∂ X 2 a 22
∂ 2U = 2 = − M 2 . g .l .Cosθ < 0 ∂θ
ce qui implique que Cos θ soit positif.
Comme on l’attendait, seule la position θ=0 est stable, c'est-à-dire celle pour laquelle la masse M2 est en bas.
113
Préparation à l’agrégation l’agrégation de mécanique mécanique : Polycopié Polycopié de Mécanique Mécanique des Solides Indéformab Indéformables les APPLICATIONS, TECHNIQUES TECHNIQUES DE CALCUL CALCUL
4.3
Cas général, méthode directe ou méthode de Liapounov Dans le cas général, la méthode décrite ci-dessus ne permet pas de déterminer les positions d’équilibre. En effet, en absence de potentiel par exemple, la méthode de Lejeune Dirichlet ne pourra pas être appliquée. On appliquera donc la méthode directe.
◊
Etat de mouvement Si l’on considère un système mécanique à n paramètres qi, l’état du mouvement de ce système est défini par une solution connue du système différentiel formé par les équations du mouvement :
= qi (t ) ∀i = 1, n
qi
-
Equilibre : C’est une solution particulière pour laquelle tous les paramètres sont constants. q& i
-
= q&&i = 0 ∀i = 1, n
Mouvement stationnaire : C’est un problème plus général pour lequel, certains paramètres du mouvement sont constants, tandis que d’autres voient leur vitesse de variation rester constante. Ce type de mouvement est fréquemment rencontré lors de l’étude du fonctionnement de machines en régime permanent. qi (t ) = cst
◊
∀i = 1, j et q& i (t ) = cst ∀i = j, n
Détermination de l’équilibre par la méthode directe Dans le cas général, on écrit les équations du mouvement et on cherche directement les positions telles que les paramètres du mouvement correspondant à l’état du mouvement étudié restent constants. &&i ∀i = 1, n par zéro, dans les équations du mouvement, et on S’il s’agit d’un équilibre on remplace q& i et q cherche une solution qi vérifiant le système d’équation restant.
◊
Stabilité au sens de Liapounov
(
)
~& (t ) deux états du mouvement distinct provenant de conditions initiales Soit (q (t ), q& (t )) et q~(t ), q (q o (t ), q& o (t )) et q~o (t ), q~& o (t ) très légèrement différentes. On dira que le mouvement est stable si :
qio − q~io < λ ∀λ tel que il existe ε tel que ~& q& io − qio < λ
qi (t ) − q~i (t ) < ε >0. q& (t ) − q~& (t ) < ε avec t >0. i i
On peut aussi exprimer la condition de stabilité de manière globale à l’aide d’une fonction écart qui doit rester petite : r 2
=
∑ (q (t ) − q~ (t )) + ∑ (q (t ) − q~ (t )) i
i
2
i
◊
&i
&
2
i
i
Théorème de Liapounov Dans le cas général, lorsqu’on a un équilibre au sens strict du terme, on étudie sa stabilité en linéarisant les équations su mouvement autour de la position d’équilibre étudiée. On aboutit alors à un système d’équations différentielles linéaires à coefficients constants, qui peut être résolu analytiquement : M .
d 2 w 2
dt
+ C
dw dt
+ K .w = F (t )
114
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
L’oscillateur linéaire L associé à ce système est alors le suivant : M .
d 2 w dt 2
+ C
dw dt
+ K .w = 0
On cherche des solutions de la forme : w(t ) = r .e ω .t où les pulsations propres
ω sont solutions de
l’équation caractéristiques du système :
∆ = D(ω ) = Det M .ω 2 + C .ω + K = 0 Théorème de Liapounov : Soit L l’oscillateur linéaire associé au système mécanique étudié et à la position d’équilibre So. Si toutes les racines de l’équation caractéristique de L ont leur partie réelle négative alors So est stable. Si l’une au moins a sa partie réelle positive alors So est instable. Si l’une au moins des solutions à sa partie réelle nulle, on ne peut pas conclure.
◊
Stabilité asymptotique En considérant les mêmes états du mouvement que dans le paragraphe précédent, on dira que le système présente une stabilité asymptotique si : t →∞ q (t ) − q~ (t ) → 0 i i → ∞ t ~ & q& i (t ) − q i (t ) → 0
◊
Stabilité orbitale Si l’on considère deux points d’origine très voisins dans un système de solides, le système reste stable si : MM * (t o )
< λ ⇒ MM * (t ) < ε ∀t
Si l’on considère par exemple deux satellites tournant sur des orbites très proches mais avec des périodes différentes. Si les deux satellites sont très voisins à l’instant to, la différence des périodes peut entraîner des écarts importants lorsque t augmente.
5
◊
VIBRATIONS AUTOUR D’UNE POSITION D’EQUILIBRE STABLE
Introduction Pour un système donné, si des positions d’équilibre stables ont été identifiées, on peut s’intéresser aux petits mouvements autour de ces positions d’équilibre. L’étude des petits mouvements se fait à l’aide des équations du mouvement linéarisées, qui ont la forme suivante dans le cas général : M (t , q o ).
d 2 w 2
dt
+ C (t , q o ).
dw dt
+ K (t , q o ).w = F (t , q o )
Lorsque les mouvements du systèmes sont de petits mouvement autour d’une position d’équilibre stable qo=qe, alors les matrices M, C et K ne dépendent plus du temps. M (q e ).
d 2 w 2
dt
+ C (q e ).
dw dt
+ K (q e ).w = F (t , q e )
115
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
Seul le second membre est susceptible de dépendre du temps. On se trouve alors dans le cas de systèmes linéaire d’équations différentielles à coefficients constant, qui peuvent être résolus de manière analytiques.
◊
Vibrations libres Lorsque le second membre de cette équation est nul, que la matrice C est nulle et que K est symétrique définie positive. Le système est un oscillateur libre : M (q e ).
d 2 w 2
dt
+ K (q e ).w = 0
Les solutions sont des combinaisons linéaires de vecteurs w(t) sinusoïdaux, appelées modes propres, qui sont les vecteurs propres de la matrice M-1K, multipliés par Cos(ωt), où ω2 est la valeur propre correspondante.
◊
Vibrations forcées Lorsque le second membre est non nul et périodique de pulsation ωf, que la matrice C est nulle et que K est symétrique définie positive, on parle alors de vibrations forcées. M (q e ).
d 2 w dt 2
+ K (q e ).w = F (t )
Les solutions de ce système d’équations, sont des combinaisons linéaires de modes propres plus une solution particulière, qui présente des singularités en ωf=ωi où les ωi sont les i pulsations propres du système.
◊
Vibrations amorties Lorsque le second membre est nul, que la matrice C est une combinaison linéaire quelconque de M et K, on parle alors de vibrations amorties. M (q e ).
d 2 w 2
dt
+ C (q e ).
dw dt
+ K (q e ).w = 0
Projetée dans la base définie parles vecteurs propres de la matrice M-1K, le système est alors complètement découplé. Et les solutions son des combinaisons linéaires de solution du type suivant : 2 2 c − d .t −t d i2 −ω i2 + b e t d i −ω i avec d i = i ε i (t ) = e i Rea e 2.m i Ceci correspond en général à une solution sinusoïdale amortie.
6
6.1
MECANIQUE DES CHOCS-PERCUSSIONS
Introduction Lorsque deux solides entrent en contact au cours de leur mouvement, il se produit ce que l’on appelle un « choc ». On s’intéresse à l’état du mouvement juste avant (instant t 1) et juste après (instant t 2) le choc, t 1 et t2 étant très proches, de l’ordre de un centième à un millième de secondes. Pendant cet intervalle de temps, la position d’un point M a très peu changé, tandis que sa vitesse a subit une forte variation. Si l’on
116
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
confond les instants t 1 et t 2, on parlera alors d’une discontinuité du vecteur vitesse. De ce fait les techniques de calcul exposées plus haut ne pourront pas être appliquées. L’analyse précise du choc, c'est-à-dire de ce qui se passe entre les instants t 1 et t 2, est un problème très complexe. Lorsque M arrive au contact avec l’obstacle, il se produit des déformations locales qui dépendent de la nature des matériaux en présence, de la vitesse d’impact, de la nature des surfaces… Elle nécessite une étude complète en mécanique des milieux continus, en viscoplasticité et en dynamique et avec établissement et perte de contact, ce qui reste un problème difficile à traiter. En mécanique des solides on cherche une modélisation simplifiée permettant, si l’on connaît l’état du mouvement avant le choc de déterminer l’état du mouvement après le choc.
6.2
Cas d’un point matériel En l’absence de choc, le principe fondamental de la dynamique s’écrit dans un repère Galiléen Rg, pour un point matériel M de masse m : r
r
m.Γ( M / Rg ) = F ext
En intégrant entre deux instants t1 et t2 , il vient : r
t =t 2
r
m.V (t 2 , M / Rg ) − m.V (t 1 , M / Rg ) =
∫
r
F ext .dt
t =t 1
Le second membre de cette équation est appelé l’impulsion intégrée du vecteur F ext pendant l’intervalle de temps (t 2-t 1). Lorsque la limite de ce second membre, lorsque l’intervalle de temps tend vers zéro, est un vecteur fini non nul, on parle alors de percussion. Ceci correspond physiquement à l’application d’une force très grande pendant un temps très court, du type « Dirac ». t =t 2
r
P =
lim (t 2−t 1)→0
∫
r
F ext .dt
t =t 1
On retrouve alors la discontinuité du vecteur vitesse, à l’instant t, que l’on note de la manière suivante : r r m.V t + , M / Rg − m.V t − , M / Rg
(
(
r
r
= [[m.V (t , M / Rg )]]t = P
Par définition P est appelée le vecteur-percussion à l’instant du choc. Sa dimension physique est M.L.T-1. Si le vecteur percussion est donné, il détermine complètement la discontinuité du vecteur vitesse lors du choc.
6.3
Cas d’un solide ou d’un système de solides Pour un système de solides Σ, le principe fondamental de la dynamique postule que, il existe au moins un repère Rg appelé repère Galiléen, le torseur dynamique de Σ dans son mouvement par rapport à Rg est égal au torseur des actions mécaniques extérieures à Σ. Ce qui se note :
{D (Σ / Rg )} = {F (Σ → Σ )} Soit m la masse et G le centre d’inertie du système matériel référentiel galiléen Rg. On a alors :
r
r
m.Γ(t , G / Rg ) {D (Σ / Rg )}= r δ G (t , Σ / Rg ) G
Σ en mouvement par rapport au
et
(Σ → Σ ) r {F (Σ → Σ )}= R M G (Σ → Σ ) G 117
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
En intégrant comme précédemment entre deux instants t1 et t2 on obtient :
t 2 r Σ → Σ R dt ( ) r r m.(V (t 2, G / Rg ) − V (t 1, G / Rg )) t 1 r = t 2 r σ G (t 2, Σ / Rg ) − σ G (t 1, Σ / Rg ) r G M G (Σ → Σ )dt G t 1
∫
∫
En supposant que les limites des intégrales des efforts et moments extérieurs ont des limites finies non nulles lorsque l’intervalle de temps tend vers zéro, on obtient l’expression de la discontinuité du vecteur quantité de mouvement et du moment cinétique. t 2 r r = Σ → Σ P lim R dt ( ) r (t 2 −t 1)→0 m.[[V (t , G / Rg )]]t t 1 = r t 2 [[σ G (t , Σ / Rg )]]t r r G M G (Σ → Σ )dt K G = (t 2−lim t 1)→0 t 1 G
∫
∫
r
m.[[V (t , G / Rg )]]t r est la discontinuité du torseur cinétique du système matériel Σ à [ [ ] ] t , / Rg ( ) Σ σ G t G
•
l’instant du choc t . r
P r est le torseur des percussions extérieures appliquées au système pendant le choc. K G G
•
On peut donc énoncer le principe fondamental des chocs : Pour tout système matériel en mouvement par rapport à un repère galiléen Rg, pendant un choc, la discontinuité du torseur cinétique du système matériel à l’instant du choc, est égal au torseur des percussions extérieures appliquées au système. Ce qui se note :
[[C (Σ / Rg )]]t = P (Σ → Σ )
◊
Remarque 1 : A partir du principe fondamental des chocs, on peut, comme pour le principe fondamental de la dynamique, démontrer le théorème des percussions réciproques.
(
P Σ1
◊
→ Σ 2 ) = −P (Σ 2 → Σ1 )
Remarque 2 : Toutes les forces qui restent bornées pendant le choc, telles que par exemple, les forces de pesanteur, correspondent à un torseur de percussion nul.
◊
Remarque 3 : Le principe fondamental des chocs est défini dans un repère galiléen, cependant lorsqu’on écrit le principe fondamental de la dynamique dans un repère R en mouvement par rapport à Rg, si les forces d’inertie d’entraînement et de Coriolis restent finies au cours du choc, elles correspondent alors à des percussions nulles. On pourra dans ce cas appliquer de la même façon le principe fondamental des chocs dans le repère non Galiléen R.
6.3.2
Percussion de liaison Considérons une liaison entre deux solides S1 et S2. Si le système subit un choc, il faut introduire à l’instant du choc une percussion de liaison :
118
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
( → S 2 ) = −P (S 2 → S 1 )
P S 1
Si la liaison est parfaite on supposera que le torseur de percussion de liaison à l’instant du choc possède les mêmes caractéristiques que le torseur de liaison habituel.
6.3.3
Choc sans frottement entre deux solides On considère deux solides S1 et S2, qui à l’instant to entrent en contact en un point I.
Le contact étant supposé sans frottement, le torseur des percussions de contact peut s’écrire de la façon suivante: r
− P .n 21 P (S 1 → S 2 )= 0 I Le principe fondamental des chocs appliqué au solide S2 permet alors d’écrire :
[[C (S 2 / Rg )]]t = P (S 1 → S 2) La valeur de P étant inconnue, ceci constitue un système de 6 équations à 7 inconnues. Comme d’habitude on écrira une loi de comportement pour le choc, qui permettra d’éliminer une inconnue. n.b. Le cas des chocs avec frottement est beaucoup plus complexe, et ne sera pas abordé ici.
◊
Définition de e L’équation supplémentaire, est donnée par la condition au point d’impact : r
r
r
n 21 ⋅ V ( I ∈ S 2 / Rg ) − V ( I ∈ S 1 / Rg ) t +
r
r
r
= −e n 21 ⋅ V ( I ∈ S 2 / Rg ) − V ( I ∈ S 1 / Rg ) t −
Qui se note aussi : r r n 21 ⋅ V t + , I , S 2 / S 1
(
r
= −e n 21 ⋅ V (t − , I , S 2 / S 1 r
où e est un coefficient qui dépend des matériaux en présence, de la nature des surfaces, de la température …
◊
Propriété de e et P : 0 ≤ e ≤ 1 et 0 ≤ P
-
Hypothèse des solides indéformables : 0 ≤ e Compte tenu de l’hypothèse des milieux solides indéformables, il ne peut pas y avoir interpénétration des deux solides au point de contact. Dans ce cas, compte tenu de la définition dee : r r n 21 ⋅ V t + , I , S 2 / S 1
r
r
= −e n 21 ⋅ V t − , I , S 2 / S 1
Le coefficient e doit donc être positif ou nul.
119
Préparation à l’agrégation de mécanique : Polycopié de Mécanique des Solides Indéformables APPLICATIONS, TECHNIQUES DE CALCUL
-
Principe des puissances virtuelles : 0 ≤ P On applique le principe des puissances virtuelles au problème du choc, et on emploie comme * torseur des vitesses virtuel {V } la restriction à chacun des solides de la discontinuité du champ des
moments du torseur des vitesses pendant le choc. On obtient alors, en tenant compte du fait que les torseurs des percussions des actions à distance sont nuls :
[[C (S 1 / Rg )]]t × {V (S 1 / Rg )}* = [[P (S 2 → S 1 )]]t × {V (S 1 / Rg )}* [[C (S 2 / Rg )]]t × {V (S 2 / Rg )}* = [[P (S 1 → S 2)]]t × {V (S 2 / Rg )}* Comme par ailleurs :
(
P S 1 → S 2
) = −P (S 2 → S 1)
On trouve en additionnant : 2
∑= [[
(
)}* = [[P (S 2 → S 1 )1]]t × ({V (S 1 / Rg )}* − {V (S 2 / Rg )}* )
)]] { (
C S i / Rg × V S i / Rg t
i 1
On également que :
[[C (S / Rg )]]t × {V (S / Rg )}* =
r
r
∫ [[V ( P / Rg )]]t dm ⋅V * ( P / Rg ) S
Donc : r
r
∫ [[V ( P / Rg )]]t dm ⋅V * ( P / Rg ) = [[P (S 2 → S 1 )]]t × ({V (S 1 / Rg )}* − {V (S 2 / Rg )}* )
S 1∪ S 2
On choisit d’abord, comme champ virtuel, la discontinuité du champ des vitesses des points de chacun des solides r
[[
r
]]
V * ( P / Rg ) = V ( P / Rg ) t Alors : r
r
r
r 2 ∫ [[V ( P / Rg )]]t dm = − P .n12 ⋅ [V (t + , I , S 1 / S 2)− V (t − , I , S 1 / S 2)]
S 1∪ S 2
D’après la loi de choc sans frottement, employée : r
r
r 2 [ [ V ( P / Rg )]]t dm = P (1 + e ).n12 ⋅ V (t − , I , S 1 / S 2) ∫
S 1∪ S 2
r r Comme les deux solides sont entrés en contact, il faut qu’avant le choc n12 ⋅ V t − , I , S 1 / S 2
(
>0
Comme e doit être positif, pour éviter l’interpénétration des surfaces et comme le premier membre est positif, on en déduit que P est nécessairement positif.
-
Principe des puissances virtuelles : e ≤ 1 On applique le principe des puissances virtuelles au problème du choc, et on emploie comme champ de vitesse virtuel, la moyenne de la vitesse avant et après le choc : r
V * ( P / Rg ) =
1
r
[V ( P / Rg )t 2
+
r
+ V ( P / Rg )t − ]
Comme : r
r
∫ [[V ( P / Rg )]]t dm ⋅V * ( P / Rg ) = [[P (S 2 → S 1 )]]t × ({V (S 1 / Rg )}* − {V (S 2 / Rg )}* )
S 1∪ S 2
120