КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ ОСАДОК ЖЕСТКОГО ШТАМПА НА ОДНОРОДНОМ ВЕСОМОМ ОСНОВАНИИ

Аннотация


Приведены результаты компьютерного моделирования процесса осадок жесткого штампа, в процессе которого учтены различные факторы, влияющие на деформируемость однородного основания. В результате установлено, что величина отношения модулей деформации материала фундамента и грунта существенно влияет на характер распределения вертикальных перемещений основания под подошвой штампа. Чем это отношение больше, тем равномернее его осадка. В первом приближении графическая зависимость S / H = f (Е о / Е ш ) может быть аппроксимирована кусочно-линейной функцией. При Е o / E ш £ 0,01 эпюра вертикальных перемещений, построенная на оси симметрии штампа, приобретает некоторую нелинейность. Осадки штампа при всех прочих равных условиях существенным образом зависят от его ширины. Эта зависимость аппроксимирована полиномом второй степени, погрешность такой аппроксимации не превышает 5 %. Зависимость осадки штампа от интенсивности внешней нагрузки оказалась линейной при всех численных значениях переменных расчетных параметров. Полученные результаты можно использовать при условии, что в каждом конкретном случае глубина развития пластических областей под краями фундамента, вычисленная на основании условия пластичности Кулона, не превышает одной четверти ширины фундамента. Все вычисления проведены при помощи компьютерной программы, имеющей государственную регистрацию, в которой формализован метод конечных элементов в упругой постановке. Размеры расчетной схемы и степень ее дискретизации выбираются таким образом, что дальнейшее их увеличение не оказывает влияния на величины напряжений и перемещений в точках расчетной схемы, расположенных в активной зоне фундамента. Расчетная схема разбита на 30 000 одинаковых треугольных элементов, сопряженных в 15 251 узлах, ширина матрицы жесткости системы - 206.

Полный текст

