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

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

Рубрика Программирование, компьютеры и кибернетика
Вид статья
Язык русский
Дата добавления 30.05.2017
Размер файла 141,7 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

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

Т.В. Бондаренко

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

Вращающееся на опорах кольцо представляет собой сложную систему с распределенными параметрами, поведение которой описывается дифференциальными уравнениями в частных производных. При вращении кольца в нем возникают колебания, которые и необходимо определить. Сложность вычисления повышает наличие нулевых ограничений на перемещение кольца в точках опор. [2] Задачу определения собственных частот колебаний конструкции можно решить путем представления решения в виде ряда Фурье с заданием нулевых перемещений конструкции в точках опор. При таком подходе достигается высокая точность результатов, но усложняется процесс решения вследствие возникновения дополнительных условий связи на решение, которые требуют учета в системе дифференциальных уравнений алгебраических уравнений и неопределенных множителей Лагранжа. дифференциальный уравнение частота колебание

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

Рассмотрим конструкцию, представляющую собой кольцо, базирующееся на двух опорах (рис.1.).

Рис. 1 Кольцо, базирующееся на двух опорах

Введем следующие обозначения:

- u, v - радиальное и тангенциальное перемещение точек кольца;

- и - система координат связанная с кольцом, определяющая угол поворота: начало отсчета угла находится в первой точке опор, отсчет угла поворота ведется против часовой стрелки;

- - угловые координаты точек опор

- - угол между опорами.

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

Точки опор разбивают координатную прямую на два участка, заключенных между ними:

- первый участок, заключенный между точками ,

- второй участок, заключенный между точками .

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

В качестве базисных функций используем линейные функции вида , при этом неизвестные коэффициенты k0, k1, которые вычисляются исходя из равенства нулю перемещений кольца в точках опор.

Зададим перемещения точек кольца в виде ряда Фурье, используя синус угла поворота точки кольца:

(1)

где:

- A - коэффициент пропорциональности, задающий амплитуду колебаний кольца пропорционально дуге участка кольца, заключенной между точками опор:

для дуги, содержащей первый участок кольца;

для дуги, содержащей второй участок кольца.

- bui, bvi - неизвестные обобщенные координаты;

- и - угол, задающий положение точки кольца.

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

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

1. для первого участка функция , тогда система уравнений имеет вид:

2. для второго участка функция , тогда система уравнений имеет вид:

Решив полученные системы уравнений, получили следующие значения коэффициентов базисных функций:

.

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

На первом участке кольца:

.

На втором участке кольца:

.

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

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

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

Для нахождения неизвестных коэффициентов на каждом из двух участков кольца необходимо решить следующую систему уравнений:

На первом участке кольца:

.

На втором участке кольца:

.

Решая полученную систему уравнений по правилу Крамера, вычислили главный определитель системы, который оказался равен нулю, а вспомогательный определитель для коэффициента k11 оказался не нулевым, значит, система решений не имеет.

Таким образом, определили, что полином второй степени нельзя использовать для замены аргумента в исходных формулах перемещений точек кольца при заданных начальных условиях.

При использовании в качестве замены аргумента полинома третьей степени имеем:

- на первом участке кольца:

.

- на втором участке кольца:

.

Для нахождения коэффициентов k0, k1, k2, k3 на каждом из двух участков кольца, заключенных между точками опор, необходимо решить следующую систему уравнений:

Главный определитель полученной системы не равен нулю, поэтому по правилу Крамера эта система имеет единственное решение.

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

На первом участке кольца:

.

На втором участке кольца:

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

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

Кинетическая энергия не вращающегося кольца вычисляется по формуле:

Потенциальная энергия вычисляется по формуле:

.

где

J - момент инерции сечения кольца

r - радиус средней линии кольца,

h - толщина кольца,

a - ширина кольца,

E - модуль упругости кольца,

с - удельная плотность кольца.

Предполагаем, что радиус кольца остается постоянным. [3]

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

Литература

1. Серов, В.В., Дуганов, В.Я. Определение собственных частот колебаний массивного вращающегося кольца [Электронный ресурс], 2006 г. Режим доступа: http://www.rusnauka.com/ONG_2006/Matemathics/17845.doc.htm (доступ свободный). Загл. с экрана. Яз. рус.

2. Полунин, А.И. Математическое моделирование динамики упругой вращающейся оболочки с опорами [Текст]: Монография / А.И. Полунин. Белгород: Изд-во БГТУ, 2008.

3. Бабаков И.М. Теория колебаний [Текст]: М.: «Наука», 1968. 560 с.

Размещено на Allbest.ru

...

