Главная > ТЕОРЕТИЧЕСКАЯ МЕХАНИКА (А.П.Маркеев)
НАПИШУ ВСЁ ЧТО ЗАДАЛИ
СЕКРЕТНЫЙ БОТ В ТЕЛЕГЕ
<< Предыдущий параграф Следующий параграф >>
Пред.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
След.
Макеты страниц

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

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

ДЛЯ СТУДЕНТОВ И ШКОЛЬНИКОВ ЕСТЬ
ZADANIA.TO

91. Уравнение движения. Определение реакций. Рассмотрим твердое тело, имеющее две неподвижные точки O и O1 (рис. 91).

Пусть F и F1 — реакции связи в точках O и O1,R — главный вектор активных сил, а MO — главный их момент относительно точки O.

Примем точку O за начало неподвижной системы координат OXYZ, ось OZ которой направим по оси OO1. С телом жестко свяжем систему координат Oxyz; ее ось Oz направлена вдоль оси OO1. Тело имеет одну степень свободы. За обобщенную координату примем угол φ, образуемый осями OX и Ox.

Чтобы получить уравнения, описывающие движение твердого тела, воспользуемся теоремами об изменении количества движения и момента количеств движения. Обе части соответствующих этим теоремам уравнения (5) п. 86 и уравнения (8) п. 87 спроектируем на оси вращающейся системы координат Охуz. Для этого воспользуемся формулой (5) п. 30 , связывающей абсолютную и относительную производные вектора. Получим:
Md~vCdt+Mω×vC=R+F+F1,d~KOdt+ω×KO=MO+OO1×F1.

Здесь M — масса тела, ω — его угловая скорость, vC — скорость центра масс тела; реакция F не вошла в уравнение (2), так как ее момент

относительно точки O равен нулю. Пусть в связанной с телом системе координат
R=RxRyRz,MO=MxMyMz,F=FxFyFz,F1=F1xF1yF1z,ω=pqr,OC=xCyCzC,KO=KxKyKz.

Очевидно, что p=0,q=0,r=φ˙, и из формул (8) п. 82 следует, что Kx=Jxzφ˙,Ky=Jyzφ˙,Kz=Jzφ˙, где Jxz, Jyz — центробежные, а Jz — осевой моменты инерции тела для точки O.

Принимая во внимание, что vC=ω×OC и обозначая буквой h расстояние между неподвижными точками тела O и O1, получим скалярную форму векторных уравнении (1) и (2) в виде
MyCφ¨MxCφ˙2=Rx+Fx+F1x,MxCφ¨MyCφ˙2=Ry+Fy+F1y,0=Rz+Fz+F1z,Jxzφ¨+Jyzφ˙2=MxhF1y,Jyzφ¨Jxzφ˙2=My+hF1x,Jzφ¨=Mz.

Последнее уравнение не содержит реакций и является дифференциальным уравнением вращения твердого тела вокруг неподвижной оси. Остальные пять уравнений служат для нахождения реакций. Последняя задача является неопределенной. Действительно, из третьего уравнения системы (3) видно, что нельзя отдельно найти продольные реакции Fz и F1z, а можно определить лишь их сумму. Эта сумма не зависит от характера вращательного движения тела. Поперечные реакции Fx,F1x,Fy,F1y находятся из первого, второго, четвертого и пятого уравнений
Рис. 92
системы (3); они зависят от вращения тела.
ПРимер 1. Равнобедренный прямоугольный треугольпиик OO1A вращается вокруг вертикальной оси, к которой он подвешен катетом

OO1=a (рис. 92). Какова должна быть угловая скорость вращения, чтобы боковое давление на нижнюю опору O равнялось нулю? Треугольник считать тонкой однородной пластинкой.

Для решения задачи воспользуемся уравнениями (3). В рассматриваемом случае
xC=0,yC=a/3,h=a,
Jxz=0, а для Jyz получаем
Jyz=yzdm=2ma20az(0zydy)dz==ma20az3dz=14ma2.

Далее,
Rx=Ry=0,Rz=mg,Mx=13mga,My=Mz=0.

Учитывая еще, что по условию задачи Fx=Fy=0, получаем уравнения (3) в виде
13maφ¨=F1x,13maφ˙2=F1y,0=mg+Fz+F1z,14ma2φ˙2=13mgaaF1y,14ma2φ¨=aF1x,Jzφ¨=0.

Из последнего уравнения следует, что φ˙=ω= const, m. е. вращение треугольника происходит с постоянной угловой скоростью.

