Теория конечных деформаций - Finite strain theory

В механика сплошной среды, то теория конечных деформаций-также называемый теория больших деформаций, или же теория больших деформаций-имеет дело с деформации в которых деформации и / или вращения достаточно велики, чтобы опровергнуть допущения, присущие теория бесконечно малых деформаций. В этом случае недеформированная и деформированная конфигурации континуума существенно различаются, что требует четкого различия между ними. Обычно это случается с эластомеры, пластически деформирующий материалы и другие жидкости и биологический мягких тканей.

Смещение

Рис. 1. Движение сплошного тела.

Смещение тела состоит из двух компонентов: жесткое тело смещение и деформация.

  • Смещение твердого тела состоит из одновременного перевод (физика) и вращение тела без изменения его формы или размера.
  • Деформация подразумевает изменение формы и / или размера тела по сравнению с исходной или недеформированной конфигурацией. в текущую или деформированную конфигурацию (Рисунок 1).

Изменение конфигурации сплошного тела можно описать поле смещения. А поле смещения это векторное поле всех векторов смещения для всех частиц в теле, что связывает деформированную конфигурацию с недеформированной конфигурацией. Расстояние между любыми двумя частицами изменяется тогда и только тогда, когда произошла деформация. Если смещение происходит без деформации, то это смещение твердого тела.

Материальные координаты (лагранжевое описание)

Смещение частиц индексируется переменной я можно выразить следующим образом. Вектор, соединяющий положения частицы в недеформированной конфигурации и деформированная конфигурация называется вектор смещения. С помощью на месте и на месте , оба из которых являются векторами от начала системы координат до каждой соответствующей точки, мы имеем Лагранжево описание вектора смещения:

Где являются ортонормированными единичные векторы которые определяют основа пространственной (лабораторной) системы координат.

Выраженное в терминах материальных координат поле смещения:

Где вектор смещения, представляющий перенос твердого тела.

В частная производная вектора смещения относительно материальных координат дает тензор градиента смещения материала . Таким образом, мы имеем

куда это тензор градиента деформации.


Пространственные координаты (описание Эйлера)

в Эйлерово описание, вектор, выходящий из частицы в недеформированной конфигурации до ее расположения в деформированной конфигурации называется вектор смещения:

Где являются единичными векторами, которые определяют основу системы координат материала (тело-рама).

Выраженное в терминах пространственных координат поле смещения:

Частная производная вектора смещения по пространственным координатам дает тензор градиента пространственного смещения . Таким образом, мы имеем

Взаимосвязь материальной и пространственной систем координат

являются направляющие косинусы между материальной и пространственной системами координат с единичными векторами и , соответственно. Таким образом

Отношения между и тогда дается

Знаю это

тогда

Объединение систем координат деформированной и недеформированной конфигураций

Обычно системы координат для деформированной и недеформированной конфигураций накладываются друг на друга, что приводит к , а направляющие косинусы принимают вид Дельты Кронекера, т.е.

Таким образом, в материальных (недеформированных) координатах смещение может быть выражено как:

А в пространственных (деформированных) координатах смещение может быть выражено как:

Тензор градиента деформации

Рис. 2. Деформация сплошного тела.

Тензор градиента деформации связана как ссылки и текущей конфигурации, как видно единичных векторов и , поэтому это двухточечный тензор.

В силу предположения о непрерывности , имеет обратное , куда это тензор градиента пространственной деформации. Затем по теорема о неявной функции,[1] в Якобиан детерминант должно быть неособый, т.е.

В тензор градиента деформации материала это тензор второго порядка который представляет градиент функции отображения или функционального отношения , который описывает движение континуума. Тензор градиента деформации материала характеризует локальную деформацию в материальной точке с вектором положения , т.е. деформацию в соседних точках, преобразованием (линейное преобразование ) элемент линии материала, исходящий из этой точки из эталонной конфигурации в текущую или деформированную конфигурацию, при условии непрерывности в функции отображения , т.е. дифференцируемая функция из и время , откуда следует, что трещины а пустоты не открываются и не закрываются во время деформации. Таким образом, мы имеем

Вектор относительного смещения

Рассмотрим частица или материальная точка с вектором положения в недеформированной конфигурации (рисунок 2). После перемещения тела новое положение частицы, обозначенное в новой конфигурации задается векторным положением . Системы координат для недеформированной и деформированной конфигурации могут быть совмещены для удобства.

