Идентификация анизотропных материалов и моделирование процессов конечного деформирования гипоупругих тел тема автореферата и диссертации по механике, 01.02.04 ВАК РФ

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

Христич Дмитрий Викторович

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

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

АВТОРЕФЕРАТ 1 1 МАР 2015

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

Тула 2015

005560298

005560298

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

Научный консультант — доктор физико-математических наук,

профессор Маркин Алексей Александрович

Официальные оппоненты: Миронов Борис Гурьевич,

доктор физико-математических наук, профессор, ФГБОУ ВПО «Чувашский государственный педагогический университет им. ИЯ. Яковлева», ректор

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

Шоркин Владимир Сергеевич, доктор физико-математических наук, профессор,

ФГБОУ ВПО «Госуниверситет - УНПК», ведущий научный сотрудник

Ведущая организация - ФГБОУ ВПО «Воронежский

государственный университет»

Защита диссертации состоится мая 2015 г. в /^^часов на заседании диссертационного совета Д 212.271.02 при ФГБОУ ВПО «Тульский государственный университет» по адресу: 300012, г. Тула, пр. Ленина, 92 (12105).

С диссертацией можно ознакомиться в библиотеке и на сайте ФГБОУ ВПО «Тульский государственный университет» (http://tsu.tula.ru).

Автореферат разослан « Я. 5"» фФрраЛЯ 2015 г.

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

Толоконников Лев Алексеевич

Актуальность темы исследования и степень её разработанности.

Структура большого числа естественных и искусственных материалов обладает элементами симметрии, которые во многом определяют их физические свойства. Проблеме описания симметрии свойств анизотропных материалов посвящены многочисленные работы как зарубежных, так и отечественных авторов: К. Трусделла, А.Е. Грина и Дж. Адкинса, A.B. Шубникова, С.Г. Лехницкого, Л.Д. Ландау и Е.М. Лифшица, В.В. Лохина и Л.И. Седова, Ю.И. Сиротина, Я. Рыхлевского, К.Ф. Черныха, Н.И. Остросаблина. О необходимости проведения испытаний по определению типа анизотропии свойств материалов говорилось в монографии А.Е. Грина и Дж. Адкинса. В работах Б.Д. Аннина и Н.И. Остросаблина, И.Ю. Цвелодуба, М.А. Hayes, J.P. Jaric, A.N. Nords предлагалось идентифицировать тип симметрии свойств материала на основе преобразований тензора упругости, отнесённого к произвольной системе координат и содержащего 21 константу материала.

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

Конструкционные материалы, обладающие анизотропией свойств, могут работать в условиях повышенных температур и больших деформаций, проявляя при этом существенную нелинейность поведения. Построению термомеханических моделей конечного деформирования анизотропных тел посвящены работы А.Е. Грина и Дж. Адкинса, В.И. Левитаса, К.Ф. Черныха, A.C. Кравчука, A.A. Маркина и М.Ю. Соколовой. В этих работах определяющие соотношения модели строятся как обобщение на случай конечных деформаций известных линейных соотношений путём использования различных энергетически сопряжённых мер напряжений и конечных деформаций. Однако в настоящее время практически отсутствуют работы, содержащие экспериментально конкретизируемые модели термоупругого поведения.

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

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

Работа выполнялась в рамках грантов РФФИ № 10-01-97501-р_центр_а «Моделирование термомеханических процессов в анизотропных средах», 12-01-31176-мол_а «Термомеханические модели процессов конечного деформирования», 13-01-97501-р_центр_а «Большие деформации и разрушение

анизотропных тел», 14-01-31138-мол_а «Построение и экспериментальное обоснование определяющих соотношений нелинейной теории упругости», 14-31-50244-мол_нр «Экспериментальная методика определения упругих свойств анизотропных материалов», ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы, проект 2012-1.4-12-000-1004-6437, госзадания № 467 Минобрнауки России.

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

Основными задачами работы являются:

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

2) разработка программы экспериментов с макрообразцами по определению типа начальной упругой анизотропии материала;

3) анализ линейных и нелинейных упругих свойств квазикристаллов;

4) постановка и решение связанных задач неизотермического конечного деформирования анизотропных тел.

Научная новизна работы:

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

2) разработана программа экспериментов, позволяющая идентифицировать тип симметрии анизотропного материала;

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

4) для аксиальных квазикристаллов, гексагональных кристаллов и трансверсально-изотропных материалов определена структура тензоров шестого ранга, характеризующих их нелинейные упругие свойства. Исследована структура этих тензоров в зависимости от порождающих элементов группы симметрии анизотропных материалов;

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

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

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

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

Методология и методы исследования. Сформулированные задачи исследования решаются на основе применения термомеханического подхода к описанию процессов деформирования анизотропных твёрдых тел в рамках гипотезы сплошности. В качестве основы для формулировки термомеханически обоснованных нелинейных определяющих соотношений для конечного равновесного деформирования анизотропных материалов используется теория процессов A.A. Ильюшина, положения которой в работах ГЛ. Бровко, В.И. Левитаса, A.A. Маркина, A.A. Рогового, П.В. Трусова были распространены на случай конечного деформирования изотропных сред.

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

На защиту выносятся следующие положения:

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

2) термомеханические модели поведения гипоупругих изотропных и анизотропных материалов при конечных деформациях;

3) программа экспериментов по идентификации типа начальной анизотропии материала;

4) структура тензоров, описывающих линейные и нелинейные термоупругие свойства квазикристаллов и кристаллических анизотропных материалов;

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

Достоверность полученных результатов подтверждается:

1) использованием в качестве методологической основы исследований теории процессов A.A. Ильюшина и термомеханического подхода, введённого в механику сплошных сред Л.И. Седовым, A.A. Ильюшиным, И.И. Гольденблатом, У. Ноллом, К. Трусделлом;

2) исследованием симметрии свойств материала на основе классических работ A.B. Шубникова, Ю.И. Сиротина, Э. Спенсера, А. Грина, В.В. Лохина и Л.И. Седова, а также известных положений теории групп;

3) использованием надёжных численных методов решения краевых задач;

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

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

Апробация работы. Основные положения и результаты работы доложены на научном семинаре по механике деформируемого твёрдого тела им. JI.A. Толоконникова (научный руководитель - профессор A.A. Маркин, г. Тула, 2014 г.), на научном семинаре кафедры теории упругости МГУ имени М.В. Ломоносова (научный руководитель - профессор И.А. Кийко, г. Москва, 2012, 2014 гг.), на научном семинаре по механике деформируемого твёрдого тела НИИ механики МГУ имени М.В. Ломоносова (научный руководитель -академик РАН И.Г. Горячева, г. Москва, 2014 г.), на научном семинаре кафедры прикладной математики факультета гражданского строительства университета Баухаус (научный руководитель - профессор К. Гюрлебек, г. Веймар, Германия, 2014 г.), на международных научных конференциях «Современные проблемы математики, механики, информатики» (г. Тула, 2011, 2012, 2013, 2014 гг.), на XVIII Зимней школе по механике сплошных сред (г Пермь 2013 г.), на X Всероссийском съезде по фундаментальным проблемам теоретической и прикладной механики (г. Нижний Новгород, 2011 г.), на международной конференции «Актуальные проблемы прикладной математики, информатики и механики» (г. Воронеж, 2012 г.).

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

Структура и объём диссертации. Работа состоит из введения, пяти глав, заключения, списка использованных источников и двух приложений. Объём диссертации - 252 страницы. Работа содержит 27 рисунков, 3 таблицы и список использованных источников из 252 наименований.

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

В первой главе приведены основные соотношения, необходимые для описания процессов конечного деформирования анизотропных тел. Деформации сплошной среды описываются тензором-аффинором. В полярное разложение аффинора деформаций Ф = 1М1 = 1*-V входят левая и правая меры искажений и = 11т и V = V7, а также ортогональный тензор поворота Я (К )• На их основе строятся тензоры деформаций Коши-Грина

1 2

£ = 'г(и -Е), Генки Г = 1п11, линеаризованный тензор деформаций

Содержание работы

В работах ГЛ. Бровко, В.И. Левитаса, A.A. Маркина, П.В. Трусова процессы конечного деформирования сплошной среды описываются с, помощью неголономных мер деформаций различных типов. В данной работе

используется одна из них, определяемая из уравнения М = ^LT1 • U + Ü • 1Г1).

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

dV

сплошной среды: М-Е = 0 = 1п-. Девиатор М связан только с

dV0

формоизменением.

Для описания напряжённого состояния сплошной среды используются тензор истинных напряжений Коши S, обобщённый тензор Коши £ = effS, «повёрнутый» обобщённый тензор Коши £R =R E R"', тензор условных

напряжений (первый тензор Пиолы-Кирхгофа) Р = (ф~') энергетический

тензор напряжений (второй тензор Пиолы-Кирхгоффа) Т = • 2 • Ф 1.

Процессам деформирования и нагружения, определённым тензорами е = siJe'eJ и Т = Ttjä'P, ставятся в соответствие их векторные образы э = эа7а и ö = aja, « = 0,1,2,3,4,5 в шестимерном декартовом пространстве A.A. Ильюшина с ортонормированным базисом /'". Базисным векторам iа пространства Е6 соответствуют тензоры I" пространства Е3 такие, что е = эа1а для любого симметричного тензора второго ранга. Тензоры I" образуют обобщённый канонический тензорный базис A.A. Ильюшина и представляются диадными разложениями

^«-^(е'г'+в^+ё3?3), 1,=-^(2ё3г3-г,г,-г2г2), 12=-^(ё2ё2-ё1ё1), i'^fä'+gV), 1*=-^(ё2ё3+ё3ё2), (1)

На его основе строится базис тензоров четвёртого ранга +

В пространстве Е6 основное термомеханическое тождество, записанное с использованием удельной свободной энергии у = U — ST, имеет вид

dy/ + SdT = —¿г • (Б. Вектор напряжений и энтропия определяются Ро

соотношениями

дщ „ дц/ ч

Основное термомеханическое соотношение в форме Гиббса в

шестимерном пространстве записывается в виде сЮ = ——э-Ыа-БсГГ, где Б -

Ро

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

дв 0 ов

Среда является изотропной, если для любого ортогонального преобразования материального базиса х',=хг<2 выполняется требование

(Г1 • Т(х,/) • О = • е(х, г) • д]^ У<2 е 8,

(4)

где 3 - функционал, связывающий напряжения, деформации и температуру, g ~ полная ортогональная группа. Если условие (4) выполняется не при всех <2, то среда называется анизотропной. Анизотропная среда обладает симметрией свойств при выполнении условия

О"1 • Т(д=,0-<2 = з{[<Г! • е(3г,г)• <2,Г(5,0,г1 уд е 8л,

гДе 8а ~ подгруппа полной ортогональной группы, называемая группой анизотропии. Разные анизотропные материалы могут обладать симметрией свойств различного типа. Тип анизотропии (симметрии свойств) среды полностью определяется группой .

В линейной теории термоупругости зависимость тензора напряжений от тензора деформаций еЛ и температуры Т имеет вид

8 = М"£Л-В-(Г-Г0), (5)

где N - начальный тензор упругости четвёртого ранга; В - тензор коэффициентов температурных напряжений второго ранга. Обратные соотношения имеют вид: еЛ =С--8 + \-(Т-Т0), где С - тензор упругих податливостей четвёртого ранга, тензор второго ранга А образован коэффициентами температурных расширений.

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

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

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

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

Выбирая различные пары энергетически сопряжённых тензоров напряжений и деформаций, связь (6) можно представить в виде T(/) = F£[E(/),7X0] или Sr(/) = Fa/[M(0,7'(0], где Fe и FM — некоторые функции. Использование в определяющих соотношениях мер деформаций б и М, инвариантных относительно жёсткого поворота, позволяет тождественно удовлетворить требованию материальной объективности. Прямая термомеханическая задача состоит в установлении связи (6).

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

= (7)

то есть e(t) = Ff1 [T(f),Г(/)] или M(t)=FÁ'/ [Ек(/),Г(0] • Установление связи (7) составляет обратную термомеханическую задачу.

В линейной теории термоупругости решениями прямой и обратной термомеханических задач являются соотношения Дюгамеля-Неймана (5) и обратные им соотношения, а в изотермических процессах - закон Гука, который в пространстве Е6 записывается в виде линейных соотношений

¿? = пэ, э = ссх, (8)

где п, с — симметричные тензоры второго ранга, являющиеся шестимерными образами тензоров N и С. Связь между тензором четвёртого ранга N в трёхмерном пространстве и его образом п задаётся соотношениями

N0ki=ßünaf>ß"< "ap=PÏNijkiPu> в которых матрицы ß° и ßl являются известными из работ A.A. Ильюшина.

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

_1

которого N ■ иа = Ла<аа, саа =(с„ ••£„) 2ва или в шестимерном пространстве п-За= ЛаЗа, причём вектор wa является шестимерным образом тензора <аа. В общем случае разложение тензора п по собственному базису представляется в виде

п =í¿ana, (9)

а=I

где п — число различных корней характеристического уравнения, базисные тензоры íia, соответствующие однократному корню Ла, имеют вид Sla = й>ата, а соответствующие корню кратности к — fïa -(baà>a + ¿>a+i<öa+i +...+ ä>a+k-i<üa+k-i • Из выражения (9) и соотношений (8) следует разложение тензора упругих податливостей по собственным базисам

С=1т-",г- (Ю)

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

Из разложений (9) и (10) следует инвариантность тензоров п и с как относительно ортогональных преобразований, связанных с выбором начальной системы координат, так и относительно преобразований вращения и отражения в собственных подпространствах.

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

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

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

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

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

PoV = - Z naß3a3ß + X Grslr) - £ Ьаэа(Т-Г0) + р0у/(°\Т), (И)

0 у=\ сг=0

2

где эа, s(/) - линеиные и квадратичные инварианты деформаций, m и и -число этих инвариантов для рассматриваемого типа анизотропного материала,

(т т т Л

щ=-сеТ0 — 1п—-—+ 1 , сс - удельная теплоёмкость при постоянных

Ч^о 'о )

деформациях.

Если считать входящие в выражение (11) коэффициенты константами материала, то в соответствии с выражениями (2) получим

т-1 - п т—1

<7= Е ПаРэр1а + X2(7%,- I Ьа1а(Т-Тп), (12)

ог,£=0 у=I а=0

1 т-\ Т

Р0а=0 'о

Соотношения (12) при бесконечно малых деформациях в изотермических процессах совпадают с законом Гука. Для изотропного материала соотношения (12) принимают вид а = /1°°э0/о + 20'э(|), где первое слагаемое определяет

гидростатические напряжения, а второе описывает напряжения при формоизменении. По аналогии с изотропным материалом назовём входящие в выражение (12) константы п"*3 модулями упругости, С - обобщёнными модулями сдвига. Константы Ьа характеризуют температурные напряжения.

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

Т = N • -е - В(Г - Г0), (13)

Ек=1Ч"М-В(Г-:Го). (14)

В изотермических процессах выражения (13), (14) принимают вид

Т = 1Ч-е, (15)

= N • -М. (16)

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

Полагая в представлении свободной энергии (11) коэффициенты , С функциями соответствующих линейных и квадратичных инвариантов: и0Г/,=па/'(эй,э1>...,з)Я_,), С =ОгЦ2г)), получим следующую термомеханическую модель анизотропного материала, удовлетворяющую предельной форме частного постулата:

ду/ ^

ОЭ

(8паР я V » ( ~ - Л

эаэр + Эр -Ъа(Т- Т0) к + 21

дЭа ) 7=1

1 т_1 Т

Ро а=0 'о

8Нг) )

Чг)' (17)

В соответствии с выражением (17) процесс в каждом неодномерном собственном подпространстве э(г) не зависит от процессов в других

собственных подпространствах: если э = э(/), то а = <?(г); если э = эаТа, то и

В работе рассмотрен вариант соотношений (17) для случая, когда п°0(эа,э/}) = П^ + Яа'>(эа + э/)), а где п«3, я", С",

С - константы:

о ^=14 2 у

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

Взаимное влияние процессов, происходящих в различных собственных подпространствах, можно учесть, полагая, что входящие в выражение для свободной энергии (11) коэффициенты пар, Сг зависят от всех инвариантов деформаций, определённых для данного типа материала: ПаР =я"/?(э0,э1,...,эт_„5(21),5(22),...,5(2л)), с =О''(Э0,Э1,...,Эш_1,5(21),^2),...,5(2л)). В этом случае напряжения определяются соотношениями

т-1 а=0

т-1 /7=0

2 — Ьа(Т — Т0) к + 2±0%) +

(18)

эаэр + Т.

ё-дОг 2 )-

\г)

т-1

2 —Т-э^Эд + 22 -^-¿оо |э(9) -

аг,/?=0 СЯ/. -' 1 4

а<»>

íff-' Г-!&(„

ч2о.,7=0 сЬр " Р г=1 '

Если рассмотреть процесс э = 5(1) = Э3Г3 в трансверсально-изотропном материале, то из соотношений (18) определяются напряжения

М- дб1 2- дв1 2 т

&(1>

ас1 2 -5(оэ(1) -

ас? 2г ас' 2т

20х + д°Х-2

аЭз2

Тэз Рз'з

которые имеют составляющую, ортогональную э(1), поэтому определяющие

соотношения (18) удовлетворяют обобщению частного постулата, но не удовлетворяют предельной форме обобщения.

Для изотропного материала соотношения (18) имеют вид

а —

00 15л00 2 8вХ 2

/п +

26' +

а«00 2

д4)

Эп+2

ас'

(0 * (19)

Формоизменение (э0=0, приводит к появлению

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

дэ0

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

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

1 т-1 _ и , т-1

А>С = ~ I саРоаор - £ - £ ааоа{Т-Т0) -с„раТ0

а,0=0 г=I о-=0

т, т т ,

—1п---+ 1

тп тп Тп

'О 'о 10

где оа, I* - линейные и квадратичные инварианты напряжений, с"^, I?, аа -

константы материала, са — удельная теплоёмкость при постоянных

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

т-1 а п т-1 „

э= Е саРац1а + Е 2Бгз{г) + X аа>а(Т-Г0), О у=1 ог=0

1 т-1 Т

8 = —^аааа+с1Т\п—.

Ро а=о 'о

Эти определяющие соотношения являются решениями обратной термомеханической задачи.

В третьей главе разработана программа идентификации типа начальной упругой анизотропии материала. Знание типа анизотропии требуется для разработки экспериментальных программ по конкретизации материальных констант и функций, которые входят в определяющие соотношения, построенные во второй главе. В работах Н.И. Остросаблина, И.Ю. Цвелодуба и других авторов идентификация типа симметрии свойств материала выполняется на основе преобразований тензора упругости, отнесённого к произвольной системе координат и содержащего 21 константу материала. В диссертации предлагается программа, которая не требует предварительного определения упругих постоянных.

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

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

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

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

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

Если е1#Е2*е3, то направления главных осей анизотропии а1, а2, а3 определяются однозначно. Канонические оси анизотропии Р, Р, Р в этом случае совпадают с главными осями анизотропии, а материал по своим свойствам может быть ромбическим, моноклинным или триклинным.

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

Если е1=е2фе3, то однозначно определяется главная поворотная ось а3, а ортогональные ей векторы а1, а2 выбираются произвольно. Базисы а1, а2, а и к , к2, Р связаны ортогональным тензором поворота <13 = со5<р(Рк1 +РР) + ьтр(РР -РР) + РР, так что а' = к'-<}3, / = 1,2,3,

а материал может быть тригональным, тетрагональным или гексагональным.

Угол <р в гексагональном материале может быть произвольным, так как упругие свойства такого материала инвариантны относительно любых поворотов вокруг вектора Р = а3.

Для определения угла <р для тригонального и тетрагонального материалов требуется два эксперимента: растяжение-сжатие вдоль векторов ах,

а2 с тензором напряжений ¡4(а1а1 - а1 а1) и сдвиг в плоскости этих векторов с тензором напряжений /5(а'а2 + а2а1), который можно осуществить, выполняя растяжение-сжатие под углом 45° к направлениям а1, а2. Обозначим измеренные в каждом из этих опытов компоненты тензора деформаций £у и е^ соответственно.

Из эксперимента на двухосное растяжение-сжатие тригональный материал можно отличить от тетрагонального и гексагонального по наличию отличных от нуля компонент деформаций и Ё2ъ- По данным этого эксперимента угол <р для тригонального материала определяется по формуле

3 гг23

Тетрагональный материал отличается от гексагонального тем, что в опыте на двухосное растяжение-сжатие для первого е!2 * О ПРИ а Для

второго £|2=0. При этом для обоих типов материалов е13=еп=0. В случае тетрагонального материала угол <р определяется по формулам

<3("" = -1агС18_ или ^">=-1 агс18_2*"(? .

4 4 *Ц'5-£12'4

Если вычисленный по указанным формулам угол <р = 0, то канонические

оси анизотропии к' совпадают с главными осями анизотропии а'.

По результатам экспериментов на двухосное растяжение-сжатие и сдвиг в плоскости векторов а1 = к1, а2 = к2 тетрагональный материал отличается от гексагонального следующим образом: если

то материал гексагональный, а если

gn ^ gl2

'4 h '

то материал является тетрагональным.

—1 —2 —3

Если Е1=е2=£3, то в качестве главных осей анизотропии а , а , а следует выбрать оси лабораторной системы координат, в которой проводился эксперимент. Материал может быть изотропным или кубическим. В кубическом материале взаимная ориентация векторных базисов а' и к' определяется ортогональным тензором поворота Q = q^k-'k1 : а' = к' • Q, / = 1,2,3. Показано, что для определения положения канонических осей анизотропии к' относительно осей лабораторной системы а' достаточно провести два эксперимента на одноосное сжатие в лабораторной системе координат. Компоненты тензора Q определяются с помощью численного решения системы уравнений.

Если одна из главных осей анизотропии совпадает с канонической осью анизотропии кубического материала, например, аъ=къ, то тензор поворота имеет вид <3 = <23 =со$д>{к1к1 + кгк2^+$т<р{кук2-А2Р) + РА3, и угол <р

определяется выражением <р(к) = — ак^ &22_ £эз, где - компоненты тензора

2 еп

деформаций при сжатии образца вдоль вектора а1.

Если по результатам описанных экспериментов тензор то

материал является кубическим. Если же <2 = Е, то главные и канонические оси анизотропии совпадают, и для идентификации типа материала требуются дополнительные эксперименты.

Для того, чтобы отличить изотропный материал от кубического, требуется провести эксперимент на двухосное растяжение-сжатие по направлениям векторов а1, а2 с тензором напряжений 14(ага1 - а2а2) и эксперимент на сдвиг в этой плоскости с тензором напряжений /5(а1а2 + а2а'). Идентификацию типа изотропного и кубического материалов в канонических

осях анизотропии можно выполнить по критерию (20): если £ц. = £12.) то

и '5

материал изотропный, а если — * , то материал является кубическим.

и (5

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

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

Идентификация ромбического, моноклинного и триклинного материалов также может быть произведена в экспериментах по кручению сплошных круговых цилиндров моментами Мх = Мхах, Мг= М2а2, Л/3=А/3а3, где Л/,, М2, М3 - величины приложенных моментов. При этом ось цилиндра направлена вдоль одного из векторов а1, а2, а3. Определение типа материала

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

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

В четвёртой главе исследуются термомеханические свойства квазикристаллов - материалов, которые не имеют периодической кристаллической структуры, но обладают дальним порядком апериодического типа. Квазикристаллы не обладают трансляционной симметрией, поэтому могут иметь поворотные оси симметрии 5, 7, 8, 10-го и более высоких порядков, недопустимые для периодически упорядоченных кристаллов.

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

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

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

гз 2лг

квазикристалла являются поворот вокруг оси к на угол —

= со5^(*'Р + РР) + 5Ш^(РР -РР) + РР (21)

и поворот

<2 = Н (22)

где тензор К = собог^РР + РР) + зтаг(РР-РР) + РР задаёт ориентацию оси третьего порядка относительно выбранной декартовой системы координат,

причём угол а определяется из геометрических соображений так, что 5 + 2у[5 . /ю-2л/5

———, Бюа = Л———, а тензор

Кз3 * = ~^ёгё1 + ё2ё2^ + ~(ё1ё2-ё2ёх} + ё3ё3 задаёт поворот на угол

вокруг поворотной оси третьего порядка.

Инвариантными относительно преобразований (21) и (22) являются базисный тензор второго ранга 1°, базисный тензор четвёртого ранга I00 и линейная комбинация тензоров 1п+122 +133+ 144 +155. Канонические представления тензоров А и С для икосаэдрических квазикристаллов записываются в виде А = а01°, С=с00100+с11^1п + 122 + 133+144+155| и

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

Аксиальные квазикристаллы имеют одну поворотную ось симметрии порядка п = 5 или п>6. В плоскостях, перпендикулярных оси симметрии, атомы расположены квазипериодически. Квазипериодические плоскости упакованы периодическим образом вдоль оси симметрии. Элементами симметрии аксиальных квазикристаллов являются поворотная ось симметрии порядка п, совпадающая с направлением вектора А3, и плоскость симметрии, перпендикулярная этой оси. Порождающими элементами группы симметрии такого квазикристалла являются ортогональный тензор поворота

+ + (23)

и 4 ' п * '

и ортогональный тензор, задающий отражение относительно плоскости симметрии,

= РР + РР - РР. (24)

Для аксиальных квазикристаллов при любом значении п = 5 или п>6 инвариантными относительно обоих преобразований (23), (24) являются базисные тензоры второго ранга 1°, I1 и тензоры четвёртого ранга I00, I01, I11,

а также линейные комбинации базисных тензоров 122 + 133, 144 +155. Канонические представления тензоров А и С для квазикристалла, имеющего поворотную ось пятого порядка и перпендикулярную ей плоскость симметрии, записываются в виде С = с00100+с01101+сп1и+с22(122 + 133)+С44(144 + 155),

А = (2ц10+я,11. Разложения тензоров А и С по каноническому тензорному базису совпадают с каноническими представлениями этих тензоров для гексагонального материала. Это значит, что в отношении термоупругих свойств аксиальные квазикристаллы ведут себя как кристаллы, относящиеся к гексагональной сингонии, или трансверсально-изотропный материал.

cosa =

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

* г. 3V

будут характеризоваться тензором шестого ранга D = —у, симметричным по

дг

первой, второй, третьей парам индексов и их перестановкам:

Dijklnm ~ Djiklmn = Ay/tow = Цijklnm ~ ^ijmnU &klijmn •

Применим метод построения инвариантных комбинаций базисных тензоров для определения структуры тензора D, описывающего нелинейные упругие свойства аксиальных квазикристаллов, гексагонального и трансверсально-изотропного материалов. Аналогично тензорному базису полусимметричных тензоров четвёртого ранга laß построим ортонормированный базис тензоров шестого ранга

\aßr = ^{\alßlr + VVlß + \ß\aV + lßl4a + Vlalß + \r\ßla)b (25)

симметричных по первой, второй, третьей парам индексов и их перестановкам.

Тензор шестого ранга D однозначно представляется разложением по базису (25) в виде D = Daprlaßr. В это разложение должны войти инвариантные

относительно преобразований (23) и (24) базисные тензоры \aßr. Инвариантными относительно преобразования Q являются базисные тензоры,

для которых выполняются условия =laßr. Непосредственная проверка

выполнения этого условия показывает, что каноническое представление тензора шестого ранга D, описывающего нелинейные упругие свойства квазикристалла, имеющего поворотную ось порядка п = 5 или п > 6 и перпендикулярную ей плоскость симметрии, имеет вид

D = D^I000 + Dml°°l + Ami0" + Aul'" + A,22(i°22 +1°33)+

+ Аж^044 +1°55) + Dl2i{lm + 1Ш) + A« (i144 +(26)

+ Z)244 (i244 " I255 - 2I345) + ü245 (2I245 + I344 - I355).

Порождающими элементами группы симметрии гексагонального

—г л

материала являются ортогональный тензор поворота вокруг оси к на угол —

и ортогональный тензор отражения относительно плоскости, перпендикулярной оси к3, <}х (24). Каноническое представление тензора шестого ранга Б, описывающего нелинейные упругие свойства гексагонального материала, имеет вид

Л = Сода!000 + АюГ1 + 4,, Г* + Ап1Ш + 2 +1°33) +

+Й044(1М4 + 1055)+А22(1122 +1133) + О144(1144+1155) + А22(1222-31233)+ (27)

+Опз (31223 -I333) + Д244 (I244 -1255 - 21345) + £>245 (21245 +1344 -1355).

Каноническое представление тензора шестого ранга Б, описывающего нелинейные упругие свойства трансверсально-изотропного материала, совпадает с представлением (26).

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

Если порождающим элементом группы симметрии анизотропного материала служит только ортогональный тензор поворота вокруг оси к3, то для квазикристаллов с поворотной осью симметрии пятого порядка каноническое представление тензора шестого ранга Б отличается от разложения (26), так как оно будет включать все двенадцать тензоров \арг и их линейных комбинаций, инвариантных относительно преобразования поворота (21):

О = Ахх>100° + Аю!*00' + ¿>01.1011 + А1II''1 + ^»022 (I022 +1°33) +

+£>044 (I044 +1055) + £»,22 (I122 +1133) + я144 (I144 + 1155) + +^224 (I224 + 21235 -1334 ) + п225 (I225 - 21234 -1335 ) +

+¿>244(I244 -1255 - 21345)+ Я245 (21245 +1344 -1355).

Для квазикристаллов с поворотной осью симметрии порядка п>6 и трансверсально-изотропного материала при отсутствии плоскости симметрии каноническое представление тензора шестого ранга Б имеет вид (26). В случае гексагонального материала при отсутствии плоскости симметрии каноническое представление тензора Б имеет вид (27). При этом наборы инвариантных базисных тензоров шестого ранга у гексагонального материала и квазикристалла с поворотной осью симметрии пятого порядка без плоскости симметрии различны.

Таким образом, тензоры второго и четвёртого рангов, описывающие линейные термоупругие свойства, имеют одинаковую структуру для гексагональных кристаллов, трансверсально-изотропных материалов и аксиальных квазикристаллов с осью симметрии порядка /7 = 5 или п> 6. Каноническое представление тензора Б, описывающего квадратичную зависимость напряжений от деформаций, для гексагональных материалов содержит двенадцать независимых компонент, а для трансверсально-изотропных материалов и аксиальных квазикристаллов с поворотной осью симметрии и-го порядка и перпендикулярной плоскостью симметрии - десять независимых компонент, причём представления тензора Б для этих материалов одинаковы. Тензор Э для трансверсально-изотропных материалов и аксиальных квазикристаллов с поворотной осью симметрии порядка п> 6 имеет одинаковый вид при наличии и при отсутствии плоскости симметрии. Для аксиальных квазикристаллов с поворотной осью симметрии 5-го порядка каноническое представление тензора О в этом случае содержит двенадцать независимых компонент, но отличается от представления тензора Б для гексагонального материала набором инвариантных базисных тензоров.

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

о

равновесия в отсчётной конфигурации V- Р + Р = б с использованием кинематических соотношений и связи между тензорами напряжений получено условие равновесного протекания процесса деформирования в отсчётной конфигурации

(и 'У -Ек • И + и-1 -¿к • Л +1Г1 -6(уУ)</К0 =

= / Р0 ■ 5ййГ0 +\Р-ШУ0, (28) ¿о Уо

где Р0 = Ид ■ Р - скорость изменения вектора внешней нагрузки, приложенной

на внешней поверхности 2Г0 с вектором единичной внешней нормали Ид, Р -скорость изменения вектора массовых сил.

Для конкретизации вариационного уравнения равновесия (28) для гипоупругого материала используются тензорно-линейные определяющие соотношения (14). Из этих соотношений при постоянных N и В следует скоростная форма записи определяющих соотношений для гипоупругого материала

±к =1Ч--М-ВГ. (29)

В работе получено уравнение теплопроводности в вариационной форме

(о '

\ (В ■ -МТ + с€рйт)8Тс!Уй = -1 щ ■ - | Л0 • УГ

К г„ кЛ

N (о \

■8 УТ

Ч У

<1Уй, (30)

где % - вектор теплового потока, Л0 - отнесённый к начальному состоянию тензор теплопроводности с положительно определённой матрицей компонент. В уравнении (30) учитывается влияние деформационных характеристик процесса на температурное поле.

Основными уравнениями системы связанной термомеханической задачи являются: вариационное условие равновесного протекания процесса деформирования (28); определяющие соотношения в форме (29); уравнение теплопроводности в вариационной форме (30); эволюционные соотношения

М = |(и_1-и+й-и-1), и • и + и ■ й = (^ у) • Фт + Ф • (V ,

адо-^вр, = \/хеК0;

начальные условия

й(х,/0) = й0(х), 2к(х,/0) = 2:(к)(х), 7\х,/0) = Г0(х) VxeK0; граничные условия статического, кинематического или смешанного типа

Р0 = Ро(х,() Ухе£р (31)

й = м*(х,/) Ухе2Г„ V/ >Г0, (32)

причём Го = 27, и £„ и Ери, Ер П Еи = 0, Ер П = 0, £я П £ри = 0;

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

Т = Т*(х,1) Ухе2Гг У/>/0, (34)

Чо=Чо(Я,0 Ухе£д (35)

О

^•Л0-УГ + «.(Г-Ги) = 0 УхеГс У/>/0, (36)

где а - коэффициент теплообмена. Поверхности Ет, £ч, £с не пересекаются: = ¿г и 2Г? II2ГС, ЕтГ[£ч=0, £тГ\£с = 0, £дГ\£с = 0. При задании

граничных условий полагаем, что функции (31Н34) являются дифференцируемыми функциями времени.

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

С использованием разработанной программы выполнен расчёт напряжённо-деформированного состояния в полом прямом круговом толстостенном (.КвИШиар = 0,75) изотропном цилиндре под действием внутреннего давления при постоянной температуре Т0= 300К без учёта

23

влияния деформаций на температурное поле. Граничные условия задачи: г = Явн: = />0«» = />«=0, Т = 0; г = Янар: Р^ = Р™ = /><« = 0, Г = 0; г = 0: аг = иф = иг -0, Т = 0; г = Ь\ иг = иф = иг = 0, 71 = 0; р - скорость изменения давления.

Расчёты проводились до деформаций 20%. В этих расчётах интерес представляет изменение формы цилиндра с ростом давления. Геометрически нелинейная постановка задачи позволяет проследить за образованием и развитием «бочкообразное™». На рисунке 1 изображены формы расчётной области (половины осевого сечения цилиндра) при различных уровнях деформаций, полученные в процессе расчёта.

а) 6) в) г)

Рисунок 1 - Формы расчётной области для цилиндра при различных деформациях: а) 5%; б) 10%; в) 15%; г) 20%

Анализ полученных результатов показывает, что линейное решение совпадает с нелинейным до уровня деформаций 1,5%. При деформации 20% наибольшая величина радиального перемещения точки на внутренней поверхности цилиндра, рассчитанная по нелинейной теории, превышает соответствующую величину, рассчитанную по линейной теории, на 17%.

Решена задача о конечном деформировании полого прямого кругового толстостенного (Лв„ / =0,5) изотропного цилиндра при следующих

граничных условиях: = Нйн: Р$г) = р1, Р0(ф) = Я0(г) = 0, Твн=500К; г = ЯШ1р :

= Ф) = роМ=0, Тнар=ЗООК; 2 = 0: Р^ = 0, «ф=«г=0,

ЧоГ) = <7оФ) = 9ог) = 0; * = /Ь(Л) = 0, Мф = «2= 0, д(0г) = <#> =0, где р -скорость изменения давления, Твп, Т - температуры окружающей среды,

