Хармониски треперник

Хармониски треперник (или хармониски осцилатор) — систем во класичната механика, на којшто, кога ќе му биде нарушена рамнотежната состојба, сè создава повратна сила, F, пропорционална на поместувањето, x:

{\vec {F}}=-k{\vec {x}}\,

каде k е позитивна непроменлива.

Ако F единствената сила која делува на системот, сситемот се нарекува едноставен хармониски треперник, и е подложен на едноставно хармониско движење: синусоидно треперење околу рамнотежна точка, со постојан замав (амплитуда) и постојана честота (која не зависи од замавот).

Ако имаме сила на триење (придушување) пропорционална на брзината, хармонискиот треперник се опишува како придушен треперник. Во зависност од коефициентот на триење, системот може:

Трепери со фрквенција помала отколку таа во непридушениот случај, и замав кој се намалува со текот на времето (подпридушен треперник).
Се враќа во рамнотежната положба, без треперење (презадушен треперник).

Граничното решение меѓу подпридушен треперник и презадушен треперник се случува при одредена вредност на коефициентот на триење, која се нарекува „критично придушување“.

Ако е присутна временско зависна сила, хармонискиот треперник се опишува како присилен треперник.

Механичките примери вклучуваат нишало (со малоаголни поместувања), маси поврзани на пружини, и акустични системи. Други подеднакви системи вклучуваат електрични хармониски треперења како што се RLC кола. Моделот на хармонискиот треперник е од важност за физиката, бидејќи секоја маса која е под дејство на силаво стабилна рамнотежа се однесува како хармониски треперник за мали вибрации. Хармониските треперници ги има насекаде во природата и се искористуваат во многу уреди направени од човекот, како што се часовници и радио кола. Тие се создавачи на сите синусоидни вибрации и бранови.

Едноставен хармониски осцилатор

Едноставен хеамониски треперник е треперник кој не е ни присилен но не е придушен. Истиот се состои од маса m, која дејствува со единствена сила, F, која ја привлекува масата во насока на точката x=0 и зависи само од местоположбата на масата во x и постојаната k. Рамнотежната состојба на силите (Втор Њутнов закон) за системот е:

F=ma=m{\frac {\mathrm {d} ^{2}x}{\mathrm {d} t^{2}}}=-kx.

Решавајќи ја оваа диференцијална равенка, се добива дека движењето е опишано со функцијата:

x(t)=A\cos \left(\omega t+\phi \right),

каде

\omega ={\sqrt {\frac {k}{m}}}={\frac {2\pi }{T}}.

Движењето е периодично, повторувајќи се на синусоиден начин со постојан замав, A. Покрај замавот, движењето на едноставниот хармониски треперник е опишан со сопствениот период T, времето потребно за едно треперење или неговата честота f = ¹⁄_T, бројот на треперења во единица време. Местоположбата во одредено време t исто тка зависи од фазата, φ, којашто ја одредува почетната точка на синусниот бран. Периодот и честотата се одредени од големината на масата m и постојаната на силата k, додека замавот и фазата се одредени од почетната местоположба и брзина.

Брзината и забрзувањето на едноставниот хармонискиот треперник, трепери со истата честота и местоположба но со променети фази. Брзината е максимална за нулта поместување, додека пак забрзувањето е со спротивна насока од поместувањето.

Потенцијалната енергија складирана во едноставниот хармониски треперник во местоположбата x е:

U={\frac {1}{2}}kx^{2}.

Придушен хармониски треперник

Однесувањето на системот во зависност од вредноста на придушувачкиот однос ζ

Придушен хармониски треперник, којшто се придушува поради силата на триење

Кај вистинските треперници, триењето, или придушувањето, го успорува движењето на системот. Поради силата на триење, брзината се намалува во однос силата на триење која делува. Додека едноставните хармониски движења треперат со обновливата сила која делува на системот, придушените хармониски движења се под дејство на сила на триење. Во многу системи кои вибрираат силата на триење F_f може да биде пропорционален со брзината v на предметот: F_f = −cv, каде c се нарекува вискозен придушен коефициент.

Рамнотежната на сила (Втор Њутнов закон) за придушени треперења е:

F=-kx-c{\frac {\mathrm {d} x}{\mathrm {d} t}}=m{\frac {\mathrm {d} ^{2}x}{\mathrm {d} t^{2}}}.

