آونگ (مکانیک)

[چپ|300px|بندانگشتی|پویانمایی یک آونگ که بردارهای [معادلات حرکت|سرعت و شتاب را نشان می‌دهد.]]
**آونگ** یا **پاندول** (به انگلیسی: Pendulum) جسمی است که از یک تکیه‌گاه ثابت آویزان شده و آزادانه تحت تأثیر نیروی گرانش به عقب و جلو نوسان می‌کند. هنگامی که یک آونگ از حالت تعادل یا استراحت خود به سمت جانبی جابجا می‌شود، تحت تأثیر نیروی بازگرداننده‌ای ناشی از گرانش قرار می‌گیرد که آن را به سمت موقعیت تعادل شتاب می‌دهد. وقتی رها می‌شود، نیروی بازگرداننده که بر جرم آونگ اثر می‌گذارد، باعث می‌شود که حول نقطه تعادل نوسان کرده و به عقب و جلو تاب بخورد. ریاضیات آونگها به‌طور کلی بسیار پیچیده است. با این حال، می‌توان فرضیات ساده‌کننده‌ای انجام داد که در مورد آونگ ساده اجازه می‌دهد معادلات حرکت برای نوسانات با زاویه کوچک به‌صورت تحلیلی حل شوند.

## آونگ گرانشی ساده
یک *آونگ گرانشی ساده* یک مدل ریاضی ایده‌آل از یک آونگ واقعی است. این مدل شامل وزنه‌ای (یا باب) در انتهای یک طناب بدون جرم است که از یک محور آویزان شده و فاقد اصطکاک است. از آنجا که در این مدل اتلاف انرژی اصطکاکی وجود ندارد، وقتی جابجایی اولیه به آن داده شود، با دامنه ثابت به عقب و جلو نوسان می‌کند. این مدل بر اساس فرضیات زیر است:
* میله یا طناب بدون جرم، غیرقابل کشش است و همواره تحت کشش باقی می‌ماند.
* وزنه (باب) یک جرم نقطه‌ای است.
* حرکت در دو بُعد رخ می‌دهد.
* حرکت انرژی خود را بر اثر اصطکاک خارجی یا مقاومت هوا از دست نمی‌دهد.
* میدان گرانشی یکنواخت است.
* تکیه‌گاه بی‌حرکت است.

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

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

### اثبات‌ها

بندانگشتی|'''شکل ۱.''' نمودار نیروی یک آونگ گرانشی ساده.|چپ|300px
شکل ۱ در سمت چپ را در نظر بگیرید که نیروهای وارد بر یک آونگ ساده را نشان می‌دهد. توجه داشته باشید که مسیر آونگ یک کمان از دایره را طی می‌کند. زاویه  بر حسب رادیان اندازه‌گیری می‌شود و این نکته برای این فرمول حیاتی است. پیکان آبی نیروی گرانش وارد بر وزنه است و پیکان‌های بنفش همان نیرو هستند که به مؤلفه‌های موازی و عمود بر حرکت لحظه‌ای وزنه تجزیه شده‌اند. جهت سرعت لحظه‌ای وزنه همواره در امتداد محور قرمز است که محور مماسی در نظر گرفته می‌شود زیرا جهت آن همیشه مماس بر دایره است. قانون دوم نیوتن را در نظر بگیرید:
F=ma
که در آن  مجموع نیروهای وارد بر جسم،  جرم و  شتاب است. معادله نیوتن را می‌توان تنها برای محور مماسی اعمال کرد. این به این دلیل است که تنها تغییرات سرعت مورد نظر است و وزنه مجبور به ماندن در مسیر دایره‌ای است. پیکان بنفش کوتاه نشان‌دهنده مؤلفه نیروی گرانش در محور مماسی است و می‌توان با استفاده از مثلثات بزرگی آن را تعیین کرد. بنابراین:
\begin{align} F &= -mg\sin\theta = ma, \qquad \text{پس} \\ a &= -g \sin\theta,\end{align}
که در آن  شتاب ناشی از گرانش در نزدیکی سطح زمین است. علامت منفی در سمت راست نشان می‌دهد که  و  همیشه در جهت مخالف یکدیگر هستند. این منطقی است زیرا وقتی آونگ بیشتر به سمت چپ تاب می‌خورد، انتظار می‌رود که به سمت راست شتاب بگیرد.

این شتاب خطی  در امتداد محور قرمز را می‌توان با فرمول‌های طول کمان به تغییرات زاویه  مرتبط کرد؛  طول کمان است:
\begin{align} s &= \ell\theta, \\ v &= \frac{ds}{dt} = \ell\frac{d\theta}{dt}, \\ a &= \frac{d^2s}{dt^2} = \ell\frac{d^2\theta}{dt^2}, \end{align}
بنابراین:
\begin{align} \ell\frac{d^2\theta}{dt^2} &= -g \sin\theta, \\ \frac{d^2\theta}{dt^2} + \frac{g}{\ell} \sin\theta &= 0. \end{align} 

معادله (۱) را می‌توان با استفاده از دو تعریف برای گشتاور به دست آورد.
\boldsymbol{\tau} = \mathbf{r} \times \mathbf{F} = \frac{d\mathbf{L}}{dt}.

ابتدا با تعریف گشتاور روی وزنه آونگ با استفاده از نیروی ناشی از گرانش شروع می‌کنیم.
\boldsymbol{ \tau } = \mathbf{l} \times \mathbf{F}_\mathrm{g} ,
که در آن  بردار طول آونگ و  نیروی ناشی از گرانش است.

