Научный журнал
Международный журнал прикладных и фундаментальных исследований
ISSN 1996-3955
ИФ РИНЦ = 0,593

ОЦЕНКА МАТЕМАТИЧЕСКОЙ ТОЧНОСТИ И ФИЗИЧЕСКОЙ ДОСТОВЕРНОСТИ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ ПЛОСКИХ НЕСТАЦИОНАРНЫХ УПРУГИХ ВОЛН НАПРЯЖЕНИЙ (ДЕЛЬТА ФУНКЦИЯ) В ПОЛУПЛОСКОСТИ

Мусаев В.К. 1
1 Московский государственный университет путей сообщения Императора Николая II
Рассмотрена некоторая информация о моделировании плоских нестационарных упругих волн напряжений в упругой полуплоскости. Для решения поставленной задачи применяются уравнения волновой теории упругости. На основе численного моделирования уравнений волновой теории упругости разработаны методика, алгоритм и комплекс программ для решения линейных нестационарных динамических задач. Основные соотношения метода конечных элементов получены с помощью динамического равновесия. Получена явная двухслойная схема. При решении сложных задач возникают проблемы оценки достоверности полученных результатов. В работе рассматривается оценка точности и достоверности результатов численного моделирования волн напряжений при распространении плоской нестационарной упругой волны в полуплоскости. В качестве воздействия применяется дельта функция или треугольный импульс. Решается система уравнений из 59048 неизвестных.
вычислительная механика
численный метод
алгоритм
комплекс программ Мусаева В.К.
нестационарные упругие волны
плоская продольная нестационарная волна
фундаментальное воздействие
распространение плоских волн
исследуемая расчетная область
дельта функция
треугольный импульс
верификация
математическая точность
физическая достоверность
точность численного метода
техносферная безопасность
1. Мусаев В.К. О достоверности компьютерного моделирования нестационарных упругих волн напряжений в деформируемых телах сложной формы // Международный журнал прикладных и фундаментальных исследований. – 2014. – № 11. – С. 10–14.
2. Мусаев В.К. Моделирование нестационарных упругих волн напряжений в деформируемых областях с помощью метода конечных элементов в перемещениях // Современные наукоемкие технологии. – 2014. – № 12 (1). – С. 28–32.
3. Мусаев В.К. Моделирование безопасности по несущей способности дымовых труб с основанием при взрыве атомной бомбы в Нагасаки // Международный журнал прикладных и фундаментальных исследований. – 2014. – № 12. – С. 198–203.
4. Мусаев В.К. Математическое моделирование поверхностных волн напряжений в задаче Лэмба при воздействии в виде дельта функции // Международный журнал прикладных и фундаментальных исследований. – 2015. – № 2–1. – С. 25–29.
5. Мусаев В.К. Определение нестационарного напряженного состояния при вертикальном сосредоточенном взрывном воздействии на набережной речного порта с незаполненным водным объектом // Современные наукоемкие технологии. – 2015. – № 2. – С. 88–92.
6. Мусаев В.К. Численное решение задачи о распространении нестационарных упругих волн напряжений в подкрепленном круглом отверстии // Современные наукоемкие технологии. – 2015. – № 2. – С. 93–97.
7. Мусаев В.К. Решение задачи о сосредоточенном упругом взрывном воздействии в виде дельта функции на набережной речного порта с заполненным водным объектом на 25 % // Современные наукоемкие технологии. – 2015. – № 4. – С. 52–55.
8. Мусаев В.К. Решение задачи о распространении плоских продольных волн в виде импульсного воздействия // Международный журнал экспериментального образования. – 2015. – № 4–2. – С. 326–330.
9. Мусаев В.К. Исследования устойчивости явной двухслойной линейной конечноэлементной схемы для внутренних узловых точек на равномерной прямоугольной сетке // Современные наукоемкие технологии. – 2015. – № 5. – С. 39–42.
10. Мусаев В.К. Численное моделирование плоских продольных волн в виде импульсного воздействия (восходящая часть – четверть круга, средняя – горизонтальная, нисходящая – линейная) в упругой полуплоскости // Международный журнал экспериментального образования. – 2015. – № 11–2. – С. 222–226.

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

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

Расчеты проводились при следующих единицах измерения: килограмм-сила (кгс); сантиметр (см); секунда (с). Для перехода в другие единицы измерения были приняты следующие допущения: 1 кгс/см2 ≈ 0,1 МПа; 1 кгс с2/см4 ≈ 109 кг/м3.

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

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

Постановка задачи с начальными и граничными условиями

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

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

Точные уравнения двумерной (плоское напряженное состояние) динамической теории упругости имеют вид

musaev01.wmf, musaev02.wmf,

musaev03.wmf,

musaev04.wmf,

musaev05.wmf,

musaev06.wmf,

musaev07.wmf, musaev08.wmf, musaev09.wmf,

musaev10.wmf, (1)

где σx, σy и τxy – компоненты тензора упругих напряжений; εx, εy и γxy – компоненты тензора упругих деформаций; u и v – составляющие вектора упругих перемещений вдоль осей OX и OY соответственно; ρ – плотность материала;

musaev11.wmf – скорость продольной упругой волны;

musaev12.wmf – скорость поперечной упругой волны;

ν – коэффициент Пуассона; E – модуль упругости; musaev13.wmf – граничный контур тела Г.

Систему (1) в области, занимаемой телом Г, следует интегрировать при начальных и граничных условиях.

