Поперечная изотропия

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

Поперечно изотропный материал является один с физическими свойствами , которые являются симметричными относительно оси, перпендикулярные к плоскости изотропии . Эта поперечная плоскость имеет бесконечное количество плоскостей симметрии, и поэтому в этой плоскости свойства материала одинаковы во всех направлениях. Следовательно, такие материалы также известны как «полярные анизотропные» материалы. В геофизике вертикально-поперечная изотропия (VTI) также известна как радиальная анизотропия.

Этот тип материала демонстрирует гексагональную симметрию (хотя технически это перестает быть верным для тензоров ранга 6 и выше), поэтому количество независимых констант в тензоре упругости (четвертого ранга) уменьшено до 5 (с 21 независимого константы в случае полностью анизотропного твердого тела ). Тензоры (второго ранга) удельного электрического сопротивления, магнитной проницаемости и т. Д. Имеют две независимые константы.

Пример трансверсально изотропных материалов [ править ]

Поперечно-изотропный эластичный материал.

Примером поперечно изотропного материала является так называемая осевая однонаправленная волокнистая композитная пластина, в которой волокна имеют круглое поперечное сечение. В однонаправленном композите плоскость, перпендикулярная направлению волокна, может рассматриваться как изотропная плоскость при больших длинах волн (низких частотах) возбуждения. На рисунке справа волокна будут выровнены по оси, перпендикулярной плоскости изотропии.${\ displaystyle x_ {2}}$

С точки зрения эффективных свойств геологические слои горных пород часто интерпретируются как поперечно-изотропные. Расчет эффективных упругих свойств таких слоев в петрологии был придуман апскейлингом Бэкуса , который описан ниже.

Матрица симметрии материала [ править ]

Матрица материала обладает симметрией относительно данного ортогонального преобразования ( ), если она не изменяется при выполнении этого преобразования. Для неизменности свойств материала при таком преобразовании нам потребуется${\ displaystyle {\ underline {\ underline {\ boldsymbol {K}}}}}$${\ displaystyle {\ boldsymbol {A}}}$

${\ displaystyle {\ boldsymbol {A}} \ cdot \ mathbf {f} = {\ boldsymbol {K}} \ cdot ({\ boldsymbol {A}} \ cdot {\ boldsymbol {d}}) \ подразумевает \ mathbf { f} = ({\ boldsymbol {A}} ^ {- 1} \ cdot {\ boldsymbol {K}} \ cdot {\ boldsymbol {A}}) \ cdot {\ boldsymbol {d}}}$

Следовательно, условие симметрии материала (используя определение ортогонального преобразования)

${\ displaystyle {\ boldsymbol {K}} = {\ boldsymbol {A}} ^ {- 1} \ cdot {\ boldsymbol {K}} \ cdot {\ boldsymbol {A}} = {\ boldsymbol {A}} ^ {T} \ cdot {\ boldsymbol {K}} \ cdot {\ boldsymbol {A}}}$

Ортогональные преобразования могут быть представлены в декартовых координатах матрицей, заданной как${\displaystyle 3\times 3}$${\displaystyle {\underline {\underline {\boldsymbol {A}}}}}$

${\displaystyle {\underline {\underline {\boldsymbol {A}}}}={\begin{bmatrix}A_{11}&A_{12}&A_{13}\\A_{21}&A_{22}&A_{23}\\A_{31}&A_{32}&A_{33}\end{bmatrix}}~.}$

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

${\displaystyle {\underline {\underline {\boldsymbol {K}}}}={\underline {\underline {{\boldsymbol {A}}^{T}}}}~{\underline {\underline {\boldsymbol {K}}}}~{\underline {\underline {\boldsymbol {A}}}}}$

Для трансверсально изотропного материала матрица имеет вид${\displaystyle {\underline {\underline {\boldsymbol {A}}}}}$

${\displaystyle {\underline {\underline {\boldsymbol {A}}}}={\begin{bmatrix}\cos \theta &\sin \theta &0\\-\sin \theta &\cos \theta &0\\0&0&1\end{bmatrix}}~.}$

