Caveats on formulating finite elasto-plasticity in curvilinear coordinates

إليك ملخص تقني مفصل للورقة البحثية بعنوان "تحذيرات حول صياغة اللدونة المرنة المحدودة في الإحداثيات المنحنية" (Caveats on formulating finite elasto-plasticity in curvilinears coordinates) من إعداد "بريت" وآخرون.

1. بيان المشكلة

تتناول الورقة فجوة حرجة في التنفيذ العددي لـ اللدونة المرنة المحدودة (finite-strain elasto-plasticity) ضمن أنظمة الإحداثيات المنحنية (وتحديداً التماثل المحوري). بينما تعتمد ميكانيكا الأوساط المتصلة بشكل أساسي على تحليل الموتر (tensor analysis)، وهو ثابت الإطار، فإن تطبيقات العناصر المحدودة (FE) تعتمد على التمثيلات المصفوفية في قواعد اختيارية للمستخدم.

• المشكلة الجوهرية: عند الانتقال من القواعد الكارتيزية (المتعامدة) إلى القواعد المنحنية (غير المتعامدة، مثل الإحداثيات الأسطوانية)، غالباً ما تفشل الصيغ القياسية أو تنتج أخطاءً. وتحديداً، تؤدي التوسعات الساذجة لتعريفات الكارتيزي لـ تدرج التشوه (deformation gradient)، والجاكوبي (Jacobian)، ومقاييس الإجهاد (stress measures) إلى عدم دقة في تحول الحجم والوصف الكينماتيكي.
• التحدي: في المسائل ذات التماثل المحوري (الشائعة في الهندسة الجيوتقنية، وتشكيل المعادن، والميكانيكا الحيوية)، تتفاقم هذه الأخطاء بسبب التشوهات الكبيرة وإدراج التأثيرات غير المرنة (اللدونة). تسلط الورقة الضوء على أن إجهاد كوشي (Cauchy stress) في اللدونة ذات الانفعال المحدود يعتمد على تغيرات التكوين بما يتجاوز الحالة المرنة الحالية، مما يتطلب تمييزاً دقيقاً بين المكونات المرنة واللدنة للجاكوبي وتدرج التشوه.
• الهدف: توفير إجراء صارم وخطوة بخطوة لصياغة تحليل العناصر المحدودة المتماثل محورياً يتعامل بشكل صحيح مع هذه التعقيدات الهندسية دون اللجوء إلى نظرية المانيفولد (manifold theory) المجردة للريمان، مع التركيز بدلاً من ذلك على التنفيذ العملي ضمن الأطر الكارتيزية القياسية باستخدام تغيير صريح للقواعد.

2. المنهجية

استخدم المؤلفون مسألة توسع تجويف (توسع أسطوانة مجوفة) كنموذج مرجعي للكشف عن التناقضات الكينماتيكية ومعالجتها. تتضمن المنهجية ما يلي:

• تحليل الموتر في القواعد المختلطة: تعمل الورقة ضمن التمثيلات الكارتيزية القياسية، ولكنها تشتق صراحة قواعد التحويل إلى القواعد المنحنية (الأسطوانية) المتوافقة (covariant) والمتغايرة (contravariant). وهي تتجنب استخدام مانيفولدات الهندسة التفاضلية لصالح تغيير صريح في قواعد الأساس.
• التعريفات الكينماتيكية:
  • تدرج التشوه ($F$): مُعرّف بصرامة كمتجه ثنائي النقاط (two-point tensor) يربط بين المنحنيات المادية. توضح الورقة أن $F$ ليس مجرد تدرج لمجال متجه في الإحداثيات المنحنية، ويتطلب موتر الإزاحة (shifter tensor - $S$) للربط بين قواعد مختلفة من التكوينات.
  • الجاكوبي ($J$): تبين الورقة أن نسبة الحجم تعتمد على محدد معاملات المترية. اشتق المؤلفون العلاقة $J = \det(F) \sqrt{\det(g)}/\sqrt{\det(G)}$، مما يثبت أن $\det(F_{curvilinear}) \neq J$ بشكل مباشر.
  • التفكيك المرن-اللدن: يتم تفكيك تدرج التشوه ضربياً ($F = F^e F^p$). تعالج الورقة صعوبة تعريف مترية التكوين الوسيط الخالي من الإجهاد. وتقترح خطية عملية للجاكوبي المرن ($J_e$) باستخدام أثر الانفعال اللوغاريتمي المرن ($\text{tr}(\epsilon^e)$)، مما يبسط الخطية المتسقة.
• مقاييس الإجهاد: تقدم الورقة إجهاد إسليف-زيتا ($\zeta = J_e \sigma$) كمقياس إجهاد مفضل لللدونة ذات الانفعال المحدود مع قواعد التدفق المرتبطة، مما يضمن الاتساق الديناميكي الحراري والحفاظ على التكوين الوسيط الخالي من الإجهاد.
• الأشكال الضعيفة والتقطيع (Discretization): تم تحويل الأشكال القوية لتوازن الزخم إلى أشكال ضعيفة لكل من صياغات "لاغرانج الكلي" (Total Lagrangian - TL) و"لاغرانج المحدث" (Updated Lagrangian - UL). ثم تم اختزالها إلى صيغ ثنائية الأبعاد متماثلة محورياً عبر تكامل الإحداثي المحيطي، مع التعامل بعناية مع حدود الحجم والسطح $2\pi r$.
• الخطية (Linearization): تم اشتقاق خطية "نيوتن-رافسون" متسقة، تأخذ في الاعتبار صراحة المشتقات المتوافقة ورموز كريستوفل (Christoffel symbols)، والتي تكون غير صفرية في الإحداثيات المنحنية.

