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

Сборник трудов конференции СПбГАСУ 2014 ч

.1.pdf
Скачиваний:
252
Добавлен:
29.03.2015
Размер:
41.4 Mб
Скачать

Раздел 3. Аналитические и численные методы исследований оснований и фундаментов

L1

 

 

dL

 

 

 

 

L0

 

 

 

 

 

;

(13)

 

2 v T

L P

L dL

0 1

 

 

 

 

2E'

 

o

 

1

 

 

 

 

 

 

 

0

 

 

 

 

Некоторые задачи, решаемые в квадратурах.

а) Р1=0; Рn=Pn(z,y); уравнение равновесия у”=P(z,y)(1+y’2)3/2. При – Р(z,y)/T=2a(y+bz+c)=2au

 

 

 

b

 

 

 

 

 

 

1

 

 

 

 

au2 D

 

 

 

 

z(u)

 

 

 

 

 

 

u

 

 

 

 

 

1

 

du D2

;

(14)

b2 1

b2 1

b2 1 au2 D1 2

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

au2 D

 

 

 

 

y(u)

 

 

 

 

 

u

 

 

 

 

 

1

du bD2

C1

;

b2

 

1

b2 1

 

b2 1 au2 D1 2

 

 

 

 

 

 

 

 

 

где D1, D2 – постоянные интегрирования, интегралы эллиптического типа.

б) Руу(z), Рzz(z) вводим N(z)=2Py/ω, M(z)=2Pz(ω):

 

z

 

 

 

 

 

 

 

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

(15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

z

 

 

 

 

z

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

exp

Ndz exp M exp

Ndz dz

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

я

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Руdz T0 sin 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

в) Оболочка с собственной массой G/p, кг/м2 при P1=0, Р’n= Р’n(z),

G=G(z):

 

 

 

 

 

 

 

 

T0 sin 0 fn dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

у

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

T 2 cos2

 

2T

f

 

sin

 

 

f

 

f

2 f 2

 

 

 

 

 

0

 

 

0

n

0

g

 

 

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

n

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fn

Pn 'dz ;

f g Gdz .

 

 

 

 

 

 

 

 

 

 

(16)

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В частном случае G=const, P’n=const=P – решение Б.И. Друзя:

z To

1 P Gsin 0

 

;

 

 

y

D

 

 

Pcos

 

 

GA

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

G

 

P Gsin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P Gsin

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Gcos

 

 

 

 

 

 

 

 

 

 

 

 

T (P2

 

G sin

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

L 2D

 

 

 

 

 

 

 

PA

 

 

 

 

,

где D

 

0

 

 

 

 

 

 

 

 

 

 

0

 

;

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

G

 

 

 

 

 

 

 

P Gsin

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

Ptg

 

G

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

arctg

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

(17)

 

 

 

 

 

P2 G2

 

 

P2

G2

 

 

 

 

 

 

 

 

 

Т=То – Gz [при любом Р’n(z,y), составляющие массы в формулах пунктов (в, д) в Р’n искусственно не введены].

г) Р1=0, Рn=Pn(L):

571

Современные геотехнологии в строительстве и их научно-техническое сопровождение

 

0

 

1

L P dL ;

z L cos L dL ;

y L sin L dL ;

(18)

 

 

Т

 

n

 

 

 

 

 

0

 

0

0

 

д) Р’n= Р’n(α), G=const, P1=0:

 

 

 

 

 

T

 

o cos ad

 

z

o

1

exp

G

 

 

;

 

 

G

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

cos ad

 

 

 

 

 

 

y To

 

G

sin ad

; Pn

G sin ; T T0

Gz . (19)

exp

 

 

 

 

 

 

 

 

 

 

 

 

е) P1=0, Рn= Рn(α) [из случае (д) не следует]:

 

 

 

 

 

 

 

 

 

 

 

0

cos ad

 

 

 

 

0

sin ad

 

 

 

 

 

z T

 

P

 

;

 

y T

 

P

 

 

.

 

(20)

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

Если представить P P P

u P

u2 ... P

 

u m ; u=sin α, то

 

 

 

n

 

 

0

 

(1)

 

(2)

 

 

 

 

 

 

 

(m)

 

 

 

 

u

 

 

 

 

 

 

u0

 

 

 

du

 

 

 

 

 

 

 

 

 

 

z T 0 du ;