где ось - ось симметрии . Матрица материала остается неизменной при повороте на любой угол вокруг оси.${\displaystyle x_{3}}$${\displaystyle \theta }$${\displaystyle x_{3}}$

В физике [ править ]

Материальные линейные соотношения в физике можно выразить в виде

${\displaystyle \mathbf {f} ={\boldsymbol {K}}\cdot \mathbf {d} }$

где - два вектора, представляющие физические величины, и - материальный тензор второго порядка. В матричной форме${\displaystyle \mathbf {d} ,\mathbf {f} }$${\displaystyle {\boldsymbol {K}}}$

${\displaystyle {\underline {\underline {\mathbf {f} }}}={\underline {\underline {\boldsymbol {K}}}}~{\underline {\underline {\mathbf {d} }}}\implies {\begin{bmatrix}f_{1}\\f_{2}\\f_{3}\end{bmatrix}}={\begin{bmatrix}K_{11}&K_{12}&K_{13}\\K_{21}&K_{22}&K_{23}\\K_{31}&K_{32}&K_{33}\end{bmatrix}}{\begin{bmatrix}d_{1}\\d_{2}\\d_{3}\end{bmatrix}}}$

Примеры физических проблем, соответствующих вышеуказанному шаблону, перечислены в таблице ниже. ^[1]

Проблема	${\displaystyle \mathbf {f} }$	${\displaystyle \mathbf {d} }$	${\displaystyle {\boldsymbol {K}}}$
Электрическая проводимость	Электрический ток ${\displaystyle \mathbf {J} }$	Электрическое поле ${\displaystyle \mathbf {E} }$	Электрическая проводимость ${\displaystyle {\boldsymbol {\sigma }}}$
Диэлектрики	Электрическое смещение ${\displaystyle \mathbf {D} }$	Электрическое поле ${\displaystyle \mathbf {E} }$	Электрическая проницаемость ${\displaystyle {\boldsymbol {\varepsilon }}}$
Магнетизм	Магнитная индукция ${\displaystyle \mathbf {B} }$	Магнитное поле ${\displaystyle \mathbf {H} }$	Магнитная проницаемость ${\displaystyle {\boldsymbol {\mu }}}$
Теплопроводность	Поток горячего воздуха ${\displaystyle \mathbf {q} }$	Температурный градиент ${\displaystyle -{\boldsymbol {\nabla }}T}$	Теплопроводность ${\displaystyle {\boldsymbol {\kappa }}}$
Диффузия	Поток частиц ${\displaystyle \mathbf {J} }$	Градиент концентрации ${\displaystyle -{\boldsymbol {\nabla }}c}$	Диффузность ${\displaystyle {\boldsymbol {D}}}$
Течение в пористой среде	Взвешенная скорость жидкости ${\displaystyle \eta _{\mu }\mathbf {v} }$	Градиент давления ${\displaystyle {\boldsymbol {\nabla }}P}$	Проницаемость для жидкости ${\displaystyle {\boldsymbol {\kappa }}}$
Эластичность	Стресс ${\displaystyle {\boldsymbol {\sigma }}}$	Штамм ${\displaystyle {\boldsymbol {\varepsilon }}}$	Жесткость ${\displaystyle \mathbf {C} }$

Использование в матрице означает, что . Использование ведет к и . Энергетические ограничения обычно требуют, а значит, и должны . Следовательно, свойства материала трансверсально изотропного материала описываются матрицей${\displaystyle \theta =\pi }$${\displaystyle {\underline {\underline {\boldsymbol {A}}}}}$${\displaystyle K_{13}=K_{31}=K_{23}=K_{32}=0}$${\displaystyle \theta ={\tfrac {\pi }{2}}}$${\displaystyle K_{11}=K_{22}}$${\displaystyle K_{12}=-K_{21}}$${\displaystyle K_{12},K_{21}\geq 0}$${\displaystyle K_{12}=K_{21}=0}$

