WWW.DISSERS.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА

   Добро пожаловать!


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

Чигинский Дмитрий Сергеевич

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

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

АВТОРЕФЕРАТ

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

Тверь 2012

Работа выполнена в ФГБОУ ВПО «Тульский государственный университет».

Научный консультант: доктор технических наук, профессор Трещев Александр Анатольевич

Официальные оппоненты: Гараников Валерий Владимирович доктор технических наук, профессор, ФГБОУ ВПО «ТвГТУ», зав. каф. «Техническая механика» Гордон Владимир Александрович доктор технических наук, профессор, ФГБОУ ВПО «Госуниверситет – УНПК», зав. каф. «Высшая математика»

Ведущая организация: ФГУП «ГНПП «Сплав» (г. Тула)

Защита состоится « 31 » мая 2012 г. в ___:___ часов на заседании диссертационного совета Д 212.262.02 при ФГБОУ ВПО «Тверской государственный технический университет» по адресу: 170026, г. Тверь, наб. Афанасия Никитина, 22, ауд. Ц-.

С диссертацией можно ознакомиться в библиотеке ФГБОУ ВПО «Тверской государственный технический университет».

Автореферат разослан «___» _________ 2012 г.

Ученый секретарь диссертационного совета Гультяев Вадим Иванович

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

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

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

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

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

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

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

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

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

• используя термодинамический потенциал Гиббса и потенциальные уравнения состояния для изотропных материалов, чувствительных к виду напряжённого состояния, предложенные в работах Матченко Н.М., Толоконникова Л.А. и Трещева А.А., получить феноменологические соотношения термоупругости для изотропных разносопротивляющихся материалов;

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

• выбрать и обосновать принятый метод решения прикладных задач, разработать алгоритм решения поставленных термомеханических задач об определении НДС тонких пластин из изотропных материалов с усложненными свойствами и программную реализацию на ЭВМ;

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

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

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

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

Методы исследования, использованные в диссертационной работе:

• общепринятые, строго обоснованные допущения и гипотезы теории расчета тонких пластин, базирующиеся на фундаментальных законах механики деформируемого твердого тела;

• итерационный метод «упругих решений», разработанный А.А. Ильюшиным, для решения нелинейных уравнений деформирования изотропных материалов.

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

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

• уравнения теории термоупругости для решения связанных задач по расчёту НДС тонких пластин в условиях термомеханического нагружения, изготовленных из разносопротивляющихся материалов, зависящими от вида напряженного состояния, с учетом влияния температуры на механические характеристики материала и напряженного состояния на распределение температуры в элементах конструкций;

• математическая модель, вариант алгоритма и программная реализация итерационного метода решения задач по расчёту НДС тонких пластин, изготовленных из материалов с «усложненными» свойствами, находящихся в условиях термомеханического нагружения;

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

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

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

Практическая значимость работы, выполненной в рамках госбюджетной НИР № 27.06 «Актуальные проблемы технологии строительных материалов и проектирования конструкций», заключается в следующих результатах:

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

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

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

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

Работа выполнена частично за счет средств гранта губернатора Тульской области в сфере науки и технологии 2008 года. Свидетельство №13-2008 от 16.01.2009.

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

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

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

• на 3-й; 7-й Международной конференции по проблемам горной промышленности, строительства и энергетики «Социально-экономические и экологические проблемы горной промышленности, строительства и энергетики» (г. Тула, ТулГУ, 2007; 2011 г.);

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

• на II-й; V-й Молодёжной научно-практической конференции студентов «Молодёжные инновации» (г. Тула, ТулГУ, 2008; 2011 г.);

• на VII-й Международной научно-технической конференции, посвященной 50летию Пензенского государственного университета архитектуры и строительства (г. Пенза, ПДЗ, 2008 г.);

• на III-й Магистерской НТК (г. Тула, ТулГУ, 2008 г.);

• на 9-й; 10-й; 12-й Международной конференции «Актуальные проблемы строительства и строительной индустрии» (г. Тула, ТулГУ, 2008; 2009; 2011 г.);

• на V-й Международной конференции «Надёжность и долговечность строительных материалов, конструкций и оснований фундаментов» (г. Волгоград, ВолгГАСУ, 2009 г.);

