Назначение оптимального режима резания

Использование электронно-вычислительных машин дает возможность определить наиболее оптимальные режимы резания. Для решения таких задач удобно воспользоваться методом линейного программирования, поскольку основные формулы теории резания выражаются степенными зависимостями, которые после логарифмирования превращаются в линейные зависимости. Одной из главных задач указанных расчетов является создание математической модели, наиболее точно описывающей основные закономерности процесса резания, которая затем реализуется на ЭВМ. Для установления математической модели составляются основные уравнения технических ограничений и одно уравнение, подлежащее оптимизации. В качестве критерия оптимальности принимается себестоимость или производительность операции. Уравнение для определения критерия оптимальности как функция элементов режимов резания называется оценочной или целевой функцией. Примем в качестве оценочной функции

f=ns=max (16.32)

и запишем уравнения технических ограничений, имеющих место при работе станка.

1. Режущая способность инструмента:

v=(CvKvD)/(TmtXvsYv)=(pDn)/1000 (16.33)

отсюда

nsYv=(318CvKv)/(DTmtXv) (16.34)

2. Эффективная мощность станка:

Ne=(Pzv)/6120£Ne ст (16.35)

Выразив скорость резания через D и n, получим:

(CPKPtXpsYppDn )/(6120×1000) £Ne ст (16.36)

отсюда

nsYp£(195×104Ne ст)/(CpKpDtXp) ( 16.37)

3. Допустимая шероховатость обработанной поверхности:

Rz=(Cr tXrsYrjZrj1Zr)/rqr£ Rz доп (16.38)

отсюда

sYr£(Rz доп rqr)/(Cr tXr(jj1)Zr) (16.39)

4. Максимальная сила, допустимая прочностью слабого звена механизма подачи станка:

Px=CPxKPtXpxsYpx £ PМП (16.40)

отсюда

sYpx£ PМП /(CPxKPtXpx ) (16.41)

5. Минимальная подача станка: S³S min

6. Максимальная подача станка: S£Sст max

7. Минимальная частота вращения шпинделя станка: n³n ст min.

8. Максимальная частота, вращения шпинделя станка: n£n ст max.

Чтобы систему ограничений и оценочную функцию привести к линейной форме, прологарифмируем полученные выражения, предварительно умножив для удобства вычисления во всех выражениях подачу на 100. В результате получим:

ln n+Yv×ln(100s)=ln((318×100YvCvKv)/(DTmtXv)) (16.42)

ln n+YPln(100s) £ln((195×104×100YpNe ст)/(CPzKPDtXp) (16.43)

Yr ln(100s) £ln((100Yrхrqr)/CrtХr(jj1)Zr) (16.44)

YPx ln(100s) £ln((100YPМП)/CPxKP tХPx (16.45)

ln(100s)³ln(100sст min) (16.46)

ln(100s)£ln(100sст max) (16.47)

ln n³ln £n ст min (16.48)

ln n£ln n ст max (16.49)

f=ln n + ln(100s)=max (16.50)

Обозначив ln n=x1; ln(100s)=x2 и правые части выражений через bi, получим математическую модель оптимального режима резания

Назначение оптимального режима резания - student2.ru
Назначение оптимального режима резания - student2.ru

Рис. 16.2. Геометрическая интерпретация математической

модели оптимального режима резания

Графическая интерпретация математической модели оптимального процесса резания представлена на рис. 16.2. Линии I-VIII изображают уравнения ограничений.

Штриховой линией, наклоненной к оси абсцисс под углом 45°, изображена оценочная функция. Область возможных решений системы очерчена выпуклым многоугольником, координаты вершин которого являются корнями совместного решения уравнений системы А. Чтобы найти оптимальное решение среди многих решений системы ограничений, необходимо среди точек многоугольника ABСDE найти такие, для которых f=max. Многоугольник ABСDE ограничен линиями минимальной частоты вращения шпинделя и минимальной подачи, линией подачи по допустимой шероховатости, линиями режущей способности инструмента и эффективной мощности станка. При перемещении прямой IX параллельно самой себе из начала координат в точку А, функция f будет расти и достигнет максимума в точке С, координаты которой X1oпт и Х2опт будут соответствовать оптимальному решению системы.

При решении задачи на ЭВМ, прежде всего определяют координаты X1 и Х2 всех точек пересечения прямых. Далее выясняют, какие из значений Х1 и Х2 удовлетворяют всем уравнениям системы. Затем определяют координаты вершины прямоугольника, для которой эта сумма имеет максимальную величину

X1oпт2опт=max (16.51)

Так как

X1oпт=ln nопт; Х2опт=ln(100sопт); (16.52)

то

n=eX1опт; sопт=0,01еХ2опт. (16.53)

Глава 17

Основные понятия теории планирования экспериментов.

Наши рекомендации