Prediction of plastic deformation zone width in SAW-welded structural steel butt joints using numerical simulations and analytical models

ملخص تقني: التنبؤ بعرض منطقة التشوه اللدن في وصلات التراكب الفولاذية الملحومة باستخدام القوس المغمور (SAW)

بيان المشكلة
يعد التنبؤ الدقيق بالتشوهات الناجمة عن اللحام أمراً بالغ الأهمية لتصميم وتصنيع الهياكل الفولاذية الضخمة، لا سيما في التطبيقات البحرية والبحرية المفتوحة. ويُعتبر التحديد الدقيق لعرض منطقة التشوه اللدن ($b_{pl}$)، المعروفة أيضاً باسم منطقة الانفعال المتأصل، شرطاً أساسياً مسبقاً لطرق التنبؤ المبسطة بالتشوه — وتحديداً تلك القائمة على نظرية الانفعال المتأصل والأحمال المكافئة. وبينما توفر تحليلات العناصر المحدودة الحرارية-المرنة-اللدنة (TEP-FEM) التقليدية دقة عالية، إلا أن تكلفتها الحسابية غالباً ما تكون باهظة للهياكل المعقدة. وفي المقابل، تعتمد النماذج التحليلية الموجودة على افتراضات مبسطة قد لا تستطيع استيعاب التفاعلات الحرارية-الميكانيكية المعقدة في اللحام بالقوس المغمور (SAW) للفولاذ الإنشائي بشكل كامل. وهناك حاجة إلى علاقات تنبؤية قوية وفعالة حسابياً تأخذ في الاعتبار معايير تكنولوجية وتصميمية محددة لتقدير $b_{pl}$ بدقة عبر نطاقات مختلفة من سماكة الألواح.

المنهجية
استخدمت الدراسة نهجاً هجيناً يجمع بين المحاكاة العددية، والتصميم الإحصائي للتجارب (DOE)، وتحليل الانحدار لتطوير نماذج تنبؤية لـ $b_{pl}$.

1. النمذجة العددية: أُجريت تحليلات حرارية عابرة باستخدام برنامج ANSYS مع نموذج مصدر حرارة غاوسي سطحي. غطت عمليات المحاكاة فئتين متميزتين من سماكة الألواح بناءً على السلوك الحراري:

   • الألواح الرقيقة ($b_{pl-tp}$): 6–26 مم، وهي تمثل هياكل السفن البحرية.
   • الألواح السميكة ($b_{pl-hp}$): 30–50 مم، وهي تمثل المنشآت البحرية المفتوحة.
     استخدم النموذج شبكة غير منتظمة مع تكثيف محلي بالقرب من اللحام لالتقاط تدرجات درجة الحرارة بدقة. تم تحديد حدود المنطقة اللدنة بـ "درجة حرارة الخضوع الأولى" ($T_{Y1}$)، وهي عتبة يبدأ عندها التشوه غير العكسي، والتي حُسبت بناءً على قوة خضوع المادة وخصائص التمدد الحراري.

2. تصميم التجارب (DOE): تم تطبيق تصميم Box–Behnken ثلاثي المستويات (BBD) بـ 15 تجربة لكل مجموعة ألواح.

   • المتغيرات المستقلة للألواح الرقيقة: مدخلات الحرارة للحام ($q_l$)، قوة خضوع المعدن الأساسي ($\sigma_{YB}$)، وسماكة اللوح ($h$).
   • المتغيرات المستقلة للألواح السميكة: مدخلات الحرارة للحام ($q_l$)، درجة حرارة البداية/التسخين المسبق ($T_0$)، وسماكة اللوح ($h$).
   • المتغير التابع: عرض منطقة التشوه اللدن ($b_{pl}$).

3. التحقق من النموذج: تم التحقق من النموذج العددي مقابل نتائج مرجعية من Grams [17] لمواد ASTM A36 وAISI 316L، حيث أظهر توافقاً بنسبة 8–10%.

4. تحليل الانحدار: أنتج التحليل الإحصائي باستخدام برنامج STATISTICA معادلات تنبؤية متعددة الحدود. وتم تقييم أهمية المتغيرات وحدود التفاعل عند مستوى ثقة 95%.

5. التحليل المقارن: تمت مقارنة نماذج الانحدار المطورة مع الصيغ التحليلية لـ Ueda وآخرون [4, 5] ومع نتائج المحاكاة العددية لتقييم دقة التنبؤ.

المساهمات والنتائج الرئيسية

