Captulo 4

Ecuaciones constitutivas
4.1. Introduccion: comportamiento uniaxial de
materiales elasticos, plasticos y viscosos
Hasta este punto hemos analizado la descripcion matematica del cambio
de forma o deformacion y de las fuerzas internas o tensiones que se desarro-
llan en un solido deformable como consecuencia de las acciones exteriores.
Este captulo trata de las ecuaciones que caracterizan la respuesta de un
determinado material frente a las cargas exteriores, es decir, de la relacion
entre tensiones y deformaciones. Estas ecuaciones se denominan ecuacio-
nes constitutivas pues describen el comportamiento macroscopico resul-
tante de la constitucion interna de un material. La respuesta real de uidos
y solidos es extremadamente compleja y variada. No solo se ve inuida por
acciones puramente mecanicas, sino tambien por cambios de temperatura o
humedad, e incluso por la historia del material, es decir por situaciones en
las que se ha encontrado con anterioridad. No resulta posible escribir un con-
junto de ecuaciones que describan el comportamiento de un material frente a
todas las acciones posibles en cualquier situacion, por ello es preferible anali-
zar lo que denominamos modelos ideales de comportamiento, asociados
a comportamientos simplicados que reproduciran con suciente precision
la respuesta de algunas sustancias frente a determinadas acciones con valo-
res comprendidos en un cierto intervalo. En principio consideraremos tres
tipos de comportamiento ideal: el modelo elastico, el modelo plastico y el
modelo viscoso. Para ilustrarlos es conveniente recurrir , en primer lugar,
al analisis del comportamiento uniaxial instantaneo del acero estructural en
traccion, que se observa en la gura y se puede obtener mediante un ensayo
normalizado en el laboratorio. En el diagrama se representa en ordenadas
la tension nominal (fuerza de traccion dividida por el area inicial de la
probeta) frente al alargamiento unitario nominal (deformacion lineal de la
muestra medida como la relacion entre el incremento de longitud entre dos
marcas y la distancia inicial entre ellas).
4-1
CAP

ITULO 4. ECUACIONES CONSTITUTIVAS 4-2

En el primer tramo OA del diagrama la respuesta es sensiblemente lineal
y proporcional, de forma que = E . La constante de proporcionalidad
es el conocido modulo de elasticidad o modulo de Young. El tramo AB
es corto y se caracteriza por la perdida de proporcionalidad. La totalidad
del recorrido OB se caracteriza por el hecho de que toda la deformacion
es recuperable: si reducimos el nivel de tension la deformacion se reduce y
el punto (, ) recorre la graca en sentido inverso por el mismo camino
inicial. A se denomina lmite de proporcionalidad. Estos dos tramos iniciales
denen lo que entendemos por respuesta elastica ideal: la deformacion
se recupera completamente una vez cesa la fuerza que la produce.
Una vez sobrepasado el punto B (lmite elastico) se produce un aumento
no controlado de la deformacion sin aumento apreciable de la fuerza (uen-
cia). Las deformaciones alcanzadas en este tramo ya no se pueden recuperar
totalmente. Si se lleva a cabo la descarga de la probeta se retrocede por una
curva que tiene una forma similar al tramo OB, de forma que al alcanzar el
origen de tensiones queda una deformacion remanente o plastica, que ya no
es recuperable. Al cargar de nuevo se regresa hasta el punto donde se haba
iniciado la descarga. Esta zona es caracterstica de la respuesta plastica.
Alcanzado un determinado punto C se produce un cierto endurecimiento
de la muestra, y nalmente la fractura de la misma.
Para explicar el modelo de comportamiento viscoso podemos recurrir
a un ensayo de compresion sobre una probeta de hormigon llevado a cabo
ahora con cargas de larga duracion. En el primer tramo del ensayo la tension
aumenta hasta un cierto valor
0
que produce una deformacion
0
. Si se
mantiene el valor de la tension y se procede a medir la evolucion de la
deformacion en el tiempo se observa que sigue aumentando progresivamente
hasta un cierto valor al que tiende asintoticamente. El comportamiento en
descarga es analogo y la recuperacion de parte de la deformacion se produce
lentamente. Este comportamiento es caracterstico de los materiales viscosos.
Los modelos ideales considerados son los siguientes:
Modelo elastico. Estudia la respuesta instantanea del material. Se
caracteriza por que la deformacion producida por una fuerza se recu-
pera totalmente cuando esta cesa.
Modelo plastico. Estudia la respuesta instantanea del material. Se
caracteriza por que la deformacion no se recupera cuando cesa la fuerza
que la produjo.
Modelo viscoso. Estudia la respuesta diferida del material. Se carac-
teriza por que la deformacion sigue aumentando a lo largo del tiempo
para un nivel constante de tension.
C. L

