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

книги из ГПНТБ / Мирцхулава, Ц. Е. Надежность гидромелиоративных сооружений

.pdf
Скачиваний:
20
Добавлен:
22.10.2023
Размер:
11.64 Mб
Скачать

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

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

Среднюю частоту срыва частиц (отказов) в единицу времени можно установить по выражению

 

v / ? / P

( R d ) d R d ;

(188)

 

 

 

 

t

 

Ф (Rd)= - 7 =

exp

Г

( R d - R d ) s

(189)

 

2°Rd

/ 2

n'Rd

 

 

 

' Проводя выкладки, аналогичные изложенным в пункте 2 на­ стоящей главы, получаем и для прогноза устойчивости размыву несвязных грунтов следующую зависимость:

 

аxv

 

 

Д

(191)

Т —

 

 

 

%

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

= 1

 

2 1 п — = In Р

+

(192)

Тн

\ т 7

 

 

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

120

Используя эту характеристику, уравнение (185) устойчивости частицы, расположенной на дне, с учетом условия работы коэф­ фициентом т можно записать в виде

- ^

д = R d,

(193)

т

а

 

где т — коэффициент условий

работы

(в данном случае он учи­

тывает влияние наносов в коллоидном состоянии на протекание процесса размыва) [76].

Выполняя выкладки, аналогичные приведенным

в работе

[76], получим

 

 

 

^дядоп —1,25 | /

0 Д 7о71н (Тч

Т,)d.

(194)

'^я.доп — (lg

у

(v4-yo)d.

(195)

Используя характеристику надежности т]н, можно установить допускаемые (неразмывающие) скорости течения водного потока для несвязных грунтов, которые обеспечат безотказную работу размываемого основания с заданной надежностью и вероятно­ стью Р jipn сроке воздействия размывающего потока Т, числе вы­

бросов v T , дисперсии

напряжения отрывка %v и сопротив-

Д

4

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

 

Обеспечить 100%-ную надежность устойчивости против размы­ ва основания теоретически невозможно. Однако при довольно вы­ сокой вероятности Р можно практически гарантировать работу основания без размыва.

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

В заключение следует отметить, что предложенная методика прогноза начала движения наносов с использованием приемов теории надежности (с вероятностной постановкой проблемы) хо­ тя и несовершенна, но позволила построить достаточно простую стохастическую модель для оценки важного этапа эрозии. Рас­ смотренная стохастическая модель хорошо объясняет устойчи­ вость и трогание наносов. Анализ полученных результатов пока­ зывает, что прогноз процесса по выведенным зависимостям зна­ чительно точнее и объективнее, чем решение этих задач, получаемое по данным экспериментальных исследований критических ско­ ростей с использованием необъективных показателей — по срыву отдельных частиц или массовому движению.

121

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

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

4. УЧЕТ ВРЕМЕНИ ПРИ ЭРОЗИОННЫХ РАСЧЕТАХ

Прогноз скорости потока, при которой начинается передвиже­ ние частиц, составляющих ложе водотока, необходим для реше­ ния многих задач, связанных с эрозионными расчетами, проекти­ рованием каналов, спрямлений, берегозащитных и руслорегули­ рующих сооружений и т. д. За последние годы предложен ряд зависимостей, позволяющих достаточно точно учитывать факторы, связанные с решением подобных вопросов [73, 76, 77, 78]. Рас­ пространение этих зависимостей на процессы, где необходим учет времени воздействия потока на размываемое основание, не обо­ сновано.

Многочисленные .эксперименты движения донных наносов, вы­ полненные различными исследователями, а также автором, с фик­ сацией процесса скоростной кино- и микрокиносъемкой позво­ лили уточнить механизм движения донных наносов. В частности, скоростная кино- и микрокиносъемка позволила уточнить, что до отрыва при определенных значениях скоростей течения час­ тица вибрирует, раскачивается и лишь при появлении актуаль­ ных пульсационных скоростей, силовое воздействие которых больше сопротивляемости частицы отрыву, срывается (перекаты­ вается). Но для этого необходимо, чтобы продолжительность воз­ действия пульсационной скорости была достаточна для перево­ рачивания частицы вокруг точки перекатывания. В противном случае частица приподнимается и вновь возвращается на место расположения [73].

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

уравнением (для частиц, в которых не

появляются

обусловлен­

ные мелкозернистостью

силы сцепления)

[76]

 

( а д + /У а )= '0 в /2Яь

(196)

т

 

 

 

122

 

Рис. 23. Схема

перекатывания

частиц

под

воздействием

 

 

 

 

руслового

потока.

 

 

 

 

 

где

 

 

п — коэффициент

перегрузки,

учитываю­

 

 

 

 

щий пульсационный характер скоро­

 

 

 

 

стей

(следовательно, и нагрузки на

 

 

 

 

частицу);

 

 

 

 

(в дан­

 

 

 

т — коэффициент условий работы

 

 

 

 

ном случае он учитывает влияние на­

 

 

 

 

носов в коллоидном состоянии);

 

 

 

К\ — коэффициент

однородности,

учиты­

/i =

6 id; /2 =

6 2 с/; /3 =

 

