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

Конспект лекций Портовые сооружения

.pdf
Скачиваний:
691
Добавлен:
16.03.2016
Размер:
7.44 Mб
Скачать

31

определять с учетом скорости ветра, а также разгона или продолжительности действия ветра по аналитическим формулам /12/:

hd

hd

0,16 Vgw2

0,16 Vgw2

1

1

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g L

 

 

 

 

 

 

3

 

 

 

 

 

 

 

1 6, 0 10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

Vw

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g t

0,635

 

 

 

 

 

3

 

 

 

 

 

1, 04

10

 

 

 

 

 

 

 

 

 

 

1

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где Vw – скорость ветра, м/с; g – ускорение свободного падения, м/с2; L – разгон, м; t – продолжительность ветра, с.

 

 

 

 

 

 

 

 

 

 

 

 

 

g L

 

g t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

значения hd и

По значениям

безразмерных параметров

 

 

 

 

и

 

 

следует вычислять

V

2

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w

 

 

w

 

 

 

 

 

 

 

принимать меньшую величину.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Средний период волн T , с, и среднюю длину d , м, следует определять по формулам:

 

 

 

g

 

d 0,625

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

V

h

 

 

 

 

 

 

 

 

 

 

g T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T 19,5

w

 

 

 

 

 

 

,

 

 

 

 

 

 

 

d

 

 

 

 

,

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

g

 

Vw

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Высоту волны i%-ой обеспеченности в системе

волн hi, м, необходимо определять по формулам: hi ki % hd ,

где ki –коэффициент перехода к высоте волны обеспеченностью i%.

Определение параметров волн в прибрежной мелководной зоне

Мелководная акватория

 

 

 

 

 

 

 

 

 

 

 

 

К мелководным акваториям

относятся водоёмы с глубиной d 0,5 d

и

уклонами дна

i 0,001 или практически горизонтальном дне.

 

 

 

 

 

 

 

 

Элементы волн в этом случае также

устанавливают по темже

графикам.

С этой целью

 

 

 

 

 

 

 

 

 

 

 

 

g L

g d

g h

g T

 

вычисляют безразмерные величины

 

 

и

 

, снимают с графика

 

и

 

 

и определяют

 

V 2

V 2

V 2

V

 

 

w

w

w

w

 

среднюю высоту h , средней период T и среднюю длину волн. Высоту волн 1% обеспеченности для

 

 

 

 

 

 

 

 

этой зоны вычисляют умножением средней высоты волн h на коэффициент k1% ,

который определяют

 

 

 

g d

 

g L

 

по графикам (рис. 1.4) в зависимости от безразмерных величин

 

и

 

,

принимая в расчете

V 2

V 2

 

 

 

w

 

w

 

минимальное значение коэффициента.

32

Зона трансформации волн

При распространении волн на открытой части моря в прибрежную зону меняются их структура и параметры. При подходе волн к прибрежной зоне под некоторым углом относительно направления изобат, что имеет место в реальных условиях морских побережий, с уменьшением глубины происходит постоянное изменение направления распространения волн, искривление линий гребней и разворачивание их в сторону берега, при этом гребни стремятся занять положение, параллельное изобатам. В результате гребень волн растягивается, и параметры волн изменяются. Это явление называется рефракцией. Объясняется тем, что мористый участок гребня обгоняет бережной, т.к. скорость распространения волн уменьшается с уменьшением глубины воды. Явление рефракции существенно влияет на изменение параметров волн.

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

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

разворачивается, и расстояние между лучами увеличивается до l lo . Вывод искомой зависимости основывается на предположении о том, что поток энергии между выделенными лучами постоянен /3, 5/:

hd 2 ld cd nd h2 l c n,

где h – высота составляющей; с – фазовая скорость волн; n – доля энергии, перемещающейся вместе с гребнем (символы с индексом «d» относятся к глубоководной зоне, символы без индекса – к мелководной). Из уравнения следует, что:

 

 

h2 h 2

k

r

k

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

t

 

 

 

 

 

 

 

 

 

 

 

где

kr

 

l

 

– коэффициент рефракции; kt

 

 

nd cd

- коэффициент трансформации.

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n c

 

 

 

 

 

 

 

 

Коэффициент kr зависит от угла d

 

 

 

 

 

 

d

kt - только от

 

 

