Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

2019

.pdf
Скачиваний:
2
Добавлен:
16.06.2024
Размер:
4.24 Mб
Скачать

Рис. 2.9. Оптимизируемая ферма

Эта задача представляется следующей математической моделью:

найти

min

 

Ni

 

li

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

F

 

 

0

 

при

 

kx

 

 

 

 

 

 

 

 

Fky

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Для данной задачи получаем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

min

 

N1

 

 

 

l1

 

 

 

N2

 

 

 

l2

 

 

 

N3

 

 

 

l3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

при

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N1 cos45 N3

cos45 P

 

N sin 45 N

2

N

3

 

sin 45 0

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(2.36)

(2.37)

(2.36*)

(2.37*)

Переходя к неотрицательным парам переменных, придем к виду:

 

min (x1 x2 ) 1,414 (x3 x4 ) (x5 x6 ) 1,414

 

 

при

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

0,707

0,707

0

0

0,707

0,707

 

x

 

10

 

 

 

 

 

 

 

 

3

 

 

0,707

0,707

1

1

0,707

0,707

 

x

0

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

5

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

6

 

 

Ответ:

(2.38)

(2.39)

x1 7,07; x6 7,07. 71

Оптимальная ферма изображена на рис. 2.10.

Рис. 2.10. Оптимальная ферма

Ni li 7,07 1,414 7,07 1,414 19,99.

Пример № 2. Найти статически определимую ферму оптимальной конфигурации. Исходная ферма представлена на рис. 2.11.

Рис. 2.11. Исходная схема фермы

Математическая модель задачи оптимизации:

 

 

при

найти

min

 

Ni

 

li

 

(2.36)

 

 

 

 

 

 

 

 

 

 

 

 

 

Ni cos i

 

 

k

 

 

 

 

Px

 

 

(2.37)

 

Ni sin i Pyk

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l1 1;

l2 1,414;

l3 1;

l4 1,414;

l5 1;

l6 1.

 

 

 

72

 

 

 

 

 

 

Уравнения равновесия узлов имеют вид:

 

 

 

 

 

 

 

N1 N2 cos45 10

 

 

 

 

 

 

 

 

 

 

 

 

 

N2 sin 45

N3

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N1 N4 cos45 0

 

 

 

 

 

(2.37**)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

4

sin 45

N

5

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

N

cos45

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

N

6

10

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После перехода к неотрицательным парам получим:

 

 

 

 

 

min f0 (x1 x2 ) 1,414(x3

x4 ) x5

x6

(2.36**)

 

 

 

 

 

1,414(x7 x8 ) x9 x10 x11 x12

 

 

 

 

 

 

 

Симплекс-таблица задачи:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

x2

x3

 

 

x4

 

x5

x6

 

x7

 

x8

x9

 

x10

x11

x12

 

 

–1

1

–0,7

 

0,7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

10

 

 

 

–0,7

 

0,7

 

–1

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

–1

1

 

 

 

 

 

 

 

 

 

 

0,7

–0,7

 

 

 

 

 

0

f(x)

 

 

 

 

 

 

 

 

 

 

 

 

–0,7

0,7

–1

 

1

 

 

20

 

 

–0,7

 

0,7

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

–1

1

10

Решение этой задачи:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x4 14,14;

x5 10,02;

 

 

x10 20,05.

 

 

2.10.Постановка задачи оптимизации балок

ирам на основе теории предельного равновесия

Будем считать материал балок и рам идеально жесткопластичным. При изгибе балки из такого материала при некоторой нагрузке в точке с наибольшим изгибающим моментом возникает пластический шарнир. Пластический шарнир – это точка оси балки, в которой возникает изгибающий момент, необходимый для полного перехода данного сечения в пластическую стадию. Этот постоянный по величине изгибающий момент называется «пластическим моментом сечения» М0. Его можно вычислить исходя из размеров сечения и предела текучести материала σТ:

M 0 Wпл т ,

где Wпл – пластический момент сопротивления. 73

Условия прочности стержня M1 Mk0 , то есть действующий изгибаю-

щий момент меньше или равен пластическому моменту сечения. Положим, что в конструкции имеется достаточное количество точек, в

которых могут одновременно образоваться пластические шарниры. Ввиду наличия пластических шарниров конструкция является статически определимой и ее прочность можно исследовать с помощью обычных линейных уравнений равновесия.

Пластические шарниры могут возникать в сечениях, примыкающих к узлам и местам защемления, а также в точках приложения сосредоточенных сил [18].

Уравнения равновесия можно подразделить на три группы: Уравнения равновесия узлов

q

Mik 0 ,

k 1

