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

Механика сплошной среды

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

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

Перечислим основные механические граничные условия в за­ дачах МСС в предположении, что граница области тела 2 задана

уравнением

Ф(дс,

t)= О в эйлеровом

пространстве и единичная

нормаль v = п к ней известна:

 

 

 

 

 

V(JC

t) = v = в1” **

на Ф(х,

0 =

0.

(12.1)

 

 

I grad Ф |

 

 

 

 

Поверхность

Ф=0

представим состоящей из

трех

частей:

2 Г 2 Р,

2 vp-

Кинематические условия на 2 Wв зависимости от свойств среды могут сводиться либо к заданию полного вектора скорости v (или

вектора перемещения и), например на неподвижной

поверхности

Ф (*)=0 — условие полного сцепления

(прилипания):

 

наФ(л:) = 0 v = 0

(или и=0),

(12.2)

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

на Ф(л;) = 0 vv = 0 (или vu=0).

(12.3)

Соотношение (12.3)

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

вие непроницаемости

на движущейся поверхности

Ф(х,

t)

нахо­

дится

так: в момент

 

t+ dt частица среды

с координатой

х

в мо­

мент t

будет иметь

 

координату x + \dt\ если

она

была

 

на по­

верхности,

т. е. Ф(х,

 

/)=0, то и останется

на ней, т. е. (&(x+vdt,

t + dt)= 0,

следовательно, после вычитания

и

деления на

dt по­

лучим

 

 

 

 

 

 

 

 

 

- ^ v + ^ = 0 , dx dt

где v — вектор скорости среды. Но

-^ -= grad® ,

и потому на основании (12.1) получаем

VV

1

дФ ___i_r,(v)

I grad Ф I

 

 

(12.4)

(12.5)

( 12.6)

где — скорость движения поверхности Ф вдоль ее нормали. Условие (12.4) или (12.6) означает, что нормальная проекция

скорости

v

среды

совпадает

с

это условие

совпадает

с (12.3),

если

Ф = Ф ( х ) = 0 . Если задано движение

каждой фикси­

рованной

точки х

поверхности

Ф ( х ,

/ ) = 0 ,

т- е- задан вектор

ско­

рости \ ф ( х , t),

то

условие

сцепления

примет

вид

v(x,

/) =

=Уф(х, t) на Ф ( х , /)= 0 , при этом автоматически

выполнено ус­

ловие (12.6)

 

на Ф =0 уф (или и= 11ф).

(12.7)

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

на Ф(лг, t) = 0 P{v) •у=р{ф)

 

P(v) • v = Sl/VjV/ = оцЩП'!

(12.8)

(пример идеальной жидкости о// ——рб;/), либо к заданию полно­ го вектора Рф* внешней силы (напряжения). Поскольку на лю­

бой площадке вектор внутреннего напряжения P(v) = Slvh то ус­ ловие в напряжениях имеет вид:

на Ф =0

P(v) =

==oiinie/=Po)

(12.9)

Смешанные условия

на 2 ир

состоят в задании

на Ф = 0 двух

(или одной) составляющих вектора v и одной (или двух) состав­

ляющей вектора P(v), а всего трех скалярных условий, т. е. час­ тично — задания условий (12.7), частично — (12.9). Если а, Ь, с — три ортогональных вектора на поверхности, то три условия должны быть такими, чтобы одно относилось к направлению а, другое — к Ь, третье — к с, т. е. не было бы двух, относящихся

к одному и тому же вектору, иначе задача МСС,

как правило,

оказывается неопределенной.

граничные условия

состоят либо

Основные температурные

в задании на поверхности температуры:

 

на Ф =0

Т=Тф (л:, 0,

(12.10)

либо в задании теплового потока:

 

на Ф = 0

qv = рФ(лг, /).

(12.11)

Условие (12.10) возможно потому, что на границе двух однотип­

ных сред, разделенных

поверхностью Ф= 0,

температуры

бывают

одинаковыми и условие

(12.10)

означает непрерывность

темпера­

туры на Ф = 0. Условие

(12.11),

если только

на поверхности Ф = 0

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

образуется

(например,

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

Ф = 0), то

правая часть

(12.11) будет состоять из теплообразова­

ния (]фт и потока тепла от внешнего тела. Поток ^фг часто счита­ ется пропорциональным разности температуры тела Т на Ф = 0 и температуры внешнего тела ТФ:

где k называется коэффициентом теплоотдачи. Если твердое тело с границей Ф = 0 сильно нагрето, то в пустоту (приближенно — в воздух) оно отдает лучистое тепло\ в условии (12.11) в этом случае q<& может быть взято согласно закону Стефана—Больц­ мана:

 

Ш '

<1212'>

где Т — температура тела на поверхности

(К), с — постоянная,

зависящая

от свойств тела. Для абсолютно черного тела

с0 =

= 5,8 Вт/м2

для «серых» тел с= гс0, причем

е для разных

тел

имеет значения от 0,96 для окисленных шероховатых черных ме­ таллов до 0,3 для блестяще полированных.

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

динамические задачи требуют кроме граничных условий поста-

о

новки еще и начальных условий. В перемещениях и=х— о(JC, t) =

= и= ф(х, t)—х

эти условия,

как

и в теоретической

механике,

имеют вид

 

о

 

 

 

 

 

i= t0,

 

v=tMAr),

 

или

и = ф 1(лг),

(12.13)

 

 

 

 

 

 

t = t0,

u=

cp2 (х),

-$-=Ч>а(х),

 

 

 

 

 

at

 

 

где ф, ф — заданные векторы (обычно

ф2 = 0). Начальное усло­

вие для температуры

 

 

 

 

 

 

t = t0, Т = х'(х) = х'(х),

(12.14)

где %— заданные функции координат.

Выше рассмотрены ус­

Условия на

поверхностях

разрывов.

ловия на поверхности 2(Ф = 0), являющейся границей

рассматри­

ваемой области

G движения

среды. В динамических

проблемах

МСС необходимо учитывать возможность возникновения поверх­ ностей разрыва внутри области G. Такие поверхности могут воз­

никать в средах, обладающих упругими свойствами

(наряду с

другими свойствами, отраженными операторами

и, QT,

(§ 12)) за счет прилагаемых в какой-то момент внезапно внеш­ них сил (F, ^(v)) или других воздействий (соответствующих пара­ метрам р). Эти поверхности движутся внутри G, выходя на гра­ ницу 2, отражаясь от нее или сообщая ее неподвижным участкам сильные толчки, иногда вызывая разрушения и т. п. Движение по­ верхности разрыва называют волной разрыва в среде. Типичными примерами являются сейсмические волны при землетрясениях, ударные волны в атмосфере от взрывов и сверхзвуковых движе­ ний летательных аппаратов.

Возможность или невозможность возникновения волн в среде полностью определяется типом присущих ей функционалов состо­

яния и QT (§ И), которые в уравнениях движения и распрост­ ранения тепла дифференцируются по координатам и времени. Но именно на поверхностях разрывов они могут терпеть разрывы, и потому дифференциальные уравнения должны пониматься в обоб­ щенном смысле или заменяться интегральными. Это означает: либо решение задачи МСС, т. е. х(х, /), v(x, t), Г, р .. надо ис­ кать во всей области G в виде обобщенных функций, либо по­ верхности разрывов выделить из G и включить в состав поверх­ ности 2, на которой записываются «граничные условия», и тогда искать в получившейся области классические решения.

Второй путь в МСС пока

является более распространенным,

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

особо уравнения на поверхностях

разрывов.

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

которой терпят разрыв тензор деформации б, температура Т и, следовательно, тензор скорости деформации V и могут иметь раз­ рыв функционалы fFiJ' уравнений состояния, т. е. тензор напряже­

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

 

 

 

 

 

д Н

 

Я(дг,

0 =

0,

grad И

 

d t

(12.15)

П = V =

u < v )

=

 

 

 

grad Н |

 

I grad Н |

 

Индексом

«1»

отмечаем состояние

среды

перед фронтом

волны,

индексом

«2» — за фронтом (т. е.

уже возмущенное волной сос­

тояние). Поверхность Н = 0 предполагается настолько гладкой, что в окрестности точки х в момент t ее можно заменить касательной плоскостью и рассмотреть малый элемент площади ее поверхности

Д2//, одинаковый в момент

t

и t + dt (dt = const).

За

время dt

в неподвижном слое объема

АV// == Д2/*- v^dt пространства наблю­

дателя произойдут следующие

изменения: значения

в

момент t

плотности, скорости, температуры, внутренней энергии,

вектора

напряжения, параметров р (следовательно, тензоров деформаций,

скорости деформаций, напряжений

)

 

f>i= p(*.

t), vl = v(x,

t), Tl= T { x ,

t),

ui=u(Jf, /),

Piv)= a'4Ar,

t)tii(x, t)eh

