WWW.BOOK.LIB-I.RU
БЕСПЛАТНАЯ  ИНТЕРНЕТ  БИБЛИОТЕКА - Электронные ресурсы
 


«УДК 539.3 В. Д. Будак1, А. Я. Григоренко2, М. Ю. Борисенко1, Е. В. Бойчук1 ОПРЕДЕЛЕНИЕ СОБСТВЕННЫХ ЧАСТОТ ЭЛЛИПТИЧЕСКОЙ ОБОЛОЧКИ ...»

УДК 539.3

В. Д. Будак1, А. Я. Григоренко2, М. Ю. Борисенко1, Е. В. Бойчук1

ОПРЕДЕЛЕНИЕ СОБСТВЕННЫХ ЧАСТОТ ЭЛЛИПТИЧЕСКОЙ ОБОЛОЧКИ

ПОСТОЯННОЙ ТОЛЩИНЫ МЕТОДОМ КОНЕЧНЫХ ЭЛЕМЕНТОВ

Исследуются свободные колебания тонкой изотропной эллиптической оболочки постоянной толщины на основе метода конечных элементов, реализованного на высокопроизводительном программном обеспечении для инженерного анализа – FEMAP. Приводится сравнительный анализ числовых результатов с полученными экспериментально.

Постановка проблемы. Во многих областях техники широкое применение находят некруговые цилиндрические оболочки. Для эксплуатации таких конструкций необходимо знать условия прочности и надежности, в частности, иметь представление о частотах и формах свободных колебаний, поскольку необходимо избегать резонансных частот, которые ведут к разрушению конструкции в целом. В связи с этим особое значение приобретают надежные методы численного и экспериментального определения динамических характеристик таких конструкций.

В данной работе исследуются свободные колебания консольно закрепленной эллиптической цилиндрической оболочки постоянной толщины. Задачам о свободных колебаниях некруговых цилиндрических оболочек посвящено небольшое количество публикаций. В работах [2, 3, 9] рассматривались свободные колебания элементов оболочечных конструкций, в частности некруговых цилиндрических оболочек, рассчитанные численными методами.

Собственные частоты и формы колебаний также определяются экспериментальными методами. Одним из таких методов является метод голографической интерферометрии, который позволяет в реальном времени наблюдать и фиксировать с большой точностью резонансные частоты и формы колебаний видимого участка исследуемой поверхности [1, 7]. Анализ отмеченных публикаций дает возможность сделать вывод, что не существует единой точки зрения относительно применения того или другого подхода к решению указанного класса задач.

Чаще всего динамику оболочечных конструкций моделируют в рамках классической теории Киргофа – Лява с использованием численных методов [4–6]. В случае переменной толщины или произвольной формы серединной поверхности перемещения точек таких оболочек описывается системой дифференциальных уравнений в частных производных с переменными коэффициентами. В этом случае не удается разделить переменные с использованием рядов Фурье в одном из координатных направлений, и решение задач о свободных колебаниях оболочек переменной толщины или произвольной формы срединной поверхности приводит к значительным трудностям вычислительного характера.

В настоящее время широкое применение для решения задач механики получило использование Computer-Aided Engineering (CAE) систем, которые могут рассчитывать конструкцию любой формы благодаря использованию метода конечных элементов. Кроме того, CAE системы дают пользователю возможность оценить поведение компьютерной модели изделия в реальных условиях эксплуатации, проверить дееспособность конструкции без значительных вложений времени и средств. Одной из таких систем является пре- и постпроцессор для выполнения инженерного анализа методом конечных элементов – FEMAP с решателем NX NASTRAN [8].

Целью настоящей работы является исследование свободных колебаний упругой изотропной цилиндрической оболочки эллиптического поперечного сечения на основе программы FEMAP и сравнение их с экспериментально полученными значениями.





