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

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

..pdf
Скачиваний:
6
Добавлен:
12.11.2023
Размер:
9.02 Mб
Скачать

П3. МЕТОДЫ ПОСТРОЕНИЯ КООРДИНАТНЫХ ФУНКЦИЙ

 

1

 

 

h1

1

 

h1

 

 

2

 

нейной полосе высотой h

 

 

E

≤ +

 

в направлении оси Е

 

со скоростью

 

 

 

 

 

 

 

2

 

 

2

 

 

 

Vf

=

V0 h0

. Необходимо в области =

3

j построить непрерывное соленои-

 

 

 

h1

j=1

дальное векторное поле v = vi ei с непрерывными компонентами его градиента, удовлетворяющее следующим граничным условиям: в Ω1 компоненты вектора

скорости v1 = 0; v2 = V0, в Ω2 нормальная к поверхностям E1 = h2 и E1 = − h2 со-

ставляющая вектора равна нулю, в Ω3 компоненты вектора v1 = 0; v2 = Vf . Ранее отмечалось, что соленоидальное векторное поле (П3.25) удобно стро-

ить с помощью функции тока (П3.27). Область движения потока удовлетворительно описывается функцией тока

 

h0 E1

 

(П3.54)

ψ = −V0 h

.

 

Действительно, в подобласти Ω1 имеем h = h0 и ψ = –V0Е1, что после дифференцирования по Е1 и Е2 по формулам (П3.27) обеспечивает выполнение граничных

 

 

 

 

h

 

 

условий v1

= 0; v2

= V0. При этом на верхней границе E1

= −

0

 

ψ = ψ+ = 0,5V0h0 =

2

 

 

 

 

 

 

 

 

h

 

 

= const; на нижней границе E1

=

0

 

ψ = ψ= –0,5V0h0 = const; на оси симмет-

2

 

 

 

 

рии (Е1 = 0) ψ = 0. Легко убедиться в том, что функция тока (П3.54) является непрерывной во всей области Ω и сохраняет постоянные значения на соответствующих линиях тока. Например, для верхней и нижней границ подобласти

Ω

достаточно в (П3.54) подставить E = −

h

и E =

h

, что приводит к значени-

 

 

2

1

2

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

ям ψ = ψ+ и ψ = ψсоответственно. В подобласти Ω , где h = h

1

и ψ = −V

h0 E1

,

h

 

 

 

 

3

 

0

 

 

 

 

 

 

 

 

 

 

1

 

также выполняются условия непрерывности функции тока (ψ = ψ+ при E1 = − h21

581

ПРИЛОЖЕНИЕ (МАТЕМАТИЧЕСКИЕ ОСНОВЫ МЕХАНИКИ СПЛОШНЫХ СРЕД)

и ψ = ψпри E1 = h21 ) и граничные условия для вектора v , так как после под-

становки последнего значения функции тока в (П3.54) получаем v1 = 0; v2 = Vf. По условию задачи текущая высота полосы h = h(Ek) в подобласти Ω2 на уча-

стке д (П3.53) должна обеспечивать конфигурацию верхней и нижней границ

в виде окружности радиуса R. Если уравнение окружности записать в эйлеровых координатах Ek:

h = h1 + 2 R2 R2 E22 ,

(П3.55)

где д E2 0 , то нетрудно убедиться в том, что непрерывная в Ω функция тока (П3.54) после подстановки ее в (П3.27) с учетом (П3.55) приводит к разрыву векторного поля v и к скачкообразному изменению компонент его градиента на уровнях E2 = −д и Е2 = 0.

Для склейки разрывного в Ω векторного поля воспользуемся методом переходных зон. Выделим в области Ωпять участков: I (EЕ2 Ен); II (Ен Е2 Ен1); III (Ен1 Е2 Еf1); IV (Еf1 Е2 Еf); V (Еf Е2 E+), где Ен = – д – δ1; Ен1 = – д + δ1; Еf1 = –δ2; Еf = δ2. Введем вспомогательное множество координат xk такое, чтобы на участках I, III, V оно совпадало с основным множеством координат

xk = Ek.

(П3.56)

В переходных зонах II и IV протяженностью 2δ1 и 2δ2 соответственно координату х1 оставим без изменения, а координату х2 представим в виде склеивающей функции

 

 

 

 

 

 

*

 

1

 

*

 

 

x1 = E1 ; x2 = 0,5(E2 + a) bg

9cos E

+

 

 

cos3E

,

(П3.57)

3

 

 

 

 

 

 

 

 

 

 

 

 

где для участка II: a = Ен1; b = 0,625; g =

2(Eн1 Eн )

;

E* =

(2 E2 EнEн1 ) ; для

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

g

 

участка IV: a = Е ; b = –0,625;

g =

2(E f 1 E f )

;

E* = (2 E2 E f E f 1 )

. Теперь

 

f1

 

 

π

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

по формулам (П3.27) с помощью (П3.51), учитывая (П3.56) и (П3.57), можно перейти к вычислению компонент вектора v

582

П3. МЕТОДЫ ПОСТРОЕНИЯ КООРДИНАТНЫХ ФУНКЦИЙ

v1 =V0

h0 E1

hf ; v2 =V0

h0

,

h2

h

где текущая высота (П3.55) принимает вид

h = h + 2 R2 R2

x2

;

1

2

 

первая производная этой величины по х2

(П3.58)

(П3.59)

h′ =

2 x2

.

(П3.60)

R2 x22

В переходных зонах (участки II и IV)

 

x2

*

*

 

 

f =

 

= 0,5 + b(9sin E

+ sin3E

).

(П3.61)

E

 

2

 

 

 

 

На участках I и V имеем f = 0, а на участке III – f = 1.

По формуле (П1.78), используя (П3.58), найдем компоненты градиента

( v ) вектора скорости v :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v1

= −

v2

=V0

h0

hf ;

 

 

 

 

 

 

 

 

 

E1

 

E2

h2

 

 

 

 

 

v

 

 

h E

 

 

 

 

2h2 f 2

 

v

 

 

1

=V

0 1

h′′ f + hf ′−

 

 

 

 

 

;

2

= 0,

 

h2

 

 

h

E

 

E

2

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

где

 

2 x

 

b

(9cos E* + 3cos3E* );

f ′ =

E22

=

 

g

 

2

 

 

 

вторая производная h по х2

 

 

h′′ =

2 R2

 

.

(R2 x22 )3

(П3.62)

(П3.63)

(П3.64)

583

ПРИЛОЖЕНИЕ (МАТЕМАТИЧЕСКИЕ ОСНОВЫ МЕХАНИКИ СПЛОШНЫХ СРЕД)

Упражнение П3.10. Показать, что вспомогательные координаты (П3.56), (П3.57) склеивают функцию тока (П3.54) на стыке участков I…V с четвертым порядком надежности

Отметим, что локальная замена координат в переходных зонах естественно приводит к некоторому изменению их границ, а значит, и к изменению всей области исследования. При решении задач МСС, в частности ОМД, метод пе! реходных зон становится эффективным, когда в первом приближении граница области исследования аппроксимируется непрерывной линией или поверхнос! тью с изломами, в которых появляются либо разрывы, либо особенности той или иной природы и того или иного порядка. Из практики известно, что в при! роде такие изломы обходятся потоками реальных материалов, и с этой точки зрения применение переходных зон в окрестности изломов является оправдан! ным. Добавим, что протяженность переходных зон (в рассмотренном приме! ре – величины 2δ1 и 2δ2) в вариационной постановке задачи при необходимос! ти можно рассматривать как варьируемые параметры.

Все вышеприведенные методы пока были связаны с локальной аппроксима! цией искомой функции в подобластях Ωj и склейкой этих аппроксимаций на границах подобластей одномерных или двухмерных областей исследования. Применение метода разделения переменных (П3.4) позволяет задачу о построе! нии одной функции в N!мерном пространстве свести к построению N функций в одномерных пространствах.