и относительной глубины

 

, коэффициент

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

d

 

 

 

 

 

угла d

 

cd

 

 

 

глубины

 

 

.

При заданном значении

 

отношение

 

 

уменьшается с

уменьшением

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

c

 

 

 

nd

относительной глубины. Отношение с уменьшением относительной глубины увеличивается, так n

как на глубокой воде с гребнем перемещается половина волновой энергии, а в очень мелководной зоне (в зоне разбивающихся валов) – вся энергия /13/. Ход изменения коэффициента kr зависит от формы изобат – при прямолинейных и вогнутых (например, в бухте) изобатах коэффициент kr уменьшается с уменьшением относительной глубины, так как волновые лучи расходятся. В случае схождения волновых лучей (например, в области мыса) коэффициент kr увеличивается с уменьшением

относительной глубины. В целом изменение общего коэффициента kr kt зависит от конкретной

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

Определение параметров волн в зоне трансформации начинается с построения диаграммы рефракции на заданном плане с изобатами. Около изобат выписываются абсолютные d и относительные

d

глубины d . Перед построением плана рефракции, полагая, что малые колебания глубины не влияют

33

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

Между сглаженными изобатами проводятся средние линии, отвечающие глубине dср 12 di di 1 .

На глубокой воде в соответствии с направлением расчетного волнения наносятся волновые лучи, расстояние между которыми выбираются в зависимости от требуемой точности построения плана рефракции. Положение лучей выбирается таким образом, чтобы после рефракции эти лучи попали в район будущего строительства порта. Практически следует брать 5-6 лучей, а расстояние между ними 1- 2 длины волны на глубокой воде с последующим сгущением. Построение ведется для каждого луча

независимо от других Поскольку влияние дна начинает сказываться на глубине d0 0,5 d , то

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

d0ср

d0

d1

 

. В точке пересечения

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

исходного

луча

со

средней

 

изобатой d0cp поворачивается

на угол

r

, значение которого

 

 

 

 

 

 

 

 

 

 

 

 

0

 

определяется уравнением:

 

 

 

 

 

 

 

 

 

i

 

 

sin ri

 

 

,

 

 

 

 

 

 

 

 

sin(

ri

 

ri 1

)

 

 

 

 

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

где ri - угол схождения или расхождения волновых лучей, т.е. угол поворота луча в точке изменения глубины, который откладывается всегда в сторону уменьшающихся глубин; r - угол рефракции, т.е.

угол между лучом и нормалью к средней изобате в точке пересечения её лучом; i - длина волны на соответствующей изобате.

i

Отношение i 1 вычисляется по уравнению:

i

 

thkidi

,

 

 

thk

d

i 1

 

i 1

 

i

 

 

где k 2 .i

34

Рефракция волн на мелководье

Нумерация изобат ведется от глубины к берегу.

Вместо уравнения можно пользоваться номограммой Снипа. Построение диаграммы рефракции в этом случае ведется в соответствии с указаниями СНиП /7/. Каждый волновой луч, начиная от

изобаты d0 0,5 d продолжают до промежуточной между первой и ближайшей к ней изобате и

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

меньших глубин на угол ri , величина которого определяется в соответствии со схемой и графиком, приведенными на рис. 6 прил. I СНиП /7/.

Схема и графики для построения плана рефракции

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

Лучи проводятся до критической глубины dcr , после которой они не меняют направление.

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

hi % kt kr kl ki hd ,

где kt - коэффициент трансформации; kr - коэффициент рефракции; kl - обобщенный коэффициент потерь; ki - коэффициент перехода к высоте волны i%-ой обеспеченности; hd - исходная средняя

высота волны на глубоководной акватории.

Коэффициент рефракции определяется по формуле:

35

kr aad ,

где ad – кратчайшее расстояние между смежными волновыми лучами в глубоководной зоне; a

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

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

Коэффициент трансформации принимается по кривой 1 (рис. 5 прил. I CНиП /7/), в зависимости

d

от относительной глубины .d

Обобщенный коэффициент потерь kl

должен определяться

по

заданным

значениям

относительной глубины и уклонам дна по табл. 1.4 прил. I СНиП /7/. При уклонах дна 0,03 и более

следует принимать значение обобщенного коэффициента потерь равным 1.

 

 

 

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

