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

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

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

Удк 539.3 перераспределение конечных деформаций, вызванное образованием концентраторов напряжений

МОСКОВСКИЙ ОРДЕНА ЛЕНИНА, ОРДЕНА ОКТЯБРЬСКОЙ РЕВОЛЮЦИИ

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

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

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

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

ВЕРШИНИН АНАТОЛИЙ ВИКТОРОВИЧ

УДК 539.3

ПЕРЕРАСПРЕДЕЛЕНИЕ КОНЕЧНЫХ ДЕФОРМАЦИЙ, ВЫЗВАННОЕ

ОБРАЗОВАНИЕМ КОНЦЕНТРАТОРОВ НАПРЯЖЕНИЙ

01.02.04 – механика деформируемого твердого тела

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

Москва – 2007 г.

2

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

Научный руководитель – доктор физико-математических наук, профессор Левин Владимир Анатольевич

Официальные оппоненты – доктор физико-математических наук, профессор Бровко Георгий Леонидович – доктор технических наук, профессор Морозов Евгений Михайлович

Ведущая организация – Институт прикладной механики РАН

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

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

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

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

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

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

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

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

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

Вклад в развитие как нелинейной теории упругости, так и эксперимента для нее внесли многие отечественные и зарубежные специалисты, в частности, Г.М.Бартенев, В.Л.Бидерман, В.Д.Бондарь, М.Ф.Бухина, И.И.Ворович, Н.В.Зволинский, Л.М.Зубов, Ю.А.Крутков, Л.И.Кутилин, А.И.Лурье, Н.Ф.Морозов, В.В.Новожилов, В.А.Пальмов, П.М.Риз, Г.Н.Савин, Л.И.Седов, Л.А.Толоконников, Т.Н.Хазанович, К.Ф.Черных, 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 и многие другие. Развитие теории наложения малых деформаций на большие началось с середины 60-х годов прошлого века, наиболее подробно вопросы теории были рассмотрены в работах киевской школы механиков под руководством А.Н. Гузя. Создание и развитие теории многократного наложения больших деформаций для тел из упругого или вязкоупругого материала было осуществлено Г.С. Тарасьевым, В.А. Левиным. В работах В.А. Левина рассмотрены также вопросы зарождения и развития дефектов в рамках механики деформируемого твердого тела при конечных деформациях, предложены нелокальные критерии прочности и модели, учитывающие возникновение и развитие зон предразрушения (совместно с Е.М. Морозовым), разработаны методы оценки эффективных характеристик пористых материалов при конечных деформациях и их наложении (совместно с В.В. Лохиным, К.М. Зингерманом). Задача наложения деформаций для двухконстантного потенциала (с учетом ряда упрощающих допущений) была рассмотрена Л.М. Нечаевым. Одним из современных подходов к решению таких задач можно считать использование метода конечных элементов, а также систем компьютерной алгебры для получения приближенных аналитических решений.

Численное решение рассматриваемых задач может быть найдено с использованием МКЭ в совокупности с методом Галёркина. Данный метод был предложен в 1915 г. Б.Г. Галёркиным как приближенный метод решения краевых задач. Ранее, в 1913г., метод применялся для решения конкретных задач теории упругости И.Г. Бубновым, в связи с чем именуется также методом Бубнова - Галёркина. Теоретическое обоснование метода принадлежит М. В. Келдышу (1942). Применение метода конечных элементов к задачам линейной и нелинейной теории упругости подробно рассмотрено в работах Л. Дж. Сегерлинда, О.Зенкевича, Одена.

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

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

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

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

– разработка программного обеспечения для реализации указанного алгоритма на базе МКЭ;

а также получение приближенного аналитического решения задачи об образовании эллиптического жесткого включения с использованием системы компьютерной алгебры «Mathematica 5.0».

Научная новизна.

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

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

Практическая значимость заключается в постановке плоских и пространственных задач об образовании и развитии концентраторов напряжений в предварительно нагруженных телах, а также программной реализации алгоритма решения указанных задач. Разработанный программный комплекс использовался при выполнении работ по гранту РФФИ (проект № 06-01-00682)1.

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

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