TОва се презапишува во обликот

{\frac {\mathrm {d} ^{2}x}{\mathrm {d} t^{2}}}+2\zeta \omega _{0}{\frac {\mathrm {d} x}{\mathrm {d} t}}+\omega _{0}^{\,2}x=0,

каде

\omega _{0}={\sqrt {\frac {k}{m}}}

се нарекува 'непридушена аголна честота на треперникот' и

\zeta ={\frac {c}{2{\sqrt {mk}}}}

се нарекува 'однос на придушување'.

Чекор на придушување на хармонискиот треперник, кривите се нацртани за три вредности на μ = ω₁ = ω₀√1−ζ². Времето е изразено во единици на времето на распаѓање τ = 1/(ζω₀).

Вредноста на односот на придушување ζ критично го определува однесувањето на системот. Тогаш придушениот хармониски треперник може да биде:

Презадушен (ζ > 1): Системот се враќа (експоненцијалното распаѓање) во почетна положба без треперење. Поголемите вредности на односот на придушување ζ се враќаат поспоро во рамнотежната состојба.
Критично придушување (ζ = 1): Системот се враќа во почетната положба што е можно побрзо без притоа да дојде до треперење. Ова е посакувано дејство кај придушните системи како што се.
Непридушени (ζ < 1): Системот трепери (со малку поразлична честота отколку во непридушениот случај) при што замавот постепено се намалува до нула. Аголната честота на непридушениот хармониски треперник е определена со:

\omega _{1}=\omega _{0}{\sqrt {1-\zeta ^{2}}}.

Q факторот на придушениот треперник сè запишува как:

Q=2\pi \times {\frac {\text{Skladirana energija}}{\text{Izgubena energija vo eden ciklus}}}.

Q е поврзана со односот на придушување преку равенството:

$Q={\frac {1}{2\zeta }}.$

Присилен хармониски треперник

Присилените хармониски треперници сè придушени треперници кои се под дејство на надворешно применета сила F(t).

Втор Њутнов закон ја добива формата

F(t)-kx-c{\frac {\mathrm {d} x}{\mathrm {d} t}}=m{\frac {\mathrm {d} ^{2}x}{\mathrm {d} t^{2}}}.

Кој обично се презапишува во формата

{\frac {\mathrm {d} ^{2}x}{\mathrm {d} t^{2}}}+2\zeta \omega _{0}{\frac {\mathrm {d} x}{\mathrm {d} t}}+\omega _{0}^{2}x={\frac {F(t)}{m}}.

Оваа равенка може да биде решена точно за секоја присилена сила, користејќи ги решенијата z(t) кои ја задоволуваат равенката во која нема присуство на сили:

{\frac {\mathrm {d} ^{2}z}{\mathrm {d} t^{2}}}+2\zeta \omega _{0}{\frac {\mathrm {d} z}{\mathrm {d} t}}+\omega _{0}^{2}z=0,

и кои можат да бидат изразени како придушени сиснусоидални осцилации,

z(t)=A\mathrm {e} ^{-\zeta \omega _{0}t}\ \sin \left({\sqrt {1-\zeta ^{2}}}\ \omega _{0}t+\phi \right),

во истиот случај каде ζ ≤ 1. ЗамавотA и фазата φ го одредуваат однесувањето потребно за да се исполнат почетните услови.

Дискретна функција

Во овој случај ζ < 1 и единечна Хевисајдова функција со услов x(0) = 0:

{F(t) \over m}={\begin{cases}\omega _{0}^{2}&t\geq 0\\0&t<0\end{cases}}

решението е:

x(t)=1-\mathrm {e} ^{-\zeta \omega _{0}t}{\frac {\sin \left({\sqrt {1-\zeta ^{2}}}\ \omega _{0}t+\varphi \right)}{\sin(\varphi )}},

со фаза φ добиена од

\cos \varphi =\zeta .\,

Времето потребно на еден треперник за да се приспособи на променетите надворешни услови е со големина одредена како τ = 1/(ζω₀). Во физиката, приспособувањето е наречено релаксаторско, и τ се нарекува време на релаксација.

Синусоидна присилна сила

Постојана состојба на промена на замавот со честота и задушување на присилен едноставен хармониски треперник.^[1]^[2]

Во случај на синусоидна присилна сила:

{\frac {\mathrm {d} ^{2}x}{\mathrm {d} t^{2}}}+2\zeta \omega _{0}{\frac {\mathrm {d} x}{\mathrm {d} t}}+\omega _{0}^{2}x={\frac {1}{m}}F_{0}\sin(\omega t),

каде $\,\!F_{0}$ е замавот на присилната сила и $\,\!\omega$ е присилената честота на синусоиден присилувачки механизам. Овој вид на систем се појавува во AC присилените RLC кола (отпорник-намотка-кондензатор) и присилени пружински системи кои имаат внатрешна механичка отпорност или пак надворешен воздушен отпор.

Општото решение е сумата на премините кое зависи од почетните услови, и од стабилната состојба, која е независна од почетните услови и зависи само од амплитудата на прислината сила $\,\!F_{0}$ , присилената честота, $\,\!\omega$ , непридушената аголна честота $\,\!\omega _{0}$ , и односот на придушување $\,\!\zeta$ .

решението на стабилна состојба е пропорционално со присилната сила со настанатата фазна промена на $\,\!\phi$ :

x(t)={\frac {F_{0}}{mZ_{m}\omega }}\sin(\omega t+\phi )

каде

Z_{m}={\sqrt {\left(2\omega _{0}\zeta \right)^{2}+{\frac {1}{\omega ^{2}}}\left(\omega _{0}^{2}-\omega ^{2}\right)^{2}}}

е апсолутната вредност на импеданца и

\phi =\arctan \left({\frac {\omega _{0}^{2}-\omega ^{2}}{2\omega \omega _{0}\zeta }}\right)

е фазата на треперењето во однос на присилната сила, ако вредноста на аркустангенсот се земе да биде меѓу -180 степени и 0.

За одредена присилена честота наречена резонанса, или резонантна честота $\,\!\omega _{r}=\omega _{0}{\sqrt {1-2\zeta ^{2}}}$ , замавот (за дадена $\,\!F_{0}$ ) е максимален. Ефектот на резонанса се случува само кога $\,\zeta <1/{\sqrt {2}}$ , осносно за значајно непридушени системи. За силно непридушени системи вредноста на замавот ќе биде доволно голем близу до резонантната честота.

Параметарски треперник

Параметарски треперник е присилен хармониски треперник при кој присилната енергија се обезбедува од различните параметри на треперникот, како што се придушната и силата на обновување. Познат пример за параметарско треперење е лулашката која ја има на детските игралишта.^[3]^[4]^[5] Лицето кое се лула на лулашка може да го засили замавот на треперењата без надворешна присилна сила (туркање), со промената на моментот на инерција на лулашката лулајќи се нанапред и наназад или седнувајќи и станувајќи, во ритам со треперењата на системот. Примери за параметарски треперници можат да имаат променливи резонантни честоти $\omega$ и придушувањето $\beta$ .

Универзална равенка на треперник

Равенката:

{\frac {\mathrm {d} ^{2}q}{\mathrm {d} \tau ^{2}}}+2\zeta {\frac {\mathrm {d} q}{\mathrm {d} \tau }}+q=0

позната како универзална равенка на треперник бидејќи сите линиски треперници од втор ред можат да се сведат во таа. Ова е направено преку бездимензионизација.

Ако прсисилната функција е f(t) = cos(ωt) = cos(ωt_cτ) = cos(ωτ), where ω = ωt_c, равенката ја добива формата

{\frac {\mathrm {d} ^{2}q}{\mathrm {d} \tau ^{2}}}+2\zeta {\frac {\mathrm {d} q}{\mathrm {d} \tau }}+q=\cos(\omega \tau ).

Решението на диференцијалната равенка содржи два дела, „премин“ и „стабилна состојба“.

Решение на премин

Решението засновано на обична диференцијална равенка е за is for преодни постојани c₁ и c₂

$q_{t}(\tau )={\begin{cases}\mathrm {e} ^{-\zeta \tau }\left(c_{1}\mathrm {e} ^{\tau {\sqrt {\zeta ^{2}-1}}}+c_{2}\mathrm {e} ^{-\tau {\sqrt {\zeta ^{2}-1}}}\right)&\zeta >1{\text{ (overdamping)}}\\\mathrm {e} ^{-\zeta \tau }(c_{1}+c_{2}\tau )=\mathrm {e} ^{-\tau }(c_{1}+c_{2}\tau )&\zeta =1{\text{ (critical damping)}}\\\mathrm {e} ^{-\zeta \tau }\left[c_{1}\cos \left({\sqrt {1-\zeta ^{2}}}\tau \right)+c_{2}\sin \left({\sqrt {1-\zeta ^{2}}}\tau \right)\right]&\zeta <1{\text{(underdamping)}}\end{cases}}$

