CC BY
141
^^МЕТОДИЧЕСКИЕ ОСНОВЫ
СОВЕРШЕНСТВОВАНИЯ ПРОЕКТИРОВАНИЯ И ПРОИЗВОДСТВА ТЕХНИЧЕСКИХ СИСТЕМ
ПРИМЕНЕНИЕ СЛО/СЛЕТЕХНОЛОГИЙ ДЛЯ ИССЛЕДОВАНИЯ РАБОТОСПОСОБНОСТИ ЭЛЛАСТИЧНОГО ОБТЮРАТОРА В УСЛОВИЯХ ИМПУЛЬСНОГО НАГРУЖЕНИЯ ВЫСОКИМ ДАВЛЕНИЕМ
Г.В. Лепеш1, Е.С. Иванова2, К.А.Егоров3
1 Санкт-Петербургский государственный экономический университет (СПбГ-
ЭУ),191023, Санкт-Петербург, ул. Садовая, д. 21 2'ъАкционерное общество «Центральный научно-исследовательский институт материалов», (АО «ЦНИИМ»), 191014, Санкт-Петербург, Ул. Парадная, 8
Проведена сравнительная расчетная оценка применимости моделей Ми-Грюнайзена и Муни-Ривлина при исследовании напряжённо-деформированного состояния обтюрирующих устройств изделий в условиях высокоскоростного импульсного нагружения. Расчеты проведены в явных Лагранжевых координатах в программном пакете ANSYS с использованием решателя Autodyn.
Ключевые слова: обтюрирующие устройства, газодинамические импульсные устройства, гиперупругие материалы, динамика нагружения, пластическая деформация, уравнение состояния Ми-Грюнайзена, модель Муни-Ривлина, напряжённо-деформированное состояние.
APPLICATION OF CAD/CAE TECHNOLOGIES FOR RESEARCH OF OPERABILITY OF THE ELASTIC OBTURATOR IN CONDITIONS OF SHOCK LOADING BY THE HIGH PRESSURE
G.V. Lepesh, E.S. Ivanova, K. A. Egorov
Sankt Petersburg State University of Economics (SPbGEU), 191023, St. Petersburg, Sadovaya, 21 Joint stock company «Central scientific-research institute of materials»,
191014, St. Petersburg, Paradnaya St.,8
It was shown the comparative estimation of Mei-Gruneisen and Mooney-Rivlin models usage adequacy by obturating devices stress-strain state researching in case of high speed pulse loading. Calculations are carried out in Lagrangian explicit coordinates by Ansys program with Autodyn solver use .
Keywords: the obturating devices, gasdynamic pulse devices, hyper elastic materials, dynamics of loading, plastic deformation, the equation of a condition of Mi-Gryunaizena, Mooney-Rivlin's model, the intense deformed state.
Одной из основных причин разгарно-эрозионного действия при работе газодинамических импульсных устройств (ГИУ) является прорыв газов между трубой и обтюрирующими устройствами изделия. Практика работы газодинамических импульсных устройств показы-
вает, что появление прорыва газов вызывает повышенный износ трубы ГИУ на том участке, где имеет место прорыв. При значительном износе происходит недопустимое уменьшение максимального давления газов и начальной скорости изделия.
1 Лепеш Григорий Васильевич - доктор технических наук, профессор, заведующий кафедрой Машины и оборудование бытового и жилищно-коммунального назначения, СПбГЭУ, тел.:+7 921 751 2829,e-mail: gregoryl@yandex. ru;
2Иванова Елена Сергеевна - кандидат технических наук, главный конструктор по направлению АО «ЦНИИМ», тел.:+7 921 3608546, e-mail: ilena-bgtu@yandex.ru;
3Егоров Константин Александрович - научный сотрудник АО «ЦНИИМ», тел.:+7 921 9825455, e-mail: konstant55 @rambler.ru
В процессе проектирования ГИУ проводится тщательная экспериментальная отработка обтюрирующих устройств изделий, которая требует учета множества взаимосвязанных факторов, происходящих в сложных напряженных условиях их функционирования. В настоящее время в качестве эксперимента прежде всего предпочтение отдают численному эксперименту, проводимому на базе популярных CAD/CAE систем, таких как Autodesk, SolidWorks, Ansys и др. [1].
Популярность этих пакетов для решения задач оценки функционирования изделий в сложных динамических условиях обусловлена предоставляемой ими возможностью средств (собственных баз и решателей), обеспечивающих возможность исследования линейных и нелинейных процессов, в том числе происходящих в средах с различными физическими свойствами. При этом задача исследователя заключается в определении границ расчетных областей с использованием CAD системы, в выборе уравнения состояния - модели поведения материалов с различными свойствами в различных условиях, выборе подходящего решателя и определения граничных для расчетной области и начальных условий для происходящих в ней процессов.
Процесс функционирования обтюрирующих устройств происходит в условиях высокой динамики нагружения и высоких давлений. При этом элементы обтюрирующих устройств, как правило, изготавливаемые из полиамидных гиперупругих материалов, испытывают конечные пластические деформации - изменяют форму в процессе выполнения своей функции при относительно высокой скорости деформирования. Для описания процессов деформирования и разрушения в подобных условиях применяют различные физические модели, учитывающие условия нагружения и поведение сплошной деформируемой среды в этих усло-виях[2]. К возможным для данного случая следует отнести модель Ми-Грюнайзена, учитывающую динамику процесса, а также модели Муни-Ривлина, Нео-гука и Огдена, применяемые традиционно для гиперупругих материалов, испытывающих конечные деформации и высокие скорости деформирования и объемные давления. В данной статье оценивается применимость моделиМи-Грюнайзена [3] дополненной уравнением вязкопластичного поведения материала - уравнением Купера-Саймондса[4] и Муни-Ривлена [3, 4].
Структура модели Ми-Грюнайзена такова, что правая часть, уравнения состояния состоит из двух слагаемых, первое из них соот-
ветствует потенциальному давлению и зависит только от изменения объёма, второе обусловлено динамической составляющей процесса -колебаниями кристаллической решётки и пропорционально энергии этих колебаний е:
Р = Р П (V) + Г [e - e0(V)],
(1)
где Р - давление, РП(У) - потенциальное давление, е - внутренняя энергия, в0(У) - внутренняя энергия при начальном объёме, V - объём, Г -коэффициент Грюнайзена, характеризующий тепловое давление со стороны колеблющихся атомов.
Г (V) = V
дР де
(2)
V ^ /к
При низких температурах вклад второго слагаемого в полное давление мал, и уравнение состояния в основном определяется первым членом - потенциальным давлением. С ростом температуры роль теплового давления возрастает. Оно может сравняться с потенциальным давлением и даже превзойти его. В этих условиях значением потенциального давления пренебрегают, и вместо него в расчёт принимается зависимость теплового давления от объёма.
Если учитывать, что деформирование обтюрирующих элементов при функционировании ГИУ носит ударно-волновой характер, тогда, для определения уравнения состояния необходимо применить уравнение ударной адиабаты Гюгонио с определением коэффициента Грюнайзена [1]:
из = С + , (3)
где и5- скорость ударного фронта; и1 - массовая скорость ударно-сжатого вещества; С -скорость звука в веществе перед фронтом; 5" -показатель адиабаты:
е-^ (4)
Здесь Г - коэффициент Грюнайзена, характеризующий тепловое давление со стороны колеблющихся атомов.
Адиабатическое движение деформируемых упругопластических сред описывается следующей системой уравнений:
ё р Л ёи „ ёЕ
f + = 0; р-1
dt
=v<»; р -л
= 5 е
V ij
;е, = 2 (VU +VU); Р = P(V,E);
сЮ.
щ
а
+ 20ГЮО1] ■ 20(еп -е8п) =...
-■ ^+Ср8«);
3
а ■-р8 + Ю , X' а ер
1] г 1] ацз 2у2 1} 1
(5)
7 ■ (а0 + Бе")
(6)
ж ■ с (/1 - 3)+с (12 - 3)+с у-I)2. (6)
Она описывает как линейную, так и объёмную деформацию; С1 и С2 константы материала, связанные с модулем сдвига, обеспечивающие минимальное отклонение кривых истинных напряжений от кривых напряжений полученных экспериментально.
где р - плотностьвещества; V - дифференциальный оператор Гамильтона; V - компонент вектора скорости; а^ - компонент тензора
напряжений; е - компонент тензора скоростей деформаций; Е - энергия единицы массы вещества; Ю - компонент девиатора напряжений;
G - модуль сдвига; X - некоторый множитель, учитывающий динамическую составляющую сжимаемого вещества; а^ -метрические коэффициенты основного базиса системы координат; р - давление; V - объём; У - динамический предел текучести; - функция Кронеке-ра.
Система уравнений (5) включает дифференциальные уравнения законов сохранения массы импульса, энерги, а также физические и кинематические уравнения, связывающее внешние нагрузки с компонентами напряжений и деформаций [2].
Уравнение Купера-Саймондса дополняют систему (5) путем учетареологического упрочненияпутем учета зависимости предела текучести от скорости деформации:
Для несжимаемых материалов
дЖ
Ж
■ с
дЖ д/„
■ С2; /1, / 2 - главные инварианты деви-
атора деформаций; /1 , /2 - главные инварианты деформации Коши-Грина, J - детерминант матрицы градиента деформации,
. 3(1 - 2ц)
а ■- - коэффициент несжимаемости
20(1 + ц)
материала, ц ■ 0,49 - коэффициент Пуассона для эластомеров, G - модуль сдвига.
Инварианты деформаци Коши-Грина:
/ ■ и2 + и;+и2;
/ 2 ■ ии2 + ии+ии2;
/3 ■ и2и22и32 ■ 1 ,
(7)
Здесь: а0 ■ 92 МПа - статический предел текучести; е - мгновеннаяскорость деформации (в начальный момент времени е = 1);Д Е,р - константы материала: В=40,08 МПа - модуль упрочнения; ^=1, р=0,001 - коэффициенты Купера-Саймондса.
В модели Муни-Ривлина учитываются свойства эластомеров испытывать относительно большие конечные деформации за счёт перераспределения объёма, т.е. эти материалы считаются практически несжимаемыми с коэффициентом Пуассона близким к 0,5 [4,5]. Основной характеристикой модели Муни-Ривлена определяется функция плотности энергии деформации. Общая и широко используемая форма функция плотности потенциальной энергии деформации Ж была предложена Рив-лином [6]:
'3 ^2^3
где и1, и2, и3 - главные значения тензора деформаций. При /3 ■ 1 материал считается несжимаемым.
Напряжения в материале:
а] ■-р8] + 2—/] - 2—1 - . (8)
] ] д/х д/2
Особенностью модели Муни-Ривлена является возможность определения реальных свойств эластомеров с любой степенью точности путем сопоставления расчетных и экспериментальных значений компонент НДС при подборе коэффициентов модели. Это является значимым основанием для ее широкого применения при расчете изделий, изготавливаемых из эластомеров.
Далее приведена сравнительная оценка НДС обтюрирующего устройства, выполненного из полимерного материала с использованием моделей Ми-Грюнайзена и Муни-Ривлина. Моделирование проводилось в программном пакете ANSYS с использованием решателя Autodyn в явных Лагранжевых координатах. Задача решалась в осесимметричной постановке. На рисунке 1 представлен общий вид расчетной модели.
На рисунке 2 представлена схема приложения граничных условий.Граничными условиями решения задачи в обоих случаях определены следующие:
- перемещения изделия, заданы соответствующими положению изделия в различные моменты времени в соответствии с кривой (рисунок 4);
- давление газов р на запоясковую часть изделия, переменное во времени в соответствии с кривой изображённой на рисунке 3.
Обтюрирующие ycmpoücmöo
Рисунок 1 - Общий вид геометрической модели
Рисунок 2 - Схема приложения граничных условий
В качестве физических моделей с, учётом особенностей моделей материалов программного пакета Ansys, определены:
- для трубы ГИУ - идеальной упругости.
- для изделия - Ильюшина А.А. (малой упругого-пластической деформации с линейным упрочнением);
- для обтюрирующего устройства (полиамид) - двух параметрическая модель Муни-Ривлина и уравнение состояния Ми-Грюнайзена (таблица 1).
Определение коэффициентов модели Муни-Ривлена производилась на основаниии
Таблица 1 - Механические свойства полиамида
экспериментальной зависимости с-в (рисунок
5).
500 400 300 200 100
(D
S
и
(D
ч «
0
0
0,005
Время, с Рисунок 3 - Кривая давления
6000 S 5000
g 4000 я
0
0,005
Время, с
Рисунок 4 - Перемещение изделия в различные моменты времени
0 0,5 1
Относительная деформация
Рисунок 5 - Зависимость напряжений от деформации с линейным упрочнением для полиамида
1,5
Модель Ro, кгс/м3 et, МПа G, МПа Kr С, м/с S Q, МПа C2, МПа D, 1/МПа
Ми-Грюнайзена 1100 0,49 92 533,6 0,87 2620 1,63 - - -
Муни-Ривлина 1100 0,49 92 533,6 - - - -37,176 147 0,00194
*Здесь: Ко - плотность материала, ц - коэффициент Пуассона, О-модуль сдвига, КГ - коэффициент Грюнайзена, С- скорость звука, 5-показатель адиабаты, С1 и С2 - постоянные материала, Б - параметр несжимаемости.
0
Результаты расчета проиллюстрированы картинами напряжений и деформаций, соответствующих времени максимального давления для моделей Ми-Грюнайзена и Муни-Ривлина (см. рисунки 6 - 7).
б
Рисунок 6 - Интенсивность напряжений обтюрирующего устройства при максимальном давлении (МПа):а - модель Ми-Грюнайзена; б -модель Муни-Ривлина
Ривлена при максимальном значении интенсивности напряжений 374 МПа максимальная деформация составляет порядка 74 %., в то время как для модели Ми-Грюнайзена при максимальном значении интенсивности напряжений 239 МПа максимальная деформация составляет порядка 25 %. Для того, чтобы выяснить причину расхождения расчётов был проведён ряд численных экспериментов по сжатию цилиндрических образцов полиамидного материала с разной скоростью деформации (рисунок 8) с использованием модели Муни-Ривлина и Ми-Грюнайзена.
Моделирование, проводилось в программном пакете ANSYS с использованием решателя Autodyn в осесимметричной постановке. Механические свойства материала соответствовали таблице 1.
Граничные условия деформирования определялись получением скоростей деформации, сопоставимых со скоростями деформации обтюрирующих устройств, при работе ГИУ.
Рисунок 7- Интенсивность деформаций обтюрирующего устройства при максимальном давлении: а - модель Ми-Грюнайзена; б - модель Муни-Ривлина
Из результатов расчётов видно, что при прочих равных условиях, для модели Муни-
Рисунок 8 - Схема приложения граничных условий образца
На рисунке 9 представлены сравнительные зависимости интенсивности напряжений от скорости деформации, полученные в результате сжатия образцов.
Из результатовследует, что напряжения, полученные с использованием модели Ми-Грюнайзена и Муни-Ривлина различаются. Причем очевидным является влияние динамики нагружения, выраженное постепенным увеличением напряжений по модели Ми-Грюнайзена при росте скорости деформации на участке до 2000 с-1 вследствие учитываемых реологических свойств материала. Далее, при увеличении скорости деформирования,напряжения, полученные с использованием модели Ми-Грюнайзена резко увеличиваются, что, очевидно, объясняется влиянием ударно-волновых процессов, также учитываемыхв модели Ми-Грюнайзена.Этим и следует объяснить существенную разницу в результатах значений НДС
обтюрирующих устройств, приведённых на рисунках 6 и 7.
Таким образом, очевидно, что при оценке напряжённо деформированного состояния обтюрирующих устройств, работающих в условиях высоких скоростей нагружения, необ-
ходимо учитывать влияние динамики процесса деформирования. Критическим значением для полиамидных устройств следует принять значение скорости деформации 2000 с-1, определяемое величиной форсирования и скоростью прохождении изделия по трубе.
Рисунок 9 - Диаграмма зависимости интенсивности напряжений (МПа) от скорости деформации
Выводы
Проведен сравнительный анализ функционирования обтюрирующих устройств ГИУ, изготовленных из полиамида, с применением для расчета НДС моделей Ми-Грюнайзена и Муни-Ривлена. В результате получены значительные расхождения результатов расчетом с использованием данных моделей. Проведенные дополнительные численные исследования показали значительное влияние на НДС динамики нагружения, учитываемой уравнением состояния Ми-Грюнайзена с учетом вязкопластичного поведения материала в соответствии с уравнением Купера-Саймондса, имеющего различный характер при скоростях деформирования до и более 2000 с-1). На основании того, что реальные условия деформирования обтюрирующих устройств происходят в условиях, когда скорости деформации могут достигать значенийв несколько десятков и сотен с-1 расчеты для них целесообразно проводить с моделью Ми-Грюнайзена с учетом влияния скорости деформации на пластическое течение материала, например по уравнению Купера-Саймондса.
Литература
1. Никишов С.Ю., Лепеш Г.В. Сокращение объема ресурсных испытаний резинометаллических амортизаторов по результатам предварительных численных исследований./ Технико-технологические проблемы сервиса. №4(26), 2013 г. С.7- 10
2.Физика взрыва / Под ред. Л. П. Орленко. — Изд. 3-е, испр. - В 2 т. Т. 2. — М.: ФИЗМАТЛИТ, 2004. -656 с.
3. Уравнения состояния твердых тел при высоких давлениях и температурах. ЖарковВ. Н., КалининВ. А. Наука», 1968.
4. Dietenberger, M. Development of a high strain-rate
dependent vehicle model / M. Dietenberger, M.Buyuk, C-D Kan. - LS-DYNA Anwenderforum,
Bamberg, 2005, B-III-1-10.
5. Черных К. Ф. Нелинейная теория упругости в машиностроительных расчётах. -Л.: Машиностроение, Ленингр. отд-ние 1986.-336 с.
6 Пальмов В. А. Определяющие уравнения термоупругих, термовязких и термопластических материалов: учеб. пособие / В. А. Пальмов. — СПб.:Изд-во Политехн. ун-та, 2009. — 138 с.
