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

Автор: Ерисов Ярослав Александрович, Гречников Федор Васильевич, Сурудин Сергей Викторович, Разживин Василий Андреевич

Журнал: Известия Самарского научного центра Российской академии наук @izvestiya-ssc

Рубрика: Машиностроение и машиноведение

Статья в выпуске: 2 т.22, 2020 года.

Бесплатный доступ

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

Еще

Предельная деформация, условие пластичности, кристаллографическая ориентация, модифицированный критерий предельного нагружения

Короткий адрес: https://sciup.org/148314221

IDR: 148314221   |   DOI: 10.37313/1990-5378-2020-22-2-118-123

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

В практике листовой штамповки стандартом для оценки условий начала локализованного образования шейки в процессе формообразования листового металла при различных видах нагружения (от одноосного растяжения, сдвига и плоской деформации до равноосного двухосного растяжения) стали кривые предельных деформаций (FLC) [1], которые впервые были использованы в работах [2,3].

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

экспериментальные методы относительно трудоемкие и дорогостоящие, так как для их проведения необходимо специализированное оборудование, большое количество испытываемого материала [27]. При этом полученные результаты не всегда удается воспроизвести.

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

Кинематические критерии [9], в отличие от энергетических, предсказывают не максимальную несущую способность заготовки, а ее предельное деформированное состояние в момент локализации деформации, когда на большей части заготовки формоизменение останавливается. Наибольшее распространение среди кинематических методов получили методы Г. Свифта и Р. Хилла [10,11]. Оценка, достоверности данных методов показала, что метод Г. Свифта целесообразно использовать для построения кривой предельной деформации только для области, где обе главных деформаций в плоскости листа принимают положительные значения, а метод Р.

Хилла, где одна из деформаций имеет отрицательные значения. Таким образом, совместное применение данных методов позволяет теоретически построить кривую предельных деформаций [12].

В работах З. Марчиняка и К. Кучинского локализация деформаций объясняется тем, что она возникает в начальной геометрической неоднородности (канавке) на поверхности листа [13,14]. В более поздних работах [15-18] модель Марчиняка-Кучинского неоднократно совершенствовалась для того, чтобы учесть различные особенности поведения материала при пластической деформации, например, анизотропию свойств. Однако, основным недостатком всех этих моделей является условие наличия дефекта в виде канавки, размеры и расположение которого значительно влияют на результаты расчетов, кроме того, в большинстве случаев, для их нахождения требуются экспериментальные данные.

В последние годы широкое распространение получил модифицированный критерий предельного нагружения, предложенный П. Хора [19-24], в соответствии с которым достижение максимального усилия в процессе нагружения является необходимым, но не достаточным условием наступления разрушения. При локализации деформации в процессе формообразования происходит изменение деформированного состояния в области шейки на плоское, что приводит к возникновению дополнительных растягивающих напряжений, как раз и являющихся причиной разрушения деформируемой заготовки. Таким образом, модифицированный критерий предельного нагружения учитывает изменение деформированного состояния при локализации деформаций, что позволяет значительно повысить точность теоретического построения кривых предельных деформаций.

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

РАЗРАБОТКА МАТЕМАТИЧЕСКОЙ МОДЕЛИ

Рассмотрим процесс двухосного растяжения прямоугольного листа из ортотропного материала в условиях плоского напряженного состояния ( О з = 0 ). Криволинейные координатные линии u и V выбраны так, чтобы они совпадали с траекториями главных нормальных напряжений О" 1 и о 2 , а также с главными осями анизотропии (рисунок 1).

Рис. 1. Элемент листовой заготовки, находящийся в двухосном напряженном состоянии

Модифицированный критерий предельного

нагружения по П. Хора [19-24] записывается следующим образом:

д о ^ d r i д т 8

~&^ d m j И^ = 1 ,

где т 8 - показатель деформационного состояния:

т 8

где 8 1 и 8 2 - главные нормальные деформации д о  д о дт8