فعلاً فقط بزرگی گشتاور روی آونگ را در نظر می‌گیریم.
|\boldsymbol{\tau}| = -mg\ell\sin\theta,
که در آن  جرم آونگ،  شتاب ناشی از گرانش،  طول آونگ و  زاویه بین بردار طول و نیروی ناشی از گرانش است.

سپس تکانه زاویه‌ای را بازنویسی می‌کنیم.
\mathbf{ L } = \mathbf{r} \times \mathbf{p} = m \mathbf{r} \times (\boldsymbol{\omega} \times \mathbf{r}) .
دوباره فقط بزرگی تکانه زاویه‌ای را در نظر بگیرید.
 |\mathbf{L}| = mr^2 \omega = m \ell^2 \frac{d\theta}{dt} .
و مشتق زمانی آن:
\frac {d}{dt}|\mathbf{L}| = m \ell^2 \frac{d^2\theta}{dt^2} ,

سپس می‌توان بزرگی‌ها را با استفاده از  مقایسه کرد:

-mg\ell\sin\theta = m \ell^2 \frac{d^2\theta}{dt^2} ,
بنابراین:
\frac{d^2\theta}{dt^2} + \frac{g}{\ell} \sin\theta = 0,
که همان نتیجه به دست آمده از طریق تحلیل نیرو است.

بندانگشتی|'''شکل ۲.''' مثلثات یک آونگ گرانشی ساده.|چپ|300px

این معادله همچنین می‌تواند از طریق اصل پایستگی انرژی مکانیکی به دست آید: هر جسمی که مسافت عمودی h را سقوط کند، انرژی جنبشی برابر با آنچه در سقوط از دست داده به دست می‌آورد. به عبارت دیگر، انرژی پتانسیل گرانشی به انرژی جنبشی تبدیل می‌شود. تغییر در انرژی پتانسیل برابر است با:
\Delta U = mgh.

تغییر در انرژی جنبشی (جسم از حالت سکون شروع به حرکت کرده) برابر است با:
\Delta K = \tfrac12 mv^2.

از آنجا که هیچ انرژی از دست نمی‌رود، افزایش در یکی باید برابر با کاهش در دیگری باشد:
\tfrac12 mv^2=mgh.

تغییر سرعت برای یک تغییر ارتفاع معین را می‌توان به صورت زیر بیان کرد:
v = \sqrt{2gh}.

با استفاده از فرمول طول کمان در بالا، این معادله را می‌توان بر حسب  بازنویسی کرد:
\begin{align} v = \ell\frac{d\theta}{dt} &= \sqrt{2gh}, \quad \text{پس} \\ \frac{d\theta}{dt} &= \frac{\sqrt{2gh}}{\ell}, \end{align}
که در آن  مسافت عمودی است که آونگ سقوط کرده است. به شکل ۲ نگاه کنید که مثلثات یک آونگ ساده را نشان می‌دهد. اگر آونگ نوسان خود را از زاویه اولیه  آغاز کند، آنگاه ، فاصله عمودی از نقطه اتکا، برابر است با:
y_0 = \ell\cos\theta_0.

به طور مشابه، برای  داریم:
y_1 = \ell\cos\theta.

سپس  تفاضل این دو است:
h = \ell\left(\cos\theta-\cos\theta_0\right).

بر حسب  داریم:

این معادله به عنوان *انتگرال اول حرکت* شناخته می‌شود؛ سرعت را بر حسب موقعیت می‌دهد و شامل یک ثابت انتگرال‌گیری مربوط به جابجایی اولیه () است. سپس، با اعمال قاعده زنجیره‌ای نسبت به زمان مشتق می‌گیریم تا شتاب به دست آید:
\begin{align}
\frac{d}{dt}\frac{d\theta}{dt} &= \frac{d}{dt}\sqrt{\frac{2g} \ell \left(\cos\theta-\cos\theta_0\right)}, \\
\frac{d^2\theta}{dt^2} & = \frac12\frac{-\frac{2g} \ell \sin\theta}{\sqrt{\frac{2g} \ell (\cos\theta-\cos\theta_0)}}\frac{d\theta}{dt} \\
& = \frac12\frac{-\frac{2g}{\ell} \sin\theta}{\sqrt{\frac{2g}{\ell} (\cos\theta-\cos\theta_0)}}\sqrt{\frac{2g}{\ell} (\cos\theta-\cos\theta_0)} = -\frac g \ell \sin\theta, \\
\frac{d^2\theta}{dt^2} &+ \frac{g}{\ell} \sin\theta = 0,
\end{align}

که همان نتیجه به دست آمده از طریق تحلیل نیرو است.

بندانگشتی|چپ|مختصات یک آونگ گرانشی ساده.|300px

 را می‌توان همچنین از طریق مکانیک لاگرانژی به دست آورد. به طور خاص، با استفاده از معادلات اویلر-لاگرانژ (یا معادلات نوع دوم لاگرانژ) با تعیین لاگرانژی سیستم (\mathcal{L})، قیدها (q) و حل دستگاه معادلات زیر:

\frac{d}{dt}\left(\frac{\partial \mathcal{L}}{\partial \dot{q_j}}\right) = \frac{\partial \mathcal{L}}{\partial q_j}.

اگر مبدأ دستگاه مختصات دکارتی به عنوان نقطه آویز (یا لولا) تعریف شود، آنگاه موقعیت وزنه برابر است با:

 x = \ell\sin{\theta},
 y = -\ell\cos{\theta},

و سرعت وزنه، که از طریق مشتق‌گیری مختصات نسبت به زمان محاسبه می‌شود (با استفاده از نشان‌گذاری نقطه برای نشان دادن مشتقات زمانی):

 \dot{x} = \ell\dot{\theta}\cos{\theta},
 \dot{y} = \ell\dot{\theta}\sin{\theta}.

