Document 123658

ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА
2010
М атем ати ка. М е хан и ка. Ин форма тик а
Вып.2(2)
УДК 539.3
Приближенная оценка глубины проникания
сваи в грунт при многократном ударе
В. Н. Аптуков 1 , А. В. Фонарев 2
Пермский государственный университет, Россия, 614990, Пермь, ул. Букирева, 15
aptukov@psu.ru; (342) 239-68-19
2
Институт механики сплошных сред УрО РАН, Россия, 614013, Пермь, ул. Академика Королева, 1
afonarev@icmm.ru; (342) 237-83-08
1
При полной постановке и решении задачи проникания недеформируемого тела при ударе по
нормали в грунт возникают значительные вычислительные затруднения. Статья посвящена
приближенному подходу к описанию этого процесса. В качестве параметров математической
модели используются данные стандартных испытаний грунтов. Приведено сравнение расчетных данных, полученных по приближенной методике, с результатами численного моделирования на основе явной лагранжевой конечно-разностной схемы. Полученные результаты
могут быть использованы при оценивании параметров оборудования, применяемого при
ударном вытрамбовывании котлованов и погружении свай.
Ключевые слова: грунт; свая; проникание; методика.
Пусть на каждую точку тела, соприкасающуюся с грунтом, действует удельное усилие
 c , направленное по нормали, и удельное уси-
1. Динамическое проникание
твердого тела в грунт
Уравнение динамического проникания
твердого тела с осью симметрии в грунт [1, 2]
имеет вид
m
d
 Fc ,
dt
лие  c , направленное по касательной к поверхности в этой точке (рис. 1). Введем две системы
координат: фиксированную Ox , связанную с
поверхностью грунта, и локальную rO1 z , связанную с движущимся телом.
Тогда сила сопротивления примет вид
(1)
где m – масса тела;  – скорость; Fc – сила
сопротивления прониканию.
Fc    (  c cos    c sin  )dS ,
(2)
S
где S – площадь поверхности тела, соприкасающаяся с грунтом; tg  f ( z ) ; f ( z ) – переменный радиус образующей тела.
Выражение (2) преобразуем к виду
x( t )
Fc  2
 (
c
  c f  ) fdz .
(3)
0
Рис. 1. Схема удельных усилий сопротивления грунта, фиксированная и локальная системы координат

Отметим, что уравнение (1) с правой
частью (3) является практически точным,
здесь не было сделано пока никаких предпо-
© В. Н. Аптуков, А. В. Фонарев, 2010
41
В. Н. Аптуков, А. В. Фонарев
ложений. Основной проблемой является задание удельных усилий  c и  c .
Очевидно, что соотношение (6) дает завышенные значения глубины проникания, так
как при оценке сопротивления грунта не учитывались трение и инерционные силы.
2. Простейшая модель проникания
тела в грунт
Максимально простой вариант выражения (3) заключается в следующем: пренебрегаем трением  c  0 , а сопротивление грунта
прониканию
полагаем
постоянным:
 c  const .
Тогда уравнение (1) преобразуется к
виду
d
m
 Sm ( x) c  f 2 ( x) c ,
dt
Рис. 2. Схема проникания тела при
постоянном сопротивлении грунта
без учета трения
(4)
где S m – площадь миделевого сечения (поперечного сечения тела на уровне поверхности
грунта).
В данном варианте текущая сила сопротивления грунта определяется средним давлением  c и площадью миделя, и она не зависит от формы поверхности тела, уже внедренного в грунт.
Представив
ускорение
в
виде
d 1 d 2 , уравнение (4) можно проинте
( )
dt 2 dx
грировать:
1
m( 2 ( x)  02 )   cV ( x),
2
,
(5)
x
3. Учет динамического сопротивления
Более точный подход к решению задачи
проникания твердого тела в грунт должен
учитывать трение  c и более сложный закон
сопротивления  c . Представим удельное сопротивление прониканию в виде суммы собственного (статического)  0 и инерционного
(динамического) сопротивления [3]:
(7)
 c   0  k 2 ,
