Finite strain homogenization of periodic rod networks with application to semi-flexible biopolymers
تقدم هذه الورقة إطار عمل للتجانس الحسابي لللانفعال المحدود بناءً على نظرية قضيب كوسيرات دقيقة هندسياً لتوصيف الاستجابات الميكانيكية غير الخطية لشبكات البوليمرات الحيوية شبه المرنة الدورية، حيث نجحت في رصد ظواهر مثل التصلب بالانفعال، والانكماش الحجمي، وتأثير بوينتينغ العكسي، مع إظهار توافق قوي مع البيانات التجريبية والقدرة على نمذجة البنى المترية المرنة.
المؤلفون الأصليون:Vinayak, Prashant K. Purohit, Ajeet Kumar
تُعد الطبيعة أستاذة في الهندسة، حيث تبني كل شيء بدءاً من الصدفة الصلبة للكائنات البحرية وصولاً إلى الأنسجة الرخوة واللينة داخل أجسادنا. هذه المواد ليست كتلًا صلبة؛ بل هي شبكات معقدة مكونة من ألياف دقيقة، تشبه إلى حد كبير شبكة مجهرية. بعض هذه الألياف صلبة، بينما البعض الآخر مرن بما يكفي للانحناء والالتواء. وعندما يحاول العلماء فهم كيفية سلوك هذه الشبكات البيولوجية، فإنهم يواجهون لغزاً صعباً: فالألياف تتحرك بطرق معقدة تتغير اعتماداً على مقدار السحب أو الضغط أو الالتواء الذي تتعرض له. وللتنبؤ بكيفية استجابة قطعة من النسيج لقوة ما، يجب على الباحثين النظر في سلوك هذه الألياف الفردية، ثم اكتشاف كيف تخلق حركتها الجماعية قوة أو نعومة المادة بأكملها. وهذا يمثل تحدياً لأن الألياف يمكن أن تنبعج، أو ترتد، أو تتمدد بطرق يصعب حسابها دون استخدام حواسيب قوية.
لقد طور فريق من الباحثين طريقة جديدة لحل هذا اللغز من خلال إنشاء نموذج افتراضي لشبكات الألياف هذه. وبدلاً من محاولة محاكاة كل ليف بمفرده في قطعة كبيرة من النسيج، وهو أمر يتطلب قدرة حوسبية هائلة، ركزوا على نمط صغير متكرر يمثل الهيكل بأكره. لقد تخيلوا مكعباً صغيراً يحتوي على عدد قليل من الألياف المرتبة في أشكال محددة، مثل نجمة ذات ثمانية أذرع أو نسخة أكثر تعقيداً ذات أربعة عشر ذراعاً. وباستخدام نهج رياضي يعامل هذه الألياف كقضبان مرنة تماماً، قاموا بمحاكاة ما يحدث عندما يتم شد هذا المكعب الصغير أو ضغطه أو قصّه. ومن خلال مراقبة كيفية تحرك وانحناء هذه الألياف القليلة داخل الصندوق الافتراضي، تمكنوا من التنبؤ بكيفية سلوك شبكة ضخمة من ألياف مماثلة في العالم الحقيقي.
وجد الباحثون أن الطريقة التي تستجيب بها هذه الشبكات للقوة تعتمد بشكل كبير على ما إذا كان يُسمح للألياف بالانحناء أو إذا كانت مجبرة على البقاء مستقيمة. فعندما سحبوا شبكتهم الافتراضية بطريقة تحاكي شد قطعة من النسيج، انحنت الألياف والتوت في البداية. ومع استمرار السحب، استقامت الألياف وبدأت في التمدد، مما جعل المادة تصبح أكثر صلابة فجأة. هذا "التصلب الناتج عن الانفعال" هو سمة شائعة في الأنسجة البيولوجية مثل الجلطات الدموية والجلد، مما يسمح لها بحماية الجسم من التمزق. كما كشف النموذج أنه مع تمدد الشبكة، فإن حجمها ينكمش بشكل ملحوظ، وهي ظاهرة تساعد في تفسير كيفية قدرة الأنسجة على الحفاظ على شكلها أثناء السحب.
وعندما قام الباحثون بضغط الشبكة بدلاً من سحبها، تغير السلوك بشكل جذري. فقد انبعجت الألياف، أو انحنت جانباً، مما تسبب في تليين المادة بدلاً من تصلبها. هذا الفقدان المفاجئ للصلابة هو ما يسمح للجلطات الدموية بالانضغاط دون أن تتفكك فوراً. واكتشف الفريق أن حجم الصندوق الافتري الذي استخدموه كان مهماً؛ حيث أظهرت الصناديق الأكبر أن الألياف تنبعج عند ضغوط أقل، مما يطابق التجارب الواقعية على الجلطات الدموية بشكل وثيق. كما اختبروا نسخة من النموذج حيث لا تكون الألياف مستقيمة بل ملفوفة مثل الزنبركات. وسمحت هذه الألياف الملفوفة للشبكة بالتمدد لمسافة أكبر بكول من التصلب، مما يحاكي سلوك بعض المواد البيولوجية التي تحتاج إلى مرونة عالية.
وفي اختبار آخر، قام الباحثون بليّ الشبكة لمحاكاة حركة القص، مثل انزلاق طبقة من مجموعة أوراق اللعب فوق أخرى. وقد لاحظوا تأثيراً غريباً حيث تحاول المادة الانكماش في اتجاه الليّ، مما يخلق قوة سحب يجب مواجهتها للحفاظ على استقرار الشكل. هذا "تأثير بوينتينغ العكسي" هو ميزة فريدة أخرى للهلاميات البيولوجية، وقد نجح النموذج في التقاطه من خلال إظهار كيفية انبعاج الألياف وإعادة ترتيب نفسها تحت تأثير الالتواء. وأظهرت عمليات المحاكاة أن انحناء وانبعاج هذه الألياف الصغيرة هو السبب الرئيسي لهذه السلوكيات غير العادية، وليس مجرد تمدد الألياف نفسها.
قارن الفريق بين نتائج الحاسوب وتجارب حقيقية على الكولاجين والفايبرين، وهما البروتينات التي تشكل الأنسجة الضامة والجلطات الدموية. وقد طابق النموذج الافتراضي البيانات التجريبية بشكل جيد جداً، مما أكد أن نهجهم يمكنه التنبؤ بدقة بكيفية تفاعل هذه المواد المعقدة مع القوى المختلفة. كما أظهروا أنه من خلال تغيير شكل الألياف في النموذج، يمكنهم تصميم مواد جديدة بخصائص محددة، مثل القدرة على التمدد لمسافات أكبر بكثير من الأنسجة الطبيعية. يوفر هذا العمل أداة قوية لفهم ميكانيكا الحياة على المستوى المجهري، ويقدم مخططاً لإنشاء مواد جديدة مستوحاة حيوياً يمكنها محاكاة مرونة وقدرة التصاميم الطبيعية على التكيف.
بيان المشكلة تُظهر شبكات البوليمرات الحيوية شبه المرنة، مثل الكولاجين والفايبرين، سلوكيات ميكانيكية معقدة وعالية اللاخطية، بما في ذلك التصلب عند الانفعال، والانكماش الحجمي تحت تأثير الشد، واللين تحت تأثير الضغط، وتأثير بوينتينج العكسي تحت تأثير القص. تنشأ هذه الظواواهر من الطبيعة الهرمية للمواد والتفاعل بين التشوه الماكروسكوبي (الجهري) وميكانيكا الألياف على المستوى الميكروي (المجهري)، وتحديداً الانحناء والبعج (buckling). وبينما توجد نماذج متصلة ظواهرية (مثل نماذج السلاسل المتوافقة واللا-متوافقة)، إلا أنها غالباً ما تفتقر إلى الرؤية حول الآليات الميكروية المحددة التي تقود هذه السلوكيات. ومن ناحية أخرى، فإن المحاكاة الحاسوبية كاملة النطاق لشبكات الألياف العشوائية مكلفة حاسوبياً. يعالج هذا العمل الحاجة إلى إطار عمل حاسوبي يلتقط التشوهات الميكروية الكبيرة وغير المتوافقة (بما في ذلك البعج) مع الحفاظ على الكفاءة الحاسوبية من خلال عملية التجانس.
المنهجية يقترح المؤلفون نهج تجانس حاسوبي للانفعال المحدود مطبق على عناصر حجم تمثيلية (RVEs) دورية مكونة من 8 سلاسل و14 سلسلة. تدمج المنهجية المكونات التالية:
النمذجة الميكروية: يتم نمذجة الألياف الفردية باستخدام نظرية قضبان كوسيرات (Cosserat) الجيومترية الدقيقة (GE). يسمح هذا الصياغة بالتقاط التشوهات الكبيرة للغاية، بما في ذلك التمدد، والانحناء، والالتواء، والقص، دون القيود المرتبطة بالتشوهات الصغيرة لنظريات عارضة أويلر-بيرنولي أو تيموشينكو.
إطار التجانس: تمت صياغة المشكلة باستخدام نهجين: الموجه بالانفعال والموجه بالإجهاد.
الموجه بالانفعال: يتم تحديد تدرج التشوه الماكروسكوبي (FM)، ويتم تقليل طاقة الـ RVE مع مراعاة القيود.
الموجه بالإجهاد: يتم تحديد إجهاد فيما الأول لكيرشوف الماكروسكوبي (PM)، مما يجعل تدرج التشوه الماكروسكوبي مجهولاً. وهذا أمر بالغ الأهمية لنمذجة الشروط الحدية الجانبية الخالية من الإجهاد (traction-free) الشائعة في التجارب.
القيود والشروط الحدية: تخضع الـ RVEs لقيود الوصلات الداخلية (صلبة أو كروية/كرة ومقبس)، وشروط حدية دورية (PBCs) مشتقة من قاعدة كوشي-بورن، وقيود عالمية لإزالة أنماط الجسم الصلب.
الحل العددي: يتم حل مشكلة تقليل الطاقة المقيدة باستخدام مخطط نيوتن-رافسون المعزز باستمرار طول القوس (arc-length continuation) لتتبع المسار. يتيح هذا تتبع الحلول ما بعد البعج دون إدخال عيوب اصطناعية.
تحليل الاستقرار: يتم تقييم استقرار الحلول التي تم الحصول عليها باستخدام معيارين: "التحميل الصلب" (اضطرابات في متغيرات التكوين فقط) و"التحميل اللين" (اضطرابات في كل من متغيرات التكوين وتدرج التشوه الماكروسكوبي). المعيار الأخير هو شرط أكثر صرامة.
المساهمات الرئيسية
استعادة ظواهر البوليمرات الحيوية الفريدة: تُظهر الدراسة أن شبكات القضبان ذات الـ 8 والـ 14 سلسلة المجانسة يمكنها استعادة ظواهر محددة لوحظت في الأنسجة البيولوجية، بما في ذلك:
التصلب عند الانفعال والانكماش الحجمي الهائل تحت تأثير الشد أحادي المحور.
اللين المفاجئ تحت تأثير الضغط.
تأثير بوينتينج العكسي (توليد إجهاد شد طبيعي تحت تأثير القص البسيط).
تحديد الآلية: حدد المؤلفون أن اللاخطية والتشوه غير المتوافق (non-affine) على المستوى الميكروي، وتحديداً انحناء وبعج الألياف، هي المحركات الأساسية لهذه الظواهر الماكروسكوبية.
تحليل ما بعد البعج الخالي من العيوب: حصل العمل على حلول ما بعد البعج باستخدام نهج تتبع مسار غير خطي دون الاعتماد على العيوب الهندسية الأولية، مع التحقق بدقة من الاستقرار عبر معايير التحميل الصلب واللين.
التحقق والتوسيع:
تم التحقق من النموذج مقابل بيانات تجريبية لكولاجين الشد أحادي المحور (Roeder et al., 2002) ونماذج تحليلية للفايبرين (Brown et al., 2009).
تمت مقارنة استجابة الضغط مع تجارب جلطات الدم (Liang et al., 2017)، والتي أظهرت توافقاً جيداً عند مراعاة حجم الـ RVE ونحافة الألياف.
تم توسيع الإطار ليشمل القضبان الحلزونية داخل وحدة الخلية ذات الـ 8 سلاسل، مما نجح في التقاط سلوك التمدد الموسع الملاحظ في المواد المبتكرة (metamaterials) المرنة.
النتائج
الشد أحادي المحور:
الموجه بالانفعال: الاستجابة تهيمن عليها عملية التمدد وهي خطية.
الموجه بالإجهاد: تُظهر الاستجابة تصلباً عند الانفعال. في البداية، يكون التشوه مهيمناً بالانحناء؛ ومع محاذاة الألياف مع اتجاه الشد، ينتقل النمط إلى التمدد. يتسبب هذا الانتقال في تأثير بواسون غير خطي وانكماش حجمي كبير.
14 سلسلة مقابل 8 سلاسل: يُظهر نموذج الـ 14 سلسلة تشعباً متميزاً حيث تبتعد (buckle) القضبان الجانبية، مما يؤدي إلى لين مفاجئ وتقليل حجمي دراماتيكي، وهو ما يتناقض مع الزيادة الحجمية الرتيبة التي يشهدها مسار الـ 8 سلاسل غير المبتعج.
الضغط:
تُظهر كل من الـ RVEs منطقة خطية تليها مرحلة ثبات (لين) بسبب بعج الألياف.
يُظهر نموذج الـ 14 سلسلة استجابة ثلاثية الخطوط: نظام خطي أولي، ثم أول مرحلة ثبات (بعج القضبان المستقيمة)، ثم ثاني مرحلة ثبات (بعج القضبان المائلة).
يعيد النموذج إنتاج سلوك الرغوة (تدرج معامل يونغ والإجهاد الحرج مع الكسر الحجمي) الملاحظ في شبكات الفايبرين.
القص البسيط:
يُظهر كلا النموذجين استجابة إجهاد-انفعال ثنائية الخط. عند الوصول إلى نقطة التشعب، تبتعد القضبان القطرية للجسم، مما يسبب ليناً في الإجهاد.
لوحظ تأثير بوينتينج عكسي كبير: مع ابتعاد القضبان وانكماشها في الاتجاه العمودي، يتولد إجهاد شد طبيعي للحفاظ على الحدود الدورية (F11=1). وتكون قيمة هذا الإجهاد الطبيعي قابلة للمقارنة مع إجهاد القص.
الأهمية والادعاءات يزعم البحث أن أهميته الرئيسية تكمن في توفير إطار عمل تجانس حاسوبي دقيق وخالٍ من العيوب يربط بين ميكانيكا الألياف الميكروية والاستجابة النسيجية الماكروسكوبية. من خلال استخدام قضبان جيومترية دقيقة، يلتقط النموذج التشوهات الكبيرة وتحولات الاستقرار (البعج) التي تعد حاسمة لسلوك البوليمرات الحيوية ولكن غالباً ما تغفلها النماذج الأبسط. يؤكد المؤلفون أن نهجهم يستعيد ظواهر فريدة مثل الانكماش الحجمي وتأثير بوينتينج العكسي حصرياً من خلال ميكانيكا بعج وانحناء الألياف، دون استدعاء معاملات ملاءمة ظواهرية. كما أن النجاح في التحقق من النموذج مقابل البيانات التجريبية للكولاجين وجلطات الدم، بالإضافة إلى القدرة على نمذجة المواد المبتكرة الحلزونية، يشير إلى فائدة الإطار لفهم الظواهر الفيزيائية الحيوية الحالية وتصميم المواد المهندسة المستوحاة حيوياً في المستقبل. يشير المؤلفون بتواضع إلى أن النموذج الحالي مرن محلي بحت، وأن العمل المستقبلي سيحتاج إلى دمج عدم المرونة، وتلامس الألياف، والتأثيرات غير الموضعية لالتقاط سلوك الشبكات الكثيفة بشكل كامل.