بنابراین، لاگرانژی برابر است با:

\begin{align}
\mathcal{L} &= E_k - E_p \\
&= \frac{1}{2} m v^2 - m g h \\
&= \frac{1}{2} m (\dot{x}^2 + \dot{y}^2) - m g \ell (1 - \cos{\theta}) \\
&= \frac{1}{2} m \ell^2 \dot{\theta}^2 - m g \ell + m g \ell \cos{\theta}.
\end{align}

**معادله اویلر-لاگرانژ** (مفرد است زیرا تنها یک قید q = \theta وجود دارد) بدین صورت است:

\begin{align}
\frac{d}{dt}\left(\frac{\partial \mathcal{L}}{\partial \dot{\theta}}\right) &= \frac{\partial \mathcal{L}}{\partial \theta} \\
\frac{d}{dt}(m \ell^2 \dot{\theta}) &= -m g \ell \sin{\theta} \\
m \ell^2 \ddot{\theta} &= -m g \ell \sin{\theta} \\
\ddot{\theta} &= -\frac{g}{\ell}\sin{\theta}. \\
\end{align}

که می‌توان آن را بازآرایی کرد تا با  که از طریق تحلیل نیرو به دست آمد، مطابقت داشته باشد:

 \frac{d^2 \theta}{d t^2} + \frac{g}{\ell}\sin{\theta} = 0.

استخراج از طریق مکانیک لاگرانژی، اگرچه برای یک آونگ ساده زیاده‌روی به نظر می‌رسد، اما برای سیستم‌های پیچیده‌تر و آشوب‌ناک، مانند آونگ‌دوبل، بسیار مفید است.

## تقریب زاویه کوچک
بندانگشتی|300px|تقریب زاویه کوچک برای تابع سینوس: برای ، می‌توان تقریب  را در نظر گرفت.
معادله دیفرانسیل ارائه شده در بالا به سادگی حل نمی‌شود و هیچ راه حلی وجود ندارد که بتوان آن را بر حسب توابع مقدماتی نوشت. با این حال، با افزودن یک محدودیت به اندازه دامنه نوسان، فرمی به دست می‌آید که حل آن به سادگی امکان‌پذیر است. اگر فرض شود که زاویه بسیار کمتر از ۱ رادیان است (اغلب کمتر از ۰٫۱ رادیان یا حدود ۶ درجه ذکر می‌شود)، یا:
\theta \ll 1,
آنگاه با جایگزینی  در  با استفاده از تقریب زاویه کوچک:
\sin\theta\approx\theta,
معادله یک نوسانگر هماهنگ به دست می‌آید:
\frac{d^2\theta}{dt^2}+\frac{g}{\ell} \theta=0.

خطای ناشی از این تقریب از مرتبه  است (بر اساس بسط تیلور برای ).

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

این حرکت یک حرکت هماهنگ ساده است که در آن  دامنه نوسان (یعنی حداکثر زاویه بین میله آونگ و خط عمود) است. دوره تناوب تقریبی حرکت برابر است با:

که به عنوان قانون کریستیان هویگنس برای دوره تناوب شناخته می‌شود. توجه داشته باشید که تحت تقریب زاویه کوچک، دوره تناوب مستقل از دامنه  است؛ این خاصیت هم‌زمانی (ایزوکرونیسم) است که گالیله کشف کرد.

## قاعده سرانگشتی برای طول آونگ
T_0 = 2\pi\sqrt{\frac \ell g} نتیجه می‌دهد:  \ell = \frac{g}{\pi^2}\frac{T_0^2} 4.

اگر از واحدهای SI استفاده شود (یعنی اندازه‌گیری بر حسب متر و ثانیه)، و با فرض اینکه اندازه‌گیری در سطح زمین انجام می‌شود، آنگاه  و  (۰٫۹۹۴ تقریب تا ۳ رقم اعشار است).

بنابراین، تقریب‌های نسبتاً معقول برای طول و دوره تناوب عبارتند از:
\begin{align} \ell &\approx \frac{T_0^2}{4}, \\ T_0 &\approx 2 \sqrt\ell \end{align}
که در آن  تعداد ثانیه‌های بین *دو* تپش (یک تپش برای هر طرف نوسان) است و  بر حسب متر اندازه‌گیری می‌شود.

## دوره تناوب با دامنه دلخواه
بندانگشتی|'''شکل ۳.''' انحراف دوره تناوب «واقعی» آونگ از تقریب زاویه کوچک. مقدار «واقعی» با ارزیابی عددی انتگرال بیضوی به دست آمده است.|چپ|300px

بندانگشتی|'''شکل ۴.''' خطاهای نسبی با استفاده از سری توانی برای دوره تناوب.|چپ|300px

[بندانگشتی|312x312px|'''شکل ۵.''' انرژی پتانسیل و [پرتره فاز یک آونگ ساده. توجه داشته باشید که محور  که زاویه است، پس از هر ۲ رادیان به خود بازمی‌گردد.]]

برای دامنه‌های فراتر از تقریب زاویه کوچک، می‌توان دوره تناوب دقیق را ابتدا با معکوس کردن معادله سرعت زاویه‌ای به دست آمده از روش انرژی () محاسبه کرد:
\frac{dt}{d\theta} = \sqrt\frac\ell{2g}\frac{1}{\sqrt{\cos\theta-\cos\theta_0}}
و سپس انتگرال‌گیری روی یک چرخه کامل،
T = t(\theta_0 \rightarrow 0 \rightarrow -\theta_0 \rightarrow 0 \rightarrow\theta_0),
یا دو برابر نیم‌چرخه
T = 2 t(\theta_0 \rightarrow 0 \rightarrow -\theta_0),
یا چهار برابر ربع‌چرخه
T = 4 t(\theta_0 \rightarrow 0),
که منجر می‌شود به:
T = 4\sqrt\frac\ell{2g}\int^{\theta_0}_0 \frac{d\theta}{\sqrt{\cos\theta-\cos\theta_0}} .