где q – количество сходящихся в i-й узел стержней.

Уравнения проекций внешних и внутренних сил на горизонтальную ось (для рам) (рис. 2.12):

 

m

 

 

 

S j Qtk 0 ,

(2.41)

 

j t

k 1

 

где t

номер яруса рамы;

 

 

τ –

количество стоек на t-м этаже;

 

m

количество этажей рамы.

 

 

Рис. 2.12. К уравнениям проекций на горизонтальную ось

74

Поперечные силы Qtk удобно выражать через концевые моменты. Уравнения равновесия ригелей (балок) (рис. 2.13):

M A (l un ) l Mn Mb (un l) Mnб,

(2.42)

где n – количество нагрузок на ригель.

При оптимизации объема материала рамы (стоимости материала) сделаем допущение о пропорциональности площади сечения стержня рамы предельному пластическому моменту этого стержня. Тогда объем (стоимость) материала рамы можно записать в виде

Z V M 0l

M 0l

2

... M 0l

(2.43)

1 1

2

n n

 

Рис. 2.13. К уравнениям равновесия балок (ригелей)

Рассмотрим численный пример.

2.11. Пример оптимизации статически неопределимой балки

Рассмотрим задачу синтеза двухпролетной балки (рис. 2.14), загруженной сосредоточенными нагрузками.

Рис. 2.14. Оптимизация сечений двухпролетной балки

75

Приняв в узле 2 одно расчетное сечение для изгибающего момента, можно исключить из рассмотрения уравнение равновесия узлов. Тогда система ограничений примет вид:

уравнения равновесия балок

0,5M

2

M

1

Pl

4;

(2.44)

 

 

 

 

 

 

 

1 1

 

 

0,5M2 M3

P2l2

4;

(2.45)

условия прочности

 

 

 

 

 

 

 

 

 

 

 

 

M1

 

 

 

M10 ,

 

(2.46)

 

 

 

 

 

 

M2

 

 

 

M10 ,

 

(2.47)

 

 

 

 

 

 

 

M2

 

 

 

M20 ,

 

(2.48)

 

 

 

 

 

 

 

 

M3

 

 

 

M20.

 

(2.49)

 

 

 

 

 

 

Целевая функция этой задачи

 

 

 

 

 

 

 

 

 

 

f (x) V 3M10 4M20.

(2.50)

Представим изгибающие моменты в сечениях парами положительных

чисел:

 

 

 

 

 

 

 

 

 

 

M1 x1 x2 ,

 

(2.51)

M2 x3 x4 ,

 

(2.52)

M3 x5 x6.

 

(2.53)

Неизвестные пластические моменты балок обозначим

 

 

 

 

x M 0 ,

 

(2.54)

 

 

 

 

7

 

1

 

 

 

 

 

x M 0.

 

(2.55)

 

 

 

 

8

 

2

 

 

Тогда задача оптимизации будет представлена в виде

 

min f (x) 3x7 4x8

(2.56)

при

 

 

 

 

 

 

 

 

 

 

0,5x3 0,5x4

x1 x2 7,5;

(2.57)

 

 

 

76

 

 

 

 

0,5x3 0,5x4 x5 x6 20;

(2.58)

x1 x2 x7 0,

(2.59)

x3

x4

x7

0,

(2.60)

x3

x4

x8

0,

(2.61)

x5

x6

x8

0.

(2.62)

Симплекс-таблица для этой задачи имеет вид: неизвестных – 8; уравнений-равенств – 2; ограничений – 6.

 

x1

x2

x3

 

x4

x5

 

x6

x7

x8

 

 

1

–1

0,5

 

–0,5

 

 

 

 

 

7,5

 

 

 

0,5

 

–0,5

1

 

–1

 

 

20

 

1

–1

 

 

 

 

 

 

–1

 

 

 

 

 

1

 

–1

 

 

 

–1

 

 

 

 

 

1

 

–1

 

 

 

 

–1

 

f(x)

 

 

 

 

 

1

 

–1

 

–1

 

 

 

 

 

 

 

 

 

3

4

 

Решение этой задачи:

 

 

 

 

 

 

 

 

 

 

 

M 0

5

x 5

 

 

 

 

 

 

 

1

 

 

7

 

 

 

 

 

 

 

 

M 0

17,5

x 17,5

 

 

 

 

 

 

2

 

 

8

 

 

 

 

 

Эпюра моментов в предельном состоянии приведена на рис. 2.14.

M1 x1 5;

M2 x3 5;

M3 x5 17,5;

 

x2 x4 x6 0.

2.12. Пример оптимизации рамы

Конфигурация рамы и схема нагрузки указаны на рис. 2.15.

