ديناميكيات الأجسام الصلبة

محرك بولتون ووات البخاري
يمكن وصف حركة كل مكون من مكونات محرك بولتون ووات البخاري (1784) بواسطة مجموعة من معادلات علم الحركة وعلم الحركة.

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

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

ديناميكا الأجسام الصلبة المستوية

إذا تحرك نظام من الجسيمات موازيًا لمستوى ثابت، يُقال إن النظام مقيد بحركة مستوية. في هذه الحالة، تتبسط قوانين نيوتن (الحركية) لنظام صلب مكون من N جسيمًا، Pᵢ ، حيث i = 1، ...، N ، لعدم وجود حركة في الاتجاه k . حدد القوة المحصلة وعزم الدوران عند نقطة مرجعية R ، للحصول على F=أنا=1شمالمأناأأنا،تي=أنا=1شمال(رأنا-R)×مأناأأنا،{\displaystyle \mathbf {F} =\sum _{i=1}^{N}m_{i}\mathbf {A} _{i},\quad \mathbf {T} =\sum _{i=1}^{N}(\mathbf {r} _{i}-\mathbf {R} )\times m_{i}\mathbf {A} _{i},}

حيث يشير r i إلى المسار المستوي لكل جسيم.

تُعطي حركية الجسم الصلب صيغة تسارع الجسيم P i بدلالة الموضع R والتسارع A للجسيم المرجعي، بالإضافة إلى متجه السرعة الزاوية ω ومتجه التسارع الزاوي α لنظام الجسيمات الصلب كما يلي: أأنا=α×(رأنا-R)+ω×(ω×(رأنا-R))+أ.{\displaystyle \mathbf {A} _{i}={\boldsymbol {\alpha }}\times (\mathbf {r} _{i}-\mathbf {R} )+{\boldsymbol {\omega }}\times ({\boldsymbol {\omega }}\times (\mathbf {r} _{i}-\mathbf {R} ))+\mathbf {A} .}

بالنسبة للأنظمة المقيدة بالحركة المستوية، يكون متجه السرعة الزاوية ومتجه التسارع الزاوي موجهين على طول المحور k عموديًا على مستوى الحركة، مما يبسط معادلة التسارع هذه. في هذه الحالة، يمكن تبسيط متجهات التسارع بإدخال متجهات الوحدة eᵢ من نقطة المرجع R إلى نقطة rᵢ ومتجهات الوحدة .تأنا=ك×هـأنا{\textstyle \mathbf {t} _{i}=\mathbf {k} \times \mathbf {e} _{i}}، لذا أأنا=α(Δرأناتأنا)-ω2(Δرأناهـأنا)+أ.{\displaystyle \mathbf {A} _{i}=\alpha (\Delta r_{i}\mathbf {t} _{i})-\omega ^{2}(\Delta r_{i}\mathbf {e} _{i})+\mathbf {A} .}

ينتج عن ذلك القوة المحصلة المؤثرة على النظام كما يلي: F=αأنا=1شمالمأنا(Δرأناتأنا)-ω2أنا=1شمالمأنا(Δرأناهـأنا)+(أنا=1شمالمأنا)أ،{\displaystyle \mathbf {F} =\alpha \sum _{i=1}^{N}m_{i}\left(\Delta r_{i}\mathbf {t} _{i}\right)-\omega ^{2}\sum _{i=1}^{N}m_{i}\left(\Delta r_{i}\mathbf {e} _{i}\right)+\left(\sum _{i=1}^{N}m_{i}\right)\mathbf {A} ,} وعزم الدوران تي=أنا=1شمال(مأناΔرأناهـأنا)×(α(Δرأناتأنا)-ω2(Δرأناهـأنا)+أ)=(أنا=1شمالمأناΔرأنا2)αك+(أنا=1شمالمأناΔرأناهـأنا)×أ،{\displaystyle {\begin{aligned}\mathbf {T} ={}&\sum _{i=1}^{N}(m_{i}\Delta r_{i}\mathbf {e} _{i})\times \left(\alpha (\Delta r_{i}\mathbf {t} _{i})-\omega ^{2}(\Delta r_{i}\mathbf {e} _{i})+\mathbf {A} \right)\\{}={}&\left(\sum _{i=1}^{N}m_{i}\Delta r_{i}^{2}\right)\alpha \mathbf {k} +\left(\sum _{i=1}^{N}m_{i}\Delta r_{i}\mathbf {e} _{i}\right)\times \mathbf {A} ,\end{aligned}}}

أينهـأنا×هـأنا=0{\textstyle \mathbf {e} _{i}\times \mathbf {e} _{i}=0}وهـأنا×تأنا=ك{\textstyle \mathbf {e} _{i}\times \mathbf {t} _{i}=\mathbf {k} }هو متجه الوحدة العمودي على المستوى لجميع الجسيمات P i .

استخدم مركز الكتلة C كنقطة مرجعية، وبذلك تتبسط معادلات قوانين نيوتن لتصبح F=مأ،تي=أناجαك،{\displaystyle \mathbf {F} =M\mathbf {A} ,\quad \mathbf {T} =I_{\textbf {C}}\alpha \mathbf {k} ,}

