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

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

Мусаев В.К. 1
1 МГМУ
Приводится некоторая информация моделирования упругих волн напряжений в бесконечной полосе при воздействии в виде треугольного импульса. Задача решается с помощью численного моделирования нестационарных динамических уравнений теории упругости. Отраженные волны от свободных поверхностей бесконечной полосы создают физическую картину стоячих волн. Для решения двумерной нестационарной динамической задачи математической теории упругости с начальными и граничными условиями используем метод конечных элементов в перемещениях. Задача решается методом сквозного счета, без выделения разрывов. Применяется однородный алгоритм. С помощью метода конечных элементов в перемещениях, линейную задачу с начальными и граничными условиями привели к линейной задаче Коши. Показаны компоненты нормальных напряжений в характерной области исследуемой задачи.
нестационарные волны
численный метод
перемещение
скорость перемещений
ускорение
напряжение
теория упругости
краевая задача
задача с начальными условиями
задача Коши
методика
алгоритм
комплекс программ
однородный алгоритм
импульсное воздействие
треугольный импульс
бесконечная полоса
стоячая волна
1. Мусаев В.К. Численное моделирование динамического напряженного состояния сооружений уравнениями двумерной теории упругости и пластичности. Автореферат диссертации на соискание ученой степени доктора технических наук по специальности 01.02.04. – М.: Совинтервод, 1993. – 46 с.
2. Мусаев В.К. Численное решение волновых задач теории упругости и пластичности // Вестник Российского университета дружбы народов. Серия прикладная математика и информатика. – 1997. – № 1. – С. 87–110.
3. Мусаев В.К. Математическое моделирование упругих волн напряжений в сложных деформируемых телах // Вестник Российского университета дружбы народов. Серия проблемы комплексной безопасности. – 2007. – № 1. – С. 62–76.
4. Мусаев В.К. Об оценке достоверности и точности численного решения нестационарных динамических задач // Вестник Российского университета дружбы народов. Серия проблемы комплексной безопасности. – 2007. – № 3. – С. 48–60.
5. Мусаев В.К. Оценка достоверности и точности результатов вычислительного эксперимента при решении задач нестационарной волновой теории упругости // Научный журнал проблем комплексной безопасности. – 2009. – № 1. – С. 55–80.
6. Мусаев В.К. О достоверности компьютерного моделирования нестационарных упругих волн напряжений в деформируемых телах сложной формы // Международный журнал прикладных и фундаментальных исследований. – 2014. – № 11. – С. 10–14.
7. Мусаев В.К. Моделирование нестационарных упругих волн напряжений в деформируемых областях с помощью метода конечных элементов в перемещениях // Современные наукоемкие технологии. – 2014. – № 12 (1). – С. 28–32.
8. Мусаев В.К. Оценка точности и достоверности численного моделирования при решении задач об отражении и интерференции нестационарных упругих волн напряжений // Успехи современного естествознания. – 2015. – № 1 (часть 7). – С. 1184–1187.
9. Musayev V.K. Modeling of non-stationary of stress waves in solid deformable bodies complex area // International Journal Of Applied And Fundamental Research. – 2014. – № 2; URL: www.science-sd.com/457-24639.
10. Musayev V.K. Estimation of accuracy of the results of numerical simulation of unsteady wave of the stress in deformable objects of complex shape // International Journal for Computational Civil and Structural Engineering. – 2015. – Volume 11, Issue 1. – P. 135–146.

О численном методе, алгоритме и комплексе программ моделирования волн напряжений

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

При динамическом и импульсном воздействии в сооружении распространяются волны напряжений.

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

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

mus01.wmf, mus02.wmf,

mus03.wmf,

mus04.wmf,

mus05.wmf, mus06.wmf,

mus07.wmf, mus08.wmf,

mus09.wmf, mus10.wmf, (1)

где σx, σy и τxy – компоненты тензора упругих напряжений;

εx, εy и γxy – компоненты тензора упругих деформаций;

u и v – составляющие вектора упругих перемещений вдоль осей ОХ и OУ соответственно;

ρ – плотность материала;

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

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

v – коэффициент Пуассона;

Е – модуль упругости;

mus13.wmf – граничный контур тела Г.

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

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

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

mus16.wmf, mus17.wmf, mus18.wmf, (2)

где mus19.wmf – матрица инерции;

mus20.wmf – матрица жесткости;

mus21.wmf – вектор узловых упругих перемещений;

mus22.wmf – вектор узловых упругих скоростей перемещений;

mus23.wmf – вектор узловых упругих ускорений;

mus24.wmf – вектор узловых упругих внешних сил.

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

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

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

mus25.wmf, mus26.wmf. (3)

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

mus27.wmf,

mus28.wmf. (4)

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

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

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

mus29.wmf mus30.wmf, (5)

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

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

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

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

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

musaev1.tif

Рис. 1. Воздействие в виде треугольного импульса

На границе пластинки mus31.wmf (рис. 2) приложено нормальное напряжение σy (рис. 1), которое при 1 ≤ n ≤ 3 (n = t/∆t) изменяется линейно от 0 до P, а при 3 ≤ n ≥ 5 от P до 0 (P = σ0, σ0 = – 0,1 МПа). Граничные условия для контуров BC и AD при t > 0 mus32.wmf. Отраженные волны от контуров BC и AD не доходят до исследуемых точек. Контур CD свободен от нагрузок. Исследуемая расчетная область имеет 4002 узловые точки. Решается система уравнений из 16008 неизвестных.

musaev2.tif

Рис. 2. Постановка задачи о стоячих волнах в бесконечной полосе

musaev3.tif

Рис. 3. Нормальное напряжение mus33.wmf во времени 0 ≤ n ≤ 100 в точке B1

musaev4.tif

Рис. 4. Нормальное напряжение mus34.wmf во времени 0 ≤ n ≤ 500 в точке B1

musaev5.tif

Рис. 5. Нормальное напряжение mus35.wmf во времени 0 ≤ n ≤ 1990 в точке B1

musaev6.tif

Рис. 6. Нормальное напряжение mus36.wmf во времени 0 ≤ n ≤ 100 в точке B1

musaev7.tif

Рис. 7. Нормальное напряжение mus37.wmf во времени 0 ≤ n ≤ 500 в точке B1

musaev8.tif

Рис. 8. Нормальное напряжение mus38.wmf во времени 0 ≤ n ≤ 1990 в точке B1

Для примера на рис. 3–5 представлено изменение нормального напряжения mus39.wmf (mus40.wmf) во времени n в точке B1.

Для примера на рис. 6–8 представлено изменение нормального напряжения mus42.wmf (mus43.wmf) во времени n в точке B1.

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


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

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

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

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