Introducci

on a la mec
anica
de los s
olidos deformables

Ignacio Romero Olleros

E.T.S. Ingenieros Industriales
cnica de Madrid
Universidad Polite

4 de Febrero, 2014

Indice

Captulo 1. Algebra
y c
alculo tensorial . . . . . . . . . . . . . . . . . . . . . . . . . . .

1. Vectores y campos vectoriales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

1.1. Componentes y cambio de base . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

1.2. Operaciones algebraicas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

1.3. Calculo vectorial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

2. Tensores y campos tensoriales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

2.1. Componentes y cambio de base . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

2.2. Operaciones algebraicas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

2.3. Tensores con propiedades especiales . . . . . . . . . . . . . . . . . . . . . . . . .

2.4. Descomposiciones de tensores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

2.5. Autovectores y autovalores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

10

2.6. Calculo tensorial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

11

3. Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

13

Captulo 2. Estudio del equilibrio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

15

1. El modelo del s
olido deformable . . . . . . . . . . . . . . . . . . . . . . . . . . . .

15

2. Fuerzas que act

uan sobre los s
olidos deformables . . . . . . . . . . . .

16

2.1. Fuerzas volumetricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

17

2.2. Fuerzas de superficie o de contacto . . . . . . . . . . . . . . . . . . . . . . . . . .

17

3. Fuerzas internas en un cuerpo deformable . . . . . . . . . . . . . . . . . .

18

4. El tensor de tensiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

21

4.1. Interpretacion fsica de las componentes del tensor de

tensiones . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

24

5. Ecuaciones de equilibrio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

25

5.1. Principio fundamental . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

25

5.2. Equilibrio de fuerzas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

27

5.3. Equilibrio de momentos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

27

6. Tensiones principales y direcciones principales de

tensi
on . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

28

7. Representaci
on gr
afica de un tensor de tensiones . . . . . . . . . . . .

29

8. Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

35

Captulo 3. An
alisis de la deformaci
on de los cuerpos . . . . . . . . . . . . . .

45

1. Introducci
on . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

45

2. Cinem
atica de un cuerpo deformable . . . . . . . . . . . . . . . . . . . . . . .

46

3. C
alculo de deformaciones. El tensor de deformaci
on
infinitesimal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

49

3.1. El tensor de deformaciones infinitesimales . . . . . . . . . . . . . . . . . . . .

49

3.2. Calculo de deformaciones longitudinales . . . . . . . . . . . . . . . . . . . . .

50

3.3. Calculo de deformaciones angulares . . . . . . . . . . . . . . . . . . . . . . . . .

52

3.4. Interpretacion geometrica de las componentes del tensor de

deformacion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

54

4. La deformaci
on volum
etrica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

56

5. Deformaciones principales y direcciones principales de

deformaci
on . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

58

6. Galgas extensom
etricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

59

7. El diagrama de Mohr de deformaciones . . . . . . . . . . . . . . . . . . . . .

60

8. Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

61

ii

Captulo 4. Elasticidad y termoelasticidad lineal . . . . . . . . . . . . . . . . . .

67

1. El concepto de elasticidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

68

2. El principio de superposici
on . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

68

3. Las constantes el
asticas de un material is
otropo . . . . . . . . . . . . .
3.1. El modulo de Young . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3.2. El coeficiente de Poisson . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

69
70
70

4. La ley de Hooke generalizada . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

71

5. Las ecuaciones de Lam

e. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

73

6. Deformaciones y tensiones proporcionales . . . . . . . . . . . . . . . . . . .

74

7. Termoelasticidad lineal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

75

8. Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

78

Captulo 5. Estudio del problema el

astico completo . . . . . . . . . . . . . . . .

81

1. Enunciado completo del problema . . . . . . . . . . . . . . . . . . . . . . . . . .

82

2. El principio de los trabajos virtuales . . . . . . . . . . . . . . . . . . . . . . . .

83

3. Las ecuaciones de Navier . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

84

4. El principio de Saint Venant . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

85

5. Estados planos de tensi

on y deformaci
on . . . . . . . . . . . . . . . . . . . .
5.1. Estados de tension plana . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.2. Estados de deformacion plana . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
5.3. El diagrama de Mohr en estados planos . . . . . . . . . . . . . . . . . . . . . .

85
86
86
87

6. Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

90

Captulo 6. La energa el
astica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

93

1. El trabajo de las fuerzas exteriores sobre un cuerpo

deformable . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

93

2. La energa el
astica de deformaci
on . . . . . . . . . . . . . . . . . . . . . . . . . .

95

iii

3. El teorema de reciprocidad de Maxwell-Betti . . . . . . . . . . . . . . .

97

4. Introducci
on a los modelos estructurales . . . . . . . . . . . . . . . . . . . .

98

5. Los coeficientes de influencia . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

98

6. Los teoremas de Castigliano . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 101

7. Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103

Captulo 7. Estudio de la finalizaci

on del comportamiento
el
astico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 109
1. Criterios de fallo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 110
2. Criterios de fluencia para materiales d
uctiles . . . . . . . . . . . . . . . . 111
2.1. El criterio de Tresca . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
2.2. El criterio de von Mises . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
3. Criterios de rotura para materiales fr
agiles . . . . . . . . . . . . . . . . . 113
3.1. El criterio de Rankine . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
3.2. El criterio de Mohr-Coulomb . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 114
3.3. El criterio de Drucker-Prager . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 115
4. Problemas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 116
Bibliografa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 119

iv

Captulo 1

Algebra
y c
alculo tensorial

Cada disciplina de la ciencia y la ingeniera tiene su lenguaje, que ha ido

evolucionando con los a
nos y que permite su descripcion de la forma mas clara y
sencilla posible. La elasticidad es una teora de campos y, como tal, se expresa mas
clara y sencillamente en el lenguaje de los vectores y los tensores. Se podra decir
que el lenguaje natural de la teora de la elasticidad, aquel en el que los conceptos
aparecen mas claramente representados, es el de los vectores y los tensores, y por ello
parece recomendable su uso. En estas paginas se recogen los conceptos mas basicos,
que son los que se emplearan en el desarrollo de la asignatura. Un desarrollo mas
completo de los conceptos del algebra y calculo tensorial se pueden encontrar, entre
muchos otros, en [4, 5, 6, 8].

1. Vectores y campos vectoriales

La definicion completa de un vector y un campo vectorial se puede consultar en
cualquier libro de algebra. En lo que sigue, llamaremos vector simplemente a un
elemento cualquiera de Rd , siendo d igual a 3 en estas notas, aunque la gran parte
de los conceptos que se presentan son validos tambien para otras dimensiones. Un
campo vectorial definido en R3 es una funcion que para todo punto en
devuelve un vector.

Notaci
on: Para diferenciar los vectores de los escalares se emplean en la literatura
distintas notaciones. As, dependiendo del libro u autor que se consulte, un mismo
vector se puede ver escrito como u,
u, ~u, u, . . . Entre todas estas utilizaremos la
primera.

1.1. Componentes y cambio de base

Sea B = {e1 , e2 , e3 } una base cartesiana de R3 . Cualquier vector v R3 se
puede expresar de la forma
v = v 1 e1 + v 2 e2 + v 3 e3 =

3
X

v i ei ,

(1.1)

i=1

y v1 , v2 , v3 se llaman las componentes de v en la base B. Las componentes de un

vector cambian seg
un la base a la cual se refieran, por lo tanto no se debe confundir
el vector mismo con su representacion.
Para expresar que una terna v1 , v2 , v3 son las componentes del vector v en la
base B escribiremos:

v1
{v}B = v2
.
(1.2)

v3 B
A menudo, cuando no hay posibilidad de confusion porque solo se ha definido una
base se emplea la notacion


v1
v1
{v} = v2
, o simplemente v = v2
.
(1.3)


v3
v3

La relacion entre las componentes de un vector, referidas a dos bases distintas

se obtiene de la siguiente manera. Sea B la base anteriormente definida y B 0 una
nueva base cartesiana formada por los vectores ortormales {e01 , e02 , e03 }. Un vector
ei cualquiera de la base B se puede expresar como suma de vectores de la base B 0
de la forma:
ei = ai1 e01 + ai2 e02 + ai3 e03 ,

siendo

aij = ei e0j .

(1.4)

Otro vector cualquiera v se puede escribir indistintamente como combinacion lineal

de los elementos de B o de los de B 0 :
v=

3
X

v i ei =

i=1

3
X

vj0 e0j .

j=1

Sustituyendo la expresion (1.4) e identificando las componentes se obtiene

vj0

3
X

aij vi .

(1.5)

i=1

Esta u
ltima relacion se puede expresar
0

a11
v1
v20
= a12
0 0
a13
v3 B

matricialmente como

a21 a31 v1
,
a22 a32 v2

v3 B
a23 a33

(1.6)

o de forma compacta
{v}B 0 = [A]T {v}B .

(1.7)

Si las bases B y B 0 son ortonormales, la matriz de cambio de base [A], es una matriz
ortogonal, es decir, que verifica [A]1 = [A]T . Mas a
un, si las dos bases tienen la
misma orientacion ([e1 e2 e3 ] = [e01 e02 e03 ], ver debajo el significado de la operacion
[]), entonces el determinante de [A] es igual a 1 por tanto es una rotacion.
Es habitual referirse a los vectores de la base cartesiana de R3 como {i, j, k} y
las componentes de un vector v en dicha base como vx , vy , vz .

1.2. Operaciones algebraicas

Los vectores de Rd poseen las operaciones vectoriales basicas de suma y
multiplicacion por un escalar. Para realizar operaciones vectoriales nos vemos

obligados a menudo a emplear las componentes de un vector, pero es importante

recalcar que el resultado es independiente de la base escogida, siempre que todos los
vectores que intervengan se expresen en la misma base. Por ejemplo, para calcular el
vector c = a + b, utilizamos las componentes de todos ellos en la base B y podemos
emplear la expresion:



c1
a1
b1
c
= a2
+ b2
.
(1.8)
2


c3 B
a3 B
b3 B
Ademas, en el espacio eucldeo se define el producto escalar de dos vectores con
la expresion
a b = a1 b1 + a2 b2 + a3 b3 .

(1.9)

El producto escalar, como el resto de operaciones de las que tratamos, es una

operacion intrnsica que no depende de la base escogida. La norma de un vector
se indicara como |a| y se define de la siguiente forma

|a| = a a .
(1.10)
Cualquier vector no nulo se puede normalizar, multiplicandose por el inverso de su
norma, y obteniendose un vector unitario. Dado un vector cualquiera a y otro
vector cualquiera unitario u, se definen la proyeccion de a sobre u y la proyeccion
de a sobre el plano normal a u como
k

a
u = a au .

au = (a u)u ,

(1.11)

Esta descomposicion es u
nica y se puede escribir a = au + a
u.
El producto vectorial de dos vectores se indica como a b y se define el
producto mixto a, b, c como
[a b c] = a b c .

(1.12)

1.3. C
alculo vectorial
En Teora de Campos se estudian los principales operadores diferenciales que
act
uan sobre los campos escalares y vectoriales. Estos son el gradiente, la divergencia
y el rotacional. Para definirlos, consideramos en esta seccion una base cartesiana

{i, j, k} y coordenadas x, y, z. Ademas usamos la notacion

,x =

,
x

,y =

,
y

,z =

,
z

(1.13)

para cualquier campo = (x, y, z).

Sean y v respectivamente un campo escalar y uno vectorial definidos sobre
R3 y diferenciables. Las operaciones diferenciales mencionadas se definen, en
componentes:

,x
{grad (x, y, z)} = ,y
,

,z
div v(x, y, z) = vx,x + vy,y + vz,z ,
(1.14)


 i

j
k 

rot v(x, y, z) =  x y z 
 vx vy
vz 
Existen numerosas relaciones entre los operadores diferenciales y teoremas
integrales que los emplean. En este curso utilizaremos dos u
nicamente: el teorema
de la divergencia y la formula de la integral por partes.

Teorema 1.1: Teorema de Gauss o de la divergencia Sea un subconjunto

de R3 y su contorno. Si v es un campo vectorial definido en se verifica:
Z
Z
div v d =
v n d ,
(1.15)

siendo n la normal saliente a .

Para el mismo dominio, si es un campo escalar, entonces
Z
Z
Z
div v d =
v n d
v grad d .

(1.16)

2. Tensores y campos tensoriales

Los tensores son objetos son objetos del algebra tan u
tiles como los vectores, y
su uso en varias ramas de la mecanica es muy habitual. Aunque existen tensores de
cualquier orden (entero), en esta asignatura emplearemos el termino tensor para

referirnos a los de segundo orden, los mas habituales en mecanica. Por ejemplo, el
tensor de inercia que se emplea en Mecanica Clasica es un tensor de segundo orden.
Los tensores se estudian a menudo en algebra bajo el nombre de
homomorfismos y son simplemente aplicaciones lineales de Rd en Rd , es decir,
funciones lineales que transforman un vector en otro. Para cualquier vector a Rd ,
un tensor T es una operacion lineal tal que T (a) es otro vector. Por sencillez, los
parentesis se eliminan y se escribe simplemente b = T a.
Un campo tensorial no es mas que una funcion que para cada punto de un
dominio devuelve un tensor. Volviendo al ejemplo de la Mecanica Clasica, el tensor
de inercia es un campo tensorial que depende del punto donde se calcule. Ademas
transforma vectores en vectores. Si el punto de evaluacion es el centro de gravedad
o un punto fijo, este tensor transforma la velocidad angular en el momento cinetico
respecto al punto.

Notaci
on: Igual que en el caso de los vectores, existe una notacion especial que
permite distinguir los tensores de segundo orden del resto de objetos (escalares,
vectores, . . . ). Tambien esta notacion depende del autor o del libro que se consulte
~~
A,
y un mismo tensor se puede escribir como A, A,
A, . . . En estas notas se empleara
la primera de ellas y se evitara la confusion entre vectores y tensores de segundo
orden empleando siempre que no se indique lo contrario letras min
usculas en el
primer caso y may
usculas en el segundo.

2.1. Componentes y cambio de base

Recordamos que un tensor es simplemente una operacion que transforma
vectores en vectores, y que es lineal. Pues bien, en particular se pueden usar
tensores y operarlos sobre los vectores de una base B. Con ello se pueden definir las
componentes de un tensor T como los nueve escalares
Tij = ei (T ej ) ,

i = 1, 2, 3

Las componentes de un tensor referidas a

matriz, y se escribe

T11 T12

[T ]B = T21 T22
T31 T32

j = 1, 2, 3 .

(1.17)

una base B se muestran en forma de

T13
T23 ,
T33 B

(1.18)

de forma analoga a la expresion en un vector columna de un vector (1.2). Como en

el caso de los vectores, la matriz de un tensor en una base cualquiera no debe de
confundirse con el tensor propiamente dicho.
La propiedad de linealidad de los tensores implica que las componentes del vector
b que resulta de la aplicacion de un tensor T sobre un vector a se pueden obtener
multiplicando la matriz [T ]B y el vector columna {a}B . Es decir, si b = T a,
entonces

T11 T12 T13
a1
b1
a
.
{b}B = [T ]B {a}B , o mas explcitamente
b2
= T21 T22 T23
2

b3 B
T31 T32 T33 B a3 B
(1.19)
Observando la definicion de las componentes de un tensor deducimos que estas
dependen de la base en la que se exprese el tensor. Para hallar la relacion entre
componentes de un mismo tensor en dos bases cartesianas distintas B y B 0 escribimos
la relacion b = T a en componentes de las dos bases.
{b}B = [T ]B {a}B ,

y {b}B 0 = [T ]B 0 {a}B 0 .

(1.20)

La expresion (1.7) relaciona las componentes de los vectores a y b en las dos bases
as que la segunda ecuacion de (1.20) se puede escribir como
[A]T {b}B = [T ]B 0 [A]T {a}B .

(1.21)

Despejando {b}B y comparando el resultado con la primera ecuacion de (1.20)

se deduce que la expresion que relaciona las componentes de T en las dos bases
consideradas es
[T ]B = [A][T ]B 0 [A]T .
(1.22)

2.2. Operaciones algebraicas

Los tensores poseen las operaciones de suma, multiplicacion y multiplicacion
por un escalar. La expresion matricial del resultado de todas estas operaciones es
la correspondiente operacion matricial operada sobre las matrices de componentes
de los tensores. Insistimos, como en el caso de los vectores, que el resultado es
independiente de la base escogida.
La traza de un tensor es la suma de los elementos de la diagonal de su matriz
de compomentes:
tr[T ] = T11 + T22 + T33 .
(1.23)
7

La traza de un tensor no depende tampoco de la base en la que se exprese su matriz

de componentes y se dice que es por tanto un invariante del tensor. Esto se puede
comprobar calculando la traza del tensor en la ecuacion (1.22). La operacion traza
es lineal as que, dado un escalar y dos tensores T , S,
tr[T ] = tr[T ] ,

tr[T + S] = tr[T ] + tr[S] .

(1.24)

El producto escalar entre tensores de orden dos se escribe con el smbolo :

y se define como la operacion que a toda pareja de tensores T , S asocia el escalar
T : S definido por
T : S = tr[S T T ] .
(1.25)
En componentes, la operacion de la doble contraccion, como tambien se conoce a
este producto escalar, es simplemente
T :S=

3 X
3
X

Tij Sij .

(1.26)

i=1 j=1

Como esta operacion define un producto escalar, tambien se puede definir una
norma asociada de tensores:

kT k = T : T .
(1.27)
El determinante de un tensor T es el escalar det(T ) que verifica
[T a T b T c] = det(T )[a b c] .

(1.28)

Ademas, se puede demostrar, que el determinante se puede obtener calculando

el determinante de la matriz de componentes del tensor, en cualquier base. El
determinante es, por tanto, otro invariante del tensor.

2.3. Tensores con propiedades especiales

Dependiendo de sus propiedades, los tensores se clasifican empleando unos
calificativos identicos a los de la clasificacion de las matrices. En primer lugar,
el tensor identidad I es el u
nico que verifica Ia = a para todo vector a. La
matriz de componentes de I, en cualquier base, es la matriz identidad. El tensor
nulo es el u
nico tensor tal que T + 0 = T , para cualquier tensor T .

Dado un tensor cualquiera T definimos el tensor traspuesto T T como aquel

que verifica
(T T a) b = a (T b) ,
(1.29)
para cualquier pareja de vectores a, b. Usando esta definicion se dice que un tensor
es simetrico si es igual a su traspuesto, y antisimetrico (o hemisimetrico) si es el
opuesto de su traspuesto. Se comprueba inmediatamente que la matriz asociada a
un tensor simetrico es simetrica, en cualquier base, y la matriz asociada a un tensor
antisimetrico es a su vez antisimetrica, tambien en cualquier base.
Los tensores antisimetricos tienen una propiedad que emplearemos mas adelante
y es que el efecto de aplicar un tensor antisimetrico W sobre un vector cualquiera
a es el mismo que el de multiplicar vectorialmente un vector w, llamado el vector
axial de W , sobre a. Es decir, que para todo a,
Wa = w a .

(1.30)

Ademas esta relacion es recproca, y por ello multiplicar vectorialmente un vector

w por otro vector cualquiera a es equivalente a multiplicar un tensor antisimetrico
W , que es u
nico, y que se llama el tensor antisimetrico asociado al vector w.
Un tensor es desviador si tiene traza nula y un tensor es esferico si es de la
forma T = pI, siendo p un escalar.

2.4. Descomposiciones de tensores

Todo tensor T se puede descomponer de forma unvoca en una parte simetrica
T s y otra antisimetrica T a de forma que T = T S + T a . Se comprueba facilmente
que cada una de estas partes son:
T s = 12 (T + T T ) ,

T a = 12 (T T T ) .

(1.31)

Ademas, todo tensor T se puede descomponer de forma unvoca en una parte

esferica y otra desviadora. La parte esferica, que denominaremos T vol tiene la misma
traza que T y la parte desviadora T des , no tiene traza. Ambas se calculan as:
T vol =

tr[T ]
I ,
3

T des = T T vol .

(1.32)

2.5. Autovectores y autovalores

Dado un tensor T cualquiera, se dice que el vector v es un autovector y su
autovalor asociado si v es unitario y
T v = v .

(1.33)

Para calcular los autovalores buscamos las soluciones no triviales de la ecuacion

(1.33) y para ello hay que resolver la ecuacion
det(T I) = 0 .

(1.34)

Esta ecuacion es un polinomio de tercer grado que tiene por expresion

3 + I1 (T )2 I2 (T ) + I3 (T ) = 0 .

(1.35)

Las funciones I1 , I2 , I3 son los llamados invariantes principales del tensor T , porque
no dependen de la base, y su expresion explcita es
I1 (T ) = tr[T ] ,

I2 (T ) = 21 (tr[T ]2 tr[T 2 ]) ,

I3 (T ) = det(T ) .

(1.36)

Como en algebra de matrices, una vez calculados los tres autovalores 1 , 2 y

3 , se calculan sus autovectores asociados buscando las bases de los espacios nulos
de los tensores
T I ,
(1.37)
que no son vacos por definicion. Cuando el tensor es simetrico, el siguiente teorema
espectral garantiza que los autovalores y autovectores cumplen una propiedades
especiales que se emplearan muy a menudo en la mecanica de solidos deformables.
Por su importancia incluimos una demostracion del teorema.

Teorema 2.2: Los tres autovalores de un tensor simetrico S son reales y sus tres
autovectores asociados forman una base ortonormal, llamada la base principal del
tensor, que denominamos B . En esta base, la expresi
on matricial del tensor es:

1 0
0
[S]B = 0 2 0
.
(1.38)
0
0 3 B

10

 n: Demostramos primero que los tres autovalores son reales. Como el

Demostracio

polinomio caracterstico de S es de tercer orden existen tres autovalores 1 , 2 , 3

que en principio pueden ser complejos. Si v es el autovector asociado a un autovalor
es el autovector conjugado entonces
de los tres y v
Sv = v
v = |v|2 .
v

(1.39)

Conjugando ambos lados de la ecuacion anterior, se obtiene

2.
v S
v = |v|

(1.40)

es decir que
Igualando las identidades de (1.39) y (1.40) concluimos que = ,
es real.
La demostracion de la segunda parte es inmediata si los autovales son distintos,
pero consideramos el caso mas general. Como antes, sean (1 , 2 , 3 ) los tres
autovalores de S (ordenados de cualquier manera) y (v1 , v2 , v3 ) sus autovectores
correspondientes. Si w es un vector ortogonal a v1 , entonces Sw es tambien
ortogonal a v1 , porque v1 Sw = Sv1 w = 1 v1 w = 0. As pues S, cuando
se restringe al subespacio de vectores ortogonales a v1 es tambien un tensor de
ese conjunto a s mismo. Por lo tanto tendra dos autovalores y autovectores, que
forzosamente deberan ser ortogonales a v1 . Tomando uno cualquiera que llamamos
2 y v2 al autovector, repetimos el mismo argumento para el subespacio de vectores
ortogonales a v1 y v2 para concluir que los tres autovectores son ortonormales.
En la base principal tenemos