Рассмотрим теперь материальную точку соседний , с вектором положения . В деформированной конфигурации эта частица имеет новое положение заданный вектором положения . Предполагая, что отрезки линии и соединение частиц и как в недеформированной, так и в деформированной конфигурации, соответственно, очень малы, то мы можем выразить их как и . Таким образом, из рисунка 2 мы имеем

куда это вектор относительного смещения, который представляет собой относительное смещение относительно в деформированном состоянии.

Приближение Тейлора

Для бесконечно малого элемента , и в предположении непрерывности поля смещения можно использовать Расширение ряда Тейлора вокруг точки пренебрегая членами более высокого порядка, чтобы аппроксимировать компоненты вектора относительного смещения для соседней частицы в качестве

Таким образом, предыдущее уравнение можно записать как

Производная по времени от градиента деформации

Расчеты, которые включают зависящую от времени деформацию тела, часто требуют вычисления производной по времени от градиента деформации. Геометрически непротиворечивое определение такой производной требует экскурсии в дифференциальная геометрия[2] но мы избегаем этих проблем в этой статье.

Производная по времени от является

куда - скорость. Производная в правой части представляет собой градиент скорости материала. Обычно это преобразовывают в пространственный градиент, т. Е.

куда это градиент пространственной скорости. Если градиент пространственной скорости постоянен, указанное выше уравнение может быть решено точно, чтобы дать

предполагая в . Есть несколько методов вычисления экспоненциальный над.

Связанные величины, часто используемые в механике сплошных сред: тензор скорости деформации и спиновый тензор определяется, соответственно, как:

Тензор скорости деформации показывает скорость растяжения элементов линии, а тензор спина указывает скорость вращения или завихренность движения.

Материальная производная по времени от обратного градиента деформации (сохраняя фиксированную исходную конфигурацию) часто требуется в анализах, которые включают конечные деформации. Эта производная

Вышеупомянутое соотношение можно проверить, взяв материальную производную по времени от и отмечая, что .

Преобразование поверхности и элемента объема

Чтобы преобразовать величины, которые определены относительно площадей в деформированной конфигурации, в те, которые относятся к площадям в эталонной конфигурации, и наоборот, мы используем соотношение Нансона, выраженное как

куда - площадь области в деформированной конфигурации, это та же область в эталонной конфигурации, и является внешней нормалью к элементу площади в текущей конфигурации, а является внешней нормалью в эталонной конфигурации, это градиент деформации, и .

Соответствующая формула преобразования элемента объема:

Полярное разложение тензора градиента деформации

Рис. 3. Представление полярного разложения градиента деформации.

Градиент деформации , как и любой обратимый тензор второго порядка, можно разложить с помощью полярное разложение теоремы, в произведение двух тензоров второго порядка (Truesdell, Noll, 1965): ортогонального тензора и положительно определенного симметричного тензора, т. е.

где тензор это собственный ортогональный тензор, т.е. и , представляющий вращение; тензор это правый тензор растяжения; и в левый тензор растяжения. Условия верно и оставили означает, что они находятся справа и слева от тензора вращения , соответственно. и оба положительно определенный, т.е. и для всех ненулевых , и симметричные тензоры, т.е. и , второго порядка.

Из этого разложения следует, что деформация линейного элемента в недеформированной конфигурации на в деформированной конфигурации, т.е. , может быть получен либо путем предварительного растяжения элемента на , т.е. с последующим вращением , т.е. ; или что то же самое, применяя жесткое вращение во-первых, т.е. с последующим растяжением , т.е. (См. Рисунок 3).

Из-за ортогональности

так что и имеют то же самое собственные значения или же основные участки, но разные собственные векторы или же основные направления и , соответственно. Основные направления связаны между собой

Это полярное разложение, уникальное как обратима с положительным определителем, является следствием сингулярное разложение.

Тензоры деформации

В механике используется несколько не зависящих от вращения тензоров деформации. В механике твердого тела наиболее популярны правые и левые тензоры деформации Коши – Грина.

Поскольку чистое вращение не должно вызывать деформаций в деформируемом теле, часто удобно использовать независимые от вращения меры деформации в механика сплошной среды. Поскольку вращение с последующим его обратным вращением не приводит к изменению () мы можем исключить поворот, умножив своим транспонировать.

Правый тензор деформации Коши – Грина

В 1839 г. Джордж Грин ввел тензор деформации, известный как правый тензор деформации Коши – Грина или же Тензор деформации Грина, определяется как:[4][5]

Физически тензор Коши – Грина дает нам квадрат локального изменения расстояний из-за деформации, т.е.

