Material ortótropo

La madera es un ejemplo de un material ortótropo: sus propiedades en tres direcciones perpendiculares (axial, radial y circunferencial) son diferentes.

Un material ortótropo u ortótropico es una material cuya microestructura local tiene dos o tres ejes ortogonales entre sí, de doble simetría rotacional, de forma que sus propiedades físicas son, en general, diferentes en las direcciones de cada uno de esos ejes. Por esta razón, los materiales ortótropos son anisótropos, ya que sus propiedades físicas dependerán de la dirección en que sean medidas. En cambio, un material isótropo exhibirá las mismas propiedades físicas en cualquier dirección. Los materiales transversalmente isótropos son un subtipo de material ortótropo.

Ortotropía general en física

Relaciones constitutivas para materiales anisótropicos

El comportamiento de los materiales se presenta en teorías físicas como relaciones constitutivas. Una gran cantidad de comportamientos físicos se pueden representar con modelos lineales de materiales que toman la forma de un tensor de segundo orden. El tensor material relaciona dos vectores y puede ser escrito como:

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

donde d , f {\displaystyle \mathbf {d} ,\mathbf {f} } son dos vectores que representan magnitudes físicas y K {\displaystyle {\boldsymbol {K}}} es el tensor material de segundo orden. Si se expresa la ecuación anterior en términos de componentes respecto de un sistema de coordenadas ortogonales resulta:

f i = K i j   d j   . {\displaystyle f_{i}=K_{ij}~d_{j}~.}

En la relación anterior se ha asumido la sumatoria sobre índices repetidos (notación de Einstein). En forma matricial se obtiene:

f _ = K _ _   d _ [ f 1 f 2 f 3 ] = [ K 11 K 12 K 13 K 21 K 22 K 23 K 31 K 32 K 33 ] [ d 1 d 2 d 3 ] {\displaystyle {\underline {\mathbf {f} }}={\underline {\underline {\boldsymbol {K}}}}~{\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}}}

Ejemplos de problemas físicos que encajan en esta descripción están listados en la lista a continuación:[1]

Problema f {\displaystyle \mathbf {f} } d {\displaystyle \mathbf {d} } K {\displaystyle {\boldsymbol {K}}}
Conducción eléctrica Corriente eléctrica
J {\displaystyle \mathbf {J} }
Campo eléctrico
E {\displaystyle \mathbf {E} }
Conductividad eléctrica
σ {\displaystyle {\boldsymbol {\sigma }}}
Dieléctricos Densidad de flujo eléctrico
D {\displaystyle \mathbf {D} }
Campo eléctrico
E {\displaystyle \mathbf {E} }
Permitividad
ε {\displaystyle {\boldsymbol {\varepsilon }}}
Magnetismo Inducción magnética
B {\displaystyle \mathbf {B} }
Excitación magnética
H {\displaystyle \mathbf {H} }
Permeabilidad magnética
μ {\displaystyle {\boldsymbol {\mu }}}
Conducción de calor Flujo de calor
q {\displaystyle \mathbf {q} }
Gradiente térmico
T {\displaystyle -{\boldsymbol {\nabla }}T}
Conductividad térmica
κ {\displaystyle {\boldsymbol {\kappa }}}
Difusión Flujo de partículas
J {\displaystyle \mathbf {J} }
Gradiente de concentración
c {\displaystyle -{\boldsymbol {\nabla }}c}
Coeficiente de difusión
D {\displaystyle {\boldsymbol {D}}}
Mecánica de sólidos tensión mecánica
σ {\displaystyle {\boldsymbol {\sigma }}}
Deformación
ε {\displaystyle {\boldsymbol {\varepsilon }}}
Relación Constitutiva
K {\displaystyle {\boldsymbol {K}}}
Flujo en medios porosos Velocidad ponderada
η μ v {\displaystyle \eta _{\mu }\mathbf {v} }
Gradiente barométrico
P {\displaystyle {\boldsymbol {\nabla }}P}
Permeabilidad
κ {\displaystyle {\boldsymbol {\kappa }}}

Condición para simetría del material

