Численное моделирование динамического контактного взаимодействия упругопластических тел тема автореферата и диссертации по механике, 01.02.04 ВАК РФ

Введение

Глава 1. АЛГОРИТМ ЧИСЛЕННОЙ РЕАЛИЗАЦИИ ГРАНИЧНЫХ УСЛОВИЙ КОНТАКТНОГО ВЗАИМОДЕЙСТВИЯ.

1.1 Формулировка граничных условий контакта с учетом сил трения

1.2 Эквивалентность квазивариационного неравенства закону трения Амонтона - Кулона.

1.3 Вычислительный алгоритм

Глава 2. ПРИМЕНЕНИЕ АЛГОРИТМА В ДИНАМИЧЕСКИХ

ЗАДАЧАХ

2.1 Модель упругопластического деформирования при конечных поворотах

2.2 Одномерная схема с реконструкцией.

2.3 Алгоритм двуциклического расщепления

2.4 Процедура корректировки напряжений

2.5 Тестовые расчеты.

2.5.1 Задача о контакте упругого слоя с жесткой плоскостью

2.5.2 Задача о контакте упругопластического слоя с жесткой плоскостью

2.5.3 Задача о вдавливании жесткого цилиндрического штампа в упругую пластину

Глава 3. РЕЗУЛЬТАТЫ ЧИСЛЕННЫХ РАСЧЕТОВ.

3.1 Механическая обработка образца

3.2 Моделирование образования стружки в процессе обработки заготовки резанием.

3.3 Волнообразование при сварке взрывом

3.3.1 Обзор основных результатов о механизме волнообразования

3.3.2 Задача о потере устойчивости слоя

3.3.3 Задача о бегущей нагрузке

3.3.4 Косое соударение пластин

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

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

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

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

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

Теория Г. Герца основана на следующих допущениях:

- поверхности взаимодействующих тел гладкие и несогласованные;

- деформации малы (а значит, применима линейная теория упругости);

- размеры области контакта малы по сравнению с размерами каждого из контактирующих тел и с радиусами кривизны их поверхностей;

- трение на поверхности контакта отсутствует.

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

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

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

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

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

Идея подхода впервые сформулирована А. Синьорини в 1933 году в работе [3] при изучении задачи о равновесии упругого тела, соприкасающегося без трения с плоским абсолютно жестким основанием (штампом), под действием заданной системы сил. Вторично и с большей полнотой Синьорини изложил свои результаты в 1959 году [4]. Исследованная им задача сводится к минимизации функционала упругой энергии J на множестве полей перемещений, удовлетворяющих в области контакта ограничению un\dQ < 0, где - плоский участок границы упругого тела Q, п - единичная внешняя нормаль к <90. Ограничение описывает условие непроникания точек тела внутрь штампа. Задача Синьорини эквивалентна вариационному неравенству J'u (и* — и) > 0, где и ж и* принадлежат множеству К = {и* : и*п\дп < 0} допустимых вариаций решения, J'u - вариация или слабая производная функционала J.

Долгое время, вплоть до широкого использования вычислительной техники в конце 50-х - начале 60-х годов, эффективных методов численного решения задач с односторонними ограничениями не существовало, и потому интерес к задаче Синьорини и ее обобщениям проявлялся в основном теоретиками. Были доказаны различные теоремы о существовании решений - Г. Фикерой, Г. Дюво и Ж.-Л. Лионсом, X. Леви и др.; обзор этих работ можно найти в книге [5]. Наиболее весомый вклад в разработку проблемы с позиций функционального анализа и теории уравнений с частными производными был внесен, по-видимому, итальянским математиком Г. Стампаккья (обзор результатов Стампаккьи насчитывает более 80 наименований; см., например, [6]). Важные результаты исследований в области теории вариационных неравенств изложены и в работах других математиков [7 - 12].

Методы теории вариационных неравенств применялись к исследованию разрешимости контактных задач механики деформируемого твердого тела с неизвестной границей в работах А. С. Кравчука [13 - 15] (без учета сил трения) и [16] (с учетом трения). Так, например, в [16] приведены постановки квазистатической и динамической контактных задач в приращениях для общего случая граничных условий на поверхности контакта, а также доказательство существования и единственности обобщенного решения. Полученные Кравчуком результаты наиболее полно изложены в монографии [17].

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

В монографии [18] собраны необходимые для обоснования численных методов теоретические результаты. Численным методам решения упру-гопластических задач на основе теории вариационных неравенств посвящены работы [19, 20], где используется метод локальных вариаций при исследовании контактной задачи о встречном прогибе двух упругих пластин, а также при решении задач упругопластического кручения призматических стержней. Влияние сил трения на поведение тел в зоне контакта в этих работах не учитывалось. Упругие и упругопластические контактные задачи с учетом трения рассматривались в [21 - 35].

Исторически первым и наиболее широко известным и используемым в расчетах является закон сухого трения Амонтона - Кулона (чаще называемый просто законом Кулона), связывающий касательные компоненты вектора усилий контактного взаимодействия с вектором скорости движения одного из контактирующих тел относительно другого. Рассмотрим, к примеру, задачу о контакте деформируемого тела с абсолютно жестким неподвижным штампом. В этом случае закон трения Амонтона - Кулона в области контакта (возможно, неизвестной заранее и определяемой в процессе решения) представляет собой совокупность следующих условий: knrl < fKnl если vT = О, ' \(Гпт\ = fWnnl И т^г = -г-^г, если VT Ф 0.

I^nrl Frl

Здесь апп и сгпт - нормальная и касательная составляющие вектора напряжений ап, действующего на площадке деформируемой поверхности' с нормалью п, / - коэффициент трения скольжения, vr - касательная составляющая вектора скорости v = du/dt, t - время или любой другой параметр, определяющий последовательность смены состояний деформируемого тела. Закон трения Кулона описывает как случай сцепления, когда vT = 0, так и случай проскальзывания, когда выполняется равенство crnTj\crnT\ = — vT/\vT\. Этот закон был подвергнут обширной экспериментальной проверке; в настоящее время хорошо изучены пределы его применимости и известны многочисленные обобщения. Несмотря на ограниченность и схематичность закона Кулона, он во многих случаях дает удовлетворительные результаты по прогнозу усилий контактного взаимодействия соприкасающихся деталей машин и элементов сооружений, а также представляет собой хороший "полигон" для отработки математического аппарата исследования и решения конкретных задач с учетом трения. Одно из обобщений закона трения Амонтона - Кулона, использованное, в частности, в работах А. А. Ильюшина [36] и С. С. Григоряна [37], представляется в следующем виде:

Wnr\ < f\&nn| и одновременно \<гпт\ < rs, если vT = 0, &пт\ = fKn\ ИЛИ \(ТПТ\ = rs, и -—т~ = --— если vT ф 0,

Тпт I 1*4 где ts - либо константа, либо функция деформаций и скоростей деформаций, имеющая смысл предела текучести материала на сдвиг.

По-видимому, первой работой по применению вариационного подхода к исследованию контактных задач теории упругости с учетом трения в строгой математической постановке была работа Г. Дюво и Ж.-Л. Лион-са, опубликованная в 1971 году; полное изложение полученных результатов с некоторыми обобщениями дано в [8]. Одним из основных предположений работ Дюво и Лионса по данной проблеме является предположение о том, что нормальное давление в зоне контакта - известная функция координат точек поверхности тела; данное предположение было снято в работе А. С. Кравчука [16]. Следует отметить, что во многих практически интересных задачах взаимное влияние нормальных и касательных напряжений невелико: можно вначале определить нормальную компоненту апп вектора напряжений на контактной поверхности, считая трение нулевым, а после этого решать задачу об определении касательной компоненты стпт, считая апп заданным. Более того, существуют случаи, когда задачи определения величин апп и стпт разделяются точно; например, так будет в контактных задачах в постановке Герца, когда контактирующие тела заменяются полупространствами, причем коэффициенты Пуассона совпадают. Еще одним весьма существенным предположением теории, развитой Г. Дюво и Ж.-Л. Лионсом, является предположение о возможности использования в квазистатических задачах закона трения Амонто-на - Кулона с перемещениями вместо скоростей. Условия, при которых такое предположение допустимо, установлены Кравчуком [16].

Оригинальный подход к решению контактных задач с учетом трения развит в ряде работ А. В. Вовкушевского и Б. А. Шойхета (см., напр., [21]): величина ипп предполагается заданной; задача сводится к минимизации функционала; варьируемые функции удовлетворяют условию непроникания; применяются теоремы о выпуклых функционалах, заданных на выпуклых множествах. Еще один подход заключается в интерпретации ограничения в контактной задаче с трением как ограничения в задаче о контакте деформируемых тел с идеальными односторонними связями, но для случая, когда контактирующие поверхности имеют микрорельеф ("зубцы"). Такая интерпретация позволяет применять для решения задачи с трением методы, которые ранее были развиты для решения контактных задач без трения, но с микрорельефом. При решении контактных задач для тел с шероховатыми поверхностями обычно деформации микронеровностей, образующих шероховатость, моделируются включением между контактирующими поверхностями некоторого упругого слоя. В статье Вовкушевского [22] аналогичный подход используется для вариационной постановки в духе задачи Синьорини с трением, причем учитывается податливость слоя по направлениям нормали и касательной. Показано, что включение упругого слоя приводит к качественным изменениям исходной задачи и может рассматриваться как способ ее регуляризации.

В цикле работ Р. В. Гольдштейна и А. А. Спектора [23 - 27] исследуются контактные задачи с сухим трением при неизвестной границе между областями проскальзывания и сцепления. Задача определения касательных напряжений и смещений сводится к эквивалентной вариационной задаче на площадке контакта, что дает возможность получить ее эффективное численное решение в трехмерном случае. Устанавливается существование решения, а в статическом случае и его единственность. Рассматриваются статическая задача о контакте двух упругих тел из одинаковых материалов и квазистатическая задача о стационарном качении упругого тела по основанию из того же материала. Исследуется также поведение во времени решения нестационарной задачи о взаимодействии движущегося упругого тела и упругого основания; для этой задачи доказывается, что поля касательных напряжений трения и скоростей проскальзывания со временем стремятся к стационарным. Строятся вариационные оценки решений задач с неизвестной площадкой контакта и задач о трещинах с частично налегающими поверхностями. Предполагается, что развитие трещины сопровождается образованием областей, где ее поверхности приходят в контакт. В неизвестных заранее зонах налегания имеется трение с коэффициентом, зависящим от нормального давления и величины относительного касательного смещения поверхностей. Рассматривается вариационная задача минимизации негладких неквадратичных функционалов, зависящих от скачков смещений в области разреза. Устанавливаются некоторые свойства решения задачи, а также качественные результаты о зависимости его интегральных характеристик от формы трещины и параметров закона трения.

В статье Ю. И. Няшина и С. А. Чернопазова [28] развивается вариационный подход к решению упругой контактной задачи с учетом трения, подчиняющегося закону Кулона. Рассматривается квазистатический процесс деформирования упругого тела. В связи с трудностями при исследовании разрешимости задач с трением, область контакта считается известной и не изменяющейся при нагружении. Исследуется влияние касательных перемещений на контактное давление. Формулируются достаточные условия, при выполнении которых это влияние отсутствует. С помощью оператора влияния перемещений на контактное давление доказывается существование и единственность обобщенного решения вариационной постановки задачи. 3. Досталь и В. Вондрак в [29] представили метод численного решения квазивариационных неравенств, описывающих равновесие упругих тел в контакте с трением. Задача сводится к последовательности хорошо обусловленных задач с заданным трением, которые переформулируются в задачи квадратичного программирования с ограничениями. Затем алгоритм решения этих задач применяется к решению результирующей контактной задачи с кулоновым сухим трением.

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

Большой класс задач составляют задачи для тел с учетом пластических деформаций; основополагающими здесь являются работы Р. Хил-ла [38], П. П. Мосолова и В. П. Мясникова [39]. Идеи и методы А. С. Кравчука в применении к упругопластическим контактным задачам используются в цикле работ В. И. Кузьменко [31 - 34], а также в статье Ю. И. Ня-шина и С. А. Чернопазова [35]. Так, например, в работах Кузьменко рассматриваются задачи определения напряженно-деформированного состояния упругопластического тела при контактном взаимодействии с жестким штампом. Исследуется квазистатический процесс деформирования при малых перемещениях и деформациях. Накладывается требование взаимного непроникания тела и штампа, а также требование отсутствия растягивающих напряжений на контактной поверхности. Для описания касательного взаимодействия тела и штампа используется закон трения Амонтона - Кулона; допускается, что шероховатость контактирующих поверхностей может быть различной в различных точках; поэтому коэффициент трения, характеризующий шероховатость поверхностей, считается функцией положения точки. С точки зрения касательного взаимодействия поверхность возможного контакта деформируемого тела и штампа можно в каждый момент времени разбить на три части: на одной из частей контакт отсутствует, на второй части происходит сцепление контактирующих поверхностей, а на третьей осуществляется контакт и взаимное скольжение поверхностей тела и штампа. Линии раздела трех указанных областей заранее неизвестны и должны быть определены при решении задачи. На остальной части поверхности деформируемого тела заданы переменные нагрузки и перемещения. Показано, что исходная контактная задача в обобщенной постановке эквивалентна некоторому квазивариационному неравенству относительно скоростей перемещений. Исследуются вопросы существования и единственности решения неравенства. Излагается методика численного решения: дискретизация квазивариационного неравенства по времени, в результате чего получается последовательность вариационных неравенств, решение которых сводится к поиску минимумов квадратичных функционалов на множествах допустимых функций; дискретизация экстремальных задач при помощи метода конечных элементов; решение возникающих конечномерных задач нелинейного программирования с использованием обобщенного метода последовательной верхней релаксации. Приводится пример решения задачи контактного взаимодействия кусочно-однородной упругопластической полосы и выпуклого штампа. Численно получены эпюры контактного давления и конфигурации пластической зоны в различные моменты времени при различных траекториях движения штампа. В процессе счета обнаружено, что, несмотря на значительное различие зон пластических деформаций в процессе деформирования, в окончательном положении штампа размеры и форма пластических областей, соответствующих различным траекториям движения штампа, достаточно близки. В работе Няшина и Чернопазова [35] рассматривается дифференциальная и вариационная постановка квазистатической задачи контактного взаимодействия упругопластического тела с жесткой опорой. В качестве определяющих соотношений приняты уравнения теории пластического течения с изотропным упрочнением, являющиеся частным вариантом теории упругопластичес-ких процессов А. А. Ильюшина. Приводится доказательство существования и единственности обобщенного решения. В [40] рассмотрены математические аспекты квазистатической контактной задачи об ударе торца прямолинейного стержня о поперечную преграду. Построена модель задачи, основанная на законе течения Прандтля - Рейсса. Приведены две вариационные постановки задачи, первая - в терминах напряжений и скоростей, вторая - в напряжениях. Применен подход, в котором использована вязкая регуляризация решения. Существование и единственность решения установлены с помощью теории эволюции решений и методов выпуклого анализа. Приведены некоторые априорные оценки приближенных решений. Кроме того, регуляризация вязкоупругопластической задачи может рассматриваться как основа численного анализа, для этого предложен алгоритм, сходимость которого подтверждена.

Под руководством Н. X. Арутюняна выполнен цикл работ по решению контактных задач для растущих деформируемых тел, обладающих сложными реологическими свойствами; результаты подытожены в монографиях [41, 42], где излагаются основы механики растущих тел применительно к задачам взаимодействия последовательно возводимых сооружений с деформируемыми основаниями, действия штампов на массивные тела изменяющейся вследствие наращивания конфигурации, контакта деформируемых тел в процессах типа напыления, осаждения, фазовых превращений, намораживания и т. п. А. А. Спектором [43], И. Г. Горячевой и М. Н. Добычиным [44], В. М. Александровым и Е. В. Коваленко [45, 46] и рядом других авторов развивается направление по учету износа контактирующих поверхностей.

Вариационный подход к решению контактных задач с неизвестными границами для оболочек и пластин предложен в работах Г. И. Львова [47] и А. М. Хлуднева [48 - 58]. В [47] с применением вариационных неравенств задача о взаимодействии тонкой пологой оболочки с абсолютно жестким штампом без учета трения в области контакта сводится к задаче минимизации функционала Лагранжа на множестве допустимых перемещений; доказываются существование и единственность решения при определенных допущениях о свойствах диаграммы деформирования. В статьях А. М. Хлуднева [48 - 54] доказана разрешимость широкого класса упругопластических краевых задач, формулируемых с использованием вариационных неравенств; подробно исследована проблема оптимального управления в контактных задачах типа Синьорини для пластин и оболочек. Так, например, в [50] рассмотрена квазистатическая задача упругопластического деформирования оболочек, а в [52] - задача о контакте упругопластической пластины с жестким штампом в точной постановке; для доказательства существования решения вариационного неравенства в [52] предложена процедура введения штрафных операторов; исследована связь между задачей о контакте упругопластической пластины со штампом и контактной задачей для жесткопластической пластины. В [55 - 58] решаются задачи деформирования упругих и неупругих тел, содержащих трещины. На берегах трещины задается краевое условие, имеющее вид неравенства и моделирующее условие взаимного непроникания берегов. Найден полный набор естественных краевых условий в виде системы равенств и неравенств, выполняющихся на берегах трещины. Предложенный подход отличается от классического, в котором краевые условия на берегах трещины имеют вид равенств.

Численному решению вариационных неравенств посвящены работы В. А. Ковтуненко [59 - 61], в которых с помощью метода штрафа дифференциальные неравенства аппроксимируются нелинейными дифференциальными уравнениями, линеаризуемыми итерационными методами. В [59] рассмотрено классическое вариационное неравенство, описывающее задачу о контакте упругой пластины с жестким препятствием, предложены итерационные методы аппроксимации данного неравенства с использованием оператора штрафа, доказаны результаты о сходимости решений и найдены оценки погрешности. Проведен ряд численных экспериментов по изучению поведения упругой квадратной пластины, защемленной по краям и находящейся под действием жесткого препятствия. Неизвестная область контакта пластины с препятствием находится после определения нормального прогиба пластины; рассматривается зависимость прогиба, контактных усилий и изгибающих моментов от нагрузки. В [60] рассмотрена система вариационных неравенств, описывающая контакт упругопластической пластины Кирхгоффа, жестко защемленной по краям, с жестким штампом при неизвестных границах области контакта и зоны пластичности. В [61] исследуется задача равновесия упругопластической балки, шарнирно закрепленной по краям, под действием жесткого штампа. Для описания поведения материала используется модель пластического течения. Итерационно построена система линейных дифференциальных уравнений и доказана сходимость ее решения к решению исходной задачи. Численные расчеты по построенной схеме проиллюстрированы на модельных примерах.

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

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

Стационарная динамическая контактная задача для абсолютно жесткого штампа, движущегося с постоянной дозвуковой скоростью вдоль границы упругой полуплоскости, в рамках классической линейной теории упругости исследована впервые JI. А. Галиным [63]; им получено точное решение этой задачи. Это одно из немногих точных решений в динамических контактных задачах с заранее неизвестной зоной контакта. В монографии [64] рассмотрено качение вязкоупругого цилиндра по вязкоупругой полуплоскости в наиболее общем случае, когда площадка контакта имеет участки с трением и сцеплением. Л. А. Галин [64] исследовал также задачу о напряженном состоянии, возникающем в упругом теле при вдавливании в него одного или нескольких движущихся с постоянной скоростью штампов произвольной формы с учетом трения Кулона. При решении контактных задач им использованы методы теории функций комплексного переменного и теории потенциалов.

Исследованию нестационарных динамических задач с изменяющейся в процессе деформирования областью контакта посвящено большое число работ. Так, например, в статье А. И. Гулидова и В. М. Фомина [65] рассматриваются контактные задачи о продольном ударе цилиндрического и усеченного конического стержней об абсолютно твердую преграду в двумерной постановке. Численно моделируется процесс отскока упруго-пластического стержня конечной длины от преграды с помощью модифицированного метода М. Л. Уилкинса. Проведенные серии расчетов при различных скоростях соударения позволили выяснить картину распространения упругопластических волн и проанализировать явление отскока. Сравнение результатов расчетов для стержней различной формы показало, что даже при малых углах конусности время отскока от преграды усеченного конического стержня существенно отличается от времени отскока цилиндрического стержня. Это явление связано с процессом кумуляции упругой волны нагрузки и разгрузки из-за изменения площади сечения и их взаимодействия с падающей пластической волной. В работе В. Г. Ковалева и В. М. Косенкова [66] изложены постановка, метод и результаты решения задачи о контактном взаимодействии вложенных друг в друга толстостенных упругопластических цилиндров конечной длины под действием импульсной нагрузки на внутреннюю поверхность внутреннего цилиндра. На основе разработанного алгоритма исследовано влияние закона нагружения, толщины внутреннего цилиндра, начальной величины и формы зазора между цилиндрами, наличия в зазоре газа на величину остаточных деформаций цилиндров. В статье Е. В. Коваленко и В. Б. Зеленцова [67] рассмотрены контактная задача об антиплоском сдвиге штампом упругого полупространства и плоская задача о вдавливании штампа в упругую полуплоскость. Для решения этих неустановившихся динамических задач применяются преобразования Лапласа -Карсона по времени и Фурье по пространственной координате; в результате задачи сводятся к интегральным уравнениям Фредгольма первого рода. Получены простые аналитические решения для больших и малых времен взаимодействия штампа с полупространством.

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

Воздействие взрывных и ударных нагрузок на металлические образцы и конструкции различной формы в большинстве случаев реализуется в условиях несимметричного нагружения. Теоретические исследования пространсвенных задач несимметричного динамического взаимодействия контактирующих тел связаны с большими трудностями. В то же время характерные особенности деформирования и разрушения твердых тел при несимметричном взаимодействии можно качественно проследить теоретически в плоской постановке. Так, в работе И. Е. Хорева, В. А. Горе льского, С. А. Зелепугина и В. Ф. Толкачева [69] проведено численное моделирование кратерообразования и разрушения плоских преград методом конечных элементов в двумерной постановке. Соударяющиеся тела описываются сжимаемой упруговязкопластической 'средой, поведение которой при динамических нагрузках характеризуется модулем сдвига, динамическим пределом текучести и вязкостью. Проведено исследование процесса взаимодействия стального плоского ударника со стальной пластиной при различных скоростях соударения. Изучены конфигурации взаимодействующих тел, поля скоростей в них, а также распределения изолиний удельного объема трещин, описывающих характер разрушения преграды в процессе проникания ударника. Сравнение с экспериментальными данными показало, что расчеты, проведенные для плоского случая, качественно верно отражают существенно трехмерный процесс взаимодействия ударников с преградами при ударе под углом.

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

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

В статьях А. И. Гулидова и др. [74 - 76] разработан численный алгоритм для решения задач высокоскоростного взаимодействия упругоплас-тических тел в двумерном случае. В качестве одного из примеров применения алгоритма рассмотрена задача о пробитии алюминиевых преград свинцовой пулей. Численно обнаружено, что процесс пробивания существенно зависит от типа преграды (монолитная преграда, преграда с контактирующими слоями, преграда с разнесенными слоями).

В статье А. Б. Киселева [77] показано, что для определения предела прочности материала можно использовать эксперименты по откольному разрушению при плоском соударении металлических пластин. Предел прочности определяется в результате сравнения данных единичного расчета процесса соударения на ЭВМ с данными эксперимента.

В работах В. В. Алехина, С. Н. Коробейникова и др. [78 - 80] представлен численный алгоритм для решения геометрически и физически нелинейных статических и динамических контактных задач методом конечных элементов. Сначала определяются геометрически взаимные проникновения контактирующих тел; затем из решения уравнений равновесия (движения) определяются контактные силы, препятствующие этим проникновениям, или методом множителей Лагранжа, или методом штрафных функций. Приводятся результаты тестовых расчетов для трехмерных контактных задач. В. В. Егуновым и А. Б. Конюховым [81] разработана методика численного моделирования квазистатического и динамического контактного взаимодействия упругопластических тел на основе метода конечных элементов, неявных схем интегрирования по времени (схемы Хаболта и Ньюмарка) и метода штрафных функций.

В динамических задачах граничные условия контакта с учетом сил трения могут быть сформулированы в виде квазивариационных неравенств. В работах А. С. Кравчука [16, 17], как уже упоминалось ранее, приведена постановка динамической контактной задачи в приращениях для общего случая граничных условий на поверхности контакта, а также доказательство существования и единственности обобщенного решения. В [82] рассмотрены модели контактного взаимодействия твердых деформируемых тел как без трения, так и с учетом трения по различным теориям. Особое внимание уделено новой модели, учитывающей силы адгезионного сцепления. Эта модель, по существу, представляет собой теорию разрушения, в том числе с учетом вязкости. Алгоритмы базируются на идеях теории оптимизации.

Один из общих подходов к исследованию динамических задач в рамках теории упругопластического течения Прандтля - Рейсса основан на формулировке определяющих соотношений необратимого деформирования в виде принципа максимума скорости диссипации энергии (см. [8, 11]). Этот подход развивается в работах Б. Д. Аннина и В. М. Садовского [83 - 85]. Для описания процесса упругопластического деформирования вводится вариационное неравенство с линейным оператором гиперболического типа. Граничные условия контактного взаимодействия тел формулируются в виде квазивариационных неравенств с односторонними ограничениями. Предложен алгоритм численного решения динамических задач, первый этап которого состоит в расчете на один шаг по времени вспомогательного решения системы дискретных уравнений, аппроксимирующих систему уравнений динамической теории упругости, а второй этап - в корректировке решения, включающей в себя процедуру вычисления проекции определенной комбинации основного и вспомогательного решений на множество допустимых вариаций. В упругопластических задачах этот алгоритм обеспечивает автоматическое выполнение в зоне необратимой деформации условия пластичности и условия положительности скорости диссипации энергии. Предложенные способы корректировки решения можно рассматривать как обобщение процедуры корректировки М. JI. Уилкинса, обладающее по сравнению с этой процедурой меньшей дополнительной (схемной) диссипацией энергии и большей точностью. В качестве примера применения алгоритмов корректировки решения к исследованию геометрически нелинейных задач в [83, 85] рассмотрена контактная задача динамического деформирования упругопластической плиты (как однородной, так и слоистой) на оправке произвольной формы. Разработана модель деформирования слоя при малых деформациях и произвольных по величине поворотах элементов. Как внутренние граничные условия на поверхности контакта соседних слоев плиты после расслоения, так и граничные условия контакта нижнего слоя с оправкой сформулированы в виде квазивариационных неравенств с ограничениями на допустимые перемещения. Трение на поверхности контакта не учитывалось. Представлены результаты расчетов динамики однородных и слоистых плит под действием импульсной нагрузки взрывного типа.

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

Данная работа посвящена численному моделированию динамического контактного взаимодействия упругопластических тел.

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

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

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

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

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

В первом параграфе описана модель упругопластического деформирования, предложенная в работе Б. Д. Аннина и В. М. Садовского [83]. Она основана на разложении тензора градиентов деформации, определяющего линейное преобразование элемента тела из отсчетного состояния в актуальное, в произведение ортогонального тензора конечного поворота и симметричного тензора, компоненты которого в предположении малости деформаций отождествляются с компонентами единичного тензора. Эта модель состоит из уравнений движения, закона Гука для упругих составляющих деформации, принципа максимума скорости диссипации энергии, описывающего процесс пластического формоизменения, и уравнения для угла поворота. Переход материала из упругого состояния в пластическое определяется условием пластичности Мизеса. Объемная деформация происходит по линейному упругому закону. Зависимость предела текучести и параметров упругости от величины механических напряжений и температуры не учитывается.

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

Во втором параграфе описана монотонная ENO-схема решения одномерных задач, являющаяся уточнением схемы распада разрыва С. К. Годунова. В этой схеме на шаге "предиктор" инварианты Римана определяются методом характеристик как кусочно-линейные сплайны. Для определения производных от инвариантов используется итерационный процесс с предельной реконструкцией, предложенный в статье В. Ф. Ка-менецкого и А. Ю. Семенова [87]. Реконструкция сводится к нахождению монотонных кусочно-линейных сплайнов, имеющих минимальные разрывы на границах соседних ячеек.

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

В четвертом параграфе описана процедура корректировки напряжений (см. [85]), которая получается путем аппроксимации определяющих соотношений необратимого пластического течения материала. Известный вариант такой процедуры предложен М. JI. Уилкинсом [89].

В. М. Садовским [85] разработана корректировка, позволяющая сохранить второй порядок точности схемы.

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

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

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

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

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

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

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

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

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

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

Результаты работы докладывались на XXXIV и XXXV Международт ных научных студенческих конференциях "Студент и научно-технический прогресс" (Новосибирск, 1996 г. и 1997 г.); конференциях молодых ученых ИВМ СО РАН (Красноярск, 1997 г. и 1999 г.); Международных конференциях "Математические модели и методы их исследования" (Красноярск, 1997 г. и 2001 г.); Втором Всероссийском семинаре "Теория сеточных методов для нелинейных краевых задач" (Казань, 1998 г.); Второй и Третьей Сибирских школах-семинарах "Математические проблемы механики сплошных сред" (Новосибирск, 1998 г. и 1999 г.); Международной конференции "Mathematics in applications", посвященной академику С. К. Годунову (Новосибирск, 1999 г.); Международной конференции "Синергетика. Самоорганизующиеся процессы в системах и технологиях" (Комсомольск-на-Амуре, 2000 г.). В целом работа докладывалась на семинарах "Математическое моделирование в механике" ИВМ СО РАН.

Исследования по теме диссертации проводились при финансовой поддержке Грантового центра при Новосибирском государственном университете (проект № 21, 1996-1997 гг.), Красноярского краевого фонда науки (грант 8F0022, 1999 г.) и Российского фонда фундаментальных исследований (гранты 97-01-00434, 99-01-00453, 01-01-06319, 01-01-00921).

Основные результаты опубликованы в работах [90 - 102].

Основные результаты диссертационной работы состоят в следующем:

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

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

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

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

ЗАКЛЮЧЕНИЕ

1. Джонсон К. Механика контактного взаимодействия. - М.: Мир, 1989. - 510 с.

2. Развитие теории контактных задач в СССР. М.: Наука, 1976. -494 с.

3. Signorini A. Sopra alcune questioni di elastostatica // Atti Soc. Ital. per Progresso delle Scienze. 1933. - P. 513-533.

4. Signorini A. Questioni di elasticita non linearizzata о semilinearizzata // Rend, di Matem. e delle sue appl. 1959. - V. 18.

5. Фикера Г. Теоремы существования в теории упругости. М.: Мир-, 1974. - 160 с.

6. Киндерлерер Д., Стампаккья Г. Введение в вариационные неравенства и их приложения. М.: Мир, 1983. - 256 с.

7. Экланд И., Темам Р. Выпуклый анализ и экстремальные задачи. -М.: Мир, 1979. 400 с.

8. Дюво Г., Лионе Ж.-Л. Неравенства в механике и физике. М.: Наука, 1980. - 383 с.

9. Главачек И., Гаслингер Я., Нечас И., Ловишек Я. Решение вариационных неравенств в механике. М.: Мир, 1986. - 270 с.

10. Байокки А., Капело А. Вариационные и квазивариационные неравенства. М.: Наука, 1988. - 488 с.

11. Панагиотопулос П. Неравенства в механике и их приложения. М.: Мир, 1989. - 492 с.

12. Бенсуссан А., Лионе Ж.-Л. Импульсное управление и квазивариационные неравенства. М.: Наука, 1987. - 597 с.

13. Кравчук А. С. К задаче Герца для линейно и нелинейно упругих тел конечных размеров // ДАН СССР. 1976. - Т. 230, № 2. - С. 308-310.

14. Кравчук А. С. Постановка задачи о контакте нескольких деформируемых тел как задачи нелинейного программирования // ПММ. -1978. Т. 42, вып. 3. - С. 466-474.

15. Кравчук А. С., Сурсяков В. А. Численное решение геометрически нелинейных контактных задач // ДАН СССР. 1981. - Т. 259, № 6. - С. 1327-1329.

16. Кравчук А. С. К теории контактных задач с учетом трения на поверхности соприкосновения // ПММ. 1980. - Т. 44, вып. 1. -С. 122-129.

17. Кравчук А. С. Вариационные и квазивариационные неравенства в механике. М.: МГАПИ, 1997. - 340 с.

18. Гловински Р., Лионе Ж.-JI., Тремольер Р. Численное исследование вариационных неравенств. М.: Мир, 1979. - 576 с.

19. Баничук Н. В., Петров В. М., Черноусько Ф. Л. Численное решение вариационных и краевых задач методом локальных вариаций // ЖВМиМФ. 1966. - Т. 6, № 6. - С. 947-961.

20. Картвелшлвили В. М. Численное решение двух контактных задач для упругих пластин // Изв. АН СССР. МТТ. 1976. - № 6. -С. 68-72.

21. Вовкушевский А. В., Шойхет Б. А. Расчет массивных гидротехнических сооружений с учетом раскрытия швов. М.: Энергоиздат, 1981. - 136 с.

22. Вовкушевский А. В. Вариационная постановка и методы решения контактной задачи с трением при учете шероховатости поверхностей // Изв. ВНИИ гидротехн. 1997. - Т. 232, № 1. - С. 151-160.

23. Спектор А. А. Вариационный метод исследования контактных задач с проскальзыванием и сцеплением // ДАН СССР. 1977. - Т. 236, № 1. - С. 39-42.

24. Гольдштейн Р. В., Спектор А. А. Вариационные оценки решений смешанных пространственных задач теории упругости с неизвестной границей // Изв. АН СССР. МТТ. 1978. - № 2. - С. 82-94.

25. Спектор А. А. Вариационные методы в пространственных задачах о нестационарном взаимодействии упругих тел с трением // ПММ. -1987. Т. 51, вып. 1. - С. 76-83.

26. Спектор А. А. Качественные свойства и поведение во времени решений пространственных контактных задач с трением // Докл. АН СССР. 1988. - Т. 301, № 6. - С. 1337-1340.

27. Гольдштейн Р. В., Спектор А. А. Вариационный метод исследования пространственных смешанных задач о плоском разрезе в упругой среде при наличии проскальзывания и сцепления его поверхностей // ПММ. 1983. - Т. 47, вып. 2. - С. 276-285.

28. Няшин Ю. И., Чернопазов С. А. Вариационный метод решения контактной задачи теории упругости с трением // ПММ. 1997. - Т. 61, № 4. - С. 692-702.

29. Dostal Z., Vondrak V. Duality based solution of contact problem with Coulomb friction // Arch. Mech. 1997. - V. 49, No. 3. - P. 453-460.

30. Дурнев В. А. Решение задачи Синьорини с несимметричным трением на границе // Тр. Псков, политехи, ин-та. 1998. - № 2. - С. 152-155.

31. Кузьменко В. И. О вариационном подходе в теории контактных задач для нелинейно-упругих тел // ПММ. 1979. - Т. 43, вып. 5. -С. 893-901.

32. Кузьменко В. И. О контактных задачах теории пластичности при сложном нагружении // ПММ. 1984. - Т. 48, вып. 3. - С. 473-481.

33. Кузьменко В. И. Контактные задачи для упругопластической полосы при сложном нагружении // Изв. АН СССР. МТТ. 1985. - № 6. -С. 128-135.

34. Кузьменко В. И. Контактные задачи теории пластичности с учетом трения на контактной поверхности // Трение и износ. 1987. - Т. 8, № 1. - С. 45-52.

35. Няшин Ю. И., Чернопазов С. А. К постановке контактной задачи упругопластичности // ПММ. 1989. - Т. 53, вып. 6. - С. 1023-1027.

36. Ильюшин А. А. Вопросы теории течения пластического вещества по поверхностям // ПММ. 1954. - Т. 18, вып. 3. - С. 265-288.

37. Григорян С. С. Новый закон трения и механизм крупномасштабных горных обвалов и оползней // ДАН СССР. 1979. - Т. 244, № 4. -С. 846-849.

38. Хилл Р. Математическая теория пластичности / Перев. с англ. М.: ИЛ, 1956. - 407 с.

39. Мосолов П. П., Мясников В. П. Механика жесткопластических сред. М.: Наука, 1981. - 208 с.

40. Shillor Meir, Sofonea Mircea. A quasistatic contact problem for an elasto-plastic rod // J. Math. Anal, and Appl. 1998. - V. 217, No. 2. - P. 579596.

41. Арутюнян H. X., Дроздов А. Д., Наумов В. Э. Механика растущих вязкоупругопластических тел. М.: Наука, 1987. - 471 с.

42. Арутюнян Н. X., Манжиров А. В., Наумов В. Э. Контактные задачи механики растущих тел. М.: Наука, 1991. - 176 с.

43. Спектор А. А. Стабилизация режимов трения при контактном взаимодействии движущегося упругого тела и основания // Трение и износ. 1990. - Т. 11, № 1. - С. 49-56.

44. Горячева И. Г., Добычин М. Н. Изнашивание неоднородно упрочненных поверхностей // Трение и износ. 1986. - Т. 77, № 6. - С. 985-992.

45. Александров В. М. Контактные задачи в трибологии // Механика и научно-технический прогресс. Т. 3. Механика деформируемого твердого тела. М.: Наука, 1988. - С. 170-180.

46. Александров В. М., Коваленко Е. В. Математические методы в контактных задачах с износом // Нелинейные модели и задачи механики деформируемого твердого тела. М.: Наука, 1984. - С. 77-89.

47. Львов Г. И. Вариационная постановка контактной задачи для линейно-упругих и физически нелинейных пологих оболочек // ПММ. -1982. Т. 46, вып. 5. - С. 841-846.

48. Хлуднев А. М. Об односторонней задаче, связанной со слабо нелинейным параболическим оператором // Сиб. мат. журн. 1978. - Т. 19, № 2. - С. 412-417.

49. Хлуднев А. М. Вариационный подход к проблеме контакта пологой оболочки с жестким телом // Дифференциальные уравнения с частными производными. Тр. семинара С.Л. Соболева. Новосибирск: ИМ СО АН СССР, 1981. - № 2. - С. 109-114.

50. Хлуднев А. М. О существовании решений в квазистатической задаче упругопластического деформирования оболочек // Сиб. мат. журн. -1984. Т. 25, № 5. - С. 168-176.

51. Хлуднев А. М. О вариационном неравенстве для оператора пологих оболочек с ограничением на границе // ПММ. 1987. - Т. 51, вып. 2.- С. 345-348.

52. Хлуднев А. М. О вариационных неравенствах в контактных пластических задачах // Дифференциальные уравнения. 1988. - Т. 24, № 9.- С. 1622-1628.

53. Хлуднев А. М. Оптимальное управление пластиной над препятствием // Сиб. мат. журн. 1990. - Т. 31, № 1. - С. 172-178.

54. Хлуднев А. М. Оптимальное управление в одномерных упругопластических моделях // ПМТФ. 1991. - № 5. - С. 112-115.

55. Хлуднев А. М. Об экстремальных формах разрезов в пластине // Изв. РАН. МТТ. 1992. - № 1. - С. 170-176.

56. Хлудпев А. М. Контактная задача для пологой оболочки с трещиной // ПММ. 1995. - Т. 59, вып. 2. - С. 318-326.

57. Хлуднев А. М. Задача о равновесии термоупругой пластины, содержащей трещину // Сиб. мат. журн. 1996. - Т. 37, № 2. - С. 452-463.

58. Хлуднев А. М. О контакте двух пластин, одна из которых содержит трещину // ПММ. 1997. - Т. 61, вып. 5. - С. 882-894.

59. Ковтуненко В. А. Метод численного решения задачи о контакте упругой пластины с препятствием // ПМТФ. 1994. - Т. 35, № 5.- С. 142-146.

60. Ковтуненко В. А. Итерационный метод решения вариационных неравенств в контактной упругопластической задаче с использованием метода штрафов // ЖВМиМФ. 1993. - Т. 33, № 9. - С. 1409-1415. ■

61. Ковтуненко В. А. Численное решение задачи о контакте упругопластической балки для модели Тимошенко // Изв. РАН. МТТ. 1996. -№ 5. - С. 79-84.

62. Ковтуненко В. А. Решение задачи о балке с разрезом // ПМТФ. -1996. Т. 37, № 4. - С. 160-166.

63. Галин JI. А. Контактные задачи теории упругости. М.: Гостех-теоретиздат, 1953. - 264 с.

64. Галин JI. А. Контактные задачи теории упругости и вязкоупругости.- М.: Наука, 1980. 304 с.

65. Гулидов А. И., Фомин В. М. Численное моделирование отскока осе-симметричных стержней от твердой преграды // ПМТФ. 1980. -№ 3. - С. 126-132.

66. Ковалев В. Г., Косенков В. М. Динамическое контактное взаимодействие упругопластических цилиндров под действием импульсной нагрузки // Изв. РАН. МТТ. 1997. - № 1. - С. 84-92.

67. Коваленко Е. В., Зеленцов В. Б. Асимптотические методы в нестационарных динамических контактных задачах для упругого полупространства // ПМТФ. 1997. - Т. 38, № 1. - С. 111-119.

68. Кондауров В. И., Петров И. Б., Холодов А. С. Численное моделирование процесса внедрения жесткого тела вращения в упругопласти-ческую преграду // ПМТФ. 1984. - № 4. - С. 132-139.

69. Хорев И. Е., Горельский В. А., Зелепугин С. А., Толкачев В. Ф. Исследование деформирования и кинетики разрушения контактирующих тел при несимметричном динамическом взаимодействии // ФГВ.- 1983. Т. 19, № 5. - С. 119-123.

70. Горельский В. А., Хорев И. Е., Югов Н. Т. Динамика трехмерного процесса несимметричного взаимодействия деформируемых тел с жесткой стенкой // ПМТФ. 1985. - № 4. - С. 112-118.

71. Горельский В. А., Зелепугин С. А., Сидоров В. Н. Численный анализ соударения разноплотных тел при ударе под углом // Изв. РАН. МТТ.- 1999. № 3. - С. 45-54.

72. Кошур В. Д., Немировский Ю. В. Континуальные и дискретные модели динамического деформирования элементов конструкций. Новосибирск: Наука, 1990. - 198 с.

73. Гулидов А. И., Шабалин И. И. Моделирование разрушенного материала дискретными частицами конечного размера // ПМТФ. 1997.- Т. 38, № 3. С. 14-19.

74. Киселев А. Б. О критерии динамического разрушения при ударном взаимодействии упругопластических тел // Вестн. МГУ. Сер. 1. Математика, механика. 1986. - № 6. - С. 46-51.

75. Алехин В. В., Коробейников С. Н. Алгоритм решения трехмерных контактных задач методом конечных элементов // Численные методы решения задач теории упругости и пластичности: Тр. 13-й Межрес-публ. конф. / Новосибирск, ИТПМ СО РАН, 1995. С. 4-12.

76. Егунов В. В., Конюхов А. Б. Численное моделирование нестационарного контактного взаимодействия упругопластических тел на основе МКЭ, неявных схем интегрирования и метода штрафных функций // Прикл. пробл. прочн. и пластич. 1997. - № 55. - С. 120-125.

77. Кравчук А. С. Применение методов и алгоритмов теории оптимизации в контактных задачах механики деформируемого твердого тела // Прикл. пробл. прочн. и пластич. 1997. - № 57. - С. 109-115.

78. Annin В. D., Sadovsky V. М. A numerical analysis of laminated elastic-plastic plates under dynamic loading // Composites Science and Technology. 1992. - V. 45. - P. 241-246.

79. Аннин Б. Д., Садовский В. М. О численной реализации вариационного неравенства в задачах динамики упругопластических тел // ЖВМиМФ. 1996. - Т. 36, № 9. - С. 177-191.

80. Садовский В. М. Разрывные решения в задачах динамики упругопластических сред. М.: Наука, 1997. - 208 с.

81. Szefer G. Dynamic contact of bodies experiencing large deformations // Acta Mech. 1997. - V. 125, No. 1. - P. 217-233.

82. Каменецкий В. Ф., Семенов А. Ю. Самосогласованное выделение разрывов при сквозных расчетах газодинамических течений // ЖВМиМФ. 1994. - Т. 34, № 10. - С. 1489-1502.

83. Фрязинов И. В. Экономичные симметризованные схемы решения краевых задач для многомерного уравнения параболического типа // ЖВМиМФ. 1968. - Т. 8, № 2. - С. 436-443.

84. Уилкинс М. JI. Расчет упругопластических течений // Вычислительные методы в гидродинамике. М.: Мир, 1967. - С. 212-263.

85. Бычек (Садовская) О. В. Численное моделирование динамического контактного взаимодействия в пространственной постановке // Материалы XXXIV Междунар. научной студ. конф. "Студент и научно-технический прогресс": Математика. Новосибирск, НГУ, 1996. -С. 13.

86. Бычек (Садовская) О. В. Моделирование динамического контактного взаимодействия // Труды семинара "Математическое моделирование в механике" ВЦК СО РАН. Деп. в ВИНИТИ. - № 446 - 1397. -С. 41-51.

87. Бычек (Садовская) О. В. Моделирование динамического контактного взаимодействия нелинейно-упругих тел // Материалы XXXV Между-нар. научной студ. конф. "Студент и научно-технический прогресс": Математика. Новосибирск, НГУ, 1997. - С. 19-20.

88. Бычек (Садовская) О. В. Моделирование динамического контактного взаимодействия нелинейно-упругих тел // Труды семинара "Математическое моделирование в механике" ВЦК СО РАН. Деп. в ВИНИТИ. - N°. 3357 - В97. - С. 76-90.

89. Бычек (Садовская) О. В., Садовский В. М. О численном исследовании динамического контактного взаимодействия деформируемых тел // Тезисы докладов Междунар. конф. "Математические модели и методы их исследования". Красноярск, 1997. - С. 50.

90. Бычек (Садовская) О. В. Моделирование динамического контактного взаимодействия нелинейно-упругих тел // Тезисы докладов конференции молодых ученых КНЦ СО РАН. Красноярск, 1997. - С. 15-16.

91. Бычек (Садовская) О. В., Садовский В. М. К исследованию динамического контактного взаимодействия деформируемых тел // ПМТФ.- 1998. Т. 39, № 4. - С. 167-173.

92. Садовская О. В. О численной реализации условий динамического контактного взаимодействия с учетом трения // Материалы Всерос. семинара "Теория сеточных методов для нелинейных краевых задач".- Казань: Унипресс, 1998. С. 63-64.

93. Садовская О. В. О численном исследовании соударения упругоплас-тических тел с учетом конечных поворотов // Динамика сплошной среды. Сб. науч. трудов. Вып. 114. Математические проблемы механики сплошных сред. Новосибирск, 1999. - С. 196-199.

94. Садовская О. В. Численное решение динамических контактных задач // Материалы конференции молодых ученых ИВМ СО РАН. -Красноярск, 1999. С. 42-54.

95. Annin В. D., Sadovskaya О. V., Sadovskii V. M. Variational approach to dynamic contact problems of elastoplasticity // International Conf. "Mathematics in applications" (abstracts). Novosibirsk, 1999. -P. 18-19.

96. Аннин Б. Д., Садовская О. В., Садовский В. М. Численное моделирование косого соударения пластин в упругопластической постановке // Физическая мезомеханика. 2000. - Т. 3, № 4. - С. 23-28.

97. Карманов В. Г. Математическое программирование. М.: Наука, 1986. - 288 с.

98. Садовский В. М. Методы решения вариационных задач механики.- Новосибирск: Изд-во СО РАН, 1998. 184 с.

99. Годунов С. К. Разностный метод численного расчета разрывных решений уравнений гидродинамики // Матем. сб. 1959. - Т. 47, вып. 3. - С. 271-306.

100. Годунов С. К., Забродин А. В., Иванов М. Я. и др. Численное решение многомерных задач газовой динамики. М.: Наука, 1976. -400 с.

101. Колган В. П. Применение принципа минимальных значений производной к построению конечно-разностных схем для расчета разрывных решений газовой динамики // Уч. зап. ЦАГИ. 1972. - Т. 3, № 6.- С. 68-77.

102. Родионов А. В. Монотонная схема второго порядка аппроксимации для сквозного расчета неравновесных течений // ЖВМиМФ. 1987.- Т. 27, № 4. С. 585-593.

103. Родионов А. В. Повышение порядка аппроксимации схемы С. К. Годунова // ЖВМиМФ. 1987. - Т. 27, № 12. - С. 1853-1860.

104. Малышев А. П. Монотонная разностная схема повышенной точности для численного моделирования волновых процессов // ЖВМиМФ.- 1996. Т. 36, № 9. - С. 155-158.

105. Колмогоров В. Л. Напряжения. Деформации. Разрушение. М.: Металлургия, 1970. - 230 с.

106. Колмогоров В. Л. Механика обработки металлов давлением. М.: Металлургия, 1986. - 688 с.

107. Тиме И. А. Сопротивление металлов и дерева резанию. Санкт-Петербург, 1870.

108. Зворыкин К. А. Работа и усилие, необходимые для отделения металлических стружек. Санкт-Петербург, 1893.

109. Брике А. А. Резание металлов (строгание). Санкт-Петербург, 1896.

110. Усачев Я. Г. Явления, происходящие при резании металлов // Изв. Петроградского политехи, ин-та. 1915. - Вып. 23.

111. Клушин М. И. Резание металлов. М.: Машгиз, 1958. - 363 с.

112. Лоладзе Т. Н. Стружкообразование при резании металлов. М.: Машгиз, 1952. 200 с.

113. Лоладзе Т. Н. Прочность и износостойкость режущего инструмента. М.: Машиностроение, 1982. - 320 с.

114. Туркович Б. Ф. О касательных напряжениях при резании металлов // Конструирование и технология машиностроения. 1970. — № 1. — С. 154-161.

115. Блэк У. Модель напряжения пластического течения при резании металлов // Конструирование и технология машиностроения. 1979.- № 4. С. 124-139.

116. Розенберг А. М., Розенберг О. А. Механика пластического деформирования в процессах резания и деформирующего протягивания. -Киев: Наукова думка, 1990. 320 с.

117. Розенберг А. М., Еремин А. Н. Элементы теории процесса резания металлов. М.: Машгиз, 1956. - 319 с.

118. Зорев Н. Н. Вопросы механики процесса резания металлов. М.: Машгиз, 1956. - 367 с.

119. Бобров В. Ф., Сидельников А. И. Особенности образования суставчатой и элементной стружек при высокой скорости резания // Вестник машиностроения. 1976. - № 7. - С. 61-66.

120. Чижов В. Н., Михайлов С. В. Исследование закономерностей завивания сливной стружки // Сб.: Оптимизация процессов резания жаро-и особопрочных материалов. Уфа, 1985. - С. 70-73.

121. Кабалдин Ю. Г., Олейников А. И., Шпилев А. М., Бурков А. А. Математическое моделирование самоорганизующихся процессов в технологических системах обработки резанием. Владивосток: Дальнаука, 2000. - 195 с.

122. Дерибас А. А., Кудинов В. М., Матвеенков Ф. И., Симонов В. А. Сварка взрывом // ФГВ. 1967. - Т. 3, № 1. - С. 111-118.

123. Allen W. A., Mapes J. М., Wilson W. G. An effect produced by oblique impact of a cylinder on a thin target //J. Appl. Phys. 1954. - V. 25, No. 5. - P. 675-676.

124. Абрахамсон Г. P. Остаточные периодические деформации поверхности под действием перемещающейся струи // Тр. Амер. об-ва инж.-мех. Сер. Е. Прикл. механика. 1961. - Т. 28, № 4. - С. 45-55.

125. Bahrani A. S., Black Т. J., Crossland В. The mechanics of wave formation in explosive welding // Proc. Roy. Soc. Ser. A. 1967. - V. 296, No. 1445. - P. 123-136.

126. Hunt J. H. Wave formation in explosive welding // The Philosoph. Mag.- 1968. V. 17, No. 148. - P. 669-680.

127. Годунов С. К., Дерибас А. А., Козин Н. С. Волнообразование при сварке взрывом // ПМТФ. 1971. № 3. - С. 63-72.

128. Дерибас А. А., Кудинов В. М., Матвеенков Ф. И. Влияние начальных параметров на процесс волнообразования при сварке металлов взрывом // ФГВ. 1967. - Т. 3, № 4. - С. 561-568.

129. Дерибас А. А., Кудинов В. М., Матвеенков Ф. И., Симонов В. А. О моделировании процесса волнообразования при сварке взрывом // ФГВ. 1968. - Т. 4, № 1. - С. 100-107.

130. Захаренко И. Д. Сварка металлов взрывом. Минск: Навука i тэх-шка, 1990. - 205 с.

131. Симонов В. А. Области сварки взрывом. Основные параметры и критерии. Новосибирск, ин-т гидродинамики им. М. А. Лаврентьева СО РАН, 1995. - 61 с.

132. Лаврентьев М. А., Ишлинский А. Ю. Динамические формы потери устойчивости упругих систем // Докл. АН СССР. 1949. - Т. 64, № 6.- С. 779-782.

133. Корнев В. М., Яковлев И. В. Модель волнообразования при сварке взрывом // ФГВ. 1984. - Т. 20, № 2. - С. 87-90.

134. Волнообразование при косых соударениях: Сб. статей. Новосибирск: Изд-во института дискрет, матем. и информатики, 2000. -221 с.1. РОССИЙСКАЯ ГОСУДАРСТВЕННА1. БИБЛИОТР.ГЛ'