Changes in vehicle's movement surface profile due to snow cover, weathering, melting, deformation of the supporting base, wheel smoothing ability and other processes

全文:

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

В статье показано приложение аналитических методов, подробно изложенных в [1-3], к исследованию изменения профиля поверхности движения, которое может быть вызвано различными причинами: ростом снежного покрова, деформацией под нагрузкой, износом, выветриванием и др. Кроме того, приведенная модель позволяет описывать траекторию центра колеса, движущегося по неровному профилю дороги при длинах неровностей, сопоставимых с радиусом колеса. Последняя задача является наименее разработанной ввиду сложности ее формализации и описания: при движении колеса по неровному профилю возникают положения, когда колесо имеет двухточечный контакт (перекатывается через яму), что приводит к потере гладкости решений. Наиболее близкими исследованиями в рамках данного аспекта изучения поверхности движения являются работы Барахтанова Л.В., Аникина А.А. [4], Кулешова А.П., Молева Ю.И., Жилейкина М.М. [5]. В работах Барахтанова Л.В. и Аникина А.А. применяется вероятностный подход к изучению заснеженных поверхностей движения; в работе Жилейкина М.М. используется кусочно-линейная аппроксимация поверхности движения. Предлагаемая в настоящей работе математическая модель носит более общий характер по сравнению с перечисленными исследованиями, а также позволяет учитывать такие эффекты, как скорость роста поверхности, неравномерность роста поверхности (в зависимости от направления нормали); изменение профиля в ветреную погоду со снегом, когда наблюдается продольное движение профиля и т.п.

Будем считать, что в начальный момент t₀ ординаты профиля дороги описываются (начальной) функцией:

$z=z\left(x,t_0\right)=z_0$.

С течением времени  происходит эволюция профиля дороги, причины которого на данном этапе не требуют конкретизации (рост снежного покрова; деформация основания; сглаживающая способность шины и пр.). Предполагается только, что изменение профиля дороги происходит перпендикулярно начальному профилю со скоростью . Причем скорость  может быть как положительная, так и отрицательная – в зависимости от физической постановки задачи (рост поверхности или ее спад). Рассмотрим произвольную точку (т.P) на начальной поверхности с координатами {x_P, z_P}, ее эволюция будет происходить по траектории, перпендикулярной – в каждый момент времени t – к текущему профилю дороги z = z(x, t) (рис. 1).

 

Рис. 1. Кинематика роста профиля поверхности

Fig. 1. Kinematics of surface profile growth

 

Будем обозначать траекторию эволюции точки Р функцией вида:

$z_p=\xi \left(x_p,t\right)$. (1)

Назовем траекторию точки «ростком» (рис. 2). Координаты точки, принадлежащей «ростку», должны удовлетворять уравнениям:

$\left.\begin{array}{c}\frac{d{x}_p}{dt}=w\times \sin \phi \\ \frac{d{z}_p}{dt}=w\times \cos \phi \end{array}\right\}$. (2)

 

Рис. 2. «Росток» в точке P

Fig. 2. «Sprout» at point P

 

При рассмотрении эволюции координат точки Р {x_P(t), z_P(t)} должно выполняться также условие:

$z_p\left(t\right)=z\left(x_p\left(t\right),t\right)=\xi \left(x_p,t\right)$. (3)

Тогда второе уравнение системы примет вид:

$\frac{\partial z}{\partial t}+\frac{\partial z}{\partial x}\frac{d{x}_p}{dt}=w\times \cos \phi$. (4)

Подставляя в последнее уравнение (4) первое уравнение системы (2), получим:

$\frac{\partial z}{\partial t}+\frac{\partial z}{\partial x}w\times \sin \phi =w\times \cos \phi$. (5)

Учтем еще две зависимости из элементарной тригонометрии, которые позволят последнее уравнение выразить от одной неизвестной переменной (замкнуть):

$\frac{\partial z}{\partial x}=-tg\phi ; \cos \phi =\frac{1}{\sqrt{1+{\left(tg\phi \right)}^2}}$. (6)

Таким образом, получаем уравнение в виде

