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

b

.pdf
Скачиваний:
22
Добавлен:
17.03.2016
Размер:
1.59 Mб
Скачать

если для любой непрерывной функции η(х) С0 , η(а) = η(b) = 0,

х [a,b]

b F(x) η(x)dx = 0 ,

a

где функция F(x) непрерывна на отрезке интегрирования [a,b] (F(x) С0), то F(x) 0 на том же отрезке.

 

у

 

 

 

 

 

 

 

Доказательство:

 

 

 

 

 

 

 

 

Предположим,

что

в

 

 

y(х0)0

 

 

 

 

 

точке

х0

 

[a,b]

 

 

 

y(х)

 

 

F(x0)

0.

Тогда

из

 

 

η(х)

 

 

 

 

 

 

 

 

 

непрерывности

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

функции

 

F(x)

0

 

 

 

 

 

 

 

х

следует, что если F(х0)

 

 

 

 

 

 

 

0, то

F(x) сохраняет

 

 

 

 

 

 

 

 

а

x1

x0

 

x2

b

 

 

знак

в

некоторой

 

 

 

 

 

 

 

 

 

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

x1

Рис. 1.4. К доказательству основной леммы

< x0 < x2. Но тогда,

выбрав функцию η(х),

 

 

вариационного

исчисления

также

сохраняющей

 

 

 

 

 

 

 

 

 

< x0 < x2

 

 

 

 

 

 

 

знак в окрестности

x1

и равной нулю вне этой окрестности (рис. 1.4), получим:

 

 

b F(x) η(x)dx = x2F(x) η(x)dx 0 ,

 

 

 

 

 

 

a

 

 

 

x1

 

 

 

 

 

так как

произведение

F( x ) η( x ) сохраняет знак на отрезке

x1 < x0 < x2 и равно нулю вне этого отрезка. Так как, предположив, что в некоторой точке x0 отрезка [а, b] F(х0) 0, мы пришли к противоречию с условиями леммы, то, следовательно, F(x) 0 на отрезке интегрирования [а, b] .

Замечание. Доказательство леммы не изменится, если функции F(x) и η(х) будут принадлежать к классу Сn, то есть быть непрерывными и иметь непрерывные производные до п-го порядка.

Основная лемма может быть доказана также для произвольной области D на плоскости, т.е., если

∫∫F(x, y) η(x, y)dA = 0 ,

F(x, y) η(x, y) C

D

 

при произвольной функции η(х,у) С (η(х,у) = 0 на границе области D), то F(x,y) 0.

Далее лемма может быть обобщена на трехмерное и любое n- мерное пространство.

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

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

вида

 

b

 

 

J[y( x )]= F[x, y(x), y (x)]dx ,

 

 

 

a

 

у(х) С1 при у(а) = А, у(b) = В. Функцию

