La resistencia de materiales clásica es una disciplina de la ingeniería mecánica,

la ingeniería estructural, la ingeniería civil y la ingeniería de materiales que

estudia la mecánica de sólidos deformables mediante modelos simplificados. La
resistencia de un elemento se define como su capacidad para resistir esfuerzos y
fuerzas aplicadas sin romperse, adquirir deformaciones permanentes o deteriorarse
de algún modo.

Un modelo de resistencia de materiales establece una relación entre las fuerzas

aplicadas, también llamadas cargas o acciones, y los esfuerzos y desplazamientos
inducidos por ellas. Generalmente las simplificaciones geométricas y las
restricciones impuestas sobre el modo de aplicación de las cargas hacen que el
campo de deformaciones y tensiones sean sencillos de calcular.

Para el diseño mecánico de elementos con geometrías complicadas la resistencia de

materiales suele ser abundante y es necesario usar técnicas basadas en la teoría de
la elasticidad o la mecánica de sólidos deformables más generales. Esos problemas
planteados en términos de tensiones y deformaciones pueden entonces ser resueltos
de forma muy aproximada con métodos numéricos como el análisis por elementos
finitos.

Índice
1 Enfoque de la resistencia de materiales
2 Hipótesis cinemática
2.1 Hipótesis cinemática en elementos lineales
2.2 Hipótesis cinemática en elementos superficiales
3 Ecuación constitutiva
4 Ecuaciones de equivalencia
4.1 Elementos lineales
4.2 Elementos bidimensionales
5 Ecuaciones de equilibrio
5.1 Ecuaciones de equilibrio en elementos lineales rectos
5.2 Ecuaciones de equilibrio en elementos planos bidimensionales
6 Relación entre esfuerzos y tensiones
6.1 Elementos lineales o unidimensionales
6.2 Elementos superficiales o bidimensionales
7 Relación entre esfuerzos y desplazamientos
8 Véase también
9 Referencias
9.1 Bibliografía
9.2 Enlaces externos
Enfoque de la resistencia de materiales
La teoría de sólidos deformables requiere generalmente trabajar con tensiones y
deformaciones. Estas magnitudes vienen dadas por campos tensoriales definidos sobre
dominios tridimensionales que satisfacen complicadas ecuaciones diferenciales.

Sin embargo, para ciertas geometrías aproximadamente unidimensionales (vigas,

pilares, celosías, arcos, etc.) o bidimensionales (placas y láminas, membranas,
etc.) el estudio puede simplificarse y se pueden analizar mediante el cálculo de
esfuerzos internos definidos sobre una línea o una superficie en lugar de tensiones
definidas sobre un dominio tridimensional. Además las deformaciones pueden
determinarse con los esfuerzos internos a través de cierta hipótesis cinemática. En
resumen, para esas geometrías todo el estudio puede reducirse al estudio de
magnitudes alternativas a deformaciones y tensiones.

El esquema teórico de un análisis de resistencia de materiales comprende:

La hipótesis cinemática establece cómo serán las deformaciones o el campo de

desplazamientos para un determinado tipo de elementos bajo cierto tipo de
solicitudes. Para piezas prismáticas las hipótesis más comunes son la hipótesis de
Bernouilli-Navier para la flexión y la hipótesis de Saint-Venant para la torsión.
La ecuación constitutiva, que establece una relación entre las deformaciones o
desplazamientos deducibles de la hipótesis cinemática y las tensiones asociadas.
Estas ecuaciones son casos particulares de las ecuaciones de Lamé-Hooke.
Las ecuaciones de equivalencia son ecuaciones en forma de integral que relacionan
las tensiones con los esfuerzos internos.
Las ecuaciones de equilibrio relacionan los esfuerzos internos con las fuerzas
exteriores.
En las aplicaciones prácticas el análisis es sencillo. Se construye un esquema
ideal de cálculo formado por elementos unidimensionales o bidimensionales, y se
aplican fórmulas preestablecidas en base al tipo de solicitación que presentan los
elementos. Esas fórmulas preestablecidas que no necesitan ser deducidas para cada
caso, se basan en el esquema de cuatro puntos anteriores. Más concretamente la
resolución práctica de un problema de resistencia de materiales sigue los
siguientes pasos:

Cálculo de esfuerzos, se plantean las ecuaciones de equilibrio y ecuaciones de

compatibilidad que sean necesarias para encontrar los esfuerzos internos en función
de las fuerzas aplicadas.
Análisis resistente, se calculan las tensiones a partir de los esfuerzos internos.
La relación entre tensiones y deformaciones depende del tipo de solicitación y de
la hipótesis cinemática asociada: flexión de Bernouilli, flexión de Timoshenko,
flexión esviada, tracción, pandeo, torsión de Coulomb, teoría de Collignon para
tensiones cortantes, etc.
Análisis de rigidez, se calculan los desplazamientos máximos a partir de las
fuerzas aplicadas o los esfuerzos internos. Para ello puede recurrirse directamente
a la forma de la hipótesis cinemática o bien a la ecuación de la curva elástica,
las fórmulas vectoriales de Navier-Bresse o los teoremas de Castigliano.
Hipótesis cinemática
La hipótesis cinemática es una especificación matemática de los desplazamientos de
un sólido deformable que permite calcular las deformaciones en función de un
conjunto de parámetros incógnita.

El concepto se usa especialmente en el cálculo de elementos lineales (por ejemplo,

vigas) y elementos bidimensionales, donde gracias a la hipótesis cinemática se
pueden obtener relaciones funcionales más simples. Así pues, gracias a la hipótesis
cinemática se pueden relacionar los desplazamientos en cualquier punto del sólido
deformable de un dominio tridimensional con los desplazamientos especificados sobre
un conjunto unidimensional o bidimensional.

Hipótesis cinemática en elementos lineales

La resistencia de materiales propone para elementos lineales o prismas mecánicos,
como las vigas y pilares, en las que el desplazamiento de cualquier punto se puede
calcular a partir de desplazamientos y giros especificados sobre el eje
baricéntrico. Eso significa que por ejemplo para calcular una viga en lugar de
espeficar los desplazamientos de cualquier punto en función de tres coordenadas,
podemos expresarlos como función de una sola coordenada sobre el eje baricéntrico,
lo cual conduce a sistemas de ecuaciones diferenciales relativamente simples.
Existen diversos tipos de hipótesis cinemáticas según el tipo de solicitación de la
viga o elemento unidimensional:

La hipótesis de Navier-Bernouilli, que se usa para elementos lineales alargados

sometidos a flexión cuando las deformaciones por cortante resultan pequeñas.
La hipótesis de Timoshenko, que se usa para los elementos lineales sometidos a
flexión en un caso totalmente general ya que no se desprecia la deformación por
cortante.
La hipótesis de Saint-Venant para la extensión, usada en piezas con esfuerzo normal
para zonas de la viga alejadas de la zona de aplicación de las cargas.
La hipótesis de Saint-Venant para la torsión se usa para piezas prismáticas
sometidas a torsión y en piezas con rigidez torsional grande.
La hipótesis de Coulomb se usa para piezas prismáticas sometidas a torsión y en
piezas con rigidez torsional grande y sección circular o tubular. Esta hipótesis
constituye una especialización del caso anterior.
Hipótesis cinemática en elementos superficiales
Para placas y láminas sometidas a flexión se usan dos hipótesis, que se pueden
poner en correspondencia con las hipótesis de vigas:

hipótesis de Love-Kirchhoff
hipótesis de Reissner-Mindlin
Ecuación constitutiva
Las ecuaciones constitutivas de la resistencia de materiales son las que explicitan
el comportamiento del material, generalmente se toman como ecuaciones constitutivas
las ecuaciones de Lamé-Hooke de la elasticidad lineal. Estas ecuaciones pueden ser
especializadas para elementos lineales y superficiales. Para elementos lineales en
el cálculo de las secciones, las tensiones sobre cualquier punto (y,z) de la
sección puedan escribirse en función de las deformaciones como:

{\displaystyle \sigma (y,z)=E\ \varepsilon (y,z)}{\displaystyle \sigma (y,z)=E\

\varepsilon (y,z)}
{\displaystyle {\begin{cases}\sigma _{xx}=\sigma &\varepsilon
_{xx}=\varepsilon \\\sigma _{yy}=0&\varepsilon _{yy}=-\nu \varepsilon \\\sigma
_{zz}=0&\varepsilon _{zz}=-\nu \varepsilon \end{cases}}}{\displaystyle
{\begin{cases}\sigma _{xx}=\sigma &\varepsilon _{xx}=\varepsilon \\\sigma
_{yy}=0&\varepsilon _{yy}=-\nu \varepsilon \\\sigma _{zz}=0&\varepsilon
_{zz}=-\nu \varepsilon \end{cases}}}

En cambio, para elementos superficiales sometidos predominantemente a flexión como

las placas la especialización de las ecuaciones de Hooke es:

{\displaystyle {\begin{bmatrix}\sigma _{xx}\\\sigma _{yy}\\\sigma

_{zz}\end{bmatrix}}={\frac {E}{1-\nu ^{2}}}{\begin{bmatrix}1&\nu &0\\\nu
&1&0\\0&0&1-\nu \end{bmatrix}}{\begin{bmatrix}\varepsilon _{xx}\\\varepsilon
_{yy}\\\varepsilon _{zz}\end{bmatrix}}}{\displaystyle {\begin{bmatrix}\sigma
_{xx}\\\sigma _{yy}\\\sigma _{zz}\end{bmatrix}}={\frac {E}{1-\nu ^{2}}}
{\begin{bmatrix}1&\nu &0\\\nu &1&0\\0&0&1-\nu \end{bmatrix}}
{\begin{bmatrix}\varepsilon _{xx}\\\varepsilon _{yy}\\\varepsilon
_{zz}\end{bmatrix}}}

Además de ecuaciones constitutivas elásticas, en el cálculo estructural varias

normativas recogen métodos de cálculo plástico donde se usan ecuaciones
constitutivas de plasticidad.

Ecuaciones de equivalencia
Las ecuaciones de equivalencia expresan los esfuerzos resultantes a partir de la
distribución de tensiones. Gracias a ese cambio es posible escribir ecuaciones de
equilibrio que relacionen directamente las fuerzas aplicadas con los esfuerzos
internos.

Elementos lineales
En elementos lineales rectos las coordenadas cartesianas para representar la
geometría y expresar tensiones y esfuerzos, se escogen normalmente con el eje X
paralelo al eje baricéntrico de la pieza, y los ejes Y y Z coincidiendo con las
direcciones principales de inercia. En ese sistema de coordenadas la relación entre
esfuerzo normal (Nx), esfuerzos cortantes (Vy, Vz), el momento torsor (Mx) y los
momentos flectores (My, Mz) es:

{\displaystyle {\begin{matrix}N_{x}=\int _{\Sigma }\sigma _{x}dydz&V_{y}=\int

_{\Sigma }\tau _{xy}dydz&V_{z}=\int _{\Sigma }\tau _{xz}dydz\\M_{x}=\int _{\Sigma }
(-\tau _{xy}z+\tau _{xz}y)dydz&M_{y}=\int _{\Sigma }z\sigma _{xx}dydz&M_{z}=\int
_{\Sigma }-y\sigma _{xx}dydz\end{matrix}}}{\displaystyle {\begin{matrix}N_{x}=\int
_{\Sigma }\sigma _{x}dydz&V_{y}=\int _{\Sigma }\tau _{xy}dydz&V_{z}=\int
_{\Sigma }\tau _{xz}dydz\\M_{x}=\int _{\Sigma }(-\tau _{xy}z+\tau
_{xz}y)dydz&M_{y}=\int _{\Sigma }z\sigma _{xx}dydz&M_{z}=\int _{\Sigma }-y\sigma
_{xx}dydz\end{matrix}}}

Ejes usuales para una pieza prismática recta, con una sección transversal recta, a
la que se refieren los esfuerzos de sección.
Donde las tensiones que aparecen son las componentes del tensor tensión para una
pieza prismática:

