X. Oliver y C. Agelet - Cuestiones y Problemas de MMC
X. Oliver y C. Agelet - Cuestiones y Problemas de MMC
X. Oliver y C. Agelet - Cuestiones y Problemas de MMC
GH
0 HF iQ L FD G H 0 H G LRV & RQ WL Q X RV
X. Oliver Olivella
y C. Agelet de Saracibar Bosch
por
compilacin:
J. Martnez Parejo
Cuestiones resueltas
1 Descripc i n del
movimiento
CUESTIONES RESUELTAS
CR 1-1 Justificar si son ciertas o falsas las siguientes afirmaciones:
a) Si el campo de velocidades es estacionario, el campo de aceleraciones tambin lo es.
b) Si el campo de velocidades es uniforme, el campo de aceleraciones es siempre nulo.
c) Si el campo de velocidades es estacionario y el medio es incompresible el campo de
aceleraciones es siempre nulo. ( Tema 5 )
Resolucin:
a) Que el campo de velocidades sea estacionario significa lo siguiente:
v (x, t )
= 0 v ( x)
t
La derivada material de la velocidad da la aceleracin y por tanto:
a(x, t ) =
v (x, t )
+ v ( x, t ) v ( x, t ) = v ( x ) v ( x )
t
La expresin que queda no depende del tiempo. As se concluye que la afirmacin del
enunciado es cierta.
b) Si el campo de velocidades es uniforme implica que no depende de la coordenada
espacial y por tanto:
v (x, t ) v ( t )
v (x, t )
v ( t )
+ v ( x , t ) v ( x, t ) =
t
t
[v(t )]ij =
Por lo tanto, la afirmacin es falsa ya que
v j (t)
x i
=0
v ( t )
no tiene porqu ser nulo.
t
Cuestiones resueltas
v (x, t )
+ v ( x, t ) v ( x, t ) = v ( x ) v ( x )
t
CR 1-2 Determinar en qu condiciones las siguientes ecuaciones pueden ser las trayectorias
de un medio continuo:
x = X t + Z(e 2 1)(1 t )
2
2
y = Y t + Z(e e )(1 t )
z = e 2 Z(1 t ) + Zt
Resolucin:
F (X, t ) es el tensor Gradiente Material de la Deformacin, definido de la siguiente manera:
x
Fij = i
X j
Se debe cumplir:
F >0
t 0 (e 2 1)(1 t )
t0
2
2
2
t + e (1 t) = t(e 1) + e > 0
t0
2
e
t <
1.157
e 2 1
GENERALIZACIN:
Para un caso ms general, donde se tenga la Forma Cannica de las Ecuaciones del Movimiento:
x = C1 t + C 3 (e 2 1)
x = x(C1 , C 2 , C 3 , t ) y = C 2 t + C 3 (e 2 e 2 )
z = e 2 C 3 (1 t)
se deber cumplir la misma condicin F > 0 . Pero ahora las expresiones son diferentes a
las anteriores.
Tomando como instante inicial un tiempo t * :
Cuestiones resueltas
X = C1t + C 3 (e 2 1)
2
Y = C 2 t + C3 (e e )
Z = e C 3 (1 t )
y despejando C i , i {1,2,3} :
1
Z
(1 e 2 )]
C1 = [X
t
1 t
1
Z
(1 e 4 )]
C 2 = [Y
t
1
t
Z
C 3 = 2
e (1 t )
t 0
Estas ecuaciones estn definidas para
. Sustituyendo en las ecuaciones de las
t 1
trayectorias del medio continuo resulta:
t
t Z
Z
2
2
x = t X t 1 t (1 e ) + 1 t * (1 e )
t
t Z
Z
(1 e 4 ) +
(1 e 4 )
y = Y
1
t
*
t
t
1
t
1 t
Z
z =
1 t
t0
x i
y su determinante, queda la condicin:
X j
F=
O lo que es lo mismo:
t2 1 t
>0
2
t 1 t
t0
t 0
t <1 , t <1
t >1 , t >1
Resolucin:
Como se da la expresin espacial de la velocidad y se pide el valor de la aceleracin en un
punto x = ( 1,1,1 ) T del espacio, no hace falta encontrar las Ecuaciones de Movimiento.
nicamente hay que aplicar:
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
donde:
Cuestiones resueltas
dv( x, t ) v (x, t )
=
+ v ( x, t ) v ( x , t )
dt
t
v = [0 , z(e t e t ) , 0]T
t
x
0
0
1
t
t
v = [x z , z(e + e ) , 0] = 0
0
0
y
t
t
1
(
e
e
)
0
z
v v = [x z , 0 , 0]T
a = [x z , z(e t e t ) , 0]
a(x = x , t = 2) = [0 , e 2 e 2 , 0]
y=
1
(Y + Z)e t + (Y Z)e t
2
z=
1
(Y + Z)e t (Y Z)e t
2
t
y z = (Y Z)e t
Y Z = ( y z )e
X=x
1
1
Y = ( y + z )e t + ( y z )e t
2
2
1
1
t
Z = ( y + z ) e ( y z )e t
2
2
y ahora la velocidad:
x( X, t ) 1
t
t
1
v ( X, t ) =
= (Y + Z)e (Y Z)e
2
t
2
1
t
t
1
(Y + Z)e + (Y Z)e
2
Cuestiones resueltas
0
v (x, t ) = z
y
(1)
v( X, t ) 1
t
t
1
= (Y + Z)e + (Y Z)e
a( X, t ) =
t
2
2
1
t
t
1
(Y + Z)e (Y Z)e
2
(2)
(3)
Resolucin:
a) Las trayectorias cumplen las siguientes ecuaciones:
dx = at x = 1 at 2 + C
1
dt
2
dy
= y y = C2e t
dt
dz = 0 z = C
3
dt
Cuestiones resueltas
x = 1 at 2 + X
2
y = Ye t
z=Z
v f = e t y
df = e t y + e t y = ( + )e t y = 0 x t
t
dt
f(x, y, z, t) e t y C = 0 e t y = C ( + )e t y = C( + )
xt
( + )e t y = 0
x t t C( + ) = 0 t
+ = 0
t := x | F(x, y, z, t) = e 2t (x 2 + y 2 ) + z 2 e t C = 0
donde C 0 es una constante, determinar (t ) sabiendo que las partculas que estn
sobre dicha superficie son siempre las mismas.
Resolucin:
F define la superficie material t := {x | F(x, y, z, t) = 0} , y por tanto deber cumplir:
donde
dF = F + v (x, t ) F = 0 x t
t
dt t
Cuestiones resueltas
F = 2e 2t (x 2 + y 2 ) 2tz 2 e t 2
t
2xe 2t
F = 2ye 2t
2
2ze t
v F = (x , y , z(t )) F = 2x 2 e 2t + 2y 2 e 2t + 2z 2 e t (t )
2
resultando finalmente:
2z 2 ((t) t)e t = 0 x t t
F(x, y, z, t) e 2t (x 2 + y 2 ) + z 2 e t C = 0
z 2 = [C e 2t (x 2 + y 2 )]e + t
(t) = t
CR 1-7
a) Obtener la divergencia en el punto (0 , 0 , 0 )T del campo vectorial definido en
coordenadas esfricas por v = ar e r .
b) Justificar que, en el caso ms general, para cualquier campo vectorial como el de la
Fig.1 ( convergente hacia el punto A ) o el de la Fig.2 (divergente desde el punto A ), el
valor de la divergencia en A es negativo y positivo, respectivamente.
Fig.1
Fig.2
Resolucin:
a) El campo de velocidades en coordenadas esfricas es el siguiente:
v r = ar
v = 0
v = 0
La definicin de divergencia en coordenadas esfricas es:
v=
1 2
1
1 v
+
+
r
(
v
)
(
v
sen
)
r
r sen
r sen
r 2 r
1
3ar 2
3
(ar
)
=
= 3a
r 2 r
r2
Cuestiones resueltas
Para el caso de la Figura 1 resulta a < 0 ya que el vector velocidad tiene diferente
sentido que las r crecientes, y por tanto, la divergencia es negativa.
Para el caso de la Figura 2 resulta a > 0 ya que el vector velocidad tiene el mismo
sentido que las r crecientes, y por tanto, la divergencia es positiva.
b) Para demostrarlo, se toma un volumen diferencial que envuelva el punto A y se razona
sobre l:
FIGURA 1:
Los vectores n y v tienen diferente sentido, ya que n es el vector unitario
exterior del volumen dV y v va hacia dentro. Por tanto su producto ser
negativo.
vdV = ( v) dV
v
v
n
<
0
dV
dS
dV
dV
( v) A < 0
FIGURA 2:
En este caso, ocurre todo lo contrario, v y n tienen el mismo sentido y, por
tanto, su producto ser positivo.
vdV = ( v) dV
dV
dV
( v ) A > 0
Problemas Resueltos
PROBLEMAS RESUELTOS
PR 1-1 Para el siguiente campo de velocidades:
v x = 0 , v y = 0 , v z = f (x, y )z
Se pide:
a) Hallar las trayectorias y las lneas de corriente indicando cmo son.
b) Dada la funcin g(x, y, z, t ) = Ln(xyz)
f
t , determinar el valor de f(x, y) para que
x
x
rea S
Resolucin:
a) Para hallar las trayectorias se debe integrar la descripcin espacial de la velocidad ya que:
dx
= v(x, t )
dt
dz
= f(C1 , C 2 ) z
dt
Ln (z )= f(C1 , C 2 ) t + k
dz
= f(C1 , C 2 )dt
z
z(t) = C 3 e f(C1 ,C 2 ) t
10
Problemas Resueltos
x(0) = C1 = X
y(0) = C 2 = Y
z(0) = C 3 = Z
x=X
y=Y
(1)
z = Ze
f ( X , Y ) t
Para hallar las lneas de corriente, se tendra que integrar el campo de velocidades respecto
, pero en el problema el campo de velocidades en descripcin espacial no depende
explcitamente del tiempo,
v
= 0 , es decir, el campo de velocidades es estacionario y, por
t
tanto, las lneas de corriente coinciden con las trayectorias. La ecuacin de las lneas de
corriente ser entonces:
x = C1
y = C2
z = C 3 X e
f C1 , C 2
z
Las trayectorias (y las lneas de corriente) son verticales ya que
slo vara la componente z , debido a que el campo de
velocidades es vertical. Como las lneas de corriente son sus
envolventes tambin deben ser verticales.
y
x
b) Recordar que (x , t ) es superficie material si, y slo si, su derivada material es nula:
d(x , t ) (x , t )
+ v (x , t ) = 0
=
dt
t
f
t=0
x
v g = v x
vy
g
x
g
g
1
xy = f(x, y)
vz = vz
= f(x, y) z
y
z
xyz
g
z
11
Problemas Resueltos
f(x, y)
dg
=
+ f(x, y) = 0
dt
x
Si se toma f(x, y) 0 :
f(x, y)
+ f(x, y) = 0
x
f(x, y)
f(x, y)
1
= f(x, y)
=1
x
f(x, y)
x
Para el caso f (x, y ) = 0 se tiene que v(x, t ) = 0 , es decir, el medio continuo est en reposo y,
por tanto, toda superficie es material.
0
x
, siendo F =
.
X
F ( X, t)
x i
.
X j
X
y
F=
X
z
x
Y
y
Y
z
Y
Z 1 0
0
y
= 0 1
0
Z
f(X, Y)t
z * * e
Z
f (X, Y )
e f (X, Y ) t
12
Problemas Resueltos
y*= Y
z* = Z e f (X,Y) 1
X 2 + Y 2 + Z 2 e 2 f(X, Y) R 2 = 0
e) La cantidad de masa por unidad de tiempo que atraviesa una superficie se define de la
siguiente forma:
S = v ndS
(2)
S
Para este problema y segn el apartado c), la densidad tiene la siguiente expresin:
=
A
e At
Se observa que la densidad no depende de los puntos del espacio, sino que nicamente
lo hace del tiempo, por lo tanto puede salir fuera de la integral.
Y si la superficie encierra un volumen V y aplicamos el Teorema de la Divergencia queda:
S = v ndS = vdV
S
(3)
13
Problemas Resueltos
S1
S1
S2
S2
S3
S3
S3
S3
v x v y v z
+
+
)dV = (0 + 0 + A)dV = AV = ALS
V
x
y
z
_____________________________________________________________________
PR 1-2 La descripcin espacial del campo de velocidades de un fluido perfecto es:
y
v(x, t ) = z e t ,
, v z
1+ t
1
e t
, siendo 0 = ctte .
2
1+ t
Se sabe tambin que la superficie (x, t ) = x z (1 + t ) e t + k = 0 es una superficie material.
14
Problemas Resueltos
Resolucin:
a) En primer lugar, hay que utilizar la informacin sobre la superficie para poder
completar el campo de velocidades espacial, ya que el enunciado no da explcitamente
la componente v z . Que dicha superficie sea material significa que todas las partculas
que la componen son siempre las mismas, y expresado matemticamente queda:
d
=
+ v = 0
dt t
= z(e t + e t (1 + t))
t
t
t
t
t
z(e + e (1 + t)) + ze (1 + t)e v z = 0 , t
T
= 1 , 0 , (1 + t)e t
v z = z
Por lo tanto:
y
v (x, t ) = ze t ,
, z
1+ t
dy
y
=
dt 1 + t
dy
dt
=
y 1+ t
ln y = ln(1 + t ) + k y = Y(1 + t )
dy
y
dy d
=
lny =
+ k y = C 2 e 1+ t
d 1 + t
y 1+ t
1+ t
(1)
15
Problemas Resueltos
x = C 3 e t + C 1
y = C2
1
e +t
z = C 3 e
(2)
v (x, t )
+ v ( x , t ) v ( x, t )
t
dv ( X, t )
dt
T
dv ( X, t )
= [0 , 0 , Ze t ]
dt
0 ze t
p (x, t ) p (x, t ) p (x, t )
p =
=
,
,
0
,
0
,
y
z
1+ t
x
( X,0)
F ( X, t )
16
Problemas Resueltos
( X,0) = 0 :
0
t
1
F ( X, t ) = 0 1 + t 0
0
e t
0
=
F = (1 + t)e t
0et
1+ t
(3)
Y=
y = Y(1 + )
1+
z = Ze z = Ze
Z = z e
x = X + Z
(4)
Para cada instante t , se puede visualizar la lnea de traza segn el parmetro , que da
la situacin espacial de las partculas tintadas.
Se puede comprobar que para t = se debe cumplir:
x = x
y = y
z = z
ya que lo que indica es que la lnea de traza est pasando por el punto de vertido.
Ahora hay que delimitar esta lnea para cada tiempo t :
En este caso, el primer punto que se tinta es el que pasa por el punto de vertido para
= t 1 , mientras que el ltimo estar pasando, en el instante = t , por el punto de
vertido.
17
Problemas Resueltos
Un extremo es = t x = x *
x = x z e t1 (t 1 t)
y = y
El otro extremo es = t 1
1+ t
1 + t1
z = z e t1 t
t1 < t < t 2
(x
, y * , z * punto de vertido
= t1
z
=t
lnea de traza
y
x
= t1
x = x z e t1 (t 1 t)
1+ t
y = y
+ t1
1
t1 t
z=z e
t > t2
(x
= t2
x = x z e t 2 (t 2 t)
1+ t
y = y
+ t2
1
t2 t
z=z e
, y * , z * punto de vertido
= t1
z
= t2
lnea de traza
( v n)dS
(v )dV
18
Problemas Resueltos
+ (v ) = 0 (v ) =
t
t
Hay que darse cuenta que dicha derivada es constante en el espacio. As, nicamente
hay que sustituir:
=
dV = V
t
t
0 te t 4
=
R 3
2
(1 + t) 3
19
Cuestiones Propuestas
CUESTIONES PROPUESTAS
CP 1-1 Definir descripcin material y espacial de un movimiento.
CP 1-2 Enunciar e interpretar fsicamente las condiciones que debe verificar la forma
cannica de la ecuacin de movimiento: x = ( X, t )
y = Y + Zt
z = Yt 2 + Z
20
Cuestiones Propuestas
a) Un observador situado en el punto del espacio (1,1,1) que midiera la densidad de las
partculas que van pasando por dicho punto.
b) Un observador que viajara con una partcula que en el instante t = 1 ocupaba la
posicin (1,1,1) , midiendo su densidad.
d (t)
= X(t) ,
dt
v Y = Ye t
, vZ = 0
xi
(1 + t)
v = [y 2 e 2t , y, 0]
se vierte un colorante desde el instante t = 1 hasta el instante t = 2 en el punto (1,1,1) .
Obtener la ecuacin de la lnea de traza en el instante t = 5 y las coordenadas de sus
T
y = C 2 e t
z = C 3 e 2 t
determinar la ecuacin de la lnea de traza, indicando los puntos inicial y final de la misma
para t = 5 .
CP 1-14 La descripcin espacial del campo de velocidades de un fluido es:
v = [ye t , ze t , 0]
y
z
v = x ,
,
1 + t 1 + t 1 + t
Razonar si la superficie t := { x | 3xt + y + z = 0} es material.
21
Cuestiones Propuestas
y
, z , demostrar que cualquiera
CP 1-16 Dado el campo de velocidades: v = x ,
2t
2t
2t
2
2
2
2
de las superficies t := { x | x + y + z = R t}, donde R es una constante, es una
superficie material.
y , z
CP 1-17 Dado el campo de velocidades: v = x ,
, y sabiendo que la superficie
2t
2t
2t
t := { x | x 2 + y 2 + z 2 = t + } es material, encontrar todos los posible valores de y .
y = C 2 e
z = C3
23
Problemas Propuestos
PROBLEMAS PROPUESTOS
PP 1-1 Para el movimiento definido por el campo de velocidades:
v x = ye t
vy = y
vz = 0
Se pide:
S
y
S
Se pide:
v y = by ;
vz =
z
t+c
24
Problemas Propuestos
A
A
S
y
h
y
x
x
Se pide:
a + x 2 dx = x a + x 2 + a Ln x + a + x 2
2
2
ij = ij E ll + 2 E ij
, cttes
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Resueltas
'HVFULSF L yQGHOD
GHIRUPDFLyQ
&8(67,21(65(68(/7$6
&5 6H WLHQH HO VLJXLHQWH WHQVRU HVSDFLDO GH GHIRUPDFLyQ SDUD HO PRYLPLHQWR GH XQ
PHGLRFRQWLQXR
0
e(x, t ) = 0
te tz
0
te tz
0
0
0 t(2e tz e t )
t =0
z Z
$
z
E
ds
dS
y Y
y
[
x X
(QODGHIRUPDFLyQGHXQVHJPHQWRVHGHILQHHOHVWLUDPLHQWRHQODGLUHFFLyQ t FRPR
=
1
= ds
1 2t e t dS
(QHVWHFDVRHOYHFWRUXQLWDULR t YDOH
2SHUDQGRUHVXOWDORVLJXLHQWH
T
t = 1 (1 1 1)
3
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Resueltas
0
te tz 1
0
1
t e t = 1 [1 1 1] 0
= 1 te t
0
0
1
3
3
te tz 0 t(2e tz e t ) 1 3
1
1
=
t =2 =
t
2
1 + te
1 + 4 e2
3
3
l AB = dS =
A
1 ds = 1 bds = 1 3
l AB = 3 + 4e 2
&5(QXQFLHUWRLQVWDQWHHOPRYLPLHQWRGHXQPHGLRFRQWLQXRYLHQHGHILQLGRSRU
x 1 = X 1 + AX 3
x 2 = X 2 AX 3
x 3 = X 3 AX1 + AX 2
+DOODUHOWHQVRUPDWHULDOGHGHIRUPDFLyQ'HPRVWUDUTXHVL A HVSHTXHxRHOPRYLPLHQWR
SXHGHHQWHQGHUVHFRPRXQDURWDFLyQGHVyOLGRUtJLGR\FDOFXODUHOYHFWRUGHURWDFLyQ
5HVROXFLyQ
(O7HQVRU0DWHULDOGH'HIRUPDFLyQ E VHUHODFLRQDFRQHO7HQVRU*UDGLHQWHGHOD'HIRUPDFLyQ F
VHJ~QODH[SUHVLyQ
E = 1 (F T F )
2
SRUORWDQWRVHGHEHHQFRQWUDUHOYDORUGH F
0 A
1
x i
1 A
F=
= 0
j A A 1
0 A
1
T
1
A
F = 0
A A 1
1 + A 2 A 2
0
T
2
2
1+ A
0
F F= A
0
0
1 + 2A 2
(OWHQVRU E UHVXOWDILQDOPHQWH
1 1 0
2
1
E = A 1 1 0
2
0 2
0
6HJ~QODWHRUtDJHQHUDOVHGHILQHHO7HQVRU*UDGLHQWH0DWHULDOGH'HVSOD]DPLHQWRVFRPR
0 0 A
J = F = 0 0 - A
A A 0
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Resueltas
(Q HO FDVR GH A SHTXHxR VH FXPSOHQ ODV FRQGLFLRQHV GH GHIRUPDFLRQHV LQILQLWHVLPDOHV
SRUORWDQWRHO7HQVRU,QILQLWHVLPDOGH5RWDFLyQYDOGUi
0 A
0
T
1
0 A
= (J J ) = 0
2
A A 0
\HO9HFWRU,QILQLWHVLPDOGH5RWDFLyQ
23 A
= 31 = A
12 0
(OPRYLPLHQWRGHODVSDUWtFXODVSXHGHGHILQLUVHFRPR
dx = F dX
GRQGH
F = + +
\DGHPiVHQSHTXHxDVGHIRUPDFLRQHVVHFXPSOHTXH
0 0 0
T
1
E = = ( J + J ) = 0 0 0
2
0 0 0
(VWRLQGLFDTXHQRH[LVWHQGHIRUPDFLRQHVSRUORWDQWRKDEUiXQJLURLQILQLWHVLPDOGHVyOLGR
UtJLGR
dx = ( + )dX
&5(QXQFLHUWRLQVWDQWHHOFDPSRGHGHVSOD]DPLHQWRVGHXQPHGLRFRQWLQXRHV
u X = (a 1 1)X
u Y = (a 2 1)Y + a 1. X
u Z = (a 3 1)Z
y = Y + u Y = a 2 Y + a 1 . X
z = Z + u Z = a 3 Z
2EWHQLGDVpVWDVVHFDOFXODHO7HQVRU*UDGLHQWHGHOD'HIRUPDFLyQ\VXGHWHUPLQDQWH
a1
F = a 1.
0
a2
0
0
0
a 3
F = a 1a 2 a 3
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Resueltas
4XHHOVyOLGRVHDLQFRPSUHVLEOHLPSOLFDTXHQRVHPRGLILFDVXYROXPHQ\SRUORWDQWR
dV = dV0 F
F = a 1a 2 a 3 = 1
4XHXQVHJPHQWRSDUDOHORDOHMH Z QRVHDODUJXHFRQOOHYDTXH
ds = dS
= 1 + 2Tz E Tz = 1 + 2E zz = 1
GRQGHHO7HQVRU0DWHULDOGH'HIRUPDFLyQYDOH
a 1 (1 + . 2 ) 1 a 1a 2.
0
T
2
1
1
E = (F F ) =
a 1a 2.
a2 1
0
2
2
0
0
a 32 1
DVtVHFRQFOX\HORVLJXLHQWH
a3 =1
1 + 2E zz = 1 E zz = 0
a 3 = 1
3RU~OWLPRVHHVWXGLDODFRQGLFLyQGHTXHHOiUHDGHXQHOHPHQWRVLWXDGRHQHOSODQR XZ QR
VHKDPRGLILFDGR
6L VH GHILQH XQ iUHD HQ HO SODQR XZ VHJ~Q ORV YHFWRUHV dX 1 \ dX 2 \ OHV DSOLFDPRV OD
GHIRUPDFLyQFRUUHVSRQGLHQWH
1
dX1 = 0 dS1
0
(OiUHDUHVXOWDQWHVHUi
0
dX 2 = 0 dS 2
1
a1
dx1 = F dX1 = a 1. dS1
0
0
dx 2 = F dX 2 = 0 dS 2
a 3
i
dx1 dx 2 = a 1
a 1.
0
k a 1 a 3.
0 = a 1a 3 dS1dS 2
a 3 0
$OLJXDODUODViUHDVLQLFLDO\ILQDOUHVXOWD
rea inicial = dS1dS 2
rea final = (a 1 a 3.
+ (a 1a 3 ) 2 dS1 dS 2
(a 1a 3) 2 + (a1a 3 ) 2 = 1 a12a 32 2 + a 12 a 32 = 1
6LVHWLHQHHQFXHQWDTXH a 32 = 1 TXHGDORVLJXLHQWH
a12 (1 + 2 ) = 1
2
a1 = 1 +
1
1 + 2
)LQDOPHQWHVHREWLHQHQODVSRVLEOHVVROXFLRQHV
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Solucin 1 : a 1 =
Solucin 2 : a 1 =
Cuestiones Resueltas
, a 2 = 1 + .2
1
1+ .
1
Solucin 3 : a 1 =
Solucin 4 : a 1 =
1+ .
1
, a2 = 1+ .2
1+ .
1
1+ .
, a3 =1
, a 3 = 1
, a2 = 1+ .2
, a3 = 1
, a 2 = 1 + .2
, a 3 = 1
&53DUDODGHIRUPDFLyQGHILQLGDSRU
x = X + mZ ,
y=Y
z = Z mX
GRQGH m HV XQD FRQVWDQWH KDOODU HO YROXPHQ GH OD GHIRUPDGD GH XQD HVIHUD GH FHQWUR
A(0,0,0) \UDGLR R
5HVROXFLyQ
(QSULPHUOXJDUVHGHEHFDOFXODUHO7HQVRU*UDGLHQWHGHOD'HIRUPDFLyQ\VXGHWHUPLQDQWH
1 0 m
x i
F=
= 0 1 0
X
j m 0 1
F =1+ m2
7HQLHQGRHQFXHQWDTXH F QRGHSHQGHGHOHVSDFLRODUHODFLyQH[LVWHQWHHQWUHYROXPHQ
LQLFLDO\ILQDOTXHGDGHODVLJXLHQWHPDQHUD
dV = F dV0
dV = F dV
V = F V0
V = (1 + m 2 ) 4 5
3
&5 (Q XQ PHGLR FRQWLQXR VH SURGXFH XQD GHIRUPDFLyQ FRQ ODV VLJXLHQWHV
FRQVHFXHQFLDVVREUHHOWULiQJXORGHODILJXUD
\
D (OVHJPHQWR OA SDVDDPHGLU (1 + p)
YHFHVVXORQJLWXGLQLFLDO
E (OiQJXOR AOB GLVPLQX\HHQ q UDGLDQHV
VXYDORULQLFLDO
F (OiUHDSDVDDVHU (1 + r ) YHFHVODLQLFLDO
G p, q, r, s << 1
[
2
6DELHQGRTXHODGHIRUPDFLyQHVXQLIRUPHTXHHOHMH z HVXQDGLUHFFLyQSULQFLSDOGHOWHQVRU
JUDGLHQWHGHODGHIRUPDFLyQTXHpVWHHVVLPpWULFR\TXHHOHVWLUDPLHQWRHQGLFKDGLUHFFLyQ
HV z = 1 + s REWHQHUHOWHQVRUGHSHTXHxDVGHIRUPDFLRQHV
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Resueltas
5HVROXFLyQ
&RPR OD GHIRUPDFLyQ HV XQLIRUPH HQWRQFHV HO 7HQVRU *UDGLHQWH GH OD 'HIRUPDFLyQ F QR
GHSHQGHUiGHODVYDULDEOHVHVSDFLDOHV&RQVHFXHQWHPHQWHHO7HQVRUGH'HIRUPDFLyQ E \ORV
HVWLUDPLHQWRV WDPSRFR
5HFRUGDUTXHVHHVWiEDMRFRQGLFLRQHVGHGHIRUPDFLRQHVLQILQLWHVLPDOHV
6HWHQGUiTXHUHODFLRQDUODORQJLWXGLQLFLDO\ILQDOGHXQVHJPHQWRSDUDOHORDOHMH x
A
= x dX = x OA inicial
O
O
x =1+ p
= (1 + p)OA inicial
OA final =
OA final
x dX
x = 1 + xx = 1 + p xx = p
ngulo inicial =
ngulo final = + - xy
- xy = xy = 2 xy = q
xy =
q
2
F HVVLPpWULFR\HOHMH z HVGLUHFFLyQSULQFLSDO
F11
F = F12
F12
F22
0
u x
u x
u x
1 + x
u x u y
u x (x, y)
y
z
0
=
=0
u y
u y
u y
z
z
u y (x, y)
0 =+ J =
1+
y
z
u z = u z = 0 u (z)
F33
u
z
u
u
x
z
z
y
1+ z
y
z
x
u
u
xz = 1 x + z = 0 xz = 0
2 z
x
u
u
xz = 1 x + z = 0 xz = 0
2 z
x
u z
= z 1
z
zz = s
z =1+ s
zz =
$GHPiVVHJ~QODH[SUHVLyQ
F = + +
Fzz = 1 + s
/DUHODFLyQH[LVWHQWHHQWUHODViUHDVHV
dA = F dA 0 F 1
3DUDFDOFXODUODLQYHUVDGH F VHXWLOL]DODVLJXLHQWHQRWDFLyQ
B11
F = B12
B12
0
B 22
0
0 1 + s
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Resueltas
GRQGHVHFDOFXODODLQYHUVDGH B \DVtVHREWLHQHODLQYHUVDGH F
B 1
C
= 11
C12
6LGHILQLPRVHOYHFWRUiUHD
C
11
C12
1
C12
F
=
C 22
C12
C 22
0
0
0
1
1+ s
0
0
dA 0 = 0 dA 0 F 1 =
0
dA 0
dA 0
1 + s
\WHQLHQGRHQFXHQWDTXH
F = tr ( ) + 1
'HVSUHFLDQGRORVWpUPLQRVGHVHJXQGRRUGHQUHVXOWD
dA = (1 + r)dA 0
yy = r p
dA = (1 + p + s + yy ) 1 dA 0
1+ s
<ILQDOPHQWHWHQLHQGRHQFXHQWDTXHHO7HQVRUGH'HIRUPDFLyQHVVLPpWULFR
p
q
=
2
0
r p 0
0
s
q
2
dx = e at dX x = e at X + C
3DUD KDOODU OD )RUPD &DQyQLFD GH ODV (FXDFLRQHV GH 0RYLPLHQWR LPSRQHPRV x = X HQ t = 0 \
TXHGD
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Resueltas
x = e at X
y = Y
z = Z
SDUDHQFRQWUDUODYHORFLGDG~QLFDPHQWHVHWLHQHTXHGHULYDU
v = dx = aXe at = ax
x dt
v y = 0
v z = 0
(O7HQVRU9HORFLGDGGH'HIRUPDFLyQVHREWLHQHGHODVLJXLHQWHIRUPD
v
(d )ij = 1 i + j
2 x j x i
a 0 0
d = 0 0 0
0 0 0
&5'HPRVWUDUTXHSDUDWRGRSURFHVRGHGHIRUPDFLyQHQHOTXHORVGHVSOD]DPLHQWRV
HYROXFLRQDQ OLQHDOPHQWH FRQ HO WLHPSR (u = (X, t ) = f (X ) t ) HO WHQVRU GH GHIRUPDFLRQHV
LQILQLWHVLPDOHV SXHGHREWHQHUVHPHGLDQWHODOLQHDOL]DFLyQHQ t = 0 GHOWHQVRUPDWHULDOGH
GHIRUPDFLRQHVGHODIRUPD = (X, t ) = E (X,0 ) t
8WLOL]DUODH[SUHVLyQ E = F T d F
5HVROXFLyQ
3DUDHOWLHPSR t = 0 VHFXPSOHTXHHO7HQVRU*UDGLHQWHGHOD'HIRUPDFLyQVHFRQYLHUWHHQOD
LGHQWLGDG (F = ) \WHQLHQGRHQFXHQWDODH[SUHVLyQGHOHQXQFLDGRTXHGD
(X(x,0 ),0 ) = D(X(x,0 ),0 )
E
SRUORWDQWRUHVXOWD
( X,0) = ( X,0)
E
3RURWUDSDUWHWHQLHQGRHQFXHQWDODH[SUHVLyQGHORVGHVSOD]DPLHQWRV u
u j
u
ij = 1 i +
2 X j X i
f j
1 f i
=
+
2 X
X i
j
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Resueltas
f
j
TXHGDXQDH[SUHVLyQOLQHDOHQ t
\GHILQLHQGRXQWHQVRU A FRPR A ij = 1 i +
2 X j X i
= At
TXHDOGHVDUUROODUUHVXOWD
( X, t ) = ( X,0) + ( X,0) t
( X,0) = 0
( X,0)
( X,0) = E
( X,0) t
( X, t ) = E
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Problemas Resueltos
352%/(0$65(68(/726
356REUHHOWHWUDHGURGHODILJXUDVHSURGXFHXQDGHIRUPDFLyQXQLIRUPH (F = ctte ) FRQ
ODVVLJXLHQWHVFRQVHFXHQFLDV
/RVSXQWRV O A \ B QRVHPXHYHQ
(OYROXPHQGHOVyOLGRSDVDDVHU " p" YHFHVHOLQLFLDO
/DORQJLWXGGHOVHJPHQWR AC SDVDDVHU p / 2 YHFHVODLQLFLDO
(OiQJXOR AOC SDVDDVHUGH 45
z
&
$
6HSLGH
D -XVWLILFDUSRUTXpQRSXHGHXWLOL]DUVHODWHRUtDGHODGHIRUPDFLyQLQILQLWHVLPDO
E 2EWHQHUHOWHQVRUJUDGLHQWHGHODGHIRUPDFLyQORVSRVLEOHVYDORUHVGH " p" \HOFDPSR
GHGHVSOD]DPLHQWRVHQVXVIRUPDVPDWHULDO\HVSDFLDO
F 'LEXMDUHOVyOLGRGHIRUPDGR
5HVROXFLyQ
D (O iQJXOR AOC SDVD GH 90 D 45 SRU OR TXH HYLGHQWHPHQWH QR VH WUDWD GH XQD
SHTXHxDGHIRUPDFLyQ$GHPiVHQSHTXHxDVGHIRUPDFLRQHVVHFXPSOHTXH << 1 , \
HQHVWHSUREOHPD = / 4 0,7854
2EVHUYDFLyQODVGHIRUPDFLRQHVQRWLHQHQGLPHQVLRQHVHQLQJHQLHUtDVHDFRVWXPEUDD
FRQVLGHUDUFRPRSHTXHxDVGHIRUPDFLRQHVFXDQGRVRQGHORUGHQGH 10 -3 - 10 -4
E ,PSRQGUHPRVXQDSRUXQDODVFRQGLFLRQHVGHOHQXQFLDGR
3DUWLHQGRGHTXH F( X, t )= ctte \VDELHQGRTXH dx = F dX VHSXHGHLQWHJUDUHVWD~OWLPD
H[SUHVLyQ\VHREWLHQH
x = dx = FdX = F dX = F X + C
F11
F = F21
F31
F12
F22
F32
F13
F23
F33
C1
C = C 2
C 3
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Problemas Resueltos
FRQORTXHUHVXOWDQ 12 LQFyJQLWDV
6HLPSRQHQDKRUDODVFRQGLFLRQHVGHOHQXQFLDGR
(OSXQWR O QRVHPXHYH
(OSXQWR$QRVHPXHYH
(OSXQWR%QRVHPXHYH
0
0
0
0 = F 0 + C C = 0
0
0
0
a
a a F11
0 = F 0 = a F
21
0
0 a F31
0
0 a F12
a = F a = a F
22
0
0 a F32
F11 = 1
F21 = 0
F31 = 0
F12 = 0
F22 = 1
F32 = 0
$JUXSDQGRWRGDODLQIRUPDFLyQUHVXOWD
1 0 F13
F = 0 1 F23
0 0 F33
Vfinal = pVinicial
/D LJXDOGDG dVf = F dV0 SHUPLWH UHODFLRQDU ORFDOPHQWH YRO~PHQHV GLIHUHQFLDOHV HQ
GLIHUHQWHVLQVWDQWHVGHWLHPSR(QHVWHFDVR F HVFRQVWDQWHSDUDFDGD t ILMDGRSRUOR
TXHVHSXHGHLQWHJUDUODH[SUHVLyQ\VDFDUHOGHWHUPLQDQWHGH F IXHUDGHODLQWHJUDO
Vf = dVf =
V
3RUORTXHVHGHEHGHLPSRQHU
l AC,final =
p
2
F dV
V0
= F dV0 = F V0
V0
F = F33 = p
l AC,inicial
&RPR F HVFRQVWDQWHODWUDQVIRUPDFLyQHVOLQHDOHVGHFLUWUDQVIRUPDUHFWDVHQUHFWDV
SRU OR TXH AC HQ OD FRQILJXUDFLyQ GHIRUPDGD WDPELpQ VHUi XQ VHJPHQWR UHFWLOtQHR
$Vt
1 0 F13 0 aF13
x C = F X C = 0 1 F23 0 = aF23
0 0 p a ap
l AC, final = l AC = (aF13 , aF23 , ap) (a,0,0) = (a(F13 1), aF23 , ap ) =
2
= (a(F13 1) ) + (aF23 ) 2 + (ap) 2 = a (F13 1) 2 + F23
+ p2 =
2
p
2
l AC =
3RUORWDQWR
2
2
(F13 1) 2 + F23
+ p 2 = p (F13 1) 2 + F23
=0
p
2
2a = p a
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Problemas Resueltos
F13 = 1 ; F23 = 0
6HWLHQHDKRUD
1 0 1
F = 0 1 0
0 0 p
SRUORTXHIDOWDHQFRQWUDUHOYDORUGH p
AOC final = 45 = / 4
1 0 1 1 1
dX = (1,0,0) dx = F dX = 0 1 0 0 = 0
0 0 p 0 0
1 0 1 0 1
( 2)
( 2)
(2)
dX = (0,0,1) dx = F dX = 0 1 0 0 = 0
0 0 p 1 p
(1)
(1)
(1)
6HJXLGDPHQWHVHLPSRQHPRVTXH
(1)
( 2)
cos (AOC final ) = cos 45 = dx(1) dx ( 2 ) = 2
2
dx dx
dx (1) = 1 ,
dx ( 2 ) = 1 + p 2
dx (1) dx ( 2 ) = 1
3RUORWDQWR
1
1 + p2
3HUR F = p > 0 \SRUWDQWR p = 1
2
1
=
p = 1
2
2
(O7HQVRU*UDGLHQWHGHOD'HIRUPDFLyQHVSRUWDQWR
F =
/DVHFXDFLRQHVGHOPRYLPLHQWRVRQ x = F X
[ 1 0 1 ; ; + =
\ = 0 1 0 < = <
] 0 0 1 = =
(OFDPSRGHGHVSOD]DPLHQWRVHQGHVFULSFLyQPDWHULDOHV
<HQGHVFULSFLyQHVSDFLDOHV
Z
U(X, t ) = x X = 0
0
z
u(x, t ) = 0
0
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Problemas Resueltos
F 5HSUHVHQWDFLyQJUiILFDGHOWHWUDHGURGHIRUPDGR
z
&
&
2 D
% %
$ $
[
35(OVyOLGRGHODILJXUDVXIUHXQDGHIRUPDFLyQXQLIRUPH
6HSLGH
D 2EWHQHUODH[SUHVLyQJHQHUDOGHODGHVFULSFLyQPDWHULDOGHOFDPSRGHGHVSOD]DPLHQWRV
U(X, t) HQIXQFLyQGHOWHQVRUJUDGLHQWHPDWHULDOGHORVGHVSOD]DPLHQWRV J
E 2EWHQHUGLFKDH[SUHVLyQVDELHQGRTXHDGHPiVVHFXPSOHQODVVLJXLHQWHVFRQGLFLRQHV
GHFRQWRUQR
U Y = U Z = 0 , X, Y, Z
UX
X =0
=0
UX
X =L
, Y, Z
DFRWDQGRORVYDORUHVVLJQLILFDWLYRV
z
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Problemas Resueltos
5HVROXFLyQ
D 8QDGHIRUPDFLyQXQLIRUPHLPSOLFD
F( X, t) = F( t ) , t , X
6H TXLHUH REWHQHU OD H[SUHVLyQ JHQHUDO GH OD GHVFULSFLyQ PDWHULDO GHO FDPSR GH
GHVSOD]DPLHQWRV U(X, t ) HQIXQFLyQGHOWHQVRU*UDGLHQWH0DWHULDOGHORV'HVSOD]DPLHQWRV
J 3RUGHILQLFLyQVHWLHQHTXH
F = + J
3RU OR WDQWR VL F HV FRQVWDQWH VHJ~Q OD SRVLFLyQ J WDPELpQ OR VHUi 7HQLHQGR HQ
FXHQWDODGHILQLFLyQGH J HLQWHJUDQGRVXH[SUHVLyQTXHGDUiORVLJXLHQWH
J=
U( X, t )
X
dU = JdX
dU = J dX
dU = JdX
U=JX+C
GRQGH C HVXQDFRQVWDQWHGHLQWHJUDFLyQ3RUORWDQWRUHVXOWD
U(X) = J X + C
E 3DUWLHQGRGHOUHVXOWDGRDQWHULRU\DSOLFDQGRODVFRQGLFLRQHVGHFRQWRUQRGHOSUREOHPD
VHREWHQGUiHOYDORUGH J \ C
&RQGLFLRQHVGHFRQWRUQR
U Y = U Z = 0 , X, Y, Z /RVSXQWRV~QLFDPHQWHVHGHVSOD]DQHQODGLUHFFLyQ;
UX
X =0
= 0 , Y, Z (OSODQR<=HQHORULJHQHVILMR
UX
X =L
, Y, Z (VWHSODQRVHPXHYHHQODGLUHFFLyQ;XQLIRUPHPHQWH
(VFULELHQGR HQ FRPSRQHQWHV HO UHVXOWDGR GHO DSDUWDGR D VH YHUiQ PHMRU ODV
HFXDFLRQHV\ODVFRQFOXVLRQHVDODVTXHVHSXHGHOOHJDU
U X = J 11 X + J 12 Y + J 13 Z + C1
U Y = J 21 X + J 22 Y + J 23 Z + C 2
U Z = J 31 X + J 32 Y + J 33 Z + C 3
'HODFRQGLFLyQGHFRQWRUQR
U Y = 0 , X, Y, Z J 21 = J 22 = J 23 = C 2 = 0
U Z = 0 , X, Y, Z J 31 = J 32 = J 33 = C 3 = 0
'HODFRQGLFLyQGHFRQWRUQR
UX
X =0
= 0 , Y, Z J 12 = J 13 = C1 = 0
'HODFRQGLFLyQGHFRQWRUQR
UX
X =L
= , Y, J 11 L = J 11 =
L
)LQDOPHQWHTXHGD
'HVFULSFLyQ GH OD GHIRUPDFLyQ
/
L
J = 0
Problemas Resueltos
0 0
0 0
0 0
0
; C = 0
0
X
L
U ( X) = J X + C = 0
F 3DUD MXVWLILFDU ORV SRVLEOHV YDORUHV SRVLWLYRV \ QHJDWLYRV TXH SXHGH WRPDU OD
FRQGLFLyQDLPSRQHUHVODVLJXLHQWH
F >0
3RUORWDQWRVHGHEHFDOFXODUHOGHWHUPLQDQWHGH F
1 +
L
F = + J = 0
0 0
1 0
0 1
F =1+ > 0
L
> L
G 3DUDFDOFXODUORVWHQVRUHVHVSDFLDO\PDWHULDOGHODVGHIRUPDFLRQHVDVtFRPRHOWHQVRU
GHGHIRUPDFLRQHVLQILQLWHVLPDOHVVHSDUWHGHVXVUHVSHFWLYDVGHILQLFLRQHV
Tensor espacial de deformaciones :
Tensor material de deformaciones :
Tensor de deformaciones infinitesimales :
e = 1 ( F T F 1 )
2
E = 1 (F T F )
2
= 1 (J T + J )
2
L + 2 L2 / 1 + L
0
e=
0 0
0 0
0 0
+ 1 2
L 2 L2
; E=
0
0 0
0 0
0 0
; = 0
0 0
0 0
0 0
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Problemas Resueltos
H *UiILFRV E XX ; e xx ; x
L
E, e,
E XX
x
1
2
= 1
L
e xx
L
1
2
E XX HVXQDSDUiERODGHJUDGRTXHSDVDSRUHORULJHQFRQXQPtQLPRHQ
/ L = 1 E XX = 1 / 2
x HVODUHFWDLGHQWLGDGSHQGLHQWH\SDVDSRUHORULJHQ
e xx WLHQHGRVDVtQWRWDVXQDYHUWLFDOHQ / L = 1 \XQDKRUL]RQWDOHQ e xx = 1 / 2
6H FRQFOX\H SRU WDQWR TXH SDUD GHIRUPDFLRQHV / L SHTXHxDV ODV WUHV IXQFLRQHV
WLHQHQXQFRPSRUWDPLHQWRPX\VLPLODU\WLHQHQODPLVPDSHQGLHQWHHQHORULJHQ(V
GHFLUVHREWHQGUtDHOPLVPRUHVXOWDGRFRQFXDOTXLHUDGHODVGHILQLFLRQHVGHWHQVRUGH
GHIRUPDFLRQHV
6LQHPEDUJRIXHUDGHHVWHGRPLQLRJUDQGHVGHIRUPDFLRQHVVHYHTXHODVWUHVFXUYDV
VRQFODUDPHQWHGLIHUHQWHV
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Propuestas
&8(67,21(635238(67$6
&3-XVWLILFDUODVLJXLHQWHDILUPDFLyQ6yORVRQItVLFDPHQWHSRVLEOHVORVPRYLPLHQWRV
GHPHGLRVFRQWLQXRVHQORVTXHHOGHWHUPLQDQWHGHOWHQVRUJUDGLHQWHGHODGHIRUPDFLyQHV
SRVLWLYR
&3$SDUWLUGHODGHILQLFLyQGHORVWHQVRUHVPDWHULDO\HVSDFLDOGHGHIRUPDFLyQGHGXFLU
MXVWLILFDQGRWRGRVORVSDVRVODVH[SUHVLRQHVTXHUHODFLRQDQORVDODUJDPLHQWRVXQLWDULRVFRQ
GLFKRVWHQVRUHV
&3 'HILQLU HVWLUDPLHQWR HQ XQ SXQWR VHJ~Q XQD GLUHFFLyQ \ GHGXFLU VX H[SUHVLyQ
PDWHPiWLFDHQIXQFLyQGH
D (OWHQVRUPDWHULDOGHGHIRUPDFLyQ
E (OWHQVRUHVSDFLDOGHGHIRUPDFLyQ
F 9HULILFDUTXHSDUDHOFDVRGHSHTXHxDVGHIRUPDFLRQHVDPEDVH[SUHVLRQHVFRLQFLGHQ
&3'HGXFLUHQIXQFLyQGHOWHQVRUHVSDFLDOGHGHIRUPDFLRQHVODH[SUHVLyQGHOiQJXOR
TXH IRUPDQ HQ OD FRQILJXUDFLyQ PDWHULDO GRV VHJPHQWRV GLIHUHQFLDOHV TXH HQ OD
FRQILJXUDFLyQ HVSDFLDO HVWiQ RULHQWDGRV VHJ~Q GRV HMHV FRRUGHQDGRV FDUWHVLDQRV
GHIRUPDFLyQILQLWD
&3'HPRVWUDUODLJXDOGDG dV = det (F ) dV0
&3 ,QWHUSUHWDU ItVLFDPHQWH HO VLJQLILFDGR GH ODV FRPSRQHQWHV GHO WHQVRU GH
GHIRUPDFLRQHV PDWHULDO PHGLDQWH HO DQiOLVLV GH OD GHIRUPDFLyQ GH XQ SDUDOHOHStSHGR
HOHPHQWDO3DUWLFXODUL]DUORDOFDVRGHSHTXHxDVGHIRUPDFLRQHV
&3 (VFULELU OD H[SUHVLyQ GH OD GHVFRPSRVLFLyQ SRODU GHO WHQVRU JUDGLHQWH GH OD
GHIRUPDFLyQ SDUD HO FDVR GH GHIRUPDFLyQ ILQLWD \ GH GHIRUPDFLRQHV LQILQLWHVLPDOHV
LQGLFDQGRHOVLJQLILFDGRGHFDGDWpUPLQRHLQWHUSUHWiQGRORGHIRUPDJUiILFD
&3 $ SDUWLU GH OD GHVFRPSRVLFLyQ SRODU HQ SHTXHxDV GHIRUPDFLRQHV MXVWLILFDU HO
PRYLPLHQWRDOUHGHGRUGHXQDSDUWtFXODFRPRXQDGHIRUPDFLyQPiVXQDURWDFLyQ
&3'HGXFLUODH[SUHVLyQGHODGHIRUPDFLyQYROXPpWULFDSDUDORVFDVRVGHGHIRUPDFLyQ
ILQLWD\GHIRUPDFLyQLQILQLWHVLPDO
&3 'HGXFLU ODV H[SUHVLRQHV GH OD GHULYDGD PDWHULDO GH WHQVRU JUDGLHQWH GH OD
GHIRUPDFLyQ\GHOWHQVRUPDWHULDOGHGHIRUPDFLyQ
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Propuestas
&3-XVWLILFDUHQTXpFRQGLFLRQHVHOWHQVRUYHORFLGDGGHGHIRUPDFLyQFRLQFLGHFRQOD
WDVDGHYDULDFLyQGHULYDGDPDWHULDOGHOWHQVRUPDWHULDOGHGHIRUPDFLyQ
&33DUDHOFDPSRGHYHORFLGDGHV v = [ax, by, x ]T VHSLGH
D 2EWHQHUODVHFXDFLRQHVGHOPRYLPLHQWRHQIRUPDFDQyQLFDWRPDQGRFRPRFRQILJXUDFLyQGH
UHIHUHQFLDODFRUUHVSRQGLHQWHD t = 1
E /DVHFXDFLRQHVGHODVOtQHDVGHFRUULHQWH
F /DHYROXFLyQGHODGHQVLGDGDORODUJRGHOWLHPSRHQIXQFLyQGHVXYDORUHQHOLQVWDQWHGH
UHIHUHQFLD
&38QDGHIRUPDFLyQYLHQHGDGDSRUODVHFXDFLRQHV
x = X + 4Z
y = Y 2Z , z = Z X
+DOODUORVWHQVRUHVHVSDFLDO\PDWHULDOGHGHIRUPDFLyQ
&33DUDORVVLJXLHQWHVFDPSRVGHYHORFLGDG
2y
v 1 = x ,
, 3z
1 + t 1 + t 1 + t
v 2 = X , 2Y , 3Z
1 + t 1 + t 1 + t
6HSLGH
D 2EWHQHUODGHVFULSFLyQPDWHULDOGH v 1 \ODGHVFULSFLyQHVSDFLDOGH v 2 FRQILJXUDFLyQGH
UHIHUHQFLD t = 0
E /DGLVWULEXFLyQGHGHQVLGDGHVHQDPERVFDVRVGHQVLGDGLQLFLDO 0
F 3DUD HO FDPSR GH YHORFLGDGHV v 1 REWHQHU ODV GHVFULSFLRQHV PDWHULDO \ HVSDFLDO GHO
FDPSRGHGHVSOD]DPLHQWRVDVtFRPRGHORVWHQVRUHVGHGHIRUPDFLyQPDWHULDO*UHHQ
/DJUDQJH\HVSDFLDO$OPDQVL
G 5HSHWLUFSDUDFRQILJXUDFLRQHVPX\SUy[LPDVDODVGHUHIHUHQFLD (t 0 )
H 'HPRVWUDUTXHHQHVWHFDVRORVGRVWHQVRUHVGHGHIRUPDFLyQFRLQFLGHQ
&30HGLDQWHWUHVH[WHQVyPHWURVVHPLGHQORVDODUJDPLHQWRVXQLWDULRVHQXQSXQWR P
GHXQDSLH]DVRPHWLGDDWHQVLyQSODQDHQWUHVGLUHFFLRQHVIRUPDQGRiQJXORVGH 60 YHU
ILJXUD6HDQGLFKRVDODUJDPLHQWRV a b \ c 'HWHUPLQDUHOWHQVRUGHGHIRUPDFLRQHVHQHO
SXQWR P
\
F
E
3
D
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Propuestas
&36HDQODVHFXDFLRQHVGHOPRYLPLHQWR
x=X
y = Y + Z2 t
z = Z + Y2t
+DOODUODORQJLWXGHQHOLQVWDQWH t = 1 GHOVHJPHQWRTXHHQODFRQILJXUDFLyQGHUHIHUHQFLD
t = 0 HVUHFWR\XQHORVSXQWRV A(0,0,0) \ B(0,1,1)
&36HDQODVHFXDFLRQHVGHPRYLPLHQWRGHXQPHGLRFRQWLQXR
x = X + Yt
y=Y
z=Z
0
0
E = te
0
0
te tY
&36HDQODVHFXDFLRQHVGHPRYLPLHQWRGHXQPHGLRFRQWLQXR
x = X + Y 2 t
y=Y
z = Z
y=Y
z = Z Xt
&DOFXODUHOiQJXORTXHIRUPDEDQHQHOLQVWDQWH t = 0 ORVVHJPHQWRVGLIHUHQFLDOHVTXHHQHO
LQVWDQWH t = t VRQSDUDOHORVDORVHMHV x z
&36HDHOPRYLPLHQWR
x = X + te tX
y=Y
z=Z
(VWXGLDU FyPR YDUtD D OR ODUJR GHO WLHPSR XQ HOHPHQWR GLIHUHQFLDO GH iUHD TXH HQ t = 0
WRPDHOYDORUGH dS 0 \VXYHFWRUQRUPDOHV n 0 = 1 (1 ,1 , 1)
3
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Propuestas
&3 &DOFXODU HO WHQVRU GH SHTXHxDV GHIRUPDFLRQHV HQ XQ SXQWR VDELHQGR TXH XQ
WULiQJXORUHFWiQJXORGHODGRV 3ds 4ds \ 5ds VLWXDGRHQGLFKRSXQWRVXIUHXQDODUJDPLHQWR
GHXQ 3% GHODKLSRWHQXVDXQDGLVPLQXFLyQGHXQ 6% HQHOiQJXORUHFWR\XQDXPHQWR
GHO 4% HQVXiUHD
&3(OPHGLRFRQWLQXRGHODILJXUDHVWiVRPHWLGRDXQHVWDGRGHGHIRUPDFLyQSODQD
XQLIRUPH
u x = u x (x, y) , u y = u y (x, y) , u z = 0 , x, y, z
'HWHUPLQDUHOWHQVRUGHSHTXHxDVGHIRUPDFLRQHVVDELHQGRTXH
\
'
(OiUHDSDVDDVHU 1 + p YHFHVVXYDORULQLFLDO
(OVHJPHQWR AB SDVDWHQHUXQDORQJLWXG
1 + q YHFHVODLQLFLDO
&
(OVHJPHQWR AC QRVHGHIRUPD
$
%
&3(OHOHPHQWRGLIHUHQFLDOGHODILJXUDVHHVWLUDHQODGLUHFFLyQGHOHMH x XQDFDQWLGDG
= tdS , (0 < t < ) 'LEXMDUORVJUiILFRV E XX t e xx t \ xx t GRQGH E xx e xx \ xx
VRQ ODV FRPSRQHQWHV xx GH ORV WHQVRUHV PDWHULDO GH GHIRUPDFLRQHV HVSDFLDO GH
GHIRUPDFLRQHV \ GH SHTXHxDV GHIRUPDFLRQHV UHVSHFWLYDPHQWH $FRWDU ORV YDORUHV
VLJQLILFDWLYRV
\
W
G6
\
[
W W
WG6
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Cuestiones Propuestas
&36HDQODVHFXDFLRQHVGHPRYLPLHQWR
x = X Z (X, t)
y = Y
z = Z + (X, t)
-XVWLILFDU TXH WRGRV ORV SODQRV SHUSHQGLFXODUHV DO HMH X DQWHV GH OD GHIRUPDFLyQ
SHUPDQHFHQ SODQRV GHVSXpV GH OD PLVPD +LSyWHVLV GH 1DYLHU%HUQRXLOOL SDUD IOH[LyQ GH
YLJDV
&3'DGRVGRVFDPSRVGHGHVSOD]DPLHQWRV U (1) (X, t ) ; U (2 ) = (X, t ) \XQWHUFHURGHOD
IRUPD U (3) = U (1) + k (t )U (2 ) GHPRVWUDU TXH VXV UHVSHFWLYRV WHQVRUHV PDWHULDOHV GH
GHIRUPDFLyQILQLWD E (1) , E ( 2 ) , E (3) HLQILQLWHVLPDO (1) , (2 ) , (3) FXPSOHQ
(3 ) = (1) + k (2 )
E (3) E (1) + kE (2 )
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Problemas Propuestos
352%/(0$635238(6726
336HFRQVLGHUDQORVVLJXLHQWHVFDPSRVGHYHORFLGDGHVHQXQPHGLRFRQWLQXR
X
v = Y
Z
I
6HSLGH
x
, v = y
z
II
D &DOFXODUORVFRUUHVSRQGLHQWHVWHQVRUHVPDWHULDOHV\HVSDFLDOHVGHGHIRUPDFLyQ E I , E II
\ e I , e II
E 'HPRVWUDU TXH HQ XQ FLHUWR HQWRUQR PX\ SUy[LPR D OD FRQILJXUDFLyQ GH UHIHUHQFLD
(t 0 , t 0 ) ORVFXDWURWHQVRUHVGHGHIRUPDFLyQGHODSDUWDGRDFRLQFLGHQ\FDOFXODU
VXYDORU-XVWLILFDUHVWDFRLQFLGHQFLD
F 2EWHQHU ORV FRUUHVSRQGLHQWHV WHQVRUHV 9HORFLGDG GH 'HIRUPDFLyQ (d I , d II ) \ ODV
GHULYDGDVPDWHULDOHVGHORVWHQVRUHVGHGHIRUPDFLyQ (E I , E II , e I , e II )
G &DOFXODUHOYDORUGHORVVHLVWHQVRUHVGHODSDUWDGRFHQXQHQWRUQRPX\SUy[LPRDOD
FRQILJXUDFLyQGHUHIHUHQFLD (t 0 , t 0 ) \SDUDODSURSLDFRQILJXUDFLyQGHUHIHUHQFLD
(t = 0 ) -XVWLILFDUORVFDVRVHQTXHFRLQFLGHQ
338QPHGLRFRQWLQXRVHPXHYHPDQWHQLHQGRXQHVWDGRGHGHIRUPDFLyQXQLIRUPHGH
WDOPRGRTXHORVWUHVHMHVFRRUGHQDGRVVRQOtQHDVPDWHULDOHV
6HSLGH
D /DH[SUHVLyQPiVJHQHUDOGHODVHFXDFLRQHVGHPRYLPLHQWREDMRHVWDVFRQGLFLRQHV
E 6DELHQGRDKRUDTXHORVHVWLUDPLHQWRVVHJ~QODVWUHVGLUHFFLRQHVFRRUGHQDGDVVRQ
LJXDOHV\TXHODGHQVLGDGHYROXFLRQDGHODIRUPD = 0 e 3at REWHQHUODGHVFULSFLyQ
HVSDFLDOGHOFDPSRGHYHORFLGDGHV\GHOWHQVRU*UDGLHQWHGHOD'HIRUPDFLyQ
F 6XSRQLHQGRTXHHOWHQVRUGHWHQVLRQHVGH&DXFK\HVGHODIRUPD
K(t)x g(x, z, t) 0
- k(t)z 0
= h(y)
0
0
0
n= 1
3
1 WLHQHODGLUHFFLyQGHOHMH x \PyGXORXQLGDGSDUDWRGRLQVWDQWH
3
GHWLHPSRREWHQHUODVIXQFLRQHV g \ h \HOSDUiPHWUR K 0 DVtFRPRODVIXHU]DVGH
YROXPHQGHVFULSFLyQHVSDFLDOTXHJHQHUDQHOPRYLPLHQWR ( Tema 4 )
1
3
336REUHODGHVFULSFLyQPDWHULDOGHXQFLHUWRFDPSRGHGHVSOD]DPLHQWRV U(X, Y, Z) VH
WLHQHODVLJXLHQWHLQIRUPDFLyQ
(VOLQHDOHQ X, Y, Z
(VDQWLVLPpWULFRUHVSHFWRDOSODQR Y = 0 HV
GHFLUVHFXPSOHTXH
U (; < = ) = U(;< = ) , ; < =
Z, w
C
Y, v
D
A
D
X, u
6HSLGH
D 'HWHUPLQDU OD H[SUHVLyQ PiV JHQHUDO GH GLFKR FDPSR SDUD TXH VH FXPSODQ ODV
FRQGLFLRQHV\
E 'HWHUPLQDU OD H[SUHVLyQ GHO FDPSR U FXDQGR DGHPiV VH FXPSOH OD FRQGLFLyQ
2EWHQHU HO WHQVRU JUDGLHQWH GH OD GHIRUPDFLyQ \ HO WHQVRU PDWHULDO GH GHIRUPDFLyQ
'LEXMDUODGHIRUPDGDGHOHOHPHQWRGHODILJXUDDFRWDQGRORVYDORUHVVLJQLILFDWLYRV
F 'HWHUPLQDUODVGLUHFFLRQHVGHILQLGDVSRUVXVYHFWRUHVXQLWDULRV T SDUD ORV FXDOHV OD
GHIRUPDFLyQVHUHGXFHDXQHVWLUDPLHQWRQRKD\URWDFLyQ
127$6HFRQVLGHUDUiQGHIRUPDFLRQHVILQLWDVQRLQILQLWHVLPDOHV
33 6REUH HO VyOLGR GH OD ILJXUD VH SURGXFH XQD GHIRUPDFLyQ XQLIRUPH WDO TXH ORV
SXQWRV A B \ C QRVHPXHYHQ6XSRQLHQGRGHIRUPDFLyQLQILQLWHVLPDOVHSLGH
D ([SUHVDU HO FDPSR GH GHVSOD]DPLHQWRV
HQIXQFLyQGHYDORUHVJHQpULFRVGHODV
GHIRUPDFLRQHV\GHODVURWDFLRQHV
E ,GHQWLILFDU ODV FRPSRQHQWHV QXODV GHO
WHQVRU GH GHIRUPDFLRQHV \ H[SUHVDU HO
YHFWRU GH URWDFLRQHV HQ IXQFLyQ GH ODV
GHIRUPDFLRQHV
D
(
D
6HVDEHDGHPiVTXH
(O VHJPHQWR AE SDVD D PHGLU
(1 + p) YHFHVVXORQJLWXGLQLFLDO
(O YROXPHQ GHO VyOLGR SDVD D VHU
(1 + q) YHFHVVXYDORULQLFLDO
(O iQJXOR VH LQFUHPHQWD HQ XQ YDORU
r HQUDGLDQHV
%DMRHVWDVFRQGLFLRQHVVHSLGH
'
D
\
&
AB = AC = AD = a
EBC =
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Problemas Propuestos
F 2EWHQHU HO WHQVRU GH GHIRUPDFLRQHV HO YHFWRU GH URWDFLyQ \ HO FDPSR GH
GHVSOD]DPLHQWRVHQIXQFLyQGH p q \ r
127$/RVYDORUHVGH p q \ r VRQSHTXHxRV\SXHGHQGHVSUHFLDUVHLQILQLWpVLPRV GH
VHJXQGRRUGHQ
336REUHHOVyOLGRGHODILJXUDVHSURGXFHXQDGHIRUPDFLyQXQLIRUPHFRQODVVLJXLHQWHV
FRQVHFXHQFLDV
z
/RV HMHV x z VRQ GRV OtQHDV
PDWHULDOHV(OSXQWR A QRVHPXHYH
(OYROXPHQGHOVyOLGRQRYDUtD
(OiQJXOR xy QRYDUtD
(O iQJXOR yz VH LQFUHPHQWD HQ r
UDGLDQHV
(O VHJPHQWR AF SDVD D PHGLU
1 + p YHFHVVXORQJLWXGLQLFLDO
(O iUHD GHO WULiQJXOR ABE SDVD D VHU
1 + q YHFHVVXYDORULQLFLDO
D
$
\]
'
[\
%
&
6HSLGH
D
[
D ([SUHVDU HO FDPSR GH GHVSOD]DPLHQWRV HQ IXQFLyQ GH YDORUHV JHQpULFRV GH ODV
GHIRUPDFLRQHV\URWDFLRQHV
E ,GHQWLILFDUODVFRPSRQHQWHVQXODVGHOWHQVRUGHGHIRUPDFLRQHV\H[SUHVDUHOYHFWRUGH
URWDFLRQHVHQIXQFLyQGHODVFRPSRQHQWHVGHOWHQVRUGHGHIRUPDFLRQHV
F 2EWHQHU HO WHQVRU GH GHIRUPDFLRQHV HO YHFWRU GH URWDFLyQ \ HO FDPSR GH
GHVSOD]DPLHQWRVHQIXQFLyQGH p q \ r
127$/RVYDORUHVGH p q \ r VRQSHTXHxRV\SXHGHQGHVSUHFLDUVHORVLQILQLWpVLPRVGH
VHJXQGRRUGHQ
336REUHHOSDUDOHOHStSHGRGHODILJXUDVHSURGXFHXQDGHIRUPDFLyQXQLIRUPHTXHOR
OOHYDDGHIRUPDUVHFRPRVHLQGLFD6HSLGH
D 2EWHQHU HO WHQVRU PDWHULDO GH GHIRUPDFLyQ \ HO WHQVRU GH GHIRUPDFLyQ LQILQLWHVLPDO
UHGRQGHDUODFLIUDGHFLPDOHQDPERVFDVRV
E 2EVHUYDQGR ORV YDORUHV REWHQLGRV SDUD ORV GRV WHQVRUHV GH GHIRUPDFLyQ DQDOL]DU \
MXVWLILFDUVLODKLSyWHVLVGHSHTXHxDVGHIRUPDFLRQHVVHUtDDGHFXDGDSDUDHVWHFDVR(VWD
KLSyWHVLVVHVXSRQGUiFLHUWDSDUDORVVLJXLHQWHVDSDUWDGRVF\G
F 2EWHQHU HO FDPSR GH GHVSOD]DPLHQWRV \ HO YHFWRU GH URWDFLyQ VDELHQGR TXH HO SXQWR
O SHUPDQHFHILMRHOSXQWR A HVWiVREUHHOHMH x \HOSXQWR B HVWiVREUHHOSODQR xy
G 2EWHQHUHOYROXPHQGHOVyOLGRGHIRUPDGR
z
'
'
&
1 = 1.005
2
2 = 0.996
2
3 = 0.998
2
D' E' = 0.996 a
&
(
F
%
2
[
E
D
F
33 (QHO SDUDOHOHStSHGR GH OD ILJXUD HO SODQR A BC D VXIUH XQ GHVSOD]DPLHQWR GH
YDORU
U x = (t ) (t )Y ;
U y = (t )X
Uz = 0
Z
D
$
&
2
'HVFULSFLyQ GH OD GHIRUPDFLyQ
Problemas Propuestos
G 5HSHWLUHODSDUWDGRFXWLOL]DQGRODWHRUtDGHGHIRUPDFLyQLQILQLWHVLPDO'HPRVWUDUTXH
FXDQGR p \ q VRQSHTXHxRVORVUHVXOWDGRVGHORVDSDUWDGRVF\GFRLQFLGHQ
z
&
&
2
AA' = p
BB' = q
%
\
$
[
CC' = q
p>0
q>0
(FXDFLR Q HVGH
FRPSDWLELOLGDG
&8(67,21(65(68(/7$6
&5 (Q XQ FLHUWR LQVWDQWH GH XQ SURFHVR GH GHIRUPDFLyQ HO WHQVRU GH SHTXHxDV
GHIRUPDFLRQHVWLHQHHOVLJXLHQWHYDORU
y
8x
2
y
=
x
2
3 x 2z 0
2
3 x 2 z
x3
2EWHQHU HO YHFWRU GH GHVSOD]DPLHQWRV \ HO WHQVRU GH URWDFLyQ VDELHQGR TXH DPERV VH
DQXODQHQHOSXQWR (0,0,0)
5HVROXFLyQ
QLFDPHQWHVHKDQGHXWLOL]DUODVHFXDFLRQHVGHO)RUPXODULR%iVLFRHLQWHJUDUODVWHQLHQGRHQ
FXHQWDTXHVHHVWiEDMRFRQGLFLRQHVGHGHIRUPDFLyQLQILQLWHVLPDO
(OWHQVRUGHSHTXHxDVGHIRUPDFLRQHV YLHQH SDUWLFXODUL]DGR SDUD XQ LQVWDQWH GHWHUPLQDGR
GHWLHPSRHVGHFLU
y
8x
2
y
*
(x, t ) = (x ) =
x
2
3 x 2z 0
2
3 x 2 z
x3
1
1
=0 ;
= 0 1 = C1 (t)
y
z
2
2
3
3
;
=0 ;
= x 2 2 = & 2 (W ) x 2 z
2
2
\
]
3 3
3
3
;
=
= 0 3 = & 3 (W ) + y
2
2
\
]
3 x 2 z
2
3y
2
0
u 1
u 1
= 2y ;
=0
y
z
u 2
u 2
;
;
=x
=0
y
z
u 3
u 3
;
;
=0
= x3
y
z
u 1 = 4x 2 y 2 + C1' (t )
u 2 = xy + C '2 (t )
u 3 = x 3 z + C 3' (t )
6H KD GH LPSRQHU DKRUD TXH HO SXQWR RULJHQ QR VH GHVSOD]D $Vt VH REWLHQH WDPELpQ
IiFLOPHQWHTXHODVFRQVWDQWHV ( C i' )VRQQXODV
4x 2 y 2
u(x ) = xy
x 3z
3
0
1
2
1
0
3 y 3 x 2 z
2
2
3
0
0
(x ) = y
2
3 2
0
0
x z
2
&5 'HWHUPLQDU OD GHVFULSFLyQ HVSDFLDO GHO FDPSR GH YHORFLGDGHV FRUUHVSRQGLHQWH DO
VLJXLHQWHWHQVRUYHORFLGDGGHGHIRUPDFLyQ
t e t x
d(x, t )= 0
0
0
t e + 1
te y + 1
0
0
d
=
= te y
= = & (W ) + te y
[
\
]
2
2
2
=0 ;
=0 ;
= 0 2 = C 2 (t )
x
y
z
3
3
3
=0 ;
=0 ;
= 0 3 = C 3 (t )
x
y
z
6HLPSRQHQDKRUDODVFRQGLFLRQHVGHFRQWRUQR
t 1 t + C1
C1 = 1
para x = 0 : 0 = 0 = C 2
C2 = 0
C =0
0 C 3
3
(OUHVXOWDGRILQDOHV
9(&725 9(/2&,'$'
te y 1
(x, t ) = 0
0
v 1
v1
v1
= te tx ;
=0 ;
= 0 v1 = C1 (t ) + e tx
x
y
z
v 2
v 2
v 2
=0 ;
=0 ;
= 2 v 2 = C 2 (t ) + 2 z
x
y
z
v 3
v 3
v 3
=0 ;
= 2 te y ;
= 0 v 3 = C 3 (t ) + 2 te y
x
y
z
6HLPSRQHQDFRQWLQXDFLyQODVFRQGLFLRQHVGHFRQWRUQR
W + &
v 0 = = &
SDUD x = 0 :
W W + &
& = t- 1
& =
& = W
/DGHVFULSFLyQHVSDFLDOGHOFDPSRGHYHORFLGDGHVTXHGDILQDOPHQWH
e tx + t 1
v(x W )= 2 z
2 te y t
&8(67,21(635238(67$6
&3 ,QGLFDU HO SDSHO TXH MXHJDQ ODV HFXDFLRQHV GH FRPSDWLELOLGDG GH XQ FDPSR GH
GHIRUPDFLRQHV HQ OD 0HFiQLFD GH 0HGLRV &RQWLQXRV \ H[SOLFDU HVTXHPiWLFDPHQWH HO
SURFHVRGHREWHQFLyQGHGLFKDVHFXDFLRQHVSDUDHOFDVRGHGHIRUPDFLyQLQILQLWHVLPDO
&3(QXQFLHUWRLQVWDQWHGHXQSURFHVRGHGHIRUPDFLyQHOWHQVRUGHGHIRUPDFLRQHV
LQILQLWHVLPDOHVWLHQHHOVLJXLHQWHYDORU
y
8x
2
y
= x
2
3 2
2 x z 0
3 2
x z
2
3
x
2EWHQHU HO YHFWRU GH GHVSOD]DPLHQWRV \ HO WHQVRU GH URWDFLyQ VDELHQGR TXH DPERV VH
DQXODQHQHOSXQWR ( 0, 0, 0)
&3+DOODUHOFDPSRGHYHORFLGDGHVHQXQLQVWDQWHVDELHQGRTXHHQGLFKRLQVWDQWH t HO
WHQVRUYHORFLGDGGHGHIRUPDFLyQHV
e x + x e x
y
d =
y x
0 3y
3y 0
&3 'HGXFLU HO FDPSR GH YHORFLGDGHV FRUUHVSRQGLHQWH DO WHQVRU YHORFLGDG GH
GHIRUPDFLyQ
0
d(x, t )= t e t y
0
te t y
0
0
0
t e t z
&3 'HWHUPLQDU OD GHVFULSFLyQ HVSDFLDO GHO FDPSR GH YHORFLGDGHV FRUUHVSRQGLHQWH DO
WHQVRUYHORFLGDGGHGHIRUPDFLyQVLJXLHQWH
0
d(x, t )= 0
t e t z
0
te t y
0
te t z
0
0
SDUD
SDUD
=
=
=Y =
]
=
W , [ \
[ ]
&36HWLHQHHOWHQVRUYHORFLGDGGHGHIRUPDFLyQVLJXLHQWH
0 y 0
d(x, t) = y 1 0
0 0 t
2EWHQHU OD IRUPD FDQyQLFD GH ODV HFXDFLRQHV GHO PRYLPLHQWR VDELHQGR TXH HQ HO SXQWR
(0,0,0) HO FDPSR GH YHORFLGDGHV v (x, t ) \ HO WHQVRU YHORFLGDG GH URWDFLyQ w (x, t ) VH
DQXODQHQWRGRLQVWDQWH
61
4 Tensiones
Cuestiones Resueltas
4 Tensione s
CUESTIONES RESUELTAS
CR 4-1 Un medio continuo se mueve con un campo de velocidades cuya descripcin
espacial es: v(x, t ) = [z , x , y]T . El tensor de tensiones de Cauchy es de la forma:
y g(x, z, t) 0
(x, t ) = h(y) z(1 + t) 0
0
0
0
= T
=
x
C
0
y
C z(1 + t) 0 = [0 0 0]
z
0
0
0
x
0 1 0
v = [z x y] = 0 0 1
y
1
0
0
z
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
62
4 Tensiones
Cuestiones Resueltas
v v = [z x
0 1 0
y] 0 0 1 = [y z x ]
1 0 0
En un cierto plano que pasa por dicho punto las tensiones normal y tangencial son y ,
respectivamente. Razonar si son posibles los siguientes valores de y :
a) = 10 , = 1 .
b) = 5 , = 4 .
c) = 3 , = 1 .
Resolucin:
Se dibujar el Crculo de Mohr en 3-D para el estado tensional que se define y para los
puntos pedidos:
Punto b
Punto a
3 = 2
2 = 5
1 = 10
Punto c
La zona sombreada es la zona donde es posible encontrar puntos que representen estados
tensionales. nicamente hay que representar cada punto (estado tensional) y ver si cae en
dicha zona.
Se concluye que ningn estado tensional es posible.
Nota: en el caso c) el signo negativo de las tensiones tangenciales no se tiene en cuenta
porque cuando se representa en el Crculo de Mohr en 3-D se considera el valor absoluto.
Cosa que no sucede con las tensiones normales (traccin = positiva).
63
4 Tensiones
Cuestiones Resueltas
max = max = b
a-b
a+b
Crculo 3
Crculo 2
Crculo 1
2
1
2000-2004 by: X. Oliver & C. Agelet de 3Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
64
4 Tensiones
Cuestiones Resueltas
El crculo nmero:
1 : Corresponde a planos paralelos a 3.
2 : Corresponde a planos paralelos a 1.
3 : Corresponde a planos paralelos a 2.
Se tienen en cuenta ahora las siguientes consideraciones:
1. x es la tensin principal mayor:
=3
=2
2 = 3
3 = 9
1 = x = 1
=3
=2
2 = x = 1
3 = 1
1 = 5
1
2
45
ESTADO I
45
ESTADO II
Resolucin:
ESTADO III
65
4 Tensiones
Cuestiones Resueltas
Para poder sumar los dos estados, las tensiones deben actuar sobre los mismos planos.
Como los dos estados presentan orientaciones diferentes, deberemos buscar las tensiones
del estado II existentes sobre los planos del estado I. Para ello, representaremos el Crculo
de Mohr del estado II:
1
= 1
= 0
= 3
= 0
Plano b:
Plano a:
45
45
> 0
< 0
Plano c:
Polo
Plano horizontal
1
Plano b
Plano a
Plano vertical
Para realizar el crculo, se representan los planos a y b, ya que conocemos todas sus
coordenadas. Adems, como pertenecen al eje de abcisas, definen el dimetro del crculo,
por lo tanto, queda determinado.
El polo se encuentra como interseccin de lneas paralelas a los planos por los puntos que
los representan. Una vez obtenido, se pasa una lnea vertical, para una plano vertical, y una
horizontal, para un plano horizontal, por el polo, encontrando as sus representaciones.
Finalmente queda:
5
2
1
7
1
1
2
ESTADO I
2
2
ESTADO II
ESTADO III
66
4 Tensiones
Cuestiones Resueltas
E
D
45
45
C
*
A
figura a
figura b
Resolucin:
Primero se tiene que representar el crculo del estado tensional:
= > 0
= +
= 0
Plano inclinado:
Plano vertical:
= > 0
Plano vertical
Plano inclinado
Polo
Ahora tan slo hay que fijarse en los planos que, pasando por el polo, cumplen la
condicin pedida:
max
max
min
El caso que provoca mximas tensiones coincide con el apartado b), ya que corresponde al
caso de tensin principal 1 .
67
4 Tensiones
Cuestiones Propuestas
CUESTIONES PROPUESTAS
CP 4-1 Enunciar los siguientes postulados de la Mecnica de Medios Continuos:
a) Postulado de Cauchy.
b) Principio de accin y reaccin.
CP 4-2 Deducir la ecuacin de Cauchy a partir del principio de conservacin de la cantidad
de movimiento.
CP 4-3 Definir el tensor de tensiones esfrico, la tensin media y el tensor desviador y
escribir las relaciones entre ellos.
CP 4-4 Justificar la relacin entre las direcciones y las tensiones principales de un tensor de
tensiones y de su parte desviadora.
CP 4-5 Hallar la relacin entre el Crculo de Mohr correspondiente a un estado tensional
plano y el de su tensor desviador.
CP 4-6 Las componentes cartesianas del tensor de tensiones de Cauchy en un punto son:
1 0 1
= 0 1 a
1 b 0
Se pide:
xy 5y
=
5y 4x
(en MPa )
68
4 Tensiones
Cuestiones Propuestas
, (a, m cttes.)
2
______________________________________________________________________
CP 4-13 El estado tensional en un punto de un medio continuo es tal que x = 1 , xz = 0 ,
xy = 0 y y z . Se sabe adems que la tensin tangencial, sobre los planos paralelos al
eje x , toma un valor mximo * = 4 . Si max es la mxima tensin tangencial actuante
sobre cualquier plano obtener todas las posibles combinaciones de y y z para que:
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
69
4 Tensiones
Cuestiones Propuestas
a) max = 6
b) max = 4
CP 4-14 Del estado tensional en un punto de un medio continuo se sabe que la mxima
tensin tangencial en planos paralelos a 1 es 1max = 2 . Obtener todos los valores de 1 ,
2 y 3 que hacen posible el estado tensional = 2 , = 2 sobre un cierto plano, en los
siguientes casos (independientes entre s):
=2.
a) La mxima tensin tangencial en planos paralelos a 2 es max
2
max
b) La mxima tensin tangencial en planos paralelos a 3 es 3 = 0 .
c) La mxima tensin tangencial en planos paralelos a 2 es max
=4.
2
a) = 4 , = 2
b) = 4 , = 1
c) = 7 , = 0
2
*
*
______________________________________________________________________
1
CP 4-17 Determinar de forma grfica todos los posibles valores
*
*
de sabiendo que en un cierto plano las tensiones normal ( )
1
*
y tangencial ( ) son:
a) = 3 , = 1
2
b) = 1 , = 1
______________________________________________________________________
CP 4-18 La figura representa el estado tensional en un punto
del que se sabe que la tensin tangencial mxima sobre
cualquier plano es max = 2 . Dibujar los correspondientes
Crculos de Mohr en 3D. A dicho estado se le aade un estado
tensional hidrosttico ( = * 1). Determinar todos los posibles
valores de * para que sea posible, en un cierto plano, el
estado tensional = 6 , = 0 .
v > 0
*
1
1
1
70
4 Tensiones
Cuestiones Propuestas
y
P
60
4
4
71
4 Tensiones
Cuestiones Propuestas
4
CP 4-27 Determinar los valores de y para los que son posibles los siguientes estados
tensionales, siendo > 0 y = 0.5 :
(a)
(b)
(c)
1
72
4 Tensiones
Cuestiones Propuestas
b) = 45
**
3
2
CP 4-31 Un estado tensional plano es tal que, en el plano donde la tensin tangencial es
mxima positiva (criterio del Crculo de Mohr), la tensin normal viene dada por:
* = 10 5 *
3
, * 0
Se pide:
a) Obtener el valor de las tensiones principales en funcin de * y * .
b) Obtener el mximo valor que puede tomar una tensin principal en estas condiciones.
CP 4-32 Dados los estados tensionales planos I y II, determinar el valor de para que en
el estado III = I + II, ocurra que:
a) No haya compresiones en ningn plano.
b) No haya tracciones en ningn plano.
c) Sea un estado de cortante puro.
4
3
3
4
ESTADO I
ESTADO II
73
4 Tensiones
Cuestiones Propuestas
3
ESTADO II
ESTADO I
CP 4-34 Calcular los posibles valores de , , y para que el estado tensional III sea
la suma de I y II. ( 0 )
4
45
45
3
ESTADO I
ESTADO II
ESTADO III
60
a
60
*b
45
2 3
3
ESTADO I
ESTADO II
3
3
74
4 Tensiones
Cuestiones Propuestas
CP 4-36 Dados los estados tensionales I y II, determinar los posibles valores de y de
para que en el estado III = I + II se verifique que la tensin principal 2 sea positiva y
forme un ngulo de 30 con el eje y .
4
30
ESTADO I
ESTADO II
CP 4-37 Determinar todos los posibles valores de * para que en el estado tensional suma
de I y II se verifique:
a)
b)
c)
d)
e)
*
*
2
1
ESTADO I
*
*
ESTADO II
75
5 Ecuaciones de conservacin-balance
Cuestiones Resueltas
5 Ecuacio n es de
conservacin-balance
CUESTIONES RESUELTAS
CR 5-1 Justificar si son ciertas o no las siguientes afirmaciones:
a) El flujo de masa a travs de una superficie material cerrada slo es nulo si el
movimiento es estacionario.
b) El flujo de masa a travs de una superficie de control cerrada es nulo si este flujo es
estacionario.
Resolucin:
a) La afirmacin es falsa porque una superficie material siempre est formada por las
mismas partculas, y por tanto, no puede ser atravesada por ninguna partcula a lo largo
del movimiento. Por este motivo, el flujo de masa a travs de una superficie material es
siempre nulo, independientemente de si el movimiento es estacionario o no.
b) La afirmacin es cierta ya que la Ecuacin de Continuidad aplicada a un flujo estacionario
implica lo siguiente:
+ (v ) = 0
t
(v ) = 0
=0
Ec. de continuidad
Flujo estacionario
(v )dV =
V
v n dS = 0
76
5 Ecuaciones de conservacin-balance
Cuestiones Resueltas
a) La afirmacin es cierta ya que una superficie material siempre est formada por las
mismas partculas, y un flujo convectivo implicara que alguna debera atravesar la
superficie.
b) La afirmacin es falsa. Sera cierta si el material fuese incompresible, pero no se cumple
en general:
v dV = v n dS = 0
= ctte v = 0
v dV = 0
V
v n dS =
, V
v dV
V
v = 0 x V
v=0
d
= 0 (el fluido es incompresible)
dt
dt t
=0
t
v = 0
CR 5-4 Un gas a presin en un recipiente de masa M sale del mismo a una velocidad
estacionaria v por un conducto de seccin S . Justificar hacia dnde tender a moverse el
conjunto y con qu aceleracin lo har.
v
S
Resolucin:
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
77
5 Ecuaciones de conservacin-balance
Cuestiones Resueltas
Suponemos por hiptesis que el peso del gas es despreciable frente al del recipiente, y que
la presin atmosfrica tambin lo es frente a la presin interior.
La fuerza que actuar sobre el gas tendr la siguiente expresin si se tiene en cuenta la 2
Expresin del Teorema del Transporte de Reynolds:
F= d
dt
v dV =
V
(v ) dV +
t
v (n v ) dS =
vve dS = v 2 S e
donde:
n es la normal unitaria exterior a la superficie que encierra al gas.
e es el vector unitario en la direccin de la velocidad.
Por otro lado, la fuerza que acta sobre el gas es:
F = Fd / g + Fpresin atm = Fd / g
Segn el Principio de Accin y Reaccin la fuerza que realiza el gas sobre el depsito ser:
Fg/d = Fd/g = v 2 S e = Ma g / d
a g/d =
v 2 S
e
M
S
v
Resolucin:
Si se tiene en cuenta la 2 Expresin del Teorema del Transporte de Reynolds y que el rgimen es
estacionario, las fuerzas que actan sobre el fluido tienen la siguiente expresin:
Fext/f
d
v dV =
(v ) dV + v (n v) dS = v (n v) dS
V
V
V
S
t
dt
78
5 Ecuaciones de conservacin-balance
Cuestiones Resueltas
Seccin S2
Superficie Slat-1
Superficie Slat-3
Seccin S4
En primer lugar, hay que notar que el vector velocidad del fluido a lo largo de las
superficies S lat-1 y S lat -3 es perpendicular a la normal unitaria exterior al volumen que
encierra al fluido:
Slat 1
v n = 0
Slat 3
Lo mismo ocurre en la pared del disco.
En las secciones S 2 y S 4 los vectores v y n no son perpendiculares, pero como existe
simetra y v es perpendicular a F , no aportan componentes de fuerza horizontales.
As que nicamente quedar:
ext/f
ext/f
v (n v) dS =
ve(e ve) dS = v S e
2
Y ahora tan slo hay que despejar el valor del mdulo de la fuerza F :
F = v 2 S + pS
79
5 Ecuaciones de conservacin-balance
Problemas Resueltos
PROBLEMAS RESUELTOS
PR 5-1 Por la tubera de la figura circula, en rgimen estacionario, un caudal de agua Q
desde el extremo A (de seccin S A ) al extremo B (de seccin S B < S A ). La tubera est
sustentada en el punto O mediante el elemento rgido P - O .
Se pide:
a) Determinar las velocidades de entrada y salida v A y v B en funcin del caudal.
b) Determinar los valores del ngulo que hacen mximo y mnimo el mdulo de la
reaccin F en O , y los correspondientes valores de dicha reaccin.
c) Lo mismo que en el apartado b) pero respecto al momento M en O .
d) Determinar la potencia W de la bomba necesaria en cada caso para proporcionar el
caudal Q .
SB
eB
vB
R
pA
vA
pB
eA
P
SA
O
M
HIPTESIS:
1) El agua es un fluido perfecto ( ij = p ij ) e incompresible.
2) Se desprecia el peso de la tubera y del agua.
Resolucin:
a) Que el agua sea un fluido incompresible significa lo siguiente:
= ctte para una misma partcula
d
=0
dt
vdV = 0 , V
(1)
Se debe ahora establecer el volumen adecuado de integracin. Para ello hay que buscar
un volumen de control tal que su contorno sea una superficie cerrada (S = V ) , para
poder aplicar el Teorema de la Divergencia:
v dV =
n v dS
, V
(2)
n v dS = 0 , V
80
5 Ecuaciones de conservacin-balance
Problemas Resueltos
SA
nv dS =
nv dS +
SA
( e A ) v A e A dS +
SB
SB
nv dS =
e B v B e B dS =
= - v A SA + v B SB = 0
v A SA = v B SB = Q
Q
SA
(3)
Q
vB =
SB
bdV +
t dS = d
dt
v dV
(4)
v dV =
t
v dV +
v(n v ) dS
(5)
SA
SA
) dS +
v(n v) dS +
v A e A (e A v A e A
SB
v(n v)dS =
SB
v B e B (e B v B e B ) dS
R = v 2A S A e A + v 2B S B e B
R = Q 2 1 e A + 1 e B
SB
SA
(6)
81
5 Ecuaciones de conservacin-balance
Problemas Resueltos
Se han de analizar ahora las diferentes fuerzas que forman parte de R . Al despreciar
las fuerzas msicas (b = 0 ) , slo se tienen en cuenta las fuerzas superficiales, es decir,
las aplicadas en el contorno del volumen de control ( S A , S B y S lat , sta ltima,
superficie lateral de las paredes):
R=
=
=
SA
SA
bdV +
t dS +
p A e A dS +
SB
t dS =
t dS +
SB
Slat
t dS
tdS =
p B (e B )dS + R t/f
Aqu R t/f representa las fuerzas ejercidas sobre el fluido por las paredes de la tubera,
que en principio desconocemos, pero que utilizando la expresin ( 6 ) pueden
obtenerse:
R t/f = R
SA
p A e A dS
SB
p B ( e B )dS
R t/f = v 2A S A e A + v 2B S B e B p A S A e A + p B S B e B
R t/f = (v 2A + p A )S A e A +(v 2B + p B )S B e B
(7)
Q2
R t/f =
+ p ASA e A +
+ p B S B e B
SA
SB
Ahora habr que pensar qu relacin tiene R t/f con nuestra incgnita F . Para ello
hay que tener en cuenta el Principio de Accin y Reaccin, y considerar la tubera y el
elemento rgido P - O como un slido nico. En estas condiciones, la fuerza que ejerce
el fluido sobre la tubera es:
R f/t = R t/f
Como es la nica accin que acta sobre el slido, y teniendo en cuenta que el peso de
la tubera es despreciable, esta fuerza deber ser compensada por la accin exterior F
para llegar al equilibrio:
R f/t + F = 0
F = R f/t = R t/f
Y por tanto, se obtiene finalmente el valor de F segn ( 7 ):
F = (v 2A + p A )S A e A +(v 2B + p B )S B e B
Q2
F =
+ p A S A e A +
+ p BS B e B
S
S
A
(8)
82
5 Ecuaciones de conservacin-balance
Problemas Resueltos
- eA
- FA eA
eB
FB eB
F mnima = Q 2 1 1 + p B S B p A S A
SB SA
Valor mximo de F :
= 3
2
F mxima = Q 2 1 1 + p B S B + p A S A
SB SA
c) Para encontrar el valor del momento M , respecto del punto O , hay que aplicar la
ecuacin de Balance del Momento de la Cantidad de Movimiento:
M liq =
r bdV +
r tdS = d
dt
r vdV
(9)
donde M liq es el momento resultante de los momentos que actan sobre el fluido.
Por otro lado, desarrollando el trmino de la derecha segn la segunda expresin del
Teorema del Transporte de Reynolds, se obtiene la siguiente expresin:
d
dt
r vdV =
t
r vdV +
(r v)(n v)dS
( 10 )
Igual que en el apartado b), se tiene un rgimen estacionario, y por tanto, la derivada
local es nula. Tambin aqu, hay que tener en cuenta que n y v son perpendiculares en
las paredes de la tubera. De esta forma, segn las expresiones ( 9 ) y ( 10 ), se obtiene:
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
83
5 Ecuaciones de conservacin-balance
M liq =
(r v)(n v)dS +
SA
SB
Problemas Resueltos
(r v)(n v)dS
( 11 )
SB
R v B (e Z )(e B v B e B )dS
M liq = v 2B RS B e Z =
Q2
Re Z
SB
( 12 )
El siguiente paso consiste en valorar las contribuciones de las fuerzas msicas, que en
este caso son nulas (b = 0 ) , y de las fuerzas superficiales:
M liq =
=
SA
r bdV +
r tdS +
=0+
SB
SB
r tdS =
r tdS +
Slat
r tdS =
r tdS =
R p B e Z dS + M t/f =
= R p B S B e Z + M t/f
donde M t / f es el momento que ejerce la tubera sobre el fluido. Para hallar su valor
se utilizar la expresin ( 12 ):
M t/f = M liq R p B S B e Z =
= v 2B RS B e Z R p B S B e Z
M t/f = RS B (v 2B + p B )e Z = R (
Q2
+ p B S B )e Z
SB
( 13 )
84
5 Ecuaciones de conservacin-balance
Problemas Resueltos
+ p B S B e Z
M = RS B (v 2B + p B )e Z = R
SB
Se puede comprobar que este resultado no depende del ngulo , y por tanto, su
mdulo tendr un valor constante.
1 v 2 dV +
2
: d dV
( 14 )
: d dV = 0
v x v x v x
x
y
z
v y v y v y
= pTr (l) = pTr
=
x
y
z
v
v z v z
z
y
z
x
v x v y v z
= p v = 0
+
+
y
z
x
= p
1 v 2 dV =
2
t
1 v 2 dV +
2
1 v 2 (n v)dS
2
SA
SA
1 v 2 (n v)dS +
1 v 2 (n v)dS =
SB 2
2
1 v 2 ( v )dS +
1 v 2 ( v )dS =
A
A
B
B
SB 2
2
85
5 Ecuaciones de conservacin-balance
Problemas Resueltos
1
1
= v 3A S A + v 3BS B
2
2
Utilizando la expresin ( 3 ), se obtiene el resultado final:
W = 1 Q 3 12 12
2
SB SA
pr , vr
pR , vR
R
r
pR
HIPTESIS:
1) Fluido perfecto: = p1
2) Suponer una distribucin radial de velocidades (v = -v r e r ) .
3) Densidad del fluido: = ctte .
4) Fuerzas msicas despreciables.
5) El problema debe resolverse utilizando exclusivamente formas globales (integrales) de
las ecuaciones de conservacin-balance.
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
86
5 Ecuaciones de conservacin-balance
Problemas Resueltos
Resolucin:
a) Para encontrar la velocidad hay que utilizar el Principio de Conservacin de la Masa, que en
forma global tiene la siguiente expresin:
dM = d
dt dt
dV = 0
V
dV = t dV +
V
v ndS
dV = 0
V t
dV =
V
v ndS = 0
V
r r
v dV
2
87
5 Ecuaciones de conservacin-balance
K=1
2
v 2r dV = 1
V
2
R& 2 (t)
Problemas Resueltos
R 4 (t)
dV =
r4
r 2 dr
R(t ) r 4
b vdV +
t vdS = d
dt
1 v 2 dV + : ddV
V
2
&
Pe = t vdS = K
V
Sabiendo que los fluidos nicamente sufren tensiones de compresin y que p es una
constante por hiptesis, se desarrolla el primer trmino:
Pe = t vdS =
V
(p r (e r )) ( v r e r )dS =
R 2 (t)
= p R R& (t)4R 2 (t) + lim{p r v r 4r 2 }= p R R& (t )4R 2 (t) + lim p r R& (t) 2 4r 2
r
r
r
2
2
&
&
Pe = p R R(t)4R (t) p R(t)R (t)4
d (R
&& 2 (t)R 3 (t))
dt
1
p R (t) = p +
2
R& (t)R 2 (t)
&& (t ) + 3 R& 2 (t )
p R (t) = p + R (t )R
88
5 Ecuaciones de conservacin-balance
donde
Entonces se tiene:
Problemas Resueltos
&& (t ) + 3 R& 2 (t )
p = R (t )R
2
d R 3 (t) = 3R 2 (t)R& (t )
dt
p d R 3 (t) = 3 d (R& 2 (t)R 3 (t))
dt
2 dt
Para encontrar el instante en el que desaparece la burbuja, nicamente hay que hallar
R(t) e imponer R(t = t final ) = 0 .
89
5 Ecuaciones de conservacin-balance
Cuestiones Propuestas
CUESTIONES PROPUESTAS
CP 5-1 Dada la descripcin espacial de una cierta propiedad por unidad de volumen de
las partculas, interpretar fsicamente el valor de:
a)
d
dt
b)
c) d
dt
d)
dV
dV
V
______________________________________________________________________
CP 5-2 Indicar e interpretar fsicamente el significado de las siguientes expresiones en la
Mecnica de Medios Continuos:
a)
(x( X, t + t) , t + t) - (x( X, t) , t)
t 0
t
lim
(x, t + t) - (x, t)
t 0
t
b) lim
c)
lim
Vt + t
(x, t + t )dV
Vt
(x, t)dV
t 0
d) lim
Vt
(x, t + t )dV
Vt
(x, t)dV
t 0
______________________________________________________________________
CP 5-3 Interpretar fsicamente los siguientes conceptos:
a) Derivada material y local de una propiedad.
b) Derivada material y local de una integral de volumen.
c) Justificar por qu motivos podra interpretarse que
v ndS
es la derivada
______________________________________________________________________
CP 5-4 Deducir las relaciones entre:
a) La derivada local y la derivada material de una magnitud .
b) La derivada local y material de una integral dV , cuando todo el flujo de la
V
90
5 Ecuaciones de conservacin-balance
Cuestiones Propuestas
91
5 Ecuaciones de conservacin-balance
Cuestiones Propuestas
______________________________________________________________________
CP 5-20 Definir:
a) Variables termodinmicas
b) Variables de estado
c) Procesos termodinmicos reversibles e irreversibles.
Indicar la utilidad del 2 principio de la termodinmica en la Mecnica de Medios
Continuos.
______________________________________________________________________
CP 5-21 Para un estado tensional de presin hidrosttica ( ij = p ij ) :
a) Demostrar que la potencia tensional es
p d
dV .
dt
V
93
5 Ecuaciones de conservacin-balance
Problemas Propuestos
PROBLEMAS PROPUESTOS
PP 5-1 El campo de velocidades de un fluido incompresible en rgimen estacionario es:
v(x ) = [- x - y; x - y; f(z)]T
Se pide:
a)
b)
c)
d)
Qa = 3Q
va
Sa = S1
Pb =P1
h
Qb = Q
vb
Sb = S1
Pb =P1
HIPTESIS:
1) Peso del fluido y de la tubera despreciables.
94
5 Ecuaciones de conservacin-balance
Problemas Propuestos
PP 5-3 En el aspersor para riego cuya planta y alzado se presentan en la figura entra a
travs de la seccin C un caudal de agua Q a una presin P y sale por los extremos A y
B a la presin atmosfrica Patm . El flujo se supone estacionario.
Se pide:
a) Las velocidades de entrada y de salida. (Justificar la frmula empleada).
b) La resultante en el punto O de las fuerzas ejercidas sobre el fluido del interior del
aspersor (fuerzas y momentos).
c) Las reacciones a ejercer sobre dicho punto O , para evitar que el aspersor se desplace
verticalmente.
d) Suponiendo que I 0 e I1 son, respectivamente, los momentos centrales de inercia del
aspersor en vaco y del aspersor lleno de agua, respecto del punto O , calcular la
aceleracin angular de giro del aspersor .
e) Si W * es la potencia de la bomba necesaria para proporcionar el caudal Q , calcular la
potencia necesaria para proporcionar un caudal 2Q .
z
Q/2
Q/2
Q/2
B
vB
vA
Q vC
PC = P
SC = S*
h/2
h/2
PB = Patm
SB = S
vB
A
PA = Patm
SA = S
vA
Q/2
ALZADO
PLANTA
HIPTESIS:
1) Fluido incompresible
2) Fluido perfecto: ij = p ij
3) Se desprecia el peso del aspersor y del agua de su interior.
4) S A = S B = S
5) S C = S*
6) m = I
95
5 Ecuaciones de conservacin-balance
Problemas Propuestos
c) La relacin entre S1 y S 2 (en funcin del parmetro ) necesaria para que la turbina
no gire.
y
R
v1
p= patm= 0
S1
1
Q
v2
S2
2
Q
S2
R
v2
HIPTESIS:
1) Fuerzas msicas despreciables.
2) Rgimen estacionario.
3) Fluido incompresible.
4) Fluido perfecto: ij = p ij
v1
S1
B
v1
z
e 3
v2
e$ 2
p = patm 0
h
h
h/2
h/2
e$ 1
x
96
5 Ecuaciones de conservacin-balance
Problemas Propuestos
HIPTESIS:
1) Rgimen estacionario
2) Fluido incompresible
3) Las presiones sobre las paredes laterales de la tubera se suponen constantes e iguales a
la presin de entrada p .
4) Fluido perfecto: ij = p ij
5) Se desprecia el peso del fluido y de la compuerta.
a
A
P2
a/2
a/2
P1
O
P2
v2
v1
P1
Detalle
HIPTESIS:
1) Se desprecian las fuerzas msicas en el fluido.
2) Fluido perfecto: ij = p ij .
Considere el anlisis por metro lineal.
97
6 Elasticidad lineal
Cuestiones Resueltas
6 Elasticid a d lineal
CUESTIONES RESUELTAS
CR 6-1 Justificar si son ciertas o no las siguientes afirmaciones:
a) En un material termoelstico, los trminos isentrpico y adiabtico son equivalentes.
b) Para materiales termoelsticos lineales, la segunda analoga trmica es aplicable siempre.
Resolucin:
a) Segn el 2 Principio de la Termodinmica, tenemos la siguiente expresin:
i
0 s&loc
= 0 s& ( 0 r q) 0
Para un material termoelstico siempre tendremos procesos reversibles y por tanto esta
desigualdad se convierte en igualdad:
i
0 s&loc
= 0 s& ( 0 r q) = 0
(1)
98
6 Elasticidad lineal
Cuestiones Resueltas
t
R
a cos
;
2
E zz = E z = E rz = 0 ; a << 1
Calcular el vector de tracciones t , en el contorno, en
E rr = E = a sin
E r =
componentes cilndricas
HIPTESIS:
1) , son las constantes de Lam.
Resolucin:
El Tensor de Deformacin coincidir con el Tensor de Deformacin Infinitesimal y tendr la
siguiente expresin en coordenadas cilndricas:
a sin
a cos
E(r, , z ) = (r, , z ) =
2
0
a cos
2
a sin
0
cos
0
2( + ) sin
= a cos
2( + ) sin
0
0
0
2 sin
El vector normal y exterior al cilindro tendr la siguiente expresin en cilndricas:
n = (1 , 0 , 0 )T
99
6 Elasticidad lineal
Cuestiones Resueltas
Resolucin:
Para resolver este problema se aplicar la 2 Analoga Trmica. Se tiene un problema I que se
descompondr en suma de dos problemas: II y III. Como no hay coacciones exteriores el
contorno u = {} , = V y t * = 0 .
PROBLEMA I:
b = 0
Acciones : t = 0 , en V
I =
Respuesta : u, ,
PROBLEMA II:
b II = b = 0
Acciones : t II = t = 0 , en V
II = 0
u II = 0 , II = 0 , II = 0
Respuesta :
PROBLEMA III:
b III = 0
Acciones : t III = 0 , en V
=
III
Respuesta :
z + 31 x + 32 y + C 3
t
CR 6-4 Un slido elstico, lineal e istropo, est sometido en todo su contorno exterior a
una presin constante,de valor , junto con un incremento de temperatura en su interior
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
100
6 Elasticidad lineal
Cuestiones Resueltas
u = 0 , en u
Respuesta :
u,,
PROBLEMA III:
Se resuelve primero el problema III, cuya solucin es trivial:
b
III
Acciones : t III
u
III
III
= 1 ( )
= n , en
Respuesta :
= 0 , en u
=
PROBLEMA II:
b
II
Acciones : t II
u
II
II
= 1 ()
= (p + )n , en
= u = 0 , en u
=0
101
6 Elasticidad lineal
Cuestiones Resueltas
( I)
+u
(II)
en
en u
Por hiptesis, se supondr que la solucin es la suma de las soluciones de los dos sistemas:
u = u (I) + u (II)
= (I) + (II)
= (I) + (II)
ECUACIONES GEOMTRICAS:
1 (u (II) + u (II) ) = 1 [(u (I) + u (II) ) + (u (I) + u (II) )] = 1 [u + u ]
ij = ij(I) + ij(II) = 1 (u i,(I)j + u (I)
j,i
j,i ) +
i, j
j,i
i, j
i, j
j,i
j,i
i, j
2
2
2
2
ECUACIONES CONSTITUTIVAS:
(I)
(II)
(II)
ij = ij(I) + ij(II) = [C ijkl (I)
]=
kl + ij ] + [C ijkl kl + ij
(II)
(I )
= C ijkl [ (I)
+ (II) ] = C ijkl kl + ij
kl + kl ] + ij [
Se comprueba, por tanto, que se cumplen todas las condiciones y que esta es, entonces, la
solucin del problema. Se concluye que el principio de superposicin tambin funciona en
termoelasticidad.
103
6 Elasticidad lineal
Problemas Resueltos
PROBLEMAS RESUELTOS
PR 6-1 Una corona cilndrica de altura h , radio interior R y radio exterior 2R , se
encuentra en el interior de un cilindro de rigidez infinita de altura h y radio 2R + a , con
a << R . Se somete la corona a un campo uniforme de temperatura . Determinar:
a) El valor de para el que las paredes laterales de la corona y el cilindro entran en
contacto.
b) Dibujar, acotando los valores significativos, la curva - , siendo el alargamiento
del radio interior de la corona. Determinar el valor de para el cual dicho radio
recupera su valor original.
c) Dibujar, acotando valores significativos, las curvas rr - , - , zz - , en los
puntos A y B .
a
2R
A
a
h
HIPTESIS:
1) Mdulo de Elasticidad: E
2) Coeficiente de Poisson: = 0
3) Coeficiente de dilatacin trmica :
4) Material elstico lineal e istropo.
5) No se considera el peso propio.
6) No se considera el rozamiento entre paredes.
2R
Resolucin:
En el problema se pueden distinguir dos fases distintas:
1 Fase
104
6 Elasticidad lineal
Problemas Resueltos
2 Fase
tanto resulta:
b III = 1 ( ) = 0
(1)
105
6 Elasticidad lineal
Problemas Resueltos
Adems, se impondr que u z (z) = 0 en todos los puntos, ya que no se dispone de ningn
dato sobre Ssuperior y S inferior de la corona. Para estas caras no se puede plantear ninguna
condicin de contorno en desplazamientos y de aqu no se pueden determinar las
constantes de integracin de u z en caso de suponer u z 0 .As que finalmente se adopta:
u(r ) = [u r (r) 0 0]
Hay que tener en cuenta que se est trabajando en coordenadas cilndricas y por tanto la
ecuacin anterior hay que adaptarla a este sistema de coordenadas. Previamente, y partir de
los datos, hay que establecer el valor de los parmetros que intervienen en el problema
(E, , = 0 ) :
=0
2( + )
E
=
=G = E
2(1 + )
2
= E = E
1 2
=
(2)
E [ 1 ( r u r )] = 0
r r r
(3)
u r = Ar + B
r
u(r ) = Ar + B 0 0
106
6 Elasticidad lineal
Problemas Resueltos
A B
r2
(r ) = 0
0
0
A + B2
r
0
0
0
(4)
Por ltimo, para obtener las tensiones se usa la ecuacin constitutiva de un material elstico
lineal istropo, particularizado con las expresiones ( 2 ):
= (Tr ())1 + 2
= E
(5)
1 Fase
(6)
CONDICIN DE CONTORNO EN r = R
Si r = R y segn las condiciones de contorno en el problema anlogo se tiene
que:
t = n = t * + n = t * + n
107
6 Elasticidad lineal
Problemas Resueltos
A
Teniendo en cuenta ( 6 ) y ( 7 ):
B
=
R2
(7)
A =
B=0
(8)
0
0
0
II = 0
0
0
0
(9)
0
0
E
0
E 0
II =
0
0
0
0
0
= 0
0
0
0
0
( 10 )
0
0 0
= 0 0
0
0 0 E
Para obtener el valor de para el cual las paredes laterales del cilindro y la corona entran
en contacto, basta imponer:
u r (r = 2R) = a
2R = a
= a
2R
( 11 )
Se volver a usar la misma geometra que en el problema inicial, pero ahora existe contacto
entre la corona y el cilindro. Es como resolver un problema aparte, con la misma geometra
pero con diferentes condiciones de contorno.
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
108
6 Elasticidad lineal
Problemas Resueltos
En este apartado se prescindir de la primera fase, pero recordando que el slido parte de
un estado que es el resultado de la anterior fase, es decir, ha tenido unos desplazamientos,
deformaciones, tensiones e incrementos de temperatura. Se utilizar la variable .
Se aplica la analoga trmica, igual que en el apartado anterior. El problema III queda igual
y por tanto tambin su resultado ( 1 ).Por tanto hay que resolver el problema II con las
mismas expresiones ( 3 ), ( 4 ) y ( 5 ):
(Aqu tambin se omiten los subndices II, aunque hay que volver a remarcar que todo est
referido al problema II y no al original)
Se tienen dos incgnitas a determinar: A , B . Por ello se deben establecer dos condiciones
de contorno. Se pueden imponer las tensiones en la pared lateral interna de la corona y los
desplazamientos de la exterior:
CONDICIN DE CONTORNO EN r = 2 R
Si r=2R y segn las condiciones de contorno u en el problema anlogo se tiene:
u r (r = 2 R) = 0
( 12 )
CONDICIN DE CONTORNO EN r = R
Si r = R se tiene segn las condiciones de contorno en el problema anlogo:
t = n = t * + n = t * + n
B
=
R2
( 13 )
Teniendo en cuenta ( 12 ) y ( 13 ):
A = 1
5
B = 4 R 2
5
( 14 )
109
6 Elasticidad lineal
Problemas Resueltos
III
III
1 (1 + 4 R 2 )
0
0
5
r2
2
1 (1 4 R ) 0
=
0
5
r2
0
0
0
1 E (1 + 4 R 2 )
0
0
2
5
2
1 E (1 4 R ) 0
0
=
5
r2
0
0
0
( 15 )
1 (1 + 4 R 2 )
0
0
2
5
2
1 (1 4 R ) 0
=
0
5
r2
0
0
0
4 E (1 + R 2 )
5
r2
=
0
( 16 )
2
4 E (1 R )
0
5
r2
0
E
Hay que tener en cuenta que no se ha considerado para nada que se ha partido del estado
final de la fase 1 y no de un estado neutro como se ha supuesto en este apartado.
En realidad se parte del estado con movimientos, deformaciones, tensiones e incremento
de temperatura correspondientes a = :
u inicial = u 1 Fase ( = )
inicial = 1 Fase ( = )
( 17 )
inicial = 1 Fase ( = )
( 18 )
110
6 Elasticidad lineal
Problemas Resueltos
u 2 Fase = u inicial + u()
2 Fase = inicial + ()
2 Fase = inicial + ()
Por lo tanto, para determinar , bastar con imponer que el movimiento del radio
interior, segn la 1 fase, sea igual pero de signo contrario al de la 2 fase, para que as el
movimiento total sea nulo.
1 fase: desplazamiento para r = R y = . Segn ( 10 ) y ( 11 ) se obtiene:
1 = u r (r = R , = ) = R = a
2
( 19 )
= 5a
6R
( 20 )
= a
2
3 R
5
111
6 Elasticidad lineal
Problemas Resueltos
rr
punto B:
rr = 0
3 E
5
punto A
rr = 1 E a
2 R
punto A
= 5 E a
6 R
punto B
=4E a
3 R
8 E
5
zz
punto B y punto A
zz = 4 E a
3 R
112
6 Elasticidad lineal
Problemas Resueltos
y
2
; yz =
x
2
HIPTESIS:
1) Propiedades del material: E , G .
x
COMENTARIO:
Se va a demostrar que la solucin de Coulomb en el caso de una pieza empotrada con un
momento torsor en el otro extremo no es una solucin simplificada como la mayora de las
de Resistencia de Materiales, sino que se trata de una solucin exacta.
Resolucin:
a) Para obtener el campo de desplazamientos debemos integrar el campo de
deformaciones dado en el problema.
Primero se encuentra el Vector de Rotacin :
1 E xz E xy
=
x
y
z
2
1 E yz E yy
=
= 0 0 = 0 1 no depende de y ( 1 1 ( y ))
y
y
z
1 E zz E yz
=
= 0 0 = 0 1 no depende de z ( 1 1 ( z ))
z
y
z
1 =
x + C1
2
113
6 Elasticidad lineal
2 E xx E xz
=
= 0 0 = 0 2 2 ( x)
x
z
x
2 E xy E yz
=0 =
y
z
x
2
2
2 E xz E zz
=
= 0 0 = 0 2 2 ( z)
z
z
x
2 =
y + C2
2
3 E xy E xx
=
= 0 0 = 0 3 3 ( x)
x
x
y
3 E yy E xy
=
= 0 0 = 0 3 3 ( y)
y
x
y
3 E yz E xz
=
= + =
z
x
y
2 2
3 = z + C3
= E xz + 2 = y y + C 2 = y +
= y + C 2
z
2
2
dz
u x = zy C 3 y + C 2 z + C1
u y
= E xy + 3 = z + C 3 u y = zx + C 3 x + 2 ( z )
x
u y
= E yy = 0
y
u y
d 2 ( z )
= E yz 1 = x + x C1 = x +
= x C1
z
2
2
dz
u y = zx + C 3 x C1 z + C 2
u z
= E xz 2 = y + y C 2 = C 2 u z = C 2 x + 3 ( y )
x
2
2
d 3 ( y )
u z
= E yz + 1 = x x + C1 = C1 =
y
2
2
dy
Problemas Resueltos
114
6 Elasticidad lineal
Problemas Resueltos
u z
= E zz = 0
z
u z = C 2 x + C1 y + C 3
z =0
= 0 , x, y C 3 y + C1 = 0 , x, y
C1 = C 3 = 0
uy
z =0
= 0 , x, y C 3 x + C 2 = 0 , x, y
uz
z =0
= 0 , x, y C 2 x + C1 y + C 3 = 0 , x, y
C 2 = 0
C 3 = C1 = C 2 = 0
Finalmente queda:
zy
u(x, y, z ) = zx
Es decir, =
u 2 = 2 z 2 (y 2 + x 2 ) = 2 z 2 r 2 u = zr
115
6 Elasticidad lineal
Problemas Resueltos
0
u(r, , z ) = rz
0
Para hallar (r, , z ) a partir de u(r, , z ) se usarn las ecuaciones que relacionan
desplazamientos y deformaciones en coordenadas cilndricas. El campo de deformaciones
resultante en cilndricas es:
0
(r, , z ) = 0
0
0
r
2
0
r
2
0
= 2G
0 Gr 0
y
GR
z = rG ,
zz = 0
GR
116
6 Elasticidad lineal
Problemas Resueltos
R = tdS
S
M = r tdS
S
( rdS)e = ( GrrdS)e
S
= G r 2 dS e z = GI 0 e z
S
M = GI 0
MT
no es ninguna
GI 0
117
6 Elasticidad lineal
Cuestiones Propuestas
CUESTIONES PROPUESTAS
CP 6-1 Obtener la ecuacin constitutiva de un material elstico lineal a partir del potencial
elstico.
CP 6-2 Justificar si es cierta o no la siguiente afirmacin: En un material termoelstico, la
potencia tensional es una diferencial exacta.
CP 6-3 Sea la expresin del potencial elstico:
u ( ) = 1 Ke 2 + G ij ij
2
E
K =
3
(
1
2)
G=
E
2(1 + )
CP 6-9 Explicar de forma esquemtica (no hace falta escribir las ecuaciones en forma
explcita) el esquema de resolucin (obtencin de , , u ) de un problema esttico,
elstico-lineal, planteado en desplazamientos, indicando cundo se utilizan las diversas
condiciones de contorno.
______________________________________________________________________
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
118
6 Elasticidad lineal
Cuestiones Propuestas
119
6 Elasticidad lineal
Problemas Propuestos
PROBLEMAS PROPUESTOS
z
R
h
HIPTESIS:
1) Propiedades del material: =
= (r ) = 0 + 1 r
PP 6-2 Un cilindro de radio R i est situado en el interior de una corona cilndrica de radio
interior R i + 2e y radio exterior R e . Entre el cilindro y la corona se encuentra una junta
elstica de radio interior R i y cuyo espesor es e . Se somete a la corona a una presin
exterior p . Se pide:
a) Obtener el campo de desplazamientos, deformaciones y tensiones en el cilindro y la
corona.
b) Dibujar los grficos U r - p , donde U r es el desplazamiento radial, y rr - p , donde rr
es la tensin radial para los puntos A , B y C de la figura.
C
Re
junta
A
Ri
DATOS:
1) R i = 1
2) R e = 2
3) = 0
cilindro
corona
e e
120
6 Elasticidad lineal
Problemas Propuestos
4) E (Mdulo de Young)
HIPTESIS:
1) Ley constitutiva de la junta elstica: p * = K * , donde p * es la presin que acta sobre
la junta, * es el acortamiento de su espesor y K es su mdulo de deformacin.
2) e << R i .
3) Se admite un comportamiento de deformacin plana en pequeas deformaciones.
PP 6-3 Un cilindro de material istropo, elstico y lineal, se encuentra apoyado sobre una
base rgida. A una distancia a muy pequea ( a << h ) de su cara superior se encuentra
otra superficie rgida. Sobre el cilindro acta una presin uniforme p en toda la superficie
lateral. Dibujar, acotando los valores significativos, las siguientes curvas:
a) Curva p - , donde es el acortamiento del radio R del cilindro.
b) Curva p - A , donde A es la tensin normal a la superficie de contacto inferior en el
punto A .
a
HIPTESIS:
1) No se considerar la accin del peso propio.
2) Coeficientes de Lam: =
3) Se considera el problema cuasiesttico.
PP 6-4 Dos cilindros macizos constituidos por materiales elsticos distintos se encuentran
superpuestos verticalmente y confinados entre dos paredes infinitamente rgidas. Los
cilindros se hallan sometidos a la accin de las presiones exteriores p y p ( p > 0 , > 0 )
tal como se indica en la figura. Se pide:
a) Obtener el campo de desplazamientos para los dos cilindros en funcin de las
constantes de integracin (a partir de hiptesis razonadas).
121
6 Elasticidad lineal
Problemas Propuestos
b) Indicar todas las condiciones de contorno a aplicar para los distintos contornos del
problema.
c) Suponiendo que el valor de la constante es tal que la seccin de contacto de los dos
cilindros no tiene desplazamiento vertical, calcular las constantes de integracin y el
valor de .
z
HIPTESIS:
1) Cilindro superior: 1 = 1
2) Cilindro inferior: 2 = 2
3) Se supone rozamiento nulo entre
cilindro-cilindro y cilindro-pared.
4) Se desprecia el peso propio.
h/2
h/2
r
R
A
R2
HIPTESIS:
1) Mdulo de elasticidad: E
2) Coeficiente de Poisson: = 0
3) R 1 = R
4) R 2 = 2 R
R1
122
6 Elasticidad lineal
Problemas Propuestos
junta
h
x
L/2
L/2
ALZADO
HIPTESIS:
1) Suponer desplazamientos de la forma: u = [ u (x ) , v(y ) , w (z ) ]T
2) Material elstico lineal.
3) = .
4) = ctte .
5) Peso del ral despreciable.
z
a
SECCIN
123
7 Elasticidad plana
Cuestiones Resueltas
7 Elasticid a d plana
CUESTIONES RESUELTAS
CR 7-1 Justificar si son ciertas o no las siguientes afirmaciones:
a) Si un estado tensional bidimensional tiene un punto singular, todas las isoclinas pasan
por dicho punto.
b) Si un estado tensional bidimensional es uniforme, todas las lneas de mxima tensin
tangencial son paralelas entre s.
Resolucin:
a) Un punto singular se define de la siguiente manera:
El estado tensional se
representa por un punto
1 = 2
=0
Por lo tanto, todas las direcciones son principales y dado un ngulo cualquiera, la
tensin principal formar un ngulo con el eje x . Por dicho punto pasar entonces
una isoclina de ngulo , y como esto se cumple para cualquier , todas las isoclinas
pasarn por dicho punto.
La afirmacin es cierta.
b) Al tener un estado uniforme de tensiones, esto implica que el Crculo de Mohr sea igual
en todos los puntos, y por tanto, los planos de mxima tensin cortante sern los
mismos en todos los puntos, as que las lneas de mxima tensin cortante sern rectas
paralelas entre s.
La afirmacin es cierta.
Direcciones de max
124
7 Elasticidad plana
Cuestiones Resueltas
CR 7-2 La placa rectangular de la figura est sometida a los siguientes estados de tensin
plana:
y
1) x = 0 ; y = b > 0 ; xy = 0
2) x = 0 ; y = 0 ; xy = my , m > 0
Dibujar para cada estado las lneas isostticas, las lneas de mxima y mnima tensin
tangencial e indicar los puntos singulares.
Resolucin:
ESTADO TENSIONAL I:
Direccin de max
2
1
1
II
Direccin de min
Lneas isostticas:
Isostticas I
Isostticas II
max
min
125
7 Elasticidad plana
Cuestiones Resueltas
2
Direccin de max
y>0
II
my
my
Direccin de min
II
y<0
my
Direccin de min
P
my
2
Direccin de max
Lneas isostticas:
Isostticas I
Isostticas II
Puntos singulares
max
min
Puntos singulares
126
7 Elasticidad plana
Cuestiones Resueltas
; xy = ax ; z = b
(a > 0 , b > 0 )
Dibujar en el plano de la figura las lneas isostticas, las lneas de mxima y mnima tensin
tangencial e indicar los puntos singulares.
Resolucin:
= ax = 0 x = 0
Puntos singulares:
Las tensiones normales ya cumplen que son iguales.
Isostticas y lneas de mxima y mnima tensin tangencial:
x>0
Direccin de max
y
2
ax
ax
II
1
Direccin de min
1
x<0
y
2
II
I
ax
x
ax
Direccin de min
P
2
Direccin de max
x<0
Isostticas I
Isostticas II
Puntos singulares
x<0
max
min
Puntos singulares
127
7 Elasticidad plana
Cuestiones Resueltas
; y = 2 x 3 3xy 2
; xy = 3x 2 y ; xz = yz = z = 0
Obtener y dibujar los puntos singulares y las redes de isostticas para cada tensin
principal.
Resolucin:
Los puntos singulares se definen segn:
x = y
=0
Se impondr la definicin para obtener la condicin que han de verificar los puntos
singulares:
x = x3 = 0
x=0
3
2
y = 2 x 3xy = 0
= 3x 2 y = 0
x = x3
y = 0
3
2
3
y = 2 x 3 xy = 2 x
x=0
dy y
=
dx
x
x 2 y 2 = C1
xy = C 2
Por tanto, las isostticas son dos familias de hiprbolas equilteras ortogonales entre s.
La recta singular x = 0 divide a la placa en dos regiones. Para identificar la familia de las 1
se tendr que tomar un punto en cada regin:
Punto (1,0) : x = 1 ; y = 1 = 2 ; xy = 0
Punto (-1,0) : x = 1 = 1 ; y = 2 ; xy = 0
Por tanto, la red de las isostticas es la siguiente:
x<0
Isostticas I
Isostticas II
Puntos singulares
x>0
Isostticas I
Isostticas II
Puntos singulares
128
7 Elasticidad plana
Cuestiones Resueltas
CR 7-5 Dibujar las lneas isostticas en la seccin transversal del cilindro de la figura.
Suponer un campo de desplazamientos de la forma:
p
u = Ar + B
r
r
=
0
u
u = 0
z
A>0 ;
B>0
Resolucin:
Se va a considerar la siguiente notacin:
f = f (r ) = u r (r )
1
= (u + u) = 0 f
2
r
0 0
0
0
= f'+ f + 2f ' = 2( + ) A 2 B
rr
r
r2
= f'+ f + 2 f = 2( + ) A + 2 B
r
r
r2
f
zz = f'+ = 2A
r
r = rz = z = 0
I =
II = rr
Isostticas I
Isostticas II
129
7 Elasticidad plana
Cuestiones Resueltas
CR 7-6 Se tiene una tubera circular de radio interior r1 = 100 mm y exterior r2 = 250 mm ,
suficientemente larga para suponer un estado de deformacin plana, sometida a una
presin interior de 1 MPa . Para una seccin normal del tubo se pide:
a) Calcular y dibujar las leyes de variacin de las tensiones principales sobre el espesor.
b) Dibujar las isostticas diferenciando las dos familias
Resolucin:
Se utilizarn coordenadas cilndricas. Considerando un campo de desplazamientos, por
simetra, uno factible sera:
u = [u r (r ) , 0 , 0]
r r r
r
r2
= 1 (u + u ) = 0
2
0
0
A + B2
r
0
0
0
= 2( + ) A 2 B
rr
r2
= 2( + ) A + 2 r 2
zz = 2A
r = rz = z = 0
r = r2
r12
1
=
A
p
rr = p
2( + ) r22 r12
2 2
rr = 0
B = 1 r1 r2 p
2 r22 r12
Y sustituyendo estos valores en las expresiones de las tensiones se obtienen los resultados
finales:
rr =
=
zz =
r12
r22 r12
r12
r22 r12
(1
(1 +
2
r22
r2
r22
r2
) p = 0.1904 119.2047
r
) p = 0.1904 + 119.2047
r
2
r1
r
p = 2 2 1 2 p = 0.3808
2
2
+ r2 r1
r2 r1
130
7 Elasticidad plana
Cuestiones Resueltas
I =
II = rr
Isostticas I
Isostticas II
1.3808
0.3808
r
100
250
-1
rr
131
7 Elasticidad plana
Cuestiones Propuestas
CUESTIONES PROPUESTAS
CP 7-1 Definir cules son las hiptesis bsicas que se realizan en los estados de tensin
plana y deformacin plana, y cules son sus consecuencias. Citar, para cada caso, un
ejemplo en el que se puedan aplicar dichas hiptesis.
______________________________________________________________________
CP 7-2 Describir el problema tipo de deformacin plana: Geometra, cargas y coacciones
actuantes, hiptesis, campo de desplazamientos, deformaciones, tensiones y relaciones
tensin-deformacin.
______________________________________________________________________
CP 7-3 Describir el problema tipo de tensin plana: Geometra, cargas y acciones
actuantes, campo de desplazamientos, deformaciones y tensiones. Poner un ejemplo de
aplicacin en Ingeniera Civil.
______________________________________________________________________
CP 7-4 Definir lneas isostticas y lneas de mxima tensin cortante. Deducir su ecuacin
diferencial y justificar el ngulo que forman entre s.
______________________________________________________________________
CP 7-5 Definir: Isoclinas, isobaras, isostticas, lneas de mxima tensin cortante, puntos
singulares y puntos neutros.
______________________________________________________________________
CP 7-6 Definir y esquematizar grficamente los siguientes conceptos en estados planos de
tensin-deformacin: lneas isostticas, lneas isoclinas, lneas isobaras, puntos singulares y
puntos neutros.
133
8 Plasticidad
Cuestiones Resueltas
8 Plasticid a d
CUESTIONES RESUELTAS
3
Resolucin:
a) La afirmacin es cierta ya que el origen representa el estado neutro y no puede estar en
estado de plastificacin, sino que est en el interior de la superficie de fluencia, lugar
geomtrico del campo elstico.
b) La afirmacin es falsa pues los tres ejes han de ser de simetra, y en la figura no lo son.
CR 8-2 Para cada uno de los siguientes casos, justificar la forma de la superficie de fluencia
en el espacio de tensiones principales:
a) f (I12 ) = 0 .
b) f (J 2 ) = 0 .
2
c) aI12 + b oct
= c , a, b, c estrictamente positivos .
Resolucin:
a) En este primer caso se tiene una condicin sobre la tensin media ya que se cumple:
I 1 = x + y + z = 3 m
Como el primer invariante est elevado al cuadrado, esto impone la distancia a planos
octadricos, tanto en un sentido como en el otro.
3
3 m
3 m
Eje de tensin
hidrosttica
134
8 Plasticidad
Cuestiones Resueltas
distancia = 3 oct
3
3 oct
2
Eje de tensin
hidrosttica
d
De donde se deducen las siguientes relaciones:
d = x = 3 oct = 3 I 1 I 1 = 3x
3
oct = R
R = y = 3 oct
3
2
aI12 + b oct
y
by
2
x
+
=1
= c 3ax +
=c
3
c
3c
3a
b
2
Esta es la expresin de una elipse en ese plano x - y . Adems, como no interviene el 3er
invariante en la definicin de la superficie de fluencia, habr simetra de revolucin
alrededor del eje hidrosttico y, por lo tanto, la rotacin de la elipse alrededor del eje x
( eje hidrosttico) acotar la superficie definitiva.
En conclusin, si se toman unos ejes x ( hidrosttico), y , z la superficie de fluencia
tendr la siguiente expresin:
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
135
8 Plasticidad
Cuestiones Resueltas
x +
+ z =1
c
3c
3c
3a
b
b
3
z y
x
Eje de tensin
hidrosttica
CR 8-3 Justificar cules de los criterios de fluencia plstica, Von-Mises, Tresca, DruckerPrager y Mhr-Coulomb, podran ser adecuados para modelizar un material en cada uno de
los siguientes casos.
a)
b)
c)
d)
Resolucin:
a) Que la plastificacin dependa de J 2 , es decir, que dependa de oct , significa que
depende nicamente de la distancia del estado tensional al eje hidrosttico. Por lo tanto,
la superficie de fluencia ser un cilindro de directriz el eje de tensiones hidrostticas. El
nico criterio con esta forma es el de Von-Mises.
b) Aqu no hay dependencia de la variable J 3 , esto implica que existe simetra respecto al
eje de tensiones hidrosttico, cosa que cumplen los criterios de Von-Mises y DruckerPrager. Pero como la superficie de fluencia depende de I1 , es decir, de la distancia al
estado neutro (1 = 2 = 3 = 0) , Von-Mises se descarta.. Queda, finalmente, el criterio
de Drucker-Prager.
c) Al depender de I1 , no cumplirn esta condicin ni Von-Mises ni Tresca. Al depender
de J 3 , tampoco lo har Drucker-Prager. Por lo tanto, nicamente queda el criterio de
Mhr-Coulomb.
d) Al no depender de I1 , aceptaremos Von-Mises y Tresca. Como depende de J 3
descartaremos Von-Mises, y por lo tanto, quedar el criterio de Tresca.
136
8 Plasticidad
Cuestiones Resueltas
c
R=
A
2
tg = c
B
1
A
2
sen =
B 1 A
= arcsen 2
B 1 A
c = B tg
137
8 Plasticidad
Cuestiones Resueltas
C
e
2
d = c e
tg
e + d
e
c = 1 1 1
e
e
tg
2 sen
e
1
sen = 2
d=
e
e
2 sen
d
+
2
c = e 1 + 1
tg 2 sen
CR 8-6 Una probeta de hormign sometida a compresin uniaxial uniforme rompe segn
un plano inclinado que forma un ngulo con el de carga. Determinar dicho ngulo en
funcin del ngulo de rozamiento interno y de la cohesin del material.
Resolucin:
En el plano de tensiones de Mohr se representar la rotura de la siguiente manera:
138
8 Plasticidad
Cuestiones Resueltas
c
P
=
4 2
x =
xy =
y = 0 ; xz = 0
= 0
yz
z = 0
M
y
2
2
e
2
R
A
1
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
139
8 Plasticidad
Cuestiones Resueltas
2
2
R = + 2 = + 4
2
2
e
2
la Figura 2.
e
G
Figura 1
Figura 2
Resolucin:
El Tensor de Tensiones para el ensayo de la figura 2 tendr la expresin:
0 0
= 0 0 m = 1 ( x + y + z ) = 0 =
3
0 0 0
donde:
J 2 = 1 ( : ) = 1 ( : ) = 2
2
2
e = 0
e
3
; e =
e
3G
140
8 Plasticidad
Cuestiones Resueltas
3 m = 3
Eje de tensin
hidrosttica
2
1
m =
+ b = 0 a + b = 0
0 0
0 0 = u
3
0 0
0
2 0
0 1 0
0 0 1
= 2 u
3
a=
2
u
3
u
3
b=
2
u
3
u
3
2 = a 1 + b
u
3 u
3
141
8 Plasticidad
Problemas Resueltos
PROBLEMAS RESUELTOS
PR 8-1 La pieza metlica PQRS tiene un espesor e y est constituida por materiales
distintos, (1) y (2), considerados elastoplsticos perfectos. Se somete la pieza a un ensayo de
corte puro en la mquina de la Figura 1 de tal forma que se produce un estado de tensiones
y deformaciones uniforme dado por:
x = y = z = 0 ; xz = yz = 0 ; xy = =
h
x = y = z = 0 ; xz = yz = 0 ; xy = 0
Cuando se ensaya por separado una pieza de cada uno de los materiales se obtiene una
curva - como la de la Figura 2. Se pide:
a) El lmite elstico que se obtendra en un ensayo de tensin uniaxial de cada material por
separado, suponiendo que obedecen al criterio de Von-Mises.
Cuando se ensaya la pieza constituida por los dos materiales se obtiene la curva P de la
Figura 3. Se pide:
b) El valor de la carga y del desplazamiento elstico Pe y e .
c) dem para los valores plsticos Pp y p .
d) Las coordenadas P de los puntos C y D de la Figura 3.
h
h/2 h/2
Marco rgido
(x)
R
(1)
(2)
(1)
(2)
Figura 1
Pp
Pe
Figura 2
3e
A
e
C
D
Figura 3
142
8 Plasticidad
Problemas Resueltos
HIPTESIS:
Material (1): G = G ; e = *
Material (2): G = G ; e = 2 *
Resolucin:
a) Con los datos anteriores se quiere obtener el lmite elstico que se obtendra en un
ensayo de tensin uniaxial de cada uno de los materiales por separado, suponiendo que
obedecen el Criterio de Von-Mises. As se sabe que la plastificacin empieza cuando:
= e
= (3J 2 ) 2
J 2 = 1 Tr ( 2 )
2
= esf
esf = m 1
m = 1 Tr ()
3
2
0 0
= 0 0 = 2 = 0
0
0 0 0
0
2
0
0 J 2 = 2 = 3
0
b) Al ensayar la pieza constituida por los dos materiales se obtiene un grfico como el de la
Figura 3. Se quiere hallar el valor de la carga Pe y desplazamiento elstico e , que
marcan el final de la fase de comportamiento elstico de la pieza. En el enunciado se
dice que al ensayar por separado cada material, el grfico resultante es el de la Figura 2,
siendo e = para (1) y e = 2 para (2). Adems, se dice que G es comn a los dos
materiales, es decir, tienen la misma pendiente en la grfica .
Ahora bien, se busca saber el comportamiento conjunto de la pieza. Es de suponer que
el comportamiento de sta ser elstico mientras los dos materiales estn en su
dominio elstico correspondiente. Por lo tanto, como el intervalo elstico del material
(1) es menor, ser este material el que defina el dominio elstico del conjunto,(hasta el
punto A de la Figura 3).
143
8 Plasticidad
Problemas Resueltos
Punto A
Punto A
1A
2A
MATERIAL (1)
MATERIAL (2)
Si se aplica equilibrio entre la fuerza P y las tensiones que tiene cada material en el
punto A resulta:
(hay que igualar fuerzas, por tanto se multiplican las tensiones por la superficie en la que
actan, considerando la unidad como magnitud perpendicular al papel)
Pe = h 1A + h 2A = h + h = h
2
2
2
2
e = 1A h = 2A h = h
G
Resumiendo:
Pe = h
e = h
G
Punto A
Punto A
Punto B
Punto B
G
1A
1B
MATERIAL (1)
G
2A
2B
MATERIAL (2)
144
8 Plasticidad
Problemas Resueltos
2
2
2
2
p = 1B h = 2B h = 2 h
G
Pp = 3 h
2
p = 2 h = 2 e
G
Punto B
Punto B
3 A
3 A
Punto D
Punto C
Punto C
Punto D
MATERIAL (1)
MATERIAL (2)
'B = 3 A = 3 h = 3 e
G
PB' = PB = 3 h
2
145
8 Plasticidad
Problemas Resueltos
PC = h 1C + h C2 = h ( ) + h 0 = h
2
2
2
2
2
C = 1C h = h
PC = h
2
C =
= e
G
Y para el punto D :
PD = h 1D + h 2D = h ( ) + h (2 ) = 3 h
2
2
2
2
2
D = 2D h = h
PD = 3 h
2
h
D =
= e
G
T
u z ] = 0
rz 0
h
146
8 Plasticidad
Problemas Resueltos
M
Material (2)
Material (1)
e
R/2
Figura 2
h
M
Mp
M2
M1
Me
Figura 1
1 2
Figura 3
HIPTESIS:
Material (1): G = G ; e = *
Material (2): G = G ; e = 2 *
Resolucin:
a) Se buscar primero el tensor de deformaciones en cilndricas (pequeas deformaciones)
a partir del campo de desplazamientos. As se obtiene el siguiente resultado:
0
= 0
0
0
r
2h
0
r
2h
0
147
8 Plasticidad
Problemas Resueltos
= 0
rr
0
Gr
h
0
0
0
Gr
h
GR
h
R
h
r
0 e
GR
h
Gr
0 e el momento actuante vale:
h
4
R
R Gr 3
Gr
r
rdrd = 2
dr = GR
(1)
0 h
0
h
2h
r(r)dS =
S
e = h
GR
148
8 Plasticidad
Problemas Resueltos
Este es el valor del ngulo de giro del extremo libre del cilindro para el que empieza la
plastificacin en los puntos exteriores de la pieza (material (1)). El momento
correspondiente se hallar sustituyendo e en la relacin ( 1 ):
4
M e = M ( e ) = GR e
2h
3
M e = R
2
el material (1) tendr una distribucin de tensiones plstica, mientras que en el (2) ser
elstica.
Para calcular este giro 1 se debe imponer
la siguiente condicin:
r = R ; = 1 =
e =
G1 R
=
2h
= 1
1 = 2 h
GR
Este es el valor mnimo de giro para el que el material (1) se encuentra plastificado en su
totalidad.
La relacin ( 1 ) entre M y aqu no es vlida ya que el material (1) se comporta
elastoplsticamente y el (2) lo hace todava elsticamente. La relacin en este apartado es
la siguiente:
M1 =
r rdrd +
R
2
R
2
0
G1 r
r
rdrd = 2
h
r 2 dr + 2G
R
2
1
h
R
2
0
r 3 dr
M 1 = 31 R 3
48
149
8 Plasticidad
Problemas Resueltos
= 2
G 2 R
= 2
2h
2 = 4 h
GR
r rdrd +
R
R
2
0
r G 2
h
r rdrd
M 2 = 17 R 3
24
= p
2
r rdrd +
R
2
R
2
0
r 2 rdrd
M p = 3 R 3
4
151
8 Plasticidad
Cuestiones Propuestas
CUESTIONES PROPUESTAS
CP 8-1 Explicar los conceptos de plano octadrico, oct y oct . Obtener las expresiones de
oct y oct en funcin de las tensiones principales.
______________________________________________________________________
CP 8-2 Justificar por qu las superficies de fluencia para materiales istropos son simtricas
respecto a tres planos en el espacio de tensiones principales.
______________________________________________________________________
CP 8-3 Justificar la siguiente afirmacin: Para una superficie de fluencia plstica de la
forma f (J 2 ) = 0 , un estado de tensiones hidrosttico nunca produce plastificacin.. A qu
superficie de fluencia conocida sera aplicable dicha afirmacin?
______________________________________________________________________
CP 8-4 Demostrar que si la superficie de fluencia plstica no depende del primer invariante
del tensor de tensiones, la deformacin plstica es incompresible.
______________________________________________________________________
CP 8-5 Indicar de forma concisa cules son las hiptesis fundamentales de la teora
matemtica de la plasticidad.
______________________________________________________________________
CP 8-6 Enunciar las hiptesis fundamentales de la teora matemtica incremental de la
plasticidad en tres dimensiones.
______________________________________________________________________
CP 8-7 Interpretacin de la condicin de consistencia en la teora de la plasticidad.
CP 8-8 Criterio de plastificacin de Tresca: definicin, expresin en funcin de las
tensiones principales, justificacin de su forma en el espacio de tensiones y particularizar al
caso de una dimensin.
CP 8-9 Justificar los motivos por los que el criterio de Tresca ( F 1 3 K ), recibe el
nombre de criterio de mxima tensin cortante.
______________________________________________________________________
CP 8-10 Demostrar que la superficie de fluencia de Von-Mises puede expresarse en
funcin de las tensiones principales como:
1
3 2 + 2 + 2 1 ( + + )2 2 = 0
u
1
2
3
1
2
3
2
3
_____________________________________________________________________
CP 8-11 Justificar si llegar a plastificar segn los criterios de Von-Mises, Mhr-Coulomb y
Drucker-Prager, un material sometido a un estado de:
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
152
8 Plasticidad
Cuestiones Propuestas
m a
3
______________________________________________________________________
CP 8-15 Indicar la forma que tendran las siguientes superficies en el espacio tridimensional
de tensiones principales:
m
a) m = ctte
c)
= ctte
oct
b) oct = ctte
______________________________________________________________________
______________________________________________________________________
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
153
8 Plasticidad
Cuestiones Propuestas
(1 3 )2 m 2 < 0
oct
K2 < 0
m
(K, m cttes.)
b
a
1
2
Interseccin con plano
octadrico por (0,0,0)
______________________________________________________________________
CP 8-21 Un material plstico perfecto se somete a un ensayo de traccin uniaxial
( 1 = , 2 = 3 = 0 ) y empieza a plastificar para un valor de la tensin = e . Si su
ngulo de rozamiento interno es de 30 , determinar el valor de la cohesin que habra que
asignar al material para poder aplicarle el criterio de Mhr-Coulomb.
______________________________________________________________________
CP 8-22 Un material elastoplstico se somete a un ensayo de corte puro (I) y a un ensayo
de traccin uniaxial (II). La plastificacin se produce, respectivamente, para = a y = b .
Determinar los valores de la cohesin y del ngulo de rozamiento interno suponiendo un
modelo de Mhr-Coulomb.
(I)
(II)
154
8 Plasticidad
Cuestiones Propuestas
______________________________________________________________________
155
8 Plasticidad
Problemas Propuestos
PROBLEMAS PROPUESTOS
PP 8-1 Se considera una esfera maciza de radio R 1 rodeada de una corona esfrica de radio
interior R 1 y radio exterior R 2 . La esfera y la corona son del mismo material y estn
inicialmente en contacto sin ejercerse ninguna presin. En un cierto instante se empieza a
calentar la esfera interior hasta un incremento de temperatura . Se pide:
a) Determinar el valor de la presin exterior que sera
necesaria aplicar a la corona para que sta no sufra
cambio de volumen (hiptesis de pequeas
deformaciones).
b) Determinar los campos de desplazamientos,
de
deformaciones y de tensiones, en la esfera y la corona
en estas condiciones.
c) Obtener el menor valor de para el cual se inicia la
plastificacin en algn punto, suponiendo las anteriores
condiciones y considerando el criterio de Von-Mises.
R2
R1
HIPTESIS:
1) Propiedades del material: E , = 0 , , f , R 1 = 1 , R 2 = 3 .
2) Se supondrn despreciables las fuerzas msicas.
3) Se conoce la solucin en desplazamientos y tensiones para una corona esfrica de radio
interior R i y exterior R e sometida a una presin interior Pi y exterior Pe es, para
=0:
u r (r )
u= 0
0
rr
= 0
0
u r = Cr +
0
0
C1
r2
C=
3
3
Pi Ri Pe Re
3
3
E Re Ri
C1 =
3
3
Pi Pe Ri Re
2 E Re 3 Ri 3
2C
rr = E C 31
r
C
= = E C + 31
PP 8-2
a) Obtener la expresin de los criterios de fluencia de Von-Mises y Tresca en funcin de
las tensiones principales y el lmite elstico e para estados tensionales de la forma:
1 > 0 , 2 > 0 , 3 = 0
156
8 Plasticidad
Problemas Propuestos
corona est sometida a una presin interior de valor P . Obtener la curva P , donde
es el desplazamiento radial del punto B de la corona.
d) Calcular los valores de e que garantizan que la plastificacin de la corona empieza
despus del aplastamiento total de la junta. Considerar el criterio de Tresca.
Junta elstica
Re
R
P
Ri
HIPTESIS:
1) Corona: E , = 0 , e .
2) Junta: P = K .
3) Fuerzas de volumen despreciables.
e e
PP 8-3 Se considera una esfera maciza de radio R 1 y material (1) rodeada de una corona
esfrica de radio exterior R 2 y material (2). Ambas esferas estn inicialmente en contacto
sin ejercerse ninguna presin. Se somete al conjunto a una presin exterior P y a un
incremento de temperatura .Se pide:
a) Determinar los valores posibles de y P (positivos o negativos) para que no se
pierda el contacto entre la esfera y la corona. Representarlos grficamente en un
diagrama P - .
b) Obtener el estado tensional de la corona y de la esfera para dichos valores.
c) En estas condiciones, obtener para cada valor de la presin P el valor de * para el
que se inicia la plastificacin en algn punto del conjunto segn los criterios de VonMises y Mhr-Coulomb. Dibujar las correspondientes grficas P - * (grficos de
interaccin).
HIPTESIS:
1) Material (1): E , = 0 , 2 , f , C, = 30 , R 1 = 1
2) Material (2): E , = 0 , , f , C , = 30 , R 2 = 2
3) Se conoce la solucin en desplazamientos y tensiones para una corona esfrica de radio
interior R i y exterior R e sometida a una presin interior Pi y exterior Pe es, para
=0:
u r (r )
3
3
3
3
u= 0
0
rr
= 0
0
u r = Cr +
0
0
C1
r
C=
Pi Ri Pe Re
E R e Ri
C1 =
Pi Pe Ri Re
2 E Re 3 Re 3
2C
rr = E C 31
r
C
= = E C + 31
157
8 Plasticidad
Problemas Propuestos
HIPTESIS:
1) Propiedades del material:
- cohesin: C
- ngulo de rozamiento interno: = 30 .
- constante trmica:
- Constantes de Lam: =
2) No se considerar el peso propio.
3) Rozamiento cilindro-suelo nulo.
A=
(1 ) p q
E
rr = = p
zz = q
,
,
B=
2 p + q
E
r = rz = z = 0
a) Demostrar que la solucin propuesta cumple todos los requisitos necesarios para ser
solucin del problema elstico lineal de la Figura 1.
b) Se consideran ahora los dos cilindros de la Figura 2, separados por una distancia muy
pequea e << H . Se somete el cilindro inferior a un incremento de temperatura
creciente 0 . Se considera el material de los dos cilindros elastoplstico con un
diagrama tensin-deformacin uniaxial como el de la Figura 3. Dibujar, acotando los
valores significativos, las curvas y zz en el punto C del cilindro inferior:
1. Antes del contacto entre los dos cilindros.
2. Despus del contacto y hasta el inicio de la plastificacin.
158
8 Plasticidad
Problemas Propuestos
z
q
e
C
Figura 1
Figura 2
10e
Figura 3
HIPTESIS:
1) Se despreciarn el rozamiento, las fuerzas msicas y las de inercia.
dx
v
=
0
w = (x )
Se pide:
a) Justificar que cualquier seccin normal al eje x permanece, despus de la deformacin,
plana y normal a la deformada del eje x (Hiptesis de Bernouilli).
b) Justificar que, si se desprecia el coeficiente de Poisson, = 0 , el campo de
desplazamientos dado es solucin del problema elstico. Obtener el valor de (x ) y los
correspondientes campos de desplazamientos, deformaciones y tensiones.
c) Obtener el mximo valor del momento M ( M e ) para el cual todos los puntos de la
viga estn en rgimen elstico y el correspondiente valor de la flecha (desplazamiento
w ) en el extremo libre de la viga.
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
159
8 Plasticidad
Problemas Propuestos
y, v
M
x, u
m
HIPTESIS:
1) No se considera el peso propio.
2) Se considerarn pequeas deformaciones.
3) Momento de inercia de la seccin respecto al eje y : I = 1 mh 3
12
Pieza infinitamente
rgida
A
L
Figura 1
160
8 Plasticidad
Problemas Propuestos
Pp
Pe
E
(1)
2 eL
E
E
e
(2)
Pq
Figura 2
HIPTESIS:
1) Mdulo de Young: E
2) Coeficiente de Poisson: = 0
3) Lmite elstico: e
4) Espesor de la pieza: b
e
L . Posteriormente se descarga nuevamente hasta P = 0 . Se pide:
E
C
e
45
45
O
10 e
A
161
Cuestiones Propuestas
9 Ecuacio n es constitutivas
en fluidos
CUESTIONES PROPUESTAS
CP 9-1 Definir los siguientes conceptos:
a) Presin hidrosttica.
b) Presin media.
c) Presin termodinmica.
Justificar en qu casos sern iguales la presin media y la presin termodinmica y cundo
lo sern las tres.
______________________________________________________________________
CP 9-2 Definir los siguientes conceptos:
a)
b)
c)
d)
e)
Fluido de Stokes.
Fluido Newtoniano.
Fluido perfecto.
Fluido incompresible.
Flujo barotrpico.
______________________________________________________________________
CP 9-3 A partir de la frmula general de la ecuacin constitutiva para fluidos viscosos,
deducir la ecuacin constitutiva para un fluido Newtoniano istropo en componentes
esfricas y desviadoras.
NOTA: Utilizar la expresin: ij kl + ( ik jl + il jk ) .
______________________________________________________________________
CP 9-4 Obtener la ecuacin constitutiva de un fluido Newtoniano, barotrpico e istropo
en componentes desviadoras y esfricas.
______________________________________________________________________
CP 9-5 Definir presin hidrosttica, presin media y presin termodinmica para un fluido
Newtoniano y obtener las relaciones entre ellas.
______________________________________________________________________
162
Cuestiones Propuestas
163
10 Mecnica de fluidos
Cuestiones Resueltas
10 Mecnic a de fluidos
CUESTIONES RESUELTAS
CR 10-1 Un fluido perfecto y barotrpico tiene una ecuacin cintica de estado definida
p
, donde k es constante. Obtener la distribucin de presiones en rgimen
k
cuasiesttico (aceleracin nula) bajo la accin del campo gravitatorio b = [0 0 - g ]T .
por = 0 +
Resolucin:
La Ecuacin Constitutiva de un fluido perfecto es:
= p1
El hecho de que sea fluido barotrpico implica que la densidad es una funcin de la
presin. Esta relacin es, precisamente, la Ecuacin Cintica del enunciado:
p
k
= (p) = 0 +
p = p(z)
p
= 0
y
dp
p
dp g
+ p = 0 g
0 + g = 0
dz
k
dz k
Solucin homognea :
z
dp g
+ p = 0 p = Ce k
dz k
Solucin particular :
p = k 0
g
z
Ce k
164
10 Mecnica de fluidos
Cuestiones Resueltas
CR 10-2 El fluido barotrpico del interior de la tubera de la figura tiene por ecuacin
cintica de estado:
p = Ln
0
( y 0 cttes.)
P2
v2
S2
Resolucin:
Segn la Ecuacin de Continuidad (global, espacial) se tiene:
d
dt
dV = 0
V
dV = t dV +
V
v ndS
dV +
v ndS = 0
dV = 0
V
v ndS = 0
= 0e
Ahora tan slo hay que sustituir esta expresin en la anterior para encontrar el resultado
que se pide:
p1
p2
0 e v1S1 = 0 e v 2 S 2
p 2 p1
v1S1
v 2S2
vS
p 2 = p 1 + ln 1 1
v 2S2
165
10 Mecnica de fluidos
Cuestiones Resueltas
CR 10-3 Calcular el valor por unidad de espesor de la fuerza horizontal F que hay que
aplicar en el punto B de la compuerta semicircular de la figura, que puede girar alrededor
de la rtula A , para que la recta AB sea vertical. La compuerta separa dos niveles distintos
de cotas h y h , de un mismo fluido.
F
R
h
A
Hiptesis:
1) Se desprecia el peso de la compuerta y la presin atmosfrica.
Resolucin:
Las nicas fuerzas que actan sobre la compuerta son las fuerzas de presin de los fluidos,
la fuerza F y la reaccin en A (componente horizontal H y componente vertical V ).
F
B
1 gh 2
2
1 g 2 h 2
2
H
Figura 1
Figura 2
Como los fluidos ejercen una presin perpendicular a la superficie de la compuerta y sta es
circular, la resultante de integrar las presiones sobre la superficie pasar por el centro de la
circunferencia que define la compuerta. As, si se plantea el equilibrio de momentos
respecto del centro de la circunferencia (ver Figura 1) resulta:
FR = HR H = F
Y si ahora se impone equilibrio de fuerzas horizontales, sabiendo que los fluidos ejercen
una presin horizontal con una distribucin triangular (ver Figura 2), se tiene:
2F + 1 (gh )(h ) = 1 (gh)h
2
2
F = 1 gh 2 (1 2 )
4
166
10 Mecnica de fluidos
Cuestiones Resueltas
45
m/2
45
Resolucin:
Si se consideran unos ejes cartesianos en el vrtice del cajn, se puede calcular su centro de
masas:
Consideramos que todo el cajn tiene densidad 1 resulta:
Peso :
W1 = 1 m(2m)1 = 1 m 2
2
Centro de gravedad :
y G1 = 2 m
3
Como se ha considerado que la punta del cajn tena densidad 1 , ahora se deber
considerar que tiene densidad 31 1 = 21 :
Peso : W2 = 1 m m(21 ) = 1 1 m 2
2 2
2
Centro de gravedad : y G 2 = 2 m = m
3 2 3
W1 y G1 + W2 y G 2
W1 + W2
=5m
9
El empuje del fluido ser equivalente al volumen de fluido desplazado y, por estar en
equilibrio, deber ser igual al peso del cajn:
WE = 1 h(2h ) 0 = h 2 0
2
3 1
h = m
2 0
W1 + W2 = 3 m 2 1
El centro de gravedad donde se aplica esta fuerza es el centro geomtrico del volumen de
fluido desplazado, es decir:
yg = 2 h = m 2 1
3
3 0
Para que el equilibrio sea estable es necesario que el centro geomtrico del fluido sea
superior al del cajn:
yg > yG
1 > 25 0
54
167
10 Mecnica de fluidos
Cuestiones Resueltas
CR 10-5 Admitiendo que el contenido de una copa de cava es un fluido perfecto y que la
presin termodinmica es igual a la hidrosttica, calcular la aceleracin con que suben las
burbujas de gas de la copa desde el fondo de la misma.
a
Hiptesis:
1) Densidad del gas: a
2) Densidad del fluido: c
3) c = a
Resolucin:
Si se admite que el contenido es un fluido perfecto, entonces se cumple que:
= p1
= p1
Segn el Teorema de Arqumedes se puede evaluar el empuje ascendente que ejerce el cava
sobre la burbuja como:
E = c Vg
Por lo tanto, la diferencia entre estas dos fuerzas ser la resultante que acte sobre la
burbuja:
F = ma = a Va
c a
g
a Va = ( c a )Vg a =
F=EW
a
a = ( 1)g
CR 10-6 Un depsito circular de gran dimetro que est lleno de agua vierte por un
pequeo orificio lateral situado a una altura H por debajo del nivel del agua en el depsito.
Si el caudal vertido es Q , obtener el dimetro D del orificio.
H
B
168
10 Mecnica de fluidos
Cuestiones Resueltas
Resolucin:
Como el depsito es de gran dimetro, se puede suponer que la diferencia de cotas H se
mantiene constante con el tiempo, por lo que se tendr rgimen estacionario.
Por otro lado, como se est tratando con agua, se est ante un fluido perfecto e
incompresible.
As, se concluye que se puede aplicar la Ecuacin de Bernouilli:
z+
p v2
+
= ctte.
g 2g
+0
g
v B = 2gH
p atm v 2B
B: 0 +
+
2g
g
A: H +
Finalmente se tiene:
Q = v B S B = 2gH D
4
D=
4Q
2gH
v2=2m/s
h2
Hiptesis:
1) Se desprecia la presin atmosfrica.
Resolucin:
Teniendo en cuenta la Ecuacin de Continuidad se puede establecer la siguiente relacin en el
volumen de control que definen las secciones consideradas:
v1 h 1 = v 2 h 2
h 2 = h1
v1
v2
h 2 = 1 h1
2
169
10 Mecnica de fluidos
Cuestiones Resueltas
1 : (H + h 1 ) + 0 +
v2
2: h2 + 0 + 2
2g
v12
2g
v 22 v 12
H = h 2 h1 +
2g
H = h 1 + 3
2g
&
S1
S2
Hiptesis:
1) Suponer que el fluido es perfecto, incompresible y que el rgimen es estacionario.
2) Se desprecia la presin atmosfrica.
3) S1 y S 2 son secciones transversales.
4) La densidad del fluido es .
Resolucin:
Como es un fluido perfecto el estado tensional ser de la forma:
= p1
v2 =
S1
S
v1 = 1 &
S2
S2
& 1
2g
S1
p = 1 & 2 gH
2 S 2
De aqu se deduce que p ha de ser constante para cualquier punto en contacto con el
mbolo ( x = H ). Entonces se obtiene que:
p = ctte. , x = H F = pS1
170
10 Mecnica de fluidos
Cuestiones Resueltas
F=
S1
S1 1& 2 gHS1
2 S 2
h/2
v(y)
h/2
Hiptesis:
1) Se desprecian las fuerzas msicas.
2) Flujo estacionario.
Resolucin:
La Ecuacin de Navier-Stokes tiene la siguiente expresin:
p + ( + ) ( v) + 2 v + b = dv
dt
v = 0
t
0
v(y)
v v = [v(y) 0 0]
y
0
0 0
0 0 = 0
0 0
Finalmente, queda:
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
dv = 0
dt
171
10 Mecnica de fluidos
Cuestiones Resueltas
p + 2 v = 0 a +
v(y) a
= y + C1
y
2 v(y)
=0
y 2
2 v(y) a
=
y 2
v(y) = a y 2 + C1 y + C 2
2
a h + C h + C = 0
1
2
2 2
2
2
a h C h + C = 0
1
2
2 2
2
2
C1 = 0
2
C = ah
2
8
Finalmente se obtiene:
2
v(y) = a y 2 h
2
4
v(y)
CR 10-10 Un depsito circular de gran dimetro y seccin S A est lleno de agua. Este
depsito vierte por un pequeo orificio lateral de seccin S B situado a una distancia H(t)
por debajo del nivel del agua en el depsito. Obtener la ecuacin diferencial
(f(H,H& ,S A , S B ) = 0) que gobierna el problema de vaciado (se supone rgimen estacionario).
SA
H(t)
B
SB
Resolucin:
Segn el Teorema de Bernouilli se verifica esta expresin para una lnea de corriente:
p v2
+
= ctte.
g 2g
Por lo tanto, aplicando esta propiedad a los puntos A y B se obtiene:
z+
p
v
p
v
zA + A + A = zB + B + B
2g
g
g
2g
172
10 Mecnica de fluidos
Cuestiones Resueltas
H(t) +
p atm v 2A
p
v2
+
= 0 + atm + B
g
2g
g
2g
2gH(t) + v 2A v 2B = 0
Basta plantear, por ltimo, las siguientes consideraciones sobre velocidades del fluido en las
secciones consideradas:
Q
&
= H
SA
&
S
H
Q
A
=
vB =
SB
S B
vA =
2
& 2 1 S A
2gH + H
S2
= 0
R(t)
Resolucin:
Como se tiene simetra esfrica, se va a considerar la velocidad en coordenadas esfricas. Se
har la hiptesis de que este campo de velocidades es de la forma:
v = (v r (r) , 0 , 0) T
1 (r 2 v ) = 0 r 2 v = A v = A
r
r
r
r 2 r
r2
r =R
dR(t)
= R (t)
dt
A = R A = R 2R
R2
Finalmente queda:
vr =
R 2 (t)R (t)
r2
173
10 Mecnica de fluidos
Cuestiones Resueltas
CR 10-12 A partir del teorema de las fuerzas vivas demostrar que para un fluido
Newtoniano incompresible que se mueve en el interior de un recipiente rgido, inmvil y
cerrado y despreciando las fuerzas msicas, se cumple:
dv 2 = 4 d : d
dt
Hiptesis:
1) Condicin de adherencia: v = 0 , en
Resolucin:
El Teorema de las Fuerzas Vivas tiene la siguiente expresin:
Pe =
b vdV + t vdS = dt
V
1 v 2 dV +
2
: d dV
V
Como no se consideran fuerzas msicas ( b = 0 ) y como el recipiente est inmvil, cosa que
implica por adherencia que la velocidad del fluido en el contorno es nula, entonces se
puede concluir:
Pe = 0
d
dt
1 v 2 dV = : d dV
V
2
A partir del Teorema del Transporte de Reynolds se puede reescribir la expresin anterior como:
1 dv 2 dV =
2 dt
V2
: d dV , V
1 dv 2 = : d
2 dt
Finalmente, se concluye:
dv 2 = 4 d : d
dt
175
10 Mecnica de fluidos
Problemas Resueltos
PROBLEMAS RESUELTOS
PR 10-1 Un flotador en forma de pirmide y de base cuadrada es utilizado para obturar el
orificio (cuadrado) de desage de un depsito que contiene agua. Se pide:
a) Obtener la resultante del empuje hidrosttico F ejercido sobre el flotador para todos
los posibles niveles de agua H del depsito, es decir:
0<Hb
Hb
a
Alzado
HIPTESIS:
1) Se desprecia la presin atmosfrica.
Planta
Resolucin:
a) El empuje hidrosttico se puede calcular como el empuje que sufrira la pirmide sin
tener en cuenta la presencia del orificio y, posteriormente, restarle la fuerza
correspondiente al empuje del agua del agujero, ya que en el problema sta no existe.
Caso 0 < H b
176
10 Mecnica de fluidos
Problemas Resueltos
Pero como el agujero est vaco, esta fuerza deber ser restada a la total ya que en
la realidad no acta:
Fresultante = Ftotal Fagujero
Fresultante = g 1 b 3 (b H )3 gHa 2
3
Caso H b
F max = F(b a )
F(b ) = gb 1 b 2 a 2
3
ga 2
b-a
H0
Donde el valor de H para F max se obtiene derivando F(H) , del primer caso, respecto
de H e igualando a cero:
F max = F(b a ) = g 1 b 3 + 2 a 3 ba 2
3
Y el valor del calado H 0 para el que la fuerza ejercida por el agua sobre el flotador vale
cero:
177
10 Mecnica de fluidos
Problemas Resueltos
3
H0 = b 2
3a
H1
H2
H
H0
PR 10-2 Sobre una placa rgida indefinida, de densidad * y espesor t , acta una fuerza
tangencial f * por unidad de superficie. La placa desliza a una velocidad v * sobre un plano
inclinado con un ngulo . Entre la placa y el plano inclinado hay dos fluidos
Newtonianos, distintos e inmiscibles, de viscosidades 1 y 2 , que se distribuyen en dos
capas del mismo espesor h . Se pide:
a) Establecer y razonar las hiptesis correspondientes a los campos de presin y de
velocidad.
b) Integrar las correspondientes ecuaciones diferenciales y obtener, salvo las constantes de
integracin, la distribucin de velocidades y presiones en cada fluido.
c) Indicar y justificar las condiciones de contorno que hay que aplicar para determinar
dichas constantes de integracin.
d) Determinar completamente los campos de velocidades y de presin, as como las
tensiones en cada fluido. Dibujar la distribucin de cada variable (velocidades, presin y
tensiones) sobre una seccin transversal como la A - A' , acotando los valores
significativos.
e) Obtener el valor de v * en funcin de f * y el caudal q que pasa por una seccin
semicircular como la B - B' .
178
10 Mecnica de fluidos
Problemas Resueltos
y
t
f*
v*
B'
A'
h
x
1
h
2
B
HIPTESIS:
1) Fluidos incompresibles.
2) Rgimen estacionario.
3) Fuerzas de masa despreciables en los fluidos.
4) Presin atmosfrica despreciable.
Resolucin:
a) Comentar primero que la dimensin z perpendicular al papel no interviene.
Se supondrn una presin y una velocidad del tipo:
p = p(x, y)
v = v (x, y) = [v x (x, y) 0 0]T
v = [v x (y) 0 0]T
0=
2vx
p
+
x
y 2
179
10 Mecnica de fluidos
Problemas Resueltos
componente y
0=
p
y
p = p(x)
f(x) = f(y)
dp
= k p(x) = kx + A
dx
d2vx k
=
v x (y) = k y 2 + By + C
2
dy 2
v x (y)
xy = yx =
kx A k y + B
0
( x, y) = y + B kx A
0
0
kx A
Donde las constantes de estas expresiones ( k , A , B , C ) sern diferentes para cada fluido.
180
10 Mecnica de fluidos
Problemas Resueltos
y =0
= p2
y =0
fluidos.
CONDICIONES DE CONTORNO PARA LAS TENSIONES
La condicin de continuidad en tensiones que hay que imponer en la interfase de los dos
lquidos no afecta a todo el tensor , sino slo al vector de tracciones t . Deber
verificarse:
t1
y =0
= t 2
y =0
n 1 = [0 1 0]T
Teniendo en cuenta que la normal n es la normal exterior: 2
T
n = [0 1 0]
y=0
= 2xy
y=0
k 1 x + A1 = g t cos , x
A 1 = g t cos
k1 = k 2 = 0
k 1 x + A1 = k 2 x + A 2 , x
A1 = A 2 = g t cos
y=0
1 B1 = 2 B 2
k 1 = k 2 = 0
v 2x (y) =
v y + 1
h 2
1+ 1
2
1
2
v y + 1
1+ 1
2
p 1 = p 2 = g t cos = ctte.
1xy = 2xy = 1
h 1 + 1
2
= ctte
181
10 Mecnica de fluidos
Problemas Resueltos
fluido
x
1 + 1 h
2
v = h 1 + 1 (f g t sin )
1 2
Para calcular el caudal que atraviesa la superficie B - B , basta con tener en cuenta que
los fluidos son incompresibles y por tanto el caudal que pasa es el mismo si
consideramos el segmento recto que une B y B :
q=
BB' curva
v ndS =
BB' recta
v ndS =
v x (y)dy
q = v h 1 +
2 1 + 2
pide:
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
182
10 Mecnica de fluidos
Problemas Resueltos
a) Determinar los campos de velocidad y presin en la zona del fluido de la Figura 2 (zona
ABCD ), salvo las constantes de integracin.
b) Indicar y justificar las condiciones de contorno necesarias para determinar el valor de
las constantes de integracin.
c) Obtener dichas constantes y las correspondientes expresiones de los campos de
velocidad y presin.
d) Obtener la expresin de las tensiones en la zona ABCD del fluido.
e) Suponiendo que el campo de tensiones y en la superficie A - A es uniforme e igual
al del punto A , demostrar que existe una relacin entre la fuerza por unidad de
longitud F , aplicada sobre el mbolo, y la velocidad de descenso del mismo, & ( t ) , de la
forma F = & ( t) . Obtener el valor de .
y
F
a
B'
h
A'
A
O
m
x
L
Figura 1
Figura 2 (detalle)
HIPTESIS:
1) Fuerzas de masa despreciable en el fluido.
2) Peso del mbolo despreciable.
3) Flujo estacionario.
4) Presin atmosfrica despreciable.
Resolucin:
a) HIPTESIS AADIDAS:
Nada depende de la direccin z perpendicular al papel, pues el problema est indefinido
en esta direccin. Considrese una situacin bidimensional, siendo, en principio por
hiptesis:
T
v = [v x (x, y ) v y (x, y ) 0]
Por otro lado, sobre las paredes AB y CD , debe cumplirse: v x = 0 , por condicin de
impenetrabilidad (un fluido no puede penetrar en un slido).
Sin embargo, es conveniente introducir otra hiptesis ms para simplificar el problema.
Se trata de una hiptesis aproximada, no exacta. Se supondr que v x = 0 en toda la zona
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
183
10 Mecnica de fluidos
Problemas Resueltos
p = p(x, y )
= 0 v y = v y (x )
p = p(y)
2 v y
p
+
0=
y
x 2
0=
2vy
p
=
y
x 2
2vy
x 2
=K
v y
x
= K x + C2
v y (x ) = 1 K x 2 + C 2 x + C 3
2
184
10 Mecnica de fluidos
Problemas Resueltos
y=m+h
= p atm = 0
v ndS
S
El mbolo baja a una velocidad & y, por tanto, el mbolo se introduce en el fluido
ganndole espacio a ste. Es decir, es como si entrase una seccin (la del mbolo) a cierta
velocidad, empujando al fluido hacia fuera. Por esta razn se puede definir como un caudal
entrante (velocidad superficie):
Q entrante = & L
En cambio, el caudal saliente de fluido a travs de los espacios entre el mbolo y las
paredes laterales se puede calcular segn la definicin general de caudal:
Q saliente = 2
Sa
v ndS = 2 v y (x )dx
0
2 v y (x )dx = & L
0
v y (x ) x = 0 = 0
v y (0 ) = C 3
C3 = 0
v y (x ) x =a = &
v y (a ) = 1 K a 2 + C 2 a
2
&
C2 = K a
a 2
3
2
a
a
2 v y (x )dx = 2 K x 2 + C 2 x dx = 2 K a + C 2 a = & L
0
0 2
2
2 3
K=
6 &
(a + L )
a3
&
C 2 = 2 + 3 L
a
a
p(y ) y = m + h = 0 K (m + h ) + C1
C1 =
6 &
(a + L )(m + h )
a3
185
10 Mecnica de fluidos
Problemas Resueltos
6 &
(a + L )(m + h y )
a3
&
v y (x ) = 33 (a + L )& x 2 + 2 + 3 L x
a
a
a
p = p(y ) =
d) Para obtener las tensiones en la zona ABCD del fluido, basta con utilizar las
expresiones correspondientes a las Componentes del Tensor de Tensiones para Fluidos
Newtonianos en Coordenadas Cartesianas. Utilizando las expresiones del apartado c) de la
velocidad y la presin se obtiene lo siguiente:
x = p ; y = p ; z = p
xy = yx =
v y
x
; xz = zx = 0 ; yz = zy = 0
p
v
y
=
x
donde
v y
x
v y
x
p
0
= & 6 x3 (a + L ) + 3 L2 + 2
a
a
a
COMENTARIO:
De hecho, la hiptesis de flujo estacionario, mientras el mbolo baja, es una hiptesis
excesiva y no rigurosa, pues en algn instante el mbolo llegar hasta el punto ms
bajo y el flujo variar. Para poder aplicar sta hiptesis, m + h debera ser una longitud
muy grande o & una velocidad muy baja.
1 = xy x=a
Criterio de
signos positivo
= y
2 = xy
x =L+ a
Se van a calcular las tensiones correspondientes para obtener las fuerzas resultantes y
aplicar equilibrio.
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
186
10 Mecnica de fluidos
Problemas Resueltos
A=
p(y ) y = m = Km + C1
* =
6
(a + L )h&
3
a
En las superficies laterales, por simetra se tiene que 1 = 2 , por lo tanto slo se
calcular 1 :
1 = 2 = xy
x =a
& L
= & 6 a3 (a + L ) + 3 L2 + 2 1 = 2 =
3 + 4
a
a a
a
a
F = &
=
2h L2
L
3 2 + 6 + 4
a a
a
187
10 Mecnica de fluidos
Cuestiones Propuestas
CUESTIONES PROPUESTAS
CP 10-1 Escribir las ecuaciones de gobierno del problema general de Mecnica de Fluidos,
enunciando cada ecuacin y caracterizando el nmero de ecuaciones y de incgnitas.
Particularizar para las ecuaciones del problema mecnico cuando ste est desacoplado del
problema trmico.
CP 10-2 Escribir las ecuaciones de gobierno de un problema de Mecnica de fluidos sealando
el nmero de ecuaciones disponibles y las incgnitas del problema. Hacer lo mismo para el
caso particular en que el fluido sea Newtoniano, indicando cundo los problemas mecnico y
trmico estn desacoplados.
CP 10-3 Deducir las siguientes ecuaciones:
a) Ecuacin fundamental de la hidrosttica (equilibrio).
b) Ecuacin de equilibrio de un fluido perfecto (ecuacin de Euler).
CP 10-4 Demostrar el teorema de Arqumedes.
CP 10-5 Deducir la expresin del Teorema de Bernouilli para un fluido perfecto,
incompresible y en rgimen estacionario, bajo fuerzas gravitatorias. Interpretar grficamente el
resultado.
CP 10-6 Deducir la ecuacin de Bernouilli para un fluido perfecto, incompresible y en rgimen
estacionario, indicando el significado de cada trmino y el significado de cada trmino.
NOTA: Utilizar la igualdad: v v = 2 v + 1 v 2
2
2 t
donde
P=
(p )dp
0
= gz
188
10 Mecnica de fluidos
Cuestiones Propuestas
Fluido incompresible.
Fluidos con viscosidad volumtrica nula.
Fluidos perfectos.
Hidrosttica.
189
10 Mecnica de fluidos
Cuestiones Propuestas
A
F
H
HIPTESIS:
1) Se desprecia el peso de la compuerta
CP 10-16 Determinar el valor mnimo de la fuerza F que hay que hacer en el punto B de la
compuerta AB de la figura articulada en el punto A, para que no se produzca prdida de agua.
c
A
b
HIPTESIS:
1) Se desprecia el peso de la compuerta.
y
R
x
O
HIPTESIS:
1) Se desprecia el peso de la compuerta.
HIPTESIS:
1) Se desprecia el peso de la compuerta.
F
O
B
190
10 Mecnica de fluidos
Cuestiones Propuestas
CP 10-19 En la pared de un depsito existe una llave de cierre que gira alrededor del punto O
de la figura. Calcular las fuerzas y el momento resultante, por unidad de espesor, que el fluido
ejerce sobre la llave.
e2
HIPTESIS:
1) Se desprecia
compuerta.
el
peso
de
e1
la
R
O
e3
CP 10-20 Calcular el peso del lastre W que hay que poner en el fondo del cajn de la figura,
de peso W , para que ste flote en equilibrio estable.
NOTA: Los pesos son por unidad de espesor.
W
W
h/2
HIPTESIS:
1) La densidad del agua es .
S1
PA
vA 0
B
S2
PB = Patm
vB
191
10 Mecnica de fluidos
Cuestiones Propuestas
y
H
x
a
R1
R2
193
10 Mecnica de fluidos
Problemas Propuestos
PROBLEMAS PROPUESTOS
PP 10-1 La figura muestra la seccin longitudinal de una canal con una compuerta por el que
circula agua en rgimen estacionario. Para secciones transversales suficientemente alejadas de la
compuerta como la 1, aguas arriba, y la 2, aguas abajo, puede considerarse un campo de
velocidades horizontal y uniforme y una distribucin de presiones hidrosttica. Se pide:
a) Obtener los valores de la velocidad en las secciones 1 y 2.
b) Calcular el valor del empuje E que ejerce el agua sobre la compuerta (por unidad de
anchura del canal).
c) Admitiendo que la distribucin de presiones p 0 sobre cada punto de la solera del canal es
igual al peso de la columna de agua por encima de la misma ( p 0 =g h ) , obtener la
distancia d a la solera del punto de aplicacin de dicho empuje.
NOTA: Todos los resultados deben darse en funcin de H 1 , H 2 y la densidad del agua
exclusivamente.
E
compuerta
1 H
1
solera
2
H2
p 0 =g h
HIPTESIS:
1) Se considera al agua como un fluido perfecto incompresible.
2) Presin atmosfrica despreciable.
3) Fuerzas de rozamiento agua-solera y agua-compuerta despreciables.
4) No se considera el peso propio de la compuerta.
194
10 Mecnica de fluidos
Problemas Propuestos
c) Las fuerzas, por unidad de longitud, F1 y F2 ejercidas sobre cada cilindro para que se
produzca el movimiento descrito.
d) La energa disipada, por unidad de tiempo y de longitud, debida a los efectos viscosos (en
funcin de las constantes de integracin).
F2
v2
R2
A
R1
F1
v1
v2
HIPTESIS:
1. Rgimen estacionario.
2. Se desprecian las fuerzas de masa y las de inercia.
3. Debido al carcter indefinido de la longitud se considerar que el flujo (y sus propiedades)
no varan en direccin longitudinal.
PP 10-3 Una placa de longitud indefinida y espesor 2a separa dos fluidos newtonianos
incompresibles que se mueven entre dos contornos de longitud indefinida situados a una
distancia h de la placa (ver figura). La placa y el contorno superior se mueven a velocidades
v/2 y v , respectivamente. Se pide:
a) Los campos de velocidades, presiones y tensiones (en funcin de las constantes de
integracin).
b) Las constantes de integracin, aplicando las condiciones de contorno pertinentes.
c) Las fuerzas por unidad de superficie F1 y F2 ejercidas sobre la placa y el contorno superior
necesarias para que se produzca el movimiento descrito.
d) La energa disipada, por unidad de tiempo y de superficie perpendicular al plano del dibujo,
debida a los efectos viscosos.
F2
v
B
h
a
a
h
F1
A
v/2
HIPTESIS:
195
10 Mecnica de fluidos
Problemas Propuestos
d
H
p=0
a
a
L
HIPTESIS:
1) El flujo se supone paralelo al plano x - y .
2) Se desprecian las fuerzas de inercia.
3) Se consideran Q , V y H por unidad de longitud en direccin z .
4) El peso de la placa y la presin atmosfrica son despreciables.
PP 10-5 Un slido se peso W desliza con una velocidad v sobre una superficie horizontal
lubricada mediante un fluido Newtoniano de viscosidad . La superficie inferior del slido
est inclinada, con una pendiente muy pequea , respecto a la horizontal. Se pide:
a) Obtener el campo de velocidades en el fluido en funcin de v , h(x) y q , donde q es el
caudal de fluido que fluye por unidad de anchura a travs de una superficie vertical entre el
suelo y el slido.
b) Obtener la distribucin de presiones en el fluido en funcin de , , q , v , h(x) , h 1 y
P0 , sabiendo que en el punto A la presin es la atmosfrica P0 .
c) Obtener el caudal q en funcin de v , h 1 y h 2 sabiendo que la presin en el punto B el
P0 .
d) Obtener la nueva expresin de la presin en funcin de , , v , h(x) , h 1 , h 2 y P0 .
e) Obtener las componentes del tensor de tensiones en funcin de , , v , h(x) , h 1 , h 2 y
P0 .
196
10 Mecnica de fluidos
Problemas Propuestos
f) Obtener la expresin del peso del slido W y de la fuerza horizontal F , necesaria para
moverlo en funcin de , , v , h(x) , h 1 y h 2 .
g) Obtener el cociente = F/W (coeficiente de rozamiento efectivo).
y
F
W
B
A
h1
h2
h(x)
L=
h 2 h1
HIPTESIS:
1) Fluido incompresible en rgimen estacionario.
2) Peso del fluido despreciable.
3) Flujo muy lento (fuerzas de inercia despreciables).
v x
0
x
g(z )]T
,a >0
Se pide:
a)
b)
c)
d)
197
10 Mecnica de fluidos
Problemas Propuestos
1
2
PP 10-8 Un disco de radio R gira con velocidad angular constante a una distancia a de una
superficie horizontal. Entre el disco y dicha superficie se encuentra un fluido de
comportamiento Newtoniano con viscosidad . Se pide:
a) La expresin del campo de velocidades en el fluido antes de aplicar las condiciones de
contorno.
b) Dicha expresin despus de aplicar las condiciones de contorno.
c) La expresin de la presin y de la tensin tangencial z .
d) El valor del momento M que se ha de aplicar en el eje del disco para mantener el
movimiento.
z
M
R
HIPTESIS:
1) Pueden despreciarse las fuerzas de inercia por ser el movimiento suficientemente lento.
2) Fluido incompresible.
3) No se considera el efecto de las paredes laterales (se desprecian los efectos del rozamiento
fluido-pared lateral).
4) Se considerar que el campo de velocidades vara linealmente con la distancia a la superficie
inferior.
2000-2004 by: X. Oliver & C. Agelet de Saracbar
E.T.S. Enginyers de Camins, Canals i Ports (Universitat Politcnica de Catalunya)
198
10 Mecnica de fluidos
Problemas Propuestos
HIPTESIS:
1) Fluido incompresible.
2) No se consideran las fuerzas de volumen.
3) Rgimen estacionario.
PP 10-10 Un cilindro macizo de longitud indefinida y radio R gira con una velocidad
constante, concntricamente con otro cilindro hueco de radio interior 2R que a su vez gira a
velocidad angular constante + a . En el espacio entre ambos cilindros hay un lquido
Newtoniano incompresible de viscosidad . Se pide:
a) Determinar el campo de velocidades del fluido
en rgimen estacionario.
b) Determinar los momentos que hay que aplicar
sobre ambos cilindros para mantenerlos a las
velocidades indicadas (indicando los sentidos).
c) Determinar el valor de la presin
termodinmica en las siguientes condiciones:
1. Los dos cilindros giran a la misma
velocidad.
2. El cilindro interior no se mueve.
3. El cilindro exterior no se mueve.
2R
+a
199
10 Mecnica de fluidos
Problemas Propuestos
HIPTESIS:
1) Movimiento estacionario.
2) Movimiento suficientemente lento para despreciar las
fuerzas de inercia.
3) Flujo incompresible.
R1
R2
HIPTESIS:
1) Fuerzas msicas despreciables
en el fluido.
2) Peso del mbolo despreciable.
3) Flujo estacionario.
4) Presin atmosfrica despreciable
B'
C
&
R2
E
R1
F
E'
D'
A'
B'
C'
a
200
10 Mecnica de fluidos
Problemas Propuestos
PP 10-13 Un fluido Newtoniano incompresible fluye por el espacio entre dos esferas,
concntricas y muy cercanas, al ser introducido por un pequeo orificio en la parte superior de
la esfera exterior. Sea P = PA PB la diferencia de presin entre los puntos A y B de la
figura y PB la presin atmosfrica. Se pide:
a) Plantear el campo de velocidades y presiones y las correspondientes ecuaciones
diferenciales y condiciones de contorno.
b) Obtener el campo de velocidades y presiones en funcin de P .
c) Calcular el caudal que atraviesa la superficie horizontal M - N en funcin de P .
A
M
kR
N
R
R/2
OBSERVACIN:
sen d=ln tg 2
1
HIPTESIS:
1) Movimiento muy lento (fuerzas de inercia despreciables).
2) Fuerzas de volumen despreciables.
3) No se considera la perturbacin del flujo
producida por el orificio de entrada.
4) El ngulo es muy pequeo (puede tomarse
tg = ).
5) Rgimen estacionario.
A
W
R
201
10 Mecnica de fluidos
Problemas Propuestos
R
v1 v1 p1 1
1 11
v2 p2
2 2
HIPTESIS:
1) Presin atmosfrica despreciable.
2) Se sugiere considerar campos de velocidades para los fluidos de la forma:
v(r , , ) = [v r
T
v ] = 0 0 Ar + B2 sin
debiendo justificarse que se cumplen todas las ecuaciones de gobierno del problema.
3) Despreciar las fuerzas de inercia y de masa en el exterior de la esfera, pero no en el interior.
202
10 Mecnica de fluidos
Problemas Propuestos
PP 10-16 Sobre un fluido perfecto que est en reposo en el instante t = 0 actan unas fuerzas
por unidad de masa b = [y 0 0]T (descripcin espacial). La distribucin (espacial) de
presiones es p(x, t ) = (t )ze -t y la distribucin espacial de la densidad es (x, t ) = (t ) .Se pide:
a) Las descripciones espaciales del campo de
aceleraciones y del campo de velocidades y las
ecuaciones de las trayectorias (en forma
cannica) y de las lneas de corriente.
b) Determinar la funcin (t ) sabiendo que en el
instante inicial (0 ) = 0 .
c) Determinar el valor del parmetro para que
el flujo sea irrotacional. Para dicho caso
determinar los potenciales de fuerzas msicas
( b = ) y de velocidades ( v = ).
Demostrar que en estas condiciones se
cumple la ecuacin:
p
+ + 1 v 2 +
=0
t
2
R
y
, x , t
203
11 Principios variacionales
Cuestiones Resueltas
11 Principio s variacionales
CUESTIONES RESUELTAS
CR 11-1 A partir de la expresin del principio de los trabajos virtuales:
V0
: dV0 = bu dV0 +
V0
t udS
u u=0 en u
W( )
W = ij ij = :
G (u )
(u)
G = t u
= b u
Ahora se puede reescribir la expresin del principio de los trabajos virtuales como:
WdV0 +
V0
V0
V0
WdV0 +
dV0 + Gd = 0 u u=0
V0
en
u u=0 en
dV0 + Gd = 0
V0
WdV0 +
V0
dV0 + Gd
se puede concluir:
U = 0
u u=0 en
donde:
2
C ijkl = W
ij kl
Si se calcula:
U(u + u) =
U(u) =
V0
V0
1 (u + u ) : C : (u + u)dV
0
2
1 (u) : C : (u)dV
0
2
V0
b udV0
b (u + u)dV
t ud
0
V0
t * (u + u)d
resulta:
1 (u) : C : (u)dV + 1 (u ) : C : (u )dV +
0
0
V0 2
2
+ 1 (u) : C : (u)dV0 b udV0 t * ud
V0 2
V0
U(u + u) U(u ) =
V0
y utilizando:
V0
: dV0 +
V0
1 (u) : C : (u)dV
0
2
V0
b udV0 t * ud
V0
: dV0
V0
b udV 0 t * ud
se tiene que:
U (u + u) U (u ) =
V0
1 (u ) : C : (u) dV
0
2
2
Y considerando que el tensor C ijkl = W es definido positivo:
ij kl
U(u + u) U(u) 0
( x x + y y + z z + xy xy + xz xz + yz yz )d
205
11 Principios variacionales
Cuestiones Resueltas
a) Tensin plana.
b) Deformacin plana.
Resolucin:
En tensin plana se cumple:
xz = yz = z = 0
W int =
( x x + y y + xy xy )d
( x x + y y + xy xy )d
207
11 Principios variacionales
Cuestiones Propuestas
CUESTIONES PROPUESTAS
CP 11-1 Indicar y justificar cul es el mbito de validez y aplicacin en Mecnica de Slidos de
los siguientes principios variacionales:
a) Principio de los Trabajos Virtuales
b) Principio de la Energa Potencial mnima.
CP 11-2 Enunciar el principio de minimizacin de la energa potencial, indicar su campo de
validez y deducirlo a partir de la expresin del Principio de los Trabajos Virtuales.
CP 11-3 Demostrar que el principio de los trabajos virtuales es una condicin necesaria para el
equilibrio de un medio continuo.
CP 11-4 Demostrar que el Principio de los Trabajos Virtuales es condicin suficiente para que
un sistema de tensiones y desplazamientos sea solucin de las ecuaciones de equilibrio.