Механика сплошной среды
..pdfМенее общим* является предположение, что u(t), ф(£), s(t) —
функции («внутренних», |
«структурных», |
«микроскопических»...) |
|||
«параметров |
состояния» |
mk = \ik(t), |
подчиняющихся эволюцион |
||
ным уравнениям |д = со(р, П (0), |
|
|
|
||
|
——k.^ —= ak (И- (Q. П (0) |
(*=0, 1, .. •), |
(10.25") |
||
|
at |
|
|
|
|
где со(|д, П) |
— выбираемые из различных |
физических, |
нередко |
интуитивных и качественных, соображений функции. Импонирует гипотеза, что скорость роста объемной плотности величины p,(f, х)
(трещин — в твердом теле или воздушных |
пузырьков — в воде |
и т. п.) является какой-то функцией самой |
p(f) и параметров |
процесса 11(f) или реакции R (t) (10.5). |
|
Соотношения (10.25") нередко дополняются гипотезой о функ ции рассеяния и работе
а»*(0=а>к(М-(0, И (0 )М 9 .
-±-b'A=Phd\ih (k=0, |
1 ,...,л ), |
(10.26) |
||
Р |
|
|
|
|
после чего из (10.14) и (10.20) |
получается уравнение |
|
||
( ~ r - + T - f -----u>k-Pk)d\ih= 0, |
(10.27') |
|||
\ д^к |
дцк |
|
} |
|
решение которого предполагается в виде |
|
|||
|
рк = - г - + т- г -— |
“V |
(Ю-27") |
|
|
д[1к |
d\ik |
|
|
Параметры типа р могут относиться к числу параметров Р(т), |
||||
Рл в (10.5). |
(10.25), |
(10.25'), |
(10.25") являются |
пара |
Все функционалы |
метрами термомеханической природы, и соответствующие им оп ределяющие уравнения приводимы к виду (10.5).
В случае (10.24) |
уравнения |
состояния |
имеют вид согласно |
|||
(10.14) |
|
|
|
|
|
|
а |
|
&i = |
л |
дг|э |
(10.27) |
|
йГ ’ |
~д¥' |
|||||
|
^ |
7 ’ |
|
Но без экспериментального изучения конкретной сплошной среды, ее структуры и свойств, которые и могут быть отражены параметрами типа (10.25'), (10.25"), нельзя указать прямых опы тов, которые позволили бы установить число и выражение пара
метров состояния через характеристики процесса |
t), T(t). |
В частности, в равновесных опытах с постоянными |
по объему и |
меняющимися во времени параметрами состояния T(t)t р«(0> б'С?г выражение (10.2) на основании (10.22) преобразуется к виду (при ш* = 0)
т. е. б' Q T имеет интегрирующий множитель Г-1:
^2
(10.28)
Этот интеграл зависит не от того, по какому закону во времени: подается тепло 6<3г и изменяются параметры рл, а только от ко нечных и начальных значений параметров. Равенство (10.28) и служит для нахождения энтропии s(T, р). Аналогичное равенство* вытекающее для рассматриваемого объема из (10.12), служит для определения внутренней энергии а(Т, р). Но для нахождения числа и физического смысла параметров р в макроопытах могут быть полезны лишь функционалы (10.6), так как параметры в силу основного постулата должны преобразовывать функциона лы процесса (10.25) в функции этих параметров.
Для необратимых процессов основное термодинамическое тож
дество (10.14) |
может быть преобразовано к общему виду, |
если |
|||||
обобщить определение процесса я(т), реакции r(t) |
и скалярного |
||||||
потенциала V* |
которые до этого принимались я = П, |
r= R , |
I |
= |
|||
Тождественные преобразования дают |
|
|
|
|
|
||
dw = du—Tdr\—ad§ = dij?,, + r\dT—ad? = |
|
|
|
||||
= dy— Tdr\ + S i = d%+ т]dT + d a i= 0 , |
|
(10.29') |
|||||
где |
|
|
|
|
|
|
|
|
% = и ~ - Т т ] , ф = u — aS, |
|
|
(10.29) |
|||
|
% = u - T i \ — o § = % - - e i = q > — Tr\. |
|
|||||
|
|
|
|
||||
Отсюда, выбирая различные допустимые пары процессов я и |
|||||||
соответствующих реакций г, получим всегда |
|
|
|
|
|||
|
w‘(it (t)) = |
V' (it)—r$) it (0 = |
0 |
|
(10.30) |
||
для произвольного л (т), |
При |
этом |
процессы я, |
реакции |
|||
г* и скалярный функционал |
V1 определяются таблицей 5 функций |
||||||
(я, г, V). |
г можно рассматривать |
|
|
|
|
|
|
Объекты л, |
как |
(4X4) |
квадратные |
||||
матрицы или |
семимерные |
(7-мерные) |
ортогональные |
векторы |
|
|
|
|
|
|
|
|
|
|
|
|
Т а б л и ц а 5 |
||
|
|
я |
|
|
|
Г |
|
V' |
|
|
V (t)= v* |
|
||
1 |
(Л. 3) |
|
|
(Т,о) |
|
и( |
|
|
|
и |
|
|||
2 |
(л. о) |
|
(Т, |
—9) |
|
¥ |
|
|
Ф= ц — а^Г |
|||||
3 |
(Г. |
9 ) |
|
|
(—л. <?) |
|
< |
|
|
*£- II |
1 |
.3 |
||
4 |
|
(Т, |
а) |
|
(—Л> - 9 ) |
|
%* |
|
|
X = tri — а?' |
||||
5 |
(и, |
1) |
|
(Г-1, —Г-%) |
|
V |
|
|
|
Л |
|
|||
|
|
|
|
|
о |
|
|
|
|
л0 |
о |
|
|
(10.31) |
|
|
|
|
Я = ( Я 0, Я0, JTJL, |
>ЛБ) = |
|
|
|||||||
|
|
|
|
О ||я"»"|| |
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
Г = |
г0 |
О |
|
= |
о |
|
,Гб) |
( т ,/1=1, 2,3) |
|
||||
|
|
(r, г0, гъ |
|
|
||||||||||
|
|
|
О |
||rmn|| |
|
|
|
|
|
|
|
|
|
|
со скалярным умножением |
|
|
|
|
|
|
|
|||||||
|
гя = |
|
0 |
0 |
|
. |
|
. |
0 |
0 |
|
|
(10.32) |
|
|
г я = г 0(я0)-+ г0я0+ ... + гБяБ= г 0(я0)‘+ гтпятГ1 |
|||||||||||||
причем |
о |
о |
|
|
|
элементы |
(скаляры) |
в |
столбцах таблицы; |
|||||
я0 г0— первые |
||||||||||||||
Лтпу ттп — компоненты тензоров, |
указанных |
в столбцах я, г; |
||||||||||||
(я0 яь |
|
Я5) |
и (г0 |
гь |
г5) |
— построенные согласно § 9 по |
||||||||
тензорам я и 7 шестимерные ортогональные векторы |
(деформаций |
|||||||||||||
и напряжений, |
и а, или наоборот — согласно таблице). |
|
||||||||||||
Например, в третьей строке таблицы: |
|
|
|
|
|
|||||||||
|
|
|
л(х) = |
(Г, |), г—г1—(— л> о); |
|
|
|
|
||||||
|
|
|
п0= Т , |
г0= — г); птп= втп, гтп= а тп) |
|
(10.33) |
||||||||
|
|
|
|
о о |
|
. |
ГтпПтп= ОтпЕтп. |
|
|
|||||
|
|
|
|
Го Ы ' = — Л Т, |
|
|
||||||||
Уравнение |
(10.30) имеет вид |
|
|
|
|
|
|
|
||||||
|
|
|
|
w* = |
'Фп (л (т) ) ~ |
г< (л) я (0 = 0 ; |
|
|
Г‘ (я)я = — л' (Я) 7 (0 + <У?п(Я) ётп (t).
Если ввести в него вместо ц рассеяние |
ад* |
и |
энтропию s |
|||||
(10.20) и функцию ф (10.23), положив |
|
|
|
|
|
|
||
г^ф ^ + Гт^ф -Е Ts, Ti\= Ts— ад* |
|
|
|
(10.35) |
||||
сохранив процесс я = (Г, <§) и введя |
реакцию rs= (—s, |
а), |
полу |
|||||
чим основное термодинамическое тождество |
(10.23) |
в виде |
|
|||||
ф* (я)—г*(я)я(/) + ад* = |
0; я = |
(71 £), |
гs = |
(— s, а). |
(10.36) |
|||
Все параметры, |
входящие |
в уравнение |
(10.30) |
и в таблицу |
||||
(я, г, V), могут быть измерены в J^-опытах; предполагается, что |
||||||||
именно по данным этих ^-опытов и построен функционал |
У'(я), |
|||||||
вследствие чего (10.30) и называется |
тождеством. Здесь |
полная |
||||||
аналогия с законом |
Ньютона |
/пг = Р (г(т)). «Образцы |
сплошной |
среды» нередко практически недоступны для .^-опытов, но неко торые закономерные проявления ее свойств в некоторых натурных условиях имеются. Тогда можно принять некоторую гипотезу о виде одного скалярного функционала У*(я), например энергии и*(я), и если из тождества (10.30) найти уравнения состояния, т. е. реакцию /*'(я), то решить соответствующую натурному явле нию задачу МСС; сравнение результатов с известными законо мерными проявлениями позволит оценить правильность гипотезы.
Этот путь научного познания древнее наших наук. |
среды, |
||||
Определение реакции |
г (я), т. е. уравнений состояния |
||||
из тождества |
(10.30) по |
одному |
заданному |
функционалу |
0 (я ) , |
следовательно, |
представляется |
необходимым, |
как и выбор |
этого |
функционала с учетом свойств довольно большого числа уже из вестных в МСС.
Вид тождества (10.36) может быть придан тождеству (10.30) применительно к любой строке таблицы (я, г, V); он отличается выделенным слагаемым ад*^50, имеющим много названий: рассе яние, диссипация, некомпенсированное тепло и т. д. Как уже ска
зано, все параметры, входящие в |
таблицу (я, |
г, |
V) и |
(10.30), |
в принципе могут быть измерены |
в J t-опытах, |
но |
этого |
нельзя |
сказать о w* = T(s—г\)^0 и об энтропии s: измерим функционал
энтропии Т1'(я), но не st(n) и ад* по отдельности. Эта неопреде ленность породила ряд гипотез о существовании вместе с (10.30) еще одного термодинамического уравнения (kd).
Целесообразно следующее доопределение функционалов рассе яния ад* и энтропии s, соответствующее сложившимся в МСС представлениям: если в конечном объеме среды возможен ста
ционарный периодический |
процесс |
П(Г, # )= П (х , |
/)=П (х, / + |
|
+ 2л/ш) с частотой со то |
рассеяние |
ад* |
и энтропия |
5 для произ |
вольного момента t и любого объема |
V с границей 2 удовлетво |
|||
ряют условиям доопределения: |
|
|
|
*4+2л/(й |
d t\p w ’dV = Q1] Q ,= |
/4-{-2я/(0 |
(— qn) d2; |
/4*f2Jl/0 |
dt\pTsdV=Q . |
С |
J dt |
jj |
|||
it |
V |
tx |
2 |
*4 |
V |
|
|
|
|
|
(10.37) |
Здесь Qi — количество тепла, выделяемое за один цикл массой объема V через 2 во внешнюю среду.
При заданном функционале ф*(я) не найдено общего решения уравнения (10.36), т. е. реакции iV(n). Но если ф'(я) имеет
дифференциал Фреше, например
|
|
|
|
|
t |
|
|
|
|
(я + 6я)— ф (я) = 6 ^ = J ^ |
J 6я (т) dr, |
(10.38) |
|||||
|
|
|
|
|
и ' |
т |
|
|
и если ф'(л) |
имеет вид |
|
|
|
|
|
|
|
|
|
У (л)= |
|
|
) Я + |
|
|
(10.39) |
причем для |
любого |
процесса |
функция фиу^О, то |
(10.36) |
имеет |
|||
частное решение, физически допустимое, |
|
|
|
|||||
|
г '(я )= ( |
) |
, |
w = ~ ^ w > 0; |
|
(10.40) |
||
|
|
Dn |
)x=t |
|
|
|
|
|
иначе |
|
|
|
|
|
|
|
|
|
\ |
DT i t |
Oi]= |
— |
\ Den It |
|
(10.41) |
|
|
|
|
P |
|
|
|||
причем необходимо еще выполнение (10.37) для s. |
(10.30) |
возни |
||||||
На основе произвола л(т) |
в тождестве в форме |
|||||||
кает и более строгий подход к определению реакции г'(я). |
|
|||||||
Поскольку тождество (10.30) |
справедливо для любого процесса |
|||||||
я(т), оно верно и для я + бл, |
и возникает вариационный принцип, |
|||||||
на основе которого можно искать решение |
(10.30) при следующих |
предположениях и обозначении со = гя.
1.Функции физически допустимых процессов я(т), реакций
г(т) с их производными |
я(т), г(т) |
и мощностей со(т) =г(т)я(т), |
которые определены на |
интервале |
/о ^ т ^ < о о , принадлежат ли |
нейному функциональному пространству Ф, состоящему из кусоч
но-непрерывных функций с разрывами в счетном |
множестве |
М0 |
|
фиксированных точек полуинтервала t0^ r < ti, |
где |
t\> t — |
на |
ибольшее значение параметра t. |
|
со(т), г(т) и |
|
Для двух любых наборов функций из Ф — я(т), |
|||
я'(т), г'(т), со'(т) — разности процессов, реакций |
и мощностей |
бя (т)= я' (т)—л (т), бг (т)= гх (я + бл)—гт (я), бсо (т)=б (г (т) я (т))=г (т) бя (т) 4- бг (т) я (т) + бг (я) бя (т)
принадлежат Ф, причем бл (x)=d/dx • бя(т).
2. Из тождества (10.30) и определения следует, что для любо
го фиксированного я(т) |
существует функционал V1(я)= V (t) е Ф, |
а для любых я(т), бя(т) |
— приращение 6V'=V* (я+ бя)—У '(я)= |
= б у (t) е Ф.
3. По определению функционала г'(я), который при заданном
я(т) равен г(т)еФ , для любых заданных процессов я(т) |
и для |
произвольных процессов бя (т) справедливо тождество |
|
6V" (0—бсо (0 = 0. |
(10.42) |
Вариационный термодинамический принцип. При заданном функционале I/', определенном на Ф, в фиксированном процессе я(т) и в произвольном ненулевом процессе бя(т) существуют первые вариации, или ненулевые линейные функционалы U(8n), Ь'(бя), такие, что
V* (я + бя)—V‘ (я)= U (бя (т)) + о (бN),
г' (яЦ-бя)—г' (я)=Ь' (бя) + о ( Ш )
по некоторой норме 8N, определенной на Ф и имеющей вид
sup |kj6n (т) | + sup 162бл (т) |, т е (/0 |
t) \M 0, |
где k\>0, ki>0 — конечные величины. |
реакции г'(я) удов |
При заданном У'(я) и 6iV-»-0 функционал |
летворяет функциональному тождеству относительно бя, -которое
соответствует |
(10.42) |
|
|
6V' (я (т))—г* (я (т)) бл (0—бг' (л (т))я (0 = 0. |
(10.43) |
||
Здесь 61^', бг' суть L', L', т. е. первые вариации |
V‘ и г'. |
|
|
Теорема . |
Если множество функций я(т) |
совпадает с прост |
|
ранством непрерывно дифференцируемых функций, то |
|
||
t |
t |
|
|
6V' = I |
W dxv (*, x) + a (0 бя (/„) = \ 8nu (x) dxv‘i (t, т)+ а‘1(t) 6п4/(*0), |
to |
Г. |
|
t |
|
бг' = J бя (x) dxR (t, x) + бл (t0) В (0, |
|
to |
t
бг'я (0 = 1 бя(/ям (0 dxR‘M (t, T) + бя;/ (/0) nhl (0 Bitkl (0, (10.44) to
где индексы i, j, k, I принимают значения 0, 1, 2, 3, вектор v(t, х) (известен при данных V\ л) и тензор-функция R(t, х) с ограни ченной вариацией по х определены с точностью до аддитивной постоянной (зависящей только от t).
Рассмотрев множество процессов 6я(т), проходящих через точ ку t0>получим из (10.43) функциональное уравнение для реакции
г '(я (т )):
v (t, т)—я (t)R(t, т)—гh(t, т)= с(0 , |
(10.45) |
где c(t) — постоянный (по т) вектор
С(t)=v(t, t)— n(t)R(t, t)— rh[t, t)=v(t, t0)~n(t)R(t, t0)— rh(t, i0)
и rh(t, T ) — вектор скачков функционала г*(я(т)), который в простейшем случае равен
Г*(*, x)=T(t)h(t, т),
причем h(t, т) = 1 при x = t и h(t, |
т) =0 при т <t. |
|
При еще более общих предположениях установлены условия, |
||
при которых 1/'(я) = \ i ‘ (n)n(t) |
и тождество (10.30) |
в виде усло |
вия ортогональности |
|
|
(И (я)—г'(я))я(/) = 0 |
(10.46) |
имеет единственным решением г'(я) = V /(я).
Изложенное выше относится к термомеханическим процессам, определяемым деформациями, напряжениями, энтропией и темпе ратурой. Если существуют другие физические поля, определяемые
параметрами рг |
(г= 1 2, |
...), включающие и структурные типа |
mk в уравнения |
вносятся |
поправки. Предполагается, что возни |
кают дополнительные силы, действующие на весь объем мыслен но выделенной частицы массы рДг/, совершающие за время Ы ра боту
и что возникает дополнительный к б'QT приток тепла pq$tAV за счет проникающего облучения или источника диссипативного ха рактера. Например, в случае электромагнитного поля с электри ческим вектором Е и магнитным Н (§ 22)
6'Лр=№рЛ = (1/4я) (EdD+ HdB), |
(10.47) |
где D В — векторы электрической и магнитной индукции, т. е. набор параметров рг состоит из компонент двух векторов: р= = (D, В); объемная плотность источника тепла рq$ представляет джоулево тепло.
Первый и второй законы термодинамики получаются в обоб щенном виде:
|
р |
01 |
01 + W» ~ div Ч + Р9р. |
|
(10.48) |
|
|
|
|
|
|
|
|
|
|
Рт ■$“ = |
—div ^+ Р^р+ рш* ’ |
|
|
|
|
|
01 |
|
|
|
|
причем входящие |
в них функции |
или функционалы |
(и, |
s, w*, |
||
qt, |
) определяются процессом |
изменения функций |
времени |
|||
&{х), |
Т(т), р(т). |
|
|
|
|
|
Интегральная форма первого и второго законов для фиксиро |
||||||
ванной движущейся массы |
с объемом Vg и поверхностью |
2 g по |
лучается из (10.48) интегрированием по объему и из определения внутренней энергии Ug и энтропии Sg рассматриваемой массы:
Ut (t)= f pudV, Sg |
рsdV |
(10.49) |
Из первого равенства (10.48) и (8.29) |
получаем |
|
dUg= J (Sild£tf—div qdt + Wg/dt + pq^dt) dV = ve
— b' A— dKg— J qnd2dM-6'lp,
или согласно (10.2')
d (Ke+ Ug) = 8 ’A + 8'Qr + б'Лр,
(10.50)
U A ^ ^ W ^ + pq^dVdt.
Здесь Kg — кинетическая энергия, б'A — работа массовых и по
верхностных сил, б'QT — приток тепла, б'Ар — приток |
энергии, |
|
связанный с параметрами р. |
с множителем |
Т~х и |
Интегрируя второе уравнение (10.48) |
||
учитывая преобразование |
|
|
§ j r d i v q d V = \ - L q n d 2 + Г |
q grad TdV, |
(10.51) |
Vg
получим интегральную форму закона баланса энтропии
d S g = - ^ - L q n d m + ^ - L (w*— -yqgradT + Pq^dVdt. (10.52)
|
|
vg |
Согласно |
закону |
Фурье —q gradT=A(gradT)2> 0 . Сумму |
w*— q grad Т/Т ^ 0 |
называют иногда объемной плотностью не |
компенсированного тепла в термомеханических процессах. След ствие (10.52) в виде
s g + £ j r |
qnd S - j jr M td V = |
£ ( —qf - - + ^ ) dV> 0 (10.53). |
* |
vg |
ve |
называют иногда неравенством Клаузиуса—Дюгема. Оно неопре деленно, если не определена энтропия S (10.49) или w*.
Отметим еще одно постулируемое фундаментальное неравенство для замкнутых по деформациям изотермических процессов [13]: работа внутренних напряжений в таких циклах неотрицательна
^ SiJ dei/ = ф |
0, |
Т = const. |
(10.54)- |
§ 11. ЗАМКНУТЫЕ СИСТЕМЫ УРАВНЕНИЙ МСС |
|
||
Для малой окрестности |
физической |
точки (частицы) |
среды |
установлены дифференциальные и интегральные уравнения сохра нения массы, импульса (уравнения движения), сохранения энер гии, баланса энтропии (уравнение притока тепла), а также урав нения, связывающие тензор напряжения и вектор теплового по тока с деформациями, температурой и немеханическими задан ными параметрами. Эти соотношения в принципе определяются,, и притом однозначно, непосредственно в «^-опытах для всех воз можных в частице процессов; поскольку все входящие в эту сис тему равенств параметры измеряются приборами и системе удов летворяют, группа параметров, названная реакцией (г), одно значно определяется группой процесса (я). Следовательно, для малой частицы решение существует: r(t) = г /(я(т)). Поэтому пе речисленная система уравнений в МСС называется замкнутой для всех внутренних точек области движения среды.
В частице состояние в любой |
момент |
однородно |
(одинаково |
|
во всех точках частицы), поэтому для нее |
все |
параметры внутри |
||
и на границе частицы одинаковы |
(деформации, |
напряжения, тем |
||
пература, градиент температуры |
и т. д.). В конечных |
областях |
движения среды, кроме задания замкнутой системы для внутрен них точек, необходимы граничные и начальные условия.
В силу основного постулата МСС решение замкнутой системы уравнений МСС для внутренних точек существует при некоторых
начальных условиях (/ = /0) |
и условиях на границе области дви |
|
жения или рассматриваемой области. |
|
|
Адекватность методов |
Лагранжа и Эйлера в МСС позволяет |
|
пользоваться любым из них. Сначала будем следовать |
методу |
|
Лагранжа и систему координат х/ в начальный момент |
времени |
|
считать декартовой ортогональной. Тензоры напряжений |
и дефор |
маций по-прежнему обозначим S, <£. Общие уравнения МСС для внутренних точек области движения следующие. ^
1.Уравнения, выражающие тензор деформаций <§ и плотность
•среды р через вектор перемещения u(x, t) или х=<р=х+и (§ 4):
|
Ъ |
gii= A U h Aki= - 1£-; рЛ =р0; |
(11.1) |
|
|
|
|
|
|
|
I= d efu ; ег/= ( £ ;/-6,/)/2; |
Y g = A = \A * \, |
|
|
Эти уравнения алгебраически выражают е»/ и р через первые про |
||||
изводные |
дх‘/дх! = A\j= 6/ + ди'/дх'. |
|
|
|
2. |
Уравнения движения, которые эквивалентны одному вектор |
|||
ному уравнению |
(§ 8): |
|
|
|
|
|
V(S‘ + p (F - * ) = 0, |
- ^ - = х , |
(11.2') |
причем F= F(x, t) — заданная массовая сила. Уравнения записы ваются также в виде
|
po( ? - X O = |
^ ( A 4 'S '" 0 . |
|
|
(Н-2) |
Структура этих трех скалярных уравнений, очевидно, такова: они |
|||||
линейно связывают между собой напряжения |
S ij |
и их |
первые |
||
производные S% со вторыми производными |
х по времени, а так |
||||
же первыми и вторыми производными х по координатам. |
|
||||
3. |
В тех случаях, когда за |
независимые искомые функции при |
|||
нимаются деформации гц или напряжения |
Sl/, |
должны |
быть |
||
удовлетворены еще условия совместности (§ 5) |
|
|
|
||
|
Rlmpq—О, |
|
|
|
|
где индексы Itnpq имеют шесть независимых комбинаций |
(1221 и |
||||
т. д.). Тензор кривизны RimPq алгебраически |
выражается |
через |
|||
е//, е//, k и 6//, kmy причем через вторые производные от е;/ — линей |
|||||
но. Эти шесть уравнений кратко запишем в виде |
|
|
|||
|
R = R (ei/t |
eijtkin)= 0. |
|
|
(11.3) |
Существенно, что если, как и сделаем, за независимую искомую функцию принять вектор х или перемещение и, то на основании (11.1) уравнения (11.3) обратятся в тождества и потому могут не рассматриваться.
Следовательно, для вектора и и тензора напряжений S полу чается только одно уравнение (11.2'), т. е. для девяти функций
х1, ха, х3 S11 S22 S33 S12= S 21 S23= S 32 S31 = S13
мы из (11.2) имеем только систему трех дифференциальных урав-