AZARO INTRODUCCI

ON A LA MEC

ANICA DE S

OLIDOS
CAP

ITULO 4. ECUACIONES CONSTITUTIVAS 4-3

4.2. Modelos elasticos: clasicacion
Dentro de los modelos elasticos, a su vez, es posible establecer una cla-
sicacion
Modelo hiperelastico Un material se denomina hiperelastico cuan-
do existe una cierta funcion W() de las componentes del tensor de
deformaciones para la que se cumple
=
W

. (4.1)
La funcion W act ua como potencial de las tensiones.
Modelo elastico Se dice que un material es elastico si existe un estado
natural en el que se encuentra libre de tensiones, y alrededor de ese
estado es posible denir una funcion uno a uno entre las componentes
del tensor de tensiones y las del tensor de deformaciones
= F(). (4.2)
Modelo hipoelastico El modelo hipoelastico es un modelo que rela-
ciona tensiones y deformaciones a nivel incremental
d = F(d, ). (4.3)
4.3. Energa de deformacion y energa
complementaria
Consideremos un solido libre de fuerzas en la conguracion de referencia.
Supongamos que sobre el act ua un sistema de fuerzas de volumen b y de
supercie

t que se aplica gradualmente hasta alcanzar su valor nal. Tras la
aplicacion de las cargas el solido se deforma hasta su conguracion actual
denida por el campo de desplazamientos d. El trabajo realizado por las
fuerzas exteriores es
T =
_
V
_
d
0
b dd dV +
_
V
_
d
0

t dd dA. (4.4)
Por el principio de conservacion de la energa, el trabajo de las fuerzas
exteriores T debe transformarse en energa durante el proceso de cambio
de forma. Si bien parte de esta energa puede perderse, por ejemplo, en la
formacion de suras o en forma de calor, otra parte de ella queda acumulada
en el material como una forma de energa potencial. La parte que se acumula
en el material se denomina energa de deformacion y se denota mediante
la letra U.
C. L

AZARO INTRODUCCI

ON A LA MEC

ANICA DE S

OLIDOS
CAP

ITULO 4. ECUACIONES CONSTITUTIVAS 4-4

Cuando la totalidad del trabajo de las fuerzas exteriores se transforma
en energa de deformacion a lo largo del proceso de cambio de forma, el
material se puede calicar de conservativo y se cumple
T = U.
El material hiperelastico ideal re une esta propiedad, y es, por tanto, capaz
de devolver integramente la energa que ha almacenado.
Mediante un smil sencillo podemos describir la energa de deformacion
como el area encerrada bajo la curva que describe el proceso carga des-
plazamiento. A partir de esta descripcion de la energa es posible introducir
una magnitud con un sentido fsico mas dicil de interpretar. Se trata de la
energa complementaria, denida como el area que queda sobre la curva
carga desplazamiento, hasta el nivel de carga alcanzado.
4.4. El modelo hiperelastico
Ya hemos introducido el material hiperelastico como el material elasti-
co conservativo, capaz de devolver la energa acumulada. Introduciremos a
continuacion el concepto de densidad de energa de deformacion.
4.4.1. Densidad de energa de deformacion
Por analoga con la denicion de un campo de fuerzas conservativo en
el que la fuerza ejercida sobre la partcula se obtiene como el gradiente de
una funcion potencial escalar en el material hiperelastico es posible denir
una funcion potencial de las variables cinematicas (deformaciones), cuya de-
rivada respecto de estas proporcione las variables estaticas (tensiones). Esta
funcion se denomina densidad de energa de deformacion y se denota
W(). Su integral en el volumen del solido es la energa de deformacion, y
sus derivadas respecto de las componentes de la deformacion proporcionan
las componentes de la tension:
=
W