контактирующей с внутренней и наружной поверхностями цилиндра. Начальная температура цилиндра Г0 = 300К.

Решения найдены в двух вариантах постановки задачи: связанной и несвязанной. Поля напряжений, деформаций и температуры не зависят от осевой координаты. Результаты решения представлены на рисунках 2-7. Сплошные линии на этих рисунках соответствуют решению связанной задачи, штриховые - несвязанной. Графики построены для моментов времени 0,4с, 0,8с, 1,2с, 1,6с, когда внутреннее давление достигает величин соответственно 0,Ю, 0,20, 0,30, 0,4й, где О — модуль сдвига материала цилиндра.

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

Рисунок 5 - Распределение температуры в цилиндре

Рисунок 4 - Распределение осевых напряжений в цилиндре

Рисунок 6 - Распределение радиальных Рисунок 7 - Распределение окружных

деформаций в цилиндре деформаций в цилиндре

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

Целью расчёта является определение напряжённо-деформированного состояния цилиндрического баллона с внутренним радиусом Лви = 0,076 м и наружным радиусом Яиар = 0,1 м под действием внутреннего давления,

находящегося в неоднородном нестационарном температурном поле. Решается связанная термомеханическая задача. Начальная температура во всех точках цилиндра Т0 = 273К. Граничные условия (31>—(36) в этой задаче имеют вид:

г = Явн: = />« = />« = 0, Т(!) = 70 + у/; г = Янар: = = = 0,

Т(1) = Т0 + уГ, 2 = 0: = 0, и9 = и2= 0, <#>=^>=^> = 0; 2 = ь, РМ = 0,

мФ = мг=0, д(0г} = = = 0, где р - скорость изменения давления, у -скорость изменения температуры.

Результаты расчётов напряжённо-деформированного состояния композитных баллонов представлены на рисунках 8-11. На этих рисунках показаны распределения напряжений, отнесённых к величине внутреннего давления, по толщине стенки баллона. Размеры баллона во всех расчётах были одинаковыми. Расчёты проводились для четырёх вариантов числа слоев, угла намотки в и скорости изменения температуры на поверхности баллона: 1) четыре двойных спиральных слоя, 6 = ±30°, у = 0,5К/с; 2) четыре двойных