$$\begin{gathered} \frac{\partial z}{\partial t}-tg\phi \times w\times \sin \phi =w\times \cos \phi ; \\ \cos \phi \frac{\partial z}{\partial t}=w{\left(\sin \phi \right)}^2+w{\left(\cos \phi \right)}^2=w; \end{gathered}$$

$\frac{\partial z}{\partial t}=\frac{w}{\cos \phi }$; (7)

$\frac{\partial z}{\partial t}=w\sqrt{1+{\left(\frac{\partial z}{\partial x}\right)}^2}$. (8)

Рассмотрим ограничение уравнения (8) на вид начального профиля опорной поверхности, т.е. будем считать, что углы φ являются малыми. В этом случае справедливы приближенные равенства:

$$\begin{gathered} \frac{\partial z}{\partial x}=-tg\phi \approx -\phi ; \\ \sqrt{1+{\left(tg\phi \right)}^2}\approx 1+\frac{1}{2}{\left(\frac{\partial z}{\partial x}\right)}^2+\dots \end{gathered}$$

Тогда уравнение (8) приведется к виду:

$\frac{\partial z}{\partial t}=w+\frac{w}{2}{\left(\frac{\partial z}{\partial x}\right)}^2$. (9)

Из уравнения видно, что одно из решений, которому это уравнение удовлетворяет, является решение вида:

$z\left(x,t\right)=wt$. (10)

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

$z^{\text{'}}\left(x,t\right)=z\left(x,t\right)-wt$.

Тогда уравнение (10) примет вид:

$\frac{\partial z^{\text{'}}}{\partial t}=\frac{w}{2}{\left(\frac{\partial z^{\text{'}}}{\partial x}\right)}^2$.

Вернемся к общему уравнению (8) и обобщим его на случай анизотропного роста, т.е. учтем зависимость скорости роста w от направления нормали. Введем новую функцию $u\left(x,t\right)=\frac{\partial z}{\partial x}$, смысл которой – тангенс угла наклона ростка. Тогда уравнение (8) представится в виде:

$$\begin{gathered} \frac{\partial z}{\partial t}=w\left(u\right)\sqrt{1+u{\left(x,t\right)}^2}, или \\ \frac{\partial z}{\partial t}=\Phi \left(u\right). \end{gathered}$$ (11)

Продифференцируем (11) по x, получим уравнение для поля наклонов u(x, t):

$\frac{\partial u}{\partial t}+C\left(u\right)\frac{\partial u}{\partial x}=0$, (12)

где $C\left(u\right)=\frac{d\Phi \left(u\right)}{du}=\frac{d}{du}\left(w\left(u\right)\sqrt{1+u^2}\right)$.

Последнее уравнение (12) является уравнением Римана или уравнением простой волны, решения которого методом характеристик приводятся в [2, 3].

Пример 1. Рассмотрим случай, когда поверхность растет с разной скоростью в различных направлениях. Анизотропию роста учтем, полагая скорость w зависящей от u.

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

$w\left(u\right)=w\cos \phi =\frac{w}{\sqrt{1+u^2}}$.

Подставив полученное выражение в уравнение (9), получим:

$\Phi \left(u\right)=w\left(u\right)\sqrt{1+u^2}=w$;

откуда следует, что уравнение роста (9) имеет тривиальное решение:

$z\left(x,t\right)=z_0+wt$,

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

Пример 2. Рассмотрим более общий случай, когда на поверхность движения осаждается частицы снега, падающие с разных направлений. Пусть функция D(φ) – диаграмма направленности, задающая интенсивность потока падающих частиц снега под углом φ к вертикальной оси z. Тогда полная интенсивность частиц, падающих на площадку поверхности, наклоненную под углом φ к оси z определяется выражением

$w\left(\phi \right)={\int }_{{\phi }_{-}}^{{\phi }_{+}}D\left(\tau \right)\cos \left(\phi -\tau \right)d\tau$. (13)

Первый пример, соответствующий вертикально падающему снегу, является частным случаем (13) с сингулярной диаграммой направленности:

$D\left(\phi \right)=w\delta \left(\phi \right)$.

