Архитектура Аудит Военная наука Иностранные языки Медицина Металлургия Метрология
Образование Политология Производство Психология Стандартизация Технологии


АНАЛИТИЧЕСКИЕ НАУЧНЫЕ ИССЛЕДОВАНИЯ НА АВТОМОБИЛЬНОМ ТРАНСПОРТЕ



Тема 9. Изучаемые вопросы: Структурная схема объекта исследования; Установление причинно-следственных связей; Составление математической модели исследуемого процесса; Написание и отладка программы; Проверка адекватности и настройка математической модели; Моделирование исследуемого процесса; Анализ результатов и формулирование выводов.

Аналитическое исследование проводится с использованием расчетных методов, с использованием ЭВМ. Аналитическое исследование, как правило, проводят параллельно с экспериментальными исследованиями. Это позволяет, во-первых, значительно сократить трудоемкость и материальные затраты на проведение исследований. Во-вторых, точность проведенных расчетов показывает, насколько глубоко и полно изучил автор объект исследования.

Типовая структурная схема процесса аналитического научного исследования в области эксплуатации автомобилей представлена на рис. 28.

1) Перед началом экспериментального исследования формулируют цель и задачи аналитического исследования.

Рис. 28. Структурная схема аналитического научного исследования

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

В качестве примера, на рис. 29 приведена структурная схема объекта исследования – процесса диагностирования автомобильной тормозной системы с пневматическим тормозным приводом. Данная схема разработана на основе структурной схемы процесса функционирования объекта исследования, представленной на рис. 1. В схеме тормозная система автомобиля, его подвеска подрессоренные и неподрессоренные массы, колеса с эластичными шинами, тормозные механизмы и пневматический тормозной привод, а также дорога, трансмиссия и ДВС представлены в виде 7-и отдельных блоков.

3) На следующем этапе аналитического исследования устанавливают причинно-следственные связи между блоками схемы. Взаимное влияние между этими блоками показано стрелками, с указанием параметров, посредством которых это влияние осуществляется.

 

Рис. 29.Структурная схема процесса диагностирования тормозной системы автомобиля с пневматическим тормозным приводом.

 

Так, например, блок «пневматический тормозной привод», имеющий давление сжатого воздуха в системе, равное Рвх, воспринимает управляющее усилие Uупр водителя на педаль тормоза. Это входные параметры блока «пневматический тормозной привод». Выходными параметрами блока «пневматический тормозной привод» являются: давление Рij сжатого воздуха в камерах тормозных механизмов; темп его изменения ; время срабатывания привода tк; рабочие ходы Lш штоков исполнительных аппаратов; управляющие функции WT отключающие подачу топлива в двигателе и WC создающие искусственное сопротивление движению автомобиля на режимах частичного торможения. Аналогично выявляют параметры взаимного влияния для каждого блока исследуемого процесса.

4) После установления причинно-следственных связей приступают к составлению математической модели исследуемого процесса. Как правило, такая модель формируется из математических описаний каждого отдельно взятого блока, входящего в состав структурной схемы исследуемого процесса. Рассмотрим пример составления математического описания для блока «колеса с эластичными шинами», входящей в состав схемы. Колеса, оснащенные эластичными шинами, преобразуют действующие на них кинематические и силовые параметры, в продольные Rxij и боковые Ryij реакции, определяя показатели процесса торможения автомобиля.

Входными функциями блока «колеса с эластичными шинами» являются: тормозные моменты МTij, развиваемые тормозными механизмами или тормозом - замедлителем, и их первые производные ; время срабатывания тормозных механизмов tc; углы поворота управляемых колес Q1 и Q2; углы увода эластичных шин dij; продольная Vx и боковая Vy составляющие вектора скорости автомобиля; нормальные реакции Rzij опорной поверхности в пятне контакта шины с дорогой; коэффициенты продольного jx и бокового jy сцепления колеса с опорной поверхностью; высоты микронеровностей дороги q(l).

Выходными функциями блока «колеса с эластичными шинами» являются: продольные Rxij и боковые Ryij реакции; угловые частоты wij вращения колес автомобиля; и моменты Mj сцепления колес с опорной поверхностью; разность тормозных сил DFT колес левой и правой стороны; силы FDш деформации шин, при движении по дорожным микронеровностям.

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

В начале составим расчетную схему (рис. 30) процесса торможения колеса с эластичной шиной и запишем уравнение динамики.

Для описания динамики процесса торможения автомобильного колеса составим уравнение моментов, относительно оси его вращения и решим

Рис. 30. Расчетная схема процесса динамики торможения автомобильного колеса с эластичной шиной

 

относительно старшей производной:

(34)