Введение Существенные статические нагрузки от инженерных сооружений обусловливают возникновение дополнительных напряжений в грунтовом массиве основания, что определяет его деформации и, как следствие, вертикальные и горизонтальные перемещения и крены сооружения. Существует и в настоящее время широко применяются, большое количество методов расчета осадок: Ф. Шлейхера [1] - для бесконечного полупространства, находящегося под действием прямоугольной загруженной площадки, метод послойного суммирования Н.Н. Маслова, Н.А. Цытовича, С.Б. Ухова [2-4], метод эквивалентного слоя [3], метод линейно-деформируемого слоя К.Е. Егорова [5; 6] и др. Базой практически для всех методов является решение Ж.В. Буссинеска [7; 8] задачи о распределении напряжений в упругом полупространстве от действия вертикальной сосредоточенной силы на граничной горизонтальной плоскости (рис. 1, а), при этом расчетная схема сводится к схемам, приведенным на рис. 1, б, в. Анализ формул для вычисления осадок говорит о том, что они учитывают вес слоя грунта, расположенного выше подошвы фундамента (обратной засыпки), вертикальные напряжения σz, возникающие от нагрузки, действующей на фундамент, модули деформации E и коэффициенты Пуассона (бокового давления ξ) грунта основания и материала фундамента, форму и размер фундамента, которые входят в том числе и в выражение для гибкости последнего, записываемое следующим образом [9]: , (1) где - соответственно коэффициенты Пуассона и модули деформации материала сваи и грунта основания; b и l - соответственно ширина фундамента и его полудлина; J - момент инерции поперечного сечения фундамента. а б в Рис. 1. Расчетная схема к задаче Буссинеска (а), имитация заглубленного прямоугольного фундамента (б) и эпюры вертикальных напряжений σz от собственного веса грунта и внешней нагрузки для определения величины сжимаемой толщи (в) Fig. 1. Calculation scheme for the Boussinesq problem (a), simulation of a buried rectangular foundation (б) and diagrams of vertical stresses σz from the own weight of the soil and external load to determine the value of the compressible thickness (в) В то же время грунт ниже подошвы фундамента считается невесомым, наличие грунта засыпки имитируется вертикальной равномерно распределенной нагрузкой интенсивности γН (см. рис. 1, б), учет жесткости фундамента значительно затрудняют расчеты в связи с большим количеством параметров в формуле (1) и, наконец, нет возможности непосредственно учитывать перераспределение давлений на основание и выравнивание осадок за счет жесткости надфундаментного строения (рис. 2), которые целесообразно учитывать в совместном расчете системы «здание (сооружение) - фундамент - основание» [10]. Факт перераспределения давления (напряжений) за счет жесткости надфундаментной части сооружения наглядно иллюстрируется результатами исследования моделей оснований фундаментов, выполненных из оптически активных материалов, поляризационно-оптическим методом [11], которые приведены на рис. 2. В настоящей работе приводятся некоторые результаты численного исследования зависимости осадок жесткого штампа, нагруженного равномерно распределенной нагрузкой, являющегося начальным этапом работы по выявлению общих закономерностей осадок жесткого фундамента и надфундаментной конструкции как единого целого. а б Рис. 2. Картины изохром, характеризующие распределение напряжений в моделях основания из оптически активного материала под жестким (а) и гибким (б) штампами (цитируется по работе [10]) Fig. 2. Isochrome patterns characterizing the distribution of stresses in base models made of optically active material under (a) rigid and (б) flexible dies (cited from [10]) Компьютерное моделирование Параметры моделирования При проведении компьютерного моделирования осадок будем использовать компьютерную программу [12], в которой для определения напряжений и перемещений формализован метод конечных элементов (МКЭ), реализующий линейно-упругую модель грунтовой среды. На рис. 3 изображен фрагмент расчетной схемы МКЭ, которая ввиду симметричности ленточного или плитного фундамента включает в себя лишь его половину. При расчете задаются тривиальные граничные условия: узлам, расположенным на вертикальных границах расчетной схемы (РС), запрещены горизонтальные перемещения, а вертикальные не ограничены. Узлам, расположенным на нижней горизонтальной границе РС, напротив, запрещены вертикальные перемещения, а горизонтальные не ограничены; нижние угловые узлы РС неподвижны. По контакту «штамп - основание» принято условие полного прилипания грунта. Считается, что толщина штампа равна глубине его заложения. Размеры расчетной схемы и степень ее дискретизации выбираются таким образом, что дальнейшее их увеличение не оказывает влияние на величины напряжений и перемещений в точках РС, расположенных в активной зоне фундамента. В нашем случае расчетная схема разбита на 30 000 одинаковых треугольных элементов, сопряженных в 15 251 узлах, ширина матрицы жесткости системы - 206. Расчетная схема составлена таким образом, что за условную единицу измерения длины принята глубина заложения фундамента Н, т.е. все геометрические размеры РС выражены в ее долях. За условную единицу объемной силы принят объемный вес грунта γо и величина интенсивности внешней нагрузки задается в долях γоН. Результаты расчетов при этом, в частности значения напряжений и перемещений, будут измеряться в безразмерных долях γоН и Н соответственно, что весьма удобно при анализе этих результатов. Перейдем теперь к определению переменных расчетных параметров и диапазонам их изменения. Рис. 3. Фрагмент схемы МКЭ для расчета осадок штампа Fig. 3. Fragment of the FEM scheme for calculating the plate settlement Из формулы (1) следует, что (2) где - некоторый безразмерный коэффициент, причем, , так как [13; 14]. Из этого соотношения следует, что . Учитывая, что для бетона [14], а среднее значение коэффициента бокового давления для глинистых грунтов ξо = 0,75 [15], то используя известное соотношение между ξо и , получим, что величина = 0,43, а среднее значение отношения . Такая величина отношения не оказывает сколько-нибудь заметного влияния на значения напряжений, а значит, и перемещений в точках активной зоны фундамента [16]. Учитывая выражение (2), можно утверждать, что при моделировании осадок ленточного и плитного фундаментов в качестве одного из параметров можно использовать величину отношения модулей деформации материала фундамента и грунта основания . Так как величины k и t индивидуальны для каждого фундамента и соответствующих грунтовых условий и величина их отношения может изменяться в широком диапазоне, то при проведении компьютерного моделирования диапазон изменения величины тоже должен быть достаточно широк. Причем - жесткость фундамента не может быть меньше жесткости основания. Условие означает, что нагрузка приложена непосредственно к поверхности грунта основания (фундамент как таковой отсутствует). Кроме отношения на величину осадок плитного фундамента будут оказывать влияние величина отношения ширины фундамента к его толщине (в нашем случае глубине заложения) и величина интенсивности внешней нагрузки, которая задается в безразмерных единицах . При дальнейшем изложении материала величину будем обозначать Для того чтобы наверняка охватить весь диапазон изменения вышеперечисленных параметров, при проведении компьютерного моделирования, принято, что им последовательно присваиваются следующие значения: 5; 10; 20; 50, и . Таким образом, выполнено 240 вычислительных операций, результаты которых приводятся ниже. Под осадкой штампа будем понимать вертикальное перемещение точки «0» подошвы штампа, лежащей на вертикальной оси симметрии фундамента (см. рис. 1, б). Анализ результатов компьютерного моделирования Влияние величины отношения Ешт / Eo на осадки штампа при фиксированных значениях параметров 2b / H и q / γoH На рис. 4 изображены картины изолиний вертикальных перемещений (осадок) в однородном основании под жестким штампом шириной 2b / H = 10 и q / γoH = 12 при величинах отношения модулей деформации материалов основания и штампа Eшт / Еo = 1,0; 10; 100; 1000, а на рис. 5 - эпюры этих перемещений (осадок), построенные при тех же условиях на вертикальной оси, проходящих через середину штампа. Рис. 4. Изолинии вертикальных перемещений (осадок) в однородном основании под жестким штампом шириной 2b / H = 10 и q / γoH = 12 при величинах отношения модулей деформации материалов основания и штампа Eшт / Еo = 1,0; 10; 100; 1000 - а, б, в, г соответственно Fig. 4. Isolines of vertical displacements (settlement) in a homogeneous base under a rigid stamp with a width of 2b / H = 10 and q / γoH = 12 at the values of the ratio of the deformation moduli of the base materials and the stamp Eшт / Еo = 1,0; 10; 100; 1000 - a, б, в, г respectively Анализ изображений на рис. 4 показывает, что с ростом численного значения величины отношения Eшт / Еo кривизна изолиний перемещений уменьшается, происходит постепенное выравнивание осадок под подошвой фундамента. Чем больше Eшт / Еo, тем равномернее осадки штампа, а их численные значения увеличиаются пропорционально росту численного значения этого отношения. Заметим, что данное обстоятельство отмечено и автором работы [10]. Анализ эпюр осадок, построенных на вертикальной оси симметрии фундамента, говорит о том, что изменение численных значений перемещений по глубине происходит по линейному закону, однако при Eшт / Еo = 1000 даже визуально можно заметить некоторую кривизну образующей линии эпюры. Рис. 5. Эпюры вертикальных перемещений (осадок) в однородном основании под жестким штампом шириной 2b / H = 10 и q / γoH = 12 при величинах отношения модулей деформации материалов основания и штампа; 1000 - а, б, в, г соответственно, построенные на вертикальной оси, проходящих через середину штампа Fig. 5. Diagrams of vertical displacements (settlement) in a homogeneous base under a rigid stamp with a width of 2b / H = 10 and q / γoH = 12 at the values of the ratio of the moduli of deformation of the materials of the base and the stamp; 1000 - a, б, в, г respectively, built on a vertical axis passing through the middle of the stamp На рис. 6 изображены графические зависимости вида S / H = f (Ео / Еш) при различных значениях величины интенсивности равномерно распределенной нагрузки q / γoH, приложенной к поверхности штампа, и фиксированных значениях величины отношения ширины штампа к глубине заложения (толщине) 2b / H. Нам представляется, что в первом приближении возможна линейная аппроксимация зависимостей вида S / H = f (Ео / Еш) на отрезках Ео / Еш Î [0,001-0,01] и Ео / Еш Î [0,01-0,05]. Таким образом, все графические аппроксимации можно представить в виде кусочно-линейных функций, имеющих перелом в точках с абсциссой Ео / Eш » 0,01. Однако из рис. 6, д, е, где приведены результаты построения этих графических зависимостей в оболочках MathCad и Excel, видно, что на отрезке Ео / Еш Î [0,01-0,05] явно наблюдается нелинейность. Поэтому при проведении дальнейших исследований необходимо провести оценку погрешности такой аппроксимации и в случае необходимости предложить другой ее вид. Рис. 6. Графические зависимости вида S / H = f (Ео / Еш) при различных значениях величины q / γoH (равные 0, 3, 12), и фиксированных значениях величины отношения ширины штампа к глубине заложения (2b / H) (а-г); графические аппроксимации тех же зависимостей в оболочках MathCad и Excel (д, е) Fig. 6. Graphic dependences of the form S / H = f (Ео / Еш) for different values of q / γoH (to 0, 3, 12), and fixed values of the ratio of the width of the stamp to the depth of placement (2b / H) (а-г); graphical approximations of the same dependencies in MathCad and Excel shells (д, е) Зависимость осадок штампа от его ширины Как видно из рис. 7, зависимость осадки штампа от его ширины имеет весьма сложный характер, и ее не просто анализировать. На первом этапе проведена аппроксимация этих зависимостей полиномом второй степени, пример такой аппроксимации приведен на рис. 7, д. Установлено, что погрешность такой аппроксимации не превышает 5 %. Отметим, что попытка аппроксимации графической зависимости S / H = f (2b / H) линейной либо кусочно-линейной функцией дают весьма существенную ошибку аппроксимации (более 50 %). Зависимость осадок штампа от интенсивности нагрузки Данные зависимости, как и следовало ожидать, являются линейными при всех численных значениях переменных расчетных параметров, рассмотренных в настоящей работе. Рис. 7. Графические зависимости вида S / H = f (2b / H) при различных значениях величины (Еш / Ео) (равные 1, 5, 10, 20, 50) и фиксированных значениях интенсивности внешней равномерно распределенной вертикальной нагрузки q = nγН, приложенной к штампу (а-г); пример полиномиальной аппроксимации (д) Fig. 7. Graphical dependencies of the form S / H = f (2b / H) for different values of the quantity (Eш / Eo) (to 1, 5, 10, 20, 50) and fixed values of the intensity of the external uniformly distributed vertical load q = nγН applied to the stamp (a-г); polynomial approximation example (д) Рис. 8. Графические зависимости вида S / H = f (q / γH) при различных значениях величины отношения модулей деформации материалов штампа и основания (Еш / Ео) (равные 1, 5, 10, 20, 50) и фиксированных значениях величины отношения ширины штампа к глубине заложения: а - 2b / H = 2, б - 2b / H = 4, в - 2b / H = 10, г - 2b / H = 20 Fig. 8. Graphical dependences of the form S / H = f (q / γH) at different values of the ratio of the deformation moduli of the stamp and base materials (Eш / Eo) (to 1, 5, 10, 20, 50) and fixed values of the ratio of the width of the stamp to the insertion depth: а - 2b / H = 2, б - 2b / H = 4, в - 2b / H = 10, г - 2b / H = 20 Возможность использования полученных результатов Графические аппроксимации полученных зависимостей могут быть использованы для предварительной оценки осадки ленточных и плитных фундаментов. Применяя различные интерполяционные методы и последовательно переходя от одной пары переменных расчетных параметров к другой, можно получить значение упругой осадки фундамента при условии, что его размеры и упругие свойства его материала соответствуют параметрам моделирования. На рис. 9 приведены изображения областей пластических деформаций, построенных на основании условия пластичности Кулона [18], в однородном основании жесткого штампа при различных значениях отношения модулей деформации материалов основания и штампа (см. рис. 9, а-в). Здесь же изображена графическая зависимость относительной глубины развития областей пластических деформаций под краями штампа единичной толщины и относительной ширины 2b / H = 10 при интенсивности внешнего воздействия q / γoH = 3 от величины отношения Eo / Eш в однородном основании, сложенном суглинком (c = 47 кПa, φ = 26o, Ео = 39 МПа, n = 0,28) (г). Рис. 9. Области пластических деформаций, построенные на основании условия пластичности Кулона, в однородном основании жесткого штампа при различных значениях отношения модулей деформации материалов основания и штампа (а-в); г - зависимость относительной глубины развития областей пластических деформаций под краями штампа единичной толщины (Н = 1) и относительной ширины 2b / H = 10 при интенсивности внешнего воздействия q / γoH = 3 от величины отношения Eo / Eш в однородном основании, сложенном суглинком (c = 47 кПa, φ = 26o, Ео = 39 МПа, n = 0,28) Fig. 9. Areas of plastic deformations, built on the basis of the Coulomb plasticity condition, in a homogeneous base of a rigid stamp at different values of the ratio of the deformation moduli of the materials of the base and stamp (a-в); г - dependence of the relative depth of development of areas of plastic deformations under the edges of a stamp of unit thickness (Н = 1) and relative width 2b / H = 10 at the intensity of external action q / γoH = 3 on the value of the ratio Eo / Eш in a homogeneous base composed of loam (c = 47 kPa, φ = 26o, Ео = 39 МPа, n = 0,28) Для обоснования возможности использования «упругих» осадок для дальнейших расчетов, необходимо для каждого конкретного варианта расчета (конкретных величин q, с, φ, Ео, n, 2b / H) построить графическую зависимость вида ΔZ / H = f (Eo / Eш). Если глубина развития пластических деформаций окажется меньше четверти ширины фундамента (ΔZ < 0,5b), то полученный результат можно считать достоверным [19]. Выводы В процессе компьютерного моделирования процесса осадок плитного фундамента рассмотрено влияние различных факторов на деформируемость однородного основания. Установлено: 1. Величина отношения модулей деформации материала фундамента и грунта существенно влияет на характер распределения вертикальных перемещений основания под подошвой фундамента. Чем это отношение больше, тем равномернее осадка фундамента. В первом приближении графическая зависимость S / H = f (Ео / Еш) может быть аппроксимирована кусочно-линейной функцией. При Еo / Eш £ 0,01 эпюра вертикальных перемещений, построенная на оси симметрии штампа, приобретает некоторую нелинейность. 2. Осадки штампа при всех прочих равных условиях существенным образом зависят от его ширины. На начальном этапе исследования эта зависимость аппроксимирована полиномом второй степени, погрешность такой аппроксимации не превышает 5 %. 3. Зависимость осадки штампа от интенсивности внешней нагрузки оказалась линейной при всех численных значениях переменных расчетных параметров. 4. Полученные результаты можно использовать при условии, что в каждом конкретном случае глубина развития пластических областей под краями фундамента, вычисленная на основании условия пластичности Кулона, не превышает одной четверти ширины фундамента.