kr при фиксированном

 

 

 

 

 

 

значении d зависит только от глубины в

расчетной точке и

генерального

направления

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

коэффициент рефракции для участков дна с прямолинейными изобатами (при 0 d

60 ,

град.)

следует вычислять по формуле:

 

 

 

 

1 C

2

sin2

d

 

 

 

 

 

kr

 

N

 

 

,

 

 

 

 

cos

2

d

 

 

 

 

 

 

 

 

 

 

 

 

 

где d – угол между направлением волн и нормалью к изобатам в расчетной точке; CN

– параметр,

определяемый по формуле.

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

d

графикам, представленным на рис.4 прил. I СНиП /7/ при заданных безразмерных параметрах

 

 

и

 

 

 

 

 

 

 

 

 

 

 

d

h1%

g T 2 , при этом период волн принимается равным периоду волн в глубоководной зоне или по

аналитическим зависимостям. При этом следует отметить, что длина волн определяется только

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

th0,5 2 d , d d

Определение параметров волн по данной методике ведется до критической глубины, то есть до глубины первого обрушения. Начиная с dcr , элементы волн определяются как в прибойной зоне. При

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

Прибойная зона

36

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

Прежде чем приступить к определению элементов волн в прибойной зоне, необходимо определить глубину первого dcr и последнего обрушения dcr ,u . Критическая глубина с учетом рефракции волн определяется методом последовательных приближений. Полагая в первом

приближении

 

hcur

1%

 

 

h

i

 

, находят по

(рис.5 прил. I СНиП /7/)

с учетом уклона дна i параметр

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

g T

 

 

 

 

g T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

cr

и глубину dcr I

первого приближения.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ход расчета можно представить в следующем виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1-ое приближение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

i

 

 

hглd ,1%

 

рис.Снипа в зависимости от уклона i dcr I .

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g T

 

 

g T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2-ое приближение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h1

 

 

 

 

 

на глубине d

 

 

I

определяется k , k , k и

h1

k

 

k

 

k

 

k

 

hd , далее

1%

 

по рис.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cr

 

 

 

t l

r

1%

 

t

 

l

 

r

1%

 

 

 

g T 2

 

 

 

 

 

Снипа в зависимости от уклона i dcr2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3-е приближение:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

на глубине dcr2 kt, kl, kr hcr2 dcr3.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Расчеты проводят до тех пор пока d i cr

d i 1cr .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Как правило, бывает достаточно 3-4 приближений. Сходимость предыдущего и последующего приближений можно считать удовлетворительной, если разность между ними не превышает 3%.

Глубину первого обрушения волн dcr , м, в прибойной зоне при постоянном уклоне i

допускается определять по формуле:

Критическая глубина, соответствующая последнему обрушению волн при переменных уклонах дна, определяется путем последовательного расчета критических глубин в прибойной зоне для участков дна с постоянными уклонами. Для определения уклонов необходимо построить профили дна по рассматриваемым волновым лучам. Первая критическая глубина, как было изложено выше, определяется методом последовательных приближений. Вторая и последующая критические глубины определяется путем умножения глубины предыдущего обрушения на коэффициент, принимаемый по таблице 1.7 в зависимости от уклона того участка, где происходит предыдущее обрушение.

 

 

 

 

 

 

 

 

 

 

 

Таблица 1.7

 

 

 

 

 

 

 

 

 

 

 

 

 

Уклон дна i

 

0,01

0,015

0,02

0,025

0,03

0,035

0,04

0,045

 

0,05

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициент ku

 

0,75

0,63

0,56

0,5

0,45

0,42

0,4

0,37

 

0,35

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Критическую глубину,

соответствующую последнему обрушению волн dcr,u

при постоянном

уклоне дна, следует определять по формуле:

dcr,u kun 1 dcr ,

37

где ku - коэффициент, принимаемый по табл. 1.7; n

– число обрушений (включая первое),

принимаемое из ряда n=2, 3, и 4 при выполнении неравенств: kun 2

0, 43 и kun 1

0, 43.

При определении глубины последнего обрушения

dcr,u

коэффициент ku

или произведение

коэффициентов не должны приниматься менее 0,35.

 

 

 

 

При уклонах дна более 0,05 следует принимать значение критической глубины dcr

dcr,u . При

