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

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

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

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

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

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

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

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

The operating characteristics of machine parts and units are defined in many respects by the physical and mechanical properties of their surface layers. Unfortunately, it is still not quite clear what parameters and mechanisms are responsible for one or another modification of surface layer properties. In this context, computer modeling techniques can be a useful tool in studying the variation of surface properties in contact interaction and run-in. Of fundamental importance is the possibility to consider the processes occurring on nanoscales and scales of individual atoms. In the work, loading conditions in plastic surface deformation was reproduced on the macroscale (traditional approach), atomic scale, and mesoscale by computer modeling with the finite element method, movable cellular automata method, and molecular dynamics method. The modeling results are in good qualitative agreement with data of experimental measurements.

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

Численно моделируется процесс ударного взаимодействия двух анизотропных тел: цилиндрического ударника и пластины. Моделирование проводится методом конечных элементов в трехмерной постановке. Исследовано влияние ориентации осей симметрии анизотропного материала на разрушение.

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

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

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

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

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

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

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

The paper is a theoretical study into the influence of the content and distribution of soft-phase inclusion agglomerates in the matrix of a ceramic composite specimen on its strength and deformation properties. The movable cellular automata method was used to simulate uniaxial compression of two-dimensional composite material specimens with an aspect ratio of 1:1. It is found that the strength and deformation properties of the generated composites decrease nonlinearly with the growing volume fraction of inclusions. The average size of inclusion agglomerates at the same volume fraction of the soft-phase particles slightly affects the strength and deformation properties of the simulated specimens. The obtained theoretical results can be used to develop new ceramic materials, such as composite ceramics with dimensions preserved at varying temperature.

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

Plastic flow localization in commercially pure titanium (VT1-0 according to the Russian classification) with the surface modified by low-energy high current electron beams has been numerically studied. The structure and mechanical properties of the modified surface layer and titanium substrate correspond to the experimentally observed ones and are taken into account explicitly as initial data of a dynamic boundary value problem. The tension of titanium structures with a modified surface layer is simulated by the finite difference method in a plane strain formulation. The dependence of the plastic strain localization parameters on the mechanical properties of structural elements in the titanium substrate has been determined.Численно исследованы процессы локализации пластического течения в технически чистом титане ВТ1-0, поверхностно модифицированном низкоэнергетическими сильноточными электронными пучками. Структура и механические свойства модифицированного поверхностного слоя и титановой подложки соответствуют экспериментально наблюдаемым и учитываются в расчетах явно в качестве начальных данных краевой динамической задачи. Растяжение структур титана с модифицированным поверхностным слоем моделируется методом конечных разностей в постановке плоской деформации. Установлена зависимость характеристик локализации пластической деформации от механических свойств элементов структуры титановой подложки.

Книга посвящена одному из современных методов исследований в материаловедении - моделированию атомной структуры дефектов. Включает в себя методически последовательное изложение взаимосвязанных разделов в физике реальных кристаллов, таких как структура и энергия точечных и планарных дефектов, структура ядра дислокаций, компьютерное моделирование взаимодействия дефектов, динамические процессы радиационных повреждений в металлах и сплавах; атомистические расчеты структуры основного состояния микрокластеров и кинетики их превращений и др. Перечисленные разделы являются достаточно устоявшимися направлениями в физике реальных кристаллов и ведены в курсы лекций для студентов физико-технического и металлургического факультетов УрФУ. Большое внимание уделяется связи результатов атомистических расчетов с особенностями поведения реальных металлов и сплавов при деформации, разрушении, радиационном повреждении, фазово-структурных превращениях. Систематизирован богатый литературный материал и результаты работы авторов по изучению дислокационной структуры, особенностей деформации и разрушения монокристаллического Ti3Al. Книга адресована научным работникам, магистрам, аспирантам и студентам старших курсов физико-технических специальностей, а также изучающим дефекты и их влияние на свойства твердых тел.

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

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

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

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

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

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

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

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

The regular grid generation method was generalized for a 3D region based on the solid mechanics equations that describe deformation of a finite volume. Test calculations of fracture of ceramic composites with inclusions are represented which are based on the developed model of quasi-brittle media with regard to damage accumulation.

