- •Г.Г КАШЕВАРОВА, Т.Б. ПЕРМЯКОВА
- •Предисловие
- •Введение
- •Общие сведения о математическом моделировании.
- •Численные методы
- •Элементы теории погрешности
- •Понятия мастера и надстройки
- •Глава 1. Основные понятия матричного исчисления. Матрицы в расчетах строительных объектов
- •1.1. Матрицы и векторы. Определения
- •1.2. Матрицы специального вида
- •1.3. Действия над матрицами
- •1.4. Нормы матрицы и вектора
- •1.5. Матрицы в задачах строительной механики
- •1.5.1. Матрицы влияния внутренних сил
- •1.5.2. Матричная форма расчета статически определимых ферм
- •1.5.3, Матричная форма метода сил
- •1.5.4. Матричная форма метода перемещений
- •1.6. Матрицы в расчетах инженерных сетей
- •1.7. Функции Excel для операций над матрицами
- •Категория: математические. Функции:
- •2.1. Системы линейных алгебраических уравнений
- •2.2.1. Метод Гаусса
- •2.2.2. Метод Гаусса для СЛАУ с ленточными матрицами
- •2.2.3. Метод прогонки
- •2.2.4. Метод (схема) Холецкого
- •2.3. Итерационные методы решения систем линейных алгебраических уравнений
- •2.3.1. Метод Якоби (простых итераций)
- •2.3.2. Метод Гаусса - Зейделя.
- •2.3.3. Условия сходимости итерационного процесса
- •2.5. Обусловленность задач и вычислений, или как узнать, что получены правильные ответы
- •2.6. Вычисление определителя
- •2.7. Вычисление обратной матрицы
- •2.8. Нахождение собственных значений и собственных векторов матриц
- •2.8.1. Вводные замечания
- •2.8.2. Методы развертывания вековых определителей
- •2.8.3. Итерационные методы определения максимального по модулю собственного значения
- •2.9.1. Реализация метода Гаусса средствами приложения Excel
- •Последовательность действий:
- •2.9.4. Реализация метода Зейделя средствами приложения Excel
- •Последовательность действий:
- •3.1. Отделение корней
- •3.2. Этап уточнения корня
- •3.2.1. Метод половинного деления (бисекций)
- •3.2.2.Метод хорд
- •3.2.3. Метод Ньютона (метод касательных)
- •3.2.4. Модифицированный метод Ньютона
- •3.3. Системы нелинейных уравнений
- •3.4. Реализация численных методов решения нелинейных уравнений средствами приложения Excel
- •3.4.1. Решение нелинейных уравнений
- •Последовательность действий
- •4.2.3. Интерполяционный полином Эрмита
- •4.2.4. Сплайн-интерполяция
- •Глава 4. Аппроксимация
- •4.1. Задача и способы аппроксимации
- •4.2. Интерполирование функций
- •4.2.1. Постановка задачи интерполирования
- •4.2.2. Интерполяционная формула Лагранжа
- •4.3. Среднеквадратичное приближение функций
- •4.3.1. Метод наименьших квадратов
- •4.3.4. Квадратичное (параболическое) приближение
- •4.3.4. Эмпирические формулы с двумя параметрами. Метод выравнивания
- •4.4. Решение задач аппроксимации с помощью электронных таблиц Excel
- •4.4.1. Построение линейной эмпирической формулы методом наименьших квадратов
- •Последовательность действий
- •Последовательность действий
- •5.1. Квадратурные формулы прямоугольников
- •5.2. Квадратурная формула трапеций
- •5.3. Квадратурная формула Симпсона
- •5.4. Реализация методов численного интегрирования средствами приложения Excel
- •Глава 6. Численные методы решения дифференциальных уравнений с начальными и краевыми условиями
- •6.1.1. Задачи Коши и краевые задачи
- •6.2.1.Классификация уравнений и типы задач
- •6.3. Численные методы решения задач Коши
- •6.3.1. Метод Эйлера
- •(геометрический метод решения задачи Коши)
- •6.4. Численные методы решения краевых задач
- •Разностная схема краевой задачи для обыкновенного дифференциального уравнения 2-го порядка
- •Конечно-разностная аппроксимация функций двух переменных
- •Сходимость метода конечных разностей
- •6.5. Вариационный подход к решению краевых задач
- •6.5.1. Основные понятия вариационного исчисления
- •6.5.2. Связь решения краевой задачи с нахождением минимума функционала
- •6.5.3. Метод Ритца
- •6.6.1. Решение задачи Коши методом Эйлера
- •Построение второй итерации
- •Последовательность действий.
- •Порядок построения графиков приближенных решений краевой задачи
- •Глава 7. Метод конечных элементов
- •7.1. Основные положения МКЭ
- •Построение расчетной модели
- •Аппроксимация искомой функции
- •Составление разрешающих уравнений
- •Решение системы линейных алгебраических уравнений
- •7.2.1. Классификация конечных элементов
- •Одномерный конечный элемент
- •Двухмерные конечные элементы.
- •элемента
- •Одномерный симплекс-элемент
- •Двухмерный треугольный симплекс-элемент
- •7.2.3. Интерполирование векторных величин
- •7.2.4. Разбиение области на конечные элементы
- •7.2.5 Нумерация узлов и элементов
- •7. 3. Основные соотношения МКЭ
- •7.3.1. Получение разрешающих уравнений на примере плоской задачи теории упругости
- •7.3.2. Примеры разрешающих уравнений в задачах расчета строительных объектов
- •7.4. Другие типы конечных элементов
- •7.4.1. Элементы Эрмита
- •7.5. Теоретическая и практическая сходимость МКЭ
- •7.6.1. Специализированные программные комплексы
- •7.6.2. Универсальные программные комплексы
- •8.1.1. Математическая модель задачи оптимизации
- •8.1. Общие сведения
- •8.1.2. Необходимые и достаточные условия экстремума функции
- •8.1.3. Классификация задач математического программирования
- •8.2. Постановка задачи оптимального проектирования
- •8.2.1. Определение входных и выходных параметров
- •8.2.2. Выбор целевой функции
- •8.2.3. Назначение ограничений
- •8.2.4. Нормирование управляемых и выходных параметров
- •8.2.5. Примеры постановок задач оптимального проектирования
- •8.3. Задачи линейного программирования
- •8.3.1. Общая постановка задачи ЛП
- •8.3.2. Геометрический смысл системы линейных неравенств с двумя неизвестными
- •8.3.4. Симплекс-метод решения задач линейного программирования
- •Задача об оптимальном плане выпуска продукции
- •Задача об оптимальном раскрое материалов (о минимизации отходов)
- •Задача о планировании смен на предприятии
- •Задача о покрытии местности при строительстве объектов
- •Транспортная задача
- •Задача о назначениях (проблема выбора)
- •8.3.6. Двойственные задачи в линейном программировании
- •8.4. Нелинейные задачи оптимизации
- •8.4.1. Выпуклые множества и выпуклые функции
- •8.4.2. Классификация численных методов решения нелинейных задач оптимизации
- •Основные этапы поиска экстремума
- •8.4.3. Численные методы одномерного поиска
- •Метод перебора или равномерного поиска
- •Метод дихотомии (или половинного деления)
- •Метод квадратичной интерполяции
- •Метод покоординатного спуска
- •Метод градиентного спуска
- •Метод наискорейшего спуска
- •Метод сопряженных градиентов
- •Метод Ньютона
- •Метод штрафных функций
- •8.5. Решение задач оптимизации с помощью электронных таблиц Excel
- •Литература
- •Оглавление
- •Численные методы решения задач строительства на ЭВМ
где
|
|
' l l |
'*12 |
|
'1/, |
|
|
г = П |
- |
- 1 г 2 \ |
г 2 2 |
г 2 /1 |
- матрица жесткости (1.30) |
||
г 2 |
|||||||
' • " Г |
|
|
|
||||
|
|
|
г п 2 |
|
Г/ш |
|
|
|
|
|
1 |
0 |
0 |
|
|
А = Л| |
А2 |
А/«]= |
0 |
1 |
0 |
= Е - единичная матрица. |
|
|
|
|
|||||
|
|
|
0 |
0 |
1 |
|
|
Следовательно, |
|
|
|
|
|||
|
|
|
г = 8"1, |
|
(1.31) |
т.е. матрица жесткости и матрица податливости взаимно обратны.
1.6. Матрицы в расчетах инженерных сетей
Матричная форма расчета является удобной и при расчете инженерных сетей. Для иллюстрации рассмотрим гидравлическую сеть, изображенную на рис. 1.6.
В случае медленных (ламинарных) течений поток Q через поперечное сечение трубы пропорционален разности давлений в начале и конце трубы. Таким образом, для элемента е6 (рис. 1.7) потоки в эту трубу в узлах 2 и 3 соответственно будут
Яг = C<V' (Рг - Рз) - |
Яз = - с<!" (Рг ~ Рз) Л 1 -32) |
где р 2 и Рз - давления в узлах 2 и 3,
Qi и Оз - потоки в тех же узлах,
с - постоянная, зависящая от свойств жидкости, диаметра и длины трубы.
Рис. 1.7. Элемент (труба) е6
В матричной форме уравнения (1.32) можно записать
|
|
|
Рг |
(1.32 о) |
|
|
|
|
Рз |
||
|
|
|
|
|
|
или |
|
|
|
|
|
кС(' |
к4 |
Рг |
|
|
|
Л22 |
Л23 |
|
(1.32 |
б) |
|
к4 |
к4 |
Рз |
|
||
.32 |
*33 |
|
|
|
Уравнение (1.32 б) - это матричное уравнение для элемента ей. Оно также может быть записано в виде расширенного матричного уравнения, включающего все узлы сети:
' 0 ■ |
'о |
0 |
0 |
|
0 |
0 |
0 |
0 |
0' ~Р\~ |
|
|
Q? |
0 |
к 4 |
к с6 |
0 |
0 |
0 |
0 |
0 |
Рг |
|
|
|
Л22 |
л 23 |
|
|
|
|
|
|
|
||
Q? |
0 |
к 4 |
к е“ |
0 |
0 |
0 |
0 |
0 |
Pi |
|
|
|
л 32 |
л зз |
|
|
|
|
|
|
|
||
0 |
0 |
0 |
0 |
|
0 |
0 |
0 |
0 |
0 |
Ра |
(1.32 в) |
0 |
0 |
0 |
0 |
|
0 |
0 |
0 |
0 |
0 |
Ръ |
|
|
|
||||||||||
0 |
0 |
0 |
0 |
|
0 |
0 |
0 |
0 |
0 |
Рб |
|
0 |
0 |
0 |
0 |
|
0 |
0 |
0 |
0 |
0 |
Рг |
|
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
0 |
Рхi . |
|
или в форме |
|
|
|
|
|
|
Q v =к' Р, |
|
(1.32 г) |
где |
р - вектор, |
компоненты |
которого рх, рг, ...,р 8 равны |
|
давлениям в узлах сети. |
|
|
||
|
Предположим, что жидкость поступает в сеть в узлах 1,2,..., |
|||
8 с |
расходами |
Rt, R2 |
Rs соответственно. |
Уравнение |
неразрывности для узла 2, например, имеет вид |
|
|||
|
*2 |
=О+О +О + 0 2е‘ +0 + Q? +0 + Q? |
(1.33) |
|
|
/= | |
|
|
|
Система уравнений типа (1.33) может быть записана следующим образом:
|
' К |
8 |
' o f |
|
|
|
R = |
*г |
вг |
1У |
(1.34) |
||
- Z |
||||||
|
|
|||||
|
|
|
|
|||
|
А . |
|
QI. |
|
|
Подстановка уравнений типа (1.32 г) в (1.34) приводит к
матричному уравнению системы
Kp=R. |
(1.35) |
Для заданных подводимых потоков узловые давления могут быть найдены путем решения уравнения (1.35). После этого расходы
через каждую трубу можно вычислить с помощью уравнений типа
(1.32).
■Пример 1.6. Пусть сеть труб, показанная на рис. 1.6, является частью системы водоснабжения поселка. Из-за повреждения дамбы давление в системе мало и поток в трубах считается ламинарным. Основываясь на этом, вычислим вытекающие потоки в узлах 6, 7 и 8 для следующих данных:
элемент es (подводящая магистраль) |
d=5 см, |
/= 1200 м; |
|
элементы eh е3, es |
|
rf=2,5 см, /= 900 м; |
|
элементы е2, |
e 7 |
d=2,5 см, |
/= 450 м; |
Давление в узле 5 |
|
|
2 кг/см2 |
Давление в узлах 6, 7, 8 |
атмосферное (1 кг/см2). |
||
Перепад давления Ар для ламинарного потока в трубе диаметром d и |
|||
длиной / вычисляется по уравнению Гагена - Пуазейля |
|
||
|
Др =32ф / d 2, |
|
(1.36) |
где р - динамическая вязкость и v - средняя скорость, с использованием согласованной системы физических единиц.
Предполагается, что вода в системе имеет постоянную температуру 16°С, динамическая вязкость равна 1,1-10-4 кгс/м2, все трубы горизонтальны и утечка в стыках пренебрежимо мала.
Из эксперимента известно, что поток в трубе ламинарный, если число Рейнольдса Re = pvd/\i, не превышает 2000. Для болыших значений числа Рейнольдса связь между перепадом давления и расходом становится нелинейной.
Принимая плотность воды 1О*3 кг/см3, проверим корректность предположения о ламинарности потоков в трубах. Расходы на выходе из узлов 6, 7 и 8 равны 830, 350 и 480 см3/С; числа Рейнольдса изменяются от 3100 в элементе е5>до 40000 в элементе е6 и, следовательно, предположение о ламинарности потока некорректно.