0
1

,
, {v2 }B = 1
{v1 }B = 0


0 B
0 B

,
{v3 }B = 0

1 B

(1.41)

por lo que la expresion matricial de S ha de ser como se indica en (1.38).

Un tensor simetrico con dos autovalores iguales se llama cilndrico, y cuando
los tres son iguales, esf
erico. En este u
ltimo caso el tensor ha de ser proporcional
al tensor identidad.

2.6. C
alculo tensorial
Como en el caso del calculo vectorial resumimos u
nicamente los elementos mas
basicos del calculo tensorial, aquellos que emplearemos en la asignatura. En primer
lugar hay que indicar que existen expresiones para el gradiente y la divergencia

11

de tensores de cualquier orden (no solo de segundo orden). Puesto que no son
necesarios no los definiremos aqu. Tan solo seran necesarios el gradiente de un
campo vectorial, que es un campo tensorial, y la divergencia de un campo tensorial,
que es un campo vectorial.
Para definir el gradiente de un campo vectorial de forma intrnsica (sin
coordenadas) se usa la misma idea que el caso de un campo escalar. En esta situacion
se define el gradiente de un campo escalar como aquel u
nico vector grad que
verifica que, para cualquier pareja de puntos Po y P
(P ) = (Po ) + grad (Po ) r + O(|r|2 ) ,

(1.42)

siendo r el vector que va del punto Po al punto P . En el caso de un campo vectorial

v su gradiente es el u
nico tensor grad v tal que
v(P ) = v(Po ) + (grad v(Po ))r + O(|r|2 ) .
En una base cartesiana las componentes del tensor grad v son

vx,x vx,y vx,z

[grad v(x, y, z)] = vy,x vy,y vy,z .
vz,x vz,y vz,z

(1.43)

(1.44)

El operador divergencia de un campo tensorial T se define de forma intrnsica

como aquel u
nico campo vectorial div T que verifica
(div T ) a = div (T T a) ,

(1.45)

para todo vector a. En coordenadas cartesianas, las componentes del vector div T
son

Txx,x + Txy,y + T xz, z

{div T (x, y, z)} = Tyx,x + Tyy,y + T yz, z
.
(1.46)

Tzx,x + Tzy,y + T zz, z

El teorema de la divergencia para campos tensoriales es practicamente
identido a (1.15), pues establece que para toda region R3 de contorno y
un campo tensorial T definido en ella, la integral de la divergencia de T es igual al
flujo saliente de T :
Z
Z
div T d =
T n d ,
(1.47)

12

siendo n la normal saliente a . Para el mismo dominio, si v es un campo vectorial,

entonces la formula de la integraci
on por partes es
Z
Z
Z
div T v d = (T n) v d
T : grad v d .
(1.48)

3. Problemas
1) Demuestra que, para cualquier tensor T , su traza verifica tr[T ] = I : T . Calcula
la traza del tensor identidad.
2) Demuestra el resultado (1.16) utilizando el teorema de la divergencia sobre el
campo vectorial v.
3) Comprueba que el tensor antisimetrico asociado a un vector de componentes
v = {vx , vy , vz }T es el que tiene por matriz de componentes

0
vz vy
[V ] = vz
0
vx .
vy vx
0
4) Si S es un tensor simetrico y W antisimetrico, comprueba que S : W = 0. si
V es un tensor esferico y E uno desviador, comprueba tambien que V : E = 0.
5) Demuestra que el producto escalar de vectores no depende de la base escogida.
6) Si u es un campo vectorial, demuestra que el vector axial de la parte
antisimetrica del tensor grad (u) coincide con 12 rot (u).
7) El tensor T en la base B = {e1 , e2 , e3 } tiene la siguiente expresion matricial

2 3 0
[T ] = 3 4 0 .
0 0 5
Encontrar
la expresion matricial
de T en la base B 0 = {e01 , e02 , e03 } si e01 =

3
3
1
1
0
0
2 e1 + 2 e2 , e2 = 2 e1 + 2 e2 y e3 = e3 . Comprobar que la traza y el
determinante no dependen de la representacion matricial.

13

14

Captulo 2

Estudio del equilibrio

1. El modelo del s
olido deformable
La materia no es contiua. Si empleamos un microscopio de suficiente resolucion
podremos apreciar como esta se compone de multitud de atomos separados entre
s, los cuales a su vez estan formados por un n
ucleo diminuto y nubes de electrones
lejanos a estos. Esta observacion es valida para cualquier tipo de cuerpo: solido,
lquido o gaseoso.
En cualquier rama de la ciencia, e ingeniera en particular, se formulan
modelos matematicos de la realidad para poder explicar su comportamiento y
poder predecir su comportamiento futuro. Todos los modelos son inexactos, pues
asumen simplificaciones de la materia para que las ecuaciones resultantes puedan
ser manejables y se puedan resolver al menos en algunos casos. Hay en cualquier

15

caso una cierta jerarqua de modelos de la realidad fsica, desde los mas sencillos e
inexactos, hasta los mas complejos y precisos.
En mecanica de solidos, que es la disciplina que nos concierne, existe el modelo
de partcula que Galileo y Newton, entre otros, introdujeron. Seg
un este modelo,
la dinamica de solidos puede estudiarse considerando que estos son puntos dotados
de masa. Un modelo de complejidad mayor es el de solido rgido, que incorpora
detalles sobre la distribucion de la masa y la orientacion de los cuerpos.
Los dos modelos indicados no describen ni la deformabilidad de los cuerpos,
ni la posibilidad de rotura/fallo, ni las diferencias entre distintos materiales, ni los
efectos de la temperatura sobre los cuerpos... Para incorporar todos estos aspectos
se formula un modelo mas complejo, llamado el modelo de solido deformable, que
sigue siendo imperfecto e inexacto, pero cuya precision a la hora de reproducir lo
que ocurre con los solidos reales es mucho mayor que la de la partcula o el solido
rgido.
En este curso describiremos el modelo de solido deformable para formados por
materiales:
continuos, por lo tanto ignorando la estructura atomica de la materia;
homogeneos, es decir indistinguibles punto a punto;
el
asticos, o sea cuya deformacion desaparece de forma instantanea cuando las
cargas se retiran;
is
otropos, que indica que las propiedades son las mismas en todas las direcciones
del espacio.
A
un con estas limitaciones, el modelo matematico resultante es muy complejo,
as que incorporamos una restriccion mas, que no tiene que ver con el material sino
con el tipo de problema que vamos a estudiar y es que solo consideramos problemas
en los que las deformaciones y sus gradientes son muy peque
nas.
Con estas restricciones, estudiaremos el comportamiento de solidos, entendidos
como subconjuntos R3 , con contorno formados por un conjunto continuo de
puntos que debido a la accion de fuerzas exteriores y/o temperatura pueden adoptar
formas distintas a la original.

2. Fuerzas que act

uan sobre los s
olidos deformables
Como la mecanica trata de las fuerzas y su efecto sobre los cuerpos, el primer
paso para describir en que consiste el solido deformable consiste en delimitar que
fuerzas vamos a considerar y cuales no.

16

En el caso de la partcula, las u

nicas fuerzas que se admiten son las fuerzas
puntuales. En el modelo del solido rgido, ademas de las primeras, se adminten
pares de fuerzas. En el modelo de solido deformable no se admite ninguna de las
dos anteriores y sin embargo se permiten dos nuevos tipos de fuerzas llamadas
fuerzas volumetricas y fuerzas de superficie.

2.1. Fuerzas volum

etricas
Las fuerzas volumetricas son fuerzas que act
uan sobre cada diferencial de
volumen del cuerpo, o equivalentemente, sobre cada diferencial de masa. El ejemplo
clasico es el de la fuerza de la gravedad, que act
ua sobre cada elemento diferencial
de volumen tirando de el hacia abajo. Matematicamente, las fuerzas volumetricas
se describen con un campo vectorial fv : R3 de forma que sobre el diferencial
de volumen en el punto P act
ua una fuerza diferencial fv dv. La resultante,
por tanto, de todas las fuerzas volumetricas que act
uan sobre un cuerpo es
Z
Rv =
fv (P ) dv ,
(2.1)

o en componentes en una base cartesiana B = {i, j, z},

Rvx Z fvx (P )
f (P ) dv .
Rvy =

vy

fvz (P )
Rvz

(2.2)

2.2. Fuerzas de superficie o de contacto

Las fuerzas de superficie, tambien llamadas fuerzas de contacto, son fuerzas
aplicadas sobre el cuerpo a traves de su contorno . Matematicamente se expresan
como un campo vectorial fs : R3 definido sobre el contorno de fuerzas por
unidad de superficie. Sobre un diferencial de area sobre el punto P act
ua
una fuerza total de valor fs dA y, por tanto, la resultante de todas las fuerzas de
superficie actuando sobre un cuerpo es:
Z
Rs =
fs (P ) dA ,
(2.3)

17

o en componentes en una base cartesiana

Rsx Z fsx (P )
R
=
fsy (P ) dA .
sy
f (P )
Rsz
sz

(2.4)

En mecanica, se denomina tensi

on a la fuerza aplicada sobre la unidad de area. A
diferencia de la presi
on, la tension tiene direccion y sentido.

Observaciones 2.1:
i. Los solidos deformables no admiten fuerzas ni pares puntuales.
ii. Las fuerzas volumetricas tienen dimensiones de F/L3 y las de superfice, de F/L2 .
iii. Las fuerzas de superficie se pueden descomponer en su componente normal y
tangencial a la superficie del cuerpo. Si consideramos un punto P , y la
normal a n en dicho punto, podemos calcular
k

fs n = (fs n)n ,

fs
n = fs fs n .

(2.5)

Es com
un estudiar solidos deformables sujetos en una parte de su contorno que
denominaremos u , de forma que = u t con u t = . En el contorno u es
solido deformable tiene sus desplazamientos impedidos y para ello la sustentacion
ejerce unas fuerzas de superficie de valor desconocido a priori que se encargan de
satisfacer dicha restriccion. Por el contrario, o bien t es una superficie libre o bien
hay fuerzas de superficie conocidas, de tal manera que los desplazamientos de sus
puntos son desconocidos cuando se plantea el problema (ver Figura 2.1).
La resultante de la fuerzas de superficie sobre u se llama la reacci
on sobre el
cuerpo y no se conoce hasta que se resuelve todo el problema de contorno.

3. Fuerzas internas en un cuerpo deformable

El concepto de fuerza interna es central para el estudio de cuerpos deformables
y es nuevo, en el sentido de que no existe en los modelos de partculas con masa o
en el de cuerpos rgidos.
Cuando se estudia un cuerpo deformable sometido a fuerzas externas se deduce
que, aunque no se puedan medir, deben de existir fuerzas en el interior del mismo.
Estas no se pueden medir porque para ello habra que partir el cuerpo, creando una

18

fs

fv

u
Figura 2.1:

El modelo de cuerpo deformable

nueva superficie externa y por tanto dejaran de ser fuerzas internas. Pero, sin duda
deben de existir para mantener la cohesion entre sus partculas y para transmitir
las fuerzas aplicadas desde el exterior, en la superficie o en el interior.
Para comprender este nuevo concepto, consideramos un cuerpo formado por una
parte roja y otra verde unidas por una superficie suave S. Cuando sometemos a
todo el cuerpo a fuerzas externas (volumetricas o de superficie) tambien aparecen
fuerzas entre las dos partes diferenciadas por sus colores. Si quitasemos la parte
verde, para que la parte roja no lo notara deberamos aplicar sobre la superficie
S algunas fuerzas que sustituyeran el efecto de la parte roja sobre la primera, y
viceversa. Estas fuerzas no se controlan desde el exterior, no son fuerzas aplicadas,
sino que aparecen en todos los cuerpos deformables, por su propia naturaleza. Lo
que es importante comprender es que la u
nica accion que la parte verde realiza
sobre la roja es la de unas fuerzas de superficie aplicadas sobre S, y que esto es
cierto para cualquier cuerpo deformable y cualquier superficie interior que queramos

19

considerar1 . Las fuerzas que se transmiten en esta superficie, por unidad de area,
reciben el nombre de tensiones internas, y dependen, en general, del punto del solido
que se investigue y del corte imaginario que se considere. El concepto de tension
interna y externa estan muy relacionados y resumimos su definicion:

Definici
on 3.2: El vector tensi
on t en un punto P del s
olido , cuando este se
corta imaginariamente con una superficie S, es el vector de fuerzas por unidad
de superficie que el resto del cuerpo realiza sobre este punto y superficie.
Como el vector tension esta siempre definido sobre una superficie de normal n,
se definen su proyeccion sobre la normal tn y sobre la superficie misma t de la
manera estandar:
k
tn = (t n)n ,
t
(2.6)
n = t tn .
Se definen las componentes intrnsecas de la tensi
on definida sobre una
superficie de normal n como
p
2
2
(2.7)
n = t n ,
| | = |t
n | = |t| n .
Notese que la componente normal n tiene signo, pero que la componente tangencial
| | siempre es positiva, o nula.
En principio, el vector t de tension en un cuerpo depende del punto P sobre el que
se eval
ue y de la superficie que haya cortado (imaginariamente) dicho cuerpo. Con
objeto de simplificar las ecuaciones de la mecanica de solidos deformables Cauchy
propuso la siguiente condicion, que ha pasado ha llamarse el Principio de Cauchy:
el vector t en un punto P que pertenece a una superficie interna solo depende
de P y de la normal n a dicha superficie en P , matematicamente:
t = t(P, n) .

(2.8)

No hace muchos a
nos se demostro que esta hipotesis no es necesaria, sino que se
puede demostrar que as ocurre siempre, y este resultado se conoce como el teorema
de Noll, su descubridor.

Esta aproximaci
on de hecho ignora el efecto de fuerzas volumetricas entre ambas partes,
que en la naturaleza son muy debiles.

20

4. El tensor de tensiones
El resultado fundamental del analisis del equilibrio en cuerpos deformables
se debe al propio Cauchy y se conoce como el teorema de Cauchy, pues tiene
demostracion.

Teorema 4.3: (Teorema de Cauchy) En un cuerpo deformable en equilibrio

existe un campo de tensores T = T (P ) tal que el campo de tensiones y el de fuerzas
de superficie se pueden expresar como
t(P, n) = T (P )T n , Si P
fs (P )

= T (P )T n , Si P t

(2.9)

El tensor T se conoce como el tensor de tensiones.

La razon por la que este resultado es tan importante es que simplifica la
dependencia funcional de la tension t en cualquier punto interior del cuerpo y sobre
cualquier superficie. De ser una dependencia no lineal t = t(P, n), esta pasa a ser
lineal en la normal n y esto tiene unas consecuencias enormes, no solo desde el punto
de vista de calculo, sino tambien en la obtencion de las ecuaciones de equilibrio.
En una base cualquiera B = {e1 , e2 , e3 }, y en particular en la base cartesiana
B c = {i, j, k}, el teorema de Cauchy se puede expresar en componentes

T11 T21 T31
t1
n1
= T12 T22 T32
,
t2
n2


t3
T13 T23 T33 B n3 B
B
.
(2.10)
Txx Tyx Tzx
tx
nx
t
= Txy Tyy Tzy
n
y
y
tz B c
Txz Tyz Tzz B c nz B c
Aunque la expresion del tensor de tensiones cambia seg
un la base, la expresion
(2.9) es valida en cualquier sistema de coordenadas. Esta formula es una expresion
intrnseca, ya que el teorema de Cauchy no hace referencia a ning
un observador ni
sistema de coordenadas.

Observaciones 4.4:
i. Las dimensiones del tensor de tensiones son de F/L2 .
ingeniera se emplean los MPa.

21

Habitualmente en

ii. El tensor de tensiones admite la siguiente descomposicion:

1
T = p1 + s , siendo p = tr[T ], s = T + p1 .
3

(2.11)

El escalar p es la presion asociada al tensor T y s es la tension desviadora.

iii. Cuando conocemos el campo de tensores T , conocemos todo el estado tensional
del cuerpo, incluyendo las tensiones en el contorno.
iv. El campo de tensiones es u
nico.
v. Cuando el campo de tensores no depende del punto, sino que es constante, se
dice que el estado tensional es homogeneo.

n: La demostraci
Demostracio
on del teorema de Cauchy emplea los argumentos

propuestos por el mismo Cauchy, usando el llamado tetrahedro de Cauchy.

Sea un tetraedro diferencial recto centrado en el punto P , con uno
de sus vertices coincidente con el centro de un sistema de coordenadas de base
B = {e1 , e2 , e3 }. La cara opuesta al origen del sistema de coordenadas tiene area
dA y normal n = n1 e1 + n2 e2 + n3 e3 . Las otras tres caras tienen areas
dA1 = n1 dA ,

dA2 = n2 dA ,

dA3 = n3 dA .

(2.12)

Llamando t a la tension sobre el area dA y t1 , t2 , t3 a las tensiones sobre las otras

tres caras se tiene que
t = t(P, n) , t1 = t(P, e1 ) , t2 = t(P, e2 ) , t3 = t(P, e3 ) ,

(2.13)

y por tanto el equilibrio de fuerzas se expresa como:

t(P, n) dA + t(P, e1 )n1 dA + t(P, e2 )n2 dA + t(P, e3 )n3 dA + fv (P ) dV = 0 .
(2.14)
Como las fuerzas volumetricas multiplican a un infinitesimo de orden superior, estas
se pueden despreciar en la suma anterior. Para continuar, tomamos el lmite en la
ecuacion anterior cuando n e1 para obtener
t(P, e1 ) = t(P, e1 ) .

(2.15)

Como este resultado es valido para cualquier base y vector e1 se concluye que
t(P, n) = t(P, n) ,

22

(2.16)

resultado conocido como el corolario de Cauchy. Utilizando este resultado en la

ecuacion (2.14) obtenemos
t(P, n) = t(P, e1 )n1 + t(P, e2 )n2 + t(P, e3 )n3 .

(2.17)

Esta relacion expresa que la dependencia del vector t en la normal n es lineal y por
lo tanto existe un tensor que denominamos T T tal que
t(P, n) = T T (P )n .

(2.18)

Si la normal n coincide con una normal a la superfice exterior del cuerpo

concluimos que
fs (P ) = T T (P )n .
(2.19)

Ejemplo 4.5: El paraleleppedo de la figura de la izquierda esta sometido a un

estado tensional que, en el sistema cartesiano indicado (con unidades de metros),
se representa con la matriz

xy
y2
xz
[T ] = y 2 xz + yz 0 MPa .
xz
0
z2
Calcular la tension normal y tangencial sobre el plano de la derecha (pasa por
tres vertices) en punto central del paraleleppedo. Datos Lx = 4 m, Ly = 5 m,
Lz = 3 m.
z
z

x
y

Consideramos tres vertices del paraleleppedo por los que pasa el plano de la
figura de la derecha. Sus vectores de posicion son rA = 4i, rB = 5j y rC = 3k.
23

El vector unitario normal al plano de la figura es

n=

(rB rA ) (rC rA )
1
(15i + 12j + 20k) .
=
|(rB rA ) (rC rA )|
769

El punto central del paraleleppedo sobre el cual se desea calcular las tensiones
es P (2, 2.5, 1.5) m y la tension en P sobre el plano es:

5
25/4
3
15
210

1
1
0
{t(P, n)} = 25/4 27/4
12 =
699
MPa .
769 20
769 90
3
0
9/4
Las componentes intrnsecas de este vector son:
n = t(P, n) n =

7047
MPa ,
769

| | =

p
|t|2 n2 = 4.86 MPa .
4

4.1. Interpretaci
on fsica de las componentes del tensor de
tensiones
Cada una de las componentes de T en una base tiene un significado especial y
en ingeniera reciben nombres que hacen referencia a su direccion, como veremos a
continuacion. Si escogemos una base cualquiera B = {e1 , e2 , e3 }, el vector tension
en un punto P que act
ua sobre una superfice de normal ei que pasa por dicho
T
punto es ti = T ei . La tensi
on normal a esta superficie es por tanto
i = ti ei = Tii
y las tensiones tangenciales o cortantes a esta superficie son por tanto
ij = ej T T ei = Tij con i 6= j .

(2.20)

En general, la componente Tij es el valor de la tension que act

ua sobre una superficie
de normal ei , en direccion ej . A la tension Tii se le da el smbolo especial i y, cuando
i 6= j, se emplea la notacion ij = Tij .
Lo anterior aplica tambien a las componentes de T en la base cartesiana
B c = {i, j, k}. En esta tenemos las tres tensiones normales
x = Txx = i T T i ,

y = Tyy = j T T j ,

24

z = Tzz = k T T k ,

(2.21)

y las tensiones tangenciales

xy = Txy = j T T i ,

xz = Txz = k T T i ,

yx = Tyx = j T T i ,

...
(2.22)

En la Figura 2.2 se pueden apreciar algunas de las componentes del tensor de

tension sobre un cubo diferencial

5. Ecuaciones de equilibrio
En esta seccion se obtienen las ecuaciones que expresan que un cuerpo
deformable esta en equilibrio.

5.1. Principio fundamental

De la misma manera que Newton establecion en su segunda ley las condiciones
necesarias y suficientes para el equilibrio de una partcula, y que luego se extienden al
caso de los solidos rgidos a
nadiendo la condicion del equilibrio de momentos, existen
un principio fundamental que generaliza los dos anteriores al caso de los cuerpos
deformables. Como tal principio no es demostrable, pero es fundamental pues en el
se basa el estudio del equilibrio de solidos, fluidos y estructuras deformables.
El principio fundamental de la est
atica de los cuerpos deformables
establece que toda region R contenida o igual a un cuerpo deformable esta en
equilibrio estatico. De forma matematica:
Z
Z
Z
fv dV +
t dA +
fs dA = 0 ,
para toda R
(2.23)
R
R\
R
Mas a
un, para que el equilibrio estatico sea completo, el momento resultante
de todas las fuerzas que act
uan sobre R tambien ha de ser nulo, as pues, para un
sistema de referencia arbitrario, podemos escribir:
Z
Z
Z
rfv dV +
rt dA+
rfs dA = 0 ,
para toda R (2.24)
R
R\
R
siendo r el vector de posicion de los puntos en R.
Las expresiones (2.23) y (2.24) tienen forma integral y su forma diferencial es
mucho mas u
til por lo que la obtenemos a continuacion. Existen dos maneras de

25

t(N, k)
t(A, i)

A
O

t(E, j )

t(O, j )
S

t(D, i)
t(S, k)
k
j
i
zz (N )

zy (N )

zx (N )
A

yz (E)

yx (O)
yy (O)

yy (E)

yx (E)

D
yz (O)
zx (S)
S

zy (S)

zz (S)

Figura 2.2:

Vectores tensi
on y sus componentes sobre un cubo diferencial

26

llegar a ellas: de forma grafica y de forma analtica. A continuacion describimos

la u
ltima de ellas y dejamos la primera para desarrollarla en clase.

5.2. Equilibrio de fuerzas

Empleando el Principio de Cauchy, las dos integrales de superficie de la
expresion (2.23) se pueden escribir como:
Z
Z
Z
t dA +
fs dA =
T n dA .
(2.25)
R\
R
R
Empleando el teorema de la divergencia sobre la region R, esta integral de superficie
se puede expresar como:
Z
Z
T
T n dA =
div T T dv .
(2.26)
R
R
Sustituyendo este resultado en la expresion (2.23), se obtiene
Z


div T T + fv dv = 0 .
R

(2.27)

Pero esta u
ltima expresion es valida para cualquier regi
on R del cuerpo, y eso
u
nicamente es posible si el integrando es identicamente cero en todo punto, es decir
div T T + fv = 0 .

(2.28)

5.3. Equilibrio de momentos

La demostracion analtica de la expresion diferencial correspondiente a la
ecuacion (2.24) requiere unas operaciones tensoriales que estan fuera del temario, as
pues solo indicamos el resultado de las mismas. Empleando el principio de Cauchy
de la misma manera que en la seccion anterior, la ecuacion integral del equilibrio de
momentos se puede escribir como:
Z
Z
r fv dV +
r (T T n) dA = 0 .
(2.29)
R
R
Usando (2.28) en la primera integral e integrando por partes la segunda se obtiene

27

(sin demostrar)
Z
Z
Z
T
r (div T ) dv +
hemi[T ] dv
r (div T T ) dv = 0 ,
R
R
R

(2.30)

siendo hemi el operador que obtiene la parte hemisimetrica de un tensor. Cancelando

los terminos con las divergencias se obtiene
Z
hemi[T ] dv = 0 ,
(2.31)
R
y como este resultado ha de ser valido para cualquier region R, concluimos que
el integrando se anula y por tanto la parte hemisimetrica de T ha de ser cero, es
decir, T es un tensor simetrico. Resumiendo, las dos ecuaciones de equilibrio de los
cuerpos deformables se pueden escribir de forma diferencial como
div T + fv = 0 ,

T = TT .

(2.32)

6. Tensiones principales y direcciones principales de

tensi
on
A partir del estado tensional de un cuerpo T = T (P ) podemos calcular el vector
tension t(P, n) que act
ua sobre cada unidad de superficie en el punto P con normal
n. La direccion del vector tension no se conoce a priori, y dependiendo del estado
tensional, puede ser cualquiera.
En particular, surge el interrogante de si, dado la tension en un punto T (P ), se
puede cortar el cuerpo por una superficie tal que el vector tension t(P, n) resulta
tener la misma direccion que la normal n. Para resolver esta cuestion, planteamos
matematicamente:
T (P )n = n .
(2.33)
La cuestion planteada es equivalente a encontrar los autovectores del tensor T , que
por ser simetrico, siempre son ortonormales. Los autovalores asociados seran reales
siempre. Los valores propios de T se conocen como las tensiones principales y
los autovalores como las direcciones principales de tensi
on.
Los tres autovectores son ortonormales y forman una base B p = {v1 , v2 , v3 }
denominada la base principal de tension en el punto P . Las tensiones principales
se suelen indicar como 1 , 2 , 3 y, por convenio, salvo que se indique lo contrario
se tomara siempre de forma que
1 2 3 .
28

(2.34)

La base principal diagonaliza el tensor de tensiones. En otras palabras, la matriz

asociada a T en la base B p es diagonal de la forma

1 0
0
(2.35)
[T ]B p = 0 2 0 .
0
0 3
Las tensiones principales y las direcciones principales de tension son, en general,
distintas en cada punto de un cuerpo deformable. Solo si el estado tensional
es homogeneo todas ellas seran las mismas en todo el cuerpo. Proporcionan
informacion muy u
til para comprender el estado tensional en un punto, y no
tiene sentido hablar de las tensiones principales de un cuerpo o de sus direcciones
principales de tension.

Observaciones 6.6:
i. Cuando dos tensiones principales son identicas se dice que el tensor de tensiones
es cilndrico. En este caso hay un direccion principal, la asociada a la tension
principal distinta, y un plano de direcciones principales. Dos vectores unitarios
cualesquiera, ortogonales entre s, forman, junto con la primera direccion
principal, la base principal de tensiones.
ii. Si las tres tensiones principales son iguales, se dice que el tensor de tensiones
es esf
erico. Cualquier vector es una direccion principal de tension y cualquier
base es una base principal.

7. Representaci
on gr
afica de un tensor de tensiones
Como un tensor (simetrico) es un objeto difcil de interpretar se han propuesto
varias representaciones graficas que proporcionan algo de informacion sobre el mismo
y permite una evaluacion cualitativa sus propiedades. La representacion mas u
til en
mecanica de solidos es el llamado diagrama de Mohr que representa graficamente
en un diagrama cartesiano todos los posibles valores de las componentes intrnsecas
de tension asociadas al tensor de tensiones en un punto. En realidad, se puede
emplear el diagrama de Mohr para representar graficamente las propiedades de
cualquier tensor simetrico, por ejemplo, el tensor de inercia.
La forma de construir el diagrama de Mohr sera la siguiente. Dado un tensor
simetrico T , escogemos un vector unitario cualquiera n1 . El vector tension cuando
el cuerpo se corta con una superficie de normal n1 es t1 = T n1 , cuyas tensiones
29

| |
C2

C1

C3
3

2
Figura 2.3:

Diagrama de Mohr

intrnsecas son n1 y | |1 . En un diagrama cartesiano se puede dibujar el punto

de coordenadas (n1 , | |1 ). Si ahora escogemos un vector unitario distinto n2
podemos, por el mismo razonamiento, calcular las tensiones intrnsecas n2 y | |2 y
dibujarlas tambien en el mismo diagrama cartesiano sobre el punto de coordenadas
(n2 , | |2 ). Este proceso se puede repetir indefinidamente seleccionando siempre
vectores unitarios y marcando en el diagrama cartesiano las tensiones intrnsecas
resultantes. El resultado de este proceso es, en el lmite, una superficie acotada por
tres crculos que tiene la forma de la Figura 2.3.
La demostracion de que la superficie de posibles componentes intrnsecas es la
indicada en la Figura 2.3 es la siguiente. En primer lugar, como | | es no negativo,
la superficie ha de estar siempre en la parte por encima del eje de abcisas. Ademas,
como n y | | son funciones continuas de T y de n, la superficie ha de ser conexa.
Mas a
un, como hay tres tensiones principales u
nicamente, la superficie intersecta el
eje de abcisas u
nicamente en tres puntos, (i , 0), i = 1, 2, 3.
Para continuar, expresamos el tensor T en su base principal y escogemos
vectores normales n que, tambien expresados en la base principal, son de la forma
n = cos v1 + cos v2 + cos v3 . Los angulos , , son los formados por n y
cada una de las direcciones principales. El vector tension expresado en esta base es
t = 1 cos v1 + 2 cos v2 + 3 cos v3 y por tanto sus componentes intrnsecas
30

son de la forma:

1 cos cos
n = 2 cos cos = 1 cos2 + 2 cos2 + 3 cos2 ,

3 cos
cos
(2.36)
| |2 = 12 cos2 (1 cos2 ) + 22 cos2 (1 cos2 ) + 32 cos2 (1 cos2 )
= 12 cos2 sen2 + 22 cos2 sen2 + 32 cos2 sen2 ,
con la restriccion 1 = cos2 + cos2 + cos2 . Si eliminamos cos de la restriccion,
podemos interpretar (2.36) como la ecuacion parametrica de una superficie en el
plano.
Para continuar definimos
b [0, /2] el angulo que satisface cos
b = | cos | y
by
analogamente
b, limitandones a partir de ahora a estudiar el diagrama de Mohr
b 2 + cos
para vectores normales unitarios de la forma n = cos
b v1 + cos v
bv3 . Si
estudiamos el caso lmite, por ejemplo,
b = /2, es decir el lugar geometrico en
el plano (n , | |) correspondiente a aquellos planos de normal perpendicular a v3 ,
b = sen2
obtenemos a partir de (2.36) que cos2
by
n = 1 cos2
b + 2 sen2
b=

2
1
(cos2
b + 1 sen2
b ) + (sen2
b + 1 cos2
b)
2
2

1 + 2 1 2
+
cos 2b
,
2
2
1 2
sen 2b
.
| |2 =
2
=

(2.37)
Interpretamos que una parte del contorno de la superficie que buscamos es un arco
1 2
2
de circunferencia centrado en ( 1 +
2 , 0) y con radio
2 . Esta semicircunferencia
la denominamos C3 , el lugar geometrico de los componentes intrnsecas de tension
en planos cuya normal es perpendicular a v3 . Repitiendo el mismo argumento, pero
escogiendo planos cuyas normales sean perpendiculares a v1 y a v2 obtenemos que el
contorno de la region admisible es la union de C3 con otras dos semicircunferencias
que llamamos C1 y C2 cuyas propiedades son analogas a las de C3 . Estas tres
circunferencias estan indicadas en la Figura 2.3.
De esta construccion se sigue que dado el diagrama de Mohr correspondiente al
estado tensional en un punto, y considerando un plano que corta al solido pasando
b 2 + cos
por dicho punto y con normal n = cos
b v1 + cos v
bv3 , las componentes
intrnsecas del vector tension se representan el el plano (n , | |) en un punto que
31

| |max

| |

3
Figura 2.4:

Representaci
on gr
afica del estado tensional sobre un plano.

debe de estar dentro de la superficie delimitada por C1 , C2 y C3 .

Mas a
un, el diagrama de Mohr permite la construccion inversa, aunque esta
no la demostramos. Dado un punto A dentro de la region comprendida entre las
tres circunferencias de Mohr, es posible determinar de forma grafica el valor de
b

b , ,
b, los angulos que forman la normal al plano cuyas componentes intrnsecas se
corresponden con el punto A. En la Figura 2.5 se indica la construccion geometrica.

Observaciones 7.7:
i. El diagrama de Mohr no permite distinguir las componentes intrnsecas de los
vectores tension correspondientes a normales que forman angulos mayores de
/2 con los ejes principales. En terminos geometricos, el diagrama de Mohr
representa las componentes intrnsecas en planos cuyas normales pertenecen al
primer octante del sistema de coordenadas principal.
ii. El diagrama de Mohr de un estado tensional cilndrico es simplemente una
semicircunferencia que corta al eje horizontal en 1 , 3 . Si el estado tensional es
esferico, el diagrama de Mohr degenera en un punto sobre el eje horinzontal de
coordenada 1 = 2 = 3 .
iii. Los estados tensiones de planos cuya normal es perpendicular la direccion
principal primera v1 se corresponden con los puntos de C1 (analogamente, para
C2 y C3 ).
iv. La mayor tension cortante | |max en un punto se corresponde con el radio de la
mayor circunferencia de Mohr. Ver Figura 2.4.

32

| |
C2

C1

C3

b
3

b
2

2 +3
2

| |
C2

C1

C3
b
3

b
2

1 +3
2

| |
C2

C1

C3

3
Figura 2.5:

b
1 +2
2

b
Obtenci
on gr
afica de los angulos
b , ,
b en el diagrama de Mohr
33

v. El modulo del vector t correspondiente a un plano cuya representacion en el

diagrama de Mohr es el punto A, el la distancia del centro de coordenadas a
dicho punto (ver Figura 2.4). Ademas, el angulo que forma el vector t y la
normal al plano es .

Ejemplo 7.8: El estado tensional en un punto de un solido deformable tiene una

expresion matricial respecto al sistema de coordadas x, y, z que es

3
0

[T ] = 3
1
0 MPa .
0
0
1

Se pide:
1. Determinar las tensiones principales y dibujar el diagrama de Mohr del estado
tensional.
2. Calcular la tension tangencial maxima en el punto e indicar el angulo que
forma la normal del plano correspondiente con los tres ejes principales de
tension.
3. Calcular el angulo que forma el eje x con cada una de las direcciones
principales de tension.

1) Los autovalores de T son las races del polinomio caracterstico

(1 ) ((1 )(1 ) 3) = 0,
es decir 1 = 2 MPa, 2 = 1 MPa, 3 = 2 MPa.
2) La tension tangencial maxima es el radio del crculo de Mohr mayor, es
decir max = 2 MPa. En la figura se calculan graficamente los angulos de la
normal al plano con mayor tension tangencial respecto de los ejes principales
( = 45o , = 90o , = 45o ).

34

| | (MPa)
max = 2 = 90

= 45
= 45
n (MPa)
-2

3) El eje x es normal al un plano sobre el cual act

ua una tension normal
n = x = 1 MPa y una tension tangencial de modulo | | = |xy | = 3 MPa.
Dibujando este punto sobre el diagrama de Mohr encontramos los angulos que
forma el eje x con las direcciones principales ( = 60o , = 90o , = 30o ).
| | (MPa)

= 90

= 60
= 30
n (MPa)
-2

8. Problemas

2.1 Razonar si las siguientes afirmaciones son verdaderas o falsas:

a) Las tensiones principales en un punto de un solido son 1 = 8 MPa, 2 = 0
35

MPa, 3 = 5 MPa. Hay alg

un plano tal que la tension normal y tangencial
en el mismo valgan n = 2 MPa y = 3 MPa? y cuando n = 5 MPa y
= 1 MPa?
b) Un solido posee un estado tensional esferico con 1 = 2 = 3 = p. Su
diagrama de Mohr corresponde a tres crculos identicos de centro en el origen
y radio p.
c) Un solido deformable cuyo tensor de tensiones sea el mismo en todo punto
no puede estar sometido a ninguna fuerza volmetrica.
d) Si tomamos g = 10 m/s2 ,
1 MPa = 1kp/cm2
e) La u
nica fuerza volumetrica que puede actuar sobre un cuerpo es el peso.

2.2 Un solido deformable en equilibrio se encuentra sometido a un campo de

tensiones cuya expresion matricial, en un sistema cartesiano de coordenadas es:

3x2 2xy
2y 2
[T ] = yx
x2 (y 2 + x2 ) .
zx
zy
(z 2 + xz)
Cuanto valen x , zx y zy ?
Cuales son las dimensiones de ?
Cual es el campo de fuerzas volumetricas que act
ua sobre el cuerpo?
Cual es la tension que se aplica desde el exterior en un punto de la superficie
(x, y, z) = (1, 1, 2) m en el que la normal tiene valor n = 113 (2i + 3j).
Cuales son sus componentes intrnsecas?
e) Encuentra las tensiones principales en el punto (x, y, z) = (1, 0, 1) m si
= 2 106 N/m4 .

a)
b)
c)
d)

2.3 Comprueba que en los siguientes casos el estado tensional propuesto verifica las
ecuaciones de equilibrio en el interior y en el contorno del cuerpo.
a) Un solido deformable sometido a presion exterior uniforme de valor p tiene
tensor de tensiones T = pI, siendo I el tensor unidad.

36

b) Un cilindro cuyo eje de revolucion coincide con k se encuentra sometido a

una presion p en sus caras planas y a presion q en su superficie lateral tiene
un estado tensional cuya expresion matricial en el sistema cartesiano es:

q 0
0
[T ]xyz = 0 q 0
0
0 p
c) Un paraleleppedo como el de la
figura tiene lados de dimensiones
a, b y c. Este cuerpo se encuentra
sometido a tensiones normales
sobre sus caras uniformemente
repartidas cuyas resultantes son,
como indica la figura, F1 , F2 , F3 .
El estado tensional del cuerpo es:
F1

0
0
bc
[T ] = 0 Fac2 0
0
0 Fab3

F3
a

F2

F1

F1

F2
F3

z
y

2.4 Un cuerpo deformable esta sometido a un estado tensional homogeneo cuya

representacion matricial en un sistema cartesiano de coordenadas es:

3/2
0
5/2

[T ] = 3/2
3/2
0 MPa .
0
0
3
1) Dibujar el diagrama de Mohr correspondiente al estado tensional.
Si cortamos el cuerpo con un plano de vector normal n tal que el vector tension t
actuando sobre dicho plano solo tenga componente tangencial:
2) Calcular el modulo de todos los posibles vectores t.
3) Indicar los angulos que forman, en cada caso, n con la segunda direccion
principal.

37

2.5 Un cilindro deformable de longitud

L y radio a esta sujeto por un extremo
y sometido a la accion de una fuerza F
normal a su superficie libre y repartida
uniformemente sobre ella.
Sabiendo que el cilindro esta en
equilibrio, postular una expresion para
el tensor de tensiones y comprobar que
verifica las ecuaciones de equilibrio en el
interior y el en contorno.
Ademas, encontrar el valor del vector tension actuando en los siguientes puntos
y planos, indicando sus componentes normal y tangencial,
a)
b)
c)
d)

Punto
Punto
Punto
Punto

A = Lj, normal n = j.
B = L2 j, normal n = j.
B, normal n = i.

B, normal n = 22 (j + k).

2.6 Si las tensiones principales en un punto de un solido elastico valen 30, 10 y

-10 MPa respectivamente, se pide:
Hallar el modulo del vector tension para la direccion en la que aparezca max .
Hallar las componentes intrnsecas de la tension en la direccion que forma 30o
con la direccion principal 1 y 60o con la direccion principal 2.

38

2.7 Un solido deformable de volumen V

se encuentra en equilibrio, sumergido en
un fluido de peso especfico y sometido
a una fuerza volumetrica. El estado
tensional del cuerpo se expresa, en el
sistema de coordenadas de la figura, con
una matriz de expresion:

z
0
0
z
0 = z [I] .
[T ] = 0
0
0
z

z
V

i) Encontrar el valor del campo de fuerzas volumetricas actuando sobre el cuerpo.

ii) Comprobar que cada punto de la superficie del cuerpo esta sometido a la accion
de la presion hidrostatica y que esta no tiene componente tangencial.
iii) Encontrar el valor de la resultante de las fuerzas ejercidas por el fluido sobre el
cuerpo, demostrando que es el vector opuesto al peso del volumen desplazado,
es decir, el principio de Arqumedes.

2.8 En la figura aparecen dibujados tres estados tensionales, con las tensiones
en MPa. Para cada uno de ellos,
a) Dibuja las componentes de los vectores tension en las caras ocultas de los
cubos diferenciales.
b) Indica la expresion del tensor de tensiones en la base cartesiana de la figura.
5
4

z
x

10

2
3

39

2.9 En la figura aparecen dibujados tres estados tensionales, con las tensiones
en MPa. Para cada uno de ellos,
a) Dibuja las componentes de los vectores tension en las caras ocultas de los
cubos diferenciales.
b) Indica la expresion del tensor de tensiones en la base cartesiana de la figura.
5
4

z
x

10

2
3

2.10
Un
paraleleppedo
deformable tiene las dimensiones
en metros que se indican en la
figura. El cuerpo esta sujeto en
el plano x = 0 y sometido a una
fuerza axial F en direccion del eje
x, aplicada uniformemente sobre
la cara x = 1 m. Cuando se aplica
una fuerza, el cuerpo se parte por el plano x + y = 0, 5. Para repararlo se emplea
un pegamento que resiste tensiones normales y tangenciales maximas de valor:
max = 40 KPa ,

max = 80 KPa .

1) Postula la expresion matricial del tensor de tensiones T (x, y, z) (no es necesario

verificar que cumple las ecuaciones de equilibrio).
2) Calcula las componentes intrnsecas de tension sobre el plano x + y = 0, 5.
3) Cual es el valor de la fuerza maxima que puede resistir el cuerpo reparado
suponiendo que la zona mas debil del mismo es la superficie pegada?

40

2.11 El solido deformable de la

figura se encuentra sometido a
un estado tensional T (x, y, z) del
cual se sabe que
z = y = xz = yz = 0 ,
x = y(2 + (1 x)) MPa ,
y que la tension xy no depende
de la variable z. Las caras del
cuerpo perpendiculares a los ejes y y z estan libres de tensiones y el cuerpo no
esta sometido a ninguna fuerza volumetrica.
1) Demuestra que xy = 21 (y 2 0.152 ).
2) Dibuja el diagrama de Mohr del estado tensional en el punto (x, y, z) = (0, 0, 0).
(nota: cotas en metros. El sistema de coordenadas esta situado en en centro de
la cara menor y sus ejes son paralelos a las aristas del solido)

2.12 Las tensiones principales en un punto son 2, 5 y 7 MPa, respectivamente.

Identifica graficamente la tension cortante maxima y mnima (en valor absoluto) de
aquellos planos en los que:
a)
b)
c)
d)

La tension normal es 6 MPa.

La tension tiene modulo 6 MPa.
La tension normal es el doble de la tension tangencial.
El angulo que forma su normal con el eje principal primero es de 30o .

41

2.13 El paraleleppedo de la
figura esta sometido a un campo
de tensiones cuya respresentacion
matricial en el sistema indicado
es:
2

x
zy 2 z 2
[T ] = zy 2 x
0 MPa .
2
z
0 xz
Si las dimensiones de la figura son metros,
1) Calcular la fuerza que se ejerce sobre la cara x = 1 m del solido.
2) Calcular la fuerza volumetrica total que se realiza sobre el cuerpo.
3) En que punto de la cara x = 1 m aparece la mayor tension tangencial?
(nota: cotas en metros)

2.14 Un cuerpo deformable ocupa un volumen V y esta sometido a un estado

tensional que, referido a una base cartesiana, es de la forma

x2 y
yz 2
2xyz z 3 /3
,
yz 2
y 3 /3
y2z
[T (x, y, z)] = K
3
2
2xyz z /3
y z
0
siendo K una constante con dimensiones de F/L5 . Comprobar que la resultante de
todas las fuerzas que act
uan sobre el contorno del cuerpo es nula.

2.15 Un punto de un cuerpo deformable esta sometido a un estado tensional cuyo

diagrama de Mohr se muestra en la figura. Sobre este diagrama se indica un punto
que se corresponde con las componentes intrnsecas de tension sobre un plano de
normal n.
b
1) Determina graficamente el valor de
b , ,
b.
2) Calcula el vector tension t sobre dicho plano.
3) Indica que angulo forma t con n.

42

| | (MPa)

/6
-10

10

20

43

30

40

n (MPa)

44

Captulo 3

An
alisis de la deformaci
on
de los cuerpos

1. Introducci
on
En el captulo 2 se ha estudiado la distribucion de fuerzas internas en los cuerpos
deformables, que es una de las diferencias fundamentales con los cuerpos rgidos.
En este captulo nos ocupamos del estudio de la deformabilidad en s misma, como
concepto geometrico. En primer lugar la propondremos una definicion matematica
de dicho concepto y nos ocuparemos, a continuacion, de calcularla exactamente en
todos los puntos de un cuerpo, suponiendo que las deformaciones son peque
nas.
La hipotesis de peque
nas deformaciones limita la aplicabilidad de los conceptos