• на VII-й Международном научном симпозиуме «Проблемы прочности, пластичности и устойчивости в механике деформируемого твёрдого тела», посвящённом 80-летию со дня рождения заслуженного деятеля науки и техники РФ профессора В. Г. Зубчанинова (г. Тверь, ТГТУ, 2010 г.);

• на VII-й Всероссийской научно-технической конференции студентов, аспирантов и молодых ученых, посвященной 50-летию первого полета человека в космос «Молодёжь и наука» (г. Красноярск, Сиб. федер. ун-т., 2011 г.);

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

В полном объеме диссертация докладывалась 14 марта 2012 года на научном семинаре по МДТТ им. Л.А. Толоконникова при ФГБОУ ВПО «Тульский государственный университет», под руководством доктора физ.-мат. наук, профессора А.А.

Маркина, на расширенном заседании кафедры «ССМиК» Тульского государственного университета 2 апреля 2012 года, а также 4 апреля 2012 года на научном семинаре по МДТТ при ФГБОУ ВПО «Тверской государственный технический университет», под руководством доктора техн. наук, профессора В.Г. Зубчанинова.

Публикации. По теме диссертационной работы опубликовано 20 печатных работы. Основное содержание диссертации отражено в 13 статьях, в том числе 3 работы в изданиях рекомендуемых ВАК РФ.

Структура и объем работы. Диссертация состоит из введения, четырех разделов, заключения, списка литературы, включающего 154 наименования, приложений. Диссертация содержит 142 страницы основного текста, в том числе 20 рисунков, 3 таблицы и приложения на 40 страницах, включающие результаты, текст программы расчёта тонких пластин и документы о внедрении. Общий объём работы — 182 страницы.

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

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

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

Первую группу составляют модели, в основу которых положена зависимость механических характеристик материала от знаков возникающих напряжений или развивающихся деформаций. В рамках данной группы рассмотрены определяющие соотношения, предложенные С.А. Амбарцумяном, Р.М. Джонсом, П.Н. Ельчаниновым, М.И. Климовым, А.Ф. Макеевым, Д.А.Р. Нельсоном, И.Г. Овчинниковым, В.В. Петровым, Б.В. Пономаревым, Г.С. Шапиро, Г.В. Бригадирова, Н.М. Матченко, J.N. Reddy, F. Tabaddor, C.W. Bert и другими.

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

Модели соотношений третьей группы строятся посредством учета взаимного влияния изменения объема и формоизменения при помощи специфического представления деформаций разрыхления как части полных деформаций или формулировкой дилатационных зависимостей. Уравнения этой группы получены в теоретических и экспериментальных исследованиях К.А. Агахи, Д.Л. Быкова, С.С. Вялова, А.И. Козачевского, В.И. Кудашова, В.Н. Кузнецова, В.П. Устинова и других.

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

Русинко, В.Г. Зубчанинова, Н.Л. Охлопкова, В.В. Гараникова, R.M. Jones, D.A.R.

Nelson и других.

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

Пономарева, В.В. Пикуля, С.П. Ти-мошенко, А.А. Хачатряна, Х.М. Муштари, С.А.

Кузнецова, К.З. Галимова, М.С. Корнишина и других.

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

Жидкова, P.E. Hart, J.N. Reddy, C.W. Bert, Y.S. Hsu и других. Приведены материалы по ее развитию, дана классификация различных теорий, отмечены их достоинства и недостатки.

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

Анализ особенностей каждого подхода для сформулированной цели исследования показал, что особенный интерес, с этой точки зрения, представляет использование теории термоупругости для изотропных материалов, рассмотренной в работах Н.М. Матченко и А.А. Трещева. Эта теория лишена большинства недостатков обнаружившихся в рассмотренных моделях. В рамках закона теплопроводности Фурье и условий динамического равновесия, с использованием термодинамического потенциала Гиббса, в ней получены основные дифференциальные уравнения разномодульной теории термоупругости: уравнение теплопроводности, включающее в связанном случае учет влияния вида напряженного состояния и уравнения динамического равновесия. Обобщение указанной теории на случай тонких пластин позволяет получить необходимые соотношения для решения задач термоупругости для пластин из разносопротивляющихся материалов.

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

Потенциал деформаций W принимается как функция характеристик нормированного пространства №2:

W = W (,, S0), (1) где , — угол и фаза напряжений; S0 — модуль вектора полного напряжения на 2 октаэдрической площадке S0 = +, = ijij 3 — средние нормальные напряжения; = SijSij 3 — касательные напряжения на октаэдрической площадке;

Sij = ij -ij — компоненты девиатора напряжений, ij — символ Кронекера.

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

2 W = b1 + b3 + b2 + b4 + b5 cos3 , (2) ( ) ( ) где b1 = 1,5(A + 2C); b2 = 1,5(A - C); b3 = 1,5(B + 4D + 2E) / 3; b4 = 4,5(B - E) / 3;

b5 = 0,75 B - 2D + 2E 2 / 3 ; = S0 и = S0 — нормированные нормальные и ( ) касательные напряжения на октаэдрической площадке; cos3 — фазовый инвариант.

Потенциал деформаций W как функции характеристик нормированного пространства № 1 принимается в виде:

W = W (I, III,S). (3) Выражение (3) в рамках квазилинейной теории, представленное в виде степенного полинома от нормированных напряжений, до третьей степени включительно, после простых преобразований, может быть использовано в виде:

2 2 W = 0,5 A + B1 1 + A + B2 2 + A + B3 3 + C + E3 + ( ) ( ) ( ) [ (4) +D +2 12 + C + E1 + D +3 23 + C + E2 + D +3 13.

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

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

Воспользовавшись теорией Матченко Н.М. и Трещёва А.А. задачу построения общей теории термоупругих разносопротивляющихся материалов, необходимо отказаться от подхода, основанного на принципе суперпозиции, и использовать понятия неравновесной термодинамики. Рассматриваются малые изменения температуры, таo кие, что / To << 1. В этом случае зависимостью механических и теплофизических характеристик материала от температуры можно пренебречь. В качестве независимых параметров состояния выбраны компоненты тензора напряжений ij и температуры T. В этом случае удобно воспользоваться термодинамическим потенциалом Гиббса:

Г = Г(ij,T ) ; (5) dГ = -eijdij - LdT, (6) где L — плотность энтропии.

В силу потенциальности соотношений (5) справедливы выражения, соответствующие формулам Кастильяно:

еij = -Г ij ; L = - Г T. (7) Рассмотрим представление термодинамического потенциала Гиббса (5) в виде:

Г = Г(,, S0,T ) (8) Представим функцию (8) в виде степенного полинома от нормы второго пространства S0 с коэффициентами разложения, зависящими от качественных инвариантов, и от температуры. Тогда, используя методику нормированного пространст ва №2 и считая, что все чисто температурные слагаемые представлены одной функцией t ( ), получим:

2 2 - = b1 + b3 + b2 + b4 + b5 cos3 + t 0 bt1 + bt2 + bt1 . (9) ( ) ( ) ( ) ( )+ Используя характеристики нормированного пространства №1, потенциал Гиббса для разносопротивляющегося материала можно записать следующим образом:

Г = Г(I, III, S,T ), (10) Окончательно для плотности энергии Гиббса применительно к характеристикам нормированного пространства напряжений №1 получим:

-Г = Гt 0 At + Bt1 1 + At + Bt2 2 + At + Bt3 3 + ( ) ( ) ( ) ( )+ 2 2 +0.5 A + B1 1 + A + B2 2 + A + B3 3 + C + E3 + (11) ( ) ( ) ( ) [ +D + 2 12 + C + E1 + D +3 23 + C + E2 + D +3 13, () () () где At, Bt — температурные константы потенциала нормированного пространства №1, вычисляемые по результатам обработки опытов, связанных с нагреванием и охлаждением одноосно растягиваемых и одноосно сжимаемых образцов из разносопротивляющегося материала.

Для полного совпадения двух форм (9) и (11) необходимо чтобы между константами этих форм имелись следующие связи:

bt1 = 3Bt ; bt2 = 3At. (12) Выражение для функции Гt 0 при относительно малых изменениях темпе( ) ратуры Т0, как и в классической теории термоупругости имеет вид:

Гt 0 = C 0 2T0. (13) ( ) ( ) где C — теплоемкость материала.

Уравнения состояния для термоупругого материала наряду с зависимостями между компонентами тензоров деформаций eij и напряжений ij должны содержать выражения для плотности энтропии L. Указанные уравнения можно получить, применив к потенциалу Гиббса (9) операцию дифференцирования (7). Очевидно, что уравнения связи между компонентами тензоров напряжений и деформаций представляют собой тензорно-нелинейную зависимость и если в этих уравнениях отделить линейную от нелинейной части, то получим:

eij = 2b2ij 3 + 2 b1 -b2 ij 3 + bt2ij 3 + Rij, (14) ( ) где Rij = Rij (,, S0, ) нелинейные слагаемые (представлены в диссертации).

L = bt1 + bt2 + bt1 + dt dT. (15) ( ) Используя формулы Кастильяно применительно к потенциалу в форме (11), придем к соотношениям, связывающим главные напряжения и деформации [77]:

еk = Ak + C n + m + Rk ; k n m ; (16) ( ) ( ) где Rk = Rk (I, III, S, ) нелинейные слагаемые (представлены в диссертации).

L = At + Bt1 1 + At + Bt2 2 + At + Bt3 3 + dГt dT, (17) ( ) ( ) ( ) В соотношениях (14) и (16) линейные члены отделены от нелинейных Rij, Rk.

( ) Здесь следует иметь в виду, что характер квазилинейных функций не позволяет непосредственно перейти от уравнений вида eij = eij km к зависимостям ij = ij ekm, ( ) ( ) так как Rij = Rij km. Поэтому обращая линейные члены уравнений (14), получим ( ) тензора напряжений:

ij = D1 + D2 eij - 3D2eij - D3 ij - Hij, (18) ( ) где e = eijij 3 ; Hij = Hij,, S0, ( ) нелинейные функции компонентов тензора напряжений, D1, D2, D3 механические константы материала.

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

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

и Трещёва А.А. имеют вид:

а) закон изменения объема e = 3K0 + 3D0 + bt1 + bt2 ; (19) ( ) б) закон изменения формы Э = 1+ tg2 2G0 + 3D0 + bt1 ( );

