Geum.ru
   Авторефераты по всем темам  >>  Авторефераты по техническим специальностям

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

УДК 539.3 ПЕКАРЬ

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

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

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

Тверь 2012

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

Официальные оппоненты: Шешенин Сергей Владимирович, доктор физико-математических наук, профессор, профессор МГУ им. М.В. Ломоносова Шеретов Юрий Владимирович, доктор физико-математических наук профессор, профессор Тверского госуниверситета Ведущая организация ООО Научно-технический центр НИИШП

Защита состоится л29 мая 2012 года в 14:00 часов на заседании диссертационного совета Д212.263.04 при Тверском государственном университете по адресу: 170002, г. Тверь, Садовый пер., 35, ауд. 2

С диссертацией можно ознакомиться в библиотеке Тверского государственного университета по адресу: 170100, г. Тверь, ул. Володарского, 44а Объявление о защите диссертации опубликовано л28 апреля 2012 года на сайте ВАК и на официальном сайте Тверского государственного университета по адресу: Автореферат разослан л28 апреля 2012 года И.О. ученого секретаря диссертационного совета, доктор физико-математических наук Зингерман К.М.

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

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

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

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

Задачам нелинейной теории упругости посвящено большое количество как отечественных, так и зарубежных работ, в частности работы Г.М. Бартенева, М.Ф. Бухиной, И.И. Блоха, Л.М. Зубова, Ю.И. Койфмана, Л.И. Кутилина, А.И. Лурье, Н.Ф. Морозова, В.В. Новожилова, В.А. Пальмова, П.М. Риза, Л.И. Седова, Г.С. Тарасьева, Л.А. Толоконникова, Т.Н. Хазановича, К.Ф. Черныха, P.J. Blats, A.E. Green, W.L. Ko, M.A. Moony, F.D. Murnaghan, W. Noll, R.S. Rivlin, L.R.G. Treloar, C. Truesdell, O. Watanabe, W. Zerna и многих других.

Более сложными являются задачи связанные с анализом поведения элементов конструкций, когда необходимо учитывать вязкоупругие процессы, происходящие в материале тела. Теории и эксперименту задач вязкоупругости посвящены работы А.А. Адамова, В.П. Матвеенко, Н.А. Труфанова, И.Н. Шардакова, Ю.Н. Работнова и его учеников, цикл работ А.А. Ильюшина и Б.Е. Победри, работы Н.Х. Арутюняна и его учеников, а также работы Р. Кристенсена. В этих работах построены модели и дано теоретическое обоснование для малых и конечных деформаций.

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

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

Для получения результатов данной работы использовался метод конечных элементов (МКЭ). Применение метода конечных элементов к задачам линейной и нелинейной теории упругости подробно рассмотрено в работах Л. Сегерлинда, О. Зенкевича, Дж. Джу и Р. Тейлора.

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

Целями работы являются:

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

- разработка детализированного алгоритма решения квазистатических задач с учетом нелинейностей высших порядков и алгоритма решения динамических задач в рамках описанной модели;

- разработка программного модуля для решения задач;

- проведение численных экспериментов для конкретных задач.

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

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

Практическая ценность работы заключается в возможности использования результатов расчета на начальных стадиях проектирования изделий из резиноподобных и полимерных материалов, а также в задачах мониторинга (например, когда в процессе эксплуатации возникает и может начать развиваться дефект). Разработанный программный модуль использовался при выполнении работ по грантам РФФИ (проекты 11-01-12043-офи-м-2011, 11-08-01284-а), госконтракту № 8757р / 14004 от 14 января 2011 г. с Фондом содействия развитию малых форм предприятий в научно-технической сфере по программе Участник Молодёжного Научно-Инновационного Конкурса (У.М.Н.И.К.) и при разработке CAE FIDESYS (Проект поддержан фондом Cколково).

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

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

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

Ч результаты численных экспериментов для конкретных задач.

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

Результаты работы полностью или частично докладывались автором и обсуждались на следующих научных семинарах и конференциях: Научные конференции Современные проблемы математики, механики, информатики в 2005, 2006, 2008, 2011 гг. (г. Тула); Научные конференции Ломоносовские чтения в 2006, 2007, 2008, 2009, 2010. гг. в МГУ им. М.В. Ломоносова (г. Москва); 7 научная конференция Актуальные проблемы состояния и развития нефтегазового комплекса России в Российском государственном университете нефти и газа им. И. М. Губкина в 2007 г. (г. Москва); Научно - практическая конференция линженерные системы - 2009, 6 - 9 апреля 2009 года, Москва;