${\displaystyle {\underline {\underline {\boldsymbol {K}}}}={\begin{bmatrix}K_{11}&0&0\\0&K_{11}&0\\0&0&K_{33}\end{bmatrix}}}$

В линейной эластичности [ править ]

Условие симметрии материала [ править ]

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

${\displaystyle {\underline {\underline {\boldsymbol {\sigma }}}}={\underline {\underline {\mathsf {C}}}}~{\underline {\underline {\boldsymbol {\varepsilon }}}}}$

или, используя обозначение Фойгта ,

${\displaystyle {\begin{bmatrix}\sigma _{1}\\\sigma _{2}\\\sigma _{3}\\\sigma _{4}\\\sigma _{5}\\\sigma _{6}\end{bmatrix}}={\begin{bmatrix}C_{11}&C_{12}&C_{13}&C_{14}&C_{15}&C_{16}\\C_{12}&C_{22}&C_{23}&C_{24}&C_{25}&C_{26}\\C_{13}&C_{23}&C_{33}&C_{34}&C_{35}&C_{36}\\C_{14}&C_{24}&C_{34}&C_{44}&C_{45}&C_{46}\\C_{15}&C_{25}&C_{35}&C_{45}&C_{55}&C_{56}\\C_{16}&C_{26}&C_{36}&C_{46}&C_{56}&C_{66}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{1}\\\varepsilon _{2}\\\varepsilon _{3}\\\varepsilon _{4}\\\varepsilon _{5}\\\varepsilon _{6}\end{bmatrix}}}$

Условием симметрии материала в линейно-упругих материалах является. ^[2]

${\displaystyle {\underline {\underline {\mathsf {C}}}}={\underline {\underline {{\mathsf {A}}_{\varepsilon }}}}^{T}~{\underline {\underline {\mathsf {C}}}}~{\underline {\underline {{\mathsf {A}}_{\varepsilon }}}}}$

где

${\displaystyle {\underline {\underline {{\mathsf {A}}_{\varepsilon }}}}={\begin{bmatrix}A_{11}^{2}&A_{12}^{2}&A_{13}^{2}&A_{12}A_{13}&A_{11}A_{13}&A_{11}A_{12}\\A_{21}^{2}&A_{22}^{2}&A_{23}^{2}&A_{22}A_{23}&A_{21}A_{23}&A_{21}A_{22}\\A_{31}^{2}&A_{32}^{2}&A_{33}^{2}&A_{32}A_{33}&A_{31}A_{33}&A_{31}A_{32}\\2A_{21}A_{31}&2A_{22}A_{32}&2A_{23}A_{33}&A_{22}A_{33}+A_{23}A_{32}&A_{21}A_{33}+A_{23}A_{31}&A_{21}A_{32}+A_{22}A_{31}\\2A_{11}A_{31}&2A_{12}A_{32}&2A_{13}A_{33}&A_{12}A_{33}+A_{13}A_{32}&A_{11}A_{33}+A_{13}A_{31}&A_{11}A_{32}+A_{12}A_{31}\\2A_{11}A_{21}&2A_{12}A_{22}&2A_{13}A_{23}&A_{12}A_{23}+A_{13}A_{22}&A_{11}A_{23}+A_{13}A_{21}&A_{11}A_{22}+A_{12}A_{21}\end{bmatrix}}}$

Тензор упругости [ править ]

Использование конкретных значений в матрице , ^[3] , можно показать , что четвертый ранг тензор упругости жесткости может быть записан в 2-индексе Фойхта обозначения в качестве матрицы${\displaystyle \theta }$${\displaystyle {\underline {\underline {\boldsymbol {A}}}}}$

${\displaystyle {\underline {\underline {\mathsf {C}}}}={\begin{bmatrix}C_{11}&C_{12}&C_{13}&0&0&0\\C_{12}&C_{11}&C_{13}&0&0&0\\C_{13}&C_{13}&C_{33}&0&0&0\\0&0&0&C_{44}&0&0\\0&0&0&0&C_{44}&0\\0&0&0&0&0&(C_{11}-C_{12})/2\end{bmatrix}}={\begin{bmatrix}C_{11}&C_{11}-2C_{66}&C_{13}&0&0&0\\C_{11}-2C_{66}&C_{11}&C_{13}&0&0&0\\C_{13}&C_{13}&C_{33}&0&0&0\\0&0&0&C_{44}&0&0\\0&0&0&0&C_{44}&0\\0&0&0&0&0&C_{66}\end{bmatrix}}.}$

