Математическая модель динамики процесса вырубки уплотнений из терморасширенного графита

Аннотация


Представлены основные этапы получения прокладок уплотнительных из терморасширенного графита. Предложена методика получения математической модели динамики процесса вырубки этих прокладок на прессовом оборудовании. Процедура получения модели состоит из двух взаимосвязанных этапов. На первом этапе с помощью эксперимента получена зависимость технологической силы процесса вырубки от перемещения рабочего органа, которая имеет ряд особенностей из-за специфических свойств терморасширенного графита. Приведены принципиальная схема эксперимента и его результаты, которые аппроксимировались в виде аналитических зависимостей на основе метода наименьших квадратов. Разработана расчетная схема динамики процесса вырубки, учитывающая динамическую взаимосвязь между двигателем, исполнительным органом в виде кривошипно-ползунного механизма пресса и обрабатываемым материалом. С учетом принятых допущений построена приведенная к кривошипному валу динамическая модель и рассчитаны основные ее геометрические, массовые, а также силовые параметры. Рассмотрены основные этапы функционирования прессового оборудования в процессе вырубки. Основные параметры сведены в соответствующую таблицу, элементы которой использованы при составлении сформированной математической модели. При этом учтены конструктивно-технологические особенности пресса и выделены этапы его работы (рабочий и холостой ход). Структура системы уравнений движения нелинейная и учитывает характеристику и вид технологической силы вырубки, аналитическое значение которой получено на основе проведенных экспериментов. Реализация решения этих уравнений осуществляется с помощью программного модуля DYNAM с использованием метода Рунге - Кутты четвертого порядка. Представлен ряд элементов модуля. В качестве примера приведена графическая иллюстрация одного из результатов расчета. По результатам проведенных исследований сделаны выводы.

Полный текст

