Курсовая работа: Суперэлементное моделирование пространственной системы "плита – грунтовое основание"

алгоритмическая модель системы.

В соответствии с этим возникают следующие задачи:

определение (формирование) системы исследования;

построение математической модели системы;

разработка алгоритм решения модели.

2. Основные понятия теории упругости

2.1 Напряжения

????? ????????????? ????????? ??????: ??????? ?????????????? ????????, ??????? ????????????? ??? ????????? ???, ?????????? ?????????????? ?? ??? ???????? ????????, ?????????? ???? ???????????? ?????????? ?? ??????????? ??????? ??????? ??, ?????????? ????? ?????, ???????? ?????? ???????? ????????????? ???? ? ?? ??????? ??????????? ??????? F.

В общем случае напряжения по сечению распределены неравномерно, чтобы определить значение напряжения в некоторой точке этой плоскости, возьмём элементарную площадку δ F в окрестности данной точки и предположим, что силы, возникающие на этой площадке, сводятся к равнодействующей δ P . Если теперь равномерно стянуть элементарную площадку δF, то в пределе получится отношение δP/δF, которое определит величину напряжения, возникающего на плоскости АВ в некоторой точке. Направление этого напряжения будет совпадать с направлением равнодействующей δP. В общем случае напряжение направлено под некоторым углом к площадке δF, на которой оно действует, и обычно раскладывается на две составляющие: нормальное напряжение, перпендикулярное к площадке δF, и касательное напряжение, действующее в плоскости площадки.

2.2 Деформации

При рассмотрении деформации в упругом теле предполагается, что

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

Рассмотрим бесконечно малый элемент dxdy вблизи точки О тела.

Можно показать, что относительное удлинение по направлению оси y задается производной.

Рассмотрим теперь изменение угла между отрезками ОА и ОВ, которые до деформирования тела были взаимно перпендикулярны. Если u и v - перемещения точки О в направлениях осей x и y, то перемещения точки А в направлении оси у и точки В в направлении оси х будут соответственно равны. Поэтому первоначально прямой угол АОВ между отрезками ОА и ОВ уменьшается на величину, которая представляет собой деформацию сдвига между осями х и у.

3. Основная концепция метода конечных элементов

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

Симплекс элементам соответствуют полиномы, содержащие константу и линейные члены:

;

Здесь коэффициентов столько сколько узлов.

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

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

Это соотношение содержит шесть коэффициентов, поэтому рассматриваемый элемент должен иметь шесть узлов.

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

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

4. Характеристики тетраэдрального элемента

4.1 Функции перемещений

На фигуре 1 изображен тетраэдральный элемент ijpm в системе координат x, y, z.

Перемещение любой точки определяется тремя компонентами u, v, w в направлениях координат x, y, z. Таким образом, вектор перемещений имеет вид

. (1)

Если для задания линейного закона изменения какой-либо величины в плоском треугольном элементе требовались три узловых значения, то в трехмерном случае необходимо задать четыре узловых значения. По аналогии с представлением (4.3) можно записать, например,

. (2)

Приравнивая эти выражения перемещением узловых точек, получаем четыре уравнения типа

и т.д. (3)

из которых определяются коэффициенты .