CC BY
22
DOI: 10.18721/JPM.10113 УДК: 539.3
Н.П. Кузнецов', Б.Е. Мельников2, А.С. Семенов2
1 Новгородский государственный университет им. Ярослава Мудрого,
г. Великий Новгород, Российская Федерация;
2 Санкт-Петербургский политехнический университет Петра Великого, г. Санкт-Петербург, Российская федерация
верификация и адаптация моделей пластичности при сложном нагружении с промежуточными полными
и частичными разгрузками
В целях более точного описания процессов сложного пассивного нагру-жения выполнены экспериментальные исследования упругопластического деформирования стальных трубчатых образцов при пропорциональном и непропорциональном (монотонном и циклическом) нагружениях, включающих частичные и промежуточные разгрузки. Проведено сравнение результатов экспериментов с прогнозами структурной (реологической) упругопластической модели и многоповерхностной теорией пластичности с одной активной поверхностью. Предложены модификации двух указанных выше моделей в целях повышения точности прогнозов.
ПЛАСТИЧНОСТЬ, ПАССИВНОЕ НАГРУЖЕНИЕ, ЭКСПЕРИМЕНТ, МОДЕЛИРОВАНИЕ.
Введение
В теории пластического течения, получившей в настоящее время наиболее широкое распространение, предполагается, что разгрузка материала происходит без изменения пластических деформаций. Однако результаты ряда исследований [1 — 19] показывают, что при пассивном нагружении (разгрузка, нагружение внутри поверхности нагружения и по касательной к ней) происходит накопление пластических деформаций. Эндохронная теория пластичности [4, 5] и обобщенная модель Прандтля (модель Мазинга) [6] обеспечивают описание пластического деформирования при пассивном нагружении. Формулировки определяющих уравнений указанных теорий ориентированы, прежде всего, на уточненное описание активного деформирования, а параметры материала определяются без учета специфики процессов при пассивном нагружении.
Экспериментальное изучение свойств поля пластических податливостей [20] послужило основой для разработки многопо-
верхностной теории пластичности с одной активной поверхностью [14 — 19]. Понятие пластических податливостей прямо или косвенно используется многоповерхностными теориями пластичности [21, 22]. Вместе с тем, накопленный экспериментальный материал по свойствам пластических податливостей материалов нельзя считать достаточным. Прежде всего, это касается деформирования при пассивных путях на-гружения. Накопление остаточных деформаций при таких нагружениях необходимо учитывать не только при расчете деформирования [14 — 19], но и при расчетах накопления повреждений [13, 23, 24], при исследовании процессов сверхпластичности [25], при оценке работоспособности и точности машинных агрегатов [26].
Настоящее исследование проведено для дальнейшего развития и обоснования концепции пластических податливостей и разработки методики выполнения расчетов, связанных со сложным переменным на-гружением, включающим промежуточные полные и частичные разгрузки.
определяющие уравнения для описания пассивного нагружения
Векторная форма представления напряжений и деформаций. При двухпараме-трическом нагружении с неизменными и совпадающими направлениями главных напряжений и деформаций пятимерные пространства девиаторов редуцируются до двумерных. Возможны различные равноправные подходы [7] для того, чтобы вводить координаты двумерного пространства, удовлетворяющие условию равенства длины вектора и интенсивности тензора. Одним из возможных вариантов задания координат и Х2 для точки, характеризующей напряженное состояние, является следующий:
1 3
^ = соэ ^ = < --(ст2 +СТ3) = -^
^ 2 Тз 2 (1)
12 = < аш^ = "у(ст2 -3з) = + 2S2),
где <31 — интенсивность напряжений по Мизесу; у — угол между вектором напряжений 2 и осью Х1; о1, о2, о3 — главные напряжения; S1, £2, $3 — главные значения де-виатора напряжений.
Аналогичные выражения вводятся и для пластических деформаций. они связывают
координаты Ер, Е2 соответствующей точки в двумерном пространстве пластических деформаций с главными значениями де-виатора пластических деформаций ер, ер следующими соотношениями:
Ер = ер соэ ф = е1,
Ер -
1
= ер эт ф = —¡= (ер
' Тз 1
(2)
где ер — интенсивность пластических деформаций, ф — угол между вектором пластических деформаций № и осью Ер.
В пространстве напряжений совокупностям напряженных состояний, объединенных общим признаком (развитие пластических деформаций заданной величины, разрушение и т. п.), соответствуют некоторые границы (гиперповерхности, в частном случае это гиперсферы в рассматриваемом векторном пространстве, линии — при двухосном нагружении). Одна из таких границ показана на рис. 1, Ь окружностью радиуса Са с центром в точке а.
Структурная (реологическая) модель. Рассмотрим структурную модель пластического деформирования материала [13, 15, 28] (рис. 1, а). Она предполагает, что в процессе нагружения граница текучести сохраняет свои форму и размер С (рис. 1, Ь).
Рис. 1. Структурная модель пластического деформирования материала (а) и схематическое представление движения поверхностей нагружения (Ь); Е1, Б2 — координаты пространства напряжений, Е1, Е2 — координаты пространства деформаций; X, У — координаты локальной системы координат, связанной с нагружением на последнем этапе; 1 — 3 — номера структурных элементов; а., р., у. — жесткости упругих элементов (координаты центров поверхностей), Са, Ср, С — пределы текучести пластических элементов (радиусы поверхностей); с — вектор
напряжений в двумерном пространстве напряжений
Ее положение определяется предысторией нагружения. Перемещение этой границы начинается при ее пересечении изнутри вектором напряжений.
Вектор приращения пластических деформаций ДЕр пропорционален проекции вектора приращения напряжений ДЕ на внешнюю нормаль к границе поверхности нагружения и направлен по этой нормали. Длины приращений рассматриваемых векторов |ДЕр = Дер и |ДЕ| = До. связаны соотношением
|ДЕР| = Н |ДЕ|СО8(ДЕлДЕр),
(3)
где На — величина модуля податливости, постоянная для всех точек границы.
При изменении напряженного состояния в пределах внутренней -окружности все элементы сухого трения неподвижны (рис. 1, Ь). При достижении границы элемент 1 начинает движение относительно элемента 2. В упругой связи между ними появляются усилия. Их можно представить натяжением упругих элементов а1 и а2. Усилия а1 и а2 являются координатами центра окружности радиуса Са в плоскости Х1, Х2 (рис. 1, Ь). Аналогично усилия в элементах в1 и Р2, у1 и у2 являются координатами центров окружностей радиусов Ср и С. Параметрам а1 и а2, Р1 и Р2, у1 и у2 соответствуют векторы а, в и у.
При нагружении по пути 1 ^ 2 (рис. 1, Ь) пластические деформации не развиваются, поскольку расстояние от текущей точки до центра окружности 5 меньше ее радиуса С :
а
8а =7 (^ -<*1)2 + (Е2 -а2)2 < С а. (4)
При дальнейшем движении вдоль прямой 1—2 при переходе через точку 2 (8а > Са) начинается движение элемента а по радиусу к -окружности, т. е. под углом
!1 - а1 а = агС^ — 1
Е2 - а2
(5)
который в общем случае не совпадает с заданным наклоном прямой, характеризующей изменение напряженного состояния (см. рис. 1, Ь). Элементы в и у остаются неподвижными пока выполняются условия
5а > Са, 8Р < С и пластические деформации соотношением
(
8у < Су. В этом случае определяются
ДЕр = Н (Е - а)
гаЛ /
1-Са
8а
(6)
При движении вдоль прямой 1 — 2 при переходе через точку 3 элемент в начинает движение (8а > Са, 8р > Ср, 8у < Су), и пластические деформации определяются выражением
( п Л
ДЕ = Н (Е - а)
пЛ 7
+ Нв(Е
в)
1 -Са
8аУ
( С Л 1-С
8в)
(7)
В общем случае, при движении всех трех окружностей 8а > Са, 8р > Ср, 8у > Су приращение пластических деформаций определяется равенством
ДЕ = Н (Е - а)
пЛ 7
(
1 -Са
8а
л
(
+ Н„(Е - в)
1 - С
8в
Л
(
+ ЩЕ - у)
С Л
(8)
которое, в отличие от предложенного в работе [13], вводится на основе дифференциальных зависимостей теории пластического течения.
Деформации на всем пути нагружения определяются суммированием приращений деформаций на каждом шаге. При разбиении пути нагружения на шаги целесообразно в качестве граничных точек выбирать те напряженные состояния, при которых в движение включается (или из него выключается) очередной элемент сухого трения структурной модели. Так, для пути, представленного на рис. 1, Ь, таковыми являются точки 2, 3 и 4.
Многоповерхностная теория с одной активной поверхностью. Математическая модель, основанная на предположении о невозможности чисто упругого деформирования как при нагрузке (увеличение о.), так и при разгрузке (уменьшение о.) была предложена в работе [20]. Авторами было принято, что соотношение (3) применимо в любой точке пути вторичного нагружения.
Утверждалось в качестве постулата, что направление вектора и модуль податливости Н являются функциями напряженного и деформированного состояний в конце первого шага нагружения (нулевого цикла) и текущей точки пути. В экспериментах с никелевыми и стальными тонкостенными трубками было установлено, что вектор ориентирован по нормали к окружности, проведенной через текущую точку и точку окончания первого нагружения, а величина Н во всех точках такой окружности примерно постоянна [14 — 20].
Пусть после нагружения по пути 0 ^ 1 (рис. 2) в точке 1 достигнуто напряженное состояние с интенсивностью напряжений <з;1. Затем, после излома траектории нагру-жения в точке 1, следует движение по пути 1 ^ 3. Направление приращения пластических деформаций совпадает с радиусом 2'2 окружности, проведенной через точки 1 и 2 из центра, лежащего на луче 0 — 1. Модуль вектора ношением
ДE р\ = Дер определяется соот-
|ДБР| = Н|Д£|со8(Д£^),
(9)
аналогичным по структуре равенству (3), где
(
Я ^ Лер
Н = 1 Н1^ , Я1 = Н1 . I Я1J Л 3
Подобный способ построения окружностей применяется и для последующих точек пути нагружения, для которых выполняется условие Як > Як-1, где к — номер точки на пути нагружения (к = 1, 2, ..., Л). Точку, после которой это условие нарушается (например, точка 3 на рис. 2), назовем точкой реверса. Строгая математическая формулировка условия появления точки реверса в пространстве девиаторов напряжений дана в работе [16]. После прохождения такой точки возникает новая поверхность, которая становится активной (подвижной), в то время как все остальные поверхности фиксируются. Центр активной окружности располагается на радиусе 33 (точка 4' на рис. 2). Направление приращения пластических деформаций Д№ совпадает с радиусом 4'4. При вычислении приращения пластической деформации можно исполь-
Рис. 2. Схематическое представление движения поверхностей нагружения (окружность и дуги) согласно многоповерхностной теории с одной активной поверхностью (ее центр расположен в точке 4'). Путь нагружения 0 ^ 1 ^ 2 ^ 3 ^ 4 в пространстве напряжений указан тонкими стрелками; в точках 2 и 4 жирными стрелками показаны направления приращений пластических деформаций; 1 и 3 — точки реверса
зовать соотношение (9) с заменой величин Я1, Н1 на Я3, Н3, соответствующих точке реверса 3.
Методика проведения опытов и образцы
При проведении экспериментов в условиях объемного напряженного состояния использовалась модернизированная установка [27], состоящая из камеры высокого давления для испытания образцов в условиях объемного напряженного состояния, из источников давления, гидравлических систем высокого и низкого давления, контрольно-измерительных приборов и различных приспособлений для обслуживания.
Для приготовления образцов использовались отрезки цельнотянутых труб длиной 4 м, из одной партии стали 1Х18Н10Т-ВД (ЭП502ВД). Размеры образцов: длина 170 мм;
номинальный наружный диаметр — 26,4 мм;
толщина стенки — 0,3 мм; отношение толщины стенки образца к его среднему диаметру — 0,011.
Указанное отношение соответствует рекомендуемым значениям и обеспечивает хорошую равномерность распределения кольцевых и радиальных напряжений по толщине стенки. Это позволяет вычислять напряжения для среднего слоя по формулам для безмоментных оболочек. Химический состав материала труб (сталь) представлен работе [18].
методика проведения опытов основана на том, что процесс нагружения рассматривается в пространстве напряжений Х1 и Х2, а процесс деформирования - в пространстве пластических деформаций Ер и Ер:
= стг - "2(сте + )> = )>
Еp = рp еp = — (р p + F p) M > 2 „J3 z eh
(10)
где ог, о0, ог - ненулевые компоненты тензора напряжений в цилиндрической системе координат, совпадающие с главными напряжениями; ер, ер, ер - главные пластические деформации.
Для вычисления длины пути пластического деформирования (параметра Одкви-ста ер ) нами использовалось соотношение
ер ={7 (й Ер )2 + (й Ер )2. (11)
Гидростатическое давление (среднее напряжение) определялось на основе зависимости
V3
(12)
а)
Ч»=0°
В основе методики нагружения лежало ступенчатое приложение нагрузок, позволявшее вести нагружение вдоль любого направления в пространстве Х1, Х2, которое характеризуется параметром вида девиато-ра у:
у = arctg (S2 / S1). (13)
Нагружение может быть реализовано в трубчатом образце (т. е. при у = 0, 1, ..., 120°). Напряженное состояние образца определяется параметром вида девиато-ра (13), интенсивностью напряжений о. и отношением инвариантов о / о..
r mean' 1
Эксперименты проводились в камере высокого давления и включали нагрузки и разгрузки образцов по прямым лучевым путям в пространстве напряжений, с варьированием параметра у. В качестве первого шага выполнялось осевое растяжение при у = 0° до некоторой выбранной интенсивности напряжений ога, затем следовала полная разгрузка, кольцевое растяжение при значениях у = 115 — 120° до приблизительно такой же интенсивности напряжений, и после этого — опять полная разгрузка до о. = 0 (см. лучи 0 - 1, 1 - 0, 0 - 4, 4 - 0 на рис. 3, b).
Переход при нагружении с луча у = 0° на луч у = 115 — 120° использовался потому, что давал возможность получать наибольшие пластические деформации, по сравнению с другими возможными реализациями лучевых путей повторного нагружения. На каждом прямолинейном участке нагруже-
Рис. 3. Пути нагружения образцов (указаны стрелками) для идентификации параметров моделей (а) и сравнение прогнозов рассматриваемых моделей (Ь); о.0 - интенсивность напряжений, соответствующая первоначальному нагружению
ния до различных значении о.0 проводились промежуточные разгрузки для измерения осевоИ (ег) и окружной (е0) деформации. Количество образцов составило 15. Уровни начальных интенсивностей напряжений о.0 изменялись от 680 до 900 МПа.
Диаграмма деформирования. Для определения ряда параметров структурной модели и многоповерхностной модели используется диаграмма деформирования (кривая зависимости интенсивности напряжений от интенсивности деформаций) при монотонном одноосном нагружении. Оценка возможности ее использования для описания сложного напряженного состояния (гипотеза единой кривой) требует проверки, которая также проводилась в рамках настоящего исследования.
Для аппроксимации диаграммы деформирования в данной работе использовалась степенная зависимость:
а,. = A(sf )т + B,
(14)
содержащая три параметра материала: A, B и т. Текущие размеры образца определяются с учетом зависимости (14).
Для получения диаграммы деформирования проводились испытания на растяжение тонкостенных трубчатых образцов на прессе Гагарина и в камере высокого давления.
В той же камере были проведены опыты по растяжению трубчатых образцов при повышенных значениях среднего напряжения о для оценки влияния отношения
mean ^ 1
о /о. на характер зависимости а, (sf). От-
mean i г г i v i '
ношение среднего напряжения к интенсивности напряжений менялось в пределах 0,330 - 0,466.
Путь нагружения в координатах Х1, Х2 представляет собой прямую линию: луч, исходящий из начала координат. Для получения зависимостей а,, (sf) и значений о, е. в момент разрушения, в процессе нагружения проводились промежуточные разгрузки для измерения деформаций. Значения параметра вида девиатора у (см. выше) меняли от 0 до 120° с шагом в 15°, отношение о /о. менялось в пределах от 0,321
mean ,
до 0,616. По каждому лучу было испытано 5 — 8 образцов.
определение параметров структурной модели. Для нахождения параметров структурной модели использовалась методика, предложенная в работе [13]. В соответствии с данной методикой, их определяли по аппроксимирующей ломаной линии, вписанной в экспериментальную диаграмму циклического деформирования при одноосном растяжении-сжатии.
Чтобы изучить поведение параметров модели при дальнейшем нагружении (после первичного растяжения), были проведены опыты как на повторно-лучевое, так и на повторно-циклическое виды нагруже-ния. При этом принималось во внимание изменение значений С, С.
Р' У
определение параметров многоповерхностной теории пластичности. Для нахождения параметров указанной теории пластичности с одной активной поверхностью нагружения нами рассматривалась степенная аппроксимация модуля пластической податливости Н в следующем виде:
[ Я V
H = сил
vai0 у
(15)
где о.0 — максимальная интенсивность напряжений, достигнутая на предыдущих этапах сложного пути нагружения; Н0 — модуль податливости максимально достигнутого уровня о.0 (Н0 = ( ер / (ст(); Я — радиус окружности равных податливостей, соответствующий текущему значению о.; С, п — параметры модели.
Податливость Н0 определялась по уравнению единой кривой (14):
Н0=тл Г- (16)
Значения Н и Я рассчитывались по результатам опытов, например, после экспериментов на чередующиеся осевые и кольцевые растяжения трубчатого образца.
Результаты экспериментального определения параметров моделей
В результате проведенных опытов на растяжение тонкостенных трубчатых образцов на прессе Гагарина и в камере высокого давления на десяти образцах была
л о.
ее с ш
СО
со CD
Й
1400 1200 1000 800 600 400 200
х о—: (к
• 0
"О"
f
0.1
0.2
0.3
Strain intensity
0.4
0.5
0.6
Рис. 4. Диаграмма деформирования стали 1Х18Н10Т-ВД
(данные от разных образцов представлены различными символами)
установлена следующая аппроксимация диаграммы деформирования (рис. 4):
ст, = 1470е°'678 + 325 (МПа). (17)
Получены следующие значения предела прочности (в истинных напряжениях): ои ~ 1200 — 1240 МПа и истинная деформация при разрыве е^« 0,5.
В опытах по растяжению трубчатых образцов при повышенных значениях среднего напряжения о было установлено,
r mean J 1
что в диапазоне исследованных о /о.
mean ,
(0,330 — 0,466) существенных отклонений как по характеру зависимости о (е), так и по значениям ои, еиг отмечено не было.
В результате проведенных опытов на 5 — 8 образцах были получены зависимости ст , (еР) и значения о,, е, в момент разрушения (как отмечалось выше, в процессе нагружения проводились промежуточные разгрузки для измерения деформаций).
На основании полученных экспериментальных данных был сделан вывод о возможности использования концепции единой кривой. Отклонения от среднего значения по напряжениям для различных путей нагружения не превышали ± 4 %. Установлено, что параметр у влияет как на значения о,, е, в момент разрушения, так и на характер разрушения.
Параметры структурной модели. В результате проведенных экспериментальных
исследовании при сложном непропорциональном нагружении были найдены параметры структурной модели при различных значениях о00 (табл. 1).
В работе [13] значения параметров модели для материала принимались посто-
Таблица 1
Значения параметров структурной модели при различных значениях с 0
Параметр Значение, МПа
о ,0 = 772 МПа ,0 о„ = 871 МПа ,0
С а 325 325
Св 711 765
С Y 770 820
E az 94 100 163 000
Eae 126 800 178 800
E а 110 400 170 900
Eez 2 960 8 560
Еве 3 610 9 230
E 3 280 8 890
E YZ 1 500 2 110
EYe 2 370 1 860
E Y 1 930 1 980
Значения параметра е0 составили 21,04 и 21,49 % при о ,0 = 772 и 871 МПа соответственно.
янными, не зависящими от истории его нагружения. Анализ экспериментальных результатов, полученных нами при повторно-лучевых нагружениях, показал, что численное значение этих параметров зависит от уровня интенсивности напряжений о0, достигаемого в базовом опыте. Поэтому опыты проводились при различных значениях о(0 (о(0 = 700 - 900 МПа). В табл. 1 даны в качестве примера соответствующие значения для двух значений оя.
Расчеты параметров Ср, С для опытов с различными значениями о0 позволили заключить о приблизительно линейной зависимости указанных параметров от вели-
чины о.п:
СТ;0(5) = 0,886 - 0,9; Ср(5)
ст
,0(33)
С
= 0,929 - 0,9;
Р(33)
С.
у(5)
С.
= 0,939 - 0,9.
у(зз)
В результате анализа поведения параметров модели при дальнейшем нагружении (после первичного растяжения), при котором, как уже отмечалось, были проведены опыты как на повторно-лучевое, так и на повторно-циклическое виды нагружения, была установлена закономерность изменения параметров, позволяющая более точно прогнозировать результаты эксперимента [29]. При этом принималось во внимание изменение значений С , С .
Р' У
Параметры многоповерхностной теории пластичности. Податливость Н0 определялась на основе уравнения (16):
Н =
1,475 (ст,.0 - 325
1470 I 1470
(МПа). (18)
Значения параметров Н и Я были рассчитаны по результатам опытов на чередующиеся осевые и кольцевые растяжения трубчатого образца.
Необходимо отметить, что в базовом опыте характер зависимости о.(е.) заметно меняется при переходе от первого кольцевого растяжения после первого осевого на-гружения к последующим чередующимся осевым и кольцевым растяжениям. Поэтому, как и в работе [13], сделано предположение об изменении свойств материала
после первых циклов сложного нагружения и принято, что первому нагружению при максимальном отличии по углу у отвечают значения параметров С = 1,6; п = 25,3, а последующим нагружениям, независимо от пути нагружения, отвечают значения С = 0,2; п = 5,9.
Учитывая характер аппроксимации зависимостей, мы пришли к выводу, что величина интенсивности напряжений о.0 при первичном растяжении в базовом опыте не влияет на значения коэффициентов А и п.
Кроме вышеизложенных, были опробованы и другие аппроксимирующие зависимости для пластических податливостей, в частности
Н = Н Ю(0'01892- 1б(Я/ст, 0 ))-04-4
Н = Н 10[(-®( Я/4 0)-0,03)/0,085-0,6]/0,2
Результаты расчетов с применением этих соотношений в целом дали худшее соответствие данным опытов, по сравнению с результатами расчетов, использующих соотношения (15).
Результаты испытаний и сопоставительный анализ моделей
Повторно-пропорциональные пути нагружения. Многоповерхностная теория пластичности с одной активной поверхностью достаточно хорошо описывает процесс на-гружения по пропорциональным путям с разгрузками и переходами на другое направление пропорционального нагружения. Структурная модель элемента материала с постоянными параметрами дает более заметные расхождения результатов с данными опытов.
Параметры структурной модели, определяемые в базовом опыте, зависят от величины интенсивности напряжений ст,. 0 при первичном нагружении образца, а уровень значений ст, при последующих циклах возрастал от цикла к циклу. Поэтому вычисления проведены с использованием эмпирической зависимости, установленной в опытах [29]:
С" = С'. I 1 + 0,35
ст
- Я
Я
(19)
где С— значение параметра С, отвечающее текущему значению наибольшей интенсивности напряжений о.. тах, индекс ' = Р, У (см. рис. 1, Ь); С' — значение параметра С], соответствующее значению предыдущей наибольшей интенсивности напряжений; Я — значение предыдущей наибольшей интенсивности напряжений (радиус гиперповерхности, соответствующей о. тах на предыдущих этапах нагруже-ния (см. рис. 3, Ь).
В табл. 2 приведены в качестве примера типичные данные, полученные нами в опытах по деформированию образцов в указанных условиях.
Структурная модель элемента материала с постоянными коэффициентами обеспечивает качественное соответствие экспериментальным данным, но количественное сопоставление приводит к выводу о большом отклонении от них. Применение коррекции параметров приводит к существенному улучшению прогноза, использующего соотношения структурной модели, и приближает по эффективности скорректированную модель к многоповерхностной теории.
Таким образом, результаты расчета деформаций по соотношениям многоповерхностной теории с одной активной поверхностью и расчета по структурной модели, использующей коррекцию параметров, дают хорошее соответствие с данными экспериментов.
Сложные пути нагружения. На рис. 5 показаны экспериментальные результаты, найденные при сложных путях нагружения для одного из таких характерных процессов (см. рис. 3, Ь); приведенные данные сопоставлены с расчетными, полученными для различных моделей материала. Параметр Одквиста на момент разрушения (на второй стадии) достиг следующих уровней: опытное значение — 0,340; прогнозируемое при использовании структурной модели с коррекцией параметров — 0,251;
прогнозируемое согласно многоповерхностной теории — 0,345.
таким образом, прогностические расчеты по обеим моделям дали удовлетворительное совпадение с данными эксперимента.
При малых деформациях несколько лучшее соответствие достигалось с использованием структурной модели материала, при развитых деформациях — модели, основанной на концепции поля пластических податливо-стей. Это объясняется тем, что первая из этих сопоставляемых моделей более точно аппроксимирует начальный участок диаграммы деформирования, но хуже отражает участок развитых деформаций. Этот недостаток можно устранить, если корректировать параметры модели с учетом повышения уровня достигаемых интенсивностей напряжений. такая коррекция существенно приближает ее прогноз к данным эксперимента (см. рис. 5).
Результаты расчетов по теории пластического течения с кинематическим упрочнением (вторая из двух рассматриваемых моделей) указывают на неадекватный (неудовлетворительный) прогностический характер ее описания, относящегося к процессу сложного циклического деформирования.
Структурная модель элемента материала в целом несколько лучше описывает процесс сложного повторно-пропорционального циклического нагружения, чем многоповерхностная теория пластичности с одной активной поверхностью. Это можно объяснить тем, что модельное описание пути этих нагружений повторяет путь нагруже-ния в базовом опыте.
однако следует отметить, что для повышения адекватности структурной модели потребовалась коррекция ее параметров. Такая мера указывает на ее большую чувствительность к условиям (по напряжениям) проведения базового эксперимента и к уровню интенсивностей напряжений, достигаемых в ходе опытов.
Многоповерхностная теория пластичности с одной активной поверхностью лучше описывает исследованные сложные циклические пути нагружения, которые могут существенно отличаться от пути нагружения в базовом опыте. Это говорит о ее универсальном характере в отношении расчетов деформирования образцов материала и элементов конструкций, подверженных сложным циклическим нагружениям.
Таблица 2
Значения параметров пластического деформирования при повторно-лучевых путях нагружения после первичного осевого растяжения
Напряжение в конце участка, МПа Наклон пути нагружения Ч>, ' AEf, % AEf, % Asf,% sf,%
Опыт Расчет Опыт Расчет Опыт Расчет Опыт Расчет
МТ С CK МТ С CK МТ С CK МТ С CK
-332 890 736 0 115 0 -0,9 2,0 -0,6 2,1 -5,6 7,5 -1,0 2,0 1,6 -0,3 1,1 -0,2 8,4 "1,2 -0,1 -0,0 1,8 2,0 1,2 2,2 0,4 8,1 2,1 2,5 1,8 3,8 1,2 3,4 10.4 18.5 2,1 4,6
-330 921 732 0 115 0 "1,1 3,6 -1,0 2,3 -3,2 6,1 -0,5 2,2 1,0 -0,1 1,3 "0,1 3,8 "1,3 0,6 "1,4 1,4 3,6 1,6 2,4 5,3 6,8 1,0 2,8 5,3 8,9 5,1 7,5 23,8 30,6 5,5 8,3
-328 957 722 0 115 0 -0,6 3,2 -0,9 2,6 "2,9 6,1 -0,4 2,3 0,8 0,0 1,1 "0,1 3,0 -1,5 0,2 -0,8 1,0 3,2 1,4 2,7 4,5 6,9 0,7 3,0 9,8 13,0 8,9 11,6 35,1 42,0 9,0 12,0
-332 991 717 0 115 0 -0,4 3,9 -0,9 2,5 "3,7 6,8 -0,4 2,3 0,5 0,2 1,0 -0,0 3,5 "1,9 0,1 "1,7 0,7 3,9 1,3 2,6 5,4 7,7 0,7 2,9 13,7 17,6 12,9 15,5 47,3 55,0 16,7 15,6
-338 1006 717 0 115 0 -0,6 2,1 -0,8 1,3 -5,0 6,7 -0,4 1,3 0,4 -0,0 0,9 -0,0 4,3 "2,1 0,1 -0,8 0,7 2,1 1,2 1,3 6,8 7,7 0,9 2,0 18,3 30,3 16,8 18,1 61,9 69,6 16,5 18,5
-343 1020 722 0 115 0 -0,5 1,9 -0,9 1,2 -5,5 7,2 -0,6 1,4 0,4 0,1 0,9 -0,0 4,3 "2,5 "0,1 -0,6 0,6 2,0 1,2 1,3 7,3 8,3 1,1 2,1 21,0 22,9 19,3 20,7 76,8 85,1 19.5 21.6
-344 1070 713 0 115 0 -0,4 5,2 -0,8 3,8 -5,5 9,7 -0,4 3,3 0,3 0,7 0,8 -0,0 3,9 "3,4 0,0 -2,0 0,5 5,2 1,2 3,9 7,0 1,0 0,7 4,1 23,4 28,6 21,8 25,7 92,1 103 22,4 26,4
-354 1110 718 0 115 0 -0,6 4,2 -0,8 3,1 "8,9 2,5 -0,8 3,1 0,2 1,6 0,7 -0,0 5,9 -4,6 -0,2 "1,1 0,6 4,5 1,1 3,2 1,0 4,1 1,4 3,9 29,3 33,8 26,8 30,0 114 128 27,8 31,7
Примечание. Первичное осевое растяжение проведено до ол = 872 МПа и = 0,245.
Принятые сокращения: МТ — многоповерхностная теория пластичности с одной активной поверхностью, С, СК — соответственно структурная модель с постоянными коэффициентами и с коррекцией параметров материала по формуле (19).
4 0.35 е?
1 3 4 \J.Zo 2 ч ч \
1 1 V * \ * 1 ) > у >
/ / г ч ^^^^ > и. 14 •«ч ц^ > ^Ч^ Ч ч Ч ► ч, /
5 и. и*™ Ч¥ о* -г
-0.08 -0.07 -0.06 -0.05 -0.04 -0.03 -0.02 -0.01
о
0.01
Рис. 5. Сравнение экспериментальных (1) и расчетных (2 — 5) результатов путей деформирования с использованием различных моделей для пути нагружения (см. рис. 3, Ь): многоповерхностной теории с одной активной поверхностью (2); структурных моделей с коррекцией (3) и без коррекции (4) параметров; теории пластического течения с кинематическим упрочнением (5)
Для многоповерхностной теории не выявлено существенной зависимости ее параметров от уровня нагружения как в базовом опыте, так и при произвольных путях нагружения. Она адекватно описывает сложное циклическое нагружение с интенсив-ностями напряжений как выше, так и ниже интенсивности напряжений, при которой определялись параметры модели. Дополнительной коррекции параметров при этом не требуется.
Заключение
Проведено сравнение прогнозов многоповерхностной теории с одной активной поверхностью и структурной модели с результатами экспериментальных исследований на стали 1Х18Н10Т-ВД при сложном непропорциональном пассивном нагруже-нии.
Многоповерхностная теория пластичности с одной активной поверхностью на-гружения обеспечивает хорошую точность прогнозов для сложных путей нагружения, с участками полной и частичной разгрузки.
Для всех исследованных путей нагру-жения результаты расчетов с использованием структурной модели при неизменных значениях параметров модели и числе элементов дает худшее соответствие данным опытов, чем результат расчетов по многоповерхностной модели. Прогнозы с применением теории пластического течения с кинематическим упрочнением значительно уступают по точности двум вышеназванным моделям.
Вариант структурной модели с переменными параметрами позволяет значительно уменьшить погрешности прогноза. Погрешности расчетов с использованием такой модели и многоповерхностной теории близки для всех исследуемых путей нагружения.
В рамках многомодельного метода [30] исследования пластического деформирования установлены области достоверного применения многоповерхностной теории с одной активной поверхностью нагружения и структурной модели материала.
Работа выполнена при финансовой поддержке гранта РФФИ № 16-08-00845-а.
СПИСОК ЛИТЕРАТУРЫ
[1] Зубчанинов В.Г. Об активных и пассивных процессах сложного нагружения-разгружения // Современные проблемы термовязкопластич-ности. Труды II школы-семинара, М.: МАМИ, 2007. С. 3-18.
[2] Зубчанинов В.Г., Гультяев В.И., Зубчанинов Д.В. Экспериментальное исследование процессов сложного нагружения материалов на многозвенных траекториях // Современные проблемы термовязкопластичности. труды II школы-семинара. М.: МАМИ, 2007. С. 19-24.
[3] Зубчанинов В.Г., Гультяев В.И., Зубча-нинов Д.В. Экспериментальное исследование стали 45 по траекториям типа «Плоский винт» // Современные проблемы прочности, пластичности и устойчивости. Тверь: ТГТУ, 2007. С. 104-109.
[4] Кадашевич Ю.И., Мосолов А.Б. Эндо-хронные теории пластичности: основные положения, перспективы развития // Известия АН СССР. Механика твердого тела. 1989. № 1. С. 161-168.
[5] Кадашевич Ю.И., Мосолов А.Б. О соотношениях эндохронной теории пластичности с «новой» мерой внутреннего времени // Технология легких сплавов. 1990. № 3. С. 32-36.
[6] Пальмов В.А. Колебания упругопласти-ческих тел. М.: Наука. 1976. 328 с.
[7] Ильюшин А.А. Пластичность. Основы общей математической теории. М.: Изд-во АН СССР, 1963. 271 с.
[8] Lissenden C.J., Colaiuta J.F., Lerch B.A. Hardening behavior of three metallic alloys under combined stresses at elevated temperature // Acta Mechanica. 2004. Vol. 169. No. 1. Pp. 53-77.
[9] Kuwabara T., Kuroda M., Tvergaard V., Nomura K. Use of abrupt strain path change for determining subsequent yield surface: experimental study with metal sheets // Acta Mater. 2000. Vol. 48. No. 9. Pp. 2071-2079.
[10] Khan A.S., Parikh Y. Large deformation in polycrystalline copper under combined tensiontorsion, loading, unloading and reloading or reverse loading: a study of two incremental theories of plasticity // International Journal of Plasticty. 1986. Vol. 2. No. 4. Pp. 379-392.
[11] Rees D.W.A. Research note on a reexamination of neutral loading experiments // International Journal of Plasticity. 1988. Vol. 4. No. 1. Pp. 91 -102.
[12] Rees D.W.A. Applications of classical plasticity theory to non-radial loading paths // Proc. R. Soc. Lond. A. 1987. Vol. 410. No. 1839. Pp. 443-475.
[13] Павлов П.А. Основы инженерных расчетов элементов машин на усталость и длительную прочность. Ленинград: Машиностроение, 1988, 251 с.
[14] Мельников Б.Е., Изотов И.Н., Кузнецов
Н.П. Пространственный упругопластический элемент для описания нестационарного деформирования // Известия ВНИИГ им. Б.Е. Веденеева. 1986. Т. 197. № 1. С. 75-79.
[15] Мельников Б.Е., Изотов И.Н., Кузнецов Н.П. Расчет и экспериментальное исследование сложных путей упругопластического деформирования // Проблемы прочности. 1990. № 8. С. 14-17.
[16] Мельников Б.Е., Семенов А.С. Многоповерхностная теория пластичности с одной активной поверхностью пластических податливо-стей // Труды ЛГТУ. 1991. № 144. С. 26-31.
[17] Мельников Б.Е., Семенов А.С. Описание эволюции поверхностей равной податливости в многоповерхностной теории пластичности // Труды СПбГТУ. 1996. № 456. С. 52-61.
[18] Изотов И.Н., Кузнецов Н.П., Мельников Б.Е., Митюков А.Г. Упруго-пластическое поведение стальных трубчатых образцов при сложном трехосном нагружении // Научно-технические ведомости СПбГТУ. 2003. № 3 (33). С. 114-125.
[19] Изотов И.Н., Мельников Б.Е., Семенов А.С., Семенов С.Г. Определение свойств пластических податливостей на никелевых образцах // Научно-технические ведомости СПбгПУ. 2009. № 2. С. 165-174.
[20] Изотов И.Н., ягн Ю.И. Изучение пластического деформирования металла с деформационной анизотропией, созданной в процессе предварительного нагружения // ДАН СССР. 1961. Т. 139. № 3. С. 576-579.
[21] Mroz Z. On the description of anisotropic workhardening // J. Mech. Phys. Solids. 1967. Vol. 15. No. 3. Pp. 163-175.
[22] Mroz Z., ^mes V.A., Zienkiewicz o.C. An anisotropic hardening model for soils and its application to cyclic teding // Int. J. for Numerical and Analytical Methods in Geomechanics. 1978. Vol. 2. No. 3. Pp. 203-221.
[23] Мельников Б.Е., Кузнецов Н.П., Митю-ков А.Г., Семенов А.С., Семенов С.Г. Применение модели энергетической оценки прочности при нескольких циклах нагружения в случае сложного неоднородного напряженного состояния // Современные проблемы ресурса материалов и конструкций. Труды III школы-семинара 2009 г. М.: МГТУ (МАМИ), 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 31th Summer School "Advanced Problems in Mechanics". June 22-July 2, 2003. St. Petersburg: Polytechnical Publishing House, 2003. Pp. 172-176.
[26] Артюх В.Г. Точность предохранителей для металлургических машин. Мариуполь: Изд-во ПГТУ, 2000. 177 с.
[27] Павлов П.А., Изотов И.Н., Кузнецов Н.П., Митюков А.Г., Симанёнок А.Н. Установка
для испытания материалов при объемном напряженном состоянии // Проблемы прочности. 1986. № 12. С. 104-107.
[28] Иванченко С.Н., Кадашевич Ю.И., По-мыткин С.П. Квазистатическая теория пластичности, учитывающая уплотнение материала // Известия РАН. Механика твердого тела. 2000. № 5. С. 21-26.
[29] Izotov I.N., Kuznetsov N.P., Melnikov B.E., Mityukov A.G., Musienko A.Yu. New variants of the multisurface theory of plasticity. Comparison with the experimental data // Proc. of SPIE. 2000. Vol. 4046. Pp. 362-367.
[30] 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. Supрl. 2. Pp. 615-616.
СВЕДЕНИЯ ОБ АВТОРАХ
КУЗНЕцоВ Николай Петрович — кандидат технических наук, доцент кафедры технологии машиностроения Новгородского государственного университета им. Ярослава Мудрого, г. Великий Новгород, Российская Федерация.
173003, Российская Федерация, г. Великий Новгород, Большая Санкт-Петербургская ул., 41.
МЕЛьНИКоВ Борис Евгеньевич — доктор технических наук, заведующий кафедрой сопротивления материалов Санкт-Петербургского политехнического университета Петра Великого, г. Санкт-Петербург, Российская Федерация.
195251, Российская Федерация, г. Санкт-Петербург, Политехническая ул., 29
СЕМЕНоВ Артем Семенович — кандидат физико-математических наук, доцент кафедры механики и процессов управления Санкт-Петербургского политехнического университета Петра Великого, г. Санкт-Петербург, Российская Федерация.
195251, Российская Федерация, г. Санкт-Петербург, Политехническая ул., 29
Kuznetsov N.P., Melnikov B.E., Semenov A.S. VERIFICATION AND ADAPTATION OF THE PLASTICITY MODELS UNDER COMPLEX VARIABLE LOADING WITH INTERMEDIATE COMPLETE AND PARTIAL UNLOADINGS.
The experimental studies of elasto-plastic deformation of tubular steel samples under proportional and non-proportional (monotonic and cyclic) loadings, including partial and intermediate loadings, have been conducted with the aim of improving the accuracy of the description of the complex passive loading processes. The plastic strain accumulation was observed in the course of tests carried out under passive loading. However, this effect turned out not to be described by the plastic flow theory. This result required the development of an alternative material model. The comparisons of experimental results with the predictions of the structural (rheological) elasto-plastic model and the multisurface theory of plasticity with one active surface were made. Modifications of the constitutive equations were proposed in order to improve the accuracy of the material response prediction.
PLASTICITY, PASSIVE LOADING, EXPERIMENT, MODELING.
REFERENCES
[1] V.G. Zubchaninov, Ob aktivnykh i passivnykh protsessakh slozhnogo nagruzheniya-razgruzheniya [On the active and passive processes of complex loading-unloading], Moscow, Sovremennyye problemy termovyazkoplastichnosti, Trudy II shkoly-seminara. (2007) 3—18.
[2] V.G. Zubchaninov, V.I. Gultyayev, D.V. Zubchaninov, Eksperimentalnoye issledovaniye protsessov slozhnogo nagruzheniya materialov na mnogozvennykh trayektoriyakh [An experimental investigation of processes of complex loading of materials in the multilink paths], Moscow, sovremennyye problemy termovyazkoplastichnosti, Trudy II shkoly-seminara. (2007) 19—24.
[3] V.G. Zubchaninov, V.I. Gultyayev, D.V. Zubchaninov, Eksperimentalnoye issledovaniye stali 45 po trayektoriyam tipa "Ploskiy vint" [An experimental investigation of steel 45 in the "flat screw" paths] Tver, Sovremennyye problemy prochnosti, plastichnosti i ustoychivosti. (2007) 104-109.
[4] Yu.I. Kadashevich, A.B. Mosolov,
Endokhronnye teorii plastichnosti: osnovnyye polozheniya, perspektivy razvitiya [Endochronic plasticity theories: the main statements, future development], Izvestiya Rossiyskoy akademii nauk. Mekhanika tverdogo tela. (1) (1989). 161-168.
[5] Yu.I. Kadashevich, A.B. Mosolov, O sootnosheniyakh endokhronnoy teorii plastichnosti s "novoy" meroy vnutrennego vremeni [On the relationships of endochronic plasticity theory with a "new" measure of internal time], Tekhnologiya legkih splavov. (3) (1990) 32 -36.
[6] V.A. Palmov, Kolebaniyauprugoplasticheskikh tel [Vibrations of elastic-plastic bodies], Moscow, Nauka, 1976.
[7] A.A. Ilyushin, Plastichnost. Osnovy obshchey matematicheskoy teorii [Plasticity. Background of the general mathematical theory], Moscow, Izd-vo AN sssR, 1963.
[8] C.J. Lissenden, J.F. Colaiuta, B.A. Lerch, Hardening behavior of three metallic alloys under combined stresses at elevated temperature, Acta Mechanica 169 (1) (2004) 53-77.
[9] T. Kuwabara, M. Kuroda, V. Tvergaard, K. Nomura, Use of abrupt strain path change for determining subsequent yield surface: experimental study with metal sheets, Acta Mater. 48 (9) (2000) 2071-2079.
[10] A.S. Khan, Y. Parikh, Large deformation in polycrystalline copper under combined tensiontorsion, loading, unloading and reloading or reverse loading: a study of two incremental theories of plasticity, International Journal of Plasticity. 2 (4)
(1986) 379-392.
[11] D.W.A. Rees, Research note on a reexamination of neutral loading experiments, International Journal of Plasticity. 4 (1) (1988) 91-102.
[12] D.W.A. Rees, Applications of classical plasticity theory to non-radial loading paths, Proc. R. Soc. Lond. A. 410 (1839) (1987) 443-475.
[13] P.A. Pavlov, Osnovy inzhenernykh raschetov elementov mashin na ustalost i dlitelnuyu prochnost [Background of fatigue and long-time strength engineering calculations of machine elements], Leningrad, Izd-vo "Mashinostroyeniye", 1988.
[14] B.E. Melnikov, I.N. Izotov, N.P. Kuznetsov, Prostranstvennyy uprugoplasticheskiy element dlya opisaniya nestatsionarnogo deformirovaniya [A spatial elastic-plastic element for description of unsteady deformation], Leningrad, Izvestiya VNIIG im. B.E. Vedeneyeva. 197 (1986) 75-79.
[15] B.E. Melnikov, I.N. Izotov, N.P. Kuznetsov, Raschet i eksperimentalnoye issledovaniye slozhnykh putey uprugoplasticheskogo deformirovaniya [Calculations and experimental investigation of complex paths of elastic-plastic deformation], Problemy prochnosti. (8) (1990) 14-17.
[16] B.E. Melnikov, A.S. Semenov, Mnogopoverkhnostnaya teoriya plastichnosti s odnoy aktivnoy poverkhnostyu plasticheskikh podatlivostey [The multisurface theory of plasticity with one active surface of plastic compliances], Trudy SPbGTU. (144) (1991) 26-31.
[17] B.E. Melnikov, A.S. Semenov, Opisaniye evolyutsii poverkhnostey ravnoy podatlivosti v mnogopoverkhnostnoy teorii plastichnosti [Description of evolution of equal-compliance surfaces in the multisurface theory of plasticity], Trudy SPbGTU. (456) (1996) 52-61.
[18] I.N. Izotov, N.P. Kuznetsov, B.E. Melnikov, A.G. Mityukov, Uprugoplasticheskoye povedeniye stalnykh trubchatykh obraztsov pri slozhnom trekhosnom nagruzhenii [Elastic-plastic behavior of steel tubular samples under the complex triaxial loading], St. Petersburg State Polytechnical University Journal, (3) (2003) 114-125.
[19] I.N. Izotov, B.E. Melnikov, A.S. Semenov, et al., Opredeleniye svoystv plasticheskikh podatlivostey na nikelevykh obraztsakh [Determination of plastic compliance properties of nickel samples], St. Petersburg State Polytechnical University Journal, (2) (2009) 165-174.
[20] I.N. Izotov, Yu.I. Yagn, Izucheniye plasticheskogo deformirovaniya metalla s deformatsionnoy anizotropiyey, sozdannoy v protsesse predvaritelnogo nagruzheniya [Study of
plastic deformation of metal with strain anisotropy induced during preliminary loading], DAN SSSR. 139 (3) (1961) 576-579.
[21] Z. Mroz, On the description of anisotropic workhardening, J. Mech. Phys. Solids. 15 (3) (1967) 163-175.
[22] Z. Mroz, V.A. Norries, o.C. Zienkiewicz,
An anisotropic hardening model for soils and its application to cyclic loading, Int. J. for Numerical and Analytical Methods in Geomechanics. 2 (3) (1978) 203-221.
[23] B.E. Melnikov, N.P. Kuznetsov, A.G. Mityukov, et al., Primeneniye modeli energeticheskoy otsenki prochnosti pri neskolkikh tsiklakh nagruzheniya v sluchaye slozhnogo neodnorodnogo napryazhennogo sostoyaniya [An application of the model for an energetic estimation of low-cycle fatigue in the case of complex inhomogeneous stress state], Sovremennyye problemy resursa materialov i konstruktsiy. Trudy III shkoly-seminara 2009 g. M. MGTU (MAMI) (2009) 140-147.
[24] B.E. Melnikov, A.S. Semenov, Energeticheskaya model nakopleniya povrezhdeniy pri malotsiklovoy ustalosti [Energetic model of damage accumulation in the low-cycle fatigue], XXI Peterburgskiye chteniya po problemam prochnosti, K 100-letiyu so dnya rozhdeniya L.M. Kachanova i Yu.N. Rabotnova Sb. Mater. (2014) 59-61.
[25] D.A. Kitaeva, Y.I. Rudaev, About kinetic
equations of dynamic superplasticity model, Proceedings of the 31th Summer School "Advanced Problems in Mechanics", June 22—July 2, 2003. St. Petersburg, Polytechnical Publishing House, (2003) 172—176.
[26] V.G. Artyukh, Tochnost predokhraniteley dlya metallurgicheskikh mashin [Safety-device accuracy for metallurgical machines], Mariupol, Izd-vo PGTU, 2000.
[27] P.A. Pavlov, I.N. Izotov, N.P. Kuznetsov, Ustanovka dlya ispytaniya materialov pri obyemnom napryazhennom sostoyanii [The testing machine in the volumetric stress], Problemy prochnosti, (12) (1986) 104—107.
[28] S.N. Ivanchenko, Yu.I. Kadashevich, S.P. Pomytkin, Kvazistaticheskaya teoriya plastichnosti, uchityvayushchaya uplotneniye materiala [Quasi-static plasticity theory considering the material consolidation], Izvestiya Rossiyskoy Akademii nauk. Mekhanika tverdogo tela. (5) (2000) 2 —26.
[29] I.N. Izotov, N.P. Kuznetsov, B.E. Melnikov, et al., New variants of the multisurface theory of plasticity, Comparison with the experimental data, Proc. of SPIE. 4046 (2000) 362—367.
[30] B.E. Melnikov, A.S. Semenov, Creation and application of hierarchical sequence of material models for numerical analysis of elasto-plastic structures, Zeitschrift fur Angewandte Mathematik und Mechanik. 76 (Supl. 2) (1996) 615—616.
THE AUTHORS
kuznetsov Nikolay P.
Yaroslav-the-Wise Novgorod State University
41 Bolshaya Sankt-Peterburgskaya St., 173003, Veliky Novgorod, Russian Federation [email protected]
melnikov Boris E.
Peter the Great St. Petersburg Polytechnic University
29 Politekhnicheskaya St., St. Petersburg, 195251, Russian Federation
semenov Artem S.
Peter the Great St. Petersburg Polytechnic University
29 Politekhnicheskaya St., St. Petersburg, 195251, Russian Federation
© Санкт-Петербургский политехнический университет Петра Великого, 2017
