- •Студентові
- •1. Цикл "Нелінійні властивості фізичного маятника" Мета циклу
- •Склад циклу
- •Загальна теоретична частина
- •Опис програмних засобів
- •1.1.2. Сценарій роботи
- •Xlabel('Час (в періодах малих власних коливань)','FontSize',12);
- •Xlabel('Кут у градусах','FontSize',12); ylabel('Швидкість','FontSize',12);
- •1.1. Залежність частоти власних коливань маятника від їх амплітуди
- •1.1.1. Мета роботи
- •1.1.3. Теоретична частина. Неізохронність власних коливань. Скелетні криві
- •Xlabel('Відносна частота власних коливань')
- •Xlabel('Відносна частота власних коливань')
- •1.1.4. Контрольні запитання
- •1.2. Амплітудно-частотні характеристики маятника
- •1.2.1. Мета роботи
- •1.2.2. Сценарій роботи
- •1.2.3. Теоретична частина. Амплітудно-частотні характеристики маятника
- •Xlabel('Відносна частота збудження'), ylabel('Амплітуда (градуси)')
- •1.2.4. Завдання
- •1.2.5. Методичні вказівки до виконання
- •1.2.6. Контрольні запитання
- •1.3. Випрямний ефект маятника при вібрації точки його підвісу
- •1.3.1. Мета роботи
- •1.3.2. Сценарій роботи
- •1.3.3. Теоретична частина. Випрямний ефект
- •1.3.4. Завдання
- •1.3.5. Контрольні запитання
- •1.4. Параметричні коливання маятника при вертикальній вібрації точки його підвісу
- •1.4.1. Мета роботи
- •1.4.2. Сценарій роботи
- •1.4.3. Теоретична частина. Параметричні коливання маятника
- •1.4.4. Завдання
- •1.4.5. Методичні вказівки до виконання
- •1.4.6. Контрольні запитання
- •2. Цикл "Динаміка гіроскопів" Вступ
- •2.1. Поводження симетричного зрівноваженого гіроскопа
- •2.1.1. Мета роботи
- •2.1.2. Сценарій роботи
- •2.1.3. Теоретична частина
- •2.1.3.1. Вільний рух сзг
- •2.1.3.2. Рух сзг під дією сталого моменту сил
- •2.1.3.3. Вплив тертя на рух сзг
- •2.1.4. Опис програмних засобів
- •2.1.5. Завдання
- •2.1.6. Методичні вказівки до виконання
- •2.1.7. Контрольні запитання
- •2.2. Поводження гіромаятника
- •2.2.1. Мета роботи
- •2.2.2. Сценарій роботи
- •2.2.3. Теоретична частина
- •2.2.3.1. Вільний рух гіромаятника
- •2.2.3.2. Вплив тертя на поводження гіромаятника
- •2.2.4. Опис програмних засобів
- •Xlabel('Час (с)','Fontsize',12)
- •2.2.5. Завдання
- •2.2.6. Методичні вказівки до виконання
- •2.2.7. Контрольні запитання
- •2.3. Поводження гіроскопа у кардановому підвісі
- •2.3.1. Мета роботи
- •2.3.2. Сценарій роботи
- •2.3.3. Теоретична частина
- •2.3.3.1. Вільний рух гкп
- •2.3.3.2. Вплив тертя на вільний рух гкп
- •2.3.3.3. Рух гкп під дією сталих моментів сил
- •2.3.3.4. Вплив тертя на рух гкп під дією сталих моментів сил
- •2.3.3.5. Рух гкп під дією гармонічно змінюваних моментів сил
- •2.3.4. Опис програмних засобів
- •If isempty(sname)
- •2.3.5. Завдання
- •2.3.6. Методичні вказівки до виконання
- •2.3.7. Контрольні запитання
- •3. Цикл "Динаміка приладів орієнтації на рухомій основі" Вступ
- •3.1. Поводження двоступеневого гірокомпаса на рухомій основі
- •3.1.1 Мета роботи
- •3.1.2. Теоретична частина
- •3.1.2.1. Вільний рух дгк
- •3.1.2.2. Рух дгк при кутовій вібрації
- •3.1.3. Опис програмного забезпечення
- •3.1.4. Завдання
- •3.1.5. Методичні вказівки до виконання
- •3.1.6. Контрольні запитання
- •3.2. Поводження гіроскопа напрямку на рухомій основі
- •3.2.1 Мета роботи
- •3.2.2 Теоретична частина
- •3.2.2.1. Гкп на обертовій основі під впливом моментів сил тертя
- •3.2.2.2. Дрейфи гкп при хитавиці основи
- •3.2.2.3. Віражна похибка гіроскопу напрямку
- •3.2.2.4. Похибка гн за хитавиці
- •3.2.3. Опис програмного забезпечення
- •Gn_po_precession
- •3.2.4. Завдання
- •3.2.5. Методичні вказівки до виконання
- •3.2.6. Контрольні запитання
- •3.3. Поводження гіромаятника на рухомій основі
- •3.3.1. Мета роботи
- •3.3.2 Теоретична частина
- •3.3.2.1. Вільний рух гіромаятника
- •3.3.2.2. Балістична похибка гм
- •3.3.3. Опис програмного забезпечення
- •3.3.4. Завдання
- •3.3.5. Методичні вказівки до виконання
- •3.3.6. Контрольні запитання
- •3.4. Поводження маятникового гірокомпаса на рухомій основі
- •3.4.1. Мета роботи
- •3.4.2 Теоретична частина
- •3.4.2.1. Вільний рух маятникового гірокомпаса
- •3.4.2.2. Балістична похибка. Умова незбурності
- •3.4.3. Опис програмного забезпечення
- •Xlabel('Час, хвилини','Fontsize',12), legend('швидкісна',' балістична',0)
- •3.4.4. Завдання
- •3.4.5. Методичні вказівки до виконання
- •3.4.6. Контрольні запитання
- •Додаток. Зразок звіту з лабораторної роботи
- •1. Мета лабораторної роботи
- •2. Лабораторна установка
- •3. Теоретичні відомості
- •4. Завдання 1
- •4.1. Результати експериментального дослідження
- •4 .2. Аналіз і висновки
- •Xlabel('Період t_0, хвилини')
- •5. Завдання 2
- •5.1. Результати експериментального дослідження
- •5.2. Аналіз і висновки
- •Загальні висновки
- •3. Цикл "Динаміка приладів орієнтації на рухомій основі" 84
- •3.1. Поводження двоступеневого гірокомпаса на рухомій основі 85
- •3.2. Поводження гіроскопа напрямку на рухомій основі 96
- •3.3. Поводження гіромаятника на рухомій основі 116
- •3.4. Поводження маятникового гірокомпаса на рухомій основі 128
1.1.3. Теоретична частина. Неізохронність власних коливань. Скелетні криві
З власного досвіду нам відомо, що період коливань гойдалки залежить від їх розмаху: якщо амплітуда розгойдування наближається до 90 градусів, період коливань помітно збільшується. Таку властивість мають лише нелінійні системи (в лінійних стаціонарних системах, поводження яких можна описати лінійними диференційними рівняннями зі сталими коефіцієнтами, частоти власних коливань залишаються тими самими за будь-яких значень амплітуди цих коливань). Її прийнято називати неізохронністю власних коливань.
Щоб вивчити теоретично цю властивість, розглянемо вільні незгасаючі коливання ФМ, коли точка його підвісу є нерухомою, а коефіцієнт демпфірування дорівнює нулю. У цих умовах рівняння (0.4) маятника набуде вигляду
. (1.1)
Спершу вивчимо поводження ФМ, користуючись
методом лінеаризації, задля чого,
припускаючи, що ФМ в умовах вільного
руху здійснює симетричні відносно
вертикалі незгасальні коливання зі
сталою амплітудою
,
наближено замінимо нелінійний член
на лінійний
. (1.2)
Величина
у загальному випадку є функцією амплітуди
коливань. Конкретний вигляд цієї функції
визначається особливостями обраного
методу лінеаризації. Якщо розглядувані
коливання є стаціонарними (тобто їхня
амплітуда не змінюється з часом), то
є сталою величиною, і відповідне
лінеаризоване рівняння
(1.3)
стає рівнянням зі сталими коефіцієнтами, розв’язок якого легко відшукується відомими способами:
. (1.4)
Значення амплітуди
і початкової фази
власних коливань визначаються (як сталі
інтегрування) початковими умовами. Не
зменшуючи узагальненості результатів,
задля спрощення вважатимемо, що у
початкову мить початкова кутова швидкість
дорівнює нулю, тобто маятник відхилено
на початковий кут
від вертикалі і відпущено без поштовху.
Тоді співвідношення (1.4) дещо спрощується
, (1.5)
причому амплітуда тепер дорівнює початковому відхиленню маятника.
Величина
(1.6)
має сенс частоти власних коливань, точніше, враховуючи безрозмірну форму первісного рівняння, – відношення частоти власних коливань маятника до значення тієї ж частоти при прагненні амплітуди коливань до нуля (малих власних коливань). Як бачимо, згідно теорії, дійсно, частота власних коливань нелінійної системи може залежати від їхньої амплітуди.
Конкретизуємо залежність (1.6).
Спочатку використаємо лінеаризацію
по кінцях. Згідно з нею, відрізок
нелінійної залежності
,
що міститься всередині інтервалу
,
замінюється на пряму, що проходить крізь
точки цього відрізку на кінцях інтервалу
(тобто, через точки
і
).
Тому у цьому випадку коефіцієнт
лінеаризації
,
який визначає тангенс куту нахилу прямої
до осі абсцис, дорівнюватиме
. (1.7)
Більш складною виявляється аналогічна формула в прямій лінеаризації:
. (1.8)
Через інтегрування частинами цей вираз для розглядуваної нелінійності (1.2) приводить до результату
. (1.9)
В теорії нелінійних коливань графік залежності частоти власних коливань від амплітуди, побудований у такий спосіб, що амплітуду відкладено вздовж осі ординат, а частоту – вдовж осі абсцис, називають скелетною кривою.
Побудуємо скелетні криві ФМ відповідно до формул (1.7) і (1.9), використовуючи можливості системи MatLAB:
a=pi/1000:pi/100:4.5;
lam1=sqrt(sin(a)./a);
lam2=sqrt(3*(sin(a)-a.*cos(a))./a.^3);
plot(lam1,a,lam2,a,'.'), grid
title('Скелетні криві ФМ')
