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

Методическое пособие 382

.pdf
Скачиваний:
4
Добавлен:
30.04.2022
Размер:
1.09 Mб
Скачать

контактной поверхности, при деформации пренебрегаем бочкообразностью.

Пусть заготовка под действием внешней силы P

получила деформацию на величину h (рис.16). Примем, что

перемещения в заготовке распределены линейно вдоль оси

z .

Тогда для произвольного сечения справедливо:

 

 

 

 

 

 

u

z

z h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Очевидно, что граничные кинематические условия по

перемещениям выполняются.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

 

 

k

 

 

h

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z=h

 

 

 

uz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

Рис.16.

Напряженное состояние при осадке цилиндрической заготовки – осесимметричное. Компонентами тензора

деформаций являются величины , , z , , z , z Величину осевой деформации можно определить

непосредственно:

z uz h

z h

Остальные деформации определим из условия несжимаемости:

81

 

 

 

 

 

 

0;

 

u

 

u

 

u

z

0

 

 

z

 

 

z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

последнюю формулу можно преобразовать к виду:

 

 

 

 

1

 

 

u

 

 

h 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрируя, получим:

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

1

h 2

f z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 h

 

,

 

 

 

 

 

 

при

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0,u 0 f z

0

 

 

 

 

 

 

 

 

 

 

 

Окончательно:

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

1 h

 

 

1 h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 h

2 h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, тогда

 

 

 

 

 

 

 

 

 

Сдвиговые деформации в площадках отсутствуют,

поскольку

 

 

напряженное состояние

 

осесимметричное,

следовательно:

 

 

 

 

z

0

. С

другой

стороны

 

 

 

 

 

 

 

 

 

u

 

 

u

z

0

 

 

 

 

 

 

 

 

 

z

 

z

 

 

 

 

 

 

 

 

 

 

,

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

Таким образом,

деформации

z

- главные.

Попутно

заметим,

 

что,

поскольку

все

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

однородное.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интенсивность

 

 

деформаций определим следующим

образом:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

i

 

 

 

 

2

z 2 z 2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3 h 2

3 h

 

2

 

h

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

2 h

2 h

 

 

 

h

82

Составляющие уравнения баланса работ11:

Ap P h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

A

s

 

 

i

dV

s

 

 

2 d dz

s

 

 

h

2 d dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

0 0

 

 

 

 

 

 

 

 

 

0 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h 2

2

 

 

 

z h

 

d

2

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2 d k

h 2

 

2 k udf 2 k u

 

h

 

2 2 d

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

2

3

2

 

 

 

 

d 3 h

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

3

 

 

 

 

 

 

 

 

 

 

 

12

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P h

 

d 2

 

 

 

 

 

 

 

 

 

 

 

d 3

 

h

 

 

 

 

 

 

 

 

 

 

s

 

 

4

h

k

12

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

приняв

 

 

трение

 

 

по

 

 

Прандлю

 

k s s ,

окончательно

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2

 

 

1

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

P s

 

 

 

 

 

 

1

 

 

s

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h

 

-

 

уже

 

известная

нам формула

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Зибеля.

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

11 Коэффициент 2 в формуле для A присутствует т.к. существует две поверхности, трение на которых одинаково.

83

Лекция 8. Использование вариационных методов при решении задач проектирования технологических процессов.

Еще одним методом, относящихся к группе энергетических, является вариационный метод. Применение вариационного метода для анализа задач обработки давлением связано с уральской школой: И.Я.Тарновский, А.А.Поздеев, В.Л.Колмогоров и др. Прежде чем перейти к изложению метода остановимся на некоторых элементах вариационного исчисления, которые будем использовать в дальнейшем.

Понятие функционала. Основная задача вариационного исчисления.

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

Функционалом называется переменная величина, зависящая от одной или нескольких функций.12 Функционал задан, если каждой функции из некоторого класса функций поставлено в соответствие некоторое число.

Одним из простейших функционалов является длина кривой, соединяющая две точки (рис.17).

Пусть y f (x) - произвольная кривая, соединяющая две точки:

x a, y A

x b, y B

12 Функция, в отличие от функционала зависит от переменных. Если отвлечься от строго математического определения функционал – это функция, от одной или нескольких функций.

84

 

 

y=f(x)

 

 

 

Таким

 

y

 

образом,

 

класс

 

 

 

 

 

 

 

кривых

 

определен

 

 

 

 

 

B

 

 

 

тем,

что

они

 

l

 

 

проходят

 

через

две

 

 

 

 

заданные

точки.

A

 

 

 

Наложим

 

на

 

 

 

 

 

функции,

 

 

 

 

 

 

 

 

 

 

определяющие

 

 

 

 

 

 

кривые,

 

 

 

x

дополнительное

 

 

 

ограничение, что они

 

a

b

 

 

 

непрерывны

и

 

 

 

 

Рис.17.

дважды

дифференцируемы.

 

Длина кривой

 

l

 

 

 

dl

 

l

 

 

dx2 dy2

b

a

dy

2

b

 

 

