حركة المقذوفات

مسارات مكافئة لنفثات المياه
مكونات السرعة الابتدائية للرمي المكافئ
تكون المسارات الباليستية مكافئة إذا كانت الجاذبية متجانسة، وإهليلجية إذا كانت شعاعية.

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

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

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

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

علم المقذوفات ( من اليونانية القديمة βάλλειν bállein ' يرمي ' ) هو علم الديناميكا الذي يتعامل مع طيران وسلوك وتأثيرات المقذوفات، وخاصة الرصاص والقنابل غير الموجهة والصواريخ وما شابه ذلك؛ علم أو فن تصميم وتسريع المقذوفات لتحقيق الأداء المطلوب. 

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

تتجاهل المعادلات الأساسية لعلم المقذوفات جميع العوامل تقريبًا باستثناء السرعة الابتدائية وزاوية الإطلاق وتسارع الجاذبية الأرضية المفترض ثباته. غالبًا ما تتطلب الحلول العملية لمسألة المقذوفات مراعاة مقاومة الهواء والرياح الجانبية وحركة الهدف وتسارع الجاذبية الأرضية المتغير مع الارتفاع، وفي مسائل مثل إطلاق صاروخ من نقطة على سطح الأرض إلى أخرى، تُؤخذ في الاعتبار المسافة بين الأفق وانحناء الأرض (سرعة دورانها المحلية).v(لأت)=ωR(لأت){\textstyle v(lat)=\أوميغا R(lat)}). لا تحتوي الحلول الرياضية التفصيلية للمشاكل العملية عادةً على حلول مغلقة الشكل ، وبالتالي تتطلب طرقًا عددية لمعالجتها.

المسار في الفراغ

في حركة المقذوفات، تكون الحركة الأفقية والحركة الرأسية مستقلتين عن بعضهما البعض؛ أي أن أياً منهما لا يؤثر على الأخرى. هذا هو مبدأ الحركة المركبة الذي وضعه غاليليو عام 1638، [ 1 ] واستخدمه لإثبات الشكل المكافئ لحركة المقذوفات. [ 2 ]

تكون المكونات الأفقية والرأسية لسرعة المقذوف مستقلة عن بعضها البعض.

المسار الباليستي هو مسار مكافئ ذو تسارع متجانس، كما هو الحال في المركبات الفضائية ذات التسارع الثابت في غياب القوى الأخرى. على الأرض، يتغير مقدار التسارع مع الارتفاع.ز(y)=ز0/(1+y/R)2{\textstyle g(y)=g_{0}/(1+y/R)^{2}}والاتجاه (الأهداف البعيدة) مع خطوط الطول والعرض على طول المسار. ينتج عن ذلك مسار إهليلجي ، وهو قريب جدًا من القطع المكافئ على نطاق صغير. مع ذلك، إذا تم قذف جسم ما واستُبدلت الأرض فجأة بثقب أسود له نفس الكتلة، فسيتضح أن المسار الباليستي جزء من مدار إهليلجي حول ذلك "الثقب الأسود"، وليس قطعًا مكافئًا يمتد إلى ما لا نهاية. عند السرعات العالية، يمكن أن يكون المسار دائريًا أيضًا (كما هو الحال في رواد الفضاء في المدار الأرضي المنخفض ، والأقمار الصناعية الثابتة بالنسبة للأرض على ارتفاع 5 أمتار ).56{\textstyle {\frac {5}{6}}}R)، مكافئ أو زائد (إلا إذا تشوه بفعل أجسام أخرى مثل القمر أو الشمس). [ 3 ]

في هذه المقالة، يتم تناول التسارع الجاذبي المتجانس(ز=ز0){\textstyle (g=g_{0})}يفترض ذلك.

تسريع

بما أن التسارع موجود فقط في الاتجاه الرأسي، فإن السرعة في الاتجاه الأفقي ثابتة، وتساويv0كوسθ{\displaystyle \mathbf {v} _{0}\cos \theta }الحركة الرأسية للمقذوف هي حركة جسيم أثناء سقوطه الحر. هنا، يكون التسارع ثابتًا، ويساوي g . [ ملاحظة 1 ] مركبات التسارع هي:

أx=0{\displaystyle a_{x}=0}،
أy=-ز{\displaystyle a_{y}=-g}.*

*يمكن الإشارة إلى تسارع المحور y أيضًا باسم قوة الأرض(-Fز/م){\textstyle (-F_{g}/m)}على الشيء (الأشياء) محل الاهتمام.

سرعة

لنفترض أن المقذوف قد أُطلق بسرعة ابتدائيةv(0)v0{\displaystyle \mathbf {v} (0)\equiv \mathbf {v} _{0}}والتي يمكن التعبير عنها كمجموع المكونات الأفقية والرأسية على النحو التالي:

v0=v0xx^+v0yy^{\displaystyle \mathbf {v} _{0}=v_{0x}\mathbf {\hat {x}} +v_{0y}\mathbf {\hat {y}} }.

المكوناتv0x{\displaystyle v_{0x}}وv0y{\displaystyle v_{0y}}يمكن إيجادها إذا كانت زاوية الإطلاق الأولية θ معروفة:

v0x=v0كوس(θ){\displaystyle v_{0x}=v_{0}\cos(\theta )}،
v0y=v0الخطيئة(θ){\displaystyle v_{0y}=v_{0}\sin(\theta )}

يظل المركب الأفقي لسرعة الجسم ثابتًا طوال حركته. أما المركب الرأسي للسرعة فيتغير خطيًا، [ ملاحظة 2 ] لأن تسارع الجاذبية الأرضية ثابت. ويمكن حساب مركبات السرعة عند أي لحظة زمنية t بتكامل التسارعين في الاتجاهين x و y ، كما يلي:

vx=v0كوس(θ){\displaystyle v_{x}=v_{0}\cos(\theta )}،
vy=v0الخطيئة(θ)-زت{\displaystyle v_{y}=v_{0}\sin(\theta )-gt}.

مقدار السرعة (وفقًا لنظرية فيثاغورس ، والمعروفة أيضًا بقانون المثلث):