Терморасширенный графит - материал с высокой стойкостью к температурным и химическим воздействиям. Обладает гибкостью, податливостью, прочностью, что делает его высокоэффективным набивочным и уплотнительным материалом, широко применяемым в различных областях промышленности [1-4]. В частности, используется для герметизации валов насосов, штоков арматуры, фланцевых соединений и тому подобного в качестве прокладок уплотнительных (ПУТГ). Схема получения ПУТГ представлена на рис. 1, где конечный продукт получается в процессе вырубки на прессовом оборудовании [5-7]. Рис. 1. Последовательность процесса изготовления ПУТГ: ИГ - интеркалированный графит; МЛГ - материал графитовый листовой Следует учесть, что механические свойства ТРГ, вызванные универсальной природой материала и методами его получения, существенно отличаются от свойств металлов. Вследствие этого технологическая сила процесса вырубки имеет ряд особенностей, основная из которых - нелинейность ее изменения, что существенно влияет на динамические процессы изготовления изделия из ПУТГ [1, 3, 4, 8]. Статья посвящена разработке методики определения динамических параметров и характеристик процесса вырубки ПУТГ на прессовом оборудовании. Цель исследований - формирование математической модели динамики процесса вырубки уплотнений из терморасширенного графита, на основе которой возможна оценка нагруженности элементов прессового оборудования. Структура работы заключалась в определении закона изменения технологической силы и изучения с учетом этого динамических явлений, сопровождающих процесс вырубки. Для определения закономерности изменения технологической силы на первом этапе проведены экспериментальные исследования с помощью испытательного комплекса Instron 3369 [3, 9, 10]. Комплекс имеет структуру, функционирующую по замкнутой системе управления (рис. 2). Принципиальная схема эксперимента представлена на рис. 3. Управление испытаниями и формирование результатов эксперимента осуществлялись с помощью программного пакета Instron Bluehill 3 [10]. В качестве испытуемого материала выбраны прокладки ПУТГ из листового графитового картона, армированного перфорированной фольгой из нержавеющей стали. Фрагмент результатов одного из испытаний представлен нелинейной диаграммой «нагрузка F (Н) - перемещение Х (м) инструмента» (рис. 4). Для проведения дальнейших исследований эта диаграмма аппроксимировалась аналитическими зависимостями. В основе приближения использовался метод наименьших квадратов по 60 точкам [11-14]. При аппроксимации комбинацией линейных и параболических зависимостей погрешности достигали величин 38 %. Однако для полинома четвертой степени погрешность оказалась значительно меньшей. В результате расчетов значение силы (H) (1) где Хϵ [0; 3], мм. Фрагмент экспериментальной и теоретической зависимостей технологической силы F представлен на рис. 5. Рис. 2. Структура комплекса Instron Рис. 3. Схема эксперимента: 1 - плата; 2 - штанцформа; 3 - образец МГЛ-2-212; 4 - лист фанеры; 5 - плита Рис. 4. Диаграмма нагружения листа МГЛ Следует отметить, что на отдельных участках максимальные погрешности достигали 8,6 % (рис. 6, а) и 4,5 % (рис. 6, б). Полученные результаты позволили перейти к этапу формирования математической модели динамики процесса вырубки. В качестве объекта Рис. 5. Экспериментальная и теоретическая диаграммы нагружения: - экспериментальные данные; - теоретические данные исследования выбран пресс, исполнительным органом которого является центральный кривошипно-ползунный механизм (КПМ). Его ползун, соединенный с пуансоном, непосредственно производит вырубку [5, 6, 7, 15]. До момента контакта инструмента с заготовкой при отсутствии технологической силы - холостой ход. Внедрение инструмента в заготовку до получения уплотнения - рабочий ход (наличие силы ). Отвод инструмента и возвращение в первоначальное положение - холостой ход. Расчетная схема основана на методиках, изложенных в работах [16-20]. В качестве звена приведения выбран кривошипный вал (рис. 7). При этом приняты следующие допущения: - звенья кривошипно-ползунного механизма абсолютно жесткие; - стол и станина пресса неподвижны; - технологическая сила представлена зависимостью (1) и имеет «скачок» по окончании процесса вырубки; Рис. 6. Сравнение диаграмм нагружения: а - первый участок; б - второй участок; - экспериментальные данные; - теоретические данные - весь процесс вырубки соответствует установившемуся движению КПМ; - связь кривошипного вала и шкива на валу ротора двигателя имитируется упругой связью с коэффициентами упругости СР и сопротивления hР, характеризующими свойства, например, ременной передачи [16, 20]. Расчетная схема процесса вырубки представлена на рис. 7. Рис. 7. Расчетная схема динамики процесса вырубки ПУТГ: а - структурная схема пресса; б - кривошипно-ползунный механизм; в - приведенная расчетная схема Здесь приняты следующие обозначения: 1 - электродвигатель; 2 - шкив ротора двигателя; 3 - упруго-демпферная связь, имитирующая ременную передачу с приведенными коэффициентом жесткости СР и коэффициентом сопротивления hР; 4 - маховик; 5 - кривошипный вал радиусом r; 6 - противовес; 7 - шатун длиной l; 8 - инструмент; Z - толщина заготовки, м; Х - вертикальное перемещение инструмента, м. Согласно схеме (см. рис. 7, б и в), JР и JК - соответственно приведенные моменты инерции ротора двигателя со шкивом на его валу и инерции кривошипа с маховиком, шатуна и ползуна с инструментом, кг·м2; φР и φК - соответственно приведенная обобщенная координата ротора двигателя и кривошипного вала, рад; МПС - приведенный момент сил от весов звеньев КПМ, технологической силы и момента сил трения МТР, Н·м; МД - приведенный момент двигателя, Н·м. Математическая модель для расчетной схемы (см. рис. 7), учитывающая динамическую взаимосвязь между двигателем, передаточным механизмом (ременной передачей) и исполнительным органом, формируется на основе уравнений движения в дифференциальной форме [18, 20, 21]: (2) Ряд основных параметров и вспомогательных коэффициентов системы (2) после преобразований в соответствии с расчетной схемой сведен в таблицу. Значения параметров и вспомогательных коэффициентов № п/п Коэффициент Значение 1 JК 2 МПС 3 МД 4 JР 5 b 6 d 7 k В таблице J0 и JШК - соответственно моменты инерции ротора двигателя и шкива на его валу, кг·м2; J - момент инерции маховика и противовеса, кг·м2; mК, mИ и mШАТ - соответственно массы кривошипа, инструмента и шатуна, кг; М - момент двигателя, Н·м; iД-К - величина передаточного отношения от двигателя на кривошипный вал; GИ и GШАТ - вес инструмента и шатуна, Н; b, d, k - вспомогательные коэффициенты для проведения расчетов. Все преобразования выполнены с учетом [6, 16-19]. Особенность математической модели (2) заключена в ее нелинейности и «кусочном» представлении силы F вида (1). Характер ее воздействия описывается функцией переключения: При этом учитывался характерный для прессового оборудования недоход δ инструмента из-за конструктивных особенностей пресс-формы [5]. Здесь H - ход инструмента, H = 2r; X - перемещение инструмента в соответствии с выбранным направлением, (см. рис. 7, а). Для реализации дальнейших исследований динамики процесса вырубки разработан программный пакет DYNAM в системе MatLab [22, 23]. Основные этапы расчетов: - формирование исходных данных на основе конструктивно-технологических свойств процесса с возможностью их корректировки; - задание начальных условий; - решение задач кинематики; - метод Рунге - Кутты четвертого порядка [24-26]; - уточнение шага интегрирования; - возможность вывода результатов в виде диаграмм и таблиц. В качестве примера на рис. 8-10 приведены результаты одного из вариантов расчета вырубки на прессе КД2322 [15]. При этом толщина заготовки Z = 3 мм, недоход инструмента δ = 1 мм, λ = 0,071, ω = 25 рад/с, Н = 50 мм. Рис. 8. Изменение технологической силы вырубки Рис. 9. Результирующие усилия вырубки Приведенные диаграммы (см. рис. 9, 10) представляют собой зависимость нагружения Q (Н) от перемещения инструмента и угла поворота кривошипного вала при непосредственном процессе вырубки. Резюмируя изложенное, можно сделать следующие выводы: - получена зависимость технологической силы процесса вырубки от перемещения рабочего органа, учитывающая особенности материала из ТРГ; - предложенная математическая модель позволяет учесть взаимосвязь между особенностями обрабатываемого материала, конструктивно-технологическими и силовыми параметрами оборудования в процессе вырубки; - математическое моделирование на основе полученных результатов позволит в дальнейшем осуществить выбор эксплуатационных параметров оборудования, обеспечивающих эффективность процесса вырубки ПУТГ. Рис. 10. Зависимость усилия вырубки от перемещения инструмента и угла поворота кривошипа

