Notas Medios

Notas de Clase: Mecánica de los Medios Continuos
Grupo de Investigación en Mecánica Aplicada
Grupo de Investigación en Mecánica Aplicada

Departamento de Ingenierı́a Civil
Escuela de Ingenierı́as
Universidad EAFIT
Medellı́n, Colombia
2020
Índice general
1. Conceptos Preliminares 5
1.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.2. Concepto de densidad generalizada . . . . . . . . . . . . . . . . . . . . 6
1.3. Conexión entre el continuo fı́sico y el continuo matemático . . . . . . . 7
1.4. Caracterı́sticas generales de un modelo del medio continuo . . . . . . . 9
1.5. Conceptos Preliminares . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.5.1. Hipótesis Básica de Continuidad . . . . . . . . . . . . . . . . . . 18
1.5.2. Partı́cula o Punto Material . . . . . . . . . . . . . . . . . . . . . 18
1.5.3. Punto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
1.5.4. Cuerpo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2. Análisis de Tensiones 21
2.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
2.2. Interacciones mecánicas-Origen del concepto de tensión . . . . . . . . . 22
2.2.1. Fuerzas de cuerpo . . . . . . . . . . . . . . . . . . . . . . . . . . 23
2.2.2. Fuerzas de superficie-El vector de tracciones . . . . . . . . . . . 24
2.2.3. Primer postulado de Cauchy (o Principio de Euler-Cauchy 1822) 28
2.2.4. Fuerzas por unidad de superficie sobre un plano con dirección
normal n̂ . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
2.2.5. Tensor de tensiones . . . . . . . . . . . . . . . . . . . . . . . . . 33
2.2.6. Transformación por rotación. Caso general . . . . . . . . . . . . 37
2.2.7. Transformación por rotación 2D. Cı́rculo de Mohr . . . . . . . . 40
2.3. Valores y direcciones principales para un estado de tensiones en 3D . . 55
2.3.1. Circulo de Mohr en 3D . . . . . . . . . . . . . . . . . . . . . . . 61
2.4. Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
3. Equilibrio 83
3.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83
3.2. Ecuaciones de Equilibrio diferencial . . . . . . . . . . . . . . . . . . . . 84
3.2.1. Coordenadas cilı́ndricas . . . . . . . . . . . . . . . . . . . . . . 90
i
ÍNDICE GENERAL ii
3.3. Soluciones clásicas para sólidos . . . . . . . . . . . . . . . . . . . . . . 98

3.3.1. Barra prismática soportada en la superficie superior y sometida
a la acción de su propio peso . . . . . . . . . . . . . . . . . . . . 98
3.3.2. Viga en voladizo . . . . . . . . . . . . . . . . . . . . . . . . . . 99
3.3.3. Cilindro sometido a presión interna y externa . . . . . . . . . . 100
3.3.4. Flamant . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101
3.4. Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
4. Análisis de Deformaciones 118

4.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 118
4.1.1. Concepto intuitivo de deformación . . . . . . . . . . . . . . . . 120
4.2. Transformaciones lineales . . . . . . . . . . . . . . . . . . . . . . . . . . 123
4.2.1. Componente rotacional . . . . . . . . . . . . . . . . . . . . . . . 125
4.2.2. Componente simétrica . . . . . . . . . . . . . . . . . . . . . . . 127
4.2.3. Componentes escalar y distorsional . . . . . . . . . . . . . . . . 129
4.3. Tensor de deformaciones unitarias . . . . . . . . . . . . . . . . . . . . . 135
4.3.1. Rotación de cuerpo rı́gido . . . . . . . . . . . . . . . . . . . . . 138
4.3.2. Componente simétrica ε . . . . . . . . . . . . . . . . . . . . . . 139
4.3.3. Componente escalar, volumétrica o isotrópica . . . . . . . . . . 140
4.3.4. Componente distorsional . . . . . . . . . . . . . . . . . . . . . . 141
4.3.5. Ejemplo: Barra sometida a su peso propio . . . . . . . . . . . . 146
4.4. Ejercicios propuestos . . . . . . . . . . . . . . . . . . . . . . . . . . . . 150
5. Elasticidad 164
5.1. Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 164
5.2. Hipótesis básicas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 165
5.3. Relación tensión-deformación . . . . . . . . . . . . . . . . . . . . . . . 166
5.4. Idealizaciones en 2D . . . . . . . . . . . . . . . . . . . . . . . . . . . . 167
5.4.1. Tensión plana . . . . . . . . . . . . . . . . . . . . . . . . . . . . 168
5.4.2. Deformación plana . . . . . . . . . . . . . . . . . . . . . . . . . 169
5.5. Ejemplos resueltos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 170
5.6. Ejercicios . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 180
Referencias 186
ÍNDICE GENERAL 1
Presentación
Gran parte de los problemas de Ingenierı́a se encuentran enmarcados en el correcto

cálculo, determinación y análisis de la distribución de desplazamientos, tensiones y
cambios de configuración para medios deformables sometidos a diferentes tipos de exci-
taciones. La mecánica del medio continuo es el modelo matemático desarrollado a partir
de principios fı́sicos introducidos en la mecánica Newtoniana que permite resolver dicho
problema. En el mejor escenario posible y dependiendo del grado de complejidad de las
geometrı́as de los medios a tratar, de las condiciones de frontera o del comportamiento
mismo de los materiales constitutivos, se pueden lograr soluciones analı́ticas a dichos
problemas. Sin embargo, en las situaciones más complejas es necesario abordarlo me-
diante el uso de herramientas computacionales o experimentales. En cualquier caso, el
Ingeniero debe tener el suficiente conocimiento y dominio de las hipótesis básicas de los
diferentes modelos matemáticos, ası́ como de las ecuaciones que gobiernan el problema
para poder ası́ resolver, plantear, interpretar y analizar la amplia gama de posibles
situaciones en el campo de la mecánica de medios continuos.
Estas notas de clase, desarrolladas por diversos docentes adscritos al Grupo de Inves-
tigación en Mecánica Aplicada, están orientadas a servir como herramienta de apoyo a
los estudiantes del curso Mecánica de los Medios Continuos del pregrado en Ingenierı́a
Civil de la Universidad EAFIT. Los contenidos de las mismas no son otra cosa que
una recopilación más o menos concienzuda de material disponible en diferentes textos
estructurado de manera que el material sea de fácil acceso a estudiantes con una forma-
ción fı́sica y matemática correspondiente al cuarto semestre de pregrado. Es importante
aclarar que estas notas no tienen el rigor en cuanto a formato y redacción de un texto
y que no buscan reemplazar la gran cantidad de excelentes contribuciones en Mecánica
de los Medios Continuos disponibles en la literatura (ver sección de Bibliografı́a).
Las notas están organizadas de manera consistente con el curso cuya estructura y con-
tenido se describe a continuación.
Objetivos generales del curso
Permitir al estudiante entender el modelo del medio continuo como una estrategia
ingenieril de solución al problema central de la mecánica en un sistema de “infinitas”
partı́culas de manera que se genere un soporte fı́sico-matemático lo suficientemente
robusto para enfrentar problemas de la mecánica aplicada.
ÍNDICE GENERAL 2
Objetivos especı́ficos
Desarrollar el concepto matemático y fı́sico de un tensor de orden 2.
Identificar las hipótesis fundamentales y que de manera conjunta con el concepto

del continuo matemático reducen los conceptos de fuerza y desplazamientos rela-
tivos entre partı́culas a los conceptos de tensión y deformación y al mismo tiempo
re-visitar los conceptos de tensión y deformación y conocer algunas de las leyes
constitutivas que gobiernan la relación entre estos.
Reducir el problema central de la mecánica para un sistema de infinitas partı́culas,

al de determinar los campos de desplazamientos, tensiones y deformaciones.
Identificar matemáticamente el problema de la determinación de desplazamientos,

tensiones y deformaciones en un medio continuo como un Problema de Valores en
la Frontera (PVF).
Metodologı́a
Exposiciones magistrales, Elaboración de Ejemplos, Lecturas asignadas. El curso se di-

vide en 2 partes. En la parte I se revisan las hipótesis básicas que sustentan el modelo y
se hace un breve repaso de álgebra vectorial. En la parte II se presentan los contenidos
teóricos del modelo. En esta se parte de una revisión de los conceptos de tensión y
deformación y obteniendo como resultado final las ecuaciones gobernantes o ecuaciones
de campo. Posteriormente se ligan las tensiones y deformaciones a través de un modelo
constitutivo o ley de Hooke e introduce la hipótesis simplificadora de suponer el equili-
brio en la configuración no-deformada del medio para llegar a la teorı́a linealizada de la
elasticidad. A la luz de esta teorı́a idealizada se estudian finalmente algunas soluciones
analı́ticas o cerradas, discutiendo brevemente los métodos de solución y concentrándose
más en el análisis y entendimiento de la solución misma.
ÍNDICE GENERAL 3
Contenido de la materia
I. Introducción
1. Presentación del problema (1 semana).
1.1 Introducción y motivación del modelo del continuo.
II. FUNDAMENTOS TEÓRICOS
2.a Análisis de tensiones. Parte 1 (4 semanas).
2.1 Concepto de tensión (Definición, Primer postulado de Cauchy, Definición

de tensor).
2.2 Ecuaciones de transformación y esfuerzos extremos en 2D y 3D. Cı́rculo

de Mohr
2.3 Solución de problemas
◦ Parcial 01.
2.b Análisis de tensiones. Parte 2 (4 semanas).
2.5 Ecuaciones de equilibrio en un medio continuo.
2.6 Estudio e interpretación de soluciones de tensiones.
◦ Parcial 02.
3. Análisis de Deformaciones (5 semanas).
3.1 Transformaciones lineales.
3.2 Definición de concepto de deformación.
3.3 Tensor de deformaciones.

ÍNDICE GENERAL 4
3.4 Teorı́a de la Elasticidad. Ley de Hooke.
3.5 Problemas.
◦ Parcial 03.
4. Proyecto docente (2 semanas).
4.1 Estudio de solución numérica.
◦ Examen Final.
Capı́tulo 1
Conceptos Preliminares
1.1. Introducción
En esta sección se presentan las hipótesis básicas que permiten construir el modelo
del medio continuo. Dicho modelo debe entenderse como una propuesta para abordar
problemas de mecánica de sistemas de infinitas partı́culas usando como herramienta
matemática las ideas de continuidad propuestas en el cálculo diferencial.
Al final de esta sección el estudiante debe estar en la capacidad de:
• Entender el paso de variables discretas a variables continuas en términos de den-

sidades.
• Identificar la diferencia entre el concepto de partı́cula en el modelo de la mecánica

Newtoniana y el concepto de punto material en el modelo del continuo.
• Identificar las funciones fuente o causa y efecto o respuesta en el modelo del

continuo.
• Identificar los lı́mites de aplicabilidad del modelo del continuo.
5
CAPÍTULO 1. CONCEPTOS PRELIMINARES 6
1.2. Concepto de densidad generalizada
A una escala de tamaño lo suficientemente pequeña la materia goza de una composición

molecular/atómica/microestructural o, en general, de carácter discreto. Lo anterior
implica que un volumen finito de material no es otra cosa que una gran colección
de partı́culas, sean estas: moléculas, átomos, partı́culas subatómicas o cualquier otro
ente discreto.
Si para esta colección se desea calcular la distribución en el tiempo de alguna función

respuesta o efecto β resultante de la consideración de alguna función causa α, será
entonces necesario resolver un sistema acoplado de ecuaciones diferenciales ordinarias
tan grande (en número de ecuaciones) como partı́culas posea la colección. Este sistema
de ecuaciones es obtenido a partir de algún principio fı́sico fundamental, conectando
las funciones α y β, y postulado a partir de observación experimental. En términos
generales es posible escribir este principio como en la ecuación (1.1),

dβ
f = α, (1.1)
dt
en la cual β es la cantidad fı́sica que experimenta cambios en el tiempo y que denomi-
naremos en lo que sigue como función flujo, mientras que α es la cantidad fı́sica que
genera los cambios y la denominaremos la función interacción.
Por ejemplo, en el problema mecánico la ecuación (1.1) corresponderı́a a la segunda ley

de Newton en cuyo caso identificamos α ≡ F~ , siendo F~ la fuerza actuante sobre una
partı́cula de masa m y β~ ≡ mV~ la cantidad de movimiento asociada con la velocidad
V~ de esta partı́cula. En este caso f estarı́a definida como f (x) = x para un argumento
x.
Por tanto, si se dispone de un postulado fı́sico fundamental como el dado en la ecua-

ción (1.1), la determinación de las funciones flujo β para cada una de las partı́culas es
simplemente una tarea operativa. Sin embargo, la gran cantidad de partı́culas de las
cuales puede estar compuesto el sistema hacen que este tratamiento “exacto” a nivel
discreto sea imposible de realizar en términos prácticos. Es decir, la posibilidad de una
solución exacta al problema se encuentra limitada por razones operacionales o compu-
tacionales. En términos prácticos es válido considerar el volumen finito de material
como un sistema de infinitas partı́culas para cuya solución se requiere tratar con infi-
nitas ecuaciones diferenciales ordinarias. Este planteamiento puede entenderse como la
pregunta o problema general que se busca resolver mediante un modelo de los medios
continuos. Este busca una solución práctica al problema de calcular una distribución
de funciones flujo, debidas a unas interacciones en un sistema de “infinitas” partı́culas.
Un tratamiento aproximado se hace posible gracias a la Hipótesis Básica de Continuidad

Matemática. Esta encuentra su motivación en el carácter continuo que soporta el espacio
matemático; entre dos puntos, sin importar que tan cercanos se encuentren, siempre
es posible encontrar otro punto. Se tiene, entonces, la definición clásica de continuidad
matemática o concepto de lı́mite y en la cual el punto podrı́a definirse a partir de una
celda volumétrica ∆V cuyo tamaño se reduce por un proceso de paso al lı́mite
lı́m ∆V = Punto . (1.2)

∆V →0
La generalización de este concepto al sistema de infinitas partı́culas se reduce a despre-

ciar el carácter discreto de la materia y asumir que estas se encuentran distribuidas en el
espacio matemático de manera continua, es decir, esta generalización equivale a afirmar
que las infinitas partı́culas ocupan un volumen finito sin dejar vacı́os ni discontinuidades
entre las mismas. En lo que sigue denominaremos esta distribución fı́sicamente ficticia
El Continuo Fı́sico. En otras palabras, la celda de volumen ∆V independientemente
de su tamaño siempre estará compuesta del mismo material ya que se ha despreciado
completamente todo ente micro-estructural.
Podemos en este punto precisar el objetivo fundamental de un modelo de los medios

continuos:
Un modelo de los medios continuos pretende estudiar la distribución espacio-

temporal de las interacciones y flujos resultantes en un sistema discreto de
infinitas partı́culas que es aproximado mediante un proceso de promediado
basado en la hipótesis básica de continuidad y por ende, convirtiendo las
variables involucradas a nivel discreto en campos o distribuciones.
1.3. Conexión entre el continuo fı́sico y el continuo

matemático
Retomando el tratamiento del problema de determinación de una función respuesta

β para un sistema de muchas partı́culas sometidas a algún tipo de excitación externa
o interacción α, es claro que, desde el punto de vista operativo, es imposible resolver
las ecuaciones resultantes (discretas) para un sistema de infinitas partı́culas, por lo que
es necesario recurrir a la hipótesis de continuidad matemática con dos consecuencias
inmediatas:
Se renuncia al concepto de partı́cula.
Se hace disponible la potente herramienta del continuo matemático o cálculo

diferencial.
Por ejemplo, en el caso de la mecánica de los medios continuos no tendrá sentido hablar
términos como “la fuerza aplicada a una partı́cula” en el sentido discreto. Más aún, en
las ecuaciones gobernantes a nivel del continuo este concepto mismo de fuerza discreta
desaparece matemáticamente y es necesario utilizar representaciones equivalentes que
tengan significado dentro del contexto de los medios continuos, es decir funciones, como
en este caso la función tensorial de tensiones.
La hipótesis de continuidad genera entonces un modelo que producirá resultados co-

rrespondientes a una aproximación con respecto a una solución discreta1 . El orden del
error en la aproximación, será siempre relativo y estará fuertemente ligado a la relación
entre el tamaño caracterı́stico del espécimen en cuestión y el tamaño caracterı́stico de
la micro-estructura que se ha despreciado. Por ejemplo, considere el caso de una viga
de concreto. En ésta, la menor dimensión caracterı́stica del espécimen corresponderı́a al
espesor h mientras que la mayor dimensión caracterı́stica de la micro-estructura seria el
diámetro del agregado grueso φ. La aproximación introducida al asumir el medio como
un continuo fı́sico dará errores despreciables siempre que h φ. En la medida en que
h → φ la contribución del agregado a la interacción mecánica es relevante y el modelo
del continuo se vuelve impreciso.
El hecho de renunciar al concepto de partı́cula y de apelar al continuo matemático

implica que toda variable de tipo discreto o asociada al concepto de partı́cula es tratada
ahora como una distribución o función. Replanteando en estos términos y por segunda
vez el objetivo fundamental de un modelo de los Medios Continuos este se deberı́a
reescribir como sigue:
Un modelo de los medios continuos pretende estudiar la variación espacio-

temporal de las distribuciones o funciones de respuesta (o efecto) debidas
a una distribución espacio-temporal de funciones fuente (o causa).
En el caso particular del problema mecánico el interés es describir la respuesta del

sistema en términos de cantidad de movimiento. En este escenario las fuentes co-
rresponden a fuerzas o mejor aún, a distribuciones de fuerzas. Más concretamente, una
vez introducida la hipótesis de continuidad, el objetivo será el de determinar funcio-
1
Si se asume que el resultado que se obtendrı́a de un modelo discreto es exacta y representa el
sistema de interés.
nes continuas representativas de la cantidad de movimiento en respuesta a funciones

representativas de las fuentes o fuerzas.
Se puede afirmar entonces que el carácter infinito del sistema de partı́culas genera
la necesidad de un modelo matemático para describir ahora la interacción2 entre 2
sistemas. Uno de los sistemas es el medio continuo en cuestión (sistema de muchas
partı́culas) que interactúa o en el caso mecánico intercambia cantidad de movimiento
con otro sistema externo. Denominaremos este (otro medio continuo) como el medio
ambiente [1].
1.4. Caracterı́sticas generales de un modelo del me-

dio continuo
Para apoyar el análisis considérese un medio infinito ocupando un volumen B(t) y

que representa en este caso el medio ambiente, figura 1.1. Este medio se encuentra
acotado por una superficie S∞ localizada a una distancia infinita de cualquier punto del
medio. Consideramos además un medio continuo (el sistema) ocupando un volumen
V (t) y localizado al interior de este espacio infinito. La interfase entre el sistema y el
medio ambiente corresponde a una superficie S(t)3 . Las preguntas que nos proponemos
responder son:
¿Cómo describir la interacción o intercambio de cualquier variable fı́sica del tipo

respuesta entre el medio (ambiente) infinito B(t) y el medio continuo V (t) debido
a una distribución de variables de tipo fuente o causa?
¿Cómo expresar esta interacción como una función continua, es decir llevada a
nivel del punto matemático como se hace en el cálculo diferencial?
Con referencia a la figura 1.1 utilizaremos la siguiente descripción de las regiones in-
volucradas. La región localizada del lado positivo del vector normal n̂(+) a S(t) se
denominará en lo que sigue la región externa. De manera consistente, la región locali-
zada del lado negativo del vector normal n̂(−) a S(t) se denominará la región interna.
Deseamos entonces caracterizar la interacción entre las regiones interna y externa de-
2
En el caso mecánico una fuerza no es otra cosa que la representación equivalente de una interacción
entre dos sistemas.
3
En la figura 1.1 la frontera del medio continuo representada por S podrı́a ser por ejemplo la de
una viga, una presa o si se quiere tan solo una frontera matemática.
bida a la existencia de una fuente localizada ya sea sobre B(t) o sobre V (t). Además
queremos describir esta interacción como una función continua.
n̂( +)
Región
n̂( −)
externa
S(t)
S∞ Región
B(t ) V (t ) interna
Figura 1.1. Medio continuo interactuando con un espacio infinito.
En el problema planteado, tanto las variables fuente como las variables causa se en-
cuentran conectadas en términos discretos a través de una ley como la indicada en
ecuación (1.1). Ya señalamos que para tratar estas variables de carácter discreto co-
mo una función, se procede partiendo del concepto de lı́mite. Para esto supondremos
que realizaremos particiones del medio continuo en celdas de volumen (cada una de las
cuales será también un sistema de infinitas partı́culas dado que se ha despreciado la
microestructura a todas las escalas) y que para eliminar la dependencia del tamaño de
la partición daremos a estas celdas un tamaño cada vez más pequeño hasta alcanzar un
estado lı́mite en el cual supondremos que cada celda se convierte en lo que denomina-
remos por ahora un punto material para diferenciarlo del punto matemático. Es decir
a la celda de volumen de la ecuación (1.3) le estamos asignando material.
lı́m ∆V = P unto material (1.3)

∆V →0
A la luz de la continuidad matemática el punto material definido en la ecua-

ción (1.3) no posee tamaño. Este sin embargo encierra o esta conformado por
infinitas partı́culas; es decir, el punto material es a la vez un medio continuo.
Matemáticamente dicho punto siempre existe, esta es la hipótesis básica de conti-

nuidad. Fı́sicamente en el proceso de cálculo del lı́mite dado en la ecuación (1.3)
eventualmente se alcanzará un tamaño finito donde la hipótesis de continuidad
comienza a colapsar y el lı́mite deja de ser válido; a partir de este tamaño el efec-
to de la micro-estructura despreciada desde el inicio comienza a ser importante
y en el lı́mite eventualmente la celda se convierte en una de las partı́culas de la
micro-estructura. De hecho, podrı́a convertirse en un vacı́o.
Si las dimensiones caracterı́sticas del espécimen son varios órdenes de magnitud

mayores que este tamaño lı́mite finito, siempre será valido trabajar con la defini-
ción matemática de punto, ecuación (1.3). Es decir si el ∆V fı́sico es menor en
varios órdenes de magnitud que todas las dimensiones importantes del sistema
de infinitas partı́culas será válida la hipótesis de continuidad fı́sica dada en la
ecuación (1.3).
Regresemos nuevamente al nivel de la partı́cula y mantengamos la designación de α pa-

ra la variable fuente o causa que aparece de forma explı́cita en asocio con una partı́cula
y la cual produce cambios de una variable tipo efecto o respuesta β. Deseamos consi-
derar esta variable como una distribución o función en el medio continuo. Denotaremos
esta función resultante como Φ. Esta puede ser de carácter escalar, vectorial o en gene-
ral tensorial. Recordemos además que las variables α y β están relacionadas mediante
una ley general o principio fundamental

dβ
f =α (1.4)
dt
El tratamiento propuesto en la ecuación (1.4) necesariamente nos lleva a tener que

promediar la variable α sobre las celdas de volumen, es decir a asociarlas al punto, vı́a
el concepto de lı́mite implicado en el cálculo diferencial. Por ejemplo, supondremos que
a una celda volumétrica finita ∆V1 le corresponde un ∆α1 , a una celda inferior ∆V2 un
∆α2 y ası́ sucesivamente de manera tal que en el estado lı́mite del punto en el sentido
de ecuación (1.3) la relación ∆α/∆V tenderá a una función continua Φ y que en el
punto material asume por lo tanto un valor constante.
∆α
lı́m =Φ (1.5)
∆V →0 ∆V
En ecuación (1.5) ∆V tiende a cero mientras que la relación ∆α/∆V se asume

aproximándose al valor dado Φ.
El lı́mite de la ecuación (1.5) no ha sido calculado a través de ningún procedi-

miento matemático explı́cito sino que su existencia se impone en el modelo como
hipótesis de partida. Es decir se asume la existencia de dicho lı́mite.
El proponer la existencia del lı́mite dado por la ecuación (1.5) ha sido motivado
por la hipótesis de continuidad fı́sica y su relación directa con el continuo ma-
temático. En estos términos la existencia de dicho lı́mite podrı́a entenderse como
un postulado fundamental sobre el cual se está construyendo el modelo del medio
continuo.
Si los resultados del modelo construido a partir del lı́mite dado en ecuación (1.5)
logran predecir las observaciones experimentales entonces el modelo resultante
habrá sido exitoso, en caso contario este debe ser enriquecido, mejorado o en el
caso extremo abandonado. Por ejemplo si Φ en ecuación (1.5) no asume un valor
constante sino una variación (por ejemplo lineal) entonces será posible capturar
fluctuaciones lineales (o localizaciones) de Φ en el punto material. Este principio
podrı́a seguirse extendiendo para considerar fluctuaciones más fuertes como por
ejemplo cuadráticas y de orden superior, hasta que eventualmente se alcanzarı́a
un nivel de refinamiento comparable en términos de precisión con un resultado
discreto. Lógicamente, en la medida en que se aumenta el orden de la fluctuación
asumida sobre el punto material aparecen cada vez más ecuaciones.
Finalmente, la nueva función Φ estará asociada con un punto material en el sentido

de la ecuación (1.3) y se tratará ahora como una cantidad de α por unidad de
volumen que varia de manera continua a lo largo de V (t), es decir como un
campo. Se dice entonces que la función Φ representa una densidad volumétrica de
α . Cuando son vistas en asocio con una superficie estas densidades se convierten
en densidades superficiales; cantidad de α por unidad de superficie.
En la figura 1.2 se esquematiza un problema para un espécimen uni-dimensional de

tamaño L. La escala de tamaño del espécimen se denominará en lo que sigue la macro-
escala. En ésta, una función del tipo Φ varı́a a lo largo del dominio matemático ocupado
por el sistema de partı́culas o puntos materiales. En un punto material cualquiera
del dominio L, la variable Φ asume un valor constante de manera consistente con el
lı́mite dado en la ecuación (1.3). Sin embargo ya es claro que fı́sicamente sobre el
punto material se presentan o existen entes de carácter discreto y sobre los cuales las
variables presentan fuertes fluctuaciones. Por ejemplo en el inserto de la figura 1.2 se
denota por ` la escala de tamaño sobre la cual se presentan las fluctuaciones. Este
dominio se denominará en lo que sigue la micro-estructura. En esta misma figura se
esquematiza también mediante la lı́nea punteada el valor promedio de las fluctuaciones
y correspondiente al valor constante asignado en el punto material de la parte superior

de la Figura. Como conclusión las fluctuaciones asociadas con la micro-estructura o
carácter discreto de la materia no serán capturadas por el modelo del continuo válido
sobre el dominio o a la escala de tamaño L.
x
L
x
l

Figura 1.2. Dominio del Medio Continuo representando la Macro-escala y la micro-
escala.
A través del concepto del lı́mite establecido en ecuación (1.3) lo que se hace entonces
es asignar al volumen o punto material matemático el valor promedio acomodado sobre
el punto o volumen material fı́sico. En estos términos es claro que la distribución de
la función Φ a nivel macro representa una aproximación con respecto a la distribución
a nivel micro. Las predicciones a nivel del continuo basado en la ecuación (1.3) serán
válidas mientras se de la condición L `.
Más aún, nos referiremos al volumen o punto material fı́sico (como el de la parte supe-
rior de la figura 1.2) como Volumen Representativo del Material (VRM). De acuerdo
con esto se tiene entonces que un continuo clásico puede entenderse como una varie-
dad diferenciable4 de puntos materiales abstractos o Volúmenes Representativos del
Material.
De acuerdo con la hipótesis de continuidad matemática estos puntos materiales a la

vez están conformados por un número infinito de puntos materiales de manera que
matemáticamente o por abstracción siempre será posible seguir descendiendo en tamaño
sin nunca encontrar los entes micro-estructurales o de tipo discreto. Sin embargo en
otros modelos del continuo (denominados como no-clásicos) el VRM puede ser dotado
de alguna capacidad de describir (aunque aún en términos del continuo) cierto tipo
de fluctuaciones micro-estructurales. Más adelante se hará referencia a este tipo de
propuestas.
Volumen Representativo del Material5
Habiendo identificado el VRM como el bloque constitutivo fundamental de un modelo

del medio continuo es útil en esta sección dar un poco más de precisión a esta definición.
Para esto se presentan algunas micro-estructuras relativas y se reitera sobre la mane-
ra de producir un continuo clásico en términos de asignaciones de variables a puntos
materiales abstractos e ideales correspondientes a valores promedios calculados sobre
puntos materiales fı́sicos e identificables en el laboratorio.
De la discusión precedente es claro que implı́cita en la estrategia utilizada en la hipóte-

sis básica de continuidad se encuentra un proceso de homogenización o transición de
menor a mayor escala de tamaño. Por homogenización nos referimos a la asignación
de propiedades fı́sicas a los puntos materiales correspondientes a un promedio de los
valores asumidos por varias partı́culas sobre un volumen finito de material.
Como ya se mencionó y de acuerdo con postulados fundamentales de la mecánica

cuántica los sistemas formados por átomos, electrones, protones, neutrones y partı́culas
subatómicas carecen de continuidad. Es decir, a este nivel no se logra una precisión
aceptable si se opera bajo la hipótesis de continuidad.
En la dirección ascendente en cuento a tamaño el siguiente nivel de fragmentación se

4
Variedad diferenciable o, en inglés, manifold. Correspondiente a una estructura topológica
5
Algunos casos de determinación del VRM se presentan en las referencias [2, 3].
identifica en la molécula. La mecánica clásica de los medios continuos ha sido construida

sobre la hipótesis de que a este nivel efectivamente se logran resultados con la precisión
adecuada y de que por encima de este nivel no existen otros entes micro-estructurales. Es
decir, para volúmenes de tamaños iguales y superiores a la molécula las aproximaciones
basadas en la hipótesis básica de continuidad usadas en la mecánica clásica permitirı́an
llegar a resultados aceptables.
Sin embargo desarrollos relativamente recientes, motivados por las nuevas tendencias
hacia la miniaturización, han revelado escalas de tamaño varios órdenes de magnitud
superior a la escala molecular donde la mecánica clásica también comienza a agotar su
nivel de precisión. Este es el caso particular de los denominados materiales con micro-
estructura. De lo anterior se desprende el hecho de que existe un volumen mı́nimo de
material por debajo del cual es impreciso abordar el problema de las infinitas partı́culas
a partir de la hipótesis básica de continuidad. Desde el punto de vista matemático dicho
volumen será siempre considerado infinitamente pequeño.
Pudiéramos pensar en esta escala como la correspondiente a la microestructura del ma-

