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

Применение осесимметричной жесткопластической безмоментной конечноэлементной модели для определения коэффициентов трения в процессах формоизменения листовых металлов Текст научной статьи по специальности «Механика и машиностроение»

CC BY
125
27
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ОСЕСИММЕТРИЧНАЯ ЖЕСТКОПЛАСТИЧЕСКАЯ БЕЗМОМЕНТНАЯ КОНЕЧНОЭЛЕМЕНТНАЯ МОДЕЛЬ / КОЭФФИЦИЕНТЫ ТРЕНИЯ / ФОРМОВКА СФЕРИЧЕСКИМ ПУАНСОНОМ / AXISYMMETRIC RIGID-PLASTIC FINITE ELEMENT MEMBRANE MODEL / FRICTION COEFFICIENTS / HEMISPHERICAL PUNCH STRETCHING

Аннотация научной статьи по механике и машиностроению, автор научной работы — Петров В. К., Михайлова В. Л., Сухомлинов Л. Г.

Излагаются результаты расчетно-экспериментальных исследований по применению осесимметричной жесткопластической безмоментной конечноэлементной модели для определения коэффициентов трения в процессах листовой штамповки. Используются данные по формовке исследуемого листового металла сферическим пуансоном. Возможности методики демонстрируются на примерах с различными условиями трения.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по механике и машиностроению , автор научной работы — Петров В. К., Михайлова В. Л., Сухомлинов Л. Г.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Application of an axisymmetric rigid-plastic finite element membrane model for evaluation of friction coefficients in sheet metal forming processes

The results of computational and experimental investigations are presented on the application of an axisymmetric rigid-plastic finite element membrane model for evaluation of friction coefficients in sheet metal forming processes. The data on hemispherical punch stretching of a sheet metal under investigation are used. The method's possibilities are demonstrated by examining cases with different friction conditions.

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

Engineering Materials, № 6, 2004, pp. 448 - 451 (на английском)

8. http://www.ifam-dd.fraunhofer.de (по состоянию на 26.03.2012)

9. Andersen О., Stephani G., Handbook of cellular metals. Solid - state and deposition methods, ed. by H. P. Degischer and B. Kriszt, Wiley-VCH Verlag GmbH, Weinheim, 2002, pp. 56 - 70 (на английском)

10. Torobin L., Method for making hollow porous microsphers, US Patent № 4671909 (на английском)

11.Ericson R. Syntactic metals: flight weight materials, Advanced materials & Processes, December, 2002, pp. 44 - 47 (на английском)

12. www.recemat.com/en/index.html (по состоянию на 26.03.2012)

13.www.msm.cam.ac.uk/mmc/people/old/dave/ (по состоянию на 26.03.2012)

14. Kendall J.M., Lee М.С., Wang T.G. Metal shell technology based upon hollow jet instability, Journal of Vacuum Science and Technology, vol. 20, April, 1982, pp. 1091 - 1093 (на английском)

15. Lee M.C. Metal shell technology. SAMPLE Journal, November/December, 1983, pp. 7 - 11 (на английском)

16. Lee M., Kendall M., Wang T. Investigation of metallic and metallic glass hollow spheres for fusion target application, Materials processing in the reduced gravity environment of space, 1982, pp. 105 - 113 (на английском)

17. Дороготовцев В., Меркульев Ю. Методы изготовления полых микросфер-микробаллонов, Препринт// Физический институт им. Лебедева, 1989, стр. 57

18. Truong V., Takakura Н., Wells J., Minemoto Т. Production of hollow spheres of eutectic tin-lead solder through a coaxial nozzle, Journal of solid mechanics and materials engineering, Vol. 4, No. 10, 2010 (на английском)

19.Bulina N.V., Gromyko A.I., Bondarenko G.V., Marachevsky A.V., Chekanova L.A., Prokofjev D.E., Churilov G.N. The physics of metals and Metallography, Vol. 102, Suppl. 1, 2006, pp. S94-S95 (на английском)

20. Eitel W., Physikalische Chemie der Silikate. 1941, S. 826 (на немецком)

21. Torobin L.B. Method and apparatus for producing hollow metal microspheres and microspher-oids, US Patent № 4415512 (на английском)

