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

Некоторые проблемы, возникающие при создании многомассовой математической модели движения автомобиля и методы их решения Текст научной статьи по специальности «Механика и машиностроение»

CC BY
151
39
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МНОГОМАССОВАЯ МОДЕЛЬ ДВИЖЕНИЯ АВТОМОБИЛЯ / УРАВНЕНИЯ ДВИЖЕНИЯ НЕПОДРЕССОРЕННЫХ МАСС АВТОМОБИЛЯ / КИНЕМАТИЧЕСКИЕ ХАРАКТЕРИСТИКИ ПОДВЕСКИ / MULTIMASS MODEL OF VEHICLE MOTION / EQUATIONS OF MOTION OF SPRUNG MASS OF A VEHICLE / KINEMATIC CHARACTERISTICS OF SUSPENSION

Аннотация научной статьи по механике и машиностроению, автор научной работы — Бахмутов С. В., Гусаков Д. Н.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Problems on Creating of a Multimass Mathematical Model of Vehicle Motion and Methods of their Solution

The paper analyzes features of vehicle multimass model acquisition, regardless of its design and building complexities. The authors propose the technique of writing the equations of motion of sprung mass in the general moving coordinate system and describe the role of kinematic characteristics of the suspension.

Текст научной работы на тему «Некоторые проблемы, возникающие при создании многомассовой математической модели движения автомобиля и методы их решения»

Раздел 1. Наземные транспортные средства энергетические установки и двигатели. место. И как показывает анализ, если руководствоваться рекомендациями работы [3], то исходное отношение натяжений должно быть т = — = 8,33 (здесь ^ и - силы натяжения

ведущей и ведомой ветвей ремня). Окончательно с некоторым запасом можно рекомендовать для передач с подпружиненными натяжными шкивами т = 8, а не т = 5, как это принято в документах 1БО [5]. Это обеспечит увеличение ресурса ремня не только без снижения тяговой способности передачи, но и с некоторым ее увеличением. Следует отметить, что, несмотря на значительное количество работ, посвященных самонатяжным ременным передачам, ремень в них, как правило, рассматривается в виде гибкой нерастяжимой нити. Как показали многолетние исследования В.К. Мартынова [6], получить современную и точную методику расчета и проектирования ременной передачи не представляется возможным. Кроме того, требуются дополнительные экспериментальные исследования для установления оптимальных норм предварительного натяжения ремня для рассматриваемого типа передач.

Литература

1. Терешкин Л.В. Приводы генераторов пассажирских вагонов. - М.: Транспорт, 1990, 152 с.

2. Самошкин С.Л., Доронин И.С., Чернышов А.А. Приводы генераторов индивидуальных систем энергоснабжения вагонов локомотивной тяги: Обзор. - М.: ЦНИИТЭИтяж-маш,1986. 44с.

3. Баловнев Н.П. Анализ методов расчета клиноременных передач. В сб. Агрегатирование и приводы сельскохозяйственных машин. М.: НПО ВИСХОМ, 1985 с. 60-71.

4. РТМ23.2.90-85. Методика проектирования и расчета клиноременных передач сельскохозяйственных машин с натяжными подпружиненными шкивами, обеспечивающими постоянство натяжения ведомой ветви. - М.: ВИСХОМ, 1984. - 34 с.

5. Пронин Б.А., Овчинникова А.А. Расчет клиноременных передач // Вестник машиностроения. 1982. № 3. с. 23-26.

6. Мартынов В.К. Разработка теории, методов расчета и проектирования современных передач трением гибкой связью. Автореферат дис. доктора техн. наук: 05.02.02 Тула. 2004 г.

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

д.т.н., проф. Бахмутов С.В., Гусаков Д.Н.

МГТУ «МАМИ»

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

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

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

ции массы, с которой жестко связана система координат, остаются постоянными.

В случае многомассовой системы возникает проблема, обусловленная вычислением моментов инерции масс, центр которых не совпадает с центром выбранной подвижной системой координат. Логично привязать подвижную систему координат к кузову автомобиля (неподрессоренной массе), в центре масс автомобиля. Таким образом, остается задача вычисления моментов инерции неподрессоренных масс. Эта задача может быть решена двумя способами. В первом случае координаты неподрессоренных масс (колес) остаются независимыми, к ним привязываются индивидуальные системы координат, а связь систем осуществляется посредством углов Брайнта. Данный вариант неудобен тем, что на каждом шагу интегрирования необходимо пересчитывать все силовые и координатные факторы, входящие в уравнения, с учетом изменяющегося взаимного положения систем координат. Кроме того, несвязанные системы при определенной погрешности интегрирования могут «разойтись» -т.е. после возвращения системы в состояние покоя взаимные линейные и угловые координаты систем могут не соответствовать начальным. Второй вариант предполагает использование единой системы координат. В этом случае силовые факторы, как внешние, так и внутренние, входят во все уравнения системы в едином виде. Однако в таком случае моменты инерции неподрессоренных масс являются переменными и должны некоторым образом пересчиты-ваться на каждом шаге интегрирования.

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

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