n i= vt{

скачком изменятся и получат значения

P2=P(-*. t + dt), v2= .

P2V)= O'/ (JC, t + dt)n1(x, t)tj. (12.16")

Подчеркнем, что n = v и iAv) (12.15) за время dt изменяются пре­ небрежимо мало.

Скоростью D(x, t) распространения поверхности разрыва или ударной волны в среде называется скорость движения поверхности разрыва вдоль ее нормали относительно невозмущенного волной вещества, т. е. D=v(v>—vvi. Изменение скорости среды после про­ хождения волны обозначим V; таким образом,

 

 

 

AV =

V= V2—vlf

D = v ^ — v 1n.

 

 

(12.17)

Условия

непрерывности

перемещения

и(х,

/)

перед фронтом

//(x, t) = 0 и

за

фронтом

(индекс

2)

Ди = и2(х,

t)—щ (х,

/) = О

можно

рассмотреть

еще

в

момент

Й-6£,

когда

точке

х на

Н (х, t) =0

можно

сопоставить

точку x'=x + vv<?>8t,

которая будет

на поверхности tf(x', t+8t)=0\ скачок перемещения и будет

 

и2(лг', / + б/)—щ (х ’, t + bt) = (A —

vu<v>+ A —

)б /= 0 ,

(12.18'

 

 

 

 

 

 

 

 

\

дх

 

 

 

dt

)

 

 

 

где Д(йи/(3...) =

((3и/(3...)2 — (ди/д...)!.

 

Пусть \

— бесконечно

малый)

вектор

на поверхности Н (х,

£)= 0 в точке

х:

 

 

 

 

 

 

 

v£=0,

u2 -Ы, /)—UI (JC + 1,

 

 

ох

1=0.

 

 

Поскольку

%— любое

направление

в

(х, t)

на #= 0, то

отсюда

следует, что существует диагональная матрица Z= diag(z1

z2

z3),

такая,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A dui

 

/

 

 

 

 

 

 

(12.18")

 

 

 

 

 

 

Д ----- =

z£v,

 

 

 

 

 

 

 

 

 

 

 

 

 

дх

 

 

 

 

 

 

 

 

 

 

откуда

(умножая на v)

находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

А

ди1

 

 

А

ди1

 

 

 

 

 

 

 

 

 

 

 

zl—= Д ----- v=

Д ------

 

 

 

 

 

 

 

 

 

 

 

 

 

дх

 

 

 

д \

 

 

 

 

 

 

 

Из (12.18')

теперь имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

гЧЛ* + Д

 

= 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

и, внося zl в формулу

(12.18"), получим

три

векторных

равен­

ства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дис

 

ди1

 

 

 

 

 

 

 

 

дх dt

или

девять скалярных

кинематических

условий

на

поверхности

разрыва Н(х, ^)=0:

 

 

 

 

 

 

 

 

 

д

ди'

_ /

dui \

_/

ди>

\ _

Vfe

Г /

ди>

\

_/ Ы

\ j _

 

dxk

\

dxk /

2

\

J l

 

|_ \

д/

/ 2

\

/ 1 J

 

 

 

 

 

IT

 

(f, * = 1 ,

2,

3).

(12.18)

 

 

 

 

 

 

 

 

 

 

 

 

Применим к объему ДVH законы сохранения (11.28) — (11.31), полагая в них t\ = t, t2 = t+dt. Высота цилиндра, представляющего объем ДVH, равна v^dt и предполагается бесконечно малой более высокого порядка сравнительно с линейным размером площади основания Д2я:

v ^ d t= o { y d 2 H), AVH= v^dtA 2H,

и потому интегралами по боковой поверхности в законах сохране­ ния можно пренебречь. На «верхнем» основании цилиндра, име­ ющем положительную нормаль v, все величины имеют индекс «1» и постоянны на всем интервале времени на нижнем (с нормалью —v) также постоянны, имеют индекс «2». Из (11.29) находим

Др-ДКн= Д (pv) n

или

(12.19)

и(г)Др=

пД (pv).

В левой части (11.30) первый интеграл равен Д(ру)ДУя, вто­ рой равен

[PIVI (Vjn)—P2v2 (v2n)] ЛД2Н;

в правой части 53<_v)= —S~2V\ cP(v)=c7:l<1v,, и потому первый ин­

теграл равен —AtP{v)dt Д2Н второй —

порядка

pFdtAVH, т. е.

высшего порядка малости. В результате

из (11.30)

получаем

o(v)A (pv)—p2v2 (v2n) -f p^j, (v1n) = —Д ^(у)

(12.20)

В (11.31) интеграл, содержащий массовую силу, также отбрасы­ вается, интеграл, содержащий параметр р, обозначим Qp-ДЕя^; тогда, аналогично предыдущему, из (11.31) получим

и('')Д Р и + т ) ] - " 4 [р (“ +

v 1 = — A{S3iv)\) -f nAq +

J

(12,21)

Если

величина №р + р^р

ограничена

внутри

объема АУя, то

й р = 0; если же она

терпит

разрыв на

фронте

волны,

то может

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

 

 

 

 

 

t±dt

 

 

 

 

 

“«“

T ib ;

t$

ЬУН( r »+ p '''>)Л,

 

<12'22>

Уравнения — два

скалярных

(12.19),

(12.21)

и одно векторное

(12.20)

— называются динамическими

условиями на

поверхности

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

образований на основании

(12.17)

приводятся к простому

виду

р2 (D—пДу) = pjD;

PjDAv= —ДсР(v);

(12.23)

Pi D

Аи + -i- (Av)21 = —cp2V Av + nAq + Ц

 

Впрочем, (12.23)

получаются из

(12.19) — (12.21), если

в них

заменить

у1-^0,

у2-^Д^,

(12.24)

 

так как физические законы в МСС не зависят от скорости пере­ носного движения Vi окрестности физической точки.

Как видно из (12.15), (12.17), функция #(х, t) удовлетворяет уравнению

—— —V igradtf^D Igradtfl.

(12.25)

dt

 

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

ства среды, т. е. видом функционалов состояния

энергии и

и др.

оператор

Нормальным для многих сред является обратимый

т.е. имеющий взаимно однозначный обратный относительно

&в лагранжевых координатах (§ 9):

Su= & v'{zmn, Т, р). ги =&7}1 (Sm Т, р). (12.26)

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

содержащие две функции

— скалярную р(р, 71 р)

и тензорную

f i j { v m n, Т , р), разрешимые однозначно

относительно

р и v mn (при

^тпФ20) !

 

 

 

 

O u= — Р(р,

Т, P ) 6 i / + / i / ( w mn, т,

Р),

 

Vmn= fmn(aii + P8ii> Т >Р)>

Р = Р(Р-

Т, Р).

Первое из представлений (12.27) функционала У

легко

пере­

писывается

в лагранжевых^ координатах

V g = Ро) и дает

вы­

ражение

Функционал 2F (12.26), вообще говоря,

также

обра­

тим. В этих нормальных случаях функции

независимы

между

собой, функции (12.27) также независимы между собой, т. е.

допускают обратные

oTj/1 ftj1

жидкостей (§ 13)

имеют место

В особом случае

идеальных

тождества

 

 

 

 

ft/ — О (I,

/= 1 , 2, 3),

(12.27')

следовательно, обратные операторы fjj1 не существуют:

аи = pblh

Oi,6,/ = —3р,

 

S*{v)= o i;viel= —pv.

(12.28)

В особом случае идеально пластических изотропных тел суще­

ствует соотношение (§ 18)

 

(о;/—обц){аи-—об,/) = const, З а = о ;/б;/,

(12.29)

и потому оператор SF не вполне обратим.

Для всех сред, принимаемых в теории как объемно-несжимае­

мые,

 

Р=Ро> A = V 7 = 1. div v = 0 , t»mn6mn = 0,

(12.30)

оператор ЯГ не вполне определен, так как вследствие (12.30) не

могут быть заданы любые процессы (8 , 7\ р) для физической частицы, а значит, и на поверхности 2 всей области тела. Интег­ рируя divv = 0 по объему области V', получим

J divvdV= f vvd2 = 0.

(12.31)

V 2

При этом обычно существует обусловленный обратный оператор

(согласованный с условием (12.30) и дающий однозначную зави­ симость 8 или V от (5, 71 р)) или обусловленный прямой ЯГ\\

(о, в, Г, Р), *{^ ,/ = 0, g = 1.

(12.32)

Основной постулат МСС, утверждающий, что внутреннее сос­ тояние малой частицы вполне определяется природой среды, за­ данием закона изменения ее границы или сил на границе и при­ тока энергии (6М, 6'QT 6'ЛР) через границу во времени относит­ ся, конечно, и ко всей области движения среды, если учесть еще работу массовых внешних сил. Это значит, что решение полной системы уравнений МСС при граничных условиях, правильно от­

ражающих приток энергии и импульса непосредственно на грани­

це, существует.

(силы тяжести,

электромаг­

Приток энергии дальнодействия

нитных сил,

) отражается работой

массовой силы

F, входящей

в уравнения движения, и функциями

параметров р,

входящими

в законы сохранения.

 

 

Движение

физической границы 2 g области в лагранжевых ко­

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

jc= <ps (x, /), х = х 2(а, Р),

(12.33)

причем х = хг(а, р) —

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

поверхности,

на которой эти точки

расположены при t= t0t а х=ф2— ее пара­

метрическое уравнение в момент t\ или заданием вектора скорос­

ти точек vr (x, t) на

поверхности

Ф(ж,

/)= 0

эйлерова

простран­

ства

t) = \z ( x ,

I)

 

 

t) = 0,

 

v (лг,

при

Ф(лг,

(12.34)

если эта поверхность известна и является физической

границей

Eg. В этом случае, как уже выяснено,

(12.33)

и есть

решение

уравнения d ^ ld t = vz(ф2, t)y и

условия

(12.33),

(12.34)

эквива­

лентны.

 

 

 

 

 

 

 

Приток тепла определяется заданием на S нормальной состав­

ляющей вектора потока тепла qv = q* или условием на S

 

 

qgrad®—

1grad Ф | = 0 ,

 

 

(12.35)

которое может быть представлено с помощью как (12.33), так и (12.34) . С точки зрения основного постулата условия (12.33) или (12.34) и (12.35) вместе с начальным условием во всей области

/ = /0 Т = Т 0 (12.35')

являются естественными для сред с обратимым функционалом ЗГ> и решение термомеханических задач существует. Условия (12.34) могут быть заданы и на произвольно движущейся заданной по­ верхности Ф(*, t)= 0.

Для обратимых операторов с точки зрения основного постула­

та задания & или 5 эквивалентны, и потому эквивалентны в смысле существования решения задания условий (12.33) или ди­ намических условий

Sl/via /= ^ )2V (х,

t)

при

х = х2(а,

Р)

(12.36)

или

 

 

 

 

 

 

о./??,с/ =

с^у1’ (х,

t)

на

Ф(лг, /) =

0

(12.37)

с сохранением условия

(12.35).

Замена (12.35)

условием

(12.10)

физически также оправдана.

Интересно сопоставление кинематических и динамических ус­

ловий

(12.34),

(12.37) и термодинамического (12.35) с условиями

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

разрыва

(12.23):

при

постоянной

«массе»

piD =

= const

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

среде

(^ =

0, cMv)= 0,

qi = 0,

£2Р= 0)

задание скорости Av тождественно заданию силы AcP(v) и на­ оборот; вместе с заданием теплового потока nAq любое из этих механических условий однозначно определяет энергию Аа.

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

Особые случаи граничных условий для сред с необратимыми

и не вполне обратимыми операторами У уточняются с учетом их свойств. В несжимаемых средах условия (12.34), (12.36) допол­ няются ограничениями (12.31) на задаваемые векторы (ф£ уг) и дополнительным заданием граничного условия на а, вытекающим из (12. 32):

°Hnini + o=ch{n)n = Nn на Ф = 0 .

(12.38)

При условиях (12.34), (12.31) дополнительно должно быть за­ дано нормальное напряжение Nn в одной точке на поверхности. При условиях (12.37) о(х, t) является новой искомой функцией, и ей соответствует уравнение несжимаемости

divv= 0 или g = l

(12.39)

и условие (12.38) на границе, вытекающее из (12.37). В идеаль­ ных жидкостях сообщаемая через границу 2 мощность механиче­ ских сил сводится только к работе нормального давления р на нормальной составляющей вектора скорости. Поэтому условия (12.34) сводятся к заданию одного скаляра v*n:

nv=vxn{x, t) на Ф = 0,

(12.40')

а условия (12.37) на основании (12.28) приводятся к виду (12.38), причем о= —р и задан скаляр Nn:

p = — Nn на Ф = 0 .

(12.40")

В случае идеально пластической среды кинематические усло­ вия (12.33), (12.34) сохраняются и решение замкнутой системы должно существовать в силу основного постулата; по той же при­ чине при условиях (12.37) статическое решение (§ 18) при про­ извольных нагрузках не существует на основании (12.29).

Закончим общими замечаниями о краевой задаче МСС. Каж­ дое из уравнений замкнутой системы и каждое из граничных и начальных условий термомеханической задачи можно рассматри­