Научная статья на тему 'Особенности напряженно-деформированного состояния тонкостенного кольцевого образца с надрезом при внецентренном растяжении'

Особенности напряженно-деформированного состояния тонкостенного кольцевого образца с надрезом при внецентренном растяжении Текст научной статьи по специальности «Физика»

CC BY
226
26
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ТОНКОСТЕННЫЙ КОЛЬЦЕВОЙ ОБРАЗЕЦ С НАДРЕЗОМ / ВНЕЦЕНТРЕННОЕ РАСТЯЖЕНИЕ / НАПРЯЖЕННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ / КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ / THIN-WALLED CIRCULAR SINGLE-EDGE-NOTCH SPECIMEN / ECCENTRIC TENSION / STRESS-STRAIN STATE / COMPUTER SIMULATION

Аннотация научной статьи по физике, автор научной работы — Коновалов В., Козлов А.В., Коновалов Д.А., Партин А.С.

Характеристики деформационного поведения и сопротивления разрушению материала оболочки тепловыделяющих элементов (твэлов) используются при проектных обоснованиях их безопасной работы. Для определения данных характеристик в АО «Институт реакторных материалов» применяют испытание на статическое внецентренное растяжение тонкостенных кольцевых образцов с концентратором напряжений в виде надреза, вырезанных из оболочечных труб. Испытания проводят через определенные промежутки времени эксплуатации твэла. Получаемые в экспериментах зависимости силы нагружения образца от величины перемещения захвата испытательной машины имеют различный вид, зависящий от накопленной поврежденности оболочки твэла. Причины поведения кривой нагружения в рассматриваемом виде испытания можно установить на основе его компьютерного моделирования. В работе приведена схема испытания. Описана постановка упругопластической задачи расчета напряженно-деформированного состояния кольцевого образца с надрезом при его внецентренном растяжении. Рассмотрено плоское деформированное состояние. Моделирование выполнено методом конечных элементов по разработанной в ИМАШ УрО РАН для рассматриваемого вида испытания специализированной компьютерной программе. Нагрузка прикладывается малыми шагами в приращениях перемещения захвата испытательной машины. На каждом шаге нагрузки расчетный алгоритм основан на принципе виртуальной мощности в приращениях. Использованы разработанные авторами определяющие соотношения для упругопластической среды с большими пластическими деформациями. Приведены результаты компьютерного моделирования процесса нагружения образца из аустенитной стали ЧС86хд. Установлено, что кривая зависимости силы растяжения образца от перемещения захвата испытательной машины состоит из двух возрастающих прямолинейных участков и соединяющей их закругленной части. Угол наклона к горизонтальной оси у первого прямолинейного участка кривой существенно больше, чем у второго. Проанализировано изменение накопленной пластической деформации в области надреза образца в разных точках кривой нагружения и объяснены особенности поведения кривой нагружения образца на разных ее участках, в том числе с учетом образования и развития макротрещины в вершине надреза образца. Представлены поля распределений показателей напряженного состояния в области надреза образца.

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — Коновалов В., Козлов А.В., Коновалов Д.А., Партин А.С.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

FEATURES OF STRESS-STRAIN STATE OF THIN-WALLED CIRCULAR SINGLE-EDGE-NOTCH SPECIMEN UNDER ECCENTRIC TENSION

Features of deformation behavior and fracture resistance of a material of fuel element claddings are used to design and ensure their safe operation. To determine these features, JSC “Institute of Nuclear Materials” applies static eccentric tension tests of thin-walled circular specimens with a stress concentrator in the form of a single-edge-notch, which are cut from shell pipes. The tests were made in certain time intervals of fuel element cladding operation. The experiments show various dependences of the specimen loading force on the value of the test machine grip displacement. These dependences are influenced by the accumulated damage of a fuel element cladding. Loading curve behavior is explained by the test simulation. The paper presents the test scheme. We formulated the elastoplastic problem of stress-strain state simulation of a single-edge-notch circular specimen under eccentric tension. The plane strain state is considered. The finite element simulation is made using the specialized computer program for the considered type of tests. The program is developed in the Institute of Engineering Science, Ural Branch of the Russian Academy of Sciences. The load is applied gradually in small increments of the test machine grip motion. At each load step, the simulation algorithm is based on the principle of virtual power. Constitutive equations developed by the authors for the elastoplastic medium with large plastic deformations are used. The simulation results of loading processes are shown for the specimens made from the austenitic steel ChS86hd. It appears that the curve describing the dependence of tensile force in the specimen on the test machine grip displacement consists of two straight line segments with a smooth transition. The first straight segment is notoriously more inclined to horizontal axis compared to the second one. The change of the accumulated plastic strain in the region of the specimen notch in different points of the loading curve was analyzed. The features of the specimen at different points of the loading curve were explained taking into account the formation and development of a macrocrack at the top of the specimen notch. Fields of stress state indicator distributions in the region of the specimen notch are shown.

Текст научной работы на тему «Особенности напряженно-деформированного состояния тонкостенного кольцевого образца с надрезом при внецентренном растяжении»

Коновалов А.В., Козлов А.В., Коновалов Д.А., Партин А.С. Особенности напряженно-деформированного состояния тонкостенного кольцевого образца с надрезом при внецентренном растяжении // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2018. - № 3. - С. 44-52. DOI: 10.15593/perm.mech/2018.3.05

Konovalov A.V., Kozlov A.V., Konovalov D.A., Partin A.S. Features of stress-strain state of thin-walled circular single-edge-notch specimen under eccentric tension. PNRPUMechanics Bulletin, 2018, no. 3, pp. 44-52. DOI: 10.15593/perm.mech/2018.3.05

ВЕСТНИК ПНИПУ. МЕХАНИКА