توجه داشته باشید که این انتگرال با نزدیک شدن  به خط عمود واگرا می‌شود:
 \lim_{\theta_0 \to \pi} T = \infty, 
بنابراین آونگی با انرژی دقیقاً کافی برای عمودی شدن، هرگز واقعاً به آنجا نخواهد رسید. (برعکس، آونگی که نزدیک به حداکثر خود است می‌تواند زمان دلخواهی طول بکشد تا پایین بیاید.)

این انتگرال را می‌توان بر حسب انتگرال‌های بیضوی به صورت زیر بازنویسی کرد:
T = 4\sqrt\frac\ell g F\left( \frac{\pi} 2, \sin \frac{\theta_0} 2\right) 
که در آن  انتگرال بیضوی ناقص نوع اول است که به صورت زیر تعریف می‌شود:
F(\varphi , k) = \int_0^\varphi \frac {du} {\sqrt{1-k^2\sin^2 u}}\,.

یا به طور مختصرتر با جانشانی:
\sin{u} = \frac{\sin\frac{\theta}{2}}{\sin\frac{\theta_0}{2}}
که  را بر حسب  بیان می‌کند،

در اینجا  انتگرال بیضوی کامل نوع اول است که به صورت زیر تعریف می‌شود:

K(k) = F \left( \frac \pi 2, k \right) = \int_0^\frac{\pi}{2} \frac{du}{\sqrt{1-k^2\sin^2 u}}\,.
برای مقایسه تقریب با راه حل کامل، دوره تناوب یک آونگ به طول ۱ متر را روی زمین () در زاویه اولیه ۱۰ درجه در نظر بگیرید:
4\sqrt{\frac{1\text{ m}}{g}}\ K\left(\sin\frac{10^\circ} {2} \right)\approx 2.0102\text{ s}.
تقریب خطی (زاویه کوچک) مقدار زیر را می‌دهد:

2\pi \sqrt{\frac{1\text{ m}}{g}} \approx 2.0064\text{ s}.

تفاوت بین این دو مقدار، کمتر از ۰٫۲٪ است که بسیار کمتر از خطای ناشی از تغییرات  در موقعیت‌های جغرافیایی مختلف است.

از اینجا به بعد راه‌های زیادی برای محاسبه انتگرال بیضوی وجود دارد.

### راه حل چندجمله‌ای لژاندر برای انتگرال بیضوی
با توجه به (معادله ۳) و راه حل چندجمله‌ای لژاندر برای انتگرال بیضوی:
K(k) =\frac{\pi}{2}\sum_{n=0}^\infty \left(\frac{(2n-1)!!}{(2n)!!}k^{n}\right)^{2}
که در آن  نشان‌دهنده دوفاکتوریل است، یک راه حل دقیق برای دوره تناوب آونگ ساده عبارت است از:
\begin{alignat}{2}
T & = 2\pi \sqrt\frac \ell g \left( 1+ \left( \frac{1}{2} \right)^2 \sin^2 \frac{\theta_0}{2} + \left( \frac{1 \cdot 3}{2 \cdot 4} \right)^2 \sin^4 \frac{\theta_0}{2} + \left( \frac {1 \cdot 3 \cdot 5}{2 \cdot 4 \cdot 6} \right)^2 \sin^6 \frac{\theta_0}{2} + \cdots \right) \\
 & = 2\pi \sqrt\frac\ell g \cdot \sum_{n=0}^\infty \left( \left ( \frac{(2n)!}{( 2^n \cdot n! )^2} \right )^2 \cdot \sin^{2 n}\frac{\theta_0}{2} \right).\end{alignat}

شکل ۴ خطاهای نسبی را با استفاده از سری توانی نشان می‌دهد.  تقریب خطی است، و  تا  به ترتیب شامل جملاتی تا توان‌های ۲ تا ۱۰ هستند.

### راه حل سری توانی برای انتگرال بیضوی
فرمول‌بندی دیگری از راه حل فوق را می‌توان یافت اگر از سری مک‌لورن زیر:
\sin \frac{\theta_0} 2 =\frac12\theta_0 - \frac{1}{48}\theta_0^3 + \frac{1}{3\,840}\theta_0^5 - \frac{1}{645\,120}\theta_0^7 + \cdots.
در راه حل چندجمله‌ای لژاندر بالا استفاده شود.
سری توانی حاصل عبارت است از:

T = 2\pi \sqrt\frac\ell g \left( 1+ \frac{1}{16}\theta_0^2 + \frac{11}{3\,072}\theta_0^4 + \frac{173}{737\,280}\theta_0^6 + \frac{22\,931}{1\,321\,205\,760}\theta_0^8 + \frac{1\,319\,183}{951\,268\,147\,200}\theta_0^{10} + \frac{233\,526\,463}{2\,009\,078\,326\,886\,400}\theta_0^{12} + \cdots \right),
کسرهای بیشتر در دانشنامه برخط دنباله‌های صحیح (OEIS) موجود است؛ صورت کسرها در دنباله  و مخرج‌ها در دنباله  قرار دارند.

### راه حل میانگین حسابی-هندسی برای انتگرال بیضوی
با توجه به (معادله ۳) و راه حل میانگین حسابی-هندسی برای انتگرال بیضوی:
K(k) = \frac {\pi}{2 M(1-k,1+k)},
که در آن  میانگین حسابی-هندسی  و  است.

