Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Расчет водонасыщенных грунтовых оснований. Монография

.pdf
Скачиваний:
0
Добавлен:
08.09.2026
Размер:
2 Мб
Скачать
     
 
 
n
i
iiii
TthTtaaтtт
0
1
10
000
1
оnо
T,a,a
01
п
а
1i
a
i
a
i
Tt
 
 
 
n
i
Tр
ii
i
е
р
aaтрт
0
1
*
1
10
р
р
р0
 
0
*
mрm
Рис. 1.7. Линейная комбинация звеньев ломаной линии
Важно отметить, что дефекту I отвечает стыковка звеньев ломаной ли­нии по ординатам, поэтому у каждой функции базиса имеется только по од­ному параметру [37].
В общем случае для звеньев ломаной линии получается следующее аналитическое выражение:
, (1.16)
.
Условие
обеспечивает горизонтальность конечного (асимпто-
тического) участка сплайна.
Специальная запись сплайна, в которой используется функция Хевисайда с разностным аргументом, приводит к тому, что параметр последующего звена
находится независимо от параметра предыдущего звена
Аргументы
сомножителей перед функциями Хевисайда совпа-
[37].
дают, поэтому можно применить теорему запаздывания и записать изобра­жение сплайна (1.16) по Лапласу-Карсону [37]:
, (1.17)
где
– точки коллокаций.
Первая точка коллокаций назначается как предельная точка
на
полуоси
. В результате все экспоненты выражения (1.17) зануля-
ются, и начальный параметр определяется выражением [37]
. (1.18)
31
Согласно теореме о пределах предельное значение оригинала
 
tm
t
рm
*
0р
рmtm
*
0
0* рmtm
 
00 tmm
0р
0
0
р
е
i
Tр
о
р
0
п
р
при
совпадает с предельным значением изображения
и наоборот (рис.1.8) [37]:
,
.
Рис. 1.8. Взаимосвязь функций оригинала и изображения
Тогда выражение (1.18) запишется в виде
при
, следова-
тельно, начальное значение искомого оригинала по методу ломаных все­гда находится точно.
В соответствии с записью сплайна в изображениях (1.17) при возникает неопределенность:
,
которую можно раскрыть путем разложения экспоненты в ряд или с помо­щью применения теоремы о пределах [37].
В связи с этим, можно сделать вывод, что точки
и
все-
гда имеют конечное значение и называются предельными, все остальные точки называются внутренними.
В главе 2 данной работы (п.2.1) будет предложен способ назначения точек коллокаций из решения трансцендентного уравнения.
Для того чтобы приближенно найти оригинал известной в изображениях ломаной линии по методу ломаных вводится система точек коллокаций и запи­сывается на этой системе условие совпадения известной функции с изображе­нием ломаной линии. В результате составляется СЛАУ, из решения которой определяются углы наклона (в радианах) звеньев ломаной линии [37].
32
Решение вязкоупругой задачи по методу ломаных также (п.1.3) на
 
tW
s
     
 
 
n
i
iiii
ТthТtaaFtF
0
1
10
 
 
 
 
n
i
ТpТp
i
ii
ee
p
aFpFtF
1
*
*
1
1
10
i
a
01
nî
aa
0
o
Т
 
tW
s
 
tF
00 tFtW
s
tFtW
s
pWtW
ss
*
0
pFtF
*
0
0* pWtW
ss
0* pFtF
 
0F
pWpF
s
*
*
 
 
0
1
10
1
*
1
FeeaFpF
n
i
ТТ
i
ii
 
 

 
   
 
 
 
.1
0
0
00
10
0
2
*
2
*
*
*
2
*
2
Ha
Hpa
s
s
e
a
П
f
e
pa
pП
pfpW
 
  
     
   
 
Ha
e
a
П
fF
0
2
2
1
0
0
000
примере определения перемещений частиц скелета грунта. Введем специа­лизированный сплайн, аппроксимирующий оригинал
который в изображениях по Лапласу-Карсону имеет вид:
где
– искомые параметры;
Очевидно, что заданная функция
;
.
и искомая аппроксимация
должны иметь общие начальные и конечны значения [37]:
;
.
Согласно теореме о пределах начальные и конечные значения задан-
ной функции и искомой аппроксимации совпадают [37]:
(1.13):
,
,
Начальный параметр
Тогда,
;
;
;
.
определим из условия совпадения
;
;
.
33
Из условия совпадения искомой аппроксимации
pF
*
pW
s
*
nj
j
j
p
 1
