Научная статья на тему 'Методика оценки прочности металлических конструкций подъемно-транспортных машин'

Методика оценки прочности металлических конструкций подъемно-транспортных машин Текст научной статьи по специальности «Строительство и архитектура»

CC BY
243
42
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
СТАЛЬНАЯ КОНСТРУКЦИЯ / STEEL STRUCTURE / СВАРНОЕ СОЕДИНЕНИЕ / ИСПЫТАНИЕ / TESTING / ПЛАСТИЧНОСТЬ / PLASTICITY / ПРОЧНОСТЬ / STRENGTH / МЕТОД КОНЕЧНЫХ ЭЛЕМЕНТОВ. / FINITE ELEMENT METHOD. / WELDING JOINT

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Плотников Дмитрий Георгиевич, Соколов Сергей Алексеевич, Боровков Алексей Иванович, Михайлов Александр Александрович

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

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

Похожие темы научных работ по строительству и архитектуре , автор научной работы — Плотников Дмитрий Георгиевич, Соколов Сергей Алексеевич, Боровков Алексей Иванович, Михайлов Александр Александрович

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

Strength estimation method for metal structures of hoisting equipment

The main problem considered in this article is prediction of strength for welding structures of hoisting equipment. Developed is a method of estimation of plastic strain limits for welding structures that takes into account all the most important factors influencing strength and plastic properties of a welding joint (such as stress concentration factors and residual stresses, environment temperature). These factors are implemented with the help of a special welded model. The static test, described in this article, was conducted at normal temperature. Analysis of the stress-strain state of the fracture area is based on the finite element method. As a result of investigation the limits of the plasticity of the welding joint were estimated.

Текст научной работы на тему «Методика оценки прочности металлических конструкций подъемно-транспортных машин»

DOI 10.8562/JEST.214.21 УДК 620.171.2

Д.Г. Плотников, С.А. Соколов, А.И. Боровков, А.А. Михайлов

МЕТОДИКА ОЦЕНКИ ПРОЧНОСТИ МЕТАЛЛИЧЕСКИХ КОНСТРУКЦИЙ ПОДЪЕМНО-ТРАНСПОРТНЫХ МАШИН

D.G. Plotnikov, S.A. Sokolov, A.I. Borovkov, A.A. Michailov

STRENGTH ESTIMATION METHOD FOR METAL STRUCTURES

OF HOISTING EQUIPMENT

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

СТАЛЬНАЯ КОНСТРУКЦИЯ; СВАРНОЕ СОЕДИНЕНИЕ; ИСПЫТАНИЕ; ПЛАСТИЧНОСТЬ; ПРОЧНОСТЬ; МЕТОД КОНЕЧНЫХ ЭЛЕМЕНТОВ.

The main problem considered in this article is prediction of strength for welding structures of hoisting equipment. Developed is a method of estimation of plastic strain limits for welding structures that takes into account all the most important factors influencing strength and plastic properties of a welding joint (such as stress concentration factors and residual stresses, environment temperature). These factors are implemented with the help of a special welded model. The static test, described in this article, was conducted at normal temperature. Analysis of the stress-strain state of the fracture area is based on the finite element method. As a result of investigation the limits of the plasticity of the welding joint were estimated.

STEEL STRUCTURE, WELDING JOINT, TESTING, PLASTICITY, STRENGTH, FINITE ELEMENT METHOD.

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

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

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

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

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

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

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

Методика экспериментального исследования

Прочность сварного узла в условиях эксплуатации зависит от совокупности факторов:

механических свойств материала в зоне сварного шва;

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

температуры эксплуатации;

скорости деформирования материала.

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

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

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

а)

бГ

15

1 10t

30t

0,7t

А-А

<j 1 И

б)

Рис. 1. Образец для испытаний: а — чертеж образца; б — схема распределения зон в образце: 1 — основной металл пластин; 2 — наплавленный металл; 3 — зона термического влияния; 4 — точка максимальной концентрации напряжений

4

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

Для экспериментального исследования использован образец, все элементы которого были вырезаны из листовой стали марки СтЗсп толщиной 20 мм. Испытания проводились на разрывной машине Instron KN с максимальным усилием 600 кН. Продольная нагрузка передавал ась на образец через хвостовые пластины, которые зажимались в захватах разрывной машины. При этом фиксировалась диаграмма растяжения в координатах «перемещение захватов — нагрузка». Испытание производилось при нормальной температуре +23 °С, скорость перемещения захватов составляла 20 мм/с. В результате нагружения образца до усилия 580 кН он деформировался с изгибом в центральном сечении при суммарном перемещении захватов 5 = 8 мм. Образец не разрушился, нагружение прекратилось из-за исчерпания возможности

