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

Курс механики сплошных сред

..pdf
Скачиваний:
19
Добавлен:
15.11.2022
Размер:
20.58 Mб
Скачать

Поясним^ это на примере. Предположим, что среда имеет контакт с подвижной стенкой, которую можно приближенно принять за плоскость *з = 0. Если уравнение подвижной стенки записать в виде x%—g(xv хг, 0 = 0. то следует предположить, что функция g и ее производные имеют порядок 0 (т)), чтобы иметь возможность при­ менить линеаризованную теорию. Краевое условие (11,55) запишется в этом случае следующим образом:

 

 

 

U*(Xn

хг,

g) — jj = 0.

 

 

 

В линеаризованной

теории

это

условие запишется проще:

 

 

 

 

иАХг, *„ 0 ) - | = 0,

 

 

(64)

т. е. оно формулируется на невозмущенной

границе

х8 = 0.

 

Замечания. Линеаризация

может быть распространена на случай

изучения системы, движение

которой мало отличается от некоторого

известного движения

(см. V.7). В рассмотренном

выше простом слу­

чае такое движение было просто покоем.

 

 

 

 

Речь здесь идет о формальной асимптотической теории, главная

идея

которой основана

на изучении

зависимости

семейства движе­

ний системы S от параметра т], когда невозмущенное движение,

соответствующее г] = 0,

известно. Ставится

задача

рассчитать

воз­

мущение для каждой

искомой

величины,

т. е. найти первый

член

асимптотического разложения

этой

величины по т]. Однако на

прак­

тике

существенного

упрощения

в

решении не

получается,

если

асимптотическое разложение производится на основе полного реше­ ния всей совокупности уравнений, зависящих от rj.

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

ными по г],

производится соответствующая подстановка в уравнения

и

краевые

условия и в каждом из равенств берутся только главные

по

параметру г\ члены. Таким образом, приходим к уравнениям и

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

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

Все сказанное выше показывает ограниченность области приме­

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

ДОПОЛНЕНИЕ

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

носа. Для изучения

последующих

глав выводы раздела V.5 не обязательны, но

их значение

в общем

случае и роль

в современных теориях механики сплошных

сред таковы,

что их

необходимо здесь упомянуть.

 

 

V.5. ИСПОЛЬЗОВАНИЕ ПЕРЕМЕННЫХ ЛАГРАНЖА

V.6 .I.

 

Формулировка закона сохранения. Вернемся к формулировке закона

сохранения,

описанного в 11,1

в переменных

Эйлера,

и выразим

его в перемен­

ных Лагранжа. Упростим

несколько задачу, заменив

объемные плотности Jli и

Ai массовыми 59,-

и Bit и предположим,

что они выражены в переменных а и t

(х = Ф (a,

t):

 

 

 

 

 

 

 

 

 

 

 

 

 

p93i(а,

0 =

с/£/(х, /),

рВ( (a, t) = Aj(x, t).

(65)

Обозначив через *£)° область, занятую отсчетной конфигурацией 5°, опреде­

ляющую

движущуюся

область

а через р0 (а) объемную

плотность в этой кон­

фигурации,

нетрудно

преобразовать тройные

интегралы из

(11,1) с

учетом закона

сохранения

массы;

получаем

 

 

 

 

 

 

 

 

 

dt jj#

dv~ dt J

Po® ‘ du°' §х>А‘

p°s ' dv°'

(66)

Здесь

du°— элементарный

объем области

кроме того, как

указано в (6)

 

 

 

р0 = р /,

y = det(F),

f ia = ^ L = ^ L .

 

(67)

 

 

 

 

 

 

 

 

 

оаа

оаа

 

 

Остается преобразовать

 

поверхностный

интеграл

в (11,1). С использованием

внешних дифференциалов подынтегральное выражение запишется [с учетом (11,35)] в виде

 

 

ш/ = у

6XQА 6хг= ЧрФ/>;

(68)

величину

(0/ = а,уЛу

можно интерпретировать как вектор

элементарных потерь

через границу д<2),

а у р = пр до— как

вектор элемента поверхности. Замена пе­

ременных

в фр производится согласно обычным правилам

 

—2 гряг?я$FГу^а$A dtf^.

С другой стороны, из уравнений (67) и (П1.44) следует выражение

F р а У р ж ~2 Bp q rF p 0 L ^q $ F гу

Л

ea0Y^ ^fl0 ^

в котором фа— вторая форма элемента

площади *

Кроме того, величины ф°а

можно трактовать как компоненты вектора

элемента

поверхности Ла do0, где л£,

в свою очередь,— компоненты единичного

вектора,

нормального к элементу по­

верхности д*2)°, площадь которого da0 и который соответствует элементу поверх­

ности д&) движущейся системы. Таким образом,

