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

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

Мусаев В.К. 1
1 МЭСИ
Приводится некоторая информация моделирования нестационарных упругих волн в полуплоскости при импульсном воздействии (первая ветвь: восходящая часть – четверть круга, нисходящая часть – четверть круга; вторая ветвь: восходящая часть – линейная, нисходящая часть – линейная). для решения поставленной задачи применяются волновые уравнения механики деформируемого твердого тела. на основе метода конечных элементов в перемещениях разработаны численный метод, алгоритм и комплекс программ для решения линейных плоских двумерных задач, которые позволяют решать задачи при нестационарных динамических воздействиях на сложные системы. При разработке комплекса программ использовался алгоритмический язык Фортран-90. Применяется квазирегулярный подход при аппроксимации исследуемой области. Приводятся нормальные напряжения в характерных точках упругой полуплоскости.
нестационарные упругие волны
волновые уравнения
динамика сплошных сред
распространение волн
волновая теория
взаимодействие с границами
полуплоскость
неотражающие граничные условия
алгоритмический язык Фортран-90
численный метод
алгоритм
комплекс программ
квазирегулярный подход
аппроксимация исследуемой области
исследуемая область
условия на фронте плоской волны
импульсное воздействие
напряжения на фронте плоской волны
1. Мусаев В.К. Решение задачи дифракции и распространения упругих волн методом конечных элементов // Строительная механика и расчет сооружений. – 1990. – № 4. – С. 74–78.
2. Мусаев В.К. Численное решение волновых задач теории упругости и пластичности // Вестник Российского университета дружбы народов. Серия прикладная математика и информатика, 1997, № 1, С. 87–110.
3. Мусаев В.К. об оценке достоверности и точности численного решения нестационарных динамических задач // Вестник Российского университета дружбы народов. Серия проблемы комплексной безопасности. – 2007. – № 3. – С. 48-60.
4. Мусаев В.К. Оценка достоверности и точности результатов вычислительного эксперимента при решении задач нестационарной волновой теории упругости // Научный журнал проблем комплексной безопасности. – 2009. – № 1. – С. 55–80.
5. Мусаев В.К. О достоверности результатов математического моделирования нестационарных волн напряжений в объектах сложной формы // Строительная механика инженерных конструкций и сооружений. – 2014. – № 3. – С. 71–76.
6. Мусаев В.К. О достоверности компьютерного моделирования нестационарных упругих волн напряжений в деформируемых телах сложной формы // Международный журнал прикладных и фундаментальных исследований. – 2014. – № 11. – С. 10–14.
7. Мусаев В.К. Моделирование нестационарных упругих волн напряжений в деформируемых областях с помощью метода конечных элементов в перемещениях // Современные наукоемкие технологии. – 2014. – № 12 (1). – С. 28–32.
8. Мусаев В.К. Оценка достоверности численного моделирования нестационарных упругих волн напряжений в областях сложной формы // Материалы Всероссийской конференции «Механика и наномеханика структурно-сложных и гетерогенных сред. Успехи, проблемы, перспективы». – М.: ИПРИМ РАН, 2015. – С. 94–103.
9. Мусаев В.К. Численное решение задачи о распространении нестационарных упругих волн напряжений в подкрепленном круглом отверстии // Современные наукоемкие технологии. – 2015, № 2. – С. 93–97.
10. Мусаев В.К. Численное моделирование вертикального сосредоточенного упругого импульсного воздействия в виде дельта функции на границе воздушной и твердой среды с полостью в виде прямоугольника (соотношение ширины к высоте один к пятнадцати) // Международный журнал прикладных и фундаментальных исследований. – 2015. – № 2 (часть 2). – С. 220–223.

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

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

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

musa1.wmf, musa2.wmf, musa3.wmf,

musa4.wmf, musa5.wmf, musa6.wmf

musa7.wmf, musa8.wmf, musa9.wmf, musa10.wmf, (1)

где musa11.wmf, musa12.wmf и musa13.wmf – компоненты тензора упругих напряжений; musa14.wmf, musa15.wmf и musa16.wmf – компоненты тензора упругих деформаций; musa17.wmf и musa18.wmf – составляющие вектора упругих перемещений вдоль осей musa19.wmf и musa20.wmf соответственно; musa21.wmf – плотность материала; musa22.wmf – скорость продольной упругой волны; musa23.wmf – скорость поперечной упругой волны; n – коэффициент Пуассона; Е – модуль упругости; musa24.wmf – граничный контур тела Г.

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

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

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

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

musa26.wmf, musa27.wmf,

musa28.wmf, (2)

где musa29.wmf – диагональная матрица инерции; musa30.wmf – матрица жесткости; musa31.wmf – вектор узловых упругих перемещений; musa32.wmf – вектор узловых упругих скоростей перемещений; musa33.wmf – вектор узловых упругих ускорений; musa34.wmf – вектор внешних узловых упругих сил.

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

musa35.wmf,

musa36.wmf. (3)

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

musa38.wmf musa39.wmf, (4)

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

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

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

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

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

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

mus1.tif

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

mus2.tif

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

На границе полуплоскости AB приложено нормальное напряжение musa41.wmf, которое при musa42.wmf (musa43.wmf) изменяется от 0 до P, при musa44.wmf изменяется от musa45.wmf до musa46.wmf, при musa47.wmf изменяется от 0 до P и при musa48.wmf изменяется от P до 0 (musa49.wmf, musa50.wmf МПа). Граничные условия для контура BCDA при musa51.wmf musa52.wmf. Отраженные волны от контура BCDA не доходят до исследуемых точек при musa53.wmf. Расчеты проведены при следующих исходных данных: musa54.wmf; ∆t = 9,263⋅10-7 с; E = 7,1⋅104 МПа; n= 0,34; r= 2,755⋅103 кг/м3; Cp= 5398 м/с; Cs= 3078 м/с. Исследуемая расчетная область имеет 62031 узловую точку. Решается система уравнений из 248124 неизвестных.

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

mus3.tif

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

mus4.tif

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

mus5.tif

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

mus6.tif

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

На фронте плоской продольной волны имеются следующие аналитические зависимости для плоского напряженного состояния musa63.wmf и musa64.wmf. Отсюда видим, что точное решение задачи соответствует воздействию musa65.wmf (рис. 2). На рис. 3–6 показано изменение нормального напряжения musa66.wmf во времени musa67.wmf в точках B1–B4: 1 – численное решение; 2 – аналитическое решение.

Для упругих нормальных напряжений musa68.wmf и musa69.wmf имеется хорошее качественное и количественное согласование с результатами точного решения.

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

Выводы

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

Решена задача о распространении плоских продольных волн в виде импульсного воздействия в упругой полуплоскости. Исследуемая расчетная область имеет 62031 узловую точку. Решается система уравнений из 248124 неизвестных. для решения поставленной задачи используется импульсное воздействие (первая ветвь: восходящая часть – четверть круга, нисходящая часть – четверть круга; вторая ветвь: восходящая часть – линейная, нисходящая часть – линейная).

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


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

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

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

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