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.

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

(1)${displaystyle {egin{cases}{cfrac {partial sigma _{xx}}{partial x}}+{cfrac {partial sigma _{xy}}{partial y}}+{cfrac {partial sigma _{xz}}{partial z}}+b_{x}=0\{cfrac {partial sigma _{yx}}{partial x}}+{cfrac {partial sigma _{yy}}{partial y}}+{cfrac {partial sigma _{yz}}{partial z}}+b_{y}=0\{cfrac {partial sigma _{zx}}{partial x}}+{cfrac {partial sigma _{zy}}{partial y}}+{cfrac {partial sigma _{zz}}{partial z}}+b_{z}=0end{cases}}Rightarrow quad {oldsymbol { abla }}cdot {oldsymbol {sigma }}+mathbf {b} =0}$

Donde ${displaystyle mathbf {b} =(b_{x},b_{y},b_{z})}$ son las fuerzas por unidad de volumen que actúan en el interior del cuerpo y ${displaystyle {oldsymbol {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)${displaystyle {egin{cases}sigma _{xx}=lambda (varepsilon _{xx}+varepsilon _{yy}+varepsilon _{zz})+2mu varepsilon _{xx}&sigma _{xy}=2mu varepsilon _{xy}\sigma _{yy}=lambda (varepsilon _{xx}+varepsilon _{yy}+varepsilon _{zz})+2mu varepsilon _{yy}&sigma _{yz}=2mu varepsilon _{yz}\sigma _{zz}=lambda (varepsilon _{xx}+varepsilon _{yy}+varepsilon _{zz})+2mu varepsilon _{zz}&sigma _{zx}=2mu varepsilon _{zx}end{cases}}Rightarrow quad sigma _{ij}=lambda delta _{ij}varepsilon _{V}+2mu varepsilon _{ij}}$

Donde se han usado los coeficientes de Lamé (λ, μ) y la deformación volumétrica ${displaystyle 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)${displaystyle {egin{cases}varepsilon _{xx}={cfrac {partial u_{x}}{partial x}}&varepsilon _{xy}={cfrac {1}{2}}left({cfrac {partial u_{x}}{partial y}}+{cfrac {partial u_{y}}{partial x}} ight)\varepsilon _{yy}={cfrac {partial u_{y}}{partial y}}&varepsilon _{yz}={cfrac {1}{2}}left({cfrac {partial u_{y}}{partial z}}+{cfrac {partial u_{z}}{partial y}} ight)\varepsilon _{zz}={cfrac {partial u_{z}}{partial z}}&varepsilon _{zx}={cfrac {1}{2}}left({cfrac {partial u_{z}}{partial x}}+{cfrac {partial u_{x}}{partial z}} ight)end{cases}}Rightarrow quad varepsilon _{ij}={cfrac {1}{2}}left({cfrac {partial u_{i}}{partial x_{j}}}+{cfrac {partial u_{j}}{partial x_{i}}} ight)}$

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 (u_x, u_y, u_z), 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)${displaystyle sigma _{ij}=lambda delta _{ij}left({frac {partial u_{x}}{partial x}}+{frac {partial u_{y}}{partial y}}+{frac {partial u_{z}}{partial z}} ight)+mu left({frac {partial u_{i}}{partial x_{j}}}+{frac {partial u_{i}}{partial x_{j}}} ight)}$

Donde ${displaystyle 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 parciales que relacionan las fuerzas de volumen (b_x, b_y, b_z) con los desplazamientos:^[1]​

(5a)${displaystyle {egin{cases}b_{x}+(lambda +mu ){frac {partial }{partial x}}( abla cdot mathbf {u} )+mu Delta u_{x}=0\b_{y}+(lambda +mu ){frac {partial }{partial y}}( abla cdot mathbf {u} )+mu Delta u_{y}=0\b_{z}+(lambda +mu ){frac {partial }{partial z}}( abla cdot mathbf {u} )+mu Delta u_{z}=0end{cases}}}$

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

(5b)${displaystyle mathbf {b} +(lambda +mu ) abla ( abla cdot mathbf {u} )+mu Delta mathbf {u} =0}$

El problema elástico lineal también puede formularse exclusivamente en términos de tensiones, aunque en ese caso la solución es única salvo un movimiento de sólido rígido para los desplazamientos. Si se reescriben las ecuaciones de compatibilidad en términos de las tensiones y en ellas se substiyen las ecuaciones de equilibrio se llegan a las ecuaciones de Michell:

(6a)${displaystyle Delta sigma _{ij}+{cfrac {1}{1+ u }}{frac {partial ^{2}sigma _{V}}{partial x_{i}partial x_{j}}}={frac { u }{1+ u }}delta _{ij}left({oldsymbol { abla }}cdot mathbf {b} ight)-left({frac {partial f_{j}}{partial x_{i}}}+{frac {partial f_{i}}{partial x_{j}}} ight)}$

Donde:

Un caso particular de las anteriores ecuaciones, cuando las fuerzas volumétricas son constantes, son las ecuaciones de Beltrami, en las que el segund miembro se anula:

(6b)${displaystyle {egin{cases}Delta sigma _{xx}+{cfrac {1}{1+ u }}{cfrac {partial ^{2}sigma _{V}}{partial x^{2}}}=0&Delta sigma _{yz}+{cfrac {1}{1+ u }}{cfrac {partial ^{2}sigma _{V}}{partial ypartial z}}=0\Delta sigma _{yy}+{cfrac {1}{1+ u }}{cfrac {partial ^{2}sigma _{V}}{partial y^{2}}}=0&Delta sigma _{xz}+{cfrac {1}{1+ u }}{cfrac {partial ^{2}sigma _{V}}{partial xpartial z}}=0\Delta sigma _{zz}+{cfrac {1}{1+ u }}{cfrac {partial ^{2}sigma _{V}}{partial z^{2}}}=0&Delta sigma _{xy}+{cfrac {1}{1+ u }}{cfrac {partial ^{2}sigma _{V}}{partial xpartial y}}=0end{cases}}}$

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:

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

Y se calcula la energía asociada al campo de desplazamientos:

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

para cualquier punto del cuerpo elástico.

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:

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

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.

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 ${displaystyle X^{J},}$ y las coordenadas sobre el cuerpo deformado como ${displaystyle x^{j},}$ de tal manera que la deformación puede representarse por el difeomorfismo:

${displaystyle phi :Omega o 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 original no deformado en notación tensorial resultan ser:

${displaystyle {oldsymbol { abla cdot sigma }}+mathbf {b} =0quad Rightarrow quad {frac {1}{sqrt {|g|}}}{frac {partial }{partial x^{k}}}left({sqrt {|g|}}sigma ^{kj} ight)={frac {partial sigma ^{kj}}{partial x^{k}}}+sigma ^{kj}{frac {partial }{partial x^{k}}}left(ln {sqrt {|g|}} ight)+b^{j}=0}$

Donde ${displaystyle |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 ${displaystyle Sigma _{R}=(S_{IJ}),}$:

(*)${displaystyle S_{IJ}=|J|{frac {partial x^{i}}{partial X^{I}}}{frac {partial x^{j}}{partial X^{J}}}sigma _{ij},qquad |J|=left|{frac {partial (X^{1},...,X^{3})}{partial (x^{1},...,x^{3})}} ight|}$

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

${displaystyle {frac {partial X^{I}}{partial x^{k}}}{frac {partial }{partial X^{I}}}left({frac {partial x^{k}}{partial X^{K}}}{frac {partial x^{j}}{partial X^{J}}}S^{KJ} ight)+{frac {1}{|J|}}{frac {partial x^{j}}{partial X^{J}}}S^{KJ}{frac {partial ln {sqrt {|g|}}}{partial X^{K}}}+b^{j}=0}$

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).

En cambio para un problema de desplazamiento puro, donde se prescriben desplazamientos dados sobre toda la frontera de cuerpo elástico está garantizada la existencia de solución. También para problemas desplazamiento-tracción en que para cada punto de la frontera del sólido está prescrito o bien su desplazamiento o bien la fuerza puntual que actuará sobre dicho punto. En todos estos casos J. Ball demostró la existencia de solución y probó que dicha solución es un mínimo del funcional de energía, bajo ciertas condiciones técnicas sobre la ecuación constitutiva del material elástico.

---