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

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

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

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

Добавлен: 10.04.2024

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

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

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

г.Интегрирование дифференциального уравнения.

Предстоит интегрировать следующее уравнение (16.37):

i ( l + N/K)dp--

 

dQ - -С4Х

(A.2f2 1) ctg Ѳ'

)Л ;

 

{ (Д40.18)(16.37)

 

2Kt

 

 

dQ/dt можно определить, дифференцируя уравнение (16.35) по t:

 

 

dQ __/ L ( M ~ N ) \ f

2t

\

 

 

 

 

(Д40.19)

~dT~~\

Я

K

 

J \р 2

 

Д/ 2 2

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ctg Ѳ' == 2 А,/(Ä,2 1).

 

 

 

 

(Д40.20)

Подставляя (Д40.20) и (Д40.19) в уравнение (Д40.18), получаем

 

Ц і + N/ K)

dp = -

 

2L

 

 

 

-J ). F (t) d t - d Q ,

 

(Д40.21)

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Mz — i

 

 

 

 

 

 

 

F(t)=- (/a—fj.a/A.2) («a—i)

 

 

 

 

 

Функцию F (г) можно привести к виду

 

 

 

 

 

 

 

 

 

 

1

■ ){(^)Ы

 

 

 

 

 

 

 

\

(/

1 —П2

\

/

1

t + |т Д )+

 

2(і—и2А2) J \\

нА

j

\

t —нА

 

 

+<“ -1>(т= г -

т£ г )}-

 

 

 

 

 

Подставив это выражение f

(t) в уравнение (Д40.21), получим

 

И і+ ІУ/К) »

- ЛИ- (

~ ^

 

) ( A T - T T T W

 

 

. /

 

 

 

1/у -і. w

 

 

1_________! _

ъ

, .

(Д40.22)

‘ \

яТС

/ѵУVЯ,,—1/А. /V) \

іі—нА

t+иА

/

 

 

 

Из уравнений (16.23), (16.24) и (16.21) имеем

І з т у г = ~(1+ У)= - ( * +М)/(М -ЛГ).

Подставив это выражение в уравнение (Д40.22), получим

i ( l + N / K ) d p = —dQ-

L (M — N)

) ( A r - T f r ) dt -

яК

L ( M + K)

\ / ___ 1_

f - j - ' P / Ä , 'jdt.

л к

J \V t —^ н

Это уравнение можно интегрировать непосредственно:

где А — постоянная интегрирования.

362


Теперь из уравнения (16.35) можно подставить выражение для самого Q. После чего получаем окончательную форму уравнения (Д40.22):

L ( M - N ) i ( l + NI K) p = - H f - c f ( l + N l K ) л К

((12/Я2_/2) С- (/2-1)

L (М + К )

In ( - S ï 7 r ) + '4-

 

+ 1п(ігг)]+- пК

w 40-23’

Известно, что в точке Q значение t бесконечно велико (рис. 16.12), а на пло­ скости р оно находится в точке L + ici. Подставив эти значения в уравнение

(Д40.23), найдем

 

A = I L (1 + NI K) + H f - — t ~ ?0- ln ( - 1 ).

 

 

 

Подставив это выражение в уравнение (Д40.23) и деля

на і (1 -f* N /К),

получаем

iL ( M - N )

 

(/2—Ц2/А.2)

 

 

 

 

 

 

 

н "(т+г)]'

 

 

p = L + icf-

n ( N + K)

ln

