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

ISSN 2618–7159
ИФ РИНЦ = 0,442

ДОСТОВЕРНОСТЬ ЧИСЛЕННОГО МЕТОДА, АЛГОРИТМА И КОМПЛЕКСА ПРОГРАММ МУСАЕВА В.К. ПРИ РЕШЕНИИ ЗАДАЧИ О РАСПРОСТРАНЕНИИ ПЛОСКИХ ПРОДОЛЬНЫХ УПРУГИХ ВОЛН (ВОСХОДЯЩАЯ ЧАСТЬ – ЛИНЕЙНАЯ, НИСХОДЯЩАЯ ЧАСТЬ – ЧЕТВЕРТЬ КРУГА) В ПОЛУПЛОСКОСТИ

Дикова Е.В. 1
1 Московский политехнический университет
Приводится некоторая информация моделирования нестационарных упругих волн в полуплоскости при импульсном воздействии (восходящая часть – линейная, нисходящая часть – четверть круга). Для решения поставленной задачи применяются волновые уравнения механики деформируемого твердого тела. На основе метода конечных элементов в перемещениях разработаны численный метод, алгоритм и комплекс программ Мусаева В.К. Для решения линейных плоских двумерных задач, которые позволяют решать задачи при нестационарных динамических воздействиях на сложные системы. При разработке комплекса программ использовался алгоритмический язык Фортран-90. Применяется квазирегулярный подход при аппроксимации исследуемой области. Приводится нормальное напряжение в характерной точке упругой полуплоскости.
достоверность
верификация
численный метод Мусаева В.К.
алгоритм Мусаева В.К.
1. Мусаев В.К. О достоверности компьютерного моделирования нестационарных упругих волн напряжений в деформируемых телах сложной формы // Международный журнал прикладных и фундаментальных исследований. – 2014. – № 11. – С. 10–14.
2. Мусаев В.К. Моделирование нестационарных упругих волн напряжений в деформируемых областях с помощью метода конечных элементов в перемещениях // Современные наукоемкие технологии. – 2014. – № 12 (1). – С. 28–32. 
3. Мусаев В.К. Оценка точности и достоверности численного моделирования при решении задач об отражении и интерференции нестационарных упругих волн напряжений // Успехи современного естествознания. – 2015. – № 1 (часть 7). – С. 1184–1187.
4. Мусаев В.К. Численное решение задачи о распространении нестационарных упругих волн напряжений в подкрепленном круглом отверстии // Современные наукоемкие технологии. – 2015. – № 2. – С. 93–97. 
5. Мусаев В.К. Решение задачи о распространении плоских продольных волн в виде импульсного воздействия // Международный журнал экспериментального образования. – 2015. – № 4 (часть 2). – С. 326–330.
6. Мусаев В.К. Моделирование нестационарных стоячих упругих волн в бесконечной полосе при воздействии в виде треугольного импульса // Международный журнал прикладных и фундаментальных исследований. – 2015. – № 11 (часть 2). – С. 248–251.
7. Мусаев В.К. Численное моделирование плоских продольных волн в виде импульсного воздействия (восходящая часть – четверть круга, средняя – горизонтальная, нисходящая – линейная) в упругой полуплоскости // Международный журнал экспериментального образования. – 2015. – № 11 (часть 2). – С. 222–226.
8. Мусаев В.К. Оценка физической достоверности и математической точности численного моделирования плоских нестационарных упругих волн напряжений в виде импульсного воздействия (функция Хевисайда) в полуплоскости // Международный журнал прикладных и фундаментальных исследований. – 2016. – № 11–1. – С. 49–52. 
9. Мусаев В.К. Оценка математической точности и физической достоверности численного моделирования плоских нестационарных упругих волн напряжений (дельта функция) в полуплоскости // Международный журнал прикладных и фундаментальных исследований. – 2016. – № 11–2. – С. 232–235.
10. Мусаев В.К. Оценка точности компьютерного моделирования плоских нестационарных волн напряжений (прямоугольный импульс) в упругой полуплоскости // Международный журнал прикладных и фундаментальных исследований. – 2016. – № 11–2. – С. 236–239. 

Постановка задачи при нестационарных волновых воздействиях

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

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

diko002.wmf, diko003.wmf,

diko004.wmf,

diko005.wmf,

diko006.wmf,

diko007.wmf,

diko008.wmf, diko009.wmf, diko010.wmf,

diko011.wmf, (1)

где diko012.wmf, diko013.wmf и diko014.wmf – компоненты тензора упругих напряжений; diko015.wmf, diko016.wmf и diko017.wmf – компоненты тензора упругих деформаций; u и v – составляющие вектора упругих перемещений вдоль осей OX и OY соответственно; diko018.wmf – плотность материала; diko019.wmf – скорость продольной упругой волны; diko020.wmf – скорость поперечной упругой волны; diko021.wmf – коэффициент Пуассона; E – модуль упругости; diko022.wmf – граничный контур тела Г.

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

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

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

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

diko024.wmf, diko025.wmf,

diko026.wmf, (2)

где diko027.wmf – диагональная матрица инерции; diko028.wmf – матрица жесткости; diko029.wmf – вектор узловых упругих перемещений; diko030.wmf – вектор узловых упругих скоростей перемещений; diko031.wmf – вектор узловых упругих ускорений; diko032.wmf – вектор внешних узловых упругих сил.

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

diko033.wmf,

diko034.wmf. (3)

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

diko036.wmf diko037.wmf, (4)

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

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

Моделирование распространения плоских продольных упругих волн в полуплоскости

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

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

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

На границе полуплоскости AB приложено нормальное напряжение diko039.wmf, которое при diko040.wmf (diko041.wmf) изменяется от 0 до P и при diko042.wmf изменяется от P до 0 (diko043.wmf, diko044.wmf МПа (–1 кгс/см2)). Граничные условия для контура BCDA при t>0 diko046.wmf. Отраженные волны от контура BCDA не доходят до исследуемых точек при diko047.wmf. Расчеты проведены при следующих исходных данных: diko048.wmf; t = 0,125*10–4 с; E = 22,366*101 МПа (22,366*102 кгс/см2); n= 0,22; r= 1,469*103 кг/м3 (1,469*10–6 кгс•с2/см4); Cp= 400 м/с; Cs= 250 м/с. Исследуемая расчетная область имеет 20862 узловых точек. Решается система уравнений из 83448 неизвестных.

На рис. 3 представлено изменение нормального напряжения diko049.wmf (diko050.wmf) во времени n в точке B1.

dik1.tif

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

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

На фронте плоской продольной волны имеются следующие аналитические зависимости для плоского напряженного состояния diko051.wmf и diko052.wmf.

Отсюда видим, что точное решение задачи соответствует воздействию diko053.wmf (рис. 2).

dik2.tif

Рис. 2. Импульсное воздействие (восходящая часть – линейная, нисходящая – четверть круга)

dik3.tif

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

Для упругих нормальных напряжений diko056.wmf и diko057.wmf имеется хорошее качественное и количественное согласование с результатами точного решения. Таким образом, можно сделать вывод, что на точность численного решения оказывает влияние аппроксимация воздействия.

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

Выводы

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

Решена задача о распространении плоских продольных волн в виде импульсного воздействия в упругой полуплоскости. Исследуемая расчетная область имеет 20862 узловых точек. Решается система уравнений из 83448 неизвестных.

Для решения поставленной задачи используется импульсное воздействие (восходящая часть – линейная, нисходящая часть – линейная).

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

Автор выражает благодарность Мусаеву В.К. за внимание к работе.


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

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

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

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