УЧЕТ ТРЕНИЯ В МАТЕМАТИЧЕСКИХ МОДЕЛЯХ ДИССИПАТИВНЫХ СИСТЕМ
Журнал: ВЕСТНИК АСТРАХАНСКОГО ГОСУДАРСТВЕННОГО ТЕХНИЧЕСКОГО УНИВЕРСИТЕТА. СЕРИЯ: УПРАВЛЕНИЕ, ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА И ИНФОРМАТИКА
№ 2
, 2022
Рубрики: МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ
Авторы:
1.
Пензенский государственный технологический университет
(кандидат технических наук, доцент; доцент кафедры автоматизации и управления)
Пенза, Пензенская область, Россия
Пенза, Пензенская область, Россия
2.
Пензенский государственный университет
(доктор технических наук, профессор; профессор кафедры теоретической и прикладной механики и графики)
Россия
Россия
Тип:
Статья
Страницы:
с 110
по 118
Статус:
Опубликован
Получено:
06.01.2022
Одобрено:
08.04.2022
Опубликовано:
22.04.2022
Классификаторы:
УДК 621.01 Теория машиностроения (машиноведение). Механика как теоретическая основа машиностроения
Язык материала:
русский
Ключевые слова:
трение, энергия, диссипативные системы, модели, лагранжиан, зависимость
Аннотация (русский):
Получение моделей механических процессов с диссипацией на основе теории Эйлера–Лагранжа имеет несомненные преимущества перед теорией Ньютона за счет меньшего размера рассматриваемого вектора переменных, входящих в уравнения. Однако вариационная теория Эйлера–Лагранжа не применима к описанию движения систем с диссипацией. Целью работы является возможность продемонстрировать использование теории Эйлера–Лагранжа применительно к диссипативным системам с различными видами трения. Математические модели систем с диссипацией построены на основе суперпозиции механического и термодинамического лагранжианов. Для получения математического описания диссипативных систем предложено использовать полевую теорию применительно к термодинамике диссипативных процессов в рамках формализма Лагранжа. В результате проведенной работы получены уравнения Эйлера–Лагранжа для моделей трения Стокса и Кулона. На основе результатов, полученных в работе, показана возможность учета диссипации энергии в формализме Лагранжа. Предложены математические модели, описывающие динамические процессы в гетерогенных структурах с трением на основе теории Эйлера–Лагранжа. Представлены математические преобразования, которые позволяют осуществлять переход от моделей, полученных на основе формализма Лагранжа, к моделям на основе механики Ньютона
Получение моделей механических процессов с диссипацией на основе теории Эйлера–Лагранжа имеет несомненные преимущества перед теорией Ньютона за счет меньшего размера рассматриваемого вектора переменных, входящих в уравнения. Однако вариационная теория Эйлера–Лагранжа не применима к описанию движения систем с диссипацией. Целью работы является возможность продемонстрировать использование теории Эйлера–Лагранжа применительно к диссипативным системам с различными видами трения. Математические модели систем с диссипацией построены на основе суперпозиции механического и термодинамического лагранжианов. Для получения математического описания диссипативных систем предложено использовать полевую теорию применительно к термодинамике диссипативных процессов в рамках формализма Лагранжа. В результате проведенной работы получены уравнения Эйлера–Лагранжа для моделей трения Стокса и Кулона. На основе результатов, полученных в работе, показана возможность учета диссипации энергии в формализме Лагранжа. Предложены математические модели, описывающие динамические процессы в гетерогенных структурах с трением на основе теории Эйлера–Лагранжа. Представлены математические преобразования, которые позволяют осуществлять переход от моделей, полученных на основе формализма Лагранжа, к моделям на основе механики Ньютона
Ключевые слова:
трение, энергия, диссипативные системы, модели, лагранжиан, зависимость
трение, энергия, диссипативные системы, модели, лагранжиан, зависимость
Введение
Силы трения, относящиеся к силам сопротивления, возникают повсеместно: это внешнее трение в движущихся соединениях деталей машин и механизмов, трение в неподвижных прессовых, заклепочных, болтовых и др. соединениях за счет микропроскальзывания при нагрузках системы
в местах контакта, это и внутреннее трение в материалах элементов системы. К силам трения относят и силы сопротивления среды, которые возникают, например, при движении конструкции в газе или жидкости.
Механическая энергия при трении не исчезает, а происходит ее диссипация – превращение в тепловую энергию. Однако адекватно описать поведение систем под действием трения – нетривиальная задача из-за различных форм проявления трения.
В каждом конкретном случае величина трения зависит от многих факторов (нагрузки, скорости движения, материалов, шероховатости, наличия смазки, температурного воздействия и др.).
Для описания динамики гетероструктур с трением используются два метода. Векторная механика базируется на законах Ньютона и двух основных векторах – «импульсе» и «силе» – и включает анализ и синтез сил и моментов. Второй метод основывается на вариационной теории Эйлера и Лагранжа, базирующейся на скалярных величинах: кинетической и потенциальной энергиях.
В векторном методе для описания движения используются «поддерживающие связи», для которых нужно придумывать различные гипотезы. В скалярном методе все действующие силы (упругости, тяжести и др.) выходят из одной скалярной величины – «силовой функции», выполняющей роль потенциальной энергии. Однако сила трения не связана
с потенциальной энергией и оказывается вне области применения вариационных принципов. Получение адекватного описания моделей механических систем с трением, с опорой на формализм Лагранжа, является актуальной целью исследования из-за явного преимущества перед теорией Ньютона.
Преимущества вариационной теории Эйлера
и Лагранжа очевидны, т. к. энергия является величиной скалярной, ее намного проще преобразовывать и дифференцировать при заменах координат, чем ускорения и силы, входящие в уравнения Ньютона и имеющие векторный характер. Особенно эти преимущества становятся очевидными по мере увеличения размерности системы. При движении одной частицы это преимущество несущественно, однако по мере увеличения числа элементов в рассматриваемой системе, где согласно теории Ньютона для каждого элемента необходимо записывать
3-мерные уравнения движения с учетом взаимных связей, это преимущество возрастает, поэтому использование теории Ньютона сильно усложняет решение задачи с трением, что непременно требует поиска возможных путей применения вариационных методов для упрощения расчетов.
Фундаментальные основы построения моделей гетероструктур с трением
Как уже было сказано, сила трения не имеет «силовой функции», а значит формализм Лагранжа не применим к вариационной теории Эйлера и Лагранжа. Но данное ограничение не является поводом отказываться от применения теории Лагранжа. Формализм Лагранжа основывается на принципе наименьшего действия Гамильтона и является уни-фицированной процедурой, применимой к каждой физической системе [1].
Лагранжева механика является интерпретацией классической механики Ньютона. Как известно, функция Лагранжа для классической механики представляет разность между кинетической и потенциальной энергией и имеет название «функция Лагранжа», или «лагранжиан» [2, 3]:
где L – функция Лагранжа; – кинетическая энергия системы; – потенциал (потенциальная энергия).
Потенциальная и кинетическая энергии содержат в себе информацию о динамике сложных механических систем. Таким образом, формализм Лагранжа представляет собой универсальное методически обобщенное средство, позволяющее описывать различные физические системы на одинаковом методическом уровне. Использование уравнений баланса позволяет определить такие универсальные физические величины, как энергия, работа, импульс. А вся информация о процессах в системе содержится в одной скалярной функции – лагранжиане , который, в свою очередь, зависит от обобщенных координат системы
и обобщенных скоростей
и является индивидуальной структурой, входящей в формализм Лагранжа [4].
Лагранжиан представляет собой интегральное ядро принципа наименьшего действия Гамильтона:
из которого вытекает система уравнений Эйлера–Лагранжа при отсутствии внешних непотенциальных сил
которые и являются фундаментальными уравнениями движения системы.
Формализм Лагранжа для диссипативных систем
Рассмотрим формализм Лагранжа для диссипативных систем.
Идея состоит в обобщении уравнения, описывающего работу силы
где dA – элементарная работа; F – сила, для которой вычисляется элементарная работа; dl – элементарное перемещение.
Смысл обобщения состоит в замене перемещения dl, вдоль которого действует сила F, на обобщенную координату dq, характеризующую пространственную напряженность (перемещение dl, объем dV, массу, заряд, энтропию dS и др.), а вместо F используем обобщенную силу Q, иллюстрирующую характер напряженности, такую, что
Везде в результате действия обобщенной силы происходит преобразование энергии:
– потенциальной
– кинетической
При отсутствии внешних сил получим равенство
(1)
При внешних непотенциальных силах уравнение (1) примет вид
Как показывают исследования [1], диссипация, связанная с внешними непотенциальными силами, может входить в формализм Лагранжа, а силы трения могут быть учтены путем введения переменной переноса.
Диссипативные свойства системы проявляются
в передаче энергии от ее нетепловых степеней свободы к тепловым. Смоделируем трение в рамках формализма Лагранжа, используя термодинамику необратимых процессов. Для этого обобщим принцип Гамильтона на полевую теорию. Представим совокупность фундаментальных полей пространственно-временных переменных системы в виде функции
Процессы в системе в интервале времени в заданном объеме V должны удовлетворять принципу наименьшего действия Гамильтона:
Лагранжиан представляет собой функцию фундаментальных полей и их первых производных по пространственным и временным координатам: , где
Уравнения Эйлера–Лагранжа примут вид
В лагранжевом формализме явление теплопроводности описывается комплексной фундаментальной тепловой переменной
где – абсолютная температура (характеризует тепловое поле); – комплексно-сопряженные функции; – тепловая фаза.
Лагранжиан теплопереноса представим через действительные переменные
где c – удельная теплота (предполагается c = const); ω – частота, в дальнейшем выпадает из уравнений; Θ – константа, имеющая размерность температуры; λ – коэффициент теплопроводности; G(T) – функция энтропии.
В случае использования идеального материала с нулевой теплопроводностью λ = 0 лагранжиан принимает вид
(2)
Наложим полевую теорию на основе модели (2) на формализм Лагранжа.
Рассмотрим одномерное жесткое тело, движущееся по плоской поверхности. Его движение происходит под действием силы трения R и обусловлено потенциалом U(x). Для простоты размеры
и свойства тела примем такими, чтобы на участке движения ∆x тело можно было принять за материальную точку, обладающую теплоизоляционными свойствами (λ = 0), а все тепло при трении передавалось бы данному телу.
Модели трения
Рассмотрим применение формализма Лагранжа в рамках моделей:
– трения Стокса
– трения Кулона
где и – силы трения Стокса и Кулона соответственно; – коэффициент трения; – скорость тела.
Трение Стокса (вязкое трение в жидкости). Уравнение движения имеет вид
(3)
где m – масса; ẍ – ускорение; x – координата тела.
Из него получаем частное уравнение баланса механической энергии
(4)
Скорость потери механической энергии в результате трения определяется мощностью силы трения. Эти потери с противоположным знаком входят в уравнение баланса тепловой энергии
(5)
где – удельная теплота; – абсолютная температура.
Сложив (4) и (5), получим
(6)
Ни одно из уравнений (3)–(6) в отдельности
не является самосопряженным, а следовательно,
не может служить отправной точкой для формализма Лагранжа.
Для возможности применения формализма Лагранжа введем дополнительную переменную переноса Ф. Лагранжиан для задачи с диссипацией принимает вид
где и φ – базовые переменные термодинамики
в рамках лагранжева формализма.
Полученный лагранжиан является суперпозицией механического лагранжиана без учета трения,
усеченного лагранжиана теплодинамики
и лагранжиана взаимодействия
связанного с трением, характеризующего энергию диссипации (первое слагаемое) и скорость диссипации (второе слагаемое).
Уравнения Эйлера–Лагранжа получаются путем вариации переменных x, T, φ, Ф и имеют вид
(7)
(8)
(9)
Решения для уравнений (7)–(9) имеют вид
(10)
где – скорость.
Продифференцировав уравнение (10), получим уравнение для скорости движения материальной точки под действием трения Стокса:
Трение Кулона. Динамика материальной точки в одномерном случае описывается лагранжианом
Уравнения Эйлера–Лагранжа примут вид
В частном случае, при U(x) = 0 и начальных условиях при получим следующие решения [4]:
(11)
Время прекращения движения .
Продифференцировав уравнение (11), получим уравнение для скорости движения материальной точки под действием трения Кулона:
(12)
Результаты моделирования
Проведем анализ полученных уравнений, построив графики зависимостей x(t), v(t), T(t) с помощью программного пакета MathCad.
Для этого зададим параметры системы:
На рис. 1 представлен скриншот, иллюстрирующий применение пакета MathCad при расчетах зависимостей x(t), v(t), T(t) при использовании моделей с трением Стокса (на скриншоте обозначены x(t), V(t), T(t)) и моделей с трением Кулона (на скриншоте обозначены x1(t), V1(t), T1(t)).
Рис. 1. Скриншот применения пакета MathCad при расчетах
Fig. 1. Screenshot of the use of the MathCad package in calculations
Проведем сравнительный анализ решений, полученных классическим методом согласно Ньютоновской механике и с использованием формализма Лагранжа для диссипативных систем.
Анализ моделей с трением Кулона. Из анализа уравнения (11) следует, что оно однозначно описывает движение тела под действием кулоновского трения.
Графики зависимости координаты и скорости от времени, полученные в MathCad согласно уравнениям (11) и (12) представлены на рис. 2.
а б
Рис. 2. Графики зависимости координаты (а) и скорости (б) от времени
при использовании модели трения Кулона в интервале 0–500 с
Fig. 2. Graphs of dependence of a coordinate (a) and speed (б) on time
when using the Coulomb friction model in the time interval 0-500 s
Для классической кулоновской силы трения уравнения зависимости координаты и скорости
от времени для равноускоренного движения имеют вид [5]
(13)
и
(14)
где – ускорение.
Проведя замену на ускорение в уравнениях (11) и (12), получим уравнения (13) и (14), отвечающие равнозамедленному движению и согласующиеся с графиками на рис. 2, что подтверждает справедливость возможности использования формализма Лагранжа для диссипативных систем.
Выбор интервала времени от 0 до 500 обусловлен остановкой тела при v(500) = 0.
Увеличение интервала времени приводит, как
и ожидается, к смене знака скорости и уменьшению координаты (рис. 3), что теоретически оправдано при переходе выделявшейся при трении тепловой энергии в кинетическую.
а б
Рис. 3. Графики зависимости скорости (а) и координаты (б) от времени
при использовании модели трения Кулона в интервале 0–1 000 с
Fig. 3. Graphs of dependence of speed (a) and a coordinate (б) on time
when using the Coulomb friction model in the time interval 0-1.000 s
В реальной ситуации такое происходить не может, т. к. энергия расходуется на работу диссипативных сил, поэтому необходимо с осторожностью подходить к моделированию, определяя интервал времени, при котором уравнения дают адекватный результат.
График зависимости температуры тела от времени (рис. 4) демонстрирует повышение температуры тела в результате трения, что не противоречит физическим законам.
Рис. 4. График зависимости температуры тела с трением Кулона от времени
Fig. 4. Graphs of dependence of body temperature with Coulomb friction on time
Анализ моделей с трением Стокса. Как известно, вязкое трение зависит от скорости движения тела (R = – μx). С уменьшением скорости трение уменьшается. Данный факт подтверждается представленными на рис. 5 графиками зависимости координаты и скорости от времени для вязкого трения.
а б
Рис. 5. Графики зависимости координаты (а) и скорости (б) от времени при использовании модели трения
Fig. 5. Graphs of a coordinate (a) and velocity (б) on time when using a friction model
Анализируя графики на рис. 5, можно прийти
к выводу, что при наличии вязкого трения характер движения оказывается более сложным, чем при трении Кулона. Движение происходит не просто
с замедлением, но и с переменным ускорением, что следует из графика 5, б, а значит, требует для описания в классической форме более сложных зависимостей, нежели зависимости (13) и (14).
Заключение
В работе показана возможность обобщения формализма Лагранжа на задачи с трением для диссипативных систем. Предложены математические модели, описывающие динамические процессы в гетерогенных структурах с трением. Результаты моделирования подтверждают схожесть полученных моделей на основе формализма Лагранжа
с моделями на основе классической ньютоновской механики.
Рассмотренные модели, характеризующие диссипативные свойства тел в формализме Лагранжа, являются упрощенными и нуждаются в совершенствовании. В реальной ситуации движение происходит при соприкосновении тел, имеющих конечные размеры, у которых тепло распределяется определенным образом между обоими телами.
Кроме того, теплопередача происходит по поверхности тел, а не по всему объему, что требует корректировки полученных моделей. Однако задача заключалась не в строгом описании моделей механических систем с диссипацией, а в демонстрации возможности учета диссипации в формализме Лагранжа. Этот формализм позволяет получить унифицированную структуру для любых технических систем и на основе вариационного принципа Гамильтона получать приближенные решения динамических систем с использованием известных пакетов программ.
1. Шарготт М. Трение в деформируемых средах. Подход на основе полевой теории к термодинамике диссипативных процессов в рамках формализма Лагранжа // Физическая мезомеханика. Томск: Ин-т физики прочности и материаловедения Сиб. отд. РАН, 2001. № 4. С. 47-57.
2. Ландау Л. Д., Лифшиц Е. М. Теоретическая физика: учеб. пособие: в 10 т. М.: Физматлит, 2007. Т. 1. Механика. 224 c.
3. Макаров П. А. О вариационных принципах механики консервативных и неконсервативных систем // Вестн. Сыктывкар. ун-та. Сер. 1: Математика. Механика. Информатика. 2017. № 2 (23). С. 46-59.
4. Антони К.-Х. Термодинамика процесса трения и лагранжев формализм: вклад в мезоскопический под-ход в теории пластичности // Физическая мезомеханика. Томск: Ин-т физики прочности и материаловедения Сиб. отд. РАН, 2001. № 4. С. 33-46.
5. Сивухин Д. В. Общий курс физики: учеб. пособие: в 5 т. М.: Физматлит, 2020. Т. 1. Механика. 560 с.