{\displaystyle [T]_{xyz}={\begin{bmatrix}\sigma _{x}&\tau _{xy}&\tau _{xz}\\\tau

_{xy}&0&0\\\tau _{xz}&0&0\end{bmatrix}}}{\displaystyle
[T]_{xyz}={\begin{bmatrix}\sigma _{x}&\tau _{xy}&\tau _{xz}\\\tau _{xy}&0&0\\\tau
_{xz}&0&0\end{bmatrix}}}

Elementos bidimensionales
Para elementos bidimensionales es común tomar un sistema de dos coordenadas
(cartesiano o curvilíneo) coincidentes con la superficie media, estando la tercera
coordenada alineada con el espesor. Para una placa plana de espesor 2t y con un
sistema de coordenadas en el que el plano XY coincide con su plano medio. Los
esfuerzos se componen de 4 esfuerzos de membrana {\displaystyle
n_{uu},n_{uv},n_{vu},n_{vv}\,}{\displaystyle n_{uu},n_{uv},n_{vu},n_{vv}\,} (o
esfuerzos axiles por unidad de área), 4 momentos flectores y 2 esfuerzos cortantes.
Los esfuerzos de membrana usando un conjunto de coordenadas ortogonales
{\displaystyle (u,v)\,}{\displaystyle (u,v)\,} sobre una lámina de Reissner-
Mindlin:

{\displaystyle {\begin{cases}n_{uu}=\int _{-t}^{+t}\left(1-{\frac {z}

{R_{u}}}\right)\sigma _{uu}\ dz&n_{vu}=\int _{-t}^{+t}\left(1-{\frac {z}
{R_{v}}}\right)\sigma _{uv}\ dz\\n_{uv}=\int _{-t}^{+t}\left(1-{\frac {z}
{R_{u}}}\right)\sigma _{uv}\ dz&n_{vv}=\int _{-t}^{+t}\left(1-{\frac {z}
{R_{v}}}\right)\sigma _{vv}\ dz\end{cases}}}{\displaystyle
{\begin{cases}n_{uu}=\int _{-t}^{+t}\left(1-{\frac {z}{R_{u}}}\right)\sigma _{uu}\
dz&n_{vu}=\int _{-t}^{+t}\left(1-{\frac {z}{R_{v}}}\right)\sigma _{uv}\
dz\\n_{uv}=\int _{-t}^{+t}\left(1-{\frac {z}{R_{u}}}\right)\sigma _{uv}\
dz&n_{vv}=\int _{-t}^{+t}\left(1-{\frac {z}{R_{v}}}\right)\sigma _{vv}\
dz\end{cases}}}

Donde {\displaystyle R_{u},R_{v}}{\displaystyle R_{u},R_{v}} son los radios de

curvatura en cada una de las direcciones coordenadas y z es la altura sobre la
superficie media de la lámina. Los esfuerzos cortantes y los momentos flectores por
unidad de área vienen dados por:

{\displaystyle {\begin{cases}v_{u}=\int _{-t}^{+t}\left(1-{\frac {z}

{R_{u}}}\right)\sigma _{uz}\ dz&v_{v}=\int _{-t}^{+t}\left(1-{\frac {z}
{R_{v}}}\right)\sigma _{vz}\ dz\\m_{uu}=\int _{-t}^{+t}\left(1-{\frac {z}
{R_{u}}}\right)\sigma _{uu}z\ dz&m_{vu}=\int _{-t}^{+t}\left(1-{\frac {z}
{R_{v}}}\right)\sigma _{uv}z\ dz\\m_{uv}=\int _{-t}^{+t}\left(1-{\frac {z}
{R_{u}}}\right)\sigma _{uv}z\ dz&m_{vv}=\int _{-t}^{+t}\left(1-{\frac {z}
{R_{v}}}\right)\sigma _{vv}z\ dz\end{cases}}}{\displaystyle
{\begin{cases}v_{u}=\int _{-t}^{+t}\left(1-{\frac {z}{R_{u}}}\right)\sigma _{uz}\
dz&v_{v}=\int _{-t}^{+t}\left(1-{\frac {z}{R_{v}}}\right)\sigma _{vz}\
dz\\m_{uu}=\int _{-t}^{+t}\left(1-{\frac {z}{R_{u}}}\right)\sigma _{uu}z\
dz&m_{vu}=\int _{-t}^{+t}\left(1-{\frac {z}{R_{v}}}\right)\sigma _{uv}z\
dz\\m_{uv}=\int _{-t}^{+t}\left(1-{\frac {z}{R_{u}}}\right)\sigma _{uv}z\
dz&m_{vv}=\int _{-t}^{+t}\left(1-{\frac {z}{R_{v}}}\right)\sigma _{vv}z\
dz\end{cases}}}

