Реферат: Исследование сопротивления вертикальным нагрузкам бипирамидальных свай
--PAGE_BREAK--При этом существенно уменьшаются затраты труда на изготовление бипирамидальной сваи, а кроме того, по сравнению, с первым вариантом, сокращается расход арматуры. Так как в случае забивной сваи необходимо обеспечить ее целостность при транспортировании и забивке.
Как показывают экспериментальные исследования, сопротивления бипирамидальных свай имеет величину равную сопротивлению пирамидальных свай тех же размеров (длина, размер поперечного сечения в голове и нижнего конца) и при одинаковых осадках. Однако удельное сопротивление бипирамидальных свай по сравнению с пирамидальными сваями в 2,0… 2,5 раза выше. То есть, расход бетона и стали также сокращается в таких же пределах.
Однако методы расчета бипирамидальных свай до настоящего времени разработаны без использования современных численных методов, что не способствует их внедрению в практику строительства.
1.2. Методы расчета сопротивления коротких забивных свай
Для разработки надежного и эффективного проектного решения свайных фундаментов необходимо знать нагрузку, которую можно передать на сваю.
На первоначальном этапе применения свайных фундаментов, когда объем их применения был сравнительно небольшой определялась несущая способность свай и в отдельных случаях свайных фундаментов путем испытаний статической нагрузкой.
В дальнейшем были проведены исследования А.А. Луга [2], В.Н. Голубков [16], [17], посвященные определению несущей способности большого числа свай и свайных фундаментов с целью обобщения статических испытаний в различных грунтовых условиях. Эти исследования были направленными на определение несущей способности грунта по боковой поверхности сваи и под острием. В результате были составлены таблицы соответствующих расчетных сопротивлений грунта, которые вошли в СНиП [18].
С развитием техники тензометрических измерений появилось значительное число работ, в которых описаны результаты исследований распределения сил трения по боковой поверхности и доля нагрузки приходящаяся на острие. (Абраменко П.Г. [19], Бартоломей А.А. [20], Бахолдин Б.В. и Игонькин Н.Т. [21], Колесник Г.С., Шахирев В.Б., Моргун А.И. [22], Таланов Г.П., Лычев П.П. [23], Mohan D., Jain G., and Kumar V. [24], Seed H.B. and Reese L.C. [25]).
Эти исследования были направлены на уточнение характера распределения сил трения по боковой поверхности так как в СНиП [18] эпюра этих сил в однородных грунтовых условиях принята треугольной с основанием на уровне острия, а также изучению закономерности распределения усилий между боковой поверхностью и острием в процессе роста нагрузки на сваю. Эти исследования положены в основу разработки теоретических методов расчета свай, которые учитывают выявления особенностей работы свай с основанием.
В настоящее время при расчете забивных свай используется методика, изложенная в СНиПе [43], методика, разработанная в Одесском инженерно-строительном институте и изложенная во Временных указаниях по проектированию и устройству фундаментов из пирамидальных свай [25], а также методы расчета с использованием численных методов [26].
1.2.1. Определение сопротивления пирамидальных свай по методу ОИСИ
Расчет пирамидальных свай по методу ОИСИ [26] выполняется с учетом следующих основных требований:
а) средний удельный вес сухого грунта (gd, ср) уплотненного при забивке пирамидальной сваи, в пределах зоны приложения должен иметь значения 16,0 — 17,5 кн/м3;
б) величина нормативной осадки пирамидальной сваи принимается в расчете равной предельно допустимой осадке Sн = 8 см, согласно СНиП [ ], для крупнопанельных и крупноблочных бескаркасных зданий;
в) объем зоны деформаций не должен превышать объема зоны уплотнения (рис. 1.1 ).
<group id="_x0000_s1026" coordorigin=«2,-1» coordsize=«19998,20003» o:allowincell=«f»><shapetype id="_x0000_t19" coordsize=«21600,21600» o:spt=«19» adj="-5898240,,,21600,21600" path=«wr-21600,,21600,43200,,,21600,21600nfewr-21600,,21600,43200,,,21600,21600l,21600nsxe» filled=«f»><path arrowok=«t» o:extrusionok=«f» gradientshapeok=«t» o:connecttype=«custom» o:connectlocs=«0,0;21600,21600;0,21600»><img width=«291» height=«213» src=«dopb24083.zip» v:shapes="_x0000_s1026 _x0000_s1027 _x0000_s1028 _x0000_s1029 _x0000_s1030 _x0000_s1031 _x0000_s1032 _x0000_s1033 _x0000_s1034 _x0000_s1035 _x0000_s1036 _x0000_s1037 _x0000_s1038 _x0000_s1039 _x0000_s1040 _x0000_s1041 _x0000_s1042 _x0000_s1043 _x0000_s1044 _x0000_s1045"> Рв
1
2
Рис. 1.1 Схема совместной работы пирамидальной сваи и грунта основания
1 — зона уплотнения основания;
2 — зона деформаций основания.
Сопротивление пирамидальной сваи определяется по формуле:
<shapetype id="_x0000_t75" coordsize=«21600,21600» o:spt=«75» o:divferrelative=«t» path=«m@4@5l@4@11@9@11@9@5xe» filled=«f» stroked=«f»><path o:extrusionok=«f» gradientshapeok=«t» o:connecttype=«rect»><lock v:ext=«edit» aspectratio=«t»><shape id="_x0000_i1025" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image002.wmz» o:><img width=«135» height=«48» src=«dopb24084.zip» v:shapes="_x0000_i1025">
где Егр.ср — значение среднего модуля объемной деформации уплотненного грунта в пределах объема зоны деформаций, которые определяются по графику Егр.ср = f(gd), (см. [26]);
Vsc — объемная осадка сваи, которая определяется по формуле:
<shape id="_x0000_i1026" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image004.wmz» o:><img width=«151» height=«43» src=«dopb24085.zip» v:shapes="_x0000_i1026">
Su — предельно-допустимая осадка Su = 8 см;
Vc — объем погруженной части сваи;
l — длина погруженной части сваи;
b — коэффициент, принимаемый равным b = 0,5;
Vac — объем зоны деформаций пирамидальной сваи, определяемый по графику Vac = f(Vsc) [26];
Fэ — эффективная площадь поперечного сечения сваи, обуславливающая ее объемную осадку, Fэ = Vc/l.
1.2.2. Определение сопротивления пирамидальных свай по СНиП
Несущую способность Fd, пирамидальной сваи с наклоном боковых граней ip> 0,25 допускается определять как сумму сил расчетных сопротивлений грунта основания на боковой поверхности сваи и под ее нижним концом по формуле[18, 27]:
<shape id="_x0000_i1027" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image006.wmz» o:><img width=«351» height=«48» src=«dopb24086.zip» v:shapes="_x0000_i1027">
где Ai — площадь боковой поверхности сваи в пределах i-го слоя грунта, м2;
a — угол наклона граней пирамидальной сваи, град.;
jli, cli — расчетные значения угла внутреннего трения, град., и сечения, кПа;
d — размер стороны нижнего конца сваи, м;
n1, n2 — коэффициенты, значения которых определяются по таблице СНиПа [ 18 ] приложения 2.
Сопротивление грунта под острием сваи P'i и по ее боковой поверхности Рi определяется по формуле:
<shape id="_x0000_i1028" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image008.wmz» o:><img width=«533» height=«80» src=«dopb24087.zip» v:shapes="_x0000_i1028">
где Ei — модуль деформации i-го слоя грунта, определяемый по результатам прессиометрических испытаний, кПа;
Vi — коэффициент Пуассона i-го слоя грунта, принимаемый в соответствии с требованиями главы СНиПа «Основания зданий и сооружений» [ 43 ];
x — коэффициент, значения которого определяются по таблице СНиПа [18].
Природное боковое давление грунта Рoi (кПа) определяют по формуле:
<shape id="_x0000_i1029" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image010.wmz» o:><img width=«151» height=«52» src=«dopb24088.zip» v:shapes="_x0000_i1029">
где gi — удельный вес грунта i-го слоя кН/м3;
hi — средняя глубина расположения i-го слоя грунта, м.
Начальное давление грунта:
Ppi = Poi(1 + sinji) + cicosji.
При отсутствии прессиометрических испытаний грунта несущая способность пирамидальной сваи определяется по формуле:
Fd = gc[RA + Shi(uifi + u0ipEiKixr)],
где gc — коэффициент условий работы сваи в грунте, принимаемый gc = 1;
R — расчетное сопротивление грунта под нижним концом сваи, кПа, принимаемое по таблице №1 [18];
A — площадь сечения конца сваи;
hi — толщина i-го слоя грунта, соприкасающегося с боковой поверхностью сваи, м;
fi — расчетное сопротивление i-го слоя грунта основания на боковой поверхности сваи, кПа, принимаемое по таблице 2 [18];
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] следует, что время, которое затрачивается ЭВМ для решения трехмерных задач МГЭ и МКЭ при одинаковой точности обычно в четыре — десять раз меньше при использовании МГЭ. Эта разница может быть гораздо ощутимее для классов задач, при решении которых использование МГЭ особенно целесообразно:
1. Системы, границы которых частично находятся в бесконечности. Поскольку процедуре решения задачи МГЭ автоматически удовлетворяет граничным условиям на бесконечности, отсутствует потребность в дискретизации этих границ. В то время как в методе граничных элементов границы в бесконечности должны быть аппроксимированы значительным количеством удаленных элементов.
2. Системы, содержащие полубесконечные области с ненагруженными участками свободной границы. В этом случае, нет нужды дискретизировать ненагруженные области, которые как правило, составляют большую часть свободной поверхности, если использовать подходящее фундаментальное решение, например решение Буссинеска или Миндлина.
Метод конечных разностей в области оснований и фундаментов нашел применение в расчетах конструкций на упругом основании.
В работе Клепикова С.Н. [40] освещено широкое применение МКР к большому классу задач по расчету конструкций на упругом основании, включая определение коэффициента на упругом основании и свайных оснований, расчет балок на изгиб и кручение, расчет перекрестных балок, рам, балок-стенок и плит, опирающихся на упругое основание произвольной жесткости.
Опубликовано сравнительно небольшое число работ, в которых используется МКР для расчета свай и свайных фундаментов. Можно отметить работу Федоровского В.Г. [28] в которой рассмотрена задача расчета сваи на действие продольной и поперечной нагрузки и выполнены расчеты с использованием метода Тейлора и метода конечных разностей (МКР).
Автор отмечает, что выбор метода расчета определяется не только его прогнозируемой точностью, но и быстродействием, простотой, необходимым объемом памяти в ЭВМ. В связи с этим в большинстве случаев можно рекомендовать МКР в задачах, где число расчетов велико, например при расчетах свайных кустов.
В настоящее время представляется маловероятным расширение области применения МКР в расчетах свай и свайных фундаментов в связи с тем, что разработаны и другие численные методы, например метод конечных элементов и метод граничных элементов, которые дают возможность более полно отразить реальные условия совместной работы свай и их оснований.
Более широкое применение в расчетах и проектировании свайных фундаментов получил МКЭ. Приложению МКЭ в анализе свай и свайных фундаментов посвящены исследования Бойко И.П. [29], Оттавиани М. [30], Петрашевича Г. [31] и др.
продолжение
--PAGE_BREAK--Применение метода конечных элементов в области проектирования свайных фундаментов позволило выявить важные особенности их взаимодействия с грунтами основания с учетом многих факторов, определяющих совместную работу.
Однако, необходимость дискретизации пространства занимаемого сваями и окружающим грунтом приводит к образованию значительного числа конечных элементов и как следствие большой системы линейных алгебраических уравнений. В связи с этим значительно возрастают затраты машинного времени, которое имеют высокую стоимость.
Так, по данным исследований Оттавиани М. [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. Структура программы
Расчет сопротивления бипирамидальных свай при действии вертикальной нагрузки реализован на алгоритмическом языке Turbo Pascal [52] с помощью программы sv63m.pas, разработанной в Винницком государственном техническом университете. Программа sv63m.pas состоит из следующих процедур:
INPUT — эта процедура считывает исходные данные: геометрические характеристики фундамента, свойства грунта, заданную осадку фундамента.
MATR — вычисляются коэффициенты влияния матрицы [K]ij и свободные коэффициенты wedi.
CAUSP — решается система линейных алгебраических уравнений, в результате определяются неизвестные значения напряжений на боковой поверхности и под нижним концом фундамента.
OUTPUT — определяются касательные напряжения по боковой поверхности фундамента и нормальные напряжения под нижним концом, а так же радиальные напряжения действующие на боковую поверхность фундамента; определяются сосредоточенные силы действующие на i-х элементах боковой поверхности (силы трения) и нижнего конца фундамента — нормальные силы, сумма соответствующих сил дает значения общего усилия по боковой поверхности и под нижним концом, а их сумма общее сопротивление фундамента.
В программе используются следующие основные переменные:
NE1 := NEA + NEB + NEC — число граничных элементов на боковой поверхности фундамента;
NN1 — число граничных узлов на боковой поверхности фундамента;
NE2 — число граничных элементов в плоскости нижнего конца фундамента;
NN2 — число граничных узлов в плоскости нижнего конца фундамента;
NE3 — число граничных элементов по окружности фундамента;
NN3 — число граничных элементов по окружности фундамента;
ls1 — длина первого (верхнего) участка фундамента;
ls2 — длина второго (среднего) участка фундамента;
ls3 — длина третьего (нижнего) участка фундамента;
ls := ls1 + ls2 +ls3 — общая длина фундамента;
E — модуль деформации грунта;
mu — коэффициент Пуассона для грунта;
ed1 — вертикальные перемещения узлов боковой поверхности фундамента;
ed2 — горизонтальные перемещения узлов боковой поверхности фундамента;
ed3 — вертикальные перемещения узлов нижнего конца фундамента;
ar1 — радиус фундамента в верхнем сечении I первого участка;
ars — радиус фундамента в нижнем сечении среднего участка;
arN — величина радиуса фундамента на уровне нижнего конца фундамента;
NE = NE1 + NE2 — число граничных элементов на поверхности фундамента;
NK1 := NE1 + 1 — номер элемента матрицы К из
NEE = 2 * NE1 — номер элемента глобальной матрицы К
NC2 := NЕЕ +1 — номер элемента глобальной матрицы К.
tga1 — тангенс угла наклона боковой поверхности (грани) среднего участка фундамента;
tga2 — тангенс угла наклона боковой поверхности нижнего участка фундамента;
NEA — число граничных элементов на первом (верхнем) участке фундамента в вытрамбованном котловане;
NEB — число граничных элементов на втором участке фундамента;
NEC — число граничных элементов на третьем (нижнем) участке фундамента;
HH1 — шаг граничных узлов на первом участке;
HH2 — шаг граничных узлов на втором участке;
HH3 — шаг граничных узлов на третьем участке;
inz [i,1], inz [i,2] — связность граничных элементов боковой поверхности фундамента;
inc [i,1], inc [i,2] — связность элементов нижнего конца фундамента;
int [i,1], int [i,2] — связность элементов окружности по боковой поверхности фундамента и в плоскости нижнего конца фундамента (в точках источников);
2.2.2. Дискретизация боковой поверхности и нижнего конца фундамента
<line id="_x0000_s1046" from=«304.8pt,-247.5pt» to=«319.25pt,-247.45pt» o:allowincell=«f» strokeweight=«1pt»><img width=«22» height=«2» src=«dopb24089.zip» v:shapes="_x0000_s1046"> 1
1
2 I
2
3
3
4
4 II
5
5
6
<group id="_x0000_s1047" coordorigin=«1,1» coordsize=«19998,19984» o:allowincell=«f»><img width=«233» height=«549» src=«dopb24090.zip» v:shapes="_x0000_s1047 _x0000_s1048 _x0000_s1049 _x0000_s1050 _x0000_s1051 _x0000_s1052 _x0000_s1053 _x0000_s1054 _x0000_s1055 _x0000_s1056 _x0000_s1057 _x0000_s1058 _x0000_s1059 _x0000_s1060 _x0000_s1061 _x0000_s1062 _x0000_s1063 _x0000_s1064 _x0000_s1065 _x0000_s1066 _x0000_s1067 _x0000_s1068 _x0000_s1069 _x0000_s1070 _x0000_s1071 _x0000_s1072 _x0000_s1073 _x0000_s1074 _x0000_s1075 _x0000_s1076 _x0000_s1077 _x0000_s1078 _x0000_s1079 _x0000_s1080 _x0000_s1081 _x0000_s1082 _x0000_s1083 _x0000_s1084 _x0000_s1085 _x0000_s1086 _x0000_s1087 _x0000_s1088 _x0000_s1089 _x0000_s1090 _x0000_s1091 _x0000_s1092 _x0000_s1093 _x0000_s1094 _x0000_s1095 _x0000_s1096 _x0000_s1097 _x0000_s1098 _x0000_s1099 _x0000_s1100 _x0000_s1101 _x0000_s1102 _x0000_s1103 _x0000_s1104"> 6
7
7
8
8
9
9 III
10
11
12
13
Рис. 2.1. Схема дискретизации боковой поверхности
фундамента в вытрамбованном котловане
t, t
1 2 3 4 5 6 (NN2)
<group id="_x0000_s1105" coordsize=«19995,20000» o:allowincell=«f»><img width=«253» height=«60» src=«dopb24091.zip» v:shapes="_x0000_s1105 _x0000_s1106 _x0000_s1107 _x0000_s1108 _x0000_s1109 _x0000_s1110 _x0000_s1111 _x0000_s1112 _x0000_s1113 _x0000_s1114 _x0000_s1115 _x0000_s1116 _x0000_s1117 _x0000_s1118 _x0000_s1119 _x0000_s1120 _x0000_s1121 _x0000_s1122"> 0 ar
1 2 3 4 5 (NE2)
Рис. 2.2. Схема дискретизации нижнего конца фундамента
По длине фундамента в вытрамбованном котловане разбивается на три участка: верхний, средний (II), нижний (III) (рис. 2.1).
Количество граничных элементов задается в пределах каждого участка соответственно: NEA, NEB, NEC. Кроме того, для каждого участка задается длина (ls1, ls2, ls3). Угол наклона боковой поверхности участков II и III задан тангенсом угла наклона (tga1 и tga2) (см. рис. 2.3).
<shape id="_x0000_i1030" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image015.wmz» o:><img width=«29» height=«56» src=«dopb24092.zip» v:shapes="_x0000_i1030">
<group id="_x0000_s1123" coordsize=«20000,20165» o:allowincell=«f»><img width=«156» height=«271» src=«dopb24093.zip» v:shapes="_x0000_s1123 _x0000_s1124 _x0000_s1125 _x0000_s1126 _x0000_s1127 _x0000_s1128 _x0000_s1129 _x0000_s1130 _x0000_s1131 _x0000_s1132 _x0000_s1133 _x0000_s1134 _x0000_s1135"> a1 <shape id="_x0000_i1031" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image018.wmz» o:><img width=«37» height=«85» src=«dopb24094.zip» v:shapes="_x0000_i1031">
a2 <shape id="_x0000_i1032" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image020.wmz» o:><img width=«45» height=«141» src=«dopb24095.zip» v:shapes="_x0000_i1032">
Рис. 2.3.
При известных длине участков и количестве граничных элементов на них определяются коэффициенты i-узлов по длине фундамента:
Z[i] = Z[i-1] + HH1 — I участок;
Z[i] = Z[i-1] + HH2 — II участок;
Z[i] = Z[i-1] + HH3 — II участок,
где <shape id="_x0000_i1033" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image022.wmz» o:><img width=«165» height=«85» src=«dopb24096.zip» v:shapes="_x0000_i1033"> - шаг граничных узлов на боковой поверхности фундамента в вытрамбованном котловане.
Узлы qi при обходе граничных элементов по окружности при заданном числе элементов NE3 и диапазона изменения угла q = 0...p определяем по формуле (см. рис. 2.4):
Ai = Ai-1 + H3,
где H3 = p/NE3 - шаг граничных узлов по окружности радиус которой, равен радиусу узла в точке приложения (j).
<group id="_x0000_s1136" coordorigin="-1,1" coordsize=«20002,19999» o:allowincell=«f»><img width=«159» height=«169» src=«dopb24097.zip» v:shapes="_x0000_s1136 _x0000_s1137 _x0000_s1138 _x0000_s1139 _x0000_s1140 _x0000_s1141 _x0000_s1142 _x0000_s1143 _x0000_s1144 _x0000_s1145 _x0000_s1146 _x0000_s1147 _x0000_s1148 _x0000_s1149 _x0000_s1150 _x0000_s1151 _x0000_s1152 _x0000_s1153"> p/2
q
p 0
Рис. 2.4.
Радиус i-го узла на боковой поверхности фундамента в вытрамбованном котловане определим при известных его значениях ar1, ars, arN и тангенсах угла наклона tga1, tga2 по формуле
I участок
ar[i]=ar1;
II участок
ar[i]=ar[i-1] — tga1 * HH2;
III участок
ar[i]=ar[i-1] — tga1 * HH3.
Координаты узлов в плоскости нижнего конца фундамента определим из следующих соотношений (см. рис. 2.5)
координат по длине фундамента Z[i]=ls;
(ls — общая длина фундамента в вытрамбованном котловане),
координат в радиальном направлении ar[i]=ar[i+1] + H2,
где H2 — шаг узлов, находящихся на нижнем конце фундамента.
ar[NE1 + 1]
ar[NE1 + 2]<group id="_x0000_s1154" coordorigin=",-1800" coordsize=«20000,22464» o:allowincell=«f»><img width=«252» height=«252» src=«dopb24098.zip» v:shapes="_x0000_s1154 _x0000_s1155 _x0000_s1156 _x0000_s1157 _x0000_s1158 _x0000_s1159 _x0000_s1160 _x0000_s1161 _x0000_s1162 _x0000_s1163 _x0000_s1164 _x0000_s1165 _x0000_s1166 _x0000_s1167 _x0000_s1168 _x0000_s1169 _x0000_s1170 _x0000_s1171 _x0000_s1172 _x0000_s1173 _x0000_s1174">
ar[NE + 1]=0
Рис. 2.5. Схема узлов на нижнем конце фундамента
В работе использовано понятие «связность элементов». Так как производится дискретизация поверхности фундамента в условиях осессимметричной задачи, то граничные элементы представлены прямыми линиями находящимися между граничными узлами и каждый граничный элемент, определяется если задать узлы которые его ограничивают (рис. 2.6).
<group id="_x0000_s1175" coordorigin=«1» coordsize=«19999,20000» o:allowincell=«f»><img width=«98» height=«92» src=«dopb24099.zip» v:shapes="_x0000_s1175 _x0000_s1176 _x0000_s1177 _x0000_s1178 _x0000_s1179"> 2
i
1
Рис. 2.6. Схема к понятию связности элементов
В данной работе для наглядности введены отдельно связности i-х элементов на боковой поверхности фундамента, в плоскости нижнего конца, и по окружности фундамента:
inz[i,1] inz[i,2],
inc[i,1] inc[i,2],
int[i,1] int[i,2],
<group id="_x0000_s1180" coordsize=«20000,20000» o:allowincell=«f»><img width=«45» height=«21» src=«dopb24100.zip» v:shapes="_x0000_s1180 _x0000_s1181 _x0000_s1182">где i — номер граничного элемента;
1, 2 — номера граничных узлов, окружающих связывающий i-й элемент (см. рис. 2.6).
2.2.3. Формирование матрицы коэффициентов влияния и свободных членов СЛАУ
При формировании коэффициентов глобальной матрицы влияния, отражающих зависимость перемещения точки наблюдения (i), когда источник возмущения находится в точке (j) используется решение Миндлина для силы приложений внутри упругого полупространства. Иногда для зависимости, когда действует единичная сила, эти решения называют фундаментальными. Для вертикальной силы Рв=1 зависимость для перемещений KW, когда точка наблюдения имеет координаты В(z,r), а источник возмущения находится на оси Z (радиальная координата равна нулю) на глубине с, запишется в виде:
<group id="_x0000_s1183" coordsize=«20001,20001» o:allowincell=«f»><img width=«291» height=«224» src=«dopb24101.zip» v:shapes="_x0000_s1183 _x0000_s1184 _x0000_s1185 _x0000_s1186 _x0000_s1187 _x0000_s1188 _x0000_s1189 _x0000_s1190 _x0000_s1191 _x0000_s1192 _x0000_s1193 _x0000_s1194 _x0000_s1195 _x0000_s1196 _x0000_s1197 _x0000_s1198">
с 0 0
r
с N
<shape id="_x0000_i1034" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image029.wmz» o:><img width=«21» height=«25» src=«dopb24102.zip» v:shapes="_x0000_i1034">\s <shape id="_x0000_i1035" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image031.wmz» o:><img width=«19» height=«20» src=«dopb24103.zip» v:shapes="_x0000_i1035">\s
Рв
x(с,0) r B(z,r)
Z
Рис. 2.7. Схема обозначений в формуле Миндлина для сосредоточенной силы Рв, приложенной внутри упругого полупространства
<shape id="_x0000_i1036" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image033.wmz» o:><img width=«415» height=«104» src=«dopb24104.zip» v:shapes="_x0000_i1036"> (2.1)
где
<shape id="_x0000_i1037" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image035.wmz» o:><img width=«141» height=«64» src=«dopb24105.zip» v:shapes="_x0000_i1037"> (2.2)
<shape id="_x0000_i1038" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image037.wmz» o:><img width=«120» height=«56» src=«dopb24106.zip» v:shapes="_x0000_i1038"> (2.3)
G — модуль сдвига грунта;
E — модуль деформации грунта;
v — коэффициент Пуассона грунта.
KW — вертикальное перемещение точки В при действии вертикальной силы Рв=1 в точке x(0, с).
Применение решения Миндлина к задаче о сопротивлении фундамента вертикальной нагрузке состоит в том, что точка приложения силы и точка наблюдения, в которой возникают вертикальные перемещения находятся на боковой поверхности или на нижнем конце. В связи с этим в формуле (2.1) выражения для R1 и R2 принимают вид:
продолжение
--PAGE_BREAK-- <shape id="_x0000_i1039" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image039.wmz» o:><img width=«140» height=«28» src=«dopb24107.zip» v:shapes="_x0000_i1039"> (2.4)
<shape id="_x0000_i1040" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image041.wmz» o:><img width=«147» height=«28» src=«dopb24108.zip» v:shapes="_x0000_i1040"> (2.5)
где
<shape id="_x0000_i1041" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image043.wmz» o:><img width=«289» height=«44» src=«dopb24109.zip» v:shapes="_x0000_i1041"> (2.6)
r — горизонтальная компонента расстояния от оси Z до точки B;
arc — горизонтальная компонента расстояния от оси Z до точки x;
r1 — горизонтальная компонента расстояния от точки В (точки наблюдения) до точки x (источник, место приложения силы);
R2 — расстояние от точки x' (фиктивный источник) до точки B;
<group id="_x0000_s1199" coordsize=«20000,20000» o:allowincell=«f»><img width=«263» height=«203» src=«dopb24110.zip» v:shapes="_x0000_s1199 _x0000_s1200 _x0000_s1201 _x0000_s1202 _x0000_s1203 _x0000_s1204 _x0000_s1205 _x0000_s1206 _x0000_s1207 _x0000_s1208 _x0000_s1209 _x0000_s1210 _x0000_s1211 _x0000_s1212 _x0000_s1213">R1 — расстояние от точки x (источник) до точки B.
x(с,arc)
<shape id="_x0000_i1042" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image046.wmz» o:><img width=«22» height=«23» src=«dopb24111.zip» v:shapes="_x0000_i1042">\s
<shape id="_x0000_i1043" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image048.wmz» o:><img width=«16» height=«14» src=«dopb24112.zip» v:shapes="_x0000_i1043">\s
q <shape id="_x0000_i1044" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image050.wmz» o:><img width=«12» height=«12» src=«dopb24113.zip» v:shapes="_x0000_i1044">\s B(z,r)
a
Рис. 2.8. Схема к определению координат точки приложения x(с,arc) и точки наблюдения B(z,r)
При определении коэффициентов влияния глобальной матрицы К учитываются различные варианты расположения источников (сил) и точек наблюдения.
<group id="_x0000_s1214" coordsize=«20000,20167» o:allowincell=«f»><img width=«177» height=«195» src=«dopb24114.zip» v:shapes="_x0000_s1214 _x0000_s1215 _x0000_s1216 _x0000_s1217 _x0000_s1218 _x0000_s1219 _x0000_s1220 _x0000_s1221 _x0000_s1222 _x0000_s1223 _x0000_s1224 _x0000_s1225 _x0000_s1226 _x0000_s1227 _x0000_s1228 _x0000_s1229 _x0000_s1230 _x0000_s1231 _x0000_s1232 _x0000_s1233 _x0000_s1234 _x0000_s1235 _x0000_s1236 _x0000_s1237 _x0000_s1238 _x0000_s1239 _x0000_s1240 _x0000_s1241 _x0000_s1242">
dc <shape id="_x0000_i1045" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image053.wmz» o:><img width=«28» height=«32» src=«dopb24115.zip» v:shapes="_x0000_i1045">
· i
Рис. 2.9. Схема к интегрированию решения Миндлина
(матрица KSS)
— источники расположены на боковой поверхности фундамента и точки наблюдения так же находятся на боковой поверхности. Для наглядности рассмотрим фундамент в вытрамбованном котловане (см. рис. 2.1) боковая поверхность которого разбита на j элементов (j=1,NE1) и имеются точки наблюдения i, находящиеся посредине граничных элементов. При вычислении коэффициента влияния входящего в матрицу [KSS]ij осуществляется интегрирование решения Миндлина по окружности находящейся на глубине с и радиусом arc и интегрирования полученных значений решения по высоте j-го элемента. Таким образом элементы подматрицы [KSS]ij определяются
<shape id="_x0000_i1046" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image055.wmz» o:><img width=«299» height=«59» src=«dopb24116.zip» v:shapes="_x0000_i1046"> (2.7)
где <shape id="_x0000_i1047" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image057.wmz» o:><img width=«308» height=«44» src=«dopb24117.zip» v:shapes="_x0000_i1047"> (2.8)
· i
<group id="_x0000_s1243" coordsize=«20000,20000» o:allowincell=«f»><img width=«147» height=«262» src=«dopb24118.zip» v:shapes="_x0000_s1243 _x0000_s1244 _x0000_s1245 _x0000_s1246 _x0000_s1247 _x0000_s1248 _x0000_s1249 _x0000_s1250 _x0000_s1251 _x0000_s1252 _x0000_s1253 _x0000_s1254 _x0000_s1255 _x0000_s1256 _x0000_s1257 _x0000_s1258 _x0000_s1259 _x0000_s1260 _x0000_s1261 _x0000_s1262 _x0000_s1263 _x0000_s1264 _x0000_s1265"> j
·
Рис. 2.10. Схема к интегрированию решения Миндлина
(матрица KBS)
— источники находятся на нижнем конце фундамента, а точки наблюдения на боковой поверхности. Количество элементов на нижнем конце j (1,NE2), а количество точек на боковой поверхности i=1,NE1. Интегрирование решения Миндлина выполняется по граничных элементам нижнего конца, представленных в виде кольца (рис. 2.10). При этом формируются коэффициенты подматрицы [KBS]ij
<shape id="_x0000_i1048" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image060.wmz» o:><img width=«307» height=«57» src=«dopb24119.zip» v:shapes="_x0000_i1048"> (2.9)
где <shape id="_x0000_i1049" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image062.wmz» o:><img width=«315» height=«41» src=«dopb24120.zip» v:shapes="_x0000_i1049"> (2.10)
r — горизонтальная компонента расстояния от оси Z до точки В;
eps — горизонтальное расстояние от оси Z до точки источника x;
de — ширина граничного элемента j нижнего конца фундамента (ширина кольца).
<group id="_x0000_s1266" coordsize=«20000,20000» o:allowincell=«f»><img width=«146» height=«269» src=«dopb24121.zip» v:shapes="_x0000_s1266 _x0000_s1267 _x0000_s1268 _x0000_s1269 _x0000_s1270 _x0000_s1271 _x0000_s1272 _x0000_s1273 _x0000_s1274 _x0000_s1275 _x0000_s1276 _x0000_s1277 _x0000_s1278 _x0000_s1279 _x0000_s1280 _x0000_s1281 _x0000_s1282 _x0000_s1283 _x0000_s1284 _x0000_s1285 _x0000_s1286 _x0000_s1287 _x0000_s1288 _x0000_s1289 _x0000_s1290">
<shape id="_x0000_i1050" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image053.wmz» o:><img width=«28» height=«29» src=«dopb24122.zip» v:shapes="_x0000_i1050">
i
· ·
Рис. 2.11. Схема к интегрированию решения Миндлина
(матрица KSB)
Если источники находятся на боковой поверхности фундамента, а точки наблюдения на нижнем конце. здесь формируются коэффициенты подматрицы [KSB]ij, i=1,NE2 j=1,NE1, которые учитывают влияние загружения боковой поверхности фундамента на перемещение элементов нижнего конца
<shape id="_x0000_i1051" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image066.wmz» o:><img width=«304» height=«57» src=«dopb24123.zip» v:shapes="_x0000_i1051"> (2.11)
где <shape id="_x0000_i1052" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image068.wmz» o:><img width=«304» height=«44» src=«dopb24124.zip» v:shapes="_x0000_i1052"> (2.12)
j (элемент j)
<group id="_x0000_s1291" coordorigin=",3" coordsize=«20000,19996» o:allowincell=«f»><img width=«166» height=«290» src=«dopb24125.zip» v:shapes="_x0000_s1291 _x0000_s1292 _x0000_s1293 _x0000_s1294 _x0000_s1295 _x0000_s1296 _x0000_s1297 _x0000_s1298 _x0000_s1299 _x0000_s1300 _x0000_s1301 _x0000_s1302 _x0000_s1303 _x0000_s1304 _x0000_s1305 _x0000_s1306 _x0000_s1307 _x0000_s1308 _x0000_s1309 _x0000_s1310 _x0000_s1311 _x0000_s1312 _x0000_s1313 _x0000_s1314 _x0000_s1315 _x0000_s1316 _x0000_s1317 _x0000_s1318">
i (точка наблюдения i)
· ·
Рис. 2.12. Схема к интегрированию решения Миндлина
матрицы (КВВ)
Последний вариант взаимодействия частей фундамента, когда источники находятся на нижнем конце фундамента, а точка наблюдения так же находится на нижнем конце фундамента.
Для вычисления коэффициентов влияния загружения элементов нижнего конца (j=1,NE2) на точки наблюдения, находящиеся посередине элементов нижнего конца, вычисляется двойной интервал
<shape id="_x0000_i1053" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image071.wmz» o:><img width=«309» height=«57» src=«dopb24126.zip» v:shapes="_x0000_i1053"> (2.13)
где <shape id="_x0000_i1054" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image062.wmz» o:><img width=«315» height=«41» src=«dopb24120.zip» v:shapes="_x0000_i1054">
Если учитываются вертикальные перемещения грунта примыкающего к поверхности фундамента, только от действия вертикальных сил, приложенных на боковой поверхности (KSS, KSB) и на нижнем конце (KBS, KBB), то глобальная матрица К имеет вид
<shape id="_x0000_i1055" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image073.wmz» o:><img width=«168» height=«84» src=«dopb24127.zip» v:shapes="_x0000_i1055"> (2.14)
Система алгебраических уравнений для определения неизвестных напряжений на боковой поверхности и под нижним концом записывается следующим образом
<shape id="_x0000_i1056" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image075.wmz» o:><img width=«229» height=«84» src=«dopb24128.zip» v:shapes="_x0000_i1056"> (2.15)
где fsb — неизвестные напряжения на поверхности фундамента;
wed — вектор-столбец единичных перемещений узлов поверхности фундамента. В случае, если принять сваю абсолютно жесткой (т. е. несжимаемой), то перемещения всех узлов будут одинаковыми. В данной работе компоненты вектора-столбца wed принимались равными осадке фундамента при которой график зависимости «нагрузки-осадки» имеет прямолинейный вид. Как показывает анализ опытных данных для призматических свай такая осадка равна 0,01 м, для пирамидальных и фундаментов в вытрамбованном котловане — 0,015..0,020 м.
Если учитывать, что на боковую поверхность фундамента действуют радиальные напряжения s2, то глобальная матрица [K] будет содержать девять подматриц и уравнение равновесия (2.15) примет вид:
<shape id="_x0000_i1057" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image077.wmz» o:><img width=«295» height=«143» src=«dopb24129.zip» v:shapes="_x0000_i1057"> (2.16)
где KRS — матрица, которая содержит коэффициенты влияния на вертикальные перемещения узлов боковой поверхности фундамента, при загружении элементов боковой поверхности радиальными напряжениями s2 (sigm2);
KSU — матрица, коэффициенты которой отражают связь между горизонтальными перемещениями узлов боковой поверхности фундамента, когда боковая поверхность загружена вертикальными напряжениями;
KRU — матрица содержащая коэффициенты влияния, которые отражают зависимость между горизонтальными перемещениями узлов боковой поверхности фундамента при загружении элементов боковой поверхности горизонтального напряжения s2;
KBU — матрица, коэффициенты которой отражают зависимость горизонтальных перемещений узлов боковой поверхности фундамента при загружении элементов нижнего конца вертикальными напряжениями s1;
KRB — матрица, коэффициенты которой отражают связь между вертикальными перемещениями узлов нижнего конца фундамента при загружении элементов боковой поверхности радиальными напряжениями s2.
{fsb} — вектор-столбец, содержащий неизвестные: касательные напряжения на боковой поверхности фундамента t, горизонтальные напряжения на боковой поверхности фундамента s2 и вертикальные напряжения на нижнем конце фундамента s1;
<shape id="_x0000_i1058" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image079.wmz» o:><img width=«125» height=«80» src=«dopb24130.zip» v:shapes="_x0000_i1058"> - вектор-столбец, содержащий заданные вертикальные перемещения узлов боковой поверхности фундамента ed1; горизонтальные перемещения узлов боковой поверхности ed2 (если свая не сжимается ed2=0); вертикальные перемещения узлов нижнего конца фундамента ed3.
Фундаментальное решение Миндлина в матрицах KRS и KRB имеет следующее выражение:
<shape id="_x0000_i1059" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image081.wmz» o:><img width=«348» height=«107» src=«dopb24131.zip» v:shapes="_x0000_i1059"> (2.17)
где
<shape id="_x0000_i1060" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image083.wmz» o:><img width=«175» height=«28» src=«dopb24132.zip» v:shapes="_x0000_i1060"> (2.19)
<shape id="_x0000_i1061" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image085.wmz» o:><img width=«177» height=«28» src=«dopb24133.zip» v:shapes="_x0000_i1061"> (2.20)
x = r×cosq — arc; (2.21)
y = -r×sinq. (2.22)
Коэффициенты матрицы KRS вычисляются с использованием фундаментального решения Миндлина KW3 и интегрирования выражения
<shape id="_x0000_i1062" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image087.wmz» o:><img width=«307» height=«57» src=«dopb24134.zip» v:shapes="_x0000_i1062"> (2.23)
где r = arz. (2.24)
Коэффициенты матрицы KRB вычисляются с использованием фундаментального решения Миндлина KW3 и интегрирования выражения
<shape id="_x0000_i1063" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image089.wmz» o:><img width=«315» height=«57» src=«dopb24135.zip» v:shapes="_x0000_i1063"> (2.25)
где <shape id="_x0000_i1064" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image068.wmz» o:><img width=«304» height=«44» src=«dopb24124.zip» v:shapes="_x0000_i1064"> (2.26)
При вычислении коэффициентов матриц KSU и KBU используется решение Миндлина
<shape id="_x0000_i1065" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image091.wmz» o:><img width=«403» height=«131» src=«dopb24136.zip» v:shapes="_x0000_i1065"> (2.27)
где R1, R2, r1 — определяются по формулам (2.4), (2.5), (2.6).
Коэффициенты матрицы KSU вычисляются интегрированием выражения
<shape id="_x0000_i1066" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image093.wmz» o:><img width=«312» height=«57» src=«dopb24137.zip» v:shapes="_x0000_i1066"> (2.28)
где <shape id="_x0000_i1067" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image095.wmz» o:><img width=«341» height=«41» src=«dopb24138.zip» v:shapes="_x0000_i1067"> (2.29)
Коэффициенты матрицы KBU равны интегралу
<shape id="_x0000_i1068" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image097.wmz» o:><img width=«317» height=«57» src=«dopb24139.zip» v:shapes="_x0000_i1068"> (2.30)
где <shape id="_x0000_i1069" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image099.wmz» o:><img width=«348» height=«41» src=«dopb24140.zip» v:shapes="_x0000_i1069"> (2.31)
Фундаментальное решение Миндлина в матрице KRU определяется формулой
<shape id="_x0000_i1070" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image101.wmz» o:><img width=«564» height=«239» src=«dopb24141.zip» v:shapes="_x0000_i1070">(2.32)
где R1, R2, x, y — определяются по формулам (2.19), (2.20), (2.21), (2.22).
Коэффициенты матрицы KRU определяются интегралом
<shape id="_x0000_i1071" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image103.wmz» o:><img width=«309» height=«59» src=«dopb24142.zip» v:shapes="_x0000_i1071"> (2.33)
где r = arz. (2.34)
2.2.4. Определение напряжений на поверхности фундамента
Когда сформирована глобальная матрица К и задан вектор-столбец
<shape id="_x0000_i1072" type="#_x0000_t75" o:ole=""><imagedata src=«1.files/image105.wmz» o:><img width=«131» height=«85» src=«dopb24143.zip» v:shapes="_x0000_i1072"> (2.35) решается система алгебраических уравнений (2.16) методом Гаусса с помощью процедуры GAUSP, в результате получим значения напряжений t и s2 в узлах боковой поверхности и напряжение s1 в узлах нижнего конца фундамента.
2.2.5. Определение общего сопротивления фундамента
Усилия на элементах боковой поверхности фундамента получим
<shape id="_x0000_i1073" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image107.wmz» o:><img width=«100» height=«49» src=«dopb24144.zip» v:shapes="_x0000_i1073"> (2.36)
а усилия на элементах нижнего конца
<shape id="_x0000_i1074" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image109.wmz» o:><img width=«107» height=«49» src=«dopb24145.zip» v:shapes="_x0000_i1074"> (2.37)
Суммарное значение силы трения определяется
<shape id="_x0000_i1075" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image111.wmz» o:><img width=«81» height=«45» src=«dopb24146.zip» v:shapes="_x0000_i1075"> (2.38)
а сила под нижним концом
<shape id="_x0000_i1076" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image113.wmz» o:><img width=«81» height=«45» src=«dopb24147.zip» v:shapes="_x0000_i1076"> (2.39)
Общее сопротивление фундамента при заданной осадке r = ed1 равно
Рс = Рб + Р0; (2.40)
Таким образом в результате применения изложенной методики расчета по методу граничных элементов с использованием решения Миндлина можно определить общее сопротивление фундамента в вытрамбованном котловане при заданной осадке.
Раздел 3. Результаты теоретических исследований сопротивления бипирамидальных свай
В данной работе согласно, описанной в разделе 2 методике, выполнены расчеты сопротивления бипирамидальных свай для грунтовых условий и типоразмеров свай по результатам исследований, представленных в работах [10 ¸ 15]. Теоретические модели взаимодействия свай в этих работах построены на основе теории проф. Голубкова В.Н. с использованием понятий зон уплотнения и деформаций. Эта теория построена на применении опытных данных, имеет полуэмпирический характер и требует дальнейшего развития.
Сравнение результатов экспериментов, выполненных в натурных условиях и расчетов с использованием метода граничных элементов позволяет оценить достоверность и надежность нового метода прогноза осадок бипирамидальных свай.
Далее кратко рассмотрены результаты полевых исследований сопротивления бипирамидальных свай выполненных в полевых условиях [10, 11].
Экспериментальные исследования сопротивления бипирамидальных свай выполнены на двух опытных площадках. Первая площадка представлена лессовидным суглинком (модуль деформации Е = 14500 кПа и коэффициент Пуассона n = 0,35). Вторая площадка представлена лессом (модуль деформации Е = 12000 кПа и коэффициент Пуассона n = 0,38).
В экспериментальных исследованиях была поставлена задача выявить рациональные соотношения между геометрическими размерами верхней части бипирамидальной сваи (оголовка) и нижней ее части (острия). В связи с этим, были испытаны четыре типоразмера свай на первой площадке (С-1, С-2, С-3, С-4) и пять типоразмеров на второй площадке (С'-1, С'-2, С'-3, С'-4, C'-4a). Кроме того на обеих площадках были испытаны статической нагрузкой забивные оголовки (С-0 и С'-0) и призматические сваи (С-пр и С'-пр). Геометрические размеры свай представлены на рис. 3.1. Как видно из рис. 3.1 бипирамидальные сваи С-1, С-2, С-3 (первая площадка) и сваи С'-1, С'-2, С'-3 (вторая площадка) имеют одинаковые размеры верхней части (оголовка), а длина нижней части для свай С-1 и С'-1 равна 0,7 м., для свай С-2 и С'-2 — 1,2 м., для свай С-3 и С'-3 — 1,7 м.
Здесь было намечено выявить влияние длины нижней части бипирамидальной сваи на работу оголовка. Предполагалось, что уплотненный грунт при забивке нижнего конца создает условия для повышения сопротивления верхней части. И как показывают опытные данные (см. таблицу) на первой площадке сопротивления сваи С-3 (р = 394 кН), выше сопротивление сваи С-1 (р = 264 кН) в 1,49 раза, а соотношения тех же показателей для свай на второй площадке составляет — 1,33. То есть при увеличении бетона на 27% имеем большие приращения сопротивления вертикальной нагрузке. В связи с этим можно считать, что среди рассмотренных типоразмеров сваи, наиболее рациональной является бипирамидальная свая С-3 для первой площадки и С'-3 для второй площадки.
На рис. 3.2 представлены значения сопротивлений бипирамидальных свай, полученные экспериментальным путем (Рэкс) и расчетом по изложенной в разделе 2 методике (Рт). Как видно из рисунка, значительная часть теоретических данных близко расположена к прямой, проходящей через начало координат и под углом 45о к осям координат (случай идеального совпадения экспериментальных и теоретических данных). Вместе с тем для сваи С'р, Спр, Спир, С'пир — теоретические данные больше экспериментальных. То есть, для этих типов свай необходимо вводить коэффициент запаса (надежности) больше единицы. Согласно данных приведенных в таблице 3.1 (колонка 7) этот коэффициент не превышает 18% (С'пр — Рэкс/Рт = 0,814). Для свай Ср, С'-4а экспериментальные данные на 20% превышают теоретические. То есть, в этом случае сопротивление свай может быть занижено по сравнению с действительным и здесь можно использовать коэффициент надежности меньше единицы, если принять такое соотношение
gн = Рт/Рэкс, (3.1)
где gн — коэффициент надежности расчета.
Для оценки влияния продольной формы свай введены коэффициенты которые определяются по формулам:
Коэффициенты остроты сваи
<shape id="_x0000_i1077" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image115.wmz» o:><img width=«93» height=«28» src=«dopb24148.zip» v:shapes="_x0000_i1077"> (3.2)
где В — размер поперечного сечения сваи в голове;
Vсв — объем погруженной части сваи;
Коэффициенты полноты сваи
<shape id="_x0000_i1078" type="#_x0000_t75" o:ole="" fillcolor=«window»><imagedata src=«1.files/image117.wmz» o:><img width=«75» height=«48» src=«dopb24149.zip» v:shapes="_x0000_i1078"> (3.3)
где L — длина заглубленной (погруженной) части сваи.
В данной работе выполнено исследование влияния коэффициента yв на сопротивление бипирамидальных, пирамидальных и призматических свай (Рт), которое определено теоретически. При этом подразумевалось, что теоретическое значение, как показывает ранее выполненный здесь анализ, отражает экспериментальные данные с точностью достаточной для практики проектирования, но имеют более плавный характер изменения по сравнению с экспериментальными данными, которые имеют разброс, обусловленный методикой испытаний (измерение осадки, нагрузка), процессом забивки, изготовления свай.
продолжение
--PAGE_BREAK--
еще рефераты
Еще работы по строительству
Реферат по строительству
Исследование несущей способности свай по результатам динамических испытаний в водонасыщенных
25 Июня 2015
Реферат по строительству
Производство линейных конструкций свай
1 Сентября 2013
Реферат по строительству
Проектирование устройства буронабивных свай
1 Сентября 2013
Реферат по строительству
Газоснабжение жилого дома
1 Сентября 2013