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

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

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

Далее рассматриваются способы реализации наиболее трудо емкого этапа - равномерного заполнения Рр - указанной проце­ дуры. Заметим, что модификации алгоритма, приведенного на рис. 4, могут различаться методами аппроксимации Рр по узло­

вым точкам, начиная с простейших линейных моделей; расчет па-

•—ijr\

раметров «Г0^

осуществляется сверткой типа (2.49) в выделен­

ной области Рс

при Ль- А., к, f,£ {1: 1_} й'нормализации вида

(2.47).

1

Равномерное заполнение области компромиссных решений узловыми точками. Рассмотрим возможность использования то­ чек, равномерно заполняющих гиперкуб в пространстве перемен­ ных проектирования [33], для равномерного размещения узлов аппроксимации в области Парето Рр (2.32). Основанный на Я П ^- поиске метод многокритериальной оптимизации [33] может ока­ заться неэффективным, если расчет конструкции для каждой точки Xq. требует значительного ресурса ЭВМ (что отличает задачи про­ ектирования строительных и других конструкций). Только малая часть точек ЛПГ - последовательностей входит в (2.32), причем

отбор F (£]e

Рр производится после генерации заданного числа

векторов

и расчетов конструкции.

Вопрос о равномерном заполнении Рр состоит в том, во-пер­ вых, как перевести свойства ЛП^-точек из Еп в Рр, во-вторых, как эффективно реализовать на ЭВМ класс соответствующих экст­

ремальных задач.

_

^

Пусть t -число критериев

 

, t , a FJ G НО; 1} - их

нормализации. Отметим два эвристических приема „равномерно­ го” заполнения Рр узловыми точками.

П о с т а н о в к а 1. В узлы располагают из условия равен­ ства ( t - 1 )-го критериального значения F * соответствующим ко­

ординатам Ху,т *€ [0; 1]^ точек Л -последовательностей

[33] для

( t - 1 )-мерного куба; ^

разыскивают из решения задачи вида

Ш . Найти: m ax {<?> ( я ) =

t f ) }

50)

 

 

у о 1 .2 ,-...,

 

П о с т а н о в к а критериев Fj , j

2 . В узловых точкам Рр отношение значений- t f равно отношению соответствующих ко­

ординат

точек ЛПт -последовательностей в t -мерном кубе.

Размещение в (2.32) узлов с указанным свойством осуществляет­ ся путем решения задач вида [14)

П2: Найти m ax { Ф ( F, <5) = m in ^ ( X ) /d j } .

<г 5 , >

Относительно коэффициентов оЛ для постановок П2 реализуются

требования

F ^jF ^.

Значения коэффициентов o(j находят

из условий нормировки o (j -

X j* * f 2 X *

Алгоритмы совместной реализации множеств экстремальных задач. Повышение эффективности решения наборов экстремаль­ ных задач (2.50), (2.51) размещения „узлов” в Рр опирается на принцип их совместной максимизации: „приписывания” группы задач 9 i вида (2.50) - (2.51) к одной „траектории подъема” (в хо­

де оптимизации одних задач Ф ^ j вводятся в область

решений,

и некоторые другие Ф^5р . Вблизи JPp траектории задач

рас­

щепляются, и точка расщепления служит стартовой для процессов максимизации % .

Алгоритм „группового подъема” составляют этапы: 1) задание

количества г „узлов” Рр

и начальной точки £

Чг» 2) форми­

рование Г задач

вида (2.50) или £2.51) и выставление началь­

ных точек их траектории подъема ( J ^ r j€ Ъх )

 

3) назна­

чение модальной (ведущей) задачи Фг в

группе; 4)

проверка

допустимости Х ^ ^ Ъ у . ,

в случае невыполнения условий пере­

ход к выбору новой точки

(п. 7); 5) вычисление в очеред­

ной точке поиска

 

значений целевых функций ^

 

и

критериев (2.50), (2,51); 6) отбор задач Ф ,

для которых

оценки

(2.50), (2.51) улучшены:

 

s - i) »

формирование

новой группы задач Ф зд , замена для задач

Ф ф

точек траекто­

рий подъема на 2

f »;

ч

если для ведущей задачи Ф>

точка

1 5ггг н е обеспечивает

увеличения оценок

(2.50)

или

(2.51), то

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

и выработки

новой точки X s+ 1{ r } . в рамках принятого метода оптимизации; если задача Фг решена, то перейти к п. 9; 8) перейти к п. 3); 9)

распечатка и запоминание параметров решения ФГ^; исключение задачй Фг из списка { г } ; если список { г } исчерпан, перейти

