книги / Теория наведенной неоднородности и ее приложения к проблеме устойчивости пластин и оболочек
..pdfКритерии устойчивости докритического процесса деформирования 241
Вводя обозначения для изгибных кривизн и обобщенных уси лий, для ДП получим:
Дх ! г ~ Т Т * ; |
N ir J Sijdxs; |
ДМу= 1 Д о !)М »;(1 1.1.8) |
|
dXjdXi |
—Ь |
—h |
|
дп = [/(aMijAXy+Nij5^ |
^ j ds . |
(11.1.9) |
В качестве критерия выступает стационарность энергетическо го эквивалента ДП ( или стационарность вариации 8*2ДП)
5(б.2ДП) = 0 .
Теперь следует использовать выражения для ДМу1 через де формации, а затем через перемещения .
Очевидно, из (11.1.9) можно получить краевую задачу, которая определяет критический параметр как собственное значение, а при ращения (ДиД Диг1, Диз1) выступают в качестве собственных форм.
Далее для болыцей наглядности можно принять допущение о малости "почти кососимметрических" интегралов по толщине (как в главе 7, (7.1.8)), тогда для моментов:
AMjj = DjjkiAXki j
и критерий имеет вид: |
|
ДП = |Д уиД х1,Д х‘+Ни~ ^ ~ ^ ] < 1 5 = 0 . |
(11.1.10) |
Поступая,как в главе 10, нормируем поле напряжений плоской пластины Ny в состоянии S„ параметром R: Nij=RNjj°. Тогда R может определяться как собственное значение краевой задачи или из усло вия стационарности (11.1.10).
Аналог отношения Релея для этой пластины также позволяет формулировать вариационную задачу для определения минималь ного значения R:
| | ( D ijkidXk|Ax|j)ds
242 Приложение вариационных критериев к проблеме устойчивости
Выражения (11.1.10), (11.1.11) принимают узнаваемую класси ческую формулировку в упругом случае, когда есть полная эквива лентность задачи в приращениях и в полных функциях, при этом Djjk| превращается в цилиндрическую жесткость и выносится за интеграл.
Без принятого выше допущения и учета всех слагаемых в вы ражении для приращений моментов имеем:
А М у = SykiAeid + DijkiA xki .
Тогда для энергетического эквивалента получим более сложное выражение;
Подставляя вместо приращений деформаций их линеаризован ные выражения через приращения перемещений, получим энерге тический эквивалент в перемещениях в виде квадратичной формы. Находя его минимум, определяем критический уровень суммарных усилий.
При всех критериальных формах, приведенных здесь, вопрос о единственности продолжения процесса деформирования из S„ в Sn+i решается из условия соответствия текущего значения параметра внешнего процесса критическому значению, полученному из реше ния (11.1.11) или (11.1.10).
11.1.2.Энергетический эквивалент продольно сжатых замкнутых цилиндрических оболочек
Рассматривается осесимметрическое деформирование круго вой цилиндрической оболочки радиуса R, продольно сжатой равно мерно распределенным по контуру усилием, принимается кинема тическая гипотеза Кирхгофа-Лява.
Тогда геометрические соотношения в приращениях имеют вид:
Аец= Дец+хз Ахи; Де12=0 ; д е22= - диз = д822 .
R
(11.1.12)
Критерии устойчивости докритического процесса деформирования 243
Пусть, наряду с исходным процессом деформирования и дегра дации, когда поперечных перемещений нет, имеет место бифурка ция процесса, и наряду с исходным продолжением {Дщ, 0, 0}, появ ляется смежное продолжение процесса, также являющееся осесим метричным {Aui + Ай/, 0, 0 + Диз1}. До состояния S„ оболочка де формировалась без прогибов.
Варьируя (линеаризуя) геометрические соотношения, записан ные для основного и смежного продолжений, получаем:
- для основного продолжения
Дец = А е п = ^ ^ ’ Д е 22 = Ае^ = 0 *,
дк\
- для смежного продолжения
Аеп= Дец+ Aeli= Дец+Деп + хзДх},;
Де22= 0+ ~Д из’ |
Де?2= ®’ |
Aen = ~ ^ U" * |
|
|
К |
|
0X1 |
Из энергетического эквивалента |
|
|
|
ДП=Г |
. 1 . 1 , |
0 Д и з ЗД и з |
d v |
AaijAejj+aij- |
|||
|
|
0Xj 0Xj |
|
получим для данного случая; |
|
|
|
АП= f Да! ^Де! 1+хз Дх| j)+ — ДсггДиз+сл |
v ' (П.1.13) |
Используя уравнения состояния, нужно перейти к перемеще ниям. Уравнения состояния для вариаций приращений исходного процесса для данного случая имеют вид:
Д а и - Ч 'и и Д е п + Ч 'и г г Д е и ; |
щ j 14 ) |
Да22 = 2211Де!1+ 4*2222Дв22 . Проинтегрируем (12.14) по толщине;
ДИп-Bi111Ле{1+ Вп22Де22 /
(1 1 .1 .1 5 )
ДЫ22= В2211Де!1+ В2222Дегг
244 Приложение вариационных критериев к проблемеустойчивости
Выразим AN22 через ANи по (11.1.15) i
ANai = aANu+ЬДеи» где введено |
(11.1.16) |
а= В ш ь=В2г!г- Ш м В щ |.
Вии Вии
Проинтегрировав в (11.1.13) по толщине, подставим в него (11.1.16). учитывая, что ANn= 0 (как в пластине, п. 11.1.1):
An=JJ A M i.A x J .+ b ^ A u ^ + N n ^ - ^ jls . (П 1.П)
Запишем выражение дм и Для этого случая: |
|
|
||
Д М п -SH II^ " ~ |
+ Sii22~9Au3 + DiinAXn |
. |
(11.1.18) |
|
|
dxi |
R |
|
|
Можно использовать упрощенные соотношения (7.1.8) главы 7. |
||||
Тогда |
|
|
|
|
|
ДМi1—Di 111 Ах}1 . |
|
(11.1.19) |
|
В этом случае для энергетического эквивалента |
1 |
|
||
ДП = [ |
|
J |
(И .1.20) |
|
Dmi(Ax‘n)2+b^j(Au!i)2+ N i ^ ~ j |
i s ., |
|||
При R-»oo (пластина при плоском напряженном состоянии, |
||||
сжатая в одном |
направлении, с суммарными усилиями только |
Nn) (11.1 -20) переходит в (11.1.10) для этого же случая нагружения. В данном случае (i 1.1.20) представляет собой вторую вариа
цию функционала ДП, полученного в главе. 7 (формула (7.1.13) с учетом всех нелинейных относительно приращений слагаемых и внешнего процесса, который здесь при варьировании "заморозился").
Критерий формулируется как условие стационарности 6ДП=0 или минимум аналога отношения Рэлея.
Критерии устойчивости докритического процесса деформирования 245
11Л.3. Энергетический эквивалент сжатых по образующей цилиндрических панелей
Рассматривается пологая цилиндрическая панель, сжатая с торцевых кромок равномерно распределенными контурными уси лиями, деформирующаяся при деградации свойств материала. Свя занный процесс описывается теорией наведенной неоднородности.
Принята кинематическая гипотеза Кирхгофа-Лява. Докритическое деформирование характеризуется отсутствием прогибов. Изу чается возможность бифуркации этого процесса, когда в критической точке возникает смежное продолжение с ненулевым приращением прогиба Айз.
Геометрические связи в приращениях для панели имеют вид:
Деу = Деу + ХзДхи» >J= 1 |
Aei3 = 0 . |
Для исходного процесса прогибов нет:
Ну 'У dxj Эх,
•Для смежного продолжения:
где
( 11.1.21)
1 _ ЭАи} . ЭАиг S12— :— + ------1
д х2 д х\
Выписываем выражение для энергетического эквивалента:
246 Приложение вариационных критериев к проблеме устойчивости
После интегрирования по толщине
АП = JJ| ANSjAe!j+ДМ уAXy+N11
Для обобщенных внутренних усилий справедливы соотноше
ния ;
ANjj Bjjk)Деы SijklAXkl
AMSj = SijkiAEICI+ DijkiAx'k|
Тогда для энергетического эквивалента
Получим выражение для энергетического эквивалента как квадратичную форму относительно приращений перемещений. Для этого после преобразований подставим выражения вариаций прира щений деформаций и изгибных кривизн через перемещения (11.1.21) в выражение для ДП (11.1.22):
ДП=JJ|в , 11I(AE! I) +4Bi 112Де'|2 Де!i+2Bi122ДЕ22Де! I+4BI212(AEI2) +
+4В,222Де12Де22 + В2222(де22)2 + Dm I(AX!I) +4Di112Дх12дх!I+
+2^Si111Ах!1Де|i+2Si 112Ах!2Aei1+S1122AX22AEI 1+2S1211Ax!jAEI2 +
+4SmiAX2,Aei2 + 2Sl222AX22AEi2+2S22IlAx!,As!2+2S22J2Ax!2AE22+
Теперь в приращениях перемещений:
Вариационные формулировки критерия бифуркации |
|
247 |
|||||||||||||
|
|
+4йш 1— |
Й 5Ди! |
|
|
5Диг |
1 |
J |
) |
&L |
|||||
лп=,^ашй5xiг ) |
+ЗЛи^ |
+2Bim |
R |
||||||||||||
|
|
|
|
5x2 |
5xi)J |
5xi |
|
|
5x2 |
|
5xi |
||||
* |
|
|
|
|
|
( 5Aul, дАиг\ (>5Диг |
|
|
|
|
|||||
|
5Aul | 5Аиг |
|
|
|
|
|
|||||||||
+ 4 В 1 2 1 2 |
5хг |
|
5xi' |
+ 4B l222\ Эхг |
5xJ |
l 5X2 R ^ I- |
|
||||||||
ЭДиг |
|
j_ |
Диз +D1111 |
|
|
|
5хДиз а2Диз |
|
|
||||||
+ В2222 |
5 х 2 |
|
R |
+4D,H2— |
--J-+ |
|
|||||||||
. |
|
|
|
|
|
|
|
5xi5x2 |
5XI |
|
|
||||
+2Di |
а2Диз 5* Дц| |
+ 4 D l2 1 2 |
|
|
|
а2Диз 5^Диз |
|
||||||||
5x2 |
|
5xi |
|
|
|
5XI5X2 |
5X2 |
|
|||||||
|
|
|
V5xi5x2/ |
|
|
|
|||||||||
+ D2222V |
5x2 |
) |
v |
Эх? |
5X1 |
|
|
5XI5X2 |
5XI |
|
|
||||
|
|
|
|
а2 диз ади! |
а2диз fад и !, адиО . |
|
|
|
|||||||
+Sll22---2— Г - +2Sl211■—2~ \ |
"T"”+ ^7“ |
+ |
|
|
|
||||||||||
|
|
|
|
5x2 |
5xi |
|
5xi |
1 5X2 |
5xi J |
|
|
|
|||
+ 4 S 122 1 |
а2 Диз f 5Auj |
5AU; +2 & а З Д ^ 1+ ^ ) + |
|
||||||||||||
|
|
aXi3X2V 5X2 |
5xi |
|
|
5X2 |
l 5X2 |
|
5xi J |
|
|||||
|
а2Диз f5AU2 i |
11. |
53Диз f 5Диг _ |
RAuS)4 |
|||||||||||
+ S 2 2 11 |
2 |
l |
5 2 |
Rлиз|+2&ш^ |
и |
||||||||||
|
5x |
|
x |
|
|
|
|
ds .
11.2. Вариационные формулировки критерия бифуркации процесса при докритическомдеформировании с изгибом
11.2.1.Критерий устойчивости для пологих оболочек
сгипотезой Кирхгофа-Лява
Рассмотрим пологие оболочки постоянной кривизны, геомет рические особенности которых обсуждаются в работе, начиная с гла вы 6, где выводилось вариационное уравнение, п.6.1, рис.67, соот ношения (6.1.1) - (6.1.7).
248 Приложение вариационных критериев к проблемеустойчивости
Согласно гипотезе Кирхгофа-Лява, для приращений деформа
ций
Деу = Деу + х3Дху, ij= 1,2; Де*з = 0 . |
( 11.2. 1) |
Для деформаций срединной поверхности Деу выполнены (7.2.2) главы 7.
Запишем энергетический эквивалент для пологой оболочки для состояния Sn, считая, что процесс деформирования и деградации имеет бифуркацию, в результате которой появляется смежное с ос новным продолжение процесса, отличающееся от основного продол жения на бесконечно малую вариацию, не нарушающую граничных условий, обозначаемого в этой главе верхним индексом
Для энергетического эквивалента
( 11.2.2)
для данного случая, после интегрирования по толщине, получим
(11.2.3)
Здесь нужно брать линеаризованные зависимости для прира щений, имеющих смысл вариаций.
Для приращений внутренних усилий используются выражения в вариациях (без неоднородности):
A N jj-B jjk iA sici+ S ijk i^ X k i ;
ДМу ~SijklДбк1 Djjkl^Xkl *
Функционал для пологой оболочки обсуждался в главе 7. Мож но сравнить (7.2.10), (7.2,12) с (11.2.3).
Подставляя в (11.2.4) линеаризованные выражения д е‘к|,
Вариационные формулировки критерия бифуркации |
249 |
I _ 1 ( Деуи - -
2 \
|
|
|
д х \д х \ - |
(11.2.5) |
3Aui, 3AuJ |
|
|
||
+ |
ЗДиз Зиз, Зиз ЗДиГ) . |
|||
----- + |
------3xj |
---------- +---------- |
-к;:" kijАйз8ij |
|
3xj |
|
3xi 3xj 3xj |
3xj |
и подставляя полученные выражения в (11.2.3)^ получаем энергети ческий эквивалент (в вариациях перемещений), минимизация кото рого дает критерий бифуркации:
= f j(Bi-»k l ^ e ij+ 2SijkiAXkiAeSj+ DijkiА Х ы AXy+ N jj ^ 3 |
Jds . |
( 11.2.6)
Иногда можно пренебрегать "почти кососимметричными" ин тегралами по толщине Sjjki, тогда (11.2.6) упрощается и форма кри терия становится более похожа на формулировку в упругости, когда рассматривается бифуркация состояния равновесия в полных функ циях (если считать тензоры Вуи, Dpi не функциями напряженнодеформированного состояния и свойств деградации, а константами материала).
Накопленный уровень усилий (напряжений) Ny на текущем ша ге инкрементальной теории должен нормироваться неизвестным па раметром, который и определяется из условия стационарности АП. Значение параметра проверяется на соответствие ведущему пара метру внешнего процесса, после чего делается вывод о том, проис ходит ли бифуркация в состоянии Sn.
1 1 .2 .2 .Критерии устойчивости для пологих оболочек.
Геометрическая модель с учетом поперечных сдвигов
В отличие от предыдущего п.11.2.1, здесь усложняются геомет рические инкрементальные соотношения, более сложный вид приоб ретают соотношения между обобщенными усилиями и деформа циями в приращениях. Эти соотношения для рассматриваемой здесь уточненной кинематической модели были получены в главе 6.
Получим выражение для энергетического эквивалента, стацио нарная точка которого определит возможность существования по-
250 Приложение вариационных критериев к проблеме устойчивости
бочного продолжения процесса деформирования и деградации в теории наведённой неоднородности.
Для энергетического эквивалента
д п = ш |
. 1 а I I |
дДиз дАиз |
dv |
(11.2.7) |
|
ДацДец+ац— -----— |
|||||
|
|
OXi |
OXj |
|
|
Для Д(у|, д е{| принимаются |
варьированные зависимости. Для |
Деу нужно принять линеаризованные (6.1.6),(6.1.7) главы 6. Нели
нейные деформации срединной поверхности (6.1.7) переходят в (11.2.5), а (6.1.6) линейны по приращениям пяти независимых пере мещений: д и‘т , т = 1,2,3, Дуь i=l,2. Напомним, что индекс "°" ис
пользуется в главе 6 для обозначения перемещений в срединной по верхности. Здесь этот индекс будет опущен, все использованные обозначения перемещений будут означать перемещения срединной поверхности. В отличие от п.11.2.1 в (11.2.7) будут входить де':,.
Подставляя в (11.2.7) |
выражения |
(6.1.6) |
и |
интегрируя в |
(11.2.7) по толщине, получим: |
|
|
|
|
ДП = 1 ДМуД8у+(дм{—ДМJ1)-— |
AQJ *)Де!з~ |
|||
,, |
«1д2Диз |
адизадиз |
( 11.2. 8) |
|
ds . |
||||
|
axjaxj |
axi |
axj |
|
Обозначение "" принято здесь для обозначения обобщенных усилий высших порядков сдвиговой модели, для них приведены оп ределения в (6.5.3).
Для обобщенных усилий в (6.3.3) - (6.3.5) получены выражения через приращения деформаций и перемещений. Эти выражения сле дует записать в вариациях, то есть принять для входящих в них при ращений деформаций линеаризованные зависимости и отбросить не однородность в (6.3.3) - (6.3.5).
Как видно, из (11.2.8) получится квадратичная форма относи тельно приращений перемещений.
Запишем частный случай, полагая нулю "почти кососиммет ричные" жесткостные интегралы по толщине в выражениях (6.3.3)- - (6.3.5), которые определены в (6.3.6) обозначениями Syu, Syu*. Sijki**. Тогда из (6.3.3) - (6.3.5) следуют