Другим способом построения функций в многомерном пространстве явля! ется метод конечных элементов (МКЭ). Суть его состоит в том, что область ис! следования Ω разбивается на «конечные элементы», т. е. на конечное количество подобластей Ωj без разрывов и пересечений, так, чтобы объединение подобластей Ωj образовывало Ω. С этой точки зрения все рассмотренные ранее методы ло! кальной аппроксимации относятся к МКЭ в одномерных областях. Для много! мерных пространств в качестве подобластей используют симплексы (многогран! ники), в вершинах которых вид локальных аппроксимаций определяется свя! зями, накладываемыми на искомую функцию.

Контрольные вопросы

1.Каким требованиям должны удовлетворять основное решение и его кор! ректировка?

2.В чем суть метода разделения переменных?

3.Как используются R!функции для построения корректирующих функций?

4.Какие комплексные функции называются аналитическими?

5.Что называется комплексным потенциалом?

584

П3. МЕТОДЫ ПОСТРОЕНИЯ КООРДИНАТНЫХ ФУНКЦИЙ

6.Какие поля скоростей называются гармоническими?

7.Что такое функция тока, консервативная функция и в каком соотношении находятся эти функции для гармонических полей скоростей?

8.Что такое конформное отображение и как это отображение выполняется?

9.Перечислите свойства интеграла К. Шварца–Э. Кристоффеля. Что означает нормировка этого интеграла?

10.Какие гармонические течения называются простейшими?

11.Чему равен вектор скорости в точке бифуркации?

12.В чем суть метода суперпозиции гармонических течений?

13.Что такое порядок надежности склейки двух функций?

14.В чем суть метода конечных элементов?

585

РЕКОМЕНДАТЕЛЬНЫЙ БИБЛИОГРАФИЧЕСКИЙ СПИСОК

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

Аркулис Г. Э., Дорогобид В. Г. Теория пластичности. – М.: Металлургия, 1987. – 352 с.

Бабич В. К., Лукашкин Н. Д., Морозов А. С. и др. Основы металлургического производства. Учебник. – М.: Металлургия, 2000. – 240 с.

Бердичевский В. Л. Вариационные принципы механики сплошной среды. – М.: Наука, 1983. – 448 с.

Богатов А. А., Мижирицкий О. И., Смирнов С. В. Ресурс пластичности при обработке давлением. – М.: Металлургия, 1984. – 144 с.

Бэкофен В. Процессы деформации. – М.: Металлургия, 1977. – 288 с.

Гантмахер Ф. Р. Теория матриц. – М.: Наука, 1967. – 576 с.

Гольденблат И. И. Некоторые вопросы механики деформируемых сред. – М.: Гостехиздат, 1955. – 272 с.

Грин А. Е., Адкинс Дж. Большие упругие деформации – М.: Мир, 1965. – 456 с.

Гун Г. Я. Теоретические основы обработки металлов давлением. М.: Металлургия, 1980. – 456 с.

Жермен П. Курс механики сплошных сред. – М.: Высшая школа, 1983. – 399 с.

Ильюшин А. А. Механика сплошной среды. Изд. 2. – М.: МГУ, 1978. – 228 с.

Качанов Л. М. Основы теории пластичности. – М.: Наука, 1969 – 420 с.

Колмогоров В. Л. Механика обработки металлов давлением. – М.: Металлургия, 1986. – 688 с.

Коппенфельс В., Штальман Ф. Практика конформных отображений. – М.: ИЛ, 1963. – 406 с.

Кристенсен Р. Введение в механику композитов. – М.: Мир, 1982. – 334.

Кучеряев Б. В. Потапков Н. А. Механика сплошных сред. Раздел: Математические методы решения задач ОМД. – М.: МИСиС, 1992. – 164 с.

Лаврентьев М. А., Шабат Б. В. Методы теории функций комплексного переменного. – М.: Наука, 1965. – 716 с.

Ломакин В. А. Теория упругости неоднородных тел. – М.: МГУ, 1976. – 386 с.

