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

Основы механики твердого деформируемого тела

.pdf
Скачиваний:
128
Добавлен:
14.05.2015
Размер:
26.55 Mб
Скачать

ГЛАВА 1. СИСТЕМЫ С ОДНОЙ СТЕПЕНЬЮ СВОБОДЫ

1.1. Собственные колебания системы без учета сил сопротивления.

Пусть на невесомой конструкции сосредоточена точечная масса m1, движение которой определяется перемещением u1(t). Кроме того, имеется заданная сила P (t), возбуждающая колебания. Предполагается, что ни сопротивление среды, ни трение в соединениях, ни так называемое внутреннее трение (т. е. рассеивание механической энергии из-за необратимых деформаций материала) не оказывают существенного влияния на колебательный процесс. Требуется описать движение данной системы.

Решение задачи начинается с того, что по направлению движения массы

прикладывается сила инерции P1

(принцип Д’Аламбера), равная произведе-

нию массы m1 на ускорение

..

, взятому с обратным знаком:

..

u1

P1 = −m1u1.

Затем используется принцип наложения, который позволяет представить искомое перемещение в виде:

u1 = u11P1 + u10P (t).

Обозначения для единичных перемещений u11 и u10 обычные. Иллюстрацией к сказанному служит рис. 1.1.

Полученное уравнение после подстановки в него

силы P1 записывается следующим образом:

 

..

+ ω

2

u1

= ω

2

u10P (t),

(1.1)

u1

 

 

где

 

 

 

1

 

 

 

 

ω2 =

 

.