La matriz del material K _ _ {\displaystyle {\underline {\underline {\boldsymbol {K}}}}} tiene simetría respecto de una transformación ortogonal ( A {\displaystyle {\boldsymbol {A}}} ) si no cambia cuando se la somete a una transformación de ese tipo. Para la invarianza de las propiedades del material cuando se le aplica una transformación de tales características se requiere que:

A f = K ( A d ) f = ( A 1 K A ) d {\displaystyle {\boldsymbol {A}}\cdot \mathbf {f} ={\boldsymbol {K}}\cdot ({\boldsymbol {A}}\cdot {\boldsymbol {d}})\implies \mathbf {f} =({\boldsymbol {A}}^{-1}\cdot {\boldsymbol {K}}\cdot {\boldsymbol {A}})\cdot {\boldsymbol {d}}}

Por lo tanto, la condición de simetría del material es (usando la definición de una transformación ortogonal)

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

Las transformaciones ortogonales se pueden representar en coordenadas cartesianas mediante una matriz A _ _ {\displaystyle {\underline {\underline {\boldsymbol {A}}}}} de 3 × 3 {\displaystyle 3\times 3} dada por

A _ _ = [ A 11 A 12 A 13 A 21 A 22 A 23 A 31 A 32 A 33 ]   . {\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}}~.}

Por lo tanto, la condición de simetría se puede escribir en forma matricial como

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

Propiedades materiales del ortótropo

Un material ortótropico tiene tres planos de simetría ortogonal. Si tomamos un sistema coordinado ortonormal de manera que los ejes coincidan con las normales por los tres planos de simetría, las matrices de transformación son:

A 1 _ _ = [ 1 0 0 0 1 0 0 0 1 ]   ;     A 2 _ _ = [ 1 0 0 0 1 0 0 0 1 ]   ;     A 3 _ _ = [ 1 0 0 0 1 0 0 0 1 ] {\displaystyle {\underline {\underline {{\boldsymbol {A}}_{1}}}}={\begin{bmatrix}-1&0&0\\0&1&0\\0&0&1\end{bmatrix}}~;~~{\underline {\underline {{\boldsymbol {A}}_{2}}}}={\begin{bmatrix}1&0&0\\0&-1&0\\0&0&1\end{bmatrix}}~;~~{\underline {\underline {{\boldsymbol {A}}_{3}}}}={\begin{bmatrix}1&0&0\\0&1&0\\0&0&-1\end{bmatrix}}}

Se puede ver que si la matriz K _ _ {\displaystyle {\underline {\underline {\boldsymbol {K}}}}} para un material es invariable ante la reflexión de dos planos ortogonales, entonces es también invariable ante la reflexión sobre un tercer plano ortogonal.

Teniendo en cuenta la reflexión A 3 _ _ {\displaystyle {\underline {\underline {{\boldsymbol {A}}_{3}}}}} sobre el plano 1 2 {\displaystyle 1-2\,} Entonces tenemos:

K _ _ = A 3 T _ _   K _ _   A 3 _ _ = [ K 11 K 12 K 13 K 21 K 22 K 23 K 31 K 32 K 33 ] {\displaystyle {\underline {\underline {\boldsymbol {K}}}}={\underline {\underline {{\boldsymbol {A}}_{3}^{T}}}}~{\underline {\underline {\boldsymbol {K}}}}~{\underline {\underline {{\boldsymbol {A}}_{3}}}}={\begin{bmatrix}K_{11}&K_{12}&-K_{13}\\K_{21}&K_{22}&-K_{23}\\-K_{31}&-K_{32}&K_{33}\end{bmatrix}}}

La relación citada anteriormente implica que K 13 = K 23 = K 31 = K 32 = 0 {\displaystyle K_{13}=K_{23}=K_{31}=K_{32}=0} . Tengamos ahora en cuenta una reflexión A 2 _ _ {\displaystyle {\underline {\underline {{\boldsymbol {A}}_{2}}}}} sobre el plano 1 3 {\displaystyle 1-3\,} Entonces tenemos:

