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

книги из ГПНТБ / Картвелишвили Н.А. Потоки в недеформируемых руслах

.pdf
Скачиваний:
7
Добавлен:
25.10.2023
Размер:
12.62 Mб
Скачать

При этом, конечно, U{Fi = UzF2- Исключая из этих уравнений U, получаем уравнение для определения глубины затопления Я:

FF*,

[ - \ ( ^ ) + Х ( Л ) + 0\

- U, }2=g{Fy

-

Я2 у2 ).

(35.16)

F-Fo

 

 

 

 

 

 

 

 

 

Если F =АНп,

 

то это уравнение имеет вид

 

 

 

 

- ^ ± І Я " { 2VgT,(VlT1-VW)

 

+ Ul

-

L/2

}2 =

 

а для русла прямоугольного сечения (/г=1)

при

t/i = t72 = 0

будет

вн(УЖ-У

я ) 2 = ( я 2

- я 2 2 )

[(^)2 -1],

 

или, полагая 2 = УЯ,

 

 

 

 

 

 

 

2 Е - WoZ4

+ 1 6 Я 2

V~H[ z3 -

(8Я; +

Я2 )

Я 2 г 2 + Я 2 = 0 .

 

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

I

а: —~и'

%

ІірЩіЩПЩііЩШіЩШШШЩШПШШІЇЩїШ

Рис. 63. Волны при частичном разрушении плотины.

тины образуется отверстие. В результате в

некоторый момент

после разрушения

будем иметь

следующую

картину

(рис.

63).

В верхнем бьефе движется волна понижения,

причем

скорость { ] '

установившегося течения на участке CD, по аналогии с (35.14),

равна

( 7 ' = - Ц Г ) + ) , ( ^ ) + ( У ,

 

(35.17)

 

 

В нижнем бьефе движется прерывная волна повышения, причем

по аналогии с (35.15),

 

 

 

 

 

F"-F2

g

У -^У*

(35.18)

 

 

Расход Q через

плотину

в общем случае зависит от глубин

Н'

и Н" слева и справа

от плотины, т. е.

 

 

 

Q = ? ( # ' ,

H")=U'F'=U"F",

 

(35.19)

причем, если истечение иезатопленное, то функция ф зависит только

от Нг

и не зависит

от Н". Исключая U' и U"

с помощью

(35.19)

из

(35.17)

и (35.18),

получаем систему

двух уравнений для

опре­

деления Н'

и Н":

 

 

 

 

 

 

 

 

 

 

^ ' { - А ( / 7 ' ) + М Л ) + с / 1 ) =

