معادله دیفرانسیل حرکت – از صفر تا صد (+ دانلود فیلم آموزش رایگان)

مجید عوض زاده
مجید عوض زاده
۱۴۰۲/۰۲/۳۰
۳۵ دقیقه
۵۱۹۵
آموزش متنی جامع
آموزش ویدئویی
در مطالب گذشته وبلاگ فرادرس مفاهیمی از فیزیک کلاسیک همچون معادله ناویر-استوکس، معادلات مومنتوم (تکانه) و معادله حرکت یک جسم در حال ارتعاش را توضیح دادیم. شاید جالب باشد بدانید که منشا تمامی این مفاهیم، قوانین نیوتن است. به طور دقیق‌تر می‌توان گفت شکل کمی این روابط بیشترین ارتباط را با قانون دوم نیوتن دارد. از این رو در این مطلب قصد داریم تا نحوه بدست آوردن معادله دیفرانسیل حرکت یک سیستم را با استفاده از این قانون، توضیح دهیم.
معادله دیفرانسیل حرکت – از صفر تا صد (+ دانلود فیلم آموزش رایگان)
997696

فیلم آموزشی معادله دیفرانسیل حرکت

دانلود ویدیو

معادله دیفرانسیل حرکت

قانون دوم نیوتن در حالت کلی ارتباط بین نیروی FF وارد به جسمی به جرم mm را که با شتاب aa در حرکت است، توضیح می‌دهد. این قانون می‌گوید شتاب جرم به صورت خطی نسبت به نیرو افزایش می‌یابد. به صورتی دقیق‌تر فرض کنید نیروی FF به جسمی به جرم mm وارد می‌شود. در این صورت می‌توان رابطه زیر را برای این جرم نوشت:

a→=F→m    ,    F→=ma→=md2r→dt2\large {\overrightarrow{a} = \frac { \overrightarrow { F } } { m } \; \; \text{,}\;\;}\kern-0.3pt{\overrightarrow{F} = m \overrightarrow { a } = m\frac { { { d ^ 2 } \overrightarrow{r } } }{ { d { t ^ 2} } } }

توجه داشته باشید که عبارت فوق، معادله‌ای برداری را نشان می‌دهد؛ به طور دقیق‌تر، دو کمیت rr و FF بردار هستند. همچنین این رابطه، برای سیستمی با جرم ثابت صدق می‌کند. در حقیقت در مواردی همچون مکانیک نسبیتی که جرم وابسته به سرعت است، باید قانون دوم نیوتن را به صورت تغییرات تکانه بر حسب زمان و به صورت زیر بیان کرد:

F→=dp→dt\large \overrightarrow { F } = \frac { { d \overrightarrow {p}}} { { d t } }

به pp، تکانه یا مومنتوم گفته می‌شود. در حالتی کلی نیروی FF وابسته به دستگاه مختصات انتخاب شده برای جسم است. در ادامه نیرو به صورت متغیر نسبت به زمان در نظر گرفته شده است.

نیروی متغیر با زمان: F=F(t)\large F = F ( t )

حالتی را در نظر بگیرید که در آن جسم حرکتی تک‌بعدی را تجربه می‌کند. در این صورت قانون دوم نیوتن را می‌توان بر حسب معادله دیفرانسیل مرتبه دوم بیان کرد. در ادامه این معادله ارائه شده است.

md2xdt2=F(t)\large m \frac{{ { d^2 }x } }{{d { t^ 2 }} } = F \left( t \right)

با انتگرال‌گیری از رابطه فوق، سرعت جسم مطابق با رابطه زیر بدست خواهد آمد.

v(t)=v0+1m∫0tF(τ)dτ\large { v \left( t \right) }={ {v_0} + \frac {1} { m }\int\limits_0^t {F\left( \tau \right)d\tau } }

برای نمونه حالتی را در نظر بگیرید که در آن جسم با سرعت اولیه v=0v=0 در زمان t=0t=0 شروع به حرکت کرده است. با انتگرال‌گیری دوباره از عبارت فوق، معادله حرکت یا همان تابع x(t)x(t) به صورت زیر بدست می‌آید.

x(t)=x0+∫0tv(τ)dτ\large x \left ( t \right ) = { x _ 0 } + \int \limits _ 0 ^ t { v \left ( \tau \right ) d \tau }

در رابطه فوق x0x_0 مکان اولیه بردار جابجایی و ττ متغیر انتگرال‌گیری است.

وابستگی نیرو به سرعت: F=F(v)\large F = F ( v )

