- •Математический факультет Кафедра вычислительной математики и программирования
- •Содержание
- •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 Метод сеток для уравнения гиперболического типа
3.3.3. Метод конечных элементов в теории упругости.
При применении метода конечных элементов для решения задач теории упругости сплошное тело рассматривают условно состоящим из конечного числа отдельных элементов, вплотную прилегающих друг к другу и шарнирно скрепленных между собой в вершинах этих элементов. Форма и размеры тела остаются неизменными. В случае полупространства или полуплоскости выделяется область где значениями напряжений или перемещений, возникающих от нагрузок, можно пренебречь. Форма элементов может быть различной и зависит от формы рассматриваемой области или ее участков. Сплошное тело, разделенное на элементы, казалось бы обладает большей податливостью, что приведет к искажению напряжения и деформаций. Для того чтобы этого не произошло, необходимо ввести определенные условия, приводящие к идентификации напряженно - деформируемого состояния тела. Это достигается требованиями выполнения условий сплошности, в частности, на границах между элементами, т.е. разделение сплошной среды на элементы в МКЭ не сопровождается ее разрезом, элементы не являются отдельными кусками, а лишь выделяются из сплошной среды. В общем случае среда может быть неоднородной по своим физико-механическим свойствам. Однако разбивку на элементы следует производить так, чтобы в пределах одного элемента участок среды можно было бы рассматривать как однородный. Причем любой другой элемент, оставаясь так же однородным, может характеризоваться свойствами, отличными от соседних элементов. Таким образом, система элементов будет в целом представлять неоднородную среду.
В данной работе будут использоваться так называемые симплекс элементы, т.е. элементы, которым соответствуют полиномы, содержащие константу и линейные члены:
0+1x+2y или 0+1x+2y+3z. (*)
Очевидно, число коэффициентов в полиномах (*) больше на единицу размерности координатного пространства, или, коэффициентов столько, сколько вершин в конечном элементе. Деформация и напряжение в любом конечном элементе выражаются через перемещения по известным формулам. В узлах элементов вводятся концентрированные силы, статически эквивалентные напряжениям на границе соответствующего элемента и внешним силам, приложенным к нему. Условие равновесия этих сил в любом узле приводит к системе алгебраических уравнений для основных неизвестных параметров.
3..4 Построение конечно - элементных соотношений
для плоской задачи линейной теории упругости
Метод конечных элементов позволяет оперировать с элементами различной формы. Однако для практических задач наиболее удобными оказались треугольные элементы, позволяющие легко сгущать сетку в местах ожидаемых высоких градиентов и удобные при оконтуривании границ рассчитываемой области.
Рассмотрим метод построения конечно-элементных соотношений для плоской задачи теории упругости. Пусть имеется двумерный симплекс-элемент с вершинами i, j, k (рис.3.4). Перемещение, как величина векторная, в каждом узле будет представлено двумя компонентами: {}T={U ,V}, где = i, j, k . Следовательно, для всего конечного элемента получим: {}еT={i ,J,k}. Перемещения внутри элемента должны однозначно определятся этими шестью компонентами. Для рассматриваемого элемента Функция перемещений представляется линейным полиномом: =1+2x+3y.
Рис. 3.4. Схема деформации треугольного элемента.
Поэтому для аппроксимации перемещений в области элемента будем иметь:
u =1+2 x+3 y, (3.6)
v =4+5 x+6 y,
где 1,2,....,6 - параметры линеаризации, постоянные для элемента.
Для компоненты U в узлах i, j, k получим:
ui =1+2 xi+3 yi,
uj=1+2 xj+3 yj, (3.7)
uk =1+2 xk+3 yk,
Решая (3.7) относительно 1,2,3 и подставляя полученные для них выражения в (3.6) получим:
U = [(ai+bi x+ciy)Ui + (aj+bj x+cjy)Uj+ (ak+bk x+cky)Uk] /(2S), (3.8)
где ai = xj yk - xk yi ; bi = yj-yk ; ci = -(xj – xk) (3.9)
Коэффициенты aj, bj, ak, bк, получаются из (3.9) круговой перестановкой индексов. Для компоненты V аналогично:
V = [(ai+bi x+ciy)Vi + (aj+bj x+cjy)Vj+ (ak+bk x+cky)Vk] /(2S), (3.10)
Деформации определяются из уравнений Коши
{}= = . (3.11)
П одставляя в уравнения (3.11) равенства (3.8) и (3.10) и производя дифференцирование, получим :
Напряжения {}= [D]*{} = [D]*[B]*{ }e (3.13)
Матрица [D] для случая плоского напряжённого состояния имеет вид:
[D]= ;
Для вывода основного уравнения МКЭ воспользуемся принципом стационарности
полной энергии системы:
( 3.14)
где П = 0,5* {}T{}dS – {}T{P} (3.15)
{P} - вектор внешних сил.
Подставим в (3.15) соотношения (3.12) и (3.13):
П = 0,5* {}T[B]T[D] [B] {}e dS - {}T{P}, (3.16)
т.к. в (3.16) векторы и матрицы от координат не зависят, то
П = 0,5*{}T[B]T[D] [B] {}e S - {}T{P}, (3.17)
Подставляя (3.17) в (3.14), получим:
S*[B]T[D] [B] {}e - {P} = 0, (3.18)
Введём обозначение: [K] = S*[B]T[D] [B], (3.19)
тогда (3.18) примет вид: [K] { }e ={ P} (3.20)
Это и есть основное уравнение МКЭ.
Если в (3.19) выполнить матричные операции, то получим:
где [K]ij = [K]jiT,
m, n = i, j, k ; = E / (1-2).
Матрица жесткости системы, определяющая жесткость конструкции в целом, строится в следующем порядке
Kij= , (3.23)
где N - количество конечных элементов дискретизованной области.
Суммирование в (3.23) производится учитывая номера конечных элементов и узлов. Покажем сказанное на примере. Пусть дискретизованная область имеет вид рис. 3.5.
Матрица жесткости для этого ансамбля(глобальная матрица)
может быть получена следующим образом:
1) формируем соответствующее матричное поле, пока не заполненное (двумерный массив);
2) рассматриваем i-й конечный элемент и для него вычисляем матричные коэффициенты по (3.22), которые прибавляем к глобальной матрице на позиции, определяемые глобальными
Рис. 3.5 – Схема номерами рассматриваемых узлов.
дискретизации Следовательно, соблюдая глобальную нумерацию узлов и
номера конечных элементов (верхние индексы), получим:
(3.24)
Ввиду того, что матрица жесткости системы имеет ленточную структуру и симметрична относительно главной диагонали, в памяти ЭВМ достаточно сформировать и хранить лишь члены из ленты по одну сторону от главной диагонали (включая и диагональные члены). После построения матрицы жесткости всей конструкции [K], составляется система линейных уравнений.
[K]{}= {P} (3.25)
где {P} и {} - векторы узловых сил и перемещений системы;
{P}= ; {}= . (3.26)
В своём первоначальном виде система (3.25) решения не имеет, т.к. матрица жёсткости [K] сингулярная – её главный определитель равен 0. Отметим, что построение основного уравнения метода конечных элементов (3.25) производилось без учёта граничных условий. Учёт граничных условий в (3.25) приводит к изменению матрицы жёсткости [K] и векторов {P} и {}. Матрица [K] уже не будет сингулярной и система (3.25) будет иметь решение.
После определения перемещения, по уравнениям (3.12) и (3.13) определяются деформации и напряжения в элементах.