Пространственная разнонаправленность климатических изменений обусловливает необходимость создания прогнозных моделей состояния растительного покрова. Предлагаемое нами решение для пространственного моделирования биопродуктивности природных экосистем создает основу для дальнейшей количественной оценки потенциального поглощения СО2, что в настоящее время считается одной из наиболее актуальных экологических задач. В предлагаемой работе с помощью методов геоанализа условий произрастания отработана оригинальная технология моделирования потенциального распространения биоценозов и их продуктивности, отражающая внутренний аттрактор развития биоценозов под влиянием локальных гидролого-климатических условий их произрастания. Методика реализована для территории юго-востока Западной Сибири в рамках общедоступной ГИС Saga на основе цифровой модели рельефа и данных климатического реанализа WoldClim 2.0. Для прогноза данных на период третьего десятилетия ХХI в. использовался региональный климатический прогноз В.В. Паромонова на основе адаптивной модели-метола экспоненциального сглаживания. Верификация результатов моделирования выполнялась на основе базы данных "Продуктивность экосистем Северной Евразии ". В результате получены пространственнораспределенные данные в виде растров с высоким пространственным разрешением для среднемноголетней потенциальной биопродуктивности по данным реанализа за 1970-2000 гг. и спрогнозированным данным на период 2021-2030 гг. Выявлена как положительная, так и отрицательная тенденция потенциальной биопродуктивности для различных природных зон юго-востока Западной Сибири, обусловленная пространственной разнонапрвленностью изменения теплоэнергетических ресурсов и атмосферных осадков по территории. В целом наблюдаемое потепление на достаточно дренированных территориях способствует увеличению биологической продуктивности агро- и биоценозов, а на переувлажненных - к ее уменьшению.

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

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

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

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

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

Проведена серия натурных экспериментов с применением созданного испытательного стенда «Трибо» в виде перемещаемой бетонной плиты, рассматриваемой как искусственное аллохтонное крыло на шероховатой плоскости сегмента Ангарского разлома в Прибайкалье. Наблюдаемые эффекты контактного взаимодействия неровностей в зоне скольжения фиксировались с применением деформометрической и динамометрической аппаратуры, а также четырех сейсмических станций «Байкал-7HR», широко используемых для регистрации землетрясений. Изучено влияние ударных воздействий на параметры сейсмических колебаний от возникающих источников сейсмических импульсов в процессе меняющегося фрикционного трения. Предложенный способ физического моделирования и полученные данные могут быть полезными для изучения режимов, интерпретации сейсмологических наблюдений, совершенствования среднесрочного прогноза горных ударов и землетрясений, а также представляют интерес для разработки новых физических моделей формирования разномасштабных очагов сейсмической диссипации в тектонических разломах.??.

A series of full-scale experiments was performed to study the influence of impact loads on the parameters of seismic vibrations initiated in variable friction. The study was conducted on a test setup Tribo which is a movable concrete slab modeling an allochthon on a rough plane of the Angara fault segment in Baikal region. Contact interactions of asperities in the slip zone were recorded using strain and load measuring equipment and four seismic stations Baykal-7HR widely used for recording earthquakes. The proposed physical modeling method and obtained results can be of interest for the development of new physical models of differently scaled sources of seismic energy dissipation in tectonic faults and can be useful for seismological studies, related data interpretation, and improvement of extended forecast of rock bursts and earthquakes.

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

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

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

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

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

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

Представлена математическая модель для описания напряженно-деформированного состояния конструкций, изготовленных из ортотропных материалов. На ее основе численно моделируется динамическое нагружение конструкций, по траекториям малой кривизны (по классификации А. А. Ильюшина теории упругопластических процессов). Пластическое деформирвание материала описывется в рамках теории течения. Данная модель описывает упругопластическое деформирование ортотропных материалов, имеющих при разрушении остаточные деформации, не более 6%. К таким материалам относятся: бетоны, углеродистые стал в определенном диапазоне низких температур, серые чугуны, а также многие сплавы. Представлены результаты численного моделирования нагружения преграды из транстропного слава Д16Т стальным ударником в диапазоне скоростей ударного нагружения (400-1200) м/с.

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

В работе рассмотрены физически обоснованные модели пластической деформации и разрушения. Иерархическое моделирование подразумевает явный учет внутренней структуры композита как возможности введения масштабного фактора, введение различных моделей механического поведения пластичных матриц и подложек, хрупких и квазихрупких включений и т.д. для описания разных физических процессов и их взаимовлияния. Проведена серия численных экспериментов по нагружению металлов и сплавов, которые используются в качестве подложки и матрицы при разработке материалов с покрытиями и композитов на металлической основе, в широком диапазоне температур и скоростей деформирования. Предложен критерий разрушения типа Хубера, который учитывает различия в критических величинах для разных типов локальных состояний: растяжение и сжатие. Исследованы процессы разрушения мезообъемов композита Al–Al2O3 и материалов, поверхностно упрочненных методами электронно-лучевой наплавки и диффузионного борирования, в том числе с градиентным подслоем. Показано, что комплексное механическое поведение композиции как целого контролируется взаимосвязанными процессами формирования полос локализованного сдвига в матрице/подложке и растрескивания включений/покрытий.

В работе проведено численное исследование деформации и разрушения мезообъема композиционного материала Al + Al2O3 при растяжении и сжатии. Иерархическое моделирование подразумевает введение различных моделей для описания механического поведения пластичной матрицы и хрупких керамических включений: упруго-пластическая формулировка с учетом деформационного упрочнения и модель растрескивания, использующая критерий разрушения типа Хубера, соответственно. Критерий учитывает различия в критических величинах для разных типов локальных состояний: растяжение и сжатие. Показано, что комплексное механическое поведение мезообъема композита как целого контролируется взаимосвязанными процессами формирования полос локализованного сдвига в матрице и растрескивания включений. Применительно к структурно-неоднородной среде предложенный критерий наибольшей интенсивности касательных напряжений позволяет правильно описать направление распространения трещин при отрыве для разных способов нагружения (растяжение, сжатие).??.

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

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

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

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

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