y T

 

 

 

 

 

 

.

 

 

 

(21)

 

P

1 u2

 

 

 

 

 

 

 

 

u Pn

 

 

 

u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

ж) Р1= Р1(α), Рn= Рn(z); при u=sin α решение следует из системы

d 2u

 

 

P

 

 

 

du

 

2

dP

 

du

 

 

1

 

 

 

 

dz2

 

 

1

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

0

;

 

 

 

 

 

 

 

 

dz

 

P

 

 

P

1 u

2

 

dz

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

n

 

 

 

dz

 

 

udu

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

C .

 

 

 

 

 

 

 

(22)

 

 

 

du

 

1 u2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В частном случае действия кулоновской силы трения Р1=fPn=tgφPn решение первого из этих уравнений: cos e fa sin( ) B A Pn z dz ,

где В и А – постоянные интегрирования; А, Рn, du/dz взяты по абсолютной величине.

з) Протягивание цилиндрической оболочки поперек твердого цилиндра с радиусом кривизны R=R(β), где β=π/2-α при действии сил трения Р1=fPn; T=T0e-fβ, Pn=T0e-fβ/R(β)зависимости Эйлера. Характерно, что падение усилия в оболочке от сил трения при протягивании не зависит от формы цилиндра R(β); последняя влияет только на контактные давления Pn.

Методы расчета цилиндрической оболочки при произвольно распределенной вдоль сечения нагрузке. Весьма распространен графоаналитический метод У. Томсона (Кельвина). Он удобен при действии одних

нормальных нагрузок

 

Рn Pn z, y, , L,... .

(23)

Используется уравнение равновесия, из которого

 

R T P ,

(24)

n

 

Идея метода – аппроксимация кривой сечения дугами окружностей с общей касательной в месте стыковки. Задают натяжение Т=const и угол α0

572

Раздел 3. Аналитические и численные методы исследований оснований и фундаментов

(рис. 1, д). По формуле (13) при α=α0, z=y=L=0 вычисляют для начала

координат Pn(0); по (14) R(0) и этим радиусом проводят дугу О1. В качестве шага назначают Δα или z, L и т. п. Для точки 1 снова вычисляют Pn(1) и Pn(2) по формулам (23) и (24); радиусом P(1) проводят дугу 1-2 и т. д. Конец

процесса определяется граничными условиями.

Аналитический метод, имеющий ту же основу, применим при действии на оболочку нормальных и касательных нагрузок. Когда в качестве переменной принят угол α, алгоритм выглядит так.

Заданы в точке О Т(0)0, α(0)0, шаг Δα; z(0)=y(0)=L(0)=0. Вычисляют по

известным зависимостям Рn(0)

Pn (0) , z(0) , y(0) , L(0) ,... ;

Р1(0)

P1 (0) , z(0) , y(0) , L(0) ,... .

 

Шаг 1.

R(0) T0 Pn(0) ;

(1)

0

;

z(1)

2sin(

2) cos 0 2 R(0) ;

 

z(1)

z(1) ; y(1) 2sin 2 sin 0 2 R(0) ; y(1) y1;

L(0) R(0) ;

 

 

 

L(1)

L(1) ; Р1(1)

P1 (1) , z(1) , y(1) , L(1) ,... ;

Рn(1) Pn (1) , z(1) , y(1) , L(1) ,... ;

 

 

 

 

 

 

 

T(1)

T0 P1(0) P1(1)

R(1) 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Шаг 2. R(i 1)

T(i 1)

Pn(i 1) ;

(i)

(i 1)

;

 

 

 

 

 

 

 

 

 

 

 

z(i) 2sin(

2) cos (i 1)

 

2 R(i 1) ;

z(i)

z(i 1) z(i) ;

 

 

 

 

 

 

 

 

 

 

 

y(i) 2sin

2 sin (i 1)

2 R(i 1) ;

y(i) y(i 1)

yi ;

L(i) R(i 1) ;

 

 

 

 

L(i)

L(i 1) L(i) ; Р1(i) P1 (i) , z(i) , y(i) , L(i) ,... ;

Рn(i)

Pn (i) , z(i) , y(i) , L(i) ,... ;

 

 

 

T(i) T(i 1) P1(i 1) P1(i) R(i 1) 2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подъем плоского груза длинной цилиндрической оболочкой

 

 

 

 

 

 