El tensor tensión de una lámina general para la que valen las hipótesis de
Reissner-Mindlin es:

{\displaystyle [T]_{uvz}={\begin{bmatrix}\sigma _{uu}&\sigma _{uv}&\sigma

_{uz}\\\sigma _{uv}&\sigma _{vv}&\sigma _{vz}\\\sigma _{uz}&\sigma
_{vz}&0\end{bmatrix}}}{\displaystyle [T]_{uvz}={\begin{bmatrix}\sigma _{uu}&\sigma
_{uv}&\sigma _{uz}\\\sigma _{uv}&\sigma _{vv}&\sigma _{vz}\\\sigma _{uz}&\sigma
_{vz}&0\end{bmatrix}}}

Un caso particular de lo anterior lo constituyen las láminas planas cuya

deformación se ajusta a la hipótesis de Love-Kirchhoff, caracterizada por que el
vector normal a la superficie media deformada coincide con la normal deformada. Esa
hipótesis es una muy buena aproximación cuando los esfuerzos cortantes son
despreciables y en ese caso los momentos flectores por unidad de área en función de
las tensiones vienen dados por:

{\displaystyle {\begin{matrix}m_{x}=\int _{-t}^{t}z\sigma _{xx}dz&m_{y}=\int _{-

t}^{t}z\sigma _{yy}dz&m_{yx}=m_{xy}=\int _{-t}^{t}z\sigma _{xy}dz\end{matrix}}}
{\displaystyle {\begin{matrix}m_{x}=\int _{-t}^{t}z\sigma _{xx}dz&m_{y}=\int _{-
t}^{t}z\sigma _{yy}dz&m_{yx}=m_{xy}=\int _{-t}^{t}z\sigma _{xy}dz\end{matrix}}}

Donde las tensiones que aparecen son las componentes del tensor tensión para una
lámina de Love-Kirchhoff:

{\displaystyle [T]_{xyz}={\begin{bmatrix}\sigma _{xx}&\sigma _{xy}&0\\\sigma

_{xy}&\sigma _{yy}&0\\0&0&0\end{bmatrix}}}{\displaystyle
[T]_{xyz}={\begin{bmatrix}\sigma _{xx}&\sigma _{xy}&0\\\sigma _{xy}&\sigma
_{yy}&0\\0&0&0\end{bmatrix}}}

Ecuaciones de equilibrio
Las ecuaciones de equilibrio de la resistencia de materiales relacionan los
esfuerzos internos con las fuerzas exteriores aplicadas. Las ecuaciones de
equilibrio para elementos lineales y elementos bidimensionales son el resultado de
escribir las ecuaciones de equilibrio elástico en términos de los esfuerzos en
lugar de las tensiones.

Las ecuaciones de equilibrio para el campo de tensiones generales de la teoría de

la elasticidad lineal:

{\displaystyle {\frac {\partial \sigma _{xx}}{\partial x}}+{\frac {\partial \sigma

_{xy}}{\partial y}}+{\frac {\partial \sigma _{xz}}{\partial z}}=b_{x}}
{\displaystyle {\frac {\partial \sigma _{xx}}{\partial x}}+{\frac {\partial \sigma
_{xy}}{\partial y}}+{\frac {\partial \sigma _{xz}}{\partial z}}=b_{x}}
{\displaystyle {\frac {\partial \sigma _{yx}}{\partial x}}+{\frac {\partial \sigma
_{yy}}{\partial y}}+{\frac {\partial \sigma _{yz}}{\partial z}}=b_{y}}
{\displaystyle {\frac {\partial \sigma _{yx}}{\partial x}}+{\frac {\partial \sigma
_{yy}}{\partial y}}+{\frac {\partial \sigma _{yz}}{\partial z}}=b_{y}}
{\displaystyle {\frac {\partial \sigma _{zx}}{\partial x}}+{\frac {\partial \sigma
_{zy}}{\partial y}}+{\frac {\partial \sigma _{zz}}{\partial z}}=b_{z}}
{\displaystyle {\frac {\partial \sigma _{zx}}{\partial x}}+{\frac {\partial \sigma
_{zy}}{\partial y}}+{\frac {\partial \sigma _{zz}}{\partial z}}=b_{z}}

Si en ellas se trata de substituir las tensiones por los esfuerzos internos, se

llega entonces a las ecuaciones de equilibrio de la resistencia de materiales. El
procedimiento, que se detalla a continuación, es ligeramente diferente para
elementos unidimensionales y bidimensionales.
Ecuaciones de equilibrio en elementos lineales rectos
En una viga recta horizontal, alineada con el eje X, y en la que las cargas son
verticales y situadas sobre el plano XY, las ecuaciones de equilibrio relacionan el
momento flector (Mz), el esfuerzo cortante (Vy) con la carga vertical (qy) y tienen
la forma:

{\displaystyle {\frac {dM_{z}}{dx}}=V_{y}\qquad \land \qquad {\frac {dV_{y}}{dx}}=-

q_{y}\qquad \Rightarrow \qquad {\frac {d^{2}M_{z}}{dx^{2}}}=-q_{y}}{\displaystyle
{\frac {dM_{z}}{dx}}=V_{y}\qquad \land \qquad {\frac {dV_{y}}{dx}}=-q_{y}\qquad
\Rightarrow \qquad {\frac {d^{2}M_{z}}{dx^{2}}}=-q_{y}}

Ecuaciones de equilibrio en elementos planos bidimensionales

Las ecuaciones de equilibrio para elementos bidimensionales (placas) en flexión
análogas a las ecuaciones de la sección anterior para elementos lineales (vigas)
relacionan los momentos por unidad de ancho (mx, my, mxy), con los esfuerzos
cortantes por unidad de ancho (vx, my) y la carga superficial vertical (qs):

{\displaystyle {\begin{matrix}{\cfrac {\partial m_{x}}{\partial x}}+{\cfrac

{\partial m_{xy}}{\partial y}}=v_{x}\\{\cfrac {\partial m_{xy}}{\partial x}}+
{\cfrac {\partial m_{y}}{\partial y}}=v_{y}\end{matrix}}\quad \land \quad {\frac
{\partial v_{x}}{\partial x}}+{\frac {\partial v_{y}}{\partial y}}=-q_{s}\qquad
\Rightarrow \qquad {\frac {\partial ^{2}m_{x}}{\partial x^{2}}}+2{\frac {\partial
^{2}m_{xy}}{\partial y\partial x}}+{\frac {\partial ^{2}m_{y}}{\partial y^{2}}}=-
q_{s}}{\displaystyle {\begin{matrix}{\cfrac {\partial m_{x}}{\partial x}}+{\cfrac
{\partial m_{xy}}{\partial y}}=v_{x}\\{\cfrac {\partial m_{xy}}{\partial x}}+
{\cfrac {\partial m_{y}}{\partial y}}=v_{y}\end{matrix}}\quad \land \quad {\frac
{\partial v_{x}}{\partial x}}+{\frac {\partial v_{y}}{\partial y}}=-q_{s}\qquad
\Rightarrow \qquad {\frac {\partial ^{2}m_{x}}{\partial x^{2}}}+2{\frac {\partial
^{2}m_{xy}}{\partial y\partial x}}+{\frac {\partial ^{2}m_{y}}{\partial y^{2}}}=-
q_{s}}