МАТЕМАТИЧЕСКОЕ И ПРОГРАММНОЕ ОБЕСПЕЧЕНИЕ ДЛЯ ИССЛЕДОВАНИЯ ПОПЕРЕЧНЫХ КОЛЕБАНИЙ ГРЕБНЫХ ВАЛОВ СУДОВ СТАТЬЯ
Журнал: ВЕСТНИК АСТРАХАНСКОГО ГОСУДАРСТВЕННОГО ТЕХНИЧЕСКОГО УНИВЕРСИТЕТА. СЕРИЯ: УПРАВЛЕНИЕ, ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА И ИНФОРМАТИКА
Том 2014 № 3
, 2014
Рубрики: МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
Авторы:
1.
Астраханский государственный технический университет
2.
Астраханский государственный технический университет
Россия
Россия
Тип:
Статья
Страницы:
с 102
по 112
Статус:
Опубликован
Получено:
05.07.2014
Одобрено:
25.07.2014
Опубликовано:
25.07.2014
Классификаторы:
ВАК 05.12.2013 Системы, сети и устройства телекоммуникаций
ВАК 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
ВАК 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
ВАК 05.13.10 Управление в социальных и экономических системах
ВАК 05.13.18 Математическое моделирование, численные методы и комплексы программ
ВАК 05.13.19 Методы и системы защиты информации, информационная безопасность
УДК 629.5.035-233.1:004.45
ГРНТИ 20.01 Общие вопросы информатики
ГРНТИ 28.01 Общие вопросы кибернетики
ГРНТИ 49.01 Общие вопросы связи
ГРНТИ 50.01 Общие вопросы автоматики и вычислительной техники
ГРНТИ 82.01 Общие вопросы организации и управления
ВАК 05.13.01 Системный анализ, управление и обработка информации (по отраслям)
ВАК 05.13.06 Автоматизация и управление технологическими процессами и производствами (по отраслям)
ВАК 05.13.10 Управление в социальных и экономических системах
ВАК 05.13.18 Математическое моделирование, численные методы и комплексы программ
ВАК 05.13.19 Методы и системы защиты информации, информационная безопасность
УДК 629.5.035-233.1:004.45
ГРНТИ 20.01 Общие вопросы информатики
ГРНТИ 28.01 Общие вопросы кибернетики
ГРНТИ 49.01 Общие вопросы связи
ГРНТИ 50.01 Общие вопросы автоматики и вычислительной техники
ГРНТИ 82.01 Общие вопросы организации и управления
Язык материала:
русский
Ключевые слова:
валопровод, поперечные колебания, собственная частота, форма колебаний, численный эксперимент
Аннотация (русский):
Отрабатывается расчетная схема для исследования поперечных колебаний гребного вала судов. Рассмотрены три варианта опирания и нагружения. Полученные решения позволяют исследовать как собственную частоту, так и форму колебаний. С целью автоматизации численных экспериментов разработаны три программы для ЭВМ, реализованные на алгоритмическом языке Maple. Проведены численные эксперименты с использованием разработанных программ для конкретных судов. Анализ результатов численных экспериментов позволил установить, что податливость подшипников гребного вала уменьшает его собственную частоту; динамические нагрузки от работающего гребного винта могут существенно изменить параметры центровки валопровода и привести к динамической неустойчивости валопровода. Полученные решения могут быть использованы для определения собственной частоты и форм колебаний конструкций, расчётные схемы которых сводятся к подобным расчетным схемам.
Отрабатывается расчетная схема для исследования поперечных колебаний гребного вала судов. Рассмотрены три варианта опирания и нагружения. Полученные решения позволяют исследовать как собственную частоту, так и форму колебаний. С целью автоматизации численных экспериментов разработаны три программы для ЭВМ, реализованные на алгоритмическом языке Maple. Проведены численные эксперименты с использованием разработанных программ для конкретных судов. Анализ результатов численных экспериментов позволил установить, что податливость подшипников гребного вала уменьшает его собственную частоту; динамические нагрузки от работающего гребного винта могут существенно изменить параметры центровки валопровода и привести к динамической неустойчивости валопровода. Полученные решения могут быть использованы для определения собственной частоты и форм колебаний конструкций, расчётные схемы которых сводятся к подобным расчетным схемам.
Ключевые слова:
валопровод, поперечные колебания, собственная частота, форма колебаний, численный эксперимент
валопровод, поперечные колебания, собственная частота, форма колебаний, численный эксперимент
Введение Судовой валопровод - часть двигательно-движительной установки судна (или корабля), предназначенная для передачи крутящего момента от двигателя к гребному винту и восприятия осевого усилия, возникающего при ее работе [1]. Система валопровода любого судна является одной из основных систем, обеспечивающих его жизнедеятельность. От качества работы валопровода зависит не только ход судна, но и условия обитания на нем. Любой «отказ» в надежной и устойчивой работе валопровода может привести к серьезным последствиям вплоть до гибели судна. Актуальность решаемой задачи Проблема обеспечения надежности работы валопровода и дейдвудных подшипников при эксплуатации судов, особенно крупнотоннажных, является одной из наиболее важных и сложных, не имеющих до сих пор удовлетворительного решения. Практически срок эксплуатации крупных судов без ремонта стал зависеть почти исключительно от технического состояния дейдвудного узла. Именно поэтому проблеме расчетов и проектирования валопроводов посвящено большое число работ как исследовательского, так и технологического характера. Однако проведенные мероприятия не устранили многие причины, снижающие устойчивость работы валопровода, а используемые расчетные схемы не позволяют изучить целый ряд явлений, имеющих место при работе валопровода. Проблема состоит в том, что валопровод при движении судна подвергается воздействию многочисленных нагрузок, стационарных и переменных. Многие из нагрузок еще недостаточно изучены. Кроме того, сама конструкция валопровода является сложной. Попытка учесть как можно точнее нагрузки, действующие на валопровод, и особенности его строения приводит к громоздким и сложным вычислениям, что затрудняет исследования. Использование средств вычислительной техники при расчете валопровода позволяет отчасти решить стоящие при его проектировании задачи. Однако внедрение ЭВМ в практику расчета валопровода сдерживается нехваткой или даже отсутствием соответствующих алгоритмов решения задач, возникающих в процессе эксплуатации валопровода, и программ для реализации разрабатываемых алгоритмов. Таким образом, разработка новых и совершенствование существующих расчетных схем валопроводов судов, реализация новых алгоритмов решения возникающих в практике задач и создание компьютерных программ актуальны в настоящее время и для практики, и с точки зрения программирования. Разработка алгоритма решаемой задачи Выбор расчётной схемы. Одной из важнейших динамических задач для любой упругой системы, в том числе и валопровода, является исследование ее собственной частоты и форм колебаний. Существуют различные методы расчета собственной частоты валопровода [2]. Однако во всех методах, с целью упрощения и сокращения объема вычислений, вводятся различные допущения и упрощения, что приводит к погрешности получаемых результатов. Использование ЭВМ позволяет оценить погрешность, вносимую различными допущениями, и построить более точные расчетные схемы. Нами разработаны алгоритмы решения некоторых задач, связанных с исследованием динамического воздействия на валопровод. С целью отработки расчетной схемы для исследования динамического воздействия на валопровод рассмотрим возможные ее варианты, отличающиеся нагружением и характером опирания. Из теории колебаний известно, что в сложных механических системах может «резонировать» любая ее часть, собственная частота колебаний которой совпадает с частотой возмущающей нагрузки. Именно поэтому при исследовании колебаний таких систем можно независимо рассматривать любую ее часть, заменяя влияние отброшенных элементов упругими связями [3]. В валопроводе крупнотоннажных судов наименьшую собственную частоту имеет кормовой участок гребного вала [4, 5], поэтому именно эта часть гребного вала интересует нас в первую очередь. В качестве расчетной схемы кормового участка гребного вала примем однопролетную балку постоянного сечения с консолью, нагруженную собственным весом q, весом гребного винта F и нагрузкой, изменяющейся с лопастной частотой (рис. 1). Рис. 1. Расчетные схемы кормового участка гребного вала: а - опоры защемляющая и шарнирная; б - опоры защемляющая и упругая; в - обе опоры шарнирные: l, L - длины отдельных участков гребного вала; q - погонный вес; m - погонная масса вала; М - масса винта; EI - изгибная жесткость сечения вала; Fasin ωt и Masin ωt - периодические нагрузки, действующие на вал; ω - круговая частота изменения нагрузок; с - жесткость упругой опоры; M0, Q0, R - реакции связей С целью оценки влияния опирания балки на ее собственную частоту в настоящей работе рассмотрены три варианта граничных условий: а - защемляющая и шарнирная опоры; б - защемляющая и упругая опоры; в - обе опоры шарнирные. Известно, что колебания системы с постоянными параметрами происходят относительно ее статического равновесия, поэтому при исследовании ее колебаний обычно статическая деформация системы под действием постоянных нагрузок не рассматривается, т. к. она не влияет на процесс колебаний. Математическое описание алгоритма решения задачи. Дифференциальное уравнение поперечных колебаний стержня имеет вид [6]: - : . (1) Решение уравнения (1) ищем по методу Фурье: , (2) где - уравнение формы изгиба балки в процессе колебаний. После подстановки выражения (2) в (1) получаем следующее дифференциальное уравнение для форм колебаний балки: , (3) где . (4) Решение уравнения (3) находим по методу начальных параметров Коши: - на участке : ; - на участке : , где - так называемые начальные параметры, т. е. соответственно прогиб, угол поворота сечения в начале координат (z = 0) и изгибающий момент и поперечная сила в этом же сечении; - система фундаментальных функций с единичной матрицей уравнения (3) аргументов и : ; ; ; . Решение рассматриваемых задач Постоянные интегрирования и неизвестную реакцию находим из граничных условий (рис. 1): - схема а: - схема б: - схема в: В результате, после соответствующих преобразований (см. [7-9]), находим: - схема а: (5) где - схема б: (6) где - схема в: (7) где Так как прогибы сечений y, углы поворота сечений , изгибающие моменты и поперечные силы Q в сечениях стержней связаны известными соотношениями: , то после определения постоянных интегрирования и реакции R можно исследовать форму колебаний стержня и усилия в сечениях в процессе колебания стержня. Условия для вычисления собственных частот колебаний получим, приравняв к нулю определители систем уравнений (5)-(7), т. е. условие для определения собственных частот колебаний имеет вид . (8) При заданных параметрах системы из условия (8) определяем коэффициент , а затем собственную частоту по формуле (см. выражение (4)): . (9) Алгоритм программы для ЭВМ для исследования поперечных колебаний стержня Численные эксперименты по исследованию поперечных колебаний гребного вала судов связаны с весьма громоздкими вычислениями. С целью автоматизации расчетов были разработаны программы для ЭВМ, которые объединены в программный комплекс для расчета поперечных колебаний гребных валов судов «Валопровод 13» [10]. Комплекс программ состоит из трех программ, каждая из которых производит вычисления по одной из расчетных схем, изображенных на рис. 1. В каждой программе рассматриваются два варианта задачи: 1) задана периодически изменяющаяся нагрузка; нужно исследовать форму колебания и динамическое нагружение стрежня; 2) заданы параметры стержня; нужно вычислить собственную частоту. Типовой вариант блок-схемы программ для ЭВМ представлен на рис. 2. Рис. 2. Блок-схема программы для ЭВМ для исследования поперечных колебаний гребного вала Все программы, входящие в пакет «Валопровод 13», реализованы в системе компьютерной математики Maple [11]. Выбор системы Maple обусловлен широкими возможностями этой программы по решению математически сложных задач, а также нашими личными предпочтениями. Система Maple быстро и эффективно выполняет символьные и численные вычисления (например, дифференцирование функций, численное и аналитическое интегрирование, вычисление пределов функций), а также решение систем линейных, нелинейных, дифференциальных уравнений, решение трансцендентных уравнений и множество других операций. При этом не только заметно сокращается рутинная часть вычислений, но и, благодаря развитым графическим возможностям системы, предоставляется возможность получить наглядный иллюстративный материал к условию и решению задачи, оформленный в виде рисунков, графиков, диаграмм, 2D- и 3D-анимированных моделей. Рассмотрим краткое описание, фрагменты листингов и примеры решения типовой программы. Краткое описание программы (схема б). Рассматривается кормовой участок гребного вала, содержащий консоль гребного вала, винт и дейдвудный пролет (см. рис. 1, б). По заданным параметрам стержня, к которым относятся: L - общая длина вала, м; l - длина участка вала, м; g - ускорение свободного падения, м/с2; EI - изгибная жесткость сечения вала, кН ∙ м2; F - вес винта, кН; q - погонный вес вала, кН/м - определяется важнейшая динамическая характеристика - собственная частота поперечных колебаний. Решение задачи сводится к трансцендентному алгебраическому уравнению: (10) Для определения корней уравнения (10) использовано приближенное численное решение при помощи Maple-команды fsolve. Поскольку уравнение (10) является трансцендентным уравнением, то команда fsolve, по умолчанию, выводит только один корень. При решении трансцендентных уравнений желательно указывать предполагаемое положение корней. С этой целью путем построения графика функции , для которой разыскивается решение, определен нужный диапазон (или начальная точка). Для определения всех корней уравнения (10) предварительно получено графическое решение - построен график функции , в котором абсциссы точек пересечения графика с осью Ох являются искомыми корнями. Таким образом, разработанная программа позволяет определить корни трансцендентного уравнения вида и затем по уравнению (9) вычислить собственную частоту поперечных колебаний гребных валов судов. Влияние расчетной схемы на величину собственной частоты поперечных колебаний рассмотрим на примере гребного вала с параметрами [1, с. 69]: F = 73,85 кН; q = 11,05 кH/м; Жесткость упругой опоры принимаем в пределах кН/м2 [12]. Пример работы программы для определения собственной частоты поперечных колебаний (расчет произведён для схемы б) 1. Очищаем внутреннюю память Maple, подключаем пакеты plots, plottools, linalg: restart; # Одновременная очистка всех переменных with(plots): with(plottools): with(linalg): # подключение специальных библиотек подпрограмм: plots - библиотека, позволяющая строить двухмерные графики plottools - библиотека для построения графических примитивов linalg - библиотека линейной алгебры 2. Вводим числовые исходные данные: L:= _: l:= _: g:= _: EI:=_: F:= _: q:= _: # общая длина вала, м # длина участка вала, м # ускорение свободного падения, м/с2 # изгибная жесткость сечения вала, кН∙м2 # вес винта, кН # погонный вес вала, кН/м 3. Производим вычисления: M := F/g; m := q/g; alpha := betta/L; betta1 := betta*l/L; betta2 := betta*(L-l)/L; lambda:=C*L^3/EI; # масса винта, кг # погонная масса вала, кг/м # коэффициенты 4. Записываем систему фундаментальных функций: K1b := (cosh(betta)+cos(betta))/2; K1b2 := (cosh(betta2)+cos(betta2))/2; K2b := (sinh(betta)+sin(betta))/2; K2b2 := (sinh(betta2)+sin(betta2))/2; K3b := (cosh(betta)-cos(betta))/2; K3b1 := (cosh(betta1)-cos(betta1))/2; K4b := (sinh(betta)-sin(betta))/2; K4b1 := (sinh(betta1)-sin(betta1))/2; K4b2 := (sinh(betta2)-sin(betta2))/2; # система фундаментальных функций 5. Записываем определитель системы уравнений - трансцендентное уравнение вида : DD := proc (betta) options operator, arrow; (K1b-lambda*K3b1*K2b2/betta^3)*(K1b-lambda*K4b1*K1b2/betta^3+M*betta*(K4b-lambda*K4b1*K4b2/betta^3)/(m*L))-(K2b-lambda*K4b1*K2b2/betta^3)*(K4b-lambda*K3b1*K1b2/betta^3+M*betta*(K3b-lambda*K3b1*K4b2/betta^3)/(m*L)) end proc; 6. Путем численного решения уравнения находим его корни (значения коэффициента при в диапазоне значений от a до b): fsolve(DD(betta), betta =a .. b); betta :=_. 7. Удостоверяемся, что найденное значение является корнем уравнения (при подстановке найденного значения должно выполняться условие ): evalf(Delta); 0. 8. Для найденного значения определяем собственную частоту поперечных колебаний гребных валов судов, 1/с: omega:=(betta/L)^2*sqrt(EI/m); evalf(omega). Результаты вычислений (для схем а, б и в) приведены в таблице. Результаты вычислений собственной частоты Схема (рис. 1) Параметры а б в 2,49 2,29 2,38 2,35 140 118 128 126 Разработанные нами компьютерные программы могут быть использованы для определения собственной частоты и форм колебаний конструкций, расчётные схемы которых сводятся к подобным расчетным схемам. Заключение Таким образом, нами предлагается алгоритм исследования собственной частоты и форм поперечных колебаний гребного вала валопроводов судов. Рассмотрены три варианта расчётных схем. С целью автоматизации расчетов при исследовании поперечных колебаний валопроводов судов разработаны алгоритмы решения задач на ЭВМ, которые реализованы в виде программ для ЭВМ на алгоритмическом языке Maple. Проведены численные эксперименты с использованием разработанных программ для конкретных судов. Анализ результатов численных экспериментов позволил установить следующее: - податливость подшипников гребного вала уменьшает его собственную частоту до по сравнению с частотой, рассчитанной без учета податливости подшипников; - динамические нагрузки, например от работающего гребного винта, могут существенно изменить параметры центровки валопровода, рассчитанные только статическим нагружением, что может привести к динамической неустойчивости валопровода; - полученные результаты соответствуют жесткости опирания балки: наибольшую жесткость имеет схема а - собственная частота поперечных колебаний наибольшая; чем больше жесткость упругой опоры, тем выше собственная частота - схема б. Разработанные компьютерные программы могут быть использованы для определения собственной частоты и форм колебаний конструкций, расчётные схемы которых сводятся к подобным расчетным схемам.
1. Подсевалов Б. В. Словарь стандартизованной терминологии в судостроении / Б. В. Подсевалов, А. П. Фомин. Л.: Судостроение, 1990. 240 с.
2. Комаров В. В. Валопроводы рыбопромысловых судов. Ч. 1. Конструкция, эксплуатация и общие вопросы проектирования / В. В. Комаров, А. С. Курылёв. Астрахань: Изд-во АГТУ, 1997. 166 с.
3. Дондошанский В. К. Расчет колебаний упругих систем на ЭВМ / В. К. Дондошанский. М.: Машиностроение, 1965. 368 с.
4. Николаев В. А. Конструирование и расчет судовых валопроводов / В. А. Николаев. Л.: Гос. союз. изд-во судостроит. пром-ти, 1958. 358 с.
5. Денисова Л. М. Сравнительный анализ собственных частот различных участков гребного вала судов / Л. М. Денисова, А. И. Миронов // Вестн. Астрахан. гос. техн. ун-та. Сер.: Морская техника и технология. 2011. № 2. С. 28-31.
6. Пановко Я. Г. Введение в теорию механических колебаний / Я. Г. Пановко. М.: Наука, 1971. 239 с.
7. Миронов А. И. К исследованию поперечных колебаний гребных валов. Часть 1 / А. И. Миронов // Вестн. Астрахан. гос. техн. ун-та. Сер.: Морская техника и технология. 2013. № 2. C. 125-130.
8. Миронов А. И. К исследованию поперечных колебаний гребных валов. Часть 2. Влияние упругой податливости подшипников на процесс колебаний вала / А. И. Миронов // Вестн. Астрахан. гос. техн. ун-та. Сер.: Морская техника и технология. 2014. № 1. С. 77-82.
9. Миронов А. И. Влияние колебаний вала на параметры «центровки» валопровода / А. И. Миронов, Е. В. Пономарёва // Вестн. Астрахан. гос. техн. ун-та. Сер.: Морская техника и технология. 2014. № 2. С. 86-94.
10. Миронов А. И. Программный комплекс для расчета поперечных колебаний гребных валов судов. Свид-во об офиц. регистрации 2014611733 / Миронов А. И., Пономарёва Е. В.; заявка № 2013660700; зарегистр. в реестре программ для ЭВМ 10.02.2014.
11. Матросов А. В. Maple 6. Решение задач высшей математики и механики / А. В. Матросов. СПб.: БХВ-Петербург, 2001. 528 с.
12. Денисова Л. М. Исследование поперечных колебаний гребных валов / Л. М. Денисова, А. И. Миронов // Вестн. Астрахан. гос. техн. ун-та. Сер.: Морская техника и технология. 2005. № 2 (25). C. 98-103.
