Volumetric Growth in Linear Elasticity Driven by an Optimality Criterion

ملخص تقني: النمو الحجمي المدفوع بالتحسين في إطار المرونة الخطية

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

الإعداد يعتمد على المرونة الخطية في بُعدين (مع إمكانية التوسع المباشر إلى ثلاثة أبعاد). يشغل الجسم مجالاً $\Omega$ بحدود مقسمة إلى أجزاء "ديريكليه" ($\partial_D \Omega$) و"نيومان" ($\partial_N \Omega$). يتم نمذجة النمو عبر تفكك إضافي لموتر الانفعال، حيث يكون إجمالي الانفعال هو مجموع الانفعال المرن وانفعال النمو غير المرن ($E_g$). وبناءً عليه، يتم تعديل معادلات التوازن بحيث يعمل $E_g$ كمصدر داخلي.

المشكلة الجوهرية هي تقليل تدريجي (incremental minimization) لمتغير الإزاحة $u^{(i)}$ وموتر النمو $E_g^{(i)}$ عند كل خطوة زمنية $i$. يخضع التحسين لثلاثة قيود رئيسية:

1. التوازن: يجب أن تتحقق معادلات توازن المرونة الخطية.
2. توازن الكتلة: يجب أن يتطابق تغير الحجم الكلي (أو المحلي) الناتج عن النمو مع زيادة محددة مسبقاً $\Gamma^{(i)}$.
3. عدم الارتداد (التراكم): يجب أن يكون النمو تراكمياً بحتًا، مما يعني أن زيادة موتر النمو يجب أن تكون موجبة شبه محددة ($\lambda_{\min}(E_g^{(i)} - E_g^{(i-1)}) \geq 0$).

يشفر الدالة الهدفية $\Phi$ الآلية الدافعة. وتستقصي الورقة معيارين محددين:

• تقليل العمل الخارجي: ويمثل تكيفاً مدفوعاً ميكانيكياً حيث يسعى الجسم لزيادة الصلابة.
• تقليل المحيط: ويمثل تطوراً مدفوعاً هندسياً (مبدأ الأيزوبيريمتريك/المحيط المتساوي).

يتم تضمين حد تنظيم تربيعي، يتناسب مع $|E_g^{(i)} - E_g^{(i-1)}|^2$، لمعاقبة الانحرافات الكبيرة عن التكوين السابق. هذا الحد، المستوحى من نظرية "دي جورجي" للحركات الدنيا، يضمن الاستمرارية الزمنية ويؤدي إلى بنية "تدفق متدرج" (gradient-flow) في الحد المستمر.

المنهجية
يستخدم المؤلفون طريقة العناصر المحدودة (FEM) لتقطيع المسألة المستمرة.

• التقطيع (Discretization): يتم تقريب مجال الإزاحة باستخدام دوال أساس مستمرة وخطية القطع، بينما يتم تقريب موتر النمو باستخدام دوال أساس ثابتة القطع على عناصر مثلثية.
• الاختزال: من خلال حل معادلات التوازن الخطية صراحةً لمجال الإزاحة بدلالة موتر النمو، يقوم المؤلفون بإلغاء متغيرات الإزاحة. يؤدي هذا إلى اختزال المسألة إلى مسألة تقليل مقيدة ذات أبعاد منتهية تتضمن فقط متغيرات النمو.
• التحسين: يتم حل المسألة المنفصلة الناتجة عددياً باستخدام أداة fmincon في برنامج MATLAB باستخدام خوارزمية البرمجة التربيعية المتتالية (SQP). كما يتم توفير التدرجات التحليلية لدالة الهدف والقيود لضمان التقارب.
• التقريب التحليلي: للكشف عن البنية الأساسية، يستنتج المؤلفون حلاً تحليلياً تقريبياً عن طريق تجاهل قيد القيمة الذاتية غير الخطي مؤقتاً (عدم الارتداد). يسمح هذا باستنتاج تعبير مغلق لزيادة النمو، مما يوضح أن التطور يتبع تدفقاً متدرجاً مسقطاً في فضاء الانفعالات غير المرنة.

النتائج الرئيسية
أُجريت تجارب عددية على ثلاث حالات نموذجية: عارضة مثبتة من الطرفين، وعارضة كابولية (كلاهما تحت حمل خارجي)، وجسم حر (تقليل المحيط).

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

الأهمية والادعاءات
تدعي الورقة القدرة على صياغة النمو الحجمي كعملية مدفوعة بالتحسين حيث يتم اختيار موتر النمو بواسطة مبدأ اختيار ميكانيكي بدلاً من قانون ظواهري.

• التحول النظري: تميز هذه الدراسة نفسها من خلال التعامل مع النمو ليس كمدخل حركي وصفي، بل كنتيجة لتقليل دالة محددة خاضعة للقيود الميكانيكية.
• الرؤية الهيكلية: يكشف النموذج أن التطور يمكن تفسيره كتدفق متدرج مسقط في فضاء الانفعالات غير المرنة. ويسلط المؤلفون الضوء على أن تقعر دالة الامتثال (compliance functional) أمر بالغ الأهمية؛ فعند فقدان التقعر، قد تظهر عدم الوحدانية والتموضع (localization).
• إثبات المفهوم: تعمل النظرية الخطية كإثبات للمفهوم، حيث تُظهر أن التحسين المقيد ميكانيكياً وحده كافٍ لتوليد مورفولوجيات (أشكال) معقدة ناتجة عن النمو (مثل تكيف الشكل وأنماط الإجهاد المتبقي) دون افتراض قوانين نمو بيولوجية محددة.
• الآفاق المستقبلية: يذكر المؤلفون أن هذا الإطار يمهد الطريق لتوسيع النهج ليشمل المرونة غير الخطية بالكامل والأنظمة الأكثر واقعية التي تتضمن أهدافاً طاقية متنافسة، رغم أنهم لا يدعون إلى حل هذه الأنظمة البيولوجية المعقدة في هذا العمل الحالي.

تخلص الورقة إلى أن الإطار المقترح ينجح في توحيد التوازن الميكانيكي، وتوازن الكتلة، وعدم الارتداد ضمن سياق متغير (variational setting)، مما يوفر أداة قوية للتنبؤ بالمورفولوجيا الناتجة عن النمو.