45

presentados en este captulo, pero no se debe de olvidar que la gran mayora de las
estructuras y de las maquinas trabajan en situaciones de peque
nas deformaciones
as pues sus aplicaciones son numerosas. El estudio de las llamandas deformaciones
finitas o grandes deformaciones queda para cursos mas avanzados.
Una vez que el concepto de deformacion este claramente definido estudiaremos,
en el captulo 4, que existe una relacion, llamada constitutiva, entre las tensiones y las
deformaciones, cerrando con ello el planteamiento del modelo del cuerpo deformable.

2. Cinem
atica de un cuerpo deformable
Un cuerpo sometido a fuerzas y acciones externas puede responder de manera
que sus puntos cambien de posicion. Dado un punto P , denominamos P 0 la
posicion en el espacio que el punto P ocupa despues de que el cuerpo sufra las
acciones exteriores. Definimos el campo de desplazamientos u : R3 como
u(P ) = rP 0 rP

(3.1)

siendo rP 0 y rP los vectores de posicion de los puntos P 0 y P respectivamente

respecto de un sistema de coordenadas cualquiera, pero fijo. El concepto de campo
de desplazamiento se estudia en el modelo de los cuerpos rgidos, pues no es exclusivo
de los cuerpos deformables.
Sin embargo, intuitivamente, se entiende que el concepto de deformacion, aunque
relacionado con el desplazamiento, no es lo mismo. Un cuerpo en el que todos
sus puntos se desplazan con un campo u constante sabemos que no se deforma.
Tampoco se deforma un cuerpo que rota alrededor de un eje o de un punto fijo.
La definicion precisa del concepto de deformacion no es evidente. Los cuerpos
tienen muchas maneras de deformarse, si se miran en su conjunto. Sin embargo,
localmente (a nivel diferencial) veremos que solo hay dos modos de deformacion. Lo
que define a la deformacion es el desplazamiento relativo, y para ello es necesario
estudiar el campo cinematico en el entorno de cada punto, porque, a diferencia de lo
que ocurre en los solidos rgidos, cada entorno se puede deformar de forma distinta.
As, en una estructura sometida a cargas se puede hablar de que una viga esta muy
deformada y otra no; en el chasis de un vehculo, una parte esta muy deformada y
otra a penas, etc.
A pesar de la aparente dificultad en describir la deformacion en toda
generalidad, argumentamos a continuacion que a nivel local solo existen dos tipos
de deformaciones. El primer tipo de deformacion que existe se refiere al cambio de

46

Q0
P0

dr 0

u(Q)
u(P )

dr

P
Figura 3.1:

Deformaci
on en el entorno de un punto de un cuerpo deformable

longitud y se conoce como la deformaci

on longitudinal unitaria de un cuerpo
en un punto P y en una direccion . Para poder definirlo, consideremos dos puntos
P y Q infinitesimalmente proximos sobre un cuerpo deformable, y llamemos dr
al vector diferencial que va desde P a Q. Cuando el cuerpo se deforma, estos dos
puntos pasan a ocupar las posiciones P 0 y Q0 , y denominamos ahora dr 0 al vector
diferencial que los une. Es el cambio relativo de longitud de un elemento diferencial
de lnea con origen en el punto P y con direccion . Su definiticion matematica es
ex (P, ) :=

| dr 0 | | dr|
.
| dr|

(3.2)

El segundo tipo de deformacion posible es el cambio de angulo. Para definirlo,

consideremos ahora que sobre el cuerpo deformable existen dos vectores diferenciales
dr 1 y dr 2 con origen en el punto P . Cuando el cuerpo se deforma, estos vectores se
transforman en nuevos vectores diferenciales dr 01 y dr 02 . Cada uno de ellos puede
haber cambiado su longitud debido a la deformacion longitudinal, pero ademas
puede haber cambiado el angulo que forman entre s. Llamando 1 y 2 a las
direcciones de los vectores dr 1 y dr 2 , se define la deformaci
on angular como
ex (P, 1 , 2 ) := cos(0 ) cos() ,

(3.3)

siendo el angulo (en radianes) que forman dr 1 , dr 2 y 0 el angulo que

forman dr 01 , dr 02 .

47

2
dr20
dr10
0
P0

dr2
P

dr1

Deformaci
on angular en el entorno de un punto de un cuerpo
deformable

Figura 3.2:

Observaciones 2.1:
i. Tanto las deformaciones longitudinales como las angulares no tienen
dimensiones.
ii. Ambas medidas son locales, es decir relativas a un punto, no a un cuerpo en su
conjunto. Ademas, estan relacionadas con direcciones. No tiene sentido hablar
de la deformacion longitudinal de un cuerpo ni de la deformacion longitudinal
de un punto. Tan solo el concepto de la deformacion longitudinal de un punto,
en una direccion tiene sentido.
iii. Las definiciones de deformacion presentadas en esta seccion son validas incluso
para deformaciones finitas y por ello se denominan exactas.

48

3. C
alculo de deformaciones. El tensor de deformaci
on
infinitesimal
El problema que pretendemos resolver en esta seccion es el siguiente: dado un
campo de desplazamientos u : R3 , cuales son las deformaciones ex y ex en
todos los puntos y todas las direcciones posibles? Este es el problema central de la
cinematica de los cuerpos deformables.
Para calcular las deformaciones en cualquier punto sera necesario determinar la
forma local del campo de desplazamientos alrededor de dicho punto. Como siempre
en teora de campos, esta informacion la recoge el gradiente:

Definici
on 3.2: Dado un campo de desplazamientos u : R3 se define el
tensor gradiente de desplazmientos u como aquel que verifica
u(P + dr) = u(P ) + u(P ) dr + O(| dr|2 ) .

(3.4)

La expresi
on en coordenadas cartesianas de la matriz asociada al tensor u es

ux,x (P ) ux,y (P ) ux,z (P )

[u(P )] = uy,x (P ) uy,y (P ) uy,z (P ) .
uz,x (P ) uz,y (P ) uz,z (P )
El gradiente de desplazamientos es tambien adimensional y, como veremos
despues, nos servira para calcular deformaciones. Para simplificar el calculo de las
mismas vamos a suponer a partir de ahora que el cuerpo al desplazarse se deforma
muy poco. La definicion precisa de que significa esto es la siguiente:

Definici
on 3.3: Se dice que un cuerpo experimenta una deformacion pequena
si kuk  1. Esto ocurre si y s
olo si todas las componentes de u son mucho
m
as peque
nas que 1.

3.1. El tensor de deformaciones infinitesimales

Cuando calculemos deformaciones comprobaremos que estas solo dependen
de la parte simetrica de u y a este objeto lo denominaremos el tensor de
deformaci
on, y juega un papel central en el modelo del solido deformable.

49

Definici
on 3.4: Dado un campo de desplazamientos u : R3 , definimos la
deformaci
on infinitesimal D como el campo de tensores simetricos


D(P ) := S u(P ) = 21 u(P ) + uT (P ) .

(3.5)

La parte de u que no esta asociada a la deformacion infinitesimal D, es decir

la parte hemisimetrica del tensor, s que esta asociada al movimiento local y recibe
la siguiente definicion:

Definici
on 3.5: La parte hemisimetrica de u(P ) es el campo tensorial de giro
infinitesimal
:= a u(P ) =

1
2

u(P ) uT (P ) .

(3.6)

Como es un tensor hemisimetrico tiene un vector axial asociado , llamado

el vector de giro infinitesimal. Este campo vectorial satisface adem
as
= 21 rot u .

(3.7)

La interpretacion geometrica completa de estos campos tensoriales es la

siguiente. Si en un punto P se escogen tres vectores diferenciales ortogonales
dr 1 , dr 2 , dr 3 , cuando el cuerpo se deforme, estos tres vectores cambian de modulo
y direccion transformandose en tres nuevos vectores infinitesimales dr 01 , dr 02 , dr 03 .
Para cada uno de ellos se puede escribir
dr 0i = (I + D + ) dr i = dr i + D dr i + dr i .

(3.8)

As pues, los tensores D y caracterizan, de forma completa, la transformacion

geometrica local, para cada entorno diferencial de los puntos del cuerpo deformable.

3.2. C
alculo de deformaciones longitudinales
Para obtener una expresion que nos permita obtener el valor de ex en funcion de
u y sus gradientes, sustituimos el desarrollo de Taylor del campo de desplazamiento
en la expresion (3.2). Sea el vector unitario en la direccion en la que queremos

50

calcular la deformacion longitudinal. Entonces,

| dr + u(Q) u(P )|
| dr + u(P ) dr|
| dr 0 |
1=
1=
1
| dr|
| dr|
| dr|
p
= | + u(P )| 1 = + 2 u(P ) + |u(P )|2 1 .

ex :=

(3.9)

La expresion para la deformacion longitudinal ex es una funcion no lineal. Sin

embargo, si las deformaciones son peque
nas podemos aproximar
p
ex + 2 u(P ) 1 ,
(3.10)

Y utilizando un desarrollo de Taylor para la funcion 1 + x obtener finalmente

ex u(P ) .

(3.11)

Si A es un tensor cualquiera, es immediato comprobar que A = AS y

concluimos

Definici
on 3.6: Se define la deformaci
on longitudinal infinitesimal en un
punto P y una direcci
on cualquiera como el escalar
(P, ) := D .

(3.12)

Observaciones 3.7:
i. La deformacion longitudinal infinitesimal es una aproximacion al la verdadera
deformacion longitudinal ex , que es mucho mas complicada de calcular. La
aproximacion es tanto mas valida cuanto mas peque
na sea la cantidad kuk.
Por tanto, solo es exacta cuando la deformacion sea infinitesimal. Para
deformaciones finitas se puede dar el caso de que un cuerpo que se mueve
rgidamente tenga deformacion no nula.
ii. La deformacion es unitaria, y por tanto adimensional.
Cuando un cuerpo se deforma, una curva material C se deforma tambien pues
cada uno de sus puntos se desplaza debido al movimiento del cuerpo. A menudo
es interesante encontrar la longitud de la curva deformada a partir de la longitud
inicial y de la deformacion longitudinal unitaria en cada punto. Si la longitud de
la curva sin deformar es L, cada punto de la curva lo denominamos P y el vector
51

tangente a la curva en P es , entonces

Z
0
L = (1 + (P, (P ))) dS .
C

(3.13)

Ejemplo 3.8:
Un cuarto de aro de radio r se deforma
seg
un el campo de desplazamientos
u(x, y, z) =

x2
i,
2

siendo x un eje del sistema cartesiano

situado en el centro del aro como se
indica en la figura. Calcular:

i) la deformacion longitudinal unitaria

en cualquier punto del aro y direccion
circunferencial.

y
x

ii) la longitud del aro deformado.

El vector tangente al aro en un punto generico es = sen i + cos j, siendo

[0, /2] el angulo que forma el vector de posicion del punto con el eje x. La
deformacion longitudinal unitaria en dicho punto y direccion es:
(P, ) = D = x sen2 .
Como x = r cos el valor de la deformacion es simplemente (P, ) =
2 r sen2 cos . La longitud del trozo de aro deformado es:
L0 =

Z
C

dS 0 =

/2

/2

(1 + (, )) dS =
=0

(1 + r sen2 cos ) dS =

=0

r + r2 .
2
3
4

3.3. C
alculo de deformaciones angulares
El calculo de la deformacion angular es, como en el caso anterior, complejo en
general. Sin embargo, bajo algunas hipotesis simplificadores dicho calculo se vuelve
mas sencillo.
52

dr20
0

dr10

dr2
= /2
P
Figura 3.3:

dr1

Interpretaci
on geometrica de la deformacion angular infinitesimal

Para empezar, u
nicamente consideraremos el caso de que 1 y 2 sean
ortonormales. A partir entonces de la expresion (3.3) de la deformacion angular
se obtiene
dr 01 dr 02
(I + u) dr 1 (I + u) dr 2 | dr 1 | | dr 2 |
=
| dr 1 | | dr 2 |
| dr 01 | | dr 02 |
| dr 01 | | dr 02 |
1
= (I + u)1 (I + u)1
(1 + 1 )(1 + 2 )

ex (1 , 2 ) := cos(0 ) =

= (1 2 + u1 2 + 1 u2 + u1 u2 ) (1 + 1 )(1 + 2 ) .
(3.14)
Si las deformaciones son peque
nas, la expresion anterior se puede simplificar
resultando
ex (1 , 2 ) u1 2 + 1 u2
(3.15)
= u1 2 + T u1 2 = 21 D2 .
Basandose en la anterior aproximacion se propone la siguiente

53

Definici
on 3.9: Se define la deformaci
on angular infinitesimal en un punto
P entre dos direcciones ortogonales 1 , 2 como

(P, 1 , 2 ) := 1 D2 .
2

(3.16)

Cuando 1 y 2 son ortonormales, la expresion (3.3) indica que ex es el coseno

del angulo 0 una vez que el cuerpo esta deformado, o equivalentemente, el seno
del angulo
1 +
2 de la Figura 3.3. Cuando las deformaciones son peque
nas
sen(
1 +
2 )
1 +
2 . En conclusion, para peque
nas deformaciones, la medida
de deformacion angular es aproximadamente igual al angulo que se cierra . Esta
aproximacion es tanto mejor cuando menor sea la deformacion y es la interpretacion
habitual en ingeniera de la medida de deformacion .

3.4. Interpretaci
on geom
etrica de las componentes del tensor de
deformaci
on
Las componentes de la matriz asociada a D en una base ortogonal cualquiera
B = {e1 , e2 , e3 } son Dij con i, j = 1, 2, 3.
A partir de las definiciones de la deformacion longitudinal y angular se
comprueba que ii = ei Dei , para cualquier i = 1, 2, 3. Es decir, la componente
Dii es la deformacion longitudinal unitaria en direccion de ei . Para cada pareja
i, j = 1, 2, 3 con i 6= j, tambien se comprueba inmediatamente que ij /2 = ei Dej ,
es decir, que la componente Dij de la matriz de deformaciones es la deformacion
angular que sufre el angulo comprendido entre ei y ej . Por estas razones es habitual
escibir la matriz [D] como

11 212 213
[D] = 221 22 223 .
(3.17)
31
32

33
2
2
En el caso de que la base ortongonal sea la cartesiana, la matriz [D] es escribe de
la siguiente forma:

xy
xz
xx
2
2

yz
[D] = yx
yy
.
(3.18)
2
2
zy
zx
zz
2
2
Las compontentes Dij , i 6= j son las medidas de deformacion que aparecen en el
tensor D. Sin embargo las deformaciones ij no aparecen tal cual en D, aunque su

54

uso es mas habitual en ingeniera que las primeras. Se suele decir que ii y ij son
medidas de deformacion ingenieriles.

Observaciones 3.10:
i. Cuando D no depende del punto de evaluacion el estado de deformacion del
cuerpo es homogeneo.
ii. Cuando un cuerpo experimenta un desplazamiento u y el tensor D es nulo en
todo punto no significa que el cuerpo no haya sufrido deformacion alguna. Esto
sera rigurosamente cierto u
nicamente si la deformacion es infinitesimal. Para
ser precisos diremos que el campo de desplazamientos es infinitesimalmente
rgido. Cuan rgido es realmente depende de la peque
nez del gradiente de
deformaciones. Por ejemplo, su un cuerpo sufre una rotacion finita alrededor de
un punto fijo sabemos que este no se ha deformado. Sin embargo se comprueba
que el campo de deformaciones D no se anula. En esta situacion el gradiente
u no es peque
no as que D no es una medida fiable de deformacion.
iii. Como las definiciones (3.5) y (3.16) son tensoriales, son validas en cualquier
sistema de referencia que se emplee, no necesariamente cartesiano.
iv. Las expresiones de cambio de base para las componentes de un tensor descritas
en el captulo 1 se aplican directamente al tensor D.

Ejemplo 3.11:
El cuadrado de la figura sufre un
estado de deformacion homogeneo con
z = xz = yz = 0. Se sabe que, al
deformarse,

a) el lado AC aumenta su tama

no un 1%.
b) el angulo ACD se cierra 2o .

c) el segmento BD se acorta un 3%.

Encontrar la expresion matricial del
tensor de deformaciones en el sistema
xy y, a partir de este, calcular cuanto
se abre el angulo EF D.

y
60
A

xB

A partir del dato del apartado a) se deduce que x = 0.01. Tambien,

directamente a partir de b) se sigue que xy = /90. La deformacion

55

longitudinal unitaria en direccion de BD se calcula de la siguiente manera

1
1

0.01 180
0 2
2

3
3
0.03 =
180
y
0
2
2
0
0
0
0
0
y, a partir de esta expresion, se obtiene y = 2.3 102 .
Por u
ltimo, para encontrar la deformacion angular del angulo EF D se utiliza la
formula

2
2
(1 , 2 ) = 1 D2
con
1 =
(i + j) , 2 =
(j i) .
2
2
2
Operando con los datos obtenidos anteriormente se calcula

1
0 1
0.01
180

2
2

2.3 102 0
(1 , 2 ) =
1 180
1
= 0.017 .

2
2
2
0
0
0
0
0
Concluimos que = 3.3 102 , es decir, que el angulo EF D se abre 3.3 102
radianes, es decir, 1.9o .
4

4. La deformaci
on volum
etrica
Otra caracterstica propia de los cuerpos deformables es que su volumen cambia
cuando estos se someten a fuerzas y/o acciones exteriores. Para cuantificar este
efecto se propone tambien una medida local de deformacion volumetrica. En un
cuerpo deformable, cada diferencial de volumen dV se transforma debido a la
deformacion en un nuevo diferencial dV 0 . Se define la deformaci
on volum
etrica
relativa como el cambio relativo local de volumen:
ex (P ) =

dV 0 dV
.
dV

(3.19)

El calculo de la deformacion volumetrica exacta ex es complicada, al igual que

ocurra con la deformacion longitudinal y la angular. Tambien como en estos
casos el calculo se puede simplificar cuando las deformaciones son infinitesimales.
Para encontrar la aproximacion necesaria se consideran tres vectores diferenciales

56

dr 1 , dr 2 , dr 3 sobre un punto P , no coplanarios. Entonces,

dV 0
[ dr 01 dr 02 dr 03 ]
1=
1
dV
[ dr 1 dr 2 dr 3 ]
[(1 + u) dr 1 (1 + u) dr 2 (1 + u) dr 3 ]
=
1
[ dr 1 dr 2 dr 3 ]
det(1 + u)[ dr 1 dr 2 dr 3 ]
1
=
[ dr 1 dr 2 dr 3 ]
= det(1 + u) 1 .

ex (P ) =

(3.20)

La funcion determinante es no lineal y se puede aproximar det(1 + u) 1 + div u

cuando la deformacion es peque
na. Por lo tanto concluimos que
ex (P ) div u(P ) .

(3.21)

As pues la divergencia del desplazamiento es una aproximacion al cambio relativo

de volumen y se define

Definici
on 4.12: La deformaci
on volum
etrica infinitesimal en un punto
P es
(P ) := div u(P ) .
Adem
as, se puede demostrar que tr[D] = div u.
Como en el caso del resto de deformaciones infinitesimales, la cantidad es una
aproximacion de la verdadera deformacion infinitesimal ex , cuya validez es tanto
mayor cuanto mas peque
na sea la deformacion.
Empleando el concepto de deformacion volumetrica se puede descomponer el
tensor D en dos contribuciones aditivas. La primera parte, la volumetrica, es un
tensor de deformacion que tiene deformacion volumetrica igual a . La segunda
parte, llamada desviadora, tiene traza nula y por lo tanto es una deformacion
responsable de cambio de forma, pero no de volumen:

D = I +e ,
3

(3.22)

y calculando e = D 3 I.
Para calcular el cambio de volumen que experimenta un cuerpo entero o una
region el mismo, basta con integrar la deformacion volumetrica en dicho dominio.

57

As pues,
Z
V =

(P ) dV .
(3.23)
R
Solo si la deformacion volumetrica es uniforme se puede escribir V = V .

Ejemplo 4.13: Un cubo deformable de lado L experimenta un campo de

desplazamientos
z y2
x3
z3
u(x, y, z) = sen( ) i + 2 j +
k,
L L
L
3L2

(3.24)

referido a un sistema de coordenadas cartesiano con origen en el centro del cubo

y ejes coordenados paralelos a las aristas. Calcular el incremento de volumen
del cubo en la deformacion.
La deformacion volumetrica relativa es
(x, y, z) = div u(x, y, z) =

z2
,
L2

(3.25)

as que el incremento total de volumen del cubo sera:

Z L/2 Z L/2 Z L/2
V =
(x, y, z) dV
Z

x=L/2 y=L/2 z=L/2

Z L/2 Z L/2
L/2

=
=

x=L/2
L3

12

y=L/2

z=L/2

z2
dx dy dz
L3

(3.26)

5. Deformaciones principales y direcciones principales

de deformaci
on
El tensor de deformacion infinitesimal D es simetrico y por tanto posee tres
autovalores reales 1 2 3 denominados las deformaciones principales. Sus
autovectores asociados {w1 , w2 , w3 } forman un base ortonormal y se conocen como
base principal de deformaci
on.

58

Observaciones 5.14:
i. Estos valores y vectores propios dependen de cada punto del cuerpo a menos
que la deformacion sea homogenea, en cuyo caso seran los mismos.
ii. La matriz asociada al tensor de deformaciones en la base principal es diagonal,
y las deformaciones principales son las componentes de la diagonal.
iii. La deformacion longitudinal unitaria en la direccion principal wi es:
(P, wi ) = wi D(P )wi = wi i wi = i .

(3.27)

iv. La deformacion angular que sufre el angulo comprendido entre las direcciones
principales i, j con i 6= j es:

(P, wi , wj ) = wi Dwj = wi j wj = 0 ,
2

(3.28)

con lo que se demuestra que en angulo entre vectores de la base principal no

cambia con la deformacion. Cuando el cuerpo se deforma, la terna principal se
transforma en otra terna de vectores ortogonales.

6. Galgas extensom
etricas
Las galgas extensometricas son unos sensores que se emplean para medir las
deformaciones longitudinales unitarias en el contorno de los cuerpos deformables.
El mecanismo de medida es resistivo: las galgas son unas resistencias muy peque
nas
que se adhieren a la superficie del cuerpo con un pegamento. Cuando el cuerpo se
deforma, la galga pegada se deforma solidariamente al cuerpo y esto hace que vare
su resistencia. Un ohmmetro calibrado relaciona el cambio de resistencia con la
deformacion experimentada en el punto donde la galga estaba pegada.
Por supuesto, para que las galgas proporcionen una medida precisa, estas no
pueden ser grandes y el pegado ha de ser muy bueno. Hay empresas especializadas
que venden galgas extensometricas de distinto tama
no y calidad, todas ellas
calibradas.
Debido al mecanismo de medida, las galgas solo pueden medir deformaciones
longitudinales y solo en la superficie de los cuerpos. Dada una direccion cualquiera
tangente a la superficie del cuerpo, la deformacion longitudinal unitaria que mide
una galga orientada seg
un una direccion es:
= D .

