Модификация дифференциальных уравнений волновой теории продольного соударения стержней

Математика и механика. Физика

УДК 621.01

МОДИФИКАЦИЯ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ ВОЛНОВОЙ ТЕОРИИ ПРОДОЛЬНОГО СОУДАРЕНИЯ СТЕРЖНЕЙ

И.А. Жуков, Л.Т. Дворников

Сибирский государственный индустриальный университет, г. Новокузнецк E-mail: zhival@yandex.ru

Показана модификация одномерной волновой теории Сен-Венана при продольном соударении стержней сложных геометрических форм с учетом деформации в радиальном направлении.

Ключевые слова:

Удар, колебания, деформация, импульс, боек, волновод.

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

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

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

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

• классический ньютоновский метод;

• метод Герца, предполагающий области контакта упругими, а тела твердыми;

• метод, предполагающий тела полностью упругими, но распространение напряжений по телам мгновенным;

• метод одномерной волны Сен-Венана;

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

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

Наиболее широкое применение при исследовании ударных систем получила одномерная волновая теория удара, сформулированная Барре де Сен-Венаном [2], нашедшего общее решение задачи в форме, допускающей практические приложения. Теория Сен-Венана построена для тонких стержней с плоскими торцами на тех допущениях, что: 1) плоские, поперечные к оси стержня, сечения остаются плоскими в процессе распространения волн продольной деформации; 2) материал стержня подчиняется закону Гука, т. е. деформации остаются в пределах упругости; 3) соприкосновение соударяющихся тел происходит в один и тот же момент времени по всей площади ударного торца.

По теории Сен-Венана процесс продольного соударения стержней полностью определяется функцией смещения ы(х,1) поперечного сечения стержня с координатой х в момент времени ¡. Ось координат поперечных сечений х совпадает с прямой, проходящей через центры сечений. Процесс распространения волн продольных колебаний в стержне описывается дифференциальным уравнением гиперболического типа

д и

X Х+сЫ

Рис. 1. Процесс распространения волн продольных колебаний в стержне переменного поперечного сечения

Выделим элемент стержня ограниченный боковой поверхностью и двумя поперечными сечениями, отстоящими на бесконечно малом расстоянии йх друг от друга, с координатами х и х+йх. Выделенный элемент будет иметь объем йУ = Б (х) йх,

и массу

йт = рйУ = р Б(х)йх.

Ускорение смещения сечения а(х,1) определятся соответственно второй производной по времени от функции смещения сечений

