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

40

Таким образом, с математической точки зрения задача об изгибе консольно защемлëнного стержня сводится к задаче Коши (2.2), (2.3), а задача об изгибе шарнирноëопртого стержня – к краевой задаче (2.2), (2.4).

2.1.Интегрирование дифференциальных уравнений n – го порядка методом понижения порядка

Если правая часть уравнения (2.1) является известной непрерывной функцией f (x) только одной переменной x или не со-

держит искомую функцию y: f (x,y,y′′, ,y(n1) ), или не содер-

жит явно независимую переменную x: f (y,y,y′′, ,y(n1) ), то для

решения уравнения (2.1) может быть применëн метод понижения порядка.

1. y(n) (x) = f (x). (2.5)

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

y(n1) = f (x)dx +C1,

y(n2) = dxf (x)dx +C1x +C2,

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

В результате общее решение уравнения будет иметь вид

y = ∫∫ ∫ f (x)dx + C1 xn1 + C2xn2 + + Cn1 x + Cn ,

(2.6)

где в правой части - n-кратный интеграл от функции f (x)

и мно-

гочлен (n-1)-ой степени от x, коэффициентами которого являются

n произвольных постоянных.

 

 

 

 

 

 

 

 

Если

для уравнения (2.5) решается задача Коши с началь-

ными условиями

y(x0 ) = y0,

y(x0 ) = y0,…, y(n1) (x0) = y0

(n1), то

частное решение уравнения может быть записано в виде

 

 

x

 

x

 

x

 

y0(n1)

 

(x x )n1 +

y0(n2)

 

(x x )n2

 

y =

dx

dx....

f (x)dx +

 

 

+...+

 

 

(n 2)!

 

 

 

 

(n 1)!

0

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x0

 

x0

 

x0

 

 

 

 

 

 

 

 

 

+ y0(x x0 )+ y0.

 

 

 

 

 

 

 

 

41

Пример. Найти максимальный прогиб шарнирно опëртого стержня при действии на него равномерно распределëнной поп е- речной нагрузки интенсивности q(x) = q0 (см. рис. 5б).

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

M (x) = q20l x q02x2 ,

где первое слагаемое – момент от опорной реакции, а второе – момент от равномерно распределëнной поперечной нагрузки. Поэтому уравнение изгиба стержня (2.2) примет вид

 

1

 

 

 

2

y′′ =

q0l x q0x

.

 

 

 

 

2

2

 

 

EI

 

Это уравнение относится к типу уравнений (2.5). Интегрируя его дважды, получим общее решение в виде

 

q0

 

 

3

 

x

4

 

 

y =

lx

 

 

+C x +C

.

 

 

 

 

 

 

 

6

12

1

2

 

2EI

 

 

Произвольные постоянные находим из граничных условий (2.4):

 

C

= −

 

l3

,

 

C

 

= 0.

 

 

 

 

12

 

 

 

 

 

 

1

 

 

 

 

 

2

 

 

 

 

 

 

 

Подставляя их в общее решение, получим

 

 

 

y =

q

l4

x 3

 

x

4

x

(2.7)

0

 

 

2

 

 

 

 

 

 

.

 

 

 

 

 

24EI

l

 

 

l

 

l

 

Максимальный прогиб стержня при действии равномерно

распределëнной нагрузки при

x = l / 2

 

будет равен

 

 

 

 

 

y

 

= −

5q

l4

 

.

 

 

 

 

 

 

 

max

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

384EI

 

 

 

2.

y

(n)

=

 

 

 

 

′′

,

,y

(n1)

).

(2.8)

 

