Приближенная оценка глубины проникания сваи в грунт при многократном ударе

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

Рубрика Математика
Вид статья
Язык русский
Дата добавления 26.04.2019
Размер файла 89,4 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

В. Н. Аптуков, А. В. Фонарев

Размещено на http://www.allbest.ru/

46

ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА

2010 Математика. Механика. Информатика Вып.2(2)

41

Пермский государственный университет

Институт механики сплошных сред УрО РАН

Приближенная оценка глубины проникания сваи в грунт при многократном ударе

В.Н. Аптуков, А.В. Фонарев

Аннотация

моделирование грунт численный лагранжевый

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

Ключевые слова: грунт; свая; проникание; методика.

Annotation

The approached estimation of penetration depth of the pile in a ground under multiply impact

V. N. Aptukov, A. V. Fonarev

From the computational point of view, the statement and solution of the problem on rigid body immersion into the ground under normal impact conditions present essential difficulties. In this paper, we develop an approximation approach to the model describing the processes of impact consolidation of foundation pits and immersion of piles. The data of standard tests are used as mathematical model parameters. A comparison is made between the results of calculation obtained by the proposed approximation method and the results of numerical modeling calculated using an explicit finite difference scheme. The results obtained in this study can be used to evaluate the parameters of equipment for impact consolidation of foundation pits and immersion of piles.

Key words: ground; pile; penetration; method.

Динамическое проникание твердого тела в грунт

Уравнение динамического проникания твердого тела с осью симметрии в грунт [1, 2] имеет вид

, (1)

где - масса тела; - скорость; - сила сопротивления прониканию.

Рис. 1 Схема удельных усилий сопротивления грунта, фиксированная и локальная системы координат

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

Тогда сила сопротивления примет вид

, (2)

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

Выражение (2) преобразуем к виду

. (3)

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

Простейшая модель проникания тела в грунт

Максимально простой вариант выражения (3) заключается в следующем: пренебрегаем трением , а сопротивление грунта прониканию полагаем постоянным: .

Тогда уравнение (1) преобразуется к виду

, (4)

где - площадь миделевого сечения (поперечного сечения тела на уровне поверхности грунта).

В данном варианте текущая сила сопротивления грунта определяется средним давлением и площадью миделя, и она не зависит от формы поверхности тела, уже внедренного в грунт.

Представив ускорение в виде , уравнение (4) можно проинтегрировать:

, (5)

где - текущий объем кратера (котлована) (рис.2); - начальная скорость удара.

Пусть тело достигло глубины и остановилось: , тогда

, (6)

где - работа, совершенная силой сопротивления грунта при образовании кратера (котлована) глубиной L.

Соотношение (6) представляет собой энергетическое равенство начальной кинетической энергии тела работе по образованию кратера. В этом смысле удельная величина сопротивления грунта прониканию является энергией, затраченной на образование единицы объема кратера .

Очевидно, что соотношение (6) дает завышенные значения глубины проникания, так как при оценке сопротивления грунта не учитывались трение и инерционные силы.

Рис. 2 Схема проникания тела при постоянном сопротивлении грунта без учета трения

Учет динамического сопротивления

Более точный подход к решению задачи проникания твердого тела в грунт должен учитывать трение и более сложный закон сопротивления . Представим удельное сопротивление прониканию в виде суммы собственного (статического) и инерционного (динамического) сопротивления [3]:

, (7)

где - коэффициент формы; - плотность грунта.

В этом случае уравнение движения примет вид

(8)

Уравнение (8) интегрируется в конечном виде для тел с простой геометрией, например для конуса или цилиндра с конической головкой. Если рассматривается задача об ударном вытрамбовывании котлована, то трением можно пренебречь, что еще больше упрощает анализ.

Задача о забивке сваи в грунт при многократном ударе

Рассмотрим задачу о забивке дизель-молотом длинной сваи в грунт при многократном ударе (рис.3). Задается масса сваи m и молота M, при n-м ударе задается высота падения молота . В ходе решения требуется определить общую глубину погружения сваи и погружение ее за один удар .

