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

книги / Оптимизация подкрепленных цилиндрических оболочек

..pdf
Скачиваний:
0
Добавлен:
12.11.2023
Размер:
10.25 Mб
Скачать

 

 

 

 

 

 

 

 

 

 

 

 

А,

Номер проекта

кН

Р

Т ,« т

еэ--ю3 кн

3,10*

кН ££-103кН.

 

<£ю?кн

ё .н

в ;

6+

46

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Н

Н

б',%

 

 

 

 

 

 

 

 

 

 

1

4

0,963

36

3,25

3,34

3,04

3,66

2,65

324,3

291

348

0,94

2

6

0,95

45

3,66

3,73

3,36

4,04

3,2

347,4

318

371

1,42

3

6

0,9999

45

4,85

5,06

4,86

5,15

15

402,3

374

422

1/4

4

4

0,95

20

2,97

3,15

2,80

3,37

2,90

319,4

290

339

3,9

5

4

0,9999

20

з,в

3,92

3,64

3,98

3,64

326,1

315

334

1,2

 

V®1 V "1

 

 

Р

 

 

Р *

Р "

6 '

GCT.

Номер проекта

 

 

 

ст

ч

ч

 

м

М

 

м

 

 

Р*

 

н

Н

 

 

 

 

 

 

 

 

м

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

0,9

9,0

0,55.

8,0

20

1.0

0,96

0,94

0,94

0,968

334

324,1

2

0.5

9J5

0,6

8,25

26

0,978

0,76

0,78

0,8

0.894

353,8

351,1

3

1,0

1,5

0,75

9,0

8

1.0

1.0

1,0

1,0

1,0

409,2

406,7

4

0,75

9,0

0,55

0,75

26

0,997

0,92

0&4

0,92

0,92

326

323,8

5

0.7

* £

0,55

9,0

20

1.0

0,96

0,98

0,90

0,992

327,2

326,6

дД ,

д а ,

еЧ

О

о:% И

6,2

12,3

321,2

*2

13.1

342,1

0

0

394,9

5,7

2,4

307,1

4,2

4,2

322

 

Таблице 18

46

4 с г*

Q X

 

гст

в°,%

хЮ3,

х ю *

 

кН

кК

0,94

3,34

15

2,6

3,73

3,87

3.05,15 15

5,5

3,18

3,26

1,4

3,97

15

ба

= 6 . 1 0 2 кН, Р = 0,95, Т = 45 лет; 3 = б0 « 6 1 0 2

кН, Р '

II

0,9999, Т = 45 лет; 4 - 6Q = 4 • 1 0 2 кН, Р = 0,95, 7 = 2 0

лет; 5 -

о

д

- 4 • 1 0 2 кН, Р = 0,9999, Г = 2 0 лет. Результаты исследований

 

показывают, что распределения весов оболочек близки к нормаль* ному закону, при уменьшении дисперсии нагрузки, а также уве* личеиии требований к оценкам Р и Т средние квадратические

отклонения уменьшаются. С увеличением требований к надежно­ сти и сроку эксплуатации {проекты 3,5) влилние неточности изго­ товления параметров оптимальных оболочек уменьшается (рас­

пределения 4 и Q& ): разброс значений Qi и

сокращается.

На рисунке указаны характеристики. Q < и Аг оптимальных обо­

лочек, обеспечивающие заданную надежность Р

и долговечность

Т- при изменяющейся нагрузке. Приводится характерный вид кри­ вых распределения Of, а также упорядоченных значений 0^/ гДе п - номер варианта случайной реализации проекта £ . При вы­ соких требованиях к Р (проекты 3,5) ограничение Q% выполнено (0 ft > 15000 кН). На рис. 23 цифрами 1,2,' 4 /- отмечены гранич­ ные значения критических нагрузок^обеспечивающие надежность

оболочек; при малых отклонениях

область экстремума „раз­

мыта”, содержит близкие значения Q 2.

 

В табл. 17 обозначено^ Q 9 значения Q\ и Q 2 »обеспечива­

