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

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

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

ББК 22.261 335

УД1С625Г04

Оптимизация подкрепленных цилиндрических оболочек / В, А. Заруцкий, Ю. М. Почтман, В. В. Скалозуб. - К.: Выща шк., 1990.-138 с. - ISBN 5-11-001853-7

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

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

Табл. 20. Ил. 21. Библиогр.: 35 наэв.

Р ецензенты : д-р техн, наук, проф. В .И,Гуляее(Киевский инженерно-строительный институт), д-р техн. наук, проф. В. С. Гуд- 'рамович (Институт технической механики АН УССР)

Редакция литературы по математике и физике Редактор Л И. Попеначвнко

з1603040000—107 61-90 - М 211(04)-90

ISBN 5-11-001853-7

© В. А. Заруцкий, Ю. М. Почтман,

 

В. В. Скалозуб, 1990

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

Если методы прямого расчета подкрепленных цилиндрических оболочек при различных воздействиях достаточно подробно отра­ жены в литературе (например, в монографиях [2 - 4]), то среди из­ вестных нам публикаций практически отсутствуют обобщающие работы, посвященные оптимальному проектированию этих слож­ ных конструкций.

Цель настоящей монографии - по возможности восполнить этот пробел.

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

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

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

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

лась следующим

образом: В. А. Заруцким написана гл. 1;

Ю. М. Почтманом -

гл. 3; В. В. Скалозубом - гл. 2, параграф 3.10

гл. 3 и приложение.

 

В каждой главе используется своя независимая система обо­ значений, так как авторы глав 1 и 3, в которых излагаются методы и результаты расчетов, использовали обозначения, принятые в работах по механике подкрепленных оболочек [2 - 4], а автор гл, 2 и приложения - обозначения, принятые в литературе по оп­ тимальному проектированию [14; 31; 35].

Глава 1. ЭЛЕМЕНТЫ ТЕОРИИ И НЕКОТОРЫЕ РАСЧЕТНЫЕ ФОРМУЛЫ

В этой главе излагаются исходные соотношения теории, а так­ же методика и отдельные результаты решения задач о собствен­ ных колебаниях, а также устойчивости многослойных ортотропных ребристых круговых замкнутых цилиндрических оболочек. Использованные ниже соотношения получены путем обобщения на многослойные ортотропные оболочки соотношений, приведен­ ных ранее в [2; 3; 4], где рассмотрены однослойные оболочки, об­ шивка которых изготовлена из изотропных материалов. Сформу­ лированы также некоторые подходы к решению задач оптимиза­ ции ребристых оболочек.

1.1. ИСХОДНЫЕ ПРЕДПОЛОЖЕНИЯ И ЗАВИСИМОСТИ

Ребристая оболочка рассматривается как система, состоящая из собственно оболочки (обшивки) и соединенных с ней по линиям контакта ребер. Принимается, что напряженно-деформированное состояние многослойной ортотропной обшивки может быть опре­ делено в рамках теории упругих тонких оболочек, основанной на гипотезах Кирхгофа - Лява, а для расчета ребер может быть ис­ пользована теория криволинейных стержней Кирхгофа - Клебша.

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

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

Рио. 1

Связь перемещений точек осей j -го продольного ребра и j i -го кольцевого ребра с перемещениями срединной поверхности об­ шивки, а также углов закручивания поперечных сечений ребер с углами поворота нормальных элементов обшивки не линиях сопряжения обшивки и ребер записывается в виде*)

V jW - Hj

(* ) = П х , j/j) + £ <pt <x . у . ),

u r j( x ) = u r ( x r y j),

I г г <*А.у),

V / , W =

где UtJf,uf- компоненты вектора перемещений срединной поверх­ ности обшивки; Uj, Vj, Ufj, V/, f KP/ “ осевое перемещение, пере­ мещение в направлении касательной к срединной поверхности об­ шивки и прогиб j -го продольного ребра, угол поворота и угол за­

кручивания поперечного сечения этого ребра; и ; , М

W,

f :