terial a través de la cual seria posible identificar de manera discreta vacı́os, impurezas,
defectos o en general heterogeneidades. Supongamos que en esta escala las heterogenei-
dades se encuentran separadas a unas distancias promedio ` (la micro-escala definida
en la sección anterior). Por ejemplo, consideremos un volumen de material de tamaño
caracterı́stico comparable a esta distancia `. Es claro que si se estuviera analizando un
espécimen con dimensiones caracterı́sticas L (la macro-escala) cercanas o comparables
con ` la validez de un tratamiento con base en la hipótesis de continuidad serı́a fácilmen-
te cuestionable pues se generarı́an fluctuaciones como las de la figura 1.2. La figura 1.3
ilustra una de estas microestructuras correspondiente a una aleación plomo/estaño en
la cual la separación de los defectos es del orden de 5.0µm. Si se estuviera analizando
vı́a la mecánica del medio continuo un espécimen con un tamaño cercano al de estas
micro-estructuras, los resultados no capturarı́an concentraciones de las variables del
tipo Φ sobre las diferentes fases, por ejemplo sobre la fase de plomo o de estaño. La
pregunta es entonces: ¿a partir de que tamaño L del espécimen es válido comenzar a
suponer que efectivamente se tiene un comportamiento representable como un medio
continuo?. En otras palabras, a partir de que tamaño L es válido afirmar que si se divide
el medio en celdas infinitesimales a las cuales se les confieren propiedades mecánicas
promedio, las predicciones del modelo estarán dentro de rangos de precisión aceptables
desde el punto de vista de la ingenierı́a?

Figura 1.3. Fotografı́a aumentada de una aleación de plomo/estaño utilizada en
aplicaciones de microelectrónica.
La respuesta depende claramente del tipo de micro-estructura inherente en cada mate-

rial. En los desarrollos tempranos de la mecánica del medio continuo el problema fue
evadido asumiendo que se estaba a escalas varias veces por encima de las distancias
inter-moleculares despreciando la existencia de entes fundamentales de tamaño supe-
rior. Aunque para la ingenierı́a clásica este modelo del medio continuo ha resultado
ser lo suficientemente preciso, el desarrollo de tecnologı́as e industrias emergentes6 ha
empezado a llevar al modelo al lı́mite de su precisión. En esta presentación se llama
la atención sobre la existencia de esta escala inferior por debajo de la cual no serı́a
válido proceder vı́a el continuo y se introduce el concepto de volumen representativo
del material como una propiedad mecánica adicional cuyo valor depende de cada ma-
terial. Claramente esta propiedad ` tiene dimensiones de longitud y varı́a de material
a material y de microestructura a microestructura.
En lo que sigue asumiremos que el sistema de infinitas partı́culas (no necesariamente

todas iguales) estará ocupando un espacio fı́sico de manera continua (en el sentido fı́sico
y matemático), o equivalentemente se asumirá continuidad de puntos y de partı́culas;
entre dos partı́culas del sistema, sin importar que tan cercanas se encuentren siempre
será posible asignar otra partı́cula. En otras palabras se le está confiriendo acá a la
partı́cula el carácter de punto material. Este espacio se subdividirá en volúmenes in-
finitesimales ∆V a cada uno de los cuales se les conferirán propiedades (mecánicas)
constantes. Nos referiremos a este volumen como Volumen Representativo del Ma-
terial (VRM)7 . Asumiremos además que siempre estaremos operando en especı́menes
con tamaños superiores a este VRM y que estaremos haciendo mediciones a escalas
igualmente superiores. La presencia de la microestructura y de un punto material fı́sico
6
Como la micro-electrónica y la nano-tecnologı́a
7
VRM: Mı́nimo volumen sobre el cual el material es estadı́sticamente homogéneo, F (x) = x
o VRM se esquematiza en la figura 1.4. Cabe anotar que la ingenierı́a moderna ha

reconocido la potencia práctica de la hipótesis del medio continuo como herramienta
de solución de problemas y ha planteado modelos que intentan ampliar las capacidades
del mismo a escalas bastante cercanas a las del VRM. Estas se han denominado como
teorı́as no-locales, de cinemática enriquecida o simplemente teorı́as no-clásicas de la
mecánica de los medios continuos [4–9].

Figura 1.4. Esquematización del Volumen Representativo del Material y de las ma-
cro y meso-escalas.
Con base en las discusiones anteriores, estos modelos no-clásicos lo que hacen es dotar
al VRM de información referente a la micro-estructura en términos de variables que
varı́an de manera continua, como por ejemplo parámetros relativos a su tamaño, que
permiten describir distribuciones no-constantes localizadas al interior del punto material
matemático lo que permite dar cuenta en términos del continuo de fluctuaciones no
capturables por una propuesta clásica.
1.5. Conceptos Preliminares
Antes de presentar el método general que se utilizará para reducir las variables dis-
cretas a distribuciones, densidades o campos es necesario introducir algunos conceptos
preliminares.
1.5.1. Hipótesis Básica de Continuidad
Aunque la descripción de los cambios de configuración a los cuales se puede someter

un medio continuo se estudiará de manera formal y detallada en secciones posteriores,
será útil hacer acá una referencia inicial a este respecto, en particular en lo concerniente
con la descripción del movimiento en el caso de un medio deformable. En este sentido
resultará útil y común comparar las configuraciones, localizaciones o lugar geométrico
de las infinitas partı́culas del sistema en dos instantes.
Un instante o configuración inicial que denotamos como correspondiente a t = 0

, donde t representa un marco de referencia temporal. Nos referiremos a esta
configuración indistintamente como el medio continuo en t = 0, la configuración de
referencia o la configuración no-deformada. Se asumirá que en esta configuración
el medio continuo se encuentra libre de interacciones con el medio ambiente.
Un instante o configuración posterior t = t. Nos referiremos a esta configuración

indistintamente como el medio continuo en t = t, la configuración actual o
la configuración deformada. Se asume que el medio continuo es llevado a esta
configuración precisamente por la acción de interacciones externas.
Al comparar estas dos configuraciones se asumirá que en ambos estados la cantidad

de materia es la misma. Más aún, siguiendo un enfoque causa-efecto supondremos que
es la acción de las interacciones externas la que produce los cambios de configuración
entre t = 0 y t = t.
La hipótesis básica de continuidad desde el punto de vista fı́sico consiste en suponer que
las infinitas partı́culas que conforman el sistema ocupan o llenan de manera continua
un espacio fı́sico determinado de volumen V y acotado por una superficie o frontera S.
Esta hipótesis fı́sica permite ahora representar todas las propiedades como funciones
continuas y posibilita el uso del cálculo como herramienta matemática. Más adelante
se identificarán algunas consecuencias de esta hipótesis básica.
1.5.2. Partı́cula o Punto Material
En el contexto de la mecánica de los medios continuos por partı́cula o punto material se

estará haciendo referencia a cantidad de materia o porción infinitesimal de un cuerpo
pero sin conceder a esta geometrı́a alguna tal y como se estipula en el lı́mite dado en
la ecuación (1.3). Tampoco se detalla ésta es un átomo, molécula o que tipo de ente
fundamental. A pesar de la carencia de geometrı́a la partı́cula tendrá siempre asociado

un volumen y una densidad de masa. Esta contradicción aparente es una consecuen-
cia del uso del concepto matemático de lı́mite. El punto material podrá experimentar
desplazamiento, velocidad, aceleración, cambios de moméntum o en general flujos e in-
teracciones. El punto material se considerará como el elemento constitutivo primario o
elemental de un modelo de los medios continuos.
1.5.3. Punto
Por punto se hará referencia a una localización geométrica fija con respecto a un sistema
de referencia. El punto como tal tampoco posee geometrı́a alguna y mucho menos tendrá
asociadas cantidades fı́sicas. Sin embargo, en un instante dado, un punto puede estar
siendo ocupado por un punto material con propiedades fı́sicas asociadas. En este caso
hablaremos de la propiedad en el punto pero teniendo en cuenta que en realidad esta
será debida a la contribución del punto material. Más aún, pueden existir interacciones
del punto con el medio externo de manera tal que cuando una partı́cula ocupa este
punto, está siente una contribución adicional o agente modificador de sus propiedades.
1.5.4. Cuerpo
Nos referimos a un cuerpo como un sistema de infinitas partı́culas (ahora puntos ma-
teriales) lo cual no constituye ninguna novedad ya que en cursos básicos de fı́sica se
habı́a trabajado con el caso de cuerpos rı́gidos. Se trataba precisamente de un sistema
de infinitas partı́culas pero en este caso el apelativo de cuerpo rı́gido hacı́a referencia
al hecho de que las distancias inter-particulares antes y después de la consideración
de las fuerzas permanecı́an constantes. Claramente esta hipótesis permitı́a abordar el
problema sin la necesidad de hacer referencia alguna al material constitutivo y permitı́a
tratarlo como una partı́cula con la adición de la interacción rotacional.
En lo que sigue se liberará la hipótesis anterior y se asumirá que como resultado de

las interacciones del sistema con el medio ambiente este puede exhibir cambios en las
distancias intra-particulares. La magnitud de estos cambios dependerá en cada caso de
la naturaleza de las fuerzas inherentes que le dan el carácter denso al sistema, es decir
enlaces electroquı́micos, electromagnéticos, etc. En otras palabras a igual intensidad de
fuerzas aplicadas a cuerpos de igual geometrı́a pero de materiales diferentes los cambios
en distancias entre partı́culas dependerán del material del cuerpo. Los cambios en las
distancias intra-partı́cula se reflejarán posiblemente no solo en cambios de posición sino
también de forma y de tamaño del cuerpo. Uno de los propósitos del modelo de la
mecánica del medio continuo es caracterizar o describir estos cambios a nivel local o del
punto material definido este en el sentido de la ?? en la sección introductoria.
Capı́tulo 2
Análisis de Tensiones
2.1. Introducción
En la sección introductoria se presentaron algunas ideas simples que permitı́an pasar

del problema discreto (modelo de partı́culas) al continuo (modelo de puntos materia-
les). En términos generales, se planteaba la conversión de las ecuaciones gobernantes
representativas de principios fı́sicos propuestos a nivel de la partı́cula a ecuaciones de
campo mediante un proceso de homogenización.
En esta sección se particulariza el problema al caso de la Mecánica Newtoniana para

sistemas de muchas partı́culas, dando como resultado el concepto de tensión el cual
resulta ser una cantidad de tipo tensorial (descrita ahora por 3 vectores indicando la
distribución de fuerzas a lo largo de 3 direcciones).
Al final de esta sección el estudiante deberá estar en la capacidad de:
• Reconocer los tipos de fuerzas presentes en el modelo de los medios continuos.
• Reconocer en el concepto de esfuerzo (o fuerza por unidad de superficie) y en

su caracter tensorial una herramienta matemática para describir los cambios de
cantidad de movimiento en el modelo del medio continuo.
• Entender el vector de tracciones como una descripción de las fuerzas por unidad
de superficie asociadas con un plano determinado.
21
CAPÍTULO 2. ANÁLISIS DE TENSIONES 22
• Determinar las componentes del vector de tracciones en asocio con un plano de-
terminado dado el tensor de esfuerzos en un punto del medio.
• Transformar de un sistema de referencia a otro las componentes del tensor de

esfuerzos.
• Evaluar condiciones de falla en términos de tensiones a través de un punto del

medio continuo.
• Identificar los invariantes del tensor de esfuerzos como las descripciones absolutas
del estado de tensiones en un punto.
• Determinar las tensiones extremas y las direcciones en que estas se presentan para
un estado de tensiones en 2D y 3D.
• Calcular los invariantes para estados de tensiones en 2D y 3D.
• Representar el tensor de tensiones en la orientación principal.
2.2. Interacciones mecánicas-Origen del concepto de

tensión
En el caso particular del problema mecánico, la ecuación (2.1) tomada de la sección

introductoria

dβ
f =α (2.1)
dt
se reduce a la segunda ley de Newton para una partı́cula de masa m, velocidad V~ ,

fuerza F~ y cantidad de movimiento P~ = mV~ . Nuevamente estipulamos que en este caso
la variable respuesta o efecto β corresponde a la cantidad de movimiento P~ = mV~ y la
variable fuente o causa α que introduce cambios temporales en β corresponde a la fuerza
F~ . La ley general de la ecuación (2.1) válida a nivel de la partı́cula se reduce entonces
a la segunda ley de Newton que relaciona los cambios en la cantidad de movimiento y
la fuerza como
dP~ dmV~
F~ = ≡ . (2.2)
dt dt
Esta ley de carácter discreto1 indica que la fuerza es el dispositivo matemático que
permite representar cambios de cantidad de movimiento entre una partı́cula y el medio
ambiente. Claramente identificamos en el concepto matemático de fuerza el punto de
partida del modelo de la Mecánica de los Medios Continuos.
2.2.1. Fuerzas de cuerpo
Nuevamente nos apoyaremos en la esquematización de la figura 2.1 en la cual se consi-

dera un medio infinito que representa el medio ambiente. Este se encuentra acotado
por una superficie S∞ encerrando el volumen B(t). Consideramos además un medio
continuo (el sistema) localizado al interior de este espacio infinito y delimitado por
una superficie S(t) que acota un volumen V (t) y que se encuentra en contacto directo
con el medio ambiente a través de su frontera S(t).
n̂( +)
Región
n̂( −)
externa
S(t)
S∞ Región
B(t ) V (t ) interna
Figura 2.1. Medio continuo interactuando con un espacio infinito.
Utilizaremos la siguiente descripción de las regiones involucradas: La región localizada

del lado positivo del vector normal n̂(+) a S(t) se denominará en lo que sigue la región
externa. De manera consistente, la región localizada del lado negativo del vector normal
n̂(−) a S(t) se denominará la región interna.
1
Discreto significa existente a nivel de la partı́cula
Se desea entonces caracterizar la interacción mecánica entre las regiones interna y ex-
terna debida a la existencia de una fuente de cantidad de movimiento pero tratada
esta como una distribución o función. Se denota allı́ esta función como f~, la cual en lo
que sigue equivale a la función general Φ y representa una densidad de fuerza (fuerza
por unidad de volumen).
Aunque esta densidad se representa en la figura como localizada en un punto se enfatiza

en que la misma puede estar distribuida en todo el volumen incluyendo la región interna
o estar localizada solo en el infinito de la región externa. En cualquier caso esta fuente
es introducida por algún agente externo.
En el primer caso de fuentes internas o una distribución interna de fuerzas por unidad
de volumen decimos que estas son representativas de interacciones a distancia. Esto
significa que el agente externo que las produce no se encuentra en la vecindad del medio
en cuestión. Estos son por ejemplo los campos gravitacionales o electromagnéticos.
2.2.2. Fuerzas de superficie-El vector de tracciones
De cualquier forma, ya sea que existen fuentes internas o externas se generara también
una interacción a través de la superficie o por contacto directo entre los dos medios.
La figura 2.2 representa el diagrama de cuerpo libre de la región interna. En éste se
muestran las interacciones ~t generadas por contacto directo a través de la superficie
S(t) y las interacciones a distancia f~, distribuidas estas sobre todo el volumen V (t).
La función vectorial ~t representa una función o distribución de fuerzas por unidad de
superficie equivalente acá a la interacción por contacto directo entre las regiones interna
y externa. Estas fuerzas en la superficie son necesarias para mantener la condición de
continuidad entre el medio ambiente y el medio continuo.
t
f
t
Figura 2.2. Diagrama de cuerpo libre de la región interna interactuando con la
región externa.
La figura 2.3 muestra el correspondiente diagrama de cuerpo libe para la región externa.
En esta figura, la función ~t nuevamente representa la fuerza por unidad de superficie
sobre la superficie S(t) de la región externa que es igual y opuesta a la función ~t actuando
sobre la superficie S(t) de la región interna (y de acuerdo con el principio de acción y
reacción). Estas fuerzas son equilibradas por una función análoga pero localizada sobre
la frontera o superficie S∞ .
Nuevamente la función ~t ha sido construida con fundamento en el continuo matemático

y su definición de estado lı́mite. En el modelo clásico de la mecánica de los medios
continuos el carácter de esta interacción por contacto directo se encuentra descrito por
el primer postulado de Cauchy. Como tal, este admite verificación experimental más
no demostración matemática. Se llegará entonces a un modelo de la mecánica de los
medios continuos basado en la definición de fuerzas superficiales de Cauchy.
S(t)
S∞ !!
B(t )
Figura 2.3. Diagrama de cuerpo libre de la región externa interactuando con la

región interna.
En algunos tratamientos es común encontrar las funciones f~ y ~t diferenciadas como

“fuerzas” externas e internas. Esto es motivado en el hecho de que las segundas apa-
recen toda vez que propongamos un nuevo diagrama de cuerpo libre haciendo cortes o
exponiendo nuevas superficies extraı́das del sistema original. Por ejemplo si del volumen
V (t) aisláramos un nuevo volumen b(t), sobre la nueva superficie expuesta aparecerı́a
entonces una nueva función ~t que efectivamente mantiene unidas las dos porciones re-
sultantes. En estos términos la distribución ~t era una fuerza interna antes de realizar
la partición y aparece ahora como una fuerza externa en el nuevo diagrama de cuerpo
libre. En cualquier caso la distribución ~t aparece para equilibrar una función fuente
externa f~. En este sentido también consideramos en este tratamiento a la función ~t
como externa.
Es válido entonces clasificar las fuerzas como externas cuando actúan sobre el sistema y
como internas cuando actúan entre partes del sistema. Sin embargo debe aclararse que
mediante una elección apropiada de un diagrama de cuerpo libre extraı́do del sistema
original, cualquier fuerza interna al sistema original puede exponerse como una fuerza
externa.
En este tratamiento denominaremos las fuentes o fuerzas f~ como fuerzas de cuerpo y

las fuerzas ~t como fuerzas por unidad de superficie o tracciones. Tratamos ambas como
externas con la salvedad de que las primeras son debidas a agentes externos a distancia
mientras que las segundas son debidas a contacto directo a través de superficies del
medio continuo.
Dependencia del vector de tracciones en la dirección del vector normal
En lo que sigue se presentará un ejemplo (ver figura 2.4) aplicado a un caso particular
que pretende servir de abrebocas o justificación para la introducción de manera formal
del primer postulado de Cauchy conducente al concepto matemático de tensión.
n̂2
P
W n̂1
V
S
Superficie no deformable.

Figura 2.4. Interacción por contacto en un punto P sobre la frontera de un medio
continuo de volumen V acotado por la superficie S.
En el ejemplo ilustrado en la figura 2.4 se extrae o aı́sla para efectos de su análisis

el volumen acotado internamente por la superficie matemática (abstracta, artificial o
virtual) S. Este volumen a la vez pertenece a un medio general libre de fuerzas en
sus fronteras laterales y soportando un peso W sobre parte de la frontera superior. En
la figura inserta y correspondiente al punto P , se realizan cortes mediante planos con
vectores normales n̂1 y n̂2 con dirección horizontal y vertical respectivamente.
La figura inserta muestra también una aproximación (si se quiere burda) a la distri-
bución de fuerzas por unidad de superficie (fuerzas de contacto o tracciones) sentida
desde dos planos pasando por el punto P y en las cercanı́as de este. Claramente, en un
lı́mite infinitesimal alrededor del punto P el vector normal a la superficie experimenta
un cambio de dirección (de n̂2 a n̂1 ). Este cambio es análogo al que se presenta en las
distribuciones de fuerza sobre los 2 planos pasando por P . Esta conexión sugiere que
la distribución de fuerzas por unidad de superficie (o interacción por contacto directo)
esta relacionada con la dirección del vector normal a la superficie. Esta dependencia de
la interacción a nivel de la superficie con el vector normal a la misma fue propuesta por
Cauchy en el siguiente postulado.
2.2.3. Primer postulado de Cauchy (o Principio de Euler-Cauchy

1822)
En dos regiones externa e interna de un medio continuo de acuerdo con una superficie
de separación S y un vector normal n̂, la región externa ejerce sobre la región interna
una fuerza (interacción) de contacto ∆F~ = ∆F~ (∆S, n̂) en cada punto P que depende
solo de la posición del punto ~xP y de la dirección del vector normal n̂ a S.
Considérese por ejemplo una serie de regiones o elementos de área alrededor de un

punto P de la frontera matemática S tal y como se muestra en la figura 2.5. En esta
se muestran las vistas frontal y transversal de diferentes elementos de área alrededor
del punto P . La lı́nea punteada representa un plano tangencial a la superficie S en
P y cuya dirección normal se representa mediante el vector n̂. Con cada elemento de
área ∆S alrededor del punto P es posible asociar una distribución de fuerzas sobre
la superficie. Denominamos esta distribución como ∆F~ . La misma puede reducirse al
sistema equivalente fuerza-momento resultante (R,~ M ~ R ) sobre el punto P .
n̂
(ΔS , ΔF )3
(ΔS , ΔF )2
(ΔS , ΔF )1 P
MR
MR
MR
Figura 2.5. Vector de tracciones en un punto del medio continuo P .
~ y momento M
Para cada elemento de área considerado, la fuerza R ~ R resultantes forman
un ángulo determinado con la dirección n̂. El postulado de Cauchy propone que si se

realiza
un proceso
continuo de reducción de los elementos de área ∆S, las cantidades
~ ~
R /∆S y MR /∆S tienden a lı́mites definidos. En el primer caso la función resultante
~t se denominará el vector de tracciones o fuerzas por unidad de superficie en
el punto P , ecuación (5.7), mientras que en el segundo la función resultante µ ~ se
denominará el vector de tracciones de par o momentos por unidad de superficie
y el cual asumiremos como cero en el estado lı́mite, ecuación (2.4)2 .
~t(n̂) ≡ ∆F~ dF~

lı́m = . (2.3)
∆S(n̂)→0 ∆S(n̂) dS
~
∆M
~ (n̂) ≡
µ lı́m = 0. (2.4)
∆S(n̂)→0 ∆S(n̂)
Mediante el proceso de lı́mite fue posible reducir la interacción a nivel de la superficie a

puntos (función continua). El vector de tracciones representa entonces una descripción
intensiva de la interacción por contacto directo entre las regiones interna y externa.
Esta interacción debe además satisfacer la tercera ley de Newton o principio de acción
y reacción.
Como consecuencia nos encontramos con que este vector de tracciones es igual (en
magnitud) y opuesto (en dirección) a la representación de la interacción sobre la región
externa. Esta condición se expresa mediante la relación de equilibrio de la ecuación (2.5)
para el punto P , figura 2.6, y el diagrama de cuerpo libre de la figura 2.3.
~t(n̂+ ) + ~t(n̂− ) = 0 (2.5)

2
El asumir este lı́mite como cero da origen a un modelo clásico de la mecánica de los medios
continuos. Podrı́a ser perfectamente válido asumir este lı́mite como diferente de cero.
−
P
t ( nˆ )
n̂( −)
n̂( +)
+
t ( nˆ )
P
Figura 2.6. Transmisión de tracciones a ambos lados de una superficie pasando por
el punto P .
2.2.4. Fuerzas por unidad de superficie sobre un plano con

dirección normal n̂
Consideremos ahora el elemento diferencial de la figura 2.7 el cual se ha construido

alrededor de un punto material P . En el elemento se han realizado cortes a través de
los 3 planos perpendiculares a los ejes cartesianos y con vectores normales −î, −ĵ y −k̂
y a través de un plano inclinado con vector normal n̂ = nx î + ny ĵ + nz k̂. Sobre este
plano también se traza un vector ν̂ = νx î + νy ĵ + νz k̂ paralelo al plano. Los vectores
de tracción sobre los planos cartesianos están dados por ~t(−î) , ~t(−ĵ) y ~t(−k̂) , mientras
que el vector de tracciones sobre el plano inclinado está dado por ~t(n) . Si se sabe que el
plano de corte tiene un vector de superficie dS~ = dS n̂, donde dS es el área de la sección
superficial, se desea plantear las ecuaciones de equilibrio en las direcciones n̂ y ν̂.
Para resolver el problema será necesario convertir los vectores de tracción sobre cada
una de las caras a vectores fuerza y posteriormente proyectar las fuerzas resultantes en
las direcciones n̂ y ν̂. Denotando las componentes escalares del vector ~t(n) a lo largo de
las direcciones n̂ y ν̂ como σnn y τnν respectivamente, este puede escribirse como;
~t(n̂) = σnn n̂ + τnν ν̂. (2.6)

