Teoría de Mindlin–Reissner

La deformación de un elemento de placa, en el que se resalta la superficie media (rojo) y la normal a la superficie (azul).

La teoría de Mindlin–Reissner (denominada también como teoría de placas Mindlin-Reisner) es parte de las teoría de placas y láminas en el estudio de mecánica de sólidos deformables. Se trata de una variante extendida de la teoría de Kirchhoff–Love aplicada a las placas delgadas que considera además los esfuerzos cortantes. La teoría fue propuesta inicialmente en el año 1951 por Raymond Mindlin.​ No obstante, una idea similar fue propuesta anteriormente por Eric Reissner en 1945.​ Ambas teorías eran un intento de explicar la deformación de la zona intermedia de una lámina delgada sujeta a flexión, en la que la normal a la superficie media permanece recta, pero no necesariamente perpendicular a la superficie. La teoría de Mindlin-Reissner se emplea en el cálculo de las deformaciones y tensiones en una placa cuyo grosor es de un décimo que las dimensiones de superficie mientras que la teoría de Kirchhoff-Love solo se aplica a placas muy delgadas.

La teoría de losas de Mindlin-Reissner se usa también para el análisis de placas laminadas pues en este tipo de placas las deformaciones transversales son muy importantes. Esta teoría también se extiende al análisis de láminas delgadas pues son muy versátiles para el análisis de placas tanto gruesas como delgadas y con materiales con distribución homogénea o materiales compuestos.

Generalidades de la teoría de Reissner-Mindlin

En esta teoría se relajan las condiciones para el movimiento en el plano de la placa, así que el movimiento en el plano $xy$ se considera en los cálculos. Estos movimientos quedan descritos de la siguiente manera:

$\theta _{x}\approx \tan \theta _{x}={\frac {\partial w(w,y)}{\partial x}}+\phi _{x}$

$\theta _{y}\approx \tan \theta _{y}={\frac {\partial w(x,y)}{\partial y}}+\phi _{y}$

Donde $\phi _{x}$ y $\phi _{y}$ representan la distorsión angular que sufre la sección transversal al no permanecer ortogonal al plano medio de la placa. Esta formulación solamente es válida para relaciones $h/L<0.2$

Hipótesis de la teoría de Reissner-Mindlin

La teoría de placas de Reissner-Mindlin tiene las siguientes suposiciones:

1. El dominio en el que se desarrolla la teoría tiene la siguiente forma:

$\Omega =\{(x,y,z)\in \mathbb {R} ^{3}|z\in \left[{\frac {-t}{2}},{\frac {t}{2}}\right],(x,y)\in A\subset \mathbb {R} ^{2}\}$

Donde $t$ representa el grosor de la placa.

1. Los esfuerzos normales en la dirección z son equivalentes a cero. Así:

$\sigma _{z}(x,y,z)=0$

Esta hipótesis está en contradicción aparente con la segunda hipótesis pero en realidad no causa problemas, pues se supone que no se aplican cargas puntuales.

1. El desplazamiento en la dirección z es equivalente a w(x,y), de esta manera el desplazamiento a lo largo del eje vertical quedaría expresado como:

$u_{z}(x,y,z)=w(x,y)$

. Esto quiere decir que el desplazamiento vertical no cambia a lo largo del grosor de la placa

1. El desplazamiento en dirección x y y depende de la rotación alrededor el eje de la siguiente manera:

$u_{x}=u=-z\theta _{x}(x,y)$

$u_{y}=v=-z\theta _{y}(x,y)$

Esto quiere decir que la sección transversal antes de la deformación permanece plana pero no ortogonal al plano medio de la viga después de la deformación. $\theta$ representa el giro de la sección normal al eje neutro de la viga, esto es $z=0$

Por simplicidad se considera que el materia es isotrópico y linealmente elástico.

Definición del campo de desplazamientos y vector de movimientos

El campo de desplazamientos queda definido de la siguiente manera:

${\textbf {u}}={\begin{bmatrix}u_{1}\\u_{2}\\u_{3}\\\end{bmatrix}}={\begin{bmatrix}u(x,y,z)\\v(x,y,z)\\w(x,y,z)\\\end{bmatrix}}={\begin{bmatrix}-z\theta _{x}\\-z\theta _{y}\\w(x,y)\\\end{bmatrix}}$

Donde: $\tan \theta _{x}\approx \theta _{x}$ y $\tan \theta _{y}\approx \theta _{y}$ Para pequeñas deformaciones, cumpliendo así las hipótesis de la teoría en cuestión. De acuerdo a todo lo formulado con anterioridad se puede definir el vector de movimientos que representa el movimiento del plano medio de la losa (es decir z=0) de la siguiente manera:

