Problema elástico

Problema elástico

Problema elástico

El problema elástico es el problema físico-matemático de encontrar los desplazamientos y las tensiones en un sólido deformable elástico, partiendo de la forma original del sólido, de las fuerzas actuantes sobre el mismo y de los desplazamientos impuestos de algunos puntos de la superficie del sólido.

El problema elástico lineal es un tipo particular de problema elástico en que tanto la ecuación constitutiva, como la relación entre deformaciones y desplazamientos viene dada por ecuaciones lineales. En este caso el problema generalmente se puede reducir a un problema mixto de Dirichlet-Von Neumann sobre un dominio del espacio que coincide con la forma del cuerpo elástico antes de la deformación.

Contenido

El problema elástico lineal

El problema elástico lineal está formado por un sistema de 15 ecuaciones diferenciales lineales, más sus condiciones de contorno. Este sistema de ecuaciones está formado por las tres ecuaciones de equilibrio que expresan que la suma de fuerzas sobre cualquier punto del cuerpo es cero:

(1) \begin{cases}
\cfrac{\part\sigma_{xx}}{\part x}+ \cfrac{\part\sigma_{xy}}{\part y}+ \cfrac{\part\sigma_{xz}}{\part z} + b_x = 0\\
\cfrac{\part\sigma_{yx}}{\part x}+ \cfrac{\part\sigma_{yy}}{\part y}+ \cfrac{\part\sigma_{yz}}{\part z} + b_y = 0\\
\cfrac{\part\sigma_{zx}}{\part x}+ \cfrac{\part\sigma_{zy}}{\part y}+ \cfrac{\part\sigma_{zz}}{\part z} + b_z = 0 \end{cases} \Rightarrow \quad \boldsymbol\nabla\cdot\boldsymbol\sigma+ \mathbf{b} = 0

Donde \mathbf{b} = (b_x,b_y,b_z) son las fuerzas por unidad de volumen que actúan en el interior del cuerpo y \boldsymbol\sigma = [\sigma_{ij}] es el tensor tensión de cuerpo. Además de las tres ecuaciones anteriores interviene las ecuaciones constitutivas que representan seis ecuaciones escalares más:

(2) \begin{cases}
\sigma_{xx} = \lambda(\varepsilon_{xx}+\varepsilon_{yy}+\varepsilon_{zz})+2\mu\varepsilon_{xx} & 
\sigma_{xy} = 2\mu \varepsilon_{xy}\\
\sigma_{yy} = \lambda(\varepsilon_{xx}+\varepsilon_{yy}+\varepsilon_{zz})+2\mu\varepsilon_{yy} & 
\sigma_{yz} = 2\mu \varepsilon_{yz}\\
\sigma_{zz} = \lambda(\varepsilon_{xx}+\varepsilon_{yy}+\varepsilon_{zz})+2\mu\varepsilon_{zz} & 
\sigma_{zx} = 2\mu \varepsilon_{zx} \end{cases}
\Rightarrow \quad \sigma_{ij} = \lambda\delta_{ij}\varepsilon_V + 2\mu\varepsilon_{ij}

Donde se han usado los coeficientes de Lamé (λ, μ) y la deformación volumétrica \varepsilon_V = \varepsilon_{xx}+\varepsilon_{yy}+\varepsilon_{zz}. Finalmente se requiere una relación entre desplazamientos y deformaciones, que en última instancia permite relacionar las fuerzas aplicadas con los desplazamientos sufridos por el cuerpo. En teoría de la elasticidad lineal estas relaciones vienen dadas por:

(3) \begin{cases}
\varepsilon_{xx} = \cfrac{\part u_x}{\part x} & 
\varepsilon_{xy}=\cfrac{1}{2}\left( \cfrac{\part u_x}{\part y}+\cfrac{\part u_y}{\part x} \right) \\
\varepsilon_{yy} = \cfrac{\part u_y}{\part y} & 
\varepsilon_{yz}=\cfrac{1}{2}\left( \cfrac{\part u_y}{\part z}+\cfrac{\part u_z}{\part y} \right) \\
\varepsilon_{zz} = \cfrac{\part u_z}{\part z} & 
\varepsilon_{zx}=\cfrac{1}{2}\left( \cfrac{\part u_z}{\part x}+\cfrac{\part u_x}{\part z} \right) \end{cases} \Rightarrow \quad \varepsilon_{ij}=\cfrac{1}{2}\left( \cfrac{\part u_i}{\part x_j}+ \cfrac{\part u_j}{\part x_i} \right)

