Научная статья на тему 'Динамика стержня при продольном ударе телом'

Динамика стержня при продольном ударе телом Текст научной статьи по специальности «Физика»

CC BY
436
99
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
СТЕРЖЕНЬ / УДАР / ТЕОРИИ СЕН-ВЕНАНА / ГЕРЦА И СИРСА / ВРЕМЯ СОУДАРЕНИЯ / УДАРНЫЙ ИМПУЛЬС / ПАРАМЕТРИЧЕСКИЙ РЕЗОНАНС / ROD / IMPACT / SAINT-VENANT THEORY / HERTZ THEORY / SEARS THEORY / IMPACT TIME / IMPACT IMPULSE / PARAMETRIC RESONANCE

Аннотация научной статьи по физике, автор научной работы — Беляев Александр Константинович, Ма Чин-чин, Морозов Никита Федорович, Товстик Петр Евгеньевич, Товстик Татьяна Петровна

Рассматривается продольный удар упругим телом по концу упругого стержня с закрепленным противоположным концом. Учитывается распространение упругих волн в стержне и локальные деформации в зоне контакта. После отскока тела стержень совершает свободные продольные колебания, которые при определенных условиях могут вызвать параметрические поперечные колебания, имеющие характер биений. В зависимости от параметров задачи теоретически и экспериментально определяются время соударения, форма ударного импульса и максимальная амплитуда поперечных колебаний при параметрическом резонансе. Обсуждаются результаты экспериментальных исследований, проведенных в Национальном Тайваньском университете. Библиогр. 12 назв. Ил. 6. Табл. 1.

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

Похожие темы научных работ по физике , автор научной работы — Беляев Александр Константинович, Ма Чин-чин, Морозов Никита Федорович, Товстик Петр Евгеньевич, Товстик Татьяна Петровна

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

DYNAMICS OF ROD UNDER AXIAL IMPACT BY A BODY

A longitudinal elastic impact a body makes on a thin rod is considered. The results of theoretical, finite element and experimental approaches to the solution of the problem are compared. The theoretical approach takes into account both the propagation of longitudinal waves in the rod and the local deformation due to the Hertz model. This approach leads to a differential equation with a delayed argument. A three-dimensional dynamic problem is considered within the framework of finite-element approach, in which the wave propagation and local deformation are automatically taken into account. A benchmark of these two approaches showed a complete qualitative and satisfactory quantitative agreement of the results in the contact force and the impact time.In the experiment, only the impact time was determined. Comparison of the measured impact time with the theoretical and finite-element results was satisfactory. Due to the fact that the tested rod was relatively short, an approximate model with two degrees of freedom was also developed to calculate the impact time.The problem of excitation of transverse oscillations after the impactor rebound off the rod is solved. For parametric resonance, the motion has the character of beats at which the energy of longitudinal oscillation is transferred into the energy of transverse oscillation and vice versa. An estimate for the maximum possible amplitude of transverse oscillation is obtained. Refs 12. Figs 6. Table 1.

Текст научной работы на тему «Динамика стержня при продольном ударе телом»

УДК 539.3 Вестник СПбГУ. Математика. Механика. Астрономия. 2017. Т. 4 (62). Вып. 3

МЯО 74М20

ДИНАМИКА СТЕРЖНЯ ПРИ ПРОДОЛЬНОМ УДАРЕ ТЕЛОМ*

А. К. Беляев1, Ч.-Ч. Ма2, Н. Ф. Морозов3,

П. Е. Товстик3, Т. П. Товстик1, А. О. Шурпатов4

1 Институт проблем машиноведения Российской Академии Наук, Российская Федерация, 199178, Санкт-Петербург, Большой пр. В. О., 61

2 Национальный Университет Тайваня,

Тайвань, 10617, Тайбэй, Рузвельт роад, Секция 4, 1

3 Санкт-Петербургский государственный университет,

Российская Федерация, 199034, Санкт-Петербург, Университетская наб., 7—9

4 Санкт-Петербургский политехнический университет Петра Великого, Российская Федерация, 195251, Санкт-Петербург, Политехническая ул., 29

