
2065
.pdf
A2
|
|
|
Начало |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
B2 |
|
|
|
|
|
|
|
|
|
B3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||||
|
Ввод |
|
|
|
Значения |
|||||||||
|
исходных |
|
|
|
|
|
параметров |
|
||||||
|
|
|
|
|
|
|||||||||
|
данных |
|
|
|
, , ,h,n |
C2 |
Расчет параметров a0,a1,q1,q2
|
|
D2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
Цикл для i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
с шагомh |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
E2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
Расчет |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
YМЛ,YМП |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
F2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
F3 |
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
Расчет |
|
|
|
|
|
|
Размер |
|
|
||||||||||||||
|
|
|
|
пятна |
|
|
|
||||||||||||||||||
|
|
Y1Л,Y1П |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
контакта 2X0 |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
G2
Запись файла данных
H1 |
|
H3 |
|
|
H2 |
Да |
|
Значение n |
i=i+h |
||
i n |
Нет
I2
Конец
Коэффициенты, характеризующие корреляционнуюфункцию случайного процесса
Коэффициенты, характеризующие рекуррентное уравнение случайного процесса
Вертикальные координаты соответственно левой и правой колеи несглаженного профиля
Вертикальные координаты соответственно левой и правой колеи сглаженного профиля
Рис. 2.16. Блок-схема алгоритма реализации математической модели микрорельефа
60

|
|
Y |
|
|
Определение FYЛ |
|
|
|
|
|
|
|
||
Определение FФY |
1 |
FФ |
|
|
|
|
|
KД |
|
|
|
|
||
|
|
|
|
Определение YR |
|
|
YРОЛЗ |
Л |
|
|
|
|
||
|
|
|
|
|
|
YРОЗ |
e |
-tгp |
1 |
КШ |
||||
|
|
|
|
|
FX |
|
|
|
|
|||||
Определение FФX |
|
|
|
|
|
|
|
|
|
Трр+1 |
р |
|||
|
sin2 |
|
|
g |
3 |
|
e |
-tгp |
1 |
КШ |
||||
|
|
|
|
|
|
|
|
|
|
Трр+1 |
р |
|||
|
|
|
|
|
|
|
|
РОЗ |
L|2 |
П |
|
|||
|
|
|
|
|
|
|
|
YРОПЗ |
|
|
|
|
||
|
|
|
|
|
sin |
Определение FПР |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||
Микрорельеф |
|
X |
|
|
Определение FXЛ |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
e-tПp |
KД |
|
|
|
|
|
yП(n) |
Y1П |
|
X+X0 |
|
|
|
|
DYOZ |
|
|
|
|
|
hШтП |
1 |
y(X)dX |
WK |
|
Y1 |
|
K2 |
hШтП |
|
|
|
|
|
||
|
|
2X0 |
|
|
|
|
|
|
|
|
||||
|
|
X-X0 |
X+X0 |
|
1/ 2 |
1-K |
|
e-tp |
|
|
РО |
|
|
|
yЛ(n) |
|
1 |
WK |
|
YРОП |
|
tgg |
|
|
|||||
|
y(X)dX |
|
|
|
|
|
1/ L3 |
|
|
|
||||
|
|
|
|
|
YR YРО |
|
|
|
|
|
||||
|
Y1Л |
2X0 |
|
|
|
YРОЛ |
|
|
|
|
||||
Y2П |
X-X0 |
|
|
|
|
|
|
|
|
|
|
|||
|
1 |
X+X0 |
WK |
1-KББ |
|
|
|
hШтЛ |
|
YРО |
|
|
|
|
|
y(X)dX |
|
|
|
|
|
|
|
||||||
|
|
2X0 |
|
|
|
|
|
|
hШтЛ |
|||||
|
Y3П |
X-X0 |
X+X0 |
|
|
YБП |
|
|
|
|
|
|
|
|
e |
-tБp |
1 |
y(X)dX |
WK |
KББ |
|
|
YБ |
|
|
|
|
|
|
|
2X0 |
|
1/ 2 |
|
|
|
|
|
|
|||||
|
Y3Л |
X-X0 |
|
|
|
YБЛ |
K |
|
|
|
|
|
|
|
|
1 |
X+X0 |
|
KББ |
|
|
|
|
|
|
||||
e |
-tБp |
y(X)dX |
K |
|
|
|
|
|
|
|
|
|
||
|
2X0 |
|
W |
|
|
|
DYOX |
|
|
|
|
|
|
|
|
|
X-X0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X+X0 |
|
|
|
|
|
|
|
|
|
|
|
Y2Л |
|
1 |
y(X)dX |
WK |
1-KББ |
|
|
K1 |
|
|
|
|
|
|
|
2X0 |
|
|
|
|
|
|
|
|
|
||||
|
X-X0 |
|
|
|
|
|
e-tЛp |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рис. 2.17. Структурная схема обобщенной математической модели автогрейдера |
|
|
|
|
|
61

