Оптимизация конструкций и технологий (стр. 1 )

ОПТИМИЗАЦИЯ КОНСТРУКЦИЙ И ТЕХНОЛОГИЙ

Требования к проектируемым зданиям и сооружениям, к организации и технологиям их возведения постоянно повышаются. Если раньше, например, к конструкциям предъявлялись лишь требования прочности, жёсткости и устойчивости, то теперь необходимо, чтобы они были предельно лёгкими и, в то же время, экономичными. В связи с этим возникают задачи оптимального проектирования конструкций и программирования производственных процессов. Однако практическое применение даже давно разработанных методов математического программирования сдерживается их сложностью, громоздкостью, большими объёмами вычислений, требованиями высокой программистской квалификации к отраслевым специалистам.

Цель данного доклада состоит в иллюстрации того, что упомянутые затруднения легко преодолеваются с помощью популярного вычислительного комплекса MATLAB. Рассматривается класс задач проектирования, организации, технологий, управления производством, квалифицируемых как задачи линейного программирования.

В среде MATLAB задачи линейного программирования решаются с помощью функции linprog. Функция linprog решает задачу линейного программирования в форме

  (min),

[x, fval] = linprog (c, A, b, Aeq, beq, lb, ub) - основной вид.

Основными входными данными linprog являются: вектор коэффициентов целевой функции c, матрица ограничений-неравенств A, вектор правых частей ограничений-неравенств b, матрица ограничений-равенств Aeq, вектор правых частей ограничений-равенств beq, вектор lb, ограничивающий план x снизу, вектор ub, ограничивающий план x сверху. На выходе функция linprog даёт оптимальный план x задачи (1) и экстремальное значение целевой функции fval.

Если какой-то из входных параметров отсутствует, на его место следует поставить квадратные скобки [], за исключением случая, когда это последний параметр в списке. Например, если нужно решить задачу без ограничений равенств, то оператор вызова функции linprog будет выглядеть так:

[x, fval] = linprog(c, A, b, [], [], lb, ub).

I. Проектирование конструкций

1.Неразрезные балки

Пример 1. Заданы схема (рис. 1) и две нагрузки для двухпролётной жёстко-пластической неразрезной балки из двух элементов постоянного прямоугольного сечения, определяемых пластическими моментами и . Требуется подобрать их значения, соответствующие минимуму массы конструкции в условиях упруго-пластической работы. Показана ориентировочная эпюра изгибающих моментов. Характерные сечения обозначены цифрами.

Решение

Массу конструкции можно принять пропорциональной предельным пластическим моментам и длинам стержней. Тогда целевая функция примет вид

,  (1)

где

,  (2)

Здесь b, h – размеры прямоугольного поперечного сечения.

Данная балка один раз статически неопределимая. Эквивалентная система метода сил имеет вид по рис.2, где лишняя неизвестная метода сил. Выразим через неё максимальные пролётные моменты

,  .  (3)

Условия прочности (пластичности) в опасных сечениях 1, 2, 3 имеют вид:

,  ,  ,  . (4)

Здесь условия прочности для сечения 3 выписаны дважды: для сечений слева и справа.

Для приведения к каноническому виду в среде Matlab двойные неравенства (4) перепишем в виде

,  ,  ,  ,  (5)

,  ,  ,    (6)

и в неравенствах (5) заменим выражениями из (3) и получим

,

.  (7)

В компьютерной программе на языке Matlab каноническую математическую модель данной задачи можно представлять в 2-х вариантах:

Целевая функция (1), равенства (3), неравенства (5), (6). Целевая функция (1), неравенства (6), (7).

Перепишем (3) в виде, совпадающем с каноническим

  (8)

Пусть параметры имеют следующие значения:

Компьютерная программа, реализующая 1 вариант математической модели, привела к результатам:

и графику для эпюры изгибающих моментов [x, fval] = linprog (f, A, b, Aeq, beq, lb, ub);

Пример 2. ( стр. 92) Заданы схема и две нагрузки для двухпролётной жёстко-пластической неразрезной балки из двух элементов постоянного прямоугольного сечения, определяемых пластическими моментами и (рис. 1). Требуется подобрать их значения, соответствующие минимуму массы конструкции в условиях упруго-пластической работы. Показана ориентировочная эпюра изгибающих моментов. Характерные сечения обозначены цифрами.