Рассматривается продольный удар упругим телом по концу упругого стержня с закрепленным противоположным концом. Учитывается распространение упругих волн в стержне и локальные деформации в зоне контакта. После отскока тела стержень совершает свободные продольные колебания, которые при определенных условиях могут вызвать параметрические поперечные колебания, имеющие характер биений. В зависимости от параметров задачи теоретически и экспериментально определяются время соударения, форма ударного импульса и максимальная амплитуда поперечных колебаний при параметрическом резонансе. Обсуждаются результаты экспериментальных исследований, проведенных в Национальном Тайваньском университете. Библиогр. 12 назв. Ил. 6. Табл. 1.

Ключевые слова: стержень, удар, теории Сен-Венана, Герца и Сирса, время соударения, ударный импульс, параметрический резонанс.

1. Введение. Рассматриваются динамические процессы, происходящие при продольном ударе телом по стержню. Задача соударения упругих тел при учете их деформации имеет богатую историю. Сен-Венан [1] с учетом распространения продольных волн рассмотрел продольное соударение стержней. Оказалось, что время соударения, найденное теоретически, плохо подтверждается экспериментами. Это связано с трудновыполнимым в эксперименте условием, при котором концы стержней были бы идеально плоскими. В то же время Герц [2], опираясь на исследования Буссинеска [3] по деформации упругого полупространства, решил задачу о прямом центральном соударении упругих тел (типа шаров). Учитываются только местные статические деформации, а распространение волн не рассматривается. При этом совпадение теоретического и экспериментального времен соударения оказалось хорошим. Сирс [4] объединил подходы Сен-Венана и Герца и рассмотрел соударение стержней с закругленными концами. При этом учитывались и местные деформации, и распространение волн. Этот подход приводит к хорошему совпадению теоретических и экспериментальных результатов и используется в многочисленных последующих работах [5, 6].

После отскока ударника стержень совершает свободные продольные колебания, которые в случае параметрического резонанса возбуждают поперечные колебания, имеющие характер биений. В работах [7-9] исследовано возбуждение поперечных колебаний, вызванных продольным ударом, причем ударная сила Р(£) была заданной функцией времени. В настоящей работе сила Р(£) находится в процессе решения задачи о продольном ударе телом по стержню. Получена оценка амплитуды биений при главных параметрических резонансах.

