Effect of plastic compressibility on the strain rate intensity factor in compression of a material layer between parallel plates

Abstract


Solutions of boundary value problems for several rigid plastic models may be singular. In particular, some components of the strain rate tensor and the quadratic invariant of the strain rate tensor may approach infinity in the vicinity of certain surfaces. Some models for predicting the high gradient of material properties near frictional surfaces in metal forming processes are based on such behavior of the velocity field. The coefficient of the leading singular term in a series expansion of the quadratic invariant of the strain rate tensor in the vicinity of frictional interfaces is called the strain rate intensity factor. The magnitude of the second invariant of the strain rate tensor is controlled by this factor that depends on geometric parameters of the boundary value problem and parameters involved in the material model. The present paper deals with the effect of plastic compressibility of the material that obeys the double shearing model on the strain rate intensity factor in compression of a plastic layer between two parallel plates. The system of equations comprising the equilibrium equations and constitutive equations is hyperbolic. It is assumed that the surface of the contact between the plates and the deforming material is an envelope of characteristics. An analytic solution of the boundary value problem is found under plane strain deformation. End effects near the free surface of the layer and its center are ignored. The dependence of the strain rate intensity factor on parameters of the boundary value problem including the parameter that controls plastic compressibility is found. In case of the plastically incompressible material, the solution coincides with the available solution.

Full Text