разрывной машины. Раскрытие стыка составило 6 мм, в основании концентратора образовались надрывы материала глубиной 4 мм.

Для определения механических свойств стали были проведены испытания стандартных цилиндрических образцов основного металла и зоны термического влияния по ГОСТ 1497—84. Испытания показали, что основной металл имеет предел текучести стт = 310 МПа, временное сопротивление ств = 500 МПа, относительное удлинение 55 = 0,3. Металл зоны термического влияния имеет предел текучести стт = 403 МПа, временное сопротивление ств = 575 МПа, относительное удлинение 55 = 0,21.

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

Рис. 2. Конечно-элементная модель: а — общий вид; б — конечно-элементная сетка в зоне вершины концентратора

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

Конечно-элементная модель строилась с учетом симметрии образца, поэтому соответствовала только его четверти (рис. 2). Свойства симметрии были реализованы с помощью граничных условий (на плоскости перемещение их = 0, на плоскости ХУ перемещение и = 0). В связи с тем, что поле напряжений и деформаций в зоне концентратора имеет весьма значительные градиенты, достоверность результатов расчета существенно зависит от качества конечно-элементной сетки. Поэтому разбивка была выполнена таким образом, чтобы в вершине трещиноподобного концентратора обеспечива-

лась регулярная сетка и размер элемента не больше 0,02а (где а — характерный размер трещины, в данном случае равный толщине листа 20 мм). Для построения модели использовались линейные элементы типа SOLID (CHEXA) с 8 узлами. Материалы основных элементов и зоны термического влияния принимались изотропными с билинейными диаграммами механических свойств [9], которые строились на основании данных, полученных при испытании цилиндрических образцов. Для обоих материалов модуль упругости принимался равным Е = 2,1-105 МПа, модуль упрочнения Ет = 1000 МПа.

Поле остаточных сварочных напряжений в образце создавалось путем моделирования термомеханического цикла сварки. Он включал нагрев материала сварного шва от температуры окружающей среды 20 °С до 1500 °С и последующее его остывание до температуры 20 °С. Коэффициент теплового расширения линейно возрастал от 11,5-10—6 при температуре 20 °С до 18,9-10-6 при температуре 1500 °С. Граничные условия определялись свойствами симметрии

образца. В результате этого расчета максимальные значения остаточных сварочных напряжений в зоне концентрации достигали 239 МПа, что составляет 80 % от предела текучести основного металла.

Моделирование захватов в процессе на-гружения модели осуществлялось с помощью граничных условий, обеспечивавших поступательное перемещение захватной части образца (иу = 0, и = 0 на поверхностях, контактирующих с захватами). При этом захваты и элементы разрывной машины, связывающие их с датчиком перемещений, предполагались абсолютно жесткими. Нагружение модели было кинематическим, то есть задавалось монотонное перемещение, и вычислялась требуемая сила. Расчет производился по стандартному алгоритму расширенного нелинейного статического анализа SOL600 в Ж Nastran.

Анализ результатов

В результате численного анализа получено поле напряжений и деформаций в модели, построен график «перемещение захватов — нагрузка», который сопоставлен с экспериментальным (рис. 3). Как видно из графиков, общая картина процесса деформирования образца удовлетворительно описывается численным расчетом. Имеется расхождение графиков в области малых деформаций. Это можно объяснить влиянием упругих деформаций захватов на начальном участке нагружения.

Локальными характеристиками напряжен -но-деформированного состояния в концентраторе являются интенсивность пластической деформации ер1 и коэффициент жесткости напряженного состояния п . Интенсивность пластических деформаций вычисляется как

V =-

42

