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
a eratur Temp
´ Presentacion
ic ul á dr
a
ras Estructu
Vi g
as
Placa s
s do i l Só PCVI y tratamiento numérico
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
VIII
PROBLEMAS RESUELTOS DE MECÁNICA DEL MEDIO CONTINUO
Abbreviations
Abreviaturas PVCI PVC MEF MEC MDF
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
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
Unidades (SI) longitud
m - metro
energía, trabajo, calor
masa
kg - kilogramo
potencia
tiempo temperatura
s - segundo K - Kelvin m s m s2
permeabilidad viscosidad dinámica flujo de masa
energía
J = Nm - Joules
conductividad térmica:
fuerza
N - Newton
frecuencia
presión, tensión
Pa ≡
velocidad aceleración
flujo de energía
N - Pascal m2
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
J = Nm - Joules J ≡ W Vatio s m2 Pa × s kg m2s J m2s W mK 1 ≡ Hz Hertz s kg m3 J m3
Prefijo
Símbolo
Potencia
kilo Mega Giga Tera
k
10 10 3
M G
10 6
T
1012
10 9
XII
PROBLEMAS RESUELTOS DE MECÁNICA DEL MEDIO CONTINUO
Nomenclatura
´ 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
b B
Tensor izquierdo de deformación de Cauchy-Green, tensor de deformación de Finger Tensor de deformación de Piola
B
Entropía creada interiormente
r r b( x , t )
b
m s2
Fuerzas másicas (por unidad de masa)
m s2
N m3
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 Diferencial de volumen m3
Ce [C ]
C in c Cv
Cp
r dA r da dV
Pa Pa Pa
Pa mol m3 m m
XIV
PROBLEMAS RESUELTOS DE MECÁNICA DEL MEDIO CONTINUO
E
E eˆ i ˆi , ˆj, kˆ
Tensor material de deformación Green-Lagrange, tensor de deformación de Green, tensor de deformación Green-St. Venant Tensor de deformación finita Euleriana o tensor de deformación de Almansi Módulo de elasticidad longitudinal o módulo de Young 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
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
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) Momento angular
nˆ
ˆ N r NO r p P
p p
Fuerza másicas por unidad de volumen
m m m m Pa
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
Primer tensor de tensiones de Piola-Kirchhoff, tensor de tensiones nominales o tensor de tensiones Pa Lagrangiano Presión media Pa Presión termodinámica Pa
NOTACIÓN
r r q( x , t ) Q Q
Flujo de calor o vector del flujo no convectivo
XV
J m2s
R S
Tensor ortogonal 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 de referencia)
r r ( x, t )
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 V ( x, t )
Pa
Pseudo vector tensión (configuración de referencia)
U&
r r u( x , t ) r r u( X , t ) 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) 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 LeviCivita
Vector posición coordenada material Vector posición coordenada espacial
m s m s m rad = ms s J =W s m m 1 K
m m
XVI
PROBLEMAS RESUELTOS DE MECÁNICA DEL MEDIO CONTINUO
εV
Deformación deformaciones)
ε
Tensor de deformación infinitesimal
η
Entropía específica
J kg K
κ
Módulo de deformación volumétrico
Pa
κ
Difusividad térmica
λ
Estiramiento
λ, µ ν
Constante 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 )
volumétrica
(para
pequeñas m
Densidad de masa en la configuración de referencia Densidad de masa en la configuración actual
m m m
m2 s m m Pa
kg m3 kg m3 kg m3 kg m3 kg m3
τ φ
Tensor de tensiones de Cauchy o tensor de tensiones Pa verdaderas Componente normal del vector tracción Pa Componente tangencial del vector tracción Pa Tensión media Pa Tensiones principales Pa Tensión normal octaédrica Pa Tensión tangencial octaédrica o tensión de corte Pa octaédrica Tensión de corte máximo Pa Tensor de tensiones de Kirchhoff Pa Ángulo de fricción interno
ψ
Energía libre de Helmholtz por unidad de masa
Ψ
Energía libre de Helmholtz por unidad de volumen J (densidad de energía) m3
Ψ (ε ) = Ψ e
Densidad de energía de deformación
ψ
Ángulo de dilatancia Tensor tasa del tensor de rotación material Tensor de vorticidad
PROBLEMAS RESUELTOS DE MECÁNICA DEL MEDIO CONTINUO
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
PROBLEMAS RESUELTOS DE MECÁNICA DEL MEDIO CONTINUO
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 Ejercicios Resueltos 1.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
)
(
)
( )
r
2
r − a
2
r 2 b cos 2 θ
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
r
donde β es el ángulo que forman los dos vectores a y b . Solución:
2
PROBLEMAS RESUELTOS DE MECÁNICA DEL MEDIO CONTINUO
2
Partiendo de la definición del módulo de un vector se cumple que: r r a+b
2
(
r r = a+b
r r
r
)⋅ (ar + br ) = ar ⋅ ar + ar ⋅ br + br ⋅ ar + br ⋅ br
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