Revisiting the Stress Field Inside an Elastic Sphere Subjected to a Concentrated Load

ملخص تقني: إعادة دراسة مجال الإجهاد داخل كرة مرنة خاضعة لحمل مركز

صياغة المسألة
تتناول هذه الدراسة مسألة القيمة الحدية لتحديد مجالات الإجهاد والإزاحة الكاملة داخل جسم صلب متجانس، متماثل المناحي، ومرن خطياً، على شكل كرة نصف قطرها $R$، خاضعة لحمل مركز عمودي على سطحها. وبينما توجد حلول كلاسيكية للمجالات اللانهائية أو شبه اللانهائية (مثل حلول كيلفين وبوسينيسك)، فإن الحلول التحليلية الدقيقة للأجسام ثلاثية الأبعاد ذات الحدود المنحنية تظل نادرة. يتضمن التكوين المحدد إجهاداً ضاغطاً مركزاً $\sigma_0$ يُطبق عند نقطة على السطح (القطب الشمالي في البداية). ويشير المؤلفون إلى أن مثل هذا الحمل المنفرد ينتهك تقنياً توازن القوى العالمي، مما يؤدي إلى حركة الجسم الجاسئ؛ ومع ذلك، يركز التحليل حصرياً على توزيع الإجهاد الداخلي، مع معاملة انتقال الجسم الجاسئ كنمط منفصل $n=1$ لا يؤثر على التشوه المرن الداخلي.

المنهجية
تمت صياغة التحليل ضمن إطار الديناميكا المرنة الخطية ثلاثية الأبعاد. استخدم المؤلفون الخطوات المنهجية التالية:

1. المعادلات الحاكمة: استُخدمت معادلة نافييه-كوشي، مع جعلها غير بعدية باستخدام الطول المميز ($R$) والزمن ($R/v_L$، حيث $v_L$ هي سرعة الموجة الطولية).
2. تمثيل الجهد: تم تفكيك مجال الإزاحة باستخدام التفكيك هلمهولتز عبر الجهود القياسية ($\phi$) والمتجهة ($\chi$). أدى هذا إلى اختزال معادلة الديناميكا المرنة المتجهة إلى معادلتين موجيتين قياسيتين تتوافقان مع الأنماط الطولية (موجات P) والمستعرضة (موجات S).
3. التوسع الطيفي: تم توسيع الحل باستخدام التوافقيات الكروية (متعددات حدود ليجاندر $P_n(\cos \theta)$) ودوال بيسل الكروية المعدلة. تم تحديد معاملات هذه التوسعات بشكل صريح من خلال فرض شروط حد الإجهاد عند السطح ($r=R$).
   ס 4. تحويل لابلاس: لحل المسألة المعتمدة على الزمن، طُبق تحويل لابلاس على الجهود. تم حل معادلات هلمهولتز المعدلة الناتجة في نطاق لابلاس، واستُخلصت المعاملات في صيغة مغلقة.
4. الحد الاستاتيكي: تم الحصول على الحل المرن الاستاتيكي بدقة كحد زمني طويل ($t \to \infty$) للصياغة الديناميكية باستخدام نظرية القيمة النهائية لتحويل لابلاس.
5. التعميم: تم تعميم الحل للحمل عند القطب ليشمل أي مواضع تحميل اختيارية $(R, \theta_0, \phi_0)$ باستخدام التحويلات الدورانية. يسمح هذا بالمعالجة المنهجية للأحمال المركزة المتعددة عبر مبدأ التراكب.
6. التقييم الديناميكي: تم تقييم تحويل لابلاس العكسي باستخدام تشويه المسار ونظرية البواقي. حدد هذا المساهمات من الأقطاب على المحور التخيلي، والتي تقابل الحل الاستاتيكي، وموجات P، وموجات S.

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

• حل تحليلي بصيغة مغلقة: تقدم الورقة تعبيرات صريحة وصيغ مغلقة لجميع مكونات موتر الإجهاد ($\sigma_{\rho\rho}, \sigma_{\theta\theta}, \sigma_{\phi\phi}, \sigma_{\rho\theta}$) ومجالات الإزاحة. وعلى عكس التمثيلات السلسلية السابقة التي غالباً ما تظل في أشكال ضمنية أو شديدة الترابط، تم تحديد جميع معاملات التوسع بشكل صريح.
• إطار موحد استاتيكي-ديناميكي: يوحد العمل بين المسألة الاستاتيكية والديناميكية، حيث يُشتق الحل الاستاتيكي كحد صارم للصياغة الديناميكية بدلاً من كونه اشتقاقاً استاتيكياً منفصلاً.
• تحليل الأنماط الصريح: يحدد التحليل صراحةً نمط التوافق الكروي $n=1$ باعتباره يمثل انتقال الجسم الجاسئ، وهو ما يمكن استبعاده عند حساب الإجهادات الداخلية. يوضح هذا نقطة غالباً ما تُعالج ضمنياً في المعالجات التحليلية الأخرى.
• فرق الإجهاد الرئيسي: تمكن التعبيرات ذات الصيغة المغلقة من الحساب المباشر للإجهادات الرئيسية وفروقها في جميع أنحاء باطن الكرة، وهي كمية بالغة الأهمية في علم الضوء المرن ثلاثي الأبعاد (photoelasticity).

النتائج

• توزيع الإجهاد الاستاتيكي: وُجد أن فرق الإجهاد الرئيسي يكون في أعلى مستوياته بالقرب من نقطة التحميل، حيث يظهر سلوكاً تباعدياً يتناسب مع $(1-\rho)^{-2}$ عند الاقتراب من السطح. على طول محور التحميل، يكون الإجهاد الشعاعي ضاغطاً، بينما تكون الإجهادات المماسية شدّية، مما يعكس التمدد الجانبي.
• انتشار الموجات الديناميكية: يلتقط الحل الديناميكي انتشار موجات P وموجات S بسرعات متميزة. تظهر النتائج ظهور:
  • موجات P أولية وموجات S تنتشر بشكل مركزي من المصدر.
  • موجات رايلي متمركزة بالقرب من السطح.
  • موجات منعكسة (موجات PP) وموجات محولة للأنماط (موجات PS) ناتجة عن تفاعلات الحدود.
  • هياكل غلاف معقدة ناتجة عن تراكب الموجات المنبعثة في أوقات مختلفة، بما يتوافق مع النظائر ثنائية الأبعاد.
• الأحمال المتعددة: تم إثبات مبدأ التراكب للحصول على تكوينات بأحمال متعددة (مثل أربعة أحمال توازن بعضها البعض)، مما يظهر مناطق تضخيم وإلغاء للإجهاد.
• التحقق: تمت مقارنة النتائج التحليلية مع محاكاة طريقة العناصر المحدودة (FEM). أظهرت المقارنة توافقاً جيداً بشكل عام، مع وجود اختلافات طفيفة بالقرب من السطح تُعزى إلى قيود دقة الشبكة في تمثيل الأحمال المركزة في FEM.

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