K _ _ = A 2 T _ _   K _ _   A 2 _ _ = [ K 11 K 12 0 K 21 K 22 0 0 0 K 33 ] {\displaystyle {\underline {\underline {\boldsymbol {K}}}}={\underline {\underline {{\boldsymbol {A}}_{2}^{T}}}}~{\underline {\underline {\boldsymbol {K}}}}~{\underline {\underline {{\boldsymbol {A}}_{2}}}}={\begin{bmatrix}K_{11}&-K_{12}&0\\-K_{21}&K_{22}&0\\0&0&K_{33}\end{bmatrix}}}

Esto implica que K 12 = K 21 = 0 {\displaystyle K_{12}=K_{21}=0} . Por lo tanto, las propiedades materiales de un material ortótropico están descritas por la matriz:

K _ _ = [ K 11 0 0 0 K 22 0 0 0 K 33 ] {\displaystyle {\underline {\underline {\boldsymbol {K}}}}={\begin{bmatrix}K_{11}&0&0\\0&K_{22}&0\\0&0&K_{33}\end{bmatrix}}}

Ortotropía en la elasticidad lineal

Ejemplos en elasticidad

Un material elástico ortótropo (lineal y homogéneo) requiere en general especificar 9 constantes elásticas, mientras que un material transversalmente isótropo sólo requiere 5 y un material puramente isótropo únicamente requiere 2.

Un ejemplo común de un material ortótropico con dos ejes de simetría es un polímero reforzado por fibras de vidrio o carbono paralelas. La resistencia y rigidez de un material compuesto de estas características generalmente serán mayores en la dirección paralela a las fibras respecto de la dirección transversal. Otro ejemplo es una membrana biológica, en la cual las propiedades en el plano de la membrana son diferentes a aquellas en la dirección perpendicular. Estos materiales son llamados a veces transversalmente isótropicos. Un ejemplo familiar de un material ortótropico con tres ejes mutuamente perpendiculares es la madera, en la que las propiedades (tales como resistencia y rigidez) a lo largo de sus fibras y en cada una de dos direcciones perpendiculares son diferentes. La ecuación de Hankinson provee una forma de cuantificar la diferencia en resistencia entre las diferentes direcciones. Otro ejemplo es un metal que ha sido laminado para producir una lámina. En él, las propiedades en la dirección de laminado y en cada una de las dos direcciones transversales serán diferentes debido a la estructura anisotrópica que se produce durante el laminado.

Es importante recordar que un material que es anisotrópico en una escala de medida puede ser isotrópico en otra, en general, mayor. Por ejemplo, la mayoría de los metales son policristalinos con granos muy pequeños. Cada uno de los granos individuales puede ser anisotrópico, pero el material como globalmente está formado por granos de orientación al azar, por lo que sus propiedades mecánicas medidas serán un promedio de las propiedades en todas las posibles orientaciones de los granos individuales.

Elasticidad anisótropica

En la elasticidad lineal, la relación entre tensión y deformación dependen del tipo de material en consideración. Esta relación es conocida como Ley de Hooke. Para materiales anisótropicos lineales, la Ley de Hooke puede ser escrita como:[2]

σ = c : ε {\displaystyle {\boldsymbol {\sigma }}={\mathsf {c}}:{\boldsymbol {\varepsilon }}}

En donde σ {\displaystyle {\boldsymbol {\sigma }}} es el tensor de la tensión, ε {\displaystyle {\boldsymbol {\varepsilon }}} es el tensor del esfuerzo, y c {\displaystyle {\mathsf {c}}} es el tensor de fuerza elástico. Si los tensores de las expresiones anteriormente citadas están descritos en términos de componentes con respecto a un sistema ortonormal coordinado podemos escribir:

σ i j = c i j k   ε k {\displaystyle \sigma _{ij}=c_{ijk\ell }~\varepsilon _{k\ell }}

En donde la suma se ha asumido sobre índices repetidos (convenio de la suma de Einstein). Desde que la tensión y el esfuerzo son simétricos, y desde que la relación de tensión-esfuerzo en elasticidad lineal puede ser derivada desde una función de esfuerzo de densidad energética, Las siguientes simetrías se mantienen por materiales elásticos lineales.

c i j k = c j i k   ,     c i j k = c i j k   ,     c i j k = c k i j   . {\displaystyle c_{ijk\ell }=c_{jik\ell }~,~~c_{ijk\ell }=c_{ij\ell k}~,~~c_{ijk\ell }=c_{k\ell ij}~.}