22.Frosch R.A. Method and apparatus for producing concentric hollow spheres, US Patent № 4279632 (на английском)

23.Beggs J.M. Method and apparatus for producing gas-filled hollow spheres, US Patent № 4344787 (на английском)

24. Hendricks C.D. Method for producing small hollow spheres, US Patent № 4133854 (на английском)

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

к.т.н. доц. Петров В.К., к.т.н. доц. Михайлова В.Л., д.т.н. проф. Сухомлинов Л.Г.

Университет машиностроения 8(495)223-05-23, доб. 1318

Аннотация. Излагаются результаты расчетно-экспериментальных исследований по применению осесимметричной жесткопластической безмоментной конечноэлементной модели для определения коэффициентов трения в процессах листовой штамповки. Используются данные по формовке исследуемого листового металла сферическим пуансоном. Возможности методики демонстрируются на примерах с различными условиями трения.

Ключевые слова: осесимметричной жесткопластическая безмоментная ко-

Раздел 2. Технология машиностроения и материалы. нечноэлементная модель, коэффициенты трения, формовка сферическим пуансоном

Эффективность применения математических моделей для исследования процессов листовой штамповки существенно зависит от достоверности принимаемых за основу исходных данных, в число которых входят параметры, характеризующие пластические свойства соответствующего листового металла, а также коэффициенты трения. Наибольшие трудности при этом связаны с достоверной оценкой коэффициентов трения. Среди используемых способов определения коэффициентов трения в зоне контакта металлического листа с инструментом можно отметить подходы, связанные с расчетно-экспериментальным исследованием обтяжки металлической полосы по цилиндрическому инструменту [1], обтяжки цилиндрической оболочки по бочкообразному инструменту [2], формовки листовой заготовки сферическим пуансоном [3]. Можно также указать на специальные устройства для измерения коэффициентов трения в операциях листовой штамповки (см., например, статью [4]). В работах [5,6] определение коэффициентов трения производится расчетно-экспериментальным путем с применением конечноэлементных моделей, которые затем могут быть использованы и для исследования операций формообразования с найденными коэффициентами трения. В настоящей статье расчетно-экспериментальное определение коэффициентов трения основывается, как и в [6], на использовании безмоментной жесткопластической конечноэлементной модели [7,8]. При этом возможности принятой методики демонстрируются на примерах с существенно различающимися условиями трения между формуемым металлическим листом и инструментом.

Основные положения используемой при исследованиях вычислительной модели состоят в следующем.

Исходим из предположения, что формуемая из листового металла осесимметричная оболочка относится к классу тонких безмоментных оболочек. Упругими деформациями на фоне больших пластических деформаций пренебрегаем, считая материал оболочки жестко-пластическим. Используем предложенный Р. Хиллом вариант теории течения для трансвер-сально изотропного материала с изотропным упрочнением. Считаем, что взаимодействие оболочки с инструментом осуществляется в соответствии с кулоновским законом трения. Меридиан срединной поверхности рассматриваемой оболочки в ее исходном недеформиро-ванном состоянии разбиваем на такое количество N участков малых размеров, чтобы в течение всего процесса деформирования допустимо было бы пренебрегать их кривизной, считая эти участки прямолинейными. Процесс формоизменения подобной безмоментной оболочеч-ной модели, состоящей из указанных N элементарных оболочек с прямолинейными образующими, рассматриваем как пошаговый, при котором переход из известного состояния в момент времени I в новое состояние, относящееся к моменту времени 1 + 1^.1, осуществляется с малыми приращениями деформаций. На указанном малом временном интервале At (шаге нагружения) формулировку задачи для принятой дискретной модели оболочки выполняем в терминах узловых перемещений (перемещений концов отрезков разбиения упомянутого меридиана оболочки) с использованием цилиндрической системы координат (X, Г, ф).

Обозначая через х, Г и X , Г координаты точек срединной поверхности дискретной модели рассматриваемой оболочки в начале и конце текущего шага нагружения, а через их, иг осевые и радиальные перемещения указанных точек на данном шаге, мы можем записать:

х*=х + их, (хог). (1)

