Научный журнал
Международный журнал прикладных и фундаментальных исследований

ISSN 1996-3955
ИФ РИНЦ = 0,580

ПРИМЕНЕНИЕ ПАКЕТОВ СИМВОЛЬНОЙ МАТЕМАТИКИ К ИССЛЕДОВАНИЮ КОНЦЕНТРАЦИИ НАПРЯЖЕНИЙ НА КОНТУРЕ ОТВЕРСТИЙ

Щукина Н.А. 1
1 ФГБОУ ВПО «Московский государственный университет экономики статистики и информатики (МЭСИ)»
Статья посвящена применению метода эффектов второго порядка для решения плоских задач нелинейной теории упругости. Данный метод использует разложение в степенные ряды по малому параметру объектов, описывающих напряженно-деформированное состояние. В рамках построенной приближенной математической модели найдено аналитическое решение задачи о концентрации напряжений на контуре квадратного отверстия, свободного от нагрузок, при равномерном растяжении на бесконечности. Данный алгоритм реализован в пакете символьной математики Maple. Показано, что величина коэффициента концентрации напряжений зависит от внешней нагрузки. При этом при определенных значениях малого параметра в области концентраторов напряжений наблюдается раздвоение максимума на два симметричных.
нелинейная теория упругости
эффекты второго порядка
приближенная математическая модель
коэффициент концентрации напряжений
система автоматизированного вычисления
1. Жуков, Б.А. Эффекты третьего порядка в исследовании концентрации напряжений около отверстий / Б.А. Жуков, Н.А. Щукина // Известия ВолгГТУ. – 2010. – Т.1. – №3. – С. 113–118.
2. Жуков, Б.А. Модель эффектов третьего порядка в статических задачах расчетов резинотехнических изделий / Б.А. Жуков, Н.А. Щукина // Изв. вузов. Сев.-Кавк. регион. Серия Естественные науки. – 2010. – № 3. – С. 24–27.
3. Жуков Б.А. Один вариант метода Синьорини при плоской деформации несжимаемого материала / Б.А. Жуков // Изв. РАН. МТТ. – 2001. – № 4. – С. 59–67.
4. Мусхелишвили Н. И. Некоторые основные задачи математической теории упругости. 5-е изд. – М.: Наука, 1966. – 707 с.

В настоящее время поиск аналитических решений задач нелинейной теории упругости является достаточно сложной задачей. Одной из сложностей при построении математических моделей таких задач является учет несжимаемости материала. Условие несжимаемости несет дополнительную информацию о геометрии деформирования. В работе [1] предложена приближенная математическая модель нелинейной теории упругости для плоской деформации однородного изотропного материала, в рамках которой условие несжимаемости выполняется автоматически. В качестве метода построения приближенной модели плоской деформации используется метод возмущений, использующий разложения по степеням малого параметра объекты, описывающие напряженно-деформированное состояние. Однако из-за громоздких вычислений разложение выше второго порядка практически не используется. Появление современных пакетов символьной математики позволяет написать программы, облегчающие манипулирование с громоздкими выражениями, описывающими эффекты второго порядка при произвольном напряженно-деформированном состоянии. Эти эффекты выделяются при ограничении разложений для радиус-вектора частиц в текущей конфигурации shcuki1.wmf и функции гидростатического давления p членами до второго порядка по малому параметру shcuki2.wmf.

Математическая постановка задачи

В рамках рассматриваемой модели можно ограничиться выражением для потенциала энергии деформации в форме, предложенной в [2]:

shcuki3.wmf,

где m, µ1 и µ2 – константы, причем m – модуль сдвига линейной теории, shcuki4.wmf – главные инварианты меры деформации Коши G.

Следуя работе [3], разложение радиус-вектора точек в плоскости, ортогональной оси OZ, в текущей конфигурации представляется в виде

shcuki5.wmf. (1)

Здесь r – радиус-вектор точек в отсчетной конфигурации, shcuki6.wmf – оператор Гамильтона в базисе исходной конфигурации, shcuki7.wmf – симплектический оператор. Дифференциальные операторы действуют только на первый множитель справа, не являющийся оператором. Точка означает скалярное произведение.

Используя разложение (1), получим представления в рамках эффектов второго порядка тензора напряжений Коши

shcuki8.wmf, (2)

где shcuki9.wmf, shcuki10.wmf,

shcuki11.wmf, shcuki12.wmf, shcuki13.wmf,

shcuki14.wmf,

shcuki15.wmf – функция, связанная с гидростатическим давлением, знак «Т» – транспонирование.

Уравнения равновесия приводятся к системам дифференциальных уравнений для эффектов первого и второго порядков соответственно:

shcuki16.wmf,

shcuki17.wmf, (3)

Условие интегрируемости системы уравнений (3) для эффектов первого, второго и третьего порядков соответственно приводит к дифференциальным уравнениям:

shcuki18.wmf, shcuki19.wmf.

Обозначим shcuki20.wmf – разложение вектора плотности внешних сил, приложенного к деформированной боковой поверхности и рассчитанного на единицу площади этой поверхности. Тогда силовые граничные условия в напряжениях запишутся в виде

shcuki21.wmf,

shcuki22.wmf, (5)

где shcuki23.wmf, shcuki24.wmf.

Если плотность внешних сил shcuki25.wmf задана в отсчетной конфигурации, а нагружение «мертвое», то shcuki26.wmf, shcuki27.wmf.