v=vx2+vy2{\displaystyle v={\sqrt {v_{x}^{2}+v_{y}^{2}}}}.

النزوح

الإزاحة وإحداثيات القذف المكافئ

في أي وقتت{\displaystyle t}، الإزاحة الأفقية والرأسية للمقذوف هي:

x=v0تكوس(θ){\displaystyle x=v_{0}t\cos(\theta )}،
y=v0تالخطيئة(θ)-12زت2{\displaystyle y=v_{0}t\sin(\theta )-{\frac {1}{2}}gt^{2}}.

مقدار الإزاحة هو:

Δر=x2+y2{\displaystyle \Delta r={\sqrt {x^{2}+y^{2}}}}.

ضع في اعتبارك المعادلات،

x=v0تكوس(θ){\displaystyle x=v_{0}t\cos(\theta )}وy=v0تالخطيئة(θ)-12زت2{\displaystyle y=v_{0}t\sin(\theta )-{\frac {1}{2}}gt^{2}}[ 4 ]

إذا تم حذف t بين هاتين المعادلتين، فسيتم الحصول على المعادلة التالية:

y=لون برونزي(θ)x-ز2v02كوس2θx2=لون برونزيθx(1-xR).{\displaystyle y=\tan(\theta )\cdot x-{\frac {g}{2v_{0}^{2}\cos ^{2}\theta }}\cdot x^{2}=\tan \theta \cdot x\left(1-{\frac {x}{R}}\right).}

هنا R هو مدى المقذوف .

بما أن g و θ و v₀ ثوابت، فإن المعادلة أعلاه تأخذ الشكل التالي :

y=أx+بx2{\displaystyle y=ax+bx^{2}}،

حيث a و b ثابتان. هذه معادلة قطع مكافئ، لذا فإن المسار قطع مكافئ. محور القطع المكافئ رأسي.

إذا كان موضع المقذوف (x,y) وزاوية الإطلاق ( θ أو α ) معروفين، فيمكن إيجاد السرعة الابتدائية عن طريق حل المعادلة القطعية المذكورة أعلاه لإيجاد v 0 :

v0=x2زxالخطيئة2θ-2yكوس2θ{\displaystyle v_{0}={\sqrt {{x^{2}g} \over {x\sin 2\theta -2y\cos ^{2}\theta }}}}.

الإزاحة في الإحداثيات القطبية

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

ر(ϕ)=2v02كوس2θ|ز|(لون برونزيθثانيةϕ-لون برونزيϕثانيةϕ){\displaystyle r(\phi )={\frac {2v_{0}^{2}\cos ^{2}\theta }{|g|}}\left(\tan \theta \sec \phi -\tan \phi \sec \phi \right)}.

في هذه المعادلة، تمثل نقطة الأصل منتصف المدى الأفقي للمقذوف، وإذا كانت الأرض مستوية، يتم رسم القوس المكافئ في هذا المدى.0ϕπ{\displaystyle 0\leq \phi \leq \pi }يمكن الحصول على هذا التعبير بتحويل المعادلة الديكارتية كما هو مذكور أعلاه بواسطةy=رالخطيئةϕ{\displaystyle y=r\sin \phi }وx=ركوسϕ{\displaystyle x=r\cos \phi }.

مدة الرحلة أو إجمالي مدة الرحلة بأكملها

يُطلق على إجمالي الزمن t الذي يبقى فيه المقذوف في الهواء اسم زمن التحليق. أما المركبة y للمسار فهي

y=y0+v0yت-12زت2{\displaystyle y=y_{0}+v_{0y}t-{\frac {1}{2}}gt^{2}}

بعد انتهاء الرحلة، يعود المقذوف إلى المحور الأفقي (المحور السيني)، لذلكy=0{\displaystyle y=0}.

y0+v0yت-12زت2=0{\displaystyle y_{0}+v_{0y}t-{\frac {1}{2}}gt^{2}=0}

حل هذه المسألة باستخدام القانون العام للمعادلة التربيعية ينتج عنه

ت=v0y±v0y2+2زy0ز{\displaystyle t={\frac {v_{0y}\pm {\sqrt {v_{0y}^{2}+2gy_{0}}}}{g}}}

لاحظ أننا نريد الجذر التربيعي الموجب في هذه الحالة. بافتراض أن المقذوف يبدأ من 0 (y0=0{\displaystyle y_{0}=0}) وباستخدامv0y=v0الخطيئة(θ){\displaystyle v_{0y}=v_{0}\sin(\theta )}وهذا يعطي:

ت=2v0الخطيئة(θ)ز{\displaystyle t={\frac {2v_{0}\sin(\theta )}{g}}}

لاحظ أننا أهملنا مقاومة الهواء على المقذوف.

إذا كانت نقطة البداية على ارتفاع y 0 بالنسبة لنقطة الاصطدام، فإن زمن الطيران هو:

ت=دvكوسθ=vالخطيئةθ+(vالخطيئةθ)2+2زy0ز{\displaystyle t={\frac {d}{v\cos \theta }}={\frac {v\sin \theta +{\sqrt {(v\sin \theta )^{2}+2gy_{0}}}}{g}}}

كما سبق، يمكن اختزال هذا التعبير ( y 0 يساوي 0) إلى

ت=vالخطيئةθ+(vالخطيئةθ)2ز=2vالخطيئةθز،{\displaystyle t={\frac {v\sin {\theta }+{\sqrt {(v\sin {\theta })^{2}}}}{g}}={\frac {2v\sin {\theta }}{g}},}=2vالخطيئة(45)|ز|=2v|ز|{\displaystyle {\frac {2v\sin {(45^{\circ })}}{|g|}}={\frac {{\sqrt {2}}v}{|g|}}}

إذا كانت θ تساوي 45 درجة.

لوvy0=0{\displaystyle v_{y0}=0}مثال على ذلك بالنسبة لجسم يسقط من حالة السكون (الرأسي) على ارتفاعy0{\displaystyle y_{0}}وعند السقوط الحر تحت تأثير الجاذبية، يصبح زمن التحليق،