Debido a las simetrías citadas anteriormente, la relación de tensión-esfuerzo en elasticidad lineal puede ser expresada en forma matriz como:

[ σ 11 σ 22 σ 33 σ 23 σ 31 σ 12 ] = [ c 1111 c 1122 c 1133 c 1123 c 1131 c 1112 c 2211 c 2222 c 2233 c 2223 c 2231 c 2212 c 3311 c 3322 c 3333 c 3323 c 3331 c 3312 c 2311 c 2322 c 2333 c 2323 c 2331 c 2312 c 3111 c 3122 c 3133 c 3123 c 3131 c 3112 c 1211 c 1222 c 1233 c 1223 c 1231 c 1212 ] [ ε 11 ε 22 ε 33 2 ε 23 2 ε 31 2 ε 12 ] {\displaystyle {\begin{bmatrix}\sigma _{11}\\\sigma _{22}\\\sigma _{33}\\\sigma _{23}\\\sigma _{31}\\\sigma _{12}\end{bmatrix}}={\begin{bmatrix}c_{1111}&c_{1122}&c_{1133}&c_{1123}&c_{1131}&c_{1112}\\c_{2211}&c_{2222}&c_{2233}&c_{2223}&c_{2231}&c_{2212}\\c_{3311}&c_{3322}&c_{3333}&c_{3323}&c_{3331}&c_{3312}\\c_{2311}&c_{2322}&c_{2333}&c_{2323}&c_{2331}&c_{2312}\\c_{3111}&c_{3122}&c_{3133}&c_{3123}&c_{3131}&c_{3112}\\c_{1211}&c_{1222}&c_{1233}&c_{1223}&c_{1231}&c_{1212}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{11}\\\varepsilon _{22}\\\varepsilon _{33}\\2\varepsilon _{23}\\2\varepsilon _{31}\\2\varepsilon _{12}\end{bmatrix}}}

Una representación alternativa en Notación Voigt es:

[ σ 1 σ 2 σ 3 σ 4 σ 5 σ 6 ] = [ 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 ] [ ε 1 ε 2 ε 3 ε 4 ε 5 ε 6 ] {\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}}}

o también

σ _ _ = C _ _   ε _ _ {\displaystyle {\underline {\underline {\boldsymbol {\sigma }}}}={\underline {\underline {\mathsf {C}}}}~{\underline {\underline {\boldsymbol {\varepsilon }}}}}

La matriz de fuerza C _ _ {\displaystyle {\underline {\underline {\mathsf {C}}}}} en la relación anterior satisface el punto de simetría.[3]

Condición para materiales simétricos

La matriz de fuerza C _ _ {\displaystyle {\underline {\underline {\mathsf {C}}}}} satisface una simetría de condición dada si no se cambia cuando está sometida a la transformación ortogonal correspondiente. La transformación ortogonal podría representar simetría con respecto a un punto, un eje, o un plano. Las transformaciones ortogonales en elasticidad lineal incluyen rotaciones y reflexiones, pero no cambios de transformación en la forma y pueden ser representadas en coordinadas ortonormales, dadas por la matriz 3 × 3 {\displaystyle 3\times 3} A _ _ {\displaystyle {\underline {\underline {\mathbf {A} }}}} dada por:

A _ _ = [ A 11 A 12 A 13 A 21 A 22 A 23 A 31 A 32 A 33 ]   . {\displaystyle {\underline {\underline {\mathbf {A} }}}={\begin{bmatrix}A_{11}&A_{12}&A_{13}\\A_{21}&A_{22}&A_{23}\\A_{31}&A_{32}&A_{33}\end{bmatrix}}~.}

En Notación Voigt, la transformación de la matriz para el tensor de fuerza puede ser expresada como una matriz: 6 × 6 {\displaystyle 6\times 6} dada por A σ _ _ {\displaystyle {\underline {\underline {{\mathsf {A}}_{\sigma }}}}} [3]

