A Hybridized Staggered Discontinuous Galerkin--Mixed Finite Element Method for Strain Gradient Elasticity
تقترح هذه الورقة طريقة عناصر منتهية مختلطة بأسلوب "غاليركين المتقطعة المتدرجة الهجينة" لمرونة تدرج الانفعال، والتي تحقق استقراراً متيناً تجاه المعلمات وتقارباً أمثلاً من خلال إعادة صياغة المسألة من الرتبة الرابعة إلى نظام من الرتبة الأولى مع فرض التماثل القوي وصياغة هجينة تقلل من الاقتران العالمي.
تخيل عالماً تتغير فيه قواعد كيفية انحناء المواد وتمددها اعتماداً على مدى صغر الحجم الذي تنظر إليه. لأكثر من قرن، اعتمد المهندسون على المرونة الكلاسيكية، وهي مجموعة من القوانين التي تعامل المواد كنسيج مستمر وسلس. تعمل هذه القوانين بشكل مثالي للجسور والمباني والآلات الكبيرة، ومع ذلك، فإنها تفشل عند تطبيقها على الهياكل المجهرية الموجودة في التكنولوجيا الحديثة، مثل الأسلاك الدقيقة في الرقائق الإلكترونية أو العوارض الرقيقة في المستشعرات الطبية. عند هذه المقاييس، تظهر المواد "تأثيرات الحجم"، حيث تسلك سلوكاً مختلفاً لمجرد أنها صغيرة. فقد يلتوي سلك رفيع بقوة أكبر مما تتنبأ به الفيزياء الكلاسيكية، أو قد يقاوم عارضة مجهرية الانحناء بطرق غير متوقعة. ولتفسير ذلك، طور العلماء "مرونة تدرج الانفعال"، وهي نظرية أكثر تقدماً تضيف "مقياس طول" محدداً إلى المعادلات. هذا المعامل الإضافي يأخذ في الاعتبار البنية الداخلية للمادة، مما يسمح للرياضيات بالتقاط تلك السلوكيات الغامضة المعتمدة على الحجم.
ومع ذلك، فإن التحدي يكمن في أن هذه المعادلات المتقدمة صعبة الحل للغاية على الحاسوب. فهي مسائل من الدرجة الرابعة، مما يعني أنها تتضمن علاقات معقدة بين كيفية تشوه المادة وكيف يتغير هذا التشوه عبر الفضاء. وغالباً ما تتعثر الطرق التقليدية لحل هذه المسائل عندما تكون المادة غير قابلة للانضغاط تقريباً (مثل المطاط) أو عندما يكون مقياس الطول صغيراً جداً، مما يؤدي إلى نتائج غير دقيقة أو ظهور "طبقات حدودية عددية" حيث يبتكر الحاسوب أخطاءً اصطناعية بالقرب من الحواف. لقد جرب الباحثون مناهج متنوعة، لكن العديد منها يتطلب ضبطاً دقيقاً للمعاملات أو يفشل في توفير طريقة مباشرة لحساب الإجهاد الكلي، وهو كمية حيوية للمهندسين الذين يصممون الأجهزة المجهرية الآمنة.
في دراسة جديدة، اقترح الباحثان بوهان يانغ وإريك ت. تشونج طريقة جديدة لمعالجة هذه المشكلة. فقد طورا طريقة هجينة تجمع بين تقنيتين رياضيتين قويتين: طريقة "غاليركين" المتدرجة غير المتصلة وطريقة العناصر المحدودة المختلطة. وبدلاً من محاولة الحل لإزاحة المادة وحدها، يقدم نهجهما عدة متغيرات جديدة لتمثيل أنواع الإجهاد المختلفة داخل المادة. حيث يعاملان الإجهاد الكلي، جنباً إلى جنب مع نسخ مقاسة من القوى الداخلية، كمتغيرات مستقلة. ومن خلال القيام بذلك، أعادا صياغة مسألة الدرجة الرابعة الصعبة إلى نظام من معادلات الدرجة الأولى الأبسط. هذا التحول يسمح لهما باستخدام أنواع محددة من المساحات الرياضية التي تتعامل طبيعياً مع تعقيد المسألة دون الحاجة إلى إصلاحات اصطناعية أو حدود جزائية.
يكمن جوهر ابتكارهما في كيفية تنظيم الحساب. فقد قسما هندسة المادة إلى شبكة من المثلثات ثم أنشآ شبكة ثانية متداخلة داخلها. وعلى هذا الهيكل المزدوج، فرضا قواعد مختلفة للاستمرارية، مما يسمح للقطع الرياضية بالتلاحم بطريقة تحافظ على التماثل الفيزيائي للإجهاد. يضمن هذا التصميم بقاء الطريقة مستقرة ودقيقة بغض النظر عن مدى صغر مقياس طول المادة أو مدى قرب المادة من كونها غير قابلة للانضغاط. ومن السمات الرئيسية لعملهما أنه يمكنهما إثبات متانة الطريقة رياضياً؛ حيث لا تتدهور الدقة عند تغير المعاملات، وهو نقطة فشل شائعة في الطرق الأخرى. علاوة على ذلك، قدما نسخة "مهجنة" من مخططهما تسمح لهما بحل معظم المتغيرات محلياً داخل كل مثلث صغير، مما يترك فقط عدداً قليلاً من المتغيرات العالمية ليتم حلها معاً. وهذا يقلل التكلفة الحسابية بشكل كبير، مما يجعل الطريقة فعالة بما يكفي للاستخدام العملي.
اختبر الباحثون طريقتهم من خلال تجارب عددية، حيث قاموا بمحاكاة سيناريوهات تتضمن حلولاً سلسة وأخرى تتميز بطبقات حدودية قوية، وهي تغيرات حادة في السلوك بالقرب من حواف المادة. وفي كل حالة، حققت طريقتهم معدلات الدقة المتوقعة وظلت مستقرة عبر نطاق واسضع من خصائص المواد، بما في ذلك القيم القصوى للصلابة وعدم القابلية للانضغاط. وخلافاً لطرق الإزاحة السابقة التي تعاني عندما يكون مقياس الطول ضئيلاً جداً، حافظ نهجهما على الدقة المثلى دون خلق أخطاء اصطناعية. كما أظهرا أن طريقتهما توفر حساباً مباشراً ودقيقاً للإجهاد الكلي، وهي نتيجة غالباً ما يصعب الحصول عليها باستخدام التقنيات الأخرى. وبينما يركز العمل الحالي على الأشكال ثنائية الأبعاد مع شروط حدودية محددة، فإن نجاح هذا النهج يشير إلى مسار واعد للمضي قدماً في نمذجة المواد المعقدة ذات المقياس المجهري. وتؤكد الدراسة أنه من خلال إعادة التفكير في كيفية هيكلة هذه المعادلات وحلها، فإنه من الممكن التغلب على الصعوبات طويلة الأمد لمرونة تدرج الانفعال، مما يوفر أداة موثوقة للجيل القادم من الهندسة المجهرية.
بيان المشكلة تتناول الورقة البحثية التقريب العددي لنموذج مرونة تدرج الانفعال (SGE)، وتحديداً الصيغة التي اقترحها "أيفانتيس" ومعاونوه. هذا النموذج هو مسألة اضطراب منفرد من الدرجة الرابعة، محكومة بمعلمة طول مادية ι∈(0,1] وثوابت لاميه λ وμ. وتطرح المعادلة الحاكمة، −(1−ι2Δ)(μΔu+(λ+μ)∇divu)=f، تحديين رئيسيين للتقطيع:
الاستمرارية عالية الرتبة: بما أنها مسألة من الدرجة الرابعة، فإن طرق العناصر المحدودة القياسية تتطلب توافقاً من نوع H2. وعلى الرغم من وجود عناصر متوافقة مع C1، إلا أنها معقدة في البناء ومكلفة حاسوبياً. أما الطرق غير المتوافقة فغالباً ما تعاني من عدم الاستقرار أو تتطلب ضبطاً دقيقاً لمعلمات الجزاء.
قوة المعلمات (Parameter robustness): المسألة هي مسألة اضطراب منفرد بالنسبة لـ ι. عندما تقترب ι من الصفر، يؤول النموذج إلى المرونة الكلاسيكية، لكن شرط ديريكليه من الرتبة العليا (∂nu=0) لا يمكن تحقيقه بواسطة الحل الكلاسيكي، مما يؤدي إلى طبقة حدودية فيزيائية بسمك O(ι). وتتعرض الطرق القائمة على الإزاحة، التي تفرض هذا الشرط بقوة على الفضاء المنفصل، لظهور "طبقة حدودية عددية" بسمك محكوم بحجم الشبكة h، مما يؤدي إلى تدهور معدلات التقارب إلى O(h1/2) عندما تكون ι≪1. بالإضافة إلى ذلك، عندما تؤول λ→∞ (حد عدم الانضغاط)، قد تعاني الطرق من قفل الحجم (volumetric locking).
إعادة صياغة من الدرجة الأولى: تمت إعادة صياغة معادلة الدرجة الرابعة كنظام من الدرجة الأولى عبر إدخال الإجهاد الكلي σ ومتغيرين مساعدين: إجهاد كوشي المقياس σ~ والإجهاد الفائق المقياس χ. ينتج عن ذلك نظام مكون من أربعة حقول يشمل الإزاحة u، والإجهاد الكلي σ، وإجهاد كوشي المقياس σ~، والإجهاد الفائق المقياس χ. ويتم فرض التماثل بقوة على جميع موترات الإجهاد الثلاثة.
الفضاءات المنفصلة:
يتم تقطيع الإزاحة u والإجهاد الكلي σ باستخدام فضاءات الدي جيه المتدرجة (SDG). تستخدم هذه الفضاءات تقسيم "ألفيلد" (Alfeld split) للشبكة الأولية لإنشاء شبكة ثنائية، مما يفرض الاستمرارية عبر الحواف الأولية لـ u والاستمرارية العمودية عبر الحواف الثنائية للإجهاد الكلي σ.
يتم تقطيع الإجهادات المقياسة σ~ وχ باستخدام أزواج رافير-توماس (RT) وفضاءات "دي جيه" (DG)، مع فرض التماثل.
الهجنة (Hybridization): لتقليل حجم النظام العالمي، يتم تخفيف الاستمرارية العمودية بين العناصر للإجهاد الكلي والإجهاد الفائق المقياس، وإعادة فرضها بشكل ضعيف عبر مضاعفات لاغرانج المعرفة على هيكل الشبكة (الحواف الثنائية). يسمح هذا بـ التكثيف الاستاتيكي المحلي، حيث يتم استبعاد درجات الحرية الداخلية محلياً، مما يترك فقط المضاعفات كمتغيرات مقترنة عالمياً.
الشروط الحدودية: يتم دمج شرط ديريكليه من الرتبة العليا بشكل طبيعي في الصيغة التباينية بدلاً من فرضه بقوة على الفضاء المنفصل، مما يتجنب الطبقة الحدودية العددية.
المساهمات الرئيسية والنتائج النظرية تقدم الورقة تحليلاً نظرياً صارماً يدعم الطريقة المقترحة:
الاستقرار: يثبت المؤلفون وجود شرط inf-sup على الفضاء المتماثل ومتراجحة كورن المنفصلة. تُستخدم هذه الأدوات لإثبات الاستقرار الجبري للمخطط.
التقارب:
تحت انتظام كافٍ، تحقق الطريقة معدلات تقارب مثالية للإزاحة والإجهادات.
والأهم من ذلك، أن تقديرات الخطأ هي موحدة في كل من ι وλ. تحتفظ الطريقة بمعدلات التقارب المثالية حتى عندما تكون ι≪1، متجنبة التدهور الذي تظهره الطرق القائمة على الإزاحة.
يتتبع التحليل صراحةً الاعتماد على انتظام الحل والمعلمات.
الهجنة: تطور الورقة المخطط الهجين، وتثبت أن المضاعفين (الذين يقربان آثار u وισ~) يحققان شروط inf-sup الخاصة بهما، مما يضمن استقرار النظام العالمي المختزل.
حد عدم الانضغاط: تم تقديم تعديل لحالة λ→∞ (حد عدم الانضغاط)، يتضمن تفكيك الإجهاد للتعامل مع قيد التقييس العالمي المطلوب لضمان جودة الحل.
النتائج العددية تم تسجيل تجربتين عدديتين للتحقق من النظرية:
الحل السلس: تؤكد حالة اختبار ذات حل سلس وخالٍ من التباعد تحقيق معدلات تقارب مثالية في كل من معايير L2 وH1 المنفصلة للإزاحة، ومعدلات مثالية للإجهاد الكلي. وتظهر النتائج قوة عبر نطاق واسع من λ (من 1 إلى 108) وι (من 1 إلى 10−8).
الطبقة الحدودية: تؤكد حالة اختبار ذات طبقة حدودية قوية (حيث يختلف الحل الدقيق لمسألة SGE عن الحل الكلاسيكي بالقرب من الحدود) أن الطريقة تحافظ على معدلات التقارب المثالية وقوة المعلمات حتى عندما تكون ι صغيرة جداً (10−6,10−8).
الأهمية تزعم الورقة أنها تقدم أول مخطط مختلط رباعي الحقول لمرونة تدرج الانفعال مدعوم بتحليل صارم. وتكمن أهميتها في:
تجنب الطبقة الحدودية العددية: من خلال دمج شرط الحدود من الرتبة العليا تباينياً بدلاً من فرضها بقوة، تحافظ الطة على معدلات تقارب مثالية لقيم ι الصغيرة، وهي مشكلة معروفة في الطرق القائمة على الإزاحة.
قوة المعلمات: المخطط مستقر ومتقارب بشكل موحد في كل من معلمة طول المقياس ι وثابت لاميه λ، مما يعالج مشكلات الاضطراب المنفرد وعدم الانضغاط في آن واحد.
التقريب المباشر للإجهاد: بخلاف الطرق القائمة على الإزاحة، يوفر هذا التشكيل المختلط تقريباً مباشراً للإجهاد الكلي، وهو كمية ذات أهمية عملية، دون الحاجة إلى معالجة لاحقة للمشتقات من الرتب العليا.
الكفاءة الحسابية: تقلل النسخة الهجينة من النظام العالمي إلى اثنين فقط من المضاعفات على هيكل الشبكة، مما يقلل التكلفة الحسابية بشكل كبير مقارنة بنظام الحقول الأربعة غير الهجين.
العمل حالياً مقتصر على بُعدين وشروط حدودية متماثلة (مشدودة من الطرفين)، مع الإشارة إلى أن التوسع إلى ثلاثة أبعاد وشروط حدودية عامة يعد من الأعمال المستقبلية.