ISSN 0130–9420. Мат. методи та фіз.-мех. поля. 2014. – 57, № 1. – С. 145-152. 145 Исходные соотношения. Рассматривается задача о собственных частотах и формах колебаний. Уравнение динамики при отсутствии внешних влияний и переменных кинематических граничных условий имеет вид && & [M ]{q (t )} + []{q (t )} + [K ]{q (t )} = 0, (1) где [M ] – матрица масс конструкции; [] – матрица коэффициентов сил вязкого демпфирования; [K ] – матрица жесткости; {q} – вектор узловых смещений; точкой и двумя точками над q в уравнении обозначены первая и вторая ее производные по времени соответственно.

В случае отсутствия демпфирования уравнение (1) примет вид && [M ]{q (t )} + [K ]{q (t )} = 0. (2) Решением матричного уравнения (2) является {q (t )} = {A} cos (wt + b), (3) где {A} – вектор амплитудных значений узловых перемещений; w = 2pf – круговая частота; b – фаза колебаний.

После прямой подстановки (3) в (2) и сокращения на cos (wt + b) получим систему алгебраических уравнений:

( - w2 [M ] + [K ]){A} = 0. (4) В этой системе ненулевые значения компонент {A} возможны лишь при условии, что det [[K ] - w2 [M ]] = 0. (5) Если квадратичные матрицы [M ] и [K ] положительно определены, что обычно так и есть для задачи линейной упругости, то уравнение (5) имеет N неотрицательных решений – собственных частот wk, причем возможны парные значения (здесь N – количество неизвестных в системе (4)).

При выполнении условия (5) одно из уравнений (4) является следствием остальных. Поэтому каждому значению wk соответствует определенное соотношение между амплитудами Аki. Другими словами, все амплитуды вектора могут быть выражены через одну из них. Соотношения между амплитудами Аki определяют k -ю собственную форму колебаний.

Из уравнений (4) следует, что все степени свободы в процессе колебаний с собственной частотой wk совершают синхронное движение. Таким образом, конфигурация конструкции не меняет своей базовой формы, а меняются только амплитуды.

Имея N значений собственных частот wk, решение системы (2) можно искать в виде линейной комбинации из N выражений (3):

N {k } cos (wk t + bk ).

{q (t )} = (6) k =1 Таким образом, изменение формы линейной упругой конструкции со временем есть линейной комбинацией всех ее собственных форм.

Поскольку значения компонент собственных векторов {k } могут быть найдены не однозначно, а с точностью до постоянного множителя, то обычно их нормируют по правилу {Ak }• [M ]{Ak } = 1. (7) Доказано, что собственные векторы {k } ортогональны относительно матриц [M ] и [K ], то есть {Ak }• [M ]{Am } = 0, {Ak }• [K ]{Am } = 0, k m. (8) Еще можно отметить, что обычно ищутся не все корни уравнения (5), а несколько (обозначим через Na ) наименьших значений, поскольку только при низших собственных частотах амплитуды колебаний имеют относительно большие значения. Иначе говоря, несколько первых собственных частот и форм колебаний достаточно для получения удовлетворительного приближения решения (6). Для этого разработано несколько алгоритмов, лучшие из них используются в расчетном модуле NX Nastran [8].

Методика решения. С помощью FEMAP была построена геометрия оболочки в виде цилиндрической поверхности эллиптического сечения с размерами: высота h = 120 мм, большая полуось a = 51.8 мм, малая полуось b = 37.295 мм.

Параметры моделирования выбраны в соответствии с формой и размерами оболочки, которая исследовалась на собственные частоты и формы колебаний методом голографической интерферометрии (рис. 1). Оболочка 1 выточена вместе с массивным телом 2 для реализации жёсткого закрепления по контуру при экспериментальных исследованиях. Поскольку метод голографической интерферометрии показал нулевые перемещения балясины, то при расчете с помощью FEMAP рационально моделировать оболочку консольно защемленной по одному из контуров.

Материалом оболочки является сталь 40Х с параРис.1 метрами: модуль Юнга E = 214 ГПа, коэффициент Пуассона n = 0.26, плотность r = 7820 кг/м. Разбивка производилась линейными четырёхугольными plate-элементами размером 1мм 1мм постоянной толщины d = 2 мм и содержала 33396 узлов и 33120 элементов (рис. 2).

