טנזור התמד

טנזור ההתמד (או טנזור האינרציה) הוא דרך נוחה וקצרה להציג את מומנטי ההתמד של הגוף. כתיבת מומנט ההתמד כטנזור מאפשרת למצוא קשר נוח בין המהירות הזוויתית של גוף לתנע הזוויתי ולאנרגיה הקינטית שלו. תנע זוויתי של גוף (אלמנט מסה) ביחס לציר מסוים מוגדר כמכפלה וקטורית של התנע הקווי במרחק מהציר: d L → = r → × d p → = r → × v → d m = d m r → × ( ω → × r → ) {\displaystyle {\vec {dL}}={\vec {r}}\times {\vec {dp}}={\vec {r}}\times {\vec {v}}dm=dm{\vec {r}}\times ({\vec {\omega }}\times {\vec {r}})} כאשר: d L → {\displaystyle {\vec {dL}}} הוא תנע זוויתי של אלמנט מסה p → {\displaystyle {\vec {p}}} התנע הקווי v → {\displaystyle {\vec {v}}} וקטור המהירות d m {\displaystyle dm} אלמנט המסה r → {\displaystyle {\vec {r}}} הוא וקטור המרחק של אלמנט המסה מציר הסיבוב. במערכת צירים קרטזית הווקטור נתון על ידי r → = ( x , y , z ) {\displaystyle {\vec {r}}=(x,y,z)} ω → {\displaystyle {\vec {\omega }}} וקטור המהירות הזוויתית של אלמנט המסה. במערכת צירים קרטזית הווקטור נתון על ידי ω → = ( ω x , ω y , ω z ) {\displaystyle {\vec {\omega }}=(\omega _{x},\omega _{y},\omega _{z})} נוכל לחשב את התנע הזוויתי בדרך של חישוב מכפלות וקטוריות בעזרת דטרמיננטות. ω → × r → = | x ^ y ^ z ^ ω x ω y ω z x y z | = ( ω y z − ω z y , ω z x − ω x z , ω x y − ω y x ) {\displaystyle {\vec {\omega }}\times {\vec {r}}={\begin{vmatrix}{\hat {x}}&{\hat {y}}&{\hat {z}}\\\omega _{x}&\omega _{y}&\omega _{z}\\x&y&z\end{vmatrix}}=(\omega _{y}z-\omega _{z}y,\omega _{z}x-\omega _{x}z,\omega _{x}y-\omega _{y}x)} r → × ( ω → × r → ) = | x ^ y ^ z ^ x y z ω y z − ω z y ω z x − ω x z ω x y − ω y x | = ( ω x y 2 − ω y x y − ω z x z + ω x z 2 , ω y z 2 − ω z y z − ω x y x + ω y x 2 , ω z x 2 − ω x z x − ω y z y + ω z y 2 ) {\displaystyle {\vec {r}}\times ({\vec {\omega }}\times {\vec {r}})={\begin{vmatrix}{\hat {x}}&{\hat {y}}&{\hat {z}}\\x&y&z\\\omega _{y}z-\omega _{z}y&\omega _{z}x-\omega _{x}z&\omega _{x}y-\omega _{y}x\end{vmatrix}}=(\omega _{x}y^{2}-\omega _{y}xy-\omega _{z}xz+\omega _{x}z^{2},\omega _{y}z^{2}-\omega _{z}yz-\omega _{x}yx+\omega _{y}x^{2},\omega _{z}x^{2}-\omega _{x}zx-\omega _{y}zy+\omega _{z}y^{2})} נוכל לפשט את הביטוי: r → × ( ω → × r → ) = ω x ( y 2 + z 2 , − y x , − z x ) + ω y ( − x y , z 2 + x 2 , − z y ) + ω z ( − x z , − y z , x 2 + y 2 ) {\displaystyle {\vec {r}}\times ({\vec {\omega }}\times {\vec {r}})=\omega _{x}(y^{2}+z^{2},-yx,-zx)+\omega _{y}(-xy,z^{2}+x^{2},-zy)+\omega _{z}(-xz,-yz,x^{2}+y^{2})} כעת כדי למצוא את התנע הזוויתי הכולל, נבצע אינטגרל מ0 עד M (כך ש-M היא המסה הכוללת של הגוף המסתובב): L → = ∫ 0 M d m r → × ( ω → × r → ) {\displaystyle {\vec {L}}=\int \limits _{0}^{M}dm{\vec {r}}\times ({\vec {\omega }}\times {\vec {r}})} , וכדי לפתור את האינטגרל יהיה עלינו רק לסכום את המכפלות הווקטוריות עבור כל רכיב בנפרד. L x = ∭ V ω x ( y 2 + z 2 ) ρ ( r ) d V − ∭ V ω y y x ρ ( r ) d V − ∭ V ω z z x ρ ( r ) d V {\displaystyle L_{x}=\iiint \limits _{V}\omega _{x}(y^{2}+z^{2})\rho (r)dV-\iiint \limits _{V}\omega _{y}yx\rho (r)dV-\iiint \limits _{V}\omega _{z}zx\rho (r)dV} L y = − ∭ V ω x x y ρ ( r ) d V + ∭ V ω y ( z 2 + x 2 ) ρ ( r ) d V − ∭ V ω z z y ρ ( r ) d V {\displaystyle L_{y}=-\iiint \limits _{V}\omega _{x}xy\rho (r)dV+\iiint \limits _{V}\omega _{y}(z^{2}+x^{2})\rho (r)dV-\iiint \limits _{V}\omega _{z}zy\rho (r)dV} L z = − ∭ V ω x x z ρ ( r ) d V − ∭ V ω y y z ρ ( r ) d V + ∭ V ω z ( x 2 + y 2 ) ρ ( r ) d V {\displaystyle L_{z}=-\iiint \limits _{V}\omega _{x}xz\rho (r)dV-\iiint \limits _{V}\omega _{y}yz\rho (r)dV+\iiint \limits _{V}\omega _{z}(x^{2}+y^{2})\rho (r)dV} V הוא פונקציית הנפח של הגוף, ו-ρ, פונקציית הצפיפות. ניתן לרשום בצורה מטריציונית: ( L x L y L z ) = ( ∭ V ( y 2 + z 2 ) ρ ( r ) d V − ∭ V y x ρ ( r ) d V − ∭ V z x ρ ( r ) d V − ∭ V x y ρ ( r ) d V ∭ V ( z 2 + x 2 ) ρ ( r ) d V − ∭ V z y ρ ( r ) d V − ∭ V x z ρ ( r ) d V − ∭ V y z ρ ( r ) d V ∭ V ( x 2 + y 2 ) ρ ( r ) d V ) ( ω x ω y ω z ) {\displaystyle {\begin{pmatrix}L_{x}\\L_{y}\\L_{z}\end{pmatrix}}={\begin{pmatrix}\iiint \limits _{V}(y^{2}+z^{2})\rho (r)dV&-\iiint \limits _{V}yx\rho (r)dV&-\iiint \limits _{V}zx\rho (r)dV\\-\iiint \limits _{V}xy\rho (r)dV&\iiint \limits _{V}(z^{2}+x^{2})\rho (r)dV&-\iiint \limits _{V}zy\rho (r)dV\\-\iiint \limits _{V}xz\rho (r)dV&-\iiint \limits _{V}yz\rho (r)dV&\iiint \limits _{V}(x^{2}+y^{2})\rho (r)dV\end{pmatrix}}{\begin{pmatrix}\omega _{x}\\\omega _{y}\\\omega _{z}\end{pmatrix}}} המטריצה באגף השמאלי, היא וקטור התנע הזוויתי, המטריצה האמצעית היא טנזור ההתמד והמטריצה הימנית היא וקטור המהירות הזוויתית. אם מדובר בגופים נקודתיים (לא גוף רציף), יש להחליף את האינטגרלים בסכומים. לפעמים גם כותבים: ( L x L y L z ) = ( I x x I x y I x z I y x I y y I y z I z x I z y I z z ) ( ω x ω y ω z ) {\displaystyle {\begin{pmatrix}L_{x}\\L_{y}\\L_{z}\end{pmatrix}}={\begin{pmatrix}I_{xx}&I_{xy}&I_{xz}\\I_{yx}&I_{yy}&I_{yz}\\I_{zx}&I_{zy}&I_{zz}\end{pmatrix}}{\begin{pmatrix}\omega _{x}\\\omega _{y}\\\omega _{z}\end{pmatrix}}} I = d e f [ I x x I x y I x z I y x I y y I y z I z x I z y I z z ] {\displaystyle \mathbf {I} {\stackrel {\mathrm {def} }{=}}{\begin{bmatrix}I_{xx}&I_{xy}&I_{xz}\\I_{yx}&I_{yy}&I_{yz}\\I_{zx}&I_{zy}&I_{zz}\end{bmatrix}}} הרכיבים על האלכסון הראשי הם מומנטי התמד של המסה ביחס לציר המערכת:   I x x {\displaystyle \ I_{xx}} - הוא מומנט ההתמד של המסה סביב ציר   X {\displaystyle \ X} ומוגדר כ: I x x   = d e f   ∑ i = 1 N m i ( y i 2 + z i 2 ) {\displaystyle I_{xx}\ {\stackrel {\mathrm {def} }{=}}\ \sum _{i=1}^{N}m_{i}(y_{i}^{2}+z_{i}^{2})\,\!}   I y y {\displaystyle \ I_{yy}} - הוא מומנט ההתמד של המסה סביב ציר   Y {\displaystyle \ Y} ומוגדר כ: I y y   = d e f   ∑ i = 1 N m i ( x i 2 + z i 2 ) {\displaystyle I_{yy}\ {\stackrel {\mathrm {def} }{=}}\ \sum _{i=1}^{N}m_{i}(x_{i}^{2}+z_{i}^{2})\,\!}   I z z {\displaystyle \ I_{zz}} - הוא מומנט ההתמד של המסה סביב ציר   Z {\displaystyle \ Z} ומוגדר כ: I z z   = d e f   ∑ i = 1 N m i ( x i 2 + y i 2 ) {\displaystyle I_{zz}\ {\stackrel {\mathrm {def} }{=}}\ \sum _{i=1}^{N}m_{i}(x_{i}^{2}+y_{i}^{2})\,\!} ושאר האיברים הם מכפלות הנקראות מכפלות התמד של המסה ביחס לזוג צירים נתון:   I x y {\displaystyle \ I_{xy}} - הוא מכפלת ההתמד של המסה ביחס לצירים   X , Y {\displaystyle \ X,Y} ומוגדר כ: I x y = I y x   = d e f   − ∑ i = 1 N m i x i y i {\displaystyle I_{xy}=I_{yx}\ {\stackrel {\mathrm {def} }{=}}\ -\sum _{i=1}^{N}m_{i}x_{i}y_{i}\,\!}   I x z {\displaystyle \ I_{xz}} - הוא מכפלת ההתמד של המסה ביחס לצירים   X , Z {\displaystyle \ X,Z} ומוגדר כ: I x z = I z x   = d e f   − ∑ i = 1 N m i x i z i {\displaystyle I_{xz}=I_{zx}\ {\stackrel {\mathrm {def} }{=}}\ -\sum _{i=1}^{N}m_{i}x_{i}z_{i}\,\!}   I y z {\displaystyle \ I_{yz}} - הוא מכפלת ההתמד של המסה ביחס לצירים   Y , Z {\displaystyle \ Y,Z} ומוגדר כ: I y z = I z y   = d e f   − ∑ i = 1 N m i y i z i {\displaystyle I_{yz}=I_{zy}\ {\stackrel {\mathrm {def} }{=}}\ -\sum _{i=1}^{N}m_{i}y_{i}z_{i}\,\!} ומתקיים   I x y = I y x , I y z = I z y , I x z = I z x {\displaystyle \ I_{xy}=I_{yx},I_{yz}=I_{zy},I_{xz}=I_{zx}} . בעזרת טנזור ההתמד נוח להציג את התנע הזוויתי של הגוף בכפל מטריצות: L = I ⋅ ω {\displaystyle L=I\cdot \,\omega } ונוח להציג את האנרגיה הקינטית של הגוף:   E k = 1 2 ω → t I ω → {\displaystyle \ E_{k}={\frac {1}{2}}\,{\vec {\omega }}^{t}I{\vec {\omega }}} . ניתן לעבור מערכות צירים בעזרת הפעלת טרנספורמציה על טנזור ההתמד ללא צורך לחשב את מומנטי ההתמד מחדש במערכת הצירים החדשה. אם   T {\displaystyle \ T} היא מטריצת סיבוב, ו   I {\displaystyle \ I} הוא טנזור ההתמד במערכת צירים ידועה, אז טנזור ההתמד במערכת הצירים המתקבלת על ידי הפעלת   T {\displaystyle \ T} נתון על ידי I n e w = T − 1 ⋅ I ⋅ T {\displaystyle I_{new}=T^{-1}\cdot I\cdot T} מערכת שבה מכפלות ההתמד של הגוף מתאפסות נקראת מערכת צירים ראשיים. ובגלל העובדה שטנזור ההתמד הוא מטריצה סימטרית ניתן ללכסן אותו על ידי מערכת צירים אורתוגונלית. זוהי מסקנה בעלת משמעות פיזיקלית חשובה כי כתוצאה מכך לכל גוף תלת־ממדי ניתן לבחור שלושה צירים ניצבים שיהוו עבורו צירים ראשיים. כיוון שהתנע הזוויתי מתקבל ממכפלת טנזור האינרציה במהירות הזוויתית, מתקבלת תופעה מפתיעה - התנע הזוויתי לא חייב להיות מקביל למהירות הזוויתית. תופעה זו גורמת לכך שגופים המסתובבים באופן חופשי, יכולים לבצע תנועה מסובכת למדי. אם, לדוגמה, נזרוק עט באוויר כך שהוא מסתובב בערך סביב צירו, נגלה כי קצוות העט "מציירים" באוויר מעגלים קטנים. תופעה זו מתקבל כיוון שחוק שימור התנע דורש כי התנע הזוויתי יישאר קבוע. אם המהירות הזוויתית אינה מקבילה לתנע הזוויתי, מוכרח להתקיים שווקטור המהירות הזוויתי יקיף את ווקטור התנע הזוויתי במעגלים. תופעה זו נקראת נקיפה (פרצסיה) ומתוארת מתמטית בעזרת משוואות אוילר תוך שימוש בטנזור ההתמד במערכת צירים ראשית. בטבלה מוצגים טנזורי מומנט ההתמד ביחס לצירים הראשיים של גופים פשוטים נוספים. להצגת הטבלה לחצו על "הצגה". Sybil P. Parker, McGraw Hill Encyclopedia of Engineering, 1983. Irving H. Shames, Engineering Mechanics, Prentic - Hill International Inc. 1970 מומנט האינרציה כטנזור (באנגלית) טנזור האינרציה (באנגלית) מומנט האינרציה (באנגלית) רשימת נוסחאות למומנטי התמד בוויקיפדיה האנגלית