М Ч

J1

м _

<fKpji - аналогичные характеристики / , -го кольцевого ребра; уг­ лы поворота нормали к срединной поверхности обшивки относи­ тельно координатных линий у и х выражаются через if, Uf следу­ ющим образом:'

•Принятое правило знаков соответствует случаю, когда ребра укрепле­ ны на внутренней поверхности обшивки. В отличие от 1.1 в 1.2 - 1.4 ui if, w, Vv - компоненты перемещений координатной поверхности и угла поворота нормали к ней.

h j'~

толщина оболочки и радиус ее срединной поверхности;

HJt

Hj1 - расстояния соответственно от оси j-ro продольного и

j-го кольцевого ребер до срединной поверхности обшивки.

нСоотношения (1.1) удобно переписать в виде

U jtX ) - U ( X .jr j) 4 .h / r i 1 X ,ft},

irJ (x) = tnx,y.)+hj f i {x,^.j,

u r j(* ) = w ( x , y j ) ,

urJ t ly ) = w ( x i t , y ) ,

где h j - 0,5 h+H jt ^ j1=0l5 h + H ;1, £ L ,U: -

координаты осей со­

ответственно j\ -го кольцевого и j -го продольного ребер.

В'дальнейшем применяется также следующая форма пред­

ставления соотношений (1.2):

 

2я *

 

 

и.. (X) = J

[а(х,у) + нj v,(Х,у)]8ц-у)йу,

vjlx) =

]

№ , у ) + Л}Уг {х,у)]8(и-и)йи,

urj(x)=

°гл т

*

J

unx,j/)8ty~yj)(ty;

 

 

гкг

 

VKPj(X)= J чг (х,у) % - у ; ) <fy; -

(13;

V yJ = /

^ [ V ( I , y) + h j( <f>2 (CC,y)] g ( X - X j j i l l ,

u rjW ) = j a ; l x , y ) S { x - X j ) d x ,

оL

где L - длина оболочки;

8 ( х - X j ), S (i/-u .) - дельта-функции Ди­

рака [3].

1

Деформации срединной поверхности обшивки (£ if L %,fc12) и

изменения ее кривизны (J ^ ,X 2) И кручения (х 12) выражаются че­ рез компоненты вектора перемещений по формулам [1]

£ e i f L .

е

= — + ^

(1 4 ,

1 ЭХ • * оу

г ' «

в г

д у 1

а2ш

:

шч. -у

„ Л

ьду_скл

г _*аг*'

*

' ^ а + г 1 *

1 t~ 1дхду

Г 'З у f a '

Деформации ребер выражаются через соответствующие ком­ поненты вектора перемещений их осей и углы закручивания попе­

речных сечений так:

 

 

 

€}

dU; .

 

® b E i,

■»

=

>

9Х* ’ 1

дХ

X

J

дХ

V

 

 

 

 

 

(1.5)

 

 

 

 

 

•ФЛ

 

 

 

 

д уг

г

Соотношения (1.2) - (1.5) далее использованы для получения по­ тенциальной и кинетической энергии системы.

1.2. ЭНЕРГИЯ СИСТЕМЫ. УРАВНЕНИЙ ДВИЖЕНИЯ

/равнения движения могут быть получены в результате приме­ нения принципа стационарности действия [21], согласно которому

t „

S IV = 0.

(1.6)

Здесь W - JL d t -действие по Гамильтону, где t ] t *- фиксирован-

ныв-моменты времени; LSK- П - функция Лагранжа, К и П - кине­ тическая и полная потенциальная энергия рассматриваемой си­ стемы; t - временная координата. Полная потенциальная энергия определяется как сумма потенциальных энергий деформации об­ шивки, продопьных и кольцевых ребер и потенциала внешних сил, приложенных как к обшивке, так и к ребрам, и записывается в виде

+ Х П.. + Г П , + А + 2 Л ,-*-!’ А;

<1.7)

1

*

у--» 1 у ,- '

Входящие в это выражение составляющие полной потенциальной энергии определяются по следующим формулам [1; 20]:

 

