Файл: Чайлдс Э. Физические основы гидрологии почв.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 10.04.2024

Просмотров: 218

Скачиваний: 5

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

ГЛАВА 16

Конформное преобразование и годограф. Роль глубины залегания водоупора

16.1. Сущность конформного преобразования

Если одну величину, например X, называемую зависимой перемен­ ной, можно выразить с помощью уравнения в виде функции другой величины, например р, называемой независимой переменной, то ее можно также выразить как функцию третьей величины ѵ, если из­ вестна зависимость между ц и ѵ. Переход от уравнения, выража­ ющего связь между X и ц, к уравнению, выражающему связь между X и V, называется преобразованием. Цель подобных преобразований всегда состоит в том, чтобы получить уравнение, поддающееся реше­ нию. С рядом таких примеров мы уже встречались, в частности в пара­ графе 12.6 [уравнение (12.26)].

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

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

связывающее одну

комплексную величину с другой, например

 

W (р) == Ф -f іф,

(16.1)

где

(16.2)

р = x-\-iz.

Мы видели, что когда W определено, это решение позволяет по­ строить на диаграмме в координатах х и z эквипотенциали и линии тока. Такую диаграмму обычно называют плоскостью р. На примерах, приведенных в главе 14, было показано, как отыскать граничные условия данной задачи, когда определено решение, и отмечалось, что обычно требуется найти решение для заданных граничных усло­ вий. Иначе говоря, мы можем сначала построить в плоскости р гра­ ницы, а затем пытаться найти решение в виде сетки эквипотенциалей и линий тока. Границы могут быть сложными и не полностью извест­ ными. Например, в случае задачи с дренажем форма поверхности уровня грунтовых вод вначале неизвестна.

Переход от одной комплексной переменной р к другой комплекс­ ной переменной, скажем а, где

* “ / ( р ) = £ + *?,

( 1 6 . 3 )

приводит к другой диаграмме с осями £ и |, которую можно рассчи­ тать, если задано р, т. е. х и г, так что для каждой точки на плоскости


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

Преобразование, выраженное уравнением (16.3), имеет то свой­ ство, что форма любого малого элемента границы остается неизмен­ ной, меняются только его размер и ориентация1. Таким образом,

 

 

 

небольшой

 

квадрат

на

 

 

 

плоскости р остается квад­

 

 

 

ратом и на плоскости о,

 

 

 

хотя размер

и

положенно

 

 

 

его

могут

измениться.

 

 

 

Следовательно,

гидроди­

 

 

 

намическая

сетка

между

 

 

 

источником и стоком в пло­

 

 

 

скости р остается гидроди­

 

 

 

намической сеткой с той же

 

 

 

формой ячейки

в

плоско­

 

 

 

сти а

и продолжает удо­

 

 

 

влетворять

закону

Дарси

 

 

 

при том же значении К.

 

 

 

Величина

угла

поворота

 

 

 

и изменение размеров

не­

 

 

 

большого

элемента могут

 

 

 

быть различны

в

разных

 

 

X

частях

диаграммы,

по­

Рис. 16.1. Гипотетическая

гидродинамиче­

этому,

чтобы

соответство­

ская сетка одного угла проводящего прямо­

вать изменившейся

форме

угольного

листа.

границ в преобразованной

Плоскость р, где

р =

ж -f- iz.

плоскости,

сама гидроди­

 

 

 

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

женная на рис. 16.2, есть

линия, параллельная оси функции

тока,

а линия тока — прямая,

параллельная оси потенциалов.

Таким

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

ред.


образом, в плоскости W гидродинамическая сетка представляет со­ бой равномерную квадратную сетку, образованную ортогональными семействами прямых. Из свойства конформных преобразований сле­ дует. что гидродинамическая сетка в плоскости р также образована семействами ортогональных кривых и представляет собой совокуп­ ность ячеек, которые при очень малых размерах тоже будут квадрат­ ными. Этот результат следует непосредственно из решения уравне­ ния Лапласа, рассмотренного в параграфе 14.5 и в Дополнении 34.

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

границы могут быть весьма

 

 

 

 

сложными,

как,

например,

 

 

 

 

в случае, описывающем тече­

 

 

 

 

ние

воды,

поступившей на

 

 

 

 

поверхность

грунтовых

вод,

 

 

 

 

к каждой

из

нескольких

 

 

 

 

равноотстоящих

дрен,

зало­

 

 

 

 

женных на

одинаковой

глу­

 

 

 

 

бине (рис. 14.1 и 15.5). Тем

 

 

 

 

не менее при каждом последу­

 

 

 

 

ющем

преобразовании

гра­

 

 

 

 

ницы

могут

становиться все

 

 

 

 

проще и проще, пока в конце

 

 

 

 