Матрица упругости и жесткости имеет 5 независимых констант, которые связаны с хорошо известными инженерными модулями упругости следующим образом. Эти инженерные модули определены экспериментально.${\displaystyle C_{ij}}$

Матрица податливости (обратная матрице упругой жесткости) имеет вид

${\displaystyle {\underline {\underline {\mathsf {C}}}}^{-1}={\frac {1}{\Delta }}{\begin{bmatrix}C_{11}C_{33}-C_{13}^{2}&C_{13}^{2}-C_{12}C_{33}&(C_{12}-C_{11})C_{13}&0&0&0\\C_{13}^{2}-C_{12}C_{33}&C_{11}C_{33}-C_{13}^{2}&(C_{12}-C_{11})C_{13}&0&0&0\\(C_{12}-C_{11})C_{13}&(C_{12}-C_{11})C_{13}&C_{11}^{2}-C_{12}^{2}&0&0&0\\0&0&0&{\frac {\Delta }{C_{44}}}&0&0\\0&0&0&0&{\frac {\Delta }{C_{44}}}&0\\0&0&0&0&0&{\frac {2\Delta }{(C_{11}-C_{12})}}\end{bmatrix}}}$

где . В инженерных обозначениях ${\displaystyle \Delta :=(C_{11}-C_{12})[(C_{11}+C_{12})C_{33}-2C_{13}C_{13}]}$

${\displaystyle {\underline {\underline {\mathsf {C}}}}^{-1}={\begin{bmatrix}{\tfrac {1}{E_{\rm {x}}}}&-{\tfrac {\nu _{\rm {yx}}}{E_{\rm {x}}}}&-{\tfrac {\nu _{\rm {zx}}}{E_{\rm {z}}}}&0&0&0\\-{\tfrac {\nu _{\rm {xy}}}{E_{\rm {x}}}}&{\tfrac {1}{E_{\rm {x}}}}&-{\tfrac {\nu _{\rm {zx}}}{E_{\rm {z}}}}&0&0&0\\-{\tfrac {\nu _{\rm {xz}}}{E_{\rm {x}}}}&-{\tfrac {\nu _{\rm {xz}}}{E_{\rm {x}}}}&{\tfrac {1}{E_{\rm {z}}}}&0&0&0\\0&0&0&{\tfrac {1}{G_{\rm {xz}}}}&0&0\\0&0&0&0&{\tfrac {1}{G_{\rm {xz}}}}&0\\0&0&0&0&0&{\tfrac {2(1+\nu _{\rm {xy}})}{E_{\rm {x}}}}\end{bmatrix}}}$

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

${\displaystyle E_{L}=E_{\rm {z}}=C_{33}-2C_{13}C_{13}/(C_{11}+C_{12})}$

Аналогично, поперечный модуль Юнга равен

${\displaystyle E_{T}=E_{\rm {x}}=E_{\rm {y}}=(C_{11}-C_{12})(C_{11}C_{33}+C_{12}C_{33}-2C_{13}C_{13})/(C_{11}C_{33}-C_{13}C_{13})}$

Inplane модуль сдвига является

${\displaystyle G_{xy}=(C_{11}-C_{12})/2=C_{66}}$

а коэффициент Пуассона для нагрузки вдоль полярной оси равен

${\displaystyle \nu _{LT}=\nu _{zx}=C_{13}/(C_{11}+C_{12})}$.

Здесь L представляет продольное (полярное) направление, а T представляет собой поперечное направление.

В геофизике [ править ]