Восемнадцатый и Девятнадцатый симпозиум Проблемы шин и резинокордных композитов в 2007 и в 2008г. соответственно (г. Москва); Всероссийской научнопрактической конференции Математика, информатика, естествознание в экономике и в обществе в 2007 году (г. Москва); 73rd EAGE Conference & Exhibition incorporating SPE EUROPEC 2011, Vienna, Austria.

Публикации. Основные результаты диссертации представлены в публикациях, 2 из которых в изданиях из списка ВАК.

Структура и объём работы. Диссертация состоит из введения, четырех глав, заключения, приложения и списка использованных источников из 1наименований. Работа изложена на 121 странице машинописного текста, содержит 33 рисунка.

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

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

В п.1.1 приводятся основные термины и обозначения, используемые в работе:

n n R - радиус-вектор частицы в -м состоянии;

i - лагранжевы (материальные) координаты частицы;

n n э - базисные векторы в -м состоянии;

i n nun R R - вектор перемещений, характеризующий переход из n предыдущего (n 1)-го состояния в последующее -е состояние;

p - градиент в базисе p-го состояния;

pp qp - аффинор деформаций, q, p I un (I un) nq1 nqхарактеризующий переход из q-го состояния в p-е;

Gq, p q, p T, p - тензорная мера деформаций, описывающая изменение q деформаций при переходе тела из состояния q в состояние p и соответствующая мере Грина ( - тензорная мера Грина);

G0,pn n n, fn, - плотность, массовая сила и давление в -м состоянии;

m n mn - относительное изменение объема при переходе из -го в -е, состояние;

0,n - тензор истинных напряжений, накопленных в теле при переходе из n начального в -е состояние;

0,n (0,nT )1 0,n 0,n - энергетический тензор напряжений, накопленных n в теле при переходе из начального в -е состояние;

nn n 1 0,n - тензор обобщенных (полных для -го состояния) 0,n 0,n напряжений, определенный в координатном базисе n-го состояния;

m n - тензор обобщенных (полных для -го состояния) напряжений, 0,n m определенный в координатном базисе произвольного -го состояния:

mn T,m n,m ;

0,n 0,n n m m m - тензор обобщенных дополнительных напряжений, q,p 0,p 0,q m определенный в координатном базисе произвольного -го состояния;

k n - граница тела в -м состоянии в координатах k -го состояния;

n k k N - нормаль к ;

n n - знак двойной скалярной свертки;

T - знак транспонирования.

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

n В п.1.3. приведены уравнения равновесия и движения на момент времени n в координатном базисе состояния на момент времени, а также уравнения n равновесия и движения на момент времени в координатном базисе состояния k (k n) на момент времени.

В п.1.4. описаны граничные условия, используемые в работе.

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

Для несжимаемого материала I (t) (G0,n : I)G0,n1 pnG0,n1, 0,n (1) где - интегральный оператор вида:

t Aexp(t) [(t)] 0 (t) )( )d K(t) . (2) , K(t t1 0 Для сжимаемого материала 0,n(t) b((G0,n : I) 3)]n(1 0,n)5/3 I (G0,n : I)G0,n1 B0,nG0,n1, (3) [1 где - интегральный оператор вида:

t Aexp(tm) [(t)] 0 (t) )( )d K(t) . (4) , K(t t1 0 Вторая глава посвящена механической и математической постановке граничных задач теории многократного наложения больших деформаций об образовании концентраторов напряжений в предварительно нагруженных телах.

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

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

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

Рис. 1.а Рис. 1.б Рис. 1.в Отметим, что форма и взаимное расположение отверстий могут задаваться как в момент образования, так и в другой, заранее заданный, момент времени 1 ' 1). И это, естественно, обуславливает математическую постановку задачи.

( В п.2.2 рассмотрено восемь случаев полной математической постановки задач об образовании концентраторов напряжений в телах из нелинейно-вязкоупругих тел: форма концентратора может быть задана как в момент образования k, так и в момент n k ; материал может моделироваться как сжимаемым (в том числе и слабосжимаемым), так и несжимаемым материалом; постановка задачи может быть как квазистатической, так и динамической.

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

Уравнение движения:

kk 2u(t) , (5) 1 0,k (t) k,n(t) k 0,n t u(t) здесь - вектор перемещений, характеризующий переход из состояния n k на момент в состояние на момент и выраженный в координатах k-го состояния.

Граничные условия:

k k k N (t) (t) P(t)(1 0,n(t)) N (t) *,n1(t) k,n1(t). (6) n 0,n n k k n (t) Начальные условия:

u u(tk ) 0; (tk ) vk, (7) t vk k 0 где - скорости, накопленные к моменту времени, они равны нулю для.

Уравнение несжимаемости:

0,n 0. (8) Определяющие соотношения для несжимаемого вязкоупругого материала (1), (2).

Геометрические соотношения:

1 0,n(t) det0,n(t), (9) 0,n(t) 0,k k,n(t) (10), k k,n I u(t). (11) Связь между плотностями в различные моменты времени k 1 0,k 0. (12) Завершает постановку задачи связь между тензором обобщенных напряжений в различных координатных базисах:

k (t) T (t) 0,k. (13) 0,n 0,n 0,k u(t) Решение этой задачи позволяет найти, в частности, вектор k перемещений из k-го состояния в n-е как функцию радиус-вектора R, т.е. в координатах k-го состояния.

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

Третья глава посвящена методам решения задач.

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

В п.3.2 обсуждаются особенности применения метода конечных элементов для несжимаемых материалов и слабосжимаемых материалов.

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

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

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

k 1 k k U U V (1 2)Ak 2Ak 1, (14) k 1 k V V (1 1)Ak 1Ak 1. (15) Схема обладает вторым порядком аппроксимации. Безусловная устойчивость данного метода гарантирована при 2 1 0.5.

В п.3.6 описывается общий алгоритм расчета интеграла (2) и (4) из вязкоупругих определяющих соотношений.

T Пусть необходимо вычислить интеграл вида:.

I (T ) K(T )g( )d Предположим, что каким-то образом нам удалось получить приближение функции K (t) ядра из (2) или (4), линейной комбинации конечного числа экспонент на [0,T ] отрезке, то есть удалось подобрать такие Ci и i, что m i K(t) e t ,t [0,T ].

C i im it Такая комбинация Ce называется рядом Прони. Тогда вычисление i iинтеграла можно производить по следующей рекуррентной формуле m (1i ) i Ik 1 : I(tk 1) , (16) C Ik g(t*k ) i i i i tk i i (tk ) Ik g( )dЧ интегралы, вычисленные на предыдущем шаге по где e g(tk 1) g(tk ) i i e (tk1tk ), а gk 1(t*) времени, а i определяются как Преимуществом использования такого метода вычисления интеграла является то, что для использования формулы необходимо хранить в памяти лишь интегралы, вычисленные на предыдущем шаге.

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

В п.3.8 описано отделение от сингулярности в интеграле из определяющих K (t) t соотношений. В нашем случае ядро в (2) имеет особенность при, а ряд t Прони ограничен при, поэтому для любых фиксированных Ci, i верно, что:

i K(t) e t t0 C , i i поэтому необходимо провести корректировку формулы (16), лотступив от точки Aexp(t) t 0. Так как K(t) , то необходимо раскладывать в ряд Прони лишь t1 1 't i функцию C 'e, тогда, с учетом отделения от t t 1. Пусть 1 i i сингулярности формула для вычисления интеграла из (2) примет вид i I(tn1) A Ci 'e( i ')(tn1tn )In Agk 1(t*)Jn, (17) i tn (1 e( i ')(tn tn1)) i ( i ')(tn ) i In g( )d e( i ')(tn tn1)I gk (t*n) где, ne ( i ') а Jn ) [ ] [, (tn1 tn )],[ ] t 1etdt,[, x] t 1etdt.

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

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

В п. 4.1.1 рассматривается задача об одноосном растяжении * (t) ( )однородной квадратной пластины (сторона 10, центр (0,0)) из 0 0 несжимаемого вязкоупругого материала. Механические свойства материала моделировались определяющими соотношениями (1) со значением констант 0.016 0.000167с 0 1МПа A 0.0135с1,,,. Задача решается при плоской деформации. Для этой задачи было найдено как численное, так и аналитическое решения и произведено их сравнение. Аналитическое решение данной задачи выглядит следующим образом:

{1 (1 xx2 1/ xx2)xx2} pn xx (18) {1 (1 xx2 1/ xx2) / xx2}xx2 pn xx {1 (1 xx2 1/ xx2)} pn zz * Из уравнений (18) аналитически была получена величина для 0,005i,i 0,1xx 1 0.1sin(50t) на отрезке [0, 0.5] в точках. Затем, для полученных граничных условий, было найдено численное решение для xx для 0.005 / j, j 1,10,100,10различных шагов по времени ( ).

Рис. 2.а Рис. 2.б На рис. 2 представлена зависимость величины относительной разности между численным и аналитическим решением Err (аналит числен ) / аналит 100 от xx xx xx времени: для а) j = 100, для б) j =1000. Как видно из рисунка 4, численное решение сходится к аналитическому решению при измельчении шага по времени. При этом максимальная относительная разность между решениями для j=1 равна 3,8464%, для j=10 равна 0.1508%, для j =100 равна 0,0138% и для j = 1000 равна 0.002224%.