вающий неоднородность наносов;

6 3 6 ? — плечи

лобового

усилия,

подъемного

 

 

 

 

усилия*-веса частицы в воде

(прихо­

 

 

 

 

дящегося на одно зерно) относитель­

Лобовое

 

 

но точки перекатывания.

[76].

результирующее усилие М.

А.

Дементьева

 

 

 

 

 

г/=

 

 

 

(197)

 

 

 

Рл УоАжСЦЙ2

;

 

 

 

равнодействующую подъемного усилия

 

 

 

 

 

 

 

 

Рв =

уоАуа2 е? 2 vl

;

 

 

 

(198)

вес

частицы

в воде

 

 

 

 

 

 

 

 

 

 

 

GB=

(уч—уо)с?3,

 

 

 

(199)

« 1 с?2;

a2d2— площади

миделя для

лобового

и подъемного усилий;

А*; Ау — коэффициенты лобового и подъемного усилий, зави­ сящие от зон движения и условий обтекания, уста­ навливаем по данным:

123

П о д с т а в л я я з н а ч е н и я ( 1 9 7 ) — (199)

в (1 9 6 ), и м е е м :

 

То щ у - ( 6 1^ 0 1 + 63X3/ 02) =

о

(уч— уо) d 8 2Ki .

( 2 0 0 )

т

zs

 

 

Следует считать, что движение частицы начинается, когда мо­ мент веса частицы в воде относительно оси перекатывания станет равным нулю. Очевидно, в этом случае при одностороннем движе­ нии потока прекращение силового воздействия потока не вызовет возвращения частицы на прежнее место, как это происходит при ее раскачивании. Частица возвращается на прежнее место, если величина силы воздействия недостаточная или продолжительность воздействия мала.

На основе опытных данных и наблюдений установлено, что шаровидная частица диаметром d выступает со дна (которое имеет наклон к горизонту а) на 0,7 d. Рассматривая устойчивость частицы, заключаем, что лобовая сила потока Рл действует па­ раллельно плоскости дна и стремится сдвинуть частицу вдоль течения, а подъемная Рв — оторвать частицу от дна и переместить ее в вертикальном направлении. Принимаем, что при предельном состоянии частицы Рв приложена в центре ее, а Рл — выше се­ редины выступа А= 0,7 d.

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

С этой

целью уравнение движения частицы (шара)

запишем с

учетом

приведенных выше зависимостей предельного

равновесия:

 

/ 5

= — (^ л /1+^ в /з)-С Е/2С03ф,

(201)

 

d№

т

 

/ — момент инерции шара относительно оси вращения (пере­ катывания) .

Интегрируя уравнение (201), имеем:

(202)

Постоянную С можно определить по начальным условиям, а именно до перекатывания частица неподвижна (рис. 23), то есть

 

 

^

= 0 .

 

 

 

 

 

dt

 

 

 

Тогда ^ = 0

и ф =

ф0 + а .

 

 

 

 

Обозначая

(Д + 1 +^в/з)

— = Р ,

из выражения

(202) полу-

чим:

 

 

т

 

 

 

r _

2 G + . .

.

2Р . . .

 

 

(203)

 

С

—- 81п(ф0+ а )

j

(фо+а).

124

П о с л е в н е с е н и я С в

(202)

и м е е м

 

— ( ф _ ф 0 _

а ) — _

_ [ s i n cp— s m ( c p 0+ a ) ]

(204)

В результате интегрирования (204) получаем

t =

 

dy

 

(205)

2Р

 

2 Gnl<2

V

( Т — ? о -

[ s i n tp — s i n ( c p o + a ) ]

 

Постоянную С можно установить также по начальным усло­ виям, когда ^ = 0 ; ф = ф о + а .

Выражение (205) HHTefpnpyeM приближенно, для этого разла­

гаем sin ф по показательному ряду при значении ф = ф0+ а . Если ограничиться первыми двумя членами, получаем

d?

2 а д

[sin (tp0+ a ) + ( t p — ср0— a ) COS ( c p o + a ) — Sin (<f>0+ a ) ] .

/

(206)

 

Для установления входящих в уравнение (205) коэффициентов лобового и подъемного усилий, а также других, используем зави­ симость для определения допускаемых (неразмывающих) скоро­ стей течения [75, 29, 48]

Ъ „ . = '

.

2

5

(207,

и ее вид с нераскрытыми коэффициентами [76,

82]

 

1,25

2 g m

( у ч — V o ) 82d K ] ;

(208)

 

V

Зто(81аЛ л-+8з“2^у)и

 

 

 

6 2 = 0,5; K i = L

 

 

Приравнивая (207) и (208), будем иметь

 

 

6iai7x+ 6 3 a27?/= 0,0733.

 

 

Когда t = 0, ф= ф0 + а ,

поэтому

С = 0.

 

 

При ф= 90°----- получаем

 

 

 

 

 

 

;

(209)

 

 

 

 

( 210)

125

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

 

 

0,3663'f4d6

-?и—а

'0,0733 f 0d3n v ^

(211)

 

-

2 ~gm

- 0,2617('c4—f0)<73 cos (<f0+«)

 

 

Если при разложении sirup ограничиться тремя членами, ре­ зультаты расчета по этой зависимости практически не будут

отличаться от

(205).

(211) относительно донной

скорости, имеем

Решая уравнение

v м

I

2 gm

0,3663

 

 

 

O,0733-f0d3n

gt*

 

 

+0,2617(Тч-

То)й?*cos (<?„+«)],

(212)

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

 

Когда уклоном дна ввиду малости можно пренебречь, для принятой высоты Выступа А= 0,7 d значение угла ср = 23035/ и принимая приближенно у чф о т = 1 ,6 5 тс/м3,

При предельном состоянии, когда ид^ = г>д ДОп,

t =

(217)

126

5. ОЦЕНКА НАДЕЖНОСТИ СОПРЯГАЮЩИХ СООРУЖЕНИИ ПРИ МЕСТНОМ РАЗМЫВЕ

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

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

Глубину местного размыва [76, 78] для связных грунтов уста­

навливают по зависимости

 

 

 

t =

(

v .

- 7. 5&0)

+ 0,25# н.б.

(218)

 

\

/ 1—0,175 ctg3

 

 

 

 

Чдоп.В

" 1

 

 

Величину

1 >д доп.в

допускаемой

неразмывающей

скорости в

воронке размыва с учетом характеристики надежности устанав­

ливаем по формуле [76]

 

 

доп.в— 1,25 \ /~

~^

т)и