концов

задача

не сведется

 

 

 

 

к течению от точечного источ­

 

 

 

 

ника к точечному стоку на

 

 

 

 

бесконечной

плоскости.

Ре­

 

 

 

 

шение для этого

случая хо­

 

 

 

 

рошо

 

известно.

Подобная

Рис.

16.2.

Конформное

преобразование

серия

преобразований

была

рис.

16.1

на плоскость

W, где W =

проведена Густафсоном

[81].

 

 

= Ф + іф.

 

Решение можно снова преоб­ разовать к плоскости р и получить его в первоначальной постановке.

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

Часто используется преобразование, при котором некоторый многоугольник одной плоскости превращается в прямую из преобра­ зованной плоскости. Особенно удобно, когда этой прямой является одна из осей. Такое преобразование называют преобразованием Шварца — Кристофеля. Оно состоит в следующем.

Если граница на плоскости р представляет собой многоугольник с вершинами в точках р1; р2, . . ., р„, внутренние углы при которых


равны соответственно a lt а 2, . . а„, то преобразование

р = С,1 j (о —

(а— (ц)(«./я-і) . . .

(а —p„)(a'l/It"1)d0 + C?

 

 

(16.4)

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

с г = р + іѵ ,

то точки, соответствующие вершинам рх, р2 и т. д. на плоскости р, являются точками р,1( ц2 и т. д. на действительной оси плоскости о. Сг и С2 — это константы, которые находят путем выбора определен­ ных соответствий между обеими плоскостями, поскольку на оси р, можно произвольно расположить до трех соответствующих точек. Если одной из них произвольно приписать значение бесконечности, то соответствующий член в уравнении (16.4) станет равен единице.

Доказательство справедливости уравнения (16.4) приведено в До­ полнении 38.

16.2. Годограф

Если частью границ задачи является поверхность уровня грунто­ вых вод, то, прежде чем приступить к преобразованию плоскости р, предстоит преодолеть существенное препятствие. Оно состоит в том, что поверхность уровня грунтовых вод является границей, форма которой вначале неизвестна, а потому для нее нельзя задать х и z, а следовательно, и р. Это препятствие преодолевают тем, что гра­ ницы изображают на диаграмме, на осях которой отложены не сами х и z, а горизонтальная и вертикальная составляющие скорости по­ тока в каждой точке, например соответственно ѵх и ѵг. Такая диа­ грамма называется годографом. На плоскости годографа, как пока­ зано в Дополнении 39, уровень грунтовых вод представится как окружность с известными центром и радиусом, хотя доля полной окружности, которую занимает сама поверхность уровня грунтовых вод, как правило, вначале неизвестна. Но и эту неопределенность можно устранить, задав часть окружности и получив затем решение, из которого будет видно, какое преобразование граничных условий необходимо при данном задании части границы. Иначе говоря, про­ извольный выбор части окружности годографа является как бы частью определения данной задачи, и эта часть раскрывается в окон­ чательном решении.

Можно указать следующие изображения границ задач на пло­ скости годографа. Линия тока, образующая угол Ѳ с горизонталью, характеризуется следующей зависимостью между ѵх и ѵг\

vxp z = tgQ.

На диаграмме годографа эта зависимость изобразится прямой с наклоном Ѳ, равным наклону линии тока на плоскости р, но проходя­


щей через начало координат. Если граничная линия тока есть пря­ мая с наклоном Ѳ, ее годографом является параллельная прямая, проходящая через начало координат.

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

=

а годограф представляет пря­ мую с наклоном я/2 + Ѳ, про­ ходящую через начало коор­ динат, т. е. перпендикулярную соответствующей границе в про­ странстве диаграммы.

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

N

может являться

потоком,

поддерживающим

испарение

с

поверхности. Будучи отри­

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

(л/2 + Ѳ),

Рис. 16.3. Уровень грунтовых вод в поч­ ве с влагопроводностью К при восходя­

щем к поверхности вертикальном по­ токе N изображается на плоскости годо­ графа ѵх + іѵг полуокружностью.

ѵ%н- к + 0K - N ) ß Y = {(Z + N ) ß ) \

(16.5)

Из этого уравнения следует, что годограф представляет собой окружность радиуса + N)J2, центр которой находится на оси ѵг, на расстоянии — 7Ѵ)/2 ниже начала координат. Как показано на рис. 16.3, эта окружность пересекает ось ѵг в точках +ІѴ и —К. На той же диаграмме показана линия, проведенная от пересечения в точке (О, —К) к точке (ѵх, vz) на окружности, которая соответ­ ствует точке на уровне грунтовых вод, где наклон равен Ѳ. Сама эта линия также образует угол Ѳ с осью ѵг.

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