ющие заданные Р и Т;

0_- математическое ожидание веса при

решении задачи (2.40);

G~t G* - -область изменения и G0 - вес

оптимальной оболочки в постановке модели нелинейного про-

граммированил (2.52); Q1

-

математическое ожидание ограниче­

ния прочности, a Q^t Q^

-

область изменения значений; Q о -

нижняя граница Q* при случайных реализациях

,* AG = о

- е ° ; д а , - а * - а ; , - д а г = а а - а - .

 

 

 

Табл. 18 содержит, оптимальные проекты оболочек с учетом

требований сортамента {ftc, Нс,

к}, а также их вес

^

и надежность pQji

 

® i

^ 1 ст> ® 2 ст " пара­

метры Qj и Qs при стандартизации

компонент

вектора Х0^

В табл. 18 указаны вероятности: выполнения ограничений прочней оти PQ , критических нагрузок потери устойчивости PQ 2, обес­

печения надежности

Р -

наименьшая оценка надежности

при случайных реализациях

 

Показатель веса решения задачи

в постановке (2.39) для значения

Р = 0 ,8 обозначен через G

Исследования показывают, что все проекты, полученные в рам­ ках НЛП или ЗСП подходов с последующей заменой оптимальных параметров значениями из сортамента, являются допустимыми; здесь относительные потери в весе менее € %. С убыванием дис­ персии нагрузки GQ разброс значений веса в и характеристик

и Q<i также уменьшается (рис. 23). Для всех рассмотренных случаев установлена близость оптимальных проектов оболочек, полученных для моделей нелинейного и стохастического програм­ мирования и с учетом дискретности сортамента, причем с возра­ станием оценок Р и Т влияние отклонений от параметров опти-

йс,Ьш-

мальной конструкции Х ° убывает и подход с позиций нелинейного программирования с последующей стандартизацией Х° становит­ ся все более оправданным.

Отметим, что с расширением условий неопределенности на этапе постановки задачи ОПК эффективность использования ком­ промисса в виде (2.52) возрастала. Количественным выражением этого факта явилось уменьшение параметра Ы= ( Fc - Г и ) / F

в области оценок характеристик FC[F F w ]; Fi - значения кри­ териев, найденные из решения задач оптимизации функции (2.52).

3 .1 0 .3 . Оптимизация цилиндрических оболочек многоцелевого назначения

Рассмотрим задачу проектирования цилиндрических оболочек длины L и радиуса У, подкрепленных продольными и поперечны­ ми ребрами, при условии, что конструкции испытывают раздель­ ное воздействие осевых сил Р^ и радиального давления Рг в за­ висимости от возможных условий эксплуатации. Здесь парамет­ ры оболочек, оптимальных для класса названных воздействий, по­ лучаются как результат некоторого компромисса и могут быть найдены в рамках многокритериального подхода, изложенного в гл. 2. Считал, что управляемыми параметрами являются Н 0 - толщина обшивки, Ьй, Ьш- толщины, высоты стрингеров и шпангоутов соответственно, а также к, к А- числа ребер каждо­ го направления, сформулируем многокритериальные задачи не­ линейного программирования согласно 2.1: построим функции критериев и ограничений (геометрических, прочностных, крити­ ческих нагрузок потери устойчивости конструкции и элементов).

Функция веса оболочки имеет вид (2.41). Выражения для оп­ ределения критических нагрузок .общей' потери устойчивости оболочки получим в рамках энергетического метода, следуя [4]. А именно, принимаются допущения, принятые в 3.10.1, учитывают­ ся дискретный характер расположения подкрепляющих элементов без эксцентриситета и продольные деформации ребер, связанные с деформацией оболочки, используется одночленная аппроксима­ ция прогибов оболочки W (как и в [2]). Потенциальная энергия си­ стемы складывается из энергии деформации обшивки, энергии, накопленной при деформировании шпангоутов и стрингеров, за вычетом работы внешних сил. Учитывая, что оболочка находится в двухосном, а ребра - в одноосном напряженном состоянии, причем продольные деформации ребер и оболочки совпадают, по­ лучаем для.обоих случаев нагружения PN и Рг выражения для

