Материал: Методы решения задач проектирования технологических процессов обработки давлением. Шагунов А.В., Корольков В.И

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

Будем рассматривать значение функционала только на

семействе кривых y x, . В этом случае функционал превращается в функцию параметра 14.

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

При 0 функция достигает экстремума,

поскольку y x, становится экстремалью f (x) . Необходимым условием экстремума функции является

равенство нулю ее первой производной:

d J f (x) y 0 d

Производная ' в точке 0 называется первой вариацией функционала J

J d

d 0

Таким образом, необходимым условием экстремали является равенство нулю первой вариации функционала15

y(x) f (x) J 0

Эйлер показал, что первая вариация функционала

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

Приводится к виду

 

 

F

 

d F

 

b

 

 

 

 

 

 

 

 

 

J

 

 

 

 

ydx

 

a

y

 

dx y'

Тогда необходимое условие стационарности функционала J 0 эквивалентно уравнению:

14 Действительно, если считать уравнение функционали y f x

известным, то функция F зависит от x,

15 Можно увидеть определенную аналогию между дифференциалом функции и первой вариацией функционала.

91

F d F 0

y dx y'

Это уравнение впервые получено Эйлером и носит его

имя.

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

– это прямая, поэтому можно легко проверить правильность полученного результата.

Как было показано ранее, эта задача сводится к отысканию экстремали функционала:

 

J b

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 y' 2

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

F

 

 

y'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

y'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 y' 2

 

 

 

 

 

 

 

 

Здесь F(x, y, y')

 

, y

,

 

 

1 y' 2

.

 

 

 

 

 

 

 

 

 

d

 

 

 

y'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Уравнение Эйлера:

dx

 

1 y' 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y'

 

 

 

 

C const

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Откуда

 

 

 

1 y' 2

 

 

 

 

 

 

 

 

.

 

Решение этого

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

обыкновенного дифференциального

уравнения

получим

следующим образом:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

C

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 y' 2

 

 

 

 

 

 

 

 

 

 

y C1x C2

 

 

 

 

 

 

 

1 C 2

 

 

 

 

 

 

 

 

,

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

y A x a ,

y B x b . Тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

A C1a C2

C1

A B

,C2

A a

 

A B

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a b

 

 

 

b

B C1b C2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

Окончательно получим уравнение экстремали

92

y A B x A a b

A B a A a b

A B x a a b

Очевидно, что это уравнение прямой, проходящей через

точки a, A и b, B , что и требовалось доказать.

В общем случае функционал может быть определен для нескольких функций:

J y1, y2 , , yn F x, y1, y2 , , yn , y1 ', y2 ', , yn ' dx

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

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

F

d F

0

 

 

 

 

yi

dx yi '

 

Поскольку рассуждения приемлемы для любой функции, то необходимым условием экстремума функционала будет система дифференциальных уравнений

Fy'i dxd Fy'i' 0, i 1,2, , n

определяющая систему экстремалей. Метод Ритца

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

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

93

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

Пусть неизвестная функция y x , являющаяся экстремалью функционала J имеет вид:

y x ai i x

i 1

Здесь i (x) - некоторые заранее заданные, т.е.

известные функции (Ритц назвал их координатными), ai - неизвестные коэффициенты, подлежащие определению.

Координатные функции i (x) должны удовлетворять следующим условиям:

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

сравнения. (В нашем примере с кривой, проходящей через две точки, это означает, что все они должны

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

граничным

условиям

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

 

 

2.Удовлетворять условиям полноты – ряд, построенный

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

Тогда функционал J y принимает вид:

J y x J

 

a

 

 

 

x

 

i

i

 

 

i 1

 

 

.

Поскольку координатные функции i (x) заранее заданы, то функционал становится функцией коэффициентов

ai 17:

16Искомыми функциями в задачах теории пластичности является функции изменения скоростей материальных точек в объеме деформируемого тела.

17Зависит от переменных, а не от функции.

94

J J ai .

Условие стационарности функционала

J 0

в этом случае преобразуется в систему бесконечного числа уравнений:

J

0

i 1,2 ,

ai

 

 

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

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

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

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

подходящих функций.

Последовательность решения задачи определения напряженно-деформированного состояния методом Ритца.

1. Выбор подходящих функций и количества членов ряда.

95