спиральных слоя, 0 = ±45°, у = 0,5К/с; 3) четыре двойных спиральных слоя, 0 = ±45°, у = 1К/с; 4) восемь двойных спиральных слоев, 9 = ±45°, у = 1К/с.

Рисунок 8 -- Распределение окружных Рисунок 9 - Распределение осевых

напряжений по толщине стенки баллона напряжений по толщине стенки баллона

Расчёты по предлагаемой методике позволяют найти касательные напряжения возникающие между слоями композита из-за стремления

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

Рисунок 10 - Распределение касательных Рисунок 11 - Распределение касательных напряжений по толщине стенки баллона напряжений по толщине стенки баллона

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

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

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

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

Схема нагружения шара на второй стадии приведена на рисунке 12: г, = 9 мм, г2 - 8,16 мм, Р = 60°. Начальная температура во всех точках Г0 = 293К.

Граничные условия (31)—(36) при рассматриваемой схеме нагружения принимают вид: хе Р = —р(п$, Тн = 2000К; дге2Г2: у„=0, Рх= 0, ГИ=293К; хе£3: р(г) = Р{(ф) = ри) =о,7;= 293К.

Результаты расчётов напряжённо-

деформированного состояния резинового шара представлены на рисунках 13-16. На этих рисунках показаны распределения напряжений, отнесённых к модулю сдвига О.

Рисунок 12-Схема Д° приложения давления (при р = 0) в шаре

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

а)