حيث M هي الكتلة الكلية و I C هي عزم القصور الذاتي حول محور عمودي على حركة النظام الصلب ويمر عبر مركز الكتلة.

جسم صلب ثلاثي الأبعاد

وصف التوجه أو الموقف

تم تطوير عدة طرق لوصف اتجاهات الجسم الصلب في ثلاثة أبعاد. ويتم تلخيصها في الأقسام التالية.

زوايا أويلر

تُنسب أول محاولة لتمثيل اتجاه ما إلى ليونارد أويلر . فقد تخيّل ثلاثة أطر مرجعية يمكن تدويرها حول بعضها البعض، وأدرك أنه بالبدء بإطار مرجعي ثابت وإجراء ثلاث دورات، يمكنه الحصول على أي إطار مرجعي آخر في الفضاء (باستخدام دورتين لتثبيت المحور الرأسي ودورة أخرى لتثبيت المحورين الآخرين). تُسمى قيم هذه الدورات الثلاث زوايا أويلر .ψ{\displaystyle \psi }يُستخدم للدلالة على التبادر،θ{\displaystyle \theta }التذبذب، وϕ{\displaystyle \phi }الدوران الذاتي.

زوايا تيت-برايان

زوايا تايت-برايان، طريقة أخرى لوصف التوجيه

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

متجه التوجيه

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

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

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

مصفوفة التوجيه

مع ظهور المصفوفات، أُعيدت صياغة نظريات أويلر. ووُصفت الدورانات بمصفوفات متعامدة تُعرف بمصفوفات الدوران أو مصفوفات جيب التمام الاتجاهي. وعند استخدامها لتمثيل اتجاه ما، تُسمى مصفوفة الدوران عادةً بمصفوفة الاتجاه أو مصفوفة الوضع.

المتجه الأويلري المذكور أعلاه هو المتجه الذاتي لمصفوفة الدوران (لمصفوفة الدوران قيمة ذاتية حقيقية وحيدة ). حاصل ضرب مصفوفتي دوران هو تركيب الدورانات. لذا، كما في السابق، يمكن تحديد الاتجاه على أنه الدوران من الإطار الأصلي للوصول إلى الإطار المطلوب وصفه.

فضاء التكوين لجسم غير متناظر في فضاء ذي n بُعد هو SO( n ) × Rⁿ . يمكن تصور التوجيه بربط أساس من متجهات المماس بالجسم. يحدد اتجاه كل متجه توجيهه.

رباعي التوجيه

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

قانون نيوتن الثاني في ثلاثة أبعاد

للنظر في ديناميكيات الأجسام الصلبة في الفضاء ثلاثي الأبعاد، يجب توسيع قانون نيوتن الثاني لتحديد العلاقة بين حركة الجسم الصلب ونظام القوى وعزوم الدوران التي تؤثر عليه.

صاغ نيوتن قانونه الثاني للجسيمات على النحو التالي: "يتناسب تغير حركة الجسم طرديًا مع القوة المؤثرة عليه، ويكون في اتجاه الخط المستقيم الذي أثرت فيه القوة." [ 3 ] ولأن نيوتن كان يشير عادةً إلى حاصل ضرب الكتلة في السرعة على أنه "حركة" الجسيم، فإن عبارة "تغير الحركة" تشير إلى حاصل ضرب الكتلة في تسارع الجسيم، ولذلك يُكتب هذا القانون عادةً على النحو التالي: F=مأ،{\displaystyle \mathbf {F} =m\mathbf {a} ,} حيث يُفترض أن F هي القوة الخارجية الوحيدة المؤثرة على الجسيم، و m هي كتلة الجسيم، و a هو متجه تسارعه. ويتم تطبيق قانون نيوتن الثاني على الأجسام الصلبة من خلال دراسة نظام صلب من الجسيمات.

نظام صلب من الجسيمات

إذا تم تجميع نظام من N جسيمًا، Pᵢ ، حيث i = 1، ...، N ، في جسم صلب، فإنه يمكن تطبيق قانون نيوتن الثاني على كل جسيم من جسيمات الجسم. إذا كانت Fᵢ هي القوة الخارجية المؤثرة على الجسيم Pᵢ ذي الكتلة mᵢ ، فإن Fأنا+ج=1شمالFأناج=مأناأأنا،أنا=1،...،شمال،{\displaystyle \mathbf {F} _{i}+\sum _{j=1}^{N}\mathbf {F} _{ij}=m_{i}\mathbf {a} _{i},\quad i=1,\ldots ,N,} حيث F ij هي القوة الداخلية للجسيم P j التي تؤثر على الجسيم P i والتي تحافظ على المسافة الثابتة بين هذين الجسيمين.

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