где k  sin 2  – коэффициент формы;  –
плотность грунта.
В этом случае уравнение движения
примет вид
x
1 d 2
m ( )  2  [ c  ( 0  k 2 ) f ] fdz (8)
2 dx
0
Уравнение (8) интегрируется в конечном виде для тел с простой геометрией,
например для конуса или цилиндра с конической головкой. Если рассматривается задача
об ударном вытрамбовывании котлована, то
трением можно пренебречь, что еще больше
упрощает анализ.
V ( x)   S m ( x)dx
0
где V (x) – текущий объем кратера (котлована) (рис.2);  0 – начальная скорость удара.
Пусть тело достигло глубины
остановилось:  ( L)  0 , тогда
1
m02  W ( L)   cV ( L) ,
2
L
и
(6)
где W ( L ) – работа, совершенная силой сопротивления грунта при образовании кратера
(котлована) глубиной L.
Соотношение (6) представляет собой
энергетическое равенство начальной кинетической энергии тела работе по образованию
кратера. В этом смысле удельная величина
сопротивления грунта прониканию  c является энергией, затраченной на образование
единицы объема кратера W / V .
4. Задача о забивке сваи в грунт
при многократном ударе
Рассмотрим задачу о забивке дизель–
молотом длинной сваи в грунт при многократном ударе (рис.3). Задается масса сваи m
и молота M, при n-м ударе задается высота
падения молота H n . В ходе решения требует-
42
Приближенная оценка глубины проникания сваи в грунт …
ся определить общую глубину погружения
сваи
n
k
x
A B B n
 x0  (0n ) 2 ) exp(C ( x  x0n )) 
C C2 C
(13)
A B B
 (  2  x)