Течение пластического материала вблизи поверхностей трения значительно отличается от течения в основном объеме. В частности, это следует из исследований процессов обработки металлов давлением [1-7]. В связи с таким различием в течении материала вблизи поверхностей трения возникает значительный градиент свойств материала, что приводит к формированию тонкого слоя материала вблизи поверхности трения, свойства которого могут оказывать положительное влияние на эксплуатационные характеристики изделий [8-16]. Аналогичные результаты получены при исследовании технологических процессов со снятием стружки [9-13]. В [14] отмечается, что при моделировании процессов резанья использование вблизи поверхностей трения тех же определяющих уравнений, что и в основном объеме материала, недопустимо. Это заключение справедливо и для моделирования процессов обработки материалов давлением. В [17, 18] для предсказания формирования микроструктуры материала в тонком слое вблизи поверхности трения в процессах обработки материалов давлением предложен подход, основанный на коэффициенте интенсивности скорости деформации. Коэффициент интенсивности скорости деформации для модели идеально жесткопластического материала введен в [19]. Этот коэффициент контролирует величину эквивалентной скорости пластической деформации вблизи поверхности максимального трения. Результат, полученный в [19], обобщен на некоторые модели вязкопластических материалов в [20, 21], на модель двойного сдвига в [22] и на модель двойного скольжения и вращения в [23]. Описание модели двойного сдвига, которая используется в настоящей статье, представлено в [13]. Коэффициент интенсивности скорости деформации зависит от параметров краевой задачи и модели материала. Обзор решений для коэффициента интенсивности скорости деформации содержится в [25]. Все эти решения получены для пластически несжимаемых материалов. Целью настоящего исследования является определение влияния пластической сжимаемости на величину коэффициента интенсивности скорости деформации на примере осадки полосы между двумя параллельными плитами в условиях плоскодеформированного состояния. Принимается модель двойного сдвига, учитывающая пластическую сжимаемость материала [26]. Постановка задачи Рассмотрим слой материала, сжимаемый двумя параллельными плитами в условиях плоскодеформированного состояния (рисунок). Материал подчиняется модели двойного сдвига [26]. Текущую толщину слоя обозначим , текущую ширину . Предполагается, что . Скорость движения каждой плиты равна U. Решение строится для мгновенного напряженно-деформированного состояния. Краевая задача имеет горизонтальную и вертикальную оси симметрии. Оси декартовой системы координат совпадают с этими осями симметрии. В связи с этим достаточно рассмотреть первый квадрант. Приближенное решение, построенное ниже, является обобщением решения Прандтля для сжатия слоя идеально жесткопластического материала между параллельными плитами [27]. В частности, краевые условия при и не могут быть удовлетворены точно. В решении Прандтля эти условия удовлетворяются в интегральном смысле. В рассматриваемом случае такой подход также возможен. Однако эти условия не оказывают никакого влияния на окончательный результат исследования и поэтому не включены в постановку краевой задачи. Остальные краевые условия выполняются точно. Рис. Геометрическая схема процесса Fig. Geometry of the process Компоненты тензора напряжения в декартовой системе координат обозначим , и . Одно из краевых условий для напряжений следует из симметрии краевой задачи. В частности, (1) при . Второе краевое условия для напряжений должно быть сформулировано при . Известно, что уравнения для напряжений при применении модели двойного сдвига являются гиперболическими [26]. Примем, что поверхность совпадает с огибающей семейства характеристик. Проекции вектора скорости на оси декартовой системы координат обозначим и . Одно из краевых условий для скоростей следует из симметрии краевой задачи. В частности, (2) при . Второе краевое условие для скоростей следует из закона движения плит: (3) при . Уравнения модели двойного сдвига включают условие пластичности Кулона-Мора и закон пластического течения. В условиях плоскодеформированного состояния эти уравнения, учитывающие пластическую сжимаемость материала, приведены в [26]. Условие пластичности имеет вид (4) Здесь p - линейный инвариант тензора напряжения; - квадратичный инвариант тензора напряжения; - угол внутреннего трения и k - коэффициент сцепления. Отметим, что в [26] предполагается, что . Общий вид условия пластичности Кулона-Мора используется в [13]. Пусть - угол между осью x и максимальным главным напряжением, отсчитываемый от оси против хода часовой стрелки. Тогда компоненты тензора напряжения выражаются как (5) Направление течения материала в первом квадранте (см. рисунок) показывает, что . В этом случае из (5) следует, что (6) Уравнения закона пластического течения имеют вид (7) Здесь обозначает полную производную по времени, а величина характеризует пластическую сжимаемость материала. Пластически несжимаемый материал получается при . В общем случае величина определяется историей деформирования. Однако в публикуемой работе, следуя [26], предполагается, что - постоянная величина. При и уравнения (7) сводятся к уравнениям ассоциированного закона течения. Уравнения (4), (5) и (7) дополняются уравнениями равновесия в форме (8) Система уравнений (4), (5) и (8) является гиперболической. Уравнения характеристик этой системы имеют вид [13, 26] (9) Таким образом, огибающая характеристик совпадает с поверхностью , если или при . Принимая во внимание (6), краевое условие для напряжений на поверхности плиты записываем в форме (10) при , где . Общее решение По аналогии с решениями аналогичных краевых задач [26-30] полагаем, что и не зависят от . Величина p исключается с помощью (4). В частности, (11) Тогда уравнения (5) принимают вид (12) Подставляя (12) в (8) и учитывая, что не зависит от , найдем (13) Исключим в этих уравнениях производную . В результате получим (14) Так как не зависит от х, то общее решение этого уравнения имеет вид (15) Здесь - произвольная функция своего аргумента. Подставим (15) в первое уравнение системы (13). В результате получим (16) Это уравнение может быть представлено в форме (17) Левая сторона этого уравнения является функцией z, а правая - функцией y. Поэтому уравнение (17) может иметь решение, только если (18) где C - постоянная величина. Общее решение уравнения (18) имеет вид , где - постоянная интегрирования. С учетом (15) это решение может быть представлено как (19) Из (17) и (18) получаем уравнение для в виде (20) Общее решение этого уравнения имеет форму (21) Здесь - постоянная интегрирования. Из (1), (5) и (6) следует, что при . Подставляя это краевое условие в (21), найдем, что . Тогда решение (21) преобразуется к виду (22) Для определения B воспользуемся краевым условием (10). Подставляя его в (21), найдем (23) При вычислении интеграла, входящего в (19), удобно перейти от интегрирования по y к интегрированию по с помощью (20). Тогда, принимая во внимание (15), (24) Здесь необходимо исключить B с использованием (23). Постоянная не может быть определена без введения краевого условия для напряжений в интегральном виде. Зависимость компонент тензора напряжения от декартовых координат вычисляется из (12) при помощи (22) и (24) в параметрическом виде. Перейдем к определению поля скорости. По аналогии с решениями аналогичных краевых задач [26-30] примем, что (25) Тогда, краевые условия (2) и (3) выполняются. Подставляя (25) во второе уравнение системы (7), получим . Интегрируя это уравнение, найдем (26) Здесь функция интегрирования не зависит от x. Подставляя (25) и (26) в первое уравнение системы (7) и используя (20) для исключения производной , получим (27) Вычислим производную . Так как не зависит от x и , то (28) Из (22) следует, что (29) Тогда (30) Подставляя (20), (25) и (30) в (28), можно убедиться, что . В этом случае уравнение (27) принимает вид (31) Интегрируя и исключая B с помощью (23), найдем (32) Здесь - постоянная интегрирования. Подставляя (32) в (26), получим (33) Уравнения (22), (25) и (33) определяют зависимость проекций вектора скорости от декартовых координат в параметрическом виде. Постоянная не может быть найдена без введения краевого условия для скоростей в интегральном виде. Коэффициент интенсивности скорости деформации Коэффициент интенсивности скорости деформации введен в [19] для случая идеально жесткопластического материала и является коэффициентом при главном сингулярном члене в разложении второго инварианта тензора скорости деформации в ряд в окрестности огибающей семейства характеристик. Это разложение имеет вид (34) при . Здесь определения для второго инварианта тензора скорости деформации выписано для случая плоской деформации; , и - компоненты тензора скорости деформации; D - коэффициент интенсивности скорости деформации и s - расстояние по нормали от огибающей. Очевидно, что при . Поскольку и ограничены, то поведение второго инварианта тензора скорости деформации вблизи огибающей полностью контролируется сдвиговой скоростью деформации. Из (20), (23), (25) и (33) находим (35) Из (10) следует, что при . Таким образом, при . Вблизи точки сдвиговая скорость деформации представляется выражением (36) при . Уравнение (20) вблизи огибающей (или ) имеет вид (37) при . Решение этого уравнения, удовлетворяющее краевому условию при , имеет форму (38) при . Подставляя (38) в (36) и исключая B при помощи (23), найдем (39) при . В рассматриваемом случае . Поэтому, сравнивая (34) и (39), получаем выражение для коэффициента интенсивности скорости деформации в виде (40) В случае идеально жесткопластического материала ( и ) это выражение совпадает с полученным в [25]. Выводы Объемная пластическая деформация влияет на коэффициент интенсивности скорости деформации. Из (40) и (41) следует, что если материал разрыхляется при пластической деформации ( ), то коэффициент интенсивности скорости деформации меньше в предположении о пластической несжимаемости материала. Если материал уплотняется ( ), то объемная пластическая деформация уменьшает коэффициент интенсивности скорости деформации. Таким образом, в соответствии с теорией [17, 18], величина влияет на формирование свойств материала в тонком слое вблизи поверхностей трения.

