УДК 517.17: 539.37
Об условиях выпуклости изотропной функции от тензора второго ранга
Владимир М. Садовский*
Институт вычислительного моделирования СО РАН, Академгородок, 50/44, Красноярск, 660036,
Россия
Получена 18.09.2010, окончательный вариант 25.10.2010, принята к печати 10.12.2010 Для скалярной функции, зависящей от инвариантов тензора второго 'ранга, получены условия выпуклости и сильной выпуклости относительно компонент этого тензора в произвольной декартовой системе координат. Показано, что если функция зависит только от четырех инвариантов: трех главных значений симметричной части тензора и модуля псевдовектора антисимметричной части, то полученные условия являются необходимыми и достаточными. Построена специальная система выпуклых инвариантов для конструирования потенциалов напряжений и деформаций в механике структурно неоднородных упругих сред, проявляющих моментные свойства.
Ключевые слова: изотропная тензорная функция, выпуклость, инварианты, нелинейная упругость, пластичность.
Выпуклые функции, определенные на пространствах тензоров напряжений, деформаций и скоростей деформации, применяются при конструировании потенциалов нелинейной теории упругости [1—3], диссипативных потенциалов и поверхностей текучести в теории вяз-копластических сред [4,5]. Для изотропных материалов, механические свойства которых не зависят от направления, такие функции — инварианты своих аргументов. В классических моделях механики деформируемых сред используются симметричные тензоры, обладающие тремя функционально независимыми инвариантами. Таковыми могут служить, например, главные значения тензора. Однако эти инварианты, вообще говоря, не являются выпуклыми функциями тензорных компонент, поэтому вопрос о выпуклости функции, зависящей от главных значений, не тривиальный.
В неклассических моделях, учитывающих микроструктуру материала, применяются несимметричные тензоры [6, 7]. В общем случае они имеют шесть функционально независимых инвариантов [8,9]. В настоящей работе строится специальная система выпуклых инвариантов для произвольного тензора второго ранга и приводятся новые или мало известные в специальной литературе условия выпуклости функции, заданной на пространстве инвариантов, относительно компонент тензора в произвольной системе координат.
Вводные понятия, случай симметричного тензора
Действительная функция f (u) от m независимых переменных u = (u 1, U2,..., um) называется выпуклой, если ее надграфик epi f = {(u, z) | f (u) ^ z} представляет собой выпуклое множество. По определению, это множество вместе с любыми своими точками (u, z) и (u, z) содержит точки отрезка (ux,zx) = A(u,z) + (1 — A)(u, z) (0 ^ A ^ 1). Выпуклая функция может принимать бесконечные значения, но тождественно равные и —ж функции в дальнейшем не рассматриваются.
*sadov@icm.krasn.ru © Siberian Federal University. All rights reserved
Для выпуклой функции множество уровня {и | /(и) ^ С} является выпуклым множеством при любом значении С.
Функция выпукла в том и только в том случае, если множество, на котором она конечна: ^т / = {и | /(и) < +<»}, выпукло и для любых элементов и и и этого множества выполняется неравенство
/(иЛ) < Л/(и) + (1 - Л)/(и) (0 < Л < 1). (1)
Первый критерий выпуклости состоит в следующем. Дифференцируемая функция /(и) выпукла в том и только в том случае, когда множество ^т / выпукло и для любых двух векторов и и и € ^т / справедливо неравенство
/(и) - /(и) > (и - и) • , (2)
в котором д/(и)/ди — вектор-градиент /(и), а точка между векторами означает скалярное произведение.
На его основе доказывается второй критерий выпуклости. Дважды непрерывно-дифференцируемая функция /(и) выпукла в том и только в том случае, когда матрица Гессе д2/(и)/ди2 неотрицательно определена в каждой точке выпуклого множества ^т /.
В качестве простого примера с помощью первого критерия можно показать, что если матрица А размерности гохи неотрицательно определена, то функция /(и) = VиАи выпукла на всем пространстве. Непосредственно из второго критерия следует, что в этом случае выпукла квадратичная форма /(и) = иАи. Для единичной матрицы А = I это означает выпуклость евклидовой нормы вектора |и| и скалярного квадрата |и|2.
Выпуклая функция /(и) называется сильно выпуклой, если вместо (1) выполняется более сильное неравенство
/(иЛ) < Л/(и) + (1 - Л)/(и) - Л(1 - Л) £ |и - и|2
с некоторой сколь угодно малой положительной постоянной £. Принимая во внимание очевидное тождество
|иЛ|2 = Л|и|2 + (1 - Л)|и|2 - Л(1 - Л)|и - и|2,
можно показать, что функция /(и) сильно выпукла в том и только в том случае, когда для некоторого £ > 0 функция /(и) - £ |и|2 выпукла. Таким образом, если /(и) дважды непрерывно-дифференцируемая функция, то для ее сильной выпуклости необходимо и достаточно, чтобы матрица д2/(и)/ди2 -2 £ I была неотрицательно определенной. Это условие обеспечивает равномерную положительную определенность матрицы Гессе для функции /(и) - важное свойство, во многих случаях гарантирующее термодинамическую корректность моделей механики.
Доказательства приведенных выше утверждений общеизвестны [10,11]. Далее в качестве вектора и будет рассматриваться вектор, составленный из компонент произвольного тензора второго ранга в трехмерном пространстве.
Условия выпуклости изотропной функции от компонент симметричного тензора получены в [12]. Пусть /(Л1, Л2, Л3) — выпуклая дважды непрерывно-дифференцируемая функция, зависящая от главных значений тензора Л = Л* (звездочка служит для обозначения сопряженного тензора), симметричная относительно аргументов:
/ (Ль Л2, Лз) = / (Л2, Ль Лз) = / (Л1, Лз, Л2).
Тогда порождаемая ею функция /(Л) выпукла относительно тензорных компонент Лij = Лji (г, 3 = 1, 2, 3) в произвольной декартовой системе координат.
Для доказательства рассматриваются три линейно независимые комбинации главных значений, предварительно перенумерованных в порядке убывания (Л1 ^ Л2 ^ Л3):
= ЛЬ М2 = Л1 +Л2, =Л1 + Л2 + Лз.
Оказывается, каждая из них является выпуклой функцией от Л^-. Выпуклость ^(Л) можно доказать с помощью известного представления для наибольшего собственного числа симметричной матрицы (принимается правило суммирования по повторяющимся индексам):
= тах Л^- щ V',
1Уг is j ,
|V| = 1
в соответствии с которым
М1(ЛЛ) = max ЛЛ Vj Vj = max (Л Л jj + (1 — Л) Л^) Vj Vj ^
|v|=1 |v| = 1
^ Л max Л j Vj Vj + (1 — Л) max Л^- Vj Vj = Л М1(Л) + (1 — Л) ^1(Л).
|v| = 1 |v| = 1
Выпуклость ^э(Л) очевидна, так как М3 = Лц + Л22 + Л33 — линейная функция. Наконец, комбинация М2 может быть представлена в виде суммы
М2 = М3 — Л3 = М3 — min Лjj Vj Vj = М3 + max {—Лгj Vj Vj},
|v| = 1 |v| = 1
и ее выпуклость следует из выпуклости каждого из слагаемых.
В силу выпуклости и симметрии f (Л1, Л2, Л3) выполняются неравенства
f > f > f дЛ1 > 5Л2 > 5Л3.
Действительно, полагая
f (Л) = f (Л Л1 + (1 — Л) Л2, Л Л2 + (1 — Л) Л1, Л3),
можно показать, что
f (Л) < Лf (Л1, Л2, Л3) + (1 — Л) f (Л2, Л1, Л3) = f (1), то есть что функция f (Л) не убывает в окрестности точки Л = 1. Таким образом,
> 0: (Л, — Л2,£ + (Л2 — Л0£ > О,
откуда df/dA.1 > ^/<9Л2, если Л1 > Л2. Если же Л1 = Л2, то производные df/dA.1 и df/дЛ2 оказываются равными в силу симметрии f (Л1, Л2, Л3) по первым двум аргументам. Аналогично доказывается неравенство df/dA.2 > df/^3.
Принимая во внимание, что Л1 = М1, Л2 = М2 —М1, Л3 = М3 —М2, получим два неравенства
f = f > 0 f = f > 0 (3)
дЛ1 ' дм2 5Л2 5Л3 ^ ' w
которые используются при проверке выпуклости f(Л). Матрица Гессе для этой функции равна
d2f = ^ B^q d2f дмр + ^ df д2Мр
dA-jj <9ЛЫ 5Лйг дмр dMq дЛjj дмр дЛjj 5Лйг'
Здесь слагаемое, отвечающее р = 3, опущено, так как мз линейно зависит от Лij и, следовательно, д2мз/дЛ^- дЛк = 0. Существование производных от главных значений Лк и, как следствие, от выпуклых инвариантов мр по компонентам тензора следует из известной теоремы матричной алгебры.
Неотрицательная определенность матриц, составленных из слагаемых первой суммы в правой части полученного равенства, доказывается непосредственно: если aij — произвольные коэффициенты, то
а д2/ дмр ^ = ^ д2/ ^
дЛи дмр дЛ^- ij 9 дмр дмд
где = дуи,р/дЛ^. Это выражение неотрицательно в силу второго критерия, записанного по отношению к функции / = /(м1, М2, Мз), выпуклость которой может быть установлена на основании неравенства (1) с учетом линейной зависимости Мк от Л1,Л2,Лз. Неотрицательная определенность матриц, отвечающих слагаемым второй суммы, следует из выпуклости М1(Л) и М2(Л), а также из неравенств (3). Таким образом, утверждение полностью доказано.
На самом деле справедливо обратное утверждение: если функция выпукла относительно компонент симметричного тензора в произвольной системе координат, то она, очевидно, выпукла в системе главных осей этого тензора относительно главных значений. Кроме того, она симметрична, поскольку в силу изотропности допускается произвольная перенумерация главных значений. Таким образом, рассматриваемые условия выпуклости являются необходимыми и достаточными.
Учитывая, что в силу инвариантности евклидовой нормы тензора
/(Л) - £ Л^Л^ = /(Л1, Л2, Лз) - £ (Л2 + Л2 +лз),
нетрудно доказать следующее утверждение.
Для того чтобы функция /(Л1, Л2, Лз) была сильно выпуклой относительно компонент симметричного тензора в произвольной декартовой системе координат, необходимо и достаточно, чтобы она была сильно выпуклой как функция трех переменных - главных значений тензора, и симметричной относительно этих переменных.
Случай несимметричного тензора
Справедливо более общее утверждение. Если изотропная дважды непрерывно-дифференцируемая скалярная функция /(Л) зависит только от четырех аргументов — трех главных значений симметричной части Ля = (Л + Л*)/2 тензора Л и модуля псевдовектора, ассоциированного с антисимметричной частью Ла = (Л - Л*)/2, то для ее выпуклости (сильной выпуклости) относительно компонент тензора Л, записанного в произвольной декартовой системе координат, необходимо и достаточно, чтобы она была симметричной относительно главных значений Ля:
/ (Л1, Л2, лз, |ла|) = / (Л2, Л1, Л|, |ла|) = / (Л1, лз, л^ , |ла|),
неубывающей по |Л°| (чтобы удовлетворяла условию д//д|Л°| ^ 4 £ |Л°| при некотором £ > 0) и выпуклой (сильно выпуклой) по совокупности аргументов.
Действительно, пусть / — выпуклая функция относительно компонент Л^-, зависящая только от указанных четырех аргументов - инвариантов тензора Л. В системе главных осей симметричной части этого тензора, которая может быть получена специальным поворотом
•Af 0 0 0 —ла
0 As л2 0, да = л3а 0
0 0 A3/ V—ла Ла
исходной координатной системы, тензоры Ля и л? принимают следующий вид:
л?
л = | 0 л? 0 | , л = | лз 0 -л?
Требуется показать, что функция /(л®,л?, А^, |ла|) является выпуклой, неубывающей по |л?| = ^(л?)2 + (л?)2 + (л?)2 и симметричной относительно л к.
Симметрия этой функции следует из ее изотропности - независимости от поворота системы координат, выпуклость следует из выпуклости / (л) на множестве тензоров частного вида
'л® -л? 0 \ л=|ла л? 0 (л? >0). 0 0 лз/
Неубывание функции по последнему аргументу доказывается методом от противного. Предположим, что в окрестности некоторой точки л®, л?, л|, |л?| = х эта функция строго убывающая. Тогда д/(л®, л?, л3,х)/д|л?| < 0. Рассмотрим следующую пару тензоров с одинаковыми системами инвариантов:
л —х 0 ( Л1 —х cos у
х As л2 0s, л = х cos у As л2
0 0 лз/ \—х sin у 0
As
(у — произвольный угол), на которых функция принимает одно и то же значение. Для этой пары
Í 0 х(1 - cos у) х sin уЛ
Л - Л = —х(1 - cos у) 0 0
у —X sin у 0 0
Учитывая, что в общем случае 4|Ла|2 = (Л32 — Л2з)2 + (Л13 — Л31)2 + (Л21 — Л12)2 д 1 / 0 Л12 — Л21 Л13 — АзЛ 1 /0 —х o^
ал =4^4 Л21— л12 0 Л23— М =2х х 0 0 дЛ 4|Л 1 \Л31 — Л13 Л32 — Л23 0 J 2х \0 0 0>
на основе второго критерия выпуклости получаем противоречие:
/I , N / д/ А/Т А ,8|Ла| д/ , 0 > i,C1(Aij — Aij) ^ = ^ J>i; — Aij) = — ^х(1 — cos > 0,
которое завершает доказательство необходимости.
Для доказательства достаточности после перенумерации главных значений симметричной части А| в порядке убывания рассмотрим выпуклые инварианты
М1 = Л1, М2 = Л1 + Л2, мз = Л1 + Л2 + Л3.
Кроме того, инвариант (Л) = |Ла | является выпуклой функцией относительно Aij- как частный случай функции вида \/uAu с неотрицательно определенной матрицей A.
Повторяя фрагмент приведенного выше доказательства для случая симметричного тензора, можно установить, что д//дм1 ^ 0 и д//дм2 ^ 0. Учитывая линейность уравнений, связывающих аргументы выпуклой функции /(А®, Л2, Л3, |Ла|) с инвариантами Mfc, можно с помощью неравенства (1) установить выпуклость функции /(мъ M2, Мз, M4):
/(мЛ,м4,мЛ, м4) = /((Л?)Л, (А2)л, (А3)л, |Ла|4) < Л/(Л1,Л2, Л3, |Ла|)+
0
+(1 - Л)/(Л1, Л2, Л^, |Ла|) = Л/(Дь Д2, Дэ, М4) + (1 - А)/(дь Д2, Дэ, М4).
По второму критерию выпуклости соответствующая матрица Гессе размерности 4 х 4 неотрицательно определена. Матрица Гессе по переменным Л^- вычисляется по формуле
д 2/ ^ дДд д2/ дДР , У^ д/ д2Др
+
дЛ^-дЛы дЛы ддрддд дЛ^- ^=3 дДр дЛ^-дЛы'
каждое из слагаемых в правой части которой определяет неотрицательно определенную матрицу. Следовательно, функция /(Л) выпукла.
Доказательство сильной выпуклости при выполнении условий, заключенных в скобки, можно получить, переходя от функции /(Л) к функции
/(Л) - £ЛуК, = /(Л1, Л2, ЛЭ, |Ла|) - £ ((ЛЦ2 + (Л2)2 + (ЛЭ)2 + 2|Ла|2).
Система выпуклых инвариантов
Изотропные функции, зависящие от инвариантов несимметричных тензоров, используются в теории моментных сред. В таких средах наряду с поступательными учитываются вращательные степени свободы частиц, поэтому тензор деформаций оказывается несимметричным. Однако доказанное утверждение практически непригодно в этой теории. Оно носит частный характер: в общем случае несимметричный тензор имеет шесть функционально независимых инвариантов Л| и Л?. Простым примером изотропной функции, которая не может быть выражена через четыре инварианта Л| и |Л°|, является определитель, описывающий изменение объема частицы:
= Л^ЛЭ + Л^Л?)2 + Л2(Л^ )2 + ЛЭ(Л? )2.
Фиксируя главные значения Л|, можно указать две системы инвариантов Л? и Л к с одинаковыми модулями |Л°| = |Лно с разными определителями.
Построим более содержательную систему выпуклых инвариантов, через которую определитель и другие характеристики деформированного состояния среды вычисляются однозначно. Пусть Д и V — ортогональный и симметричный тензоры, участвующие в полярном разложении Кэли:
Л = Д • V, Д • Д* = I, V = V*, V > 0.
Тензор V имеет три функционально независимых инварианта - главные значения V ^ V ^ Vэ, с помощью которых описывается деформация среды. Наибольшее из них неотрицательно в силу условия det V ^ 0. У тензора Д инвариантов тоже три, но особый интерес представляет только один из них — угол поворота частицы.
Покажем, что следующие комбинации инвариантов деформации:
Д1 = V!, Д2 = 7 V2 + VI, Дэ = ^ V2 + VI + (4)
являются выпуклыми. Действительно, так как Л* • Л = V2, то
Д1 = тах .
кИ1
Входящее сюда подкоренное выражение представляет собой неотрицательную квадратичную форму относительно компонент тензора Л, поэтому функция ^Л^Лу щи^ выпукла и, как следствие, выпукла функция Д1. В этом легко убедиться, проверяя неравенство (1).
На основании формулы дэ = \/Лг.Л. доказывается выпуклость дэ. Наконец Д2 можно записать в виде
Д2 ^ V Дэ - VI = Л.Лг. - тт ЛыЛк2щи. = тах ^Лг.Лг. - ЛьЛу^^. (5)
V у кИ1 кИ1
Из цепочки уравнений и неравенств
ЛкгЛк.^ < тахЛ^Лд.= V2 < V!2 + V2 + VI = Л.Лг.
кМ
вытекает, что подкоренное выражение в правой части (5) также представляет собой неотрицательную квадратичную форму относительно Л.. Отсюда следует выпуклость Д2.
В дополнение к системе (4) в качестве выпуклого инварианта, характеризующего поворот, примем Д4 = |Л°|. Этот инвариант тождественно равен нулю при отсутствии поворота частицы, когда Д = I, и отличен от нуля при Д = I.
Справедливо следующее утверждение. Если функция /(д1, Д2, дэ, Д4) выпукла и ее первые производные д//ддк по всем аргументам неотрицательны, то она выпукла относительно компонент тензора Л. Если, кроме того, выпуклой является функция /(Д1, Д2, Дэ, Д4) -е Дэ и выполняется неравенство д//ддэ ^ 2 е дэ при некотором достаточно малом значении е > 0, то функция / является сильно выпуклой относительно компонент тензора Л.
Заметим, что аналогичным способом можно выписать выпуклые инварианты несимметричного тензора, основываясь на левом разложении Кэли в произведение симметричного и ортогонального тензоров: Л = Ш • Д, и что вопрос о построении полной системы выпуклых инвариантов, состоящей из шести функционально независимых комбинаций, остается открытым.
В заключение приведем специальную систему условий, при выполнении которых функция, зависящая только от трех инвариантов У1, и Уэ несимметричного тензора Л, заведомо не является выпуклой. Принципиально то, что этим условиям автоматически удовлетворяет потенциал напряжений в теории гиперупругости [3] как функция от тензора дисторсии: /, = 0, если ^ = = V! = 1 (упругая энергия равна нулю в естественном, неде-формированном состоянии среды) и /^э) > 0, если (V - 1)2 + (V - 1)2 + (^э - 1)2 > 0 (упругая энергия положительна в деформированном состоянии).
От противного, допустим, что хотя бы одна из таких функций выпукла. Тогда для нее множество уровня
{Л | /(Л) < 0}
является выпуклым множеством. Так как /(Л) ^ 0, то это множество представимо в виде {Л | /(Л) = 0}. Из приведенной системы условий следует, что оно совпадает с группой ортогональных тензоров. Противоречие в том, что такая группа не образует выпуклое множество, поскольку очевидно, например, что выпуклая комбинация двух ортогональных тензоров Д = I и
/0 -1 0\ 1 1 /1/2 -1/2 0\ Д= 1 0 0 , равная - Д + - Д =1/2 1/2 0 , 0 0 1 2 2 0 0 1
не является ортогональным тензором.
Работа выполнена при финансовой поддержке РФФИ (код проекта 11-01-00053) и Междисциплинарного интеграционного проекта Сибирского отделения РАН № 40.
Список литературы
[1] Л.И.Седов, Механика сплошной среды, Т. 2, М., Наука, 1994.
[2] Ю.Н.Работнов, Механика деформируемого твердого тела, М., Наука, 1979.
[3] С.К.Годунов, Е.И.Роменский, Элементы механики сплошных сред и законы сохранения, Новосибирск, Научная книга, 1998.
[4] П.П.Мосолов, В.П.Мясников, Механика жесткопластических сред, М., Наука, 1981.
[5] Г.И.Быковцев, Д.Д.Ивлев, Теория пластичности, Владивосток, Дальнаука, 1998.
[6] В.А.Пальмов, Основные уравнения теории несимметричной упругости, Прикл. матем. и мех., 28(1964), вып. 3, 401-408.
[7] О.В.Садовская, В.М.Садовский, Математическое моделирование в задачах механики сыпучих сред, М., Физматлит, 2008.
[8] Э.Спенсер, Теория инвариантов, М., Мир, 1974.
[9] П.А.Жилин, Модифицированная теория симметрии тензоров и тензорных инвариантов, Известия ВУЗов, Северо-Кавказский регион, Естественные науки. Нелинейные проблемы механики сплошных сред, 2003, 176-195.
[10] Р.Т.Рокафеллар, Выпуклый анализ, М., Мир, 1973.
[11] В.Г.Карманов, Математическое программирование, М., Наука, 1980.
[12] W.H.Yang, A useful theorem for constructing convex yield functions, Trans. ASME. J. Appl. Mech., 47(1980), № 2, 301-305.
Conditions for Convexity of the Isotropic Function of the Second-rank Tensor
Vladimir M. Sadovskii
For a scalar function, depending on the invariants of the second-rank tensor, condition of convexity and strong con-vexity are obtained, with respect to the components of this tensor in an arbitrary Cartesian coordinate system. It is shown that if a function depends only on the four invariants: three principal values of the symmetric part of a tensor and modulus of pseudovector of the antisymmetric part, these conditions are necessary and sufficient. A special system of convex invariants is suggested to construct potentials for the stresses and strains in the mechanics of structurally inhomogeneous elastic media, exhibiting moment properties..
Keywords: isotropic tensor function, convexity, invariants, nonlinear elasticity, plasticity.