(4.5)
U =
_
V
W() dV. (4.6)
Por otro lado, la densidad de energa de deformacion se puede interpretar
como la suma de las areas bajo cada una de las curvas de las componentes
de la tension respecto de las de la deformacion
W() =
_

0
dW =
_

0
W

: d =
_

0
: d.
C. L

AZARO INTRODUCCI

ON A LA MEC

ANICA DE S

OLIDOS
CAP

ITULO 4. ECUACIONES CONSTITUTIVAS 4-5

En notacion clasica,
W() =
_

x
0

x
d
x
+
_

y
0

y
d
y
+
_

z
0

z
d
z
+ 2
_

xy
0

xy
d
xy
+ 2
_

xz
0

xz
d
xz
+ 2
_

yz
0

yz
d
yz
. (4.7)
Haciendo el cambio de variable
ij
= 2
ij
se obtiene la expresion de la
densidad en funcion de las distorsiones angulares
W() =
_

x
0

x
d
x
+
_

y
0

y
d
y
+
_

z
0

z
d
z
+
_

xy
0

xy
d
xy
+
_

xz
0

xz
d
xz
+
_

yz
0

yz
d
yz
. (4.8)
4.4.2. Densidad de energa complementaria
4.5. El modelo elastico lineal
4.5.1. Material anisotropo
La relacion lineal entre el tensor de tensiones y el de deformaciones se
expresa as
= C. (4.9)
En componentes

ij
= C
ijkl

kl
(4.10)
En esta ecuacion, C es un tensor de cuarto orden denominado tensor de
elasticidades. Tiene por tanto 3
4
= 81 componentes. No es difcil demostrar
que la simetra de y implica las siguientes simetras de C:
C
ijkl
= C
jikl
C
ijkl
= C
ijlk
C
ijkl
= C
klij
.
Estas simetras reducen de 81 a 21 las componentes distintas del tensor.
Tambien se puede obtener este mismo resultado expresando en forma vec-
torial la ecuacion (4.9):
= C. (4.11)
En esta nueva ecuacion, es un vector que contiene las seis componentes
independientes del tensor de tensiones
= {
x

y

z

xy

xz

yz
}, (4.12)
C. L

AZARO INTRODUCCI

ON A LA MEC

ANICA DE S

OLIDOS
CAP

ITULO 4. ECUACIONES CONSTITUTIVAS 4-6

es el vector que contiene los alargamientos unitarios en las direcciones
coordenadas y las distorsiones angulares entre ellas
= {
x

y

z

xy

xz

yz
}, (4.13)
y C es la denominada matriz constitutiva, que es una matriz cuadrada
simetrica de 6 6 elementos, de los cuales, por tanto, solo 21 son indepen-
dientes.
4.5.2. Material isotropo: ecuaciones de Lame, ley de Hooke
generalizada
El tensor C debe ser un tensor isotropo de cuarto orden. A partir de
esta condicion y utilizando sus condiciones de simetra, es posible demos-
trar que el n umero de parametros independientes en C se reduce a 2, que
denominaremos y , y que la ecuacion (4.9) se transforma en
= e I + 2 , (4.14)
expresion que constituye las denominadas ecuaciones de Lame en la que
interviene la deformacion volumetrica introducida en el captulo de cinemati-
ca. La forma clasica de las ecuaciones de Lame es

x
= e + 2
x