(1.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m1u11

 

Линейное дифференциальное уравнение (1.1) должно интегрироваться с учетом двух начальных условий. Если возмущающая сила P (t) отсутствует, то система совершает свободные (собственные) колебания. Чтобы возбудить их, надо вывести систему из состояния покоя, а затем предоставить ее самой себе. В этом случае уравнение (1.1) становится однородным:

..

2

u1

= 0,

(1.1a)

u1

+ ω

а его решением являются гармонические колебания

u1 = A sin ωt + B cos ωt.

(1.3)

582

Часть V

Следовательно, величина (1.2) представляет собой квадрат круговой частоты собственных колебаний системы с одной степенью свободы.

Постоянные A и B, определяющие форму колебаний (1.3), находят по начальным условиям задачи. Систему можно вывести из равновесия, если отвести массу m1 на расстояние u0, а затем плавно отпустить, либо сообщить неподвижной массе импульс, т. е. задать начальную скорость v0. Задание начальных условий и состоит в указании величин u0 и v0 в некоторый начальный момент времени t0:

A sin ωt0 + B cos ωt0 = u0,

ωA cos ωt0 − ωB sin ωt0 = v0.

Отсюда следует

 

 

 

 

 

 

 

 

 

A = u0 sin ωt0 +

v0

 

 

B = u0 cos ωt0

v0

 

(1.4)

 

cos ωt0,

 

 

 

sin ωt0.

ω

 

ω

Если же t0 = 0, то

A = v0/ω,

B = u0.

 

 

 

(1.4a)

 

 

 

 

Теперь можно найти силу инерции:

 

 

 

 

 

..

= m1ω

2

(A sin ωt + B cos ωt).

 

(1.5)

P1 = −m1u1

 

 

Таким образом, сила P1 меняется во времени по тому же закону, что и

перемещение (1.3), т. е.

 

 

 

 

 

 

 

 

 

 

 

 

P1 = m1ω2u1.

 

 

 

(1.5a)

Интерес представляет амплитудное значение этой силы, определяемое в тот

.

.

= 0:

момент времени, когда P 1

= u1

A cos ωt − B sin ωt = 0.

Значит, tg ωt = A/B, а так как

cos ωt = (1 + tg2ωt)1/2, sin ωt = cos ωt

·

tg ωt,

 

 

то (см. формулы (1.3) и (1.5))

max |u1| = A2 + B2, max |P1| = m1ω2 A2 + B2.

С учетом зависимостей (1.4) отсюда следует:

max |u1| =

 

 

 

 

 

,

max |P1

| = m1

ω

 

 

 

 

 

 

 

ω2

 

+ u0

 

ω2

 

+ u0.

 

 

 

v0

2

2

 

 

 

 

2

 

 

v0

2

2

 

(1.6)

(1.6a)

Глава 1

583

Формулы (1.6a) используются при построении эпюр динамических усилий. На рис. 1.2a изображена невесомая консоль, несущая на свободном

конце массу m. Так как m1 = m и

l3 u11 = (M 1, M 1) = 3EI ,

то, согласно формуле (1.2),

ω2 = 3mlEI3 .

На конструкцию действуют сила веса P = mg и сила инерции P1. Опасны два способа приложения этих сил, указанные на рис. 1.2b. Пусть, для определенности, v0 = = 0, а величина u0 задана. Тогда (см. формулы (1.6a))

max |P1| = m1ω2u0.

При ω2u0 > g огибающая эпюра изгибающих моментов имеет вид, представленный на рис. 1.2c.

Следующий пример демонстрирует ход решения задачи в случае, когда система с одной степенью свободы имеет несколько точечных масс (рис. 1.3a). Состояние системы характеризуется углом ϕ поворота жесткого диска. Последний крепится к земле при помощи шарнира A и пружины B заданной жесткости r. Жесткость r – это реактивная сила, которая возникает в пружине при изменении ее длины на единицу. На рис. 1.3b показано одно из возможных положений системы в процессе колебаний. Силы инерции P1, P2, P3 направлены в сторону перемещений u1, u2, u3 соответственно. Очевидно (u1 = 2, u2 = u3 = ),

.. ..

.. .. .. ..

Pi = −miui, u1

= 2lϕ, u2 = u3 = lϕ,

m1 = m2 = 2m, m3 = m, uB = 2lϕ.

Уравнение движения составляется как условие равновесия ΣMA = 0. При этом учитывается, что величины ui пренебрежимо малы по сравнению с размером l (рис. 1.3c):

ΣMA = 0 : P1 · 2l + P2 · l + P3 · l − VB · 2l = 0

584

 

 

 

 

 

Часть V

или

 

..

2..

 

 

2..

2

2

ϕ = 0.

8ml ϕ − 2ml ϕ − ml ϕ − 4rl

Стало быть,

..

 

4r

 

 

 

 

 

 

 

 

 

ϕ +

 

 

ϕ = 0

 

 

 

11m

 

 

и из сравнения этой записи с формулой (1.1a) следует, что

ω2 = 4r/11m.

1.2. Вынужденные колебания системы без учета сил сопротивления. Колебания, протекающие при наличии возмущающей силы P (t), называют вынужденными. Задача о вынужденных колебаниях сводится к интегрированию неоднородного дифференциального уравнения (1.1) по заданной функции P (t). Общее решение (1.3) указанного уравнения уже имеется, а частное решение можно получить, например, методом вариации постоянных. Этот метод хорош тем, что позволяет при любом воздействии P (t) записать частное решение уравнения (1.1) в квадратурах. Однако для тех возмущающих воздействий, которые рассматриваются в настоящем пособии, прибегать к готовым квадратурным формулам нет необходимости: непосредственное интегрирование уравнения (1.1) оказывается более простой операцией. Приводимые ниже решения ряда задач подтверждают сказанное.

1.2.1. Колебания системы при внезапно появляющейся силе. Пусть в момент времени t0 = 0 к системе прикладывается сила P , которая затем остается постоянной. В этом случае постоянной будет и правая часть уравнения (1.1), так что его решение очевидно:

u1 = A sin ωt + B cos ωt + u10P.

(1.7)

Следует отметить, что как раз на величину u10P переместилась бы точка сосредоточения массы, если бы сила P прикладывалась достаточно медленно, т. е. статически. Поэтому перемещение u10P называют статическим

Глава 1

585

перемещением и обозначают через uст:

 

uст = u10P.

(1.8)

Пусть при t = 0 система находится в покое. Тогда из равенств (1.7) и (1.8) следует:

1) u0

= 0 : B + uст = 0, 2) v0

= 0 : ωA = 0;

 

 

(1.7a)

u1 = uст(1 cos ωt).

Величину, заключенную в круглые скобки, называют динамическим коэффициентом и обозначают через µ:

µ = 1 cos ωt.

Этот коэффициент показывает, во сколько раз динамическое перемещение отличается от статического:

u1 = µuст.

(1.7b)

Интерес представляет максимальное значение функции µ(t). Оно достигается в тот момент времени, когда cos ωt = 1, и равно 2. Графики функций P (t) и u1(t) приведены на рис. 1.4.

1.2.2. Колебания при линейно возрастающей силе. Пусть сила P (t) линейно зависит от времени:

P (t) = βt (0 ≤ t),

где β – константа. Эта зависимость приближенно описывает нарастание давления взрывной волны. Если в некоторый фиксированный момент времени t = τ сила P (t) достигает значения P , то β = P/τ и

P (t) = Pτ t.

При такой силе уравнение (1.1) имеет решение

u1 = A sin ωt + B cos ωt +

1

uстt,

(1.9)

 

 

τ

 