j
s
j
pWpF
*
*
 
j
s
п
i
ТpТp
j
i
pWee
p
aF
ijij
*
1
1
1
10
 
 
 
j
j
s
п
i
ТpТp
i
p
F
pW
eea
ijij
 
 
1
0
*
1
1
0p
t
 
1
0
0
*
1
1
F
pW
ТТa
s
п
i
iii
tWpW
ss
0
*
 
 
 
 
1
0
0
................................................................
1
0
*
1
1
*
1
1
F
pW
ТТa
p
F
pW
eea
s
п
i
iii
j
j
s
п
i
ТpТp
i
ijij
i
a
 
tF
pW
s
*
i
Т
j
p
p
*
pm
*
p
*
pm
*
с заданной функ-
цией
точки
на системе точек коллокаций
составим СЛАУ [37]:
,
,
.
Последнее уравнение СЛАУ записывается специальным образом для
в изображениях или
в оригинале [37]:
,
.
В результате получаем:
ются искомые параметры
описывающее искомый оригинал для функции ного решения при одном и том же порядке СЛАУ сильно зависит от выбора
временных точек
СЛАУ. Известная в изображениях функция цией специального вида ров, подбирая которые можно аппроксимировать функцию того, функция
параметров ее оригинал можно представить в явном виде [37].
.
Решая систему n линейных уравнений с n неизвестными, определя-
, и получается выражение
, приближенно
. Точность получен-
и точек коллокаций
.
Таким образом, метод ломаных сводится к составлению и решению
аппроксимируется функ-
, зависящая от переменной p и от парамет-
должна иметь такой вид, чтобы при любом количестве
34
. Кроме
Основная сложность метода ломаных состоит именно в назначении
i
Т
j
p
2
2500см
5.2:1
5.3:1
s
W
l
W
s
l
Ha
za
s
о
i
s
ee
Eа
zzW
i
2
2
2
)(
Ha
l
o
l
e
E
h
zW
2
1)0(
точек коллокаций, так как чем точнее будут они найдены, тем точнее будет значение искомой функции. Более подробно вопрос о назначении времен­ных точек
и точек коллокаций
будет рассматриваться в п.2.1 [37].
1.5. Определение механических характеристик водонасыщенных грун-
тов с учетом давлений воды в порах
Любые строительные объекты строятся на грунтовом основании или располагаются в толще грунта. Под влиянием приложенных нагрузок грунты основания деформируются, и сооружения испытывают вертикальные и гори­зонтальные перемещения. Эти перемещения не происходят сразу же и пол­ностью по приложению нагрузки, а длятся некоторое время [37].
Таким образом, возникает необходимость прогнозирования величины пе­ремещений, а также продолжительности их протекания во времени, для чего необходимо определять соответствующие механические характеристики [37].
Согласно [28-30] механические характеристики определяют на осно­вании исследований с использованием нескольких приборов: одноосного сжатия, одноплоскостного среза, компрессионного сжатия и трехосного осесимметричного сжатия (стабилометр). Испытаниям подвергают образцы с ненарушенной структурой (в виде паст заданной плотности и влажности) с отношением диаметра d к высоте h 1:2 или 1:1,5. Недостатком таких схем приборов является малая высота и малый диаметр испытуемых образцов, вследствие чего, двухфазный образец переходит в однофазный, так как по­ровая вода не зажимается частицами грунта и полностью сбрасывает с себя нагрузку. Следовательно, испытания образцов стандартных размеров [1, 28­30] не могут характеризовать систему в целом и отображают ее в плане то­чечного определения.
Для более достоверного представления механических свойств двухфаз­ной системы предлагается [4, 58] испытывать образцы с площадью попереч-
ного сечения до
и с отношением диаметра d к высоте h
, удаленные от дневной поверхности и расположенные в зоне, где боль-
шую часть нагрузки воспринимает на себя поровая вода.
Двухфазная система моделируется упругим вариантом кинематиче­ской модели [53], который отвечает стабилизированному состоянию грунта.
На основании решения дифференциального уравнения (1.9) с задан­ными граничными условиями (1.10) определяются перемещения частиц
твердой
и жидкой фазы
, напряжения в скелете грунта
и в жидкой
фазе [37]:
;
35
или
;
i
za
oi
s
ezz
2)(
;
i
za
oi
l
ezz
2
1)(
1-
ма
2
i
zz
i
l
ex
za
oi
l
zzezz
i
2
1)(
 
 
o
i
l
ex
i
zz
z
a
)(
1ln
1
2
)(1МПа
Е
s
i
zz
)()(
2
2
2
i
s
ex
Ha
za
s
о
i
s
zzWee
aЕ
zzW
i
s
Е
1
Ha
o
i
s
ex
s
e
azzW
E
2
1
1
2
l
Е
1
0z
01)0(
2
zWe
E
h
zW
l
ex
Ha
l
o
l
Ha
o
l
ex
l
eh
zW
E
2
1
01
1-
ма
2
)(1МПа
Е
s
.
Универсальный параметр кинематической модели
определя-
ется [4] из условия совпадения теоретического значения порового давления в сечении
с экспериментальным значением той же величины при ста-
билизированном состоянии образца [37]:
.
Индекс «ex» отвечает слову «experiment - эксперимент».
После преобразований получается:
.
Параметр
определяется из условия совпадения теоретиче-
ского значения перемещения поперечного сечения образца (скелета) в сече­нии
с экспериментальным значением той же величины [37]:
.
В результате преобразований записывается выражение для определе-
ния механической характеристики
[37]:
.
Параметр
чения перемещения жидкой фазы в сечении
определяется из условия совпадения теоретического зна-
с экспериментальным зна-
чением той же величины при стабилизированном состоянии образца [37]:
.
Окончательная формула для этого параметра имеет вид [37]:
Исследования [4, 58] показали, что значения параметров не зависят от ме­ста расположения датчиков в образце и применяемой методики их определения.
Изменение параметров
и
во времени будет описано в п.2.3.
36
2. ПРИМЕНЕНИЕ МЕТОДА ЛОМАНЫХ
р
*
 
t
р
*
 
tm
рm
*
р
*
рm
*
ii
рртрр
**
ni ,...,0
ni
i
i
Т
 1
ni
i
i
р
 1
ni
i
i
Т
 1
Т
о
Т
0
о
Т
п
Т
пп
ТТ 1,01
12
1,0ппТТ
23
1,0ппТТ
ni
i
i
Т
 1
ni
i
i
р
 1
К ОБРАБОТКЕ ЭКСПЕРИМЕНТА
2.1. Способы назначения точек коллокаций
Первоначально для сокращения записей вводятся следующие обозна­чения [37]:
– известное в изображениях решение вязкоупругой задачи для
фиксированной точки тела;
– искомый оригинал решения
и
– оригинал и изображение стандартной функции метода
;
ломаных, определяемые соответственно по формулам (1.16) и (1.17).
В п.1.4 отмечалось, что для определения оригинала по методу ломаных известную в изображениях функцию
циального вида
(1.17) с использованием метода коллокаций [37]:
аппроксимируют функцией спе-
,
.
Процесс нахождения приближенного оригинала сводится к экспонен­циальной аппроксимации функции (1.17), недостаток которой заключается в том, что после некоторого (оптимального) количества членов ряда, при дальнейшем увеличении числа слагаемых точность аппроксимации будет ухудшаться, и функция (1.16) начнет осциллировать между точками колло­каций. Кроме того, точность аппроксимации при одном и том же порядке
СЛАУ существенно зависит от выбора узлов
и точек коллокаций
, что будет показано в п.2.2, 2.4 [37].
Для составления СЛАУ по методу ломаных первоначально необхо­димо задаться системой временных точек
, образующих сильно воз-
растающую последовательность. Это связано с тем, что функции, характе­ризующие вязкоупругие свойства меняются быстро на начальном участке, и по мере возрастания времени скорость их изменения убывает и стремится к нулю при
Начальная точка ; конечная точка
[37].
соответствует моменту загружения, то есть
берется равной значению времени t, при котором
можно считать, что все механические характеристики уже достигли своих асимптотических значений, последующие точки должны отличаться друг от друга на порядок [37]:
Рассмотрим алгоритм назначения временных точек локаций
,
метода ломаных, предложенный Шалабодовым В. И [90].
,
, и т.д.
и точек кол-
37
Известная зависимость между безразмерными аргументами t и р:
 
р
dtteрt
1
0
*
i
i
T
р
1
1,...,2,1 ni
0
1
о
о
Т
р
01
п
п
T
р
о
р
0
п
р
о
р
0
п
р
iiTT ;
 
ii
TthTthtf
ii
TрTр
eeрf
*
0*
i
Tр
i
р
eТрf
2
lnln
ii
TT
р
.
В монографии Работнова Ю.Н. [71] отмечается, что если бы аргу­менты t и р были размерными, то их размерности были бы взаимно обрат-
ными. Формально перенесем (обоснование приведено ниже) это соображе­ние на численные соотношения фиксированных значений [37]:
,
, (2.1)
что приводит к назначению точек коллокаций по узлам ломаной линии.
Первая точка назначается для начального момента времени:
,
конечная – для времени, при котором механические характеристики дости­гают своих асимптотических значений:
.
Точки
и
называются предельными, они являются условием
того, что начальное и конечное значения искомой аппроксимации и задан­ной функции будут совпадать. В связи с этим, можно сделать вывод, что рассматриваемые способы назначения точек коллокаций разработаны для определения внутренних точек, а две предельные точки
и
всегда имеют конечное значение [37].
Для обоснования взаимообратной связи узлов сплайна и точек колло­каций необходимо заключить узел сплайна в малый отрезок
и
на этом отрезке ввести ступенчатую функцию
которой по теореме запаздывания отвечает изображение в виде разности экспонент
Необходимое условие экстремума
где определяется точка экстремума
На основании известного разложения логарифма
,
.
,
.
38
xx 1ln
xx 1ln
имеем
i
ii
T
TT
lnln
i
ii
T
TT
lnln
i
i
i
i
i
i
T
T
T
T
T
р
1
2
lnln
lim
0
i
t
i
р
1i
T
i
T
ii
TT ;
1
tT 0
 
ii
TthTthtf
1
ii
TрTр
eeрf
1
*
что позволяет перейти к пределу
,
;
,
.
Исследование связи
и
узлов сплайна. Для этого на отрезке
и
, основанной на двух соседних значениях
полуоси
необ-
ходимо ввести ступенчатую функцию (рис.2.1):
,
и на основании теоремы запаздывания можно переписать ее в изображениях по Лапласу-Карсону (рис.2.2):
.
Рис. 2.1. График ступенчатой функции
Рис. 2.2. Изображение ступенчатой функции
39
Условие экстремума путем приравнивания производной от изображе-
0
1
/
*
ii
Tр
i
Tр
ii
eTeTрf
ii
Tр
i
Tр
ii
eTeT
1
iiiii
TрTTрT
lnln
1
i
рр
1
1
lnln
ii
ii
i
TT
TT
р
1...,3,2 ni
ii
TT ;
1
ni
i
i
Т
 1
11
ijijijij
TpTpTpTp
eeeem
10 m
ijij
TpTp
ii
eea
1
11
ijij
TpTp
ii
eea
1т
ния к нулю:
Отсюда
Логарифмируя, получается:
.
.
.
Вводится обозначение для точки экстремума
и выражение для
ее определения:
,
. (2.2)
Основной особенностью способа назначений точек коллокаций по формуле (2.2) является то, что он отражает информацию оригинала на всем
отрезке
, а не только в окрестности точки
, как это было в пер-
вом случае. Однако рассмотренный прием не учитывает влияние соседних участков [37].
Предложенный Шалабодовым В. И алгоритм назначения точек колло­каций на основании формул (2.1) и (2.2) позволяет описывать монотонные функции. В случае немонотонных функций рассмотренный алгоритм можно применять только для СЛАУ малого порядка (до пяти), что будет показано в п.2.3 [37].
Автором предложен новый способ назначения точек, который позво­ляет описать немонотонные функции ломаной линией по методу ломаных с количеством звеньев больше пяти.
Осцилляция функций появляется в результате неточности аппрокси­мации, что вызвано накоплением погрешностей в процессе расчета. Чтобы повысить точность приближенного решения необходимо улучшить обу­словленность матрицы СЛАУ, для этого точки коллокаций предлагается находить из решения трансцендентного уравнения [37]:
использование которого приводит к тому, что элементы матрицы СЛАУ, стоящие на главной диагонали
получаются больше следующих за ними в строке элементов
в m
раз.
,
, (2.3)
40
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]