Mathematical Modeling of Non-Stationary Elastic Stress Waves (Transient Process) at Influence (Vertical Concentrated as a Triangular Pulse) on the Surface of a Half-Plane (Lamb’s Problem)
В.К. Мусаев
DOI: https://doi.org/10.46698/vnc.2020.32.87.010
2020-12-21
Геология и геофизика Юга России
Abstract:Рассматривается задача о численном моделировании продольных, поперечных и поверхностных волн на свободной поверхности упругой полуплоскости при воздействии вертикальной сосредоточенной нагрузки в виде треугольного импульса. Полученные результаты исследуемой задачи являются актуальными, так как позволяют выявить типы волн на поверхности упругой полуплоскости, которые применяются в сейсмологии и сейсмостойкости геообъектов. Цель работы. Для оценки несущей способности геообъектов нужна информация о напряженном состоянии. Поэтому получение информации об амплитудах и формах продольных, поперечных и релеевских волн в задаче Лэмба в виде контурных напряжений является актуальной фундаментальной научной задачей. Методика. Для решения нестационарной динамической задачи теории упругости с начальными и граничными условиями используется метод конечных элементов в перемещениях. С помощью метода конечных элементов в перемещениях, линейную задачу с начальными и граничными условиями привели к линейной задаче Коши. Предложен квазирегулярный подход к решению системы линейных обыкновенных дифференциальных уравнений второго порядка в перемещениях с начальными условиями и к аппроксимации исследуемой области. Методика основывается на схемах: точка, линия и плоскость. Исследуемая область разбивается по пространственным переменным на треугольные и прямоугольные конечные элементы первого порядка. По временной переменной исследуемая область разбивается на линейные конечные элементы с двумя узловыми точками. При разработке комплекса программ использовался алгоритмический язык Фортран-90 наблюдается динамический процесс в виде стоячих волн. Результаты. Приводится некоторая информация о численном моделировании упругих волн напряжений в упругой полуплоскости при сосредоточенном волновом воздействии в виде треугольного импульса (дельта-функции). Исследуемая расчетная область имеет 12008001 узловых точек. Решается система уравнений из 48032004 неизвестных. Показано изменение упругого контурного напряжения на свободной поверхности полуплоскости в разных точках. Амплитуда поверхностных волн Релея существенно больше амплитуд продольных, поперечных и других волн при сосредоточенном вертикальном воздействии в виде треугольного импульса на поверхности упругой полуплоскости. После поверхностных волн Релея наблюдается динамический процесс в виде стоячих волн The problem of numerical modeling of longitudinal, transverse and surface waves on the free surface of an elastic half-plane under the influence of a vertical concentrated load in the form of a triangular pulse is considered. The obtained results of the problem under study are relevant, since they allow us to identify the types of waves on the surface of an elastic half-plane that are used in seismology and seismic stability of geo objects. The aim. To assess the load-bearing capacity of geo objects, you need information about the stress state. Therefore, obtaining information about the amplitudes and shapes of longitudinal, transverse, and Rayleigh waves in the lamb problem in the form of contour stresses is an urgent fundamental scientific task. Method. To solve a non-stationary dynamic problem of elasticity theory with initial and boundary conditions, the finite element method in displacements is used. Using the finite element method in displacements, the linear problem with initial and boundary conditions was led to the linear Cauchy problem. A quasi-regular approach to solving a system of second-order linear ordinary differential equations in displacements with initial conditions and to approximating the domain under study is proposed. The method is based on the following diagrams: point, line, and plane. The study area is divided by spatial variables into triangular and rectangular finite elements of the first order. According to the time variable, the study area is divided into linear finite elements with two nodal points. The Fortran-90 algorithmic language was used in the development of the software package. Results. Some information is provided on numerical modeling of elastic stress waves in an elastic half-plane under concentrated wave action in the form of a triangular pulse (Delta-function). The estimated area under study has 12008001 nodal points. A system of equations consisting of 48032004 unknowns is solved. The change in the elastic contour stress on the free surface of the half-plane at different points is shown. The amplitude of surface Rayleigh waves is significantly greater than the amplitudes of longitudinal, transverse, and other waves when a concentrated vertical action is performed in the form of a triangular pulse on the surface of an elastic half-plane. After surface Rayleigh waves, a dynamic process is observed in the form of standing waves