این منجر به یک فرمول جایگزین و با همگرایی سریع‌تر برای دوره تناوب می‌شود:
T = \frac{2\pi}{M\left(1, \cos\frac{\theta_0} 2 \right)} \sqrt\frac\ell g.

اولین تکرار (ایترِیشن) این الگوریتم نتیجه می‌دهد:
T_1 = \frac{2T_0}{1 + \cos \frac{\theta_0}{2}}. 

این تقریب دارای خطای نسبی کمتر از ۱٪ برای زوایای تا ۹۶٫۱۱ درجه است. از آنجا که \frac{1}{2}\left(1+\cos\left(\frac{\theta_0}{2}\right)\right) = \cos^2 \frac{\theta_0}{4},، این عبارت را می‌توان به صورت خلاصه‌تر نوشت:
T_1 = T_0 \sec^2 \frac{\theta_0}{4}.

بسط مرتبه دوم \sec^2(\theta_0/4) به T \approx T_0 \left(1 + \frac{\theta_0^2}{16} \right) کاهش می‌یابد.

تکرار دوم این الگوریتم نتیجه می‌دهد:
T_2 = \frac{4T_0}{1 + \cos \frac{\theta_0}{2} + 2\sqrt{\cos \frac{\theta_0}{2}}} = \frac{4T_0}{\left(1 + \sqrt{\cos \frac{\theta_0}{2}} \right)^2}. 

این تقریب دوم دارای خطای نسبی کمتر از ۱٪ برای زوایای تا ۱۶۳٫۱۰ درجه است.

## فرمول‌های تقریبی برای دوره تناوب آونگ غیرخطی
اگرچه دوره تناوب دقیق T را می‌توان برای هر دامنه محدود  با محاسبه انتگرال بیضوی کامل مربوطه K(k) (که در آن k \equiv \sin(\theta_0/2)) تعیین کرد، اما در کاربردها اغلب از این کار اجتناب می‌شود زیرا بیان این انتگرال به صورت یک فرم بسته بر حسب توابع مقدماتی امکان‌پذیر نیست. این امر راه را برای تحقیق بر روی فرمول‌های تقریبی ساده برای افزایش دوره تناوب آونگ با افزایش دامنه باز کرده است (که در آزمایشگاه‌های فیزیک مقدماتی، مکانیک کلاسیک، الکترومغناطیس، آکوستیک، الکترونیک، ابررسانایی و غیره مفید است). فرمول‌های تقریبی یافت شده توسط نویسندگان مختلف را می‌توان به صورت زیر طبقه‌بندی کرد:

* فرمول‌های «زاویه نه چندان بزرگ»، یعنی آنهایی که تخمین‌های خوبی برای دامنه‌های زیر \pi/2 رادیان (یک حد طبیعی برای یک وزنه در انتهای یک ریسمان انعطاف‌پذیر) ارائه می‌دهند، اگرچه انحراف نسبت به دوره دقیق به صورت یکنواخت با دامنه افزایش می‌یابد و برای دامنه‌های نزدیک به \pi رادیان نامناسب هستند. یکی از ساده‌ترین فرمول‌های یافت شده در متون، فرمول زیر توسط لیما (۲۰۰۶) است: T \approx -\,T_0 \, \frac{\ln{a}}{1-a}، که در آن a \equiv \cos{(\theta_0/2)}.
* فرمول‌های «زاویه بسیار بزرگ»، یعنی آنهایی که دوره دقیق را به صورت مجانبی برای دامنه‌های نزدیک به \pi رادیان تقریب می‌زنند، با خطایی که برای دامنه‌های کوچک‌تر به صورت یکنواخت افزایش می‌یابد (یعنی برای دامنه‌های کوچک نامناسب هستند). یکی از بهترین این فرمول‌ها، فرمول کرومر (Cromer) است: T \approx \frac{2}{\pi}\,T_0\,\ln{(4/a)}.

البته، افزایش T با دامنه زمانی که \pi/2 باشد آشکارتر است، همانطور که در آزمایش‌های بسیاری با استفاده از میله صلب یا دیسک مشاهده شده است. از آنجا که زمان‌سنج‌ها و حسگرهای دقیق در حال حاضر حتی در آزمایشگاه‌های فیزیک مقدماتی موجود هستند، خطاهای تجربی یافت شده در آزمایش‌های «زاویه بسیار بزرگ» به اندازه کافی کوچک هستند که با دوره دقیق مقایسه شوند، و توافق بسیار خوبی بین تئوری و آزمایش‌هایی که در آن اصطکاک ناچیز است، یافت شده است. از آنجا که این فعالیت توسط بسیاری از مربیان تشویق شده است، جستجو برای یک فرمول تقریبی ساده برای دوره تناوب آونگ که برای تمام دامنه‌های ممکن معتبر باشد و بتوان داده‌های تجربی را با آن مقایسه کرد، انجام شد. در سال ۲۰۰۸، لیما یک فرمول میانگین وزنی با این ویژگی استخراج کرد:
T \approx \frac{r\,a^2\,T_\text{Lima} +k^2\,T_\text{Cromer}}{r\,a^2+k^2},
که در آن r = 7.17 است، که حداکثر خطای تنها ۰٫۶٪ (در \theta_0 = 95^\circ) را نشان می‌دهد.

## جابجایی زاویه‌ای با دامنه دلخواه
بسط سری فوریه برای \theta(t) به صورت زیر داده می‌شود:


\theta(t) = 8\sum_{n \geq 1\text{ فرد}}\frac{(-1)^{\left\lfloor{n/2}\right\rfloor}}{n}\frac{q^{n/2}}{1+q^{n}}\cos(n\omega t)


