книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред
.pdf§ 2.4. Вариационные постановки задач |
81 |
|
J <х(е(и)) • • |
е(6и) dV —J pl\ (е(6и)) dV = |
А е, |
V |
V |
(2.4.60) |
-J l l(s(u))5pdV = 0. v
Вариационная постановка квазистатической задачи МДТТ для несжимае мой среды заключается в нахождении системы функций (и,р), удовлетворя ющих системе (2.4.60).
2 .4 .7 . Вариационная постановка задачи в напряжениях
Рассмотрим квазистатическую задачу (2.3.17) в напряжениях. Будем на зывать поле тензора сг(х) статически допустимым, если оно является непрерывно-дифференцируемым в У и удовлетворяет статическим граничным условиям (2.3.17г) и уравнениям равновесия (2.3.17а). Если же поле сг(х) удовлетворяет всем уравнениям (2.3.17), т. е. кроме того еще и уравнениям совместности (2.3.176), выраженным через напряжения, то такое поле будем называть действительным.
Запишем функционал, называемый кастильянианом:
/С (сг) |
Z dV, |
(2.4.61) |
V
где Z = р £(сг,во) — потенциал (2.1.48). Тогда, если рассматривать дважды непрерывно-дифференцируемое поле сг(х), имеет место следующая теорема. Теорема 2.4.6 (вариационный принцип Кастильяно). Пусть тензорный оператор (2.3.17в) является потенциальным:
в = П (а) = dZjdcr, |
(2.4.62) |
тогда среди всех статически допустимых полей сг(х) действительное поле отличается тем, что для него и только для него кастильяниан имеет стационарное значение:
6JC{cr) = 0. |
(2.4.63) |
Подставляя (2.4.61) в (2.4.63), получаем вариационное уравнение для данного принципа:
8К,= |
SZ dV = |
' dZ |
•• 5а dV = |
s(a) • • да dV = 0. |
(2.4.64) |
|
|
да |
|
|
|
V |
|
V |
|
V |
|
▼ Пусть статически допустимое поле <т(х) является действительным, тогда оно удовлетворяет всей системе (2.3.17), т. е. существует поле деформа ции е(а(х)), удовлетворяющее уравнениям совместности, а, следовательно, существует и вектор перемещений и(х), удовлетворяющий этим уравнениям (в соответствии с теоремой 2.3.2):
82 |
Глава 2. Упругие среды с малыми деформациями |
е(а) = def и.
Рассмотрим дополнительно еще два произвольных статически допустимых поля Ti(x) и Т 2(х), удовлетворяющих уравнениям
|
V • Т а + р f —0, х G П; |
|
(2.4.65) |
|
|
И *T Q| ^ —t пе> |
^ — 1, 2. |
|
(2.4.66) |
Домножим (2.4.65) на вектор и и проинтегрируем по V: |
|
|||
11• V • Т а dV -j- |
р f • u dV = |
0. |
(2.4.67) |
|
У |
У |
|
|
|
Преобразуя первый интеграл, согласно формуле (2.4.3), получаем |
|
|||
T Q • • е(а ) с/П -Ь |
tne *и <in + р f • u (in = |
0, се = 1, 2. |
(2.4.68) |
|
у |
у |
|
|
|
Если определить вариацию 5а как разность двух статически допустимых полей 6<т= Т х - Т 2, то, вычитая уравнение (2.4.68) при се = 2 из (2.4.68) при се = 1, получаем вариационное уравнение (2.4.64).
Докажем теорему в обратную сторону. Пусть выполняется уравнение (2.4.64) для статически допустимого поля а(х). Покажем, что при этом выполняются уравнение совместности (2.3.176) и определяющие соотношения (2.3.17в).
Поскольку поле а(х) удовлетворяет уравнениям (2.3.17а) и (2.3.17г) по
условию, то можно образовать новый функционал |
|
К! {а) = /С(<т) + w • (’V • о- + °р f) dV 9 |
(2.4.69) |
v |
|
вариация которого, очевидно, также равна нулю: 5К! = 0, a w(x) — некоторое достаточно гладкое векторное поле, определенное в П.
Преобразуя второй интеграл в (2.4.69) по формуле (2.4.3), получаем
К! = К,+ |
• w - |
• • def w dV + р f • w d n |
(2.4.70) |
|
у |
у |
|
Вычислим вариацию этого функционала, принимая во внимание, что второй и четвертый интегралы в правой части (2.4.70) не содержат а и, следовательно, их вариация равна нулю:
0 = 5К! = |
е(а) • • да dV — 5а • • def w dV = |
(е(а) —def w) • • 5а dV. |
у |
у |
у |
|
|
(2.4.71) |
§ 2.4. Вариационные постановки задач |
83 |
Отсюда следует, что, в силу произвольности 5сг, должно выполняться соотно
шение |
|
е(<т) = def w. |
(2.4.72) |
Иначе говоря, существует векторное поле w(e(cr(x))), удовлетворяющее урав нениям (2.4.72) с заданным тензором в, который находят по тензору сг с помощью определяющих соотношений в(сг) = П(сг). Тогда w(x) — истинное поле перемещений, поскольку, в силу теоремы 2.3.2, из уравнений (2.4.72) следует, что выполняются уравнения совместности (2.3.176).
Таким образом, найдены поля е(сг) и u = w(e(cr(x))), удовлетворяющие всей системе (2.3.17), значит сг(х) и есть действительное поле перемеще ний. А
Вариационная постановка задачи в напряжениях (2.3.17) заключается в нахождении статически допустимого поля <т(х), удовлетворяющего вариаци онному уравнению (2.4.64). Следует отметить, что эта постановка крайне редко используется при решении задач МСС, поскольку ее решение ищется
вклассе статически допустимых функций. Иными словами, если в вариацион ных принципах Лагранжа и Хеллингера — Рейсснера решение и подбирается из класса функций, удовлетворяющих лишь условиям (2.4.1а) на границе, то
впринципе Кастильяно — в классе функций сг, удовлетворяющих уравне ниям равновесия и статическим граничным условиям (2.4.65), (2.4.66), что является существенно более сложной задачей.
2.4.8. Условие устойчивости Адамара
Важную роль в механике деформируемых сред играет теорема единствен ности решения квазистатической задачи. Для ее формулировки необходимо ввести специальный класс твердых сред (по-прежнему рассматриваем изотер мические процессы, хотя для неизотермических можно ввести аналогичное определение).
Определение 2.4.1. Твердую среду с малыми деформациями, описываемую ф-моделъю (2.1.40а), называют с в е р х у с т о й ч и в о й по А дамару, если для любых несовпадающих пар полей вДх) и е2(х), определенных в области V, выполняется неравенство Адамара (условие эллиптичности):
((т(в2) - <x(ei)) • • (е2 - e\)dV > 0. |
(2.4.73) |
у
В случае, если имеет место нестрогое неравенство (2.4.73), то среду называют устойчивой по Адамару.
Примером сред, сверхустойчивых по Адамару, являются линейно-упругие среды (см. § 2.6).
84 |
Глава 2. Упругие среды с малыми деформациями |
Определение 2.4.2. Твердую среду с малыми деформациями, описыва емую Q-моделью (2.1.46), (2.1.47), называют с в е р х у с т о й ч и в о й по А д а м а р у , если для любых несовпадающих пар полей сгДх) и сг2(х), определенных в области V (произвольной), выполняется неравенство
(cr2 - c ri) -- (е2(ег2) - £ \ { a \ ) ) d V > 0. |
(2.4.74) |
v
Для идеальных твердых сред условие устойчивости Адамара (2.4.73) удобно записать в иной форме.
Определение 2.4.3. Идеальную твердую среду с малыми деформациями, для которой ф-моделъ определяющих соотношений (2.1.40) удовлетворяет условию
h .. М Д .. h > mh • • h |
Vh, |
(2.4.75) |
де |
|
|
где т — некоторое положительное число, |
т > 0; h |
— произвольный |
ненулевой симметричный тензор второго ранга, называют и д е а л ь н о й у п р у г о й с ре д о й с п о л о ж и т е л ь н ы м к а с а т е л ь н ы м м о д ул е м .
Здесь да/де — тензор четвертого ранга, представляющий собой производную тензорной функции по тензорному аргументу (см. т. 1, § 4.13) (предполагаем, что такая производная существует).
С учетом свойства потенциальности (2.1.40д), условие (2.4.75) можно
также представить в виде |
|
д2Ч! |
(2.4.76) |
h • • - —— • • h ^ mh • • h Vh, |
|
де де |
|
где использован тензор второй производной от скалярной функции тензорного аргумента (см. т. 1, §4.13).
Определение 2.4.4. Идеальную твердую среду с малыми деформациями, для которой (-модель (2.1.46)—(2.1.48) удовлетворяет условию
—ph • • |
д ^ • • h = h • • |
• • h |
^ mh • • h |
Vh, |
(2.4.77) |
|
да да |
da |
|
|
|
где m — некоторое положительное число, m > |
0; h — |
произволь |
|||
ный ненулевой |
симметричный |
тензор |
второго |
ранга, |
называют |
и д е а л ь н о й у п р у г о й с р ед ой с п о л о ж и т е л ь н о й к а с а т е л ь н о й п о д а т л и в о с т ь ю .
Для дальнейшего введем некоторую вспомогательную функцию
№ = o-(e(ui + £w )) • • e(w ) dV, |
1, |
(2.4.78) |
v
где w (x ) — некоторое фиксированное поле; u i( x ) — поле перемещений. В силу линейности оператора малых деформаций, имеем
e (u i + £w) = e (u i) + £e(w ). |
(2.4.79) |
§ 2.4. Вариационные постановки задач |
85 |
Поскольку по предположению тензорная функция а{е) является диффе
ренцируемой, то существует и производная |
|
|
|
||
/ ч |
да |
9e(ui + £w) |
|
|
|
</Щ) = £(W) |
. . - |
--------- --------dV — |
|
|
|
V |
|
r)rr |
|
|
|
|
|
|
|
(2.4.80) |
|
|
|
e(w) • • — (e(ui + £w)) • • e(w) dV, |
|||
|
|
у |
|
|
|
которую предполагаем непрерывной функцией |
на всем отрезке [0, 1]. Для |
||||
такой функции <р(£) имеет место представление |
|
|
|||
|
|
¥>(0 = ¥>(0) + |
Г € [0 ,1 ], |
|
|
и, следовательно, |
|
|
|
|
|
|
^ (1) - ^ ( 0) = ^ (Г ), |
o ^ |
r ^ 1. |
(2.4.81) |
Выберем в качестве функции w разность двух полей перемещений w =
= U2 —ui, тогда из (2.4.78) имеем |
|
|
|
|
Щ ) = |
(т (£2) • • (е2 - £ \ ) dV, |
<р(0) = |
<r(ei) •• (e2 - e i ) dV. |
(2.4.82) |
У |
|
V |
|
|
Подставляя выражения (2.4.82) и (2.4.80) в (2.4.81), получаем |
|
|||
(а(е2) - |
<т(е1)) • • (е2 - ei) dV = |
e(u2 - |
uj) • • — (e(uj + |
|
v |
|
v |
|
|
+ <^(u2 —ui))) • • s i (112 —ui) dV ^ m |
e(w) • • e(w) dV > 0. |
(2.4.83) |
||
|
|
|
у |
|
Таким образом, доказана следующая теорема.
Теорема 2.4.7. Для идеально-упругой среды, для которой функция (2.4.78) является непрерывно-дифференцируемой на [0, 1], условие устой чивости Адамара (2.4.73) выполняется тогда и только тогда, когда среда обладает положительным касательным модулем, т. е. выполняется нера венство (2.4.74).
2.4.9. Теорема единственности для квазистатических задач
Теорема 2.4.8. Для твердой среды с малыми деформациями, описываемой ф-моделъю и сверхустойчивой по Адамару, решение квазистатической задачи (2.2.9) является единственным.
▼ Пусть имеется два различных решения задачи (2.2.9): иДх) и U2(x), х G V. Им соответствуют разные поля деформаций вДиД = def ui и £2(щ) = = def U2 и поля напряжений cr\ = <r(ei), а 2 — <х(£2)-
86 Глава 2. Упругие среды с малыми деформациями
Поскольку каждый набор (ua , s a, cra) при а = 1,2 является решением, то
он удовлетворяет интегральному соотношению (2.4.4): |
|
|
|
сга • • s(w) dV = |
w • t ne dT, + p f • w сПЛ a |
= 1, 2, |
(2.4.84) |
у |
у |
|
|
для всякого w, удовлетворяющего (2.4.1). |
|
|
|
Вычитая соотношение (2.4.75) при се = 1 из (2.4.75) |
при а = 2, получаем |
||
тождество Кирхгофа: |
|
|
|
(<х(е2) —<x(ei)) **£(w ) dV = 0. |
|
(2.4.85) |
у
Выберем в качестве поля w разность решений: w = 112 —ui. Очевидно, что при этом условие (2.4.1) будет выполнено, тогда, в силу линейности оператора def, получаем
e(w) = е(u2 - |
ui) = e(u2) - e(ui) = е2 - £\. |
(2.4.86) |
Подставляя (2.4.86) в (2.4.85), приходим к соотношению |
|
|
((т(в2) - |
<x(ei)) • • (е2 - ei) dV = 0, |
(2.4.87) |
у |
|
|
которое, очевидно, противоречит условию (2.4.73) сверхустойчивости среды по Адамару и может выполняться, только если
e(u2) = e(ui) Vx G V. |
(2.4.88) |
Из (2.4.88) следует, что при этом должны совпадать и поля напряжений сг2 и (Ть а поля перемещений и2(х) и ui(x), в соответствии с формулой Чезаро (2.3.49), могут отличаться только на смещение твердого тела как жесткого целого за счет выбора векторов UQ и OJQ. Однако поскольку оба
поля и2 и ui |
удовлетворяют граничным условиям (2.2.9д), они совпадают на |
|
поверхности |
и имеет место тождество |
|
0 = U 2 ls„ _ U I IE U = U 0 2 - u01 + ( ^ 02- ^ 01) |
X ( х - х 0) | Ец V x e E u , (2 .4 .8 9 ) |
|
которое справедливо только, если UQ2 = UQI и CJQ2 = ^01 • Таким образом, поля |
||
перемещений и2(х) и ui(x) тоже совпадают. |
А |
Теорема 2.4.9. Для несжимаемой твердой среды с малыми деформация ми, описываемой ф-моделъю (2.2.12) и сверхустойчивой по Адамару, т. е. удовлетворяющей условию (2.4.73), решение квазистатической задачи
(2.2.13) является единственным. |
|
|
|
▼ Пусть |
имеется два различных решения задачи |
(2.2.13): иа (х), |
_ра (х) |
(а = 1, 2). |
Тогда тождество Кирхгофа (2.4.85) для |
несжимаемой |
среды |
принимает вид
|
§ 2.4. Вариационные постановки задач |
|
87 |
|
(^ (е 2) - |
^ (e i)) •• e(w) dV |
2 |
dV = О. |
(2.4.90) |
(.Р -p i)/i(e (w )) |
|
|||
v |
|
v |
|
|
Если выбрать в качестве поля w разность двух решений: w = 112 —ui, то
в силу несжимаемости среды |
|
I i(e(w)) = Ii(e(u2)) - Ii(e(ui)) = О, |
|
и тождество (2.4.90) примет вид |
|
(F(e2) - F ( e x)) • • (в2 - вД dV = 0. |
(2.4.91) |
у
В то же время неравенство Адамара (2.4.73), если е2 и е\ не произвольны, а являются решениями задачи (2.2.13), указывает на то, что должно выпол няться иное соотношение:
р Ч е 2) - ^Р{в\) - (pa - Р1)Е) • • (е2 - |
eO dV = |
V |
|
( £ (е2) - |
:F(ei)) ■■{e2 - e x) d V > 0. (2.4.92) |
|
V |
|
Следовательно, соотношение (2.4.91) может выполняться, |
только если е2 = |
|
= вь а, следовательно, и |
= Т ( е i). Тогда соотношение (2.4.90) прини |
|
мает вид |
|
|
|
(Р2 -P i)^i(^(w )) dV = 0 |
(2.4.93) |
|
у |
|
и имеет место для произвольной скалярной функции Ji(e(w)), что возможно, только если р2 = р\.
Дальнейший ход доказательства о единственности поля перемещений и(х) такой же, как и в теореме 2.4.8. А
2.4.10. Необходимые условия разрешимости задачи в напряжениях
Отметим, что при доказательстве теоремы единственности существенную роль играли граничные условия для перемещений на части Т,и. Если же на всей поверхности тела £ заданы только статические граничные условия, то приходим к постановке (2.3.17) задачи в напряжениях. Очевидно, что для нее тоже будет иметь место теорема единственности, но только для напряжений и деформаций, а перемещения будут определяться по формуле Чезаро только с точностью до движения тела как жесткого целого (векторы UQ и и) останутся неопределенными). Для исключения жесткого движения к постановке в напряжениях следует добавить два условия для перемещений — обычно это условия закрепления какой-либо точки XQ G Е:
Глава 2. Упругие среды с малыми деформациями
х = XQ: и = О, V х и = 0. |
(2.4.94) |
Кроме того, задачу в напряжениях (2.3.17) или соответствующую ей задачу (2.3.22) можно разрешить только при выполнении определенных условий для входных данных задачи — векторов f и t ne.
Теорема 2.4.10. Необходимыми условиями разрешимости квазистатической задачи (2.3.17) в напряжениях являются условия обращения в нуль главного вектора V и главного момента М всех действующих на тело сил:
|
V = О, |
М = 0, |
(2.4.95) |
где |
|
|
|
V = p i d V + t ne dE, |
ЛА — |
p x x f dV + X X tne dE. |
(2.4.96) |
у |
У |
|
|
▼Действительно, пусть задача (2.3.17) имеет решение. Тогда, проинте грировав уравнение равновесия (2.3.17а) и использовав формулу Гаусса — Остроградского (т. 1, (3.5.13)), а также граничное условие (2.3.17г), получим
О |
V • О- dV + p i d V = |
t ne dE + pfdV . |
(2.4.97) |
у |
у |
У |
|
Если же предварительно векторно умножить (2.3.17а) на х, а затем
проинтегрировать по V, |
имеем |
|
|
|
О = (V • а) х х dV + |
р f х х |
dV = |
V • (cr х х) dV + |
р f х х dV = |
у |
у |
у |
у |
|
|
|
|
t ne х х dH + p f |
х х dV. (2.4.98) |
|
|
|
у |
|
Здесь использована формула ковариантного дифференцирования (см. т. 1, упр. 3 к § 2.4). А
Очевидно, что условия (2.4.95) необходимы и для разрешимости квазистатической задачи (2.2.9) в перемещениях.
Теорема 2.4.11. Среднее значение тензора напряжений, определенное как
<*■) |
1 |
a dV |
(2.4.99) |
|
V |
||||
|
|
|
v
для задачи в напряжениях (2.3.17) зависит только от массовых и поверх ностных сил и может быть вычислено до решения краевой задачи (2.3.17):
(<г) |
1 |
Р (f ® х + х ® f) dV + |
1 |
(tn (g) х + х (g) t n) dE. (2.4.100) |
|
W |
|
W |
X |
|
v |
|
|
▼Действительно, домножая уравнение равновесия (2.3.17а) тензорно на
|
§ 2.5. Нелинейно-упругие среды с малыми деформациями |
89 |
||||
вектор х, с учетом формулы (т. 1, (2.4.27)) получаем |
|
|||||
О = |
(V • сг) 0 х dV + p f (8) x |
= |
V • (cr 0 x) dy — (7 • V (8) x dV+ |
|||
у |
у |
|
У |
|
У |
|
|
+ °pi (g) х dV = |
t n (8) x dH — |
cr • E |
+ p f ^ x d V , |
(2.4.101) |
|
|
у |
|
У |
|
У |
|
или |
|
|
|
|
|
|
|
стdV |
р f (8) х йУ + t n (8) х d£. |
(2.4.102) |
|||
|
У |
У |
|
£ |
|
|
Симметризуя эту формулу, действительно приходим к (2.4.100). А
Упражнение к § 2.4
Упражнение 1. Доказать теорему 2.4.5.
§ 2.5. Нелинейно-упругие среды с малыми деформациями
2.5.1.Особенности принципа материальной симметрии
вслучае малых деформаций
До сих пор не рассматривался детальный вид определяющих соотноше ний — достаточно было иметь их общий вид в ^-форме (2.1.39) или в ("-форме (2.1.45), однако при переходе к решению конкретных задач необходим «кон кретный» вид зависимости функций ф, у и сг от тензора е.
Из общей теории определяющих соотношений, изложенной в т. 2, гл. 3, известно, что зависимости ф ~ е и сг ~ е должны записываться с помощью инвариантов тензора е. Такой способ записи определяющих соотношений обеспечивает их корректность — определяющие соотношения согласуются с
принципом материальной симметрии (при преобразовании отсчетной кон- |
|
о |
* |
фигурации JC в другую отсчетную |
конфигурацию JC определяющие соот |
ношения не изменяются, если это ^-преобразование принадлежит группе симметрии Gs данной твердой среды). Однако, если рассматриваются малые
деформации твердых сред, то приходим к противоречию. Действительно,
о *
iT-преобразования отсчетной конфигурации Н: 1C —►/С для основных групп
Gs твердых тел (группы изотропии, трансверсальной изотропии, ортотропии и другие группы кроме тождественной) предполагают наличие конечных
поворотов, а постулаты теории малых деформаций (см. п. 2.1.2) «разрешают»
о *
лишь малые повороты при преобразованиях: /С —►/С и /С —►/С. Таким образом,
90 |
Глава 2. Упругие среды с малыми деформациями |
если при переходе из /С в /С осуществляется конечный поворот, а из /С в /С —
о *
малый, то в итоге получаем конечный поворот JC—►JC —►/С.
Из-за этого противоречия формально мы «не имеем права» даже про верять определяющие соотношения (2.1.39) и (2.1.45) на согласованность с принципом материальной симметрии, хотя интуитивно ясно, что такая согла сованность должна быть. Выход из этого противоречия таков: при проверке определяющих соотношений на соблюдение принципа материальной симмет рии (и только при ней) вместо допущения (2.1.2) для градиента деформации
F используют соотношение |
|
F = F • Q T, |
(2.5.1) |
где Q — ортогональный тензор из некоторой группы Gs, а один из тензоров |
||
* |
|
|
F или F удовлетворяет соотношению |
|
|
F = Е + AF, |
| | A F | | < 1 , |
(2.5.2а) |
ИЛИ |
|
|
F = Е + AF, |
|| AF ||< 1. |
(2.5.26) |
Тогда определяющие соотношения (2.1.39), (2.1.40) или (2.1.45) должны конструироваться таким образом, чтобы выполнялся принцип материальной симметрии (в указанных выше ограничениях):
ф(е, |
в) = ф(е, в), |
(2.5.3а) |
^ Цв, |
0) = fj(e, в), |
(2.5.36) |
< ст= ;F(e, в), |
(2.5.3в) |
|
w*(e, |
в) = w*(e, в) |
(2.5.3г) |
(функция рассеивания й;*, как и ранее, — это единое обозначение, а не зна-
*>i< >i<
чение в конфигурации /С ), где е и сг — тензор малых*деформаций и тензор
напряжений, введенные для отсчетной конфигурации /С, отличающейся от /С iT-преобразованием с ортогональным тензором Q е Gs.
Если рассматриваем только идеальные среды, то вместо определяющих соотношений (2.1.39) имеем соотношения (2.1.40), для которых принцип материальной симметрии сводится к одному-единственному требованию (см. т. 2, теорема 3.8.1) для потенциала:
Ф(е, в) = Ф(е, в). |
(2.5.4) |
Аналогичное требование накладывается и на определяющие соотношения в С-форме (2.1.46), (2.1.47):
Z(a, в) = Z(a, в). |
(2.5.5) |