*(х, о =ддг1. д/

Относительная деформация е

е( x, t) =-

д x

ст( x, t) = E -е( x, t) = E

ди

F

д x ди

д x

Второй закон Ньютона для выделенного элемента запишется в виде

dm  a(x, t) = F(x + dx, t) - F(x, t),

д2u

P^ S (x)-

P^ S

дt2 dx

д2 и

дt2

дx

д/2 дх2

где а=лЕ/р - скорость распространения упругих волн продольной деформации в стержне с плотностью материала р и модулем упругости Е.

Как видно, ур. (1) учитывает только свойства материалов соударяющихся тел и никак не учитывает их геометрию. Рассмотрим процесс распространения волн продольных колебаний в стержне переменного поперечного сечения (рис. 1).

Б(х) х+сЫ

Подставляя выражение усилия (2) в (3), получаем модификацию уравнения продольных колебаний стержней с учетом их геометрии

чд2и

Ур. (4) для стержней переменного поперечного сечения показано в работе Я.Г. Пановко [3].

Выполняя операции дифференцирования, можно представить ур. (4) в другой форме

2 д2и(х,/) 2 1 йБ(х) ди(х,/) д2и(х,/) „ ,гч

а-+ а--— —--^^

дх Б (х) йх дх д/

Для стержня постоянного поперечного сечения ур. (5) преобразуется в уравнение (1).

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

Начальные условия следующие:

- в момент начала взаимодействия смещения сечений стержня м>(х^) и бойка ы(х,1) равны нулю:

w( х,0) = 0, и (х,0) = 0;

- скорость смещения, определяемая частной производной по времени, для бойка равна его предударной скорости У0, а для стержня равна нулю:

3w( x,0) 3t

= V0,

ди( x,0) 3t

= 0.

Граничные условия, определяющие состояние концов бойка и стержня, следующие:

- в процессе взаимодействия смещения на границе бойка и стержня равны:

^(0, г) = и (0,г);

- в процессе взаимодействия силы взаимодействия на границе бойка и стержня равны:

?(0) ^(0,т) = „ ди(0,т);

^(0) Г - ^0 : ;

дх дх

- неударный торец бойка свободен от деформаций:

д^

дх

= 0;

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

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

lim ^UlXA = о.

х дх

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

Однако экспериментальные опыты свидетельствуют о некорректности допущений, принятых в одномерной волновой теории удара, и невозможности их осуществления на практике. В ударных системах сложной геометрической формы предположение, что в процессе деформирования плоские поперечные сечения остаются плоскими, весьма сомнительно, т. к. совместно с продольными колебаниями возникают и поперечные. Эти факты свидетельствуют о необходимости рассмотрения деформации стержней в радиальном направлении. Впервые на необходимость более детального рассмотрения задачи о продольном ударе указал Т. Юнг [4], введя ограничение на предударную скорость бойка и принимая в учет влияние его массы. В 1876 г. Л. Похгаммером была опубликована работа [5], в которой он рассматривает задачу о продольных колебаниях цилиндрического стержня кругового сечения, учитывая при этом поперечные перемещения.

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

Вывод модифицированного уравнения продольных колебаний в цилиндрической системе координат показан в работе [6]. Однако использование цилиндрической системы координат накладывает ограничения на возможность расчета ударных

Рис. 2. Колебания в стержне при ударе

Напряженно-деформированное состояние и кинематика стержня при ударе (рис. 2) определяются функцией смещения и(х,{) сечения вдоль оси, как и по теории Сен-Венана. Относительная продольная деформация элемента определяется выражением

ди

ех =-.

дх

Согласно закону Гука, при продольной деформации тела возникают относительные деформации ег в поперечном направлении

ди

ег = -/ех = -/л-= —--,

дх д-

где /л - коэффициент Пуассона; иг(х,г,$ - функция радиального перемещения.

Тогда функция радиального перемещения полностью определяется функцией осевого перемещения сечения и расстоянием г от точки сечения до оси х

ди

иг (х, r, t) = -¡r -

дх

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

б| Ldt = 0,

где Ь=(Т—П) - функция Лагранжа, определяемая разностью кинетической Т и потенциальной П энергий системы.

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

Таким образом, функция Лагранжа определяется выражением

ь = тх + Г, - .

Тогда

I

(T + T -П)dt = °

Кинетическая энергия выделенного элемента с функцией площади поперечного сечения £(х) в осевом направлении определяется выражением

г.=и 2 ^ ]'=2 м ^дт11]2 ^

1 х

Аналогично получается выражение кинетической энергии в радиальном направлении

\ 2

диг

T =

И 1

l S

dm

д t

2J/2p ( J г2dS) (^дгÏ dX'

2, si) I дхдt J

Интеграл

J r 2dS =J p

s

представляет собой полярный момент инерции сечения.

Следовательно

T = 2 Jp/2 J р

( д2и ( x, t ) ^2 v dxdt у

dX'

Для определения потенциальной энергии выделенного элемента учитывается, что при продольном деформировании возникают только нормальные напряжения стх(х,0=ст(х,0, параллельные продольной оси, которые в пределах упругости определяются законом Гука

ди

ст(x, t) = Esx(x, t) = E -

дх

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

=

J J 2ст

дх

С учетом полученных соотношений (7)-(9) определяется выражение (6) принципа Гамильто-на-Остроградского для продольных колебаний выделенного элемента с учетом поправки Релея

2

1 Л ди

äJJ

2 ps

ï

dt

+тРЛ jp

д2 и ( х, t) dxdt

-1 es (x) дМ

2 I дх

dxdt = О.

Задача вывода уравнения продольных колебаний сводится к исследованию функционала действия по Гамильтону-Остроградскому на экстремум.

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

д

дх

д у

дt I ду J дхдt

кдуг у

=

где ег =

ди

v = -

vr=

д 2и

дх dt dхdt

Вычисляя производные, входящие в условие экстремальности (10), получаем

д дх

ES

ди ( х, t ) дх

±(pS ( х) ^^ V

dt Г dt J

dхdt

p/л2 J p

д2и ( х, t ) dхdt

= О.

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

£ дх [х (х) Ёдт (х) ^+

д

д3и

+P"Tx I JP(X) дхд(2

= 0'

Если пренебречь деформацией в радиальном направлении, приняв /=0, получается уравнение (4).

Приведенное модифицированное уравнение (11) свидетельствует о том, что поправка Релея приводит к появлению члена более высокого порядка малости. Поскольку поправка Релея зависит от коэффициента Пуассона /л и полярного момента инерции Jp(x), можно сделать вывод, что решение зависит от большего числа механических параметров и особенностей геометрии стержней, чем по одномерной теории Сен-Венана. Таким образом, предложенная авторами модификация может стать фундаментальной основой при решении проблем, связанных с продольным ударом.

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

1. Александров Е.В., Соколинский В.Б. Прикладная теория и расчеты ударных систем. - М.: Наука, 1969. - 201 с.

2. Clebs^ A. Theorie de l'elasticite des corps solides / V.f. Saint-Venant. - Paris: Dunod, 1883. - 980 p.

3. Пановоко Я.Г. Основы прикладной теории колебаний и удара. 7 Изд. 3-е., доп. и переработ. - Л.: Машиностроение, 1976. -

С. 133.

4. Thomas Y. A course of lectures on natural philosophy and the mechanical arts. - London, 1807.

5. Pochhammer L. Über die Fortpflanzungsgeschwindigkeiten kleiner Schwingungen in einem unbegrentzen isotropen Kreiszylinder // J. f. d. reine und angew. - 1876. - Math., 81. - S. 324.

Мясников А.А. Модифицированное уравнение продольных колебаний стержней переменного поперечного сечения в цилиндрической системе координат // Матер. VII научно-практ. конф. по проблемам машиностроения, металлургических и горных машин / Под ред. проф. Л.Т. Дворникова. - Новокузнецк: СибГГМА, 1998. - С. 70-79.

Пат. 2090753 РФ. МПК6 Е21С 5/02. Ударная бурильная машина / Л.Т. Дворников, Е.Ф. Губанов. - Приоритет от 08.06.1993, опубл. 20.09.97, Бюл. № 26.

Поступила 09.06.2008 г.

УДК 621.01

ЗАДАЧА СТРУКТУРНОГО СИНТЕЗА ПРОСТРАНСТВЕННЫХ АССУРОВЫХ ТРЕХЗВЕННЫХ МЕХАНИЗМОВ

Л.Т. Дворников, М.Г. Попугаев

Сибирский государственный индустриальный университет, г. Новокузнецк E-mail: fdba@yandex.ru

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

Ключевые слова:

Трехзвенные механизмы, однозвенные группы Ассура, кинематические пары.

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

В трехзвенных механизмах подвижными являются два звена, а третьим звеном считается стойка, относительно которой рассматривается движение. Так как ведущим звеном в них будет или кривошип или ползун, обладающие подвижностью Ж=1, то при создании механизма с той же подвижностью, что и у ведущего звена, т. е. с Ж=1, вторым звеном, должна являться система, обладающая нулевой подвижностью, так называемая группа Ассура. Если найти структуры всех возможных однозвенных групп Ассура, то тем самым будет решена задача о полном составе трехзвенных ассуровых механизмов.

Обратимся к поиску структур однозвенных пространственных групп Ассура. Известно [1], что все пространственные кинематические цепи описываются структурной формулой А.П. Малышева, которая записывается в виде уравнения

№ = 6п - 5р5 - 4рА - 3р3 - 2р2 - р..

В (1) обозначены: Ж - подвижность кинематической цепи (число ее обобщенных координат); п -число подвижных звеньев цепи; р1 - число используемых в цепи кинематических пар, где I - класс пары, соответственно: 1=5 - пятый, /=4 - четвертый и т. д.

Однозвенным пространственным группам Ассура соответствуют следующие параметры: Ж=0 и п=1, тогда формула (1) принимает вид

5ръ + 4р4 + 3р3 + 2 р 2 + р, = 6.

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

р = т + (т - 1)п/-1 +... + п +... + 2п2 + п,,

п = 1 + пт-1 +... + п1 +... + п 2 + п..

В

р = р5 + р4 + рз + р2 + рV