авторефераты диссертаций БЕСПЛАТНАЯ  БИБЛИОТЕКА

АВТОРЕФЕРАТЫ КАНДИДАТСКИХ, ДОКТОРСКИХ ДИССЕРТАЦИЙ

<< ГЛАВНАЯ
АГРОИНЖЕНЕРИЯ
АСТРОНОМИЯ
БЕЗОПАСНОСТЬ
БИОЛОГИЯ
ЗЕМЛЯ
ИНФОРМАТИКА
ИСКУССТВОВЕДЕНИЕ
ИСТОРИЯ
КУЛЬТУРОЛОГИЯ
МАШИНОСТРОЕНИЕ
МЕДИЦИНА
МЕТАЛЛУРГИЯ
МЕХАНИКА
ПЕДАГОГИКА
ПОЛИТИКА
ПРИБОРОСТРОЕНИЕ
ПРОДОВОЛЬСТВИЕ
ПСИХОЛОГИЯ
РАДИОТЕХНИКА
СЕЛЬСКОЕ ХОЗЯЙСТВО
СОЦИОЛОГИЯ
СТРОИТЕЛЬСТВО
ТЕХНИЧЕСКИЕ НАУКИ
ТРАНСПОРТ
ФАРМАЦЕВТИКА
ФИЗИКА
ФИЗИОЛОГИЯ
ФИЛОЛОГИЯ
ФИЛОСОФИЯ
ХИМИЯ
ЭКОНОМИКА
ЭЛЕКТРОТЕХНИКА
ЭНЕРГЕТИКА
ЮРИСПРУДЕНЦИЯ
ЯЗЫКОЗНАНИЕ
РАЗНОЕ
КОНТАКТЫ

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

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ

имени М.В. ЛОМОНОСОВА

МЕХАНИКО-МАТЕМАТИЧЕСКИЙ ФАКУЛЬТЕТ

На правах рукописи

УДК 539.3

Сержкин Алексей Александрович

ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ СОУДАРЕНИЯ

ЦИЛИНДРА С НЕДЕФОРМИРУЕМОЙ ПРЕГРАДОЙ

МЕТОДОМ РАЗДЕЛЕНИЯ ПО ФИЗИЧЕСКИМ

ПРОЦЕССАМ НА ПОДВИЖНЫХ ЭЙЛЕРОВЫХ

СЕТКАХ

Специальность 01.02.04 – механика деформируемого твердого тела

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Научный руководитель:

доктор физико математических наук, профессор Киселев Алексей Борисович Москва -

Работа выполнена на кафедре газовой и волновой динамики механико математического факультета Московского государственного университета имени М.В. Ломоносова и в Центре фундаментальных и прикладных исследований ВНИИ автоматики имени Н.Л. Духова.

Научный руководитель: Доктор физико-математических наук, профессор Киселев Алексей Борисович

Официальные оппоненты: Доктор физико-математических наук, старший научный сотрудник Пшеничнов Сергей Геннадиевич Доктор физико-математических наук, старший научный сотрудник Рыбакин Борис Петрович

Ведущая организация: Московский государственный Университет леса

Защита состоится «17» мая 2013 г. в 16 часов на заседании диссертационного совета Д 501.001.91 при Московском государственном университете имени М.В. Ломоносова по адресу: 119992, Москва, Ленинские горы, Главное здание МГУ, Механико-математический факультет, аудитория 16-10.

С диссертацией можно ознакомиться в библиотеке механико математического факультета МГУ (Главное здание, 14 этаж) Автореферат разослан «15» апреля 2013 г.

Ученый секретарь диссертационного совета Д 501.001.91, доктор физико-математических наук, профессор С.В. Шешенин

Общая характеристика работы

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

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

Исследование характера деформирования стержня в зоне 1.

ударяемого торца, а также деформирования самой плоскости торца.

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

Продолжительность взаимодействия стержня и преграды. Анализ 2.

критериев отскока.

Влияние трения между стержнем и преградой на конечную 3.

форму стержня.

Особенности и причины образования ступенчатой формы 4.

боковой поверхности стержня.

Определение основных закономерностей разрушения стержня.

5.

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

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

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

Работа выполнена при поддержке РФФИ (гранты № 09-01-00144а, № 12-01-00425а), а также программы «Взаимодействие жидких, твердых и газообразных сред" (лаборатория динамики деформируемых сред механико математического факультета МГУ имени М.В. Ломоносова;

номер госрегистрации № 01.9.80 005453).

Апробация работы.

Основные результаты диссертации докладывались на следующих научных конференциях и семинарах:

Ломоносовские чтения МГУ. Москва, ноябрь 2011.

1.

Ломоносовские чтения МГУ. Москва, апрель 2012.

2.

Ломоносовские чтения МГУ. Москва, апрель 2013.

3.

Международная научная конференция студентов, аспирантов и 4.

молодых учных «Ломоносов-2013», апрель 2013.

XI Забабахинские научные чтения. Снежинск. 16-20 апреля 5.

2012.

6. Advanced Problems in Mechanics. St. Petersburg. July 2-8, 2012.

7. European Congress on Computational Methods in Applying Science and Engineering (ECCOMAS 2012). September 10-14, 2012.

V-VII научно-технические конференциии молодых ученых 8.

ВНИИ автоматики имени Н.Л. Духова (март 2011, март 2012, март 2013).

Научно- исследовательский семинар кафедры газовой и волновой 9.

динамики механико-математического факультета МГУ под рук. академика РАН Р.И. Нигматулина.

10. Научно-исследовательский семинар кафедры механики композитов под рук. профессора Б.Е. Победри.

11. Научно-исследовательский семинар кафедры теории пластичности механико-математического факультета МГУ под рук. член корр. РАН Е.В. Ломакина.

12. Научно-исследовательский семинар кафедры теории упругости механико-математического факультета МГУ под рук. профессора И.А.

Кийко.

Публикации. По работе имеется 5 публикаций, в том числе две статьи в журналах из перечня ВАК.

Личный вклад автора состоял в участии в разработке численного метода и создании программного комплекса «ТИС», участия в постановке задачи соударения цилиндра с недеформируемой преградой, в адаптации комплекса для расчета рассматриваемой задачи, в проведении расчетов и анализе их результатов.

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

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

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

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

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

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

( vz ) ( vr ) vr t z r r ( vz ) ( vz ) ( vr vz ) vr vz zz zr zr t z r r z r r vr ( vr ) ( vr ) ( vr vz ) rz rr rr 0 (1) t r z r z r r ( e) ( ev z ) ( evr ) evr ( zz vz rz vr ) t z r r z ( zr vz rr vr ) rz vz vr v rr r r r а также уравнений упругопластического течения Прандтля-Рейса:

S zz S zz S zz S zz v r vz vr vz vr S zz t z r r r z vz vr vr vz 2 S zz 3 z r r z S rr S rr S rr S rr vr vz vr vz vr S zr t z r r r z 2 vz vr vr vr (2) 2 S rr 3 z r r r Srz Srz S rz Srz vr S zz Srr vz vr vz vr t z r r 2 r z vz vr Srz r z где 2Y 3W,Q 2Y 2 (3) 2Y 0, Q vr vr vz vz vz W 2 S rr S rz S zz (4) r z r r z Здесь используются следующие обозначения: - модуль сдвига, Y предел текучести материала при простом растяжении, плотность, vz2 vr (vz, vr ) - вектор скорости, e - плотность полной энергии v (на единицу массы), плотность внутренней энергии,,,, zz rr zr компоненты тензора напряжений, который раскладывается на шаровую и девиаторную чаcти p S zz, p Srr, zz rr pS, S zr, p ( ) / 3.

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

vr vz vr B H * * t z r r 1 H( ) H( ) 0 (5) 4 0 Y0 ln Y0 ln (6), 3 vr Su Su * * vz vr C Su H Su 0 (7) t z r r (1 )(1 ) (1 )(1 ) Система (1-7) замыкается с помощью уравнения для шаровой части тензора напряжений:

v p K ln ( ) B ln(1 ) (8) c 0 Здесь K объемный модуль материала, коэффициент объемного V расширения, теплоемкость при постоянных деформациях, c начальная плотность материала.

