مبدأ العمل الأقل في الميكانيكا التحليلية

الخلفية




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

مشكلة في حركة الكرة

ليس بهذه البساطة


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

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

الصورة

كلا النوعين من حركة الكرة ممكنان ، لكن الإجراء في الحالة الثانية سيظهر أكثر ، لأنه طوال هذا الوقت ستتحرك الكرة بطاقة حركية غير صفرية.

كيف تحفظ مبدأ العمل الأقل بحيث يكون عادلاً في مثل هذه المواقف؟ سنتحدث عن هذا في المرة القادمة.

فما هي المشكلة في رأيي؟

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

ستتحدث هذه المقالة عن كيفية بناء الميكانيكا على أساس PND ، وستحاول أن توضح للقارئ أن المشكلة التي طرحها مؤلف المنشور المذكور مفقودة.

1. تعريف عمل هاميلتون. مبدأ العمل الأقل


يسمى عمل هاميلتون وظيفية

S = i n t l i m i t s t 2 t 1 \ ، L l e f t ( m a t h b f q ( t ) ، d o t m a t h b f q ( t ) r i g ح ر ) \ ، د ر       


أين

L يسار( mathbfq(t)، dot mathbfq(t) right)=T left( mathbfq(t)، dot mathbfq(t) right) Pi( mathbfq)


هي وظيفة Lagrange لبعض النظام الميكانيكي حيث (حذف الحجج أدناه) T هي الطاقة الحركية للنظام ؛ P - طاقته الكامنة ؛ ف (ر) هو ناقل الإحداثيات المعممة لهذا النظام ، وهي دالة للوقت. ويعتقد أنه يتم إصلاح الغرز الزمنية t 1 و t 2 .

لماذا الوظيفة لا تعمل؟ لأن الوظيفة ، بحكم تعريفها ، هي قاعدة يتم بموجبها ربط رقم من مجال التعريف (وسيطة دالة) برقم آخر من مجال القيم. الوظيفة تختلف من حيث أن الوسيطة ليست رقمًا ، ولكنها دالة كاملة. في هذه الحالة ، هذا هو قانون حركة النظام الميكانيكي q (t) ، المحدد على الأقل خلال الفترة الزمنية بين t 1 و t 2 .

سمحت لنا أعمال طويلة الأمد (وهذا ما يُطرح بشكل معتدل!) أعمال علماء الميكانيكا (بما في ذلك ليونارد أويلر المثير للإعجاب) بصياغة

مبدأ العمل الأقل:
نظام ميكانيكي تم تحديد وظيفة Lagrange له L يسار( mathbfq(t)، dot mathbfq(t) right) ، يتحرك بطريقة تجعل قانون حركته q (t) يوفر الحد الأدنى للوظيفة

S= int limitst2t1L left( mathbfq(t)، dot mathbfq(t) right)dt to min


دعا عمل هاميلتون.
بالفعل من تعريف PND ذاته ، فإنه يتبع حقيقة أن هذا المبدأ يؤدي إلى معادلات الحركة فقط لفئة محدودة من الأنظمة الميكانيكية. من اجل ماذا؟ ودعنا نكتشف ذلك.

2. حدود تطبيق مبدأ العمل الأقل. بعض التعاريف للأصغر


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

قوة العمل أثناء التنقل
ضع في اعتبارك نقطة تتحرك على طول المسار AB ، الذي يتم تطبيق القوة عليه  vecF . يتم تحديد الإزاحة الصغيرة بشكل لا نهائي لنقطة على طول المسار بواسطة المتجه d vecs المماس الموجه إلى المسار.

العمل الأولي للقوة  vecF على التحرك d vecs تسمى القيمة العددية تساوي

dA= vecF cdotd vecs


ثم ، فإن العمل الكامل للقوة على تحريك النقطة على طول المسار AB هو تكامل منحني

A= int limitsAB vecF cdotd vecs



نقاط الطاقة الحركية
الطاقة الحركية للنقطة T هي العمل الذي يجب أن تكمله القوى المطبقة على نقطة الكتلة m من أجل تحويل النقطة من الحركة إلى الراحة بسرعة  vecv
نحسب الطاقة الحركية وفقًا لهذا التعريف. دع النقطة تبدأ بالانتقال من حالة الراحة تحت تصرف القوى المطبقة عليها. في الجزء من المسار AB ، يكتسب السرعة  vecv . نحسب العمل الذي قامت به القوات المطبقة على النقطة ، والتي ، وفقًا لمبدأ استقلالية القوات ، نستبدل الناتج  vecF