k̂ (nˆ )
t
µˆ ( -iˆ )
t n̂
τnµ
σ nn
( -jˆ )
t ĵ
ˆ)
iˆ t
( -k
Figura 2.7. Elemento diferencial cortado con un plano de vector normal n̂ tras lo
cual se expone el vector de tracciones ~t(n) .
P
Planteando primero Fn = 0 se tiene;
σnn dS − [~t(−î) · n̂]dSx − [~t(−ĵ) · n̂]dSy − [~t(−k̂) · n̂]dSz = 0.
Ahora, expresando los vectores de tracción ~t(−î) , ~t(−ĵ) y ~t(−k̂) en términos de sus com-
ponentes escalares en las direcciones î, ĵ y k̂ es posible escribir la expresión anterior
como:
σnn dS − (σxx nx + τxy ny + τxz nz )dSx − (τyx nx + σyy ny

+τyz nz )dSy − (τzx nx + τzy ny + σzz nz )dSz = 0
finalmente, escribiendo las áreas de los planos cartesianos dSx , dSy y dSz en términos
del área dS del plano inclinado de acuerdo con:
dSx = dSnx
dSy = dSny
dSz = dSnz
se tiene que la fuerza por unidad de superficie en la dirección normal al plano de corte
está dada por:
σnn = σxx n2x + τxy nx ny + τxz nx nz + τyx nx ny

(2.7)
+σyy n2y + τyz nz ny + τzx nx nz + τzy ny nz + σzz n2z .
P
Procediendo de manera similar planteamos Fν = 0 como;
τnν dS − [~t(−î) · ν̂]dSx − [~t(−ĵ) · ν̂]dSy − [~t(−k̂) · ν̂]dSz = 0.
y expresando nuevamente los vectores de tracción en términos de sus componentes

escalares obtenemos;
τnν dS − (σxx νx + τxy νy + τxz νz )dSx − (τyx νx + σyy νy

+τyz νz )dSy − (τzx νx + τzy νy + σzz νz )dSz = 0.
Luego de usar nuevamente las expresiones para las áreas de las caras asociadas a los
planos cartesianos en términos del área del plano inclinado de acuerdo con:
dSx = dSnx
dSy = dSny
dSz = dSnz
se tiene que la fuerza por unidad de superficie en la dirección paralela al plano de corte
está dada por:
τnν = σxx νx nx + τxy νy nx + τxz νz nx + τyx νx ny

(2.8)
+σyy νy ny + τyz νz ny + τzx νx nz + τzy νy nz + σzz νz nz .
2.2.5. Tensor de tensiones
De las ecuación (2.7) y ecuación (2.8) es posible concluir que para determinar comple-
tamente el vector de tracciones ~t(n) en el punto material P sobre una superficie con
dirección normal n̂ es necesario conocer los vectores de tracción ~t(−î) , ~t(−ĵ) y ~t(−k̂)
asociados con 3 direcciones perpendiculares.3
Estos 3 vectores de tracción (o equivalentemente sus respectivas componentes escalares)

forman una cantidad tensorial (tensor de orden 2) denominada el tensor de tensiones.
Estos permiten definir de manera completa el estado de fuerzas por unidad de superficie
en asocio con cualquier dirección de un punto material del medio continuo. Es común
almacenar o representar este tensor en una matriz de 3 × 3 en la cual cada vector fila
almacena las componentes escalares de los vectores de tracción.
 
σxx τxy τxz
 τyx σyy τyz 
τzx τzy σzz
La figura 2.8 muestra estos vectores y sus componentes escalares. En dicha figura y en
la representación matricial del tensor de tensiones los ı́ndices i y j en un término tı́pico
σij hacen referencia al plano sobre el que está aplicado el vector de tracciones y a la
dirección sobre la cual se está proyectando el mismo.
3
Es posible demostrar que en realidad basta con conocer los 3 vectores de tracción sobre planos aso-
ciados con 3 direcciones no co-lineales y no necesariamente ortogonales. Sin embargo por conveniencia
se eligen las 3 direcciones cartesianas.
(kˆ)
z k̂
t
ˆ
( j)
t
z
ĵ
σxx
(iˆ) y
t σ zz τ xy
î
τ zy
τ xz
τ zx τ yz
x τ yx
σ yy τ xz σ yy
y
τ xy τ yx
τ yz
σ xx
τ zx
τ zy
x σ zz
Figura 2.8. Vectores de tracción actuando sobre los planos de un elemento diferencia
en un sistema de referencia Cartesiano.
En esta representación también se han utilizado los simbolos σ y τ para diferenciar las
componentes normales (almacenadas en la diagonal de la matriz) y paralelas al plano
(almacenadas por fuera de la diagonal de la matriz). Por ejemplo la componente τxy hace
referencia a la proyección en la dirección ĵ (ı́ndice y) del vector de tracciones asociado
al plano con normal î (ı́ndice x). Similarmente, la componente τzx hace referencia a la
proyección en la dirección î (ı́ndice x) del vector de tracciones asociado al plano con
normal k̂ (ı́ndice z).
Cálculo del vector de tracciones sobre un plano con vector normal n̂ (Formula
de Cauchy).
En la sección anterior se derivaron las ecuación (2.7) y ecuación (2.8) las cuales per-
mitı́an calcular las componentes σnn y τnν que definı́an el vector de tracciones ~t(n) en el
plano inclinado con vector normal n̂, mediante las 9 componentes escalares del tensor
de tensiones (o esfuerzos). Nótese que estas componentes permiten expresar el vector
~t(n) en el sistema de referencia n̂ − ν̂ de acuerdo con:
~t(n̂) = σnn n̂ + τnν ν̂.
Aunque desde el punto de vista fı́sico las componentes σnn y τnν son bastante útiles
toda vez que dan cuenta de lo que experimenta el medio a lo largo del plano en térmios
de tensiones normales y paralelas a este, en ocasiones resulta útil o inclusive necesario el
determinar directamente el vector ~t(n) en el sistema de referencia x − y − z en el cual se
conoce el tensor de tensiones σij . Este cálculo corresponde simplemente a la proyección
en la dirección n̂ de cada uno de los vectores ~t(î) , ~t(ĵ) y ~t(k̂) que definen el estado de
tensiones en el punto P como se ilustra a continuación.
(n̂) (n̂)
Primero escribimos el vector ~t(n) en términos de sus componentes escalares tx , ty y
(n̂)
tz como:
~t(n̂) = t(n̂) (n̂) (n̂)

x î + ty ĵ + tz k̂.
Estas componentes se obtienen tras plantear y resolver las ecuaciones de equilibrio del
elemento mostrado en la figura 2.7 en las direcciones î, ĵ y k̂ respectivamente. Las
ecuaciones resultantes pueden escribirse de la forma:
t(n̂)
x = t
~(î) · n̂
t(n̂)
y = t
~(ĵ) · n̂
t(n̂)
z = t
~(k̂) · n̂
o en la siguiente ecuación matricial:
 
(n̂)   
 tx
 
 σxx τxy τxz  nx 
(n̂)
ty =  τyx σyy τyz  ny . (2.9)
 t(n̂)
 
 τzx τzy σzz

nz

z
La ecuación (2.9) correspondiente al producto entre el tensor de tensiones y el vector

normal a una plano (o superficie de corte) se conoce como la fórmula de Cauchy. Ma-
temáticamente ésta no es otra cosa que la proyección del tensor σij a lo largo de la
dirección n̂. De manera condensada la ecuación (2.9) se puede escribir como:
t(n) = [σ] {n}

(2.10)
Ejemplo 1: Tracciones sobre un plano arbitrario a partir de σij
El estado de tensiones sobre un punto material P de un medio continuo está dado por
σij . Determinar el vector de tracciones sobre un plano T definido por los puntos A, B
y C. Adicionalmente determinar las componentes de la tensión en la dirección normal
y paralela al plano.
Solución
Primero necesitamos determinar la ecuación del plano definido por los puntos A, B y
C y la dirección del vector normal. Recodando que la ecuación del plano esta dada por
ax + by + cz = d
donde los coeficientes a, b y c son las componentes del vector normal n̂. Se tiene entonces
que:
~rAC × ~rAB
n̂ =
k~rAC × ~rAB k
luego:
a ≡ nx
b ≡ ny
c ≡ nz
mientras que d = axA + byA + czA . Ahora es posible determinar el vector de tracciones
sobre una superficie con vector normal n̂ usando la ecuación (2.9):
 
(n̂)   
 tx
 
 σxx τxy τxz  nx 
(n̂)
ty =  τyx σyy τyz  ny .
 t(n̂)
 
 τzx τzy σzz

nz

z
Para determinar las componentes normal y paralela se proyecta el vector ~t(n̂) a lo largo
(n̂) (n̂) (n̂)
de n̂. Si se denotan las componentes de ~t(n̂) por tx , ty y tz se tiene que:
σnn = ~t(n̂) · n̂ ≡ t(n̂) (n̂) (n̂)

x nx + ty ny + tz nz
mientras que la tensión paralela al plano de corte se determina como:
2
τnν = ~t(n̂) · ~t(n̂) − σnn
2
2.2.6. Transformación por rotación. Caso general
Al igual que en el caso vectorial, en varias aplicaciones será conveniente conocer la

descripción del tensor de tensiones bajo la rotación de su sistema de referencia, o más
interesante aún, deseamos saber como varı́an las tensiones en las diferentes direcciones
de un punto material del medio continuo. Este problema corresponde al problema de
transformación bajo rotación o la determinación de las componentes del tensor de
tensiones en un sistema de referencia primado x0 − y 0 − z 0 a partir de las componentes
en un sistema de referencia no primado x − y − z, ver figura 2.9.
z
z
z'
σ zz
z'z
'
τ z'y' y'
σ τ y'z' y'y
'
τ zy σ
τ zx τ yz
τ z'x'
τ xz σ yy τ x'z' τ y'x'
y τ x'y' y
τ xy τ yx
σ xx '
σ x'x
x x
x'
Figura 2.9. Representación de tensor de tensiones en un sistema x − y − z y en un

sistema alternativo x0 − y 0 − z 0 rotado con respecto al primero.
En la figura se hace la descripción esquemática del tensor en el sistema primado y no

primado respectivamente. En la sección anterior se resolvió el problema de la determina-
ción del vector de tracciones en asocio con una dirección normal n̂ (ver ecuación (2.9)).
Para resolver el problema de transformación bajo rotación basta con reconocer las di-
recciones iˆ0 , jˆ0 y k̂ 0 del nuevo sistema de referencia como vectores normales a lo largo
0 0 0
de los cuales se desea determinar los vectores de tracciones ~t(i ) , ~t(j ) y ~t(k ) . En la fi-
gura 2.10 se presentan en lı́nea púrpura las direcciones asociadas con el sistema de
referencia x − y − z, mientras que en lı́nea negra se presentan las
direcciones asociadas con el sistema de referencia primado x0 − y 0 − z 0 . Denotaremos

por σ al tensor de tensiones en el primer sistema de referencia y por σ 0 el del segundo.
Figura 2.10. Representación de tensor de tensiones en sistema x−y −z y x0 −y 0 −z 0 .
Utilizando la ecuación (2.9), o equivalentemente la ecuación (2.10) determinemos el

vector de tracciones en caras con vectores normales iˆ0 , jˆ0 y k̂ 0 de acuerdo con:
0
{t(i ) } = [σ] {i0 }
0
{t(j ) } = [σ] {j 0 } (2.11)
(k0 )
{t } = [σ] {k 0 }
Ahora, estos vectores deben ser expresados en términos de componentes escalares en el
sistema de referencia x0 − y 0 − z 0 para lo cual proyectamos cada uno de ellos a lo largo
de las direcciones iˆ0 , jˆ0 y k̂ 0 como sigue
~t(i0 ) = [~t(i0 ) · î0 ]î0 + [~t(i0 ) · ĵ 0 ]ĵ 0 + [~t(i0 ) · k̂ 0 ]k̂ 0
~t(j 0 ) = [~t(j 0 ) · î0 ]î0 + [~t(j 0 ) · ĵ 0 ]ĵ 0 + [~t(j 0 ) · k̂ 0 ]k̂ 0 (2.12)
(k0 ) (k0 ) (k0 ) (k0 )
~t = [~t · î0 ]î0 + [~t · ĵ 0 ]ĵ 0 + [~t · k̂ 0 ]k̂ 0
Estas expresiones equivalen a:

~t(i0 ) = σx0 x0 î0 + τx0 y0 ĵ 0 + τx0 z0 k̂ 0
~t(j 0 ) = τy0 x0 î0 + σy0 y0 ĵ 0 + τy0 z0 k̂ 0
~t(k0 ) = τz0 x0 î0 + τz0 y0 ĵ 0 + σz0 z0 k̂ 0 .
Finalmente si escribimos los vectores iˆ0 , jˆ0 y k̂ 0 en el sistema de referencia x − y − z en

términos de sus cosenos directores como:
î0 = Cx0 x î + Cx0 y ĵ + Cx0 z k̂

ĵ 0 = Cy0 x î + Cy0 y ĵ + Cy0 z k̂ (2.13)
k̂ 0 = Cz0 x î + Cz0 y ĵ + Cz0 z k̂
Substituyendo las ecuación (2.11) y ecuación (2.13) en la ecuación (2.12) se tiene:
[σ 0 ] = [C][σ][C]T (2.14)
donde [C], que es conocida como la matriz de cosenos que relacionan los dos sistemas
y está dada por:
 
C x0 x C x0 y C x0 z
[C] = Cy0 x
 Cy 0 y Cy 0 z 
Cz 0 x Cz 0 y Cz 0 z
2.2.7. Transformación por rotación 2D. Cı́rculo de Mohr
En muchas aplicaciones prácticas los problemas de tensiones pueden ser estudiados

mediante un análisis en dos dimensiones. Cuando éste es el caso la transformación por
rotación puede ser desarrollada a partir de una herramienta gráfica conocida como el
cı́rculo de Mohr. En lo que sigue se deducen las expresiones para pasar de un tensor
con representación en un sistema x − y a su respresentación en un sistema primado
x0 − y 0 . En la figura 2.11 se muestra la descripción esquemática del tensor en el sistema
primado y no primado respectivamente.
Figura 2.11. Izquierda: Tensor [σ] en el sistema de referencia x − y, Derecha: Tensor

[σ 0 ] en el sistema de referencia x0 − y 0 . El sistema de referencia x0 − y 0 se encuentra
rotado un ángulo θ respecto a x − y.
Desarrollando de la forma expandida la ley de transformación por rotación presentada

en ecuación (2.14), para el caso 2D se tiene:
0 0

σxx τxy Cθ Sθ σxx τxy Cθ −Sθ
0 0 = (2.15)
τxy σyy −Sθ Cθ τxy σyy Sθ Cθ
Desarrollando la multiplicación matricial y escribiendo las tres expresiones que definen

0 0 0
los tres escalares σxx , σyy y τxy :
0
σxx = σxx Cθ2 + σyy Sθ2 + 2τxy Sθ Cθ
0
σyy = σxx Sθ2 + σyy Cθ2 − 2τxy Sθ Cθ (2.16)
0
τxy = (σyy − σxx )Sθ Cθ + τxy (Cθ2 − Sθ2 )
Sumando la primera y la segunda en la expresión anterior se tiene:

0 0
σxx + σyy = σxx + σyy (2.17)
La ecuación (2.17) muestra que la suma de las componentes normales (elementos sobre
la diagonal) es invariante (no cambia) bajo la rotación del sistema de referencia. Elimi-
0
nando σyy , se tiene:
0
σxx − (σxx + σyy ) = −(σxx Sθ2 + σyy Cθ2 − 2τxy Sθ Cθ ). (2.18)
Utilizando las relaciones trigonométricas:
1
sen2 θ = [1 − cos(2θ)]
2
1
cos2 θ = [1 + cos(2θ)]
2
cos θ − sen2 θ = cos(2θ)
2
1
sen θ cos θ = sen(2θ)
2
la ecuación (2.18) se puede reescribir como:

0 1 − C2θ 1 + C2θ
σxx − (σxx + σyy ) = − σxx + σyy − τxy S2θ
2 2
o agrupando términos:
0 σxx + σyy σxx − σyy

σxx − = C2θ + τxy S2θ (2.19)
2 2
Reemplazando las identidades la tercer expresión de la ecuación (2.15) quedarı́a:
0 σxx − σyy
τxy =− S2θ + τxy C2θ (2.20)
2
Para terminar se suman los cuadrados de las ecuación (2.19) y ecuación (2.20)
2 2
0 σxx + σyy 0 2
σxx − σyy 2
σxx − + τxy = + τxy (2.21)
2 2
Comparando la ecuación (2.21) con la ecuación de un cı́rculo, r2 = x2 + y 2 , se puede

concluir que ésta representa la ecuación de un cı́rculo con radio R y centro C en el

sistema de referencia σ − τ .
2
2 σxx − σyy 2
R = + τxy (2.22)
2
(2.23)

σxx + σyy
C= ,0 (2.24)
2
Reescribiendo la ecuación (2.21) se tiene una ecuación describe todos los vectores de
tracción asociados con las infinitas direcciones del espacio para un tensor dado de acuer-
do con:
0 20
2
(σxx − C) + τxy = R2 (2.25)
Las componentes del tensor en el espacio σ vs τ corresponden a 2 puntos A y B dados
por;
A = (σxx , −τxy )
(2.26)
B = (σyy , +τxy )
El punto A corresponde a la componente vectorial del tensor asociada con la dirección
î mientras que el punto B corresponde a la componente vectorial del tensor asociada
con la dirección ĵ . Nótese que en la ecuación (2.26) el signo (−) que antecede a la
componente τxy del punto A es artificial. Esto equivale a distorsionar la componente
vectorial asociada con la dirección î. En la figura 2.12 se representan los puntos A y B
sobre el cı́rculo. En este caso se asume que σyy > σxx y τxy > 0.
Figura 2.12. Cı́rculo de Mohr para para σyy > σxx y τxy > 0.
En conclusión, dado el tensor de tensiones σij en el sistema de referencia x − y es

posible trazar el cı́rculo correspondente al estado de esfuerzos en el punto mataerial y a
partir de este se pueden determinar sus componentes en cualquier otra dirección. Esta
aparece en el cı́rculo con una orientación 2θ con respecto a la dirección original debido
a la parametrización de las ecuaciones en la derivación misma de la ecuación del cı́rculo.
En consecuencia, un angulo de θ en el problema fı́sico aparece como 2θ en el cı́rculo de
Mohr.
Un problema de especial interés, y facilmente soluble con el cı́rculo de Mohr, es el de

determinar las direcciones en las cuales se presentan las máximas y mı́nimas tensiones
normales y paralelas. Estas aparecen en el cı́rculo de la figura 2.12 rotuladas como
λ1 = σmáx , λ2 = σmı́n y τ = τmáx .
Ejemplo 2: Cı́rculo de Mohr
La figura 2.13 muestra dos estados de tensiones para el mismo punto de un medio
continuo.
Y
10 MPa Y
10 MPa
10 MPa
50 MPa 20 MPa
=30º X =60º X
(a) Estado de tensión 1 (b) Estado de tensión 2
Figura 2.13. Estado de tensiones en un punto.
Para cada estado de carga mostrado en la figura 2.13a y en la figura 2.13b, respectiva-
mente.
• Determinar los esfuerzos normales y cortantes sobre los planos mostrados con li-
nea punteada.
Solución:
Inicialmente determinemos el tensor de tensiones para los estados 1 y 2 mostrados

en la figura 2.13a y la figura 2.13b respectivamente.

0.0 −50.0 20.0 −10.0
[σ1 ] = [σ2 ] =
−50.0 10.0 −10.0 −10.0
Aunque hay diversos caminos para llegar a una solución, para dar una mayor
ilustración del manejo del cı́rculo de Mohr haremos uso de esta herramienta 4 .
En la figura 2.14a y en la figura 2.14b se muestra el cı́rculo de Mohr para ambos
estados de tensiones. Dicho estado se indica mediante la lı́nea X − Y .
4
Es buen ejercicio cotejar estos resultados haciendo uso del manejo matricial de la transformación.
X X'
A'
X
X'
A' 60º 60º
a a
B'
Y' Y
B'
Y'
Y
Figura 2.14. Cı́rculo de Mohr con tensor en sistema original x − y y sistema rotado
x0 − y 0
Para encontrar las tensiones normales y tangenciales solicitadas hacemos una ro-
tación de 60◦ en el cı́rculo de Mohr (30◦ en el problema fı́sico) en sentido contrario
a las manecillas del reloj. Esta rotación se indica en el cı́rculo de Mohr de la figu-
ra 2.14 mediante las lı́neas punteadas color rojo X 0 − Y 0 que unen los puntos A0
y B 0 . Si se lee la información de los puntos A0 y B 0 se obtienen los tensores en el
sismema x0 − y 0 . El estado de tensiones en el sistema rotado serı́a entonces dado
por:
0
τxy = R sen(β)
0
σxx = C − R cos(β) (2.27)
0
σyy = C + R cos(β)
Donde R y C son el radio y el centro del cı́rculo de Mohr respectivamente.

◦ 50
Para el estado 1 R = 50.25, C = 5.0 y β = α − 60 . α = arctan
5
◦ ◦ 10
Para el estado 2 R = 18.03, C = 5.0 y β = 180 − (α + 60 ). α = arctan
15
Los tensores en el sistema rotado serı́an:

0 −40.80 −20.67 3.83 −17.99
[σ1 ] = [σ20 ] =
−20.67 50.80 −17.99 6.16
A partir de los tensores anteriores se concluye que las tensiones normal y tan-
gencial al plano solicitado son σ = −40.80 y τ = −20.67 para el estado 1 y
σ = 3.83 y τ = −17.99 para el estado 2 respectivamente. En la figura 2.15a y en
la figura 2.15b se representa el cubo de tensiones para cada estado.
Y Y
Y' Y'
50.80 6.16
20.67 40.80 17.99 3.83
X' X'
20.67 17.99
30º 30º
X X
Figura 2.15. Cubo de Tensiones en sistema rotado x0 − y 0
• Para cada uno de los estados de tensiones mostrados determinar los valores pro-
pios y las direcciones principales.
Solución:
Los valores principales no son otra cosa que los valores máximos axiales (nor-
males) a los cuales está sometido el punto. Este estado se logra precisamente
cuando el valor de las tensiones tangenciales es nulo. En la figura 2.16a y en la
figura 2.16b se muestra en el cı́rculo de Mohr, el sistema original x−y y el sistema
principal xp − yp , mientras que en la figura 2.17a y en la figura 2.17b se muestra
el cubo de tensiones para el estado principal.
X
84.29º
Ap Bp
33.69º Ap Xp
Xp Yp Bp
Yp
Figura 2.16. Cı́rculo de Mohr con valores máximos normales. Puntos Ap , Bp .
Y 16.84º
Y Yp
Yp
55.05 45.05
13.03
Xp
42.14º 30º
X
X 23.03
Xp
Figura 2.17. Cubo de tensiones en sistema rotado xp − yp .
Nótese que para encontrar el estado

principal 1 rotamos desde el sistema original
50
xy un ángulo α = arctan en sentido contrario a las manecillas del reloj,
5
10
mientras que para el estado principal 2 rotamos un ángulo α = arctan en
15
sentido de las manecillas 5 . Por similitud a la ecuación (2.27) las ecuaciones para
5
La rotación puede ser hecha en cualquier sentido; lo importante es no perder de vista la ubicación
encontrar los valores se podrı́an escribir como:

p
τxy =R
p
σxx =C ±R (2.28)
p
σyy = C ∓ R
Los estados de tensiones principales quedarı́an escritos como:

p −45.25 0.0 p 23.03 0.0
[σ1 ] = [σ2 ] = (2.29)
0.0 55.25 0.0 −13.03
• Asuma que ambos estados de tensiones ocurren al mismo tiempo y determine las
direcciones principales y valores propios para el estado de tensiones resultantes.
Solución: Para hacer superposición de los estados de tensiones ambos deben

de estar en el mismo sistema de referencia. En nuestro caso se tiene que inicia-
lemente ambos están escritos en el mismo sistema de referencia xy, por lo que
bastará hacer la suma componente a componente de ambos tensores.

0.0 −50.0 20.0 −10.0 20.0 −60.0
[σt ] = + = (2.30)
−50.0 10.0 −10.0 −10.0 −60.0 0.0
En la figura 2.18a se presenta el cı́rculo de Mohr con el estado original producto

de la superposición (linea X − Y ) y el estado principal (linea Xp − Yp ) mientras
que en la figura 2.18b se muestra el cubo de tensiones del estado principal.
de los ejes y tener en cuenta que estas rotaciones coinciden en sentido con el problema fı́sico pero que
corresponden a valores de ángulo doble en el cı́rculo
Y
40.26º
X Yp
50.83
80.53º
Bp Ap
Yp Xp
X
70.83
Xp
Y
(a) Circulo de Mohr (b) Cubo de tensiones
Figura 2.18. Estado de tension principal superpuesto
El tensor escrito en direcciones principales es:

p 70.83 0.0
[σ2 ] = (2.31)
0.0 −50.83
Ejemplo 3: Planos de falla
Un punto material en un medio continuo está sometido al estado de tensiones mostrado

en la figura 2.19. Si las lı́neas punteadas corresponden a planos de falla del material y
sabiendo que los dos planos son perpendiculares entre sı́, se pide:
(a)
Figura 2.19. Estado de tensiones
• Calcular las tensiones normal y tangencial al plano 1.
Solución:
Escribiendo el estado de tensiones como:

T 0
[σ] =
0 0
Para el plano 1 el vector normal es nˆ1 = [− sen θ, cos θ], y la tensión en el plano
se calcula como:

T 0 − sen θ
t1 = [σ]nˆ1 =
0 0 cos θ

−T sen θ
t1 =
0
La tensión normal se calcula proyectando la tensión t1 sobre el vector nˆ1 ,


−T sen θ − sen θ
σn = t1 · nˆ1 = · ≡ T sen2 θ
0 cos θ
El esfuerzo tangencial se puede calcular proyectando la tensión t1 sobre un vector
unitario inscrito en el plano 1. Dando como resultado:
τ = T sen θ cos θ
• Calcular las tensiones normal y tangencial al plano 2.
Solución:
Siguiendo un proceso análogo al anterior, se llega a los siguientes resultados:
σn = T cos2 θ
τ = T sen θ cos θ
• Suponiendo que la falla va a ocurrir en alguno o en ambos planos, y que el ma-

terial es infinitamente resistente a esfuerzos normales, pero no soporta esfuerzos
cortantes mayores a τadm = T /4, calcule el ángulo θ al cual se produce la falla e
indique en cuál plano se presenta la misma.
Solución:
Graficando el estado de tensiones como un cı́rculo de Mohr (figura 2.20), se ubica

τadm = T /4. Luego se llega a la relación:
T /4 1 π
sen 2θ = = →θ=
T /2 2 12
Figura 2.20. Circulo de Mohr correspondiente al Estado de tensiones de la figura 2.19
La falla se da en cualquier plano, ya que el valor del cortante es igual es ambos

por ser perpendiculares.
Ejemplo 4: Refuerzo
Un punto material de un medio continuo se encuentra sometido al estado de tensiones

correpondiente al cı́rculo de Mohr de la figura 2.21 y en el cual τR es la resistencia al
corte del material.
Si se requiere aplicar una tracción S = 3τR , proponga una medida de mejoramiento del
material, en términos de una aplicación adicional de tensiones, que garantice que éste
no falle por cortante.
Del cı́rculo de Mohr, se puede ver que la única manera de evitar que el material falle,
es llevándolo a un estado de esfuerzos que no alcance τR . En otras palabras, se necesita
reducir el radio del cı́rculo de Mohr. La reducción de radio necesaria se muestra en la
figura 2.22.
Figura 2.21. Estado de tensiones
(a)
Figura 2.22. Cambio en el estado de tensiones
De la figura 2.22 se ve el radio inicial R = S/2. El radio del nuevo cı́rculo debe ser
R = S/3 para que el material no falle por cortante. Entonces
S σxx − σyy S
R= = → σyy =
3 2 3
Figura 2.23. Cambio en el estado de tensiones
En la figura 2.23 se ilustra le cambio en las configuraciones de tensiones.
Ejemplo 5: Problema de Cauchy usando Python
En el enlace: https://github.com/AppliedMechanics-EAFIT/Notas-MMC/blob/master/
notebooks/Ej3_Cauchy.ipynb
2.3. Valores y direcciones principales para un esta-

do de tensiones en 3D
En la sección anterior, para el caso de estados de tensiones en 2 dimensiones (2D),

identificamos direcciones o planos trazados por un punto material para los cuales el
vector de tracciones solo tenı́a componente en la dirección normal al plano. Mas aún,
se demostró que estas componenetes correspondı́an a las tensiones máximas y mı́nimas
que experimentaba el punto material. En esta sección se aborda el mismo problema
para el caso de estados de tensiones en 3D.
k̂ (nˆ )
t
µˆ ( -iˆ )
t n̂
τnµ
σ nn
( -jˆ )
t ĵ
ˆ)
iˆ t
( -k
Figura 2.24. Elemento diferencial cortado con un plano de vector normal n̂ tras lo
cual se expone el vector de tracciones ~t(n) .
Para esto consideremos nuevamente el tetrahedro de Cauchy (figura 2.24) representan-

do el estado de tensiones en un punto material del medio continuo. Sobre este se traza
un plano con vector normal externo n̂. Adicionalmente, sobre este plano también tra-
zamos un vector µ̂ contenido en el plano inclinado. De otro lado, el vector de tracciones
asociado con el plano inclinado se denota como ~t(n̂) . Nótese que el vector de tracciones
y el vector normal n̂ no son colineales.
En la ecuación (2.32) el vector de tracciones se escribe en términos de sus componentes

normal y paralela al plano, denotadas como σnn y τnµ respectivamente.
~t(n̂) = σnn n̂ + τnν ν̂. (2.32)
De la representación anterior claramente se identifica que la máxima componente normal

σnn se presenta cuando las direcciones de ~t(n̂) y n̂ coincidan, en cuyo caso podemos
escribir;
~t(n̂) = σ̃n̂ (2.33)
donde σ̃ representa dicho valor máximo. De otro lado, tenemos que el vector ~t(n̂) también
puede ser encontrado a partir del tensor de tensiones y de la dirección n̂ a través de la

fórmula de Cauchy escrita en notación matricial como;
(n̂)
t = [σ] {n̂} . (2.34)
Comparando la ecuación (2.33) con la ecuación (2.34) podemos escribir;
[σ − σ̃I] {n̂} = {0} (2.35)
el cual corresponde a un problema de autovalores para la matriz de esfuerzos. Este

problema tiene una solución no trivial (n̂ 6= 0) cuando;
det [σ − σ̃I] = 0.
En foma expandida escribimos;
 
σxx − σ̃ τxy τxz
det  τyx σyy − σ̃ τyz  = 0
τzx τzy σzz − σ̃
y evaluando el determinante se tiene;
σ̃ 3 − (σxx + σyy + σzz )σ̃ 2 + (σxx σyy + σxx σzz + σzz σyy − τxy
2 2
− τxz 2
− τyz )σ̃−
2 2 2
(σxx σyy σzz + 2τxy τxz τyz − σxx τyz − σyy τxz − σzz τxy ) = 0. (2.36)
• La ecuación (2.36) corresponde a la ecuación caracterı́stica del problema de au-

tovalores dado en ecuación (2.35).
• Dado que el tensor de tensiones es simétrico y real su ecuación caracterı́stica tiene

3 raices reales distintas denotadas como σ̃1 , σ̃2 y σ̃3 y tales que σ̃1 ≥ σ̃2 ≥ σ̃3 .
• Los valores principales σ̃1 , σ̃2 y σ̃3 corresponden a las tensiones máxima, interme-
dia y mı́nima que experimenta el punto material.
• Para cada raı́z (tensión extrema) de la ecuación caracterı́stica existe una dirección
asociada n̂ la cual se encuentra usando la ecuación (2.35) y la condición n2x + n2y +
n2z = 1.
De otro lado, considerando que le ecuación (2.36) da como solución las tensiones prin-
cipales σ̃1 , σ̃2 y σ̃3 necesariamente se tiene que los coeficientes deben ser independientes
del sistema de referencia en que se tenga el tensor de tensiones. De lo anterior se tiene
que la ecuación caracterı́stica se puede escribir como;
σ̃ 3 − Iσ σ̃ 2 + IIσ σ̃ − IIIσ = 0
donde;
Iσ = σxx + σyy + σzz

2 2 2
IIσ = σxx σyy + σxx σzz + σzz σyy − τxy − τxz − τyz
2 2 2
IIIσ = σxx σyy σzz + 2τxy τxz τyz − σxx τyz − σyy τxz − σzz τxy
son constantes independientes del sistema de referencia6 .
Consideremos el siguiente estado de tensiones con componentes correspondientes a un

sistema de referencia x − y − z;
 
200 100 300
σ =  100 0 0 .
300 0 0
Se pide determinar las tensiones principales y sus direcciones asociadas.
La ecuación caracterı́stica esta dada por;
σ̃ 3 − 200σ̃ 2 − 100000σ̃ = 0
6
Estas constantes se denominan los invariantes del tensor de tensiones y juegan un papel análogo
a la magnitud de un vector la cual permanece constante independientemente del sistema de referencia
en el que se escriba el vector
donde identificamos los 3 invariantes como;
Iσ = 200
IIσ = −100000
IIIσ = 0.
Resolviendo la ecuación caracterı́stica determinamos las tensiones principales σ̃1 = 432,
σ̃2 = 0 y σ̃3 = −232.
La dirección asociada con la tesion σ̃1 = 432 resulta de resolver;
    
−232 100 300  nx   0 
 100 −432 0  ny = 0
300 0 −432 nz 0
   
n2x + n2y + n2z = 1.
De la segunda y tercera ecuación se tiene;
100
ny = nx
432
300
nz = nx .
432
Usando estas en n2x + n2y + n2z = 1 se tiene que;
nx = ±0.807
mientras que ny = 0.187 y nz = 0.560. Por lo tanto la dirección asociada con la tensión
principal σ̃1 = 432 corresponde al vector n̂1 = 0.807î + 0.187ĵ + 0.560k̂.
Similarmente, resolviendo para σ̃2 = 0 se tiene;
    
200 100 300  nx   0 
 100 0 0  ny = 0
300 0 0 nz 0
   
De la segunda (o la tercera) nx = 0 mientras que de la primera ny = −3nz . Finalmente

de la condición n2x + n2y + n2z = 1 se tiene que la dirección asociada con la tensión
principal σ̃2 = 0 corresponde a n̂2 = −948ĵ + 0.316k̂. Finalmente para la tensión
principal σ̃3 = −232 resolvemos;
    
432 100 300  nx   0 
 100 232 0  ny = 0 .
300 0 232 nz 0
   
100 300
De la segunda y la tercera se tiene que ny = − 231 nx y nz = − 231 nx mientras que de la
condición n2x + n2y + n2z = 1 se tiene que la dirección asociada con la tensión principal
σ̃3 = −232 corresponde a n̂3 = 0.592î − 0.255ĵ − 0.766k̂.
Nótese que en el paso final del cálculo de las direcciones principales n̂1 , n̂2 y n̂3 se uso
la condicı́on n2x + n2y + n2z = 1 para determinar una de las componentes a partir de las
cuales era posible determinar las 2 restantes.
Por ejemplo en el caso de σ̃1 = 432 esta condición arrojó como resultado nx = ±0.807.
Para determinar las componentes ny y nz usamos las relaciones ny = ± 100 n y nz =
432 x
300
± 432 nx sobre las cuales trasladamos los 2 posibles signos de nx . Lo anterior indica que
las direcciones principales corresponden a vectores con direcciones opuestas lo cual es
equivalente a observar de manera simultanea las tracciones sobre la dirección positiva
y negativa del punto material. Por ejemplo, en el caso de la tensión principal σ̃1 = 432
obtenemos las direcciones n̂1 = 0.807î + 0.187ĵ + 0.560k̂ y su opuesta −n̂1 = −0.807î −
0.187ĵ − 0.560k̂.
Tarea
Comprobar la ortogonalidad de las direcciones principales. Transformar el tensor de

tensiones al sistema de referencia determinado por las direcciones principales y verificar
que el tensor se diagonaliza en dicha dirección.
2.3.1. Circulo de Mohr en 3D
Con los tres valores principales σ̃1 , σ̃2 y σ̃3 es posible realizar una representación plana
del problema de tensión en 3D tal y como se muestra en figura 2.25. Estos circulos
tienen centro y radio dados por;
σ̃1 + σ̃3
σcI =
2
σ̃1 + σ̃2
σcII =
2
σ̃2 + σ̃3
σcIII =
2
σ̃1 − σ̃3
RI =
2
σ̃1 − σ̃2
RII =
2
σ̃2 − σ̃3
RIII =
2
Es claro que para poder graficar el cı́rculo de Mohr en 3D es necesario resolver primero el
problema de tensiones y direcciones principales. Aunque una vez disponible, el cı́rculo
se puede usar para localizar el estado de tensiones en otras direcciones (intersección
de los 3 radios) por el momento solo se hará uso del mismo para calcular el cortante
máximo.
R1
R3 R2
σ~3 C3 σ~2 C2 C1 σ~1 σ
Figura 2.25. Cı́rculo de Mohr
De figura 2.25 es importante resaltar que el esfuerzo cortante máximo está dado por:
σ̃1 − σ̃3
τmax =
2
2.4. Ejercicios
1. La figura 2.26 muestra el estado de tensión en un punto P del Medio Continuo

para varias direcciones: 7
7
Tomado del problema 3.12 en Chávez, E. Problemas resueltos de mecánica del medio Continuo
2
y
6
x
Figura 2.26. Estado de tensiones en el punto en diferentes direcciones.
a) Obtener el estado de tensiones [σ] en el punto P del Medio Continuo.
b) Representar el estado de tensiones mediante el cı́rculo de Mohr.
c) Hallar las tensiones normales máximas y la dirección en la que se presentan.
d ) Utilizando el cı́rculo de Mohr, hallar las tensiones cortantes máximas y la

