CC BY
18Научно-технический и производственный журнал
-------ЖИЛИЩНОЕ ---
СТРОИТЕЛЬСТВО
Подземное строительство
УДК 624.044.2
К.Г. ШАШКИН, канд. техн. наук, ООО «ПИГеореконструкция» (Санкт-Петербург); В.А. ШАШКИН, инженер, Санкт-Петербургский государственный архитектурно-строительный университет
Решение задач предельного равновесия с использованием метода конечных элементов
Рассмотрены вопросы перехода задач деформирования основания в задачи предельного равновесия. Приведено сравнение задач, решенных с применением метода конечных элементов с численными и аналитическими решениями теории предельного равновесия, выделены основные факторы, влияющие на вид графика «нагрузка-осадка» для штампа на идеально упругопластическом основании.
Ключевые слова: напряженно-деформированное состояние, предельное равновесие, метод конечных элементов.
Решение задач предельного равновесия для определения несущей способности различного типа фундаментов имеет длинную историю. Точные решения для случая невесомого грунта были получены в начале XX в. Л. Прандт-лем (1921 г.) и Г. Рейсснером (1924 г.). Для практического случая основания с ненулевым удельным весом, углом внутреннего трения и сцеплением К. Терцаги (1943 г.) сформулировал классическую трехчленную формулу:
Ри, кПа 40000
35000 30000 25000 20000 15000 10000 5000
0
q
0 1 2 3 4 5 6
Рис. 1. Зависимостьри(4) при ф=45о и с=0по аппроксимации [3]и по точному решению [4]
N
400
У
350 300 250 200 150 100 50 0
50
(1)
в которой коэффициенты N и N определяются в зависимости от значений угла внутреннего трения. По решению Л. Прандтля:
4 1 —вшф
Определение коэффициента N оказывается существенно более сложным. Методику решения данной задачи для плоской постановки сформулировал В.В. Соколовский (1942 г.). В.Г. Березанцев обобщил данную задачу для случая осевой симметрии. Вопросам оценки несущей способности фундамента с учетом собственного веса грунта посвящены работы Г. Мейергофа (1951 г.), Й. Хансена (1970 г.), А. Весича (1975 г.). Основополагающие теоретические решения были получены Д. Друкером (1952 г.). Для случая идеально упругопластического материала с ассоциированным законом пластического течения им были сформулированы теоремы о верхней и нижней границах, в пределах которых должно находиться значение предельной нагрузки.
Рис. 2. Сравнение величин коэффициента N, вычисленного по различным решениям
Рис. 3. Сравнение положения предельных поверхностей: 1 — по условию Кулона—Мора; 2 — по выражению дит=Мр (4); 3 — по модели Друкера—Прагера
а
а
Ф
а
Подземное строительство
ц м .1
Научно-технический и производственный журнал
100
300
500
700
900
Р, кПа 1100
0 L-0
10
ф
30
40
50
Рис. 6. Сравнение областей сдвига (максимальные перемещения) при решении упругопластических задач и по теории предельного равновесия при ф =20°
S, м
Рис. 4. Сравнение графиков «нагрузка-осадка» при использовании условия (4), модели Друкера—Прагера (описанный конус вокруг пирамиды Кулона) и классической модели с критерием Кулона—Мора
100 90 80 70 60 50 40 30 20 10
d, м
Рис. 7. Зависимость предельного давления от глубины заложения фундамента й шириной Ь=2м при с=20кПа; ф=0о
Рис. 5. Сравнение коэффициента Ny по результатам решения упругопластических задач в программе FEM models и по решениям С. Мартина, С. Слоана и В.Г. Федоровского
В последние годы решением задач предельного равновесия успешно занимались А.М. Караулов (2002 г.) [1], В.Г. Федоровский (2005 г.) [2], К.В. Королев (2010 г.) [3] и др.
Следует отметить, что коэффициенты N, приведенные в отечественных нормативных документах и в [3], отражают решение задачи при наличии пригрузки и угла внутреннего трения при условии:
,= g + c-c&|>>1 уЬ
При приближении q' к нулю зависимость p(q) приобретает нелинейный характер, в результате чего, как отмечается в [2, 3], формула (1) становится принципиально неверной (рис. 1). Поэтому коэффициенты Ny, приведенные в [3] и в отечественных нормах, отвечают линейному участку зависимости pu(q) и несправедливы для случая штампа на поверхности или на небольшой глубине в сыпучей среде.
Задача определения несущей способности штампа на поверхности сыпучей среды (при q1, стремящейся к нулю) долгое время являлась недоступной для точного теоретического решения. Ввиду сложностей, связанных с учетом соб-
Рис. 8. Характер областей сдвига (максимальные перемещения на стадии потери устойчивости) при различном заглублении фундамента шириной Ь=2м при с=20кПа; ф=0о: а — 1 м; б — 5 м
ственного веса, получено большое число примерных значений коэффициента отличия между которыми зачастую значительны, особенно для углов внутреннего трения, превышающих 30о. С. Слоаном и др. [5] путем решения задачи предельного равновесия методом конечных элементов были получены значения коэффициента N по нижнему и верхнему пределам, отличающимся не более чем на 3,4% друг от друга. Величины полученные Мартином [4], находятся в границах, полученных Слоаном. В статье В.Г. Федоровского [2] приведено точное решение более общей задачи с учетом наклонного действия нагрузки на штамп. В частном случае для идеально шероховатого штампа и вертикальной нагрузки коэффициенты N практически точно совпадают с результатами, приведенными в статье С. Мартина, и лежат в границах, приведенных в статье С. Слоана.
Таким образом, анализ данных работ позволяет утверждать, что для коэффициента N в случае фундамента на по-
Научно-технический и производственный журнал
Подземное строительство
Теория Угол внутреннего трения ф
5o 10o 15o 20o 25o 30o 35o 40o 45o
СП 50-101-2004* 0,4 1,2 2,7 5,76 11,74 24,78 55 132,02 355,22
Eurocode 7 (Vesic) 0,1 0,52 1,58 3,94 9 20,1 45,24 106,04 267,74
Слоан верхняя граница 0,12 0,456 1,238 2,961 6,738 15,237 35,649 88,39 240,88
нижняя граница 0,115 0,434 1,178 2,822 6,431 14,567 33,951 83,327 224,945
Мартин, Федоровский 0,11 0,43 1,18 2,84 6,49 14,76 34,48 85,58 234,3
* Коэффициент по СП 50-101-2004 приведен в соответствие с формулой (1).
верхности сыпучего основания наконец найдено точное решение. Авторы [2, 4, 5] различным путем пришли практически к одинаковым значениям коэффициента которые существенно меньше значений, приведенных в отечественных нормах. Причиной расхождения является то, что для специфической задачи о незаглубленном штампе в несвязной среде коэффициент N строго говоря, не является независимым коэффициентом в формуле (1), а сам является функцией от величины пригрузки, поскольку функция ра(д") становится нелинейной (рис. 1).
Коэффициенты по работам С. Мартина, С. Слоана, В.Г. Федоровского в сравнении с данными отечественных и зарубежных норм приведены в таблице и на рис. 2. Учитывая наличие точных решений для N при д'<1, представляется актуальным ввести ограничение на применение обычных коэффициентов N при д'<1. Для этого участка целесообразно использовать более точную аппроксимацию нелинейного участка зависимости ра(д") по сравнению с формулой К. Терцаги (1).
Факторы, влияющие на вид графика «нагрузка-осадка» при решении задач о вдавливании штампа:
- модуль деформации грунта и размер расчетной области;
- величина коэффициента Пуассона;
- вид деформирования (дренированная или недрениро-ванная постановка задачи);
- вид пластического течения: ассоциированный или неас-социированный, например равнообъемный, закон пластического течения;
- начальное напряженное состояние, обычно характеризуемое коэффициентом бокового давления К0 (те программы, которые не позволяют регулировать коэффициент К0 на стадии природного напряженного состояния, непригодны для решения геотехнических задач);
- особенности поверхности текучести.
На последнем факторе стоит остановиться особо. Не секрет, что в различных конечно-элементных программах используются модели с разными поверхностями текучести, близкими к классическому критерию Кулона-Мора. Например, в модели Друкера-Прагера вместо шестигранной пирамиды Кулона-Мора используется коническая поверхность (рис. 3).
Результаты решения конечно-элементных задач показывают, что коническая поверхность, описанная вокруг пирамиды Кулона-Мора, приводит к существенно большей несущей способности, чем по классическому критерию Кулона-Мора (рис. 4). С другой стороны, коническая поверхность, вписанная в пирамиду Кулона-Мора, не соответствует критерию прочности при ст2=ст3, т. е. при обычном напряженном состоянии в трехосных испытаниях. В результате пользователь модели Друкера-Прагера оказывается перед дилеммой, чаще всего им не осознаваемой: либо существенно завысить несущую способность основания по сравнению с классическими решениями, либо пренебречь возможностью моделирования трехосных испытаний.
Другим примером отклонения от поверхности Кулона-Мора являются модели типа Cam Clay. Данные модели традиционно строятся в координатах р=(ст1+ст2+ст3)/3 и q=(a1 — ст3)/2. В этих осях можно провести линию:
,,пД q,m=Mp. (4)
м= 3sin<|)
При З-вшф эта линия соответствует критерию Кулона-Мора, но при одном условии: ст2=ст3. В целом же критерий (4) дает пирамиду, несколько отличающуюся от пирамиды Кулона-Мора (рис. 5). В результате решения идеально упругопластической задачи с критерием (4) мы получаем существенно иной о график «нагрузка-осадка» и другую несущую способность фундамента (рис. 6). Поэтому при использовании моделей типа Cam Clay с критерием (4) следует иметь в виду, что несущая способность фундамента
Р, кПа
)0
S, м
Рис. 9. Зависимость осадки S от нагрузки р при различном заглублении фундамента
100
80
60
40
20
0
Рис. 10. Зависимость предельного давления от эксцентриситета приложения нагрузки е для фундамента шириной Ь=2 м при с=20кПа, ф=0о
Ри, кПа
е
Подземное строительство
Ц M .1
Научно-технический и производственный журнал
е=0,1 м
е=0,4 м
J
е=0,6 м
е=1 м
Рис. 11. Характер областей сдвига (максимальные перемещения на стадии потери устойчивости) при различном эксцентриситете e для фундамента шириной b=2м при c=20кПа; ф=0о
получается несколько выше, чем по решениям теории предельного равновесия.
Таким образом, следует обратить внимание специалистов на множество факторов, от которых зависит график «нагрузка-осадка», даже при использовании простейшей идеально упругопластической модели. При использовании более сложных моделей количество факторов многократно увеличивается.
Сравнение коэффициента N по результатам решения упругопластических задач и по теории предельного равновесия. Проведенные исследования позволили установить, что проблемы с получением предельной несущей способности, описанные в [3], связаны не с недостатками МКЭ, а с особенностями решения данной задачи при использовании конечных элементов 1-го порядка. Использование этих элементов для решения геотехнических задач возможно только при использовании специальной техники интегрирования вектора сил. В противном случае элементы 1-го порядка дают некорректные значения предельной нагрузки и некорректные деформации в допредельном состоянии. Кроме того, решение упругопластических задач, переходящих в задачи предельного равновесия, сопряжено с рядом дополнительных трудностей. В частности, при угле внутреннего трения выше 30о для получения корректного решения приходится существенно увеличивать точность решения задачи. На рис. 7 приведено сравнение результатов расчета коэффициентов N по программе FEM models (использование выражения (3) для элементов 1-го порядка) и по приведенным в п. 1 решениям теории предельного равновесия. Как видно из сравнения, метод конечных элементов позволяет получить решения, достаточно близкие к значениям, полученным С. Слоаном, С. Мартином и В.Г. Федоровским [2, 4, 5], чего не удается достичь при использовании программы PLAXIS в связи с особенностями работы алгоритма автовыбора шага нагружения.
На рис. 8 приведено сравнение областей предельного состояния, полученных по решению упругопластической задачи методом конечных элементов и по программе С. Мартина [4], демонстрирующее их удовлетворительное совпадение.
Оценка несущей способности заглубленного фундамента и фундамента с эксцентрично приложенной нагрузкой выполнялась для идеально связной среды (с нулевым углом внутреннего трения). В соответствии с классическими решениями учет заглубления фундамента производится через добавление по поверхности нагрузки q, равной давлению грунта выше подошвы заглубленного фунда-
мента. Однако в реальности при увеличении заглубления начинает играть роль не только пригрузка основания рядом с фундаментом, но и прочность грунта выше подошвы штампа, что приводит к увеличению предельной нагрузки по сравнению с теоретическим решением (рис. 9).
Как видно из рис. 8, а, при небольшой глубине фундамента происходит выпор по форме, близкой к незаглуб-ленному фундаменту. При увеличении заглубления фундамента область выпора становится менее определенной (рис. 8, б). Кроме того, из рис. 9 видно, что при увеличении заглубления штампа увеличивается зона линейного поведения на графике «нагрузка-осадка».
Учет эксцентриситета приложения нагрузки в отечественных нормах выполняется путем уменьшения ширины рассчитываемого фундамента b'=b—2e, где e - эксцентриситет приложения нагрузки. Данный подход также вполне оправдан при небольшом эксцентриситете (рис. 10). Однако при приближении точки приложения нагрузки к краю фундамента такой подход приводит к снижению до нуля несущей способности. В то же время очевидно, что при расположении нагрузки на краю фундамента несущая способность фундамента должна быть сильно сниженной, но ненулевой. Конечно-элементное решение показывает именно такую зависимость (рис. 10). Формы потери устойчивости при эксцентричном приложении нагрузки приведены на рис. 11.
Для фундаментов на несвязном основании при небольшом заглублении подошвы представляется целесообразным внести коррективы в нормативные документы на основании точных решений, полученных С. Слоаном, С. Мартином и В.Г. Федоровским. Необходимо также внести ограничения на использование обычных (заложенных в нормативные документы) коэффициентов N при небольшом заглублении фундаментов в несвязном грунте или откорректировать традиционно применяемую трехчленную формулу и применить более точную аппроксимацию для малозаглуб-ленных фундаментов.
Решение задач предельного равновесия методом конечных элементов позволяет получить коэффициенты Ny, близкие к решениям С. Слоана, С. Мартина и В.Г. Федоровского.
Использование метода конечных элементов позволяет уточнить особенности потери устойчивости основания для ряда задач, которые имеют только приближенное решение по теории предельного равновесия, в частности для задач с заглубленным штампом, с эксцентриситетом приложения нагрузки или с критериями прочности, отличающимися от критерия Кулона-Мора.
Список литературы
1. КарауловА.М. Несущая способность основания осесиммет-ричных фундаментов. Новосибирск: СГУПС, 2002. 104 с.
2. Федоровский В.Г. Несущая способность сыпучего основания ленточного фундамента при действии наклонной внецентренной нагрузки. // Основания фундаментов и механика грунтов. 2005. № 4. С. 1-7.
3. Королев К.В. Плоская задача теории предельного равновесия грунтов. Новосибирск: СГУПС, 2010. 250 с.
4. Martin C.M. Exact bearing capacity calculations using the method of characteristics // Proc. 112-th Int. Conf. IACMAG, Turin. 2005. Рр. 441-450.
5. Hjiaj M., Lyamin A.V., Sloan S.W. Numerical limit analysis solutionsfor the bearing capacity factor Ny // International Journal of Solids and Structures 42. 2005. Pp. 1618-1704.
