Управление деформацией сдвига и ее роль в формировании напряженного и деформированного состояний (НДС) Текст научной статьи по специальности «Физика»

УДК 621.771 (0.75.8)

DOI: 10.14529/теМ 70312

УПРАВЛЕНИЕ ДЕФОРМАЦИЕЙ СДВИГА И ЕЕ РОЛЬ В ФОРМИРОВАНИИ НАПРЯЖЕННОГО И ДЕФОРМИРОВАННОГО СОСТОЯНИЙ (НДС)

Н.В. Судаков

г. Челябинск, Россия

На основе единой математической модели обобщены варианты осадки прямоугольной заготовки с различными схемами действия напряжений контактного трения. Методом минимизации невязок при решении переопределенной системы уравнений определены поля напряжений и скоростей. Показана возможность эффективного влияния на НДС путем целенаправленного управления схемой действия контактного трения, что может оказаться полезным при снижении усилия деформации, формировании структуры материала, уплотнении некомпактных материалов и при решении каких-либо других специальных задач.

Впервые показано влияние крутящего момента от сил трения при несимметричной осадке на характер напряженного и деформированного состояний в объеме заготовки. При этом, в частности, установлено наличие зон затрудненной деформации и их смещение к краям заготовки по мере увеличения степени несимметрии схемы действия внешнего трения. Показано, что влияние крутящего момента на НДС растет с увеличением коэффициента трения и отношения продольных размеров заготовки к ее высотным размерам.

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

Сделаны рекомендации по использованию метода минимизации невязок при решении переопределенной системы уравнений в задачах по анализу НДС при обработке материалов давлением. Проиллюстрирована его эффективность, сочетающая простоту с достаточно высокой точностью по сравнению с вариационным методом.

Ключевые слова: осадка; деформация; напряжения; трение; момент крутящий; сила; симметричная; несимметричная; эпюра; минимизация; невязка; переопределенная; система; уравнения.

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

обоснования их эффективности уделяется недостаточное внимание. Настоящая работа ставит своей целью в определенной мере дополнить имеющуюся информацию по затронутой проблеме.

В рамках единой математической модели покажем влияние схемы действия напряжений трения на НДС в объеме тела. Поставленную задачу решим на примере сжатия заготовки из идеально пластичного материала в условиях плоской деформации при различных схемах действия напряжений трения (рис. 1). Напряжения трения представим линейными зависимостями от координаты х:

х1 ( х) := ЦХ

Т9 ( X

(х) := цх

у + ь (У- 1) у-х (у-1)

(1)

где ц - коэффициент трения; х^ - сопротивление материала пластической деформации

У

V

О

рр^п^щщ

X

а)

тх(х)

у-

+ ь

У

б)

гДх)

-у-

А ъ

г2(х)

У

V

V

в)

Рис. 1. Некоторые из возможных вариантов эпюр трения при осадке заготовки прямоугольного сечения: а - симметричная осадка (у = 0); б - несимметричная осадка (у = 0,5); в - предельный

случай несимметричной осадки (у = 1)

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

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

В основу разработки математической модели положен метод минимизации невязок при решении переопределенной системы уравнений, предложенный в работе [2].

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

из замкнутой системы уравнений МСС, включая граничные условия, условие несжимаемости и другие возможные ограничения.

Касательные напряжения аху зададим

степенным полиномом, ограничившись членами четвертой степени

I 2 . I 2 . 4 .

а ^ = ax + a2 x + a3 xy + a4 y + a5 x +

3 2 2 3 4

+ a6 xy + a7 x y + a8 x y + a9 y ,

(2)

где а1 - параметры, подлежащие определению.

При симметричной и несимметричной осадке должны выполняться условия а х = а |_х

'+ y

и а

..+ x И - y

= а |_ - условия асимметрии. Усло-

&1+У

вия асимметрии будет выполняться, если суммарная по х и у степень отдельных членов ряда

2 з

(2) будет четной (ху, х , ху и т. д.). По этой причине из ряда (2) исключены члены с нечетной суммарной по х и у степенью. Удовлетворив граничные условия для аху :

аху|у = И _Т1(х) = 0 и аху|у=_И _т2(х) = 0 (3)