1 p

qp 2 2 p q 1

p

qp

 

 

b

 

 

 

b 2

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

1

 

 

; t

;

(25)

 

 

 

 

 

 

 

 

 

p

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

E1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q qPR R L ; P

; T tPR ; B bR ;

 

1

E / E

2

,

 

 

 

 

 

 

 

 

 

 

 

0 0 k

 

 

R0

 

 

 

0

0

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Q – подъемная сила; Р – давление; R0 – исходный радиус сечения; Lk – расчетная длина оболочки; Т – кольцевые натяжения; В – высота оболочки под грузом; ν,ν1, Е1, Е2 – коэффициенты Пуассона и модули упругости ортотропного материала в кольцевом и продольном направлениях. Для оболочек с тюбичными оконечностями длиною L Lk L 2,8B. Для нерастяжи-

мой оболочки B 2R0 Q PLk ; Т РВ2.

На основании проведенного анализа автором будут разработаны алгоритмы инженерного метода расчета оболочек, подверженных внешним статистическим и динамическим нагрузкам, которые в дальнейшем будут учтены в результате, экспериментальных исследований.

573

Современные геотехнологии в строительстве и их научно-техническое сопровождение

Литература

1.Применение грунтонаполняемых оболочечных конструкций для укрепления оснований в сложных геологических условиях / Буняев М.С./ Дефекты зданий и сооружений. Усиление строительных конструкций // Материалы XVII научно-методической конференции ВИТУ (21 марта 2013 года) / ВИ(ИТ) ВАМТО (ВИТУ). – Спб., 2013. – 204с.

2.Кашарина Т. П. Совершенствование конструкций, методов научного обоснования, проектирования и технологии возведения облегченных гидротехнических сооруже-

ний / Т. П. Кашарина. – М.,2000. – 40с.

3.Кашарина Т. П. Мягкие гидросооружения на малых реках и каналах / Т. П. Кашарина. – М.: Мелиоратизация и водное хозяйство,1997. – 56с.

4.Магула В. Э. Судовые эластичные конструкции / В. Э. Магула – Л. : Судо-

строение, 1978. – 263 с.

5.Шихирин В.Н., Ионова В.Ф., Шальнев О.В., Котляренко В.И. Эластичные механизмы и конструкции. Монография. – Иркутск: Издательство ИрГТУ, 2006. – 286с.

6.Дыба В.П. Определение формы абсолютно гибкой подпорной стенки и напряженного состояния грунтовой засыпки / В.П. Дыба, В.В. Лифанов // Изв. Сев.-Кавк. науч. центра высш. шк. техн. Науки. – 1983. – №3. – с. 76-78.

7.Х. Гибкие бункера. – М.: Стройиздат, 1980. – 168 с.

УДК 539.3

Д.А. Бекмирзаев, Н.А. Нишонов, Т. Юлдашев

(ИМиСС АН РУз, г. Ташкент)

МЕТОД КОНЕЧНЫХ РАЗНОСТЕЙ В СЕЙСМОДИНАМИКЕ ТРУБОПРОВОДОВ, ВЗАИМОДЕЙСТВУЮШИХ С ГРУНТОМ

1. Как показывают последствия прошедших землетрясений в мире подземные трубопроводы в большинстве случаев разрушаются в резких поворотах, присоединения трубопроводов к массивным узлам, в прохождении трубопроводов через разломы и др. Больше всего влияют свойства грунтовых условий. Чем мягче грунт, тем больше повреждений. Если считать, что около 1/3 расходы на строительство приходить на системы жизнеобеспечения (данные ООН), то естественно вопросы исследования напряженнодеформированного состояния подземных трубопроводов при сейсмических нагружениях являются актуальной.

В стенах института Механики и сейсмостойкости сооружений им. М.Т. Уразбаева АН РУз проведены экспериментальные и теоретические исследования сейсмодинамики подземных трубопроводов. В частности установлено наличие относительного движения трубы в грунтах при сейсмических нагружениях. Изучены экспериментальные значения коэффициент равномерного сдвига трубопроводов относительно грунта [1, 2, 3].

Составлены дифференциальные уравнения движения сложных систем подземных трубопроводов, сформулированы граничные и начальные условия

574

Раздел 3. Аналитические и численные методы исследований оснований и фундаментов

иусловия стыковок трубопроводов на сложных узлах, как в случае жестких