* Работа выполнена при поддержке РФФИ (гранты №16.51.52025 МНТ-а, № 16.01.00580-а). (¡5 Санкт-Петербургский государственный университет, 2017

Обсуждаются результаты экспериментальных исследований, проведенных в Национальном Тайваньском университете.

2. Продольный удар телом по стержню. Рассматривается упругий стержень длиной Ь, правый конец которого закреплен (рис. 1). К левому концу в начальный момент времени со скоростью уо подлетает ударник массой т, в результате чего в системе возникает контактное взаимодействие. Ударник считается абсолютно твердым телом, однако учитывается местное смятие а в зоне контакта. По стержню со скоростью а распространяются плоские продольные волны. Сила тяжести не учитывается. Определяется величина контактной силы Р(£) и время соударения.

Рис. 1. Удар телом по стержню.

Это классическая задача удара тела по стержню. Запишем условие совпадения координат точек контакта [5, 6, 10]:

иа-а = и, иа = г>о1--— / [1-в)Р{в)с1в, и= [ Р{в)У{1 - в)с1в. (2.1)

т 1 о 3 о

Здесь ио — перемещение ударника, и — перемещение левого конца стержня,

У (Ь + 4Ь/с) = У (*),

(2.2)

_ _а_ ( 1, 0 < I < 2Ь/с, ( ] ~ ЕЕ { -1, 2Ь/с < * < 4Ь/с,

Р = Нс?!\ +

V 4 Е 4Ео ) У

где У(£) —реакция стержня на приложенный к нему единичный импульс, контактная сила Р(£) связана с местным смятием а по формуле Герца [2], Ео,Е и vo, V — модули Юнга и коэффициенты Пуассона ударника и стержня соответственно, Е — площадь поперечного сечения стержня, с = \fWJ~p — скорость звука в стержне, р — плотность материала стержня, к — сумма кривизн соударяющихся тел в точке контакта (предполагается, что в зоне контакта тела осесимметричны). Если тело ударяет по плоскому концу стержня и материал соударяющихся тел одинаковый, то

Р = ка3/2,

к

2 Ет/Ё

3(1 — г/2)'

(2.3)

где Е — радиус кривизны ударника в точке контакта.

В результате (2.1) становится нелинейным интегральным уравнением относительно силы Р (£) вида

— [ Р{6){1-6)с16+ [ Р{в)У{1-т ]о ./о

+

к У

2/3

уоЬ.

(2.4)

Удар заканчивается при обращении силы Р(£) в нуль.

V

Перейдем к безразмерным переменным, считая длину стержня и скорость звука в нем, равными единице, и отнеся силу P к ЕЕ. Тогда уравнение (2.4) примет вид

£ / P(Ti)(T-Ti)dn + / P(Ti)Y(T—Ti)dn +гР(т)2/3 = «ст, Y(T) = ( —1)[т/2], (2.5) J о ./с

где

Р Р at ^ щ М 1 /£Р\2/3

Здесь через [z] обозначена целая часть числа z, параметр £ равен отношению массы стержня M = pFL и массы ударника m, параметр r учитывает местное смятие.

Уравнение (2.5) содержит три безразмерных параметра £, V0 и r. При r = 0 местное смятие не учитывается, уравнение (2.5) становится линейным, время удара зависит только от £, а сила Р(т) зависит от £ и линейно зависит от Vo. При r = 0 уравнение (2.5) имеет аналитическое решение, построенное, например, в [1, 6].

При r > 0 путем дополнительного изменения масштаба три параметра £, Vo, r можно свести к двум. Положим Р = r3Р, Vo = r3«0. Тогда в уравнении (2.5) величины Р и Vo заменятся на Р и го = «о/r3, а r заменится на 1, т. е. в уравнении (2.5) останется два параметра V0 и £.

Вернемся к уравнению (2.5). Дифференцирование по т дает [6]

^ = Р(пУ1т1-Р(т) +2Р(т-2)-2Р(т-4) + .. )j . (2.7)

Уравнение (2.7) интегрируем при условии Р(т) = 0 при т < 0. Чтобы уйти от тривиального решения Р(т) = 0 при численном интегрировании, следует задать Р(0) малым положительным числом, например 10-10. Интегрируем уравнение (2.7) последовательно на промежутках (0, 2), (2,4), (4, 6) и т. д., пользуясь тем, что на предыдущем промежутке значение функции Р(т) уже найдено. При этом учитываем непрерывность функции Р при т = 2, 4,... Заметим, что в пакете Wolfram Mathematica-9 возможно непосредственное интегрирование уравнения (2.7) без разбиения времени на интервалы.

В зоне контакта характер распространения плоской продольной волны нарушается. В работах [4, 6] влияние краевого эффекта в зоне контакта учитывается в рамках предположения о том, что плоская волна деформации распространяется по стержню, начиная с некоторого удаления от его конца. Это приводит к тому, что слагаемое Р(т) в уравнении (2.7) заменяется на Р(т — S), где S — величина порядка толщины стержня [6]. Нам представляется, что такой подход нуждается в уточнении, в основном, из-за того, что не учитываются силы инерции участка (0, S) стержня. В дальнейшем считаем S = 0.

Выше функция Y(t) и уравнение (2.7) получены на основании анализа распространения волн по стержню. Те же результаты могут быть получены методом Фурье. Рассмотрим краевую задачу (в переменных (2.6))

сРи д2и ди

дт2 dx2 ' dx

= — Р(т), U(1,T )=0 (2.8)

x=0

с нулевыми начальными условиями. Ее решение имеет вид u(x, т) ^^ ifk(т) cos Vkx, Vk = (k — 0.5)n,

k= 1

d2

Vk . 2

dT 2

+ v2vk = 2P(T ), Vk (0) = 0,

dT

(2.9)

T=0

2 />т -

<Рк(т) = — Р(т1) Зш(1Ук(т — Т\))с1т1. ^к ] 0

Теперь перемещение левого конца стержня записывается в виде

((t)= u(0,t)=/ Р(п)У(т - Ti)dTi, У(т) = 2^

k=1

sin Vk T Vk

(2.10)

Ряд У(т) медленно сходится к разрывной функции У(т). На рис.2 показана сумма 20 первых слагаемых при 0 < т < 6.

Рис. 2. Функция Y(t).

3. Метод конечных элементов. Сравнение результатов. Рассмотрим продольный удар телом с радиусом закругления Д по плоскому концу стержня из того

же материала. Пусть L = 0.5 м, F = 5 • 10

-5 „2

R = 0.01 м, m = 0.5 кг, vo = 1 м/с, 7.8 • 103 кг/м3. Находим c = 5189 м/с, = v0/c = 0.0001927, С = M/m = 0.39,

Е = Е0 = 2.1 • 1011 Н/м2, V = = 0.3, р = р0 М = 0.195 кг и безразмерные параметры у0 г = 0.0155, входящие в уравнение (2.7).

Зависимость Р(т) в переменных (2.6) представлена на рис. 3. Кривая Р(т) имеет три максимума. Продолжительность удара примерно в 6 раз превосходит время пробега волны по длине стержня.

Для ударника и стержня с теми же параметрами были проведены расчеты для задачи в трехмерной постановке по методу конечных элементов (МКЭ) [10]. Было использовано 3000000 элементов, число степеней свободы около 1000000. Подробности выбора элементов приводятся в [10]. На рис. 4 представлено сравнение результатов расчетов по одномерной модели (2.7) и по МКЭ. Сравнение кривых на рис. 4 говорит о том, что одномерная модель дает более короткое время удара и большую величину максимальной ударной силы, чем МКЭ в трехмерной постановке.

4. Послеударные колебания стержня. Пусть отскок ударника произошел в момент времени т = Т. При т > Т стержень совершает свободные продольные колебания, которые при определенных условиях могут вызвать изгибные колебания.

0

Рис. 3. Зависимость Р(т) в безразмер- Рис. 4- Сравнение зависимостей Р(4) по

ных переменных. уравнению (2.7) — кривая 1, и по МКЭ — кри-

вая 2.

Запишем в переменных (2.6) систему уравнений, описывающих связанные продольно-поперечные колебания [8],

дт2

сРч дх2

d2w d дт2 дх

dw дх

+ м2

dw д2w дх дх2

дх4 0

0, —

ди дх

0, и(1,т ) = 0,

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

x=0

£ =

ди

дх'

д2г

дх2

= 0, х = 0,1.

(4.1)

(4.2)

Здесь и(х,т) и w(x,t)—продольное и поперечное перемещения оси стержня, fj, = \JJ/(FL2) — малый параметр толщины стержня, J — момент инерции поперечного сечения. Уравнение (4.2) записано для шарнирно опертых концов стержня. Начальные условия для функции и(х, т) при т = T совпадают с ее значениями в момент окончания удара. Имея в виду возможность появления параметрических колебаний, начальные условия для w^, т) берем ненулевыми исчезающе малыми.

В уравнении (4.2) е(х,т) —осевая деформация стержня. Пренебрегая в уравнении (4.1) нелинейным слагаемым, с учетом (2.9) находим выражение для перемещения точек стержня в виде ряда Фурье по системе собственных функций cos Vkх с безразмерными частотами Vk:

и(х, т)

w 2 I »

/ P(ri)sin(i/fc(r - ti))(1ti cosvkx = k=i VkJ 0

w w

= /,(bk cos Vk т + Ck sin Vk т) cos VkX = N. ak COs(vk т + ak ) cos Vk х, (4.3)

k=1

k=i

2

bk = ——I P(T)sm(vkT)dT, Vk 0

2

ck = — I P(t) cos(vkT)dT, Vk 0

ak

= л/b2k + ■

k

Следовательно, функция м(ж, т), а вместе с ней и функция е(ж, т) являются периодическими функциями т с периодом 4. Коэффициенты Ьк и Ск могут быть вычислены при интегрировании уравнения (2.7). Ниже при рассмотрении параметрических ре-зонансов используется амплитуда колебаний по первой форме ах. Начальная фаза «1 не используется, ибо за счет сдвига начала отсчета времени при послеударных колебаниях можно считать ах =0.

Условием возникновения параметрических резонансов для уравнения (4.2) является равенство [7-9, 11]

|wm ± wn| = kVi, m, n, k = 1, 2,. .., Vi = n/2, = n , (4.4)

w=

2

где wm — частоты поперечных колебаний по формам w(x) = sin mnx. Точнее, резо-нансы могут возбуждаться в окрестностях значений параметров, удовлетворяющих соотношению (4.4). Для главных параметрических резонансов, возбуждающихся, как правило, с наибольшей интенсивностью, будет

m = n, k = 1, 4npmm2 = 1, 1m = 4nm2, l = 1/p. (4.5)

Здесь введен большой параметр длины l, измеряемый в долях радиуса инерции поперечного сечения стержня. Формула (4.5) определяет критические длины стержня 1m, в окрестностях которых возбуждаются главные параметрические резонансы.

В линейном приближении при параметрических резонансах амплитуда поперечных колебаний неограниченно возрастает, поэтому в уравнение (4.1) введено нелинейное слагаемое, описывающее влияние поперечных колебаний на продольные. При этом колебания при параметрическом резонансе имеют характер биений с последовательной передачей энергии продольных колебаний поперечным и обратно [7-9]. С использованием двухмасштабных асимптотических разложений [12] при интегрировании системы (4.1), (4.2) для главных параметрических резонансов получена оценка амплитуды прогиба [8, 9]:

maxw(x,T) < 2ai, (4.6)

x ,т

где ai — амплитуда перемещения по первой форме продольных колебаний, определенная по формуле (4.3).

Обратимся к рассмотренному в п. 2 примеру. Параметр длины равен I = L\JF/ J, и при L = 0.5 м, F = 5 • 10-5 м2 для круглого стержня находим l = 250.7. Ближайшие критические длины равны I4 = 201.1, I5 = 314.2, поэтому для круглого стержня главный параметрический резонанс не возбуждается. Однако при тех же значениях L и F для стержня другого поперечного сечения (например, прямоугольного) возможно равенство l = I5, при этом резонанс возбуждается по форме w = wo(t)sin5nx. При интегрировании системы (2.7) по формулам (4.3) находим амплитуду продольных колебаний ai, а затем по формуле (4.6) — амплитуду биений max wo(t) = 4.27 • 10-5. Учитывая, что за единицу длины принято L = 0.5 м, находим, что амплитуда биений равна 0.0214 мм.

5. Сравнение теоретических и экспериментальных результатов. В Национальном Тайваньском университете была проведена серия экспериментов по определению времени соударения при продольном ударе с использованием пьезоэлектрического эффекта. Рассматривается ударник в виде стального шарика с теми же параметрами Eo, ро, vo, что и в примере п. 3, масса шарика m принимает 6 значений, приведенных в таблице. Был взят стержень с параметрами: длина L = 0.0301 м, площадь поперечного сечения F = 10-4п м2, модуль Юнга E = 7.342 • 1010 Н/м2, коэффициент Пуассона v = 0.34, плотность р = 2696.6 кг/м3. Скорость удара v0 = 0.3225 м/с. Находим скорость звука в стержне с = 5218 м/с и его массу M = 25.5 г.

В столбцах t3KC и ¿МКЭ таблицы в безразмерном виде приведены времена соударения, найденные экспериментально и по методу конечных элементов. В качестве единицы времени принято время пробега волны по длине стержня. Время в секундах получается умножением на L/c = 5.77 • 10-6с.

В столбце t(2.7) приведено время, полученное в результате решения уравнения (2.7), причем г»о = v0/с = 0.0000618, а безразмерные коэффициенты £ и r для разных значений m приведены в таблице.

Сравнение теоретических и экспериментальных значений времени удара

т (г) ^ЭКС ¿MIO ¿(2.7) tapp t г

0.13 6.31 5.41 4.66 4.70 196.15 2.477

0.44 8.12 10.53 6.99 7.04 57.95 2.163

0.71 9.06 13.54 8.23 8.25 35.91 2.051

1.04 10.41 14.15 9.32 9.33 24.52 1.966

3.51 14.73 18.41 14.04 14.05 7.26 1.718

8.46 19.26 21.10 18.84 18.84 3.01 1.556

На рис. 5 показан график функции Р(т) для т = 8.46 г. Для других значений т графики по форме такие же.

410

2-10

5 10 15

Рис. 5. Функция Р(т) для параметров при эксперименте.

Кривая на рис. 5 существенно отличается от кривых на рис. 2 и 3 тем, что не проявляется распространение волн по стержню. Она похожа на график функции, описывающей ударную силу при соударении двух шаров.

—► D, D

щ и

Рис. 6. Приближенная модель с двумя массами.

Это обстоятельство для описания результатов экспериментов позволяет приближенно заменить рассматриваемую систему системой с двумя степенями свободы (рис. 6), в которой Di —нелинейная пружина, моделирующая местное смятие согласно формуле (2.3), D —линейная пружина жесткостью EF/L, равной жесткости стержня на растяжение, а масса второго тела принимается равной половине массы стержня (mi = M/2). Вместо уравнения (2.7) приходим к уравнению

S = ¿Р(Т)1/3 (ío-íjTpínJdn^jTpínJoos^-n))^) . (5.1)

Время удара tapp, полученное в результате решения уравнения (5.1) при тех же начальных условиях, что и для уравнения (2.7), помещено в столбце tapp таблицы. Видим, что значения tapp и t(2.7) близки друг к другу, причем разница убывает вместе с ростом времени удара. Разумеется, для других значений параметров задачи, для которых безразмерное время удара является более коротким, модель, показанная на рис. 6, не применима.

Причины различия результатов в столбцах ¿экс, ¿МКЭ и í(2.r) таблицы подлежат дальнейшему изучению, причем одной из причин расхождения является влияние краевого эффекта, которое для короткого стержня относительно велико.

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

При эксперименте определялось только время соударения. Сравнение с теоретическим и конечно-элементным результатами было удовлетворительным. В связи с тем, что стержень был относительно коротким, для вычисления времени соударения была использована также приближенная модель с двумя степенями свободы.

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

Литература

1. Saint-Venant A. Sur le choc longitudinal de deux barres elastiques // J. de Math. (Liouville) Ser. 2. T. 12. 18б7.

2. Hertz H. Über die Berührung fester elastischer Körper // Za. f. Math. (Crelle). 1881. Bd 92.

3. Boussinesq J. Applications des potentials. Paris, 1885.

4. Sears J. E. On the longitudinal impact of metal rods with rounded ends // Proc. Cambridge Phil. Soc. 1908. Vol. 14.

5. Бидерман В. Л. Теория удара. M., 1952.

6. Зегжда С. А. Соударение упругих тел. СПб.: Изд-во С.-Петерб. ун-та, 1997. 316 с.

7. Морозов Н. Ф., Товстик П. Е. Динамика стержня при продольном ударе // Вестн. С.-Петерб. ун-та. Сер. 1. 2009. Вып. 2. С. 105-111.

8. Беляев А. К., Морозов Н. Ф., Товстик П. Е., Товстик Т. П. Биения в задаче о продольном ударе по тонкому стержню // Изв. РАН. МТТ. 2015. № 4. С. 104-117.

9. Беляев А. К., Морозов Н. Ф., Товстик П. Е., Товстик Т. П. Параметрические резонансы в задаче о продольном ударе по тонкому стержню // Вестн. С.-Петерб. ун-та. Сер. 1. 201б. Т. 3(б1), вып. 1. С. 77-94.

10. Беляев А. К., Ма Ч.-Ч., Шурпатов А. О. Полуаналитическое, конечно-элементное и экспериментальное определение контактной силы осевого соударения стержня и ударника // Науч. техн. ведомости СПбГТУ. Физ.-мат. науки. (в печати).

11. Якубович В. А., Старжинский В. М. Линейные дифференциальные уравнения с периодическими коэффициентами и их приложения. M.: Наука, 1972.

12. Боголюбов Н. Н., Митропольский Ю. А. Асимптотические методы в теории нелинейных колебаний. M.: Наука, 19б9.

Статья поступила в редакцию 22 февраля 2017 г.; рекомендована в печать 30 марта 2017 г.

Сведения об авторах

Беляев Александр Константинович — доктор физико-математических наук, профессор; [email protected]

Ма Чин-Чин — профессор; [email protected]

Морозов Никита Федорович — академик РАН, доктор физико-математических наук, профессор; [email protected]

Товстик Петр Евгеньевич — доктор физико-математических наук, профессор; [email protected] Товстик Татьяна Петровна — кандидат физико-математических наук, старший научный сотрудник; [email protected]

Шурпатов Анатолий Олегович — аспирант; [email protected]

DYNAMICS OF ROD UNDER AXIAL IMPACT BY A BODY

Alexander K. Belyaev1, Chien-Ching Ma2, Nikita F. Morozov3, Petr E. Tovstik3, Tatiana P. Tovstik1, Anatoly O. Shurpatov4

1 Institute for Problems in Mechanical Engineering RAS, Bolshoy pr. V.O., 61, St. Petersburg, 199178, Russian Federation; [email protected], [email protected]

2 National Taiwan University, No. 1, Sec. 4, Roosevelt Rd., Taipei 10617, Taiwan (R.O.C.); [email protected]

3 St. Petersburg State University, Universitetskaya nab., 7—9, St. Petersburg, 199034, Russian Federation; [email protected], [email protected]

4 Peter the Great St. Petersburg Polytechnic University, Polytechnicheskaya ul., 29, St. Petersburg, 195251, Russian Federation; [email protected]

A longitudinal elastic impact a body makes on a thin rod is considered. The results of theoretical, finite element and experimental approaches to the solution of the problem are compared. The theoretical approach takes into account both the propagation of longitudinal waves in the rod and the local deformation due to the Hertz model. This approach leads to a differential equation with a delayed argument. A three-dimensional dynamic problem is considered within the framework of finite-element approach, in which the wave propagation and local deformation are automatically taken into account. A benchmark of these two approaches showed a complete qualitative and satisfactory quantitative agreement of the results in the contact force and the impact time.

In the experiment, only the impact time was determined. Comparison of the measured impact time with the theoretical and finite-element results was satisfactory. Due to the fact that the tested rod was relatively short, an approximate model with two degrees of freedom was also developed to calculate the impact time.

The problem of excitation of transverse oscillations after the impactor rebound off the rod is solved. For parametric resonance, the motion has the character of beats at which the energy of longitudinal oscillation is transferred into the energy of transverse oscillation and vice versa. An estimate for the maximum possible amplitude of transverse oscillation is obtained. Refs 12. Figs 6. Table 1.

Keywords: rod, impact, Saint-Venant theory, Hertz theory, Sears theory, impact time, impact impulse, parametric resonance.

References

1. Saint-Venant A., "Sur le choc longitudinal de deux barres elastiques", J. de Math. (Liouville). Ser. 2 12 (1867).

2. Hertz H., "Über die Berührung fester elastischer Körper", Za. f. Math. (Crelle) 92 (1881) [in German].

3. Boussinesq J., Applications des potentials (Paris, 1885).

4. Sears J. E., "On the longitudinal impact of metal rods with rounded ends", Proc. Cambridge Phil. Soc. 14 (1908).

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

5. Biderman V.L., Theory of impact (Mashgiz, Moscow, 1952) [in Russian].

6. Zegzhda S.A., Collision of elastic bodies (Publishers of St. Petersburg University, 1997) [in Russian].

7. Morozov N. F., Tovstik P.E., "Rod dynamics under axial impact", Vestnik St. Petersburg University, Series 1, issue 2, 105—111 (2009) [in Russian].

8. Belyaev A. K., Morozov N. F., Tovstik P. E., Tovstik T. P., "Beating in the problem of axial impact on a thin rod", Mechanics of Solids 50(4) (2015).

9. Belyaev A. K., Morozov N. F., Tovstik P. E., Tovstik T. P., "Parametric Resonances in the Problem of Longitudinal Impact on a Thin Rod", Vestnik St. Petersburg University. Mathematics 49(1), 53—67 (2016).

10. Belyaev A. K., Ma Chien-Ching, Shurpatov A. O., "Semi-analytic, finite-element and experimantal determination of contact force of axial collision of rod and impactor", St. Petersburg Polytechnic University Journal: Physics and Mathematics (to appear in 2017) [in Russian].

11. Yakubovich V. A., Starzhinski V. M., Linear differential equations with periodic coefficients (Wiley, New York, 1975).

12. Bogoliubov N. N., Mitropolsky Yu. A., Asymptotic methods in the theory of non-linear oscillations (Gordon and Breach, New York, 1961).

Для цитирования: Беляев А. К., Ма Ч.-Ч., Морозов Н.Ф., Товстик П. Е., Товстик Т. П., Шурпатов А. О. Динамика стержня при продольном ударе телом // Вестник СПбГУ. Математика. Механика. Астрономия. 2017. Т. 4 (62). Вып. 3. С. 506-515. DOI: 10.21638/11701/spbu01.2017.312

For citation: Belyaev A. K., Ma Ch.-Ch., Morozov N. F., Tovstik P. E., Tovstik T. P., Shurpatov A. O. Dynamics of rod under axial impact by a body. Vestnik SPbSU. Mathematics. Mechanics. Astronomy, 2017, vol. 4(62), issue 3, pp. 506-515. DOI: 10.21638/11701/spbu01.2017.312

О международной научной конференции по механике

«ВОСЬМЫЕ ПОЛЯХОВСКИЕ ЧТЕНИЯ»

Международная научная конференция по механике «Восьмые Поляховские чтения» пройдет с 30 января по 2 февраля 2018 года в Санкт-Петербургском государственном университете и посвящается памяти В.В.Белецкого — выдающегося ученого, специалиста в области небесной механики, теории вращательных движений искусственных и естественных небесных тел. Владимир Васильевич был членом научного комитета «Поляховских чтений», внёс большой вклад в развитие секции «Динамика космического полёта».

История конференции восходит к 1997 г., когда она впервые была проведена в формате Всероссийской конференции в связи с 90-летием выдающегося ученого, профессора СПбГУ Н. Н. Поляхова (1906—1987). С 2003 года «Поляховские чтения» приобрели статус международной конференции и по традиции проводятся в Санкт-Петербурге один раз в три года. Тематика «Поляховских чтений» охватывает практически все фундаментальные направления современной механики и продолжает расширяться в соответствии с развитием науки.

Научные направления конференции

I. Теоретическая и прикладная механика. VI. Устойчивость и стабилизация

II. Динамика естественных механических и электромеханических

и искусственных небесных тел. систем.

III. Гидрогазодинамика. VII. Биомеханика.

IV. Физико-химическая аэромеханика. VIII. Физическая механика. V. Механика деформируемого твердого тела. IX. История механики.

Рабочие языки конференции — русский и английский.

Сборник тезисов докладов (на русском или на английском языках) будет опубликован в электронном виде к началу конференции. Труды конференции будут опубликованы в электронном виде (на английском языке) в периодическом издании AIP Conference Proceedings, входящем в Scopus и WoS Core Collection.

Адрес сайта конференции — https://events.spbu.ru/events/polyakhov_readings.

Электронный адрес Оргкомитета конференции — [email protected].

Председатель конференции: Никита Федорович Морозов (академик РАН).

Сопредседатели конференции: Г.А.Леонов (член-корр. РАН), Е. В.Кустова (СПбГУ), М.П. Юшков (СПбГУ).

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