нашли а2, а3, а4, а5, а7, а8.При этом использовали метод неопределенных коэффициентов (группировка членов уравнения по степеням независимой переменной и приравнивание нулю выражений при разных степенях независимой переменной). В результате получили

а^(х^ ^ ^ a9). (4)

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

Эа

Эx Эа

xy

xy

dy Эа

+ -

yy

=0;

= 0

(5)

(6)

дх ду

определили нормальные напряжения ахх и ауу с точностью до неизвестных функций интегрирования f (у) и f (х), которые представили отрезками степенных рядов:

напряжений ауу , действующих на контакте с

инструментом; Mх - крутящий момент от напряжений трения.

В итоге задача свелась к определению восьми констант (a1, a6, a9, a10, a11, a1з, a14, a15) при наличии неудовлетворенного условия пластичности:

аXX(х, у, a1, a6, a9, a10, °11) -

ayy(x? y, a6, a!3? a14? a15)

- 2т _ 1 -■

аxy (x' y5 a6' a9)2

= 0.

(11)

f(y)=a10+ an y2 + a12 y4;

\ 2 4

f (x) = a13 + a14x + a15x .

(7)

(8)

Константу a7 нашли из условия равенства нулю интегральной суммы продольных напряжений на границах заготовки (х = Ь и х = - Ь ):

h

¡ахх (Ь, y, ^ a6, a9, a10, a11, ^2)Ф = 0. (9)

- h

Константу a12 определили из условия равенства нулю суммы моментов внешних сил, действующих на заготовку,

о

Myy = i aHy=h -xdx;

-b

b

Mт = i x1(x) • hdx,

(10)

-b

где Myy - крутящий момент от нормальных

Для удобства дальнейших выкладок представим уравнение (11) следующим образом: П( x, y, a1, a6, a9, aю, an, a^, a14, a15) = 0. (12) В соответствии с методом, предложенным в работе [2], тело заготовки разбили равномерной прямоугольной сеткой, включающей 25 узлов с конкретными координатами, для которых записали уравнение (12). В результате решалась система из 25 уравнений при восьми неизвестных параметрах. Переопределенная система уравнений приближенно решалась относительно параметров ai с помощью программы, минимизирующей невязки. Такую программу, в частности, реализует встроенная в Mathcad функция Minerr. Ниже предлагается вниманию блок программы, ответственный за минимизацию невязок при решении переопределенной системы уравнений (рис. 2), а также фрагменты результатов

Рис. 2. Фрагмент программы, минимизирующей невязки

т

расчета, выполненного для условий: к = 50 мм; Ь = 100 мм; ц = 0,4; х, = 300 МПа (рис. 3-6).

При этом левый ряд рисунков характеризует НДС при симметричной осадке (у = 0), а правый - при несимметричной осадке (у = 1). Для других возможных вариантов несимметричной осадки НДС будет иметь промежуточный характер между представленными на рисунках крайними случаями. На рис. 3 показаны эпюры относительных нормальных напряжений (ауу / а^) при симметричной осадке ( у = 0 ) и несимметричных вариантах осадки. Полученные результаты подтверждают имеющиеся представления о симметричной осадке. Что касается несимметричных вариантов, то в литературе отсутствует информация о полях напряжений, корректно учитывающих специфику процессов. В частности, в известных решениях не учитывается крутящий момент, обусловленный внешними силами трения. Необходимость уравновешивания крутящего момента приводит к несимметричному распределению нормальных напряжений на противоположных контактных поверхностях. Анализ показал, что с увеличением отношения Ь / к влияние крутящего момента на асимметрию напряженного состояния ослабевает. Следует также заметить, что относительное среднее усилие деформации (Рс / ах) при у = 0, у = 0,5 и у = 1 для указанных выше условий расчета соответственно составило: 1,26; 1,18 и 1,02, а при Ь / к = 8: 1,75; 1,44 и 1,08. Эти результаты также не противоречат имеющимся представлениям.

