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

Прикладная теория колебаний

..pdf
Скачиваний:
23
Добавлен:
15.11.2022
Размер:
3.77 Mб
Скачать

Поскольку скорость в момент сжатия равна 0, K2 = 0, и тогда

A = −K = −

Gv2

.

(4.37)

1

2g

 

 

 

 

 

Сила, действующая на бульдозер,

 

 

Pд = cλд,

 

(4.38)

где c – жесткость конструкции и преграды; λд – величина суммарного сжатия ножа и тупикового ограждения.

Учитывая, что сила, действующая на бульдозер, будет совершать отрицательную работу, получим:

 

P

λ

д

 

P2

 

 

A = −

д

 

= −

д

,

(4.39)

2

 

2с

 

 

 

 

 

или, приравнивая (4.37) и (4.39), получим:

Gv2 = − Pд2 . 2g 2с

Сила после наезда

P = v Gc .

(4.40)

д g

Например, при весе машины G = 60 т и v = 0,1 м/с сжатие от Р = 1 т равно 1 см, следовательно:

 

с = 1 т/см = 100 т/м,

P = v Gc

= 0,1

60 100 2,5 т.

д

g

 

9,81

 

 

91

Стр. 91

ЭБ ПНИПУ (elib.pstu.ru)

4.5. Определение допустимой скорости транспортировки

В практике эксплуатации строительных и дорожных машин часто встречаются задачи по определению допустимых с точки зрения динамических нагрузок на конструкцию параметров движения транспортных средств. Рассмотрим определение допустимой скорости движения прицепа и собственных частот колебаний конструкции [9] .

Для прицепа массой m, двигающегося со скоростью v (точка 0 не имеет вертикальных перемещений) (рис. 4.10), определим скорость установившегося движения системы, когда амплитуда колебаний достигает максимальной величины.

Рис. 4.10. Расчетная модель прицепа

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

h= h0 1cos 2lπx ,

1

где x = vt; h – высота неровности; h0 – амплитуда неровностей; l – длина волны.

Момент инерции массы прицепа относительно точки 0 равен J0. Трение вязкое, коэффициент трения α. Масса колес при-

92

Стр. 92

ЭБ ПНИПУ (elib.pstu.ru)

цепа по сравнению с массой прицепа m мала, ею можно пренебречь.

Воспользуемся уравнением Лагранжа второго рода. Определим кинетическую и потенциальную энергию и диссипативную функцию Рэлея:

 

 

 

 

K =

mv2

+

J

ϕ2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

;

 

 

 

 

 

 

 

 

(4.41)

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

W =

 

c(lϕ−h)

2

;

 

 

 

 

 

 

 

 

 

(4.42)

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

R =

 

α(lϕ−h)2

.

 

 

 

 

 

 

 

 

 

(4.43)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставив выражения (4.41)–(4.43) в уравнение Лагранжа:

 

d

 

K

 

 

K

 

 

W

 

 

R

 

 

 

 

 

 

 

 

∂ϕ

∂ϕ

= −

∂ϕ

∂ϕ

,

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J0ϕ+αl(lϕ−h) +cl(lϕ−h) = 0

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ+ 2ξϕ+

 

 

 

 

 

 

 

 

 

2πvt

 

 

 

 

 

p2h

 

p02ϕ = bsin

 

 

 

 

 

 

−γ

+

 

0

0

,

(4.44)

 

l1

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b =

2πvh0αl

2 + p4

 

h0

 

2

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l1J0

 

 

 

 

 

0

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

tg γ =

p02l1

; 2ξ = αl2 ; p

2 = cl2 .

 

 

 

 

 

 

4ξπv

 

 

 

 

 

 

J0

 

 

 

 

 

 

0

 

 

 

J0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение уравнения (4.44) можно представить в виде

 

 

2πvt

 

 

 

 

 

 

 

 

 

 

2πvt

 

 

 

+ D.

(4.45)

ϕ = Asin

l1

 

−γ

+ B cos

l1

−γ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

93

Стр. 93

ЭБ ПНИПУ (elib.pstu.ru)

Очевидно, что D = hl0 , а для определения постоянных A

и B получим систему уравнений:

 

2

 

2πv 2

 

 

 

 

 

2πv

 

B

 

 

 

 

 

p

 

 

 

 

 

 

 

A 2ξ

 

 

 

 

 

 

= b,

 

 

 

0

 

 

l

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

2ξ

 

2πv

+

2

 

2πv

 

 

=

 

 

 

 

 

 

 

A

p0

 

l1

 

 

B

0.

 

 

 

 

 

l1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2πv

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b p0

l1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

(4.46)

 

 

 

 

 

 

 

2 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

2πv

+ 4ξ

2

 

2πv 2

 

 

 

 

p0

l1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2πv

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2ξb

l1

 

 

 

 

 

 

 

 

 

 

 

 

 

B = −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(4.47)

 

 

2

2πv 2 2

 

 

 

 

2

2πv

2

 

 

 

p0

 

l1

 

 

+4ξ

 

 

 

l1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение уравнения движения (4.44) можно также пред-

ставить в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2πvt

 

 

 

 

 

h

 

 

 

 

 

 

 

 

 

 

