ui - сума размеров сторон i-го сечения сваи, м;

u0i - сумма размеров сторон i-го поперечного сечения сваи, которое имеет наклон к оси сваи;

ip - наклон боковых граней сваи в долях единицы;

i £ 0,025

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

Ki - коэффициент, зависящий от вида грунта и принимаемый по таблице 4 СНиПа [18];

xr - реологический коэффициент, принимаемый xr = 0,8.

При расчете пирамидальных свай по СНиП, надо определить по таблице R = f(H,JL) и fi = f(H,JL).

Кроме того, коэффициенты Ki, в зависимости от вида грунта принимаются:

Ki = 0,5 (пески и супеси); Ki = 0,6 - (суглинки);

Ki = 0,7 (глины с JL = 0,18); Ki = 0,8 (глины с JL = 0,25).

Для песков Ki принимается в зависимости от крупности. При расчете несущей способности свай по СНиП надо определить R и fi, которые зависят от физических показателей JL и крупности песка, но не зависят от механических показателей E,V.

1.3. Применение численных методов для расчета свай и свайных фундаментов

Теоретические методы для прогноза поведения прогноза поведения свай и свайных фундаментов развивались на основе использования решений Мелана для плоской задачи и решения Миндлина в случае пространственной задачи. Этот подход использовали в своих исследованиях Абраменко П.Г. [19], Барвашов В.А. [9], Бартоломей А.А.[10,11], Бенерджи П. и Батерфилд Р. [26] и другие.

Бартоломей А.А. [11] на основании многочисленных экспериментальных исследований предложил методику расчета осадки ленточных свайных фундаментов. Для решения задачи использована формула Горбунова-Посадова для вертикальной составляющей перемещения в случае плоской задачи при загружении основания вертикальными силами Р, приложенными на глубине Z. Формула была получена на основании фундаментального решения Е. Мелана для плоской задачи.

При решении задачи приняты следующие допущения:

1) грунт - линейно-деформируемая среда;

2) сваи и грунт в межсвайном пространстве рассматриваются как единый массив;

3) нагрузка от сваи на грунт передается через боковую поверхность сваи и массивы грунта и в плоскости нижних концов свай;

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

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

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

В настоящее время наиболее широкое распространение получили такие численные методы: метод конечных разностей (МКР), метод конечных элементов (МКЭ), метод граничных элементов (МГЭ), а также различные их модификации, включая комбинированные, то соединяющие в различном объеме выше перечисленные пути решения одной задачи, но для различных областей исследуемой среды.

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

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

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

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

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

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

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

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

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

Из выполненных различными авторами исследований [7, 12] следует, что время, которое затрачивается ЭВМ для решения трехмерных задач МГЭ и МКЭ при одинаковой точности обычно в четыре - десять раз меньше при использовании МГЭ. Эта разница может быть гораздо ощутимее для классов задач, при решении которых использование МГЭ особенно целесообразно: