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

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

Мусаев В.К. 1
1 Московский государственный машиностроительный университет (МГМУ)
Приводится некоторая информация моделирования безопасности некоторых задач при нестационарных волновых сейсмических воздействиях с помощью метода конечных элементов. Применяются технические средства в виде вертикальных полостей для увеличения безопасности объекта. Для решения поставленных задач применяются волновые уравнения механики деформируемого твердого тела. На основе метода конечных элементов в перемещениях разработаны методика, алгоритм и комплекс программ для решения линейных динамических задач теории упругости. Основные соотношения метода конечных элементов получены с помощью принципа возможных перемещений. Линейная динамическая задача с начальными и граничными условиями с помощью метода конечных элементов в перемещениях приведена к системе линейных обыкновенных дифференциальных уравнений с начальными условиями. Задача с начальными условиями с помощью конечноэлементного варианта метода Галеркина приведена к явной двухслойной схеме. Рассмотрена постановка некоторых задачи с полостью и без полости в полуплоскости при воздействии в виде функции Хевисайда. Исследуемая расчетная область имеет 14762 узловых точек и 14520 конечных элементов. Решается система уравнений из 59048 неизвестных.
моделирование
математическое моделирование
численный метод
алгоритм
комплекс программ
метод Мусаева В.К.
нестационарные упругие волны
динамика сплошных сред
волновая теория сейсмической безопасности
конечноэлементный метод Галеркина
распространение волн
вертикальные прямоугольные полости
исследуемая расчетная область
функция Хевисайда
1. Мусаев В.К. Определение качества сооружений в детерминированной постановке с помощью математического мониторинга // Вестник Российского университета дружбы народов. Серия проблемы комплексной безопасности. – 2005. – № 1. – С. 42–47.
2. Мусаев В.К. Об управлении упругим волновым напряженным состоянием сооружений с помощью полостей в виде прямоугольников при сейсмических воздействиях // Вестник Российского университета дружбы народов. Серия проблемы комплексной безопасности. – 2007. – № 2. – С. 6–17.
3. Мусаев В.К. Численное моделирование упругих сейсмических волн напряжений в сложных деформируемых телах // Вестник Российского университета дружбы народов. Серия проблемы комплексной безопасности. – 2007. – № 4. – С. 6–22.
4. Мусаев В.К. Математическое моделирование сейсмической безопасности сооружений неглубокого заложения с помощью вертикальных полостей // Вестник Российского университета дружбы народов. Серия проблемы комплексной безопасности. – 2008. – № 1. – С. 6–14.
5. Мусаев В.К. Об устойчивости одномерной явной двухслойной конечноэлементной линейной схемы в перемещениях для внутренних узловых точек на равномерной линейной сетке // Вестник Российского университета дружбы народов. Серия проблемы комплексной безопасности. – 2008. – № 1. – С. 57–60.
6. Мусаев В.К. О численном моделировании безопасности системы «сооружение-грунт-воздух» при сейсмических воздействиях // Двойные технологии. – 2013. – № 3. – С. 47–49.
7. Мусаев В.К. О достоверности результатов математического моделирования нестационарных волн напряжений в объектах сложной формы // Строительная механика инженерных конструкций и сооружений. – 2014. – № 3. – С. 71–76.
8. Мусаев В.К. О достоверности компьютерного моделирования нестационарных упругих волн напряжений в деформируемых телах сложной формы // Международный журнал прикладных и фундаментальных исследований. – 2014. – № 11. – С. 10–14.
9. Мусаев В.К. Численное моделирование нестационарных упругих волн напряжений в некоторых задачах методического характера // Международный журнал экспериментального образования. – 2015. – № 11 (часть 2). – С. 227–230.
10. Мусаев В.К. Моделирование нестационарных процессов в геообъектах с помощью волновой теории сейсмической безопасности // Международный журнал экспериментального образования. – 2015. – № 12. – С. 347–352.

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

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

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

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

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

mus_201.wmf,

mus_202.wmf, mus_203.wmf,

mus_204.wmf,

mus_205.wmf, mus_206.wmf,

mus_207.wmf, mus_208.wmf, mus_209.wmf,

mus_210.wmf, (1)

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

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

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

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

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

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

mus_215.wmf, mus_216.wmf, mus_217.wmf, (2)

где mus_218.wmf – диагональная матрица инерции; mus_219.wmf – матрица жесткости; mus_220.wmf – вектор узловых упругих перемещений; mus_221.wmf – вектор узловых упругих скоростей перемещений; mus_222.wmf – вектор узловых упругих ускорений; mus_223.wmf – вектор внешних узловых упругих сил.

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

