1612042555-2f73713c5a14598eee34448f2d1c6445 (542294), страница 22
Текст из файла (страница 22)
А. (НГУ)ЛЕКЦИЯ 14Новосибирск, 2017 г.14 / 21Далее рассмотрим сначала случай, когда все корни λσ уравнениячастот – различные: λσ 6= λρ , ∀ σ 6= ρ (тогда различны и ūσ ).Т.к. исходная система дифференциальных уравнений движения –линейная, то линейная комбинация её частных решений спостоянными коэффициентами Cσ – является решением системы:q̄ =nXCσ ūσ sin(ωσ t + ασ )σ=1⇔q̄ =nXCσ q̄ σ(∗)σ=1при произвольных постоянных Cσ и ασ является решениемисходной системы. Покажем, что эта формулаохватывает все движения системы, т.е. является общим решением.Т.е. удовлетворяет не только дифференциальным уравнениямдвижения, что мы уже показали, но и любым начальнымусловиям (т.е.
Cσ и ασ (σ = 1, . . . , n) можно так выбрать).Батяев Е. А. (НГУ)ЛЕКЦИЯ 14Новосибирск, 2017 г.15 / 21Предварительно установим следующее свойство:Свойство 1.Aūσ ūρ = 0,σ 6= ρбилинейная форма с матрицей A для различных амплитудныхвекторов обращается в ноль.Действительно, имеем для различных σ 6= ρ, т.е. по предположениюλσ 6= λρ , равенства:(C − λσ A)ūσ = 0 | · ūρC ūσ · ūρ = λσ Aūσ · ūρ⇒(C − λρ A)ūρ = 0 | · ūσC ūρ · ūσ = λρ Aūρ · ūσУчитывая доказанное ранее свойство симметричных матриц онезависимости соответствующих им билинейных форм от порядкаумножения векторов:Aū · v̄ = Av̄ · ū,C ū · v̄ = C v̄ · ūПосле вычитания имеем:(λσ − λρ )Aūσ · ūρ = 0откуда в виду λσ 6= λρ получим требуемое свойство.Батяев Е.
А. (НГУ)ЛЕКЦИЯ 14Новосибирск, 2017 г.16 / 21Замечание: билинейная форма – если ū, v̄ – разные: Aū · v̄квадратичная форма – если ū, ū – одинаковые: Aū · ūСвойство 2. Теперь покажем, что в рассматриваемом случае,когда корни уравнения частот различные, n амплитудных векторовūσ – линейно независимы.Действительно, пусть b̄ =nPσ=1Bσ ūσ = 0Покажем что тогда все коэффициенты Bσ = 0 (σ = 1, .
. . , n)При любом фиксированном ρ будем иметь:Ã n!nXXAūρ ·b̄ = 0 ⇒ Aūρ ·Bσ ūσ =Bσ Aūρ ·ūσ = Bρ Aūρ ·ūρ = 0σ=1σ=1Но квадратичная форма Aūρ · ūρ > 0, значит Bρ = 0 для всякого ρ.Что доказывает независимость ūσ и соответствующих им решенийq̄ σ = ūσ sin(ωσ t + ασ )Батяев Е. А. (НГУ)ЛЕКЦИЯ 14Новосибирск, 2017 г.17 / 21Найдём теперь в формуле (∗) значения 2n произвольных постоянныхCσ и ασ (σ = 1, .
. . , n) так, чтобы удовлетворялись произвольныенаперёд заданные начальные условияqσ (0) = qσ0 ,q̇σ (0) = q̇σ0(σ = 1, . . . , n)Из формулы (∗) находим:nXCσ sin ασ ūσq̄ 0 =q̄˙ 0 =σ=1⇔nXq̄(0) = q̄ 0 ,˙q̄(0)= q̄˙ 0Cσ ωσ cos ασ ūσσ=1В силу линейной независимости ūσ из данной системы алгебраическихлинейных уравнений однозначно определяются произведения Cσ sin ασи Cσ ωσ cos ασ , а следовательно, т.к.
ωσ 6= 0, однозначно определяетсяCσ и ασ (ασ определяется с точностью до слагаемого, кратного 2π).Таким образом, при отсутствии в уравнении частот кратных корней,формула (∗) охватывает все движения систем, которые представляютсобой малые колебания. В случае, когда уравнение частот имееткратные (одинаковые) корни, оказывается возможным также найти nнезависимых решений в виде колебаний: ūσ sin(ωσ t + ασ ) ипредставить общее решение системы в форме (∗).Батяев Е. А. (НГУ)ЛЕКЦИЯ 14Новосибирск, 2017 г.18 / 21Замечание: Частоты ωσ – не зависят от начальных данных иполностью определяются матрицами A и C, т.е.
свойствами системы.Это свойство собственных частот малых колебаний называется— изохронность.В заключение отметим, что колебания вида q̄ σ = Cσ ūσ sin(ωσ t + ασ ),из которых складывается произвольное колебание системы называютсяq̄ σ = Cσ ūσ sin(ωσ t + ασ ) — главные колебанияВектор ūσ — амплитудный вектор σ-го главного колебания.Т.е.
в σ-ом главном колебании все обобщённые координаты совершаютгармоничные колебания с одной и той же частотой ωσ ! Отношениеамплитуд колебаний отдельных обобщённых координат определяетсяотношением соответствующих компонент амплитудных векторов.Решение (∗) иногда удобнее представить в другом, альтернативномвиде (без начальных фаз ασ , но с косинусами и синусами):q̄ =nX³´ūσ Cσ(s) sin(ωσ t) + Cσ(c) cos(ωσ t)σ=1Батяев Е. А. (НГУ)ЛЕКЦИЯ 14Новосибирск, 2017 г.19 / 21Отметим, что если в n-мерном координатном пространстве {q1 , . . . , qn }ввести A-метрику (евклидову структуру при помощи удвоеннойкинетической энергии 2T ) — под квадратом длины вектора ūпонимать величину квадратичной формыAū · ū =nXaσρ uσ uρσ,ρ=1т.е.
принять за скалярное произведение векторов ū и v̄ величинуnX(Aū, v̄) =aσρ uσ vρ – билинейную форму, где A – положительноσ,ρ=1определённая матрица, т.е. квадратичная форма всегда положительна:nX(Aū, ū) =aσρ uσ uρ > 0,σ,ρ=1тогда равенство Aūσ ūρ = 0 выражает собой свойствоортогональности амплитудных векторов, соответствующих различнымкорням уравнения частот, в A-метрике.Батяев Е.
А. (НГУ)ЛЕКЦИЯ 14Новосибирск, 2017 г.20 / 21Учитывая однородность уравнений для определения амплитудныхвекторов (C − λσ A)ūσ = 0, умножение на скаляр которого,ничего не меняет, замечаем, что амплитудные вектора ūσопределяются с точностью до постоянного множителя. А этотмножитель, в свою очередь, можно подобрать так, чтобы½0, σ 6= ρ(Aūσ , ūσ ) = 1, т.е. (Aūσ , ūρ ) = δσρ =1, σ = ρгде δσρ – символ Кронекера. Таким образом, можно сделатьамплитудные вектора ūσ — ортонормированными в A-метрике.Замечание: В силу свойства C ū · v̄ = λAū · v̄, одновременно сравенством Aūσ ūρ = 0 имеет место и равенствоC ūσ ūρ = 0(σ 6= ρ)т.е.
ортогональность амплитудных векторов в C-метрике.Батяев Е. А. (НГУ)ЛЕКЦИЯ 14Новосибирск, 2017 г.21 / 21ТЕОРЕТИЧЕСКАЯ МЕХАНИКА2 СЕМЕСТРЛЕКЦИЯ 15НОРМАЛЬНЫЕ КООРДИНАТЫВЛИЯНИЕ ПЕРИОДИЧЕСКИХ ИДИССИПАТИВНЫХ ВНЕШНИХ СИЛ НАКОЛЕБАНИЯ КОНСЕРВАТИВНЫХ СИСТЕМЛектор: Батяев Евгений АлександровичБатяев Е. А. (НГУ)ЛЕКЦИЯ 15Новосибирск, 2017 г.1 / 24Составление общего решения системы дифференциальных уравнений малыхколебанийA ¨q̄ + C q̄ = 0не зависящего от наличия кратных корней уравнения частотdet(C − λA) = 0основывается на приведении квадратичных форм к сумме квадратов.Из линейной алгебры известно, что две квадратичные формыnnXXAq̄ · q̄ =aσρ qσ qρ ,C q̄ · q̄ =cσρ qσ qρσ,ρ=1σ,ρ=1из которых хотя бы одна, например Aq̄ · q̄, является положительноопределенной (у нас они обе положительно определены), всегда можнонеособенным (невырожденным) преобразованием (заменой) переменныхq̄ = U θ̄ =nXūσ θσ⇔qρ =σ=1det U = det{ū1 , .
. . , ūn } 6= 0,привести к сумме квадратов:nXuρσ θσθ̄ = (θ1 , . . . , θn ) − новые переменныеAq̄ · q̄ =nXσ=1Батяев Е. А. (НГУ)(ρ = 1, . . . , n)σ=1ЛЕКЦИЯ 15θσ2 ,C q̄ · q̄ =nXλσ θσ2σ=1Новосибирск, 2017 г.2 / 24Aq̄ · q̄ =nXθσ2 ,C q̄ · q̄ =σ=1nXλσ θσ2σ=1Причем во второй сумме λσ – это корни уравнения частот. В самом деле:nX∗Aq̄ · q̄ = AU θ̄ · U θ̄ = U AU θ̄ · θ̄ =θσ2⇒U ∗AU = I = diag{1}C q̄ · q̄ = CU θ̄ · U θ̄ = U ∗CU θ̄ · θ̄ =σ=1nXλσ θσ2⇒U ∗CU = diag{λσ }σ=1U ∗ – сопряженная к U матрица (в R сопряжение = транспонированию). Тогдаdet(U ∗CU − λU ∗AU ) = det(diag{λσ } − λI) = det(diag{λσ − λ}) =nY(λσ − λ)σ=1с другой стороныdet(U ∗CU − λU ∗AU ) = det(U ∗(C − λA)U ) = det U ∗ · det(C − λA) · det UУчитывая, что det U 6= 0 и det U ∗ 6= 0, имеемnYdet(C − λA) =(λσ − λ) = 0σ=1Батяев Е.
А. (НГУ)ЛЕКЦИЯ 15Новосибирск, 2017 г.3 / 24Поэтому величины λσ , с точностью до порядка следования, являютсякорнями уравнения частот и однозначно определяются исходнымиквадратичными формами, т.е. не зависят от выбора заменыпеременных q̄ = U θ̄. Причем, как мы показывали, все λσ как корниуравнения частот – вещественны и положительны.Относительно векторов ūσ , из которых в виде столбцов составленаматрица преобразования U , можно сказать, что они – линейнонезависимы, поскольку соответствующий определитель отличен отнуля: det U = det kuσρ knσ,ρ=1 6= 0. Более того, оказывается ониортонормированы в A-метрике:Aūσ · ūρ = δσρ(∀ σ, ρ = 1, .
. . , n)В самом деле:ÃAq̄ · q̄ = AnXσ=1Батяев Е. А. (НГУ)! nnnXXXūσ θσ ·ūρ θρ =Aūσ · ūρ · θσ θρ =θσ2ρ=1σ,ρ=1ЛЕКЦИЯ 15σ=1Новосибирск, 2017 г.4 / 24Нормальные (главные) координатыВернемся к преобразованию. Поскольку оно линейное, то обобщённыескорости q̄˙ и θ̄˙ будут связаны аналогичным соотношением:nX˙q̄ = U θ̄˙ =ūσ θ̇σÃпоэтому в первой из формσ=1nXAq̄ · q̄ =!θσ2можно заменить q̄ и θσ ,σ=1соответственно на q̄˙ и θ̇σ . В результате получим следующие выражениядля кинетической и потенциальной энергий в новых переменных:nn11X 2T = Aq̄˙ · q̄˙ =θ̇σ2211XΠ = C q̄ · q̄ =λσ θσ222σ=1σ=1Обобщённые координаты (переменные) θσ называются —нормальными или главными координатами,в которых T и Π принимают канонический вид.Батяев Е. А.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