Об авторах

Л. А Бартоломей

Тюменский индустриальный университет

О. А Богомолова

Волгоградский государственный технический университет

В. Д Гейдт

Тюменский индустриальный университет

А. В Гейдт

Тюменский индустриальный университет

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

  1. Schleicher F. Zur Theorie des Baurgundes. - Der Bauinginieur. Berlin, 1926. - H. 48-49.
  2. Маслов Н.Н. Длительная устойчивость и деформации смещения подпорных сооружений. - М.: Энергия, 1968.
  3. Цытович Н.А. Механика грунтов. - М.: Высшая школа, 1983. - 228 с.
  4. Механика грунтов, основания и фундаменты / С.Б. Ухов, В.В. Семенов, В.В. Знаменский [и др.]. - М.: Высшая школа, 2007. - 566 с.
  5. Егоров К.Е. К вопросу о деформации основания конечной толщины // Механика грунтов: сб. тр. НИИОСП. - № 34. - М.: Госстройиздат, 1958.
  6. Егоров К.Е. К вопросу расчета основания под фундаментом с подошвой кольцевой формы // Механика грунтов: сб. тр. НИИОСП. - № 34. - М.: Госстройиздат, 1958.
  7. Boussinesq J. Applications der patentiels a l’etude de l’eqiuilibre et du mouvements de Solidos, elastiques. - Paris, 1885.
  8. Boussinesq J. Essai sur la theorie des eaux courantes. - Paris, 1877. - 680 с.
  9. Основания, фундаменты и подземные сооружения. Справочник проектировщика / М.И. Горбунов-Посадов, В.А. Ильичев, В.И. Крутов и др.; под общ. ред. Е.А. Сорочана и Ю.Г. Трофименкова. - М.: Стройиздат, 1985. - 480 с.
  10. Кушнер С.Г. Расчет деформаций оснований зданий и сооружений. - Запорожье, 2008. - 496 с.
  11. Александров А.Я., Ахметзянов М.Х. Поляризационно-оптические методы механики деформируемого тела. - М.: Наука, 1973. - 576 с.
  12. FEA: свидетельство о государственной регистрации программы для ЭВМ № 2015617889 / О.А. Богомолова и др. Зарегистрировано 23 июля 2015 г.
  13. Каган А.А. Расчетные показатели физико-механических свойств грунтов. - Ленинград: Стройиздат, 1973. - 144 с.
  14. СНиП 2.02.01.-83. Бетонные и железобетонные конструкции / Госстрой СССР. - М.: Стройиздат, 1985.
  15. Месчян С.Р. Начальная и длительная прочность глинистых грунтов. - М.: Недра, 1978. - 207 с.
  16. Амусин Б.З., Фадеев А.Б. Метод конечных элементов при решении задач горной геомеханики. - М.: Недра, 1975. - 144 с.
  17. Сопоставление результатов расчета несущей способности двухслойного основания заглубленного ленточного фундамента различными способами / О.А. Богомолова, А.И. Вайнгольц, О.В. Ермаков [и др.] // Вестник Пермского национального исследовательского политехнического университета. Строительство и архитектура. - 2014. - № 2. - С. 106-116.
  18. Coulomb C. Application des riles de maximus et minimis a quelques problemes de statique relatifs a L`architecture // Memories de savants strangers de L`Academlie des sciences de Paris. - 1773. - 233 p.
  19. Флорин В.А. Основы механики грунтов. - Л.: Госстройиздат, 1959. - Т. I. - 360 c.

Статистика

Просмотры

Аннотация - 50

PDF (Russian) - 35

Ссылки

  • Ссылки не определены.

© Бартоломей Л.А., Богомолова О.А., Гейдт В.Д., Гейдт А.В., 2022

Creative Commons License
Эта статья доступна по лицензии Creative Commons Attribution-NonCommercial 4.0 International License.

Данный сайт использует cookie-файлы

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

О куки-файлах