- •Г.Г КАШЕВАРОВА, Т.Б. ПЕРМЯКОВА
- •Предисловие
- •Введение
- •Общие сведения о математическом моделировании.
- •Численные методы
- •Элементы теории погрешности
- •Понятия мастера и надстройки
- •Глава 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
- •Литература
- •Оглавление
- •Численные методы решения задач строительства на ЭВМ
Задача об оптимальном раскрое материалов (о минимизации отходов)
В строительном деле достаточно часто приходится решать проблему раскроя полуфабрикатов или минимизации отходов производства.
Постановка задачи. Из имеющихся заготовок в виде досок длиной D каждая требуется получить 6, частей длиной L, (z=l,2,...,m). Имеется несколько вариантов раскроя Vj (/=1,2,...,п) каждой доски L (план раскроя).
При каждому-м варианте раскроя получается ау частей длиной lj. (При этом axjlx + a2JI2 +... + a -lт < D . Это условие, наложенное
на коэффициенты, содержится в определении «вариант раскроя» и не является условием оптимизации).
Требуется так распилить доски, чтобы было как можно меньше отходов, или требуемое количество частей должно быть получено из минимального количества заготовок.
В качестве независимых параметров выбираем Xj - число досок (заготовок), которое намечено распилить согласно у-му варианту.
Запись математической модели задачи:
минимизировать
|
|
|
л |
(8.67) |
Zmin(X) = *| + *2 + - + Х„ = |
||||
|
|
|
7=1 |
|
при ограничениях: |
|
|
|
|
а\\Х\ |
а12*2 + |
|
а1пХп —^1» |
|
а 21Х 1 + |
0 22 х 2 + |
••• + |
а 2пХ п — ^ 2 > |
(8.68) |
|
|
|
|
|
Я ,„ 1 * 1 |
а т 2 Х 2 + • • • + а т п Х п ~ |
’ |
||
|
Xj;> 0 |
при |
j=l,2,...,n. |
(8.69) |
Задачи об оптимальном раскрое довольно разнообразны и в качестве примера рассмотрим еще один вариант такой задачи.
■Пример 8.7. При серийном производстве некоторого изделия из полос проката длиной 5000мм необходимо вырезать 3 вида заготовок.
Номер, длина и количество заготовок: №1 длина 1655мм, 1шт.
№2 длина 1050мм, 5шт. №3 длина 210мм, 1шт.
Требуется составить оптимальный план раскроя, чтобы получить комплект заготовок для 12 изделий и израсходовать при этом минимальное количество полос.
Решение
Оптимизацию будем производить, исходя из минимальных отходов полос проката при их раскрое.
1) Составим таблицу-карту раскроя:
Способ |
Количество заготовок |
Полезно |
Длина |
Количе |
||
раскроя |
|
длиной(мм) |
используемая |
отходов |
ство |
|
|
|
|
|
длина (мм) |
(мм) |
полос |
|
1655 |
1050 |
210 |
|
|
|
1 |
3 |
0 |
0 |
4965 |
35 |
Х \ |
2 |
2 |
1 |
1 |
4570 |
430 |
х 2 |
X |
||||||
3 |
1 |
3 |
0 |
4805 |
195 |
* 3 |
4 |
0 |
4 |
1 |
4410 |
590 |
х4 |
Таким образом, получилось четыре способа раскроя полос.
2)В качестве проектных параметров возьмем: - количество прокатных полос, раскроенных /-способом.
3)Определим длину отходов при каждом способе раскроя.
4)В качестве функции цели примем суммарную длину отходов, которая должна быть минимальной:
Zn,in = 35xi + 430х2+ 195х3+ 590х4.
5) Запишем ограничения. Для 12 изделий необходимо заготовок: №1 - 12 шт.
№2 - 12 х 5 шт. №3 - 12 шт.
Ограничения записываем исходя из условий, что количество заготовок для 12 изделий должно быть не меньше соответственно 12, 60,
12:
Зх, + 2 х 2 + х 3 + 0 * х 4 >12
О х, + х2 +Зх3 + 4 х 4 >60
по количеству заготовок
х2 + х4 > 12
х,>0, / = 1,2,3,4
Задача о планировании смен на предприятии
Полученная математическая модель задачи об оптимальном раскрое материалов (8.67) - (8.69) может быть использована также при решении задач о планировании смен на предприятии.
Только в этом случае содержательный смысл параметров будет другой, а именно:
Vj (/=1,2,...,;/) - возможные в течение дня смены; |
|
|
Li (/=1,2,.. .,m) - |
определенное время дня. |
|
ctij=1, если У} |
предусматривает работу во время |
в |
противном случае <7,7=0;
bi~ число работников, требующихся в момент времени xj - количество работников смены Vj.
При планировании строительства производственных, социальных и др. объектов всегда стоит задача правильного их размещения, то есть возникает так называемая задача о покрытии. Эта задача является частным случаем рассмотренных выше задач о раскрое и планировании (в ней b,=1; а^е {0,1} ).
Задача о покрытии местности при строительстве объектов
Постановка задачи. Для строительства гидроэлектростанций могут быть использованы пять мест К, (/-1,2,...,5). Каждая ГЭС К, могла бы обслужить (яг//=1) или не обслужить {ау=0) некоторые из четырех регионов Li (/=1,2,...,4). После сооружения минимального количества электростанций каждый из четырех регионов должен быть «покрыт» хотя бы один раз.
В качестве управляемых переменных примем xj (/=1,2,...,5), причем xj=1, если в месте Vj построена ГЭС, и х/=0 в противном случае.
Математическая модель задачи получается такая же, как и в задаче оптимального раскроя (8.67) - (8.69), но с условием, что bf=1 (/=1,2,...,4) и с дополнительным условием <а/=0 или х/=\» при 7=1,2,-..,5).
Возможны следующие варианты задачи:
1) Заданы расходы на строительство VJy равные су, и ищется наиболее дешевое «покрытие» всех местностей.
Целевая функция в этом случае Zmin = с]х] +с2х2 +... + с5х5; 2) Можно потребовать также Л-кратное покрытие. В этом
случае 6, = к.
Транспортная задача
Это еще одна типичная и очень важная задача линейного программирования.
Постановка задачи. В п пунктах отправления (склады, заводы и т.п.) АиА2у..уАп имеется (или производится) некоторый однородный продукт в количествах аиа2у..,а„. Необходимо доставить этот продукт в т пунктов назначения B{yB2y..JBmв количествах b{,b2j..ybm.
Стоимость перевозки единицы груза из пункта At в пункт Bj равна Су. Заметим, что стоимость - это условное понятие, которое может означатьрасстояние, тариф, расход топлива, время и т.д.
Количество перевозимого груза из пункта А{ в пункт Bj обозначим через Ху (/=1,2,..,л; Задача заключается в определении таких величин Ху , при которых стоимость перевозок будет минимальной.
Условия задачи можно записать компактно в виде таблицы 8.1. (двойной матрицы):
|
|
|
|
Таблица 8.1 |
V |
в, |
Ъх |
Ъг |
Ьщ |
Л / |
|
|
|
|
|
а\ |
С\\ |
С\г |
С\т |
|
*11 |
|
*12 |
*1/и |
|
а2 |
сг\ |
^22 |
Clm |
|
*21 |
|
*22 |
*2ш |
|
ап |
Сп\ |
с п2 |
Спт |
|
х„\ |
|
Х„2 |
*/ш/ |
Совокупность т х п чисел Ху , т.е. матрицу X=[xv] будем называть
матрицей перевозок или планом перевозки, а матрицу С=[с,/ ] - матрицей транспортных издержек (затрат).
Сформулируем математическую модель транспортной задачи.
Транспортная задача заключается в отыскании среди допустимых планов оптимального, т.е. такого, по которому общая
стоимость перевозок минимальна, т.е
т п
Z min = Z |
Z C i/ •* !/ |
( 8 '7 ° ) |
/=1 |
у =1 |
|
План является допустимым, если числа Ху удовлетворяют
следующим естественным ограничениям: |
|
|
||||
|
хп + х12 +... + xin = а, |
(/ = 1,2,.., т) 1 |
|
|||
|
x \ j + x 2j |
+ - + x mj = |
b i |
U = U2,.., и)} |
|
|
Поскольку грузы предполагается перевозить только в одном |
||||||
направлении - |
из пунктов отправления в пункты назначения, то на |
|||||
переменные |
накладываются |
условия |
неотрицательности |
|||
переменных: |
Ху >0 |
(/ = 1,2,..,/и; |
|
у =1,2,..,и), |
(8.72) |
|
|
|
в которых первые т равенств означают, что из каждого пункта производства Aj вывозится весь произведенный продукт. Последние (и) равенства означают, что каждый пункт потребления полностью удовлетворяется.
Транспортную задачу в приведенной формулировке (8.70) - (8.72) называют закрытой (или замкнутой) транспортной задачей, в отличие от открытой, в которой
пт
! > |
. ' Х > - |
<8-73> |
/=1 |
J = 1 |
|
■ Пример 8.8. Построить математическую модель транспортной задачи.
На трех цементных заводах производится цемент одной и той же марки в количествах соответственно 100, 130, 170 тонн. Цемент следует доставить на четыре ЖБК, потребляющих его соответственно в количествах 150, 120, 80, 50 тонн. Стоимости (у.е.) перевозок одной тонны продукта с /-го (/=1,2,3) завода нау-й (/=1,2,3,4) ЖБК приведены в табл. 6.2. Спланировать перевозки так, чтобы их стоимость была минимальной.
|
|
|
|
|
|
Таблица 8.2 |
Цемент |
Стоимость перевозки |
_ |
Объем |
|||
(У-е.) |
|
^ |
|
|
производ |
|
ные |
|
|
—-""'""'Количество . |
|
ства (т) |
|
заводы |
__перевозимого продукта (т) |
ЖБК-4 |
|
|||
|
ЖБК-1 |
ЖБК-2 |
ЖБК-З |
|
||
№1 |
3 |
|
5 |
7 |
11 |
100 |
|
Хц |
Х12 |
Х13 |
X14 |
||
№2 |
1 |
|
4 |
6 |
3 |
130 |
|
Х2| |
Х22 |
Х23 |
Х24 |
||
№3 |
5 |
|
8 |
12 |
7 |
170 |
|
Х31 |
Хз2 |
Хзз |
х34 |
||
Объем |
150 |
|
120 |
80 |
50 |
|
пот |
|
|
|
|
|
|
ребления |
|
|
|
|
|
|
Проектные |
параметры: |
^ —количество (т) |
продукта, |
|||
перевозимого с /-го (/=1,2,3) завода нау-й (/-1,2,3,4) ЖБК. |
|
Тогда целевая функция (общая стоимость перевозок) имеет вид:
Zmi„ =3*11+5*12+7*13+1 1*14+*21+4*22+6*23+3*24+5*31+8*32+12*33+7*34.