иподатливых стыковках [3, 4].

Ранее также были получены отдельные численные результаты напряженно – деформированного состояния подземных трубопроводов. Ныне продолжается изучения влияния водонасышенности грунтовых условий, а также сформулированы отдельные нелинейные задачи сейсмодинамики подземных трубопроводов.

Вданной статье сделан некоторый подход использованию метода конечных разностей к исследованию напряженно деформированного состояния подземных трубопроводов при сейсмических нагружениях.

Даны алгоритмы и программа расчета, приведены отдельные численные результаты, построены графики изменения напряжений и перемещений сечений трубопроводов в зависимости от граничных условий и влияния грунтовых условий.

Здесь главным образом обращается внимания к вопросам продольного колебания подземных трубопроводов при сейсмических воздействиях.

Вконечном счете проводимые исследования будут полезным для дополнения нормативных материалов по сейсмическому строительству подземных трубопроводов.

2. Дифференциальное уравнение продольного движения подземных трубопроводов при сейсмических воздействиях запишется [5]:

 

 

2

u

(x,t)

aT2

2

u

(x,t)

b2

u

(x,t)

2

u

0

(x,t)

,

(1)

 

 

 

t 2

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

t 2

 

 

где

a2

EF ; b2

DH kx . k

x

– коэффициент равномерного сдвига трубо-

 

T

m

m

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

провода относительно грунта [1, 2, 3], EF – продольная жесткость, DН

– на-

ружной диаметр, m – масса трубы единицы длины, u0 (x,t) – движения грунта при землетрясениях, u (x,t) – относительное перемещение сечение трубы относительно грунта при сейсмических нагружениях. Значение kx – хорошо

изучено, оно сходно коэффициенту равномерного сдвига фундамента Сх – относительно грунта [6].

В качестве граничных условий принимаем, трубопровод защемленным

левым концом, а правый конец свободен:

 

u

0,t 0

 

 

 

 

 

 

 

u

 

 

 

 

 

 

 

(2)

EF

 

 

 

 

0

 

 

 

 

 

 

 

x

 

 

 

 

x l

 

 

 

 

x l

 

 

 

 

 

 

 

 

Начальными условиями будет

 

 

 

 

u

(x,0)

u

(x,0) 0 ,

(3)

575

Современные геотехнологии в строительстве и их научно-техническое сопровождение

 

 

Для

 

проведения

 

 

 

расчетов

 

принимаем

два

вида

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сpt

 

 

 

 

 

t

 

 

x

 

 

 

x,t

16

 

 

2

 

 

 

 

u

0

(x,t) Ae

 

C p sin t

 

,

u

0

0

,

 

. Здесь

A

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

Cp

 

 

 

 

0

 

 

 

 

 

 

амплитуда колебания грунта; – коэффициент затухания;

– частота коле-

баний; C p – скорость продольный волны в грунте.

 

 

 

 

 

 

 

 

 

 

Выводим безразмерные параметры:

u

DН u ,

 

 

 

 

x l x . Тогда

t

t0

t ,

дифференциальное уравнение (1) получает вид:

 

 

 

 

 

 

 

 

 

 

2u

 

2u

 

 

 

 

 

 

 

 

 

 

 

 

 

(4)

t 2

x2

b1u b2u0 .

 

 

 

 

 

 

 

 

 

 

 

 

 

b2l

2

 

 

 

 

 

 

l 2

 

 

Здесь введены следующие обозначения:

 

b

, b

 

 

 

. Граничные

 

a2

 

a2 D

 

 

 

 

 

 

 

 

 

1

 

 

 

2

 

H

условия (2) получает вид

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

DH u 0,t 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

DH u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(5)

EF

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x l

 

 

 

 

 

 

 

 

 

 

l

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для решения краевой задачи (4)–(5) с нулевыми начальными условиями применяется метода конечных разностей второго порядка точности. Полученная система алгебраических уравнений решается по явной схеме,

с удовлетворением условия Куранта h4 . При этом уравнение (4) получает вид:

j 1

 

2

j

 

2

 

2

 

j

 

2

j

 

j 1

 

2

j

ui

 

 

ui 1

(2 2

 

 

 

b1 )ui

 

 

 

ui 1

ui

 

 

 

h2

h2

 

 

h2

 

 

b2u0,i .

Краевые и начальные условия принимают вид:

u0j 0

 

 

 

 

 

 

 

DH (uij 1 uij 1 )

 

 

 

 

 

EF

 

 

0

 

 

 

 

 

 

 

 

 

 

x l

l

2h

 

 

 

 

x l

 

 

 

 

 

 

 

 

 

(6)

(7)

u

 

 

t 0

 

 

(8)

t 0 ui,0 ; u

ui,0

Напряжение подземного трубопровода вычисляется по формуле

 

E D u

E D

1

uij 1 uij 1

(9)

 

 

l x

 

l

2h

 

Здесь uij искомые относительные перемещения, h – шаг по координате, l

длина трубопровода, – шаг по времени.

Составлен компьютерный алгоритм и реализована программа на «Borland Delphi 7» на ПЭВМ.

576

Раздел 3. Аналитические и численные методы исследований оснований и фундаментов

Задача 1. Подземный трубопровод левый конец защемлен, а правый конец свободен. Действует сейсмическая нагрузка, изменяющаяся во време-

ни u0 (t) Asin t .

Механические и геометрические параметры приняты в следующем ви-

де.

E 2 106 кГ/см2;

 

F 2 R см2;

0,5 см;

DH 27 см;

R 13,5 см;

8 10 6 кГс2/см4;

 

 

m F кГс2/см2;

T 0,05 с;

 

l 10000 см;

A 2 см;

 

 

 

 

T

 

 

2

 

 

t t

 

 

(t)

 

 

 

3

 

 

с;

 

 

при u

0

.

 

 

kx 2,5 кГ/см ; t0

2

T

. u0 (t)

 

 

при 0

 

 

 

 

 

 

 

 

 

 

t t0

 

 

 

 

Результатырешениякраевойзадачиприводитсяввидеграфикарис. 1 и 2.

 

 

 

 

 

u0=Asinωt

 

 

 

 

 

 

 

 

u0=Asinωt

 

 

 

 

 

 

 

 

t=0.0015

t=0.00225

 

t=0.003

 

t=0.00375

 

 

t=0.0015

t=0.00225

 

t=0.003

t=0.00375

 

 

 

0,06

 

u, см

 

 

 

 

 

 

 

 

 

180

σ, кГ/см2

 

 

 

 

 

 

 

0,05

 

 

 

 

 

 

 

 

 

 

 

130

 

 

 

 

 

 

 

 

0,04

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,03

 

 

 

 

 

 

 

 

 

 

 

80

 

 

 

 

 

 

 

 

0,02

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,01

 

 

 

 

 

 

 

 

 

 

 

30

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x,

 

 

0

 

 

 

 

 

 

 

 

 

 

x,

-20 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20

40

60

80

100

 

 

 

 

0

20

40

60

 

 

80

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 1. Изменение перемещения по длине

 

Рис. 2. Изменение напряжения по длине

 

 

 

 

 

трубопровода

 

 

 

 

 

 

 

трубопровода

 

 

 

 

 

Задача 2. Подземный трубопровод левый конец защемлен, а правый

конец свободен. Действует сейсмическая нагрузка изменяющая по закону

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

(x,t)

 

 

 

 

C

 

50000 см/с, остальные исходные данные при-

0

Asin t

 

 

p

 

 

 

 

 

C p

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

няты как в задаче 1. Результаты решения краевой задачи приводится в виде

графика рис. 3 и 4.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u0=Asinω(t-x/cp)

 

 

 

 

 

 

 

u0=Asinω(t-x/cp)

 

 

 

 

 

 

 

t=0.0015

t=0.00225

 

t=0.003

 

t=0.00375

 

 

t=0.0015

t=0.00225

 

t=0.003

t=0.00375

 

 

 

 

 

 

 

 

180 σ, кГ/см2

 

 

 

 

 

 

 

0,05 u, см

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,04

 

 

 

 

 

 

 

 

 

 

 

150

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,03

 

 

 

 

 

 

 

 

 

 

 

120

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

90

 

 

 

 

 

 

 

 

0,02

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

60

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,01

 

 

 

 

 

 

 

 

 

 

 

30

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

100x,

0

 

 

 

 

 

x,

 

 

0

 

20

40

 

60

 

 

 

80

-30 0

20

40

60

80

100

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 3. Изменение перемещения по длине

 

Рис. 4. Изменение напряжения по длине

 

 

 

 

 

