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

книги из ГПНТБ / Прессование алюминиевых сплавов. Математическое моделирование и оптимизация

.pdf
Скачиваний:
15
Добавлен:
25.10.2023
Размер:
13.98 Mб
Скачать

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

Главные направления девиатора скоростей деформа­ ций совпадают с главными направлениями тензора.

Характеристическое уравнение имеет вид:

\Ѣк — ^ifcl = 0 или

X3— т]11 А, — г|III

О,

 

 

поскольку его первый инвариант равен нулю:

т1и Л12 _|_

Л22 Лгз _j_ Лзз

Лзі

гІ2і 'Пгг

Лзг Лзз

Ліз

Л и

И ( 5 п - У

2 + ( ^

- У

2 + (5эз-

+ (> (% + & + &)];

 

 

лІП = ЬъІ-

 

 

 

Введем величину

 

 

H = 2 V W \

 

 

(11-24)

которую назовем интенсивностью скоростей деформации сдвига. Эта величина, как и степень деформации сдви­ га Г, определяется зависимостью

і

 

Г = j HÄ,

(ІІ-25)

о

 

где интегрирование по времени t ведется вдоль траекто­ рии материальной частицы, будет часто встречаться в дальнейшем.

Криволинейные координаты

Пусть в области D заданы три однозначные и непрерывно диф­ ференцируемые функции:

а = а (х і , х2, х3), ß = ß (хі, х2, х3), у = у (лц, х2, х3),

удовлетворяющие условию

о,

D (x i, х2, х3)

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

50

ную тройку чисел (ос, ß , у) — криволинейные координаты этой точки. Условие вида:

а (а'і , х 2, х 3) — const

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

а ( * і і х2, Хз) — const, ß ( X i , х2, х3) = const,

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

Каждая точка из D может быть представлена как точка пересе­ чения трех координатных линий или трех координатных поверх­ ностей.

Будем в дальнейшем рассматривать только ортогональные кри­ волинейные системы координат. Элемент длины дуги между двумя «соседними точками» определяется формулой [10]:

ds2 = (d x j2+

(dx2)2+

(dx3)2= H l(d a)*+ Я2

(cfß)2 +

Я 2 (dy)*, (11-26)

где

 

 

 

 

 

 

 

 

 

К

дхі \ г , / ам

2 1

дх3

н\ =

 

+

дх2\ 2

лГ + 1 7

+

да

ар )

І-ІГI +

 

 

 

а р /

+

дхз

Н\-

дхіУ ,

Ідх2у

дх3\*

 

 

a ß )'■

ду +

( l ? j +

ду

 

 

 

 

 

 

 

 

Н- - коэффициенты Ляме.

Скорости, деформации в криволинейной системе координат могут быть найдены по следующим зависимостям [11]:

 

 

dva

1

 

dHa

 

ѴУ

 

дНа

 

 

 

 

ЭР +

 

ду

 

 

da

1 Яа Яр

 

На Н у

 

 

dvß

1

 

"V

Jd h .

ѵа

 

аяр

’ßß:

Но

 

Я«Яр

da

HßHv

ду

 

=ѴѴ

1

dvy

1

 

 

дНу

 

 

дНу

ң

dy

Ha Hy

да

 

ЯрЯѵ

ар

 

 

 

 

 

 

 

 

/»Э

 

 

 

а

 

(11-27)

=aß :

2

\H ß

д

 

 

 

( ѵ а \-\

К

да

І^Яр

У

к

K

/ J

 

 

 

 

 

=ßv:

J _

^

4

 

1 - Ч

+

Яр

_a_

 

EL

2

 

я ѵ aY

H ß

 

[Яр

ap

Я„

 

 

 

J _

 

 

 

 

 

Л

 

і

 

 

=VV '

2

H y dy

 

 

 

Яа

да

 

HV / J

 

 

 

 

 

4

51

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

Рассмотрим частный случай плоскопараллельного движения сплошной среды, когда поле скоростей имеет следующий вид:

Ѵ1 = ѵі (*1, *2),

ѵ2 = Щ(Д, х2), ѵ3 - 0.

При этом все линии тока параллельны фиксирован­

ной плоскости

(Х\, х2).

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

 

Если среда несжимаема, компоненты скорости удов­

летворяют уравнению несжимаемости:

dvj

dv2_ _ Q

дхі

дх2

Введем в рассмотрение функцию тока -ф (агі, х2); ком­ поненты вектора скорости определим зависимостями:

і>2 = —

(11-28)

 

0*1

Легко видеть, что при любом выборе функции тока условия несжимаемости удовлетворяются тождественно.

Подставим зависимости (П-28) в дифференциальное уравнение линии тока:

dx\ _ dx2 dty гД

дх2 дху

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

= 0

или

d\|>(хи х2) = 0.

Следовательно, вдоль линии тока выполняется ус­ ловие

ф(дГі, лг3) = const,

т.е. функция тока сохраняет постоянное значение.

52

Возьмем две точки Р\ и Р2 в области D; пусть apt и фг — соответствующие значения функции тока.

Подсчитаем поток вектора скорости через дугу С, со­ единяющую эти точки. Интеграл по поверхности (11-12) от нормальной составляющей скорости для двумерного движения сводится к интегралу по контуру С:

Я =

j vn ds.

 

 

 

 

 

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

 

к

dxi

 

 

 

 

1 Г

 

 

 

 

 

 

поэтому

йф dxi

dxz

Ѵп = »A + щк

дХі

ds

дхг

ds

 

 

Следовательно, поток равен

 

 

q =

j dty = ф2 — Фі-

 

 

(11-29) .

 

с

 

 

 

 

 

 

 

 

4. НАПРЯЖЕННОЕ СОСТОЯНИЕ

Внешние силы

Пусть нам дано некоторое тело D, ограниченное по­ верхностью 5.

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

К основной особенности этих сил следует отнести то, что они непрерывно распределены по объему или поверх­ ности тела.

Плотность

Возьмем внутри рассматриваемого тела произволь­ ную точку М и выделим в окрестности этой точки эле­ ментарный объем Аw с массой Ат.

53

Составим отношение Ат/Aw и устремим размеры эле­ мента к нулю. Величина

р(М)— lim

дю - о А ш

называется плотностью среды в точке М.

Массовые внешние силы

Пусть равнодействующая внешних сил, действующих

на выделенный элемент, равна ДФ.

Составим отношение ДФ/Ат и вычислим предел

р

ДФ

1

ДФ

F = п т — = — п т — .

Д т -> 0

A ffi

р Д ш - і- о Д а У

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

Поверхностные внешние силы

Как мы уже отмечали, на поверхности 5 действуют внешние поверхностные силы. Выделим элемент поверх­

ности с единичной нормалью п. На этот элемент дейст­

вуют силы, их равнодействующая равна АР. Устремим размеры элемента к нулю и вычислим предел

ап— lim — .

Ä S -I-OA S

Векторная величина ап (по индексу п суммирования не производится) характеризует интенсивность прило­ женных к поверхности внешних сил и называется векто­ ром поверхностных напряжений.

Внутренние напряжения

Вернемся к первоначально рассматриваемому телу D и мысленно рассечем его поверхностью 2 на две части D\ и D2 (рис. 5).

Отбросим часть D2. Оставшаяся часть тела будет оставаться в равновесии, если на поверхности 2 действу­ ет некоторая система сил.- Назовем их внутренними си­ лами.

54

Возьмем на поверхности 2 точку М и выделим в ее

окрестности элементарную площадку Д2 с нормалью п. Пусть равнодействующая сил, приложенных к пло­

щадке Д2, равна АР.

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

.. Д Р

а" = lim —

Д 2^0 Д2

действующих в точке М на площадку с нор­

малью п.

осо­

 

Существенной

 

бенностью здесь

явля­

 

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

 

ентация поверхности 2

Рис. 5. Тело под действием внешних

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

нор­

сил

мали П.

Выберем в связи с этим три площадки, проходящие

через точку М и нормальные ортам е».

—У

На площадке, перпендикулярной орту ей действует

вектор напряжений а 1с компонентами

fall I СТ12> °1 з)-

Аналогично на двух других площадках действуют векторы напряжений

° 2 — (°21> ^22> ®2.з) И CF3 — f a e i , СГ32, СТ33).

Составим матрицу

°12 <Аз

° 2 1

° 2 2

° 2 3

° 3 1

&32

° 3 3

и рассмотрим физический смысл ее компонентов. Диагональные элементы матрицы сгц, а22, сгзз назы­

ваются нормальными компонентами напряжений, по­ скольку они представляют собой проекции векторов на­ пряжений на нормаль к площадкам.

Боковые элементы матрицы называются касательны­ ми компонентами напряжений. Они являются проекция­ ми векторов на плоскости площадок.

55

Напряжения на наклонной площадке

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

произвольной площадке.

Проведем

через точку

О площадку

с нормалью

п и построим элементар­ ный тетраэдр, образован­ ный пересечением этой площадки с плоскостями, перпендикулярными к ко­ ординатным осям (рис. 6).

Устремим размеры тетраэдра к нулю, при этом он будет стягивать­

Рис. 6. Элементарный тетраэдр ся в точку 0.

В результате получа­ ем формулу Коши

о” =

О1/(•

 

(ІІ-ЗО)

или

в координатной записи

 

^itk

®ik4*

 

 

Вспоминая правило умножения матриц, найдем (И-31)

 

/ ° 1 1

° 1 2

° 1 3

K l ° я2 ° лз) — ( 4 4 4 ) 1 °2 1

° 2 2

° 2 3

 

 

^ 3 2

с*зз

Таким образом, формула (ІІ-ЗО) каждому вектору п

ставит в соответствие вектор напряжений о", действую­

щий на площадке, ортогональной к нормали п, по пра­ вилу умножения вектора на матрицу напряжений. Матрица (ста) является тензором

<*11 <*12 <*13

<*21 <*22 <*23

<*31 <*32 <*33

который называется тензором напряжений,

56

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

Итак, мы ввели в рассмотрение новый тензор второй валентности — тензор напряжений ||a,ftli.

В дальнейшем будет показано (см. с. 62—63), что этот тензор симметричен.

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

а Ч 2 + о 11 Я — стш = 0.

(11-32)

к 3 —

Корни

>з

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

В новой координатной системе тензор напряжений принимает вид

 

Ol

0

0

т„

0

о2

0

 

0

0

Оз

При этом касательные компоненты напряжений об­ ращаются в нуль.

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

Коэффициенты характеристического уравнения обра­ зуют систему инвариантов

0 = а 11 +

<*22 +

<?33 —

+

<*2 +

°з\

 

 

Оц

012

а22 ^23

 

° 3 3

°\31

CTl 0*2 ~Ь

dg [ OsOp,

CT21

0-2І +

0^32 0*3

+

а 13

<Ді

III I ,

°= К-fei -

Девиатор напряжений

Тензор напряжений можно представить в виде суммы девиатора и шарового тензора:

To = Do+ % !•

Здесь величина

СТср — ~ (<Тц + ст22+ °3з)

(ІІ-ЗЗ)

57

называется средним напряжением.

 

Компоненты девиатора равны

 

^ik

' ^ik

^cp

 

 

 

 

 

т.

e.

 

 

 

 

 

 

 

Ой

er,,,

a12

 

or13

 

 

 

° 2 1

° 2 2

°cp

° 2 3

 

 

 

<%

° 3 2

a 3 3

°cp

 

 

Характеристическому

уравнению

девиатора напря­

жений

 

 

 

 

 

 

Я3 + 5” Я - 5 іп = 0

 

 

 

 

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

 

 

S11

S 12

S22

S23

S 33

S 31

7 - \(°П — 022)2-f

 

 

+

 

+

 

 

 

S 21

S 22

S32

S 33

S13

S11

О

 

 

 

 

 

 

 

ь°зі)];

 

Величина

 

 

 

 

 

Т = V И

 

 

 

 

 

(П*34)

называется интенсивностью касательных напряжений.

5. ОСНОВНЫЕ УРАВНЕНИЯ МЕХАНИКИ СПЛОШНЫХ СРЕД

Уравнения неразрывности

Выведем уравнение неразрывности сплошной среды. Мысленно вырежем в последней объем W, ограниченный материальной поверхностью 2 (рис. 7).

Масса, заключенная внутри этой поверхности, оста­ ется неизменной во времени и равна

т= рdW = const,

отсюда

pm _ Q dt

58

или

(Н-35)

w

Вспоминая формулу (II-18), определяющую полную производную от интеграла во времени, получаем

Рис. 7. Материальный объем W

 

где vn— ViU — нормальная составляющая

скорости на

поверхности 2; k — компоненты нормали

п.

Воспользуемся формулой Гаусса — Остроградского и преобразуем интеграл по поверхности 2 в интеграл по объему:

W

Отсюда уравнение (II-38) принимает вид:

Из леммы (см. с. 28), поскольку объем произволен, следует, что

(II-36)

59

Соседние файлы в папке книги из ГПНТБ