Инварианты часто используются в выражениях для функции плотности энергии деформации. Наиболее часто используемые инварианты находятся

куда являются коэффициентами растяжения для единичных волокон, которые изначально ориентированы вдоль направлений собственных векторов правого (эталонного) тензора растяжения (обычно они не выровнены по трем осям систем координат).

Тензор деформации пальца

В ИЮПАК рекомендует[5] что обратный к правому тензору деформации Коши – Грина (называемый в этом документе тензором Коши), т. е., , называться Тензор пальца. Однако эта номенклатура не является общепринятой в прикладной механике.

Левый тензор деформации Коши – Грина или Фингера

Изменение порядка умножения в формуле для правого тензора деформации Грина – Коши приводит к левый тензор деформации Коши – Грина который определяется как:

Левый тензор деформации Коши – Грина часто называют Тензор деформации пальца, названный в честь Йозеф Фингер (1894).[5][6][7]

Инварианты также используются в выражениях для функции плотности энергии деформации. Обычные инварианты определяются как

куда - определитель градиента деформации.

Для несжимаемых материалов используется несколько иной набор инвариантов:

Тензор деформации Коши

Ранее в 1828 г.[8] Огюстен Луи Коши ввел тензор деформации, определяемый как обратный к левому тензору деформации Коши – Грина, . Этот тензор также получил название Тензор Пиолы[5] и Тензор пальца[9] в литературе по реологии и гидродинамике.

Спектральное представление

Если есть три различных основных участка , то спектральные разложения из и дан кем-то

Более того,

Заметьте, что

Следовательно, из единственности спектрального разложения также следует, что . Левая полоса () также называется тензор пространственного растяжения при этом правая растяжка () называется тензор растяжения материала.

Эффект действующий на растянуть вектор на и повернуть его в новую ориентацию , т.е.

В том же духе,

Производные от стрейч

Производные растяжения относительно правого тензора деформации Коши – Грина используются для вывода соотношений напряжение-деформация многих твердых тел, в частности гиперупругие материалы. Эти производные

и из наблюдений следует, что

Физическая интерпретация тензоров деформации

Позволять - декартова система координат, определенная на недеформированном теле, и пусть - другая система, заданная на деформируемом теле. Пусть кривая в недеформированном теле параметризовать с помощью . Его образ в деформированном теле .

Недеформированная длина кривой определяется выражением

После деформации длина становится равной

Отметим, что правый тензор деформации Коши – Грина определяется как

Следовательно,

что указывает на то, что изменения длины характеризуются .

Тензоры конечных деформаций

Концепция чего-либо напряжение используется для оценки того, насколько данное смещение локально отличается от смещения твердого тела.[1][10] Одной из таких деформаций при больших деформациях является Лагранжев тензор конечных деформаций, также называемый Тензор деформации Грина-Лагранжа или же Зеленый - тензор деформации Сен-Венана, определяется как

или как функция тензора градиента смещения

или же

Тензор деформации Грина-Лагранжа является мерой того, насколько отличается от .

В Тензор конечных деформаций Эйлера-Альманси, относящийся к деформированной конфигурации, т.е. эйлерово описание, определяется как

или как функция градиентов смещения мы имеем

Seth–Hill family of generalized strain tensors

B. R. Seth от Индийский технологический институт Харагпур was the first to show that the Green and Almansi strain tensors are special cases of a more general strain measure.[11][12] The idea was further expanded upon by Родни Хилл в 1968 г.[13] The Seth–Hill family of strain measures (also called Doyle-Ericksen tensors)[14] можно выразить как

For different values of у нас есть:

The second-order approximation of these tensors is

куда is the infinitesimal strain tensor.

Многие другие определения тензоров допустимы при условии, что все они удовлетворяют условиям, которые:[15]

  • исчезает при всех движениях твердого тела
  • зависимость на тензоре градиента смещения непрерывна, непрерывно дифференцируема и монотонна
  • также желательно, чтобы сводится к тензору бесконечно малых деформаций как норма

Примером может служить набор тензоров

которые не принадлежат классу Сета – Хилла, но имеют то же приближение 2-го порядка, что и меры Сета – Хилла при для любого значения .[16]

Коэффициент растяжения

В коэффициент растяжения - мера растяжения или нормальной деформации дифференциального линейного элемента, которая может быть определена либо в недеформированной конфигурации, либо в деформированной конфигурации.