можем

окончательно написать:

<>)i ==aip<?P=biaK'

 

(69)

полагая

 

 

^ &ip~ ^ ia Fpa* ^ict. == ^

*

(70)

где Gaj— элементы матрицы, обратной матрице градиента F.

Итак, матрица bia позволяет рассчитать ю,-, используя вектор ф£ элемента поверхности отсчетной конфигурации по формуле, аналогичной той, которая слу­ жит для расчета со,- по вектору ф^ элементарной поверхности в текущей конфи­ гурации.

С учетом введенных обозначений закон сохранения (11,1) записывается в окон­ чательном виде

^ p05|du° + J Ь{а Па d o °= ^ р0В,- dv°.

(71)

аа° а

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

 

 

 

 

 

Ро -gjf

 

a =

 

Ро^/»

 

 

 

(72)

так как

субстанциональная

производная

в

лагранжевых

переменных совпадает

с частной

производной

по

времени,

а р0

от

 

времени

не зависит. Заметим,

что

 

»

 

 

 

 

 

 

 

 

гт

 

 

 

 

bia, p =

, причем bia

выражены в переменных Лагранжа.

 

 

V.5.2

Приложение

к

механике

сложных сред. Уравнения

в переменных

Лаг­

ранжа. К вектору

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

 

напряжения

т м о ж н о

применить резуль­

таты (69)

и (70). Для

этого

следует

ввести

матрицу

tia,

называемую матрицей

Пиола— Лагранжа,

которая

дает

 

возможность

по-новому

представить тензор на­

пряжений:

 

 

 

 

 

 

 

 

 

 

 

 

 

(73)

 

 

 

 

 

Ту ==»О/уфу — hatya у

 

 

 

 

 

 

 

J(Jij =

tiaFia> tia — JQijGaj\

 

 

 

(74)

последние формулы можно переписать

в матричном виде

 

 

 

 

 

 

 

 

/ 2 =

 

TFT, Т = / 2 б т .

 

 

 

(75)

Напомним, что матрица Т **

позволяет выразить элементарное напряжение т

через вектор ф° элемента поверхности отсчетной конфигурации.

 

Из уравнения

(72)

 

следуют

уравнения

сохранения

количества движения

в лагранжевых переменных. Прежде всего из уравнения (111,1) получаем

 

 

 

 

 

 

Л2у ,

 

 

 

 

 

 

 

 

 

 

 

 

 

Р° ^

= Pogi + tia, a.

 

 

 

(76)

так как в переменных Лагранжа

у ,=

d2Xj

gi

 

обозначает

здесь массовую плот-

dt2

 

* Напомним, что аа — декартовы координаты отсчетной конфигурации.

 

** Символ "I" обозначает здесь матрицу, составленную из элементов

а не

матрицу-столбец, относящуюся к вектору Т.

 

 

 

 

 

 

 

 

ность заданных внешних сил (/,• =

 

 

 

 

 

 

 

= р?/). Далее, формула (111,2) выра­

 

 

 

 

 

 

 

жает, как видно, симметрию тен­

 

 

 

 

 

 

 

зора напряжении,

т. е.

 

 

 

 

 

 

 

 

 

UaFia =

tjaFia

 

(77)

 

 

 

 

 

 

 

На самом же деле матрица

 

 

 

 

 

 

 

Пиола — Лагранжа

не является

 

 

 

 

 

 

 

симметричной, что ведет к опреде­

 

 

 

 

 

 

 

ленным неудобствам

при

се

ис­

 

Рис.

4.

Тензор напряжений:

пользовании. Кроме того, как и

 

матрица F,

матрица

Пиола — Лаг­

п и п ° - единичные

векторы; ф и т получаются из

ранжа неудобна тем, что ее тензор­

ф° и т° путем переноса.

Название

ный характер распознается нелегко.

Отображение

Матрицы и

Последнее

следует

хотя

бы

из

ф — > X

 

компоненты

Тензор Коши

присутствия

латинского

и грече­

 

2,

о|7

ф — ►т

 

Т . П а

Матрица Пиола -

ского индексов: первый обозначает

•Ф*1°

 

 

 

Лагранжа

операцию в Т (Mf) — линейном ка­

 

S. %3

Тензор Пиола —

сательном

пространстве частицы М

 

 

 

 

Кирхгофа

в текущей конфигурации, а второй

 

 

 

 

 

 

 

относится к операции

в векторном

касательном

пространстве отсчетной

конфигурации.

 

 

 

 

 

 

Эти

неудобства могут

быть

устранены,

если рассмотреть в векторном прост­

ранстве

Т (М°),

касательном в точке М° (положение частицы в отсчетной конфи­

гурации), элементарный вектор /°, который

в результате

конвективного

переноса

переходит в вектор т в текущей

конфигурации. Согласно

(5)

имеем

 

 

 

 

 

 

 

X{ = F ta1a, Х° = Ga (T,-,

 

 

 

 

 

(78)

и, следовательно, в силу (73) и (74):

 

 

 

 

 

 

 

 

 

 

Та =

CWfP<P(3= *7GaiO/yGp/<pp= sapq)p-

 

 

 

 

(79)

Положим в матричном обозначении

 

 

 

 

 

 

 

 

 

 

 

 

S = GT = JG2GT ;

 

 

 

 

 

(80)

-заметим также,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J2 =

FSFT ,

Т = FS.

 

 

 

 

 

(81)

Матрица S

составленная

из

компонентов sap, представляет в пространстве

Т (М°)

симметричный тензор, называемый тензором Пиола— Кирхгофа и выражаю­

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

Точно так же может быть преобразовано уравнение энергии. Элементарное количество теплоты, полученное областью *2), запишется [согласно (IV.2.8)) в виде

если

положить

gdo = — <7,71,-do = — <7/ф/ = —(Оафа,

(82)

 

 

 

 

 

 

 

(ba = JqiGai

Jqt = toaFia-

(83)

Вектор со

является

вектором

потока теплоты Пиола— Кирхгофа.

ных

В переменных Лагранжа уравнение сохранения энергии в Частных производ­

будет

иметь вид

 

 

 

 

 

Ро "dГ

( в

)

i)*a ~

a i-r ,

в котором все величины

предполагаются выраженными в

переменных Лагранжа

и, следовательно,

полная

производная

фактически является Частной. Как и в

разделе IV.3.3, это

уравнение может

быть представлено

в бо^ее простом виде,

Ро d<

dT F‘a _

“ + Л

f84)

V.5.3. Сопряженное представление напряжений и скоростей деформаций. Заметим, что понятие массовой энергии внутренних сил, введенное в IV.3.3, при­ водит к произведению матриц

P~'°ijDij = Paltia

* 5 * -р ,- Ч ч >

 

(85)

В самом деле, последнее равенство следует из уравнений

 

^ice

 

 

== ~7£

/а)*

 

 

р dFia 1

/

dF ia ,

c

1

d /v

x

sa P ^ f l - ^ = 2 - S a p ( ^ / 0

--- h ^ .’a “ d/

SaP

d7

 

последнее вытекает из симметрии sap.

dp"1 dL

Принято говорить, что матрицы D, —jj- ,

, представляющие собой ско­

рости деформаций, сопряжены соответственно матрицам 2, 7\ S, которые опре­ деляют три способа описания напряжений. Примечания в IV.2.3 показывают преимущества этого определения.

V.6. ПРОИЗВОДНЫЕ ПО ВРЕМЕНИ ВЕКТОРА И ТЕНЗОРА

Выше было введено понятие полной производной. Фактически она представляет собой производную, полученную для частицы при ее переносе. Если, например, изучается вектор X (/), заданный в векторном пространстве Т (Mf)t то имеем

 

 

 

 

“ .

i .

' i

+

t

a

,

 

 

 

 

 

 

 

 

, »

о

 

 

 

 

 

Здесь

следует иметь

в виду,

что для

вычисления

правой части

в пространстве

Т (М*) взят вектор,

который

эквиполлентен

вектору

X (t + s)

из Т (A4t+S).

Но такой вектор не может характеризовать изменение X относительно самой среды,

Поясним это на следующем простом

примере. Пусть S — абсолютное твердое тело,

движущееся

в системе

отсчета

тогда

полная

производная от любого вектора,

постоянного

в 5,

не будет равна

нулю

в

а

будет равна его абсолютной ско­

рости

в системе

С физической же точки зрения является важным найти изме­

нение

изучаемой

величины

при

рассмотрении

явления

в самой

среде, так как

именно такие производные могут играть определенную роль в формулировке законов поведения.

V.6.I. Конвективная производная. Второй путь заключается в обращении к понятию конвективной производной, которая в случае вектора X определяется по формуле

 

 

DCX =

lim

- (< ^ -i (t, t + s) {X (t + s ) ) - X (0).

(86)

 

 

 

s- 0

S

 

где

— отображение,

обратное с£Г, введенному в V.1.1.

в базисе

Согласно (5)

и (11) компоненты отображения оТ -1 (X (^ + ь)} в (86)

равны

 

 

Gik(t, t + s)Xk (t + s).

 

 

 

 

 

 

Так

как G/^/7^y=

6/y, то

после дифференцирования в соответствии с

результа­

тами,

полученными в V.3, имеем

 

±0„ ч.

Теперь можно представить правую часть (86) в явной форме, пользуясь не-

X k (t + s) = Xk ( / ) + s ^ - * ( / ) + . . . ,

@ife(^» ^ + s)= $/ft— sUit*+•••

Обозначая составляющую вектора DCX/ по направлению в,- через DCX, можем написать

о Л = ^ - (/', Л -

<87)

Точно так же можно рассуждать и в отношении тензора второго ранга. Например, конвективная производная тензора напряжений запишется в виде

 

 

боц

Ui,k°k/ U/,k°ik'

(88)

 

 

=

Формулы

(87) и (88) определяют конвективные производные вектора и тен­

зора второго

ранга

соответственно. Несмотря на важность этого

понятия, на

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

не будем, так

как для правильного понимания

и примене­

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

V.6.2. Производная относительно собственного вращения или производная Яуманна. Третий путь, который позволяет ввести производную по времени, обхо­

дящий трудность, отмеченную в связи

с определением

конвективной

производной

и тесно

связанный с деформацией среды,

заключается

в том, что конвективный

перенос

(М, t', /) заменяется переносом

в собственном вращении 31 (М,

Z',

t),

введенным в V.I.3. Производная

в

собственном

вращении

вектора

X (t)

будет,

таким образом, задаваться

формулой

[аналогичной (86)]

 

 

 

 

 

 

D jX = Urn

1

(Я "1 (М,

t,

l + s) {jr(/ +

s ) } - X (0).

 

(89)

 

S -*• 0

S

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Составляющая 31"1 {X ( t s ) }

по

 

направлению

б/ будет

 

 

 

 

что следует из R~1 = RT с

Rki(t> ^+

s) (^fe (t +

s),

 

 

 

 

 

 

 

учетом

 

определения

(23).

С

другой

стороны,

со­

гласно

(32)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0 =

^ft/ = — &ik-

 

 

 

 

 

 

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

относительно собственного

вращения

вектора и тензора:

 

 

 

 

 

 

 

 

DjX/~ W

 

QikXk,

 

 

 

 

 

 

 

(90)

 

 

 

 

do,-у

 

 

 

 

 

 

 

 

 

 

(91)

 

DJatJ=

~di----- ®ikPk/®/k°ik-

 

 

 

 

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

* Строго говоря, введенная конвективная производная является контравариантной конвективной производной. Последовательное развитие понятия конвек­ тивной производной требует обращения к другим производным вектора или тен­ зора, в частности к ковариантной конвективной производной (см. задачу 25). Заметим, наконец, что введенное понятие конвективной производной — это произ­ водная Ли, ассоциированная с группой П, зависящей непрерывно от некоторого параметра и определенной на пространственно-временном многообразии, задавае­ мом движущейся системой S, поле скоростей U которого является ее бесконечно малым генератором.

также производной ЯуМанна) эти члены вычисляются через компоненты скоростей

вращения. Заметим, наконец, что выражение (90) может быть переписано с уче­ том (IV, 20):

Если же сплошная среда— абсолютно твердое тело 5, то последняя формула совпадает со следующей хорошо известной в общей механике:

в которой

о — вектор

скоростей вращения абсолютно твердого тела, движущегося

в системе

отсчета £j{,\

последний член является, как

известно, переносной ско­

ростью

вектора X.

 

 

 

 

 

V.7. ОБЩИЙ ФОРМАЛИЗМ ЛИНЕАРИЗАЦИИ

В

V.4.5

показано, каким образом осуществляется

линеаризация в том слу­

чае, когда конфигурации движущейся системы остаются

близкими к фиксирован­

ной заданной

конфигурации. В механике сплошных сред и, в частности, в гидро­

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

состояние

системы,

предполагаются известными.

 

 

Само

собой разумеется, что все замечания

V.4.5 относительно формального

характера

линеаризованной теории

остаются в силе в более

общем случае, кото­

рый рассматривается

ниже.

рассмотрения двух примеров, в которых пред­

Здесь

ставится задача — после

ставляется

интересным проследить

применение

линеаризованной теории

сформу­

лировать логически

естественные

правила математического

формализма

линеари­

зации и привести несколько приложений.

V.7.I. Типичные примеры, а) Линеаризованная аэродинамика. Рассмотрим следующее невозмущенное движение: идеальная жидкость находится в установив­ шемся движении с постоянной скоростью параллельно оси хи занимая все про­ странство вне жесткой неподвижной пластины с произвольными очертаниями,

лежащей в плоскости

*3 = 0.

Допустим

временно,

что такое движение возможно

и что при этом выполняются

условия,

 

налагаемые

присутствием

пластины

(нор­

мальная к пластине составляющая скорости равна нулю).

 

 

 

Представим теперь

себе жесткое препятствие А (крыло самолета), геометриче­

ская

форма которого

близка

к форме

пластины

Р (рис. 5). Для конкретности

будет

считать, что

все

точки

тела

А расположены

близко к Р и что во всякой

точке

поверхости

препятствия

дА,

в которой существует касательная

плоскость,

угол между последней и поверхностью

Р очень мал. Образно говоря,

тело

А по­

лучено малой деформацией и

путем

«покрытия»

пластины Р.

Поместим

наше

препятствие в упомянутое равномерное течение так, чтобы все касательные пло­ скости были близки к плоскости *з = 0; и будем считать, что вокруг А установи­ лось новое стационарное течение. Это течение не будет уже равномерным, но можно считать, что оно мало отличается от начального равномерного невозму­ щенного течения. Вследствие этого можно ожидать, что линеаризация относительно невозмущенного движения позволит получить полезную информацию.

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

. Вебозмущенное течение

Возмущенное течение

Рис. 5. Пример задачи линеаризованной аэродинамики.

Для упрощения взяты плоские течения

б) Устойчивость ваданного течения. Рассмотрим стационарное течение вязкой жидкости между двумя соосными круговыми цилиндрами с осью 0х3, в котором линии тока являются кругами с центрами на оси 0х3. Установлено, что течение такого вида существует, и все условия задачи при этом выполнены. Именно это течение выберем за невозмущенное. Пусть в некоторый заданный момент времени /0 течение претерпело малое возмущение; можно считать, что в последующие моменты (если, по крайней мере, t 10 достаточно мало) траектории частиц нового течения будут мало отличаться от траекторий невозмущенного движения. Решение линеа­ ризованной задачи позволяет получить информацию об устойчивости начального течения. В самом деле, можно искать ответ на вопрос о том, затухает или растет

возмущение, определяемое линеаризованным решением.

V.7.2. Определение линеаризованной задачи, соответствующей некоторой фиксированной задаче. В этом разделе изложим, каким образом по определению

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

а) Уравнения.

Обозначим

условно

через

fb (х, t) функции, известные

в рас­

сматриваемой задаче,

через

Ха (х.

t) — функции,

являющиеся

ее

решением

(т. е.

искомые величины);

и наконец,

(§а {/*,»

Ха} = 0

будет обозначать

уравнения, свя­

зывающие неизвестные с заданными величинами в каждой точке области

занятой системой S

в момент t.

Индексы а и b

могут принимать

конечное

число

значений.

Обозначим

через fb и Xl

заданные

величины и решение невозмущен­

ного уравнения в области So, удовлетворяющие в ней уравнению

$ а {/?» Х2} = 0,

Введем теперь

функции

fb (х,

t)

и

Ха (*,

t), определенные в области 5$ для

любого момента времени,

и

рассмотрим

выражение

 

 

 

 

 

 

 

 

 

$сс

 

 

Х2+ Т|Ха),

 

 

 

 

равное нулю при т] = 0.

Пусть

для

малых т| можно написать главную

часть дан­

ного выражения как

 

По определению, уравнения

 

 

 

 

 

 

 

 

 

 

 

 

X®, 7ь. Х а) = 0,

 

 

 

(92)

очевидно,

линейные

и

однородные

относительно и Ха, составляют

линеаризо-

ванные_ уравнения,

соответствующие

данному

невозмущенному

движению.

Вели­

чины fb всегда будут считаться заданными, а Ха—неизвестными линеаризованной

задачи. По определению, г\Ха характеризует возмущение, a Xa-\-i\Xa— линеари­ зованное возмущенное движение,

б) Граничные условия. Рассмотренное начальное движение удовлетворяет не только системе уравнений ^ а = 0, но также и некоторым граничным условиям на поверхностях в любой момент времени. Одна из таких поверхностей (границ) 2* будет представлена параметрически уравнениями x ,= g /(£ b £2, t).

Граничные условия можно записать символически:

2fi{Xat Fd, Gc (£/)} = 0.

Операторы jgp определены на

2*; Ха (£ь £2, /)— значения неизвестных на 2*;

величины Fd (h, 12, t) известны;

Gc — некоторые функции, определяемые геомет-

*з{*2. Ъ. осШ = о,

где

= g °(£ 1» ?2» 0

определяет

граничную

поверхность

2$,

которая

совпадает

с Г в невозмущенном движении.

 

 

 

 

 

 

 

 

 

 

 

 

Введем,

как и ранее,

фукции Ха (gb

g*

/),

Fd (glf g2, 0. g/ (Si,

S2,

0

и рас­

смотрим

выражение

S^{X°a+ 7]Xat F°d+r\Fd,

 

(#? +

Л&/)}, равное

нулю

при

Л = 0.

Предположим

далее,

что

главная

часть

в разложении

по

i] равна т)Za.

По определению, уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

03

№ .

F \ g\. Ха, Fd, 8 i\ = 0

 

 

_

 

 

(93)

выражают

условия,

которым должны

удовлетворять

неизвестные

Ха линеаризо­

ванной

задачи (остальные

аргументы,

как

правило,

известны)

на

граничной по­

верхности

2о области

SQ. Функции g p — аффинные функции

переменной

Ха.

 

Итак,

решение

линеаризованной

задачи

заключается

в

том,

 

чтобы

найти

в So неизвестные функции Xat удовлетворяющие в So уравнениям (92), а на гра­ нице 2о — граничным условиям типа (93).

V.7.3. Движения, близкие к равномерным. Проиллюстрируем изложенные выше формальные положения на примере V .7.1, а (линеаризованная аэродинамика). Здесь компоненты поля скоростей невозмущенного движения

Ui = V0, U2 = 0, U5 = 0;

возмущенного движения

^1 = Уо + Л^1* ^а = Л^2» ^з = Л^з-

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

* '= 1 r + U‘U‘f

Удерживая в линейной теории только члены с г], находим составляющие ускорения

Заметим, что и в общем случае для нахождения полной производной некото­ рой величины, равной 0, или постоянной в невозмущенном движении необходимо

применить к этой последней оператор

Точно так же поступим со всеми

членами

уравнений. В качестве

примера

рассмотрим уравнение неразрывности, общий вид которого

 

^ + ( Р ^ , ) . , = 0.

 

(94)

Пусть ро— объемная плотность жидкости в невозмущен ном движении,

которую

будем считать постоянной. Объемную

плотность в возмущенном движении запишем

какро + ЛР» и после линеаризации уравнение (94)

примет вид

 

% +У°Ё^ + р° и‘-1=0^

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

*з = Л£(*ь *2, 0.

близкой по форме к плоскости *3«=0. Граничное условие, естественно связанное с уравнением неразрывности (нормальная составляющая скорости жидкости или

газа относительно стенки), запишется в общем случае следующим образом (II, 55):

-jj- (л5—*з)=о,

или в развернутой форме

 

 

 

п { - ^ + [ У о +

Л*М *1. -<2- ni)]

(*и

— j\U3(xi,xi t ^g) = 0.

Удерживая

только первый член в разложении по т], получаем граничное

усло­

вие линеаризованной задачи:

 

 

 

 

и а (*х, Jt„ 0 ) = § ( х и

X», 0 +

V « -|£ (*i. Xt, t).

(95)

Из этого уравнения можно получить информацию (в данном случае величину

компонента U3 возмущенной скорости) на поверхности дг3 = 0, т. е. на границе соответствующей границе 2* в невозмущенном движении.

ГЛАВА VI

ЗАКОНЫ ПОВЕДЕНИЯ (ЧИСТО МЕХАНИЧЕСКИЕ ЯВЛЕНИЯ)

СОДЕРЖАНИЕ ГЛАВЫ

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

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

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

В последующих главах это ограничение будет снято после вве­