Пред.
След.
Макеты страниц
Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике Рассмотрим теперь нелинейную колебательную систему, у которой некоторые параметры, например масса системы, жесткость, коэффициент трения и др., медленно изменяются со временем («медленно» по отношению к естественной единице времени – периоду собственных колебаний). В этом случае мы приходим к рассмотрению следующего нелинейного дифференциального уравнения с медленно меняющимися коэффициентами: в котором, как и ранее, $\varepsilon$-малый положительный параметр, $\tau=\varepsilon t-$ «медленное» время. Построение приближенных решений уравнения (8.1) не вызывает дополнительных принцишиальных затруднений и может быть осуществлено с помощью изложенного нами выше асимптотического метода. Отметим, что для построения асимптотических рядов необходимо, чтобы коэффициенты уравнения (8.1) $m(\tau), k(\tau)$, а также $f\left(\tau, x, \frac{d x}{d t}\right)$ имели достаточное число производных по $\tau$ для всех конечных значений $\tau$, кроме того, для любых $\tau$ на интервале $0 \leqslant \tau \leqslant L$ (только на этом интервале времени $\left(0, \frac{L}{\varepsilon}\right)$ мы будем рассматривать колебательный процесс, описываемый дифференциальным уравнением (8.1)) $k(\tau) Соверпенно очевидно, что если в коэффициентах уравнения (8.1) считать $\tau$ некоторым постоянным параметром, то мы получаем уравнение, рассмотренное нами в § 1, для которого построены асимптотические формулы, учитывающие все те дополнительные явленй, которые возникают в колебательной системе из-за нелинейного возмущения. При $\tau=\varepsilon t$, т. е. если в исследуемой нелинейной колебательной системе некоторые параметры будут изменяться со временем, хотя и медленно, естественно ожидать еще некоторых дополнительных изменений в решении, не наблюдаемых в колебательных системах с постоянной массой, жесткостью и т. Д. Такими дополнительными явлениями будут, например, зависимость «собственной» частоты вт «медленного» времени и др. Учитывая все это, естественно искать общее решение уравнения 8.1) в виде разложения: в котором $u_{1}(\tau, a, \psi), u_{2}(\tau, a, \psi), \ldots$ являются периодическими функциями угла $\psi$ с периодом $2 \pi$, а величины $a$ и $\psi$ как функции времени определяются уже следующими дифференциальными уравнениями: где $\omega(\tau)=\sqrt{\frac{k(\tau)}{m(\tau)}}$ – «собственная» частота рассматривасмой колебательной системы. Таким образом, как и в § 1, задача построения асимптотических приближенных репений уравнения (8.1) сводится к определению выражений для функций и последующему интегрированию системы уравнений (8.3), определяющей амплитуду и полную фазу колебания. Как и в § 1, для однозначного определения функций, стоящих в правых частях уравнений (8.3), необходимо наложить дополнительные условия на функции $u_{1}(\tau, a, \psi), u_{2}(\tau, a, \psi), \ldots$ В качестве таких условий принимаем опять-таки условия (1.8), которые должны в данном случае выполняться для любых г на интервале $0 \leqslant \tau \leqslant L$. После сделанных предварительных замечаний приступим к определению функций (8.4). Для этого дифференцируем правую часть выражения (8.2) с учетом уравнений (8.3) и. подставляем в уравнение (8.1), правую часть которого раскладываем в ряд Тейлора. Приравнивая после этого коәффициенты при одинаковых степенях $\varepsilon$, получаем следующие уравнения: где введены обозначения: Далее, из уравнения (8.5), учитывая условия отсутствия первой гармоники в функции $u_{1}(\tau, a, \psi)$, находим: Но условие отсутствия первой гармоники в (8.9) дает нам следующие выражения для $A_{1}(\tau, a)$ и $B_{1}(\tau, a)$ : Таким образом, в первом приближении асимптотическое решение уравнения (8.1) следует искать в форме где $a$ и должны быть определены из уравнений первого приближения: Заметим, что по аналотии с $\$ 7$ уравнения первого приближения (8.12) можно записать в виде: где $\delta_{e}(\tau, a)$ и $\omega_{e}(\tau, a)$ являются соответственно әквивалентным декрементом затухания и эквивалентной частотой, причем они отличаются от приведенных в § 7 аналогичных выражений только наличием «медленного» времени и определяются выражениями: Для построения второго приближения необходимо найти из условия отсутствия первой гармоники в функции $u_{2}(\tau, a, \psi)$ выражения для $A_{2}(\tau, a)$ п $B_{2}(\tau, a)$. и следовательно, асимптотическое решение уравнения (8.1) во втором приближении будет: где $a$ и $\psi$ должны быть определены из уравнений второго приближения: Здесь $A_{1}(\tau, a)$ и $B_{1}(\tau, a)$ определяются выражениями (8.10), $A_{2}(\tau, a)$ и $B_{2}(\tau, a)$ – выражениями (8.15), а $u_{1}(\tau, a, \psi)$-согласно формуле (8.9). Заметим, что во всех полученных формулах при интегрировании по ечитаем $a$ и $\tau$ постоянными параметрами. Сопоставляя полученные выражения для первого и второго приближения с результатами § 1 , убеждаемся, что общая схема построения решений для случая колебательной системы, описываемой уравнением (8.1), будет такая же, как и в случае системы, рассматривавшейся в § 1. Полученные нами уравнения первого приближения отличаются от уравнений (1.24) наличием «медленного» времени и дополнительного слагаемого Таким образом, в первом приближении медленная изменяемость массы и коэффициента упругости, кроме нарушения гармоничности колебаний, вводит дополнительные «силы трения», знак которых будет зависеть от того, каким образом изменяются параметры исследуемой колебательной системы. Остановимся еще на следующем обстоятельстве. При построении асимштотических приближенных решений мы сводим интегрирование одного уравнения второго порядка (8.1) к интегрированию двух уравнений первого порядка (8.12) или (8.17), которые во многих случаях не могут быть проинтегрированы в әлементарных функциях, и поэтому их решения приходится искать при помощи численных методов. Используя численные методы, можно было бы проинтегрировать непосредственно уравнение (8.1), однако это сложная задача, требующая чрезвычайно много времени и в большинстве случаев трудно выполнимая в связи с возможностью накопления больпой систематической ошибки. Численное же интегрирование уравнений первого (или второго) приближения не представляет затруднений ввиду того, что в этих уравнениях переменными являются амплитуда и фаза, а не сама осциллирующая функция $x$. Для получения полной картины процесса здесь достаточно вычислить небольшое количество точек, расположенных на сравнительно «гладкой» кривой, что существенно упрощает численное интегрирование, в то время как при интегрировании непосредственно уравнения (8.1) нам пришлось бы находить не огибающую, а непосредственно синусоиду. Заметим, что для составления как уравнений первого приближения (8.12), так и второго приближения (8.17) нет необходимости пользоваться формулами (8.10) или (8.15). Мы можем для составления этих уравнений, так же как и в случае уравнения (1.1), воспользоваться непосредственно уравнениями гармонического баланса, которые в натем случае имеют вид Подставляя в подынтегральные выражения значения $x, \frac{d x}{d t}, \frac{d^{2} x}{d t^{2}}$, найденные из (8.2) с точностью до величин первого порядка малости, с учетом, разумеется, того, что $a$ и $\psi$ являются функциями времени, удовлетворяющими уравнениям (8.3), и производя интегрирование, получим для $A_{1}(\tau, a)$ и $B_{1}(\tau, a)$ выражения (8.10). Учитывая при подстановке $x, \frac{d x}{d t}, \frac{d^{2} x}{d t^{2}}$ также величины, пропорциональные $\varepsilon^{2}$, получим для $A_{2}(\tau, a)$ и $B_{2}(\tau, a)$ выражения (8.15). Заметим, что при вычислении интегралов (8.18) $\tau$, как и $a$, считаем некоторым постоянным параметром и интегрирование ведем только по $\psi$. Рассмотрим еще некоторые частные случаи уравнения (8.1), для которых уравнения первого приближения (8.12) принимают совсем простую.форму. Допустим, что в исследуемой колебательной системе отсутствует трение, и уравнение, описывающее движение, имеет следующий вид: В этом случае уравнениями первого приближения будут: Система уравнений (8.20) может быть проинтегрирована до конда. Действительно, из первого уравнения находим: где $a_{0}$ – начальное значение амплитуды при $t=0$. Подставляя найденное значение амплитуды во второе уравнение системы (8.20), получим: где Таким образом, в первом приближении колебания, описываемые уравнением (8.19), будут «синусоидальными» с амплитудой, обратно пропорциональной $\sqrt{m(\tau) \omega(\tau)}$, и фазой, изменяющейся согласно формуле (8.22). В качестве второго частного случая рассмотрим систему с медленно меняющейся массой $m(\tau)$, находящуюся под воздействием линейной упругой силы $k(\tau) x$ с медленно меняющимся коэффициентом упругости и нелинейного трения, зависящего от скорости и \”медленного» времени. Колебания этой системы описываются уравнением В этом случае уравнения первого приближения принимают следующий вид: Из второго уравнения сразу находим закон изменения полной фазы колебания Из системы (8.25) следует, что частота колебания, ошисываемого уравнением (8.24), в первсм приближении не зависит от амплитуды, а зависит только от характера медленного изменения коэффициентов $m(\tau)$ и $k(\tau)$. В качестве примера рассмотрим колебания математического маятника постоянной массы (в случае переменной массы дополнительных затруднений не возникает) при наличии малого затухания, пропорционального первой степени скорости, и медленного изменения длины. $К$ подобной схеме приводятся многочисленные задачи практики. Обозначая $\theta$ угол отклонения маятника от вертикального положения, $g$ – ускорение силы тяжести, $m$ массу маятника, $l=l(t)$ – медленно изменяющуюся длину, $2 n$-коэффициент трения, получим дифференциальное уравнение Для небольших отклонений мы можем $\sin \theta$ заменить первыми двумя членами разложения в степенной ряд, после чего уравнение (8.27) можно записать следующим образом: в котором В первом приближении согласно (8.11) и (8.12) имеем: где $a$ и должны быть определены из системы уравнений первого приближения: где $\omega(\tau)=\sqrt{\frac{g}{l(\tau)}}$ Подставляя это значение $a$ во второе уравнение системы (8.31), найдем: Формулы (8.32) и (8.33) дают возможность построить график зависимости амплитуды и фазы от времени при медленном изменении длины маятника по произвольному закону. где $\lambda=\frac{2 n}{m l}, \varphi$ – начальное значение фазы. Допустим теперь, что длина маятника изменяется по линейному закону $l(\tau)=l_{0}+l_{1} \tau ; l_{0}$ – значение длины при $t=0$, $l_{1}$ – скорость изменения длины маятника (для малого интервала времени всегда можно предположить с достаточной степен:ю точности, что длина изменяется по линейному закону). В этом случае для амплитуды и фазы имеем выражения: Согласно формуле (8.35) амплитуда колебаний при медленном изменении длины маятника будет изменяться не по экспоненциальному закону, как при обычном линейном трении, а обратно пропорционально некоторой степенной функции времени. При этом очевидно, что при $n<0, l_{1}>0$ и $\left|\frac{n}{m l_{1} \varepsilon}\right|<\frac{3}{4}$, а также при $n>0$ и $l_{1}>0$ колебания за+ тухают. Таким образом, медленное увеличение дливы маятника, как и следовало ожидать, способствует затуханию колебаний. Если $l_{1}<0, n>0$ и $\left|\frac{n}{l_{1} m_{s}}\right|<\frac{3}{4}$, то амплитуда возрастает, а при $\left|\frac{n}{l_{1} m_{\varepsilon}}\right|>\frac{3}{4}-$ убывает. При $l_{1}<0$ и $n<0$ амплитуда возрастает. При отсутствии затухания $(n=0)$ амплитуда колебаний возрастает с уменышением длины и убывает с увеличением длины. Аналогичный анализ может быть проведен и для частоты колебаний. Так, например, при отсутствии затухания мгновенная частота уменьшается при увеличении длины маятника и увеличивается при уменьшении длины. Подсчитаем для данного примера колебаний маятника с медленно меняющейся длиной второе приближение. Согласно формулам (8.16), (8.9) и (8.17) после ряда выкладок имеем: где $a$ и $\psi$ должны быть определены из системы уравнений второго приближения: которые тоже могут быть проинтегрированы до конца. Так, из первого уравнения системы (8.38) получаем следующее соотношение между $a$ и $t$ : после чего можем проинтегрировать также и второе уравнение системы (8.38).
|
1 |
Оглавление
|