Решението на премин е независно од присилната функција.

Решение на стабилна состојба

Со примена на „ методот на комплексни променливи“ со решавање на помошната равенка прикажана подолу најде реалниот дел од решението:

{\frac {\mathrm {d} ^{2}q}{\mathrm {d} \tau ^{2}}}+2\zeta {\frac {\mathrm {d} q}{\mathrm {d} \tau }}+q=\cos(\omega \tau )+\mathrm {i} \sin(\omega \tau )=\mathrm {e} ^{\mathrm {i} \omega \tau }.

Претпоставеното решение е со облик:

\,\!q_{s}(\tau )=A\mathrm {e} ^{\mathrm {i} (\omega \tau +\phi )}.

Изводите од нулти до втори ред сè:

q_{s}=A\mathrm {e} ^{\mathrm {i} (\omega \tau +\phi )},\ {\frac {\mathrm {d} q_{s}}{\mathrm {d} \tau }}=\mathrm {i} \omega A\mathrm {e} ^{\mathrm {i} (\omega \tau +\phi )},\ {\frac {\mathrm {d} ^{2}q_{s}}{\mathrm {d} \tau ^{2}}}=-\omega ^{2}A\mathrm {e} ^{\mathrm {i} (\omega \tau +\phi )}.

Заменувајќи ги овие записи во диференцијалната равенка се добива

\,\!-\omega ^{2}A\mathrm {e} ^{\mathrm {i} (\omega \tau +\phi )}+2\zeta \mathrm {i} \omega A\mathrm {e} ^{\mathrm {i} (\omega \tau +\phi )}+A\mathrm {e} ^{\mathrm {i} (\omega \tau +\phi )}=(-\omega ^{2}A\,+\,2\zeta \mathrm {i} \omega A\,+\,A)\mathrm {e} ^{\mathrm {i} (\omega \tau +\phi )}=\mathrm {e} ^{\mathrm {i} \omega \tau }.

Делејќи со експоненционалниот поим од лево се добива:

\,\!-\omega ^{2}A+2\zeta \mathrm {i} \omega A+A=\mathrm {e} ^{-\mathrm {i} \phi }=\cos \phi -\mathrm {i} \sin \phi .

Пресметувањето на реалниот и комплексниот дел се сведува на две независни равенки:

A(1-\omega ^{2})=\cos \phi \qquad 2\zeta \omega A=-\sin \phi .

Замавен дел

Цртеж на честотниот дел на идеален хармониски треперник.

Квадрирајќи ги двете страни и собирајќи ги равенките се добива:

\left.{\begin{array}{rcl}A^{2}(1-\omega ^{2})^{2}&=&\cos ^{2}\phi \\[6pt](2\zeta \omega A)^{2}&=&\sin ^{2}\phi \end{array}}\right\}\Rightarrow A^{2}[(1-\omega ^{2})^{2}+(2\zeta \omega )^{2}]=1.

следи,

A=A(\zeta ,\omega )={\text{sign}}\left({\frac {-\sin \phi }{2\zeta \omega }}\right){\frac {1}{\sqrt {(1-\omega ^{2})^{2}+(2\zeta \omega )^{2}}}}.

Споредувајќи го овој резултат со теорискиот дел за резонанса. Оваа функција на замавот е особено важна за анализа и разбирање на фрквентниот запис на системите од втор ред.

Фазен дел

Решението за φ, се делат двете равенки и се добива:

\tan \phi =-{\frac {2\zeta \omega }{1-\omega ^{2}}}={\frac {2\zeta \omega }{\omega ^{2}-1}}\Rightarrow \phi \equiv \phi (\zeta ,\omega )=\arctan \left({\frac {2\zeta \omega }{\omega ^{2}-1}}\right).

Оваа фазна функција е особено важна за анализа и разбирање на фрквентниот запис на системите од втор ред.

Целосно решение

Со комбинирање на деловите од замавното и фазното решение:

\,\!q_{s}(\tau )=A(\zeta ,\omega )\cos(\omega \tau +\phi (\zeta ,\omega ))=A\cos(\omega \tau +\phi ).

