Математическая модель режимов лазерного резания древесины березы
Знакомство с особенностями построения математической модели глубины лазерного резания березы. Общая характеристика лазерно-гравировального станка VL-4060, анализ основных функций. Рассмотрение уровней и интервалов варьирования переменных факторов.
Рубрика | Математика |
Вид | статья |
Язык | русский |
Дата добавления | 29.11.2018 |
Размер файла | 666,2 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
Размещено на http://www.allbest.ru/
Математическая модель режимов лазерного резания древесины березы
Исследования процессов лазерного резания является актуальной задачей, решение которой позволит инженерам-технологам увеличить производительность лазерных станков путем использования оптимальных режимов резания.
Объект исследования - процесс лазерного резания древесины березы.
Цель исследования - нахождение зависимости глубины резания от мощности лазерного излучателя и скорости перемещения лазерного луча.
Экспериментальная установка включает в себя лазерно-гравировальный станок VL-4060, воздушный компрессор, ПК. ПО LaserCut 5.1, электронный штангенциркуль ЦИ 0-150 с ценой деления 0,01мм, набора щупов ГОСТ 882-75, поверхностный влагомер Hydro CONTROL, весы CAPTO ГОСМ СЕ623-С.
Толщина заготовок древесины березы измерялась электронным штангенциркулем ЦИ 0-150 мм с ценой деления 0,01, с цифровым индикатором «ЭНКОР». Скорость движения технологической головки устанавливалась и контролировалась при помощи ПО LaserCut 5.1, глубина резания измерялась при помощи набора щупов ГОСТ 882-75 и штангенциркуля ЦИ 0-150 мм, с ценой деления 0,01. Заготовки древесины березы представляют собой бруски с размерами 110х20х10 мм в количестве 15 шт. Образцы не имели пороков, были высушены и имели влажность 8-10 %, что проверялось контактным влагомером Hydro CONTROL. После взвешивания на весах весы CAPTO ГОСМ СЕ623-С для эксперимента были выбраны образцы с одинаковой массой и одинаковой плотностью.
Заготовки укладывались на ячеистый стол станка, положение фокальной f плоскости было настроено на верхнюю поверхность (Рисунок 1) и в течение проведения эксперимента не изменялось. Все резы были сделаны поперек волокон древесины при одном направлении движения лазерной головки.
Рисунок 1. Положение фокальной плоскости: 1-ячеистый стол лазерно-гравировального станка, 2- бруски древесины березы
Было проведено 750 резов с различными скоростно-мощностными режимами. Каждый режим был реализован 30 раз. В ходе проведения эксперимента скорость ступенчато менялась и составляла: 10, 12, 18, 41, 90 мм/с. Мощность ЛИ в ходе проведения эксперимента изменялась ступенчато и соответствовала: 21, 24, 30, 35, 38 процентам от максимальной мощности станка, составляющей 65 Вт. Обоснование указанного диапазона мощностей приведено в [5]. Каждому из пяти значений мощности ЛИ соответствовали пять значений скорости.
Измерения глубины реза образцов выполнялись с помощью щупа и штангенциркуля. В полученный рез, глубиной h, помещался щуп высотой Н как показано на рисунке 2 и штангенциркулем измерялось расстояние Sизм от верхней поверхности брусков древесины березы до верхней поверхности щупа, а также высота брусков древесины березы.
Рис.2
Вычисление глубины реза образцов выполнялось с помощью формулы 1.
h=(k+H)-Sизм
При воздействии лазерного излучения на древесный материал глубина резания зависит от параметров: мощности лазерного излучателя, времени воздействия на материал. Формула глубины резания в общем случае известна уже длительное время (Формула 2).
где Р - мощность излучения, Вт; - радиус фокусировки, мм; - скорость резки, мм/с; - плотность вещества, кг/м3; с - теплоемкость, Дж/К; Ти - температура испарения материала, °C; - удельная энергия испарения, кДж/кг.
Однако некоторые коэффициенты, используемые в случае резания древесины, на данный момент не известны. А, именно коэффициенты Ти - температура испарения материала, Lи - удельная энергия испарения.
Поскольку глубина резания зависит обратно пропорционально от скорости резания, что неудобно для построения регрессионной модели. Решили использовать другую величину - время воздействия лазерного луча на материал, которое напрямую зависит от скорости.
Постоянными факторами являются: влажность древесины 8-10 %, масса и плотность образцов.
Переменными факторами являются: время воздействия лазерного луча на материал, которое варьировалось от 2,2 ? t ? 20 мс, мощность лазерного излучения, которое варьировалось от 22 ? Р ? 38 %.
Входными факторами являются: мощность лазерного излучения, время воздействия лазерного луча на материал.
Выходным фактором является глубина резания.
Оценочным показателем является глубина резания. Предполагаемая величина изменения ожидается в диапазоне от 0 ? у ? 5 мм.
Определение верхнего, нижнего и основного уровней факторов в натуральных и нормализованных обозначениях и интервалы их варьирования приведены в таблицах (1) и (2).
математический лазерный древесина станок
Таблица 1 Уровни и интервалы варьирования переменных факторов
Факторы |
Шаг варьирования |
Уровни варьирования факторов |
|||
Нижний |
Основной |
Верхний |
|||
Мощность лазерного излучения (Х1), % |
8,5 |
22 |
29,5 |
38 |
|
1 |
-1 |
0 |
+1 |
||
Время воздействия лазерного луча (Х2), мс |
8,888 |
2,222 |
11,111 |
20 |
|
1 |
-1 |
0 |
+1 |
Таблица 2 Рабочая матрица планирования ПФЭ 2К (по плану Коно)
Номер опыта |
Натуральные значения факторов |
Кодированные значения факторов |
Значение выходной величины |
|||
Х1 |
Х2 |
У |
||||
1 |
24,7 |
20 |
1 |
1 |
У1 |
|
2 |
13,65 |
20 |
-1 |
1 |
У2 |
|
3 |
24,7 |
2,(2) |
1 |
-1 |
У3 |
|
4 |
13,65 |
2,(2) |
-1 |
-1 |
У4 |
|
5 |
24,7 |
11,(1) |
1 |
0 |
У5 |
|
6 |
13,65 |
11,(1) |
-1 |
0 |
У6 |
|
7 |
19,5 |
20 |
0 |
1 |
У7 |
|
8 |
19,5 |
2,(2) |
0 |
-1 |
У8 |
|
9 |
19,5 |
11,(1) |
0 |
0 |
У9 |
По результатам экспериментальных исследований глубины резания в зависимости от мощности и времени воздействия лазерного излучения на древесину березы составили таблицу (3).
Таблица 3 Результаты экспериментальных исследований глубины резания в зависимости от мощности и времени воздействия лазерного излучения
Номер опыта |
Натуральные значения факторов |
Кодированные значения факторов |
Значение выходной величины |
|||
Х1 |
Х2 |
У |
||||
1 |
24,7 |
20 |
1 |
1 |
1,66676 |
|
2 |
13,65 |
20 |
-1 |
1 |
0,45235 |
|
3 |
24,7 |
2,(2) |
1 |
-1 |
0,29176 |
|
4 |
13,65 |
2,(2) |
-1 |
-1 |
0,04735 |
|
5 |
24,7 |
11,(1) |
1 |
0 |
1,15029 |
|
6 |
13,65 |
11,(1) |
-1 |
0 |
0,20971 |
|
7 |
19,5 |
20 |
0 |
1 |
1,35735 |
|
8 |
19,5 |
2,(2) |
0 |
-1 |
0,23324 |
|
9 |
19,5 |
11,(1) |
0 |
0 |
0,95293 |
Для построения математической модели (уравнения регрессии) глубины лазерного резания березы в зависимости от мощности и времени воздействия лазерного излучения определили коэффициенты b0, b1, b2, , b11, b22, b12.
b0 = 0,912;
b1 = 0,4;
b2 = 0, 484;
b11 = -0,211
b22 = -0,96
b12 = 0,243
Методом подстановки получили экспериментально уравнение регрессии по расчету глубины лазерного резания древесины березы на лазерно-гравировальном станке и представлено в виде математической модели в кодированном виде (3):
y=0,912 + 0,4х1 + 0,484х2 - 0,211х12 - 0,96х22 = 0,243х1х2
Как видно из уравнения (3) наибольшее влияние на выходную величину (глубину) оказывает время воздействия лазерного излучения на древесину (фактор х2), имеющую наибольшее значение коэффициента. Это возможно связано с тем, что при увеличении времени усиливаются экзотермические реакции окисления древесины. При увеличении мощности лазерного излучения (х1) значение глубины увеличивается, что и ожидалось.
После провели перевод регрессионного уравнения (3) из кодированного вида в натуральный вид (4).
Соотношение времени воздействия (мс) к скорости перемещения лазерного луча рассчитывали по формуле (5).
где 0,2 мм - диаметр лазерного пятна на материале; 1000 - переводной коэффициент, переводящий в мс; V - скорость перемещения лазерного луча.
В итоге методом упрощения математическая модель приняла следующий вид (6).
где H - глубина реза, мм; P - мощность ЛИ, Вт; V - скорость перемещения лазерного луча, мм/с.Экспериментальная формула 6 справедлива для определения глубины лазерного резания березы с влажностью 6-8% для станка VL4060 внутри интервала варьирования переменных факторов, приведенных в таблице 1.
Литература
1.Григорьянц А.Г., Шиганов И.Н., Мисюров А.И. «Технологические процессы лазерной обработки» - Москва: МГТУ им. Н.Э. Баумана, 2006
2.Григорьянц А.Г. «Основы лазерной обработки материалов» - Москва, 1989
3.Веденов А.А., Гладуш Г.Г. «Физические процессы при лазерной обработке материалов» - Москва: Энергоатомиздат, 1985
4.Либенсон М.Н, Яковлев Е.Б., Шандыбина Г.Д. «Взаимодействие лазерного излучения с веществом. Конспект лекций» - Санкт-Петербург, 2005
5.Кузнецов А.И., Шимон Е.В., Новоселова М.В., Николаева Т.Ю. Определение действительного диапазона изменения мощности лазерного излучения станка vl4060. // Деревообработка: оборудование, менеджмент XXI века. Труды VII международного евразийского симпозиума / Под научной ред. В.Г. Новоселова, 2012.
Размещено на Allbest.ru
...Подобные документы
Знакомство с особенностями построения математических моделей задач линейного программирования. Характеристика проблем составления математической модели двойственной задачи, обзор дополнительных переменных. Рассмотрение основанных функций новых переменных.
задача [656,1 K], добавлен 01.06.2016Проведение аналитического обзора литературы по биологическому действию бетулина. Клинические испытания в отношении В-клеток поджелудочной железы на модели аллоксанового диабета у мышей. Апробация фитопрепарата Ве-Тритерпеноид из коры березы повислой.
статья [19,5 K], добавлен 23.07.2012Проектирование математической модели. Описание игры в крестики-нолики. Модель логической игры на основе булевой алгебры. Цифровые электронные устройства и разработка их математической модели. Игровой пульт, игровой контроллер, строка игрового поля.
курсовая работа [128,6 K], добавлен 28.06.2011Синтез оптимального управления при осуществлении разворота. Разработка математической модели беспилотных летательных аппаратов. Кинематические уравнения движения центра масс. Разработка алгоритма оптимального управления, результаты моделирования.
курсовая работа [775,3 K], добавлен 16.07.2015Словесная, математическая постановка исходной задачи. Исследование математической задачи на корректность. Применение метода экспертных оценок и парных сравнений основных объективных, субъективных факторов, послуживших причиной к поступлению учиться в МАИ.
курсовая работа [145,1 K], добавлен 19.12.2009Динамическая модель как теоретическая конструкция, описывающая изменение состояний объекта. Характеристика основных подходов к построению: оптимизационный, описательный. Рассмотрение способов построения математических моделей дискретных объектов.
контрольная работа [769,7 K], добавлен 31.01.2013Литералы рассуждения и вопрос об их отрицаниях. Математическая модель отрицания для рассуждения, содержащего связную совокупность суждений. Отрицания в математической логике и дополнения в алгебре множеств. Интерпретации формул математической логики.
контрольная работа [40,8 K], добавлен 03.09.2010Знакомство с особенностями возникновения тригонометрии, рассмотрение этапов развития. Анализ способов решения треугольников, основанных на зависимостях между сторонами и углами треугольника. Характеристика аналитической теории тригонометрических функций.
презентация [654,4 K], добавлен 24.06.2014Создание математической модели движения шарика, подброшенного вертикально вверх, от начала падения до удара о землю. Компьютерная реализация математической модели в среде электронных таблиц. Определение влияния изменения скорости на дальность падения.
контрольная работа [1,7 M], добавлен 09.03.2016Предмет, методы и понятия математической статистики, ее взаимосвязь с теорией вероятности. Основные понятия выборочного метода. Характеристика эмпирической функции распределения. Понятие гистограммы, принцип ее построения. Выборочное распределение.
учебное пособие [279,6 K], добавлен 24.04.2009Общая характеристика математической модели радиотехнического сигнала. Значение спектрального разложения функций в радиотехнике. Работа вещественных одномерных детерминированных сигналов и система синусоидальных и косинусоидальных гармонических функций.
курсовая работа [1,0 M], добавлен 13.08.2011Теория игр - математическая теория конфликтных ситуаций. Разработка математической модели игры двух лиц с нулевой суммой, ее реализация в виде программных кодов. Метод решения задачи. Входные и выходные данные. Программа, руководство пользователя.
курсовая работа [318,4 K], добавлен 17.08.2013Понятие функции нескольких переменных. Аргументы, частное значение и область применения функции. Рассмотрение функции двух и трех переменных. Предел функции нескольких переменных, теорема. Главная сущность непрерывности функции нескольких переменных.
реферат [86,3 K], добавлен 30.10.2010Составление математической модели для предприятия, характеризующей выручку предприятия "АВС" в зависимости от капиталовложений (млн. руб.) за последние 10 лет. Расчет поля корреляции, параметров линейной регрессии. Сводная таблица расчетов и вычислений.
курсовая работа [862,4 K], добавлен 06.05.2009Метод интервалов как один из важнейших методов математической деятельности, связанный с вопросами нахождения нулей функции или промежутков ее знак постоянства для неравенства. Алгоритм решения дробно-рационального неравенства методом интервалов.
курсовая работа [630,7 K], добавлен 12.04.2015Метод планирования второго порядка на примере В3-плана. Получение и исследование математической модели объекта в виде полинома второго порядка. Статистический анализ полученного уравнения и построение поверхностей отклика. Расчет коэффициентов регрессии.
курсовая работа [128,4 K], добавлен 18.11.2010Деятельность при решении задач складывается из умственных действий и осуществляется эффективно, если первоначально она происходит на основе внешних действий с предметами. Главная проблема - дети не могут перейти от текста задачи к математической модели.
дипломная работа [79,2 K], добавлен 24.06.2008Свободное падение тела с учетом сопротивления среды. Зависимость перемещения и скорости падения от времени. Формулировка математической модели и ее описание. Описание программы исследования с помощью пакета Simulink. Решение задачи программным путем.
курсовая работа [1,2 M], добавлен 21.03.2011Изучение физического процесса как объекта моделирования. Описание констант и параметров, переменных, используемых в физическом процессе. Схема алгоритма математической модели, обеспечивающая вычисление заданных зависимостей физического процесса.
курсовая работа [434,5 K], добавлен 21.05.2022Решение систем уравнений по правилу Крамера, матричным способом, с использованием метода Гаусса. Графическое решение задачи линейного программирования. Составление математической модели закрытой транспортной задачи, решение задачи средствами Excel.
контрольная работа [551,9 K], добавлен 27.08.2009