ϕ = ϕ0 sin

 

 

+

0

 

,

 

 

 

 

 

(4.48)

 

l1

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ϕ0 = A2 + B2 =

 

 

 

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

2πv

4ξ

2

2πv 2

 

 

 

 

 

 

 

 

p0

l1

 

 

 

+

 

 

l1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

94

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Стр. 94

ЭБ ПНИПУ (elib.pstu.ru)

5. СТАТИСТИЧЕСКИЕ МЕТОДЫ ИССЛЕДОВАНИЯ ДЕЙСТВИЯ ДИНАМИЧЕСКИХ НАГРУЗОК

НА КОНСТРУКЦИЮ СДМ

Динамические нагрузки, действующие на конструкцию СДМ, зависят от большого количества факторов, сочетание которых носит случайный характер. Например, в турбулентной атмосфере точная форма порыва ветра всегда является в какойто мере неопределенной [2, 5, 9].

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

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

Изучением закономерностей случайных процессов занимается математическая статистика. Рассмотрим основные понятия метода математической статистики, которые используют при анализе динамических процессов.

5.1. Понятие вероятности. Функции распределения вероятностей

Рассмотрим простой опыт с бросанием монеты. Пусть благоприятным событием A будет выпадение герба. В первых n экспериментах событие A произошло ν раз.

Число ν, показывающее, сколько раз произошло данное событие A, называется частотой. Частотностью события A

называется отношение ν/n. Очевидно, что 0 ≤ nv ≤ 1.

97

Стр. 97

ЭБ ПНИПУ (elib.pstu.ru)

Эксперименты показывают, что при возрастании числа опытов частотность ν/n определенного события A стремится к некоторому более или менее постоянному значению, что отражено на рис. 5.1.

Рис. 5.1. Зависимость частотности опытов от числа опытов

Устойчивость частотности для большого объема испытаний, повторяемых при одинаковых условиях, известна давно. Имеются основания допустить, что для любого события A, связанного со случайным экспериментом, можно определить такое число P, при котором частотность появления события A при большом повторении опыта приблизительно будет равна P. Это число P называется вероятностью события A. Ясно, что

0 ≤ P ≤ 1.

Если событие A является невозможным, то Р = 0. Если же A достоверное событие, то Р = 1.

В качестве примера случайной величины рассмотрим усилие N, измеряемое в элементе конструкции.

Пусть при первом измерении получено значение усилия N1, при втором N2 и т.д., при n-м измерении Nn. Будем откладывать на графике вдоль оси абсцисс величину N, а вдоль оси ординат вероятность того, что усилие не превышает значения N = x. Эту вероятность приближенно можно определить при помощи суммирования частотностей всех измеренных значений усилий, меньших или равных х. Обозначим ее Р(N x). Полученные значения вносим в таблицу и строим по ним график.

98

Стр. 98

ЭБ ПНИПУ (elib.pstu.ru)

Рис. 5.2. Кривая интегрального распределения вероятностей

Кривая F(x) = Р(N x) называется кривой интегрального распределения вероятностей

(рис. 5.2).

Интегральная функция распределения F(x) есть неубывающая функция x, и для нее справедливы соотношения:

0 F(x) 1, F(–) = 0, F(+) = 1.

Производная от функции распределения

dF(x) = ω

dx

(x), (5.1)

если она существует, называется плотностью вероятности или функцией распределения случайной величины.

Распределение бывает двух типов – дискретное и непрерывное.

Случайная величина x называется величиной дискретного типа, если она может принимать конечное или бесконечное счетное множество значений, т.е. такое множество значений, элементы которого могут быть занумерованы в каком-нибудь порядке и выписаны в последовательность х1, х2, …, х.

Распределение дискретной случайной величины х будет полностью определено, если указать для любого значения v(v 1, 2, …) вероятность Рv того, что х принимает значение х0. Так как вероятность того, что события х1, х2, …, х произойдут, равна 1, то

n

Pv =1.

v=1

При этом интегральная функция распределения F(x) задается соотношением

99

Стр. 99

ЭБ ПНИПУ (elib.pstu.ru)

F(x) = P(xv x) = Pv ,

(5.2)

xv x

 

где суммирование распространено на все значения индекса v,

для которых x0 x. Таким

образом, P(x) – ступенчатая

 

функция, равная последова-

 

тельной величине на лю-

 

бом интервале без точки x0

 

и имеющая скачок Pv в каж-

 

дой точке xv . Распределение

Рис. 5.3. Дискретное распреде-

дискретного типа графически

можно представить ступенча-

ление случайной величины

той кривой (рис. 5.3).

 

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

Для непрерывной величины х функция распределения Р(х)

непрерывна, плотность вероятности

ω(x) = F(x)

существует

для всех значений х.

 

 

В этом случае

 

 

F(x) = P(x0 < x) = x

ω(x)dx.

(5.3)

−∞

 

 

Распределение непрерывной случайной переменной величины можно представить графиком функции распределения или графиком плотности вероятности данного распределения ω(x)

Рис. 5.4. Плотность вероятности (рис. 5.4).

100

Стр. 100

ЭБ ПНИПУ (elib.pstu.ru)