В геофизике распространено предположение, что горные образования земной коры являются локально полярными анизотропными (трансверсально изотропными); это простейший случай, представляющий геофизический интерес. Апскейлинг Backus ^[4] часто используется для определения эффективных поперечно-изотропных упругих постоянных слоистых сред для длинноволновых сейсмических волн.

Предположения, сделанные в приближении Бэкуса, следующие:

• Все материалы линейно эластичны.
• Отсутствуют источники собственного рассеивания энергии (например, трение)
• Действительно в пределе бесконечной длины волны, следовательно, хорошие результаты только в том случае, если толщина слоя намного меньше длины волны.
• Статистика распределения упругих свойств слоев стационарна, т. Е. Нет коррелированного тренда по этим свойствам.

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

• квази- P-волна ( направление поляризации почти совпадает с направлением распространения)
• квази- S волна
• S-волна (поляризованная ортогональная квази-S-волне, оси симметрии и направлению распространения).

Решения задач распространения волн в таких средах могут быть построены из этих плоских волн с использованием синтеза Фурье .

Апскейлинг Backus (длинноволновое приближение) [ править ]

Слоистая модель однородного и изотропного материала может быть увеличена до поперечно-изотропной среды, предложенной Бакусом. ^[4]

Бэкус представил эквивалентную теорию среды, гетерогенную среду можно заменить однородной, которая предсказывает распространение волн в реальной среде. ^[5] Бэкус показал, что наслоение в масштабе, намного меньшем длины волны, оказывает влияние и что ряд изотропных слоев может быть заменен однородной трансверсально-изотропной средой, которая ведет себя точно так же, как реальная среда при статической нагрузке в бесконечный предел длины волны.

Если каждый слой описывается 5 трансверсально-изотропными параметрами , задавая матрицу${\displaystyle i}$${\displaystyle (a_{i},b_{i},c_{i},d_{i},e_{i})}$

${\displaystyle {\underline {\underline {{\mathsf {C}}_{i}}}}={\begin{bmatrix}a_{i}&a_{i}-2e_{i}&b_{i}&0&0&0\\a_{i}-2e_{i}&a_{i}&b_{i}&0&0&0\\b_{i}&b_{i}&c_{i}&0&0&0\\0&0&0&d_{i}&0&0\\0&0&0&0&d_{i}&0\\0&0&0&0&0&e_{i}\\\end{bmatrix}}}$

Модули упругости для эффективной среды будут

${\displaystyle {\underline {\underline {{\mathsf {C}}_{\mathrm {eff} }}}}={\begin{bmatrix}A&A-2E&B&0&0&0\\A-2E&A&B&0&0&0\\B&B&C&0&0&0\\0&0&0&D&0&0\\0&0&0&0&D&0\\0&0&0&0&0&E\end{bmatrix}}}$

где

${\displaystyle {\begin{aligned}A&=\langle a-b^{2}c^{-1}\rangle +\langle c^{-1}\rangle ^{-1}\langle bc^{-1}\rangle ^{2}\\B&=\langle c^{-1}\rangle ^{-1}\langle bc^{-1}\rangle \\C&=\langle c^{-1}\rangle ^{-1}\\D&=\langle d^{-1}\rangle ^{-1}\\E&=\langle e\rangle \\\end{aligned}}}$

${\displaystyle \langle \cdot \rangle }$ обозначает средневзвешенный объем по всем слоям.

Сюда входят изотропные слои, так как слой изотропен, если , и .${\displaystyle b_{i}=a_{i}-2e_{i}}$${\displaystyle a_{i}=c_{i}}$${\displaystyle d_{i}=e_{i}}$

Приближение коротких и средних длин волн [ править ]

Решения задач распространения волн в линейных упругих трансверсально изотропных средах могут быть построены путем наложения решений для квази-P-волны, квази-S-волны и S-волны, поляризованной ортогональной квази-S-волне. Однако уравнения для углового изменения скорости являются алгебраически сложными, и скорости плоских волн являются функциями угла распространения . ^[6] Зависящие от направления волновые скорости для упругих волн, проходящих через материал, могут быть найдены с помощью уравнения Кристоффеля и определены как ^[7]${\displaystyle \theta }$