• تطوير المعادلات التنبؤية: استخلصت الدراسة نموذجين انحدار متميزين بصدقية إحصائية عالية ($R^2_{adj} \approx 0.97$ للألواح الرقيقة و $\approx 0.98$ للسميكة).

  • نموذج اللوح الرقيق ($b_{pl-tp}$): حدد مدخلات الحرارة للحام ($q_l$)، وسماكة اللوح ($h$)، وقوة الخضوع ($\sigma_{YB}$) كعوامل مهمة، بما في ذلك حدود التفاعل ($q_l \cdot \sigma_{YB}$، $q_l \cdot h$، $\sigma_{YB} \cdot h$) وحد تربيعي للسماكة ($h^2$).
  • نموذج اللوح السميك ($b_{pl-hp}$): حدد مدخلات الحرارة للحام ($q_l$) وسماكة اللوح ($h$) كعوامل مهيمنة، مع حد تفاعل مهم ($q_l \cdot h$) وحد تربيعي لمدخلات الحرارة ($q_l^2$). ومن الجدير بالذكر أن درجة حرارة التسخين المسبق ($T_0$) وُجد أنها غير مهمة إحصائياً ضمن النطاق المختبر (100–200 درجة مئوية)، وتم استبعادها من المعادلة النهائية.

• العوامل المهيمنة: أكد التحليل أن مدخلات الحرارة للحام هي المعلم الأكثر تأثيراً في التحكم في عرض المنطقة اللدنة في كلتا المجموعتين، يليه سماكة اللوح. وكان التفاعل بين مدخلات الحرارة والسماكة حاسماً في كلا النموذجين.

• الدقة والتحقق:

  • أظهرت نماذج الانحدار توافقاً ممتازاً مع بيانات المحاكاة العددية. بالنسبة للألواح الرقيقة، كانت 12 حالة من أصل 15 ضمن خطأ تنبؤ يتراوح بين 0–5 مم؛ وبالنسبة للألواح السميكة، كانت جميع الحالات الـ 15 ضمن هذا النطاق.
  • عند المقارنة بالصيغ التحليلية لـ Ueda وآخرون، أظهرت نماذج الانحدار المطورة أخطاء تنبؤ متوسطة أقل بكثير. أظهرت نماذج Ueda تباينات أكبر، خاصة للألواح الرقيقة تحت ظروف مدخلات حرارة عالية (أخطاء تصل إلى 32 مم في بعض الحالات)، بينما حافظت النماذج المقترحة على أخطاء تقل عموماً عن 10 مم.

• انتقال النظام الحراري: بررت الدراسة فيزيائياً تصنيف الألواح إلى فئات رقيقة (6–26 مم) وسميكة (30–50 مم) من خلال ملاحظة انتقال متميز في توزيع درجات الحرارة وآليات تدفق الحرارة بين سمك 26 مم و30 مم.

الأهمية والادعاءات
يزعم البحث أن المنهجية المقترحة توفر بديلاً موثوقاً وفعالاً حسابياً لمحاكاة العناصر المحدودة (TEP-FEM) الكاملة لتقدير عرض منطقة التشوه اللدن في وصلات التراكب الفولاذية الملحومة بـ SAW. ومن خلال دمج تصميم التجارب (DOE) مع المحاكاة العددية، أثبت المؤلفون ضرورة وجود نماذج تنبؤية منفصلة للألواح الرقيقة والسميكة بسبب اختلاف السلوكيات الحرارية.

ويؤكد المؤلفون أن هذه العلاقات التنبؤية ذات قيمة خاصة لـ:

1. التنبؤ المبسط بالتشوه: توفير معلمات مدخلات دقيقة (تحديداً عرض منطقة الانفعال المتأصل) لطرق التنبؤ بالتشوه القائمة على الأحمال المكافئة والانفعال المتأصل، والتي تُستخدم على نطاق واسع في التصميم الأولي للهياكل البحرية والمحيطية.
2. التطبيق الهندسي: تقديم أداة توازن بين الدقة والكفاءة، بحيث تكون مناسبة لمراحل التصميم والتصنيع حيث يلزم التقييم السريع.
3. تحسين العمليات: تمكين المهندسين من فهم حساسية عرض المنطقة اللدنة للمعلمات الرئيسية مثل مدخلات الحرارة وسماكة اللوح دون اللجوء إلى عمليات محاكاة مكلفة.

وتخلص الدراسة إلى أنه بينما تعد النماذج قوية ضمن النطاق التجريبي الذي تم استقصاؤه، يجب تجنب الاستقراء خارج النطاقات المحددة لمدخلات الحرارة وسماكة اللوح للحفاظ على دقة التنبؤ. ولا يدعي العمل استبدال FEM التفصيلي للتصميم النهائي، بل يهدف ليكون تقريباً موثقاً ومبنياً على الفيزياء للمعلم الوسيط الحرج وهو عرض المنطقة اللدنة.