Книги / Теоретические основы гидравлики. Ртищева А.С. 2007 г
.pdf30
Отношение силы к площади поверхности, на которую она действует,
называется напряжением.
Напряжение называется нормальным σ, если оно направлено по нормали к анализируемой поверхности.
Напряжение называется касательным τ, если оно направлено по касательной к анализируемой поверхности.
Нормальное напряжение вдоль оси x (действующее на площадку, перпендикулярную оси x) обозначают σx. Аналогично σy, σz.
Касательные напряжения записываются с двумя индексами, например, τxy, первый индекс указывает, какой оси перпендикулярна выделенная площадка, а второй – в направлении какой координатной оси действует рассматриваемое касательное напряжение.
Таким образом, в направлении каждой координатной оси действуют нормальная и две касательных составляющих напряжения. Так в направлении осей x, y, z действуют следующие составляющие напряжения соответственно:
|
σx |
τ yx |
τzx |
. |
(3.40) |
|
τxy |
σ y |
τzy |
||
|
|
||||
|
τxz |
τ yz |
σ z |
|
|
Выражение (3.40) представляет собой тензор напряжений. |
(3.41) |
||||
τxy =τ yx ; τxz =τzx ; τzy =τ yz . |
|||||
Таким образом, тензор напряжений симметричен относительно главной |
|||||
диагонали |
|
|
|
|
|
|
σx |
τxy |
τxz |
. |
(3.42) |
|
τxy |
σ y |
τ yz |
||
|
|
||||
|
τxz |
τ yz |
σz |
|
|
В невязкой жидкости касательные напряжения равны нулю |
(3.43) |
||||
τxy =τ yz =τxz = 0. |
|||||
Также уравнение (3.43) справедливо для вязкой, но покоящейся жидкости. |
|||||
Для этих случаев также верно |
|
|
|
|
|
σ x =σ y =σz = −p , |
(3.44) |
||||
где p – давление, действующее на выделенный элемент противоположно по знаку напряжениям σ, поскольку давление создает напряжение сжатия, а
положительными принято считать растягивающие напряжения. |
|
||
Из формулы (3.44) следует, что |
(3.45) |
||
p = − |
1 |
(σ x +σ y +σ z ). |
|
|
|
||
3 |
|
||
Таким образом, в вязкой покоящейся или невязкой (покоящейся или движущейся) жидкости силы трения не возникают, напряжения определяются
31
гидростатическим давлением. Для этих случаев тензор напряжений принимает вид
|
− p |
0 |
0 |
|
|
|
|
(3.46) |
|
|
|
||||||
|
0 |
− p |
0 |
|
|
|
. |
|
|
0 |
0 |
− p |
|
|
|
|
|
Нормальные напряжения в движущейся вязкой жидкости можно выразить |
||||||||
зависимостями: |
= −p + ∆py ; σ z = −p + ∆pz ; |
(3.47) |
||||||
σ x = −p + ∆px ; σ y |
||||||||
где ∆px, ∆py, ∆pz – дополнительные давления в направлении координатных осей x, y, z, обусловленное влиянием вязкости.
Таким образом, закон трения Ньютона (1.4) для плоского потока можно обобщить на пространственное течение.
Для ньютоновской жидкости, у которой связь между напряжением и скоростью деформации линейна, такой обобщенный закон трения можно представить в виде
|
2 |
r |
∂w |
x |
|
|
∂w |
|
|
|
|
∂wy |
|
|
|||||
σx = −p − |
|
µdivW + 2µ |
|
|
|
|
x |
+ |
|
|
|
|
|
; |
|
||||
3 |
∂x |
|
; τxy = 2µεxy = µ |
∂y |
∂x |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
2 |
r |
∂wy |
|
|
∂w |
x |
|
|
|
∂w |
z |
|
|
|
|
(3.48) |
||
σ y = −p − |
|
µdivW + 2µ |
|
|
|
; τxz = 2µεxz = µ |
|
|
+ |
|
|
; |
|||||||
3 |
∂y |
|
∂z |
|
|
|
|||||||||||||
|
|
|
|
|
|
|
∂x |
|
|
|
|||||||||
|
2 |
r |
∂w |
z |
|
∂wy |
|
|
|
∂w |
|
|
|
|
|
|
|||
σz = −p − |
|
µdivW + 2µ |
|
|
|
|
|
|
|
|
|
|
z |
|
|
|
|||
3 |
∂z |
|
; τ yz = 2µεyz = µ |
∂z |
+ |
∂y |
. |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
r |
∂wx |
|
∂wy |
|
∂wz |
|
||
divW = |
|
+ |
|
|
+ |
|
. |
(3.49) |
∂x |
|
∂y |
∂z |
|||||
Уравнения (3.48) носят название закон трения Стокса. |
|
|||||||
Для несжимаемой жидкости |
divW = 0 , поэтому выражения |
(3.48) |
||||||
упрощаются.
3.2.2. Дифференциальное уравнение неразрывности
Выделим в движущемся потоке элементарный объем в форме прямоугольного параллелепипеда (рис. 3.6) и определим изменение массы dm жидкости в выделенном объеме за элементарный промежуток времени dτ. Это изменение массы определяется разностью между втекающей и вытекающей массой жидкости через грани элементарного объема. Поскольку объем выделенного элемента V = dxdydz остается неизменным с течением времени, то
изменение массы жидкости dm может быть обусловлено лишь изменением ее
плотности dρ |
|
d(ρV ) |
|
|
|
|
dm |
|
dρ |
dρ |
|
|
|
dτ |
= |
dτ |
=V dτ |
= dτ |
dxdydz . |
(3.50) |
32
Рис. 3.6. Выделение элементарного объема в движущемся потоке
Определим массу жидкости, втекающую в выделенный объем за единицу времени.
Жидкость втекает через грани ABFE, AEHD и EFGH в следующих количествах: через грань ABFE − ρwx dydz ; через грань AEHD − ρwy dxdz ; через
грань EFGH − ρwz dxdy . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
Жидкость |
вытекает |
через |
грани DCGF, |
|
BCGF и ABCD в следующих |
|||||||||||||||||||||||||||||||
количествах: через грань DCGH − |
|
ρwx |
+ |
∂(ρw |
x |
) |
|
|
|
|
|
через грань BCGF − |
|||||||||||||||||||||||||
|
|
∂x |
|
dx dydz ; |
|||||||||||||||||||||||||||||||||
|
|
∂(ρwy ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
) |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂(ρw |
z |
|
|
||||
|
ρwy + |
|
|
dy dxdz |
; через грань ABCD − |
ρwz |
+ |
|
|
|
dz dxdy . |
|
|||||||||||||||||||||||||
∂y |
|
∂z |
|
|
|
|
|||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Просуммировав количества втекающей и вытекающей жидкости по всем |
||||||||||||||||||||||||||||||||||||
граням, найдем изменение массы жидкости в выделенном объеме: |
|
||||||||||||||||||||||||||||||||||||
|
|
dρ |
|
|
|
|
∂(ρw |
x |
) |
|
|
∂(ρwy ) |
|
∂(ρw |
z |
) |
|
|
|
|
|
|
|||||||||||||||
|
|
|
dxdydz = − |
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
dxdydz ; |
|
(3.51) |
|||||||||||
|
|
dτ |
∂x |
|
|
|
|
|
∂y |
|
|
|
∂z |
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
dρ |
∂(ρw |
) |
|
|
|
|
∂(ρwy ) |
|
∂(ρw |
) |
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
dτ + |
|
∂x x |
|
+ |
|
|
|
|
|
+ |
|
|
∂z z |
|
= 0 . |
|
|
|
|
(3.52) |
||||||||||||
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
Уравнение |
(3.52) |
|
называется |
|
|
|
дифференциальным |
уравнением |
||||||||||||||||||||||||||||
неразрывности или сплошности. Оно также может быть записано в виде |
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
dρ |
+ div(ρWr)= 0. |
|
|
|
|
|
|
|
|
|
|
|
(3.53) |
||||||||||||||||
|
|
|
|
|
|
|
|
dτ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Для потоков несжимаемой жидкости (стационарных и нестационарных) |
||||||||||||||||||||||||||||||||||||
уравнение неразрывности примет вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
∂wx |
+ |
|
|
∂wy |
+ |
∂wz |
|
= 0 ; |
|
|
|
|
|
|
|
|
|
|
(3.54) |
|||||||||||
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
dy |
|
|
|
dz |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
divW = 0 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.55) |
|||||||||||
33
Для стационарных потоков газа (сжимаемой жидкости) уравнение
неразрывности примет вид |
|
|
∂(ρwy ) |
|
|
|
|
|
|
||
∂(ρw |
) |
|
|
|
∂(ρw |
) |
|
|
|
||
∂x x |
|
+ |
|
|
+ |
∂z z |
|
= 0 |
; |
(3.56) |
|
|
|
∂y |
|
|
|||||||
|
|
|
div(ρW )= 0 . |
|
|
|
(3.57) |
||||
Необходимо отметить, что существуют и другие способы вывода уравнения неразрывности.
3.2.3.Дифференциальные уравнения переноса количества движения. Уравнения Эйлера и Навье-Стокса
Рассмотрим движение сплошной среды, предполагая скорость, плотность, напряжения и массовые силы непрерывными функциями времени и координат. В декартовой системе координат выделим элемент в виде прямоугольного параллелепипеда (рис. 3.6). К выделенному объему применяется закон сохранения количества движения, в соответствии с которым, изменение за определенный промежуток времени количества движения жидкости в элементарном объеме, равно импульсу внешних сил, действующих на этот объем.
Изменение количества движения в выделенном объеме происходит за счет изменения плотности жидкости и ее скорости, так и за счет разницы между втекающим и вытекающим количеством движения через границы этого объема.
Уравнение движения удобно рассматривать в проекциях на координатные оси x, y, z. Для получения проекций уравнения движения на каждую координатную ось применим закон сохранения количества движения к граням выделенного объема, перпендикулярным соответствующей оси.
Рассмотрим сначала грань ABFE, перпендикулярную оси x, и определим
составляющие потока количества движения через |
нее за время dτ |
в |
направлении координатных осей x, y, z: |
|
|
I xx (x)= ρwx wx dydzdτ =ixxdydzdτ ; |
|
|
I xy (x)= ρwx wy dydzdτ = ixy dydzdτ ; |
(3.58) |
|
I xz (x)= ρwx wz dydzdτ = ixz dydzdτ .
где ixx , ixy , ixz − составляющие потока импульса через грань, имеющую
площадь, равную 1 и перпендикулярную оси x, за единицу времени в направлении координатных осей x, y, z соответственно.
34
Потоки количества движения через грань DCGH можно выразить
зависимостями: |
|
∂(ρw |
w |
|
) |
|
|
∂i |
|
|
|
||
I xx (x + dx) |
|
ρwx wx + |
|
|
|
|
|
|
|||||
= |
x |
|
x |
|
|
dx dydzdτ = ixx + |
|
xx |
dx dydzdτ ; |
|
|||
|
|
|
∂x |
|
|
|
|
|
|
∂x |
|
|
|
|
|
|
∂(ρwx wy ) |
|
|
∂ixy |
|
|
|||||
I xy (x + dx)= ρwx wy + |
|
|
|
|
|
|
|
|
|
|
(3.59) |
||
|
|
|
|
|
|
|
|
||||||
∂x |
|
|
|
|
dx dydzdτ = ixy + |
∂x |
dx dydzdτ ; |
||||||
|
|
|
|
|
|
|
|
|
|
|
|||
I xz (x + dx) |
|
ρwx wz + |
∂(ρw |
w |
|
) |
|
|
∂i |
|
|
|
|
= |
x |
|
z |
|
|
dx dydzdτ = ixz + |
|
xz |
dx dydzdτ . |
|
|||
|
|
|
∂x |
|
|
|
|
|
|
∂x |
|
|
|
Разница |
|
между |
втекающим |
(положительным) и вытекающим |
|||||||||
(отрицательным) количеством движения через грани, перпендикулярные оси x, в направлениях осей x, y, z определяется как
∆I xx = Ixx (x)− Ixx (x + dx)= − |
∂(ρwx wx )dxdydzdτ ; |
|
|
|
∂x |
|
|
∆I xy = I xy (x)− I xy (x + dx)= − |
∂(ρwx wy ) |
dxdydzdτ ; |
(3.60) |
|
|||
|
∂x |
|
|
∆I xz = Ixz (x)− I xz (x + dx)= − |
∂(ρwx wz )dxdydzdτ . |
|
|
|
∂x |
|
|
Аналогично определяются изменения составляющих количества движения через грани перпендикулярные осям y и z:
∆I yx = I yx (y)− I yx (y + dy)= −∂(ρwy wx )dxdydzdτ ;
∂y
∆I yy = I yy (y)− I yy (y + dy)= −∂(ρwy wy )dxdydzdτ ;
∂y
∆I yz = I yz (y)− I yz (y + dy)= −∂(ρwy wz )dxdydzdτ .
∂y
∆I zx = I zx (z)− I zx (z + dz)= −∂(ρwz wx )dxdydzdτ ;
∂z
∆I zy = I zy (z)− I zy (z + dz)= −∂(ρwz wy )dxdydzdτ ;
∂z
∆I zz = I zz (z)− I zz (z + dz)= −∂(ρwz wz )dxdydzdτ .
∂z
(3.61)
(3.62)
35
Изменение за время dτ количеств движения в выделенном объеме ∆I x , ∆I y , ∆I z в направлении осей x, y, z, происходящее за счет изменения плотности жидкости и ее скорости, выразим соотношениями вида
∆I x = ∂(∂ρτwx )dxdydzdτ ;
∆I y = |
∂(ρwy ) |
dxdydzdτ ; |
(3.63) |
||
∂τ |
|
||||
|
|
|
|||
∆Iz = |
∂(ρwz )dxdydzdτ . |
|
|||
|
∂τ |
|
|
||
На выделенный объем, в общем случае, действуют внешние массовые (гравитационные или инерционные) и поверхностные силы. Найдем проекции импульса массовых сил на оси x, y, z за время dτ .
Для внешней массовой силы R = i X + jY + kZ , отнесенной к единице объема, имеем
IRx = Xdxdydzdτ ; |
IRy =Ydxdydzdτ ; |
IRz = Zdxdydzdτ . |
(3.64) |
|||||||
Проанализируем поверхностные силы. На грань ABFE, перпендикулярную |
||||||||||
оси x, действует внешняя поверхностная сила |
Px (x), а на противоположную |
|||||||||
грань DCGH (и в противоположном направлении) – сила Px (x + dx). |
|
|||||||||
r |
r |
|
|
r |
r |
|
|
∂pr |
|
|
P |
(x)= p |
|
dydz ; |
P |
(x + dx)= p |
|
+ |
x |
dx dydz . |
(3.65) |
|
|
∂x |
||||||||
x |
|
x |
|
x |
|
x |
|
|
||
Вуравнении (3.65) prx − вектор результирующего напряжения,
действующего на площадку, перпендикулярную оси x.
Вектор импульса результирующей поверхностной силы,rдействующей на перпендикулярные оси x грани выделенного объема, I Px определяется
выражением |
∂px |
|
|
|
IrPx = − |
dxdydzdτ . |
(3.66) |
||
|
||||
|
∂x |
|||
Заметим, что направление результирующей поверхностной силы, действующей на перпендикулярные оси x грани, так же, как и направление ее импульса, не совпадает с осью x.
Аналогично получаются выражения для импульсов результирующей поверхностной силы, действующей на грани, перпендикулярные осям y и z:
r |
|
∂py |
|
|
||
I Py |
= − |
|
|
dxdydzdτ . |
(3.67) |
|
∂y |
||||||
|
|
|
||||
IrPz |
= − |
∂pz |
|
dxdydzdτ . |
(3.68) |
|
∂z |
||||||
|
|
|
||||
36
Векторы напряжений px , py , pz , действующих на перпендикулярные осям x, y, z соответственно грани, можно выразить через составляющие тензора
напряжений |
r |
|
|
r |
σ x + jτxy + kτxz ; |
|
|
px = i |
|
||
r |
r |
τxy + jσ y + kτ yz ; |
|
py |
= i |
(3.69) |
|
r |
r |
τxz + jτ yz + kσ z . |
|
pz |
= i |
|
|
Суммарный импульс поверхностных сил
IrP = − irr ∂∂σxxk ∂τxz∂x
|
∂τxy |
|
|
∂τ |
|
|
|
r ∂τxy |
|||
+ |
|
|
|
+ |
|
xz |
+ |
|
|
||
|
|
|
|
|
|||||||
|
∂y |
∂z |
|
j |
∂x |
||||||
|
|
|
|
|
|
|
|||||
|
|
∂τ |
yz |
|
∂σ |
|
|
|
|
||
+ |
|
|
+ |
|
|
z |
|
|
|
||
|
|
|
|
|
|
|
|
||||
|
|
∂y |
|
|
∂z |
|
|
|
|
||
|
|
|
|
|
|
|
|
||||
|
∂σ |
y |
|
∂τ |
yz |
|
|
|
|
|
+ |
|
+ |
|
|
+ |
|
|
|
||
|
|
|
|
|
|
|||||
|
∂y |
∂z |
|
|
|
|
||||
|
|
|
|
× dxdydzdτ . |
(3.70) |
|||||
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
Таким образом, для различных составляющих закона сохранения
количества движения, получаем |
|
|
|
|
|
|
|
|
∂τxy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
− (∆I xx + ∆I yx + ∆I zx )= I Rx |
∂σ |
|
x |
|
|
|
|
|
|
|
|
|
∂τ |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
xz |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
∆I x |
− |
|
|
|
|
+ |
|
|
|
|
|
+ |
|
|
|
|
|
|
dxdydzdτ ; |
|
|
|
|||||||||||||||||||||||||||||
∂x |
|
|
∂y |
|
|
|
∂z |
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
− (∆I xy + ∆I yy + ∆I zy )= I Ry |
|
∂τ |
xy |
|
|
|
|
|
∂σ |
y |
|
|
|
|
|
|
∂τ |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
yz |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
∆I y |
− |
|
|
|
|
+ |
|
|
|
|
+ |
|
|
|
|
|
|
dxdydzdτ ; |
|
|
|
||||||||||||||||||||||||||||||
∂x |
|
|
∂y |
|
|
∂z |
|
|
|
(3.71) |
||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
− (∆I xz + ∆I yz + ∆I zz )= I Rz |
|
∂τ |
xz |
|
|
|
|
|
∂τ yz |
|
|
|
|
|
|
∂σ |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
∆I z |
− |
|
|
|
|
+ |
|
|
|
|
|
+ |
|
|
|
|
|
|
dxdydzdτ . |
|
|
|
|
||||||||||||||||||||||||||||
∂x |
|
|
∂y |
|
|
|
∂z |
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
Преобразуя эти уравнения получаем |
|
|
уравнения движения в напряжениях: |
||||||||||||||||||||||||||||||||||||||||||||||||
∂(ρw |
|
|
) |
|
|
∂(ρw |
w |
|
|
) |
|
∂(ρwx wy ) |
|
∂(ρw |
|
w |
|
|
) |
|
|
|
|
|
|
|
|
∂σ |
|
|
|
|
|
∂τxy |
|
∂τ |
|
|
|
|
||||||||||||
|
|
x |
|
|
+ |
|
x |
|
x |
|
|
+ |
|
|
|
|
+ |
|
x |
|
|
z |
|
|
|
= X |
+ |
|
|
|
|
x |
+ |
|
|
|
+ |
|
xz ; |
|
||||||||||||
|
∂τ |
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|||||||||||||||||||||||||||||
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
∂z |
|
|
||||||||||
|
∂(ρwy ) |
+ |
|
∂(ρwx wy ) |
+ |
∂(ρwy wy ) |
+ |
∂(ρwy wz ) |
|
=Y + |
∂τxy |
|
+ |
∂σ y |
+ |
∂τ yz |
; |
(3.72) |
||||||||||||||||||||||||||||||||||
|
∂τ |
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
∂y |
|
|
∂z |
|
|
|
|
|
|
|
|
|
∂x |
|
|
|
∂y |
|
∂z |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
∂(ρw |
|
) |
|
|
∂(ρw |
|
w |
|
) |
|
|
∂(ρwy wz ) |
|
∂(ρw |
|
w |
|
) |
|
|
|
|
|
|
|
∂τ |
|
|
|
|
|
|
∂τ yz |
|
∂σ |
|
|
|
|
||||||||||||
|
|
|
z |
|
|
+ |
|
x |
|
|
z |
|
|
+ |
|
|
|
|
+ |
|
z |
|
|
z |
|
|
|
= Z |
+ |
|
|
|
xz |
|
+ |
|
|
|
+ |
|
|
z . |
|
|||||||||
|
∂τ |
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
∂y |
∂z |
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
|||||||||||
37
|
Левая часть уравнений (3.72) может быть преобразована как |
|
|
||||||||||||||||||||||||||||||||||||||||||
∂(ρw |
|
) |
|
|
∂(ρw |
w |
|
) |
|
|
∂(ρwx wy ) |
|
|
∂(ρw |
w |
) |
|
|
∂w |
|
|
∂ρ |
|
∂w |
|
|
|||||||||||||||||||
|
x |
|
|
+ |
|
|
|
|
x |
|
|
x |
|
+ |
|
|
|
|
|
|
|
|
|
+ |
|
|
|
x |
z |
|
= ρ |
|
x + wx |
|
|
|
+ ρwx |
|
x + |
|
|||||
∂τ |
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
∂τ |
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
∂(ρwy ) |
|
∂z |
|
|
|
|
∂τ |
|
|
∂x |
|
||||||||||||||||
|
∂(ρw |
) |
|
|
|
|
|
|
∂w |
|
|
|
|
|
|
|
|
|
∂w |
|
|
|
∂(ρw |
) |
|
|
|
|
|
||||||||||||||||
+ wx |
|
|
|
|
|
x |
|
+ ρwy |
|
|
|
x |
+ wx |
|
|
|
|
|
|
|
|
+ ρwz |
|
|
x + wx |
|
z |
|
= |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
∂z |
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
. |
(3.73) |
|||||||||||
|
∂w |
|
|
|
|
|
∂w |
|
|
|
|
|
|
|
∂w |
|
|
|
|
|
|
|
|
∂w |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
x |
+ w |
|
x |
+ w |
|
x |
|
+ w |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
= ρ |
|
|
|
|
x |
|
|
y |
|
|
z |
|
|
|
x + |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
∂τ |
|
|
|
∂x |
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
∂ρ |
|
∂(ρw |
|
) |
|
|
∂(ρwy ) |
|
|
∂(ρw |
|
) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
+ wx |
|
|
|
|
+ |
|
|
|
|
x |
|
|
+ |
|
|
|
|
|
|
+ |
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
∂τ |
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Аналогично преобразовываются и два других уравнения.
В соответствии с уравнением неразрывности (3.52) последняя сумма в выражении (3.73), стоящая в квадратных скобках, равна нулю. Обозначим
Dwx |
= |
∂wx |
+ w |
|
∂wx |
+ w |
|
∂wx |
+ w |
|
∂wx |
, |
(3.74) |
|
|
|
|
|
|||||||||
dτ |
∂τ |
x ∂x |
y ∂y |
z ∂z |
|||||||||
где Dwdτx − субстанциальная производная.
Таким образом, получаем уравнения движения в форме уравнений Навье-
Стокса:
|
|
Dw |
x |
|
|
|
|
|
|
|
|
|
∂p ∂ |
|
|
∂w |
|
|
|
|
2 |
|
|
|
r |
|
|
|
∂ |
|
|
|
∂w |
|
|
|
∂wy |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
= X − |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
x |
|
|
|
divW |
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
||||||||||||||
ρ |
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
µ |
|
|
|
− |
|
+ |
|
|
|
|
|
|
|
|
+ |
|
|
|
+ |
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
µ |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
dτ |
|
|
|
|
|
∂w |
|
|
|
∂x |
|
|
|
|
∂x |
|
|
|
|
∂x |
|
|
|
3 |
|
|
|
|
|
|
|
∂y |
|
|
∂y |
|
|
|
|
dx |
; |
|
||||||||||||||||||
|
|
∂ |
|
|
|
|
|
x |
|
|
∂w |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
+ |
|
|
|
µ |
|
|
|
|
+ |
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
∂z |
|
|
|
|
|
|
|
dx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
Dwy |
|
|
|
|
|
|
|
|
|
∂p ∂ |
|
∂w |
|
|
|
|
∂wy |
|
∂ |
|
|
|
|
|
|
∂wy |
|
|
|
|
2 |
|
|
|
r |
|
|
|||||||||||||||||||||||
ρ |
|
|
|
|
|
|
|
=Y − |
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|
x |
+ |
|
|
|
|
|
+ |
|
|
|
|
|
2 |
|
|
|
|
− |
|
|
divW |
|
+ |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
dτ |
|
|
|
∂y |
|
∂x |
µ |
∂y |
|
dx |
|
|
|
|
µ |
|
∂y |
|
|
|
3 |
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
; |
(3.75) |
||||||||||||||||||||||
|
|
∂ |
|
|
|
|
∂wy |
|
|
|
∂w |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
+ |
|
|
|
|
µ |
|
|
|
|
|
|
+ |
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
dy |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
Dw |
z |
|
|
|
|
|
|
|
|
∂p ∂ |
|
∂w |
|
|
|
∂w |
|
|
|
|
∂ |
|
∂wy |
|
|
|
∂w |
|
|
|
|
|
||||||||||||||||||||||||||
|
|
ρ |
|
|
|
|
|
|
|
= Z |
|
|
|
|
|
|
|
|
|
|
|
µ |
|
|
|
x |
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z |
+ |
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
|
+ |
|
|
|
|
|
|
+ |
|
|
|
+ |
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
dτ |
|
|
|
∂z |
|
|
|
|
|
∂z |
|
dx |
|
|
|
µ |
∂z |
|
|
dy |
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
. |
|
|||||||||||||||||||||||||||
|
|
|
|
∂ |
|
|
|
|
|
|
|
∂w |
|
|
|
|
2 |
|
|
|
r |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
+ |
|
|
|
|
|
µ |
2 |
|
|
|
|
|
− |
|
|
divW |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
∂z |
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
38
Система уравнений (3.75) описывает движение вязкой сжимаемой
жидкости (газа). Для несжимаемой жидкости divW = 0 , поэтому уравнения Навье-Стокса упрощаются
|
|
|
|
|
Dwx |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|||||||||
|
|
ρ |
|
|
= X − ∂p |
+ µ |
∂ |
wx |
+ ∂ |
wx + ∂ |
|
wx |
; |
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
dτ |
|
|
|
|
|
|
∂x |
|
|
|
|
∂x |
2 |
|
|
|
∂y |
2 |
|
|
|
|
|
|
∂z |
2 |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
Dwy |
|
|
|
|
|
|
∂p |
|
|
|
∂2 wy |
|
∂2 wy |
|
|
|
|
|
∂2 wy |
|
|
|
|
|||||||||||||||||||||||
|
|
ρ |
|
|
|
|
|
|
=Y − |
|
|
+ µ |
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
|
; |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
dτ |
|
|
|
|
|
|
∂y |
|
|
|
|
∂x |
|
|
|
∂y |
|
|
|
|
|
|
∂z |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
Dwz |
|
|
|
|
|
|
∂p |
|
|
|
2 |
wz |
|
|
2 |
wz |
|
|
|
|
|
|
2 |
wz |
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
= Z − |
|
|
|
∂ |
|
∂ |
|
|
|
|
|
∂ |
|
|
|
|
|||||||||||||||||||||||||||||
|
|
ρ |
|
|
|
|
|
|
|
+ µ |
|
|
|
|
2 |
|
+ |
|
|
|
|
|
2 |
+ |
|
|
|
|
|
2 |
|
. |
|
|
|||||||||||||||||
|
|
|
|
dτ |
∂z |
|
∂x |
|
|
∂y |
|
∂z |
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
Для идеальной среды µ = 0 из уравнений (3.76) получим |
|||||||||||||||||||||||||||||||||||||||||||||||||||
ρ |
Dw |
x |
|
|
|
= X |
− |
∂p |
; ρ |
Dwy |
|
=Y − |
∂p |
; |
|
|
ρ |
Dw |
z |
= Z − |
∂p |
. |
|||||||||||||||||||||||||||||
|
dτ |
|
|
|
|
|
∂x |
|
dτ |
∂y |
|
|
|
dτ |
|
∂z |
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
или |
|
∂wx |
|
|
|
|
|
|
|
∂wx |
|
|
|
|
|
|
∂wx |
|
|
|
|
|
|
∂wx |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
1 ∂p |
|
|
||||||||||
|
|
|
|
|
+ wx |
|
|
+ wy |
|
|
+ wz |
|
|
= |
|
|
X |
− |
|
; |
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
∂x |
|
∂y |
∂z |
|
ρ |
|
|
ρ ∂x |
|
||||||||||||||||||||||||||||||||||||||
|
|
∂τ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
∂wy |
|
|
|
+ wx |
|
∂wy |
|
+ wy |
|
∂wy |
|
+ wz |
∂wy |
= |
1 |
|
Y |
− |
|
1 |
|
∂p |
; |
|
|||||||||||||||||||||||||
|
|
∂τ |
|
|
|
|
|
∂x |
|
|
∂y |
|
∂z |
ρ |
|
|
ρ ∂y |
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
∂wz |
+ wx |
∂wz |
|
+ wy |
|
∂wz |
|
+ wz |
∂wz |
|
= |
1 |
|
|
Z |
− |
|
|
1 |
∂p . |
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
ρ |
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
∂τ |
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
|
∂y |
|
|
|
|
|
|
∂z |
|
|
|
|
|
|
|
|
ρ ∂z |
|
|
||||||||||||||||
(3.76)
(3.77)
(3.78)
Эти уравнения называются уравнениями Эйлера и описывают движение идеальной сплошной среды.
3.2.4. Дифференциальное уравнение энергии
Рассмотрим движение сплошной среды, полагая скорость, плотность, напряжения и массовые силы непрерывными функциями времени и координат. В декартовой системе координат выделим элемент в виде прямоугольного параллелепипеда (рис. 3.6). К выделенному объему применяется закон сохранения полной энергии E, включающий внутреннюю u и кинетическую
W 2
2
определенный промежуток времени количества энергии жидкости в элементарном объеме равно работе L внешних сил, действующих на этот объем, подводимому (отводимому) количеству теплоты Qτ и выделившемуся количеству теплоты при наличии внутренних источников тепла QVτ :
ρ |
DE dxdydz = L + Q |
+ Q |
. |
|
|
|
dτ |
τ |
V τ |
|
(3.79) |
|
|
|
|
||
39 |
|
Работу внешних массовых сил выразим формулой вида |
|
Lм = (Xwx +Ywy + Zwz )dxdydzdτ . |
(3.80) |
Определим работу нормальных и касательных сил, действующих на грани выделенного объема. Действие на грань ABFE нормальных напряжений
обуславливает работу |
(x)= −σx wx dydzdτ . |
|
|
||||||
|
Lσx |
|
(3.81) |
||||||
На грань противоположную грань DCGH: |
|
||||||||
L |
(x + dx)= |
σ |
w |
x |
+ |
∂(σx wx ) |
dx dydzdτ . |
|
|
|
|
||||||||
σx |
|
x |
|
|
∂x |
|
(3.82) |
||
|
|
|
|
|
|
|
|||
Работа сжимающих напряжений считается отрицательной, а растягивающих – положительной.
Результирующая работа нормальных напряжений, действующих на грани,
перпендикулярные оси x: |
|
∂(σx wx ) |
|
|
||
L |
= L |
(x)+ L |
(x + dx)= |
dxdydzdτ . |
|
|
|
|
|||||
σx |
σx |
σx |
|
∂x |
(3.83) |
|
|
|
|
|
|||
Аналогично определяется результирующая работа касательных напряжений, действующих на грани, перпендикулярные оси x:
L |
= |
|
∂(τxy wy ) |
dxdydzdτ ; |
(3.84) |
|
|
|
|||||
τxy |
|
|
∂x |
|||
|
|
|
|
|
|
|
L |
= |
|
∂(τxz wz ) |
dxdydzdτ . |
(3.85) |
|
|
|
|||||
τxz |
|
|
∂x |
|||
|
|
|
|
|||
Аналогичные рассуждения проводятся применительно к граням, перпендикулярным осям y и z.
Таким образом, работа поверхностных сил будет определяться выражением вида
|
∂ |
(σ |
|
w |
|
+ |
τ |
|
w |
|
+τ |
|
w |
|
)+ |
|
|
∂ |
|
(τ |
|
w |
|
+σ |
|
w |
|
+τ |
|
w |
|
)+ |
|
||||
|
|
∂x |
x |
x |
xy |
y |
xz |
z |
|
|
|
|
xy |
x |
y |
y |
yz |
z |
|
||||||||||||||||||
L |
|
|
|
|
|
|
|
|
|
|
∂y |
|
|
|
|
|
|
|
|
||||||||||||||||||
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dxdydzdτ . |
(3.86) |
|
п |
|
|
∂ |
(τxz wx |
+τ yz wy +σz wz ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Работа L внешних сил определяется как |
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L = Lм + Lп . |
|
|
|
|
|
|
|
|
|
(3.87) |
||||||||||
|
Количество теплоты, подведенное через к выделенному объему, |
||||||||||||||||||||||||||||||||||||
перпендикулярно оси x, через грань ABFE: |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Qτx (x) |
= qx dydzdτ , |
|
|
|
|
|
|
|
|
(3.88) |
||||||||||||
где qx − плотность теплового потока, направленная вдоль оси x, Вт/м2. |
|
||||||||||||||||||||||||||||||||||||
|
Количество теплоты, отведенное через грань DCGH: |
|
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂q |
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Q |
(x + dx)= q |
|
+ |
|
|
dx dydzdτ . |
|
|
|
|
(3.89) |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
τx |
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
||||
