- •Математический факультет Кафедра вычислительной математики и программирования
- •Содержание
- •2 Разностные методы решения задач математической физики….…..20
- •Разностные схемы для линейных дифференциальных уравнений
- •1 Математические модели физических систем
- •1.1 Физические системы и методы их исследования
- •Основные определения и понятия
- •Системы деформируемых твёрдых тел
- •1.1.3 Подходы к исследованию систем деформируемых твёрдых тел
- •1.1.4 Средства исследования состояния систем деформируемых твёрдых тел
- •Разностные схемы для основных уравнений математической физики
- •1.2.1 Основные уравнения математической физики и их классификация
- •1.2.2 Разностная аппроксимация основных производных. По определению производная функции одной переменной записывается в виде
- •Это значит, что аппроксимация
- •Общий принцип метода сеток
- •Сходимость, аппроксимация и устойчивость разностных схем
- •2 Разностные методы решения задач математической физики
- •2.1 Разностные схемы для линейных дифференциальных уравнений
- •2.1.1 Краевые задачи для линейных дифференциальных уравнений
- •2.1.2 О задаче Дирихле
- •2.1.3 Разностные аппроксимации уравнений эллиптического типа
- •2.1.4 Разностная аппроксимация граничных условий
- •2.2 Построение разностной аппроксимации задачи Дирихле для
- •2.2.1 Разностная схема задачи Дирихле для уравнения Пуассона.
- •2.2.2 Устойчивость и сходимость разностной задачи Дирихле
- •2.3. Разрешимость разностных уравнений и способы их решения. Правило Рунге
- •2.4. Разностные схемы для линейных дифференциальных
- •2.5. Разностные схемы для линейных дифференциальных уравнений
- •2.5.1. Общая постановка задачи.
- •3. Метод конечных элементов и суперэлементов
- •3.1. Физические предпосылки методов
- •3.2 Деформации твёрдых тел
- •3. 2.1 Линейно-упругие деформации твёрдых тел
- •3.2.2. Нелинейно-упругие деформации твёрдых тел
- •3.3. Прикладные вопросы метода конечных элементов
- •3.3.1. Основная концепция метода конечных элементов
- •3.3.2 Дискретизация и оптимальная нумерация узлов области
- •3.3.3. Метод конечных элементов в теории упругости.
- •3..4 Построение конечно - элементных соотношений
- •3.5. Аналитический алгоритм метода конечных элементов для исследования
- •Свойства матрицы жёсткости
- •3.7 Аналитический алгоритм метода суперэлементов
- •Методы исследования нелинейных математических
- •4.1 Общие предпосылки методов численного исследования нелинейных
- •4.2 Итерационные численные методы
- •4.3 Безитерационные методы исследования нелинейных математических моделей деформируемых твёрдых тел и их систем
- •4.4 Сравнительный анализ эффективности методов исследования
- •5 Основы методологии компьютерного объектно-
- •5.1 Основные понятия и определения
- •5.2 Построение виртуальной физической модели системы
- •5.3 Исходные данные при компьютерном моделировании
- •5.4 Деформации грунтовых оснований при устройстве отдельных
- •6 Технология компьютерного объектно-
- •6.1 Технологические этапы компьютерного объектно-ориентированного
- •6.2 Интерфейс визуального ввода данных
- •7 Программное обеспечение компьютерного
- •7.1 Структура программного обеспечения объектно-ориентированного моделирования системы деформируемых твёрдых тел
- •7.2 Технология объектно-ориентированного программирования
- •Объектно-ориентированного моделирования физической системы
- •7.3 Программный комплекс «Энергия - 2д » моделирования
- •7.4 Программный комплекс «Энергия - ос » моделирования
- •7.5 Программный комплекс «Энергия - 3д » моделирования
- •8 Верификация технологии и программного
- •8.1 Методологические аспекты верификации
- •8.2 Верификация технологии и программного обеспечения
- •8.3 Верификация технологии и программного обеспечения
- •8.4 Верификация технологии и программного обеспечения
- •Численные методы математической физики
- •1. Метод сеток для задачи дирихле
- •1.1. Основы метода сеток
- •1.3.Варианты задания
- •Лабораторная работа №2
- •Дополнение к лабораторной работе №2
- •Лабораторная работа №3 Метод сеток для уравнения параболического типа.
- •Лабораторная работа №4 Метод сеток для уравнения гиперболического типа
4.3 Безитерационные методы исследования нелинейных математических моделей деформируемых твёрдых тел и их систем
Рассмотренные и другие методы решения нелинейных краевых задач теории упругости являются итерационными и основаны на идее решения нелинейных задач методами линейной теории упругости. Реализация такого подхода сводится к построению алгоритма изменения модуля упругости элементов деформируемого тела или к изменению вектора сил. Эти методы имеют слабую сходимость или вообще не приводят к искомому решению (метод переменной жёсткости).
В настоящей работе предлагаются оригинальные двухэтапные методы решения краевых задач нелинейной теории упругости [2, 11, 12]. В этих методах рассматривается однородное нелинейно-деформируемое твёрдое тело объёма V c границей Г и системой граничных условий, соответствующих первой или третьей краевым задачам теории упругости. Уравнение закона деформирования твёрдого тела в общем случае представим виде:
, (4.6)
где f (**) – функция модуля упругости при нелинейном деформировании;
**- совокупность параметров, определяющих значение модуля упругости.
Начальное значение модуля упругости Е = Е0 .
Твёрдому телу объёма V c границей Г и законом деформирования
(4.7)
поставим в соответствие геометрически тождественное гипотетическое линейно-упругое тело с законом деформирования
. (4.8)
Модуль упругости подлежит определению и должен быть таким, чтобы при тождественных граничных условиях для обоих тел их смещения совпадали. В соответствии с принципом возможных перемещений для всякой сплошной среды
П dV – W) = 0, (4.9)
где П - потенциал деформации; W - работа внешних сил.
Для единичного элемента
(4.10)
где ε - средняя деформация; k – коэффициент объемного сжатия.
Метод энергетической линеаризации
В разработанном методе энергетической линеаризации краевой задаче нелинейной теории упругости телу объёма V c границей Г и законом деформирования ставится в соответствие геометрически тождественное гипотетическое линейно-упругое тело с законом деформирования
. (4.11)
Модуль деформации Ег подлежит определению и должен быть таким, чтобы при тождественных граничных условиях для обоих тел их смещения совпадали.
В силу поставленного условия для модуля деформации Ег и тождественности граничных условий для рассматриваемых твёрдых тел можно утверждать, что работы внешних сил на смещениях для исследуемого нелинейно – деформируемого и гипотетического линейно-упругого тел будут равны, т.е. будем иметь Wn = Wr, тогда
,
Пn - энергия деформации нелинейно-упруго
Пr - энергия деформации гипотетического линейно - упругого тела; индексы “n” и “r”– признаки нелинейно-упругого и гипотетического линейно - упругого тел.
Полученное выражение представим в следующей форме
, (4.12)
где F( ) = .
В соответствии с основной леммой вариационного исчисления из (4.12) следует:
F( ) - / 2 = 0. (4.13)
При законе деформирования в виде степенной функции из (4.13) получим:
(4.14)
Подставим (4.14) в уравнение состояния (11) и, учитывая закон Гука
,
решим его относительно i, получим:
Если закон деформирования имеет вид
,
то определяется из уравнения: E0 - 2B/(1+m) = , в этом случае
E r = / = E0 - 2B / (1+m) .
Для полученных значений Еr решается линейная задача. Полученное решение будет являться и решением исходной нелинейной задачи.