xy
=
xy
(4.15a)

y
= e + 2
y

xz
=
xz
(4.15b)

z
= e + 2
z

yz
=
yz
, (4.15c)
y su expresion en la forma = C resulta
_

xy

xz

yz
_

_
=
_

_
+ 2 0 0 0
+ 2 0 0 0
+ 2 0 0 0
0 0 0 0 0
0 0 0 0 0
0 0 0 0 0
_

_
_

xy

xz

yz
_

_
. (4.16)
Las constantes y son los denominados parametros de Lame. Es usual
emplear la notacion alternativa G para el parametro .
Por otra parte, sumando las tres ecuaciones de Lame que proporcionan
las tensiones normales se obtiene el siguiente resultado
I

= (3 + 2) e. (4.17)
Sustituyendo en la expresion (4.14) se puede despejar la relacion inversa,
que resulta
=
1
2
_


3 + 2
I

.
C. L

AZARO INTRODUCCI

ON A LA MEC

ANICA DE S

OLIDOS
CAP

ITULO 4. ECUACIONES CONSTITUTIVAS 4-7

Si introducimos dos nuevos parametros, E y , denidos como
E =
(3 + 2)
+
=

2 ( +)
, (4.18)
con las expresiones inversas
=
E
(1 +)(1 2)
=
E
2 (1 +)
, (4.19)
es sencillo demostrar que la expresion anterior se transforma en
=
1 +
E


E
I

I. (4.20)
Esta es la expresion compacta de la Ley de Hooke generalizada; en
notacion clasica

x
=
1
E
_

x
(
y
+
z
)
_

xy
=
2 (1 +)
E

xy
(4.21a)

y
=
1
E
_

y
(
x
+
z
)
_

xz
=
2 (1 +)
E

xz
(4.21b)

z
=
1
E
_

z
(
x
+
y
)
_

yz
=
2 (1 +)
E

xy
. (4.21c)
Las constantes E y se denominan modulo de elasticidad y coeciente
de Poisson respectivamente. La expresion de la ley de Hooke generalizada
en la forma = C
1
es la siguiente
_

xy

xz

yz
_

_
=
_

_
1
E

E
0 0 0

E
1
E

E
0 0 0

E
1
E
0 0 0
0 0 0
2(1+)
E
0 0
0 0 0 0
2(1+)
E
0
0 0 0 0 0
2(1+)
E
_

_
_

xy

xz

yz
_

_
. (4.22)
4.5.3. Interpretaci on fsica de las constantes elasticas
del material de Hooke
Modulo de elasticidad
Si consideramos de nuevo el ensayo de traccion uniaxial, es facil notar
que en este ensayo
y
=
z
= 0, y ademas
xy
=
xz
=
yz
= 0.
Sustituyendo estos valores en la primera ecuacion de la Ley de Hooke
resulta

x
= E
x
,
relacion que permite una interpretacion sencilla del modulo de elasti-
cidad.
C. L

AZARO INTRODUCCI

ON A LA MEC

ANICA DE S

OLIDOS
CAP

ITULO 4. ECUACIONES CONSTITUTIVAS 4-8

Modulo de Poisson
Con las condiciones anteriores es posible calcular las deformaciones
y
y
z
en el ensayo de traccion uniaxial. El resultado es

y
=

E

x

z
=

E

x
.
Estas expresiones indican que en el ensayo uniaxial de traccion, las
dimensiones de la seccion transversal se reducen proporcionalmente al
valor del coeciente de Poisson.
Modulo de elasticidad transversal
A partir del ensayo de corte directo en la direccion X sobre caras
Z = const., en el que la unica componente no nula de la tension es
xz
se deduce a partir de la ley de Hooke

xz
=
xz
,
que muestra como = G esta relacionado con la rigidez frente a la
distorsion del material.
Modulo de deformacion volumetrica
La ecuacion (4.17) establece una relacion entre el primer invariante de
la tension y la deformacion volumetrica. Teniendo en cuenta que la
tension normal media es un tercio del primer invariante se cumple

