Multiscale Hyperelastic Characterization of Polyurethane-SiC Nanocomposites under Large Compressive Deformation: A Molecular Dynamics and Constitutive Modeling Approach
تدمج هذه الدراسة محاكاة الديناميات الجزيئية مع نمذجة "أوجدن" فائقة المرونة المُحسّنة لتوصيف السلوك الانضغاطي غير الخطي المعتمد على معدل الانفعال لمركبات "البولي يوريثين-سيليكون كاربيد" النانوية، مما ينجح في الربط بين الآليات النانوية والموصوفات التأسيسية على مستوى المتصل لاستخدامات هندسية متقدمة.
المؤلفون الأصليون:Sy-Ngoc Nguyen, Hongdeok Kim, Joonmyung Choi
تخيل أنك تحاول بناء بدلة البطل الخارق المثالية. تريدها أن تكون ناعمة ومطاطية مثل الرباط المطاطي لتتمكن من التحرك بحرية، ولكن قوية بما يكفي لإيقاف رصاصة مسرعة. هذا هو حلم علم المواد: ابتكار "المواد النانوية المركبة"، وهي تشبه خلط ذرات صغيرة فائقة القوة من الصخر في مادة بلاستيكية لزجة. الهدف هو الحصول على أفضل ما في العالمين. لفهم كيفية سلوك هذه المواد، ينظر العلماء إليها غالبًا تحت المجهر، ولكنهم يحتاجون أحيانًا إلى التكبير أكثر من ذلك، وصولاً إلى مستوى الذرات الفردية. وهنا يأتي دور "الديناميكا الجزيئية". فكر في الأمر كأنه فيلم سينمائي فائق القوة على الكمبيوتر يعرض رقصة الذرات بالحركة البطيئة، مما يسمح للعلماء بمشاهدة كيف تنضغط المادة، وتتمدد، وتعود لشكلها الأصم دون كسر عينة حقيقية أبدًا. فكرة رئيسية أخرى هي "المرونة الفائقة". بينما يعود الرباط المطاطي البسيط لشكلة بسهولة، فإن المواد ذات المرونة الفائقة هي "بهلوانيات" هذا العالم؛ إذ يمكن ضغطها لتصبح كرة صغيرة أو تمديدها لضعف حجمها ومع ذلك تعود إلى شكلها الأصلي. إن فهم كيفية قيامها بذلك بدقة أمر بالغ الأهمية لتصميم كل شيء، بدءًا من مصدات السيارات التي تمتص الصدمات وصولًا إلى الغرسات الطبية التي تتحرك مع جسمك.
تتعمق هذه الورقة البحثية في دراسة مادة بطل خارق محددة: مزيج من "البولي يوريثين" (نوع من البلاستيك) و"كربيد السيليكون" (SiC)، وهي جسيمات نانوية صلبة كالصخر. أراد الباحثون معرفة ما يحدث عندما يصطدم هذا النوع من المواد بقوة هائلة، عبر ضغطها بنسبة هائلة تصل إلى 75% — تخيل ضغط قطعة من "المارشميلو" حتى تصبح مسطحة كالفطيرة. استخدموا فيلمهم الحاسوبي (الديناميكا الجزيئية) لمحاكاة هذا السحق بسرعات هائلة، وتحديدًا بمعدلات انفعال تبلغ 109 و 1010 و 1011 في الثانية. هذه السرعات تفوق ملايين المرات أي شيء يمكن للإنسان القيام به في المختبر، لكنها ضرورية لرؤية كيفية تفاعل الذرات في أجزاء من الثانية الضئيلة التي يمكن للحاسوب التعامل معها.
اكتشف الفريق أن إضافة جسيمات "كربيد السيليكون" النانوية تعمل مثل إضافة مراسي صغيرة غير مرئية داخل البلاستيك الناعم. فعندما تُضغط المادة بلطف، يقوم البلاستيك بمعظم العمل. ولكن بمجرد أن يصبح الضغط عميقًا حقًا (بعد تجاوز 60% من الانضغاط)، تبدأ الجسيمات الصلبة في الاصطدام ببعضها البعض، مما يؤدي إلى حبس سلاسل البلاستيك في مكانها وجعل المادة تصبح فجأة أكثر صلابة وقوة. وكلما زادت سرعة الضغط، زادت قوة مقاومة المادة، لأن سلاسل البلاستيك لم يكن لديها الوقت الكافي للتحرك والابتعاد عن الطريق.
ولفهم هذه الحركات اللزجة والمجنونة، حاول العلماء مطابقة البيانات مع معادلات رياضية، تمامًا مثل محاولة وصف رحلة "الأفعوانية" بمعادلة واحدة. لقد اختبروا عدة معادلات مختلفة، بما في ذلك بعض المعادلات البسيطة ونموذجًا أكثر تعقيدًا يسمى "نموذج أودجن" (Ogden model). ووجدوا أن المعادلات البسيطة فشلت فشلًا ذريعًا؛ إذ كانت إما تظن أن المادة ستكون لينة جدًا أو تتوقع انفجار الضغط بشكل هائل بينما كانت المادة في الواقع تزداد صلابة فقط. ومع ذلك، كان "نموذج أودجن" المعقد ذو الحدين متطابقًا تمامًا. لقد نجح في رصد الرحلة بأكملها: الضغط الأولي السهل، ثم القسم الأوسط حيث تزداد الصعوبة، وصولًا إلى الارتفاع الدراماتيكي النهائي في الضغط عندما تصطدم الجسيمات النانوية ببعضها البعض.
تشير الدراسة إلى أنه باستخدام هذا النموذج الرياضي المحدد، يمكن للمهندسين الآن التنبؤ بدقة بكيفية سلوك هذه المادة في سيناريوهات العالم الحقيقي، مثل حماية هيكل خفيف الوزن من تأثير اصطدام أو تصميم غرسة طبية حيوية يمكنها تحمل حركات الجسم. ويشعر الباحثون بالثقة في أن طريقتهم تعمل لأن الرياضيات طابقت محاكاتهم الحاسوبية بدقة عالية جدًا، حيث وصلت نسبة الخطأ إلى 0.67% في أفضل الحالات. كما لاحظوا أنه بينما تصبح المادة صلبة للغاية عند السرعات العالية، فإن قدرتها على تغيير الحجم (كيف تتقلص عرضيًا عند الضغط) تظل ثابتة تقريبًا، مما يعني أنها تظل "إسفنجة" موثوقة وغير قابلة للانضغاط حتى تحت الضغط الشديد. هذا العمل يسد الفجوة بين عالم الذرات الصغير والفوضوي وعالم الهندسة الكبير والعملي، مقدمًا مخططًا جديدًا لبناء مواد تجمع بين النعومة والقوة التي لا تُقهر.
بيان المشكلة تعد متراكبات النانو من البولي يوريثين وكربيد السيليكون (PU-SiC) واعدة للتطبيقات الهندسية المتقدمة نظرًا لخصائصها الميكانيكية المعززة. ومع ذلك، فإن توصيف سلوكها الميكانيكي تحت تشوهات الضغط الكبيرة يفرض تحديات كبيرة. إن التفاعل بين المرونة الفوقية (Hyperelasticity) والاعتماد على معدل الانفعال، لا سيما عند الكسور الحجمية المنخفضة لكربيد السيليكون (SiC)، لا يزال غير مستكشف بشكل كافٍ. تفشل النماذج المرنة الخطية التقليدية في التقاط العلاقات المعقدة وغير الخطية للإجهاد والانفعال المتأصلة في هذه المواد. علاوة على على ذلك، غالبًا ما تعاني استراتيجيات مطابقة المرونة الفوقية التقليدية من أجل تحقيق معاملات مستقرة وذات معنى فيزيائي عبر نطاقات تشوه واسعة، خاصة في المتراكبات النانوية الهجينة (لينة-صلبة) حيث يؤدي التباين بين المصفوفة البوليمرية والجسيمات النانوية الصلبة إلى عدم خطية شديدة. بالإضافة إلى ذلك، تفتقر الطرق التجريبية غالبًا إلى الدقة اللازمة لتصور توزيعات الإجهاد والآليات على مقياس النانو التي تقود الاستجابة المادية.
المنهجية تستخدم هذه الدراسة إطار عمل متعدد المقاييس يدمج محاكاة الديناميكا الجزيئية (MD) مع النمذجة التأسيسية للمرونة الفوقية على مستوى المتصل (Continuum level).
محاكاة الديناميكا الجزيئية (MD): أُجريت محاكاة الديناميكا الجزيئية لجميع الذرات للتحقيق في متراكبات (PU-SiC) تحت ضغط أحادي المحور. تم نمذجة ثلاثة أنظمة مواد: بولي يوريثين نقي (0% SiC)، وNC1 (1.95% SiC)، وNC2 (3.9% SiC). غطت المحاكاة انفعالات هندسية ضاغطة تصل إلى 75% (ϵ=−0.75) عند معدلات انفعال تبلغ 109 و1010 و1011 s−1. استخدمت النماذج جسيمات نانوية كروية من كربيد السيليكون (بقطر 2.5 نانومتر) داخل مصفوفة بولي يوريثين مكونة من ثنائي ميثيلين ثنائي فينيل ثنائي أيزوسياناتات وبولي تتتراميثيلين إيثر جليكول.
توليد البيانات: ولدت المحاكاة منحنيات إجهاد-انفعال عالية الدقة ونسب الانفعال المستعرض إلى المحوري. تم وضع تركيز خاص على تباعد الإجهاد الملاحظ عند انفعالات الضغط العالية، والذي يُعزى إلى التفاعلات بين الجسيمات النانوية وصورها الدورية (محاكاة التلامس بين الجسيمات النانوية).
النمذجة التأسيسية: تم مطابقة البيانات المستمدة من الديناميكا الجزيئية مع دالات طاقة الانفعال (SEFs). قيمت الدراسة نماذج "نيو-هوكيان" القابلة للانضغاط (CnH)، ونموذج "أوجدن" أحادي الحد (Ogd1)، ونموذج "أوجدن" ثنائي الحد (Ogd2).
التحسين: استُخدمت خوارزمية التحسين البديلة (SOA) لتحديد معاملات المرونة الفوقية بقوة. نهج هذا الأسلوب تقليل جذر متوسط مربع الخطأ (RMSE) بين الإجهادات الهندسية المستمدة من الديناميكا الجزيئية وتنبؤات النموذج مع فرض قيود فيزيائية تتعلق بمعاملات القص والحجم الأولية. عالج التحسين تحديدًا تحدي مطابقة الارتفاعات الحادة في الإجهاد الناتجة عن تفاعلات الجسيمات النانوية.
المساهمات الرئيسية
إطار عمل متعدد المقاييس: تؤسس هذه الدراسة رابطًا كميًا بين آليات التشوه على المقياس الذري (MD) والأوصاف التأسيسية الماكروسكوبية (نماذج المرونة الفوقية) لمتراكبات (PU-SiC) تحت الضغط الشديد.
خوارزمية التحسين: تم تطوير طريقة تحديد قوية قائمة على التحسين لالتقاط الاستجابات الميكانيكية غير الخطية، بما في ذلك تباعد الإجهاد الناتج عن تلامس الجسيمات النانوية عند الانفعالات العالية، والتي غالبًا ما تفشل نماذج المطابقة التقليدية في تمثيلها بدقة.
اختيار النموذج: تم تحديد نموذج "أوجدن" ثنائي الحد (Ogd2) كأكثر دالة طاقة انفعال (SEF) ملاءمة من خلال مقارنة معيارية منهجية. نجح النموذج في إعادة إنتاج استجابة الإجهاد-الانفعال كاملة النطاقات، من النظام الإنتروبي اللين الأولي إلى التصلب الحاد عند الانفعال الكبير، متفوقًا على نموذجي CnH وOgd1.
التحليل البارامتري: توفر الدراسة معاملات المرونة الفوقية المحسنة كدوال للكسر الحجمي لكربيد السيليكون ومعدل الانفعال، مما يقدم وصفًا كميًاًا لتطور الصلابة والقابلية للانضغاط.
النتائج
آلية التعزيز: تعمل جسيمات كربيد السيليمان النانوية على تعزيز صلابة مصفوفة البولي يوريثين بشكل كبير، لا سيما عند الانفعالات الكبيرة (ϵ<−0.6). تعمل الجسيمات النانوية كنقاط تثبيت لسلاسل البولي يوريثين، مما يحد من الحركة ويحسن نقل الإجهاد. يتكثف هذا التأثير مع زيادة محتوى كربيد السيليكون (NC2 > NC1 > Pure PU).
الاعتماد على معدل الانفعال: تؤدي معدلات الانفعال الأعلى (1011 s−1) إلى مستويات إجهاد مرتفعة بسبب محدودية الوقت اللازم لارتخاء السلاسل وتفككها، مما يزيد فعليًا من الصلابة الظاهرية ويدفع المادة نحو نظام يشبه الحالة الزجاجية. يتضخم هذا التأثير في المتراكبات النانوية ذات المحتوى الأعلى من كربيد السيليكون.
الاستجابة المستعرضة: تظل نسبة الانفعال المستعرض إلى المحوري (نسبة بواسون) غير متأثرة بشكل كبير بالكسور الحجمية لكربيد السيليكون حتى 3.9%، وهو ما يتوافق مع الطبيعة شبه غير القابلة للانضغاط لمصفوفة البولي يوريثين. تختلف النسبة قليلاً مع معدل الانفعال ولكنها تظل قريبة من 0.35.
تباعد الإجهاد: عند انفعالات الضغط التي تتجاوز حوالي 60% (تحديدًا ϵ≈−0.60 لـ NC1 وϵ≈−0.55 لـ NC2)، يحدث تباعد حاد في الإجهاد. يُعزى ذلك إلى التفاعل بين الجسيمات النانوية وصورها الدورية، مما يحاكي شبكة تلامس تشبه التغلغل (Percolation-like contact network). نموذج Ogd2 المقترح، الذي تمت مطابقته عبر SOA، يلتقط هذا السلوك بقيم RMSE منخفضة (تتراوح بين 0.67% إلى 5.63%).
أداء النموذج: أظهر نموذج Ogd2 دقة فائقة عبر جميع تركيبات المواد ومعدلات الانفعال. في المقابل، بالغ نموذج CnH بشكل كبير في تقدير الإجهاد عند الانفعالات الكبيرة، بينما قلل نموذج Ogd1 من تقدير ارتفاع الإجهاد المدفوع بتفاعلات الجسيمات النانوية.
الأهمية تزعم الورقة أن هذا العمل يوفر منهجية تنبؤية لتصميم وتحسين متراكبات (PU-SiC). ومن خلال الربط بين الآليات النانوية والنماذج التأسيسية الماكروسكوبية، تقدم الدراسة رؤى أساسية ذات صلة بالتطبيقات التي تتطلب هياكل خفيفة الوزن عالية الأداء، والطلاءات المقاومة للصدمات، والغرسات الطبية الحيوية. إن القدرة على نمذجة التصلب غير الخطي والحساسية لمعدل الانفعال لهذه المواد عند التشوهات الكبيرة تعالج فجوة حرجة في مناهج التوصيف الحالية، مما يسمح بتنبؤات هندسية أكثر موثوقية للمواد النانوية اللينة تحت ظروف التحميل القصوى.