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

книги из ГПНТБ / Жуков А.В. Колебания лесотранспортных машин

.pdf
Скачиваний:
11
Добавлен:
27.10.2023
Размер:
20.06 Mб
Скачать

При обработке методом пересечений

м а т е м а т и ч е с к о е

о ж и д а н и е

высоты

неровностей

определяется по формуле

M[H]=I,Jll

P(HT

) Н T,

где />(//,)

— частота

появления

не­

ровностей высотой HT.

 

 

 

 

 

F(Hh

 

 

 

 

 

F(H)

 

W

 

 

 

 

 

 

 

 

 

 

 

 

 

0,8-

 

 

 

 

 

 

 

 

о.г-

 

 

 

 

 

 

0.5 Л

 

 

 

 

 

 

ОАЛ

 

 

 

 

 

 

 

 

0.1-

 

 

 

 

 

 

о.гл

 

 

 

 

 

 

 

 

о

 

10

го

н,см

 

Рис. 23.

Распределение высот

дорожных

неровностей

 

 

 

опытного участка дороги.

 

 

Д и с п е р с и я

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

D [ H ] характеризует

ее

рассеивание относительно центра группирования. Она представ­ ляет собой математическое ожидание соответствующей центри­

рованной

величины, т. е. D [ H ]

= 2 " = і

Р{Ні ) • {HT — MHT ) 2 ,

где

(Н;—МНі)

— центрированная

высота

неровности. Кроме

того,

рассеивание случайной величины характеризуется средним квад­

ратичным

значением,

в нашем случае — с р е д н е й

к в а д р а ­

т и ч н о й

в ы с о т о й

н е р о в н о с т и о„ , которая равна

У D [ H ] .

Величина з н

удобна при сравнении различных участ­

ков дорог, ее размерность совпадает с размерностью

случайной

величины.

Исследования показывают, что .кривые распределения высот дорожных неровностей, построенные по экспериментальным дан­ ным (рис. 23), близки к кривым нормального закона распреде­ ления.

В е р о я т н о с т ь

случайного

 

совокупного

 

расхождения

Р( х 2 )

между наблюденными и

выравнивающими

частотами

нормального

распределения близка

к единице

 

(более 0,95).

Поэтому выравнивание с помощью нормального

распределения

можно считать

хорошо

согласующимся с данными

наблюдений.

Для

стационарного

случайного

процесса

 

корреляционная

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

t

первого

аргумента

на оси

абсцисс, а только от промежутка

t

между

первым и

вторым

аргументами, т. е. Rx{t,

t-\-t)—R ( т ) .

 

стационарного

Следовательно, корреляционная

функция

процесса есть функция не двух, а одного аргумента. Это обстоя­

тельство сильно облегчает операции

над стационарными

случай­

ными функциями.

 

 

 

 

 

 

 

Основные свойства

корреляционной функции

следующие:

1)

симметричность

(R

( x)=R(—t

) ) ; 2) при

х = 0

R(0) =

= £ > [ / / ! ( 0 ] > 0 ; 3) / ? ( 0 ) > | / ? ( т )

| ;

4) на величину / ? ( т ) три

х=0

влияет высота

выступов

и глубина исследуемой

кривой;

5) с увеличением длины неровностей исследуемой кривой точка пересечения корреляционной функции с осью абсцисс удаляется

от начала

координат

(т. е. увеличивается в р е м я к о р р е л я ­

ц и о н н о й

с в я з и

to); 6) если в исследуемую кривую входит

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

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

р ( т ) = Я ( т ) / # 0 .

Чтобы определить нормированную корреляционную функцию,

если скорость движения

отлична от

1 м/с, достаточно

разделить

значение аргумента t

для каждого

значения р ( t )

на величи­

ну этой скорости, выраженной в метрах в секунду, оставив зна­ чения р ( х ) без изменения (рис. 24).

В процессе движения транспортная машина кроме продоль­ ных совершает и поперечно-угловые колебания, которые возни­