Исключив из второго и четвертого уравнений величину F1y, придем к соотношению, определяющему величину угловой скорости. Окончательно найдем, что
ω=2g/a.
92. Условия, при которых динамические реакции равны статическим. Если в первом, втором, четвертом и пятом уравнениях системы (3) положить φ˙=0,φ¨=0, то получим систему уравнений для определения поперечных статических реакций. Если же тело вращается, то либо φ˙, либо φ¨, либо та и другая из этих величин не будут равными нулю. Поэтому левые части упомянутых уравнений в общем случае не будут тождественно равными нулю во все время движения и, следовательно, динамические реакции отличаются от статических.

Найдем условия, при которых динамические реакции равны статическим. Приравнивая нулю левые части первого, второго, четвертого и пятого уравнений системы (3), получим следующие две пары равенств:
yCφ¨+xCφ˙2=0,yCφ˙2+xCφ¨=0.Jxzφ¨Jyzφ˙2=0,Jxzφ˙2+Jyzφ¨=0.

Равенства (4) и (5) можно рассматривать как однородные системы линейных уравнений соответственно относительно xC,yC и Jxz,Jyz. Определители этих систем одинаковы и равны величине φ¨2+φ˙4. Если тело вращается, то эта величина не может быть тождественно равной нулю. Поэтому равенства (4) и (5) удовлетворяются только при выполнении условий
xC=yC=0,Jxz=Jyz=0.

Таким образом, динамические реакции при вращении твердого тела вокруг неподвижной оси равны статическим тогда и только тогда, когда ось вращения является главной центральной осью инерции тела.
93. Уравнение движения физического маятника. Физическим маятником называется твердое тело, которое может вращаться вокруг неподвижной горизонтальной оси под действием силы тяжести. Выберем неподвижную систему координат OXYZ так, чтобы ее ось OZ совпадала с осью вращения маятника, а ось OY была направлена вертикально вниз. Связанную с маятником систему координат Oxyz выберем так, чтобы центр масс маятника лежал на оси Oy, а оси Oz и OZ совпадали. Тогда если a — расстояние от центра тяжести до оси вращения, то Mz=mgasinφ, и из последнего уравнения системы (3) получим дифференциальное уравнение движения физического маятника в виде
φ¨+mgaJzsinφ=0.

Сравнивая это уравнение с уравнением движения плоского математического маятника, которое задается равенством (6) п. 57 , находим, что физический маятник будет колебаться по такому же закону, что и математический маятник длиной
l=Jzma.

Величину l, определяемую по формуле (7), называют приведенной длиной физического маятника.
94. Фазовая плоскость для уравнения движения маятника. Для выяснения общих свойств движения систем с одной степенью

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

Будем считать, что правая часть этого уравнения удовлетворяет условию существования и единственности его решения.

Можно считать, что уравнение (8) описывает движение в системе с одной степенью свободы, в которой x играет роль обобщенной координаты, а кинетическая и потенциальная энергия определены равенствами
T=12x˙2,Π=f(x)dx.

Полная механическая энергия E=T+Π постоянна во все время движения, т. е. уравнение (8) имеет первый интеграл
E(x,x˙)12x˙2+Π(x)=h= const. 

Уравнение (8) эквивалентно системе двух уравнений
x˙=y,y˙=f(x).

Плоскость с координатами x,y называется фазовой плоскостью уравнения (8). Точки фазовой плоскости называются фазовыми точками. В каждой точке плоскости, где определена функция f(x), система (10) задает вектор с компонентами x˙,y˙; этот вектор называется ϕa зовой скоростью. Решение системы (10) задает движение фазовой точки по фазовой плоскости, причем скорость движения фазовой точки равна фазовой скорости в том месте плоскости, где в данный момент находится точка. Кривая, которую описывает фазовая точка, называется фазовой кривой. В частных случаях фазовая кривая может состоять из одной точки. Такие точки называются положениями равновесия. Вектор фазовой скорости в положении равновесия равен нулю.

Интеграл (9) позволяет легко находить фазовые кривые. На каждой фазовой кривой значение полной механической энергии E постоянно, поэтому каждая фазовая кривая целиком принадлежит одному уровню энергии E(x,x˙)=h.
Запишем интеграл (9) в виде
x˙2=2(hΠ(x)).

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

