ності
, n-параметр опору середовища
F (t)=F0? ? 2 (sin? T + 2sin 10? T), F (t) - збурювальна сила
Закон має вигляд:
(2.2)
тут - амплітуда сили, що обурює
?- Частота збуджуючої сили.
.4 Алгоритм розрахунку базової моделі і проведення досліджень
Так
Рисунок 2.2 - Графічна схема алгоритму
На даній графічній схемі (рис. 2.4) представлено короткий опис рішення задачі в системі MathCAD. Першим пунктом графічної схеми є введення вихідних даних: m-маса вантажу;
m - маса вантажу;
с - жорсткість пружин;
?- Коефіцієнт демпфірування;
?- Частота збуджуючої сили.
Далі записуємо диференціальне рівняння.
Наступним пунктом є рішення диференціального рівняння при заданих початкових умовах. Будуємо графіки переміщення, швидкості і прискорення.
Змінюючи величину частоти коливань проводимо 10 дослідів. Далі будуємо зведений графік залежності мінімуму переміщення від маси.
Побудувавши зведений графік, необхідно вказати явище резонансу.
Останнім пунктом даної графічної схеми є визначення діапозона значень частоти.
Графічна схема вирішення даної задачі наведена нижче.
matlab моделювання оператор матричний
. Опис реалізації в пакеті MathCAD
.1 Опис реалізації базової моделі
Документ MathCAD наведений у додатку A. У пакеті MathCAD здійснюється наступне:
Завдання 1.1
) Введення вихідних даних для вирішення завдання: з-жорсткість пружини. - маса вантажу; ?- Коефіцієнт демпфірування;
? - частота збуджуючої сили.
Числові значення вихідних даних кожного з параметрів беруться з п.2.2 (згідно виданим 14-му варіанту 2 завдання на курсову роботу)
) Завданням вектора початкових умов у, вектора-функції, що містить похідні D (t, y) без урахування сили, що обурює і рішення системи диференціальних рівнянь за допомогою функції rkfixed при значеннях граничних точок інтервалу, на якому шукається рішення диференціальних рівнянь від 0 до 5 в 1000 точках (не рахуючи початкової точки), в яких шукається наближене рішення. Результатом є матриця z значень часу (стовпець Z <1>), переміщення (стовпець Z <2>) і швидкості (стовпець Z <3>), що містить 1001 рядок (рахуючи початкову точку) (див. Додаток А).
) визначення (використовуючи матрицю z) значенні прискорення в 1001 точці (див. Додаток А).
) Побудова графіків функцій переміщення (див. додаток А рис.1.1), швидкості (див. додаток А рис.1.2) і прискорення (див. додаток А рис.1.3) залежно від часу (см . Додаток А).
Завдання 1.2
) Задаємо матрицю? 1 (значення частоти) з 10 значеннями змінних.
) Задаємо вектор початкових умов y, век...