که در آن q نوم بیضوی (elliptic nome) است، q=\exp\left({-\pi K\bigl(\sqrt{\textstyle 1-k^2}\bigr) \big/ K(k)}\right), k=\sin (\theta_0/2), و \omega=2\pi/T بسامد زاویه‌ای است.

اگر تعریف کنیم:
\varepsilon=\frac12\cdot \frac{1-\sqrt{\cos(\theta_0/2)}}{1+\sqrt{\cos(\theta_0/2)}} 
q را می‌توان با استفاده از بسط زیر تقریب زد:
q = \varepsilon + 2\varepsilon^5 + 15\varepsilon^{9} + 150\varepsilon^{13} + 1707\varepsilon^{17} + 20910\varepsilon^{21} + \cdots 
(ببینید ). توجه داشته باشید که \varepsilon  برای \theta_0 است، بنابراین این تقریب حتی برای دامنه‌های بزرگ نیز قابل استفاده است.
به طور معادل، زاویه را می‌توان بر حسب تابع بیضوی ژاکوبی \operatorname{cd} با مدول k بیان کرد:
\theta(t)=2\arcsin\left(k\operatorname{cd}\left(\sqrt{\frac{g}{\ell}}t;k\right)\right),\quad k=\sin\frac{\theta_0}{2}.

برای x کوچک، داریم \sin x\approx x، \arcsin x\approx x و \operatorname{cd}(t;0)=\cos t، بنابراین راه حل به خوبی توسط راه حل ارائه شده در بخش «تقریب زاویه کوچک» تخمین زده می‌شود.

## مثال‌ها
انیمیشن‌های زیر حرکت یک آونگ ساده (بدون اصطکاک) را با مقادیر فزاینده جابجایی اولیه وزنه، یا به طور معادل سرعت اولیه فزاینده، نشان می‌دهند. نمودار کوچک بالای هر آونگ، نمودار صفحه فاز مربوطه است؛ محور افقی جابجایی و محور عمودی سرعت است. با سرعت اولیه به اندازه کافی زیاد، آونگ دیگر به عقب و جلو نوسان نمی‌کند بلکه به طور کامل دور نقطه اتکا می‌چرخد.


File:Pendulum_0deg.gif|زاویه اولیه ۰ درجه، تعادل پایدار
File:Pendulum_45deg.gif|زاویه اولیه ۴۵ درجه
File:Pendulum_90deg.gif|زاویه اولیه ۹۰ درجه
File:Pendulum_135deg.gif|زاویه اولیه ۱۳۵ درجه
File:Pendulum_170deg.gif|زاویه اولیه ۱۷۰ درجه
File:Pendulum_180deg.gif|زاویه اولیه ۱۸۰ درجه، تعادل ناپایدار
File:Pendulum_190deg.gif|آونگ با انرژی که به سختی برای یک چرخش کامل کافی است
File:Pendulum_220deg.gif|آونگ با انرژی کافی برای یک چرخش کامل


## آونگ مرکب
یک **آونگ مرکب** (یا **آونگ فیزیکی**) آونگی است که در آن میله بدون جرم نیست و ممکن است دارای ابعاد گسترده باشد؛ یعنی یک جسم صلب با شکل دلخواه که حول یک نقطه اتکا O نوسان می‌کند. در این حالت دوره تناوب آونگ به گشتاور لختی آن I_O حول نقطه اتکا بستگی دارد.

معادله گشتاور بیان می‌کند:
\tau = I \alpha
که در آن:
: \alpha شتاب زاویه‌ای است.
: \tau گشتاور نیرو است.

گشتاور توسط نیروی گرانش ایجاد می‌شود، بنابراین:
\tau = - m g r_\oplus \sin\theta
که در آن:
* m جرم کل جسم صلب (میله و وزنه) است.
* r_\oplus فاصله نقطه اتکا تا مرکز جرم سیستم است.
* \theta زاویه نسبت به خط عمود است.

بنابراین، تحت تقریب زاویه کوچک، \sin\theta\approx\theta (یا به طور معادل زمانی که \theta_\mathrm{max}\ll 1)،
\alpha = \ddot{\theta} = \frac{mgr_\oplus}{I_O}\sin\theta \approx -\frac{mgr_\oplus}{I_O}\theta
که در آن I_O گشتاور لختی جسم حول نقطه اتکا O است.

عبارت برای \alpha همان فرم آونگ ساده معمولی را دارد و دوره تناوبی برابر با مقدار زیر می‌دهد:
T = 2 \pi \sqrt{\frac{I_O}{mgr_\oplus}}

و بسامدی برابر با:
f = \frac{1}{T} = \frac{1}{2\pi} \sqrt{\frac{mgr_\oplus}{I_O}}

اگر زاویه اولیه در نظر گرفته شود (برای دامنه‌های بزرگ)، آنگاه عبارت برای \alpha تبدیل می‌شود به:
\alpha = \ddot{\theta} = -\frac{mgr_\oplus}{I_O}\sin\theta
و دوره تناوبی برابر با مقدار زیر می‌دهد:
T = 4 \operatorname{K}\left(\sin^2\frac{\theta_\mathrm{max}}{2}\right) \sqrt{\frac{I_O}{mgr_\oplus}}
که در آن \theta_\mathrm{max} حداکثر زاویه نوسان (نسبت به خط عمود) و \operatorname{K}(k) انتگرال بیضوی کامل نوع اول است.

یک مفهوم مهم **طول معادل**، \ell^\mathrm{eq} است، که طول یک آونگ ساده است که همان بسامد زاویه‌ای \omega_0 آونگ مرکب را دارد:
 {\omega_0}^2 = \frac{g}{\ell^\mathrm{eq}} := \frac{mgr_\oplus}{I_O} \implies \ell^\mathrm{eq} = \frac{I_O}{mr_\oplus} 