ت=2y0ز{\displaystyle t={\sqrt {\frac {2y_{0}}{g}}}}

هذه النتيجة مستقلة عن أي سرعة أفقية ابتدائية، كما هو موضح في الرسوم المتحركة أعلاه.

زمن الرحلة إلى موقع الهدف

كما هو موضح أعلاه في قسم الإزاحة ، فإن السرعة الأفقية والرأسية للمقذوف مستقلة عن بعضها البعض.

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

x=v0تكوس(θ){\displaystyle x=v_{0}t\cos(\theta )}

xت=v0كوس(θ){\displaystyle {\frac {x}{t}}=v_{0}\cos(\theta )}

ت=xv0كوس(θ){\displaystyle t={\frac {x}{v_{0}\cos(\theta )}}}

ستعطي هذه المعادلة إجمالي الوقت t الذي يجب أن يقطعه المقذوف للوصول إلى الإزاحة الأفقية للهدف، مع إهمال مقاومة الهواء.

أقصى ارتفاع للمقذوف

أقصى ارتفاع للمقذوف

يُعرف أقصى ارتفاع سيبلغه الجسم بذروة حركته. وسيستمر ازدياد الارتفاع حتىvy=0{\displaystyle v_{y}=0}، إنه،

0=v0الخطيئة(θ)-زتح{\displaystyle 0=v_{0}\sin(\theta )-gt_{h}}.

الوقت اللازم للوصول إلى أقصى ارتفاع (h):

تح=v0الخطيئة(θ)|ز|{\textstyle t_{h}={\frac {v_{0}\sin(\theta )}{|g|}}}.

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

لحساب الإزاحة الرأسية لأقصى ارتفاع للمقذوف:

ح=v0تحالخطيئة(θ)-12زتح2{\displaystyle h=v_{0}t_{h}\sin(\theta )-{\frac {1}{2}}gt_{h}^{2}}
ح=v02الخطيئة2(θ)2|ز|{\displaystyle h={\frac {v_{0}^{2}\sin ^{2}(\theta )}{2|g|}}}

يتم الحصول على أقصى ارتفاع يمكن الوصول إليه عند θ = 90°:

حمأx=v022|ز|{\displaystyle h_{\mathrm {max} }={\frac {v_{0}^{2}}{2|g|}}}

إذا كان موضع المقذوف (x,y) وزاوية الإطلاق (θ) معروفين، فيمكن إيجاد أقصى ارتفاع عن طريق حل المعادلة التالية لإيجاد قيمة h:

ح=(xلون برونزيθ)24(xلون برونزيθ-y).{\displaystyle h={\frac {(x\tan \theta )^{2}}{4(x\tan \theta -y)}}.}

زاوية الارتفاع (φ) عند أقصى ارتفاع تُعطى بالمعادلة التالية:

ϕ=دالة الظل العكسيلون برونزيθ2{\displaystyle \phi =\arctan {\tan \theta \over 2}}

العلاقة بين المدى الأفقي والارتفاع الأقصى

العلاقة بين المدى d على المستوى الأفقي والارتفاع الأقصى h الذي تم بلوغه عندتد2{\displaystyle {\frac {t_{d}}{2}}}يكون:

ح=دلون برونزيθ4{\displaystyle h={\frac {d\tan \theta }{4}}}

أقصى مسافة للمقذوف

أقصى مسافة للمقذوف

لا يعتمد مدى المقذوف وارتفاعه الأقصى على كتلته. لذا، يتساوى المدى والارتفاع الأقصى لجميع الأجسام التي تُقذف بنفس السرعة والاتجاه. المدى الأفقي (d) للمقذوف هو المسافة الأفقية التي يقطعها عند عودته إلى ارتفاعه الابتدائي.y=0{\textstyle y=0}).

0=v0تدالخطيئة(θ)-12زتد2{\displaystyle 0=v_{0}t_{d}\sin(\theta )-{\frac {1}{2}}gt_{d}^{2}}.

حان وقت الوصول إلى الأرض:

تد=2v0الخطيئة(θ)|ز|{\displaystyle t_{d}={\frac {2v_{0}\sin(\theta )}{|g|}}}.

من الإزاحة الأفقية، أقصى مسافة للمقذوف: [ 5 ]

د=v0تدكوس(θ){\displaystyle d=v_{0}t_{d}\cos(\theta )}،

لذا [ ملاحظة 3 ]

د=v02|ز|الخطيئة(2θ).{\displaystyle d={\frac {v_{0}^{2}}{|g|}}\sin(2\theta ).}

لاحظ أن قيمة d تصل إلى أقصى قيمة لها عندما

الخطيئة(2θ)=1،{\displaystyle \sin(2\theta )=1,}

وهو ما يتوافق بالضرورة مع2θ=90{\textstyle 2\theta =90^{\circ }}، أوθ=45{\textstyle \theta =45^{\circ }}.

مسارات مقذوفات أُطلقت بزوايا ارتفاع مختلفة ولكن بنفس السرعة 10 م/ث، في فراغ ومجال جاذبية منتظم لأسفل مقداره 10 م/ث² . تقع النقاط على فترات زمنية قدرها 0.05 ثانية، ويتناسب طول ذيل المقذوف طرديًا مع سرعته. t = الزمن منذ الإطلاق، T = زمن التحليق، R = المدى، H = أعلى نقطة في المسار (مُشار إليها بالأسهم).

المسافة الأفقية الكلية (d) المقطوعة. [ 6 ]

د=vكوسθ|ز|(vالخطيئةθ+(vالخطيئةθ)2+2زy0){\displaystyle d={\frac {v\cos \theta }{|g|}}\left(v\sin \theta +{\sqrt {(v\sin \theta )^{2}+2gy_{0}}}\right)}

عندما يكون السطح مستوياً (الارتفاع الأولي للجسم يساوي صفرًا)، تكون المسافة المقطوعة: [ 7 ]

د=v2الخطيئة(2θ)|ز|{\displaystyle d={\frac {v^{2}\sin(2\theta )}{|g|}}}

