книги / Оптимизация подкрепленных цилиндрических оболочек
..pdfББК 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) |
|
|
|
|
|
V»•ФЛ |
|
|
|
|
д уг ‘ |
г |
Соотношения (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+
* _ |
к |