CC BY
10
#9(726) 2020
ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
3
Машиностроение и машиноведение
УДК 531.8:519.2 doi: 10.18698/0536-1044-2020-9-3-8
Численный метод расчета нелинейной упругой характеристики продольно-поперечного преобразователя
П.А. Пья, Ю.В. Григорьев
МГТУ им. Н.Э. Баумана
A Numerical Method for Calculating the Nonlinear Elastic Characteristic of Longitudinal-Transverse Transducers
Pyae Phyo Aung, Y.V. Grigoriev
Bauman Moscow State Technical University
Рассмотрен численный метод, позволяющий исследовать напряженно-деформированное состояние продольно-поперечного преобразователя и получать его нелинейную упругую характеристику. Предложена математическая модель, в которой использовано численное решение краевой задачи на базе уравнений плоского криволинейного стержня в программном пакете MATLAB. Напряженно-деформированное состояние и нелинейная упругая характеристика системы получены методом последовательных нагружений, основанных на линеаризованных уравнениях криволинейного стержня. Разработанную модель можно рассматривать как начальное приближение к решению пространственной задачи продольно-крутильного преобразователя.
Ключевые слова: метод последовательных нагружений, нелинейная задача, линеаризованная система уравнений, продольно-поперечный преобразователь, криволинейный стержень, нелинейная упругая характеристика
This paper presents a numerical method for studying the stress-strain state of longitudinal-transverse transducers and obtaining their nonlinear elastic characteristic. The authors propose a mathematical model that uses a direct numerical solution of the boundary value problem based on the plain curved rod equations in MATLAB. The stress-strain state and the nonlinear elastic characteristic of the system are obtained using a method of successive loading based on linearized equations of the curved rod. The proposed model can be considered as an initial approximation to the solution of the spatial problem of the longitudi-nal-torsional transducer.
Keywords: successive loading method, nonlinear problem, linearized system of equations, longitudinal-transverse transducer, plain curved rod, nonlinear elastic characteristic
4
ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
#9(726) 2020
Рассмотрим численный метод, позволяющий исследовать напряженно-деформированное состояние продольно-поперечного преобразователя и получать его нелинейную упругую характеристику.
Цель работы — разработка математической модели для расчета напряженно-деформированного состояния плоского криволинейного стержня при больших перемещениях.
Для плоского криволинейного стержня (рис. 1) определим зависимость поперечного перемещения (направленного вдоль оси х1) их1 верхнего конца 2 стержня от продольного перемещения (сонаправленного с осью х 2) их 2 его нижнего конца 1.
Считаем, что нижний конец стержня перемещается строго вертикально вдоль оси х2, а его верхний конец смещается в направляющих строго горизонтально параллельно оси Х1 [1-3]. Также принимаем, что концы жестко защемлены, т. е. не поворачиваются относительно точек крепления. Массой стержня пренебрегаем.
Заданы следующие параметры: высота Н, начальный угол Р и угол разворота дуги ф.
Углы между локальным единичным вектором сонаправленным с глобальной осью х1, и естественным единичным вектором е1, направленным по касательной к дуге окружности, определяются следующими выражениями:
• в начале стержня
п а
«о =--Р;
• в конце стержня
х2 ш 2/ \ai ш \
т ш
h 1 еУ
V 1 J-77» \а0 V
*1
е2
Ф
Рис. 1. Схема для определения перемещений плоского криволинейного стержня
К
Высота
«i = ^ - (Р+Ф).
H = R sin(P + ф) - R sin Р,
(1)
где Я — радиус окружности.
Из соотношения (1) выражаем радиус окружности
Н
R = -
(2)
sin^+ф) -sinР
Запишем систему уравнений для плоского криволинейного стрежня [1, 4, 5] как частный случай общих уравнений [6], полагая
K10 = K20 = 0; Q3 = 0; Mi = M2 = 0;
#1 =#2 = 0; из = 0,
где K10, K20 — начальные кривизны оси стержня относительно главных осей; Q3 — поперечная сила; M1 и M2 — крутящий и изгибающий момент; #1, #2 — углы поворота поперечного сечения относительно главных осей; и3 — линейное перемещение.
При этом полагаем, что P0 = AP = 0 и T0 = AT = 0 [6], где P0 и AP — внешняя сосредоточенная сила и ее малое приращение; T0 и AT — внешний сосредоточенный момент и его малое приращение.
Окончательно получаем систему из шести уравнений
dQ1
de dQ2 de dM 3
de
- K 30Q2 = 0; + K 30Q1 = 0; ■ + Q2 = 0;
4M 3 = 0;
(3)
de A3_
du1 v n
—--K 30U2 = 0;
de
du2 de
+ K 30U1 -#3 = 0,
где и Q2 — внутренние силы, направленные вдоль осей х1 и х 2; е — безразмерная осевая координата; К30 — начальная кривизна оси стержня; М3 — изгибающий момент; — угол поворота поперечного сечения стержня; А33 — изгибная жесткость стержня; и1 и и2 —
#9(726) 2020
ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
5
продольное и поперечное перемещение поперечного сечения стержня.
Одним из методов решения подобных нелинейных задач является метод последовательных нагружений [3, 5, 7], в котором каждое дискретное увеличение нагрузки описывается линейными (линеаризованными) уравнениями. Для каждого т-го шага внешняя сила р(т) = кр, где т — номер шага нагружения; к — параметр, определяющий часть общей нагрузки на каждом этапе, к = 1/п (п — количество шагов нагружения).
Система уравнений для первого шага нагружения совпадает с системой уравнений (3):
de
--K 3OQ2(1) = 0;
dQ
,(i)
de
dM3(1)
- + K 3OQI(1) = 0;
de
d#3(1)__i
de A33
+ Q2(1) = 0;
3 - — M з(1) = 0;
dui(1) de du2(1) de
- K 3OU2(1) = 0;
+ K3OUI(1) -#З(1) = 0.
Для каждого последующего т-го шага нагружения система уравнений будет выглядеть следующим образом [6, 8, 9]:
dO^(т)
-— к з(т-1)д2(т) - д2(т-1) АК з(т) = 0;
de
dQ
(m)
de
dM з(т)
de
d#3(m)__i
de A33
- + K 3(m-1)Qi(m) + Qi(m-1)AK з(т) = 0;
+ Q2(m) = 0;
3 —— M 3(m) = 0;
dui(m)
de du2(m) de
- K 3(m-1)u2(m) = 0;
+ K3(m-1)ui(m) -#3(m) = 0,
где
m-1 m-1 M (k)
K 3(m-1) = K 30 + E AK 3(k) = K 30 + E
k=2 k =2 A33
M, (m) m-1
AK 3(m) = Q(m-1) = E AQj (k).
A33 k=2
Рассмотрим граничные условия для очередного шага решаемой задачи.
При безразмерной осевой координате е = 0.
Так как нижний конец стержня жестко защемлен, угол поворота в начале участка равен нулю.
В начале координат задано перемещение вдоль оси х 2. Для решения задачи его необходимо спроецировать на естественные оси [10-12].
Матрица поворота для перехода от декартовой системы координат к естественной имеет вид
n cos ф sin ф
Цф) =
-Sin ф cos ф
С: I=-1-} •
Тогда граничные условия будут выглядеть следующим образом:
= 0; u1 =Au0sin а0; u2 =Au0cos а0,
где Au0 — малое приращение перемещения на конце стержня.
При безразмерной осевой координате е = 1.
Верхний конец стержня жестко защемлен, поэтому угол поворота в конце участка равен 0.
Верхний конец стержня может перемещаться только горизонтально в направляющих, поэтому можно записать следующее условие:
uxi = 0.
Сила, направленная вдоль оси х1, также должна быть равна нулю, поскольку верхний конец стержня может свободно перемещаться в направлении оси х1:
Qxi = 0.
Эти два граничных условия записаны в декартовой системе координат. Спроецировав их на естественные оси в конце стержня, получаем следующие алгебраические выражения [13-15]:
• для перемещения
ux1 = u1 sin a1 + u2 cos a1;
• для силы
Qx1 = Q1 sin a1 - Q2 cos a1.
Тогда граничные условия будут выглядеть следующим образом:
= 0; u1sin a1 + u2 cos a1 = 0;
Q1 sin a1 - Q2 cos a1 = 0.
6
ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
#9(726) 2020
их2, мм
5
0 I1 , , , , , ,
0 10 20 30 40 50 ихЪ мм
Рис. 2. Графики деформирования оси стержня: 1 — исходная форма оси стержня; 2 — форма оси стержня на последнем шаге нагружения
продольно-поперечного преобразователя
Полученные граничные условия должны выполняться на каждом шаге решения методом последовательных нагружений в программном комплексе МАТЬАБ, так как углы а0 и а1 остаются неизменными на протяжении всего решения задачи [15-17].
В качестве примера приведем результаты расчета плоского криволинейного стержня со следующими входными параметрами: H = 50 мм; Р = л/6; ф = л/6; ширина сечения Ь = 5 мм; высота сечения h = 1 мм; модуль Юнга E = 2 -105МПа.
Радиус кривизны рассчитываемого криволинейного стержня, вычисленный по формуле (2), Я = 0,1366 м.
В ходе решения задачи в математическом пакете МАТЬАБ получены графики деформирования оси стержня на каждом шаге нагруже-ния (рис. 2).
Также определена зависимость поперечного смещения верхнего конца стержня их\ от продольного перемещения его нижнего конца и0 (рис. 3), называемая нелинейной упругой характеристикой продольно-поперечного преобразователя.
Выводы
1. Методом последовательных нагружений в математическом комплексе МАТЬАБ проведен численный расчет напряженно-деформированного состояния плоского криволинейного стержня при больших перемещениях.
2. На основе полученных результатов построена нелинейная упругая характеристика продольно-поперечного преобразователя.
Литература
[1] Жилин П.А. Прикладная механика. Теория тонких упругих стержней. Санкт-Петер-
бург, Изд-во СПБГПУ, 2007. 101 с.
[2] Кудряшев С.Б. Решение задачи анализа и синтеза продольно-крутильного ультразву-
кового волновода. Молодой исследователь Дона, 2018, № 4(13), с. 74-76. URL: http://mid-journal.ru/upload/iblock/214/15_760_Kudryashev_74_76.pdf
[3] Grigoryev Y.V., Pyae Phyo Aung. Waveguide Analysis for Ultrasonic Medical Instruments
(Comparison of Different Methods). MATEC Web of Conferences, 2018, vol. 226, no. 04007, doi: https://doi.org/10.1051/matecconf/201822604007
[4] Бидерман В.Л. Механика тонкостенных конструкций. Статика. Москва, URSS, 2017.
485 с.
[5] Левин В.Е., Пустовой Н.В. Механика деформирования криволинейных стержней. Ново-
сибирск, Изд-во НГТУ, 2008. 208 с.
[6] Светлицкий В.А., Нарайкин О.С. Упругие элементы машин. Москва, Машиностроение,
1989. 380 с.
#9(726) 2020
ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
7
[7] Наумов А.М. Применение метода последовательных нагружений при решении задач
механики плоских стержней. Известия высших учебных заведений. Машиностроение, 2016, № 12, с. 33-42, doi: 10.18698/0536-1044-2016-12-33-42
[8] Гаврюшин С.С., Барышникова О.О., Борискин О.Ф. Численный анализ элементов кон-
струкций машин и приборов. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2014. 479 с.
[9] Елисеев В.В. Механика упругих тел. Санкт-Петербург, Изд-во СПБГПУ, 2003. 336 с.
[10] Khaled Alnefaie. Lateral and longitudinal vibration of a rotating flexible beam coupled with torsional vibration of a flexible shaft. International journal of mechanical and mechatronics engineering, 2013, vol. 7, no. 5, pp. 884-891.
[11] Svetlitsky V.A. Engineering Vibration Analysis: Worked Problems. B. 1. Berlin, Springer, 2004. 316 p.
[12] Svetlitsky V.A. Engineering Vibration Analysis: Worked Problems. B. 2. Berlin, Springer, 2010. 343 p.
[13] Феодосьев В.И. Сопротивление материалов. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2010. 590 с.
[14] Cardoni A., Harkness P., Lucas M. Ultrasonic Rock Sampling Using Longitudinal-Torsional Vibrations. Physics Procedia, 2010, no. 3, pp. 125-134, doi: 10.1016/j.phpro.2010.01.018
[15] Zhou G., Zhang Y., Zhang B. The Complex-mode Vibration of Ultrasonic Vibration Systems. Ultrasonics, 2002, vol. 40, pp. 907-911, doi: 10.1016/s0041-624x(02)00224-x
[16] Svetlitsky V.A. Statics of Rods. Berlin, Springer, 2004. 308 p.
[17] Svetlitsky V.A. Dynamics of Rods. Berlin, Springer, 2004. 340 p.
References
[1] Zhilin P.A. Prikladnaya mekhanika. Teoriya tonkikh uprugikh sterzhney [Applied mechanics.
Theory of thin elastic rods]. Sankt-Petersburg, SPBGPU publ., 2007. 101 p.
[2] Kudryashev S.B. Solution to the problem of analysis and synthesis of longitudinal-torsion ul-
trasonic waveguide. Molodoy issledovatel' Dona, 2018, no. 4(13), pp. 74-76 (in Russ.). Available at: http://mid-journal.ru/upload/iblock/214/15_760_Kudryashev_74_76.pdf
[3] Grigoryev Y.V., Pyae Phyo Aung. Waveguide Analysis for Ultrasonic Medical Instruments
(Comparison of Different Methods). MATEC Web of Conferences, 2018, vol. 226, no. 04007, doi: https://doi.org/10.1051/matecconf/201822604007
[4] Biderman V.L. Mekhanika tonkostennykh konstruktsiy. Statika [Mechanics of thin-walled
structures. Statics]. Moscow, URSS publ., 2017. 485 p.
[5] Levin V.E., Pustovoy N.V. Mekhanika deformirovaniya krivolineynykh sterzhney [The me-
chanics of deformation of curved rods]. Novosibirsk, NSTU publ., 2008. 208 p.
[6] Svetlitskiy V.A., Naraykin O.S. Uprugiye elementy mashin [Elastic elements of machines].
Moscow, Mashinostroyeniye publ., 1989. 380 p.
[7] Naumov A.M. The Application of the Method of Successive Loadings when Solving Prob-
lems of Mechanics of Flat Rods. Proceedings of Higher Educational Institutions. Machine Building, 2016, no. 12, pp. 33-42 (in Russ.), doi: 10.18698/0536-1044-2016-12-33-42
[8] Gavryushin S.S., Baryshnikova O.O., Boriskin O.F. Chislennyy analiz elementov konstruktsiy
mashin i priborov [Numerical analysis of structural elements of machines and devices]. Moscow, Bauman Press, 2014. 479 p.
[9] Eliseyev V.V. Mekhanika uprugikh tel [The mechanics of elastic bodies]. Sankt-Petersburg,
SPBGPU publ., 2003. 336 p.
[10] Khaled Alnefaie. Lateral and longitudinal vibration of a rotating flexible beam coupled with torsional vibration of a flexible shaft. International journal of mechanical and mechatronics engineering, 2013, vol. 7, no. 5, pp. 884-891.
[11] Svetlitsky V.A. Engineering Vibration Analysis: Worked Problems. B. 1. Berlin, Springer, 2004. 316 p.
[12] Svetlitsky V.A. Engineering Vibration Analysis: Worked Problems. B. 2. Berlin, Springer, 2010. 343 p.
[13] Feodos'yev V.I. Soprotivleniye materialov [Strength of materials]. Moscow, Bauman Press, 2010. 590 p.
8
ИЗВЕСТИЯ ВЫСШИХ УЧЕБНЫХ ЗАВЕДЕНИЙ. МАШИНОСТРОЕНИЕ
#9(726) 2020
[14] Cardoni A., Harkness P., Lucas M. Ultrasonic Rock Sampling Using Longitudinal-Torsional Vibrations. Physics Procedia, 2010, no. 3, pp. 125-134, doi: 10.1016/j.phpro.2010.01.018
[15] Zhou G., Zhang Y., Zhang B. The Complex-mode Vibration of Ultrasonic Vibration Systems. Ultrasonics, 2002, vol. 40, pp. 907-911, doi: 10.1016/s0041-624x(02)00224-x
[16] Svetlitsky V.A. Statics of Rods. Berlin, Springer, 2004. 308 p.
[17] Svetlitsky V.A. Dynamics of Rods. Berlin, Springer, 2004. 340 p.
Информация об авторах
ПЬЯ Пьо Аунг — аспирант кафедры «Прикладная механика». МГТУ им. Н.Э. Баумана (105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1, e-mail: [email protected]).
ГРИГОРЬЕВ Юрий Всеволодович — кандидат технических наук, доцент кафедры «Прикладная механика». МГТУ им. Н.Э. Баумана (105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1, e-mail: [email protected]).
Статья поступила в редакцию 08.04.2020 Information about the authors
PYAE Phyo Aung — Postgraduate, Department of Applied Mechanics. Bauman Moscow State Technical University (105005, Moscow, Russian Federation, 2nd Baumanskaya St., Bldg.5, Block 1, e-mail: [email protected]).
GRIGORIEV Yuri Vselovolodovich — Candidate of Science (Eng.), Department of Applied Mechanics. Bauman Moscow State Technical University (105005, Moscow, Russian Federation, 2nd Baumanskaya St., Bldg.5, Block 1, e-mail: [email protected]).
Просьба ссылаться на эту статью следующим образом:
Пья П.А., Григорьев Ю.В. Численный метод расчета нелинейной упругой характеристики продольно-поперечного преобразователя. Известия высших учебных заведений. Машиностроение, 2020, № 9, с. 3-8, doi: 10.18698/0536-1044-2020-9-3-8
Please cite this article in English as: Pyae Phyo Aung, Grigoriev Y.V. A Numerical Method for Calculating the Nonlinear Elastic Characteristic of Longitudinal-Transverse Transducers. BMSTU Journal of Mechanical Engineering, 2020, no. 9, pp. 3-8, doi: 10.18698/0536-1044-2020-9-3-8
В Издательстве МГТУ им. Н.Э. Баумана вышло в свет учебное пособие К.Ю. Вуколова, А.М. Зимина, В.И. Тройнова
«Спектроскопия водородной плазмы»
В учебном пособии приведены теоретические основы атомной и молекулярной спектроскопии водородной плазмы, методики спектроскопических измерений и описание аппаратуры, используемой для диагностики высокотемпературной плазмы.
Предназначено для формирования компетенций, предусмотренных основными профессиональными образовательными программами по направлениям подготовки 16.04.02 «Высокотехнологические плазменные и энергетические установки» (магистратура) и 03.06.01 «Физика и астрономия» (аспирантура, уровень подготовки кадров высшей квалификации). Пособие также может быть использовано при подготовке специалистов, обучающихся по специальности 24.05.02-05 «Проектирование авиационных и ракетных двигателей» (специализация «Проектирование электроракетных двигателей»).
По вопросам приобретения обращайтесь:
105005, Москва, 2-я Бауманская ул., д. 5, стр. 1. Тел.: +7 499 263-60-45, факс: +7 499 261-45-97; [email protected]; www.baumanpress.ru