C C
C
2  ( 
и погружение ее за один удар Ln .
откуда с учетом условия остановки сваи
 ( x  xkn )  0 можно получить соотношение
n
для определения положения сваи xk после nго удара:
( AC  B  BCx0n  C 2 (0n )2 ) exp(C ( xkn  x0n )) 
(14)
 AC  B  BCxkn
Соотношение (14) не позволяет полуn
чить конечное выражение для xk , поэтому
упростим уравнение (11), учитывая влияние
инерционной составляющей сопротивления в
среднем,
d 2
1
( )  A  Bx  C (0n )2 .
dx
2
Рис. 3. Схема динамической забивки длинной
сваи (параметры – слева) и разностная сетка
для расчета процесса проникания численным
методом [4, 5] (справа)
Уравнение (15) легко интегрируется,
что позволяет получить достаточно простое и
физически ясное решение для приращения
глубины проникания сваи при n-м ударе:
Скорость падения молота в момент соприкосновения со сваей без учета сил трения
в системе при n-м ударе определяется по высоте его падения:
(9)
U 0n  2 gH n .
Ln  xkn  x0n 

Из закона сохранения импульса (предполагается, что коэффициент восстановления
удара равен нулю – неупругий удар) следует
начальная скорость внедрения сваи при n-м
ударе составляет
 0n 
U 0n M .
(m  M )
(15)
(m  M )(0n ) 2
(16)
1
n 2
n
4 R( R( 0  k  (0 ) )  x0 c )
2
При первых ударах для малых значений
( x01  0 ) основное сопротивление оказывает лобовая коническая часть сваи. С ростом
глубины проникания основную роль начинает
играть сопротивление трения на боковой поверхности сваи, которое для длинных свай
может на порядок превышать лобовое сопротивление.
Очевидно, что для одинаковых начальных
скоростей сваи 0n  const с ростом глубины ее
n
0
x
(10)
Будем полагать удельное усилие трения
постоянным,  c  const , а головную часть сваи
– конической с углом 2 . Тогда закон движения сваи (8) можно записать в форме линейного неоднородного дифференциального уравнения относительно квадрата скорости:
n
проникания x0 (увеличением количества ударов
n) приращение глубины проникания Ln при n-м
ударе падает по гиперболическому закону.
d
(11)
(  2 )  A  Bx  C 2 ,
dx
4R 2 0
4R c
4R 2k
, (12)
A
, B
,C
mM
mM
mM
5. Примеры расчета
Рассмотрим некоторые численные примеры расчета по предлагаемой методике. На
рис. 4 (группа кривых 2) показаны зависимости общей осадки сваи от количества ударов
для следующих параметров задачи: R = 10 см;
где R – радиус цилиндрической части сваи.
Решение уравнения (11) с начальными
условиями  ( x  x0n )  0n имеет вид
43
В. Н. Аптуков, А. В. Фонарев
 = 30°; m  M = 30 кН;  0 = 0,5 МПа;
 грунта = 2,0 г/см3.
женной методике с результатами численного
моделирования, требующего больших вычислительных затрат.
Параметры расчета по приближенной
методике: вес сваи 20 кН; вес копра 10 кН;
диаметр сваи 40 см; удельный вес грунта 1,8
гс/см3; начальная скорость удара 3 м/с; удельное сопротивление грунта (статическое)  0 =
Там же показаны зависимости осадки
сваи за один удар от номера удара (группа
кривых 1).
Полученные зависимости качественно
хорошо согласуются с результатами расчетов
по программному комплексу [4, 5].
0,17 МПа; удельное трение  c = 0,065 МПа.
Для расчетов по программе PILE
GROUND были приняты следующие параметры грунта: сцепление c0 = 2 кгс/см2; угол
внутреннего трения  = 150; коэффициент
Пуассона  = 0,33; начальный объемный модуль K0 = 59 кгс/см2.
6. Заключение
Предложена инженерная методика расчета проникания сваи в грунт при многократном ударе, позволяющая с минимальными
затратами оценить влияние основных параметров процесса на динамику ударного погружения.
При глубинах более 3–4 м величина
осадки сваи за один удар сильно уменьшается
и не превышает 5–10 см. Показано значительное влияние трения на затраты энергии при
ударной забивке длинных свай на глубину 5–
10 м, когда количество ударов может достигать нескольких сотен.
Рис. 4. Зависимость параметров осадки сваи от
количества ударов:  c = 0,05 МПа (квадрат);  c
= 0,1 МПа (треугольник)
Список литературы
1. Аптуков В.Н., Мурзакаев Р.Т., Фонарев А.В.
Прикладная теория проникания. М.: Наука.
1992. 104 c.
2. Аптуков В.Н., Бартоломей А.А., Фонарев
А.В., Ирундин С.В. Моделирование процесса
ударного вытрамбовывания котлованов //
Основания, фундаменты и механика грунтов. 2000. № 3. С.11–15.
3. Баллистические установки и их применение
в экспериментальных исследованиях / под
ред. Н.А.Златина и Г.И.Мишина. М.: Наука.
1974. 344 с.
4. Аптуков В.Н., Фонарев А.В. Численное моделирование процессов ударного и взрывного деформирования элементов конструкций и грунтов: учеб. пос. / Перм. гос. ун-т.
Пермь, 2009. 221 с.
Рис. 5. Зависимость общей осадки сваи от количества ударов: численное решение двухмерной
динамической задачи (квадрат); приближенная
методика (линия)
Постановка задачи численного моделирования в рамках двумерной нестационарной
упругопластической контактной задачи приведена в [4], численная схема реализована в
виде пакета прикладных программ PILE
GROUND [5].
На рис. 5 показаны сравнительные результаты расчета по предлагаемой прибли-
44
Приближенная оценка глубины проникания сваи в грунт …
5. Бартоломей А.А. Омельчак И.М., Фонарев А.В. Математическое моделирование динамики погружения свай // Тр.
междунар. конф. по проблемам свайного фундаментостроения. М., 1989.
С.28–36.
The approached estimation of penetration depth
of the pile in a ground under multiply impact
V. N. Aptukov 1 , A. V. Fonarev 2
1
Perm State University, Russia, 614990, Perm, Bukireva st., 15
aptukov@psu.ru; (342) 239-68-19
2
ICMM of Ural Branch RAS, Russia , 614013, Perm, Ac. Koroleva st., 1
afonarev@icmm.ru; (342) 237-83-08
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.
45