
2065
.pdfПо остальным параметрам а2, а3 … аm частные производные имеют аналогичный вид.
n |
yi |
f xi, a1, a2... am fa1 xi, a1, a2... am |
0, |
|
|||||||||||||||
|
|
||||||||||||||||||
i 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
. . . . . . . . . . . . . . . . . . . |
(3.8) |
||||||||||||||||||
n |
y |
|
f x |
|
, a , a |
|
... a |
|
f |
x |
, a , a |
|
... a |
|
|
0. |
|
||
|
|
|
|
|
|
|
|
||||||||||||
i 1 |
|
i |
|
i |
1 |
2 |
|
m |
am |
i |
1 |
|
2 |
|
|
m |
|
|
|
|
Решение этой системы относительно а1, а2 … аm дало искомые |
||||||||||||||||||
наилучшие значения числовых параметров. |
|
|
|
|
|
|
|||||||||||||
|
Регрессионные зависимости оцениваются мерой достоверности |
||||||||||||||||||
R2, которая находится в пределах [53]: |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
0 ≤ R2 ≤ |
1. |
|
|
|
|
|
|
(3.9) |
||
|
При R2 = 0 величины, для которых определяются уравнения |
||||||||||||||||||
регрессии, являются |
независимыми; |
при |
R2 |
= |
1 имеет |
место |
функциональная (а не статистическая) зависимость. Принято считать допустимым R2 ≥ 0,7 [53].
Программный продукт SigmaPlot позволяет находить уравнение регрессии для построенной поверхности y=f(x,z).
Аппроксимируемые зависимости соответствуют рабочему процессу полутяжелого автогрейдера ГС 18.07. Значения варьируемых геометрических параметров берутся как в сторону увеличения, так и уменьшения номинальных (табл. 3.2).
|
|
Таблица 3.2. |
Значения варьируемых геометрических параметров |
||
|
|
|
Номинальное значение параметра |
|
Интервал варьирования |
L 6,2м |
|
5,2 L 7,2м |
LБ 1,5м |
|
1 LБ 2м |
K 0,4 |
|
0,3 K 0,7 |
Полученные зависимости |
Ky f (L, LБ , K) представлены в виде |
совокупности поверхностей (рис. 3.28). На рис. 3.29, 3.30 представлены сечения этих поверхностей, позволяющие более подробно увидеть, как отдельные параметры влияют на критерий эффективности.
90

|
|
Ky |
|
|
|
|
|
|
|
|
L=7,2м |
6.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L=6,7м |
6.0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L=6,2м |
5.5 |
|
|
|
|
|
|
|
|
|
L=5,7м |
5.0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
L=5,2м |
4.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4.0 |
|
|
|
|
|
|
|
|
|
|
3.5 |
|
|
|
|
|
|
|
|
|
|
3.0 |
|
|
|
|
|
|
|
|
|
|
2.50.650.60 |
|
|
|
|
|
|
|
2.0 |
|
|
|
0.550.500.450.40 |
|
|
|
|
|
1.6 |
1.8 |
|
|
|
|
|
|
|
1.4 |
|
|||
|
|
K |
|
0.35 |
|
|
1.2 |
|
,м |
|
|
|
|
|
|
0.30 |
1.0 |
|
Lб |
|
|
Рис. 3.28. Графики зависимости коэффициента сглаживания в продольной |
||||||||||
плоскости от длины базы, длины базы балансира и коэффициента базы |
||||||||||
|
|
|
автогрейдера |
|
|
|
|
|||
Ky |
|
|
|
|
|
|
|
|
|
|
6.5 |
|
|
|
|
|
|
|
|
|
|
L=7,2м |
|
|
|
|
|
|
L=6,2м |
L=6,7м |
||
6 |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
5.5 |
|
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
4.5 |
|
|
L=5,2м |
|
|
|
|
|
||
L=5,7м |
|
|
|
|
|
|
||||
4 |
|
|
|
|
|
|
|
|
|
Lб, м |
1 |
|
1.25 |
|
|
1.5 |
|
|
|
1.75 |
2 |
Рис. 3.29. Графики зависимости коэффициента сглаживания в продольной |
||||||||||
плоскости от межосевого расстояния колес балансирной тележки и длинны базы |
||||||||||
|
|
|
автогрейдера при K=0,6 |
|
|
91