A2
Начало
B2
Ввод исходных данных
C2
Расчет постоянных параметров
D2
Цикл для t=0 с шагом t
E2
Мат.модель микрорельефа
F2
Мат.модель грунта
B3
Значения параметров
подсистем
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
hШтП,hШтЛ |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
G2 |
|
|
|
|
|
|
|
|
|
G3 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
Мат.модель |
|
|
|
|
|
|
Мат.модель |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
автогрейдера |
|
|
|
|
|
|
СУРО |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
H2 |
|
YРОП,YРОЛ |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
Запись файла |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
данных:YРО, РО |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
I1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Нет |
|
|
|
I3 |
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Время |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
I2Конец |
|
|
|
|
t=t+ t |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
моделирования |
|
|
|
|
пути? |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
J2 |
|
|
Да |
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
Конец |
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Параметрымашины, обрабатываемого грунта, микрорельефа,СУРО
Коэффициенты, запаздывания
Координатырежущей кромки отвала
Стабилизация положения РО
Рис. 2.18. Блок-схема алгоритма реализации рабочего процесса автогрейдера
62
2.6. Математическая модель рабочего процесса автогрейдера
Для реализации поставленной в работе цели необходимо построить обобщенную математическую модель автогрейдера. Как показал анализ предшествующих исследований, а также анализ рабочего процесса автогрейдера, проведенные в главе 1, обобщенная математическая модель автогрейдера включает в себя математические модели следующих подсистем: базовая машина, СУ РО в поперечной плоскости и продольной плоскостях, микрорельеф обрабатываемой поверхности и учет влияния реакции грунта на вертикальные перемещения автогрейдера. Эти сложные подсистемы описаны рядом уравнений, на основе которых построены обобщенные структурные схемы. Обозначения на рисунках соответствуют обозначениям, используемым ранее при описании каждой подсистемы в отдельности.
На основе обобщенной математической модели автогрейдера была составлена структурная схема, представленная на рис. 2.17.
На рис. 2.18 представлен блок-схема алгоритма, описывающего рабочий процесс автогрейдера. Данный алгоритм используется как составная часть САПР структуры автогрейдера.
2.7. Методика экспериментальных исследований
Комплексный метод исследований предполагает проведение экспериментальных исследований, основными задачами которых являются: подтверждение адекватности математической модели объекта исследования; определение численных значений параметров, входящих в математические модели объекта; подтверждение работоспособности и эффективности технического решения, внедренного в производство [86,80].
При проведении экспериментальных исследований возможно использование двух методов: активного и пассивного [1].
Пассивный эксперимент предусматривает наблюдение за выходными координатами без вмешательства в процесс функционирования системы [1]. Результаты пассивного эксперимента используются для подтверждения адекватности математических моделей, проверки работоспособности конструкторских разработок, определения эффективности предложенных инженерных решений.
63

Активный эксперимент предусматривает формирование детерминированных воздействий на подсистемы и сложную динамическую систему в целом. Выходные координаты и переходные процессы дает информацию о свойствах объекта.
Вданной работе целесообразно использовать методологию как пассивного, так и активного эксперимента, в силу того, что внешние возмущающие воздействия имеют стохастическую природу, а при подтверждении адекватности моделей целесообразно формировать детерминированные воздействия.
Вкачестве оценки достоверности экспериментальных данных использовалась доверительная вероятность – Px. В данной работе
достаточной считалась доверительная вероятность Px 0,95. Исключение грубых ошибок измерений проводилось в
соответствии с ГОСТ 11.002-73. Для этого по данным упорядоченной выборки х1 х2 … хn вычислялось среднее арифметическое выборки [50]:
|
|
|
|
|
1 n |
|
|
(2.47) |
|
|
mx |
|
|
|
xi |
|
|
||
|
|
|
|
|
|||||
и дисперсия: |
|
|
|
|
n i 1 |
|
|
|
|
|
|
1 |
|
n |
|
|
|
||
|
|
|
|
|
|
|
|||
Dx |
|
|
|
xi |
mx |
2, |
(2.48) |
||
|
|
|
|||||||
|
|
n 1i 1 |
|
|
|
где хi – значение i-го наблюдения; n – число наблюдений.
Грубой считается ошибка, при которой выполняется неравенство:
xi mx |
tk, , |
(2.49) |
|
x |
|||
|
|
где tk, – табличное значение параметра распределения Стьюдента при степени свободы k = n - 1 и доверительной вероятности рд = 1 - у; где
у – уровень значимости ( у = 1 - рд).
С учетом принятого в работе допущения о нормальном законе распределения результатов наблюдений, истинное значение математического ожидания с вероятностью pд лежит в пределах:
mx |
|
t |
x |
mx mx |
t |
|
x |
, |
(2.50) |
||||
|
|
|
|
|
|
|
|
|
|||||
|
n |
|
|
n |
|||||||||
|
|
|
|
|
|
|
|
|
|
где t – величина, характеризующая для нормального закона распределения число средних квадратических отклонений, откладываемых от центра рассеивания в обе стороны для того, чтобы вероятность попадания в полученный участок была равна pд.
64

