COMPUTER SIMULATION OF RIGID PLATE SETTLEMENT ON A HOMOGENEOUS WEIGHT BASE
- Authors: Bartolomey L.A1, Bogomolova O.A2, Geidt V.D1, Geidt A.V1
- Affiliations:
- Industrial University of Tyumen
- Volgograd State Technical University
- Issue: Vol 13, No 2 (2022)
- Pages: 5-17
- Section: ARTICLES
- URL: https://ered.pstu.ru/index.php/CG/article/view/3511
- DOI: https://doi.org/10.15593/2224-9826/2022.2.01
- Cite item
Abstract
The paper presents the results of computer simulation of the process of hard stamp upsetting, during which various factors affecting the deformability of a homogeneous base are taken into account. As a result, it was found that the value of the ratio of the deformation moduli of the foundation material and soil significantly affects the nature of the distribution of vertical displacements of the base under the sole of the stamp. The larger this ratio, the more uniform its settlement. In the first approximation, the graphical dependence S / H = f (Е о / Е ш ) can be approximated by a piecewise linear function. At Е o / E ш £ 0,01 the diagram of vertical displacements, built on the axis of symmetry of the stamp, acquires some non-linearity. Draft of the stamp, all other things being equal, significantly depend on its width. This dependence is approximated by a polynomial of the second degree; the error of such an approximation does not exceed 5 %. The dependence of the die settling on the intensity of the external load turned out to be linear for all numerical values of the variable design parameters. The results obtained can be used provided that in each specific case the depth of development of plastic regions under the edges of the foundation, calculated on the basis of the Coulomb plasticity condition, does not exceed one quarter of the width of the foundation. All calculations were carried out using a state-registered computer program that formalizes the finite element method in an elastic formulation. The dimensions of the design scheme and the degree of its discretization are chosen in such a way that their further increase does not affect the magnitude of stresses and displacements at the points of the design scheme located in the active zone of the foundation. The design scheme is divided into 30 000 identical triangular elements conjugated at 15 251 nodes, the width of the system stiffness matrix is - 206.
Full Text
Введение Существенные статические нагрузки от инженерных сооружений обусловливают возникновение дополнительных напряжений в грунтовом массиве основания, что определяет его деформации и, как следствие, вертикальные и горизонтальные перемещения и крены сооружения. Существует и в настоящее время широко применяются, большое количество методов расчета осадок: Ф. Шлейхера [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. Полученные результаты можно использовать при условии, что в каждом конкретном случае глубина развития пластических областей под краями фундамента, вычисленная на основании условия пластичности Кулона, не превышает одной четверти ширины фундамента.About the authors
L. A Bartolomey
Industrial University of Tyumen
O. A Bogomolova
Volgograd State Technical University
V. D Geidt
Industrial University of Tyumen
A. V Geidt
Industrial University of Tyumen
References
- Schleicher F. Zur Theorie des Baurgundes. - Der Bauinginieur. Berlin, 1926. - H. 48-49.
- Маслов Н.Н. Длительная устойчивость и деформации смещения подпорных сооружений. - М.: Энергия, 1968.
- Цытович Н.А. Механика грунтов. - М.: Высшая школа, 1983. - 228 с.
- Механика грунтов, основания и фундаменты / С.Б. Ухов, В.В. Семенов, В.В. Знаменский [и др.]. - М.: Высшая школа, 2007. - 566 с.
- Егоров К.Е. К вопросу о деформации основания конечной толщины // Механика грунтов: сб. тр. НИИОСП. - № 34. - М.: Госстройиздат, 1958.
- Егоров К.Е. К вопросу расчета основания под фундаментом с подошвой кольцевой формы // Механика грунтов: сб. тр. НИИОСП. - № 34. - М.: Госстройиздат, 1958.
- Boussinesq J. Applications der patentiels a l’etude de l’eqiuilibre et du mouvements de Solidos, elastiques. - Paris, 1885.
- Boussinesq J. Essai sur la theorie des eaux courantes. - Paris, 1877. - 680 с.
- Основания, фундаменты и подземные сооружения. Справочник проектировщика / М.И. Горбунов-Посадов, В.А. Ильичев, В.И. Крутов и др.; под общ. ред. Е.А. Сорочана и Ю.Г. Трофименкова. - М.: Стройиздат, 1985. - 480 с.
- Кушнер С.Г. Расчет деформаций оснований зданий и сооружений. - Запорожье, 2008. - 496 с.
- Александров А.Я., Ахметзянов М.Х. Поляризационно-оптические методы механики деформируемого тела. - М.: Наука, 1973. - 576 с.
- FEA: свидетельство о государственной регистрации программы для ЭВМ № 2015617889 / О.А. Богомолова и др. Зарегистрировано 23 июля 2015 г.
- Каган А.А. Расчетные показатели физико-механических свойств грунтов. - Ленинград: Стройиздат, 1973. - 144 с.
- СНиП 2.02.01.-83. Бетонные и железобетонные конструкции / Госстрой СССР. - М.: Стройиздат, 1985.
- Месчян С.Р. Начальная и длительная прочность глинистых грунтов. - М.: Недра, 1978. - 207 с.
- Амусин Б.З., Фадеев А.Б. Метод конечных элементов при решении задач горной геомеханики. - М.: Недра, 1975. - 144 с.
- Сопоставление результатов расчета несущей способности двухслойного основания заглубленного ленточного фундамента различными способами / О.А. Богомолова, А.И. Вайнгольц, О.В. Ермаков [и др.] // Вестник Пермского национального исследовательского политехнического университета. Строительство и архитектура. - 2014. - № 2. - С. 106-116.
- 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.
- Флорин В.А. Основы механики грунтов. - Л.: Госстройиздат, 1959. - Т. I. - 360 c.
Statistics
Views
Abstract - 87
PDF (Russian) - 51
Refbacks
- There are currently no refbacks.