Степень растяжения дифференциального элемента (Рисунок) в направлении единичного вектора в материальной точке в недеформированной конфигурации определяется как

куда - деформированная величина дифференциального элемента .

Аналогичным образом степень растяжения дифференциального элемента (Рисунок), в направлении единичного вектора в материальной точке в деформированной конфигурации определяется как

Нормальное напряжение в любом направлении может быть выражена как функция степени растяжения,

Это уравнение подразумевает, что нормальная деформация равна нулю, то есть деформации нет, когда растяжение равно единице. Некоторые материалы, такие как эластометры, могут выдерживать коэффициенты растяжения 3 или 4 до того, как они выйдут из строя, тогда как традиционные конструкционные материалы, такие как бетон или сталь, терпят неудачу при гораздо более низких коэффициентах растяжения, возможно, порядка 1,1 (ссылка?)

Физическая интерпретация тензора конечных деформаций

Диагональные компоненты тензора конечных деформаций Лагранжа связаны с нормальной деформацией, например

куда нормальная деформация или инженерная деформация в направлении .

Недиагональные компоненты тензора конечных деформаций Лагранжа связаны с деформацией сдвига, например

куда это изменение угла между двумя линейными элементами, которые изначально были перпендикулярны направлениям и , соответственно.

При определенных обстоятельствах, т.е. при малых смещениях и малых скоростях смещений, компоненты лагранжевого тензора конечных деформаций могут быть аппроксимированы компонентами тензор бесконечно малых деформаций

Тензоры деформации в конвективных криволинейных координатах

Представление тензоров деформации в криволинейные координаты полезен для решения многих задач механики сплошных сред, таких как нелинейные теории оболочек и большие пластические деформации. Позволять обозначают функцию, с помощью которой вектор положения в пространстве строится из координат . Координаты называются "конвекционными", если они соответствуют взаимно однозначному отображению в и из лагранжевых частиц в сплошном теле. Если координатная сетка «нарисована» на теле в его начальной конфигурации, то эта сетка будет деформироваться и течь вместе с движением материала, чтобы оставаться нарисованными на тех же частицах материала в деформированной конфигурации, так что линии сетки пересекаются на одной и той же частице материала. в любой конфигурации. Касательный вектор к кривой линии деформированной координатной сетки в дан кем-то

Три касательных вектора в точке образуют локальную основу. Эти векторы связаны с векторами обратного базиса соотношением

Определим тензорное поле второго порядка (также называемый метрический тензор ) с компонентами

В Символы Кристоффеля первого рода можно выразить как

Чтобы увидеть, как символы Кристоффеля связаны с правым тензором деформации Коши – Грина, давайте аналогичным образом определим два основания: уже упомянутое, касательное к линиям деформированной сетки, а другое - к линиям недеформированной сетки. А именно,

Градиент деформации в криволинейных координатах

Используя определение градиент векторного поля в криволинейных координатах градиент деформации можно записать как

Правый тензор Коши – Грина в криволинейных координатах

Правый тензор деформации Коши – Грина имеет вид

Если мы выразим в терминах компонентов по отношению к базису {} у нас есть

Следовательно,

и соответствующий символ Кристоффеля первого типа может быть записан в следующей форме.

Некоторые соотношения мер деформации и символов Кристоффеля

Рассмотрим взаимно однозначное отображение из к и предположим, что существуют два положительно определенных симметричных тензорных поля второго порядка и это удовлетворяет

Потом,

Отмечая, что

и у нас есть

Определять

Следовательно

Определять

потом

Определим символы Кристоффеля второго типа как

потом

Следовательно,

Из обратимости отображения следует, что

Мы также можем сформулировать аналогичный результат в терминах производных по . Следовательно,

Условия совместимости

Проблема совместимости в механике сплошных сред включает определение допустимых однозначных непрерывных полей на телах. Эти допустимые условия оставляют тело без нефизических зазоров или перекрытий после деформации. Большинство таких условий применимо к односвязным телам. Для внутренних границ многосвязных тел требуются дополнительные условия.

Совместимость градиента деформации

Необходимые и достаточные условия существования совместимого поля над односвязным телом

Совместимость правого тензора деформации Коши – Грина.

Необходимые и достаточные условия существования совместимого поля над односвязным телом

Мы можем показать, что это смешанные компоненты Тензор кривизны Римана – Кристоффеля. Следовательно, необходимые условия для -совместимость заключаются в том, что кривизна деформации Римана – Кристоффеля равна нулю.

Совместимость левого тензора деформации Коши – Грина.