dirección en la que se presentan.
e) Encontrar el estado de tensiones en un sistema de referencia cartesiano que

se encuentra rotado 300 respecto al sistema mostrado en la figura.
f ) En la figura 2.26 uno de los valores de corte presentado no es correcto. ¿Diga

cual es y por que?
2. Las componentes del tensor de tensiones es un punto P son:

 
1 2 3
[σ] = 2
 4 6 
3 6 1
Encontrar:
a) El vector de tracción t en P para un plano normal al eje x1 ;
b) El vector de tracción t en P para un plano cuya normal es √1 (1, −1, 2);

6
c) El vector de tracción t en P para el plano 2x1 − 2x2 − x3 = 0;
3. La figura 2.27 muestra dos estados de tensiones para el mismo punto de un medio
continuo. 8
10 MPa 10 MPa
10 MPa
θ 20 MPa
θ =30° 50 MPa
θ=60°
a) Determinar los esfuerzos normales sobre los planos mostrados con linea pun-
teada en las figuras para cada uno de los estados.
b) Determinar los esfuerzos cortantes sobre los planos mostrados con linea pun-
teada en las figuras para cada uno de los estados.
c) Para cada uno de los estados de tensiones mostrados determinar los valores
propios y las direcciones principales.
d ) Asuma que ambos estados de tensiones ocurren al mismo tiempo y determi-

ne las direcciones principales y valores propios para el estado de tensiones
resultantes.
e) Para el estado de tensiones resultante de sumar los estados mostrados, de-

terminar el estado de tensiones en un sistema de referencia que se encuentra
a 900 del mostrado en la figura.
f ) Para el estado de tensiones resultante de sumar los estados mostrados, de-

g) Para el estado de tensiones resultante de sumar los estados mostrados, de-

8
Tomado del problema 4.8 en Reddy, J.N (2010). An introduction to continuum mechanics
h) Para el estado de tensiones resultante de sumar los estados mostrados, de-

4. La figura 2.28 muestra un cuerpo sometido a diferentes tensiones en la superficie

del mismo. Escribir el vector de tensiones en cada una de la superficies del medio
continuo. 9
3 kN/m2 2 kN/m
H G D C
F E
5 kN/m2
2.5 kN/m2
A B
Figura 2.28. Medio Continuo sometido a diferentes tensiones en su superficie.
5. En la figura 2.29 se dan las tensiones ~t1 y ~t2 en dos planos no perpendiculares:
t2
t1
60º 60º
de tensiones [σ] en el sistema x − y (horizontal y vertical),

a) Hallar el tensor √
si t1 = 4î + 2 3ĵ y ~
~ t2 = −2î + 0ĵ .
9
6. La figura 2.30 muestra dos estados de tensiones para el mismo punto de un medio
continuo.
Figura 2.30. Estados de tensiones en un punto.
a) Para cada uno de los estados de tensiones mostrados determinar los valores
propios y las direcciones principales.
b) Asuma que ambos estados de tensiones ocurren al mismo tiempo y determi-

ne las direcciones principales y valores propios para el estado de tensiones
resultantes.
c) Para el estado de tensiones resultante de sumar los estados mostrados, de-

d ) Para el estado de tensiones resultante de sumar los estados mostrados, de-

e) Para el estado de tensiones resultante de sumar los estados mostrados, de-

f ) Para el estado de tensiones resultante de sumar los estados mostrados, de-

7. El tensor de tensiones en el punto P está definido como

 
8 −4 1
[σ] =  −4 3 0.5 
1 0.5 2
Calcular el vector de tracción en el punto P según la normal del plano ABC (ver
figura figura 2.31) y descompongalo en su componente normal y tangencial
x3
C(0,0,5)
x2
B(0,2,0)
x1 A(3,0,0)
Figura 2.31
8. Escribir el tensor de tensiones y hacer la representación en el cı́rculo de Mohr

para los siguientes caso:
a) Caso unidimensional, estado de carga de tracción.
b) Caso unidimensional, estado de carga de compresión.
c) Caso bidimensional, estado de carga de tracción.
d ) Caso bidimensional, estado de corte puro.
9. Para un prisma rectangular con sección transversal A, calcular el tensor de esfuer-

zos de Cauchy en un punto en el interior si se somete a una fuerza f distribuida
uniformemente en la parte superior (ver figura 2.32). Para A = 1 m2 y f = 1000
kgf /m2 calcular el valor y la dirección de los cortantes máximos.
x2
x1
Figura 2.32
10. En la figura 2.33 se dan las tensiones para un mismo punto asociados a dos es-
tados I y II. Determinar el estado de tensiones III si el estado I y II ocurren al
mismo tiempo. 10
5
1
1 3
1
+ =
2 45º 45º
σ τ
Estado I Estado II Estado III

11. La figura 2.34a corresponde al montaje de la prueba de laboratorio de compresión

simple usada para determinar la resistencia del concreto.
10
Tomado del ejemplo 4.4 en Olivella, X; Saracı́bar, C (2000). Mecánica de Medios Continuos para
ingenieros
b
c
y
(a) Ensayo de compresión (b) Sección en planta
simple de lado b, ensayo con-
finado
Figura 2.34. Columna en ensayo de compresión
a) Hallar la tensión normal σy mı́nima que produce la falla a cortante en el

modelo presentado en la figura 2.34a, sabiendo que la capaciad máxima del
material ante esfuerzos de corte es τmax = A.
b) Chequear el equilibrio diferencial (local) de la probeta del ensayo.
c) Si ahora se hace el mismo ensayo pero con confinamiento perimetral como

el mostrado en la figura 2.34b, determine el valor minimo del esfuerzo de
confinamiento σc para obtener un un esfuerzo axial máximo σy = 10A y un
cortante máximo τmax = A.
12. En la figura 2.35a se muestra un bloque macizo de arcilla que se extrajo tempo-
ralmente de un muro estructural.Por efecto de funcionamiento del muro, todos los
puntos al interior del bloque se ven sometidos al estado de esfuerzos presentado
en la figura 2.35b.
z
50 5
15 5
5 5
20 5 5 5 20 5
y 10 5
6 5
15 5 y
x
x
(a) Bloque macizo de arcilla ex- (b) Esfuerzos en cualquier punto
traido del muro estructural del bloque (kgf /cm2 ).
Figura 2.35
Para examinar el bloque, este fue cortado como se muestra en la figura 2.36a , y
para colocarlo nuevamente en el muro, las dos piezas en que fue dividido se unirán
con un pegamento, que con el tiempo incrementa su resistencia ante esfuerzos
tangenciales y esfuerzos de tracción, tal como se muestra en la figura 2.36b
7
10
2
4
z
5
6
3
x
(a) Bloque seccionado. (b) Resistencia del pegamento.
Figura 2.36
a) Determine cuál es el tiempo mı́nimo que se debe esperar tras aplicar el

pegamento para poder colocar el bloque nuevamente en el muro sin que
falle.
13. La figura 2.37 muestra el estado de tensiones en un punto al interior de un Medio

Continuo. Sobre la cara asociada al eje y 0 se dan las tensiones σy0 y0 y τy0 x0 y sobre
la cara asociada al eje x0 se da el vector de tracciones expresado en el sistema de
referencia x − y.
Figura 2.37. Estado de tensiones en un punto al interior de un medio continuo.
a) Escribir el tensor de tensiones en el sistema de referencia x0 − y 0 .
b) Escribir el tensor de tensiones en el sistema de referencia x − y.
14. A continuación se presenta el tensor de esfuerzos para un punto de un medio

continuo  de referencia x − y − z.
 en el sistema
1 4 5
[σ] =  4 2 6 
5 6 3
Escriba el tensor de esfuerzos σ 0 asociado al mismo punto del medio conti-

nuo pero en el sistema de referencia cuyos ejes coordenados x0 − y 0 − z 0 están
descritos por los siguientes vectores unitarios:
êx0 = cos(π)î + cos(π/2)ĵ + cos(π/2)k̂
êy0 = cos(π/2)î + cos(π)ĵ + cos(π/2)k̂
êz0 = cos(π/2)î + cos(π/2)ĵ + cos(2π)k̂
15. La figura 2.38 muestra el estado de tensiones en un punto al interior de un Medio

Continuo en el sistema de referencia x − y. Adicionalmente en ésta figura se mues-
tra el sistema de referencia x0 − y 0 que se encuentra rotado y desplazado respecto
al sistema x − y. ¿Cuál de las siguientes respuestas representa correctamente el

tensor en el sistema de referencia x0 − y 0 ?.

a b
a) [σ] =
b −c

−c −b
b) [σ] =
−b a

−c − w −b − (w + z)
c) [σ] =
−b − (w + z) a−z

a−w b − (w + z)
d ) [σ] =
b − (w + z) −c − z

−c b
e) [σ] =
b a
f ) Ninguno de los tensores presentados en las opciones de la (a) a la (g).
c
b
b
a a X
b
z
b
c
Y'
w
X'
16. Se quiere someter una columna a una tensión de compresión uniforme de 100 kgf /cm2
en su superficie superior, ver la punto 16, y se dispone sólo de cinco materiales
diferentes cuyas capacidades máximas se relacionan en el cuadro cuadro 2.1:
Y
Y
yy
yy
X X
yy
(a) Columna sometida a una ten- (b) Tensor de tensiones en cual-

sión de compresión σyy . quier punto al interior de la co-
lumna.
Figura 2.39. Columna sometida a carga vertical en la parte superior
Material 1 Material 2 Material 3 Material 4 Material 5

Capacidad Máxima a 40 60 40 160 50
Cortante
Tracción
Compresión
Cuadro 2.1. Esfuerzos máximos soportados por los materiles en kgf /cm2 .
Determine con cual o cuales materiales se puede construir la columna sin que se
produzca la falla.
17. En las figura 2.40a y figura 2.40b se muestra una viga empotrada en su extremo
izquierdo y sometida a cargas distribuidas, a, y b, respectivamente, sobre las caras
indicadas.
Y Y
a b
H/2 a H/2 b
X X
L L
(a) Viga sometida a carga distri- (b) Viga sometida a carga distri-
buida a. buida b.
Figura 2.40
Se quiere estudiar los esfuerzos en el punto de coordenadas (L, H/2), para lo cual
se pide:
a) Para cada uno de las vigas y sobre los cortes mostrados en la figura 2.41a
indicar la dirección y el valor del esfuerzo normal y tangencial sobre las 6
caras presentadas.
b) Para cada uno de las vigas calcule los vectores de tracciones en el sistema de
referencia XY sobre las 6 caras indicadas, e ilústrelos completamente en la
figura 2.41b
Y Y
X X
(a) Punto (L, H/2). Presentar es- (b) Punto (L, H/2). Presentar
fuerzos normales y tangenciales. vectores de tracciones en x − y.
Figura 2.41
18. Si el tensor de esfuerzos en un punto P , en el sistema de referencia X, Y, Z está

definidido por:
kgf kgf kgf
σxx = 80 ; σyy = 50 ; σzz = 60
cm2 cm2 cm2
kgf kgf kgf

τxy = τyx = 40 ; τxz = τzx = 30 ; τyz = τzy = 10
cm2 cm2 cm2
Y las coordenadas de los puntos A, B y C en el sistema de referencia X, Y, Z

son (1, 0, 0) , (0, 2, 0) y (0, 0, 3) respectivamente. Se pide lo siguiente:
a) Calcular el vector de tracciones ~t, que actúa en la cara del punto P , que es
paralela al plano que contiene los puntos A , B y C.
b) Calcular el esfuerzo normal al que está sometida la cara del punto P , que es
paralela al plano que contiene los puntos A , B y C.
c) Calcular la magnitud del esfuerzo tangencial al que está sometida la cara del
punto P , que es paralela al plano que contiene los puntos A , B y C.
d ) Encontrar la matriz de transformación del sistema de referencia X − Y − Z

al sistema de referencia X 0 − Y 0 − Z 0 para el cual el eje X 0 es paralelo al
vector V~AB y el eje Y 0 es normal al vector normal al plano que contiene los
puntos A, B y C.
19. El estado de esfuerzos en un punto de un medio continuo está dado por:
K(t)
σxx (t) = 0.0 ; σyy (t) = K(t) ; τxy (t) = τyx (t) = −
2.0
Donde, K(t) es un parámetro que varı́a en función del tiempo t, tal como se
muestra en la gráfica y tabla presentadas en la figura 2.42.
200
150
Tiempo K
K (t) [tonf/m²]
t₀ = 0.0 s 0.0 tonf/m² 100

t₁ = 10.0 s 150.0 tonf/m²
50
t₂ = 20.0 s -‐100.0 tonf/m²
t₃ = 30.0 s 0.0 tonf/m² 0
t₄ = 40.0 s 200.0 tonf/m² 0.0 10.0 20.0 30.0 40.0 50.0
-‐50
t₅ = 50.0 s 0.0 tonf/m²
-‐100
Tiempo [s]
Figura 2.42. Variación del parametro K
a) Si el material es infinitamente resistente ante esfuerzos normales (axiales),

kgf
pero no soporta esfuerzos cortantes mayores o iguales a 12 2 determine el
cm
instante en que se presenta la falla en caso de que se presente. Si no hay falla
indı́quelo claramente.
kgf
b) Si el material únicamente soporta esfuerzos tangenciales menores a 15 ,
cm2
kgf
esfuerzos de compresión menores a 3 , y esfuerzos de tracción menores
cm2
kgf
a 20 determine el instante en que se presenta la falla en caso de que se
cm2
presente. Si no hay falla indı́quelo claramente.
20. En la figura 2.43a se presentan tres planos asociados a un punto P de un Medio

Continuo.
(a) (b)
Figura 2.43. Planos de corte en el medio continuo
√
→ 3 2 √
Si se tiene conocen tres vectores de tracciones dados por: t 1 = î + 4 2ĵ +
√ 2
√ → 2 √ √ →
5 2k̂ , t 2 = î + 2 2ĵ − 2k̂, t 3 = 4î + 1ĵ + 2k̂; determine el tensor de
2
tensiones [σ], en el sistema de referencia x − y − z, que proyectado sobre los
→ →
planos con vectores normales n̂1 , n̂2 y n̂3 genera los vectores de tensiones t 1 , t 2
→
y t 3 . Los vectores n̂1 y n̂2 se encuentran contenidos en el plano z − y, tal como
se muestra en la figura 2.43b y n̂3 = 1î.
21. El tensor de tensiones en el punto P está definido como

 
F G −H
[σ] =  G 3F 2G 
−H 2G −F
Calcular el vector de tracción en el punto P según la normal al plano definido por
los puntos A, C y el origen del sistema de referencia (ver figura 2.44). Descomponer
en componentes tangenciales y normal a dicho plano.
C(0,0,4)
A(5,0,0) B(0,6,0) y
x
Figura 2.44
22. Para el medio continuo mostrado en la figura 2.45 se sabe que el tensor de esfuerzos
asociado a un sistema de referencia xy y debido a la acción de una carga P1
(figura 2.45a) está dado por [σ]; y el tensor de esfuerzos asociado a un sistema
de referencia x0 y 0 y debido a la acción de una carga P2 (figura 2.45b) está dado
por un tensor [σ 0 ]. Ambos sistemas de referencia comparten el origen, pero se
encuentran rotados uno respecto al otro.
P1 P2 y'
x Punto A (1/ , 2/ ) x' Punto A (x', y')
(a) Carga P1 y punto en el siste- (b) Carga P2 y punto en el siste-

ma de referencia xy ma de referencia x0 y 0 .
Figura 2.45. Medio Continuo
Se tiene que [σ], en el sistema de referencia xy y debido a la carga P1 , es igual a
x3 yx2
 
−
 (x2 + y 2 )2
2P1  (x2 + y 2 )2 
[σ] =  
π  yx2 xy 2


−
(x2 + y 2 )2 (x2 + y 2 )2
Y [σ 0 ], en el sistema de referencia x0 y 0 y debido a la carga P2 , es igual a

 
0.5 1.5
[σ 0 ] =  
1.5 4.5
Si se quiere estudiar el punto A y se sabe que el valor de la carga P1 es P1 = 25.0,

y las coordenadas en el sistema de referencia xy del punto a estudiar corresponden
1.0 2.0
ax= ,y= , entonces:
π π
a) Encuentre los esfuerzos axiales máximos asociados a [σ] y [σ 0 ] de forma in-

dependiente
b) Determine el ángulo entre el sistema de referencia xy y el sistema de re-

ferencia primado x0 y 0 (ángulo entre el eje x y x0 ). Se sabe que el ángulo
conformado entre la dirección dónde se presenta el esfuerzo máximo a com-
presión del tensor [σ], y la dirección asociada al esfuerzo máximo a tracción
del tensor [σ 0 ] es 0.0◦ .
c) Si se considera que el estado de esfuerzos en un medio continuo está descrito

por la superposición de los efectos de las cargas P1 y P2 , tensores [σ] y [σ 0 ],
determine el esfuerzo tangencial (cortante) máximo al que está sometido el
punto.
23. Un medio continuo es sometido al estado de cargas externas mostrado en la fi-

gura 5.8. Estas cargas producen un estado de esfuerzos en todos los puntos al
interior del cuerpo representado por la figura 2.46b.
P
z
12-P
6
12-P y y
6
x x
(a) Elemento sometido a cargas (b) Tensor de esfuezos
perimetrales
Si el criterio de falla del material está dado por: τ = σ + 2.0 (kgf /cm2 ) y la
carga externa P es aumentada mediante la relación P = Kt + 6.0, con K > 0 y
P es dado en (kgf /cm2 ) . Determine:
a) El valor máximo de la carga P que puede ser aplicada antes que se produzca
la falla en el elemento.
b) El tiempo t en el que se presenta la falla.
c) La dirección del plano en que se produce la falla del elemento.
d ) El valor y la dirección del cortante máximo al que está sometido el punto en

el momento de la falla.
24. Dado el siguiente tensor de esfuerzos asociado a un punto de un medio continuo:

 
5 0 0
[σ] =  0 4 0 
0 0 3
a) Encuentre las direcciones principales del tensor.

b) Encontrar el valor cortante máximo.
c) Encontrar las direcciones en las cuales se presenta el cortante máximo.
d ) Determinar el valor del cortante mı́nimo y direcciones en las cuales se pre-

senta el cortante mı́nimo.
25. El estado de tensión en un punto P del Medio Continuo se da en la figura 2.47. 11
x3
1
4
1
4 σ22
1 1
x2
x1
Figura 2.47. Estado de tensiones en el punto.
Se pide:
a) Determinar el valor de la componente σ22 del tensor de tensiones para que

exista al menos un plano que pase por P que esté libre de tensiones.
b) Determinar el vector normal a dicho plano.
c) Determinar los esfuerzos axiales máximos correspondientes al estado de ten-

siones presentado.
11
Tomado del ejemplo 3.5 en Olivella, X; Saracı́bar, C (2000). Mecánica de Medios Continuos para
ingenieros
d ) Determinar las direcciones principales del tensor de esfuerzo.
e) Determinar la matriz de transformación entre el sistema de referencia x1 −

x2 − x3 y el que es colineal con las direcciones pricipales.
f ) Determinar el máximo valor del esfuerzo cortante.
26. A continuación se presenta el tensor de esfuerzos para un mismo punto de un

medio continuo en los sistemas de referencia x − y − z y x0 − y 0 − z 0 .
   
1/2 1 3/2 A−2 0 0
[σ] =  1 2 3  [σ 0 ] =  0 A 0 
3/2 3 1/2 0 0 5
a) Determine el valor de A.
b) ¿Cuál es el valor del esfuerzo máximo a compresión que se presenta en el

punto?. Determine la dirección con respecto al sistema de referencia x-y-z en
donde dicho esfuerzo se presenta.
c) ¿Cuál es el valor del esfuerzo máximo a tracción que se presenta en el punto?.

Determine la dirección con respecto al sistema de referencia x-y-z en donde
dicho esfuerzo se presenta.
d ) Determine el vector de tracciones respecto al sistema de referencia x-y-z en

la cara en la que se presenta el esfuerzo máximo a compresión.
e) Determine el vector de tracciones respecto al sistema de referencia x-y-z en

la cara en la que se presenta el esfuerzo máximo a tracción.
f ) Escriba el tensor en el sistema de referencia conformado por las direcciones

principales.
g) Cual es la magnitud del máximo esfuerzo cortante que puede experimentar

el punto.
Capı́tulo 3
Equilibrio
3.1. Introducción
En esta sección se formularán las ecuaciones de equilibrio para un medio continuo, ha-
ciendo uso del principio de conservación de moméntum y de momento de moméntum.
Estas darán como resultado un sistema de ecuaciones diferenciales parciales en las com-
ponentes del tensor de tensiones y demostrarán la simetrı́a del tensor de tensiones. Una
vez planteadas las ecuaciones, la sección concluye con el estudio de algunas soluciones1
clásicas para el caso de sólidos.
• Identificar el dominio del problema y su frontera (o superficie).
• Reconocer la diferencia entre equilibrio local y global en un medio continuo.
• Identificar condiciones de frontera de tracciones para diferentes problemas parti-

culares.
• Recuperar y/o verificar condiciones de frontera.

1
Por soluciones nos referimos a funciones de las componentes del tensor de esfuerzos que satisfacen
las ecuaciones de equilibrio y las condiciones de frontera.
83
CAPÍTULO 3. EQUILIBRIO 84
• Entender los aspectos matemáticos fundamentales del problema de valores en la

frontera de la mecánica de los medios continuos en sistemas de referencia carte-
sianos y cilı́ndricos.
• Verificar si una función de tensiones es solución a un problema particular de la

mecánica de los medios continuos.
• Analizar,desde el punto de vista de tensiones, algunas soluciones clásicas de la

mecánica de los medios continuos.
3.2. Ecuaciones de Equilibrio diferencial
En las secciones anteriores se describió el modelo del medio continuo en términos de

fuerzas. El paso del modelo discreto al modelo continuo, vı́a la hipótesis fundamental
de continuidad, se tradujo en un tratamiento de las fuerzas ya no como vectores sino
que estas tuvieron que ser tratadas en términos de funciones tensoriales con significado
de densidades de fuerza o fuerzas por unidad de superficie. Es ası́ como ahora en el
modelo continuo el tensor de tensiones se convierte en el dispositivo matemático que
produce cambios en la cantidad de movimiento en el sistema de infinitas partı́culas.
Nuestro siguiente objetivo es entonces el de proponer las ecuaciones que gobiernan

el flujo de las densidades superficiales de fuerza σ, desde sus valores conocidos en la
frontera en términos del vector de tracciones ~t(n) , hasta cualquier punto del medio
continuo. La determinación de estas tensiones no solamente es un problema con su
propio valor asociado, sino que también es el camino necesario si se desean conocer los
cambios en la cantidad de movimiento del sistema.
Es esperable además, si se sabe que el modelo matemático utilizado para el tratamiento

del sistema de muchas (infinitas) partı́culas ha sido construido a partir de las leyes de la
Mecánica Newtoniana, que las ecuaciones resultantes deban estar estrechamente ligadas
a estas (¿O son las mismas?). Acá partimos entonces de los principios fundamentales
de la mecánica Newtoniana.
Considerese un elemento diferencial arbitrario tomado de un medio continuo y de di-

mensiones dx × dy × dz como el que se muestra en la figura 3.12
2
Debe notarse que este bloque de tensiones corresponde a un elemento del continuo de tamaño
diferencial, mientras que en las secciones anteriores hacı́amos uso de un elemento similar pero que no
Denotemos las áreas superficiales de las caras de dicho elemento y su volumen como;
dSx = dydz ,
dSz = dxdy ,
dSy = dxdz ,
dV = dxdydz .
En lo que sigue nos proponemos plantear las condiciones de equilibrio del elemento
como si se tratara de un cuerpo rı́gido. Para esto supondremos además que el elemento
se encuentra sometido a la acción de fuerzas de cuerpo B~ = Bx ı̂ + By ̂ + Bz k̂.
σxx
τ xy
σ zz + dσ zz τ xz
τ zy + dτ zy
τ zx + dτ zx
τ yz + dτ yz dz
τ yx Bz
τ xz + dτ xz
By σ yy + dσ yy
σ yy Bx y
τ yx + dτ yx
τ xy + dτ xy
τ yz σ xx + dσ xx
dx
τ zx
τ zy
σzz
x dy
Figura 3.1. Tensiones en un elemento diferencial de tamaño dx dy dz.
Dado que el bloque en cuestión tiene tamaño (ası́ sea diferencial) es esperable que las
tensiones sobre una cara y otra ya no permanezcan constantes sino que estas experi-
menten variaciones también diferenciales. Por ejemplo, si la tensión normal sobre la
cara con vector −î y localizada sobre el punto O es σxx entonces se tiene que la tensión
tenı́a ningún tamaño asociado y que solo se convertı́a en una forma práctica de representar el estado
de tensiones en un punto material en un sistema de referencia cartesiano.
Por ejemplo en aquel caso si cortábamos un punto material con un plano con dirección normal n̂ y
con vector de tracciones asociado ~t(n̂) tenı́amos que este vector era igual y opuesto al asociado al corte
con un plano con dirección −n̂. De esta forma se decidió representar el tensor de esfuerzos para un
punto del continuo mediante los tres vectores de tracción asociados con 3 direcciones perpendiculares
ı̂, ̂ y k̂ y sus direcciones opuestas −ı̂, −̂ y −k̂.
∂σxx
normal sobre el punto xO + dx está dada por σxx + ∂x
dx tal y como se muestra en la
figura 3.2.
Con esta consideración de por medio planteamos ahora la condición de equilibrio en la

dirección x: X
Fx = 0 .
Y teniendo cuidado de convertir las tensiones a fuerzas tras considerar el área de cada
cara. Apoyándonos nuevamente en la figura 3.2 en la cual se muestran todas las tensiones
en la dirección x tenemos que

∂σxx ∂τyx
σxx + dx dSx − σxx dSx − τyx dSy + τyx + dy dSy
∂x ∂y

∂τzx
− τzx dSz + τzx + dy dSz + Bx dV = 0
∂z
cancelando los términos iguales y opuestos y eliminando el diferencial de volumen común
a todos los términos se tiene
∂σxx + ∂τyx dydS ∂τzx dydS
dxdS
x y +
z + Bx

dV
= 0.
∂x ∂y ∂z
z
σ xx
τ zx
τ zx + !!!!!!!
! dz
!"
τ yx τ yx
τ yx + !!!!!!! dz
!dy
Bx !"
y
σxx dx
σ xx + !!!!!!! dx
!
!"
τ zx
x dy
Figura 3.2. Variación de las tensiones que producen fuerzas en la dirección x sobre
un elemento diferencial de dimensiones dx, dy, dz.
∂σxx ∂τyx ∂τzx

+ + + Bx = 0. (3.1)
∂x ∂y ∂z
La ecuación (3.1) es una ecuación diferencial parcial que controla la variación de las
componentes σxx , τyx y τzx .
Tarea
P
Usando
P las condiciones de equilibrio en las direcciones y y z (dadas por Fy = 0 y
Fz = 0) para los elementos diferenciales mostrados en las figura 3.3a y figura 3.3b
demostrar que estas resultan en
∂τxy ∂σyy ∂τzy
+ + + By = 0
∂x ∂y ∂z
∂τxz ∂τyz ∂σzz
+ + + Bz = 0 .
∂x ∂y ∂z
z
z
σ zz
σ zz + !!!!!!! τ xz
τ xy !dz
!"
τ zy
τ zy + !!!!!!! dz
!
!" Bz τ
σ yy dy τ yz + !!!!!!!yz! dy dz
σ yy + !!!!!!! !
dz !"
σ yy By !"
τ yz τ xz
τ xz +!!!!!!! y
y !"
!dx
τ xy
τ xy + !!!!!!!
! dx dx
!" dx
τ zy
x dy x dy
σ zz
(a) Fuerzas en y (b) Fuerzas en z
Figura 3.3. Variación de las tensiones que producen fuerzas en la dirección y y z