(<2 —О

 

 

 

 

 

L

 

 

 

 

 

 

 

iL(Hf + iC) ,

/ f + n A

\

 

 

 

(Д40.24)

 

 

 

n ( N + K )

^ / —[i/X

J -

 

 

 

 

 

 

 

 

 

 

Напомним, что, согласно уравнению (16.21),

 

 

 

 

 

поэтому

 

 

( K + N ) / ( M - N ) = у,

 

 

 

 

 

 

(M + K) / ( N + K) = l + i/y.

 

 

 

 

 

 

 

 

 

 

Из уравнения (16.26) определяется ß:

 

 

 

 

 

 

 

 

 

 

p / X =

l +

ß .

 

 

 

 

 

Подставив эти значения в уравнение (Д40.24), получим после некоторых

алгебраических преобразований:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

г + 1

\-1

{(Д40.25)

 

 

 

 

 

 

 

 

 

 

J J - 1

(16-38)

W можно выразить из уравнений (16.27) и (16.35) для Q:

 

 

 

W —Q^- i NpJK = H f - i - C f ( l Jr N J K

) -

L (M — N )

ln

(t + l + P) (t —1 —p)

 

 

 

 

 

 

л К

 

 

(/+ 1)(/-1)

 

, iN

'

л

^ln

f+

1 -Р

I

ln

t+1

ß)]•

 

 

к

L“ '

V“'

1 + ^ß

 

г + 1 +

 

 

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

 

 

 

 

 

 

 

 

W = d> + i ^ =

H f + C f + ^ ^ - + - ^ -

j^CAf—JV) ln ( t - i ) - M l n

( f —

1

— ß ) -f-

+

l V l n (

« +

( M - N ) (K — N)

(

t + 1

 

 

(Д40-26)

l +

ß )

 

 

 

z + l + ß

 

 

(16.39)

 

 

 

 

(K + N)

 

 

 

 

д. Условия на границах в плоскости р.

Рассмотрев плоскость t на рис. 16.12, видим, что вдоль границы QR вели­ чины t действительные, положительные и превышают 1 + ß. Следовательно,

все логарифмы в уравнениях (Д40.25)—(Д40.26) — это логарифмы

действитель­

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

нивая мнимые части, найдем

t

( 1 6 . 4 5 )

nz / L = яс,/L + ln ( * + | + j) HIM t + 1 + ß)•

*+1