трубопровода

 

 

 

 

 

 

 

трубопровода

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

577

 

 

 

 

 

 

 

Современные геотехнологии в строительстве и их научно-техническое сопровождение

На основании полученных графиков можно заключить, что в первой задаче перемещения в разных временах изменяются до 50 м далее остается постоянными, а напряжение за переделом 50 м принимает значение нуль. Во второй задаче пределы изменения перемещений и напряжение доходит до 60 м. перемещения в данном времени остаются постоянными, а напряжения принимают значения нуль.

 

kx=2.5

kx=3.5

 

kx=4.5

kx=5.5

0,05

um ax, см

 

 

 

 

0,04

 

 

 

 

 

0,03

 

 

 

 

 

0,02

 

 

 

 

 

0,01

 

 

 

 

 

0

 

 

 

 

x,

0

20

40

60

80

100

Рис. 5. Изменение перемещения по длине

трубопровода t 0.00375,

 

 

 

kx=2.5

kx=3.5

kx=4.5

 

kx=5.5

160

 

σm ax, кГ/см2

 

 

 

 

120

 

 

 

 

 

 

80

 

 

 

 

 

 

40

 

 

 

 

 

 

0

 

 

 

 

 

x,

-40

0

20

40

60

80

100

Рис. 6. Изменение напряжения по длине

 

трубопровода, t 0.00375,

 

 

 

 

x

 

 

 

 

x

 

u

 

 

 

u

 

 

 

 

 

 

 

 

0

 

 

 

 

0

 

 

 

 

 

 

C p

 

 

 

C p

На рис. 5 и 6 приведены графики изменения перемещений и напряжений при разных значениях kx.– взаимодействия трубы с грунтом по длине трубопровода при t=0.00375.

Рассмотрим расчет полимерных трубопроводов что ныне очень актуально. При этом решается следующее дифференциальное уравнение

 

2

u

E

2

u

 

3

u

 

 

 

2R

k x (х)(

u

 

u

0 ) ,

(10)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

x 2

x 2

 

 

 

(R 2

r 2 )

t

t

 

 

 

 

 

 

 

 

 

 

 

где – плотность материала, Е – модуль упругости, – коэффициент вязко-

сти, kх (x) – коэффициент продольного взаимодействия трубопровода с грун-

товой средой, R – внешний радиус трубопровода, r – внутренний радиус трубопровода, u0 – смещение грунта. Дифференциальное уравнение (10) решается методом конечных разностей второго порядка точности. В результате полученная система алгебраических уравнений решается по неявной схеме. При этом используется метод прогонки. После определения перемещения будет возможность определения НДС полимерных трубопроводов при сейсмических нагружениях.

Задача 3. В качестве граничных условий принимаем (2). Начальные условия взяты как в (3). Геометрические, физические и механические пара-

метры полимерной трубы и грунта следующие: E=5·107кГ/м2, ρ=94 кГ·с24,

L=20 м, Rн=0.110 м, RВ=0,1037 м, Т=0,3 с, η=7,7·104 кГ·с/м2, а0=5,2·10-3 м.

Движение почвы задается импульсом и бегущим по направлению оси трубы:

578

Раздел 3. Аналитические и численные методы исследований оснований и фундаментов

 

u0 x,t

16

0

2

Сpt

, где Ср – скорость продольных сейсмических

04

,

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

волн в грунте, как приведено в работе [3].

 

 

 

 

 

 

 

 

 

 

 

Результаты численной реализации приведены на рис. 7–10 в виде рас-

пределения перемещений и продольных напряжений по безразмерной коор-

динате.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1000000 u (106)

 

 

 

 

 

1200000

u (106)

 

 

 

 

 

 

 

 

 

 

900000

 

 

 

 

 

1000000

 

 

 

 

 

 

 

 

 

 

 

800000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

700000

 

 

 

 

 

800000

 

 

 

 

 

 

 

 

 

 

 

600000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

500000

 

 

 

 

 

600000

 

 

 

 

 

 

 

 

 

 

 

400000

 

 

 

 

 

400000

 

 

 

 

 

 

 

 

 

 

 

300000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

200000

 

 

 

 

 

200000

 

 

 

 

 

 

 

 

 

 

 

100000

 

 

 

 

t

0

 

 

 

 

 

 

 

 

 

 

х

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

2

 