59

(3.29)

Mediante la colocacion de varias galgas superpuestas sobre un mismo punto

(formando rosetas) se pueden medir tambien deformaciones angulares. Para
demostrar esto supongamos que escogemos un sistema de coordadas cartesiano
(x, y, z) en un punto de la superficie, siendo el eje z normal a la superficie. Si
se pegan tres galgas sobre el punto, formando, respectivamente, angulos , y
con el eje x, tenemos que cada una de ellas mide una deformacion longitudinal
unitaria distinta:
= D ,

= D ,

= D ,

(3.30)

siendo , y los vectores unitarios en direccion de las galgas. Expresando el

tensor de deformaciones D y los vectores de direccion de las galgas en el sistema
escogido (x, y, z) obtenemos

xy
xz
cos x
cos
2
2

yz
= sen xy
y
sen = x cos2 + xy sen cos + y sen2 ,
2

2xz yz

0
z
0
2

xy
xz
cos
cos x
2
2
yz

sen = x cos2 + xy sen cos + y sen2 ,

= sen xy
y
2

2xz yz

0
z
0
2

xy
xz
cos x
cos
2
2
yz

= sen xy

sen = x cos2 + xy sen cos + y sen2 .

y
2

2xz yz

0
0
z
2
2
(3.31)
En las ecuaciones anteriores hay tres incognitas x , y , xy , tres datos geometricos
, y y tres medidas experimentales , y . A partir de los datos podemos
encontrar las deformaciones longitudinales unitarias en las direcciones x e y y la
deformacion angular.

7. El diagrama de Mohr de deformaciones

Como se indico en el captulo 2, el diagrama de Mohr se puede emplear para
representar graficamente algunas de la propiedades de cualquier tensor simetrico y,
en particular, el tensor de deformaciones. El diagrama de Mohr de deformaciones
tiene pues un aspecto identico al de tensiones, si se reinterpreta, para cada direccion
, la componente sobre el eje de abcisas como la deformacion longitudinal unitaria
() y aquella sobre el eje de ordenadas como la maxima deformacion angular
/2(, ) de todos los angulos rectos con un lado igual a .

60

Todas las propiedades geometricas del diagrama de Mohr que se enunciaron para
el tensor de tensiones tienen su correspondencia en el tensor de deformaciones. Por
ejemplo, el radio del crculo de Mohr mayor es la mitad de la maxima deformacion
angular que se da entre todos las parejas de vectores ortonormales que coinciden
en el punto. En cualquier caso, la utilidad practica del diagrama de Mohr de las
deformaciones es mucho menor que en caso de las tensiones.

8. Problemas

3.1 Un solido se deforma y el desplazamiento de un punto cualquiera (x, y, z) del

mismo en un sistema de coordenadas {O, i, j, k} viene dado por la expresion:
u(x, y, z) = 3 103 e

x+y
5

z
i + 2 103 sen( )k
2

(m) .

Se pide:
1) Hallar la representacion matricial del tensor de deformaciones infinitesimales D
en el sistema de referencia {O, i, j, k}.
b en un punto cualquiera y su vector axial
2) Hallar la matriz de giro infinitesimal
1
asociado . Comprobar que = 2 rot u.
3) Considerese un vector diferencial
de modulo dS sobre el punto (x, y, z) =
(1, 2, 1) m con direccion (0, 3/2, 1/2). Calcular la expresion del vector
despues de la deformacion, indicando su punto de aplicacion, direccion y modulo.
4) Hallar las deformaciones principales en el punto O y dibujar el diagrama de
Mohr correspondiente al estado de deformacion en dicho punto.
5) Calcular la deformacion volumetrica en O. Comprobar que es un invariante y
que por tanto
xx (O) + yy (O) + zz (O) = 1 (O) + 2 (O) + 3 (O) = div u .

3.2 Una barra recta de longitud L, alineada con el eje x, sufre una deformacion que
se describe con un campo de desplazamiento



x
x 
u(x, y, z) = 100L sen
x i + 100L 1 cos
j.
100L
100L
61

Calcular la longitud de la barra deformada.

3.3 Un solido de seccion romboidal tiene

una profundidad de 1 m. El campo de
desplazamientos del cuerpo es:
 x
y 
3x
u(x, y, z) =
+
j+
k.
400 150
100
i) Calcular la expresion matricial del tensor
de
deformaciones infinitesimales en el sistema
de referencia de la figura. Comprobar que
es homogeneo.
ii) Dibujar la seccion z = 0 deformada empleando
una escala adecuada para que se puedan
apreciar las deformaciones.
iii) Calcular las deformaciones longitudinales unitarias en un punto generico en las
direcciones i, j y k. Comprobar que estas son independientes de las coordenadas
del punto.
iv) Calcular la deformacion que sufre el angulo formado por los vectores i y j en el
origen de coordenadas.
v) Interpretar graficamente los resultados de iii) y iv)
vi) Calcular el volumen del cuerpo deformado.
vii) Calcular la longitud de los lados deformados.

62

3.4 Una hormiga se encuentra en el centro de un

disco deformable de radio r = 2m, donde se define un
sistema de coordendas cartesiano. El disco se deforma
y sus puntos experimentan un desplazamiento

u(x, y, z) = (3x + 4y) 102 i + (2x 3y) 102 j .

Sabiendo que la hormiga se mueve a una velocidad

constante de 2 m/min, cual es el tiempo mnimo que
la hormiga empleara en alcanzar el borde del disco
deformado?que direccion ha de seguir para ello?

3.5 La deformacion de un paraleleppedo solido de lados a, b, c metros esta descrita

por el vector de desplazamientos:
u(x, y, z) = (3x + 2y) 103 i + (6x + 2z) 103 k (m)
referido a un sistema de coordenadas cartesiano (O, x, y, z).
b
a

F
c

z
x

Se pide
1) Determina la expresion matricial de los tensores de deformaciones infinitesimales
y de giro infinitesimal en el sistema de referencia (O, x, y, z).
2) Calcula el volumen del paraleleppedo deformado.
del paraleleppedo deformado.
3) Calcula la longitud de la diagonal OF

63

3.6 Un cuerpo experimenta el desplazamiento correspondiente a una rotacion

de angulo alrededor del eje z cartesiano. Por ello, un punto de coordenadas
cartesianas (x, y, z) pasa, despues de la deformacion, a ocupar la posicion
(x cos + y sen , x sen + y cos , z) .
Calcula el vector de desplazamiento y la expresion matricial del tensor de
deformaciones en coordenadas cartesianas. Comprueba que la deformacion no es
nula a pesar de que el movimiento que el cuerpo experimenta es de solido rgido.
Explica razonadamente como es esto posible.

3.7 El triangulo isosceles de la figura

corresponde a la seccion z = 0 de un
prisma triangular de altura 2 m que sufre
una deformacion
u(x, y, z) = [(2x + y + K ey )i + yj + (x + y + z)k]103 (m)
siendo K una constante sin dimensiones.
1) Calcula el angulo formado por los lados OA y OB despues de la deformacion.
2) Calcula la longitud del lado AB tambien despues de la deformacion (en este
u
ltimo caso tomar K = 0 y suponer que la longitud de los catetos es de 1 m).
3) Calcula el volumen del prisma deformado suponiendo tambien que K = 0.

64

3.8 El paraleleppedo de la figura tiene un

campo de desplazamientos
u(x, y, z) = A(x+y)i+B(y+z)j+C(z+x)k (m)
siendo A, B y C tres constantes adimensionales.
Encontar su valor sabiendo que:
i) Con la deformacion, el paraleleppedo aumenta su tama
no un 2%.
ii) La diagonal de la cara x = 0 que divide en dos el angulo entre los ejes y, z reduce
su longitud un 1, 5%.
iii) El angulo formado por los vectores i y j en el origen de coordenadas forman
89o despues de la deformacion.

3.9 Una esfera de radio r = 2 m sufre una deformacion al deformarse seg

un un campo
2
2
de desplazamientos u(x, y, z) = (5x + 2y) 10 i + (2x + 3y) 10 j + 4 102 zk.
1) Demostrar que la deformacion es homogenea.
2) Calcular el volumen de la esfera deformada.
3) Calcular la longitud de los radios mas largos y mas cortos despues de la
deformacion.

3.10 Un cubo de lado l = 2 m esta colocado en el centro de un sistema de

coordenadas cartesiano con sus aristas paralelas a los ejes de coordenas. Si el
desplazamiento del cubo es:
u(x, y, z) = (6x3 + z) 102 i + (7y + x) 102 j + (4z + x2 ) 102 k .
1) Calcula el volumen del cubo deformado.
2) Calcula la longitud de la arista que conecta los puntos (1, 1, 1) m y (1, 1, 1)m,
una vez que el cuerpo se deforma.

65

66

Captulo 4

Elasticidad y
termoelasticidad lineal

En los captulos 2 y 3 se han estudiado las dos caractersticas fundamentales de

los cuerpos deformables, a saber, la existencia de fuerzas internas que garantizan su
equilibrio y su deformabilidad. En este captulo se estudia como no se trata de dos
cuestiones aisladas sino que, al contrario, son causa y efecto de una misma realidad.
La relacion entre tension y deformacion es el problema fundamental de la
mecanica de cuerpos deformables, y el mas complejo. Debido a esto u
ltimo nos
concentramos en este captulo en dos casos muy sencillos, el de la elasticidad y
la termoelasticidad lineal. A pesar de su sencillez estos dos modelos constitutivos
permiten resolver una gran cantidad de problemas de ingeniera. Quedan para otros
cursos mas avanzados el estudio de la viscoelasticidad, la plasticidad, etc

67

1. El concepto de elasticidad
El terminoelasticidad se emplea coloquialmente y posee el siguiente significado
preciso en el contexto de la mecanica

Definici
on 1.1: Se dice que un punto material P se comporta el
asticamente
cuando el tensor de tensi
on T (P ) es u
nicamente funci
on del valor en ese instante
de la deformaci
on D(P ). Supondremos por simplificar que si la deformaci
on se
anula, tambien lo hace la tensi
on.
Bajo esta definicion los efectos de la historia y de la velocidad de deformacion
son irrelevantes, fenomenos por otra parte fundamentales para el estudio de la
plasticidad y la viscoelasticidad.

Observaciones 1.2:
i. No existe ning
un material elastico en la naturaleza. Lo que s existen
son materiales que se comportan elasticamente en un rango de deformacion,
normalmente peque
no. Este rango es el habitual para el funcionamiento de las
maquinas y las estructuras por lo que basta para muchos calculos en ingeniera
mecanica.
ii. En un ciclo de deformacion la tension incial y final en un punto material ha de
ser la misma.
iii. La elasticidad no esta asocidad a la linealidad de la respuesta.
materiales que se comportan elasticamente pero de forma no lineal.

Existen

iv. La definicion del concepto de elasticidad es local por lo que permite considerar
cuerpos en los que una region se comporte linealmente y otra no.

2. El principio de superposici
on
Combinando las ecuaciones del equilibrio, la definicion de la deformacion y una
relacion constitutiva de la forma T = f (D) se define de forma completa un problema
de contorno que, bajo ciertas hipotesis en el modelo constitutivo, se puede demostrar
que tiene solucion y u
nica ademas. Sin embargo, que esta solucion exista no significa
que se pueda calcular en la practica pues el sistema resultante de ecuaciones en
derivadas parciales es bastante complejo.
Para facilitar la resolucion de problemas elasticos se adopta muy a menudo
una hipotesis de calculo que de hecho se verifica de forma muy precisa cuando

68

las deformaciones son peque

nas. Para describir esta hipotesis consideramos que
la solicitacion sobre un cuerpo se compone de las fuerzas de volumen sobre y de
superficie sobre t y y esto lo indicamos con el par (fv , fs ). Esta solicitacion produce
en el cuerpo una repuesta que consta de un desplazamiento, una deformacion y
una tension que escribimos como una terna (u, D, T ). Los dos primeros campos
estan relacionados por la definicion de la deformacion, los dos u
ltimos por la ley
constitutiva y la tension ademas ha de estar en equilibrio con la solicitacion.

Definici
on 2.3: Se dice que un cuerpo satisface el principio de superposici
on
cuando la respuesta a una solicitaci
on combinada (fv 1 +fv 2 , fs 1 +fs 2 ) es la suma
1
1
1
2
de las respuestas (u , D , T ) y (u , D 2 , T 2 ), siendo cada una de estas la debida
a la acci
on independiente de la primera y segunda solitaci
on, respectivamente.
Como las ecuaciones del equilibrio son lineales y tambien lo es la relacion entre
el desplazamiento y la deformacion, el principio de superposicion es equivalente a
suponer que el modelo constitutivo es tambien lineal. La relacion T = f (D) puede
ser por tanto descrita por una funcion f tensorial de variable tensorial y lineal. Ese
tipo de funciones se conocen como tensores de cuarto orden.
En general, la definicion de una relacion lineal entre D y T , o dicho de
otra manera un tensor de cuarto orden entre tensores simetricos, require la
definicion de 6 6 = 36 constantes elasticas. Sin embargo en este curso solo
estudiaremos materiales isotropos y para este caso solo son necesarias dos constantes
independientes.
El principio de superposicion se supondra siempre cierto a partir de ahora,
como la hipotesis de peque
nas deformaciones aunque hay que subrayar que no es
rigurosamente cierto.

3. Las constantes el
asticas de un material is
otropo
Cuando se estudia la respuesa elastica de un material isotropo bastan dos
constantes para definirla completamente. Para identificar estas dos propiedades,
puesto que son locales, se establecen ensayos mecanicos muy sencillos en los que los
estados de tension y deformacion sean homogeneos, a partir de los cuales se pueden
sacar conclusiones de la respuesta puntual.
Mencionamos otra vez que existen infinitas parejas de constantes que definen de
forma unvoca el comportamiento elastico de un material isotropo. Escogemos el

69

modulo de Young y el coeficiente de Poisson, definidos a continuacion, por ser las

dos constantes que tienen un significado fsico mas claro.
Para describir matematicamente la respuesta elastica de un material isotropo
se estudia el estado de tension homogenea mas sencillo que existe, esto es, el de
traccion uniaxial. Para ello se toma una barra de material sujeta por un extremo
y se somete a una fuerza de traccion por el extremo libre. Ignorando los efectos
cercanos a la sustentacion, el estado de tension y deformacion de todos los puntos
de esta barra es

x 0 0
x 0 0
(4.1)
[D] = 0 y 0 ,
[T ] = 0 0 0 ,
0 0 z
0 0 0
siendo (x, y, z) un sistema de coordenadas con el eje x paralelo al eje del cilindro. En
la barra se observa que la dimension axial aumenta, que las dimensiones trasversales
disminuyen y que los angulos relativos no cambian.

3.1. El m
odulo de Young
El modulo de Young E es la pendiente de la curva (, ) en un ensayo uniaxial
de traccion/compresion. Si denominamos x el eje paralelo a la direccion del ensayo,
se verifica
x
x =
.
(4.2)
E
Las dimensiones de E son F/L2 . El acero, por ejemplo, posee un modulo de Young
de valor Eacero = 210 GPa.

3.2. El coeficiente de Poisson

El coeficiente de Poisson relaciona la deformacion axial con el acortamiento
trasversal en un ensayo uniaxial de traccion. Se define este escalar mediante la
relacion
y
z
= = ,
x
x
siendo valida la segunda igualdad debido a la hipotesis de isotropa. Siendo un
cociente entre deformaciones, el coeficiente de Poisson es adimensional. En el acero,
por ejemplo, este coeficiente tiene el valor acero = 0.3.

70

4. La ley de Hooke generalizada

En la seccion anterior se han definido el modulo de Young y el coeficiente de
Poisson, como constantes elasticas que describen completamente la deformacion en
un ensayo de traccion uniaxial (o compresion). A continuacion demostramos que
a partir de estos datos podemos encontrar la relacion tension-deformacion en un
punto sometido a un estado tensional arbitrario, encontrando la llamada ley de
Hooke generalizda, pues extiende al continuo la relacion elastica de los resortes.
Para obtener dicha relacion partimos primero de la siguiente observacion: el
estado tensional mas complejo que puede ejercerse sobre un diferencial de volumen
es un estado triaxial de traccion/compresion. Efectivamente, cualquier estado
tensional puede expresarse, en la base principal de tension, como un estado triaxial
de traccion/compresion. En dicho estado las tres tensiones normales se denominan
1 , 2 , 3 y coinciden con las tensiones principales.
Aceptando como valido el principio de superposicion, estudiamos la respuesta
al estado triaxial de tension superponiendo tres estados de traccion/compresion
uniaxial. Comenzando por la traccion/compresion sobre un plano perpendicular a
la direccion principal primera, el estado de tension y deformacion correspondiente
es:

1 0 0
0
0
E
[D (1) ] = 0 E1
[T (1) ] = 0 0 0 ,
0 .
(4.3)
0 0 0
0
0
E1
Estudiando a continuacion un estado de traccion/compresion uniaxial en la direccion
principal de tension segunda obtenemos un nuevo estado de deformacion

0 0 0
E2 0
0
2
[T (2) ] = 0 2 0 ,
[D (1) ] = 0
0 .
(4.4)
E
0 0 0
0
0 E2
Finalmente, considerando el tercer estado de tension posible se obtiene que la tension
y deformacion son

0 0 0
E3
0
0
[T (3) ] = 0 0 0 ,
[D (3) ] = 0
E3 0 .
(4.5)
3
0
0
0 0 3
E
Por el principio de superposicion, la deformacion debida a un estado tensional

71

T = T (1) + T (2) + T (3) es la suma D = D (1) + D (2) + D (3) , o en forma de matriz:

1

2
3
0
0
E E E
2
1
3
.
[D] =
0
(4.6)
0
E E E
1
2
3

0
0
E
E
E
La primera conclusion que se obtiene de (4.6) es que, en un material elastico
isotropo, las bases principales de tension y deformacion coinciden. Sobre todo,
esta expresion indica la relacion mas general posible entre tension y deformacion
de un material de estas caractersticas cuando estas dos cantidades se expresan en
componentes de la base principal. Para hallar la expresion intrnseca, valida para
cualquier sistema de coordenadas, no necesariamente cartesiano reformulamos la
anterior expresion de la siguiente manera:

1 0
0
1 0 0

1+
[D] =
0 2 0 (1 + 2 + 3 ) 0 1 0
E
E
(4.7)
0
0 3
0 0 1
1+

=
[T ] tr[T ] [I] .
E
E
Esta u
ltima expresion depende solo de operadores intrnsecos, pues en ning
un lugar
se hace referencia a componentes o sistemas de coordendas, as que se puede formular
de manera completamente general la siguiente ley de Hooke generalizada:
D=

1+

T tr[T ] I .
E
E

(4.8)

Para la resolucion de problemas resulta u

til recoger la expresion en componentes
E
cartesianas de (4.8). Definimos para ello el m
odulo de cortante G = 2(1+)
y
escribimos
xy

x
(y + z ) ,
xy =
,
x =
E
E
G
y

xz
y =
(z + x ) ,
xz =
,
(4.9)
E
E
G
yz
z

z =
(x + y ) ,
yz =
.
E
E
G
En estas expresiones se puede leer un resultado adicional importante: en los
materiales elasticos isotropos las tensiones normales solo producen deformaciones
longitudinales (en todos las tres direcciones debido al efecto Poisson) y las tensiones
cortantes solo produce deformaciones angulares (cada tension cortante produce una
deformacion angular desacoplada del resto).

72

5. Las ecuaciones de Lam

e
La ecuacion (4.8) permite calcular la deformacion D en funcion de la tension T
y en esta seccion invertimos esta expresion para encontrar una formula de la tension
en funcion de la deformacion. Para ello, comenzamos amplicando el operador tr[] a
ambos lados de la igualdad (4.8) resultando en
1+

tr[T ] tr[T ] tr[I]

E

E

1+
3
tr[T ]
=
E
E
1 2
=
tr[T ] .
E

tr[D] =

(4.10)

En el captulo 3 se escogio el smbolo para indicar la traza de la deformacion, as

pues
E
tr[T ] =
.
(4.11)
1 2
Sustituyendo este u
ltimo resultdado en la ecuacion (4.8) obtenemos
D=

1+
E
T
I .
E
E 1 2

(4.12)

Despejando el tensor de tension de esta expresion se obtiene

T =

E
E
D+
I .
1+
(1 + )(1 2)

(4.13)

Para poder escribir esta expresion de forma mas compacta definimos el primer y
segundo coeficiente de Lame
=

E
,
(1 + )(1 2)

E
.
2(1 + )

(4.14)

Ambos coefiecientes de Lame tienen dimensiones de F/L2 , como el modulo de

Young, puesto que son rigideces. Como el segundo coeficiente de Lame es igual
al modulo de cortante G escribimos la expresion (4.13) finalmente
T = 2GD + tr[D]bI .

(4.15)

Esta u
ltima expresion se conoce como la ecuacion de Lame y permite obtener la
tension a partir de la deformacion. Como se trata de una ecuacion intrnseca es

73

valida en cualquier sistema de coordendas. En particular, si se expresan todos los

tensores en coordendas cartesianas se obtiene
x = 2G x + ,

xy = G xy ,

y = 2G y + ,

xz = G xz ,

z = 2G z + ,

yz = G yz ,

(4.16)

= x + y + z .

6. Deformaciones y tensiones proporcionales

En un ensayo uniaxial, una tension de traccion provoca una deformacion en
direccion de las tensiones aplicadas. En general esto no es as y un punto sometido
a un estado tensional T experimenta una deformacion D que no es proporcional a la
tension, es decir, D 6= T , para ning
un escalar . Por ejemplo, un punto sometido a
traccion uniaxial sufre deformaciones en las direcciones perpendiculares a al traccion
aplicada debidas al efecto Poisson. Sin embargo, un punto sometido a un estado
de tension de cortante puro solo experimenta deformacion angular y se comprueba
facilmente que D = (2G)1 T . Pretendemos estudiar a continuacion cuantos casos
existen de solicitaciones que provocan estados de deformacion proporcionales a estos.

Teorema 6.4: En un solido elastico isotropo solo los estados de tension esfericos
y los puramente desviadores causan estados de deformaci
on proporcionales a ellos
mismos. En el primer caso, cuando T es esferico,
D = (3)1 T ,

(4.17)

y en el segundo, cuando T es desviador,

D = (2)1 T .

(4.18)

n: Supongamos que en para un estado tensional T , la deformaci

Demostracio
on

provocada en un punto es tal que D = T . Entonces, por las ecuaciones de Lame,

T = 2T +

tr[T ]
I .
3

(4.19)

tr[T ]
I .
3

(4.20)

Esta ecuacion se puede reescribir como

