УДК 531
ОБ ИСПОЛЬЗОВАНИИ ТЕНЗОРА ЛОГАРИФМИЧЕСКОЙ ДЕФОРМАЦИИ
О. С. Садаков1, А. О. Щербакова2
Отказ от гипотезы малости перемещений в механике деформируемого тела приводит к проблемам, связанным с геометрической нелинейностью, которые до сих пор не решены до конца, несмотря на свою актуальность для практических приложений. Одна из таких проблем относится к понятию логарифмической деформации. Скалярная логарифмическая деформации в ряде случаев полезна, но тензор логарифмической деформации обычно лишен смысла.
Ключевые слова: геометрическая нелинейность; конечные деформации; логарифмическая деформация; тензор Генки; инженерная деформация; простой сдвиг.
Введение
Тензор логарифмической деформации в настоящее время является одной из наиболее популярных мер деформации тела при решении геометрически нелинейных задач. Он представляет собой обобщение скалярной логарифмической деформации
1 1 ( l ^ '= [— dx = ln
x
l0
V l0 J
(І)
введенной Людвигом в 1909 г., перенесенной на область тензоров [1]. Идея Людвига состояла в том, чтобы сравнивать текущую длину растягиваемого стержня не с начальной длиной l0, как при использовании обычной инженерной деформации
£ = ^ = ^, (2)
11 l0 l0
а с длиной стержня в предыдущем состоянии. Логарифмическая деформация (1) однозначно выражается через инженерную деформацию
e = ln (1 + £) (3)
и наоборот, инженерная выражается через логарифмическую
£ = exp(e) -1. (4)
Таким образом, зная одну, легко получить другую. Для величины e Людвиг в свое время ввел термин effective specific strain - особая эффективная деформация, а несколько позже Надаи предложил новое название - натуральная деформация [2]. В настоящее время широкое распространение получили термины: истинная деформация и логарифмическая деформация.
В 1928 г. Генки [3] предложил распространить идею Людвига на случай трехмерной деформации, введя по аналогии с выражением (3) тензор логарифмической деформации, иногда называемый его именем - тензором Генки
e = ln(/ + £) . (5)
Здесь и ниже тензоры и векторы обозначены жирными символами - для того, чтобы отличать их от скаляров. В выражении (5) £ - это тензор инженерной деформации, связывающий векторы l0 волокон тела в начальном состоянии (до деформирования) с векторами l в деформированном состоянии
Ml = е-10, Ml = l -10, (6)
1 Садаков Олег Сергеевич - доктор технических наук, профессор, кафедра прикладной механики, динамики и прочности машин, Южно-Уральский государственный университет.
2 Щербакова Алла Олеговна - кандидат технических наук, доцент, кафедра прикладной механики, динамики и прочности машин, Южно-Уральский государственный университет.
E-mail: [email protected]
Садаков О. С., Об использовании тензора
Щербакова А.О. логарифмической деформации
I - единичный тензор. Сумму тензора є и единичного тензора обычно называют тензором удлинений (иначе, тензор растяжения или тензор коэффициентов длины)
и = I + є. (7)
В 1956 г. Дойл и Эриксон [4] ввели понятие обобщенной деформации
Е п =1 (ип -1). (8)
п
Здесь разные значения п соответствовали разным мерам деформации. В частности, при п = 0 у Дойла и Эриксона получился тензор логарифмической деформации.
Таким образом, как и в случае одноосной деформации, в общем случае имеются две альтернативные меры деформации: инженерная и логарифмическая. И, как показывает практика, логарифмическая деформация более популярна среди инженеров-механиков (при решении геометрически нелинейных задач). Это, по-видимому, объясняется тем, что использование логарифма позволяет заменять произведения суммами, упрощая решение задач. Например, чтобы найти относительное изменение объема элемента твердого тела, можно вместо традиционного выражения
^ = (1 + Єі)(1 + є2)(1 + Єз) -1, (9)
Уо
использовать более простое:
1п
V Уо )
= Є + е2 + е3 . (10)
Условие несжимаемости среды при этом принимает вид
е1 + е2 + е3 = ^ (11)
то есть тензор логарифмической деформации напоминает девиатор. Удобно также разделять логарифмическую деформацию е на упругую составляющую и неупругую:
Є = Єе + Єр . (12)
Эта сумма напоминает геометрически линейный подход.
При решении инженерных задач тензор логарифмической деформации связывают с некоторыми мерами напряжений, например, Кирхгофа [5] или Коши [6, 7]. Однако достоверность полученных результатов при этом вызывает серьезные сомнения [8]. Ниже будет показано, что использование тензора логарифмической деформации имеет смысл только в главных осях.
Вычисление тензора логарифмической деформации
Известно [9], что для получения какой-либо функции двухвалентного тензора Т достаточно знать его главные числа {Л} и главные направления, которые обычно задают единичными векторами { сі} - главными векторами. Особенность главных векторов состоит в том, что при ска-
лярном умножении на тензор Т они не поворачиваются:
Т• С- =Кс. (13)
Соответствующее характеристическое уравнение
Т -Л' I = 0 (14)
позволяет найти главные числа и векторы тензора Т. Для вычисления функции тензора необходимо записать его в главных координатах, применив эту функцию к главным значениям. У симметричного тензора все главные числа вещественны, а главные направления взаимно ортогональны. Следовательно, натуральный логарифм тензора и определяется в трехмерном пространстве суммой трех слагаемых
1п и = 1п(мі )сі сі, (15)
где і =1...3, иі - это главные числа тензора и, а сі - векторы декартового базиса, совпадающие с главными направлениями тензора и . Отметим, что, согласно выражениям (7) и (13), тензоры и и а имеют одинаковые главные направления, а их главные числа отличаются на единицу
иі =єі +1. (16)
В матричной форме выражение (16) имеет вид
[e ] =
ln u1 0 0
0 ln u2 0
0 0 ln u
(17)
4з_
При заниси матричных выражений следует номнить о том, в каком именно базисе это сделано. Нанример, выражение (17) содержит координаты тензора, занисанного в декартовом базисе, нанравления которого совнадают с главными нанравлениями этого тензора. В другом базисе координаты будут другими.
Деформации и жесткие повороты
При неоднородном деформированном состоянии тела обычно рассматривают деформирование элементарных объемов, материальные волокна которых бесконечно малы. Тогда деформированное состояние каждого из таких объемов можно нринять однородным. В этом случае волокна,
нрямые до деформирования, нри деформировании не искривляются, и их удобно обозначать век-
торами. При деформировании элементарного объема его главные волокна (волокна, онисываемые главными векторами тензора деформации) не новорачиваются, а только удлиняются или укорачиваются. Однако деформирование соседних элементарных объемов снособно вызывать жесткий новорот носледнего, нри этом его главные волокна новорачиваются. Онератором линейной связи между векторами элементарных волокон в начальном и деформированном состоянии в этом случае служит тензор дисторсии F
l = F • 1о, (18)
включающий в себя деформацию элементарного объема и его жесткий новорот
F = R • U = V • R . (19)
В этом выражении, которое обычно называют нолярным разложением тензора дисторсии, U и V -это соответственно нравый и левый тензоры удлинений (симметричные тензоры):
U = Єu +1, V = Єv +1, (20)
выраженные соответственно через нравый eU и левый eV тензор инженерной деформации, а R -тензор жесткого новорота на угол p нротив часовой стрелки вокруг оси, заданной единичным вектором n
= nn + (I -nn)cosp + Э • n sinp. (21)
Здесь Э = -5укeiejeк - это тензор Леви-Чевитты (5jk - символ Веблена). На нлоскости, нернен-дикулярной вектору n, в декартовом базисе матрица этого тензора имеет довольно нростой вид
cosp - sinp
[ R] =
(22)
sinp cosp
-1 T
Поскольку тензор жесткого новорота является ортогональным тензором (то есть R = R ), связь между левым и нравым тензором удлинений, следующая из нолярного разложения дисторсии (19), линейна
V = R • U • Rт. (2З)
Это выражение означает, что левый тензор удлинений нолучается новоротом нравого тензора на угол жесткого новорота.
Полярная декомнозиция тензора дисторсии основана на следующих соображениях. На рис. 1 ноказана схема дисторсии элементарного объема тела (для наглядности рассмотрена двумерная задача), в результате которой материальные волокна вытянулись (либо укоротились) и новерну-лись. Произвольный вектор l0, отвечающий некоторому множеству материальных волокон в начальном состоянии, изменил свою длину и нанравление. Сейчас это вектор lF = F • l0, отвечающий тому же множеству материальных волокон в деформированном состоянии. Этот вектор можно нредставить двумя снособами. В нервом снособе сначала рассматривают деформацию элементарного объема, нолучая нри этом вектор lu = U • l0 , а затем его жесткий новорот
If = R • lu = • іо . (24)
F
Об использовании тензора логарифмической деформации
В результате тензор днсторсни выражается через правый тензор растяжений Г = Я • и . Во втором способе сначала рассматривают жесткий поворот, получая вектор 1Я = Я • 10, а затем деформацию
ІР = V • 1К =У_^Я • /0. (25)
\| '
Рис. 1. Дисторсия элементарного объема
Тензор дисторсии при этом выражается через правый тензор растяжений Г = V • Я.
По аналогии с введением двух тензоров удлинений - правого и левого, вводят два тензора логарифмической деформации. Взятие логарифмической функции тензоров и и V, дает соответственно правый и левый тензоры логарифмической деформации:
еи = 1пи , eV = ^. (26)
Смысл координат
Обычно тензоры логарифмической деформации записывают в главных координатах, так как именно эта форма записи позволяет извлечь из них полезную информацию в виде логарифмических деформаций главных волокон. Например, координаты тензора
еи = 1п«хСхСх + 1пиУСУСУ , (27)
записанного в декартовом базисе { сх, су} (см. рис. 1), векторы которого направлены вдоль главных осей хи и уи этого тензора, представляют логарифмические деформации главных волокон. Этим, собственно, и ограничивается полезная с точки зрения механики информация, содержащаяся в этом тензоре. Получить какую-либо информацию о деформировании других волокон с помощью этого тензора невозможно.
Что касается тензора е^,, то он согласно выражениям (23) и (26) получается жестким поворотом правого тензора удлинений
eV = 1п(Я • и • Я т).
Следовательно, его главные координаты (х^ иуг на рис. 1) совпадают с главными координатами тензора еи (хи иуи на рис. 1). По сути, тензор eV содержит ту же информацию, что и тензор еи . Однако для получения этой информации необходимо не только направить базисные векторы вдоль главных осей, но и следить за их жестким поворотом, что усложняет координатную и матричную форму записи тензорных выражений [10].
В отличие от тензора логарифмической деформации (левого или правого), тензор инженерной деформации содержат существенно больше полезной информации о деформировании тела. Во-первых, согласно выражению (6), его проекция на вектор, задающий некоторое материальное волокно, дает вектор изменения этого волокна, определяющий его удлинение (укорочение) и поворот. Во-вторых, столбцы матрицы тензора е представляют изменения базисных векторов, причем базисные векторы могут и не совпадать с главными направлениями тензора е . В-третьих, диагональные элементы матрицы [ є ] представляют линейные деформации волокон в направле-
нии базисных векторов, а недиагональные - угловые деформации. И в-четвертых, тензор инженерной деформации, позволяя найти линейную деформацию любого волокна, дает с помощью выражения (3) возможность расчета скалярной логарифмической деформации
МІ
е = 1п (1 + є) , є = —, МІ = І —10 = І /о + є • /о | -1
І
0 .
(28)
Простой сдвиг
Одной из ключевых тестовых геометрически нелинейных задач является задача о простом сдвиге треугольной конструкции, схема которой показана на рис. 2, а. На рис. 2, а в начальном состоянии конструкция изображена прерывистой линией, а в деформированном - сплошной; оси х и у - произвольные оси декартового базиса; хг и уг - главные оси деформации. Задача о простом сдвиге имеет ряд существенных достоинств. Во-первых, деформация конструкции при простом сдвиге однородна; во-вторых, дисторсия содержит как деформацию, так и жесткий поворот; в-третьих, дисторсия в любой момент времени определяется только одним скалярным параметром - смещением д.
а) схема простого сдвига
б) первые главные деформации
в) вторые главные деформации
г) деформации горизонтальных волокон
д) деформации вертикальных волокон
Рис. 2. Простой сдвиг
е) сдвиговые деформации
На примере простого сдвига удобно демонстрировать различие между главными инженерными деформациями є1 и є2, а также соответствующими главными логарифмическими деформациями е1 = 1п(є1 +1) и е2 = 1п(є2 +1) в зависимости от смещения д. На рис. 2, а и 2, б главные инженерные деформации изображены сплошными линиями, а главные логарифмические - прерывистыми. При деформациях, превышающих 100 %, инженерные деформации могут отличаться от логарифмических в несколько раз. Однако при деформациях менее 10 % (что соответствует в данном примере смещению д = 0,2) различие не превышает 5 %, что легко объясняется использованием логарифмической функции.
Иногда тензоры логарифмической деформации записывают не в главных, а в произвольных координатах, применяя к выражениям (26) операцию перехода к новому базису
[е '] = [# ][е][н ]т. (29)
Об использовании тензора логарифмической деформации
II ЄУх ЄУху II 1 1 е хх е х у 1 —[ *г
■0 ху £ •0 ху £
[ к]. (30)
Здесь матрица [е'] содержит координаты тензора логарифмической деформации в новом базисе, матрица [е] - координаты тензора логарифмической деформации в старом базисе, а [Н] - матрица перехода к новому базису. Отметим, что в случае, когда тензоры логарифмической деформации записаны в произвольном базисе, смысл координат этих тензоров как мер деформации исчезает. На рис. 2 г, д и е показаны изменения координат тензоров инженерной и логарифмической деформации при простом сдвиге в осях х и у:
"1п(е1 +1) 0
0 1п(е2 +1)
Как и в случае главных деформаций, здесь различие между инженерными и логарифмическими деформациями возрастает с ростом д. Однако, в отличие от главных деформаций, даже при малых д это различие может быть значительным. Например, при д = 0,2 разница между инженерной и логарифмической деформацией в направлении оси х составляет порядка 30 % (еГх при этом равно 1,5 %), а в направлении оси у - порядка 100 % (еГу при этом равно 0,5 %). Различие между сдвиговыми деформациями начинает превышать 5 % при д > 1,18 (Х¥ху при этом равно 51 %).
В расчетах кинетики деформирования, когда главные направления деформаций поворачиваются в процессе деформирования, в разные моменты процесса деформирования тензоры логарифмической деформации (правый и левый) дают информацию о разных волокнах тела. На рис. 3 показано, что разным смещениям д отвечают разные пары главных направлений (оси ху и уу). Значит, волокна, информацию о которых содержат тензоры логарифмической деформации, также будут различаться в зависимости от д.
Рис. 3. Однородная дисторсия квадратного элемента при простом сдвиге
Логарифмическая деформация произвольного волокна легко вычисляется с помощью тензора инженерной деформации (28), тогда как попытка использовать вместо него тензор логарифмической деформации приводит к ошибочному результату. Этот факт проиллюстрирован на рис. 4 на примере волокна а , заданного в начальном состоянии в координатах х и у вектором
1
На рис. 4 изображен график изменения функции линейной деформации этого волокна в зависимости от параметра д (сплошная линия)
[а ] =
а
график скалярной логарифмической деформации (прерывистая линия)
Єу — 1п(єг + 1) ,
и величины
-0,5\
Рис. 4. Деформации волокна а при простом сдвиге
(1) _
1п
а + Єу ■ а
а
+1
вычисленной по аналогии с линейной деформацией, когда вместо тензора инженерной деформации используется тензор логарифмической деформации (пунктирная линия). Ошибочность использования тензора логарифмической деформации становится очевидна при д = 2, когда волокно а зеркально отражается относительно оси у, не изменив длины. Его деформация (как инже-
а + Єу • а - а
єу —
е
нерная, так и логарифмическая) при этом должна быть равна нулю. Однако величина eV(/) при q = 2 нулю не равна, она составляет -0,272.
Выводы
1. Скалярная логарифмическая деформация связана однозначной зависимостью с инженерной деформацией, то есть одна легко выражается через другую. Скалярную логарифмическую деформацию удобно использовать в тех случаях, когда в выражениях логарифмируется произведение инженерных деформаций. Тогда это произведение можно заменить суммой логарифмических деформаций и выражение при этом станет линейным.
2. Тензор логарифмических деформаций (правый и левый) содержит информацию только о логарифмических деформациях главных волокон. Если записать его в произвольном базисе, то координаты этого тензора с точки зрения деформирования объекта теряют смысл.
3. В расчетах кинетики деформирования, когда главные направления деформаций поворачиваются в процессе деформирования, тензор логарифмической деформации в разные моменты процесса деформирования дает информацию о разных волокнах тела.
4. Тензор инженерной деформации в отличие от тензора логарифмической деформации содержит следующую информацию о деформировании объекта.
• Его проекция на вектор, задающий некоторое материальное волокно, дает вектор изменения этого волокна, определяющий его удлинение (укорочение) и поворот.
• Столбцы его матрицы представляют изменения базисных векторов, причем базисные векторы могут и не совпадать с главными направлениями деформации.
• Диагональные элементы его матрицы представляют линейные деформации волокон в направлении базисных векторов, а недиагональные - угловые деформации.
• Тензор инженерной деформации, позволяя найти линейную деформацию любого волокна, дает возможность расчета скалярной логарифмической деформации.
Литература
1. Ludwik, P. Elemente der Technologischen Mechanik / P. Ludwik. - Berlin: Applied Mechanics, Verlag von J. Springer, 1909. - 57 p.
2. Nadai, A. Plastic Behavior of Metals in the Strain-Hardening Range. Part I / A. Nadai // J. Appl. Phys. - 1937. - Vol. 8. - P. 205-213.
3. Hencky, H. Uber die Form des Elastizitatsgesetzes bei ideal elastischen Stoffen / H. Hencky // Zeit. Tech. Phys. - 1928. - Vol. 9. - P. 215-220.
4. Doyle, T.C. Nonlinear elasticity / T.C. Doyle, J.L. Ericksen // Advances in applied mechanics IV. - New York: Academic Press, Inc. - 1959. - P. 53-115.
5. Аннин, Б.Д. Допустимые формы упругих законов деформирования в определенных соотношениях упруго-пластичности / Б.Д. Аннин, С.Н. Коробейников // Сибирский журнал индустриальной математики. - 1998. - Т. 1, № 1. - C. 21-34.
6. ANSYS - а general purpose finite element program. Rev. 5.0. / Houston (PA): Swanson analysis system Inc. - 1996. - 510 р.
7. Горлач, Б.А. Конечные деформации в задачах формования неупругих тел / Б.А. Горлач, Е.А. Ефимов // Математическое моделирование систем и процессов. - 1992. - № 1. - C. 67-81.
8. Садаков, О.С. Конечные деформации в механике деформируемого твердого тела / О.С. Са-даков // Вестник ЮУрГУ. Серия «Математика, физика, химия». - 2005. - Вып. 6. - № 6(46) -C. 114-121.
9. Ланкастер, П. Теория матриц / П. Ланкастер. - М.: Наука. - 1973. - 280 c.
10. Truesdell, C. The Classical Field Theories / C. Truesdell, R. Toupin; S. Flugge (ed) // Encyclopedia of Physics. - Berlin: Springer-Verlag, 1960. - Vol. III/1. - P. 226-793.
Поступила в редакцию 5 марта 2014 г.
Об использовании тензора логарифмической деформации
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics” ________________2014, vol. 6, no. 3, pp. 78-85
ON THE USE OF LOGARITHMIC STRAIN TENSOR
O.S. Sadakov', A.O. Scherbakova2
The rejection of the small displacements hypothesis in the mechanics of deformable solids leads to some problems connected with geometrical nonlinearity. At present, these problems are still not completely solved despite their relevance to practical applications. One of these problems is the concept of logarithmic strain. The scalar logarithmic strain is effective in some cases but the logarithmic strain tensor is usually useless.
Keywords: geometrical nonlinearity; finite strains; logarithmic strain; Hencky strain tensor; engineering strain; simple shear.
References
1. Ludwik P. Elemente der Technologischen Mechanik. Berlin: Applied Mechanics, Verlag von J. Springer, 1909. 57 p.
2. Nadai A. Plastic Behavior of Metals in the Strain-Hardening Range. Part I. J. Appl. Phys. 1937. Vol. 8. pp. 205-213.
3. Hencky H. Uber die Form des Elastizitatsgesetzes bei ideal elastischen Stoffen. Zeit. Tech. Phys. 1928. Vol. 9. pp. 215-220.
4. Doyle T.C., Ericksen J.L. Nonlinear elasticity. Advances in applied mechanics IV. New York: Academic Press, Inc. 1959. pp. 53-115.
5. Annin B.D., Korobeynikov S.N. Sibirskiy zhurnal industrial’noy matematiki. 1998. Vol. 1, no. 1. pp. 21-34. (in Russ.).
6. ANSYS - a general purpose finite element program. Rev. 5.0. Houston (PA), Swanson analysis system Inc, 1996. 510 p.
7. Gorlach B.A., Efimov E.A. Matematicheskoe modelirovanie sistem i protsessov. 1992. no. 1. pp. 67-81.
8. Sadakov O.S. Konechnye deformatsii v mekhanike deformiruemogo tverdogo tela (Finite deformations in mechanics of deformable solids). Vestnik YuUrGU. Seriya «Matematika, fizika, khimiya». 2005. Issue 6. no. 6(46). pp. 114-121. (in Russ.).
9. Lankaster P. Teoriya matrits (Theory of matrices). Moscow, Nauka Publ, 1973. 280 p.
10. Truesdell C., Toupin R., Flugge S. (ed) The Classical Field Theories. Encyclopedia of Physics. Berlin, Springer-Verlag, 1960. Vol. III/1. pp. 226-793.
Received 5 March 2014
1 Sadakov Oleg Sergeyevich is Dr. Sc. (Engineering), Professor, Department of Applied Mechanics, Dynamics and Strength of Machines, South Ural State University.
2 Scherbakova Alla Olegovna is Cand. Sc. (Engineering), Associate Professor, Department of Applied Mechanics, Dynamics and Strength of Machines, South Ural State University.
E-mail: [email protected]