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

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

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

Уравнения состояния с учетом неравномерного

91

4.4. Уравнения состояния с учетом неравномерного, нестационарного температурного поля

Полагаем, что при нагреве тепловое расширение етувеличивает среднее упругое удлинение, связанное со средним напряжением за­ коном Гука:

(4.4.1)

eT= a (T)T(x,y,z,t).

(4.4.2)

Здесь а(Т) — коэффициент теплового расширения; T(x,y,z,t) — распределение температуры по объему тела в момент времени t; е" - полная средняя деформация, являющаяся суммой температурного расширения и упругой деформации от напряжений.

При условии несжимаемости материала К=оо для компонент пол­ ной деформации выполняется:

(4.4.3)

Отсюда получим:

е," = eij+ ET6ij

(4.4.4)

Уравнения состояния материала при этом имеют вид:

Получим уравнения состояния в приращениях:

92

Линеаризованные уравнения состояния

Представляя приращение ДА. в виде суммы ДА? и ДАЬ получим:

Д\б=а,ДДе|-Дет5(|);

(4.4.7)

дх, - » ( е : - е «) ;

(4.4.8)

_ 2 С Й -Е с*

'

(4.4.9)

и “ _ ьс

2

sii

J

CTi

 

 

 

Тогда

 

 

 

 

As,j= ^ „ ( д е и

- Д

^ ^ ^ s , j .

П А Ю )

где Ар определяется выражением (4.1.15).

С учетом сжимаемости материала для компонент полной де­

формации справедливо:

 

 

 

” е88«=е+Ет

 

 

(4-4.11)

Представляя приращение ДА в виде суммы дх? и дх* для дх?>

получим:

 

 

 

ДА?= ДАу(де|]-Дет8^

 

(4.4.12)

Здесь да|| определяется выражением (4.2.10), а приращение

дхГ ‘ выражением (4.2.11).

 

 

 

Тогда будем иметь:

 

 

 

Asy= Х.')к|(де"г(де+ДЕ,)бк|)+

syД °+ ^

s,j . (4.4.13)

Связь между приращениями

компонент тензора

напряжений

и тензора деформаций с учетом

неравномерного нестационарного

температурного поля будет:

 

 

 

Уравнения состояния с учетом сжимаемости

93

(4.4.14)

В выражения для x.*jki и &к* входят, в данном случае, парамет­

ры E J и Ек. получающие возмущение вследствие изменения темпе­ ратуры во времени.

Аналогичным образом получаем связь между приращениями компонент тензора напряжений и тензора деформаций с учетом не­ равномерного нестационарного поля при изменении объемного мо­ дуля упругости:

Летij= ^k](Aek|-Дет5к,)+ ^ j^ A o 5ij+ ^ S ij+ ^ ~ - r a 5ij •

(4.4.15) В данном случае объемный модуль К получает возмущение в

том числе и вследствие изменения температуры во времени.

4.5.Уравнения состояния с учетом сжимаемости при плоском напряженном состоянии

Рассмотрим тонкостенную конструкцию в декартовой систе­ ме координат Х1Х2Х3, причем ось хз направим перпендикулярно сре­ динной поверхности. Из условия обобщенного плоского напряжен­ ного состояния следует, что о 3з=0, е зз* 0 - Среднее напряжение в этом случае равно

Pl 1+022

(4.5.1)

Для е3з из условия плоского напряженного состояния следует:

e.i3= " 7^ e~(ei 1+622) •

(4.5.2)

1-Рс

 

94

Линеаризованные уравнения состояния

 

Учитывая (4.5.2), получим для е; выражение :

 

 

_ i / f t

eH+en+^ e u e 22+— (е?2+е?з+е2п)

(4.5.3)

 

 

|5"

^

 

 

ft

ft

.

 

Здесь введены следующие обозначения:

 

 

 

 

 

f t =1- f t + ^c'

 

 

 

 

 

ц2=4цс- ^ - 1 ;

 

(4.5.4)

 

 

 

f t =3(b

ft)J-

 

 

При построении

соотношений в приращениях,

при условии,

что ДК = 0, будем исходить из уравнения:

 

 

 

Деи

Доп

ЗК-Х

е

 

(4.5.5)

 

"

Ж

Дсг5'г

 

 

 

 

 

Представляя, как и ранее, ДХ в виде суммы двух слагаемых ДХГ и дх,**Фи выражая их через возмущенные значения параметров диаграммы деформирования Ек* и ЕЛ получим:

