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

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

Терауд Валентин Викторович

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕХНОЛОГИЧЕСКИХ ПРОЦЕССОВ, ПРОТЕКАЮЩИХ В УСЛОВИЯХ ПОЛЗУЧЕСТИ

05.13.18 — Математическое моделирование, численные методы и комплексы программ

(технические науки)

Автореферат

диссертации на соискание ученой степени кандидата технических наук

1 9 МАЙ 2011

Москва 2011

4847650

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

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

доктор физ.-мат. наук, профессор ЛОКОЩЕНКО Александр Михайлович

Официальные оппоненты:

доктор физ.-мат. наук, профессор РЕВИЗНИКОВ Дмитрий Леонидович доктор технических наук, профессор РОМАНОВ Константин Игоревич

Ведущая организация:

Санкт - Петербургский государственный университет, математике - механический факультет

Защита состоится «16» июня 2011 г. в 14 часов, в аудитории 1804 на заседании диссертационного совета Д212.129.03 при ГОУ Московском государственном индустриальном университете по адресу 115280, г. Москва, ул. Автозаводская, д. 16, ауд.1804.

к(

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

Автореферат разослан «р»мая 2011 г.

Ученый секретарь диссертационного совета к.т.н., доцент

Кузнецов А.В.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

Для достижения поставленных целей необходимо решить следующие задачи:

1) формулировка новой математической модели, характеризующей высокотемпературные свойства материала, используемого в процессах осадки и штамповки;

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

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

4) проведение серии натурных экспериментов по осадке цилиндров для сравнения с результатами вычислений по различным моделям;

5) проведение сравнительного исследования новой математической модели с известными моделями и проведениыми экспериментами на основе задач осадки сплошных и полых цилиндров;

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

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

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

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

Научная новизна fliiccqrrainiOHHOH работы состоит в следующем:

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

2. Создан программный комплекс, предназначенный для моделирования осадки цилиндров и деформации мембран.

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

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

На защиту выносятся:

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

2. Математические модели осадки сплошных и полых круговых цилиндров.

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

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

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

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

осаживания цилиндров. Создана непротиворечивая модель деформирования мембраны под действием поперечного давления при различных ограничениях на основе предложенной математической модели. Разработанный комплекс компьютерных программ позволяет по заданным исходным данным моделировать осадку круговых цилиндров и деформирование мембран. Разработанный программный комплекс т позволяет в процессе натурных экспериментов получать деформационные характеристики цилиндра в различные моменты времени по фотографиям осаживаемого цилиндра внутри высокотемпературной печи. Полученные результаты работы внедрены в ФГУП НПЦ газотурбостроения «Салют» и ФГУП ЦИАМ им. П.И.Барапова.

Апробация работы. Основные результаты работы докладывались на 21 конференции, в том числе:

• на конференциях молодых ученых НИИ механики МГУ им М.В.Ломоносова (ежегодно в течение 2006 г. - 2010 г.);

• на Всероссийской конференции молодых ученых (с международным участием) «Неравновесные процессы в сплошных средах» (ИМСС УрО РАН, Пермь, 2007 г.);

• на научно-теоретической конференции молодых ученых, специалистов, аспирантов и студентов "Будущее авиационной науки" (ЦАГИ, г. Жуковский, 2008 г.);

• на Всероссийской конференции "Успехи механики сплошных сред" (ПАПУ ДВО РАН, Владивосток, 2009 г.);

• на 37-ой конференции по механике сплошной среды (8о1МесЬ2010), (Варшава, Польша, 2010 г.);

• на конференции "Актуальные проблемы механики сплошной среды", (Дилижан, Армения, 2010 г);

• на кафедре общей и прикладной математики МГИУ (2011 г.).

Впелренне результатов работы. Результаты работы использованы в ФГУП «НПЦ газотурбостроения «Салют», что подтверждается актом внедрения результатов диссертационной работы. Результаты диссертационной работы используются при расчёте силовых параметров ковки - штамповки крупных поковок на гидравлических прессах, скорости деформации которых составляют 10"2 - 10"3 с"1.

Результаты работы также использованы в ФГУП ЦИАМ им. П.И.Баранова при моделировании сопротивления ползучести и длительной прочности материала при одноосном и сложном напряженных состояниях, что подтверждено актом внедрения.

Публикации. По теме диссертации опубликованы 23 работы, в том числе 4 статьи в журналах, рекомендованных ВАК РФ, 8 статей в других изданиях и тезисы докладов на 11 конференциях.

