מטוטלת מתמטית

מתוך ויקיפדיה, האנציקלופדיה החופשית
קפיצה אל: ניווט, חיפוש
Simple pendulum.jpg

מטוטלת מתמטית (נקראת גם מטוטלת פשוטה) היא מטוטלת שמורכבת מגוף בעל ממדים קטנים, התלוי על חוט שמסתו ומידת ההתארכות שלו בזמן התנודות ניתנים להזנחה. בנוסף, זווית התנודה של המטוטלת קטנה יחסית. זהו מודל פיזיקלי, שאינו קיים באופן מושלם במציאות, אך בזכות הקירוב ניתן לתאר את תנועת הגוף באופן פשוט. תחת קירוב זה מטוטלת מתמטית היא סוג של אוסצילטור הרמוני, ולכן מהווה מודל לתופעות פיזיקליות רבות.

לאחר הסטת הגוף מנקודת שיווי המשקל ושיחרורו, הגוף יבצע תנודות הרמוניות במישור אנכי סביב נקודת שיווי המשקל-לאורך קשת מעגל. זמן המחזור של תנודות מטוטלת מתמטית אינו תלוי במסת המשקולת ובמשרעת (זווית) התנודות. זוהי תכונה חשובה, שכן היא מאפשרת למדוד מרווחי זמן. בעבר (עד המצאת השעון החשמלי) השתמשו בני האדם בשעון המונע על ידי מטוטלת.


ניתוח מתמטי[עריכת קוד מקור | עריכה]

אנימציה של תנועת מטוטלת, המראה את השתנות וקטורי המהירות (V) והתאוצה (A).
ניתוח הכוחות על מטוטלת מתמטית

ננתח את תנועת המטוטלת, בקירוב בו החוט חסר מסה ואורכו קבוע, מסת המטוטלת נקודתית, וזווית התנודה קטנה. נסמן:

  •  \ l - אורך החוט
  •  \ m - מסת המשקולת
  •  \ g - תאוצת הכובד
  •  \ \theta - הזווית מהאנך.

ננתח את המומנטים הפועלים על המטוטלת, יחסית לנקודת התלייה. מכיוון שהחוט מחובר לנקודת התלייה, הוא אינו מפעיל מומנט. לכן, המומנט היחיד הפועל על המטוטלת הוא המומנט שמפעיל כוח הכובד, וגודלו  \ -mgl \sin \theta . (הסימן שלילי כיוון שזהו כוח מחזיר, המנוגד לכיוון ההעתק).

מומנט ההתמד של המערכת הוא פשוט I = ml^2, ולכן מתקיים

 \ -mgl \sin \theta = ml^2 \ddot \theta .

באמצעות קירוב זוויות קטנות,  \ \sin \theta \approx \theta , נקבל משוואה של אוסצילטור הרמוני:

\ \ddot \theta = - \frac{g}{l} \theta

הפתרון הכללי של משוואה זו הוא:

 \ \theta = \theta_0 \sin (\omega t + \phi),

כאשר  \omega = \sqrt {\frac{g}{l}}, ו  \ \theta_0 , \phi נקבעים על ידי תנאי ההתחלה. זוהי פונקציה מחזורית, בתדירות  f= \frac{1}{2\pi} \sqrt{ \frac{g}{l}}, ובזמן מחזור  \ T=2 \pi \sqrt{ \frac{l}{g}}.

המרחק מנקודת שיווי המשקל הוא \ x= l \theta (שוב, בקירוב של זוויות קטנות), ולכן גם המרחק מקיים תנודה הרמונית:  \ x = x_0 \sin (\omega t + \phi).

באופן כללי, ניתן לפתור את המשוואה פתרון אנליטי גם ללא הקירוב של זוויות קטנות, בעזרת אינטגרל אליפטי.

במערכת זו, כאשר H הוא המילטוניאן וכאשר H=\frac{g}{l} שרטוט הפאזה יתן הפרדה בין שני אזורים. הפרדה זו נקראת ספרטריקס.

