- •Г.Г КАШЕВАРОВА, Т.Б. ПЕРМЯКОВА
- •Предисловие
- •Введение
- •Общие сведения о математическом моделировании.
- •Численные методы
- •Элементы теории погрешности
- •Понятия мастера и надстройки
- •Глава 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
- •Литература
- •Оглавление
- •Численные методы решения задач строительства на ЭВМ
Специальные средства предусмотрены для создания расчетных моделей, поверхность которых описывается аналитически. Для формирования произвольных сеток на плоскости используется автоматическая триангуляция, с помощью которой сетка может быть нанесена на любую область расчетной схемы.
Для того чтобы обеспечить инженеру работу с расчетной схемой в привычной среде, используются разбивочные (координационные) оси. На этих осях могут выполняться операции по созданию схемы и ее фрагментации. Они могут быть показаны на всей схеме или на любом ее фрагменте.
7.6.2.Универсальные программные комплексы
Ун и в е р с а л ь н ы е программные комплексы общего назначения позволяют исследовать широкий класс конструкций и физических явлений при различных условиях нагружения и свойствах материалов. Их отличительными признаками являются: а) сложная логическая структура; б) длительное время разработки; в) обширная библиотека конечных элементов; г) модульный принцип построения; д) наличие базы данных и СУБД; е) большой объем входной и выходной информации и др.
Примеры универсальных программных комплексов, которые используют для расчета строительных объектов: ANSYS, NASTRAN, COSMOS, STARK, COSAR и др. Рассмотрим один из них:
■УНИВЕРСАЛЬНЫЙ ПРОГРАММНЫЙ КОМПЛЕКС - ANSYS
[49]- это программа конечно-элементного анализа, разработанная доктором Джоном Свенсоном и сотрудниками фирмы Swanson Analysis Systems Incorporated (SASI) в 1970 году. С того времени программа стала необходимым инструментом инженерного проектирования и проведения анализа разработок методом конечных элементов. В настоящее время
программа ANSYS доступна для работы на всех компьютерах - от персональных компьютеров (PC) до суперЭВМ.
ANSYS - это ПК, который ставит своей основной целью предоставление возможности выполнения практически любого вида анализа с помощью метода конечных элементов в интересах фактически
любой отрасли производства: строительной, автомобильной, космической, машиностроительной, и др. Универсальность программы заключается в том, что ее можно использовать при решении задач во всех
технических дисциплинах: прочности, механики, теплофизики, гидромеханики и др. Возможности программы ANSYS реализуются благодаря отдельным программам, составляющим полный пакет ANSYS. Так, например,
•программа ANSYS/LinearPlus обеспечивает решение задач линейной статики (с нелинейными элементами зазора и большими прогибами балок и оболочек), динамики и устойчивости..
•ANSYS/Thermal содержит возможности теплового анализа.
•ANSYS/Multiphysic - представляет собой программное средство анализа, позволяющее проводить расчетные исследования не только в таких отдельных областях знания, как прочность, распространение тепла, механика жидкостей и газов или электромагнетизм, но и решать связанные задачи.
•ANSYS/Mechanical™ предоставляет широкие возможности для решения сложных задач прочности конструкций, теплопередачи и акустики. Эта программа является мощным инструментом для определения перемещений, напряжений, усилий, температур и давлений, а также других важных параметров.
•ANSYS/Structural™ выполняет сложный прочностной анализ конструкций с учетом разнообразных нелинейностей, среди которых геометрическая и физическая нелинейности, нелинейное поведение конечных элементов и потеря устойчивости. Используется для точного моделирования поведения больших и сложных расчетных моделей.
•ANSYS/FLOTRAN™ - это удобное программное средство для решения разнообразных задач гидродинамики, включая ламинарное и турбулентное течение несжимаемой или сжимаемой жидкости.
Впрограмме ANSYS можно работать в интерактивном или пакетном режиме. Полностью интерактивная графика (т.е. средства и системы ввода, отображения и редактирования изображений) является составной частью программы ANSYS.
Все функции, выполняемые программой ANSYS, объединены в группы, которые называются процессорами. Программа имеет один
препроцессор, один процессор решения, два постпроцессора и несколько вспомогательных процессоров, включая оптимизатор.
Препроцессор используется для создания конечно-элементной модели и выбора опций для выполнения процесса решения.
Процессор решения используется для приложения нагрузок и граничных условий, а затем для определения отклика модели.
С помощью постпроцессора пользователь обращается к результатам решения для оценки поведения расчетной модели, а также для проведения дополнительных вычислений, представляющих интерес.
Средства препроцессорного твердотельного моделирования программы ANSYS позволяют иметь дело непосредственно с геометрической моделью, не обращаясь к специфическим объектам (узлам и элементам) конечно-элементной модели. Для того чтобы облегчить генерацию модели, программа отделяет фазу задания геометрии и граничных условий от построения сетки конечных элементов. Сначала пользователь описывает геометрию (и граничные условия, если есть желание) твердотельной модели. Если делать это в интерактивном режиме, то возможна удобная проверка входных данных. Программа строит сетку для получившейся модели, что определяет местоположение узлов и связность элементов. После того как построена твердотельная модель, ее конечно-элементный аналог (т.е. сетка узлов и элементов) может быть создана всего лишь одним обращением к меню программы. В программе ANSYS предусмотрено четыре способа генерации сетки: монтаж из простых составных частей, независимый способ, метод экструзии и адаптивное построение.
Средства твердотельного моделирования включают в себя представление геометрии, основанное на использовании сплайновой технологии, геометрических примитивов и операций булевой алгебры. При желании можно отказаться от услуг твердотельного моделировщика и построить конечно-элементную модель непосредственным заданием узлов, элементов и граничных условий.
Пользователь получает результаты анализа на стадии решения, Эта фаза применения программы ANSYS состоит в задании вида анализа и его опций, нагрузок и шага решения и заканчивается запуском на счет конечно-элементной задачи. Выбранный вид анализа указывает программе, какие разрешающие уравнения следует использовать для решения данной задачи. Каждая категория расчетов включает несколько их отдельных типов, например, статический и динамический типы прочностных расчетов. Выбором опций можно дополнительно определить особенности проводимого анализа.
Заданные нагрузки и ограничения определяют граничные условия для расчетной модели. К нагрузкам относятся ограничения степеней свободы, сосредоточенные, распределенные, объемные и инерционные усилия. Конкретный вид нагрузок зависит от вида проводимого анализа (например, приложенная в точке нагрузка может быть сосредоточенной
силой при прочностном анализе или тепловым потоком при расчете теплопередачи).
На стадии получения решения имеется возможность изменить свойства материала и атрибуты конечного элемента. Пользователь имеет возможность выбирать те конечные элементы, которые могут менять степень полинома, тем самым уменьшая общее время решения.
После того как все соответствующие параметры заданы, может быть выполнено и само решение. Пользователь поручает программе решить определяющие уравнения и получить результаты для выбранного вида анализа. В вычислительном отношении это самая интенсивная часть анализа, не нуждающаяся, однако, во вмешательстве пользователя. Она требует самых значительных затрат компьютерного времени и минимальных затрат времени пользователя. Для того чтобы получить решение за минимальное время, программа ANSYS переупорядочивает расположение элементов и узлов.
Все типы расчетов, выполняемые программой ANSYS, основаны на классических инженерных представлениях и концепциях. При помощи численных методов эти концепции формулируются в виде матричных уравнений, которые наиболее пригодны для конечно-элементных приложений.
Совокупность дискретных областей (элементов), связанных между собой в конечном числе точек (узлов), представляет собой математическую модель системы, поведение которой нужно анализировать. Основными неизвестными являются степени свободы узлов конечно-элементной модели. К степеням свободы относятся перемещения, повороты, температуры, давления, скорости, потенциалы электрических или магнитных полей; их конкретное содержание определяется типом элемента, который связан с данным узлом. В соответствии со степенями свободы для каждого элемента модели формируются матрицы масс, жесткости (или теплопроводности) и сопротивления (или удельной теплоемкости). Эти матрицы приводят к системам совместных уравнений, которые обрабатываются так называемыми “решателями” По умолчанию для решения систем уравнений используется фронтальный решатель. Решатель этого типа одновременно формирует и решает общую для нескольких элементов матрицу жесткости, состоящую из индивидуальных матриц элементов. Эта процедура последовательно продвигается через всю модель, элемент за элементом, вводя уравнения, соответствующие степеням свободы отдельного элемента. В это же время определяются основные неизвестные
и исключаются (на основе метода Гаусса) из общей матрицы. С помощью алгоритмов прямого действия —типа фронтального решателя - получают “точное” решение системы, а итеративные решатели дают сходящееся от итерации к итерации, приближенное решение.
С помощью постпроцессорных средств программы пользователь имеет возможность легко обратиться к результатам решения и интерпретировать их нужным образом, используя обширный набор команд и функций дружественного интерфейса. Результаты решения включают значения перемещений, температур, напряжений, деформаций, скоростей и тепловых потоков. Итогом работы программы на постпроцессорной стадии является графическое и/или табличное представление результатов.
Глава 8. Численные методы оптимизации
Необходимость принятия наилучших решений так же стара, как само человечество.
Термин «оптимизация» в литературе обозначает процесс отыскания наилучшего или оптимального решения. На практике под оптимизацией понимается скорее стремление к совершенству, которое, возможно, не будет достигнуто по ряду объективных причин.
С оптимизацией явно или неявно мы встречаемся в любой сфере человеческой деятельности, в том числе и в строительстве. Методы оптимизации позволяют при решении организационно- управленческих задач получить наилучшее распределение ресурсов или оптимальные объемы выпуска продукции с целью получения максимальной прибыли; в задачах технического проектирования - выбрать наилучший вариант конструкции или сооружения и т.д.
Проблема оптимизации имеет два основных аспекта:
1)нужно поставить задачу или разработать ее
математическую модель в виде математических выражений, устанавливающих основные взаимосвязи оптимизируемых параметров, формализовав при этом понятие «оптимальный», придав количественные оценки возможным вариантам, количественный смысл словам «лучше», «хуже»;
2)нужно решить задачу, уже имеющую математическую формулировку, выбрав подходящий вычислительный процесс (чаще всего - численный метод), учитывающий все условия, содержащиеся в математической модели.
Постановка задачи оптимизации является наиболее сложным и ответственным этапом. Она осуществляется с учетом
назначения реального объекта, целей проектирования (управления) и конкретных условий реализации задачи. Далее на примерах наиболее типичных задач экономики, управления и проектирования мы рассмотрим приемы и особенности