Общие условия достаточности для левого тензора деформации Коши – Грина в трехмерном пространстве не известны. Условия совместимости для двумерных поля были обнаружены Джанет Блюм.[17][18]

Смотрите также

Рекомендации

  1. ^ а б Люблинер, Якоб (2008). Теория пластичности (PDF) (Пересмотренная ред.). Dover Publications. ISBN  978-0-486-46290-5. Архивировано из оригинал (PDF) 31 марта 2010 г.
  2. ^ А. Явари, Дж. Э. Марсден и М. Ортис, О пространственных и материальных ковариантных законах баланса упругости, Journal of Mathematical Physics, 47, 2006, 042903; С. 1–53.
  3. ^ Оуэнс, Эдуардо де Соуза Нето, Джордже Перич, Дэвид (2008). Вычислительные методы пластичности: теория и приложения. Чичестер, Западный Суссекс, Великобритания: Wiley. п. 65. ISBN  978-0-470-69452-7.
  4. ^ В ИЮПАК рекомендует называть этот тензор тензором деформации Коши.
  5. ^ а б c d A. Kaye, R. F. T. Stepto, W. J. Work, J. V. Aleman (Испания), A. Ya. Малкин (1998). «Определение терминов, относящихся к не предельным механическим свойствам полимеров». Pure Appl. Chem. 70 (3): 701–754. Дои:10.1351 / pac199870030701.CS1 maint: несколько имен: список авторов (связь)
  6. ^ Эдуардо Н. Дворкин, Марсела Б. Гольдшмит, 2006 г. Нелинейная континуа, п. 25, Springer ISBN  3-540-24985-0.
  7. ^ В ИЮПАК рекомендует называть этот тензор тензором деформации Грина.
  8. ^ Йирасек, Милан; Бажант, З. П. (2002) Неупругий расчет конструкций, Wiley, стр. 463 ISBN  0-471-98716-6
  9. ^ Дж. Н. Редди, Дэвид К. Гартлинг (2000) Метод конечных элементов в теплообмене и гидродинамике, п. 317, CRC Press ISBN  1-4200-8598-0.
  10. ^ Белычко, Тед; Лю, Винг Кам; Моран, Брайан (2000). Нелинейные конечные элементы для сплошных сред и структур. (перепечатка с исправлениями, 2006 г.). John Wiley & Sons Ltd., стр. 92–94. ISBN  978-0-471-98773-4.
  11. ^ Сет, Б. Р. (1961), «Обобщенная мера деформации с приложениями к физическим задачам», Краткий технический отчет MRC № 248, Центр математических исследований, Армия США, Университет Висконсина: 1–18
  12. ^ Сет, Б. Р. (1962), "Обобщенная мера деформации с приложениями к физическим проблемам", Симпозиум IUTAM по эффектам второго порядка в упругости, пластичности и механике жидкости, Хайфа, 1962 г.
  13. ^ Хилл Р. (1968), «О материальных неравенствах для простых материалов - I», Журнал механики и физики твердого тела, 16 (4): 229–242, Bibcode:1968JMPSo..16..229H, Дои:10.1016/0022-5096(68)90031-8
  14. ^ T.C. Дойл и Дж.Л. Эриксен (1956). «Нелинейная эластичность». Успехи прикладной механики 4, 53–115.
  15. ^ З.П. Бажант и Л. Седолин (1991). Устойчивость конструкций. Теории упругости, неупругости, разрушения и повреждения. Oxford Univ. Press, New York (2-е изд. Dover Publ., New York 2003; 3-е изд., World Scientific 2010).
  16. ^ З.П. Бажант (1998). "Легко вычисляемые тензоры с симметричной обратной аппроксимацией конечной деформации Генки и ее скорости." Журнал материалов технологии ASME, 120 (апрель), 131–136.
  17. ^ Блюм, Дж. А. (1989). «Условия совместимости для левого поля деформаций Коши – Грина». Журнал эластичности. 21 (3): 271–308. Дои:10.1007 / BF00045780. S2CID  54889553.
  18. ^ Ачарья, А. (1999). "Об условиях совместимости левого деформационного поля Коши – Грина в трех измерениях" (PDF). Журнал эластичности. 56 (2): 95–105. Дои:10.1023 / А: 1007653400249. S2CID  116767781.

дальнейшее чтение

  • Макоско, К. В. (1994). Реология: принципы, измерения и приложения. Издатели ВЧ. ISBN  1-56081-579-5.

внешняя ссылка