زمانی که جسمی در یک سیالِ گازی یا مایع به حرکت در می‌آید، نیرویی را از جانب سیال حس می‌کند که مهندسان آن را درگ یا پسا می‌نامند. در سرعت‌های پایین اندازه این نیرو وابسته به سرعت جسم در سیال است. اندازه این نیرو را می‌توان در حالت کلی به صورت زیر بیان کرد:

F→=–kv→\large \overrightarrow { F } = – k \overrightarrow{v}

ضریب kk وابسته به ویسکوزیته η\eta است. در حالتی ویژه اگر جسم به صورت کروی شکل باشد، آن‌گاه نیروی درگ مطابق با قانون استوکس، به صورت زیر قابل محاسبه است.

F→=–6πηRv→\large \overrightarrow { F } = – 6 \pi \eta R \overrightarrow { v }

در رابطه فوق RR، شعاع توپ و ηη ویسکوزیته محیط یا همان سیالی است که جسم در آن در حال حرکت است. در حالت حرکتِ تک‌بعدی در سیال، قانون دوم نیوتن را می‌توان به صورت دیفرانسیلی نیز بیان کرد:

md2xdt2=mdvdt=–kv\large m \frac { { { d ^ 2 } x }} {{d {t ^ 2 }}} = m \frac { { d v } } { { d t }} = – k v

با انتگرال‌گیری از رابطه فوق و فرض سرعت اولیه v0v_0 در زمان t0t_0، داریم:

dvv=–kmdt    ⇒∫v0vduu=–km∫0tdτ\large { \frac{{ d v } }{ v} = – \frac { k } { m } d t \;\;}\Rightarrow { \int \limits_{{ v _0 } } ^v { \frac{{ d u } }{ u } } = – \frac { k} { m}\int\limits_0^t { d \tau } }

توجه داشته باشید که در عبارت فوق، uu و vv متغیر‌های انتگرال‌گیری هستند. تغییرات سرعت از vv تا v0v_0 و تغییرات زمانی از 00 تا tt در نظر گرفته شده‌اند. بنابراین حاصل انتگرال فوق به صورت زیر در می‌آید.

ln⁡v–ln⁡v0=–kmt    ⇒ln⁡vv0=–kmt    ⇒v(t)=v0e–kmt\large \begin{align*} { \ln v – \ln { v _ 0 } = – \frac { k }{ m } t \;\;}\Rightarrow {\ln \frac{v}{{{v_0}}} = – \frac { k } { m } t \;\;}\Rightarrow { v \left( t \right) = { v _ 0} { e ^{ – {\large\frac { k } { m } \normalsize}t}} } \end{align*}

از این رو اگر نیروی درگ متناسب با سرعت جسم باشد، سرعت جسم به صورت نمایی کم خواهد شد. اگر برای بار دوم از عبارت بدست آمده، انتگرال بگیریم، معادله حرکت جسم به صورت زیر بدست خواهد آمد.

x(t)=x0+∫0tv(τ)dτ=x0+∫0tv0e–kmτdτ=x0–mv0k(e–kmt–1)=x0+mv0k(1–e–kmt)\large \begin {align*} x \left( t \right) & = {x_0} + \int\limits_0^t {v\left( \tau \right)d\tau } \\ & = { { x _ 0 } + \int \limits _ 0 ^ t { { v _0 } { e ^ { – {\Large \frac{k}{m}\normalsize}\tau }}d\tau } } \\ & = { { x _ 0 } – \frac{{m{v_0}}}{k}\left( {{e^{ – {\Large\frac{k}{m}\normalsize}t}} – 1} \right) } \\ & = {{x_0} + \frac{{m{v_0}}}{k}\left( {1 – { e ^ { – { \large\frac { k } {m } \normalsize} t } } } \right) } \end {align*}

رابطه فوق نشان می‌دهد که وابستگی بین کمیت‌های xx و tt در حرکت جسم تا ایستادن، وابسته به مومنتوم اولیه جسم یا همان mv0mv_0 است. جالب است بدانید با افزایش سرعت جسم، فیزیک فرآیند نیز تغییر می‌کند. در حقیقت در سرعت بالا انرژی جنبشی جسم نه تنها صرف اصطکاک بین لایه‌ها می‌شود، بلکه منجر به حرکت لایه‌ها در جلوی جسم نیز خواهد شد. در این حالت، نیروی وارد به جسم وابسته به توان دوم سرعت خواهد بود. در حقیقت اندازه نیرو در این حالت برابر است با:

F=–μρSv2\large F = – \mu \rho S { v ^ 2 }