A σ _ _ = [ A 11 2 A 12 2 A 13 2 2 A 12 A 13 2 A 11 A 13 2 A 11 A 12 A 21 2 A 22 2 A 23 2 2 A 22 A 23 2 A 21 A 23 2 A 21 A 22 A 31 2 A 32 2 A 33 2 2 A 32 A 33 2 A 31 A 33 2 A 31 A 32 A 21 A 31 A 22 A 32 A 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 A 11 A 31 A 12 A 32 A 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 A 11 A 21 A 12 A 22 A 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 ] {\displaystyle {\underline {\underline {{\mathsf {A}}_{\sigma }}}}={\begin{bmatrix}A_{11}^{2}&A_{12}^{2}&A_{13}^{2}&2A_{12}A_{13}&2A_{11}A_{13}&2A_{11}A_{12}\\A_{21}^{2}&A_{22}^{2}&A_{23}^{2}&2A_{22}A_{23}&2A_{21}A_{23}&2A_{21}A_{22}\\A_{31}^{2}&A_{32}^{2}&A_{33}^{2}&2A_{32}A_{33}&2A_{31}A_{33}&2A_{31}A_{32}\\A_{21}A_{31}&A_{22}A_{32}&A_{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}\\A_{11}A_{31}&A_{12}A_{32}&A_{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}\\A_{11}A_{21}&A_{12}A_{22}&A_{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}}}

La transformación para el tensor de esfuerzo tiene una ligera forma diferente debido a la elección de notación. Esta transformación de matriz es:

A ε _ _ = [ 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 2 A 21 A 31 2 A 22 A 32 2 A 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 2 A 11 A 31 2 A 12 A 32 2 A 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 2 A 11 A 21 2 A 12 A 22 2 A 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 ] {\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}}}

Se puede mostrar que: A ε _ _ T = A σ _ _ 1 {\displaystyle {\underline {\underline {{\mathsf {A}}_{\varepsilon }}}}^{T}={\underline {\underline {{\mathsf {A}}_{\sigma }}}}^{-1}} .

Las propiedades elásticas de un continuo son invariables ante una transformación A _ _ {\displaystyle {\underline {\underline {\mathbf {A} }}}} si y sólo si:[3]

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

Matrices de rigidez y flexibilidad para materiales ortótropicos

Para un material elástico ortótropico existen tres planos de simetría ortogonales entre sí. Por esa razón, puede escogerse una base ortonormal, de tal manera que cada uno de los vectores que forman dicha base sea el vector normal a uno de los planos de simetría, las matrices de transformación serían:

A 1 _ _ = [ 1 0 0 0 1 0 0 0 1 ]   ;     A 2 _ _ = [ 1 0 0 0 1 0 0 0 1 ]   ;     A 3 _ _ = [ 1 0 0 0 1 0 0 0 1 ] {\displaystyle {\underline {\underline {\mathbf {A} _{1}}}}={\begin{bmatrix}-1&0&0\\0&1&0\\0&0&1\end{bmatrix}}~;~~{\underline {\underline {\mathbf {A} _{2}}}}={\begin{bmatrix}1&0&0\\0&-1&0\\0&0&1\end{bmatrix}}~;~~{\underline {\underline {\mathbf {A} _{3}}}}={\begin{bmatrix}1&0&0\\0&1&0\\0&0&-1\end{bmatrix}}}

Se puede demostrar que si la matriz C _ _ {\displaystyle {\underline {\underline {\mathsf {C}}}}} para un material elástico es invariante frente a una reflexión respecto de dos planos ortogonales, entonces también es invariante frente a una reflexión respecto del tercer plano ortogonal. Si consideraramos otra reflexión A 3 _ _ {\displaystyle {\underline {\underline {\mathbf {A} _{3}}}}} respecto del plano 1 2 {\displaystyle 1-2} entonces tendríamos:

A ε _ _ = [ 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 ] {\displaystyle {\underline {\underline {{\mathsf {A}}_{\varepsilon }}}}={\begin{bmatrix}1&0&0&0&0&0\\0&1&0&0&0&0\\0&0&1&0&0&0\\0&0&0&-1&0&0\\0&0&0&0&-1&0\\0&0&0&0&0&1\end{bmatrix}}}