кают из-за того, что правые и левые колеса

машины

неодновре­

менно наезжают

на выступы

или впадины,

т. е. микропрофили

левой и правой

колеи дороги

не совпадают.

Отсюда

ясно, что

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

воздействие, как и в продольной плоскости,

носит случайный

характер.

 

 

Для исследования влияния дороги

на

поперечно-углшые

колебания машины строят г р а ф и к и

п р е в ы ш е н и й одной

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

г>(т)=АіЄ-аМ+А2е-"М

c o s p - : ;

р( х) = е - а 1 т 1 COS Р т .

Р(П

1.0-

Рис. 24. Нормированные корреляционные функции участ­ ков дорог:

 

1, 3,

5 расчетные

кривые;

2, 4, 6 — экспериментальные

 

 

 

 

 

 

 

кривые.

 

 

 

 

 

Входящий

в уравнения

коэффициент

а характеризует

б ы-

с т р о т у

з а т у х а н и я

корреляционной функции,

а

^ — ее

к о л е б а т е л ь н ы е

с в о й с т в а

(см. рис. 24).

 

 

 

При малых значениях

 

а мииропрофиль дороги

по

харак­

теру ближе к периодическим

колебаниям частоты, равной

 

р ,со

случайной

амплитудой и фазой.

При больших его

значениях

преобладания тех или иных частот не наблюдается.

 

 

 

Корреляционные

функции являются

в р е м е н н ы м и

ста­

тистическими характеристиками. Для статистического исследо­

вания

динамической

системы нужны не временные, а ч а с т о т ­

н ы е

характеристики

воздействия. Обычно пользуются энерге-

тическими спектрами (спектральная плотность воздействия), которые являются изображением Фурье корреляционной функ­ ции.

Э н е р г е т и ч е с к и й с п е к т р есть спектр корре­ ляционной функции Ф ( ы ) случайного стационарного процесса,

который

пропорционален

квадрату

амплитуд.

Для получения

такого

опектра достаточно

взять

от нее интеграл

Фурье

[52, 53,

54], т.

е.

 

 

 

 

 

 

 

 

 

 

Ф ( ш ) = ~ | _ ~ Я Ы е - ' ^ с » .

 

(15)

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

зоваться

выражением

 

 

 

 

 

 

 

 

 

R (т) =

| _ * Ф

(ш)

е-'^сіш.

 

 

Проинтегрировав

выражение

(15) подстановкой в него зна­

чения

корреляционной

функции,

получим

расчетную

формулу

для вычисления значений спектральной плотности.

 

Если

корреляционная

функция

задана

формулами:

 

 

 

 

Р , ( х ) = Л е—М;

 

 

 

 

 

р2 ( т ) =А e-*W cos$x ;

 

 

 

 

 

Р з ( х) ге

-*'W+A2e-a'\'-\.tos$i,

 

 

выражения спектральных

плотностей в преобразованном виде

[53]будут соответственно:

ф, ( ш ) = е д Л а / [ Ц а 2 + о ) 2 ) ] ;

Ф - н ^ [ , +

, ( : + Р

)

2 + а

г +

( : _ Р )

г

] ;

о б )

Ф

1 ;

ті( а 2 ! + ш2)

^

2тг

[ а 2

2 + ( о

+

{))2

"t"

Величина нормированной спектральной плотности опреде­

ляется по формуле

р ( ш ) = Ф

)/#(0). Для вычисления спект­

ральных

плотностей при скоростях,

отличных от скорости 1 м/с,

необходимо

в их

выражения

подставить

значения

к о э ф ф и ­

ц и е н т о в к о р р е л я ц и о н н о й

с в я з и

а и

р

для нужной

скорости

движения.

 

 

 

 

 

 

 

 

Графики Ф (со) являются дорожными спектрами и показы­ вают, на какой частоте и с какой интенсивностью дорога воз­ буждает колебания движущейся транспортной машины (рис. 25). Как видно из графиков, дорога возбуждает колебания не одной