Об авторах

А. Е Кобитянский

Пермский национальный исследовательский политехнический университет

В. С Белобородов

Пермский национальный исследовательский политехнический университет

Список литературы

  1. Материалы графитовые листовые (МГЛ) и прокладки из терморасширенного графита (ТРГ) [Электронный ресурс]. - URL: http://yrel77.ru/?page_id=426 (дата обращения: 02.05.2020).
  2. Белова М.Ю., Исаев О.Ю., Смирнов Д.В. Разработка первого межгосударственного стандарта на материалы из терморасширенного графита // Актуальные проблемы порошкового материаловедения: материалы междунар. науч.-техн. конф., посвященной 85-летию со дня рождения академика В.Н. Анциферова, г. Пермь, 26-28 ноября 2018 г. - Пермь, 2018. - С. 391-394.
  3. Кобитянский А.Е., Плюснина М.А. Определение ударной нагрузки при моделировании динамики процесса вырубки уплотнения из терморасширенного графита // Электрофизические методы обработки в современной промышленности: материалы II Междунар. науч.-практ. конф. молодых ученых, аспирантов и студентов, г. Пермь, 18-20 декабря 2018 г. - Пермь, 2019. - С. 250-253.
  4. Караваев Д.М. Композиционные материалы на основе терморасширенного графита для эксплуатации при температурах до 500 °С: дис. … канд. техн. наук: 05.16.06. - Пермь: Изд-во Перм. нац. исслед. политехн. ун-та, 2016. - 16 с.
  5. Бурдуковский В.Г., Инатович Ю.В. Оборудование кузнечно-штамповочных цехов. Кривошипные машины: учеб. пособие. - Екатеринбург: Изд-во Урал. ун-та, 2018. - 168 с.
  6. Кузнечно-штамповочное оборудование: учеб. для вузов / под ред. Л.И. Живова. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2006. - 560 с.
  7. Иванов Ю.В. Научное обоснование и реализация технических решений для улучшения виброакустических характеристик кузнечно-прессовых машин и агрегатов: дис. … д-ра техн. наук: 05.02.09 / ИжГТУ им. М.Т. Калашникова. - Ижевск, 2017. - 329 с.
  8. Механические свойства композиционного материала на основе терморасширенного графита / Д.М. Караваев, А.М. Ханов, А.И. Дегтярев, Л.Е. Макарова, Д.В. Смирнов, О.Ю. Исаев // Известия Самарского научного центра Российской академии наук. - 2012. - Т. 14, № 1 (2). - С. 562-564.
  9. Универсальная электромеханическая испытательная машина instron 3369 [Электронный ресурс]. - URL: http://ckp.tsu.ru/about/equipment/7/21/ (дата обращения: 03.05.2020).
  10. Instron introduces the newest version of bluehill software. - URL: https://www.instron.ru/-/media/literature-library/press-releases/2010/05/bh3-press-release.pdf (accessed 03 May 2020).
  11. Линник Ю.В. Метод наименьших квадратов и основы математико-статистической теории обработки наблюдений. - 2-е изд., доп. и испр. - М.: Гос. изд-во физ.-мат. лит., 1962. - 349 с.
  12. Губанов B.C. Обобщенный метод наименьших квадратов. Теория и применение в астрометрии. - СПб.: Наука, 1997. - 318 с.
  13. Мусатов М.В., Львов А.А. Анализ моделей метода наименьших квадратов и методов получения оценок // Вестник Саратовского государственного технического университета. - 2009. - Т. 4, № 2 (43). - С. 137-140.
  14. Пономарев И.В., Славский В.В. О геометрической интерпретации метода наименьших квадратов // Известия Алтайского государственного университета. - 2012. - № 1-1 (73). - С. 119-121.
  15. КД2322. Пресс однокривошипный простого действия открытый наклоняемый [Электронный ресурс]. - URL: http://stanki-katalog.ru/sprav_kd2322.htm (дата обращения: 24.09.2020).
  16. Теория механизмов и машин: учеб. для вузов / К.В. Фролов [и др.]; под ред. К.В. Фролова. - М.: Высшая школа, 2002. - 496 с.
  17. Механика машин: учеб. пособие для втузов / под ред. проф. Г.А. Смирнова. - М.: Высшая школа, 1996. - 511 с.
  18. Теория механизмов и машин: учеб. пособие / М.З. Коловский, А.Н. Евграфов, Ю.А. Семенов, А.В. Слоущ. - М.: Академия, 2006. - 560 с.
  19. Левитский Н.И. Теория механизмов и машин. - М., 1979. - 576 с.
  20. Зубов В.А., Лисицын В.Д. Динамика механизмов с упругими звеньями. - Л.: Изд-во ЛПИ, 1971. - 168 с.
  21. Коловский М.З. Динамика машин. - Л.: Машиностроение, 1989. - 262 с.
  22. Дьяконов В.П. MATLAB. Полный самоучитель. - М.: ДМК Пресс, 2012. - 768 с.
  23. Ревинская О.Г. Основы программирования в MatLab: учеб. пособие. - СПб.: БХВ-Петербург, 2016. - 208 с.
  24. Демченко В.В. Метод Рунге - Кутты решения задачи Коши для обыкновенных дифференциальных уравнений первого порядка: учеб.-метод. пособие по курсу «Вычислительная математика» / МФТИ. - М., 2004. - 20 с.
  25. Мышенков В.И., Мышенков Е.В. Численные методы. Ч. 2. Численное решение обыкновенных дифференциальных уравнений: учеб. пособие для студ. спец. 073000 / МГУЛ. - М., 2005. - 109 с.
  26. Демидович Б.П., Марон И.А., Шувалова Э.З. Численные методы анализа. Приближение функций, дифференциальные и интегральные уравнения: учеб. пособие. - 5-е изд., стер. - СПб.: Лань, 2010. - 400 с.

Статистика

Просмотры

Аннотация - 91

PDF (Russian) - 32

Ссылки

  • Ссылки не определены.

Данный сайт использует cookie-файлы

Продолжая использовать наш сайт, вы даете согласие на обработку файлов cookie, которые обеспечивают правильную работу сайта.

О куки-файлах