ВведениеСилы трения, относящиеся к силам сопротивления, возникают повсеместно: это внешнее трение в движущихся соединениях деталей машин и механизмов, трение в неподвижных прессовых, заклепочных, болтовых и др. соединениях за счет микропроскальзывания при нагрузках системы в местах контакта, это и внутреннее трение в материалах элементов системы. К силам трения относят и силы сопротивления среды, которые возникают, например, при движении конструкции в газе или жидкости. Механическая энергия при трении не исчезает, а происходит ее диссипация – превращение в тепловую энергию. Однако адекватно описать поведение систем под действием трения – нетривиальная задача из-за различных форм проявления трения.В каждом конкретном случае величина трения зависит от многих факторов (нагрузки, скорости движения, материалов, шероховатости, наличия смазки, температурного воздействия и др.).Для описания динамики гетероструктур с трением используются два метода. Векторная механика базируется на законах Ньютона и двух основных векторах – «импульсе» и «силе» – и включает анализ и синтез сил и моментов. Второй метод основывается на вариационной теории Эйлера и Лагранжа, базирующейся на скалярных величинах: кинетической и потенциальной энергиях.В векторном методе для описания движения используются «поддерживающие связи», для которых нужно придумывать различные гипотезы. В скалярном методе все действующие силы (упругости, тяжести и др.) выходят из одной скалярной величины – «силовой функции», выполняющей роль потенциальной энергии. Однако сила трения не связана с потенциальной энергией и оказывается вне области применения вариационных принципов. Получение адекватного описания моделей механических систем с трением, с опорой на формализм Лагранжа, является актуальной целью исследования из-за явного преимущества перед теорией Ньютона.Преимущества вариационной теории Эйлера и Лагранжа очевидны, т. к. энергия является величиной скалярной, ее намного проще преобразовывать и дифференцировать при заменах координат, чем ускорения и силы, входящие в уравнения Ньютона и имеющие векторный характер. Особенно эти преимущества становятся очевидными по мере увеличения размерности системы. При движении одной частицы это преимущество несущественно, однако по мере увеличения числа элементов в рассматриваемой системе, где согласно теории Ньютона для каждого элемента необходимо записывать 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 timewhen 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 timewhen 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). ЗаключениеВ работе показана возможность обобщения формализма Лагранжа на задачи с трением для диссипативных систем. Предложены математические модели, описывающие динамические процессы в гетерогенных структурах с трением. Результаты моделирования подтверждают схожесть полученных моделей на основе формализма Лагранжа с моделями на основе классической ньютоновской механики. Рассмотренные модели, характеризующие диссипативные свойства тел в формализме Лагранжа, являются упрощенными и нуждаются в совершенствовании. В реальной ситуации движение происходит при соприкосновении тел, имеющих конечные размеры, у которых тепло распределяется определенным образом между обоими телами. Кроме того, теплопередача происходит по поверхности тел, а не по всему объему, что требует корректировки полученных моделей. Однако задача заключалась не в строгом описании моделей механических систем с диссипацией, а в демонстрации возможности учета диссипации в формализме Лагранжа. Этот формализм позволяет получить унифицированную структуру для любых технических систем и на основе вариационного принципа Гамильтона получать приближенные решения динамических систем с использованием известных пакетов программ.