CC BY
11
УДК 62-462:539.41 DOI 10.18698/0536-1044-2017-10-10-17
Разработка и обоснование инженерной методики расчета напряженного состояния винтовых трубок, нагруженных внутренним давлением
А.Г. Сорокина
МГТУ им. Н.Э. Баумана, 105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1
The Development and Validation of Engineering Methods of Calculating the Stress State of Spiral Tubes Loaded by Internal Pressure
A.G. Sorokina
BMSTU, 105005, Moscow, Russian Federation, 2nd Baumanskaya St., Bldg. 5, Block 1 e-mail: [email protected]
Разработана упрощенная методика расчета напряжений в винтовых трубках, нагруженных внутренним давлением. Показано, что в большинстве случаев координатные линии на срединной поверхности трубки можно приближенно рассматривать как линии кривизны. Безмоментные уравнения равновесия в частных производных вследствие геликоидальной симметрии удалось свести к одному обыкновенному дифференциальному уравнению для мембранного усилия, действующего вдоль винтовых линий. Оказалось, что решение дифференциального уравнения практически совпадает с результатом применения котельных формул, в которых используется нормальный радиус кривизны винтовой линии. Инженерная методика продемонстрирована на винтовой трубке с заданным соотношением размеров.
Ключевые слова: винтовая трубка, геликоидальная симметрия, безмоментная теория, уравнения равновесия, котельные формулы, метод наименьших квадратов
A simplified method of calculating stresses in spiral tubes under internal pressure was developed. It was shown that in most cases the coordinate lines on the middle surface of the tube could be approximated as lines of curvature. The membrane equilibrium equations in partial derivatives due to the helical symmetry were reduced to one ordinary differential equation for the membrane force acting along the helix. It transpired that the solution to the differential equation almost coincided with the result of the application of boiler formulas that used the normal radius of curvature of the helix. The engineering method was demonstrated on a spiral tube with a given ratio of sizes.
Keywords: spiral tube, helical symmetry, membrane theory, equilibrium equations, boiler formulas, least square method
Винтовые трубки (ВТ) применяют в энергетическом машиностроении, химической промышленности, приборостроении, медицинском оборудовании и других отраслях техники. На рис. 1 показана ВТ c образующим контуром в виде окружности.
Такие трубки без ограничения общности можно отнести к классу геликоидально симметричных оболочек (ГСО), срединная поверхность которых образована винтовым движением замкнутого контура. Основным видом нагрузки для ВТ является внутреннее давление
Рис. 1. ВТ с образующим контуром в виде окружности
или разность давлений снаружи и внутри ВТ. Очевидно, что для трубок с внутренним давлением напряженное состояние является близким к безмоментному. Однако в отличие от цилиндрических оболочек, нагруженных внутренним давлением, для которых вычисление напряжений в материале трубки сводится к применению двух простейших соотношений (котельные формулы), расчет напряженного состояния ГСО даже в безмоментном случае — нетривиальная задача.
Основная трудность расчета ГСО заключается в том, что винтовые линии не являются ортогональными к линиям замкнутого контура, вследствие чего необходимо использовать косоугольные гауссовы координаты. Применение косоугольных координат приводит к довольно громоздким системам уравнений, которые можно решить только приближенными или численными методами [1-14]. Обзор различных подходов к расчету оболочек циклического строения, в том числе и ГСО, выполнен в работе [15].
Альтернативой косоугольным координатам служит переход к так называемому нормальному контуру, т. е. такому контуру, который ортогонален винтовым линиям. Предыдущие исследования автора показали, что нормальный контур весьма удобен для ГСО открытого профиля, таких как ленточные пружины. Использование нормального контура в случае малых перемещений приводит к уравнениям теории оболочек в ортогональных гауссовых координатах, не совпадающих с линиями кривизны. Такие уравнения применены, например, в работе [16]. Для витых и винтовых трубок нормальный контур, как правило, является разомкнутым, что порождает неудобства, например, при по-
становке граничных условий. Таким образом, как использование косоугольных координат, так и переход к нормальному контуру, не приводят к удобным методикам расчета замкнутых ГСО.
Цель работы — разработка и обоснование инженерной методики расчета напряжений, возникающих под действием внутреннего давления в ВТ. Предлагается модифицировать котельные формулы таким образом, чтобы в измененном виде они стали применимы для расчета напряжений в ВТ.
Геометрия срединной поверхности ГСО и ВТ.
Срединная поверхность ГСО (рис. 2) задается векторно-матричным соотношением
г(«, ф) = ^ -^^(ф^а), (1)
2к
где г (а, ф) — радиус-вектор срединной поверхности ГСО; а — гауссова координата контура; k — орт оси г; Т — период винтовой поверхности; L — матрица поворота вокруг оси г на угол ф; в(а) — радиус-вектор контура.
Выражение (1) означает, что геликоидальная поверхность образована поворотом контура вокруг оси г на угол ф с одновременным переносом контура вдоль той же оси, т. е. его винтовым движением. Для перехода к проекциям достаточно конкретизировать в выражении (1) матрицу поворота и радиус-вектор в(а):
^ (ф)] =
С08 ф 8Ш ф 0
- 8Ш ф С08 ф 0
в(а) =
ах (а) ау (а) аг (а)
Рис. 2. Геликоидальная поверхность и ее контур
Я
а\
—ГХ
Еч
^ I ^"
О х
Рис. 3. Параметры ВТ
Для ВТ (рис. 3) контур является окружностью. Ее плоскость перпендикулярна винтовой линии, поэтому контур можно представить в параметрическом виде как уравнение окружности в плоскости хОг, которая дополнительно повернута вокруг оси х на угол у:
"10 0 К Я-0 0
a(a) =
0
С08 у 81п у
0
- 81п у С08 у
ь 81п а 0
' со8 а
Л
у = агйд
Т
2пЯ'
где у и Я — угол подъема и радиус винтовой линии, проходящей через центры окружностей; Ь — радиус окружности; а — полярный угол в плоскости хОг.
Согласно работе [17], параметры Ламе А, В и единичные векторы касательные к ко-
ординатным линиям, имеют вид
А = ,
В =
Эг ^ t ЭГ
Эа) ' 1 АЭа'
ЭГ V t ЭГ
Эф) ' 2 ВЭф
Отметим, что параметры Ламе не зависят от гауссовой координаты ф в силу геликоидальной симметрии. Параметр А совпадает с радиусом окружности (А = Ь), для параметра В получается довольно громоздкое выражение, поэтому оно здесь не представлено.
Синус угла X между координатными линиями и нормаль п к поверхности определяются по формулам
81п %= ^ Xt2
П :
^ X t2 81п X
Так как выражения получаются довольно громоздкими, вместо формул далее будут представлены графики для конкретного соотношения размеров
Я =-
Т = 4Ь.
(2)
На рис. 4 показана ВТ с параметрами (2).
Визуально линии контура на рис. 4 практически перпендикулярны винтовым линиям. Это подтверждается графиком синуса угла между координатными линиями (рис. 5).
На представленном графике отличие 81пх от единицы не превышает 10 %. Для ВТ с меньшим, чем в выражении (2) отношением Ь/Я отличие 81пх от единицы будет еще незначительнее. В связи с этим можно сделать первое упрощающее предположение: координатную
сетку, образованную окружностями контура и винтовыми линиями, приближенно можно рассматривать как ортогональную.
Согласно работе [17], кривизны координатных линий 1/Rl, 1/R2 и кручение поверхности 1АК12 можно найти по формулам
1 Э2г
П;
R[ A2дa2 _ д^ R2 _ В2дФ2 1 д2r
п;
(3)
R12 АВдадф
п,
где индекс « » означает, что данная кривизна не обязательно является главной (у главной кривизны штрих не указывается [17]).
Результаты расчета по формулам (3) представлены в виде графиков (рис. 6).
Как видно из полученных графиков, кручение поверхности 1Л^2 заметно меньше кривизны линии контура 1/RÍ, поэтому в выражениях, содержащих обе эти величины, кручением поверхности можно пренебречь (1/R12 = 0). Из этого следует второе упрощающее предположение: координатные линии приближенно можно рассматривать как линии кривизны.
Уравнения равновесия безмоментной ВТ и их решение. В трубках, нагруженных внутренним давлением, в основном реализуется безмомент-
вшх 1,0 р=
0,8
1 2 3 4 5 6 а
Рис. 5. График синуса угла между координатными линиями
—, —, —
6 а
Рис. 6. Графики кривизн b/RÍ (1), b/R2 (2) и кручения b/R12 (3) винтовой поверхности
ное напряженное состояние, поэтому справедливы обычные уравнения безмоментной теории оболочек [17]:
А(ТВ) + -^(А2)-|ВГ2 + ABql _0; да Адф ' да
А) + -^(В2)-дАТ1 + ABq2 _ 0; (4) дф Вда ' дф
Т1 х Т2 — ^ —_ qз, Rl ^
где Т1, Т2 и Б — усилия соответственно окружное, направленное вдоль винтовых линий и сдвигающее; qь q2, qз — распределенные нагрузки.
Учтем, что в силу геликоидальной симметрии все величины не меняются по винтовой координате ф, а также то, что оболочка нагружена только давлением и кривизны координатных линий можно рассматривать как главные:
ql _ q2 _ 0; qз _ р; _ _ _
Rl Rl Ь R2 R2
где р — внутреннее давление в ВТ.
Тогда система (4) существенно упрощается:
йа^ 1 ) йа 2
д
Вда
(В2)_ 0;
Т Т2 — + — _ р. Ь R2
Из второго уравнения следует Б = С/В2, где С — постоянная интегрирования, но так как крутящего момента в сечении нет, то Б = 0. Оставшиеся соотношения несложно сводятся к одному дифференциальному уравнению для Т2
Ь | р - А 1 В1 - йВт2 _ 0.
йа ^ ^ R2) ) й а
(5)
После определения Т2 усилие Т1 можно найти по формуле
Т1 _ь ( р - I
Некоторым осложнением расчета является то, что кривизна 1^2 = 1^2 обращается в нуль во внутренних точках интервала 0 < а < 2л (см. рис. 6). В связи с этим численное интегри-
рование уравнения (5) невозможно. Вместо численного интегрирования применен метод наименьших квадратов. Усилие T2 представлено отрезком тригонометрического ряда
T (а) = ¿ T2(c) cos ka -+ ¿ T2k sin ka, (6)
k=0 k=1
где n — максимальный номер гармоники; индекс «(с)» относится к слагаемым, симметричным относительно начала отсчета угла а; индекс «(s)» — к кососимметричным слагаемым.
Подстановка разложения (6) в уравнение (5) дает функцию невязки дифференциального уравнения равновесия
Оценка сходимости численной процедуры
f (a)=da
p - Ti 1 в) - dB.T¡.
R2) ) da
(7)
Невязка уравнения равновесия оболочки в сечении, совпадающем с контуром, определяется выражением
2л
g = рлЬ2 - J T2Ada. (8)
0
Общая невязка всей задачи, составленная из невязок (7) и (8), имеет вид
J =
f „ 12 2?Гf (a)12
pb2
da,
(9)
0 рЬ2
где множители в знаменателе добавлены для обезразмеривания невязок.
Интегралы в выражениях (8) и (9) вычислялись как конечные суммы. Минимизация невязки (9) выполнялась стандартной процедурой FindMinimum пакета Wolfram Mathematica [18]. В результате получены коэффициенты разложения (6) и построены графики мембранных усилий, представленные на рис. 7.
Из графиков видно, что напряжения в ВТ на 40 % больше, чем в прямой трубке такого же
Рис. 7. Распределение мембранных усилий T1/(pb) (1) и T2/(pb) (2) вдоль контура ВТ
Количество гармоник n Количество точек интегрирования • f T1 1 min I — I 1 pb) f T11 max I — I 1 pb)
20 12 0,838096 0,999821
20 24 0,705377 0,858443
20 48 0,822827 1,42188
20 96 0,822126 1,42408
20 120 0,821949 1,42463
20 162 0,821748 1,42526
4 162 0,821703 1,42501
диаметра, и опасная точка находится на наименьшем расстоянии от оси 2.
Для оценки сходимости вычислительного процесса рассмотрено различное количество гармоник в разложении (6) и точек интегрирования на контуре (см. таблицу).
Для построения графиков усилий Т1, Т2 были выбраны параметры, при превышении которых три значащие цифры в результатах не меняются: п = 20, количество точек интегрирования, равное 162. Как видно из таблицы, количество гармоник и точек интегрирования можно брать существенно меньшим, но вычисления происходят довольно быстро, поэтому данные параметры назначались «с запасом».
На рис. 7 усилие Т2 практически не отличается от постоянной, равной рЬ/2. В связи с этим можно принять третье упрощающее предположение: усилие ^ приближенно можно рассматривать как постоянное и определять его из уравнения равновесия торца по котельной формуле, т. е. Х2 ~pb/2. Тогда расчет упрощается полностью и сводится к инженерной методике определения напряжений в ВТ:
а! = TL , pbf 1 -_b_ h h f 2R2
а2 = * « ^, 2 h 2h
(10)
где а1 — окружное напряжение; Н — толщина ВТ; а2 — напряжение, действующее вдоль винтовых линий.
Таким образом, предлагаемую инженерную методику можно рассматривать как формальное распространение котельных формул на ВТ. Однако следует учитывать, что соотношения (10) обоснованы гораздо более серьезными
соображениями, представленными выше, чем традиционные, принимаемые при выводе котельных формул.
Для вычисления нормального радиуса кривизны винтовой линии в соотношении (10) достаточно использовать теорему Менье о связи кривизны нормального и наклонного сечений [17], из которой следует
1 _ sin а R2 R + b sin а
(
cos
arctg
T
\
2n(R + bsin а)
(11)
где c помощью арктангенса можно определить угол подъема винтовой линии для текущего радиуса поверхности, равного R + b sin а. Если в формуле (11) выполнить тригонометрические преобразования, то она несколько упрощается:
1 4л2 (R + b sin a)sin а
(12)
R2 T2 + 4л2 (R + bsin а)2
Мембранные усилия, найденные по инженерной методике, т. е. по формулам (10)—(12), сопоставлялись с более точными данными, рассчитанными выше путем минимизации невязки (9). Результат сопоставления показан на рис. 8, где T* — усилие, определенное по инженерным формулам (10)-(12).
Как видно из полученного графика, погрешность инженерной методики для ВТ с соотношением размеров R : b = 5:3, T: b = 4 составляет менее 3 %.
Известно, что безмоментное решение несправедливо вблизи линий, на которых радиус кривизны R1 или R2 обращается в бесконечность, поэтому в узких зонах, примыкающих к винтовым линиям а = 0 и а = л, возникают дополнительные моментные напряжения Однако указанные зоны невелики, и при постоянном во времени давлении и пластичном материале локальные всплески напряжений прак-
6 а
Рис. 8. Оценка погрешности приближенных формул (10)-(12)
тически не сказываются на несущей способности оболочек [17].
В связи с тем, что введенные выше гипотезы вносят некоторые погрешности в расчет напряжений, при использовании предложенной инженерной методики рекомендуется назначать коэффициент запаса на 10... 15 % большим, чем при обычном применении без-моментной теории оболочек.
Выводы
1. Предложена инженерная методика приближенного расчета напряжений в ВТ, нагруженных внутренним давлением, которая сводится к использованию выражений (10), формально совпадающих с котельными формулами.
2. Методика обоснована подробным анализом геометрических соотношений и уравнений равновесия, применяемых в безмоментной теории оболочек в криволинейных гауссовых координатах.
3. Представленный пример демонстрирует надежность методики даже для ВТ с малым отношением радиуса винтовой линии к радиусу окружности контура.
Литература
[1] Krivoshapko S.N. Geometry and strength of general helicoidal shells. Applied Mechanics Re-
views, 1999, vol. 52, no. 5, pp. 161-175.
[2] Krivoshapko S.N. Static analysis of shells with developable middle surfaces. Applied Mechan-
ics Reviews, 1998, vol. 51, no. 12, pt. 1, pp. 731-746.
[3] Халаби С.М. Моментная теория тонких винтовых псевдо-торсовых оболочек. Строи-
тельная механика инженерных конструкций и сооружений, 2001, № 10, c. 61-67.
[4] Рынковская М.И. К вопросу о расчете на прочность тонких линейчатых винтовых обо-
лочек. Строительная механика инженерных конструкций и сооружений, 2015, № 6, c. 13-15.
[5] Рынковская М.И. К вопросу расчета прямых геликоидальных оболочек по методу
В.Г. Рекача. Строительная механика инженерных конструкций и сооружений, 2006, № 2, c. 63-66.
[6] Рынковская М.И. Применение и расчет геликоидальных оболочек в архитектуре и
строительстве. Вестник Российского университета дружбы народов. Серия: Инженерные исследования, 2012, № 4, c. 84-90.
[7] Туликова Е.М. Расчет тонких упругих оболочек в форме длинного косого геликоида.
Строительная механика инженерных конструкций и сооружений, 2015, № 3, c. 23-27.
[8] Тупикова Е.М. Вычисление перемещений и усилий в железобетонных пандусах в форме
косого геликоида. Строительная механика и расчет сооружений, 2015, № 4, c. 24-28.
[9] Тупикова Е.М. Полуаналитический расчет оболочки в форме длинного пологого косого
геликоида в неортогональной несопряженной системе координат по моментной теории. Строительная механика инженерных конструкций и сооружений, 2016, № 3, c. 3-8.
[10] Тупикова Е.М. Анализ метода В.Г. Рекача для расчета напряженно-деформированного состояния оболочки в форме длинного пологого косого геликоида. Строительная механика и расчет сооружений, 2016, № 1, с. 14-20.
[11] Savicevic S., Janjic М., Vukcevic M., Sibalic N. Stress research of helicoidal shell shape elements. Machines, technologies, materials, 2013, is. 10. URL: http://www.mech-ing.com/ journal/Archive/2013/10/42_Savicevic_mtm13.pdf (дата обращения 20 мая 2017).
[12] Savicevic S. A Development of Automatized Projection of Construction Elements of Helical Shell Shape. PhD dissertation, Podgorica, Faculty of Mechanical Engineering, 2001.
[13] Krivoshapko S.N., Gbaguidi A.G. Two methods of analysis of thin elastic open helicoidal shells. International Journal of Research and Reviews in Applied Sciences, 2012, vol. 12, no. 3, pp. 382-390.
[14] Knabel J., Lewinski T. Selected equilibrium problem of thin elastic helicoidal shells. Archives of Civil Engineering, 1999, vol. 45(2), pp. 245-257.
[15] Krivoshapko S.N., Christian A. Bock Hyeng. Static and dynamic analysis of thin-walled cyclic shells. International Journal of Modern Engineering Research, 2012, vol. 2, is. 5, pp. 3502-3508.
[16] Шевелев Л.П., Корихин Н.В., Головин А.И. Состояния поля напряжений в геликоидальной оболочке. Строительство уникальных зданий и сооружений, 2014, № 2(17), c. 25-38.
[17] Бидерман В.Л. Механика тонкостенных конструкций: статика. Москва, URSS, 2017. 496 c.
[18] Дьяконов В.П. Mathematica 5.1/5.2/6. Программирование и математические вычисления. Москва, ДМК-Пресс, 2008. 574 с.
References
[1] Krivoshapko S.N. Geometry and strength of general helicoidal shells. Applied Mechanics Re-
views, 1999, vol. 52, no. 5, pp. 161-175.
[2] Krivoshapko S.N. Static analysis of shells with developable middle surfaces. Applied Mechan-
ics Reviews, 1998, vol. 51, no. 12, pt. 1, pp. 731-746.
[3] Khalabi S.M. Momentnaia teoriia tonkikh vintovykh psevdo-torsovykh obolochek [Moment
theory of thin helical pseudo-tarcowie shells]. Stroitel'naia mekhanika inzhenernykh kon-struktsii i sooruzhenii [Structural Mechanics of Engineering Constructions and Buildings]. 2001, no. 10, pp. 61-67.
[4] Rynkovskaia M.I. K voprosu o raschete na prochnost' tonkikh lineichatykh vintovykh ob-
olochek [On problem of strength analysis of thin linear helicoidal shells]. Stroitel'naia mek-hanika inzhenernykh konstruktsii i sooruzhenii [Structural Mechanics of Engineering Constructions and Buildings]. 2015, no. 6, pp. 13-15.
[5] Rynkovskaia M.I. K voprosu rascheta priamykh gelikoidal'nykh obolochek po metodu
V.G. Rekacha [To the question of the calculation of the direct helicoidal shells by the method of V.G. Rekach]. Stroitel'naia mekhanika inzhenernykh konstruktsii i sooruzhenii [Structural Mechanics of Engineering Constructions and Buildings]. 2006, no. 2, pp. 63-66.
[6] Rynkovskaia M.I. Primenenie i raschet gelikoidal'nykh obolochek v arkhitekture i stroitel'stve
[On application and analysis of helicoidal shells in architecture and civil engineering]. Vestnik Rossiiskogo universiteta druzhby narodov. Seriia: Inzhenernye issledovaniia [RUDN journal of engineering researches]. 2012, no. 4, pp. 84-90.
[7] Tupikova E.M. Raschet tonkikh uprugikh obolochek v forme dlinnogo kosogo gelikoida
[Analysis of the thin elastic shells in the form of long oblique helicoid]. Stroitel'naia mek-hanika inzhenernykh konstruktsii i sooruzhenii [Structural Mechanics of Engineering Constructions and Buildings]. 2015, no. 3, pp. 23-27.
[8] Tupikova E.M. Vychislenie peremeshchenii i usilii v zhelezobetonnykh pandusakh v forme
kosogo gelikoida [Calculation of displacements and forces in reinforced concrete ramps in the form of skew helicoid]. Stroitel'naia mekhanika i raschet sooruzhenii [Structural Mechanics and Analysis of Constructions]. 2015, no. 4, pp. 24-28.
[9] Tupikova E.M. Poluanaliticheskii raschet obolochki v forme dlinnogo pologogo kosogo ge-
likoida v neortogonal'noi nesopriazhennoi sisteme koordinat po momentnoi teorii [Emi-analytical analysis of a long shallow oblique helicoidal shell in a non-orthogonal non-conjugate coordinate system]. Stroitel'naia mekhanika inzhenernykh konstruktsii i sooruzhenii [Structural Mechanics of Engineering Constructions and Buildings]. 2016, no. 3, pp. 3-8.
[10] Tupikova E.M. Analiz metoda V.G. Rekacha dlia rascheta napriazhenno-deformirovannogo sostoianiia obolochki v forme dlinnogo pologogo kosogo gelikoida [Rekatch's method of stress-strain analysis of the shell of long shallow oblique helicoid form]. Stroitel'naia mekhanika i raschet sooruzhenii [Structural Mechanics and Analysis of Constructions]. 2016, no. 1, pp. 14-20
[11] Savicevic S., Janjic M., Vukcevic M., Sibalic N. Stress research of helicoidal shell shape elements. Machines, technologies, materials, 2013, is. 10. Available at: http://www.mech-ing.com/journal/Archive/2013/10/42_Savicevic_mtm13.pdf (accessed 20 May 2017).
[12] Savicevic S. A Development of Automatized Projection of Construction Elements of Helical Shell Shape. PhD dissertation, Podgorica, Faculty of Mechanical Engineering, 2001.
[13] Krivoshapko S.N., Gbaguidi A.G. Two methods of analysis of thin elastic open helicoidal shells. International Journal of Research and Reviews in Applied Sciences, 2012, vol. 12, no. 3, pp. 382-390.
[14] Knabel J., Lewinski T. Selected equilibrium problem of thin elastic helicoidal shells. Archives of Civil Engineering, 1999, vol. 45(2), pp. 245-257.
[15] Krivoshapko S.N., Christian A. Bock Hyeng. Static and dynamic analysis of thin-walled cyclic shells. International Journal of Modern Engineering Research, 2012, vol. 2, is. 5, pp. 3502-3508.
[16] Shevelev L.P., Korikhin N.V., Golovin A.I. Sostoianiia polia napriazhenii v gelikoidal'noi obolochke [Field stresses state in the helicoidal shell]. Stroitel'stvo unikal'nykh zdanii i sooruzhenii [Construction of Unique Buildings and Structures]. 2014, no. 2(17), pp. 25-38.
[17] Biderman V.L. Mekhanika tonkostennykh konstruktsii: statika [Mechanics of thin-walled structures: statics]. Moscow, URSS publ., 2017. 496 p.
[18] D'iakonov V.P. Mathematica 5.1/5.2/6. Programmirovanie i matematicheskie vychisleniia
[Mathematica 5.1/5.2/6. Programming and mathematical calculations]. Moscow, DMK-Press publ., 2008. 574 p.
Информация об авторе
СОРОКИНА Анастасия Германовна (Москва) — старший преподаватель кафедры «Прикладная механика». МГТУ им. Н.Э. Баумана (105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1, e-mail: [email protected]).
Статья поступила в редакцию 15.06.2017 Information about the author
SOROKINA Anastasia Germanovna (Moscow) — Senior Lecturer, Department of Applied Mechanics. Bauman Moscow State Technical University (105005, Moscow, Russian Federation, 2nd Baumanskaya St., Bldg. 5, Block 1, e-mail: [email protected]).