1. При данном h фазовые кривые могут располагаться только в той части фазовой плоскости, где выполняется неравенство Π(x)h. Эту часть плоскости называют областью возможности движения. Неравенство Π(x)h следует из того, что для реальных движений правая часть формулы (11) не может быть отрицательной, так как ее левая часть есть квадрат вещественной величины.
2. Как видно из системы (10), положения равновесия лежат на оси x фазовой плоскости, причем в положении равновесия x=x, где x критическая точка потенциальной энергии, т. е. такая точка, в которой dΠ/dx=0.
3. Если x=x является точкой локального минимума функции Π(x), причем d2Π/dx2>0 при x=x, то точка (x,0) на фазовой плоскости будет особой точкой типа центр для системы (10). Если же x=x — точка локального максимума и в ней d2Π/dx2<0, то (x,0) — особая точка типа седло.
4. Совокупность фазовых кривых симметрична относительно оси x. Это свойство вытекает из формулы (11) ввиду четности входящих в нее величин относительно x˙.
5. В точках оси x, отличных от положений равновесия, фазовые кривые ортогональны оси x. Это свойство сразу видно из системы (10), так как в этих точках x˙=0,y˙=f(x)eq0.
Перечисленные свойства позволяют по виду графика функции Π(x) сразу делать выводы о характере движения, описываемого уравнением (8). На рис. 93 для примера показаны график потенциальной энергии и соответствующие фазовые кривые. Направления движения фазовой точки показаны стрелками. При h=h1 есть положение равновесия типа центр. Это положение равновесия окружено замкнутыми фазовыми кривыми. При h>h3 фазовые кривые разомкнуты. На уровне энергии h=h3 есть положение равновесия типа седло. На уровне h=h3 расположена также фазовая кривая, начинающаяся вблизи седловой точки и при t входящая в эту точку. Эта кривая отделяет область, где расположены замкнутые кривые, окружающие точку типа центр, от разомкнутых фазовых кривых, соответствующих значениям h, большим h3. Такие кривые, которые раз-

деляют области с различным характером поведения фазовых кривых, называют сепаратрисами.

Построим фазовую плоскость для дифференциального уравнения (6), описывающего движение маятника. Для кинетической и потенциальной энергии маятника имеем выражения
T=12Jzφ˙2,Π=mgacosφ.

Если положить ω02=g/l,Π=ω02cosφ, то интеграл энергии T+Π= const запишется в виде
12φ˙2+Π=h= const. 

Рис. 94
График функции Π(φ) и фазовые кривые представлены на рис. 94. Картина фазовых кривых периодична по φ с периодом 2π. При h<ω02 движение невозможно. При h=ω02 маятник находится в положении равновесия, когда его центр масс занимает самое низкое из возможных положений. На фазовой плоскости φ,φ˙ этому положению равновесия соответствуют точки, в которых φ=2kπ(k=0,±1,±2,), а φ˙=0. Это точки типа центр. Они окружены замкнутыми фазовыми кривыми, соответствующими колебаниям маятника. Колебательным движениям маятника соответствуют значения h, удовлетворяющие неравенству ω02<h<ω02.

При h=ω02 возможны два типа движений. Один соответствует положению равновесия маятника, когда его центр масс занимает наивысшее возможное положение. Этому равновесию на фазовой плоскости соответствуют точки φ=π+2kπ(k=0,±1,±2,),φ˙=0. Это точки типа седло. Для другого типа движений при h=ω02 центр масс

маятника асимптотически при t стремится занять наивысшее положение. Асимптотическим движениям на фазовой плоскости соответствуют кривые, соединяющие точки типа седло. Эти кривые являются сепаратрисами.

При h>ω02 движение маятника будет вращательным. Для этого движения абсолютная величина угла φ монотонно возрастает. Этим движениям на фазовой плоскости соответствуют незамкнутые кривые. Сепаратрисы разделяют области колебательных и вращательных движений.
95. Некоторые сведения из теории эллиптических интегралов и эллиптических функций Якоби. В этой главе и в некоторых других разделах книги будут использоваться так называемые эллиптические интегралы и эллиптические функции. Дадим здесь необходимые определения и понятия.
Интеграл
u=F(φ,k)=0φdx1k2sin2x

называется эллиптическим интегралом первого рода. Величина k называется модулем эллиптического интеграла. Обычно считается, что k удовлетворяет неравенствам 0k<1.
Интеграл
E(φ,k)=0φ1k2sin2xdx

называется эллиптическим интегралом второго рода.
Величина
K(k)=F(π2,k)=0π2dx1k2sin2x

называется полным эллиптическим интегралом первого рода, а величина
E(k)=E(π2,k)=0π21k2sin2xdx
— полным эллиптическим интегралом второго рода.