(1 2 )T =

74

Para que esta ecuacion se cumpla para alg

un escalar solo existen dos posibilidades:
tr
[T ]
o bien T = 3 I, o bien ambos lados de la igualdad se anulan. En el primer caso
la tension es esferica y se cumple
(1 2 )

tr[T ]
tr[T ]
I =
I 1 2 = 3 = (3)1 .
3
3

(4.21)

En el segundo caso la tension es desviadora, tr[T ] = 0, y el parentesis en la

ecuacion (4.20) debe de anularse, para lo cual es necesario que = (2)1 .

7. Termoelasticidad lineal
As como la relacion entre las tensiones y deformaciones se observa
cotidianamente, tambien se aprecia en multitud de situaciones que los campos
de temperatura y tension/deformacion estan acoplados. As pues, si se calienta
un cuerpo este se deforma y a veces aparecen en el tensiones. Mas a
un, en
ciertos materiales se observa que incluso una deformacion elastica produce cambios
de temperatura (el llamado efecto Gough-Joule). Este problema acoplado es en
general muy complejo, pero si el acoplamiento en u
nicamente en un sentido (la
temperatura produce deformaciones pero no viceversa) su formulacion es sencilla.
Debido a su importancia consideramos en esta seccion las relaciones constitutivas
de la termoelasticidad lineal.
Se comprueba experimentalmente que un cuerpo isotropo, homogeneo y libre de
coacciones (u = ), cuando se calienta uniformemente se deforma sin que aparezcan
tensiones. Esta deformacion de origen puramente termico es u
nicamente volumetrica
y proporcional al incremento termico y a un coeficiente de dilatacion termica que
indicamos con el smbolo y con dimensiones de temperatura inversa. Llamando
Dter a las deformaciones termicas se cumple por tanto
Dter = T I ,

(4.22)

siendo T el salto termico respecto a una temperatura en la que no existen

deformaciones termicas. Admitiendo el principio de superposicion, podemos
formular una ley de Hooke generalizada con efectos termicos de la forma
D = Dmec + Dter =

1+

T tr[T ] I + T I .
E
E
75

(4.23)

La deformacion tiene por tanto dos componentes: una mecanica y otra termica. Esta
relacion, valida en cualquier sistema de coordenadas, tiene la siguiente expresion en
componentes cartesianas:
x

(y + z ) + T ,
E
E
y

y =
(z + x ) + T ,
E
E
z

z =
(x + y ) + T ,
E
E

x =

xy
,
G
xz
=
,
G
yz
=
,
G

xy =
xz
yz

(4.24)

La relacion de Hooke (4.23) se puede invertir para obtener las ecuaciones de

Lame con efecto de la temperatura. Como en la seccion 5, para despejar la tension
de las ley de Hooke aplicamos el operador traza a ambos lados de la identidad (4.23)
y obtenemos


1+
tr[D] =
tr[T ]
tr[T ] + T tr[I]
E
E
(4.25)
1 2
=
tr[T ] + 3T .
E
As pues, la traza de la tension es
tr[T ] =

E
3 E

T .
1 2
1 2

(4.26)

Sustituyendo este resultado en (4.23) obtenemos finalmente las ecuaciones de Lame

con efecto de la temperatura:
T = 2GD + I T I ,

(4.27)

siendo la constante

E
.
(4.28)
1 2
El tensor T I se conoce con el nombre de la tension de origen termico, as pues
=

T = Tmec + Tter ,

(4.29)

es decir, que en un solido elastico sometido a deformacion y a cambio de temperatura

las tensiones tiene dos componentes, una de origen puramente mecanico y otro
termico. Notese que un un cuerpo que no se puede deformar libremente, si se somete
a un salto termico, desarrolla tensiones de origen termico no nulas. Estas tensiones
pueden ser muy grandes en elementos de maquinas sometidos a altas temperatura
de funcionamente si estos no se dise
nan cuidadosamente.
76

Ejemplo 7.5: Un cilindro de goma con diametro d = 20 mm y longitud L = 200

mm se aloja en una cavidad cilndrica de diametro D = 20, 05 mm. Sobre el
cilindro se coloca un piston rgido.
1. Calcular el estado tensional y de deformacion en el cilindro si el piston lo
comprime con una fuerza total de 1000 N (notese que el estado de tension es
cilndrico). Indicar la longitud del cilindro deformado.
2. Calcular otra vez el estado tensional y de deformacion en el cilindro si,
manteniendo fijo el piston, el conjunto se calienta 140 o C.
Datos: suponer que las paredes del cilindro y el piston son infinitamente rgidas
y que no ejercen ning
un rozamiento sobre el cilindro. Constantes del material
del cilindro: E = 500 MPa, = 0.48, = 20 106 (o C)1 .

El estado tensional y de deformacion sera, en los dos casos, cilndrico y

homomegeo (debido a la simetra del la geometra y cargas, a la ausencia de
fuerzas volumetricas y de rozamiento). Si escogemos un sistema de coordenadas
cartesiano tal que el eje x coincida con el eje del cilindro, la expresion matricial
del la tension y deformacion en este sistema sera

x 0 0
x 0 0
[D] = 0 e 0
[T ] = 0 q 0 ,
0 0 e
0 0 q
Estudiamos ahora por separado los dos casos de carga:
1) Cuando se comprime el cilindro no se sabe si este contacta con la cavidad.
Suponiendo que no contacta, la u
nica componente no nula del tensor de tensiones
es x y su valor es x = P/A siendo A = 4 202 mm2 el area de la seccion. La
longitud del diametro deformado bajo esta hipotesis sera:
d0 = (1 + e)d = (1 +

x )d = (1 + 3, 06 103 )d = 20, 06 mm .
E

Este resultado contradice la hipotesis de que el cilindro no toca la cavidad,

as que no puede ser cierta. Si el contacto ocurre, entonces la tension q ha
de ser no nula, y la deformacion en direccion radial e es conocida y de valor
e = 0, 05/20 = 2, 5 103 . El resto de componentes de los tensores de tension y
deformacion se calculan a partir de las ecuaciones de Lame o de la ley de Hooke

77

generalizada:
x = P/A = 3, 18 MPa ,
1
(E e + x ) = 0, 53 MPa ,
q=
1
x
2q = 5, 34 103 .
x =
E
La longitud del cilindro deformado es: L0 = (1 + x )L = 198, 93 mm.
2) Suponemos, como en el caso anterior, que al dilatarse el cilindro, este no
toca con las paredes de la cavidad. En este caso, obtenemos en primer lugar
la tension en direccion axial en el cilindro a partir de la condicion de que la
longitud de este no vara:
x =

x
+ T = 0 = x = E T = 1.4 MPa .
E

A partir de esta calculamos el diametro deformado:

d0 = (1 + x ) d = (1

x + T )d = 20, 08 mm ,
E

que es contrario a la hipotesis. Por lo tanto se puede garantizar que habra

contacto entre cilindro y cavidad y que, como en el caso primero, e = 2, 5 103 .
Planteamos las ecuaciones de Hooke para las deformaciones en direccion x y
radial y obtenemos
x
q
2 + T (= 0) ,
E
E
1
x
e=
q
+ T (= 2, 5 103 ) ,
E
E

x =

y resolvermos las dos incognitas q y x que resultan tener valores x = 14, 44

MPa y q = 13, 89 MPa.
4

8. Problemas

4.1 En un cuerpo deformable con un estado tensional T se define la presion media

78

en cada punto como el escalar

1
pm = tr[T ] .
3
Demuestra que cuando el solido es elastico lineal se verifica la ecuacion constitutiva
pm = ,
siendo una constante del material denominada rigidez volumetrica y de valor
= + 23 G y la deformacion volumetrica infinitesimal.

4.2 Un solido elastico lineal con modulo de Young E = 20 GPa y coeficiente de

Poisson = 0, 4 se encuentra sometido a una presion hidrostatica de 100 MPa.
1) Postular la expresion del tensor de tensiones, razonando, sin operar, por que es
valido.
2) Si el vol
umen del solido antes de ser sometido a presion era de 102 m3 , Cual
sera su volumen bajo presion?
3) Justificar por que la respuesta a la pregunta anterior es independiente de la
forma del solido.
4) En el caso de que = 0, 5, Cual sera la deformacion volumetrica?
5) No existen materiales con coeficiente de Poisson mayor que 0,5. Intenta razonar,
a partir del ensayo descrito en este problema que ocurra al someter a presion
uniforme un material que tuviera coeficiente de Poisson 0,6. Es fsicamente
posible?

4.3 Las matrices de componentes en coordenadas esfericas del tensor de tensiones T

y del tensor de deformacion infinitesimal D son:

rr
r /2 r /2
rr r r

/2 .
[T ] = r , [D] = r /2
r
r /2 /2

Obtener la ley de Hooke en coordenadas esfericas.

79

4.4 Un cuerpo elastico deformable de

dimensiones 200 100 50 mm3 se
aloja, como se indica en la figura,
en una cavidad rgida y lisa de
dimensiones 202 100 50 mm3 .
Sobre la cara libre se aplica una
presion exterior de valor p = 500
N/mm2 . Si el cuerpo tiene modulo de
Young E = 104 N/mm2 y coeficiente
de Poisson = 0, 3, hallar todas
las componentes de los tensores de
tension y deformacion expresados en
el sistema de coordenadas (O, x, y, z)
de la figura.

4.5 Un cilindro de acero y radio 2 cm tiene un extremo

pegado a una pared rgida y otro extremo pegado a
una placa tambien considerada rgida. El cilindro se
encuentra en el interior de un cilindro hueco, tambien
pegado a la pared y la placa. Este segundo cilindro
es de bronce, y sus radios exterior e interior tiene
dimensiones de 4 y 3 cm, respectivamente. Ambos
cilindros tiene longitud L = 60 cm.
La placa esta unida a una segunda pared mediante
un resorte elastico de constante K = 5 108 N/m,
inicialmente sin deformar.
Si la camisa y el cilindro sufren un incremento de
temperatura T = 20o C, encontrar las tensiones
y deformaciones en los dos cilindros as como la
elongacion del muelle.
Datos: Ea = 210 GPa, Eb = 100 GPa, a =
1, 2 105 K1 , b = 1, 8 105 K1 .

80

Captulo 5

Estudio del problema

el
astico completo

En los captulos anteriores se han presentado los principales elementos del

modelo del solido deformable, incluyendo un caso particular de modelo constitutivo,
el elastico, que sirve de conexion entre la deformacion y la tension. Con todos estos
ingredientes se puede finalmente plantear el problema completo de un cuerpo (solido,
elastico) deformable. Matematicamente, este se expresa como un problema de
valores de contorno, mediante ecuaciones en derivadas parciales o en forma integral.
En este captulo formulamos este problema y discutimos algunas de sus propiedades
mas importantes.
El estudio completo del problema elastico se puede encontrar en [1], o en el
tratado clasico de [3].

81

fs

fv

u
Figura 5.1:

El cuerpo deformable

1. Enunciado completo del problema

Un cuerpo elastico deformable es un dominio R3 con contorno = u t .
En u el cuerpo esta sujeto, y en t hay unas fuerzas de superficie fs conocidas. Todo
el cuerpo esta sometido a fuerzas volumetricas fv . Si el cuerpo esta en equilibrio, es
isotropo y elastico, y solo se consideran peque
nas deformaciones, el desplazamiento
u : R3 , la deformacion D y la tension T satisfacen el siguiente problema de
valores de contorno:
en ,

div T + fs = 0

sobre t ,

T n = fs
T

=T ,
sobre u ,

u=0
S

D = u,
T = tr[D]I + 2G D .

82

(5.1)

Este sistema de ecuaciones en derivadas parciales es el objeto de la teora de

la elasticidad clasica. De hecho, simplemente reemplazando la u
ltima de estas
ecuaciones por una relacion constitutiva mas compleja se define el problema de
la mecanica de solidos deformables en peque
nas deformaciones.

2. El principio de los trabajos virtuales

El problema de contorno definido por las ecuaciones (5.1) se puede formular
de manera equivalente (con alguna salvedad en la regularidad de las soluciones) en
forma integral. Esta nueva formulacion, conocida con el nombre del principio de
los trabajos virtuales es la base del calculo matricial de estructuras, y de metodos
numericos como, por ejemplo, el metodo de los elementos finitos.
Antes de enunciar este principio, se define un desplazamiento virtual como
un campo u : R3 cualquiera que satisface u = 0 en u . Cualquier campo
vectorial definido de esta manera es un desplazamiento virtual valido y es, por tanto,
completamentamente independiente del desplazamiento real u.

Teorema 2.1: Un cuerpo deformable sujeto en u , sometido a fuerzas

volumetricas fv : R3 y fuerzas de superficie fs : t R3 est
a en equilibrio
bajo el campo de desplazamiento u : R3 si, para cualquier desplazamiento
virtual u se verifica
Z
Z
Z
T : D[u] dV =
fv u dV +
fs u dS ,
(5.2)

siendo D[u] = S [u], y T = f (D) con D = S [u].

n: Partiendo de la ecuaci
Demostracio
on diferencial del equilibrio interno se puede

deducir que, para un desplazamiento virtual cualquiera:

Z
0=
(div T + fv ) u dV .

Integrando por partes esta integral, usando T n = fs sobre t y u = 0 sobre u se

obtiene (5.2).

83

Observaciones 2.2:
i. El principio de los trabajos virtuales se ha demostrado a partir del principio
de Cauchy. Alternativamente, el principio de los trabajos virtuales se podra
postular y, a partir de el, obtener el principio de Cauchy. La situacion es analoga
a la que aparece en la mecanica clasica, en la cual las leyes de Newton se pueden
derivar del principio de Hamilton o viceversa.
ii. El principio de los trabajos virtuales es la formulacion integral (o variacional,
o debil, como tambien se conoce) de la ecuacion diferencial del equilibrio de
fuerzas. El resto de ecuaciones de (5.1) se imponen exactamente en cada punto.
Existen formulaciones variacionales mas generales que las del principio de los
trabajos virtuales, que imponen de forma debil el resto de las ecuaciones del
problema elastico (vease, por ejemplo, [9, 2]).

3. Las ecuaciones de Navier

En las secciones 1 y 2 se ha planteado, de dos maneras alternativas, el problema
de contorno de la mecanica de los cuerpos deformables con peque
nas deformaciones.
En ambas formulaciones, aparecen como incognita los campos de desplazamiento
u, de deformacion D y de tension T . Para la resolucion analtica de algunos
problemas resulta u
til plantear el problema de contorno u
nicamente en funci
on
del campo de desplazamientos. Cuando esto se lleva a cabo para las ecuaciones de la
elasticidad lineal se obtienen unas formulas muy compactas que reciben en nombre
de ecuaciones de Navier.
Para obtener dichas ecuaciones, basta con sustituir la expresion de la tension T
en funcion de la deformacion y esta del desplazamiento u resultando en


div div u 1 + 2G S u + fv = 0 ,


(5.3)
div [u]1 + 2G S u n = fs ,
en t ,
u=0,

en u .

Simplificando la primera de estas ecuaciones mediante las relaciones

div (div (u) 1) = grad (div [u]) ,
div (grad u) = 4u ,
T

div (grad u) = grad (div (u)) ,

84

(5.4)

se demuestra inmediatamente que (5.3) se puede escribir como

( + G) grad (div (u)) + G 4u + fv = 0 ,


div [u]1 + 2G S [u] n = fs ,

en t ,

u=0,

en u .

(5.5)

4. El principio de Saint Venant

La experiencia indica que para el estudio de la solucion a un problema de
un cuerpo deformable los detalles exactos de como estan aplicadas las fuerzas de
superficie no son muy relevantes. Por ejemplo, cuando se realiza un ensayo de
traccion, la forma de las mordazas de la maquina de traccion, aunque no puede ser
completamente aleatoria, no afecta el resultado de los ensayos. Lo mismo ocurre
con las tensiones en el terreno bajo una zapata, o en un piston cuando esta sometido
a las presiones de los gases en el interior de un cilindro.
El principio de Saint Venant establece que los campos de desplazamiento,
deformacion y tension debidos a dos distribuciones de fuerzas de superficie
estaticamente equivalentes son iguales lejos de la zona de aplicacion. Esta definicion
deja sin definir cuan lejos los efectos de los detalles en la aplicacion de las fuerzas
dejan de ser perceptibles, as que resulta un poco imprecisa. Como regla general,
se puede estimar que esta distancia es igual a la dimension caracterstica de la zona
de aplicacion de las cargas. En cualquier caso, su aceptacion es fundamental en
ingeniera y siempre lo daremos como valido.
El principo de Saint Venant data de 1855, aunque con el tiempo se ha demostrado
que no es un principio como tal sino que puede ser demostrado. Parte de la dificultad
en demostrarlo radica en que su definicion, como se comento anteriormente, es algo
imprecisa. La definicion precisa de este principio se suele asociar a [7].

5. Estados planos de tensi

on y deformaci
on
El tratamiento analtico de los problemas de cuerpos deformables es, en general,
muy complicado. Existen dos casos particulares que simplifican mucho la descripcion
matematica del problema y que, ademas, son muy habituales. Estos son los casos
de tension y deformacion plana en los que algunas de las componentes de tensor de
tension o deformacion son nulas en todos los puntos del cuerpo. Como se vera a
continuacion, esto es el resultado de geometras y cargas muy particulares.

85

5.1. Estados de tensi

on plana
Definici
on 5.3: Se dice que un cuerpo se encuentra en un estado plano de
tensi
on cuando existe un sistema de
de tensiones en todo punto del cuerpo

1 (x1 , x2 )

[T ] = 21 (x1 , x2 )
0

coordenadas (x1 , x2 , x3 ) tal que el tensor

tiene la expresi
on

12 (x1 , x2 ) 0
2 (x1 , x2 ) 0 .
(5.6)
0
0

Este estado de tension aparece, de forma muy aproximada, en cuerpos planos,

muy delgados con cargas de superficie y volumen contenidas en dicho plano. Notese
que en toda esta seccion la notacion 1 , 2 , 3 se refiere a las tensiones normales a
los planos coordenados, y no a las tensiones principales.
El tensor de deformacion en estados planos de tension tiene por expresion
1 2

12
(x1 , x2 )
(x1 , x2 )
0
E
G
1+

2 1
.
[D] = [
T tr[T ]1] = G21 (x1 , x2 )
(x1 , x2 )
0
E
E
E

0
0
E (1 + 2 )(x1 , x2 )
(5.7)
Notese que la deformacion 3 no se anula.

5.2. Estados de deformaci

on plana
Definici
on 5.4: Se dice que un cuerpo se encuentra en un estado plano de
deformaci
on cuando existe un sistema de coordenadas (x1 , x2 , x3 ) tal que el
tensor de deformaci
on en todo punto del cuerpo tiene la expresi
on

1 (x1 , x2 ) 212 (x1 , x2 ) 0

[D] = 221 (x1 , x2 ) 2 (x1 , x2 ) 0 .
(5.8)
0
0
0
Este estado de deformacion aparece, de forma muy aproximada, en cuerpos
con simetra axial y cargas ortogonales a dicho eje de simetra, que ha de coincidir
con el eje x3 del sistema de referencia indicado anteriormente. Los cuerpos que se
86

encuentran en un estado plano de deformacion tiene un campo de desplazamientos

que, empleando el sistema de referencia cartesiano que se menciona, verifica
u = u(x1 , x2 ) ,

u e3 = 0 .

(5.9)

En estos estados de deformacion, el tensor de tensiones tiene por expresion matricial

en la base {e1 , e2 , e3 }
[T ] = [2G D + tr[D]1]

2G1 (x1 , x2 ) + (x1 , x2 )

G21 (x1 , x2 )
=
0

0
.
0
(x1 , x2 )
(5.10)
Notese que, en general, la componente 3 no se anula. De hecho, podemos escribir
G12 (x1 , x2 )
2G2 (x1 , x2 ) + (x1 , x2 )
0

1 + 2 = 2G(1 + 2 ) + 2 = 2( + G) .
Como + G =

2 ,

(5.11)

la tension en la direccion x3 se puede expresar tambien como

3 = = 2( + G) = (1 + 2 ) .

(5.12)

5.3. El diagrama de Mohr en estados planos

En un estado plano, la direccion que hemos denominado x3 es principal y
la tension asociada es una tension principal (que se anula en el caso de tension
plana). En el plano x1 x2 , existen dos tensiones principales que llamamos I , II con
sus direcciones principales correspondientes vI , vII . Notese que no se correspoden
necesariamente con las tensiones principales 1 , 2 porque puede ser que 3 sea la
mayor de las tres tensiones principales o la intermedia.
Para continuar, y por simplificar la notacion, supongamos que el sistema
coordenado x1 , x2 , x3 es el cartesiano x, y, z. Entonces, en cualquiera de los dos
tipos de estados planos las tensiones principales I y II son las races del polinomio
caracterstico





xy
1 0 
2
0 =  x

= 2 (x + y ) + x y xy
.
(5.13)
xy y
0 1 

87

II

I +II
2

Figura 5.2:

Diagrama de Mohr para estados planos

Estas
se pueden escribir de forma explcita como
s

x y 2
x + y
2 .

+ xy
=
2
2

(5.14)

Consideremos ahora las componentes intrnsecas de la tension t = T n sobre

planos de normal n contenida en el plano xy, es decir, tal que n k = 0. En la base
principal B P = {vI , vII , k} este vector se puede escribir como n = cos vI +sen vII ,
as pues la tension normal sobre dicho plano es

I 0
0 cos
cos
n = n T n = sen 0 II 0 sen

0
0
0 z
0
(5.15)
= I cos2 + II sen2
I + II
I II
=
+
cos(2) .
2
2
Para calcular la componente tangencial definimos el vector unitario m = nk. Este
vector define, s
olo para problemas planos la u
nica direccion tangencial posible sobre

88

el plano de normal n donde puede haber tension tangencial. Este vector ademas
tiene expresion en la base principal m = sen vI cos vII as que podemos definir
la tension tangencial m como la proyeccion t m y calcularla explcitamente de la
siguiente manera

I 0
0 cos
sen
m = m T n = cos 0 II 0 sen

0
0
0 z
0
(5.16)
= I sen cos II sen cos
I II
=
sen(2) .
2
A partir de las expresiones (5.15) y (5.16) se interpreta que las componentes
intrnsecas (n , m ) de la tension en estados planos recorren un circunferencia en
el plano como se indica en la Figura 5.2. A diferencia del diagrama de Mohr para
estados de tension tridimensionales, en el caso plano tiene sentido representar un
crculo completo, puesto que en este caso la tensiones tangenciales m s que pueden
ser negativas.

89

6. Problemas

B
45

5.1 Se colocan tres galgas extensometricas sobre la

superficie de un cuerpo deformable como se indica
en la figura. Si las galgas miden:
3

A = 10

B = 2 10

C = 3 10

y se sabe que el solido esta en un estado de tension

