Тензор упругости

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

Определение и общий вид

В линейной упругой среде связь между компонентами тензора напряжений \(\sigma_{ij}\) и тензора деформаций \(\varepsilon_{kl}\) записывается в виде обобщённого закона Гука:

\[ \sigma_{ij} = C_{ijkl} \varepsilon_{kl} \]

где \(C_{ijkl}\) — компоненты тензора упругости. В трёхмерном пространстве тензор четвёртого ранга имеет \(3^4 = 81\) компоненту. Однако благодаря симметриям тензоров напряжений и деформаций, а также фундаментальным термодинамическим принципам, число независимых компонент значительно сокращается.

Симметрии тензора упругости

1. Симметрия по индексам напряжений: \(\sigma_{ij} = \sigma_{ji}\), следовательно \(C_{ijkl} = C_{jikl}\).
2. Симметрия по индексам деформаций: \(\varepsilon_{kl} = \varepsilon_{lk}\), следовательно \(C_{ijkl} = C_{ijlk}\).
3. Термодинамическая симметрия (симметрия главных значений): Из существования упругого потенциала следует, что \(C_{ijkl} = C_{klij}\).

Эти симметрии уменьшают количество независимых компонент до 21. Такой тензор описывает анизотропное тело общего вида (триклинную сингонию).

Матричная запись (нотация Фойгта)

Для удобства вычислений тензор упругости четвёртого ранга часто представляют в виде матрицы \(6 \times 6\) (матрицы жёсткости) с использованием нотации Фойгта. В этой нотации пары индексов заменяются одним индексом по правилу:

Тогда закон Гука записывается как:

\[ \sigma_i = C_{ij} \varepsilon_j \]

где \(i, j = 1, \dots, 6\), а \(C_{ij}\) — симметричная матрица \(6 \times 6\), содержащая 21 независимый элемент.

Влияние кристаллической симметрии

Симметрия кристаллической решётки материала накладывает дополнительные ограничения на вид тензора упругости. Чем выше симметрия, тем меньше независимых констант.

Изотропное тело

Для изотропных материалов (стекло, поликристаллы без текстуры, аморфные тела) упругие свойства не зависят от направления. В этом случае тензор упругости описывается всего двумя независимыми константами, которые обычно выбирают как модуль Юнга \(E\) и коэффициент Пуассона \(\nu\), или как модуль сдвига \(G\) и модуль объёмного сжатия \(K\).

В нотации Фойгта для изотропного тела матрица жёсткости имеет вид:

\[ C = \begin{pmatrix} \lambda + 2\mu & \lambda & \lambda & 0 & 0 & 0 \\ \lambda & \lambda + 2\mu & \lambda & 0 & 0 & 0 \\ \lambda & \lambda & \lambda + 2\mu & 0 & 0 & 0 \\ 0 & 0 & 0 & \mu & 0 & 0 \\ 0 & 0 & 0 & 0 & \mu & 0 \\ 0 & 0 & 0 & 0 & 0 & \mu \end{pmatrix} \]

где \(\lambda\) и \(\mu\) — параметры Ламе. Связь между ними и инженерными константами:

• \(\mu = G = \frac{E}{2(1+\nu)}\)
• \(\lambda = \frac{\nu E}{(1+\nu)(1-2\nu)}\)

Кристаллы кубической сингонии

Кристаллы с кубической решёткой (например, медь, алюминий, кремний) имеют три независимые константы: \(C_{11}\), \(C_{12}\) и \(C_{44}\). Матрица жёсткости:

\[ C = \begin{pmatrix} C_{11} & C_{12} & C_{12} & 0 & 0 & 0 \\ C_{12} & C_{11} & C_{12} & 0 & 0 & 0 \\ C_{12} & C_{12} & C_{11} & 0 & 0 & 0 \\ 0 & 0 & 0 & C_{44} & 0 & 0 \\ 0 & 0 & 0 & 0 & C_{44} & 0 \\ 0 & 0 & 0 & 0 & 0 & C_{44} \end{pmatrix} \]

Другие сингонии

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

• Гексагональная (трансверсально-изотропная): 5 констант.
• Тетрагональная: 6 или 7 констант.
• Ромбическая: 9 констант.
• Моноклинная: 13 констант.
• Триклинная: 21 константа.

Тензор податливости

Наряду с тензором упругости \(C_{ijkl}\) существует обратный ему тензор — тензор податливости \(S_{ijkl}\), который связывает деформации с напряжениями:

\[ \varepsilon_{ij} = S_{ijkl} \sigma_{kl} \]

Матрица податливости \(S\) является обратной к матрице жёсткости \(C\) в нотации Фойгта. Компоненты тензора податливости имеют размерность, обратную размерности напряжения (Па⁻¹).

Физический смысл компонент

Компоненты тензора упругости имеют физический смысл. Например, в изотропном материале:

• \(C_{11} = \lambda + 2\mu\) — это напряжение, необходимое для создания единичной продольной деформации при отсутствии поперечной деформации.
• \(C_{44} = \mu\) — это напряжение сдвига, необходимое для создания единичной деформации сдвига.

В кристаллах компоненты \(C_{ij}\) напрямую связаны со скоростью распространения упругих волн (акустических фононов) в различных кристаллографических направлениях. Например, скорость продольной волны в направлении [100] в кубическом кристалле равна \(\sqrt{C_{11}/\rho}\), где \(\rho\) — плотность.

Применение

Тензор упругости является ключевым понятием в:

1. Механике сплошных сред и теории упругости: Используется для расчёта напряжённо-деформированного состояния конструкций и деталей.
2. Физике твёрдого тела: Необходим для описания упругих свойств кристаллов, фазовых переходов, акустических свойств.
3. Сейсмологии и геофизике: Позволяет моделировать распространение сейсмических волн в земной коре, определять состав и структуру горных пород.
4. Материаловедении: Используется для характеризации композиционных материалов, полимеров, биологических тканей.
5. Микроэлектронике: При расчёте механических напряжений в полупроводниковых структурах и микроэлектромеханических системах (МЭМС).

Методы определения

Экспериментальное определение компонент тензора упругости обычно проводится с помощью:

• Ультразвуковых методов: Измерение скоростей продольных и поперечных акустических волн в образцах различной ориентации.
• Резонансных методов: Определение собственных частот колебаний образцов (например, метод резонансной ультразвуковой спектроскопии).
• Статических методов: Испытания на растяжение, сжатие и сдвиг образцов, вырезанных вдоль различных кристаллографических направлений.
• Методов рассеяния света (Бриллюэновское рассеяние): Для изучения упругих свойств на микро- и наноуровне.
• Методов компьютерного моделирования: Атомистическое моделирование (метод молекулярной динамики) и расчёты из первых принципов (ab initio) позволяют предсказывать упругие константы для новых материалов.

Ограничения модели

Тензор упругости описывает только линейное упругое поведение материала. Он не применим для:

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

Источники

1. Ландау Л. Д., Лифшиц Е. М. Теоретическая физика. Том VII. Теория упругости. — М.: Наука, 1987.
2. Новацкий В. Теория упругости. — М.: Мир, 1975.
3. Сиротин Ю. И., Шаскольская М. П. Основы кристаллофизики. — М.: Наука, 1979.
4. Тимошенко С. П., Гудьер Дж. Теория упругости. — М.: Наука, 1975.
5. Nye J. F. Physical Properties of Crystals: Their Representation by Tensors and Matrices. — Oxford University Press, 1985.