Выражения (1) – (5) представляют постановку граничной задачи нелинейной теории упругости. Так как уравнения (4) для эффектов первого и второго порядков приводят к бигармоническим уравнениям относительно искомых функций shcuki28.wmf и shcuki29.wmf, то поиск аналитического решения краевых задач для эффектов первого и второго порядка можно проводить единым образом в автоматическом режиме.

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

shcuki30.wmf,

где shcuki31.wmf.

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

Перейдем к комплексным переменным shcuki32.wmf, shcuki33.wmf и введем комплексные потенциалы по формуле Гурса:

shcuki34.wmf,

shcuki35.wmf.

Для данных потенциалов уравнения (4) выполняются тождественно, а решением системы уравнений равновесия (3) являются функции:

shcuki36.wmf,

shcuki37.wmf,

где

shcuki38.wmf, shcuki39.wmf, shcuki40.wmf, shcuki41.wmf.

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

shcuki42.wmf и shcuki43.wmf

для эффектов первого и второго порядков соответственно.

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

Задача о концентрации напряжений около квадратного отверстия при равномерном растяжении на бесконечности

Без потери общности будем считать радиус отверстия равным единице. Контур отверстия свободен от напряжений. На бесконечности приложена равномерная радиальная нагрузка интенсивности p, Н/м2. Универсальным силовым параметром, имеющим ту же размерность, является модуль сдвига линейной теории упругости shcuki44.wmf. Поэтому естественно возникает безразмерный малый параметр вида shcuki45.wmf.

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

shcuki50.wmf,

где shcuki51.wmf – тангенциальная компонента «плоской» части тензора истинных напряжений Коши shcuki52.wmf. Представление этой компоненты в рамках приближенной теории имеет вид shcuki53.wmf, откуда, вспоминая, что shcuki54.wmf, получим выражение для коэффициента концентрации в рамках приближенной теории:

shcuki55.wmf.

Конформное отображение внешности квадрата на внешность круга shcuki56.wmf задается формулой shcuki57.wmf. Вершинам квадратного контура соответствуют значения полярного угла shcuki58.wmf, shcuki59.wmf и shcuki60.wmf. На бесконечности приложена равномерная нагрузка интенсивности shcuki61.wmf, где shcuki62.wmf – малый параметр. Тогда на бесконечности вектор внешних сил shcuki63.wmf на площадках с нормальным вектором shcuki64.wmf имеет разложение shcuki65.wmf, shcuki66.wmf, а на площадке с нормальным вектором shcuki67.wmf разложение – shcuki68.wmf, shcuki69.wmf.

Для квадратичного приближения разложение коэффициента концентрации напряжений по малому параметру shcuki70.wmf в вершинах квадратного отверстия принимает вид shcuki71.wmf. Таким образом, наблюдаем уменьшение значения тангенциальных напряжений с увеличением величины интенсивности приложенной внешней нагрузки shcuki72.wmf. Графики распределения относительных тангенциальных напряжений на внутреннем контуре отверстия, свободном от нагрузок, вычисленных в рамках эффектов первого и второго порядков при различных значениях параметра shcuki73.wmf, изображены на рис. 1 и 2.

schuk1.tif

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

schuk2.tif

Рис. 2. Распределение относительных тангенциальных напряжений на контуре отверстия (развертка) при shcuki75.wmf. Тонкая линия соответствует решению в рамках линейной теории, толстая линия соответствует квадратичному решению

Как видно из графиков, в квадратичном приближении при различных значениях параметра shcuki76.wmf максимум величины коэффициента концентрации напряжений k достигается при различных значениях j, при этом его величина, по сравнению с классическим решением линейной теории, уменьшается. Заметим, что при учете внешней нагрузки p в разложении коэффициента концентрации напряжений наблюдаем эффект раздвоения одного максимума в вершине угловой точки на два симметричных максимума в окрестности данной вершины. При этом в самой вершине достигается минимум.

Исследовав функцию shcuki77.wmf на экстремум, получаем, что значениям

shcuki78.wmf

соответствует один экстремум: максимум достигается в вершинах квадратного контура. При shcuki79.wmf каждому значению shcuki80.wmf соответствует три точки экстремума: в вершинах квадрата достигается минимум, а двум другим корням производной shcuki81.wmf соответствуют максимумы.

Максимальные значения коэффициента концентрации напряжений shcuki82.wmf, вычисленные с точностью 10–6, и соответствующие им значения полярного угла shcuki83.wmf в зависимости от величины малого параметра h записаны в таблице.

Экстремальные значения коэффициента концентрации shcuki84.wmf и соответствующие им значения полярного угла shcuki85.wmf

shcuki86.wmf

shcuki87.wmf

shcuki88.wmf

0.1

8.4929885

0

0.2

8.0174237

0

0.3

7.5418589

0

0.364

7.2374980

0.0031250

0.4

7.0911554

0.0236328

0.5

6.8714057

0.0437500

0.6

6.8041820

0.0500000

0.7

6.8603106

0.0633789

0.8

6.9618454

0.0696289

0.9

7.1029659

0.0750000

Выводы

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

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


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

Щукина Н.А. ПРИМЕНЕНИЕ ПАКЕТОВ СИМВОЛЬНОЙ МАТЕМАТИКИ К ИССЛЕДОВАНИЮ КОНЦЕНТРАЦИИ НАПРЯЖЕНИЙ НА КОНТУРЕ ОТВЕРСТИЙ // Международный журнал прикладных и фундаментальных исследований. – 2015. – № 6-2. – С. 228-233;
URL: http://www.applied-research.ru/ru/article/view?id=6878 (дата обращения: 19.04.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074