Для среднеквадратического отклонения можно записать [86,110]:
x |
|
t |
x |
|
x x |
t |
|
x |
|
. |
(2.51) |
|
|
|
|
|
|
|
|
|
|||||
|
2n |
|
2n |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
2.8. Оценка адекватности математической модели автогрейдера
Адекватность модели, как способность отображать с необходимой точностью характеристики исследуемого объекта при изменении его параметров и внешних воздействий, является одним из подтверждений работоспособности систем автоматизированного проектирования, правомерности использования предложенных методик проектирования и выводов, полученных при исследовании математической модели. Адекватность подтверждается сравнением результатов, полученных расчетным путем с помощью математической модели, с экспериментальными данными, расхождение между которыми для решения поставленных в работе задач не должно превышать 15% .
В качестве объекта исследования был выбран автогрейдер ДЗ-143. Параметры переходных процессов являются основными показателями, характеризующими динамические свойства автогрейдера [68]. Переходные процессы регистрировались при сбрасывании колес автогрейдера, предварительно поднятых с
помощью домкрата на 0,05 м.
Адекватность статической модели автогрейдера, реализуемой на ЭВМ в соответствии с алгоритмом, блок-схема которого представлена на рис. 2.4, оценивалась по вертикальному перемещению правой крайней точки режущей кромки отвала Yроп при подаче детерминированного воздействия под элементы ходового оборудования.
Задача оценки адекватности в статическом и динамическом режимах решалась путем регистрации перемещений основных элементов базовой машины относительно поверхности обрабатываемого грунта. Перемещения элементов с помощью потенциометрических датчиков преобразовывались в изменения электрического напряжения, которые фиксировались на фотоленте магнитоэлектрическим осциллографом К12-22. Схема размещения датчиков перемещения на автогрейдере приведена на рис. 2.19.
65

D5
D3
D1 |
D6 |
|
D4
D2
Рис. 2.19. Схема размещения датчиков перемещений на автогрейдере
В качестве датчиков использовались проволочные потенциометры I класса точности типа ПЛП кругового вращения, питаемые стабилизированным постоянным напряжением. При колебаниях автогрейдера ось потенциометра поворачивается через шкив нитью, один конец которой закреплен на неподвижной опоре, а второй через пружину - на автогрейдере.
Непосредственное измерение угловых перемещений отвала технически трудно осуществимо. Поэтому угловые перемещения вычислялись по разности вертикальных перемещений правой и левой сторон отвала.
На рис. 2.20 в качестве примера приведен график вертикального перемещения правой крайней точки режущей кромки отвала Yроп при подъеме правого переднего колеса балансирной тележки автогрейдера. При расчете геометрические параметры автогрейдера принимались соответствующими автогрейдеру ДЗ-143. Угол захвата отвала =45o. Высота подъема изменялась от 0 до 0,1 м.
Анализ расчетной зависимости показал, что она достаточно близко совпадает с экспериментальной, полученной в результате натурных испытаний. Относительная ошибка соответствия математической модели ее оригиналу не превышает 10% .
66

