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

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

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

Исследуем теперь влияние эксцентриситета расположения ребер на поведение оптимальных оболочек, используя методику работы [15], основанную на использовании теории конструктивно-, ортотропных оболочек.

Задача НЛП (с учетом эксцентричности расположения ребер)

будет иметь вид

 

 

 

 

m in G (X) = х . *

,

(3.13)

 

 

1

x h

 

при условиях:: * 2 * ° :

<$ С П г о ,

 

Здесь

=

 

xt =b . i s =o,QShc; i 4 = o,osJ;

 

 

 

 

(3.14)

♦ . и - М

. - г Ц ^ * ? ) - ;

(3.15)

’ ■ ‘ " ' " • • • ' e

y

 

(316)

* > W ■ В ? ' ’ l { [ « . s , (Ь С , г , ! 1 . V . 4 ) f

- W 4

 

 

 

 

а »

% ^ ^ Х3 '0 ,0 В } (P? i X ) ^ x 1 O.U5,

11 is»;

л хг

£

г* £ .

A

___ u.

<00

1 -

В

ж к ы ‘>

 

 

N

т

» n 2 - 12(1-/U2j

N '

2

 

£ ]Г

 

£ \

д -

* * / *

_ £ .

4

s O j r R

- 0 .= — 1

______ £ .

 

J = v W > " ' 5' (' - ^ N' 4

 

■ * W ) « '

«

л .4, j

_

 

n%_________ £ •

£ - { £ E . f .

 

4 = ? $ - . ^ ' с = е т (*« * « * . ) ;

аМ)-/и*)***3

где Zc > 0 при наружном расположении стрингеров и z c < 0 - при внутреннем их расположении.

Выражение (3.14) - ограничение по критической силе потери устойчивости отдельного стрингера, ограничение (3.15) исключает потерю устойчивости панели, заключенной между стрингерами и шарнирно опертой по четырем сторонам. Выражения (3.16) и (3.17) - ограничения (с точностью до множителя У ц ) по критической на­ грузке потери устойчивости оболочки по осесимметричной и не­ осесимметричной формам соответственно. Выражение (3.16) - условие прочности и, наконец, (3,19) - геометрические ограниче­ ния на размеры.

Для численного анализа рассмотрим оптимальное проектиро­ вание шарнирно опертой оболочки при следующих исходных дан­ ных: R = 1,5 м = 8 *101 Н/м2; в0 - 157 Н/м2; р = 0,3. Расчет проводился для различных значений осевой сжимающей нагрузки А/ и отношений L / R , На варьируемые параметры накладывались о г р а н и ч е н и я : 3- 10 э м ;0 £ Лс^ 1,0•10~2 м; 0 « I $ as 2л^м .Р езупьтаты приведены в табл. 4; их анализ показывает, что при осевом сжатии таких оболочек расположение подкрепле­ ния практически не влияет на параметры оптимального проекта (через к * и к в табл. 4 обозначено соответственно количество стрингеров при наружном и внутреннем их расположении).

 

 

 

г с => 0

 

 

<

 

z z < О

и

L,

 

 

 

 

 

 

 

 

М

8 -io_f

Чг-Ю’ *

к с ю г,

е •10~2, h+

 

£ -ю "Л

 

 

G

<7

м

 

м

м

' м

м

t i

1

1,5

0,066

0,05

0,438*

3,424

94,2

10

0,066

0,05

0,295

3

1,5

0,198

0,169

0,495

3,985

67,3

14

0,198

0,182

0,520

5

1,5

0,330

0,328

0,494

3,97

94,2

10

0,331

0,321

0,550

7

1,5

0,458

0,438

0,440

4,295

94,2

10

0,465

0,380>

0,664

9

1,5

0,606

0,515

0,738

5,263

42,8

22

0,501

0,508

0,725

1

3,0

0,066

0,05

0,252

3,429

55,4

17

0,066

0,05

0,382

3

3,0

0,204