פתרון מדויק[עריכת קוד מקור | עריכה]

הסטייה של זמן המחזור בקירוב הזוויות הקטנות מזמן המחזור המדויק. ניתן לראות כי עד לזווית מפתח של 30°, מתקבל דיוק טוב מאוד - פחות מ 2%.

כאמור, ניתן למצוא את זמן המחזור של מטוטלת פשוטה גם ללא ההנחה של זוויות קטנות. יש להדגיש כי פתרון זה אינו אנליטי (כלומר, לא ניתן להציגו כהרכבה של פונקציות אלמנטריות: פולינומים, אקספוננטים ופונקציות טריגונומטריות). הפתרון מערב אינטגרל אליפטי, שאותו יש לחשב באופן נומרי.

משימור אנרגיה מתקבל מיידית כי \ v=\sqrt{2g\Delta h}. כאשר \ \Delta h הוא ההפרש בין גובה המטוטלת לגובהה המקסימלי. אנו יודעים כי  v=l\frac{d\theta}{dt} ולכן

\frac{d\theta}{dt}=\frac{\sqrt{2g\Delta h}}{l}.

מצד שני, משיקול גאומטרי ניתן לראות ש \ \Delta h=l(\cos\theta-\cos\theta_0), כאשר  \ \theta_0 היא הזווית שבה המטוטלת נמצאת בגובהה המקסימלי. מתקבלת המשוואה הדיפרנציאלית

 \frac{d\theta}{dt}=\sqrt{\frac{2g}{l}(\cos\theta-\cos\theta_0)}.

זוהי משוואה פרידה, וצריך לבצע את האינטגרל

\int dt=\int\frac{d\theta}{\sqrt{\frac{2g}{l}(\cos\theta-\cos\theta_0)}}.

ניתן לקבל את זמן המחזור על ידי אינטגרציה על רבע מחזור והכפלה ב-4:

 T=4\int_{\theta_0}^0\frac{d\theta}{\sqrt{\frac{2g}{l}(\cos\theta-\cos\theta_0)}}

אינטגרל זה אינו פתיר אנליטית, אך ניתן להביעו באמצעות  \ F(k,\phi) האינטגרל האליפטי של לז'נדר מהסוג הראשון, כך:

 T=4\sqrt{\frac{l}{g}}F(\sin\frac{\theta_0}{2},\frac{\pi}{2}) .

הרחבה - ניתוח האופנים העצמיים למטוטלת בעלת N מתנדים צמודים[עריכת קוד מקור | עריכה]

Postscript-viewer-shaded.png ערך מורחב – מטוטלות מתמטיות צמודות

מטוטלת בעלת N מתנדים צמודים, היא בעצם מספר מטוטלות הקשורות זו לזו כאשר החוט של הראשונה יוצא מסוף המטוטלת הקודמת. מספר האופנים העצמיים של מטוטלות מתמטיות צמודות, הוא בדיוק כמספר המתנדים הצמודים שבה. האופנים העצמיים מצביעים על אורך פרקי הזמן שבהם המערכת חוזרת לאותו מבנה מרחבי, ללא חשיבות לזווית ביחס לאנך לרצפה, אלא רק לזוויות בתוך המערכת. עם הוספת מתנד צמוד למטוטלת, האופנים העצמיים שלה מתפלגים באופן הבא: מספר האופנים העצמיים של המערכת גדל והפרש הערכים בין אחד למשנהו קטן. בנוסף לכך, הם הולכים ומתבדרים.

ניתוח מתמטי עבור שני מתנדים צמודים[עריכת קוד מקור | עריכה]

DoublePendulum.JPG

על מנת לחשב את משוואות התנועה עבור מטוטלת בעלת שני גופים, נשתמש במשוואות אוילר-לגראנז'. לפי משוואות אלו:

 \ T = \frac{1}{2} m_1 v_1^2 + \frac{1}{2} m_2 v_2^2