Рис.1. Расчетная схема трехосного автомобиля и внешние силовые факторы,

действующие на него.

Для переноса моментов и произведений инерции в параллельную систему координат из

I, = х' I х . + т х (у12 +

I, = у I + т х(х'2 + %12 у

I , = I + т х (у12 + х'2

Рх,у, х, у , = Р + т х х'х у' ху

Р, , х % = Р + т х х'х?

(1.1.а) (1. 1.б) (1. 1в) (1. 2.а) (1. 2.в)

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

представлена при помощи коэффициентов (оси х, у, % - первоначальные, а х , у , %' -повернутые):

Таблица 1.

Новые координаты Первоначальные координаты

х у 2

х' а1 в 71

У а2 в 72

2Х а3 в 73

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

х' = а1 х х + Д х У + ^ х % (координаты у' и 2' вычисляются аналогично)

Обратный переход осуществляется при помощи коэффициентов в вертикальных строках:

х = а, х х'+о, х у'+а3 х г' . , >

1 2 ^ 3 (и т.д.)

Моменты инерции в повернутой системе координат в общем виде выразятся: I , =а2 х I + в2 х I + г2 х I - 2 хДхг.х Р -

х 1 х у /1 % #^1/1

-2х а х в хР -2хахг.хР

1 1 ху 1 1 х%

I , = а22 х I + в х I + 722 х I - 2 х Д хг2 х Р -

у 2 х ~ 2 у '2 % ~ 2 ' 2

- 2 ха2 х в х Р - 2 ха2 хг2 х Р

2*2 ху 2 ' 2

^ = Оз2 х 4 +вз2 х Iy +Г32 х / - 2 хДз хГз х Ру2 -

- 2 ха3 х в х Р - 2 ха3 ху3 х Р

3/^3 ху 3*3 хг

Произведения инерции относительно повернутых осей:

Рху =-а1 ха2 х ^ -Д х в2 х Ь -П хП х ^ + Рху х(а хД +а2 хв1 ) + + Рхг х(а1 х ^2 +а2 х71 )+ Руг х (Д ^2 + Д х ) Рх^ = -а1 ха3 х Ix -Д хв3 х ^ -Г1 хГ3 х ^ + Рху х(а1 хв3 +а3 хД ) +

+ Р„ х(а1 х^3 +а3 хГ1 )+ Р х(Д х^3 + Д х ) Ру, =-а2 ха3 х Ix -Д хД3 х ^ -72 х73 х ^ + Рху х(а2 хв3 +а3 хД ) + + Рхг х(а2 х ^3 +а3 х ^2 )+ Р х(в2 х Уъ + Д ^2 )

(2.1.а) (2.1.б) (2.1.в)

(2.2. а) (2.2.б) (2.2.в)

принятые обозначения см. в приложении 1 к статье.

7

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

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

йг

( дт Л

дТ ди дЯ _

= О,

(3)

) дцг дцг дС

Согласно выбранному в модели расположению осей координат (рис. 1.), кинетическая энергия подрессоренной массы:

Т = — ■ V дт, 2 ^ 1

( + г,-ц/■ у, }+(Гу + ( - х, + р■ у1У

■!>■г, + (■ х1 У +

(4)

С учетом шести обобщенных координат ( С = х, С = у, Сз = г, С = р, С = Г, Сб = ( ), характеризующих положение жесткого тела в пространстве, левые части уравнений Лагран-жа для подрессоренной массы будут иметь вид:

йг йг

( дТ Л

дТ„

дС1 ) дС

= тп

{дТ Л

дТп

дС[ 2 ) дС2

= т,

о

й

2 ) ( дТ0 Л

дТо

йг ^ дС3 ) дс3 Г дТ0 Л дТ0

= т,

о

Гх-у ■Гу + Г ■ V-((2 + Г2 )■ Хо + Г р-()■ уо +

+ ( + (/■ р\го _

Гу -р■ Г +(■ Гх +((+р■ г) ■ Хо -(у2 + р2) ■ Уо +

+ (( Г - р)го

К - Г ■ Гх + р Ку +(рр ( - Г') ■ Хо +(р + ^ (У уо -

-(2 + Г2 У го

1 о

!х о ■ р + - г V + р ■ Гу ))• уо +(-К - УуГх ))о- го +

йг ^д<С4 ) дс4

+ (Р ■( - Г)Рхуо -(Р■ Г + (УРхго +((2 - Г 2 )руго + ^ ( (о - /уо )

(5)

йг

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ГдТо Л

дТо

Л о ■ Г + ( К - р Ку - Г )) хо + ( -(■ Гу + Г-Уг)•

то го -

\дЧ5 ) дС5

((■ Г + рУРху о +( 2 2 )-Рхг о + (( р-()-Руг о + р^ (■ ((хо - /го )

й )~1ТГ = 1г о (+( К + Г.у - р Кг)) хо +( ку - Г К2 - Гх)) уо +

+ ( 2 - р 2\Ру о +((■ Г - рУРг о -(((■ р + Г)Руго + р■ (/у о - ^х о )

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

дС, ч

(б)

т i _ ^ ,

кС1С1 =

дС,

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

ям в пространстве.

Абсолютные скорости элементарной неподрессоренной массы осей х, у , ^ найдем в соответствии с общими формулами (4):

5тп-

по направлению

X; = ¥х + г-(( + кг2т + кр2-р + к22-2 )-? + к2? V - кр?- р - к?- г)--(( 0 -к?-Г - кр-р - к?-2)+ к--Уг + к™-Г + крх-р

У, = УУ-(кРР-Р - к7-у; - к7-г)-(2^0 +кГ-г+кР2-р+к22-2)+

+ ( + к2? ■ уг - кр? - р - к? - г) (0 + кгх •г + крх- р + к2х - 2)- к] -Г -

(7)

- кру-р - к"-У,

2'; = к22 -Уж - Г- (+к- т + крх-р + к 2х- 2)+ ку • р - к]Р - У, - к]р - г) --( ; 0 -к?-г - кру-р - к2* -2)+ к- -г

Конечные формулы, определяющие законы движения подрессоренных и неподрессо-ренных масс по каждой из координат, соответствующих степеням свободы, получаются путем дифференцирования общих формул (4) по соответствующим обобщенным координатам.

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

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

Рис. 2. Силовые факторы, возникающие при разрыве подвески.

Выводы

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

Приложение 1.

Обозначения и индексация переменных, входящих в уравнения: 1х, 1у, Iг - моменты инерции в первоначальной системе координат;

1Х,, 1у,, 1г, - моменты инерции в новой системе координат;

х', у , г' - соответствующие смещения осей новой системы координат относительно первоначальной;

Р, Руг, Рг - произведения инерции в первоначальной системе координат;

Рх'у' , Ру'г' , Рх'г' - произведения инерции в новой системе координат;

Т - кинетическая энергия системы; и - потенциальная энергия системы; Я - диссипативная функция Релея; О, - обобщенная сила;

С, - обобщенная координата;

С - обобщенная скорость;

1Хо, 1уо, 12о и Руо, Рхго, Руго - моменты инерции и произведения инерции подрессоренной массы относительно осей х, у , г; хо , уо , го - координаты центра подрессоренной массы;

Гх, Гу, Гг, р , Г, ( - линейные и угловые скорости по соответствующим осям координат;

С, - возможное перемещение подрессоренной массы;

- возможное перемещение неподрессоренной массы; 1 - номер неподрессоренной массы;

Р№ - сила аэродинамического сопротивления;

Яг1, Яп1, Яг1 - продольная (боковая, вертикальная) реакция 1 -го колеса;

- эквивалентный силовой фактор, вводимый при разрыве подвески, действующий в

направлении соответствующей оси; Оа - вес автомобиля;

Ма - обозначение центра масс автомобиля;

Ыц - крутящий момент 1 -го колеса;

Мд1 - момент сопротивления повороту 1 -го колеса;

Ыср1 - момент сопротивления развалу ] -го колеса.

Математическая модель качения эластичного колеса по деформируемому

грунту

к.т.н., с.н.с. Чистов М.П., Наумов А.Н. МГТУ «МАМИ», 21 НИИИАТМО РФ

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

i Надоели баннеры? Вы всегда можете отключить рекламу.