Решение задачи ударного взаимодействия упругого тела и пластинки Уфлянда-Миндлина с помощью лучевого метода Текст научной статьи по специальности «Физика»

УДК 539.3

А. А. Локтев, Д. А. Локтев

РЕШЕНИЕ ЗАДАЧИ УДАРНОГО ВЗАИМОДЕЙСТВИЯ УПРУГОГО ТЕЛА И ПЛАСТИНКИ УФЛЯНДА-МИНДЛИНА С ПОМОЩЬЮ ЛУЧЕВОГО МЕТОДА

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

E-mail: prtlokt@yandex.ru

Ключевые слова: динамический контакт, пластинка Уфлянда-Миндлина, лучевой метод.

Задачи ударного взаимодействия твердых тел и пластинок рассматривались неоднократно [1-9]. Поскольку любая задача динамического контакта состоит из задачи местного смятия и задачи деформирования взаимодействующих тел вне области контакта, то существует несколько подходов к решению задач такого класса. В работах [1-4] авторы пошли по пути усложнения моделей местного смятия в зоне контакта. В [3-7] применялся подход, основанный на использовании волновых уравнений Уфлянда-Миндлина, описывающих динамическое поведение пластинки, которые позволяют предполагать, что деформация материала мишени вне области контакта происходит за счет распространения с конечной скоростью нестационарной волны поперечного сдвига. Для описания основных неизвестных динамического контакта в [6] использовалось их разложение в ряд по функциям Бесселя. Степенные ряды по поверхностной координате и времени для искомых величин использованы в работах [3-5, 7]. В [4, 5] использованы одночленные лучевые разложения, которыми описывались динамические характеристики до окончания процесса взаимодействия. В работе [7] волновой подход обобщен на мишень в виде термоупругой ортотропной пластинки, обладающей криволинейной анизотропией. Усложненная модель мишени в виде многослойной пластинки из композиционных материалов рассматривалась в [8], где исследовано влияние толщин различных слоев на конечные характеристики динамического контакта.

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

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

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

Движение ударника и движение контактной области пластинки после начала взаимодействия описываются уравнениями

т (а + гю) = — Р (г), ркпг^п) = 2пт0Цг |Г=Г0 + Р (г), (1)

где а — сжатие упругого элемента при ударе; г — прогиб пластинки; р — плотность материала пластики; Р(г) — контактная сила; к — толщина пластинки; Qr — перерезывающая погонная сила на границе контактной области.

Контактная сила для упругого ударника имеет вид [9]

Р (г) = Ех (а — г), (2)

где Ег — модуль упругости буфера.

Метод решения. Вне контактной области за фронтом волновой поверхности искомая функция представляется в виде лучевого ряда

Рис. 1. Схема динамического контакта тела и пластинки посредством буфера

по координате и времени [3-5]:

Z(r,t) = ± i [Z,w ] t=,./G (t - ^) ' H (t - , (3)

k=0

^v/aj-ki

где ] = [дк2/дЬк] — скачки производных к-го порядка по времени £ от искомой функции Z на волновой поверхности £, при t = (г — г0)/С; г0 — начальный радиус; Н(£) — единичная функция Хевисайда; С — нормальная скорость распространения фронта волны.

Уравнения, описывающие динамическое поведение упругой изотропной пластинки в полярной системе координат с учетом инерции вращения нормали к серединной плоскости и деформации поперечного сдвига, имеют вид [6]

V — 40 + ОМ + о = Р23Д; Оо + - = р'Ш (4)

г от 12 от г

• /ад. + д.\ . /д. + ö/3r

Mr = D — + а— ; Mv = D--+ а—~

yd- r / у r dr

(dW • )

Qr = Kuh I — - д I,

(5)

где г и — полярные радиус и угол; Мг и М^ — изгибающие моменты; Яг — перерезывающая сила; Д — угловая скорость вращения нормали к срединной поверхности пластинки в направлении г; Ш = гЬ — скорость прогиба; К = п2/12; ц — модуль сдвига; а — коэффициент Пуассона; Б = Е (1 — а2 £ — время; точка над величинами

означает производную по времени.

Для определения переменных коэффициентов ] необходимо записать уравнения движения точек пластинки (4), (5) в скачках и использовать условие совместности для перехода от производных по поверхностной координате к производным по времени на волновой поверхности, после чего определяются скорость сдвиговой волны и динамическое условие совместности для перерезывающей силы и скорости прогиба [5]:

С = /Кц/р, Яг = —рС'Ш. (6)

Подставляя в уравнения (1) выражения (2) и (6), получаем систему интегро-дифференциальных уравнений относительно а и ь:

m (ä + W^j = -Ei (а - w);

, 2- 2n-oKphW (7)

phnr0W =---f--+ Ei (а - w).

G

Используя преобразование Лапласа, запишем систему (7) в виде

mp2 (¿5 + w) = —Ei (¿5 — w) + mVo;

2nroKuhp _ (8)

pknrffiw =--—-w + El (a — w).

G

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

Решая эту систему, получаем выражения

Р2 + - V

w = —-—m + ; (9)

2 Ei 2 Ei p2 — p2 —

mm

_ (p2 + B0p + A0) v0

« = -7-^^, (10)

p p3 + Bo p2 + Cop + — Bo

m

El 2G(2) / 2 1

где x = 1/ti; Ao = ——2; Bo =-; Co = Ei ——2 +--

phnr 0 ro \phnr2 m/

Для обратного перехода в пространство оригиналов необходимо разложить выражения (9) и (10) на простые дроби. Для этого нужно решить характеристическое уравнение

E

p3 + Bo p2 + Cop + E Bo = 0, (11)

m

которое может иметь два комплексно сопряженных и один действительный корень или три действительных корня.

Если уравнение (11) имеет действительные корни aii, ai2, ai3, то выражение для а имеет вид

- Ai + Bi + Ci + Di (12)

a =--1---1---1--, (12)

p— aii p — ai2 p— ai3 p

где

Vo (a?i + Бфи + Ao) Vo (a?? + Бф^ + Ao)

Ai = 7-r-:-г-; Б1 =

(aii - ai?) (aii - ai3) an' (ai? - an) (ai2 - ai3) ai2'

Vo (a?3 + Boai3 + Ao) AoVo Ci = 7-w-ч-; Di =

(а1э — ац) (а1э — а12) а1з' а11а12а1з

Выражение для а в пространстве оригиналов запишем как

а = А1 ехр (а11г) + В1 ехр (а12г) + С1 ехр (а13г) + Б1. (13) Если уравнение (11) имеет один действительный корень а16 и два

комплексно сопряженных ai4, ai5, то

_ A2P + B2 C2 , D2 а = 2 . /—^ +-+ —, (14)

p2 + a' p + b' p - ai6 p

где

2

„ _ Vo (a24 + Boai4 + Ao) Vo K5 + Boai5 + Ao)

A2 =

(а14 — а15) (ам — а1б) а14 (а14 — а15) (а15 — а1б) а15

В = Уо (а^4 + Во а 14 + Ар) а ^ . с = Уо (а^ + Во а^ + Ар) (а14 — а1б) а14 ' (а26 + а'а^ + Ь') аи

= — тт~°. а' = (Во + а1б); Ь' = [Со + а^а']. Ь'а1б

Выражение для а в пространстве оригиналов примет вид а = ехр(— 1/2а'*) <^'-1/2 вт £1/2^ [2В2 — а'А2] +

+ А2 } + С ехр (а1б£) + £2, (15)

где £' = 4Ь' — а'2.

При действительных корнях (11) функция г!) представляется как

Аз Вз Сз Бз 2 тУо Ь =-+-+-+ — + —2-, (16)

р — а11 р — а12 р — а13 р тр2 — Е1

где

A = Vo (mali + Ei) (a2i + Ao + Boan) ; (aii - ai2) (аи - aiß) an (ma^ - Ei)'

Vo (ma22 + Ei) (a^ + Ao + Boai2)

B = ^o (mai2 + Ei) (a2_

(ai2 - aii) (ai2 - a^) ai2 (ma^ - Ei)'

Vo (ma2s + Ei) + Ao + Boa^) n AoVo

C3 = --—-;-;-Ö-D3 =

(а1з — аи) (а1з — а12) а1з (та^ — Е1)' апа12а1з'

Переходя к пространству оригиналов, получаем

г = Аз ехр (ац£) + Вз ехр (а^£) + Сз ехр (а1з£) +

+ Бз + 2тГ—/ Уо вЬ| У*|. (17)

тт

В случае двух комплексно сопряженных корней характеристического уравнения выражение (9) примет вид

А4р + Вл С4 Бл 2тУо ,10Ч

Ь = т4-п; +-+ — + —2—, (18)

р2 + ар + Ь р — а16 р тр2 — Е1

где

Vo (ma^ + Ei) (a?4 + Aq + Еоам)

A4 =

(ai4 - ai5) (ai4 - ais) ai4 (maf4 - Ei)

vq (ma^ + Ei) (a?5 + aQ + Boa^)

(ai5 - aw) (ai5 - a^) aVj (ma^ - Ei)' B = Vo (maÍ4 + Ei) (a24 + Aq + Воам) ;

(ai4 - ai6) a 14 (maf4 - Ei)

C = Vo (mai6 + Ei) (a^6 + Aq + Bpa^); D = AiVq 4 (af6 + a'ai6 + b') ai6 (maf6 - Ei) ' 4 b'ai6

Совершая обратное преобразование Лапласа, находим

w = exp (-a't/2) x

X {^'-i/2 [2B4 - a'A4] + A4 cos Q e'i/2t) } +

+ C4 exp (ai6t) + D4 + 2mC0-i/2Vo sh (coi/2t) . (19)

После подстановки величин w и а в (2) можно получить выражение для контактной силы P(t). В случае действительных корней уравнения (11) она может быть представлена в виде

P (t) = Ei [(Ai - A3) exp (aiit) +

+ (Bi - B3) exp (ai2t) + (Ci - C3) exp (ai3t) + + (Di - D3) - 2m (Ei/m)-i/2 Vo sh ((Ei/m)i/2 t) . (20) При наличии комплексно сопряженных корней (11) получим

P (t) =

=Ei

exp (-a't/2) {e'-i/2 sinQ^t) í2 (B-B4) +a' (A2-A4)] -

— (А4 — А2) ] + (С2 — С4) exp (а1бг) + (А, — А4) —

—2т (Е1/т)-1/2 V эЬ ((Е1/т)1/2 г) . (21)

В работе [5] задача об ударе шаром по пластинке была решена без учета инерции контактной области (р = 0), что связано со значительными трудностями при решении полного интегро-дифференциального

уравнения, полученного из системы (1). Для определения вклада инерции контактного пятна в конечные динамические характеристики рассмотрим частный случай решенной задачи при р = 0.

Характеристическое уравнение (11) преобразуется в следующее:

p2 + 2C0 p + Eim-1 = 0, (22)

где CO =

Ei

MB0

Корни этого уравнения

«19,20 = -C0 (C0)2 - Ei/m. (23)

Рис. 2. Зависимость контактной силы от времени для разных подходов к описанию процесса взаимодействия

Выражения для динамического прогиба, контактной силы получаются из выражений (12)-(21) при р = 0.

Численные исследования. Исследуем зависимость контактной силы от свойств ударника и мишени. Параметры взаимодействия приняты следующими: т = 10 г, к = 13 мм, У0 = 10 м/с, Е\ = 250 кН/м. С использованием соотношений (20) и (21) на рис. 2 построены зависимости контактной силы от времени: кривая 1 построена по результатам работы [9] без учета волновых процессов в плите; кривая 2 соответствует подходу, примененному в данной работе; кривая 3 получена по результатам работы [7], где неизвестные величины были представлены в виде многочленных лучевых разложений и при этом использовалось дополнительное граничное условие на границе области контакта; кривая 4 получена без учета инерции контактной области при прогибе. Из рис. 2 следует, что учет волновых явлений приводит к снижению максимального значения контактной силы по сравнению с силой, полученной при использовании неволнового уравнения колебаний плиты; без учета инерции контактной области контактная сила уменьшается на 8 %, что может считаться существенной величиной при решении практических задач. На рис. 3 приведены зависимости контактной силы от времени для разных материалов мишени; видно уменьшение максимальной контактной силы и увеличение времени взаимодействия для материала с меньшим модулем сдвига. На рис. 4 показана зависимость максимальной контактной силы от толщины стальной и алюминиевой пластинок и начальной скорости удара. Из рис.4 следует, что значение максимальной контактной силы увеличивается с увеличением начальной скорости удара и толщины мишени. Можно также отметить, что при одинаковых значениях начальной скорости удара и толщины пластинки максимум

Рис. 3. Зависимость контактной силы от времени для различных материалов мишени:

1 — сталь; 2 — медь; 3 — алюминий

Рис. 4. Зависимость максимальной контактной силы от безразмерной толщины пластинки для различных значений начальной скорости удара, м^ (цифры у кривых): сплошные кривые соответствуют стальной пластине, штриховые — алюминиевой

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

Выводы. 1.Учет волновых процессов уменьшает максимальное значение контактной силы. 2. Необходимо учитывать инерционный член в уравнении движения контактной области в задачах, в которых погрешность не должна превышать 8%. З.При увеличении модуля сдвига контактная сила увеличивается, а время контакта уменьшается. 4. При увеличении толщины мишени и скорости удара максимальная контактная сила увеличивается.

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

1. Phillips J. W., C a l v i t H. H. Impact of a rigid sphere on a viscoelastic plate // Trans. of the ASME, Journal of Applied Mechanics. - 1967. - Vol. 34, № 4. -P. 873-878.

2. Александров В. М., Ромалис Б. Л. Контактные задачи в машиностроении - М.: Машиностроение, 1986. - 362 с.

3. Россихин Ю. А., Шитикова М. В., Локтев А. А. Удар шара о нелинейно упругий буфер, установленный на плите перекрытия // Изв. вузов. Строительство. - 2004. - № 11. - С. 16-22.

4. Локтев А. А. Упругопластическая модель взаимодействия цилиндрического ударника и пластинки // Письма в журнал технической физики. - 2007. - Т. 33, вып. 16.-С. 72-77.

5. RossikhinYu. A., ShitikovaM. V. A ray method of solving problems connected with a shock interaction // Acta Mechanica. - 1994. - Vol. 102, nos. 1-4. -P. 103-121.

6. Филиппов А. П. Поперечный упругий удар тяжелым телом по круглой плите // Механика твердого тела. - 1971. - № 6. - С. 102-109.

7. Л о к т е в А. А. Динамический контакт ударника и упругой ортотропной пластинки при наличии распространяющихся термоупругих волн // Прикладная математика и механика. - 2008. - Т. 72, вып. 4. - С. 652-658.

8. Chen P., X i o n g J., Shen Z. Thickness effect on the contact behavior of a composite laminate indented by a rigid sphere // Mechanics of Materials. - 2008. -V. 40. - P. 183-194.

9. C o n w a y H. D., L e e H. C. Impact of an indenter on a large plate // Trans. of the ASME. Journal of Applied Mechanics. - 1970. - Vol. 37, no. 1. - P. 234-235.

Статья поступила в редакцию 18.11.2009

Алексей Алексеевич Локтев окончил Воронежский государственный архитектурно-строительный университет в 2001 г. Канд. физ.-мат. наук, зав. кафедрой Московской финансово-юридической академии. Автор около 60 научных работ в области механики контактного взаимодействия, теории упругости, волновых процессов в твердых телах, математического моделирования в механике и информационных технологий.

A.A. Loktev graduated from the Voronezh State University for Architecture and Civil Engineering in 2001. Ph. D. (Phys.-Math.), head of department of the Moscow Academy of Finance and Law. Author of about 60 publications in the field of mechanics of contact interaction, theory of elasticity, wave processes in solid bodies, mathematical simulation in mechanics and information technologies.

Даниил Алексеевич Локтев — студент МГТУ им. Н.Э. Баумана. Автор семи научных работ в области механики контактного взаимодействия, теория упругости, волновых процессов в телах, вычислительной техники и информатики.

D.A. Loktev — student of the Bauman Moscow State Technical University. Author of 7 publications in the field of mechanics of contact interaction, theory of elasticity, wave processes in bodies, computing hardware and information technologies.