sobre un elemento diferencial de dimensiones dx × dy × dz.
Por otroPlado, usando ahora las condiciones de equilibrio rotacional se tiene de la con-
dición Mz = 0;

1 ∂τyx 1 1 ∂τxy 1
τyx dSy dy + τyx + dy dSy dy − τxy dSx dx − τxy + dx dSx dx = 0
2 ∂y 2 2 ∂x 2
1 ∂τyx 2 1 ∂τxy 2

τyx dxdzdy + dy dxdz − τxy dxdzdy − dydzdx = 0
2 ∂y
2 ∂x
resultando en la siguiente igualdad entre las tensiones cortantes;
τxy = τyx .
P P
De manera similar planteando Mx = 0 y My = 0 resulta en;
τxz = τzx
τyz = τzy .
Las condiciones de equilibrio traslacional y rotacional se resumen en las ecuación (3.2)

y ecuación (3.3). Las primeras corresponden a 3 ecuaciones diferenciales parciales y las
segundas a la relación de simetrı́a del tensor. Nótese que en total se tienen 6 ecuaciones
en 9 incógnitas3 lo cual significa que desde el punto de vista de tensiones el problema
es estáticamente indeterminado.

+ + + Bx = 0
∂x ∂y ∂z
+ + + By = 0 (3.2)
∂x ∂y ∂z
+ + + Bz = 0
∂x ∂y ∂z
τxy = τyx
τxz = τzx (3.3)
τyz = τzy
En matemáticas un problema definido por un conjunto de ecuaciones que gobiernan

el comportamiento de una función (ecuaciones gobernantes), el dominio de la mismma
3
Las 9 componentes del tensor de tensiones
y las condiciones que la función satisface en la superficie o frontera del dominio, se

conoce como un problema de valores en la frontera (PVF). En nuestro caso particular
el problema esta definido de la siguiente forma:
• La función en cuestión es el tensor de tensiones.
• Las ecuaciones gobernantes son las ecuación (3.2) y ecuación (3.3).
• El dominio será por ejemplo el depósito de suelo, le presa, la viga o la cuña

autosoportada (ver ejemplo).
• Las condiciones de frontera serán el vector de tracciones especificado o conocido

sobre la frontera el cual esta relacionado con el tensor de tensiones mediante la
fórmula de Cauchy tn = σ · n.
Aunque debido a la indeterminación estática no es posible aplicar ningún método formal

de solución al PVF del medio continuo si podemos afirmar que para un problema en
partı́cular la solución del problema debe satisfacer las ecuación (3.2) y ecuación (3.3) y
las condiciones de frontera.
Se dice que una función es solución de un problema de valores en la fron-

tera cuando esta satisface sus ecuaciones gobernantes y las condiciones de
frontera.
Tarea
La ecuación
d du
EA =0
dx dx
gobierna el desplazamiento axial u(x) de una barra de sección transversal variable A(x)
y modulo de elasticidad E. La barra esta sometida a las condiciones de frontera dadas
por

du
EA =F
dx x=L
u|x=0 = 0 .
Si el área de la sección transversal esta dada por
A(x) = A0 (2 − x/L)
demostrar que la siguiente función es solución del problema de valores en la frontera.

FL 2
u= ln
EA0 2 − x/L
3.2.1. Coordenadas cilı́ndricas
Figura 3.4. Tensiones en un elemento diferencial de tamaño drdθdz.
∂σrr 1 ∂τθr ∂τzr σrr − σθθ