About the authors

E. A Lyamina

Ishlinsky Institute for Problems in Mechanics RAS; Samara National Research University

O. V Novozhilova

Bauman Moscow State Technical University

References

  1. Griffiths B.J. Mechanisms of white layer generation with reference to machining and deformation processes // Trans. ASME J. Trib. - 1987. - Vol. 109 - P. 525-530. doi: 10.1115/1.3261495
  2. Трунина Т.А., Коковихин Е.А. Формирование мелкодисперсной структуры в поверхностных слоях стали при комбинированной обработке с применением гидропрессования // Пробл. машиностр. и надежн. машин. - 2008. - № 2. - C. 71-74.
  3. Alexandrov S., Jeng Y.-R., Hwang Y.-M. Generation of a fine grain layer in the vicinity of frictional interfaces in direct extrusion of AZ31 alloy // ASME J. Manuf. Sci. Eng. - 2015. - Vol. 137, no. 5 - Paper 051003. doi: 10.1115/1.4030267
  4. Generation of a layer of severe plastic deformation near friction surfaces in upsetting of steel specimens / S. Alexandrov, L. Sidjanin, D. Vilotic, D. Movrin, L. Lang // Metals. - 2018. - Vol. 8. - Paper 71. doi: 10.3390/met8010071
  5. Sanabria V., Mueller S., Reimers W. Microstructure evolution of friction boundary layer during extrusion of AA 6060 // Procedia Eng. - 2014. - Vol. 81. - P. 586-591. doi: 10.1016/j.proeng.2014.10.044
  6. Effect of frictional conditions on the generation of fine grain layers in drawing of thin steel wires / A. Stolyarov, M. Polyakova, G.Atangulova, S. Alexandrov, L. Lang // Metals. - 2019. - Vol. 9, no. 8. - Article 819. doi: 10.3390/met9080819
  7. Effect of die angle and frictional conditions on fine grain layer generation in multi-pass drawing of high carbon steel wire / A. Stolyarov, M. Polyakova, G. Atangulova, S. Alexandrov // Metals. - 2020. - Vol. 10. - Article 1462. doi: 10.3390/met10111462
  8. Griffiths B.J., Furze D.C. Tribological advantages of white layers produced by machining // ASME Journal of Tribology. - 1987. - Vol. 109, no. 2. - P. 338-342. doi: 10.1115/1.3261363
  9. Umbrello D., Jawahir I.S. Numerical modeling of the influence of process parameters and workpiece hardness o white layer formation in AISI 52100 steel // Int. J. Adv. Manuf. Technol. - 2009. - Vol. 44. - P. 955-968.
  10. Investigations of white layer formation during machining of powder metallurgical Ni-based ME 16 superalloy / S.C. Veldhuis, G.K.Dosbaeva, A. Elfizy, G.S. Fox-Rabinovich, T. Wagg // J. Mater. Eng. Perform. - 2010. - Vol. 19. - P. 1031-1036. doi: 10.1007/s11665-009-9567-7
  11. Modeling of white and dark layer formation in hard machining of AISI 52100 bearing steel / D. Umbrello, A.D. Jayal, S. Caruso, O.W. Dillon, I.S. Jawahir // Mach. Sci. Technol. - 2010. - Vol. 14. - P. 128-147. doi: 10.1080/10910340903586525
  12. Cho D.-H., Lee S.-A., Lee Y.-Z. Mechanical properties and wear behavior of the white layer // Tribol. Lett. - 2012. - Vol. 45. - P. 123-129. doi: 10.1007/s11249-011-9869-4
  13. Experimental research material characteristics effect on white layers formation in grinding of hardened steel / X. Huang, Z. Zhou, Y. Ren, C. Mao, W. Li // Int. J. Adv. Manuf. Technol. - 2013. - Vol. 66. - P. 1555-1561. doi: 10.1007/s00170-012-4439-y
  14. Jaspers S.P.F.C., Dautzenberg J. Material behaviour in conditions similar to metal cutting: Flow stress in the primary shear zone // J. Mater. Process. Technol. - 2002. - Vol. 122. - P. 322-330. doi: 10.1016/S0924-0136(01)01228-6
  15. Kim Y.-T., Ikeda K. Flow behavior of the billet surface layer in porthole die extrusion of aluminum // Metallurg. Mater. Trans. - 2000. - Vol. 31A. - P. 1635-1643. doi: 10.1007/s11661-000-0173-4
  16. Choi Y. Influence of a white layer on the performance of hard machined surfaces in rolling contact // Proc. Inst. Mech. Eng., Part B. - 2010. - Vol. 224, no. 8. - P. 1207-1215. doi: 10.1243/09544054JEM1847
  17. Гольдштейн Р.В., Александров С.Е. Подход к предсказанию формирования микроструктуры материала вблизи поверхностей трения при развитых пластических деформациях // Физическая мезомеханика. - 2014. - Т. 17, № 5. - C. 15-20. doi: 10.24411/1683-805X-2014-00046
  18. Александров С.Е., Гольдштейн Р.В. К построению определяющих уравнений в тонком слое материала вблизи поверхностей трения в процессах обработки материалов давлением // Докл. РАН. - 2015. - Т. 460, № 3. - С. 283-285. doi: 10.7868/S0869565215030093
  19. Alexandrov S., Richmond O. Singular plastic flow fields near surfaces of maximum friction stress // Int. J. Non-Linear Mech. - 2001. - Vol. 36, no. 1. - P. 1-11. doi: 10.1016/S0020-7462(99)00075-X
  20. Alexandrov S., Mustafa Y. Quasi-static axially symmetric viscoplastic flows near very rough walls // Appl. Math. Model. - 2015. - Vol. 39, no. 15. - P. 4599-4606. doi: 10.1016/j.apm.2014.12.052
  21. Alexandrov S., Mustafa Y. Singular solutions in viscoplasticity under plane strain conditions // Meccanica. - 2013. - Vol. 48, no. 9. - P. 2203-2208. doi: 10.1007/s11012-013-9736-5
  22. Александров С.Е., Лямина Е.А. Сингулярные решения при плоском пластическом течении материалов, чувствительных к среднему напряжению // Докл. РАН. - 2002. - Т. 383, № 4. - С. 492-495.
  23. Wang Y., Alexandrov S., Lyamina E. Solution behavior in the vicinity of characteristic envelopes for the double slip and rotation model // Appl. Sci. - 2020. - Vol. 10. - Article 3220. doi: 10.3390/app10093220
  24. Spencer A.J.M. A theory of the kinematics of ideal soils under plane strain conditions // J. Mech. Phys. Solids - 1964. - Vol. 12. - P. 337-351. doi: 10.1016/0022-5096(64)90029-8
  25. Alexandrov S. The strain rate intensity factor and its applications: a review // Mater. Sci. Forum. - 2009. - Vol. 623. - P. 1-20. doi: 10.4028/www.scientific.net/MSF.623.1
  26. Spencer A.J.M. Compression and shear of a layer of granular material // J. Eng. Math. - 2005. - Vol. 52. - P. 251-264. doi: 10.1007/s10665-004-5662-9
  27. Хилл Р. Математическая теория пластичности. - М.: Гостехиздат, 1956. - 408 c.
  28. Marshall E.A. The compression of a slab of ideal soil between rough plates // Acta Mech. - 1967. - Vol. 3. - P. 82-92. doi: 10.1007/BF01453708
  29. Collins I.F., Meguid S.A. On the influence of hardening and anisotropy on the plane-strain compression of thin metal strip// Trans ASME J. Appl. Mech. - 1977. - Vol. 44. - P. 271-278. doi: 10.1115/1.3424037
  30. Nepershin R.I. Non-isothermal plane plastic flow of a thin layer compressed by flat rigid dies // Int. J. Mech. Sciences. - 1997. - Vol. 39, no. 8. - P. 899-912. doi: 10.1016/S0020-7403(96)00095-1

Statistics

Views

Abstract - 222

PDF (Russian) - 114

Cited-By


PlumX


Copyright (c) 2021 Lyamina E.A., Novozhilova O.V.

Creative Commons License
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies