Teoría de placas y láminas

De Wikipedia, la enciclopedia libre
Saltar a: navegación, búsqueda
Flexión de una placa circular empotrada en su contorno bajo la acción de una carga vertical distribuida uniformemente. La mitad izquierda muestra la forma deformada y la mitad derecha muestra la forma no deformada. La simulación mediante elementos finitos fue llevada a cabo mediante el software Ansys.

En ingeniería estructural, las placas y las láminas son elementos estructurales que geométricamente se pueden aproximar por una superficie bidimensional y que trabajan predominantemente a flexión. Estructuralmente la diferencia entre placas y láminas está en la curvatura. Las placas son elementos cuya superficie media es plana, mientras que las láminas son superficies curvadas en el espacio tridimensional (como lás cúpulas, las conchas o las paredes de depósitos).

Constructivamente son sólidos deformables en los que existe una superficie media (que es la que se considera aproxima a la placa o lámina), a la que se añade un cierto espesor constante por encima y por debajo del plano medio. El hecho de que este espesor es pequeño comparado con las dimensiones de la lámina y a su vez pequeño comparado con los radios de curvatura de la superficie, es lo que permite reducir el cálculo de placas y láminas reales a elementos idealizados bidimensionales.

Cálculo de placas[editar]

Hipótesis de Reissner-Mindlin[editar]

Deformación transversal de una placa en la hipótesis de Reissner-Mindlin donde θi y dw/dxi no tienen necesariamente que coincidir.

Las hipótesis de Reissner-Mindlin son un conjunto de hipótesis cinemáticas sobre cómo se deforma una placa o lámina bajo flexión que permiten relacionar los desplazamientos con las deformaciones. Una vez obtenidas las deformaciones la aplicación rutinaria de las ecuaciones de la elasticidad permite encontrar las tensiones, y encontrar la ecuación que relaciona desplazamientos con las fuerzas exteriores.

Las hipótesis de Reissner-Mindlin para el cálculo elástico de placas y láminas son:

  1. El material de la placa es elástico lineal.
  2. El desplazamiento vertical para los puntos del plano medio no depende de z: uz(x, y, z) = w(x, y).
  3. Los puntos del plano medio sólo sufren desplazamiento vertical: ux(x, y,0) = 0, uy(x, y,0) = 0.
  4. La tensión perpendicular al plano medio se anula: σzz= 0.

Como consecuencia los desplazamientos horizontales sólo se dan fuera del plano medio y sólo se producen por giro del segmento perpendicular al plano medio. Como consecuencia de las hipótesis de Reissner-Mindlin los desplazamientos pueden escribirse como:

 \begin{cases}
    u_x(x,y,z) = -z\theta_x(x,y) \\
    u_y(x,y,z) = -z\theta_y(x,y) \\
    u_z(x,y,z) = w(x,y) 
\end{cases}


Hipótesis de Love-Kirchhoff[editar]

En las placas en que se desprecia la deformación por cortante, puede suponerse adecuadamente una hipótesis adicional conocida como hipótesis de Love-Kirchhoff. Esta hipótesis dice que:

5. \theta_x(x,y) = \frac{\partial w}{\partial x} \qquad
\theta_y(x,y) = \frac{\partial w}{\partial y}


Esta hipótesis es análoga a la hipótesis de Navier-Bernoulli para vigas. De hecho existe un paralelo entre los modelos de vigas y de placas. El modelo de placa de Reissner-Mindlin es el el equivalente de la viga de Timoshenko, mientras que el modelo de placa de Love-Kirchhoff es el equivalente de la viga de Euler-Bernoulli.

Las hipótesis de Reissner-Mindlin combinada con la hipótesis de Love-Kirchhoff proporcionan una hipótesis cinemática para los desplazamientos. A partir de esos desplazamientos pueden calcularse fácilmente las deformaciones para una placa delgada:


\varepsilon_{xx} = \frac{\partial u_x}{\partial x} = -z\frac{\partial^2 w}{\partial x^2} \qquad
\varepsilon_{yy} = \frac{\partial u_y}{\partial y} = -z\frac{\partial^2 w}{\partial y^2} \qquad
\varepsilon_{xy} = \frac{1}{2}\left ( \frac{\partial u_x}{\partial y}+\frac{\partial u_y}{\partial x} \right ) = -z\frac{\partial^2 w}{\partial x \partial y}


En función de esas deformaciones las tensiones se calculan trivialmente a partir de las ecuaciones de Lamé-Hooke que generalizan la ley de Hooke para sólidos deformables.

Ecuación de Lagrange para placas delgadas[editar]

Para una placa plana de espesor constante en la que sean válidas las hipótesis de Reissner-Mindlin y Love-Kircchoff el descenso vertical en cada punto bajo la acción de las cargas apoyadas sobre ella viene dada por:

(1)\Delta \Delta w(x,y) = \frac{q(x,y)}{D}

Donde w(x, y) es la flecha vertical o descenso vertical de la placa en el punto de coordenadas (x, y), q(x, y) es la carga por unidad de área en el mismo punto, el operador laplaciano se define, en coordenadas cartesianas, por la siguiente suma de operadores:

\Delta = \frac {\partial^2}{\partial x^2} +
\frac {\partial^2}{\partial y^2}, \qquad \qquad \Rightarrow \Delta\Delta w = 
\frac {\part^4 w}{\part x^4} + \frac {\part^4 w}{\part x^2y^2} + \frac {\part^4 w}{\part y^4}

Y finalmente la constante D es la rigidez flexional de placas y viene dada en función del espesor de la placa (h), el módulo de Young (E), el coeficiente de Poisson (ν):

D = \frac{Eh^3}{12(1-\nu^2)}

Es interesante notar que la ecuación (1) es el análogo de la ecuación de la elástica para vigas. Para placas de espesor no constante, análogamente al caso de la ecuación de la elástica para vigas, la flecha y la carga aplicada están relacionadas por la ecuación:

(2)\Delta \left (D \Delta w(x,y) \right ) = q(x,y)

Donde ahora la rigidez flexional D es función una D(x, y) que depende del punto concreto de placa.

La resolución de la ecuación (1) en general no es trivial y requiere tanto el uso de coordenadas adecuadas (para placas rectangulares se emplean coordenadas cartesianas, pero para placas circulares se emplean coordenadas cilíndricas) como la elección de algún método adecuado de integración. Entre los más sencillos están el método de Navier-Kirchhof[1] y el método de Levy[2] que se basan en series obtenidas mediante el método de separación de variables.

Cálculo de tensiones en placas delgadas[editar]

En una lámina sometida fundamentalmente a flexión en la que se desprecia la deformación por cortante, o lámina de Love-Kirchhof, los esfuerzos internos se caracterizan por dos momentos flectores m_x, m_y\; según dos direcciones mutuamente perpendiculares y un esfuerzo de torsión m_{xy}. Estos esfuerzos están directamente relacionados con la flecha vertical w(x, y) en cada punto por:

\begin{cases} 
m_x = -D\left[\cfrac{\part^2 w}{\part x^2}+ \nu \cfrac{\part^2 w}{\part y^2}\right] & 
m_{xy} = -D(1-\nu) \left[\cfrac{\part^2 w}{\part y\part x}\right]\\
m_y = -D\left[\nu \cfrac{\part^2 w}{\part x^2}+ \cfrac{\part^2 w}{\part y^2}\right] \end{cases}

Donde:

\nu\,, es el coeficiente de Poisson del material de la placa.
D = Eh^3/12(1-\nu^2)\;, es la rigidez en flexión de la placa, siendo:
E\; el módulo de Young del material de la placa, y h el espesor de la placa.

Las tensiones sobre una placa son directamente calculabes a partir de los esfuerzos anteriores:

\begin{cases} 
\sigma_{xx} = \cfrac{12z}{h^3}m_y(x,y) & \sigma_{xy} = \cfrac{12z}{h^3}(1-\nu)^2m_{xy}(x,y)\\
\sigma_{yy} = \cfrac{12z}{h^3}m_x(x,y) & \sigma_{xz} = \sigma_{yz} = \sigma_{zz} = 0 \end{cases}