Formulación de Navier en desplazamientos

Esta formulación evita el problema de integrar las ecuaciones de compatibilidad ya que plantea directamente el problema en función de los tres campos de desplazamientos incógnita (ux, uy, uz), con lo cual las ecuaciones de compatibilidad se cumplen automáticamente y no son necesarias. Matemáticamente basta substituir las ecuaciones (3) en las ecuaciones (2) y substituir nuevamente el resultado en las ecuaciones (1). Este esquema, reduce enormemente el tamaño del sistema de ecuaciones diferenciales que hay que resolver. Para un sólido elástico lineal e isótropo las tensiones en función de los campos de desplazamiento vienen dados por:

(4) \sigma_{ij} = \lambda\delta_{ij}\left(\frac{\part u_x}{\part x}+ \frac{\part u_y}{\part y}+ \frac{\part u_z}{\part z} \right) + \mu \left(\frac{\part u_i}{\part x_j}+\frac{\part u_i}{\part x_j} \right)

Donde \lambda, \mu\; son los llamados coeficientes de Lamé que caracterizan el comportamiento elástico del material. Si las expresiones (4) se introducen en las ecuaciones de equilibrio (1) para puntos interiores del sólido se llega a un sistema de tres ecuaciones en derivadas paraciales que relacionan las fuerzas de volumen (bx, by, bz) con los desplazamientos:[1]

(5a) \begin{cases}
 b_x + (\lambda+\mu)\frac{\part}{\part x}(\nabla\cdot\mathbf{u})+\mu\Delta u_x = 0\\
 b_y + (\lambda+\mu)\frac{\part}{\part y}(\nabla\cdot\mathbf{u})+\mu\Delta u_y = 0\\
 b_z + (\lambda+\mu)\frac{\part}{\part z}(\nabla\cdot\mathbf{u})+\mu\Delta u_z = 0
 \end{cases}

Las tres ecuaciones anteriores pueden reunirse en una única ecuación vectorial como:

(5b) \mathbf{b} + (\lambda+\mu)\nabla(\nabla\cdot\mathbf{u})+\mu\Delta \mathbf{u} = 0

Formulación de Michell-Beltrami en tensiones

Formulación de Patnaik en deformaciones

Problema elástico lineal dependiente del tiempo

Unicidad del problema elástico lineal

El problema elástico lineal definido por sus 15 ecuaciones de gobierno dada por (1), (2) y (3) más las condiciones de contorno tiene solución única, tanto en su versión dependiente del tiempo con en su versión independiente del tiempo, siempre y cuando el tensor de constantes elásticas sea definido positivo, cosa que sucede para todos los materiales conocidos.

La demostración matemática de este hecho se realiza suponiendo que existen dos soluciones:

\mathbf{u}^{(1)}(x,y,z), \qquad \mathbf{u}^{(2)}(x,y,z)

que satisfacen las ecuaciones del problema y las condiciones de contorno. A continuación se define el funcional de potencia mecánica:

\dot{E}(u):=\frac{1}{2}\int_V \left( \rho||\dot\mathbf{u}||^2 +\sum_{i,j} \sigma_{ij}(\mathbf{u})\dot\epsilon_{ij}(\mathbf{u}) \right)dV

Y se calcula la energía asociada al campode desplazamientos:

\tilde{\mathbf{u}} := \mathbf{u}^{(1)} - \mathbf{u}^{(2)}

Se puede probar fácilmente que \dot{E}(\tilde{\mathbf{u}}) = 0 lo cual dado el carácter definido postivo del tensor de constantes elásticas implica necesariamente que:

\tilde{\mathbf{u}} =0 \Rightarrow \qquad \mathbf{u}^{(1)} = \mathbf{u}^{(2)}

para cualquier punto del cuerpo elástico.

Casos particulares

Si bien no se conoce un procedimiento general para hallar una solución al problema elástico en su forma más general. Existen numerosas soluciones o esquemas de solución útiles cuando el problema elástico aparece en ciertas geometrías o formas simplificadas. Típicamente el problema elástico resulta más sencillo cuando:

  • Se tiene un estado elástico de elasticidad plana, es decir, cuando en el cálculo de las tensiones y/o las deformaciones puede reducirse el problema elástico a un cierto tipo de problema bidimensional.
  • Existe simetría axial de algún tipo.
  • En cierto tipo de sólidos elásticos semi-infinitos, lo cual permite dar resupuestas aproximadas a muchos problemas de ingeniería geotécnica cuando el terreno se modeliza como un sólido elástico semi-infinito.

El problema elástico no-lineal

El problema elástico no lineal plantea dificultades adicionales serias respecto al problema no lineal entre ellas:

  1. El sistema de ecuaciones es no lineal, en este caso tanto las ecuaciones de compatibilidad como las ecuaciones que relacionan desplazamientos y deformaciones son no lineales.
  2. La forma geométrica del cuerpo antes y después de la deformación no coinciden, lo cual lleva a la dificultad adicional de describir el dominio sobre el que están definidas las ecuaciones diferenciales del problema elástico no lineal.

Estas dos dificultades generalmente son inseparables: sólo las tres ecuaciones de equilibrio son lineales si se estudia el problema sobre la geometría del cuerpo una vez deformada, aunque esta forma deformada a priori no es conocida tal como expresa la segunda dificultad.

Tensor tensión de Piola-Kirchhoff

Una posibilidad es tratar de resolver el problema elástico teniendo en cuenta que existe un difeomorfismo entre la forma del cuerpo una vez deformado y la forma del cuerpo antes de la deformación. Siguiendo la convención usual se designan las coordenadas sobre el cuerpo antes de deformar mediante X^J\, y las coordenadas sobre el cuerpo deformado como x^j\, de tal manera que la deformación puede representarse por el difeomorfismo:

\phi:\Omega \to \Omega', \qquad (x^1,x^2,x^3)= \phi(X^1,X^2,X^3)

Mediante este difeomorfismo parte puede tratar de escribirse las ecuaciones del problema en lugar sobre el dominio ocupado por el cuerpo una vez deformado, que a priori es desconocido, sobre el dominio antes de la deformación. Las ecuaciones de equilibrio sobre el dominio orginal no deformado en notación tensorial resultan ser:

\boldsymbol{\nabla\cdot\sigma} = \mathbf{b} \quad \Rightarrow \quad 
\frac{1}{\sqrt{|g|}} \frac{\part }{\part x^k} \left(\sqrt{|g|} \sigma^{kj} \right) =
\frac{\part \sigma^{kj}}{\part x^k} +
\sigma^{kj}\frac{\part }{\part x^k} \left(\ln \sqrt{|g|} \right)= b^j

Donde |g|\, es el determinante del tensor métrico que coincide con el cuadrado del jacobiano respecto a las coordenadas cartesianas. Para poder escribir las ecuaciones de equilibrio sobre el dominio de referencia sin deformar podemos usar el segundo tensor de Piola-Kirchhoff \Sigma_R = (S_{IJ})\,:

(*) S_{IJ}= |J| \frac{\part x^i}{\part X^I} \frac{\part x^j}{\part X^J} \sigma_{ij}, \qquad  
|J| = \left| \frac{\part(X^1,...,X^3)}{\part(x^1,...,x^3)}\right |

Reescribiendo (*) esa encuación gracias al tensor de Piola-Kirchhoff tenemos que las ecuaciones de equilibrio toman la forma:

\frac{\part X^I}{\part x^k} \frac{\part}{\part X^I}
\left(\frac{\part x^k}{\part X^K} \frac{\part x^j}{\part X^J} S^{KJ}\right) +
\frac{1}{|J|} \frac{\part x^j}{\part X^J} S^{KJ} \frac{\part \ln \sqrt{|g|}}{\part X^K} = b^j

Existencia y unicidad del problema elástico no lineal

Tal como muestran numerosos ejemplos de no linealidad geométrica e inestabilidad elástica como la abolladura o la inestabilidad de arcos, en general para un valor dado de las cargas existen varias configuraciones deformadas posibles, es decir, en esos casos no existe una solución única del problema elástico compatible con las condiciones de contorno y las fuerzas aplicadas.

