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

книги / Механика сплошной среды. Т.4 Основы механики твёрдых сред

.pdf
Скачиваний:
19
Добавлен:
19.11.2023
Размер:
41.7 Mб
Скачать

§ 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Ь = 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Ф(е)

= ^ • • = е

де

. . 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):

• • SedV = 0.

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 и являются независимыми, то вместо соотно­ шения (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), в силу произвольности функций ,

и 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) можно представить в виде системы