Placas rectangulares[editar]

Para una placa rectangular de dimensiones a y b con carga uniforme q (por unidad de superfice) y simplemente apoyada en sus extremos la deflexión vertical \scriptstyle w(x,y) de la misma viene dada por:

w(x,y) = \frac{16q}{\pi^6 D}\sum_{m=1,3,5,\dots}^\infty\sum_{n=1,3,5,\dots}^\infty
\frac{1}{mn}\frac{\sin\left(\frac{m\pi x}{a}\right)\sin\left(\frac{n\pi y}{b}\right)}{\left[ \left( \frac{m}{a}\right)^2+ \left(\frac{n}{b}\right)^2 \right]^2}

con:

D=Eh^3/[12(1-\nu)]\,, rigidez flexional de placas
h << \min(a,b)\,, grosor de la placa.
E, \nu\,, módulo de Young y coeficiente de Poisson del material de la placa.

La anterior serie converge muy rápidamente por lo que se obtiene una muy buena aproximación tomando sólo los 3 o 6 primeros términos, además puede demostrarse que la flecha máxima cumple:

w_\max = \frac{16q}{\pi^6 D}\sum_{m=1,3,5,\dots}^\infty\sum_{n=1,3,5,\dots}^\infty
\frac{(-1)^{\frac{m+n}{2}-1}}{mn\left[ \left( \frac{m}{a}\right)^2+ \left(\frac{n}{b}\right)^2 \right]^2} \le \frac{16q}{\pi^6 D}\left( \frac{a^2b^2}{a^2 + b^2} \right)^2 =
\frac{192}{\pi^6}\frac{q a^4}{Eh^3}(1-\nu^2)\frac{\lambda^2}{(1+\lambda^2)^2}

donde \lambda = a/b\,.

Cálculo de láminas[editar]

Una lámina es un elemento estructural bidimensional curvado. Si las placas se tratan análogamente a las vigas rectas, las láminas son el análogo bidimensional de los arcos.

Usando coordenadas curvilíneas ortogonales sobre la superficie \scriptstyle (\alpha, \beta)\, se pueden escribir las ecuaciones de equilibrio para los esfuerzos internos para una lámina de Reisner-Mindlin como:[3] [4]

\begin{cases}
-\cfrac{U'_v}{UV}n_{uv}+\cfrac{V'_u}{UV}n_{vv}+\cfrac{q_u}{R_u}-\cfrac{1}{UV}
\left[(Vn_{uu})'_u + (Un_{vu})'_v \right] = p_u \\
+\cfrac{U'_v}{UV}n_{uu}-\cfrac{V'_u}{UV}n_{vu}+\cfrac{q_v}{R_v}-\cfrac{1}{UV}
\left[(Vn_{uv})'_u + (Un_{vv})'_v \right] = p_v  \\
-\cfrac{n_{uu}}{R_u} - \cfrac{n_{vv}}{R_v} -\cfrac{1}{UV}
\left[(Vq_u)'_u + (Uq_v)'_v \right] = p_\eta \\
-\cfrac{U'_v}{UV}m_{uv}+\cfrac{V'_u}{UV}m_{vv} + q_u -\cfrac{1}{UV}
\left[(Vm_{uu})'_u + (Um_{vu})'_v \right] = m_u \\
+\cfrac{U'_v}{UV}m_{uu}-\cfrac{V'_u}{UV}m_{vu}+ q_v -\cfrac{1}{UV}
\left[(Vm_{uv})'_u + (Um_{vv})'_v \right] = m_v
\end{cases}

Donde:

{'}_u,\ {'}_v, indican las derivadas parciales respecto a las coordenadas u, v.
U = \|{\mathbf{r}'}_u\| es el módulo del vector tangente asociado a la coordenada u.
V = \|{\mathbf{r}'}_v\| es el módulo del vector tangente asociado a la coordenada v.
R_u, R_v\, son los radios de curvatura según las direcciones de las líneas coordenadas.
p_u, p_v, p_\eta\, son las fuerzas por unidad de área en cada punto de la lámina.
m_u, m_v\, son las momentos por unidad de área en cada punto de la lámina.
n_{uu}, n_{uv}, n_{vu}, n_{vv}\, son los esfuerzos de membrana.
q_u, q_v\, son los esfuerzos cortantes de la placa.
m_{uu}, m_{vv}\, son los momentos flectores de la placa.
m_{uv}, m_{vu}\, son los momentos torsores de la placa.

Cúpula bajo su peso propio[editar]

Como ejemplo de las anteriores ecuaciones podemos considerar una cúpula en forma de casquete esférico sometida a su propio peso. Cada punto de la cúpula bidimensional se puede parametrizar mediante las coordenadas (u,v)\ [= (\theta, \phi)]\,:

\mathbf{r} = R_C(\sin\theta\cos\phi, \sin\theta\sin\phi, \cos\theta), \qquad
\mathbf{p} = (p_\theta, p_\phi, p_\eta) = -p(\sin\theta, 0, \cos\theta)

Con lo cual tenemos los factores geométricos siguientes:

U=\|\mathbf{r}_\theta\| =R_C, \qquad
V=\|\mathbf{r}_\phi\| = R_C\sin\theta

Y por tanto las ecuaciones anteriores quedarán reducidas a:

\begin{cases}
+n_{\phi\phi}\cos\theta + q_\theta - (n_{\theta\theta})'_\theta - n_{\theta\theta}\tan\theta = -R_Cp\sin\theta \\
-n_{\phi\theta}\cos\theta +\cfrac{q_\phi}{\sin\theta}- (n_{\theta\phi})'_\theta - n_{\theta\phi}\tan\theta = 0  \\
-n_{\theta\theta} - \cfrac{n_{\phi\phi}}{\sin\theta} -(q_\theta)'_\theta - q_\theta\tan\theta = -R_Cp\cos\theta \\
+m_{\phi\phi}\cos\theta + R_Cq_\theta -(m_{\theta\theta})'_\theta - m_{\theta\theta}\tan\theta = 0 \\
-m_{\phi\theta}\cos\theta + R_Cq_\phi -(m_{\theta\phi})'_\theta - m_{\theta\phi}\tan\theta = 0 \end{cases}

Lámina axisimétrica[editar]

El caso general de una lámina general requiere usar coordenadas curvilíneas generales para parametrizar su superficie, eso conduce a ecuaciones de gobierno que son ecuaciones en derivadas parciales cuya integración es complicada. Sin embargo muchos casos de interés involucran láminas con simetría axial o de revolución, con cargas que también respetan la simetría axial. En esos casos la geometría de la superficie puede parametrizarse mediante una coordenada (que da su "perfil" radial), y las ecuaciones de gobierno en ese caso involucran derivadas respecto a una única coordenada, y por tanto son un sistema de ecuaciones diferenciales ordinarias. Gracias a eso el comportamiento real puede estimarse mediante métodos clásicos, ya que resulta factible integrar en algunos casos las ecuaciones de gobierno. Esto contrata con el caso general para el cual no se conocen las soluciones analíticas de las ecuaciones de gobierno, por lo que en el caso general el comportamiento sólo puede investigarse buscando soluciones numéricas aproximadas a las ecuaciones de gobierno.

El caso más simple de teoría de láminas axisimétricas es la teoría estática especial de Cosserat que describe el comportamiento de placas axisimétricas susceptibles de sufrir flexión en su superficie media, extensión de la misma y cortante en el espesor.

Referencia[editar]

Bibliografía[editar]

  • Antman, Stuart S. (1995). «X. Axisymmetric Equilibria of Cosserat Shells». Nonlinear Problems of Elasticity (libro). Applied Mathematical Sciences (en inglés) 107. Nueva York: Springer-Verlag. ISBN 0-387-94199-1. 
  • Langhaar, H. L. (1962). Energy Methods in Applied Mechanics. Wiley. ISBN 978-0-89464-364-4. .
  • Washizu, K. (1974). Variational methods in Elasticity and Plasticity. Pergamon Press. ISBN 978-0-08-026723-4. 

Enlaces externos[editar]