Summary

בהתבסס ידע ענן FE סימולציה של תהליכים פורמינג פח

Published: December 13, 2016
doi:

Summary

The following paper presents a novel FE simulation technique (KBC-FE), which reduces computational cost by performing simulations on a cloud computing environment, through the application of individual modules. Moreover, it establishes a seamless collaborative network between world leading scientists, enabling the integration of cutting edge knowledge modules into FE simulations.

Abstract

The use of Finite Element (FE) simulation software to adequately predict the outcome of sheet metal forming processes is crucial to enhancing the efficiency and lowering the development time of such processes, whilst reducing costs involved in trial-and-error prototyping. Recent focus on the substitution of steel components with aluminum alloy alternatives in the automotive and aerospace sectors has increased the need to simulate the forming behavior of such alloys for ever more complex component geometries. However these alloys, and in particular their high strength variants, exhibit limited formability at room temperature, and high temperature manufacturing technologies have been developed to form them. Consequently, advanced constitutive models are required to reflect the associated temperature and strain rate effects. Simulating such behavior is computationally very expensive using conventional FE simulation techniques.

This paper presents a novel Knowledge Based Cloud FE (KBC-FE) simulation technique that combines advanced material and friction models with conventional FE simulations in an efficient manner thus enhancing the capability of commercial simulation software packages. The application of these methods is demonstrated through two example case studies, namely: the prediction of a material’s forming limit under hot stamping conditions, and the tool life prediction under multi-cycle loading conditions.

Introduction

Finite Element (FE) simulations have become a powerful tool for optimizing process parameters in the metal forming industry. The reliability of FE simulation results is dependent on the accuracy of the material definition, input in the form of flow stress data or constitutive equations, and the assignment of the boundary conditions, such as the friction coefficient and the heat transfer coefficient. In the past few years, advanced FE simulations have been developed via the implementation of user-defined subroutines, which have significantly broadened the capability of FE software.

The use of such advanced FE simulations in the design of forming processes for structural components has been investigated by both the aviation and automotive industries, with the intention of producing lightweight structures that reduces operating costs and CO2 emissions. Particular focus has been placed on the replacement of steel components with lower density materials, such as aluminum alloys and magnesium alloys. However, these alloys, especially the stronger variants, offer limited formability at room temperature and thus complex-shaped components cannot be manufactured using the conventional cold stamping process. Therefore, advanced high temperature forming technologies, such as warm aluminum forming 1-4, hot stamping of aluminum alloys 5-9 and hot stamping of high strength steels 10, have been developed over the past decades to enable complex-shaped components to be formed. In general, high temperature forming processes involve significant temperature variations, strain rate and loading path changes 11, which would, for instance, cause inevitable viscoplastic and loading history dependent responses from the work piece materials. These are intrinsic features of high temperature forming processes and may be difficult to represent using conventional FE simulation techniques. Another desirable feature would be the ability to predict the tool life over multiple forming cycles in such processes, since they require low friction characteristics achieved through coatings that degrade with each forming operation. To represent all these features via the implementation of user-defined subroutines would be computationally very expensive. Moreover, the development and implementation of multiple subroutines would require excessive multi-disciplinary knowledge from an engineer conducting the simulations.

In the present work, a novel Knowledge Based Cloud FE (KBC-FE) simulation technique is proposed, based on the application of modules on a cloud computing environment, that enables an efficient and effective method of modeling advanced forming features in conjunction with conventional FE simulations. In this technique, data from the FE software is processed at each cloud module, and then imported back into the FE software in the relevant consistent format, for further processing and analysis. The development of these modules and their implementation in the KBC-FE is detailed.

Protocol

