Interaction between a crack and punches in an orthotropic strip
تتقصى هذه الورقة عامل شدة الإجهاد لصدع غريفيث في شريط متباين الخواص يخضع لضغط من كبسات على كلا الوجهين، وذلك باستخدام تحويلات فوريه لاشتقاق معادلات تكاملية منفردة، والتي يتم حلها بعد ذلك باستخدام كثيرات حدود تشيبيشيف لتوفير تقريبات تحليلية ونتائج عددية رسومية.
في عالم الهندسة، نادراً ما تكون المواد مثالية. فحتى أقوى المعادن والمركبات تحتوي على عيوب ضئيلة وغير مرئية يمكن أن تتحول إلى كسور كارثية تحت الضغط. كان هذا الواقع مفهوماً لأول مرة منذ أكثر من قرن على يد آلان أرنولد غريفيث، وهو مهندس أدرك أن قوة المادة لا تتحدد بالروابط بين ذراتها، بل بوجود شقوق مجهرية. وقد اكتشف أن هذه الشقوق تعمل كمراكز لتركيز الإجهاد، مما يضاعف القوة المسلطة على المادة حتى تنكسر. وللتنبؤ بموعد حدوث هذا الفشل، يستخدم العلماء قيمة تسمى "عامل شدة الإجهاد". فكر في هذه القيمة كجهاز قياس يقيس مدى تراكم الإجهاد عند طرف الشق تماماً؛ فإذا ارتفع هذا القياس كثيراً، سيتوسع الشق وتنهار البنية. ويعد هذا المفهوم حيوياً لتصميم كل شيء، بدءاً من أجنحة الطائرات وصولاً إلى المركبات الفضائية، حيث تمتلك المواد غالباً خصائص مختلفة اعتماداً على الاتجاه الذي تُقاس فيه، وهي خاصية تُعرف باسم "التباين الاتجاهي" (orthotropy).
لقد تصدى فريق من الباحثين مؤخراً لمتغير معقد لهذه المشكلة: ماذا يحدث عندما يوجد شق داخل شريط من هذه المادة، وهذا الشريط يتعرض في الوقت نفسه للضغط بواسطة أدوات صلبة، أو "مكابس"، على سطحيه العلوي والسفلي؟ يحاكي هذا السيناريو مواقف من العالم الحقيقي حيث قد يكون المكون تحت تأثير شد ناتج عن حمل ما، بينما يتعرض في الوقت ذاته لضغط من الآلات المحيطة. ركز الباحثون على شريط من مادة يحتوي على شق مركزي يمتد أفقياً. وأرادوا فهم كيف سيغير التفاعل بين الشق والضغط الخارجي من المكابس "عامل شدة الإجهاد" عند طرف الشق. ولحل هذه المعضلة، استخدموا تقنيات رياضية متقدمة لترجمة المشكلة الفيزيائية إلى مجموعة من المعادلات التي تصف القوى والحركات داخل المادة. ثم طبقوا طريقة تتضمن سلاسل رياضية خاصة لإيجال حل تقريبي، مما سمح لهم بحساب عامل شدة الإجهاد لسيناريوهات مختلفة دون الحاجة إلى إجراء تجارب فيزيائية.
بحثت الدراسة تحديداً في نوعين من المواد: البيريليوم، وهو معدن خفيف الوزن يُستخدم في تطبيقات الفضاء والطاقة النووية، ومركب "الصلب-مايلر" (Steel-Mylar)، الذي يوفر قوة عالية ومقاومة للصدمات. وقد قام الباحثون بمحاكاة كيفية تغير عامل شدة الإجهاد عبر تغيير عمق الشريط ونسبة قوة فتح الشق إلى قوة الضغط بواسطة المكابس. ووجدوا أن سلوك عامل شدة الإجهاد يعتمد بشكل كبير على المادة والتوازن بين القوى. فبالنسبة للبيريليوم، عندما كانت قوة الضغط مساوية تقريباً لقوة الفتح، انخفض عامل شدة الإجهاد في البداية مع زيادة عمق الشريط، ولكنه بدأ في الارتفاع مرة أخرى. ومع ذلك، عندما كانت قوة الضغط أضعف، انخفض عامل شدة الإجهاد بسرعة ثم استقر. وفي المقابل، بالنسبة لمركب الصلب-مايلر، تسببت قوة الضغط القوية في ارتفاع عامل شدة الإجهاد في البداية قبل أن ينخفض مع زيادة عمق الشريط، بينما أدت قوى الضغط الأضعف إلى تناقص مستمر.
تقدم هذه النتائج صورة أوضح لكيفية سلوك الشقوق في البيئات المعقدة. لاحظ الباحثون أنه مع زيادة عمق الشريط، يميل عامل شدة الإجهاد إلى الاستقرار عند القيمة واحد، مما يؤكد أن المادة تتصرف ككتلة غير محدودة عندما تكون الحدود بعيدة بما يكفي. والأهم من ذلك، أظهرت الدراسة أن زيادة قوة الضغط من المكابس تقلل عموماً من عامل شدة الإجهاد، مما يجعل من الصعب فعلياً فتح الشق ونموه. وهذا يشير إلى أنه في تكوينات معينة، يمكن أن يعمل الضغط الخارجي كآلية لوقف انتشار الشق، مما قد يمنع انهيار الهيكل. وتوفر هذه النتائج إطاراً رياضياً يمكن للمهندسين استخدامه للتنبؤ بشكل أفضل بسلامة المكونات المصنوعة من هذه المواد المتخصصة، مما يضمن إدارة التوازن الدقيق بين الشد والضغط بشكل صحيح لمنع الكسور غير المتوقعة.
ملخص تقني: التفاعل بين صدع ومثاقب في شريط متباين الخواص
بيان المشكلة يتناول هذا البحث مسألة مرونة مستوية (plane elastostatic) تتضمن صدع غريفث (Griffith crack) يقع عند المستوى المتوسط لشريط متباين الخواص (orthotropic) لانهائي السمك 2h. يتعرض الشريط لتحميل حدودي عبر مثاقب (punches) على كلا الوجهين. وتحديداً، يتم فتح الصدع (المعرف بـ ∣x∣≤1,y=0) بواسطة شد سطحي عمودي −p(x)، بينما تتعرض حدود الشريط (y=±h) لمثاقب تسلط شداً سطحياً عمودياً −q(x) على المنطقة ∣x∣<1. تهدف الدراسة إلى تحديد عامل كثافة الإجهاد (SIF) من النوع الأول (Mode-I) عند أطراف الصدع، وتحليل كيفية تغير التفاعل بين الصدع المغمور والمثاقب الخارجية مع عمق الشريط ونسبة الأحمال المطبقة. يفترض التحليل ظروف الانفعال المستوي (plane strain) ويتجاهل الانفعال اللدن عند أطراف الصدع، ملتزماً بافتراض ميكانيكا الأوساط المتصلة بأن حجم منطقة العمليات (process zone) ضئيل مقارنة بالأبعاد الكلية.
المنهجية يستخدم المؤلفون تقنية تحويل فورييه لتقليص معادلات الحركة التفاضلية الحاكمة للوسط متباين الخواص إلى نظام من المعادلات التكاملية المنفردة المتزامنة.
الصياغة الرياضية: تم حل معادلات الإزاحة باستخدام تمثيلات تكاملية تتضمن دوالاً تعسفية A(s,y) و B(s,y). تم تطبيق الشروط الحدودية على أوجه الصدع وأسطح الشريط لاشتقاق العلاقات بين هذه الدوال.
الاختزال إلى معادلات تكاملية: تم تحويل المسألة إلى زوج من المعادلات التكاملية المنفردة المترابطة ذات نوع كوشي (Cauchy-type singularities). تتضمن هذه المعادلات دوالاً مجهولة f1(t) و f2(t)، والتي ترتبط بانقطاعات الإزاحة عبر الصدع.
الحل العددي: تم حل المعادلات التكاملية المنفردة باستخدام طريقة حل متسلسلة تقريبية. تم توسيع الدوال المجهولة باستخدام كثيرات حدود تشيبيشيف من النوع الأول (Tn). ومن خلال الاستفادة من خصائص التعامد لكثيرات حدود تشيبيشيف من النوع الثاني (Un)، تم اختزال المعادلات التكاملية إلى نظام جبري لانهائي من المعادلات.
حساب عامل كثافة الإجهاد (SIF): تم اشتقاق عامل كثافة الإجهاد (KI) تحليلياً من معاملات الحل المتسلسل، وتحديداً لطرف الصدع عند x=1.
المساهمات والنتائج الرئيسية يقدم البحث تعبيراً تحليلياً تقريبياً لعامل كثافة الإجهاد (SIF) ويعرض نتائج عددية لمادتين متباينتي الخواص هما: البيريليوم ومركب الفولاذ-مايلر (Steel-Mylar). استُخدمت الثوابت المرنة لهذه المواد لإنشاء تمثيلات رسومية لعامل كثافة الإجهاد المعياري (KI/p) مقابل عمق الشريط (h) لمختلف نسب قوة المثقب إلى حمل فتح الصدع (q/p).
البيريليوم: بالنسبة لنسب q/p القريبة من 1، ينخفض عامل كثافة الإجهاد المعياري في البداية ثم يزداد مع زيادة عمق الشريط h. أما بالنسبة للنسب الأقل (q/p=0.2,0.4,0.6)، فإن عامل كثافة الإجهة ينخفض بسرعة قبل أن يستقر مع زيادة h.
مركب الفولاذ-مايلر: بالنسبة لنسب q/p القريبة من 1، يزداد عامل كثافة الإجهاد المعياري في البداية ثم ينخفض مع زيادة h. أما بالنسبة للنسب الأقل، فإن عامل كثافة الإجهاد ينخفض بثبات مع زيادة h.
الاتجاه العام: في جميع الحالات، عندما يصبح عمق الشريط كبيراً بما يكفي، يقترب عامل كثافة الإجهاد المعياري (KI/p) من الواحد الصحيح. وهذا يؤكد التوقع النظري بأنه في حالة وجود صدع في وسط متباين الخواص لانهائي، فإن التأثيرات الحدودية للشريط تتلاشى، ويتقارب عامل كثافة الإجهاد مع القيمة الخاصة بصدع معزول في مجال لانهائي.
التفسير الفيزيائي: يُعزى انخفاض عامل كثافة الإجهاد مع زيادة العمق إلى إمكانية إغلاق الصدع تماماً بواسطة القيود الحدودية، مما يقلل من نزعة الكسر. وعلى العكس من ذلك، فإن زيادة KI/p مع زيادة q تُفسر بأن الضغط يعارض فتح الصدع.
الأهمية والنطاق يضع المؤلفون هذا العمل ضمن سياق ميكانيكا الكسر للمواد غير المتجانسة (anisotropic)، والتي تزداد أهميتها في تطبيقات الطيران والهندسة عالية السرعة (مثل الصواريخ وشفرات التوربينات). تستهدف الدراسة تحديداً الفجوة في الأدبيات المتعلقة بالأوساط متباينة الخواص التي تتضمن كلاً من الصدوع والمثاقب، وهي توليفة لوحظ أنها أقل استكشافاً من مسائل الصدع القياسية.
تكمن أهمية البحث في تطوير نموذج لفهم سلوك عامل كثافة الإجهاد (SIF) وإمكانية توقف الصدع أو انتشاره نتيجة التفاعل بين الصدع المغمور والمثاقب السطحية. ويذكر المؤلفون أن النموذج المقترح يهدف إلى أن يكون مفيداً للباحثين والمهندسين العاملين مع المواد متباينة الخواص ومتجانسة الخواص. وتختتم الورقة بتواضع، مشيرة إلى أن الدراسات المستقبلية ستوسع هذا العمل ليشمل الصدوع الواقعة على خط واحد (collinear cracks) في المواد متباينة الخواص والموضوعة بين شرائط مرنة مختلفة تحت ظروف تحميل متنوعة (شد مفاجئ، تحميل نقطي، أو تحميل حراري-ميكانيكي).