A novel mixed-variable beam element and its application on the free vibration of Timoshenko beam
تقترح هذه الورقة عنصر عناصر محدودة مختلطة المتغيرات جديداً يتكون من عقدتين و8 درجات حرية لعوارض تيموشينكوو، يدمج القوى الداخلية كدرجات حرية عقدية، مما يظهر دقة حسابية فائقة ودقة مباشرة في حساب القوى الداخلية مقارنة بالعناصر التقليدية في تحليل الاهتزاز الحر.
المؤلفون الأصليون:Lv Ding, Xue Xia, Siyuan Bao, S D Wang
في عالم الهندسة، تُعد العوارض (الكمرات) هي خيول العمل الصامتة التي تدعم الجسور وناطحات السحاب وأجنحة الطائرات. ولضمان عدم فشل هذه الهياكل، يجب على المهندسين التنبؤ بكيفية اهتزازها عند تعرضها للرياح أو حركة المرور أو الزلازل. وعلى مدى أكثر من قرن، وجهت نظريتان رئيسيتان هذه التنبؤات. تفترض النظرية الأقدم والأبسط أن المقطع العرضي للعارضة يظل مسطحاً تماماً ومتعامداً مع طولها أثناء انحنائها، متجاهلة حقيقة أن المادة نفسها يمكن أن تنزلق داخلياً أو تتعرض للقص. يعمل هذا النموذج بشكل جيد مع العوارض الطويلة والنحيفة جداً، لكنه يصبح غير دقيق بالنسبة للعوارض الأقصر والأكثر سمكاً حيث تبرز أهمية الانزلاق الداخلي. أما النظرية الأكثر تقدماً، والمعروفة باسم نظرية "تيموشينكو" للعوارض، فتصحح ذلك من خلال الإقرار بأن المقطع العرضي للعارضة يمكن أن يدور بشكل مستقل عن انحنائها، مما يأخذ في الاعتبار كلاً من انزلاق الطبقات ووزن المادة نفسها. ورغم أن هذه النظرية أكثر دقة، إلا أن تحويلها إلى نموذج حاسوبي للهياكل الواقعية كان يمثل صراعاً تاريخياً؛ إذ تتطلب الطرق الحاسوبية القياسية عدداً مفرطاً من الأجزاء الصغيرة للغاية للحصول على إجابة دقيقة، أو أنها تنتج نتائج تبتعد عن الواقع كلما زادت سرعة الاهتزاز وأصبح أكثر تعقيداً.
لقد طور فريق من الباحثين من جامعة سوزهو للعلوم والتكنولوجيا وجامعة ولاية أوكلاهوما طريقة جديدة لنمذجة هذه العوارض تتجاوز هذه القيود التقليدية. فقد ابتكروا وحدة بناء متخصصة للمحاكاة الحاسوبية، وهي عنصر رقمي يلتقط فيزياء العارضة بشكل أكثر اكتمالاً من النسخ السابقة. في النماذج الحاسوبية القياسية، يتم تعريف جزء العارضة عادةً فقط من خلال مقدار حركتها للأعلى والأسفل ومدى التوائها. ومع ذلك، يضيف العنصر الجديد معلومتين حاسمتين لكل نقطة اتصال: القوى الداخلية التي تدفع وتسحب العارضة، وتحديداً قوة القص وعزم الانحناء. ومن خلال تضمين هذه القوى مباشرة في تعريف النموذج، وجد الباحثون أن بإمكانهم وصف سلوك العارضة بدقة أكبر بكثير دون الحاجة إلى تقطيع الهيكل إلى آلاف القطع الصغيرة.
اختبر الباحثون هذا العنصر الجديد عبر محاكاة الاهتزاز الحر للعوارض تحت ثلاثة ظروف شائعة: عارضة مدعومة من كلا الطرفين، وعارضة مثبتة بإحكام من كلا الطرفين، وعارضة مثبتة من طرف واحد بينما يتدلى الطرف الآخر بحرية. وقارنوا نتائجهم بالحلول الرياضية المعروفة ومحاكاة البرامج التجارية عالية المستوى. كانت النتائج مذهلة؛ فعند حساب الترددات الطبيعية — وهي المعدلات المحددة التي تهتز بها العارضة طبيعياً — أنتج العنصر الجديد أخطاءً كانت عادةً نصف أخطاء النماذج التقليدية أو أقل. فعلى سبيل المثال، في العارضة المدعومة ببساطة، أظهرت الطريقة التقليدية أخطاءً تتجاوز الثلاثة بالمائة في أنماط الاهتزاز العليا، بينما حافظت الطريقة الجديدة على أخطاء أقل بك بنسبة كبيرة من الواحد بالمائة، حتى عند استخدام شبكة خشنة نسبياً مكونة من ستة عشر جزءاً فقط. وقد ظل هذا المستوى من الدقة ثابتاً سواء كانت العارضة منتظمة السمك أو مستدقة (متغيرة الشكل من طرف لآخر).
وبعيداً عن مجرد حساب سرعة اهتزاز العارضة، فإن الطريقة الجديدة تتفوق أيضاً في التنبؤ بالقوى الداخلية داخل الهيكل. ففي العديد من النهج التقليدية، يتعين على المهندسين إجراء حسابات إضافية، غالباً ما تكون غير دقيقة، لتحديد قوى القص والانحناء بعد حل عملية الاهتزاز الرئيسية. ولأن العنصر الجديد يعامل هذه القوى كأجزاء أساسية من النموذج منذ البداية، فإن الحاسوب يحلها مباشرة جنباً إلى جنب مع الحركة. وقد أظهرت الدراسة أن توزيع هذه القوى الداخلية طابق الحلول التحليلية الدقيقة بشكل شبه تام، حتى بالنسبة للاهتزازات المعقدة وعالية التردد. وهذا يعني أنه يمكن للمهندسين الآن رؤية مكان تركز الإجهاد الأعلى داخل العارضة بتفاصيل كان من الصعب تحقيقها سابقاً دون جهد حاسوبي هائل.
إن تداعيات هذا العمل كبيرة جداً لتصميم الهياكل الحديثة خفيفة الوزن. فمع توجه الهندسة نحو مواد أخف وأكثر تعقيداً، تتقلص هوامش الخطأ، وتزداد الحاجة إلى أدوات محاكاة عالية الدقة وفعالة. يوفر هذا العنصر الجديد وسيلة لتحقيق هذه الدقة دون التكلفة الحسابية الثقيلة المتمثلة في تكرار النموذج إلى آلاف الأجزاء الصغيرة. كما أنه يوفر معياراً موثوقاً لتقييم الطرق العددية الأخرى، ويقترح مساراً مستقبلياً حيث يمكن تحليل "التوائم الرقمية" لجسورنا وطائراتنا بدقة تماثل الواقع الفيزي الذي تمثله. لقد أثبت الباحثون أنه من خلال إعادة التفكير في المتغيرات الأساسية المستخدمة لوصف العارضة، من الممكن حل مشكلات طويلة الأمد في ديناميكا الهياكل باستخدام نهج أبسط وأكثر متانة.
ملخص تقني: عنصر عارضة هجين جديد متغير المتغيرات للاهتزاز الحر لعوارض تيموشينكو
بيان المشكلة غالبًا ما تعاني عناصر عوارض تيموشينكو التقليدية، المستخدمة على نطاق واسع في تحليل المنشآت حيث يكون تشوه القص والقصور الذاتي الدوراني مهمين، من قيود في الدقة الحسابية والكفاءة. تعتمد الصيغ التقليدية عادةً على درجتين من درجات الحرية (DOFs) لكل عقدة: الإزاحة المستعرضة والدوران المقطعي. وتحدد الورقة البحثية أن هذه العناصر القياسية غالبًا ما تفشل في مراعاة تأثيرات الشروط الحدية لقوة القص وعزم الانحناء ضمن بناء المجال المتغير بشكل كافٍ. ونتيجة لذلك، قد تتطلب غالبًا شبكات دقيقة للغاية لتحقيق دقة مقبولة، لا سيما لأنماط الاهتزاز ذات الرتب العليا، أو قد تظهر قفل القص (shear locking) أو دقة غير مثالية في توزيع القوى الداخلية. وبينما تم استكشاف استراتيجيات استيفاء متنوعة (خطية، تربيعية، تكعيبية) ومبادئ تباين مختلطة، لا تزال هناك حاجة إلى عنصر يقدم في آن واحد دقة عالية، وكفاءة، وإمكانية الوصول المباشر إلى القوى الداخلية دون التعقيد المرتبط بحل المعادلات التفاضلية المتجانسة أو استخدام تعبيرات حقول إزاحة شديدة التعقيد.
المنهجية يقترح المؤلفون عنصر عارضة تيموشينكو مبتكرًا يتكون من عقدتين و8 درجات حرية (DOF)، يُطلق عليه اسم "عنصر عارضة تيموشينكو الاندماجي" (fusion Timoshenko beam element). يكمن الابتكار الجوهري في بناء درجات الحرية العقدية وحقول الإزاحة:
درجات الحرية العقدية: على عكس العناصر التقليدية، تتضمن درجات الحرية العقدية لهذا العنصر الجديد ليس فقط الإزاحة المستعرضة (w) والدوران المقطعي (ϕ)، بل وأيضًا مشتقاتهما الأولى، والتي تقابل فيزيائيًا قوة القص (Q) وعزم الانحناء (M). وبالتالي، تمتلك كل عقدة أربع درجات حرية معممة: wi,ϕi,Qi,Mi.
بناء حقل الإزاحة: استنادًا إلى نظرية عارضة تيموشينكو، حيث تكون الإزاحة المستعرضة والدوران المقطعي مستقلين متبادلاً ولكن مرتبطين عبر تشوه القص، اشتق المؤلفون تعبيرين جديدين ومتسقين لحقل الإزاحة. تم بناؤهما باستخدام دوال استيفاء هيرميت التكعيبية (cubic Hermite interpolation functions).
يتم التعبير عن حقل الإزاحة المستعرضة كدالة للإزاحات الدورانية، والدورانات، والقوى الداخلية العقدية.
يتم صياغة حقل الدوران المقطعي بشكل مماثل لضمان الاتساق مع القوى الداخلية.
نهج الاشتقاق: تم اشتقاق صيغة العنصر باستخدام مبدأ تباين من نوع الإزاحة (مبدأ هاميلتون)، وليس مبدأ تباين مختلط. لا يتم استيفاء القوى الداخلية بشكل مستقل؛ بل يتم اشتقاقها وتحويلها من العلاقات التكوينية الدقيقة بين القوى الداخلية ومشتقات الإزاحة. وهذا ما يميز العنصر المقترح عن العناصر المختلطة أو الهجينة الكلاسيكية (مثل عناصر Taylor-Filippou أو Pian).
الصياغة: تمت صياغة كثافة طاقة الانفعال والطاقة الحركية مع مراعاة تشوهات الانحناء والقص. تم اشتقاق مصفوفة الصلابة للعنصر (Ke) ومصفوفة الكتلة (Me) عن طريق تكامل تعبيرات الطاقة هذه عبر مجال العنصر، بما يتناسب مع العوارض ذات المقطع العرضي المنتظم أو الخطي المتغير. ويتم حل النظام العالمي كمسألة قيم ذاتية قياسية لتحديد الترددات الطبيعية وأشكال الأنماط.
المساهمات الرئيسية
صيغة عنصر مبتكرة: تطوير عنصر عارضة تيموشينكو بـ 8 درجات حرية يدمج القوى الداخلية (القص والعزم) مباشرة في درجات الحرية العقدية.
استيفاء متسق: اقتراح تعبيرين متكاملين لحقول الإزاحة يعتمدان على استيفاء هيرميت التكعيبي، يعملان على تنسيق الإزاحة المستعرضة والدوران مع الوفاء بالقوانين التكوينية. la
حساب مباشر للقوى الداخلية: ميزة كبيرة تتمثل في الحصول على قوى القص وعزوم الانحناء مباشرة عند العقد أثناء عملية الحل، مما يلغي الحاجة إلى عمليات الاشتقاق اللاحقة التي غالبًا ما تسبب أخطاء في العناصر التقليدية.
الكفاءة والدقة: تهدف الطريقة إلى حل عدم الكفاءة الحسابية للعناصر التقليدية عند الأنماط العليا مع الحفاظ على دقة عالية.
النتائج العددية تحقق المؤلفون من صحة العنصر المقترح من خلال تحليلات الاهتزاز الحر لعوارض تيموشينكو تحت ثلاثة ظروف حدية: بسيطة الارتكاز (S-S)، مثبتة من الطرفين (C-C)، ومثبتة من طرف وحرة من الآخر (C-F). تمت المقارنة مع الحلول التحليلية، وعناصر تيموشينكو التقليدية، وبرمجيات العناصر المحدودة التجارية (ABAQUS).
الترددات الطبيعية:
بالنسبة للعوارض المنتظمة، أظهر العنصر المقترح أخطاء نسبية أقل بكثير مقارنة بالعناصر التقليدية. وبالنسبة للأنماط عالية الرتبة، وُجد أن خطأ العنصر المقترح هو نصف أو أقل من خطأ العنصر التقليدي ذي رتبة الاستيفاء نفسها.
في حالات S-S و C-C و C-F، حافظ العنصر المقترح على أخطاء نسبية أقل من 1% حتى مع الشبكات الخشنة (مثل 16 عنصرًا)، بينما تجاوزت العناصر التقليدية غالبًا نسبة خطأ 2-4% للأنماط العليا.
بالنسبة للعوارض ذات المقطع العرضي المتغير الخطي (العوارض المستدقة)، أظهر العنصر توافقًا ممتازًا مع نتائج ABAQUS، حيث ظلت الأخطاء النسبية عمومًا أقل من 0.2% لأول خمسة أنماط، حتى مع نسب الاستدقاق الحادة.
أشكال الأنماط والقوى الداخلية:
تطابقت أشكال أنماط الإزاحة والدوران المستخرجة من العنصر المقترح بشكل وثيق مع الحلول التحليلية، حتى مع عدد قليل من العناصر.
والأهم من ذلك، كانت توزيعات قوة القص وعزم الانحناء المستمدة من العنصر المقترح متسقة للغاية مع الملامح التحليلية. وقد وقعت القيم العقدية بدقة على الأغلفة التحليلية، مما يؤكد قدرة العنصر على التقاط القيم القصوى للقوى الداخلية ودقة الطور دون الحاجة إلى تكثيف الشبكة بشكل مفرط.
الأهمية والادعاءات تزعم الورقة أن عنصر عارضة تيموشينكو الاندماجي المقترح يوفر منهجية عددية عالية الكفاءة والدقة لتحليل اهتزاز منشآت تيموشينكو. وتكمن أهميته الأساسية في قدرته على:
تقليل الاعتماد على الشبكة: تحقيق دقة فائقة بعدد أقل من العناصر مقارنة بالصيغ التقليدية، مما يقلل التكلفة الحسابية.
تعزيز الدقة: توفير أخطاء أقل بكثير للترددات الطبيعية وتوزيعات القوى الداخلية عالية الرتبة.
المخرجات المباشرة: تقديم القوى الداخلية العقدية (القص والعزم) مباشرة كجزء من الحل، مما يضمن دقة عالية في تحليل الإجهادات دون الحاجة إلى معالجة لاحقة إضافية.
المعايرة المرجعية: يشير المؤلفون إلى أن النتائج المتوقعة بواسطة هذا العنصر يمكن أن تعمل كمراجع معيارية لتقييم أداء الطرق العددية الأخرى.
تخلص الدراسة إلى أن هذا النهج يعالج قيود العناصر الموجودة من خلال دمج شروط القوة في المجال المتغير، مما يوفر أداة قوية للتطبيقات الهندسية الحديثة التي تتطلب دقة وكفاءة عاليتين.