Для перевірки вірності визначення коефіцієнту
рекомендується підрахувати значення співмножника в рішенні при =5.0,861 = 4,31с: .Таке значення співмножника (наближене до нуля) в рішенні
підтверджує факт, що вільні коливання системи на цей момент часу практично згасають; значить коефіцієнт знайдено вірно.3. Визначення амплітудних- та фазово-частотних характеристик системи.
Шляхом виведення, за допомогою ЕОМ, для заданої механічної системи з параметрами
= 0,256т; = 13,6кН.м –1; = 0,456кН.с.м–1 получимо (шляхом введення на друкарський пристрій – принтер) амплітудно- та фазово-частотніх характеристики системи та приведемо їх на рис.2 і рис.3 (відповідно).4. Розкладання функції F(t) в ряд Фур’є та визначення параметрів гармонік збурюючої сили.
Розкладемо функцію
в ряд Фур’є: , ( )де
- номер гармоніки, а - число гармонік в розкладенні.Визначимо (за допомогою ЕОМ) параметри гармонік: амплітуди
, частоти та початкової фази .Для заданої сили “прямокутного” типу з параметрами
кН, значення параметрів гармонік наведені у табл.1.Таблиця 1.
Номер гармоніки, | , кН | , | , рад. |
1 | 0,764 | 2 | 0 |
2 | 0,255 | 6 | 0 |
3 | 0,153 | 10 | 0 |
4 | 0,109 | 14 | 0 |
5 | 0,085 | 18 | 0 |
5. Дослідження вимушених коливань механічної системи.
5.1. Визначення (за допомогою ЕОМ) “точного” рішення диференціального рівняння. Аналіз рішення.
Визначимо за допомогою ЕОМ “точне” рішення
диференціального рівняння для випадку, коли сила представлена однією гармонікою ( =1). Два графіка функцій для відповідних випадків виводяться на екран ЕОМ. Перед виводом графіків на друкарський пристрій їх треба “промасштабувати”, тобто получити рішення на заданому відрізку інтегрування 0 , де рекомендується задавати рівним 8 10 . На рис.4 приведені вказані графіки функцій для заданої механічної системи. Лінія 1 відображає “точне” рішення, а лінія 2 – рішення у випадку = 1 (тобто, коли ).Із графіків видно, що функції
получаються періодичними, тобто рух механічної системи получається періодичним-коливальним. І в першому, і в другому випадку при с процеси получилися явно усталеними, тобто без вільних коливань. Але в “точному” рішенні навіть при явно виражені дві частоти – одна дорівнює (див. лінію 2 для випадку = 1; це частота першої гармоніки), а друга частота – втричі більша (це частота другої гармоніки, див. табл. 1). Рішення, що відповідають лініям 1 і 2, значно відрізняються одне від одного. Наприклад, в момент часу с ( =4,31с) значення м (точне рішення) і м (випадок = 1).5.2. Підбір (за допомогою ЕОМ) раціональної кількості гармонік в розкладанні функції .
Визначимо (за допомогою ЕОМ) функції
для випадків = 2, 3. На рис. 5, а, б лініями 1 показані графіки для “точного” рішення, а лініями 2 – графіки тих же функцій для випадків = 2; = 3 (відповідно). Із рисунків випливає, що при = 2 графік відрізняється від “точного” рішення, але в значно меншій степені, ніж при = 1. А у випадку = 3 графік практично не відрізняється від “точного” рішення.Значення відповідних функції при
с становлять м ( = 2) і м ( = 3), тобто при = 2 різниця в значеннях відповідає D = 5,7%, а при = 3 - D = 3,7%.За одержаним результатам можна зробити висновок, що для отримання рішення
з 5% точністю достатньо взяти кількість гармонік = 3 в розкладенні збурюючої сили в ряд Фур’є.