يُمكن تبسيط معادلات القوى هذه بشكلٍ كبير من خلال إدخال محصلة القوة وعزم الدوران المؤثرين على النظام الصلب. تُحسب محصلة القوة وعزم الدوران باختيار إحدى الجسيمات في النظام كنقطة مرجعية، R ، حيث تُطبق كل قوة خارجية مع إضافة عزم الدوران المرتبط بها. تُعطى محصلة القوة F وعزم الدوران T بالصيغ التالية: F=أنا=1شمالFأنا،تي=أنا=1شمال(Rأنا-R)×Fأنا،{\displaystyle \mathbf {F} =\sum _{i=1}^{N}\mathbf {F} _{i},\quad \mathbf {T} =\sum _{i=1}^{N}(\mathbf {R} _{i}-\mathbf {R} )\times \mathbf {F} _{i},} حيث R i هو المتجه الذي يحدد موضع الجسيم P i .

يتحد قانون نيوتن الثاني للجسيمات مع هذه الصيغ لحساب القوة المحصلة وعزم الدوران ليعطي، F=أنا=1شمالمأناأأنا،تي=أنا=1شمال(Rأنا-R)×(مأناأأنا)،{\displaystyle \mathbf {F} =\sum _{i=1}^{N}m_{i}\mathbf {a} _{i},\quad \mathbf {T} =\sum _{i=1}^{N}(\mathbf {R} _{i}-\mathbf {R} )\times (m_{i}\mathbf {a} _{i}),} حيث تتلاشى القوى الداخلية F<sub> ij</sub> في أزواج. تعطي حركية الجسم الصلب صيغة تسارع الجسيم P <sub>i</sub> بدلالة الموضع R والتسارع a للجسيم المرجعي، بالإضافة إلى متجه السرعة الزاوية ω ومتجه التسارع الزاوي α لنظام الجسيمات الصلب كما يلي: أأنا=α×(Rأنا-R)+ω×(ω×(Rأنا-R))+أ.{\displaystyle \mathbf {a} _{i}=\alpha \times (\mathbf {R} _{i}-\mathbf {R} )+\omega \times (\omega \times (\mathbf {R} _{i}-\mathbf {R} ))+\mathbf {a} .}

خصائص الكتلة

تُمثَّل خصائص كتلة الجسم الصلب بمركز كتلته ومصفوفة القصور الذاتي . اختر نقطة مرجعية R بحيث تحقق الشرط التالي: أنا=1شمالمأنا(Rأنا-R)=0،{\displaystyle \sum _{i=1}^{N}m_{i}(\mathbf {R} _{i}-\mathbf {R} )=0,}

ثم يُعرف باسم مركز كتلة النظام.

تُعرَّف مصفوفة القصور الذاتي [IR ] للنظام بالنسبة إلى نقطة المرجع R كما يلي:[أناR]=أنا=1شمالمأنا(أنا(SأناتيSأنا)-SأناSأناتي)،{\displaystyle [I_{R}]=\sum _{i=1}^{N}m_{i}\left(\mathbf {I} \left(\mathbf {S} _{i}^{\textsf {T}}\mathbf {S} _{i}\right)-\mathbf {S} _{i}\mathbf {S} _{i}^{\textsf {T}}\right),}

أينSأنا{\displaystyle \mathbf {S} _{i}}هو متجه العمود R iR ؛Sأناتي{\displaystyle \mathbf {S} _{i}^{\textsf {T}}}وهي منقولتها، وأنا{\displaystyle \mathbf {I} }هي مصفوفة الوحدة 3×3 .

SأناتيSأنا{\displaystyle \mathbf {S} _{i}^{\textsf {T}}\mathbf {S} _{i}}هو حاصل الضرب القياسي لـSأنا{\displaystyle \mathbf {S} _{i}}مع نفسه، بينماSأناSأناتي{\displaystyle \mathbf {S} _{i}\mathbf {S} _{i}^{\textsf {T}}}هو حاصل الضرب الموتري لـSأنا{\displaystyle \mathbf {S} _{i}}مع نفسه.

معادلات القوة والعزم

باستخدام مركز الكتلة ومصفوفة القصور الذاتي، تأخذ معادلات القوة وعزم الدوران لجسم صلب واحد الشكل التالي: F=مأ،تي=[أناR]α+ω×[أناR]ω،{\displaystyle \mathbf {F} =m\mathbf {a} ,\quad \mathbf {T} =[I_{R}]\alpha +\omega \times [I_{R}]\omega ,} وتُعرف باسم قانون نيوتن الثاني للحركة للأجسام الصلبة.

تُصاغ ديناميكيات نظام مترابط من الأجسام الصلبة، حيث B <sub>i</sub> , j = 1, ..., M ، عن طريق عزل كل جسم صلب وإدخال قوى التفاعل. وتؤدي محصلة القوى الخارجية وقوى التفاعل المؤثرة على كل جسم إلى معادلات القوة والعزم. Fج=مجأج،تيج=[أناR]جαج+ωج×[أناR]جωج،ج=1،...،م.{\displaystyle \mathbf {F} _{j}=m_{j}\mathbf {a} _{j},\quad \mathbf {T} _{j}=[I_{R}]_{j}\alpha _{j}+\omega _{j}\times [I_{R}]_{j}\omega _{j},\quad j=1,\ldots ,M.}

ينتج عن صياغة نيوتن 6 م معادلة تحدد ديناميكيات نظام مكون من م أجسام صلبة. [ 4 ]