4

6

8

10

12

14

16

18

20

-100000 0

 

500

1000

1500

200

 

 

 

 

 

 

 

1-j=500 2-j=1000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 7. Изменение перемещений по време-

Рис. 8. Изменение перемещений по коорди-

 

 

 

ни при i=50

 

 

 

нате 1 – j=500, 2 – j=1000

 

 

 

30000 σ (кГ/м2)

 

 

 

 

 

100000 σ (кГ/м2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

80000

 

 

 

 

 

 

 

 

 

 

 

20000

 

 

 

 

 

60000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10000

 

 

 

 

 

40000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

20000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

0

 

 

 

 

t

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-20000 0

2

4

 

6

8

10

12

14

16

18

20

0

 

500

1000

1500

200

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-10000

 

 

 

 

 

-40000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-60000

 

 

 

 

 

 

 

 

 

 

 

-20000

 

 

 

 

 

-80000

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

-100000

 

 

 

 

 

 

 

 

 

 

 

-30000

 

 

 

 

 

-120000

 

 

 

 

j=1106

j=1098

j=500

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рис. 9. Изменение напряжений

 

Рис. 10. Изменение напряжений по коорди-

 

 

по времени при i =25

 

 

нате при j=500, j=1098, j=1106

 

 

Импульсивная нагрузка действует в интервале t=0-1000, с возрастанием нагрузки соответственно увеличиваются продольные смещения, (в рассмотренном сечении i=50), см. рис. 7. С уменьшением нагрузки соответственно уменьшаются перемещения до нуля (t=1000), рис. 7. После снятия нагрузки колебания затухают по времени (t≥1000).

На рис. 8 показано изменение смещения по длине трубопровода при конкретном времени j=500 и j=1000. Из рис. 8 видно, что при х=0-2 смещение резко увеличивается, а потом стабилизируется.

На рис. 9 приводятся изменения продольного напряжения по времени в сечении трубопровода i=25. При этом максимальная амплитуда продольного напряжения достигается значения до σ=±25·103 кГ/м2. При t≥1000 возникают напряжения с высокой частотой и максимальной амплитудой. Далее с течением времени уменьшаются амплитуда напряжения (см. рис. 9).

579

Современные геотехнологии в строительстве и их научно-техническое сопровождение

На рис. 10 показано изменение напряжений по длине трубопровода. При воздействии нагрузки максимальные значения напряжения возникает около защемленного контура при j=500 (σ=-10·104кГ/м2). После снятия нагрузки j=1098 и j=1106 максимальные значения напряжений возникают около защемленного края. В дальнейшем напряжение уменьшается до нуля.

В настоящее время применение полимерных труб в сейсмических районах стало актуальным, проводимые, в частности нами, исследования позволяют создать рекомендации по их расчету на сейсмостойкость.

Расчеты показали, что во многих случаях укладки трубопроводов при установлении подвижности концов трубопроводов, влиянием вторых концов можно пренебречь.

Таким образом, при нормировании сейсмических нагрузок на концевых сечениях достаточно рассмотрение напряженно – деформированного состояния трубопроводов около узла или вблизи узла.

Это позволит существенно упростить расчеты и дать рекомендации программ расчета в зависимости от грунтовых условий, глубины укладки и геометрических параметров рассчитываемых трубопроводов.

Литература

1.Аронов Р.И. Исследование условий взаимодействия трубы и грунта при продольных перемещениях трубопровода // Труды ВНИИСтройнефть.– М., 1953. – вып.V. –

С. 14 – 45.

2.Емельянов Л.М. О продольных напряжениях в подземных газопроводных трубах // Труды ВНИИ природного газа: Вопросы добычи, транспорта и переработки природ-

ных газов. – М., 1951. – С. 177 – 213.

3.Рашидов Т.Р. Динамическая теория сейсмостойкости сложных систем подземных сооружений. – Ташкент: Фан, 1973. – 180 с.

4.Рашидов Т.Р., Хожметов Г.Х. Сейсмостойкость подземных трубопроводов. –

Ташкент: Фан, 1985. – 152 с.

5.Уразбаев М. Т. Сейсмостойкость упругих и гидроупругих систем. Ташкент, 1966. 254 с.

6.Баркан Д.Д. Динамика оснований и фундаментов. – М.: Стройвоенмориздат, 1948. – 412 с.

580

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]