F[x,y(x),y'(x)]

будем

считать дважды дифференцируемой.

 

 

 

Рассмотрим ε-окрестность функции у(х), положив

 

y1 (x)= y(x,α)= y(x)+α η(x);

y1( x ) = y (x,α)= y ( x ) +α

η (x);

α η( x ) =δ y( x ) и ϕ(x,α)= J[y(x,α)]= b F[x, y(x,α), y(x,α)]dx ;

η(a)=η(b)= 0 .

a

Получим вариацию функционала, как производную функционала по параметру α, положив α = 0 (см. формулу (1.1.6)). Так как интеграл имеет конечные пределы и функция

непрерывна и имеет непрерывные

F(x, y(x,α), y (x,α))

производные, можно проводить дифференцирование под знаком интеграла. При этом учитываем, что функция F[x,y(x,α),y'(x,α)] является сложной функцией аргументов х, α, y(x,α), y'(x,α).

ϕ( x,α )

 

 

= F[x, y(x,α), y (x,α)]dx =

 

 

 

 

b

 

 

 

α

 

a

 

α

 

 

 

 

 

 

b

 

y

 

 

y

b

= Fy

 

 

+

Fy

 

dx = [Fy η + Fyη]dx ,

 

α

a

 

 

 

α

a

где Fy =

F[x, y(x), y(x)]

 

,

Fy

=

F[x, y(x), y(x)]

.

 

 

 

 

 

 

y

 

 

 

 

 

y

 

 

 

 

Проинтегрируем по частям второе слагаемое в полученном

выражении

 

 

 

 

dF

 

 

 

 

dF

 

 

 

b

 

b

b

y

 

b

y

 

 

Fyη(x)dx = Fyη(x)

a

 

η(x)dx = −

 

η(x)dx .

 

 

 

 

a

 

 

 

a

dx

 

 

a

dx

Подставляя

полученное выражение в производную по α

функции ϕ(х,α)

и вычисляя ее при

 

α = 0, получим вариацию

функционала

 

 

 

b

 

 

 

 

 

 

 

 

 

 

 

J[y(x,α)]

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

F[x, y(x), y (x)]dx ;

 

 

 

a

δ J[y(x)]=

ϕ(x,α)

 

b

 

 

= Fy

 

α

 

α=0 a

 

 

dFdxy

 

b

dF

y

 

 

 

η(x)dx = Fy

 

 

δy dx .

dx

 

a

 

 

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

δ J[y(x)]= b

Fy

dFy

η(x)dx .

(1.2.4)

dx

a

 

 

 

Так как η(x) является произвольной функцией, то на основании основной леммы вариационного исчисления получаем необходимое условие экстремума функционала в виде:

Fy

dFy

= 0 .

(1.2.5)

dx

 

 

 

 

Формула (1.2.5) была получена Эйлером и носит его имя.

В формуле (1.2.5) выражение

 

dFy

представляет полную

 

dx

 

 

 

 

производную по х и вычисляется по формуле:

dFy

=

d F[x, y(x), y(x)]

=

dx

dx

y

 

 

=

2 F

 

+

2 F

y′ +

2 F

 

y′′ = F

xy

+ F

yy

y′ + F ′ ′ y′′. 1.2.6)

 

 

 

 

xy

 

 

yy

 

 

 

y y

 

 

 

 

 

 

Формула Эйлера (1.2.5) является необходимым условием

существования

 

 

 

экстремума

 

функционала

J[y(x,α)]= b F[x, y(x), y(x)]dx .

a

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

δ2J[y(x)] = δ{δ J[y(x)]} > 0,

имаксимум функционала, если вторая вариация меньше нуля

δ2J[y(x)] = δ{δ J[y(x)]} < 0.

Если вторая вариация функционала равна нулю

δ 2J[y(x)] = δ{δ J[y(x)]} = 0,

 

 

то говорят, что на функции

у(х)

функционал

достигает

стационарного значения.

 

 

 

 

 

 

J[y0(x)] в

Для доказательства

представим

функционал

ε-окрестности функции

y(x),

для

которой

первая

вариация

функционала равна нулю

-

δ J[y0(x)]

=

0, в виде первых трех

членов разложения функционала в ряд Тейлора:

 

 

J[y(x)]= J[y0 (x)]+ δJ[y0 (x)]+δ 2 J[y0 (x)].

 

Так как δ J[y0(x)] = 0, то будем иметь

 

 

 

δ 2 J[y(x)]= J[y(x)]J[y0 (x)],

 

 

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

 

 

 

 

 

 

 

1) если δ 2J[y0(x)] > 0 , то

J[y(x)] > J[y0(x)]

в ε-окрестности

функции y0(x) - функционал достигает минимума

на кривой

y0(x);

 

 

 

 

 

 

 

2) если δ 2J[y0(x)] < 0 , то J[y(x)] < J[y0(x)]

в ε-окрестности

функции y0(x) - функционал

достигает

максимума

на кривой

y0(x);

 

 

 

 

 

 

 

3) если δ 2J[y0(x)] = 0 , то

J[y(x)] = J[y0(x)]

в ε-окрестности

функции y0(x) - признак стационарности функционала.

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

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

Без доказательства приведем некоторые обобщения формулы

Эйлера.

 

 

 

 

 

 

 

 

1.

Для функционала

 

 

 

 

 

 

 

 

 

b

 

 

 

 

(n)

( x ))(x)x ,

 

 

 

′′

 

 

 

 

 

 

 