в котором статическое перемещение uст определяется формулой (1.8). Если в начальный момент времени t = 0 система находилась в состоянии покоя, то перемещение u0 и скорость v0 равняются нулю. Тогда

B = 0, A = uτ стω

586

 

 

 

 

Часть V

и

uст

 

 

 

u1 =

(ωt − sin ωt).

(1.9a)

τ ω

Из этой записи видно, что

 

 

 

 

 

µ =

ωt − sin ωt

.

(1.10)

 

 

 

τ ω

 

Графики функций P (t) и u1(t) даны на рис. 1.5.

Наибольшие напряжения в конструкции возникают в тот момент времени, когда динамический коэффициент (1.10) максимален. Если нагрузка действует в интервале времени [0, t ], то надо найти такое время t1 из этого интервала, при котором функция (1.10) достигает глобального максимума.

Так как

.

 

1

cos ωt

..

ω sin ωt

 

 

 

µ(t) =

 

τ

, µ(t) =

 

,

 

 

 

 

 

 

 

 

 

 

τ

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

.

..

 

 

 

 

 

 

µ(t) и

µ(t) одновременно обращаются в нуль при ωt = 0,

π, 2π, . . . Сказанное означает, что все стационарные точки – это точки перегиба, а потому функция µ(t) монотонно возрастает по t. Стало быть, наиболее опасен момент времени t1 = t .

1.2.3. Колебания при гармонической возмущающей силе. Гармоническая возмущающая сила возникает, в частности, при вращении неуравновешенных частей электродвигателей. Так, если ось вращения ротора массой m0 находится от оси, проходящей через центр массы m0, на расстоянии e, то возникнет центробежная сила (рис. 1.6)

Q = m0v2 . e

Линейная скорость v движения центра масс может быть выражена через угловую скорость θ вращения ротора: v = θe. В свою очередь (см. формулу (3) на с. 580), θ = πn0/30, так что

π2n2m e

Q = 0 0 . (1.11) 900

В ряде случаев расчетные усилия в конструкции обуславливаются действием лишь вертикальной составляющей

QV = Q sin θt

Глава 1

 

 

 

 

 

 

587

силы Q. При P (t) = QV уравнение (1.1) принимает вид

 

..

+ ω

2

u1

= ω

2

u10Q sin θt,

(1.12)

u1

 

 

а его решение записывается следующим образом:

 

u1 = A sin ωt + B cos ωt + C sin θt.

(1.13)

Константы A и B связаны с общим решением. Они определяются начальными условиями задачи, тогда как постоянная C, характеризующая частное решение

u1 = C sin θt,

должна удовлетворять исходному уравнению (1.12):

−θ2C sin θt + ω2C sin θt = ω2u10Q sin θt;

C =

 

ω2

 

u10Q.

 

ω2

θ2

 

 

 

 

 

 

Пусть

 

 

 

 

1

 

uст = u10Q,

µ =

 

.

 

 

1

(θ/ω)2

 

 

 

 

 

 

Тогда

C = µuст.

(1.14)

(1.15)

(1.16)

Если двигатель включается в момент времени t0 = 0, то, согласно формуле (1.14),

u1(0) = 0, u˙ 1(0) = Cθ.

Эти значения принимаются в качестве начальных перемещения u0 и скорости v0 для определения констант A и B (см. равенство (1.13)):

1) u0 = 0 : B = 0; 2) v0 = : = + Cθ A = 0.

Следовательно,

u1 = C sin θt,

(1.17)

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

Величина C представляет собой наибольшее динамическое перемещение точечной массы, которое, согласно зависимостям (1.15) и (1.16), равно статическому перемещению, умноженному на динамический коэффициент. Из записи (1.17) следует, что угловая скорость θ есть не что иное, как круго-

588

Часть V

вая частота вынужденных колебаний, обусловленных гармонической силой. Характер изменения динамического коэффициента при наращивании оборотов двигателя представлен на рис. 1.7. Если θ = n0 = 0, то µ = 1: при неработающем двигателе нет и динамического эффекта. По мере увеличения частоты θ в диапазоне от 0 до ω коэффициент µ возрастает. При

θ = ω наступает резонанс. Теоретически перемещение u1 становится неограниченным: µ → ∞. На самом деле величина u1 конечна (см. далее п. 1.5), но она столь велика, что появляется угроза выхода конструкции из строя. Вот почему правилами проектирования строительных конструкций верхнее значение коэффициента µ ограничено числом 2: µ ≤2. Из формулы (1.15)2 следует, что это ограничение можно записать так:

θ ≤ 0, 7ω.

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