№ 3, 2018 PNRPU MECHANICS BULLETIN

http://vestnik.pstu.ru/mechanics/about/inf/

Б01: 10.15593/регш.шееЬ/2018.3.05 УДК 539.3

ОСОБЕННОСТИ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ТОНКОСТЕННОГО КОЛЬЦЕВОГО ОБРАЗЦА С НАДРЕЗОМ ПРИ ВНЕЦЕНТРЕННОМ РАСТЯЖЕНИИН^^Н!

Л О 11

А.В. Коновалов1, А.В. Козлов2, Д.А. Коновалов1, А.С. Партин1

''Институт машиноведения УрО РАН (ИМАШ УрО РАН), Екатеринбург, Россия

2АО «Институт реакторных материалов» (АО «ИРМ»), Заречный, Свердловская область, Россия

О СТАТЬЕ АННОТАЦИЯ

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

В работе приведена схема испытания. Описана постановка упругопластической задачи расчета напряженно-деформированного состояния кольцевого образца с надрезом при его внецентренном растяжении. Рассмотрено плоское деформированное состояние. Моделирование выполнено методом конечных элементов по разработанной в ИМАШ УрО РАН для рассматриваемого вида испытания специализированной компьютерной программе. Нагрузка прикладывается малыми шагами в приращениях перемещения захвата испытательной машины. На каждом шаге нагрузки расчетный алгоритм основан на принципе виртуальной мощности в приращениях. Использованы разработанные авторами определяющие соотношения для упругопластической среды с большими пластическими деформациями.

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

©ПНИПУ

© Коновалов Анатолий Владимирович - д.т.н., проф., e-mail: avk@imach.uran.ru, : 0000-0001 -9131 -8636.

Козлов Александр Владимирович - д.т.н., e-mail: kozlov_av@irmatom.ru, : 0000-0003-4138-2036.

Коновалов Дмитрий Анатольевич - к.т.н., н.с., e-mail: satterkein@yandex.ru, : 0000-0002-5306-7765.

Партин Александр Сергеевич - к.т.н., с.н.с., e-mail: lmd@imach.uran.ru, : 0000-0002-8359-2579.

Anatoly V. Konovalov - Doctor of Technical Sciences, Professor, e-mail: avk@imach.uran.ru, : 0000-0001 -9131 -8636.

Alexander V. Kozlov - Doctor of Technical Sciences, e-mail: kozlov_av@irmatom.ru, : 0000-0003-4138-2036.

Dmitry A. Konovalov - CSc in Technical Sciences, Researcher, e-mail: satterkein@yandex.ru, : 0000-0002-5306-7765.

Alexander S. Partin - CSc in Technical Sciences, Senior Researcher, e-mail: lmd@imach.uran.ru, ORCiD: 0000-0002-8359-2579.

I —^—— ——I Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International | (cc) CD (j9 I License (CC BY-NC 4.0)

I^^^QjHU^H This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)

Получена: 08 июня 2018 г. Принята: 03 сентября 2018 г. Опубликована: 08 октября 20'8 г.

Ключевые слова:

тонкостенный кольцевой образец с надрезом, внецентренное растяжение, напряженно-деформированное состояние, компьютерное моделирование.

FEATURES OF STRESS-STRAIN STATE OF THIN-WALLED CIRCULAR SINGLE-EDGE-NOTCH SPECIMEN UNDER ECCENTRIC TENSION

A.V. Konovalov1, A.V. Kozlov2, D.A. Konovalov1, A.S. Partin1

Institute of Engineering Science Ural Branch of Russian Academy of Sciences (IES UB RAS), Ekaterinburg, Russian Federation Joint stock company "Institute of Nuclear Materials" (JSC "INM"), Zarechny, Sverdlovsk region, Russian Federation

ARTICLE INFO

ABSTRACT

Received: 08 June 2018 Accepted: 03 September 2018 Published: 08 October 2018

Keywords:

thin-walled circular single-edge-notch specimen, eccentric tension, stressstrain state, computer simulation

Features of deformation behavior and fracture resistance of a material of fuel element claddings are used to design and ensure their safe operation. To determine these features, JSC "Institute of Nuclear Materials" applies static eccentric tension tests of thin-walled circular specimens with a stress concentrator in the form of a single-edge-notch, which are cut from shell pipes. The tests were made in certain time intervals of fuel element cladding operation. The experiments show various dependences of the specimen loading force on the value of the test machine grip displacement. These dependences are influenced by the accumulated damage of a fuel element cladding. Loading curve behavior is explained by the test simulation. The paper presents the test scheme. We formulated the elastoplastic problem of stress-strain state simulation of a single-edge-notch circular specimen under eccentric tension. The plane strain state is considered. The finite element simulation is made using the specialized computer program for the considered type of tests. The program is developed in the Institute of Engineering Science, Ural Branch of the Russian Academy of Sciences. The load is applied gradually in small increments of the test machine grip motion. At each load step, the simulation algorithm is based on the principle of virtual power. Constitutive equations developed by the authors for the elastoplastic medium with large plastic deformations are used. The simulation results of loading processes are shown for the specimens made from the austenitic steel ChS86hd. It appears that the curve describing the dependence of tensile force in the specimen on the test machine grip displacement consists of two straight line segments with a smooth transition. The first straight segment is notoriously more inclined to horizontal axis compared to the second one. The change of the accumulated plastic strain in the region of the specimen notch in different points of the loading curve was analyzed. The features of the specimen at different points of the loading curve were explained taking into account the formation and development of a macrocrack at the top of the specimen notch. Fields of stress state indicator distributions in the region of the specimen notch are shown.

©PNRPU

Введение