Задавалось консольное закрепление по нижнему контуру оболочки.

–  –  –

Мода 6: f = 3860 Гц Мода 7: f = 3863 Гц Рис. 5 Выводы. Моделирование оболочки в виде консольно защемленной по одному из контуров цилиндрической поверхности эллиптического сечения достаточно надежно отображает условия проводимого эксперимента.

Экспериментально полученные частоты и частоты, рассчитанные с помощью NASTRAN, отличаются не более чем на 13 %, для некоторых значений различие составляет меньше 1 %.

Численным моделированием можно получить полный спектр форм колебаний в отличие от экспериментального метода, где результаты зависят от условия возбуждения, точности измерительной аппаратуры и чистоты проведения эксперимента.

Достаточно высокое совпадение экспериментальных и теоретических частот указывает на высокую достоверность обеих методик исследования.

1. Григоренко А. Я., Золотой Ю. Г., Пригода А. П., Жук И. Ю., Хоришко В. В., Овчаренко А. В. Экспериментальное исследование собственных колебаний толстостенной цилиндрической оболочки методом голографической интерферометрии // Мат. методи та фіз.-мех. поля. – 2012. – 55, № 3. – С. 93–98.

То же: Grigorenko A. Ya., Zolotoi Yu. G., Prigoda A. P., Zhuk I. Yu., Khorishko V. V., Ovcharenko A. V. Experimental investigation of natural vibrations of a thick-walled cylindrical shell by the method of holographic interferometry // J.

Math. Sci. – 2013. – 194, No. 3. – P. 239–244.

2. Григоренко А. Я., Пузырев С. В., Волчек Е. А. Исследование свободных колебаний некруговых цилиндрических оболочек с помощью метода сплайн-коллокации // Мат. методи та фіз.-мех. поля. – 2011. – 54, № 3. – С. 60–69.

То же: Grigorenko A. Ya., Puzyrev S. V., Volchek E. A. Investigation of free vibrations of noncircular cylindrical shells by the spline-collocation method // J.

Math. Sci. – 2012. – 185, No. 6. – P. 824–836.

3. Григоренко Я. М., Беспалова Е. И., Китайгородский А. Б., Шинкарь А. И. Свободные колебания элементов оболочечных конструкций. – Киев: Наук. думка, 1986. – 171 с.

4. Григоренко Я. М., Василенко А. Т. Теория оболочек переменной жесткости. – Киев: Наук. думка, 1981. – 544 с. – (Методы расчета оболочек: В 5 т. – Т. 4.)

5. Григоренко Я. М., Мукоед А. П. Решение нелинейных задач теории оболочек на ЭВМ. – Киев: Вища шк., 1983. – 286 с.

6. Григоренко Я. М., Панкратова Н. Д. Обчислювальні методи в задачах прикладної математики. – Київ: Либідь, 1995. – 272 с.

7. Макаева Р. Х., Царева А. М., Каримов А. Х. Определение вибрационных характеристик дисков постоянной толщины экспериментально-расчетным методом с применением голографической интерферометрии // Голография в России и за рубежом. Наука и практика: Сб. тр. IV-й Междунар. науч.-практич. конф. – Москва, 2007. – С. 116–118.

8. Рудаков К. Н. FEMAP 10.2.0. Геометрическое и конечно-элементное моделирование конструкций. – Киев: Нац. техн. ун-т Украины «КПИ», 2011. – 317 с.

9. Sabih I. Hayek, Boisvert Jeffrey E. Vibration of elliptic cylindrical shells: higher order shell theory // J. Acoust. Soc. Am. – 2010. – 128, No. 3. – P. 1063–1072.

ВИЗНАЧЕННЯ ВЛАСНИХ ЧАСТОТ ЕЛІПТИЧНОЇ ОБОЛОНКИ ПОСТІЙНОЇ ТОВЩИНИ