где wк – угловая скорость вращения колеса [с-1]; Мт – тормозной момент [Н× м]; Мf – момент сопротивления качению [Н× м]; Rx – реализованная касательная реакция (тормозная сила) [Н]; rко – радиус качения колеса в ведомом режиме (силовой радиус) [м]; Jк – момент инерции колеса [кГ× м2].

Тормозной момент в математической модели может быть задан в виде линейной функции:

Мт = Кт × t (35)

где Кт– темп нарастания тормозного момента [Н× м/с]; t – координата времени [с].

Момент сопротивления качению определяется по выражению:

Мf = Rz × f(v) × rко (36)

где Rz– нормальная реакция от действия нагрузки на колесо Gк [Н]; f(v) – коэффициент сопротивления качению, зависящий от скорости автомобиля V.

Если скорость движения автомобиля задана в размерности [км/час], то коэффициент сопротивления качению f(v) рассчитывают по формуле:

(37)

где f0– значение коэффициента сопротивления качению при нулевом значении скорости V.

Для вычисления реализованной касательной реакции (тормозной силы) Rx воспользуемся формулой:

Rx = Rz × f(s) × j mах (38)

где j mах – максимальное значение коэффициента сцепления колеса с опорной поверхностью дороги в области критического проскальзывания колеса; f(s) – функция проскальзывания (нормированная f(S) - диаграмма).

Математическое описание нормированной f(s) – диаграммы представлено в виде функции:

(39)

где a1 и b1 коэффициенты, которые определяют вид функции f(s); S – коэффициент проскальзывания пятна контакта шины относительно опорной поверхности дороги:

(40)

Для того чтобы вычислить коэффициенты a1 и b1 функции f(s) используют экспериментальные j (s) – диаграммы автомобильных шин, внешний вид которых представлен на рис. 31, а).

Для этого сначала на каждой j (s) – диаграмме определяют максимальное

Рис. 31 Экспериментальные характеристики автомобильной шины а) - j (s) – диаграммы; б) нормированные  
 

значение коэффициента сцепления j mах колеса с опорной поверхностью дороги в области критического проскальзывания. Затем каждую точку j (s) – диаграммы делят на j mах. Частные от деления являются ординатами графика нормированной функции проскальзывания f(s), представленной на рис. 31, б).

На графиках (рис. 31, б) нормированной функции проскальзывания f(s) определяют два параметра: коэффициент «жесткости» проскальзывания h s и коэффициент fd снижения фрикционных свойств шины в блоке.

Численное значение коэффициента «жесткости» проскальзывания h s определяется из графика (рис. 31, б) в начальной области устойчивого торможения колеса при S®0 как отношение:

(41)

Численное значение коэффициента fd снижения фрикционных свойств шины в блокетоже измеряется на графике (рис. 29, б) при S=1.

Затем рассчитывается коэффициент Кz используя уравнение:

; (42)

Далее определяем коэффициента b11 по формуле:

(43)

Находим коэффициент b1 из уравнения:

; (44)

И в заключении, рассчитываем коэффициент a1 из выражения:

(45)

Таким образом, уравнения 20¸ 31 составляют математическое описание одного из блоков исследуемого процесса, схема которого приведена на рис. 29. Аналогично разрабатывают математические описания остальных шести блоков, входящих в состав структурной схемы (рис. 29).

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

Часто в математических моделях приходится решать дифференциальные уравнения (например, уравнение 20). При этом используют численные методы: метод Эйлера, метод Рунге-Кута и пр. Поэтому программу пишут на одном из алгоритмических языков, позволяющих решать такие задачи математической модели. Такими языками являются: MATCHKAD, BASIC, TURBO-BASIC, PASKAL, DELFI, C, C+ и т.п.

а) В начале алгоритма расчета прописывают процедуру ввода исходных данных. Для нашего примера необходимо ввести в программу такие исходные данные, как: V – начальная скорость автомобиля [м/с]; Jк – момент инерции колеса [кГ× м2]; Мт max – максимальное значение тормозного момента; rко – радиус качения колеса в ведомом режиме (силовой радиус) [м]; Кт– темп нарастания тормозного момента [Н× м/с]; f0– значение начального коэффициента сопротивления качению; j mах – максимальное значение коэффициента сцепления; значение коэффициента «жесткости» проскальзывания h s; значение коэффициента fd снижения фрикционных свойств шины в блоке; шаг интегрирования дифференциального уравнения Δ t [с]; время расчета исследуемого процесса – tп [с].

б) Затем выполняют расчет исходных параметров модели, таких как:

начальное значение угловой скорости колеса: wк = ; численные значения коэффициентов для расчета нормированной функции проскальзывания Кz - (по формуле 28); b11- (по формуле 29); b1 - (по формуле 30); a1 - (по формуле 31).