plana, siendo z el eje de tension nula. Calcular el
tensor de deformacion completo en el punto en el
que las galgas realizan las mediciones (E = 20000
Kp/mm2 y = 0, 35).

y
x

5.2 Dados los estados tensionales I y II correspondientes a estados planos de tension,

1) Considerar el estado III que resulta de sumar las tensiones que crean los estados
I y II. Dibujar el diagrama de Mohr correspodiente a este tercer estado.
2) Determinar de forma grafica el valor de para que el estado III sea un estado
de cortante puro.
3) Determinar de forma grafica el valor mnimo de para que en el estado III no
haya compresion en ning
un plano.
4) Determinar de forma grafica el valor maximo de para que en el estado III no
haya traccion en ning
un plano.
5) Resuelve analticamente las tres preguntas anteriores.

90

(NOTA: las tensiones en el estado I estan expresadas en MPa).

4

45

3
4

4
3
3

4
Estado I

Estado II

5.3 Un punto de un cuerpo tiene un estado

tensional plano cuya representacion grafica se
adjunta.

5
A

1) Dibuja el diagrama de Mohr del estado

tensional.
2) Identifica sobre la circunferencia de Mohr el
estado tensional de las caras A y B.
3) Calcula a partir de la figura el valor de las
tensiones principales.
4) Indica el angulo (y el sentido) que forma la
normal nA con la direccion principal primera.
91

B
1
(Tensiones en MPa)

5.4 La figura indica el estado tensional plano de

un punto en un cuerpo deformable.

1. Halla el valor de la tension normal sabiendo

que la tension cortante maxima en ese punto
es de 5 MPa.
2. Dibuja el crculo de Mohr correspondiente al
estado de tension resultante.

I
II

3. Identifica, sobre el crculo, el estado tensional

de la cara I y de la cara II.

5.5 Calcula, a partir del diagrama de Mohr, el

valor de las tensiones normales y tangenciales
sobre cada una de las caras del triangulo equilatero
de la figura, sabiendo que esta en equilibrio y que
las tensiones de la figura estan expresadas en MPa.

92

(Tensiones en MPa)

Captulo 6

La energa el
astica

En el captulo 5 se cierra la definicion de las ecuaciones que describen el problema

de los cuerpos elasticos. En este captulo no se a
nade por tanto ning
un elemento
que no se encuentre contenido, de alguna manera, en las ecuaciones del captulo
anterior. Sin embargo, la identificacion de la expresion de la energa abre la puerta
a teoremas extremadamente u
tiles para la resolucion de problemas en Resistencia
de Materiales y Calculo de Estructuras.

1. El trabajo de las fuerzas exteriores sobre un cuerpo

deformable
El concepto de trabajo es central en mecanica. Para un cuerpo deformable
cualquiera, es mas facil definir el trabajo realizado durante un incremento
(diferencial) del campo de desplazamientos. Como existen dos tipos de fuerzas que

93

se pueden aplicar sobre un cuerpo deformable, a saber, volumetricas y de superficie,

este incremento de trabajo es
Z
Z
fs du dS .
(6.1)
dWext =
fv du dV +

Si lo que se desea es calcular el trabajo total realizdao durante un proceso de

carga completo sobre un cuerpo deformable la definicion debe de integrar (6.1).

Definici
on 1.1: Sea un cuerpo deformable sometido a un proceso de carga fuerzas
externas volumetricas y de superficie (fv (t), fs (t)) definidas en el intervalo
(t) entonces
t [0, T ]. Si en cada instante el campo de desplazamientos es u
el trabajo de las fuerzas externas durante dicha deformaci
on se define mediante
la integral:

Z T Z
Z
0
0

dV +
dA dt .
Wext =
fv u
fs u
(6.2)
0

0 (t) =
siendo u

d
(t).
dt u

Esta definicion es completamente general y revela que la definicion de trabajo

depende del proceso de carga, no solo de los valores iniciales y finales de fv y fs .
Esto complica enormemente el calculo del trabajo, porque require la determinacion
de la solucion u en cada instante de carga. Este es un resultado general que se
estudia con detalle en termodinamica.
Para solidos deformables elasticos lineales el siguiente resultado simplifica la
definicion general del trabajo:

Teorema 1.2: (T. de Clapeyron) El trabajo realizado por un sistema de cargas

(fv (t), fs (t)) definidas en el intervalo t [0, T ], con (fv (0), fs (0)) = (0, 0), sobre un
cuerpo el
astico lineal es
Z
Z
Wext = 12
fv u dV + 12
fs u dA ,
(6.3)

(T ) el campo de desplazamiento al final del

siendo fv = fv (T ), fs = fs (T ) y u = u
proceso de carga.

94

La importancia de este resultado es que, para calcular el trabajo total del

proceso, solo es necesario encontrar el desplazamiento en un estado de carga, el
final. Ademas, revela que el trabajo externo de cualquier proceso de carga solo
depende del valor final del las fuerzas volumetricas y de superficie.
n: Para demostrar este resultado supondremos que fv (t) = p(t)fv y
Demostracio

fs (t) = p(t)fs , siendo p(t) una funcion escalar que satisface p(0) = 0 y p(T ) = 1.
Esta simplificacion supone que las tres componentes de las fuerzas volumetricas y las
tres de las fuerzas de superficie varan de la misma manera. Esta hipotesis simplifica
la demostracion pero el mismo argumento que sigue se podra usar con seis factores
de proporcionalidad independientes.
Si u es el campo de desplazamientos cuando sobre el cuerpo act
ua el sistema de
(t) = p(t)u, debido a la linealidad del problema.
fuerzas (fv , fs ), en cada instante u
Utilizando la definicion 6.2 del trabajo se sigue que

Z T Z
Z
0
0
Wext =
p(t) p (t)fv u dV +
p(t) p (t)fs u dA dt
0

Z

fv u dV +

fs u dA
t

Z
Z

fs u dA
t

Z
fv u dV +

p(t) p0 (t) dt

Z
fv u dV +

Z

h

p2 (t) iT
2 0

(6.4)


fs u dA 12 ,

2. La energa el
astica de deformaci
on
De la misma manera que un resorte elastico, al deformarse, almacena energa,
tambien los cuerpos deformables lo hacen. Esta energa se llama la energa el
astica
de deformaci
on. Un sistema elastico es, por definicion, un sistema deformable
en el que todo el trabajo que realizamos para deformar el cuerpo se recupera al
descargar el sistema de las fuerzas aplicadas. Para calcular la energa elastica
de deformacion en un estado cualquiera bastara por tanto con calcular el trabajo
empleado en deformar el cuerpo hasta dicho punto.
Esta definicion no es valida para sistemas mecanicos mas generales, en los
que se pueda dar disipacion energetica por mecanismos internos (plasticidad,

95

viscoelasticidad, da
no, ...). En el caso que nos ocupa, calculamos el trabajo realizado
por todas las fuerzas exteriores aplicadas sobre un cuerpo deformable que es:
Z
Z
fs u dS .
(6.5)
Wext = 12
fv u dV + 12

Por el principio de los trabajos virtuales, si escogemos el desplazamiento virtual u

igual al desplazamiento real u obtenemos
Z
Z
Z
fs u dS =
T : D dV ,
(6.6)
fv u dV +

Comparando (6.5) y (6.6) se concluye que

Z
1
Wext = 2
T : D dV .

(6.7)

La u
ltima expresion permite calcular la energa en un estado de deformacion
empleando u
nicamente variables del interior del cuerpo y se conoce como Wint la
energa elastica almacenada.
La densidad de energa elastica almacenada se puede escribir de varias maneras.
Si expresamos todo en funcion de la deformacion, obtendremos
T : D = ( tr[D]1 + 2 G D) : D = tr[D]2 + 2 G D : D .
Cuando empleamos esta expresion definimos la funcion
Z


2
1
Wint (D) =

tr[D]
+
2
G
D
:
D
dV ,
2

(6.8)

(6.9)

y a esta funcion, con estos argmentos, se la conoce como la energa interna o la

energa de deformacion.
Cuando la densidad de energa elastica almacenada se escribe en funcion de la
tension,


1+
1+

T :D=T :
T tr[T ]1 =
T : T tr[T ]2
E
E
E
E
se puede expresar toda la energa elastica de la forma


Z
1+

1
Wint
(T ) =
T
:
T

tr[T
]
2
E
E

96

(6.10)

y se conoce como la energa interna complementaria.

Los valores numericos de la energa interna y de la energa interna
complementaria (y tambien el del trabajo de las fuerzas externas) coindicen cuando
la tension y la deformacion estan relacionadas (mediante la ecuacion de Lame). Sin
embargo, las expresiones funcionales son muy distintas y sus propiedades tambien,
como veremos. Ademas, en casos no elasticos, la energa elastica complementaria
tambien puede definirse de forma consistente, y en ellos los valores de ambas energas
no coinciden.

3. El teorema de reciprocidad de Maxwell-Betti

El teorema de reciprocidad de Maxwell-Betti es un resultado muy u
til de
aplicacion en problemas elasticos u
nicamente. Establece lo siguiente:

Teorema 3.3: Sea un cuerpo elastico sometido a dos sistemas de fuerzas que
denominamos (fs 1 , fv 1 ) y (fs 2 , fv 2 ). Si sus desplazamientos, deformaciones y
tensiones son, respectivamente, (u1 , D 1 , T 1 ) y (u2 , D 2 , T 2 ), se verifica
Z
Z
Z
Z
1
2
1
2
2
1
1
1
1
1
fv u dV + 2
fs u dS = 2
fv u dV + 2
fs 2 u1 dS . (6.11)
2

Expresado de otra manera, el teorema de Maxwell-Betti establece que el trabajo

realizado por el sistema de fuerzas (fs 1 , fv 1 ) sobre el campo de desplazamientos u2
es el mismo que el trabajo realizado por el sistema de fuerzas (fs 2 , fv 2 ) sobre el
campo de desplazmientos u1 .
n: La demostraci
Demostracio
on de este resultado es sencilla a partir del teorema

de los trabajos virtuales. Para presentarla, nos basamos en los materiales elasticos
isotropos, los u
nicos que hemos estudiado, sabiendo que el resultado es valido en
condiciones mas generales. Para estos materiales se cumple:


T 1 : D 2 = tr[D 1 ]1 + 2 G D 1 : D 2
= tr[D 1 ]tr[D 2 ] + 2 G D 1 : D 2
= D 1 : ( tr[D 2 ]1) + D 1 : (2 G D 2 )


= D 1 : tr[D 2 ]1 + 2 G D 2
= D1 : T 2 .

97

(6.12)

Para demostrar el teorema de Maxwell-Betti, escribimos el principio de los trabajos

virtuales para el primer sistema de fuerzas
Z
Z
Z
fs 1 u dS .
(6.13)
T 1 : D dV =
fv 1 u dV +

En el principio de los trabajos virtuales podemos escoger como desplazamiento

virtual el campo de desplazmiento real correspondiente al sistema segundo de
fuerzas, con lo que obtenemos
Z
Z
Z
1
2
1
2
fs 1 u2 dS .
(6.14)
T : D dV =
fv u dV +

Podemos repetir este argumento cambiando los papeles de los sistemas de fuerzas y
obtendremos:
Z
Z
Z
2
1
2
1
T : D dV =
fv u dV +
fs 2 u1 dS .
(6.15)

Comparando estas dos ecuaciones y empleando el resultado (6.12) queda demostrado

el teorema.

4. Introducci
on a los modelos estructurales
Los modelos estruturales son simplificaciones del modelo del solido deformable
que sirven para estudiar la mecanica de solidos deformables con geometras sencillas.
En todos estos modelos (la barra, la viga, la placa, etc) el concepto de punto material
desaparece y se estudian de forma conjunta grupos de puntos como por ejemplo
los contenidos en una misma seccion recta. Esta simplificacion introduce errores
y aproximaciones, pero permite resolver analticamente muchsimos problemas
sencillos, pero de gran utilidad practica.
Cada tipo de modelo estructural, como en el estudiado en el captulo 2, define
unas ecuaciones de equilibrio, unos tipos de fuerzas admisibles, una cinematica y
unas relaciones constitutivas.

5. Los coeficientes de influencia

Consideremos un modelo estructural sobre el cual seleccionamos un n
umero
finito N de puntos denominados P1 , P2 , . . . , PN . Sobre cada uno de estos puntos
98

se definen vectores unitarios i , i = 1, 2, . . . , N , y tambien sobre ellos fuerzas y

momentos concentrados Fi , Mj de forma que
Fi = Fi i ,

y Mj = Mj j ,

con 1 i, j N, i 6= j .

(6.16)

Si se quiere estudiar un sistema de fuerzas con dos o mas de ellas aplicadas sobre
un mismo punto, las definiciones anteriores siguen siendo validas si se consideran
dos puntos superpuestos.

Definici
on 5.4: Se llama desplazamiento efectivo o eficaz i en un estado
determinado de carga y deformaci
on a la proyecci
on del desplazamiento del
punto Pi seg
un la direcci
on (y sentido) del vector unitario i :
i = u(Pi ) i .

(6.17)

El desplazamiento efectivo tiene por tanto dimensiones de longitud y es la

componente del desplazamiento del punto i conjugada con la fuerza Fi , es decir,
la parte del desplazamiento que produce trabajo con la fuerza Fi . El trabajo por
tanto efectuado por la fuerza Fi es
Wi = 12 Fi u(Pi ) = 21 Fi i u(Pi ) = 12 Fi i .

(6.18)

As pues, la definicion de desplazamientos efectivos permite ignorar el uso de

vectores, en lo que respecta al calculo de trabajos y energas.
En el caso general, cuando act
uan N fuerzas, el trabajo total de las fuerzas
externas se puede escribir como
Wext =

N
X

1
2 Fi i

(6.19)

i=1

Si sobre un punto Pi existe un momento concentrado aplicado en lugar de una

fuerza, entonces se llama desplazamiento efectivo i en un estado determinado
de carga a la proyeccion del giro del punto Pi seg
un la direccion (y sentido) del
vector unitario i .
i = (Pi ) i .
(6.20)
El desplazamiento efectivo tiene en este segundo caso dimensiones de angulo y
es el la componente del giro en Pi conjugada al par Mi , es decir, la parte del giro
que produce trabajo con la fuerza Mi .
99

Si consideramos el modelo anterior, cuando u

nicamente se aplica una fuerza
Fi sobre el punto Pi (todas las demas fuerzas son cero), el punto Pj sufre un
desplazamiento efectivo que denominamos ji (con dimensiones de longitud o giro).
De la misma manera, si estudiamos la situacion ahora en la que la u
nica fuerza
sobre el sistema sea la fuerza Fj , el punto Pi sufre un desplazamiento efectivo que
indicamos como ij . El teorema de reciprocidad de Maxwell-Betti establece
que si el modelo estructural es elastico y lineal entonces se verifica la relacion
1
2 Fi

ij = 12 Fj ji ,

(6.21)

es decir, el trabajo realizado por la fuerza generalizada Fi sobre el desplazamiento

generalizado debido u
nicamente a Fj es igual a su recproco.

Definici
on 5.5: El coeficiente de influencia ij es el valor del desplazamiento
efectivo en el punto Pi cuando la fuerza o momento aplicado sobre el punto Pj
es unitario. Por la linealidad de la respuesta del sistema, los coeficientes de
influencia se calculan como
ij
ij =
(6.22)
Fj
Las dimensiones del los coeficientes de influencia son:
i) longitud/fuerza, si Fi , Fj son ambos fuerzas puntuales.
ii) longitud/(fuerzalongitud) si Fi es una fuerza y Fj es un momento.
iii) angulo/fuerza, si Fi es un par y Fj una fuerza.
iv) angulo/(fuerzalongitud), si Fi , Fj son ambos momentos puntuales.
Por el principio de superposicion, el desplazamiento efectivo total i en el punto
Pi cuando act
uan las fuerzas (generalizadas) F1 , F2 , . . . , FN es
i = i1 + i2 + . . . = i1 F1 + i2 F2 + . . . =

N
X

ik Fk

(6.23)

k=1
se puede expresar en funci
La energa interna complementaria Wint
on de los
coeficientes de influencia y de las fuerzas proyectadas, pues en un modelo estructural
=W
elastico Wint
ext y se puede escribir:

Wint
=

1
2

N
X

Fi i =

i=1

1
2

N
X
i,j=1

100

ij Fi Fj .

(6.24)

Esta u
ltima expresion es una forma cuadratica que se puede indicar con matrices
de la siguiente manera:

11 12 . . . 1N F1
F1

F2 21 22 . . . 2N F2

1
1

Wint = 2

..
.
.
.
.
..
..
..
.. .. = 2 { F } [ ] { F } . (6.25)

.
.

FN
FN
N 1 N 2 . . . N N
Como la forma cuadratica es definida positiva, los terminos de la diagonal de la
matriz de coeficientes de influencia o matrix de flexibilidad han de ser
positivos, es decir, en toda estructura ii > 0. Este resultado tiene implicaciones
fsicas: cuando una estructura se solicita u
nicamente con una fuerza (generalizada),
el desplazamiento (generalizado) efectivo del punto de aplicacion ha de ser positivo.
Si definimos la matriz de rigidez k como la inversa de la matriz de coeficientes
de influencia, la energa interna de deformacion se puede expresar tambien de forma
matricial:

k11 k12 . . . k1N 1

1

2 k21 k22 . . . k2N 2

1

Wint = 2
..
= 12 { } [ k ] { } .
..
..
..
..
..

.
.
.
.

N
N
kN 1 kN 2 . . . kN N
(6.26)
La matriz k, por la inversa de una matriz simetrica, definida positiva, tambien
tiene estas dos propiedades.

6. Los teoremas de Castigliano

Los dos teoremas de Castigliano estan basados en las expresion de la energa
interna y de la energa interna complementaria para modelos estructurales. En
resistencia de materiales y en calculo de estructuras se emplean constantemente,
sobre todo el segundo teorema, para resolver problemas hiperestaticos. Su enunciado
es el siguiente:

Teorema 6.6: (Primer teorema de Castigliano) En un modelo estructural

en equilibrio, la derivada parcial de la energa Wint (1 , 2 , . . .) respecto de un
desplazamiento efectivo i es igual al m
odulo de la fuerza generlizada aplicada sobre
el punto i:
Wint (1 , 2 , . . .)
Fi =
.
(6.27)
i
101

Y recprocamente:

Teorema 6.7: (Segundo teorema de Castigliano) En un modelo estructural

(F , F , . . .)
en equilibrio, la derivada parcial de la energa complementaria Wint
1
2
respecto del m
odulo de una fuerza generalizda Fi es igual al desplazamiento efectivo
conjugado con dicha fuerza:
i =

(F , F , . . .)
Wint
1
2
.
Fi

(6.28)

n: Las demostraciones de estos dos teoremas se obienen f

Demostracio
acilmente a

partir de las expresiones de la energa interna y la energa interna complementaria en

funcion de los coeficientes de rigidez y de influencia. En el primer caso, calculamos

Wint (1 , 2 , . . .) =
i
i

N
X

1
2

kab a b =

a,b=1

N
X

kib b ,

(6.29)

b=1

que coincide con la expresion del modulo de la fuerza generalizda Fa , lo que

demuestra el primer teorema. De forma analoga, derivando la expresion de la energa
complementaria se obtiene

Wint
(F1 , F2 , . . .) =
Fi
Fi

1
2

N
X
a,b=1

ab Fa Fb =

N
X

ib Fb ,

(6.30)

b=1

que coincide ahora con el desplazamiento efectivo i , demostrando el segundo

teorema.
En la aplicacion de estos teoremas se comprueba que la energa interna y su
complementaria, aunque de valor identico en el caso elastico, no son intercambiables.
Mas concretamente, si utilizamos la energa interna Wint en el segundo teorema
de Castigliano, por ejemplo, en vez de la energa complementaria obtendramos el
resultado
Wint (1 , 2 , . . .)
i =
=0,
(6.31)
Fi
que es absurdo.

102

7. Problemas
6.1 Se propone el siguiente ejemplo para comprobar de forma sencilla la validez del
teorema de reciprocidad. Dos muelles elasticos se colocan en serie tal y como indica
la figura.
1 = 10 N/mm

2 = 20 N/mm
1

Se consideran dos casos de carga:

1) Se aplica una fuerza F1 = 1 N sobre el punto 1 en direccion del eje x y llamamos
11 al desplazamiento del punto 1 y 21 al desplazamiento del punto 2.
2) En segundo lugar se considera el estado de deformacion debido a una fuerza
F2 = 1 N sobre el punto 2, y se llama 21 , 22 a los desplazamientos sobre los
puntos 1 y 2, respectivamente, en este segundo caso.
Comprobar que F1 12 = F2 21 .

6.2 En una estructura se escogen dos puntos A y B sobre los cuales se aplican
fuerzas y se miden desplazamientos efectivos. Sobre la estructura sin deformar se
aplica, en primer lugar, una fuerza FA en el punto A y denominamos estado 1
(E1) al resultante. Sin quitar la fuerza ya aplicada, se a
nade otra fuerza FB sobre
el punto B resultando el estado 2 (E2). La figura muestra la evolucion de los
desplazamientos efectivos en A y B en funcion del valor de las fuerzas aplicadas.
FA (kN)
5

E1

FB (kN)

E2

E2

E1
2

A (m)

103

B (m)

Se pide:
1.- Calcular, empleando el teorema de reciprocidad, el valor de .
2.- Dibujar sobre las misma figuras las curvas fuerzas-desplazamiento que
resultaran de quitar primero la fuerza FA y despues la fuerza FB .

6.3 Una viga empotrada como la de la

figura tiene longitud L. En extremo
2 se aplica una fuerza puntual P y
en el punto medio 1 un momento
concentrado M . El desplazamiento
vertical (en direccion del eje y)
cuando u
nicamente se aplica la fuerza
puntual es:
vP (x) =

P x2
(3L x) ,
6K

L/2

L/2

P
M
x
1

(6.32)

siendo K una constante de rigidez de la viga con dimensiones de N m2 . Cuando

solo se aplica el momento concentrado, el desplazamiento vertical es:

M x2

si x L/2

2K
vM (x) =
(6.33)

M L (2x L/2) si x L/2

4K
Se pide:
1) Calcular 12 , el giro que se produce en el punto 1 cuando se aplica una fuerza
de valor F = 6 N en el punto 2. (NOTA: el giro, cuando los desplazamientos
son peque
nos, se puede aproximar por la tangente, es decir (x) = v 0 (x)).
2) Calcular 21 , el desplazamiento vertical que se produce en el punto 2 cuando se
aplica un par de valor M = 5 Nm en el punto 1.
3) Comprobar el teorema de Maxwell-Betti con los dos resultados anteriores.
4) Calcular los cuatro coeficientes de influencia de la estructura, indicando sus
dimensiones.
104

6.4 Una estructura elastica y plana, de la