m
=
3 + 2
3
e.
al factor de la deformacion volumetrica lo denominamos K, modulo
de deformacion volumetrica, por tanto
K =
3 + 2
3
=
E
3 (1 2)
. (4.23)
4.5.4. Energa de deformacion en el material lineal. Formula
de Clapeyron
Volviendo al material elastico lineal en su caso mas general, si la re-
lacion entre las componentes de la tension y de la deformacion es lineal,
entonces es posible simplicar la expresion clasica de la densidad de energa
de deformacion
W() =
1
2
_

x
+
y

y
+
z

z
+
xy

xy
+
xz

xz
+
yz

yz
_
. (4.24)
Esta expresion se denomina formula de Clapeyron. Si empleamos la nota-
cion vectorial introducida antes podemos calcular la densidad de energa de
deformacion como forma cuadratica de las componentes de la deformacion
W() =
1
2

T
=
1
2

T
C . (4.25)
C. L

AZARO INTRODUCCI

ON A LA MEC

ANICA DE S

OLIDOS
CAP

ITULO 4. ECUACIONES CONSTITUTIVAS 4-9

Finalmente, la energa de deformacion se obtiene integrando la densidad de
energa en el volumen del solido
U =
1
2
_
V

T
C dV. (4.26)
4.6. Lmite del modelo elastico: criterios de
plasticacion
En esta seccion se ha seguido la referencia [1].
4.6.1. Criterios para materiales metalicos
Criterio de Tresca
Se basa en limitar el diametro de la mayor circunferencia de Mohr al
valor
e0
g(
1
,
2
,
3
) = [(
1

2
)
2

2
e0
] [(
1

3
)
2

2
e0
] [(
2

3
)
2

2
e0
] <= 0.
(4.27)
Los estados tensionales posibles expresados como puntos (
1
,
2
,
3
) se en-
cuentran comprendidos en el interior de un prisma hexagonal cuyo eje viene
denido por la direccion {1 1 1} y cuyo lado es
_
2
3

e0
. La plasticacion se
produce cuando se alcanza la supercie del prisma.
Criterio de von Mises
En este caso los estados posibles se encuentran en el interior de un cilin-
dro con la misma orientacion y el mismo radio que el prisma de Tresca. Su
expresion es
g(
1
,
2
,
3
) = (
1

2
)
2
+ (
1

3
)
2
+ (
2

3
)
2
2
2
e0
= 0. (4.28)
4.6.2. Criterios para materiales geologicos
Criterio de MohrCoulomb
La plasitcacion se alcanza cuando la mayor circunferencia de Mohr es
tangente a la recta
= c + tan , (4.29)
donde c es la cohesion y el angulo de rozamiento interno. Esta condicion
se puede expresar de la forma
g(
1
,
2
,
3
) = [(
1
K
2
C)(
2
K
1
C)] [(
1
K
3
C)(
3
K
1
C)]
[(
2
K
3
C)(
2
K
3
C)] <= 0, (4.30)
C. L

AZARO INTRODUCCI

ON A LA MEC

ANICA DE S

OLIDOS
CAP

ITULO 4. ECUACIONES CONSTITUTIVAS 4-10

con
K =
1 + sin
1 sin
(4.31)
y
C = c
2 cos
1 sin
. (4.32)
La supercie de plasticacion es, en este caso, un cono de base hexagonal.
Criterio de DruckerPrager
Corresponde a una supercie conica de ecuacion
g(
1
,
2
,
3
) = (
1

2
)
2
+(
1

3
)
2
+(
2

3
)
2
A
2
(I

+B)
2
= 0. (4.33)
Bibliografa
[1] R. Irles Mas, Mec anica de Medios Continuos Para Ingenieros Geologos,
Publicaciones de la Universidad de Alicante, 2004.
C. L

AZARO INTRODUCCI

ON A LA MEC

ANICA DE S

OLIDOS