La necesidad de que C _ _ = A ε _ _ T   C _ _   A ε _ _ {\displaystyle {\underline {\underline {\mathsf {C}}}}={\underline {\underline {{\mathsf {A}}_{\varepsilon }}}}^{T}~{\underline {\underline {\mathsf {C}}}}~{\underline {\underline {{\mathsf {A}}_{\varepsilon }}}}} implica que:[3]

[ 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 ] = [ 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 ] {\displaystyle {\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}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}}}

Esto se puede cumplir únicamente si:

C 14 = C 15 = C 24 = C 25 = C 34 = C 35 = C 46 = C 56 = 0   . {\displaystyle C_{14}=C_{15}=C_{24}=C_{25}=C_{34}=C_{35}=C_{46}=C_{56}=0~.}

Consideremos ahora la reflexión A 2 _ _ {\displaystyle {\underline {\underline {\mathbf {A} _{2}}}}} respecto del plano 1 3 {\displaystyle 1-3} . En este caso:

A ε _ _ = [ 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 ] {\displaystyle {\underline {\underline {{\mathsf {A}}_{\varepsilon }}}}={\begin{bmatrix}1&0&0&0&0&0\\0&1&0&0&0&0\\0&0&1&0&0&0\\0&0&0&-1&0&0\\0&0&0&0&1&0\\0&0&0&0&0&-1\end{bmatrix}}}

Usando de nuevo la condición de invarianza, se obtiene un requerimiento adicional:

C 16 = C 26 = C 36 = C 45 = 0 {\displaystyle C_{16}=C_{26}=C_{36}=C_{45}=0}

No se puede obtener información adicional porque la reflexión respecto del tercer plano de simetría no es independiente de las reflexiones respecto de los plano que ya hemos considerado. Por lo tanto, la matriz de rigidez o tensor de constantes elásticas de un material elástico ortótropo lineal puede ser escrita como:

C _ _ = [ C 11 C 12 C 13 0 0 0 C 12 C 22 C 23 0 0 0 C 13 C 23 C 33 0 0 0 0 0 0 C 44 0 0 0 0 0 0 C 55 0 0 0 0 0 0 C 66 ] {\displaystyle {\underline {\underline {\mathsf {C}}}}={\begin{bmatrix}C_{11}&C_{12}&C_{13}&0&0&0\\C_{12}&C_{22}&C_{23}&0&0&0\\C_{13}&C_{23}&C_{33}&0&0&0\\0&0&0&C_{44}&0&0\\0&0&0&0&C_{55}&0\\0&0&0&0&0&C_{66}\end{bmatrix}}}

La inversa de esta matriz de rigidez, denominada matriz de flexibilidad o matriz de docilidad, generalmente se escribe como:[4]

S _ _ = [ 1 E 1 ν 12 E 1 ν 13 E 1 0 0 0 ν 21 E 2 1 E 2 ν 23 E 2 0 0 0 ν 31 E 3 ν 32 E 3 1 E 3 0 0 0 0 0 0 1 G 23 0 0 0 0 0 0 1 G 31 0 0 0 0 0 0 1 G 12 ] {\displaystyle {\underline {\underline {\mathsf {S}}}}={\begin{bmatrix}{\tfrac {1}{E_{\rm {1}}}}&-{\tfrac {\nu _{\rm {12}}}{E_{\rm {1}}}}&-{\tfrac {\nu _{\rm {13}}}{E_{\rm {1}}}}&0&0&0\\-{\tfrac {\nu _{\rm {21}}}{E_{\rm {2}}}}&{\tfrac {1}{E_{\rm {2}}}}&-{\tfrac {\nu _{\rm {23}}}{E_{\rm {2}}}}&0&0&0\\-{\tfrac {\nu _{\rm {31}}}{E_{\rm {3}}}}&-{\tfrac {\nu _{\rm {32}}}{E_{\rm {3}}}}&{\tfrac {1}{E_{\rm {3}}}}&0&0&0\\0&0&0&{\tfrac {1}{G_{\rm {23}}}}&0&0\\0&0&0&0&{\tfrac {1}{G_{\rm {31}}}}&0\\0&0&0&0&0&{\tfrac {1}{G_{\rm {12}}}}\\\end{bmatrix}}}