3. المساهمات الرئيسية

1. توضيح موتر الإزاحة (Shifter Tensor): تُعرف الورقة صراحة دور موتر الإزاحة في رسم خرائط مكونات الإزاحة بين القواعد المنحنية المرجعية والحالية، مما يصحح مصدراً شائعاً للخطأ حيث تُعامل الإزاحات بسذاجة كفروق متجهات بسيطة.
2. صياغة الجاكوبي الصحيحة: توضح الورقة أن محدد تدرج التشوه في المكونات المنحنية لا يساوي نسبة الحجم ($J$) ما لم يتم تصحيحه بمحددات المترية. وتوفر الصيغة الصحيحة لـ $J$ وتفكيكه المرن-اللدن.
3. استراتيجية الخطية العملية: يقترح المؤلفون خطية محددة للجاكوبي المرن ($J_e = \exp(\text{tr}(\epsilon^e))$) تقلل من المصطلحات الإضافية وتتكامل بسلاسة مع البرامج الفرعية التكوينية القياسية، مما يتجنب الحاجة إلى حساب المشتقات المعقدة لتدرج التشوه اللدن.
4. سير العمل الخوارزمي: توفر الورقة خوارزميات مفصلة (الخوارزمية 1 والخوارزمية 2) لكل من تنفيذات العناصر المحدودة (UL و TL) المتماثلة محورياً، مع توضيح المواضع التي يجب فيها إدراج المصطلحات المنحنية (رموز كريستوفل، تحديثات المترية) مقارنة بكود كارتيزي قياسي.
5. التحقق: تتحقق الورقة من الصياغة عبر مقارنة كود العناصر المحدودة ثنائي الأبعاد المتماثل محورياً مع تنفيذ كارتيزي ثلاثي الأبعاد كامل لنفس المسألة.

4. النتائج

تحقق المؤلفون من صياغتهم باستخدام نموذج أسطوانة سميكة الجدار تعرضت لتوسع شعاعي داخلي مع إزاحة تتغير خطياً على طول الارتفاع.

• المقارنة: تمت مقارنة نتائج الصياغة المقترحة (ثنائية الأبعاد متماثلة محورياً) مع تحليل عناصر محدودة ثلاثي الأبعاد باستخدام عناصر هيكساهدري (hexahedral) ذات 20 عقدة.
• الدقة: أظهرت المقارنة عدم وجود فرق ملموس بين النتائج ثلاثية الأبعاد وثنائية الأبعاد المتماثلة محورياً بالنسبة للمتغيرات الرئيسية:
  • ضغط كوشي السالب.
  • إجهاد كوشي الانحرافي.
  • الجاكوبي المرن واللدن.
• التقارب: أظهر مخطط التكرار "نيوتن-رافسون" تقارباً قوياً (عادةً 3-4 تكرارات لكل خطوة زمنية) مع معدلات تقارب متطابقة لكل من الصياغات الثلاثية والثنائية الأبعاد.
• الكفاءة الحسابية: حققت الصياغة المتماثلة محورياً دقة تشبه الدقة ثلاثية الأبعاد بتكلفة حسابية أقل بكثير (تقطيع ثنائي الأبعاد مقابل ثلاثي الأبعاد)، مما يؤكد جدوى الطريقة المقترحة.

5. الأهمية

تعتبر هذه الورقة مهمة لمجتمع الميكانيكا الحسابية لعدة أسباب:

• منع الأخطاء: تحدد وتصحح الأخطاء الجوهرية والدقيقة التي تظهر عند تنفيذ اللدونة ذات الانفعال المحدود في الإحداثيات المنحنية. يمكن للتعريفات غير الصحيحة للكميات الكينماتيكية أن تؤدي إلى مصفوفات صلابة مماسية غير متسقة، أو عدم التقارب، أو حالات إجهاد غير صحيحة فيزيائياً.
• الربط بين النظرية والتطبيق: من خلال تجنب نظرية المانيفولد المجردة والتركيز على تحولات القواعد الصريحة، توفر الورقة نهجاً "كتيبياً" (cookbook approach) يسمح لمطوري العناصر المحدودة بتنفيذ حلول متماثلة محورياً قوية دون الحاجة إلى خبرة عميقة في الهندسة التفاضلية.
• تمكين المحاكاة المعقدة: يتيح التعامل الصارم مع الجاكوبي ومقاييس الإجهاد المحاكاة الموثوقة للمسائل المتماثلة محورياً المعقدة التي تتضمن تشوهات لدنة كبيرة، مثل اختراق التربة (CPT)، وحفر الأنفاق، وتشكيل المعادن، حيث تكون تغيرات الحجم والتأثيرات متباينة الخواص أمراً حاسماً.
• التقنين: تضع الورقة نظاماً وإجراءً واضحاً ومتسقاً للتعامل مع "المزاح" (shifter)، ومعاملات المترية، والمشتقات المتوافقة في أكواد العناصر المحددة المتماثلة محورياً، مما يعزز قابلية إعادة الإنتاج والدقة في الأبحاث المستقبلية.

باختصار، توفر الورقة التصحيحات الرياضية والخوارزمية اللازمة لضمان أن تحليلات العناصر المحددة للمسائل اللدنة المرنة المتماثلة محورياً تحت انفعالات كبيرة هي دقيقة رياضياً وقوية عددياً.