При эксплуатации в реакторах на тепловых и быстрых нейтронах происходит распухание и охрупчивание материалов внутриреакторных конструкций, в частности аустенитных сталей [1, 2], что негативно влияет на их работоспособность [3]. Под действием реакторного облучения изменяются также трещиностойкость [4], ползучесть, усталостная прочность [5, 6] и другие механические характеристики [7]. Наиболее сильным изменениям подвергаются оболочки тепловыделяющих элементов (твэлов) ядерных реакторов [8], испытывающие наибольшую дозовую нагрузку. Оболочки твэлов являются тонкостенными трубами малого диаметра, которые изготавливаются из различных сплавов: алюминия, циркония, сталей [9-11], в частности из аустенитной стали ЧС86хд [12]. Для подтверждения безопасности эксплуатации твэлов реакторов на быстрых нейтронах проводятся их послереакторные исследования [13], важной составной частью которых являются механические испытания изготовленных из оболочек твэлов образцов различных видов [14-16]. В АО «ИРМ» разработана и используется методика испытаний трубчатых образцов внутренним давлением твердого пластичного

заполнителя [17]. В ИМАШ УрО РАН совместно с АО «ИРМ» была разработана методика определения кривой упрочнения оболочек твэлов [18]. Получаемые в этих испытаниях характеристики констатируют механическое состояние оболочки, но не позволяют прогнозировать процесс деформации и разрушения полномасштабной оболочки в реальных условиях эксплуатации. Для этого надо определить характеристики сопротивления разрушению при наличии дефектов, являющихся концентраторами напряжений и деформаций.

Для определения характеристик сопротивления разрушению оболочек твэлов в АО «ИРМ» разрабатывается методика испытаний на статическое внецентренное растяжение тонкостенных кольцевых образцов с концентратором напряжений в виде надреза, вырезанных из оболочечных труб. Надрез представляет собой трещину с радиусом в ее вершине. Испытания проводят через определенные промежутки времени эксплуатации твэ-ла. Получаемые в экспериментах зависимости силы на-гружения Е образца от величины перемещения и захвата испытательной машины имеют различный вид (рис. 1), зависящий от накопленной поврежденности оболочки твэла. Кривые имеют начальный линейный участок изменения силы нагружения, после которого происходит

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

Для данного вида испытания классические характеристики сопротивления разрушению, разработанные в механике распространения трещины, такие как критический коэффициент интенсивности напряжений К, сила сопротивления раскрытию трещины ^, ./-интеграл, критическое раскрытие трещины 8С [19-21], неприемлемы [22], так как надрез отличается по форме от острой трещины, образец тонкостенный, пластическая зона у вершины надреза, как будет показано ниже, значительна, напряженно-деформированное состояние в области надреза в образце имеет сложный вид и существенно изменяется с изменением геометрии образца в процессе нагружения.

Рис. 1. Примеры экспериментальных кривых внецентренного растяжения тонкостенного кольцевого образца с надрезом из стали ЧС86хд Fig. 1. Examples of experimental curves for eccentric tension of a thin-walled circular single-edge-notch specimen made from steel ChS86hd

Характеристики сопротивления разрушению оболочек твэлов следует разрабатывать с позиций накопления по-врежденности, основываясь на аппарате механики повре-жденности [23-25], согласно которому вводится мера по-врежденности материала, как частный вариант скалярная величина ю , принимающая значение от 0 до 1. При значении ю , равном 1, в материальной точке среды образуется макротрещина. Экспериментально доказано [26-28], что пластичность металлов и сплавов зависит от безразмерных характеристик напряженного состояния, построенных на базе трех инвариантов тензора напряжений, например показателя напряженного состояния ст0/Т и параметра Ло-

де-Надаи = 2--1, где ст0 - среднее нормальное

напряжение; Т - интенсивность касательных напряжений, ^, ст2, ст3 - главные напряжения.

В литературе предложено много феноменологических математических моделей накопления поврежденно-сти [24, 29-32]. Детальный и достаточно полный обзор таких моделей сделан в работе [33]. Аргументами феноменологических моделей накопления поврежденности в материальной точке среды являются либо величина предельной пластичности как функция от показателей напряженного состояния, либо некоторая величина, зависящая от инвариантов тензора напряжений. Для получения информации о предельной пластичности проводят несколько видов испытаний, в которых варьируются значения выбранных характеристик напряженного состояния [26-28]. В работе [34] представлен обзор подходов и методов получения данной информации, а также описан разработанный комплекс испытаний, позволяющий без применения техники высоких давлений определять предельную пластичность металлов и сплавов в условиях горячей деформации в широком интервале изменений показателей напряженного состояния.

Объектом исследования данной работы являются оболочки твэлов реакторов на быстрых нейтронах из стали ЧС86хд. Целью работы является анализ изменения напряженно-деформированного состояния тонкостенного кольцевого образца с надрезом из стали ЧС86хд при его внецентренном растяжении на основе конечно-элементного моделирования. Результаты данного анализа могут служить основой для выбора схем проведения испытаний на предельную пластичность рассматриваемой стали. Проведенный информационный поиск в электронных библиотеках ScienceDirect и eLIBRARY с глубиной поиска 25 лет не выявил примеров моделирования рассматриваемого процесса.

1. Постановка задачи

На рис. 2 представлена расчетная схема испытания. В силу симметрии схемы нагружения рассмотрена половина тонкостенного кольцевого образца. Нагружающая сила Е приложена в точке Б, перемещение и данной точки равно половине суммарного перемещения и захватов испытательной машины. Здесь Я1 и Я2 - внешний и внутренний радиусы кольцевого образца; Ж -ширина стенки кольца; а - ширина стенки кольца в месте надреза; С0 - начальное (до нагружения) расстояние от оси симметрии до точки Б. Образец имеет угловой надрез с радиусом в вершине надреза.

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

СТ1

