Научный журнал
Международный журнал экспериментального образования
ISSN 2618–7159
ИФ РИНЦ = 0,425

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

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

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

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

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

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

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

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

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

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

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

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

Musayev05.wmf

Musayev06.wmf Musayev07.wmf

Musayev08.wmf

Musayev09.wmf Musayev10.wmf

Musayev11.wmf Musayev12.wmf

Musayev13.wmf Musayev14.wmf (1)

где σx, σy и τxy – компоненты тензора упругих напряжений; εx, εy и γxy – компоненты тензора упругих деформаций; u и v – составляющие вектора упругих перемещений вдоль осей OX и OY соответственно; ρ – плотность материала; Musayev15.wmf – скорость продольной упругой волны; Musayev16.wmf – скорость поперечной упругой волны; ν – коэффициент Пуассона; E – модуль упругости; Musayev17.wmf – граничный контур тела Γ.

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

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

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

Musayev18.wmf Musayev19.wmf Musayev20.wmf (2)

где Musayev21.wmf – диагональная матрица инерции; Musayev22.wmf – матрица жесткости; Musayev23.wmf – вектор узловых упругих перемещений; Musayev24.wmf – вектор узловых упругих скоростей перемещений; Musayev25.wmf – вектор узловых упругих ускорений; Musayev26.wmf – вектор внешних узловых упругих сил.

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

Musayev27.wmf

Musayev28.wmf (3)

Шаг по временной переменной координате Δt выбирается из следующего соотношения

Musayev29.wmf (i = 1, 2, 3, ...), (4)

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

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

pic_6.tif

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

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

pic_7.tif

Рис. 2. Изменение упругого контурного напряжения Musayev30.wmf во времени t/Δt в точке A1: 1 – в задаче без полости; 2 – в задаче с полостью (соотношение ширины к высоте один к двенадцати)

pic_8.tif

Рис. 3. Изменение упругого контурного напряжения Musayev31.wmf во времени t/Δt в точке A2: 1 – в задаче без полости; 2 – в задаче с полостью (соотношение ширины к высоте один к двенадцати)

pic_9.tif

Рис. 4. Изменение упругого контурного напряжения Musayev32.wmf во времени t/Δt в точке A3: 1 – в задаче без полости; 2 – в задаче с полостью (соотношение ширины к высоте один к двенадцати)

pic_10.tif

Рис. 5. Изменение упругого контурного напряжения Musayev33.wmf во времени t/Δt в точке A4: 1 – в задаче без полости; 2 – в задаче с полостью (соотношение ширины к высоте один к двенадцати)

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

В точке F перпендикулярно свободной поверхности ABEFG приложено сосредоточенное нормальное напряжение σy (рис. 1), которое при 0 ≤ n ≤ 10 (n = t/Δt) изменяется линейно от 0 до P, а при 10 ≤ n ≤ 20 от P до 0 (P = σ0, σ0 = –1 МПа (–1 кгс/см2)).

Граничные условия для контура GHIA при t > 0 Musayev34.wmf. Отраженные волны от контура GHIA не доходят до исследуемых точек при 0 ≤ n ≤ 200. Контур ABCDEFG свободен от нагрузок, кроме точки F, где приложено сосредоточенное упругое нормальное напряжение σy.

Расчеты проведены при следующих исходных данных: 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 неизвестных.

Результаты расчетов для контурного напряжения Musayev35.wmf Musayev36.wmf во времени n получены в точках A1–A4 (рис. 1), находящихся на свободной поверхности упругой полуплоскости.

На рис. 2–5 приведены контурные напряжения Musayev37.wmf во времени n в точках A1–A4: 1 – в задаче без полости; 2 – в задаче с полостью (соотношение ширины к высоте один к двенадцати).

Выводы

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

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

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

4. Рассмотрена задача о воздействии сосредоточенной взрывной волны на свободной поверхности упругой полуплоскости с полостью (соотношение ширины к высоте один к двенадцати).

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


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

Мусаев В.К. МОДЕЛИРОВАНИЕ НЕСТАЦИОНАРНЫХ УПРУГИХ ВОЛН НАПРЯЖЕНИЙ В УПРУГОЙ ПОЛУПЛОСКОСТИ С ПОЛОСТЬЮ (СООТНОШЕНИЕ ШИРИНЫ К ВЫСОТЕ ОДИН К ДВЕНАДЦАТИ) ПРИ СОСРЕДОТОЧЕННОМ ВЗРЫВНОМ ВОЗДЕЙСТВИИ // Международный журнал экспериментального образования. – 2016. – № 12-1. – С. 139-142;
URL: https://expeducation.ru/ru/article/view?id=10843 (дата обращения: 28.03.2024).

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

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