какой-либо частоты, а целый спектр колебаний. Ширина диапа­ зона частот зависит от скорости движения и качества дороги.

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

Рис. 25.

Энергетический спектр воздействия участка

 

дороги:

/ — г> = 10

км/ч; 2 — и = 15 км/ч; 3 — t) = 30 км/ч; 4 с = 60км/ч .

из анализа приведенньюс уравнений, наибольшее значение спект­ ральной плотности будет при ю = Р . Это значит, что наиболь­ шая величина 'возмущения наблюдается при этой частоте. Фи­ зически это объясняется тем, что коэффициент корреляционной связи р одновременно является и частотой главной гармоники, которая входит и в выражение корреляционной функции и в состав случайной функции, описывающей дорожный микропро­ филь. Поэтому при этой частоте наблюдается всплеск спект­ ральной плотности, причем количество таких всплесков будет равно количеству членов выражения спектральной плотности, содержащих коэффициенты р.

Из выражений (16) и (17) видно, что гари « = 0 значения спектральной плотности не равны нулю, т. е. они не проходят

через

начало координат. Из уравнения

(16), если

ш = 0 ,

имеем

из

уравнения

Ф 2 ( 0 ) = = : Я ( 0 ) Л а / [ * ( а Ч - Р 2 )],

 

 

 

(17) —

 

 

 

 

 

 

 

Фэ (0) =

Д (0)Ліа,/(я а 2 ,) +R(0)A2

«2 /[ г. ( а2

2 4-р2 ) ] .

 

 

Величина

максимума

зависит от

соотношения

параметров

а

и

р , от соотношения

постоянных коэффициентов

А\

и А2 и

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

При уменьшении а убывание кривых корреляционных фун­ кций будет более медленным, соответственно в спектре случай­ ной функции большой удельный вес приобретают малые часто­ ты. В этом случае кривая спектральной плотности вытяпивается вверх, одновременно сжимаясь с боков. Если а увеличивается, преобладание малых частот становится все менее выраженным.

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

спектральной

плотности.

 

 

 

Установлено,

что высота

неровностей

в одинаковой степени

увеличивает

или

уменьшает

ординаты

спектра,

влияние же

длины неровностей более сложно.

 

 

Энергетические спектры воздействия зависят только от ка­

чества дороги и скорости движения. Ими можно

пользоваться

три проектировании любой транспортной

машины.

 

Глава Ш

МЕТОДЫ ИССЛЕДОВАНИЯ

1

Свободные колебания

При отсутствии возмущающей силы,, если состояние равно­ весия упругой системы было нарушено и затем она предостав­ лена самой себе, система совершает с в о б о д н ы е к о л е б а - н и я.

Рассмотрим простое стандартное уравнение [55], которое применимо ко многим механическим линейным системам с од­ ной степенью свободы. Уравнение имеет вид

 

 

х-{-р2х=0,

 

 

(18)

где

х — перемещение;

 

 

 

 

P—Vcjm

— частота свободных колебаний

жесткость

 

 

упругого элемента, т — масса).

 

 

 

Ч а с т о т а

колебаний р зависит от параметров

системы, но

не зависит от начальных условий. Она является

важным показа­

телем, характеризующим колебательные свойства системы.

 

Уравнение

(18) имеет решение ;в виде х—Сх

sin

pt-\-C2smpt,

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

характеризующиеся

н а ч а л ь н ы м с м е щ е н и е м

Хо

и

на ­

ч а л ь н о й с к о р о с т ь ю v0,

т. е. х=х0,

v = v0

при

ґ = 0 .

Так

как Cj Vo/p; С20,

решение

уравнения

(18)

можно

перепи­

сать в виде

x=xocosp^+ (v0sinpt) I p.

Это уравнение соответствует простой системе с одной степенью свободы. С увеличением числа степеней свободы решение задачи усложняется.

Рассмотрим свободные колебания механической системы с двумя степенями свободы. Возьмем для примера транспортную систему — трелевочный лесной трактор, перемещающий пакет деревьев, подвешенный за один конец. Упрощенная расчетная схема поперечно-угловых колебаний приведена на рис. 26, а.