(20) в) уравнение, определяющее фазовую характеристику tg = 3b5sin3 2G0 + D0 + bt1 S(3 );

(21) г) закон изменения плотности энтропии L = bt1 + bt2 + bt1 + С0 Т0, (22) ( ) Анализируя выражения законов (19)(22), приходим к выводу о том, что полученные определяющие соотношения для квазилинейных термоупругих материалов учитывают дилатационные свойства. С другой стороны, как следует из уравнений (20) и (21), изменение формы зависит от средних напряжений. Однако, если касательные напряжения = 0, то формоизменения не происходит Э = 0, т. е. при гидростатическом и температурном напряжении не возникает деформаций сдвига. Наличие подобных свойств у разносопротивляющихся материалов подтверждается в ряде теоретических и экспериментальных исследований.

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

Для вычисления констант потенциалов деформаций (9) и (11) используются результатами опытов на одноосное растяжение и одноосное сжатие. В работах Матченко Н.М. и Трещёва А.А. были аппроксимированы экспериментальные зависимости линейными диаграммами по способу наименьших квадратов на весь диапазон деформирования образцов до разрушения и вычислены константы потенциалов деформаций (9) и (11) с учетом трех вариантов априорных предпосылок. Для используемых при апробации материалов — графита марки АРВ механические константы материала представлены в табл. 1. Причем для констант Е, b3, b4, b5 приведены значения, соответствующие второму варианту априорных предпосылок, применяемые к условиям одноосных или плоских напряженных состояний при изгибе тонких пластин.

Таблица 1. Механические константы материала - - - Технические, МПа Потенциала (11), МПа-1 Потенциала (9), МПа-Материал Е+ Е- + - А·10-5 В·10-6 С·10-6 D·10-7 E·10-6 b 1·10-5 b 2·10-5 b 3·10-6 b 4·10-5 b 5·10-Графит АРВ 3750 6130 0,200 0,35 21,490 51,767 -55,215 18,815 0 15,670 40,517 51,349 13,449 29,3Параметром, характеризующим реакцию тела на связанные механические и температурные воздействия можно назвать коэффициент температурного расширения, значение которого изменяется по всему объёму пластинки и в процессе установления температуры, в зависимости от времени. Искомый параметр является переменным коэффициентом при функции температуры в уравнениях состояния (14):

ij = bt1ij 3 + bt2ij 3. (23) В отличие от традиционного представления постоянного коэффициента температурного расширения, выражение (23) включает слагаемое, учитывающее влияние напряжённого состояния на температурные деформации.