Решението на равенката на универзалниот треперник е суперпозиција (збир) на решенијата на премин и стабилна состојба

\,\!q(\tau )=q_{t}(\tau )+q_{s}(\tau ).

Истоветни системи

Хармониските третерници се присутни во голем дел на области на инженерството и се истоветни во смисла на нивните математички модели кои се идентични (Погледај равенка на универзален треперник од погоре). Подолу има табела која ги прикажува истоветните записи на четири хармониски треперници во механиката и електрониката. Ако истоветните параметри во самата линија на табелата им се придодадат бројчено еднакви вредности, однесувањето на треперниците&mdash, нивната резултантна бранова форма, резонантна фрквенција, придушниот факторитн.tc.— се исти.

Транслационо механички	Вртливо механички	Сериски сврзани RLC кола	Паралелно сврзани RLC кола
Местоположба $x\,$	Агол $\theta \,\!$	Полнеж $q\,$	тек $\phi \,$
Брзина ${\frac {\mathrm {d} x}{\mathrm {d} t}}\,$	Аголна брзина ${\frac {\mathrm {d} \theta }{\mathrm {d} t}}\,$	Струја ${\frac {\mathrm {d} q}{\mathrm {d} t}}\,$	Напон ${\frac {\mathrm {d} \phi }{\mathrm {d} t}}\,$
Маса $M\,$	Момент на инерција $I\,$	Индукција $L\,$	Капацитет $C\,$
Постојана на пружината $K\,$	Постојана на вртењето $\mu \,$	Еластичност $1/C\,$	Сусцептанса $1/L\,$
Придушување $\gamma \,$	Триење при вртењето $\Gamma \,$	Отпор $R\,$	Спроводливост $G=1/R\,$
Присилна сила $F(t)\,$	Присилен момент на сила $\tau (t)\,$	Напон $e\,$	Струја $i\,$
Непридушена резонантна честота $f_{n}\,$ :
${\frac {1}{2\pi }}{\sqrt {\frac {K}{M}}}\,$	${\frac {1}{2\pi }}{\sqrt {\frac {\mu }{I}}}\,$	${\frac {1}{2\pi }}{\sqrt {\frac {1}{LC}}}\,$	${\frac {1}{2\pi }}{\sqrt {\frac {1}{LC}}}\,$
Диференцијална равенка:
$M{\ddot {x}}+\gamma {\dot {x}}+Kx=F\,$	$I{\ddot {\theta }}+\Gamma {\dot {\theta }}+\mu \theta =\tau \,$	$L{\ddot {q}}+R{\dot {q}}+q/C=e\,$	$C{\ddot {\phi }}+G{\dot {\phi }}+\phi /L=i\,$

Примена на конзервативна сила

Проблемот на едноставниот хармониски треперник често се сретнува во физиката, бидејќи маса во рамнотежа под влијание на некаква конзервативна сила, во границите на малите движења, се однесува како едноставен хармониски треперник.

Конзервативна сила е онаа која има функција на потенцијална енергија. Функцијата на потенцијалната енергија на хармонискиот треперник е:

V(x)={\frac {1}{2}}kx^{2}

Земајќи преодна потенцијална енергетска функција $V(x)$ , може да се развие во Тејлоров ред како $x$ околу енергетски минимум ( $x=x_{0}$ ) за да се измоделира однесувањето на малите растројувања од рамнотежната состојба.

V(x)=V(x_{0})+(x-x_{0})V'(x_{0})+{\frac {1}{2}}(x-x_{0})^{2}V^{(2)}(x_{0})+O(x-x_{0})^{3}

Бидејќи $V(x_{0})$ е минимум, првиот извод се проценува кога $x_{0}$ е нула, па линискиот запис се отстранува:

V(x)=V(x_{0})+{\frac {1}{2}}(x-x_{0})^{2}V^{(2)}(x_{0})+O(x-x_{0})^{3}

Постојаната V(x₀) е продна и мора да се отстрани, и координатна трансформација овозможува да се довбие записот на едноставниот хармониски треперник:

V(x)\approx {\frac {1}{2}}x^{2}V^{(2)}(0)={\frac {1}{2}}kx^{2}

Така, со преодна функција на потенцијална енергија $V(x)$ со постоечки втор извод, може да се уопотреби као решение за едноставниот хармониски треперник и приближно решение за малите растројувања околу рамнотежната точка.

