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

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

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

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

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

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

В работе исследуются процессы деформации и разрушения стали с боридным покрытием. Краевая динамическая задача решается численно методом конечных разностей в постановке плоской деформации. Геометрия границы раздела «покрытие - подложка» соответствует экспериментально наблюдаемой и задается в расчетах явно. Для описания механической реакции стальной основы используется упругопластическая модель изотропно упрочняющегося материала, которая включает соотношения для описания медленных течений. Исследованы особенности локализации пластической деформации и разрушения в условиях распространения полосы Чернова-Людерса в стальной подложке при растяжении.