الدوران في ثلاثة أبعاد

قد يُظهر الجسم الدوار، سواء كان تحت تأثير عزم الدوران أم لا، سلوكيات التبادر والترنح . المعادلة الأساسية التي تصف سلوك الجسم الصلب الدوار هي معادلة أويلر للحركة . τ=دلدت=دلدت+ω×ل=د(أناω)دت+ω×أناω=أناα+ω×أناω{\displaystyle {\boldsymbol {\tau }}={\frac {D\mathbf {L} }{Dt}}={\frac {d\mathbf {L} }{dt}}+{\boldsymbol {\omega }}\times \mathbf {L} ={\frac {d(I{\boldsymbol {\omega }})}{dt}}+{\boldsymbol {\omega }}\times {I{\boldsymbol {\omega }}}=I{\boldsymbol {\alpha }}+{\boldsymbol {\omega }}\times {I{\boldsymbol {\omega }}}} حيث أن المتجهات الزائفة τ و L هي، على التوالي، عزم الدوران على الجسم وعزمه الزاوي ، والكمية القياسية I هي عزم قصوره الذاتي ، والمتجه ω هو سرعته الزاوية، والمتجه α هو تسارعه الزاوي، و D هو التفاضل في إطار مرجعي قصوري و d هو التفاضل في إطار مرجعي نسبي ثابت مع الجسم.

تمت مناقشة حل هذه المعادلة عندما لا يكون هناك عزم دوران مطبق في مقالتي معادلة أويلر للحركة وقطع بوانسو الإهليلجي .

يستنتج من معادلة أويلر أن عزم الدوران τ المطبق عموديًا على محور الدوران، وبالتالي عموديًا على L ، ينتج عنه دوران حول محور عمودي على كل من τ و L. تُسمى هذه الحركة بالترنح . تُعطى السرعة الزاوية للترنح ΩP بالضرب الاتجاهي التالي :τ=ΩP×ل.{\displaystyle {\boldsymbol {\tau }}={\boldsymbol {\Omega }}_{\mathrm {P} }\times \mathbf {L} .}

ترنح الجيروسكوب

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

تحت تأثير عزم دوران ثابت مقداره τ ، فإن سرعة التبادر Ω P تتناسب عكسياً مع L ، وهو مقدار زخمها الزاوي: τ=ΩPلالخطيئةθ،{\displaystyle \tau ={\mathit {\Omega }}_{\mathrm {P} }L\sin \theta ,} حيث θ هي الزاوية بين المتجهين Ω P و L. بالتالي، إذا تباطأ دوران الجسم (على سبيل المثال، بسبب الاحتكاك)، فإن زخمه الزاوي يقل، وبالتالي يزداد معدل التبادر. يستمر هذا حتى يعجز الجهاز عن الدوران بسرعة كافية لتحمل وزنه، وعندها يتوقف عن التبادر ويسقط عن دعامته، ويعود ذلك في الغالب إلى أن الاحتكاك الذي يعيق التبادر يُسبب تبادرًا آخر يؤدي إلى السقوط.

بحسب الاصطلاح، فإن هذه المتجهات الثلاثة - عزم الدوران، والدوران، والترنح - كلها موجهة بالنسبة لبعضها البعض وفقًا لقاعدة اليد اليمنى .

العمل الافتراضي للقوى المؤثرة على جسم صلب

يتم الحصول على صياغة بديلة لديناميكيات الجسم الصلب والتي تحتوي على عدد من الميزات الملائمة من خلال النظر في العمل الافتراضي للقوى المؤثرة على الجسم الصلب.

The virtual work of forces acting at various points on a single rigid body can be calculated using the velocities of their point of application and the resultant force and torque. To see this, let the forces F1, F2 ... Fn act on the points R1, R2 ... Rn in a rigid body.

The trajectories of Ri, i = 1, ..., n are defined by the movement of the rigid body. The velocity of the points Ri along their trajectories are Vi=ω×(RiR)+V,{\displaystyle \mathbf {V} _{i}={\boldsymbol {\omega }}\times (\mathbf {R} _{i}-\mathbf {R} )+\mathbf {V} ,} where ω is the angular velocity vector of the body.

Virtual work

Work is computed from the dot product of each force with the displacement of its point of contact δW=i=1nFiδri.{\displaystyle \delta W=\sum _{i=1}^{n}\mathbf {F} _{i}\cdot \delta \mathbf {r} _{i}.} If the trajectory of a rigid body is defined by a set of generalized coordinatesqj, j = 1, ..., m, then the virtual displacements δri are given by δri=j=1mriqjδqj=j=1mViq˙jδqj.{\displaystyle \delta \mathbf {r} _{i}=\sum _{j=1}^{m}{\frac {\partial \mathbf {r} _{i}}{\partial q_{j}}}\delta q_{j}=\sum _{j=1}^{m}{\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}_{j}}}\delta q_{j}.} The virtual work of this system of forces acting on the body in terms of the generalized coordinates becomes δW=F1(j=1mV1q˙jδqj)++Fn(j=1mVnq˙jδqj){\displaystyle \delta W=\mathbf {F} _{1}\cdot \left(\sum _{j=1}^{m}{\frac {\partial \mathbf {V} _{1}}{\partial {\dot {q}}_{j}}}\delta q_{j}\right)+\dots +\mathbf {F} _{n}\cdot \left(\sum _{j=1}^{m}{\frac {\partial \mathbf {V} _{n}}{\partial {\dot {q}}_{j}}}\delta q_{j}\right)}

