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

книги из ГПНТБ / Гоберман Л.А. Прикладная механика колесных машин

.pdf
Скачиваний:
54
Добавлен:
24.10.2023
Размер:
9.73 Mб
Скачать

Используя далее уравнения Лагранжа и выражение (Х.1), приходим к следующим дифференциальным уравнениям:

 

d 2x

_ J ___ in

 

 

d z .

 

 

 

~JF

 

m T +

m 0

 

d t

W Sin(P

 

 

 

 

- f (/„ Sin Ф +

2 COS ф)

Ь

 

(X.3a)

 

 

 

 

 

 

 

 

 

 

 

 

c ß y

 

1

 

 

JL

d z

 

 

 

ЧІ2

 

mT +,

m0-!q

W C0S(P

 

 

 

d t

 

 

 

 

 

+ (^o cos ф — г sin ф)

J;

 

(X.36)

 

d2z

-

 

(n

mn

d 2x

5іпф

d 2y

 

 

 

1

ж

^ С ° 8ф

 

 

d t 2

m 0

(Ѵг

 

 

 

 

 

 

 

 

 

 

 

 

I d

 

/ d y y i ) .

 

(Х.Зв)

 

 

 

 

+

dt2

2 V d t )

J / ’

 

d 2ф

^ + «

0 (/2

 

( Q r - mn (/0 sin ф - f

2 cos ф)

d 2x

~dF -

+ 22) 1

Чф

 

 

 

 

H i 2

 

(l0cos ф -

г sin Ф)

d 2y

+

d 2z

+ 2г -§- • -g-

}• (Х.Зг)

В этих

уравнениях Q*,

 

Qz

и

— обобщенные силы на

соответствующих

обобщенных перемещениях

системы

 

 

 

 

дП

Qz

 

дП

 

дП

 

 

Qy = -

ду ;

 

d z

 

дц

 

Величина Qx численно равна суммарной касательной реакции поверхности качения на колесах машины.

Одним из параметров, определяющих положение системы, яв­ ляется указанная на рис. 57 величина w lt которую можно найти из соотношения

^ = х - ( т ^ Ѵ +у1ёф)-

(Х-4)

Для определения положения системы по параметру wt нас будет интересовать только знак этой величины: при wt 0 задние опоры еще не сорвались с препятствия, при > > 0 задние опоры сорвались с препятствия.

Положение системы определяется также наличием или поте­

рей

контакта

ее опор с

поверхностью качения.

Если

б1<;0

и

б2 <

0 , то это

означает,

что и передние и задние

опоры системы

контактируются с поверхностью качения; если

> 0

и б2 >

0 ,

то те и другие опоры не имеют контакта с поверхностью качения. Как видно из уравнений для öj и б2, знаки этих величин опре­

278

деляются

знаками

их

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

через w2

и w3.

w2 =

У +

h '— L x sin cp — h4 cos ф rKcos ф ;

 

(Х.5)

 

 

 

w3 = у

+ L 2 sin ф

cos ф rKcos ф ;

 

(Х.6 )

тогда

 

 

 

 

 

 

 

 

 

 

 

бх = - 3 - ; б2 = - 3 _ .

 

 

 

 

 

1

COS ф ’

й

COS ф

 

 

Очевидно,

что при

w2 < 0,

w3 <

0 передние и задние опоры

системы

находятся

в

контакте

с

поверхностью качения;

при

ш2 >■ О,

w3 >> 0 их контакт потерян.

 

 

После срыва задних опор с препятствия их деформация опре­

деляется

выражением

 

 

 

 

 

 

 

 

 

б -

Ш4

 

 

 

 

 

 

 

£

COS ф

 

 

 

где

 

 

 

 

 

 

 

 

 

W,

+

h) + L 2 sin ф — h4 cos ф — rKcos ф.

 

(X.7)

По аналогии с предыдущим

при

w4 •< 0 задние опоры

после

срыва их с препятствия имеют контакте основанием, при

>• 0

не имеют контакта.

Таким образом, параметры wt (і = 1ч-4) вполне определяют положение исследуемой системы. Назовем эти величины пара­ метрами положения системы. В зависимости от знака этих пара­ метров возможны следующие ситуации (рис. 57).

До срыва задних опор с препятствия:

 

1 ) Wo <

0 ;

w3 <

0 — передние и задние опоры имеют контакт

2 )

w2 >

0 ;

w3 <

 

с поверхностью качения;