כאשר: T - האנרגיה הקינטית האצורה במערכת.

U, אנרגיית הגובה האצורה במערכת היא:  \ U = -g( R_1 \cos\phi (m_1+m_2)+ m_2 R_2 \cos \psi )

פונקציית הלגראנז' L של המערכת היא:  \ L = T - U ; <נחשב עתה את משוואת אוילר-לגראנז' לפי הזווית \ \phi : \ \frac{d}{dt} \frac {\partial L}{\partial \dot \phi} = \frac {\partial L}{\partial \phi} לאחר חישוב הנגזרות החלקיות בכל אגף מתקבלת המשוואה:  \ (I): m_1 R_1^2 \ddot \phi + m_2 R_2^2 \ddot \phi + m_2 R_1 R_2 \ddot \psi \cos (\phi-\psi)- m_2 R_1 R_2 \dot \psi \sin(\phi-\psi)(\dot \phi - \dot \psi)  \ = -m_2 R_1 R_2 \dot \phi \dot \psi \sin(\phi-\psi)- gR_1 \sin \phi (m_1+m_2 )

באותה צורה, מחישוב משוואת אוילר-לגראנז' לפי הזווית \ \psi מתקבלת המשוואה:

\ (II): m_2 R_2^2 \ddot \psi + m_2 R_1 R_2 \ddot \phi \cos(\phi-\psi)- m_2 R_1 R_2 \dot \phi \sin(\phi-\psi)(\dot \phi - \dot \psi)  \ = m_2 R_1 R_2 \dot \phi \dot \psi \sin(\phi-\psi)- g m_2 R_2 \sin \psi

לאחר קירוב לזווית קטנות וארגון המשוואות, מקבלים מ-(I):

\ (2,1): R_1 \ddot \phi + \frac{m_2}{m_1 + m_2} R_2 \ddot \psi = -g\psi

הערה: משמעות הסימונים למשוואות (n,k) הינו: המשוואה ה-k עבור n מתנדים צמודים. מ-(II) מקבלים:

\ (2,2): R_1 \ddot \phi + R_2 \ddot \psi = -g\psi

ניתוח מתמטי עבור שלושה מתנדים צמודים[עריכת קוד מקור | עריכה]

TriplePendulum.JPG

מחישוב משוואות אוילר לגראנז' בצורה זהה לזו שפורטה לעיל, קירוב לזוויות קטנות, וארגון המשוואות מתקבלות המשוואה הבאות:

\ (3,1): R_1 \ddot \phi + \frac {m_2+ m_3}{m_1+m_2+m_3} R_2 \ddot \psi+\frac {m_3}{m_1+m_2+m_3} R_3 \ddot \theta = -g\phi

\ (3,2): R_1 \ddot \phi+ R_2 \ddot \psi + \frac {m_3}{m_2+m_3} R_3 \ddot \theta = -g\psi

\ (3,3): R_1 \ddot \phi + R_2 \ddot \psi + R_3 \ddot \theta = -g\theta

הכללה ל-N מתנדים צמודים[עריכת קוד מקור | עריכה]

בעזרת שש המשוואות שפותחו לעיל, ננסח עתה נוסחה כללית עבור מספר (k) של מתנד במטוטלת בעלת N מתנדים:  \ (N,k) \sum_{i=1}^K R_i \ddot \phi_i + \sum_{i=k+1}^N \ddot \phi_i R_i \frac {\sum_{j=i}^N m_i}{\sum_{j = k}^N m_i} = -g\phi_k

כאשר: \ \phi _i היא הזווית ה-i מלמעלה

הערה: בנוסחה לעיל, כאשר משתמשים במשוואה עבור k=n יש להשמיט בחישוב את האיבר האחרון (המשוואה אמורה להתחיל מ k+1 ולעלות עד N דבר אשר אינו אפשרי במקרה זה).

ראו גם[עריכת קוד מקור | עריכה]

קישורים חיצוניים[עריכת קוד מקור | עריכה]