Рис. 2. Расчетная схема нагружения кольцевого образца с надрезом в испытании на внецентренное растяжение Fig. 2. The loading scheme of a thin-walled circular single-edge-notch specimen in the eccentric tension test

Все расчеты напряженно-деформированного состояния кольцевого образца при его внецентренном растяжении выполнены методом конечных элементов с использованием компьютерной программы, разработанной в ИМАШ УрО РАН. Рассматривалась упруго-пластическая постановка задачи для плоского деформированного состояния.

Нагрузку задавали мелкими шагами. На каждом малом шаге приращения нагрузки на интервале времени [t, t + At] алгоритм вычислений основан на принципе

виртуальной мощности в приращениях, уравнение которого имеет вид

J о••VhdV-J P • hd E-

+м-а«р • ь-^о. (1)

dt V у е )

Здесь с - тензор напряжений Коши; Р - плотность поверхностных сил; Ь - вариация кинематически допустимых полей скоростей перемещений; V - набла-оператор; У ,2 - объем и поверхность деформируемого тела в актуальной конфигурации соответственно; -У , - 2 - элементы объема и площади поверхности деформируемого тела в актуальной конфигурации соответственно; точкой и двумя точками обозначены соответственно скалярное и двойное скалярное произведение тензоров.

Дифференцируя по времени интегралы в уравнении (1) по подвижным объему и поверхности и разделив полученное равенство на At, имеем

J («т + V • V - о • Vv) • • VhrfF -

V

- J [P + (V • v)P - (n • Vv • n)P] • ЫЕ + J о ••VhdV -J P • hd eV 0,

+ 1 d At dt v v

(2)

где V - вектор скорости перемещения материальной точки; Vv - градиент скорости перемещений; п - единичная внешняя нормаль в точках поверхности Е , точкой сверху обозначена полная производная по времени.

Равенство, аналогичное равенству (2), для решения краевой задачи, названное обновляемой Лагранжевой формулировкой (updated Lagrangian formulation), получено с использованием второго тензора напряжений Пиола-Кирхгоффа в работе [35]. Для первого тензора напряжений Пиола-Кирхгоффа оно приведено в работе [36]. Формула для скоростной составляющей равенства (1) для тензора напряжений Коши выведена в работе [37].

В равенстве (2) использовали следующие определяющие соотношения, полученные в работе [38]:

o = X©I + 2(X© + |a)D-Vvo-o-Vvr-J|a^S, (3)

S = о -ct0I

ст0 = K©, K = Х + уъ ц,

(4)

D = 0,5(Vv + VV); условие текучести Мизеса

Y(S, к) = Т2 - к2 = 0.

Здесь X, |i - коэффициенты Ламе; I - единичный тензор; D - тензор скоростей деформаций; S - девиатор тензора напряжений; ст0 - среднее нормальное напряжение, K - объемный модуль упругости; © - относительное изменение индивидуального объема бесконечно малой частицы среды; Т = *J0,5S --S - интенсивность касательных напряжений; к - напряжение текучести при сдвиге; J - числовой коэффициент, принимающий

значения J =0 при < 0 или S • -S < 0 при '1'= 0 и J =\ при Т= 0. S-S >0; х - параметр упрочнения,

X = ер , еР - компоненты тензора скоростей пластических деформаций в декартовой системе координат. Равенство (3) получено c использованием корота-ционной производной Коттер и Ривлина к тензору напряжений Коши, применение которой обосновано в работе [38].

На каждом малом шаге нагрузки к равенству (2) с подставленными в него определяющими соотношениями (3) применяли конечно-элементную аппроксимацию [39] и решали полученную систему линейных алгебраических уравнений, в которой неизвестными являлись компоненты вектора скорости в узлах сетки конечных элементов. В расчетах использовали четырехугольные изопараметрические элементы, образующие регулярную сетку. Скорость параметра упрочнения % в равенстве (3) определяли итерационно методом начальных деформаций [40, 41]. На каждой такой итерации условие текучести Мизеса (4) удовлетворяли также итерационно методом Ньютона с использованием обобщенного трапецеидального правила (generalized trapezoidal rule) [42].

Шаг нагрузки в виде малого перемещения Ли рассчитывался в программе из условия обеспечения в начальный момент нагружения максимальной по всем конечным элементам величины интенсивности напря-

V

E

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

2. Результаты компьютерного моделирования

В расчетах для материала образца из стали ЧС86хд задавали модуль Юнга и коэффициент Пуассона, равными соответственно 187 ГПА и 0,2. Напряжение текучести при нагружении вычисляли по формуле