критических

усилий

общей

потери

устойчивости конструкции

в виде

 

 

 

 

 

 

 

 

?л “ [ A

t f W

/ • > ( » “ r S » N с ш )

Фдо

^ '0

Мс

у

'2С

у

''1 iu

X

' 2 ш Х '

 

 

 

+

 

 

С ^ ] А и 2 Г уm) :

 

 

 

 

(3.12b;

 

 

 

 

 

 

 

 

 

 

\

=[

 

 

 

* 72cr< , X

С

гхш+ y L rx ) + &

ry +

 

 

 

 

 

 

 

+*JLr* V « +*e

 

 

 

 

 

 

 

Здесь обозначено:

, Vp

 

безразмерные характеристики

напряжений,

Г

^

^

/ l x

^ ' ) ,

?^ - / >

р/( х ,.£ ,) ;

 

 

 

0 г Ьш/1 ;

<^т п =('/яг + лг) |//п 2л

FCl

Fo5-^

- площади стрингера

и обшивки соответственно; г с= 2

соз'яА ;

= Z sin* /*£ : с

 

 

 

 

 

 

 

 

 

 

 

 

 

* L

У

ЫЛ

 

1

 

га шД

C0S* M

;

^ = . ^ ш

г тсС( -

f =

f k ;

e = 8 ( V U

;

X

' ^

o

i

Ёщ~Еш/Е0; t * x t i

v "

 

 

 

 

 

 

 

 

?K * m4

* «

/t’Tl'-hofli

<£ - « Ч

1 ^

1 !r(','lo)l (1 + t y li

 

V ,'l“

2<r t,- ^

4

C

/( * 5

L r/n 2) ; ^

= n %

I “p/ [ h

‘ lr ( 1 ♦

* 4 й ;

 

- " ^ O

n V - O

 

] ;

i p = t m 2[ ( « * - t ) t 4 r '

+ t V

^ n/(;H (f-^}(n2- i;].; m=Bir<(i/Tig ■

 

(3.124J

£

-

^

4

£

A

*

^

V

H

K

I ^

) ]

;

f c

i c Fc%

](* b l% ,J ;

C

 

 

 

 

 

Л Л

) ;

 

 

 

 

 

 

фт п \,

 

C

= W * « * £ „ , ! £ / [ £ Л

» ! - 1K 1

 

) ];

Ф

^

' Л

8)* ;

 

£

3

f c Fc m * <V

l >

' - V

 

- 1

) Фт

] ;

V

s(' " l '

y rt*J‘ i

ЛИ « Ч ^ „ 1 ;

величины, имеющие обозначение „с ” , относятся к стрингерам, а „ш " - к шпангоутам; „кр ” соответствует кручению, а „и з" - из­ гибу; *7?0- число полуволн по образующей.

При осевом сжатии возможны следующие частные случаи де­ формации оболочки [2] (для £ш = £0 расстояние шпангоута от торца равно расстоянию между шпангоутами): 0 ,1 ,2 ,4> 5 ,8 . При­ чем второй и восьмой случаи реализуются при неосесимметрич­

ной деформации, когда аппроксимация прогибов представима:

W = /sin (m o < )sin (n jfyтогда в (3.123) необходимо заменить Гу на Гуш и наоборот. Согласно [2], происходит следующая деформа­ ция ребер: О - ребра одновременно изгибаются и закручиваются;

1 -

стрингеры изгибаются; 2 - стрингеры только закручиваются;

4 -

шпангоуты закручиваются; 5 - стрингеры изгибаются, шпан­

гоуты закручиваются; 8 - стрингеры и шпангоуты только закру­ чиваются.

При нагружении равномерным давлением учтем следующие частные случаи деформации конструкции:

.

-1 .

Р *

Р

Р,

V 1 р

*г1

0)Vp = ,т А

+г 1\

+!гсИ“ Г Г1ш+-Т- t j *

м>

 

 

К .-И

 

'

 

 

 

 

 

£

' ЗШ

'

 

-I

1 г Р l Р Л +1 р

V I Р V. 1, Р .

р<1)

 

 

 

+- h

 

V T

 

 

 

 

 

 

 

,

Г р ,

Р

 

 

P i

Л и Р п .

Ч?р. =

+ М г и + г Л с + 7,г^2с/)+ 72

Зш

44J

 

 

 

 

 

 

 

гЬ Г

U -

(3.125)

Если критические усилия выступают в качестве критериев за­ дачи оптимального проектирования, то используются выражения

для у

к , причем берется минимум no (nt0,n) и видам де-

РЫ

\ г

формации (о). Ограничения по критическим напряжениям общей потери устойчивости оболочки определяются соотношениями

diQ{i) * $ (™

, С ц » " V

0;

Qe{PH ' pJ -

(3 ,гб)

в которых о = {о, 1,2,4, б, 8 } для Q = Я .

ио = { 0 , 1 ,4 ,5 } для Q =

= рг } GQ=(jQl{Et)\

пРи

Pw

6Q = N/(2wrh0+ * £ ) ;

при

Q * ^ X = P r l ( \ l r + h < F m / l L r ) l

Кроме общей необходимо также исключить возможность мест­ ной потери устойчивости элементов конструкции - шпангоутов, стрингеров. Следуя [21]г считаем расчетной моделью стрингера пластинку, свободную по одной и шарнирно опертую по трем сто­ ронам, тогда ограничение значения критического напряжения по­ тери устойчивости

а,оэг

 

(3.127)

(t-v>c)

*

N

Ограничения по максимальным напряжениям для шпангоутов при осевом сжатии и стрингеров при радиальном давлении можно

получить из соотношений (3.127) соответствующей заменой пра­ вой части на , Q е {Рц, ?г }, и величины х* на x j,

Предполагаем, что выполнение других ограничений местной потери устойчивости элементов конструкции обеспечивается за очет надлежащего выбора на основании исследований, выполнен­ ных в [2 1 ], геометрических характеристик подкреплений и других искомых параметров.

Ограничение по прочности на торцах оболочки при Q = PN имеет вид

?3e(i)=PN/(2^v'<Fc)<teA], »к»>

где [бд}~ допустимое напряжение. Формируя ограничение проч­ ности при Q = РГ1 учтем, что радиальная нагрузка Рг при отсут­ ствии изгиба оболочки эквивалентна действию сжимающих уси­ лий Ру =РГ Г [ 2 ; 4], тогда

9*6т =pr r / ( V

° - 129->

ГД* Ь = ЬШЛШ ^ .

Задание геометрических ограничений искомых параметров

£ 7

$ x L $

i = 1,7; х 7 щх 7 - const,

(3.130)

завершает формирование модели оптимального проектирования многофункциональных ребристых цилиндрических оболочек (2.41), (3.123) - (3.130). Для анализа поставленной частично целочислен­ ной задачи НЛП использовался метод случайного поиска, описан­ ный в гл. 2. При минимизации (3.126) по (/п0,л) использовался перебор по сетке с адаптацией области дт п ,как указано в 2.4.

Для корректности модели задачи оптимального проектирова­ ния вводились геометрические ограничения: а) „сохранение'* ре­ бер; б) „сохранение” допустимых соотношений высоты hp и тол­ щины ftp элементов подкрепления [2 ] в виде

а) К РЭ

s d 3 ,

(3.131)

где dv d2, d3

- заданные постоянные, а индекс

ре{с,ш }. При

проектировании d ^ выбирались равными d ^ = 2 ; ц 2 = 0 ,0 2 ;

d 3 = 0,1.

В табл. 19, 20 приводятся результаты оптимизации ребристых оболочек длиной L = 0,8 м, полученные методами случайного по­

иска, когда Е = 0,686 1 0 е кН/м

; Р = 0,3; [0ГД ] = 0^27 * 1 0 6 кН/м .

В.этих таблицах ojoaHaneHo:

Ц)=г!Ц Q^ (X1 0 Н) - осевая

нагрузка; Qp (* 10 Па) г нагрузка радиального давления; V 1 2 мэ) - объем оболочек; 4 C^t* 1 0 кН/м2) - напряжения на торцах оболочки при осевой нагрузке ( 3 . 1 2 8 ) ; ( * 106 кН/м2) - напряжения (3.129) при радиальном давлении; oQp (* 104 кН/м2) - критическая нагрузка потери устойчивости оболочки при радиаль­

ном давлении (3.125); ограничение (3.126) для случая осевого сжа­

тия QN было неактивным;

критические напряжения местной

потери устойчивости: дс (* 10

кН /ма) - стрингеров (3.127); д ш

(Х10* кН/м2) - шпангоутов .(3.127). В табл. 19, 20 указаны также

параметры волнообразования ( т н -

число полуволн по образу­

ющей при Оц, а для Qp trip - 1,

rip - число волн в окружном

направлении при 0 ^

и Q „ соответственно) и номера форм поте­

ри устойчивости

для 0 ^ и Jp (3.125) для Qp ). Характеристи­

ки форм потери устойчивости указываются независимо от того, какой вид нагрузки прикладывался и какое ограничение было ли­ митирующим.

Таблица 19

Номер проекта

 

N

° Р

 

9QN

ч 9Q1

9 Ш

 

 

У "

Q

 

У- ю“2

 

 

 

 

 

1

 

0,6

167

0.5640

0,27

0,242

 

0,273

0,085

2

 

0,5

-

8,8

0,5633

0,44

0,12

8,83

0,114

0,059

3

 

0,6

60

2,95

0,5632

0,19

0,004

2,95

0,20

0,15

4

 

0,76

ЗОВ

-

6,923

0,27

0,205

5,1

0,3.

0,084

5

 

0,76

-

5

0,923

0,204

0,01

0,165

0,073

6

 

0,76

92,4

1,9

0,920

0,27

0,003

1,91

0,278

0,12

7

 

0,5

126

6.5

0,601

0,245

0,006

6,68

0,25

0,267

В

 

0,6

200

10

0,963

0,269

0,081

11

0,270

0,18

9

'

0,6

-

5

0,56В

0,28

1 0,007

5,1

0,14

0,011

10

 

0,6

5

0,470

0,40

0,009

5,14

0,27

0,012

11

 

0,25

60

0,345

0,45

0,003

66,5

0,124

0,035

12

 

0,25

30

0,354

0,91

0,0016

зо;в

0,1

0,02

12В

Количество стрингеров при радиальном давлении и для мно­ гофункциональных конструкций определялось, как правило, ниж­ ним’ допустимым значением к * = кт -1П = 4. Оптимальные обо­ лочки разыскивались при следующих ограничениях на число ре­ бер: 4 ^ к ^ 50; 2 4 4 40. Число шпангоутов для QN также было минимально допустимым: ft*= k itn^n = Z

Табл. 19, 20 содержат и параметры оптимальных многоцеле­ вых (многофункциональных) оболочек (проекты 1-6) и оболочек, предназначенных для учета неопределенности условий функцио­ нирования и оценки несущей способности (проекты 7 ,8), а также характеристики локальных и глобальных экстремумов для одно­ критериальных задач. Многоэкстремальность решения связана с активностью ограничения (3.126) критических сил потери устой­ чивости для Qp (проекты 9-12).

В задачах (2.43)-(2.45) требовалось найти параметры много­ функциональной оболочки заданного объема (веса) / 0 , выдержи­

вающей наибольшие нагрузки

Qy

и

Qpf удовлетворяющие кри­

терию (2.44). При заданном объеме

 

= 0,564 • 102 м3 оптималь­

ные одноцелевые оболочки длиной L = 0,8 м и ф = 0,5 выдержива­

ют максимальные осевые нагрузки

Од, = 1,87 • 103 кН (проект 1)

и радиальное давление Q р = 0,88 • 102 кН/м2 (проект 2). Много­

функциональная оболочка того же объема

t

оптимальная по

(2.44), имеет критические нагрузки:

 

= 0 ,6 -103 кН и

Qp =

= 0,295 кН /м 2^(проект 3),

при

этом

значение

максимина

(2.44)

С^1 = 0,32. При ф = 0,75 и

 

= 0,923 • 10~2 м3 соответствующие

характеристики оптимальных оболочек были следующие:

=

= 0,308 *104 кН (проект 4);

Qp

= 0,5 - 102 кН/м2;

= 0,924 х

х i0 3 Н; Qp = 1,9 • 104 Па}, «4 = 0,3 (проект 6).

 

1

 

При решении задачи оптимизации оболочки в постановке (2.46), как многокритериальной задачи выбора запаса несущей способ­ ности в условиях неопределенности, задавался начальный объем

материала V = 0,8 10“ 2 м 3; ф - 0,5;

L= 0,8 м. При этом много­

целевая

оболочка имела

критические

нагрузки {Q ^

= 125 х

^

Цръ = 6,25-10

 

 

^

х 10 Н;

Па (проект 7), В предположении о воз­

можности отклонения амплитуд нагрузок допускается увеличение объема конструкции на 20 %. Тогда оптимальная по (2.46) оболоч­

ка при

=

А 9 = .0,5 имеет следующие характеристики: Q°

=

= 2 - 1 0 3 кН;

Qp

= 10г кН/м2 (проект 8). Увеличение объема кон­

струкции

на

2и%

позволило поднять критические нагрузки

на

60 %. Как и ранее, для ф = 0,5 в оптимальных решениях (проек­

ты 7,8) активны ограничения прочности

(3.128) и критических

нагрузок устойчивости (3.126).

 

Проекты 9,10 и 11,12 табл. 20 показывают многоэкстремальный характер задачи оптимального проектирования. Для проектов 9,10 ограничение (3.126) для Qp выполняется как равенство ( Qp ~ ) для различных по объему вариантов параметров оболочек;

в зо>)е активности ограничения устойчивости (3.126) критические нагрузки потери устойчивости могут-быть увеличены без увеличе­ ния объема оболочки (проект. 11 - Qp * QQ = 0,60 ■1(г кН/м*,

V = 0,346• 10* м3; проект 12 г Qp « q'Q

=*0,3 • 103 кН/м2, V* =

= 0,354- to-3 MV

1

Остановимся на вопросе выбора эапаоа несущей способности многоцелевой оболочки посредством критерия (2.46) в условиях неполной информации. Сравним результаты оптимального по (2.46) проектирования для различных значений коэффициентов

выделения ресуроов С; 0> = Сг Gopt и Мд С1 = = 1,2 (проект в. 0opt соответствует проекту 7, табл. 17) Здесь

( С г

= 1,1) получено: Q®1 = 1,5* 103кН; Qp1 = 0,76-103 кН/м3;

^ ° 1 = (Р/ “ Рз* У рз1 * °Я Ч> =

*о = 0,224-10-* м; Л с =

= 0,463 * 10“ * м; Ьш= 0,623 м; &с = 0,463 • 10“ 3 м;

Ьщ= 0,385 х

х 10**м; к = 6;

9; т - 1; п - 3; Ты ~ 0. Для проекта

8 ^ =

= 0,6; при вычислении d принималось ft

= Qw = 1,26 •

103 кН;

Ра

= 0Ы - 0,75 • 10 кН/м j соответствующие относительные

значения оС'^°/А6 = ^°/(С~1) равны о£0 = 3,

= 2. По­

скольку увеличение объема многоцелевой оболочки вызвано рас­ ширением условий неопределенности, обусловленной возможны­

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

f t ,, то

,0

.Oi

7 О

T QI

 

L

пол^ение«

> d 1

и «

^

свидетельствует о повышении

эффективности многокритериального выбора на основе (2.46) при расширении условий неопределенности (С, > С2),

Соседние файлы в папке книги