Потеря устойчивости однородного стержня при продольном ударе о стержень, взаимодействующий с жесткой преградой Текст научной статьи по специальности «Физика»

Механика

УДК 622.233.6

ПОТЕРЯ УСТОЙЧИВОСТИ ОДНОРОДНОГО СТЕРЖНЯ ПРИ ПРОДОЛЬНОМ УДАРЕ О СТЕРЖЕНЬ,

ВЗАИМОДЕЙСТВУЮЩИЙ С ЖЕСТКОЙ ПРЕГРАДОЙ

А.А. Битюрин

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

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

Введение

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

Впервые задача о динамической устойчивости упругих стержней при мгновенном продольном действии сжимающей нагрузки, которая в дальнейшем сохраняет своё постоянное значение, была решена М.А. Лаврентьевым и А.Ю. Ишлинским в работе [7]. Позднее проводилось теоретическое исследование потери устойчивости мгновенно сжатого бесконечного упругого стержня методом малых возмущений [11], а также исследование потери устойчивости полубесконечного упругого стержня [5, 8]. В работах [9, 10] Малым В.И предложен подход для исследования потери устойчивости стержней при больших и малых прогибах на основе асимптотического метода линеаризации на неустановившейся стадии выпучивания при ударе полубесконечного упругого стержня телом большой массы, движущимся с постоянной скоростью.

Выполнение поставленных задач сводилось либо к аналитическому, либо к численному решениям линеаризованного дифференциального уравнения четвёртого порядка, а также системы дифференциальных уравнений движения стержня с начальной прогибью. Недостатком такого подхода является очевидная трудность математического решения, сопряженная с громоздкими математическими выкладками, а также удовлетворительная сходимость результатов математического моделирования и экспериментальных данных лишь в узком интервале изменения той или иной характеристики, связанная с накоплением погрешности вычислений в ходе проведения математических операций. Задача о потере устойчивости неоднородных ступенчатых стержней при ударе в перечисленных работах, а также в работах [2, 5, 12, 14] не рассматривалась.

Постановка задачи

Рассматривается математическая модель продольного упругого центрального удара однородного стержня меньшего поперечного сечения, движущегося со скоростью У0, по однородно-

му стержню, взаимодействующему с жесткой преградой. Масса первого стержня т], длина /,, масса и длина второго соответственно т2 и /2, общая длина обоих стержней 1=1х+12 (рис. 1). Оба стержня состоят из одного материала.

У,э

~>’ 2

1

т2

I — 1Х +/2

Рис. 1. Схема удара однородных стержней при неудерживающих связях

Движение поперечных сечений соударяемых стержней описывается волновыми уравнениями:

Т2 ~т дТ2 _и’ и;

дх а о/

д\(х,1) 1 д2и2(х,1) _Л

Я 2 2 “ и’ 1 ~ 1 ’ '

ох а от