T=A= int limitsAB vecF cdotd vecs


وفقًا لقانون نيوتن الثاني

 vecF=m veca=m fracd vecvdt


ثم

T= int limitsAB vecF cdotd vecs=m int limitsAB fracd vecvdt cdotd vecs=m int limitsAB vecv cdotd vecv


نحن نحسب المنتج الحجمي القائم بصرامة تحت العلامة المتكاملة ، والتي نميز من خلالها المنتج الحجمي لمتجه السرعة في حد ذاته

 fracddt( vecv cdot vecv)= fracd vecvdt cdot vecv+ vecv cdot fracd vecvdt=2 vecv cdot fracd vecvdt quad(1)


من ناحية أخرى

 vecv cdot vecv=v2


التفريق بين هذه المساواة فيما يتعلق بالوقت لدينا

 fracddt( vecv cdot vecv)= fracddt(v2)=2v fracdvdt quad(2)


بمقارنة (1) و (2) ، نستنتج ذلك

 vecv cdotd vecv=vdv


ثم نحسب العمل بهدوء ، ونكشف عن التكامل المنحني الخطي من خلال واحد محدد ، مع أخذ معامل السرعة للنقطة في بداية ونهاية المسار كحدود

T=m int limitsAB vecv cdotd vecv=m int limitsv0vdv= fracmv22



القوى المحافظة ونقاط الطاقة المحتملة
ضع في اعتبارك القوة التي تعمل على نقطة ، بحيث يعتمد حجم واتجاه هذه القوة فقط على موقع النقطة في الفضاء

 vecF= vecF(x،y،z) quad(3)


دع النقطة تتحرك في الفضاء على طول المسار التعسفي AB. نحسب العمل الذي ستقوم به القوة (3)

A= int limitsAB vecF cdotd vecs= int limitsAB left(Fxdx+Fydy+Fzdz right)


نظرًا لأن إسقاط القوة على محور الإحداثيات يعتمد فقط على هذه الإحداثيات نفسها ، يمكن للمرء دائمًا العثور عليها

U=U(x،y،z)


مثل هذا

Fx= frac جزئيU جزئيX، quadFy= frac جزئيU جزئيy، quadFz= frac جزئيU جزئيz


بعد ذلك ، يتم تحويل تعبير العمل إلى

A= int limitsAB left( frac جزئيU جزئيxdx+ frac جزئيU جزئيydy+ frac جزئيU جزئيzdz right)= int limitsUBUAdU=UBUA


أين UA، ،UB هي قيم الدالة U (x، y، z) عند النقطتين A و B ، على التوالي. وبالتالي ، لا يعتمد عمل القوة التي نفكر فيها على مسار النقطة ، ولكن يتم تحديده فقط من خلال قيم الدالة U في بداية ونهاية المسار. تسمى هذه القوة قوة محافظة ، وتسمى الوظيفة المقابلة U (x ، y ، z) دالة القوة. من الواضح ،  vecF= nablaU ، وكذلك المساواة في الصفر لعمل القوة المحافظة عند التحرك على طريق مغلق. يقال أيضًا أن الدالة U (x، y، z) تحدد مجال القوة في الفضاء.
الطاقة الكامنة  Pi= Pi(x،y،z) تسمى النقاط في الفضاء مع مجال قوة معين عمل القوى الخارجية المطبقة عليه ، والتي تؤديها عند تحريك النقطة إلى موقع في الفضاء المحدد بواسطة الإحداثيات (س ، ص ، ض) من بعض الموقع التعسفي المحدد كنقطة مرجعية لمستوى الطاقة المحتملة .
دعونا نختار نقطة اعتباطية O تقع بين النقطتين A و B على مسار النقطة التي تم النظر فيها سابقًا. نفترض أن الطاقة الكامنة تساوي صفر عند النقطة O. ثم حسب التعريف

 PiA=(UAUO)


هي الطاقة الكامنة للنقطة في الموضع A و

 PiB=(UBUO)


- الطاقة الكامنة للنقطة في الموضع ب. وبالنظر إلى كل ما سبق ، فإننا نحسب مرة أخرى عمل القوى المحتملة على الانتقال من النقطة أ إلى النقطة ب