б)

в)

Рисунок 13 - Распределения напряжений в шаре при р = 0: а) радиальные напряжения; б) окружные напряжения; в) осевые напряжения

Рисунок 14 - Распределения напряжений в шаре при р = 0,4960: а) радиальные напряжения; б) окружные напряжения; в) осевые напряжения

Рисунок 15 - Распределения напряжений в шаре при р = 0,8320 : а) радиальные напряжения; б) окружные напряжения; в) осевые напряжения

а) б) в)

Рисунок 16-Распределения напряжений в шаре при р = \,П6С\ а) радиальные напряжения; б) окружные напряжения; в) осевые напряжения

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

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

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

Основные результаты работы:

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

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

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

3. Разработана программа экспериментов по определению типа начальной анизотропии материала.

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

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

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

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

Содержание диссертации отражено в 42 публикациях автора, основными из которых являются следующие:

1. Маркин, A.A. Процессы упругопластического конечного деформирования [Текст] / A.A. Маркин, М.Ю. Соколова, Д.В. Христич. - Тула: Иэд-во ТулГУ, 2011. - 374 с.

2. Маркин, A.A. Постулат А.А.Ильюшина для анизотропных материалов и вариант определяющих соотношений [Текст! / A.A. Маркин, М.Ю. Соколова, Д.В. Христич // Известия РАН. Механика твёрдого тела. -2011. - № 1. - С. 38-45.

3. Соколова, М.Ю. Модель упругопластического деформирования нелинейных анизотропных материалов |Текст] / М.Ю. Соколова, Д.В. Христич И Известия ТулГУ. Естественные науки. Вып. 2. Часть 2. - Тула: Изд-во ТулГУ, 2013. - С. 239-250.