к п. 10, в противном случае разыскивается траектория с наиболь­ шим числом задач Ф ^ ) и перехоходят к п. 3); 10) конец.

Варианты алгоритма могут различаться как способом выбора модальной задачи в группе г , (п. 3), так и методом оптимизации (п. 7). Индекс г , может выбираться: произвольно; по назначению

ЛПР; из у с л

о в и я |

Последнее

способствует

либо сохранению группы

на очередном шаге,

либо окончанию решения задачи Ф .

 

 

• i

 

Расчеты на ЭВМ свидетельствуют о достаточной эффективности постановок (2.50) и (2.51) при решении наборбв этих „узловых" за*

задач алгоритмом „группового подъема". Процедура проектиро­ вания, представленная на рис 4, использована в 3.10.

Аксиоматизированные методы решения многокритериальных задач проектирования конструкций. Применимость аксиоматиче­ ского подхода к выбору принципа оптимальности в некоторых классах МКЗ особенно ценна, так как позволяет прогнозировать характеристики компромиссных решений. За счет появления дополнительных средств априорного анализа МКЗ удается наибо­ лее точно назначить параметры известного принципа оптимально­ сти {Р у Ш }. Данный подход не заменяет, а точнее раскрывает смысл содержательного анализа МКЗ, рекомендуя рациональную форму обобщения информации о задаче.

В работе [26] обоснована возможность аксиоматизированного выбора игрового принципа компромисса вида

,

-1

min

f~ -

FjW ~ FJ

1

opt { f

< « } = rnai

i F. (I ) =

 

 

f

 

XeD

ieCislj]

*

F.

- F.

J

 

1

7

 

J

1

(2.52)

а МКЗ проектирования конструкций; приводятся признаки, классы соответствующих задач ОПК и разработаны человеко-машинные

процедуры, использующие „механизмы" прогноза оценок^ { Fj } В (2.52) обозначено / } * * = m ag Fj[X), a F * = m in Fj (Х р );р Ф j.

При этом выбор компромисса (£.52) связывается не только с ак­ сиоматикой (в данном случае с единственностью решения , его оптимальностью по Парето, независимостью от масштабов измеоения /}(£), симметричностью относительно je [1; C j ] : FL = - F j, L ,i £ f 1i l j ) , но и с выявлением физических, природных, •свойств исходной проблемы проектирования, обеспечивающих реализуемость соответствующего решения.

Остановимся на формальном выделении класса многокритери­ альных задач ОПК, для которого возможно аксиоматическое обос­ нование компромиссного критерия вида (2.52). Назовем такие за­ дачи МКЗ с совокупно противоречивыми критериями (ЗСПК). Пусть

 

 

 

 

Fj{k)U CfkU F. ( * ) * * .

Ket/{p, q};

je C irk ]}*^ (2.53)

Xs

 

P (Ft

- вектор критериальных

где Qj - знак пустого множества, a. с

 

ограничений. Критерии

F p (X ),F ^ff) являются противоречивыми

(компромиссными), если при заданном

 

- таком, что решение

задачи вида

=

Fp (X)

 

 

F

 

 

XeFp (X)

(2 54)

существуй г, ограничение дин значений