AAB=AAO+AOB=UOUA+UBUO=(UOUA)(UOUB)= PiA PiB


وبالتالي ، فإن عمل القوى المحافظة يساوي التغيير في الطاقة الكامنة لنقطة مأخوذة بعلامة معاكسة

AAB= PiA PiB=( PiB PiA)= Delta Pi


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

3. مفهوم التغيرات في الإحداثيات المعممة. بيان المشكلة التباينية


لذلك ، نعتبر الآن نظامًا ميكانيكيًا يتحرك تحت تأثير القوى المحتملة ، ويتم تحديد موقعه بشكل فريد بواسطة ناقل الإحداثيات المعممة

 mathbfq= يسار[q1،q2، dots،qs right]T quad(4)


حيث s هو عدد درجات الحرية لنظام معين.

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


الشكل 1. الحركة الفعلية والدائرية لنظام ميكانيكي

في الشكل ، التبعية qi(t) يصور منحنى أحمر. نختار اثنتين من الدعامات الثابتة العشوائية t 1 و t 2 ، ونضع t 2 > t 1 . موضع النظام  mathbfq1= mathbfq(t1) نتفق على استدعاء الموقف الأولي للنظام ، و  mathbfq2= mathbfq(t2) - الموضع النهائي للنظام.

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

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

 deltaqi(t)=qi(t)qi(t)، quad foralli= overline1،s quad(5)


ستسمى الاختلافات المتزامنة في الإحداثيات المعممة

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

في الواقع ، نحن نغير قانون الحركة ، والذي بموجبه ينتقل النظام من الحالة الأولية إلى الحالة النهائية. يتم تحديد الحالة الأولية والنهائية من خلال قانون الحركة الفعلي ، لكنني أؤكد مرة أخرى أننا لا نعرف قيمها المحددة ويمكن أن تكون ممكنة بشكل حركي ، نحن نعتقد فقط أنها موجودة وأن النظام مضمون للانتقال من موقع إلى آخر! في الموضع الأولي والنهائي للنظام ، لا نغير قانون الحركة ، وبالتالي ، فإن الاختلافات في الإحداثيات المعممة في الموضع الأولي والنهائي تساوي الصفر

 deltaqi(t1)= deltaqi(t2)=0، quad foralli= overline1،s quad(6)


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

 deltaS= delta int limitst2t1L(q1، dots،qs، dotq1، dots، dotqs)dt=0 quad(7)



4. حل المشكلة التباينية. معادلات لاغرانج من النوع الثاني


دعونا نحل مشكلتنا التباينية ، التي نحسب لها التغيير الكامل للعمل الوظيفي ونعادلها بالصفر

\ start {align} \ delta S = & \ int \ limits_ {t_1} ^ {t2} \، L (q_1 + \ delta q_1، \ dots، q_s + \ delta q_s، \، \ dot q_1 + \ delta \ dot q_1 ، \ dots ، \ dot q_s + \ delta \ dot q_s) \ ، dt - \\ & - \ int \ limits_ {t_1} ^ {t2} \ ، L (q_1 ، \ dots ، q_s ، \ ، \ dot q_1 ، \ dots ، \ dot q_s) \ ، dt = 0 \ end {align}


دعنا ندفع كل شيء تحت تكامل واحد ، وبما أن جميع العمليات على الكميات المتناهية الصغر صالحة للتغيير ، فإننا نحول هذا التمساح إلى الشكل

 int limitst2t1 left[ sum limitssi=1 frac جزئيL جزئيqi deltaqi+ sum limitssi=1 frac جزئيL جزئي dotqi delta dotqi right]dt=0 quad(8)


بناء على تعريف السرعة المعممة

 delta dotqi= fracd( deltaqi)dt


ثم يتم تحويل التعبير (8) إلى الشكل

 int limitst2t1 left[ sum limitssi=1 frac جزئيL جزئيqi deltaqidt+ sum limitssi=1 frac جزئيL جزئي dotqiد( deltaqi) right]=0


يتم دمج المصطلح الثاني في أجزاء

 sum limitssi=1 frac جزئيL جزئي dotqi deltaqi|t2t1+ int limitst2t1 left[ sum limitssi=1 frac جزئيL جزئيqi deltaqi sum limitssi=1 fracddt اليسار( frac جزئيL جزئي نقطةqi يمين) deltaqi right]dt=0 quad(10)


