Научный журнал
Международный журнал прикладных и фундаментальных исследований

ISSN 1996-3955
ИФ РИНЦ = 0,686

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

Мусаев В.К. 1
1 МЭСИ
Рассматривается задача о численном моделировании продольных, поперечных и поверхностных волн на свободной поверхности упругой полуплоскости. Приводится изменение упругого контурного напряжения на свободной поверхности полуплоскости. Для решения двумерной нестационарной динамической задачи математической теории упругости с начальными и граничными условиями используем метод конечных элементов в перемещениях. Задача решается методом сквозного счета, без выделения разрывов. Применяется однородный алгоритм. С помощью метода конечных элементов в перемещениях, линейную задачу с начальными и граничными условиями привели к линейной задаче Коши. Приводится некоторая информация о численном моделировании упругих волн напряжений в упругой полуплоскости при сосредоточенном взрывном воздействии в виде дельта функции. Показано изменение упругого контурного напряжения на свободной поверхности полуплоскости.
нестационарные волны
численный метод
перемещение
скорость перемещений
ускорение
напряжение
теория упругости
краевая задача
задача с начальными условиями
задача Коши
методика
алгоритм
однородный алгоритм
комплекс программ
продольная волна
поперечная волна
коническая волна
волна Релея
поверхностная волна
задача Лэмба
упругая полуплоскость
напряжения на свободной поверхности
1. Мусаев В.К. Решение задачи дифракции и распространения упругих волн методом конечных элементов // Строительная механика и расчет сооружений. – 1990. – № 4. – С. 74–78.
2. Мусаев В.К. Численное решение волновых задач теории упругости и пластичности // Вестник Российского университета дружбы народов. Серия прикладная математика и информатика. – 1997. – № 1. – С. 87–110.
3. Мусаев В.К. Об оценке достоверности и точности численного решения нестационарных динамических задач // Вестник Российского университета дружбы народов. Серия проблемы комплексной безопасности. – 2007. – № 3. – С. 48–60.
4. Мусаев В.К. Оценка достоверности и точности результатов вычислительного эксперимента при решении задач нестационарной волновой теории упругости // Научный журнал проблем комплексной безопасности. – 2009. – № 1. – С. 55–80.
5. Мусаев В.К. О моделировании отражения упругих волн напряжений от свободной поверхности деформируемой области // Двойные технологии. – 2012. – № 4. – С. 61–64.
6. Мусаев В.К. О моделировании безопасности технических объектов от взрывных воздействий // Стратегическая стабильность. – 2013. – № 1. – С. 69–72.
7. Мусаев В.К. Математическое моделирование волн напряжений в задаче Лэмба // Информационно-телекоммуникационные технологии и математическое моделирование высокотехнологичных систем. Материалы Всероссийской конференции с международным участием. – М.: РУДН, 2014. – С. 266–269.
8. Мусаев В.К. Математическое моделирование интерференции нестационарных упругих волн напряжений в виде треугольного импульса от свободной поверхности пластинки / В.К. Мусаев, С.В. Ситник, А.А. Тарасенко, В.Г. Ситник, М.В. Зюбина // Современные проблемы науки и образования. – 2014. – № 4; URL: www.science-education.ru/118-14118 (дата обращения: 21.09.2014).
9. Мусаев В.К. Математическое моделирование отражения нестационарных упругих волн напряжений в виде треугольного импульса от свободной поверхности пластинки / В.К. Мусаев, С.В. Ситник, А.А. Тарасенко, В.Г. Ситник, М.В. Зюбина // Фундаментальные исследования. – 2014. – № 9 (часть 7). – С. 1466–1470; URL: www.rae.ru/fs/?section=content&op=show_article&article_id=10004353 (дата обращения: 21.09.2014).
10. Мусаев В.К. О достоверности компьютерного моделирования нестационарных упругих волн напряжений в деформируемых телах сложной формы // Международный журнал прикладных и фундаментальных исследований. – 2014. – № 11 – С. 10–14.

Постановка задачи

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

usaev1.tif

Рис. 1. Некоторое тело Г в прямоугольной декартовой системе координат XOY

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

musaev2.wmf, musaev3.wmf,

musaev4.wmf,

musaev5.wmf,

musaev6.wmf,

musaev7.wmf,

musaev8.wmf, musaev9.wmf, musaev10.wmf,

musaev11.wmf, (1)

