Development of Design-Oriented Bond Strength Equations for FRP Bars Embedded in Concrete Exposed to Elevated Temperatures Using Explainable Machine Learning
تُطوّر هذه الدراسة نموذج تعلم آلي هجين (WOA–ETR) للتنبؤ بدقة بقوة تماسك قضبان البوليمر المعزز بالألياف (FRP) في الخرسانة تحت درجات حرارة مرتفعة، وتستفيد من تحليل SHAP لتحديد العوامل المؤثرة الرئيسية، وتستخلص معادلات تصميمية قوية وقابلة للتفسير الفيزيائي بناءً على هذه الرؤى.
تُعد الخرسانة العمود الفقري للبنية التحتية الحديثة، لكن قضبان الفولاذ المخفية بداخلها يمكن أن تصدأ وتضعف بمرور الوقت، خاصة في البيئات القاسية أو المالحة أو الرطبة. ولحل هذه المشكلة، اتجه المهندسون إلى قضبان البوليمر المعززة بالألياف، وهي نوع من التعزيز المصنوع من ألياف قوية محبوسة داخل راتينج بلاستيكي. هذه القضبان لا تصدأ، وهي أخف وزناً من الفولاذ، ولا توصل الكهرباء، مما يجعلها مثالية للجسور والمباني في الظروف المسببة للتآكل. ومع ذلك، فإن هذه القضبان البلاستيكية لها نقطة ضعف حرجة: فهي حساسة للحرارة. فعند تعرضها لدرجات حرارة عالية، يمكن للراتينج الذي يربط الألياف ببعضها أن يلين ويتدهور، تماماً كما يفقد شمع الشمعة شكله تحت أشعة الشمس. هذا الليونة تضعف التماسك الحيوي، أو "الارتباط"، بين القضيب والخرسانة المحيطة به. وإذا فشل هذا الارتباط، يمكن للهيكل أن يفقد قوته، ومع ذلك، كان التنبؤ بدقة بمقدار قوة التماسك المتبقية بعد التعرض لحريق أو حرارة لغزاً صعباً للمهندسين.
لقد وضع فريق من الباحثين حداً لهذا اللغز من خلال الجمع بين البيانات الهندسية التقليدية والتعلم الحاسوبي المتقدم. قاموا بجمع مجموعة ضخمة مكونة من 132 تجربة واقعية حيث تم سحب هذه القضبان البกว่า من كتل خرسانية بعد تسخينها إلى درجات حرارة متفاوتة، تراوحت من درجة حرارة الغرفة حتى 600 درجة مئوية. كان الهدف هو إيجاد طريقة موثوقة للتنبؤ بمدى قوة التماسك في ظل هذه الظروف. وبدلاً من الاعتماد على معادلات بسيطة غالباً ما تغفل عن تعقيد الموقف، استخدم الباحثون نموذجاً حاسوبياً متطوراً. لقد دربوا هذا النموذج باستخدام تقنية تسمى "انحدار الأشجار الإضافية" (extra trees regression)، والتي تبني العديد من أشجار اتخاذ القرار لإيجاد الأنماط في البيانات، ثم قاموا بتنقيحه باستخدام "خوارزمية تحسين الحوت" (whale optimization algorithm). تعمل هذه الخطوة الثانية كمحرك بحث ذكي، حيث تقوم باستمرار بتعديل إعدادات النموذج لإيجاد الطريقة الأكثر دقة للتنبؤ بالنتائج، محاكيةً بذلك الطريقة التي تصطاد بها حيتان الأحدب في مجموعات منسقة.
أثبت النموذج الحاسوبي دقة ملحوظة، حيث تنبأ بشكل صحيح بقوة التماسك في 94.2 بالمائة من حالات الاختبار. والأهم من ذلك، لم يرد الباحثون مجرد "صندوق أسود" يعطي إجابات دون تفسير؛ بل استخدموا أدوات خاصة للنظر داخل النموذج وفهم أي العوامل كانت الأكثر أهمية. واكتشفوا أن درجة الحرارة التي تعرض لها القضيب كانت العامل الأهم على الإطلاق، يليه عن قرب ملمس سطح القضيب وكمية الخرسانة التي تغطي القضيب. وبينما تؤكد الحكمة الهندسية التقليدية غالباً على قوة الخرسانة نفسها، أظهرت هذه الدراسة أن هذه القضبان تحديداً في الحرارة العالية، كانت درجة الحرارة وملمس سطح القضيب أكثر حسماً. كما كشف النموذج أن طول القضيب المغمر في الخرسانة يغير كيفية سلوك التماسك بطرق مختلفة، مما يشير إلى أن معادلة واحدة لا يمكنها وصف كل حالة.
ولجعل هذه النتائج مفيدة للتصميم في العالم الحقيقي، ترجم الباحثون رؤى الحاسوب إلى معادلات واضحة وعملية. لقد أدركوا أنه نظراً لأن التماسك يسلك سلوكاً مختلفاً اعتماداً على طول الجزء المدفون من القضيب في الخرسانة، فقد احتاجوا إلى تقسيم البيانات إلى ثلاث مجموعات: أطوال انغمار قصيرة، ومتوسطة، وطويلة. وبالنسبة لأقصر القضبان، قاموا بتقسيم البيانات بشكل أكبر بناءً على ملمس سطح القضيب، مثل ما إذا كان مغطى بالرمل أو محززاً. ومن خلال إنشاء معادلات منفصلة لكل سيناريو من هذه السيناريوهات المحددة، أنتجوا صيغاً كانت غير منحازة تقريباً، مما يعني أن القيم المتوقعة كانت مطابقة تماماً للنتائج التجريبية الفعلية. توفر هذه المعادلات الجديدة للمهندسين أداة شفافة وموثوقة لتصميم هياكل أكثر أماناً يمكنها تحمل درجات الحرارة العالية، مما يسد الفجوة بين تعلم الآلة المعقد وسلامة البناء العملية.
ملخص تقني: تطوير معادلات تصميمية لقوة التماسك لأسياخ الـ FRP المغروسة في الخرسانة المعرضة لدرجات حرارة مرتفعة باستخدام التعلم الآلي القابل للتفسير
بيان المشكلة توفر أسياخ البوليمر المعزز بالألياف (FRP) مزايا متانة كبيرة مقارنة بتسليح الفولاذ في البيئات القاسية؛ ومع ذلك، فإن تطبيقها محدود بسبب الحساسية الحرارية لمصفوفة الراتنج البوليمري العضوي (عادةً الإيبوكسي أو فينيل إستر). فعند التعرض لدرجات حرارة مرتفعة، يخضع الراتنج لعملية انتقال زجاجي وتحلل حراري، مما يضعف بشكل كبير التماسك بين سيخ الـ FRP والخرسانة المحيطة به. وبينما تحاول العلاقات التجريبية الحالية وصف سلوك التماسك المعتمد على درجة الحرارة، إلا أنها غالبًا ما تعاني من محدودية التطبيق، وانخفاض دقة التنبؤ، وعدم القدرة على استيعاب التفاعلات المعقدة وغير الخطية بين المعاملات المؤثرة مثل درجة حرارة التعرض، والبروفيل السطحي، وطول الانغراس، وخصائص الخرسانة. علاوة على ذلك، فإن نماذج التعلم الآلي (ML) الحالية غالبًا ما تعمل كـ "صناديق سوداء"، وتفتقر إلى التفسير الفيزيائي المطلوب للدمج المباشر في أكواد التصميم الهندسي.
المنهجية تقترح هذه الدراسة إطار عمل متكامل يجمع بين التعلم الآلي المتقدم والذكاء الاصطناعي القابل للتفسير (XAI) لتطوير معادلات تصميمية قابلة للتفسير فيزيائيًا. تسير المنهجية عبر ثلاث مراحل رئيسية:
تجميع وقبل معالجة قاعدة البيانات: تم تجميع قاعدة بيانات شاملة مكونة من 132 عينة من اختبارات السحب (pull-out) من الأدبيات الموجودة. تتضمن مجموعة البيانات متغيرات مثل مقاومة الضغط للخرسانة (fc′)، وقطر سيخ الـ FRP (db)، ونسبة الغطاء إلى قطر السيخ (c/db)، ونسبة طول الانغراس إلى قطر السيخ (le/db)، ونوع سيخ الـ FRP (GFRP، BFRP)، والبروفيل السطحي (مغطى بالرمل، مضلع، ملفوف بالألياف، إلخ)، ودرجة حرارة التعرض (T). تم ترميز المتغيرات الفئوية، وتم تطبيع البيانات (normalization).
نمذجة هجينة للتعلم الآلي: تم تطوير نموذج هجين يجمع بين خوارزمية تحسين الحوت وريجرسون الأشجار الإضافية (WOA–ETR). استُخدمت خوارزمية تحسين الحوت (WOA)، وهي خوارزمية استدلالية (metaheuristic)، لضبط المعلمات الفائقة (hyperparameters) لنموذج ريجرسون الأشجار الإضافية (ETR) تلقائيًا. تم تدريب النموذج والتحقق من صحته باستخدام 10 جولات مستقلة مع تحقق تبادلي خماسي الطيات (5-fold cross-validation) لضمان المتانة. تم تقييم الأداء باستخدام R2، وRMSE، وMAE، وMAPE.
التفسير واستخراج المعادلات: للتغلب على طبيعة "الصندوق الأسود" للتعلم الآلي، تم استخدام تفسيرات شابلّي المضافة (SHAP) ومخططات الاعتماد الجزئي (PDP) لتحديد أهمية الميزات وتصور التأثيرات الهامشية للمتغيرات المدخلة. بناءً على هذه الرؤى، تم تصنيف قاعدة البيانات بشكل منهجي إلى مجموعات فرعية وفقًا لأنظمة طول الانغراس والبروفيلات السطحية. بعد ذلك، تم تطبيق ريجرسون المربعات الصغرى غير الخطي على هذه المجموعات الفرعية لاشتقاق معادلات تنبؤية مبسطة وقابلة للتفسير فيزيائيًا.
المساهمات الرئيسية
نموذج WOA–ETR الهجين: تطوير نموذج هجين جديد يتفوق على خوارزميات التعلم الآلي التقليدية (بما في ذلك RF، وBR، وGB، وXGBoost، وETR المستقل) في التنبؤ بقوة التماسك تحت درجات حرارة مرتفعة.
إطار الذكاء الاصطناعي القابل للتفسير: تطبيق تحليلات SHAP وPDP على مشكلة قوة التماسك، مما وفر تسلسلًا هرميًا واضحًا للعوامل المؤثرة وكشف عن العلاقات غير الخطية التي تغفلها الصيغ التجريبية التقليدية.
معادلات موجهة للتصميم: اشتقاق سلسلة من المعادلات التنبؤية المحددة بناءً على التصنيف المدفوع بالتعلم الآلي. وبخلاف المعادلة العالمية الواحدة، تأخذ هذه المعادلات في الاعتبار آليات نقل التماسك المتميزة في أنظمة الانغراس القصيرة والمتوسطة والطويلة، وكذلك البروفيلات السطحية المحددة.
التحقق الصارم: تستخدم الدراسة تقنيات تحقق واسعة النطاق، بما في ذلك 10 جولات مستقلة، والتحقق المتبادل "ترك واحد خارج" (LOOCV)، وإعادة أخذ العينات بطريقة البوتستراب (bootstrap resampling)، لتقييم استقرار النموذج، وقدرة التعميم، وعدم يقين المعلمات.
النتائج
الأداء التنبئي: حقق نموذج WOA–ETR دقة عالية على مجموعة بيانات الاختبار بـ R2 قدره 0.942، وRMSE قدره 1.536 ميجا باسكال، وMAPE قدره 16.812%. كما أظهر استقرارًا وقدرة تعميم فائقتين مقارنة بالنماذج المرجعية.
أهمية الميزات: حدد تحليل SHAP أن درجة حرارة التعرض هي العامل الأكثر تأثيرًا، تليها البروفيل السطحي، ونسبة الغطاء إلى قطر السيخ، ونسبة طول الانغراسا. ومن الجدير بالذكر أن مقاومة ضغط الخرسانة ونوع السيخ أظهرا أهمية نسبية أقل في هذا السياق عالي الحرارة.
السلوك غير الخطي: كشف تحليل PDP أن تدهور قوة التماسك غير خطي للغاية بالنسبة لدرجة الحرارة، حيث يحدث انخفاض سريع بين درجة الحرارة المحيطة وحوالي 350 درجة مئوية، وبعد ذلك يستقر المنحنى. وبالمثل، أظهرت نسبة طول الانغراس ثلاث مناطق استجابة متميزة (مستقرة، انتقالية، وهضبة)، مما يبرر نهج التصنيف.
دقة المعادلات: قدمت المعادلات التنبؤية المشتقة تنبؤات غير منحازة تقريبًا، حيث تراوحت متوسط النسب التجريبية إلى المتوقعة بين 1.02 و1.06. تراوحت قيم R2 للمعادلات المشتقة بين 0.716 و0.986. وأظهرت المعادلات الخاصة بأطوال الانغراس المتوسطة والطويلة أداءً استثنائيًا (R2>0.96) وتغيرًا منخفضًا.
المتانة: أكدت تحليلات البوتستراب وLOOCV أن معادلات الانغراس المتوسط/الطويل والأسياخ المضلعة ذات الطلاء الرملي (SP=5) تمتلك استقرارًا عاليًا في المعلمات. وأظهرت معادلات المجموعات الفرعية الأخرى (مثل الأسياخ المغطاة بالرمل) عدم يقين أعلى في المعلمات بسبب صغر حجم العينات، لكنها حافظت على دقة تنبؤية مقبولة.
الأهمية والادعاءات تزعم الورقة أن أهميتها الأساسية تكمن في سد الفجوة بين نماذج التعلم الآلي عالية الدقة "الصندوق الأسود" والتصميم الهندسي العملي. فمن خلال استخدام التعلم الآلي ليس فقط للتنبؤ ولكن كأداة للاكتشاف لفهم الآليات الحاكمة، نجحت الدراسة في ترجمة أنماط البيانات المعقدة إلى معادلات تحليلية مبسطة وشفافة.
ويؤكد المؤلفون أن إطار العمل المقترح يقدم بديلًا أكثر تكيفًا ودقة للصيغ التجريبية وأكواد التصميم الحالية (مثل ACI 440.1R)، والتي غالبًا ما تعتمد على ظروف تجريبية محدودة. وتشدد الدراسة على أنه بينما يعمل نموذج WOA–ETR كمقدر قوي، فإن المعادلات المشتقة توفر أداة فعالة حاسوبيًا وقابلة للتفسير فيزيائيًا للمهندسين لتقدير قوة التماسك تحت درجات الحرارة المرتفعة. ومع ذلك، يشير المؤلفون بتواضع إلى أن موثوقية بعض المعادلات مقيدة بالعدد المحدود من العينات في مجموعات بيانات فرعية محددة (خاصة لأطوال الانغراس الطويلة وبروفيلات سطحية معينة)، مما يشير إلى أن الحصول على المزيد من البيانات التجريبية سيكون مفيدًا لتنقيح هذه الصيغ وتوسيع نطاق تطبيقها.