الجبر الخطي وتطبيقاته في الهندسة

الجبر الخطي: لغة الآلات الحديثة

تخيّل ذراعاً روبوتية بستة مفاصل تلحم هيكل سيارة. كل مفصل يدور بزاوية معينة، والنتيجة هي موضع واتجاه رأس اللحام في الفضاء ثلاثي الأبعاد. كيف يحسب الحاسوب الزوايا الست اللازمة لوضع رأس اللحام في نقطة محددة؟ الجواب: الجبر الخطي (Linear Algebra) — المصفوفات والمتجهات والتحويلات الخطية.

هذا ليس رياضيات نظرية — هذا هو الأساس الذي تعمل به أنظمة التحكم والروبوتات وتحليل الإجهاد في كل مصنع حديث.

المتجهات: الاتجاه والمقدار معاً

المتجه (Vector) هو كمية لها مقدار واتجاه. القوة متجه (10 نيوتن باتجاه اليمين)، السرعة متجه، الإزاحة متجه.

في الفضاء ثلاثي الأبعاد، كل متجه يُكتب بثلاثة مركبات:

v = [vx, vy, vz]

عمليات أساسية:

• الجمع: u + v = [ux+vx, uy+vy, uz+vz] — لجمع القوى المؤثرة على قطعة
• الجداء النقطي: u · v = ux×vx + uy×vy + uz×vz — يُعطي العمل المبذول بواسطة قوة
• الجداء المتقاطع: u × v — يُعطي متجهاً عمودياً على المستوى، يُستخدم في حساب عزم الدوران

تطبيق صناعي: في تحليل القوى على عنصر إنشائي (كمرة أو عمود)، كل قوة تُمثَّل بمتجه. مجموع المتجهات يُعطي القوة المحصّلة — إذا كانت صفراً، الجسم في توازن استاتيكي.

المصفوفات: تنظيم البيانات والتحويلات

المصفوفة (Matrix) هي جدول من الأرقام مُرتّب في صفوف وأعمدة. في الهندسة تُستخدم لأمرين رئيسيين:

تمثيل أنظمة المعادلات

تصوّر دائرة كهربائية بثلاث حلقات. قوانين كيرشهوف تُعطي:

10×I₁ + 5×I₂ + 0×I₃ = 24
5×I₁ - 8×I₂ + 3×I₃ = 0
0×I₁ + 3×I₂ - 6×I₃ = -12

بصيغة المصفوفات: Ax = b

A = | 10   5   0 |    x = | I₁ |    b = |  24 |
    |  5  -8   3 |        | I₂ |        |   0 |
    |  0   3  -6 |        | I₃ |        | -12 |

الحل: x = A⁻¹ × b — مصفوفة واحدة معكوسة تُعطي كل التيارات دفعة واحدة.

التحويلات الخطية

المصفوفات تُمثّل أيضاً تحويلات هندسية: الدوران، التحجيم، الانعكاس، القص.

مصفوفة الدوران في المستوى بزاوية θ:

R(θ) = | cos(θ)  -sin(θ) |
       | sin(θ)   cos(θ) |

تصوّر قطعة معدنية على طاولة ماكينة CNC. لتدويرها 45 درجة حول المحور z، نضرب كل إحداثية في مصفوفة الدوران. هذا بالضبط ما يفعله متحكم CNC آلاف المرات في الثانية.

حل أنظمة المعادلات الخطية

طريقة غاوس (Gaussian Elimination)

أكثر الطرق شيوعاً. تحوّل المصفوفة المُوسَّعة [A|b] إلى شكل مثلثي علوي باستخدام عمليات الصفوف الأولية:

1. استبدال صف بمضاعف صف آخر مطروح منه
2. تبديل صفين
3. ضرب صف بثابت غير صفري

ثم نحل بالتعويض الخلفي (Back Substitution) من أسفل لأعلى.

طريقة LU Decomposition

تُفكّك المصفوفة A إلى حاصل ضرب مصفوفتين:

A = L × U

حيث L = مثلثية سفلى، U = مثلثية عليا. الميزة: إذا تغيّر طرف المعادلة b (لكن A بقيت نفسها)، نحل مباشرة بدون إعادة التفكيك. هذا شائع في المحاكاة الصناعية حيث نحل نفس النظام بعشرات قيم الحمل المختلفة.

الطرق التكرارية (Iterative Methods)

لأنظمة ضخمة (آلاف أو ملايين المعادلات) — مثل تحليل العناصر المحدودة (FEA):

