Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Оптимизация технологических процессов механической обработки

..pdf
Скачиваний:
41
Добавлен:
15.11.2022
Размер:
8.74 Mб
Скачать

ного штучного времени /шгр

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 — Р12) 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 постоянна на ней и поэтому любую ее точку можно принять за минимум, в частности точку пересе­ чения кривой с ребром многогранника, т. е. выполнение равенства Р12 = «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' =

К

 

№ Ч г'= К г,

(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

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

W

 

 

 

 

 

 

 

 

 

 

 

 

 

После введения переменных

хг =

In

tu

х2 = In /2» •••»

=

In

оценочной функции у = In

(У и обозначения

С/ =

In /?у — In

X

*д/ — */0

_

/,

^д<

*|Д

 

 

-----1-----—

"Hi

 

 

 

 

 

 

 

 

 

 

 

 

Для оптимизации полученной математической модели могут быть при­ менены численные методы решения задач линейного программирова­ ния. Математически эту задачу оптимизации можно сформулировать следующим образом. Пусть хс — /-я управляемая переменная (/ = 1, 2, 3, п). Требуется определить такие значения х£, чтобы выражение
п
5] PyiXi в зависимости от содержания задачи (вида оптимизируемого
{=\
эксплуатационного свойства) было максимизировано или минимизи­ ровано с учетом наложенных ограничений на изменение управляемых переменных х£.
В линейном программировании любая задача максимизации может быть сведена к эквивалентной задаче минимизации (и наоборот), если одновременно с изменением «знака» оптимизации произвести изменение знаков перед всеми коэффициентами в выражении для целевой функ-
п

получим систему уравнений типа

 

П

 

С ,= Е р й

(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, ..., т) преобра­