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

Adhesive contact with exponential adhesive interaction is simulated with the use of the method of dimensionality reduction. The developed procedure is illustrated with an example of adhesion of a cylindrical punch and an elastic half space. However, it is general and can be used for any form of interaction potential and any form of indenter.

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

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.

A multilevel approach is used to numerically investigate physical and mechanical properties of titanium-based bcc alloys and their behavior under conditions identical to selective laser sintering. Elastic properties of beta-Ti-Nb alloy are calculated within the first principles approach. An algorithm is proposed and tested to optimize the calculations and reduce their number by more than 5 times. A molecular dynamics method is employed to study structural changes of titanium and niobium powder particles during sintering and to calculate adhesion characteristics of nanoparticles of the produced alloy depending on the external action. The simulation results are in good agreement with the known experimental data and can be used as input data both for numerical models of a higher spatial scale and for the optimization of production parameters of titanium alloys by additive technologies.В работе с использованием многоуровневого подхода проведено численное исследование физико-механических свойств ОЦК-сплавов на основе титана и изучено их поведение в условиях идентичных селективному лазерному спеканию. В рамках первопринципного подхода вычислены упругие свойства бета-сплава Ti-Nb. Предложен и апробирован алгоритм оптимизации расчетов, позволивший сократить объем вычислений более чем в 5 раз. С использованием метода молекулярной динамики исследовано изменение структуры частиц порошка титана и ниобия в процессе плавления и рассчитаны адгезионные характеристики наночастиц образовавшегося сплава в зависимости от условий внешнего воздействия. Результаты моделирования хорошо коррелируют с известными экспериментальными данными и могут быть использованы в качестве входных данных для численных моделей большего пространственного масштаба, а также для оптимизации параметров процесса получения титановых сплавов с использованием аддитивных технологий.

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

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

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

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

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

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

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

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

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

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

The paper studies tribological aspects of nanostructuring burnishing of steels. The efficiency of the process in improving the tribological properties of steels is assessed as regards the choice of an indenter material and lubricant-coolant reasoning from the friction coefficient at the “indenter-treated part” contact and from the absence of adhesive bond and fatigue microcracks. It is shown that synthetic diamond and dense boron nitride are promising indenter materials for nanostructuring burnishing of corrosion-resistant 20X13 steel and cement 20X steel. It is demonstrated that nanostructuring burnishing increases the wear resistance of structural steels under abrasive action and sliding friction in different media (lubricant, water, air, and argon) due to suppressed processes of microcutting, plastic edging, fatigue and oxidation wear, and adhesive bonding.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Проведено численное исследование деформации и разрушения мезообъема композиционного материала на металлической основе. Иерархическое моделирование подразумевает 1) явный учет внутренней структуры композита, как возможности введения масштабного фактора, 2) введение различных моделей для описания механического поведения пластичной матрицы и хрупких керамических включений, дающих возможность описания различных физических процессов и их взаимовлияния, а именно: пластического течения с учетом деформационного упрочнения и разрушения. Критерий разрушения типа Губера учитывает различия в критических величинах для разных типов локальных состояний - растяжения и сжатия. Показано, что комплексное механическое поведение мезообъема композита как целого контролируется взаимосвязанными процессами формирования полос локализованного сдвига в матрице и растрескивания включений. Установлено, что даже в условиях одноосного внешнего сжатия в композиционном материале могут возникать локальные области растяжения, величина напряжений в которых имеет существенное значение: порядка 30% от величины сжимающих нагрузок для исследуемого типа структуры. Показано, что незначительная модификация критерия наибольшей интенсивности касательных напряжений применительно к структурно-неоднородной среде позволяет правильно описать направление распространения трещин при отрыве для разных способов нагружения (растяжение, сжатие).

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

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

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

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

Численно моделируется воздействие импульса сжатия на элемент ракетного двигателя на твердом топливе. Анализируется напряженно-деформируемое состояние топлива и оболочки.

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

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