Scientific journal
International Journal of Applied and fundamental research
ISSN 1996-3955
ИФ РИНЦ = 0,564

ON THE RELIABILITY OF COMPUTER SIMULATION OF NON-STATIONARY ELASTIC STRESS WAVES IN DEFORMABLE BODIES OF COMPLEX SHAPE

Musayev V.K. 1
1 MESI
1138 KB
Provides information about the numerical simulation of non-stationary elastic waves in complex deformable areas. Problems are solved using numerical modeling of wave equations of the elasticity theory. Based on the finite element method in the movements of the developed algorithm and software package for solving linear flat two-dimensional problems, which allow to solve complex problems under non-stationary dynamic impacts. When developing complex programs used algorithmic language Fortran-90. The study area is divided by the considered variables on the final elements of the first order. This paper considers the evaluation of the accuracy and reliability of results of numerical modeling of stress waves in regions of complex shape. Comparison with experimental, analytical and numerical methods.
mathematical modeling
dual voltage
round hole
photoproject supported the hole
the wave theory of elasticity dynamic elasticity theory
diffraction
stress concentration
displacement
velocity
displacement
acceleration
seismic impact
the Heaviside function
finite element method
complex programs
anchor point
an explicit two-layer scheme

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

Некоторые вопросы в области постановки, разработки методики, алгоритма и достоверности результатов численного моделирования нестационарных динамических задач рассмотрены в следующих работах [3–10].

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

Для решения двумерной плоской динамической задачи теории упругости с начальными и граничными условиями – используем метод конечных элементов в перемещениях [3–5, 8–10].

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

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

mus02.wmf, mus03.wmf, mus04.wmf, (1)

где: mus05.wmf – матрица инерции; mus06.wmf – матрица жесткости; mus07.wmf – вектор узловых упругих перемещений; mus08.wmf – вектор узловых упругих скоростей перемещений; mus09.wmf – вектор узловых упругих ускорений; mus10.wmf – вектор узловых упругих внешних сил.

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

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

mus11.wmf, mus12.wmf. (2)

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

mus13.wmf,

mus14.wmf. (3)

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

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

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

mus16.wmfmus17.wmf, (4)

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

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

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

Решение методических задач

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

mus1а.tif

Рис. 1. Изменение упругого контурного напряжения mus18.wmf в точке 1 во времени mus19.wmf на контуре свободного круглого отверстия при воздействии плоской продольной упругой волны типа функции Хевисайда

mus2а.tif а)

mus3а.tifб)

Рис. 2. Экспериментальное воздействие mus20.wmf во времени mus21.wmf, полученное методом динамической фотоупругости: а – фотограмма картин полос; б – экспериментальное воздействие, принятое при численном решении методом конечных элементов в перемещениях

mus4а.tif а)

mus5а.tifб)

Рис. 3. Изменение упругого контурного напряжения mus22.wmf в точке 1 во времени mus23.wmf на контуре свободного круглого отверстия при воздействии mus24.wmf: а – фотограмма картин полос; б: 1 – экспериментальные результаты, полученные методом динамической фотоупругости; 2 – результаты численного решения, полученные методом конечных элементов в перемещениях

В сечении на расстоянии 1,9H (рис. 1) при 0 ≤n ≤ 10 (n = t/Δt) скорость упругого перемещения mus25.wmf изменяется линейно от 0 до P (P = σ0/(ρCp) (σ0 = – 0,1 МПа (– 1 кгс/см2)), а при n ≥ 10 mus27.wmf Контур круглого отверстия ABCD предполагается свободным от нагрузок при t > 0. Граничные условия для контура EFGH при t > 0 mus29.wmfОтраженные волны от контура EFGH не доходят до исследуемых точек при 0 ≤ n ≤ 260. Исследуемая расчетная область имеет 1536 узловых точек. Контур круглого отверстия аппроксимирован 28 узловыми точками.

На рис. 1 показано изменение упругого контурного напряжения mus31.wmf(mus32.wmfв точке 1 во времени mus33.wmf(mus34.wmf 1 – результаты аналитического решения [1, 2]; 2 – результаты численного решения, полученные методом конечных элементов в перемещениях [4, 6–8].

Расхождение для максимального упругого контурного напряжения составляет 6 %.

На рис. 2 показано экспериментальное воздействие σ01 во времени mus36.wmf, полученное методом динамической фотоупругости: а – фотограмма картин полос; б – экспериментальное воздействие, принятое при численном решении методом конечных элементов в перемещениях.

На рис. 3 показано изменение упругого контурного напряжения mus37.wmfв точке 1 во времени mus38.wmf при воздействии σ01: а – фотограмма картин полос; б: 1 – экспериментальные результаты, полученные методом динамической фотоупругости [4, 6–8]; 2 – результаты численного решения, полученные методом конечных элементов в перемещениях [4, 6–8].

Расхождение для максимального упругого контурного напряжения составляет 2 %.

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

mus6а.tif

Рис. 4. Изменение упругого контурного напряжения mus39.wmf в точке 1 во времени mus40.wmf на внутреннем контуре подкрепленного круглого отверстия при воздействии плоской продольной упругой волны типа функции Хевисайда

Начальные условия приняты нулевыми. В сечении на расстоянии 1,6H (рис. 4) при 0 ≤ n1 ≤ 10 (n1 = t/Δt1) скорость упругого перемещения mus41.wmf изменяется линейно от 0 до P1 = σ0/(ρ2Cp2), а при n1 ≥ 10 mus42.wmf

Внутренний контур подкрепленного отверстия ABCD предполагается свободным от нагрузок при t > 0. На границе подкрепления и среды EFGH приняты условия непрерывности перемещений. Граничные условия для контура IJKL при t > 0 mus46.wmf Отраженные волны от контура IJKL не доходят до исследуемых точек при 0 ≤ n1 ≤ 540 (mus48.wmf– подкрепление; mus49.wmf – среда). Исследуемая расчетная область имеет 1536 узловых точек. Внутренний контур подкрепления аппроксимирован 28 узловыми точками. По толщине подкрепление аппроксимировано двумя узловыми точками.

На рис. 4 показано изменение контурного напряжения mus50.wmf в точке 1 во времени mus51.wmf (mus52.wmf 1 – результаты аналитического решения [2]; 2 – результаты численного решения, полученные методом конечных элементов в перемещениях [4, 8].

Расхождение для максимального упругого контурного напряжения составляет 12 %.

Выводы

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

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

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