Метод определения особых траекторий колебания груза 2d-пружинного маятника Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Семкив О.М.

Разработан метод определения средствами дифференциального исчисления и компьютерной графики геометрической формы особой траектории перемещения по плоскости груза 2d-пружинного маятника в зависимости от массы груза, начальной длины пружины в ненагруженном состоянии, коэффициента жесткости пружины и начальных условий возникновения колебаний.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Семкив О.М.

THE METHOD OF DETERMINING SPECIAL LOAD OSCILLATION TRAJECTORIES OF THE 2D-SPRING PENDULUM

A method for determining a geometric shape of a special trajectory of movement of the 2D-spring pendulum weight on the path by the means of differential calculus and computer graphics, depending on the load weight, the initial length of the spring in the unloaded condition, the coefficient of spring rigidity and the initial conditions of oscillation occurrence is developed.

Текст научной работы на тему «Метод определения особых траекторий колебания груза 2d-пружинного маятника»

МЕТОД ВИЗНАЧЕННЯ ОСОБЛИВИХ ТРАСКТОР1Й КОЛИВАНЬ ВАНТАЖУ 2D-ПРУЖИННОГО МАЯТНИКА

О.М. Семк1в, к.т.н., ст.н.сп., Нацюнальний ун1верситет цивильного захисту УкраУни, м. Харкчв

Анотаця. Розроблено метод визначення засобами диференщального числення та комп'ютерног графти геометричног форми особливог траекторИ перемщення по площит ва-нтажу 2d-пружинного маятника залежно eid маси вантажу, початковог довжини пружини у ненавантаженому стат, коефщента жорсткостi пружини та початкових умов виникнення коливань.

Ключов1 слова: 2d-пружинний маятник, ттегральна крива, фазова траекторiя, коефщент жорсткостi пружини, комп'ютерна графта.

МЕТОД ОПРЕДЕЛЕНИЯ ОСОБЫХ ТРАЕКТОРИЙ КОЛЕБАНИЯ ГРУЗА 2D-ПРУЖИННОГО МАЯТНИКА

О.М. Семкив, к.т.н., ст.н.с., Национальный университет гражданской защиты Украины, г. Харьков

Аннотация. Разработан метод определения средствами дифференциального исчисления и компьютерной графики геометрической формы особой траектории перемещения по плоскости груза 2d-пружинного маятника в зависимости от массы груза, начальной длины пружины в ненагруженном состоянии, коэффициента жесткости пружины и начальных условий возникновения колебаний.

Ключевые слова: 2d-пружинный маятник, интегральная кривая, фазовая траектория, коэффициент жесткости пружины, компьютерная графика.

THE METHOD OF DETERMINING SPECIAL LOAD OSCILLATION TRAJECTORIES OF THE 2D-SPRING PENDULUM

O. Semkiv, Ph. D. (Eng.), Sr. Researcher, National University of Civil Protection of Ukraine

Abstract. A method for determining a geometric shape of a special trajectory of movement of the 2D-spring pendulum weight on the path by the means of differential calculus and computer graphics, depending on the load weight, the initial length of the spring in the unloaded condition, the coefficient of spring rigidity and the initial conditions of oscillation occurrence is developed.

Key words: 2D-spring pendulum, integral curve, phase trajectory, spring rigidity coefficient, computer graphics.

Пружинним маятником називають коливаль-ну систему з мехашчною пружиною (або за-мшником И функцш - деяким еластичним матерiалом), що закршлена одним кшцем

стащонарно, а на шшому кшщ знаходиться вантаж певно! маси. При цьому заведено вважати, що вантаж коливасться виключно завдяки зусиллям пружини - як у бш И стис-нення, так i у бш И розтягнення, i що конструктивно забезпечусться поперечне «незги-

нання» пружини. Звичайно пружину розгля-дають за умови нерухомостi И ош. Тобто зазначене «одновимiрне» коливання пружини здшснюеться переважно завдяки еластич-ним властивостям матерiалу, з якого и виготовлено.

Але конструкторiв може зацiкавити [1-3] i «двовимiрне» коливання тша пружини у вертикально площинi хОу навколо стацiонарно закршленого кiнця, подiбно до коливання традицшного математичного маятника (рис. 1).

Рис. 1. Схема 2а?-пружинного маятника