Fq(X) обращаете/* в ра­

венство: Fq {£ р ) «

Cq

 

 

О п р е д е л е н и е .

Критерии { / 7

f совокупно противоречи­

вы, если для любых /э,q e t

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

сти (2.64). Иначе говоря, в

FPqlX ) имеют место соотношения типа

когда

(Рро(*) являются

некоторыми монотонно

убывающими функциями.

Доказательство существования и единственности решений ЗСПК, соответствующих моделям (2.53), (2.54) и найденных по­ средством компромисса (2.52), имеет идейную близость с работой

[14]. При этом решение $ Q^ является оптимальным

по

Парето,

в также симметрично: ^ c(JT0p ^ ) = F.^ №cpf),

Щ '

Из вы­

ражения (2.47) для нормализации частных критериев и равенств

Fp~ Fq получается условие инвариантности решения к ли­

нейным преобразованиям Fj (£)•

Помимо содержательных классов ЗСПК (например, оптималь­ ное проектирование многофункциональных ребристых оболочек в условиях неопределенности: МКЗ с критериями „вес” , „надеж­ ность” , „долговечность” и др. [26]) можно указать и формальные признаки характеристик задач ОПК, удовлетворяющих отношениям совокупной противоречивости (2.53), (2.54). Связи, требуемые для ЗСПК, обнаруживаются между значениями целевой функции и ак­ тивными ограничениями в однокритериальных задачах ОПК. Пере­ вод величин, представленных этими ограничениями, в состав ком­ понентов векторного критерия (2.2) приводит к модели ЗСПК, Та­ ким образом, совокупная противоречивость некоторых характе­ ристик задач ОПК носит локальный характер, Зависит от области значений^оценок характеристик конструкций и в той или другой степени присуща всем задачам ОПК.

Аксиоматический подход к решению МКЗ оптимизации подкреп­ ленных цилиндрических оболочек рассмотрен в 3.10.

В дайной главе на основе результатов, полученных в гл. 1, 2, рассматриваются основные постановки задач, относящихся к оп­ тимальному проектированию широкого класса подкрепленных много- и однослойных цилиндрических оболочек. Учитывается дис­ кретный характер подкрепляющего силового набора. Анализиру­ ются постановки задач, позволяющие отразить широту и многооб­ разие многокритериального подхода к исследованию проблем, возникающих при разработке моделей оптимального проектиро­ вания по нескольким критериям качества конструкций.

Так как оптимизационные задачи сводятся в основном к зада­ чам нелинейного программирования, принципиальные особенно­ сти которых подробно описаны □ 2.1, 2.4, то при их решении ис­ пользуются эффективные алгоритмы метода случайного поиска [17, 31]. Численные результаты получены в основном с помощью ЭВМ типа ЕС (язык программирования ФОРТРАН). Для всех типов рассмотренных оболочек детально анализируется поведение оп­ тимальных конструкций и даются конкретные рекомендации для их проектирования. Все обозначения приведены в гл. 1 и 2.

3.1. ОПТИМАЛЬНОЕ ПРОЕКТИРОВАНИЕ ОДНОСЛОЙНЫХ -ОБОЛОЧЕК ПРИ ОСЕВОМ СЖАТИИ

Рассмотрим шарнирно опертую по краям круговую цилиндриче­ скую оболочку радиуса г и длины L, подкрепленную продольными ребрами жесткости (стрингерами) и сжатую заданной осевой на­ грузкой М, которая осуществляется в виде равномерно распреде­ ленных по торцам оболочки сжимающих напряжений р (рис. 5). Критические напряжения определяются по линейной теории, в уп­ ругой области работы материала, характеристики которого (мо­ дуль упругости £, предел пропорциональности ff7 и коэффициент Пуассона/U ) известны. Стрингеры представляют собой тонкостен­ ные стержни открытого профиля (в виде прямоугольной полосы или равнобокого уголка), расположенные симметрично относи­ тельно срединной поверхности оболочки. Требуется найти такие значения толщины обшивки h, толщины ребра Ъ(высота ребра Н является зависимой величиной и определяется по норматив­

ным значениям

обеспечивающим местную устойчивость

н

стрингера), а также количество ребер к , чтобы при заданной на­ грузке оболочка имела минимальный объем и одновременно, вы­ полнялись условия прочности и устойчивости. Таким образом, нужно найти минимум функции:

l/=[2irrh +kF)Lt

(3.1)

где F - площадь поперечного сечения ребра, при следующих ус­ ловиях [2]:

(Z w rh + kF)G T ^ Ы;

(3.2)

2 7 r£ y ft2(1+

(3 . 3 )

2nrh

 

Условие (ЗЛ) - условие прочности, а условие (3.3) - ограничение по критической нагрузке потери устойчивости. На основании (1.46) и в предположении, что ребра работают только на изгиб в радиаль­ ной плоскости и кручение, а также что для расчета обшивки при­ менима техническая теория оболочек, после соответствующих преобразований для параметра у получаем такие соотношения:

1) общий случай деформации при потере устойчивости

1

 

* 6 7 у -------- ]

у

,

— « ( У ,2(1 v y

'

1*С(Р/27ггЛ)[Уз(1У)

2h‘rz

2rfh }

Гг

где

* / * ) ] ,

(ЬЛ)

/ г / , ;

2) первый частный случай деформации

 

 

 

 

 

 

 

V

 

-

 

I1(7ГНГЛ

 

 

 

 

 

 

 

(3-5)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

т2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A 'h

 

 

 

 

 

 

 

 

 

'З) второй частный случай деформации

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

А П У )

 

 

 

 

 

 

 

 

 

 

|2+— — ^ в а л

и

c * i ,

 

Ь = Vl2(l-/12J

 

JTrEh

 

 

 

 

 

(3.S)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

р

I

 

л С м ^ в й - У ) ,

 

 

 

/i2 ( i- /u * ;

*

-------- ------------ и*------------- J 1 е ш с > 1

 

I -

■ с г

 

тгг £ Н ь

 

 

 

 

 

ГАв

с -

 

 

 

 

г *

 

 

 

 

 

 

 

 

 

 

 

 

/T h " /« ( 1 .- //* )

 

 

 

 

 

 

 

 

Используя обозначения

 

 

 

 

 

 

 

 

 

 

« , = Л ;

 

 

г 8=& ;

* ,

= * ;

A , = t ir r L ;

 

 

 

 

,

 

 

 

 

 

 

 

27Г/*й_

 

 

GL

 

 

/12=ЦгЛ-1); И3=АА; V

- f i - 1-;

й2= Т (гЛ' 1);

о

1.

D

_

2 ic F ,

-n _

2AJ 1

,

_

A

_

-тч

_

^

 

N Л‘

 

 

 

 

 

 

гяг '

гя-г'

V

/3(1^

 

IH JrB j

 

 

 

 

я -

г Л2(Д-1)

2Л- 1

-

-

1

Д =

 

 

 

 

 

 

 

 

 

 

 

Ь

»

Л. Л ZZ

_______

3/.¥(1+|М )'

Be= i» 3

[

10

 

« р - н ) 1'

?" у5й У )

's=

 

 

24i!Dp’ ^ “Vfwyi'

 

 

AV

I )

 

 

2В°'

 

 

 

 

5ЛГ2 ’

 

 

 

 

г 1 Д !^ г О :.^ 1 .р = — • Ъ =2 < Д в =

 

ЗД,

 

" a

r d + f i ) 1

 

6W1 + U)

1

« 8 ’

л

т

1

«

 

V - l b - V , 2 D- D « b ^ E E i,

 

 

V

. 6^

1Л?

 

Z03-

-°18

 

fr(f+fl)

 

 

 

 

3

_

 

,

 

IPS '\/F

 

 

 

 

 

 

 

 

2>« " г * ( 1 + /|)

'

^1 2 (1 -jU*)

 

 

 

 

 

и подставляя их в (3.1) - (3.6), получаем две задачи (соответствен­ но для подкрепления уголком и полосой) частично целочисленно­ го НЛП (x s может принимать только целые значения).

Стрингеры из равнобоких уголков. Найти неотрицательные зна­

чениях,,:^ и £ Jfкоторые минимизируют функцию

 

У

-

А

М

и удовлетворяют ограничениям:

 

 

ei V V

l x j * 1'r

 

(38)

х л

Величина7 принимает наименьшее из следующих значений:

5> ? г = A

,

+ 1

/ <1+Д , * * * з I х 1). «сли

;

И

Д Д 1 + Г я 4 х * / £ г )

 

 

 

9 '

12 2 3 '

1 ' '

 

 

 

^

а .

т 4х.5

 

 

 

W

!

 

)•

е с л и с > 1 '

 

 

 

 

 

 

где с = у>х3 \ ^ .

Стрингеры из прямоугольной полосы. Найти неотрицательные значения x v x t и которые минимизируют функцию

У = А ^ * А $ х 1 х 3

0.10)

и удовлетворяют ограничениям

13.11)

(3.12)

где ч также принимается наименьшим из следующих значений:

1 + ®17г г

од

8)1?3 | л ,з<! + V f t !X0 /XIX1 ecjw <c> 1 .

Поставленные задачи нелинейного программирования (3.7) - (3.12) решались методом случайного поиска по методике, описан­ ной в 2.32.

В качестве примера здесь рассмотрим задачу отыскания опти­ мальных параметров цилиндрической стрингерной .оболочки при следующих данных: г = 0,5 м; L - 0,387 м; Е =. 6,87 • 101 Н/м ; IУт = 148 • 106 Н/ма; N= 25,4 • IQ4 Н; Л = 16,34. Стрингеры выпол­ нены: а) из разнобоких уголков; б) из прямоугольной полосы. Огра­ ничения (в метрах) на варьируемые параметры принимаются сле­ дующими: 0,001 0,005; 0 ,0 0 1 ^ 6 ^ 0,005; 4 ^ fc < 50. Учиты­ вается целочисленный характер переменной к (количество ребер).

Результаты решения приведены в табл. 2. Видно, что для оболо­ чек со стрингерами в виде уголков критические напряжения ми­ нимальны в общем случае деформации, а для оболочек со стрин­ герами в виде полос - в первом частном случае. Сравнение ре­ зультатов, полученных для оболочек, подкрепленных стрингерами из полосы,и стрингерами из уголка, позволяет сделать вывод, что для стрингер-полосы определяющими являются частные случаи деформации. При этом толщина стенки оболочки уменьшается, а количество стрингеров увеличивается (по сравнению со стрин­ гер-уголком). Объем оболочки с подкреплением в виде уголка зна­ чительно ниже объема оболочки, подкрепленной продольными ребрами в виде полосы.

В связи с этим значительный интерес представляет исследова­ ние поведения параметров оболочки при изменении осевой сжи­ мающей силы N. Так, при А/ = /^А /о (/Уо ~ 2 6 ,*М 0 ^) рассмотрим оболочки, нагруженные силами, являющимися безразмерными ве­ личинами М, = 0,5 - 4, для двух случаев: а) в качестве стрингера принимается уголок; б) в качестве стрингера принимается полоса.

В табл. 3 приведены значения оптимальных параметров оболо­ чек, объемы, параметры критических напряжений, значения кри­ тических и нормальных напряжений, полученных для указанных выше значений сжимающей силы. На основе этих данных можно сделать вывод, что с увеличением сжимающей силы: 1) намечает­ ся тенденция к увеличению толщины стенки оптимальной оболоч­ ки; 2) растет толщина стрингера; 3) количество ребер несколько уменьшается; 4) объем оболочки растет, но при небольших нагруз­ ках (NQ - 0,5 - 0,8) объем оболочки, подкрепленной стрингеромуголком, выше, чем при тех же нагрузках для оболочек, подкреп­ ленных стрингером-полосой. При увеличении нагрузки объемы оболочек для обоих случаев подкрепления выравниваются. Кроме

того, следует отметить, что определяющими при выборе оптималь ных параметров оболочки явились общий и первый частный слу. чаи деформации. На рис. 6 показана зависимость объема оболочки от значения сжимающей силы для подкрепления в виде уголка и полосы. Хорошо видно, что при небольших нагрузках оболочки, подкрепленные стрингерами в виде полосы, легче, чем оболочки, подкрепленные стрингерами в виде уголка. В дальнейшем с ро­ стом нагрузки это различие в весе исчезает. Рис. 7 демонстрирует изменение толщин стрингеров в виде полосы и уголка в зависи­ мости от нагрузки А/0,

 

 

 

 

 

 

 

 

 

 

Таблица2

 

 

 

Оптимальные

 

 

Параметры

 

 

 

 

 

параметры

п о

’,

критического

 

Тип

 

 

 

 

 

напряжения

 

 

 

 

 

 

СМ3

 

подкрепления

Л,

Лр,

К '

 

 

 

 

 

 

 

 

мм

 

МИ

ММ

 

 

 

h

и

 

 

 

 

 

 

 

 

 

 

 

 

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

0,28

1,26

34

0,961

 

4,99

3,81

3,86

Стрингер-уголок

0,35

1,24

26

0,917

 

2,35

2,36

2,52

 

 

 

 

 

 

 

 

 

 

Таблица 3

 

л,

V

К

К-ю 5.

 

 

 

 

б га5G

ГО"6

Стрингер

ИИ

ММ

шт

см3 '

1

*1

 

h

fy ,»

Я

/м?

*1

 

 

 

 

 

 

 

 

 

 

 

 

0,5

0,16

1

39

 

674,7

4,203

3,49

 

13,4

76,7

73,0

0.6

0,16

1

39

 

698,3

3,64

3,62

 

10,8

90,5

87,5

0.8

0,196

1,06

40

 

793,7

3.87

3,97

 

12,71

102,5

96

1.0

0,351

1,24

26

 

917,6

2,35

2,36

 

2,51

120,6

123

1,2

0,297

1,24

32

969

3,13

3,27

 

5,65

124

127

1,4

0,312

1,30

32

1043

3,70

3,64

 

8,27

157

137

2,0

0,311

1,55

34

1380

5,13

4,00

 

13,83

170

145

3,0

0,401

2,22

26

2071

8,33

3.14

 

15,48

173

147

4,0

0,482

2,83

22

2758

11,47

2,72

 

17,29

180

148

0,5

0,171

1,00

43

479

5,97

4,55

 

7,74

107

107,0

0,6

0,174

1,00

46

502

5,87

5,15

 

9,038

123

122,0

0.8

0,226

1.07

41

573

4,761

4,61

 

5,264

143

143,0

1,0

0,286

1,26

43

691

4,99

3,81

 

3,860

147

147,5

1,2

0,325

1,24

45

830

4,39

6,39

 

5,303

196

147,5

1,4

0,198

1,63

43

966

17,44

5,81

 

38,57

158

148,0

2,0

0,509

2,46

20

1382

10,93

2,46

 

3,55

172

147,5

3,0

0,359

3,28

24

2068

39,31

3,35

 

33,2

165

148,0

4,0

0,416

4,34

. 19

2763

60,81

2,59

 

37,32

143

148,0

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