книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред
.pdf§ 2.4. Вариационные постановки задач |
71 |
Тогда для ди будет выполнено интегральное соотношение (2.4.4):
сг(е(и)) • • 6s(u) dV = |
pi -d u d V + t ne-5udS . |
(2.4.8) |
V |
У |
|
Уравнение (2.4.8) называют вариационным уравнением для квазистатической задачи МДТТ. Запишем его в несколько ином виде.
Обозначим работу внешних поверхностных и массовых сил на переме щениях и:
4 е |
4 е |
_и 4 е |
4 е |
— t п е |
’ U d Y l, А т |
— |
pi • и dV, |
(2.4.9) |
л. |
— Лу* |
л т, |
л ^ |
у
а потенциальную энергию твердого тела:
П = |
(2.4.10) |
у
где Ф = рф(е,в) — потенциал (2.1.40д), и запишем функционал
L(u) = П - А е, |
(2.4.11) |
называемый лагранжианом. Тогда если рассматривать кинематически допу стимое поле и(х) в классе дважды непрерывно-дифференцируемых функций, то имеет место следующая теорема.
Теорема 2.4.1 (вариационный принцип Лагранжа). Пусть тензорный оператор а(е) для твердого тела является потенциальным, т. е. удовле творяет условию (2.1.406):
<у = <х(е) = дЯ//де, |
(2.4.12) |
тогда среди всех кинематически допустимых полей и(х) действительное поле отличается тем, что для него и только для него лагранжиан L имеет стационарное значение:
5L(и) = 0. |
(2.4.13) |
Замечание 2.4.1. В условия теоремы 2.4.1 неявно входит еще и определе ние вариации функционала 6L — достаточно простое, если оператор <х(е) является просто тензорной функцией <т(е), и совсем нетривиальное, если это оператор интегрального типа (см. т. 1, §4.15). Далее будем полагать, что определяющие соотношения являются тензорными функциями.
▼ Пусть выполнено уравнение (2.4.13) для поля и, тогда, подставляя выражения (2.4.9)—(2.4.11) в (2.4.13), имеем
5Ь = 5Ф dV - pi •SudV - tnp • Su dH = 0. |
(2.4.14) |
уу
72 |
Глава 2. Упругие среды с малыми деформациями |
Учитывая замечание 2.4.1 и принимая во внимание потенциальность опе ратора (2.4.12), получаем
5Ф = — • • 5е = сг(в(и)) • • 5е(и). |
(2.4.15) |
Подставляя (2.4.15) в (2.4.14), приходим к вариационному уравнению (2.4.8). Тогда, сделав все выкладки в обратном порядке, снова получим интегральное соотношение (2.4.2):
(п • сг — t ne) • 5и (П2 — 5и • (V • ст + pi) dV = О, |
(2.4.16) |
V |
|
из которого, в силу произвольности функций 5и, следует уравнение равнове сия (2.2.9а).
Остальные соотношения (2.2.9б)-(2.2.9д), очевидно, при этом тоже вы полняются, т. е. рассматриваемое поле и — это действительное поле переме щений.
Докажем теорему в обратную сторону. Пусть и — действительное поле перемещений, т. е. удовлетворяет всей системе (2.2.9). Тогда, повторив все выкладки (2.4.1)—(2.4.8), снова приходим к вариационному уравнению (2.4.8) для и. Принимая во внимание потенциальность оператора (2.4.12) и соот ношение (2.4.15), получаем, что из формулы (2.4.8) следует соотношение (2.4.14), которое и означает, что лагранжиан (2.4.11) имеет стационарное значение для поля и. А
Вариационная постановка квазпстатпческой задачи (2.2.9) заключается в нахождении кинематически допустимого поля перемещений и(х), удовле творяющего вариационному уравнению (2.4.8).
2.4.3. Вариационный принцип Хеллингера — Рейсснера
Кроме вариационного принципа Лагранжа, существуют и другие вариа ционные принципы. Рассмотрим наиболее широко применяемый на практике такой принцип.
Будем исходить из вариационного уравнения (2.4.8). Обозначим несколько иначе дифференциальный оператор (2.4.5):
е = L(u), |
(2.4.17) |
L(u) = def u |
(2.4.18) |
и добавим нулевое слагаемое в (2.4.8):
|
§ 2.4. Вариационные постановки задач |
73 |
||
а(е) • • |
SL(u) dV + Scr(e) |
(L(u) - e) d V - |
|
|
v |
v |
|
|
|
|
|
• d u d V - |
t ne • 5u |
= 0. (2.4.19) |
|
|
у |
|
|
Примем, что перемещения u и деформации е являются независимыми искомыми функциями, т. е. полагаем, что соотношение Коши (2.4.17) уже может не иметь места. Тогда вариационное уравнение (2.4.19) превращается
вновое самостоятельное уравнение, а вариационный принцип формулируется
ввиде следующей теоремы (если и рассматривают в классе дважды, а в —
вклассе один раз непрерывно-дифференцируемых функций).
Теорема 2.4.2 (принцип Хеллингера — Рейсснера). Пусть определяющие соотношения сг(е) являются Н-потенциальными, т. е. существует такая скалярная функция Ф(е), для которой
дсг(е) _ дФ(е) де де
тогда среди всех пар полей (и,в), где и — кинематически допустимое поле, е — некоторое произвольное симметричное тензорное поле, действи тельные поля и и е (т. е. удовлетворяющие всем соотношениям (2.2.9)) отличаются тем, что для них и только для них имеет стационарное значение следующий функционал:
8J{и, е) = О, |
(2.4.20) |
где
J (и, е) = (с(е) • • L(u) - Ф(е)) dV |
pi • и dV |
t ne • и |
(2.4.21) |
у |
у |
|
|
▼Пусть выполнено уравнение (2.4.20) для полей и и в, тогда, вычисляя вариацию от функции Ф
8Ф(е) |
= ^ • • 8е = е • |
де |
. . Se = е • • Scr(e) |
(2.4.22) |
|
де |
|
|
|
и подставляя ее в (2.4.21) и (2.4.20), получаем уравнение |
|
|||
8J = (сг(е) • • |
5L(u) + Scr(e) • • |
L(u) —е • • dcr(e)) dV SAe |
0, (2.4.23) |
|
v |
|
|
|
|
которое в точности совпадает с (2.4.19).
Примем во внимание, что 5L(u) = L(5u), и преобразуем первый интеграл в (2.4.19) в соответствии с формулами (2.4.3):
74 |
Глава 2. Упругие среды с малыми деформациями |
|||
a --S L (u ) dV = |
СV • (Su • <т) - |
5u • V • <т) dV = |
|
|
у |
У |
|
|
|
|
|
|
п • сг • SudU |
5u- V -a dV. (2.4.24) |
|
|
|
у |
|
Подставляя (2.4.24) в (2.4.19), приведем это уравнение к виду |
||||
(V |
• сг + pi) • Su + |
(n • сг —t ne) • Su dH + (L(u) - |
e) - ~ - - S e d V = 0. |
|
|
|
|
|
OS |
V |
|
TI<J |
V |
(2.4.25) |
|
|
|
|
В силу произвольности полей Su и Ss, отсюда следует, что выполняются уравнение равновесия (2.2.9а), граничные условия (2.2.9г) и соотношения Коши (2.1.75): е = L(u), остальные соотношения в (2.2.9) выполнены по условию теоремы. Следовательно, поля (и, в) — это действительные поля перемещений и деформаций.
Докажем теорему в обратную сторону. Если и и е — действительные поля, удовлетворяющие всей системе (2.2.9), то можно образовать интегральное соотношение (2.4.25), а затем, повторяя все выкладки в обратном порядке, вновь получить уравнение (2.4.20). А
Отметим, что поскольку в вариационном уравнении Хеллингера — Рейсснера (2.4.19) вариации Su и 5е являются независимыми, то вместо соотно шения (2.4.19) можно использовать систему вариационных уравнений:
'J V ( e ) - SL(u) dV = SAe(u),
(2.4.26)
J ( L ( u ) - £) .A - W ,y
Тогда вариационная постановка Хеллингера — Рейсснера для квазистатической задачи (2.2.9) заключается в нахождении системы из кинематически допустимых полей перемещений и(х) и поля тензора е(х), удовлетворяющих системе вариационных уравнений (2.4.26).
2.4.4. Вариационная постановка динамических задач МДТТ
Для динамической задачи (2.2.5) также можно сформулировать вариаци онную постановку.
Введем три специальных класса полей: 1) класс переменных векторных полей w(x, t) (х G V U X, t > 0), непрерывно-дифференцируемых один раз по х и два раза по t и удовлетворяющих граничному условию (2.4.1); 2) и 3) классы векторных полей у(х) и z(x), определенных в области V.
Тогда, умножая уравнение движения (2.2.5) на векторное поле w и ин тегрируя его по V, а также умножая начальные условия (2.2.5е) и (2.2.5д)
§ 2.4. Вариационные постановки задач |
75 |
соответственно на поля у и z и интегрируя по V, получаем следующее интегральное тождество:
w • (V |
д2' |
У • (u - и0) dV + |
z • (v —Vo) |
= 0. |
• а + р (f - —у ) ) d V + |
||||
у |
дГ |
у |
у |
|
|
|
(2.4.27)
Здесь введены обозначения для перемещения и скорости при t = 0:
u(x) = u(x, 0), v(x) = ^ ( x , 0). (2.4.28)
Преобразовав первый в (2.4.27) интеграл по формулам (2.4.3)-(2.4.5),
имеем |
|
|
|
|
|
|
|
о |
д2и ч |
7ТГ |
tne ‘ W (П2 |
(T (S (U )) • • e(w) dV + |
|
||
P w |
• (f - -7^2) |
d V + |
|
||||
v |
|
|
|
|
V |
|
|
|
|
|
+ |
у • (u - u0) dV + |
z . (v —vo) dV = 0. |
(2.4.29) |
|
|
|
|
у |
|
|
у |
|
Выберем поля w, у |
и z следующим образом: |
|
|||||
|
|
w = 5и, |
у = т |
5и, |
z = п 5v. |
(2.4.30) |
Здесь ди — разность двух кинематически допустимых переменных векторных полей при t > 0; ди — разность двух допустимых полей начальных перемеще ний при t = 0; 5v — разность двух допустимых полей начальных скоростей при t = 0.
Допустимыми полями начальных перемещений и скоростей называем про извольные стационарные векторные поля в V, не зависящие от t, а действи тельными полями начальных перемещений и скоростей называем такие поля и и v, которые удовлетворяют условиям (2.2.5д) и (2.2.5е), т. е. совпадают с
и0 и v0.
Ненулевые константы т и п выбираем из условия, чтобы все слагаемые в нижеследующем уравнении (2.4.31) имели одинаковую физическую размер ность.
Тогда, подставляя (2.4.30) в (2.4.29), получаем следующее вариационное уравнение для динамической задачи (2.2.5):
<х(е(и)) • • 5е(u) dV + |
г ~ д2и |
pi • 6и dV + t ne • 5и |
р—FT • 6и dV = |
||
у |
dt2 |
У |
У |
||
+ т |
(и —UQ) • Su dV + n (v —vo) • 5v dV = 0. (2.4.31) |
|
у |
|
у |
76 |
Глава 2. Упругие среды с малыми деформациями |
Обратим внимание на тот факт, что если во втором слагаемом заменить вариацию ди на оператор дифференцирования du/dt, то это слагаемое будет представлять собой полную производную по t от кинетической энергии К (т. 2, (2.4.1)):
К = 2 °p\dn/dt\2dv, |
о d |
ди |
dV = |
>du д2и |
-dT = -xJ dt |
dt |
d v ■ (2-4 -з2> |
||
|
р— |
|
|
|
V |
V |
|
|
V |
поскольку при малых деформациях область V неподвижна.
Введем «почти кинетическую энергию» К ' — функционал, вариация ко торого формально совпадает с производной от кинетической энергии (т. 2,
(2.4.1)), если оператор d/dt заменить на вариацию 5, т. е. |
|
|
6К' = |
5и dV, |
(2.4.33) |
|
у |
|
при этом сам функционал К ' явного вида не имеет.
Иногда вектор f' = d2u /d t2 условно рассматривают как вектор массовых сил инерции, который так же, как и f, не варьируют. Тогда для К ' можно ввести явное выражение, однако вариация такого функционала SK' будет совпадать с (2.4.33) лишь условно, т. е. с условием, что варьируется поле и, но не варьируется d2u /d t2.
Введем также начальное значение потенциальной кинетической энергии:
П0(й) = у |й - и 0|2 сП/, |
K Q(V ) = | |v - v0|2 dV, |
(2.4.34) |
v |
V |
|
тогда можно сформулировать следующий вариационный принцип, полагая, что и — дважды непрерывно-дифференцируемая функция по х и t.
Теорема 2.4.3 (вариационный принцип Даламбера). Пусть тензорный оператор <т(е) является потенциальным, т. е. удовлетворяет условию (2.4.12), тогда среди всех допустимых систем полей (u(x, t), u(x), v(x)) действительные системы полей отличаются тем, что для них и только для них следующий функционал имеет стационарное значение:
SD(u, u, v) = О, |
(2.4.35) |
где |
|
D = К ' + П - (KQ + П0) - Ае. |
(2.4.36) |
Т Для доказательства достаточно вычислить вариацию |
|
SB = 8К ' + 8118KQ - <5П0 - 8Ае |
(2.4.37) |
§ 2.4. Вариационные постановки задач |
77 |
и подставить выражения (2.4.9), (2.4.10), (2.4.33), (2.4.34) в (2.4.37), а затем (2.4.37) в (2.4.35), в результате действительно придем к вариационному уравнению (2.4.31).
Если выполнено соотношение (2.4.35), то имеет место и (2.4.31), которое, в силу независимости вариаций 5и, 5и и 5v, эквивалентно системе трех
вариационных уравнений: |
|
|
' |
гл2 |
(2.4.38) |
<fe(u) dV + |
°р2Л . 5udV = SAe, |
|
J |
дt2 |
|
V |
|
|
(u —uo) • <5u dV = 0, |
(2.4.39) |
|
r |
|
|
(v —VQ) ‘ 5v dV = 0. |
(2.4.40) |
|
у |
|
|
Преобразовав уравнение (2.4.38) с помощью формул (2.4.2)-(2.4.5) в об ратном порядке, получим
|
dV - |
(n • сг — t ne) • 6u dH = 0. |
(2.4.41) |
Из |
уравнений (2.4.39)-(2.4.41), в силу произвольности функций 5и, 5и |
||
и 5v, |
следует, что выполняются все уравнения системы (2.2.5), т. е. поля и, |
||
о |
о |
|
|
=VQ являются действительными.
Вобратную сторону теорема будет доказана, если провести все рассужде ния от (2.2.5) и (2.4.39)-(2.4.41) до (2.4.37), (2.4.35) в обратном порядке. А
Вариационная постановка динамической задачи МДТТ (2.2.5) заключа ется в нахождении системы функций (и, и, v), удовлетворяющих системе вариационных уравнений (2.4.38)-(2.4.40).
При численном решении динамической задачи МДТТ, например методом конечных элементов (МКЭ), система (2.4.38)-(2.4.40) сводится к системе
обыкновенных дифференциальных уравнений второго порядка с начальными условиями.и = UQ, v
2.4.5. Вариационная постановка задачи теплопроводности твердого тела
Сформулируем вариационную постановку задачи теплопроводности (2.2.8). Введем специальные классы:
1)скалярных переменных полей $(х,£), определенных в T U S Vt > 0, непрерывно-дифференцируемых один раз по х и по t и удовлетворяющих
78 |
Глава 2. Упругие среды с малыми деформациями |
|
|
|
нулевому граничному условию на части £# поверхности £ |
области V: |
|
|
i? Е. |
0; |
(2.4.42) |
2) стационарных скалярных полей /3(х), определенных в области V . Тогда, домножая уравнение теплопроводности (2.2.8) на функцию д и
интегрируя по V, а затем вычитая из него уравнение (2.2.8г), умноженное на
/3 и проинтегрированное по V, получаем |
|
|
|||
' ° дв q , |
i?V • (Л • V0) dV - |
pqmP dV - |
(в-e0)pdv 0, |
(2.4.43) |
|
pc*;— i? dV - |
|||||
v |
v |
|
V |
V |
|
где в = 0(x, 0). |
|
|
|
|
|
Преобразуем второе слагаемое: |
|
|
|
||
i?V • (Л • V0) dV |
V • (i?A • V0) dH - Vi? • Л • V0 dH = |
|
|||
v |
|
у |
у |
|
|
|
|
|
gei? dH |
Vi? • Л • V0 dK |
(2.4.44) |
|
|
|
|
у |
|
Здесь использовано граничное условие (2.2.86) и учтено условие (2.4.42) для функции 1?.
Назовем нестационарное поле 0(х, £) (х е У, t > 0) допустимым, если оно является один раз непрерывно-дифференцируемым по х в У U S и по t и удовлетворяет граничному условию (2.2.8в):
0 |Ев = 0е. |
(2.4.45) |
О |
|
Допустимым стационарным полем 0(х) (х е Н) назовем произвольное
скалярное поле в Н.
Выбирая i? и (3 в виде |
|
|
= 50, |
/3 = /?о50 |
(2.4.46) |
и подставляя в (2.4.43), с учетом (2.4.45) получаем вариационное уравнение для нестационарной задачи теплопроводности:
° дв 7Тг |
V0 • Л • V50 dH + |
ge50 dV — |
рсу—50 (ЗП + |
||
v |
v |
|
|
pqmd0 dV - A) (i0 - в0)5в dV = 0. (2.4.47) |
|
|
у |
у |
Здесь /?о — ненулевая размерная константа.
§ 2.4. Вариационные постановки задач |
79 |
Если обозначить |
|
|
QE = - qe6 dV, |
Qm = p qm6 dV, |
|
|
V |
|
m = ( i / 2) V0 • Л • V0 dV, |
U0{6) = (/30/2) (в - во)2 dV, |
(2.4.48) |
у |
у |
|
а также определить функционал U{6), вариация которого имеет вид |
|
|
8U(6) |
|
(2.4.49) |
у |
|
то можно сформулировать следующий вариационный принцип, рассматривая допустимые поля в в классе дважды непрерывно-дифференцируемых функций
по х.
о
Теорема 2.4.4. Среди всех допустимых полей (0(х, t), в(х)) действитель ные поля, т. е. удовлетворяющие всей системе (2.2.8), отличаются тем, что для них и только для них следующий функционал имеет стационар
ное значение: |
_ |
0 |
|
|
|
5Ь(в, |
в) = |
0, |
(2.4.50) |
где |
в) = U(в) - и0(°в) + Л(0) - |
Q sW - 0т(в). |
(2.4.51) |
|
Ь(в, |
▼Для доказательства достаточно подставить формулы (2.4.48), (2.4.49) и (2.4.51) в (2.4.50) и убедиться в том, что уравнение (2.4.50) совпадает с (2.4.47). Дальнейший ход доказательства аналогичен доказательству теоре
мы 2.4.1. А
о
В силу произвольности вариаций 86 и 86, из (2.4.47) получаем вариацион
ную постановку нестационарной задачи теплопроводности (2.2.8) относитель
но
но системы функций (6(x,t),6(x)): |
|
|
|
||
' J pCv^SedV = |
- J V 0 • Л • VSOdV - J qe86dv + J °pqmS9dV, |
|
|||
< V |
dt |
V |
Eg |
У |
(2.4.52) |
j(6 - |
60)86dV = 0. |
|
|
||
|
|
|
|||
<y |
|
|
|
|
|
2.4.6. Вариационная постановка квазистатической задачи |
|||||
|
|
для несжимаемой среды |
|
||
Сформулируем |
вариационную |
постановку |
квазистатической задачи |
(2.2.13) для несжимаемой среды.
Введем кроме векторного поля w(x), удовлетворяющего (2.4.1), скалярное поле у(х), определенное в V U £. Домножая уравнение равновесия (2.2.13а)
80 |
Глава 2. Упругие среды с малыми деформациями |
на w, а условие несжимаемости (2.2.136) на у, интегрируя получившееся выражение по У и складывая их друг с другом, получаем интегральное соотношение
w(Vp + V • сг(е)) dV + |
°pi • w dV + |
Vh(e) dV = 0. |
(2.4.53) |
у |
у |
У |
|
Преобразовав первое слагаемое в (2.4.53) по формулам (2.4.1)—(2.4.5), находим
(<х(е(и)) —рЕ) • • e(w) dV |
yh(e(u))dV = |
|
v |
V |
|
|
°pi • w dV + |
t ne • w dE. (2.4.54) |
Выбирая в качестве w и у: |
у |
|
|
|
|
w = Su и у = 6p, |
(2.4.55) |
где dp — вариация допустимых полей гидростатических напряжений, т. е. определенных в П и И, но не обязательно удовлетворяющих системе (2.2.13), получаем вариационное уравнение
((т(в(и)) —_рЕ) • • 6s(u) dV ■ |
Ii {s{u))6pdV = Ae(u). |
(2.4.56) |
У |
У |
|
Если ввести гидростатическую потенциальную энергию |
|
|
Пр — |
Ph(e) dV |
(2.4.57) |
v
и лагранжиан
L \u ,p ) = U + Up - A e, |
(2.4.58) |
где П и Ае определяют по (2.4.9) и (2.4.10), то можно сформулировать вариационный принцип.
Теорема 2.4.5. Пусть тензорный оператор а(е) потенциален, т. е. удо влетворяет (2.2.12), тогда среди всех систем (и,_р) кинематически допу стимых полей перемещений и допустимых полей гидростатических на пряжений действительные поля, т. е. удовлетворяющие системе (2.2.13), отличаются тем, что для них и только для них лагранжиан V имеет стационарное значение:
5L'(и , р) = 0. |
(2.4.59) |
▼Доказательство теоремы оставим в качестве упр. 1 к § 2.4. |
А |
В силу независимости др и ди, уравнение (2.4.56) можно представить в виде системы