Апробация работы. Основные положения диссертационной работы докладывались и обсуждались на международной научной конференции, посвященной 80-летию со дня рождения профессора Л.А.Толоконникова в 2003 г. (г. Тула);

на Международной научной конференции «Современные проблемы математики, механики, информатики» в 2004, 2005 гг. (г. Тула);

на Международной конференции «Фундаментальные проблемы разработки нефтегазовых месторождений, добычи и транспортировки углеводородного сырья» в 2004 г. (г. Москва);

на пятнадцатом, шестнадцатом и семнадцатом симпозиумах «Проблемы шин и резинокордных композитов» в 2004, 2005, 2006 гг. (г. Москва);

на 6-й научно-технической конференции «Актуальные проблемы состояния и развития нефтегазового комплекса России» в 2005 г.

(г. Москва);

на научных конференциях «Ломоносовские чтения» в 2004, и 2006 гг. в МГУ им. М.В. Ломоносова;

на научных семинарах кафедры «Вычислительная механика» (под руководством чл.-корр. РАН А.В. Забродина);

на шестом Всероссийском семинаре “Сеточные методы для краевых задач и приложения” в 2005 г. (г. Казань);

на седьмом международном конгрессе по вычислительной механике в 2006 г. (Лос Анджелес, США).

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

Результаты работы использовались при написании книги В.А. Левин, В.В. Калинин, К.М. Зингерман, А.В.Вершинин “Развитие дефектов при конечных деформациях. Компьютерное и физическое моделирование”. М. ФИЗМАТЛИТ. 2007. 392с. (в печати) Структура и объём работы. Диссертация состоит из введения, трех глав, заключения, приложения и списка литературы. Работа изложена на страницах машинописного текста, содержит 65 рисунков, список использованных источников из 117 наименований.

КРАТКОЕ ИЗЛОЖЕНИЕ ДИССЕРТАЦИИ Во введении обоснована актуальность темы диссертационной работы, сформулированы цели и задачи данной работы, приведена аннотация содержания диссертации.

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

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

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

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

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

k k = эi i – оператор градиента;

q n n m n = I + up = I up m,n q, p = q,k k, p, – аффинор p = m +1 p = m + деформаций ( m n );

E q, p = ( m, p T, p m,q T,q ) – тензор деформаций, описывающий m m m изменение деформаций при переходе тела из состояния q в состояние p и отнесенный к координатному базису m -го состояния;

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

Fq, p = T, p q, p – тензорная мера деформаций, описывающая q изменение деформаций при переходе тела из состояния q в состояние p и соответствующая мере Фингера ( F0,1 – тензорная мера Фингера);

n, f n – плотность и массовая сила в n -м состоянии;

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

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

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

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

m n 0,n = T,m 0,n n,m ;

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

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

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

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

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

(Индекс над символом, кроме э и R, указывает номер состояния, в координатном базисе которого вычисляется данная величина).

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

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

Сжимаемый материал (на примере материала Мурнагана) 0 0 0 0,n = E 0,n I I + 2G E 0,n + 3C3 E 0,n I I + 0 0 0 +C4 E 0,n I I + 2C4 E 0,n I E 0,n + 3C5 E 0,n, (1) Несжимаемый материал (на примере материала Муни) µ {(1 + ) F } 0,n = + (1 ) ( F0,n I ) F0,n F0,n p0,n I, (2) 0, n В п.1.4. приведены уравнения равновесия и граничные условия в координатном базисе текущего состояния, а также для тела, находящегося в n -м состоянии, в координатах произвольного k -го состояния ( k n ).

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

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

(Рис. 1а) - начальное состояние. В теле отсутствуют напряжения и деформации (Рис. 1б) – первое промежуточное состояние (Рис. 1в) – второе промежуточное состояние Затем в теле намечается некоторая замкнутая поверхность (рис. 1б), и удаляется часть тела, ограниченная этой поверхностью. Действие удаленной части тела на оставшуюся заменяется (по принципу освобождаемости от связей) силами, распределенными по этой поверхности. Далее эти силы, перешедшие в разряд внешних, «мгновенно» изменяются на большую величину, например, уменьшаются до нуля (под термином «мгновенно изменяются» при рассмотрении статических задач не следует понимать, что данное приложение (снятие) нагрузок приводит к деформированию тела в динамическом режиме). Тело, приобретая (теряя) большие дополнительные деформации и напряжения (по крайней мере, в окрестности вновь образованной граничной поверхности), переходит во второе промежуточное состояние (рис. 1в). Такое нагружение можно продолжить и дальше.

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