■[уч Т о )^+ 1,25(С " + Р в + Р ц ) ].

у

0,3т0

у

где т)н— характеристика надежности, учитывающая'как пульса-

ционный

характер

скорости, так

и

неоднородность

сцепления

грунтов

размываемого

основания;

 

g — ускорение

свободного

падения;

 

 

 

т — коэффициент условий

работы;

 

 

 

уо— плотность

воды;

 

 

 

 

 

уч— плотность

агрегатов;

 

 

 

 

d — средний диаметр агрегатов;-

 

 

 

С”— нормативная усталостная прочность на

разрыв;

 

Р о— давление, оказываемое

струей на

подлежащие отрыву

агрегаты и слагающие дно воронки;

 

во­

Ри— пригружающее действие, обусловленное глубиной

ды в воронке размыва.

 

основания

гид­

Для сохранения

устойчивости

подмываемого

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

Я > / ,

(219)

127

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

t — глубина размыва, обусловленная как сопротивляемо­ стью грунта основания размыву, так и скоростью потом­ ка во входном сечении.

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

ется следующим образом:

 

P ^ P ( H ^ t ) .

(220)

Вероятность Р непревышения фактического

размыва зависит

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

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

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

как известно,

что величина максимального

значения возрастает

с увеличением

длины отрезка (вероятность

превышения данного

значения переменной зависит от числа наблюдений за этой пере­ менной) .

Для получения интересующего нас ряда (или наиболее полной характеристики случайной величины — закона распределения) рас­ четных величин глубины местного размыва можно использовать современные методы теории вероятностей и математической статис­ тики (закон распределения экстремальных величин Гумбеля [22] и др.).

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

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

Задачи о «выбросах» случайной функции за некоторый уро­ вень, который позволяет находить максимальное значение орди­ наты, за определенный промежуток времени наиболее просто мож­

128

но решить при подчиненности процесса нормальному распределе­ нию.

Если данные наблюдений отклоняются от нормального про­ цесса, для нахождения наибольших значений и вероятности их появления можно использовать методы статистики экстремальных значений [122]. Этот прием позволяет по ограниченному ряду максимальных значений, полученному из наблюдений, оценить возможную максимальную величину при последующих наблюде­ ниях за более длительный срок службы сооружения. Рассматри­ ваемый процесс, как_известно, полностью харак4еризует матема­

тическое ожидание t, дисперсия Dt и корреляционная функция Rt(т). Эти параметры не изменяются с течением времени.

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

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

пт

Р = —— е~пТ,

(221)

га !

 

где п — среднее число выбросов от предела Я

кривой процесса

t(T)-

 

m — число превышений предела Я за время Т.

Когда отсутствует превышение, то есть т = 0, нет и отказа,

тогда зависимость примет вид

 

р = е ~ пТ.

(222)

Если допустить, что процесс нормальный, среднее число вы­ бросов от предела Я можно установить по формуле Райса [156]

 

 

п = п

ехр

-(H-t)2

 

 

(223)

 

 

2D,

 

 

 

 

 

 

 

 

 

 

 

Внося

значения

п по формуле

(222)

в

(223),

получим

 

H = t

У 2РД1п (Й 7)-1 п (—1пЯ)]

 

(224)

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H=tr\n,

 

 

 

 

 

 

У2Р,[\п (пТ)— 1 п ( —

I n Р)\

 

(225)

 

 

■— 1 +

 

7

 

 

 

 

 

 

 

 

 

 

 

где г)н — характеристика надежности,

зависящая

от

дисперсии

глубины

местного

размыва

(которая, в свою очередь,

зависит от

изменения скорости ниспадающей

струи

во входном

сечении, а

9 З а к а з 6767

129

Соседние файлы в папке книги из ГПНТБ