1.2.4. Колебания при кинематическом возмущении системы. Примером динамического кинематического воздействия служат колебания основания сооружения при землетрясении. Пусть U0(t) – заданное смещение некоторой опоры системы с одной степенью свободы, а u10 – перемещение по направлению движения массы m1 от воздействия U0 = 1. Поскольку масса m1 совершает также и собственные колебания u1(t), то колебательный процесс поддерживается силой

P

(t) =

m

u¨

1

m

u

U¨

0

,

(1.18)

1

 

1

 

1

10

 

 

 

а так как u1(t) = u11P1(t), то (см. обозначение (1.2))

u¨1

+ ω2u1 = −u10U¨0.

(1.19)

¨

для силы инерции, обусловленной кинемати-

При обозначении P0 = −m1U0

ческим воздействием, уравнению (1.19) можно придать вид, мало чем отличающийся от традиционной записи (1.1):

 

u

 

u¨1 + ω2u1 =

10

P0(t).

(1.19a)

 

 

m1

 

Глава 1

589

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

о зависимости ускорений ¨ оснований сооружений от времени. Графики за-

U

висимостей ¨ ( ), называемые акселерограммами, имеют вид, представлен-

U t

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

то получится некая усредненная диаграмма ¨ ( ), ко-

U t

торую можно аппроксимировать какой-либо гладкой кривой (рис. 1.8b). Зачастую в качестве аппроксимирующей кривой берется периодическая функция

U¨ (t) = ate−bt sin θt,

(1.20)

где a, b, θ – константы, определяемые по усредненной акселерограмме.

Получить решение уравнения (1.19) с правой частью, порождаемой функцией (1.20), в квадратурах не удается. Поэтому правую часть зависимости (1.20) заменяют отрезком ряда:

U¨

n

(t) = Ak sin kθt

 

k

 

=1

и сводят решение исходной задачи к решению серии задач с гармонической возмущающей силой (при k = 1, 2, . . . , n). После объединения всех n решений вычисляется сила инерции (1.18) и строятся динамические эпюры усилий. Ясно, что эта задача довольно громоздка даже в том нетипичном для строительной практики случае, когда система имеет только одну степень свободы.

1.2.5. Колебания при импульсивном воздействии. Импульсивным называют воздействие, которое прикладывается к системе за очень малый промежуток времени [0, τ ], а затем исчезает. Характерной особенностью указанного воздействия является то, что закон изменения силы P (t) на интервале [0, τ ], а зачастую и сама длительность τ приложения этой силы остаются неизвестными. Интегральная же характеристика

τ

 

 

S =

P (t)dt

(1.21)

0

590

Часть V

кратковременной нагрузки, именуемая импульсом силы, может быть установлена по количеству движения тела, оказывающего воздействие на систему. Однако прежде, чем воспользоваться импульсом S для исследования ее (системы) колебаний, полезно на примерах убедиться в том, что, не зная его, решить задачу о движении системы при действии произвольной кратковременной силы не удастся. Для этого достаточно задаться двумя различными законами изменения силы P (t) на одном промежутке времени, получить основные характеристики движения системы для каждой из нагрузок и сравнить результаты.

Итак, пусть сначала сила P на интервале [0, τ ] не меняется (рис. 1.9a). Решение уравнения (1.1) можно записать, отталкиваясь от формул (1.7a) и (1.3):

u1(t) =

uст(1 cos ωt), 0 ≤ t ≤ τ ;

(1.22)

 

 

A sin ωt + B cos ωt, τ < t.

 

 

 

 

 

Таким образом, при t [0, τ ] масса m1 движется по закону, отвечающему внезапно приложенной возмущающей силе, а после снятия нагрузки система совершает свободные гармонические колебания. При t = τ из верхней строки (1.22) следует

u1(τ ) = uст(1cos ωτ ), u˙ 1(τ ) = ωuст sin ωτ.

Величины u1(τ ) и u˙ 1(τ ) принимаются в качестве начальных параметров u0 и v0 для определения констант A и B, характеризующих движение массы

на этапе τ < t:

 

A sin ωτ + B cos ωτ = uст(1 cos ωτ ),

A cos ωτ − B sin ωτ = uст sin ωτ.

 

 

Отсюда можно установить, что

 

A = uст sin ωτ, B = uст(cos ωτ − 1).

(1.23)

График функции (1.22) изображен на рис. 1.9b. Штриховой линией показан закон колебаний при внезапно приложенной силе P , которая в интервале [0, ∞] не меняется (ср. с рис. 1.4).

Максимальное значение перемещения u1 можно найти при помощи первого из соотношений (1.6), если подставить в него константы A и B по формулам (1.23):

 

 

 

ωτ

max u1 = uст 2(1 cos ωτ ) = 2uст sin

 

.

2