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

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

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

11 Толщину экрана /ри с. 5 .3 / представим подинойом

*^Ot^Z+^itc^4c^S » /5 *0 7 /

коэффициенты которого определяется выражениями

 

 

 

 

 

 

 

 

Сдо® Wg.

 

 

 

 

 

 

 

 

 

 

 

 

С2о=(Н1 -НгУ<£н,

 

 

 

 

 

 

 

 

 

 

Саг~(Нз~Н2)1 сСгв1

г

 

 

 

 

 

Спк[1Н^И2У(НгНг) ^ ' ( И 3% )2

^ :]1и иd 2j \

 

(dL1ir<£ec)-

координаты

 

*on eifa

A i t

i

- 1 ,4 .

 

 

 

Дая вычисления векторов^чузстштальнооти

ij

наобхсдр^

i(o

определить_производные dJi/dZ

. Производные

от

толзин

h4

 

у . . . ,

hfi

оуперелементов евдана

по ущ ш влашы парамatw

рам

в га с л я в т по состоюениям

dtCpj , г

. dCgt

 

• гг

dhp

Г dCoo . dCia г 2 .

/2

J dHj

* T t i f * "

Щ ~ ***

йН+ « W

W

W

Щ

‘ \ о ,

j- & B .

 

 

 

 

 

 

 

 

где

 

c/ii4 , с / ? *

- коордлнатм центра

суперэхемент*

5 р

 

/5

23/.

 

 

 

 

 

 

 

 

 

 

т Матрица производныхdc/dH от коэффюрентов Cm(Coo,CtO'Cei> Cj2) полинома /5.87/ по уэловда втчэним толщин айрана

Н"(Н,. Ht.Hs.Hif"*” ввд

 

0

i

0

0

tic _

и *

- U u

0

0

 

 

 

йй

0

 

U U

о

 

 

f<4

4 l(d» d^ f -Л.

С. Толщина на борте /рис.

5 . 4 / представлена о тд ел ьн о

e

подобластях

Sif .

 

 

 

1 *3

 

A* /---- :------ A |

 

' U

r

 

A10 1

A 7

 

j o

g

 

As

AIff < 4

 

ftie. 5.4.

 

1. В подсылаете Q j

толщина ведается

полиномом

Л * CQO* C10cLf * CQ2 cLl

/5.88/

коэффициенты которого находят по значениям толщины Н5 . Н6

Н7 *HlQ

Соо~"ь,

Счо~ Ws H6)Jc(19 ,

Сог~1^т“Haifahi 2 С1г~[№ю'Нб)~(Нз~Н6)--££-(H7 -H6)J £ ■ ]l(cCno’d*iah

*- 1.52 -

*7

Производные от толщин

hp

супврэлеиентов

Sp, Sp £ Q-i*

борта по управляющим параметрам

определяют из выражений

 

Hi,

 

 

I

 

dH.

 

 

 

 

 

 

 

 

* Щ Л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

[ о

, \

jc {l,2 .3 ,i.b } .

*

 

 

 

 

 

'/

 

 

Sp

Шт-

где o t f *

. 0^2» -

координаты центра суперэлемента

 

тжца производных

d C .

от C-[c0o,Ci0,ln .CnYдо

Ня[Н5,Нд,Н7Щ

имеет

ввд

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

1

 

 

0

 

 

0

 

 

 

CL15

 

 

 

JL

 

0

 

 

0

 

dc

 

 

 

 

d < 5

 

1

 

 

 

 

 

0

 

 

-

 

'

 

 

 

0

 

dH

 

 

 

dm

d-27

 

d-if

 

 

 

 

 

1

t , L

dtfo)jj'j J*

\

1 -Л-

 

 

4

 

 

 

d<5

dir

 

 

 

ditod]

 

d^d2w l

 

£&g

 

 

 

 

 

 

 

if

 

В подобласти[

толщина задается полиномом

 

h-Coo+Ciod-i * cot(^£“d-zaf * Ciid^dg-cLgg)2,

 

/5 .6 9 /

коэффициенты которого вычисляют по знамениям толщин

W7 , / / j .

Нд , НЮ

 

 

 

 

 

 

 

 

 

 

 

 

 

Coo ~ Мд •

 

 

 

 

 

 

 

 

 

 

 

 

 

Сю~(Нд~Нд)/с(.м,

,

 

 

 

 

 

 

 

 

 

С’в 2 = ^ 7 * ' ^

) / с^г7

 

 

 

, *

 

 

 

 

 

С12ж[(Ню'Н8ННд-

 

 

~ ^ 7

~ft fa ir ]//Io £ f» e tf л?),

 

где

d l7 ~(cC2 7-dfp)

<^гющ(^-210^

2$)

 

 

Jm

 

Производные от hp

 

 

 

« о

 

 

 

суперэленонтов i p . 5 p ° c Q

 

по Н/

определяются соотношениями

 

 

 

 

 

 

 

dCooасоо Qcic^f-CIQ*

 

dc,ас9*//

~j

12. dCjg .

iy

/

■ it

 

dHj

d H j 1*

dHjP2* ^

9'

dH/ *< *™ ~ **1

dhp^J a n J

 

_

 

 

 

 

 

/ «

К

» .

 

 

Щ 1 0 ,

 

 

 

 

 

 

 

 

j = 1,6

 

 

 

 

 

 

 

 

 

 

- IS3 -

йС

ж

 

 

0

 

1

0

 

0

 

 

0

 

1

4

 

0

 

 

 

<к#

k ig

 

 

 

 

 

 

 

 

4 г

 

1

0

 

0

 

 

" /»

 

 

°-1Г

tktf

 

 

 

1

IT

4 + Ь ш А **\и

 

 

4

JC W

^

\

« м ? » '

л ,

110<к

L

 

 

 

>(С оо,С ю ,

Н =(Н 7, Н в. Н а ,Н ,о )Т.

 

Рис.5.5.

l.fo a p m яа коиуов, жл хе кик ■ на борте, представлена одаедано в водобшсмх Qi1/, Q f /рас, 5 .6 / е помощью пода-

нфмов, а ха л о тн ц х^б .в в /. /5 .8 9 /.

-

UpoEM W ie dkJ iH j,j- G 6 ,

, Р » ж ДОП,«ДО KMX

 

 

 

 

• деоцесое

опяшм-

 

Оболочка нахсщмтел под внешни давлением интенсивности

10"^ кг/шй н карактермвуется фиаичеокшк параметрами:

f

- 62.40 и * « г|

d -

0,3i

G'-EKZtM)) ; К1- 5/6

/Wtfepxai оболочка «отроги»*/,

Вибиралх оледупщие раамарн

баллона /мц/: Я*

649;

- П9б; Йгот 945;

24,7;

В0 =

172.5; H0

138,0 . В локальной « с т е л е

координат

 

«очки

В ,,

В2,Ан, В5 , Af g

имеют координате:

В / /- 1 5 ,

-5 0 /,

В2 / - 3 ,

-3 1 /,

Аде/О, О/,

В3 /- 1 0 ,

4 5 /,

А д е/-24,4,

8 0 /.

Точки

А5 , Аде, Aj,A16tAetA15 и

Аде ,

Л /Э , A f *

лежат

в соответствующих горияонтальник плоскостях.

 

 

 

Ьго. 5

.7 .

 

 

 

 

 

 

 

 

Значения толщин

Щ /мм/

в точках

Aj

.

/ ■ 1,16

в

на­

тальной конструкции выбирали равндаи: 9 ,2 4

; 8

,5 ; 9 ,7 6 ;

10,

27;

12,87; 9 ,5 ; 9 ,5 , 9 ,6 ;

12,81;

12,94; 9 ,5 ;

6

,1 7

;

6 ,1 7 ;

6 ,1 7 ;

 

9 ,5 ;

9 ,5 ;

Объем материала

/м асса экрана и борта/ до

оппаш аа-

цнн

J j -

З .Н Т О 5 мм5.

 

 

 

 

 

 

 

 

 

 

На рис. 5 .6 , 6 .7

кравши

I

показаны опоры напряжения

(5у/

на внешней поверхности оболочки в сечениях

AgA^AgA 3

/малое

сечение/ и

А2 А3 А5 А$

/большое

сечение/. Б малой сечении

баллона возняквот максимальные растягивавшие напряжения.

 

i e l , 4

 

В роди управляющих параметров выступали толщины

 

,

экрана, а также толщины

//£- , £ -

1.20 айрана и борта,

на

кото­

рые накладывались двухсторонние

ограничения: 5 ,0

«S

15 ,0 ;

i

ш1 ,3 ,

4, . . . » 10; 5

£

 

Ю. Прочностные ограничения

учитывали в

подобластях

<Qy

, Q 2 н выбирали в виде двух­

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

 

 

 

 

 

[ б

- ] < 6

ц < [ 6 J ,

[6 - ]< 6 (t*[ 6 t]

i=(,2.

/5 .9 0 /

 

Пределы изменения напряжений:

£ б - ]=

-3,0 кгЛп-1 '

 

 

[ б у » 0 .9' кг/мм ,

 

 

 

 

 

 

 

 

 

 

Задачи оптимизации массы решали методом внешних штрафных

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

Т а б л и ц а 5. 1

п/п

« <

Н± -

' 6,02

IU 1 Jo/lO1-.

К

 

 

в|09

6 ,6 5

2.63

 

16

I

~ 5 ,4 4

132

2

7,92

5,04

7

55

5,12

2,62

254

16

3

8,22

5 36

6,37

5,35

2,62

366

Гб

П р и м е ч а н и е : 1 - результаты, полученные методом Ялетчера-Пауэлла * применением <ю1фяжешшх уравнений; 2 - е при­ менением кокесщо-разностных аппроксимаций; о - методом помока По деформированному многограннику.

В первом случав управление осуществлялось толщинами

,

L « 1 ,4

экрана цри нажженных /равных начальным/ толщинах бор­

та. В табл. 5 .1

приведены полученные в результате

опгомигации

 

толщины вкрзна

Н* , I ■ 1 ,4 , значение оптимизируемого функцио­

нала do

/массы бадлона/,

количество вычислений

К целевой

 

.функпия,

ц>оцонтннЙ вьогрш

d *

 

 

Пожученные различили методами оптимальные проекты равно-

 

Дэнчы, хстг значения опташльных управляющих параметров в них

 

иотчэжмг. отличаются. Однако с точки зрения вычислительных

 

затрат, мотод Флетчера-Пауэлла с применением сопряженных ураз?-

 

ц«НкЯ для

гррзделвния «феноводнях намного эффективнее.

 

На рис. 5 .6 -5 .7 кривые 2 отракмот эщры напряжения 5# соответственно в мадам и больше» сечениях оптимальной конструк­

ции

В точке

Ду

напряжения в оптимальных 1фоектах достигают

своего додуста! ого

значения.

 

 

 

 

 

Чмеленндаи экспериментами установлено, что уменьшение тол­

щин

Hg I

HJQ

борта приводит к значительному уменьше­

нию напряжений в

экране

баллона. Кривые 3 /рио.

5 .6 , 5 .7 /

иллюстрируют егшры напряжения

бу

^соответствую щих сечениях

.для

вначений

толщин

- 5 ,0 ;

t =

5,10 /толщины экрана рав­

ны начальниц/. При этом напряжения на

борте находятся в допус­

тимых пределах.

Объем данного

баллона

J0= 2,24*10э мм3 , т .е ,

на 28 %легче исходного.

 

 

 

, • . . > Н10

 

Поэтому в случае управления толщинами

за начальное приближение выбирали конструкция со следующим рас­

пределением

толщин:

Н<* 9,24; //2 »

8 ,5 ; Н$» 9 ,7 6 ; Н±« J0 ,2 7 ;

И5 « . . . =

5 ,0 .

Для достижения

оптимального проекта мето­

дом сопряженных градиентов с применением сопряженных уравнений потребовалось 160 перерасчетов баллона. При этом 14 раз вычисля­ ли градиент оптимизируемой функции. В результате получали сле­

дующее распределение толщин:

Wy «= 6 ,3 9 ;

5, 0;

Нэ -

5,7 2 ;

^ 4

- 6 ,0 ;

т . . . в Н10= 5 ,0 , Толщины борта баллона

в

про­

цессе

оптимизации не изменялись; толщины экрана сильно умень­

шались / в центре

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

рсое?

го нижнего предела/. Объел полученной конструкция

1 ,6 4

*1

0 ° мм ,

т .е . на 47 %легче исходного

баллона. Следовательно,

подбирая

одновременно толщину экрана и борта баллона, можно значительно

понизить его вес.

+

Кривые 4 /рис. 5 .6 ,

5 .7 / иллюстрируют этюры напряжения 6Г,и

в малом и большом сечении баллона. В указанных сечениях нормальные напряжения находятся внутри допустимой области /5 .5 5 /, Одна­

ко вблизи диагонального сечения вдоль линии

А+, А5напря -

женин достигают своего верхнего предела.

 

Пример 5 .2 .

Рассмотрим

баллон кикескспа

/рис. 5 .2 /,

находящийся

под

воздействием

внешнего давления. Исходные данные

конструкции

такие же,

как

в предыдущем примере. Задача оптими­

зации

состоит в

том,

чтобы подобрать толщину экрана бадлонр,

объем

которого не превигоал бы заданного значения. о це;ак> .<инж-

миэации возникающих в нем напряжений

 

 

 

 

Jc^min

max о{в),

 

 

 

 

2 е в ,

Q.QZi‘

 

 

J,№1* J ° .

nt -{2:

H , Hj, H4] T, Hi* Hi * Hi*, £ - / 7 },

 

t(S)-max {<?« ,G*t,6 u ,6« },

Hi-5,0,

i - i j ; Ht-19,0, 1‘ Ш ;Н *г -Ю,1?*3,НЮ*

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

ностных аппроксимаций.

 

 

 

Принимаем з а

начальное приближено исходную конструкцию

2 ^

■ /9 .2 4 ;

8 ,5 ; 9,7 6 ; 1 0 ,2 7 /. Однако методом Флетчера-Пауэл­

ла не удается

уменьшить напряжения в баллоне. Поэтому

выби­

раем конструкций оо следующим распределением толщин:

 

2 ^ -

/1 2 ,0 ; 9 ,9 ;

12 ,0 ; 1 2 ,0 / /недопустимая по объему конструк­

ция/. В результате

оптимизация в этом

случае максимальные .растя­

гивающие напряжения снижаются с 0,74

кг/мм^ до 0,60 кг/мы _на

19 %. При атом

подучаем следующее распределение толщины;

2 * -

■ 10,52; 9,7 0 ;

8 ,2 6 ; 9,96 . По сравнению с исходной конструк­

цией,

толщина в малом сечении и в центре баллона увеличилась,

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

Отметим, что как с применением сопряженных уравнений,

так

ж с помощью хонечно-разноотных аппроксимаций производных,

зау­

чены. одинаковые оптимальные проекты. В первом случае для дости­ жения оптимального проекта провели 37 перерасчетов конструкции 3 вычисления градиента/; во втором - 51 расчет конструкции.

Пример 5 ,3 , С целью рационального распределения материа­ ла и повшения несущей способности тройникового соединения /пример 3 .9 / рассмотрим задачу построения дискретно-равнопроч­ ной оболочечной конструкции. В роли 6 нщ? выбирали максималь­ ное по модулю значение напряжений на суперэлементе. Управление осуществляли толщиной трубы и патрубка. Для решения этой задачи

потребовалось

7 перерасчетов тротилового соединения. Напряжен­

ное

состояние

в полученной конструкции изображено на рис. 3 .25 -

3 .2 7

штриховой линией.

Заметай, что повторные перерасчеты конструкции проводили по сокращенной схеме, описанной в параграфе 5 .2 . Вследствие в того затраты машинного времени на построение дискретно-равно­ прочного тройникового соединения составили половину еатрат на один первый прямой расчет.

В полученной дискрето-равнопрочной конструкции максималь­ ные напряжения снизились на 60 %. При этом толщина патрубка увеличилась на 94 %, а толщина трубы уменьшилась на 8 %. В результате объем конструкции получился на 4 ,4 %меньше, чем исходной.

 

 

 

 

СПИСОК ЛИТЕРАТУРЫ

 

 

 

 

 

 

 

1. А б о в с к и й

Н.П., А н д р е е в

 

Н.П., Д о р у -

г а

 

А.П. Вариационные принципы теории упругости и теории обо­

лочек. М.,

Т97В. - 287 с.

 

 

 

 

 

 

 

 

 

2. А м б а р ц у м я н

С» А. Общая теория

анизотропный

ободочек» М., 1974. - 44В с .

 

 

 

 

 

 

 

 

3.

Б а з а р а

М., Ше т т и

К. Нелинейное программиро­

вание. Теория и алгоритмы. М., 1982.

-

534 с.

 

 

 

 

4. Б а л а н д и н

В. А. . Ш а п о ш н и к о в

Н.Н. Расчет

осесимметричных круглых пластин и оболочек по методу конечных

элементов с учетом поперечного сдвига / / Сб. научн.

тр. ВНИИХМ.

1978.

№ 79.

С .18-24.

 

 

 

 

 

 

 

 

 

 

5.

Б а

н и ч у к

Н.В. Оптимизация форм упругих

тел. М.,

I960.

-

255

с.

 

 

 

 

 

 

 

 

 

 

6.

Б а

н и ч у к

Н.В. Современные проблемы оптимизации

конструкций

/ / Иэв. АН СССР. Механика

твердого

тела.

1982. Ж 2.

С. 110-124.

 

 

 

 

 

 

 

 

 

 

 

7 .

Б а н и ч у к

Н .В .,

К о б е л е в

З.В. Некоторые

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

ободочек

вращения / / Изв.

АЛ

АрмССР.

Механика.

1983. Т. 86. № 2.

С .10-17.

 

 

 

 

Б а

н и ч у к

Н.В.

Введение

в

оптимизацию конструкций.

М.,

1985. -

302 с»

 

 

 

 

 

 

 

 

 

 

9 .

 

Б а т е

К .,

В и л с с

н

Е. Численные методы анализа

иметод конечных элементов. М., 1982. - 447 с.

10.Б е р к о в Н.А. Исследование напряженного состояния соединений цилиндрических оболочек: Автореф. дис. . . . канд. техн. наук. М., 1982. - 24 с.

11.

Б о г о м о л о в

С. И. . Ж у р а в л е в а

А.М. Коле­

бания сложных механических систем. Харьков, 1978.

- 136 с .

12.

Б о г о м о л о в

С.И .,

Л у ц е н к о

С. С. . Н а з а ­

р е ч к о

С.А.

0

применений суперпараметрического

оболочечно­

го конечного элемента :t расчету колебаний лопаток

турбгмащии / /

Пробл.

прочности,

£982.

№ 6.

С.71-74.

 

 

13.

Б о

г о

с л о

б

с

к а

я

Н.М., 3 у е в

Б .К ., К а п у с ­

т и н

С, Л.

Применение простейших сдвиговых моделей ККЗ для

расчета

осесимметричных

 

оболочек / / Пршсл. пробл. прочгюсти

и пластичности.

1930.

Зып.

£5. С.95-99.

 

 

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