2

 

 

 

 

1

 

 

 

dx

1 y' dx

 

 

dx

 

a

 

 

 

 

 

 

 

 

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

Как видно из этих примеров, функционалы могут быть определены кратными интегралами.

Когда функция y f (x) , от которой зависит функционал, непрерывна и дважды дифференцируема, функционал может быть записан в следующем общем виде:

J J y x abF x, y, y' dx

Всегда, когда функция f (x) принимает определенный вид (соответствует определенному классу), функционал принимает то или иное значение – число, характеризующее

данную функцию f (x) .

85

В

примере

с

длиной

кривой

J l ,

F(x, y, y') 1 y' 2 :

J ab1 y' 2 dx

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

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

Вариационный метод как частный случай энергетического метода.

Теперь можно ответить на вопрос, каким образом используется вариационное исчисление в теории решения задач при проектирование процессов обработки давлением.

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

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

86

С одним из таких функционалов мы уже встречались ранее – это функционал полной мощности кинематически возможного поля скоростей идеально жестко-пластичного тела.

J s i*dV k v*df pi vi*dF

V

f

Fp

Это выражение полностью удовлетворяет определению

функционала. Функцией,

от

которой зависит функционал

полной мощности, является функция изменения скорости v по объему деформируемого тела (интенсивность скоростей

деформации i может быть определена по известному полю скоростей), независимой переменной – координаты материальных точек. Класс функций также переделен – они должны удовлетворять граничным условиям и условиям непрерывности – т.е. быть кинематически возможными.

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

Таким образом, задача отыскания поля скоростей является вариационной задачей и может быть сформулирована следующим образом:

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

Метод решения задач математической физики, который состоит в том, что интегрирование дифференциальных уравнений при заданных краевых условиях заменяется отысканием функции реализующей минимум

87

соответствующего функционала называется энергетическим методом.

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

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

при которых функционал полной мощности стационарен.

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

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

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

0.5 ij ij dV piui dF

 

 

V

FP

Здесь Fp - это часть поверхности деформируемого тела

объемом V

,

на котором заданы удельные внешние силы pi

i x, y, z

;

ui

- кинематически возможное поле

88

перемещений, ij - поле деформаций Коши, соответствующее

кинематически возможному полю перемещений; ij - поле напряжений, удовлетворяющее обобщенному закону Гука. Этот принцип используют для вариационной формулировки метода конечных элементов.

Используют также функционал Маркова для жесткопластического тела (с учетом упрочнения), определяющий полную мощность для непрерывного поля скоростей13:

J s i dV pi vi dF

V FP

Здесь vi - кинематически возможное поле скоростей,

ij - поле скоростей деформаций, соответствующее

кинематически возможному полю скоростей; ij - поле напряжений, удовлетворяющее физическим уравнениям связи напряженного и деформированного состояний для кинематически возможного поля скоростей.

Вариация функционала и ее свойства Рассмотрим вариационную задачу отыскания

экстремума функционала

J y x abF x, y x , y' x dx

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

y(a) A, y(b) B .

Напомним, что экстремалью называют функцию, сообщающую экстремум функционалу. Функции, из которых выбирается экстремаль, называются функциями сравнения.

13 В зарубежной литературе часто вместо терминов «полная мощность» и «полная энергия» используют термины «дополнительная энергия (работа)» и «дополнительная мощность».

89

Пусть:

 

 

 

 

 

 

 

 

 

 

 

 

f (x)

- экстремаль функционала J ,

 

 

 

 

 

y

 

 

 

f(x)

 

 

 

 

f (x)

- функция

B

y

 

 

 

 

 

 

 

 

 

сравнения,

 

 

 

 

 

 

 

бесконечно

близкая

 

 

 

 

f(x)

 

 

 

 

 

dy

 

 

к экстремали.

 

A

 

 

dx

 

 

 

 

 

 

 

 

Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

разность

 

между

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

функцией сравнения

 

 

 

 

 

 

 

 

и

экстремалью

при

 

 

 

 

 

 

x

 

одном

и

том

же

 

 

 

 

 

 

 

значении

аргумента

 

a

 

 

 

b

 

 

 

 

 

 

 

 

x

 

 

называется

 

 

 

 

 

 

 

 

 

 

Рис.18.

 

 

 

 

 

вариацией

 

y

 

 

 

 

 

 

 

 

 

 

функции f (x) :

 

 

 

 

 

 

 

 

 

 

 

y f (x) f (x)

 

 

 

 

 

 

 

 

 

Вариацию y не следует путать с дифференциалом dy .

Дифференциал dy -

бесконечно малое приращение функции,

обусловленное бесконечно малым приращением аргумента

dx . Вариация же – бесконечно малое приращение функции

f (x) , которое создает новую функцию f x y y .

 

 

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

 

Выделим

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

семейство

 

функций

сравнения, бесконечно близких к экстремали.

 

 

 

 

 

y x, f x y

 

 

 

 

 

 

 

 

В

это

семейство

входи

и

сама

экстремаль:

при

0 y x,0 f x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

90