Рассматривается задача нормального и косого взаимодействия удлиненного стального изотропного цилиндрического ударника с ортотропной пластиной в диапазоне скоростей удара до 2000 м/с в случае, когда вектор скорости ударника не совпадает с его продольной осью. Исследовано влияние угла нутации на разрушение пластины и на поведение ударника в процессе внедрения. Задача решается чилсленно, методом конечных элементов в трехмерной постановке. Поведение материал ударника описывается упругопластической модель., поведение материала преграды описывается в рамках упруго-хрупкой модели.

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

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

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

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

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

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

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

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

Описана методика количественного определения неравновесности микроструктуры материалов по величине коэффициента детектирования Кdet, который рассчитывается программным обеспечением методом дифракции обратнорассеянных электронов. разработанная методика опробована на примере титана технической чистоты ВТ1-0 с крупнозернистой и ультрамелкозернистой структурой, полученной с использованием методов интенсивной пластической деформации по различным технологическим схемам.

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

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

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

Методами дифракции обратнорассеянных электронов проведено исследование изменения микроструктуры приповерхностного слоя никелида титана после импульсных воздействий на поверхность образцов сплава потоками ионов кремния средних энергий. Показано, что после ионно-пучкового облучения поверхности образцов наблюдается изменение и фрагментация структуры приповерхностного слоя на глубину 5-15 мкм, меньшую среднего размера исходного зерна исследуемого сплава. Установлено, что характерными особенностями слоя с фрагментированной структурой являются присутствие в нем мартенситной фазы В19' и высокая концентрация межфазных и внутрифазовых границ раздела, линейные размеры фрагментов превышают 1 мкм, измельчение структуры слоя под облученной поверхностью неоднородно и зависит от кристаллографической ориентации исходного зерна. Высказано предположение, что одной из причин интенсивной фрагментации отдельных зерен исходной фазы В2 после ионно-пучковой обработки является близость ориентации основных систем скольжения направлению воздействия ионными пучками. Возможно, это привело к более раннему, по сравнению с остальными зернами, запуску процесса пластической деформации таких зерен и, как результат, частичной фрагментации их структуры.

We study the problem of wear of a rotationally symmetric profile subjected to oscillations with small amplitude. Under these conditions, sliding occurs at the boundary of the contact area while the inner parts of the contact area may still stick. In a recent paper, Dimaki with colleagues proposed a numerically exact simulation procedure based on the method of dimensionality reduction (MDR). This drastically reduced the simulation time compared with conventional finite element simulations. The proposed simulation procedure requires carrying out the direct and the inverse MDR transformations in each time step. This is the main time consuming operation in the proposed method. However, solutions obtained with this method showed a remarkable simplicity of the development of wear profiles in the MDR space. In the present paper, we utilize these results to formulate an approximate model, in which the wear is simulated directly in the onedimensional space without using integral transformations. This speeds up the simulations of wear by further several orders of magnitude.

We study the problem of wear of a rotationally symmetric profile subjected to oscillations with small amplitude. Under these conditions, sliding occurs at the boundary of the contact area while the inner parts of the contact area may still stick. In a recent paper, Dimaki with colleagues proposed a numerically exact simulation procedure based on the method of dimensionality reduction (MDR). This drastically reduced the simulation time compared with conventional finite element simulations. The proposed simulation procedure requires carrying out the direct and the inverse MDR transformations in each time step. This is the main time consuming operation in the proposed method. However, solutions obtained with this method showed a remarkable simplicity of the development of wear profiles in the MDR space. In the present paper, we utilize these results to formulate an approximate model, in which the wear is simulated directly in the one-dimensional space without using integral transformations. This speeds up the simulations of wear by further several orders of magnitude.

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

The paper reports on theoretical study to elucidate the influence of geometric (width) and mechanical characteristics of phase interfaces on strength, ultimate strain, and fracture energy of metal ceramic composites. The study was performed by computer simulation with the movable cellular automaton method and a well-developed mesoscale structural composite model that takes explicit account of wide transition zones between reinforcing inclusions and the matrix. It is shown that the formation of relatively wide “ceramic inclusions-binder” interfaces with gradual variation in mechanical properties allows a considerable increase in the mechanical properties of the composite. Of great significance is not only the interface width but also the gradient of mechanical properties in the transition zone. The presence of defects and inclusions of nano- and atomic scales in interface regions can increase internal stresses in these regions, induce a steep gradient of mechanical properties in them, and hence decrease strain characteristics and fracture energy of the composite.

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

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