הנדסת מכונות מתקדמת: מדריך מעשי לשיטת האלמנטים הסופיים (FEA) עם MATLAB

מבוא לשיטת האלמנטים הסופיים בהנדסת מכונות

שיטת האלמנטים הסופיים (FEA) היא כלי חישובי רב-עוצמה המאפשר למהנדסים לנתח התנהגות של מבנים מורכבים בתנאי עומס שונים. בקורס MECH0059 Advanced Computer Applications in Engineering, סטודנטים להנדסת מכונות נדרשים לפתח תוכנית MATLAB המבוססת על אלמנטים מרובעים בעלי ארבעה צמתים (4-noded quadrilateral elements) לצורך חישוב עיבור ומאמצים בלוח דק. מדריך זה נועד להסביר את העקרונות המרכזיים, תוך שימת דגש על יישומים מעשיים והבנת התהליך.

עקרונות יסוד של שיטת האלמנטים הסופיים

שיטת האלמנטים הסופיים מבוססת על פירוק מבנה רציף למספר סופי של אלמנטים קטנים, המחוברים זה לזה בנקודות הנקראות צמתים (nodes). עבור כל אלמנט מחשבים את מטריצת הקשיחות (stiffness matrix), המתרגמת עיבור לכוחות. בהנחה של מאמץ מישורי (plane stress) עם מודול אלסטי E=40 GPa, יחס פואסון ν=0.3 ועובי לוח 2 מ"מ, נדרש לחשב את מטריצת הקשיחות האלמנטרית. דוגמה טרנדית: בדומה לאופן שבו רשת נוירונים בבינה מלאכותית מפצלת תמונה לחלקיקים קטנים לצורך זיהוי עצמים, גם FEA מפצלת מבנה לאלמנטים קטנים לניתוח מדויק.

בניית תוכנית MATLAB לחישוב מטריצת קשיחות

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

function k = elementStiffness(E, nu, t, coords)
% coords: 4x2 matrix [x1 y1; x2 y2; x3 y3; x4 y4]
% E: Young's modulus, nu: Poisson's ratio, t: thickness
% Returns 8x8 stiffness matrix

% Jacobian matrix and integration points (2x2 Gauss quadrature)
% ... (code for shape functions, B matrix, D matrix)
k = zeros(8,8);
for i = 1:2
  for j = 1:2
    xi = [-1/sqrt(3), 1/sqrt(3)];
    eta = [-1/sqrt(3), 1/sqrt(3)];
    [B, detJ] = computeB(xi(i), eta(j), coords);
    D = (E/(1-nu^2)) * [1 nu 0; nu 1 0; 0 0 (1-nu)/2];
    k = k + B' * D * B * detJ * t;
  end
end
end

הקוד משתמש באינטגרציה נומרית בשיטת גאוס-לגנדר (2x2 Gauss quadrature) לחישוב מדויק של מטריצת הקשיחות. יש לוודא שהקואורדינטות של הצמתים מוזנות בסדר הנכון (בכיוון השעון או נגדו).

ניתוח מבנה עם 2 ו-4 אלמנטים

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

בחינת השפעת מודול האלסטיות ותנאי העומס

בחלק השלישי, יש להריץ את התוכנית עבור שלושה ערכי מודול אלסטי שונים (למשל 30, 50, 60 GPa) ועבור שלושה כיווני עומס שונים (זוויות θ שונות, למשל 45°, 60°, 90°). יש להשוות את תוצאות העיבור והמאמצים ולנתח את המגמות. לדוגמה, הגדלת מודול האלסטיות תקטין את העיבור עבור אותו עומס, בעוד ששינוי זווית העומס תשנה את התפלגות המאמצים. ניתוח זה חיוני להבנת השפעת תכונות החומר ותנאי הגבול על התנהגות המבנה.

אימות התוצאות באמצעות תוכנה מסחרית (ANSYS)

כדי לאמת את נכונות התוכנית, יש למודל את אותו מבנה בתוכנת ANSYS Mechanical APDL ולהשוות את התוצאות. ההשוואה צריכה לכלול את העיבורים והמאמצים בצמתים קריטיים. סטייה קטנה (פחות מ-5%) מעידה על תוכנית מדויקת. דוגמה טרנדית: בדומה לאימות מודל בינה מלאכותית באמצעות מערך בדיקה (test set) לפני פריסה בייצור, גם כאן יש לאמת את הקוד מול תוכנה מוכרת.

הכנת מצגת PowerPoint לסיכום העבודה

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

טיפים להצלחה בוויווה (Viva)

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

סיכום

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