0 ij = ε 0 0 0 ij = µ X 2X 0 0 µ Xl 2 0 3 − 0 µ X 22 X3 −µ 0 l 32 X 2 −µl 2 X 2 − − 20 3 3 µ Xµ l X 2 l2 3 − µ 2lX32 2 X l 2X l 2 3
∫
(
µ X 0 2X ε 3 2 0 − µ l ij = 0 X 2 − 0 ε 3 0 µ 0 X 2 µ X 0 ij 2= µ X3 l 2X 00 − µ 0l 2 2 X 3 X 2 l µ X 2 l 2 3 − µ2 XX 0 0 X 2 − 0 −µ l 22 X3332 3 µ X l 2 X 2 2 −µ 2 3 l 0 3 X 2 −µ 2 − µµ l 2 3 X l X 2 l2 2X 23 X −µ 2 l 2 3 X ll 2 3 2X l 2 3
ij
∫
∫
ε
∫
0
2 2
= 0
3
0
∫
(
0 −µ X 2 3 −µ 2 X l 2X l 2 3
∫
∫
0
[ ∫
−µ l X 2 3 l2
)]
∫
µ X 0 2X 2 3 l −µ 0 X 2 −µ 3 X 2 l2 3 −µ 2 X l 2X l 2 3
µ X 0 2X 2 3 −µ l 0 X 2 −µ 3 X 2 l2 3 −µ 2 X l 2X l 2 3
3
2
X l 2X
[
[
]
PROBLEMAS RESUELTOS DE
[
∫
(
ε
X 2
0
]
)]
)]
[
2 −µ l X 2 3 l2
µ X
)]
−µ
−µ
(σ
3
−ν
∫
∫
II
MECÁNICA DEL MEDIO CONTINUO
Nomenclature
III
Problemas Resueltos de ´ Mecanica del Medio Continuo EDUARDO WALTER VIEIRA CHAVES
IV
MECÁNICA DEL MEDIO CONTINUO
Presentación
fus
ió n
s Sue lo
-di
Hi
Fl uid os
jo Flu
Mec . de
Co nv ec ció n
T é rm i
co
´ Presentacion
ic ul á dr
a
ras Estructu
Vi g
as
Placa s
s do i l Só PVCI y Estrategias de Solución
Mov. Sólido Rígido
Ecuaciones Constitutivas
Ecuaciones Fundamentales de MMC Tensiones Cinemática del continuo
Tensores
VI
MECÁNICA DEL MEDIO CONTINUO
´ para el Alumno Guia 1) NO SE MEMORIZA EJERCICIO. 2) Una vez que la teoría haya sido estudiada, intentar resolver los ejercicios sin mirar la solución. Es importante que el alumno ante un nuevo problema desarrolle la habilidad de dar la solución al problema con los conocimientos adquiridos. 3) Tener en cuenta que, en general, un ejercicio es un caso particular de la teoría. Es muy importante saber reconocer cuando estamos ante una aproximación del caso general. 4) A veces, la solución de un ejercicio se puede obtener por varios caminos. Una vez resuelto el ejercicio, intentar verificar si existe otra forma de resolverlo. 5) Cuidado, puede haber erratas, seáis críticos...
Contenido
Contenido ABREVIATURAS ............................................................................................................................................... IX OPERADORES .................................................................................................................................................... X UNIDADES (SI) ................................................................................................................................................ XI NOTACIÓN .................................................................................................................................................. XIII FÓRMULAS ÚTILES ................................................................................................................................... XVII 1 TENSORES ......................................................................................................................... 1 1.1 VECTORES, NOTACIÓN INDICIAL ..................................................................................................... 1 1.2 OPERACIONES CON TENSORES DE ORDEN SUPERIOR ............................................................... 14 1.3 TRANSPUESTA ...................................................................................................................................... 20 1.3.1 Simetría y Antisimetría ........................................................................................................... 20 1.4 COFACTOR. ADJUNTA. TRAZA. TENSORES PARTICULARES. DETERMINANTE ....................... 26 1.5 DESCOMPOSICIÓN ADITIVA DE TENSORES .................................................................................. 41 1.6 LEY DE TRANSFORMACIÓN. INVARIANTES. .................................................................................. 42 1.7 AUTOVALORES, AUTOVECTORES Y TRANSFORMACIONES ORTOGONALES ........................... 49 1.8 REPRESENTACIÓN ESPECTRAL ........................................................................................................ 58 1.9 TEOREMA DE CAYLEY-HAMILTON ................................................................................................. 63 1.10 TENSORES ISÓTROPOS Y ANISÓTROPOS ...................................................................................... 82 1.11 DESCOMPOSICIÓN POLAR............................................................................................................... 91 1.12 TENSOR ESFÉRICO Y DESVIADOR ................................................................................................. 91 1.13 OTROS ................................................................................................................................................. 92 1.14 NOTACIÓN DE VOIGT ..................................................................................................................... 92 1.15 CAMPO DE TENSORES....................................................................................................................100 1.16 TEOREMA CON INTEGRALES ........................................................................................................122 2 CINEMÁTICA DEL CONTINUO ................................................................................. 153 2.1 DESCRIPCIÓN DEL MOVIMIENTO, DERIVADA MATERIAL, VELOCIDAD, ACELERACIÓN.................................................................................................................................153 2.2 TENSORES DE DEFORMACIÓN FINITA, DEFORMACIÓN HOMOGÉNEA ..............................181 2.3 DESCOMPOSICIÓN POLAR DEL GRADIENTE DE DEFORMACIÓN ..........................................225 2.4 DEFORMACIÓN INFINITESIMAL ....................................................................................................244 3 TENSIONES ................................................................................................................... 257 3.1 FUERZA, TENSORES DE TENSIONES, VECTOR TENSIÓN..........................................................257 3.2 ECUACIÓN DE EQUILIBRIO, TENSIONES Y DIRECCIONES PRINCIPALES ..............................268 3.3 OTRAS MEDIDAS DE TENSIÓN.......................................................................................................278 3.4 MÁXIMA TENSIÓN DE CORTE, CÍRCULO DE MOHR ..................................................................279 3.5 PARTICULARIDADES DEL TENSOR DE TENSIONES ....................................................................286 3.6 ESTADO TENSIONAL EN DOS DIMENSIONES .............................................................................302 3.7 TENSIONES EN COORDENADAS CILÍNDRICAS Y ESFÉRICAS ...................................................309
VIII
MECÁNICA DEL MEDIO CONTINUO: PROBLEMAS RESUELTOS
4 LEYES FUNDAMENTALES DE LA MECÁNICA DEL MEDIO CONTINUO..........313 4.1 INTRODUCCIÓN A PROBLEMAS DE FLUJO ................................................................................... 326 4.2 INTRODUCCIÓN A MOVIMIENTO DE SÓLIDO RÍGIDO ............................................................. 334 5 INTRODUCCIÓN A: ECUACIONES CONSTITUTIVAS, PVCI, Y ESTRATEGIAS DE SOLUCIÓN DEL PVCI ..................................................................................................371 6 ELASTICIDAD LINEAL..................................................................................................461 6.1 ELASTICIDAD TRIDIMENSIONAL ................................................................................................... 461 6.2 ELASTICIDAD BIDIMENSIONAL ..................................................................................................... 500 6.3 INTRODUCCIÓN A ELEMENTOS ESTRUCTURALES 1D ............................................................. 534 6.4 TORSIÓN ............................................................................................................................................. 549 6.5 ENERGÍA DE DEFORMACIÓN PARA ELEMENTOS 1D ............................................................... 565 6 BIBLIOGRAFÍA .............................................................................................................. 585
Abreviaturas
Abreviaturas PVCI PVC MEF MEC MDF MMC sii
Problema de Valor de Contorno Inicial Problema de Valor de Contorno Método de los Elementos Finitos Método de los Elementos de Contorno Método de las Diferencias Finitas Mecánica del Medio Continuo si y solo si
Latin i.e. et al. e.g. etc. Q.E.D. v., vs. viz.
id est et alii exempli gratia et cetera Quod Erat Demonstrandum versus vidilicet
es decir y otros por ejemplo y así sucesivamente lo que se quería demostrar versus a saber
- nu - ómicron - pi - theta - ro (rho) - sigma - tau - ypsilon - sigma - omega - xi - psi - dseta
Operadores
Operadores 〈•〉 =
• +• 2
• Tr (•)
(•) T (•) −1
paréntesis de MacAuley norma Euclidiana de • traza de (•) transpuesta de (•)
(•) esf
inversa de (•) inversa de la transpuesta de (•) parte simétrica de (•) parte antisimétrica de (•) parte esférica de (•) o parte hidrostática
(•) dev
parte desviadora de (•)
•
módulo de • salto de • producto escalar determinante de (•) Cofactor de • ; adjunta de (•) traza de (•) doble producto escalar operador diferencial escalar (Laplaciano) producto tensorial gradiente de • divergencia de • producto vectorial Primer, segundo y tercer invariantes del tensor •
(•) −T (•) sym (•) anti
[[•]] ⋅
det(•) ≡ • cof (•) adj(• ) Tr (•) : ∇2 ⊗ ∇ • ≡ grad(•) ∇ ⋅ • ≡ div (•) ∧ I • , II • , III • D• ≡ •& Dt r • •ˆ 1 I I sym ≡ I
Derivada material de • Vector Vector unitario (versor) Tensor identidad de segundo orden Tensor identidad de cuarto orden Parte simétrica del tensor identidad de cuarto orden
Unidades (SI)
Unidades (SI) m - metro kg - kilogramo s - segundo K - Kelvin
m s m s2 J = Nm - Joules N - Newton N Pa ≡ 2 - Pascal m W mK 1 ≡ Hz Hertz s
corriente eléctrica cantidad de sustancia intensidad luminosa
energía, trabajo, calor potencia permeabilidad viscosidad dinámica flujo de masa flujo de energía densidad de masa densidad de energía
Prefijo pico nano micro mili centi deci
Símbolo
Potencia
p
10 10 −12
η µ
m c d
10
−9
10 −6 10 −3 10 −2 10
Prefijo kilo Mega Giga Tera
A - ampere mol - mol cd - candela
J = Nm - Joules J ≡ W Vatio s m2 Pa × s kg m2s J m2s kg m3 J m3
Símbolo
Potencia
k
10 10 3
M G
10 6
T
1012
10 9
XII
MECÁNICA DEL MEDIO CONTINUO: PROBLEMAS RESUELTOS
Constantes Fisicas ´ Constante de gravitación Universal de Newton: G = 6.67384 × 10 −11 Velocidad de la luz en el vacío: c = 299 792 458
m m ≈ 300 000 000 s s
m3 kg s 2
Notación
´ Notacion r r r r A( X , t ) ≡ a ( X , t ) Aceleración (configuración de referencia)
A
Matriz de transformación de base
r r a ( x, t )
Aceleración (configuración actual)
B0 B ∂B
Medio continuo en la configuración de referencia - t = 0 Medio continuo en la configuración actual - t Contorno de B (frontera)
r r b( x , t ) b
m s2
N m3 Tensor izquierdo de deformación de Cauchy-Green, m 2 tensor de deformación de Finger m2
Fuerzas másicas (por unidad de masa)
B
Tensor de deformación de Piola
B
Entropía creada interiormente
b
m s2
J sK J Manantial de entropía local por unidad de masa y por unidad de tiempo kg s K
c
Tensor constitutivo elástico Matriz elástica (notación de Voigt) Tensor constitutivo inelástico Tensor de deformación de Cauchy Calor específico a volumen constante Calor específico a presión constante Cohesión
cc
Concentración
C
Tensor derecho de deformación de Cauchy-Green
DV
Deformación volumétrica
D
Tensor velocidad de deformación o tensor tasa de deformación o tensor tasa de deformación Euleriana o tensor estiramiento Diferencial de área en la configuración de referencia m2 Diferencial de área en la configuración actual m2
Ce [C ]
in
C c Cv
Cp
r dA r da
Pa Pa Pa
Pa mol m3 m3 m3
MECÁNICA DEL MEDIO CONTINUO: PROBLEMAS RESUELTOS
XIV
Diferencial de volumen m3 Tensor material de deformación Green-Lagrange, tensor m 2 de deformación de Green m2
dV E
ˆi , ˆj, kˆ
Tensor de deformación finita Euleriana o tensor de m 2 deformación de Almansi m2 Módulo de elasticidad longitudinal o módulo de Young Pa Base Cartesiana en notación simbólica Base Cartesiana
F
Gradiente de deformación
G H
Módulo de elasticidad transversal Tensor de deformación de Biot
H
Entropía total
r HO
Momento angular
J
Determinante del Jacobiano
e E eˆ i
r
J ( X , t) r
Tensor gradiente material de los desplazamientos
j ( x, t )
Tensor gradiente espacial de los desplazamientos
r J
Tensor de difusividad
K
Tensor de conductividad térmica
K
Energía cinética
r L
Cantidad de movimiento lineal
l
Tensor gradiente espacial de velocidad
m M
ˆ N
Masa total Tensor de tensiones de Mandel Vector unitario normal a una superficie (configuración actual) Vector unitario normal a una superficie (configuración de referencia)
r r p = ρb
Fuerza másicas por unidad de volumen
nˆ
m m Pa J K kgm 2 = Js s m3 m3 m m m m mol m2s W J = mK smK J kg m s m sm kg Pa
N m3
p p
Primer tensor de tensiones de Piola-Kirchhoff, tensor Pa de tensiones nominales Presión media Pa Presión termodinámica Pa
r r q( x , t )
Flujo de calor o vector del flujo no convectivo
Q
Tensor ortogonal
P
J m2s
NOTACIÓN
Q r r ( x, t )
XV
R S
Potencia calorífica J Función escalar que describe en forma espacial el calor J generado por las fuentes internas por unidad de masa kg s Tensor ortogonal de la descomposición polar Segundo tensor de tensiones de Piola-Kirchhoff Pa
s
r
Flujo de entropía
T r ˆ r t (n) ( x , t , nˆ ) r (Nˆ ) t0 r T ( x, t ) t t0 ≡ t = 0
Tensor de tensiones de Biot Vector tracción (configuración actual)
Pa
Temperatura Tiempo Tiempo inicial
K s s J =W s J kg
Tasa de la energía interna
u
Energía interna específica
r U( X , t r V ( x, t )
Pa
Pseudo vector tensión (configuración de referencia)
U&
r r u( x , t ) r r u( X , t )
J kg s m 2 Pa
Vector desplazamiento (Euleriana) m Vector desplazamiento (Lagrangiana) m Tensor derecho de estiramiento, o tensor de estiramiento Lagrangiano, o tensor de estiramiento material Tensor izquierdo de estiramiento, o tensor de estiramiento Euleriano, o tensor de estiramiento espacial
r r r r V ( X , t ) ≡ v ( X , t ) Velocidad (configuración de referencia)
m s m s m rad = ms s J =W s m m 1 K
r r v ( x, t )
Velocidad (configuración actual)
W
Tensor spin o tensor velocidad de rotación
w int r X r x
Potencia tensorial
α
Coeficiente de expansión térmica
δ ij ε1 , ε 2 , ε 3
Delta de Kronecker Deformaciones principales
ε
Alargamiento unitario
ijk
Símbolo de permutación, componentes del tensor Levi-Civita
εV
Deformación deformaciones)
ε
Tensor de deformación infinitesimal
Vector posición coordenada material Vector posición coordenada espacial
volumétrica
m m
(para
pequeñas m 3 m3 m m
MECÁNICA DEL MEDIO CONTINUO: PROBLEMAS RESUELTOS
XVI
η
Entropía específica
J kg K
κ
Módulo de deformación volumétrica
Pa
κ
Difusividad térmica
λ
Estiramiento
λ, µ ν
Constantes de Lamé Coeficiente de Poisson
ρ
Densidad de masa
ρS
Densidad de masa de la solución
ρf
Densidad de masa del fluido r
ρ 0 ( x, t ) r
ρ ( x, t ) σ r σN r σS σm σ1 , σ 2 , σ 3 r σ oct r τ oct τ max
τ φ ψ Ψ
Densidad de masa en la configuración de referencia Densidad de masa en la configuración actual Tensor de tensiones de Cauchy o tensor de tensiones verdaderas Componente normal del vector tracción Componente tangencial del vector tracción Tensión media Tensiones principales Tensión normal octaédrica Tensión tangencial octaédrica o tensión de corte octaédrica Tensión de corte máximo Tensor de tensiones de Kirchhoff Ángulo de fricción interno Energía libre de Helmholtz específica (por unidad de masa) Densidad de energía libre de Helmholtz (por unidad de volumen)
Ψ (ε ) = Ψ e
Densidad de energía de deformación
ψ
Ángulo de dilatancia Tensor tasa del tensor de rotación material
Ω
m2 s m m Pa
kg m3 kg m3 kg m3 kg m3 kg m3 Pa Pa Pa Pa Pa Pa Pa Pa Pa
Lista de identidades trigonométricas http://en.wikipedia.org/wiki/Trigonometric_identity
Algunas Expansión en Serie ∂f 1 ∂2 f 1 ∂3 f 2 ( x − a) + ( x − a ) + ( x − a ) 3 + L (serie de Taylor) 2 3 ∂x 2! ∂x 3! ∂x n n − ( 1 ) x 2 + L ; ( x < 1) (serie binomial) (1 + x) n = 1 + nx + 2! 1 1 exp x = 1 + x + x 2 + x 3 + L 2! 3! 1 1 Ln(1 + x) = x − x 3 + x 5 − L 3! 5! 1 2 1 4 cos( x) = 1 − x + x − L 2! 4! 1 3 1 5 sin( x) = x − x + x − L 3! 5! 1 1 cosh( x) = 1 + x 2 + x 4 + L 2! 4! 1 3 1 5 sinh( x) = x + x + x + L 3! 5! 1 1 π tan( x) = x + x 3 + x 5 + L x < 3 15 2 f ( x) = f (a ) +
Algunas Derivadas d d x d 1 d 1 (exp x ) = exp x ; (a ) = Ln(a) a x ; [Ln( x)] = ; [log a ( x)] = dx dx dx x dx xLn(a ) d 1 ∂f ( x) [Ln( f ( x))] = dx f ( x) ∂x donde e ≡ exp es el exponencial y Ln es el logaritmo natural, donde se cumple que: Ln(exp x ) = x
and
d [sin( x)] = cos( x) ; dx 1 d [arcsin( x)] = dx 1 − x2
exp Ln( x ) = x
d d [cos( x)] = − sin( x) ; [tan( x)] = sec 2 ( x) dx dx d d −1 1 ; [arccos(x)] = ; [arctan(x)] = 2 dx dx 1 + x2 1− x
Lista de derivadas http://en.wikipedia.org/wiki/List_of_derivatives
FÓRMULAS ÚTILES
XIX
Algunas Integrales
∫ exp dx = exp x
x
1
∫ x dx = Ln( x) ∫ ∫
∫
;
∂f ( x) exp f ( x ) dx = exp f ( x ) ∂x
∫ Ln( x)dx = xLn( x) − x + C
;
u = sin −1 + C a a −u du 1 u = tan −1 + C 2 2 a a +u a
∫u
du
2
2
du 2
u −a
2
=
u 1 sec −1 + C a a
Lista de integrales http://en.wikipedia.org/wiki/List_of_integrals
Solución de Funciones Función cuadrática ax 2 + bx + c = 0
solución →
x=
− b ± b 2 − 4ac 2a
( a ≠ 0)
Regla de Ruffini http://en.wikipedia.org/wiki/Ruffini%27s_rule
Expresiones relacionadas con el círculo: Ecuación del círculo: ( x1 − a ) 2 + ( x2 − b) 2 = r 2 Área del círculo: A = πr 2 Longitud de la circunferencia: C = 2πr x2 r θ b
a
x1
XX
MECÁNICA DEL MEDIO CONTINUO: PROBLEMAS RESUELTOS
Expresiones relacionadas con la elipse: x2 r x
b
f2
θ
f1
x1
b
a
r
Ecuación de la elipse: x = r = a 2 − b2 a2 Área de la elipse: A = πab
Excentricidad: e =
a
p 1 + e cos θ
;
0 < e < 1 , donde a 2 =
p2 se cumple. (1 − e 2 ) 2
1 Tensores La notación indicial fue introducida por ‘Einstein (1916, sec. 5), who later jested to a friend, "I have made a great discovery in mathematics; I have suppressed the summation sign every time that the summation must be made over an index which occurs twice..." (Kollros 1956; Pais 1982, p. 216). ‘ Ref. (Wolfram MathWorld (Einstein Summation))
1.1 Vectores, Notación Indicial Ejemplo 1.1
r
r
Probar que si a y b son vectores se cumple que:
(ar ∧ br )⋅ (ar ∧ br ) = (ar ⋅ ar )(br ⋅ br ) − (ar ⋅ br )
2
Solución:
)
(
(ar ∧ br )⋅ (ar ∧ br )
2 r r 2 r r = a ∧ b = a b sin θ r 2 r 2 r 2 r 2 r 2 r = a b sin 2 θ = a b 1 − cos 2 θ = a b 2 r 2 r 2 r r r 2 r 2 r r 2 = a b − a b cos θ = a b − a ⋅ b r r r r r r 2 = (a ⋅ a ) b ⋅ b − a ⋅ b r r r 2 r r r 2 donde hemos considerado que a ⋅ a = a y b ⋅ b = b .
(
( ) ( )
Ejemplo 1.2 r
r
r
)
(
)
( )
2
r − a
2
r b
2
cos 2 θ
r
Probar que: si c = a + b , el módulo de c puede ser expresado a través de la siguiente relación: r c =
r a
2
r r r + 2 a b cos β + b r
2
r
donde β es el ángulo que forman los dos vectores a y b . Solución: Partiendo de la definición del módulo de un vector se cumple que: r r a+b
2
Universidad de Castilla- La Mancha Ciudad Real - España
(
r r = a+b
)⋅ (ar + br ) = ar ⋅ ar + ar ⋅ br + br ⋅ ar + br ⋅ br Draft
Por: Eduardo W. V. Chaves (2014)
MECÁNICA DEL MEDIO CONTINUO: PROBLEMAS RESUELTOS
2
r r
r
r r
r
2
r r
2
Teniendo en cuenta que a ⋅ a = a , b ⋅ b = b
r r
y que a ⋅ b = b ⋅ a (conmutativo), concluimos
que: r r a+b
r r r r r r r r = a⋅a + a⋅b + b ⋅ a + b ⋅b
2
r = a r = a
2
r r r + 2a ⋅ b + b
2
r r r + 2 a b cos β + b
2
β
r
r a
demostración que a − b =
2
r r a+b =
r a
2
r r r + 2 a b cos β + b
r r r − 2 a b cos β + b r
r
r
2
2
r
r
r r a+b
2
2
2
. Luego es de fácil
. r r
r
NOTA: Partiendo de la expresión a + b = a + 2a ⋅ b + b valor a + b
r a
2
con lo cual demostramos que r
r b
2
2
podemos concluir que el
será máximo cuando β = 0º resultando que
r = a
2
r r r + 2a ⋅ b + b
r = a
2
r r r +2a b + b
(
r r = a + b
2
β = 0º 2
r b
r a
)
r r r r a+b = a + b
2
r
r
r
r
Luego para cualquier otro valor de 0º < β ≤ 180 º el valor a + b será menor que a + b . r
r
r
r
luego, a + b ≤ a + b :
r b
r r r r r c = a+b ≤ a + b
r r r c = a+b r b
r a
r
r
r
r
r
r
De forma análoga se puede demostrar que a ≤ c + b y b ≤ a + c que es la conocida desigualdad triangular, donde se cumple que: a
c