Далi зазначену коливальну конструкцiю бу-демо називати 2а?-пружинним маятником (тобто такою, яка здшснюе коливання у площиш). Траeкторieю коливання вантажу назвемо слщ на площинi, одержаний в результат його перемiщення. Особливою трае-кторieю вважатимемо таку лiнiю, геометрич-на форма яко! вiдрiзнятиметься деякою зако-номiрнiстю порiвняно з траекторiями перемiщення вантажу в результатi хаотичних коливань 2а?-пружинного маятника.

Доцiльнiсть дослiдження 2^-пружинних мая-тниюв виникла у зв'язку з виявленими мож-ливостями !х «нестандартних» використань як у теоретичному плаш, так i на практицi [1-5]. Дшсно, вiдомi зв'язки пружинних мая-тниюв з фазовими траекторiями динамiчних систем на поверхш тора, а також iз теорiею математичних бiльярдiв [7]. У робот [10] наведено розрахунки коливань 2^-пружин-ного маятника, але без аналiзу впливу його

параметрiв. Крiм того, диференцiальнi рiв-няння коливань 2^-пружинних маятникiв по-дiбнi до рiвнянь задач «хижак-жертва» [6], що вiдкривае напрям дослiджень. При зазна-чених коливаннях цiкавiсть викликае геомет-рична форма траектори перемщення по площинi хОу вантажу [8-10]. Ця траекторiя iлюструе розв'язок вщповщних диференща-льних рiвнянь (1), що описують коливання 2^-пружинного маятника. Одержат геомет-ричнi форми траекторiй перемщення по площиш вантажу допоможуть шюструвати розв'язки певних задач (наприклад, задачi «хижак-жертва» [6]), i !х розгляд дозволить аналiзувати (подiбно до того, як у теори коливань застосовують фiгури Лiссажу) характер розв'язюв у сумiжних за змютом задачах.

Мета i постановка завдання

Метою роботи е розробка методу визначення засобами диференщального числення та комп'ютерно! графши геометрично! форми особливо! траектори перемiщення по площи-ш вантажу 2^-пружинного маятника залежно вщ його параметрiв, де головним параметром буде значення коефщента жорсткостi k пружини. На характер коливань 2^-пружинного маятника впливатимуть такi параметри: маса вантажу, початкова довжина пружини у не-навантаженому сташ, коефiцiент жорсткосп пружини i початковi умови шщдавання коливань - такi як початковий кут вщхилення осi пружини, початкова швидюсть кута вщ-хилення 2^-пружинного маятника та швид-юсть початкового подовження тiла пружини. Для практичних впроваджень необхiдно роз-робити метод визначення набору значень цих параметрiв, якi б забезпечили особливi траектори перемiщення вантажу.

Моделювання коливань 2^-пружинного маятника

Описувати коливання 2^-пружинного маятника будемо [9] за допомогою системи дифе-ренцiальних рiвнянь

