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

книги из ГПНТБ / Барский И.Б. Динамика трактора

.pdf
Скачиваний:
21
Добавлен:
25.10.2023
Размер:
14.66 Mб
Скачать

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

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

корреляционные функции. Для

воздействий

qx(l*) и qy{l*)

вза­

имные корреляционные функции имеют вид

 

 

Rxy(l*)

= \\m

-±-

\ qx(l)qy(l

+

l*)dl;

 

 

 

- L ,

 

}

(53)

 

 

 

 

 

Ryx(l*)

= \im

- J -

Г qy(l)qx(l

+

l*)dl.

 

 

 

 

 

 

)

 

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

VRA0)Ry(0) '

Для определения авто- и взаимно корреляционной функции по известной реализации случайного процесса интегралы (52) и

(53) заменяются конечными

суммами

 

 

 

JV-ц

 

Rx(\i)zzz

7

ЗД+ц!

 

N-\x

j£j

 

N-ti

W-H

[*

£

где u = — ; p. = 0, 1, 2, ...; N — число интервалов; v

= — ; v = 1,

Д

Д

2, ... ; здесь Д длина разбиения пройденного пути

(шаг).

130

Интервал L 0 = N& определяется из соотношения L 0 = Ю / т а х ,

где / т а х — максимальная волна функции неровности, т. е. мак­ симальный отрезок оси / между двумя соседними нулями функ­ ции.

Шаг Д может быть либо определен по формуле Д =

где

/ mm — минимальная длина волны неровности, либо непосредст­ венно по записи случайного процесса, так чтобы функция мало изменялась на интервале разбиения.

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

Обычно для

 

аппроксимации

корреляционных

функций

не­

ровностей и сил рабочих

сопротивлений

их

сперва нормируют,

т. е. делят на

максимальное значение

ординаты

Rx(0),

Ry{0),

Rxy(0),

а затем

подбирают

функциональную зависимость вида

Р(/*) = ^

 

п

 

 

 

т

 

 

 

=

У ^ А

^

cos ру* +

^

Л

sin р,| /• |,

 

 

 

i = l

 

 

i=n+l

 

 

 

где Ai,

a,i, Pi — неопределенные

коэффициенты.

 

 

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

р(/*) = Де - 0 '"* 1 cos ру* + А2е~а^

cos р2 /*

(54)

показано в работе [31].

 

 

Часто в первом приближении можно положить

 

Л1 = р0 (при /* = 0); А2 = 0; а 2 =

0; р2 =

0.

Если обозначить /* средний период прохождения через нули (по переменной /*) корреляционной фукции р(/*), то коэффици­ енты ai и Pi можно определить по следующим зависимостям:

a1 0 = - L _ l n

Е2

| / , |

Pi(при / = / , )

Рю -

— •

9*

131

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

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

Спектр показывает, какого рода колебания преобладают в данном процессе, какова его внутренняя структура.

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

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

 

 

 

оо

 

 

 

 

5((») = 2 j#(/*)cosco/*d/*,

 

(55)

где со = 2nf Ус-

 

о

 

 

 

 

 

 

 

 

Можно получить и обратную

зависимость

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

функции от спектральной

плотности

 

 

 

 

 

оо

 

 

 

 

#(/*) = - L

('S(co)cosco/*dco.

 

(56)

 

 

л

J

 

 

 

При /* =

0 имеем

 

о

 

 

 

 

 

00

 

 

 

 

 

 

 

 

 

#(0) = />= — Г S(<o)rf(B.

 

 

 

 

 

 

о

 

 

Формулы

 

 

л J

 

обратным

(55) и (56) называют

также прямым и

преобразованием Фурье.

 

 

 

 

 

Пользуясь преобразованием Фурье, можно

найти

спектраль­

ную плотность для нормированной

функции, заданной

формулой

(54). Не приводя преобразований, запишем результат

 

 

5(со) = 2/?(0)

\

а , ( ( о 2 + а? + р?)

 

 

 

 

 

( а > 2 - а 2 - в 2 ) 2 + 4

 

 

 

+ А

а 2

( с о 2 + а2 + Р2)

 

(57)

 

 

 

 

 

( а , 2 - а 2 - В 2 ) 2 + 4а2 2 со2

132

В дальнейшем понадобится также выражение для спектраль­ ной плотности ускорения, т. е. второй производной функции не­ ровности. Спектральная плотность второй производной от функ­ ции равна спектральной плотности самой функции, умноженной на со4:

5 у с к И = co4S(co). При этом необходимо обеспечить условие

оо

 

J SycK(co)d(u < оо.

(58)

о

 

Приведем формулы прямого и обратного преобразования Фурье для взаимных корреляционных функций и спектральных плотностей:

 

Sxy(v>)=

j

Rxu{l*)er№dl*\

 

 

 

 

— о о

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

Rxy(l*) =

^ -

J4 ,((D) e / M '*do .

 

 

 

 

— оо

 

 

 

 

 

Такой результат получен потому, что взаимная

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

ная функция — нечетная

функция,

а

взаимная

спектральная

плотность — комплексная

функция.

 

 

 

 

 

Если случайный

процесс

образуется

суммированием двух

(в общем случае и большим

числом)

стационарных и стацио­

нарно связанных случайных процессов х и у

 

 

 

z(l*)

=

x(l*) +

y(l%

 

 

 

то корреляционная функция процесса z(t)

равна

 

=

Rx{i*)

+ ВД*) + RxyV*) +

Ryxin

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

 

 

 

 

 

 

 

Sz(co) = Sx(tit) + Sy(a>) + Sxy{(o)

+

Syx((o).

 

Для некоррелированных процессов

*,(/*) =/?*(/*) + З Д * ) ;

Sz(et) = Sx(m) + Sy(ta).

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

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

133

Установим количественную связь между двумя способами

аппроксимации.

 

 

 

Для этого

используем метод

неканонических

разложений

стационарных

случайных функций [38], в соответствии с которым

случайную функцию q ( t ) можно представить

в виде гармоники

к \ sin озМ- %2 cos со^, где случайные

величины

к \ , К2,

со независи­

мы и удовлетворяют условиям

 

 

 

 

 

>

 

(59)

 

М[Я?] = М[Я|] = Я,(0), J

 

 

где символ М обозначает математическое ожидание случайной величины.

Законы распределения вероятностей случайных величин Xi и

Х2 произвольны.

Плотность

распределения случайной величины

со определяется

формулой

 

 

 

 

1

S(a>)

i

ии

(60)

 

 

/(СО):

 

J P(T)COS atxdx.

 

RAO)

 

Зная распределение частот гармоники, можно найти распре­ деление длин волн неровностей, если воспользоваться зависи­ мостью при v = 1 м/с

со = 2 л Л(о,

где / — длина волны неровности.

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

« 0 =

Тогда

dh f(h(l)).

dl

 

 

 

 

 

f(0 = ^

f

(

^

 

 

Плотность распределения

высот

неровностей

Н определяет­

ся через функцию распределения длин /

 

/(Я)

= I

dg(H)

ПёШ)],

 

 

 

dH

 

 

где g(H) — функция, обратная

по отношению к

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

ной зависимости между высотами и длинами

 

 

Н = а + Ы,

 

8(H)-

Н а I .

 

134

Тогда

При а = 0 имеем

Таким образом удалось связать два способа обработки мик­ ропрофилей полей.

Обратимся к конкретному анализу экспериментальных дан­ ных.

В табл. И приведены результаты обработки микропрофилей с помощью теории случайных функций. Измерения микропрофи­ лей во всех случаях проводились на участках длиной не менее 50—100 м.

Оценим высоты неровностей сельскохозяйственных полей по

классификации, данной в работе (40}.

 

 

 

 

К малоизношенным покрытиям ( У R(0)

^

1,5

см)

относят­

ся шесть участков, к сильноизношенным ( У R(0)

= 1,5 ч- 3,0 см)

двенадцать участков и к разбитым (j/i/?(0)

>

3,0

см)

—десять

участков.

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

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

По данным табл. 11 были построены в логарифмических ко­

ординатах спектральные плотности

ускорений

жесткого

катка,

движущегося по неровностям

(для

краткости: ускорений,

созда­

ваемых неровностями),

для

скорости

движения

трактора

v =

= 1 м/с. Затем набор

спектральных

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

вался как совокупность реализации случайных функций угловой скорости со. Для характеристики этих функций построены рас­ пределения ординат в трех сечениях, соответствующих to = 10; 25; 50 1/с. Во всех сечениях закон распределения функции S (со) оказался логарифмически нормальным с практически постоян­ ным коэффициентом вариации для всех сечений и математиче­ ским ожиданием и дисперсией, зависящими от переменной <о. Следовательно, функция S(co) является нестационарной случай­ ной функцией.

135

 

 

 

 

 

 

 

 

 

+ А2е~а^1*^

cos ру*, v = 1 м/с)

 

Таблица 11

 

Характеристики неровностей

(р = Ахе~ai1'*'

 

 

Автор

 

 

Д о р о г а (фон ,

рельеф)

 

 

YR(0), СМ

 

 

а,, 1/с

а 2 , 1/с

Pi. l.c

Я. М. Певзнер,

Булыжная, со впадинами и буграми

. . . .

2,50-3,28

0,85

0,15

0,50

0,20

2,00

А. А. Тихонов

Булыжная,

удовлетворительного качества . . .

1,35-2,24

1,00

0,45

»

 

 

 

 

 

 

 

 

 

0,85

0,15

0,20

0,05

0,60

»

 

 

 

 

 

 

 

 

0,50-1,24

1,0

0,15

Н. М. Антышев

Стерня озимой пшеницы, по направлению уборки

2,40

1,0

0,42

0,29

»

Стерня

озимой

пшеницы, против

направления

3,50

.—

1,0

0,53

0,33

 

 

 

 

 

 

 

 

 

 

Стерня

озимой

пшеницы,

против

направления

3,26

0,9

0,1

0,70

0,20

1,57

 

 

 

 

 

 

 

 

 

»

Лущеная

стерня

озимой

пшеницы,

поперек

на-

2,74

0,95

0,05

0,50

0,30

1,18

 

 

 

 

 

 

 

 

 

*

Вспаханная

стерня озимой

пшеницы,

поперек

на­

 

 

 

 

 

 

правления

уборки, дно

борозды

 

 

2,42

0,7

0,3

0,65

3,20

1,57

 

 

 

*

Вспаханное

поле,

поперек

 

направления вспашки,

4,09

0,9

0,1

0,50

0,40

6,48

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

»

Паровое

поле, поперек направления предшествую-

3,65

0,6

0,4

0,60

0,75

1,57

 

 

 

 

 

 

 

 

 

»Кукурузное поле, вторая междурядная продель­

ная

культивация, вдоль рядка

1,97

— 1,0

3,44

2,09

Кукурузное поле, вторая междурядная продоль­

 

 

 

 

 

ная

культивация, вдоль рядка, по колее трак-

3,36

1,0

 

3,20

1,57

 

 

 

»Вспаханное поле, имитация поперечной культива­

ции кукурузы, по колее трактора

3,10

0,75 0,25 0,86

0,30

2,36

Н. М. Антышев

Стерня

кукурузы

после

уборки

на

силос,

вдоль

2,28

1,0

0,71

0,785

»

Стерня

кукурузы

после

уборки

на

силос,

поперек

3,22

0,70

0,3

0,25

0,60

1,57

 

 

 

 

 

 

 

 

 

 

 

 

 

»

Грунтовая

полевая,

по

колее

 

 

 

 

2,12

1,0

0,58

0,63

»

 

 

 

 

 

 

 

 

 

 

 

 

0,815

1,0

0,13

1,05

А. А. Силаев

Полевая

и

вспаханный

луг, перпендикулярно бо-

15,20

—•

1,0

0,01-0,11

0,025-0,14

 

 

 

 

 

 

 

 

 

 

 

 

 

Б. Н. Кириенко

 

 

 

 

 

 

 

 

 

 

 

 

1,36-2,78

1,0

0,19—0,78

0,49—0,75

»

Стерня,

вдоль

направления

вспашки

 

1,81—2,02

1,0

0,18-0,47

0,33-0,425

»

Зяблевая вспашка, поперек направления борозд 2,45-2,80

-

1,0

0,75-3,60

0,43-0,69

»

Пропашное

поле,

занятое кукурузой

 

2,72-3,35

1,0

2,66—4,30

8,90

»

Пропашное

поле,

занятое

 

картофелем, вдоль

2,06 -2,23

1,0

0,55-0,61

0,42-0,76

 

 

 

 

 

 

 

 

 

 

 

 

 

В. Я. Анилович

 

 

 

 

 

 

 

 

 

 

 

 

5,18

1,0

1,30

 

Поперек

пахоты

при сплошной

культивации * .

4,20

1,0

0,29

5,40

 

Поперек борозд поля при уборке

пшеницы * . .

5,45

1,0

7,50

53,50

М. И. Ландсман

Хлопковое

поле

после

посева

 

 

 

 

1,29

1,0

1,20

А. И. Корсун

Хлопковое

поле

после

первого

полива

. . . .

1,35

0,5

0,5

1,20

1,20

2,10

»

Хлопковое

поле

после

второго

полива

. . . .

1,48

1,0

0,60

0,36

Дю Ин Ю

Картофельное

поле,

вдоль

борозд *

 

3,10

1,0

23,00

60,00

»

Картофельное

поле,

поперек

борозд * . . . .

2,80

1,00

42,00

60,00

»

 

 

 

 

 

 

 

 

 

 

 

 

3,46

1,00

22,00

100,00

»

 

 

 

 

 

 

 

 

 

 

 

 

4,60

1,00

16,00

81,00

»

Грунтовая

без

местных

неровностей * . . . .

4,55

1,00

21,00

71,00

»

Стерня

вдоль

борозд *

 

 

 

 

 

 

3,46

1,00

 

30,00

60,00

»

Стерня

поперек

борозд *

 

 

 

 

 

3,46

1,00

20,00

50,00

^

* Приведены данные д л я ускорений, создаваемых неровностями; е д и н и ц а измерения УR(0) м/с2 .

Такую функцию удобно представить в виде произведения 5 (со) = S0(co)p,

где So (со) —математическое ожидание S(co); р — случайная ве­ личина. Математическое ожидание параметра р равно единице.

Дисперсия р постоянна и равна коэффициенту вариации функции S (со).

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

W(p) =

W(SoP)S0.

 

 

На рис. 80 дан график

функции W(p).

Функция W(p)

не за­

висит от угловой скорости

о). Плотность

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

W(p)

и(Ро)

 

 

 

 

Ро

 

Р'о

 

 

0

0,2

0,4

0,6

0,8

р

Рис. 80.

Плотность

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

параметра

р

 

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

Это свидетельствует о том, что среднее значение функции 5(со), равное S0(co) и соответствующее М(р) = 1, имеет малую вероятность появления (0,2). В то же время наиболее вероят­ ным значением параметра р будет 0,05. Этот результат ставит под сомнение возможность ограничиваться в расчетах подрессорных систем лишь средним значением спектральной плотности; необ­ ходимо принимать во внимание функцию распределения вероят­ ностей и определять наиболее вероятное значение показателя плавности хода трактора.

Для расчета колебаний тракторов небходимо также найти аналитическую формулу функции So (со). График аппроксимирую­

щего

выражения, удовлетворяющего

формуле

(57)

при А2 = 0

и А\ = 1, приведен на рис. 81, а. Предложенная

аппроксимация

функции

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

участке

угловой

скорости

(0—40)

1/с хорошо

совпадает

с функцией 50 (со),

а

при

со >

> 40

1/с

она описывает ее в среднем. Коэффициенты

математи­

ческого

ожидания

функции

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

плотности

равны

R(0)

= 26,5 м2 4 ; р = 39 1/с;

а = 8,9

1/с.

 

 

 

 

 

138

Удобство представления функции S (<о) в виде произведения 5о(со)р состоит в том, что найденная совокупность случайных функций обеспечивает аналитическую аппроксимацию, которая удовлетворяет условию существования функций спектральной плотности.

Коэффициенты R(0),

а, р получены для скорости

v = 1 м/с.

Для любой другой скорости движения машины

[32,

2] R\{^)

=

= 26,5и2 м2 4 ; «i = av

1/с; В =рЧ>

1/с. Предложенная

аппрокси­

мация удобна при расчете ускорений подрессоренных

масс,

по­

скольку при расчетах

возможно

применение

табулированных

S0(U),M'/CS

 

S„'(CJJ,CMZ/C

 

 

 

О

20

40

60ц1/с

0

4

8 со, 1/с

 

 

ч)

 

 

 

6)

Рис. 81. График функции S0(co) для расчета:

а — ускорений подрессоренных масс; б — перемещений

функций. Для того чтобы оценить порядок коэффициентов полу­ ченного аппроксимирующего выражения спектральной плотно­ сти, по данным табл. 9 построены корреляционные зависимости между R(0), а, В (рис. 82, а). На эти графики нанесены точки, соответствующие коэффициентам аппроксимирующего выраже­ ния функции спектральной плотности.

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

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

139

При расчетах подвесок тракторов возникает необходимость также рассчитывать и перемещения подрессоренных масс.

Расчеты, аналогичные вышеприведенным для ускорений воз­

действий, позволили

получить графики плотности

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

параметра р 0 '

(см. рис. 80), обобщенной

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

сти неровностей для

расчета

перемещений

подрессоренных масс

S 0 ' ( M )

(см.

рис.

81, б)

и корреляционных

зависимостей

(рис. 82,

б).

 

 

 

 

 

ос, 1/с

 

 

ос, 1/с

 

 

О

20

40

60

80 Д 1/с о

г

4

6

8 р, 1/с

0

10

20

30

40 а, 1/с 0

2

4

6

8 «, 1/с

 

 

 

а)

 

 

 

6)

 

Рис. 82. Корреляционные зависимости между коэффициентами аппроксимирующих выражений спектральной плотности для расчета:

а — ускорений подрессоренных масс; б — перемещений

График обобщенной спектральной плотности аппроксимиру­ ется общей формулой (91) при следующих значениях парамет­ ров:

Я(0) = 8,56 см2 ; а = 0,995 1/м; 6 = 0; Л 2 = 0; Л 1 = = 1 .

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

Представляет интерес выяснить физический смысл и некото­ рые количественные зависимости для коэффициентов а, р\ R(0)

140

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