П

о

*1(Г(С-!£в +2С

«

££+С ег +С

)йхйЦ +

 

 

 

2 * { у ii

1

1 г

гг г

ев

< г '

J

 

 

 

 

о о L гКг

 

 

 

 

о

о

 

 

 

 

+г /

j

‘V

! +zv

s v

V

’ +

 

 

LZJtr

 

О

О

 

 

 

 

 

 

 

i

1

[К.А V

V

£i V

 

W

**А V

 

 

 

^

 

4

i V

‘ +rw **+iX

4 & - * » & * * *

‘’•‘и

 

 

 

 

0

 

 

 

 

 

 

 

 

 

n J'i=

г

l

 

 

 

 

 

£

^

KP; , 17 , ^ T: <t1D)

 

 

 

 

 

f ^ * V

«y \ ) w d t d y ;

 

 

«.a.

Ai = -

FjC x

 

 

 

 

 

 

 

(1. 12)

Г 2- J ttjVj+XfjtfjXtx;

 

 

 

 

 

Fj Gu

v rr

 

 

 

 

 

 

( w

 

 

 

 

-

P

Ч

’ 1

 

'1 4

 

 

 

 

 

 

 

 

 

 

-------------------- -

где в дополнение к приведенным в 1.1 приняты такие обозначения:

 

M+N

 

 

 

 

 

 

С

= Г ,

8** (6.

-S.

.),

 

 

 

4 M+N <

0

2

 

*.

 

 

 

0,',, И

у

W

1-

(11*>

 

1 " +NA

 

 

 

 

 

^

^ | , в

; ’ [«;а- б ^ - ^ ( в ^ Ч - ' )+зд(5« г М ]-

* .* , = 1 ,2 ,8 ;

5gg - толщина 3j -слоя; М - число слоев, отсчитываемых От внут­ ренней поверхности обшивки до координатной поверхности; N - число остальных слоев; Д - расстояние вдоль нормали от коорди­

натной поверхности до внутренней поверхности обшивки; h = М+ N

~ 2 8S - толщина многослойной обшивки;

 

 

в, _ 6 .

"

b f * < * “ И ? * ,* * “

" "

 

 

9

( fj* t E^21 $1* G $t - модули упругости и коэффициенты Пуас­

сона материала зг-го слоя соответственно в направлениях осей х и у и модуль сдвига этого слоя); Л у; Л ?}~ площади поперечных сечений, их моменты инерции при изгибе в радиальной плоскости, в плоскости, касательной к срединной по­ верхности обшивки, и при кручении /-го продольного и го коль­ цевого ребер соответственно; E j, Gj, E r , - модули упруго­

сти и сдвига материала этих ребер; А, ^ - числа продольных и кольцевых ребер; $Xf Gyt6j,6j- докритические напряжения в обшивке и ребрах. При решении задач устойчивости, рассматрива­ емых в этой книге, докритическое состояние принимается безмоментным; принимается также, что Gj - 6 г ,а -&у

Кинетическая энергия рассматриваемой ребристой оболочки определяется как сумма кинетических энергий обшивки, меридио­ нальных и кольцевых ребер:

 

к

(115)

 

К=К0+£ К,+ 2 К;

Здесь

1 J ,'1 1

 

Pj FJ f

r 9tii. 2

^ 0 l { ^ 4 ] d y , (1.18)

4 = ¥ J

К аг>

 

 

- плотность материала s2-ro слоя обшивки;

Pj, Pj- плотность

материалаj

и jA-го ребер;

 

J

V “ +V i ;

Д5г - расстояние от координатной поверхности до з2-го сЛ0Я* Подставляя (1.8) и (1.5) в (1.6), после элементарных преобразо­

ваний, с учетом соотношений (1.2) -

(1.5), нетрудно получить сле­

дующее вариационное уравнение:

2fjy

-

l i j r r

J{H [Хх6и 1- f y i r + X , Sui)dxay - j [ T x bu+Sx ?d+

* _

к

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