ПРИМЕНЕНИЕ МЕТОДА КОНЕЧНЫХ РАЗНОСТЕЙ ДЛЯ РЕШЕНИЯ ДИНАМИЧЕСКИХ ЗАДАЧ

Научная статья
DOI:
https://doi.org/10.23670/IRJ.2017.60.094
Выпуск: № 6 (60), 2017
Опубликована:
2017/06/19
PDF

Литвинова Э.В.

ORCID: 0000–0001–5549–5627, кандидат технических наук,

Академия строительства и архитектуры

ФГАОУ ВО «КФУ имени В.И. Вернадского»

ПРИМЕНЕНИЕ МЕТОДА КОНЕЧНЫХ РАЗНОСТЕЙ ДЛЯ РЕШЕНИЯ ДИНАМИЧЕСКИХ ЗАДАЧ

Аннотация

Предлагается применение метода конечных разностей для решения динамических задач строительной механики на примере динамического изгиба жестко защемленной призматической балки под действием внезапно приложенной равномерно распределенной нагрузки, неизменной во времени. В методе конечных разностей область непрерывного изменения аргумента заменяется конечным (дискретным) множеством узлов, называемым сеткой. Метод конечных разностей – сеточный метод, основанный на замене производных разностными отношениями, прост и удобен для вычислений.

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

Litvinova E.V.

ORCID: 0000–0001–5549–5627, PhD in Engineering, Academy of Construction and Architecture V.I. Vernadsky Crimean Federal University

APPLICATION OF FINITE DIFFERENCES METHOD FOR DYNAMIC PROBLEMS SOLUTION

Abstract

The article proposes the application of finite differences method for solving the dynamic problems of building mechanics with the help of dynamic bending of a rigidly clamped prismatic beam under the impact of a suddenly applied uniformly distributed load, unchanged in time. In the method of finite differences, the domain of continuous change of the argument is replaced by a finite (discrete) set of nodes, called a grid. Finite differences method – the grid method, based on the replacement of derivatives by differencing ratios, is simple and convenient for calculations.

Keywords: finite difference method, grid functions, differencing scheme, dynamic load application.

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

Рассмотрим один из методов дискретизации аналитических моделей функционирования строительных конструкций – метод конечных разностей (МКР).

Пусть на отрезке image001 требуется найти решение некоторого уравнения

image002

где L – некоторый дифференциальный оператор.

Будем предполагать, что решение image003 уравнения (1) имеет необходимое число производных. Для его решения применим МКР [1 – 3].

В МКР область непрерывного изменения аргумента заменяется конечным (дискретным) множеством узлов, называемым сеткой. Сетку вводят следующим образом. Отрезок изменения аргумента image004 разбивают на n равных частей длиной image005 каждая. Сами точки image006 называются узлами сетки. На сетке вместо функции непрерывного аргумента image003 рассматриваются функции дискретного аргумента image007, определяемые в узлах сетки и называемые сеточными функциями. Производные заменяются (аппроксимируются) соответствующими разностными отношениями, то есть линейной комбинацией значений сеточной функции в нескольких узлах [1, 3].

Необходимо найти разностный аналог производной функции image003 в касательной к кривой в точке. Это можно сделать с помощью замены истинной производной следующими разностными отношениями [3]:

image008

Первое выражение в (2) определяет левую одностороннюю разностную производную, второе выражение – правую одностороннюю разностную производную, третье – двусторонняя разностная производная.

Очевидно, располагая выражением для первой производной, можно выписать разностные аналоги и для производных высших порядков [3]:

image009

image010

image011

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

Часто при решении задач с помощью МКР (главным образом при аппроксимации граничных условий) используют правые и левые разности, аналогичные по точности аппроксимации центральным разностям. Приведем без вывода наиболее характерные из них [3]:

image012

image013

image014

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

image015

Рис. 1 – Схема крепления и нагружения балки

Применение МКР для решения динамических задач покажем на примере динамического изгиба жестко защемленной призматической балки (рис. 1), под действием внезапно приложенной равномерно распределенной нагрузки, неизменной во времени image016 Здесь image017 – единичная функция Хэвисайда

image018

Схема задачи, показанная на рис. 1, отличается в динамическом приложении нагрузки. До приложения нагрузки балка неподвижна и не имеет прогиба.

Динамика системы с внезапным приложением нагрузки описывается дифференциальным уравнением

image019

Здесь: image020 – погонная масса; ρ – плотность материала; А – площадь поперечного сечения; Е – модуль Юнга; J – момент инерции поперечного сечения; w – величина прогиба балки; x – горизонтальная координата вдоль оси балки.

Если предположить, что колебания системы происходят при наличии сил сопротивления, пропорциональных скорости движения поперечных сечений балки, то уравнение (3) можно записать так

image021

где α – коэффициент сопротивления. Уравнение (3) необходимо дополнить начальными и граничными условиями:

image022

image023

