Применение экстраполяционных и итерационных методов в инженерных теплофизических расчетах

Тип работы:
Реферат
Предмет:
ТЕХНИЧЕСКИЕ НАУКИ


Узнать стоимость

Детальная информация о работе

Выдержка из работы

УДК 621. 9
А. А. Липатов, Ю. Л. Чигиринский
ПРИМЕНЕНИЕ ЭКСТРАПОЛЯЦИОННЫХ И ИТЕРАЦИОННЫХ МЕТОДОВ В ИНЖЕНЕРНЫХ ТЕПЛОФИЗИЧЕСКИХ РАСЧЕТАХ
Волгоградский государственный технический университет
leeandrej@yandex. ru- julio-tchigirinsky@yandex. ru
Рассмотрено усовершенствование методики расчета температуры на передней и задней поверхностях режущего инструмента профессора А. Н. Резникова, заключающееся в учете зависимости теплопроводности обрабатываемого материала от температуры путем реализации итерационного алгоритма. Дополнительно введено экспериментальное определение сил, действующих на площадке износа, методом экстраполяции на нулевую площадку износа. Представлено сравнение результатов расчетов для случаев линейной и нелинейной экстраполяции. Показано, что учет изменения теплопроводности и температуропроводности аустенитной стали с ростом температуры приводит к существенной коррекции результатов по сравнению с известной методикой: повышение теплопроводности и температуропроводности сопровождается снижением температуры на передней поверхности и существенным повышением — на площадке износа.
Ключевые слова: расчет температуры, итерация, зависимость от температуры, теплопроводность, экстраполяция на нулевую площадку износа, силы резания, точение, аустенитная сталь, передняя поверхность, задняя поверхность, площадка износа.
A. A. Lipatov, Ju. L. Tchigirinsky
APPLICATION EXTRAPOLATION AND ITERATIVE TECHNIQUE AT ENGINEERING THERMAL DESIGN
Volgograd State Technical University
The improvement of a technique temperature calculation of on fore and back surfaces of a cutter of prof. A. N. Reznikov, gathered in the account of a relationship of heat conductivity of a fashioned stuff from temperature by realization of iterative algorithm is viewed. Is extra the developmental determination of the forces, wears acting to a platform by a method of extrapolation on zero wear platform is introduced. Comparison of results of calculations for cases of streak and nonlinear extrapolation is presented. It is shown, that the account of a modification of heat conductivity and temperature conductivity austenitic to steel with growth of temperature leads to essential correction of results in comparison with a known technique: the raise of heat conductivity and temperature conductivity is accompanied by decrease in temperature on a fore surface and an essential raise — on a platform of a wear.
Key words: temperature calculation, iteration, temperature relationship, thermal conductivity coefficient, extrapolation on zero wear platform, cutting force, turning, austenitic steel, a fore surface, a back surface, wear platform.
Температура резания является важнейшим В настоящее время общепринятой для раз-фактором, определяющим вид контактного дельного определения температур на передней
взаимодействия и интенсивность износа. Со- поверхности и площадке износа является мето-
вершенствованию методики расчета темпера- дика профессора А. Н. Резникова. Подробное ее
туры на режущих поверхностях лезвийного ин- описание приведено в изданном достаточным
струмента и посвящена настоящая работа. тиражом первоисточнике [4]. Можно лишь от-
метить, что алгоритм расчета учитывает взаимный нагрев контактных площадок, а также их нагрев за счет теплоты, выделяющейся в зоне стружкообразования. Многие исходные данные для расчета (усадка стружки, длина контакта по передней поверхности, силы резания и др.) должны быть получены экспериментально. Эта методика и была принята за базовую.
Вместе с тем известный метод не лишен недостатков. Во-первых, базовая методика предполагает использование экспериментальных данных о составляющих силы резания (вертикальной Рг и горизонтальной Р^,), действующих на инструмент в целом. Распределение же их по контактным площадкам (передней поверхности и площадке износа задней поверхности) производится априори: исходя из предела прочности обрабатываемого материала и ширины площадки износа кз задаются силы (нормальная Ркп и тангенциальная Р^Т) на последней, затем вычисляются силы (также нормальная и тангенциальная) на передней поверхности (для нулевого переднего угла у это сводится к вычитанию: Руп Р2 Phт, Рут Рху Ркп).
Во-вторых, в расчете используются средние значения теплопроводности X и температуропроводности ю обрабатываемого материала в предполагаемом диапазоне значений температуры 0 (введение зависимостей Х (0) и ю (0) непосредственно в дифференциальные уравнения теплопроводности вызвало бы значительные трудности при их решении).
Для устранения первого из указанных недостатков авторами был разработан способ экспериментального выделения сил Ро и Рхуо, действующих на переднюю поверхность, из суммарного усилия резания путем экстраполяции экспериментально полученных зависимостей Р2(Из) и Рху (Из) на нулевую площадку износа [2]. После этого для конкретных значений Нз вычитанием определяются усилия на самой площадке износа: РЫ = Рху — Рхуо, Ры = Ре ~ Рео (силы Руп и Ру: на передней поверхности предполагаются неизменными и при у = 0 равными соответственно Ро и Рхуо).
При выявлении характера зависимостей Рг (кз) и Рху (кз) были использованы различные (линейная, степенная, экспонента и др.) функции. Статистическая обработка [3] результатов показала, что наилучшим образом (по критерию минимизации стандартного отклонения) зависимости Ре (Нз) и Рху (Из) могут быть представлены степенной функцией вида
Р = Л (Из + В) к + С с показателем степени К & lt- 1 (т. е. по мере увеличения площадки износа рост сил несколько замедляется). Для сравнения в расчетах была использована и реализуемая простейшими стандартными программами линейная экстраполяция (отклонение нелинейных степенных графиков от прямых не превышало 2−2,5%).
Для разрешения второй проблемы было решено, сохранив математический базис метода А. Н. Резникова, использовать для учета зависимостей теплофизических параметров стали от температуры процесс итерационного приближения. Алгоритм строился следующим образом. Для предварительного расчета температур в зоне стружкообразования 0с, на передней поверхности инструмента 0п и площадке износа задней поверхности 0з использовались рекомендуемые [3] средние значения тепло- и температуропроводности обрабатываемого материала. Затем производился их перерасчет для соответствующих зон деформации согласно полученным в первом приближении значениям температур (разумеется, на этой стадии величины X и ю для разных зон деформации становились различными). С использованием полученных таким образом значений теплофизических параметров обрабатываемой стали для зоны стружкообразования Хс и юс, передней поверхности Хп и юп и площадки износа Хз и юз вновь рассчитывались соответствующие температуры. По новым значениям 0с, 0п и 0з опять рассчитывались величины теплопроводности и температуропроводности для всех зон деформации и т. д. Ряды получаемых при этом значений 0с, 0п и 0з достаточно быстро сходились (для вычисления температур с точностью до 1 °C требовалось не более 6−7 итераций).
С целью проведения сравнительного анализа разработанная программа позволяла рассчитывать значения температур как по усовершенствованной, так и по базовой [4] методике. Естественно, обеспечивалась возможность задействовать два описанных выше усовершенствования по отдельности (в том числе использовать результаты как линейной, так и нелинейной экстраполяции сил).
Для построения зависимостей теплофизических параметров обрабатываемого материала от температуры Х (0) и ю (0) использовались экспериментальные точки, взятые из целого ряда литературных источников. Статистическая обработка справочных данных показала, что наи-
лучшим образом (с минимальным стандартным отклонением) указанные зависимости могут быть представлены степенными функциями. В результате были получены следующие эмпирические выражения для зависимостей теплофизических параметров стали 12Х18Н10Т от температуры:
Х = 0,972(0 + 965)0'036 -1,210,
ю = 0,089(0 + 23)0,038 -0,063.
Отличие рассчитанных по этим формулам величин, от рекомендуемых [4] средних значений (X = 0,054 кал/см-с-град, ю = 0,05 см2/с) в диапазоне температур 150−1000 °С оказалось достаточно существенным: для X — от -23 до +24%, для ю — от -9,5 до +5,6%.
Получение необходимых для расчетов исходных экспериментальных данных осуществлялось при продольном точении стали 12Х18Н10Т резцами из твердого сплава ТТ20К9 со скоростями резания 0,5−2,0 м/с при подаче 0,3 мм/об и глубине резания 1,5 мм. Геометрия резцов была следующей: передний угол у = 0, главный задний угол, а = 10, главный угол в плане ф = 45 (все эти параметры были учтены непосредственно в программе расчета). Величина площадки износа Нз варьировалась в пределах до 0,4 мм (методика экспериментов по измерению сил резания инструментом с различными значениями Нз приведена в [2]).
Наиболее характерные результаты для базовой и различных вариантов предложенной методики расчетов приведены в табл. 1 и 2.
Как следует из табл. 1, экспериментальное разделение действующих на передней и задней
поверхности сил приводит к незначительному снижению расчетной температуры передней поверхности и гораздо более существенному повышению температуры площадки износа, так как определенные из опыта силы на площадке оказались значительно больше заданных априори. При этом температура 0з оказалась заметно выше в случае использования более точной нелинейной экстраполяции. Учет зависимостей Х (0) и ю (0) также приводит в большинстве случаев к возрастанию расчетной температуры на площадке износа: из-за меньшей (чем на передней поверхности) температуры на площадке износа теплопроводность и температуропроводность стали оказываются меньше рекомендуемых средних значений, что затрудняет сток тепла в заготовку. Более высокие температуры передней поверхности приводят, наоборот, к возрастанию величин X и ю для контактных объемов стружки, что приводит к интенсификации стока тепла в обрабатываемый материал и снижению расчетных контактных температур (за исключением скорости резания 0,5 м/с, при которой температура 0п еще невелика). В целом использование предложенной методики по сравнению с базовой фиксирует снижение температуры передней поверхности инструмента и рост (весьма существенный) температуры на площадке износа. Этот результат лучше согласуется с выявленными ранее [1] особенностями контактного взаимодействия при обработке аусте-нитной стали (исчезновение застойных явлений на площадке износа при меньших скоростях резания, чем на передней поверхности).
Таблица 1
Расчетные температуры на контактных поверхностях инструмента при обработке с различными скоростями (ширина площадки износа ^ = 0,2 мм)
Методика расчета Скорость резания V, м/с
0,5 0,75 1,0 1,5 2,0
Расчетные температуры на передней поверхности (0п) и площадке износа (0з)
0п 0з 0п 0з 0п 0з 0п 0з 0п 0з
Базовая [3] 660 363 762 386 825 395 927 415 1024 442
Экспериментальное определение сил
(линейная экстраполяция), X и ю —
константы 649 430 750 461 813 473 915 498 1010 538
Экспериментальное определение сил
(нелинейная экстраполяция), X и ю —
константы 644 473 744 514 805 534 906 578 1000 628
Экспериментальное определение сил
(линейная экстраполяция), учет зави-
симостей Х (0) и ю (0) 658 456 742 482 790 492 865 512 935 546
Экспериментальное определение сил
(нелинейная экстраполяция), учет за-
висимостей Х (0) и ю (0)* 653 495 736 530 783 546 858 581 926 621
* Окончательный вариант предложенной методики.
В табл. 2 приведены данные об изменении расчетной температуры на площадке износа в зависимости от самого износа (величины Кз) при обработке со скоростью V =1,5 м/с. Помимо уже выявленных при анализе табл. 1 закономерностей, данные табл. 2 показывают достаточно очевидную тенденцию роста температуры 0з по мере износа инструмента.
Также обращает на себя внимание следующее. Представление зависимостей Р2(Кз) и Рху (К) степенными функциями (как было указано выше, более точное, с меньшим стандартным отклонением) фиксирует существенное изменение характера роста 0з при износе инструмента. Уже при минимальном износе 0,05 мм температура площадки достигает достаточно высоких значений, приближающихся к 500 °C. При дальнейшем возрастании Кз рост 0з замедляется. Для расчетов с использованием линейного представления зависимостей Р2(Кз) и РхуЮ фиксируется более равномерный рост 0з при увеличении износа Кз. Это также может быть объяснено присущим аустенитной стали характером зависимостей теплофизических параметров X и ю от температуры.
Таблица 2
Расчетные температуры на площадке износа задней поверхности в зависимости от величины hз- v =1,5 м/с
Методика расчета Величина площадки износа, Из, мм
0,05 0,1 0,2 0,3 0,4
Базовая [3] 338 369 413 448 479
Экспериментальное определение сил (линейная экстраполяция), X и ю — константы 386 431 498 549 592
Экспериментальное определение сил (нелинейная экстраполяция), X и ю -константы 472 530 578 600 612
Экспериментальное определение сил (линейная экстраполяция), учет зависимостей X (0) и ю (0) 416 456 512 554 588
Экспериментальное определение сил (нелинейная экстраполяция), учет зависимостей X (0) и ю (0)* 496 543 581 596 604
* Окончательный вариант предложенной методики.
При малых величинах Кз предполагаемая невысокая температура на площадке износа (а следовательно, и малая теплопроводность стали) привела бы к снижению стока тепла в заготовку, что является фактором, приводящим к росту величины 0з. Поэтому температура площадки износа оказывается выше. По мере увеличения Кз и соответствующего повышения температуры возрастает и теплопроводность контактных объемов стали, а следовательно и сток тепла в заготовку. В результате рост величины 0з замедляется, что и фиксируется результатами расчетов. В данном случае (в сравнении с линейной экстраполяцией) описанный выше механизм действует не через учет зависимостей X и ю от температуры, а через характер зависимостей изменения сил на площадке износа — ускоренный рост при малых Кз и замедление роста при дальнейшем износе. При этом небольшое отклонение от линейности зависимостей Р2(К) и Р, у (Кз) (т. е. суммарных сил) сопровождается для малых Кз весьма существенным (в процентном отношении) увеличением сил на площадке износа.
В заключение можно отметить, что полученные с помощью предложенной модернизированной методики расчетные данные хорошо согласуются с выявленными ранее результатами при обработке аустенитных сталей и свидетельствуют о ее корректности.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Липатов, А. А. Влияние характера зависимости теплопроводности обрабатываемого материала от температуры на закономерности контактного взаимодействия и тип стружкообразования / А. А. Липатов // СТИН. — 2006. -№ 8. — С. 37−40.
2. Липатов, А. А. Методика определения сил на задней поверхности режущего инструмента / А. А. Липатов, Ю. Л. Чигиринский, С. И. Кормилицын // СТИН. — 2010. -№ 8. — С. 6−8.
3. Никифоров, С. К. Методика отбраковки измерений с аномальными значениями среднеквадратической ошибки / С. К. Никифоров, А. А. Степченко, Н. П. Алаев // Известия ВолгГТУ. Серия «Актуальные проблемы управления, вычислительной техники и информатики в технических системах». Вып. 5: межвуз. сб. науч. ст. / ВолгГТУ. -Волгоград, 2008. — № 8. — С. 45−48.
4. Резников, А. Н. Теплофизика процессов механической обработки материалов / А. Н. Резников. — М.: Машиностроение, 1981. — 279 с.

ПоказатьСвернуть
Заполнить форму текущей работой