(

(Д40.27)


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

x = L,

как это и следует из того, что известно о границах на плоскости р. Аналогично, из действительных частей уравнения (Д40.26) имеем

яФIL =-- (Hf -\- c f ) - \ - ~ £(Af — N) ln (t —1)М ln (t 1 — ß) +

+ A l n ( « + i +

(MßHN) ( K — N)

(

t + i

)

I

(Д40.28)

 

K + N

V

і + l + ß

J

I

(16.46)

Поскольку Ф есть сумма Я и z, Я можно найти, вычитая z из Ф с исполь­ зованием уравнений (Д40.28) и (Д40.27). В результате

M — N

In (т+г) +

М + К

( 1 + Ü J L ) .

{ (Д40-29)

л Н / А = nHfJL - К

К In

 

(16.47)

На границе DP величины t действительные и положительные, но заключены между нулем и единицей. Следовательно, как t — 1, так и t — 1 — ß — отри­

цательные числа, поэтому можно написать:

In

( * - l - ß

C+ln (

l+ ß-* N

Vf + i + ß

l + ß + < J 9

ln (t — 1) = гя + ln (1 t)

и

ln (t 1 — ß) = г'л + ln (1 + ß—t),

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

.i* /L = n Cf/L +

l n ( | ± i - p f ) +

( | ) l n (

1 + t

\

/

(Д40.30)

ß + 1+ f

/ ’

\

(16-42)

л Ф / Ь = ~

(H[+cf) + ~

[ ( М - А )

l n ( l - / ) - M l n ( ß + l - / )

+

 

+ l V l n ( ß + l + /)

,

( М - - N) (К — N)

/'

1 + t

 

 

(Д40.31)

1

K + N

 

1П1 ß + 1 + 1

 

 

(16.43)

 

 

 

 

 

 

 

n H/ L = nH[ JL+

M — N

ln

(

\_4_ Æf +

Jf

/• ß + 1 + i

 

f

(Д40.32)

K

1 + t ) '

к

 

1 П 1. ß + 1 - t

 

\

(16.44)

 

 

 

\

 

 

Известно, что в точке D на этой границе,

а именно на оси дрены, z = О,

a t = 1. Подставляя эти величины в уравнение (Д40.30), получаем важный ре­

зультат:

ЛCf

 

 

+ (ѵ ) ln ß+ 2

= 0,

 

 

~т~ +

1п- ß+ 2

 

 

или

 

 

 

 

J (Д40.33)

ncfIL

= ln (l +

2 /ß) + (2 /v) ln (l +

ß/2 ).

1

(16.40)

 

 

 

 

 

Комбинируя это равенство с уравнением (16.32), получаем дополнительное

выражение для Ь/, а именно

 

 

 

{^

 

n b f l L = ln (1 + 2 /ß) + (2/ Y) ln (i Yqqr) -

1 6 .4 1 )

Наконец, на границе

DS величины t действительные

и положительные

и заключены между 1 и 1 +

ß. Поэтому требуется только преобразовать лога­

рифм t 1 — ß к виду іл

+

ln (1 + ß — t). Результатом на этот раз явятся

уравнения (16.48)—(16.50).

 

 

 

 

 

 


ГЛАВА 17

Некоторые трехмерные дренажные гидродинамические сетки

17Л. Приближение Дюпюи—Форхаймера для трехмерных задач

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

трем прямоугольным координатам. Точной

теории, сопоставимой

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

в этом случае не суще­

ствует, однако, обратясь к приближенной

теории Дюпюи — Фор­

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

Как и в случае двухмерных задач, основное допущение теории Дюпюи — Форхаймера состоит в том, что поток ограничивается гори­ зонтальным направлением, и, следовательно, дрены доходят до водо­ упора, а эквипотенциальные поверхности вертикальны, и значения потенциала на них соответствуют высотам, на которых они пересе­ кают уровень грунтовых вод. Горизонтальное направление потока будет различным в разных точках. Горизонтальный поток имеет составляющие в направлении двух горизонтальных прямоугольных осей координат х и у. В соответствии с законом Дарси эти составля­ ющие скорости ѵх и ѵу пропорциональны составляющим наклона уровня грунтовых вод dZJdx и dZ/dy, поскольку, в соответствии с до­ пущением Дюпюи — Форхаймера, последние являются составля­ ющими градиента потенциала.

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

Qx = —KZôy(dZ/dx),

тогда как в параллельном сечении на координате х + ôx величина потока равна

Qx+bx = - К ôУ {Z оэгідх) + {д [Z (dZ/dx)]/dx} ôx}.


Первое из этих двух уравнений описывает приток в призму, а второе — отток из нее, поэтому скорость изменения влагозапаса

(dS]dt)x, связанного с изменчивостью расхода в

направлении х ,

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

 

dSldt)x — КЬу {д [Z (dZ/dx)]/dx}Ьх.

(17.1а)

Точно так же для изменения расхода в направлении у имеем

(dS/dt)y = KÔx{d [Z(dZ/dy)]/dy}ôy.

(17.16)

Далее, при наличии на поверхности осадков, имеющих интенсив­ ность q, связанная с ним скорость влагозапаса (dS/dt)qбудет равна

(dS/dt)q — qbxdy.

(17.1в)

Очевидно, что результирующая скорость влагозапаса,

связанная

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

dS/dt = (ôx by) [q+ K{d [Z (,3Z/dx)]ldx} + K{d[Z (dZ/dy)]/ dy}]. (17.2)

В частном случае, когда устанавливается стационарное состояние, скорость изменения влагозапаса равна нулю, и уравнение (17.2) после подстановки l]^dZ2]dx вместо Z (dZ]dx) и такой же подстановки для направления у принимает вид

â*Z2ïdx* + d*ZP/dy* - 2 q /K .

(17.3)

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

17.2. Круговая периферическая дрена

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

nr*q = - 2 nrZK (dZ/dr),

(17.4)

где q — интенсивность осадков, a Z — высота уровня грунтовых вод в точке г. Если периферическая дрена представляет собой канаву, отстоящую на расстояние R от центра, причем вода в ней стоит на