Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

מערך בעיות — אופני תנודה עצמיים ושרשרת N המסות

מערך בעיות: אופני תנודה עצמיים ושרשרת N המסות

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

לאורך כל הדף, שרשרת פירושה NN מסות שוות mm המחוברות בקפיצים זהים בקשיחות ksk_s במרווח aa, כשהקצוות קבועים אלא אם נאמר אחרת, כך ש

ωp=2ksmsin⁡ ⁣(pπ2(N+1)),ap(j)∝sin⁡ ⁣(pπjN+1),c=aksm.\omega_p = 2\sqrt{\frac{k_s}{m}}\sin\!\left(\frac{p\pi}{2(N+1)}\right), \qquad a_p(j) \propto \sin\!\left(\frac{p\pi j}{N+1}\right), \qquad c = a\sqrt{\frac{k_s}{m}} .

בעיה 1 — אימות הניחוש

(א) רשמו את משוואת התנועה של מסה jj בשרשרת עם קצוות קבועים, ואמרו במדויק איזה תפקיד ממלאים התנאים x0=xN+1=0x_0 = x_{N+1} = 0 — הם אינם משוואות תנועה.

(ב) הציבו xj=a(j)cos⁡ωtx_j = a(j)\cos\omega t וצמצמו את הבעיה לנסיגה ב-a(j)a(j) בלבד.

(ג) הראו ש-ap(j)=sin⁡(pπj/(N+1))a_p(j) = \sin(p\pi j/(N+1)) מקיים את שני תנאי הקצה לכל pp שלם, ושהוא מקיים את הנסיגה הפנימית עבור ה-ωp\omega_p שצוטט למעלה. (רמז: sin⁡(θ+φ)+sin⁡(θ−φ)=2sin⁡θcos⁡φ\sin(\theta + \varphi) + \sin(\theta - \varphi) = 2\sin\theta\cos\varphi.)

(ד) הראו ש-pp ו-p+2(N+1)p + 2(N+1) מתארים את אותה תנועה, ושה-p=N+1p = N+1 אינו מתאר תנועה כלל. הסיקו שיש בדיוק NN אופנים ולא יותר, ואמרו מדוע המספר הזה היה חייב לצאת שווה למספר המסות.

בעיה 2 — שלוש מסות ביד, צורות תחילה

קחו N=3N = 3.

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

(ב) רשמו את MM ואת KK ופתחו את det⁡(K−ω2M)=0\det(K - \omega^2 M) = 0. תקבלו משוואה ממעלה שלישית; פרקו אותה לגורמים תוך שימוש בצורה ששרטטתם לאופן האמצעי כרמז לאחד השורשים.

(ג) הראו ששלושת השורשים הם ω=ks/m {2−2, 2, 2+2}\omega = \sqrt{k_s/m}\,\{\sqrt{2 - \sqrt2},\ \sqrt2,\ \sqrt{2 + \sqrt2}\}, ובדקו אותם מול הצורה הסגורה.

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

בעיה 3 — השרשרת החופשית משני הקצוות

הסירו את שני הקירות, ונותרות NN מסות המחוברות ב-N−1N-1 קפיצים.

(א) רשמו את KK והראו שסכום כל שורה הוא אפס. הסיקו ש-a=(1,1,…,1)\mathbf{a} = (1,1,\dots,1) הוא אופן עם ω=0\omega = 0.

(ב) הראו שהקואורדינטה המודאלית של אופן האפס פרופורציונלית לתנע הכולל חלקי המסה הכוללת, ושמשוואת התנועה שלה, q¨0=0\ddot q_0 = 0, היא שימור התנע בניסוח אחר.

(ג) הסבירו מדוע נוסחת אנרגיית האופן 12(q˙p2+ωp2qp2)\tfrac12(\dot q_p^2 + \omega_p^2 q_p^2) עדיין נכונה לאופן הזה ומצטמצמת ל-12q˙02\tfrac12 \dot q_0^2.

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

בעיה 4 — קריאת גרף נפיצה

נמסר לכם גרף של ω\omega מול kk עבור שרשרת לא מוכרת. זו קשת סינוס, אפס ב-k=0k=0, שמגיעה למקסימום של 6.0×10136.0\times10^{13} רדיאן לשנייה ב-k=1.6×1010k = 1.6\times10^{10} מ−1^{-1}.

(א) מצאו את המרווח aa ואת הצירוף ks/m\sqrt{k_s/m}.

(ב) מצאו את מהירות הגל באורכי גל ארוכים cc, והשוו אותה למהירות הקול במוצק טיפוסי, כמה אלפי מטרים לשנייה. הגיבו.

