книги / Оптимизация подкрепленных цилиндрических оболочек
..pdfИсследуем теперь влияние эксцентриситета расположения ребер на поведение оптимальных оболочек, используя методику работы [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~№ |
К |
|
|
||
г |
|
2± |
|
^11 |
** |
<1>л |
|||||||
3 |
|
1 |
|
|
|
|
'2 |
|
3 |
Р г |
|||
« |
|
- |
|
|
|
|
|
|
|
1 |
|
2 |
|
где Ф |
|
отношения высоты |
поперечного сечения профиля |
||||||||||
к его толщине для указанных выше элементов подкрепления. |
|||||||||||||
|
|
S |
f n |
|
|
_ |
s F |
|
|
* 5 |
|
|
|
|
г - |
|
|
|
1П |
|
|
|
|||||
|
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 * |
|
условия прочности |
|
|
|
|
|
|
|
||
<г |
-----------------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 =