or collecting the coefficients of δqjδW=(i=1nFiViq˙1)δq1++(1=1nFiViq˙m)δqm.{\displaystyle \delta W=\left(\sum _{i=1}^{n}\mathbf {F} _{i}\cdot {\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}_{1}}}\right)\delta q_{1}+\dots +\left(\sum _{1=1}^{n}\mathbf {F} _{i}\cdot {\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}_{m}}}\right)\delta q_{m}.}

Generalized forces

For simplicity consider a trajectory of a rigid body that is specified by a single generalized coordinate q, such as a rotation angle, then the formula becomes δW=(i=1nFiViq˙)δq=(i=1nFi(ω×(RiR)+V)q˙)δq.{\displaystyle \delta W=\left(\sum _{i=1}^{n}\mathbf {F} _{i}\cdot {\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}}}\right)\delta q=\left(\sum _{i=1}^{n}\mathbf {F} _{i}\cdot {\frac {\partial ({\boldsymbol {\omega }}\times (\mathbf {R} _{i}-\mathbf {R} )+\mathbf {V} )}{\partial {\dot {q}}}}\right)\delta q.}

Introduce the resultant force F and torque T so this equation takes the form δW=(FVq˙+Tωq˙)δq.{\displaystyle \delta W=\left(\mathbf {F} \cdot {\frac {\partial \mathbf {V} }{\partial {\dot {q}}}}+\mathbf {T} \cdot {\frac {\partial {\boldsymbol {\omega }}}{\partial {\dot {q}}}}\right)\delta q.}

The quantity Q defined by Q=FVq˙+Tωq˙,{\displaystyle Q=\mathbf {F} \cdot {\frac {\partial \mathbf {V} }{\partial {\dot {q}}}}+\mathbf {T} \cdot {\frac {\partial {\boldsymbol {\omega }}}{\partial {\dot {q}}}},}

is known as the generalized force associated with the virtual displacement δq. This formula generalizes to the movement of a rigid body defined by more than one generalized coordinate, that is δW=j=1mQjδqj,{\displaystyle \delta W=\sum _{j=1}^{m}Q_{j}\delta q_{j},} where Qj=FVq˙j+Tωq˙j,j=1,,m.{\displaystyle Q_{j}=\mathbf {F} \cdot {\frac {\partial \mathbf {V} }{\partial {\dot {q}}_{j}}}+\mathbf {T} \cdot {\frac {\partial {\boldsymbol {\omega }}}{\partial {\dot {q}}_{j}}},\quad j=1,\ldots ,m.}

It is useful to note that conservative forces such as gravity and spring forces are derivable from a potential function V(q1, ..., qn), known as a potential energy. In this case the generalized forces are given by Qj=Vqj,j=1,,m.{\displaystyle Q_{j}=-{\frac {\partial V}{\partial q_{j}}},\quad j=1,\ldots ,m.}

D'Alembert's form of the principle of virtual work

The equations of motion for a mechanical system of rigid bodies can be determined using D'Alembert's form of the principle of virtual work. The principle of virtual work is used to study the static equilibrium of a system of rigid bodies; however, by introducing acceleration terms in Newton's laws this approach is generalized to define dynamic equilibrium.

Static equilibrium

يُعرَّف التوازن السكوني لنظام ميكانيكي من الأجسام الصلبة بشرط أن يكون الشغل الافتراضي للقوى المؤثرة صفرًا لأي إزاحة افتراضية للنظام. يُعرف هذا بمبدأ الشغل الافتراضي. [ 5 ] وهذا يُكافئ شرط أن تكون القوى المعممة لأي إزاحة افتراضية صفرًا، أي Q i = 0.

لنفترض أن نظامًا ميكانيكيًا يتكون من n جسمًا صلبًا، Bᵢ ، حيث i = 1، ...، n ، ولنفترض أن محصلة القوى المؤثرة على كل جسم هي أزواج القوة والعزم، Fᵢ و Tᵢ، حيث i = 1 ، ... ، n . لاحظ أن هذه القوى المؤثرة لا تشمل قوى رد الفعل عند نقاط اتصال الأجسام. أخيرًا، افترض أن السرعة Vᵢ والسرعات الزاوية ωᵢ ، حيث i = 1، ...، n ، لكل جسم صلب، تُحدد بإحداثي عام واحد q. يُقال إن نظام الأجسام الصلبة هذا له درجة حرية واحدة .

