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

Abstract

Актуальность обусловлена необходимостью предварительного численного моделирования оценки плотности лёсса, уплотняемого глубинными взрывами, с целью снижения экономических затрат при проведении производственных работ. Устранение просадочности грунта реализуется на этапах проектирования и строительства зданий и сооружений для обеспечения их безопасной длительной эксплуатации. Лёссовые просадочные грунты широко распространены в мире, часто в Европе и Азии. Лёссам характерны низкая плотность скелета грунта, высокая пылеватость, макропористость. Цель. Методом конечных разностей численно оценить плотность грунта после уплотнения глубинными взрывами. Методы. Математическое и численное моделирование; метод конечных разностей для оценки плотности грунта; вычислительный эксперимент. Результаты и выводы. Построены дискретные линейные динамические системы на основе применения метода конечных разностей для решения начально-граничных задач в рамках математического моделирования уплотнения просадочных лёссовых грунтов. Основу математической модели образует параболическое дифференциальное уравнение. Для решения использована шеститочечная симметричная схема Кранка-Николсона. Полученные динамические системы учитывают начальное и граничные условия, входные воздействия коэффициента диффузии грунта, мощности заряда взрывчатого вещества, вектор горизонтального распространения газа. Результатом решения систем являются значения плотности грунта в узлах сеточной функции. Оценка плотности грунта, найденная методом конечных разностей, имеет второй порядок точности по пространственным координатам и по времени. Построенные шеститочечные схемы абсолютно устойчивы. Реализован вычислительный эксперимент по предложенной методике оценки плотности грунта, который показал адекватность опытным данным натурного строительного объекта

Relevance. The need for preliminary numerical modeling of the density assessment of loess compacted by deep blasts in order to reduce economic costs during production work. Soil subsidence was eliminated at the stages of design and construction of buildings and structures to ensure their safe, long-term operation. Loess subsidence soils are widespread throughout the world, often in Europe and Asia. Loess is characterized by low density of the soil skeleton, high dust content, and macroporosity. Aim. To numerically estimate the density of soil after compaction by deep blasting using the finite difference method. Methods. Finite difference method; computational experiment; numerical modeling. Results and conclusions. To solve initial boundary problems within the framework of mathematical modeling of loess compaction by deep explosions, the finite difference method was used. The model is based on a parabolic differential equation. A six-point symmetrical Crank-Nicholson circuit was used. The constructed linear systems of equations take into account the initial and boundary conditions, the input effects of the soil diffusion coefficient, the power of the explosive charge, and the vector of horizontal gas distribution. The result of solving the systems are the values of soil density at the nodes of the grid function. The estimate of soil density found by the finite difference method has second order accuracy in spatial coordinates and time. The constructed six-point circuits are absolutely stable. A computational experiment was implemented using the proposed method for assessing soil density, which showed adequacy to the experimental data of a natural construction site