
- •Розділ 1. Загальні питання теорії моделювання
- •1.1. Поняття моделі в прикладному дослідженні
- •1.2. Моделювання і математика
- •1.3. Обчислювальний експеримент
- •Розділ 2. Математичні моделі інженерних задач
- •2.1. Лінійні алгебраїчні моделі
- •2.1.1. Задача про розподіл токів у електричному ланцюзі
- •2.1.2. Системи лінійних алгебраїчних рівнянь
- •2.1.3. Коректність, стійкість, достовірність
- •2.1.4. Методи рішення лінійних систем на еом
- •2.1.5. Вправи
- •2.2. Моделі, що зводяться до задач про власні значення
- •2.2.1. Власні вектори і власні значення матриць
- •2.2.2. Методи рішення задач про власні значення на еом
- •2.2.3. Вправи
- •2.3. Нелінійні алгебраїчні та трансцендентні моделі
- •2.3.1. Задача про розмір критичної сили для стрижня
- •2 .3.2. Корені нелінійних рівнянь і систем
- •2.3.3. Методи рішення нелінійних рівнянь і систем
- •2.3.4. Вправи
- •2.4. Моделі, що зводяться до диференціальних рівнянь.
- •2.4.1. Математичні моделі механічних систем
- •2.4.2. Математичні моделі для електричних ланцюгів
- •2.4.3. Моделі систем автоматичного регулювання
- •2.4.4. Математичні моделі вагонів
- •2.4.5. Рішення звичайних диференціальних рівнянь
- •Метод Ейлера
- •Метод Рунге-Кутта
- •2.4.7. Вправи
2.4.7. Вправи
Вправа 8. З задачами про удар і коливання часто доводить зштовхуватися в аерокосмічній промисловості і на транспорті, де є численні джерела порушення коливань. Усунення ударних і вібраційних навантажень має винятково велике значення для забезпечення нормальної роботи приладів і систем керування і створення комфортних умов для екіпажа.
Звичайно для захисту від надмірних вібрацій у конструкції транспортного засобу вводять пружні опори постачені пристроями, що забезпечують деяке демпфування коливань. Такі опори різко зменшують частоти власних коливань конструкції, забезпечуючи їхню істотну відмінність від частот силових чинників , що обурюють. Таке рішення ефективно як засіб захисту від стаціонарних коливань, проте у випадку ударних навантажень податливість опор може призвести до неприпустимо великих зсувів.
Відомо, що від цієї хиби вільні системи підвіски, у яких використовуються пружини із симетричною нелінійною характеристикою, жорсткість котрих пргресивно збільшується при великих відхиленнях від ”робочої точки”.
Пристрій, показаний на малюнку , складається з маси m , пов'язаної з жорсткою стінкою через пружину постійної жорсткості k, демпфер із коэфициентом демпфировання c і пружиною з нелінійною характеристикою, що створює силу опору, рівну добутку постійної k* на зсув у третьому ступені.
x
с
m
к
к*
Тоді рух розглянутої системи описується лінійним диференціальним рівнянням :
mx + сx + kx +k*x3 = 0 ,
і традиційні ”точні” методи не дозволяють знайти залежність зсуву х від часу. Тому для рішення зазначеного рівняння прийдеться удатися до чисельного методу.
Нехай параметри системи имеют наступні значения:
K= 2,0 н / cм ; k* = 2,0 н/см3 с = 0,15 н c/ см ; m = 1кг, а початкові умови задані у виді
x(0) = 10 см , x ‘ (0) = 0 .
Проведемо моделювання руху розглянутої механічної системи в інтервалі часу 0 t 1 c.
Щоб вирішити цю задачу одношаговим чисельним методом, необхідно звести диференціальне рівняння другого порядку до двох диференціальних рівнянь першого порядку. У результаті отримаєм:
Щоб довжина виражалися в сантиметрах, змінимо розмірність маси m =0,01 Н с2/див. Вирішити задачу методом Рунге-Кутта.
Вправа 9. Тіло, зображене на малюнку, рухається по плоскій поверхні з тертям, що обумовлює демпфирування коливань. Його маса m=4,5 кг, жорсткість пружини k=175 H/м, коефіцієнт тертя f=0,3. Розрахувати рух тіла в інтервалі часу 0<=t<=2 із початковими умовами x(0)=7,5см, x ‘ (0)=0. Розрахунок виконати по методу Рунге-Кутта, результат навести графічно.
X(t)
K m
Вправа 10. Диференціальне рівняння згибної лінії бруска з постійним поперечним перетином має вигляд
,
при
початкових умовах
y(0)=0,
y’(0)=0
і
значеннями L=1м,
PL2/EJ=2.
Вправа 11. Матеріальна точка масою m закріплена на кінцях однорідного стрижня і чинить вільні коливання. Промоделювати рух цієї системи і виконати розрахунок для повного циклу коливань.
Примітка. Скористатися наступним рівнянням руху, де незалежної перемінної є час:
Обчислити задачу методом Рунге-Кутта.
11Студентам впевнитися у наведеному результаті самостійно.