Лурье А. И. Теория упругости. – М.: Наука, 1970. – 940 с.

Лыков А. В. Теория теплопроводности. – М.: Высшая школа, 1967. – 480 с.

Малинин Н. Н. Прикладная теория пластичности и ползучести. – М.: Машиностроение, 1968. – 400 с.

Мейз Дж. Теория и задачи механики сплошных сред. – М.: Мир, 1974. – 320 с.

Мусхелишвили Н. И. Курс аналитической геометрии. – М.: Высшая школа, 1967. – 656 с.

Мусхелишвили Н. И. Некоторые основные задачи математической теории упругости. – М.: Наука, 1966. – 708 с.

Победря Б. Е. Механика композиционных материалов. – М.: МГУ, 1984. – 336 с.

Прагер В. Введение в механику сплошных сред. – М.: ИЛ, 1963. – 311 с.

586

РЕКОМЕНДАТЕЛЬНЫЙ БИБЛИОГРАФИЧЕСКИЙ СПИСОК

Работнов Ю. Н. Механика деформируемого твердого тела. – М.: Наука, 1979. – 744 с.

Рвачев В. Л. Методы алгебры логики в математической физике. – Киев: Наукова думка, 1974. – 258 с.

Седов Л. И. Механика сплошной среды. Т. 1, 2. – М.: Наука, 1976.

Степанский Л. Г. Расчеты процессов обработки металлов давлением. – М.: Машиностроение, 1977. – 424 с.

Стечкин С. Б., Субботин Ю. Н. Сплайны в вычислительной математике. – М.: Наука, 1976. – 248 с.

Томсен Э., Янг Ч., Кобаяши Ш. Механика пластических деформаций при обработке металлов. – М.: Машиностроение, 1969. – 503 с.

Фавар Ж. Курс локальной дифференциальной геометрии. – М.: ИЛ, 1969. – 559 с.

Эльсгольц Л.Э. Дифференциальные уравнения и вариационное исчисление. – М.: Наука, 1965. – 424 с.

587

ПРЕДМЕТНЫЙ УКАЗАТЕЛЬ

Предметный указатель

Аксиома памяти 221

– размерности 220 Аксиомы Нолла У. 219

Альтернативная (кососимметричная, антисимметричная) часть тензора ранга n 331

Альтернирование тензора ранга n 331 Антисимметричная (альтернативная, кососимметричная) часть тензора ранга n 460

Анизотропия 18 Армирование 120

Базис взаимный 487

– пространства 533

– ортонормированный 536 Баланс мощности 194,196

работы 197

тепла 202

Безмоментная теория напряжений 211 Брахистохрона 361

Валентность тензора 448 Вариационное исчисление 361 Вариация функции 539

функционала 539 Варьирование функции 539 Вектор единичный 443

– деформации 220

искажения (дисторции) 28, 85

нормированный 443

нулевой 443

перемещения 27

скорости перемещения 32

теплового потока 202 Векторная линия 66

Вилка оценки эффектиных модулей 280

Вихреисточник 572

Вихрь (ротор) порядка n 487

полный 487

тензора ранга m (левый, правый) 486 Взаимно однозначное соответствие 18, 562 Взаимная ортогональность векторов 444

Воздействие внешнее 19

термомеханическое 22

Волокно 24

материальное 26

направленное 24 Вычитание тензоров 455

Вязкопластичный потенциал i-го композита 279

– эффективного модуля 279

Гармоническое течение 568 Гетерогенность (неоднородность) 16 Гидростатическое давление 153 Гипотеза единой кривой 206 Главная диагональ матрицы 446 Главные значения тензора 463 Голономные связи 553 Гомеоморфизм 17

Гомогенность (однородность) тела 16

макроскопическая 22

Градиент деформации материальный 27

– пространственный 28

порядка n тензора ранга m (левый, правый) 485 Граница внешнего воздействия 20

многосвязная 20

односвязная 20

свободная 20

тела 19

588

ПРЕДМЕТНЫЙ УКАЗАТЕЛЬ

Граничная точка 19 Группы преобразований симметрии 211