4. Соколова, М.Ю. О симметрии термоупругих свойств квазикристаллов [Текст] / М.Ю. Соколова, Д.В. Христич // Прикладная математика и механика. - 2014. - Т. 78. - Вып. 5. - С. 728-734.

5. Соколова, М.Ю. Описание конечных деформаций твёрдых тел в отсчётной конфигурации [Текст] / М.Ю. Соколова, Д.В. Христич // Прикладная механика и техническая физика. - 2012. - Т. 53 - № 2 - С 156-166.

6. Христич, Д.В. Аналитическое определение симметрии свойств квазикристаллов |Текст] / Д.В. Христич // Известия ТулГУ. Естественные науки. Вып. 1. - Тула: Изд-во ТулГУ, 2012. - С. 81-S8.

7. Христич, Д.В. Варианты нелинейной связи между напряжениями и деформациями в анизотропных материалах [Текст] / Д.В. Христич //

Известия ТулГУ. Естественные науки. Вып. 1. Часть 1. - Тула: Изд-во ТулГУ, 2014. - С. 216-224.

8. Христич, Д.В. К вопросу об определении главных осей анизотропии материала [Текст] / Д.В. Христич // Известия ТулГУ. Естественные науки. Вып. 2. - Тула: Изд-во ТулГУ, 2014. - С. 203-213.

9. Христич, Д.В. Компьютерное моделирование экспериментов по определению типа начальной анизотропии упругих материалов [Текст] / Д.В. Христич // Известия ТулГУ. Естественные науки. Вып. 4. - Тула: Изд-во ТулГУ, 2014. - С. 110-119.

10. Христич, Д.В. Критерий экспериментальной идентификации гексагонального, тригонального и тетрагонального материалов [Текст] / Д.В. Христич // Вестник КГТУ им. А.Н.Туполева. - 2013. - № 2. - С. 67-72.

11. Христич, Д.В. Критерий экспериментальной идентификации изотропного и кубического материалов [Текст] / Д.В. Христич И Известия ТулГУ. Естественные науки. Вып. 3. - Тула: Изд-во ТулГУ, 2012. - С. 110118.

12. Христич, Д.В. Критерий экспериментальной идентификации ромбического, моноклинного и триклинного материалов [Текст] / Д.В. Христич // Известия ТулГУ. Естественные науки. Вып. 3. - Тула: Изд-во ТулГУ, 2013. - С. 166-178.

13. Христич, Д.В. Моделирование процесса конечного деформирования анизотропных тел [Текст] / Д.В. Христич // Вычислительная механика сплошных сред. - 2013. — Т. 6, № 4. - С. 410-419.

14. Христич, Д.В. Нелинейные упругие свойства анизотропных кристаллов и аксиальных квазикристаллов [Текст] / Д.В. Христич // Известия ТулГУ. Естественные науки. Вып. 1. - Тула: Изд-во ТулГУ, 2013. -С. 111-122.

15. Христич, Д.В. Постановка задачи конечного деформирования анизотропных тел в терминах начальной конфигурации [Текст] / Д.В. Христич, Ю.В. Астапов, Л.В. Глаголев // Известия ТулГУ. Естественные науки. Вып. 3. - Тула: Изд-во ТулГУ, 2014. - С. 148-157.

16. Христич, Д.В. Программа экспериментов по определению главных осей анизотропии материала [Текст] / Д.В. Христич, P.A. Каюмов, И.З. Мухамедова // Известия КГ АСУ. - 2012. -Ä3 (21). - С. 216-224.

17. Христич, Д.В. Решение краевых задач нелинейной термоупругости [Текст] / Д.В. Христич, М.Ю. Соколова // Известия ТулГУ. Естественные науки. Вып. 1. - Тула: Изд-во ТулГУ, 2010. - С. 123-136.

18. Соколова, М.Ю. Нелинейные задачи анизотропной термоупругости [Текст] / М.Ю. Соколова, Д.В. Христич // IX Всероссийский съезд по теоретической и прикладной механике. Аннотации докладов. Т. III (Нижний Новгород, 22-28 августа 2006). - Нижний Новгород: Изд-во Нижегородского госуниверситета им. Н.И.Лобачевского, 2006. - С. 197-198.

Христич Д.В.

Автореферат

Изд.лицЛР № 020300 от 12.02.97. Подписано в печать 11.02.2015 Формат бумаги 60x84 Бумага офсетная. Усл.печ. л. 1,9 Уч.изд. л. 1,6 Тираж 100 экз. Заказ 003 Тульский государственный университет 300012, г. Тула, просп.Ленина, 92. Отпечатано в Издательстве ТулГУ 300012, г. Тула, проспЛенина, 95.