(ג) גל עם k=2×109k = 2\times 10^{9} מ−1^{-1}: מצאו את ω\omega מן היחס המדויק ומן הקירוב הליניארי, ותנו את שגיאת הקירוב באחוזים.

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

בעיה 5 — השרשרת שפרטו עליה כטור פורייה

שרשרת של N=20N = 20 נמשכת למשולש סימטרי — מסה jj מוזזת ב-Amin⁡(j,N+1−j)/N+12A\min(j, N+1-j)/\tfrac{N+1}{2} — ומשוחררת ממנוחה.

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

(ב) מקדמי פורייה של גל משולש יורדים כ-1/p21/p^2 עם סימן מתחלף. אמתו נומרית ש- q1/q3≈−9q_1/q_3 \approx -9 ו-q1/q5≈25q_1/q_5 \approx 25, וזהו את הסטיות הקטנות מ--9 ומ-25 המדויקים — הנמדדות -9.07 ו-25.6 — כתוצאה של דגימת משולש ב-20 נקודות ולא ברציפות.

(ג) חשבו את שיעור האנרגיה הכוללת שנמצא ביסוד. קשרו את תשובתכם לחידה שבראש עמוד המודול.

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

בעיה 6 — גבול הרצף, נמדד

בעזרת coupled.chain_mode_frequencies ו-coupled.chain_continuum_frequencies:

(א) עבור p=1p = 1, חשבו את השגיאה היחסית מול ערך הרצף עבור N∈{10,20,40,80,160}N \in \{10, 20, 40, 80, 160\} והתאימו את סדרה בצירים לוגריתמיים. דווחו את המעריך.

(ב) חזרו על ההתאמה כשאתם משתמשים ב-NN כפרמטר העידון במקום ב-N+1N+1. תקבלו בערך 1.95. הסבירו איזה מהם נכון ומדוע, בלי להיעזר בשאלה איזה מהם קרוב יותר ל-2.

(ג) הניחו את השגיאה החזויה (pπ)2/(24(N+1)2)(p\pi)^2/(24(N+1)^2) מעל המדידות שלכם. דווחו את הפער היחסי הגדול ביותר על פני הרשת שלכם.

(ד) חזרו על (א) עבור p=1,2,4,8,16p = 1, 2, 4, 8, 16 ב-N=100N = 100 קבוע והתאימו את המעריך ב-pp. ציינו את ה-pp הגדול ביותר שעבורו הייתם מוכנים לכנות שרשרת בת 100 מסות “מיתר עד כדי 1%”, ותנו את המספר שמצדיק את הבחירה.

בעיה 7 — אתגר: פגם שלוכד אופן

בנו שרשרת קבועה של N=41N = 41 והחליפו את המסה המרכזית ב-αm\alpha m.

(א) עבור α=0.3\alpha = 0.3, פתרו את האופנים והראו שהתדירות הגבוהה ביותר נמצאת מעל ראש הסרט של השרשרת המושלמת 2ks/m2\sqrt{k_s/m}. דווחו את שני המספרים.

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

(ג) הסבירו מדוע אופן שמעל הסרט אינו יכול להתפשט, ומדוע זו אותה טענה כמו “הוא מקומי”. ההסבר שלכם צריך להשתמש ב-ω(k)\omega(k) ולומר מה kk היה צריך להיות.

(ד) סרקו את α\alpha מ-1 כלפי מטה ומצאו היכן האופן המקומי נפרד מן הסרט. אחר כך חזרו עם α>1\alpha > 1: דווחו אם פגם כבד מייצר אופן מקומי מעל הסרט, מתחתיו, או כלל לא, והסבירו את חוסר הסימטריה.

בעיה 8 — אתגר: פרמי–פסטה–אולם–צינגו

הוסיפו איבר קובי חלש לכל קפיץ, כך שהכוח בין שכנים הופך ל-ks δ (1+β δ)k_s\,\delta\,(1 + \beta\,\delta) כאשר δ\delta היא ההתארכות ו-β\beta קטן.

(א) שנו אינטגרטור ורלה-מהירות כך שיטפל בכוח הזה. אמתו שב-β=0\beta = 0 הוא משחזר את coupled.simulate_coupled.

(ב) התניעו שרשרת של N=32N = 32 באופן הנמוך שלה ועקבו אחר אנרגיות האופנים לאורך ריצה ארוכה, כשאתם מטילים על צורות האופנים הליניאריות לכל אורכה. אמרו מדוע השימוש בצורות הליניאריות לגיטימי כאן ומה היה נשבר אילו β\beta היה גדול.

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

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