уклонах дна 0,001 i 0,01 необходимо принимать:

 

 

 

 

dcr,u 0, 43 dcr ,

 

Высоту

волны в

прибойной зоне hsur1% , м, необходимо определять для заданных уклонов дна i по графикам СНиПа. При

этом по безразмерной величине

dsur

 

принимается значение

hsur1%

и соответственно определяется

 

 

 

 

 

 

g

 

2

 

 

d

T

hsur1% .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Длина волны в прибойной зоне

 

sur ,

м, определяется по верхней огибающей кривой на рис.

СНиПа в зависимости от параметра

 

dsur

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Высоту волны i %-ой обеспеченности в прибойной зоне, м, необходимо определять по формуле:

 

hsur1% ki hsur1% ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициент

перехода к высоте заданной обеспеченности ki

следует принимать по таблице 1.8.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 1.8

 

 

 

 

 

 

 

 

 

 

 

Обеспеченность высоты волны i%

 

0,1

1,0

 

2,0

5,0

13,0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Коэффициент ki

 

 

 

 

 

 

 

 

 

 

1,15

1,0

 

0,95

0,9

0,85

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Расчет параметров волн на защищенных акваториях порта

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

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

Волна, проникающая на огражденную территорию через ворота порта, подвергается ряду изменений. Прежде всего, происходит дифракция волн при выходе в порт, т.к. волна, огибая в плане всякое препятствие, имеет свойство разворачиваться и менять свое направление. При этом количество энергии, приходящееся на участок гребня (фронта) волны, равный ширине ворот, распределяется на большую длину ее гребня. Волна как бы растягивается при непрерывном удлинении гребня (фронта) волны. Это неминуемо приводит к уменьшению волновой удельной энергии (приходящейся на единицу

длины гребня волны), и к уменьшению высоты волн на огражденной акватории ( E h2 ). При этом уменьшение величины удельной волновой энергии, а, следовательно, и высоты волн происходит неравномерно. Линия, соединяющая самые высокие точки гребней волн, называется главным волновым лучом.

Помимо дифракции на волновой режим в порту при изменяющихся глубинах влияет и рефракция волн, также вызывающая их поворот и уменьшение энергии. Кроме того, волна затухает и

38

потому, что происходит гашение ее энергии в результате трения на границах водной среды (у дна, у стенок оградительных и причальных сооружений и т.п.). Наконец, весьма значительное влияние на режим в порту оказывает отражение волн от сооружений и берегов. Отраженные волны, двигаясь навстречу подходящим волнам, интерферируют с ними и создают сложную картину, которую обычно называют толчеей.

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

Измеряя высоту волн в воротах порта hвх и в разных точках на акватории hi , вычисляют значение

hi ki . Получаемый в ряде точек безразмерный коэффициент ki , называемый коэффициентом

hвх

остаточной высоты, дает возможность построить изолинии остаточной волны и выделить те районы акватории, которые удовлетворяют поставленным требованиям. Совокупность значений ki на

акватории порта иногда называется паспортом порта. Коэффициент ki в общем случае можно записать в виде:

ki kdif kr kl kt ko ,

где kdif – коэффициент дифракции волн; kr – коэффициент рефракции волн; kl – коэффициент потерь (потери энергии по длине, главным образом диссипация (рассеивание), в зависимости от глубины воды и уклона дна и изменяющийся с их уменьшением от 1,0 до 0,66); kt – коэффициент

трансформации: kt f (d ) (зависит от глубины, трения о дно); ko – коэффициент интерференции

(зависит от отражающей способности сооружений, ограничивающих акваторию, в основном - причальных).

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

Теоретическим исследованиям дифракции волн посвящены работы Ю.З. Алешкова, А.М. Жуковца, Ю.М. Крылова, Д.Д. Лаппо, И.А. Степанова и др. /5, 17, 18, 19, 20/. Однако многие теоретические решения, при всей их внешней стройности изложения, использовать для практических целей пока не представляется возможным, так как точность вычислений по ним недостаточна. Это объясняется тем, что дифракция волн жидкости является более сложным явлением, чем дифракция световых, звуковых и электромагнитных волн, аналогию с которой используют некоторые авторы. В этой связи особо хотелось бы отметить теоретические и экспериментальные исследования В.К. Завьялова /17/, по результатам которых были разработаны рекомендации по расчету дифракции волн, включенные в СНиП /7/.