0,190

0,226

3,012

49,6

19

0,197

0,186

0,302

5

3,0

0,331

0,320

0,551

4,29

235,6

4

0,331

0,321

0,550

7

3,0

0,461

0,439

0,440

4,444

94,2

10

0,459

0,449

0,352

9

3,0

0,590

0,505

0,726

5,241

44,8

21

.0,592

0,563

0,430

1

4,5

0,066

0,05

0,339

2,753

58,9

16

0,069

0,054

0,380

3

4,5

0,204

0,192

0,294

3,221

85,6

11

0,197

0,182

0,361

5

4,5

0,331

0,320

0,551

4,312

235,6

4

0,331

0,321

0,550

7

4,5

0,459

0,452

0,315

4,082

188,5

5

0,462

0,441

0,440

 

______________I'

0,591

0,506

0,725

5,247

44,8

2*1

0,593

0,508

0,725

 

 

 

 

 

 

 

 

 

 

32. ОПТИМИЗАЦИЯ ПАРАМЕТРОВ ОБОЛОЧЕК С БОЛЬШИМИ ПРЯМОУГОЛЬНЫМИ ВЫРЕЗАМИ ПРИ ОСЕВОМ СЖАТИИ

Ребристые цилиндрические оболочки, являющиеся во многих случаях основными несущими элементами тонкостенных конст­ рукций, часто имеют различные вырезы. Это определяет актуаль­ ность их оптимального проектирования по весовому критерию при удовлетворении конструктивным и технологическим требованиям. В данном параграфе рассмотрены продольно сжатые ребристые цилиндрические оболочки с большими прямоугольными вырезами, подкрепленными: а) замкнутыми рамами; б) окаймляющими стрин­ герами и шпангоутами повышенной жесткости.

Ниже приводится методика расчета на устойчивость [6; 9], не­ обходимая для последующего синтеза соответствующих оптими­ зационных моделей.

3 .2 .1 . Расчетные формулы для определения критических напряжений

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

Используются следующие условные обозначения: Л,Г/2 - тол­

щина, радиус, и длина оболочки;

 

- ширина и длина отвер­

стия; G1

- шаг расстановки стрингеров и шпангоутов; f t 7,7кр -

площадь,

момент

инерции

при

изгибе

и кручении

стрингеров;

-

момент инерции

при

изгибе

и кручении

шпангоутов;

то же для вертикальных элементов рам, подкрепляющих

отверстия; 7|П,2|*Р

- по же для кольцевых рам, подкрепляющих

отверстия; Et f i -

модуль упругости материала и

коэффициент

Пуассона; т , п - число полуволн по образующей и число полных волн вдоль окружности; к, кь Нг - числа стрингеров, шпангоутов, отверстий; J2 - моменты инерции при изгибе и кручении дополнительных шпангоутов по кольцевым краям отверстий;

F

ез

вз1*

£3.

• Г = ---------

; с< = -------- ;

( 3 ----------

оГ ~ ------

2TTrh

27ГГд

i j r r d

* I D

 

КР

 

r j. J

 

, ш В L • Д

- i l l ----- .

 

Г, in

'

т

-

ц г)

 

m frr

t

-

b -

 

m =

 

L

 

 

Г

 

Потеря устойчивости рас­

 

сматривается

в

линейной

 

постановке,

дополнительный

 

прогиб задается выражением

 

 

Г

cos

„ ■Н •

 

 

 

 

(3.20)

 

Параметр

позволяет пере­

 

носить

начало

координат

 

в окружном

 

направлении.

 

Принимается, что

докритиче-

 

ское напряжение безмомент-

 

но, а напряжения в обшивке

 

ах = - р

 

= г = о, при

 

этом в стрингерах действуют

 

усилия

/Vc = р £ , а в стойках

 

подкрепляющей

отверстие

 

рамы усилия

 

 

 

 

Np=-o.sN c l^ -- n )- o ,5 p b h .