+ + + + Br = 0
∂r r ∂θ ∂z r
∂τrθ 1 ∂σθθ ∂τzθ 2τrθ
+ + + + Bθ = 0 (3.4)
∂r r ∂θ ∂z r
∂τrz 1 ∂τθz ∂σzz τrz
+ + + + Bz = 0
∂r r ∂θ ∂z r
Ejemplo 1: Cuña autosoportada
Considere la cuña doble de lado ` y ángulo interno 2φ mostrada en la figura 3.5. Se

asume que esta se encuentra contenida en el plano X − Y , con condiciones de carga
idealizables mediante un estado plano de tensiones (es decir 2D). La cuña se encuentra
cargada por tracciones uniformes de intensidad S aplicadas sobre sus 4 caras de tal
manera que esta se encuentra auto-equilibrada. Determinar la solución al problema de
valores en la frontera proponiendo una solución y verificando las condiciones
de frontera.
l l
S
x
S
S
l l
Figura 3.5. 2D Self-equilibrated wedge.
En el caso de estados de tensiones idealizables en 2 dimensiones y contenidos en el plano

x − y nada es función de z y las ecuaciones de equlibrio ecuación (3.2) y ecuación (3.3)
se reducen a
∂σxx ∂τxy
+ =0
∂x ∂y
(3.5)
∂τxy ∂σyy
+ =0
∂x ∂y
mientras que el equilibrio rotacional se reduce a
τxy = τyx .
Para proponer una solución planteamos las ecuaciones de equilibrio global de la cuña,
lo cual es equivalente a asumir que las solución es constante. Se tiene entonces que
X
Fx = 0 −→ −`S cos φ + σxx ` sen φ = 0
X
Fy = 0 −→ −`S sen φ − σyy ` cos φ = 0
de donde resulta la siguiente función para el campo de tensiones
σxx = S cot φ
σyy = −S tan φ (3.6)
τxy = 0
con la condición τxy = 0 debida a la simetrı́a del problema.
Para verificar que esta solución satisface las ecuaciones gobernantes reemplazamos la
ecuación (3.6) en ecuación (3.5), notando que en este problema en particular no hay
fuerzas de cuerpo (B~ = 0) resultando en;
∂(S cot φ)
=0
∂x
y
∂(−S tan φ)
= 0.
∂y
Sin embargo, existen infinitas funciones σ = σ(x, y) que satisfacen las ecuación (3.5)
y para garantizar que la solución propuesta es la solución correcta al problema es
necesario verificar que esta satisface además las condiciones de frontera del problema.
Sean n̂1 ,n̂2 ,n̂3 ,n̂4 los vectores normales externos a las caras expuestas de la cuña y
escritos en el sistema de referencia x − y como;
n̂1 = − sen φêx + cos φêy

n̂2 = − sen φêx − cos φêy
n̂3 = + sen φêx + cos φêy
n̂4 = + sen φêx − cos φêy .
Para calcular las componentes del vector de tracciones sobre cada una de las caras
usamos la fórmula de Cauchy escrita como
ti = σij n̂j .
Luego, sobre la cara con vector normal n̂1 tenemos
tx = −S cos φ
ty = −S sen φ ,
similarmente sobre la cara con vector normal n̂2
tx = −S cos φ
ty = +S sen φ ,
sobre la cara con vector normal n̂3
tx = +S cos φ
ty = −S sen φ ,
y finalmente, sobre la cara con vector normal n̂4
tx = +S cos φ
ty = +S sen φ ,
los cuales coinciden con las tracciones especificadas como condiciones de frontera del
problema.
Ejemplo 2: Empujes sobre presa de concreto
Si se sabe que la función de tensiones dada en la ecuación (3.7) es solución al problema

de valores en la frontera del medio continuo sobre el dominio mostrado en la figura 3.6,
en donde γ es una constante.4 Se pide:
σxx = γx − 2γy
σyy = −γx (3.7)
τxy = −γy
4
Un script de Python con las funciones de esfuerzo está disponible en: https://github.com/
casierraa/MMC/blob/master/dam.py.
Figura 3.6. Dominio de validez de la solución.
• Identificar la frontera del problema.
• Verificar que la función dada en la ecuación (3.7) efectivamente satisface las con-
diciones de equilibrio local.
• Determinar las condiciones de frontera sobre el dominio (graficar).
• Verificar que la solución satisface condiciones de equilibrio global.
Solución
Sustituyendo la solución dada en ecuación (3.7) en las ecuaciones de equilibrio local y

considerando que Bx = 0 y By = 0 se verifica que la solución satisface equilibrio
∂(γx − 2γy) ∂(−γy)

+ + Bx = 0
∂x ∂y
γ−γ =0
∂(−γy) ∂(−γx)
+ + By = 0
∂x ∂y
0 + 0 = 0.
Sean los vectores normales a las caras vertical, horizontal e inclinada n̂1 = −̂, n̂2 = ı̂
y n̂3 = − sen φ ı̂ + cos φ ̂. En la superficie con normal n̂1 (y = 0) se tiene que;
σxx = γx
σyy = −γx
y
τxy = 0 ,
luego,
tx γx 0 0 0
= = .
ty 0 −γx −1 γx
Esta distribución se muestra en la figura 3.8
!!H
Figura 3.7. Distribución de tracciones sobre la cara vertical.
Similarmente, sobre la superficie con normal n̂2 = ı̂ (x = H) se tiene

σxx = γH − 2γy
σyy = −γH
y
τxy = −γy ,
luego,
tx γH − 2γy −γy 1 γH − 2γy
= = .
ty −γy −γH 0 −γy
La distribución de la componente tx sobre esta cara, corresponde a la tensión normal

σxx tal y como se muestra en la figura 3.7.
!!H
!!H
X
Figura 3.8. Distribución de tensiones normales sobre la cara horizontal.
Similarmente, la componente ty corresponde a la componente τxy . Sobre esta cara se

tiene en y = 0 que ty = 0, mientras que en y = H ty = −γH. Esta distribución se
muestra en la figura 3.9
!!H
Figura 3.9. Distribución de tensiones cortantes sobre la cara horizontal.
En la figura 3.10 se presentan las distribuciones de tensiones sobre toda la presa en

forma de contornos y curvas de nivel de isotensiones (igual valor de tensión).
σxx σyy
0.0
1.00
0.5 0.75
0.50
1.0 0.25
0.00
τxy τmax
0.0 0.25
0.50
0.5 0.75
1.00
1.0
0.0 0.5 1.00.0 0.5 1.0
Figura 3.10. Variación de las tensiones que producen fuerzas en la dirección x y y
sobre la presa.
Para verificar el equilibrio global consideremos el diagrama de cuerpo libre mostrado

en la figura 3.11
1
FY
H/3
H/6 H/6
1
Fx Fx
2
Figura 3.11. Diagrama de cuerpo libre de la presa.
donde
γH 2 γH 2
Fx1 = , Fx2 = − ,
4 4
luego, X
Fx = 0 ,
y
γH 2
Fy1 =
2
luego, X
Fy = 0
similarmente,
X 4 1 5
MA = γH 3 + γH 3 − γH 3 = 0 ,
24 24 24
con lo que se verifica que la solución satisface equilibrio global.
Tarea
Para el problema de la presa;
• Demostrar que la solución dada solo es valida para φ = π/4
• Si la presa es de concreto, el cual presenta baja resistencia a la tracción determinar

la zona de la presa que se debe reforzar para que esta no falle. Nota: Identificar la
superficie donde las tensiones normales (σxx ) cambian de compresión a tracción.
3.3. Soluciones clásicas para sólidos
En esta sección se presentan algunas soluciones clásicas en términos de tensiones. Se

invita al estudiante a realizar el análisis de las mismas usando razonamientos similares a
los presentados en los 2 ejemplos anteriores y de acuerdo con los objetivos de aprendizaje
declarados al inicio del capitulo.
3.3.1. Barra prismática soportada en la superficie superior y

sometida a la acción de su propio peso
La barra prismatica mostrada en la figura 3.25


Figura 3.12. Barra prismática sometida a la acción de peso propio.
se encuentra atada sobre la superficie superior y esta sometida a la acción de su propio

peso. La función dada en la ecuación (3.8) es la solución del problema,
σzz = γz
σxx = σyy = 0 (3.8)
τxy = τxz = τyz = 0
3.3.2. Viga en voladizo
La figura 3.13 muestra una viga en voladizo sometida a una carga P en uno de sus
extremos:
H L
X
H/2
Z
0.5
Y
Figura 3.13. Viga en voladizo.
El tensor solución se da en la punto 1.

P
σxx = − (XY )
I
σyy = σzz = 0
(3.9)
P H2
τxy = τyx = − ( − Y 2)
2I 4
τxz = 0
3.3.3. Cilindro sometido a presión interna y externa
Para este problema se requieren las ecuaciones de equilibrio local en coordenadas polares
dadas por:
∂σrr 1 ∂τθr σrr − σθθ
+ + + Br = 0
∂r r ∂θ r (3.10)
∂τθr 1 ∂σθθ 2τrθ
+ + + Bθ = 0
∂r r ∂θ r
La figura 3.16 muestra un cilindro de radio interno a y radio externo b sometido a

presiones internas y externas pa y pb respectivamente
pb
pa b
x
a
Figura 3.14. Cilindro sometido a presión interna y externa.
El tensor solución se da en la punto 1

2 2

b
r2
−1 1 − ar2
σrr = − b2 pa − a2
pb
2 − 1 1 − 2
2a b
(3.11)
b a2
r2
+1 1 + r2
σθθ = b2 pa − a2
pb
a 2 − 1 1 − b 2
3.3.4. Flamant
La figura 3.15 muestra un semi-espacio sometido a una carga lineal superficial P .

P
y
r
x
Figura 3.15. Medio continuo sometido a una carga P.
El tensor solución se da en la ecuación (3.12)
2P
σrr = − Cosθ
πr
σθθ = 0 (3.12)
τrθ = 0
Ejemplo 3: aplicación de solución para cilindro sometido a pre-

sión interna y externa
1. La solución para un cilindro de radio interno a y radio externo b sometido a

presiones internas y externas pa y pb respectivamente (ver figura 3.16), está dado
por:

b2 a2 b2 a2
r2
−1 1− r2 r2
+1 1+ r2
σrr = − b2
pa − a2
pb σθθ = b2
pa − a2
pb τrθ = 0
a2
−1 1− b2 a2
−1 1− b2
pb
pa b
x
a
Figura 3.16. Cilindro sometido a presión interna pa y presión externa pb
Ejemplo 4: Aplicación de soluciones
En el enlace https://github.com/AppliedMechanics-EAFIT/Notas-MMC/blob/master/notebooks/
PROBLEMA_4_16.ipynb
3.4. Ejercicios
1. Si el tensor de esfuerzos en cualquier punto de la cuña presentada en la figura 3.17 es:

−Scotφ 0
[σ] =
0 Stanφ
S
x
S
l
Figura 3.17. Cuña de espesor e sometida a tensiones tangenciales constantes (S) en

dos caras.
P
a) Verificar el equilibrio global ( F = 0.0).
b) Verificar el equilibrio a nivel diferencial.
c) ¿Es posible encontrar esfuerzos cortantes τxy al interior de la cuña mayores a S?.
Responder sı́ o no y justificar su respuesta.
d ) Calcule el vector de tensiones en cada cara de la cuña.
2. Encontrar las tensiones normales máximas y las direcciones principales de un estado

tensional que está dado por el tensor de tensiones siguiente:
 
1 1 0
[σ] =  1 4 0 
0 0 2
3. En la figura 3.18 se muestra un elemento con densidad ρ, de sección circular, sopor-

tado de su extremo superior y está sometido solamente a la acción de su peso propio.
Se indica el tensor de esfuerzos en el sistema coordenado x − y para los puntos a y b
mostrados en la figura.

0 0 0 0
Punto a: [σ] = Punto b: [σ] =
0 2γL/3 0 γL/3
X Gravedad (g)
= *g
L/3
a
L L/3
2r
b
L/3
2r
Y
Figura 3.18. Barra colgada
Si el tensor de esfuerzos varı́a de forma lineal a lo largo del eje y, se pide:

a) El tensor de esfuerzos para cualquier punto al interior del elemento.
b) ¿El elemento se encuentra en equilibrio a nivel diferencial?. Justifique su respuesta
matemáticamente.
c) Si el elemento es infinitamente resistente ante esfuerzos de tracción y de compre-
sión, pero su capacidad máxima ante esfuerzos cortantes está dada por τmax .
Determine:
1) ¿Cuál es la longitud máxima posible del elemento?.
2) ¿Cuál es la fuerza en el soporte superior cuando se produce la falla por corte?.
3) ¿Cuál es la tensión en la sección transversal inferior del elemento (y = L)
cuando se produce la falla por corte?.
d ) Si el elemento es infinitamente resistente ante esfuerzos cortantes y de compresión,
pero su capacidad máxima ante esfuerzos de tracción está dada por σtrac .
Determine:
cuando se produce la falla por tracción?.
3) ¿Cuál es la fuerza en el soporte superior cuando se produce la falla por trac-
ción?.
e) Si el elemento es infinitamente resistente ante esfuerzos cortantes y de tracción,
pero su capacidad máxima ante esfuerzos de compresión está dada por σcomp .
Determine:
cuando se produce la falla por compresión?.
3) ¿Cuál es la fuerza en el soporte superior cuando se produce la falla por com-

presión?.
4. El tensor de tensiones de un medio continuo se describe con la siguiente diádica:
−2x2 yx 4z 2 + yx 2z 2 yx − xz
 
[σ] =  4z 2 + yx yz + 4y 2 −yz − xyz 

2
2z yx − xz −yz − xyz 4z 3
a) Calcular la fuerza de cuerpo necesaria para mantener el medio continuo en equi-

librio estático.
5. Si en un medio continuo, el campo de tensiones está dado por el tensor:
x2 y (1 − y 2 )x
 
0
[σ] =  (1 − y 2 )x (y 3 − 3y)/3 0 
0 0 2z 2
Determinar:
a) La fuerza de cuerpo necesaria para mantener el medio continuo en equilibrio estáti-
co.
b) Las tensiones principales en el punto P (a, 0, 0).
c) La tensión tangencial (corte) máxima en el punto P .
6. El campo de tensión de un medio continuo está dado por: 5
 
0 0 2x2
[σ] =  0 1 4x1 
2x2 4x1 1
donde x1 y x2 son coordenadas cartesianas.

a) Despreciando las fuerzas de cuerpo ¿está el cuerpo en equilibrio?
b) Determinar el vector tensión que actúa en un punto (1, 2, 3) según el plano x1 +
x2 + x3 = 6;
c) Determinar la proyección del vector tensión según la dirección normal y tangencial
al plano x1 + x2 + x3 = 6.
7. El campo del tensor de tensiones de Cauchy viene representado por sus componentes
como:
x21 x2 (a2 − x22 )x1
 
0
[σ] = k  (a2 − x22 )x1 31 (x32 − 3a2 x2 ) 0  ,
0 0 2ax23
donde k y a son constantes. Encontrar el campo de fuerzas de cuerpo b (por unidad de
volumen) necesario para que el sistema esté en equilibrio.
5
8. En la figura 3.19a se muestra la solución fundamental para semi-espacio sometido a una

carga lineal superficial, P . Los esfuerzos al interior del suelo están dados por:
2P 2P 2P
σxx = − Cos3 θ σyy = − Sen2 θCosθ τxy = − SenθCos2 θ
πr πr πr
P
y
y
r
H
x
x
L1 L2
(a) Solución fundamental carga su- (b) Estructuras y tuberı́a. L1 =

perficial 6m, L2 = 10m, H = 12m
Figura 3.19. Estructuras sobre la superficie
Sı́ se desea instalar una tuberı́a, perpendicular al plano mostrado (XY ), de dı́ametro
despreciable en el punto indicado en la figura 3.19b, y se sabe que el suelo está sometido
a la acción de dos estructuras que le transmiten de forma independiente una carga lineal
superficial, P = 100 T on/m. Determine:
a) ¿Cuál es el esfuerzo mı́nimo de compresión que debe ser capaz de soportar la
tuberı́a para que no se dañe?.
b) ¿Cuál es el esfuerzo mı́nimo de cortante que debe ser capaz de soportar la tuberı́a
para que no se dañe?.
9. El campo de tensión de un medio continuo está representado por:

 
a a a
[σ] =  a a a 
a a a
donde a es una constante.
a) Determinar los valores principales.
b) Determinar las direcciones principales.
c) Determinar el valor del esfuerzo cortante máximo.
10. En la figura 3.20 se muestra un cuerpo cuyo tensor de esfuerzos [σ] en el sistema de
referencia xy está dado por:
q 12 q
σxx = (2x3 y − 4xy 3 + c2 xy) − px σyy = − 3 (4c3 x − 2xy 3 + 6c2 xy)
8c3 5 8c
q 6
τxy = τyx = 3 [3x2 (c2 − y 2 ) − (c4 − y 4 ) + c2 (c2 − y 2 )]
8c 5
Donde q, p y c son constantes.
c c
y
x
Figura 3.20. Medio continuo
En las siguientes preguntas seleccione la opción correcta y justifique matemáticamente

su selección:
a) El cuerpo se encuentra en equilibrio a nivel diferencial si:
1) La fuerza de cuerpo en la dirección y es p
2) No hay fuerzas de cuerpo
3) La fuerza de cuerpo en la dirección x es q
4) La fuerza de cuerpo en la dirección x es px.
5) La fuerza de cuerpo en la dirección x es p.
b) El vector de tracciones sobre la cara y = c en el sistema de referencia xy está dado
por:
" q 3 c − 4xc3 + 12 c3 x) − px
#
(2x
1) [t] = 8c3 5
−qx

−qx
2) [t] =
0

0
3) [t] =
−qx
" q 3 c − 4xc3 + 12 c3 x) − px
#
− (2x
4) [t] = 8c3 5
qx
5) Ninguna de las anteriores.
c) La magnitud de la fuerza en direccion y que actua en la cara y = c es:
1) F = qH 2 .
qH 2
2) F = .
2
3) F = 2qc.
qH
4) F = .
2
5) Ninguna de las anteriores.

d ) Si el tensor [σ] se emplea como la solución de esfuerzos para la presa mostrada en
la figura 3.21, seleccione cual de las siguientes afirmaciones es correcta.
c c
y
Presa
almacenado H
Suelo de soporte x
Figura 3.21. Presa
1) El tensor satisface plenamente las condiciones de frontera en todas las caras

(y = c, y = -c, x = 0, x = H).
2) El tensor satisface plenamente las condiciones de frontera en las caras y = c

y x = 0.
3) El tensor no satisface las condiciones de frontera en la cara y = -c.
4) El tensor no satisface las condiciones de frontera en la cara x = 0.
5) Ninguna de las afirmaciones es correcta.
e) Aceptando como solución el tensor presentado en el numeral [1] y considerando
que la presa está construida con suelo compactado cuya resistencia ante esfuerzo
kgf kgf kgf
cortante es de 0.5 2 ; si H = 3.0 m, c = 1.1 m, q = 1000 3
y p = 2000
cm m m3
se puede afirmar:
1) La presa falları́a ante esfuerzos de corte en el plano xy porque el cortante en
kgf
el punto de coordenadas (3.0 m,1.1 m) es mayor a 0.5
cm3
2) La presa falları́a ante esfuerzos de corte en el plano xy porque la magnitud
del cortante máximo en el punto de coordenadas (0.0 m,-1.1 m) es mayor a
kgf
0.5 .
cm2
3) La presa no falları́a ante esfuerzos de corte en el plano xy porque la magnitud
del cortante máximo en el punto de coordenadas (3.0 m,0.0 m) es menor a
kgf
0.5 .
cm2
4) La presa falları́a ante esfuerzos de corte en el plano xy porque la magnitud
del cortante máximo en el punto de coordenadas (3.0 m,-1.1 m) es mayor a
kgf
0.5
cm2
11. En la figura 3.22 se muestra un cuerpo de sección rectangular de ancho unitario, altura
H, soportado en un extremo (X = 0), y sometido en su otro extremo (X = L) a la
acción de una fuerza P en la dirección Y. Si el tensor de esfuerzos en el sistema coor-
denado XY Z está dado por:
P
σxx = − (XY ) σyy = σzz = 0
I
P H2
τxy = τyx = − ( − Y 2) τxz = 0
2I 4
Donde I es el momento de inecia de la sección respecto a el eje Z y P es la fuerza.
H L
X
H/2
Z
0.5
Y
Figura 3.22. Viga en volodizo
a) Verifique equilibrio diferencial.

b) Verifique las condiciones de frontera.
c) ¿Las condiciones de frontera coinciden con las presentadas en el enunciado del
problema?. De ser afirmativa o negativa su respuesta, justifı́quela.
d ) Verifique equilibrio global.
e) ¿Cual es el valor del esfuerzo cortante máximo y en donde se presenta?.
f ) ¿Cual es el valor del esfuerzo normal máximo y en donde se presenta?.
12. En la figura 3.23a se muestra una estructura especial actuando sobre un cuerpo de
concreto de gran dimensión. El estado de esfuerzos al interior del cuerpo debidos a la
acción de la estructura, en un sistema coordenado cilı́ndrico (r, θ, z), están dados por:
2Q Q√
10mmσrr = − Cosθ; σθθ = σzz = 0; τrθ = τrz = 0; τθz = 3
πr πr
Q es una constante positiva, r es la distancia desde la ubicación de la estructura al
punto de evaluación y θ es el ángulo medido desde la horizontal en sentido antihorario.
2h
x x
h
r
y y
h
(a) Solución fundamental (b) Estructuras y cable
Si se desea instalar un cable (de √

dı́ametro despreciable) perpendicular al plano mostrado
(XY ) a una profundidad h = 4 2 y se sabe que el cuerpo está sometido a la la acción
de dos estructuras iguales, tal como se muestra en la figura 3.23b. Determine:
a) Si el cable es infinitamente resistente ante esfuerzos normales pero τcrit es el

máximo esfuerzo de corte que resiste, ¿cuál es el máximo valor de la carga Q que
pueden transmitir las estructuras para que el cable no falle?.
b) Si τcrit = 10, ¿cuál es el valor de Q?.
c) Si Q = 5, ¿cuál es el valor mı́nimo de τcrit para que el cable no falle por corte?.
13. En las figura 3.23c y figura 3.23d se muestra las soluciones para una masa de suelo
sometido a una carga superficial P y una carga distribuida uniforme superficial w .
Los esfuerzos al interior del suelo debidos a la carga P y w, en un sistema coordenado
cilı́ndrico (r, θ, z), están dados:
   
cos θ 0 0 π + 2θ − sen (2θ) 1 − cos (2θ) 0
2P  w
[σ]P = − 0 0 0  [σ]w = − 1 − cos (2θ) π + 2θ + sen (2θ) 0 
πr π
0 0 0 0 0 0
Donde P y w son las cargas, θ es el ángulo medido desde la vertical y positivo en sentido
antihorario.
P W
y y
r r
x x
(c) Solución carga vertical P (d) carga uniforme superficial w
Figura 3.23. Soluciones fundamentales
Si a la masa de suelo se transmiten de manera simúltanea las cargas P = 100 T on/m y

w = (100 + tiempo) T on/m2 , y se quiere instalar un sistema de acueducto conformado
por dos tuberı́as A y B (de dı́ametro despreciable) perpendicular al plano mostrado
(XY ), en donde la profundad para la tuberı́a A es HA = 1.0/π m y para la tuberı́a B
es HB = 0.75HA tal y como se muestra en la figura 3.24. Determine:
P HA W
y
HB
x B HA
A
Figura 3.24. Cargas y tuberı́as
a) Determine cuál tuberı́a falla primero y en que tiempo se produce la falla, si se sabe
que el material de la tuberı́a A es infinitamente resistente a esfuerzos cortantes
T on
pero no soporta esfuerzos axiales mayores o iguales a 210, 5 2 y el material de la
m
tuberı́a B es infinitamente resistente a esfuerzos axiales pero no sorporta esfuerzos
T on
cortantes mayores o iguales a 100, 25 2 .
m
14. En la figura 3.25 se presenta una barra de radio C y longitud L. El estado de esfuerzos
del medio continuo en el sistema cartesiano x, y, z está dado por el tensor:
σxx = 0; σyy = 0; σzz = 0;

(3.13)
τxy = 0; τxz = ωy; τzy = −ωx
C x
y L
z
Donde ω es una constante.
a) Determine el valor del máximo esfuerzo de compresión al que está sometido la

barra y la ubicación (coodenadas de los puntos) donde se presentan esos valores
máximos.
b) ¿Describa que problema resuelve el estado de esfuerzos en a barra? Sugerencia:

para facilitar el análisis estudie las caras con vector normal en dirección z.
c) Si el cuerpo tiene longitud L = 50cm, ω = 5000 T onf /m3 , además es infinitamente
resistente ante esfuerzos de tracción y de compresión, pero su capacidad máxima
ante esfuerzos tangenciales es τmax = 88 kgf /cm2 . Determine el valor máximo del
radio C que podrı́a tener el cuerpo.
15. Si a la barra del punto punto 14 adicionamente le quisieramos considerar el efecto del
peso propio γ y sabemos que el tensor de esfuerzos en el sistema coordenado x, y, z en
cualquier punto de los planos z = L/3 y z = 2L/3 (ver figura 3.26) debido solo a este
efecto está dado por:
   
0 0 0 0 0 0
P untos(z=L/3) : [σ] =  0 0 0  P untos(z=2L/3) : [σ] =  0 0 0 
0 0 2γL/3 0 0 γL/3
γ = ρg, ρ es la densidad y g la gravedad.
C x
L/3
y L/3 L
z
L/3
Figura 3.26. Función peso propio
Además dicho esfuerzo es descrito por una función lineal en z.
a) Si se tiene en cuenta que estado final de esfuerzos es la suma de estados indivi-

duales, determine cuál es el tensor de esfuerzos del problema que considera con-
juntamente el tensor de esfuerzos del punto 15 y el peso propio.
b) Diga si el cuerpo está o no en equilibrio diferencial (local) estático. Sı́ es del caso,
de una propuesta para lograr dicho equilibrio. Justifique su respuesta matemáti-
camente
16. El estado de esfuerzos en un punto de un medio continuo está dado por:

σxx (t) = K(t); σyy (t) = −6K(t); σzz (t) = −K(t);
τxy (t) = 0; τxz (t) = 0; τyz (t) = 0
Donde, K(t) es un parámetro que varı́a en función del tiempo t, tal como se mues-
tra en la gráfica y tabla presentadas en la figura 3.27.
4
t (s) K(t) (kgf/cm2 )
2
K(t) (kgf/cm²)
0.0 2.0 0
5.0 2.0 2
10.0 4.0 4
15.0 0.0 6
20.0 -6.0 0 5 10 15 20 25
25.0 0.0 t (s)
Figura 3.27. Variación del parámetro K

a) Si el material es infinitamente resistente ante esfuerzos normales (axiales), pero
kgf
no soporta esfuerzos cortantes mayores o iguales a 13 2 determine el instante
cm
en que se presenta la falla en caso de que se presente. Si no hay falla indı́quelo
claramente.
17. En la figura 3.28 se presenta la sección transversal de un cilindro. Los esfuerzos en el

dominio están representados por:
a2 b2 (po − pi ) pi a2 − po b2 a2 b2 (po − pi ) pi a2 − po b2
σr = + σθ = − + τrθ = 0
(b2 − a2 )r2 b2 − a2 (b2 − a2 )r2 b2 − a2
b
x
a
Figura 3.28. Sección transversal cilindro
a) Verifique el equilibrio diferencial.

b) ¿Cuales son las cargas externas a las cuales está sometido el cilindro?. Indı́quelas
gráficamente.
c) ¿Cual es el esfuerzo axial máximo y dónde se presenta?.
d ) ¿Cual es el esfuezo de corte máximo y dónde se presenta?.
18. En la√figura√3.29 muestra el estado de tensiones en un punto de coordenadas (X, Y ) =

(55/ 2, 55/ 2) al interior de un Medio Continuo. Sobre una de las caras asociadas al
eje Y 0 se dan las tensiones σY 0 Y 0 y τY 0 X 0 y sobre una de las caras asociadas al eje X 0 se
da el vector de tracciones expresado en el sistema de referencia XY .
Y
Y
X'
Y'
(r, )
X
X
a) Escribir el tensor de tensiones en el sistema de referencia polar r, θ.

19. Los tensores

   
2P cos θ 2Q sen θ
[σ(r, θ)]P = − (2α + sen(2α))r
0 
[σ(r, θ)]Q = − (2α − sen(2α))r
0 
0 0 0 0
representan el estado de tensiones en un punto de coordenadas (r, θ) al interior de las

cuñas de semi-ángulo interno α (medido en radianes) mostradas en las figura 3.30a y
figura 3.30b, resultantes de la aplicación de cargas puntuales P y Q respectivamente.
Y Q Y
α α α α
r r
θ θ
X X
(a) Cuña de semi-ángulo interno (b) Cuña de semi-ángulo interno

α sometida a una carga vertical α sometida a una carga vertical
P. Q.
Figura 3.30. Cuña de semi-ángulo interno α.
Usando dichas soluciones se pide determinar el estado de tensiones para un punto P de

coordenadas (r, θ) localizado al interior de un depósito de suelo sometido a una carga
puntual de magnitud F como se muestra en la punto 4.
F
β
Y
r
θ
P(r,θ )
Figura 3.31. Depósito de suelo sometido a una carga puntual F .

Capı́tulo 4
Análisis de Deformaciones
4.1. Introducción
• Reconocer la diferencia entre desplazamientos absolutos y relativos entre los pun-

tos materiales.
• Entender los cambios geométricos impartidos al aplicar una transformación lineal

sobre un campo vectorial.
• Identificar los efectos fı́sicos contenidos en el tensor de deformaciones infinitesi-

males.
• Calcular el tensor de deformaciones infinitesimales a partir del campo de despla-

zamientos relativos.
• Reconocer la diferencia entre la configuración deformada de un medio continuo y

la deformada del punto material
Del análisis de las fuerzas “internas” se concluyó —usando las 3 leyes de Newton sobre
un volumen arbitrario de un medio continuo— que estas corresponden a nivel local a
una función (continua) de carácter tensorial y que satisface el problema de valores en la
frontera especificado por las ecuaciones de equilibrio (leyes generalizadas de Newton) y
118
CAPÍTULO 4. ANÁLISIS DE DEFORMACIONES 119
la asignación de las tracciones en la superficie del medio (ecuación (4.1), ecuación (4.2)):

+ + + Bx = 0
∂x ∂y ∂z
+ + + By = 0 (4.1)
∂x ∂y ∂z
+ + + Bz = 0 ,
∂x ∂y ∂z
y
τxy = τyx
τxz = τzx (4.2)
τyz = τzy .
Se dispone de 6 ecuaciones en 9 incógnitas controlando la distribución espacial de las

fuerzas internas σ sobre los diferentes puntos materiales del medio continuo. Dichas
fuerzas internas son generadas por fuerzas de superficie (tracciones) ~t(n̂) o por interac-
~ Sin embargo, en este punto el problema se hace
ciones a distancia o fuerzas de cuerpo B.
(matemáticamente) indeterminado ya que se dispone de más incógnitas que ecuaciones.
La indeterminación desde el punto de vista de fuerzas no debe ser para nada sorpresiva,
ya que si éstas se consideran como las “funciones causa” aún falta involucrar los efectos
o “funciones respuesta”. Estas funciones corresponden precisamente a los cambios (ci-
nemáticos) o de configuración, experimentados por las infinitas partı́culas una vez son
sometidas a interacciones externas.
Claramente se identifica la primera necesidad de estudiar los cambios de configuración

como la de tratar de establecer una conexión entre estos cambios y las fuerzas internas
σ y que permita a la postre eliminar la indeterminación matemática del problema de
valores en la frontera (PVF) (especificado por las ecuación (4.1) y ecuación (4.2)).
Aunque es claro que las nuevas configuraciones aparecen debidas a los desplazamientos
a los que son sometidas las diferentes partı́culas, a la hora de establecer la conexión
entre éstas y las fuerzas internas σ es importante considerar que son en realidad los
desplazamientos relativos quienes generan estas fuerzas internas. La aparición de des-
plazamientos relativos necesariamente se reflejará en términos de cambios de forma del
medio continuo.
4.1.1. Concepto intuitivo de deformación
La figura 4.1 muestra, en sombreado gris, una viga en voladizo en su configuración

original (antes de aplicar las cargas externas) y en lı́nea punteada su configuración final
(o deformada)1 resultante tras aplicar alguna distribución de cargas externas en forma
de tracciones.
Como resultado de la aplicación de las cargas externas los puntos materiales que con-
forman la viga experimentan desplazamientos. En lo que sigue denominaremos campo
de desplazamientos a la variación espacial de estos, y lo denotaremos como ~u(~x). Tal
y como se aprecia de manera intuitiva, los desplazamientos varı́an para los diferentes
puntos de la viga generándose desplazamientos relativos.
Figura 4.1. Configuración original y deformada de una viga en voladizo.
Por ejemplo, es evidente que los desplazamientos sobre el empotramiento son ”nu-
los”mientras que los puntos cercanos al extremo derecho de la viga experimentan los
máximos desplazamientos. Sin embargo, sabemos que las fuerzas internas (o tensiones)
están ligadas al desplazamiento relativo entre los puntos materiales. El propósito de
esta sección es hacer una descripción local (o diferencial) de los desplazamientos rela-
tivos entre puntos materiales y que pueda ser conectada posteriormente con el tensor
de tensiones a través de propiedades de los materiales . Con esto, no solo se conseguirá
completar las ecuaciones para el problema de valores en la frontera del medio continuo,
sino que además será posible describir la nueva configuración del medio, resultante de
la aplicación de interacciones externas tanto a nivel local (diferencial) como global.
Para aclarar, nuevamente de manera intuitiva, el carácter local de los cambios relativos
la figura 4.2 muestra la configuración original (figura 4.2a) y deformada (figura 4.2b)
1
En la configuración deformada los puntos materiales o elementos constitutivos de la viga estarán
sometidos a un estado de tensiones que satisface las ecuaciones de equilibrio local ( a nivel diferencial)
y global.
de la viga en voladizo. Como elemento auxiliar de análisis esta vez se ha dibujado una
rejilla de elementos perfectamente cuadrados sobre la configuración no-deformada. Si
suponemos que estos elementos son de lado h es claro que en el modelo del continuo
el punto material corresponde a uno de estos en el estado limite cuando h → 0. Es
evidente de la figura 4.2b como los diferentes elementos cambian su forma una vez la
viga pasa de la configuración descargada a la cargada. Mas aún es evidente de la misma
figura como estos cambios varı́an de elemento a elemento.
(a) Configuración original con ele-

mentos diferenciales no deformados
(b) Configuración final con elementos

diferenciales deformados
Figura 4.2. Comparción entre configuración original y deformada con elementos

diferenciales.
La figura 4.3 muestra nuevamente ambas configuraciones y sobre estas se han seleccio-
nado para efectos ilustrativos dos elementos localizados sobre el empotramiento y sobre
el extremo libre del voladizo. En la configuración original ambos elementos son idénticos
y corresponden a cuadrados de lado h. En la configuración deformada, mostrada con
lı́neas punteadas, se resaltan nuevamente estos 2 elementos de donde se pueden hacer
las siguientes observaciones:
El elemento sobre la parte superior del empotramiento experimenta cambios de tamaño

y de forma a través de alargamientos en la dirección longitudinal, mientras que el
elemento sobre la parte inferior experimenta acortamientos igualmente en la dirección
longitudinal. Estos cambios suceden a pesar de que dichos elementos se encuentran
sobre el empotramiento donde los desplazamientos como tal son nulos. De otro lado, el
elemento sobre el extremo libre de la viga experimenta pequeños cambios de forma (no
apreciables en el dibujo) y experimentan las mayores rotaciones como si fuera un cuerpo
rı́gido. Los pequeños cambios de forma son equivalentes a pequeños desplazamientos
relativos aunque estos puntos experimentan los mayores desplazamientos absolutos.
Figura 4.3. Elementos con máxima deformación y máxima rotación.
En resumen, el propósito de esta sección es hacer una descripción local (o diferencial) de

los desplazamientos relativos entre los puntos materiales de un medio continuo cuando
este es sometido a la acción de fuerzas externas. En general, como se verá mas adelan-
te, estos cambios se describirán en términos de cambios de tamaño (o magnitud) y de
orientación de las ”fibras materiales”que emanan de la partı́cula en las diferentes direc-
ciones tal y como se esquematiza en la figura 4.4. Esta muestra un punto material de
un medio continuo conjuntamente con su vecindad matemática representada por fibras
imaginarias de material, mostradas en lı́neas punteadas. En la figura el punto material
en estudio se muestra con relleno de color negro y se denomina P . Adicionalmente se
muestra un punto material arbitrario (sin relleno) y denominado P 0 .En el lado izquier-
do se muestra la localización de ambos puntos materiales en un instante de tiempo
t = 0, correspondiente a la configuración original o descargada. Tras la aplicación de
las cargas ambos puntos materiales se desplazan a las posiciones señaladas como Q y Q0
respectivamente y mostradas en lado derecho de la figura correspondiente a un instante
de tiempo t = t. Como resultado de los desplazamientos experimentados por los puntos
materiales la fibra material que une los mismos experimenta cambios de tamaño y de
orientación.
t =0 t =t
P'
Q Q'
P
U ( A, t )
A
x
Figura 4.4. Fibras materiales en la configuración original y deformada.
En lo que sigue estudiaremos la relación entre 2 fibras arbitrarias como las mostradas
en la figura desde 2 puntos de vista. Inicialmente se estudiarán los cambios de magnitud
y orientación experimentados por un vector posición tras operarlo con una matriz con-
servando la propiedad de lı́nealidad. Posteriormente estudiaremos el problema a nivel
infinitesimal tras analizar un par de puntos, unidos por una fibra material de tamaño di-
ferenial. Dicho analisis permitirá encontrar el tensor graiente de desplazamientos como
medida local (o diferencial) de la deformación2 .
4.2. Transformaciones lineales
Con el propósito de conceptualizar el proceso deformación o cambios de configuración

a nivel local o diferencial es conveniente apoyarnos como herramienta de análisis en las
transformaciones lineales3 . En términos simples estas corresponden a la traformación
de un vector en otro conservando la propiedad de lı́nealidad. Esta idea se ilustra en
la ecuación (4.3) en la que ~r = rx î + ry ĵ representa el vector posición de un punto
P . Como resultado de alguna acción externa, representada por un tensor A este punto
material experimenta un desplazamiento ∆ ~ = ∆x î + ∆y ĵ con lo que que el vector ~r se
2
Aunque hasta el momento el tratamiento del medio continuo no ha especificado ninguna diferencia
entre medios sólidos, liquidos o gases el tratamiento cinematico que se presenta si esta restringido al
caso de pequeñas deformaciones
3
Esta sección esta basada en los textos: Principles of Solid Mechanics. Rowland Richards Jr. CRC
Press, 2001 y Introduction to the Mechanics of Continuous Medium. Lawrence E Malvern. Prentice
Hall, 1969.
transforma en el vector ρ~ = ρx î + ρy ĵ. Esta relación puede escribirse como:

ρ~ = A · ~r (4.3)
z
Q(x,y,z)
Δ
P(x,y,z)
ρ
r
o y
Figura 4.5. Transformación del vector ~r en el vector ρ

~.
El vector ~r al transformarse en el vector ρ~ experimenta 2 tipos de cambios:
i Cambios de magnitud (o de tamaño).
ii Cambios de orientación (dirección).
La relación ecuación (4.3) entre ~r y ρ~ puede escribirse en forma explicita como:

ρx axx axy rx
= (4.4)
ρy ayx ayy ry
La ecuación (4.4) indica que la matriz A transforma el vector ~r en el vector ρ~ y por ende
ésta contiene toda la información referente a los cambios de magnitud y dirección que
experimentan las infinitas direcciones que emanan del punto O. Una transformación
como la descrita en la ecuación (4.4) se denomina lineal si satisface la condición dada
por:
A · (~r1 + ~r2 ) = A · ~r1 + A · ~r2 (4.5)
la cual implica que lineas rectas y paralelas antes de la transformación permanecen

rectas y paralelas después de la transformación o en otras palabras que conservan la ley
del paralelogramo. Esto se ilustra en la figura 4.6:
z c'
a'
r2 )
[A] r1
](
r1
[A
[A]r2 d'
o
r1 y
r2 a
r1 r
b
2
x
c
Figura 4.6. Definición de linealidad.
Alternativamente, es posible escribir la transformación en términos del vector de des-

~ que experimenta la cabeza del vector ~r como sigue:
plazamientos ∆
∆~ = ρ~ − ~r
∆~ = A · ~r − ~r
ρ~ = [A − I] · ~r
∆~ = D · ~r
y de manera explı́cita escribimos

∆x dxx dxy rx
= (4.6)
∆y dyx dyy ry
La matriz D se denomina la matriz de transformación de desplazamientos. En lo que

sigue estudiaremos transformaciones como las dadas en la ecuación (4.6) identificando
los efectos que dicha matriz impone sobre las diferentes direcciones que emanan de un
punto material.
4.2.1. Componente rotacional
A continuación nos disponemos a estudiar el efecto de los diferentes términos en la

matriz de transformación de desplazamientos D. Abordemos inicialmente los cambios
de orientación impartidos por D en los diferentes vectores. Considerando la figura 4.7,

en esta el vector original ~r se transforma en el vector ρ tras realizar una rotación α.
r α
ρ
x
Figura 4.7. Rotación a través de un ángulo α del vector ~r en el vector ρ
~.
Relacionando las componentes rectangulares de ambos vectores se tiene que:
ρx = rx + αry
y
ρy = ry − αrx .
De manera explicita, pero en términos de los desplazamientos de la cabeza del vector

posición ~r lo anterior resulta en:

∆x 0 α rx
= (4.7)
∆y −α 0 ry
De la ecuación (4.7) se concluye que la componente rotacional del tensor de transfor-

mación de desplazamientos D es una matriz con ceros en la diagonal y términos iguales
pero con signo opuesto a uno y otro lado de la diagonal. En el lenguaje del análisis
tensorial, extendido en este caso a matrices, a estos tensores se les conoce como tenso-
res anti-simétricos. En conclusión la componente rotacional de la transformación lineal
corresponde a un tensor anti-simétrico.
4.2.2. Componente simétrica
Como ya ha sido identificado en los planteamientos anteriores, la matriz de transforma-

ción de desplazamientos D imparte en las diferentes direcciones cambios de dirección y
de magnitud. El haber identificado el efecto rotacional como la componente anti-simétri-
ca de D permite concluir que el término restante, y simétrico por demás, contendrá los
efectos asociados a los cambios de tamaño o magnitud de los diferentes vectores. De
dicha conclusión resulta natural resolver el tensor D precisamente en una componente
simétrica y en una componente anti-simétrica como:
dxy +dyx dxy −dyx
dxx dxy dxx 2
0 2
= dyx +dxy + dyx −dxy . (4.8)
dyx dyy 2
d yy 2
0
En esta partición es evidente que el segundo término del lado derecho corresponde efec-
tivamente a la componente rotacional del tensor (componente anti-simétrica), mientras
que el primer término necesariamente contendrá la información relativa a los cambios
de tamaño. Nótese además que este primer término es simétrico con respecto a la dia-
gonal. De forma general escribimos esta partición de la matriz de transformación de
desplazamientos D como:
D =ε+ω (4.9)
donde ε contiene los efectos de cambio de tamaño y ω contiene los efectos de cambio
de orientación.
Para estudiar los efectos de las diferentes componentes del tensor de transformación de
desplazamientos D usaremos el cuadrado unitario de la figura 4.8.
B(0,1) C(0,1)
A(0,1) x
Figura 4.8. Cuadrado unitario a ser transformado por el tensor D.
Como método de análisis aplicaremos las diferentes componentes del tensor D sobre los
vectores posición de los puntos A, B y C dados por:

~rA = 1.0î
~rB = 1.0ĵ
~rC = 1.0î + 1.0ĵ
y asumiremos que los puntos O, A, B y C permanecen unidos por lı́neas rectas de
acuerdo con la propiedad de linealidad. Posteriormente compararemos la configuración
original y deformada del elemento.
Consideremos inicialmente la componente simétrica de la ecuación (4.9). Por razones

que se harán evidentes más adelante denotemos los elementos de esta componente como:

εxx γ̄
ε=
γ̄ εyy
Para determinar la configuración deformada aplicamos la expresión:

∆x εxx γ̄ rx
=
∆y γ̄ εyy ry
a los vectores posición ~rA , ~rB y ~rC . Esto da como resultado los vectores de desplaza-
miento:
∆~ A = εxx ı̂ + γ̄̂
~ B = γ̄ı̂ + εyy ̂
∆
~ C = (εxx + γ̄)ı̂ + (εyy + γ̄)̂
∆
y la configuración deformada mostrada en lı́nea punteada en la figura 4.9.
y
C'
B'
B C
A'
A x
Figura 4.9. Configuraciones original y deformada de un cuadrado unitario debido a

la componente simétrica ε del tensor de transformación de desplazamientos D.
4.2.3. Componentes escalar y distorsional
La componente simétrica del tensor de transformación de desplazamientos es análoga

al tensor de tensiones estudiado en secciones anteriores. De hecho, dado su carácter
simétrico este puede ser estudiado en términos del circulo de Mohr. Identificando esta
conexión es posible resolver esta componente en términos adicionales que permiten
ahondar en el entendimiento fı́sico del análisis de deformaciones. En particular considere
la siguiente descomposición:
εxx +εyy εxx −εyy
εxx γ̄ 2
0 2
0 0 γ̄
= εxx +εyy + εyy −εxx + (4.10)
γ̄ εyy 0 2
0 2
γ̄ 0
en la cual se ha separado la componente simétrica en un término escalar (correspon-

diente al centro del cı́rculo de Mohr) y en 2 términos con información distorsional
controlando el radio del cı́rculo. Escribamos la descomposición de ε como:
ε = pI + S + C
εxx +εyy
donde se reconoce p = 2
con lo que el primer término se escribe como:
εxx +εyy
2
0
pI = εxx +εyy .
0 2
εxx −εyy
Similarmente, haciendo s = 2
escribimos:
εxx −εyy
0 s 0
S= 2
εyy −εxx ≡
0 2
0 −s
y finalmente
0 γ̄
C= .
γ̄ 0
Para entender el efecto fı́sico de cada uno de estos términos aplicaremos nuevamente la
transformación al cuadrado unitario previamente definido.
La componente escalar pI genera los siguientes vectores de desplazamientos tras aplicar

la ecuación (4.6):
~ A = pı̂
∆
~ B = p̂
∆
~ C = pı̂ + p̂
∆
resultando en el cambio de forma mostrado en la figura 4.21
B' C'
B CC
A x
A'
Figura 4.10. Configuraciones original y deformada de un cuadrado unitario debido

a la componente escalar de la componente simétrica del tensor de transformación de
desplazamientos.
De manera similar, la primera componente distorsional S genera los siguientes vectores

de desplazamientos:
~ A = sı̂
∆
~ B = −s̂
∆
~ C = sı̂ − s̂
∆
resultando en el cambio de forma mostrado en la figura 4.11.

B'
B C
C'
A A' x
Figura 4.11. Configuraciones original y deformada de un cuadrado unitario debi-

do a la primera componente distorsional de la componente simétrica del tensor de
transformación de desplazamientos.
Luego, aplicando la segunda componente distrosional C da como resultado los siguientes

vectores de desplazamientos:
~ A = γ̂
∆
~ B = γı̂
∆
~ C = γı̂ + γ̂
∆
y en los cambios de configuración mostrados en la figura 4.12
C'
B'
B
C
A'
A x
Figura 4.12. Configuraciones original y deformada de un cuadrado unitario debi-

do a la segunda componente distorsional de la componente simétrica del tensor de
transformación de desplazamientos.
Finalmente, y en aras de la consistencia apliquemos la componente rotacional ω (ver

ecuación (4.9)) al cuadrado unitario obteniendo como resultado los siguientes vectores
de desplazamiento:
~ A = −α̂
∆
~ B = αı̂
∆
~ C = αı̂ − α̂
∆
y el cambio de configuración mostrado en la figura 4.13
B B' C
C'
A
x
A'
Figura 4.13. Configuraciones original y deformada de un cuadrado unitario debido
a la componente anti-simétrica del tensor de transformación de desplazamientos.
Ejemplo
Un tensor dado D permite calcular los desplazamientos ∆ de los puntos en un espacio

(en 2-dimensiones) y cuya posición original se encuentra dada por el vector posición ~r
de acuerdo con
~ = D · ~r .
∆
El este tensor resulta de superponer (o sumar) los 2 efectos mostrados en la figura 4.14a
y figura 4.14b en la cual las lı́neas continuas representan la configuración original y las
lı́neas punteadas la configuración final.
Y
Y
B '(0, 2)
B '(0.5,1)
B(0,1)
B(0,1)
C (1,0)
X X
C '(2,0)
C (1,0)
C '(1, −0.5)
(a) Cambios debidos a la componente (b) Cambios debidos a la componente

simétrica anti-simétrica
Figura 4.14. Transformación lı́neal de elemento triangular.
Se pide:
(i) Determinar el tensor D.
(ii) Determinar las componentes simétrica ε y anti-simétrica ω del tensor D.
(iii) Determinar las direcciones y valores principales de la componente simétrica

ε.
En la primera transformación se tienen los siguientes vectores de posición originales

~rB = 1.0̂
~rC = 1.0ı̂
y los siguientes vectores de desplazamientos para los puntos B y C respectivamente;
~ B = 1.0̂
∆
~ C = 1.0ı̂
∆
planteando para cada punto la ecuación:

∆x εxx γ̄ rx
=
∆y γ̄ εyy ry
permite determinar;
εxx γ̄ 1.0 0.0
ε≡ = .
γ̄ εyy 0.0 1.0
De manera análoga, en el caso de la segunda transformación se tiene que

0 α 0.0 0.5
ω≡ =
−α 0 −0.5 0.0
por lo que el tensor resultante corresponde a:

1.0 0.0 0.0 0.5 1.0 0.5
ε+ω ≡ + ≡ .
0.0 1.0 −0.5 0.0 −0.5 1.0
Problema propuesto
La figura 4.15 corresponde a la configuración original de un medio continuo. Si sobre

este se aplica el tensor de transformación de desplazamientos:

dxx dxy
D=
dyx dyy
determinar la configuración deformada del medio indicando claramente el efecto de cada

una de las componentes definidas de acuerdo con:
D = ε + ω ≡ pI + S + C + ω
E(0,1) B(1,1)
A(−1,1)
G(−1,0) H (1,0)
x
D(−1,−1) F(0,−1) C(1,−1)
Figura 4.15. Transformación al espacio completo.

4.3. Tensor de deformaciones unitarias
Procedamos ahora a comparar lo que sucede con una fibra materiales infinitesimal tras
haber sido sometida al proceso de deformación. Ésta se muestra en la figura 4.16 conec-
tando los puntos materiales que ocupan las posiciones P y Q antes de la deformación y
las posiciones p y q en la configuración deformada. En la configuración original la fibra
P Q esta dada por el vector d~r expresado en términos de magnitud y dirección como:
d~r = n̂ dS
donde
dx dy
n̂ = ı̂ + ̂ .
dS dS
!
dx
Q
!
dr
P p
Figura 4.16. Fibras materiales emanando de un punto material P en las configura-

ciones original y deformada.
En la configuración deformada la fibra P Q se convierte en la fibra pq dada por el

vector d~x. Ambas se comparan en la figura 4.17 en la que también se identifica el
desplazamiento ~uP del punto material P y el desplazamiento ~uQ del punto material Q.
q
!
uQ
Q !
dx
!
dr !
du
! p
P uP
Figura 4.17. Desplazamiento relativo entre 2 puntos materiales infinitesimalmente
cercanos.
Escribiendo el desplazamiento relativo d~u en términos de los vectores fibra original y

deformado d~r y d~x se tiene:
d~u = d~x − d~r
o usando una relación como:
d~x = A · d~r
donde A transforma d~r en d~x se tiene:
d~u = [A − I] · d~r
o, finalmente,
d~u = J · d~r . (4.11)
En la expresión ecuación (4.11) el tensor J 4 permite obtener el desplazamiento re-

lativo d~u del punto material Q con respecto a P . Escribiendo d~u en términos de sus
componentes rectangulares como
d~u = du ı̂ + dv ̂
permite escribir la ecuación (4.11) en forma explicita de acuerdo con:

du Jxx Jxy dx
= .
dv Jyx Jyy dy
4
El tensor J es el análogo al tensor D visto en la sección de trasformaciones lineales
Reconociendo los términos del tensor J como las variaciones de las componentes rec-
tangulares du y dv a lo largo de distancias diferenciales dx y dy permite escribir:
∂u ∂u
du = Jxx dx + Jxy dy ≡ dx + dy
∂x ∂y
∂v ∂v
dv = Jyx dx + Jyy dy ≡ dx + dy
∂x ∂x
o en forma explicita:  
∂u ∂u
du ∂y  dx
= ∂x

dv ∂v ∂v dy
∂x ∂y
y en donde el tensor J corresponde al gradiente de desplazamientos. Normalizando el

vector de desplazamientos relativos por la magnitud original dS de la fibra, permite
obtener el desplazamiento relativo unitario como una transformación
    
 du  ∂u ∂u  dx 
dS
=  ∂x ∂y  dS
 dv  ∂v ∂v  dy 
dS ∂x ∂y dS
o
d~u
= J · n̂. (4.12)
dS
Para identificar los diferentes efectos contenidos en el gradiente de desplazamientos este

se descompone en términos de sus componentes simétrica y anti-simétrica como:
J =ε+ω
o en forma explicita:
     
∂u ∂u ∂u 1 ∂u ∂v 1 ∂u ∂v
∂x 2 ∂y
+ ∂x
0 2 ∂y
− ∂x
 ∂x ∂y 
= + 
∂v ∂v 1 ∂u ∂v ∂v
∂x ∂y 2 ∂y
+ ∂x ∂y
− 12 ∂u
∂y
− ∂v
∂x
0
Los efectos fı́sicos de las diferentes componentes del tensor J serán estudiados ana-
lizando los cambios impuestos por este sobre el elemento diferencial mostrado en la
figura 4.18. En particular, se calcularán los desplazamientos de los puntos A, B y C
con vectores posición:
~rA = dx ı̂
~rB = dy ̂
~rC = dy ı̂ + dx ̂ .
B C
dy
dx A x
Figura 4.18. Elemento diferencial de tamaño dx × dy.
4.3.1. Rotación de cuerpo rı́gido
Estudiando inicialmente la componente anti-simétrica ω del tensor J aplicando la trans-

formación sobre el elemento diferencial se tiene respectivamente para los desplazamien-
tos relativos de los puntos A, B y C:
       
1 ∂u ∂v
du 0 − dx  0 
= 2 ∂y ∂x  ≡
dv  − 12 ∂u ∂v
− ∂x 0  0  − 1 ∂u − ∂v dx
A ∂y 2 ∂y ∂x
   
du  1 ∂u − ∂v dy 
2 ∂y ∂x
=
dv   0 
B
   
du  1 ∂u − ∂v dy 
2 ∂y ∂x
=
dv  − 1 ∂u − ∂v dx
C 2 ∂y ∂x
dando como resultado la configuración deformada de la figura 4.19.
Rigurosamente, el elemento experimenta un pequeño cambio de volumen (área) con-

siderado como de segundo orden. Claramente esta componente no involucra fuerzas
internas o tensiones. De otro lado se debe identificar que la rotación del elemento esta
dada por:
1 ∂u ∂v
ω= − .
2 ∂y ∂x
1 ⎛ ∂u ∂v ⎞
-
2 ⎜⎝ ∂y ∂x ⎟⎠
B B' C
C'
dy
dx A
x
1 ⎛ ∂u ∂v ⎞
-
2 ⎜⎝ ∂y ∂x ⎟⎠
A'
Figura 4.19. Rotación de cuerpo rı́gido de un elemento diferencial.
4.3.2. Componente simétrica ε
Para el punto A
       
∂u 1 ∂u ∂v ∂u
du
∂x 2 ∂y
+ ∂x
dx  dx 
=  ≡ ∂x .
dv  1 ∂u
+ ∂v ∂v  0   1 ∂u + ∂v dy 
A 2 ∂y ∂x ∂y 2 ∂y ∂x
Similarmente, para los puntos B y C,

       
du  1 ∂u + ∂v dy  du  ∂u dx + 1 ∂u + ∂v dy 
2 ∂y ∂x ∂x 2 ∂y ∂x
= , = ,
dv   ∂v
dy  dv   1 ∂u + ∂v dx + ∂v dy 
B ∂y 2 ∂y
C ∂x ∂y
resultando en la configuración mostrada en la figura 4.20.

y C'
B'
B C
dy
A'
dx A x
Figura 4.20. Cambio de forma en el elemento diferencial debido a la componente

simétrica del gradiente de desplazamientos.
Ahora, estudiemos de manera individual los efectos de los distintos componentes del
tensor simétrico ε.
4.3.3. Componente escalar, volumétrica o isotrópica
      
1 ∂u ∂v dx  1 ∂u + ∂v dx
du
2 ∂x
+ ∂y
0 2 ∂x ∂y
=  ≡
dv  1 ∂u ∂v
0 2 ∂x
+ ∂y 0  0 
A
 
0

du
 
=
dv B  1 ∂u + ∂v dy 
2 ∂x ∂y
   
du  1 ∂u + ∂v dx
2 ∂x ∂y
=
dv   1 ∂u + ∂v dy 
C 2 ∂x ∂y
B'
C'
1 ⎛ ∂u ∂v ⎞
+ dy
2 ⎜⎝ ∂x ∂y ⎟⎠ C
B
dy
A A'
dx x
1 ⎛ ∂u ∂v ⎞
+ dx
2 ⎜⎝ ∂x ∂y ⎟⎠
Figura 4.21. Cambio de forma en el elemento diferencial debido a la componente

escalar del gradiente de desplazamientos.
4.3.4. Componente distorsional
      
1 ∂u ∂v dx  1 ∂u − ∂v dx
du
2 ∂x
− ∂y
0 2 ∂x ∂y
=  ≡
dv  0 − 12 ∂u
∂x
∂v
− ∂y 0  0 
 A  
du  0 
=
dv  − 1 ∂u − ∂v dy 
B 2 ∂x ∂y
   
du  1 ∂u − ∂v dx 
2 ∂x ∂y
=
dv  − 1 ∂u − ∂v dy 
C 2 ∂x ∂y
B C
1 ⎛ ∂v ∂u ⎞
- dy C'
2 ⎜⎝ ∂y ∂x ⎟⎠ B'
dy
A A'
dx x
1 ⎛ ∂u ∂v ⎞
- dx
2 ⎜⎝ ∂x ∂y ⎟⎠
Figura 4.22. Cambio de forma en el elemento diferencial debido a la primera com-

ponente distorsional del tensor de deformaciones unitarias.
       
1 ∂u ∂v
du 0 2 ∂y
+ ∂x
dx  0 
=  ≡
dv  1 ∂u ∂v
+ ∂x 0  0   1 ∂u + ∂v dx
A 2 ∂y 2 ∂y ∂x
   
du  1 ∂u + ∂v dy 
2 ∂y ∂x
=
dv   0 
B
   
du  1 ∂u + ∂v dy 
2 ∂y ∂x
=
dv   1 ∂u + ∂v dx
C 2 ∂y ∂x
y
1 ⎛ ∂u ∂v ⎞
+
2 ⎜⎝ ∂y ∂x ⎟⎠ C'
B B' C
dy
A' 1 ⎛ ∂u ∂v ⎞
+
2 ⎜⎝ ∂y ∂x ⎟⎠
dx A x
Figura 4.23. Cambio de forma en el elemento diferencial debido a la segunda com-

ponente distorsional del tensor de deformaciones unitarias.
Ahora, en el argot ingenieril es común a cada término de la componente simétrica

del tensor gradiente de desplazamientos asignarle una nomenclatura particular. Para
generalizar la definición consideremos el tensor ε expresado en coordenadas cartesianas.
 

∂u 1 ∂u ∂v 1 ∂u
+ ∂w

 ∂x 2 ∂y
+ ∂x 2 ∂z ∂x 
 
 
[ε] =  1 ∂v ∂u
 2 ∂x + ∂y
∂v 1 ∂v
+ ∂w ,
∂y 2 ∂z ∂y 
 
 
1 ∂w ∂u 1 ∂w ∂v ∂w

2 ∂x
+ ∂z 2 ∂y
+ ∂z ∂z
o, escrito de otra forma:  

1 1
 εxx γxy γxz 
 2 2 
1 1 
[ε] = 
 γxy εyy γyz 
2 2 

1 1 
γxz γyz εzz
2 2
En donde εxx , εyy y εzz son conocidas como deformaciones axiales o normales y γxy ,γxz
y γyz como deformaciones angulares o deformaciones por cortante.
Por otro lado, la componente rotacional podrı́a escribirse

 
0 ωxy ωxz
[ω] = −ωxy 0 ωyz  ,
−ωxz −ωyz 0
en donde ωxy , ωxz y ωyz representan la rotación de cuerpo rı́gido respecto a los ejes z,
y y x respectivamente.
Para terminar, el tensor ε en coordenadas cilı́ndricas está definido como:

 
1 1
 εrr γrθ γrz 
 2 2 
1 1 
[ε] = 
 γrθ εθθ ,
γθz 
2 2 
1 1 
γrz γθz εzz
2 2
en donde
∂u u 1 ∂v ∂w
εrr = , εθθ = + , εzz = ,
∂r r r ∂θ ∂z
1 ∂u ∂v v ∂v 1 ∂w ∂u ∂w
γrθ = + − , γθz = + , γrz = + .
r ∂θ ∂r r ∂z r ∂θ ∂z ∂r
Ejemplo
El campo de desplazamientos para la cuña auto-soportada mostrada en la figura 4.24

esta dado por:
S S
u= K1 (ν, φ)(x − ` cos φ) v = − K2 (ν, φ)y
E E
donde u y v corresponden a los desplazamientos en las direcciones horizontal y vertical
respectivamente y K1 (ν, φ) y K2 (ν, φ) son constantes que dependen del semi-ángulo de
la cuña y del material. E es una constante que depende del material.
l l
S
x
S
S
l l
Figura 4.24. Cuña auto-soportada.
• Esbozar la configuración deformada (final) de la cuña
• Los tensores gradiente de desplazamientos J, deformaciones unitarias y rotacio-

