Будем рассматривать значение функционала только на
семействе кривых 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