يُعطى الشغل الافتراضي للقوى والعزوم، F i و T i ، المطبقة على هذا النظام ذي درجة الحرية الواحدة، بالعلاقة التالية:دلتادبليو=أنا=1ن(FأناVأناq˙+تيأناωأناq˙)دلتاq=سؤالدلتاq،{\displaystyle \delta W=\sum _{i=1}^{n}\left(\mathbf {F} _{i}\cdot {\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}}}+\mathbf {T} _{i}\cdot {\frac {\partial {\boldsymbol {\omega }}_{i}}{\partial {\dot {q}}}}\right)\delta q=Q\delta q,} أين سؤال=أنا=1ن(FأناVأناq˙+تيأناωأناq˙)،{\displaystyle Q=\sum _{i=1}^{n}\left(\mathbf {F} _{i}\cdot {\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}}}+\mathbf {T} _{i}\cdot {\frac {\partial {\boldsymbol {\omega }}_{i}}{\partial {\dot {q}}}}\right),} هي القوة المعممة المؤثرة على هذا النظام ذي درجة الحرية الواحدة.

إذا تم تعريف النظام الميكانيكي بواسطة m إحداثيات معممة، q j ، j = 1، ...، m ، فإن النظام لديه m درجات حرية ويتم إعطاء الشغل الافتراضي بواسطة، دلتادبليو=ج=1مسؤالجدلتاqج،{\displaystyle \delta W=\sum _{j=1}^{m}Q_{j}\delta q_{j},} أين سؤالج=أنا=1ن(FأناVأناq˙ج+تيأناωأناq˙ج)،ج=1،...،م.{\displaystyle Q_{j}=\sum _{i=1}^{n}\left(\mathbf {F} _{i}\cdot {\frac {\partial \mathbf {V} _{i}}{\partial {\dot {q}}_{j}}}+\mathbf {T} _{i}\cdot {\frac {\partial {\boldsymbol {\omega }}_{i}}{\partial {\dot {q}}_{j}}}\right),\quad j=1,\ldots ,m.} هي القوة المعممة المرتبطة بالإحداثي المعمم q j . ينص مبدأ الشغل الافتراضي على أن التوازن السكوني يحدث عندما تكون هذه القوى المعممة المؤثرة على النظام صفرًا، أي سؤالج=0،ج=1،...،م.{\displaystyle Q_{j}=0,\quad j=1,\ldots ,m.}

تحدد هذه المعادلات m التوازن الساكن لنظام الأجسام الصلبة.

قوى القصور الذاتي المعممة

لنفترض جسمًا صلبًا واحدًا يتحرك تحت تأثير قوة محصلة F وعزم دوران T ، وله درجة حرية واحدة مُعرَّفة بالإحداثي المعمم q . بافتراض أن نقطة مرجع القوة المحصلة وعزم الدوران هي مركز كتلة الجسم، فإن قوة القصور الذاتي المعممة Q* المرتبطة بالإحداثي المعمم q تُعطى بالعلاقة التالية: سؤال*=-(مأ)Vq˙-([أناR]α+ω×[أناR]ω)ωq˙.{\displaystyle Q^{*}=-(M\mathbf {A} )\cdot {\frac {\partial \mathbf {V} }{\partial {\dot {q}}}}-\left([I_{R}]{\boldsymbol {\alpha }}+{\boldsymbol {\omega }}\times [I_{R}]{\boldsymbol {\omega }}\right)\cdot {\frac {\partial {\boldsymbol {\omega }}}{\partial {\dot {q}}}}.}

يمكن حساب قوة القصور الذاتي هذه من الطاقة الحركية للجسم الصلب، تي=12مVV+12ω[أناR]ω،{\displaystyle T={\tfrac {1}{2}}M\mathbf {V} \cdot \mathbf {V} +{\tfrac {1}{2}}{\boldsymbol {\omega }}\cdot [I_{R}]{\boldsymbol {\omega }},} باستخدام الصيغة سؤال*=-(ددتتيq˙-تيq).{\displaystyle Q^{*}=-\left({\frac {d}{dt}}{\frac {\partial T}{\partial {\dot {q}}}}-{\frac {\partial T}{\partial q}}\right).}

يمتلك نظام مكون من n جسمًا صلبًا بإحداثيات معممة m طاقة حركية تي=أنا=1ن(12مVأناVأنا+12ωأنا[أناR]ωأنا)،{\displaystyle T=\sum _{i=1}^{n}\left({\tfrac {1}{2}}M\mathbf {V} _{i}\cdot \mathbf {V} _{i}+{\tfrac {1}{2}}{\boldsymbol {\omega }}_{i}\cdot [I_{R}]{\boldsymbol {\omega }}_{i}\right),} والتي يمكن استخدامها لحساب قوى القصور الذاتي المعممة m [ 6 ]سؤالج*=-(ددتتيq˙ج-تيqج)،ج=1،...،م.{\displaystyle Q_{j}^{*}=-\left({\frac {d}{dt}}{\frac {\partial T}{\partial {\dot {q}}_{j}}}-{\frac {\partial T}{\partial q_{j}}}\right),\quad j=1,\ldots ,m.}

التوازن الديناميكي