وبالتالي، تُحقق أقصى مسافة عندما تكون الزاوية θ تساوي 45 درجة. وهذه المسافة هي:

دمأx=v2|ز|{\displaystyle d_{\mathrm {max} }={\frac {v^{2}}{|g|}}}

تطبيق نظرية الشغل والطاقة

وفقًا لنظرية الشغل والطاقة، فإن المركبة الرأسية للسرعة هي:

vy2=(v0الخطيئةθ)2-2زy{\displaystyle v_{y}^{2}=(v_{0}\sin \theta )^{2}-2gy}.

تتجاهل هذه الصيغ مقاومة الهواء الديناميكية وتفترض أيضًا أن منطقة الهبوط تقع على ارتفاع موحد 0.

زاوية الوصول

"زاوية الوصول" هي الزاوية ( θ ) التي يجب إطلاق المقذوف عندها من أجل قطع مسافة d ، بالنظر إلى السرعة الابتدائية v .

الخطيئة(2θ)=زدv2{\displaystyle \sin(2\theta )={\frac {gd}{v^{2}}}}

هناك حلان:

θ=12دالة الجيب العكسية(زدv2){\displaystyle \theta ={\frac {1}{2}}\arcsin \left({\frac {gd}{v^{2}}}\right)}(مسار ضحل)

ولأنالخطيئة(2θ)=كوس(2θ-90){\displaystyle \sin(2\theta )=\cos(2\theta -90^{\circ })}،

θ=45+12أركوس(زدv2){\displaystyle \theta =45^{\circ }+{\frac {1}{2}}\arccos \left({\frac {gd}{v^{2}}}\right)}(مسار حاد)

الزاوية θ المطلوبة للوصول إلى الإحداثيات ( x , y )

مسار مقذوف في الفراغ لزوايا إطلاق مختلفة. سرعة الإطلاق ثابتة لجميع الزوايا، 50  م/ث، وتسارع الجاذبية الأرضية 10  م/ ث² .

لإصابة هدف على مسافة x وارتفاع y عند إطلاقه من (0,0) وبسرعة ابتدائية فإن زاوية (زوايا) الإطلاق المطلوبة θ هي:

θ=دالة الظل العكسي(v2±v4-ز(زx2+2yv2)زx){\displaystyle \theta =\arctan {\left({\frac {v^{2}\pm {\sqrt {v^{4}-g(gx^{2}+2yv^{2})}}}{gx}}\right)}}

يمثل الجذران للمعادلة زاويتي الإطلاق المحتملتين، ما لم تكونا تخيليتين، وفي هذه الحالة لن تكون السرعة الابتدائية كافية للوصول إلى النقطة ( س ، ص ) المختارة. تسمح هذه الصيغة بإيجاد زاوية الإطلاق المطلوبة دون التقيد بـy=0{\textstyle y=0}.

يمكن أيضًا التساؤل عن زاوية الإطلاق التي تسمح بأقل سرعة إطلاق ممكنة. يحدث هذا عندما تتساوى الحلول المذكورة أعلاه، مما يعني أن المقدار الموجود أسفل علامة الجذر التربيعي يساوي صفرًا. يتطلب هذا، tan θ = v² / gx، حل معادلة تربيعية لـv2{\displaystyle v^{2}}[ 8 ] ونجد

v2/ز=(2زy±ز4y2+4(1x)2)/2ز=y±y2+x2.{\displaystyle v^{2}/g=(2gy\pm g{\sqrt {4y^{2}+4(1x)^{2}}})/2g=y\pm {\sqrt {y^{2}+x^{2}}}.}

هذا يعطي

θ=دالة الظل العكسي(y/x+y2/x2+1).{\displaystyle \theta =\arctan \left(y/x+{\sqrt {y^{2}/x^{2}+1}}\right).}

إذا رمزنا للزاوية التي ظلها y/x بالرمز α ، [ 9 ] فإن

لون برونزيθ=الخطيئةα+1كوسα،{\displaystyle \tan \theta ={\frac {\sin \alpha +1}{\cos \alpha }},}مقلوبه:
لون برونزي(π/2-θ)=كوسαالخطيئةα+1،{\displaystyle \tan(\pi /2-\theta )={\frac {\cos \alpha }{\sin \alpha +1}},}[ 10 ]
كوس2(π/2-θ)=12(الخطيئةα+1){\displaystyle \cos ^{2}(\pi /2-\theta )={\frac {1}{2}}(\sin \alpha +1)}
2كوس2(π/2-θ)-1=كوس(π/2-α){\displaystyle 2\cos ^{2}(\pi /2-\theta )-1=\cos(\pi /2-\alpha )}

وهذا يعني

θ=π/2-12(π/2-α).{\displaystyle \theta =\pi /2-{\frac {1}{2}}(\pi /2-\alpha ).}

بمعنى آخر، يجب أن يكون الإطلاق بزاوية(π/2+α)/2{\textstyle (\pi /2+\alpha )/2}في منتصف المسافة بين الهدف وسمت الرأس (المتجه المعاكس للجاذبية).

إجمالي طول مسار الرحلة

يُعطى طول القوس المكافئ الذي يرسمه المقذوف، L ، بافتراض أن ارتفاع الإطلاق والهبوط متساويان (لا توجد مقاومة هواء)، بالصيغة التالية:

ل=v022ز(2الخطيئةθ+كوس2θln1+الخطيئةθ1-الخطيئةθ)=v02ز(الخطيئةθ+كوس2θtanh-1(الخطيئةθ)){\displaystyle L={\frac {v_{0}^{2}}{2g}}\left(2\sin \theta +\cos ^{2}\theta \cdot \ln {\frac {1+\sin \theta }{1-\sin \theta }}\right)={\frac {v_{0}^{2}}{g}}\left(\sin \theta +\cos ^{2}\theta \cdot \tanh ^{-1}(\sin \theta )\right)}

