Материал: 864

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

данными: при t 0 (в начальный момент)

v v0 . Подставляя началь-

ные данные в общее решение (5.1), имеем

 

v0 0 C1, или C1

v0 .

Таким образом, уравнение (5.1) принимает вид

v at v0 .

(5.2)

Уравнение (5.2) представляет закон движения в виде v f (t).

Так как v dS , то решение задачи Коши (5.2) можно преобразовать к dt

виду

 

dS

 

at v ,

 

 

 

 

 

 

 

 

 

dt

 

 

0

 

 

 

или

 

 

 

 

 

 

dS atdt v0dt.

 

 

(5.3)

Интегрируя равенство (5.3), получаем общее решение задачи:

S

 

1

at2 v t C

2

.

(5.4)

 

 

 

 

2

0

 

 

 

 

 

 

 

 

Для определения постоянной интегрирования C2 используем начальные условия: при t=0 S S0 . Подставляя эти значения в равенство (5.4), получаем

S0 0 0 C2,

или C2=S0. Таким образом, уравнение движения (5.4) принимает вид

S

1

at2

v t S

0

.

 

(5.5)

 

 

2

 

0

 

 

 

 

 

 

 

 

 

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

S f (t).

положить a g;

 

 

Если в уравнениях (5.2) и (5.5)

v0 0;

S0 0;

Sh,то получаем закон свободного падения тела в пустоте:

vgt; h 1 gt2 .

2

В разделе 4 пособия была рассмотрена нормальная система дифференциальных уравнений (4.1). С механической точки зрения эта система задает поле скоростей движений, где t – время; x1,x2,...,xn прямоугольные координаты точки фазового пространства.

Решение системы есть функции

51

x1 1(t,t0,x10,x20,...,xn0);

........................................ (5.6)

xn n(t,t0,x10,x20,...,xn0),

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

(x10,x20,...,xn0). При такой интерпретации система (4.1) называется динамической, а каждое ее решение (5.6) – движением. Кривая, описываемая точкой при движении, называется траекторией движения.

Уравнения (5.6) задают траекторию движения в параметрическом виде, t – параметр. Отметим, что при n = 2 имеем движение в плоскости (x1,x2), при n = 3 имеем движение в пространстве

(x1,x2,x3).

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

Пример 37. Задача о движении материальной точки под действием силы F, рассматриваемая в теоретической механике, сводится к системе трех дифференциальных уравнений второго порядка. Обозначим через Fx,Fy,Fz проекции силы F на оси координат. Положе-

ние точки в любой момент времени t определяется ее координатами x, y, z, которые являются функциями от времени t; проекции вектора

скорости на оси координат есть dx,dy ,dz (в теоретической механике dt dt dt

их обычно обозначают через x, y,z).

Пусть сила F (а следовательно, и ее проекции) зависит от времени t, положения x, y,z точки и скорости движения точки. Тогда математической моделью поступательного движения материальной точки массы m, согласно закону Ньютона, является система трех дифференциальных уравнений второго порядка:

 

m

d2x

 

 

F

(t,x, y,z,

 

 

dx

 

 

dy

 

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

,

 

 

 

);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt2

 

 

x

 

 

 

dt

 

 

dt

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d2 y

 

 

 

 

 

 

dx

 

 

dy

 

 

 

dz

 

m

 

 

 

 

 

 

Fy

(t,x, y,z,

 

 

,

 

 

 

 

,

 

 

 

 

);

(5.7)

 

 

dt2

 

dt

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

m

d2z

 

 

F (t,x, y,z,

dx

 

,

dy

 

,

dz

).

 

 

 

 

 

 

 

 

 

 

 

 

dt

2

 

 

z

 

 

 

dt

 

 

dt

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

52

Решением этой системы являются функции x(t,C10,C20); y(t,C10,C20); z(t,C10,C20), которые дают траекторию движения точки в пространстве R3, константы C10,C20 определяются из начальных условий.

Очевидно, что в случае плоского движения, когда траекторией является кривая в плоскости XOY, имеем систему

 

 

 

mx

Fx(t,x, y,x, y);

(5.8)

 

 

my

Fy(t,x, y,x, y).

 

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

Получим систему вида (5.8) для решения конкретной задачи: тело брошено под углом к горизонту и движется в среде, сопротивление которой пропорционально скорости v тела; найти траекторию движения.