Ky |
|
|
|
|
6.5 |
|
L=6,7м |
|
|
6 |
L=7,2м |
|
|
|
5.5 |
|
|
|
|
|
|
|
|
|
5 |
|
|
|
|
4.5 |
|
|
|
L=5,2м |
4 |
|
|
|
|
|
|
|
L=5,7м |
|
3.5 |
|
|
|
|
3 |
|
|
|
L=6,2м |
|
|
|
|
|
2.5 |
|
|
|
K |
0.3 |
0.4 |
0.5 |
0.6 |
0.7 |
Рис. 3.30. Графики зависимости коэффициента сглаживания в продольной |
||||
плоскости от длины базы и коэффициента базы автогрейдера при LБ=1,5м |
Полученные зависимости K f (L, LБ , K) представлены в виде совокупности поверхностей (рис. 3.31).
|
|
|
L=7,2м |
K |
|
20 |
||
|
||
L=6,7м |
|
|
|
18 |
|
L=6,2м |
|
|
L=5,7м |
16 |
|
|
||
L=5,2м |
14 |
|
|
||
|
12 |
|
|
10 |
|
|
8 |
60.65 |
|
0.600.55 |
|
0.50 |
|
K |
0.450.400.350.30 1.0 |
2.0
1.8
|
1.4 |
1.6 |
|
1.2 |
,м |
||
|
|||
|
б |
||
|
L |
|
Рис. 3.31. Графики зависимости коэффициента сглаживания в поперечной плоскости от длины базы, длины базы балансира и коэффициента базы автогрейдера
92

На рис. 3.32, 3.33 представлены сечения этих поверхностей, |
|||||
позволяющие более подробно увидеть, как отдельные параметры |
|||||
влияют на критерий эффективности. |
|
|
|||
Kγ |
|
|
|
|
|
15 |
|
|
|
|
|
L=7,2 |
L=6,7 |
L=6,2м |
|
L=5,7 |
|
14 |
|
|
|
|
|
13 |
|
|
|
|
|
|
|
L=5,2 |
|
|
|
12 |
|
|
|
|
|
11 |
|
|
|
Lб, м |
|
1 |
1.25 |
1.5 |
1.75 |
||
2 |
|||||
Рис. 3.32. Графики зависимости коэффициента сглаживания в поперечной плоскости |
|||||
от межосевого расстояния колес балансирной тележки и длинны базы автогрейдера |
|||||
|
|
при K=0,5 |
|
|
Kγ |
|
|
|
|
19 |
|
|
|
|
18 |
|
L=6,7 |
|
|
17 |
|
|
|
|
16 |
L=7,2 |
|
|
|
15 |
|
|
|
|
14 |
|
|
|
|
13 |
|
|
|
|
12 |
|
|
|
L=5,2 |
11 |
|
|
|
|
|
|
|
м |
|
10 |
|
|
L=6,2 |
|
9 |
|
|
L=5,7м |
|
8 |
|
|
м |
|
7 |
|
|
|
K |
0.3 |
0.4 |
0.5 |
0.6 |
0.7 |
Рис. 3.33. Графики зависимости коэффициента сглаживания в поперечной |
||||
плоскости от длины базы и коэффициента базы автогрейдера при LБ=1,25 м |
Анализируя полученные графические зависимости, можно сделать вывод о том, что изменение варьируемых геометрических параметров автогрейдера в указанных интервалах оказывает различное влияние на коэффициенты сглаживания обрабатываемой поверхности.
93
Наиболее существенное влияние оказывают значения коэффициента базы машины и длина базы. С увеличением значений коэффициента базы машины и длина базы автогрейдера улучшается его планирующая способность. Межосевое расстояние колес балансирной тележки не значительно влияет на планирующую способность в обеих плоскостях.
Результаты аппроксимации зависимостей критериев эффективности от варьируемых параметров приведены в табл. 3.3.
Таблица 3.3
Уравнения регрессии, аппроксимирующие зависимости критериев эффективности от основных геометрических параметров автогрейдера
Критерий |
Длина |
|
Уравнение регрессии |
|
Досто- |
эффек- |
базы |
|
|
верность |
|
тивности |
L |
|
|
|
R2 |
|
5,2 |
Ky=1.332+0.353*Lб+3.376*K-0.113*Lб2+2.586*K2 |
0.9910 |
||
|
|
|
|
||
|
5,7 |
Ky=1.337+0.623*Lб+2.33*K-0.21*Lб2+4.554*K2 |
0.9944 |
||
|
|
|
|
||
Ky |
6,2 |
Ky=0.036+1.149*Lб+6.538*K-0.371*Lб2+0.918*K2 |
0.9937 |
||
|
|
|
|
||
|
6,7 |
Ky=-1.208+1.363*Lб+11.316*K-0.423*Lб2-3.469*K2 |
0.9882 |
||
|
|
|
|
||
|
7,2 |
Ky=-2.651+1.670*Lб+16.985*K-0.523*Lб2-8.847*K2 |
0.9826 |
||
|
|
|
|
||
|
5,2 |
K =2.287+0.183*Lб+17.997*K-0.031*Lб2+0.454*K2 |
0.9874 |
||
|
|
|
|
||
|
5,7 |
K =1.590+0.632*Lб+18.475*K-0.142*Lб2+3.138*K2 |
0.9846 |
||
|
|
|
|
|
|
K |
6,2 |
|
2 |
2 |
0.9843 |
|
K =-0.785+1.859*Lб+24.544*K-0.5*Lб -0.341*K |
|
|
||
|
6,7 |
K =-2.673+1.398*Lб+35.705*K-0.411*Lб2-9.914*K2 |
0.9727 |
||
|
|
|
|
|
|
|
7,2 |
|
K =-5.705+2.974*Lб+44.626*K-0.882*Lб2- |
|
0.9584 |
|
|
17.848*K2 |
|
94