А .и

.с .

 

(4.5.6)

ДХй

 

 

д х Г ^ г + т ' Л Д

-

(4.5.7)

 

 

Ео

D

 

У_ Ес

Ес

Ес*~Ес Ес

 

1+^с

(1+Цс)2Ео

 

/ - E I - E C

Ек~Ес Ее .

 

l+ ft

(l+Hc)2ЕоV

 

В =

 

 

 

+е|5).

=^-М с)м се?+(2цс-1)(е?1+е!2)+ 2(2-цс)е)1е22.б (1 -Цс)(е?г+^

Уравнения состояния с учетом сжимаемости

95

С=

ei1+Ц2ем) Aei1+ (2ц.

P2ei

1) Дегг+2p3(ei2Aei2+ei3Aei3+e23Дегз);

 

D = 2( l - H Jef

-B v 0^

.

х' Ео

Подставляя (4.5.6) и (4.5.7) в (4.5.5), после преобразования получим:

 

 

+

 

+ ^ ’ SS

(4.5.8)

Здесь Хм

представляют следующие выражения:

 

Л.ГТ= ~~ (^Pjei 1+ 1^

22); Хп= "^(2м-[е22+Р2е11)’

 

X*2 = A.2*

= ^ 3e j2 ^

Л.1з = Хз|в Ц-3емр*

А.гз^Хз^Рзегз” ’

• (4.5.9)

Вводя обозначение

 

 

 

 

 

a.№,= «>*Sj.+ ^ S i j

 

(4.5.10)

получим:

 

 

 

 

 

 

Asij= Xji Д§к1 +

SijAo + ~ sij

(4.5.11)

Переходя от компонент приращений девиаторов напряжений и деформаций к компонентам приращений тензоров напряжений и деформаций, получим:

Д сгв«х5,Д еы + ' ^ Л с 8» + ^ Ч « • (4.5.12)

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

аС -

96 Линеаризованные уравнения состояния

Учет температурных деформаций приводит,

по аналогии с

(4.4.14), к следующему выражению:

 

 

 

Д с г > "к,(леГ,-де'6 к|) +

Д<т8у +

sij •

(4.f 13)

Здесь так же,

как и в (4.4.14),

в выражения

для x.Jki и Д*.Г

входят параметры

е * и Ек, получающие возмущение

вследствие

изменения температуры.

4.6.Уравнения состояния с учетом сжимаемости

иизменения объемного модуля упругости при плоском напряженном состоянии

Полученные соотношения можно обобщить для случая изме­ нения объемного модуля упругости материала, то есть для ДК * О

в (2.1.15). Используем полученные ранее выражения для АХ

(4.3.4)

и Дрс (4.3.5), учитывающие изменение К.

 

Вычислим приращение интенсивности деформаций в

(4.3.4)

и (4.3.5) для случая плоского напряженного состояния при условии, что

Де; = ЕсК*-К.Ес В_ + _С

(4.6.1)

е; 6К2 D* D*

Здесь В и С имеют тот же смысл, что и в выражениях (4.5.6) и (4.5.7), а выражение для D* имеет следующий вид:

0 * = 2 (Н » !)е? + В & ^

(4.6.2)

Представляя Д*. в виде суммы двух слагаемых д ^ Ф* и дх,**\ в данном случае будем иметь:

Уравнения состояния с учетом сжимаемости

91

 

 

С

 

(4.6.3)

АХ?, =у

 

6К2

D

(4.6.4)

 

Здесь

*Ес~Ес + __Ес__ ЕсК ~ЕсК

1+Мс (1+цс)2 6К2

Подставляя (4.6.3) и (4.6.4) в выражение (4.3.1), после преобра­ зования получим:

 

LX.ki Ski _

ч +

^

о В , • (4-6.5)

 

з к х

 

 

 

з к 2

Здесь

представляет следующее выражение:

 

A.iT= |4(2p,ei 1+2р2е2г);

Хп ~ ^ *(-И-

i

eiО’

 

Л,Гг*= A.2I* М^зе^ ;

 

^13 = ^31

 

(4.6.6)

 

 

 

 

Хгз ~Л.32 0 *Кзегз

Вводя обозначение

(4.6.7)

Xijkl" ^8ik8jl+ ^ SU’

получим:

98

Линеаризованные уравнения состояния

 

 

 

Тогда компоненты приращений тензоров напряжений

и де­

формаций будут связаны зависимостью:

 

 

 

Atfip А.SI Деи +

Да5у+ ~ ~ ц +х ~ г г о5ч

(4-6-9)

Учет температурного воздействия приводит к выражению:

 

Доу=Хр(декг AeT6ki)+

Дстбу+ ~~~ Sij+A. -" —

а §у .

(4.6.10)

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

4.7. Эквивалентные формулировки уравнений состояния. Уравнения состояния в инкрементальной форме

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

Процесс деформирования твердого тела будем описывать по­ следовательностью состояний равновесия S0, Si..... Sn, Sn+i3..., Sr, где So - исходное состояние твердого тела. Обычно в нем нет напря­ жений и деформаций, нет и приложенных внешних воздействий; Sr конечное состояние, при котором требуется определить напря­ женно-деформированное состояние, исследовать устойчивость или решить другие задачи; Sn-произвольное промежуточное равновесное состояние. В состоянии Sn известны все параметры внешнего и внутреннего процесса: уровень приложенных нагрузок, уровень дег­ радации физико-механических свойств материала, тензоры напряже­ ний и деформаций в точках тела, вектор перемещения каждой точки

 

Эквивалентные формулировки уравнений состояния

99

объема.

Эти параметры процесса,

или переменные состояния,

из­

вестны

не только в состоянии равновесия S„, но и во всех прой­

денных за историю процесса состояниях Si,...,Sn-i.

Каждое из со­

стояний Sn будем иногда называть шагом процесса.

 

 

Построение инкрементальной теории состоит в

формулирова­

нии инкрементальных уравнений

(уравнений в приращениях или

скоростях), позволяющих получить решение для приращений на ша­ ге Sn+i, чтобы с помощью переменных состояния Sn и найденных приращений шага S„+i получить переменные состояния S„+i.

При формулировании теории будем исходить из того, что шаги S„ и Sn+i бесконечно близки друг к другу, а состояний при переходе из So в Sr бесконечно много. Следует отличать построение инкре­ ментальной теории от организации дискретного алгоритма расчета, когда следует интересоваться еще и проблемами сходимости, точно­ сти, устойчивости алгоритма. Алгоритм строится на основе сфор­ мулированной теории, где решены принципиальные вопросы обос­ нования ее соотношений, математически промоделированы физиче­ ские процессы, определены области применимости моделей.

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

Таким образом, в дальнейшем считаем, что шаги процесса не­ прерывно следуют друг за другом, а величины приращений, отде­ ляющих состояние S„ от состояния Sn+i, настолько малы, насколько это требуется для математических построений.

Под процессом деформирования So,..., S„, S„+i,...,Sf будем по­ нимать единый взаимосвязанный процесс деформирования и изме­ нения физико-механических свойств материала в точке тела, зави­ сящий от напряженно-деформированного состояния, которое, в свою очередь, зависит от изменяющихся свойств материала. В со­ стоянии So будем считать материал однородным, изотропным и фи­ зически нелинейным, то есть зависимость компонент тензора на­ пряжений от компонент тензора деформаций может быть отражена следующей формулой:

о# - огу(ец,е|2 ,...,езз)

ixj - 1,2,3 .

(4.7.1)

100

Линеаризованные уравнения состояния

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

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

Aaij = Аоу(Деп,...,Аев,е|1,...,езз,оц,...сз?,...).

(4.7.2)

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

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

A sy = УрА ек1 +фу или А ец =FijkiA ski+fij .

(4.7.3)

Здесь Ч'р, фу, Fyki, fjj являются функциями переменных состояний. Приращения в инкрементальных теориях отсчитываются от текущих значений переменных состояния.

В предыдущем параграфе обсуждался аппарат теории наведен­ ной неоднородности, с помощью которого учитываются изменения, происходящие в материале нагруженной, конструкции в результате внешних воздействий. Таким аппаратом являются функции дегра­ дации Ф\ (i=l,2,...,m), вводимые в диаграмму деформирования мате­ риала <Ti(ej, Ф|,..., Фт). Для функций деградации имеется система оп­ ределяющих уравнений в приращениях, опирающаяся на экспери­ ментальные данные о понижении показателя объективной прочно­ сти, о смещении кривой длительной прочности, об изменении диа­ грамм деформирования материала в результате увеличения плотно­ сти микроповреждений.