При малых k интегралы (15), (16) представляются в виде сходящихся рядов по степеням k :
K=π2(1+14k2+964k4+),E=π2(114k2364k4).

Из (13) и (14) можно получить следующее выражение для производных эллиптических интегралов по модулю k :
F(φ,k)k=1k2(E(φ,k)k2F(φ,k)kksinφcosφ1k2sin2φ),E(φ,k)k=E(φ,k)F(φ,k)k,

где k2=1k2,k дополнительный модуль.
Если в равнствах (19) и (20) положить φ=π2, то получим производные по k от полных эллиптических интегралов (15) и (16):
dKdk=E(k)k2K(k)kk2,dEdk=E(k)K(k)k.

Функция, являющаяся результатом обращения эллиптического интеграла первого рода, называется амплитудой и обозначается так:
φ=amu.

Функции z=sn(u,k) (эллиптический синус) и z=cn(u,k) (эллиптический косинус) определяются так:
z=sn(u,k)=sinφ=sinamu и z=cn(u,k)=cosφ=cosamu.

Так как sinφ и cosφ имеют период 2π по φ, то согласно (13) и (15), эллиптические синус и косинус имеют по u период, равный 4K(k). Функция дельта амплитуды z=dn(u,k) определяется так:
z=dn(u,k)=dφdu=1k2sin2φ=1k2sn2(u,k).

Функция дельта амплитуды имеет период 2K(k) по u. Функции φ=amu,z=sn(u,k), z=cn(u,k),z=dn(u,k) аналитичны относительно k и при k0 стремятся соответственно к функциям φ=u,z=sinu, z=cosu,z=1.

Эллиптические функции Якоби удовлетворяют следующим легко проверяемым тождествам:
Рис. 95
sn2u+cn2u=1,dn2u+k2sn2u=1.

Справедливы следующие формулы дифференцирования эллиптических функций:
ddusnu=cnudnu,dducnu=snudnu,ddudnu=k2snucnu.

Графики эллиптических функций Якоби представлены на рис. 95. 96. Интегрирование уравнения движения маятника. Рассмотрим три случая в соответствии с возможными значениями константы h в интеграле (12).
1. ω02<h<ω02. В этом случае, как показано в п. 94 , маятник совершает колебания. Пусть β — максимальный угол, на который отклоняется маятник от своего вертикального положения, соответствующего значению φ, равному нулю. Тогда h=ω02cosβ и интеграл (12) запишется в виде
φ˙2=2ω02(cosφcosβ).

Положим k1=sin(β/2) и сделаем замену переменных
sin(φ/2)=k1sinψ.

Тогда интеграл энергии (27) примет следующий вид:
ψ˙2=ω02(1k12sin2ψ).

Если принять, что при t=0φ=0, то отсюда получаем
ω0t=0ψdx1k12sin2x=F(ψ,k1),
т. е. ψ=am(ω0t). Поэтому из (23) и (28) имеем окончательно
φ=2arcsin(k1snω0t).

Функция φ периодична по t с периодом τ, вычисляемым, согласно п. 95 , по формуле
τ=4K(k1)/ω0.

Воспользовавшись разложением (17), получим, что при небольшом значении максимального угла отклонения β маятника от вертикали имеем приближенное значение периода
τ=2πl/g,

что совпадает с известным значением периода малых колебаний маятника. При учете двух первых членов разложения τ в ряд по β получаем более точное значение периода
τ=2πlg(1+β216).

Если βπ, то k11 и период колебаний τ неограниченно возрастает.
2. h>ω02. В этом случае маятник находится в режиме вращения. Пусть при t=0φ=0,φ˙=φ˙0. Тогда h=12φ˙02ω02 и интеграл (12) запишется в виде
φ˙2=φ˙02(1k22sin2φ2),

где введено обозначение
k22=4ω02φ˙02.

Так как h>ω02, то φ˙02>4ω02 и, следовательно, k22<1. Из (35) имеем
φ˙02t=F(φ2,k2)=0φ/2dx1k22sin2x.

Следовательно,
φ=2am(φ˙0t/2).

Если начальная угловая скорость велика, т. е. φ˙02ω02, то приближенно φ=φ˙0t и вращение маятника мало отличается от равномерного.
3. h=ω02. Этот случай соответствует асимптотическим движениям маятника. Интеграл энергии (12) в этом случае дает соотношение
φ˙2=4ω02cos2(φ/2).

Если при t=0φ=0 и φ˙>0, то отсюда после интегрирования получаем
φ=π+4arctg(eω0t).

1
Оглавление
email@scask.ru