Пределы интегрирования в (13) учитывают эффекты затенения рассматриваемого элемента профиля опорной поверхности z = z(x, t) остальной ее частью, т.е. являются ее функциями:

${\phi }_{+}={\phi }_{+}\left(z\right), {\phi }_{-}={\phi }_{-}\left(z\right)$.

Иллюстрация, поясняющая зависимость пределов интегрирования от вида опорной поверхности в окрестности рассматриваемой точки, показана на рис. 3.

 

Рис. 3. Зависимость пределов интегрирования диаграммы направленности осаждаемых частиц снега на опорную поверхность

Fig. 3. Dependence of the integration limits of the directional diagram of snow particles deposited on the supporting surface

 

На рис. 3 для определенности показана область рассматриваемой точки профиля поверхности, с углом нормали к вертикальной линии φ > 0. При этом пределы интегрирования определяются из условия возможности попадания частиц снега на поверхность: слева – угол падения ограничен касательной к поверхности; справа – горизонтальной прямой, проходящей через рассматриваемую точку.

В рамках данного примера рассмотрим изотропный поток D(φ) = const и отсутствие затенения (φ₊ = π / 2, φ_– = φ – π / 2):

$w\left(\phi \right)={\int }_{\phi -\pi /2}^{\pi /2}D\cos \left(\phi -\tau \right)d\tau =D\left(1+\cos \phi \right)$.

Пример 3. Рассмотрим случай, когда снегопад сопровождается ветром, а снег падает под углом φ₀ ≠ 0 к зениту. Начальный профиль дороги z = z(x, t₀) = z₀ имеет ограниченные углы наклона нормалей:

$\left|\phi -{\phi }_0\right|<\frac{\pi }{2}$.

Определим эволюцию профиля заснеженной поверхности движения. В этом случае диаграмма направленности имеет вид:

$D\left(\tau \right)=w\delta \left(\tau -{\phi }_0\right)$.

Ограничение на углы нормалей гарантирует, что для любой точки исходного профиля z(x, t₀) δ-функция не выйдет за пределы интегрирования. Следовательно, скорость прилипания снега на участок поверхности дороги, нормаль к которому составляет угол φ с вертикалью, определяется выражением:

$w=w\cos \left(\phi -{\phi }_0\right)=w\left(\cos \phi \cos {\phi }_0+\sin \phi \sin {\phi }_0\right)$.

Подставив это выражение в (13), перейдя от φ к u и введя обозначения w₁ = cosφ₀; w₂ = sinφ₀, запишем уравнение в терминах z(x, t):

$\frac{\partial z}{\partial t}=w_1+w_{2}u или \frac{\partial z}{\partial t}-w_1\frac{\partial z}{\partial x}=w_2, z\left(x,t_0\right)=z_0$.

Последнее уравнение является линейным и имеет решение в виде:

$z\left(x,t\right)=z_0\left(x-w_{2}t\right)+w_{1}t$.

Вид решения показывает, что профиль заснеженного основания не меняет формы, а смещается в сторону ветра (рис. 4).

 

Рис. 4. Продольное смещение профиля опорной поверхности при снегопаде, сопровождающемся ветром

Fig. 4. Longitudinal displacement of the supporting surface profile during snowfall accompanied by wind

 

Пример 4. Рассмотрим случай роста с изотропной скоростью, т.е. скоростью не зависящей от угла φ наклона нормали w(u) = w = const. В этом случае $C\left(u\right)=\frac{wu}{\sqrt{1+u^2}}=w\sin \phi$, ростки являются прямыми и скорость роста поверхности вдоль прямого ростка является постоянной. Рост профиля поверхности при этом задается параметрически [1]:

$x=y+\frac{u_0\left(y\right)}{\sqrt{1+u_0^2\left(y\right)}}wt, z=z_0\left(y\right)+\frac{wt}{\sqrt{1+u_0^2\left(y\right)}}$. (14)

Если рассмотреть в качестве исходного профиля кривой z = z(x, t) синусоиду, получим:

$z_0\left(x\right)=z_0\cos kx; u_0\left(x\right)=z_{0}k\sin kx$.

