Оптимизация технологических процессов механической обработки
..pdfного штучного времени /шгр
F = х - J ^ + O - X ) 4?=£-, |
|
||
|
*шт.ср |
^оп.ср |
|
где /шт.ср. Соп.ср — значения |
интегральных |
средних производитель |
|
ности (штучного времени) и себестоимости, |
определяемых в |
интерва |
|
лах ДЛЯ скорости (Omin, «maxi. ПОДЭЧИ lSm|n, |
Smax] И ГЛубИНЫ |
рвЗЭНИЯ |
|
^шах]. Эта зависимость |
может быть преобразована к виду (3.84). |
Область технических ограничений в рассматриваемом случае представ ляет криволинейный многогранник, заданный системой неравенств:
по стойкости инструмента
|
|
|
|
|
|
(3.96) |
|
по мощности главного привода станка |
|
|
|
||||
v n*+V z/*z ^ |
/V„11 • 102 • 60 |
|
(3.97) |
||||
|
СЛ |
|
|
||||
|
|
|
|
|
|
||
по прочности державки и жесткости режущего инструмента |
|
||||||
v n*sy4 Xz |
min |
4В, К |
. |
500 |
Г |
13.98> |
|
1вр^г"г'Л |
’ |
СгУ.Р |
|||||
по жесткости детали |
|
|
|||||
|
|
|
|
|
|
||
v |
у |
^ 246fe«Djp6 |
|
(3.99) |
|||
vsyvt у iC ---------- - |
|
L3CPA РУ
по прочности механизма подач
v~nx$yxtxx < |
||
|
СхК |
|
по шероховатости поверхности (для |
||
точения) |
|
|
о.оов/г'^г0*8 |
(3.101) |
|
(ффО0,35 |
||
|
||
по максимальной подаче, допустимой ки |
||
нематикой станка, |
|
|
s < s max; |
(3.102) |
по максимальной скорости, допустимой ки нематикой станка,
o < jiD « max/1000. |
(3.103) |
Найдем точку в области допустимых ре шений, в которой целевая функция F (v, s, f) будет минимальная. При общем под ходе необходимо рассмотреть четыре типа точек А, Б, В, Г (рис. 34).
(3.100)
Рис. 34. Графическое пред ставление задачи определения точек в области технических ограничений, в которых целе вая функция F (v, s, t) при нимает минимальное значе ние (при непрерывных зна
чениях v и s).
Точка А находится внутри многогранника решений; Б — на одной из его криволинейных граней; В — на одном из его ребер; Г — в од ной из вершин многогранника. Далее покажем, что в условиях рас сматриваемой задачи отыскание минимума оценочной функции F (v, s; t) может вестись только в точках В и Г
С л у ч а й 1. Допустим, что целевая функция F (v, s, t) достигав feTминимума в одной из внутренних точек многогранника типа А. Тог да в этой точке должны одновременно выполняться равенства
дF/dv = 0, дFids = 0, dF/dt = 0.
Вычислим частные производные функции
В7Л + JL 1—B'lT-, dTIdv -
|
+ *+£-£]■ |
(3.104) |
|
Используя |
известную зависимость |
для стойкости инструмента |
|
|
Т = ClJmf~Xv,mv~1,ms~yv/m9 |
(3.105) |
|
определяем |
значение |
|
|
dTIdv = — — |
= (— Mm) Ш |
(3.106) |
|
|
m |
|
|
Подставив значение dTfdv в выражение (3.104), получим |
|
||
|
-5Г = - ^ [ Л + |
4 - ( 1 - 1//п).- |
(3.107) |
|
|
Анализ полученного выражения показывает, чтоpF/dv — 0 эквивалент но равенству
А + -у -(1 — Mm) = 0. |
( 3 .1 0 8 ) |
Вычисляя аналогично частные производные dF/ds и dF/dt, получаем еше два равенства
< 4 + f - ( l |
____m J = |
о, |
(3 .1 0 9 ) |
||
А |
+ 4 - t l |
____ V ) = |
° . |
(3 .1 1 0 ) |
|
|
|
|
m t |
|
|
Очевидно, что если xv< l, |
yv < |
1 |
(что всегда выполняется в реальных |
||
задачах), то равенства ( 3 .1 0 8 ) , |
(3 .1 0 9 ) не |
могут одновременно иметь |
места. Таким образом, случай 1 (нахождение минимума целевой функ ции F (v, s, f) в точке А) может быть исключен из дальнейшего рассмот рения.
С л у ч а й 2. Допустим, что целевая функция F (a, s, f) достига ет минимума на поверхности одной из криволинейных граней много гранника в точке типа 5 .
Анализ вида технических ограничений (3 .9 6 ) — (3 .1 0 3 ) показыва ет, что все грани рассматриваемого многогранника могут быть описа ны в общем виде параметрическим уравнением поверхности
Л У = /с , |
( 3 . 1 1 1 ) |
где р, q, г, К ~~ коэффициенты и постоянные факторы, определяющие^ вид технического ограничения.
Показатели степеней р, q, г могут в отдельных случаях равняться нулю. При q Ф 0, что выполняется для всех ограничений, кроме огра ничения (3.103), уравнение (3.111) может быть переписано в виде
s = / ( / ; “, |
(3.112). |
где р = — p!q\ а = — r/q. |
F (v, s, |
Подставляя значение s в формулу (3.84) для целевой функции |
t), получаем функцию F, двух переменных и и /, т. е. производим су
жение функций F на поверхности s = |
|
Предварительно выража |
||||
ем стойкость инструмента Т в зависимости (3.105) |
через |
два парамет |
||||
ра — и и I |
|
|
|
|
|
|
|
Т = cl/mr gv,mv~'/m (K1vat*)-Vv"n = К2Г а't r p\ |
(3.113>- |
||||
где К2 = |
К Г Уо,ЩС1/т\ а, = |
Лг+jyjL ; |
р, = |
1+ -fc - . |
|
|
T o r l a |
[А + |
|
" 7 |
W |
м ‘ + В‘Ы |
|
где |
Ах = |
А1Къ Вг = В'1КХК* |
|
(3.114)- |
Для того чтобы минимум функции Z7! достигался в одной из точек, внутри грани s = (т. е. не на ребре), необходимым условием яв ляется выполнение равенств
dFJdv = 0, dFJdt = 0.
Проверим возможность выполнения этого условия. Вычислим эти частные производные, обозначив для краткости 1 + а = а 2, 1 + Р =-
= Рг:
Е± |
^ 2+1 |
[Аг + |
В1/а,УР|] + |
|
1 — |
|
ди |
|
|
|
|
|
|
|
= |
м , |
+ В. (1 - |
х ) |
'“‘Л |
(3-115)- |
Аналогично вычисляя dFjdt и приравнивая обе частные |
производные |
|||||
нулю, записываем систему уравнений |
|
|
|
|||
|
^1 + |
В, (1 — Р1/Р2) f'V*' |
= 0, |
|
||
|
А, + В, (1 — аг/а2) |
|
= 0 . |
(3.116). |
||
Вычитая из |
первого уравнения второе, |
получаем |
|
|||
|
|
|
|
= 0 . |
(3.117) |
Последнее равенство возможно лишь при (V02 = aJ a2-
Таким образом, лишь при случайном совпадении коэффициентов. ах/а2 и Pi/p2 система уравнения может иметь общее решение. Уточним,
что в этом случае уравнения, входящие в (3.116), полностью совпада ют и из них можно найти решение
f V ' = Q , |
(3.118) |
где Q — положительный корень уравнения (3.117). Тогда точка ми нимума оценочной функции должна отыскиваться при решении систе мы уравнений
/ V |
1* |
= Q , |
(3.119) |
s = |
К |
/ Л |
|
Так как система состоит из двух уравнений с тремя неизвестными, ее решение будет представляться в виде некоторой кривой, лежащей на поверхности и обладающей тем свойством, что dF^fdv = dF^ldt = 0 на этой кривой. Это означает, что функция F постоянна на ней и поэтому любую ее точку можно принять за минимум, в частности точку пересе чения кривой с ребром многогранника, т. е. выполнение равенства Р1/Р2 = «i/a2 при случайном совпадении коэффициентов приводит к случаю отыскания минимума целевой функции в точке типа В на ребре многогранника. Поэтому случай 2 (нахождение минимума целевой функции F (и, s, 0 в точке типа Б) может также быть исключен из даль нейшего рассмотрения.
Таким образом, оптимизация режимов механической обработки для непрерывных значений параметров и, s и t сводится к определению минимума целевой функции в точках типа В или Г Первоначально отыскивается минимум целевой функции на ребре многогранника.
Рассмотрим ребро, образованное пересечением двух поверхностей,
описываемых параметрическими уравнениями |
|
||
xf'sq'f' |
= |
Kv |
(3.120) |
vp*sq4rt |
= |
К2- |
(3.121) |
Это кривая, которую можно задать в параметрической форме относи тельно t формулами
|
|
v = |
Qity' , |
|
|
(3.122) |
|
|
s = |
Q2/v% |
|
|
(3.123) |
где |
Yi = |
r#\ -V '» |
. v |
Г2Р1 |
rt^2 . |
|
PiQi — Pi'h ’ Y] |
Q\P* |
Ч0Р1 |
|
|||
|
|
4-< |
|
|
Pi |
|
|
|
Р\Яг~РгЯ\ |
; Q* = |
/(’1 |
\ 4 iP t— QzPi |
|
|
|
I |
|
|
|
|
Подставляя в формулу (3.84) для |
целевой |
функции F (v, s, |
t) выра |
|||
жения переменных и и s и вводя ранее принятые |
обозначения, |
записы |
||||
ваем |
|
|
|
|
|
|
|
Л. (0 = |
[Аг + Я Л |
|
(3.124) |
||
где г|) = Yi + Y2 + |
I’- * = |
|
У2Уу- . |
|
|
|
Найдем минимум целевой функции F2 (f) с помощью производной
~т~ = ~ |
|
|
х |
х |
[фЛ, + |
В,(ф — к) Г]. |
(3.125) |
Приравняв dFJdt = 0, |
получим |
уравнение, из |
которого определим |
значение t |
|
|
|
|
t = [ |
I1'* |
(3.126) |
После определения t по зависимостям (3.122) и (3.123) определяются значения v и s.
Из всех полученных точек для рассмотренных ребер отбрасываются те, которые не принадлежат многограннику ограничений. Последнее положение приводит к необходимости дополнительно определять зна чение функции F (о, s, t) в вершинах многогранника, что упрощает по строение общего алгоритма оптимизации режимов обработки.
Определение точки вершины многогранника (точки типа Г) произ
водится на основе решения системы |
из трех |
уравнений |
tf ls9'f' = |
К1» |
|
№ Ч г'= К г, |
(3.127) |
|
t/V 7 r’ = |
К3. |
|
Для упрощения решения уравнения (3.127) приводятся к линейному виду логарифмированием
Pixi + ЯЛ + ГЛ = bv |
|
||
р2хi + |
ЯЛ + |
г2хз = Ьг, |
(3.128) |
Psxi + |
Язх1 + |
гзхз = Ь3, |
|
где хг = In о; х2 = In s; x3 = |
ln/; 6г = In /Сг; 6a = lnffs; |
63 = 1пЛГ3. |
Система линейных уравнений (3.128) может быть решена одним из известных способов. В рассматриваемом алгоритме эта задача реша лась по правилу Крамера.
Общий алгоритм оптимизации режимов механической обработки для непрерывных значений параметров о, s и t показан на рис. 35.
7. Параметрическая оптимизация технологических процессов при обеспечении
эксплуатационных свойств деталей
При параметрической оптимизации технологических процессов и, в частности, при расчете режимов резания, как было показано выше, преимущественно используются в качестве оценочной функции максимальная производительность, минимальная себестоимость техно логической операции или их компромиссное значение. Однако, как пока зал проведенный анализ, в отдельных случаях целесообразнее исполь зовать другие оценочные функции. Так, при технологическом обеспе-
Г |
Н А Ч А Л О |
3 |
|
|
Ввод данных
Формирование конс тант, входящих в целе - бую функцию и техни ческие ограничения
С о ■ ........
Формирование массива
у (L ) = 0_______
г 4 ---------
1 10ю
с
1L - 0
□« ______
Г™ __
О
ГH-Z.+/
г- 9 ------------ ------------------
Вычисление v(L), s(L ), T(L) -точки минимума оценочной функции на пересечении 1 -го и (1-1)-го ограничений
чении заданных характеристик качества поверхности (микрогеомет рии, волнистости, микротвердости, остаточных напряжений и струк туры поверхностного слоя) с целью повышения эксплуатационных свойств деталей машин более полным критерием оптимальности будет
Рис. 35. Блок-схема алгоритма оптимизации режимов механической обработки для непрерывных значений и, s, L
служить функциональная зависимость эксплуатационного свойства от технологических факторов обработки. Это вызвано тем, что на ха рактеристики качества поверхности оказывает влияние не только ско рость резания и подача, которые определяют критерий оптимальности по технологической себестоимости, но и ряд других факторов (геомет рия и материал инструмента, давление и ток при ЭМО и т. д.).
Для упрощения построений математической модели технологиче скую наследственность целесообразнее учитывать не через одноимен ные характеристики качества поверхности, а через некоторые комплекс ные выражения, в наибольшей степени оказывающие влияние на
рассматриваемый процесс обработки. В настоящей работе для этих це
лей использовался безразмерный комплекс Д = Rm>x/pbt/v, где р, Ь, v — характеристики шероховатости соответственно радиус выступов и параметры кривой опорной поверхности.
Учитывая, что характеристики качества поверхности формируют ся главным образом на окончательной операции, первоначально ограничим рассматриваемую модель рамками одной операции «р, с вы дачей в качестве выходного параметра значения безразмерного ком плекса Д пред. которое может быть использовано в дальнейшем при оп тимизации на уровнях предварительной обработки. Для упрощения описаний математических зависимостей индексы обозначений операций
([) опущены. В этом случае ограничения для рассматриваемой модели можно выразить в виде системы уравнений
Ri = |
|
. . • tnln, |
|
|
R2 = |
k2 An’pV N *” .. |
(3.129) |
||
|
|
|
|
|
|
• |
• • • |
Апл |
|
D — Ь Д*тО Ал|Л п2 |
||||
Km — |
Km^npe/ , |
f2 |
• • • l n |
, |
где R lt R2, .... Rm — характеристики качества поверхности, значения которых необходимо получить на данной операции; klt k2, .... km; k10, •••» *mu — коэффициенты регрессии, постоянные для данной опера ции; tu t2.......tn— технологические факторы обработки.
Аналогично может быть выражена зависимость для критерия оп
тимальности — эксплуатационного |
свойства |
(износостойкости, кон |
||||
тактной жесткости, усталостной |
прочности |
и т. д.) |
|
|||
f I |
— £> Л ки0 |
/*«1/*и2 |
ftu n |
(3.130) |
||
и |
— киа пред£1 |
*2 |
• • Л п |
• |
В выражениях (3.129) и (3.130), кроме технологических факторов tu
t2, ..., tn, неизвестным является значение Д пред, которое можно |
полу |
чить через коэффициенты технологической наследственности ад и для |
|
рассматриваемого метода обработки из зависимости |
|
Д = а А р е д - |
(3.131) |
Учитывая, что коэффициент ад экспериментально может быть получен в виде
|
ад = кьфЧ***. . . |
|
(3.132) |
|
после преобразований |
получаем |
1 |
|
|
|
|
|
(3.133) |
|
л '»“ |
" |
( « А ,,* » 4 |
ьь |
|
ж г |
|
|||
Подставим значение Д„ред в выражения (3.129) и (3.130) |
||||
|
|
feio |
|
fkln |
«1 = *х kЛ/*Д1<*Д2 |
t^ti" |
|||
д 1 |
2 |
|
|
-^?2 — ^2 *А**Л1^Л2
л 1 2
*" “ *»/ ^ ? 5 Г — X 4l'i
/*Лп Ьл |
Мп |
(3.134) |
1п |
I |
|
л |
|
|
femO |
|
|
д : ? * г1*%‘"1 •■■<” .
,,
ftuQ
U - k |
I_______-_______Y*~ /W«2 |
(3.135) |
|
u |
&/*Д1/*Д2 |
/*Дл I |
|
|
л A1 ’2 |
|
|
После преобразования этих выражений получим следующую систему уравнений:
(3.136)
, А ,-ТГ^ |
lu |
кА \кт О \ |
( и |
кАпктО \ |
R„ |
|
|
|
|
и выражение для оценочной функции |
|
|
||
, . .— |
/. |
яал«о \ |
/ .. |
*Дп*ио;\ |
и- *■(£)*'.* |
|
|
(3.137) |
|
|
|
|
|
Приводим полученные выражения к линейному виду bjO
1л Rj = In |
k A \k jO |
\ KJl |
|
|
bA |
+ ^ / n |
^A*^/0 |
) l n < „ , |
|
|
ьд |
ku0 |
|
4
(3.138)
r f H |
Л01 |
* д А о |
\ |
|
+ |
|
|||
6A |
/ |
|
|
||||||
|
|
|
|
||||||
|
|
|
|
|
|
|
|||
4 " |
|
^An^uO |
„ . |
|
|
|
(3.139) |
||
|
1 |
6Д |
W |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
После введения переменных |
хг = |
In |
tu |
х2 = In /2» •••» |
= |
In |
|||
оценочной функции у = In |
(У и обозначения |
С/ = |
In /?у — In |
X |
|||||
*д/ — */0 |
„ |
_ |
/, |
^д< |
*|Д |
|
|
||
-----1-----— |
"Hi |
|
|
|
*Д |
’ |
|
||
|
|
|
|
|
|
|
|
получим систему уравнений типа |
|
П |
|
С ,= Е р й |
(3.140) |
{=1 |
|
и уравнения оценочной функции |
|
п |
|
У = Ру + (S PytXf |
(3.141) |
ции [11]. Так, максимизация У1р У[Х1 эквивалентна минимизации f=i
п
£ (—pyt) х(. Кроме того, любое линейное уравнение, а также любую си- i=1
стему линейных уравнений можно представить в виде некоторой сово купности неравенств с помощью одного дополнительного ограниче ния, т. е. систему уравнений
|
t p / i X ^ C , |
|
(3.142) |
|
можно записать в виде |
|
|
|
|
Е |
PitXi < С/, |
2 a fa < |
Р, |
(3.143) |
1=1 |
|
(=[ |
|
|
т |
т |
|
|
|
где а , = Ц р н , р = — |
2 С / . |
|
|
|
/=I |
/=! |
|
фигурирующей |
в той |
Когда некоторое i значение переменной хь |
или иной линейной модели, не ограничено в знаке, в процессе нахож дения численного решения оказывается полезным провести переход от переменных, не имеющих ограничения в знаке, к неотрицательным переменным. Так, если переменная х£для i = 1, 2, ..., k ^ п не огра-
п
ничена в знаке, то ограничения 2 Pnxi = С/ (j' = 1, 2, ..., т) преобра