3.3.2 Алгоритм решения задачи безусловной оптимизации градиентным методом
Одна из основных целей проектирования заключается в оптимизации решений, т.е. в достижении заданных характеристик при наименьших затратах или наилучших характеристик проектируемых систем при ограниченных затратах имеющихся ресурсов [15].
Сущность оптимизации сводится к отысканию при наложенных ограничениях таких значений переменных х1, х2, …, хn, которые дают минимум (или максимум) целевой функции [15]
f (X) f x1, x2,..., xn . |
(3.10) |
Общая задача оптимизации может быть сформулирована в следующем виде. Необходимо найти значения переменных х1, х2, …, хn, при которых целевая функция f принимает экстремальное значение с учетом функциональных ограничений (равенств) и граничных условий (неравенств).
Задачи оптимизации с точки зрения методов решения делятся на два класса [15, 59, 86]:
-задачи безусловной оптимизации;
-задачи условной оптимизации.
Задача безусловной оптимизации представляет собой поиск оптимума целевой функции без всяких дополнительных условий
[15, 59, 86]:
f(x) → min(max). |
(3.11) |
Такие задачи на практике встречаются крайне редко, но метод их решения служит основой для решения практических задач оптимизации.
Задача условной оптимизации в общем виде записывается как
[15, 59, 86]:
F = f(xj) → min;
gi(xj) ≤ bi;
dj ≤ xj ≤ Dj; (3.12) i = 1…m;
j = 1…n.
95