1. פיתוח של טמפרטורה גבוהה יצירת מודל חיזוי הגבל לייזר לחתוך את דגימות לבדיקות formability מבין הסדינים AA6082 סגסוגת אלומיניום (עובי 1.5 מ"מ) לתוך הגיאומטריות נבחרים 12. Etch דפוס רשת, מורכב נקודות עגולות 0.75 מ"מ קוטר עם מרווח קבוע של 1 מ"מ, על פני השטח של הדגימות באמצעות שיטה אלקטרוליטי 13. ידני להחיל גרפיט גריז כחומר סיכה בצד הלא חרוט. הרכב את אסדת כיפת המבחן במכבש הידראולי גבוה שיעור 12. השתמש מכונת בדיקה אוניברסלית הידראולי 250 kN. חמם את אסדת כיפת מבחן לטמפרטורת בדיקות ולהגדיר את האגרוף במהירות מרגשת מתמיד. אז ליזום את הבדיקה. הערה: טמפרטורות הבדיקה הן 300, 400, ו -450 ° C, בהתאמה. מהירויות הבדיקות לכלול 75, 250, ו -400 מ"מ / s. עצור את הבדיקה בבית המופע הראשון של מזמוזים. הערה: עיתונות Strokדואר (כלומר, גובה דגימה סופי) מוגדר כך מזמוזים רק הוא ציין על הדגימה נוצרה. מדוד את גובה הדגימה הסופי באמצעות מד גובה, ולחשב את זני שיעורי זן מרביים (שיעור השינוי של זן ביחס לזמן) באמצעות 3D אופטי להרכיב מערכת ניתוח. לנתח את השינויים ריווח לרשת כדי לחשב את הזנים בכל נקודה של הדגימה נוצרה. ודא 3D האופטי להרכיב מערכת ניתוח כולל מצלמה, דגימה נוצרה, וברים בקנה מידת כיול 14. הערה: הדגימה מועברת במרכז פטיפון וסגורה עם ברי הסולם, ועמדות ביחס שלהם נשארו קבועות למשך הניתוח. הגדר את המצלמה בגובה קבוע (למשל, 50 ס"מ) וזווית (למשל, 30, 50, או 70 מעלות) כדי הדגימה, ולצלם מעל סיבוב מלא (360 מעלות) של פטיפון, במרווחים של 15 ° . הערה: preseעבודת NT, שלושה סטים של תמונות נרכשו מגובה רב מצלמה מרובה וזוויות כדי למפות את הזנים מעל הדגימה כולו 15. טענת את התמונות לתוך תוכנת הניתוח להרכיב 3D האופטי, והמשך לחשב את הזנים. עושה זאת על ידי לחיצה על פונקצית 'האליפסות צרורות מחשוב ", אשר מזהה את נקודות רשת, ואחריו לחיצה על" נקודות ורשת 3D מחשוב' הפונקציה אשר בונה את הרשת. הערה: חישוב הזנים ולדמיין אותו במצב ההערכה. פלט ההפצות להתאמץ כדי לקבוע את זני גבול עבור כל דגימה המבוססת על ISO 12004 16, ואת עלילת תרשימי הגבול להרכיב עבור מהירויות ויוצרים שונות וטמפרטורות טביעה. כיול מודל חומר AA6082 בטמפרטורות שונות 300 כדי 500 מעלות צלזיוס שיעורי זן מן 0.1 כדי 10 שניות -1. הערה: מודל החומר והקבוע שלה עבור AA6082מפורטים הפניה 17. ליישם ולאחד את פונקצית תשואת איזוטרופי Hosford 18, מרצ'יניאק-קוצ'ינסקי (ח"כ) התאוריה 19 ומודל החומר בשלב 1.12 לתוך אלגוריתם אינטגרציה כדי לגבש מודל חיזוי גבול הטביעה. הערה: המודל מתואר התייחסות 11. כיול ולאמת את המודל שפותח על צעד 1.13 באמצעות ובין התוצאות שהושגו בשלב 1.11. לחזות את הגבולות ויוצרים באמצעות המודל אומת 11 משלב 1.14. הערה: איור 1 מציג את תחזיות המודל המתקבל בטמפרטורות שונות, במהירות ויצרו של 250 מ"מ / s, או באופן שקול, שיעור זן של 6.26 s -1. פיתוח 2. מודל אינטראקטיבי חיכוך / תלבש בצע בדיקות כדור-על-דיסק עבור מצופה דגימות (דיסק) כן ניטריד טיטניום (פח) ציפוי על פלדת נושאותדיסק GCr15 באמצעות קשת לקתודה מקרטעת magnetron בתדר באמצע, עם הפרמטרים בתצהיר שניתן ביחס 20. באמצעות מיקרוסקופ אלקטרונים סורק (SEM), להשיג הטופוגרפיה משטח / חתך של המדגם מצופה. מדוד את עובי הציפוי פח באמצעות תמונות SEM ידי השוואת הטופוגרפיה (בהירות חוזה) של חומרי בסיס וציפוי. הערה: הפרוצדורות ניתן למצוא התייחסות 20. השתמש-filometer פרו משטח אור לבן הבין-ferometric להשיג את חספוס פני השטח של המדגם. מניח את המדגם תחת העדשה ולהתאים את מיקרוסקופ כדי להשיג מבנה שטח ברור. להאיר את המדגם ולהתאים את הזוויות של x ו- y הצירים להתבונן רצועות הפרעה ברורות (אשר ניתן לנטר מהמסך). גדר עומק ברוטו בתוכנה ולהתחיל מדידה. באופן אוטומטי לסרוק את השטח מדגם ולחשב את חספוס פני השטח. להעריך את חוזק חסיד of המדגם באמצעות בודק מיקרו שריטות. החלת עומס הגדלה (מקסימום 50 N) ובמרחק מאפס (5 מ"מ מקסימום) על ציפוי הפח. קבע את הַכשָׁלָה עומס קריטי של הציפוי ולקבל את עקומות מיקרו שריטות 20. להעריך קשיות של המדגם באמצעות indenter קשיות. הפעל עומס סטטי של 20 N על המדגם עבור 15 שניות. מדוד אלכסון הרושם שנעשה על ידי indenter, ולאחר מכן לקבל את ערכי הקשיות מן הבוחן. כדור-על-דיסק לבצע בדיקות על tribometer בסביבת הסביבה (טמפרטורה 25 ° C, לחות 30%). השתמש 6 מ"מ קוטר WC-6 הכדור% (מיקרו-קשיות 1,780 HV, כוח שחיקה 1,380 N / ס"מ, מודולוס אלסטיות 71 GPA) כמקבילה כנגד דיסקת מצופה. כוונו את מהירות הזזה ביחס 5 מ"מ / s. הפעל עומס רגיל של 200 נ התחל הערכים המנועים ולהקליט חיכוך באמצעות tribometer. תפסיקו את הבדיקה ב 180, 350 ים, 400 s, ו -450 s, בהתאמה, כדי לנתח את המסלול ללבוש באמצעות oמיקרוסקופ ptical 20. מדוד את הטופוגרפיה של השטח שחוק באמצעות profilometer משטח interferometric אור לבן לאחר בדיקה. חזור על בדיקות (שלב 2.1.6) עם המון נורמלי שונים (300 N, 400 N). קבע את האבולוציה של מקדם החיכוך עד ההתמוטטות של הציפוי הקשה, המאופיינת בעלייה חדה של מקדם החיכוך מגרש את האבולוציה של מקדם חיכוך נגד הזמן אחרי הקלטת ערכי החיכוך בשלב 2.1.6. הערה: האבולוציה של מקדם החיכוך מוצגת התייחסות 20. להעריך את האבולוציה של מקדם החיכוך מבחינת התנהגות ללבוש ואת המנגנונים הקשורים. הערה: האבולוציה של חיכוך מאופיינת לשלושה שלבים שונים: (i) בשלב חיכוך נמוך, (ii) בשלב חיכוך חריש, וכן (iii) פירוט ציפוי במה 20,21. להעריך את בגדי היםטאטעס ב 180 s ידי לקטוע את הבדיקה באופן ידני, ולאחר מכן לנתח את המסלול ללבוש באמצעות מיקרוסקופ אופטי. הערה: שלב זה הוא לחקור את הפסולה ללבוש לשלב חיכוך הנמוך כמתוארת בשלב 2.2.2. חזור על שלב 2.2.3 ב 350 s, 400 s, ו -450 s, בהתאמה. לפתח את המודל חיכוך אינטראקטיבי לאפיין את μ מקדם חיכוך הכולל על ידי שילוב של α חיכוך μ הראשוני עם חיכוך החריש של חלקיקי חומרת μ Pc (כפי שמוצג במשוואה. (1)) 20. (1) מערבב את חיכוך החריש בין הכדור לבין המצע (μ Ps) עם עובי הציפוי המיידי (ח), אשר תמחיש את העלייה החדה מושרת פירוט ציפוי של חיכוך חריש μ Pc (Eq. (2)). הערה: במקרה זה, μ Pc שווה μ Ps כאשר עובי הציפוי שנותר הוא אפס (המציין את הפירוט המלא של הציפוי הקשיח). (2) שם λ 1 ו λ 2 פרמטרי מודל הציגו לייצג את המשמעות הפיסית של התהליך ללבוש. λ 1 מתאר את השפעתו של חלקיקים ללבוש גדולים בפח, ו λ 2 מייצג את עוצמת השפעת חיכוך החריש, המאופיינת במדרון של מקדם החיכוך. השתמש באלגוריתם אינטגרציה המבוססת זמן כדי להשיג את ההתפתחות של עובי הציפוי נותר ומדגים הבלאי המצטבר בתנאי קשר שונים. עדכון עובי הציפוי בכל לולאת חישוב לפי משוואה. (3). (3) <img alt="משוואה 3" src="/files/ftp_upload/53957/53957eq3.jpg"/> כאשר h 0 הוא עובי הציפוי הראשוני הוא השיעור ללבוש תלוי הזמן של הציפוי. שינוי החוק בגדי Archard 22 (Eq. (4)) וליישם אותו במודל הנוכחי. (4) כאשר K הוא מקדם שחיקה, P הוא לחץ הקשר, V היא מהירות זזה, ו- C H הוא הקשיות בשילוב של הציפוי ואת המצע. השתמש במודל של קורסונסקי לחשב את הקשיות המשולבות (Eq. (5)). (5) כאשר H s הוא הקשיות של המצע, α הוא יחס הקשיות בין ציפוי מצע β הוא מקדם ההשפעה של העובי. לייצג את עומס פרמטרים תלוי λ 1 ו- K על ידי כוח lמשוואות aw. (6) (7) שם κ λ1, κ K, Ν λ1 ו Ν K הם קבועים חומר הקשור בהתפתחות של חיכוך 20. התאם את מודל החיכוך האינטרקטיווי את תוצאות הניסוי באמצעות אלגוריתם אינטגרציה שפותח בקבוצה 'המחברים מנת לקבוע את הפרמטרים במודל. 3. סיפורי סימולציה KBC-FE KBC-FE במקרה סימולציה מחקר 1: חיזוי של גבול להרכיב בתנאים בהבלטה חמה צור שם פרויקט סימולציה חדש בתוכנת סימולציה FE. בחר את התהליך כמו 'בול חם טביעה' ואת הסוג הפותר כמו 'PAM-AutoStamp' כאשרשמירת הפרוייקט. ייבא את למות הדלת הפנימית על ידי לחיצה על 'CAD כלים ייבוא' ולאחר מכן 'ייבוא transfe r' הדלת הפנימית קובץ הגיאומטריה 'IGS' לתוך ממשק גרפי תוכנת סימולציה FE. בחר את אסטרטגיית "להרכיב חם 'עבור meshing של כלים. שם האובייקט המיובא כמו'תמותי'. חזור על שלב 3.1.2 ו 'יבוא' האובייקטים של פאנץ' Blankholder, בהתאמה. הקישו על 'בלנק' תחת לשונית 'הגדרת'. לחץ על 'הוסף ריק' ב 'עורך בלנק ", ולהגדיר את" אובייקט חדש' כמו 'בלנק'. לאחר מכן בחר את הסוג כמו 'Surface בלנק'. בחר 'מתאר' עבור סוג ההגדרה ולייבא את b הצורה הריקהלחיצה y על 'ייבוא מקובץ CAD'. גדר 'חידוד' כמו 'רמה שהוטלה' ובחר רמה 1 תחת 'אפשרויות Mesh'. בטל 'meshing האוטומטי' ולהגדיר 'גודל Mesh' עד 4 מ"מ. הגדרת מאפייני החומר 'עורך בלנק'. לחץ על 'טען חומר' תחת לשונית 'חומר'. בחר את 'AA6082' (יחידה: מ"מ · ק"ג · ms · C) חומר כמו תכונות החומר. הגדר את 'בכיוון מתגלגל' ל 'x = 1'. הגדר את 'העובי בלנק' עד 2 מ"מ, ו 'הטמפרטורה הראשונית' ריק ל -490 מעלות צלזיוס. הערה: תכונות החומר ודגם חומר מתוארות ביחס 17. לחץ על 'תהליך </strong> 'תחת' 'הכרטיסיות ובחר "הגדרת סמל +' כדי לטעון מאקרו חדש. דפדף אל 'Hotforming בול ' ובחר 'HF_Validation_DoubleAction_GPa.ksa'. בתיבת הדו-שיח 'התאמה אישית', להפעיל את בלנק, Die, פונץ ', ו Blankholder. תחת הלשונית 'שלבים', להפעיל כוח המשיכה, אחזקות, Stamping, מרווה. כוונו את כל הפרמטרים 'האובייקטים מייחסים' תחת לשונית 'הגדרה' להתכתב עם הגדרת הניסוי בפועל (כוח החזקה ריק = 50 kN, להרכיב מהירות = 250 מ"מ / s, מקדם חיכוך = 0.1, העברת חום מקדם 23 בתור פונקציה של לחץ פער ומגע). לחץ על 'בדוק' סמל לבדוק את הגדרת סימולציה ולהבטיח שלא חלו טעויות ההגדרות הנ"ל. יש ללחוץ על הצלמית 'חישוב' כדי להפעיל את סיניסוח מתואר. הערה: התוכנה רושמת 11 מדינות במהלך הסימולציה במחשב מארח. לאחר השלמת הדמיה, להתבונן בתוצאות הסימולציה בממשק גרפי תוכנת סימולציה FE, והמשך להקליט "תסריט" בשל מעשה ייצוא ערכי גובה, כלומר, המתח העיקרי (קרום), זן קטין (קרום), וטמפרטורה של כל האלמנטים הריקים, עבור מדינת סימולציה שצוינה. לחץ על 'שיא' וערכי קונטור יצוא ידני. לחץ על 'עצור' כדי לעצור את ההקלטה. שמור את התסריט כדי לחזור על אותה הפעולה עבור כל 11 מדינות הסימולציה. לחץ על 'הפעל' סמל לטעון את התסריט, לחץ על 'האם כל "לייצא את ערכי הגובה. הערה: כל קונטור / מדינה בודדת, התוכנה באופן אוטומטי מייצאת הערכים בקבצים 'ASCII' תחת 'major_strain_statenאדמדם ',' minor_strain_statenumber ', ו' temperature_statenumber ', בהתאמה. שמור את כל הקבצים לייצא למחשב ענן. הפעל את 'מודל חיזוי מזמוז' (קוד מודול כלומר, ענן) יחד עם כל הקבצים מיוצאים במחשב הענן. לחזות את ההתפרצות של מזמוז באמצעות יצירת מודל חיזוי גבול במחשב הענן. הערה: זה מודל 11 נותן למשתמשים את האפשרות להפעיל את מודל החיזוי על כמרכיב בודד או כל האלמנטים של ריק. ידני קלט את פרטי הסימולציה / פרמטרים "מודל חיזוי המזמוז '. הזן את מספר המדינות בהדמיה (מדינת 11), שבץ המוחלט של תהליך ההטבעה (157 מ"מ), רוקע מהירות (250 מ"מ / s), זן מגוון של עניין (קריטריון לבחירת האלמנט, למשל, זן> 0.2) ו כל האלמנטים. הערה: straבטווח מגביל את האלמנטים עבורו מזמוזים יכולים להתבצע על ידי קביעת קריטריון אלמנט, למשל, רק את הרכיבים עם זן מהותי אחרון יותר מ -0.2 נבחרו לצורך הערכה נוספת במודול. לאחר השלמת חישוב מודול במחשב ענן, לשמור באופן אוטומטי את כל הנתונים (מזמוז תוצאות חיזוי) לקבצים 'ASCII' מעוצב. טען את המצב הסופי של תוצאות הסימולציה FE. בכרטיסייה 'קווי המתאר', לחץ על 'ייבוא' ולאחר מכן 'ערכים סקלר'. בחר את הקובץ 'ASCII' המתקבל בשלב לעיל. הצג את תוצאות חיזוי מזמוז בתוכנת סימולצית FE. KBC-FE במקרה סימולציה מחקר 2: חיזוי חיי הכלי בתנאי העמסה רב-מחזור צור שם פרויקט סימולציה חדש בתוכנת סימולציה FE. בחר את procESS כמו 'רגילה הטבעה' ואת סוג פותר כמו 'PAM-AutoStamp' בעת שמירת הפרויקט. ייבא את הגיאומטריה למות על ידי לחיצה על 'CAD כלי היבוא' ולאחר מכן 'יבוא העברה' קובץ גיאומטרית צורת U למות 'IGS' לתוך ממשק גרפי תוכנת סימולצית FE. בחר את האסטרטגיה 'אימות' עבור meshing של כלים. שם האובייקט המיובא כמו'תמותי'. חזור על שלב 3.2.2 לייבא את האובייקטים של פאנץ' Blankholder, בהתאמה. הקש על 'בלנק' תחת לשונית 'הגדרה'. 'הוסף ריק' ב 'העורך בלנק', הגדר את 't objec החדש' כמו 'בלנק', ולאחר מכן בחר את הסוג כמו 'Surface בלנק'. בחר "ארבעה פויnts 'עבור סוג ההגדרה ולהגדיר את הגודל הריק 120 × 80 מ"מ 2. גדר 'חידוד' כמו 'שהוטלה רמה': רמת 1 תחת 'אפשרויות Mesh'. בטל 'meshing האוטומטי' ולהגדיר 'גודל Mesh' ל -1.5 מ"מ. הגדרת מאפייני החומר 'עורך בלנק'. הקש על 'טען חומר' תחת לשונית 'החומר'. בחר את 'AA5754-H111' (יחידה: מ"מ · ק"ג · ms · C) חומר כמו תכונות החומר. הגדר את 'בכיוון מתגלגל' ל 'x = 1'. הגדר את 'עובי בלנק' ל -1.5 מ"מ. לחץ על 'תהליך' תחת לשונית 'הגדרה' ובחר בסמל '+' כדי לטעון מאקרו חדש. דפדף אל' בול היתכנות' ובחר 'SingleActioin_GPa.ksa'. בתיבת הדו-שיח 'התאמה אישית', להפעיל את בלנק, Die, פונץ ', ו Blankholder. תחת 'שלבים', להפעיל כוח משיכה, אחזקות, ומחתים. הגדר את כל 'הפרמטרים' בסימולציה להתכתב עם הגדרת הניסוי בפועל (כוחות החזקה ריקה = 5, 20, 50 kN, בהתאמה, מהירות להרכיב = 250 מ"מ / s, מקדם חיכוך = 0.17). 'בדוק' סט-אפ סימולציה ולהבטיח שלא חלו טעויות ההגדרות הנ"ל. לחץ על סמל 'חישוב' ולהתחיל את 'חישוב' עבור סימולציה כיפוף 11 מדינות צורת U במחשב המארח. לאחר השלים ההדמיה, יצוא לתאם "נתונים ונתונים 'לחץ קשר' אוטומטי עבור יצירת העבודהכלים (אגרוף, למות ומחזיק ריק) כקבצים 'ASCII' (לפי שלבים 3.1.11 ו 3.1.12). שמור את כל הקבצים לייצא למחשב ענן. הפעל את 'מודול חיזוי חיי הכלי' יחד עם כל הקבצים מיוצאים במחשב הענן. ידני קלט להרכיב פרמטרי 'מודול כלי חי החיזוי ". הזן את הפרמטרים הבאים: מספר המדינות (מדינה 11), שבץ כולל (70 מ"מ), רוקע מהירות (250 מ"מ / sec) ואת מקדם חיכוך ראשוני (0.17). בחר בכלי (אגרוף, למות, או בעל ריק), ולאחר מכן להתחיל את החישוב עבור אלמנט אחד או את כל האלמנטים. לאחר השלמת חישוב מודול במחשב ענן, לשמור באופן אוטומטי את כל הנתונים (כולל עובי הציפוי שנותר מיידי מקדם חיכוך) לקבצים 'ASCII' מעוצב. לטעון ולהציג את עובי הציפוי הנותרים frictiעל מקדם עבור האלמנטים הרלוונטיים תוכנת סימולציה FE (לפי שלב 3.1.17).

Representative Results

סימולציה KBC-FE לחיזוי מזמוז בתהליך ביול חם, השימוש ריק צורה אופטימיזציה לא רק לחסוך עלות החומר, אלא גם לעזור כדי לצמצם את נוכחותם של פגמים, כגון מזמוז, פיצוח, וקמטים. הצורה הריקה הראשונית משפיעה על זרימת החומר באופן משמעותית במהלך גיבוש, ולכן עיצוב הגיוני של הצורה הריקה הוא קריטי להצלחת התהליך בהבלטה החם ואיכות המוצרים הסופיים. כדי להפחית את המאמצים של ניסויי הניסוי וטעייה כדי לקבוע הגיאומטריה הריקה האופטימלית, סימולצית KBC-FE הוכחה להיות שיטה יעילה ויעילה במיוחד למזעור באזורים עם מזמוז. באמצעות טכניקה זו, כל סימולציה לוקחת כ 2 שעות, ואילו חישוב מודול ענן במקביל לחיזוי מזמוז יושלם בתוך 4 שעות. <p class="jove_content" fo:keep-together.within-page = "1"> איור 4 מראה את האבולוציה של הצורה הריקה המשמשת ביול החם, דוגמא הרכיב פנימי דלת רכב. הצורה הריקה הראשונית, אמצה מתהליך הביול קר קונבנציונלי, שמשה לראשונה בסימולצית KBC-FE. תוצאות הניסוי באיור 4 (א) מראים כי כישלון גדול (פיצוח או מזמוז) באזורים גלויים לאחר בהבלטה חמה. לאחר איטרציה אחת של אופטימיזציה הצורה הריקה, ניתן לראות באיור 4 (ב) כי פנל מוצלח כמעט מלא נוצר עם מזמוזים הרבה פחות, בהשוואה לשימוש בצורה הריקה הראשונית. ניתן לראות כי יש עדיין אינדיקציה מזמוז על הכיסים, בפינה הימנית העליונה ופינות שמאל של החלונית. לאחר אופטימיזציה נוספת באיור 4 (ג), בצורה הריקה אופטימיזציה לבסוף הושגה ללא מזמוזים גלויים בלוח. הצורה הריקה אופטימיזציה נקבעה על ידי סימולצית KBC-FE אומתה באופן ניסיוני דרך בהבלטה חמהניסויים שנערכו על קו ייצור אוטומטי לחלוטין המוצע על ידי יצרן מערכת ייצור. סימולצית KBC-FE עבור חיזוי כלי חיים סימולציות FE קונבנציונליות של תהליכי מתכת יוצר מבוצעות עבור מחזור אחד. עם זאת, בסביבת ייצור, מחזוריים להרכיב מספר מבוצעים על כלי מסוים, שבו הוא נמצא כי עלייה במספר המחזורים להרכיב התוצאה היא וריאציה גדלה בין המרכיבים יוצר. וריאציה זו במהלך טעינה לכלי רב-מחזור היא התוצאה של שינוי טופוגרפית משטח. לדוגמא, הטעינה רב-המחזור של יצירת כלים עם ציפוי פונקציונלי תוביל להפחתת עובי ציפוי עקב הבלאי. יתר על כן, הפירוט של הציפוי גם יושפע להרכיב פרמטרים, כגון העומס / לחץ, ויצר מהירויות, וכו טכניקת KBC-FE מאפשרתסימולציה של מתכת ויוצרי תהליכים בתנאי העמסה רב-מחזור, שהוא חיוני עבור חיזוי חי השירות של יצירת כלים עם ציפוי תפקודי מתקדם. כדי לחקור את ההשפעות של כוח החזקה ריק על חיי הכלי, ערכי כוח ההחזקה ריקים של 5, 20, ו -50 kN נבדקו במהירות ויוצרים קבועה של 250 מ"מ / s. איור 5 מציג את התפלגות עובי ציפוי כולים נותרת עם כוחות החזקה ריקים שונים אחרי 300 מחזורי טביעה. זה עולה בבירור כי עובי הציפוי נותר פוחת עם גידול הכח מחזיק הריק. איור 6 מציג את ציפוי לחץ ואת הזמן שנותר הפצת עובי עם כוחות החזקה ריקים של 5, 20, ו -50 kN, בהתאמה, לאורך מרחק curvilinear של הקובייה לאחר 300 מחזורי טביעה. מאז באזור AB מייצג המזון למותבאזור ראנס במהלך תהליך כיפוף צורת U, הלחץ והמרחק ללבוש יחסית באזור זה, היו הרבה יותר גבוהים מאשר באזורים אחרים של הקובייה. כתוצאה מכך, שחיקה של הציפוי בעיקר התרחשה בתחום זה. ישנם שני ערכי השיא של הפחתת עובי הציפוי ב 20 kN 50 kN שמתאימים לשני פסגות בלחץ. בינתיים, עובי הציפוי נותר פוחת עם הגידול של כוח החזקה ריק. עובי הציפוי נותר הנמוך ביותר עם כוחות החזקה ריקים של 5, 20, ו -50 kN, היה 0.905, 0.570, ו 0.403 מיקרון, בהתאמה, שבו עובי הציפוי הראשוני היה 2.1 מיקרון. איור 1: השוואה בין זני גבול הניסיוניים וחזה להרכיב בטמפרטורות שונות. זני הגבול להרכיב להגדיל עקב עלייה בטמפרטורה, במהירות קבועה של 250 מ"מ/ s, או באופן שקול, שיעור זן של 6.26 s -1. אנא לחץ כאן כדי לצפות בגרסה גדולה יותר של דמות זו. איור 2: תרשים סכמטי לסימולצית FE ענן מבוסס ידע של תהליך גיבוש מתכת. תוכנת סימולציה FE מסחרי, משמש להרצת הסימולציה ולייצא את התוצאות הנדרשות עבור מודולים בודדים. המודולים, למשל, formability, העברת חום, חוזק יוצרי פוסט (מיקרו), חיזוי חיי הכלי, עיצוב כלי, וכו ', לעבוד במקביל ובאופן עצמאי בענן, ולכן המאפשר שילוב של חיתוך ידע קצה ממקורות מרובים לתוך סימולציות FE . אנא CLאיכס כאן כדי לצפות בגרסה גדולה יותר של דמות זו. איור 3: גיאומטריה של חתיכת העבודה וכלי צורת U כיפוף סימולציה. הכלים, כלומר, אגרוף, בעל ריק ולמות, הם המודל באמצעות אלמנטים נוקשה. אלמנטי מעטפת משמשים חתיכת עבודת אלמנטים (ריק). אנא לחץ כאן כדי לצפות בגרסה גדולה יותר של דמות זו. איור 4: אבולוציה של צורה ריקה עבור ביול חם של פנל דלת פנימי (המוצג סימולצית FE). שמאל: הנתונים במסגרות ירוקות מייצגים צורות ריקות בכל שלב אופטימיזציה, ואלה באדומיםמסגרות מתאימות בצורה הריקה לפני הייעול שלה. מימין: מזמוז תוצאות חיזוי בכל שלב אופטימיזציה. (א) תוצאות ראשוניות עם כישלון גדול (פיצוח / מזמוז יוצגו בצבע אדום), (ב) צמצום כישלון עם מזמוזים כמה לאחר השלב הראשון של אופטימיזציה, (ג) צורה ריקה אופטימיזציה סופי ללא מזמוזים גלויים. אנא לחץ כאן כדי לצפות בגרסה גדולה יותר של דמות זו. איור 5: התפלגות עובי ציפוי הנותרת (המוצג סימולצית FE) עם כוחות החזקה ריקים: (א) 5 kN, (ב) 20 kN, ו (ג) 50 kN, אחרי 300 מחזורים להרכיב במהירות ביול קבועה של 250 מ"מ / s. אנאלחץ כאן כדי לצפות בגרסה גדולה יותר של דמות זו. איור 6: חיזוי של מגע עם שימוש בלחץ ועובי ציפוי הנותרים עם כוחות החזקה ריקים: (א) 5 kN, (ב) 20 kN, ו (ג) 50 kN, לאורך מרחק curvilinear של הקובייה במהירות ביול קבועה של 250 מ"מ / s. אנא לחץ כאן כדי לצפות בגרסה גדולה יותר של דמות זו.

Discussion

טכניקת סימולצית KBC-FE מאפשרת סימולציות מתקדמות להתנהל מחוץ לאתר באמצעות מודולים ייעודיים. זה יכול לפעול מודולים פונקציונליים על סביבת ענן, המקשרים את בלוטות מתחום התמחות אחרת, על מנת להבטיח כי סימולציות תהליך מתנהלות בצורה מדויקת ככל האפשר. ההיבטים הקריטיים סימולצית KBC-FE עשויים להיות כרוכים עצמאיים של קודי FE, יעילות של החישוב, ודיוק של מודולים הפונקציונליים. המימוש של כל פונקציה מתקדמת מודול היה להסתמך על פיתוח מודל חדש ו / או טכניקה ניסויית רומן. לדוגמא, מודול הגבול להרכיב מפותח המבוססת על מודל החיזוי ויוצר גבול המאוחד החדש 11, ואת תחזית מודול חי חיכוך כלי פותחתה כיום על ידי יישום מודל החיכוך האינטרקטיווי 20. טכניקת סימולצית KBC-FE מציעה גם את הפונקציה של חישוב סלקטיבי, כלומר, רק אלמנטים למלא את הבחירהקריטריונים נבחרו לצורך הערכה נוספת של מודולים בודדים. למשל, מודול חיזוי חיי כלי בוחר באופן אוטומטי את האלמנטים שלשם הציפוי הקשיח נוטה התמוטטות, לפי דירוג שיעור הבלאי של כל האלמנטים ב -1 להרכיב מחזור, ולכן בדרך כלל פחות מ -1% של האלמנטים ייבחרו נוסף ערכות חיי כלי בתנאי העמסה רב-מחזור. במחקר הנוכחי, ניבוי חיי כלי לאחר 300 מחזורי טביעה יכול להסתיים בתוך 5 דקות.

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

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

כטכנולוגיה בעתיד, סימולצית KBC-FE תסתמך על פיתוח חבילות תוכנת סימולצית FE מבוססת אינטרנט ייעודי ויציב, אשר ידרוש רווחי מאוד, אבל לגמרי מודל עסקי שונה שיוקם על ידי מפתחי תוכנה. בנוסף, רשת פנימית ייעודית צריכה להיות בנויה בתוך המפלגות השיתופיות על מנת להבטיח אבטחת מידע והאמינות המלאה של המערכת התעשייתית. </p>

Disclosures

The authors have nothing to disclose.

Acknowledgements

The financial support from Innovate UK, Ultra-light Car Bodies (UlCab, reference 101568) and Make it lighter, with less (LightBlank, reference 131818) are gratefully acknowledged. The research leading to these results has received funding from the European Union’s Seventh Framework Program (FP7/2007-2013) under grant agreement No. 604240, project title ‘An industrial system enabling the use of a patented, lab-proven materials processing technology for Low Cost forming of Lightweight structures for transportation industries (LoCoLite)’. Significant support was also received from the AVIC Centre for Structural Design and Manufacture at Imperial College London, which is funded by Aviation Industry Corporation of China (AVIC).

Materials

AA6082-T6 AMAG Material
AA5754-H111 AMAG Material
 1000 kN high-speed press ESH Forming press
ARGUS GOM Optical forming analysis
PAM-STAMP 2015 ESI FE simulation software
Matlab  MathWorks Numerical calculation software
Gleeble 3800 DSI Uniaxial tensile test
High Temperature Tribometer (THT) Anton Paar Friction property test
NewViewTM 7100 ZYGO Surface profilometer
 Magnetron sputtering equipment Coating deposition
Microhardness tester Wolpert Wilson Instruments
Nano-hardness indenter  MTS

References

  1. Miller, W. S., et al. Recent development in aluminium alloys for the automotive industry. Mater. Sci. Eng. A. 280 (1), 37-49 (2000).
  2. Bolt, P. J., Lamboo, N. A. P. M., Rozier, P. J. C. M. Feasibility of warm drawing of aluminium products. J. Mater. Process. Tech. 115 (1), 118-121 (2001).
  3. Li, D., Ghosh, A., et al. Effects of temperature and blank holding force on biaxial forming behavior of aluminum sheet alloys. J. Mater. Eng. Perform. 13 (3), 348-360 (2004).
  4. Toros, S., Ozturk, F., Kacar, I. Review of warm forming of aluminum-magnesium alloys. J. Mater. Process. Tech. 207 (1-3), 1-12 (2008).
  5. Wang, L., Strangwood, M., Balint, D., Lin, J., Dean, T. A. Formability and failure mechanisms of AA2024 under hot forming conditions. Mater. Sci. Eng. A. 528 (6), 2648-2656 (2011).
  6. Wang, L., et al. . TTP2013 Tools and Technologies for Processing Ultra High Strength Materials. , (2013).
  7. El Fakir, O., et al. Numerical study of the solution heat treatment, forming, and in-die quenching (HFQ) process on AA5754. Int. J. Mach. Tool. Manu. 87, 39-48 (2014).
  8. Raugei, M., El Fakir, O., Wang, L., Lin, J., Morrey, D. Life cycle assessment of the potential environmental benefits of a novel hot forming process in automotive manufacturing. J. Clean. Prod. 83, 80-86 (2014).
  9. Liu, J., Gao, H., Fakir, O. E., Wang, L., Lin, J. HFQ forming of AA6082 tailor welded blanks. MATEC Web of Conferences. 21 (05006), (2015).
  10. Karbasian, H., Tekkaya, A. E. A review on hot stamping. J. Mater. Process. Tech. 210 (15), 2103-2118 (2010).
  11. El Fakir, O., Wang, L., Balint, D., Dear, J. P., Lin, J. Predicting Effect of Temperature Strain Rate and Strain Path Changes on Forming Limit of Lightweight Sheet Metal Alloys. Procedia Eng. 81, 736-741 (2014).
  12. Shi, Z., et al. . the 3rd International Conference on New Forming Technology. , 100-104 (2012).
  13. Mohamed, M. S., Foster, A. D., Lin, J., Balint, D. S., Dean, T. A. Investigation of deformation and failure features in hot stamping of AA6082: Experimentation and modelling. Int. J. Mach. Tool. Manu. 53 (1), 27-38 (2012).
  14. Hosford, W. F. Comments on anisotropic yield criteria. Int. J. Mech. Sci. 27 (7), 423-427 (1985).
  15. Marciniak, Z., Kuczyński, K. Limit strains in the processes of stretch-forming sheet metal. Int. J. Mech. Sci. 9 (9), 609-620 (1967).
  16. Ma, G., Wang, L., Gao, H., Zhang, J., Reddyhoff, T. The friction coefficient evolution of a TiN coated contact during sliding wear. Appl. Surf. Sci. 345, 109-115 (2015).
  17. Põdra, P., Andersson, S. Simulating sliding wear with finite element method. Tribol. Int. 32 (2), 71-81 (1999).
  18. Archard, J. F. Contact and Rubbing of Flat Surfaces. J. Appl. Phys. 24 (8), 981-988 (1953).
  19. Liu, X., et al. Determination of the interfacial heat transfer coefficient in the hot stamping of AA7075. MATEC Web of Conferences. 21 (05003), (2015).

Play Video

Cite This Article
Zhou, D., Yuan, X., Gao, H., Wang, A., Liu, J., El Fakir, O., Politis, D. J., Wang, L., Lin, J. Knowledge Based Cloud FE Simulation of Sheet Metal Forming Processes. J. Vis. Exp. (118), e53957, doi:10.3791/53957 (2016).

View Video