Уравнение теплопроводности получено используя закон теплопроводности Фурье, при предположении, что скорость распространения тепла бесконечно велика:

qi = -T,i, (24) где qi компоненты вектора теплового потока; коэффициент теплопроводности.

Для построения уравнений притока тепла используется аналогичная работам Матченко Н.М. и Трещёва А.А. методика. Локальная формула приращения энтропии:

ТL,t = -qi,i +U, (25) где U удельная мощность источников тепла.

Подставим уравнение (24) в (25) получим зависимости приращения энтропии от температуры:

ТL,t = Т,ii +U. (26) Дифференциальное уравнение теплопроводности получено при рассмотрении совместно с (26) выражения (15), продифференцированного по времени и умноженного на температуру.

3At,t +BtS,t T + CT,t -Т,ii -U = 0. (27) ( ) Считая изменения температуры малыми, в выражении температуры можно пренебречь членом T0 << 1. После линеаризации и замены переменных получим уравнение притока тепла:

3At,t +BtS,t T0 +U = 0. (28) ,ii -C0,t -( ) Выражение (28) в отличие от классического уравнения теплопроводности остается нелинейным, несмотря на проведенную линеаризацию по температуре. Кроме традиционного слагаемого, определяющего связанность полей деформаций и температуры, в уравнение (28) входит нелинейная компонента, учитывающая влияние вида напряженного состояния вида напряженного состояния на процесс теплопроводности.

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

0 ,ii -µ,t - ui,it = N -U, (29) где µ = C - 3AtT0D3 ; = AtT0 D1 - 2D2 ; N = T0 -At ijHij,t +BtS,t .

( ) ( ) Очевидно, что в уравнении теплопроводности (29) фигурируют дополнительные слагаемые, учитывающие связность задачи и влияние чувствительности вида напряженного состояния на теплофизические характеристики материалов.

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

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

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

Рассматривается упругое равновесие тонкой пластины толщиной h из разносопротивляющегося материала под действием поперечной равномерно распределенной нагрузки с интенсивностью q. Положение любой точки определено в ортогональной системе координат xk k =1, 2, 3. При этом плоскость, образованную осями x1 и x( ) совмещается со срединной поверхностью пластинки в недеформированном состоянии, а ось x3 ориентируется перпендикулярно этой поверхности в направлении прогиба.

Рассматриваются достаточно тонкие пластины, такие, чтобы применение гипотез Кирхгофа-Лява не вызывало возражений. Деформированное состояние пластин определяется компонентами перемещений точек срединной поверхности u1, u2, w.

Компоненты тензора деформаций выражаются через параметры деформации ij и кривизны ij срединной поверхности i, j = 1, 2 :

( ) eij = ij + x3ij. (30) Рассматриваются геометрически линейные задачи изгиба тонких пластин в декартовой и цилиндрической системе координат. Применительно к конкретным условиям задачи деформации и кривизны срединной поверхности запишем следующим образом, для декартовой системы координат:

ij = ui, j + u 2, ij = -w,ij. (31) ( ) j,i Для цилиндрической системы координат выражения (30) примут вид:

er = r + zr, e = + z. (32) где r, — кривизны в радиальном и окружном направлениях; r, — деформации срединной поверхности.

В цилиндрической системе координат деформации и кривизны срединной поверхности записываются следующим образом:

r = u,r, = u r, r = -w,rr, = - w,r r, (33) где w,r — угол поворота нормали к срединной поверхности пластины; — радиальu, w ные перемещения и прогибы срединной плоскости.

Усилия в срединной поверхности пластин определяются в виде:

h / 2 h / Nij = ijd x3, Mij = ij x3 d x3, (34) -h / 2 -h / где Nij – продольные усилия; Mij – изгибающие и крутящие моменты.

Для этих усилий имеем уравнения статики (с учетом продольных усилий в срединной плоскости):

N11,1 + N12,2 = 0, N12,1 + N22,2 = 0, (35) M11,11 + 2M12,12 + M22,22 = -q - N11w,11 -2N12w,12 -N22w,.

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