Данная динамическая система имеет две степени свободы,

которые характеризуются

обобщенными

координатами

<р и р.

Потенциальная энергия

выражается как

сумма П =

П 1 +іП 2 ,

причепи

 

 

 

П ^ С р Д / 2 ! ^ 0,5 GtfK <f2 ;

U2=Gn h = 0,5 Ga L ? 2 — 0,5DG„ p2,

где h — вертикальное перемещение точки 0 3 при отклонении •подрессоренной массы на угол ср и подвеса пачки на угол Р;

ср — жесткость рессоры;

п— число осей;

h

— расстояние между рессорами;

G

— вес подрессоренной массы трактора;

Нк

— расстояние от оси крена до центра тяжести подрес­

 

соренной .массы трактора;

G n

вес части пакета, приходящийся на трактор;

L

— расстояние от центра крена до точки подвеса 02.

Рис. 26. Схема поперечно-угловых колебаний трактора при полуподве­ шенном способе трелевки леса (а) и схема отклонений подрессоренных

масс (б).

Общее выражение для потенциальной энергии будет

П = 0 , 5 (0,5 n/Vp — GHK — GaL)

ср2 +0,5 Ga

D p2.

 

Кинетическая энергия равна сумме кинетических

энергий

первого Ті и второго Т2

тел:

 

 

 

 

 

Г = 0 , 5 h

?2 +0,5 / 0 2 п р ' 2 + 0 , 5

G'„ v\

/g,

 

 

где / 0 2 п — момент инерции

пачки относительно

оси,

проходя­

щей через

точку

подвеса и центр

касания

второго

конца иачки с землей;

 

Іт

— момент инерции относительно оси крена;

 

vy

— скорость перемещения центра тяжести пачки;

v2y^v20/9

(v0

— скорость точки

0 3

пакета);

 

 

G'п

— полный вес пачки.

обоб­

Подставив

значение скорости v2^, выраженной через

щенные

координаты (ірис. 26,6),

и

выполнив некоторые

преоб­

разования, получим

 

 

 

T = 0,5[/r +

( D + L ) 2 G / n / ( 9 g ) ] ? 2

+ 0