Расчет дифракции волн в порту при защищенности акватории порта одиночным молом

При ограждении порта одиночным молом волна проникает в порт, огибая голову мола, и при определении интенсивности волнения на акватории порта необходимо рассматривать задачу дифракции волн одиночным молом. В основу решения этой задачи многие авторы /17, 19/ используют гипотезу Т. Юнга о дифракции, как результате передачи энергии вдоль фронтов волн, положение которых на акватории определяется методом Гюйгенса-Френеля. Количественное выражение поперечной передачи энергии условно принято в виде потока энергии q или количества энергии, проходящего через поперечное сечение отдельной волны за единицу времени. Согласно Т. Юнгу поток энергии

39

E

пропорционален длине волны , групповой скорости u и градиенту энергии на участке фронта l

волны l , т.е.:

 

 

 

q u

E

u gradl dE,

Подобное выражение

 

l

 

 

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

/17/.

Выражение для коэффициента дифракции акватории, огражденной одиночным молом в этом случае имеет вид:

 

 

 

 

 

 

 

1 2, 3

 

r 1,7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

kdif ,s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

r 1,7

 

r

0,67

 

 

 

2,5

 

 

 

0,17

r

 

 

 

1

2, 3

 

 

1,1

 

 

 

th

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tg

при 0;

 

 

 

 

 

 

 

 

 

 

 

где

f

при 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Полученное выражение представлено (рис.7 прил. I СНиП /7/) в виде номограммы, которая значительно упрощает и ускоряет выполнение расчетов.

В соответствии с рекомендациями СНиП /7/ высоту дифрагированной волны hdif , м, на огражденной акватории необходимо определять по формуле:

hdif kdif hi ,

где kdif – коэффициент дифракции волн; hi - высота исходной волны i%-ной обеспеченности (в

расчетную обеспеченность высот волн в системе при определении защищенности портовых акваторий необходимо принимать равной 5%).

В качестве расчетной длины принимается исходная длина на входе в акваторию. Коэффициент дифракции волн kdif для акватории, огражденной одиночным молом, следует

принимать в соответствии со схемой и графиками согласно штриховой линии со стрелками или рассчитывать по формуле. Для этого необходимо определить следующие параметры:

- , град - угол между границей волновой тени (ГВТ) и линией, соединяющей голову мола с

расчетной точкой. В зоне волновой тени угол считается положительным, вне зоны тени –

отрицательным; при этом следует отметить, что границей волновой тени (ГВТ) является линия, проведенная через голову мола параллельно исходному лучу;

r

- - относительное расстояние от головы мола до расчетной точки,

где r –расстояние от головы мола до рассматриваемой точки, м; - длина волны, м;

- , град - угол между границей волновой тени и линией мола, проходящей по внутренней грани оградительного сооружения.

40

Расчет дифракции волн в порту при защищенности акватории порта сходящимися молами

Решение задачи дифракции волн сходящимися молами производилось с использованием гипотезы Т. Юнга и баланса энергии, а также экспериментально /17/. Фронт волн на акватории определялся методом Гюйгенса-Френеля. Ось координат l =0 проходит через голову правого мола параллельного лучу исходных волн. На линии входа (начальное условие):

 

r

 

l

 

 

 

l

 

 

 

 

 

 

 

E(

 

;

 

)

 

 

Eo

i

 

,

 

 

 

 

 

 

 

 

 

r

0

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

B

 

l

0

 

 

 

 

i

 

,

 

 

 

 

 

 

 

 

где В – расчетная ширина входа на акваторию. Функция распределения энергии на акватории определялась суммированием основного и отраженных от каждого мола потоков волновой энергии, т.е.:

 

qотр.1

 

 

 

 

 

 

 

 

qc qo

lп1

 

r

1

 

 

 

 

 

 

 

qотр.2

 

 

 

 

 

0,

 

 

 

 

 

 

lп2

 

r

2

 

Местоположение

 

 

 

 

 

 

главного луча волнения, которым является линия, проходящая через точки фронтов волн, имеющих максимальные высоты, определялось исследованием на «максимум» функции энергии волн на акватории.

Полученные расчетно-эмпирические зависимости также были представлены графически и в виде номограмм (рис.8 и рис.9 прил. I СНиП /7/).