J[y(x)]= F(x, y(x), y (x), y

 

(x),..., y

 

 

 

a

 

 

 

 

 

 

 

при

у(а) = уа,

у(b) = уb,

у'(а)

 

= у'a,

у'(b)

= у'b, …,

y( n1 )( a ) = ya( n1 ) ,

y( n1 )( b ) = yb( n1 ) ,

 

 

экстремум

функционала

реализуется для функции у(х), отвечающей условию:

 

F

d

 

F

+

d 2

 

F

d 3

 

F

+ ... + (1)n

 

d n

 

 

F

 

= 0 .

(1.2.7)

y

 

 

dx2

 

dx3

 

dxn y(n)

 

 

dx y

 

y′′

y′′′

 

 

 

 

 

 

2. Для функционала от функции нескольких переменных,

например, в

трехмерном пространстве

х, у,

z

 

для

взаимно

независимых функций f(x, у, z), ϕ (x, у, z), ψ(x, у, z) С

 

 

J = ∫∫∫F( x, y,z, f ,ϕ,ψ , f x ,ϕx ,ψ x , f y ,ϕy ,ψ y , f z ,ϕz ,ψ z )d,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

f x = f ,

f y

 

= f ,

f z =

f , ϕx =

ϕ

,.....,

ψ z

=

ψ

.

 

 

z

 

 

 

 

x

 

 

y

 

 

 

z

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

На границе области

 

 

 

 

 

 

 

 

 

 

 

 

 

 

интегрирования функции

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f, ϕ, ψ удовлетворяют

 

 

 

 

 

 

 

 

 

 

 

 

 

 

граничным условиям

f = fs, ϕ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=ϕs,

ψ = ψs.

 

 

Так как функции f(x, у, z),

ϕ(x, у, z),

ψ(x, у, z) - независимы, то

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

δ J f

= 0 ,

δ Jϕ = 0

 

 

δ Jψ

= 0 ,

 

 

 

(1.2.8)

которым соответствуют формулы Эйлера:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

δ J f = 0

F

 

F

 

 

F

 

 

F

;

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x f x

 

y f y

 

z f z

 

δ Jϕ = 0

F

∂ ∂F

 

∂ ∂F

 

 

 

 

F

;

(1.2.9)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ

x ϕx

 

y ϕy

 

 

 

 

 

 

 

 

 

 

 

z ϕz

 

δJψ = 0

F

F

F

 

F

.

 

ψ

x ψ x

 

 

 

 

 

 

 

 

 

y ψ y

 

 

 

 

z ψ z

 

Частные производные по х, у, z являются полными частными производными и вычисляются по формулам:

∂ ∂F

=

 

2 F

+

 

2 F

 

f

+

2 F

2 f

+

 

2 F 2 f

+

 

 

2 F

2 f

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x f x

xf x

 

f f x x

f x2

x2

f y f x xy

f z f x xz

 

 

 

 

 

 

 

 

 

 

 

 

∂ ∂F

 

=

 

2 F

 

+

 

2 F

 

f

+

 

2 F 2 f

 

+

2 F

2 f

+

 

 

2 F

 

2 f

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y f y

 

yf y

 

 

f f y y

 

f x f y xy

f y2

y 2

 

 

f z f y yz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ ∂F

=

 

2 F

+

 

2 F

 

f

+

 

2 F 2 f

 

+

 

2 F 2 f

 

+

2 F

 

2 f

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z f z

 

zf z

f f x z

f x f z xz

 

f y f z

 

yz

 

f z2

 

z 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично могут быть получены формулы для полных частых производных подынтегральной функции F при вариациях функционала по аргументам ϕ(x, у, z), ψ(x, у, z).

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

1.3. Задачи на экстремум функционалов.

Рассмотрим несколько задач на минимум (экстремум) функционалов и использование формулы Эйлера.

Пример 1. Определить функцию у(х), на которой достигается экстремум функционала

π

2

 

2

 

 

π

 

 

 

J1[y( x )]= (y

y

)dx , у(0) = 0,

(1.3.1)

 

 

y

=1.

0

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

Для решения задачи

 

 

 

 

используем условия Эйлера

 

 

 

 

(1.2.5) минимума функционала

 

 

 

 

δJ1

[y(x)]= 0 Fy

d

Fy= 0

 

 

 

 

 

 

 

 

 

 

 

 

.

dx

 

 

 

 

 

 

 

 

 

Для данного функционала имеем

F(x, y, y)= y 2 y2 .

 

 

Тогда получаем

Fy = y (y 2 y2 )= 2y ; Fy= y(y 2 y2 )= −2y;

 

d

Fy= −2

dy

= −2y′′ .

 

 

dx

dx

 

 

 

 

 

Подставляя результаты в формулу Эйлера, получим

 

 

 

2у + 2у" = 0 или у" + у = 0.

(1.3.2)

Таким

образом, условие Эйлера минимума функционала

J1[y(x)]

привело

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

уравнению второго порядка с постоянными коэффициентами

(1.3.2)

и граничными

условиями, известными из постановки

задачи

у(0) = 0,

 

π

=1.

y

 

 

 

 

2

 

Составляя характеристическое

уравнение

к2 + 1 = 0,

имеющее два чисто мнимых корня

k1,2 = ±i,

получим общее

решение дифференциального уравнения (1.2.2) в виде:

у(х) = С1 sin x + C2 cos x.

Из граничных условий определяем

 

 

 

у(0) = С2 = 0,

 

π

= C1

=1 .

 

y

 

 

 

 

 

2

 

 

 

Окончательно получаем, что функционал J1[y(x)] достигает

экстремума на кривой

у(х) = sin x.

 

 

 

 

 

Пример 2. Определить функцию на которой достигается

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

 

 

 

 

 

 

J1[y( x )]= 2 (y2 + y2 )dx ,

 

у(0) = 0,

у(2) = 1. (1.3.3)

0

 

F(x,y,y')

 

 

Подынтегральная

функция

отличается

от

предыдущего примера лишь знаком перед

у',

поэтому из условия

Эйлера получим дифференциальное уравнение

y y′ = 0 , общим

решением которого является функция

 

 

 

 

у(х) = С1 ех+ С2 е .

 

 

 

 

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

 

 

 

у(0) = C1 + С2 = 0,

С2 = -C1;

 

 

у(2) = С1 е2+ С2 е -2 = С1 (е2 - е -2 ) = 2 С1 sh2 = 1,

C =

1

.

 

 

 

1

2sh2

 

 

 

Откуда получаем функцию y(x), на которой

функционал

J2[y(x)] достигает экстремума

 

 

 

 

y( x ) = 2sh1 2 (ex ex )= shshx2 .

Пример 3. Определить функцию, на которой достигается экстремум функционала

J3 [y( x )]= 2 (y2

+ y sin x)dx , у(0) = 0,

у(π) = 1.

 

 

(1.3.4)

Решение:

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dFy

 

 

 

 

 

 

 

F = у

2

+ у sinx,

Fy = sin х,

Fy= 2y',

 

 

 

= 2y

′′

;

 

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

Fy

dFy

= sin x 2y′′ = 0 ,

y′′

=

1

 

sin x .

 

 

 

 

 

 

 

dx

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрируя

дважды

найденное

из уравнения

 

 

 

Эйлера

выражение и используя граничные условия, получим:

 

 

 

 

 

 

 

 

 

 

y(x)= C x + C

2

1

sin x ,

у(0) = С2 = 0,

у(π) = С1 π = 1,

C

 

=

 

1

.

2

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

π

 

Окончательно получаем, что экстремум рассматриваемого функционала реализуется на кривой

y(x)= πx 12 sin x .

Пример 4. Найти линию минимальной длины, соединяющей две точки на плоскости с координатами х1, у1 и х2, у2.

Хотя ответ на этот вопрос ясен, получим его, используя формулу Эйлера.

Расстояние между двумя точками плоскости определяется функционалом:

L[y( (x)]= x2

1 + y, у(х1) = у1, у(х2) = у2.

(1.3.5)

x1

 

 

Так как функция F не зависит от у(х) , то Fy = 0, и формула Эйлера запишется в виде:

dFdxy= Fyyy′′ = 0 ,

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

у" = 0

или

Fyy = 0.

Решением первого дифференциального уравнения будет функция у(х) = С1х + С2, которая определяет прямую линию на плоскости. Из граничных условий получим:

y(x1) = C1x1 + С2 = y1 и y(x2) = C2x2 + С2 = y2.

Решая полученную систему уравнений, имеем

C

=

y2 y1

,

C

2

= y

y2 y1

x ,

 

 

1

 

x2 x1

 

1

x2

x1

1

 

 

 

 

 

 

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

y( x ) = y1 + y2 y1 (x x1 ) x2 x1

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

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

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

Пример 5. Рассмотрим задачу о минимальной поверхности вращения.

Найти кривую у(х), проходящую через точки у(х0)= у0, у(х1)= у1, которая при вращении вокруг оси х образует поверхность минимальной площади.

Площадь поверхности вращения определяется интегралом

x

 

 

S = 1

y( x ) 1 + y2 (x)dx .

(1.3.6)

x0

 

 

 

Находим минимум

 

функционала (1.3.6), используя

 

условие Эйлера.

Предварительно заметим, что, так как F = y( x )

1 + y2 ( x )

не зависит напрямую от х (F = F[y(x),y'(x)]), то в формуле (1.2.6) частная производная по х равна нулю, и, следовательно, условие Эйлера принимает вид:

F

y

d

F

= F

y

F

yy

y′− F ′ ′ y′′ = 0 .

(1.3.7)

dx

 

 

y

 

 

y y

 

Если выражение (1.3.7) домножить на у', то оно становится полной производной по х

(Fy FyyFyyy′′) y′ = dxd (F Fyy)= 0

. 1.3.8)

В справедливости полученного выражения можно убедится непосредственным дифференцированием:

d

(F Fyy)

=

(F Fyy)

y′ +

(F Fyy)

y′′ =

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

y

 

 

 

 

 

 

= Fy y′ − Fyyy2 + Fyy′′ − Fyyyy′′ − Fyy′′ =

 

 

= (F

y

F

yy

y′ − F ′ ′ y′′) y.

 

 

 

 

 

 

 

 

 

 

 

 

 

y y

 

 

 

 

 

 

 

 

 

 

Интегрируя уравнение (1.3.8), получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F Fyy′ = C 1.

 

 

 

 

(1.3.9)

 

Так как

 

F

=

 

 

1 + y

2

 

=

yy

, то из (1.3.9) имеем:

 

 

y

y

 

 

1 + y2

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y 1+ y2

 

yy

 

y′ = y (1+ y2 )yy2

=

y

 

= C1 .

 

 

 

 

 

 

1+ y2

 

 

 

 

1+ y2

 

 

1+ y2

 

 

Откуда получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y′ =

y

2 C12

 

 

 

 

 

C dy

 

 

 

 

 

 

 

 

 

C1

 

,

 

или

 

1

 

= dx .

 

 

 

 

 

 

 

 

 

 

y2 C12

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Интегрирование последнего выражения дает

 

 

 

y + y2

C 2

x + C2

= C1

ln

 

1

C1

 

 

 

 

 

или

 

 

 

x+C2

 

 

 

 

 

 

 

 

 

 

y + y2 C 2

= C e C1 .

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x+C2

 

 

Умножая последнее соотношение на

 

y

2

C

2

C1

,

 

y

 

 

 

e

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= C e

x+C2

 

 

 

 

 

 

 

 

 

 

y y2 C 2

 

C1 .

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Складывая полученное соотношение с предыдущим,

получим решение в явном виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(x)= C ch

x + C2

.

 

 

 

 

(1.3.10)

 

 

 

 

 

 

 

1

 

C1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

Пример 6. Получить кривую, минимизирующую функционал:

J[y( x )]= l

 

EI z ( x )

y′′2

( x ) q( x ) y( x ) dx ,

(1.3.11)

 

0

2

 

 

 

у(0) = 0,

у′′(0) = 0, у( l ) = 0,

у′′( l ) = 0,

 

где у(х) - прогиб балки; ЕIz - изгибная жесткость балки (Iz - момент инерции поперечного сечения балки, может быть переменным); q(x) - интенсивность распределенной нагрузки.

Функционал (1.3.11) представляет собой полную энергию деформации изгиба балки, а минимизация этого функционала при заданных граничных условиях, как следует из принципа Лагранжа, который будет рассмотрен ниже, дает форму прогиба балки от действия распределенной нагрузки q(x). Приведенные граничные условия соответствуют шарнирному опиранию балки.

 

EI

z

( x )

y′′2

 

зависит от у"(х),

Функция

F =

 

 

( x ) q( x ) y( x )

 

2

 

 

 

 

 

и условие Эйлера определяется формулой (1.2.7):

 

 

 

 

 

F

d

 

 

F

 

+

d

2 F

= 0 ;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

d x y

d x2 y′′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

= −q( x ) ;

 

 

 

 

F

 

= 0 ;

 

 

 

 

 

 

 

F

= EI z ( x ) y′′ ;

y

 

 

 

 

y

 

 

 

 

 

 

y′′

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2 F

 

 

d 2

 

 

 

 

 

d 2 y

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

EI

z

( x )

 

 

 

 

 

 

= (EI

z

( x ) y′′)

;

 

 

 

2

 

 

 

 

 

2

 

 

 

 

2

 

 

dx

 

y′′

 

 

dx

 

 

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и, следовательно, условие минимума функционала получим в виде:

d 2

 

 

 

d 2 y

 

 

 

 

EI

z

( x )

 

 

 

q( x ) = 0 .

(1.3.12)

 

2

 

2

dx

 

 

dx

 

 

 

 

 

 

 

 

 

 

 

Таким образом, условию Эйлера минимума функционала (1.3.12) соответствует дифференциальное уравнение изгиба балки.

В случае балки постоянной жесткости это уравнение запишется в виде:

y IV q( x ) = 0 , EI z

интегрируя которое при q = const , получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

2

3

 

 

qx4

 

 

 

 

 

 

у(х) =

С0 +

С1 х + С2 х

 

+ С2 х

 

 

.

 

 

 

 

24EI z

Удовлетворяем граничные условия задачи при

х = 0:

 

 

 

 

 

 

 

у(0) = С0 = 0;

 

 

 

у′′(0) = 2С2 = 0.

 

 

 

 

 

Учитывая вычисленные коэффициенты, удовлетворяем

граничные условия при

х = l:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y(l)= C1 l + C3

3

 

ql 4

 

 

 

 

′′

 

 

 

ql 2

 

 

l

 

 

 

 

y

l

 

 

 

 

 

 

 

;

 

 

 

= 0 ;

24EI z

 

(l)= 6C3

 

2EI z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

 

 

ql

 

 

 

C

 

 

ql 3

 

 

 

 

ql 3

 

 

ql 3

 

 

3

=

 

 

 

,

 

= −

 

 

 

 

+

 

 

= −

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

12EI z

 

1

 

12EI z

 

 

24EI z

 

24EI z

 

 

 

 

 

 

 

 

 

 

 

 

 

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

y( x ) = − 24qEI z (l 3 x 2lx3 + x45 ),

что соответствует точному решению сопротивления материалов. В частности, при х = 0,5l получим прогиб в середине балки

l

 

q

 

3 l

 

l 3

 

l 4

 

 

5

 

ql 4

y

 

 

= −

 

l

 

 

2l

 

 

 

+

 

 

 

 

= −

 

 

 

.

2

24EI z

2

2

2

384

EI z

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

1.4.Изопериметрическая задача вариационного исчисления.

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

Среди кривых, удовлетворяющих условиям у(а) = А, у(b) = В, на которых функционал

J1[y( x )]= b G(x, y, y)dx

(1.4.1)

a

 

принимает заданное значение l, найти ту, для которой другой функционал

J[y( x )]= b F(x, y, y)dx

(1.4.2)

a

 

достигает экстремума.

Решение этой задачи определяется теоремой [16,17], которая здесь приводится без доказательства.

Теорема. Если функционал J[y( x )]= b F(x, y, y)dx достигает

a

экстремума на кривой у = у(х), у(а) = А, у(b) = В, а функционал J1 [y(x)] на этой кривой удовлетворяет условию

J1 [y(x)]= b G(x, y, y)dx = l , и кривая у = у(х) не является

a

экстремалью функционала J1[y(x)], то существует такая постоянная λ, что эта кривая является экстремалью функционала

J0 [y( x )]= b (F + λ G)dx .

(1.4.3)

a

 

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

Решение, удовлетворяющее условиям Эйлера для функционала (1.4.3), даст семейство кривых, определяемых параметром λ, из которых искомое решение должно удовлетворять условию наложенному на функционал (1.4.1).

Пример. Рассмотрим задачу о наибольшей площади, описанную выше. Найти кривую в верхней полуплоскости,

проходящую через точки

 

(-a,0) и

(a,0),

имеющую заданную

длину

2l (l

> а) и охватывающую вместе с отрезком

[-а, а]

максимальную площадь.

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Ищем кривую у = у(х), для которой

 

 

 

 

 

 

 

 

J1 [y(x)]= a

1 + y2 dx = 2l , у(-a) = у(a) = 0,

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

а интеграл

(функционал)

J[y(x)]= a y(x)dx

 

 

принимает

максимальное значение.

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Составим функционал

 

 

 

 

 

 

 

 

 

 

 

 

J

 

[y(x)]= J[y(x)]+ λ J

 

[y(x)]=

a

+ λ

1 +

y

2

 

 

(1.4.4)

0

1

y

 

 

dx .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

Тогда F = y + λ 1 + y2 ;

Fy =1;

Fy

= λ

 

 

y.

 

 

 

 

 

 

 

 

 

 

 

 

1 + y2

 

Составим условие Эйлера

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 λ

d

y

= 0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

1 + y2

 

 

 

 

 

Из этого условия, интегрируя его, находим

 

 

 

 

 

 

 

 

x λ

y

= C1 или

 

λ

y

= x1

C1 .

 

 

 

 

 

 

 

 

 

 

 

 

1+ y2

 

 

 

 

 

1 + y2

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]