${\textbf {u}}={\begin{bmatrix}w\\\theta _{x}\\\theta _{y}\\\end{bmatrix}}$

Con toda la formulación que se tiene se puede observar $\theta _{x}$ y $\theta _{y}$ no dependen solamente de la deflexion, por lo tanto se puede tratar como una variable independiente.

Definición de las curvaturas de la losa

La matriz de curvaturas y giro de la losa se definen de manera matricial como:

$\kappa ={\begin{bmatrix}\kappa _{x}\\\kappa _{y}\\\kappa _{xy}\\\end{bmatrix}}={\begin{bmatrix}{\frac {\partial \theta _{x}(x,y)}{\partial x}}\\{\frac {\partial \theta _{y}(x,y)}{\partial y}}\\{\frac {\partial \theta _{x}(x,y)}{\partial x}}+{\frac {\partial \theta _{y}(x,y)}{\partial y}}\\\end{bmatrix}}$

Cabe anotar que $\kappa _{i}>0$ representa una concavidad positiva (estiramiento del plano inferior de la placa), $\kappa _{i}<0$ una concavidad negativa (estiramiento del plano inferior de la placa) y $\kappa _{i}=0$ que no existe deformación en la viga (véase Concavidad).

Definición del campo de deformaciones

Las deformaciones según la teoría de Mindlin-Reissner estarían dadas de la siguiente manera:

$\varepsilon _{x}={\frac {\partial u}{\partial x}}=-z{\frac {\partial \theta _{x}}{\partial x}}$

$\varepsilon _{y}={\frac {\partial v}{\partial y}}=-z{\frac {\partial \theta _{y}}{\partial y}}$

$\varepsilon _{z}={\frac {\partial w}{\partial z}}=0$

se tiene que $\varepsilon _{z}=0$ debido a que el desplazamiento vertical no depende de la posición del punto a lo largo del plano transversal de la viga.

Las deformaciones angulares deben cumplir con las hipótesis planteadas anteriormente de la siguiente manera:

$\gamma _{xy}={\frac {\partial u}{\partial y}}+{\frac {\partial v}{\partial x}}={\frac {\partial (-z\theta _{x})}{\partial y}}+{\frac {\partial (-z\theta _{y})}{\partial x}}=-2z{\frac {\partial ^{2}w}{\partial x\partial y}}-z\left({\frac {\partial \phi _{y}}{\partial x}}+{\frac {\partial \phi _{x}}{\partial y}}\right)$

$\gamma _{xz}={\frac {\partial u}{\partial z}}+{\frac {\partial w}{\partial x}}=-\phi _{x}$

$\gamma _{yz}={\frac {\partial v}{\partial z}}+{\frac {\partial w}{\partial y}}=-\phi _{y}$

Por lo tanto el campo de deformaciones queda definido de la siguiente manera:

$\varepsilon ={\begin{bmatrix}\varepsilon _{x}\\\varepsilon _{y}\\\gamma _{xy}\\\gamma _{xz}\\\gamma _{yz}\\\end{bmatrix}}={\begin{bmatrix}\varepsilon _{f}\\\varepsilon _{c}\end{bmatrix}}$

Donde $\varepsilon _{f}={\begin{bmatrix}\varepsilon _{x},\varepsilon _{y},\gamma _{xy}\end{bmatrix}}^{T}$ representan el vector de deformaciones por flexion y $\varepsilon _{c}={\begin{bmatrix}\gamma _{xz},\gamma _{yz}\end{bmatrix}}^{T}$ el vector de deformaciones por cortante.

Definición del campo de esfuerzos de la losa

Los esfuerzos que se tienen en cuenta según la teoría de Mindlin-Reissner son los siguientes:

$\sigma ={\begin{bmatrix}\sigma _{x}\\\sigma _{y}\\\sigma _{z}\\\tau _{xy}\\\tau _{xz}\\\tau _{yz}\\\end{bmatrix}}={\begin{bmatrix}\sigma _{x}\\\sigma _{y}\\\tau _{xy}\\\tau _{xz}\\\tau _{yz}\\\end{bmatrix}}={\begin{bmatrix}\sigma _{f}\\\sigma _{c}\end{bmatrix}}$

Donde $\sigma _{f}={\begin{bmatrix}\sigma _{x},\sigma _{y},\tau _{xy}\end{bmatrix}}^{T}$ representan el vector de esfuerzos debidos a flexion y $\sigma _{c}={\begin{bmatrix}\tau _{xz},\tau _{yz}\end{bmatrix}}^{T}$ el vector de esfuerzos debidos a cortante. Se toma a $\sigma _{z}=0$ según la hipótesis 2 de la teoría de losas de Mindlin-Reissner.