Es más, puede incluso suceder que para ciertas configuraciones ni siquiera existe equilibrio compatible con ciertos valores de la deformación, como sucede en la inestabilidad de arcos. Esta dificultad se resuelve si considera el problema elástico general dependiente del tiempo, entonces siempre existe solución (aunque ésta puede no ser única si la solución corresponde a una configuración estática de equilibrio).

Referencias

  1. Ortiz Berrocal, 1998, p. 152.
  • Ortiz Berrocal, L., Elasticidad, McGraw-Hill, 1998, ISBN 84-481-2046-9.
  • Timoshenko, S.P. y Godier J.N., Theory of elasticity, McGraw-Hill, 1951.
  • S. N. Patnaik y D. A. Hopslins: "Stress Formulation in Three Dimensional Elasticity", NASA-TP-2001-210515, Septiembre 2001 [1].
  • S. N. Patnaik y S. Pai: "Boundary compatibility condition and rotation in elasticity", International Journal of Physical Sciences, Vol. 1 (2), pp. 081-084, Octubre 2006.
Obtenido de "Problema el%C3%A1stico"

Wikimedia foundation. 2010.

Игры ⚽ Поможем сделать НИР

Mira otros diccionarios:

  • Muelle elástico — Saltar a navegación, búsqueda Para otros usos de este término, véase Muelle. «Resorte» redirige aquí. Para otras acepciones, véase Resorte (desambiguación). Muelles de tracción. Se conoce como …   Wikipedia Español

  • Elasticidad (mecánica de sólidos) — Una varilla elástica vibrando, es un ejemplo de sistema donde la energía potencial elástica se transforma en energía cinética y viceversa. En física e ingeniería, el término elasticidad designa la propiedad mecánica de ciertos materiales de… …   Wikipedia Español

  • Formulación débil de una ecuación diferencial — La formulación débil (o formulación variacional) de un problema definido mediante ecuaciones diferenciales es una forma alternativa en que dichas ecuaciones se escriben en forma integral, dando lugar a ecuaciones tratables mediante los métodos… …   Wikipedia Español

  • Elasticidad plana — La elasticidad plana se refiere al estudio de soluciones particulares del problema elástico general y al estudio del conjunto de aplicaciones técnicas en que aparecen dichos estados elásticos de tensión deformación reducibles a problemas planos o …   Wikipedia Español

  • Mecánica de sólidos deformables — La mecánica de los sólidos deformables estudia el comportamiento de los cuerpos sólidos deformables ante diferentes tipos de situaciones como la aplicación de cargas o efectos térmicos. Estos comportamientos, más complejos que el de los sólidos… …   Wikipedia Español

  • Mecánica de la fractura — Dentro de las ramas de la mecánica, la mecánica de la fractura está teniendo un enorme auge en la actualidad al estudiar los procesos catastróficos de rotura en estructuras: aviones, puentes... Mecánica de fractura ó Mecánica de la fractura es… …   Wikipedia Español

  • Método de los elementos finitos — Solución de MEF en 2D para una configuración de un magnetostato, (las líneas muestran la dirección de la densidad de flujo calculada, y el color, su magnitud) …   Wikipedia Español

  • Resistencia de materiales — Saltar a navegación, búsqueda Para otros usos de este término, véase Resistencia. La resistencia de materiales clásica es una disciplina de la ingeniería mecánica y la ingeniería estructural que estudia los sólidos deformables mediante modelos… …   Wikipedia Español

  • Lagrangiano — En física, un lagrangiano es una función matemática a partir de la cual se pueden obtener la evolución temporal, las leyes de conservación y otras propiedades importantes de un sistema físico. De hecho, en física moderna el lagrangiano se… …   Wikipedia Español

  • Ecuación de compatibilidad — Saltar a navegación, búsqueda Una ecuación de compatibilidad es una ecuación adicional a un problema mecánico de equilibrio necesaria para asegurar que la solución buscada es compatible con las condiciones de contorno o para poder asegurar la… …   Wikipedia Español

Compartir el artículo y extractos

Link directo
Do a right-click on the link above
and select “Copy Link”