h D1u1,11 +0,5 D1 - D2 u2,12 +0,5 D1 + D2 u1,22 = g1; (36) ( ) ( ) h D1u2,22 +0,5 D1 - D2 u1,12 +0,5 D1 + D2 u2,11 = g2; (37) ( ) ( ) w,1111 +2w,1122 +w,2222 D1h3 12 = q + g3; (38) ( ) 0 ,33 -( - 3 AtT0D3 ,t - AtT0 D1 - 2 D2 u1,1t +u2,2t = g4. (39) C ) ( )( ) где gi = gi,,S0, ( ) нелинейные функции (i = 1,2,3), в развернутом виде представлены в диссертации.

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

Начальные и граничные условия для температуры, заключаются в том, что в начальный момент времени пластина имеет температуру T0, и нагревается с двух сторон с некоторым перепадом температуры T=(T1-T2), на верхней и на нижней поверхностях пластин поддерживается постоянная температура T1 и T2, соответственно.

Разрешающие уравнения (36)—(39) вместе с начальными и граничными условиями образуют замкнутую систему, описывающую связанную задачу термоупругости для тонких прямоугольных пластин, выполненных из изотропных разносопротивляющихся материалов.

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

Усилия в срединной поверхности пластин определяются аналогично (34).

Условия равновесия принимаются с учетом влияния продольных усилий:

N + Nr - N r = 0, Mr,r + Mr - M r = Qr (40) ( ) ( ) r,r где Nr, N — продольные усилия в радиальном и окружном направлениях; Mr, M — изгибающие моменты.

Величина поперечной силы Qr в случае равномерно распределенной нагрузки определяется следующим образом:

Qr = -N w,r - qr 2 (41) r где q — интенсивность поперечной равномерно распределенной нагрузки.

Подставляя выражения (34), записанные для цилиндрической системы координат в (40), а также учитывая (18), (32) и (33), можно получить систему уравнений равновесия пластины формы в перемещениях (2 уравнения). Полученные уравнения равновесия и уравнение притока тепла (29) образуют полную систему дифференциальных уравнений описывающих напряженное состояние круглых пластин в условиях термомеханического нагружения. После выполнения преобразований, разрешающая система дифференциальных уравнений в частных производных приобретает вид:

* hrD1 rr +u,r r - u r2 g1 ; (42) (u, )= * (43) (w, + w,rr r - w,r r2)D h3 12 = qr 2 + g2;

rrr 0 0 * ,zz -( - 3 AtT0D3 ,t -AtT0 D1 - 2 D2 u,rt = g3, (44) C ) ( ) * где gi = gi,,S0, ( ) нелинейные функции (i = 1,2,3), в развернутом виде представлены в диссертации.

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

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

Ограничимся рассмотрением пластин в линейной геометрической постановке.

С целью сравнения, получаемых решений проводилось сравнение результатов расчета по рассмотренной выше модели (обозначение на рисунках — РС) с результатами расчета на основе физических соотношений Дюгамеля-Неймана, а именно не связанная термоупругость без учёта разносопротивляемости с усреднёнными деформационными характеристиками при растяжении и сжатии (КН), с деформационными характеристиками принятыми как при растяжении (П+) и при сжатии (П-), а также с результатами расчетов, полученными по моделям Толоконникова-Матченко и С.А.

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

Разбиение пластины, для использования метода конечных разностей, проводилось с постоянным шагом по трём направлениям с постоянным шагом hi (i=1,2,3) на N-1 участков. Дискретизация по времени также проводилась с постоянным шагом, равным h4. Разностные аналоги разрешающих дифференциальных уравнений получены путем замены производных конечными разностями: первого порядка «правыми», «левыми», «центральными»; «центральными» второго порядка и бигармоническими четвёртого порядка. Для аппроксимации производных по времени применяется неявная конечно-разностная схема. Значения функций в законтурных и граничных точках определяются из граничных и начальных условий. Численные исследования показали устойчивость принятой разностной схемы для различных сеток. Исследование сходимости показали достаточную точность принятого шага разбиения сетки.

Уравнения (36)(39) для прямоугольных и (42)(44) для круглых пластин представлены в форме (все нелинейные члены выписаны в правых частях), которая удобна для применения метода «упругих решений». В процессе, прослеживается влияние температуры на механические характеристики материалов и напряженного состояния на распределение температуры по толщине и в плане пластинки. Решение системы дифференциальных уравнений на каждой итерации по методу «упругих решений» находится частное решение системы разрешающих уравнений. Алгоритм решения рассмотренной задачи реализован в интерактивной среде для инженерных расчетов MATLAB (MathWorks Inc., США).