أينv0{\displaystyle v_{0}}هي السرعة الابتدائية،θ{\displaystyle \theta }زاوية الإطلاق وز{\displaystyle g}يمثل التسارع الناتج عن الجاذبية قيمة موجبة. ويمكن الحصول على هذا التعبير بتقييم تكامل طول القوس للقطع المكافئ للارتفاع والمسافة بين حدود الإزاحة الابتدائية والنهائية (أي بين الصفر والمدى الأفقي للمقذوف) بحيث:

ل=0رأنزهـ1+(دyدx)2دx=0v02الخطيئة(2θ)/ز1+(لون برونزيθ-زv02كوس2θx)2دx.{\displaystyle L=\int _{0}^{\mathrm {range} }{\sqrt {1+\left({\frac {\mathrm {d} y}{\mathrm {d} x}}\right)^{2}}}\,\mathrm {d} x=\int _{0}^{v_{0}^{2}\sin(2\theta )/g}{\sqrt {1+\left(\tan \theta -{g \over {v_{0}^{2}\cos ^{2}\theta }}x\right)^{2}}}\,\mathrm {d} x.}

إذا كان زمن الرحلة هو t ،

ل=0تvx2+vy2دت=02v0الخطيئةθ/ز(زت)2-2زv0الخطيئةθت+v02دت.{\displaystyle L=\int _{0}^{t}{\sqrt {v_{x}^{2}+v_{y}^{2}}}\,\mathrm {d} t=\int _{0}^{2v_{0}\sin \theta /g}{\sqrt {(gt)^{2}-2gv_{0}\sin \theta t+v_{0}^{2}}}\,\mathrm {d} t.}

المسار في الهواء

مسارات كتلة مقذوفة بزاوية 70 درجة:  بدون مقاومة (قطع مكافئ) مع مقاومة ستوكس مع مقاومة نيوتن  

تُنشئ مقاومة الهواء قوة (بالنسبة للمقذوفات المتناظرة) تكون دائماً موجهة ضد اتجاه الحركة في الوسط المحيط ولها مقدار يعتمد على السرعة المطلقة:Fأأنار=-و(v)v^{\displaystyle \mathbf {F_{air}} =-f(v)\cdot \mathbf {\hat {v}} }تعتمد قوة الاحتكاك خطيًا على السرعة (و(v)v{\displaystyle f(v)\propto v}) عند السرعات المنخفضة جدًا ( مقاومة ستوكس ) والتربيعية (و(v)v2{\displaystyle f(v)\propto v^{2}}عند السرعات العالية ( مقاومة نيوتن ). [ 11 ] يتم تحديد الانتقال بين هذه السلوكيات بواسطة رقم رينولدز ، الذي يعتمد على سرعة الجسم وحجمه وكثافته.ρ{\textstyle \rho }واللزوجة الديناميكيةη{\textstyle \eta }بالنسبة لأرقام رينولدز الأقل من 1 تقريبًا، تكون العلاقة خطية، أما فوق 1000 (في حالة التدفق المضطرب ) فتصبح تربيعية. في الهواء، الذي يتمتع بلزوجة حركيةη/ρ{\displaystyle \eta /\rho }حوالي 0.15  سم 2 /ثانية، وهذا يعني أن قوة السحب تصبح تربيعية في v عندما يكون حاصل ضرب سرعة الجسم وقطره أكبر من حوالي 0.015 م 2 /ثانية، وهو ما يحدث عادة بالنسبة للمقذوفات.

  • سحب ستوكس:Fأأنار=-كSتoكهـsv{\displaystyle \mathbf {F_{air}} =-k_{\mathrm {Stokes} }\cdot \mathbf {v} \qquad }Rهـ1{\textstyle Re\lesssim 1})
  • نيوتن دراغ:Fأأنار=-ك|v|v{\displaystyle \mathbf {F_{air}} =-k\,|\mathbf {v} |\cdot \mathbf {v} \qquad }Rهـ1000{\textstyle Re\gtrsim 1000})
مخطط الجسم الحر لجسم لا تؤثر عليه سوى الجاذبية ومقاومة الهواء.

يوضح مخطط الجسم الحر على اليمين مقذوفًا يتعرض لمقاومة الهواء وتأثيرات الجاذبية. ويُفترض هنا أن مقاومة الهواء تكون في الاتجاه المعاكس لسرعة المقذوف.Fأأنار=-و(v)v^{\displaystyle \mathbf {F_{\mathrm {air} }} =-f(v)\cdot \mathbf {\hat {v}} }

مسار مقذوف مع مقاومة ستوكس

سحب ستوكس، حيثFأأنارv{\displaystyle \mathbf {F_{air}} \propto \mathbf {v} }ينطبق هذا فقط عند السرعات المنخفضة جدًا في الهواء، وبالتالي فهو ليس الحالة النموذجية للمقذوفات. ومع ذلك، فإن الاعتماد الخطي لـFأأنار{\displaystyle F_{\mathrm {air} }}علىv{\displaystyle v}ينتج عن ذلك معادلة تفاضلية بسيطة للغاية للحركة

ددت(vxvy)=(-μvx-ز-μvy){\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}{\begin{pmatrix}v_{x}\\v_{y}\end{pmatrix}}={\begin{pmatrix}-\mu \,v_{x}\\-g-\mu \,v_{y}\end{pmatrix}}}

حيث تصبح المركبتان الديكارتيتان مستقلتين تمامًا، وبالتالي يصبح حلها أسهل. [ 12 ] هنا،v0{\displaystyle v_{0}}،vx{\displaystyle v_{x}}وvy{\displaystyle v_{y}}سيتم استخدام الرموز للدلالة على السرعة الابتدائية، والسرعة في اتجاه المحور والسرعة في اتجاه المحور y ، على التوالي. وسيتم استخدام الرمز m للدلالة على كتلة المقذوف.μ:=ك/م{\displaystyle \mu :=k/m} . بالنسبة للاشتقاق، الحالة الوحيدة هي حيث0oθ180o{\textstyle 0^{o}\leq \theta \leq 180^{o}}يتم أخذ ذلك في الاعتبار. مرة أخرى، يتم إطلاق المقذوف من نقطة الأصل (0,0).

