CC BY
16
ЧЕБЫШЕВСКИЙ СБОРНИК Том 18 Выпуск 3
УДК 539.3 Б01 10.22405/2226-8383-2017-18-3-254-278
ТОЧНЫЕ РЕШЕНИЯ ЗАДАЧ ТЕОРИИ
МНОГОКРАТНОГО НАЛОЖЕНИЯ БОЛЬШИХ ДЕФОРМАЦИЙ ДЛЯ ТЕЛ, ОБРАЗОВАННЫХ ПОСЛЕДОВАТЕЛЬНЫМ СОЕДИНЕНИЕМ ДЕФОРМИРОВАННЫХ
__и
ЧАСТЕЙ1
К. М. Зингерман2 (г. Тверь), Л. М. Зубов3 (г. Ростов-на-Дону)
Аннотация
В статье приведены и развиты разработанные совместно с профессором МГУ им. М.В. Ломоносова В.Ан. Левиным подходы к точному аналитическому решению задач о больших деформациях составных изделий (тел) из несжимаемых изотропных нелинейно-упругих материалов, части которых предварительно деформированы. Решение этих задач представляет интерес при анализе напряжений в элементах конструкций, изготавливаемых из предварительно нагруженных частей. Результаты могут быть использованы для тестирования промышленного программного обеспечения, предназначенного для численного моделирования аддитивных технологий.
Постановка задач осуществляется на основе теории наложения больших деформаций и в рамках этой теории может быть сформулирована следующим образом. Части изделия, первоначально не связанные между собой, подвергаются начальному деформированию и переходят в промежуточное состояние. Затем эти части соединяются между собой. Соединение происходит по некоторым поверхностям, общим для каждой пары соединяемых частей. Далее тело, составленное из нескольких частей, деформируется как единое целое под действием приложенной к нему дополнительной нагрузки и переходит в конечное состояние. Предполагается, что на поверхностях, по которым соединены части тела, выполняются условия идеального контакта, т.е. векторы перемещений в соединяемых частях изделия на этих поверхностях совпадают.
1 Работа выполнена при частичной поддержке Министерства образования и науки России в рамках базовой части Государственного задания (проект № 9.7446.2017/8.9).
2 Зингерман Константин Моисеевич, заведующий кафедрой вычислительной математики Тверского государственного университет, Zingerman.KM@tversu.ru
3 Зубов Леонид Михайлович, профессор кафедры теории упругости Института математики, механики и компьютерных наук ЮФУ, lmzubov@sfedu.ru
Точные решения для изотропных несжимаемых материалов найдены с использованием известных универсальных решений и могут быть рассмотрены как обобщение этих решений на случай наложения больших деформаций.
В статье детально рассмотрены следующие задачи:
— задача о напряженно-деформированном состоянии в двух полых круговых упругих цилиндрах (трубах), один из которых был предварительно деформирован и вставлен в другой цилиндр (задача Ламе-Гадолина);
— задача о кручении составного цилиндра;
— задача о больших деформациях изгиба составного бруса, состоящего из нескольких предварительно деформированных частей (слоев). Приведены математические постановки этих задач, методы и некоторые результаты их решения. Исследовано влияние предварительных деформаций на напряженно-деформированное состояние, анализируются нелинейные эффекты.
Ключевые слова: нелинейная теория упругости, наложение больших деформаций, предварительно деформированные тела, точные аналитические решения, аддитивные технологии.
Библиография: 25 названий.
EXACT SOLUTIONS OF PROBLEMS OF THE
THEORY OF REPEATED SUPERPOSITION OF LARGE STRAINS FOR BODIES CREATED BY SUCCESSIVE JUNCTION OF STRAINED
PARTS
K. M. Zingerman (Tver), L. M. Zubov (Rostov-na-Donu)
Abstract
Large strains of composite solids made of incompressible isotropic nonlinear-elastic materials are analyzed for the case in which the parts of these solids are preliminarily strained. The approaches to exact analytical solutions of these problems are given and developed in cooperation with V.An. Levin. He is a professor at the Lomonosov Moscow University. The solution of these problems is useful for stress analysis in members containing preliminarily stressed parts. The results can be used for the verification of industrial software for numerical modeling of additive technologies.
The problems are formulated using the theory of repeated superposition of large strains. Within the framework of this theory these problems can be formulated as follows. Parts of a member, which are initially separated from one another, are subjected to initial strain and passes to the intermediate state. Then these parts are joined with one another. The joint is performed by some surfaces that are common for each pair of connected parts. Then
the body, which is composed of some parts, is strained as a whole due to additional loading. The body passes to the final state. It is assumed that the ideal contact conditions are satisfied over the joint surfaces. In other words, the displacement vector in the joined parts is continuous over these surfaces.
The exact solutions for isotropic incompressible materials are obtained using known universal solutions and can be considered as generalizations of these solutions for superimposed large strains.
The following problems are considered in detail:
— the problem of stress and strain state in two hollow circular elastic cylinders (tubes) one of which is preliminarily strained and inserted into another cylinder (the Lame-Gadolin problem);
— the problem of torsion of a composite cylinder;
— the problem of large bending strains of a composite beam consisting of some preliminarily strained parts (layers). The mathematical statements of these problems are given, the methods of solution are presented, and some results of solution are shown. The impact of preliminary strains on the state of stresses and strains is investigated, and nonlinear effects are analyzed.
Keywords: nonlinear theory of elasticity, superposition of large strains, prestrained bodies, exact analytical solutions, additive technologies.
Bibliography: 25 titles.
1. Введение
В статье приведены и развиты разработанные совместно е профессором МГУ им. М.В. Ломоносова В.Ан. Левиным подходы к точному аналитическому решению задач о больших деформациях составных изделий (тел) из несжимаемых изотропных нелинейно-упругих материалов, части которых предварительно деформированы. Решение этих задач представляет интерес при анализе напряжений в элементах конструкций, изготавливаемых из предварительно нагруженных частей. Результаты могут быть использованы для тестирования промышленного программного обеспечения, предназначенного для численного моделирования аддитивных технологий.
Постановка задач осуществляется на основе теории наложения больших деформаций [1, 2, 3, 4, 5] и в рамках этой теории может быть сформулирована следующим образом. Части изделия, первоначально не связанные между собой, подвергаются начальному деформированию и переходят в промежуточное состояние. Затем эти части соединяются между собой. Соединение происходит по некоторым поверхностям, общим для каждой пары соединяемых частей. Далее тело, составленное из нескольких частей, деформируется как единое целое под действием приложенной к нему дополнительной нагрузки и переходит в конечное состояние. Предполагается, что на поверхностях, по которым соединены части тела, выполняются условия идеального
контакта, т.е. векторы перемещений в соединяемых частях изделия на этих поверхностях совпадают.
Точные решения для изотропных несжимаемых материалов найдены с использованием известных универсальных решений [6, 7, 8] и могут быть рассмотрены как обобщение этих решений на случай наложения больших деформаций,
2. Задача Ламе-Гадолина
Рассматривается для случая конечной плоской деформации задача о напряженно-деформированном состоянии в двух полых круговых упругих цилиндрах (трубах), один из которых был предварительно деформирован и вставлен в другой цилиндр. Эта задача может быть рассмотрена как обобщение задачи Гадолина [10] на случай наложения конечных деформаций [2]. Отметим, что решение этой задачи при малых деформациях хорошо известно, Это решение приведено, например, в [11].
Постановка задачи в общем случае следующая. Первый цилиндр, к внутренней поверхности которого может быть приложено давление, деформируется под действием этого давления и некоторой нагрузки, приложенной к внешней поверхности, таким образом, что его границы сохраняют круговую форму, а внешний радиус уменьшается на заданную величину. Затем этот цилиндр вставляется во второй (недеформированный) цилиндр, внутренняя поверхность которого совпадает с внешней поверхностью первого цилиндра после предварительной деформации. Цилиндры находятся в промежуточном состоянии (для второго цилиндра оно совпадает с недеформированным), Затем нагрузка, приложенная к внешней поверхности первого цилиндра, «мгновенно» убывает до нуля. Цилиндры приходят во взаимодействие. Вследствие этого в них возникают дополнительные деформации и соответствующие им напряжения. Затем возможно изменение давления, приложенного к внутренней поверхности первого цилиндра. Это также приводит к возникновению дополнительных деформаций и напряжений. Границы цилиндров меняются, и они переходят в конечное (текущее) состояние.
Результаты решения этой задачи могут быть использованы при расчете остаточных напряжений, возникающих при запрессовке полых цилиндров, и при исследовании их возможного разрушения под действием этих напряжений на основе различных критериев [14], Результаты могут быть также применены для тестирования САЕ-систем на основе метода конечных элементов, учитывающих эффекты, связанные с перераспределением больших деформаций,
2.1. Математическая постановка задачи.
В данной задаче можно различать три состояния системы цилиндрических тел: начальное, промежуточное и конечное. Для обозначения этих со-
стояний далее используются индексы 0, 1 и 2 соответственно.
Математическая постановка задачи в координатах промежуточного состояния в обозначениях теории наложения больших деформаций [2] имеет следующий вид:
— для промежуточного состояния (т.е. для определения начальной деформации первого цилиндра) записывается уравнение равновесия
V ■ &
0,1
0,
граничные условия
N ■ °
0,1
— РN (г 1 - Го)|г1
иг,
Г1
и кинематическое соотношение
1
Ф
о,1 —
(Уго)
(1)
(2)
(3)
— для конечного состояния (т.е. для задачи о совместном деформировании цилиндров после того, как один из них вставлен в другой) записывается уравнение равновесия
У" [(Ф1,2) ■ *0,2] — 0
(4)
граничные условия
N ■ а
0,2
—РN N■а
г1
и кинематические соотношения
0,2
г2 1 3
11 Ф1,2 = VГ2 — I + У И2, Фо,2
Фо,1 ■ Ф12.
(5)
(6)
Постановка задачи для обоих состояний включает также определяющие соотношения для материала Бартенева-Хазановича [9]
2^Го,„1/2 — Ро,п1, Го,„
Ф0,га ' Фо,га,
П — 1, 2
(7)
и условия несжимаемости
Ф
к,п
1,
П
1,2, к<п.
(8)
Кроме того, в постановку задачи входит условие непрерывности вектора
2
г2 и вектора нормальных напряжений N ■ сто,2 на границе Г2 между цилиндрами в текущем состоянии. Здесь использованы следующие обозначения: &о,п _ тензор истинных полных для п-го состояния напряжений; Фк,п _ аффинор деформации при переходе из сосотяпия к в состояние щ Ео,га — мера деформации, соответствующая мере Фингера при переходе из состояния к
1
1
2
0
в состояние щ ип — вектор перемещений при переходе из (п — 1)-го в п-е состояние; Г^ и Г — внутренняя и внешняя граница первого цилиндра в состоянии соответственно (Г2 при к = 1, 2 является и внутренней границей второго цилиндра в соответствующем состоянии); Г3 — внешняя граница второго цилиндра в состоявии; иг — заданный вектор смещения внешней границы первого цилиндра при его начальной деформации; Р — давление, приложенное к внутренней границе первого цилиндра.
Отметим, что вследствие осевой симметрии задачи вектор единичной нор-
2 1
мали к границам цилиндров не меняется при деформации, т.е. N = N. Это позволяет переписать граничные условия (5) и условие непрерывности вектора нормальных напряжений в координатах промежуточного состояния.
Рассмотрим решение задачи в цилиндрической системе координат (г, в, г). Пусть го — радиальная координата частицы в начальном состоянии, г — в промежуточном, р — в конечном. Обозначим через Ь внутренний радиус первого цилиндра в промежуточном состоянии. Его внешний радиус в этом состоянии без ограничения общности можно считать равным 1, Пусть внешний радиус первого цилиндра в начальном состоянии равен 1 + е, а внешний радиус второго цилиндра до деформации — с.
2.2. Нахождение решения.
Для решения применяется известный подход, основанный на использовании универсального решения задачи об осесимметричной деформации полого цилиндра из несжимаемого изотропного материала при больших деформациях [7]. Отличие состоит в том, что в данном случае деформирование осуществляется в два этапа, поэтому необходимо сначала определить начальные деформации, а затем — дополнительные.
Интегрирование условий несжимаемости (8) с учетом осевой симметрии задачи (и в предположении о том, что движение частиц происходит только в радиальном направлении) позволяет получить выражения для го и р через г с точностью до постоянной. После подстановки полученных выражений в кинематические и определяющие соотношения и затем в уравнения равновесия получаются линейные дифференциальные уравнения для множителей Лагранжа ро,п (п = 1, 2). Решение этих уравнений позволяет определить ро,п и затем — напряжения. Постоянные определяются из граничных условий. Сначала определим начальную деформацию первого цилиндра. Условие
несжимаемости (8) при п = 1 может быть записано в виде — —0 = 1, откуда
г аг
го = у/г2ГС. (9)
Диагональные компоненты аффинора начальных деформаций с учетом (9) определяются по формулам
/ т ч / 1 л/Г2 + С ,т . Г Г ,т\ 1
(фоД- =( Ат) =—— • = - = ^^, (Фо,1)„ = 1,
а недиагональные компоненты равны нулю.
Подстановка этих компонент в определяющие соотношения (7) дает:
+С 2pr ,
=---+ Po,l, = r^^ + Po,i. (10)
r v r + С
Наконец, после подстановки выражений (10) в уравнение равновесия
^ [Ki)„] + 1 [Ki)„ - Ki^] = о
получается дифференциальное уравнение для p0 i:
j- (po,i) = 0,
dr
откуда p0)i = В, В — постоянная. Подстановка этого выражения в формулы (10) позволяет получить окончательные выражения для начальных напряжений, Постоянные С и В определяются из граничных условий (2), которые могут быть записаны в виде
roUi = ! + ^ (^o,i)rr|r=6 = -Р.
С
С =2 е + е2. (11)
Найдем теперь дополнительные деформации. Характеристики напряженно-деформированного состояния в первом цилиндре будем помечать верхним индексом int, а во втором цилиндре — индексом ext. Индекс 0,2, указывающий номера состояний (0 — начальное, 2 — конечное), будем для краткости опускать. Тензор дополнительных деформаций Ф^2 будем обозначать через Ф, чтобы отличать его от тензора полных деформаций Ф0,2 = Ф,
л г 1 т 1 р dp
Условие несжимаемости det Ф i 2 = 1 может быть записано в виде — — = 1,
г dr
откуда _
" р = V^2 + Q, (12)
Q — постоянная, Поскольку должно выполняться условие непрерывности р на границе между цилиндрами, значение этой постоянной будет одним и тем же для обоих цилиндров.
Ненулевые компоненты аффинора дополнительных деформаций с учетом (12) определятся следующим образом:
т dp г т р л/г2 + Q т , .
Фгг = ~т = = - =--, = 1. 13
dr у/г2 + Q г г
Для ненулевых компонент аффинора полных деформаций в первом цилиндре имеют место следующие выражения:
т int / т ч ^ Vr2 + 2 е + е2 т int , т . Jr2 + Q
Ф^ = (^0,i)rr$rr = -r-n-- , ФЩ = (^0,i)ööФвв - V
гт.—:—> ^ üü аа-^аа , 9 --
л/г2 + Q Vr2 + 2 е + е2
^int = 1
ZZ
Во втором цилиндре начальные деформации отсутствуют, поэтому полные деформации в нем равны дополнительным.
После подстановки компонент этих тензоров в определяющие соотношения (7) имеем:
int 2 ¡1 -у2 + 2 е + е2 int int 2 ¡i л/г2 + Q ^ ( ,
arr =-, 2 -+ Р , = / 2 2 + Р , (14)
\Jr2 + Q Vr2 + 2 е + е2
ext 2 ßr , ext ext 2 ß \Zf2 + Q , ext /ir\
arr = / 2 , „ + P , = -"- + P . (15)
\Л2 + V Г
Уравнение равновесия (4) для внутреннего цилиндра может быть записано в виде
d 1
jr: №/^rr] + - [с*/Фгг - /Ф*] = о.
Аналогично записывается уравнение равновесия для внешнего цилиндра. Подставляя в эти уравнения выражения (13), (14), (15), получаем дифференциальные уравнения для множителей Лагранжа
d d
-т- (Pext) = 0, -т- (pint) = 0, dr аг
откуда
pint = Bint, pext = Bext, (16)
Вint и Bext — постоянные. Подстановка этих выражений в формулы (14), (15) позволяет получить окончательные выражения для напряжений. Постоянные Q, Вint и Bext определяются из граничных условий (5) и условия непрерывности нормальных напряжений на границе между цилиндрами, С учетом сделанного выше замечания о неизменности единичых нормалей к цилиндрическим поверхностям при деформации эти условия могут быть записаны в виде
int I _ _р ^ext I _ о л-ext I _ л-int I
U rr \r=b Г , u rr |r=c 0, u rr |r=l u rr | r=l.
После подстановки в эти условия выражений (14), (15), (16) они запишутся в виде
2 » + 62 + ßint = -р, + в ext = 0, (17)
2 » + 62 + ßint - - В ext = 0. (18)
/Г+Q -TTQ к J
Выражая ßint и Bext из уравнений (17) и подставляя полученные выражения в (18), получаем уравнение
-1 + 2 е + е2 -b2 + 2 е + е2 1
2 ß
+
V- + Q V^Tq V-+Q л/с^Тя
р,
из которого может быть найдена постоянная Q. Решив это нелинейное уравнение численно, можно найти Ви ВеХ из уравнений (17),
2.3. Результаты числовых расчетов.
Рассмотрим некоторые числовые результаты. Отметим, что в данной задаче более естественно задавать в качестве исходного данного внутренний радиус первого цилиндра не в промежуточном, а в начальном состоянии. Этот радиус (обозначим его Ьо) связан с ради усом Ь в промежуточном состоянии соотношением
Ь = л/Ъв2 - 2 е - е2,
которое можно получить, подставив г = Ь и го = Ьо в (9) с учетом (11).
Для тонкостенных цилиндров наблюдается немонотонная зависимость напряжений в первом цилиндре от параметра е, т.е. при некотором значении е напряжения с дальнейшим ростом е начинают убывать. В качестве примера рассмотрим случай с = 1.1, Ьо = (1 + е)/1.1. Для этого случая па рис. 1 дана зависимость напряжения агг на границе между цилиндрами от е, а на рис. 2 показаны зависимости напряжения аее на границах цилиндров от е. Сплошная линия на этом рисунке соответствует напряжению на внутренней границе первого цилиндра, пунктирная линия с длинными штрихами — напряжению на его внешней границе (эти два напряжения взяты с противоположным знаком), штрихпунктирная линия — напряжению на внутренней границе второго цилиндра, пунктирная линия с короткими штрихами — напряжению на его внешней границе.
[ЛЬ]
Рис. 1: Зависимость напряжения агг на границе между цилиндрами от е для материала Бартенева-Хазановича при с = 1.1, Ьо = (1 + е)/1.1. Напряжение взято с противоположным знаком
На всех рисунках напряжения отнесены к модулю р.
[ЛЬ]
2.5 т
1.5 т й 0.5-
Рис. 2: Зависимости напряжения аее на границах цилиндров от е для материала Бартенева-Хазановича при с = 1.1 Ьо = (1 + е)/1.1
3. Кручение составного цилиндра с предварительно деформированным включением
В данном параграфе приведено решение статической задачи о кручении составного цилиндра из несжимаемых нелинейно-упругих материалов. Цилиндр содержит центральное круговое цилиндрическое включение, которое предварительно скручено и растянуто (или сжато) вдоль оси и скреплено с ненапряженным внешним полым цилиндром. Точное аналитическое решение этой задачи для произвольных нелинейно-упругих изотропных несжимаемых материалов получено [12, 13] на основе универсального решения для несжимаемого материала — решения задачи о кручении, растяжении и изменении диаметра круглого цилиндра [15, 6, 7, 16]. Детальный анализ решения выполнен для случая, когда цилиндры изготовлены из материалов типа Муни
Постановка задачи в рамках теории наложения больших деформаций может быть сформулирована следующим образом. В начальном состоянии во внутреннем цилиндре отсутствуют напряжения и деформации. Этот цилиндр подвергается начальному деформированию — скручиванию и осевому растяжению или сжатию, которое вследствие несжимаемости материала цилиндра сопровождается изменением его диаметра. В результате этой деформации цилиндр переходит в промежуточное состояние, в котором он находится в равновесии. Далее этот цилиндр вставляется в полый недеформированный цилиндр, внутренний радиус которого совпадает с радиусом первого (сплошного) цилиндра. Цилиндры скрепляются и подвергаются совместной дополнительной деформации — скручиванию и осевому растяжению или сжатию,
0
0.2 0.4 0.6 0.8 1 1.2 1.4
[17]
которое, как и ранее, сопровождается изменением диаметров цилиндров. Составной цилиндр переходит в конечное состояние, в котором он находится в равновесии. Боковая поверхность составного цилиндра в этом состоянии считается свободной от нагрузок, а к его основаниям приложены осевая сила и крутящий момент (рис, 3),
Рис, 3: Схема нагружения составного цилиндра, ф — осевая сила, С — крутящий момент.
Решение этой задачи позволяет, в частности, исследовать эффект Пойн-тинга, который состоит в изменении длины цилиндрического тела при кручении [18]; исследование этого эффекта для однородного цилиндра из несжимаемого материала выполнено в [19],
3.1. Основные обозначения.
Для краткости будем опускать в этом параграфе нижние индексы у тензоров, описывающих переход из начального в конечное состояние, <у = ^о,2 _ тензор напряжений Коши, Е = Е0,2 — мера деформации Фингера, Ф = Ф0,2 — аффинор полной деформации,
Ф;пй = Ф0,1 _ аффинор начальной деформации, соответствующей переходу из начального состояния в промежуточное,
Фааа = Ф1,2 _ аффинор дополнительной деформации, т.е. деформации при переходе из промежуточного состояния в конечное; 1тА = А : I — след тензора второго ранга А, И. — радиус-вектор частицы в конечном состоянии,
в
в
xs (s = 1, 2, 3) — декартовы координаты частиц в начальном (недеформи-рованном) состоянии, х'т (т =1, 2, 3) — в промежуточном, Х^ (к = 1, 2, 3) — в конечном,
ifc — координатные орты декартовых координат,
р — давление в несжимаемом теле, не выражаемое через деформацию, р, в, ( — цилиндрические координаты частиц внутреннего цилиндра в естественном (недеформированном) состоянии;
г, <р, z — цилиндрические координаты частиц внутреннего цилиндра в промежуточном (предварительно напряженном) состоянии и частиц внешнего цилиндра в естественном (недеформированном) состоянии;
R, Ф, Z — цилиндрические координаты частиц цилиндров в конечном состоянии.
Указанные координаты связаны с декартовыми координатами соотношениями
х1 = р cos в, х2 = р sin в, х3 = (
х[ = Г cos tp, х'2 = г sin tp, х'3 = Z,
X1 = R cos Ф, Х2 = R sin Ф, Х3 = Z.
ep, eg, e^, er, e^, ez, вд, e^ e^ — орты, касательные координатным линиям цилиндрических координат:
ep = i1 cos в + i2 sin в, eg = —i1 sin в + i2 cos в, e^ = i3; er = i1 cos ip + i2 sin >p, ev = —i1 sin <p + i2 cos >p, ez = i3; eR = i1 cos Ф + i2 sin Ф, = —i1 sin Ф + i2 cos Ф, e^ = i3.
3.2. Начальная деформация внутреннего цилиндра.
Начальная (предварительная) деформация включения представляет собой осевое растяжение-сжатие и кручение и задается отображением [16, 6, 7, 15]
г = г(р), <р = в + фо(, z = АоС (19)
фо = const, А0 = const
Тензор отвечающий начальной деформации (19), имеет следующее
выражение
dv т
^mit = eP er + - eg ev + фorec ev + Aoec ez (20)
В (19), (20) ф0 a A0 — угол закручивания и кратность осевого удлинения, определяющие начальную деформацию включения. На основании (20) из условия несжимаемости
det Vinit = -—Ао = 1
находим функцию г(р):
г(р) = А0" 2 Р (21)
С учетом (21) выражение (20) принимает вид
i
Фтй — (eper + ee e^) + фогес e^ + Аовс ez (22)
3.3. Деформации и напряжения в составном цилиндре.
Рассмотрим задачу о кручении и растяжении-сжатии составного цилиндра, Цилиндрический стержень, подвергнутый описанной выше начальной деформации, вставляется в полый круговой цилиндр, находящийся в естественном, ненапряженном состоянии. После скрепления центрального включения с внешним цилиндром составной цилиндр подвергается кручению и растяжению-сжатию как одно целое силами, приложенными к торцам цилиндра, Внешний цилиндр в ненапряженном состоянии занимает область 0 ^ z ^ I, 0 ^ <р ^ 2ж, г1 ^ г ^ г0. Отсюда находится радиус включения:
0 ^ р ^ Л02г1. За отечетную конфигурацию составного цилиндра примем такую, в которой точки тела имеют цилиндрические координаты г, <р, z. Эта конфигурация является естественной для внешней цилиндрической оболочки и предварительно напряженной для центрального включения. Решение задачи о совместном кручении и осевом растяжении-сжатии составного цилиндра ищем в виде
R — R(r), Ф — <р + фz, Z — Xz, 0 ^ г ^ го (23)
где ф и Л — угол закручивания и кратность удлинения составного цилиндра (ф, X — const). Условие несжимаемости позволяет определить функцию R(r):
R(r) — Л" 2 г (24)
При помощи (23), (24) вычисляется аффинор дополнительной деформации Фас1с1 i
*add — А"1 (ereR + е^еф) + ф Rezеф + Aezez (25)
Аффинор Ф полной деформации для внутреннего цилиндра представляется в виде
дХ дх'
Ф — ФШ ■ Фас1d; Фаас1 — im ifc, Ф;пй — —mis im (26)
Для внешнего цилиндра начальная деформация отсутствует, поэтому для Ф—Ф
Обозначив через R1
и Ro соответственно радиус внутреннего цилиндра и внешний радиус внешнего цилиндра в конечном (деформированном) состоянии, на основании (24) имеем
Ri — R(n) — X" 2 п, Ro — R(ro) — X" 2 ro. (27)
Определяющее соотношение для произвольного изотропного несжимаемого материала имеет вид [6]
* = ЫЛ, ^)Е - Х2(Л, ^)Е-1 - р1, (28)
здесь Е — мера деформации Фингера:
Е = Фт ■ Ф, (29)
дШ дШ
Х1 = 2ж, Х2 = 2ж,
Ш (/1, /2) — функция удельной потенциальной энергии упругого материала, /1; 12 — инварианты деформации: 11 = tгЕ, 12 = tr (Е-1), Далее будем обозначать функции хь Х2; характеризующие механические свойства материала внутреннего цилиндра, через х!и х2; чтобы отличить их от функций и характеризующих механические свойства материала внешнего цилиндра.
Учитывая (29), (26), (25), можно получить выражения меры деформации Е Е-1
цилиндре 0 ^ Я < Я1 и во внешнем полом цилиндре Я1 < Я ^ Я0. Далее, используя (28), можно получить выражение тензора напряжений Коши в составном цилиндре как функции координат Д, Ф, Z в конечном состоянии. Это выражение подставляется в уравнение равновесия, что позволяет найти функцию р = р(Я). Зная эту функцию, можно найти окончательные выражения для напряжений [13], Эти выражения здесь не приводятся ввиду их громоздкости, С использованием этих выражений определяются продольная сила ^ и крутящий момент С, приложенные к каждому основанию цилиндра:
До До
я = 2п аг ШЯ, С = 2п тzфR2dR (30)
3.4. Результаты численных расчетов.
При описании численных результатов будем использовать следующие представления для постоянных материалов:
+ р) р(1 - р) о_ М1 + по _ М1 - М
и =—4—, ^ = —^—, Со =-^-, Су¡2 =-^-,(31)
здесь и р — модули сдвига материалов внутреннего и внешнего цилиндров соответственно, р0ъ ¡3 — безразмерные постоянные. Отметим, что условиями положительности потенциальной энергии деформации являются неравенства [6]
0, ро > 0, -1 < р < 1, -1 < ро < 1.
Численные расчеты были выполнены при следующих значениях параметров: модуль сдвига включения р0 = 2^, р1 = р = -0.24 г0 = 1.75г1, угол начального закручивания ф0 = (г1)-1, кратность начального осевого удлинения
А0 = 2.
Была исследована зависимость угла закручивания -ф от кратности осевого удлинения А при полной деформации для случая, когда осевая сила равна 0, Эта зависимость приведена на рис, 4 для разных значений А0, Угол закручивания отнесен к величине (п)-1.
Из рисунка можно видеть, что в составном цилиндре имеет место эффект Пойнтинга, т.е. кручение цилиндра сопровождается его осевым растяжением или сжатием. Если А0 < 1, то при сравнительно малых по модулю углах закручивания имеет место осевое сжатие (А < 1), а в остальных случаях — растяжение. Если же А0 > 1, то при любых углах закручивания имеет место осевое растяжение.
Рис, 4: Зависимость угла закручивания -ф от кратности осевого удлинения А для разных значений А0 (эти значения указаны под рисунком).
Как отмечено в [20, 21], знак эффекта Пойнтинга в нелинейно-упругих материалах зависит от того, выполняются ли для этих материалов дополнительные неравенства, в частности, так называемое эмпирическое неравенство, Указано, что если материал несжимаемый и удовлетворяет эмпирическому неравенству, то для такого материала возможен только положительный эффект Пойнтинга, т.е. при кручении однородного стержня, изготовленного из такого материала, длина стержня увеличивается. Для материала типа Муни эмпирические неравенства совпадают с условиями положительности потенциальной энергии деформации, т.е. ^ > 0 —1 < ¡3 < 1, Таким образом, при выполнении этих неравенств в стержне из материала типа Муни возможен только положительный эффект Пойнтинга, т.е. А > 1, А для
составных предварительно нагруженных стержней из материалов типа Му-ни, как видно из рис, 4, возможно как удлинение, так и укорочение. Это нетривиальный эффект, вызванный предварительной деформацией внутренней части стержня,
4. Изгиб составного бруса с предварительно деформированными слоями
В данном параграфе решена плоская задача о больших деформациях изгиба составного бруса, состоящего из нескольких предварительно деформированных частей (слоев). Получено [23] точное аналитическое решение этой задачи для случая, когда брус изготовлен из изотропных нелинейно-упругих несжимаемых материалов. Решение является универсальным для данного класса материалов и обобщает известное универсальное решение для однородного бруса [6, 7, 16],
Механическая постановка задачи следующая. Брусья в форме прямоугольного параллелепипеда из несжимаемого материала, находящиеся в ненапряженном состоянии, подвергаются начальной деформации — растяжению или сжатию — и переходят в промежуточное состояние. Начальная деформация каждого бруса однородна, но различна для разных брусьев. Для некоторых брусьев начальная деформация может отсутствовать. Далее эти брусья соединяются между собой (рис, 5), Предполагается, что на гранях, по которым происходит соединение, выполняются условия идеального контакта, т.е. совпадают перемещения соединяемых брусьев, К торцам составного бруса прикладывается изгибающий момент. Далее составной брус деформируется как единое целое. Под действием изгибающего момента и вследствие начальной деформации слоев происходит его изгиб. Брус переходит в конечное состояние,
4.1. Основные уравнения постановки задачи
Используются те же основные обозначения, что и в предыдущем параграфе.
Постановка задачи осуществляется в координатах промежуточного состояния и включает следующие уравнения.
Уравнение равновесия (массовые силы отсутствуют)
11
V ■ Б = 0. (32)
к
Здесь Б — несимметричный тензор напряжений Пиолы, определенный в конфигурации к-то состояния. Этот тензор выражается через тензор истинных напряжений а следующим образом:
Б =(ае1 ¥аМ) (Ф^) 1 ■ а.
(33)
Упругий потенциал W и определяющие соотношения для несжимаемых изотропных нелинейно-упругих материалов |6|:
W = W (Д, /2) , (34)
D = [Xi (h, /2) + hX2 (h, /2)]^fnit ■ Ф - X2 (h, /2) ■ G ■ Ф - p (Ф^)-1,
(35)
dW dW 1
Xi = 2 X2 = 2 -h = tr G, /2 = tr (G ). (36)
0I1 0I2
Здесь /1; I2 — инварианты деформации, хъ Х2 _ функции отклика материала, р — множитель Лагранжа, определяющий давление в несжимаемом теле.
Условия несжимаемости материала:
det Фinit = det Фааа = det Ф = 1. (37)
Кинематические соотношения:
G = Ф ■ Фт, Ф = Ф^ ■ Фааа. (38)
Рассмотрим плоскую деформацию прямоугольного составного бруса, занимающего в промежуточном состоянии область 0 ^ х ^ /, 0 ^ у ^ h (рис. 5). Через ж, у, г обозначены декартовы координаты частицы бруса в промежуточном состоянии. Размер по координате г, один и тот же для всех слоев, значения не имеет. Брус состоит из п слоев, к-й слой занимает в промежуточном состоянии область hk-1 ^ у ^ hk (0 = h0 ^ h1 ^ ... ^ hn-1 ^ hn = h). Начальная (предварительная) деформация fc-ro слоя (к = 1,...,п) относительно начального ненапряженного состояния материала однородна и задается следующим аффинором деформации (A(k) = const):
Фтй = A(fc) i1 i1 + A-1 i2 i2 + i313, (39)
где ij — координатные орты.
I
Рис. 5: Форма составного бруса в промежуточном состоянии
Если к-й слой бруса не был предварительно напряжен, то промежуточное состояние этого слоя совпадает с начальным, А(&) = 1 и Фщ^ = I-
Материалы слоев составного бруса могут различаться. Будем далее обозначать функцию потенциальной энергии деформации &-го слоя {Кк-1 ^ у < < Кк) через Ш(к\ а соответствующие функции отклика через х!^ и
Л*) = 2 ^^ Л*) = 2 ^^ ' '
Х1 =2 011 'Х2 =2 д12 ■
Отметим, что для расматриваемого класса задач 11 = 12 вследствие условия несжимаемости и условия плоской деформации.
4.2. Нахождение решения задачи.
Дополнительную деформацию составного бруса будем разыскивать в виде, соответствующем нелинейному чистому изгибу |22|
X = р (у) sin fix, Y = р (у) cos fix, Z = z, fi = const. (40)
где X,Y,Z — декартовы координаты точек бруса в конечном состоянии. При деформации (40) каждая прямая у = const превращается в дугу окружности с радиусом р (у). Геометрический смысл параметра fi состоит в том, что fix — угол поворота сечения х = const вокруг оси, ортогональной плоскости бруса. Форма составного бруса в конечном состоянии показана па рис. 6. Отметим, что при fi > 0 брус изгибается выпуклостью вверх, а при fi < 0 — выпуклостью вниз.
Рис. 6: Форма составного бруса в конечном состоянии
Аффинор дополнительной деформации, отвечающий отображению (40), имеет вид
п • dp.
Фааа = fip ii ei + —e2 + 13 13, ^
в! = 11 cos fix — i2 sin fi x, e2 = 11 sin fix + i2 cos fi x.
Для к-го слоя бруса в соответствии с формулами (38), (39), (41)
d п
Ф = \(k)fip 1i ei + A7fc1-^ 12 в2 + 1з 1з. (42)
Вычисляемая па основании (38), (41), (42) мера деформации Коши-Грипа для &-го слоя бруса (в области < у < К&) имеет следующий вид:
G = \к)2/3У i1 i1 + A(fc)-^d^J i2 i2 + 1з i3
Эта мера деформации не зависит от координаты ж, что обеспечивает сведение исходной двумерной задачи к одномерной.
Функция р (у) находится из условия несжимаемости
^ = 1. (44)
Уравнение (44) имеет следующее решение:
р2 = р0 + 2/3-1 у, ро = р(0), (45)
где р0 — неизвестная постоянная.
Подстановка полученных выражений в определяющие соотношения позволяет найти напряжения. Подставив выражения для напряжений в урав-
нахождения постоянных интегрирования используются граничные условия
^22)(0) = 0, ^2?(h) = 0, (46)
обозначающие отсутствие нагружения оснований составного бруса, и условия сопряжения на поверхностях раздела слоев
^22-1)(hk) = D(2)(hk), к = 1,...,п - 1. (47)
Для того, чтобы удовлетворить краевым условиям на концах бруса х = 0, х = I в интегральном смысле, определим главный вектор f и главный момент M сил, действующих в сечении х = const при деформации вида (40), Из уравнений равновесия и граничных условий (46), (47) вытекает, что f=0
висит от точки приведения (полюса). Выбирая за полюс центр окружностей,
= const
hfc
n S
M = M i3, M = M(3) = ^ D^kkpdy. (48)
k=1.
d 2
i^J i2 i2 + 1з 1З. (43)
hk-
Согласно (48) изгибающий момент одинаков во всех сечениях бруса. Это означает, что реализация деформации изгиба (40) требует приложения к концам бруса только изгибающих моментов, величина которых определяется формулой (48), Функция М (¡3) является основной характеристикой, определяющей сопротивление составного бруса при изгибе.
4.3. Результаты численных расчетов
Приведем некоторые результаты численных расчетов для материала Тре-лоара (неогуковекого материала) [6], когда Xi ) = Ук = const = 0, к = 1,..., п.
На рис, 7, 8 даны результаты расчетов для трехслойного бруса (п = 3) при h2 = 3hi; h3 = h = 5hi; у2 = у3 = А(i) = А(3) = 1 (т.е. нижний и верхний слои предварительно не деформированы и имеют одинаковые свойства, а предварительно деформирован только средний слой; кратность А(2) начального удлинения этого слоя в направлении оси х будем обозначать через А),
На рис, 7 показаны зависимости изгибающего момента М от параметра ¡3 для разных значений начального удлинения А, Эти зависимости монотонны и асимптотически стремятся к некоторым предельным значениям, зависящим от А, при ¡3 ^
М Цй?
30
20
10
/
t
/ / .
//'
^у'-р-
^30" -40 -
0.25.....0.5--1.0---1.5--2.0
-1-
0.5
М
Рис. 7: Зависимости изгибающего момента М от параметра @ для разных значений кратности начального удлинения А для материала Трелоара (значения А указаны под рисунком)
На рис. 8 приведена зависимость параметра ¡3 от кратности начального удлинения А для случая М = 0, т.е. при отсутствии внешних нагрузок. Можно видеть, что если средний слой бруса предварительно растянут в направлении оси х (т.е. при А > 1), то брус изгибается выпуклостью вверх, а если этот слой предварительно сжат (А < 1), то брус изгибается выпуклостью вниз.
Отметим также, что зависимость А(3) при отсутствии внешних нагрузок немонотонна и имеет точку максимума при А ~ 1.75, Это означает, что в случае предварительного растяжения среднего слоя существует предельный относительный угол поворота поперечных сечений бруса при отсутствии внешних нагрузок.
Рис, 8: Зависимость параметра /3 от кратности начального удлинения А для материала Трелоара при отсутствии изгибающего момента
Авторы благодарят профессора В.Ан. Левина (МГУ им, М.В, Ломоносова) за совместное обсуждение подходов к постановке задач и результатов их решения,
СПИСОК ЦИТИРОВАННОЙ ЛИТЕРАТУРЫ
1, Levin V.A. Theory of repeated superposition of large deformations. Elastic and viseoelastie bodies// International Journal of Solids and Structures, 1998, V. 35. P. 2585-2600.
2, Левин В.А. Многократное наложение больших деформаций в упругих и вязкоупругих телах. М,: Наука. Физматлит, 1999. 224 с.
3. Левин В,А,, Зингерман К.М, Плоские задачи многократного наложения больших деформаций. Методы решения. М,: Физматлит, 2002. 272 с.
4. V.A. Levin, К.М, Zingerman, A.V. Vershinin, E.I. Freiman, A.V. Yangirova, Numerical analysis of the stress concentration near holes originating in
previously loaded viseoelastie bodies at finite strains// International Journal of Solids and Structures. 2013. Vol. 50. Issue 20-21. P. 3119-3135.
5. Levin V.A., Vershinin A.V. Non-stationary plane problem of the successive origination of stress concentrators in a loaded body. Finite deformations and their superposition// Communications in Numerical Methods in Engineering. 2008. V. 24. P. 2229-2239.
6. Лурье А.И Нелинейная теория упругости. М,: Наука, 1980. 512 с.
7. Трусделл К. Первоначальный курс рациональной механики сплошных сред. М,: Мир, 1975. — 592 с.
8. Rivlin R.S. Large elastic deformations of isotropic materials. VI. Further results in the theory of torsion, shear, and flexure// Philos. Trans. Roy. Soc. London, 1949. A252. P. 173-195.
9. Бартенев Г.М., Хазанович Т.Н. О законе высокоэластичных деформаций сеточных полимеров // Высокомолекулярные соединения. 1960. Т. 2, № 1. С.20-28.
10. Гадолин А.В. Теория орудий, скрепленных обручами. Артиллерийский журнал. 1861, JV2 12.
11. Седов Л.И. Механика сплошной среды. Т. 2. М,: Наука, 1994. 560 с.
12. Левин В.А., Зубов Л.М., Зингерман К.М. Кручение составного нелинейно упругого цилиндра с предварительно напряжённым включением / / Докл. РАН. 2013. Т. 453, № 5. С. 507-510.
13. Levin V.A., Zubov L. М,, Zingerman К. М. The torsion of a composite, nonlinear-elastic cylinder with an inclusion having initial large strains // Int. J. Solids a. Struct. 2014. V. 51. № 6. P. 1403-1409.
14. Левин В.А., Морозов E.M., Матвиенко Ю.Г. (под редакцией В.А. Левина) Избранные нелинейные задачи механики разрушения. М,: Физматлит, 2004. - 407 с.
15. Грин А., Адкинс Дж. Большие упругие деформации и нелинейная механика сплошной среды. М,: Мир, 1955. 445 с.
16. Rivlin R.S. Large elastic deformations of isotropic materials // Philos. Trans. Roy. Soc. London, 1948. A240. P. 459-508.
17. Moony M.A. Theory of large elastic deformation // Journal of Applied Physics. 1940. № 11. P. 582-592.
18. Povnting J. H. On pressure perpendicular to the shear planes in finite pure shears, and on the lengthening of loaded wires when twisted// Proc. R. Soc. Lond. A 1909 82, p. 546-559.
19. Зубов . I.M. О прямом и обратном эффекте Пойнтинга в упругих цилиндрах // Докл. РАН. 2001. Т. 380, № 2. С. 194-196.
20. Mihai L.A., Goriely A. Positive or negative Povnting effect? The role of adscititious inequalities in hvperelastie materials // Proe. E. Soe. A. 2011. V. 467. P. 3633-3646.
21. Moon H,, Truesdell C. Interpretation of adscititious inequalities through the effects pure shear stress produces upon an isotropic elastic solid // Arch. Ration. Mech. Anal. 1974. V. 55. P. 1—17.
22. Зеленина А. А, Зубов Л.М. Нелинейная теория чистого изгиба призматических упругих тел// Прикл. мат. и мех. 2000. Т. 64. Вып. 3. С. 416-424.
23. Левин В.А., Зубов Л.М., Зингерман К.М. Точное решение задачи о нелинейном изгибе составного бруса с предварительно деформированным слоем при конечных деформациях // Докл. РАН. 2015. Т. 460, JV2 2. С. 155-158.
24. Levin V.A., Zubov L. М,, Zingerman К. М. An exact solution for the problem of flexure of a composite beam with preliminarily strained layers under large strains // Int. J. Solids a. Struct. 2015. V. 67-68. P. 244—249.
25. Levin V.A,, Zubov L.M., Zingerman K.M. An exact solution for the problem of flexure of a composite beam with preliminarily strained layers under large strains. Part 2. Solution for different types of incompressible materials // International Journal of Solids and Structures. 2016. V.100-101. P.558-565.
REFERENCES
1. Levin V.A. Theory of repeated superposition of large deformations. Elastic and viseoelastie bodies. International Journal of Solids and Structures. 1998. V. 35. P. 2585-2600.
2. Levin, VA 1999, Mnogokratnoe nalozhenie bol'shih, deformaciy v uprugih i vyazkouprugih telah [Repeated superposition of large deformationsin elastic and viseoelastie bodies], Fizmatlit, Moscow.
3. Levin, VA & Zingerman, KM 2002, Ploskie zadaehi mnogokratnogo nalozheniya bol'shih, deformaciy. Metody resheniya. [Plane problems of repeated superposition of large strains. Methods of solution] Fizmatlit, Moscow.
4. Levin, V.A,, Zingerman, K.M., Vershinin, A.V,, Freiman, E.I., Yangirova, A.V, 2013. "Numerical analysis of the stress concentration near holes originating in previously loaded viseoelastie bodies at finite strains" Int. J. Solids and Structures. Vol. 50. Issue 20-21. P. 3119-3135.
5. Levin V.A., Vershinin A.V. 2008, "Non-stationary plane problem of the successive origination of stress concentrators in a loaded body. Finite deformations and their superposition" Communications in Numerical Methods in Engineering. Vol. 24. P. 2229-2239.
6. Lurie, AI. 1980, Nelineynaya teoriya uprugosti [Nonlinear theory of elasticity] Nauka, Moscow.
7. Truesdell, C. 1972, A First Cource in Rational Continuum Mechanics. Baltimore, Maryland, The Johns Hopkins University.
8. Rivlin R.S. 1949. "Large elastic deformations of isotropic materials. VI. Further results in the theory of torsion, shear, and flexure." Philos. Trans. Roy. Soc. London , Vol. A252. P. 173-195.
9. Bartenev, G.M., Khazanovich, T.N. 1960. "The law of high-elastic strains of mesh polymers". Vysokomolekulyarnye soedineniya. Vol. 2. .V" 1. P. 20-28.
10. Gadolin, A.V., 1861. "Theory of cannons fastened by bands." Artilleriyskiy zhurnal. 1861, № 12.
11. Sedov, LI 1994, Mekhanika Sploshnoy Sredy [Mechanics of Continuous Medium] Vol. 2. Nauka, Moscow.
12. Levin, V.A., Zubov, L.M. & Zingerman, K.M. 2013 "Torsion of a composite nonlinear elastic cylinder with a prestressed inclusion," Doklady Physics, Vol. 58, No. 12, pp. 540-543.
13. Levin V.A., Zubov L.M. & Zingerman K.M. 2014. "The torsion of a composite, nonlinear-elastic cylinder with an inclusion having initial large strains" Int. J. Solids a. Struct. Vol. 51. № 6. P. 1403-1409.
14. Levin, VA, Morozov, EM & Matvienko Yu.G. (Ed. by V.A. Levin). 2004, Izbrannye Nelineinye Zadachi Mehaniki Razrusheniya [Selected Nonlinear Problems of Fracture Mechanics] Fizmatlit, Moscow.
15. Green, A.E., & Adkins, J.E., 1960. Large Elastic Deformations and Non-linear Continuum Mechanics. Oxford, Clarendon Press.
16. Rivlin R.S. 1948. "Large elastic deformations of isotropic materials" Philos. Trans. Roy. Soc. London, Vol. A240. P. 459-508.
17. Moony, M.A. 1940. "Theory of large elastic deformation," Journal of Applied Physics. № 11. P. 582-592.
18. Povnting J. H. 1909, "On pressure perpendicular to the shear planes in finite pure shears, and on the lengthening of loaded wires when twisted," Proc. R. Soc. Lond. A. Vol. 82, pp. 546-559.
19. Zubov, L.M., 2001, "Direct and inverse Povnting effects in elastic cylinders," Doklady Physics Vol, 46, pp. 675—677,
20. Mihai, L.A., Gorielv, A,, 2011, "Positive or negative Povnting effect? The role of adscititious inequalities in hvperelastie materials," Proc. R. Soc. A. V. 467. P. 3633-3646.
21. Moon, H,, Truesdell, C. 1974. "Interpretation of adscititious inequalities through the effects pure shear stress produces upon an isotropic elastic solid," Arch. Ration. Mech. Anal. Vol. 55. P. 1—17.
22. Zelenina, A.A., & Zubov, L.M., 2000. "The non-linear theory of the pure bending of prismatic elastic solids," J. Appl. Math. Mech. Vol. 64, pp. 399406.
23. Levin, V.A,, Zubov L.M., Zingerman, K.M., 2015. "Exact solution of the nonlinear bending problem for a composite beam containing a prestressed layer at large strains," Doklady Physics. Vol, 60, No, 1, P. 24-27,
24. Levin V.A,, Zubov L.M., & Zingerman, K.M, 2015, "An exact solution for the problem of flexure of a composite beam with preliminarily strained layers under large strains," Int. J. Solids a. Struct. Vol. 67—68. P. 244—249.
25. Levin V.A,, Zubov L.M., & Zingerman, K.M, 2016. "An exact solution for the problem of flexure of a composite beam with preliminarily strained layers under large strains. Part 2. Solution for different types of incompressible materials," International Journal of Solids and Structures. Vol.100-101. P.558-565.
Тверской государственный университет.
Южный федеральный университет.
Получено 10.07.2017
принято в печать 14.09.2017