mus_224.wmf,

mus_225.wmf. (3)

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

mus_226.wmf mus_227.wmf, (4)

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

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

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

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

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

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

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

musaev1.tif

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

От точки B параллельно свободной поверхности ABC приложено нормальное напряжение σx, которое при 0 ≤ n ≤ 10 (mus_228.wmf) изменяется линейно от 0 до P, а при n ≥ 10 равно P (P = σ0, σ0 = 0,1 МПа (1 кгс/см2)). Граничные условия для контура CDEA при t > 0 mus_229.wmf. Отраженные волны от контура CDEA не доходят до исследуемых точек при 0 ≤ n ≤ 200. Контур ABC свободен от нагрузок, кроме точки B.

musaev2.tif

Рис. 2. Воздействие типа функции Хевисайда

Расчеты проведены при следующих исходных данных: mus_230.wmf; ?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 м/с. Исследуемая расчетная область имеет 14762 узловых точек и 14520 конечных элементов. Решается система уравнений из 59048 неизвестных. В характерных областях исследуемой задачи получены контурные напряжения и компоненты тензора напряжений.

2. Рассмотрим задачу о воздействии плоской продольной сейсмической волны (рис. 2) параллельной свободной поверхности упругой полуплоскости с полостью (соотношение ширины к высоте один к четырем) (рис. 3). От точки F параллельно свободной поверхности ABEFG приложено нормальное напряжение σx, которое при 0 ≤ n ≤ 10 (mus_231.wmf) изменяется линейно от 0 до P, а при n ≥ 10 равно P (P = σ0, σ0 = 0,1 МПа (1 кгс/см2)).

Граничные условия для контура GHIA при t > 0 mus_232.wmf. Отраженные волны от контура GHIA не доходят до исследуемых точек при 0 ≤ n ≤ 200. Контур ABCDEFG свободен от нагрузок, кроме точки F. Расчеты проведены при следующих исходных данных: mus_233.wmf; ?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 м/с. Исследуемая расчетная область имеет 14762 узловых точек и 14516 конечных элементов. Решается система уравнений из 59048 неизвестных. В характерных областях исследуемой задачи получены контурные напряжения и компоненты тензора напряжений.

3. Рассмотрим задачу о воздействии плоской продольной сейсмической волны (рис. 2) параллельной свободной поверхности упругой полуплоскости с полостью (соотношение ширины к высоте один к восьми) (рис. 4). От точки F параллельно свободной поверхности ABEFG приложено нормальное напряжение σx, которое при 0 ≤ n ≤ 10 (mus_234.wmf) изменяется линейно от 0 до P, а при n ≥ 10 равно P (P = σ0, σ0 = 0,1 МПа (1 кгс/см2)). Граничные условия для контура GHIA при t > 0 mus_235.wmf. Отраженные волны от контура GHIA не доходят до исследуемых точек при 0 ≤ n ≤ 200. Контур ABCDEFG свободен от нагрузок, кроме точки F.

musaev3.tif

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

musaev4.tif

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

Расчеты проведены при следующих исходных данных: mus_236.wmf; ?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 м/с. Исследуемая расчетная область имеет 14762 узловых точек и 14512 конечных элементов. Решается система уравнений из 59048 неизвестных.

4. Рассмотрим задачу о воздействии плоской продольной сейсмической волны (рис. 2) параллельной свободной поверхности упругой полуплоскости с полостью (соотношение ширины к высоте один к двенадцати) (рис. 5). От точки F параллельно свободной поверхности ABEFG приложено нормальное напряжение σx, которое при 0 ≤ n ≤ 10 (mus_237.wmf) изменяется линейно от 0 до P, а при n ≥ 10 равно P (P = σ0, σ0 = 0,1 МПа (1 кгс/см2)). Граничные условия для контура GHIA при t > 0 mus_238.wmf. Отраженные волны от контура GHIA не доходят до исследуемых точек при 0 ≤ n ≤ 200. Контур ABCDEFG свободен от нагрузок, кроме точки F. Расчеты проведены при следующих исходных данных: mus_239.wmf; ?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 м/с. Исследуемая расчетная область имеет 14762 узловых точек и 14508 конечных элементов. Решается система уравнений из 59048 неизвестных.

musaev5.tif

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

Выводы

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


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

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

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

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