cual se dibuja solo una parte, se encuentra
sometida a la accion de dos fuerzas puntuales
PA , PB aplicadas sobre los puntos A y B,
respectivamente, con la direccion y sentido
indicados en la figura.
Cuando se aplica una carga puntual PA =
3 103 N, el punto A se desplaza 3 mm hacia
la izquierda y el punto B, 2 mm hacia la
izquierda y 1 mm hacia abajo. Cuando se
aplica u
nicamente una fuerza PB = 7 103 N,
el punto B tiene un desplazamiento efectivo
cuyo modulo es 4 mm.

PA

A
PB
30o

1) En que direccion y sentido es el desplazamiento del punto B en este u

ltimo
caso de carga? Por que?
2) Encontrar los 4 coeficientes de influencia del sistema, indicando sus dimensiones.
3) Cuanto vale (y en que sentido) el desplazamiento horizontal del punto A cuando
se aplica una fuerza sobre el de valor PA = 2103 N y otra de valor PB = 3103
N sobre el punto B?

6.5 Sobre un punto de una estructura se

aplican dos fuerzas como se indica en la
figura. Se conocen los siguientes coeficientes
de influencia:
11 = 2 102 mm/N ,
22 = 3 102 mm/N ,
12 = 102 mm/N .
1. Encontrar el desplazamientos efectivo en
la direccion de cada una de las fuerzas.
2. Encontrar graficamente la posicion del
punto deformado despues de la aplicacion
de las dos fuerzas (usa la escala que creas
mas conveniente).

105

6.6 Se consideran dos muelles elasticos de constantes k1 , k2 siendo el primero mas

rgido que el segundo. Ambos muelles tiene un extremo sujeto y el otro libre.
1) Si se aplica una fuerza de traccion F igual a ambos muelles, Cual almacena
mas energa elastica?
2) Si el extremo libre de cada uno de los muelles se desplaza una distancia , cual
almacena mas energa elastica?
3) Como se almacena mas energa elastica en un muelle, sometiendolo a una fuerza
de traccion o a una de compresion? (iguales en modulo).
4) Como colocaras los muelles (en serie o en paralelo) para que almacenaran la
mayor cantidad posible de energa cuando se aplica una u
nica fuerza F en el
extremo libre? Y si se aplica un desplazamiento sobre el extremo libre?

6.7 Demuestra que la energa elastica almacenada en un cuerpo homogeneo sometido

a una presion hidrostatica uniforme p es
W =

V 2
p ,
2K

siendo K la rigidez volumetrica del material y V el volumen del cuerpo.

6.8 En la estructura de la figura

se conocen los siguientes coeficientes
de influencia:
11 = 2 10
3

21 = 10

P1
60o

mm/N ,

mm/N .

Calcular el vector de desplazamiento

del punto P1 cuando se aplican
cargas de valor F1 = 1000 N y
F2 = 2000 N, sabiendo que la componente
horizontal de dicho desplazamiento
es ux (P1 ) = 3 mm.

106

F1
x
P2

F2

6.9 Una estructura elastica esta sometida a dos fuerzas puntuales F1 , F2 . Sobre
cada una de los puntos, y en la direccion de la fuerza correspondiente, se definen
desplazamientos efectivos 1 , 2 . La matriz de coeficientes de influencia es:


4
[] =
1

1
1


mm/kN .

1) Indicar la expresion analtica de la energa potencial interna complementaria

U (F1 , F2 ), como una funcion u
nicamente de las fuerzas aplicadas.
2) Encontrar los desplazamiento efectivos 1 , 2 que aparecen cuando se aplican
fuerzas F1 = 3 KN, F2 = 5 KN.
3) Calcular el trabajo Wext realizado por las fuerzas externas indicadas.
4) Comprobar que el trabajo externo es igual a la energa potencial elastica
almacenada en el cuerpo.

6.10 Considera un paraleleppedo de lados Lx , Ly , Lz paralelos a los ejes

coordenados. Si las caras perpendiculares al eje x estan sometidas a una fuerza
por unidad de superficie de valor f , demuestra que la energa de deformacion del
paraleleppedo es
Wint = 21 K 2 ,
siendo = EA/Lx , A = Ly Lz y = f A/K.

107

6.11 En una estructura como la de la

figura act
ua una fuerza puntual y un
momento concentrado, en la direccion y
sentido indicados. Cuando act
ua una
fuerza de valor F = 104 N, se sabe que el
desplazamiento efectivo bajo la fuerza es
de 5mm y el giro efectivo bajo el momento
es de 103 rad, que es en sentido opuesto
al momento de la figura. Cuando se aplica
u
nicamente un par de valor M = 1000
N m se sabe que el giro bajo M es de
102 rad. Determinar el valor de los
cuatro coeficientes de influencia indicando
claramente sus unidades.

108

Captulo 7

Estudio de la finalizaci
on
del comportamiento
el
astico

Hasta ahora en el curso nos hemos limitado a estudiar el comportamiento elastico

de los materiales, obteniendo resultados energeticos y teoremas cuya aplicacion se
limita a dicho rango. Todos los materiales poseen un rango elastico limitado, y
por tanto, para estar seguros de que bajo las solicitaciones consideradas el material
sigue siendo elastico estudiamos en este captulo como detectar cuando se agota
dicho regimen. El objetivo no es, por tanto, el estudio de la ineasticidad, o de la
fractura, que son el objeto de otras materias relacionadas. En este captulo nos
limitaremos a estudiar modelos que nos informen de la validez de las hipotesis de la

109

elasticidad y, mas a
un, que nos proporcionen algo de informacion sobre cuanto de
lejos esta cada punto material del final del rango elastico. Todos estos objetivos se
estudian distinguiendo los materiales d
uctiles de los fragiles porque no hay un u
nico
modelo que se adapte bien a todos los tipos de materiales.

1. Criterios de fallo
Los solidos salen del regimen de comportamiento elastico por motivos muy
distintos, dependiendo de la microestructura de los materiales que los constituyen.
Por ejemplo, los metales dejan de ser elasticos cuando plastifican debido a la
nucleacion de dislocaciones en la red cristalina de cada grano. Los polmeros tambien
salen del regimen elastico, pero en este caso se debe a desenrollamiento de cadenas
polimericas. Por u
ltimo, los materiales ceramicos o el hormigon dejan de ser elasticos
debido a la aparicion de microfisuras. Por unificar conceptos, llamaremos fallo a
la finalizacion del comportamiento elastico de un material, independientemente del
micromecanismo responsable del mismo.
Un criterio de fallo es un modelo matematico que intenta explicar cuando se
inicia el fallo de un punto material a partir del estado de tensiones y/o deformaciones
del mismo. Aunque estan inspirados en la micromecanica de los materiales, los
criterios de fallo son solo formulas sencillas que, con uno o varios parametros, ajustan
los resultados experimentales de la mejor forma posible. No hay ning
un criterio de
fallo exacto para todo estado tensional T .
En este curso estudiaremos criterios de fallo de la forma f (T ) 0. Cuando
la funcion f (T ) es negativa, el punto se encuentra en regimen elastico. Cuando
f (T ) = 0, el criterio predice que se produce el fallo. Lo que ocurre si f > 0 no tiene
interes porque el criterio no proporciona entonces informacion u
til. Cuando el valor
de f (T ) es negativo, su modulo indica, la distancia que esta el punto del fallo.
Aunque no lo definamos con precision, si f (T1 ) < f (T2 ), entonces el estado T1 esta
mas lejos del fallo que el estado T2 .
Por simplificar mas a
un los criterios de fallo, nos basaremos en el ensayo de
traccion/compresion pura para definir los criterios de fallo. En un material d
uctil,
sabemos que el fallo plastico ocurre cuando la tension alcanza el lmite elastico e ;
en cambio, un material fragil falla cuando la tension alcanza el valor r , la tension
de rotura. Si definimos la tension u
ltima u al lmite elastico, si el material es d
uctil,
o la tension de rotura, si el material es fragil, consideraremos en este curso criterios

110

de fallo siempre de la forma:

f (T ) = eq (T ) u ,

(7.1)

siendo eq un escalar que denominamos la tensi

on equivalente y que siempre ha
de definirse de acuerdo a un criterio de fallo.
Para cuantificar la severidad de un estado tensional respecto de un criterio de
fallo, se define el coeficiente de seguridad del estado tensional T respecto al
criterio de fallo f como el escalar n tal que
f (nT ) = 0 .

(7.2)

De acuerdo a las dos definiciones anteriores, la tension equivalente eq (T ) es aquella

tension que en un ensayo de traccion/compresion pura tendra el mismo coeficiente
de seguridad que T .
Un criterio de fallo no puede depender de T de cualquier manera. Para que este
sea fsicamente correcto, por ejemplo, no puede ser una funcion de las componentes
de la matriz asociada a T que dependa del sistema de coordenadas escogido.
Expresado de otra manera, la funcion f solo puede depender de invariantes de
T y si ademas, es isotropa, no puede depender de ninguna direccion. Existen
infinitos invariantes del tensor tension, pero solo se pueden escoger tres que
sean funcionalmente independientes. Tpicamente se escogen, bien los invariantes
principales descritos en el captulo 1, o bien las tres tensiones principales. Por
unificar conceptos utilizaremos siempre estas tres u
ltimas y, abusando de la notacion,
escribiremos:
f (T ) = f (1 , 2 , 3 ) = eq (1 , 2 , 3 ) u .
(7.3)

2. Criterios de fluencia para materiales d

uctiles
Independientemente de los micromecanismos responsables de la finalizacion del
comportamiento elastico en los materiales d
uctiles, estos se caracterizan por una
rama plastica muy larga hasta el fallo definitivo. Por ello, todos los criterios de fallo
de materiales d
uctiles se llaman criterios de fluencia.
Entre los materiales d
uctiles, los mas comunes son los metales. Existen varios
criterios para modelar su fallo y a continuacion describimos los dos mas habituales.

111

2.1. El criterio de Tresca

El criterio de Tresca (1814-1885) se basa en una serie de experimentos llevados
a cabo entre 1864 y 1873 por dicho ingeniero frances. En ellos, Tresca estudio
la deformacion plastica y el punzonamiento de placas y cilindros de plomo, cobre,
parafina, hielo, ... Los informes de estos experimentos fueron, durante 80 a
nos, los
mas completos sobre el tema de plasticidad. En ellos se describen, por primera
vez, el regimen elastico, el endurecimiento plastico y la fluencia de los metales.
Sobre este u
ltimo aspecto, ademas de identificar por vez primera que los metales
fluyen como lquidos, de forma isocorica, demostro que esto ocurre siempre bajo un
estado tensional en el que la tension tangencial maxima tiene un valor caracterstico,
constante para cada material. Como en un ensayo de traccion pura la tension
tangencial maxima toma el valor /2 propuso la siguiente funcion de fluencia:
T resca
fT resca (1 , 2 , 3 ) = eq
(1 , 2 , 3 )e ,

T resca
eq
(1 , 2 , 3 ) = 1 3 .
(7.4)

2.2. El criterio de von Mises

El segundo criterio de fluencia que consideramos fue formulado por Maxwell
hacia 1865, pero se suele atribuir a von Mises (18831953) que lo publico en 1913.
La motivacion para este criterio tambien se encuentra en el comportamiento de los
metales y expresa matematicamente que la plasticidad ocurre cuando la energa de
distorsion alcanza un umbral caracterstico del material.
La energa de distorsion es la energa que tiene la parte desvidora de la tension
definida como s = T + pm 1 con pm = e 13 tr[T ]. En un ensayo de traccion pura,
el valor de esta energa cuando se alcanza el lmite elastico se puede calcular y
es (1 + )/(3E)e2 . Calculando tambien esta energa en funcion de las tensiones
principales se puede establecer la siguiente funcion de fluencia:
vM
fvM (1 , 2 , 3 ) = eq
(1 , 2 , 3 ) e ,

siendo la tension equivalente respecto al criterio de von Mises igual a

q
vM
eq
(1 , 2 , 3 ) = 12 [(1 2 )2 + (2 3 )2 + (3 1 )2 ] .

112

(7.5)

(7.6)

Ejemplo 2.1: Un punto de un cuerpo deformable ductil esta sometido a un estado

tensional cuya matriz asiociada, en un sistema de referencia cartesiano, es

10 10 0
20 0 MPa .
(7.7)
[T ] = 10
0
0
15
Calcular la tension equivalente en el punto seg
un los criteriosd de Tresca y von
Mises. Si se sabe que el lmite elastico del material es e = 80 MPa, calcular
ademas el factor de seguridad del estado tensional anterior seg
un cada uno de
los dos criterios indicados.
Las tensiones principlales de este estado tensional son

1 = 15 + 5 5 MPa ,
2 = 15 MPa ,
3 = 15 5 5 MPa .

(7.8)

Las tensiones equivalentes seg

un los criterios de Tresca y von Mises son:

vM
Tr
eq
= 19.37 MPa .
(7.9)
eq
= 10 5 = 22.36 MPa ,
En cada caso, el factor de seguridad es
nT r =

80
= 3.58 ,
22.36

nvM =

80
= 4.13 .
19.37

(7.10)

3. Criterios de rotura para materiales fr

agiles
Los materiales fragiles fallan de forma s
ubita, sin aparente fluencia, y por ello los
criterios de fallo se denominan criterios de rotura. Ademas, otra caracterstica
que distingue los materiales fragiles de los d
uctiles es su habitual anisotropa pues
resisten mucho mas a compresion que a traccion.

3.1. El criterio de Rankine

El criterio de Rankine predice que un punto material falla cuando, bien la tension
principal mayor 1 alcanza la tension de rotura a traccion rt , o bien la menor tension
principal 3 alcanza la tension de rotura a compresion rc .
113

| |

(0, C)

rc

rt

Figura 7.1:

Representaci
on gr
afica del criterio de Mohr-Coulomb.

Matematicament el criterio de Rankine se puede expresar como

Rankine
fRankine (1 , 2 , 3 ) = eq
(1 , 2 , 3 ) rt ,
Rankine
eq
(1 , 2 , 3 ) = max(1 , rt 3 |rc |)

(7.11)

3.2. El criterio de Mohr-Coulomb

La motivacion para el criterio de Mohr-Coulomb surge de la observacion
experimental que indica la resistencia al cortante de ciertos materiales es sensible
a la presion media. Este tipo de comportamiento se asemeja a la ley de Coulomb
de la friccion y fue Mohr en 1882 quien noto que la condicion de fallo en este caso
coincide con el instante en el que el mayor crculo de Mohr es tangente a una recta,
denominada la recta caracterstica del material.
Como se puede apreciar en la Figura 7.1, la recta caracterstica intersecta el
eje vertical del diagrama de Mohr en el punto (0, C). La constante C del material
indica su resistencia a cortante cuando la tension normal es nula y se llama por
ello la cohesi
on del mismo. El angulo determina cuanto crece la resistencia a la
cortadura en funcion de la tension normal. Por analoga con la ley de Coulomb del
rozamiento, esta constante material se llama el
angulo de fricci
on del material.
En la figura Figura 7.1 se observa que los estado tensionales correspondientes a
los estados de tension y compresion pura en el punto de rotura son tangentes a la

114

recta caracterstica del material. Por tanto se puede escribir:

sen =

rt /2
,
H rt /2

sen =

y tambien

|rc |/2
.
H + |rc |/2

(7.12)

Igualando ambas expresiones del seno del angulo de friccion se obtiene que
H=

rt
,
1k

con k =

rt
.
|rc |

(7.13)

Una vez obtenida el valor de la tension H para la cual el material no resiste ning
un
esfuerzo tangencial, se puede despejar el valor del seno del angulo de friccion como:
sen =

1k
.
1+k

(7.14)

Por u
ltimo, y tambien a partir de la Figura 7.1, se puede escribir que, en
cualquier estado de fallo se ha de verificar:
sen =

1 3
2
3
1 +
2

(7.15)

Y sustituyendo los valores de sen y H obtenidos, respectivamente, en (7.14)

y (7.13) resulta que en cualquier estado de fallo:
1 k3 rt = 0 .

(7.16)

Concluimos que la funcion de fluencia para el criterio de Mohr-Coulomb se puede

escribir como:
MC
fM C (1 , 2 , 3 ) = eq
(1 , 2 , 3 ) rt ,

En el caso en el rt = rc

MC
eq
(1 , 2 , 3 ) = 1 k3 .
(7.17)
el criterio de Mohr-Coulomb coincide con el de Tresca.

3.3. El criterio de Drucker-Prager

De la misma manera que el criterio de Mohr-Coulomb generaliza el criterio de
Tresca introduciendo una dependencia de la resistencia con la presion, Drucker y
Prager en 1950 propuesieron una extension del criterio de von Mises para capturar
el mimo efecto. La funcion de fallo en este caso es de la forma
fDP (1 , 2 , 3 ) = kdev T k p K ,

115

(7.18)

Ejemplo 3.2: Un solido esta sometido a una solicitacion de forma que en un punto
el estado tensional se puede expresar, en una base cartesiana, como

10 10 0
[T ] = 10 15 0 MPa .
0
0
2

(7.19)

La tension de rotura a traccion del material es rt = 10 MPa y la de compresion

es rc = 40 MPa. Calcular la tension equivalente en el punto seg
un los criterios
de Rankine y de Mohr y los factores de seguridad en cada caso. Dibujar el
diagrama de Mohr del estado tensional en el punto y los diagramas de los estados
tensionales cuando la tension es T 0 = nT , siendo n cada uno de los coeficientes
de seguridad previamente calculados.
4

4. Problemas
7.1 Suponiendo que et < |ec |, determinar en funcion de et y de k = et /|ec | el
punto de interseccion de la recta caracterstica de un material con el eje horizontal y
el angulo que forma con el mismo. (Pista: considera un estado esferico de tension
y encuentra el valor de la presion cuando se produce el fallo).
Demostrar ademas que un punto material falla, de acuerdo con el criterio de
Mohr-Coulomb, cuando existe un plano en el que la tension tangencial alcanza el
valor
= C n tan .
(NOTA: C y se llaman, respectivamente, la cohesion y el angulo de friccion del
material)

7.2 Razona las siguientes cuestiones:

1) Un material es tal que, cuando se ensaya mediante un ensayo a traccion pura,
el criterio de Rankine, el de Tresca y el de von Mises predicen el mismo valor
para la tension normal u
ltima. Sin embargo, si el ensayo es de cortante, cada
uno de los modelos predice una tension u
ltima distinta. Como es eso posible?
Cual de todas las predicciones es la correcta?
2) Un solido se somete a un estado hidrostatico de tension T = pI. Cual es la
presion u
ltima seg
un el criterio de von Mises? Comenta la respuesta.

116

3) Sea f (1 , 2 , 3 ) = fb(1 , 2 , 3 ) C una funcion de fluencia. Podemos elegir un

valor de la constante C de forma que el criterio prediga exactamente la tension
de plastificacion en el ensayo de traccion y de forma aproximada la tension de
plastificacion en el ensayo de cortante puro, o viceversa. En general, estos dos
valores de la constante C seran distintos. Cual es mejor? Cual es mas exacto?.
4) Por que un criterio de fluencia de la forma f (T ) = xx yy es inaceptable?
5) Si tuvieras un material especfico con el que construir una pieza mecanica, o una
estructura, Cual criterio de fallo escogeras y por que?

7.3 Un punto de un solido se encuentra en un estado de tension plana descrito por

el esquema de la figura ( > 0).
1. Dibuja el diagrama de Mohr.
2. Encuentra la expresion del tensor de tensiones
en el sistema de coordenadas xyz.
3. Si el material tiene un lmite elastico e ,
encontrar el factor de seguridad de dicho
estado tensional seg
un los criterios de Tresca y
von Mises.
4. Si el lmite de rotura a compresion es rc =
1, 2r , cual es el factor de seguridad seg
un los
criterios de Mohr y de Rankine?

7.4 Un punto de un solido deformable se encuentra sometido a un estado de tension

que, en un sistema de coordenadas dado se expresa como:

30 5 0
[T ] = 5 12 0 MPa
(7.20)
0
0 6
1) Dibuja los crculos de Mohr de este estado tensional.
2) Encuentra la tension equivalente al estado tensional (7.20) seg
un los criterios de
Rankine, Tresca, von Mises y Mohr simplificado (rt /rc = 0, 8).
117

3) Si el lmite elastico a traccion del material es de 50 MPa, calcula los coeficientes

de seguridad del estado tensional seg
un los mismos criterios que en el apartado
anterior.
4) Dibuja los tres crculos de Mohr en el estado tensional n T , siendo n el factor
de seguridad obtenido con el criterio de Tresca. Repite el dibujo siendo n ahora
el factor de seguridad que resulta de emplear el criterio de Mohr simplificado.
Comprueba en este caso que el mayor de los crculos de Mohr es tangente a la
recta caracterstica del material.

7.5 El fallo de un punto material se estudia con el criterio de Mohr. Se sabe que la
cohesion de dicho material es C = 12 MPa y el angulo de friccion es = 30o . Si
las tensiones principales en ese punto son 1 = 5 MPa, 2 = 3 MPa y 3 = 0 MPa,
encontrar graficamente el coeficiente de seguridad en el punto respecto al criterio de
fallo indicado.

7.6 El estado tensional en un punto es

plano y su diagrama de Mohr aparece en el
la figura. Calcular, sabiendo que el lmite
elastico es e = 180 MPa:
1) La tension equivalente seg
un el
criterio de von Mises y el coeficiente
de seguridad asociado.
2) La tension equivalente seg
un el
criterio de Tresca y el coeficiente de
seguridad asociado.

(MPa)

n (MPa)
30

40

7.7 Se desea emplear el criterio de Mohr-Coulomb para estudiar la rotura de un

material fragil. Se sabe que el angulo de friccion de dicho material es = 30o y
que, al someter una probeta a un ensayo de cortante puro, la tension tangencial en
el instante del fallo es r = 10 MPa. Determinar la cohesion del material.

118

Bibliografa

[1] M. E. Gurtin. The linear theory of elasticity. VIa/2, 1972.

[2] K D Hjelmstad. Structural mechanics. Springer Science+Business Media, second
edition, 2005.
[3] A E H Love. A treatise on the mathematical theory of elasticity. 1927.
[4] L E Malvern. Introduction to the Mechanics of a Continuum Medium.
Prentice-Hall, Englewood Cliffs, N.J., 1969.
[5] G E Mase. Continuum Mechanics. McGraw-Hill, 1970.
[6] X Oliver and C Agelet de Saracibar.
ingenieros. Ediciones UPC, 2000.

Mec
anica de medios continuos para

[7] R Stenberg. On Saint-Venants principle. Q J Appl Mech, 11:393402, 1954.

[8] A M Stuart and O Gonzalez. A first course in continuum mechanics. Cambridge
University Press, 2008.
[9] Ky
uichiro Washizu. Variational Methods in Elasticity and Plasticity. Pergamon,
1982.

119