в) Затем последовательно решают основные уравнения модели с 26 по 20. В программу расчета вводят ограничения. Например: Если Мт ³ Мт max, то переменной Мт присваивать значение Мт max.

г) Для интегрирования дифференциального уравнения (20) численным методом Эйлера, в программу расчета вводят оператор интегрирования:

wк (i) = wк (i-1) × Δ t (46)

где wк (i) – значение угловой скорости колеса на i-м шаге интегрирования;

wк (i-1) – значение угловой скорости колеса на предыдущем, (i-1)-м шаге интегрирования; - угловое ускорение колеса, определяемое по формуле (20); Δ t – шаг интегрирования (для обеспечение приемлемой точности расчетов Δ t =0, 001с. Не более! ).

В программе расчета угловой скорости колеса обязательно вводят ограничение: Если wк £ 0, то переменной wк присваивать значение 0.

Для расчета времени исследуемого процесса ti вводят оператор:

ti = t(i-1) + Δ t (47)

где t(i-1) – значение времени на предыдущем шаге интегрирования.

д) Как правило, результаты расчета на модели необходимо визуализировать либо в цифровом, либо в графическом виде. Поэтому в программу записывают операторы вывода на экран монитора либо цифровой информации, либо графиков функций, например: wк = f(t); Мт = f(t); Rx = f(t), или Rx = f(s).

е) Решение основных уравнений модели (26 ¸ 20 и 32¸ 33) рекомендуется решать в цикле. Цикл последовательного решения уравнений модели можно задать либо специальными операторами (предусмотренными в языках программирования) либо написанием в программе условного оператора, например: Если время счета t £ tп то следовать к первому (по порядку записи в программе) уравнению модели, в противном случае – прекратить расчеты.

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

6) Для проверки адекватности математической модели выполняют предварительные расчеты исследуемого процесса. Это необходимо для того, чтобы количественно оценить погрешности расчетов с результатами эксперимента.

7) Проверку адекватности математической модели выполняют на основе статистических методов.

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

, (48)

где: а и b - параметры модели, Е - отклонение от линии регрессии.

На первом этапе с использованием пакета MIKROSOFT EXCEL рассчитывается значение E по формуле:

, (49)

где: y- данные, полученные на модели; yie - экспериментальные данные;

n - степень свободы.

Критерий значимости корреляционного коэффициента рассчитывается с использованием выражения:

, (50)

где: r -корреляционный коэффициент.

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

, (51)

Ковариацией двух случайных величин называется математическое ожидание произведения отклонений X и Y от своих математических ожиданий:

, (52)

где X и Y - случайные параметры; , - математическое ожидание параметров X и Y.

Математически, коэффициент корреляции может быть представлен как:

, (53)

, (54)

где: X и Y - значения, полученные в ходе эксперимента и на модели;

D(X) и D(Y) - дисперсии значений полученных в ходе эксперимента и на модели X и Y;

cov[X, Y] - ковариациязначений полученных в ходе эксперимента и на модели X и Y.

Для расчета коэффициента корреляции наиболее удобна следующая формула:

. (55)

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

Для оценки адекватности модели также используют F-критерий. Величина F определяется с использованием программы MIKROSOFT EXCEL, в результате проведения дисперсионного анализа. На основании сравнения расчетного значения критерия Фишера с его критическим значением, делается заключение о значимости модели и ее адекватности результатам эксперимента.

8) Если по результатам расчетов математическая модель исследуемого процесса признаются не адекватной, то производят настройку математической модели с целью её уточнения или (и) дополнения.

Так, например, в вышеописанной математической модели процесса торможения автомобильного колеса уточнению подлежат такие параметры как: Jк – момент инерции колеса [кГ× м2]; Кт– темп нарастания тормозного момента [Н× м/с]; rко – радиус качения колеса в ведомом режиме (силовой радиус) [м]; f(v) – коэффициент сопротивления качению; j mах – максимальное значение коэффициента сцепления, а также значения коэффициентов h s и fd определяющих вид функции проскальзывания f(s).

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

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

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

11) На основании выявленных зависимостей делают выводы о тенденции изменения, например, параметров X объекта исследования при изменении параметров U. Выводы делают о каждой выявленной закономерности.

12) По итогам аналитического исследования составляют соответствующий документ – отчет, заключение, диссертацию, в котором (в зависимости от вида документа) описывают цель, задачи исследования, используемый математический аппарат, методики, а также результаты исследования и выводы.


Поделиться:



Популярное:

Последнее изменение этой страницы: 2016-05-29; Просмотров: 1102; Нарушение авторского права страницы


lektsia.com 2007 - 2024 год. Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав! (0.045 с.)
Главная | Случайная страница | Обратная связь