nes ω.
• Determinar las direcciones de máximo cambio de tamaño y máxima rotación.
Para esbozar la configuración deformada de la cuña es preciso evaluar el campo de

desplazamientos en puntos de la misma. Debido a la linealidad de las funciones del
campo de desplazamientos u y v y por la simetrı́a del problema miramos lo que pasa
en cada una de las esquinas 5 . Es decir, evaluamos los desplazamientos en los puntos
A = (0, 0), B = (` cos φ, ` sen φ), C = (` cos φ, −` sen φ) y D = (2` cos φ, 0). En la
figura 4.25 se muestra esquemáticamente en linea punteada la configuración deformada
de la cuña. Nótese como los puntos ubicados en el eje de simetrı́a A − D solo tienen
desplazamiento en x, u, mientras que los puntos en el eje de simetrı́a B − C solo en y,
v. El punto central del medio continuo no se desplaza.
5
Para funciones de desplazamiento más complejas puede ser preciso usar métodos más elaborados
para la evaluación de dichas funciones.
a A D d
c
C
Figura 4.25. Configuración final Cuña auto-soportada.
Por otro lado, el tensor gradiente de desplazamientos está dado por:

 
S K1 (ν, φ) 0
J=  
E 0 −K2 (ν, φ)
y, en este caso, se tiene que ε ≡ J.
Luego, las rotaciones de cuerpo rı́gido son nulas
ω = 0.
Nótese que el tensor ε ya está escrito en sus valores principales, lo que quiere decir que
las deformaciones axiales máximas están dadas por: ES K1 (ν, φ) como el estiramiento
máximo, y ES K2 (ν, φ) como el acortamiento máximo.
4.3.5. Ejemplo: Barra sometida a su peso propio
En la figura 4.26 se muestra una barra de sección circular de radio R, altura H, fija
en su extremo inferior y que está sometida solamente a la acción de su propio peso. Se
indica el campo de desplazamientos en el sistema coordenado xyz.

γ γ γ 2 γν 2 γ 2
u=ν zx , v=ν zy , w=− z − (x + y 2 ) + H ,
E E 2E 2E 2E
donde u, v, w son los desplazamientos asociados a los ejes x, y, z, respectivamente. E, γ
y ν son constantes del material.
R X
H
Y Z
Figura 4.26. Barra apoyada
• Esboce la configuración deformada de la barra.
Solución:
Para esbozar una configuración deformada lo hacemos a partir del campo de

desplazamientos u, v, w. Para realizarla hagamos un corte transversal (plano xy)
y un corte longitudinal (plano xz) y evaluemos que pasa con los desplazamientos
en algunos puntos de la frontera.

γ
x = 0, y = 0, z = 0. =⇒ u = 0, v = 0, H w=
2E
γν γ 2
x = ±R, y = 0, z = 0. =⇒ u = 0, v = 0, w = − R2 + H
2E 2E
γν γ 2
x = 0, y = ±R, z = 0. =⇒ u = 0, v = 0, w = − R2 + H
2E 2E
x = 0, y = 0, z = H. =⇒ u = 0, v = 0, w = 0
γ γν
x = ±R, y = 0, z = H. =⇒ u = ±ν HR, v = 0, w = − R2
E 2E
γ γν
x = 0, y = ±R, z = H. =⇒ u = 0, v = ±ν HR, w = − R2
E 2E
En la figura figura 4.27 aparecen la configuración inicial (lı́nea punteada) y final

del elemento.
X Y
R
X
Z
Figura 4.27. Cortes con configuración final
• Calcule el tensor gradiente de desplazamientos y descompóngalo en [ε] (simétrico)

y [ω] (asimétrico).
Solución:
El tensor gradiente de desplazamientos está dado por:

 
νz 0 νx
γ
[D] =  0 νz νy  .
E
−νx −νy −z
El tensor ε y ω están dados por:

   
νz 0 0 0 0 νx
γ  γ
[ε] = 0 νz 0  , [ω] =  0 0 νy  .
E E
0 0 −z −νx −νy 0
• ¿Cuáles son los valores máximos y mı́nimos de las deformaciones axiales y en que
puntos se presentan?
Solución:
γ γ
εxx = ν z =⇒ εmı́n = 0 (z = 0), εmáx = ν H (z = H)
E E
γ γ
εyy = ν z =⇒ εmı́n = 0 (z = 0), εmáx = ν H (z = H)
E E
γ γ
εzz = − z =⇒ εmı́n = 0 (z = 0), εmáx = − H (z = H)
E E
• ¿Cuáles son los valores máximos y mı́nimos de los desplazamientos u, v y w y en

que puntos se presentan?
Solución:
γ γ
u = ν zx =⇒ umı́n = 0 (x = 0), umáx = ν HR (z = H, x = |R|)

E γE
γ
v = ν zy =⇒ vmı́n = 0 (y = 0), vmáx = ν HR (z = H, y = |R|)

E E
Para el caso del desplazamiento en z los máximos dependen de la relación entre

H 2 y R2 = x2 + y 2 . Supongamos H > R.
γ 2
wmı́n = 0 (x = y = 0, z = H), wmáx = H (x = y = z = 0)
2E
• ¿Existen puntos exentos de rotación de cuerpo rı́gido?, responda sı́ o no y justifique

su respuesta.
Solución:
Sı́, si existen puntos exentos de rotación. La lı́nea longitudinal cuando x = y = 0

no experimenta rotación de cuerpo rı́gido.
• ¿Es posible encontrar puntos sometidos a deformación angular?, responda sı́ o no

y justifique su respuesta.
Solución:
Sı́ es posible. Como puede verse en el tensor ε los valores principales son diferentes,
por lo es posible encontrar alguna dirección sometida a deformación angular.
Tarea
Describa la configuración deformada de las partı́culas.
4.4. Ejercicios propuestos
1. En la figura 4.28 se muestra la configuración inicial y final en donde √

las coordena-
√
das puntos
√ son A,
√ B, C y D son A = (1, 1), B = (−1, −1), C = (− 2/2, 2/2),
D = ( 2/2, − 2/2). Sabiendo que ~r es el vector de posición original, ρ~ es el
~ es vector de desplazamiento
vector de posición final y M
~ = D~r.
a) Determinar la matriz de desplazamientos lineales [D] tal que M
b) Determinar la matriz de transformación lineal [A] tal que ρ~ = A~r

c) Descomponer [D] en componente simétrica y componente asimétrica. Mos-

trar gráficamente lo efectos.
d ) Usando el circulo de Mohr, determinar las direcciones principales y los valores

propios de la matriz de transformación lineal simétrica.
Configuración inicial
Configuración final
Figura 4.28. Representación gráfica de una transformación lineal.
2. En la figura figura 4.29 se muestra la configuración inicial y final (linea punteada)

del un triángulo de lado unitario formado por los puntos A, B, C. Para pasar de
la configuración inicial a la final se aplica una trasformación lineal de desplaza-
mientos [D], tal que {∆} = [D] · {r}, donde ∆ es el vector desplazamiento, y r es
el vector posición inicial.
B '(1,2)
Y
B(0,1)
A(0,0)
C(1,0) C '(2,0) X

Figura 4.29. Configuraciones inicial y deformada
a) Determinar la transformación lineal [D], y descomponerla en sus componente

simétrica [ε] y asimétrica [ω].
b) Determinar las direcciones principales de la transformación [D]. Es posible

usar el cı́rculo de Mohr.
c) Mostrar matemática y gráficamente el efecto independiente de las diferentes

componentes del tensor lineal de desplazamientos. Las componentes a con-
siderar son la isotrópica, las distorsionales y la rotación de cuerpo rı́gido, tal
que [D] = p[I] + [S] + [C] + [ω]
d ) Escribir en sus direcciones principales la componente simétrica [ε] de la trans-

formación [D].
3. La figura figura 4.30 muestra un rectángulo infinitesimal de lado b y altura h en

la configuración sin deformar el cual es deformado por alguna acción externa y
dicha deformación se describe como la superposición lineal de varias deformaciones
individuales, las cuales se dan a continuación:

d 0
a) ε1 =
0 d

e 0
b) ε2 =
0 −e

0 f
c) ε3 =
f 0
Figura 4.30. Rectángulo en la configuración no deformada.
Se pide determinar lo siguiente:
a) Calcular los desplazamientos de los puntos A, B y C para cada uno de los

estados de deformación y dibujar las configuraciones deformadas.
b) Si se sabe que las dimensiones del rectángulo b = 1.00 y h = 1.00 y los valores
de las transformaciones son d = 0.01, e = 0.005 y f = 0.0025. Calcular las
direcciones principales y los valores propios para la superposición lineal de
las deformaciones.
4. En la figura figura 4.31 se muestra una barra de sección circular de radio R, altura
H, soportada de su extremo superior y que está sometido solamente a la acción
de su peso propio. El campo de desplazamientos u, v, w en el sistema coordenado
xyz de para la barra está dado por:
γ γ γ 2 γν 2 γ 2
u = −ν zx , v = −ν zy , w= z + (x + y 2 ) − H ,
E E 2E 2E 2E
donde E, γ y ν son constantes del material (positivas).
H
Z
Figura 4.31. Barra colgada
a) Esboce la configuracióon deformada de la barra.
b) Calcule el tensor de deformaciones y descompóngalo en su componente simétri-

ca: [ε] y asimétrica [ω]
c) Cuáles son los valores máximos y mı́nimos de las deformaciones axiales, εxx ,
εyy , εzz , y en que puntos se presentan.
d ) Cuáles son los valores máximos y mı́nimos de los desplazamientos: u, v, w y

en que puntos se presentan.
e) ¿Existen puntos exentos de rotación de cuerpo rı́gido?, responda sı́ o no y

justifique su respuesta.
f ) ¿Es posible encontrar puntos sometidos a distorsión angular?, responda sı́ o

no y justifique su respuesta.
5. En la figura figura 4.32 se muestra un medio continuo de sección rectangular de

ancho unitario, altura H, soportado en un extremo (X = 0), y sometido en su
otro extremo (X = L) a la acción de un momento flector M alrededor del eje Z.
Si el campo de desplazamientos en el sistema coordenado XY Z está dado por:
M M νM
u=− XY , v= (X 2 + ν(Y 2 − Z 2 )) , w= YZ,
EI 2EI EI
donde u, v, w son los desplazamientos asociados a los ejes X, Y, Z respectivamente.

E y ν son constantes (positivas) del material e I es el momento de inercia alrededor
del eje Z de la sección transversal.
H L
X
H/2
Z
0.5
Y
Figura 4.32. Viga sometida a momento flector
a) En la figura figura 4.33 se muestran esbozos de posibles configuraciones de-

formadas (sección achurada) de la viga mostrada en la figura . Determine en
cada caso cual es la configuración correcta.
• Para la sección transversal en X = L.
• Para una sección longitudinal cuando Y > 0.
• Para una sección longitudinal asociada al plano (Y = −H/2)
• Para una sección longitudinal asociada al plano (Z = 0)

Y Y
Z Z Z Z Z Z
Z Z
a Y b Y c d Y e Y f Y g h Y
Z Z Z
X X X
i j k
Z Z Z
X X X
l m n
Z Z
X X X
o p Y q
t
X X X
Y r Y s Y
Figura 4.33. Posibles configuraciones deformadas de sección transversal
b) Calcule el tensor de deformaciones [ε] y el de rotación [ω]
c) Si el valor de la distorsión angular máxima es γ = γmáx determine cual es el

momento máximo posible, Mmáx que podrı́a ser aplicado.
d ) ¿Cuáles son las coordenadas x, y, z de los puntos que experimentan mayores

rotaciones de cuerpo rı́gido?
e) ¿Cuáles son los puntos que experimentan la menor rotación de cuerpo rı́gido?
f ) Cuales son los valores máximos y mı́nimos de las deformaciones axiales, εxx ,
εyx , εzz y en que puntos se presentan.
g) ¿Cuáles son los valores máximos y mı́nimos de los desplazamientos, u, v, w,

y en qué puntos se presentan?
h) ¿Existen puntos exentos de rotación de cuerpo rı́gido?, responda sı́ o no y

justifique su respuesta.
i ) ¿Es posible encontrar puntos sometidos a distorsión angular?, responda sı́ o

6. En la figura figura 4.34 se muestra la solución de un semi-espacio sometido a una

carga lineal superficial P . Los desplazamientos al interior del suelo debidos a la
carga P están dados por:
2P (1 − ν)P
ur (r, θ) = − ln r cos θ − θ sen θ + B cos θ
πE πE
2νP 2P (1 − ν)P (1 − ν)P
uθ (r, θ) = sen θ + ln r sen θ − θ cos θ + sen θ − B sen θ
πE πE πE πE
w(r, θ) = 0
donde E es el módulo de elasticidad, ν es la relación de Poisson del material y B

es una constante positiva.
P
y
r
Figura 4.34. Solución carga vertical
Si en un suelo sometido a la acción de dos estructuras que le transmiten una

carga lineal superficial, P y 3P , respectivamente tal como se muestra en la figura
figura 4.35 se desea instalar una tuberı́a, perpendicular al plano mostrado (XY ),
de diámetro despreciable, en el punto indicado en la figura figura 4.35.
P 3P
y
d
Figura 4.35. Estructuras y tuberı́a.
a) Determine cuál es el valor de la deformación axial máxima.
b) Si la máxima deformación angular a la que puede ser sometida la tuberı́a sin

fallar, es igual a: γmáx . ¿Cuál es el máximo valor de la carga distribuida P
que pueden transmitir las estructuras para que la tuberı́a no falle?
c) Esboce la configuración deformada de la partı́cula correspondiente al punto

donde se ubica la tuberı́a. Solo considere el plano xy.
7. En la figura figura 4.36 se muestra una barra de radio R y longitud L. El campo

de desplazamientos en el sistema coordenado xyz, está dado por:
u = −θyz , v = θxz , w = 0.
donde u, v, w son los desplazamientos asociados a los ejes x, y, z respectivamente.

θ es una constante positiva.
R X
H
Y Z
a) Dibuje la configuración deformada de la barra.

b) Calcule el tensor de deformaciones, [D], y determine: tensor de deformación

[ε] (simétrico) y tensor de rotación [ω] (asimétrico).
c) Ilustre la deformación de la partı́cula en el punto de coordenadas y = R y

x = 0. Para la ilustración use cuadrados de tamaño diferencial contenidos en
los planos xy, xz, yz de forma independiente. En la ilustración debe incluirse
los efectos de cuerpo rı́gido.
d ) ¿Cuáles son los puntos que experimentan mayores rotaciones de cuerpo rı́gi-
do?
e) ¿Cuáles son los valores máximos y mı́nimos de los desplazamientos, u, v, w

y en que puntos se presentan:
f ) ¿Es posible encontrar puntos sometidos a distorsión angular?, responda sı́ o

8. En la figura figura 4.37 se presenta una cuña de espesor t sometida a la acción de

una carga (en el plano XY ) distribuida sobre su perı́metro.
Y
l
Figura 4.37. Cuña de espesor t
El campo de desplazamientos de la cuña esta dado por:

2B 2B
u(x, y) = (cot φ + ν tan φ)(x − l cos φ) , v(x, y) = − (tan φ + ν cot φ)y ,
E E
donde B, E y ν son constantes positivas.
Además, se sabe que:

E E E
σxx = (εxx + νεyy ) , σyy = (εyy + νεxx ) , τxy = γxy .
(1 − ν 2 ) (1 − ν 2 ) 2(1 + ν)
a) Dibuje la configuración deformada de la cuña.
b) Ilustre la deformación de la partı́cula en el punto de coordenadas X = 0 y

Y = 0. En la ilustración deben incluirse los efectos de cuerpo rı́gido.
c) Si ν = 0.25 y φ = 53.13◦ , determine la reacciones (Fuerzas) en el soporte .
9. Si el campo de desplazamientos para la barra mostrada en la 4.38 está dado por:

1 − ν2 1+ν
u (x) = − σ0 x , v = 0, w (z) = ν σ0 z .
E E
t
y
H
L
z H/2
x
t/2
Figura 4.38. Viga de Longitud = L, Altura = H y Espesor = t. El origen del sistema

coordenado coincidiendo con el centroide de la barra.
a) Explique a que acciones externas está sometida la barra.
b) Dibuje la configuración deformada de la barra.
c) ¿Hay partı́culas sometidas a rotaciones de cuerpo rı́gido? (responda sı́ o no

y justifique su respuesta).
d ) ¿Hay partı́culas que NO experimenten distorsiones angulares? (responda sı́

o no y justifique su respuesta).
e) En que coordenadas está la partı́cula que sufre mayor deformación axial y

cuál es el valor de esta deformación.
f ) En que coordenadas está la partı́cula que sufre mayor distorsión angular y

cuál es el valor de esta deformación.
g) Encuentre los desplazamientos y deformaciones del punto de coordenadas

(0, 0, 0).
10. Para un campo de desplazamientos dado por:
u = a(x2 − 5y 2 )ı̂ + (2axy)̂
a) Determinar el tensor de deformación;
b) Obtener las deformaciones principales.
11. En la figura figura 4.39 se presenta una barra sobre la cual se aplica una fuerza en
dirección x Fx , que genera un esfuerzo constante σx sobre la sección transversal
donde está aplicada. La barra inicialmente está separada de una lámina delgada
de acero una distancia d tal y como se muestra en figura 4.39b. El campo de
desplazamientos para la barra mostrada está dado por:
1 1
u (x) = xσx , v = 0, w (z) = − z(H 2 − y 2 )σx ,
E1 E2
donde E1 = 200000 kgf/cm2 , y E2 = 1000000 kgf/cm4 , L = 1000 cm, H = 10 cm.

σx es el esfuerzo en dirección x el cual tiene una variación en el tiempo de acuerdo
a como lo muestra la figura 4.39c.
x 2000
Fx Lámina 1800
x 1600
d
A 1400
H H
σx (kgf/cm2 )
H H 1200
L
L 1000
2L 800
z
600
y z
L 400
L
200
00 100 200 300 400 500
2H 2H
H H Tiempo (s)
(a) Configuración (b) Vista plano (c) Curva de esfuerzo.
inicial de la barra xz. Separación
en 3D. Aplicación inicial barra -
de la fuerza lámina, d = 2.0
cm. φ = 45◦ .
Figura 4.39. Barra sometida a fuerza axial Fx
a) Para el punto A de la barra con coordenadas (x, y, z) = (L, 0, B) grafique

los desplazamientos u y w en función de tiempo.
b) Dibuje la configuración deformada de la partı́cula para el punto de coorde-

nadas (x, y, z) = (0, 0, 0)
c) Si el valor de σx = 50 kgf/cm2 y se estudian solo los puntos de la cara z =

1.0cm determinar el valor de la distorsión angular máxima y las coordenadas
x, y, z de los puntos en que se presenta.
d ) Si el esfuerzo σx se incrementa hasta que el punto A con coordenadas (x, y, z) =

(L, 0, B) toque la lámina, determine el tiempo t en el que la barra toca la
lámina. ¿Cuáles serı́an las coordenadas finales x, z del punto A en ese ins-
tante?
12. En un ensayo de un elemento de concreto se mide la deformación en un mismo

punto mediante tres galgas de deformación. Las galgas de deformación se disponen
en un arreglo tal y como se muestra en la figura 4.40.
y y x'
y''
0.0004
y' x''
0.0002
60º 60º
x x
- 0.0003
Figura 4.40. Galgas de deformación
Si las galgas solo registran la deformación axial y las mediciones de laboratorio

reportaron: εxx = −0.0003; εx0 x0 = 0.0004; εy00 y00 = 0.0002. Se solicita
a) Tensor de deformaciones en el sistema xy, x0 y 0 y x00 y 00 .
b) Distorsión angular máxima.
c) Deformación axial máxima.
d ) Configuración deformada de la partı́cula.

Capı́tulo 5
Elasticidad
5.1. Introducción
Al final de esta sección el estudiante deberá estar en la capacidad de
• Entender las hipótesis básicas de la elasticidad lineal.
• Reconocer la relación de tensión-deformación para materiales elásticos lineales.
• Identificar los efectos fı́sicos relacionados con las constantes de material como el
módulo de Young y el coeficiente de Poisson.
• Identificar las consideraciones implicadas en las idealizaciones de esfuerzo plano

y deformación plana.
• Operar en problemas de elasticidad simples.
Como se discutió en el capı́tulo 5 las condiciones de equilibrio en términos de tensiones

se reducen a 3 ecuaciones diferenciales parciales en las 9 componentes del tensor de
tensiones y a 3 relaciones de simetrı́a resultantes del equilibrio rotacional del punto
material: en total estas correspondı́an a 6 ecuaciones en 9 incógnitas. Al introducir los
cambios cinemáticos en términos del tensor de deformaciones unitarias se involucran en
164
CAPÍTULO 5. ELASTICIDAD 165
el problema las siguientes 6 ecuaciones adicionales y 9 incógnitas más

∂u ∂v ∂w
εxx = , εyy = , εzz = ,
∂x ∂y ∂z

1 ∂u ∂v 1 ∂u ∂w 1 ∂v ∂w
εxy = + , εxz = + , εyz = + ,
2 ∂y ∂x 2 ∂z ∂x 2 ∂z ∂y
con lo que se tienen 12 ecuaciones en 18 incógnitas.
Aunque todos los problemas abordados hasta el momento han estado relacionados con
sólidos, también es cierto que en ningún momento hemos hecho distinción alguna de un
material a otro. Es decir las soluciones en términos de tensiones para problemas como el
de la cuña autosoportada son igualmente válidas independientemente del material del
que este fabricada la cuña. En esta sección se presentan relaciones entre las tensiones
y las deformaciones las cuales cambiarán necesariamente de un material a otro. Por
ejemplo, refiriéndonos al problema de la cuña, es natural esperar que si 2 cuñas con
iguales geometrı́as y sometidas a los mismos cortantes de magnitud S, pero fabricadas
de 2 materiales diferentes, por ejemplo acero y madera, éstas se deformen en magnitudes
diferentes. De forma general podemos escribir estas relaciones como
    

 σxx 
 C11 C12 C13 C14 C15 C16   εxx 
σ C C C C C C ε
   
yy 21 22 23 24 25 26 yy

 
   
 

     
σzz C C C C C C ε
    
31 32 33 34 35 36  zz

=  . (5.1)

 τyz 
 C41 C42 C43 C44 C45 C46   
 εyz 
τxz  C51 C52 C53 C54 C55 C56   ε 
    
 xz 

   
  
τxy C61 C62 C63 C64 C65 C66 εxy
  
Con la ecuación (5.1) se tienen ahora 18 ecuaciones en 18 incógnitas con lo cual la

solubilidad o no del problema depende ahora de las condiciones especificadas en la
frontera.
5.2. Hipótesis básicas
En esta sección se particulariza el modelo del medio continuo al modelo más simple: un
sólido elástico lineal. Las hipótesis básicas para la teorı́a de elasticidad lineal son:
El material es elástico. Esto implica que el cuerpo retorna a su estado inicial

cuando las fuerzas externas dejan de actuar sobre él. Esto implica que el trabajo
realizado sobre el cuerpo se almacena como energı́a de deformación y que no hay

disipación de energı́a en este proceso.
Los materiales son isótropos. Esto implica que las propiedades de los mate-
riales no dependen de la dirección en la que se midan. Por ejemplo, la madera
no es un material isótropo, pues las propiedades mecánicas son diferentes en la
dirección de la fibra que perpendicular a esta.
Los materiales son homogéneos. Esto implica que las propiedades de los
materiales son iguales en todos los puntos del medio.
Los materiales son lineales. Un material es lineal si las deformaciones pue-

den escribirse como una relación lineal de las tensiones, es decir, si satisface la
ecuación (5.1).
Las deformaciones son pequeñas. Esto implica que las ecuaciones de equilibrio
son válidas en la configuración inicial del medio.
Estas consideraciones hacen que la ecuación (5.1) se simplifique enormemente, y que la

cantidad de constantes involucradas se reduzcan de 36 a 2.
5.3. Relación tensión-deformación
Considerando las hipótesis planteadas en el párrafo anterior las relaciones de la ecua-

ción (5.1) se reducen a las relaciones tensión - deformación dadas en la ecuación (5.2):
E νE
σxx = εxx + (εxx + εyy + εzz ) ,
(1 + ν) (1 + ν)(1 − 2ν)
E νE
σyy = εyy + (εxx + εyy + εzz ) ,
(1 + ν) (1 + ν)(1 − 2ν)
E νE (5.2)
σzz = εzz + (εxx + εyy + εzz ) ,
(1 + ν) (1 + ν)(1 − 2ν)
τxy = Gγxy ,
τxz = Gγxz ,
τzy = Gγyz ,
y las relaciones deformación - tensión dadas en la ecuación (5.3):

1
εxx = [σxx − ν(σyy + σzz )] ,
E
1
εyy = [σyy − ν(σxx + σzz )] ,
E
1
εzz = [σzz − ν(σxx + σyy )] ,
E (5.3)
τxy
γxy = ,
G
τxz
γxz = ,
G
τzy
γzy = .
G
En la ecuación (5.2) y en la ecuación (5.3), conocidas comúnmente como la Ley de

Hooke, las constantes E (módulo elástico o de módulo de Young) , G (módulo elástico
de cortante) y ν (relación de Poisson) son propias de cada material. El valor de dichas
constantes se puede determinar mediante ensayos de laboratorio de tensión simple para
el módulo elástico E y relación de Poisson ν y de torsión pura para el módulo de
cortante G. Por otro lado es posible hallar la siguiente relación entre constantes
E
G= .
2(1 + ν)
5.4. Idealizaciones en 2D
En general, es difı́cil resolver problemas en 3D. Por esta razón, existen varias idealiza-
ciones que consideran el medio continuo como si este fuera plano. Estas simplificaciones
consideran que:
• el medio está acotado por dos superficies planas paralelas al plano xy;
• la sección transversal del medio es uniforme en la dirección z;
• la carga se distribuye uniformemente en la dirección z; y
• no hay tracciones tangenciales en las superficies planas que acotan el medio.

Debemos resaltar que la tercera restricción implica que no habrá torsión. Estas restric-
ciones no implican, sin embargo, que el medio continuo debe ser delgado.
Las variables involucradas en las idealizaciones planas son reducidas a 2 desplazamientos

y 3 esfuerzos/deformaciones. Especı́ficamente:
Desplazamientos: u, v;
Deformaciones: εxx , εyy , γxy ;
Esfuerzos: σxx , σyy , τxy ;
en coordenadas cartesianas y
Desplazamientos: ur , uθ ;
Deformaciones: εrr , εθθ , γrθ ;
Esfuerzos: σrr , σθθ , τrθ ;
en coordenadas polares.
5.4.1. Tensión plana
Decimos que un medio continuo está en tensión plana (o esfuerzo plano) si el estado
de esfuerzos es tal que una de sus componentes principales es cero en todo el medio.
Podemos elegir esta dirección principal alineada con z, lo que implica σzz = σxz =
σyz = 0 en todo el cuerpo. Esta idealización es común para placas delgadas con cargas
paralelas a las caras. En algunos casos, se puede considerar como tensión plana placas
levemente curvadas. Este es el caso de un cilindro con pared delgada bajo la acción de
una presión, como en un tanque.
Las relaciones deformación-tensión en este caso se reducen a

ν
εzz = − (σxx + σyy )
E
1
εxx = (σxx − νσyy )
E (5.4)
1
εyy = (σyy − νσxx )
E
τxy
γxy =
G
Note que εzz 6= 0, sin embargo, puede prescindirse de él para el análisis, pues puede
obtenerse a partir de los esfuerzos/deformaciones en el plano.
5.4.2. Deformación plana
Decimos que un medio continuo está en deformación plana si el estado de deformaciones

es tal que una de sus componentes principales es cero en todo el medio. Podemos elegir
esta dirección principal alineada con z, lo que implica εzz = εxz = εyz = 0 en todo el
cuerpo. Esta idealización es común para cuerpos que tienen una dimensión muy grande
comparada las demás.
Las relaciones de deformación esfuerzo deformación para las componentes axiales serı́an
entonces
1
εxx = [σxx − ν(σyy + σzz )]
E
1
εyy = [σyy − ν(σxx + σzz )]
E
1
0 = [σzz − ν(σxx + σyy )]
E
de donde
σzz = ν(σxx + σyy )
Eliminando σzz se tiene que:

1 − ν2 ν(1 + ν)
εxx = σxx − σyy
E E
1 − ν2 ν(1 + ν) (5.5)
εyy = σyy − σxx
E E
τxy
γxy =
G
Reescritura como tensión plana
Es posible escribir las formulaciones de tensión plana y deformación plana de manera

equivalente. En el caso de deformación plana, si se introduce los siguientes módulos
equivalente
E
Ē =
1 − ν2
y
ν
ν̄ = ,
1−ν
podemos escribir las relaciones de deformación-tensión como
ν̄
εzz = − (σxx + σyy )
Ē
1
εxx = (σxx − ν̄σyy )
Ē
1
εyy = (σyy − ν̄σxx )
Ē
τxy
γxy = ,
G
note que el módulo de cortante G no se ve afectado por esta transformación.
5.5. Ejemplos resueltos
Ejemplo 1: Cuña autosoportada
Regresando al ejemplo de la cuña auto-soportada (figura 5.1)

l l
S
x
S
S
l l
Figura 5.1. Cuña auto-soportada.
Usando argumentos de equilibrio determinamos que el campo de tensiones solución del

problema estaba dado por:
σxx = S cot φ
σyy = −S tan φ (5.6)
τxy = 0
Usando las relaciones tensión-deformación para tensión plana (o deformación plana con
módulos efectivos) se tiene que:
1
εxx = [S cot φ + νS tan φ]
E
1
εyy = − [S tan φ + νS cot φ]
E
Recordando las relaciones:

{d~u} = [ε]{d~x} + [ω]{d~x
las cuales para el caso plano se simplifican a
γxy ωxy
du = εxx dx + dy − dy
2 2
γxy ωxy
dv = εyy dy + dx + dx
2 2
y reconociendo que por la simetrı́a del problema las componentes rotacionales son nulas
se tiene que
Z
1 1
u= [S cot φ + νS tan φ] dx ≡ [S cot φ + νS tan φ]x + A ,
E E
y Z
1 1
v=− [S tan φ + ν cot φ] dy ≡ − [S tan φ + νS cot φ]y + B .
E E
donde A y B son constantes de integración.
Ejemplo 2: Barra sometida a carga axial
En la figura 5.2a se muestra una barra de sección rectangular con a = 40 cm, b = 50 cm

y longitud inicial Lc = 500 cm, que es sometida a un ensayo de compresión confinada
(está contenida en un recipiente indeformable), mientras que en la figura 5.2b se muestra
una barra de sección circular con diámetro d = 50 cm y longitud inicial Lt = 497 cm
sobre la cual se realiza un ensayo de tracción simple. Las constantes del material, tanto
para la barra rectangular como para la barra circular, son E = 150000 kgf/cm2 y
ν = 0.20.
Ambos ensayos se realizan con su extremo fijo en el plano z = 0 y están cargados con
un sistema que no transmite esfuerzos cortantes y que hace que el esfuerzo axial se
distribuya uniformemente sobre la superficie transversal de las barras. Los esfuerzos
que se inducen a las barras varı́an con el tiempo conforme a lo mostrado en la figura
figura 5.2c
525 Esfuerzo barra circular

z
350
σz (kgf/cm2 )
F z
F 175
00 350 700 1050 1400

Tiempo (s)
0 Esfuerzo barra rectangular
250
Lc
Lt
σz (kgf/cm2 )
500
750
1000
a b d
x y x y 12500 200 400 600
Tiempo (s)
800 1000 1200 1400
(a) Barra (b) Barra (c) Curvas de esfuerzo

compresión tracción
Figura 5.2. Esquema ensayos y curvas de aplicación carga
Determine:
¿Cuál es el tiempo t en el cual las dos barras tienen la misma longitud?
Solución:
La barra circular está solamente sometida a esfuerzo axial σzz , por lo que la
deformación es
σzzt
εzzt =
E
y la longitud final serı́a:
Lfcir = Linicir (1 + εzzt )
σzzt (5.7)
Lfcir = Linicir 1 +
E
Por otro lado, la barra cuadrada está sometida a esfuerzo axial σxx , σyy y σzz ,
pero solo se deforma en el eje z que es la deformación de interés. Para este caso
entonces,
σzzc
εzzc =
(2µ + λ)
donde
E νE
µ= y λ= .
2(1 + ν) (1 + ν)(1 − 2ν)
Llamando A = (2µ + λ) se tiene que la longitud final para la barra cuadrada serı́a:
Lfcua = Linicua (1 + εzzc )

σzzc (5.8)
Lfcua = Linicua 1 +
A
Resta, entonces, igualar las longitudes y despejar el tiempo. Para hacer esto,
es preciso tener en cuenta barra circular es cargada solo hasta que el esfuerzo
es σ = 350 kgf/cm2 por lo que habrá que considerar ambas condiciones. Por
simplicidad inicialmente determinemos que pasarı́a si t > 350 s, es decir en donde
σzzt = 350 kgf/cm2 . Reemplazando el valor de los esfuerzos en función del tiempo
en la ecuación (5.7) y la ecuación (5.8)

350
Lfcir = 497 1 + (5.9)
E

t
Lfcua = 500 1 − (5.10)
A
donde A ' 1.1111E.
Igualando la ecuación (5.9) y la ecuación (5.10) y despejando t se obtiene
t = 613.4 s .
En la barra circular, ¿cuáles son los desplazamientos y las deformaciones en el

punto de coordenadas (x, y, z) = (20, 15, 450) cm para el tiempo t = 300 s?
Las deformaciones asociadas a la barra son:

σzzt
εxxt = −ν
E
σzzt
εyyt = −ν (5.11)
E
σzzt
εzzt =
E
reemplazando σzzt por 300
εxxt = −0.0004
εyyt = −0.0004 (5.12)
εzzt = 0.002
Por otro lado para determinar los desplazamientos integramos la deformación. Por
ser deformación constante la integración se obtiene multiplicando la deformación
por las distancias hasta el punto de evaluación.
u = −0.0004(20 cm) = −0.008 cm

v = −0.0004(15 cm) = −0.006 cm (5.13)
w = 0.002(450 cm) = 0.9 cm .
¿Cuál es el valor del máximo esfuerzo cortante en un punto al interior de la barra

rectangular durante el ensayo?
Los esfuerzos para la barra cuadrada están dados por
σxxc = λεzzc
σyyc = λεzzc (5.14)
σzzc = (2µ + λ)εzzc
El corte máximo se da entonces

2µ + λ − λ
τmáx = εzzc
2
llamando de nuevo A = (2µ + λ),

µ
τmáx = σzzc
A
para σzzc = 1000 kgf/cm2 ,
µ
τmáx = 1000 kgf/cm2 .
A
Ejemplo 3: Cilindro sometido a torsión
Se tiene un cilindro sometido a torsión con un giro por unidad de longitud 2β. El campo
de desplazamientos para este problema está dado por
u = −2βzy , v = 2βzx , w = 0,
en donde 2β = 2θ/h, representa el giro por unidad de longitud y está dado en rad/m
si los desplazamientos están en m.
Figura 5.3. (Izquierda) Esquema de un cilindro sometido a una torsión en la parte

superior y empotrado en la parte inferior. (Derecha) Configuración deformada para
el cilindro sometido a torsión. Se dibuja una rejilla sobre el sólido por claridad.
Si el cilindro debe ir al interior de una cavidad cilı́ndrica, es decir, el radio de la

cavidad debe ser mayor o igual al del cilindro deformado √ ¿Cuál es el radio interno
mı́nimo de la cavidad si h = 1 m, R = 10 cm y 2θ = 21/10 rad?
Para resolver este problema debemos saber cuál es el máximo radio que el cilindro
va a tener luego de deformarse, para ello necesitamos la región con mayores des-
plazamientos. Para saber en dónde ocurren los mayores desplazamientos podemos
analizar las ecuaciones de desplazamientos para ver que esto pasa en z = h. El
campo de desplazamiento allı́ serı́a (sin considerar el desplazamiento en z que es
nulo):
u = −2θy, v = 2θx .
Alternativamente, podrı́amos saber que el mayor desplazamiento pasa en z = h

analizando la configuración deformada que se presenta en los anexos.
Adicionalmente, sabemos que el nuevo radio del cilindro estarı́a dado por la posi-
ción de los puntos ubicados en la circunferencia x2 + y 2 = R2 . Estos puntos tienen
como ecuaciones paramétricas
x = R cos α ,
y = R sen α ,
en donde α es un parámetro. Tenemos, pues, que la posición de los nuevos puntos
de la circunferencia estarı́an dados por
xdeformado = R cos α − 2Rθ sen α ,
ydeformado = R sen α + 2Rθ cos α .
Usando el teorema de Pitágoras, podemos calcular la magnitud del nuevo radio

2
Rdeformado = x2deformado + ydeformado
2
= [R cos α − 2Rθ sen α]2 + [R sen α + 2Rθ cos α]2

= 4R2 θ2 sen2 α + 4R2 θ2 cos2 α + R2 sen2 α + R2 cos2 α
= 4R2 θ2 [sen2 α + cos2 α] + R2 [sen2 α + cos2 α]
= 4R2 θ2 + R2
= R2 (1 + 4θ2 ) .
que es equivalente a √
Rdeformado = R 1 + 4θ2 ,
o
Rdeformado √
= 1 + 4θ2 .
R
También se pudo haber analizado cómo se movı́a el punto (R, 0), que en la confi-
guración deformada estarı́a en (R, 2θR). Lo que nos permitirı́a llegar a la misma
conclusión, pero de manera más simple.
√ √
Ahora, si 2θ = 21/10, entonces θ = 21/20 y tendrı́amos que
r
Rdeformado 21
= 1+4×
R 400
r
21
= 1+
100
r
121
=
100
11
= ,
10
Lo que nos permite concluir que el radio máximo es 11 cm.
Si el cilindro está hecho de hormigón, material para el cual se tiene que γmáx =
10−4 ¿Cuál es el valor máximo de θ posible antes que el material falle si R = 1 m
y h = 10 m?
El tensor de deformaciones unitarias en este caso es

 
0 0 −βy
[ε] =  0 0 βx  .
−βy βx 0
Para calcular γmáx necesitamos encontrar los valores principales de [ε]. Estos son
la solución del polinomio caracterı́stico
λ3 − λβ 2 (x2 + y 2 ) = 0 ,
o
λ3 − λβ 2 r2 = 0 .
Este polinomio puede factorizarse como
λ[λ2 − β 2 r2 ] = 0 ,
las raı́ces son, entonces,
λ1 = |βr|, λ2 = 0, λ3 = −|βr| ,
y γmáx = |βr| − (−|βr|) = 2|βr|. El valor máximo es cuando r = R y tenemos que
γmáx = 2βR ,
en donde no tuvimos en cuenta el signo de β, pues este sólo expresa si la torsión

se hizo en sentido horario o anti-horario.
Despejando, obtenemos que

γmáx h
θ= ,
2 R
y sustituyendo los valores dados
10−4

10
θ= = 5 × 10−4 .
2 1
Si el cilindro está hecho de acero 1020 con una resistencia al corte τcr = 280 MPa,
un módulo de Young de 186 GPa y un coeficiente de Poisson de 0.29 ¿Cuál es el
valor de par de torsión M que hace fallar al material y a qué ángulo de torsión
θmáx sucede si R = 20 mm y h = 1 m? Use la expresión
Iz
M= τmáx ,
R
para calcular el par de torsión, en donde Iz es el segundo momento de área para
la sección transversal.
πR4
Tenemos que τmáx ≤ τcr , y además Iz = , luego
2
πR3
M= τcr
2
π × (20 × 10−3 m)3 × 280 × 106 Pa
=
2
= 1120π N m
≈ 3518.58 N m .
Procedemos a calcular el ángulo de torsión al que sucede la falla. A partir de las

relaciones de deformación-tensión, tenemos que el esfuerzo está dado por
 
0 0 −y
[σ] = βG  0 0 x  ,
−y x 0
y el esfuerzo cortante máximo estarı́a dado por τmáx = 2βGR, o

θRE
τmáx = ,
(1 + ν)h
y resolviendo esta ecuación para θ, obtenemos
(1 + ν)hτcr
θmáx =
ER
(1 + 0.29) m × 280 × 106 Pa
=
186 × 109 Pa × 20 × 10−3 m
301
= rad
3100
≈ 0.097 rad .
5.6. Ejercicios
1. Para un campo de desplazamientos dado por:

u = a(x2 − 5y 2 )ı̂ + (2axy)̂
a) Determinar el tensor de deformación
b) Obtener las deformaciones principales
c) Obtener los esfuerzos principales
2. En la figura 5.4 se muestra una columna prismática, de sección cuadrada (lado

a), con una tracción axial igual a P . Si el material del que está hecha la pieza
tiene propiedades ν y E, encontrar los campos de esfuerzos y deformaciones.
x3
h
x2
Figura 5.4. Columna con carga axial.
3. En la figura 5.5 se muestra la misma columna del ejercicio anterior, pero en

este caso se quiere estudiar la deformación causado por su propio peso. Si las
componentes del tensor de esfuerzos son:
σ33 = ρg(h − x3 ) ,
σ11 = σ22 = 0 ,
σ12 = σ13 = σ23 = 0 .
a) Calcular los campos de deformaciones si las propiedades del material son ν

y E.
x3
h
x2
Figura 5.5. Columna bajo la acción de su propio peso.
4. En la figura 5.6 se muestra una viga en voladizo de sección unitaria con una carga
distribuida en su extremo de forma parabólica.
P c2 y2

t= 1 − 2 ĵ
2I c
.La solución para esfuerzos es de la forma
P
σxx = − xy ,
I
σyy = 0 ,
P
σxy = − (c2 − y 2 ) .
2I
Y para el campo de desplazamientos es:

ν P y3
2
x 2 c2 P y

1 l
ux = − + − − ,
G E 6I E E G 2I
νP xy 2 P x3 P l2 x P l3
uy = + − + ,
2EI 6EI 2EI 3EI
en donde G = E/(2 + 2ν) es el módulo de cortante e I es el momento de inercia.
a) Determinar las condiciones de frontera del problema.
b) Verificar que las soluciones de esfuerzos y desplazamientos satisfacen las con-

diciones de frontera.
c) Encontrar la región de la viga que está a compresión y la región tensión, y

por tanto el eje neutro. Y determinar la ecuación de desplazamientos para
este eje neutro. Comparar con la teorı́a clásica de vigas.
1 l
c x
c
P
y
Figura 5.6. Viga en voladizo con una carga distribuida en su extremo.
5. En la figura 5.7 se muestra un paralelepı́pedo de dimensiones a, b , c, constituido

por un material homogéneo elástico y lineal se aloja en una cavidad de la misma
forma y dimensiones, cuyas paredes son de un material lo suficientemente rı́gido
para poderlo suponer indeformable. Sobre la abertura de la cavidad de dimensio-
nes a × b y a través de una placa rı́gida de peso y rozamiento despreciables se
aplica, perpendicularmente a ella, una fuerza F que comprime al bloque elástico.
1
Calcular:
a) Las fuerzas laterales ejercidas por las paredes de la cavidad sobre el parale-
lepı́pedo;
b) La variación de altura experimentada por el mismo.
6. En un punto de un suelo que podemos considerar como un sólido elástico lineal

se conoce la deformación volumétrica εV = −2 × 10−3 , la deformación tangencial
εxy = −3 × 10−3 y la deformación εxx = 0 . El suelo está sometido a un estado de
deformación plana en el plano xy. Se pide:
a) Componentes cartesianas del tensor de deformación.
b) Suponiendo que las constante elásticas son E = 50 MPa , ν = 1/4 , obtener

las componentes del tensor de tensiones y sus valores principales. Obtener
asimismo las direcciones en las que las tensiones normales y tangenciales son
máximas o mı́nimas y sus valores.
1
Tomado del ejemplo 3.5 en [10].
F
z
y
c
b
Figura 5.7. Paralelepı́pedo
7. Un sólido se halla sometido a deformación plana, siendo las componentes del

tensor de deformación en un punto:
 
−2 3 0
→(2)
ε =  3 −10 0  × 10−3
0 0 0
Considere que el sólido tiene un comportamiento elástico lineal e isótropo, definido
por módulo elástico de Young E = 10 MPa y coeficiente de Poisson ν = 0.25 . Se
pide:
a) Obtener las deformaciones principales y las direcciones en que se producen;
b) Obtener las componentes del tensor de tensiones de Cauchy;
c) Obtener las máximas y mı́nimas tensiones normales;
d ) Se sabe que el material rompe cuando en algún plano se alcanza una tensión
tangencial que supere 40 kPa. Verificar si se produce la rotura.
8. En la figura 5.8a se muestra una barra, de sección rectangular, a = 40 cm, b = 50

cm y longitud inicial Lc = 500 cm, que es sometida a un ensayo de compresión
confinada (está contenida en un recipiente indeformable). En la figura 5.8b se
muestra una barra, de sección circular, de diámetro d = 50 cm y longitud inicial

Lt = 497 cm sobre la cual se realiza un ensayo de tracción simple. Las constantes
del material, tanto para la barra rectangular como para la barra circular, son
E = 150000 kgf/cm2 y ν = 0.20.
Ambos ensayos se realizan con su extremo fijo en el plano z = 0 y están cargados

con un sistema que no transmite esfuerzos cortantes y que hace que el esfuerzo
axial se distribuya uniformemente sobre la superficie transversal de las barras. Los
esfuerzos que se inducen a las barras varı́an con el tiempo, conforme a lo mostrado
en la figura 5.8c
525 Esfuerzo barra circular
z 350
F z σz (kgf/cm2 )
175
F
00 350 700 1050 1400

Tiempo (s)
0 Esfuerzo barra rectangular
250
Lc
σz (kgf/cm2 )
Lt 500
750
1000
a b d
x y x y 12500 200 400 600
Tiempo (s)
800 1000 1200 1400
(a) Barra com- (b) Barra trac- (c) Curvas de esfuerzo

presión ción
Figura 5.8. Esquema ensayos y curvas de aplicación carga
a) ¿Cuál es el tiempo t en el cual las dos barras tienen la misma longitud?
b) En la barra circular, ¿cuáles son los desplazamientos u, v y w y las deforma-

ciones εxx εyy y εzz en el punto de coordenadas (x, y, z) = (20, 15, 450) cm
para el tiempo t = 300 s?
c) ¿Cuál es el valor del máximo esfuerzo cortante en un punto al interior de la

barra rectangular durante el ensayo?
d ) ¿Cuál es el área mı́nima que alcanza la sección transversal de la barra circular

durante el ensayo?
e) ¿En la barra rectangular, en el tiempo t = 1000 s cuál son las fuerzas Fx y

Fy que le transmite la barra a las paredes del recipiente?
9. En la figura 5.9a se presenta una placa en forma de rombo sometida a la acción

de una carga, en el plano XY , que se encuentra distribuida sobre su perı́metro y
cuyo campo de desplazamientos está dado por:
A A 2Aν
u (x, y, z) = (1 − ν) x, v (x, y, z) = (1 − ν) y, w (z) = − z
E E E
con A una constante positiva, E el módulo de elasticidad del material y ν la
relación de Poisson.
y
y
l l
x
x
l l
(a) Placa (b) Esquema para dibujar las

cargas axiales y tangenciales
Determine la deformación unitaria axial máxima ε y la deformación angular

máxima γ que se presentan en la placa ¿En qué puntos se presentan?.
Dibujar la configuración deformada de la placa en el plano xy ¿Cuáles son

los puntos que experimentan los mı́nimos desplazamientos?.

l
Ilustre la deformación de la partı́cula con coordenadas , 0, 0 en el plano
4
xy. Utilice cuadrados de tamaño diferencial.
Determine y represente gráficamente, sobre la figura (b), la tensión normal

y tangencial a la cara de la placa, contenida en el cuadrante I del plano XY .
Bibliografı́a
[1] Juan H. Cadavid, R.: Mecánica del medio continuo: una iniciación. Fondo Editorial
Universidad EAFIT, 2009.
[2] Bonda, NR y IC Noyan: Effect of the specimen size in predicting the mechanical
properties of PbSn solder alloys. IEEE Transactions on Components, Packaging,
and Manufacturing Technology: Part A, 19(2):208–212, 1996.
[3] Bonda, NR y IC Noyan: Deformation inhomogeneity and representative volume in

Pb/Sn solder alloys. Metallurgical Transactions A, 23(2):479–484, 1992.
[4] Voigt, Woldemar: Theoretische studien über die elasticitätsverhältnisse der krysta-
lle. Königliche Gesellschaft der Wissenschaften zu Göttingen, 1887.
[5] Toupin, Richard A.: Elastic materials with couple-stresses. Archive for Rational
Mechanics and Analysis, 11(1):385–414, 1962.
[6] Mindlin, Raymond David: Micro-structure in linear elasticity. Archive for Rational
Mechanics and Analysis, 16(1):51–78, 1964.
[7] Eringen, A Cemal: Linear theory of micropolar elasticity. Journal of Mathematics

and Mechanics, páginas 909–923, 1966.
[8] Cosserat, Eugene y François Cosserat: Théorie des corps déformables. A. Hermann
et fils,, 1909.
[9] Aero, E. L. y E. V. Kuvshinskii: Fundamental equations of the theory of elastic

media with rotationally interacting particles. Soviet Physics-Solid State, 2(7):1272–
1281, 1961.
187
BIBLIOGRAFÍA 188
[10] Chaves, Eduardo: Problemas resueltos de Mecánica del Medio Continuo. CIMNE,
2014.
[11] Timoshenko, Stephen y J.N. Goodier: Theory of Elasticity. McGraw-Hill, 3a edi-

ción, 1970.
[12] Reddy, Junuthula Narasimha: An introduction to continuum mechanics. Cambrid-

ge university press, 2013.
[13] Mase, G Thomas, Ronald E Smelser y George E Mase: Continuum mechanics for
engineers. CRC press, 2009.
[14] Olivella, Xavier Oliver y Carlos Agelet de Saracı́bar Bosch: Mecánica de medios
continuos para ingenieros, volumen 92. Universidad Politécnica de Catalunya,
2002.
[15] Spiegel, Murray R Murray R y cols.: Teorı́a y problemas de análisis vectorial y una
introducción al análisis tensorial. Mc Gray HIll, 1969.

Notas Medios

Cargado por

Copyright:

Formatos disponibles

Notas Medios

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Notas Medios

Cargado por

Copyright:

Formatos disponibles

Notas de Clase: Mecánica de los Medios Continuos

Grupo de Investigación en Mecánica Aplicada

Grupo de Investigación en Mecánica Aplicada

3.3. Soluciones clásicas para sólidos . . . . . . . . . . . . . . . . . . . . . . 98

4. Análisis de Deformaciones 118

Gran parte de los problemas de Ingenierı́a se encuentran enmarcados en el correcto

Objetivos generales del curso

Desarrollar el concepto matemático y fı́sico de un tensor de orden 2.

Identificar las hipótesis fundamentales y que de manera conjunta con el concepto

Reducir el problema central de la mecánica para un sistema de infinitas partı́culas,

Identificar matemáticamente el problema de la determinación de desplazamientos,

Exposiciones magistrales, Elaboración de Ejemplos, Lecturas asignadas. El curso se di-

1. Presentación del problema (1 semana).

1.1 Introducción y motivación del modelo del continuo.

II. FUNDAMENTOS TEÓRICOS

2.a Análisis de tensiones. Parte 1 (4 semanas).

2.1 Concepto de tensión (Definición, Primer postulado de Cauchy, Definición

2.2 Ecuaciones de transformación y esfuerzos extremos en 2D y 3D. Cı́rculo

2.3 Solución de problemas

2.b Análisis de tensiones. Parte 2 (4 semanas).

2.5 Ecuaciones de equilibrio en un medio continuo.

2.6 Estudio e interpretación de soluciones de tensiones.

3. Análisis de Deformaciones (5 semanas).

3.1 Transformaciones lineales.

3.2 Definición de concepto de deformación.

3.3 Tensor de deformaciones.

3.4 Teorı́a de la Elasticidad. Ley de Hooke.

4. Proyecto docente (2 semanas).

4.1 Estudio de solución numérica.

Al final de esta sección el estudiante debe estar en la capacidad de:

• Entender el paso de variables discretas a variables continuas en términos de den-

• Identificar la diferencia entre el concepto de partı́cula en el modelo de la mecánica

• Identificar las funciones fuente o causa y efecto o respuesta en el modelo del

• Identificar los lı́mites de aplicabilidad del modelo del continuo.

1.2. Concepto de densidad generalizada

A una escala de tamaño lo suficientemente pequeña la materia goza de una composición

Si para esta colección se desea calcular la distribución en el tiempo de alguna función

Por ejemplo, en el problema mecánico la ecuación (1.1) corresponderı́a a la segunda ley

Por tanto, si se dispone de un postulado fı́sico fundamental como el dado en la ecua-

Un tratamiento aproximado se hace posible gracias a la Hipótesis Básica de Continuidad

lı́m ∆V = Punto . (1.2)

La generalización de este concepto al sistema de infinitas partı́culas se reduce a despre-

Podemos en este punto precisar el objetivo fundamental de un modelo de los medios

Un modelo de los medios continuos pretende estudiar la distribución espacio-

1.3. Conexión entre el continuo fı́sico y el continuo

Retomando el tratamiento del problema de determinación de una función respuesta

Se renuncia al concepto de partı́cula.

Se hace disponible la potente herramienta del continuo matemático o cálculo

La hipótesis de continuidad genera entonces un modelo que producirá resultados co-

El hecho de renunciar al concepto de partı́cula y de apelar al continuo matemático

Un modelo de los medios continuos pretende estudiar la variación espacio-

En el caso particular del problema mecánico el interés es describir la respuesta del

nes continuas representativas de la cantidad de movimiento en respuesta a funciones

1.4. Caracterı́sticas generales de un modelo del me-

Para apoyar el análisis considérese un medio infinito ocupando un volumen B(t) y

¿Cómo describir la interacción o intercambio de cualquier variable fı́sica del tipo

Figura 1.1. Medio continuo interactuando con un espacio infinito.

lı́m ∆V = P unto material (1.3)

A la luz de la continuidad matemática el punto material definido en la ecua-