Menu
Home
People
Places
Arts
History
Plants & Animals
Science
Life & Culture
Technology
Reference.org
Orthotropic material
open-in-new
<h2 id="orthotropy-in-physics">Orthotropy in physics</h2> <h3>Anisotropic material relations</h3> <p>Material behavior is represented in physical theories by <a href="/facts/Constitutive_relation/IXTjEAbL">constitutive relations</a>. A large class of physical behaviors can be represented by linear material models that take the form of a second-order <a href="/facts/Tensor/pNjFo5tV">tensor</a>. The material tensor provides a relation between two <a href="/facts/Euclidean_vector/bCLrdksy">vectors</a> and can be written as </p> f = K ⋅ d {\displaystyle \mathbf {f} ={\boldsymbol {K}}\cdot \mathbf {d} } <p>where d , f {\displaystyle \mathbf {d} ,\mathbf {f} } are two vectors representing physical quantities and K {\displaystyle {\boldsymbol {K}}} is the second-order material tensor. If we express the above equation in terms of components with respect to an <a href="/facts/Orthonormal/JfgL1nAh">orthonormal</a> <a href="/facts/Coordinate_system/Mtd0q4SP">coordinate system</a>, we can write </p> f i = K i j d j . {\displaystyle f_{i}=K_{ij}~d_{j}~.} <p><a href="/facts/Einstein_notation/UXwS43wU">Summation over repeated indices</a> has been assumed in the above relation. In matrix form we have </p> 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}}} <p>Examples of physical problems that fit the above template are listed in the table below.<a class="footnote-ref" id="fnref:2" href="#fn:2"><sup>2</sup></a> </p> <table><tbody><tr><th>Problem</th><th> f {\displaystyle \mathbf {f} } </th><th> d {\displaystyle \mathbf {d} } </th><th> K {\displaystyle {\boldsymbol {K}}} </th></tr><tr><td><a href="/facts/Electrical_conduction/akyCKvMy">Electrical conduction</a></td><td><a href="/facts/Electrical_current/Dc1xdnUX">Electrical current</a> J {\displaystyle \mathbf {J} } </td><td><a href="/facts/Electric_field/tt4MooBs">Electric field</a> E {\displaystyle \mathbf {E} } </td><td><a href="/facts/Electrical_conductivity/akyCKvMy">Electrical conductivity</a> σ {\displaystyle {\boldsymbol {\sigma }}} </td></tr><tr><td><a href="/facts/Dielectric/PnVt9P2G">Dielectrics</a></td><td><a href="/facts/Electrical_displacement/YERMcaRg">Electrical displacement</a> D {\displaystyle \mathbf {D} } </td><td><a href="/facts/Electric_field/tt4MooBs">Electric field</a> E {\displaystyle \mathbf {E} } </td><td><a href="/facts/Electric_permittivity/Sd3WtgT1">Electric permittivity</a> ε {\displaystyle {\boldsymbol {\varepsilon }}} </td></tr><tr><td><a href="/facts/Magnetism/eo75vApl">Magnetism</a></td><td><a href="/facts/Electromagnetic_induction/T36JjtgY">Magnetic induction</a> B {\displaystyle \mathbf {B} } </td><td><a href="/facts/Magnetic_field/Ta8Ev588">Magnetic field</a> H {\displaystyle \mathbf {H} } </td><td><a href="/facts/Magnetic_permeability/GA7CcRLU">Magnetic permeability</a> μ {\displaystyle {\boldsymbol {\mu }}} </td></tr><tr><td><a href="/facts/Thermal_conduction/5vRktgcS">Thermal conduction</a></td><td><a href="/facts/Heat_flux/YjF54Tcf">Heat flux</a> q {\displaystyle \mathbf {q} } </td><td><a href="/facts/Temperature_gradient/YhWj21zy">Temperature gradient</a> − ∇ T {\displaystyle -{\boldsymbol {\nabla }}T} </td><td><a href="/facts/Thermal_conductivity/qWnSG3nQ">Thermal conductivity</a> κ {\displaystyle {\boldsymbol {\kappa }}} </td></tr><tr><td><a href="/facts/Diffusion/KbAxroBa">Diffusion</a></td><td>Particle <a href="/facts/Flux/Aw3P5Bvb">flux</a> J {\displaystyle \mathbf {J} } </td><td><a href="/facts/Concentration_gradient/lcIz28qY">Concentration gradient</a> − ∇ c {\displaystyle -{\boldsymbol {\nabla }}c} </td><td><a href="/facts/Mass_diffusivity/GSngNem0">Diffusivity</a> D {\displaystyle {\boldsymbol {D}}} </td></tr><tr><td><a href="/facts/Fluid_dynamics/y7uzSCfg">Flow</a> in <a href="/facts/Porous_media/KTMRISJf">porous media</a></td><td>Weighted fluid <a href="/facts/Velocity/Q3o1LVDe">velocity</a> η μ v {\displaystyle \eta _{\mu }\mathbf {v} } </td><td><a href="/facts/Pressure_gradient/v7BmlAAS">Pressure gradient</a> ∇ P {\displaystyle {\boldsymbol {\nabla }}P} </td><td><a href="/facts/Fluid_permeability/ja7dRwgU">Fluid permeability</a> κ {\displaystyle {\boldsymbol {\kappa }}} </td></tr></tbody></table> <h3>Condition for material symmetry</h3> <p>The material matrix K _ _ {\displaystyle {\underline {\underline {\boldsymbol {K}}}}} has a symmetry with respect to a given <a href="/facts/Orthogonal_transformation/s7LGBDFX">orthogonal transformation</a> ( A {\displaystyle {\boldsymbol {A}}} ) if it does not change when subjected to that transformation. For invariance of the material properties under such a transformation we require </p> 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}}} <p>Hence the condition for material symmetry is (using the definition of an orthogonal transformation) </p> 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}}} <p>Orthogonal transformations can be represented in Cartesian coordinates by a 3 × 3 {\displaystyle 3\times 3} matrix A _ _ {\displaystyle {\underline {\underline {\boldsymbol {A}}}}} given by </p> 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}}~.} <p>Therefore, the symmetry condition can be written in matrix form as </p> K _ _ = A T _ _ K _ _ A _ _ {\displaystyle {\underline {\underline {\boldsymbol {K}}}}={\underline {\underline {{\boldsymbol {A}}^{T}}}}~{\underline {\underline {\boldsymbol {K}}}}~{\underline {\underline {\boldsymbol {A}}}}} <h3>Orthotropic material properties</h3> <p>An orthotropic material has three <a href="/facts/Orthogonal/JVIjYPog">orthogonal</a> <a href="/facts/Plane_of_symmetry/eSeR8nbD">symmetry planes</a>. If we choose an orthonormal coordinate system such that the axes coincide with the normals to the three symmetry planes, the transformation matrices are </p> 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}}} <p>It can be shown that if the matrix K _ _ {\displaystyle {\underline {\underline {\boldsymbol {K}}}}} for a material is invariant under reflection about two orthogonal planes then it is also invariant under reflection about the third orthogonal plane. </p><p>Consider the reflection A 3 _ _ {\displaystyle {\underline {\underline {{\boldsymbol {A}}_{3}}}}} about the 1 − 2 {\displaystyle 1-2\,} plane. Then we have </p> 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}}} <p>The above relation implies that K 13 = K 23 = K 31 = K 32 = 0 {\displaystyle K_{13}=K_{23}=K_{31}=K_{32}=0} . Next consider a reflection A 2 _ _ {\displaystyle {\underline {\underline {{\boldsymbol {A}}_{2}}}}} about the 1 − 3 {\displaystyle 1-3\,} plane. We then have </p> 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}}} <p>That implies that K 12 = K 21 = 0 {\displaystyle K_{12}=K_{21}=0} . Therefore, the material properties of an orthotropic material are described by the matrix </p> <blockquote> 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}}} </blockquote> <h2 id="orthotropy-in-linear-elasticity">Orthotropy in linear elasticity</h2> <h3>Anisotropic elasticity</h3> <p>In <a href="/facts/Linear_elasticity/pPLIvs7a">linear elasticity</a>, the relation between <a href="/facts/Stress_(physics)/ErAXDbbd">stress</a> and <a href="/facts/Infinitesimal_strain_theory/j4CJFdsU">strain</a> depend on the type of material under consideration. This relation is known as <a href="/facts/Hooke%2527s_law/7hTUcKUm">Hooke's law</a>. For anisotropic materials Hooke's law can be written as<a class="footnote-ref" id="fnref:3" href="#fn:3"><sup>3</sup></a> </p> σ = c ⋅ ε {\displaystyle {\boldsymbol {\sigma }}={\mathsf {c}}\cdot {\boldsymbol {\varepsilon }}} <p>where σ {\displaystyle {\boldsymbol {\sigma }}} is the stress <a href="/facts/Tensor/pNjFo5tV">tensor</a>, ε {\displaystyle {\boldsymbol {\varepsilon }}} is the strain tensor, and c {\displaystyle {\mathsf {c}}} is the elastic <a href="/facts/Stiffness_tensor/7hTUcKUm">stiffness tensor</a>. If the tensors in the above expression are described in terms of components with respect to an <a href="/facts/Orthonormal/JfgL1nAh">orthonormal</a> <a href="/facts/Coordinate_system/Mtd0q4SP">coordinate system</a> we can write </p> σ i j = c i j k ℓ ε k ℓ {\displaystyle \sigma _{ij}=c_{ijk\ell }~\varepsilon _{k\ell }} <p>where summation has been assumed over repeated indices. Since the stress and strain tensors are <a href="/facts/Symmetric_tensor/CWKBfgnv">symmetric</a>, and since the stress-strain relation in linear elasticity can be derived from a <a href="/facts/Strain_energy_density_function/PhKsKtGB">strain energy density function</a>, the following symmetries hold for linear elastic materials </p> 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}~.} <p>Because of the above symmetries, the stress-strain relation for linear elastic materials can be expressed in matrix form as </p> [ σ 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}}} <p>An alternative representation in <a href="/facts/Voigt_notation/EbE6JTGA">Voigt notation</a> is </p> [ σ 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}}} <p>or </p> σ _ _ = C _ _ ε _ _ {\displaystyle {\underline {\underline {\boldsymbol {\sigma }}}}={\underline {\underline {\mathsf {C}}}}~{\underline {\underline {\boldsymbol {\varepsilon }}}}} <p>The <a href="/facts/Stiffness_matrix/kjsEp3rj">stiffness matrix</a> C _ _ {\displaystyle {\underline {\underline {\mathsf {C}}}}} in the above relation satisfies <a href="/facts/Point_symmetry/UzswAm8N">point symmetry</a>.<a class="footnote-ref" id="fnref:4" href="#fn:4"><sup>4</sup></a> </p> <h3>Condition for material symmetry</h3> <p>The stiffness matrix C _ _ {\displaystyle {\underline {\underline {\mathsf {C}}}}} satisfies a given symmetry condition if it does not change when subjected to the corresponding <a href="/facts/Orthogonal_transformation/s7LGBDFX">orthogonal transformation</a>. The orthogonal transformation may represent symmetry with respect to a <a href="/facts/Point_symmetry/UzswAm8N">point</a>, an <a href="/facts/Axis_of_symmetry/bpvJqcnS">axis</a>, or a <a href="/facts/Plane_of_symmetry/eSeR8nbD">plane</a>. Orthogonal transformations in linear elasticity include rotations and reflections, but not shape changing transformations and can be represented, in orthonormal coordinates, by a 3 × 3 {\displaystyle 3\times 3} matrix A _ _ {\displaystyle {\underline {\underline {\mathbf {A} }}}} given by </p> 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}}~.} <p>In Voigt notation, the transformation matrix for the stress tensor can be expressed as a 6 × 6 {\displaystyle 6\times 6} matrix A σ _ _ {\displaystyle {\underline {\underline {{\mathsf {A}}_{\sigma }}}}} given by<a class="footnote-ref" id="fnref:5" href="#fn:5"><sup>5</sup></a> </p> 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}}} <p>The transformation for the strain tensor has a slightly different form because of the choice of notation. This transformation matrix is </p> 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}}} <p>It can be shown that A ε _ _ T = A σ _ _ − 1 {\displaystyle {\underline {\underline {{\mathsf {A}}_{\varepsilon }}}}^{T}={\underline {\underline {{\mathsf {A}}_{\sigma }}}}^{-1}} . </p> <blockquote> <p>The elastic properties of a continuum are invariant under an orthogonal transformation A _ _ {\displaystyle {\underline {\underline {\mathbf {A} }}}} if and only if<a class="footnote-ref" id="fnref:6" href="#fn:6"><sup>6</sup></a> </p> 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 }}}}} </blockquote> <h3>Stiffness and compliance matrices in orthotropic elasticity</h3> <p>An orthotropic elastic material has three <a href="/facts/Orthogonal/JVIjYPog">orthogonal</a> <a href="/facts/Plane_of_symmetry/eSeR8nbD">symmetry planes</a>. If we choose an orthonormal coordinate system such that the axes coincide with the normals to the three symmetry planes, the transformation matrices are </p> 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}}} <p>We can show that if the matrix C _ _ {\displaystyle {\underline {\underline {\mathsf {C}}}}} for a linear elastic material is invariant under reflection about two orthogonal planes then it is also invariant under reflection about the third orthogonal plane. </p><p>If we consider the reflection A 3 _ _ {\displaystyle {\underline {\underline {\mathbf {A} _{3}}}}} about the 1 − 2 {\displaystyle 1-2\,} plane, then we have </p> 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}}} <p>Then the requirement 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 }}}}} implies that<a class="footnote-ref" id="fnref:7" href="#fn:7"><sup>7</sup></a> </p> [ 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}}} <p>The above requirement can be satisfied only if </p> 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~.} <p>Let us next consider the reflection A 2 _ _ {\displaystyle {\underline {\underline {\mathbf {A} _{2}}}}} about the 1 − 3 {\displaystyle 1-3\,} plane. In that case </p> 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}}} <p>Using the invariance condition again, we get the additional requirement that </p> C 16 = C 26 = C 36 = C 45 = 0 . {\displaystyle C_{16}=C_{26}=C_{36}=C_{45}=0~.} <p>No further information can be obtained because the reflection about third symmetry plane is not independent of reflections about the planes that we have already considered. Therefore, the stiffness matrix of an orthotropic linear elastic material can be written as </p> <blockquote> 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}}} </blockquote> <p>The inverse of this matrix is commonly written as<a class="footnote-ref" id="fnref:8" href="#fn:8"><sup>8</sup></a> </p> S _ _ = [ 1 E 1 − ν 21 E 2 − ν 31 E 3 0 0 0 − ν 12 E 1 1 E 2 − ν 32 E 3 0 0 0 − ν 13 E 1 − ν 23 E 2 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 {21}}}{E_{\rm {2}}}}&-{\tfrac {\nu _{\rm {31}}}{E_{\rm {3}}}}&0&0&0\\-{\tfrac {\nu _{\rm {12}}}{E_{\rm {1}}}}&{\tfrac {1}{E_{\rm {2}}}}&-{\tfrac {\nu _{\rm {32}}}{E_{\rm {3}}}}&0&0&0\\-{\tfrac {\nu _{\rm {13}}}{E_{\rm {1}}}}&-{\tfrac {\nu _{\rm {23}}}{E_{\rm {2}}}}&{\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}}} <p>where E i {\displaystyle {E}_{\rm {i}}\,} is the <a href="/facts/Young%2527s_modulus/MV5JMCIM">Young's modulus</a> along axis i {\displaystyle i} , G i j {\displaystyle G_{\rm {ij}}\,} is the <a href="/facts/Shear_modulus/9e9IvJnz">shear modulus</a> in direction j {\displaystyle j} on the plane whose normal is in direction i {\displaystyle i} , and ν i j {\displaystyle \nu _{\rm {ij}}\,} is the <a href="/facts/Poisson%2527s_ratio/ppIsaujA">Poisson's ratio</a> that corresponds to a contraction in direction j {\displaystyle j} when an extension is applied in direction i {\displaystyle i} . </p> <h3>Bounds on the moduli of orthotropic elastic materials</h3> <p>The strain-stress relation for orthotropic linear elastic materials can be written in Voigt notation as </p> ε _ _ = S _ _ σ _ _ {\displaystyle {\underline {\underline {\boldsymbol {\varepsilon }}}}={\underline {\underline {\mathsf {S}}}}~{\underline {\underline {\boldsymbol {\sigma }}}}} <p>where the compliance matrix S _ _ {\displaystyle {\underline {\underline {\mathsf {S}}}}} is given by </p> 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}}} <p>The compliance matrix is <a href="/facts/Symmetric_matrix/08CGwjNl">symmetric</a> and must be <a href="/facts/Positive-definite_matrix/Hbr6AuPS">positive definite</a> for the <a href="/facts/Strain_energy_density_function/PhKsKtGB">strain energy density</a> to be positive. This implies from <a href="/facts/Sylvester%2527s_criterion/eSGHgbkS">Sylvester's criterion</a> that all the principal <a href="/facts/Minor_(linear_algebra)/CMudF5Eh">minors</a> of the matrix are positive,<a class="footnote-ref" id="fnref:9" href="#fn:9"><sup>9</sup></a> i.e., </p> Δ k := det ( S k _ _ ) > 0 {\displaystyle \Delta _{k}:=\det({\underline {\underline {{\mathsf {S}}_{k}}}})>0} <p>where S k _ _ {\displaystyle {\underline {\underline {{\mathsf {S}}_{k}}}}} is the k × k {\displaystyle k\times k} principal <a href="/facts/Submatrix/qa8FI4ko">submatrix</a> of S _ _ {\displaystyle {\underline {\underline {\mathsf {S}}}}} . </p><p>Then, </p> Δ 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}}} <p>We can show that this set of conditions implies that<a class="footnote-ref" id="fnref:10" href="#fn:10"><sup>10</sup></a> </p> 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} <p>or </p> 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} <p>However, no similar lower bounds can be placed on the values of the Poisson's ratios ν i j {\displaystyle \nu _{ij}} .<a class="footnote-ref" id="fnref:11" href="#fn:11"><sup>11</sup></a> </p> <h2 id="see-also">See also</h2> <ul><li><a href="/facts/Anisotropy/YMqTQWVB">Anisotropy</a></li> <li><a href="/facts/Clinotropic_material/SCgBWdPR">Clinotropic material</a></li> <li><a href="/facts/Stress_(mechanics)/ErAXDbbd">Stress (mechanics)</a></li> <li><a href="/facts/Infinitesimal_strain_theory/j4CJFdsU">Infinitesimal strain theory</a></li> <li><a href="/facts/Finite_strain_theory/rEjpC9oD">Finite strain theory</a></li> <li><a href="/facts/Hooke%2527s_law/7hTUcKUm">Hooke's law</a></li></ul> <h2 id="further-reading">Further reading</h2> <ul><li><a href="http://www.oofem.org/resources/doc/matlibmanual/html/node6.html">Orthotropy modeling equations</a> from <a href="/facts/OOFEM/I75SE3Vd">OOFEM</a> Matlib manual section.</li> <li><a href="http://www.efunda.com/formulae/solid_mechanics/mat_mechanics/hooke_orthotropic.cfm">Hooke's law for orthotropic materials</a></li></ul> <h2 id="references">References</h2> <ol> <li id="fn:1"><p>Geraldes DM et al, 2014, A comparative study of orthotropic and isotropic bone adaptation in the femur, International Journal for Numerical Methods in Biomedical Engineering, Volume 30, Issue 9, pages 873–889, DOI: 10.1002/cnm.2633, http://onlinelibrary.wiley.com/wol1/doi/10.1002/cnm.2633/full <a href="/wiki/International_Journal_for_Numerical_Methods_in_Biomedical_Engineering" target="_blank">/wiki/International_Journal_for_Numerical_Methods_in_Biomedical_Engineering</a> <a href="#fnref:1" class="footnote-back-ref">↩</a></p></li> <li id="fn:2"><p>Milton, G. W., 2002, The Theory of Composites, Cambridge University Press. <a href="#fnref:2" class="footnote-back-ref">↩</a></p></li> <li id="fn:3"><p>Lekhnitskii, S. G., 1963, Theory of Elasticity of an Anisotropic Elastic Body, Holden-Day Inc. <a href="#fnref:3" class="footnote-back-ref">↩</a></p></li> <li id="fn:4"><p>Slawinski, M. A., 2010, Waves and Rays in Elastic Continua: 2nd Ed., World Scientific. [1] <a href="https://web.archive.org/web/20090210192845/http://samizdat.mines.edu/wavesandrays/WavesAndRays.pdf" target="_blank">https://web.archive.org/web/20090210192845/http://samizdat.mines.edu/wavesandrays/WavesAndRays.pdf</a> <a href="#fnref:4" class="footnote-back-ref">↩</a></p></li> <li id="fn:5"><p>Slawinski, M. A., 2010, Waves and Rays in Elastic Continua: 2nd Ed., World Scientific. [1] <a href="https://web.archive.org/web/20090210192845/http://samizdat.mines.edu/wavesandrays/WavesAndRays.pdf" target="_blank">https://web.archive.org/web/20090210192845/http://samizdat.mines.edu/wavesandrays/WavesAndRays.pdf</a> <a href="#fnref:5" class="footnote-back-ref">↩</a></p></li> <li id="fn:6"><p>Slawinski, M. A., 2010, Waves and Rays in Elastic Continua: 2nd Ed., World Scientific. [1] <a href="https://web.archive.org/web/20090210192845/http://samizdat.mines.edu/wavesandrays/WavesAndRays.pdf" target="_blank">https://web.archive.org/web/20090210192845/http://samizdat.mines.edu/wavesandrays/WavesAndRays.pdf</a> <a href="#fnref:6" class="footnote-back-ref">↩</a></p></li> <li id="fn:7"><p>Slawinski, M. A., 2010, Waves and Rays in Elastic Continua: 2nd Ed., World Scientific. [1] <a href="https://web.archive.org/web/20090210192845/http://samizdat.mines.edu/wavesandrays/WavesAndRays.pdf" target="_blank">https://web.archive.org/web/20090210192845/http://samizdat.mines.edu/wavesandrays/WavesAndRays.pdf</a> <a href="#fnref:7" class="footnote-back-ref">↩</a></p></li> <li id="fn:8"><p>Boresi, A. P, Schmidt, R. J. and Sidebottom, O. M., 1993, Advanced Mechanics of Materials, Wiley. <a href="#fnref:8" class="footnote-back-ref">↩</a></p></li> <li id="fn:9"><p>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. <a href="#fnref:9" class="footnote-back-ref">↩</a></p></li> <li id="fn:10"><p>Ting, T. C. T. (1996), "Positive definiteness of anisotropic elastic constants", Mathematics & Mechanics of Solids, 1 (3): 301–314, doi:10.1177/108128659600100302, S2CID 122747373. <a href="/wiki/Mathematics_%26_Mechanics_of_Solids" target="_blank">/wiki/Mathematics_%26_Mechanics_of_Solids</a> <a href="#fnref:10" class="footnote-back-ref">↩</a></p></li> <li id="fn:11"><p>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. <a href="#fnref:11" class="footnote-back-ref">↩</a></p></li> </ol>