Рис. 3 Схема динамической забивки длинной сваи (параметры - слева) и разностная сетка для расчета процесса проникания численным методом [4, 5] (справа)

Скорость падения молота в момент соприкосновения со сваей без учета сил трения в системе при n-м ударе определяется по высоте его падения:

. (9)

Из закона сохранения импульса (предполагается, что коэффициент восстановления удара равен нулю - неупругий удар) следует начальная скорость внедрения сваи при n-м ударе составляет

. (10)

Будем полагать удельное усилие трения постоянным, , а головную часть сваи - конической с углом . Тогда закон движения сваи (8) можно записать в форме линейного неоднородного дифференциального уравнения относительно квадрата скорости:

, (11)

, (12)

где - радиус цилиндрической части сваи.

Решение уравнения (11) с начальными условиями имеет вид

(13)

откуда с учетом условия остановки сваи можно получить соотношение для определения положения сваи после n-го удара:

(14)

Соотношение (14) не позволяет получить конечное выражение для , поэтому упростим уравнение (11), учитывая влияние инерционной составляющей сопротивления в среднем,

. (15)

Уравнение (15) легко интегрируется, что позволяет получить достаточно простое и физически ясное решение для приращения глубины проникания сваи при n-м ударе:

(16)

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

Очевидно, что для одинаковых начальных скоростей сваи с ростом глубины ее проникания (увеличением количества ударов n) приращение глубины проникания при n-м ударе падает по гиперболическому закону.

Примеры расчета

Рассмотрим некоторые численные примеры расчета по предлагаемой методике. На рис. 4 (группа кривых 2) показаны зависимости общей осадки сваи от количества ударов для следующих параметров задачи: = 10 см; = 30°; = 30 кН; = 0,5 МПа; = 2,0 г/см3.

Там же показаны зависимости осадки сваи за один удар от номера удара (группа кривых 1).

Полученные зависимости качественно хорошо согласуются с результатами расчетов по программному комплексу [4, 5].

Рис. 4 Зависимость параметров осадки сваи от количества ударов: = 0,05 МПа (квадрат); = 0,1 МПа (треугольник)

Рис. 5 Зависимость общей осадки сваи от количества ударов: численное решение двухмерной динамической задачи (квадрат); приближенная методика (линия)

Постановка задачи численного моделирования в рамках двумерной нестационарной упругопластической контактной задачи приведена в [4], численная схема реализована в виде пакета прикладных программ PILE GROUND [5].

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

Параметры расчета по приближенной методике: вес сваи 20 кН; вес копра 10 кН; диаметр сваи 40 см; удельный вес грунта 1,8 гс/см3; начальная скорость удара 3 м/с; удельное сопротивление грунта (статическое) = 0,17 МПа; удельное трение = 0,065 МПа.

Для расчетов по программе PILE GROUND были приняты следующие параметры грунта: сцепление = 2 кгс/см2; угол внутреннего трения = 150; коэффициент Пуассона = 0,33; начальный объемный модуль = 59 кгс/см2.

Заключение

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

При глубинах более 3-4 м величина осадки сваи за один удар сильно уменьшается и не превышает 5-10 см. Показано значительное влияние трения на затраты энергии при ударной забивке длинных свай на глубину 5-10 м, когда количество ударов может достигать нескольких сотен.

Список литературы

Аптуков В.Н., Мурзакаев Р.Т., Фонарев А.В. Прикладная теория проникания. М.: Наука. 1992. 104 c.

Аптуков В.Н., Бартоломей А.А., Фонарев А.В., Ирундин С.В. Моделирование процесса ударного вытрамбовывания котлованов // Основания, фундаменты и механика грунтов. 2000. № 3. С.11-15.

Баллистические установки и их применение в экспериментальных исследованиях / под ред. Н.А.Златина и Г.И.Мишина. М.: Наука. 1974. 344 с.

Аптуков В.Н., Фонарев А.В. Численное моделирование процессов ударного и взрывного деформирования элементов конструкций и грунтов: учеб. пос. / Перм. гос. ун-т. Пермь, 2009. 221 с.

