Remerciez-le!

Remerciez @Admin pour avoir partagé cet document gratuitement, de la manière la plus simple, en partageant sur les réseaux sociaux.

TS - TP Physique n°11 Eric DAINI – Lycée Paul Cézanne – Aix en Provence - http:

TS - TP Physique n°11 Eric DAINI – Lycée Paul Cézanne – Aix en Provence - http://labotp.org vlim(exp) = 0,94 m.s-1 τexp = 0,085 s RESOLUTION D'UNE EQUATION PAR LA METHODE D'EULER (Correction) Objectifs • Mettre en oeuvre la méthode d'Euler pour déterminer une solution numérique d'une équation différentielle. • Trouver le modèle le plus pertinent pour l'expression de la force de frottement fluide. I. EXPERIENCE • Une bille d'acier tombe verticalement dans de l'huile (voir photo). Une chronophotographie de l'expérience et un pointage Aviméca ont permis d'obtenir le graphe v(t) ci-contre. • La vitesse limite vlim et le temps caractéristique τ τ τ τ de la chute sont rappelés sur le graphe. • Les caractéristiques de la bille d'acier sont : - Masse : m = 4,08 g; Diamètre : d = 10 mm; Volume V. - Masse volumique de la bille : ρ ρ ρ ρ = 7 796 kg.m-3 - Masse volumique de l'huile : ρ ρ ρ ρh = 920 kg.m-3. - Intensité de la pesanteur: g = 9,80 m.s-2 II MODELISATION DE LA CHUTE • Selon le type de force de frottement que l'on considère, l'équation différentielle du mouvement sur la vitesse s'écrit : Cas n°1: f = k1.v B = v A + dt dv 1 (1) Cas n°2 : f = k2.v2 B = v A + dt dv 2 2 (2) avec A = m k et B =       ρ ρ − h 1 g 1) Valeur de la constante B avec 3 chiffres significatifs: B =       ρ ρ − h 1 g = 920 9,80 1 7796   × −     = 8,64 m.s-2. 2) Si v = vlim = 0,94 m.s-1 = cte alors dv dt = 0 : dans l'équation (1): A1.vlim = B ⇔ 1 lim B 8,64 A v 0,94 = = = 9,19 s-1 dans l'équation (2): A2.v²lim = B ⇔ 2 lim B 8,64 A v² (0,94)² = = = 9,78 m-1 3) Expressions numériques des deux équations différentielles: dv 9,19.v 8,64 dt + = (1) 2 dv 9,78.v 8,64 dt + = (2) III RESOLUTION NUMERIQUE D'UNE EQUATION DIFFERENTIELLE SELON LA METHODE D'EULER 1) La méthode d'Euler • Principe : pour des intervalles de temps ∆t "petits" on peut faire l'approximation suivante: a(t) = dt dv ≈ v v(t t) v(t) t t ∆ + ∆ − = ∆ ∆ soit : v(t+∆t) = v(t) + a(t) × × × × ∆t. ∆ ∆ ∆ ∆t est appelé pas d’itération de la méthode d’Euler. L’expression de a(t) dépend du modèle utilisé. Dans la suite on choisit : ∆ ∆ ∆ ∆t = 0,010 s. vlim(exp) = 0,94 m.s-1 τexp= 0,085 s vlim τ vlim τ TS - TP Physique n°11 Eric DAINI – Lycée Paul Cézanne – Aix en Provence - http://labotp.org 2) Calculs "à la main" • Cas n° 1 : modélisation en f = k 1.v: La bille est lâchée sans vitesse initiale à la date t0 = 0,000 s . a) a1(t) = dt dv = B - A1.v = 8,64 – 9,19 × × × × v(t) b) 3) Résolution avec un tableur • Les calculs fastidieux "à la main" peuvent être automatisés avec un ordinateur et un tableur, tel qu'Excel . Cas n°1 : modélisation de f en f = k 1.v a) Formule de calcul pour calculer t1 dans la cellule A10: = A9+B$3 formule de calcul pour calculer v(t1) dans la cellule B10: = B9+(B$5-B$4*B9)*B$3 b) Graphe de v(t) obtenu par la méthode d'Euler et les points expérimentaux: c) Le modèle en f = k1.v n'est pas satisfaisant car le graphe v(t) obtenu avec la méthode d'Euler est situé au–dessous des points expérimentaux. t(s) v(t) en (m.s-1 ) a1(t) en (m.s-2 ) t0 = 0,000 8,640 0,000 t1 = to+ ∆t = 0,010 7,846 0,086 t2 = t1+ ∆t = 0,020 7,125 0,165 t3 =t2 + ∆t = 0,030 6,470 0,236 t4 =t3 + ∆t = 0,040 5,876 0,301 Cas n° 2 : modélisation en f = k 2.v² a) a2(t) = dt dv = B - A1.v = 8,64 – 9,78 × × × × v(t) b) t(s) v(t) en (m.s-1 ) a2(t) en (m.s-2 ) t0 = 0,000 0,000 0,000 t1 = to+ ∆t = 0,010 0,086 0,086 t2 = t1+ ∆t = 0,020 0,172 0,172 t3 =t2 + ∆t = 0,030 0,256 0,256 t4 =t3 + ∆t = 0,040 0,336 0,336 TS - TP Physique n°11 Eric DAINI – Lycée Paul Cézanne – Aix en Provence - http://labotp.org Cas n°2 : modélisation de f en f = k 2.v2 a) formule de calcul pour calculer t1 dans la cellule A10: = A9+B$3 formule de calcul pour calculer v(t1) dans la cellule B10: =B9+(B$5-B$4*B9^2)*B$3 b) graphe de v(t) obtenu par la méthode d'Euler et les points expérimentaux. c) Le modèle f = k2.v2 n'est pas satisfaisant car le graphe v(t) obtenu avec la méthode d'Euler passe au-dessus de certains points expérimentaux. Cas n°3 : modélisation de f en f = k 3 . vn • Les deux modélisations précédentes n'étant pas satisfaisantes, on envisage une troisième possibilité : f = k3.vn. L'équation différentielle correspondante s'écrit : Cas n°3 : f = k3.vn B = v A + dt dv n 3 a) Pour le modèle f = k1.v, les points expérimentaux sont au-dessus du graphe v(t) tandis que pour le modèle f = k2.v² certains points expérimentaux sont au-dessous du graphe v(t). On peut donc conclure que n est compris entre 1 et 2 dans le modèle: f = k3.vn b) Pour v = vlim = 0,94 m.s-1 = cte alors dv dt = 0 : équation (3): A3.vn lim = B ⇔ 3 n lim B A v = c) formule de calcul pour A3 dans la cellule B4: = B$5/0,94^B6 formule de calcul pour calculer t1 dans la cellule A11: = A10+B$3 formule de calcul pour calculer v(t1) dans la cellule B11: = B10+(B$5-B$4*B10^B$6)*B$3 d) La valeur de n avec le curseur qui donne la modélisation la plus satisfaisante est: n = 1,5 e) Graphe de v(t) obtenu par la méthode d'Euler et les points expérimentaux avec n = 1,5. TS - TP Physique n°11 Eric DAINI – Lycée Paul Cézanne – Aix en Provence - http://labotp.org f) La modélisation de f la plus satisfaisante est: f = k3.v1,5 avec: k3 = m.A3 = 4,08.10-3 × 9,48 = 3,87.10-2 (S.I) finalement: f = 3,87.10-2 × × × × v1,5 g) Lorsque le régime permanent est atteint, v = vlim = 0,94 m.s-1 donc: f = 3,87.10-2 × 0,94 = 3,5.10-2 N h) Poids de la bille: P = m.g = 4,08.10-3 × 9,80 = 4,0.10-2 N f + PA en régime permanent Poussé d'Archimède: PA = ρh × V × g = 920 × (4π × (5,0.10-3)3 / 3) × 9,80 = 4,7.10-3 N f + PA = 3,5.10-2 + 4,7.10-3 = 3,97.10-2 N ≈ 4,0.10-2 N = P Le résultat pouvait être prévu car en régime permanent le vecteur vitesse est constant donc: P + f + PA = 0 En projection suivant un axe vertical orienté vers le bas: P – f – PA = 0 ⇔ P = f + PA. 4) Influence du pas ∆ ∆ ∆ ∆t sur la résolution de l'équation différentielle par la méthode d'Euler a) Graphes de v(t) pour ∆ ∆ ∆ ∆t = 0,010 s et pour ∆ ∆ ∆ ∆t = 0,100 s: b) Pour ∆t > 0,100 s, le graphe v(t) s'éloigne nettement des points expérimentaux. Pour être applicable, la méthode d'Euler doit utiliser un pas d'intégration ∆t "pas trop grand". c) Pour une bonne modélisation ∆t = 0,010 s. La constante de temps est : τ τ τ τ = 0,085 s On a donc: τ / ∆t = 8,5 ⇔ ∆t = τ / 8,5 ≈ τ / 10. Une condition d'utilisation de la méthode d'Euler sur ∆ ∆ ∆ ∆t est donc: ∆ ∆ ∆ ∆t ≤ τ τ τ τ / 10. TS - TP Physique n°11 Eric DAINI – Lycée Paul Cézanne – Aix en Provence - http://labotp.org IV APPLICATION DE LA METHODE D'EULER A LA CHARGE D'UN CONDENSATEUR • On considère un circuit (RC) série alimenté par un générateur de tension constante E. 1) Schéma du circuit (RC) et équation différentielle sur la tension uc(t) aux bornes du condensateur 2) Exprimons uc(t+∆t): Pour des intervalles de temps ∆t "petits" on peut faire l'approximation suivante: c du dt ≈ c u t ∆ ∆ On a: c du dt ≈ c u t ∆ ∆ = E τ - c u τ or ∆uc = uc(t+∆t) – uc(t) donc il vient: uc(t + ∆t) = uc(t) + c E u −     τ   × ∆t • Dans la cellule B10 on tape: =B9+H$5*(H$3-B9)/H$4 • Sachant qu'une solution de l'équation différentielle est uc(t) = E.(1 – exp(-t/τ)), uploads/Sante/ 10-euler-correction.pdf