f (x,y ,y

 

 

Уравнение (2.8) не содержит искомой функции y(x).

Рассмотрим процедуру интегрирования уравнения данного типа на примере уравнения второго порядка

y

′′

= f (x,y ).

(2.9)

 

 

 

42

 

Понижение порядка достигается подстановкой

(2.10)

y (x) = z(x).

 

 

Тогда

 

(2.11)

y (x) = z (x),

′′

 

и уравнение (2.9) приводится к уравнению первого

порядка от-

носительно функции z(x):

 

 

 

(2.12)

z (x) = f (x,z).

Интегрируя уравнение (2.12), находим его общий интеграл в

виде

z(x) = ϕ(x,C1 ),

(2.13)

где С1 - произвольная постоянная. Далее в (2.13) заменяем левую

часть согласно (2.10). и вновь получаем уравнение первого порядка относительно искомой функции y:

y (x) = ϕ(x,C1 ),

(2.14)

 

Интегрируя уравнение (2.14), находим общее решение исходного уравнения (2.9) в виде

y(x) = ξ(x,C1 ,C2 ).

(2.15)

Замечание. Если уравнение (2.9) не содержит ни искомой функции y, ни её производных до ( k- 1) - го порядка включительно, то есть имеет вид

y(n) = f (x, y(k ), y(k +1), , y(n1) ),

то его порядок может быть понижен сразу на k единиц подстановкой

y(k ) = z(x).

Пример. Решить уравнение x3 yI V + x2 y′′′ =1.

 

Это уравнение не

содержит явно искомой функции y(x) и

её первых производных

y ,

y

′′

и относится ко второму из рас-

 

 

 

 

смотренных нами типов. Применяя подстановку

(2.16)

 

y

 

(x) = z(x),

 

 

′′′

 

 

получаем линейное неоднородное уравнение первого порядка

43

x3 z′+ x2z =1

или

1

1

 

 

z′+

 

 

x z =

 

.

(2.17)

x3

Интегрируем его методом вариации произвольной постоянной. Сначала решаем соответствующее однородное уравнение

z′+

1

 

x z = 0.

(2.18)

Разделяем в нëм переменные:

dzz = −dxx .

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

ln z = −ln x + lnC1 = ln Cx1 .

Следовательно, общее решение уравнения (2.18) будет

z =

C1 .

(2.19)

 

x

 

Далее ищем решение уравнения (2.17) в форме, аналогичной по структуре выражению (2.19), но произвольную постоянную в (2.19) заменяем неизвестной функцией v(x):

 

 

 

 

 

z =

v(x)

.

 

 

(2.20)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

Подставляя (2.20) в (2.17), получим

 

 

 

 

 

 

 

 

 

 

1

 

 

 

v (x)x v(x) + v(x)

=

,

 

dv(x)

 

1

x2

 

x2

 

x3

 

откуда следует

=

и,

после разделения переменных,

dx

x2

 

 

 

 

 

 

 

 

 

dv(x) = dxx2 .

Интегрируя это уравнение, находим v(x) = −1/ x +C1. Поэтому согласно (2.20) имеем

z = −

1

+ C1 .

(2.21)

x2

 

x

 

44

Заменяя в выражении (2.21) z по формуле (2.16), приходим к уравнению третьего порядка относительно искомой функции y(x):

y

′′′

1

C1

 

 

 

+ x .

 

= − x2

(2.22)

 

Уравнение (2.22) содержит

в правой

части известную

функцию от x и относится к первому из рассмотренных нами типов (см. (2.5)). Интегрируя его последовательно три раза, окончательно получим общее решение исходного уравнения, содержащее 4 произвольных постоянных :

 

 

 

 

y′′ = 1 + C1 ln x + C2,

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

y′ = ln x +C1x(lnx 1)+C2x +C3 ,

 

 

 

 

 

 

 

x2

 

3

 

 

 

 

x2

 

 

 

 

y = x(lnx 1) + C

 

 

ln x

 

 

+ C

 

 

+ C

x + C

.

 

 

 

2

2 2

 

 

 

 

1 2

 

 

 

 

 

 

3

4

 

3.

y

(n)

=

f (y,

 

′′

 

y

(n1)

).

 

 

 

(2.23)

 

y , y

, ,

 

 

 

 

Это уравнение не содержит явно независимой переменной x. В частном случае уравнение второго порядка данного типа будет

y

′′

= f (y,y ).

(2.24)

 

 

Понижение порядка достигается подстановкой

(2.25)

y (x) = z(y).

 

 

 

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

 

′′

dz dy

dz

 

y

(x) = dy dx

= dy z .

(2.26)

Поэтому уравнение (2.24) приводится к уравнению первого порядка относительно функции z(y):

z dz

= f (y,z).

(2.27)

dy

 

 

Интегрируя уравнение (2.27), находим его общее решение в виде

z(y) =ϕ(y,C1 ),

(2.28)

где C1 - произвольная постоянная. Далее в (2.28) заменяем левую часть согласно (2.25) и вновь приходим к уравнению первого по-

45

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

функции y : y (x) =ϕ (y,C1 ), следовательно,

 

 

 

 

dy

= dx.

(2.29)

ϕ(y,C )

 

 

1

 

 

Интегрируя уравнение (2.29), окончательно получим общий интеграл исходного уравнения (2.24) в виде

 

 

 

 

 

y(x) = ξ(x,C1 ,C2 ).

 

(2.30)

Пример. Решить задачу Коши:

 

 

 

 

 

 

 

 

 

16y

3

y

′′

= 4 y

4

1,

y(0) =

 

2

 

 

 

2

 

 

 

 

2

,

2 .

 

 

 

y (0) =

Уравнение не содержит явно независимой переменной x и потому относится к третьему из перечисленных типов. Принимаем y(x) = z(y), тогда с учетом (2.26) получим уравнение первого

порядка с разделяющимися переменными

16y3 z dydz = 4 y4 1.

Приводим его к виду zdz = 4 y4 3 1 dy. Интегрируя, получим

16y

 

z2

 

y2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

y2

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

=

 

+

 

 

 

+C

 

,

z = ±

 

 

 

 

 

+

 

 

 

 

 

 

+C

, где C

= 2C

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

8

 

16y2

1

 

 

 

 

 

 

 

4

 

 

 

 

16y2

 

 

 

1

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Заменяя z

 

на y

согласно (2.25),

 

 

 

вновь получаем уравнение

первого порядка

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y′ = ±

4 y4

+16C y2 +

1

 

.

 

 

 

(2.31)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4 y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Определим произвольную постоянную C1

из начальных ус-

ловий. При x = 0 равенство (2.31) примет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

1

+16

1 C +1

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= ±

 

 

 

 

 

 

 

 

 

2

1

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда получим 2 = ±2+8C1 , C1 = 14 . Подставляя C1 в (2.31), приходим к уравнению

 

 

 

 

 

 

 

 

 

46

 

 

 

 

 

 

 

2y2 +1

 

y

 

 

4 y4

+ 4 y2 +1

 

 

 

=

 

 

 

 

 

= 4 y (знак минус отброшен, так как при

 

 

4 y

 

 

 

этом не удовлетворяются начальные условия).

Разделяя переменные, имеем

4 ydy

 

= dx , следовательно, после

2y2 +1

 

 

интегрирования получим ln(2y2 +1) = x +ln C2, 2y2 +1 = C2ex . Из первого начального условия находим C2 = 2.

В результате решение задачи Коши запишется в виде y = ex 12.

Пример. Решить задачу Коши: y

′′

= 2y

3

,

y(0) =1,

 

 

 

y (0) =1.

Полагая, как и выше,

 

 

y

′′

dz

получим

y (x) = z(y),

 

(x) = dy z ,

уравнение первого порядка

z dz = 2y3

.

Разделяя переменные и

 

dy

 

 

 

 

 

 

 

 

 

 

интегрируя, получим z = ±

 

 

 

 

 

 

y′ =

 

(знак

y4 +C1

.

 

 

 

 

 

y4 +C1

минус отброшен, так как при этом не удовлетворяются начальные условия).

Используя начальные условия, находим C1 = 0. Поэтому y′ = y2 , и уравнение легко интегрируется: 1y = x +C2. Далее из первого

начального условия находим C2 = −1. Поэтому решение задачи Коши запишется в виде

y =11 x .

Замечание. При решении задачи Коши часто целесообразно (как в приведëнных примерах) определять значения произвол ь- ных постоянных в процессе решения, а не после нахождения общего решения исходного уравнения, так как при произвольных значениях C1,C2,...интегрирование может оказаться весьма слож-

ным, а то и вообще невозможным. Действительно, например, в последней задаче решение уравнения y′ = y4 +C1 приводит по-