В п. 4.1.2 решается серия задач о квазистатической деформации квадратной пластинки (сторона 10, центр (0,0)) с круговым отверстием (радиус 1, центр (0,0)) под воздействием растягивающих напряжений / 0 C, где yy C 0.01, 0.02, 0.05, 0.1. Механические свойства материала моделировались 0 1МПа определяющими соотношениями (1) со значением констант, 0.016 0.000167с A 0.0135с1,,. Задача решается при плоской деформации.

На рис. 3 представлена форма пластины и распределение компоненты yy yy тензора полных истинных напряжений для растяжения на границе 0.1 (рис.

1 0 0, рис. 3.б:

3.а: ).

Для данной задачи было проведено сравнение с численно-аналитическим решением, полученным К.М. Зингерманом с учетом нелинейностей только второго порядка. Сравнивались значения перемещений ux в точке пересечения x границы отверстия и оси.

Рис. 3.а Рис. 3.б 2 ux ux На рис. 4 изображена относительная разница D между решением max(ux ) ux ux, полученным с учетом нелинейностей только 2-го порядка и решением , полученным с учетом нелинейностей высших порядков.

Рис. 4.а Рис. 4.б Как видно из графиков, на рис. 4.a при малых деформациях (около 5%) учет нелинейностей выше 2-го несущественен ( D 1%). При больших деформациях (более 5%) разница между решениями начинает возрастать. В частности, при деформациях около 15% D 7%, а при деформациях в 30% (рис.4.б) D 32%. Этот факт показывает необходимость учета нелинейностей высших порядков при решении задач такого типа.

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

В п.4.2.3 решается серия задач о динамическом деформировании квадратной пластинки (сторона 10, центр (0,0)) с эллиптическим отверстием (большая полуось 0.6, эксцентриситет 3, центр (0,0)) под воздействием растягивающих напряжений yy 0.1 и дальнейшем динамическом образовании (t=3.8) эллиптического отверстия (большая полуось 0.6, экcцентриситет 3, центр 0.6,0.6 + i*0.2), где i=0, 2, 4, 6. Механические свойства материала моделировались определяющими 0.00 1МПа соотношениями (1) со значением констант,, 0.001с A 0.0139с1,. Задача решается при плоской деформации.

Рис. 5.а. Рис. 5.б.

На рис. 5 представлена форма пластины и распределение компоненты тензора полных истинных напряжений для различных моментов времени yy (5.a.: 0 0, 5.б.: 2 7 ) Рис. 6.а Рис. 6.б На рис. 6.а и рис. 6.б показаны зависимости компоненты тензора полных yy истинных напряжений для 2-х точек первой полости (6.а. - правой, 6.б. - левой).

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

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

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

3. Проведены численные эксперименты и получены следующие эффекты:

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

Публикации по теме диссертации 1. Левин В.А., Калинин В.В., Вершинин А.В., Пекарь Г.Е. Решение плоской задачи о концентраторе произвольной формы образованном в нагруженном теле. Конечные деформации // Известия Тульского госуниверситета. Серия "Дифференциальные уравнения и прикладные задачи". Ч Тула: изд-во ТуГУ, 2006. Ч Т. 12., Вып. 1. Ч С. 167-172.

2. Пекарь Г.Е. Об одной нестационарной задаче образования концентратора напряжений в нагруженном теле из несжимаемого вязкоупругого материала. Конечные деформации и их перераспределение // Вестник Нижегородского университета им. Н.И. Лобачевского. Ч Нижний Новгород, 2011. Ч № 4 (4) Ч С. 1689Ц163. Левин В. А., Вершинин А. В., Пекарь Г.Е., Саяхова Л. Ф., Труфен К.

Н., Филипенко Е. В., Яковлев М. Я. Использование нелокального критерия прочности в задачах теории многократного наложения больших деформаций // Ломоносовские чтения. Тезисы докладов научной конференции. Ч М., 2006. Ч С. 45.

4. Левин В. А., Пекарь Г. Е., Филипенко Е. В., Яковлев М. Я. Плоская задача об образовании полости произвольной формы в нагруженном теле из нелинейно-упругого материала. Конечные деформации // Ломоносовские чтения.

Тезисы докладов научной конференции. Ч М., 2007. Ч С. 106Ц107.

