Рис. 5. Расчетный участок трубопровода
h
= 1
м, то общее число элементов массива
W
(
i
)
должно быть не меньше
протяженности участка
НК
в метрах.
2. Формируется массив чисел
DW
(
i
)
для описания продольных смеще-
ний грунта
w
гр
в узловых точках конечно-разностной сетки.
3. Формируется массив
Q
(
i
)
для описания продольных реакций грунта
q
z
в узловых точках конечно-разностной сетки.
4. Задаются исходные данные: координаты точек
Н
,
А
,
В
,
K
;
w
гр
(
А
) =
=
DW
(
i
A
)
;
w
гр
(
B
) =
DW
(
i
B
)
;
h
=
HZ
;
q
пр
=
QPR
;
W
пр
=
WPR
;
δ
т
=
DT
;
D
;
E
. Все координаты должны быть больше единицы. Рекомен-
дуется принять
z
H
= 10
, а координаты остальных точек выбрать по их
расстояниям от точки
Н
(в метрах).
5. Определяются номера узлов характерных точек:
i
Н
= 10;
i
А
=
i
H
+
НА
;
i
В
=
i
H
+
НВ
;
i
К
=
i
H
+
НК
.
6. Определяется продольное смещение грунта в узлах конечно-разностной
сетки:
DW
(
i
) = 0
, если
i < i
А
;
DW
(
i
) =
DW
(
i
A
) +
DW
(
i
B
)
−
DW
(
i
A
)
i
B
−
i
A
(
i
−
i
A
)
, если
i
A
6
i
6
i
B
;
DW
(
i
) = 0
, если
i > i
B
.
Здесь смещения заданы линейно, но допустимы и другие аппроксимации.
7. Составляется цикл, состоящий из пунктов 7–9 с изменением индекса
i
от
i
H
до
i
K
с шагом 1.
8. Определяется реакция грунта
Q
(
i
) =
−
(
QPR/WPR
)[
W
(
i
)
−
−
DW
(
i
)]
; если полученное значение
Q
(
i
)
больше
+
Q
пр
, то принимает-
ся
Q
(
i
) = +
Q
пр
; если же
Q
(
i
)
<
−
Q
пр
, то принимается
Q
(
i
) =
−
Q
пр
.
9. Вычисляется комплексный коэффициент
R
=
HED
=
h
2
Eδ
т
.
10. Вычисляется
W
(
i
) = 0
,
5 [
W
(
i
−
1) +
W
(
i
+ 1) +
RQ
(
i
)]
.
11. Проверяется, есть ли изменение в массиве
W
(
i
)
в пределах
i
H
6
i
6
i
K
.
Если оно произошло по сравнению с предыдущим циклом, то повторяются
п. 6–9.
Сравнение последующего состояния с предыдущим проводится по сумме
абсолютных значений
W
(
i
)
на всем участке
SUMMA
=
i
к
X
i
н
|
W
(
i
)
|
.
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2011. № 4
115
