FEM1_006 שיטת אלמנטים סופיים בתלת-ממד

מבוא

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

כאשר הוא השארית ו- היא פונקציית מבחן.

תפוח אדמה הנדסי. (Zohdi, 2018).

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

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

מוביל ל-,

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

אשר, מכיוון שההטרחה , מוביל ל-

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

למה ? זה לא אמור להיות ?

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

לפיכך, נוכל לרשום:

אם נחליט להגביל את הבחירות שלנו של לאלה כך ש-, יש לנו,
כאשר הוא התזוזה המוחלת על , עבור אלסטיות ליניארית עם עיבור אינפיניטסימלי, הבעיה בצורתה החלשה היא:

מצאו , , כך ש-, :

כאשר על .

כמו בניסוח החד-ממדי, זה נקרא “צורה חלשה” מכיוון שהיא אינה דורשת את הגזירות של המאמץ . במילים אחרות, דרישות הגזירּות נחלשו. נדגיש שאם נבדוק את הפתרון עם כל פונקציות המבחן האפשריות של חלקות מספקת, אז הפתרון החלש שקול לפתרון החזק. יתר על כן, בהינתן שהפתרון האמיתי חלק מספיק, הצורות החלשה והחזקה הן שקולות.

מרחבים סובולב הילברטיים

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

באותו אופן כמו בחד-ממד, אנו מגדירים את כמרחב של פונקציות וקטוריות שרכיביהם ב- , כלומר:

לפיכך, אנו מנסחים את הבעיה החלשה בצורה הבאה:

מצאו , , כך ש-, :

קירוב אלמנטים סופיים (FEM)

נוח לרשום את התבנית הבילינארית בצורה מטריציונית באופן הבא:

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

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

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

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

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

או , כאשר, לדוגמה¹:

נוח לרשום:

וכן . אם נבחר את עם אותו בסיס, אך עם צירוף לינארי שונה , אז נוכל לכתוב:

מכיוון ש- שרירותי (כלומר ), נקבל:

באופן מפורש, היא מערכת המשוואות שיש לפתור.

טרנספורמציות גלובליות/מקומיות

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

וכן

בפירוק החישובים לאלמנטים, , כאשר (מספר האלמנטים), כאשר

וכן

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

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

אלמנט סופי דו-ממדי. (Zohdi, 2018).

יצירת רשת ופונקציות קישוריות

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

מספר אלמנט	צומת גלובלי 1	צומת גלובלי 2	צומת גלובלי 3	צומת גלובלי 4
1	1	2	3	4
2	2	3	4	5
3	7	8	9	10
4	6	7	8	9

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

מגבלות על אלמנטים

נזכיר שבחד-ממד יש לנו חישוב מהסוג הבא:

ברור שיעקוביאן אפס יוביל לבעיות ולאינטגרלים שעלולים להיות סינגולריים. בחד-ממד, היה קל להימנע מכך מכיוון שהצמתים ממוספרים באופן סדרתי, וכל עוד הצמתים אינם מתלכדים, הדבר לא יקרה, מאחר ש-. עם זאת, ברור ש- אינו פיזיקלי, מכיוון שמשמעותו היא שצמתים שכנים ממופים ‘הפוך’ (אחד דרך השני). יעקוביאנים שליליים יכולים גם להוביל למטריצות קשיחות לא מוגדרות (indefinite). כמו בניסוחים חד-ממדיים, בניסוחים דו-ממדיים ותלת-ממדיים, יש לוודא ש- בכל התחום. ישנן שתי דרכים שבהן יכולים להיווצר יעקוביאנים לא-חיוביים: (1) האלמנטים אינם קמורים (non-convex) ו-(2) מספור הצמתים שגוי, מה שגורם לאלמנטים ‘להימשך הפוך’. עלינו להבטיח ש-, מכיוון של- יש משמעות פיזיקלית: הוא היחס בין הנפח הדיפרנציאלי של אלמנט המסטר לנפח הדיפרנציאלי של האלמנט הסופי. אם הצמתים ממוספרים נכון כדי להבטיח שהם לא ‘נמשכים הפוך’ (לדוגמה, ראו את האיור הבא) ושהאלמנטים קמורים, אז .

דוגמה לרשת ממופה עבור רצועה חצי מעגלית.

אלמנטים טובים ורעים: דוגמאות

נתבונן באלמנט לינארי דו-ממדי. פונקציות הצורה הן:

פונקציות המיפוי הן:

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

באופן מפורש:

אלמנט לינארי דו-ממדי ודוגמאות למיפוי. (Zohdi, 2018).

עבור ארבעת המקרים באיור לעיל, יש לנו:

• מקרה 1: אלמנט זה קביל, מאחר ש- בכל האלמנט. היעקוביאן קבוע במקרה זה.
• מקרה 2: אלמנט זה אינו קביל, מאחר ש- בכל האלמנט. הבעיה המהותית היא שהצמתים ממוספרים באופן שגוי, מה שהופך את האלמנט ל-‘הפוך’.
• מקרה 3: אלמנט זה קביל, מאחר ש- בכל האלמנט. בעוד שהיעקוביאן אינו קבוע בכל תחום האלמנט, הוא חיובי וחסום.
• מקרה 4: אלמנט זה אינו קביל, מאחר ש- באזורים מסוימים של האלמנט. למרות שהיעקוביאן חיובי בחלקים מסוימים של התחום, יעקוביאן שלילי בחלקים אחרים יכול לגרום לבעיות, כגון סינגולריות פוטנציאלית במטריצת הקשיחות.
• עבור אלמנטים לינאריים, המדד המרכזי לאלמנט בעייתי הוא אי-קמירות (non-convexity) של האלמנט (גם אם הוא ממוספר נכון).

הערות שוליים

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