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

книги / Теория наведенной неоднородности и ее приложения к проблеме устойчивости пластин и оболочек

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

Критерии устойчивости докритического процесса деформирования 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) следуют