CC BY
29
DOI 10.5862/JEST.254.16 УДК 539.3:539.214
Б.Е. Мельников, А.С. Семенов ИССЛЕДОВАНИЕ ПОЛЯ ПЛАСТИчЕСКИХ ПОДАТЛИВОСТЕй
при сложных пассивных путях нагружения
B.E. Melnikov, A.S. Semenov
STUDY OF PLASTIC COMPLIANCE FIELD UNDER COMPLEx PASSIVE LOADING
Представлены результаты экспериментальных исследований упругопластического деформирования при сложных переменных пассивных путях нагружения образцов из технически чистого никеля. Произведен анализ погрешностей эксперимента и предложены методы повышения точности испытаний. Исследованы свойства поля пластических податливостей в условиях пассивного нагружения. Для всех рассмотренных путей нагружения геометрические места равных пластических модулей податливости близки к окружностям. На основе полученных результатов предложены определяющие уравнения многоповерхностной теории пластичности с одной активной поверхностью равных пластических податливостей. Предложен метод определения констант модели, характеризующих эволюцию центра активной поверхности нагружения, по данным опытов на сложное непропорциональное нагружение, состоящее из первичного нагружения, разгрузки, вторичного нагружения в другом направлении и разгрузки. Представлены результаты верификации введенных определяющих уравнений на основе испытаний тонкостенных трубчатых образцов при сложных переменных (монотонное и циклическое) нагружениях, включающих частичные и промежуточные разгрузки.
ПЛАСТИЧЕСКОЕ ДЕФОРМИРОВАНИЕ МЕТАЛЛОВ; ПОЛЕ ПЛАСТИЧЕСКИХ ПОДАТЛИВОСТЕЙ; СЛОЖНОЕ ПАССИВНОЕ НАГРУЖЕНИЕ; ЭКСПЕРИМЕНТ; МОДЕЛИРОВАНИЕ.
The paper presents the findings of experimental studies on elastic-plastic deformation under complex-variable passive loading of commercially pure nickel samples. The analysis of experimental errors is performed and methods to improve test accuracy are proposed. The properties of the field of plastic compliances are investigated under passive loading. For all loading paths the shapes of equal plastic compliance are close to circles. On the basis of the obtained results, the constitutive equations of multisurface plasticity theory with one active surface of plastic compliances are suggested. A method for determining the model constants, characterizing the evolution of the center of the active surface is proposed based on the experiments under nonproportional complex loading, consisting of primary loading, unloading, secondary loading in a different direction and unloading. The verification of the proposed constitutive equations is performed on the basis of the experiments on thin-walled tubular specimens under complex (monotonic and cyclic) loadings, including partial and intermediate unloading.
PLASTIC DEFORMATION OF METALS; PLASTIC COMPLIANCE FIELD; COMPLEX PASSIVE LOADING; EXPERIMENT, MODELING.
Введение
Анализ усталостной прочности элементов конструкций, работающих в условиях сложного непропорционального циклического воздействия, требует использования в расчетах уточненных моделей неупругого деформирования. Получившая широкое распространение на практике теория пластического течения не учитывает накопление пластических деформаций при пассивном нагружении, предполагая, что разгрузка, нагружение внутри поверхности нагружения и по касательной к ней являются чисто упругими процессами. Однако в ряде случаев пассивного нагружения наблюдается накопление пластических деформаций. Поэтому необходимо разработать альтернативные модели пластичности, основой для которых могут служить результаты ряда экспериментальных и теоретических исследований [1—17] пассивного нагружения. В теории упруго-пластических процессов А.А. Ильюшина [3], эндохронной теории пластичности [4, 5] и обобщенной модели Прандтля (модель Мазинга) [6] учитывается пластическое деформирование при пассивном нагружении, однако определяющие уравнения и постоянные задаются исходя из режима активного нагружения без учета специфики процессов при пассивном нагружении.
Экспериментальное изучение свойств поля пластических податливостей [18, 19] послужило основой для разработки многоповерхностной теории пластичности с одной активной поверхностью [12—17]. Понятие пластических податливостей прямо или косвенно используется многоповерхностными теориями пластичности [20, 21]. Однако свойства пластических податливостей экспериментально исследованы слабо. Прежде всего это касается деформирования при пассивных путях нагружения. Учет накопления остаточных деформаций при таких нагружениях актуален не только при расчете деформирования [12—17, 22], но и при расчете накопления повреждений [11, 23, 24], анализе сверхпластичности [25], оценке надежности элементов машин и конструкций [26—28].
Целями работы были углубленное исследование процессов сложного пассивного нагружения и совершенствование определяющих уравнений многоповерхностной теории пластичности с од-
ной активной поверхностью на основе экспериментального изучения упругопластического деформирования трубчатых образцов из технически чистого никеля при непропорциональном сложном (монотонное и циклическое) нагружении, включающем частичные и промежуточные нагрузки.
Материал и образцы для испытаний
Материал (технически чистый никель, примеси - Si-0,068 %; Fe-0,025 %; Си-0,02 %) для проведения опытов выбран в связи с тем, что обладает сравнительно высоким модулем упругости — Е = 2,0740й Па, большой величиной отношения временного сопротивления к пределу текучести и малой ползучестью при нормальной температуре. Для данного материала установлено [17], что угол между плоскостью нагружения и плоскостью деформирования составляет около 5 %, а ошибка в вычислении приращения пластических деформаций как проекции на эту плоскость не превосходит 0,4 %. В работе [29] показано, что погрешность между удельными работами пластической деформации, определенными с учетом и без учета отклонения от подобия девиаторов, максимально достигает 1,73 %. Это весьма важно для последующего расчета поврежденности с использованием энергетических моделей накопления повреждений.
Образцы, на которых проводились испытания, вырезаны из холоднотянутых труб длиною около метра; наружные диаметры 8,06—8,10 мм, толщина стенки 0,18—0,20 мм. Это позволило создать практически однородное напряженно-деформированное состояние. По толщине стенки располагались не менее 6—7 зерен. Образцы браковались, если отклонения в замерах толщин стенок в десяти контролируемых точках составляли более чем ±0,005 мм. Несмотря на некоторый разброс в величинах толщин стенок вдоль окружности сечения, средние значения для рассматриваемого сечения отличались менее чем на ±0,001 мм. Если погонные веса вырезанных из одной трубы образцов отличались более чем на 1%, то образцы также браковались. По весу, диаметру и плотности (последняя определялась взвешиванием в дистиллированной воде и воздухе и равнялась 8,78 • 103 кг/м3) образца вычислялась
осредненная толщина стенки. При расхождениях (они были менее ±0,002 мм) с замерами в расчетах использовалась толщина, найденная по результатам взвешивания.
Образцы отжигали (одновременно для одной серии опытов) при 860 °С с последующим охлаждением в печи. На основании металлографического анализа установлено, что в продольном и поперечном шлифах зерна (средний поперечный размер 0,033 мм) — равноосные полиэдры со следами двойников. Испытания проводились через 20 дней после отжига. Каждый опыт проводился непрерывно.
Для крепления в установке на концы образца после отжига напаивались две втулки. При пайке во избежание перегрева рабочая часть образца омывалась водой. Перед отжигом, после него и после пайки проверялась прямолинейность образцов.
Установка для испытаний детально описана в [17].
Оценка погрешностей эксперимента
При анализе экспериментальных данных упруго-пластического деформирования материала рационально использовать векторное представление девиаторов напряжений и деформаций, введенное для случая плоской деформации Л. Прандтлем и позднее развитое В. Прагером, А.А. Ильюшиным [3] и др.
В совмещенном пространстве девиаторов напряжений = а, |з =л/3т) и пластических деформаций (эр = ер ,Эр =ур /л/3) для построения полей пластических податливостей h использо-
валась зависимость
Н =
8 Эр
где обозначено:
8Е| со8(а-Р)'
(1)
8 Эр =
^8ер)2 + 3(8ур)2; 18|^(8о)2 + 3(8т)2;
%/38т
а = аг^
8ур л/38е
в = аг^ -
(2)
8Ер ' ~ 8а в которых 8ст, 8х — изменения нормального и касательного напряжений за ступень нагруже-ния; 8е р и 8у р — изменения пластической осевой и сдвиговой деформаций.
Диаграмма строилась по средним значениям ст, х на ступени вторичного нагружения. Вычисление напряжений и деформаций проводилось в соответствии с методикой, предложенной в [17].
Расчеты выполнялись с учетом изменения размеров образцов в ходе опыта. Погрешности, получаемые с использованием соотношений (1), (2), соответствуют следующему разбросу данных:
в положении точек на поверхности напряжений разброс по ст доходил до ±2,5 • 106 Па, по - до ±3 • 106 Па;
в интенсивности пластических деформаций при вторичных нагружениях в точках, по которым строились геометрические места равных h, при h = 10-11 Па-1 погрешность в 8е р составляла ±0,001 %;
в направлении 8Э р при |8Э р| = 0,01 % погрешность доходила до ±5°; при |8Э р| = 0,05% — до ±1°.
В соответствии с методикой, изложенной в [17], учитывались погрешности, вызванные как погрешностями в измерениях, так и применением приближенных соотношений.
Методология проведения эксперимента
Методология включала одновременный отжиг 10 образцов; испытания начались через 20 дней и продолжались две недели. С каждым образцом испытания проводились непрерывно в течение 8—10 часов. Для контроля работы приборов и проверки отсутствия изгиба образца после каждой ступени нагружения снимались отсчеты. Определялись напряжения и деформации. Методология эксперимента подробно описана в [17].
Подготовительные опыты
Для изучения однородности образцов и начальной анизотропии материала были поставлены эксперименты при пропорциональном нагружении. Установлено, что при аппроксимации данных опытов уравнением
Ер = аак, (3)
где Ер, а1 — интенсивности пластических деформаций и напряжений, их отклонения не превосходят 2 % (к = 5,6; а = 1,25 • 1014 Па-5-6). Приняты модули упругости Е = 2,07 • 1011 Па и сдвига G = 0,75 • 1011 Па. При шести вариантах
пропорционального нагружения расхождение по величинам напряжений менее 3 %.
Диаграммы деформирования при нагруже-нии осевой силой и различными по знаку крутящими моментами практически совпали. Отклонения опытных точек от кривой, построенной по зависимости
л с р
Н = = ако,к-1, (4)
йо1 1
не превышают величину ±5 % по напряжениям. Геометрические места равных Дер (для значений 0,006; 0,012; 0,018; 0,036 %) и равных к (10-11; 2 • 10-11; 3 • 10-11; 5 • 10-10; 10-9 Па-1) могут отклоняться на такую же величину от соответствующих окружностей.
На основании данных по изменению скорости ползучести проведены эксперименты для первичного нагружения при нулевой выдержке, 15-минутной и 16-часовой выдержке. После этого образцы разгружались и нагружались по различным программам. Результаты испытаний для второго и третьего вариантов выдержки оказались близки. Сопротивление образца, не подвергавшегося выдержке, было приблизительно на 20 % ниже.
Величины выдержек выбирались по скорости ползучести, т. е. при интенсивности деформации е р0 < 1% применялась выдержка 15 мин., при ер = (3—4)% — 30 мин. и т.д.
Исследование поля пластических податливостей
При активном нагружении растягивающей силой и крутящим моментом в различных вариантах пассивного нагружения в ходе 16 опытов установлено, что значение к и направление 5Эр практически не зависят от 5£: наибольший поворот 5Эр , соответствующий изменению направления 5£ на 90°, составляет ± 2° от среднего направления.
Эволюция поверхности упругого состояния по допуску 0,05 % на образцах получивших пластическую деформацию 1,1; 1,9; 2,2 и 3,5 %, соответствует характерным свойствам, отмеченным в [30].
Для изучения влияния характера пути пассивного нагружения 10 образцов были подвергнуты растяжению, 3 образца — кручению и 4 образца первично нагружались по пути о = \/3т .
Далее задавались различные пути пассивного нагружения. Результаты опытов на 10 образцах при различных путях пассивного нагружения показали, что величины вторичных пластических деформаций, направления векторов их приращений, величины к в первом приближении определяются первичным нагружением и текущим напряженным состоянием.
По достижении заданного напряженного состояния осуществлялась 15-минутная выдержка, после этого образец частично разгружался. При некотором, предварительно намеченном для данного образца, напряженном состоянии разгрузка прекращалась и начиналось вторичное нагружение. На плоскости напряжений вторичный путь составлял прямой угол с путем первичного нагружения.
На каждой ступени нагружения определялись величины, необходимые для расчета характеристик анизотропии. Вторичное нагружение доводилось до величины интенсивности напряжений, близкой к о.0, достигнутой при первичном на-гружении. Затем обследованный вторичный путь проходился в обратном направлении до выхода на путь первичного нагружения и производилось нагружение образца по первичному пути до получения следующего значения о.0. Опять давалась 15-минутная выдержка под полной нагрузкой. В такой же последовательности производилась разгрузка и вторичное нагружение образца. Достижение на том же образце нового напряженного состояния и последующее вторичное нагружение производились так же, только время выдержки образца при полной нагрузке было увеличено до 20 минут.
По вторичным путям, составляющим прямой угол с направлением пути первичного нагружения, было испытано 8 образцов и 2 образца нагружались по «косым» вторичным путям. Результаты всех этих испытаний позволили построить искомые участки геометрических мест равных к и равных Дер.
В описанных выше испытаниях как при первичном, так и при вторичном нагружении одновременно изменяли осевую силу и крутящий момент.
В основных опытах по изучению деформационной анизотропии реализованы 16 вариантов
нагружения, испытано 54 образца, пройдено 118 вторичных путей нагружения.
Опыты поставлены для трех вариантов ин-тенсивностей напряжений первичного нагружения — 103, 150, 306 МПа, достигавшихся различными путями.
Построены геометрические места равных пластических податливостей h для первичных нагружений соответственно осевой силой, крутящим моментом, по пути а = >/3т и по двух-звенному пути.
Получены геометрические места равных h для первичных нагружений осевой силой (рис. 1), крутящим моментом и по двухзвенному пути. Пути вторичного нагружения показаны тонкими линиями. Во всех случаях пропорционального нагружения центры окружностей лежат на линии первичного нагружения.
На рис. 2 изображены опытные точки и окружности, соответствующие равным значениям h, при различных величинах первичного нагружения, что демонстрирует влияние достигнутого при первичном нагружении о.0 на размеры и положение окружностей равных h. Центры окружностей располагаются на отрезке о-о".
Во всех изученных случаях деформационная анизотропия проявлялась при первых пластиче-
Т3т , МПа
л/3т, МПа
о, МПа
Рис. 1. Геометрические места равных пластических податливостей h для первичных нагружений осевой силой. Для всех кривых одинаковое значение о.0 = 103 МПа и различные h
200
о, МПа
Рис. 2. Геометрические места равных пластических податливостей h для первичных нагружений осевой силой. Для всех кривых одинаковое значение h = 10-10 Па-1 и различные о.0
ских деформациях (около 0,2 %) и не исчезала при больших значениях е Геометрические места равных h аппроксимировались окружностями, положения центров и радиусы которых находились по методу наименьших квадратов. При этом отклонение не превосходило 2—3 % величины радиуса, т.е. имело тот же порядок, что и разброс, обусловленный неоднородностью образцов.
Во всех случаях пропорционального нагруже-ния эти центры лежат на линии первичного на-гружения. При сложных путях первичного нагружения они смещались к большей проекции Э0Р.
Форма и свойства поверхностей равных остаточных деформаций соответствуют приведенным в работах [30, 31] и многочисленным экспериментальным данным для поликристаллических металлов.
При вторичных нагружениях, не выходящих за пределы круга Мизеса, соответствующего последней точке первичного нагружения, поле равных h в первом приближении однозначно определяется историей первичного нагружения.
Все геометрические места равных h в диапазоне 10-11 Па-1 < h < 10-10 Па-1 близки к окружностям. Центры этих окружностей смещены относительно начала координат. Направление
вектора смещения центра соответствует отношению компонент первичной пластической деформации, модуль этого вектора зависит только от значения к, отвечающего данной окружности, и не зависит от истории нагружения. Радиусы окружностей не зависят от характера пути первичного нагружения и определяются соответствующим значением к и интенсивностью напряжений о достигнутой при первичном нагружении. Чтобы определить размеры и положение этих окружностей, кроме величин к и о , необходимо располагать значениями пяти констант материала, определяемых по результатам трех опытов.
При к < 1-10-11 Па-1 геометрические места равных к приближенно можно считать окружностями, проходящими через последнюю точку пути первичного нагружения. Центры этих окружностей располагаются на той же прямой, что и центры окружностей, отвечающих большим значениям к.
Величины вторичных пластических деформаций модулей податливости и направления векторов приращений вторичных пластических деформаций в первом приближении однозначно определяются лишь историей первичного на-гружения и текущим напряженным состоянием при вторичном нагружении. Следовательно, деформационную анизотропию, созданную данным первичным нагружением, можно характеризовать одним полем к и одним полем направлений 5Эр.
Из-за понижения точности измерений при к < 10-11 Па-1 зависимости, характеризующие поля, могут быть менее точными.
Расчет процессов деформирования проводился с помощью соотношений [17]
Ро =
А £,
0 .
к с
к = ак
ю
Я + В
о I о + В
А Ло_ к
Б- В
С о к-1
Базовые эксперименты и методика определения постоянных
Значения постоянных в соотношениях (5), (6) — а = 3,15-10-46 Па-56; к = 5,6; А = 3,610-4; В = 3,2-109 Па; D = 4,12-109 Па — получены осреднением данных ряда опытов. Величины а и к могут быть определены из опыта на линейное растяжение (в результате аппроксимации диаграммы испытаний прямой линией, построенной в координатах ^о—^ер). Погрешность в вычислениях этих постоянных определяется неоднородностью материала, которая дает отклонения не более ±2 %. Среднеквадратические погрешности а и к равны соответственно 4 % и 0,3 % при о = 107 Па.
Для определения константы А необходимо знать величину хотя бы одного смещения центра окружности равной пластической податливости. Такая окружность может быть построена по данным двух опытов на вторичное нагружение, первичные нагружения которых были одинаковыми. Наибольшая погрешность определения А будет при использовании окружностей, отвечающих меньшим значениям к, т.к. р в этом случае велико.
Для определения постоянных В и Б необходимы значения минимум двух радиусов —Я1 и Я^, соответствующих разным к. Модуль вектора смещения центра поверхности |р| определяется величиной модуля пластической податливости к. Сопоставление величин |р| позволяет установить зависимость
Р =
к
(7)
(5)
(6)
где ро , р^ — проекции вектора смещения поверхности равных пластических податливостей; а, к, А, В, D — постоянные материала.
Направление вектора р определяется соотношением компонент первичных пластических деформаций (5).
Для установления зависимости радиуса окружности равной пластической податливости от о и для определения постоянных в координатах о — Я были нанесены точки, соответствующие изучаемым значениям о и к (рис. 3). Точки, отвечающие одинаковым значениям к, располагаются весьма близко от прямых, пересекающихся в точке С.
Я, МПа
300-
200 -
100 ■
-100 -
М0^ Па1
2110"^ Па"'
-о—й = З'Ю"5 Па"1
—*—к = 5'10^ Па1
10' 10"" Па"'
к-1
>10
получаем
Ъ о + В . Б
к-.1— + В ! ак
К = °о±В (к-1А + в) - В; Б У ак
(
к = ак
Я + В
\°1 о + В
\к-1
Б - В
С помощью соотношений (7) и (11) найдены значения |р| и Я для окружностей, отвечающих значениям h, 10-11Па-1, равным 1, 2, 3, 5 и 10. Из (10) следует
Ъ о + в
Я1 + в
Ъ о + В
^ + В -. (13)
° к+ в ° к -¡¡^ + в
\ак у ак
Отсюда могут быть найдены константы В и Б:
300
о.„, МПа
¡0'
Рис. 3. Зависимость радиуса окружности равных пластических податливостей от интенсивности напряжений первичного нагружения
Ее координаты соответствуют величине В, которую можно рассматривать как некоторую постоянную материала, определяющую развитие деформационной анизотропии. Для сопоставления поведения материала при вторичных нагруже-ниях с его поведением при первичном нагружении (рис. 4) были нанесены точки, соответствующие величинам радиусов Я0 окружностей равных h для исходного материала, которые получены из (4):
1
В =
А к-1 - А к-\\И2 2 V ак 1 V ак
(А - А)+
Б = (Ъ 0 + в)
А - А
(14)
(8)
Экспериментально определенные точки располагались на вертикальной прямой с абсциссой о* = 9,2-105 Па. Тогда для любого луча
+ В Я + В
(9)
о + В Яо + В Полагая о* + В = Б, после подстановки в (8)
Я + В
(10)
Откуда могут быть получены выражения для Я и к
(11)
(12)
При этом чем больше разность (Я:—Я2) радиусов, тем меньше погрешность В и Б. Поэтому при постановке опытов следует выбирать возможно большее о , а окружности равной пластической податливости к строить для наименьшего и наибольшего достоверно определяемых значений к.
По описанной методике из данных двух опытов (о.п = 3,42 • 108 Па, к = 10-11 Па-1 и к = 10-10 4 ¡0 ' '1 2
Па-1) определены В = 3,32-108 Па, Б=4,24 • 108 Па, а по осредненным данным всех опытов — В = 3,2 • 108 Па, Б=4,12 • 108 Па.
Отличия в результатах расчетов по соотношениям (4), (5) с использованием этих значений оказались невелики.
Таким образом, для определения констант А, В и Б достаточно иметь две окружности равных к. Чтобы построить такую окружность, центр которой лежит на векторе первичной пластической деформации, необходимо знать две точки окружности, достаточно удаленные друг от друга. Окружности, отвечающие о0 = 3,42 • 108 Па, были построены по двум точкам. Средние точки (см. рис. 2) оказались близки к этим окружностям. Также не дали значительных отклонений от этих окружностей точки, полученные в ряде проверочных опытов (см. рис. 4).
Тэт, МПа
/г=10"10Па"1
ТЭтМПа
о, МПа
Рис. 4. Геометрические места равных пластических
податливостей к для различных первичных нагружений (1 — осевое растяжение; 2 — крутящий момент; 3 — нагружение о = л/3т )
Исследован вопрос о направлении векторов 5Э по отношению к нормали к поверхностям равных бер. Проведены испытания первично-растянутых образцов (см. рис. 4), первичноза-крученных крутящим моментом, а также нагруженных первично либо пропорциональными, либо сложными путями до о = (соответственно точки 1, 2, 3).
Построенные по приближенным формулам (7), (11) окружности соответствуют следующим значениям о.0 , 108 Па: 1; 1,27; 1,47; 1,77; 2,19; 2,62; 3 и 3,35.
Координаты экспериментальных точек преобразовывались так, чтобы направление вектора первичной пластической деформации Э0Р повернулось до совпадения с осью о (при нагруже-нии крутящим моментом — на 90°; по пути о = л/3т — на 45° и т. д.).
На рис. 5 показаны геометрические места равных бе? (кривая 1— 0,006 %, 2 — 0,012 %, 3 — 0,018 %, 4 — 0,036 %) при первичном нагру-жении растягивающей силой до о0 = 1,48 • 108 Па. Направления вторичных нагружений изображены штрих-пунктирными линиями. Для ряда путей вторичных нагружений стрелками показаны направления векторов 5ЭР.
Отклонения 5Э от нормалей увеличиваются при удалении от точки окончания первичного
о МПа
Рис. 5. Геометрические места равных 8ер . (1 — 0,006 %; 2 — 0,012 %; 3 — 0,018 %; 4 — 0,036 %) при первичном нагружении растягивающей силой (о0 = 148 МПа; е0 = 1,4 %)
нагружения. В среднем эти отклонения составляли 5°. Практическое использование геометрических мест равных бер для определения направления векторов 5ЭР затруднено, так как формы этих поверхностей существенно меняются в зависимости от величины бер [30].
Отклонения от нормалей к поверхностям равных к составляло 6—7°, т. е. немного больше отклонений от поверхностей равных бер .
Определяющие уравнения многоповерхностной теории при одной активной поверхности
Рассмотрим вариант определяющих уравнений, учитывающих описанные выше особенности упруго-пластического деформирования при сложном пассивном нагружении:
й о
= 4 Б--(йе-йе р); (15)
.Р_,д/ (о л)
й е р = к-
д О
д/ (О А )
до
-й о,
(16)
где е и ер — тензоры полной и пластической деформации; /(оА) = /к+1(оА) — соответствует активной (к+1)-й поверхности; оА = о — о0 — тензор активных напряжений; о0 — центр (к + 1)-й поверхности равных податливостей, характеризующий ее смещение как жесткого целого; к — модуль пластической податливости (к + 1)-й
поверхности; 4Б — тензор модулей упругости четвертого ранга.
Использование окружностей в качестве двумерных сечений поверхности равной податливости позволяет ввести в рассмотрение гиперсферу Мизеса
ЯоА ) =/(нА) = ^ (3/28 А 8 А ), (17)
где sA — девиатор тензора активных напряжений оА. В этом случае закон течения (16) принимает вид
1 р 3 8А ® 8А ■■ds й е р =-к—-А-.
(18)
А А
Уравнения для величин, характеризующих положение центров поверхностей So., необходимые для вычисления йеР по (18), вводятся на основе экспериментально наблюдаемой закономерности, устанавливающей, что центр поверхности б р. + 1 находится на прямой, соединяющей б р. и точку реверса б [14, 15]:
Ъ0] + 1 = *0] + й] ^К, ^0]
(19)
^я^ 8я, ^а)"(8д+1 вя,)
ач =-
Сч' Уяя+1 Яя/.
(в я ^Ъдч ) "(в я+1 в я )
, ]=1,..., к, (20)
где Бд. — девиатор напряжений, соответствующий реверсу нагрузки, под которым понимаем точку излома пути нагружения с последующим приращением напряжений, направленным внутрь поверхности равной податливости. Условие возникновения нового, у-го реверса можно определить следующим образом:
до
■йо < 0.
(21)
Возникновение нового, у-го реверса ст0/ обусловливает появление новой (/+1)-й поверхности. Для вычисления местоположения центра активной поверхности в (20) вместо + 1 следует использовать значение девиатора напряжений б.
Уравнения эволюции внутренних переменных Бр. (19), (20) имеют голономный рекуррентный характер, что позволяет учесть предысторию
процесса нагружения без введения дифференциальных уравнений.
Модуль пластической податливости к, фигурирующий в уравнениях (16) и (18), является функцией активных напряжений и напряжений в точке реверса наибольшей поверхности за всю историю нагружения Зд.. Сопоставление результатов расчетов и опытов показало, что хорошую точность обеспечивает степенная аппроксимация (12).
Предложенный вариант теории является обобщением на трехмерный случай рассмотренных выше особенностей поведения материала при сложном нагружении, а также удовлетворяет термодинамическим ограничениям [32]. Особенности интегрирования введенных определяющих уравнений (15), (18)—(20) рассмотрены в [15].
Сравнение расчетных и экспериментальных путей деформирования
Пути первичного и вторичного нагружения в пространстве напряжений показаны на рис. 6 сплошной линией. На рис. 7 представлены пути деформирования, отвечающие путям вторичного нагружения. Деформации, вычисленные с помощью многоповерхностной теории, показаны на рис. 7 пунктирной линией, с использованием теории пластического течения — штрих-пунктирной линией, экспериментальные данные маркированы сплошной линией. Отличие от экспериментальных данных деформаций, рас-читанных по многоповерхностной теории, не превышает 0,009 %, а по теории пластического течения — 0,09 % для ех, 0,125 % для еху.
На рис. 6 также показаны пути нагружения, полученные по теории течения (штрих-пунктирная линия) и по исследуемой теории (пунктирная линия), отвечающие наблюдаемому в опыте (сплошная линия) пути вторичного деформирования. Наибольшее отклонение от действительного пути нагружения составляет 1,08 • 107 Па для рассматриваемой теории и 3,33 • 107 Па для теории течения.
n/3t , МПа
* * \ 1 ...... 1 -i
k К - 11 А 1 11 } * ( - - -4
и 1-1-т--1-1-1-1-1-
-50 0 50 100 о, МПа
Рис. 6. Пути нагружения: 1 — эксперимент (сплошная линия); 2 — многоповерхностная теория (пунктирная линия); 3 — теория течения (штрих-пунктирная линия)
0,3-
0,2-
0,1 -
0,0
-0,2 -0,1 0,0 е, %
Рис. 7. Пути деформирования: 1 — эксперимент (сплошная линия); 2 — многоповерхностная теория (пунктирная линия); 3 — теория течения (штрих-пунктирная линия)
Заключение
Для всех рассмотренных путей нагружения геометрические места равных пластических модулей податливости к близки к окружностям, центры которых лежат на прямой, соединяющей исходную точку первоначального нагружения и точку реверса. Приращение пластической деформации направлено по нормали к этим окружностям.
Радиусы окружностей, соответствующих геометрическим местам равных пластических по-датливостей не зависят от характера пути первичного нагружения и определяются значением к и интенсивностью напряжений о , до которой доводилось первичное нагружение.
Выполнена систематическая оценка влияния различных источников погрешности (неточности в геометрии образцов, погрешности измерительной аппаратуры, статистического разброса характеристик и микроструктуры материала образцов) на точность результатов при проведении экспериментов в условиях сложного пассивного нагружения.
Для определения констант модели, характеризующих эволюцию центра активной поверхности нагружения, предложен метод, использующий результаты опытов на сложное непропорциональное нагружение, состоящее из двух участков на-гружений.
Сравнение прогнозов многоповерхностной теории при одной активной поверхности с данными опытов на сложное пассивное нагружение продемонстрировало их хорошее соответствие (отличие менее 5 % для рассматриваемых траекторий деформирования).
Работа выполнена при финансовой поддержке гранта Российского Фонда Фундаментальных исследований (РФФИ) № 16-08-00845-а.
Y/лА%
■ ч. "Н.
ji ЛЯ \
JJ JB L if *
т
СПИСОК ЛИТЕРАТУРЫ
1. Зубчанинов В.Г., Гультяев В.И., Зубчанинов Д.В.
Экспериментальное исследование процессов сложного нагружения материалов на многозвенных траекториях // Современные проблемы термовязко-пластичности. Труды II школы-семинара. М., 2007. С. 19-24.
2. Зубчанинов В.Г. Об активных и пассивных процессах сложного нагружения-разгружения // Совре-
менные проблемы термовязкопластичности. Труды II школы-семинара. М., 2007. С. 3—18.
3. Ильюшин А.А. Пластичность. Основы общей математической теории. М.: Изд-во АН СССР, 1963. 271 с.
4. Кадашевич Ю.И., Мосолов А.Б. Эндохронные теории пластичности: основные положения, перспективы развития // Известия Российской академии наук. Механика твердого тела. 1989. № 1. С. 161.
5. Кадашевич Ю.И., Мосолов А.Б. О соотношениях эндохронной теории пластичности с «новой» мерой внутреннего времени // Прикладная математика и механика. 1990. Т. 54, № 4. С. 689.
6. Пальмов В.А. Колебания упруго-пластических тел. М.: Наука, 1976. 328 с.
7. Kuwabara Т., Ishiki M., Kuroda M., Takahashi S. Yield locus and work hardening behavior of a thin-walled steel tube subjected to combined tension-internal pressure // J. Phys. IV France 105. 2003. P. 347-354.
8. Smith S., Almroth B.O. An experimental investigation of plastic flow under biaxial stress // Experimental mechanics. 1970. № 10. P. 217-224.
9. Khan A.S., Parikh Y. Large deformation in poly-crystalline copper under combined tension-torsion, loading, unloading and reloading or reverse loading: a study of two incremental theories of plasticity // International Journal of Plasticty. 1986. Vol. 2. P. 379-392.
10. Rees D.W.A. Research note on a reexamination of neutral loading experiments // International Journal of Plasticity. 1988. Vol. 4. P. 91-102.
11. Павлов П.А. Основы инженерных расчетов элементов машин на усталость и длительную прочность. Л: Машиностроение, 1988. 251 с.
12. Мельников Б.Е., Изотов И.Н., Кузнецов Н.П. Пространственный упругопластический элемент для описания нестационарного деформирования // Известия ВНИИГ им. Б.Е. Веденеева. 1986. Т. 197. С. 75-79.
13. Мельников Б.Е., Изотов И.Н., Кузнецов Н.П.
Расчет и экспериментальное исследование сложных путей упругопластического деформирования // Проблемы прочности. 1990. № 8. С. 14-17.
14. Мельников Б.Е., Семенов А.С. Многоповерхностная теория пластичности с одной активной поверхностью пластических податливостей // Труды СПбГТУ. 1991. № 144. С. 26-31.
15. Мельников Б.Е., Семенов А.С. Описание эволюции поверхностей равной податливости в многоповерхностной теории пластичности // Труды СПбГТУ. 1996. № 456. С. 52-61.
16. Изотов И.Н., Кузнецов Н.П., Мельников Б.Е., Митюков А.Г. Упругопластическое поведение стальных трубчатых образцов при сложном трехосном на-гружении // Научно-технические ведомости СПбГТУ, 2003. № 3(33). С. 114-125.
17. Изотов И.Н., Мельников Б.Е., Семенов А.С., Семенов С.Г. Определение свойств пластических по-датливостей на никелевых образцах // Научно-технические ведомости СПбГПУ, 2009. № 2(78). С. 165-174.
18. Изотов И.Н., Ягн Ю.И. Изучение пластического деформирования металла с деформационной анизотропией, созданной в процессе предварительного нагружения // ДАН СССР. 1961. Т. 139, № 2(78). С. 576-579.
19. Изотов И.Н. Изучение анизотропии металла созданной некоторыми видами пластического де-
формирования: Автореф. ... дисс. на соиск. уч. ст. канд. тех. наук. Ленинград. 1962.
20. Mroz Z. On the destription of anisotropic workhardening // J. Mech. Phys. Solids. 1967. Vol. 15. № 3. P. 163-175.
21. Mroz Z., V.A., Zienkiewicz O.C. An anisotropic hardening model for soils and its application to cyclic leading // Int. J. for Numerical and analytical Methods in Geomechanics. 1978. Vol. 2. № 3. P. 203-221.
22. Melnikov B.E., Semenov A.S. Creation and application of hierarchical sequence of material models for numerical analysis of elasto-plastic structures // Zeitschrift fur Angewandte Mathematik und Mechanik. 1996. Vol. 76. № Suppl. 2. P. 615-616.
23. Мельников Б.Е., Кузнецов Н.П., Митюков А.Г., Семенов А.С., Семенов С.Г. Применение модели энергетической оценки прочности при нескольких циклах нагружения в случае сложного неоднород-но-го напряженного состояния // Современные проблемы ресурса материалов и конструкций. Труды III школы-семинара М.: Изд-во МГТУ (МАМИ), 2009. С. 140-147.
24. Мельников Б.Е., Семенов А.С. Энергетическая модель накопления повреждений при малоцикловой усталости // Сб. матер.: XXI Петербургские чтения по проблемам прочности. К 100-летию со дня рождения Л.М. Качанова и Ю.Н. Работнова. 2014. С. 59-61.
25. Kitaeva D.A., Rudaev Y.I. About kinetic equations of dynamic superplasticity model // Proceedings of the XXXI Summer School «Advanced Problems in Mechanics». 2003. St. Petersburg: Polytechnical Publishing House, 2003. P. 172-176.
26. Semenov A., Semenov S., Nazarenko A., Getsov L. Computer simulation of fatigue, creep and thermal-fatigue cracks propagation in gas-turbine blades // Materi-ali in Tehnologije. 2012. Vol. 46. № 3. P. 197-203.
27. Артюх В.Г. Точность предохранителей для металлургических машин. Мариуполь: Изд-во ПГТУ, 2000. 177 с.
28. Melnikov B.E., Semenov A.S. Multimodel analysis as strategy of reliable description of elastic-plastic deformation processes // Proc. of IV International Conference on Computational Plasticity. Fundamentals and Applications. 1995. P. 181-189.
29. Касап К. Учет отклонения от подобия девиа-торов при определении работы пластической деформации // Periodica Polytechnica Mechanical engineering. 1971. Т. 13. № 1. С. 17-25.
30. Васин Р.А. Определяющие соотношения теории пластичности // Итоги науки и техники. Серия: Механика деформируемого тела. 1990. Т. 21. С. 3-75.
31. Ягн Ю.И., Шишмарев О.А. Некоторые результаты исследования границ упругого состояния пластически растянутых образцов никеля // Доклады АН СССР. 1958. Т. 119, № 1. С. 46-48.
32. Семенов А.С. Анализ термодинамической состоятельности многоповерхностной теории пластичности с одной активной поверхностью // На-
учно-технические ведомости Санкт-Петербургского государственного политехнического университета. 2003. № 3(33). С. 103-113.
REFERENCES
1. Zubchaninov v.G., Gultyayev v.l., Zubchaninov
D.v. Eksperimentalnoye issledovaniye protsessov slozh-nogo nagruzheniya materialov na mnogozvennykh trayektoriyakh. Sovremennyye problemy termovyazkoplas-tichnosti. Trudy IIshkoly-seminara. M., 2007. S. 19—24.
2. Zubchaninov v.G. Ob aktivnykh i passivnykh protsessakh slozhnogo nagruzheniya-razgruzheniya. Sovremennyye problemy termovyazkoplastichnosti. Trudy II shkoly-seminara. M., 2007. S. 3—18.
3. Ilyushin A.A. Plastichnost. Osnovy obshchey matematicheskoy teorii. M.: Izd-vo AN SSSR, 1963. 271 s. (rus.)
4. Kadashevich Yu.i., Mosolov A.b. Endokhronnyye teorii plastichnosti: osnovnyye polozheniya, perspektivy razvitiya. Izvestiya Rossiyskoy akademii nauk. Mekhanika tverdogo tela. 1989. № 1. S. 161. (rus.)
5. Kadashevich Yu.I., Mosolov A.b. O sootnosheni-yakh endokhronnoy teorii plastichnosti s «novoy» meroy vnutrennego vremeni. Prikladnaya matematika i mekhanika. 1990. Vol. 54. № 4. S. 689. (rus.)
6. Palmov v.A. Kolebaniya uprugo-plasticheskikh tel. M.: Nauka. 1976. 328 s. (rus.)
7. Kuwabara T., ishiki M., Kuroda M., Takahashi s. Yield locus and work hardening behavior of a thin-walled steel tube subjected to combined tension-internal pressure. J. Phys. IV France 105. 2003. P. 347-354.
8. smith s., Almroth B.o. An experimental investigation of plastic flow under biaxial stress. Experimental Mechanics. 1970. № 10. P. 217-224.
9. Khan A.s., Parikh Y. Large deformation in poly-crystalline copper under combined tension-torsion, loading, unloading and reloading or reverse loading: a study of two incremental theories of plasticity. International Journal of Plasticty. 1986. Vol. 2. P. 379-392.
10. Rees D.W.A. Research note on a reexamination of neutral loading experiments. International Journal of Plasticity. 1988. Vol. 4. P. 91-102.
11. Pavlov P.A. Osnovy inzhenernykh raschetov ele-mentov mashin na ustalost i dlitelnuyu prochnost. L: Mashinostroyeniye, 1988, 251 s. (rus.)
12. Melnikov B.E., izotov i.N., Kuznetsov N.P. Pros-transtvennyy uprugoplasticheskiy element dlya opisaniya nestatsionarnogo deformirovaniya. Izvestiya VNIIG im. B.Ye. Vedeneyeva, 1986. T.197. S. 75-79. (rus.)
13. Melnikov B.E., izotov i.N., Kuznetsov N.P. Ra-schet i eksperimentalnoye issledovaniye slozhnykh putey uprugoplasticheskogo deformirovaniya. Problemy proch-nosti. 1990. №8. S. 14-17. (rus.)
14. Melnikov B.E., semenov A.s. Mnogopoverkh-nostnaya teoriya plastichnosti s odnoy aktivnoy poverkh-nostyu plasticheskikh podatlivostey. Trudy SPbGTU. 1991. №144. S. 26-31. (rus.)
15. Melnikov B.E., semenov A.s. Opisaniye evolyutsii poverkhnostey ravnoy podatlivosti v mnogopoverkhnost-noy teorii plastichnosti. Trudy SPbGTU. 1996. № 456. S. 52-61. (rus.)
16. izotov i.N., Kuznetsov N.P., Melnikov B.E., Mi-tyukov A.G. Uprugoplasticheskoye povedeniye stalnykh trubchatykh obraztsov pri slozhnom trekhosnom nagru-zhenii. Nauchno-tekhnicheskiye vedomosti SPbGTU. 2003. № 33. S. 114-125. (rus.)
17. izotov i.N., Melnikov B.E., semenov A.s., semenov s.G. Opredeleniye svoystv plasticheskikh podatlivostey na nikelevykh obraztsakh. Nauchno-tekhnicheskiye vedomosti SPbGPU. 2009. № 2(78). S. 165-174. (rus.)
18. izotov i.N., Yagn Yu.I. Izucheniye plastichesk-ogo deformirovaniya metalla s deformatsionnoy aniz-otropiyey, sozdannoy v protsesse predvaritelnogo nagruzheniya. Doklady AN SSSR. 1961. T. 139. № 3. S. 576-579. (rus.)
19. izotov i.N. Izucheniye anizotropii metalla soz-dannoy nekotorymi vidami plasticheskogo de-formirovaniya: Avtoref. diss. ... na soisk. uch. st. kand. tekh. nauk. Leningrad. 1962. (rus.)
20. Mroz Z. On the destription of anisotropic work-hardening. J. Mech. Phys. Solids. 1967. Vol. 15, № 3. P. 163-175.
21. Mroz Z., Norries v.A., Zienkiewicz o.c. An aniso-tropic hardening model for soils and its application to cyclic loading. Int. J. for Numerical and analytical Methods in Geomechanics. 1978. Vol. 2, № 3. P. 203-221.
22. Melnikov B.E., semenov A.s. Creation and application of hierarchical sequence of material models for numerical analysis of elasto-plastic structures. Zeitschrift fur Angewandte Mathematik und Mechanik. 1996. Vol. 76. № Suppl. 2. P. 615-616.
23. Melnikov B.E., Kuznetsov N.P., Mityukov A.G., semenov A.s., semenov s.G. Primeneniye modeli ener-geticheskoy otsenki prochnosti pri neskolkikh tsiklakh nagruzheniya v sluchaye slozhnogo neodnorodnogo napryazhennogo sostoyaniya. Sovremennyye problemy resursa materialov i konstruktsiy. Trudy III shkoly-semina-ra. M.: Izd-vo MGTU (MAMI), 2009. S. 140-147. (rus.)
24. Melnikov B.E., semenov A.s. Energeticheskaya model nakopleniya povrezhdeniy pri malotsiklovoy ust-alosti. Sb. mater.: XXI Peterburgskiye chteniya po prob-lemam prochnosti. K 100-letiyu so dnya rozhdeniya L.M. Kachanova i Yu.N. Rabotnova. 2014. S. 59-61. (rus.)
25. Kitaeva D.A., rudaev Y.I. About kinetic equations of dynamic superplasticity model. Proceedings of the XXXI Summer School«Advanced Problems in Mechanics». 2003. St. Petersburg: Polytechnical Publishing House, 2003. P. 172-176.
26. Semenov A., Semenov S., Nazarenko A., Getsov L.
Computer simulation of fatigue, creep and thermal-fatigue cracks propagation in gas-turbine blades. Materiali in Tehnologije. 2012. Vol. 46, № 3. P. 197-203.
27. Artyukh V.G. Tochnost predokhraniteley dlya metallurgicheskikh mashin: Mariupol: Izd-vo PGTU, 2000. 177 s. (rus.)
28. Melnikov B.E., Semenov A.S. Multimodel analysis as strategy of reliable description of elastic-plastic deformation processes. Proc. of IVInternational Conference on Computational Plasticity. Fundamentals and Applications. 1995. P. 181-189.
29. Kasap K. Uchet otkloneniya ot podobiya deviato-rov pri opredelenii raboty plasticheskoy deformatsii. Peri-
odica Polytechnica Mechanical engineering. 1971. Vol. 13, № 1. P. 17-25.
30. Vasin R.A. Opredelyayushchiye sootnosheniya teorii plastichnosti. Itogi nauki i tekhniki. Seriya mekhan-ika deformiruyemogo tela. 1990. Vol. 21. S. 3-75. (rus.)
31. Yagn Yu.I., Shishmarev O.A. Nekotoryye rezultaty issledovaniya granits uprugogo sostoyaniya plasticheski rastyanutykh obraztsov nikelya. Doklady AN SSSR. 1958. Vol. 119, № 1. S. 46-48. (rus.)
32. Semenov A.S. Analiz termodinamicheskoy sostoyatelnosti mnogopoverkhnostnoy teorii plastichnosti s odnoy aktivnoy poverkhnostyu. Nauchno-tekhnicheskiye vedomosti Sankt-Peterburgskogo gosudarstvennogo politekhnicheskogo universiteta. 2003. № 3(33). S. 103113. (rus.)
СВЕДЕНИЯ ОБ ABTOPAX/AUTHORS
МЕЛЬНИКОВ Борис Евгеньевич — доктор технических наук заведующий кафедрой Санкт-Петербургского политехнического университета Петра Великого. 195251, Россия, г. Санкт-Петербург, Политехническая ул., 29. E-mail: [email protected]
MELNIKOV Boris Е. — Peter the Great St. Petersburg Polytechnic University. 29 Politechnicheskaya St., St. Petersburg, 195251, Russia. E-mail: [email protected]
СЕМЕНОВ Артем Семенович — кандидат физико-математических наук доцент Санкт-Петербургского политехнического университета Петра Великого. 195251, Россия, г. Санкт-Петербург, Политехническая ул., 29. E-mail: [email protected]
SEMENOV Artem S. - Peter the Great St. Petersburg Polytechnic University. 29 Politechnicheskaya St., St. Petersburg, 195251, Russia. E-mail: [email protected]
© Санкт-Петербургский политехнический университет Петра Великого, 2016