5. Левин В.А., Пекарь Г.Е. О модели образования дефекта в нагруженном теле из вязкоупругого материала. Конечные деформации // Проблемы шин и резинокордных композитов (Материалы восемнадцатого симпозиума). Ч М.: ООО Научно-технический центр НИИШП, 2007. Ч Т. 2, С. 42-6. Пекарь Г.Е., Левин В.А. Вариант алгоритма решения на базе мкэ плоской задачи об образовании концентратора напряжений произвольной формы в предварительно нагруженном теле из вязкоупругого несжимаемого материала.

Конечные деформации // Труды всероссийской научно-практической конференции Математика, информатика, естествознание в экономике и обществе. Секция физики, техники и электроники. Ч М., 2007. Ч С. 26.

7. Левин В.А., Калинин В.В., Агапов Н.А., Кукушкин А.В., Саяхова Л.Ф., Труфен К.Н., Пекарь Г.Е. Моделирование взаимовлияния нановключений в рамках механики деформируемого твердого тела с помощью пакета ABAQUS // Тезисы доклада "Современные проблемы математики, механики, информатики".

Тула, 2005. Ч С. 227.

8. Левин В.А., Вершинин А.В., Пекарь Г.Е., Саяхова Л.Ф., Труфен К.Н.

Некоторые возможности использования многофункционального специализированного программного комплекса Наложение // Ломоносовские чтения. Тезисы докладов научной конференции. Ч М., 2007. Ч С. 48.

9. Вершинин А. В., Калинин В. В., Пекарь Г. Е.Некоторые результаты решения плоской задачи о концентраторе напряжений произвольной формы, образованном в нагруженном теле. Конечные деформации // Материалы международной научной конференции "Современные проблемы математики, механики, информатики". Ч Тула, 2006. Ч С. 106Ц107.

10. Левин В. А., Вершинин А. В., Пекарь Г. Е.К решению плоской задачи о принудительном образовании дефекта в нагруженном вязкоупругом теле.Конечные деформации.//Материалы международной научной конференции "Современные проблемы математики, механики, информатики". Ч Тула, 2008. Ч С. 249Ц251.

11. Левин В.А., Вершинин А.В., Пекарь Г.Е., Фрейман Е.И. О разработке совместимого специализированного программного комплекса "НАЛОЖЕНИЕ", предназначенного для учета изменения нагрузок, дефектов и свойств материала в процессе нагружения при больших деформациях // Инженерные системы Ч 2009.

материалы международной научно-практической конференции. Ч М., 2009. Ч Т 1. Ч C. 125Ц132.

12. Левин В.А., Пекарь Г.Е. Алгоритм и вариант программной реализации задачи об учете перераспределения в теле конечных деформаций (для различных типов определяющих соотношений) // Проблемы шин и резинокордных композитов (Материалы девятнадцатого симпозиума). Ч М.: ООО Научнотехнический центр НИИШП, 2008. Ч Т. 2, С. 67 - 13. Левин В.А., Пекарь Г.Е. Оценка взаимодействия и взаимовлияния концентраторов напряжений, последовательно образуемых в теле из вязкоупругого материала. Конечные деформации // Ломоносовские чтения.

Тезисы докладов научной конференции. Ч М., 2008. Ч С. 117Ц118.

14. Левин В.А., Пекарь Г.Е. Некоторые результаты решения задач об образовании концентраторов напряжений различной формы в телах из вязкоупругого материала. Конечные деформации // Ломоносовские чтения.

Тезисы докладов научной конференции. Ч М., 2009. Ч С. 103.

15. Левин В.А., Пекарь Г.Е. Разработка программного модуля для решения задач вязкоупругости. Слабосингулярные ядра // Ломоносовские чтения.

Тезисы докладов научной конференции. Ч М., 2010. Ч С. 124Ц125.

16. Пекарь Г.Е. Разработка модуля CAE Fidesys для вычисления интегралов свертки. Решение задачи об образовании отверстия в пластине из нелинейно вязкоупругого материала // Материалы международной научной конференции "Современные проблемы математики, механики, информатики". Ч Тула, 2011. Ч C. 175.

17. M. Charara, A. Vershinin, D. Sabitov, G. Pekar. Solving the 3D Acoustic Wave Equation with Higher-order Mass-lumped Tetrahedral Finite Elements // 73rd European Association of Geoscientists and Engineers Conference and Exhibition 2011.

Ч Vienna, Austria, 2011. Ч P. 41Ц45.

Заказ № 369-i/04/2012 Подписано в печать 27.04.12 Тираж 100 экз. Усл. п.л. 0,ООО УЦифровичокФ, тел. (495) 649-83-www.cfr.ru ; e-mail:zak@cfr.ru    Авторефераты по всем темам  >>  Авторефераты по техническим специальностям