, 5 ( / ° 2 П + G'n D2/(9g))

р 2 +

+[DG'n(l+L)/(9g)]V<?.

Применив уравнение Лагранжа, (будем иметь систему урав­

нений:

 

 

 

 

 

 

 

[/т+ G'n(D+L)2 /(9g)

 

] ip+ [DG'n

(1 +1)/ (9g) ] p +

(0,5С р м/2 , -

 

— GHK — GUL) ? =

0;

 

[DG'n(l+L)/(9g)]

?'+,[/°2п

+

GnD2/

( 9 g ) ] p + G „

D?=0,

или

 

 

 

 

 

 

 

 

a u cp+a 1 2 p - f cn'f = 0 ;

 

 

ai2 'f'+a2 2 P+C22p =

0,

(19)

где a „ = / T + G / n ( ^ + ^ ) 2

/

(9g);

 

 

cn=Q£cvnl2x—GH-GnL;

al2=DG'n(\+L)/(9g);

 

 

cu=0;

 

(20)

a2 2

= / ° 2 n + G ' n D 2 / ( 9 g ) ;

 

c22=GnD.

 

Уравнение чіастот будет иметь вид

 

 

{(0,bcpnPl-GHK-GnL)-k*[l1+G'a(D+L)>/(9g)])

 

 

[ G n D —

-

£ 2 ( / ° 2 n + G ' n D 2 / ( 9 g ) ) ] -

[ D G ' n ( l + L ) 2 / ( 9 g ) ] 2 - ^ = 0 .

Квадраты п а р ц и а л ь н ы х

ч а с т о т равны:

 

n2 i =

cn/a u ;

 

n22 =

c22/a22

 

(пі — частота свободных колебаний системы, если представить, что масса хлыстов сосредоточена в точке 0 2 ; п2 — частота сво­ бодных колебаний физического маятника, образуемого пачкой, когда точка 02 неподвижна).

Введя обозначение

X2=[DG'n(l+L)4(9g)]2/{[Ir+G'n(D+L)2li9g)]

• (/02п+

+ G'n D2 /,(9g))),

получим

(1 + Х 2 ) ^ 4 - ( n 2 i + n 2 2 ) f c 2 + n V * 2 2 = 0 ,

(21)

откуда

находим квадраты ч а с т о т

г л а в н ы х

 

к о л е б а н и й

системы:

 

 

 

 

 

 

 

 

 

 

 

 

 

k2i = 0&[n*l+n*2

K ( n V H 2 i ) 2 - H

X2 "2 i"2 2]/(1 — х 2

) ;

 

^ 2 2 = 0 , 5 [ п 2 1 + п 2 2

V(n22-n2,)2+4

 

х2 «2 іл2

2 ]/(1 — X2 )-

(22)

К о э ф ф и ц и е н т ы

ф о і р м

главных

-колебаний

равны:

 

 

C i 2

 

k2i

_ ft _ а 1 2

 

^ 2 2

 

 

 

, 9 о ч

 

^ =

~

*

^ 2

р -

'

Р 2 —

~ 2

р

- -

 

1<")

 

 

#11

 

М і - К 2

 

 

 

 

Я 1 —

й

2

 

 

Уравнение

(21) должно иметь два положительных

корня k2i

и k22,

причем

должны

соблюдаться

условия: 0<.k2i^

 

сии;

С22/Я.22 <

^ 2 2 < +

° ° -

 

 

 

 

 

 

 

 

 

 

 

 

іПірсоставлении

системы

(19)

необходимо

 

соблюдать

условия:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

О ц > 0 ; а 2 2 > 0 ;

 

а п а 2 2 — а 2 і 2 > 0 ;

 

 

 

 

 

с ц > 0 ;

с 2

2 > 0 ;

 

Сцс2 2

с 2

1 2 > 0 . '

 

 

 

Коэффициенты форм главных колебаний связывают между

собой обобщенные координаты: cpU) =

fjjfi(0

;

<р<2) =

Р 2 { П 2 ) .

 

Из

формул

(23)

видно,

что Рі >0; Рг <0,

поэтому при

главном колебании низшей частоты k\

знаки

<р и $ одинаковы,

а при колебаниях

высшей частоты — различны.

 

 

 

 

Примем, что координаты

tp и Р совершают

гармонические

колебания одинаковой частоты k, одинаковой или прямо проти­

воположной фазы, отличающиеся друг от друга

только ампли­

тудами. Тогда

с? =Axsm(kt-\-b\)\$

=A2sm(kt-\-o.2).

 

Общее

решение

системы (19) представится в виде:

<? =

Сі (с2 2 £2 ia2 2 )

s i n (kxt +

— C 2 (cx2

— k22aX2)

sin (k2t-\- a2 );

 

 

 

 

 

 

 

(24)

P =

— Ci(c2i

k2xan)

s i n (kit-}-0Li)+C2{cn

k22an)

s i n ( k 2 t + a 2 ) .

 

Уравнения

(24) содержат постоянные C\, C2,

ax и a2 , кото­

рые должны быть определены из начальных условий. При г = 0 задаются значения обобщенных координат и скоростей: ср = <р0;

В результате решения системы, состоящей из четырех урав­

нений с четырьмя неизвестными

(Си С2 , а! и а 2 ) , получим:

. e i

=

atctg

 

ti^+fo^b;

 

 

 

 

CCS4 +Po S2

 

=

arCtg

(УО S3+fto Si) &2.

a2

 

: —:

,

 

 

 

 

CPS3+P0S4

 

С

?os

* + Po S 2

 

 

1

(S1S41— S3S2) - Sin a,

'

Соседние файлы в папке книги из ГПНТБ