Динамическая реакция пластины на действие движущейся нагрузки

УДК 539.3

Динамическая реакция пластины на действие движущейся нагрузки

Александр Н.Блинов

Институт математики, Сибирский федеральный университет, Свободный 79, Красноярск, 660041,

Россия

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

Ключевые слова: установившееся движение, упругое основание, резонанс.

Введение

Знание реакции упругого слоя на действие движущейся нагрузки имеет большое практическое значение [1]. Например, при некоторой скорости нагрузки, называемой критической, возникает явление резонанса, при котором прогиб упругого слоя возрастает, что может привести к разрушению пластины.

В настоящей работе рассматривается упругий слой на упругом основании, по которому движется с постоянной скоростью нагрузка. В первой части работы изучается одномерная постановка задачи: бесконечная упругая балка на винклеровом основании под действием постоянной нагрузки, движущейся равномерно, т.е. с постоянной скоростью. Получено точное значение прогиба балки и найдена критическая скорость. Во второй части работы исследуется двумерная постановка задачи. В качестве модели для упругого слоя взята система уравнений Миндлина [2], а для описания упругого основания — уравнения теории упругости в случае плоского деформированного состояния. Условия сопряжения между ними — скользящий контакт. Получено решение для сосредоточенной нагрузки, как функции x — vt. Условие возникновения резонанса сводится к нелинейному алгебраическому уравнению. Наименьший положительный корень уравнения и является критической скоростью.

1. Реакция балки на винклеровом основании на действие движущейся нагрузки

Груз весом P равномерно движется со скоростью v вдоль бесконечной балки, лежащей на сплошном однородном основании. Картина изгиба оси балки остается стационарной равномерно движущейся со скоростью движения груза. Такое явление называют бегущей изгиб-ной волной. Вертикальная координата груза остается неизменной, и поэтому вертикальное ускорение груза равно нулю, а давление груза на балку равно весу P. Основание балки

© Siberian Federal University. All rights reserved

будем считать линейно деформируемым и следующим гипотезе Винклера: Q = — ку, где Q — интенсивность реакции основания, у — прогиб; к — коэффициент пропорциональности, характеризующий жестокость основания, называемый коэффициентом постели. Будем считать, что упругое основание лишено свойств инерции и демпфирования. Таким образом, упругость — единственное физическое свойство, которым обладает данная модель упругого основания. Пусть х — абсцисса текущего сечения балки, отсчитываемая от некоторого неподвижного начала координат, £ — время, Е.1 — жесткость поперечного сечения балки при изгибе. Дифференциальное уравнение изгиба балки записывается в форме

д4у д2у

ЕЛ3Х4 = — ку'

Здесь правая часть представляет собой интенсивность нагрузки в точке х в момент времени £ и состоит из двух членов: инерционной нагрузки (т — масса единицы длины балки) и реакции упругого основания.

Запишем уравнение (1) в форме

0 + 241 + Ь2у = 0,

где

a =

m ,2 k

2EJ 62 = EJ.

Решение этого уравнения ищем в виде y = f (x — vt). Здесь аргумент x — vt представляет собой абсциссу текущего сечения балки, отсчитываемую от подвижного начала координат, который совмещен с грузом.

После подстановки уравнение в частных производных (2) переходит в следующее обыкновенное дифференциальное уравнение:

fIV + 2av2f'' + b2f = 0.

Уравнение (4) представляет собой дифференциальное уравнение бегущей волны. Решая его, получим:

f = e-a«

где С = x — vt, а величины a и в определяются через коэффициенты уравнения (4) по формулам:

a = в = УЦО^ .

Будем считать, что решение (5) относится к той части бегущей волны, которая расположена впереди движущегося груза (т.е. при С > 0). Для другой части волны (при С < 0) уравнение изогнутой оси имеет ту же форму:

f1 = e-a«

но с другими значениями констант.

Найдем восемь постоянных C¿ и D¿ из (5) и (7). Для этого следует воспользоваться условиями:

{на бесконечности f = 0 при С = <х, fi =0 при С =

под грузом f (0) = fi (0), f' (0) = f1 (0), f'' (0) = fl' (0),

Л"

Из (8) находим:

В1 = В2 = Сз = С4 = 0.

Р Р

С1 = — Dз =--„ ^ /—::;-—г— , Со = — Dл =--

2Е7в (а2 + в2)' 2 4 2Е7а (а2 + в2) '

Таким образом, бегущая волна имеет совершенно симметричную форму, и для анализа полученных результатов достаточно рассмотреть, например, часть волны при £ ^ 0:

1 = - 2Е7а

Наибольший интерес представляет значение прогиба под грузом, т.е. при £ = 0 :

Р

1 (0) = -2ЕЗа. (а2 + в2) . Подставим сюда а, в из (6), получим

Р

I (0) = -

— ау2 )

При у = 0, т.е. в случае неподвижной нагрузки, получим известный результат для прогиба (обозначаемого 1с) под действием силы, приложенной к бесконечно длиной балке на упругом основании. При скорости у > 0 прогиб под грузом будет превосходить величину 1с. Отношение м = I/1с называется коэффициентом динамичности, и он равен

2у/кЕ.

Из этой формулы видно влияние всех параметров системы. С ростом скорости у динамический коэффициент увеличивается и при скорости

2л/кШ укр = \ -

V т

обращается в бесконечность. Найденное значение скорости является критическим.

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

2. Реакция пластины на упругом основании на действие движущейся нагрузки

Упругий слой толщины 2к лежит на поверхности полубесконечной упругой среды. Слой и полупространство однородные, изотропные и линейно-упругие. По слою движется нормальная нагрузка Е (х, £) с постоянной скоростью у. Рассматривается плоское деформированное

2

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

Перемещения и напряжения в полупространстве можно записать в следующем виде:

ду дф ду дф

дх ду ' ду дх '

Тх _ 2(1 — 0) дV

М

+

20

д2у _ 2_д2ф

1 — 20 дх2 1 — 20 ду2 дхду

