CC BY
51
Технология
УДК 621.753
АЛГОРИТМ РАСЧЕТА ОПЕРАЦИОННЫХ РАЗМЕРОВ ПРИ ОБРАБОТКЕ ПОВЕРХНОСТЕЙ ВРАЩЕНИЯ
М.М. Тверской
THE ALGORITHM OF CALCULATION OF THE DIMENSIONS FOR MACHINING OPERATIONS OF SURFACE OF REVOLUTION
M.M. Tverskoy
Алгоритм расчета операционных размеров при многооперационной обработке поверхности вращения основан на составлении и решении векторной размерной цепи для каждой из операций, начиная с последней. Многозвенная векторная размерная цепь приводится к трехзвенной линейной размерной цепи, замыкающим звеном которой является припуск. Составляющими звеньями этой цепи являются радиусы поверхностей, получаемых на двух смежных переходах, и результирующий эксцентриситет между этими поверхностями. Радиус поверхности, получаемой на предшествующем переходе, неизвестен и подлежит расчету. Приводится алгоритм вычисления вероятностных числовых характеристик результирующего эксцентриситета и последующего расчета операционного размера. Расчет по описанному алгоритму продолжается до определения операционного размера заготовительной операции.
Ключевые слова: операционный размер, размерная цепь, эксцентриситет, припуск.
The algorithm of calculation of the dimensions for machining operations during multi-operational machining of surface of revolution is based on making-up and solving of vector dimensional chain for each operation starting with the last one. Multi-component vector dimensional chain is transformed to three-component linear dimensional chain, in which allowance performs as closed component. Partial components of this chain are radiuses of the surfaces, which were made during two neighboring operations, and the resulting eccentricity between these surfaces. The radius of the surface, which was made during previous operation, is unknown and has to be calculated. The algorithm of calculation of probabilistic numerical characteristics of the resulting eccentricity and the following calculation of the dimension for machining operation is shown. The calculations according to the algorithm described are going on until the dimensions for producing the metal bar are found.
Keywords: dimension for machining operation, dimensional chain, eccentricity, allowance.
В работе [1] изложена методика вероятностного расчета размерных цепей, содержащих звенья-эксцентриситеты. В качестве математической модели эксцентриситета принят двух- или трехмерный случайный вектор, алгебраические проекции которого на оси координат есть случайные величины, имеющие нормальное распределение (Гаусса) с нулевым математическим ожиданием
и равными дисперсиями а2х = а2 = a2z = а2. Величина ст в этом случае называется радиальным
средним квадратическим отклонением (СКО).
Рассмотрим двухмерный случай (рис. 1). Производится окончательная обработка цилиндрической поверхности заданного радиуса г. Для расчета размера R, выдерживаемого на предшествующей операции (переходе) составляется векторная размерная цепь
Z = R-r-Eb, (1)
—У
где Z - вектор припуска, нормальный к поверхности, получаемой в результате обработки на рас-сматриваемом технологическом переходе; г - радиус-вектор этой поверхности; R - радиус -
вектор поверхности, получаемой на предыдущем переходе, он подлежит определению; Еа -
-» —>
вектор эксцентриситета поверхности с радиусом г относительно оси поверхности с радиусом 7?.
—>
Вектор Еа есть разность векторов, имеющих начало в точке О, соответствующей оси исходной базовой поверхности:
Еа =Er~ER.
(2)
Учитывая, что каждый из векторов Ег и
—>
Ей в общем случае есть векторная сумма нескольких слагаемых (см. рис. 1), выражение (2) запишем в следующем виде:
Еа =Е\ + Ег + Еъ—Еа—Еъ. (3)
Суммирование независимых случайных —>
векторов Еі в правой части выражения (3) производим, определяя числовые характеристики суммы их алгебраических проекций на оси прямоугольной системы координат. При этом в силу устойчивости закона нормального распределения [2] алгебраические проекции —»
суммарного вектора Еа , рассматриваемые как случайные величины, будут иметь также нормальное распределение с математическим ожиданием, равным нулю, и дисперсиями, равными
сумме дисперсий алгебраических проекций суммируемых векторов Еі :
(4)
;=i
где а, - радиальные СКО суммируемых векторов; п - число суммируемых векторов в правой части (3); <7Д - радиальное СКО суммарного вектора.
Поскольку суммируемые векторы Ei имеют круговое нормальное распределение
—>
(oxi = сyi = Gj), то такое же распределение имеет и суммарный вектор Еа, у которого
СХА = °УА = °А ■
—>
В этом случае модуль вектора Еа распределен по закону Рэлея (или по закону Максвелла для трехмерного вектора).
Поле рассеивания скалярной случайной величины £, =
, подчиняющейся закону Рэлея,
(5)
принимается равным со, = 2aEl4l, где оЕ. - СКО случайной величины Et.
Величина СКО ст£ связана с радиальным СКО ст, (СКО алгебраических проекций случай-
ного вектора Еі ) следующей зависимостью
°Еі
= сг J2-
п
(6)
Учитывая, что со,- =£',тах ~ЕШп , при ЕШа =0 получаем со,- =£',тах . Тогда из соотношений (5) и (6) получаем
(7)
2,
Подставляя соотношение (7) в выражение (4), получаем 1
(8)
Величина стд является искомым (и единственным) параметром закона Рэлея, которому подчиняется распределение случайной величины ЕА.
Для нахождения СКО этой случайной величины вновь используем соотношение (6):
°ЕА =0,189■ (9)
Соотношение (9) позволяет вычислять СКО величины эксцентриситета (ЕА) поверхности с
радиусом г относительно оси поверхности с радиусом R, т. е. СКО модуля вектора Еа , входящего в уравнения (1) и (2). Исходными являются величины Eijmx, регламентируемые в конструкторской и технологической документации.
Находим математическое ожидание случайной величины 2?д, используя соотношения закона Рэлея:
mEA~°EA,L =°ЕАлН^-=:1’913(УЕА- (10)
я/2. _ 1 л
12-л/2 °Ел\4-л
Возвратимся к рис. 1 и уравнению (1), из которых следует, что величина припуска (модуля
—>
вектора Z) различна в разных точках обрабатываемой поверхности с радиусом г. Направление —^
вектора Еа определяет ту точку (В) на обрабатываемой поверхности, где припуск имеет минимальную величину. Это справедливо для любой реализации процесса обработки. Для указанного направления все векторы, входящие в уравнение (1), коллинеарны, поэтому это векторное уравнение вырождается в алгебраическое:
Z = R - г - ЕА . (11)
Уравнение (11) - это уравнение линейной размерной цепи (рис. 2). Рассматривая все величины, входящие в это уравнение, как случайные, находим математическое ожидание замыкающего звена - припуска
mZ~mR~mr~mEA-
Среднее квадратическое отклонение величины припуска
W + 4a , (13)
где aR = ^ ТD ; стг = Td ; TD, Td - заданные допуски на диаметры соответствующих поверхностей.
Для определения математического ожидания mR неизвестного составляющего звена R необходимо знать величину математического ожидания mz замыкающего звена - припуска Z. Однако регламентируемым значением этого звена является его минимальное значение Zmin [1].
Случайная величина Z является суммой трех случайных величин, входящих в правую часть уравнения (11). Плотности распределения радиусов R и г описываются законом Гаусса, плотность распределения эксцентриситета ЕА - законом Рэлея. Разность радиусов (R- г) распределена также по закону Гаусса [2], тогда закон распределения случайной величины Z есть композиция
cjz = V°
законов Гаусса и Рэлея. Вид кривой плотности распределения/(г) зависит [3] от величины отношения средних квадратических отклонений
1 =
’Е А
VCT«+(Tr
(14)
При ст£Д = 0, т. е. при Еа = const имеем композицию двух законов Гаусса, которая снова дает закон Г аусса со средним квадратическим отклонением
°z =л/стк+стг ■ (15)
С вероятностью 99,73 % минимальное значение припуска будет на расстоянии Зст2 от его математического ожидания.
Если = стг = 0, т. е. R = const, г = const, то распределение плотности замыкающего звена -припуска будет описываться законом Рэлея. Поле рассеивания припуска с вероятностью 99,76 % принимаем равным 2 аЕд V7 = 5,2915а£д (см. выражение (5)). Математическое ожидание припуска отстоит от минимального значения припуска на расстоянии
f I-----л
a = 2aF yfl -
ТПр — (Уг
V 4-я
= 3,378a£ =3,378ctz.
(16)
Таким образом, в общем случае (см. рис. 2)
УЯг — ^ггап ко.
где к находится в интервале (3,0; 3,378).
При реальных (для рассматриваемых условий) значениях отношения / = (0______________4) величина к
изменяется от 3,0 до 3,25 [3]. Подставляя (16) в (12), получаем
(17)
mR
— тг + + Zmjn+& + СТГ +ст£д .
Для определения величины к может быть использована линейная аппроксимация к = 3+0,06 I, (18)
где I определяется отношением (14).
После вычисления по формулам (17), (18) значения математического ожидания тК искомого радиуса /? находим предельные значения операционного размера
0ЯГ = 2'%±^ГО. (19)
Рассмотрим пример. На чертеже детали задан размер 0 40 ± 0,06, что соответствует г = 20 ± 0,03; тг= 20; Td = 0,12; ог = 0,01. В технологической карте заданы отклонения ±0,12 на операционный размер 0D. Тогда TD= 0,24; aR = 0,02. Расчетное значение [1] минимального припуска Zmin=0,2. Величины допускаемых эксцентриситетов: £1пш =0,2; Е2тах = 0,15; £3пах=0,1; Е4тах = 0,2; £5шх =0,15 (см. рис. 1). По формуле (9) вычисляем СКО суммарного эксцентриситета
о£Д = 0,189-у/0,22 +0Д52 +0,12 +0,22 +0Д52 - 0,069443.
По формуле (10) подсчитываем математическое ожидание суммарного эксцентриситета.
mEh = 0,069443^^— » 0,1328.
По формуле (13) находим СКО припуска az = д/0,022 +0,012 + 0,0694432 = 0,07295.
По формуле (14) находим отношение СКО:
, 3.1.
V0,022 +0,012
По формуле (18) находим к = 3 + 0,06x3,1 = 3,186; koz = 0,2324.
По формуле (17) рассчитываем математическое ожидание mR= 20 + 0,1328 + 0,2 + 0,2324 к 20,5652.
По формуле (19) получаем операционный размер 0 41,1304 ±0,12.
Полученный размер округляем в большую сторону до того знака после запятой, который имеют отклонения, и окончательно получаем D = 0 41,14 ±0,012.
Если поверхность с радиусом г обрабатывается за несколько операций (переходов), то полученный размер D является исходным при расчете операционного размера для предшествующей операции (перехода). Расчет по изложенному выше алгоритму продолжается до получения операционного размера для заготовительной операции (например, штамповки).
Если обрабатываемая поверхность вращения является внутренней (отверстием), то алгоритм расчета операционных размеров остается прежним, но в размерной цепи (11) известным является звено R, а определяемым - звено г.
Литература
1. Тверской, М.М. Расчет размерных цепей со звеньями-эксцентриситетами / М.М. Тверской//Вестник ЮУрГУ. Серия «Машиностроение». - 2005. - Вып. 6. - № 1 (41). - С. 168-180.
2. Вентцель, Е.С. Теория вероятностей/Е.С. Вентцель. - М.: Наука, 1969. - 576 с.
3. Точность производства в машиностроении и приборостроении / под ред. А.Н. Гаврилова. -М.: Машиностроение, 1973. — 567 с.
Поступила в редакцию 22 февраля 2011 г.
Тверской Михаил Михайлович. Доктор технических наук, профессор, заведующий кафедрой «Автоматизация механосборочного производства», Южно-Уральский государственный университет (г. Челябинск). Область научных интересов - технология и автоматизация машиностроения. Тел.: (351)267-92-12.
Mikhail М. Tverskoy. The Doctor of engineering science, the head of “Mechano-Assembley Automation” department of South Ural State University, Chelyabinsk, Russia. Professional interests: technology and automation of engineering industry. Tel.: (351) 267-92-12.