Подобные документы

  • Решение дифференциальных уравнений с использованием классических алгоритмов численных методов Эйлера и Рунге-Кутта 4-го порядка. Команды, используемые при решении обыкновенных дифференциальных уравнений в системе вычислений. Результат работы программы.

    курсовая работа [226,6 K], добавлен 05.04.2013

  • Численные методы решения задач. Решение алгебраических и трансцендентных уравнений. Уточнение корня по методу половинного деления. Решение систем линейных уравнений методом итераций. Методы решения дифференциальных уравнений. Решение транспортной задачи.

    курсовая работа [149,7 K], добавлен 16.11.2008

  • Суть метода Рунге-Кутта и его свойства. Решение дифференциальных уравнений первого порядка. Вычислительный блок Given/Odesolve. Встроенные функции rkfixed, Rkadapt, Bulstoer. Решения линейных алгебраических уравнений в среде MathCad и Microsoft Excel.

    курсовая работа [1,1 M], добавлен 02.06.2014

  • Применение итерационных методов численного решения системы линейных алгебраических уравнений при вычислении на ЭВМ. Математические и алгоритмические основы решения задачи, метод Гаусса. Функциональные модели и блок-схемы, программная реализация решения.

    курсовая работа [527,5 K], добавлен 25.01.2010

  • Требования к языкам программирования, их эффективность, лаконичность, ясность, реальные возможности. Создание языка С#. Применение систем линейных алгебраических уравнений для практических задач, сущность и особенности метода Крамера для их решения.

    курсовая работа [118,1 K], добавлен 13.11.2009

  • Системы линейных алгебраических уравнений. Код программы для решения систем линейных алгебраических уравнений. Математические и алгоритмические основы решения задачи методом Гаусса. Программная реализация решения. Алгоритмы запоминания коэффициентов.

    лабораторная работа [23,5 K], добавлен 23.09.2014

  • Численные методы решения нелинейных уравнений, систем линейных и нелинейных алгебраических уравнений, дифференциальных уравнений, определенных интегралов. Методы аппроксимации дискретных функций и методы решения задач линейного программирования.

    методичка [185,7 K], добавлен 18.12.2014

  • Итерационные методы решения нелинейных уравнений, системы линейных алгебраических уравнений (СЛАУ). Решение нелинейных уравнений методом интерполирования. Программная реализация итерационных методов решения СЛАУ. Практическое применение метода Эйлера.

    курсовая работа [1,6 M], добавлен 20.01.2010

  • Решение системы обыкновенных дифференциальных уравнений в программе Matlab. Применение метода Рунге–Кутты. Априорный выбор шага интегрирования. Построение трехмерного графика движения точки в декартовой системе координат и создание видеофайла формата AVI.

    контрольная работа [602,8 K], добавлен 04.05.2015

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

    курсовая работа [939,4 K], добавлен 16.01.2014

  • Использование MS Excel для математических расчетов. Описание численных методов решения системы линейных алгебраических уравнений. Решение систем линейных алгебраических уравнений с методами Крамера и Зейделя и с помощью табличного процессора MS Excel.

    курсовая работа [1,6 M], добавлен 14.02.2021

  • Решение системы линейных алгебраических уравнений методом Гаусса с выборкой ведущего элемента. Изучение особенности программной реализации алгоритма, составленной средствами разработки Microsoft Visual Studio. Проведение сложения и умножения двух матриц.

    курсовая работа [702,6 K], добавлен 22.03.2015

  • Алгоритм решения систем линейных уравнений методом Гаусса, его этапы. Система уравнений для определения коэффициентов сплайна, представляющая собой частный случай систем линейных алгебраических уравнений. Программная реализация, тестовый пример.

    курсовая работа [431,8 K], добавлен 15.06.2013

  • Решение линейных дифференциальных уравнений численными и символьными методами в рамках пакета компьютерной математики MathCAD. Сравнения результов решений и применение их при исследовании функционирования автоматических систем и электрических агрегатов.

    контрольная работа [51,5 K], добавлен 07.05.2009

  • Рассмотрение двух способов решения систем линейных алгебраических уравнений: точечные и приближенные. Использование при программировании метода Гаусса с выбором главного элемента в матрице и принципа Зейделя. Применение простой итерации решения уравнения.

    курсовая работа [879,8 K], добавлен 05.06.2012

  • Численные методы решения задачи Коши для обыкновенных дифференциальных уравнений: Эйлера, Рунге-Кутта, Адамса и Рунге. Техники приближенного решения данных уравнений: метод конечных разностей, разностной прогонки, коллокаций; анализ результатов.

    курсовая работа [532,9 K], добавлен 14.01.2014

  • Сферы использования компьютеров, сущность и языки программирования. Применение модифицированного метода Гаусса и расширенной матрицы для решения системы линейных алгебраических уравнений (СЛАУ). Разработка программы, системные требования для ее работы.

    курсовая работа [657,1 K], добавлен 09.01.2014

  • Постановка задачи, математические и алгоритмические основы решения системы линейных алгебраических уравнений. Решение системы данных уравнений методом Гаусса с выбором главного элемента по столбцу. Функциональные модели и блок-схемы решения задачи.

    курсовая работа [428,9 K], добавлен 25.01.2010

  • Решение уравнения методом половинного деления. Программа в Matlab для уравнения (x-2)cos(x)=1. Решение нелинейных уравнений методом Ньютона. Интерполяция заданной функции. Решение системы линейных алгебраических и обыкновенных дифференциальных уравнений.

    курсовая работа [1,4 M], добавлен 15.08.2012

  • Возможности математического пакета MathCad в среде Windows 98 для использования матричной алгебры и решения системы линейных алгебраических уравнений. Методы решения систем линейных алгебраических уравнений. Сравнение метода Гаусса с методом MathCad.

    практическая работа [62,6 K], добавлен 05.12.2009

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.