книги / Расчет и оптимизация оболочек с резными срединными поверхностями
..pdf11 Толщину экрана /ри с. 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. |
|
|
||||||
|
|
3» |
Б а |
н и ч у к |
Н.В. |
Введение |
в |
оптимизацию конструкций. |
||||
М., |
1985. - |
302 с» |
|
|
|
|
|
|
|
|
||
|
|
9 . |
|
Б а т е |
К ., |
В и л с с |
н |
Е. Численные методы анализа |
иметод конечных элементов. М., 1982. - 447 с.
10.Б е р к о в Н.А. Исследование напряженного состояния соединений цилиндрических оболочек: Автореф. дис. . . . канд. техн. наук. М., 1982. - 24 с.
11. |
Б о г о м о л о в |
С. И. . Ж у р а в л е в а |
А.М. Коле |
||||||||
бания сложных механических систем. Харьков, 1978. |
- 136 с . |
||||||||||
12. |
Б о г о м о л о в |
С.И ., |
Л у ц е н к о |
С. С. . Н а з а |
|||||||
р е ч к о |
С.А. |
0 |
применений суперпараметрического |
оболочечно |
|||||||
го конечного элемента :t расчету колебаний лопаток |
турбгмащии / / |
||||||||||
Пробл. |
прочности, |
£982. |
№ 6. |
С.71-74. |
|
|
|||||
13. |
Б о |
г о |
с л о |
б |
с |
к а |
я |
Н.М., 3 у е в |
Б .К ., К а п у с |
||
т и н |
С, Л. |
Применение простейших сдвиговых моделей ККЗ для |
|||||||||
расчета |
осесимметричных |
|
оболочек / / Пршсл. пробл. прочгюсти |
||||||||
и пластичности. |
1930. |
Зып. |
£5. С.95-99. |
|
|