0 — передние опоры потеряли контакт, зад­

 

 

 

 

 

ние опоры находятся в контакте с по­

3)

w2 0 ;

w3 >

 

верхностью качения;

0

— задние опоры потеряли контакт, перед­

 

 

 

 

 

ние опоры находятся в контакте с по­

4)

щ2 >

0;

w3 >

 

верхностью качения;

0 — передние и задние опоры потеряли кон­

 

 

 

 

 

такт с поверхностью качения.

После срыва задних опор с препятствия:

 

5)

w4 > 0 ;

w2 >

0 — задние и передние опоры не имеют кон-

и>4 >

0 ;

w2 0

такта с поверхностью качения;

6 )

— задние

опоры не

имеют контакта, пе-

 

 

 

 

 

рідние

опоры — в

контакте с поверх­

7)

w4 <

0;

w2 <

0

ностью качения;

 

передние и задние опоры имеют контакт

споверхностью качения;

8)w4 < 0 ; w2 > 0 ' задние опоры имеют контакт с поверх­

ностью качения, передние колеса не имеют контакта.

279

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

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

Рис. 57. Возможные положения машины при съезде с препятствия

Для системы с передними и задними ведущими колесами (все опоры — движители) ,

Qx = сг (бг + б2) у;

для системы с передними ведущими и задними ведомыми колесами

Qx = ст(бху — б2/к);

для системы с задними ведущими и передними ведомыми колесами

Qx = сг (б2у — бJ K);

для системы со всеми ведомыми колесами

 

Qx = — сг (бг +

б2) /к,

где у — коэффициент, зависящий от

величины тягового усилия

на движителях; у — срСІ1£;

 

Фсц — коэффициент

сцепления движителей с поверхностью ка­

чения;

определяющий

использование силы сцеп­

I — коэффициент,

ления; /к — коэффициент сопротивления качению.

280

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

(/к = 0 ; Фсц = 0 ).

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

Примем следующие обозначения:

 

FK. 1 — коэффициент связи передних опор

(FKX — 0; fK>ФсцЕ)‘>

FK2 — коэффициент связи задних опор

(FK2 = 0 ; /к‘> фсцІ);

Кг =

K2=

>;

II

со

cr

COS ф

 

 

Cr

показатели, зависящие от жесткости опор;

cos2 cp

Cr

 

 

cos3 q>

)

 

 

 

D x =

l0 sin cp + z cos cp; D3 =

+ h) sin ф —Lx;

D 2

== l0 cos ф — г sin ф; D4

— г/ sin ф + L2.

С учетом этих

обозначений

 

 

 

<Эхі =

К і Щ Р К 1;

Q x2 = K ^ g F K z ,

Qx3 =

KiW^FKz-,

 

 

Qyo = — (m, + m0) g\

 

Qyl =

K-2W2\

Qy2—

 

Qy3 — ---

 

 

Qq,o = — m0gDs -f MpK;

 

Qfpi —

KsD3w2',

Qp2 =

KsDjWg,

Qrpi =

K3D3w3,

где Л4рк— реактивный

момент

от двигателя;

вводят в расчет

при исследовании бокового опрокидывания системы. Обобщенная сила

Qz = —с0 (г — hr) — /n0g cos ф.

Значения остальных обобщенных сил для различных ситуаций приведены в табл. 16.

Решение системы дифференциальных уравнений на ЭВМ.

Исходную систему дифференциальных уравнений преобразуем следующим образом.

Обозначаем

d x = vx-

d y

d z

d(p

CO

d t

d t

~ d f ~ V z ’ ~ d f

 

457-

281

282

Ситуации и обобщенные силы системы

 

 

­

 

 

 

 

 

 

 

ситу№ ации

Параметры положения системы

 

 

 

1

w4 <

0;

w2 <

0;

w3 <

0

 

2

w4 <7 0;

w2 >

0;

<

0

 

3

w4 < 0 ;

ay2 <

0;

wg >

0

 

4

ffi'i <

0;

a>2 >

0;

<jy3 >

0

 

5

w4 >

0;

>

0;

w4 >

0

 

6

1% >

0;

да, <

0;

w4 >

0

 

7

> 0 ;

ш2 <

0;

<

0

 

8

> 0 ; w2 > 0; ш4 < 0

Qx

QXi + Qx2

qX2

Qx>

0

0

*3*1 + ^*3

Обобщенные силы

Qy

Qy, +

+ Qy*

Qy, + Qy*

Qyo + Qyi

Qyо

Qy.

Qyо + ^

Q</0 +

+ <3,3

Qy, + Qya

Т а б л и ц а 16

^cp

<2фо +^4)3 + ^Фз