در رابطه فوق، μμ، ضریب نسبت، SS، سطح مقطع جسمِ در حال حرکت و ρρ چگالی سیال است. رفتاری که در بالا توضیح داده شده، رفتاری غیر خطی محسوب می‌شود. این رفتار در شرایطی رخ می‌دهد که عدد رینولدز کمتر از  100100 باشد.

Re→=ρvLη>100\large \overrightarrow {\text{Re}} = \frac {{\rho vL}} {\eta } > 100

در رابطه فوق ReRe نشان دهنده عدد بی‌بعد رینولدز است. هم‌چنین ηη، ویسکوزیته و LL، طول مشخصه جسم است. برای نمونه در حالتی که توپی در هوا در حال حرکت باشد، طول مشخصه برابر با شعاع توپ در نظر گرفته می‌شود.

حال می‌خواهیم معادله حرکت را در حالتی بدست آوریم که نیروی وارد به جسم با توان دوم سرعت تغییر می‌کند. در این حالت معادله دیفرانسیل حاکم بر حرکت جسم به صورت زیر قابل بیان است.

md2xdt2=mdvdt=–μρSv2\large { m \frac{ { { d^ 2 } x } } {{d { t ^ 2 } }} = m\frac { {d v } }{ { d t }} }={ – \mu \rho S { v ^ 2}}

با انتگرال‌گیری از رابطه فوق داریم:

md2xdt2=mdvdt=–μρSv2\large { m \frac { { {d ^ 2 } x } } {{ d { t ^2 } } } = m \frac { { d v} } { { d t } } }={ – \mu \rho S { v ^ 2 } }

توجه داشته باشید که در رابطه بالا uu و ττ متغیر‌های انتگرال‌گیری هستند. سرعتِ uu از مقدار v0v_0 شروع شده و تا مقدار vv کاهش می‌یابد. هم‌چنین این فرآیند از زمان 00 تا  tt رخ می‌دهد. بنابراین با قرار دادن این مقادیر به عنوان بازه‌های انتگرال‌، به رابطه زیر خواهیم رسید.

−(1v–1v0)=–μρSmt⇒1v=1v0+μρSmt⇒v(t)=11v0+μρSmt=v01+μρSv0mt\large \begin {align*} -\left( {\frac{1}{v} – \frac { 1} { { { v _ 0 }} } } \right) & = – \frac{{\mu \rho S } }{ m } t \\ & \Rightarrow \frac{1}{v} = \frac { 1}{ { { v _0 }} } + \frac { { \mu \rho S } } { m } t \\ & \Rightarrow \large v \left( t \right) = \frac { 1} { { \frac{1} { {{ v _ 0 }}} + \frac { { \mu \rho S} } {m } t }} = \frac { { {v _ 0} } } {{1 + \frac { { \mu \rho S{v_ 0 } } } { m } t } } \end {align*}

با انتگرال‌گیریِ دوباره از رابطه فوق، شکل نهایی معادله حرکت به صورت زیر بدست خواهد آمد.

$$ \large \begin {align*} \require{cancel} x\left( t \right) & = \int\limits_0^t {\frac{ { {v _ 0 } } } { {1 + \frac{{\mu \rho S {v _ 0 } }} {m }\tau }}d\tau } \\ & = \int\limits_0^t {\frac{{\cancel { v _0 } } } { { 1 + \frac { { \mu \rho S {v_0 } } } {m } \tau }}\frac{{d\left( {1 + \frac{{\mu \rho S { v _ 0} } } { m } \tau } \right) } }{ {\frac{{\mu \rho S\cancel{v_0}}}{m}}}} \\ & = {\frac { m }{ { \mu \rho S } }\int\limits _ 0 ^ t {\frac{{d\left( {1 + \frac { {\mu \rho S { v_ 0 } } } { m }\tau } \right ) } } { {1 + \frac{{\mu \rho S { v _0 }} } { m } \tau } } } } \\ & = {\frac{m}{{\mu \rho S } }\left[ {\left. {\ln \left( {1 + \frac { { \mu \rho S { v _0 } } }{m}\tau } \right)} \right|_0^t} \right] } \\ & = { \frac { m } { { \mu \rho S } } \ln \left( {1 + \frac{{\mu \rho S { v _ 0 } } } {m } t } \right) } \end {align*} $$

مطالب مجله فرادرس را در جستجوی گوگل سریع‌تر ببینید

بر اساس رای ۱۶ نفر
آیا این مطلب برای شما مفید بود؟
اگر پرسشی درباره این مطلب دارید، آن را با ما مطرح کنید.
منابع:
Math24
PDF
مطالب مرتبط
نظر شما چیست؟

نشانی ایمیل شما منتشر نخواهد شد. بخش‌های موردنیاز علامت‌گذاری شده‌اند *