Стойки будем проектировать одинакового сечения. Тогда объем материала рамы

V 2 4M10 6M20 ,

(2.63)

77

 

или

f (x) V 8M10 6M20.

(2.63*)

Рис. 2.15. Схема оптимизируемой портальной рамы

Составим уравнения равновесия отдельных частей рамы: а) уравнения равновесия узлов

M2 M3 0,

M5 M6 0;

б) уравнение проекций на горизонтальную ось

M1 M2 M6 M7 10; 4 4

в) уравнение равновесия для ригеля

0,5M3 M4 0,5M5 30.

Условия прочности стержней:

M1

 

 

M10

или

M10 M1 M10 ,

 

M2

 

 

M20

или

M10

M2

M10 ,

 

 

M3

 

 

M20

или

M20

M3 M20 ,

 

 

M4

 

 

M20

или

M20

M4

M20 ,

 

 

78

(2.64)

(2.65)

(2.66)

(2.67)

(2.68–69)

(2.70–71)

(2.72–73)

(2.74–75)

M5

 

 

 

M20

или

M20

M5

M20 ,

(2.76–77)

 

M6

 

 

 

M10

или

M10

M6

M10 ,

(2.78–79)

 

 

M7

 

 

 

M10

или

M10

M7

M10.

(2.80–81)

 

 

 

Количество неравенств можно в два раза уменьшить, добавив в каждое соответствующее Mi0 . Из (2.68–2.81) получим:

 

 

M1 M10 2M10

 

(2.82)

 

 

M2 M10 2M10

 

(2.83)

 

 

M3 M20 2M20

 

(2.84)

 

 

M4 M20 2M20

 

(2.85)

 

 

M5 M20 2M20

 

(2.86)

 

 

M6 M10 2M10

 

(2.87)

 

 

M7 M10 2M10.

 

(2.88)

Введем обозначение

 

 

 

 

 

 

 

x M

i

M 0 ,

i 1,2,...,n;

k 1,2,...,m.

(2.89)

i

k

 

 

 

 

 

 

 

x

 

M 0 ,

 

(2.90)

 

 

n k

 

k

 

 

где n – количество сечений;

 

 

 

 

 

 

m – количество различных сечений.

 

 

 

Тогда задача линейного программирования примет вид

 

 

 

min f (x) 8x8

6x9

 

(2.91)

при

 

 

 

 

 

 

 

 

 

x2 x3 x8 x9 0

 

(2.92)

 

 

x5 x6 x8 x9 0

 

(2.93)

 

 

x1 x2 x6 x7 40

(2.94)

 

 

79

 

 

 

 

 

 

 

 

x3 2x4 x5

60

 

 

(2.95)

 

 

 

 

 

x1 2x8 0

 

 

(2.96)

 

 

 

 

 

x2 2x8 0

 

 

(2.97)

 

 

 

 

 

x3 2x9 0

 

 

(2.98)

 

 

 

 

 

x4 2x9 0

 

 

(2.99)

 

 

 

 

 

x5 2x9 0

 

 

(2.100)

 

 

 

 

 

x6 2x8 0

 

 

(2.101)

 

 

 

 

 

x7 2x8 0.

 

 

(2.102)

Данная задача заносится в следующую симплекс-таблицу:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x1

x2

x3

x4

 

x5

 

x6

x7

x8

x9

 

 

 

1

–1

 

 

 

 

 

 

–1

1

 

 

 

 

 

 

 

1

 

–1

 

–1

1

 

 

–1

1

 

 

 

 

 

–1

1

 

 

40

 

 

 

–1

2

 

–1

 

 

 

 

 

60

 

1

 

 

 

 

 

 

 

 

–2

 

 

 

 

1

 

 

 

 

 

 

 

–2

 

 

 

 

 

1

 

 

 

 

 

 

 

–2

 

 

 

 

 

1

 

 

 

 

 

 

–2

 

 

 

 

 

 

 

1

 

 

 

 

–2

 

 

 

 

 

 

 

 

 

1

 

–2

 

 

f(x)

 

 

 

 

 

 

 

 

1

–2

 

 

 

 

 

 

 

 

 

 

 

8

6

 

Неизвестных – 9. Ограничений – 11. Равенств – 4.

В незаполненных клетках таблицы – нули. Решение данной задачи дает:

x1 0;

x2 6,666; x3

6,666;

x4

33,333;

x5 0;

x6 0;

 

x7 33,333;

x8 16,666;

x9 16,666.

 

Итак, сечение рамы: M10 16,666;

M20

16,666 .

 

 

 

 

80

 

 

 

 

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]