- •Математический факультет Кафедра вычислительной математики и программирования
- •Содержание
- •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 Метод сеток для уравнения гиперболического типа
8 Верификация технологии и программного
ОБЕСПЕЧЕНИЯ КОМПЬЮТЕРНОГО МОДЕЛИРОВАНИЯ
СОСТОЯНИЯ СИСТЕМ ДЕФОРМИРУЕМЫХ ТВЁРДЫХ ТЕЛ
8.1 Методологические аспекты верификации
Верификация разработанных программных комплексов «Энергия-09» и технологии компьютерного объектно-ориентированного моделирования состояния систем деформируемых твёрдых тел проводилась методом вычислительного эксперимента. Для оценки адекватности построенной модели системы и конкретного программного комплекса по определению состояния исследуемой системы использовались соответствующие реальные задачи механики грунтов и строительной механики, для которых есть теоретические и (или) экспериментальные данные. Все задачи сводились к третьей краевой задаче математической физики и решались в линейной и нелинейной постановке методом конечных элементов в сочетании с методом энергетической линеаризации.
8.2 Верификация технологии и программного обеспечения
моделирования двумерных систем механики грунтов
Физическая постановка задачи
В качестве модельной задачи был принят фундамент из большеразмерных коробчатых плит на нелинейно-деформируемом грунтовом основании. Необходимо исследовать и оценить осадку плиты. Предполагается, что по всей контактной поверхности имеется полное сцепление материала плиты с грунтом. Это значит, что вся полость коробчатой плиты должна быть заполнена грунтом. На верхнюю плоскость фундаментной плиты действует нормальная внешняя нагрузка. Грунт, находящийся в полости, будет уплотняться, в основном, за счёт сжатия, его горизонтальные перемещения в поперечном направлении будут ограничены боковыми стенками полости. Продольные перемещения в силу геометрических размеров плиты можно принять нулевыми. Эти обстоятельства позволяют решение поставленной пространственной задачи свести к решению соответствующей плоской задачи. При достижении определённой степени уплотнения грунт внутри полости будет передавать нагрузку на нижележащие слои грунтового основания. При этом плитный коробчатый фундамент с уплотнённым в полости грунтом, возможно, будет работать как единый конструктивный элемент системы «фундамент - грунтовое основание».
В настоящей работе при компьютерном моделировании приняты следующие физико-механические характеристики элементов рассматриваемой системы «Плитный коробчатый фундамент – грунтовое основание». Модуль упругости для фундамента E = 36000 МПа (360000 кг/см²), для грунта E=36 МПа (360 кг/см²); коэффициент Пуассона для фундамента µ= 0,02, для грунта µ= 0,2.
Поскольку исследуется эффективность фундаментов из большеразмерных коробчатых плит, то базовой задачей для сравнительного анализа принят сплошной плитный фундамент, имеющий внешние геометрические и физико-механические характеристики аналогичные исследуемому типу фундамента, рисунок 8.1. Материалоёмкость базового и исследуемого типов фундамента будет различной, что и определит эффективность фундаментов из большеразмерных коробчатых плит [20, 30, 34, 35, 37, 39].
Компьютерное моделирование
Модельная задача №1 (Плитный сплошной фундамент – базовая задача)
Исходные данные: векторы шагов дискретизации
hx = {70, 50, 40, 40, 40, 50, 50, 40, 40, 40, 50, 70},
hy = {40, 40, 50, 50, 50, 50, 50, 60, 60, 80};
Eгр = 36 МПа, гр = 0.2; Епл = 36000 МПа, пл = 0,02,
расчётная область: 580 х 530 см; плита - 260 х 230 см;
нагрузка Р = 16,2 кН (1620 кгс), {P} = {135, 270, 270, 270, 270, 270, 135}.
Модельная задача №2 (Плитный коробчатый фундамент)
Исходные данные: векторы шагов дискретизации
hx = {70,50,40,40,40,50,50,40,40,40,50,70}; hy = {40,40,50,50,50,50.50,60,60,80};
Eгр = 36 МПа , гр = ,.2; Епл = 36000 МПа; пл = 0,02;
нагрузка Р = 16,2 кН (1620 кгс); {P} = {135, 270, 270, 270, 270, 270, 135}.
расчётная область: - 580 х 530 см; плита - 260 х 230 см; полость - 180 х 150 см.
Рисунок 8.1 – Постановка базовой задачи и дискретизация расчётной области
Анализ результатов компьютерного моделирования осадки плитного
коробчатого фундамента на лёссовидном грунтовом основании
Поставленная задача решалась при условии линейного и нелинейного деформирования элементов исходной системы. На рисунках 8.1, 8.2, 8.3 срединная часть поперечного сечения плиты показана узлами 44 – 48 и 57 – 61. Центральная часть полости плиты представлена узлами: 45-47, 58-60.
По значению перемещений в центральных узлах конструкции можно определить степень уплотнения грунта в полости коробчатой плиты, что показано на рисунке 8.2 и в таблице 8.1, где размеры поперечного сечения плиты приняты равными 260 х 30 см.
Рисунок 8.2 – Коробчатый фундамент и дискретизация расчётной области
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
19 |
20 |
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 |
30 |
31 |
32 |
33 |
34 |
35 |
36 |
37 |
38 |
39 |
40 |
41 |
42 |
43 |
44 |
45 |
46 |
47 |
48 |
49 |
50 |
51 |
52 |
53 |
54 |
55 |
56 |
57 |
58 |
59 |
60 |
61 |
62 |
63 |
64 |
65 |
66 |
67 |
68 |
69 |
70 |
71 |
72 |
73 |
74 |
75 |
76 |
77 |
78 |
79 |
80 |
81 |
82 |
83 |
84 |
85 |
86 |
87 |
88 |
89 |
90 |
91 |
92 |
93 |
94 |
95 |
96 |
97 |
98 |
99 |
100 |
101 |
102 |
103 |
104 |
105 |
106 |
107 |
108 |
109 |
110 |
111 |
112 |
113 |
114 |
115 |
116 |
117 |
118 |
119 |
120 |
121 |
122 |
123 |
124 |
125 |
126 |
127 |
128 |
129 |
130 |
131 |
132 |
133 |
134 |
135 |
136 |
137 |
138 |
139 |
140 |
141 |
142 |
143 |
Рисунок 8.3 – Схема дискретизации расчётной области
Как видно из таблицы, для внутренних узлов сплошного и коробчатого плитных фундаментов характерны незначительные горизонтальные перемещения. Сплошной плитный фундамент перемещается как твёрдое тело, при заданных условиях его осадка составила 1,64 см. При тех же исходных условиях осадка коробчатого плитного фундамента 1,7см. Это увеличение осадки произошло за счёт уплотнения грунта внутри полости коробчатого фундамента, в верхней части полости уплотнение несколько большее, чем уплотнение в нижней части, что показано в таблице 8.1 Значение осадок при других размерах поперечного сечения сплошной и коробчатой плит и их экономическая эффективность приведены в таблице 8.2.
Таблица 8.1 - Значения горизонтальных и вертикальных составляющих перемещений
в узлах деформируемой области плитных фундаментов (в см.)
№№ узлов |
Плитный сплошной фундамент |
Плитный коробчатый фундамент |
||||||
Линейное решение |
Нелинейное решение |
Линейное решение |
Нелинейное решение |
|||||
u |
v |
u |
v |
u |
v |
u |
v |
|
6 |
- 0,002 |
1,64 |
- 0,003 |
2,62 |
- 0,003 |
1,700 |
- 0,003 |
2,74 |
7 |
- 0,002 |
1,64 |
- 0,002 |
2,62 |
- 0,005 |
1,707 |
- 0,0006 |
2,75 |
8 |
- 0,003 |
1,64 |
- 0,002 |
2,63 |
- 0,008 |
1,71 |
- 0,004 |
2,76 |
32 |
- 0,006 |
1,64 |
- 0,005 |
2,62 |
- 0,015 |
1,7 |
- 0,017 |
2,74 |
33 |
- 0,006 |
1,64 |
- 0,006 |
2,62 |
- 0,013 |
1,7 |
- 0,014 |
2,75 |
34 |
- 0,006 |
1,64 |
- 0,006 |
2,63 |
- 0,011 |
1,71 |
- 0,012 |
2,76 |
45 |
- 0,008 |
1,63 |
- 0,01 |
2,62 |
- 0,035 |
1,64 |
- 0,05 |
2,65 |
46 |
- 0,008 |
1,64 |
- 0,01 |
2,62 |
- 0,014 |
1,62 |
- 0,019 |
2,61 |
47 |
- 0,008 |
1,64 |
- 0,01 |
2,62 |
- 0,011 |
1,65 |
- 0,016 |
2,65 |
58 |
- 0,001 |
1,63 |
- 0,015 |
2,62 |
- 0,045 |
1,57 |
- 0,064 |
2,51 |
59 |
- 0,001 |
1,64 |
- 0,015 |
2,62 |
- 0,005 |
1,51 |
- 0,006 |
2,42 |
60 |
- 0,001 |
1,64 |
- 0,015 |
2,63 |
- 0,02 |
1,56 |
- 0,033 |
2,51 |
Таблица 8.2 – Осадки сплошных и коробчатых фундаментных плит (в см.)
hпс |
230 |
180 |
130 |
80 |
hп |
150 |
100 |
50 |
0 |
Sлпс |
1,64 |
1,9 |
2,2 |
2,6 |
Sнпс |
2,62 |
3,2 |
3,9 |
4,8 |
Sлпк |
1,7 |
1,96 |
2,25 |
- |
Sнпк |
2,75 |
3,3 |
4,0 |
- |
С % |
45 |
38,5 |
26,6 |
- |
hпс – толщина сплошной фундаментной плиты;
hп - высота полости коробчатой фундаментной плиты;
Sлпс , Sнпс - осадки сплошной фундаментной плиты при линейном и нелинейном
деформировании основания;
Sлпк ,- осадки коробчатой фундаментной плиты при линейном и нелинейном
деформировании основания;
С, % - экономическая эффективность коробчатой фундаментной плиты.
В целом при размерах поперечного сечения плиты равными 260 х 230 см несущая способность коробчатого плитного фундамента уменьшилась на 3,6 % при упругом деформировании и на 4,6 % при неупругом деформировании грунтового основания. Эти данные говорят о сопоставимой несущей способности рассматриваемых типов фундаментов. Экономический эффект будет определяться размерами полости коробчатого плитного фундамента и может быть значительным, что показано в таблице 8.2.
Исследовалось состояние грунта в области выреза и в области, примыкающей к боковым стенкам плитного фундамента. Выборка из приведенных таблиц для грунта, отстоящего от боковой поверхности конструкции фундаментной плиты на один шаг дискретизации (40 см) и для центра внутренней полости коробчатого фундамента представлены в таблицах 8.2 и 8.3. Рассмотрены варианты устройства сплошного плитного фундамента из бетона (колонка «бетон»), полость коробчатого фундамента заполнена грунтом (колонка «грунт»). В обоих случаях значения перемещений сопоставимы, а в области, примыкающей к боковым стенкам плитного фундамента, наблюдается телескопический сдвиг грунта.
Таблица 8.3 – Горизонтальные и вертикальные перемещения грунта
в точках, отстоящих от боковой поверхности плиты на один шаг дискретизации
-
№№
узлов
бетон
грунт
u
v
u
v
3
16
29
42
55
0,28
0,05
0,02
0,01
0,04
1,21
1,23
1,21
1,18
1,14
0,27
0,05
0,03
0,08
0,12
1,25
1,27
1,25
1,22
1,18
Таблица 8.4 – Горизонтальные и вертикальные перемещения грунта в
точках полости коробчатого фундамента по оси симметрии
-
№№
узлов
бетон
грунт
u
v
u
v
7
20
33
46
59
0,03
0,015
0,001
0,02
0,04
1,7
1,7
1,7
1,7
1,7
0,003
0,02
0,04
0,03
0,002
1,82
1,82
1,82
1,72
1,57
Анализ полученных результатов позволяет сделать следующие выводы:
1 Грунт в области выреза находится в уплотнённом состоянии. Степень уплотнения грунта наибольшая вблизи верхней плиты и наименьшая снизу.
2 Перемещения грунта в области выреза незначительно отличаются от перемещений плитного фундамента.
3 Плитный фундамент вместе с уплотнённым грунтом в области выреза может рассматриваться как единая конструкция.
4 Несущая способность плитного и коробчатого фундаментов при одинаковых исходных данных отличаются незначительно.
Общий вывод
Из сделанных выводов следует: факт близкого соответствия решений для плитного и коробчатого фундаментов при одинаковых исходных данных может быть предпосылкой для разработки инженерного метода определения осадки плитного коробчатого фундамента.