Scientific journal
Modern high technologies
ISSN 1812-7320
"Перечень" ВАК
ИФ РИНЦ = 0,940

MATHEMATICAL MODELING OF ELASTIC-PLASTIC STATE OF A CYLINDRICAL DOMAIN

Artemov M.A. 1 Baranovskiy E.S. 1 Potapov N.S. 1
1 Voronezh State University
The paper studies mathematical models of plastic-rigid and elastic-plastic bodies. Forms of notation for von Mises and Tresca yield criterions are compared in the case of plane strain state of plastic-rigid body. The case of existence of the solution of axisymmetric plane strain problem during the implementation of full plasticity regime for Tresca’s yield condition in the framework of theory of plastic flow is analyzed. The correct form of notation of plasticity regimes for general form of piecewise-linear plasticity conditions is considered. Elastic and plastic compressibility is taken into account in present mathematical model. An algorithm allowing to determine the conditions, under which plastic zone corresponding to the selected plasticity regime arises on the boundary of axisymmetric cylindrical region, when choosing a piecewise-linear plasticity conditions in the general form, is proposed. Expressions for determination of stress and strain for certain plasticity regimes, when one and two regimes of plasticity are realized in a plastic zone, are derived.
mathematical modeling
elastic-plastic body
plastic-rigid body
plane strain
piecewise-linear plasticity conditions
theory of plastic flow

Вопросам математического моделирования объектов, проявляющих упругопластические свойства, посвящено большое количество работ. Однако некоторые известные математические модели и алгоритмы решения задач ставят ряд вопросов. Так, осесимметричная задача для плоского деформированного состояния для сжимаемого упругопластического тела вызывает постоянный интерес и ряд дискуссий [12], обусловленных поиском правильного решения при рассмотрении кусочно-линейных условий пластичности. В работах [1, 2] обсуждались альтернативные формы записи кусочно-линейных условий пластичности. Общие соотношения при выполнении условия полной пластичности приведены в [8]. Решение конкретных задач при учете упругой и пластической сжимаемости для кусочно-линейных условий пластичности даны в [3, 5, 6]. В работе [4] обсуждалась статическая и кинематическая определимость для пластического тела. Ошибки, связанные с некорректной формой записи кусочно-линейных условий пластичности, обсуждались в [2, 7].

Модели жесткопластического тела

Для плоского деформированного состояния при выборе математической модели изотропного жесткопластического тела, когда пластический потенциал является четной функцией тензора напряжений, осевое напряжение σz является средним арифметическим двух других главных напряжений. Если функция пластичности и пластический потенциал совпадают, то для условия пластичности Треска осевое напряжение остается неопределенным [9]

artemov01.wmf

σmin ≤ σz ≤ σmax.

При этом задача плоской деформации будет локально статически определимой для компонент тензора напряжений σx, σy, σxy. Если выбирается условие пластичности Треска, то определение всех компонент тензора напряжений при плоской деформации возможно, когда вводится некоторое дополнительное предположение. Так, для осесимметричного состояния при выборе режима полной пластичности

artemov02.wmf (1)

в области a ≤ r ≤ b, на границе которой задано давление artemov03.wmf, имеет место

artemov04.wmf σθ = 2k + σr.

Здесь (r, θ, z) – цилиндрической системы координат.

Если, например, выбирается условие пластичности Мизеса, то в случае плоской деформации имеем [11]

artemov05.wmf artemov06.wmf

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

Ниже рассматривается плоское деформированное состояние, когда выбирается модель упругопластического тела.

Модель упругопластического тела

Для сжимаемого упругого тела в задаче Ламе осевое напряжение не всегда может быть средним [3]. Такая же особенность может иметь место и для упругопластического тела. Для упругого состояния цилиндрической области a ≤ r ≤ b, нагружаемой внешним давлением pb и внутренним давлением pa [11],

artemov08.wmf artemov09.wmf artemov10.wmf

artemov11.wmf (2)

Все функции пластичности вида artemov12.wmf (s – девиатор напряжений), для напряженного состояния (2) принимают наибольшее значение на границе r = a, поэтому пластическая зона будет зарождаться на этой границе.

Рассмотрим i-й режим кусочно-линейного условия пластичности общего вида (корректная форма записи для кусочно-линейного условия пластичности)

artemov13.wmf (3)

Условия, при выполнении которых на границе r = a при зарождении пластической зоны будет выполняться режим (3), определяются по следующему алгоритму: в соотношениях (3) компоненты напряжений заменяются с учетом формул (2) и r = a. Далее, используя неравенства f i–1 ≤ k и f i+1 ≤ k, находятся допустимые границы изменения одного из давлений: pa или pb. Реализуя указанный алгоритм, приходим к следующим формулам:

artemov14.wmf

Величины pb_max и pb_min являются наибольшим и наименьшим значениями pb, определяемым из равенств

artemov15.wmf

artemov16.wmf

которые получаются из равенств fi = fi–1, fi = fi+1 после подстановки в них (2).

Аналогично определяется диапазон допустимых значений давлений для pa (p a_min ≤ pa ≤ p a_max). Значения pa_max и pa_min определяются как наибольшее и наименьшее из следующих выражений:

artemov17.wmf

когда j = i – 1 и j = i + 1. В этом случае давление на внешней границе

artemov18.wmf

Для кусочно-линейных условий пластичности частного вида эти формулы были получены ранее в [10].

Пример реализации алгоритма определения режима пластичности

Наиболее известное решение задачи об определении напряженного и деформированного состояния для осесимметричной цилиндрической трубы a ≤ r ≤ b, на боковых границах которой заданы давления pa, pb или перемещения ua, ub, относится к режиму пластичности Треска

artemov19.wmf (4)

где k = const; κ = sign(σθ – σr). Из уравнения равновесия

artemov20.wmf

и системы (4), когда σr(r = a) = –pa, решая начальную задачу, находим окружное и радиальное напряжения:

artemov21.wmf

artemov22.wmf (5)

Из ассоциированного с условием (4) закона пластического течения, принимая гипотезу о естественном состоянии, следуют пропорции

artemov23.wmf

для пластических деформаций. Поскольку для осевой деформации

artemov24.wmf

то из закона Гука, учитывая формулы (5), находим, что

artemov25.wmf (6)

где ν – коэффициент Пуассона.

Из закона Гука, учитывая формулы (5), (6), определяем упругие деформации в пластической зоне

artemov26.wmf

artemov27.wmf

artemov28.wmf (7)

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

artemov29.wmf (8)

Учитывая формулы (7) и решая уравнение (8), находим перемещения и полные деформации в пластической зоне:

artemov30.wmf

artemov31.wmf

artemov32.wmf (9)

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

artemov33.wmf (10)

Из того, что на упругопластической границе r = C пластические деформации равны нулю, следует равенство C = 2κ(1 – ν2)c2.

Применяя предложенный выше алгоритм, находим, что режим (4) может зарождаться на границе r = a, если

artemov34.wmf artemov35.wmf

Выполнение одного режима пластичности

Решение (5)–(7), (9), (10) справедливо в области a ≤ r ≤ c пластического состояния, если радиус упругопластической границы

artemov36.wmf

Граница r1 находится из условия σr = σz. На границе r = r1 для напряжений имеет место равенство

artemov37.wmf

Из этой формулы видно, что безразмерное радиальное напряжение σr/k зависит только от значения коэффициента Пуассона и всегда будет положительным (растягивающие усилия).

Если c ≤ r1, то в пластической зоне реализуется только режим (4), когда

artemov38.wmf

Радиус упругопластической границы определяется из уравнения

artemov39.wmf

которое для выбранного алгоритма решения задачи, учитывая (5), находится из граничного условия artemov40.wmf. Для режима (4) из формул (5), (6) следует, что для радиуса упругопластической границы условие a ≤ c = r1 ≤ b выполняется, если

artemov41.wmf

В упругой зоне c ≤ r ≤ b коэффициенты A, B, входящие в формулы (2), можно найти, используя условие непрерывности напряжений на границе r = c. Так что в области c ≤ r ≤ b

artemov42.wmf

artemov43.wmf

σr = ν(σr + σθ).

Выполнение двух режимов пластичности

Из формул (5), (6) при r > r1 следует, что artemov44.wmf, поэтому в области пластического состояния a ≤ r ≤ c (r1 ≤ c) на границе r = r1 произойдет переход к режиму пластичности

artemov45.wmf (11)

Для этого режима задача статически неопределима. Поэтому определение напряженного состояния требует совместного рассмотрения всех уравнений (основных и определяющих).

Для режима (11) напряжения определяются по формулам [6]

artemov46.wmf

artemov47.wmf

artemov48.wmf (12)

Из условия непрерывности напряжений на границе r = r1 из формул (5) и (12) находим

artemov49.wmf

artemov50.wmf

artemov51.wmf

Если радиус упругопластической границы r1 < c, то в пластической зоне реализуются два режима: (4) и (11). В этом случае в упругой зоне c ≤ r ≤ b

artemov52.wmf

artemov53.wmf

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

artemov54.wmf

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

Выводы

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