с р ( Я ' ,

И"),

 

 

7J^{-^r[-X(F')

 

+ X(Fl) + (Jl} -

 

U2\=g(F'y'-F2y2).

 

При незатопленном

истечении (через

водослив

или через

отвер­

стие

в плотине), когда

ф не зависит

от Н", решается сначала пер­

вое

из этих уравнений

и найденное

значение Н'

подставляется во

второе уравнение, из которого определяется Н". При затопленном истечении функция ф обычно такова, что первое уравнение разре­ шимо в явной форме относительно Н". Подставляя Н" во второе уравнение, получаем одно уравнение относительно Н'.

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

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

параграфе

(средняя скорость

остается

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

постоянной

глубине),

нельзя определить

момент

(или место)

разрушения

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

36. Об определении исходных характеристик речного русла

для гидравлических расчетов

Для расчета волнового режима необходимо предварительно за­ даться характеристиками русла, т. е. функциями Хо, Іи Х2, цо, j-ii, v0 , vi в уравнениях (33.8) и (33.9). Если известны такие характе­ ристики русла, как F, В, К, то построение указанных функций сво­ дится к подбору таких коэффициентов в аппроксимирующих вы­ ражениях

Х = л 0 + Х 1 / г - | - л 2 А 2 ,

которые давали бы наилучшее (например, в смысле способа наи­ меньших квадратов) приближение к кривым K(h), \x(h), v(h) для различных створов. По этим данным нетрудно построить искомые графики зависимости Ко, Ki, Ко, LI0-, І-Ч, vo, vi от s.

Однако вопрос о том, что именно представляют собой харак­ теристики F, В, К для речного русла неправильной формы, не так прост. Он связан с выводом уравнений движения потока в русле неправильной формы и не может быть решен вполне строго из-за того, что гидравлическая идеализация для таких русел более груба, чем для русел с плоскостью симметрии, течение в которых имеет меньше особенностей, не охватываемых гидравлической идеализа­ цией. В главе IV на форму русла было наложено ограничение, со­ стоявшее в том, что русло имеет вертикальную плоскость симмет­ рии. В таких руслах динамическая ось потока и линия центров тяжести сечений лежат в указанной плоскости симметрии. В даль­

D

Рис. 64. Фарватер и динамическая ось в русле симметричной не из­ меняющейся по длине формы.

нейшем уравнения, полученные для симметричных русел, перено­

сились, без специальных

оговорок,

на речные русла неправильной

формы

(в тех случаях,

когда не

учитывалось влияние кривизны

струй

в вертикальной плоскости).

Такой перенос означает прежде

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

Рассмотрим русло, образованное при движении некоторой сим­ метричной кривой (образующей, см. рис. 15) по другой кривой (направляющей) так, что плоскость образующей всегда нормальна к направляющей, ось симметрии образующей всегда лежит в верти­ кальной плоскости, а вершина скользит по направляющей. При этом предполагается, что радиус кривизны направляющей всегда больше, чем В/2. Если динамическая ось потока в таком русле ле­ жала бы на поверхности, описываемой осью симметрии образую­ щей (в плане линия АВ на рис. 64, т. е. фарватер), то очевидно, распространяя на такой поток уравнения, выведенные для потока с вертикальной плоскостью симметрии без членов, зависящих от кривизны линий тока осредненных скоростей, мы не сделаем ошибки ни в уравнении неразрывности, ни в динамическом урав­ нении. Справедливость этого утверждения в пределах, в которых приемлемы предположения об отсутствии влияния сил, вызывае­ мых кривизной струй, и о горизонтальности сечения свободной по­ верхности потока плоскостью образующей, вытекает из теорем Гульдена об объеме и поверхности тел вращения. Течение в таком русле будет (в пределах действия указанных предположений) та­ ким же, каким оно было бы в русле той же длины и уклона с на­ правляющей, расположенной в вертикальной плоскости.

Но в действительности динамическая ось потока не лежит на поверхности, образуемой осью симметрии образующей, она от­ клоняется в сторону вогнутого берега (рис. 64, линия CD). Поэтому возникает вопрос, как правильнее отсчитывать длины: по линии ЛЯ или по линии CD, а значит и определять сечение потока, которое должно быть нормально к линии отсчета длин. В случае правиль­ ного русла результаты расчета мало зависят от того, какой из этих линий отдать предпочтение.

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

Рис.

65.

Фарватер и

стрежень

""-«у^ j

^

'

 

в

речном

русле.

 

 

 

/ — меженный фарватер,

2 — стре­

 

 

 

жень

в

паводок,

3. 4 — границы

 

 

 

 

пойменного

русла .

 

 

 

ках с широкими поймами стрежень и фарватер в межень направ­

лены

в общем

одинаково (рис. 65, линия

1), в паводок же

фарва­

тер

остается

прежним

(если пренебречь

деформациями

русла),

а стрежень

(линия 2) сильно спрямляется. Поэтому если в

паводок

отсчитывать

расстояния

по фарватеру, а сечения брать нормально

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

Рассмотрим два варианта описания паводка уравнениями гид­ равлической идеализации:

1) расстояния отмеряются по стрежню

дНс

__ Q c

, Uc

дЦс

| 1

dUz

дх*

I<1

S

' <?*c

~t~ g

dtc

dxc ' dtc

2) расстояния отмеряются по фарватеру

 

дНф

Ql

иф

диф

і

1*

дхф

~

g '

< * с ^

g

 

 

дхф *

діф

 

 

16*

 

 

 

 

243.

Каждому из этих вариантов соответствуют свои характеристики русла.

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

НфнНс, Qb=mQQc,

mQ

 

і(j,=т І 4,

І*Ф=от/4.

 

где т — масштабные

коэффициенты.

Подставляя эти значения

в уравнения второго варианта, получаем:

 

 

mm*

дН,

т г т

о

<?с ,

 

с

< * с

тнт\

К\

1

 

<7С

_ dUc

• "lx'nQ

. dUc

 

g

' дхс

' mtmH

'

dt^ '

dxc 1 HI^HQ

Для подобия нужно, чтобы эти уравнения совпали с уравнени­ ями первого варианта, т. е. чтобы было

пцтх = тхт\ = т\ ^ mxmQ = тхт2н ^

Так как расходы и наполнения

русла в обоих вариантах оди­

наковы, то mq=mH = l и, следовательно,

т н

т-к

щ

или

 

 

m~\\mx,

mt=mx,

mK=Vmx.

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

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

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

Практические приемы подготовки исходных данных для расче­ тов течений в речных руслах по гидрометрическим и топографиче­ ским данным подробно рассматриваются в работе [21].

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

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

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

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

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

мени, сил и средств, чем обычные изыскательские работы.

Послед­

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

участках

очень большой протяженности.

 

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

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

рактеристика

возможного пути подхода к ее решению.

 

Напишем

уравнения (33.8) и (33.9), причем

не будем

учиты­

вать для простоты инерционные члены

 

 

X, (s) Л + 2Хо (s)

+ Х 2 (s)

2Х, (s) +

(s) и*=0,

(36.1)

 

 

ТЯ- + Ъ®-Ъ-+ът-£-=0.

 

(36.2)

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

уровня по длине

L бьефа за время

Т получена

функция

h(s, х)

в виде семейства кривых или в виде таблиц, или, наконец, в виде

аналитической аппроксимации, а также функция и (0, х).

Требуется

по этим данным определить коэффициенты

%о, /ч, /U, vo, vi уравне­

ний (36.1) и (36.2).

 

 

 

Из уравнения (36.2)

имеем

 

 

u(s, , ) = " ( 0 , ^

J К ( / ) + ^(/)Л(Л

т)] d h £ х ) dl.

(36.3)

Обозначим через Ф результат подстановки u{s, т) из формулы (36.3) в левую часть (36.1):

Ф=Х,(5 )Л.(5 >

,) + X2 (s)/f-(5 , t ) + d h ( s

d s z ) +

+ 2[X0 (s)+X1 (5)A(5, -с)] {«(0, ^ ) - | [ v o ( 0 +

vi(/)A(/, *)] X

х _ ^ о . Л | + Х о ( 5 ) | „ ( о ,

т ) - | к ( о + ^ ( О л ( ^ .

 

 

 

(36.4)

Если «(О, т)

и A(s, т)

заданы экспериментально, то Ф не будет

тождественным

нулем ни при каких Хо, Xi, Xz, vo, Vi. Поэтому все

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

і

і

 

j

j<I>2dsrfc = min.

(36.5)

о о

Задача отыскания коэффициентов уравнений по решению этих уравнений некорректна и поэтому на искомые функции Хо, Хх, Xz, va, vi следует наложить дополнительные условия, вытекающие из смысла задачи. Так как нас интересуют сглаженные осредненные характеристики русла, то необходимо, чтобы эти функции были аналитичны. Можно пойти еще дальше и искать их не в классе аналитических функций, а в более узком классе полиномов от s некоторой степени N, т. е. полагать

X J ( s ) = a i 0 + a J 1 s + . . . -r-aiNsN,

/ = 0 , 1, 2,

 

vt(s)=bl0+bns+

. . . +biNsN,

i=0, 1.

(36.6)

Тогда все дело сводится к определению коэффициентов

ац и

Ьц, при которых выполняется условие

(36.5), т. е. к решению си­

стемы алгебраических

уравнений:

 

 

і і

datj

и о

і і

о о

В частном случае, когда речное русло аппроксимируется руслом призматическим (что очень часто дает хороший результат), т. е.

когда Ко, Ki, Ко, vo, vi не функции s, а постоянные величины, эта система будет

И « - £ - * * - < > .

| І Ф - £ - * * - О .

о о

о о

1

где в первом уравнении і = 0, 1, 2, а

во втором — і = 0, 1. Выписы­

вать в развернутом виде эти уравнения мы не будем.

Слабые места определения характеристик русла по наблюде­ ниям над неустановившимся движением состоят главным образом

вследующем.

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

разных средних расходах и рассматривать Ко, Ki, К2, vo, vi как функ­ ции s только в пределах каждого данного опыта. Вообще же эти функции зависят еще и от среднего расхода или среднего уровня при производстве опыта.

2.Нельзя указать какой-либо критерий для определения сте­ пени N аппроксимирующих полиномов (36.6), эту степень нужно находить пробами.

3.Вычисления характеристик связаны с дифференцированием экспериментально определяемой функции h(s, т) и поэтому не мо­ гут быть особенно точными.

Глава IX

Н Е О Д Н О М Е Р Н Ы Е И Д Е А Л И З А Ц И И Р У С Л О В Ы Х П О Т О К О В

37. Двумерная плановая

идеализация

 

Уравнение (15.3), в котором мы

отбросим члены,

связанные

с компонентами скорости и ускорения

по направлению

Хз, и урав­

нение (15.14) после несложных преобразований дают следующую

систему уравнений двумерного

потока:

 

 

 

 

 

 

 

,

.

1

дН

\

.*

.

1

dUx

а — 1

Ui

dH

 

м + е з ( / і Ч

 

;

) = У7 1 і +1

g

 

dt

g

Н

1

dt

г

 

 

 

 

дхх

)

 

 

 

 

 

+т(-

 

Ux

dUx

. U2

 

dUx

\

a

U\

QL

 

(37.1)

 

L\

dxx

1

I.2

 

dx2

j

g

L\L2

dxx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(r

\

1

d H

\

. *

і

1

dUо

a —I

U2

 

дН .

 

Ч + Ч

{ ^

+ - ^

-dxYJ—Л

+

^ -

д

Г

g

7f ' ~Ш~ " T "

,

 

 

 

_dU2_

.

Щ_

 

dU2

\

a_

U\

dLx

 

,„ »

^

g

\

L2

" dx2

" f

I ,

'

длг, J

£

' I , Z . 2

' d x 2

'

v3'"*'

Вместо Ui и U2

введем

новые

неизвестные

qi = LzHUi

и q2 =

--LiHUz- Уравнение неразрывности

 

в новых

переменных

будет

 

^ Г

+

-

Й

- +

^

Т = °

.

 

07.4)

Имея далее

 

 

 

 

 

 

 

 

 

 

dUx

 

1

(

1

dqx

4\

dH УI

 

dt

 

L2

\

H

dt

 

 

dt

>1 >

 

dUx

 

1

dqx

 

Я\

 

dH

1i

dL2

 

dxx = "Zo7HК дхх

 

H

 

dxx

L2

dxx

 

dUx

1

(

dqx

 

4\

 

dH

4x

dL2

 

dx2

L2H

\ dx2

 

H

 

dx2

L2

dx2

 

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

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