YРОП,м
0.1
0.08
Теоретическая
0.06
0.04
Экспериментальная
0.02
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Y2П,м |
|
0.01 |
0.02 |
0.03 |
0.04 |
0.05 |
0.06 |
0.07 |
0.08 |
0.09 |
|
||||||||||||||||||||
0 |
0.1 |
Рис. 2.20. Вертикальное перемещение крайней правой точки режущей кромки отвала при подъеме правого переднего колеса балансирной тележки
Адекватность динамического математического описания автогрейдера была оценена сравнительным анализом поперечных угловых колебаний отвала.
Угол перекоса отвала является важной координатой, определяющей поперечный уклон дорожного полотна. На рис. 2.21 приведен фрагмент переходного процесса колебаний отвала автогрейдера в поперечной плоскости: расчетный переходный процесс, полученный при математическом моделировании с исходными данными, соответствующими автогрейдеру ДЗ-143 при=45o и переходный процесс, полученный экспериментальным путем при единичном ступенчатом воздействии высотой 0,05м под задним колесом балансирной тележки, что в рабочих условиях соответствовало бы съезду заднего колеса с неровности 0,05м.
Анализ переходных процессов показал, что они имеют сложную динамическую природу и, как указано в работе [65], могут быть аппроксимированы выражением
|
k T2 |
4 2 |
|
|
|
2 t |
|
2 |
|
|
|
t 0 0 |
|
|
e |
k t |
|
|
|
|
|
||
2 |
|
sin |
T |
arctg k T |
|
, (2.52) |
|||||
|
|
|
|
|
|
|
|
|
|
где 0 - начальное значение угла наклона; T - период колебаний РО в
поперечной плоскости; k - коэффициент затухания колебаний. Сравнительный анализ теоретической и экспериментальной
кривых, полученных при одинаковых исходных данных, показал, что
67

относительная погрешность по амплитуде колебаний <12% , по периоду колебаний <13% , по коэффициенту затухания <10% .
Достигнутая степень адекватности математической модели достаточна для решения поставленных в работе задач.
Степень расхождения зависит от точности параметров расчетной схемы и технического состояния реальной машины [110].
γ,рад |
|
|
|
|
|
|
|
|
|
|
0.01 |
|
|
|
экспериментальная |
|
|
|
|
||
0.008 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
0.006 |
|
|
|
|
|
|
|
|
|
|
0.004 |
|
|
|
|
|
|
|
|
|
|
0.002 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
-0.002 |
|
|
|
|
|
|
|
|
|
|
-0.004 |
|
|
|
|
|
|
|
|
|
|
-0.006 |
|
|
|
|
теоретическая |
|
|
|
|
|
-0.008 |
|
|
|
|
|
|
|
t,c |
||
-0.01 |
|
|
|
|
|
|
|
|
|
|
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
|
Рис. 2.21. Фрагмент переходного процесса колебаний |
|
|
|
||||||
|
|
отвала автогрейдера в поперечной плоскости |
|
|
|
68
3. ТЕОРЕТИЧЕСКИЕ ОСНОВЫ ИНЖЕНЕРНОЙ МЕТОДИКИ ПРОЕКТИРОВАНИЯ АВТОГРЕЙДЕРА
3.1. Методика теоретических исследований математической модели автогрейдера, выполняющего планировочные работы
При изучении сложных динамических систем широко применяется методология системного анализа. В основе системного анализа лежит понятие системы, под которой подразумевается множество компонентов, обладающих определенными свойствами с фиксированными между ними взаимосвязями [71].
В данной работе системный анализ применяется при исследовании автогрейдера. Машина рассматривается как сложная система, состоящая из отдельных взаимосвязанных друг с другом подсистем. Каждая выделенная подсистема имеет свои определенные свойства и законы функционирования, которые описываются математическими моделями. Математические модели подсистем в совокупности образуют сложную математическую модель автогрейдера.
Решение поставленных в работе задач проводилось на основе базовых положений методологии системного анализа, состоящих в следующем [86]:
–целостность системы при ее взаимодействии с внешней средой;
–структурность системы в виде совокупности связей и их причинности;
–иерархичность системы, т.е. каждый автогрейдер, представленный как система, делится на подсистемы, которые, в свою очередь, делятся на другие подсистемы и т.д. до получения неделимых элементов системы.
Решение задач с применением системного анализа проводится в соответствии со следующими этапами [71]:
1) постановка задачи – определяют объект, цели и задачи исследования, а также критерии для изучения и управления объектом.
2) анализ решаемой задачи – очерчиваются границы изучаемой системы и определяется ее структура; объекты и процессы, имеющие отношения к поставленной цели, разбиваются на изучаемую систему
ивнешнюю среду. Затем выделяют отдельные составляющие части системы, устанавливают взаимодействие между ними и внешней
69