Relación entre esfuerzos y deformaciones en la losa

Ahora, considérese una losa con un material homogéneo y en condiciones isotermales. La matriz constitutiva D usada en esta teoría corresponde a una combinación de la matriz constitutiva de tensión plana en 2D y la matriz de corte escrita a continuación:

$D={\begin{bmatrix}D_{f},0\\0,D_{c}\\\end{bmatrix}}$

En donde la matriz $D_{f}$ se denomina matriz constitutiva por flexión y está definida como:

$D_{f}={\frac {E}{1-\nu ^{2}}}{\begin{bmatrix}1,\nu ,0\\\nu ,1,0\\0,0,{\frac {1-\nu }{2}}\\\end{bmatrix}}$

y la matriz $D_{s}$ se denomina matriz constitutiva de cortante y está definida como:

$D_{c}=G{\begin{bmatrix}1,0\\0,1\end{bmatrix}}$

Ahora, al igual que en la Teoría de Vigas de Timoshenko los esfuerzos $\tau {xz}$ y $\tau _{yz}$ no tienen la distribución parabólica que estos tienen en la teoría, por lo tanto se trabaja con un $\alpha$ denominado coeficiente de distorsión lateral que hace que tanto la distribución teórica de $\tau _{xz}$ y $\tau _{yz}$ como la distribución hallada mediante la teoría de Mindlin-Reissner tengan el mismo trabajo de deformación. Así, para losas el coeficiente de distorsión está determinado por 5/6 ($\alpha =5/6$). Así, la matriz constitutiva por cortante queda definida de la siguiente manera:

$D_{s}=G\alpha {\begin{bmatrix}1,0\\0,1\end{bmatrix}}$

Por lo tanto, el campo de esfuerzos quedaría descrita de la siguiente manera de acuerdo a la ley de elasticidad de Hooke:

$\varepsilon =D\sigma$

Aplicando la teoría de losas tratada a la ley de elasticidad de Hooke se obtiene el siguiente resultado:

$\varepsilon =D\sigma ={\begin{bmatrix}D_{f},0\\0,D_{c}\\\end{bmatrix}}\sigma ={\begin{bmatrix}D_{f},0\\0,D_{c}\\\end{bmatrix}}{\begin{bmatrix}\sigma _{x}\\\sigma _{y}\\\tau _{xy}\\\tau _{xz}\\\tau _{yz}\\\end{bmatrix}}$

Por lo tanto la relación entre esfuerzos y deformaciones hacen que se puedan clasificar los esfuerzos de la siguiente manera: existen unos esfuerzos de flexión y unas deformaciones por flexión que están descritos de la siguiente manera:

$\varepsilon _{f}=D_{f}\sigma _{f}$

desglosando la ecuación se tiene que:

${\begin{bmatrix}\varepsilon _{x}\\\varepsilon _{y}\\\gamma _{xy}\end{bmatrix}}={\frac {E}{1-\nu ^{2}}}{\begin{bmatrix}1,\nu ,0\\\nu ,1,0\\0,0,{\frac {1-\nu }{2}}\\\end{bmatrix}}{\begin{bmatrix}\sigma _{x}\\\sigma _{y}\\\sigma _{xy}\end{bmatrix}}$

Y de la misma manera existen los esfuerzos por cortante y las deformaciones transversales por cortante que están relacionados de la siguiente manera:

$\varepsilon _{c}=D_{c}\sigma _{c}$

Y de manera matricial quedarían descritas de la siguiente manera:

${\begin{bmatrix}\gamma _{xz}\\\gamma _{xz}\end{bmatrix}}=G\alpha {\begin{bmatrix}1,0\\0,1\end{bmatrix}}{\begin{bmatrix}\tau _{xz}\\\tau _{yz}\end{bmatrix}}$

Por lo tanto las deformaciones se puede escribir de la siguiente manera:

$\varepsilon ={\begin{bmatrix}\varepsilon _{f}\\\varepsilon _{c}\end{bmatrix}}$

Y la relación entre esfuerzos y deformaciones estaría descrita a continuación:

$\sigma =D\varepsilon ={\begin{bmatrix}D_{f},0\\0,D_{c}\end{bmatrix}}{\begin{bmatrix}\varepsilon _{f}\\\varepsilon _{c}\end{bmatrix}}$

Vector de deformaciones generalizadas

El vector de deformaciones $\varepsilon$ se puede escribir de la siguiente manera:

${\begin{bmatrix}\varepsilon _{x}\\\varepsilon _{y}\\\tau _{xy}\\\tau _{xz}\\\tau _{yz}\end{bmatrix}}={\begin{bmatrix}-z{\frac {\partial \theta _{x}}{\partial x}}\\-z{\frac {\partial \theta _{y}}{\partial y}}\\-2z{\frac {\partial ^{2}w}{\partial x\partial y}}-z\left({\frac {\partial \phi _{y}}{\partial x}}+{\frac {\partial \phi _{x}}{\partial y}}\right)\end{bmatrix}}=-z{\begin{bmatrix}{\frac {\partial \theta _{x}}{\partial x}}\\{\frac {\partial \theta _{y}}{\partial y}}\\2{\frac {\partial ^{2}w}{\partial x\partial y}}-\left({\frac {\partial \phi _{y}}{\partial x}}+{\frac {\partial \phi _{x}}{\partial y}}\right)\end{bmatrix}}=-z{\hat {\varepsilon _{f}}}$

Donde ${\hat {\varepsilon _{f}}}$ representa el vector de deformaciones por flexión generalizadas

Vector de esfuerzos generalizados

Los esfuerzos se pueden expresar a través de un vector de esfuerzos generalizados ${\hat {\sigma }}$ de la siguiente manera:

${\hat {\sigma }}={\begin{bmatrix}{\hat {\sigma _{f}}}\\{\hat {\sigma _{c}}}\end{bmatrix}}={\begin{bmatrix}M_{x}\\M_{y}\\M_{xy}\\Q_{x}\\Q_{y}\end{bmatrix}}=\int _{-t/2}^{t/2}{\begin{bmatrix}-z\sigma _{x}\\-z\sigma _{y}\\-z\tau _{xy}\\\tau _{xz}\\\tau _{yz}\end{bmatrix}}dz=\int _{-t/2}^{t/2}{\begin{bmatrix}z\sigma _{f}\\\sigma _{c}\end{bmatrix}}dz$

Pero se tiene que $\sigma _{f}=D_{f}\varepsilon _{f}=-zD_{f}{\hat {\varepsilon _{f}}}$ y $\sigma _{c}={\hat {\sigma _{c}}}=D_{c}{\hat {\varepsilon _{c}}}$, reemplazando estos valores en la anterior ecuación se tiene que:

${\hat {\sigma }}=\int _{-t/2}^{t/2}{\begin{bmatrix}zD_{f}z{\hat {\varepsilon _{f}}}\\D_{c}{\hat {\varepsilon _{c}}}\end{bmatrix}}dz={\begin{bmatrix}\int _{-t/2}^{t/2}z^{2}dzD_{f}{\hat {\varepsilon _{f}}}\\\int _{-t/2}^{t/2}dzD_{c}{\hat {\varepsilon _{c}}}\end{bmatrix}}={\begin{bmatrix}{\frac {t^{3}}{12}}D_{f}{\hat {\varepsilon _{f}}}\\tD_{c}{\hat {\varepsilon _{c}}}\end{bmatrix}}={\begin{bmatrix}{\hat {D_{f}}}{\hat {\varepsilon _{f}}}\\{\hat {D_{c}}}{\hat {\varepsilon _{c}}}\end{bmatrix}}$

Las matrices ${\hat {D_{f}}}={\frac {t^{3}}{12}}D_{f}$ y ${\hat {D_{c}}}=tD_{c}$ se denominan matriz constitutiva generalizada de flexión y matriz constitutiva generalizada de cortante respectivamente.

Relación entre esfuerzos generalizados y deformaciones generalizadas

Así, los esfuerzos generalizados y las deformaciones generalizadas quedarían de la siguiente manera:

${\hat {\sigma }}={\begin{bmatrix}{\hat {\sigma _{f}}}\\{\hat {\sigma _{c}}}\end{bmatrix}}={\begin{bmatrix}{\hat {D_{f}}},0\\0,{\hat {D_{c}}}\end{bmatrix}}={\begin{bmatrix}{\hat {\varepsilon _{f}}}\\{\hat {\varepsilon _{c}}}\end{bmatrix}}={\hat {D}}{\hat {\varepsilon }}$

Bibliografía

Hughes, Thomas J.R (1987). «5». The Finite Element Method (en inglés). Dover Publications. pp. 311. (requiere registro). 

Reddy, J.N (2006). «6». And introduction to the finite element method (en inglés) (3 edición). Mc Graw Hill. pp. 311. (requiere registro). 

Oñate, Eugenio (213). «6». Structural Analysis with the Finite Element Method Linear Statics (en inglés) (1 edición). Springer. 

Véase también

• Teoría de Vigas de Timoshenko - Teoría de vigas que tiene en cuenta las deformaciones por cortante.
• Ley de Hooke
• [1]
• [2]
• [3]

• Datos: Q3564249