donde E i {\displaystyle {E}_{\rm {i}}\,} es el módulo de Young a lo largo del eje i {\displaystyle i} , G i j {\displaystyle G_{\rm {ij}}\,} es el módulo de corte en la dirección j {\displaystyle j} en el plano cuya normal está en la dirección i {\displaystyle i} , y ν i j {\displaystyle \nu _{\rm {ij}}\,} es el coeficiente de Poisson que corresponde a una contracción en la dirección j {\displaystyle j} cuando se aplica una elongación en la dirección i {\displaystyle i} .

Usando las constantes que aparecen en esa última matriz, la matriz de rigidez o matriz de constantes elásticas C _ _ {\displaystyle {\underline {\underline {\mathsf {C}}}}} se puede escribir explícitamente como:

C _ _ = 1 Δ ( ( 1 ν 23 ν 32 ) E 1 ( ν 12 + ν 13 ν 32 ) E 2 ( ν 13 + ν 12 ν 23 ) E 3 0 0 0 ( ν 21 + ν 23 ν 31 ) E 1 ( 1 ν 31 ν 13 ) E 2 ( ν 23 + ν 21 ν 13 ) E 3 0 0 0 ( ν 31 + ν 32 ν 21 ) E 1 ( ν 32 + ν 31 ν 12 ) E 2 ( 1 ν 12 ν 21 ) E 3 0 0 0 0 0 0 G 12 0 0 0 0 0 0 G 31 0 0 0 0 0 0 G 23 ) {\displaystyle {\underline {\underline {\mathsf {C}}}}={\frac {1}{\Delta }}{\begin{pmatrix}(1-\nu _{23}\nu _{32})E_{1}&(\nu _{12}+\nu _{13}\nu _{32})E_{2}&(\nu _{13}+\nu _{12}\nu _{23})E_{3}&0&0&0\\(\nu _{21}+\nu _{23}\nu _{31})E_{1}&(1-\nu _{31}\nu _{13})E_{2}&(\nu _{23}+\nu _{21}\nu _{13})E_{3}&0&0&0\\(\nu _{31}+\nu _{32}\nu _{21})E_{1}&(\nu _{32}+\nu _{31}\nu _{12})E_{2}&(1-\nu _{12}\nu _{21})E_{3}&0&0&0\\0&0&0&G_{12}&0&0\\0&0&0&0&G_{31}&0\\0&0&0&0&0&G_{23}\end{pmatrix}}}

Además la magnitud Δ {\displaystyle \Delta } viene dada por la relación:

Δ := 1 ν 12 ν 21 ν 13 ν 31 ν 23 ν 32 ν 13 ν 32 ν 21 ν 31 ν 12 ν 23 {\displaystyle \Delta :=1-\nu _{12}\nu _{21}-\nu _{13}\nu _{31}-\nu _{23}\nu _{32}-\nu _{13}\nu _{32}\nu _{21}-\nu _{31}\nu _{12}\nu _{23}}

Cotas para los módulos elásticos en materiales ortótropos

La relación deformación-tensión para un material lineal isótropo se puede escribir usando la notación de Voigt como:

ε = S : σ {\displaystyle {\boldsymbol {\varepsilon }}={\mathsf {S}}:{\boldsymbol {\sigma }}}

donde la matriz de flexibilidad S _ _ {\displaystyle {\underline {\underline {\mathsf {S}}}}} viene dada por:

S _ _ = [ S 11 S 12 S 13 0 0 0 S 12 S 22 S 23 0 0 0 S 13 S 23 S 33 0 0 0 0 0 0 S 44 0 0 0 0 0 0 S 55 0 0 0 0 0 0 S 66 ] {\displaystyle {\underline {\underline {\mathsf {S}}}}={\begin{bmatrix}S_{11}&S_{12}&S_{13}&0&0&0\\S_{12}&S_{22}&S_{23}&0&0&0\\S_{13}&S_{23}&S_{33}&0&0&0\\0&0&0&S_{44}&0&0\\0&0&0&0&S_{55}&0\\0&0&0&0&0&S_{66}\end{bmatrix}}}