Бартоломей А.А. Омельчак И.М., Фонарев А.В. Математическое моделирование динамики погружения свай // Тр. междунар. конф. по проблемам свайного фундаментостроения. М., 1989. С.28-36.

Размещено на Allbest.ru

...

Подобные документы

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

    презентация [240,6 K], добавлен 18.04.2013

  • Задачи для обыкновенных дифференциальных уравнений. Квадратурные формулы. Теоретические основы метода сеток для решения задачи Коши. Погрешность аппроксимации, устойчивость, основная теорема метода сеток. Схема предиктор-корректор 2-го порядка.

    реферат [47,4 K], добавлен 07.12.2013

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

    курсовая работа [120,8 K], добавлен 27.01.2014

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

    контрольная работа [550,9 K], добавлен 12.12.2013

  • Общая характеристика параболических дифференциальных уравнений на примере уравнения теплопроводности. Основные определения и конечно-разностные схемы. Решение дифференциальных уравнений параболического типа методом сеток или методом конечных разностей.

    контрольная работа [835,6 K], добавлен 27.04.2011

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

    курсовая работа [781,4 K], добавлен 28.12.2012

  • Основные положения теории математического моделирования. Структура математической модели. Линейные и нелинейные деформационные процессы в твердых телах. Методика исследования математической модели сваи сложной конфигурации методом конечных элементов.

    курсовая работа [997,2 K], добавлен 21.01.2014

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

    дипломная работа [5,7 M], добавлен 25.05.2014

  • Дифференциальное уравнение первого порядка, разрешенное относительно производной. Применение рекуррентного соотношения. Техника применения метода Эйлера для численного решения уравнения первого порядка. Численные методы, пригодные для решения задачи Коши.

    реферат [183,1 K], добавлен 24.08.2015

  • Использование метода конечных разностей для решения краевой задачи уравнений с частными производными эллиптического типа. Графическое определение распространения тепла методом конечно-разностных аппроксимаций производных с применением пакета Mathlab.

    курсовая работа [1,0 M], добавлен 06.07.2011

  • Общая постановка задачи решения обыкновенных дифференциальных уравнений, особенности использования метода Адамса в данном процессе. Решение системы обыкновенных дифференциальных уравнений методом Адамса и точным методом, сравнение полученных результатов.

    курсовая работа [673,6 K], добавлен 27.04.2011

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

    дипломная работа [460,8 K], добавлен 04.05.2011

  • Особенности решения линейных и нелинейных уравнений. Характеристика и практическое применение и различных методов при решении уравнений. Сущность многочлена Лагранжа и обратного интерполирования. Сравнение численного дифференцирования и интегрирования.

    курсовая работа [799,6 K], добавлен 20.01.2010

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

    статья [41,4 K], добавлен 17.10.2012

  • Понятие и оценка необходимости в статистической обработке психологических данных. Методика и основные этапы математической обработки полученных данных, его критерии и параметры: признаки и переменные, шкалы измерения, анализ и оценка уровня значимости.

    презентация [443,1 K], добавлен 28.02.2014

  • Последовательность решения линейной краевой задачи. Особенности метода прогонки. Алгоритм метода конечных разностей: построение сетки в заданной области, замена дифференциального оператора. Решение СЛАУ методом Гаусса, конечно-разностные уравнения.

    контрольная работа [366,5 K], добавлен 28.07.2013

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

    курсовая работа [596,2 K], добавлен 27.04.2011

  • Исследование числовых рядов на сходимость. Область сходимости для разных степенных рядов. Разложение функции в ряд Тейлора. Нормы сеточной функции. Исследование устойчивости разностной схемы для однородного уравнения. Совокупность разностных уравнений.

    курсовая работа [586,9 K], добавлен 19.04.2011

  • Многошаговые методы и их построение. Вычисление интеграла. Формула для определения неизвестного значения сеточной функции. Запись разностной схемы четвертого порядка. Сущность методов Адамса, Милна, прогноза и коррекции. Оценка точности вычислений.

    презентация [162,9 K], добавлен 18.04.2013

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

    диссертация [2,8 M], добавлен 19.06.2015

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.