Разработка методики и алгоритма

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

Принимая во внимание определение матрицы жесткости, вектора инерции и вектора внешних сил для тела Г, записываем приближенное значение уравнения движения в теории упругости

musaev14.wmf, musaev15.wmf, musaev16.wmf, (2)

где musaev17.wmf – диагональная матрица инерции; musaev18.wmf – матрица жесткости; musaev19.wmf – вектор узловых упругих перемещений; musaev20.wmf – вектор узловых упругих скоростей перемещений; musaev21.wmf – вектор узловых упругих ускорений; musaev22.wmf – вектор внешних узловых упругих сил.

Соотношение (2) система линейных обыкновенных дифференциальных уравнений второго порядка в перемещениях с начальными условиями.

Таким образом, с помощью метода конечных элементов в перемещениях, линейную задачу с начальными и граничными условиями (1) привели к линейной задаче Коши (2).

Для интегрирования уравнения (2) конечноэлементным вариантом метода Галеркина приведем его к следующему виду

musaev23.wmf, musaev24.wmf. (3)

Интегрируя по временной координате соотношение (3) с помощью конечноэлементного варианта метода Галеркина, получим двумерную явную двухслойную конечноэлементную линейную схему в перемещениях для внутренних и граничных узловых точек

musaev25.wmf,

musaev26.wmf, (4)

где Δt – шаг по временной координате.

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

Общая теория численных уравнений математической физики требует для этого наложение определенных условий на отношение шагов по временной координате Δt и по пространственным координатам, а именно

musaev27.wmf musaev28.wmf, (5)

где Δl – длина стороны конечного элемента.

Для исследуемой области, состоящей из материалов с разными физическими свойствами, выбирается минимальный шаг по временной координате (5).

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

Численное моделирование импульсного воздействия в упругой полуплоскости

Рассмотрим задачу о воздействии плоской продольной волны в виде дельта функции (рис. 2) на упругую полуплоскость (рис. 1).

На границе полуплоскости AB приложено нормальное напряжение σy, которое при 0 ≤ n ≤ 10 (n = t/Δt) изменяется линейно от 0 до P, а при 10 ≤ n ≤ 20 от P до 0 (P = σ0, σ0 = – 0,1 МПа (– 1 кгс/см2)). Граничные условия для контура BCDA при t > 0 musaev29.wmf. Отраженные волны от контура BCDA не доходят до исследуемых точек при 0 ≤ n ≤ 100. Расчеты проведены при следующих исходных данных: H = Δx = Δy; Δt = 1,393•10-6 с; E = 3,15•104 МПа (3,15•105 кгс/см2); ν = 0,2; ρ = 0,255•104 кг/м3 (0,255•10-5 кгс•с2/см4); Cp = 3587 м/с; Cs = 2269 м/с. Решается система уравнений из 59048 неизвестных.

mus1.tif

Рис. 1. Постановка задачи о распространении плоских продольных нестационарных упругих волн в полуплоскости

mus2.tif

Рис. 2. Воздействие в виде дельта функции (треугольный импульс)

На рис. 3-6 представлено изменение нормального напряжения musaev30.wmf (musaev31.wmf) во времени n в точках B1–B4: 1 – численное решение; 2 – аналитическое решение.

mus3.tif

Рис. 3. Изменение упругого нормального напряжения musaev32.wmf во времени t/Δt в точке B1: 1 – численное решение; 2 – аналитическое решение

mus4.tif

Рис. 4. Изменение упругого нормального напряжения musaev33.wmf во времени t/Δt в точке B2: 1 – численное решение; 2 – аналитическое решение

mus5.tif

Рис. 5. Изменение упругого нормального напряжения musaev34.wmf во времени t/Δt в точке B3: 1 – численное решение; 2 – аналитическое решение

mus6.tif

Рис. 6. Изменение упругого нормального напряжения musaev35.wmf во времени t/Δt в точке B4: 1 – численное решение; 2 – аналитическое решение

В данном случае можно использовать условия на фронте плоской волны. Предположим, что от некоторых точек упругой среды производится какое-то возмущение. Тогда из этих точек во все стороны начинают излучаться волны. На некотором расстоянии от центра возмущения рассматриваемые волны можно представить как плоские. Тогда все частицы движутся параллельно направлению распространения волны. Такие волны принято считать плоскими. На фронте плоской продольной волны имеется следующая аналитическая зависимость для плоского напряженного состояния musaev36.wmf. Отсюда видим, что точное решение задачи соответствует воздействию σ0 (рис. 2). Для нормального напряжения σy имеется хорошее качественное и количественное совпадение с результатом аналитического решения. На основании проведенных исследований можно сделать вывод о физической достоверности результатов численного решения задач при распространении нестационарных упругих волн в деформируемых телах.

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


Библиографическая ссылка

Мусаев В.К. ОЦЕНКА МАТЕМАТИЧЕСКОЙ ТОЧНОСТИ И ФИЗИЧЕСКОЙ ДОСТОВЕРНОСТИ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ ПЛОСКИХ НЕСТАЦИОНАРНЫХ УПРУГИХ ВОЛН НАПРЯЖЕНИЙ (ДЕЛЬТА ФУНКЦИЯ) В ПОЛУПЛОСКОСТИ // Международный журнал прикладных и фундаментальных исследований. – 2016. – № 11-2. – С. 232-235;
URL: https://applied-research.ru/ru/article/view?id=10470 (дата обращения: 29.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674