ينص مبدأ دالمبير للعمل الافتراضي على أن نظامًا من الأجسام الصلبة يكون في حالة اتزان ديناميكي عندما يكون العمل الافتراضي الناتج عن مجموع القوى المؤثرة وقوى القصور الذاتي مساويًا للصفر لأي إزاحة افتراضية للنظام. وبالتالي، فإن الاتزان الديناميكي لنظام مكون من n جسمًا صلبًا ذي m إحداثيات معممة يتطلب ما يلي: دلتادبليو=(سؤال1+سؤال1*)دلتاq1++(سؤالم+سؤالم*)دلتاqم=0،{\displaystyle \delta W=\left(Q_{1}+Q_{1}^{*}\right)\delta q_{1}+\dots +\left(Q_{m}+Q_{m}^{*}\right)\delta q_{m}=0,} لأي مجموعة من الإزاحات الافتراضية δq j . ينتج عن هذا الشرط m معادلة، سؤالج+سؤالج*=0،ج=1،...،م،{\displaystyle Q_{j}+Q_{j}^{*}=0,\quad j=1,\ldots ,m,} والتي يمكن كتابتها أيضاً على النحو التالي ددتتيq˙ج-تيqج=سؤالج،ج=1،...،م.{\displaystyle {\frac {d}{dt}}{\frac {\partial T}{\partial {\dot {q}}_{j}}}-{\frac {\partial T}{\partial q_{j}}}=Q_{j},\quad j=1,\ldots ,m.} والنتيجة هي مجموعة من m معادلات الحركة التي تحدد ديناميكيات نظام الجسم الصلب.

معادلات لاغرانج

إذا كانت القوى المعممة Q j قابلة للاشتقاق من طاقة الوضع V ( q 1 , ..., q m ) ، فإن معادلات الحركة هذه تأخذ الشكل التالي ددتتيq˙ج-تيqج=-Vqج،ج=1،...،م.{\displaystyle {\frac {d}{dt}}{\frac {\partial T}{\partial {\dot {q}}_{j}}}-{\frac {\partial T}{\partial q_{j}}}=-{\frac {\partial V}{\partial q_{j}}},\quad j=1,\ldots ,m.}

في هذه الحالة، نُدخل دالة لاغرانج ، L = TV ، فتصبح معادلات الحركة هذه ددتلq˙ج-لqج=0ج=1،...،م.{\displaystyle {\frac {d}{dt}}{\frac {\partial L}{\partial {\dot {q}}_{j}}}-{\frac {\partial L}{\partial q_{j}}}=0\quad j=1,\ldots ,m.} تُعرف هذه باسم معادلات لاغرانج للحركة .

الزخم الخطي والزاوي

نظام الجسيمات

يُصاغ الزخم الخطي والزاوي لنظام جسيمات صلب بقياس موضع وسرعة الجسيمات بالنسبة لمركز الكتلة. لنفترض أن نظام الجسيمات Pᵢ ، حيث i = 1، ...، n ، يقع عند الإحداثيات rᵢ والسرعات vᵢ . اختر نقطة مرجعية R واحسب متجهي الموضع والسرعة النسبيين. رأنا=(رأنا-R)+R،vأنا=ددت(رأنا-R)+V.{\displaystyle \mathbf {r} _{i}=\left(\mathbf {r} _{i}-\mathbf {R} \right)+\mathbf {R} ,\quad \mathbf {v} _{i}={\frac {d}{dt}}(\mathbf {r} _{i}-\mathbf {R} )+\mathbf {V} .}

متجهات الزخم الخطي والزاوي الكلي بالنسبة إلى نقطة المرجع R هي ص=ددت(أنا=1نمأنا(رأنا-R))+(أنا=1نمأنا)V،{\displaystyle \mathbf {p} ={\frac {d}{dt}}\left(\sum _{i=1}^{n}m_{i}\left(\mathbf {r} _{i}-\mathbf {R} \right)\right)+\left(\sum _{i=1}^{n}m_{i}\right)\mathbf {V} ,} و ل=أنا=1نمأنا(رأنا-R)×ددت(رأنا-R)+(أنا=1نمأنا(رأنا-R))×V.{\displaystyle \mathbf {L} =\sum _{i=1}^{n}m_{i}\left(\mathbf {r} _{i}-\mathbf {R} \right)\times {\frac {d}{dt}}\left(\mathbf {r} _{i}-\mathbf {R} \right)+\left(\sum _{i=1}^{n}m_{i}\left(\mathbf {r} _{i}-\mathbf {R} \right)\right)\times \mathbf {V} .}

إذا تم اختيار R كمركز الكتلة، فإن هذه المعادلات تتبسط إلى ص=مV،ل=أنا=1نمأنا(رأنا-R)×ددت(رأنا-R).{\displaystyle \mathbf {p} =M\mathbf {V} ,\quad \mathbf {L} =\sum _{i=1}^{n}m_{i}\left(\mathbf {r} _{i}-\mathbf {R} \right)\times {\frac {d}{dt}}\left(\mathbf {r} _{i}-\mathbf {R} \right).}

نظام صلب من الجسيمات