بناء على الشرط (7) لدينا

 sum limitssi=1 frac جزئيL جزئي dotqi deltaqi|t2t1=0


ثم نحصل على المعادلة

 int limitst2t1 left[ sum limitssi=1 left( frac جزئيL جزئيqi fracddt left( frac جزئيL جزئي نقطةqi يمين) يمين) deltaqi right]dt=0


بالنسبة لحدود التكامل التعسفي ، يتم ضمان اختفاء جزء لا يتجزأ من خلال اختفاء التكامل

 sum limitssi=1 left[ frac جزئيL جزئيqi fracddt left( frac جزئيL جزئي dotqi right) right] deltaqi=0 quad(11)


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

 frac جزئيL جزئيqi fracddt left( frac جزئيL جزئي dotqi يمين)=0، quad foralli= overline1،s


لا أحد يزعجنا في ضرب كل من المعادلات في (-1) والحصول على تدوين أكثر دراية

 fracddt left( frac جزئيL جزئي نقطةqi يمين) frac جزئيL جزئيqi=0، quad foralli= overline1،s quad(12)



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

هذا كل شيء ، ينتهي مبدأ العمل الأقل هناك ، وتبدأ نظرية المعادلات التفاضلية العادية ، والتي ، على وجه الخصوص ، تنص على أن حل المعادلة (12) هو دالة متجهة للشكل

 mathbfq= mathbfq(t،C1،C2، dots،C2s)


حيث C 1 ، ... ، C 2s هي ثوابت تكاملية تعسفية.

بهذه الطريقة
PND هو مبدأ أساسي يسمح للمرء بالحصول على معادلات الحركة لنظام يتم تعريف وظيفة Lagrange من أجله
نقطة! في مشاكل الميكانيكا التحليلية ، لم تعد الحسابات المذكورة أعلاه بحاجة إلى القيام بها ، يكفي استخدام نتائجها (12). الوظيفة التي تحقق المعادلة (12) هي قانون الحركة لنظام يرضي PND.

5. المشكلة مع الكرة والجدار


الآن نعود إلى المهمة التي بدأ بها كل شيء - حول الحركة أحادية البعد للكرة بالقرب من جدار مرن تمامًا. بالطبع ، بالنسبة لهذه المشكلة ، يمكن للمرء الحصول على معادلات تفاضلية للحركة. نظرًا لأن هذه معادلات تفاضلية للحركة ، فإنني أؤكد ذلك ، أي من حلولها توفر الحد الأدنى من الإجراءات الوظيفية ، مما يعني تنفيذ PND! يمكن تمثيل الحل العام لمعادلات الحركة للكرة في شكل ما يسمى صورة طورية للنظام الميكانيكي قيد النظر. صورة طور هذه


الشكل 2. صورة طور النظام في مشكلة الكرة

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

هذا هو حل عام لمعادلة حركة الكرة. يوفر كل مسار من هذه المراحل الحد الأدنى من الوظيفة الوظيفية ، والتي تتبع مباشرة من الحسابات التي تم إجراؤها أعلاه.

ماذا يفعل مؤلف المهمة؟ يقول: هنا الكرة في حالة استراحة ، وللفترات الزمنية من t A إلى t B تكون الحركة صفر. إذا تم دفع الكرة ضد الجدار ، فستكون الحركة في نفس الفترة الزمنية أكبر ، حيث أن الكرة لديها طاقة حركية غير صفرية وغير متغيرة. ولكن لماذا تتحرك الكرة نحو الجدار ، لأن الراحة ستكون أقل؟ لذا فإن PND يواجه مشاكل ولا يعمل! لكننا بالتأكيد سنحل هذا في المقالة التالية.

ما يقوله الكاتب هراء. لماذا؟ نعم ، لأنه يقارن الإجراءات على الفروع المختلفة لنفس مسار المرحلة الحقيقية! في هذه الأثناء ، عند تطبيق PND ، تتم مقارنة الإجراء على المسار الفعلي وعلى العديد من مسارات التقاطع الدائري. , , !

? . . , . . . . . . ? ,

? , , t A t B , , . «» . .

. , , , t B , . , , t B , . , . ? , , . ? ربما. ? , , , .

هذا كل شيء. . .


— . , , , (12). , , (12) . .

.

Source: https://habr.com/ru/post/ar425771/


All Articles