Рис. 8

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

метричная формы выпучивания; число волн п связано с числом

£

отверстий соотношением п * л 1

где гц = 1 ,2 , 3 ,... при чет­

ном и п 2,4 ,6,... при нечетном

 

 

Возможны следующие случаи деформации [6]:

 

1) стойки и ригели ра^ы изгибаются и закручиваются

 

= — г - ; ;

т

=г,ц, 6, ...;

(ЗЛ1)

2) стойки рамы только закручиваются, ригели только изгибаются

и

s

 

(3-22)

m

 

3) стойки рамы только изгибаются, ригели только закручива­ ются

d, * ~JT (« г 1), л? = 1,3,5,

(3.23)

4) стойки и ригели рамы только изгибаются.

d

1/

; т » 1,3,5... ;

(3.24)

*

 

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

Полная энергия системы;

Э - U

о

+ U

+ U

+ U

р

+ Тп+ Т + Т ,

 

с

ш

 

о

1

р’

где UQl Uc, иш, Up

- потенциальная энергия деформации соот­

ветственно обшивки, стрингеров, шпангоутов, рамы; Г0/Т1, Тр - работы внешних сил, приложенных к обшивке, стрингерам и стой­ кам.

Условие Л3 = 0 приводит к следующему выражению для па­ раметра критических напряжений:

я А0+ А с + А ш+ Af ,

Выражения для А0 и BQ имеют вид

4 т г п *

 

 

/72

tjmhn2)2

( т К п 1)1* М

t(m\nz)

_

-чЛ171

 

 

m L i

( - я

 

sin

 

irm

 

/•

»

 

 

 

_

М ) Я<

 

 

П В

К sin

—j r ~ cos 2 n d ^ ,

£ =

 

 

if r n

а формулы для A p,5p зависят от слрая деформации (3,21)—(3^4). Ниже приведены выражения Ар и ир для первого случал дефор­ мации (3.21):

г

27ГГ

.♦ V e

, I

t + ty 3

А Р~ ‘*С1 к г [

^

J*n—

 

 

 

 

 

 

 

, 8

(1

 

 

Va)

 

 

+ l j + ^ ] - £ l

 

;

 

 