в<2

+ о

Сфо + <3фі

*Зфо

*Зфо

Сф„ + <гф1

Q<po + ^ф. + *Зф3

•Зфо + ^Фз

и приводим уравнения (Х.З) к виду

 

 

dv

~ ^

0 1

^

A n ,-^ - + A03^ f -

1

dt

л ° 2

 

dt

 

 

 

dvy _ р

 

1

р

da)

I

duy

 

 

В01 “г ^*02

~dt

В°3 ~W

 

 

Öl

 

^ 02

(X.8 )

dVz __I

r>

 

 

rfco

.

r

dvx

. r düy

- М — Ьплoi -b(^ 0 2

~äf +

b03

 

Ar ^Oi ■

 

 

 

 

 

 

 

 

 

dv,

 

~ W ~ ° 01 + D ° 2 X F +

D os~4t~ + D °4 ~W ’

 

где

1

 

 

 

 

 

 

 

 

 

 

(Qx m0u>2D22m0vz(£>cos cp);

Аoi m0+ mT

 

 

A,0 2

 

 

m0Dl

 

 

 

 

 

Щ+ mT

'

 

 

 

 

 

 

 

 

 

юз

 

 

 

m„

 

 

 

 

 

------r ---- sin cp;

 

 

 

 

 

m0-\- mr

 

 

 

1

 

(Qy +

rn0(ü2D1-f 2m0ozco sin

cp);

B0i m0 + mx

 

y

 

 

 

 

 

 

B02

 

mn

rnT ■D2;

 

 

m0+

 

 

R

—_

 

mo

 

 

 

 

08

 

m0+ mT cos cp;

 

 

COl = —

{Qz A-m0zu)2)-

 

Cq2

Ab0>

°C,0 3

 

•sin cp;

C0i — cos cp;

 

A n

 

 

 

 

2/n0ZHK(ö);

 

 

A n —

 

 

 

m0Di

 

 

 

 

^

—mo(/o + z2)

 

 

 

 

 

 

 

А>з —

 

 

—mo( / 0 + z2)

 

 

 

 

 

 

 

 

A i —

 

 

 

m0l

 

 

 

 

 

 

’rao(^o + z2)

 

 

 

 

 

 

 

 

Из уравнений (X.8 )

исключаем

 

 

 

 

düz .

 

dvx

_

 

_

rfco

 

 

~~5Г ’

 

" Ж ’

Ж ’

~ЗГ’

 

283

Начальные условия, соответствующие начальному моменту срыва передних опор с препятствия, находим из общих выражений для обобщенных сил Qz, Qy и (2Ф, которые в этом случае прирав­ ниваем нулю, а обобщенные перемещения принимаем соответ­ ственно равными ги, ун и срн. Начальное значение обобщенного перемещения х = хя находим из геометрических соотношений. В результате получаем следующие соотношения, из которых определяем начальные значения указанных обобщенных пере­ мещений системы:

m0g cos фн -f с0 (гя 1г2) = 0 ;

(тт+ m0) g cos2 срн + сг(уя —

sin срн — (hx + rK) cos cpH] +

 

 

+ cr[Ун + L2sin фя — (hx+

rK) cos cpH] = 0 ;

m0gl0cos cpH— m0^i>sin <p„ +

{cr [yu Lxsin <p —

(Ai + Л<) cos фн] (уя sin cpH— Li) +

cr [ya -I- L%sin cpH—

 

 

— (Ai 4- rK) cos фн] (уя sin Фн + L2)) = 0;

 

 

Хн = — Ljcos фи +

(ун __ i Lsln фи).

Из этих

выражений

находим

 

 

 

z„ =

1

tn.cs

 

 

 

h2---- -^соБфн;

 

 

 

 

С0

 

Уп = 1

у-7- sin ф„ 4- (hi 4- rK) cos фн — (Jüi±^S-) COS фн;

 

 

F = m0g cos3 ф(/„ cos фн — zHsin фн) +

+

сг [УоLi sin фн(hi +

гк) cos фн] (уя sin фн — Lx) f

+

сг [уо+ L2sin фя — (hi 4- rK) cos фн] (уяsin фн 4- L2);

хя = — Li cos фн 4- (уя Li sin фн) tg фн.

Для решения дифференциальных уравнений на ЭВМ теперь остается установить условия опрокидывания системы, выбрать метод интегрирования и на основе этого составить блок-схему программы.

Условия опрокидывания машины могут быть определены по

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

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

Х1С— X

costp

+ A) tg ф

 

 

284

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

: X

COS ф

+ A) tg ф

 

 

Этими выражениями определяются необходимые условия опрокидывания.

Однако при исследовании процесса опрокидывания на ЭВМ

представляет

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

изменения функций ф (і)

и ф (t), т. е.

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

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

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

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

можно проследить и по кривой ф (/).

В качестве условия опрокидывания машины можно также при­ нять условие I ф I ^ 1,5.

При решении уравнений на ЭВМ «Минск-22» была принята про­ грамма с переменным шагом интегрирования, в основу которой положен метод Рунге—Кутта.

Результаты, полученные при решении данной динамической модели на ЭВМ, позволяют сделать следующие выводы:

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

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

устойчивой;

вероятность опрокидывания машины возникает

наибольшая

в тех случаях,

когда при съезде с препятствия осуществляется

ее торможение или выключение двигателя фрикционом; наимень­

шая опасность

опрокидывания — при съезде машины с раз­

гоном;

высота препятствия уменьшается с уменьшением

критическая

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

Эти выводы были проверены экспериментально и получили свое подтверждение.

285

48. ОБОБЩЕННЫЕ МЕХАНИКО-МАТЕМАТИЧЕСКИЕ МОДЕЛИ ОПРОКИДЫВАНИЯ СЛОЖНЫХ КОЛЕСНЫХ СИСТЕМ

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

Движение таких систем может быть описано обобщенными ма­ тематическими моделями, которые получаются на основе учета шести обобщенных перемещений — трех поступательных и трех вращательных (ф— относительно вертикальной оси Cz; %— отно­ сительно продольной оси Су и ф—относительно поперечной оси Сх).

Неголономность системы в обобщенных математических моде­ лях учитывается уравнениями кинематических связей колеса с поверхностью качения; эти уравнения для колеса с пневмати­ ческой шиной были получены нами в разделе 15.

На рис. 58—61 представлены расчетные схемы для исследова­ ния процесса опрокидывания: погрузчика на базе трактора с шар­ нирно-сочлененной рамой (рис. 58); тягача с шарнирно-сочленен­ ной рамой и прицепа с передними управляемыми колесами (рис. 59); того же тягача и прицепа с поворотной передней тележ­ кой (рис. 60); того же тягача и прицепа с поворотной передней тележкой и подрессоренными передними и задними колесами

(рис. 61).

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

Исследуемую систему представим состоящей из семи масс:

тх — массы балансирной (грузовой) части трактора;

т2 — массы небалансирной (моторной) части трак­ тора;

т3, /л4, тъ, me — массы колес;

т 7 — суммарной массы навесного оборудования и груза, приведенной к центру захвата (ковша).

Обозначим также:

Jcixli Jayl и Jcizl — моменты инерции грузовой части трак­ тора соответственно относительно осей *і, Ух и Zj, проходящих через центр тяжести Сх этой части машины;

JiXl, J2 z, — моменты инерции моторной части трак­ тора относительно осей х х и zx;

286

Л Уг — момент

инерции моторной

части трак­

тора относительно оси у 2,

совпадающей

с осью горизонтального шарнира рамы;

J іХі, J іУі, JiZl — моменты инерции колес

= 3, 4, 5, 6)

относительно осей х г, у ІУ гг;

J с , х 0 J c i X l , J c b Xl , J c , Xt — моменты

инерции

колес

относительно

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

центры масс соответствующих колес;

гк — радиус вращения

колеса;

 

сп— жесткость навесного оборудования, при­ веденная к центру захата (ковша).

Упругие свойства колесных шин в радиальном направлении будем учитывать жесткостью сші (і = 3 - f - 6) и коэффициентом демп-

Рис. 58. Расчетная схема погрузчика с шарнирно-сочлененной рамой

фирования km, а в касательном и боковом направлениях — перемещениями т](. и | г, как это было принято в разделе 15.

Процесс опрокидывания машины с шарнирно-сочлененной ра­ мой, имеющей вертикальный и горизонтальный шарниры, удобно расчленить на два последовательных этапа. Первый этап — от­ носительный поворот вокруг продольных осей балансирной и не­ балансирной частей машины на углы уѵ1 и у2 до момента, пока

287

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