В систему уравнений (3.12) входят три составляющие:
-целевая функция показывает, в каком смысле решение должно быть оптимальным, то есть наилучшим, при этом возможны три вида назначения целевой функции: максимизация, минимизация, назначение заданного значения.
-ограничения устанавливают зависимости между переменными.
-граничные условия показывают, в каких пределах могут быть значения искомых переменных в оптимальном решении.
Метод множителей Лагранжа применим при наличии функциональных ограничений вида [15, 59]
fj = fj (x1, x2,…, xn) = 0, |
(3.13) |
где j = 1, 2,…, m.
Для целевой функции Z (x1, x2,…, xn) справедливо уравнение
[15, 59]
dZ |
Z |
dx |
Z |
|
|
dx |
|
|
... |
Z |
dx |
|
0 |
|||||
x |
|
|
|
|
|
n |
||||||||||||
|
1 |
|
x |
2 |
|
|
|
|
2 |
|
|
x |
n |
|
|
|||
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
или |
dZ Z dxi |
0. |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|||||||||||||
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
i 1 |
xi |
|
|
|
|
|
|
|
||||||
Продифференцировав равенство (3.14), получим [15, 59] |
||||||||||||||||||
|
|
df1 |
|
n |
|
f |
1 |
dxi |
|
|
|
|
|
|||||
|
|
|
|
|
|
0 |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
i 1 xi |
|
|
|
|
|
|
|
|
|
|||||
|
|
........................... . |
|
|
|
|
||||||||||||
|
|
dfm |
|
n |
|
fm |
|
|
dxi |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
i 1 xi |
|
|
|
|
|
|
|
|
(3.14)
(3.15)
(3.16)
Каждое из полученных m уравнений теперь умножим на пока еще неизвестный параметр λ, называемый множителем Лагранжа [15, 59]:
df f1 |
dx 0 |
|||||
|
n |
|
|
|
|
|
1 1 |
|
1 |
|
i |
|
|
i 1 |
xi |
|||||
|
|
|
|
|||
df f2 |
|
|||||
dx 0 |
||||||
|
n |
|
|
|
|
|
2 |
2 |
i 1 |
2 xi |
i |
|
|
. |
(3.17) |
|||
..................................... |
|
|
||||||||||
|
|
df |
|
n |
|
f |
m |
dx |
|
|
|
|
|
|
|
|
|
|
|
0 |
|
||||
m |
m |
m xi |
i |
|
||||||||
|
|
i 1 |
|
|
|
|
96
Сложив уравнения (3.17) и уравнение (3.15), получим [15, 59]:
n |
|
Z |
|
f1 |
|
f2 |
|
fm |
|
|
|
|
|
|
1 |
2 |
... m |
|
0. |
(3.18) |
|||||
|
|
|
|
|||||||||
|
xi |
xi |
xi |
|
dxi |
|||||||
i 1 |
|
|
|
|
xi |
|
|
Поскольку все параметры xi независимы, то для того, чтобы это уравнение удовлетворялось, достаточно, чтобы каждый из n членов равнялся нулю [15, 59].
Таким образом, получаем n уравнений [14, 59]
Z |
1 |
f1 |
2 |
f2 |
... m |
fm |
0. |
(3.19) |
|
x |
x |
x |
x |
||||||
|
|
|
|
|
|||||
i |
|
i |
|
i |
|
i |
|
|
Кроме того, имеется еще m уравнений (3.13), определяющих ограничения [15, 59].
Решение системы m + n уравнений и дает искомое оптимальное решение [15, 59].
Таким образом, задача оптимизации стала безусловной и свелась к нахождению экстремума целевой функции.
Задача безусловной оптимизации решена методом первого порядка – градиентным. Данный метод обеспечивает необходимую точность.
Алгоритм метода заключается в следующем (рис. 3.34) [86]:
1. Задать x0, 0<ε<1, ε1 > 0, ε2 > 0, M – предельное число итераций.
Найти |
|
градиент |
|
функции |
в |
произвольной |
||
точке f(x) |
|
f (x) |
|
f (x) |
T |
|
|
|
|
,..., |
|
|
|
||||
x |
x |
|
|
|
||||
|
n |
. |
|
|
||||
|
|
1 |
|
|
|
|
|
2.Положить k = 0.
3.Вычислить f(xk).
4.Проверить выполнение критерия окончания ║ f(xk)║≤ ε1: а) если критерий выполнен, то расчет окончен и x* = xk;
б) в противном случае перейти к пункту 5.
5.Проверить выполнение неравенства k ≥ M:
а) если неравенство выполнено, то расчет окончен и x* = xk; б) если нет, перейти к пункту 6.
6.Задать величину шага tk.
7.Вычислить xk+1= xk- tk f(xk).
8.Проверить выполнение условия f(xk+1) - f(xk) < 0:
а) если условие выполнено, то перейти к пункту 9;
97