موارد زیر را در نظر بگیرید:
* آونگ ساده حالت خاصی است که در آن تمام جرم در وزنه متمرکز شده و در فاصله \ell از نقطه اتکا نوسان می‌کند. بنابراین، r_\oplus=\ell و I_O=m\ell^2 است، پس عبارت کاهش می‌یابد به: {\omega_0}^2 = \frac{mgr_\oplus}{I_O}=\frac{mg\ell}{m\ell^2}=\frac{g}{\ell}. توجه کنید که \ell^\mathrm{eq}=\ell، همانطور که انتظار می‌رفت (تعریف طول معادل).
* یک میله همگن با جرم m و طول \ell که از انتهای خود آویزان است دارای r_\oplus=\frac{1}{2}\ell و I_O=\frac{1}{3}m\ell^2 است، بنابراین عبارت کاهش می‌یابد به: {\omega_0}^2 = \frac{mgr_\oplus}{I_O}=\frac{mg\,\frac{1}{2}\ell}{\frac{1}{3}m\ell^2}=\frac{g}{\frac{2}{3}\ell}. توجه کنید که \ell^\mathrm{eq}=\frac{2}{3}\ell؛ یک میله همگن طوری نوسان می‌کند که گویی یک آونگ ساده با طولی برابر با دو سوم طول خودش است.
* یک آونگ ساده سنگین: ترکیبی از یک میله همگن با جرم m_\mathrm{rod} و طول \ell که از انتهای خود آویزان است و یک وزنه m_\mathrm{bob} در انتهای دیگر. در این حالت سیستم دارای جرم کل m_\mathrm{bob}+m_\mathrm{rod} است و سایر پارامترها عبارتند از m r_\oplus=m_\mathrm{bob}\ell+m_\mathrm{rod}\frac{\ell}{2} (طبق تعریف مرکز جرم) و I_O=m_\mathrm{bob}\ell^2+\frac{1}{3}m_\mathrm{rod}\ell^2، بنابراین عبارت کاهش می‌یابد به:
{\omega_0}^2 = \frac{mgr_\oplus}{I_O}=\frac{\left(m_\mathrm{bob}\ell+m_\mathrm{rod}\frac{\ell}{2}\right)g}{m_\mathrm{bob}\ell^2+\frac{1}{3}m_\mathrm{rod}\ell^2} = \frac{g}{\ell} \frac{m_\mathrm{bob}+\frac{m_\mathrm{rod}}{2}}{m_\mathrm{bob}+\frac{m_\mathrm{rod}}{3}} = \frac{g}{\ell} \frac{1+\frac{m_\mathrm{rod}}{2m_\mathrm{bob}}}{1+\frac{m_\mathrm{rod}}{3m_\mathrm{bob}}} 
که در آن \ell^\mathrm{eq} = \ell \frac{1+\frac{m_\mathrm{rod}}{3m_\mathrm{bob}}} {1+\frac{m_\mathrm{rod}}{2m_\mathrm{bob}}} . توجه کنید که این فرمول‌ها را می‌توان تنها با در نظر گرفتن جرم میله یا وزنه به عنوان صفر، به دو مورد قبلی تعمیم داد. همچنین توجه کنید که فرمول به جرم وزنه و میله به طور جداگانه بستگی ندارد، بلکه به نسبت آن‌ها، \frac{m_\mathrm{rod}}{m_\mathrm{bob}} بستگی دارد. یک تقریب می‌تواند برای \frac{m_\mathrm{rod}}{m_\mathrm{bob}}\ll 1 ساخته شود:

{\omega_0}^2 \approx \frac{g}{\ell} \left( 1+\frac{1}{6}\frac{m_\mathrm{rod}}{m_\mathrm{bob}}+\cdots\right)

توجه کنید که این چقدر شبیه به بسامد زاویه‌ای در یک سیستم فنر-جرم با جرم مؤثر است.

## آونگ میرا و واداشته
بحث بالا فقط بر روی وزنه آونگی تمرکز داشت که تنها تحت تأثیر نیروی گرانش بود. فرض کنید یک نیروی میرایی (استهلاک)، مثلاً مقاومت هوا، و همچنین یک نیروی محرک (رانش) سینوسی بر جسم وارد شود. این سیستم یک نوسانگر هماهنگ میرا و واداشته (Driven) است و رفتاری آشوب‌ناک دارد.

معادله (۱) می‌تواند به صورت زیر نوشته شود:
m\ell^2 \frac{d^2 \theta}{dt^2}  = -mg\ell \sin \theta

(به اثبات گشتاور برای معادله (۱) در بالا مراجعه کنید).

یک جمله میرایی (استهلاک) و یک جمله نیروی محرک (رانش) می‌توانند به سمت راست اضافه شوند تا رابطه زیر به دست آید:

m\ell^2 \frac{d^2 \theta}{dt^2} = -mg\ell\sin \theta - b\frac{d\theta}{dt} + a\cos(\Omega t)

که در آن فرض شده است میرایی مستقیماً متناسب با سرعت زاویه‌ای است (این فرض برای مقاومت هوا با سرعت پایین صحیح است، همچنین ببینید نیروی پسار). a و b ثابت‌هایی هستند که به ترتیب دامنه نیروی محرک و درجه میرایی را تعیین می‌کنند. \Omega بسامد زاویه‌ای نوسانات محرک است.

با تقسیم طرفین بر m\ell^2 داریم:

\frac{d^2 \theta}{dt^2} + \frac{b}{m\ell^2}\frac{d\theta}{dt} + \frac{g}{\ell}{\sin \theta} - \frac{a}{m\ell^2}\cos (\Omega t) = 0.

برای یک آونگ فیزیکی (مرکب):

\frac{d^2 \theta}{dt^2} + \frac{b}{I}\frac{d\theta}{dt} + \frac{mgr_\oplus}{I}{\sin \theta} - \frac{a}{I}\cos (\Omega t) = 0.

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

## تعبیر فیزیکی دوره تناوب موهومی
تابع بیضوی ژاکوبی که موقعیت یک آونگ را به عنوان تابعی از زمان بیان می‌کند، یک تابع دوچندان متناوب با یک دوره تناوب حقیقی و یک دوره تناوب موهومی است. دوره تناوب حقیقی، طبیعتاً زمانی است که طول می‌کشد تا آونگ یک چرخه کامل را طی کند. پل اپل یک تعبیر فیزیکی برای دوره تناوب موهومی ارائه کرد: اگر  حداکثر زاویه یک آونگ باشد و  حداکثر زاویه آونگ دیگر، آنگاه دوره تناوب حقیقی هر یک، برابر با بزرگی دوره تناوب موهومی دیگری است.

## آونگ‌های جفت‌سده

**آونگ‌های جفت‌شده** می‌توانند بر حرکت یکدیگر تأثیر بگذارند، چه از طریق اتصال مستقیم (مانند فنری که وزنه‌ها را به هم متصل می‌کند) و چه از طریق حرکات در ساختار تکیه‌گاه (مانند سطح میز). معادلات حرکت برای دو آونگ ساده یکسان که توسط یک فنر متصل به وزنه‌ها جفت شده‌اند، با استفاده از مکانیک لاگرانژی قابل استخراج است.

انرژی جنبشی سیستم عبارت است از:
E_\text{K}=\frac{1}{2}mL^2\left(\dot\theta_1^2+\dot\theta_2^2\right)
که در آن m جرم وزنه‌ها، L طول ریسمان‌ها و \theta_1، \theta_2 جابجایی‌های زاویه‌ای دو وزنه از وضعیت تعادل هستند.

انرژی پتانسیل سیستم عبارت است از:
E_\text{p}=mgL(2-\cos\theta_1-\cos\theta_2)+\frac{1}{2}kL^2(\theta_2-\theta_1)^2

که در آن g شتاب گرانشی و k ثابت فنر است. جابجایی L(\theta_2-\theta_1) فنر از موقعیت تعادل خود، با فرض تقریب زاویه کوچک در نظر گرفته شده است.

لاگرانژی سیستم به صورت زیر خواهد بود:
\mathcal{L}=\frac{1}{2}mL^2\left(\dot\theta_1^2+\dot\theta_2^2\right)-mgL(2-\cos\theta_1-\cos\theta_2)-\frac{1}{2} k L^2(\theta_2-\theta_1)^2
که منجر به مجموعه معادلات دیفرانسیل جفت‌شده زیر می‌شود:
\begin{align}
\ddot\theta_1+\frac{g}{L}\sin\theta_1+\frac{k}{m}(\theta_1-\theta_2)&=0 \\
\ddot\theta_2+\frac{g}{L}\sin\theta_2-\frac{k}{m}(\theta_1-\theta_2)&=0
\end{align}

با جمع و تفریق این دو معادله و اعمال تقریب زاویه کوچک، دو معادله نوسانگر هماهنگ بر حسب متغیرهای \theta_1+\theta_2 و \theta_1-\theta_2 به دست می‌آید:
\begin{align}
\ddot\theta_1+\ddot\theta_2+\frac{g}{L}(\theta_1+\theta_2)&=0 \\
\ddot\theta_1-\ddot\theta_2+\left(\frac{g}{L}+2\frac{k}{m}\right)(\theta_1-\theta_2)&=0
\end{align}
که پاسخ‌های متناظر آن‌ها عبارتند از:
\begin{align}
\theta_1+\theta_2&=A\cos(\omega_1t+\alpha) \\
\theta_1-\theta_2&=B\cos(\omega_2t+\beta)
\end{align}
که در آن:
\begin{align}
\omega_1&=\sqrt{\frac{g}{L}} \\
\omega_2&=\sqrt{\frac{g}{L}+2\frac{k}{m}}
\end{align}

و A، B، \alpha، \beta ثابت‌های انتگرال‌گیری هستند.

با بازنویسی راه حل‌ها تنها بر حسب \theta_1 و \theta_2 داریم:
\begin{align}
\theta_1&=\frac{1}{2}A\cos(\omega_1t+\alpha)+\frac{1}{2}B\cos(\omega_2t+\beta) \\
\theta_2&=\frac{1}{2}A\cos(\omega_1t+\alpha)-\frac{1}{2}B\cos(\omega_2t+\beta)
\end{align}

اگر به وزنه‌ها ضربه (فشار) اولیه وارد نشود، شرط \dot\theta_1(0)=\dot\theta_2(0)=0 ایجاب می‌کند که \alpha=\beta=0 باشد، که (پس از کمی جابجایی) نتیجه می‌دهد:
\begin{align}
A&=\theta_1(0)+\theta_2(0)\\
B&=\theta_1(0)-\theta_2(0)
\end{align}

## جستارهای وابسته
* آونگ دوتایی
* آونگ کاپیتسا

## پانویس

## برای مطالعه بیشتر
*
*
*

## پیوند به بیرون
*مقاله مث‌ورلد درباره تابع ماتیوس

رده:معادلات دیفرانسیل
رده:آونگ‌ها
رده:سامانه‌های دینامیکی
رده:ریاضی فیزیک