La matriz de complianza es una matriz simétrica y debe ser definida positiva para que la función de densidad de energía de deformación sea positiva. Esto implica a partir del criterio de Sylvester que todos los menores principales de la matriz son positivos,[5]​ es decir:: Δ k := det ( S k _ _ ) > 0 {\displaystyle \Delta _{k}:=\det({\underline {\underline {{\mathsf {S}}_{k}}}})>0} donde S k _ _ {\displaystyle {\underline {\underline {{\mathsf {S}}_{k}}}}} es la submatriz principal de k × k {\displaystyle k\times k} de S _ _ {\displaystyle {\underline {\underline {\mathsf {S}}}}} .

Entonces: Δ 1 > 0 S 11 > 0 Δ 2 > 0 S 11 S 22 S 12 2 > 0 Δ 3 > 0 ( S 11 S 22 S 12 2 ) S 33 S 11 S 23 2 + 2 S 12 S 23 S 13 S 22 S 13 2 > 0 Δ 4 > 0 S 44 Δ 3 > 0 S 44 > 0 Δ 5 > 0 S 44 S 55 Δ 3 > 0 S 55 > 0 Δ 6 > 0 S 44 S 55 S 66 Δ 3 > 0 S 66 > 0 {\displaystyle {\begin{aligned}\Delta _{1}>0&\implies \quad S_{11}>0\\\Delta _{2}>0&\implies \quad S_{11}S_{22}-S_{12}^{2}>0\\\Delta _{3}>0&\implies \quad (S_{11}S_{22}-S_{12}^{2})S_{33}-S_{11}S_{23}^{2}+2S_{12}S_{23}S_{13}-S_{22}S_{13}^{2}>0\\\Delta _{4}>0&\implies \quad S_{44}\Delta _{3}>0\implies S_{44}>0\\\Delta _{5}>0&\implies \quad S_{44}S_{55}\Delta _{3}>0\implies S_{55}>0\\\Delta _{6}>0&\implies \quad S_{44}S_{55}S_{66}\Delta _{3}>0\implies S_{66}>0\end{aligned}}} Se puede demostrar que este conjunto de condiciones implica que:[6] S 11 > 0   ,     S 22 > 0   ,     S 33 > 0   ,     S 44 > 0   ,     S 55 > 0   ,     S 66 > 0 {\displaystyle S_{11}>0~,~~S_{22}>0~,~~S_{33}>0~,~~S_{44}>0~,~~S_{55}>0~,~~S_{66}>0} o: E 1 > 0 , E 2 > 0 , E 3 > 0 , G 12 > 0 , G 23 > 0 , G 13 > 0 {\displaystyle E_{1}>0,E_{2}>0,E_{3}>0,G_{12}>0,G_{23}>0,G_{13}>0} Sin embargo, no se pueden establecer cotas inferiores similares para los valores de los coeficientes de Poisson ν i j {\displaystyle \nu _{ij}} .[5]

Véase también

Referencias

  1. Milton, G. W., 2002, The Theory of Composites, Cambridge University Press.
  2. Lekhnitskii, S. G., 1963, Theory of Elasticity of an Anisotropic Elastic Body, Holden-Day Inc.
  3. a b c d Slawinski, M. A., 2010, Waves and Rays in Elastic Continua: 2nd Ed., World Scientific. https://web.archive.org/web/20110718145226/http://samizdat.mines.edu/wavesandrays/WavesAndRays.pdf
  4. Boresi, A. P, Schmidt, R. J. and Sidebottom, O. M., 1993, Advanced Mechanics of Materials, Wiley.
  5. a b Ting, T. C. T. and Chen, T., 2005, Poisson's ratio for anisotropic elastic materials can have no bounds,, Q. J. Mech. Appl. Math., 58(1), pp. 73-82.
  6. Ting, T. C. T., 1996, Positive definiteness of anisotropic elastic constants,, Mathematics and Mechanics of Solids, 1, pp. 301-314.

Enlaces externos

  • Orthotropy modeling equations from OOFEM Matlib manual section.
  • Hooke's law for orthotropic materials
Control de autoridades
  • Proyectos Wikimedia
  • Wd Datos: Q901006
  • Wd Datos: Q901006