б) если нет, то вернуться к пункту 7, положив tk tk .
2
9. Проверить выполнение условий ║xk+1- xk ║<ε2 , │f( xk+1)-f( xk)│<
ε2 :
а) если оба условия выполнены при текущем значении k и k=k-1, то расчет окончен, x* = xk+1;
б) если хотя бы одно из условий не выполнено, положить k=k+1 и перейти к пункту 3.
A1 |
|
|
|
Начало |
|
|
|
B1 |
B2 |
|
2 |
Ввод |
Значения |
|
|
исходных |
x0 , M,0< <1, |
H1 |
|
данных |
1 |
>0, |
|
|
>0, 2 |
|
|
C1 |
|
|
Циклдля tk |
|
|
|
|
Нахождение |
|
|
I1 |
f(x) |
|
|
D1 |
|
|
|
|
|
|
|
Расчетx k+1 |
|
||
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
Цикл для k=0 |
|
|
|
3 |
|
|
|
||||
|
|
|
|
|
|||||||
c шагом1 |
|
|
|
|
|
||||||
|
|
|
J1 |
Нет |
|||||||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
f(x k+)1-f(xk)<0? |
|
|
E1 |
|
|
|
|
|
|
k |
|
Да |
|
|
Расчет f(x ) |
|
K1 |
|
|
|
|
|
|
Нет |
|
|
|
|
|
|f(xk+)1-f(xk)|< 2 |
|
Да |
F1 |
|
F2 |
||xk+-1 xk||< 2 ? |
|
k |
1 |
Да |
|
||
|
|| f(x )|| 1? |
1 |
L1 |
|
|
|
|
|
|
|
|
|
|
Нет |
G2 |
x*=xk+1 |
|
Да |
G1 |
|
|
|
|
k M? |
M |
|
|
||
|
|
M1 |
|
||
|
|
|
|
|
|
|
|
Нет |
|
Конец |
|
|
|
|
|
|
|
1 |
|
2 |
|
|
|
H2
tk tk=2
K2
k=k+1 3
Рис. 3.34. Алгоритм решения задачи безусловной оптимизации градиентным методом
98
3.3.3.Выбор оптимальных значений геометрических параметров автогрейдера
Выбор оптимальных значений параметров производился градиентным методом.
В результате подстановки целевых функций для каждого значения длины базы автогрейдера были получены оптимальные значения коэффициента базы машины и межосевого расстояния колес балансирной тележки, которые представлены в табл. 3.4.
Таблица 3.4
Cравнение результатов оптимизации и действительных геометрических параметров автогрейдера
Исходный |
|
Найденные рациональные и |
||
параметр |
Автогрейдер |
действительные значения |
||
|
|
|
|
|
L,м |
|
K |
LБ ,м |
|
5,2 |
New Holland F106.6 |
0,38 |
1,25 |
|
Теоретический |
0,67 |
1,09 |
||
|
||||
5,7 |
Komatsu GD510R-1 |
0,44 |
1,53 |
|
Теоретический |
0,65 |
1,06 |
||
|
||||
6,2 |
ГС 18.07 |
0,40 |
1,5 |
|
Теоретический |
0,63 |
1,43 |
||
|
||||
6,7 |
Volvo G990 |
0,43 |
1,59 |
|
Теоретический |
0,68 |
1,04 |
||
|
||||
7,2 |
Komatsu GD825A-2 |
0,44 |
1,84 |
|
Теоретический |
0,61 |
1,37 |
||
|
Исходя из данных, приведенных в табл. 3.4, можно сделать вывод о неоптимальности значений рассмотренных геометрических параметров у существующих автогрейдеров.
3.4. Тягово-сцепной расчет
Как известно, тягово-сцепной расчет автогрейдеров заключается в определении основных параметров, в первую очередь полностью отвечающих требованиям тягового режима работы при выполнении различных технологических операций [20,21,57,58,78,81,102,104,105].
При возведении земляного дорожного полотна на один проход резания требуется 2—4 прохода по разравниванию и планировке грунта. Поэтому резание грунта автогрейдером необходимо
99