МЕТОДОМ СКІНЧЕННИХ ЕЛЕМЕНТІВ

Досліджуються вільні коливання тонкої ізотропної еліптичної оболонки постійної товщини на основі методу скінченних елементів, реалізованого на високопродуктивному програмному забезпеченні для інженерного аналізу – FEMAP. Проводиться порівняльний аналіз чисельних результатів з отриманими експериментально.

DETERMINATION OF NATURAL FREQUENCIES OF AN ELLIPTIC SHELL OF CONSTANT

THICKNESS BY THE FINITE ELEMENT METHOD






Похожие работы:

«Сообщение о существенном факте "Об отдельных решениях, принятых советом директоров эмитента 1. Общие сведения 1.1. Полное фирменное наименование Открытое акционерное общество эмитента "Центр по перево...»

«Інструкція з експлуатації Пральна машина Зміст Установка, 2-3 Розпакування та установлення по рівню Українська, 1 Підключення до електричної мережі та водопостачання IWC 7085 Перший цикл прання Технічні дані Опис пральної машини та початок прання, 4-5 Пульт управління Індикаторні лампи Початок прання Цикли прання, 6 Таблиця циклів прання На...»

«1 ОБЩИЕ ПОЛОЖЕНИЯ Примерная основная образовательная программа начального общего образования разработана в соответствии с требованиями федерального государственного образовательного стандарта начального общего образования (далее — Стандарт) к структуре основной образовательной программы, определяет цель, задачи, планируемые результаты...»

«КУРБАНОВ К.К. ПРОБЛЕМЫ ОРГАНИЗАЦИИ ЭФФЕКТИВНОГО ФУНКЦИОНИРОВАНИЯ АПК РЕСПУБЛИКИ ДАГЕСТАН 1 Статья посвящена проблемам организации эффективного функционирования АПК Республики Дагестан. Рассматрив...»

«Первые шаги регистрация на программы АССА Издано Представительством АССА в Украине ул. Красноармейская, 9/2, оф. 77, Киев, 01004, Украина тел: (+38 044) 498 34 50; факс: (+38 044) 498 34 51 www.accaglobal.com...»

«Открытое акционерное общество Иркутское авиационное производственное объединение ИНН 3807002509 1. Вид ценных бумаг (облигации), идентификационные признаки выпуска Облигаций. Процентные документарные неконвертируемые облигации на предъявителя серии А02 с обязательным централизованным хранением Вид ценных бумаг: об...»

«Газета №2 группы "Улыбка" декабрь, январь, февраль 2016 -2017 г. "Сорока белобока"Поздравляем с днм рожденья: В декабре Ларина Илью В январе: Крюкову Милану, Юрчак Арину, Тарабрина Володю! В феврале: Хомутова Дмитрия! Пословицы о зиме. Зима – не лето, в шубу одета. Снегирь прилетит о зиме известит. Год кончается, а зима начинаетс...»

«НАЦИОНАЛЬНЫЙ ЦЕНТР ОБУЧЕНИЯ НАВЫКАМ ОКАЗАНИЯ ПЕРВОЙ ПОМОЩИ "ШКОЛА БУБНОВА" В.Г. Бубнов АВТОМОБИЛЬНАЯ АПТЕЧКА Оптимальная комплектация аптечки водителя для эффективного оказания первой помощи в дорожно-транспортных происшествиях Москва ПОЯСНИТЕЛЬНАЯ ЗАПИСКА Прежде чем приступить...»

«УДК 62.531.7 Ю. М. Зубарев, д.т.н., проф., С. В. Косаревский, к.т.н., доц., СанктПетербургский институт машиностроения E-mail: texmash@zavod-vtuz.ru Координатные измерения на производстве Статья посвящена применении информационных технологий на производстве, в частность при проведении процесаа измерения. Анализи...»








 
2017 www.book.lib-i.ru - «Бесплатная электронная библиотека - электронные ресурсы»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.