Elementos de La Mecanica Del Medio Continuo - B5
Elementos de La Mecanica Del Medio Continuo - B5
Elementos de La Mecanica Del Medio Continuo - B5
Índice .................................................................................................... 1
Capítulo 1 Conceptos generales ........................................................... 7
1.1. Concepto de Mecánica del Medio Continuo ........................ 7
1.2. Matrices................................................................................ 8
1.3. Operadores especiales ........................................................ 11
1.3.1. Delta de Kronecker ..................................................... 11
1.3.2. Operador permutación ................................................ 11
1.4. Notación indicial ................................................................ 12
1.5. Vectores ............................................................................. 16
1.5.1. Norma de un vector y producto escalar ...................... 17
1.5.2. Producto vectorial ....................................................... 18
1.5.3. Notación matricial ...................................................... 19
1.5.4. Transformación de sistema coordenado ..................... 20
1.6. Tensores ............................................................................. 22
1.6.1. Producto diádico y producto tensorial ........................ 22
1.6.2. Definición de tensor de segundo orden ...................... 25
1.6.3. Multiplicación de un tensor por un escalar y suma de
tensores 26
1.6.4. Tensores característicos de segundo orden ................. 26
1.6.5. Tensores de orden n .................................................... 28
1.6.6. Productos entre tensores ............................................. 29
1.6.7. Notación matricial aplicada a tensores ....................... 32
1.6.8. Notación de Voigt ....................................................... 35
1.6.9. Valores y direcciones principales de un tensor de
segundo orden .................................................................................. 37
1.6.10. Invariantes de un tensor de segundo orden ................. 41
1.6.11. Obtención de los valores y direcciones principales de
un tensor simétrico de segundo orden .............................................. 42
1.6.12. Parte esférica y desviadora de un tensor ..................... 46
1.7. Calculo tensorial ................................................................ 49
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 2
1.2. Matrices
pijk i j j k k i
~ 1
(1.4)
2
siendo i, j y k las variables enteras cuyo valor es 1, 2 o 3. La
permutación par e impar establece dos formas de ordenar los tres índices
como se puede esquematizar en la siguiente figura.
1 1
3 2 3 2
bi Aij c j (1.12)
ai 3
a
i (1.17)
xi i 1 xi
1.5. Vectores
x2 x2
a2 e 2
e2
e2 a ai ei
e1 x1 x1
e 3 e1 a1e1
a3e3
e3
x3 (a) x3 (b)
Figura 1.2. Vectores: (a) vectores base de un sistema coordenado, (b)
descomposición en componentes rectangulares de un vector a.
a b a b cos (1.23)
ei e j ij (1.26)
a b a b sin (1.27)
ab ~
pijk ei a j bk (1.29)
x2
a b
b bi ei
x1
a ai ei
x3
Figura 1.3. Producto vectorial.
e1 e2 e3 , e2 e3 e1 , e3 e1 e2 (1.30)
a a1 a2 a3
T
(1.31)
x2 x2 x2
x1
a2 e 2 a ai ei
a ai ei a 2 e2 a1e1
e2 x1
x1 e2
a1e1 e1
a3e3 e 3 e1 e3
x3 a3 e3
x3 (a) (b)
x3
Figura 1.4. Vector a definido en dos sistemas coordenados: (a) sistema coordenado
original, (b) sistema coordenado nuevo.
M ij ei e j (1.33)
ei M ije j (1.34)
e j M ij ei (1.35)
ai M ij a j (1.36)
1.6. Tensores
a b a b a b (1.42)
a b c a b a c . b c a b a c a (1.43)
En consecuencia (a b)T b a .
ai M ip a p , bi M iqbq (1.45)
x2
x2 x2
e2
e2 x1
e1
e1
e3 x1 x1
x3 e3
(a) x3 (b)
x3
Figura 1.5. Vectores base: (a) sistema coordenado original, (b) sistema coordenado
nuevo.
Apq M pi M qj Aij
(1.50)
M pi Aij M Tjq
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 26
Aij 1
2
A A A A
ij ji
1
2 ij ji
(1.53)
1 ij (ei e j ) (1.54)
det A 16 ~
pijk ~
prst Air Ajs Akt (1.56)
c a a 1 a
Dado que el tensor unidad es simétrico, es decir que
ij ji , se demuestra que a 1 1 a como lo indica el
problema anterior.
Para la segunda demostración se define un tensor de
segundo orden D Dij ei e j tal que:
D B 1 Dij ei e j Bik ei e k kj e k e j Bik kj ei e k e k e j
Dij ei e j Bik kj ei e j Dij Bik kj
D11 B1111 B12 21 B13 31 B11
D12 B1112 B12 22 B13 32 B12
D33 B3113 B32 23 B33 33 B33
D B B 1 B
a Ma (1.64)
B MBM T
(1.65)
Producto ci ~
c ab pijk a j bk
vectorial
Producto
tensorial
C a b Cij aib j C abT
Producto
contraído tensor - c Ab ci Aijb j c Ab
vector
Producto
contraído vector - c bA c j bi Aij c AT b
tensor
Producto
contraído vector - aCb aiCij b j a T Cb
tensor - vector
Producto
contraído tensor - C AB Cij Aik Bkj C AB
tensor
Producto
contraído tensor – C A BT Cij Aik B jk C ABT
tensor trans.
Producto doble
contraído tensor - A:B Aij Bij
tensor
Producto doble
contraído tensor 4 Cij Aijkl Bkl
C A:B
orden – tensor 2
orden
Traza de un ij Aij Aii
trA 1 : A
tensor
Determinante de
det A ~
pijk A1i A2 j A3k det A
un tensor
Gradiente de un
campo escalar
a grad ai xi a
Tensor derivada o
gradiente de un B a a Bij ai x j B aT
campo vectorial
Gradiente B sa 1 a j ai
simétrico de un Bij 1
aT 12 aT
vector 1
2 a a 2 xi x j
2
T a
Divergencia de div a a ai xi
campo vectorial
A : B Aij Bij
(1.71)
A11B11 A22B22 A33B33 2 A23B23 2 A13B13 2 A12B12
Los tensores A y B se pueden representar mediante las matrices
columna A y B de la forma:
A A11 A22 A33 A23 A13 A12
T
(1.72)
B B11 B22 B33 B23 B13 B12
T
B11
B
22
B
A11 A22 A33 A23 A13 A12 33
2 B23 (1.73)
2 B13
2 B12
A11B11 A22 B22 A33B33 2 A23B23 2 A13B13 2 A12 B12
Es habitual que las entidades asociadas a la cinemática como el campo
de la deformación incluyan un factor de 2 en los tres últimos términos de
la matriz columna como se explicará más adelante.
Por otro lado, las componentes de un tensor simétrico de cuarto orden
C se pueden expresar mediante una matriz cuadrada tamaño 6. Si de
nuevo A y B son tensores simétricos de segundo orden, la operación
tensorial A C : B , escrita en notación indicial como Aij Cijkl Bkl ,
también puede expresarse matricialmente en notación de Voigt de la
Conceptos generales 37
A x x (1.76)
A .
Se construye la matriz P (de 3 por 3), ubicando un vector propio
normalizado transpuesto x (de 1 por 3) en cada una de sus filas, así:
( r )T
Conceptos generales 39
x(1)T
P x( 2)T (1.82)
x(3)T
P , se escogen los vectores propios x(1) y x( 2) ortogonales entre sí,
manteniendo la validez de las ecuaciones (1.80) a (1.85).
Si los tres valores propios son iguales, es decir si (1) ( 2) (3) ,
cualquier vector corresponde a uno de los vectores propios de la matriz
A . Para conformar una matriz de transformación ortogonal P , los
vectores propios se escogen como cualquier grupo de tres vectores
unitarios mutuamente ortogonales.
Sea A un tensor de segundo orden que se expresa mediante sus
componentes Aij en la base e i en la matriz A Aij o mediante sus
componentes Aij en la base ei , conformando la matriz A Aij .
A partir de la ecuación característica de la matriz A expresada de la
forma det A I 0 y recordando que A M A M y
T
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 40
det M det A I det M 0
T
det M det M 1
T
se concluye que:
donde es un valor propio de la matriz A . Los valores propios de
la matriz de componentes del tensor A son cantidades intrínsecas e
independientes del sistema coordenado de referencia.
Si A es un tensor simétrico sus valores propios serán tres números
reales A(1) , A( 2) , A(3) denominados componentes principales o valores
principales de A. Si además A(1) , A( 2) , A(3) son positivos, A se denomina
tensor definido positivo. Los vectores propios del tensor A indicados
como x(1) , x( 2) , x(3) se denominan direcciones principales.
Si los valores principales son distintos las direcciones principales son
mutuamente ortogonales. En cambio, si dos o tres valores principales
son iguales, las direcciones principales dadas por los vectores
x(1) , x( 2) , x(3) no estarán definidas de forma única.
Si dos valores principales son iguales, por ejemplo si A(1) A( 2) A(3) ,
el vector propio x(3) está definido, mientras que los vectores x(1) y x( 2)
corresponden a cualquier vector unitario ortogonal a x(3) . Asimismo,
cuando los tres valores principales son iguales, es decir si
A(1) A( 2) A(3) , todo vector unitario es un vector propio del tensor A.
Con el fin de conformar un sistema coordenado con las direcciones
principales, se establece que los vectores propios x(1) , x( 2) , x(3) sean
vectores unitarios mutuamente ortogonales. Por lo tanto, el sistema
coordenado principal xˆ1 , xˆ2 , xˆ3 se define como a aquel sistema
coordenado cartesiano derecho conformado por tres ejes denominados
ejes principales, los cuales son paralelos a tres vectores unitarios
llamados vectores base principal eˆ1 , eˆ 2 , eˆ 3 .
Conceptos generales 41
trA 1 : A (1.87)
A3 I1A2 I A I 31 0 (1.94)
3 I1 2 I I 3 0 (1.96)
donde:
I1 A11 A22 A33
I A11A22 A11A33 A22 A33 A122 A132 A23
2
(1.97)
I 3 A11A22 A33 2 A12 A13 A23 A11A23
2
A22 A132 A33 A122
De acuerdo con el método de Cardano, las tres raíces reales de la
ecuación cúbica mostrada anteriormente son:
( a ) 2S cos 13 13 I1
(b ) 2S cos 13 23 13 I1 (1.98)
( c ) 2S cos 13 43 13 I1
siendo,
R 13 I12 I 2
Q 13 I1I 2 I 3 272 I13
T 1
27
R3
12
S 13 R
12
Q
arccos
2T
Los valores principales (1) , ( 2) , (3) corresponden con las raíces
obtenidas ( a ) , (b ) , ( c ) reorganizadas de mayor a menor tal forma que
(1) ( 2) (3) , es decir:
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 44
siendo,
( A22 ( p ) ) A23
1( p ) det
A23 ( A33 ( p ) )
A12 A23
(2p ) det
A13 ( A33 ( p ) )
A ( A22 ( p ) )
3( p ) det 12
A13 A23
2
( p ) 1( p ) (2p ) 3( p )
2
2
1 1 3
A 1 1 1
3 1 3
(22)
arccos 0.179853 rad 0.059951 rad
2 125 3
Las tres raíces reales de la ecuación cúbica mostrada
anteriormente se obtienen de la forma:
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 46
( a ) 2 5 cos0.059951 13 3 5.46410
(b ) 2 5 cos 0.059951 23 13 3 1.46410
( c ) 2 5 cos 0.059951 43 13 3 1.00000
Estos últimos resultados corresponden a valores principales
diferentes del tensor simétrico de segundo orden A.
Ordenados de mayor a menor se pueden indicar como:
(1) 5.46410 , ( 2) 1.00000 , (3) 1.46410
Para obtener la dirección principal asociada con el valor
(1) 5.46410 se calculan los siguientes parámetros:
1(1) 14.9282 , 1( 2) 5.46410 , 1(3) 20.3923
(1) 25.8564
x2
ei
xi
( x1 , x2 , x3 )
x1
x3
n niei ei
xi
(1.107)
n1 n2 n3
x1 x2 x3
Éste resultado indica el cambio de en la dirección n.
A A A A A A
divA 11 21 31 e 1 12 22 32 e 2
x1 x 2 x 3 x1 x 2 x 3
A A A
13 23 33 e 3
x1 x 2 x 3
ai
gra a a ei e j
x j
div a dV n adS
V V
a (1.112)
V xii dV V ni ai dS
donde n es un vector normal unitario a la superficie V en la posición
x.
a
n
dS
x2 ( x1 , x2 , x3 )
x V
e2 V
e3 e1 x1
x3
div A dV n A dS
V V
Aij (1.113)
x
V i
dV ni Aij dS
V
Capítulo 2
Descripción del esfuerzo
El esfuerzo es una medida de la intensidad de la fuerza expresada como la
fuerza distribuida por unidad de área, ocurrente dentro de un cuerpo o
sobre su contorno. En la mecánica del medio continuo un cuerpo se
considera libre de esfuerzos, si las únicas acciones presentes son las
fuerzas interatómicas necesarias para mantener unidas las partículas del
mismo. El interés de este capítulo es definir y describir el esfuerzo en
una partícula considerando su representación tensorial (Lai & Saibel
1965; Spencer 1990; Mase & Mase 1999; Holzapfel 2000; Oliver &
Agelet 2002; Reddy 2004; Kojic & Bathe 2005; Wu 2005).
t dS
n
dS
x2
dV
V
e2 V
e1 x1
e3 bdV
x3
P( x1 , x2 , x3 )
V V S
t (n )
n
dS -n
x2 S
e2 V S
t (n )
e3 e1 x1
x3
Figura 2.3. Vector tracción en un punto P actuando sobre una superficie de normal
n.
t (i ) ij e j
t (1) 11e1 12e 2 13e3
(2.3)
t ( 2 ) 21e1 22e 2 23e3
t (3) 31e1 32e 2 33e3
Recordando que el producto punto entre dos vectores es igual a la
magnitud de la proyección de un vector sobre el otro y multiplicando por
e j en ambos lados de la ecuación anterior, se tiene que:
ij t (i ) e j (2.4)
x2 x2 x2
e2 t ( 2) e2
t (1) e2
t ( 3)
e1 x1 e3 e1 x1 e1 x1
e3 P P e3 P
x3 x3 x3
x2 x2 x2
t 2(1)e 2 ( 2)
t (1) t ( 2) t e 2 2
t ( 3)
t 2(3)e 2
(1)
x1 t1( 2)e1 x1 t1(3)e1 x1
t e
3 3 (1) ( 3)
t e
1 1
t e 3 3
t3( 2)e3
x3 x3 x3
(a) (b) (c)
Figura 2.4. Componentes del vector tracción asociado a la superficie de normal (a)
e1, (b) e2, y (c) e3.
respectivamente.
x3
T
t (1) S1 n
R
e1 x2
t ( n ) S
t ( 2) S 2
e2 P
S
t (3) S3 e3
x1
Q
Figura 2.5. Fuerzas actuantes sobre un tetraedro elemental
Las fuerzas actuantes en las cuatro caras del tetraedro son las
siguientes: sobre la cara PRT actúa t (1)S1 , sobre PQT actúa t ( 2)S2 ,
sobre PQR actúa t (3)S3 y sobre la cara inclinada QTR actúa t ( n )S .
Adicionalmente se ejercen fuerzas de cuerpo en el volumen del tetraedro
V iguales a bV .
De acuerdo con la segunda ley de Newton se establece el equilibrio de
fuerzas del tetraedro elemental como la sumatoria de valores promedio de
las fuerzas actuantes igualadas a la fuerza inercial aV así,
V
t ( n ) t (1) n1 t ( 2) n2 t (3) n3 (b a) (2.7)
S
Si el tetraedro elemental se reduce de tamaño manteniendo su forma el
factor V S 0 y en consecuencia:
t ( n ) t (i ) ni (2.9)
t ( n ) ni t (i ) ni ij e j (2.10)
t (j n ) ni ij (2.11)
pq M pi ij M qj (2.15)
ij ei e j ij ei e j (2.16)
Descripción del movimiento y de la deformación 63
e1 13 (e1 2e 2 2e3 )
e2 13 (2e1 e 2 2e3 )
e3 13 (2e1 2e 2 e3 )
x3
e3
x3 e1 e2
e3 33
e1 e2 31 32
13e3
13 23
(1)
t 11
11e1 12e 2 12 21 22
x1 x1 x2
x2
(a) (b)
Figura 2.6. Componentes de esfuerzo: (a) componentes del vector tracción
actuando sobre la superficie e1, (b) representación de las componentes de esfuerzo
sobre las caras positivas de un cubo infinitesimal (se omite el vector base que
define la dirección).
t1( n ) 11 21 31 n1
(n) (2.20)
t2 12 22 32 n2
t3( n ) 13 23 33 n3
x3 x3
e3
e2
e1
33
32
31
23
22
13 21 x2
12
x1 11
x2
x1
Figura 2.7. Representación de las componentes de esfuerzo sobre las caras
positivas de un cubo infinitesimal en un sistema coordenado de vectores base ei .
t dS bdV 0
(2.21)
V V
n t* x V (2.22)
t*
n
V
dS
x2
x dV
e2 x
V
e3
x3 O e1 x1 b
n
V
i ij dS b j dV 0
V
(2.23)
bdV 0
V
(2.27)
b 0 X V (2.29)
11 21 31
b1 0
x1 x2 x3
12 22 32
b2 0 (2.30)
x1 x2 x3
13 23 33
b3 0
x1 x2 x3
La segunda condición establece que el momento resultante alrededor
de un punto O es igual a cero. Siendo x el vector posición de cada fuerza
aplicada con respecto a O, como lo indica la Figura 2.8, se tiene una
ecuación de equilibrio de la forma,
x t dS x bdV 0
*
(2.31)
V V
p p jpq x pbq dV 0
~ x n dS ~
jpq p r rq (2.32)
V V
~ x x b dV 0
p
V
jpq
xr
p rq p q
(2.34)
Descripción del movimiento y de la deformación 71
~ x x b dV 0
p
V
jpq
xr
p rq p q
(2.35)
~ x
p jpq p rq x p rq x pbq 0 (2.37)
xr xr
Factorizando x p del segundo y tercer término y recordando que
x p xr pr y que pr rq pq , se obtiene que:
~
p jpq pq x p rq bq 0 (2.38)
xr
El término entre paréntesis de la ecuación anterior corresponde a la
parte izquierda de la ecuación local de equilibrio dada por la fuerza
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 72
pq qp (2.40)
t (n ) t ( p)
n( p)
n
P P
(a) (b)
Figura 2.9. Vector tracción en P : (a) asociado a una superficie arbitraria de
normal n, (b) asociado a una superficie especial cuya normal n(p) es paralela al
mismo vector tracción
13
23 33 ( p )
Tal determinante da como resultado una ecuación cúbica para ( p )
denominada ecuación característica del tensor de esfuerzos de la forma:
x̂3 x3
ê 3
ê 2
ê1
( 3)
( 2)
x̂2
(1)
x1 x2
x̂1
Figura 2.10. Valores principales de esfuerzos y direcciones principales en un punto
P.
Descripción del movimiento y de la deformación 77
n n n
( p) 2
1
( p) 2
2
( p) 2
3 1
La dirección n(3) n1(3)e1 n2(3)e2 n3(3)e3 del esfuerzo
principal (3) 2 MPa , se obtiene de la solución del
sistema de ecuaciones simultáneas conformado por las
ecuaciones escalares linealmente independientes de la
ecuación homogénea tensorial ji ji ( p ) ni( p ) 0 para j
igual a 2 y 3, y por la norma del vector unitario de
componentes ni( p ) , así:
4n1(3) 11 6 n2(3) 2 n3(3) 0
2n1(3) 2n2(3) 14 6 n3(3) 0
n n n
( 3) 2
1
( 3) 2
2
( 3) 2
3 1
Después de resolver el sistema de ecuaciones anterior se
obtienen las siguientes componentes del vector direccional
n (3) asociado al esfuerzo principal (3) :
2 2 1
n1(3) , n2(3) , n3(3)
3 3 3
Al calcular la dirección n(1) n1(1)e1 n2(1)e2 n3(1)e3 del
esfuerzo principal (1) 5 MPa se observa que las tres
Descripción del movimiento y de la deformación 81
t ( n) N n S s (2.56)
t (n ) ss t (n )
ss t (n )
s
n s
P P n Nn P n Nn
n nˆieˆ i (2.61)
t ( n ) nˆiˆ ij eˆ j
(2.62)
nˆ1 (1)eˆ1 nˆ2 ( 2)eˆ 2 nˆ3 (3)eˆ 3
x̂2
t (n )
Nn
n x̂1
ê 2
ss s
P ê1
ê 3
x̂3
Figura 2.12. Vector tracción en P asociado a una superficie arbitraria de normal n.
Componentes normal y tangencial definidas en un sistema coordenado de vectores
base principal.
N
2nˆ 2 ( 2) (3) 0 (2.68)
nˆ 2
se obtienen las siguientes soluciones. Si (1) (3) 0 y
(2) (3) 0 , las componentes del vector n corresponden a
nˆ1 0, nˆ 2 0, nˆ3 1. En cambio, si (1) (3) 0 y
(2) (3) 0 , las componentes nˆ1 , nˆ3 no están definidas y nˆ 2 0 . Si
(1) (3) 0 y (2) (3) 0 , las componentes nˆ 2 , nˆ3 no están
definidas y nˆ1 0 . En cualquiera de estos casos la componente de
esfuerzo normal N adquiere un valor estacionario que coincide con el
esfuerzo principal de (3) .
De la misma manera como se obtuvo la componente de esfuerzo
normal de la Ecuación (2.66) en términos de las componentes nˆ1 , nˆ2 del
vector n (caso III), se calcula N en función de nˆ1 , nˆ3 (caso II) y en
función de nˆ2 , nˆ3 (caso I).
I II III IV V VI
n̂1 1 0 0 1
2
0 1
2
n̂2 0 1 0 0 1
2 1
2
n̂3 0 0 1 1
2 1
2
0
(1)n n
n x̂1 x̂1
x̂1 n
( 3) n
n ê1 I n ê 2 II n ê3 III
x̂3 x̂3 x̂3
x̂2 x̂2 1
2 23s x̂2
1
2 13s
s n
x̂1 s x̂1 x̂1
n n s
1
2 12s
s 1
2 (eˆ1 eˆ 3 ) s 1
2 (eˆ 2 eˆ 3 ) s 1
2 (eˆ1 eˆ 2 )
x̂2 IV x̂2 V
x̂2
VI
1
2 23s n
1
2 12s
s s
n x̂1
x̂1 x̂1
s n
1
2 13s
n 1
(eˆ 2 eˆ 3 ) n 1
(eˆ1 eˆ 2 )
n 1
2 (eˆ1 eˆ 3 ) x̂3 2
x̂3
2
x̂3 s 1
(eˆ1 eˆ 2 )
s 1
2 (eˆ1 eˆ 3 ) s 1
2 (eˆ 2 eˆ 3 ) 2
Figura 2.13. Planos donde actúan las componentes máximo o mínimo del
esfuerzo normal y cortante.
S2 nˆ12 (21) nˆ 22 (22) nˆ 32 (23) nˆ12 (1) nˆ 22 ( 2) nˆ 32 (3)
2
(2.69)
S2
nˆ2 2 3 2 3 2 13 nˆ12 2 2 3 nˆ22 0 (2.72)
nˆ2
donde,
13 (1) (3) , 23 ( 2) (3)
12 (1) ( 2)
(2.73)
t (n) N n S s
t1( n ) N nˆ1 S sˆ1
1
2 (1) 12 ( (1) (3) ) 1
2 12 ( (1) ( 3) ) sˆ1 sˆ1 1
2
t 2( n ) N nˆ 2 S sˆ2
(2.78)
0 12 ( (1) (3) ) 0 12 ( (1) (3) ) sˆ2 sˆ2 0
t 3( n ) N nˆ 3 S sˆ3
1
2 (3) 12 ( (1) (3) ) 1
2 12 ( (1) ( 3) ) sˆ3 sˆ3 1
2
s 1
2
eˆ 1 eˆ 3
Mediante el mismo procedimiento se puede establecer la dirección del
esfuerzo cortante máximo en los casos IV, V y VI. La Figura 2.13
muestra los planos para los cuales actúa un esfuerzo cortante máximo,
indicando el vector normal al plano n y el vector unitario s que lleva la
dirección de Smax .
N
C S2 R 2
2
donde R min R R max (2.85)
S
R max
N
min
R
C
Figura 2.14. Valores factibles de componentes de esfuerzo normal y cortante
considerando el rango de variación de n̂12 .
S
R1max
R3max
R2max
R2min
R1
min
R3min N
( 3) ( 2 ) (1)
C1
C2
C3
R3min 1
( 2)
0 nˆ32 1 C3 1
(1) (2) 2 (1)
( 2)
n̂ 3
( 3)
2 max 1
R 3 2 (1)
S
Smax 12 ( (1) (3) )
N
Nmin (3) ( 2) max
(1)
N
n 1
(eˆ1 eˆ 3 ) n 1
(eˆ 2 eˆ 3 ) n 1
2 (eˆ1 eˆ 2 )
x̂3
2
x̂3
2
x̂3
x̂2 x̂2 x̂2
n ê1 n ê 2 n ê3
x̂3 x̂3 x̂3
Figura 2.16. Esfuerzos máximos y mínimos: (a) ubicación en los círculos de Mohr
de los valores máximos de esfuerzos normales y cortantes, (b) plano sobre el cual
actúan tales esfuerzos.
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 92
M 13 tr 13 ii 13 ii (2.87)
det dev
ji ji ( p ) 0
dev
11 M (dev
p) 12 13
det 12 22 M ( p)
dev
23
0
(2.93)
13 23 33 M ( p)
dev
J
dev 3
( p) 1
dev 2
( p) J 2 (dev
p) J3 0 (2.94)
x̂2 t (n ) ~
Nn
s ~s ê 2 x̂1
~
n
ê1
ê 3 ~
n 1
eˆ1 eˆ 2 eˆ 3
3
x̂3
Figura 2.17. Vector tracción actuando sobre un plano octaédrico.
1
3
(1)
2 1
3
( 2)
2 1
3
( 3)
2
(2.96)
13 (21) (22) (23)
Asimismo, la componente de esfuerzo normal definida en la Ecuación
(2.64) para dicho plano es:
N n
~ nˆ 2 nˆ 2 nˆ 2
1 (1) 2 ( 2) 3 ( 3)
3 (1) ( 2) (3)
1
(2.97)
x2
x1
x3
x2 x2
(a) 22 x1
12
11 x1
(b)
x3
x3
(c)
Figura 2.18. Sólido en condición plana de esfuerzos en el plano x1x2: (a) viga
sometida a flexión y cortante, (b) lámina sometida a tensión, (c) estado de
esfuerzos en un elemento infinitesimal.
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 96
x2 x2
x2
e2 e2 x1
e2
e1 e1
e3 x1
e3 e3 e1 x1
x3 x3 , x3
(a) (b)
Figura 2.19. Transformación entre sistemas coordenados cartesianos para un
estado plano de esfuerzos definido por: (a) vectores base e i y (b) vectores base ei .
13 23 33 0
Recordando las relaciones trigonométricas,
cos 2 12 (1 cos 2 )
sin 2 12 (1 cos 2 )
(2.105)
cos 2 cos sin
2 2
(1) 0 0
ˆ 0 ( 2) 0 (2.110)
0 0 0
(1) , ( 2) 12 11 22 11 22 2 4122 (2.111)
S S
N N
( 3) ( 2) (1) ( 3) ( 2) (1)
S
(a) ( 3) 0 (b) ( 2) 0
N
( 3) ( 2) (1)
(c) (1) 0
Figura 2.20. Círculos de Mohr para estado plano de esfuerzos.
x̂2 x2
x2
22 ( 2)
12 (1) x̂1
11
x1
x1
x3 x3 , xˆ3 ( 3) 0
(a) (b)
2 12
tan 2 (2.120)
11 22
Si se construye un triángulo rectángulo cuyo cateto adyacente a un
ángulo 2 es 12 (11 22 ) y el cateto opuesto es 12 , se obtiene que:
R ( 12 ( 11 22 )) 2 122
1
( 11 22 )
cos 2 2
, sin 2 12
R R
Reemplazando las expresiones anteriores en la ecuación (2.106) se
obtienen los esfuerzos normales máximo y mínimo, los cuales son iguales
a los esfuerzos principales contenidos en el plano x1 x2 , presentados en las
dos primeras expresiones de la ecuación (2.109), es decir,
tan 2 s
11 22
(2.123)
2 12
Si se construye un triángulo rectángulo cuyo cateto adyacente a un
ángulo 2 s es 12 y el cateto opuesto es 12 (11 22 ) , se obtiene que:
Descripción del movimiento y de la deformación 103
R ( 12 ( 11 22 )) 2 122
12 1
( 11 22 )
cos 2 s , sin 2 s 2
R R
Reemplazando las expresiones anteriores en la tercera expresión de la
ecuación (2.106), se obtiene el esfuerzo cortante máximo y mínimo
contenido en el plano, de la forma:
12max 12 ( (b ) ( c ) ) (2.125)
S
11,12
N
( 2) (1)
22 ,12
(a)
S
Smax 12 ( (1) (3) )
N
( 3) ( 2) (1)
0 0 0
0 0 0
(2.127)
0 0 0
0 0
0 0
(2.128)
0 0 0
M 0 x2
I 0 0
0 0 0 (2.129)
0 0 0
0 Gx3 Gx2
Gx 0 0
3 (2.130)
Gx2 0 0
T (2.133)
M pi M qj (i j )
Tpqij M pi e p ei
M pi M qj M pj M qi (i j )
A partir de la Ecuación (2.104), en condición plana de esfuerzos la
relación de cambio de sistema coordenado se puede escribir como:
11 cos sin 2 2 sin cos 11
2
sin 2 cos 2
2 sin cos 22
22 (2.135)
12 sin cos sin cos cos 2 sin 2 12
11
t1( n ) n1 0 0 0 n3 n2 22
{t ( n ) } [H ( n ) ]{σ} t2( n ) 0 n2 0 n3 0 n1 33 (2.137)
t3( n ) 0
0 n3 n2 n1 0 23
13
12
En particular para la condición plana de esfuerzos la ecuación anterior
se simplifica de la forma:
t1( n ) n1 0 n2 11
(n)
n1 22
(2.138)
t 2 0 n2
12
0 100 x1 100 x2
100 x1 0 0 2
kN
m
100 x2 0 0
x2
x3
r=1
x1
es decir, t ( n) 100 x1x2e1 x12e2 x1x2e3 .
Descripción del movimiento y de la deformación 111
x2
x2 n
n 1 2 e1 3 2 e 2 n
x r=1
P x1 x2
x1
x3 x1
x3
x 1 2 e1 3 2 e2 3e3 x12
(a) (b)
x2 x2
n
r=1 n
x2 x2 e 2
x1 x3 x e x1
x3 x1 1 1
Específicamente en el punto P 1 2 , 3 2 ,3 , el vector
tracción actuando en la superficie de normal
n 1 2 e1 3 2 e 2 será igual a:
0 50 50 3 1 2 25 3
t
(n )
50 0
0 3 2 25
50 3 0 0 0 25 3
t ( n) 25 3e1 e 2 3e 3 kN m 2
La componente de esfuerzo normal es igual a N t (n) n .
En particular el esfuerzo normal a la superficie lateral del
cilindro cuya normal n 1 2 e1 3 2 e 2 , en el punto
P 1 2 , 3 2 ,3 corresponde a:
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 112
N nT t ( n )
25 3
N 1 2
3 2 0 25 25 3 kN m 2
25 3
La magnitud de la componente de esfuerzo normal a la
superficie lateral del cilindro es N t (n) n , y se puede
plantear en notación matricial de la forma:
N nT t ( n )
100 x1 x2
N x1 x2 0 100 x12 200 x12 x2
100 x1 x2
x
1
3
2e1 e2 x3e3 m
n(1)
1
e1 x1e2 x2e3 para (1) 100
2
n ( 3)
1
e1 x1e2 x2e3 para (3) 100
2
0 0 0
0 0 A Bx22 2
kN
m
0 A Bx22 Cx2 x3
h
x3
x1
h
0 0 0 0 0
t
(n)
0 0 A Bh 2 1 0
0 A Bh 2 Chx3 0 A Bh 2
es decir, t ( n) A Bh 2 e3 . Para la cara x2 h el vector
tracción, considerando que su vector normal es n e2 , es
igual a t ( n) A Bh 2 e3 . Si el vector tracción es nulo,
A Bh 2 0 A Bh 2
La componente de esfuerzo actuando sobre la cara x3 0
paralela al eje x2 es igual a ns n s , siendo n e3 y
s e2 . En consecuencia la componente ns se calcula
matricialmente de la forma:
ns sT T n
0 0 0 0
ns 0 1 0 0 0 A Bx22 0 A Bx22
0 A Bx22 0 1
h a h a B
W nsdx1dx2 A Bx dx1dx2 2a Ax2 x23
2
2
h a h a
3 h
Bh 2
W 4ah A
3
Capítulo 3
Descripción del movimiento
y la deformación
En este capítulo se presentan los conceptos generales relacionados con el
movimiento de las partículas como posición, punto material y espacial,
configuración, coordenadas materiales y espaciales, descripción
lagrangeana y euleriana, entre otras. Después se describe el
desplazamiento y la deformación en un cuerpo, indicando diferentes
medidas de deformación y profundizando específicamente en la
deformación infinitesimal (Spencer 1990; Bonet & Wood 1997;
Belytschko, Liu et al. 2000; Oliver & Agelet 2002; Reddy 2004; Wu
2005).
configuración de referencia
t 0
x2 P ( X1, X 2 , X 3 )
0
X
e2 u( X, t ) P ( x1 , x2 , x3 )
t
e3 x( X, t )
e1 x1
t t
x3 configuración actual
t 0
x2 P ( X1, X 2 , X 3 )
t t1
XP
e2 P ( x1 , x2 , x3 ) t t 2
e3 P ( x1 , x2 , x3 )
e1 x1
x3
Figura 3.2. Trayectoria de una partícula P de un cuerpo.
u XP , t x XP , t XP (3.8)
1
Las palabras distorsión y alargamiento han sido traducidas del vocablo inglés
deformation.
2
La palabra deformación es la traducción del vocablo inglés strain.
Descripción del movimiento y de la deformación 125
l L (3.12)
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 126
t 0
x
L
t t
x
l
Figura 3.4. Línea material definida entre dos partículas de un cuerpo en un espacio
unidimensional
2 3 n
La deformación logarítmica escalar permite obtener la
deformación en un espacio tridimensional, sin embargo tal
proceso de generalización es complejo y computacionalmente
costoso.
l 2 L2 1 2
G
2 L2
1
2
(3.15)
l 2 L2 1 1
A 1 2 (3.16)
2l 2
2
La Figura 3.5 compara a la deformación ingenieril con las demás
medidas de deformación graficando las deformaciones en función de
1.
0.10
E , L ,G , A Ingenieril
Green
0.08
0.06
Logarítmica
0.04
Almansi
0.02
0.00
0.00 0.02 0.04 0.06 0.08 0.10
1
1.0
E , L ,G , A
0.8 Green
Ingenieril
0.6
Logarítmica
0.4
Almansi
0.2
0.0
0.0 0.2 0.4 0.6 0.8 1.0
1
Figura 3.5. Medidas de deformación normal unidimensional
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 128
x P x XP , t (3.17)
x P l d x XP Ld, t (3.18)
t 0
x2
Q
P
Ld 0 t t
Q
e2
P ld
XP
e3 t
x1 xP
e1
x3
Figura 3.6. Movimiento en una línea material definida entre las partículas P y Q.
3
A pesar de lo indicado en las notas anteriores, en español el tensor gradiente de
deformación corresponde al vocablo inglés deformation gradient tensor.
Descripción del movimiento y de la deformación 129
f i
f i ( zr ) f i ( r ) ( z j j )
z j
z r r
(3.19)
( z j j ) 2 fi
2
( z j j )n n fi
2! z 2j n! z nj
z r r z r r
xi X rP L d r , t xi X rP , t L d j xi
X j
X ,t
r
xi
xiP l di xiP L d j (3.20)
X j
X ,t
r
xi
l di L d j X r 0
X j
X ,t
r
ei (3.22)
X i
El tensor de segundo orden F obtenido del producto tensorial entre el
vector posición x y el operador nabla material se denomina tensor
gradiente material de deformación o simplemente tensor gradiente de
deformación. Las componentes de dicho tensor Fij corresponden a los
términos xi X j , por lo tanto F se puede expresar en notación compacta
e indicial de la forma:
F x
xi (3.23)
Fij
X j
xi X j xi
ik
X j xk xk (3.26)
Fij F 1
jk ik
l d F L d , L d F 1 l d
(3.30)
ld i LFij d j , Ldi Fij1ld j
l 2 ld ld ld
2
E 1
2
C~ 1 F F 1
1
2
T
(3.35)
Eij 1
2
C~ ij ij F F 1
2 ki kj ij (3.36)
~
E 12 1 B1 12 1 FT F1 (3.42)
~
Eij 12 ij Bij1 1
2 ij Fki1Fkj1 (3.43)
(3.44)
2 1 2d E d
Se puede demostrar que el tensor E es simétrico evaluando el tensor
transpuesto de la forma:
~ T
~
T
ET 12 1 B 1 12 1 B T 12 1 F 1 F T T
1
2
1 F T
F 1 E
d 2 1 2 1
1
3
4
1
4 2 1
3 4 3
d 3 1 34 1
4 2 1 2 1
Siendo d un vector unitario, la norma del vector d es la
tasa de estiramiento:
7
1
1 2 2 2 13
d 2 2 1 2 1
2
4 3 6
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 136
J X, t u X, t
(3.45)
xX, t X F 1
y en notación indicial las componentes corresponden a:
ui
J ij X r , t
X j
(3.46)
x X i
i Fij ij
X j X j
F 1 J (3.47)
E 1
2
1 J 1 J 1 1 1 1 J J
T 1
2
T
1 JT J 1
J J
(3.48)
E 1
2
T
J J
T
1 u u u uk
Eij i j k (3.49)
2 X j X i X i X j
ui
J ij xr , t
x j
(3.51)
x X
i i ij Fij1
x j x j
1 u u u u
Eij i j k k (3.53)
2 x j xi xi x j
E J J J J
1
2
T T
0 A 0 0 0 A A2 0 0
12 0 0 A 12 A 0 0 12 0 A2 0
A 0 0 0 A 0 0 0 A2
En consecuencia,
A2 A A
E 12 A A2 A
A A A2
Para calcular el tensor de deformación de Almansi es
necesario escribir las ecuaciones de movimiento inversas.
A partir de las ecuaciones de movimiento, se despejan las
coordenadas materiales en función de las espaciales,
obteniendo:
X1
1
1 A3
x1 Ax2 A2 x3
X2
1
1 A 3
A2 x1 x2 Ax3
X3
1
1 A 3
Ax1 A2 x2 x3
De las expresiones anteriores se evalúa el campo de
desplazamientos en descripción espacial ui xi X i xr ,
como:
u1
A
1 A3
A2 x1 x2 Ax3
u2
A
1 A3
Ax1 A2 x2 x3
u3
A
1 A3
x1 Ax2 A2 x3
Por lo tanto, el gradiente de desplazamiento en descripción
espacial J ij ui x j se puede escribir en notación
matricial de la forma:
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 140
A3 A A2
J 1 3 A2 A3 A
1 A
A A2 A3
Las componentes del tensor de deformación de Almansi se
obtienen de la forma:
E 12 J J T J T J
A3 A A2 A3 A2 A
1 2 1
A A3 A A3 A2
2 1 A3
A
A2 A3
2 1 A3 A
A2 A A3
A6 A4 A2 A A A A A A
5 4 3 5 4 3
1
A A A
5 4 3
A A A A A A
6 4 2 5 4 3
2 1 A3
2
A5 A4 A3 A A A A A A
5 4 3 6 4 2
A A 0
y en consecuencia,
n
J 1 1 F 1 1 F 1
1 r r 1
r 0 (3.55)
F 1 F 1 F 1 1 F 1
2 3 n n 1
J J J 2 J 3 1 J n 1
n
(3.56)
EE 1
2
J J J J
T 1
2
T
1 u u 1 ui u j (3.60)
Eij Eij i j
2 X j X i 2 x j xi
12 J JT (3.61)
1 u u
ij i j (3.62)
2 X j X i
12 u u s u (3.63)
12 J JT (3.66)
1 u u
ij i j (3.67)
2 x j xi
12 u u su (3.68)
X 2 , x2 t t X 2 , x2
Px1 , x2 , x3
Px1 , x2 , x3 x2
r P X 1 , X 2 , X 3
e2 P X 1 , X 2 , X 3 X2
r
t 0
0
X 3 , x3 e1 X 1 , x1 X 3 , x3 x1 X 1 X 1 , x1
(a) (b)
Figura 3.7. Rotación de un cuerpo rígido: (a) esquema general, (b) coordenadas
materiales y espaciales de una partícula P.
u1 X 1 cos X 2 sin X 1
u2 X 1 sin X 2 cos X 2 (3.71)
u3 0
Después de calcular las derivadas de las componentes de
desplazamiento con respecto a las coordenadas materiales se obtiene la
matriz de componentes de indicada a continuación.
cos 1 0 0
0 cos 1 0
(3.72)
0 0 0
2 d FT F d d 1 J 1 J d
T
d 1 1 1 J J T 1 J T J d (3.73)
1 d J JT JT J d
Despreciando los términos de segundo orden del gradiente material
del desplazamiento, es decir, J T J 0 , la ecuación anterior se reduce a:
1 d J J T d (3.74)
1
n
1k 2k ! k 1 12 18 2
1n 2n! n (3.75)
k 0 1 2k k! 4 1 2nn!2 4n
2 k
1 12 d J J T d 18 d J J T d
2
(3.76)
1 12 d J J T d (3.77)
E 12 d J J T d (3.78)
1
2
J J
T
(3.79)
E d d (3.80)
E di ij d j
E d111d1 d112d 2 d113d3 d 2 21d1 d 2 22d 2 d 2 23d3 (3.81)
d 3 31d1 d3 32d 2 d 3 33d3
La operación tensorial mostrada en las ecuaciones anteriores se puede
expresar en notación matricial de la forma:
E d d
T
31 32 33 d3
3.11.2. Distorsión
Sean Q y R dos partículas ubicadas en un entorno diferencial de la
partícula P dentro de un cuerpo material, para las cuales L(Q ) d (Q ) y
L( R ) d ( R ) son los vectores definidos por las líneas materiales PQ y PR
respectivamente en la configuración de referencia. Como lo muestra la
Figura 3.8 se considera que L(Q ) d (Q ) y L( R ) d ( R ) conforman un ángulo
recto. En cambio, en la configuración actual las líneas materiales PQ y
PR definen a los nuevos vectores l (Q )d(Q ) y l ( R )d( R ) respectivamente, los
cuales conforman un ángulo .
Aplicando la definición de tensor gradiente de deformación F dada en
la Ecuación (3.30) a las líneas materiales PQ y PR se tiene que:
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 148
l (Q )d (Q ) F L(Q ) d (Q )
(3.83)
l ( R )d ( R ) F L( R ) d ( R )
L(Q ) d (Q ) F 1 l (Q ) d(Q )
(3.84)
L( R ) d ( R ) F 1 l ( R ) d( R )
t 0
X 2 , x2
R
L( R ) d ( R ) 2 Q
P (Q )
L d (Q )
0 t t
R
e2 l ( R)
d ( R)
Q
P
X
e3 P l (Q )d (Q ) t
xP
e1 X 1 , x1
X 3 , x3
Figura 3.8. Configuración actual y de referencia de dos líneas materiales en un
medio continuo.
l (Q )l ( R )d(Q ) d( R ) L(Q ) L( R ) d (Q ) 1 J 1 J d ( R )
T
(3.86)
L(Q ) L( R ) d (Q ) 1 J J T J T J d ( R )
Despreciando los términos de segundo orden del gradiente material
del desplazamiento, es decir, J T J 0 , la ecuación anterior se reduce a:
l (Q)l ( R)d(Q) d( R) L(Q) L( R) d (Q) 1 J JT d ( R) (3.87)
Descripción del movimiento y de la deformación 149
d (Q ) 1 2 d ( R )
1
cos (3.89)
( R )
(Q )
(Q ) d (Q ) 1 2 d (Q )
(3.90)
( R ) d ( R ) 1 2 d ( R )
Por lo tanto el coseno del ángulo formado por las líneas materiales PQ
y PR en la configuración actual corresponde a:
d (Q ) 1 2 d ( R )
cos
d (Q ) 1 2 d (Q ) d ( R ) 1 2 d ( R )
(3.91)
d (Q ) d ( R ) d (Q ) 2 d ( R )
1 d (Q ) 2 d (Q ) 1 d ( R ) 2 d ( R )
Como lo muestra la Figura 3.9, ahora se considera que L(Q ) d (Q ) y
L( R ) d ( R ) están contenidos en el plano X 1 X 2 y conforman un ángulo recto
entre los dos. En cambio, en la configuración actual las líneas materiales
PQ y PR definen a los nuevos vectores l (Q )d(Q ) y l ( R )d( R ) respectivamente,
también contenidos en el plano x1 x2 , los cuales conforman un ángulo
12 . Si las líneas materiales PQ y PR son paralelas a los ejes
coordenados X 1 y X 2 respectivamente, es decir el vector direccional
d (Q ) e1 y d ( R ) e2 , se cumple que:
d (Q ) d ( R ) e1 e 2 0 , d (Q ) d (Q ) 1 , d ( R ) d ( R ) 1
d (Q ) d (Q ) e1 e1 11 , d ( R ) d ( R ) e 2 e 2 22 (3.92)
d (Q ) d ( R ) e1 e 2 12
Sustituyendo los resultados anteriores en la Ecuación (3.91), se tiene
que:
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 150
212
cos 12 (3.93)
1 211 1 2 22
Dado que cos 12 sin 2 12 , el ángulo 12 expresando en la
ecuación anterior es igual a:
212
12 arcsin (3.94)
2 1 211 1 2 22
Rx, t
u1
RX,0
X 2 , x2
l ( R) ( R) d( R)
X 2
L( R )
L( R ) Qx, t
( R)
12 l (Q ) 12 d (Q )
d u2
(Q )
d (Q )
d (R) Px1 , x2 , x3 , t L(Q ) Px1 , x2 , x3 , t X 1
d (Q ) L(Q ) QX,0
(b)
e2 P X 1 , X 2 , X 3
X 3 , x3 e1 X 1 , x1
(a)
Figura 3.9. Configuración actual y de referencia de dos líneas materiales en un
medio continuo considerando la distorsión en el plano x1x2: (a) esquema general,
(b) elemento diferencial.
x̂2
R
X 2 , x2
t 0
ê 2
R P ê
1
(R) ê 3 Q
d S x̂1
d (S )
P x̂3
S
d (Q ) Q 0 t t
R
e2
d (R)
d (S ) P
d (Q )
Q t
e3 S
e1 X 1 , x1
X 3 , x3
Figura 3.10. Configuración actual y de referencia de tres líneas materiales
mutuamente ortogonales en un medio continuo paralelas a las direcciones
principales de deformación.
X 2 , x2
t 0
dX( 2)
P dX(1)
dX ( 3)
0 t t
dx ( 2)
dx (1)
e2 P
XP
e3 t
e1 xP dx (3)
X 1 , x1
X 3 , x3
Figura 3.11. Cambio de un volumen infinitesimal de un cuerpo material.
x x x
dVt ~
pijk i dX r(1) j dX s( 2) k dX t(3)
X r X s X t
(3.118)
x x j xk
dVt ~
pijk i dX r(1) dX s( 2) dX t(3)
X r X s X t
De acuerdo con la Ecuación (3.113) las componentes de los vectores
paralelos a los ejes coordenados que determinan los lados del volumen en
la configuración de referencia son:
dX dX
(1) (1)
1 dX 2(1) dX 3(1) dX
T (1)
0 0
T
dX dX
( 2) ( 2)
1 dX 2( 2) dX 3( 2) 0
T
dX ( 2) 0
T
(3.119)
dX dX
( 3) ( 3)
1 dX 2(3) dX 3(3) 0
T
0 dX (3)
T
tr 12 trJ trJT trJ (3.127)
vol tr
(3.128)
ii 11 22 33
x2
x1 x2
x2 x3 11 (b)
x1
x1 x2
x3
1
2 12 (d)
22 x1
x3
(a) x3 (c)
Figura 3.12. Sólido en condición plana de deformaciones en el plano x1x2: (a) corte
de espesor unitario de un muro de contención; estado de deformaciones de un
elemento infinitesimal: (b) deformación longitudinal en x1, (c) deformación
longitudinal en x2, (b) deformación angular en el plano x1x2.
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 160
M 11 e1 e1 cos M 12 e1 e2 sin M 13 e1 e3 0
M 21 e2 e1 sin M 22 e2 e2 cos M 23 e2 e3 0 (3.131)
M 31 e3 e1 0 M 32 e3 e 2 0 M 33 e3 e3 1
(1) , ( 2) 1
2
11 22 11 22 2 4 122 (3.138)
y en notación ingenieril a:
(1) , ( 2) 1
2
11 22 11 22 2 122 (3.139)
1
2
11, 12 12
( 2) (1)
22 , 12 12
211 2 33 2 13 2 22 2 33 2 23
, (3.147)
x32 x12 x1x3 x32 x22 x2x3
Ahora, derivando la componente 11 con respecto a x2 y a x3 se tiene
que:
211 2 u1
x2x3 x2x3 x1
(3.148)
2 22 2 23 2 13 2 12
2
x1x3 x1x2 x22 x2x3
(3.150)
2 33 2 23 2 13 2 12
2
x1x2 x1x3 x2x3 x32
La siguiente expresión general en notación indicial corresponde a un
conjunto de 81 ecuaciones escalares de las cuales solo 6 son diferentes no
triviales y coinciden con las ecuaciones de compatibilidad (3.146),
(3.147), (3.149) y (3.150).
2 ij 2 kl 2 ik 2 jl
0 (3.151)
xk xl xi x j x j xl xi xk
u 1
2
T
12 (11 11 22 22 33 33 2 23 23 213 13 212 12 )
Las componentes de deformación infinitesimal en un sistema
coordenado de vectores base ei se pueden obtener de la siguiente
operación matricial:
T (3.155)
b c
45º
P x1
a
Figura 3.14. Deformímetros eléctricos en un punto P.
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 168
c d d 0 12 12 22 1
2 23 2
T 1 1 1
2 2
12 13 2 23
1
33 0
obteniendo como resultado:
c 12 11 22 12
En consecuencia la componente de deformación angular
ingenieril en el plano x1 x2 es igual a:
12 2 c a b
11 1
2
12 0
12 12 22
0
0 0 0
(1) , ( 2) 1
2
11 22 11 22 2 122
Sustituyendo a 11 a , 22 b y 12 2 c a b de
la solución del problema anterior se tiene:
(1) , ( 2) 12 a b a b 2 4 c2 4 c a b a b 2
de donde se llega a:
(1) , ( 2)
a b
a c 2 b c 2
2 2
Capítulo 4
Leyes de conservación y
principio de los trabajos
virtuales
Las leyes de conservación permiten establecer que una cantidad física en
un medio continuo se conserva en el tiempo, independientemente del
comportamiento del material del cual está compuesto y del orden de
magnitud de sus deformaciones. A continuación se describen las leyes
de conservación de magnitudes como la masa, la cantidad de movimiento
o momentum lineal4, el momento de la cantidad de movimiento o
momentum angular5 y la energía. Asimismo, al final de este capítulo se
presenta brevemente el principio de los trabajos virtuales el cual permite
establecer el equilibrio de un cuerpo a partir del trabajo interno y externo
asociado a una variación del desplazamiento (Timoshenko & Young
1965; Lanczos 1970; Spencer 1990; Oñate 1995; Popov 1998; Mase &
Mase 1999; Oliver & Agelet 2002; Shabana 2008).
4
En Mecánica clásica, cantidad de movimiento es la traducción del vocablo inglés
linear momentum.
5
El momento de la cantidad de movimiento es la traducción del vocablo inglés
angular momentum.
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 172
dΨx, t Ψx, t
vx, t Ψx, t (4.4)
dt t
El segundo término del lado derecho de la igualdad anterior que
corresponde con la diferencia entre la derivada material y la derivada
local de una propiedad se denomina derivada convectiva y está definida
como:
dΨx, t Ψx, t
vx, t Ψx, t (4.5)
dt t
(t ) x, t dV
t V
(4.8)
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 174
(t ) d x, t dV
dt Vt
(4.10)
t X, t det FX, t dV dt x, t det Fx, t dV
(t ) d
0 0
V0 V0
det F det F dV0
d d
V0
dt dt
x, t dV
Vt
(4.20)
x, t dV 0
d
dt Vt
(4.21)
m v
dv d
f ma m (4.25)
dt dt
Sustituyendo la cantidad de movimiento de la Ecuación (4.24) en la
Ecuación (4.25), se obtiene:
dp d
dt dt V
f vdV (4.26)
por lo tanto,
d
dt V
vdV t dS bdV
V V
(4.28)
vdV b dV
d
dt V
(4.29)
V
b v dV 0
V
(4.31)
b v 0 x Vt (4.32)
b 0 x Vt (4.34)
x mv
dv d
x ma x m (4.36)
dt dt
Reemplazando la cantidad de movimiento de la Ecuación (4.24) en la
Ecuación (4.25) se obtiene:
d
dt V
x ma x vdV (4.37)
x t dS x bdV
x ma
V
V (4.38)
producidopor fuerzas de superficie producidopor fuerzas de cuerpo
d
dt V
~p jpq x p vq dV ~p jpq x ptqdS ~p jpq x pbq dV
V V
(4.40)
p jpq x p nr rqdS p jpq x x p rq dV
~ ~
(4.42)
V V r
Por lo tanto,
~ x x b d x v dV 0
p
V
jpq
xr
p rq p q
dt
p q
(4.45)
~ dv
p jpq pq x p rq bq q v p vq 0 (4.48)
xr dt
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 182
K 12 v v dV
V
(4.52)
K 12 vi vi dV
V
K 12 v v dV 12 v v dV
d d
(4.53)
dt V V
dt
(4.58)
P jivi bi vi dV ji vi vi ij bi vi dV
x j
V x j
V
x j
V
u v : dV Q
El tensor de segundo orden obtenido del producto diádico v se
puede escribir como la suma entre su parte simétrica y su parte
antisimétrica, así:
v 1
2
v v 12 v v
vi 1 vi v j 1 vi v j (4.62)
x j 2 x j xi 2 x j xi
D s v 12 v v W a v 12 v v
1 v v 1 v v (4.63)
Dij i j Wij i j
2 x j xi 2 x j xi
v : D : W : D : : D (4.64)
u : DdV Q
V
u D dV Q
(4.65)
ij ij
V
u : DdV 0
V
u D dV 0
(4.66)
ij ij
V
D su 1
2
u u (4.68)
u : (4.69)
: D dV t v dS b v dV
(4.70)
V V V
W f u (4.71)
u(x) 0 x u (4.72)
: u dV t udS b udV
s
(4.75)
V t V
: u dV t u dS b u dV a u dV
s
V t V V (4.76)
u(x) 0 x u
La expresión anterior se aplica a problemas dinámicos y se denomina
principio de d’Alembert.
Las componentes con respecto a la base e i del tensor de segundo
orden resultante del gradiente simétrico del desplazamiento virtual se
presentan en notación indicial como:
1 ui u j
su ei e j (4.77)
2 x j xi
Leyes de conservación y principio de los trabajos virtuales 189
donde los tres últimos términos están multiplicados por 2, con el fin
de asegurar la equivalencia entre la operación tensorial : su y la
operación matricial u, como se explicó en el Apartado 1.6.8.
T
siendo,
11 22 33 23 13 12 T
t t
1 t2 t3
T
(4.80)
b b1 b2 b3 T
u u1 u2 u3 T
Capítulo 5
Modelo constitutivo de un
material elástico lineal
En los capítulos anteriores se ha descrito de forma independiente el
campo de los esfuerzos y de las deformaciones de un cuerpo. En este
capítulo se relacionan tales entidades mecánicas del material a través de
los modelos constitutivos, considerando específicamente un
comportamiento lineal y elástico del material (Lubliner 1990; Lemaitre
1992; Maugin 1992; Simó & Hughes 1998; Runesson 1999; Holzapfel
2000; William 2002; Kachanov 2004; Kojic & Bathe 2005).
El primer apartado indica brevemente las características generales de
los modelos constitutivos y su clasificación actual según algunos autores
(Runesson, Ottosen et al. 1991; William 2002). Las siguientes secciones
describen en detalle el modelo constitutivo para materiales elásticos
lineales isótropos, definido por la generalización de la ley de Hooke.
Al final del capítulo se presenta la relación entre el campo de los
esfuerzos y el campo de las deformaciones para el estado plano de
esfuerzos y el estado plano de deformaciones, obtenida como una
simplificación del problema tridimensional.
(a) (b) (c)
Figura 5.1. Relación uniaxial esfuerzo-deformación en carga y descarga: (a)
material elástico lineal, (b) material elástico no lineal, (c) material inelástico.
6
Los modelos constitutivos independientes de su historia corresponden al vocablo
inglés rate-independent models.
7
Los modelos constitutivos independientes de su historia corresponden al vocablo
inglés rate-dependent models.
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 196
N P A
P
N N A
L
L E L L
E E
(a) (b) (c)
Figura 5.2. Curva esfuerzo - deformación: (a) material con ablandamiento, (b)
material con endurecimiento, (c) esquema de un ensayo de tracción uniaxial.
N
E
1 G
E 1
(a) (b)
Figura 5.3. Ley de Hooke en el campo unidimensional: (a) relación esfuerzo
normal vs. Deformación longitudinal, (b) relación esfuerzo cortante vs.
deformación angular.
x2
22
x1 22 22 22
x3
1
11
33
1 33
11 33 11
11
1 33
(a) (b) (c)
11
1
11 22 33
E
22 11 22 33
1
(5.26)
E
33 11 22 33
1
E
Despejando las componentes de esfuerzo normal de la expresión
anterior se obtiene las siguientes ecuaciones:
11
E
(1 )11 22 33
(1 )(1 2 )
22
E
11 (1 ) 22 33 (5.27)
(1 )(1 2 )
33
E
11 22 (1 ) 33
(1 )(1 2 )
Por otro lado, la deformación angular ingenieril en el plano x1 x2 es
proporcional al esfuerzo cortante en dicho plano de la forma 12 12 G .
De igual manera, se obtienen las siguientes relaciones en los planos x1 x3
y x2 x3 .
12 13 23
12 , 13 , 23 (5.28)
G G G
Despejando las componentes de esfuerzo cortante se tiene que:
12 G 12 , 13 G 13 , 23 G 23 (5.29)
Las ecuaciones escalares (5.27) y (5.29) indican la relación
constitutiva de un material elástico lineal isótropo.
Dado que el módulo de elasticidad a cortante se puede expresar en
términos de E y como lo indica la Ecuación (5.37), la relación entre el
esfuerzo y la deformación se puede escribir en notación tensorial indicial
y compacta así:
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 206
E E
ij kk ij ij
(1 )(1 2 ) (1 )
(5.30)
E E
tr 1
(1 )(1 2 ) (1 )
Asimismo, la ecuación constitutiva inversa en notación tensorial
indicial y compacta corresponde a:
1
ij kk ij ij
E E
(5.31)
1
tr 1
E E
11
1
11 22 33 1 (5.34)
E E
Sea un estado de deformación angular pura en el plano x1 x2 asociado
al estado plano de esfuerzo cortante puro de la forma:
12 ; 11 22 33 13 23 0 (5.35)
Aplicando la Ecuación (3.134), la deformación longitudinal con
respecto a un sistema coordenado x1 x2 , el cual conforma un ángulo de 45
grados entre los ejes x1 y x1 son:
11 12 12 12 (5.36)
Modelo constitutivo de un material elástico lineal 207
11 ( 2 ) 0 0 0 11
( 2 ) 0 0 0
22 22
33 ( 2 ) 0 0 0 33
23 0 0 23
13 0 13
12 sim 12
(5.41)
Asimismo, la ecuación constitutiva se puede expresar en función del
módulo de Young y la relación de Poisson obteniendo:
11 (1 ) 0 0 0 11
(1 ) 0 0 0
22 22
33 E (1 ) 0 0 0 33
23 (1 )(1 2 )
1
2 (1 2 ) 0 0 23
13 2 (1 2 )
1
0 13
12 sim 2 (1 2 )
12
1
(5.42)
Al despejar el tensor de la deformación de la Ecuación (5.24) se tiene
la relación constitutiva inversa, presentada en notación tensorial como:
ij Cijkl
1
kl
(5.43)
C1 :
Las ecuaciones escalares (5.26) y (5.28) definen a las 6 componentes
de deformación en términos de las componentes de esfuerzo.
En notación de Voigt, despejando la matriz columna de deformación
de la Ecuación (5.18) se tiene que:
C1 (5.44)
A partir de las ecuaciones escalares indicadas anteriormente, la
ecuación constitutiva inversa se puede escribir en notación de Voigt y en
términos de E y de la forma:
Modelo constitutivo de un material elástico lineal 209
11 1 0 0 0 11
1 0 0 0 22
22
33 1 1 0 0 0 33
23 E 2(1 ) 0 0 23
13 2(1 ) 0 13
12 sim 2(1 ) 12
(5.45)
11
E
11 22 , 22
E
11 22 (5.48)
1 2
1 2
La componente de esfuerzo cortante 12 está dada por la primera
ecuación escalar (5.29), de la forma:
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 210
E
12 12 (5.49)
2(1 )
Sustituyendo las componentes de esfuerzo en las ecuaciones escalares
(5.26) y (5.28) se tiene que:
11
1
11 22 , 22 1 11 22
E E
33 11 22 (5.50)
E
2(1 )
13 0 , 23 0 , 12 12
E
En notación matricial las componentes del tensor de deformación se
presentan de la forma:
11 22 (1 ) 12 0
1
(1 ) 12 11 22 0 (5.51)
E
0 0 ( 11 22 )
11
E
(1 )11 22
(1 )(1 2 )
22
E
11 (1 ) 22 (5.57)
(1 )(1 2 )
E
33 11 22
(1 )(1 2 )
La componente de esfuerzo cortante 12 está dada por la primera
ecuación escalar (5.29), de la forma:
E
12 12 (5.58)
2(1 )
Por lo tanto las componentes del tensor de esfuerzos en un estado
plano de deformaciones corresponden a:
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 212
(1 ) 11 22 (1 2 ) 12 0
E (1 2 ) 12 11 (1 ) 22 0
(1 )(1 2 )
0 0 ( 11 22 )
(5.59)
Dado que la componente de esfuerzo normal 33 indicada en la
Ecuación (5.56) depende de las otras componentes de esfuerzo normal, la
matriz columna que representa el esfuerzo contiene solo las componentes
en el plano x1 x2 como lo muestra la Ecuación (5.53). Asimismo, la
matriz columna de las componentes de deformación en notación de Voigt
está indicada en la Ecuación (5.52).
De acuerdo con las expresiones (5.57) y (5.58), la relación constitutiva
de un material elástico lineal en condición plana de deformaciones escrita
en notación de Voigt es de la forma:
11 (1 ) 0 11
E
22 (1 )(1 2 ) (1 ) 0 22
12 0 0 1
2
(1 2 ) 12
(5.60)
Después de evaluados los coeficientes de la matriz columna de
esfuerzos, la componente 33 se calcula como lo muestra la Ecuación
(5.56).
Capítulo 6
Planteamiento del problema
elástico lineal
Hasta el momento se han presentado los elementos de la mecánica del
medio continuo indispensables para resolver un problema mecánico, en
particular si el material es elástico, lineal e isótropo. En este capítulo se
busca plantear el problema general de elasticidad lineal relacionando las
condiciones de equilibrio, la ecuación constitutiva del material y la
descripción del esfuerzo y la deformación (Timoshenko & Goodier 1970;
Ortiz 1998; Mase & Mase 1999; Oliver & Agelet 2002).
x V
t (x, t )
n(x) n (x, t ) t (x, t )
n(x)
t 0
V dS u(x,0) u (x)
x2 V u (x,0) u (x)
V u (x)
(x,0) u
e2
t t
Vu
e3 e1 x1
x3 x Vu b(x, t )
u(x, t ) u (x, t )
(x, t ) b(x, t ) u
(x, t ) x Vt
ij (6.1)
b j uj
xi
Tal expresión tensorial contiene a 3 ecuaciones escalares para
j 1,2,3 . Es importante recordar que la ley de conservación de la
cantidad de movimiento tiene validez independientemente de la
cinemática o del modelo constitutivo que describe al cuerpo.
Si la aceleración del sólido es nula, la expresión anterior se transforma
en la condición de equilibrio estático definida como:
(x, t ) b(x, t ) 0 x Vt
ij (6.2)
b j 0
xi
uk u u
ij ij i j (6.9)
xk x j xi
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 218
uk 2ui 2u j
ij 2 b j uj
xi xk xi x j xi
uk u u
n j ni i j t j x V (6.12)
xk x
j xi
La condición de borde que establece el valor del desplazamiento en
una parte específica del contorno se expresa de la misma forma como lo
indica la Ecuación (6.6).
ij
E
(1 )(1 2 )
kk ij (1 2 ) ij (1 ) 0 ij
(6.17)
E
tr 1 (1 2 ) (1 ) 0 1
(1 )(1 2 )
Por otro lado, se considera que la transferencia de calor en un sólido
elástico está controlada por la ley de Fourier definida como:
qi
xi (6.18)
(x, t ) q(x, t )
donde es el coeficiente de conductividad térmica y q(x, t ) es un
campo vectorial que representa el flujo de calor por unidad de área y por
unidad de tiempo, el cual se puede expresar como:
qi
c
xi t (6.19)
q c
Planteamiento del problema elástico lineal 221
11
E
11 22
(1 )(1 )
22
E
11 22 (6.22)
(1 )(1 )
E
12 12
2(1 )
La relación entre la deformación infinitesimal y el desplazamiento
corresponde a 3 expresiones escalares independientes de la forma:
u1 u2 u1 u2
11 , 22 , 12 (6.23)
x1 x2 x2 x1
A partir de las ecuaciones anteriores, se pueden obtener las
componentes independientes del desplazamiento u1 ,u2 mediante la
solución de un sistema de dos ecuaciones diferenciales parciales de la
forma:
E 2 u1 2u2
(1 )(1 ) x12 x1x 2
E 2 u1 2u2
b1 0
2(1 ) x 22 x1x 2
(6.24)
E 2 u1 2u2
(1 )(1 ) x1x 2 x 22
E 2 u1 2u2
b2 0
2(1 ) x1x 2 x12
Tales expresiones se denominan ecuaciones de Navier en condición
plana de esfuerzos.
11
E
(1 )11 22
(1 )(1 2 )
22
E
11 (1 ) 22 (6.25)
(1 )(1 2 )
E
12 12
2(1 )
En cambio, las componentes de esfuerzo cortante 13, 23 son iguales
a cero y la componente de esfuerzo normal en la dirección x3 depende de
las otras dos componentes de esfuerzo normal de la forma
33 (11 22 ) .
Para un problema estático en condición plana de deformaciones las
ecuaciones escalares de equilibrio y la relación deformación infinitesimal
– desplazamiento corresponde a las ecuaciones (6.21) y (6.23),
respectivamente. La ecuación escalar de equilibrio omitida en la
expresión (6.21) establece que la componente de esfuerzo 33 es
constante con respecto a la dirección x3 .
A partir de las ecuaciones anteriores, se pueden obtener las
componentes independientes del desplazamiento u1 ,u2 mediante la
solución de un sistema de dos ecuaciones diferenciales parciales de la
forma:
E 2 u1 2u2
(1 )
(1 )(1 2 ) x12 x1x 2
E 2 u1 2u2
b1 0
2(1 ) x 22 x1x 2
(6.26)
E 2 u1 2u2
(1 )
(1 )(1 2 ) x1x 2 x 22
E 2 u1 2u2
b2 0
2(1 ) x1x 2 x12
t t1e1 t 2e 2
n e1
n e2
x2
x1
b
n e1 h n e 2
L
x3
Figura 6.2. Placa rectangular en condición plana de esfuerzos
t1e1 t2e 2 e1 σ x x1 L h x2 h
(6.35)
t1 11 t2 12
8 0, 9 0, 10 0 (6.39)
6 6 x1 2 7 x2 2 10
2
11 2 3 x12 6 4 x1 x2 12 5 x22
x22
(6.41)
8 x1 6 9 x2 2 12
2
12 3 2 x12 4 3 x1 x2 3 4 x22
x1 x2
2 7 x1 2 8 x2 11
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 228
3 3( 1 5) (6.43)
x2
x1
b
P h
L
x3
Figura 6.3. Viga en voladizo sometida a flexión y cortante
Las diferentes condiciones de borde de este problema se indican a
continuación.
Sobre la cara izquierda en el extremo libre donde n e1 , la fuerza
de superficie solo tiene componente en dirección x2 y es equivalente a la
acción de una fuerza puntual P, de tal forma que:
t n σ t2e 2 e1 σ x x1 0 h x2 h
t1 11 0, t2 12 (6.45)
P h t2b dx2 h 12b dx2
h h
que:
P
11 4h2 (6.47)
2bh
Otras condiciones de borde se presentan en las caras superior e
inferior, en las cuales no hay fuerzas de superficie aplicadas. Siendo
n e2 en tales caras, se tiene que:
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 230
t n σ 0 e2 σ x x2 h 0 x1 L
(6.48)
t1 21 0, t2 22 0
Cuando 22 0 x x2 h 0 x1 L y sustituyendo en la
Ecuación (6.41), se tiene que:
7 0, 10 3h 2
Sustituyendo los resultados anteriores en la Ecuación (6.49) se obtiene
que:
2 h 6 0 6 2 h
2h 6 0 6 2h (6.51)
2 0, 6 0 1 0
Cuando 12 0 x x2 h 0 x1 L se tiene que:
8 0, 11 3 4 h 2
Después de sustituir la Ecuación (6.47) en la última expresión de
(6.53), se obtiene:
3P P
11 , 4 (6.54)
4bh 4bh3
Sobre la cara derecha siendo n e1 , la fuerza reactiva de superficie es
equivalente a la acción de una fuerza puntual –P y de un momento flector
–PL de tal forma que:
Planteamiento del problema elástico lineal 231
t n σ t1e1 t2e 2 e1 σ x x1 L h x2 h
t1 11, t2 12
(6.55)
P h t2b dx2 h 12b dx2
h h
7 2 L (6.56)
que:
PL
9 4L (6.57)
4bh3
Después de sustituir la Ecuación (6.54) y en la última expresión, y de
sustituir este resultado en la Ecuación (6.46), se obtiene:
9 0 , 5 0 (6.58)
3 0 , 10 0 (6.59)
11
3P
2bh3
x1x2 , 22 0 , 12
4bh3
3P 2
h x22 (6.60)
P
b2 b2
P
Figura 6.4. Lámina sometida a una fuerza distribuida en un área que tiende a cero.
P P P
b2 b2 b2 b2 b2 b2
b4 (c)
b2
b
(b) med
med
(a) max 1.378 med max 2.575 med
med P A
max 1.027 med
t(I ) t ( II )
t tM t tM
V (I ) (II )
h( I )
V h ( II )
P d h (I )
P d h(II )
V ( I ) (x P , t M ) V ( II ) (x P , tM )
( I ) (x P , tM ) ( II ) (x P , tM )
uc (x, t ) ut (7.1)
ut
uc 0 ut
(b)
uc
(a)
( 2)
u u
u (1)
u
u u
(a) (b)
Figura 7.2. Criterios de fallo material: (a) curva esfuerzo – deformación de
tracción y compresión uniaxial, (b) dominio elástico del criterio de fractura de
Rankine para un estado plano de esfuerzos.
u (1) u , u ( 2) u (7.3)
( 2)
u
u u
u
u (1)
( 3)
u
Figura 7.3. Criterios de fallo material. Dominio elástico del criterio de fractura de
Rankine para un estado tridimensional de esfuerzos.
SY
max
12 Y (7.5)
( 2)
Y Y
Y (1)
Y
Y Y
(a) (b)
Figura 7.4. Criterios de fallo material: (a) curva esfuerzo normal – deformación
longitudinal de tracción y compresión uniaxial, (b) dominio elástico del criterio de
fluencia de Tresca para un estado plano de esfuerzos.
1
2 ( 2) (3) 12 Y , 1
2 (1) (3) 12 Y
(7.6)
1
2 (1) ( 2) 12 Y
( 2) (1) ( 2) (3)
( 2)
(1) Y von Mises
( 3) (a) Y Y
Y Y
( 3) Y Tresca (1)
(b)
Figura 7.5. Criterios de fallo material. Dominio elástico del criterio de fluencia de
Tresca para un estado tridimensional de esfuerzos: (a) vista general, (b) plano
normal al eje hidrostático.
2E
(1) (22) (23)
(7.11)
uˆ
1
2E
E
(21) (22) (23) (1) ( 2) ( 2) (3) (1) (3)
[ ] [ esf ] [ dev ]
(1) 0 0 M 0 0
0 ( 2) 0 0 M 0
0 0 (3) 0 0 M
(7.12)
(1) M 0 0
0 ( 2) M 0
0 0 (3) M
siendo M 13 ( (1) ( 2) (3) ) .
De acuerdo a lo anterior, la densidad de energía de deformación se
puede descomponer de la forma uˆ uˆ ext ( ijesf ) uˆ dis( ijdev ) . La densidad
de energía de deformación por extensión û ext corresponde a la energía
producida por el tensor esférico de esfuerzos, la cual está asociada al
cambio de volumen del sólido. Asimismo, la densidad de energía de
deformación por distorsión û dis es la energía generada por el tensor
desviador de esfuerzos, la cual está asociada con el cambio de forma del
sólido.
Sustituyendo en la Ecuación (7.11) a las componentes del tensor
esférico de esfuerzos dadas en la Ecuación (7.12), se obtiene la siguiente
expresión de la densidad de energía de deformación por extensión:
1 2
uˆ ext (1) (2) (3) 2 (7.13)
6E
Asimismo, la densidad de energía de deformación por distorsión,
indicada a continuación, se calculó reemplazando en la Ecuación (7.11) a
las componentes del tensor desviador de esfuerzos de la Ecuación (7.12).
uˆ dis
1
6E
( (1) ( 2) )2 ( ( 2) (3) )2 ( (1) (3) )2 (7.14)
vm 1
2
( (1) ( 2) )2 ( ( 2) (3) )2 ( (1) (3) )2 (7.18)
vm Y (7.19)
vm 3
2
oct (7.20)
oct 3
2
Y (7.21)
J 2 16 ( (1) ( 2) )2 ( ( 2) (3) )2 ( (1) (3) )2 (7.22)
3J 2 Y (7.23)
( 2) (1) ( 2) (3)
( 2)
(1) Y
( 3) (a)
Y Y
( 3) (1)
(b)
Figura 7.6. Criterios de fallo material. Dominio elástico del criterio de fluencia de
von Mises para un estado tridimensional de esfuerzos: (a) vista general, (b) plano
normal al eje hidrostático.
Criterios de fallo material 245
( 2)
Y
Y (1)
Y
Y
Figura 7.7. Criterios de fallo material. Dominio elástico del criterio de fluencia de
von Mises para un estado plano de esfuerzos
( 2)
( 2) Y (1)
Y
Tresca von Mises
von Mises
Y (1)
Tresca
Y
Rankine
Y
Y Rankine (1) ( 2)
(a) (b)
Figura 7.8. Comparación entre los criterios de Rankine, Tresca y von Mises: (a)
dominio elástico en condición plana de esfuerzos, (b) detalle de un estado de
esfuerzo cortante puro.
( 2)
ut ut
uc
ut (1)
uc uc
(a) (b)
Figura 7.9. Criterios de fallo material: (a) curva esfuerzo normal – deformación
longitudinal de tracción y compresión uniaxial, (b) dominio elástico del criterio de
Rankine para un estado plano de esfuerzos con tracción y compresión diferenciada.
uc (1) ut , uc ( 2) ut (7.26)
m m 1 (7.28)
a , b , m uc
m 1 m 1 ut
De manera similar al criterio de Tresca, el criterio de Mohr –
Coulomb establece que un punto material tiene comportamiento elástico
si el esfuerzo cortante máximo es menor que el esfuerzo cortante máximo
obtenido de la envolvente de fallo, es decir,
1
2
(1) ( 2) a ut 12 b (1) ( 2) (7.29)
S
1
2 ( (1) ( 2) )
ensayo de compresión c~ envolvente de falla
~
Smax c~ N tan
1
2 uc ~
1
2 ut
uc ( 2) (1) ut N
ensayo de tracción
1
2 ( (1) ( 2) )
1
2 uc 1
2 ut
Figura 7.10. Criterios de fallo material: Círculos de Mohr y envolvente de fallo.
Mohr - Coulomb
(1) ( 3) (1)
( 3)
(a) (b)
Figura 7.11. Criterios de fallo material. Dominio elástico del criterio de Mohr –
Coulomb para un estado tridimensional de esfuerzos: (a) vista general, (b) plano
normal al eje hidrostático.
( 2) m1 (1) ut ut
uc
ut (1)
( 2) m( (1) ut )
uc
Figura 7.12. Criterios de fallo material: dominio elástico del criterio de Mohr-
Coulomb para un estado plano de esfuerzos con tracción y compresión
diferenciada.
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 250
dp I1 J 2 (7.33)
dp
max
( 1
3
) ut (7.35)
dp
max
( 1
3
) uc (7.36)
dp dp
max
(7.39)
es decir,
(1) ( 2) (3)
(
1 (7.40)
1
6 (1) ( 2) ) 2 ( ( 2) (3) ) 2 ( (1) (3) ) 2 2
dp
max
( 2)
Drucker – Prager a partir ( 2)
de ensayos uniaxiales
Mohr –
(1) ( 2) (3) Coulomb
(1) ( 3) (1)
Drucker – Prager a partir (b)
( 3) (a)
de ensayos biaxiales
Figura 7.13. Criterios de fallo material. Dominio elástico del criterio de Drucker –
Prager para un estado tridimensional de esfuerzos: (a) vista general, (b) plano
normal al eje hidrostático.
Elementos de la mecánica del medio continuo para cuerpos sólidos. Volumen 1 252
(1) ( 2)
1
1
3
2
(1) (1) ( 2) (22) 2
dp
max
(7.41)
( 2)
Drucker-Prager a ( 2) Mohr-
partir de ensayos
Mohr-Coulomb ut uniaxiales y m>3 ut Coulomb
uc uc
ut (1) ut (1)
uc
Drucker-Prager a
partir de ensayos uc
(a) uniaxiales y m<3 (b)
( 2)
Mohr-
Coulomb ubt
(1)
ubc
(c)
Drucker-Prager a partir
de ensayos biaxiales
Figura 7.14. Criterios de fallo material: dominio elástico del criterio de Drucker –
Prager, para un estado plano de esfuerzos con tracción y compresión diferenciada:
(a) a partir de ensayos uniaxiales con relación entre resistencias menor a 3, (b) a
partir de ensayos uniaxiales con relación entre resistencias mayor a 3 y (c) a partir
de ensayos biaxiales.
dp
max
(2 1
3
) ubt (7.42)
dp
max
(2 1
3
) ubc (7.43)
T T
x3
u3 ( x1 , x2 )
x3
L
x2
x3
x1
u2
u1
Figura 8.2. Alabeo y giro de la sección transversal de una barra sometida a torsión
pura.
Torsión pura en barras prismáticas 257
x2
0 x3 L
T T
O x3
x3
L
x3 (8.2)
x2
B
B
u2
u1
A
A x2
O
x1 x1
O
(a) (b)
11 11 22 33 211 0
22 11 22 33 2 22 0
33 11 22 33 2 33 0
23 2 23 23 x1 (8.7)
x2
13 213 13 x2
x1
12 212 12 0
Con el resultado anterior se demuestra una de las consideraciones
planteada por el método de Saint-Venant, la cual establece la presencia
exclusiva de los esfuerzos cortantes 23 y 13 .
x2 t 0 n
n n
ds n
dx2 n
ds
O x1
dx1
V
(a) (b)
x2
z0
z x2 x1
z ( x1 , x2 )
z x1 membrana
z0 z0
x1
Fq Fq
presión uniforme q
Figura 8.6. Analogía de la membrana
x2
23
13
x1
x2 dA
O x1
Figura 8.7. Momento interno producido por los esfuerzos cortantes en un elemento
diferencial de área sobre la sección transversal de la barra.
h dv v dh hv (8.21)
x2
x1
x2
13
23
T x1
(a) O
x3
(b)
Figura 8.8. Barra de sección transversal circular sometida a torsión pura: (a)
esquema general, (b) sección transversal y distribución de esfuerzos cortantes.
T J (8.29)
x2
x1
x2
13
23
T x1
(a)
x3 O
(b)
Figura 8.9. Barra de sección transversal rectangular sometida a torsión pura: (a)
esquema general, (b) sección transversal y distribución de esfuerzos cortantes.
n 1
2 Z n n2 2 8
2 Z n (1) 2
0 n 1,3,5, , (8.35)
x12 b n
La ecuación anterior representa un sistema de ecuaciones diferenciales
no homogéneas en función de x1 , cuya solución general es de la forma:
n 1
n x1 n x1 8b2
Z n c1 sinh c2 cosh (1) 2
0 (8.36)
b b n3 3
Las constantes de integración de la función Z n se obtienen aplicando
las siguientes condiciones de borde. La simetría de la superficie alabeada
con respecto al centro de torsión O requiere que c1 0 y un valor de cero
de la función de Pradtl para x1 a / 2 , establece que:
n 1
8b 2
c2 (1) 2
na (8.37)
n3 3 cosh
b
En consecuencia, la función Z n y la función de Prandtl son iguales a:
nx1
n 1 cosh
8b
2
b 0
Z n (1) 2
1 (8.38)
n3 3 na
cosh
2b
nx1
n 1 cosh
8b 2 1 b n x2
(1) 2
1 cos (8.39)
3 n 1,3,5 n3 na b
cosh
2b
Tarea
Referencias