epl - ep2 ) +(ep2 -

ep 2" ep3

)2+(■

ep3 - epl

Для анализа сходимости решения и оценки предельной пластической деформации в сечении А-А (рис. l, в) был проведен анализ трех вариантов моделей с минимальными размерами элементов 0,25 мм, 0,5 мм и l мм. По результатам этих расчетов построена зависимость максимальных пластических деформаций, соответствующих перемещению захватов на 8 мм, от размера элемента (рис. 3, б). Деформации усреднялись в объеме куба со стороной l мм. Расчетное значение пластических деформаций определялось с помощью экстраполяции на размер элемента, равного нулю. Таким образом, было получено значение предельной пластической деформации eplc = 0,44. Размер пластической зоны при максимальной нагрузке составил гпл = = 13,8 мм. Максимальное значение коэффициента жесткости напряженного состояния в процессе деформирования с увеличением нагрузки изменяется от 2,5 до 1,7 .

Таким образом, методика решения вопроса о возможности использования конкретной стали с определенной сварочной технологией для изготовления сварных конструкций включает следующие мероприятия. Необходимо изготовить серию образцов указанной конфигурации из проката наибольшей толщины с применением тех же сварочных процедур. Провести их испытания при минимальной эксплуатационной температуре по ГОСТ 15150. По результатам испытаний рассчитать значение предельных пластических деформаций в образце. Полученные данные позволят применить для сварных конструкций деформационный критерий прочности в форме [6]

V

^ ePicln ,

где ер1, ер2, ер3 — значения главных пластических деформаций в элементе.

Коэффициент жесткости напряженного состояния определяется как п =—^ [6], где а! —

первое главное напряжение; аут — эквивалентное напряжение, вычисленное согласно теории прочности Мизеса — Генки.

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

а)

Усилие, кН

500 400 300 200 100 0

*

/

1

7

б)

Пластические деформации

6 8 Перемещение захватов, мм

0,8 Размер элемента, мм

Рис. З.Результаты расчета: а — диаграмма нагружения образца, экспериментальная и расчетная; б — зависимость расчетной максимальной пластической деформации от размера элемента; в — деформированная конечно-элементная модель (--эксперимент; —--расчет)

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

Результаты проведенного эксперимента сопоставлены с оценкой предельных деформаций по методике Н.А. Махутова. Согласно [5] предельное значение интенсивности пластической деформации может быть представлена в виде ек = ек- Бе(ц), где ек- — предельная пластическая деформация, определенная при одноосном растяжении образца цилиндрической формы; Бе — коэффициент снижения предельной пластической деформации; п — коэффициент жесткости напряженного состояния.

На основании работы [5] построена приближенная зависимость коэффициента снижения предельной пластической деформации от коэффициента жесткости напряженного со-

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

1 2

А (п)=-V (1)

п2

Предельное значение интенсивности пластической деформации для испытанного образца можно оценить по методике [5], используя значение коэффициента жесткости, полученное из численного анализа, и значение предельной пластической деформации ек-, определенное по результатам испытания цилиндрических образцов. В результате вычислений получаем ек = 0,38. Таким образом, аналитическая оценка предель-

ной пластической деформации в образце весьма близка к максимальному значению интенсивности пластической деформации, полученному с помощью конечно-элементного расчета — ек = 0,44, то есть отличается на 14 %. Расхождение можно объяснить тем, что разрывы, вероятно, появились в концентраторе не при максимальном перемещении захватов, а несколько раньше, что невозможно зафиксировать при большой скорости нагружения.

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

Выводы и рекомендации

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

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

СПИСОК ЛИТЕРАТУРЫ

1. РД 24.090.52—90. Подъемно-транспортные машины. Материалы для сварных металлических конструкций / ВНИИПТМАШ. М., 1990.

2. РД 22—16—2005. Машины грузоподъемные. Выбор материала для изготовления, ремонта и реконструкции сварных стальных конструкций / ФГУП СКТБ БК. 2005.

3. Federation Europeenne de la manutention. Section I: FEM 1.001. Rules for the design of hoisting appliances. 3rd Edition. 10.1998

4. ISO 20332. Cranes — Proof of competence of steel structures. First edition 2008—12—15.

5. Махутов Н.А. Деформационные критерии разрушения и расчет элементов конструкций на прочность. М.: Машиностроение, 1981. 272 с.

6. Копельман Л.А. Основы теории прочности сварных конструкций: Учебное пособие. 2-е изд., испр.

СПб.: Издательство «Лань», 2010. 457 с.

7. Соколов С.А. Металлические конструкции подъемно-транспортных машин. СПб.: Политехника, 2005. 423 с.

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

8. Одесский П.Д., Ведяков И.И., Горпинченко В.М.

Предотвращение хрупких разрушений металлических строительных конструкций. М: Изд-во СП Интермет Инжиниринг, 1998. 220 с.

9. Ullner С., Brüser S., Subaric-Leitis A., Reth-meier M. Determination of local stress-strain properties of resistance spot-welded of advanced high-strength steels using the instrumented indentation test./ - Journal of Materials Science. Vol. 47. Issue 3. 2012. P. 15041513.

10. ГОСТ 32578—2013. Краны грузоподъемные. Металлические конструкции. Требования к материалам, 2013.