Структура и объем работы Диссертационная работа состоит из введения, 4-х глав, заключения и списка литературы из 152 наименований и 1-го приложения. Работа изложена на 210 страницах машинописного текста и содержит 106 рисунков и 25 таблиц.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

способов моделирования деформирования твердых тел и методов решения полученных соотношений и систем уравнений. Особо выделяются методы, основанные на конечно-элементном подходе. Приводится обзор известных работ по рассмотренным технологическим задачам - осадке цилиндров и деформировании мембран. В последние годы важные результаты получены С.Е.Александровым, A.JI. Воронцовым, И.В. Добровым, A.M. Локощенко, H.H. Малининым, К.И. Романовым, О.В. Сосниным, СП. Яковлевым, С.С. Яковлевым, К.И. Романовым, А. IГадай, F. К. G. Odqvist, и другими учеными. В данном направлении известны классические результаты многих исследователей. Приведено описание применяемого конечно-элементного пакета LS-DYNA.

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

- степенная связь <ти, р1{ и t (модель I)

Следует отметить, что модель (3) введена в данной работе впервые. Она, с одной стороны, описывает неустановившуюся стадию ползучести материала и, с другой стороны, учитывает уровень интенсивности действующих напряжений <г„ относительно предела кратковременной прочности материала сгь при температуре испытаний. Модель (3) с помощью дробно-степенной функции от ег„ позволяет повысить уровень нелинейности зависимости />„(ст„, 0 по сравнению со стандартной степенной зависимостью (2).

В настоящей работе для выполнения компьютерного моделирования осадки цилиндров используется программный комплекс ЬБ-ОУМА. В настоящее время 15'-Л5 является эффективной программой для решения задач соударения, взрыва, разрушения, обработки металлов давлением и других. В 1£-ВКУЛ реализованы эффективные способы решения перечисленных задач. Приводится описание явного и неявного методов решения, обосновывается актуальность применения неявного метода при моделирования задач ползучести. Описана применяемая процедура решения по р-схеме Ныомарка.

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

Рассматривается задача об осесимметричной осадке кругового цилиндра высотой 2Я0 и радиусом Л0 между двумя абсолютно жесткими плитами с кинематической и

(1)

(2)

(3)

силовой программами нагружения. В первом случае нагружение осуществляется сближением плит с не зависящей от времени / скоростью 2и'(г) = 2м>0, во втором — посредством постоянной сжимающей силы Р(/) = Р„, приложенной к плитам.

Моделирование производится двумя методами: А - приближенное аналитическое решение и В - трехмерное численное решение, основанное на МКЭ.

Постановка задачи А (без учета возможного учета бочкообразования) строится для моделей материала (2) и (3), решение основывается на общих уравнениях деформируемого твердого тела и математических преобразованиях. Для двух программ нагружения получены соотношения высоты Н(1), скорости осаживания >с(0, силы осаживания Р(1) и энергии деформирования У(1).

Постановка задачи В (с учетом возможности бочкообразования) моделируется на основе конечно-элементного подхода с помощью пакета программ Ьв-ОУЭДА. При учете условий симметрии рассматривается 1/8 часть цилиндра (см. рис. 2.1, в дальнейшем называемая просто цилиндром). Уравнения равновесия внутри цилиндра удовлетворялись в ослабленной формулировке в соответствии с принципом равенства виртуальных работ внутренних и внешних сил. Дискретизация по пространственным координатам х1={х,у, г) рис ""ЩЦ

осуществлялась посредством 8-узловых или 20-

узловых гексаэдрических элементов. Суммарное количество конечных элементов цилиндра и давящей плиты составляет около 3500 штук. Для уменьшения вычислительной погрешности счета при обработке контакта, размеры конечных элементов цилиндра и контактирующей плиты имели один порядок и составляли 10 элементов по высоте, 15 по радиусу и 30 элементов по стороне плиты.

Дискретизация по времени выполнялась по неявной схеме. Модель материала (1) имеется в базовом пакете Г^-Г)У!\7Л. Для расчета с использованием остальных моделей ((2) и (3)) были написаны программы на языке БотЦап, которые интегрировались в конечно-элементный комплекс Ь8-ПУ1ЧА средствами пользовательского программирования.

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

= К )+2е(ДвГ " Л4") (' = 1- 3),

^•=^- + «346^ ».у = 1-3),

Л 'Я

шш

После этого компоненты девиатора тензора напряжений ст* выбранной модели:

масштабировались по

^''(яГ'-^'ГГ

для модели (1), для модели (2),

"1+ М

, и = Л™модели (3),

Вычисление проводилось на персональном компьютере Intel Core 2 !.86GHz с 3ГБ ОЗУ, и в зависимости от конкретной задачи оно занимало ох 5 до 80 мин для каждого комплекта исходных данных. При этих вычислениях использовалось распараллеливание счета на несколько процессоров при помощи SMP, позволившее в среднем в 1.5 раза сократить время счета.

Проведено моделирование осадки цилиндра при различных значениях геометрических параметров: i?o = 19-5 мм и 39 мм, Но — 9.5 мм, 19.5 мм и 39 мм. Цилиндр осаживался до уровня (Но - H\)IHq = 0.4 в течение /, = 60, 740 и 2000 сек по двум программам нагружения: 1) постоянная скорость штампа и'(/) = w0; 2) постоянная сжимающая сила P(t) = Р0 - 57 - 75 кН, в зависимости от модели и коэффициента трения. Контактные условия между штампом и цилиндром моделировались законом Кулона с задаванием полного и неполного проскальзывания (ц = 0 и 0.5), а также при полном контактном прилипании ( ц=ц *).

Пример получаемых деформационных сеток представлен на рисунке 2.2. где изображена RZ-проекция конечно-элементной модели при различных значениях |л ((б) - ц — 0 и (в) - (1 = 0.5) и ((а) -1 = 0 и (б), (в) -1 = 740 сек) для (3) модели материала. Следует отметить, что применение различных моделей (1) — (3) не приводит к существенному изменению геометрии конечной сетки. При наличии полного контактного прилипания наблюдается интенсивный переход материала цилиндра с боковой поверхности на контактную, которая увеличивается только за счет перетекания материала.

На основе данной задачи в работе проведено сравнение трех рассмотренных моделей материала. Проводилось сравнение высот цилиндра 2H(t), бочкообразное™ K(f) = (Лщзх -ЛшпУЛшах. силы осаживания P(t), энергии деформирования V(f) и других, в зависимости от времени /, коэффициента трения //, программы нагружения и других параметров. На рис. 2.3 приведены зависимости H(t) при различных моделях материала (1) - (3) (кривые 1 - 3 соответственно), осаживаемого по силовой программе P(t)=P0 при отсутствии трения (2R0= 39 мм). Модели материала (!) и (2) дают схожие кривые (немного различающиеся в середине процесса осадки и совпадающие в начале и конце процесса). Модель материала (3) показывает несколько большую осадку за то же время, но, как показало экспериментальное исследование (см. главу 3), более близкую к данным, полученным в натурных экспериментах, по сравнению с моделями (1) - (2).

о

200 401) 600 800

I, сек Рос. 2.3.

О 100 21» 300 400 500 600 700

и сек Рис. 2.4.

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

Моделирование полого цилиндра по силовой и кинематической программам деформирования проведено при разных коэффициентах трения при использование различных моделей материала. Было получено, чго при увеличении коэффициента трения сила осаживания уменьшается до 20% при использован™ модели (1), в то же время для (3) модели такая зависимость не наблюдается.

На рис. 2.4 изображены зависимости высоты цилиндра 2Н{() при отсутствии трения для трех моделей, из которых следует получаемый различный характер осаживания цилиндра. Модели (2) - (2) характеризуются быстрым укорочением цилиндра в кратковременный начальный период осаживания и гораздо более медленным в оставшееся время. В свою очередь модель (3) не проявляет резкого изменения скорости сШ/<Ь (по сравнению с моделями (1) - (2)), при этом осаживание производится в условиях, приближенных к постоянной скорости деформирования. Следует отметить, что проведенные эксперименты показывают, что резкие изменения скорости осаживания, наблюдающиеся при использовании моделей (1) и (2), не происходят. Модель материала (3) наиболее правильно описывает экспериментальные данные (см. главу 3).

При моделировании вычислялась энергия деформирования Г(1), соответствующая различным программам осаживания. Данные вычисления показали, что энергия деформирования V.,., затрачиваемая на осадку цилиндра по кинематической программе, меньше, чем энергия деформирования К,,, накопленная по силовой программе. Разность -Уи.) возрастает с увеличением коэффициента трения ц и отношения Д,/#0 независимо от типа модели. В рассматриваемом примере эта разница достигала 20 %.

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

Оценка адекватности численного моделирования производилась путем его сравнения с точным аналитическим решением при отсутствии трения. Данное сравнение показало полное совпадение зависимостей Щг) , Р{{) , К(/), и'(/) и У(г) для обеих программ погружения и грех рассмотренных моделей. Кроме того, зависимости м>(/) и Н(1) для силовой программы нагружеиия сравнивались с соответствующими экспериментальными данными. Сходимость численной процедуры счета обусловлена в основном сходимостью применяемого метода Ныотона-Рафсона при решении системы нелинейных

дифференциальных уравнений МКЭ. Исследования показали, что при использовании модели (3) скорость сходимости вырастает в несколько раз.

Оценка гогрешности численного моделирования производилась с помощью метода Рунге-Ромберга при использовании моделей (2) и (3) и уменьшении шага по времени и уменьшении характерного размера элемента в 2 раза. Новая модель (3) имеет меньшую погрешность при вычислениях.

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

В третьей главе разработан метод бесконтактных измерений, позволяющий исключить проблемы измерений при высокотемпературных испытаниях путем замены реального нагретого образца на аналогичную ему компьютерную модель (числовой образ). Разработаны методологическая, физическая и программная базы. Метод был апробирован при проведении натурных экспериментов по осадке алюминиевых цилиндров различной геометрии при различных параметрах, позволивших подтвердить достоверность проведенного моделирования осаживания цилиндров (см главу 2). Для обработки полученной графической информации разработан пакет программ АзЬгп. Для определения материальных констант разработаны метод и компьютерная программа для проведения одноосных экспериментов на ползучесть при квазипостоянных напряжениях, обеспечиваемых кусочнопостояными растягивающими силами. По расчетам компьютерной программы проведены натурные эксперименты на одноосную ползучесть, из которых определялись константы материала в моделях (1) - (3). Проведен сравнительный анализ при использовании моделей (1) - (3) с данными, полученными в натурных экспериментах. Применяются метод твердотельного моделирования на основе пакета 8оН(1'\Уогк5, метод наименьших квадратов и ряд математических методов расчета.

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

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

распознавания границ и реперных точек образца. Вычислитель использует распознанные данные фотографий редактором для вычисления числовой метрической геометрии и параметров осаживания цилиндра. Из полученных метрических данных при помощи постобработчика могут быть созданы графические зависимости вычисленных величин. Постобработчик может визуализировать двумерные графические данные в виде полей на поверхности цилиндра и в виде видео процесса осадки. На рис. 3.2 показано главное окно программы автоматического редактора. Цифрами обозначено 1 - основное изображение цилиндра (увеличенное), 2 - меню параметров настройки. 3 - графики цветности и яркости выбранного горизонтального сечения, 4 - график выбранного вертикального сечения (аналогично 3), 5 - миниатюра всего кадра.

Рис. 3.2.

Данный комплекс последовательно выполняет следующие действия: распознает цилиндр на кадре, выделяет его границы, распознает реперные линии и точки их пересечения, вычисляет "пиксельные" координаты границ и реперов, преобразует "пиксельную" геометрию в метрическую и вычисляет деформационные характеристики. Границы цилиндра на боковых сторонах аппроксимировались полиномами VII степени, на торцевых сторонах прямыми линиями. По изменениям расстояний между реперными линиями вычисляются компоненты вектора перемещений и затем тензора деформаций на внешней поверхности деформированного цилиндра. Всё распознавание происходит в автоматическом режиме, с поддержкой пакетной обработки множества фотографий и ручной правки. Разработанная система измерений позволяет вычислять 15 параметров: высоту и диаметр цилиндра в любом поперечном сечении Н{г, /), £>(г, г), координаты и скорости перемещений произвольных точек на поверхности цилиндра Х{х.у\1) и У(х, V,;). бочкообразности К((), деформации на видимой поверхности цилиндра р(х,у,1) и другие параметры. Комплекс АвШт позволяет получать видеоизображение процесса осадки.

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

Разработанный метод бесконтактных измерений и пакет программ Аькт апробировался на экспериментальном исследовании осадки цилиндров из алюминиевого сплава нри высокой температуре. Были использованы цилиндры с различными геометрическими размерами (2Я0, 2Й0): (39, 39), (78, 39), (19.5, 39) и (39, 78), подготовленные специальным образом. Были проведены 23 испытания по силовой программе нагружения с различными значениями сжимающей силы. Измерение геометрических параметров цилиндра производилось разработанным оптическим методом, полученные данные обрабатывались пакетом АякШ.

На основе проведенного экспериментального исследования осадки цилиндров были сделаны уточнения в конечно-элементной модели. Использование вычисленных материальных констант из испытаний на одноосную ползучесть позволило провести моделирование осадки круговых цилиндров с использованием моделей (1) - (3) для алюминиевого сплава. На основе результатов экспериментов проведено верифицирование результатов моделирования по зависимостям Н(1), р:(1), р0(:, /), р.(г. /).

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

Была проведена серия из 39 натурных экспериментов на образцах из алюминиевого сплава при высокой температуре, на основе разработанного метода одноосного квазипостоянного нагружения, который позволил получить относительную погрешность растягивающего напряжения Аег/сг3 не более 0.3%. Из этих экспериментов определялись материальные константы в моделях (1) - (3): А, и, I, В, т, у, С, аь, /. на основе метода наименьших квадратов.

На рис. 3.3 и 3.4 в качестве примера приведено сравнение двух экспериментов (кривые 4, 5) и результатов, полученных численным моделированием при использовании моделей (1) - (3) (кривые 1, 2, 3). При рассмотрении зависимости Н(1) во времени все три модели хорошо согласуются с экспериментальными данными до уровня осадки порядка 20%. При дальнейшем увеличении осадки наилучшее согласование экспериментальных и теоретических результатов осуществляется при использовании модели (3).

Сравнение экспериментальной и теоретических зависимостей окружной деформации рд(!) в центральной части цилиндра (рис. 3.4) показывает, что в первой половине процесса

осадки следует отдать предпочтение модели (3), а при последующей осадке - модели (2).

Для оценки погрешности разработанного оптического метода измерений были проведены несколько экспериментов, в процессе которых параллельно с оптической системой измерения производилось прямое измерение высоты образца и максимального диаметра в процессе осадки. Результаты сравнения показали, что предложенный оптический метод позволяет получить значения геометрических параметров цилиндра с высокой точностью. При исключении высокочастотного шума из данных можно получить характерную попэешность порядка 0.3 мм. После испытаний геометрия каждого образца измерялась микрометром и сравнивалась с конечной геометрией, полученной оптическим методом. Для примера в табл. 3.1 приведены значения распознанных высот цилиндра 2Н2 и измеренной после испытания 2Н\. Было получено, что среднее типичное отклонение составило величину около 0.16 мм. Такое же измерение производилось для диаметра в средней части образца.

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

Рассмотрим деформирование длинной узкой прямоугольной мембраны (рис. 4.1) шириной 2а и начальной толщиной Н0, закрепленной вдоль длинных сторон и

нагруженной равномерным поперечным давлением q.

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

На первой стадии мембрана, плоская в начальном состоянии, под действием давления д мгновенно упруго деформируется, приобретая форму незамкнутой круговой цилиндрической оболочки с центральным углом 2а, (рис. 4.!). Решение строится на основе закона Гука. Проводится численный эксперимент, показывается

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

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

2Н, [мм] ] 2Н2 [мм]

2Н„=19.5 мм, К0=19.5 мм

12,0 11,8

11,6 11,6

11,0 10,7

10,7 10.8

2Но=39 мм, Я0=19.5 мм

22,0 22,2

22,3 22,2

24,5 24,5

24,3 24,3

2Н0=39 мм, Й^ЗЗ ММ

26,5 27.6

24.3 24,3

Табл. 3.1.

первой стадии. На основе предложенной модели (3) выводятся основные соотношения, зависящие от угла раствора мембраны. Зависимость угла раствора мембраны а(Г) от безразмерного времени 1 принимает вид:

-1

1

da- а> <sí°2-

2 sin*(а)

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

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

ав = 2/л/з , ав=агв1аь.

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

На третьей стадии ползучесть мембраны происходит внутри жесткой матрицы при наличии трения о ее стенки (закон Кулона). Задача решается итерационным методом. Конечноразностные соотношения принимают вид:

= dsf + 4ds™

КГ

H„H

dst

H™=H>

2 Я0Я;

(

(dt^)í

aPt

, a ctga

1-a-ctgff

1 - a ■ ctge

HJ1..

2 H.Hl

—з---WiW IM +-ж--

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

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

II,

U

0.8

40

...

f i

.............. - /' ^

1

80 120 Ж Ж Ь, ММ

Рис. 4.3

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

В качестве примера на рис. 4.3 приведены полученные в результате моделирования величины толщины мембраны Н и ее ширины Ъ конце каждой стадии (стадии 1, 2, 3) деформирования. Начальная толщина мембраны Но - 2 мм, начальная ширина 2а = 100 мм. Отрезок прямой 2 показывает, что мембрана в конце второй стадии имеет ширину 117 мм и толщину 1.7 мм, что обусловлено ее поперечным растяжением. В конце третьей стадии (кривая 3) толщина мембраны не является постоянной величиной по ее ширине, что обусловлено наличием трения о стенки матрицы.

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ И ВЫВОДЫ

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

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

3. Предложен метод определения минимальной энергии деформирования при осадке цилиндров.

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

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

ОСНОВНЫЕ ПОЛОЖЕНИЯ ДИССЕРТАЦИОННОЙ РАБОТЫ ОПУБЛИКОВАНЫ В СЛЕДУЮЩИХ НАУЧНЫХ ТРУДАХ

Статьи в журналах, рекомендованных ВАК РФ:

1. Локощенко A.M., Демин В.А., Носов В.В. (Терауд В.В.) Осадка кругового цилиндра в условиях установившейся ползучести. // Известия ВУЗов. Машиностроение. 2007. №4. С. 3-10.

2. Локощенко A.M., Моссаковский П.А., Терауд. В.В. Исследование осадки круговых цилиндров при ползучести с учетом и без учета бочкообразования. // Вычислительная механика сплошных сред. 2010. Т.З. №1. С. 52-62.

3. Локощенко A.M., Терауд. В.В. Экспериментальное подтверждение моделирования осадки цилиндров при ползучести. // Машиностроение и инженерное образование. 2011.№1. С. 24-28.

4. Терауд В.В.. Моделирование высокотемпературной осадки цилиндров при различных моделях материала. // Известия Самарского гос. тех. ун-та. Серия: Физико-математические науки. 2011. №1. С. 283-290.

Статьи в других изданиях:

5. Терауд В.В. Экспериментально-теоретическое исследование высокотемпературной осадки цилиндров // Труды конференции-конкурса молодых ученых НИИ механики МГУ им. М.В. Ломоносова (8-10.10.2008). Изд-во МГУ. М. 2009. С. 215-220.

6. Терауд В.В. Инновационная многокомпонентная система бесконтактных измерений при высоко температурных испытаниях // Труды конференции-конкурса молодых ученых НИИ механики МГУ им. М.В. Ломоносова (14-16.10.2009). Изд-во МГУ. М. 2010. С. 307-317.

7. Lokoshchenko A.M., Teraud V.V. Experimental and theoretical study of the high temperature creep shortening for a circular cylinder // 37й1 Solid mechanics conference (6.09-10.09.2010, Warsaw, Poland). Inst, of Fundam. Technol, Research. Wars. 2010. P. 194-195.

Тезисы докладов на конференциях:

8. Локощенко A.M., Терауд В.В. Деформирование мембраны под действием поперечного давления (учет мгновенного, свободного и стесненного деформирования) // Научная конф. "Ломоносовские чтения" (апрель 2009). Секция механики. Тезисы докладов. Изд-во МГУ. М. 2009. С 106.

9. Локощенко A.M., Моссаковский П.А., Терауд В.В. Двумерное решение осадки сплошных цилиндров при ползучести. // Успехи механики сплошных сред (29.0505.10. 2009, Владивосток). Тезисы Всероссийской конференции, приуроченной к 70-летию академика В.А. Левина. Владивосток.: Дальнаука. 2009. С. 89 - 90.

Терауд Валентин Викторович

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕХНОЛОГИЧЕСКИХ ПРОЦЕССОВ, ПРОТЕКАЮЩИХ В УСЛОВИЯХ ПОЛЗУЧЕСТИ

Автореферат

Подписано в печать 05.05.11 Формат бумаги 60x84/16 Усл. печ. л. 1,0. Уч.-изд. л. 1,0. Тираж 100. Заказ № 137

Издательство МГИУ, 115280, Москва, Автозаводская, 16 www.izdat.msiu.ru; e-mail: izdat@msiu.ru; тел. (495) 620-39-90

ОГЛАВЛЕНИЕ

ВВЕДЕНИЕ

ГЛАВА 1. Введение в ползучесть. Описание моделей. Обзор литературы по тематике. Обзор конечно-элементного пакета моделирования ЬБ-Оупа.

1.1. Введение в ползучесть.

1.2. Базовые модели ползучести, при одноосном растяжении.

1.3. Описание моделей, используемых в диссертации.

1.4. Обзор методов решения и способов моделирования.

1.5. Обзор рассмотренных технологических задач.

1.6. Описание пакета программ Ьз-ТУупа.

ГЛАВА 2. Моделирование осадки сплошных и полых цилиндров при различных подходах к решению, различных краевых условиях, различных моделях материалов. Исследование различий моделей материала. Проведение численного эксперимента. Разработка пакета программ для моделирования осадки цилиндров.

2.1. Математическая формулировка задачи.

2.2. Аналитическое решение без учета бочкообразования.

2.3. Численный подход (МКЭ).

2.4. Численное моделирование. Обсуждение результатов исследования.

2.5. Оценка погрешности и сходимости.

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

2.7. Поиск оптимальной программы нагружения.

2.8. Результаты решения прикладной задачи.

2.9. Экспериментальная верификация моделирования.

2.10. Пакет программ для моделирования осадки цилиндров.

2.11. Выводы по второй главе.

ГЛАВА 3. Разработка методологической, физической и программной баз для бесконтактного метода измерений. Разработка программного комплекса Авкпп. Разработка метода и компьютерной программы для проведения экспериментов при квазистатическом напряжении. Сравнение экспериментальных данных с проведенным моделированием.

3.1. Аппаратно - программный метод бесконтактных измерений.

3.2. Программная среда Askim.

3.2.1. Ручной редактор.

3.2.2. Автоматический редактор.

3.2.3. Вычислитель.

3.2.4. Постобработчик.

3.3. Аппаратная реализация метода.

3.4. Оценка точности и погрешности разработанного бесконтактного метода.

3.5. Эксперименты по осадке цилиндров.

3.6. Метод квазипостоянного нагружения. Программа Samlot.

3.7. Определение физических констант в моделях материала при помощи метода квазипостоянного нагружения.

3.8. Экспериментальная проверка проведенного моделирования.

3.9. Выводы по третьей главе.

ГЛАВА 4. Моделирование деформирования мембраны на основе предложенной новой модели. Проведение численного эксперимента. Разработка пакета программ для моделирования деформирования мембраны.

4.1. Введение.

4.2. Математическая формулировка задачи.

4.3. Стадия упругого деформирования.

4.4. Численный эксперимент.

4.5. Стадия свободного деформирования при ползучести.

4.6. Потеря устойчивости.

4.7. Численный эксперимент.

4.8. Стадия стесненного деформирования в клиновидной матрице при ползучести.

4.9. Потеря устойчивости и разрушение.

4.10. Численный эксперимент.

4.11. Пакет программ для моделирования деформации мембраны.

4.12. Программа упругого деформирования (upr m). 178 4.13 Программа свободного деформирования (svb m).

4.14. Программа стесненного деформирования (klin m).

4.15. Объединительный модуль (3state).

4.16. Выводы по четвертой главе. 184 ЗАКЛЮЧЕНИЕ 185 СПИСОК ЛИТЕРАТУРЫ 186 ПРИЛОЖЕНИЕ

5.1. Метод секущих.

5.2. Метод интегрирования Симпсона.

5.3. Метод интегрирования Гаусса.

5.4. Численная реализация моделей.

5.5. Программа Samlot.

5.6. Модуль чтения метаданных.

5.7. Модуль изображения графиков.

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

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

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

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

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

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

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

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

Для достижения поставленных целей необходимо решить следующие задачи:

1) формулировка новой математической модели, характеризующей высокотемпературные свойства материала, используемого в процессах осадки и штамповки;

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

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

4) проведение серии натурных экспериментов по осадке цилиндров при различных исходных значениях параметров, для сравнения с результатами вычислений по различным моделям;

5) проведение сравнительного исследования новой математической модели с известными моделями и проведенными экспериментами на основе задач осадки сплошных и полых цилиндров;

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

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

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

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

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

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

2. Создан программный комплекс, предназначенный для моделирования осадки цилиндров и деформации мембран.

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

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

На защиту выносятся:

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

2. Математические модели осадки сплошных и полых круговых цилиндров.

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

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

5. Бесконтактный метод измерений и программный комплекс АбШгп для обработки фотографической информации, полученной в натурных экспериментах.

Практическая значимость. Предложена новая математическая модель для описания неустановившейся ползучести, учитывающая уровень действующего напряжения относительно предела прочности материала. Разработана модель для описания осадки сплошных и полых цилиндров. Получено условие энергетической эффективности осаживания цилиндров. Разработана непротиворечивая модель деформирования мембраны под действием поперечного давления при различных ограничениях на основе предложенной математической модели. Разработанный комплекс компьютерных программ позволяет по заданным исходным данным моделировать осадку круговых цилиндров и деформирование мембран. Разработанный программный комплекс Азкгт позволяет в процессе натурных экспериментов по фотографиям осаживаемого цилиндра внутри высокотемпературной печи получать деформационные характеристики цилиндра в различные моменты времени. Полученные результаты работы внедрены в ФГУП НПЦ газотурбостроения «Салют» и ФГУП ЦИАМ им. П.И.Баранова.

Апробация работы. Основные результаты работы докладывались на 21 конференции, в том числе:

• на конференциях молодых ученых НИИ механики МГУ им М.В.Ломоносова (ежегодно в течение 2006 г. — 2010 г. );

• на всероссийской конференции молодых ученых (с международным участием) «Неравновесные процессы в сплошных средах» (ИМСС УрО РАН, Пермь, 2007 г.);

• на научно-теоретической конференции молодых ученых, специалистов, аспирантов и студентов "Будущее авиационной науки" (ЦАГИ, г. Жуковский, 2008 г.);

• на Всероссийской конференции "Успехи механики сплошных сред" (ИАПУ ДВО РАН, Владивосток, 2009 г.);

• на 37-ой конференции по механике сплошной среды (8о1МесЬ2010), (Варшава, Польша, 2010 г.);

• на конференции "Актуальные проблемы механики сплошной среды", (Дилижан, Армения, 2010 г);

• на кафедре общей и прикладной математики МГИУ (2011 г.).

Внедрение результатов работы. Результаты работы использованы в ФГУП «НПЦ газотурбостроения «Салют», что подтверждается актом внедрения результатов диссертационной работы. Результаты диссертационной работы используются при расчёте силовых параметров ковки - штамповки крупных поковок на гидравлических прессах, скорости деформации которых составляют 10'2 — 10"3 с"1.

Результаты также работы использованы в ФГУП ЦИАМ им. П.И.Баранова при моделировании сопротивления ползучести и длительной прочности материала при одноосном и сложном напряженных состояниях, что подтверждено актом внедрения.

Публикации. По теме диссертации опубликованы 23 работы, в том числе 4 статьи в журналах, рекомендованных ВАК РФ, 8 статей в других изданиях и тезисы докладов на 11 конференциях.

Структура и объем работы Диссертационная работа состоит из введения, 4-х глав, заключения и списка литературы из 152 наименований и 1-го приложения. Работа изложена на 210 страницах машинописного текста и содержит 106 рисунков и 25 таблиц.

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

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

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

3. Предложен метод определения минимальной энергии деформирования при осадке цилиндров.

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

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

ЗАКЛЮЧЕНИЕ.

1. Абрамов A.A., Юхно Л.Ф. Один метод исключения для линейных задач // Ж. вычисл. матем. и матем. физ. 1998. Т. 38 . №4. С. 547-556.

2. Азаренок Б.Н. О построении структурированных сеток в двумерных невыпуклых областях с помощью отображений // Журнал вычислительной математики и математической физики. 2009. Т. 49. № 5. С. 826-839.

3. Александров С.Е. О разрывных полях скоростей при произвольной деформации идеального жесткопластического тела // Докл. РАН. 1992. Т. 324. №4. С. 769-771.

4. Александров С.Е. Поле скорости вблизи поверхности их разрыва при произвольном течении идеального жесткопластического материала. Механика твердого тела. 1995. №5. С. 116-122.

5. Александров С.Е., Вилотич Д., Гольдштейн Р.В., Чиканова H.H. К определению диаграммы пластичности // Механика твердого тела. 1999. №4. С. 141-148.

6. Александров С.Е., Гольдштейн Р.В. Вариационный принцип в теории пластического течения с учетом образования трещин отрыва // Докл. РАН. 1995. Т. 342. № 4. С. 484-486.

7. Александров С.Е., Данилов B.JL, Чиканова H.H. О зоне торможения при моделировании осесимметричных процессов обработки металлов давлением в условиях ползучести // Механика твердого тела. 2000. №1. С. 149-151.

8. Александров С.Е., Друянов Б.А. О типе уравнений двумерных течений. идеальной жесткопластической среды при условии текучести Драккера. //Механика твердого тела. 1989. №6, С. 155-158.

9. Александров С.Е., Друянов Б.А. Об условиях трения пластических тел // Механика твердого тела. 1992. №4. С. 116-122.

10. Андреев A.B. Развитие методов прямого численного решения одномерных интегродифференциальных уравнений механики // Механика твёрдого тела. 2007. №2. С. 50-65.

11. Базаренко H.A. Взаимодействие полого цилиндра конечной длины и плиты с цилиндрической полостью с жесткими вкладышами // Прикладная математика и механика. 2010. Т.74. Вып.З. С. 455 468.

12. Базаренко H.A. Контактная задача для полого и сплошного цилиндров со свободными от напряжений торцами // Прикладная математика и механика. 2008. Т. 72. №2. С.328-341.14.17.18,19.20,2124