- •Г.Г КАШЕВАРОВА, Т.Б. ПЕРМЯКОВА
- •Предисловие
- •Введение
- •Общие сведения о математическом моделировании.
- •Численные методы
- •Элементы теории погрешности
- •Понятия мастера и надстройки
- •Глава 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
- •Литература
- •Оглавление
- •Численные методы решения задач строительства на ЭВМ
В настоящее время существует большое разнообразие в формулировках МКЭ (метод Галеркина, метод наименьших квадратов, метод глобального баланса и др.[26]). Классический подход - вариационная формулировка метода, в соответствии с которой минимизируется функционал, связанный с физическим смыслом решаемой задачи, и тогда МКЭ можно рассматривать как один из вариантов метода Ритца, в котором используются
специфические локальные координатные функции
7.1.Основные положения МКЭ
Вобщем случае алгоритм МКЭ состоит из следующих этапов.
1.Построение расчетной модели для анализа конструкции или физического явления.
2.Разбиение системы на элементы (дискретизация области).
3.Аппроксимация искомой функции на каждом конечном элементе полиномом.
4.Составление разрешающих уравнений МКЭ.
5.Решение системы линейных алгебраических уравнений относительно неизвестных узловых значений искомой функции.
Построение расчетной модели
Расчетный анализ любой конструкции начинается с построения ее расчетной схемы, встречающейся в учебниках по строительной механике и сопротивлению материалов [11, 35], т.е. с попытки установить, что именно в рассматриваемом случае является существенным, а чем можно пренебречь.
При создании расчетной схемы отображается геометрия системы, условия закрепления, принимается решение о том, какой расчет будет выполняться: линейный или нелинейный', одномерный, плоский или пространственный; следует ли учитывать силы инерции и выполнять динамический расчет или можно ограничиться статическим анализом.
После того как расчетная схема установлена, наступает период ее детального описания в форме, пригодной для выполнения расчетного анализа. МКЭ требует дискретизации расчетной схемы. В связи с этим в дополнение к общепринятому содержанию понятия расчетной схемы введем понятие расчетной модели системы
Расчетной моделью системы (сооружения) будем называть совокупность отдельных конструктивных элементов, соединенных конечным числом шарнирных и жестких узлов и нагруженных внешней нагрузкой (это могут быть силы, перемещения, температура и др.).
Это определение полностью соответствует расчетной схеме стержневых систем. Для более сложных конструктивных элементов (пластины, оболочки, массивные тела) оно также подходит, только при этом усложняются исходные матрицы, определяющие упругие свойства элементов.
При создании расчетной модели реальной конструкции необходимо тщательно следить за тем, какие внутренние усилия могут возникать в применяемых конечных элементах и как эти
усилия согласуются друг с другом в общих узлах
Такие проблемы чаще всего возникают в местах сопряжения конструктивных элементов, имеющих разную размерность, например, при сопряжении одномерных стержней с двухмерными пластинами или с трехмерными с массивами и т.п.
■ В качестве примера рассмотрим конструкцию (рис.7.1) [48], в которой в месте примыкания ригеля к массивной стене при действии на ригель вертикальных и горизонтальных нагрузок, должен возникать
изгибающий момент. Моделируя ригель стержневым элементом
(каждый узел которого имеет 3 степени свободы), а стену - плоскими элементами «балки-стенки» (узлы которой имеют по 2 степени свободы), момент в ригеле в месте примыкания будет отсутствовать (рис.7.1,а).
Для учета момента в точке примыкания необходимо запретить угловое перемещение путем налоэюения соответствующей связи. Это соответствует предположению о том, что изгибная жесткость балкистенки на несколько порядков превышает изгибную жесткость ригеля.
собой только в узлах расчетной схемы и граничат друг с другом без зазоров и наложений.
Дискретизация области (тела) включает задание числа, размеров и формы конечных элементов, которые используются для построения дискретной модели реального тела, а также нумерацию узлов и нумерацию конечных элементов.
Дискретизация одномерного тела почти тривиальна, так как она сводится только к делению отрезка на более короткие участки.
Для двухмерных областей наиболее часто используются элементы в форме треугольников и четырехугольников, причем внутри области они могут быть использованы одновременно. Элементы могут иметь как прямо-, так и криволинейные границы, что позволяет с достаточной степенью точности аппроксимировать границу любой формы.
Для трехмерных областей наиболее употребимы элементы в форме тетраэдра и параллелепипеда, которые также могут иметь прямолинейные и криволинейные границы.
Аппроксимация искомой функции
Непрерывная функция, определяемая в каждой конкретной задаче (перемещение, температура и т.п.), может задаваться произвольным образом, но чаще всего для этих целей используются полиномы, которые подбираются так, чтобы обеспечить непрерывность искомой функции в узлах на границах элементов. Например, значение непрерывной функции ср(с) в
произвольной точке |
е-vo конечного |
элемента |
можно |
аппроксимировать полиномом |
|
|
|
|
<р(с) =А {C)R + A0 |
|
(7.1) |
где А(е) - вектор-строка коэффициентов полинома; А0 - свободный член; R = R(x,y, z) - вектор координат в рассматриваемой точке.
Задача данного этапа далее заключается в определении неизвестного вектора А(е) и свободного члена А0. Для этого, используя условие непрерывности функции в узлах, коэффициенты
—(to полинома выражают через вектор узловых значений функции Ф и
координаты узлов, в результате чего получают аппроксимирующую функцию в виде