(Ту _ 2(1 — 0) д2у 20 д2 у 2д2ф

М 1 — 20 ду2 1 — 20 дх2 дхду '

ху М

д2

у д2ф д2ф

дхду дх2 ду2 Волновые потенциалы у (х, у,Ь), ф (х, у, Ь) удовлетворяют уравнениям

1 дV д2ф + д2ф _ 1 д2ф

д2у д2у

дх2 + V ~ с2 дЬ2 ' дх2 ' ду2 ~ С2 дЬ2 '

с2 = 2м (1 — 0) с2 = М С1= Р (1 — 20) , С2 = ^

где р — плотность; м — модуль сдвига; 0 — коэффициент Пуассона; С; — скорость продольных волн; Ст — скорость поперечных волн.

Используя подстановку г = х — юЬ, получим, что у, ф удовлетворяют уравнениям

2 &21 + = 0 ^2 + = 0 дг2 ду2 ' дг2 ду2 '

1 — в2

Я

1

1 — 20 в2 в2 = -2(1 — 0) ^ в С2 '

Ограничимся рассмотрением скоростей движения нагрузки, меньших скорости распространения поперечных волн в полупространстве, т.е. в < 1. В этом случае уравнения будут эллиптическими. Для пластины под действием поперечных сил Миндлин [2] вывел систему приближенных уравнений, в которых учитывается влияние инерции вращения и поперечного сдвига. С учетом подстановки г = х — Ы она имеет следующий вид:

2 к 2\ (2и 0

Т—07 — 2р1Ну ) (2 — т = 0,

{2хМ\Ь — 2р\Нги1

(2Ш , (Н

-ГТ — 2ХМ1к^— Р = Л аг2 аг

N -кт = 0,

4 М1к3 2 ,о Л ¿2Н , ((Ш ,

зГ—07 — зР^ ) (Н —

где и, Ш — перемещения срединной поверхности слоя; Н — угол поворота поперечного сечения; х = 0,845 — сдвиговой коэффициент Тимошенко; р, т — нормальное и касательное напряжения, действующие на границе раздела слоя и полупространства; Е — внешняя нагрузка на верхней поверхности слоя.

2

в

Предполагается, что решение затухает в глубине полупространства, поэтому граничные условия при y = —ж будут иметь вид

* = 0, ф = 0. (11) На верхней поверхности слоя при y = h

p = —F, т = 0.

На границе раздела слоя и полупространства рассматривается условие сопряжения в виде скользящего контакта при y = — h:

p = ay, W = v, т = 0, axy = 0.

Внешнюю нагрузку возьмем в виде F = Foeikr, где k — волновое число, eikr = cos kr+i sin kr. Решение будем искать в следующем виде [4]:

*=W*eikr, ф=Wф1кг,

где *, ф — функции только от y,

U =-^т Ueikr, W = Weikr, k k

т = те1кг, р = ре1кт, n = Кеш',

где и, Ш, т, р, n — константы.

Подставляя (13) в уравнения для потенциалов, получим

др - = 0, дУ2 - к2°2ф = 0'

Решение этих уравнений имеет вид

р = С\ек,1

где С\, С2, С3, С4 — постоянные.

Учитывая краевые условия (11), находим

С3 = С4 = 0.

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

* =_ciek9

-1- ciekq(y+h)+ikr, ф = —

Из краевых условий (12) на границе раздела при y = — h получим:

2iqC1 — C2 — C2s2 = 0,

2C2is — ^ Ci + Ciq2 = 1,

Ciq + C2i = W Ц, т = 0.

Подставив (13) в уравнения Миндлина (10), получим

2xhiW = К + §3-31 - 2 M_V - (2XMih - 2Plhv2) AW - - _ = F-, U = 0.

Из уравнений (18) и (19) находим константы

^ F^ 2iq Fo — ^lFo

C = —г, C2 =-^ —-, W =— I q

1 A , 2 1 + s2 A , A 1 q

м 2 h A A 1 _ mF- c

2q

1 + s2 0, U = 0,

где

A

2q

1 + s2

- qj B - C, B =

4x2h2 D :

1 + s2 1 - 2$

D = + -

1 - 2$ k

i ' Q1 Q oPlh§v2" k 3 1 - $1 3 Ml

2xh 4 kh3 2 , § 2 .. 1 — pih3v2 —.

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

1

F = — coskrdk = б (r), п J

где б (r) — дельта-функция Дирака.

Решение для прогиба пластины W, который нас и интересует, имеет вид

1 Г 1_

W =- -W cos krdk,

nMi J k

где

Обозначим

Тогда

w (k) = F0 (' -

L = - q), M = 2

E

2h2

pih2v2

A=

3(1 - $i) 3mi

k3LME - k2CE + k (LMx - x2hL) - Cx X + k2E

- (x + k2E) MiLF-

1W (k) = _

k V ; (k3MLE - k2CE + k (M - xh) xL - Cx) k

и

где

Условие возникновения резонанса соответствует кратному вещественному корню при разложении на элементарные множители функции (22). Из формул Кардано следует, что это условие эквивалентно равенству нулю детерминанта кубического уравнения, стоящего в знаменателе этой функции:

3 )3+

= _1 (_С_У , (М - хЬ) X

Р 3 V Мь) + МЕ ' = А (У + 1 (М - хЬ) X Сх

Ц 27 \ыь) +3 ыь МЕ МЬЕ .

Поскольку в уравнении (23) р и ц — являются функциями от V — скорости движения нагрузки, то оно и определяет величину этой скорости. Критической скоростью и будет наименьший положительный корень уравнения (23).

Список литературы

[1] Д.Ахенбах, Движущаяся нагрузка, приложенная к пластине на упругом полупространстве, Прикладная механика, (1967), №4, 83-88.

[2] J.Midlin, Influence of rotatory inertia and shear on flexural matrons of isotropic elastic plate, J. Appl. Mechanics, Trans. ASME, 73(1951), 31-38.

[3] Я.Г.Пановко, И.И.Губанова, Устойчивость и колебания упругих систем, М., Наука, 1964.

[4] А.И.Лурье, Теория упругости, М., Наука, 1970.

[5] В.Новацкий, Теория упругости, М., Мир, 1975.

Dynamic Response of Plates on the Effect of a Moving Load

Alexander N.Blinov

The problem of stationary motion of the load on the surface of plate lying on an elastic foundation is considered. From a practical point of view of a large interest is the speed at which a phenomenon of resonance may be observed. The solution of this problem described in our paper can be used as basis for practical experiments.

Keywords: elastic plate, moving load, resonance.