x(ت)=vx0μ(1-هـ-μت){\displaystyle x(t)={\frac {v_{x0}}{\mu }}\left(1-e^{-\mu t}\right)}(1ب)
y(ت)=-زμت+1μ(vy0+زμ)(1-هـ-μت){\displaystyle y(t)=-{\frac {g}{\mu }}t+{\frac {1}{\mu }}\left(v_{y0}+{\frac {g}{\mu }}\right)(1-e^{-\mu t})}(3ب)
ت=1μ(1+μزvy0+دبليو(-(1+μزvy0)هـ-(1+μزvy0))){\displaystyle t={\frac {1}{\mu }}\left(1+{\frac {\mu }{g}}v_{y0}+W{\bigl (}-(1+{\frac {\mu }{g}}v_{y0})e^{-(1+{\frac {\mu }{g}}v_{y0})}{\bigr )}\right)}.

مسار مقذوف مع مقاومة نيوتن

مسارات قافز مظلي في الهواء مع مقاومة نيوتن: الطيران المظلي؛ قاذفة قنابل بالقرب من الهدف (فوقه).v.د=103.ν=150{\textstyle v.d=10^{3}.\nu =150}سم² / ثانية؟

إن الحالة الأكثر شيوعًا لمقاومة الهواء ، في حالة أرقام رينولدز التي تزيد عن 1000 تقريبًا، هي مقاومة نيوتن بقوة تتناسب مع مربع السرعة.Fأأنار=-كv2{\displaystyle F_{\mathrm {air} }=-kv^{2}}في الهواء، الذي تبلغ لزوجته الحركية حوالي 0.15  سم 2 /ثانية، هذا يعني أن حاصل ضرب سرعة الجسم وقطره يجب أن يكون أكبر من حوالي 0.015 م 2 /ثانية.

لسوء الحظ، لا يمكن حل معادلات الحركة تحليليًا بسهولة في هذه الحالة. لذلك، سيتم دراسة حل عددي.

تم وضع الافتراضات التالية:

Fد=-12جρأvv{\displaystyle \mathbf {F_{D}} =-{\tfrac {1}{2}}c\rho A\,v\,\mathbf {v} }
أين:

حالات خاصة

على الرغم من أن الحالة العامة لمقذوف مع مقاومة نيوتن لا يمكن حلها تحليليًا، إلا أنه يمكن حل بعض الحالات الخاصة. هنا نرمز إلى السرعة النهائية في السقوط الحر بـv=ز/μ{\textstyle v_{\infty }={\sqrt {g/\mu }}}وثابت زمن الاستقرار المميزتو1/(vμ)=1/زμ{\textstyle t_{f}\equiv 1/(v_{\infty }\mu )=1/{\sqrt {g\mu }}}(بعدز{\displaystyle g}[م/ث² ] ،μ{\displaystyle \mu }[1/م])

  • حركة شبه أفقية : في حالة كون الحركة شبه أفقية،|vx||vy|{\displaystyle |v_{x}|\gg |v_{y}|}مثل الرصاصة الطائرة. يكون لمركبة السرعة الرأسية تأثير ضئيل للغاية على الحركة الأفقية. في هذه الحالة: [ 14 ]
v˙x(ت)=-μvx2(ت){\displaystyle {\dot {v}}_{x}(t)=-\mu \,v_{x}^{2}(t)}
vx(ت)=11/vx،0+μت{\displaystyle v_{x}(t)={\frac {1}{1/v_{x,0}+\mu \,t}}}
x(ت)=1μln(1+μvx،0ت){\displaystyle x(t)={\frac {1}{\mu }}\ln(1+\mu \,v_{x,0}\cdot t)}
وينطبق النمط نفسه على الحركة مع الاحتكاك على طول خط في أي اتجاه، عندما تكون الجاذبية ضئيلة (صغيرة نسبيًا).ز{\displaystyle g}وينطبق ذلك أيضًا عندما يتم منع الحركة الرأسية، كما هو الحال بالنسبة لسيارة متحركة مع إيقاف تشغيل محركها.
  • الحركة الرأسية للأعلى : [ 14 ]
v˙y(ت)=-ز-μvy2(ت){\displaystyle {\dot {v}}_{y}(t)=-g-\mu \,v_{y}^{2}(t)}
vy(ت)=vلون برونزيتصهـأك-تتو{\displaystyle v_{y}(t)=v_{\infty }\tan {\frac {t_{\mathrm {peak} }-t}{t_{f}}}}
y(ت)=yصهـأك+1μln(كوستصهـأك-تتو){\displaystyle y(t)=y_{\mathrm {peak} }+{\frac {1}{\mu }}\ln {\bigl (}\cos {\frac {t_{\mathrm {peak} }-t}{t_{f}}}{\bigr )}}
هنا
vزμ{\displaystyle v_{\infty }\equiv {\sqrt {\frac {g}{\mu }}}}وتو=1μز،{\displaystyle t_{f}={\frac {1}{\sqrt {\mu g}}},}
تصهـأكتودالة الظل العكسيvy،0v=1μزدالة الظل العكسي(μزvy،0)،{\displaystyle t_{\mathrm {peak} }\equiv t_{f}\arctan {\frac {v_{y,0}}{v_{\infty }}}={\frac {1}{\sqrt {\mu g}}}\arctan {\left({\sqrt {\frac {\mu }{g}}}v_{y,0}\right)},}
و
yصهـأك-1μlnكوستصهـأكتو=12μln(1+μزvy،02){\displaystyle y_{\mathrm {peak} }\equiv -{\frac {1}{\mu }}\ln {\cos {\frac {t_{\mathrm {peak} }}{t_{f}}}}={\frac {1}{2\mu }}\ln {{\bigl (}1+{\frac {\mu }{g}}v_{y,0}^{2}{\bigr )}}}
أينvy،0{\displaystyle v_{y,0}}هي السرعة الابتدائية الصاعدة عندت=0{\displaystyle t=0}والوضع الابتدائي هوy(0)=0{\displaystyle y(0)=0}.
لا يمكن للمقذوف أن يرتفع لمسافة أطول منترأناsهـ=π2تو{\displaystyle t_{\mathrm {rise} }={\frac {\pi }{2}}t_{f}} في الاتجاه الرأسي، عندما يصل إلى الذروة (0 م، ذروة y ) عند 0 م/ث.
  • الحركة الرأسية إلى الأسفل : [ 14 ]
