دليل تطبيق طريقة العناصر المحدودة في MATLAB للهندسة الميكانيكية

مقدمة في طريقة العناصر المحدودة وتطبيقاتها في الهندسة الميكانيكية

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

نظرية العناصر الرباعية رباعية العقد في تحليل الإجهاد المستوي

العنصر الرباعي رباعي العقد (Q4) هو عنصر ثنائي الأبعاد يستخدم لتمثيل الأشكال الهندسية غير المنتظمة. لكل عقدة درجتا حرية (إزاحة في اتجاه x وy)، مما يعطي المجموع 8 درجات حرية للعنصر. مصفوفة الصلابة للعنصر تُحسب باستخدام التكامل العددي (مثل طريقة Gauss-Legendre) على منطقة العنصر. المعادلة الأساسية هي: K = ∫ B^T D B t dA، حيث B هي مصفوفة الإجهاد-الإزاحة، D هي مصفوفة الخواص المادية، وt هي سمك الصفيحة.

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

لنبدأ بكتابة دالة MATLAB تحسب مصفوفة الصلابة لعنصر رباعي رباعي العقد. سنفترض أن معامل المرونة E = 40 GPa، ونسبة بواسون ν = 0.3، وسمك الصفيحة t = 2 mm. سنستخدم التكامل Gaussian بنقطتين في كل اتجاه (2×2). إليك الكود الأساسي:

function K = elementStiffnessQ4(E, nu, t, coords)
% coords: 4x2 matrix [x1 y1; x2 y2; x3 y3; x4 y4]
% Gauss points and weights
gp = [-1/sqrt(3), 1/sqrt(3)];
w = [1, 1];
K = zeros(8,8);
for i = 1:2
  for j = 1:2
    xi = gp(i); eta = gp(j);
    % Shape functions and derivatives
    N = 0.25*[(1-xi)*(1-eta); (1+xi)*(1-eta); (1+xi)*(1+eta); (1-xi)*(1+eta)];
    dN_dxi = 0.25*[-(1-eta), (1-eta), (1+eta), -(1+eta);
                   -(1-xi), -(1+xi), (1+xi), (1-xi)];
    % Jacobian matrix
    J = dN_dxi * coords;
    detJ = det(J);
    invJ = inv(J);
    dN_dx = invJ * dN_dxi;
    % B matrix (3x8)
    B = zeros(3,8);
    for k = 1:4
      B(1,2*k-1) = dN_dx(1,k);
      B(2,2*k) = dN_dx(2,k);
      B(3,2*k-1) = dN_dx(2,k);
      B(3,2*k) = dN_dx(1,k);
    end
    % Material matrix D (plane stress)
    D = E/(1-nu^2)*[1, nu, 0; nu, 1, 0; 0, 0, (1-nu)/2];
    K = K + B' * D * B * detJ * t * w(i) * w(j);
  end
end
end

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

تحليل هيكل صفيحة رقيقة باستخدام عنصرين وأربعة عناصر

بعد بناء دالة العنصر، ننتقل إلى تجميع المصفوفة الكلية للهيكل. في هذا المثال، سنحلل صفيحة مستطيلة مقسمة إلى عنصرين أو أربعة عناصر. سنطبق قيوداً عند الحافة السفلية (إزاحة صفر) وحمولة موزعة على الحافة العلوية. نستخدم برنامج MATLAB لحل النظام الخطي K * u = F والحصول على الإزاحات. ثم نحسب الانفعالات باستخدام العلاقة ε = B * u. النتائج الأولية تظهر أن العناصر الأربعة تعطي دقة أعلى، خاصة في مناطق تركيز الإجهاد.

دراسة تأثير تغيير الخواص المادية وظروف التحميل

لإكمال التحليل، نقوم بتغيير معامل المرونة E إلى ثلاث قيم مختلفة (مثل 30 GPa، 50 GPa، 70 GPa) ونكرر الحسابات. نلاحظ أن الإزاحة تتناسب عكسياً مع E. كما نغير زاوية التحميل من 90° إلى 60° و45° و30°. هذه التغييرات تؤثر على توزيع الإجهاد وتظهر أهمية فهم اتجاه التحميل في التصميم. يمكن مقارنة النتائج مع برنامج تجاري مثل ANSYS للتحقق من صحة البرنامج.

التحقق من صحة البرنامج باستخدام ANSYS Mechanical APDL

للتأكد من دقة برنامجنا، نقوم بنمذجة نفس الهيكل في ANSYS باستخدام عناصر PLANE182 (عناصر رباعية رباعية العقد). بمقارنة الإزاحات القصوى، نجد أن الفرق أقل من 2%، مما يؤكد صحة الخوارزمية. هذا التمرين يعزز فهم الطالب لكيفية عمل برامج العناصر المحدودة التجارية وأهمية التحقق اليدوي.

نصائح للعرض التقديمي (PowerPoint) والمناقشة الشفوية

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

خاتمة

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