– уравнение равновесия:

k k k k k k 0,n 0,n ln (1 + 0,k ) + 0,n : ( k,n ) k,n k k ( T k,n 1 ) T k,n 0,n = (3) – уравнение несжимаемости (для несжимаемых материалов):

0,n = 0, (4) – граничные условия:

k (n) k k N n 0,n = QNn (5) k k n k k (n) n (n) |do| Q N n = (1 + 0,k ) P N n k,n n k n |do| n (n) P N n - вектор истинных напряжений на элементарной площадке n n n n d o = N n | d o |. В частности, в случае отсутствия нагрузок на границе k k концентратора получаем N n 0,n = 0.

k n 0,n = 0,n, (6) где 0,n – тензор истинных напряжений на бесконечности, определенный в постановке задачи.

В постановку задачи также входят уравнения, связывающие тензор истинных напряжений 0,n с аффинором деформаций 0,n. Для различных материалов они имеют вид (1), (2).

Завершают постановку задачи геометрические соотношения:

1 + 0,n = det 0,n, (7) 0,n = 0,k k,n, (8) n k u.

k,n = I + (9) i i = k + Решение этой задачи позволяет, в частности, найти суммарный вектор n k ui перемещений из k -го состояния в n -е как функцию радиус-вектора R, i = k + т.е. в координатах k -го состояния.

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

Постановка задачи включает:

– уравнение равновесия n 0,n = 0, (10) – уравнение несжимаемости (для несжимаемых материалов) 0,n = 0, (11) – граничные условия n (n) n N n 0,n = P N n, n (12) n n В частности, в случае отсутствия нагрузок на границе концентратора n получаем N n 0,n = 0.

n n 0,n = 0,n, (13) – геометрические соотношения 1 + 0,n = det 0,n, (14) nn = I ui.

0,n (15) i = Уравнения, связывающие тензор истинных напряжений 0,n с аффинором деформаций 0,n, – соотношения (1), (2).

Решение этой задачи позволяет, в частности, найти суммарный вектор n u перемещений из начального в n -е состояние как функцию радиус i i = n вектора R, т.е. в координатах n -го состояния.

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

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

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

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

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

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

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

Применение метода Синьорини (метода последовательных приближений) и системы компьютерной алгебры «Mathematica 5.0» к исследуемым задачам позволяет найти их приближенное аналитическое решение. При этом решение исходной нелинейной задачи сводится к решению бесконечной последовательности линеаризованных задач.

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

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

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

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

v N i ( 0,n )dS = n N i 0,n d ( Ni ) 0,n dS = 0;

N i - функции формы.

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

В п.2.4 рассмотрены математическая постановка динамической задачи теории наложения и ее решение с использованием МКЭ. Постановка задачи включает:

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

d 2u n n [(1 + 0,n ) 0,n +1 (t ) n,n +1 (t )] n 2 = (16) dt Начальные условия:

du u (0) = un ;

(0) = 0 (17) dt Граничные условия:

0,n +1 (t ) = 0,n +1 (t ) (18) n n n = P(t ) (1 + 0,n +1 (t )) N n +1 (t ) *,n +11 (t ) n,n +11 (t ) (19) N n +1 (t ) 0,n +1 (t ) n k n+1 ( t ) Определяющие соотношения для потенциала Мурнагана:

n 0,n +1 (t ) = (1 + 0,n +1 (t )) *,n +11 (t ) 0,n +1 (t ) 1n +1 (t ) (20) n n, 0 0 0 0,n +1 (t ) = ( E 0,n +1 (t ) : I ) I + 2G E 0,n +1 (t ) + 3C3 ( E 0,n +1 (t ) : I ) 2 I + (21) 0 0 0 +C4 ( E 0,n +1 (t ) 2 : I ) I + 2C4 ( E 0,n +1 (t ) : I ) E 0, n +1 (t ) + 3C5 ( E 0,n +1 (t )) 0,n +1 (t ) = (1 + 0,n +1 (t )) (t ) 0, n +1 (t ) 0,n +1 (t ) * (22) 0, n + E0,n +1 (t ) = ( 0,n +1 (t ) * n +1 (t ) I ) (23) 0, Геометрические соотношения:

1 + 0,n +1 (t ) = det 0,n +1 (t ) (24) 0,n +1 (t ) = 0, n n,n +1 (t ) (25) n n,n +1 (t ) = I + u (t ) (26) n = (27) 0,n Система уравнений Галёркина, соответствующая приведенной дифференциальной постановке, выглядит следующим образом:

r€ d 2u r € M 2 + K (u ) f = 0 (28) dt r € Здесь u - вектор неизвестных узловых перемещений;

M = N T Nd, (29) r r r r n € € € € K (u ) = N T [(1 + 0, n (u )) 1 0,n +1 (u ) n, n +1 (u )]d, (30) r N nd f= T (31) Решение полученной системы ищется по методу SS22.

В п.2.5 рассматривается применение метода Синьорини к решению задачи об образовании жесткого эллиптического включения в нагруженном теле. Решение исходной нелинейной задачи ищется в перемещениях в виде бесконечной суммы u2 = u2 + u2 +..., где u2i ) ~ q i +1.

(0) (1) ( (32) Вектор u2i ) (i = 1,2,...) называют поправкой от учета эффектов (i + 1) -го ( порядка для перемещений при переходе из 1-го во 2-е состояние. При таком подходе выделяется безразмерная величина q, определяемая следующим | P|, где max = max 0,2ij. Далее в подобном виде образом: q = max max, G G i, j представляются и другие характеристики напряженно-деформированного состояния. Подставив эти разложения в уравнения, входящие в постановку задачи, и сгруппировав члены одинакового порядка, можно получить бесконечную последовательность систем линеаризованных краевых задач для расчета напряжений и деформаций в теле в конечном состоянии, при последовательном решении которых можно найти сначала u2 и 0,2 (0) (0) (нулевое приближение), затем u2 и 0,2 (первое приближение) и т.д. В (1) (1) настоящей работе выполнены расчеты для нулевого и первого приближений.

Далее рассматривается применение метода Колосова-Мусхелишвили к решению линеаризованной задачи упругости. При этом вводится тензор S, соответствующий тензору напряжений линейной упругости. Для сжимаемого материала S = L1[u ], для несжимаемого S = L2 [u, p]. Тогда векторы u, f, N и тензор S могут быть представлены в координатной форме следующим образом:

u = u1e1 + u2e2 + u3e3, (33) f = f1e1 + f 2 e2, (34) N = N1e1 + N 2 e2, (35) S = S11e1e1 + S12 e1e2 + S21e2 e1 + S22 e2 e2 + S33e3e3. (36) В рассмотрение также вводятся:

– комплексные переменные z = x1 + ix2, z = x1 ix2, (37) – функции этих переменных:

w( z, z ) = u1 + iu2, F ( z, z ) = ( f1 + if 2 ), (38) N ( z, z ) = N1 + iN 2, (39) – комбинации компонент некоторого тензора T второго ранга (в декартовой системе координат):

TI = T11 + T22 + i (T12 T21 ), TII = T11 T22 + i (T12 + T21 ) (40) Для случая плоской деформации и сжимаемого материала уравнения равновесия и граничные условия линеаризованной задачи могут быть записаны в комплексной форме следующим образом:

S II S I + = 2F, (41) z z “условие прилипания” wвкл = wматр (42) =, S II =, (43) SI I II w w w SI = 2 ( + G ) +, S II = 4G. (44) z z z Решение линеаризованной краевой задачи ищется в виде:

w = wн + wодн, S = Sн + Sодн, (45) p = pн + pодн, 33 = 33 н + 33 одн, (46) где wн, Sн, pн, 33 н – некоторое частное решение линеаризованной задачи, wодн, Sодн, pодн, 33 одн – решение линеаризованной задачи для однородной системы уравнений.

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

Для случая плоской деформации сжимаемого материала эти формулы имеют вид:

1 ( + 3G ) Fdzdz ( + G ) Fdzdz, wн = (47) 4G ( + 2G ) w w w S I н = 2 ( + G ) н + н, S II н = 4G н. (48) z z z Рассматривается подход к решению линеаризованной краевой задачи для однородной системы уравнений с использованием комплексных потенциалов Колосова-Мусхелишвили ( z ) и ( z ) :

( z )dz z ( z ) ( z )dz, wодн = (49) 2G S I одн = 2 ( z ) + ( z ), S II одн = 2 ( z ) + z ( z ). (50) В третьей главе рассматриваются результаты решения задач, постановка и методы решения которых, приведены во второй главе.

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

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

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

В п.3.1 рассматриваются примеры решенных плоских статических и динамических задач.

В 3.1.1 рассмотрена задача об одновременном образовании двух эллиптических отверстий (отношение полуосей 1:5). На приведенных ниже рисунках для случая начального всестороннего растяжения при / G = 2.097, С3/G=-0.0689, С4/G=-0.37464, С5/G=0.33709;

0,111 /G=p, 0,122 /G=p;

p/G=0. приведено распределение главной компоненты тензора полных истинных напряжений 0,21 /G. Свойства материала матрицы описываются потенциалом Мурнагана. Отношение большой полуоси отверстия к стороне матрицы 1:20.

Схема нагружения Распределение первой главной компоненты тензора истинных напряжений между полостями ( 0,21 /G) В 3.1.3 рассмотрена задача об одновременном образовании двух круговых отверстий. На приведенных ниже рисунках для случая начального одноосного растяжения при = -0.263905, 0,111 = 0;

0,122 = p;

p/ µ =0. приведено распределение главной компоненты тензора полных истинных напряжений 0,21 /G. Свойства материала матрицы описываются потенциалом Муни. Отношение большой полуоси отверстия к стороне матрицы 1:20.

Распределение первой главной компоненты тензора истинных напряжений между полостями ( 0,21 /G) 0. 0. 0. on G 0. 0. 0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 x R В 3.1.5 приведено приближенное аналитическое решение методом Синьорини задачи об образовании в бесконечно протяженном теле жесткого эллиптического включения. На приведенных ниже рисунках для случая p p начальных нагрузок на бесконечности 0,11,1 = p, 0,11,2 =, 0,12,2 = ;

2 p = со следующими значениями параметров: µ = 2.51228, = 2.63905;

µ 11 приведено распределение и матрицы при удалении от µ matrix µ matrix включения. Свойства материала матрицы описываются потенциалом Муни.

S11 S ЂЂЂ ЂЂ ЂЂЂ ЂЂЂ ЂЂ ЂЂЂ m m 0. 0. 0.15 0. 0. 0. 0. y y 0.5 1.5 2 2.5 0.5 1.5 2 2.5 -0.05 0. В 3.1.7 рассмотрена динамическая задача об одновременном образовании (с учетом динамических эффектов) в предварительно нагруженном теле двух круговых отверстий. На приведенных ниже рисунках для случая начального всестороннего растяжения при / G = 2.097, С3/G= 0,1 /G=p, 0,1 /G=p;

p/G=0. 0.0689, С4/G=-0.37464, С5/G=0.33709;

11 приведено распределение первой главной компоненты тензора полных истинных напряжений 0,21 /G в различные моменты времени, начиная с момента образования. Свойства материала матрицы описываются потенциалом Мурнагана.

В п.3.2 рассматриваются примеры решенных трехмерных задач. В 3.2. рассматривается задача об образовании шарового включения с последующим образованием внутренних слоев. Ниже приведены результаты решения задачи для случая образования трех слоев (трех « вложенных» включений) при начальном одноосном растяжении p=G/20. Радиусы включений заданы в момент их образования R3=0.8R2, R2=0.8R1. Отношение R1 к стороне матрицы принято равным 1:20. Приведены результаты решения (поле распределения полных истинных напряжений 0,n1 / G в плоскости xOy) для случая, когда материал тела (матрицы) - материал Мурнагана с константами / G = 2.097, С3/G=-0.0689, С4/G=-0.37464, С5/G=0.33709, материал первого включения материал Мурнагана /G1=2.24;

С3/G1=-1.96, С4/G1=3.61, С5/G1=-11.13 и G/G1=0.39. Материал второго включения соответствует материалу первого при пористости 0.06, а третьего включения соответствует материалу первого при пористости 0.18.

Схема нагружения изменение 0,n1 / G вдоль оси Ox поле распределения 0,n1 / G в плоскости xOy В п.3.3 анализируется влияние нелинейных эффектов и зависимость получаемого решения от размеров конечноэлементной сетки. Приводится сравнение решений, полученных разными методами. В 3.3.3 приводится сравнение линейных и нелинейных решений. Рассматривается влияние как физической, так и геометрической нелинейностей.

Ниже приведены результаты решения задачи о растяжении квадратной матрицы с круговым отверстием при одноосном растяжении p=G/20. На графике представлено изменение полных истинных напряжений 0,n1 / G вдоль оси Ox для случая, когда материал тела (матрицы) - материал Мурнагана с константами / G = 2.24, С3/G=-1.96, С4/G=3.61, С5/G=-11.13.

0. 0. 0. 0. on G _Guk on 0. G _Linon G 0. 0. 0. 0 1 2 3 4 5 6 7 8 9 x x1 x,, RRR -.-.-.- – линейное решение (закон Гука и малые деформации) - - - – нелинейное решение (закон Гука и конечные деформации) _ – нелинейное решение (потенциал Мурнагана и конечные деформации) Таким образом, разница между линейным и нелинейным решениями в точке максимальной концентрации напряжений составляет порядка 30%.

Кроме того, отметим смещение точки максимальной концентрации напряжений от вершины носика полости.

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

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

получены приближенные аналитические решения плоской задачи о жестком эллиптическом (в момент образования) включении для разных типов нелинейно-упругих материалов. Разработан алгоритм и программное обеспечение для решения этой задачи на базе системы компьютерной алгебры «Mathematica 5.0» с использованием метода Синьорини.

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

Публикации по теме диссертации 1. В. А. Левин, И. А. Мишин, А. В. Вершинин “Плоская задача об образовании включения в упругом нагруженном теле. Конечные деформации”, Вестник Московского Университета, серия 1, Математика. Механика. 2006 №1 с. 56-59.

2. В. А. Левин, А. В. Вершинин “О приближенном аналитическом решении плоской задачи о жестком эллиптическом включении, возникающем в нагруженном теле. Конечные деформации”, Известия ТулГУ, серия “Дифференциальные уравнения и прикладные задачи”, 2004, с. 151-157.

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

Материалы Шестого Всероссийского семинара. – Казань: Казанский государственный университет, 2005. - с. 157- 4. Василевский Ю., Вершинин А., Данилов А., Пленкин А. “Технология построения тетраэдральных сеток для областей, заданных в САПР” Матричные методы и технологии решения больших задач – Москва:

Институт вычислительной математики РАН, 2005. – c. 21-32.

5. В. А. Левин, И. А. Мишин, А. В. Вершинин “Решение задачи об упругом или жестком включении, возникающем в предварительно нагруженном теле из высокоэластичного материала, с помощью программного комплекса “Наложение”. 15-й симпозиум “Проблемы шин и резинокордных композитов”, г. Москва, 18-22 октября 2004г., с.33.

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

Москва, МГУ, тезисы научной конференции “Ломоносовские чтения 2004”, c. 41.

7. В. А. Левин, И. А. Артюхин, А. В. Вершинин, Д. А. Улькин “Вариант модели для описания напряженно деформированного состояния тела из высокоэластичного материала после фазового перехода. Конечные деформации”, Тезисы докладов международной научной конференции, посвященной 80-летию со дня рождения профессора Л.А.Толоконникова, стр. 184-185, Россия, Тула, 18-21 ноября 2003 г.

8. В. А. Левин, А. В. Вершинин “К оценке напряженного состояния вблизи жесткого эллиптического включения, возникшего в нагруженном нелинейно-упругом массиве. Конечные деформации”, материалы Международной конференции «Фундаментальные проблемы разработки нефтегазовых месторождений, добычи и транспортировки углеводородного сырья» г. Москва, 24-26 ноября 2004г., стр. 174.

9. В. А. Левин, А. В. Вершинин “О приближенном аналитическом решении плоской задачи о жестком эллиптическом включении, возникающем в нагруженном теле. Конечные деформации”, сборник Международной научной конференции "Современные проблемы математики, механики, информатики", Россия, Тула, 17-19 ноября г., стр. 107-108.

10. В. А. Левин, И. А. Мишин, Е. А. Рыбалка, А. В. Вершинин “К расчету напряженного состояния массива при образовании в нем полости или включения с помощью авторского программного комплекса для аналитических вычислений на ЭВМ “Наложение”. Конечные деформации”, 6-я научно-техническая конференция, посвященная 75 летию Российского государственного университета нефти и газа им. И.

М. Губкина, г. Москва, 26-27 января 2005 г., стр. 414-415.

11. А. В. Вершинин “Об использовании МКЭ для решения задач теории многократного наложения больших деформаций”, тезисы научной конференции “Ломоносовские чтения 2005”, 19 апреля 2005 г., г.Москва, МГУ, c. 43.

12. В.А. Левин, А.В. Вершинин “Плоская задача о разгрузке тела из нелинейно-упругого материала после возникновения в нем зоны предразрушения. Конечные деформации”. Сборник тезисов Международной научной конференции “Современные проблемы математики, механики, информатики”, Россия, Тула, 22-26 ноября г., стр. 223-224.

13. В.А. Левин, А.В. Вершинин, Е.И. Фрейман, Е.Д. Комолова “К оценке микронапряжений при моделировании свойств наноматериалов в рамках механики деформируемого твердого тела при конечных деформациях”. Сборник тезисов Международной научной конференции “Современные проблемы математики, механики, информатики”, Россия, Тула, 22-26 ноября 2005 г., стр. 224-225.

14. В.А. Левин, К.М. Зингерман, Е.В. Рыбалка, А.В. Вершинин “Сравнение конечно-элементного и приближенного аналитического решения некоторых плоских задач теории многократного наложения больших деформаций”. Сборник тезисов Международной научной конференции “Современные проблемы математики, механики, информатики”, Россия, Тула, 22-26 ноября 2005 г., стр. 225-226.

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

Москва, МГУ, с. 44-45.

16. В.А. Левин, А.В. Вершинин, Пекарь Г.Е., Саяхова Л.Ф., Труфен К.Н., Яковлев М.А. “Использование нелокального критерия прочности в задачах теории многократного наложения больших деформаций.

Решение модельных задач с помощью программного комплекса "Наложение". Тезисы научной конференции “Ломоносовские чтения 2006”, 18 апреля 2006 г., г. Москва, МГУ, с. 45.

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

Тезисы научной конференции “Ломоносовские чтения 2006”, 18 апреля 2006 г., г. Москва, МГУ, с.43-44.

18. Левин В.А., Вершинин А.В. “Программный комплекс “Наложение” на базе МКЭ”. Тезисы научной конференции “Ломоносовские чтения” 2006, 18 апреля 2006 г., г. Москва, МГУ, с. 44.

19. Левин В.А., Зингерман К.М., Вершинин А.В. “Программный комплекс для численного решения динамических задач об образовании повреждений в нагруженных нелинейно-упругих телах при наложении больших деформаций”. // Проблемы шин и резинокордных композитов.

Материалы 17-го симпозиума. - М.: Научно-технический центр "НИИШП", 2006. Т.1. С. 214-217.

20. В.А. Левин, В.В. Калинин, К.М. Зингерман, А.В. Вершинин Развитие дефектов при конечных деформациях. Компьютерное и физическое моделирование. М. ФИЗМАТЛИТ. 2007. 392с. (в печати) 21. В.А. Левин, В.В. Калинин, А.В. Вершинин, Г.Е. Пекарь Решение плоской задачи о концентраторе напряжений произвольной формы, образованном в нагруженном теле. Конечные деформации. Известия Тульского госуниверситета. Серия "Дифференциальные уравнения и прикладные задачи". Том 12. Вып. 1. 2006г. с. 167- 22. В.А. Левин, А.В. Вершинин Авторский программный комплекс наложение и его использование для решения плоских и пространственных статических и динамических задач для предварительно нагруженных тел при конечных деформациях.

Известия Тульского госуниверситета. Серия "Дифференциальные уравнения и прикладные задачи". Том 12. Вып. 1. 2006г. с. 145-

 




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

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