Сборник трудов конференции СПбГАСУ 2014 ч
.1.pdfРаздел 3. Аналитические и численные методы исследований оснований и фундаментов
Рис. 1. Схема напряженного (a) и деформированного (b) состояния части трубы
Данное выражение приводится к виду |
|
|||
d r |
|
r |
0 , |
(3) |
dr |
|
r |
|
|
Рассмотрим деформации элемента А1В1С1Д1(рис. 1b). В виду симметрии увеличение внутреннего давления приведет к радиальному перемещению всех точек цилиндра. Так, точки А1 и Д1 сместятся в радиальном направлении на величину (u) в положения А2 и Д2, а точки В1 и С1 – на величину u+du
в положения В2 и С2.
Тогда относительная радиальная деформация грани А1В1 будет равна
r |
А1В1 |
АВ |
|
ВВ1 АА1 |
(u du) u |
du |
, |
(4) |
||
|
|
|
||||||||
|
АВ |
АВ |
dr |
dr |
|
|
||||
Относительная окружная деформация грани А1Д1 равна |
|
|
||||||||
|
|
|
А2 Д2 А1 Д1 |
(r u) d r d u |
, |
(5) |
||||
|
|
|
||||||||
|
|
|
А1 Д1 |
r d |
r |
|
|
В соответствии с законом Гука для плоского напряженного состояния можно записать
r |
|
E |
( r ), |
(6) |
|
(1 |
2 ) |
||||
|
|
|
С учетом выше изложенного получим линейное дифференциальное уравнение 2-го порядка с переменными коэффициентами относительно радиальных перемещений, т. е.
511
Современные геотехнологии в строительстве и их научно-техническое сопровождение
d 2u |
|
1 du |
|
u |
0, |
|
dr2 |
r dr |
r2 |
||||
|
|
|
Общее решение данного уравнения сводится к виду u A r B / r,
(7)
(8)
где А и В – постоянные интегрирования, вычисляемые с учетом граничных условий.
Постоянные интегрирования А и В вычисляются при совместном решении следующих уравнений
при r rн, |
r |
|
E |
|
|
|
А(1 ) В |
1 |
|
н, |
|||||||
|
|
|
|
|
r2 |
|
|||||||||||
1 2 |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
н |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(9) |
при r rв, |
r |
|
E |
|
|
|
А(1 ) В |
1 |
|
в |
|||||||
|
|
|
|
r2 |
|
||||||||||||
1 2 |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
в |
|
|
|
Решение данных уравнений дает следующие значения |
|||||||||||||||||
|
|
|
|
1 |
|
2 |
|
|
2 |
н |
|
|
|||||
|
А |
|
rв |
в rн |
|
|
|||||||||||
|
|
E |
|
|
|
|
|||||||||||
|
|
|
|
r 2 |
r 2 |
|
|
|
, |
||||||||
|
|
|
|
|
|
|
|
|
|
н |
|
в |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(10) |
|
|
|
|
1 |
|
2 |
2 |
|
|
|
|
|
|
||||
|
В |
rв |
rн ( в н ) |
|
|||||||||||||
|
|
E |
|
|
r 2 |
r 2 |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
н |
|
в |
|
|
|
|
Для возможности использования решения данной задачи применительно к прессиометрическим испытаниям необходимо принять значение rн=∞ (бесконечное распространение грунта) и вычислить пределы
lim ( в rв2 н rн2 ) /(rн2 rв2 ) н,
rН
(11)
lim (rв2 rн2 )( в н) /(rн2 rв2 ) rв2 ( в н)
rН
Тогда радиальные перемещения можно вычислить по формуле
u |
rв |
(1 )( в н) н(1 2 ) , |
(12) |
|
E |
||||
|
|
|
Зная горизонтальные перемещения стенок камеры прессиометра (увеличение радиуса скважины) и давление внутри камеры σв и снаружи σн (давление грунта) модуль деформации грунта можно определить по следующей формуле
E |
(1 )( в н) н(1 |
2 ) |
, |
(13) |
u |
|
|||
|
|
|
|
512
Раздел 3. Аналитические и численные методы исследований оснований и фундаментов
В качестве наружного давления Рн принимали величину природного давления на глубине расположения прессиометра. При горизонтальной поверхности грунта вертикальные и горизонтальные напряжения от собственного веса грунта будут равны
zy zg h, |
(14) |
где γ – удельный вес грунта в пределах глубины расположения (h) зонда.
Всоответствии с решением Буссинеска о единой силе, приложенной
кповерхности грунта, горизонтальные напряжения равны
|
3P y2 z |
|
1 |
2 |
1 |
|
|
|
|
2R z y2 |
|
|
z |
|
||||||||
y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
||
R5 |
|
3 |
|
|
|
|
R z 2 |
R3 |
|
|
|
|||||||||||
|
2 |
|
|
R R z |
|
|
|
|
R3 |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(15) |
|
3P x2 z |
|
1 |
2 |
1 |
|
|
|
2R z y2 |
|
|
z |
|
|||||||||
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
R5 |
|
3 |
|
|
|
|
R z 2 |
R3 |
|
|
|
||||||||||
|
2 |
|
|
R R z |
|
|
R3 |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
R |
x2 y2 |
z2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Горизонтальные напряжения σх и σу от собственного веса грунта при его горизонтальной поверхности должны быть равными, что имеет место только при μ=0,5, т. е. в дальнейших вычислениях необходимо учитывать только первое слагаемое значение горизонтальных напряжений, например, σу можно получить двойным интегрированием выражения (15) в бесконечных пределах
|
|
3Pz |
( y )2 |
|
d |
|
|
Pz |
|
|
|
|
|||
y1 |
|
|
|
|
|
5 / 2 |
|
|
|
|
|
|
, |
(16) |
|
|
|
|
(x |
2 |
z |
2 |
) |
||||||||
|
|
2 ( y )2 x2 |
z2 |
|
|
|
|
|
|||||||
|
|
Pz |
d |
P h, |
|
|
|
|
|
|
|
|
|||
y2 |
|
|
|
|
|
|
|
|
|
|
|
(17) |
|||
(x )2 z2 |
|
|
|
|
|
|
|
Полученная формула учитывает коэффициент Пуассона грунта и глубину погружения прессиометра путем приложения внешнего давления, равного природному. Так как природное давление увеличивается пропорционально глубине погружения прессиометра, что ведет к снижению горизонтальных (радиальных) перемещений, а, следовательно, и увеличение модуля деформации грунта. Следовательно, в однородных грунтах модуль деформации грунта должен возрастать с увеличением глубины скважины (расположение прессиометра). В существующем ГОСТе эта зависимость обратная.
Литература
1. Методические рекомендации по определению деформационных и прочностных свойств глинистых пород методом прессиометрии. ВСЕГИНГЕО. М.1971.
513
Современные геотехнологии в строительстве и их научно-техническое сопровождение
2.ГОСТ 20276-99. Грунты. Методы полевого определения характеристик прочности и деформируемости. М.
3.ГОСТ 20276-85. Грунты. Методы полевого определения характеристик деформируемости
УДК 624.131.43
С.П. Дорошенко, А.А. Коршунов, А.Л. Невзоров
(САФУ, г. Архангельск)
ОЦЕНКА УСТОЙЧИВОСТИ ДАМБЫ С УЧЕТОМ ВЕРОЯТНОСТНОГО РАСПРЕДЕЛЕНИЯ ФИЗИКОМЕХАНИЧЕСКИХ СВОЙСТВ ГРУНТА
По результатам обработки данных лабораторных испытаний грунтов определяют нормативные и расчетные значения характеристик выделенных инженерно-геологических элементов. Причем для большинства характеристик за расчетное допускается принимать нормативное, то есть среднее значение. Это вполне приемлемо для небольших зданий и сооружений. Для гидротехнических сооружений, размеры которых достигают сотен метров, такой подход может привести к ошибкам при оценке надежности, так как кроме неоднородности основания, имеет место неоднородность грунтов карьера, сказывается и особенность укладки и уплотнения грунтов в теле сооружения. Нормы [5] при анализе устойчивости гидротехнических сооружений рекомендуют учитывать вероятностное распределение характеристик грунтов.
При расчете коэффициента запаса устойчивости дамб необходимо выявлять факторы, оказывающие наибольшее влияние на эту величину, а также оценивать вероятность появления ее критических значений.
Считается общепринятым, что большинство физико-механических характеристик грунтов имеет нормальный закон распределения.[2] Однако, в этом случае, особенно при средних значениях характеристик близких к нулю, например удельного сцепления или коэффициента фильтрации, допускается появление отрицательных значений расчетных параметров, не имеющие физического смысла.
А.Н. Колмогоров еще в середине XX века доказал [4], что размеры частиц при дроблении подчиняются логарифмически-нормальному закону распределения. Исходя из предположения о схожести искусственного измельчения и естественного выветривания горных пород, можно предположить, что логарифмически-нормальный закон распределения допускается применять и для анализа физико-механических свойств грунтов.
Сотрудниками кафедры инженерной геологии, оснований и фундаментов САФУ с 2004 года осуществлялся геотехконтроль за процессом возведения трехъярусной ограждающей дамбы хвостохранилища на одном из предприятий Архангельской области. Поперечный разрез дамбы показан на рис. 1.
514
Раздел 3. Аналитические и численные методы исследований оснований и фундаментов
Рис.1. Конструкция ограждающей дамбы
В данной статье представлены результаты оценки коэффициента запаса устойчивости дамбы с учетом вероятностного анализа свойств песчаноглинистых отложений (хвостов), намытых в экран верхового откоса (рис. 2). Обычно влиянию состава и свойств пляжных отложений на устойчивость дамб накопителей промышленных отходов должного внимания не уделяется.
Из отложений пляжа в ходе геотехконтроля было отобрано более 700 образцов, определены их состав, физические свойства, водопроницаемость и деформаци- онно-прочностные характеристики (табл. 1). Плотность частиц грунта в вероятностном анализе не рассматривалась из-за небольшого интервала варьирования.
Рис. 2. Намыв хвостовых отложений
Согласно [1], если физическая величина подчиняется логарифмическинормальному распределению, то ее логарифм является нормально распределенным, поэтому для всех физико-механических характеристик вычислялся натуральный логарифм, а далее находились основные показатели совокупностей: математическое ожидание и стандартное отклонение (табл. 1).
В отличие от хвостовых отложений 1 и 2 ярусов, для которых была известна совокупность значений прочностных характеристик c и , на 3 ярусе
использовалась совокупность значений предельного сопротивления сдвигу при различных нормальных напряжениях, а при обработке данных трехосных испытаний анализировался параметр N, который, согласно [3], зависит от
эффективных напряжений 1' и 3' . Результаты представлены в табл.2.
515
Современные геотехнологии в строительстве и их научно-техническое сопровождение
Таблица 1
Объемы лабораторных испытаний и основные показатели совокупностей физикомеханических свойств хвостовых отложений
Характеристика грунта |
Число |
Математическое |
Стандартное |
|
испытаний |
ожидание |
отклонение |
||
|
Пляж 1 |
яруса |
|
|
Плотность |
|
0,595 |
0,069 |
|
Плотность сухого грунта d |
211 |
0,372 |
0,041 |
|
Коэффициент пористости e |
|
-0,152 |
0,089 |
|
Коэффициент фильтрации k f |
19 |
-2,581 |
1,030 |
|
Модуль деформации E |
22 |
2,385 |
0,1681 |
|
Угол внутреннего трения |
22* |
3,348 |
0,065 |
|
Удельное сцепление c |
2,324 |
0,198 |
||
|
||||
|
Пляж 2 |
яруса |
|
|
Плотность |
224 |
0,587 |
0,073 |
|
Плотность сухого грунта d |
0,377 |
0,066 |
||
|
||||
Коэффициент пористости e |
121 |
-0,184 |
0,158 |
|
Коэффициент фильтрации k f |
40 |
-1,281 |
1,171 |
|
Модуль деформации E |
48 |
1,413 |
0,542 |
|
Угол внутреннего трения |
38* |
3,371 |
0,112 |
|
Удельное сцепление c |
2,090 |
0,864 |
||
|
||||
|
Пляж 3 |
яруса |
|
|
Плотность |
147 |
0,617 |
0,066 |
|
Плотность сухого грунта d |
0,389 |
0,060 |
||
|
||||
Коэффициент пористости e |
78 |
-0,159 |
0,136 |
|
Коэффициент фильтрации k f |
26 |
-2,105 |
1,002 |
|
Модуль деформации E |
52 |
0,951 |
0,470 |
|
Угол внутреннего трения |
30*+38** |
- |
- |
|
Удельное сцепление c |
- |
- |
||
|
Примечание: * – одноплоскостной срез;** – трехосное сжатие.
|
|
Таблица 2 |
Основные показатели для совокупностей прочностных свойств 3 яруса |
||
|
|
Стандартное откло- |
Характеристика грунта |
Математическое ожи- |
|
дание |
нение |
|
Касательные напряжения , 50 кПа |
-3,226 |
0,168 |
Касательные напряжения , 100 кПа |
-2,639 |
0,176 |
Касательные напряжения , 200 кПа |
-2,068 |
0,157 |
Параметр N |
1,317 |
0,165 |
516
Раздел 3. Аналитические и численные методы исследований оснований и фундаментов
Следующим шагом анализа, было определение интервала варьирования значений физико-механических свойств хвостовых отложений. В одномерном случае функция плотности вероятности для нормального распределения имеет вид [1]:
|
|
1 |
e |
(x )2 |
|
f (x) |
|
2 2 ; |
|||
|
2 |
||||
|
|
|
где x – случайная величина; – математическое ожидание; – стандарт-
ное отклонение.
Как пример, на рис. 3 показана функция f(x) для логарифма плотности хвостовых отложений пляжа пляжа первого яруса. От среднего значения откладывалось одно значение , что соответствует вероятности 68,26 %.
Вертикальными линиями показан взятый для анализа интервал значений. Следует отметить, что возможность появления значения больше верхнего и меньше нижне-
го предела составляет 15,87 %. |
|
|
Очевидно, что в зависимости от |
|
|
типа проектируемо го сооружения |
|
|
при проектировании |
необходимо |
|
расширять диапазон |
вероятности |
|
на одну-две величины . |
Рис. 3. Распределение логарифма плотности |
|
|
|
хвостовых отложений 1 яруса |
В результате обработки данных были получены интервалы варьирования значений всех характеристик хвостовых отложений (табл. 3–7).
Таблица 3
Интервалы варьирования значений удельного веса и удельного веса в водонасыщенном состоянии
Ярус дамбы |
|
γunsat, γsat, кН/м3 |
|
|
0 |
5 |
10 |
15 |
20 |
Первый
Второй
Третий
517
Современные геотехнологии в строительстве и их научно-техническое сопровождение
Таблица 4
Интервалы варьирования значений модуля деформации |
|
|
|
|
|||||||||
Ярус дамбы |
|
|
|
|
|
E, МПа |
|
|
|
|
|
|
|
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
Первый |
|
|
|
|
|
|
|
|
|
|
|
|
|
Второй
Третий
Таблица 5
Интервалы варьирования значений угла внутреннего трения
Ярус дамбы |
|
|
|
|
|
, град. |
|
|
|
|
|
||
0 |
3 |
6 |
9 |
12 |
15 |
18 |
21 |
24 |
27 |
30 |
33 |
36 |
39 |
Первый
Второй
Третий
Таблица 6
Интервалы варьирования значений удельного сцепления
Ярус дамбы |
|
|
|
|
с, кПа |
|
|
|
|
0 |
2 |
4 |
6 |
8 |
10 |
12 |
14 |
16 |
18 |
Первый
Второй
Третий
Таблица 7
Интервалы варьирования значений коэффициента фильтрации
Ярус дамбы |
|
|
|
kf, м/сут |
|
|
|
|
|
0 |
0,1 |
0,2 |
0,3 |
0,4 |
0,5 |
0,6 |
0,7 |
0,8 |
0,9 |
Первый
Второй
Третий
С учетом интервалов варьирования характеристик выполнялось моделирование дамбы в программно-вычислительном комплексе "Plaxis2D 2013". Использовались модели Мора-Кулона и упрочняющегося грунта (Hardening soil).
Согласно [6], в первой модели кроме физических свойств и коэффициента фильтрации используются следующие параметры: угол внутреннего трения , удельное сцепление c, модуль деформации E, коэффициент Пуас-
сона , а также угол дилатансии , который принимался по зависимости,
518
Раздел 3. Аналитические и численные методы исследований оснований и фундаментов
приведенной в руководстве 30 , при отрицательных значениях задавали 0 . Дополнительные характеристики грунтов для второй модели бы-
ли вычислены согласно рекомендациям [6].
Для расчета водопроницаемости при установившемся и неустановившемся режимах фильтрации использовалась модель Ван-Генухтена [6].
Алгоритм вычислений коэффициента запаса устойчивости основывался на процедуре пошагового изменения характеристик грунтов до момента отказа сооружения. Моделировалось поярусное возведение дамбы и поэтапное заполнение хвостохранилища.
На первом этапе моделирования задавались нормативные (средние) значения всех физико-механических характеристик и определялся коэффициент запаса устойчивости. Далее, по результатам вероятностного анализа, изменялась одна из характеристик – задавались ее минимальное и максимальное значения и вычислялись соответствующие им коэффициенты запаса устойчивости. Остальные свойства при этом оставались неизменными.
Результаты моделирования с использованием модели Мора-Кулона показали, что значительное влияние на устойчивость сооружения оказывает коэффициент фильтрации хвостовых отложений (табл. 8). Влияние остальных физико-механических свойств пляжных хвостовых отложений было не столь велико. Следует заметить, что согласно [6] в расчетах задают, как правило, только среднее значение коэффициента фильтрации.
Таблица 8
Влияние коэффициента фильтрации на устойчивость сооружения
Ярус |
|
Коэффициент запаса устойчивости |
||
дамбы |
при значении коэффициента фильтрации |
|||
|
минимальном |
|
среднем |
максимальном |
Первый |
1,418 |
|
1,455 |
1,474 |
|
|
|
|
|
Второй |
1,385 |
|
1,455 |
1,497 |
|
|
|
|
|
Третий |
1,427 |
|
1,455 |
1,490 |
|
|
|
|
|
При наихудшем и наилучшем распределении всех физикомеханических свойств отложений пляжа коэффициент запаса устойчивости варьировался в следующих пределах 1,310<kst<1,546.
По тому же алгоритму проводилось моделирование с использованием модели упрочняющегося грунта (Hardening soil). Вычисления показали, что при варьировании физико-механических свойств хвостовых отложений коэффициент запаса устойчивости практически не отличается от результатов, полученных с использованием первой модели. Минимальное значение коэффициента запаса устойчивости равнялось 1,316, максимальное 1,547.
519
Современные геотехнологии в строительстве и их научно-техническое сопровождение
Следует заметить, что при численном моделировании учитывалась возможность как общей потери устойчивости трехъярусной дамбы, так и локального обрушения низового откоса на одном из ярусов (рис. 4 и 5).
Рис. 4. Призма обрушения при локальной потере устойчивости
Рис. 5. Призма обрушения при общей потере устойчивости
Анализируя результаты работы, можно сказать, что при расчете гидротехнических сооружений, необходимо учитывать вероятностное распределение физико-механических свойств грунтов. Такой подход дает более полное представление о надежности сооружения. В работе анализировались только хвостовые отложения намытые в экран ограждающей дамбы, однако аналогичный подход должен использоваться для основания и тела сооружения. Кроме того тело дамбы следует зонировать в пространстве, учитывая изменение физико-механических характеристик грунтов.
Литература
1.Айвазян С. А. и др. Прикладная статистика: Основы моделирования и первичная обработка данных. Справочное изд. С. А. Айвазян, И. С. Енюков, Л. Д. Мешалкин. — М.: Финансы и статистика, 1983. — 471 с.
2.Гольдштейн М. Н. Механические свойства грунтов. – Изд. 2-е, переработ. – Москва: Стройиздат, 1971. – 368 с.
3.ГОСТ 12248-2010 Грунты. Методы лабораторного определения характеристик прочности и деформируемости. МНТКС, 2011.
4.Колмогоров А.Н. О логарифмически-нормальном законе распределения размеров частиц при дроблении. Докл. СССР, 1941, т. 31, №2, с. 99.
5.СП 23.13330.2011 Основания гидротехнических сооружений. Актуализированная редакция СНиП 2.02.02-85 М.: ФГУП ЦПП, 2011.
6.User Guide Plaxis 2D 2012, Plaxisbv, Netherlands, 2012.
520