Примери

Математичко нишало

Математичкото нишало опишува едноставно хармониско движење при што имаме непридушено движење и мали замави.

Ако имаме непридушено движење и мали замави, диференцијалната равенка која го опишува математичкото нишало е:

{\mathrm {d} ^{2}\theta  \over \mathrm {d} t^{2}}+{g \over \ell }\theta =0.

Решението на оваа равенка е:

\theta (t)=\theta _{0}\cos \left({\sqrt {g \over \ell }}t\right)\quad \quad \quad \quad |\theta _{0}|\ll 1

каде $\theta _{0}$ е најголемиот агол кој го постигнува нишалото.Периодот, времето потребно за да се исполни едно треперење, запишан со $2\pi$ поделен со времето преку косинусот ( ${\sqrt {g \over \ell }}$ here).

T_{0}=2\pi {\sqrt {\ell  \over g}}\quad \quad \quad \quad |\theta _{0}|\ll 1.

Систем пружина/маса

Системот на маса и пружина во рамнотежа (A), натисната (B) и извлечена (C) положба.

Кога пружината е истегната или натисната под дејство на маса, во пружината се јавува еластична сила. Хуковиот закон ја дава поврзаноста меѓу силата која ја поседува пружината кога истата е натисната или развлечена за извесна должина:

F\left(t\right)=-kx\left(t\right)

каде F е силата, k е постојаната на пружината, и x е поместувањето на масата во однос на рамнотежната состојба. Знакот минус во равенката означува дека пружината се придвижува во спротивна насока на поместувањето и се спречува масата да отиде во бесконечноста.

Соо употреба маетодот на урамнотежување на сили или пак еннергетски метод, може постојано да се покаже дека движењето на системот се определува преку следнава равенка:

F(t)=-kx(t)=m{\frac {\mathrm {d} ^{2}}{\mathrm {d} {t}^{2}}}x\left(t\right)=ma.

...последниот дел е Втор Њутнов закоз за движењето.

Ако почетното поместување е A, и немаме почетна брзина, решението на оваа равенка ќе се запише како:

x\left(t\right)=A\cos \left({\sqrt {k \over m}}t\right).

Ако се земе дека имаме идеална пружина без маса , $m$ е масата на крајот од пружината. Ако пак пружината има маса, нејзината делотворна маса може да биде придодадена во $m$ .

Поврзано

Наводи

↑ Ogata, Katsuhiko (2004). System dynamics (4. изд.). Upper Saddle River, NJ: Pearson Education. ISBN 9780131247147.
↑ Ajoy Ghatak (2005). Optics, 3E (3. изд.). Tata McGraw-Hill. стр. 6.10. ISBN 978-0-07-058583-6.
↑ Case, William. „Two ways of driving a child's swing“. Архивирано од изворникот на 2013-09-25. Посетено на 27 November 2011.
↑ doi:10.1119/1.18209
Овој навод ќе се дополни автоматски во текот на следните неколку минути. Можете да го прескокнете редот или да го проширите рачно
↑ Roura, P.; Gonzalez, J.A. (2010). „Towards a more realistic description of swing pumping due to the exchange of angular momentum“. European Journal of Physics. 31 (5): 1195–1207. Bibcode:2010EJPh...31.1195R. doi:10.1088/0143-0807/31/5/020.

[1] Ogata, Katsuhiko (2004). System dynamics (4. изд.). Upper Saddle River, NJ: Pearson Education. ISBN 9780131247147.

[2] Ajoy Ghatak (2005). Optics, 3E (3. изд.). Tata McGraw-Hill. стр. 6.10. ISBN 978-0-07-058583-6.

[Case-3] Case, William. „Two ways of driving a child's swing“. Архивирано од изворникот на 2013-09-25. Посетено на 27 November 2011.

[Case96-4] :10.1119/1.18209
Овој навод ќе се дополни автоматски во текот на следните неколку минути. Можете да го прескокнете редот или да го проширите рачно

[Roura-5] Roura, P.; Gonzalez, J.A. (2010). „Towards a more realistic description of swing pumping due to the exchange of angular momentum“. European Journal of Physics. 31 (5): 1195–1207. Bibcode:2010EJPh...31.1195R. doi:10.1088/0143-0807/31/5/020.

[1]

[2]

[3]

[4]

[5]