Для демонстрации возможностей предлагаемой теории решалась задача со следующими исходными данными: толщина пластины h=200 мм, размеры квадратной пластины a=b=2000 мм, размеры круглой пластины r=1000 мм, материал графит АРВ, жестко защемлена (шарнирно опёрта) по контуру; пластина нагружалась равномерно распределенной нагрузкой q=0…13 кПа; также осуществлялся нагрев поверхности пластины. Постоянная температура на нижней поверхности пластины T2=50 C, на верхней — T1=10 C. Начальная температура пластины T0=20 C. Механические характеристики материала: модуль упругости E+=3750 МПа; E-=6130 МПа; коэффициент Пуассона +=0,2, -=0,35; плотность =1700 кг/м; коэффициенты линейного теплового расширения t+=4.0·10-6 К; t-/ t+=1,5; коэффициент теплопроводности =150 Вт/(м·К); теплоемкость материала С=500 Дж/(кг·К).

На рис. 1—4 представлены некоторые полученные результаты расчёта круглой шарнирно опёртой и жёстко защемлённой пластин.

20% -----РС -РС Г+ Г+ Г20% -6 ГКН КН ---10 -5 0 0 100 200 300 400 5Напряжения, МПа r Радиус пластины, мм Рис. 1. Распределение напряжений r Рис. 2. Распределение напряжений r в центре шарнирно опёртой пластины в шарнирно опёртой пластине сверху Прогиб w жёстко защемлённой пластины 4 --РС -Г+ Г-КН --РС 10% Г+ --8 Г- 20% КН --0 100 200 300 400 50 100 200 300 400 5Радиус пластины, мм Радиус пластины, мм Рис. 3. Распределение напряжений r Рис. 4. Распределение прогиба w в жёстко защемлённой пластине снизу в жёстко защемлённой пластине Аналогичные и другие результаты, полученные для прямоугольных (квадратных) и круглых жёстко защемлённых и шарнирно опёртых пластин, представлены в диссертации.

В заключении приведены основные результаты и выводы по работе.

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

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ 1. Впервые получено решение связанной термомеханической задачи изгиба тонких пластин из изотропных разносопротивляющихся материалов. Получено новое решение научно-технической задачи механики деформированного твёрдого тела, учитывающее влияния текущего состояния НДС на величины температурных характеристик материалов и влияние температуры на напряжённое состояние. Подтверждено наличие известных фактов и обнаружены новые эффекты деформирования. Показано, что совместный учет влияния НДС на величины коэффициентов температурного расширения и эффекта разносопротивляемости, позволяют получить уточнение расчетных результатов по сравнению с существующими моделями.

2. В рамках нормированных пространств напряжений, предложенных в работах Л.А. Толоконникова, Н.М. Матченко и А.А. Трещева, проведены аналогичные исслеr Толщина, мм Напряжения, МПа r Прогиб w, мм Напряжения, МПа дования для модификации термоупругого потенциала Гиббса. Получены разрешающие уравнения термоупругости для тонких пластин из изотропных разносопротивляющихся материалов в декартовой и цилиндрической системах координат при малых перемещениях. Разработана математическая модель решения задачи.

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

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

ОСНОВНЫЕ ПОЛОЖЕНИЯ ДИССЕРТАЦИИ ОПУБЛИКОВАНЫ В РАБОТАХ (публикации в изданиях, включенных в перечень ВАК российских рецензируемых научных журналов, в которых должны быть опубликованы основные научные результаты диссертаций на соискание ученой степени кандидата наук):

1. Чигинский Д. С. Решение задачи об изгибе тонкой прямоугольной пластины из разносопротивляющихся материалов в условиях термомеханического нагружения / Д. С. Чигинский, В. Г. Теличко, А. А. Петров // Известия ТулГУ. Технические науки. Вып. 1. В 2 ч. Ч. 2. — Тула: Изд-во ТулГУ, 2009. — с. 114—120.