Два слагаемых,---и 8, описывают д81   дт8 д81

приращение напряжений, вызванное упрочнением и изменением деформированного состояния на плоское соответственно. В процессе пластической деформации происходит увеличение напряжений из-за упрочнения материала, одна-д ко, скорость упрочнения —L уменьшается. При д81

д

условии —L ° 1 , отношение деформаций не д 8 1

изменяется - т 8 = const , а значит т 8 = 0 .

_ д

Как только о , то деформированное со- д8 1     1

стояние т 8 изменится. Очевидно, что влияние изменения отношения деформаций будет зави-

сеть от используемого условия пластичности.

Воспользуемся условием пластичности, разработанным в работе [25], которые для случая

плоского напряженного состояния в главных осях анизотропии запишется следующим образом:

° , = ^ V( ^ 12 + П 31 ) - 2 П 12 то + ( ^ 12 + П 23 ) та°"1; (3)

где т ° = О "1 / Г - показатель напряженного состояния; Пу — обобщенные показатели анизотропии:

П = 1 -

15 ( A '- 1 ) 3 + 2 A

А, + А у

-А k

A - параметр анизотропии кристаллической решетки:

S I 111 - SH22

2 S 2323

S jkl — константы кристаллической решетки; A t - ориентационные факторы кристаллографиче

ской ориентировки:

A _ hik2 + k2li +12hi i 2222

( h i + ki + l i )

5 o 1         2 «                             tTj

_--------- ^L n 12 m 8 + (n 12 + 7 23 ) J— . (14)

d 8 1    7 12 7 23 7 31 N                        8 i

Дифференцируя уравнения (10) и (2), распишем второе слагаемое в критерии П. Хора (1): Ю. d т _ 2 °.m 2 (15) 5 m 8 88 1      8 i ( П 12 + П 31 ) m 8 + 2 7 12 m 8 + ( 7 12 + 7 23 )

h i , k i , l i - индексы Миллера, определяющие i -е направление в кристалле относительно системы координат, связанной с листом.

Выражения для интенсивности деформаций и главных деформаций имеют вид:

Окончательно, подставляя выражения (14), (15) и (10) в критерий (1), определим предельную интенсивность деформаций:

кр _ яУ ( 7 12 + 7 31 ) m 2 + 2 7 12 m e + ( 7 12 + 7 23 ) г

+•

V 7 12 7 23 7 31 N

ПиПаП зх Nm j

1 (16)

8 i = V2

V( 7 12 + 7 31 ) - 2 7 12 m g + ( 7 12 + 7 23 ) m o

[ ( 7 12 + П з1 ) m ^ + 2 7 12 m E + ( П 12 + П 23 ) ] [ 7 12 P + ( 7 12 + П 23 ) ]

( 7 12 + П 31 ) - 7 12 m o

8 1 ; (7)

е _ 1           (712 + 731)- 712mo

  • 81    _~r^ I                                        28i

    V2 V( 7 12 + 7 31 ) - 2 7 12 m o + ( П 12 + 7 23 ) mo

_ _ 1           (712 + 723 )mo - 712_

  • 82          I

    V2 V( 7 12 + 7 31 ) - 2 7 12 m o + ( П 12 + 7 23 ) mo

Подставляя (16) в (11) с учетом (2), определим предельные деформации по П. Хора:

е кр _ n +

1        2                     (17)

+______________________ 7 12 7 23 7 31 Nme ______________________,

[( 7 12 + 7 31 ) m 2 + 2 7 12 m e +( 7 12 + 7 23 ) ] [ 7 12 Р +( 7 12 + 7 23 ) ]

екр _ m ee кр .

Подставляя уравнения (8) в соотношение (2) получим уравнения связи между показателями напряженного m o и деформированного т 8 состояний:

т _ ( « 12 + « 23 ) m o - « 12 (9) 8 ( « 12 + « 31 ) - « 12 m o

Расчет т 8 в уравнениях (17) осуществляется по зависимости (9) для заданных значе-

ний та . Выражения (17) позволяют построить кривую предельную деформации на отрезке

m o G [ 0,1 ]

с учетом кристаллографической

ориентации структуры материала.

РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ

Используя зависимость (9), преобразуем выражения (3) и (7) к виду:

4 П 12 П 23 П 31N V( 7 12 + П 31 ) m 8 + 2 7 12 m 8 + ( П 12 + 7 23 ) oi _ ----- г-------------------;---------;----------o 1 ,

V 2                П 12 m 8 + ( П 12 + 7 23 )

8 i _ I            V( 7 12 + 7 31 ) m 8 + 2 7 12 m 8 + ( П 12 + П 23 ) 8 1 .

V 7 12 7 23 7 31 N

(11) Распишем первое слагаемое в критерии

П. Хора (1) следующим образом:

d o i   d o i d o,- д8

---1 _--- 1--- L---L,       (12) 5 8 1    d o i- d8 i 5 8 1

Примем, что упрочнение материала описывается степенной функцией:

o i _ K 8 i ,           (13)

где K и « - константы анизотропного материала.

После дифференцирования выражений для интенсивности напряжений (10), интенсивности деформаций (11) и степенной функции упрочнения (13) первое слагаемое в критерии П. Хора примет вид:

Для листовых материалов характерны определенные преимущественные кристаллографические ориентировки. Так после прокатки сплавов с кубической решеткой преобладают следующие ориентировки деформационного типа: {112}<111> (медь), {110}<112> (латунь), {123}<634> (S) и {100}<011> («куб на ребре»). В отожженных листах преимущественными являются ориентировки рекристаллизационного типа: {100}<001> (кубическая) и {110}<001> (Госса). При этом доля конкретных ориентировок определяется режимами прокатки и термической обработки.

Рассмотрим лист из меди, для которого компоненты тензора податливости S ji равны: 5 1'111 _ 15.0 ТПа-1; S 1122 _- 6.30 ТПа-1 and S 2323 _ 3.33 ТПа-1, тогда по формуле (5): A '_ 3.203 . Для указанных идеальных кристаллографических ориентировок по формулам (4) и (6) рассчитаны обобщенные показатели анизотропии (таблица 1).

На рисунке 2 показаны кривые предельных деформаций для различных кристаллографических ориентировок (таблица 1) при « _ 0,35 . Из рисунка 2 видно, что при т 8 0 ориентировкой, обеспечивающие наивысшие деформационные возможности, является {110}<001>, но при т 8 0

Таблица 1. Ориентационные факторы и обобщенные показатели анизотропии идеальных кристаллографических ориентировок

Кристаллографическая ориентировка

Ориентационные факторы

Обобщенные показатели анизотропии

Наименование

{ hkl }< uvw >

А 1

Д 2

Д 3

П 12

п 23

П 31

Медь

{112}<111>

0,333

0,250

0,250

0,533

1,116

0,533

Латунь

{110}<112>

0,250

0,333

0,250

0,533

0,533

1,116

S

{123}<634>

0,281

0,278

0,250

0,617

0,835

0,814

«Куб на ребре»

{100}<011>

0,250

0,250

0,0

-0,054

1,703

1,703

Кубическая

{100}<001>

0,0

0,0

0,0

1,703

1,703

1,703

Госса

{110}<001>

0,0

0,250

0,250

1,703

-0,054

1,703

Изотропный случай

0,20

0,20

0,20

1,0

1,0

1,0

данная ориентировка показывает самые низкие предельные деформации. Аналогичная картина наблюдается и для некоторых других рассматриваемых ориентировок. Например, ориентировка {100}<011> обеспечивает наивысшие деформационные возможности при m g >  0 , но при m g 0 - все наоборот. Кривые для изотропного случая и кубической ориентировки совпадают.

В отличии от метода Р. Хилла и Г. Свифта [26], метод П. Хора позволяет построить обе ветви кривой предельной деформации, при этом в области m g <  0 учитывается влияние кристаллографической ориентации структуры.

Сравнивая кривые предельных деформаций для различных ориентировок, можно увидеть, что дополнительные напряжения, возникающие вследствие изменения деформированного состояния на плоское, могут приводить как к увеличению (например, ориентировка {100}<011>, правая ветвь кривой), так и уменьшению предельных деформаций ({110}<001>) по сравнению со случаем, когда учитывается только упрочнение. Влияние дополнительных напряжений определяется не только кристаллографией структуры материала, но и деформированным состоянием.

а

б

Рис. 2. Кривые предельных деформаций по методу П. Хора для идеальных кристаллографических ориентировок деформационного (а) и рекристаллизационного (б) типов

ЗАКЛЮЧЕНИЕ

В данной работе были построены кривые предельных деформаций материала на основе критерия П. Хора с учетом условия пластичности, учитывающего кристаллографию структуры материала. Однозначно назвать ориентировку, обеспечивающую наилучшие деформационные возможности невозможно, так как предельные деформации зависят от напряженно – деформированного состояния вызванному, конкретной технологической операцией листовой штамповки. Так при линейном растяжении наилучшие показатели у ориентировки деформационного типа, при двухстороннем растяжении – у ориентировок рекристаллизационного типа. Таким образом, при производстве листовых заготовок необходимо назначить согласованные режимы прокатки и промежуточной термической обработки (отжига), обеспечивающие формирование заданных идеальных кристаллографических ориентировок, необходимых для повышения предельных деформаций в конкретных операциях листовой штамповки.

Исследование выполнено в рамках госзада-ния (тема FSSS-2020-0016).

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

  • Banabic, D. Formability of Metallic Materials. Plastic Anisotropy, Formability Testing, Forming Limits / D. Banabic, H.J. Bunge, K. Pohlandt, A.E. Tekkaya. -Berlin: Springer, 2000. - 335 p.
  • Keeler, S.P. Plastic instability and fracture in sheet stretched over rigid punches / S.P. Keeler, W.A. Backhoffen // ASM Trans. O. - 1964. - Vol. 56. - P. 25-48.
  • Goodwin, G.M. Application of strain analysis to sheet metal forming problems in the press shop / G.M. Goodwin // J. SAE Tech. - 1968. - Art. No. 680093.
  • Чумадин, А.С. Теория и расчеты процессов листовой штамповки (для инженеров) / Чумадин А.С. -Москва: МАТИ, 2014. - 215 с.
  • Storakes, B. Plastic and visco-plastic instability of a thin tube under internal pressure, torsion and axial tension/ B. Storakes // JJVS. - 1968. - Vol. 10, №6. - P. 510-529.
  • Малинин, Н.Н. Устойчивость двухосного пластического растяжения анизотропных листов и цилиндрических оболочек / Н.Н. Малинин // Известия АН СССР. Механика твердого тела. - 1971. - №2. - С. 115-118.
  • Korhonen, A.S. On the theories of sheet metal necking and forming limits / A.S. Korhonen // Journal of Engineering Materials and Technology. - 1978. -Vol. 100. - P. 303-309.
  • Дель, Г.Д. Технологическая механика / Г.Д. Дель. -М.: Машиностроение, 1978. - 174 с.
  • Stoughton, T.B. Review of theoretical models of the strain-based FLD and their relevance to the stress-based / T.B. Stoughton, X. Zhu // International Journal of Plasticity. - 2004. - Vol. 20. - P. 1463-1486.
  • Swift, H.W. Plastic instability under plane stress / H.W. Swift // Journal of the Mechanics and Physics of Solids. - 1952. - Vol. 1. - P. 1-18.
  • Hill, R. On discontinuous plastic states, with special reference to localized necking in thin sheets / R. Hill // Journal of the Mechanics and Physics of Solids. -1952. - Vol. 1. - P. 19-30.
  • Morales-Palma, D. On the Use of Maximum Force Criteria to Predict Localised Necking in Metal Sheets under Stretch-Bending / D. Morales-Palma, A.J. Martínez-Donaire, C. Vallellano // Metals. - 2017. -Vol. 7. - Art. No. 469.
  • Marchiniak, Z. Limit strains in the processes of stretch-forming sheet metal / Z. Marchiniak, K. Kuczynski // Int. J. Mech. Sci. - 1967. - Vol. 9. - P. 609-620.
  • Marchiniak, Z. Influence of the plastic properties of material on the forming limit diagram for sheet metal in tension / Z. Marchiniak, K. Kuczynski, T. Pokora // Int. J. Mech. Sci. - 1973. - Vol. 15. - P. 789-805.
  • Xu, S. Prediction of forming limit curves of sheet metals using Hill's 1993 user-friendly yield criterion of anisotropic materials / S. Xu, K.J.Weinmann // International Journal of Mechanical Sciences. - 1998. - Vol. 40(9). - P. 913-925.
  • Yoshida, K. Path-dependence of the forming limit stresses in a sheet metal / K. Yoshida, T. Kuwabara, M. Kuroda // International Journal of Plasticity. - 2007. -Vol. 23(3). - P. 361-384.
  • Allwood, J.M. Generalised forming limit diagrams showing increased forming limits with non-planar stress states / J.M. Allwood, D.R.Shouler // International Journal of Plasticity. - 2009. - Vol. 25(7). P. 1207-1230.
  • Hu, Q. Prediction of forming limits for anisotropic materials with nonlinear strain paths by an instability approach / O. Hu, L. Zhang, O. Ouyang, X. Li, X. Zhu, J. Chen // International Journal of Plasticity. - 2018. - Vol. 103. - P. 143-167.
  • Hora, P. A prediction method for ductile sheet metal failure using FE-simulation / P. Hora, L. Tong, J. Reissner // Proceedings of NUMISHEET'96, Dearborn, 1996. P. 252-256.
  • Hora, P. Mathematical prediction of FLC using macroscopic instability criteria combined with micro structural crack propagation models / P. Hora, L. Tong, J. Reissner // Proceedings of Plasticity conference, Ouebec, 2003. P. 364-366.
  • Hora, P. Numerical prediction of FLC using the enhanced modified maximum force criterion (eMMFC) / P. Hora, L. Tong // Proceedings of FLC-Zurich'06, 2006. P. 31-36.
  • Krauer, J. Forming limits prediction of metastable materials with temperature and strain induced martensite transformation / J. Krauer, P. Hora, L. Tong // Proceedings of the 9th International Conference on Numerical Methods in Industrial Forming Processes (NUMIFORM 2007), Porto, 2007. P. 1263-1268.
  • Hora, P. Theoretical prediction of the influence of curvature and thickness on the enhanced modified maximum force criterion / P. Hora, L. Tong // Proceedings of the 7th International Conference and Workshop on Numerical Simulation of 3D Sheet Metal Forming Processes (NUMISHEET 2008), Interlaken, 2008. P. 205-210.
  • Hora, P. Modified maximum force criterion, a model for the theoretical prediction of forming limit curves / P. Hora, L. Tong, B. Berisha // International Journal of Material Forming. - 2013. - Vol. 6(2). - P. 267-279.
  • Ерисов Я.А., Гречников Ф.В., Сурудин С.В. Критерий пластичности анизотропной среды с учетом кристаллографии структуры и его экспериментальная проверка // Вестник Магнитогорского государственного технического университета им. Г.Л. Носова. 2016. Т.14. №4. С. 42- 49.
  • Сурудин С.В., Ерисов Я.А., Разживин В.А., Петров И.Н. Построение кривых предельных деформаций по методу Хилла - Свифта с учетом кристаллографии структуры // Известия Тульского государственного университета. Технические науки. - 2019. - № 5. - С. 256-262.
  • ISO 12004-2:2008 Metallic materials - Sheet and strip - Determination of forming-limit curves -Part 2: Determination of forming-limit curves in the laboratory.
Еще
Статья научная