де 0(0 - кут вщхилення осi маятника вщ вер-тикалi; к - коефщент жорсткостi пружини; т - маса маятника; L(t) - функщя змiни до-вжини пружини; L0 - початкова довжина маятника. Крiм того, необхiдно задати умови: початковий кут 0(0) вщхилення тша маятника; початкову швидкiсть ^0(0) вiдхилення тiла пружини; початкову швидюсть DL(0) подовження тiла пружини, а також межi часу iнтегрування.

Фшсуемо значення всiх параметрiв, що впли-вають на характер коливань, ^м одного (по-значимо його як р), який назвемо головним i який буде змшюватися у наперед визначених межах.

Для певного значення р систему рiвнянь (1) розв'яжемо чисельним методом Рунге-Кутти. Це дозволяе у фазовому просторi одержати послiдовнiсть точок з координатами

(t), L (t), (> або |0 ((), 0 ((), (]•, як визна-

чатимуть певну iнтегральну криву розв'язку системи (1) (тут крапкою позначено похiдну за (). Ортогональною проекщею штегрально! криво! з фазового простору на фазову пло-щину буде фазова траекторш [11, 12].

На рис. 2 наведено зображення штегрально! криво! та фазово! траекторi!, якi в загальному випадку матимуть хаотичний характер.

За змши значення головного параметра р мають змшюватися як форма штегрально! криво!, так i форма фазово! траекторий Необ-хiдно визначити таке значення р, за якого фазова траекторiя набуде вигляду регулярно! криво!.

Взаемопроещюючими точками штегрально! криво! назвемо таю !! точки, ортогональш проекцi! яких на фазову площину матимуть однаковi (або близькi вiдносно точностi об-числень) координати.

Щц критичним значенням параметра р вва-жатимемо таке, коли штегральна крива набу-де форми, за якою !! точки стануть взаемо-проецiюючими. Наочно критичне значення параметра виявляеться у тому, що фазова траекторiя матиме вигляд регулярно! криво! (рис. 3).

Рис. 2. 1нтегральна крива та фазова траекто-рiя в загальному випадку

Рис. 3. 1нтегральна крива та фазова траекто-рiя у випадку критичного значення головного параметра р

За результатами вщтворення коливань 2^-пружинного маятника критичне значення кО вiдiб'еться в особливiй траектори форми коливання вантажу на площиш хОу.

У загальному випадку пропонуеться такий алгоритм пошуку особливих траекторш.

Крок 1. Обираемо поточне значення р i шляхом розв'язання системи рiвнянь (1) одержу-емо поточкове подання штегрально! криво! системи (1).

Крок 2. Будуемо поточкове зображення фа-зово1 траекторiï як проекци на фазову пло-щину iнтегральноï кривоï.

Крок 3. Здiйснюемо ощнку кiлькостi точок на фазовiй площиш, за умови ïx можливих сумщень при проецiюваннi; ця кiлькiсть точок вщповщатиме поточному значенню р.

Крок 4. Починаючи з кроку 2 виконуемо у цикл обчислення для ушх необxiдниx зна-чень головного параметра р.

Крок 5. Залежно вщ значення параметра р формуемо функцда Np(p) кшькосп точок на фазовiй площинi, як визначають зображення фазовоï траекторiï, та обчислюемо екстрему-ми цiеï функци.

Крок 6. За визначеними критичними значен-нями параметра р перевiряемо «дда» 2^-маятника шляхом побудови засобами комп'ютерно1' графiки слiду вiд коливання вантажу.

Далi узагальнений алгоритм детально проь люструемо на прикладi побудови слщу вiд коливання вантажу 2^-маятника залежно вiд головного параметра k - коефiцiента жорст-костi пружини.

Оскшьки розв'язок диференцiальниx рiвнянь одержано чисельним методом, то штегральна крива складатиметься з окремих точок. На фазовш площиш проекци цих точок можуть сшвпадати, тому фазову траекторда в момент виявлення критичного значення k0 зобража-тиме менша кшьюсть точок. На базi розроб-леного алгоритму для середовища Maple бу-ло складено програму обчислення кшькосп точок, якi утворюють зображення фазово1' траекторiï на фазовiй площинi за умови вра-хування можливих сумщень точок проекци.

Для цього було залучено команди пакета ImageTools для формування та аналiзу точ-кових зображень проекцш на фазовiй площиш 0цiнка виражаеться значеннями побу-довано1' функцiï Np(k). Критичне значення k0 визначаеться мiнiмальною кiлькiстю точок фазово1' траекторiï, що на графiку функ ци Np(k) мае проявлятися як екстремум (мшмум).

У результатi виконаних обчислень одержано графш Np(k) функци змiни кiлькостi точок-проекцiй на фазовш площиш (рис. 4). Ця фу-нкщя мае глобальний мшмум при k = 18,12, а також множину локальних мшмушв, серед яких два будуть яскраво вираженими при k = 23,0 i k = 28,84.

Наприклад, розв'яжемо систему рiвнянь (1) з параметрами m = 1; L0 = 1 i початковими умовами: початковий кут 9(0) = 0 вщхилен-ня; початкова швидюсть D9(0) = 0,5 вщхи-лення тша пружини; початкова швидюсть DL(0) = 0 подовження тша пружини; час t штегрування обрано у межах 0-8 л. Параметр k змшювався на iнтервалi 15 < k < 35, який було роздшено на S=250 частин. Значення вшх обраних параметрiв взято в умовних величинах.

Для визначення взаемопроецдаючих точок штегрально1' криво1' було побудовано посль довнiсть анiмацiйниx зображень фазових траекторш на фазовiй площиш, здшснених засобами комп'ютерно1' графiки.

Наочно процес вибору критичного значення параметра мае вигляд на фазовш площиш «як наведення на рiзкiсть» зображення, утво-реного хаотично переплетеними фазовими траекторiями (порiвняти рис. 2 i 3).

S = 250; kmin = 15; 130001,

Рис. 4. Графiк функци вщ k змiни кiлькостi точок на фазовш площиш

На рис. 5-7 наведено вщповщш цим випад-кам штегральш кривi та 1'х фазовi траекторiï. На кожному з рисунюв 5-7 зображено вщпо-вщно до знайденого значення k фазову траекторда: