Механика столкновений твердых тел: влияние трения и адгезии. Ii. Численное моделирование Текст научной статьи по специальности «Физика»

Ляшенко Я.А., Виллерт Э., Попов В.Л. Механика столкновений твердых тел: влияние трения и адгезии. II. Численное моделирование // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2018. - № 4. -С. 178-199. DOI: 10.15593/perm.mech/2018.4.17

Lyashenko I.A., Willert E., Popov V.L. Mechanics of collisions of solids: influence of friction and adhesion. II. Numerical modeling. PNRPUMechanics Bulletin, 2018, no. 4, pp. 178-199. DOI: 10.15593/perm.mech/2018.4.17

ВЕСТНИК ПНИПУ. МЕХАНИКА

№ 4, 2018 PNRPU MECHANICS BULLETIN

http://vestnik.pstu.ru/mechanics/about/inf/

Б01: 10.15593/регш.шееЬ/2018.4.17 УДК: 531.133.1; 531.311; 539.312

МЕХАНИКА СТОЛКНОВЕНИЙ ТВЕРДЫХ ТЕЛ: ВЛИЯНИЕ ТРЕНИЯ И АДГЕЗИИ. II. ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ

Я.А. Ляшенко12, Э. Виллерт1, В.Л. Попов134

1Берлинский технический университет, Берлин, Германия 2Сумский государственный университет, Сумы, Украина

Национальный исследовательский Томский государственный университет, Томск, Россия Национальный исследовательский Томский политехнический университет, Томск, Россия

О СТАТЬЕ

АННОТАЦИЯ

Получена: 05 июня 2018 г. Принята: 19 июня 2018 г. Опубликована: 28 декабря 2018 г.

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

столкновение, нормальный и тангенциальный контакт, адгезия, трение, скольжение, коэффициент восстановления, метод редукции размерности.

Данная работа представляет собой вторую часть обзора, посвященного физике двухчастичных столкновений твердых тел. В первой части описаны теоретические и экспериментальные работы по неупругим столкновениям твёрдых тел, в которых диссипация энергии обусловлена наличием внутреннего или внешнего трения, пластичности, адгезии, а также других каналов затухания. В настоящей части обзора мы уделили основное внимание случаю столкновений упругих частиц. Приведены результаты математического и численного моделирования, полученные авторами с применением метода редукции размерности, который позволяет описание трехмерного контакта свести к моделированию контактных процессов в эквивалентном одномерном пространстве. Рассмотрено три случая: столкновение между телами без проскальзывания (эквивалент бесконечного коэффициента трения), с конечным коэффициентом трения в зоне контакта и с наличием адгезионного взаимодействия, которое описывается в пределе ЖР Во всех рассматриваемых случаях идентифицированы определяющие процесс безразмерные переменные, и численно определены функции, задающие зависимости между этими переменными. Используя указанные безразмерные переменные, возможно рассчитать компоненты скоростей и циклическую частоту вращения, а значит, и траекторию шарика после столкновения, если известны аналогичные параметры до столкновения. Теоретические результаты сравниваются с опубликованными экспериментальными данными и показывают хорошее совпадение с известными экспериментами.

©ПНИПУ

© Ляшенко Яков Александрович - дф-мн, проф., н.с., e-mail: [email protected], 0000-0001 -7511 -3163 Виллерт Эмануэль - н.с., e-mail: [email protected]), ¡D 0000-0001 -7535-7301 Попов Валентин Леонидович (Берлин, Германия) - д.ф.-м.н., проф., зав. каф, e-mail: [email protected], 0000-0003-0506-3804

lakov A. Lyashenko - Doctor of Physical and Mathematical Sciences, Professor, Scientific Researcher, e-mail: [email protected], 0000-0001-7511-3163

Emanuel Willert - Researcher, e-mail: [email protected], iD 0000-0001-7535-7301 Valentin L. Popov - Doctor of Physical and Mathematical Sciences, Professor, Head of Department. e-mail: [email protected], 0000-0003-0506-3804

Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)

This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)

MECHANICS OF COLLISIONS OF SOLIDS: INFLUENCE OF FRICTION AND ADHESION. II NUMERICAL MODELING

I.A. Lyashenko12, E. Willert1, V.L. Popov1'3'4

''Berlin University of Technology, Berlin, Germany 2Sumy State University, Sumy, Ukraine

3National Research Tomsk State University, Tomsk, Russian Federation 4National Research Tomsk Polytechnic University, Tomsk, Russian Federation

ARTICLE INFO

ABSTRACT

Received: 05 June 2018 Accepted: 19 June 2018 Published: 28 December 2018

Keywords:

impact, normal and tangential contact, adhesion, friction, partial sliding, restitution coefficient, method of dimensionality reduction

This paper is the second part of the review on the physics of two-particle collisions of solids. The first one describes theoretical and experimental works on collisions of elastic and elastic-plastic solid bodies in the case when the energy dissipation is caused by the inner or interface friction, plasticity, adhesion, or some other damping mechanisms. In this second part, we concentrate on the case of collisions of elastic particles. Results of the analytical and numerical modelling within the dimensionality reduction (MDR) method are presented. MDR allows describing a three-dimensional contact by a reduction to an equivalent problem in a one-dimensional space. We consider three cases: a collision between bodies without slipping (equivalent to an infinite coefficient of friction), with a finite coefficient of friction in the interface, and with the presence of adhesion, which is described in the JKR limit. In all the considered cases, proper dimensionless variables have been identified which fully characterize the collision process, and the functions determining the transformation of the state before the impact into the state after the impact has been determined numerically. By using these dimensionless variables, it is possible to calculate the velocity components and the angular velocity, and hence the trajectory of the ball after the collision, in case if analogous parameters are known before the collision. Theoretical results are compared with the available results, the experimental data show agree well with the experiments.

©PNRPU

Введение

Физические процессы, происходящие при столкновениях твердых частиц, представляют существенный интерес для различных технологических процессов, которые связаны с динамикой гранулированных сред [1-5]. Сложность теоретического описания процесса столкновения связана прежде всего с тем, что даже в простейшем случае упругого столкновения двух сферических частиц динамика их взаимодействия является очень сложной [6], поскольку может реализоваться множество различных случаев. Это частичное проскальзывание либо полное скольжение во время контакта, а также отсутствие скольжения на протяжении различных фаз удара [7]. На сегодняшний день не существует полной теории либо эффективных численных методов, которые позволяли бы описывать все процессы, происходящие при столкновениях. В предыдущей работе [8] нами был проведен обзор существующих теоретических и экспериментальных работ по изучению столкновений в случаях наличия трения в зоне контакта, адгезионного взаимодействия, а также пластической деформации. В этой части обзора мы проводим численное моделирование процесса удара в рамках метода редукции размерности [9], который позволяет описать многие экспериментальные результаты других авторов, освещенных в [8].

В своих недавних работах В.Л. Попов и другие показали, что теория Герца-Миндлина может быть точно вос-

произведена с помощью замены реального трехмерного контакта [10] контактом между модифицированным профилем с линейным упругим основанием, состоящим из независимых пружин [8, 11-14]. Этот метод можно использовать не только для точного описания нормального и тангенциального контакта при произвольной истории нагрузки, но и с хорошей точностью для контактов качения [15]. Указанный подход, получивший название метода редукции размерности (MDR), значительно упрощает проблему описания контакта, открывая при этом новые возможности для построения аналитических моделей и численного анализа динамически изменяющихся нормальных и тангенциальных контактов. Помимо этого, в рамках метода MDR возможно также описание контакта с адгезией, что позволяет провести моделирование столкновения тел, между которыми существует адгезионное взаимодействие. Такие ситуации нередко возникают в различных технологических процессах, к примеру при газодинамическом напылении [16], когда важно, чтобы после столкновения частица оставалась на поверхности. Настоящий обзор посвящен описанию упругого столкновения сферических тел под углом к упругому полупространству, которое основывается на результате численных экспериментов в рамках метода редукции размерности. Как мы уже указывали, эта часть обзора (первую часть см. в [8]) посвящена непосредственно моделированию в рамках метода MDR и обобщению полученных результатов, а также их сравнению с известными экспериментами.

Работа построена следующим образом. В подразд. 1 описана постановка задачи, приведена схема рассматриваемой системы и классическое аналитическое решение. Подразд. 2 посвящен обзору результатов, полученных авторами непосредственно с помощью моделирования в рамках метода редукции размерности. Приведены основные предположения и допущения, в рамках которых работает указанный подход. Отдельно рассмотрены три случая: столкновение без проскальзывания (эквивалент бесконечного коэффициента трения), столкновение с конечным коэффициентом трения и столкновение при наличии адгезионного взаимодействия между поверхностями контактирующих тел. Во всех случаях приведено аналитическое исследование линеаризованной модели с постоянной контактной жесткостью, в результате которого определены функции, полностью описывающие столкновение. С помощью моделирования в широком диапазоне параметров показано, что найденные функции полностью задают параметры столкновения при реальном контакте между сферическим телом и упругой полуплоскостью. Проведен численный расчет указанных функций, а в случаях без проскальзывания и наличия адгезионного взаимодействия в актуальной области параметров записаны аналитические аппроксимации. Показано, что использование этих аппроксимаций позволяет найти любую величину после столкновения: угловую скорость вращения, нормальную и тангенциальную компоненту скоростей, угол отражения, коэффициенты восстановления и т.д. Для конечного коэффициента трения и случая с адгезией проведено сравнение с экспериментами и показано, что теоретические результаты описывают экспериментальные с высокой точностью. Более того, в ситуации с конечным коэффициентом трения теория в том числе описывает область малых углов падения меньше 30°, что является ее несомненным преимуществом по сравнению с другими подходами.

1. Схематическое представление задачи и ее классическое аналитическое решение

В настоящей работе мы ограничиваемся рассмотрением столкновения частиц с полупространством. Эта базовая задача может быть легко обобщена на случай столкновения двух конечных частиц. Во избежание излишнего усложнения мы отказались от явного рассмотрения механизма сведения общей задачи столкновений к задаче столкновения с полупространством. Детали этой редукции можно найти в классических учебниках по теории столкновений, например в [17]. Рассмотрим столкновение упругого шара массой т и радиусом Я с упругим полупространством, как это показано на рис. 1.

Пусть модули упругости шара и полупространства Е и Е2, коэффициенты Пуассона ^ и у2 , сдвиговые модули ^ и 02 , а плотности материалов шара и полупространства р; и р2. Начальная скорость шара V до

акта столкновения раскладывается на вертикальную уг0 и горизонтальную ух0 компоненты, начальная циклическая частота вращения шара ю0. После акта столкновения часть кинетической энергии передается телу, с которым шар контактировал, и в результате эти величины становятся равными V, —V, V, ю. Зная величины компонент скоростей, легко определить угол скольжения для падающего тела а , а также угол скольжения (3 , с которым шар отлетает от поверхности.

Рис. 1. Схематическое представление столкновения сферического тела с полупространством под углом к поверхности Fig. 1. Schematic representation of the collision of a spherical body with a half-space at an angle to the surface

Существует классическое решение этой задачи. Обозначим через F и F силы, действующие на тело

во время контакта со стороны полупространства. Уравнения движения шара запишутся в виде [8], [18]

m(v2 -vz0) = ~\'qFz(t')dt\ m(vx -vx0) = -\'qFx(t')dt I (ю-ю0) = -R\[Fx (t')dt

(1.1) (1.2) (1.3)

где г - время столкновения, а момент инерции сферы I = (2/5)тЯ2. Система уравнений (1.1)—(1.3) позволяет определить все параметры после столкновения (скорость V, угол скольжения (3 , угловую скорость со) как функции аналогичных параметров до столкновения. Однако для этого необходимо знать временные зависимости тангенциальной Е (0 и нормальной Е (г) сил, определение которых всегда представляет нетривиальную задачу. При использовании условия качения для тангенциальной компоненты скорости отражения V

v + raR - 0

(1.4)

система (1.1)—(1.3) решается аналитически точно, что приводит к известному результату [19]:

vx vx0 R®0>

_ 2

5 v„

га -—га„---—

7 0 7 R

(1.5)

(1.6)

который мы в дальнейшем будем обозначать как «классическое» решение. Изменение кинетической энергии в результате столкновения составит

АЕ = от(уX - ^20)/2 +1 -®о)/2 =

(1.7)

Из выражения (1.7) видно, что удар не является абсолютно упругим, поскольку энергия не сохраняется (ее изменение отрицательно). Отметим, что решения (1.5), (1.6) не содержат материальных параметров типа модулей сдвига, а зависят только от начальных условий, что, как показано далее, имеет место только для узкого диапазона параметров и в большинстве практических случаев несправедливо [8].

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

еп / ^0 >

= V: / ^0,

(1.8) (1.9)

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

2. Моделирование упругого удара сферического тела с упругим полупространством

Несмотря на обширное количество экспериментальных результатов и теоретических моделей, столкновение сферического тела под углом к поверхности на сегодняшний день представляет нерешенную задачу, поскольку не существует единой теории, позволяющей полностью описать процессы, происходящие в зоне контакта. В недавних работах авторов настоящего обзора [20-22] в безразмерных переменных решена задача по столкновению упругого тела с упругим полупространством (см. рис. 1) в таких случаях: столкновение без проскальзывания (эквивалент бесконечной силы трения), столкновение с конечным коэффициентом трения, а также столкновение при наличии адгезионного взаимодействия между поверхностями контактирующих тел. Причем задачи решены таким образом, что, зная начальные параметры до столкновения (скорость, угол падения и циклическую частоту вращения), аналогичные величины после столкновения определяются однозначным образом. Этот раздел посвящен обзору результатов, полученных в указанных работах.

2.1. Основные предположения

Рассмотрим систему, показанную на рис. 1. Основным моментом при построении математической модели здесь является адекватное описание нормального и тан-

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

1 - 2у1 _ 1 - 2У2 С2

(2.1)

Будем действовать в приближении полупространства [23], в котором градиент поверхности всегда намного меньше 1, что означает малость глубины индентирования d по сравнению с радиусом контакта а. Кроме того, все размеры, связанные с контактом, считаются пренебрежимо малыми по сравнению с радиусом сферы Я . Если нормальная и тангенциальная компоненты скорости уг и V, то характерное отклонение центра масс индентируемого тела имеет порядок величины « d ■ V / V . Для «типичного» столкновения с отношением V / V порядка единицы тангенциальное перемещение ^ « d, а это намного меньше, чем радиус контакта.

Адгезию в контакте будем описывать в рамках JKR-теории, предполагая, что условие «прилипания» выполняется по всему контакту. Как показано в работе [24], при межмолекулярном притяжении предположение «прилипания» физически более обоснованно, чем предположение отсутствия сдвиговых напряжений в контакте, обычно используемое в теории ЖК Тем не менее в случае упруго подобных тел, удовлетворяющих уравнению (2.1), нет никакого различия между нормальным контактом без трения и контактом без проскальзывания. Даже в случае материалов, которые не удовлетворяют уравнению (2.1), теория может быть использована в очень хорошем приближении, поскольку разница получаемых результатов в адгезионном контакте со скольжением и без него по величине имеет тот же порядок, что и различие между нормальной жесткостью цилиндрического индентора с трением и без такового, для коэффициентов Пуассона, находящихся в пределах от 0,3 до 0,5, не превышает 3 % [24]. Отметим, что столкновение рассматривается в квазистатическом случае, когда скорости воздействия намного меньше, чем скорость звука в упругом полупространстве. В этом случае эффекты, вызванные инерцией полупространства, незначительны. Помимо всего сказанного, мы пренебрегаем потерей кинетической энергии при столкновении за счет упругих волн, распространяющихся в полупространстве [25, 26].

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

вращения К = / (МЯ2) [7] и угол Ф между нормалью к полупространству и направлением между точкой

контакта и центром масс параболического индентора. Здесь мы рассматриваем случай однородного шара, для которого К = у!0,4 и Ф = 0, поскольку именно этот случай представляет наибольший интерес в различных приложениях.

2.2. Удар без проскальзывания (бесконечный коэффициент трения)

В работе [7] показано, что в зависимости от угла падения возможно различное поведение в области между поверхностями контактирующих тел: полное проскальзывание, частичное проскальзывание, либо качение с преобладанием «прилипания» (эти три зоны показаны на рис. 8 в работе [8]. Описанные ситуации соответствуют конечному коэффициенту трения ц, который существенно отражается на поведении системы. Прежде чем вводить в рассмотрение коэффициент трения, опишем задачу без проскальзывания, решение которой представлено в работе [20]. Классическое решение этой задачи дается уравнениями (1.5)-(1.7). Покажем здесь, что это решение ввиду использованного приближения (1.4) выполняется только при определенном соотношении между материальными параметрами контактирующих тел, а в произвольном случае является неверным. Отметим, что ситуация без проскальзывания рассматривалась ранее в работе [27], однако в ней учитывались неупругие эффекты, мы же ограничимся областью упругого удара.

2.2.1. Случай постоянной контактной жесткости

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

Обозначим нормальное и тангенциальное смещение центра тяжести сферического тела как и2 (направлено вниз) и и^ (направлено вправо по ходу движения), нормальную и тангенциальную силы как ^ (направление совпадает с направлением и2) и ^ (направление совпадает с их). Тогда уравнения движения запишутся в виде (ср. с (1.1) - (1.3))

тй = —FM

тй = —F'

/ф = -FfR,

(2.2)

(2.3)

(2.4)

где ф - угол поворота сферы, изменение которого задает угловую скорость через связь ф = со. Нормальная и тангенциальная силы в рассматриваемом случае определятся как

F,, = ku

N z z

Ff = kx (Ux + Яф),

(2.5)

(2.6)

где к и кх - нормальная и тангенциальная жесткости гипотетической пружины. В уравнении (2.6) мы учли тот факт, что тангенциальное перемещение точки контакта является суммой перемещений центра масс и перемещения за счет вращения несжимаемой сферы. Решая систему уравнений (2.2)-(2.6) при начальных условиях в начальный момент контакта и (0) = 0, й, (0) = у 0, иг (0) = 0 , иг (0) = уг0, ф(0) = 0 , ф(0) = со0, мы получим

"z(0 = vzo cos(roz0,

(2.7)

5 2о

+ J (Vx0 + (2.8)

5vri 5

— co0---H--

7 7 R J 1R

Ф« = | | + ^(vl0+iiG>0)cos(tV), (2.9)

где частота % /т и частота % = ^7кх / (2т) . Время контакта г определяется уравнением (2.7) и вычисляется как = л / % = л^т / кг . Скорости, в последний момент столкновения соответствующие времени г , соответственно вычисляются так:

где

uz{ti) = -vz„

• / ч - 2тл 7 К

Vx = ux (ti ) = Vx +- V C0s Л - —

CD = ф(/,.) = CD Л--FCOS

7 R

V = vxo + >

(2.10)

(2.11)

(2.12)

а величины ух и и представляют классическое решение (1.5), (1.6). Согласно полученным соотношениям изменение энергии во время столкновения примет вид

№ = -m(vxo + ^®o)2sin2 nl7

(2.13)

Отметим, что в случае кг = 14кх решение задачи с постоянной жесткостью (2.10), (2.11) и (2.13) совпадают с классическим решением (1.5), (1.6), (1.7). Из выражений (2.10), (2.11) и (2.13) следует, что комбинации

7 Ух- Ух 7 кх

—х-- cos nl—-

2 V У 2 к

7 R(га-га)

V

7AE

- cos

7к-2 к

m(v-0 + r®g)

■ - - sin

7 к^ 2 к

(2.14)

(2.15)

(2.16)

зависят только от отношения тангенциальной жесткости к нормальной кх / к2. Поскольку такое отношение не

зависит от размера контакта и, следовательно, остается неизменным в процессе индентирования, это позволяет предположить, что структура уравнений (2.14)-(2.16) справедлива для более общего случая контакта произвольной формы. Для произвольного осесимметричного тела отношение тангенциальной жесткости к нормальной постоянно и вычисляется как кх / к2 = О* / Е* [28],

где Е * определяется выражением

J_

E*

1 -v? 1 -v2

E,

E

(2.17)

G*

справедливо аналогичное соотношение в виде

G

2-v, 2-v.

4G 4G2

(2.18)

В работах [29, 30] показано, что выражения (2.17), (2.18) справедливы в том числе и для шероховатых фрактальных поверхностей, которые имеют практически все природные объекты, формирующиеся под действием шумов (флуктуаций) и других случайных сил [31-34]. Таким образом, можно предполагать, что соотношения (2.14)-(2.16) выполняются как для гладких профилей, так и для шероховатых, для которых точная зависимость может быть заменена другой, учитываю-

щей тип шероховатости. Таким образом, мы приходим к гипотезе, что в общем случае соотношения (2.14)-(2.16) должны быть определены как

7 v- - v-2

7 R(ra-ra)

V 5

7AE

V

- Р(у),

m(v-0 + Rra0) где введена функция

--1 + (Р(у) )2,

7 G *

1 У 2 E • '

(2.19)

(2.20)

(2.21)

В следующих двух подразделах мы докажем правомерность гипотезы (2.19), (2.20), а также найдем форму функции P( у).

2.2.2. Методика моделирования

Для моделирования нормального и тангенциального контакта при столкновении будем использовать метод редукции размерности (Method of Dimensionality Reduction - MDR) [9]. В рамках MDR необходимо сделать следующее. Сначала произвести замену трехмерного упругого полупространства на линейное упругое основание, состоящее из массива независимых пружин с расстоянием Ах между ними. При этом для каждой пружины тангенциальная и нормальная жесткость рассчитываются как

Акг - E'Ах, Акх - G" Ах.

На следующем шаге трехмерный профиль z - f (r) трансформируется в плоский профиль g(х) согласно правилу

f(r)

g(х) -|х||-Т== dr.

0 V х - r

Этот этап схематически показан на рис. 2.

(2.22)

а б

Рис. 2. МБЯ-трансформирование оригинального трехмерного профиля в одномерный (а) и замена упругого полупространства дискретным упругим основанием (б)

Fig. 2. MDR transformation of the original three-dimensional profile into one-dimensional (a) and replacement of the elastic half-space by a discrete elastic foundation (b)

В данном подразделе мы рассматриваем случай без учета адгезии, поэтому показанная на рисунке величина А^тах (а) = 0. В [8, 13] показано, что при вдавливании МОЯ трансформированного профиля в упругое основание и его последующем произвольном движении в нормальном и тангенциальном направлениях зависимости обеих компонент сил от перемещений в соответствующих направлениях для одномерного представления задачи полностью воспроизводят зависимости для начальной трехмерной конфигурации. Решение контактной задачи в рамках метода MDR имеет такую же точность, как и решение Каттанео [35] и Миндлина [28]: в случае произвольного значения коэффициента Пуассона присутствует ошибка, которая, однако, является несущественной [36].

Для сферы радиусом Я профиль в окрестности контакта имеет вид

г2

/ (г) = — •

2 Я

Этой функции согласно (2.22) соответствует одномерный профиль

g ( x) = x2/ R.

(2.23)

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

йг(х1)=иг~8(х1)- (2.24)

Контактный радиус а определяется согласно условию

йг (a) =uz - g(a) = 0.

(2.25)

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

дйх ) = dих + ТМср.

(2.26)

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

Fn =Е* Аг^МДХ,.),

cont

Ff = G*AxJ^ux(xi),

(2.27) (2.28)

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

2.2.3. Результаты моделирования

При моделировании процесса столкновения уравнения (2.2)-(2.4) решались численно методом Эйлера с шагом интегрирования по времени Аг = 10—8 с, а тангенциальная ^ и нормальная ^ силы определялись

по методике, подробно описанной в подразд. 2.2.2. Результаты моделирования показаны на рис. 3.

Здесь показана рассчитанная функция Р(у) (2.19),

полностью определяющая поведение системы, а также безразмерное изменение энергии при столкновении, задающееся выражением (2.20). Обе величины представлены как функции параметра у (2.21). Если контактирующие тела имеют одинаковые коэффициенты Пуассона ^ = у2 = V, то выполняется равенство О*/ Е*= 2(1 — V) / (2 — V). Из термодинамического критерия устойчивости следует соотношение —1 1 / 2 [37]. Таким образом, для изотропных тел выполняется неравенство 2 / 3 < О* / Е* < 4 / 3, которое приводит к актуальному диапазону безразмерного параметра у :

1,52 <у< 2,16.

(2.29)

Однако для анизотропных (например, ортотропных) тел эффективное отношение О * / Е * может находиться в более широком диапазоне, чем (2.29), поэтому при моделировании мы не ограничивались указанным интервалом. В результате проведенного моделирования выяснено, что все данные ложатся на единственную кривую. Таким образом гипотезы (2.19), (2.20) подтверждены.

В области у < 3 численно найденная функция Р(у) с высокой точностью аппроксимируется выражением

Р (у) = —1 + 2ехр(—ау)соэ2 [к (у — Ь)] (2.30)

с параметрами а = 0,195 , Ь = 0,061, к = 1,685 . Данные, рассчитанные по выражению (2.30), показаны на рис. 3 ромбами.

Коэффициент детерминации для аппроксимации (2.30) и численно найденной функции Р(у) Я2 > 0,9995 . Воспользовавшись выражением (2.30) для практически важного случая V = 1 / 3, мы имеем Р ~ 0,197, а для несжимаемой среды (V = 1 / 2) получим Р ~ —0,087 . В актуальном диапазоне (2.29) с гораздо более высокой точностью можно использовать аппроксимацию

P(y) = 33,405 - 85,392у + + 75,914у2 - 28,195у3 + 3,727у4,

(2.31)

которая для величин v = 1 / 3 и v = 1 / 2 дает значения Р « 0,193 и Р «—0,1. Из рис. 3 следует, что классическое решение (1.5), (1.6) в актуальном диапазоне пара-

метров (2.29) является справедливым только для единственного значения параметра у.

Отметим, что функция Р(у) не определена, если + Я%) = 0. В этом случае имеем полностью нормальный контакт, поэтому, что очевидно, ух = и и = %.

Рис. 3. Зависимости Р(у) (2.19) и -1 + [Р(у)]2 (2.20) как функции параметра у (2.21). Сплошная кривая представляет собой результат большого количества численных экспериментов с различным набором начальных условий, радиусов, упругих модулей и массы сферического тела. Вне зависимости от значений указанных величин все данные ложатся на единственную кривую, которая однозначно связывает кинематические параметры после столкновения с аналогичными

параметрами до столкновения. Классическое решение (1.5) показано горизонтальной штриховой линией. Интервал (2.29) показан вертикальными штриховыми линиями Fig. 3. Dependences Р(у) (2.19) and -1 + [Р(у)]2 (2.20) as functions of the parameter у (2.21). The solid curve is the result of a large number of numerical experiments with a different set of initial conditions, radii, elastic moduli, and the mass of the spherical body. Regardless of the values of these quantities, all data located on a single curve that uniquely relates the kinematic parameters after collision with similar parameters bevor collision. The classical solution (1.5) is shown by a horizontal dashed line. The interval (2.29) is shown by vertical dashed lines

2.3. Удар с конечным коэффициентом трения

В подразд. 2.2 проведен анализ столкновения упругого сферического тела с упругим полупространством в случае бесконечного коэффициента трения, когда по всей области контакта наблюдается «прилипание» без проскальзывания. Это предположение существенно сужает область применения полученных результатов, поскольку достоверно известно, что в определенные моменты контакта наблюдаются области проскальзывания. Например, это показано на диаграмме на рис. 8 в работе [8]. Используемый нами для определения нормальной и тангенциальной контактных сил метод MDR позволяет проводить моделирование с произвольным коэффициентом трения. Проведем в этом подразделе анализ, позволяющий обобщить полученные ранее результаты на более широкий класс рассматриваемых систем.

2.3.1. Случай постоянной контактной жесткости

Пусть между контактирующими поверхностями существует трение, характеризуемое коэффициентом трения ц. Проведем здесь анализ, аналогичный осуществленному в подразд. 2.1.1. Полученные ранее результаты (2.10), (2.11), (2.13) справедливы в случае, когда в системе с конечным коэффициентом трения ц выполняется условие «прилипания», при котором во всей зоне контакта отсутствует проскальзывание. Это условие записывается в виде

Vx0 + Яюо

|sin(ra/)| < ц|sin(raz0|• (2.32)

k„ га.

Уравнение (2.32) имеет нетривиальное решение:

% = у%, а>ас2 = 2у2 / 7, (2.33) задающееся параметром

а = Ц^г0 / I ^0 + Я®0| (2.34)

и натуральным числом у (натуральное значение этот параметр имеет только в рассматриваемом случае постоянной жесткости [22]). При выполнении условия (2.33) реализуется случай, описывающийся уравнениями (2.10), (2.11), (2.13). Рассмотрим ситуацию с произвольным коэффициентом трения ц. Условие «прилипания» задается комбинацией параметров у (2.21) и а (2.34). Это позволяет предположить, что гипотезы (2.19), (2.20) действуют и в рассматриваемом случае, однако теперь необходимо рассматривать функцию двух параметров Р(у, а). Далее мы докажем эту гипотезу численным моделированием и определим вид функции Р(у, а), но прежде получим

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

При выполнении неравенства

а < ас2 = 2у2 / 7

(2.35)

пружина будет скользить с начала контакта. В этом случае сила трения определится как

р/ = ц^« 0 + я^Х

где - знаковая функция, определяющаяся стан-

дартным образом и введенная для учета направления действия силы. Начальная фаза проскальзывания закончится в момент времени г = ^, которое является корнем уравнения

1 - "Г I созКО = 1 -.

у

7а

Поскольку контакт существует до момента времени г = л / юг, пружина будет проскальзывать на всем его протяжении, если ^ > я/юг, т.е. при выполнении неравенства

2у2

7(2у2 -1)

(2.36)

В этом случае горизонтальная скорость после столкновения может быть рассчитана как

V = ^0 — — 1р/& = ^0 — — J рм& вБИ^ + Я®0) = т т

= Ка — 2^0 ^Ке + я®0)- (2.37)

Выражение (2.37) может быть переписано в безразмерном виде как

Р = 1 — 7а (2.38)

Отметим, что решение (2.38) будет справедливо также и для нелинейного контакта в случае скольжения по всей зоне контакта на протяжении воздействия, поскольку зависимость нормальной силы р от глубины индентирования в (2.37) не была задана.

2.3.2. Методика моделирования

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

I G*mv(X,.)|.

(2.39)

Если условие (2.39) выполняется, пружина будет скользить. Для учета скольжения на каждом шаге к горизонтальному смещению пружин прибавляется инкремент (2.26), после чего для всех пружин проверяется условие (2.39). Для тех пружин, для которых это условие выполняется, устанавливается новая величина смещения

(2.40)

где знак совпадает со знаком их (х) в вьфажении (2.39), т.е. до присваивания величине их(х) нового значения. После этого нормальная р, и тангенциальная р силы рассчитываются, как и ранее, по формулам (2.27), (2.28). Поскольку р и р теперь определены, решая

численно систему уравнений (2.2)-(2.4), мы можем определить временную эволюцию интересующих нас величин. Как и ранее, сосредоточимся на нахождении функции Р( у, о) (2.19), которая позволяет рассчитать полный

набор кинетических параметров после столкновения в зависимости от их значений до момента контакта.

2.3.3. Результаты моделирования

В разд. 4 работы [8] на основе имеющихся литературных данных описано столкновение сферического тела с упругим полупространством под углом. В частности, кратко описаны результаты работы [7], где показано, что в зависимости от угла падения возможна реализация различных ситуаций - полное проскальзывание в зоне контакта (область III), преобладание «прилипания» (область I), либо частичное проскальзывание (область II). На рис. 8 в работе [8] эти области показаны римскими цифрами. Ситуация с бесконечным коэффициентом трения, когда проскальзывание отсутствует на всех этапах столкновения, рассмотрена нами в предлагаемой работе в под-разд. 2.2, где в актуальном диапазоне (2.29) функция Р( у) описывается аппроксимацией (2.31) и не зависит от параметра о (2.34). При скольжении по всей зоне контакта реализуется соотношение (2.38).

Прежде чем начинать моделирование процесса столкновения, определим, пользуясь выражениями (3.6)-(3.12) из работы [8] области I, II и III, как функции параметров у (2.19) и о (2.34). При этом условие реализации области I (преобладает "прилипание") запишется как а > ас2. Интервал, соответствующий области II (частичное проскальзывание) в рассматриваемом случае принимает вид аЛ < а < ас 2. И последняя область III (полное проскальзывание) наблюдается при выполнении неравенства а < ас1. Критические значения ас1, ас2 определяются полученными ранее формулами (2.36), (2.35). Отметим, что в случае у = 1 обе критические величины совпадают (ас1 =ас2), поэтому описанное выше справедливо только при у > 1, когда ас2 > ас1. На рис. 4 показана диаграмма, на которой показаны обсуждаемые выше области I, II и III.

Рис. 4. Диаграмма режимов столкновения с конечным коэффициентом трения. Области I, II, III описаны выше и соответствуют аналогичным областям, показанным на рис. 8 в работе [8]. На вставке к рисунку показана актуальная область (2.29)

Fig. 4. Diagram of collision modes with a finite coefficient of friction. Areas I, II, III are described above and correspond to similar areas shown in Fig. 8 in [8]. The inset in the figure shows the actual area (2.29)

На осях диаграммы отложены безразмерные переменные а и у, набор которых полностью определяет

рассматриваемую задачу. Увеличенный фрагмент диаграммы, показанный на вставке к рисунку для актуального диапазона параметров (2.29), полностью определяет возможные режимы в этом диапазоне, поскольку показывает поведение для всех значений а.

На рис. 5 показана полученная в результате численного эксперимента зависимость P( у, а), где результаты многочисленных компьютерных экспериментов ложатся на единственную поверхность, что делает ее универсальной.

Из рисунка видно, что при больших значениях параметра а (2.34) поведение описывается кривой, показанной на верхней панели рис. 3, который построен для случая с бесконечным коэффициентом трения, соответствующим отсутствию проскальзывания в зоне контакта. Однако, как видно из рис. 5, поведение в области малых значений параметра а существенно отличается, что требует проведения дополнительного анализа.

На рис. 6, а приведены зависимости P(у, а) как функции параметра а (2.34) для 17 значений у (2.21) из актуального диапазона (2.29). Из рисунка видно, что для всех кривых при малых а наблюдается универсальное поведение, не зависящее от величины у . Здесь все кривые ложатся на единственную прямую, задающуюся уравнением (2.38). Это связано с тем, что такие значения параметра а отвечают низкому трению, когда осуществляется проскальзывание на протяжении всего контакта (область III на диаграмме (см. рис. 4)). При увеличении значения а функция P(у, а) отклоняется от этой прямой, причем, как видно на вставке к рис. 5, которая представляет его увеличенный фрагмент, где и начинаются указанные отклонения, с увеличением ве-

личины параметра у отклонение начинается раньше. Это согласуется с диаграммой, представленной на рис. 4, где с ростом у область проскальзывания III сужается, поскольку уменьшается критическое значение ас1 (2.36). После отклонения от линейного поведения все кривые ведут себя одинаково: функция Р( у, а) сначала продолжает убывать, а затем при достижении минимума возрастает (возможно и последующее убывание функции), и с превышением параметром а критического значения выходит на постоянное значение, соответствующее полному прилипанию, когда Р(у, а) больше не зависит от а и описывается аппроксимацией (2.31). Описанное поведение при больших а хорошо прослеживается на трехмерной диаграмме на рис. 5. Стоит отметить, что условие а > ас2 означает преобладание "прилипания" в зоне контакта, а полное прилипание наступает при еще большем увеличении а. Поэтому при превышении а > а и дальнейшем его увеличении некоторое время наблюдается зависимость функции Р(у, а) от параметра а. К примеру, на кривых на рис. 6, а при малых значениях параметра у во всем приведенном диапазоне а > а наблюдается увеличение Р(у, а), а кривые, соответствующие большим у , при а > ас2 показывают уменьшение значений Р(у, а). С еще большим ростом а все же устанавливается постоянное значение Р(у, а), описываемое аппроксимацией (2.31), что легко видеть на рис. 6, б. Однако, из рисунка следует, что даже при значении а = 12 зависимость Р( у, а) хоть и незначительно, но все же отклоняется от функции (2.31).

2.3.4. Сравнение с экспериментом

На рис. 7 показаны экспериментальные результаты исследования соударения под углом сферического тела с плоскостью, полученные в работе [38].

Рис. 5. Численно рассчитанная зависимость Р(у,а) (2.19) как функция параметров у (2.21) и а (2.34). Полученная поверхность является результатом большого количества численных экспериментов с различным набором начальных условий, радиусов, упругих модулей, и массы сферического тела

Fig. 5. Numerically calculated dependence P(y,a) (2.19) as a function of parameters y (2.21) and a (2.34). The resulting surface is the result of a large number of numerical experiments with a different set of initial conditions, radii, elastic moduli, and the mass of a spherical body

а б

Рис. 6. Зависимости Р(у,а) как функции параметра а (2.34) при фиксированных значениях у (2.21) (а). Величина у изменяется от 1,52 до 2,16 с интервалом 0,04. Направление увеличения у показано на рисунке стрелкой. Штриховой линией показана критическая величина ас2

(2.35). Сплошная прямая на вставке к рисунку иллюстрирует зависимость (2.38); зависимости Р(у, а) как функции параметра у (2.21) при фиксированных значениях а (2.34) (б). Нижние 8 кривых построены при увеличении а от 0,3 до 1,0 с шагом 0,1. Далее а увеличивается

от 1 до 12 с шагом 1,0. Направление увеличения а показано на рисунке стрелкой. Штриховой линией показана функция (2.31) Fig. 6. The dependence Р(у, а) as a function of the parameter а (2.34) for fixed values of the parameter у (2.21) (a). The value of у varies from 1.52 to 2.16 with a step of 0.04. The direction of increasing of у is shown in the figure with the arrow. The dashed line shows the critical value ас2 (2.35). The solid straight line in the inset to the figure illustrates the dependence (2.38); the dependence Р(у, а) as a function of the parameter у (2.21) for fixed values а (2.34) (b). The lower 8 curves are plotted with an increase from 0.3 to 1.0 in 0.1 increments. Then а increases from 1 to 12 in step 1.0. The direction of increasing of а is shown in the figure by the arrow. The dashed line shows the function (2.31)

а б

Рис. 7. Тангенциальный коэффициент восстановления et (1.9) (а) и угловая скорость -ю (б) как функции угла падения (90 -а°) в области упругого удара [38]. Кривая 1 на панели (а) получена по формуле (2.42) при ц = 0,092 , en = 1,0 . Сплошные кривые

на обоих рисунках рассчитаны способом, подробно описанным в подразд. 2.3.4 Fig. 7. The tangential restitution coefficient (1.9) (a) and the angular velocity -ю (b) as a function of the angle of incidence (90-а°) in the region of elastic impact [38]. Curve 1 on the panel (a) is obtained from (2.42) at the parameters ц = 0,092 , en = 1,0 . The solid curves in both figures are calculated in the manner detailed in subsection 2.3.4

В работе изучается отражение шарика, выполненного из оксида алюминия, диаметром 5 мм, падающего под углом на упругую поверхность при нулевой начальной частоте вращения % = 0 . На рис. 7 показаны результаты экспериментов по столкновению шарика с отражающей поверхностью в области упругого удара

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

В случае, показанном на рис. 7, нормальный коэффициент восстановления ел (1.8) составлял от 0,976 до

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

0,986, что говорит о полностью упругом ударе. Из рис. 7, а следует, что минимальное значение тангенциального коэффициента восстановления е( « 0,6 реализуется при значении угла падения порядка 20°, а при полном скольжении, когда скорость имеет только тангенциальную составляющую (а° = 0°), коэффициент

восстановления г « 12. Каждая точка на рис. 7 получена усреднением минимум 10 экспериментов с фиксированными значениями начальной скорости V и угла падения (90 — а°). На рис. 7, б показана зависимость циклической частоты вращения шарика после столкновения от угла падения, которая определяется как [38, 39] (см. обозначения на рис. 1)

5 у0 сое а°

ю = —0-(1 — г )•

2 Я '

(2.41)

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

справедливо соотношение между нормальной р и тангенциальной р силами в виде

р = цр,

из чего легко получить выражение для тангенциального коэффициента восстановления [40]:

е, =1-Ц(1 + е„)*8а%

(2.42)

где для эксперимента, показанного на рис. 7, при углах падения, находящихся в диапазоне между 90 и 30°, коэффициент трения ц« 0,092. Коэффициент г, рассчитанный по формуле (2.42), соответствующий значениям ц = 0,092 и г = 1,0, показан на рис. 7, а сплошной лини-ей3, откуда видно, что приближение несжимаемого тела хорошо работает в данном случае для углов падения (90° —а°) > 30°. Это связано с тем, что в зоне контакта

2Подобное поведение наблюдалось в эксперименте, описанном в работе [41].

3Несмотря на то, что в эксперименте коэффициент г принимал значения от 0,976 до 0,986, а его среднее значение составляло 0,980, при расчете зависимости (2.42) выбрано значение еп = 1, соответствующее упругому нормальному столкновению, поскольку далее будет показано, что в области применимости соотношения (2.42) результаты моделирования, проведенного в подразд. 2.3, показанные на рис. 7 сплошными непронумерованными линиями, идеально описываются формулой (2.42).

наблюдаются области проскальзывания и «прилипания». Однако при больших значениях углов падения (когда нормальная скорость существенно меньше тангенциальной ) преобладает зона проскальзывания, и эта ситуация хорошо описывается аппроксимацией (2.42).

Поскольку данные получены в области отсутствия пластических эффектов, мы можем провести сравнение этих результатов с результатами моделирования в рамках метода MDR, проведенного в подразд. 2.3.3. В указанном подразделе рассчитана функция Р(у, а) (2.19), которая позволяет рассчитать все параметры столкновения. Причем переменная у (2.21) учитывает параметры соударяющихся тел, а а (2.34) описывает влияние коэффициента трения ц. На рис. 7 приведены зависимости тангенциального коэффициента восстановления г = V / ^о (19) и циклической скорости ю, с которой сфера отскакивает от поверхности, при нулевом начальном вращении ю0 = 0. Выразим эти параметры

через функцию Р(у, а). Воспользовавшись соотношениями (2.19), (2.20), (1.5), (1.6), (2.12) при ю0 = 0, получим

2Р(у, а) + 5

7

ю =

5у0 сое а [Р(у, а) -1]

7 Я

(2.43)

(2.44)

где в эксперименте, описанном в работе [38], радиус сферы Я = 2,5 • 103 м, а начальная скорость V = 3,9 м/с. Угол а° в формуле (2.44) - это угол скольжения, показанный на рис. 1. В работе [38] изменялся угол падения, который составляет (90 — а°) и отложен по оси абсцисс на рис. 7. Величина а в формулах (2.43), (2.44) зависит от угла скольжения а°. Так, в случае ю0 = 0 параметр а определится согласно (2.34) как

а =

(2.45)

где в эксперименте, представленном на рис. 7, коэффициент трения ц = 0,092. Таким образом, используя формулы (2.43), (2.44) и варьируя угол а°, можно получить зависимости, представленные на рис. 7, если знать вид функции Р(у, а), показанной на рис. 5 и 6. Перед началом соответствующих расчетов необходимо определить параметр у (2.21), который остается постоянным для серии экспериментов, результаты которых показаны на рис. 7. В работе [38] не приведены параметры используемых материалов, однако указаны сами материалы. Так, в качестве отражающей поверхности использовалось натриево-кальциевое стекло, а отражаемая частица была изготовлена из оксида алюминия. Более подробных данных в работе нет, а приведенных

недостаточно для точного определения параметров, поэтому возьмем типичные для этих материалов вели-чины4. Так, для используемого стекла Е2 = 7,3 • 1010 Па, О = 3,0 • 1010 Па, у2 = 0,217 . Для оксида алюминия Е = 3,7• 1011 Па, О = 1,53• 1011 Па, 0.25 и плотность р = 3970 кг/м3, что соответствует материалу, содержащему 96 % Л1203 при комнатной температуре. Указанный набор параметров позволяет найти Е* и О* (формулы (2.17), (2.18)), Е* ж 6,42-101» Па, О* « 5,64-1010 Па, что согласно формуле (2.21) приводит к значению у «1,754.

На рис. 7, а и б сплошными линиями показаны рассчитанные по формулам (2.43), (2.44) тангенциальная компонента коэффициента восстановления в и циклическая частота и, взятая с обратным знаком, при у = 1,754 и изменении а согласно выражению (2.45). Из рисунков следует, что полученная нами зависимость с большой точностью описывает экспериментальные данные во всем диапазоне углов падения, что является несомненным преимуществом по сравнению с аппроксимацией (2.42), которая показана на рис. 7, а кривой 1 и описывает поведение системы при углах падения (90° - а) > 30° . В то же время полученные нами результаты и аппроксимация (2.42) идеально совпадают в области применимости (2.42). Преимущество нашего подхода состоит еще и в том, что в нем начальная циклическая частота вращения сферической частицы может принимать любые отличные от нуля значения. Также важен и тот факт, что исследование проведено в безразмерных величинах у (2.21) и а (2.34), которые включают в себя все параметры системы. Поэтому в области упругого удара подход на основе метода MDR, описанный в данной работе, позволяет описать любой эксперимент. Более того, как это предложено в работе [42], в случае неизвестной природы контактирующих тел результаты экспериментов по столкновению частиц можно использовать для определения модуля упругости материала, для чего можно использовать соотношения (2.43), (2.44), поскольку вид функции Р(у, а) известен.

2.4. Удар с адгезией

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

4При вариации параметров для различных типов стекол и про-

центного состава Л1203 получаемые результаты изменяются несуще-

ственно.

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

В экспериментальной работе [43] показано, что нормальный и тангенциальный вклады в силы адгезии могут зависеть от направления движения частицы. В предлагаемом нами подходе мы считаем, что реализуются условия, при которых влияние тангенциального нагружения на адгезию в нормальном направлении отсутствует. В этом случае адгезионное взаимодействие может быть описано в рамках метода редукции размерности [8, 13, 14]. В [44] изучалось влияние адгезии на столкновения частиц, находящихся в жидкости с использованием теории JKR [45]. Авторы указанной работы предложили комплексное решение задачи в безразмерном виде для широкого диапазона параметров в случае чисто нормального контакта. В другой работе [46] при использовании метода дискретных элементов (DEM) изучено взаимодействие между полимерными частицами, которые моделируются как агломерации упругих микроэлементов, связанных с помощью адгезионных сил. Проведено исследование длительности столкновений, а также процессов агломерации и распада полимерных частиц во время столкновения. Задача столкновения для нормального адгезионного контакта также изучалась аналитически и экспериментально в работе [47]. В этой работе применялись адгезионные модели JKR [45], DMT [48] и Маугиса [49]. Довольно полная аналитическая модель нормального столкновения с адгезией в рамках JKR-подхода предложена в работе [50]. В [51] проведено моделирование нормального адгезионного контакта с помощью методов молекулярной динамики. И хотя существуют теоретические работы, в той или иной мере описывающие адгезионный контакт при столкновении под углом [52, 54], этот вопрос ввиду своей сложности до сих пор остается открытым. Проведем здесь соответствующий анализ в рамках метода MDR, следуя идеологии под-разд. 2.2, 2.3. Будем считать, что в случае наличия адгезии между контактирующими поверхностями по всей области контакта наблюдается зона прилипания, что эквивалентно бесконечному коэффициенту трения. В работе [21] на основании анализа, представленного в подразд. 2.2, проведен учет адгезии в рамках метода MDR [9, 13]. Приведем в этом подразделе основные результаты, полученные в указанной работе.

2.4.1. Случай постоянной контактной жесткости

В рассматриваемом случае справедливы уравнения (2.5)-(2.9). Будем характеризовать адгезионное взаимодействие между контактирующими телами отрицательным значением «глубины проникновения» -dc, при котором адгезионный контакт теряет устойчивость и происходит отрыв сферы от пружины. Момент време-

ни ^, при котором происходит отрыв, легко получить из уравнения для вертикального движения (2.7):

и (и) = Vо эт(юД.) / ю2 = —dc,

откуда следует решение:

t; =— arcsm Ю

-d ю.

1

ю

л + arcsin

( d ю V

. (2.46)

Скорости в последний момент столкновения принимают вид

Z Z

'ю / V,

7 ^ 2 к

)2 =-Vz„V1 -ß2 ,

f ( dm Y|

л + arcsin

V V Vz o ))

= F + 2Г cos[9(y,ß)] /7,

(2.48)

® = Ф(А) = ю + V сое [ф(у, Р)], (2.49) IК

где введены функции р = dcюz / и ф(у,Р) =

= у(л + агс5т Р), а также используются величины у

(2.21), V; (1.5), ю (1.6) и V (2.12). Определим потери

энергии при контакте. В рассматриваемом случае удобно отдельно вычислять потери кинетической энергии при вертикальном ДЕг и горизонтальном ЛЕх движении:

Щ =- m (Vxo + R®o)2sin2 [9(y,ß)], (2.50)

AE = -mv\ß2 /2 = -E nß2,

z zzOr '

(2.51)

где Е - начальная кинетическая энергия, связанная с нормальным движением. Уравнения (2.47)-(2.51) могут быть переписаны в виде

vz /vz0 = -/1-^

AEz / Ezo =-ß2,

7 v - V 7 Я(ю-ю)

(2.52)

(2.53)

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

Р(у, Р). Дополнительно к этим двум гипотезам в случае

адгезионного контакта добавляются еще две, связанные с диссипацией нормальной компоненты кинетической энергии:

Vz / Vzo =-F(ß),

AEZ /o =-1 + (F(ß))2

(2.56)

(2.57)

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

d =

/ т Л 2/5

15 mVzo ^

16 E*R1'2

а максимальное расстояние между сферой и полупространством в момент отрыва в адгезионном контакте согласно JKR-теории [45] - как

А 3

d = 4

'n2RAy ^1/3

где Ay - поверхностная энергия. В итоге приходим к соотношению

d

ß = :i£. = л2/3-| -

3 (16Y'5 R8/15Ay2

d

4115

E* 4/15m2/5v:5 •

Для удобства опустим константу и будем в дальнейшем использовать параметр Р в виде

ß =

/ л ч Л 2/15

' R4Ay5 ^

с*2 3 6

E mvzo у

(2.58)

2 К 5 V

7AE„

= cos[ф(y,ß)], (2.54) 2.4.2. Методика моделирования

= -sin2 [ф(у,Ю]. (2.55)

т( ^ + Яю0)

Правые части выражений (2.52)-(2.55) являются функциями параметра у = 7£./2£г =4тОГп2Е" (2.21) и отношения смещений

р = < ю / V0 = ^ /d,

где d = / ю имеет смысл максимальной глубины проникновения (в случае постоянной жесткости максимальная глубина проникновения одинакова при учете

Если глубина индентирования профиля в упругое основание равна d, то смещения индивидуальных пружин в области индентирования определяются согласно выражению (ср. с (2.23), (2.24))

uz (x) = d - g(x) = d - x2 / R.

(2.59)

В адгезионном контакте контактный радиус a увеличивается, и выражение (2.25) теряет применимость. В этом случае натяжение граничных пружин может быть найдено по «правилу Гесса» (Heß rule) [9]. Согласно этому правилу, условие равновесия граничных пружин (при x = +a) может быть записано как

V

Al = A/max (a) = V2™Ду / E

(2.60)

где величина A/max (a) показана на рис. 2. Комбинируя (2.59) и (2.60), получим

иг (a) = d - a2/ R = -Al^ (a) = —/ 2rcaAy / E*, (2.61)

откуда

d = а2/ Я 2гоаАу / Е*. (2.62)

Нормальная сила при этом определится как сумма сил от всех сжатых (в области индентирования) и растянутых (в адгезионной области) пружин:

a a /

F (a) = J u (x)dx = 2E* JI d

R

dx =

4E*a ' 3R

8лл3 Е*Ду.

(2.63)

Представим аналитические и численные результаты в безразмерных переменных:

а = а/а0, Г_=Г_/Г0, d = dld0,

где Е0, а и ^ - критические значения силы, контактного радиуса и абсолютного значения глубины инден-тирования в момент отрыва параболического профиля от упругого полупространства при условиях «контролируемой силы» [9],

3

F0 a0 =

9tcR2 Ду 8E*

do =

3л2 RA у2

"б4Е"

В безразмерных единицах соотношения (2.62) и (2.63) запишутся как

= Зй2 -4д

F =аъ-2аъ

(2.64)

(2.65)

Выражения (2.64), (2.65) воспроизводят классическое решение, полученное в [45], и могут быть использованы для тестирования численного алгоритма. Так, на рис. 8 показаны результаты тестирования численного алгоритма для нормального контакта, из чего видно, что численные результаты совпадают с аналитическим решением с большой точностью.

Рис. 8. Нормальная сила F. как функция глубины

индентирования d для нормального контакта с адгезией. Сплошная кривая показывает аналитическое решение (2.64), (2.65). Кружками показаны результаты численного моделирования, полученные в рамках метода MDR

Fig. 8. Normal force F_ as a function of the indentation depth d lor normal contact with adhesion. The solid curve shows the analytical solution (2.64), (2.65). The circles show the results of numerical simulation obtained in the framework of the MDR method В итоге численное моделирование сводится к следующим шагам. На каждом временном шаге рассчитываются тангенциальные и нормальные натяжения всех индивидуальных пружин. В предположении отсутствия проскальзывания смещения пружин с координатами x рассчитываются по формулам (2.24) и (2.26). При этом пружины, находящиеся в контакте, определяются условием (2.61). После этого рассчитывается новый контактный радиус a. Полные контактные силы Fw и F

на каждом шаге вычисляются как суммы (2.27), (2.28). На последнем этапе временного шага численно методом Эйлера решаются уравнения движения (2.2)-(2.4). В модели могут изменяться параметры R, Е*, О*, m, Aу, а также обе компоненты начальной скорости vl0, vz0 и начальная циклическая частота вращения частицы ю0.

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

в подразд. 2.3. Случай бесконечного коэффициента трения, рассмотренный в подразд. 2.2, может быть описан в рамках используемой здесь процедуры при ß = 0, что

соответствует значению Ay = 0. В работе [7] вводится в рассмотрение безразмерная тангенциальная сила

- _ 3(1 -\>)НХ

V ,, , „2 ,

4CRlG\i'

(2.66)

которая представляется как функция безразмерного времени I = I / х . где время

т = 2,9432С2R / v ,

(2.67)

0

-a

соответствует длительности упругого удара без адгезии [53]. Также в [7] вводится параметр

у = 2(1 ,

Ц(2 -v)vz 0

(2.68)

значение которого определяет режимы удара. Константа С в выражениях (2.66), (2.67) имеет вид [7]

C =

f15m(1 -v)vz20^ 16GR3

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

Кривые 1, 2 и 3 на верхней панели рис. 9 соответствуют различным значениям параметра у < 1, при которых преобладает режим «прилипания». Кружками здесь показаны результаты, полученные в работе [7] для случая без адгезии. Отличия с рассчитанными кривыми связано с тем, что мы проводили моделирование с бесконечным коэффициентом трения без проскальзывания, в то время как при конечном ц в диапазоне значений у < 1, хотя и преобладает режим «прилипания», все же есть моменты времени, в которые осуществляется скольжение. Поэтому на зависимостях р(Г), полученных в работе [7], сила р меньше по абсолютному значению, чем в ситуации без проскальзывания. При использовании процедуры (2.39), (2.40) и выражений (2.66), (2.68) мы получим идеальное совпадение результатов, однако в рамках данного подраздела нас интересует случай без проскальзывания, поскольку рассматривается влияние адгезии.

Рис. 9. Зависимости безразмерной тангенциальной силы Fy (2.66) от безразмерного времени t = t / т (2.67), полученные при параметрах R = 0,02 м, Е' = 1010 Па, р = 1000 кг/м3, ю0 = 0 рад/с, vz0 = 0,01 м/с, v = 0,3 . Верхняя панель (случай без адгезии): кривые 1-3 отвечают значениям у = 0,2; 0,5; 1,0 , кружками показаны данные из работы [7]. Нижняя панель (случай с адгезией) построена для у = 0,5 и поверхностных энергий Ду = 1,0; 5,0; 10,0; 17,0; 24,0 Дж/м2. Штриховая кривая на обеих панелях соответствует случаю без адгезии при \|/ = 0,5 Fig. 9. Dependences of the dimensionless tangential force Fx (2.66) on the dimensionless time t=t/x (2.67), obtained with the parameters R = 0.02 m, E* = 1010 Pa, р = 1000 kg/m3, ю0 = 0 rad/s, vz0 = 0.01 m/s, v = 0.3 . Top panel (case without adhesion): curves 1-3 correspond to values у = 0.2,0.5,1.0 , circles show data from [7]. The lower panel (the case with adhesion) is constructed for у = 0.5 and surface energies Ду = 1.0,5.0,10.0,17.0, 24.0 J/m2 . The dashed curve on both panels corresponds to the case without adhesion at у = 0.5

Нижняя панель на рис. 9 показывает результаты моделирования для случая с адгезией при различных значениях поверхностной энергии Ду. Из рисунка следует, что с уменьшением значения Ду поведение системы стремится к безадгезионному случаю. Помимо этого, рост величины Ду приводит к увеличению продолжительности контакта, а также частоты и амплитуды колебаний тангенциальной силы Fx . Окончания всех зависимостей соответствуют моменту времени t перед отрывом профиля от упругого основания. Стоит обратить внимание, что в случае без адгезии в момент отрыва

реализуется значение тангенциальной силы Fx = 0, в то время как в адгезионном случае в конце контакта /*" Ф 0. Значение нормальной силы F в момент отрыва при Ду Ф 0 также ненулевое (см. рис. 8).

2.4.3. Результаты моделирования

На рис. 10 показаны зависимости нормального коэффициента восстановления vz/vz0 (2.56), а также относительной потери кинетической энергии ДЕг / Ez0 (2.57) от параметра Р (2.58), описывающего адгезионные свойства системы.

Рис. 10. Зависимости нормального коэффициента восстановления vz / vz0 (2.56) (сплошная линия) и относительной потери энергии за счет движения в нормальном направлении ДЕ / Ez0 (2.57) (пунктирная линия) как функции параметра Р

(2.58). Столкновения моделировались при различных наборах начальных условий и параметров модели, таких как упругие модули, масса и поверхностная энергия. Все данные ложатся на единственную кривую. Аппроксимация (2.69) показана

ромбами

Fig. 10. Dependences of the normal restitution coefficient vz / vz0 (2.56) (solid line) and the relative energy loss due to motion in the normal direction ДЕ / Ez0 (2.57) (dashed line) as a function of

the parameter P (2.58). Collisions were modeled for different sets of initial conditions and model parameters, such as elastic moduli, mass and surface energy. All data located on a single curve. The approximation (2.69) is shown by rhombuses

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

ДЕ

Е,

= -0,0019^ + 0,0413Р-2,7194Р2 -

-24,2169Р3 + 30,073р4 -23,5737р5,

(2.69)

из которого вертикальная компонента скорости после столкновения уг определяется из (2.56), (2.57),

а непосредственно функция Е (Р) (2.56) запишется как

F (P) = -v, / Vо ЧДЕ: / Ez0 +1.

(2.70)

Коэффициент детерминации для аппроксимации (2.69) и численно определенной кривой на всем показанном диапазоне Я2 > 0,99999997 (рис. 10), что гово-

рит о высокой скоррелированности выражения (2.69) с численными результатами. Ранее в работе [55] аналитически было получено точное выражение

АЕ- = -^(12 + )р5/2 ж -14,184Р5/2, (2.71)

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

Из расчета, показанного на рис. 10, следует, что критическое значение параметра Р - это Рс ж 0,34582 . Эту величину также можно найти аналитически точно из соотношения (2.71) при АЕг / Ег0 = -1, что для критического значения Рс приводит к следующему результату

Рс =

(12 + 21/3)2/5 л2/3

0.3461590663,

(2.72)

который хорошо совпадает с наиденным ранее численным значением.

В случае Р > Рс условие уг = 0 достигается до момента отрыва, поэтому сфера «прилипает» к поверхности. Аналогичный эффект наблюдается в рамках линейной модели (2.46), когда частица не отскакивает от поверхности при выполнении условия (^ / /т > 1. Поэтому будем исследовать область «отскока», в которой р < рс.

На рис. 11 показана полная зависимость Р(у,Р). Как и в предыдущих случаях, поверхность является результатом большого количества компьютерных экспериментов с различными наборами начальных значений и материальных параметров. В актуальном интервале значений у (2.29) показанная зависимость с достаточной точностью аппроксимируется выражением

г{ у, р)« а + а у+а у2 + а у3 + а у4

(2.73)

Рис. 11. Верхняя панель: трехмерная зависимость Р(у,Р) (2.19) как функция параметров у (2.21) и Р (2.58); нижняя панель:

двумерная зависимость Р(у,Р) (2.19), на которой большие значения P отвечают более светлым участкам Fig. 11. Top panel: three-dimensional dependence Р(у,Р) (2.19) as a function of the parameters у (2.21) and Р (2.58); lower panel: two-dimensional dependence Р(у,Р) (2.19), in which large values correspond to lighter areas

с коэффициентами а , которые являются функциями параметра Р [21]:

а = 43,139 -1171,292Р + 26272,84Р2 -425164,358Р3 +

+ 2019578,141Р4 - 2730636,482Р5, а = -107,979 + 2330,636Р -59003,404Р2 + 973668,408Р3 -

- 4542494,427Р4 + 6011602,434Р5,

а = 95,513 -1659,826Р + 49929,478Р2 -830694,428Р3 +

+ 3787497,609Р4 -4889171,548Р5, а = -35,686 + 507,716Р-18909,789Р2 + 312224,871Р3 -

- 1384531,709Р4 +1737264,447Р5,

а = 4,786 - 57,27Р + 2693,537Р2 - 43475,09Р3 + + 186768,848Р4 - 227031,806Р5.

На рис. 12 сплошными кривыми показаны численно полученные результаты, а символами - аппроксимация (2.73).

Рис. 12. Сравнение численных результатов (сплошные линии) с аппроксимацией (2.73) для Р(у,Р) как функции у (2.21) при фиксированных значениях параметра Р (2.58),

показанных на рисунке Fig. 12. Comparison of the numerical results (solid lines) with the approximation (2.73) for the function у (2.21) for fixed

values of the parameter Р (2.58), which are shown in the figure

Из рисунка видно, что функция (2.73) с высокой точностью описывает зависимость Р(у,Р) от обоих параметров. Наибольшие отклонения выражения (2.73) от численных результатов наблюдаются на краях при Р и 0 и Р и Р. Это, однако, не имеет большого значения, поскольку случай Р = 0 описывается кривыми, показанными на рис. 3; в области Р< 0,34 полученные данные, как это следует из рис. 12, хорошо описываются зависимостью (2.73), а узкий интервал 0,34 < Р < Р и 0,346 не представляет большого интереса, поскольку находится близко к значению Р, при котором тело «прилипает» к поверхности. Кроме того, при необходимости всегда можно провести дополнительное моделирование согласно процедуре, подробно описанной в подразд. 2.2.2, 2.3.2 и 2.4.2.

Полученные зависимости (2.69) (или точное выражение (2.71)) и (2.73) позволяют вычислить все параметры столкновения в актуальной области значений у (2.29). Например, найдем правило вычисления углов скольжения Р° как функции угла скольжения а° (см. рис. 1). Согласно рис. 1, имеем tga° = vz0 / vj0 и tgP° =| vz | /vx. Из выражений (2.19), (2.56) получим

tgP° =

7F(P)

2

2Р(у,р) + 5 R«c [ Р(у,Р) +1]

2tg a° vz0

Y1

(2.74)

где Е(Р) определяется согласно (2.70) по известной зависимости (2.69) (либо аналогичной точной (2.71)). Из (2.74) следует, что при ю0 Ф 0 угол скольжения Р° зависит явным образом от начальной нормальной компоненты скорости . При ю0 = 0 выражение (2.74) упрощается и принимает вид

tgP° =

7F (p)tg a°

2Р(у,Р) + 5.

(2.75)

Заметим, что в случае отсутствия адгезии между контактирующими телами (Ду = 0) в уравнениях (2.74), (2.75)

следует положить Е(Р) = 1. Формула (2.75) может использоваться для объяснения серии экспериментов, проведенных в работе [38], в которой изучалось падение сферического тела под углом к отражающей поверхности при нулевой начальной частоте вращения ю0 = 0 .

2.4.4. Сравнение с экспериментом

Насколько нам известно, на сегодняшний день не существует исчерпывающих аналитических, численных либо экспериментальных результатов, которые описывали бы столкновение под углом при наличии адгезионных сил между контактирующими поверхностями. К примеру, в работе [52] проводится моделирование методом конечных элементов столкновения под углом с учетом адгезии, однако рассматривается упругоплас-тический удар, нас же интересует исключительно упругая область. Молекулярно-динамические модели хотя и учитывают адгезию естественным образом, но не в состоянии описать поведение систем на макроскопическом уровне [55]. Тем не менее существует большое количество работ, которые описывают нормальные столкновения частиц с учетом адгезии. На рис. 2 и 3 в первой части обзора [8] приведены экспериментальные результаты из этих работ, также в [8] кратко описан теоретический подход к этому вопросу. Поэтому ограничимся сравнением полученных нами результатов с имеющимися данными по нормальному столкновению частиц. В работе [57] описаны результаты эксперимента по отражению частицы, выполненной из полистирола, от гладкой поверхности кварца. На рис. 13 для сравнения с нашей теорией показаны эти же результаты, но в области упругого удара, для которого характерным является возрастание коэффициента восстановления е с ростом нормальной скорости vz0.

Фактически, сплошной линией на рис. 13 показана функция Е(Р) (2.56), которая определяется из аппроксимации (2.69) при использовании (2.70). Параметр Р (2.58) содержит информацию о параметрах системы,

которые представлены в подписях к рис. 13. Эффективный модуль упругости Е * определяется по формуле (2.17), где для полистирола Е = 3 ■ 109 Па, 0,34 , а для кварца Е2 = 7,3 ■ 1010 Па, у2 = 0,17.

Рис. 13. Сравнение экспериментальных результатов, полученных в работе [57] (показаны символами), с численными результатами (сплошная линия) при параметрах R = 0,635 -10-6 м, р = 1000 кг/м3, Е*= 3,25 -109 Па, Ду = 0,15 Дж/м2. Данные, показанные кривой, нормированы на величину 0,96 для исключения влияния неупругих эффектов в материале,

которые не учитываются в модели Fig. 13. Comparison of the experimental results obtained in [57] (shown by symbols) with numerical results (solid line) at the parameters R = 0.635 -10-6 m, р = 1000 kg/m3, E* = 3.25 -109 Pa, Ду = 0.15 J/m2. The data shown by the curve is normalized to the value 0.96 to exclude the effect of inelastic effects in the material that are not taken into account in the model

Отметим, что согласно теоретическим представлениям о нормальном столкновении с адгезией (см. формулу (3.24) и рис. 5 в [8] коэффициент восстановления е при нормальном столкновении в упругой области с ростом начальной скорости v должен асимптотически приближаться к единице, чего не происходит в эксперименте, результаты которого показаны на рис. 13. Это связано с тем, что при уже сравнительно малых скоростях в соударяющихся с поверхностью частицах начинаются неупругие эффекты, связанные с пластической деформацией и вязкой диссипацией. Для исключения этих эффектов, которые не учитываются в нашей модели, полученная в результате моделирования зависимость e(vz0) нормирована на максимальное значение 0,96. Из рис. 13 следует, что при исключении неупругих эффектов результаты моделирования с высокой точностью совпадают с результатами реального эксперимента. Заключение

В первой части обзора [8] проведен анализ теоретических и экспериментальных работ различных авторов по физике столкновений упругих сферических тел с упругой плоскостью. Показано, что в зависимости от наличия сил трения, адгезии, вязких и пластических эффектов, шероховатостей поверхностей и т.п. возможно принципиально различное поведение. В представленной (второй) части работы показаны результаты, полученные авторами настоящего обзора. Так, рассмотрены вопросы упругого удара при наличии трения между поверхностями контактирующих тел, а также случай их адгезионного взаимодействия. Отдельно выделена ситуация отсутствия проскальзывания между контактирующими телами, что эквивалентно бесконечному коэффициенту трения ц . Во всех случаях определена функция P, которая является функцией безразмерных параметров, включающих в себя тангенциальную и нормальную компоненты начальных скоростей, начальную циклическую частоту вращения сферического тела, его массу, модули упругости и сдвига контактирующих тел, коэффициент трения в случае его конечного значения, а при учете адгезии также и поверхностную энергию. Таким образом, учтены основные параметры, задающие поведение частиц при столкновении. Модель, предложенная в работе, позволяет определить кинетические параметры после столкновения, такие как нормальная и тангенциальная компоненты скорости отражения частицы, а также циклическая частота ее вращения в зависимости от значений аналогичных параметров до столкновения. Все величины, измеряемые экспериментально, такие как нормальный и тангенциальный коэффициенты восстановления, угол отражения, потери кинетической энергии и т.п., выражаются через численно найденную функцию P. Кроме этого, подробно описана сама процедура моделирования, позволяющая находить не только параметры после столкновения, а и определять изменяющиеся во времени величины на протяжении контакта, такие как глубина индентирования сферы в упругое полупространство и ее горизонтальное смещение, нормальная и тангенциальная компоненты контактных сил, временные зависимости скоростей и частоты вращения частицы и т.д. Во всех рассмотренных случаях удар считается упругим, т.е. пластические эффекты в модели не учитывают-

Библиографический список

1. Dynamics of drag and force distributions for projectile impact in a granular medium / M.P. Ciamarra, A.H. Lara, A.T. Lee, D.I. Goldman, I. Vishik, H.L. Swinney // Phys. Rev. Lett. - 2004. - Vol. 92. - No. 19. -P. 194301.

2. Jop P., Forterre Y., Pouliquen O. A constitutive law for dense granular flows // Nature. - 2006. - Vol. 441. - No. 7094. - P. 727-730.

3. Model for collisions in granular gases / N.V. Brilliantov, F. Spahn, J.-M. Hertzsch, T. Pöschel // Phys. Rev. E. - 1996. -Vol. 53. - No. 5. - P. 5382-5392.

ся. Проведено сравнение с экспериментальными результатами. Так, в случае конечного коэффициента трения для упругого удара, когда нормальный коэффициент восстановления равен единице, теоретически рассчитаны зависимости тангенциальной компоненты коэффициента восстановления и циклической частоты вращения частицы после столкновения. Сравнение указанных зависимостей с полученными экспериментально показывает хорошее согласие. При учете адгезии проведено сравнение экспериментально определенного коэффициента восстановления при нормальном столкновении с теоретической кривой, что также показывает совпадение с высокой точностью. Поскольку теория учитывает весь необходимый набор материальных параметров контактирующих тел и начальных условий, она может использоваться для объяснения и предсказания результатов экспериментов по упругому удару однородных сферических тел под углом к поверхности, когда ненулевые значения принимают как нормальная, так и тангенциальная компоненты скоростей для описанных ситуаций, т.е. столкновение частиц без проскальзывания в зоне контакта, столкновение с конечной силой трения и столкновение при учете адгезионного взаимодействия между контактирующими частицами. Отметим, что задача чисто нормального соударения является гораздо более простой, чем косой удар, и она была решена нами полуаналитически в недавней работе [56] в полном диапазоне изменений параметра Тейбора, задающего адгезионные свойства контакта (см. также работу Мауги-са [49]).

Благодарности

Публикация содержит результаты исследований, проведенных при финансовой поддержке Немецкого научно-исследовательского сообщества (DFG), проект PO 810-55-1.

Acknowledgements

The publication contains the results of researchs conducted with the financial support of the German Research Community (DFG), Project ID - PO 810-55-1.

4. Bernard B. Impacts in mechanical systems: analysis and modelling. - Berlin, New York: Springer, 2000. - 278 p.

5. Attractive particle interaction forces and packing density of fine glass powders / E.J.R. Parteli, J. Schmidt, C. Blumel, K.-E. Wirth, W. Peukert, T. Poschel // Sci. Rep. - 2014. - Vol. 4. - P. 6227 (7 pp.).

6. Родионов А.И., Матвеев К.А. К динамике удара абсолютно твердого шара по упругому полупространству // Научный вестник НГТУ. - 2012. - № 1 (46). - С. 93-108.

7. Maw N., Barber J.R., Fawcett J.N. The oblique impact of elastic spheres // Wear. - 1976. - Vol. 38. - No. 1. - P. 101-114.

8. Ляшенко Я.А., Виллерт Э., Попов В.Л. Механика столкновений твердых тел: влияние трения и адгезии. I. Обзор экспериментальных и теоретических работ // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2018. - № 2. - С. 44-61. DOI: 10.15593/perm.mech/2018.2.05

9. Popov V.L., Heß M. Method of dimensionality reduction in contact mechanics and friction. - Berlin: Springer, 2014. - 265 p.

10. Jia Y.-B. Three-dimensional impact: energy-based modeling of tangential compliance // The International Journal of Robotics Research. - 2013. - Vol. 32. - No. 1. - P. 56-83.

11. Popov V.L., Psakhie S.G., Numerical simulation methods in tribology: possibilities and limitations // Tribology International. - 2007. - Vol. 40. - No. 6. - P. 916-923.

12. Popov V.L. Method of reduction of dimensionality in contact and friction mechanics: A linkage between micro and macro scales // Friction. - 2013. - Vol. 1. - No. 1. - P. 41-62.

13. Popov V.L., Heß M. Methode der Dimensionsreduktion in Kontaktmechanik und Reibung. Eine Berechnungsmethode im Mikro- und Makrobereich. - Berlin: Springer, 2013. - 267 p.

14. Popov V.L., Hess M. Method of dimensionality reduction in contact mechanics and friction: a users handbook. I. Axial-ly-symmetric contacts // Facta Universitatis, series Mechanical Engineering. - 2014. - Vol. 12. - No. 1. - P. 1-14.

15. Wetter R., Popov V.L. Shakedown limits for an oscillating, elastic rolling contact with Coulomb friction // Int. J. of Solids and Structures. - 2014. - Vol. 51. - No. 5. - P. 930-935.

16. Клинков С.В., Косарев В.Ф. Моделирование адгезионного взаимодействия частиц с преградой при газодинамическом напылении // Физ. мезомех. - 2002. - Т. 5, № 3. - С. 27-35.

17. Stronge W. J. Impact Mechanics. - Cambridge: Cambridge University Press, 2004. - 280 p.

18. Раус Э.Дж. Динамика системы твердых тел. Т. 1. -М.: Наука, 1983. - 464 с.

19. Hauger W., Schnell W., Gross D. Technische Mechanik. Bd 3: Kinetik 7. - Berlin: Springer, 2002. - 267 p.

20. Lyashenko I.A., Popov V.L. Impact of an elastic sphere with an elastic half space revisited: Numerical analysis based on the method of dimensionality reduction // Sci. Rep. - 2015. -Vol. 5. - P. 8479 (5 p.).

21. Lyashenko I.A., Willert E., Popov V.L. Adhesive impact of an elastic sphere with an elastic half space: Numerical analysis based on the method of dimensionality reduction // Mech. Mater. -2016. - Vol. 92. - P. 155-163.

22. Willert E., Popov V.L. Impact of an elastic sphere with an elastic half space with a constant coefficient of friction: Numerical analysis based on the method of dimensionality reduction // ZAMM - Journal of Applied Mathematics and Mechanics / Zeitschrift für Angewandte Mathematik und Mechanik. - 2016. - Vol. 96. - P. 108-91095.

23. Popov V.L. Contact mechanics and friction. Physical principles and applications. 2nd Edition. - Berlin: Springer, 2017. - 362 p.

24. Influence of complete sticking on the indentation of a rigid cone into an elastic half-space in the presence of molecular adhesion / F.M. Borodich, B.A. Galanov, Yu.I. Prostov, M.M. Suarez-Alvarez // J. Appl. Math. Mech. - 2012. - Vol. 76. - No. 5. - P. 590-596.

25. Argatov I.I. Asymptotic modeling of the impact of a spherical indenter on an elastic half-space // Int. J. Solids Struct. -2008. - Vol. 45. - No. 18-19. - P. 5035-5048.

26. Hunter S.C. Energy absorbed by elastic waves during impact // J. Mech. Phys. Solids. - 1957. - Vol. 5. - No. 3. - P. 162-171.

27. Broomhead D.S., Gutkin E. The dynamics of billiards with no-slip collisions // Physica D. - 1993. - Vol. 67. - No. 1-3. -P. 188-197.

28. Mindlin R.D. Compliance of elastic bodies in contact // ASME J. Applied Mech. - 1949. - Vol. 16. - P. 259-268.

29. Campana C., Persson B.N.J., Muser M.H. Transverse and normal interfacial stiffness of solids with randomly rough surfaces // J. of Phys.: Condens. Matter. - 2011. - Vol. 23. -No. 8. - P. 085001 (9 pp.).

30. Grzemba B., Pohrt R., Teidelt E. Maximum micro-slip in tangential contact of randomly rough self-affine surfaces // Wear. - 2014. - Vol. 309. - No. 1-2. - P. 256-258.

31. Persson B.N.J. On the fractal dimension of rough surfaces // Tribol. Lett. - 2014. - Vol. 54. - No. 1. - Р. 99-106.

32. Pogrebnjak A.D., Borisyuk V.N., Bagdasaryan A.A. Numerical analysis of the morphological and phase changes in the TiN/Al2O3 coating under high current electron beam modification // Cond. Matt. Phys. - 2013. - Vol. 16. - No. 3. - P. 33803: 1-8.

33. Lyashenko I.A., Borysiuk V.N., Manko N.N. Statistical analysis of self-similar behaviour in the shear induced melting model // Cond. Matt. Phys. - 2014. - Vol. 17. - No. 2. - Р. 23003: 1-11.

34. Khomenko A.V., Lyashenko I.A., Borisyuk, V.N. Mul-tifractal analysis of stress time series during ultrathin lubricant film melting // Fluct. Noise Lett. - 2010. - Vol. 9. - No. 1. - P. 19-35.

35. Cattaneo C. Sul contatto di due corpi elastici: distribuzione locale degli sforzi // Rendiconti dell'Accademia nazionale dei Lincei. -1938. - Vol. 27. - No. 6. - P. 342-348; 434-436; 474-478.

36. Munisamy R.L., Hills D.A., Nowell D. Static axisym-metric hertzian contacts subject to shearing forces // ASME J. Appl. Mech. - 1994. - Vol. 61. - No. 2. - P. 278-283.

37. Ландау Л.Д., Лифшиц Е.М. Теория упругости. - М.: Наука, 1987. - 246 с.

38. Gorham D.A., Kharaz A.H. The measurement of particle rebound characteristics // Powder Technol. - 2000. - Vol. 112. -No. 3. - P. 193-202.

39. Ning Z. Elasto-plastic impact of fine particles and fragmentation of small agglomerates: PhD Thesis. - Birmingham, Aston University, 1995. - 269 p.

40. Sondergaard R., Chaney K., Brennen C.E. Measurements of solid spheres bouncing off flat plates // ASME J. Appl. Mech. - 1990. - Vol. 57. - No. 3. - P. 694-699.

41. Великанов Н.Л., Наумов В.А., Корягин С.И. Упругий удар тела о наклонную поверхность // Вестник Балтийского федерального университета им. И. Канта. -2013. - Вып. 10. - С. 36-42.

42. Лапшин В.Л., Рудых А.В., Глухов А.В. Разработка ударного способа определения модуля упругости материала // Системы. Методы. Технологии. - 2015. - № 2(26). - С. 37-43.

43. Waters J.F., Guduru P.R. Mode-mixity-dependent adhesive contact of a sphere on a plane surface // Proc. R. Soc. A. -2010. - Vol. 466. - P. 1303-1325.

44. Kosinski P., Hoffmann A.C. Extended hard-sphere model and collisions of cohesive particles // Phys. Rev. E. - 2011. -Vol. 84. - No. 3. - P. 031303 (13 pp.).

45. Johnson K.L., Kendall K., Roberts A.D. Surface energy and the contact of elastic solids // Proc. R. Soc. London A. - 1971. -Vol. 324. - No. 1558. - P. 301-313.

46. Lednikova B., Kosek J. The effects of adhesive forces on the collision behavior of polyolefin particles // Powder Technol. -2013. - Vol. 243. - P. 27-39.

47. Andres R.P. Inelastic energy transfer in particle surface collisions // Aerosol Sci. Technol. - 1995. - Vol. 23. - No. 1. - P. 40-50.

48. Derjaguin B.V., Muller V.M., Toporov Yu.P. Effect of contact deformations on the adhesion of particles // J. Colloid In-terf. Sci. - 1975. - Vol. 53. - No. 2. - P. 314-326.

49. Maugis D. Adhesion of spheres: the JKR-DMT transition using a dugdale model // J. Colloid Interf. Sci. - 1992. -Vol. 150. - No. 1. - P. 243-269.

50. Collision dynamics of granular particles with adhesion / N.V. Brilliantov, N. Albers, F. Spahn, T. Pöschel // Phys. Rev. E. -2007. - Vol. 76. - No. 5. - P. 051302 (12 pp.).

51. Negative normal restitution coefficient found in simulation of nanocluster collisions / K. Saitoh, A. Bodrova, H. Hayaka-wa, N.V. Brilliantov // Phys. Rev. Lett. - 2010. - Vol. 105. -No. 23. - P. 238001 (4 pp.).

52. Katta R.R., Polycarpou A.A., Hanchi J.V. High velocity oblique impact and coefficient of restitution for head disk interface operational shock // ASME J. Tribol. - 2009. - Vol. 131. -No. 2. - P. 021903 (9 pp.).

53. Hertz H.J. Ueber die berührung fester elastischer körper // J. für die Reine Angew. Math. - 1882. - Vol. 92. - P. 156-171.

References

1. Ciamarra M.P, Lara A.H., Lee A.T., Goldman D.I, Vishik I., Swinney H.L. Dynamics of drag and force distributions for projectile impact in a granular medium. Phys. Rev. Lett, 2004, vol. 92, no. 19, pp. 194301.

2. Jop P., Forterre Y., Pouliquen O. A constitutive law for dense granular flows. Nature, 2006, vol. 441, no. 7094, pp. 727-730.

3. Brilliantov N.V., Spahn F., Hertzsch J.-M., Pöschel T. Model for collisions in granular gases. Phys. Rev. E, 1996, vol. 53, no. 5, pp. 5382-5392.

4. Bernard B. Impacts in mechanical systems: analysis and modelling. Berlin, New York: Springer, 2005, 278 p.

5. Parteli E.J.R., Schmidt J., Blümel C., Wiith K.-E., Peukert W., Pöschel T. Attractive particle interaction forces and packing density of fine glass powders. Sci. Rep, 2014, vol. 4, p. 6227 (7 pp.).

6. Rodionov A.I., Matveev K.A. K dinamike udara absolyutno tverdogo shara po uprugomu poluprostranstvu [To the dynamics of the impact of an absolutely rigid ball on an elastic half-space]. Scientific Bulletin of NSTU, 2012, no. 1 (46), pp. 93-108.

7. Maw N., Barber J.R., Fawcett J.N. The oblique impact of elastic spheres. Wear, 1976, vol. 38, no. 1, pp. 101-114.

8. Lyashenko I.A., Willert E., Popov V.L. Mechanics of collisions of solids: influence of friction and adhesion. I. Review of experimental and theoretical works. PNRPU Mechanics Bulletin, 2018, no. 2, pp. 44-61. DOI: 10.15593/perm.mech/2018.2.05

9. Popov V.L., Heß M. Method of dimensionality reduction in contact mechanics and friction. Berlin, Springer, 2014, 265 p.

10. Jia Y.-B. Three-dimensional impact: energy-based modeling of tangential compliance. The International Journal of Robotics Research, 2013, vol. 32, no. 1, pp. 56-83.

11. Popov V.L., Psakhie S.G., Numerical simulation methods in tribology: possibilities and limitations. Tribology International, 2007, vol. 40, no. 6, pp. 916-923.

12. Popov V.L. Method of reduction of dimensionality in contact and friction mechanics: A linkage between micro and macro scales. Friction, 2013, vol. 1, no. 1, pp. 41-62.

13. Popov V.L., Heß M. Methode der Dimensionsreduktion in Kontaktmechanik und Reibung. Eine Berechnungsmethode im Mikro- und Makrobereich, Berlin: Springer, 2013, 267 p.

14. Popov V.L., Hess M. Method of dimensionality reduction in contact mechanics and friction: a users handbook. I. Axi-ally-symmetric contacts. Facta Universitatis, series Mechanical Engineering, 2014, vol. 12, no. 1, pp. 1-14.

15. Wetter R., Popov V.L. Shakedown limits for an oscillating, elastic rolling contact with Coulomb friction. Int. J. of Solids and Structures, 2014, vol. 51, no. 5, pp. 930-935.

16. Klinkov S.V., Kosarev V.F. Modeling of particle-substrate adhesive interaction in cold spraying. Fizicheskaja Mesomechanika, 2002, vol. 5, no. 3, pp. 27-35.

54. Brach R.M., Dunn P.F. A mathematical model of the impact and adhesion of microsphers // Aerosol Sci. Technol. -1992. - Vol. 16. - No. 1. - P. 51-64.

55. Thornton C., Ning Z. A theoretical model for the stick/bounce behaviour of adhesive, elastic-plastic spheres // Powder Technol. - 1998. - Vol. 99. - No. 2. - P. 154-162.

56. Willert E., Lyashenko I.A., Popov V.L. Influence of the Tabor parameter on the adhesive normal impact of spheres in Maugis-Dugdale approximation // Comp. Part. Mech. - 2018. -Vol. 5. - No. 3. - P. 313-318.

57. Dahneke B.E. Particle bounce or capture-search for an adequate theory: I. Conservation-of-energy model for a simple collision process // Aerosol Sci. Technol. - 1995. - Vol. 23. -No. 1. - P. 25-39.

17. Stronge W. J. Impact Mechanics. Cambridge, Cambridge University Press, 2004, 280 p.

18. Routh E.J. Dynamics of a system of rigid bodies. Part 1. London, Macmillan and Co. 1955, 448 p.

19. Hauger W., Schnell W., Gross D. Technische Mechanic. Bd 3: Kinetik 7, Berlin: Springer, 2002, 267 p.

20. Lyashenko I.A., Popov V.L. Impact of an elastic sphere with an elastic half space revisited: Numerical analysis based on the method of dimensionality reduction. Sci. Rep., 2015, vol. 5, pp. 8479

(5 pp.).

21. Lyashenko I.A., Willert E., Popov V.L. Adhesive impact of an elastic sphere with an elastic half space: Numerical analysis based on the method of dimensionality reduction. Mech. Mater, 2016, vol. 92, pp. 155-163.

22. Willert E., Popov V.L. Impact of an elastic sphere with an elastic half space with a constant coefficient of friction: Numerical analysis based on the method of dimensionality reduction. ZAMM - Journal of Applied Mathematics and Mechanics, 2016, vol. 96, pp. 1089-1095.

23. Popov V.L. Contact mechanics and friction. Physical principles and applications, Berlin, Springer, 2010, 362 p.

24. Borodich F.M., Galanov B.A., Prostov Yu. I., Suarez-Alvarez M.M. Influence of complete sticking on the indentation of a rigid cone into an elastic half-space in the presence of molecular adhesion. J. Appl. Math. Mech., 2012, vol. 76, no. 5, pp. 590-596.

25. Argatov I.I. Asymptotic modeling of the impact of a spherical indenter on an elastic half-space. Int. J. Solids Struct., 2008, vol. 45, no. 18-19, pp. 5035-5048.

26. Hunter S.C. Energy absorbed by elastic waves during impact. J. Mech. Phys. Solids, 1957, vol. 5, no. 3, pp. 162-171.

27. Broomhead D.S., Gutkin E. The dynamics of billiards with no-slip collisions. Physica D, 1993, vol. 67, no. 1-3, pp. 188-197.

28. Mindlin R.D. Compliance of elastic bodies in contact. ASME J. Applied Mech., 1949, vol. 16, pp. 259-268.

29. Campanâ C., Persson B.N.J., Müser M.H. Transverse and normal interfacial stiffness of solids with randomly rough surfaces. J. of Phys.: Condens. Matter., 2011, vol. 23, no. 8, pp. 085001 (9 pp.).

30. Grzemba B., Pohrt R., Teidelt E. Maximum micro-slip in tangential contact of randomly rough self-affine surfaces. Wear, 2014, vol. 309, no. 1-2, pp. 256-258.

31. Persson B.N.J. On the fractal dimension of rough surfaces. Tribol. Lett, 2014, vol. 54, no. 1, pp. 99-106.

32. Pogrebnjak A.D., Borisyuk V.N., Bagdasaryan A.A. Numerical analysis of the morphological and phase changes in the TiN/Al2O3 coating under high current electron beam modification. Cond. Matt. Phys., 2013, vol. 16, no. 3, pp. 33803: 1-8.

33. Lyashenko I.A., Borysiuk V.N., Manko N.N. Statistical analysis of self-similar behaviour in the shear induced melting model. Cond. Matt. Phys, 2014, vol. 17, no. 2, pp. 23003: 1-11.

34. Khomenko A.V., Lyashenko I.A., Borisyuk, V.N. Mul-tifractal analysis of stress time series during ultrathin lubricant film melting. Fluct. Noise Lett, 2010, vol. 9, no. 1, pp. 19-35.

35. Cattaneo C. Sul contatto di due corpi elastici: distribuzione locale degli sforzi. Rendiconti dell'Accademia nazionale dei Lincei, 1938, vol. 27, no. 6, pp. 342-348; 434-436; 474-478.

36. Munisamy R.L., Hills D.A., Nowell D. Static axisym-metric hertzian contacts subject to shearing forces. ASME J. Appl. Mech., 1994, vol. 61, no. 2, pp. 278-283.

37. Landau L.D., Lifshitz E.M. Theory of elasticity. Oxford, Pergamon Press, 1986, 195 p.

38. Gorham D.A., Kharaz A.H. The measurement of particle rebound characteristics. Powder Technol., 2000, vol. 112, no. 3, pp. 193-202.

39. Ning Z. Elasto-plastic impact of fine particles and fragmentation of small agglomerates: PhD Thesis. Birmingham, Aston University, 1995, 269 p.

40. Sondergaard R., Chaney K., Brennen C.E. Measurements of solid spheres bouncing off flat plates. ASME J. Appl. Mech., 1990, vol. 57, no. 3, pp. 694-699.

41. Velikanov N.L., Naumov V.A., Koryagin S.I. Uprugiy udar tela o naklonnuyu poverhnost [Elastic impact of the body on an inclined surface]. Bulletin of the Baltic Federal University named I. Kant, 2013, no. 10, pp. 36-42.

42. Lapshin V.L., Rudyih A.V., Gluhov A.V. Razrabotka udarnogo sposoba opredeleniya modulya uprugosti materiala [Development of an impact method for determining the modulus of elasticity of a material]. Sytemi. Metodi. Tekhnologii, 2015, no. 2(26), pp. 37-43.

43. Waters J.F., Guduru P.R. Mode-mixity-dependent adhesive contact of a sphere on a plane surface. Proc. R. Soc. A, 2010, vol. 466, pp. 1303-1325.

44. Kosinski P., Hoffmann A.C. Extended hard-sphere model and collisions of cohesive particles. Phys. Rev. E, 2011, vol. 84, no. 3, pp. 031303 (13 pp.).

45. Johnson K.L., Kendall K., Roberts A.D. Surface energy and the contact of elastic solids. Proc. R. Soc. London A, 1971, vol. 324, no. 1558, pp. 301-313.

46. Lednikova B., Kosek J. The effects of adhesive forces on the collision behavior of polyolefin particles. Powder Technol, 2013, vol. 243, pp. 27-39.

47. Andres R.P. Inelastic energy transfer in particle surface collisions. AerosolSci. Technol., 1995, vol. 23, no. 1, pp. 40-50.

48. Derjaguin B.V., Muller V.M., Toporov Yu.P. Effect of contact deformations on the adhesion of particles. J. Colloid Interf. Sci, 1975, vol. 53, no. 2, pp. 314-326.

49. Maugis D. Adhesion of spheres: the JKR-DMT transition using a dugdale model. J. Colloid Interf. Sci., 1992, vol. 150, no. 1, pp. 243-269.

50. Brilhantov N.V., Albers N., Spahn F., Poschel T. Collision dynamics of granular particles with adhesion. Phys. Rev. E, 2007, vol. 76, no. 5, pp. 051302 (12 pp.).

51. Saitoh K., Bodrova A., Hayakawa H., Brilliantov N.V. Negative normal restitution coefficient found in simulation of nanocluster collisions. Phys. Rev. Lett., 2010, vol. 105, no. 23, pp. 238001 (4 pp.).

52. Katta R.R., Polycarpou A.A., Hanchi J.V. High velocity oblique impact and coefficient of restitution for head disk interface operational shock. ASME J. Tribol., 2009, vol. 131, no. 2, pp. 021903 (9 pp.).

53. Hertz H. J. Ueber die beruhrung fester elastischer korper. J. fur die reine Angew. Math., 1882, vol. 92, pp. 156-171.

54. Brach R.M., Dunn P.F. A mathematical model of the impact and adhesion of microsphere. Aerosol Sci. Technol., 1992, vol. 16, no. 1, pp. 51-64.

55. Thornton C., Ning Z. A theoretical model for the stick/bounce behaviour of adhesive, elastic-plastic spheres. Powder Technol., 1998, vol. 99, no. 2, pp. 154-162.

56. Willert E., Lyashenko I.A., Popov V.L. Influence of the Tabor parameter on the adhesive normal impact of spheres in Maugis-Dugdale approximation. Comp. Part. Mech., 2018, vol. 5, no. 3, pp. 313-318.

57. Dahneke B.E. Particle bounce or capture-search for an adequate theory: I. Conservation-of-energy model for a simple collision process. Aerosol Sci. Technol., 1995, vol. 23, no. 1, pp. 25-39.