Давление гидростатическое 153 Движение нестационарное 22

поступательное 24

стационарное 22

термомеханическое 18 Движения описание 31

– лагранжево (материальное) 31

– эйлерово (пространственное) 32 Двойникование 271 Двойственная вариационная задача 555 Девиатор 460

деформации 56

– направляющий 220

скоростей деформаций 88

напряжений 153

Действительное (реальное) поле (Р-поле) скоростей 75 Деформация главная 55

горячая 257

малая 52

– линейная 53

– сдвиговая 53

конечная 34

– линейная тензора Грина А. Е. 34

Коши О. 42

Лагранжа Ж. 36

Эйлера Л. 44

– угловая (сдвиговая) тензора

 

 

Грина А. Е. 35

Коши О. 43

Лагранжа Ж. 40

Эйлера Л. 44

объемная 60

окрестности 25

остаточная 254

относительная 34

пластическая металла 202, 254

средняя 55

холодная 257

упругая 202 Диагонализация тензора 463 Диада 457

Дивергенция порядка n полная 486

тензора ранга m (левая, правая) 485 R-дизъюнкция 532

Диполь 572 Дисторция (искажение) окрестности частицы 27 Длина вектора 442

Единственность решения задач 532

Жесткое вращение 53, 56 Жонглирование индексами 33

Закон Гука Р. обобщенный 208

– для изотропных сред 169

движения 26

парности касательных напряжений 192

преобразования компонент вектора 447

сохранения вещества 26

– тепловой энергии (баланс тепла) 202

текучести ассоциированный 266

трения Зибеля Э. 334

– Кармана Т. 334

– Кулона Г.–Амантона С. А. 335

Фурье Ж. 202

Звезда 20

Идеальная жидкость 563 Изменение объема 56

формы 56 Изопериметрические задачи 553

ограничения 553

589

ПРЕДМЕТНЫЙ УКАЗАТЕЛЬ

Изотропия 18 Инвариант 447

тензора (первый, второй, третий) 462

Инверсия 211 Индексация Миллера Х. 164

Интеграл Шварца К.–Кристоффе- ля Э. 569 Интенсивность касательных напряжений 153

сдвиговых деформаций 57

– скоростей деформаций 90 Искажение (дисторция) окрестности частицы 24 Испытание с постоянной скоростью деформации 82

– – – – деформирования 82

стандартное 253

Кинематика 24 Кинематически возможное поле (КВполе) скоростей 75

– разрывное 75 Комплексный аргумент 562 Композит 18 Компоненты матрицы боковые 446

– диагональные 446

вектора ковариантные 488

– контрвариантные 487

метрического тензора 488 Конвективные слагаемые 27 Конвективный перенос тепла 204 Конечные деформации 33, 41

– линейные тензора Грина А. Е. 34

Коши О. 42

Лагранжа Ж. 36

Эйлера Л. 44

– угловые (сдвиговые) тензора Грина А.Е. 35

– – – – Коши О. 43

– – – – – Лагранжа Ж. 36

– – – – Эйлера Л. 44 Константа интеграла Шварца К.–Кри- стоффеля Э. 569

упругая Ламе Г. 217

Континуум 18 Контур текучести 266

Конформное отображение 568 R-конъюнкция 560 Координаты материальные 26

обобщенные 52

пространственные 27 Корректировка основного решения 74, 559 Корректность математической постановки задач 532

Кососимметричная (альтернативная, антисимметричная) часть тензора 331 Коэффициент вытяжки 113, 353

вязкой податливости 209

вязкости 208

– объемной 241

изменения длины волокна 32

податливости 208

Пуассона С. Д. 252

температуропроводности 203

теплопроводности 203

трения по напряжению пластического сдвига 334

– – нормальному давлению 334 Коэффициенты Фурье Ж. 536 Краевая задача ТП 24

– Коши О. 338

– Римана Г. 339

– смешанная 339

Краевой задачи постановка 22,239

кинематическая 239

математическая 22

статическая 343

– решение 22

590