На рис. 4 показаны карты линий уровня относительного значения среднего нормального напряжения (Ка = а / ах); (здесь а - среднее

нормальное напряжение, а а^ = л/3 • х^). Имеет место изменение неоднородности напряженного состояния по мере увеличения степени не симметрии (увеличение у от нуля до единицы). Так, при у = 0 максимальное отклонение Ка от среднего значения составило 22 %, а при у = 1 - 79 %. Таким образом, для условий расчета (Ь / к = 2) неравномерность напряженного состояния оказывается не в пользу несимметричной осадки. При увеличении отношения Ь / к картина меняется в пользу несимметричной осадки, например, при Ь / к = 10 соответственно имеем 56 и 9,5 %. Это можно объяснить

тем, что с увеличением Ь / к увеличивается плечо крутящего момента нормальных напряжений по сравнению с плечом напряжений трения, что и способствует более эффективному уравновешиванию момента от напряжений трения. Заметим, что математические модели, не учитывающие крутящие моменты, всегда приводили к более однородному напряженному состоянию при несимметричной осадке. Таким образом, влияние схемы действия сил трения (деформации сдвига) на качественную картину напряженного состояния достаточно очевидно.

Поле скоростей при известном напряженном состоянии определили из физических уравнений [3]:

Sx

H

%

1 %

- + С;

с уу + с '

(13)

oxy

С = Т.-,

. H

где ^ - компоненты тензора скоростей деформаций сдвига, связанные со скоростями течения &г- (&х, &у) известными геометрическими уравнениями ^ = 0,5(&г- ^ + &1); Н -интенсивность скоростей деформаций сдвига

(Н = 2^& + ^ ).

После определения ^ убедились, что условие несжимаемости (Е)хх + уу = 0) выполняется автоматически. Функцию интегрирования f (х) при решении геометрического

уравнения относительно скорости &у(х, у) нашли, удовлетворив граничное условие & у (х, к) = -¥ (здесь V- скорость инструмента). Условие &у (х, -к) = V выполнилось автоматически. Функцию интегрирования f (у) при решении геометрического уравнения относительно &х (х, у) задали отрезком степенного ряда с неизвестными коэффициентами, которые определялись из равенства ^ - = 0,

т. е. из равенства ^, полученной из физического уравнения, значению , полученной

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

С xx = Т

Рис. 3. Эпюры нормальных напряжений на контакте с инструментом

шт(Ко) = -0.828 Ка тт(Ко) = -0.924

теап(Ка) = -0.706 теап(Ка) =-0.516

тая(Ка) = -0.55 тая(Ка) = -0.243

Рис. 4. Карты линий уровня средних нормальных напряжений

Рис. 5. Карты линий уровня интенсивностей скоростей деформации сдвига

На рис. 5 представлены карты линий уровня относительного значения интенсивности скоростей деформаций сдвига Н/Н0 (здесь Н0 - интенсивность скоростей деформаций сдвига при однородной деформации). Темные участки на рис. 5 показывают наличие зон

затрудненной деформации (интенсивность скоростей деформаций сдвига меньше, чем при однородной деформации). По мере увеличения степени не симметрии (увеличение у) области затрудненной деформации смещаются к краям заготовки (см. рис. 5), а не исчезают

Рис. 6. Векторные диаграммы поля скоростей

вообще, как это показывают решения, не учитывающие наличие крутящих моментов.

Векторные диаграммы полей скоростей на рис. 6 иллюстрируют характер течения металла при симметричной и несимметричной осадке.

Заключение

1. На основе единой математической модели обобщены различные варианты несимметричной, включающей симметричную как частный случай, осадки заготовки прямоугольного сечения из идеально пластичного материала для условий плоской деформации.

2. Методом минимизации невязки при решении переопределенной системы уравнений [2] определены поля напряжений и скоростей для возможных вариантов несимметричной осадки. Максимальная ошибка при приближенном удовлетворении условия пластичности составила при у = 0 порядка ±0,2%, при у = 0,5 - ±5%, а при у = 1 -±10%. Можно предположить, что величина ошибки зависит от числа членов ряда уравнения (2), степени переопределения системы уравнений и каких-либо других факторов. Этот вопрос требует изучения.

Ранее подобной сложности задача могла быть решена только вариационным методом [3], существенно более трудоемким в реализации. Автор выражает надежду, что метод минимизации невязки при решении переопре-

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

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

