Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
диф 1.pdf
Скачиваний:
121
Добавлен:
05.06.2015
Размер:
2.19 Mб
Скачать

94

P = nπ 2.

EI l

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

P = n2π 2EI .

(2.97)

l2

 

Число n =1,2,...в данной задаче имеет смысл числа полуволн си-

нусоиды (2.94), по которой происходит искривление оси стержня.

Из формулы (2.97) следует, что существует множество значений сжимающей силы P, при которых возможны различные искривлëнные формы равновесия продольно сжатого стержня. Наименьшее значение сжимающей силы, соответствующее n = 1,

называется критической (эйлеровой) силой и равно

P = π 2EI . l2

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

y(0)y(1) = 0,

y (1) = 0.

Общее решение уравнения согласно (2.92) имеет вид y = C1 cos λx +C2 sin λx.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вычисляем

 

 

λ sin

 

λx +C2 λ cos

λx,

y (x) = −C1

 

 

 

 

y(0) = C1,

y(1) = C1 cos

 

 

+C2 sin

 

 

,

 

λ

λ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ sin

λ +C2

 

 

λ cos λ.

 

y (1) = −C1

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

λ) C1 sin λ C2 = 0,

(1cos

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

λ C1 + λ cos λ C2 = 0.

λ sin

 

95

Ненулевое решение еë существует тогда и только тогда, к о- гда определитель системы равен нулю:

1cos λ

sin λ

= 0.

λ sin λ

λ cos λ

Раскрывая определитель, получим уравнение относительно параметра λ :

λ (cosλ 1) = 0,

корни которого являются собственными значениями задачи:

λ = 0,

λ = (2πn)2,

n =1,2,...

Легко убедиться,

что при λ = 0

y = const. Подставляя

λ = 2πn, n =1,2,..., в систему линейных однородных алгебраических уравнений относительно C1 и C2 , находим C2 = 0. Соот-

ветствующие собственные функции с точностью до множителя будут y = cos 2πnx.

2.9. Дифференциальные уравнения с переменными коэффициентами

Если коэффициенты ai (x), (i =1,2,...,n) линейного неоднородного уравнения

y(n) + a1 (x) y(n1) + a2(x) y(n2) + an1 (x) y′+ an (x) y = f (x) (2.98)

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

Одним из простейших линейных дифференциальных уравнений с переменными коэффициентами является уравнение Эйлера.

96

2.9.1. Уравнение Эйлера

Уравнением Эйлера называется линейное уравнение вида

(ax +b)n y(n) + a (ax +b)n1 y(n1) + a

(ax +b)n2

y(n2) + +

 

1

2

 

 

(2.99)

+ an1 (ax

+b)y′+ an y = f (x),

 

 

 

где a, b, a1, a2, , an1, an - константы. В наиболее

распространëнном случае,

при a =1, b = 0

уравнение Эйлера

имеет вид

 

 

 

 

xn y(n) + a1xn1 y(n1) + a2xn2 y(n2) + + an1xy′+ an y = f (x). (2.100)

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

Для уравнения (2.100), полагая

 

 

 

 

 

 

 

 

 

 

 

 

 

x = et , (x > 0) t = ln x,

 

 

 

 

(2.101)

находим производные разного порядка от функции

y по новой

переменной

t:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y′ =

dy

=

dy dt

= e

t

dy

,

 

 

y′′ = e

2t d 2 y

dy

,

 

 

 

 

 

 

 

 

 

2

 

 

dx

 

dt dx

 

 

 

dt

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

y′′′ = e

3t d 3 y

3

d 2 y

+ 2

dy

 

 

 

 

(2.102)

 

 

 

 

 

 

3

 

2

 

,

 

 

 

 

 

 

 

 

 

 

 

dt

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

…………………………………

и подставляя (2.101) и (2.102) в (2.100), получим уравнение с постоянными коэффициентами относительно функции y(t).

В общем случае, для уравнения (2.99) формулы перехода (2.102) примут вид:

 

 

 

 

 

t

 

 

 

 

 

 

y

== ae

t

dy

 

ax +b = e ,

t = ln(ax +b),

 

 

dt ,

 

 

 

 

 

2

e

2t d 2 y

dy

,

3

e

3t d 3 y

3

d 2 y

+ 2

dy

y′′ = a

 

 

2

 

y′′′ = a

 

 

 

 

3

 

 

2

,

 

 

 

dt

 

 

 

 

 

 

 

dt

 

dt

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

dt

…………………………………..

97

 

Пример. Решить уравнение

 

x2 y′′+3xy′+ y = 2+ x.

(2.103)

Уравнение (2.103) есть уравнение Эйлера второго порядка. Применим замену независимой переменной и производных по формулам (2.101), (2.102). Тогда уравнение примет вид

e

2t

e

2t d 2 y

dy

+3e

t

e

t dy

+ y = 2

+e

t

 

 

 

2

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

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

d 2 y

+ 2dy

+ y = 2+ et .

(2.104)

dt 2

dt

 

 

Его общее решение

y = y0 + y1 + y2,

где y0 - общее решение соответствующего однородного уравнения:

y0 = C1et +C2 t et ,

y1 - частное решение неоднородного уравнения

d 2 y

+ 2

dy

+ y = 2

,

dt

2

dt

 

 

 

y2 - частное решение неоднородного уравнения

d 2 y

+ 2dy

+ y = et .

dt2

dt

 

Эти решения находятся методом подбора (см. выше) и имеют вид

y1 = 2, y2 = 0,25et .

Поэтому общее решение уравнения (2.104) будет

y = C et +C

2

t et + 2+0,25et .

(2.105)

1

 

 

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

 

 

98

 

 

y = C 1 +C

 

ln x

+ 2+0,25x.

 

 

 

 

1 x

2 x

Замечание. Решение однородного уравнения Эйлера мож-

но искать в виде

 

 

 

 

y = ekt = (et )k = xk .

Тогда y′ = kxk1, y′′ = k(k 1)xk2,...Подстановка этих соотно-

шений в

исходное дифференциальное уравнение (2.100) (при

f (x) = 0)

приводит его сразу к соответствующему характеристи-

ческому уравнению относительно k.

Простому корню k1 характеристического уравнения соответствует частное решение xk1 , а m-кратному действительному

корню k - m линейно независимых решений вида:

xk1 , xk1 ln x,

1

 

xk1 (lnx)2,..., xk1 (ln x)m1.

 

Если коэффициенты дифференциального уравнения действительны, а характеристическое уравнение имеет комплексные сопряженные корни кратности m, то уравнение Эйлера будет иметь 2m частных линейно независимых решений вида:

xα cos(β ln x),

xα sin (β ln x),

xα ln xcos(β ln x),

xα ln x sin (β ln x),

……………………..

……………………..

xα (lnx)m1 cos(β ln x),

xα (lnx)m1 sin (β ln x).

В общем случае решения уравнения с переменными коэффициентами (2.98) строятся приближёнными методами. Одним из эффективных приближённых методов является метод степенных рядов [5,9,10].

Рассмотрим два варианта применения этого метода:

1)метод степенных рядов в форме метода неопределённых коэффициентов;

2)разложение решения задачи Коши в степенной ряд с помощью формулы Тейлора16.

16Т е й л о р Брук (18.08.1685 – 29.12.1731) – английский математик.