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

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

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

§ 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)

у

у

 

 

 

Если определить вариацию как разность двух статически допустимых полей 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)