Тогда параметрическое уравнение линии z = z(x, t) имеет вид:

$z=\mu +\frac{\epsilon \tau \sin \mu }{\sqrt{1+{\epsilon }^2{\sin }^2\mu }}, \eta =\epsilon \cos \mu +\frac{\tau }{\sqrt{1+{\epsilon }^2{\sin }^2\mu }}$,

где введены безразмерные переменные kx = z; ky = μ; ckt = τ; kz = η и безразмерный параметр kz₀ = ε.

Профиль опорного основания z = z(x, t), построенный по приведенным выше параметрическим формулам (при ε = 0.5, для различных t), показан на рис. 5.

 

Рис. 5. Рост профиля поверхности z = z(x, t)

Fig. 5. Surface profile growth z = z(x, t)

 

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

Пример 5. Приведенные в предыдущем примере параметрические зависимости, описывающие рост поверхности, применимы и для описания выветриваемых или тающих поверхностей движения, если применить отрицательную скорость роста w < 0. Профили тающей поверхности приведены на рис. 6. С ростом параметра (времени), вершины начинают заостряться, а впадины, наоборот – уплощаются.

 

Рис. 6. Таяние профиля поверхности z = z(x, t)

Fig. 6. Surface profile melting z = z(x, t)

 

Рассмотренные примеры показывают универсальность математической модели (12), позволяющей описывать широкий круг явлений, связанных с изменением профиля опорного основания транспортных средств.

Особо следует отметить пример 4: одной из важных интерпретаций решений (14) графиков рис. 5 является их соответствие кривой, описываемой центром колеса при качении по неровному профилю дороги [6]. Однако данная интерпретация осложняется нарушением гладкости решений уравнения (12): при малом радиусе колеса  относительно кривизны профиля опорной поверхности, решения (14) соответстуют уравнениям эквидистанты к кривой, заданной параметрически, которые выводятся из элементарных соображений дифференциальной геометрии; при увеличеннии радиуса колеса (или увеличении кривизны профиля) возникает так называемая градиентная катастрофа, которая характеризуется потерей гладкости и однозначности решения. Градиентная катастрофа, в интерпретации качения колеса по неровному профилю, означает критическое значение радиуса колеса , когда возникает двухточечный контакт при перекатывании через яму (лакуну). При значениях радиуса колеса R^*, больших критического значения, имеет место неоднозначность в решении уравнения (12), из-за появления «ласточкиных хвостов» (рис. 7).

 

Рис. 7. Потеря однозначности решений уравнения (12) при качении колеса по неровному профилю

Fig. 7. Loss of uniqueness of solutions of equation (12) when a wheel rolls on an uneven profile

 

С математической точки зрения, градиентная катастрофа наступает в момент потери однозначности эйлерового описания координат. Как следует из (1), на ростке z_p = ξ(x_p, t) имеет место постоянство лагранжевой координаты (т.е. росток является траекторией точки с координатой (x_p, z_p)). Обычно расмотрение изменения профиля проводиться в неподвижной системе координат, в частности, в конкретных фиксированных точках с декартовыми координатами. Такой более естественный способ наблюдения суть эйлерово описание.

Связь эйлеровых и лагранжевых координат дается уравнением, вытекающим из (3):

$z_p\left(t\right)=z\left(x_p\left(t\right),t\right)$.

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

 

Рис. 8. Возникновение неоднозначности эйлерова и лагранжева описаний

Fig. 8. The emergence of ambiguity of Eulerian and Lagrangian descriptions

 

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

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

Например, в случае качения колеса по неровному профилю в области многозначности решения, получить нужное обобщенное решение можно, накладывая дополнительные ограничения: 1) требование однозначности решения; 2) требование максимальности решения.

Рассмотрим, следствием какого интегрального инварианта является дифференциальное уравнение (12), описывающее, в том числе, траекторию центра колеса, катящегося по неровному профилю дороги. Для этого умножим уравнение (12):

$\frac{\partial u}{\partial t}+C\left(u\right)\frac{\partial u}{\partial x}=0$

