книги / Теория наведенной неоднородности и ее приложения к проблеме устойчивости пластин и оболочек
..pdfУравнения состояния с учетом неравномерного |
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, Ф|,..., Фт). Для функций деградации имеется система оп ределяющих уравнений в приращениях, опирающаяся на экспери ментальные данные о понижении показателя объективной прочно сти, о смещении кривой длительной прочности, об изменении диа грамм деформирования материала в результате увеличения плотно сти микроповреждений.