2. Чигинский Д. С. Связанная задача термомеханического изгиба тонких прямоугольных пластин из изотропных разносопротивляющихся материалов / Д. С. Чигинский, А. А. Трещёв, В. Г. Теличко // Известия ТулГУ. Технические науки. Вып. 2. Проблемы специального машиностроения. — Тула: Изд-во ТулГУ, 2011. — с. 494—502.

3. Чигинский Д. С. Анализ определяющих соотношений для нелинейных изотропных разносопротивляющихся материалов в задачах термоупругости / Д. С. Чигинский, А. А. Трещёв, В. Г. Теличко, // Известия ТулГУ. Технические науки. Вып. 2. Проблемы специального машиностроения. — Тула: Изд-во ТулГУ, 2011. — с. 547—555.

(публикации в остальных изданиях):

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

Изд-во ТулГУ, 2007. – с. 288–294.

5. Чигинский, Д. С. Расчёт напряжённо-деформированного состояния круглых пластин из материалов с усложнёнными свойствами в условиях термомеханического нагружения / Д. С. Чигинский, А. А. Петров, В. Г. Теличко // Материалы МНК «Современные проблемы математики, механики, информатики». – Тула: ТулГУ, 2007. – с.

196–197.

6. Чигинский, Д. С. Решение связанной задачи термоупругости для тонких прямоугольных пластин из графита / Д. С. Чигинский, А. А. Петров, В. Г. Теличко // Материалы международной научной конференции «Современные проблемы математики, механики, информатики». – Тула: Изд-во ТулГУ, 2007. – с. 198–199.

7. Чигинский, Д. С. К решению задачи об изгибе тонкой круглой пластины из разносопротивляющихся материалов в условиях термомеханического нагружения / Д. С. Чигинский, А. А. Петров, В. Г. Теличко // Сб. статей VII Международной НТК, посвященной 50-летию ПГУАС. – Пенза: ПДЗ, 2008. – с. 153–156.

8. Чигинский, Д. С. Характерные результаты расчета НДС тонкой прямоугольной пластины из разносопротивляющегося материала в условиях термомеханического нагружения / Д. С. Чигинский, В. Г. Теличко // Надёжность и долговечность строительных материалов, конструкций и оснований фундаментов. Материалы V Международной конференции, г. Волгоград, 23–24 апреля 2009 г. В 3-х ч. Ч. III. – Волгоград:

ВолгГАСУ, 2009. – с. 242–247.

9. Чигинский, Д. С. Решение связанной задачи о поперечном изгибе тонкой прямоугольной пластины из графита в условиях термомеханического нагружения / Д. С.

Чигинский, А. А. Трещёв, В. Г. Теличко // Актуальные проблемы прикладной математики, информатики и механики. – Воронеж: ВГУ, 2009. – Ч. 2. – с. 204–209.

10. Чигинский, Д. С. Термомеханический изгиб прямоугольных пластин из разносопротивляющихся материалов с учетом связанности задачи / Д. С. Чигинский, А. А.

Трещёв, В. Г. Теличко // Исследования и инновационные разработки РААСН. – Москва-Иваново: РААСН-ИГАСУ, 2010. – Т. 2. – с. 216–221.

11. Чигинский, Д. С. Связанная термомеханическая задача изгиба тонких прямоугольных пластин из изотропных материалов с усложнёнными свойствами / Д. С. Чигинский // Молодёжь и наука: сборник материалов VII Всероссийской НТК студентов, аспирантов и молодых ученых, посвященной 50-летию первого полета человека в космос / отв. ред. О. А. Краев. – Красноярск: Сиб. федер. ун-т., 2011. – с. 37–40.

12. Чигинский, Д. С. Связанные термомеханические задачи изгиба тонких прямоугольных пластин из изотропных разносопротивляющихся материалов / Д. С. Чигинский, В. Г. Теличко, А. А. Трещёв // Материалы международной научной конференции «Современные проблемы математики, механики, информатики». – Тула: Изд-во ТулГУ, 2011. – с. 230–235.

13. Чигинский, Д. С. Связанные задачи термомеханического изгиба тонких прямоугольных пластин из изотропных разносопротивляющихся материалов / Д. С. Чигинский // 7-я Международная конференция по проблемам горной промышленности, строительства и энергетики. Материалы конференции. Т. 2. – Тула: Изд-во ТулГУ, 2011. – с. 103–111.







© 2011 www.dissers.ru - «Бесплатная электронная библиотека»

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