УДК 534.1:674.05
С. С. Г а в р ю ш и н, В. С. Прокопов, М. А. Б л о х и н
РАЗРАБОТКА МЕТОДИКИ ЧИСЛЕННОГО АНАЛИЗА ДИНАМИЧЕСКИХ ХАРАКТЕРИСТИК МНОГОПИЛЬНОГО СТАНКА С КРУГОВЫМ ПОСТУПАТЕЛЬНЫМ ДВИЖЕНИЕМ ДЕРЕВОРЕЖУЩИХ ПОЛОТЕН
Рассмотрена комплексная математическая модель, предназначенная для анализа деревообрабатывающего оборудования принципиально нового класса - многопильных станков с круговым поступательным движением дереворежущих полотен. Проведены численные исследования влияния конструктивных и технологических параметров на функциональные характеристики оборудования. Предложена методика определения диапазона рациональных значений параметров для обеспечения безопасной и эффективной работы, основанная на декомпозиции механической системы с поэтапным расчетом основных узлов станка. Для исследования прочностных и динамических характеристик пильных модуля и блока используется метод конечных элементов, а для анализа движения пильных полотен — качественные методы исследования движения нелинейных систем.
E-mail: [email protected]; [email protected] [email protected]
Ключевые слова: лесопильное оборудование, деревообработка, многопильный станок, коленчатый станок, пильное полотно, динамика, прочность, нелинейные колебания, резонансный режим.
Одной из актуальных задач для перехода к высокотехнологичной инновационной экономике является разработка в области создания деревообрабатывающих станков новых технических решений, сохраняющих преимущества и исключающих недостатки существующего парка оборудования. Для решения этой задачи предложено [1, 2] деревообрабатывающее оборудование принципиально нового класса — многопильный станок с круговым поступательным движением полосовых пильных полотен (рис. 1).
Простота и надежность конструкции предлагаемого станка позволяет обеспечить высокие функциональные характеристики, среди которых следует особо выделить: улучшение качества обработанных поверхностей, снижение энергопотребления; относительно малую массу и динамическую сбалансированность основных узлов; повышенную мобильность оборудования. Особенностью такого оборудования является передача движения с ведущего на ведомый вал непосредственно через пильные полотна, которые выполняют роль гибких связей и включены в состав пильных блоков. В состав блоков также входят эксцентрики, упругие элементы крепления и корпусные детали.
Рис. 1. Принципиальная схема (а) и опытный образец (б) многопильного станка "Шершень М2001" с круговым поступательным движением полосовых пильных полотен
Все попытки создания аналогичного оборудования опытно-эмпирическим путем окончились неудачей [3]. Во время тестовой эксплуатации опытно-промышленных образцов, как правило, возникали режимы работы, сопровождавшиеся поперечным изгибом и скручиванием полотен, что приводило к формированию криволинейных пропилов, излому зубьев, задиру обрабатываемых поверхностей и в ряде случаев к разрушению пильных полотен. Такие режимы наблюдались как при распиловке древесины, так и на холостом ходу. Причины возникновения подобных режимов объясняются сложным динамическим поведением системы, возможностью появления при определенных условиях как обычных, так и параметрических резонансов, а также существенной нелинейностью ее отдельных элементов (конструкция, близкая по схеме к исследуемой, рассматривается в работах Е.К.Зеемана и В.И.Арнольда [4]).
Решение задачи синтеза конструкции невозможно без комплексного исследования динамического поведения системы с использованием приема декомпозиции конструкции и последующего анализа основных узлов и элементов станка в целях отстройки от нежелательных резонансных режимов.
В доступной авторам настоящей статьи литературе в силу новизны конструкции не удалось найти приемлемую методику, позволяющую провести всестороннее исследование данной задачи. Авторская методика предполагает использование приема декомпозиции конструкции, что позволит в простейшем случае свести проблему к двум основным задачам: рассмотрению поведения дереворежущих полотен и исследованию динамики системы двух валов, соединенных гибкими связями. Важно отметить, что дереворежущие полотна находятся в условиях
а б
Рис. 2. Схема (а) и расчетная схема (б) пильного модуля:
1 и 6 — верхняя и нижняя корпусные детали, 2 и 7 — верхний (ведомый) и нижний (ведущий) валы, 3 — эксцентрик, 4 — упругие элементы крепления пильного полотна, 5 — полосовое пильное полотно
предварительного натяжения с эксцентриситетом, а также испытывают воздействие сил резания и сил инерции, возникающих вследствие вращения валов. При анализе поведения пильных полотен используются основные соотношения механики деформируемого твердого тела применительно к тонкостенным конструкциям [5, 6]. Численная реализация решений проводится методом конечных элементов [7-9]. Используемый подход позволяет максимально адекватно описать условия закрепления и нагружения рабочих органов, автоматизировать проведение расчетов и выполнить анализ результатов расчета для различных вариантов исполнения. В соответствии с используемой конструкторской классификацией, система двух вращающихся валов, соединенных гибкими связями, образует пильный блок. Описание движения пильного блока можно свести к исследованию системы нелинейных дифференциальных уравнений, для решения которых используются методы численного интегрирования по времени, спектральный анализ реализаций и качественные методы исследования нелинейных дифференциальных уравнений (отображение Пуанкаре и отображение сдвига) [10-12].
Введем следующие обозначения (рис.2, б): (£) = £ — угол поворота ведущего (нижнего) вала; — угловая скорость вращения ведущего вала, = (£) + Л^ — угол поворота ведомого вала; Ь0 — расстояние между осями валов, Ь' = Ь0 + ЛЬ — расстояние между осями вращения корпусных деталей при повороте ведомого вала относительно ведущего; е — эксцентриситет оси вращения корпусной детали. Запишем выражение для ЛЬ — изменения расстояния между
осями вращения корпусных деталей при взаимном повороте валов:
AL(A<p,<pB(t)) = ((e(cos(^(í) + А() - cos((,(t))))2+
i
+ (Lq + e(sin(^(t) + А() - sin^t))))2) 2 - Lo. (1)
Расстояние h от центра вращения верхнего вала до линии действия упругой силы, возникающей при деформировании приведенного упругого элемента на AL, запишется как
h(A(,( (t)) =
|e(L0 cos(((t) + A() + e sin A()|
(2)
Учитывая геометрические соотношения (1) и (2), уравнение движения пильного модуля для принятой расчетной схемы (рис. 2, б) можно представить в виде
IА ( + KAL(A(, ()h(A(, (в) + aA ( = 0, (3)
где K — приведенная жесткость полотна и упругих элементов; I — момент инерции участка верхнего вала с эксцентриком; a — коэффициент демпфирования.
В случае колебаний с малой амплитудой уравнение движения преобразуется к виду
IА j + Ke2A( cos2 (в + aAj = 0. (4)
Частота малых собственных колебаний данной системы при фиксированном угловом положении ведущего вала (в составляетр0((в) =
Ke2 cos2 (в 0 „ r ( ---. В этой системе есть особые точки ^угловое поло-
п п \
жение ведущего вала, равное +— и —— I, в которых частота соб-
2 2 /
ственных колебаний с малой амплитудой равна нулю, что приводит к необходимости рассмотрения нелинейных уравнений движения без ограничения амплитудных отклонений. На рис. 3 приведены фазовые траектории системы при кинематическом возбуждении колебаний
^(t) = Av cos(w^t),
где A^ — амплитуда и — частота возбуждения.
Проведенные исследования показывают, что в системе (3) могут реализовываться нелинейные режимы движения, а также возможно возникновение бифуркаций, что требует использования более детального анализа в целях выявления всех возможных типов движения.
Рис.3. Вид фазовых траекторий при изменении амплитуды внешнего воздействия (а, б) (стрелкой показана эволюция траекторий при увеличении амплитуды) и частоты внешнего воздействия (в, г) (стрелкой показана эволюция траектории при увеличении )
Вместе с тем, отметим, что система (3) не описывает наличие предварительного натяжения полотен. Учет начального натяжения требует более детальной расчетной схемы. Рассмотрим простейшую расчетную модель, описывающую систему валов с шестью упругими связями (рис. 4).
Уравнение движения пильного блока для принятой расчетной схемы (см. рис. 4) без учета крутильной и изгибной податливости валов примет вид
Др+KДЬ(Др, рв)Л(Др, Рв)+KДЬ (Др, Рв+3) Др, Рв+|) + +KДЬ (Др, Рв) h (Др, Рв) +KДЬ(Др, Рв)Н(Др, Рв) +
+ KДЬ (Др, Рв + 4Р\н(Др, Рв + ^
р 4п \
3 J \ Г
■
+
+ KДЬ (Др, Рв + у МДр, Рв + у) + аДР = 0, (5)
7 (к
-) Л( Дс^в + т
где /^ — суммарный момент инерции верхнего (ведомого) вала с установленными эксцентриками; К — приведенная жесткость полотна и упругих элементов в пильном модуле.
Рис. 4. Упрощенная (а) и расчетная (б) схемы пильного блока с шестью модулями
Рассмотрев малые отклонения углового положения ведомого вала Д^ и проведя тригонометрические преобразования уравнения (5), получим линейное уравнение
1^Дф + 3 Ке2 Д^ + аДф = 0. (6)
Отметим, что при наличии трех и более упругих связей между валами (при условии сохранения эксцентриками взаимного углового положения) уравнение движения становится линейным и перестает зависеть от углового положения ведущего вала. Частота собственных крутильных колебаний ведомого вала как жесткого целого относительно ведущего вала, возникших за счет деформирования упругих
/зке2 ^ „
элементов, равна р0 = \ —^—. Данный результат подтверждается
V 1е
анализом нелинейного уравнения (5). Спектральный анализ реализаций показывает наличие только одной гармоники с частотой р0 в угловом перемещении ведомого вала для широкого спектра управляющих конструктивных параметров и начальных условий.
Следует отметить, что при учете возможного разброса механических характеристик упругих элементов, погрешностей изготовления элементов станка и погрешностей сборки была обнаружена возможность реализации в системе квазипериодических движений в присутствии несоизмеримых частот (рис. 5). Этот факт свидетельствует о том, что при невыполнении определенных требований к качеству деталей возможно появление резонансных режимов в рабочей области.
Рис. 5. Отображение Пуанкаре для квазипериодического движения с несоизмеримыми частотами (а) и спектр Фурье реализации колебательного движения системы (б)
Рис. 6. Конечно-элементная модель пильного полотна, установленного в модуль
Важной задачей является исследование динамических характеристик пильных полотен в зависимости от конструктивных параметров, поскольку известно, что совпадение частот вращения валов с частотами собственных колебаний полотна может привести к возникновению резонансного режима работы. Напомним, что для исключения возможности статической потери устойчивости пильные полотна при установке предварительно растягиваются. Причем для стабилизации режущей кромки полотна нагружаются не "центрально", а с некоторым эксцентриситетом, при этом передняя режущая кромка полотна оказывается в условиях более интенсивного растяжения, чем задняя. Конечно-элементная модель пильного модуля приведена на рис. 6 и состоит из стержневых, оболочечных и специальных конечных элементов (используется программный комплекс АРМ WinMachine). Упругие элементы, обеспечивающие поджатие пильного полотна, моделируются стержневыми конечными элементами третьего порядка, построенными по гипотезе Кирхгофа, имеющими 2 узла и 12 степеней свободы. Пильное полотно и корпусные детали моделируются плоскими обо-лочечными треугольными и четырехугольными элементами второго порядка, построенными по дискретной гипотезе Кирхгофа, имеющими 24 и 18 степеней свободы соответственно. Штифты крепления, подшипники и эксцентрики моделируются специальными конечными элементами типа жесткая вставка, имеющими 2 узла. Рассматривались несколько вариантов исполнения пильных полотен, отличающихся длиной, шириной и толщиной рабочей части.
В результате анализа динамических характеристик пильного модуля были получены зависимости спектра частот собственных коле-
О 909 1817 2726 3635 Н 0 901,4 1802,8 2704,2 3605,6^, Н
в г
Рис. 7. Зависимость спектра частот от силы натяжения при:
свободной длине полотна 652 (а), 465 мм (б... г); шириной 70 (а), 80 (б, в), толщиной 1,4 (а, в); 1,3 (б), 1,2 мм (г)
баний от силы предварительного натяжения для различных вариантов исполнения пильных полотен (рис. 7).
В несбалансированной системе в процессе движения на полотно могут действовать значительные инерционные силы, влияние которых минимизируется посредством специальной балансировки пильного модуля — приведением центров масс деталей, составляющих пильный модуль, к осям вращения корпусных деталей.
Изгибная податливость валов, а также наличие погрешностей изготовления деталей станка и разброс механических характеристик упругих элементов поджатия полотен приводят к различным натяжениям пильных модулей в зависимости от их расположения в пильном блоке. Сравнение зависимостей спектра частот собственных колебаний пильных модуля и блока в сборе показано на рис. 8.
По результатам исследований можно выделить возможные причины возникновения резонансных режимов работы:
• наличие частоты собственных колебаний пильного полотна в рабочем диапазоне частот вращения валов (до 50 Гц или 3000 об/мин);
Рис. 8. Зависимость спектра частот от сил натяжения пильных модуля (а) и блока (б)
Рис. 9. Области допустимых значений частоты вращения валов и сил натяжения полотна и области, в которых возможно возбуждение двух форм колебаний на одной частоте
• наличие частоты собственных колебаний ведомого вала как жесткого целого за счет податливости упругих элементов натяжения в рабочем диапазоне частот вращения валов;
• отсутствие или недостаточная точность балансировки пильного модуля, что может привести к появлению моментов сил инерции, изгибающих полотно в своей плоскости, достаточных для потери устойчивости плоской формы изгиба;
• совпадение частоты собственных колебаний ведомого вала как жесткого целого за счет податливости упругих элементов натяжения с одной из частот собственных колебаний пильного полотна.
По результатам расчета динамических характеристик пильных полотен можно выделить допустимые и нежелательные области изменения сил натяжения и частоты вращения валов для каждого варианта исполнения (рис. 9).
На основе полученных результатов предложена методика выбора рациональных параметров принципиально нового дереворежущего
Рис. 10. Блок-схема методики выбора рациональных параметров многопильного станка
оборудования (рис. 10). Данная методика была использована для синтеза опытно-промышленного образца многопильного станка с круговым поступательным движением полотен "Шершень М2001".
СПИСОК ЛИТЕРАТУРЫ
1. Бл охин М. А. Расчет и проектирование многопильного станка с круговым поступательным движением полосовых пильных полотен: Дисс.... канд. техн. наук. -М., 2005. - 151 с.
2. П а т е н т № 2131806 РФ, 6 В27 В 3/00. Пильный модуль, пильный блок и устройство для распиловки / М.А. Блохин.
3. А.С. 370026 СССР, МКИ В27В3/00. Лесопильная рама.
4. А р н о л ь д В. И. Теория катастроф, 3-е изд., доп. - М.: Наука, 1990. - 128 с.
5. Белкин А. Е., Гаврюшин С. С. Расчет пластин методом конечных элементов: Учеб. пособие. - М.: Изд-во МГТУ им. Н.Э. Баумана, 2008. - 232 с.
6. К о н д р а т ю к А. А., Ш и л ь к о В. К. Оценка напряженного состояния ленточных пил // Изв. Томского политехн. ун-та. 2004. - Т. 307, № 2. - С. 138-142.
7. Бате К., В и л с о н Е. Численные методы анализа и метод конечных элементов: Пер. с англ. - М.: Стройиздат, 1982. - 448 с.
8. Зенкевич О. Метод конечных элементов в технике / Пер. с англ.; Под ред. Б.Е. Победри. - М.: Мир, 1975. - 543 с.
9. Курков С.В. Метод конечных элементов в задачах динамики механизмов и приводов. - СПб.: Политехника, 1991. - 223 с.
10. Б у т е н и н Н. В., Н е й м а р к Ю. И., Ф у ф а е в Н. Л. Введение в теорию нелинейных колебаний. - М.: Наука, 1985. - 385 с.
11. Старжинский В. М. Прикладные методы нелинейных колебаний. - М.: Наука, 1977. - 256 с.
12. Awrej cewicz Jan, Krysko Vadim A. Chaos in structural mechanics // Springer-Verlag Berlin Heidelberg. - 2008. - 426 p.
13. Воскресенский C. А. Резание древесины. - М.-Л. : Гослесбумиздат, 1955. - 190 с.
14. Г а в р ю ш и н С. С., Б а р ы ш н и к о в а О. О., Б о р и с к и н О. Ф. Численные методы в проектировании гибких упругих элементов. - Калуга, 2001. -198 с.
15. Григолюк Э. И., Шалашилин В. И. Проблемы нелинейного деформирования: Метод продолжения решения по параметру в нелинейных задачах механики твердого деформируемого тела. - М.: Наука, 1988. - 232 с.
16. П а т е н т № 3929048 США, В27В 3/00. Многорядная пила с возвратно-поступательным движением полотна.
Статья поступила в редакцию 3.11.2010
Сергей Сергеевич Гаврюшин родился в 1949 г., окончил МВТУ им. Н.Э. Баумана в 1972 г. Д-р техн. наук, профессор, первый проректор - проректор по научной работе МГТУ им. Н.Э. Баумана. Автор более 200 научных работ в области механики деформируемого твердого тела, численных методов анализа инженерных задач, нелинейного поведения тонкостенных конструкций, биомеханики.
S.S. Gavryushin (b. 1949) graduated from the Bauman Moscow Higher Technical School in 1972. D. Sc. (Eng.), professor, first vice-rector on scientific work of the Bauman Moscow State Technical University. Author of more than 200 publications, specialist in the field of mechanics of deformable solid body, numerical methods for analysis of engineering problems, nonlinear behavior of thin-walled constructions, biomechanics.
Михаил Анатольевич Блохин родился в 1949 г., окончил МВТУ им. Н.Э. Баумана в 1974 г. Канд. техн. наук, генеральный директор ООО "Магеллан 2001". Автор шести научных работ и трех патентов в области динамики и прочности.
M.A. Blokhin (b. 1949) graduated from the Bauman Moscow Higher Technical School in 1974. Ph. D. (Eng.), general director of Limited Liability Company "Magellan 2001". Author of 6 publications and 3 patents for inventions, specialist in the field of dynamics and strength, wood processing.
Владимир Сергеевич Прокопов родился в 1980 г., окончил МГТУ им. Н.Э. Баумана в 2004 г. Руководитель направления прочностного анализа ООО НТЦ "АПМ". Автор пяти научных работ в области динамики и прочности, численных методов анализа инженерных задач (МКЭ).
V.S. Prokopov (b. 1980) graduated from the Bauman Moscow State Technical University in 2004. Chief on strength analysis of the Limited Liability Company "NTTs "APM". Author of 5 publications, specialist in the field of dynamics and strength, numerical methods for analysis of engineering problems.
Михаил Иванович Осипов (1938 г.-2010 г.)
Ушел из жизни Михаил Иванович Осипов — замечательный педагог, крупный ученый, вся творческая жизнь которого была связана с Московским государственным техническим университетом им. Н.Э. Баумана и кафедрой "Газотурбинные и нетрадиционные энергоустановки", которой он руководил многие годы.
Михаил Иванович, еще будучи студентом, обратил на себя внимание кафедры и ее руководителя Владимира Васильевича Уварова своей колоссальной работоспособностью, постоянным желанием узнать и освоить возможно больше нового. После успешного окончания в 1963 г. МВТУ им. Н.Э. Баумана он был оставлен на кафедре и активно включился в масштабные научно-исследовательские работы, проводимые в проблемной лаборатории "Газотурбостроение", спектр которых простирался от газотурбинных установок для первого отечественного газотурбовоза, сверхмощных газовых турбин, турбовентиляторных двигателей большой степени двухконтурности — до новых направлений, связанных с использованием МГД-генераторов. В 1973 г. им была защищена кандидатская диссертация. Михаил Иванович прошел на кафедре все ступени роста от инженера-конструктора до профессора, заведующего кафедрой. Одновременно много энергии им было вложено в развитие факультета "Энергомашиностроение": с 1978 г. по 1984 г. он работал заместителем декана, а затем по 1997 г.
был проректором НУК "Э" (деканом факультета "Энергомашиностроение"). С 1998 г. Михаил Иванович бессменно руководил родной ему кафедрой.
Будучи ищущим и инициативным ученым, Михаил Иванович не только всячески содействовал развитию новых направлений энергетики, связанных с использованием нетрадиционных и возобновляемых источников энергии, но и явился инициатором создания такой специальности на кафедре. В связи с этим с 1999 г. кафедра газовых турбин получила новое название — "Газотурбинные и нетрадиционные энергоустановки".
При его участии и под его руководством были выполнены многие исследования и разработки, внесшие большой вклад в развитие отечественной энергетики. Признанием значимости этих работ являются многочисленные награды, которых он был удостоен.
Активная энергия Михаила Ивановича распространялась за пределы нашей страны: слушателями его докладов и лекций были аудитории Великобритании, Чехословакии, Финляндии, Германии, Южной Кореи, Швеции, США, Франции. С 1993 г. Михаил Иванович являлся президентом совета Восточно-Европейского отделения Международной энергетической ассоциации и членом международного жюри по премиям за лучшие работы в области энергетики.
Все кто с ним общался, всегда поражались его эрудиции, его энциклопедическим знаниям. Он был членом ученых советов МГТУ им. Н.Э. Баумана, НУК "Э", совета РАН по прямому преобразованию энергии, Национального комитета РАН по тепломассообмену, членом редакционных коллегий известных научно-технических журналов, таких как "Вестник МГТУ им. Н.Э. Баумана. Серия "Машиностроение", "Проблемы энергетики", международный журнал "Аэромеханика и теплофизика".
Михаил Иванович был замечательным педагогом, любившим студентов, но и требовательным, не прощавшим халтуры. Но с самой высокой шкалой требований он подходил к себе. Его авторитет и среди студентов, и среди коллег был неизменно высок.
Человек высокой культуры, он был простым в общении, чутким и внимательным к окружающим, притягивающим к себе людей. Он умел создавать хороший микроклимат в коллективе.
Память о Михаиле Ивановиче Осипове останется на долгие годы в сердцах всех знавших его.