Сетевое издание

Современные проблемы науки и образования

ISSN 2070-7428

"Перечень" ВАК

ИФ РИНЦ = 1,006

• Авторы
• Резюме
• Файлы
• Ключевые слова
• Литература

Шиверский А.В. 1 Горохов А.Ю. 1 Максимов П.В. 1 Блинов А.В. 1

---

1 ФГБОУ ВПО «Пермский национальный исследовательский политехнический университет»

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

Статья в формате PDF

436 KB

анизотропия

многооператорная задача

вязкоупругость

метод конечных элементов

1. Думанский А.М. Обращение определяющих уравнений наследственно-упругой ортотропной среды при плоском напряженном состоянии // Механика композиционных материалов и конструкций. – 2000. – Т. 6, № 4. – С. 482–493.

2. Лехницкий С.Г. Теория упругости анизотропного тела. – М.: Наука, 1977. – 416 с.

3. Малый В.И., Труфанов Н.А. Метод квазиконстантных операторов в теории вязкоупругости кусочно-однородных материалов // Деформирование и разрушение структурно неоднородных материалов. – Свердловск: УрО РАН, 1989. – С. 78–85.

4. Методы прикладной вязкоупругости / Адамов А.А. [и др.]. – Екатеринбург: УрО РАН, 2003. – 411 с.

5. Плуме Э.З. Сравнительный анализ ползучести однонаправленных композитов, армированных волокнами различного типа // Механика композиционных материалов. – 1985. – № 3. – С. 431–436.

6. Работнов Ю.Н. Элементы наследственной механики твердых тел. – М.: Наука. – 1977. – С. 37–43.

7. Труфанов Н.А., Черепанова Л.М. Применение частичных аппроксимаций для повышения точности метода квазиконстантных операторов // Вестник ПГТУ. Прикладная математика и механика. – 2000. – № 1. – С 70–76.

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

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

Рассматриваются следующие варианты нагружения: сдвиг усилиями, параллельными главным направлениям (задача 1); касательные усилия, равномерно распределенные по кромке отверстия (задача 2); внутренне давление по кромке отверстия (задача 3) (рис. 1).

Рис. 1. Схема нагружения пластинки: а) сдвиг усилиями, параллельными главным направлениям (задача 1); б) касательные усилия, равномерно распределенные по кромке отверстия (задача 2); в) внутренне давление по кромке отверстия (задача 3)

Постановка задачи в рамках линейной теории вязкоупругости:

· уравнения равновесия: , ;

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

· физические соотношения в условиях плоского напряженного состояния (ПНС):

(1)

где – время; , , , – интегральные операторы вида

(2)

Граничные условия при различных видах внешнего нагружения будут иметь вид:

· сдвиг усилиями, параллельными главным направлениям (задача 1)

, .

· касательные усилия, равномерно распределенные по кромке отверстия (задача 2)

,

· внутренне давление по кромке отверстия (задача 3)

, , ,

.

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

В упругой постановке рассмотренные задачи имеют аналитические решения, приведенные в [2]. Построены конечно-элементные модели задач и выполнена проверка их точности и сходимости. Сравнение численного решения с аналитическим производилось по окружному напряжению возле кромки отверстия. Численные расчеты проведены в пакете ANSYS с использованием следующих упругих характеристик материала пластинки: Па, Па, , Па. Внешнее усилие принято постоянными , где – функция Хевисайда. Результаты приведены в табл.1.

Таблица 1

Оценка численных моделей задач

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

В [5] приведены данные об аппроксимации кривых ползучести однонаправленного органопластика. Ядра функций ползучести приняты в виде суммы экспонент. Физические соотношения для органопластика [5] содержат четыре независимых оператора вида (2). Обращение определяющих уравнений вязкоупругой среды произведено в предположении, что все операторы являются квазиконстантными. В этом случае для любого момента времени справедливо равенство , где – тензор постоянных упругости, а – тензор постоянных податливости. Значения показатели квазиконстантности операторов приведены в табл. 2.

Таблица 2

Показатели квазиконстантности органопластика и стеклотекстолита

Характеристики ортотропной пластинки получены в виде пакета из двух ортогонально армированных однонаправленным композитом слоев: , где , – тензора постоянных упругости ортотропной пластинки и однонаправленного органопластика; – тензор постоянных упругости однонаправленного органопластика при повороте на .

Для стеклотекстолита из работы [1] определяющие уравнения содержат три независимых оператора: , и , где – дробно-экспоненциальная функция Работнова. Дальнейшие расчеты проводились в предположении, что операторы являются квазиконстантными. Значения показателей квазиконстантности определены по формуле из работы [4] и приведены в табл. 2:

,

где , , ,

, .

Вычисление значений функций и произведено в пакете MatLab. При больших значениях аргумента использованы асимптотические разложения и связь функции с функцией Миттаг—Лаффлера [6].

Поставленные краевые задачи 1, 2, 3 для материала органопластика и стеклотекстолита решены при постоянной внешней нагрузке приближенным методом квазиконстантных операторов. Произведена оценка погрешности решения задач указанным методом. Оценка погрешности решения для материала, описываемого несколькими независимыми операторами, производится по формуле [4]:

, где – число независимых операторов.

Для каждого из слагаемых погрешности выполняется оценка , где – показатели квазиконстантности независимых операторов.

Исходя из приведенных в таблице 2 показателей квазиконстантности имеем следующие оценки погрешности: органопластик – ; стеклотекстолит – .

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

На рисунке 2 в полярной системе приведены графики изменения величин для трех различных моментов времени. Здесь – окружное напряжение возле кромки отверстия, а коэффициент «0,5» имитирует отверстие радиуса 0,5, аналогично графикам, предложенным в работе [2].

Рис. 2. Эволюция окружных напряжений у кромки отверстия.

Задача 1: а) органопластик, б) стеклотекстолит. Задача 2: в) органопластик, г) стеклотекстолит. Задача 3: д) органопластик, е) стеклотекстолит

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

Рис. 3. Изменение коэффициента концентрации напряжений во времени возле кромки отверстия: а) органопластик; б) стеклотекстолит

Количественное изменение коэффициентов концентрации, приведенных на рис.3, представлено в табл. 3.

Таблица 3

Изменение коэффициента концентрации напряжений во времени

Библиографическая ссылка