ep=0,5C 3 [ | ,* , ( ^ -

 

 

 

2 Г

mZ

nB

г г

m l.

п д

 

 

ЭSHI

Stn‘

5 т

sin

;

^2 л

 

 

г г

 

2 г

1 Г

'

пб

 

 

i r

С = t r

Li 6 3 m

Kp

cr

 

f ns

Зтг(1 -/Цг ) 1

дг

 

 

BEJ,П

 

£2ia

 

<l>2

- углы

перекоса соответственно

вертикального и гори-,

эонтального элементов рамы.

Выражения для остальных параметров, входящих в (3.25), а так­ же для А р и Вр в остальных случаях деформации ввиду их гро­ моздкости не приводятся.

Относительные жесткости вертикальных, горизонтальных эле­ ментов рамы и дополнительных шпангоутов определяются по фор­ мулам

1

j 4 U В Л/

Л/

Л

f ■*

У л. лI _д_

п

2 к г

 

£-2

 

 

а ;

Р

„ = - £ 7

- --------- *

7

 

 

 

 

 

 

 

 

Г *

 

* „ ь е

*

! . ? * $ $

A

,

et' - ш ) '£ .

p

J L (a te )

1/1

1

 

^ «< 1

 

V

1'

"in

3

F>

 

1

 

 

Bs

p in t

п

А

-

я(1~№

К

 

 

г

 

 

^11

**

<1>л

3

 

1

 

 

 

 

'2

 

3

Р г

«

 

-

 

 

 

 

 

 

 

1

 

2

где Ф

 

отношения высоты

поперечного сечения профиля

к его толщине для указанных выше элементов подкрепления.

 

 

S

f n

 

 

_

s F

 

 

* 5

 

 

 

г -

 

 

 

 

 

 

 

1r r S

9

л -

iS

 

r, = L 8

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

jrr

 

 

 

 

(3.27)

 

 

rf

 

4

 

 

 

 

 

 

 

 

 

 

 

CL1

* 2ff

 

 

 

 

где t n ,f:ir t,

a $ ;

 

 

 

 

 

площади поперечного сечения вертикальных го­

ризонтальных элементов м|»до. ajajoxe дополнительных шпангоу­ тов; коэффициенты В, К, г , - относительные парамет­ ры изгибной жесткости, жесткости при кручении и площади попе­ речного сечения для рассматриваемого профиля подкрепления[6].

П о д к р е п л е н и е о т в е р с т и я с п л о ш н ы м и с т р и н ­ г е р а м и и ш п а н г о у т а м и . Часто подкрепления отверстий осуществляются вертикальными окаймляющими ребрами, длина которых равна длине оболочки; поперечное сечение ребер посто­ янно по длине. В сочетании с сильными шпангоутами, установлен­ ными на уровне края отверстий, эти стрингеры образуют систему подкрепления, выполняющую функции описанной выше рамы. Ме­ тодика расчета при таком варианте подкрепления остается в ос­

новном той же. Отличаются только выражения для Ир и 0р

кото­

рые приведены ниже

 

/ L - в о < Л

т Ч

(Ш)

п " г

12(1-/иу ;(1+5»)+2W

* - [ #

+1 Ы

+ Ч ( % ^ < ) ( U S , ) ,

(3.29)

,

я, n В

cosin d t

 

-ffle ^ = ('1 )

cos-y

 

Указанная система подкрепления эффективнее по сравнению с подкреплением рамами, так как мощные окаймляющие стринге­ ры связаны не только с обшивкой, но и с граничными диафрагма­ ми. Формулы (3.28), (3.29) учитывают совместность деформаций обшивки и ребер не только по прогибу, но и по углам поворота и закручивания.

3 .2 .2 . Определение оптимальных параметров

Рассмотрим шарнирно опертую по краям круговую ребристую цилиндрическую оболочку с заданными геометрическими разме­

рами Lf Q.VQ/ (-2 (рис. 8). Осевая нагрузка р -— равномерно рас­

пределена по торцам оболочки. Обшивка ослаблена большими прямоугольными вырезами, которые расположены регулярно в ок­ ружном направлении и симметрично относительно краев. Вырезы подкреплены по контуру рамами постоянной жесткости; в уровне криволинейных краев установлены окаймляющие шпангоуты. Обо­ лочка выполнена из материала с заданными физическими харак­ теристиками Е, (*•

Для анализа влияния отдельных факторов на оптимальное ре­ шение удобнее пользоваться такой системой параметров, в кото­ рой каждый параметр является независимой величиной.- Введем следующие безразмерные параметры, характеризующие вес, раз­ меры и число подкрепляющих оболочку и отверстия ребер: й} -

ЛК- отношение объема ребер к объему оболочки; о71 =

= % /Уш ^отношение объема стрингеров к объему шпангоутов;

t - n f f - относительная толщина обшивки; 0 ,-d jh - относительное расстояние между стрингерами; 0.^=а * /к - относительное рас­

стояние между шпангоутами; f t fa, отношение площади по­ перечного сечения вертикальных и горизонтальных элементов рамы, а также дополнительных шпангоутов к соответствующей площади сечения обшивки. Принято, что элементы подкрепления выполнены из тонкостенных профилей типа уголок, двутавр, тавр, швеллер или Z . Будем считать, что потеря устойчивости оболочки с отверстиями при недеформируемом контуре дополнительных подкрепляющих элементов происходит так же, как в ребристой оболочке без отверстий, усиленной рядовыми ребрами жесткости. Поэтому для описания выпучивания, в котором участвуют только

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

Методика, описанная в 1.4, позволяет рассмотреть такую де­ формацию, при которой теряет устойчивость лишь панель между отверстиями. В этом случае полагаем в выражении (3.20), что дли­

на выпучивающегося участка равна L * и, следовательно, м -

иг

Выражения для А0, BQ, нетрудно получить, приняв L-L

исоответственно изменив характеристики жесткости шпангоутов

икольцевых элементов рам, т. е. заменив L 2 на L в выражениях

для d 1t h -Таким образом, учитываются все характерные случаи потери устойчивости ребристой цилиндриче­ ской оболочки с отверстиями.

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

а) для оболочки с отверстиями

 

/Укр = у E t f i + f K ) 2 ? rr/i;

(З.Э0)

'б) для оболочки без отверстий, подкрепленной сеткой ребер

Nf,p = , t E f ( 1 ' * f * ) 2 r r W l.

(3.31)

Для расчета на устойчивость требуется при определении парамет­ ра критических напряжений fy.) в каждом из возможных слу­ чаев деформации при потере устойчивости выполнять минимиза­ цию по параметрам волнообразования Ш и п . Минимизация осу­ ществляется путем перебора т и п по алгоритму, описанному в 2.3.2., с вычислением на каждом г-м шаге JL . (f y j) для соответст­ вующего случая деформации и последующей выборкой из всех

VSj(V ij) минимального $£?иг(р™1п),

которое

принимается. в

расчет.

_ _

^

^

Обозначим независимые параметры СО, 0 )

^ ,

0 . соот­

ветственно через x it хг, т 3,

Тогда задача определе­

ния оптимального подкрепления формулируется как задача НЛП. В качестве функции цели и ограничений; зависящих от вектора

X

з д выбраны:

 

 

объем конструкции

 

*0 = I r { 2 J r r if 1+ и х 1х г ^ ~ ( х ^ х 7 * 1 i l l

и х

 

Кt, FпI 2

 

 

я т т L + x « 1 _ i B k 2 Г t *

 

Xjm in ^ Xj * Xjm a x>

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

 

 

 

 

 

 

 

-----------------2яггд?г ( - ---------------------<

-

 

«

&r i

 

 

0

JC£ ))

 

 

 

 

 

e0lFcQ+a/a)-dxr )

 

 

 

 

 

iF n

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

Н”

= 2Ю 1 ™ Е хг 1и

'

* ) * / / ,

 

3=

1, 2, , 7

(5.32)

r tKP = 2 jr fm in f l i

X

X

) * N

 

t ' -

1,2, - • • 1

(555)

1

*

>

 

 

 

 

i + x t

J

 

 

 

где 6"0 - напряжения в обшивке; (Гр - напряжения в стойках рамы при действии усилия Л/р - N£ (В /а + 1 )-Q ,6 fi0 B h .

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

£ = / ( X J - ► m i n ,

X e F ^

(аз4)

при ограничениях по прочности:

 

д о

?(12) = 0-о-ат *о, дг [Х) = Вр-6т^а,

по устойчивости:

 

 

J 4 ( i )

= N - N j p « Q

(15 6 )

и конструктивно-дахнологических ограничениях:

 

7 = 1 , 2 , . . . , 7 .

(3.37)

Особенностью рассматриваемого класса задач частично це­ лочисленного НЛП является то, что оптимизируемые параметры одновременно являются функциями • целочисленных параметров ролнообразования в окружном и осевом направлениях. Следова­ тельно, функции, стоящие в ограничениях, в области изменения оптимизируемых параметров терпят разрывы и применение поис­ ковых алгоритмов, предназначенных для минимизации непрерыв­ ных функций, "эляетсл невозможным. При изучении особенностей задачи (3.34) - (3.37) был обнаружен ряд локальных минимумов на границе области поиска.

Для численной иллюстрации рассмотрим оптимальное проекти­ рование ребристой цилиндрической оболочки (рис. 8) при следу­ ющих значениях исходных данных: L = 40 см; Г = 20. см; к= 4; Lg =

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