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

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

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

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

The paper studies how the stress state of the interface between structural elements in a block-structured medium affects its deformation response to dynamic loading. It is shown that the normalized shear stress and mean stress are the major factors that determine the deformation response of the interface. We propose to describe the dependence of the value of induced irreversible displacement at the interface on the normalized shear stress using a logistic function. The central point of this function is the point of transition from the quasi-elastic to quasi-plastic stage of the interface shear deformation. The obtained empirical dependences are important for understanding the mechanism of irreversible strain accumulation in fault zone fragments and, particularly, for the development of an earlier proposed approach to estimate the characteristic level of active shear stresses in separate tectonic fault regions.

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