لتخصيص هذه الصيغ لجسم صلب، افترض أن الجسيمات متصلة ببعضها البعض بشكل صلب بحيث يتم تحديد مواقع P i ، حيث i=1,...,n بواسطة الإحداثيات r i والسرعات v i . حدد نقطة مرجعية R واحسب متجهات الموضع والسرعة النسبية. رأنا=(رأنا-R)+R،vأنا=ω×(رأنا-R)+V،{\displaystyle \mathbf {r} _{i}=(\mathbf {r} _{i}-\mathbf {R} )+\mathbf {R} ,\quad \mathbf {v} _{i}=\omega \times (\mathbf {r} _{i}-\mathbf {R} )+\mathbf {V} ,} حيث ω هي السرعة الزاوية للنظام. [ 7 ] [ 8 ] [ 9 ]

الزخم الخطي والزخم الزاوي لهذا النظام الصلب، مقاسين بالنسبة إلى مركز الكتلة هما: ص=(أنا=1نمأنا)V،ل=أنا=1نمأنا(رأنا-R)×vأنا=أنا=1نمأنا(رأنا-R)×(ω×(رأنا-R)).{\displaystyle \mathbf {p} =\left(\sum _{i=1}^{n}m_{i}\right)\mathbf {V} ,\quad \mathbf {L} =\sum _{i=1}^{n}m_{i}(\mathbf {r} _{i}-\mathbf {R} )\times \mathbf {v} _{i}=\sum _{i=1}^{n}m_{i}(\mathbf {r} _{i}-\mathbf {R} )\times (\omega \times (\mathbf {r} _{i}-\mathbf {R} )).}

تُبسط هذه المعادلات لتصبح، ص=مV،ل=[أناR]ω،{\displaystyle \mathbf {p} =M\mathbf {V} ,\quad \mathbf {L} =[I_{R}]\omega ,} حيث M هي الكتلة الكلية للنظام و[IR ] هي مصفوفة عزم القصور الذاتي المعرفة بواسطة [أناR]=-أنا=1نمأنا[رأنا-R][رأنا-R]،{\displaystyle [I_{R}]=-\sum _{i=1}^{n}m_{i}[r_{i}-R][r_{i}-R],} حيث [r i − R] هي المصفوفة المتناظرة المائلة التي تم إنشاؤها من المتجه r iR.

التطبيقات

  • تحليل الأنظمة الروبوتية
  • التحليل البيوميكانيكي للحيوانات أو البشر أو الأنظمة الشبيهة بالبشر
  • تحليل الأجسام الفضائية
  • فهم الحركات الغريبة للأجسام الصلبة. [ 10 ]
  • تصميم وتطوير أجهزة الاستشعار القائمة على الديناميكيات، مثل أجهزة الاستشعار الجيروسكوبية.
  • تصميم وتطوير تطبيقات متنوعة لتحسين استقرار السيارات.
  • تحسين رسومات ألعاب الفيديو التي تتضمن أجسامًا صلبة

انظر أيضاً

مراجع

  1. ب. بول، علم الحركة وديناميكيات الآلات المستوية، برنتيس هول، نيوجيرسي، 1979
  2. LW Tsai، تحليل الروبوت: ميكانيكا المعالجات التسلسلية والمتوازية، جون وايلي، نيويورك، 1999.
  3. موسوعة بريتانيكا، قوانين نيوتن للحركة .
  4. KJ Waldron و GL Kinzel، علم الحركة والديناميكا، وتصميم الآلات ، الطبعة الثانية، جون وايلي وأولاده، 2004.
  5. توربي، بروس (1984). "أساليب الطاقة". الديناميكا المتقدمة للمهندسين . سلسلة HRW في الهندسة الميكانيكية. الولايات المتحدة الأمريكية: دار نشر كلية سي بي إس. ISBN 0-03-063366-4.
  6. TR Kane and DA Levinson, Dynamics, Theory and Applications , McGraw-Hill, NY, 2005.
  7. ماريون، جيه بي؛ ثورنتون، إس تي (1995). الديناميكا الكلاسيكية للأنظمة والجسيمات (الطبعة الرابعة ). طومسون. ISBN  0-03-097302-3..
  8. سايمون، ك. ر. (1971). الميكانيكا ( الطبعة الثالثة). أديسون-ويسلي. رقم ISBN  0-201-07392-7..
  9. تيننباوم، ر. أ. (2004). أساسيات الديناميكا التطبيقية . سبرينغر. ISBN 0-387-00887-X..
  10. غوميز، ر. و.؛ هيرنانديز-غوميز، ج. ج.؛ ماركينا، ف. (25 يوليو 2012). "أسطوانة قافزة على مستوى مائل" . المجلة الأوروبية للفيزياء 33 ( 5). IOP: 1359–1365 . arXiv : 1204.0600 . Bibcode : 2012EJPh...33.1359G . doi : 10.1088/0143-0807/33/5/1359 . S2CID 55442794. تاريخ الاسترجاع: 25 أبريل 2016 . 

للمزيد من القراءة

  • إي. ليمانيس (1965). المسألة العامة لحركة الأجسام الصلبة المقترنة حول نقطة ثابتة. ( سبرينغر ، نيويورك).
  • WB Heard (2006). ميكانيكا الأجسام الصلبة: الرياضيات والفيزياء والتطبيقات. ( Wiley-VCH ).