REFERENCES

1. RD 24.090.52—90. Podyemno-transportnyye mash-iny. Materialy dlya svarnykh metallicheskikh konstruktsiy/ [RD 24.090.52-90. Hoisting equipment. Materials for welding metal structures]./ VNIIPTMASh. M., 1990. (rus.)

2. RD 22—16—2005. Mashiny gruzopodyemnyye. Vy-bor materiala dlya izgotovleniya, remonta i rekonstruktsii svarnykh stalnykh konstruktsiy/[RD 22-16-2005. Hoisting equipment. Choice of materials for repair and reconstruc-

tion of welding metal structures] / FGUP SKTB BK, 2005. (rus.)

3. Federation Europeenne de la manutention. Section I: FEM 1.001. Rules for the design of hoisting appliances. 3rd Edition. 10.1998. (rus.)

4. ISO 20332. Cranes - Proof of competence of steel structures. First edition 2008—12—15.

5. Makhutov N.A. Deformatsionnyye kriterii razrush-eniya i raschet elementov konstruktsiy na prochnost. [De-

formation fracture criterion and structure analysis of metal structures]. M.: Mashinostroyeniye, 1981. 272 s. (rus.)

6. Kopelman L.A. Osnovy teorii prochnosti svarnykh konstruktsiy: Uchebnoye posobiye. 2-ye izd., ispr. [Foundations of the theory of strength for welding structures].: SPb. Izd-vo «Lan», 2010. 457 s. (rus.)

7. Sokolov S.A. Metallicheskiye konstruktsii podyem-no-transportnykh mashin. [Metal structures of hoisting equipment]. SPb.: Politekhnika, 2005. 423 s. (rus.)

8. Odesskiy P.D. , Vedyakov I.I., Gorpinchenko V.M. Predotvrashcheniye khrupkikh razrusheniy metalli-

cheskikh stroitelnykh konstruktsiy / [Prevention of brittle fracture of metal structures]. M.: Izd-vo SP Inter-met Inzhiniring, 1998. 220 s. (rus.)

9. Ullner S. , Bruser S, Subaric-Leitis A., Reth-meier M. Determination of local stress-strain properties of resistance spot-welded of advanced high-strength steels using the instrumented indentation test. Journal of Materials Science. Vol. 47. Issue 3. 2012. P. 1504-1513 (rus.)

10. GOST 32578—2013. Krany gruzopodyemnyye. Metallicheskiye konstruktsii. Trebovaniya k materialam/ [GOST 32578-2013 Cranes. Metal structures. Requirements for materials]. 2013. (rus.)

СВЕДЕНИЯ ОБ АВТОРАХ

ПЛОТНИКОВ Дмитрий Георгиевич — аспирант Санкт-Петербургского политехнического университета Петра Великого. 195251, Россия, г. Санкт-Петербург, Политехническая ул., 29. E-mail: [email protected]

СОКОЛОВ Сергей Алексеевич — доктор технических наук профессор Санкт-Петербургского политехнического университета Петра Великого 195251, Россия, г. Санкт-Петербург, Политехническая ул., 29. E-mail: [email protected]

БОРОВКОВ Алексей Иванович — кандидат технических наук профессор проректор по перспективным проектам Санкт-Петербургского политехнического университета Петра Великого. 195251, Россия, г. Санкт-Петербург, Политехническая ул., 29. E-mail: [email protected]

МИХАЙЛОВ Александр Александрович — ассистент Санкт-Петербургского политехнического университета Петра Великого. 195251, Россия, г. Санкт-Петербург, Политехническая ул., 29. E-mail: [email protected]

AUTHORS

PLOTNIKOV Dmitry G. — Peter the Great St. Petersburg Polytechnic University. 29 Politechnicheskaya St., St. Petersburg, 195251, Russia. E-mail: [email protected]

SOKOLOV Sergey A. — Peter the Great St. Petersburg Polytechnic University. 29 Politechnicheskaya St., St. Petersburg, 195251, Russia. E-mail: [email protected]

BOROVKOV Alexey I. — Peter the Great St. Petersburg Polytechnic University. 29 Politechnicheskaya St., St. Petersburg, 195251, Russia. E-mail: [email protected]

MICHAILOV Alexandr A. — Peter the Great St. Petersburg Polytechnic University. 29 Politechnicheskaya St., St. Petersburg, 195251, Russia. E-mail: [email protected]

© Санкт-Петербургский политехнический университет Петра Великого, 2015

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