Координаты и перемещения, относящиеся к узлам дискретной модели, будем обозначать как Х1., г{, Х1. , Г{ , ы'х, и'г (/ = 1,2,..., А^ + 1). С учетом этого, принимая во внимание связи (1), для значений я и длины прямолинейной образующей / -ого элемента (/ = 1,2,...,Л^) в начале и конце шага нагружения можно записать:

[х1+1-х1У + {г1+1 -гг)\

I- /- (2)

**=- ху+(С - гу = ^(х1+1 - х,+и1;1 - иу+(г1+1 - гг + и1;1 - иу .

Здесь и далее индекс / у величин, относящихся к / -ому элементу, для простоты опускается.

Приращение деформации (малое относительное удлинение) элемента модели на шаге нагружения в меридиональном направлении (8^ ) определяем с учетом (2) по схеме:

я* я я*2 s2 1

я 2-я 2-я (3)

+2 • (гм -гг)- (и;+1 - и].) + (и*;1 - и'У + К+1 - и'У ]

Отмечаем, что выражение (3) является квадратично нелинейным по отношению к узловым перемещениям на шаге нагружения. Оно справедливо для любых значений этих перемещений при условии, что | « 1.

Малое относительное удлинение на шаге нагружения окружного материального волокна в середине элемента (приращение деформации е в окружном направлении в середине

элемента) определяем с учетом (1) по схеме:

еФ = (г* -г)/г = йг/г, (4)

г = (гг+ гм)/2, г* = (г* + г*+1)/2, йг = (и\ + и'г+1)/2. (5)

Приращение деформации 8Л по толщине элемента на шаге нагружения определяем исходя из условия несжимаемости, а именно:

Ч=- (е,+еф). (6)

Основываясь на принятом варианте теории течения, соотношения между устанавливающимися в элементах модели на шаге нагружения напряжениями и соответствующими приращениями деформаций записываем в виде:

2-(Я+ 2) 3-(1 + 2-Д)'ДЙ

у/Я + 2 Г, „ ч2 ч2 „ , ч2п1/2 (7)

„ ; •—■-[(1+д)-е,+д-Еф], <=>Ф),

Ле = л —

3 (1 + 2-Я)

о = Ф(е + Ае),

где: и С7ф - меридиональное и окружное напряжения, СТ - эквивалентное напряжение,

Дё - приращение накопленной эквивалентной деформации,

8и8 + Де = 8 - величины накопленной эквивалентной деформации в начале и конце шага нагружения,

Я и С7 = Ф(е) - экспериментально определяемые параметр нормальной анизотропии и

функция упрочнения листового металла. Считаем, что все нагрузки, действующие на рассматриваемую дискретную модель формуемой оболочки, сведены к обобщенным узловым силам 2л • и 271-Рг' (/= 1,2,...,Л/Г +1) в радиальном и окружном направлениях. Исходим из принципа возможных перемещений, утверждающего, что в состоянии равновесия работа приложенных к деформируемому телу сил на вариациях возможных перемещений равна работе напряжений на

соответствующих вариациях деформаций. В результате для рассматриваемой модели оболочки на шаге нагружения получаем вариационное уравнение вида:

ЛГ+1

2тг • • 8е, + аф • 5вф) -гЬ = 2тг • ■ 5< + ^' • Ъи\.). (8)

1=\

Суммирование в левой части уравнения (8) осуществляется по всем элементам модели.

С использованием связей (3)-(7) вариационное уравнение (8) формулируется исключительно в терминах узловых перемещений. Приравнивая коэффициенты при одноименных вариациях узловых перемещений в левой и правой части такого уравнения, получаем для / -ого узла дискретной модели оболочки на шаге нагружения уравнения равновесия вида:

/,;(ад ) = /•;.

Ц{и1,и1) = К\ (9)

0' = /-1, /, / + 1),

где: Ьх, Ц. - нелинейные алгебраические операторы, отражающие геометрическую и физическую нелинейность рассматриваемой задачи на шаге нагружения.

В случае закрепленного / -ого узла относящаяся к этому узлу пара силовых уравнений (9) заменяется кинематическими условиями вида:

и'= О, и/=0. (10)

Если 7 -ый узел находится в свободном состоянии (не контактирует с инструментом), в соответствующих ему уравнениях (9) следует положить:

^'=0, ^'=0. (11)

Для узла, находящегося в контакте с инструментом, контур которого на шаге нагружения описывается уравнением ^(х, г) = 0, должно выполняться кинематическое условие ви-

4>(х1+и'х,г1+и'г) = 0 (12)

и силовое условие (соответствующее кулоновскому закону трения) вида:

ы

к

А и'

(13)

где: |И - коэффициент трения,

Аиге - относительное перемещение узла вдоль поверхности инструмента на шаге нагружения,

' и I'],' - тангенциальная и нормальная по отношению к контуру инструмента компоненты вектора узловой силы. Обозначая как ех, в'г и пх, п'г компоненты векторов касательной и нормали к контуру

инструмента в текущем положении 7-ого узла, силы Ре\ Рп\ входящие в уравнение (13),

можно выразить через силовые факторы Рх', I']' в виде:

(14)

т-" г г т—. г 7—>г ^ '

к< =П • Г +П • Г .

п X X г г

Использование записи (14) предполагает, что входящие в нее факторы /'_',', I']' посредством связей (9) явно выражены через искомые узловые перемещения на шаге нагружения.

Что касается перемещения Аи'е, входящего в формулировку уравнения (13), то для него имеет место зависимость вида:

(15)

где: А11 - заданное перемещение инструмента в направлении оси х на шаге нагружения.

Замечаем, что физические соотношения (7) и силовые соотношения (13) имеют особен-

ности при Дв = О

и

Аи1

0. в целях устранения этих особенностей предусмотрена регу-

ляризация указанных соотношений путем введения в пренебрежимо малых окрестностях

значении Дв = 0

и

Аи\

= 0 соответствующих фиктивных линейных относительно Дв и

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Аи\

участков. Подобный прием [7, 8] позволяет напряжениям , С7ф и силам при не-

обходимости снижаться вплоть до нулевых значений, обеспечивая выполнение условий равновесия на шаге нагружения.

Решение физически и геометрически нелинейной контактной задачи для описанной дискретной модели оболочки на шаге нагружения сводится посредством итерационной процедуры к решению последовательности линейных задач. При этом линеаризация исходной нелинейной системы уравнений на шаге нагружения в рамках такой процедуры осуществляется с использованием методов Ньютона и переменных параметров. Итерационные уточнения продолжаются до достижения заданной относительной точности (80Т) по перемещениям. Решение соответствующей системы линейных алгебраических уравнений выполняется по методу Гаусса.

Перейдем теперь к заявленному исследованию с применением описанной модели.

В качестве базового эксперимента при оценке коэффициента трения между формуемым металлическим листом и инструментом использовалась формовка круглой заготовки (из данного материала) сферическим пуансоном (рисунок 1).

Д.

-?

А .А, к и

Р I и

Рисунок 1 - Схема формовки сферическим пуансоном

Конкретные исследования выполнялись применительно к образцам из листовой стали 08КП толщиной И = 0,75 мм . Из испытаний на одноосное растяжение была установлена кривая упрочнения данного материала в виде С7 = А ■ б" (где: =541 МПа, и =0,22) и определено значение коэффициента нормальной анизотропии в виде Я = 1,26. Эксперименты по формовке производились с использованием пуансона диаметром Пп = 60 мм и матрицы, диаметр отверстия которой составлял величину Пм = 72 мм, а радиус скругления рабочей кромки -величину гм - 3 мм . Фланец круглой заготовки из исследуемого материала удерживался от перемещений в процессе формовки прижимным кольцом с рифтом диаметром 1)р =110 мм .

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

Раздел 2. Технология машиностроения и материалы, рошок канифоли. Отметим, что соответствующие этим условиям экспериментальные данные на обсуждаемых ниже графиках изображаются символами •, ■, А.

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

При проведении расчетов было принято N = 220 , AU = 0,5Dn /3600 и S0T = 0,001.

Сплошные кривые на рисунке 2, помеченные цифрами 1, 2, 3, представляют собой полученные расчетным путем зависимости нагрузки Р на пуансоне от перемещения U пуансона в процессе формовки данного листового металла при значениях коэффициента трения |Ы = 0,14; 0,3; 0,5 (соответственно). Видно, что на начальном этапе формоизменения (10 мм <U <20 мм) величина коэффициента трения практически не влияет на силовую характеристику процесса формовки. Это подтверждают и представленные здесь же экспериментальные результаты, относящиеся к существенно разным условиям трения на поверхности пуансона. Однако момент разрыва (значение перемещения U*, при котором сила Р достигает своего предельного значения) существенным образом зависит от значения коэффициента трения |Ll.

Р, кН 30 ■ 25 ■ 20 15 ■ 10 ■ 5 ■ 0

10 15 20 25 U, мм

Рисунок 2 - Графики зависимости силы, приложенной к пуансону, от его перемещения при различных значениях коэффициента трения

Кроме указанных трех кривых, расчетным путем были получены также аналогичные силовые характеристики для целого ряда других значений коэффициента трения JLL, и для каждого из таких значений определены деформации оболочки в момент достижения пре-

\ 2

J

У

Раздел 2. Технология машиностроения и материалы, дельного состояния.

Рисунок 3 - График зависимости предельного значения меридиональной деформации в полюсе формуемой оболочки от коэффициента трения

На рисунке 3 представлена полученная расчетным путем зависимость меридиональной логарифмической деформации (еДО)) в полюсе формуемой оболочки в предельном состоянии от коэффициента трения (LX на поверхности пуансона (сплошная кривая).

Здесь же представлены экспериментальные результаты (полученные с доведением формуемой оболочки до разрыва), относящиеся к трем упоминавшимся случаям покрытия поверхности пуансона. Сопоставление указанных расчетных и экспериментальных результатов позволяет заключить, что значения коэффициента трения при покрытии поверхности пуансона полиэтиленовой пленкой, минеральным маслом и порошком канифоли могут быть оценены величинами 0,14; 0,3 и 0,5 (соответственно). Представленное на рисунках 4, 5, 6 сопоставление расчетных и (обозначенных точками) экспериментальных результатов по полной картине распределения деформаций (на момент разрыва оболочки) в зоне контакта оболочки с пуансоном применительно к упомянутым трем случаям покрытия его поверхности подтверждает такую оценку. Здесь 8ф - логарифмические деформации оболочки в окружном направлении. Цифрами 1, 2, 3 помечены расчетные результаты, полученные при значениях перемещения С/[мм] пуансона вида (27,8; 28,4; 28,8) в случае JLX = 0,14 (рисунок 4); вида (20,0; 20,4; 20,8) в случае JLL = 0,3 (рисунок 5) и вида (16,0; 16,4; 16,6) в случае JLL = 0,5 (рисунок 6).

3

/V2

1,2,3

[1 1 3 2 3 3 3 43

Рисунок 4 - Распределение деформаций вдоль радиуса заготовки к моменту разрыва формуемой оболочки в случае нанесения полиэтиленовой пленки на поверхность

пуансона

Из указанных рисунков видно, что результаты, помеченные цифрой 1, соответствуют моменту достижения формуемой оболочкой предельного состояния, поскольку дальнейшее незначительное продвижение пуансона (варианты 2 и 3) приводит к локализации деформации. Об этом можно судить по наблюдаемому катастрофическому росту пиков кривых распределения деформаций е3 с одновременным прекращением роста деформаций в зоне кон-

такта формуемой оболочки с пуансоном.

0J 0:6 0:5 ОД 0,3 0,2 0:1 0 "0,1

3

X"2

[ 1 1, 2, 3 0 2 ] 3 0 £<Р 4Ь

г, мм

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

0.5 ОД 0,3 0,2 0,1 о

■0.1 -I

3

2

L

■ 1, 2, 3

II 1 : 2 ] 3 3 43

1.UU

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

В качестве общего вывода по изложенной статье отметим, что представленные результаты расчетно-экспериментальных исследований позволили продемонстрировать возможности применения осесимметричной жесткопластической безмоментной конечноэлементной модели для определения коэффициентов трения в процессах формоизменения листовых металлов. Описанная методика включает испытания с доведением исследуемого листового металла до разрыва в процессе формовки сферическим пуансоном и сопоставление замеренных деформаций в зоне полюса отформованной оболочки с тем, что дает вычислительная модель формуемой оболочки в момент достижения предельного состояния.

Литература

1. Pearce R. Some effects of friction in punch-stretching// Dev. Draw. Metals. Proc. Int. Conf. / London. 11-13 May. 1983. / London. 1983. P. 249-254.

2. Одинг С.С., Максименков В.И. Определение коэффициента трения при формировании осесимметричной оболочки. //Известия вузов. Машиностроение. 1984. № 8. с. 23-27.

3. Ведмедь Ю.П., Юдович С.З., Фишман И.М. Определение коэффициента трения между пуансоном и заготовкой при листовой штамповке. // Кузнечно-штамповочное производство. 1972. №4. с. 25-26.

4. Fratini L., Lo Casto S., Lo Valvo E. A technical note on an experimental device to measure friction coefficient in sheet metal forming. // Journal of Materials Processing Technology. 2006. V. 172. № l.P. 16-21.

5. Wang N.-M., Budiansky R. Analysis of sheet metal stamping by a finite-element method. // Trans. ASME. J. Appl. Mech. 1978. V. 45. № 1. P. 73-82.

6. Матвеев А.Д., Сухомлинов Л.Г., Жиляев С.Д., Соловей Е.М. Расчетно-экспериментальная методика определения коэффициента трения для предельного формоизменения при листовой формовке. // Вопросы исследования прочности деталей машин. Сборник научных

трудов кафедры «Прикладная механика» МИП. М. 1994. с. 104-107.

7. Сухомлинов Л.Г., Энгельсберг В.К. Конечноэлементная система автоматизированного расчета напряженно-деформированного состояния тонких оболочек в процессах осесим-метричного формоизменения под действием жестких штампов. // Известия вузов. Машиностроение. 1989. № 3. с. 66-71.

8. Sukhomlinov L.G., Engelsberg V.K., Davydov V.N. A finite element membrane model for the analysis of axisymmetric sheet metal forming processes. // Int. J. Mech. Sei. 1992. V. 34. № 3. P. 179-193.

Оценка безотказности токарных станков с ЧПУ по параметрам точности обработки

к.т.н. доц. Авдеев В.Б.

Университет машиностроения 8 (495) 223-05-23, assi(aimami.ru

Аннотация. В статье рассмотрен метод определения точностной надёжности

токарных станков с ЧПУ с учётом взаимного влияния процессов, снижающих точность обработки.

Ключевые слова: металлообрабатываюгцне станки, параметрическая надежность, безотказность.

Для металлорежущих станков наиболее характерны параметрические отказы. Существующие методы оценки параметрической надёжности экспериментальным путём можно разбить на две группы: статистические методы и методы, исследующие формирование отказа на основе анализа физической картины его возникновения. Трудность применения статистических методов связана с необходимостью иметь большое количество объектов испытаний и длительностью испытаний. В связи с этим распространение получила вторая группа методов, общая методология которых разработана проф. A.C. Прониковым [1]. Применение данного подхода к конкретном типу оборудования требует дополнительных разработок и исследований с учётом конструкции и условий эксплуатации станков этого типа.

Важным при оценке надежности станка является выбор параметра, по которому будет оцениваться уровень надежности. Наиболее полно характеризовать качество станка может точность обрабатываемых деталей. Выбор этого критерия качества обусловливает комплексный подход к оценке надежности, под которым понимается, во-первых, рассмотрение технологической системы как единого целого, т.е. учёт влияния на точность деталей не только станка, но и приспособления, инструмента, детали; во-вторых, рассмотрение одновременного воздействия на формирование точности всех процессов и явлений, возникающих при работе станка. Комплексность метода требует и позволяет учесть взаимное влияние процессов, протекающих в различных элементах технологической системы.

Целью данной работы является разработка метода определения безотказности станков токарной группы с ЧПУ по параметрам точности обрабатываемых деталей с учетом взаимного влияния ухудшающих точность процессов.

Для достижения поставленной цели необходимо было решить следующие задачи:

а) определить, как происходит формирование точности обработки на станке данной конструкции и как точность изменяется во времени;

б) составить математическую модель, описывающую процесс изменения точности во времени;

в) проверить степень адекватности модели и реального процесса обработки деталей на станке;

г) произвести оценку надежности станка по выбранному параметру.

В качестве объекта исследования взят станок с ЧПУ токарной группы среднего типоразмера. Поскольку исследование станка во всём диапазоне технологических возможностей

i Надоели баннеры? Вы всегда можете отключить рекламу.