к = а0(1 + "х)"2,

где " = 439 МПа; " =14,86; " = 0,216.

В расчетах задали следующие применяемые в экспериментах размеры (в мм) кольцевого образца: Я-[ = 3,05, Я2 = 3,45, С0 = 1,6, Ж = 0,4, а = 0,28, ширина кольца 4,15 мм, радиус в вершине надреза 0,03 мм, угол надреза 30°.

На рис. 3 приведена рассчитанная кривая изменения силы растяжения Е кольца от перемещения захвата и. Кривая имеет два возрастающих прямолинейных участка и соединяющую их закругленную часть. Угол наклона к горизонтальной оси у первого прямолинейного участка кривой существенно больше, чем у второго. На данной кривой выделили пять точек, в которых проанализировали напряженно-деформированное состояние кольцевого образца при его внецентренном растяжении. Перемещение захвата и в этих точках равняется соответственно 0,1, 0,26, 0,51, 1,2 и 2 мм. На рис. 4 для этих точек представлены в области надреза образца распределения значений меры пластиче-

!

ской деформации %(/) = |х(т)с/т, которая является

о

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

Рис. 3. Рассчитанное изменение силы растяжения F кольцевого образца от перемещения захвата u Fig. 3. The calculated change of the tensile force F in the circular specimen dependent on the grip displacement u

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

Как видно из рис. 4 и 5, деформированное и напряженное состояния в зоне надреза кольца существенно неоднородны. Самые большие значения степени пластической деформации сдвига % имеют место непосредственно в вершине надреза (рис. 5, а). В нижней половине по толщине стенки кольца преобладают растягивающие напряжения, в верхней - сжимающие. Перед радиусом выреза на некотором расстоянии от него выделяется область с наибольшими растягивающими напряжениями и положительными значениями показателя напряженного состояния ст0/Т, лежащими в пределах 2,14-2,75 при u = 2 (рис. 5, б). При охвате

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

%

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

0.000 -0.027 0.027-0.054 0.051-0.081 0.081 -0.107 0.107-0.134 0.134-0.161 0.161 -0.188 0.188-0.215

Рис. 4. Распределение степени пластической деформации сдвига х при перемещениях захвата u, равных 0,26 мм (а) и 0,51 мм (б) (соответствует точкам 2 и 3 на рис. 3) Fig. 4. The distribution of plastic shear strain х at the grip displacement u, equal to 0.26 mm (a) and 0.51 mm (b). This corresponds to points 2 and 3 in Fig. 3

V

%

WO- 0.258

ao/T

-2 12 - -1 51

99 --0 74

Рис. 5. Распределение степени пластической деформации сдвига х (а), показателя напряженного состояния ст0/Т (б) и параметра Надаи-Лоде (в) в области пластической деформации при u = 2 мм (точка 5 на рис. 3) Fig. 5. Distribution of plastic shear strain х (а), stress state indicator а0/Т (b) and Nadai-Lode stress state parameter (c) in plastic deformation region at u=2 mm (point 5 in Fig. 3)

напряженное состояние в этой зоне становится близким к состоянию при чистом сдвиге, так как значения параметра Лоде-Надаи близко к нулю (см. рис. 5, в).

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

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

Заключение

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

Анализ выполненных расчетов показал следующее:

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

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

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

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

- перед радиусом выреза на некотором расстоянии от него выделяется область с наибольшими растягивающими напряжениями и положительными значениями показателя напряженного состояния а0/Т, что определяет рост поврежденности материала образца в этой части и образование и развитие макротрещины;

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

б

а

в

48904853485353485348535323

- отклонение после перегиба от монотонного возрастания и понижение экспериментальной кривой нагрузки свидетельствует о начале образования макротрещины в зоне его надреза и ее развитии;

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

Библиографический список

1. К вопросу о радиационном распухании и радиационном охрупчивании аустенитных сталей. Ч. I. Экспериментальные результаты / Б.З Марголин [и др.] // Вопросы материаловедения. - 2009. - № 2 (58). - С. 89-98.

2. Neustroev V.S., Garner F.A. Very high swelling and embrittlement observed in a Fe-18Cr-10Ni-Ti hexagonal fuel wrapper irradiated in the BOR-6O fast reactor // Fusion Materials. -2007. -Vol. 43 - P. 109-122.

3. Debusson P.H., Seran J.L., Soulat P. Irradiation embrittlement of reactor internals materials // Effect of Irradiation on Water Reactor Internals. - Paris, 1997. - ECSC-EEC-EAEC: Brussels-Luxembourg. -1997. - Vol. 3. - P. 30-34.

4. Bates F., Gilbert E.R. Effect of stress on swelling in 316 stainless steel // J. Nucl. Mat. - 1978. - Vol. 71 - P. 286-292.

5. Creep-fatigue evolution of normalized and tempered modified 9Cr-1Mo / K. Aoto, R. Komine, F. Ueno, H. Kavasaki, Y. Wada // Nuclear Engineering and Design. - 1994. - Vol. 153. -No. 1. - P. 97-110.

6. Creep-fatigue evolution method for Mod.9Cr-1Mo weldment / T. Asayana, S. Hasebe, Y. Heakawa, Y. Wada // Proceeding of 12lthSMiRT. - 1993. - Vol. 1. - P. 123.

7. Evaluation procedures for irradiation effects and sodium environmental effects for the structural design of Japanese fast breeder reactors / T. Asayana, Y. Abe, N. Miyaji, M. Koi, E. Yoshida // ASME J. Pressure Vessel Technol. - 2001. - Vol. 123. - P. 49.

8. Garner F.A. Irradiation Performance of Cladding and Structural Steels in Liquid Metal Reactors // Material Science and Technology: A Comprehensive Treatment. - 1994. - Vol. 10. -P. 419-543.

9. Герасимов В.В., Монахов А.С. Материалы ядерной техники: учебник для вузов - 2-е изд., перераб. и доп. - М.: Энергоатомиздат, 1982. - 288 с.

10. Бойко В.И., Кошелев Ф.П. Ядерные технологии в различных сферах человеческой деятельности: учеб. пособие. -Томск: Изд-во ТПУ, 2006. - 342 с.

11. Требования к конструкционным материалам. Тепловыделяющие элементы ядерных реакторов [Электронный ресурс]. -URL: http://leg.co.ua/arhiv/generaciya/teplovydelyayuschie-elementy-yadernyh-reaktorov-26.html (дата обращения: 07.06.2018).

12. Разработка конструкционной стали для твэлов и ТВС быстрых натриевых реакторов / А.В. Целищев, В.С. Агеев, Ю.П. Буданов [и др.] // Атомная энергия. - 2010. - Т. 108, № 4. - С. 217-221.

13. Оборудование и методики послереакторных исследований в комплексе защитных камер АО «ИРМ» / Л.П. Синельников [и др.] // Атомная энергия. - 2016. - Т. 121, № 4. -С. 187-194.

14. СТП 086-288-99. Оболочки облученных и необлучен-ных твэлов. Методика выполнения измерений характеристик механических свойств при растяжении кольцевых образцов из

Благодарности

Работа выполнена по бюджетной теме № 03912016-0001 (г/р АААА-А18-118020790140-5).

Acknowledgments

The work was performed according to the state assignment No. 0391-2016-0001 (s/r AAAA-A18-118020790140-5)

тонкостенных труб из стали и сплавов / ГНЦ РФ «НИИАР». -М., 1999. - 18 с.

15. ОИ 001.325-2006. Методика выполнения испытаний для определения характеристик механических свойств при растяжении образцов труб в поперечном направлении из сплавов циркония / АО «ВНИИНМ». - М., 2006. - 30 с.

16. Grigoriev V., Jacobsson R., Schrire D. Advanced techniques for mechanical testing of irradiated cladding materials // International Atomic Energy Agency, available at: http://www.iaea.org/ inis/collection/NCLCollectionStore/_Public/33/018/33018365.pdf (accessed 07 Juny 2018).

17. Свидетельство на полезную модель № 26129. Устройство для испытаний трубчатых образцов / Козлов А.В., Брюшко-ва С.В., Ковалев И.Н., Кузнецов Ю.Г. Зарегистр. Гос. реестре полезных моделей Рос. Федерации. - М., 10 ноября 2000 г.

18. Методика определения кривой упрочнения материала оболочек твэлов / А.В. Коновалов, Д.И. Вичужанин, А.С. Пар-тин, А.В. Козлов // Заводская лаборатория. Диагностика материалов. - 2017. - Т. 83, № 7. - C. 58-61.

19. Разрушение. Т.2. Математические основы теории разрушения / под ред. Г. Либовица. - М.: Мир, 1975. - 766 с.

20. Партон В.З., Морозов Е.М. Механика упругопласти-ческого разрушения. - М.: Наука, 1985. - 504 с.

21. Херцберг Р.В. Деформация и механика разрушения конструкционных материалов. - М.: Металлургия, 1989. - 576 с.

22. Работнов Ю.Н. Механика деформируемого твердого тела. - М.: Наука, 1988. - 712 с.

23. Качанов Л.М. Основы механики разрушения. - М.: Наука, 1974. - 312 с.

24. Колмогоров В.Л. Напряжения, деформации, разрушение. - М.: Металлургия, 1970. - 229 с.

25. Sumio Murakami. Continuum Damage Mechanics -A Continuum Mechanics Approach to the Analysis of Damage and Fracture. - Dordrecht, Heidelberg, London, New York: Springer, 2012. - 402 p.

26. Bridgman P.W. Studies in Large Plastic Flow and Fracture with Special Emphasis on the Effects of Hydrostatic Pressure. -New York: McGraw-Hill, 1952. - 362 p.

27. Some problems of large plastic deformation of metals at high pressure / B.I. Beresnev [et al.]. - New York: Pergamon press, 1963. - 79 p.

28. Колмогоров В.Л., Шишминцев В.Ф. Зависимость пластичности стали от гидростатического давления // Физика металлов и металловедение. - 1966. - Т. 21, № 6. - С. 910-914.

29. Смирнов-Аляев Г.А., Розенберг В.М. Теория пластических деформаций металлов: механика конечного формоизменения. - М.; Л.: Машгиз. Ленингр. отд-е, 1956. - 367 с.

30. McClintock F.A. A criterion of ductile fracture by the growth of holes // Journal of Applied Mechanics. - 1968. -Vol. 35. - P. 363-371. DOI: org/10.1115/1.3601204

31. Rice J.R., Tracey D.M. On the ductile enlargement of voids in triaxial stress fields // Journal of the Mechanics and Physics of Solids. - 1969. - Vol. 17. - P. 201-217. DOI: org/10.1016/0022-5096(69)90033-7

32. Hancock J.W., Mackenzie A.C. On the mechanisms of ductile failure in high-strength steels subjected to multi-axial stress-states // Journal ofthe Mechanics and Physics of Solids. - 1976. - Vol. 24. - No. 2-3. - P. 147-160. DOI: org/10.1016/0022-5096(76)90024-7

33. Trong Son Cao. Models for ductile damage and fracture prediction in cold bulk metal forming processes: a review // International Journal of Material Forming. - 2017. - Vol. 10. - No. 2. -P. 139-171. DOI 10.1007/s12289-015-1262-7

34. Смирнов С.В., Вичужанин Д.И., Нестеренко А.В. Комплекс испытаний для исследования влияния напряженного состояния на предельную пластичность металла при повышенной температуре // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2015. -№ 3. - С. 146-164. DOI: 10.15593/perm.mech/2015.3.11

35. Bathe K.J., Ramm E., Wilson E.L. Finite element formulations for large deformation dynamic analysis // International Journal for Numerical Methods in Engineering. - 1975. - Vol. 9. -P. 353-386.

References

1. Margolin B.Z., Kursevich I.P., Sorokin A.A., Lapin A.N., Kokhanov V.I., Neustroev V.S. K voprosu o radiatsionnom raspukhanii i radiatsionnom okhrupchivanii austenitnykh stalei. Chast' I. Eksperimental'nye rezul'taty [On the issue of radiation swelling and radiation embrittlement of austenitic steels. Part I. Experimental results]. Voprosy materialovedeniya - Material science issues, 2009, vol. 2, no. 58, pp.89-98.

2. Neustroev V.S., Garner F.A. Very high swelling and embrittlement observed in a Fe-18Cr-10Ni-Ti hexagonal fuel wrapper irradiated in the BOR-60 fast reactor. Fusion Materials, 2007, vol. 43, pp. 109-122.

3. Debusson P.H., Seran J.L., Soulat P. Irradiation embrittlement of reactor internals materials. "Effect of irradiation on water reactor internals", Paris, 1997. - ECSC-EEC-EAEC: Brussels-Luxembourg, 1997, vol. 3, pp. 30-34.

4. Bates F., Gilbert E.R. Effect of stress on swellingin 316 stainless steel. J.Nucl. Mat., 1978, vol. 71, pp. 286-292.

5. Aoto K., Komine R., Ueno F., Kavasaki H., Wada Y. Creep-fatigue evolution of normalized and tempered modified 9Cr-1Mo. Nuclear Engineering and Design, 1994, vol. 153, no. 1, pp. 97-110.

6. Asayana T., Hasebe S., Heakawa Y., Wada Y. Creep-fatigue evolution method for Mod.9Cr-1Mo weldment. Proceeding of 12lthSMiRT, 1993, vol. 1, p. 123.

7. Asayana T., Abe Y., Miyaji N., Koi M., Yoshida E. Evaluation procedures for irradiation effects and sodium environmental effects for the structural design of Japanese fast breeder reactors. ASME J. Pressure Vessel Technol., 2001, vol. 123, p. 49.

8. Garner F.A. Irradiation Performance of Cladding and Structural Steels in Liquid Metal Reactors. Material Science and Technology: A Comprehensive Treatment, 1994, vol. 10, pp. 419-543.

9. Gerasimov V.V., Monahov A.S. Materialy jadernoj tehniki: uchebnik dlja vuzov [Materials of nuclear technology]. Moscow, Jenergoatomizdat, 1982, 288 p.

10. Bojko V.I., Koshelev F.P. Jadernye tehnologii v razlichnyh sferah chelovecheskoj dejatel'nosti [Nuclear technologies in various spheres of human activity]. Tomsk, Izd-vo TPU, 2006, 342 p.

11. Trebovanija k konstrukcionnym materialam -Teplovydeljajushhie jelementy jadernyh reaktorov [Requirements

36. Cheng J.-H., Kikuchi N. An analysis of metal forming processes using large deformation elastic-plastic formulations // Computer Methods in Applied Mechanics and Engineering. -1985. - Vol. 49. - No. 1. - P. 71 - 108.

37. Поздеев А.А., Трусов П.В., Няшин Ю.И. Большие уп-ругопластические деформации: теория, алгоритмы, приложения. - М.: Наука, 1986. - 232 с.

38. Коновалов А. В. Определяющие соотношения для уп-ругопластической среды при больших пластических деформациях // Изв. РАН. Механика твердого тела. - 1997. - № 5. -С. 139-147.

39. Зенкевич О. Метод конечных элементов в технике. -М.: Мир, 1975. - 539 с.

40. Ильюшин А.А. Пластичность. - М.: Гостехиздат, 1948. - 376 с.

41. Биргер И.А. Некоторые общие методы решения задач теории пластичности // ПММ. - 1951. - Т. XV. - Вып. 6. -С. 765-770.

42. Ortiz M., Popov E.P. Accuracy and stability of integration algorithms for elastoplastic constitutive relations // International Journal for Numerical Methods in Engineering. - 1985. - Vol. 21. -P. 1561-1576.

for structural materials - Fuel element claddings of nuclear reactor], available at: http://leg.co.ua/arhiv/generaciya/teplovydelya-yuschie-elementy-yadernyh-reaktorov-26.html.

12. Celishhev A.V., Ageev V.S., Budanov Ju.P., et al. Razrabotka konstrukcionnoj stali dlja tvjelov i TVS bystryh natrievyh reaktorov [Development of structural steel for fuel element claddings and fuel assemblies of fast sodium reactors]. Atomnaja jenergija - Nuclear energy, 2010, vol. 108, no. 4, pp. 217-221.

13. Sinel'nikov L.P., Averin S.A., Kozlov A.V., Shcherbakov E.N., Kinev E.A. Oborudovanie i metodiki poslereaktornykh issledovanii v komplekse zashchitnykh kamer AO «IRM» [Equipment and methods of post-reactor research in the complex of protective chambers of JSC «IRM»]. Atomnaya energiya - Nuclear energy, 2016, vol. 121, no. 4, pp. 187-194.

14. STP 086-288-99. Obolochki obluchennykh i neobluchennykh tvelov. Metodika vypolneniya izmerenii kharakteristik mekhanicheskikh svoistv pri rastyazhenii kol'tsevykh obraztsov iz tonkostennykh trub iz stali i splavov [Shells of irradiated and unirradiated fuel rods. Methodology of measurement characteristics of mechanical tensile properties of ring specimens made of thin-walled tubes and pipes of steel and alloys]. GNTs RF «NIIAR», 1999, 18 p.

15. OI 001.325-2006. Metodika vypolneniya ispytanii dlya opredeleniya kharakteristik mekhanicheskikh svoistv pri rastyazhenii obraztsov trub v poperechnom napravlenii iz splavov tsirkoniya [Test procedure for determining the characteristics of mechanical properties of tensile pipe samples in the transverse direction of zirconium alloys]. AO «VNIINM», 2006, 30 p.

16. Grigoriev V., Jacobsson R., Schrire D. Advanced techniques for mechanical testing of irradiated cladding materials, International Atomic Energy Agency. available at http://www.iaea.org/ inis/collection/NCLCollectionStore/_Public/33/018/33018365.pdf

17. Kozlov A.V., Bryushkova S.V., Kovalev I.N., Kuznetsov Yu.G. Ustroistvo dlya ispytanii trubchatykh obraztsov. Svidetel'stvo na poleznuyu model' №26129, Zaregistrirovan v Gosudarstvennom reestre polznykh modelei Rossiiskoi federatsii [Device for testing tubular samples. Certificate for utility model

No. 26129, Registered in The state register of utility models of the Russian Federation]. Moscow, 10 November, 2000.

18. Konovalov A.V. Vichuzhanin D.I. Partin A.S., Kozlov A.V. Metodika opredeleniya krivoi uprochneniya materiala obolochek tvelov [Determination of True Stress-Strain (Hardening) Curve for the Fuel Rod Material]. Zavodskaya laboratoriya. Diagnostika materialov - Industrial laboratory. Diagnostics of materials, 2017, vol. 83, no. 7, pp. 58-61.

19. Fracture. Volume 2. Mathematical fundamentals. Edited by H. Liebowitz. Moscow, Mir, 1975, 766 p.

20. Parton V.Z., Morozov E.M. Mekhanika uprugoplastiches-kogo razrushenia [Mechanics of elastic-plastic fracture]. Moscow, Nauka, 1985, 504 p.

21. Hertzberg R.V. Deformatsiya i mekhanika razrusheniya konstruktsionnykh materialov [Deformation and fracture mechanics of engineering materials]. Moscow, Metallurgiya, 1989, 576 p.

22. Rabotnov Yu. N. Mekhanika deformiruemogo tverdogo tela [Mechanics of a deformable solid]. Moscow, Nauka, 1988, 712 p.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

23. Kachanov L.M. Osnovi mekhaniki razrushenia [Foundations of fracture mechanics]. Moscow, Nauka, 1974, 312 p.

24. Kolmogorov V.L. Napriazheniia, deformatsii, razrushenie [Stresses, strains, destruction]. Moscow, Metallurgiia, 1970, 229 p.

25. Sumio Murakami. Continuum Damage Mechanics - A Continuum Mechanics Approach to the Analysis of Damage and Fracture. -Dordrecht, Heidelberg, London, New York, Springer, 2012, 402 p.

26. Bridgman P.W. Studies in Large Plastic Flow and Fracture with Special Emphasis on the Effects of Hydrostatic Pressure. New York, McGraw-Hill, 1952, 362 p.

27. Some problems of large plastic deformation of metals at high pressure / B.I. Beresnev [et al.]. - New York: Pergamon press, 1963. - 79 p.

28. Kolmogorov V.L., Shishmintsev V.F. Ductility of steel as a function of hydrostatic pressure. Physics of Metals and Metallography, 1966, vol. 21, no. 6, pp. 97-101.

29. Smirnov-Aliaev G.A., Rozenberg V.M. Teoriia plasticheskikh deformatsii metallov: mekhanika konechnogo formoizmeneniia [The theory of metal deformation: the final forming mechanics]. Moscow; Leningrad, Mashgiz, 1956, 367 p.

30. McClintock F.A. A criterion of ductile fracture by the growth of holes, Journal of Applied Mechanics,1968, vol. 35, pp. 363-371. doi: org/10.1115/1.3601204

31. Rice J.R., Tracey D.M. On the ductile enlargement of voids in triaxial stress fields. Journal of the Mechanics and Phys-

ics of Solids, 1969, vol. 17, pp. 201-217. doi: org/10.1016/0022-5096(69)90033-7

32. Hancock J.W., Mackenzie A.C. On the mechanisms of ductile failure in high-strength steels subjected to multi-axial stress-states, Journal of the Mechanics and Physics of Solids, 1976, vol. 24 (2-3), pp. 147-160. doi: org/10.1016/0022-5096(76)90024-7

33. Trong Son Cao. Models for ductile damage and fracture prediction in cold bulk metal forming processes: a review, International Journal of Material Forming, 2017, vol. 10, no. 2, pp. 139-171. DOI 10.1007/s12289-015-1262-7

34. Smirnov S.V., Vichuzhanin D.I., Nesterenko A.V. A set of tests for studying the effect of the stress state on ultimate metal plasticity at high temperature, PNRPU Mechanics Bulletin, 2015, no. 3. pр. 146-164. doi: 10.15593/perm.mech/2015.3.11

35. Bathe K.J., Ramm E., Wilson E.L. Finite element formulations for large deformation dynamic analysis, International journal for numerical methods in engineering, 1975, vol. 9, pp. 353-386.

36. Cheng J.H., Kikuchi N. An analysis of metal forming processes using large deformation elastic-plastic formulations, Computer methods in applied mechanics and engineering, 1985, vol. 49, no. 1, pp. 71-108.

37. Pozdeev A.A., Trusov P.V., Njashin Ju.I. Bol'shie uprugoplasticheskie deformacii: teorija, algoritmy, prilozhenija [Large elastic-plastic deformations: theory, algorithms, applications]. Moscow, Nauka, 1986, 232 p.

38. Konovalov A.V. Opredeljajushhie sootnoshenija dlja uprugoplasticheskoj sredy pri bol'shih plasticheskih deformacijah [Constitutive equations for elastic-plastic medium under large plastic deformations]. Izvestija RAN. Mehanika tverdogo tela -News of RAS. Solid mechanics, 1997, no. 5, pp. 139-147.

39. Zenkevich O. Metod konechnyh jelementov v tehnike [Finite element method in engineering]. Moscow, Mir, 1975, 539 p.

40. Il'jushin A.A. Plastichnost' [Plasticity]. Moscow, Gostehizdat, 1948, 376 p.

41. Birger I.A. Nekotorye obshhie metody reshenija zadach teorii plastichnosti [Some general methods for solving problems of plasticity theory]. PMM, 1951, vol. XV, no. 6, pp. 765-770.

42. Ortiz M., Popov E.P. Accuracy and stability of integration algorithms for elastoplastic constitutive relations, International journal for numerical methods in engineering, 1985, vol. 21, pp. 1561-1576.

i Надоели баннеры? Вы всегда можете отключить рекламу.