Решение. В любой точке M(x, y) траектории на тело действуют две силы: тяжести P mg и сопротивления среды F kv(рис. 1).

Y

F M(x,y)

α

P

X

.

Рис. 1

Найдем проекции на оси координат равнодействующей этих сил

Fp :

 

 

 

^

OX);

Fpx

Pcos(P,

OX) F cos(F,

 

 

 

 

 

 

(5.9)

 

 

 

 

 

Fpy

Pcos(P,

OY) F cos(F,

 

OY).

 

 

53

Так как cos(P, OX) 0; cos(F, OX) dx ; ds

cos(F, OY) dy , система (5.9) принимает вид ds

dx

;Fpx kv

 

ds

Fpy

mg kv

dy

.

 

 

 

ds

cos(P, OY) 1;

(5.10)

Поскольку v ds , то систему (5.10) можно записать в виде dt

dx

;Fpx k

 

 

dy

dt

Fpy

k

mg.

 

 

 

dt

Тогда, согласно второму закону Ньютона [см. формулу (5.8)], имеем дифференциальные уравнения движения

 

mx kx;

 

 

(5.11)

my ky

mg.

Каждое из дифференциальных уравнений, входящих в систему (5.11), содержит только одну неизвестную функцию, поэтому интегрирование системы сводится к интегрированию по отдельности каждого из уравнений системы. Первое уравнение системы есть линейное однородное дифференциальное уравнение с постоянными коэффициентами, второе – линейное неоднородное дифференциальное уравнение с постоянными коэффициентами. Уравнения такого типа были рассмотрены в подразделе 3.3, поэтому приводим решение уравнений

системы (5.11) без особых комментариев.

m

d2x

k

dx

0, следова-

dt2

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

x C1x1

C2x2; со-

ставляем характеристическое уравнение:

 

 

 

 

 

 

 

 

 

k

 

m 2 k 0 (m k) 0 0;

k

 

 

 

 

 

 

t

 

 

 

 

 

 

x 1;x

2

e m .

 

 

 

1

 

 

2

 

 

m

1

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

Тогда функция x C C

 

есть

общее решение первого

e m

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

уравнения системы (5.11).

Решаем второе уравнение системы (5.11):

54

m

d2 y

k

dy

mg

d2 y

 

k dy

g.

 

 

 

 

 

 

 

 

dt2

dt

dt2

m dt

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

y ;

 

 

 

 

 

 

 

 

 

t

 

ищем в виде y

y

 

 

 

 

 

 

 

 

 

 

y

C

 

C

e m

;

 

 

 

 

 

 

 

 

 

3

 

 

 

4

 

 

 

 

Общее решение уравнения

y At;(y ) A;(y ) 0.

Подставляем в уравнение y ,(y ) ,(y ) и находим значение A:

 

k

A g A

gm

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

y

gm

t; общее решение

 

m

 

 

 

 

k

 

k

 

 

 

 

 

 

 

 

 

k

t

 

 

 

 

 

 

 

 

 

 

второго уравнения системы (5.11) есть функция y C3 C4e m

 

gmt. Таким образом, получили общее решение системы (5.11): k

 

 

 

 

 

 

k

t

 

 

x(t) C C

 

 

 

;

 

 

e m

 

 

1

2

 

 

 

 

 

 

 

(5.12)

 

k

 

 

 

gm

 

 

 

 

t

 

 

 

t.

 

 

 

 

 

m

 

 

 

 

 

y(t) C3 C4e

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

Найдем частное решение системы (5.11), используя начальные условия: x(0) 0; y(0) 0;x (0) v0 cos ; y (0) v0 sin .

Подставим первые два из начальных условий в общее решение (5.12), получим систему

0 C1 C2;

(5.13)

0 C3 C4.

Из системы (5.12) имеем, дифференцируя оба уравнения систе-

мы:

k

dx k C2e m ;

dt m t

 

 

 

 

k

 

 

 

 

dy

 

k

 

t

 

gm

 

 

 

 

 

 

C4e m

 

 

.

 

m

 

dt

 

 

 

 

 

k

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

 

v0 cos

 

k

C2;

 

 

 

 

 

m

 

 

k

 

 

(5.14)

v0 sin

 

C4

 

gm

.

m

 

 

 

 

 

 

 

k

Из систем (5.13) и (5.14) находим произвольные постоянные:

55