• طريقة يعقوبي (Jacobi): تُحدّث كل متغير بناءً على القيم السابقة
• طريقة غاوس-زايدل: تستخدم القيم المُحدَّثة فوراً — تتقارب أسرع

القيم الذاتية والمتجهات الذاتية

تصوّر عموداً فولاذياً تحت حمل ضغط. عند حمل معين يبدأ بالانبعاج (Buckling). هذا الحمل الحرج هو قيمة ذاتية (Eigenvalue) لمعادلة الاستقرار، وشكل الانبعاج هو المتجه الذاتي (Eigenvector).

التعريف الرياضي:

A × v = λ × v

المصفوفة A تُحوّل المتجه v إلى مضاعف لنفسه — لا يتغير اتجاهه، فقط مقداره يتغيّر بمعامل λ.

إيجاد القيم الذاتية:

det(A - λI) = 0

هذه المعادلة المميزة (Characteristic Equation) تُعطي متعددة حدود درجتها n (حجم المصفوفة). جذورها هي القيم الذاتية.

مثال 2×2:

A = | 4  1 |
    | 2  3 |

det(A - λI) = (4-λ)(3-λ) - 2 = λ² - 7λ + 10 = 0
λ₁ = 5,  λ₂ = 2

تطبيقات القيم الذاتية في الصناعة

تطبيق: تحليل الإجهاد (Stress Analysis)

في تحليل إجهاد عنصر ميكانيكي، حالة الإجهاد عند نقطة تُمثَّل بـ مصفوفة الإجهاد:

σ = | σxx  τxy  τxz |
    | τxy  σyy  τyz |
    | τxz  τyz  σzz |

هذه مصفوفة متماثلة 3×3. قيمها الذاتية هي الإجهادات الرئيسية (Principal Stresses) — وهي أكبر وأصغر إجهاد يتعرض له العنصر. المتجهات الذاتية تُعطي الاتجاهات الرئيسية.

لماذا هذا مهم؟ لأن معايير الانهيار (مثل معيار فون ميزس) تعتمد على الإجهادات الرئيسية:

σ_VM = √(0.5 × ((σ₁-σ₂)² + (σ₂-σ₃)² + (σ₃-σ₁)²))

إذا تجاوز σ_VM إجهاد الخضوع، القطعة ستنهار.

تطبيق: أنظمة التحكم

نظام تحكم خطي يُوصف بمعادلة الحالة:

dx/dt = A × x + B × u
y = C × x

حيث x = متجه الحالة، u = المدخلات، y = المخرجات.

استقرار النظام يعتمد كلياً على القيم الذاتية للمصفوفة A:

• إذا كانت كل القيم الذاتية لها جزء حقيقي سالب → النظام مستقر
• إذا كانت أي قيمة ذاتية لها جزء حقيقي موجب → النظام غير مستقر

مهندس التحكم يُصمّم متحكماً يُغيّر القيم الذاتية (يُسمى Pole Placement) لضمان الاستقرار.

تطبيق: حركيات الروبوت (Robotics Kinematics)

ذراع روبوتية بـ n مفصل. كل مفصل يُوصف بمصفوفة تحويل 4×4 (تحتوي دوران + إزاحة). الموضع النهائي لرأس الذراع:

T_total = T₁ × T₂ × T₃ × ... × Tₙ

حيث كل Tᵢ هي مصفوفة التحويل المتجانسة (Homogeneous Transformation Matrix):

Tᵢ = | R₃ₓ₃  d₃ₓ₁ |
     | 0₁ₓ₃    1  |

R = مصفوفة دوران 3×3، d = متجه إزاحة.

الحركية الأمامية (Forward Kinematics): معرفة الزوايا → حساب موضع الرأس (ضرب مصفوفات). الحركية العكسية (Inverse Kinematics): معرفة الموضع المطلوب → حساب الزوايا (أصعب بكثير — يحتاج حل أنظمة معادلات غير خطية غالباً).

مصفوفة يعقوبي (Jacobian Matrix) تربط بين سرعات المفاصل وسرعة الرأس:

v = J(θ) × dθ/dt

J هي مصفوفة 6×n — أداة أساسية لتخطيط المسارات وتجنب المواضع المفردة (Singularities) حيث يفقد الروبوت درجة حرية.

ملخص عملي: متى تحتاج ماذا؟

الجبر الخطي ليس ترفاً أكاديمياً — هو الأداة الحسابية الأساسية التي تعمل بها المعدات الصناعية الحديثة. من متحكم CNC بسيط إلى ذراع روبوتية بستة محاور — كلها تُنفّذ عمليات مصفوفات في كل لحظة.