Применим для решения краевой задачи (4) – (6) МКР. Выберем на плоскости image024 ( image025 – координата, image026 – время) прямоугольную равномерную сетку с шагами image027 и image028. Искомые сеточные функции image029 (j – номер шага по координате x, i – номер шага по времени t, image030 – итоговая функция прогиба балки). Для аппроксимации производных, входящих в уравнение (4), а также граничных и начальных условий применим формулы центральных конечных разностей:

 image031

Разделим (4) на массу m

image032

Подставим выражения (7) и граничные условия (5) и (6) в уравнение (8). Получим систему алгебраических уравнений относительно сеточных функцийimage029  внутри области image033 в момент времени i

image034

image035

Рис. 2 – Двумерная разностная сетка

Как видно из выражения (9), в точке на плоскости image036 с координатами (j, i) ее прогиб формируется с помощью семи узловых точек, образующих на плоскости «крест» (рис. 2). При этом значение узловой переменной на новом шаге по времени image037 может быть явно выражено из уравнения (9). Поэтому разностную схему (9), построенную с помощью соотношений (7) называют явной разностной схемой типа «крест»

image038

Здесь

image039

Граничные условия (6) не зависят от времени и могут быть записаны следующим образом:

image040                 (11)

Из начальных условий (5) имеем

image041

Для замыкания системы нам необходимо получить значения узловых перемещений image029 на первом шаге по времени (i = 1). Они могут быть найдены при помощи разложения в ряд Тейлора с учетом начальных условий (10)

image042

Для определения image043 воспользуемся уравнением (8). Тогда

image044

Отсюда следует, что image045

Таким образом, имеем разностное уравнение (10), граничные условия (11) и начальные условия для первого и второго слоя разностной схемы

image046

которые аппроксимируют краевую задачу (8) с граничными и начальными условиями с точностью до величины  image047

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

Кроме того, не будем учитывать внутреннее сопротивление конструкции. Тогда уравнение движения балки (4) примет вид

image048

image049

Здесь image050

Зададим исходные данные: image051 – погонная нагрузка, image052 − длина балки, image053 − модуль упругости, image054 − допускаемое напряжение.

Максимальный изгибающий момент в месте жесткой заделки балки image055 Допускаемый момент сопротивления для данных условий составляет image056 Пусть балка имеет двутавровое сечение. Тогда из таблицы сортаментов выбираем двутавр № 18, для которого: image057 − площадь поперечного сечения. Погонная масса балки при значении плотности материала image058 составляет image059 Прогиб балки при статической нагрузке в фиксированных сечениях составляет:

image060

Перейдем к динамическому расчету. Разобьем балку по длине на пять равных отрезков 04-07-2017 12-58-15

Применим МКР к уравнению (12), получим четыре уравнения:

image062

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

Вполне естественно, что максимальный изгибающий момент в заделке балки должен быть в два раза больше, т. е. image064 image065

С целью обеспечения прочности балки при заданной динамической нагрузке необходимо определить допускаемый момент сопротивления сечения image066image067. Из сортамента выбираем двутавр № 22а: image068 При таких условиях решение динамической задачи дает:  image069

При этом  image070

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

Рассматриваемая схема проста и удобна для вычислений.

 

Список литературы / References

  1. Курамшина Р.П. Численные методы в строительстве и их реализация: учебное пособие / Р.П. Курамшина. – Братск: ГОУВПО «Братский государственный университет», 2010. – 104 с.
  2. Чемодуров В.Т. Методы теории планирования эксперимента в решении технических задач: монография / В.Т. Чемодуров, В.В. Жигна. – Симферополь: ИТ «АРИАЛ», 2012. – 110 с.
  3. Чемодуров В.Т. Численные методы в строительстве: учебное пособие / В.Т. Чемодуров, М.С. Сеитжелилов. – Симферополь: ИТ «АРИАЛ», 2016. – 112 с.

Список литературы на английском языке / References in English

  1. Kuramshina R.P. Chislennye metody v stroitel'stve i ih realizacija: uchebnoe posobie [Numeral methods in building and their realization: train aid] / R.P. Kuramshina. – Bratsk: GOUVPO «Bratskij gosudarstvennyj universitet» [Bratsk State University], 2010. – 104 p. [in Russian]
  2. Chemodurov V.T. Metody teorii planirovanija jeksperimenta v reshenii tehnicheskih zadach: monografija [Methods of theory of planning of experiment are in the decision of technical tasks: monograph] / V.T. Chemodurov, V.V. Zhigna. – Simferopol': IT «ARIAL», 2012. – 110 p. [in Russian]
  3. Chemodurov V.T. Chislennye metody v stroitel'stve (dlja studentov vseh special'nostej ochnoj i zaochnoj form obuchenija): uchebnoe posobie [Numeral methods are in building: train aid] / V.T. Chemodurov, M.S. Seitzhelilov. – Simferopol': IT «ARIAL», 2016. – 112 p. [in Russian]