${\displaystyle {\begin{aligned}V_{qP}(\theta )&={\sqrt {\frac {C_{11}\sin ^{2}(\theta )+C_{33}\cos ^{2}(\theta )+C_{44}+{\sqrt {M(\theta )}}}{2\rho }}}\\V_{qS}(\theta )&={\sqrt {\frac {C_{11}\sin ^{2}(\theta )+C_{33}\cos ^{2}(\theta )+C_{44}-{\sqrt {M(\theta )}}}{2\rho }}}\\V_{S}&={\sqrt {\frac {C_{66}\sin ^{2}(\theta )+C_{44}\cos ^{2}(\theta )}{\rho }}}\\M(\theta )&=\left[\left(C_{11}-C_{44}\right)\sin ^{2}(\theta )-\left(C_{33}-C_{44}\right)\cos ^{2}(\theta )\right]^{2}+\left(C_{13}+C_{44}\right)^{2}\sin ^{2}(2\theta )\\\end{aligned}}}$

где - угол между осью симметрии и направлением распространения волны, - массовая плотность, - элементы матрицы упругой жесткости . Параметры Томсена используются для упрощения этих выражений и облегчения их понимания.${\displaystyle {\begin{aligned}\theta \end{aligned}}}$${\displaystyle \rho }$${\displaystyle C_{ij}}$

Параметры Томсена [ править ]

Параметры Томсена ^[8] представляют собой безразмерные комбинации модулей упругости , характеризующие трансверсально-изотропные материалы, которые встречаются, например, в геофизике . В терминах компонентов матрицы упругой жесткости эти параметры определяются как:

${\displaystyle {\begin{aligned}\epsilon &={\frac {C_{11}-C_{33}}{2C_{33}}}\\\delta &={\frac {(C_{13}+C_{44})^{2}-(C_{33}-C_{44})^{2}}{2C_{33}(C_{33}-C_{44})}}\\\gamma &={\frac {C_{66}-C_{44}}{2C_{44}}}\end{aligned}}}$

где индекс 3 указывает ось симметрии ( ). Эти параметры вместе с соответствующими скоростями продольных и поперечных волн можно использовать для характеристики распространения волн через слабоанизотропные слоистые среды. Эмпирически параметры Томсена для большинства слоистых горных пород намного ниже 1.${\displaystyle \mathbf {e} _{3}}$

Это имя относится к Леону Томсену, профессору геофизики из Хьюстонского университета , который предложил эти параметры в своей статье 1986 года «Слабая упругая анизотропия».

Упрощенные выражения для волновых скоростей [ править ]

В геофизике анизотропия упругих свойств обычно мала, и в этом случае . Когда точные выражения для скоростей волн, приведенные выше, линеаризуются в этих малых величинах, они упрощаются до${\displaystyle \delta ,\gamma ,\epsilon \ll 1}$

${\displaystyle {\begin{aligned}V_{qP}(\theta )&\approx V_{P0}(1+\delta \sin ^{2}\theta \cos ^{2}\theta +\epsilon \sin ^{4}\theta )\\V_{qS}(\theta )&\approx V_{S0}\left[1+\left({\frac {V_{P0}}{V_{S0}}}\right)^{2}(\epsilon -\delta )\sin ^{2}\theta \cos ^{2}\theta \right]\\V_{S}(\theta )&\approx V_{S0}(1+\gamma \sin ^{2}\theta )\end{aligned}}}$

где

${\displaystyle V_{P0}={\sqrt {C_{33}/\rho }}~;~~V_{S0}={\sqrt {C_{44}/\rho }}}$

- скорости продольных и поперечных волн в направлении оси симметрии ( ) (в геофизике это обычно, но не всегда, вертикальное направление). Обратите внимание, что можно дополнительно линеаризовать, но это не приводит к дальнейшему упрощению.${\displaystyle \mathbf {e} _{3}}$${\displaystyle \delta }$

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

См. Также [ править ]

• Закон Гука
• Линейная эластичность
• Ортотропный материал

Ссылки [ править ]