نمای شبیه متلب از مسیر حرکت
لینکهای فیزیکی ضخیمتر همراه با مسیر بسیار نازک نقطه انتهایی.
نمایش زنده دینامیک آشوبی
شبیهساز آشوب آونگ دوگانه
این شبیهساز تحت وب نشان میدهد که چگونه یک سامانه مکانیکی کاملاً قطعی میتواند در عمل غیرقابل پیشبینی شود، وقتی حرکت آن نسبت به شرایط اولیه بسیار حساس باشد.
نمای مکانیزم شبیه CAD
مقایسه آونگ اصلی و آونگ اندکی آشفتهشده برای نمایش اثر پروانهای.
زمان شبیهسازی
00:00.00
مدل اصلی
\[ \begin{alignedat}{2} \dot{\theta}_1 &= \omega_1,\qquad& \dot{\omega}_1 &= f_1(\theta_1,\theta_2,\omega_1,\omega_2),\\ \dot{\theta}_2 &= \omega_2,\qquad& \dot{\omega}_2 &= f_2(\theta_1,\theta_2,\omega_1,\omega_2) \end{alignedat} \]
آماده است. شرایط اولیه را وارد کنید و دکمه اجرای شبیهسازی را بزنید.
آشوب چیست؟
آشوب به معنی تصادفی بودن نیست. یک سامانه آشوبی همچنان قطعی است؛ یعنی آینده آن توسط معادلات مشخص و دقیق تولید میشود. با این حال، اختلافهای بسیار کوچک در شرایط اولیه میتوانند با گذشت زمان به اختلافهای بزرگ تبدیل شوند. از آنجا که هر اندازهگیری واقعی دقت محدودی دارد، پیشبینی بلندمدت چنین سامانهای حتی با دانستن معادلات حاکم، در عمل بسیار دشوار یا ناممکن میشود.
آونگ دوگانه یکی از نمونههای کلاسیک آشوب است، زیرا فقط دو لینک صلب و دو زاویه دارد، اما کوپل غیرخطی میان این دو لینک آنقدر قوی است که میتواند حرکت نامنظم، غیرتناوبی و بسیار حساس نسبت به شرایط اولیه ایجاد کند.
اثر پروانهای در این شبیهساز
اثر پروانهای یعنی یک اغتشاش بسیار کوچک در حالت اولیه یک سامانه دینامیکی غیرخطی میتواند در نهایت اختلافی بزرگ و قابل مشاهده ایجاد کند. در این شبیهساز، مسیر اصلی و مسیر آشفتهشده تقریباً از شرایط اولیه یکسان شروع میشوند. پنل CAD مانند هر دو حالت آونگ دوگانه را نمایش میدهد، در حالی که پنل شبیه متلب مسیر نقطه انتهایی لینک دوم را برجسته میکند.
در ابتدای حرکت، دو پاسخ ممکن است تقریباً یکسان به نظر برسند. اما چون آونگ دوگانه یک سامانه مکانیکی کوپلشده و شدیداً غیرخطی است، همان اختلاف کوچک بهتدریج تقویت میشود. پس از زمان کافی، مسیرهای نقطه انتهایی بهوضوح از هم جدا میشوند و حرکت دیگر شبیه یک نوسان ساده یا هارمونیک نیست. در عوض، نقطه انتهایی ناحیهای نامنظم از صفحه را پر میکند که نشانه بصری رفتار آشوبی است.
این حساسیت فقط به اغتشاش کوچک مصنوعی اضافهشده توسط شبیهساز محدود نیست. حتی تغییرهای کوچک در شرایط اولیه واقعی یا پارامترهای قابل تنظیم نیز میتوانند حرکت آینده را کاملاً تغییر دهند. برای نمونه، تغییر جزئی زاویه لینک اول، زاویه لینک دوم، سرعت زاویهای اولیه هر لینک، طول لینکها یا جرم لینکها میتواند جملههای شتاب در معادلات حاکم را تغییر دهد. با ادامه شبیهسازی، این تغییرهای کوچک جمع میشوند و ممکن است به توالی کاملاً متفاوتی از تاب خوردن، چرخش، وارونگی و مسیرهای نقطه انتهایی منجر شوند.
بنابراین آونگ دوگانه یک نمونه بسیار مناسب از آشوب قطعی است. حرکت تصادفی نیست؛ بلکه از معادلات دیفرانسیل قطعی تولید میشود. با این حال، چون حالت آینده به شرایط اولیه و پارامترهای سیستم بسیار حساس است، پیشبینی بلندمدت آن بسیار دشوار میشود. دو شبیهسازی که تقریباً از یک نقطه شروع میکنند، پس از زمان کافی میتوانند به حرکتهایی برسند که از نظر فیزیکی کاملاً متفاوت به نظر میرسند.
در عمل، تغییر پارامترهایی مانند طول لینک، جرم لینک، زاویه اولیه یا سرعت زاویهای اولیه فقط باعث کندتر یا سریعتر شدن همان حرکت قبلی نمیشود؛ بلکه میتواند رفتار کیفی سیستم را تغییر دهد. یک مجموعه پارامتر ممکن است برای مدتی حرکت نرم و نسبتاً منظم تولید کند، در حالی که مجموعهای بسیار نزدیک به آن ممکن است به سرعت حرکتهای چرخشی، حلقهای و مسیرهای کاملاً نامنظم ایجاد کند.
مختصات آونگ دوگانه
فرض کنید \(\theta_1\) و \(\theta_2\) موقعیتهای زاویهای مطلق لینک اول و لینک دوم باشند که از راستای قائم اندازهگیری میشوند. مختصات نقاط انتهایی به صورت زیر است:
\[ x_1 = L_1 \sin(\theta_1),\qquad y_1 = L_1 \cos(\theta_1) \] \[ x_2 = x_1 + L_2 \sin(\theta_2),\qquad y_2 = y_1 + L_2 \cos(\theta_2). \]
نمایش گرافیکی، این مختصات فیزیکی را به مختصات بوم رسم تبدیل میکند و نقطه اتصال اصلی را نزدیک مرکز بالایی هر پنل قرار میدهد. در پنل سمت چپ، مسیر نقطه انتهایی با خط بسیار نازک رسم میشود و خود دو لینک نیز با خط ضخیمتر نمایش داده میشوند تا حرکت، معنی فیزیکی خود را حفظ کند.
معادلات دینامیکی استفادهشده در حلگر جاوااسکریپت
برای جرمهای \(m_1,m_2\)، طولهای \(L_1,L_2\)، گرانش \(g\)، زاویههای \(\theta_1,\theta_2\)، سرعتهای زاویهای \(\omega_1,\omega_2\)، و \(\Delta = \theta_1 - \theta_2\)، حلگر عددی مدل کوپلشده استاندارد آونگ دوگانه را محاسبه میکند:
\[ \dot{\omega}_1 = \frac{-g(2m_1+m_2)\sin\theta_1 - m_2g\sin(\theta_1-2\theta_2) - 2m_2\sin\Delta\left(\omega_2^2L_2 + \omega_1^2L_1\cos\Delta\right)} {L_1\left(2m_1+m_2-m_2\cos(2\Delta)\right)} \]
\[ \dot{\omega}_2 = \frac{2\sin\Delta\left(\omega_1^2L_1(m_1+m_2)+g(m_1+m_2)\cos\theta_1+\omega_2^2L_2m_2\cos\Delta\right)} {L_2\left(2m_1+m_2-m_2\cos(2\Delta)\right)}. \]
شبیهسازی این معادلات را با روش رانگ-کوتای مرتبه چهارم انتگرالگیری میکند. این روش نسبت به روش ساده اویلر حرکت نرمتری تولید میکند و واگرایی آشوبی را بدون تخریب سریع ظاهر مکانیکی انیمیشن، قابل مشاهدهتر میسازد.