на величину $\frac{dC\left(u\right)}{dt}$, тогда уравнение допускает вид:

$\frac{\partial v}{\partial t}+v\frac{\partial v}{\partial x}=0$, (15)

где $v\left(x,t\right)=C\left(u\left(x,t\right)\right)$.

Уравнение (15) является уравнением Римана в стандартной форме, которое выражает в дифференциальной форме более общий интегральный закон сохранения. Определить закон сохранения дифференциального уравнения – значит найти инвариант этого уравнения. Для этого представим (15) в дивергентной форме:

$\frac{\partial v}{\partial t}+\frac{\partial }{\partial x}\left(\frac{v^2}{2}\right)=0$.

Интегрируя последнее уравнение по x в бесконечных пределах, считая v(x, t) стремящемся к нулю при x → ±∞, придем к выражению:

$\frac{d}{dt}{\int }_{-\infty }^{+\infty }v\left(x,t\right)dx=0$.

Откуда следует, что инвариантом является выражение:

$I={\int }_{-\infty }^{+\infty }v\left(x,t\right)dx={\int }_{-\infty }^{+\infty }v\left(x\mathrm{,0}\right)dx=const$. (16)

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

 

Рис. 9. Структура разрывов гладкости траекторий центра колеса в области многозначности

Fig. 9. Structure of discontinuities of smoothness of wheel center trajectories in the domain of ambiguity

 

На рис. 9 показаны биссектрисы, соответствующие профилю опорной поверхности. Согласно определению [6], биссектриссой L(M) множества M (область ограниченная сверху профилем опорной поверхности) называется совокупность всех точек из дополнения множества M до плоскости (R² \ M), которая имеет не менее двух проекций на M. Как следует из определения, биссектриссы профиля являются теми кривыми, на которых траектории центра колеса (при различных радиусах колеса R) терпят разрыв гладкости (оставаясь непрерывными). На рис. 9 показана качественная картина поведения биссекторных кривых (перевернутое дерево): биссектрисы с ростом радиуса катящегося колеса попарно сливаются между собой и в пределе при R → ∞ вырождаются в одну единственную биссектрису. Предельный случай R → ∞ можно переформулировать как качение колеса с конечным радиусом, но по микропрофилю, неровностями которого можно пренебречь.

Рассматривают другой тип решений нелинейных дифференциальных уравнений первого порядка, так называемые «слабые» решения [7]. Слабыми решениями называют функции, удовлетворяющие упомянутым уравнениям всюду, где существую производные данных функций, и имеющие разрывы (самих функций или их производных) вдоль кривых (биссекторных кривых рис. 9), которые определяются из дополнительных соображений.

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

作者简介

P. Dmitriev

R.E. Alekseev Nizhny Novgorod State Technical University

编辑信件的主要联系方式.
Email: pavel-9524755555@yandex.ru
ORCID iD: 0009-0007-6344-963X

доцент, канд. техн. наук

俄罗斯联邦, Nizhny Novgorod

A. Anikin

R.E. Alekseev Nizhny Novgorod State Technical University

Email: anikin.zvm@mail.ru
ORCID iD: 0000-0003-0368-4199

начальник НИЛ транспортных машин и транспортных технологических комплексов, д-р техн. наук

俄罗斯联邦, Nizhny Novgorod

Yu. Molev

R.E. Alekseev Nizhny Novgorod State Technical University

Email: moleff@yandex.ru
ORCID iD: 0000-0002-0429-4590

профессор, д-р техн. наук

俄罗斯联邦, Nizhny Novgorod

U. Vakhidov

R.E. Alekseev Nizhny Novgorod State Technical University

Email: umar-vahidov@mail.ru
ORCID iD: 0000-0003-4109-8406

заведующий кафедрой, д-р техн. наук

俄罗斯联邦, Nizhny Novgorod

S. Korostelev

I.I. Polzunov Altai State Technical University

Email: korsan73@mail.ru
ORCID iD: 0009-0003-1382-5589

заведующий кафедрой, доцент, д-р техн. наук

俄罗斯联邦, Barnaul

参考

补充文件