Метод граничных элементов нашел применение в задачах связанных с теорией потенциала, теорией упругости, пластичности, вязкопластичности, вопросах теории теплопроводности, а также в расчетах изгибов тонких упругих пластин, колебаний деформируемых тел, распространения волн в средах, динамики жидкости.
Метод граничных элементов также может быть использован в сочетании с другими численными методами, такими как методы конечных элементов или конечных разностей, т. е. в смешанных формулировках. Соответствующие комбинированные решения почти неограниченно расширяют область применения методов, потому, что метод граничных элементов обладает четко выраженными преимуществами для областей больших размеров, в то время как методы конечных элементов являются удобным средством включения в такие системы объектов быстрого изменения свойств.
Из выполненных различными авторами исследований [7, 12] следует, что время, которое затрачивается ЭВМ для решения трехмерных задач МГЭ и МКЭ при одинаковой точности обычно в четыре - десять раз меньше при использовании МГЭ. Эта разница может быть гораздо ощутимее для классов задач, при решении которых использование МГЭ особенно целесообразно:
1. Системы, границы которых частично находятся в бесконечности. Поскольку процедуре решения задачи МГЭ автоматически удовлетворяет граничным условиям на бесконечности, отсутствует потребность в дискретизации этих границ. В то время как в методе граничных элементов границы в бесконечности должны быть аппроксимированы значительным количеством удаленных элементов.
2. Системы, содержащие полубесконечные области с ненагруженными участками свободной границы. В этом случае, нет нужды дискретизировать ненагруженные области, которые как правило, составляют большую часть свободной поверхности, если использовать подходящее фундаментальное решение, например решение Буссинеска или Миндлина.
Метод конечных разностей в области оснований и фундаментов нашел применение в расчетах конструкций на упругом основании.
В работе Клепикова С.Н. [40] освещено широкое применение МКР к большому классу задач по расчету конструкций на упругом основании, включая определение коэффициента на упругом основании и свайных оснований, расчет балок на изгиб и кручение, расчет перекрестных балок, рам, балок-стенок и плит, опирающихся на упругое основание произвольной жесткости.
Опубликовано сравнительно небольшое число работ, в которых используется МКР для расчета свай и свайных фундаментов. Можно отметить работу Федоровского В.Г. [28] в которой рассмотрена задача расчета сваи на действие продольной и поперечной нагрузки и выполнены расчеты с использованием метода Тейлора и метода конечных разностей (МКР).
Автор отмечает, что выбор метода расчета определяется не только его прогнозируемой точностью, но и быстродействием, простотой, необходимым объемом памяти в ЭВМ. В связи с этим в большинстве случаев можно рекомендовать МКР в задачах, где число расчетов велико, например при расчетах свайных кустов.
В настоящее время представляется маловероятным расширение области применения МКР в расчетах свай и свайных фундаментов в связи с тем, что разработаны и другие численные методы, например метод конечных элементов и метод граничных элементов, которые дают возможность более полно отразить реальные условия совместной работы свай и их оснований.
Более широкое применение в расчетах и проектировании свайных фундаментов получил МКЭ. Приложению МКЭ в анализе свай и свайных фундаментов посвящены исследования Бойко И.П. [29], Оттавиани М. [30], Петрашевича Г. [31] и др.
Применение метода конечных элементов в области проектирования свайных фундаментов позволило выявить важные особенности их взаимодействия с грунтами основания с учетом многих факторов, определяющих совместную работу.
Однако, необходимость дискретизации пространства занимаемого сваями и окружающим грунтом приводит к образованию значительного числа конечных элементов и как следствие большой системы линейных алгебраических уравнений. В связи с этим значительно возрастают затраты машинного времени, которое имеют высокую стоимость.
Так, по данным исследований Оттавиани М. [30] на расчет одиночной сваи было затрачено менее одной минуты машинного времени, а куста из 2´2 свай - требовалось около 200 минут, а куста из 3´3 свай -250 минут машинного времени.
Кроме того, в связи с особенностями дискретизации исследуемой трехмерной области при использовании МКЭ (большое число узлов, элементов, значительный объем данных о начальном и краевых условиях задачи) требуются увеличенные затраты времени на подготовку и ввод исходных данных в ЭВМ.
В настоящее время для анализа свайных фундаментов более конкурентоспособным является МГЭ. Так при его реализации в линейной постановке задачи требуется дискретизация только границы исследуемой области, то есть боковой поверхности сваи и плоскости подошвы ее нижнего конца.
В работах Р. Батерфилда и П.К. Бенарджи [26] рассмотрено поведение абсолютно-жестких и сжимаемых свай погруженных в однородное линейно деформируемое полупространство. Анализ поведения выполнен с использованием МГЭ. В качестве фундаментального решения использовано решение Миндлина о сосредоточенной силе приложенной внутри упругого полупространства.
Результаты исследований представлены в виде графиков, которые показывают влияние относительного заглубления сваи, отношения модуля упругости сваи к модулю сдвига грунта, влияние уширения нижнего конца сваи на перемещение одиночных свай при действии вертикальной нагрузки. Кроме того, исследовалось влияние относительно заглубления свай, расстояния между сваями и отношения модуля упругости сваи к модулю сдвига грунта полупространств на несущую способность кустов свай типичной конфигурации (2´2, 3´3 и 5 свай). Сравнение результатов своих теоретических исследований и экспериментальных исследований в полевых условиях других авторов позволило сделать вывод о том, что данные об осадке одиночной сваи можно экстраполировать на поведение группы (куста) свай.
В работе Швеца А.В. и др. [32] представлен метод определения вязкоупругого напряженно-деформированного состояния (НДС) в активной зоне биконической сваи. Рассмотрена осесимметричная задача линейной теории вязкоупругости, которая решалась методом интегральных преобразований. Отмечается, что решение упругой задачи может быть реализовано как в аналитическом виде, так и численно. В данной работе упругое решение построено методом конечных элементов. Проведены расчеты НДС в активной зоне биконической сваи. Однако, результаты расчетов и их анализ в работе не приведены.
Таким образом, анализ показывает, что численные методы, особенно МКЭ и МГЭ находят применение в исследованиях сложных явлений совместной работы свай и свайных фундаментов.
В связи с возможностью и необходимостью применения ЭВМ в расчетах оснований и фундаментов, актуальным вопросом является разработка методик использования МГЭ, который имеет значительные преимущества по сравнению с другими численными методами, особенно для областей с бесконечными границами. Учитывая, что МГЭ еще не использовался в расчетах ленточных свайных фундаментов, его применение здесь будет способствовать более полному исследованию важной в практичном отношении проблемы.
На основании анализа состояния вопроса применения коротких свай в промышленном и гражданском строительстве намечены такие задачи:
1. Разработка методики расчета бипирамидальных свай по деформациям основания с применением метода граничных элементов.
2. Анализ результатов экспериментальных данных сопротивлений бипирамидальных свай вертикальным нагрузкам.
3. Выполнение расчетов сопротивления бипирамидальных свай на ЭВМ с использованием метода граничных элементов.
4. Сравнение теоретических и экспериментальных данных сопротивления бипирамидальных свай.
Раздел 2. Применение МГЭ в расчетах сопротивления
бипирамидальных свай
2.1. Алгоритм определения сопротивления бипирамидальных свай вертикальным нагрузкам с использованием МГЭ
Алгоритм расчета свай с применением МГЭ состоит из следующих основных этапов:
- дискретизация (разбивка) поверхности фундамента в вытрамбованном котловане (боковой поверхности и нижнего конца);
- определение коэффициентов матриц влияния сил действующих на поверхности фундамента на точки (узлы) дискретизации с использованием фундаментального решения Миндлина [41];
- формирование глобальной матрицы коэффициентов влияния и свободных членов (использования граничных условий);
- решение системы линейных алгебраических уравнений т. е. боковой поверхности и в плоскости нижнего конца фундамента;
- определение сопротивления грунта на боковые поверхности и под нижним концом фундамента в вытрамбованном котловане, а так же общего сопротивления фундамента при заданной осадке.
2.2. Расчет бипирамидальных свайна ЭВМ
2.2.1. Структура программы
Расчет сопротивления бипирамидальных свай при действии вертикальной нагрузки реализован на алгоритмическом языке TurboPascal [52] с помощью программы sv63m.pas, разработанной в Винницком государственном техническом университете. Программа sv63m.pas состоит из следующих процедур:
INPUT - эта процедура считывает исходные данные: геометрические характеристики фундамента, свойства грунта, заданную осадку фундамента.
MATR - вычисляются коэффициенты влияния матрицы [K]ij и свободные коэффициенты wedi.
CAUSP - решается система линейных алгебраических уравнений, в результате определяются неизвестные значения напряжений на боковой поверхности и под нижним концом фундамента.