v˙y(ت)=-ز+μvy2(ت){\displaystyle {\dot {v}}_{y}(t)=-g+\mu \,v_{y}^{2}(t)}
vy(ت)=-vtanhت-تصهـأكتو{\displaystyle v_{y}(t)=-v_{\infty }\tanh {\frac {t-t_{\mathrm {peak} }}{t_{f}}}}
y(ت)=yصهـأك-1μln(ضرب بالعصا(ت-تصهـأكتو)){\displaystyle y(t)=y_{\mathrm {peak} }-{\frac {1}{\mu }}\ln(\cosh {\biggl (}{\frac {t-t_{\mathrm {peak} }}{t_{f}}}{\biggr )})}باستخدام الدوال الزائدية
بعد مرور بعض الوقتتو{\displaystyle t_{f}}عند y=0، يصل المقذوف إلى سرعة نهائية تقريبًا.-v{\displaystyle -v_{\infty }}.

الحل العددي

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

ددت(xyvxvy)=(vxvy-μvxvx2+vy2-ز-μvyvx2+vy2){\displaystyle {\frac {\mathrm {d} }{\mathrm {d} t}}{\begin{pmatrix}x\\y\\v_{x}\\v_{y}\end{pmatrix}}={\begin{pmatrix}v_{x}\\v_{y}\\-\mu \,v_{x}{\sqrt {v_{x}^{2}+v_{y}^{2}}}\\-g-\mu \,v_{y}{\sqrt {v_{x}^{2}+v_{y}^{2}}}\end{pmatrix}}}.

يسمح هذا النهج أيضًا بإضافة تأثيرات معامل السحب المعتمد على السرعة، وكثافة الهواء المعتمدة على الارتفاع (في المنتج).ج(v)ρ(y){\displaystyle c(v)\rho (y)}) ومجال الجاذبية المعتمد على الموقعز(y)=ز0/(1+y/R)2{\textstyle g(y)=g_{0}/(1+y/R)^{2}}(متىyR:زز0/(1+2y/R)ز0(1-2y/R){\textstyle y\ll R:g\lesssim g_{0}/(1+2y/R)\approx g_{0}(1-2y/R)}(انخفاض خطي).

مسار مرتفع

مسارات إطلاق الصواريخ الباليستية الكورية الشمالية هواسونغ-14 ، وهواسونغ-15 ، وهواسونغ-17

يُعدّ المسار المرتفع حالة خاصة من المسارات الباليستية للصواريخ، وهو مسارٌ يكون فيه أوج الصاروخ أعلى من مسار الطاقة الأدنى للوصول إلى نفس المدى. بعبارة أخرى، يسافر الصاروخ على ارتفاع أعلى، وبالتالي يستهلك طاقة أكبر للوصول إلى نقطة الهبوط نفسها. قد يُلجأ إلى هذا المسار لأسباب مختلفة، مثل زيادة المسافة إلى الأفق لتحسين مدى الرؤية/الاتصال، أو لتغيير زاوية اصطدام الصاروخ بالأرض. تُستخدم المسارات المرتفعة أحيانًا في كلٍّ من هندسة الصواريخ ورحلات الفضاء . [ 15 ]

حركة المقذوفات على نطاق كوكبي

مسار المقذوف حول كوكب، مقارنة بالحركة في مجال جاذبية منتظم

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

يجب تعديل معلمات المسار من قيم مجال الجاذبية المنتظم المذكورة أعلاه. يُفترض أن نصف قطر الأرض هو R ، و g هي جاذبية السطح القياسية.v~:=v/Rز{\displaystyle {\tilde {v}}:=v/{\sqrt {Rg}}}[ 16 ]​

المدى الكلي بين لحظة الإطلاق ولحظة الاصطدام:

د=v2الخطيئة(2θ)ز/1-(2-v~2)v~2كوس2θ{\displaystyle d={\frac {v^{2}\sin(2\theta )}{g}}{\Big /}{\sqrt {1-\left(2-{\tilde {v}}^{2}\right){\tilde {v}}^{2}\cos ^{2}\theta }}} (حيث زاوية الإطلاق) θ=12أركوس(v~2/(2-v~2)){\displaystyle \theta ={\tfrac {1}{2}}\arccos \left({\tilde {v}}^{2}/(2-{\tilde {v}}^{2})\right)})

أقصى مدى للمقذوف عند زاوية الإطلاق المثلى θ=45 درجة :

دمأx=v2ز/(1-12v~2){\displaystyle d_{\mathrm {max} }={\frac {v^{2}}{g}}{\big /}\left(1-{\tfrac {1}{2}}{\tilde {v}}^{2}\right)}     معv<Rز{\displaystyle v<{\sqrt {Rg}}}، أول سرعة كونية

أقصى ارتفاع يمكن أن يصل إليه المقذوف فوق سطح الكوكب :

ح=v2الخطيئة2θز/(1-v~2+1-(2-v~2)v~2كوس2θ){\displaystyle h={\frac {v^{2}\sin ^{2}\theta }{g}}{\Big /}\left(1-{\tilde {v}}^{2}+{\sqrt {1-\left(2-{\tilde {v}}^{2}\right){\tilde {v}}^{2}\cos ^{2}\theta }}\right)}

أقصى ارتفاع للمقذوف عند إطلاقه عمودياً (θ=90{\displaystyle \theta =90^{\circ }}):

حمأx=v22ز/(1-12v~2){\displaystyle h_{\mathrm {max} }={\frac {v^{2}}{2g}}{\big /}\left(1-{\tfrac {1}{2}}{\tilde {v}}^{2}\right)}     معv<2Rز{\displaystyle v<{\sqrt {2Rg}}}، السرعة الكونية الثانية أو سرعة الإفلات ،