где м,(х,0, и2(х,1)~ продольное перемещение поперечного сечения стержней 1 и 2, х- координата сечения, I - время, а - скорость распространения продольной волны деформации.

Начальные условия определяют состояние стержней перед их соударением:

при ,=,-=0 ^-‘оКу,.^'оКоаи2(х,1„) н 0 ы 0 дх д( дх

Краевые условия определяют отсутствие силы в сечении х = 0 и равенство нулю скорости сечения х = I при взаимодействии стержня 2 с жесткой преградой:

д»1(0.0 = 0 ди2(1,1 ) = 0

дх 5/

а также определяют равенство сил и условия сопряжения стержней в сечениях х = 1Х при непосредственном их взаимодействии:

= (5)

дх дх дх

дщ(}\^ )_ди2(1и0 с дщ(1х,( )

& 5/ ’ дх

либо отсутствие сил в ударных сечениях стержней, если их взаимодействие отсутствует:

д3£\:1\ = 0> (]112,(1\-0 = р; если ^ /)—¿¿2(/[, о <0, (7)

ох ох

где £ - модуль упругости первого рода, и _ площади поперечных сечений соответственно первого и второго стержней, а - скорость распространения продольной волны деформации.

Метод решения

Решение дифференциальных уравнений (1), (2) реализуется методом Даламбера в виде [1]:

И](х,0 = /1 (^-х) + <рх{ей + х), 0<х</], (8)

м2(х,0 = У2(о^-д') + 9?2(°^ + д:)» 1\^х<1х+12, (9)

= ~/{{а1-х) +(р{{М+х), (10)

дх

vx(x,t)=^~9-= a[f{(at-x)+<p{(at+x)], (11)

ot

e2{x,t) = du*M = -f2(at-x) + (p2(at + x) , (12)

dx

V2(x,t)=^ll= a[fHat-x)+(p’2{at + x)]. (13)

at

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

f’(at-x) = f'(at-x)l—; <p'(at + х) = <p'(at + х)I—, деформацию в сечении и его скорость

а а

s(x,t) = - f\at - х) + ф'{м + х), v(x,t) = v(x,t) = f'(at - х) + (p'(at + х).

VQ

Пример 1

Пусть в момент времени Т = tj построена диаграмма относительной продольной деформации ё для первого однородного участка произвольной многоступенчатой стержневой системы при продольном ударе (рис. 2, а).

/,

■if.....................................'................................../

Т- t.

ь,

1 h

р . m

/

jf

/

a)

P

-t "it '/

A

6)

Рис. 2. а) диаграмма относительной продольной деформации для одного из участков произвольной многоступенчатой стержневой системы; б) эквивалентная схема нагружения рассматриваемого однородного участка стержневой системы

Аналогичным путем можно построить диаграмму для любого однородного участка [15]. Величину продольной силы Р1 на первом однородном участке, в соответствии с законом Гука [13] и используя данные диаграммы, представим в виде

а

(14)

Для определения эквивалентной продольной нагрузки, приложенной в торце первого однородного участка (рис. 2, б), можно воспользоваться зависимостью (15), представленной в [13]:

Р =

ЭКв

V2

чАу

\2

V. А у

+...+Р

V2

чАу

+ Л

Ґ, л2 Ел.

чАу

(15)

В формулах (14)и(15)а- скорость распространения продольной волны деформации в данном материале стержней, Е - модуль упругости, А - площадь поперечного сечения рассматриваемого однородного участка, ё1 - относительная продольная деформация на /'-м отрезке рассматриваемого однородного участка многоступенчатого стержня, ЬХ,Ь2,Ь1,Ьп — координаты приложения продольных сил Р1 на этом участке.

Для тонких однородных участков (/ > 15с1, где I и с! - соответственно длина и диаметр рассматриваемого участка или стержня) при расчете на устойчивость применима формула Эйлера [6, 13]

Ркр=7Г2Е^/{М1)2, (16)

где Jmin - минимальный момент инерции поперечного сечения относительно главных центральных осей инерции, ¡л - коэффициент приведения длины, зависящий от схемы закрепления стержня [13], / - длина рассматриваемого однородного участка. В данном случае это будет длина /,, поскольку I обозначена общая длина ступенчатого стержня. Если рассматриваемый однородный участок жестко связан со вторым, то коэффициент приведения длины для первого однородного участка /и принимается равным 2 [13]. Аналогичным путем можно рассматривать любой однородный участок произвольной многоступенчатой стержневой системы и в зависимости от его связей с другими участками или с жесткой преградой /л будет принимать значения от 0,7 до 2 [6, 13].

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

Для выполнения условия сохранения устойчивого состояния многоступенчатой стержневой системы необходимо, чтобы Рэкв< Ркр , т. е. используя (14)—(16) и, учитывая, что для круглого

поперечного сечения Jm■m -

ЖСІ

~64~

Уо<КР

п2с12а

16(М)

Ґ г \2

V А

получим

\2

+

(¿2-Ю

V

А

+ (ё3-ё2)

/

ґьУ

(Ь-Ъ) Т +(ёп~Ъ)

V Ау

г А л2

к

V А у

-і

.(17)

Формула (17) применима для расчета критической предударной скорости,многоступенчатой стержневой системы с любым числом однородных участков.

Пример 2

Рассмотрим удар о жесткую преграду однородных стержней с длинами 1Х = 12 =0,5/. Соот-

А

ношение площадей поперечных сечений первого и второго стержней А =— = 0,5 (рис. 3).

А,

о

ч

1 2 '"'1 X

,

./ / т - /

и

Ш7

!1 к+Ь=1

Рис. 3. Схема удара однородного стержня меньшего поперечного сечения о стержень, взаимодействующий с жесткой преградой

0,251

0,51

0,751

0,51/а

1/а

1,51/а

21/а

2,51/а

31/а

/О =0,5 ФО =0,5

¿О =0

v0= 1

/Г 0,5 Ф1= -0,17

= -0,17 -0.17 i2 = 0

-0,34

Sj = -0,67 V|= 0,33

""/l=-0,17

<pi=-0,4

/1 = Ф1=

fl=

v,=

-0,4

-0,4

8Г -0,23 . vj= -0,57

I4

= 0 • -0,8

... ' /Г--0.4

Is Ф1= -0,4

:Я=_:М_ Ii=q_

По

II,

/о = ° ФО =0

/2= 0,34 Ф2=0

£2 = О $2 = О

II 2 "

S2 = -0,34 V2= 0,34

- II3

’/¿=0,34 Ф2= -0,34

s2= -0,68 v2=0

Л = -0,22 Ф2= -0,34 II4

£2= -0,12 v^= -0,56

'"?2= -0,22 I! 5 92=-0,22

н„

/1=-0,22' Ф2=-0,22 62=0

§2=0 V2= -0,44

v2= -0,44

Рис. 4. Поле состояний при ударе однородного стержня меньшего поперечного сечения о стержень,

взаимодействующий с жесткой преградой

Применим метод характеристик для построения поля состояний (рис. 4). Области состояний 10 - 15, По - Пб с соответствующими значениями /\at-x ), ф'(си + х ), ё(х,(), у(х,1) определяют параметры прямых и обратных волн деформаций, продольную деформацию и скорость поперечных сечений. Длительность состояния для произвольного сечения определяется разностью ординат /, которые имеют точки наклонных линий для рассматриваемого сечения.

По данным поля состояний (рис. 4) строим диаграмму относительной продольной деформации ё (рис. 5).

0,51/а

1/а

1,51/а

21/а

2,51/а

і

1 2

1“':!:' - -

1

0 0,251 0, 51 0,751 1

е~ 0

Г

-0,67 ^ 0 -0,34

еГ

-0,68

-0,23 І1 0 -0,12

Рис. 5. Диаграмма относительной продольной деформации

Анализируя диаграмму, можно сделать вывод, что максимальное значение относительной продольной деформации на первом стержне имеет место при I=1,0//а (=-0,67 ). На втором стержне ётт =-0,68 при t-l,5l/a . Следовательно, в данные моменты времени при достижении предударной скорости некоторого критического значения стержневая система может потерять устойчивость. Для расчета критической предударной скорости используем полученную формулу (17), в которую подставляем одно единственное значение относительной продольной деформации для каждого стержня: сначала £тах=-0,67 для первого, затем £тах=-0,68 для второго, получив два значения критической скорости. Учитывая коэффициент приведения длины ц=\, соответствующий неудерживающим связям в сечении их взаимодействия и в сечении, взаимодействующим с жесткой преградой, получим соответствующие значения предударной

скорости ¥кр =0,37

п2й2а

. ж2сі2а

2 для первого и Укр =0,36 для второго стержней. Анализируя по-

лученный результат, заключаем, что потеря устойчивости рассматриваемой стержневой системы наступит при /=1,0На, поскольку в данный момент времени достигается критическое значение предударной скорости для первого стержня.

Заключение

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

Литература

1. Алимов, О.Д. Распространение волн деформаций в ударных системах / О.Д. Алимов, В.К. Манжосов, В.Э. Еремьянц. - М.: Наука, 1985. - 354 с.

2. Алфутов, Н.А. Основы расчета на устойчивость упругих систем / Н.А. Алтуфов. - М.: Машиностроение, 1978. - 312 с.

3. Битюрин, А.А. Моделирование продольного удара однородных стержней при неудерживающих связях / А.А. Битюрин, В.К. Манжосов // Вестник УлГТУ. - 2005. - № 3. - С. 23-25.

4. Вольмир, А.С. Исследование процесса выпучивания стержней при ударе / А.С. Вольмир, И.Г. Кильдибеков // Доклады АН СССР. - 1966. - Т. 2. - Вып. 10. - С. 10-17.

5. Вольмир, А.С. Устойчивость упругих систем / А.С. Вольмир. - М.: ГИТТЛ, 1962.- 880 с.

6. Дарков, А.В. Сопротивление материалов / А.В. Дарков, Г.С. Шпиро. - М.: Высшая школа, 2003.-641 с.

7. Лаврентьев, М.А. Динамические формы потери устойчивости упругих систем / М.А. Лаврентьев, А.Ю. Ишлинский // Доклады АН СССР. - 1949. - Т. 65, № 6. - С. 42-45.

8. Малый, В.И. Длинноволновое приближение в задачах о потере устойчивости при ударе / В.И. Малый // Известия АН СССР. МТТ. - 1972. - № 4. - С. 138-144.

9. Малый, В.И. Выпучивание стержня при продольном ударе. Малые прогибы / В.И. Малый // Известия АН СССР. МТТ. - 1973. - № 4. - С. 181-186.

10. Малый, В.И. Выпучивание стержня при продольном ударе. Большие прогибы / В.И. Малый // Известия АН СССР. МТТ. - 1975. - № 1. - С. 52-61.

11. Малышев, Б.М. Устойчивость стержня при ударном сжатии / Б.М. Малышев // Известия АН СССР. МТТ. - 1966. - № 4. - С. 137-142.

12. Пановко, Я.Г. Устойчивость и колебания упругих систем / Я.Г. Пановко, И.И. Губанова. -М.: Наука, 1987.-352 с.

13. Писаренко, Г.С. Справочник по сопротивлению материалов / Г.С. Писаренко, А.П. Яковлев, В.В. Матвеев // Киев: Наукова думка, 1989. - 732 с.

14. Тимошенко, С.П. Устойчивость стержней, пластин и оболочек / С.П. Тимошенко // М.: Наука, 1974. - 808 с.

15. Битюрин, А.А. Продольный удар неоднородного стержня о жесткую преграду / А.А. Битюрин, В.К. Манжосов // Ульяновск: Изд-во УлГТУ, 2009. - 164 с.

Поступила в редакцию 21 февраля 2010 г.

STABILITY LOSS OF THE HOMOGENEOUS CORE AT LONGITUDINAL IMPACT ON THE CORE, INTERACTING WITH THE HARD BARRIER

Mathematical modeling of longitudinal elastic central impact of a homogeneous core of a smaller cross-section on the homogeneous core integrating with a hard barrier with unilateral constraints in contact and impact sections is carried out. The rate of critical preimpact speed leading to loss of stability in the considered core system is simulated.

Keywords: longitudinal impact, stability, deformation, mathematical modeling, speed, hard barrier, core.

Bityurin Anatoly Aleksandrovich is a doctoral candidate, Ulyanovsk State Technical University.

Битюрин Анатолий Александрович - докторант, Ульяновский государственный технический университет.

e-mail: denjgy0706@yandex.ru