CC BY
25
ВА ГРЕЧИШНИКОВ, ВМЖОЕАРЕВ, КАШПЕЙШЩ АМЦЫБУЛЬСЖИЙ
чистной. Такая схема подобна групповой, применяемой в протяжках для обработки шлицевых отверстий [1].
Успешные производственные испытания сверл разработанной конструкции показали следующие их основные достоинства:
1) высокая прочность взаимозаменяемых резцов, простота их заточки с большими запасами на переточку обеспечивают большой срок их службы и высокую надеж-' ность;
2) наличие большого пространства для размещения стружки предотвращает ее паке= тирование и позволяет работать с повышенными подачами;
3) возможность регулировки сверл по диаметру сокращает число их типоразмеров;
4) для повышения скорости резания возможно использование резцов с напайными пластинами из твердых сплавов марок T5K1Q, ВК8 или ТТ10К8Б.
СПИСОК ЛИТЕРАТУРЫ
.. 1 . Протяжки для обработки отверстий/ Д.К.Маргулис, М.М.Тверской, В.Н.Ашихмин и др . -М. Машиностроение, 1986. -232 с.
Томский политехнический университет
УДК 621.951.41
В.А.ГРЕЧИШНИКОВ, В.ИЖОКАРЕВ, Е.А.КОПЕЙКИН, А.ЮДЫБУЛЬСКИЙ
РАСЧЕТ НА КРУЧЕНИЕ СТЕБЛЯ СВЕРЛА ДЛЯ ГЛУБОКОГО ВИБРАЦИОННОГО СВЕРЛЕНИЯ МЕТОДОМ КОНЕЧНЫХ ЭЛЕМЕНТОВ
Рассмотрена методика расчета на кручение стебля сверла для глубоких отверстий методом конечных элементов.
При сверлении глубоких отверстий, в особенности малых диаметров, наибольшую трудность представляет сверление вязких материалов. Образующаяся при этом сливная стружка имеет форму гофрированной ленты большой длины, поэтому возможно ее пакетирование и защемление в стружечной канавке сверла. В результате на поверхности отверстий появляются задиры, а защемление стружки может привести к поломке инструмента.
Для повышения стабильности работы технологической системы и производительности обработки в таких случаях применяют механическое или кинематическое дробление стружки, а также используют адаптивное глубокое сверление [1]. Увеличение жесткости стебля сверла также в значительной степени способствует повышению стабильности работы инструмента и производительности обработки. Однако оценить напряженно-деформированное состояние сверла для глубокого сверления в его различных сечениях, оснащенного каналами для подвода смазочно-охлаждающей жидкости (СОЖ) и имеющего, как правило, сложный профиль, классическими теоретическими методами в ряде случаев крайне трудно.
В таких случаях весьма эффективны численные методы расчета и, в частности, метод конечных элементов (МКЭ). Являясь приближенным методом, МКЭ может обеспечивать приемлемую точность расчетов при соответствующем разбиении рассматриваемой области на конечные элементьт Значительная трудоемкость расчетов вызывает необходимость составления вычислительных программ и использования вычислительной техники. В то же время алгоритм МКЭ носит достаточно стандартный и во многом универсальный характер, что значительно облегчает решение задач на ЭВМ [2, 3].
Рассмотрим последовательность решения плоской (двумерной) задачи кручения стебля сверла при разбиении сечения сверла на конечные элементы треугольной формы.
Профиль поперечного сечения стебля сверла для глубокого вибрационного сверления с прямыми стружечными канавками представлен на рис. 1,а [4]. С учетом симметрии профиля на конечные элементы достаточно разбить лишь четвертую часть профиля (рис. 1,6).
14 20 26 32
б
Рис. 1. Поперечное сечение сверла для глубокого вибрационного сверления: а - общий вид; б - разбиение % части сечения на конечные элементы
Примем постоянными и известными модуль сдвига материала б (МПа), относительный угол закручивания 0 (рад/1000мм), размеры сечения стебля сверла. Неизвестными считаем крутящий момент М (Н-мм), обеспечивающий заданный угол закручивания5 касательные напряжения в сечении т2Х? т2у (МПа).
Для решения задачи используем мембранную аналогию Прандтля, согласно которой вводимую функцию напряжений ср можно представить как некоторую упругую поверхность, охватывающую поперечное сечение стержня (рис. 2). Здесь крутящий момент пропорционален объему, охватываемому этой поверхностью, а касательные напряжения в произвольной точке сечения можно определить по формулам [2, 5]
т т
ах'
(1)
Дифференциальное уравнение для функции напряжений можно представить в виде
Э2ср Э2ср
йх2 ду2
- + 209 = 0.
(2)
В.А.ГРЕЧИШНИКОВ, ВЖКОКАРЕВ, КАЖОЛЕЙКИН, АШЦЫБУЛЬСКИЙ
Граничными являются условия: на наружном контуре <р = 0; (3)
на внутреннем контуре ср = const.
Поскольку в стебле сверла имеется отверстие для подачи СОЖ, необходимо часть воображаемой мембраны заменить горизонтальной неупругой поверхностью (пластиной) с углом наклона, равным нулю [5].
Представим функцию напряжений ср в виде линейной функции, для каждого, треугольного элемента
ф = щ + а2х + «Зу, (4)
где а15а2,а3 - коэффициенты.
Формулу (4) можно представить в преобразованном виде для произвольного элемента с узлами i, j, m, обозначенными в направлении Рис.2. Функция напряжений В виде упругой ПО- против часовой стрелки (в противном случае
верхности для сечения стержня с отверстием г
- г г значения площади элементов будут отрица-
тельными)
Ф = N.Oi + N/Dj + NmOm = [N] {Ф}, (5)
где Ni?NJ?Nm- вводимые функции формы соответственно в узлах i, j, m; Ф19Ф-,Фт - значения функции напряжений соответственно в узлах i,j,m; [N] и {Ф} - соответственно вектор-строка узловых значений функции формы и вектор-столбец узловых значений функции напряжений
Функция формы, например, для узла i
Ni=^(ai+biX + Ciy)?
где А - площадь треугольника ijm; щ3Ь13щ - коэффициенты.
ai = xjym-xmyj;
bi=yj-ym;
3=xm-xf
Площадь треугольного элемента может быть определена по формуле
(6)
(7)
А = -2
1
х>
х-
(8)
Уi Уз
Хш Уш
Суммарное значение функции ф с учетом множества конечных элементов в системе
9 = 1ф(е\ (9)
е=1
где ф(е)- значение функции в элементе е; Е-число элементов в системе.
РАСЧЕТ НА КРУЧЕНИЕ СТЕБЛЯ СВЕРЛА ДЛЯ ГЛУБОКОГО ВИБРАЦИОННОГО СВЕРЛЕНИЯ
МЕТОДОМ КОНЕЧНЫХ ЭЛЕМЕНТОВ
_211
При реализации МКЭ применяем минимизацию некоторого функционала связанного с выражением (2), при вариационном рассмотрении задачи [2]
^ ¿([кСе)] {Ф} + {Г(е)» - 0 5 (Ю)
О е=1
где [к(е)]- матрица жесткости элемента е; {Ф}~ вектор узловых значений функции напряжений; {£ ^} - вектор нагрузки для элемента е.
Глобальная матрица жесткости элементов
[К] = Х[к(е)]. (п)
е=1
Глобальный вектор нагрузки
Тогда
* (12)
е=1
[К]{Ф} = {Р}. (13)
Матрица жесткости элемента е
[к(е)]= | [В(е)]т[Б(е)][В(е)]аУ, (14)
у(с)
где [В(е)]т- матрица, транспонированная к матрице [В(е)]; [В(е)]- матрица, полученная дифференцированием функций формы МГе); [В(е)]- матрица упругих характеристик материала (для данного случая, согласно [ 2], [0(е)]=[1]).
Для определения матрицы градиентов [В(е)] продифференцируем функцию формы [ф(е) ] по х и у и представим в виде
'ъ<в) ь2е) ... ь«У
[В(е)]= 1 / 1 ] 2А(е)
г(е) г(е) (е)
2 ••• Сп J
(15)
где и - число узловых точек в системе.
Коэффициенты, составляющие матрицу преимущественно имеют нулевые
значения, за исключением тех, которые определяют узлы элемента е. Например, для первого элемента (рис. 2) ненулевыми являются коэффициенты Ь1?Ь2,Ь5,с1?с2,с5, Коэффициенты для первого элемента определяем по формулам
Ь1 =у2-у5; с,=х5-х2;
ь2 =у5-У1; с2 = х1-х5; -(1б)
Ь5=У1-У2; с5=х2-х1.
Принимая единичной толщину рассматриваемой области и учитывая, что [В(е)]=[1], приводим выражение (14) к виду
[к(е)] = [В(е)]т[В(е)]А(е). ■ (17)
Вектор нагрузки для элемента е
В.А.ГРЕЧИШНИКОВ, В. И. КО КАР ЕВ, Е.А.КОПЕЙКИН, А.Ю.ЦЫБУЛЬСКИЙ
{f(e)}= f2G(e)0[N(e) ]TdV. (18)
v(e)
Например, для первого элемента с учетом того, что сумма функций формы для любой точки элемента равна единице [2], вектор нагрузки
ff(1)} = [2G(1)6[N(1)]TdV = -G(i)9A(i)(l 1 0 0 1 0 ... 0)т. (19)
yh 3
После определения [к(е)] и {fve)} для всех элементов рассматриваемой области производим вычисление глобальных матрицы жесткости и вектора нагрузки по формулам (11) и (12), преобразуем систему уравнений (13) с учетом граничных условий (3) и решаем ее относительно {Ф}.
Преобразование системы уравнений с учетом первого граничного условия (ср = 0 на наружном контуре) не вызывает затруднений, поскольку для этого необходимо лишь исключить из этой системы уравнения и элементы столбцов глобальной матрицы жесткости, определяющие узлы сечения, расположенные на наружном контуре (1, 5, 9, 15, 21, 27, 30-35).
Для преобразования системы уравнений с учетом второго граничного условия (ф = const на внутреннем контуре) необходимо определить значение функции напряжений ф по контуру отверстия. Данная задача решалась следующим образом: дополнительно разбивалось на конечные элементы сечение отверстия, определялись значения функции напряжений в узлах для сплошного сечения (с одним граничным условием ф = 0 на наружном контуре), находилось среднее значение функции в узлах, расположенных по внутреннему контуру, а затем вновь составлялась и решалась система уравнений относительно {Ф} при известных граничных условиях (3).
В дальнейшем вычислялись напряжения tzx и т2ув соответствии с (1), которые считаются постоянными в пределах одного элемента
(^)Т=[В«]{Ф>. (20)
дх ду
Крутящий момент, воспринимаемый элементом е
Т(е) = 2 J ф(е)(1А - 2 |[М(е)]{Ф}ёА =2А(е)[Ф]т(М<е) N(2e) ... Mff . (21)
А(«> А<е>
Для первого элемента, например, Т® =2А(1)[Ф]Т(Ц1} - ^)=^А(1)[Ф]Т(1 10 0 10... 0)т A(1)(®j +Ф5).(22)
Крутящий момент, воспринимаемый рассматриваемой областью сечения (без учета влияния на момент наличия отверстия), равен удвоенному объему, охватываемому воображаемой упругой поверхностью
т = Ет(е)- • (23)
е=1
Суммарный момент, воспринимаемый сечением стебля сверла, с учетом того, что разбиению подвергалась лишь четвертая его часть, а также с учетом влияния на величину момента наличия отверстия
РАСЧЕТ НА КРУЧЕНИЕ СТЕБЛЯ СВЕРЛА ДЛЯ ГЛУБОКОГО ВИБРАЦИОННОГО СВЕРЛЕНИЯ
МЕ ТО ДОМ КОНЕ ЧНЫХ ЭЛЕМЕНТОВ
М = 4Т + 2тхг2Ф0ТВ5 (24)
где г-радиус отверстия ; Фотв - значение функции напряжений на границе отверстия.
Площадь отверстия сложного профиля определяется как сумма площадей конечных элементов, на которые разбивается сечение отверстия. В связи с аппроксимацией функции напряжений значение крутящего момента, полученное в результате численного эксперимента, отличается от точного решения в меньшую сторону, причем точность расчета повышается при использовании более мелкой сетки разбиения области на конечные элементы.
Крутильная жесткость стебля сверла постоянного сечения
С = М/0. (25)
При расчете стебля сверла с сечением, показанным на рис. 1, диаметром 7,2 мм с отверстием диаметром 2,0 мм и минимальной толщиной стенки 1,25 мм при 0=80000 МПа, 9=1 град/1000мм, при разбиении области на 130 элементов крутящий момент ориентировочно составил М=196 Н-мм, крутильная жесткость С=11,26 Н-м2. Максимальные касательные напряжения, зафиксированные в одном из элементов, т2Х =0,0125 МПа, г^ -5,8 МПа.
Экспериментальные исследования на кручение стебля сверла (сталь 20Х) с указанными конструктивными параметрами показали, что значение жесткости, полученное расчетом по МКЭ (11,26 Н-м2) при числе элементов п = 130, меньше среднего экспериментального значения (12,12 Н-м2) на 7,1 %.
Расчет стебля сверла на продольный изгиб может быть выполнен при использовании конечно-элементной модели, если момент инерции сечения стебля сверла относительно оси у представить как сумму моментов инерции конечных элементов.
Использование ЭВМ значительно сокращает трудоемкость расчетов, позволяет оперативно произвести анализ напряженного состояния стебля сверла сложного сечения, сравнивать различные конструкции на жесткость, производить оптимизационное проектирование.
СПИСОК ЛИТЕРАТУРЫ
1. Обработка глубоких отверстий / Н.Ф.Уткин, Ю .И. К и ж н я е в, С.К.Плужников и др.; Под общ. ред. Н.Ф. У тки на.-Л.: Машиностроение. 1988.- 269с.
2. Сегер линд Л. Применение метода конечных элементов. ~М.:Мир, 1979. - 388 с.
3. Розин Л . А. Задачи теории упругости и численные методы их решения.- Спб.: Изд-во СпбГТУ, 1998.- 531с. ' '
4. Прогрессивные режущие инструменты и режимы резания - металлов: Справочник/ В И.Баранчиков, А.В.Жаринов, Н.Д.Юдина и др. ; Под общ. ред. В И.Баранчикова.-М.: Машиностроение, 1990.-400с.
5. Тимошенко С.П., Гудьер Д ж. Теория упругости. - М: Наука, 1979. -531 с.
Московский государственный университет «СТАНКИН»