مدة الرحلة:

ت=2vالخطيئةθز12-v~2(1+12-v~2v~الخطيئةθدالة الجيب العكسية2-v~2v~الخطيئةθ1-(2-v~2)v~2كوس2θ){\displaystyle t={\frac {2v\sin \theta }{g}}\cdot {\frac {1}{2-{\tilde {v}}^{2}}}\left(1+{\frac {1}{{\sqrt {2-{\tilde {v}}^{2}}}\,{\tilde {v}}\sin \theta }}\arcsin {\frac {{\sqrt {2-{\tilde {v}}^{2}}}\,{\tilde {v}}\sin \theta }{\sqrt {1-\left(2-{\tilde {v}}^{2}\right){\tilde {v}}^{2}\cos ^{2}\theta }}}\right)}

انظر أيضاً

ملحوظات

  1. g هي تسارع الجاذبية الأرضية .9.81م/s2{\displaystyle 9.81\,\mathrm {m/s^{2}} }(بالقرب من سطح الأرض).
  2. يتناقص عندما يتحرك الجسم لأعلى، ويزداد عندما يتحرك لأسفل
  3. 2الخطيئة(α)كوس(α)=الخطيئة(2α){\displaystyle 2\cdot \sin(\alpha )\cdot \cos(\alpha )=\sin(2\alpha )}

مراجع

  1. غاليليو غاليلي، علمان جديدان ، ليدن، 1638، ص 249
  2. نولتي، ديفيد د. ، غاليليو غير المقيد (مطبعة جامعة أكسفورد، 2018) ص 39-63.
  3. 1 2 بوركو، ليور م.؛ برايس، ريتشارد هـ. (2005-06-01). "المسار الباليستي: قطع مكافئ، قطع ناقص، أم ماذا؟" . المجلة الأمريكية للفيزياء . 73 (6): 516-520 . arXiv : physics/0310049 . doi : 10.1119/1.1866097 . ISSN 0002-9505 . 
  4. ستيوارت، جيمس؛ كليج، دان؛ واتسون، سليم (2021). حساب التفاضل والتكامل: الدوال المتسامية المبكرة ( الطبعة التاسعة). بوسطن، ماساتشوستس: سينجج. ص 919. ISBN   978-1-337-61392-7.
  5. موبس، ويليام؛ لينغ، صموئيل جيه؛ ساني، جيف (19-09-2016). "4.3 حركة المقذوفات - فيزياء الجامعة، المجلد 1 | أوبن ستاكس" . openstax.org . تم الاطلاع عليه بتاريخ 04-06-2026 .
  6. غالانت، جوزيف (2012). ممارسة الفيزياء باستخدام دفتر الملاحظات العلمي: منهج لحل المشكلات . جون وايلي وأولاده . ص 132. ISBN  978-1-119-94194-1.مقتطف من الصفحة 132. لاحظ أنه تم استبدال قيمة yy 0 في المصدر بقيمة y 0 في المقالة.
  7. تاتوم (2019). الميكانيكا الكلاسيكية (ملف PDF) . ص. الفصل 7. 
  8. مع V=v 2 ، V 2 - 2gy.V - (gx) 2 = 0:V1،2=(2زy±4ز2y2+4(زx)2)/2{\displaystyle V_{1,2}=(2gy\pm {\sqrt {4g^{2}y^{2}+4(gx)^{2}}})/2}
  9. تأن(θ)=y/x+(y/x)2+1{\textstyle tan(\theta )=y/x+{\sqrt {(y/x)^{2}+1}}}أينy/x=sأنان(α)/جos(α):sأنان2α+جos2α=1=(تأن2α+1)جos2α{\textstyle y/x=sin(\alpha )/cos(\alpha ):sin^{2}\alpha +cos^{2}\alpha =1=(tan^{2}\alpha +1)cos^{2}\alpha }
  10. 1/لون برونزي(π/2-θ)=كوس(π/2-θ)/الخطيئة(π/2-θ)=الخطيئةθ/كوسθ{\displaystyle 1/\tan(\pi /2-\theta )=\cos(\pi /2-\theta )/\sin(\pi /2-\theta )=\sin \theta /\cos \theta }=(الخطيئةα+1)/كوسα=لون برونزيα+1/كوسα{\displaystyle =(\sin \alpha +1)/\cos \alpha =\tan \alpha +1/\cos \alpha }=لون برونزيθ{\displaystyle =\tan \theta }
  11. ستيفن تي. ثورنتون؛ جيري بي. ماريون (2007). الديناميكا الكلاسيكية للجسيمات والأنظمة . بروكس/كول. ص 59. ISBN  978-0-495-55610-7.
  12. أتام ب. آريا؛ أتام باركاش آريا (سبتمبر 1997). مقدمة في الميكانيكا الكلاسيكية . برنتيس هول إنترناشونال. ص 227. ISBN  978-0-13-906686-3.
  13. ريجينالد كريستيان، برناردو؛ خوسيه بيريكو، إسغيرا؛ جازمين داي، فاليخوس؛ جيف جيرارد، كندا (2015). "حركة المقذوفات المتأثرة بالرياح". المجلة الأوروبية للفيزياء . 36 (2) 025016. Bibcode : 2015EJPh...36b5016B . doi : 10.1088/0143-0807/36/2/025016 . S2CID 119601402 . 
  14. 1 2 3 والتر غرينر (2004). الميكانيكا الكلاسيكية: الجسيمات النقطية والنسبية . سبرينغر ساينس آند بيزنس ميديا. ص 181. ISBN  0-387-95586-0.
  15. الدفاع الصاروخي الباليستي، مسرد المصطلحات، الإصدار 3.0 ، وزارة الدفاع الأمريكية ، يونيو 1997.
  16. lgeorge2. "الفصل 2 - هندسة المدار" .{{cite web}}: صيانة CS1: الأسماء الرقمية: قائمة المؤلفين ( رابط )