где musaev12.wmf, musaev13.wmf и musaev14.wmf – компоненты тензора упругих напряжений; musaev15.wmf, musaev16.wmf и musaev17.wmf – компоненты тензора упругих деформаций; u и v – cоставляющие вектора упругих перемещений вдоль осей OX и OY соответственно; r – плотность материала; musaev18.wmf – скорость продольной упругой волны; musaev19.wmf – скорость поперечной упругой волны; n – коэффициент Пуассона; E – модуль упругости; musaev20.wmf – граничный контур тела Г.

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

Начальные условия в области Г зададим в виде

musaev21.wmf, musaev22.wmf, musaev23.wmf,

musaev24.wmf, musaev25.wmf, (2)

где musaev26.wmf, musaev27.wmf, musaev28.wmf и musaev29.wmf – заданные в области Г функции.

Граничные условия зададим в виде:

составляющих компонентов тензора упругих напряжений на границе S1

musaev30.wmf, musaev31.wmf,

musaev32.wmf; (3)

составляющих компонентов вектора упругих перемещений на границе S2

musaev33.wmf, musaev34.wmf, musaev35.wmf, (4)

где l и m – направляющие косинусы; musaev36.wmf, musaev37.wmf, musaev38.wmf и musaev39.wmf – заданные на границе S функции.

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

Метод решения двумерной плоской динамической задачи теории упругости

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

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

musaev40.wmf, musaev41.wmf, musaev42.wmf, (5)

где musaev43.wmf – матрица инерции; musaev44.wmf – матрица жесткости; musaev45.wmf – вектор узловых упругих перемещений; musaev46.wmf – вектор узловых упругих скоростей перемещений; musaev47.wmf – вектор узловых упругих ускорений; musaev48.wmf – вектор узловых упругих внешних сил.

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

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

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

musaev49.wmf, musaev50.wmf. (6)

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

musaev51.wmf,

musaev52.wmf, (7)

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

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

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

musaev56.wmf, musaev57.wmf, (8)

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

Результаты численного эксперимента показали, что при k = 0,5 обеспечивается устойчивость двумерной явной двухслойной конечноэлементной линейной схемы в перемещениях для внутренних и граничных узловых точек на квазирегулярных сетках.

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

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

Решение задачи о сосредоточенном упругом взрывном воздействии

В упругой полуплоскости от сосредоточенного воздействия распространяются продольные, поперечные, рэлеевские и конические волны. Они распространяются с разной скоростью. В работе [7] приведена информация о численном решении задачи Лэмба. Рассмотрим задачу о воздействии сосредоточенной волны в виде дельта функции (рис. 3) перпендикулярной свободной поверхности упругой полуплоскости (рис. 2). В точке B перпендикулярно свободной поверхности musaev59.wmf приложено упругое нормальное напряжение musaev60.wmf (рис. 3), которое при musaev61.wmf musaev62.wmf изменяется линейно от 0 до P, а при musaev63.wmf от P до 0 (musaev64.wmf, musaev65.wmf МПа (–1 кгс/см2)).

usaev2.tif

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

usaev3.tif

Рис. 3. Воздействие типа дельта функции

usaev4.tif

Рис. 4. Изменение упругого контурного напряжения musaev66.wmf во времени musaev67.wmf в точке A1

usaev5.tif

Рис. 5. Изменение упругого контурного напряжения musaev68.wmf во времени musaev69.wmf в точке A2

usaev6.tif

Рис. 6. Изменение упругого контурного напряжения musaev70.wmf во времени musaev71.wmf в точке A3

Граничные условия для контура CDEA при musaev72.wmf musaev73.wmf. Отраженные волны от контура CDEA не доходят до исследуемых точек при musaev74.wmf. Контур ABC свободен от нагрузок, кроме точки B, где приложено сосредоточенное упругое нормальное напряжение musaev75.wmf.

Расчеты проведены при следующих исходных данных: musaev76.wmf; musaev77.wmf = 1,393⋅10-6 с; E = 3,15⋅10 4 МПа (3,15⋅10 5 кгс/см2); musaev78.wmf= 0,2; musaev79.wmf= 0,255⋅104 кг/м3 (0,255⋅10-5 кгс с2/см4); musaev80.wmf= 3587 м/с; musaev81.wmf= 2269 м/с. Решается система уравнений из 48032004 неизвестных.

На рис. 4–6 показано изменение упругого контурного напряжения musaev82.wmf (musaev83.wmf) во времени n в точках A1–A3 (рис. 2), находящихся на свободной поверхности упругой полуплоскости (расстояние между точками: A1 и A2 равно H; A2 и A3 равно H).

Вывод

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


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

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

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

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