Литература

1. Выдрин, В.Н. Принципиальные и теоретические основы нового процесса «прокатка-волочение» / В.Н. Выдрин, Л.М. Агеев // Теория и технология прокатки. - Челябинск: ЧПИ, 1971. - № 76. - С. 3-21.

2. Судаков, Н.В. Дискретно-локальный метод минимизации невязок при приближенном решении уравнений в задачах обработки металлов давлением / Н.В. Судаков // Наука ЮУрГУ: материалы 66-й научной конференции Секции технических наук. - 2014. - С. 16691675. - кир:/^расе. 8ти.ги/хт1ш/ЫШгеат/ каМк/0001.74/4356/18^/.

3. Колмогоров, В.Л. Механика обработки металлов давлением: учеб. для вузов / В.Л. Колмогоров. - М.: Металлургия, 1986. - 688 с.

Судаков Николай Владимирович, д-р техн. наук, профессор, г. Челябинск; nv-sudakov@ yandex.ru.

Поступила в редакцию 15 мая 2017 г.

DOI: 10.14529/met170312

CONTROLLING OF SHEAR DEFORMATION

AND ITS ROLE IN FORMING A STRESS-STRAIN STATE

N.V. Sudakov, [email protected] Chelyabinsk, Russian Federation

Basing on a common mathematical model the variants of compressing rectangular blank with different schemes of contact friction forces are generalized. With a method of minimizing discrepancy for solving overdetermined system of equations stress and velocity field is determined. Potential for effective control of stress-strain state by targeted control of scheme of contact friction action is shown, that may be useful for decreasing deformation forces, forming material structure, compressing non-compact materials and for solving some special applications.

Torque of friction influence on character of stress-strain and deformed state in non-symmetrical compressing in blank volume is shown for the first time. Existence of hardened deformation zones and their displacement to the edges of a blank while increasing a non-symmetricity of scheme of external friction is determined. The fact that influence of torque on stress-strain state is growing with increase of friction constant and height to length ratio of a blank is shown.

Comparative evaluation of deformation heterogeneity in symmetrical and non-symmetrical compression modes is made. It depends on size ratio of a blank and is not one-valued.

Some recommendation is made on use of method of discrepancy minimization in solving overdetermined system of equations for stress-strain state in metal forming. Its effectiveness, simplicity and high accuracy in comparison to a variational method is shown.

Keywords: compression; deformation; stress-strain state; friction; torque; force; symmetrical; asymmetric; epure; minimization; discrepancy; overdetermined; system; equations.

References

1. Vydrin V.N, Ageev L.M. [Principal and Theoretical Foundations of a New Process of "Rolling-Drawing"]. Teoriya i tekhnologiya prokatki [Theory and Technology of Rolling]. No. 76. Chelyabinsk, ChPI Publ., 1971, pp. 3-21. (in Russ.)

2. Sudakov N.V. [Discrete-Local Method for Minimizing Residuals for Approximate Solution of Equations in Problems of Metal Forming]. Nauka YuUrGU. Materialy 66-y nauchnoy konferentsii [Science of SUSU. Materials of the 66th Scientific Conference]. Chelyabinsk, SUSU Publ., 2014, pp. 1669-1675. (in Russ.)

3. Kolmogorov V.L. Mekhanika obrabotki metallov davleniem [Mechanics of Metal Forming]. Moscow, Metallurgiya Publ., 1986. 688 p.

Received 15 May 2017

ОБРАЗЕЦ ЦИТИРОВАНИЯ

Судаков, Н.В. Управление деформацией сдвига и ее роль в формировании напряженного и деформированного состояний (НДС) / Н.В. Судаков // Вестник ЮУрГУ. Серия «Металлургия». - 2017. - Т. 17, № 3. -С. 99-105. DOI: 10.14529/теЙ70312

FOR CITATION

Sudakov N.V. Controlling of Shear Deformation and Its Role in Forming a Stress-Strain State. Bulletin of the South Ural State University. Ser. Metallurgy, 2017, vol. 17, no. 3, pp. 99-105. (in Russ.) DOI: 10.14529/metl70312