В качестве критерия начала макроразрушения используется энтропийный критерий достижения удельной диссипацией предельного значения (Киселев А.Б., Юмашев М.В. Деформирование и разрушение при ударном нагружении. Модель повреждаемой термоупругопластической среды // ПМТФ. – 1990. - № 5).

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

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

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

Для простоты рассматривается система уравнений без учета параметров поврежденности (1-2). Она может быть записана в следующей, дивергентой форме:

fk q H (9) t xk vk vk vz vz kz vk vr vr kr fk vk E v q E kj j Где - вектор состояния среды, а - вектор потока S zz S zz vk S rz S rz vk S rr S rr vk в направлении координаты k, k z, r. Правая часть системы (9) z, r, j представляет собой недифференциальные члены уравнений системы (1)-(7).

vr r vz vr rz r r vr rr r r r vr e rz v z vr rr vr H (10) r r S zz vr vz vr 2 vz vr vr vz S zz 2 S r r z 3 z r r z S rz vr S zz S rr vz vr vz vr S rz r 2 r z r z S rr vr vz vr 2 vz vr vr vr S zr 2 S rr r r z 3 z r r r Основная идея численного метода состоит в расщеплении определяющих соотношений математической модели в соответствии с основными физическими процессами. Использование этого подхода сводит решение к решению набора более простых задач без потери точности. А именно, расщепим решение системы (9) на два шага, условно назовем их «гидродинамический» и «упругопластический». На «гидродинамическом»

шаге будем решать следующую систему уравнений:

fk q 0 (11) t xk Она представляет систему уравнений гидродинамики с учетом наличия фиксированных в каждой лагранжевой точке значений девиатора тензора напряжений. Ее решение строится на подвижной эйлеровой расчетной сетке с помощью метода конечного объема. Аппроксимация численных потоков производится по методу Русанова, а так же, с целью повышения точности схемы, используется метод подсеточного линейного восполнения MUSCL.

На втором шаге будем решать систему уравнений:

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

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

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

Для проверки метода приводится сравнение полученных численных и экспериментальных результатов для стержней из алюминия l0 0.04694 м, r0 0.00391м, V0 275м / с и l0 0.02347 м, r0 0.00391м, 373м / с. (Wilkins M.L., Guinan M.W. Impact of cylinders on a rigid V boundary // J. Appl. Phys. 1973. V. 44. № 3. Pp. 1200-1216.). Из рис.2 видно, что в момент окончания пластического деформирования размеры и форма стержня, полученные численно, совпадают с размерами и формой стержня, полученными в результате эксперимента.

Кроме того, наблюдается сходство с численным решением данной задачи, полученным с использованием метода, основанного на эйлеровом подходе с использованием схемы ENO (Essentially Non-Oscillatory).

(Udaykumar H.S., Tran L., Belk D.M., Vanden K.J. An Eulerian method for computation of multimaterial impact with ENO shock-capturing and sharp interfaces // Journal of Computational Physics. 2003. V. 186. Pp. 136–177).

Рис. 2. Сравнение численных результатов (сверху) и экспериментальных данных (снизу).

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

Рис. 3. Контуры конечной формы стержня при расчетах с различными коэффициентами трения.

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

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

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

Также в данном разделе приводится анализ длительности процесса соударения. Установлено, что момент отскока стержня от преграды совпадает с моментом достижения кинетической энергией своего первого локального максимума после минимального значения (рис. 6). Зависимость времени отскока от скорости соударения имеет ступенчатый вид.

Ступенчатость вызвана особенностями взаимодействия упругих и пластических волн в осевом и радиальном направлениях.

Рис. 6. Зависимость безразмерных значений времени отскока (красная кривая), времени достижения кинетической энергией минимального значения (зеленая кривая) и следующего за ним локального максимума (синяя кривая) от скорости соударения В данном разделе также приводится следующее исследование взаимодействия стержня с преградой. Отслеживались площадь области контакта и сила взаимодействия. Было рассмотрено три варианта начальной скорости стержня - v0 20 м / с, v0 60 м / с, v0 273 м / с. Таким значениям скоростей будут соответствовать различные режимы деформирования. В первом случае, когда v0 20 м / с будет наблюдаться чисто упругий режим.

60 м / с часть стержня деформируется в Во втором случае, когда v упругопластическом режиме, но интенсивность пластических деформаций невысока. И, наконец, при скорости v0 273 м / с стержень будет деформироваться в упругопластическом режиме, причем наблюдается довольно высокая интенсивность пластических деформаций. На рис. приведены указанные зависимости.

Рис. 7. Величина силы взаимодействия (слева) и площади контакта стержня с преградой (справа) для начальных значений.

20 м / с (синяя кривая), v0 60 м / с (красная кривая), v 273 м / с (серая кривая).

v Как видно, при высоких скоростях взаимодействия происходит частичный отрыв поверхности ударяемого торца от преграды. Период колебаний силы в начале взаимодействия соответствует удвоенному времени прохода ударной волны от оси симметрии до свободной поверхности. Кроме того отслеживалось максимальное значение скорости точки отрыва. В случае скорости соударения v0 20 м / с, когда наблюдается упругий режим деформирования, максимальная скорость точки отрыва составляет vb 10526 м / с, что почти вдвое превышает скорость звука для алюминия c 5330 м / с. Это связано с тем, что деформации малы, и фронт волны разгрузки в продольном направлении сильно не искажается, а остается практически плоским, и отрыв всех точек поверхности торца происходит почти одновременно.

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

273 м / с, максимальная скорость точки В третьем случае, когда v 4353 м / с. При этом данная отрыва получается значительно ниже, vb величина достигается в момент, когда происходит окончательный отрыв стержня. Рассматриваемые скорости в случаях частичного отрыва значительно ниже.

В четвертой главе приводятся численные решения задачи с учетом повреждаемости.

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

(Киселев А.Б., Юмашев М.В. Деформирование и разрушение при ударном нагружении. Модель повреждаемой термоупругопластической среды // ПМТФ. – 1990. - № 5). Был рассмотрен случай столкновения стержня с начальными размерами l 0 0.04694 м, r0 0.00782 м со скоростью 750 м / с. В этом случае критерий разрушения выполняется v одновременно в двух областях в районе преграды - около оси симметрии и в районе свободной поверхности (рис. 8).

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

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

Математическое моделирование фрагментации тонкостенных конструкций и компактных элементов при взрывном нагружении и ударном взаимодействии // Математическое моделирование. - 2012. - Т. 24. - № 2. – С. 33-66).

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

2 Rk ( D*) N (13) Здесь R - средний радиус торообразной области «лепесткового»

разрушения, коэффициент 0 k 1 определяет часть запасенной энергии, идущую на образование свободной поверхности, обычно его значение k 0.5, - плотность внутренней энергии в разрушаемой области, D * величина предельной удельной диссипации, - количество энергии, необходимое для образования единицы свободной поверхности (так называемая ударная вязкость материала);

[Z] – целая часть числа Z.

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

NAv r 3 V (14) NAS r 2 kU Коэффициенты Av и As определяют пропорциональность объема и площади поверхности фрагмента соответствующей степени характерного радиуса, их отношение является функцией формы. Для шарообразных Av фрагментов эта функция максимальна,, для кубических фрагментов As Av Av 1, для тетраэдрических. Здесь U - упругая ( D)V As 6 As энергия, запасенная в области объемом V.

Для выбранного случая начальных условий получается, что число "лепестков" равно 9, а число фрагментов около оси симметрии Av 10 3. Для шарообразных фрагментов эта величина порядка N As 1000.

Основные результаты и выводы Метод разделения по физическим процессам позволяет с 1.

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

При увеличении силы трения между стержнем и преградой 2.

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

Конечная форма боковой поверхности стержней при 3.

пластическом деформировании имеет ступенчатый вид, характеризуемый возникновением «вторичной области деформации». Факт образования которой определяется соотношением длины и радиуса стержня и не зависит от скорости соударения.

Критерием отскока, или окончанием процесса пластического 4.

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

Зависимость продолжительности контакта с преградой от 5.

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

В процессе деформирования возможен частичный отрыв 6.

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

При больших скоростях соударения разрушение стержня может 7.

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

Список публикаций Киселев А.Б., Серёжкин А.А. Численное исследование 1.

особенностей процесса соударения упругопластического цилиндра с недеформируемой преградой // ПММ. – 2013. – Т. 77.

Меньшов И.С., Мищенко А.В., Серёжкин А.А. Численное 2.

моделирование упругопластических течений методом Годунова на подвижных эйлеровых сетках. // Математическое моделирование. – 2013. – Т.

25.

3. Menshov I., Mischenko A., Serezhkin A. An eulerian Godunov-type scheme for calculation of the elastic-plastic flow equations with moving grids // Europ. Congress on Comput. Methods in Appl. Sc. and Eng. (ECCOMAS 2012). J.

Eberhardsteiner et. al. (eds.). Vienna, Austria, September 10-14, 2012. CD format, 2012. Article 2164. 20 p.

Мищенко А.В., Серёжкин А.А., Меньшов И.С., Киселев А.Б. Метод 4.

разделения по физическим процессам для моделирования деформирования и разрушения твердых тел. // Забабахинские научные чтения: сб. материалов XI Межд. конф. 16-20 апреля 2012. – Снежинск: Изд-во РФЯЦ-ВНИИТФ, 2012. С. 306.

5. Киселев А.Б., Меньшов И.С., Серёжкин А.А. Вычислительный комплекс «ТИС»: тестирование на задачах гидродинамики // Ломоносовские чтения. Научная конф. Секция механики. Апрель 2012 года. Тезисы докладов. - М.: Изд-во Моск. ун-та, 2012. - С. 90-91.



 




 
2013 www.netess.ru - «Бесплатная библиотека авторефератов кандидатских и докторских диссертаций»

Материалы этого сайта размещены для ознакомления, все права принадлежат их авторам.
Если Вы не согласны с тем, что Ваш материал размещён на этом сайте, пожалуйста, напишите нам, мы в течении 1-2 рабочих дней удалим его.