Problemas resueltos de Mecánica de Fluidos
Capítulo 1: Fluidostática
Marcos Lema & Anne Gosset
6 de abril de 2021
1. Una presa de anchura b forma un embalse de agua de profundidad H.
Calcular la distribución de presiones sobre la presa, la fuerza resul-
tante y el momento que ha de soportar la pared. Determinar donde
se localiza el centro de presiones, de tal forma que el momento que
ocasiona la fuerza resultante sobre este punto sea igual al calculado
previamente.
Solución
a) Distribución de presiones
La distribución de presiones permite obtener el valor de la presión en función de las
variables del problema. Para ello tenemos que integrar la ecuación fundamental de la
fluidostática:
∇P = ρf~m
En este caso, f~m = −g k,
b por lo tanto:
!
∂b ∂ b ∂
i+ j + kb P = −ρg kb
∂x ∂y ∂z
1
Separando cada una de las componentes de la ecuación para posteriormente integrar,
obtenemos:
∂P
=0
∂x
∂P
=0
∂y
∂P
= −ρg
∂z
Las dos primeras integrales son directas de tal forma que: Px = cte. y Py = cte. Este
resultado indica que la presión no cambia cuando varía la coordenada x e y. Para
integrar la ecuación en derivadas parciales de la componente z, procedemos primero a
hacer el cambio a derivadas totales:
∂P
dP = dz
∂z
Por lo que la ecuación a integrar es:
Z P2 Z z2
dP = −ρg dz
P1 z1
Donde P1 , P2 , z1 , z2 son los límites de integración que definiremos posteriormente. El
resultado de la integral es:
P2 − P1 = −ρg (z2 − z1 ) + C
Es fácil deducir que la constante de integración es 0. Además, como la presión en
la superficie libre es igual a la atmosférica, podemos sustituir uno de los límites de
integración por los siguientes valores:
P1 = Pa y z1 = H
Eliminando los subíndices del otro límite de integración y reordenando la cuación
resultante, obtenemos la distribución de presiones buscada.
P = Pa + ρg (H − z)
En esta ecuación, la presión tiene una parte constante, debido a la presión atmosférica,
y una parte variable debido al líquido y que depende la profundidad definida por la
coordenada z. Si representamos gráficamente dicha evolución del valor de la presión,
obtenemos los prismas de presiones que se muestran en la siguiente figura. En dicha
figura también se ha incluído el prisma de presiones en la cara seca de la presa, la cual
está únicamente bajo el efecto de la presión atmosférica.
2
b) Fuerza resultante
La fuerza resultante debido la presión que el agua ejerce sobre la presa se obtiene
integrando la siguiente ecuación:
dF~ = −P dA n
b
En la ecuación anterior, P toma el valor de la distribución de presiones obtenida ante-
riormente, donde también se tiene que incluir el valor de la presión en la cara seca de
la pared..
dF~ = − [Pa + ρg(H − z)] dA(−bi) − Pa dA(bi)
La fuerza resultante debido a la presión atmosférica a uno y otro lado de la pared se
anulan. Sabiendo que dA = dzdy, por lo que la ecuación a integrar es:
Z F Z H Z b
dF = ρg (H − z)dz dy
0 0
" #H " #b
z2
F = ρg (Hz − y
2 0 0
1
F~ = ρgbH 2 (bi)
2
c) Momento
El cálculo del momento se realiza con la siguiente ecuación vectorial:
~ = ~r × dF~
dM
donde ~r es el vector distancia medido desde el centro de momentos, que en el presente
caso coincide con el origen de coordenadas. Por lo tanto:
~r = z kb
3
La expresión del diferencial de fuerza es la misma que la obtenida previamente:
dF = ρg(H − z)dzdy(bi)
Para el cálculo del producto vectorial, calculamos el vector unitario resultante, el cual
se obtienes resolviendo el siguiente tensor:
bi
kb j
b
kb × bi = 0 0 1 = bj
1 0 0
Por lo tanto:
~ = ~r × dF~ = ρg(Hz − z 2 )dzdy(bj)
dM
Integramos ahora para calcular el momento:
Z M Z H Z b
2
dM = ρg (Hz − z )dz dy
0 0
" #H " #b
Hz 2 z 3
M = ρg − y
2 3 0 0
Una vez sustituidos los límites de integración, se obtiene el valor del momento:
3
~ = ρgb H bj
M
6
Para obtener el centro de presiones,ZP , suponemos que la fuerza resultante se concentra
en dicho punto, de tal forma que el momento es el mismo. Por lo tanto:
H
FR ZP = M → ZP =
3
El centro de presiones, ZP , es un punto ficticio que nos permite calcular de forma
sencilla el momento que produce la fuerza resultante. Dicho punto está situado a H/3
del fondo, tal como se indica en la figura.
4
2. Calcular el peso que ha de tener un semicilindro hueco, de radio R
y longitud b, de tal forma que cuando está flotando en agua la línea
de flotación coincide con el borde superior, tal como se muestra en
la figura.
Solución
Para la resolución de este poblema se van a plantear dos métodos distintos.
a) Integrar la ecuación de fuerza.
~ = −P dA n
dF b
A la vista de los datos indicados en la figura, sabemos que:
P = ρ g r senθ
dA = b r dθ (En coordenadas cilíndricas)
~ = −b ρ g r2 senθ dθ n
dF b
Proyectando F~ en su componente horizontal y vertical, el ser el semicilindro una figura
simétrica la componente horizontal se anula, por lo que únicamente integraremos la
componente vertical: Z Z 0
dF = b ρ g r2 sen2 θ dθ
−π
Teniendo en cuenta que:
1 cos(2θ)
sen2 θ = −
2 2
5
! 0 0 !
Z 0
1 cos(2θ) 1 1 1
Z
2
dF = b ρ g r − dθ = b ρ g r2 θ − sen (2θ)
−π 2 2 2 −π 2 2 −π
Sustituyendo los límites de integración, obtenemos:
π
Fr = b ρ g r 2
2
Al ser el semicilindro una figura simétrica con igual distribución de presiones a ambos
lados, el centro de presiones estaría situado sobre el eje de simetría y la fuerza resultante
tiene la dirección y sentido del vector unitario (k),
b tal como se indica en la figura.
π
F~r = b ρ g r2 kb
2
b) Aplicando el principio de Arquímedes.
Una forma sencilla de resolver este problema es haciendo uso del principio de Ar-
químedes, donde únicamente tenemos que calcular el volumen de agua desplazado por
el semicilindro, y multiplicarlo por la densidad del agua y la gravedad para obtener el
peso. Dicho valor es igual a la fuerza vertical que experimenta el cilindro:
π
V = r2 b
2
π
Fv = V · ρ g = b ρ g r 2
2
Que es el mismo valor que hemos obtenido resolviendo la ecuación integral.
Como el enunciado nos pide calcular el el peso que ha de tener el semicilindro, este
tiene que ser igual a la fuerza vertical ejercida por el agua, F~v . Si el peso fuese menor,
el semicilindro emergería, desplazando un volumen menor de agua, lo que resultaría
en una menor fuerza vertical que igualaría el nuevo peso. Si el peso fuese mayor, el
semicilindro se hundiría.
6
3. La compuerta de la figura, de forma triangular y ancho b, cierra un
depósito de agua y puede pivotar en el punto donde se ha situado el
origen. Calcular la fuerza que es necesario aplicar en el punto A, para
que la compuerta no se abra bajo la influencia del agua contenida en
el depósito.
Solución
Para la resolución de este poblema se van a plantear dos métodos distintos.
a) Prisma de presiones sobre la compuerta.
Teniendo en cuenta que el origen está en la superficie libre y que la única fuerza
másica es debida a la gravedad, la distribución de presiones resulta ser:
P = Pa − ρ g z
El prisma de presiones aparece representado en la figura, normal a la superficie. Para
el cálculo de la fuerza eliminamos Pa , al estar todas las caras de la compuerta bajo la
influencia de la presión atmosférica. La fuerza resultante la calculamos con el volumen
del prisma de presiones:
1
F~r = − ρ g H b Ln
b
2
7
H
L=
cos(45)
donde L es la longitud del plano inclinado y n b es el vector unitario normal a dicho
plano. El centro de presiones está situado a L/3 medido desde A, que es el punto de
máxima presión en el agua.
1 1 H
Zp = L =
3 3 cos(45)
Para calcular FA planteamos un balance de momentos, para lo que descomponemos la
fuerza resultante sobre la compuerta en su componente vertical y horizontal:
1
F~h = F~r sen(45)bi = ρ g b H 2bi
2
1
F~v = F~r cos(45)kb = ρ g b H 2 kb
2
El balance de momentos resultanto respecto al punto de pivote sería:
1 2
FA H + Fv · L cos(45) = Fh · L cos(45)
3 3
1 1 1 2
FA H + ρ g b H 2 · H = ρ g b H 2 · H
2 3 2 3
Obteniendo:
1
F~A = ρ g H 2 b(−bi)
6
b) Proyectanto el valor de la presión.
Otra forma de resolver este tipo de problemas es proyectando el valor de la presión en
la superficie en contacto con el agua sobre la dirección vertical y horizontal, tal como
se muestra en la figura. De esta forma, se puede obtener fácilmente la fuerza resultante
en la proyección vertical y horizontal calculando el volumen de ambos prismas de pre-
siones. El centro de presiones en la proyección horizontal estaría situado a 2/3H del
punto de pivote y el de la proyección vertical a 1/3H.
1 2
F~h = ρ g H 2 b bi → Zp = H
2 3
1 2
F~v = ρ g H 2 b kb → Zp = H
2 3
Aplicamos ahora el balance de momentos para despejar el valor de FA :
1 2
FA H + Fv · H = Fh · H
3 3
1
F~A = ρ g H 2 b(−bi)
6
8
4. Un ascensor de base 1x1 m2 está inundado de agua hasta una altura
H = 2m. Tanto las paredes como la puerta del ascensor son estancas,
por lo que no hay fugas de agua al exterior, donde la presión es la
atmosférica. La puerta tiene unas dimensiones h = 1, 8 m de alto por
b = 1 m de ancho, y se abre gracias a unas bisagras instaladas en un
lateral de la misma. Si inicialmente el ascensor está detenido, se pide
calcular:
a) La distribución de presiones en el volumen fluido y la fuerza resultante sobre el fondo
del ascensor.
b) La fuerza F que hay que aplicar a la puerta, en el punto que se indica en la figura, para
evitar que gire con respecto a las bisagras, y mantenerla así cerrada.
c) En un momento dado, el ascensor empieza a subir con una aceleración de 5m/s2 :
Calcular la nueva distribución de presiones y la fuerza resultante sobre el fondo.
d) Debido al peso del conjunto, el cable del ascensor se rompe, y éste cae libremente:
Calcular la distribución de presiones y la fuerza resultante sobre el fondo.
e) Calcular la fuerza F para mantener la puerta cerrada en este caso.
9
Solución
a) Situando el sistema de referencia en el fondo del ascensor, al existir únicamente la fuerza
másica debido a la gravedad, la distribución de presiones se puede calcular fácilmente
como:
P = Pa + ρ g (H − z)
Para el cálculo de la fuerza resultante sobre el fondo del ascensor, como la distribución
de presiones para z = 0 es constante y el area del fondo es A = b2 , obtenemos:
P = Pa + ρ g H
F~r = P · A = ρ g H b2 (−k)
b
Nótese que para el cálculo de la fuerza no hemos tenido en cuenta la presión atmosférica,
al anularse el valor de esta con la presión atmosférica a la que están sometidas las caras
exteriores del ascensor.
b) En este caso, la distribución de presiones sobre la puerta ya no es constante y, además,
la altura de la puerta es menor que la altura de agua en el interior del ascensor. Por la
tanto, la representación gráfica de la distribución de presiones sobre la puerta tendrá
una parte constante y otra variable, tal como se indica a continuación:
En la figura anterior ya no hemos representado la presión atmosférica, porque a la hora
de calcular la fuerza resultante sobre la puerta tampoco la tendremos en cuenta por la
misma razón explicada previamente. La fuerza resultante la calculamos con el volumen
del prisma de presiones.
F~cte = ρ g b h (H − h)bi
1
F~var = ρ g b h2 )bi
2
El centro de presiones estará situado a una altura de h/2 para Fcte y de h/3 para
la Fvar , ambos a una distancia b/2 del punto de pivote que son las bisagras. Como
10
la fuerza F se aplica a una distancia b del punto de pivote, el balance de momento
necesario para obtener el valor F se plantea como:
F b = (Fcte + Fvar )b/2
1 1
F~ = ρ g b h (H − h) + ρ g b h2 (−bi)
2 4
c) En el vector fuerzas másicas, la gravedad y la acelaración tiene igual dirección y sentido,
ya que el vector aceleración apunta en sentido contrario al movimiento del ascensor.
f~m = (g + a)(−k)
b
La nueva distribución de presiones sería:
P = Pa + ρ (g + a) (H − z)
Una vez más, la distribución de presiones sobre el fondo es constante y la presión
atmosférica se anula al haber la misma presión en el exterior del ascensor. Por lo
tanto:
F~r = P · A = ρ (g + a) H b2 (−k)
b
d) En este caso, la aceleración del ascensor coincide con el valor de la gravedad, por lo
que el vector fuerzas másicas es nulo:
f~m = (g − g) kb = 0
Por lo tanto, la distribución de presiones es constante e igual a:
P = Pa
Lo que significa que la fuerza resultante sobre el fondo del ascensor es nula.
e) De nuevo, si la distribución de presiones en el líquido contenido en el ascensor es igual
a la presión exterior, la fuerza necesaria para mantener la puerta cerrada es nula.
11
5. Un tanque cilíndrico de radio R está lleno de agua hasta una altura
H. Se pide calcular la distribución de presiones y la fuerza resultante
en el fondo cuando el tanque gira a una velocidad ángular ω.
Solución
Para obtener la distribución de presiones es necesario integrar la ecuación fundamental
de la fluidostática, ∇P = ρ f~m , sabiendo que en la presente configuración, el vector fuerzas
másicas depende de la gravedad y de la aceleración centrífuga. En coordenadas cilíndricas,
tenemos:
f~m = g(−k)
b + ω 2 r rb
Integrando la ecuación fundamental de la fluidostática, obtenemos:
Z P2 Z z2 Z r2
d P = −ρ g dz + ρ ω 2 r dr
P1 z1 r1
!
2 r2 2 r12
P2 − P1 = −ρ g(z2 − z1 ) + ρ ω −
2 2
La distribución de presiones nos permite calcular la presión en función de las variables del
problema, que en el presente caso son z y r. Por lo tanto, en la ecuación anterior necesitamos
sustituir uno de los límites de integración por condiciones de contorno conocidas. Si asig-
namos estos valores conocidos a las variables P1 , r1 y z1 , observamos que no existe ningún
punto del volumen fluido donde conozcamos las tres variables. Por ejemplo, sabemos que en
la superficie libre la presión es la atmosférica y conocemos los valores de r en el centro y en
el borde del tanque, pero no conocemos los valores de z en ningún caso, ya que no sabemos
cuando se eleva la superficie libre en los borde ni cuando se deprime en el centro. Por eso
necesitamos plantear una nueva ecuación, que será la conservación de la masa, ya que la
cantidad de agua cuando el tanque está en reposo y cuando gira es la misma, al no haber
desbordamiento.
Para ello, primero asignamos como valores de P1 , r1 y z1 los correspondientes al pun-
to de la superficie libre situado en el eje de simetría donde, tal como hemos dicho, solo
12
desconocemos el valor de z. Por lo tanto:
P1 = Pa ; r1 = 0 y z1 = z0
Reescribimos la distribución de presiones con estos valores y cambiamos los valores de P2 ,
r2 y z2 por los mas genéricos P , r y z.
!
2 r2
P = Pa − ρ g(z − z0 ) + ρ ω
2
Cuando el tanque está en reposo, el volumen de agua en el tanque es:
V = Π R2 H
Cuando el tanque gira, tenemos que integrar el volumen de agua desde la superficie libre
hasta el fondo del tanque. En coordenadas cilíndrica, el diferencial de volumen se expresa
como:
dV = 2Π z r dr
Donde z indica, para cada dr, la distancia desde el fondo del tanque hasta la superficie libre,
por lo que necesitamos conocer la curva que forma dicha superficie. Esto se puede obtener
sustituyento P = Pa en la distribución de presiones, ya que esta es la presión que tenemos
en toda la superficie libre. !
2
2 r
0 = −ρ g(z − z0 ) + ρ ω
2
ω 2 r2
z = z0 +
2g
Sustituyendo z en la ecuación de dV e integrando, obtenemos:
!
Z Z R
ω 2 r2
dV = 2 Π z0 + r dr
0 2g
!
R2 ω 2 R4
V = 2 Π z0 +
2 8g
Igualando al volumen inicial obtenemos el valor de z0 que buscábamos:
ω 2 R2
z0 = H −
4g
Por lo que la distribución de presiones en el tanque es:
!
ω 2 R2 r2
P = Pa − ρ g(z − H + ) + ρ ω2
4g 2
Para obtener la fuerza resultante en el fondo, integramos la ecuación del diferencial de fuerza:
~ = −P dA n
dF b
13
Sabiendo que la presión en el fondo del tanque la obtenemos con la ecuación de la distribución
de presiones cuando z = 0.
!
ω 2 R2 r2
P = Pa + ρ g(H − ) + ρ ω2
4g 2
Por lo tanto, en coordenadas cilíndricas:
Z R" !#
Z
ω 2 R2 r2
dF = Pa + ρ g(H − ) + ρ ω2 2 Π r dr
0 4g 2
" ! #
ω 2 R2 R4
F~ = Pa + ρ g(H − ) Π R2 + ρ ω 2 Π (−zb)
4g 4
6. Un depósito de base cuadrada de lado L, está lleno de agua hasta una
altura H. Cuando se somete a una aceleración, la inclinación de la
superficie libre es tal que el fondo queda parcialmente al descubierto,
tal como se indica en la figura. Bajo estas condiciones, se pide:
a) La distribución de presiones en el volumen fluido.
b) El centro de presiones en la pared opuesta a la dirección de la aceleración.
c) La fuerza resultante sobre el fondo del depósito.
Solución
a) Para obetner la distribución de presiones, integramos la ecuación fundamental de la
fluidostática teniendo en cuenta que el vector fuerzas másicas es el siguiente:
f~m = −abi − g kb
Por lo tanto, la ecuación a integrar es:
14
!
∂b ∂ b
i + k P = −ρabi − ρg kb
∂x ∂z
Integramos por separado cada una de las componentes de la ecuación anterior:
∂P
= −ρa
∂x
∂P
= −ρg
∂z
Haciendo el cambio a derivadas totales, obtenemos las dos ecuaciones a integrar:
Z P2 Z x2
dP = −ρa dz
P1 x1
Z P2 Z z2
dP = −ρg dz
P1 z1
Donde P1 , P2 , x1 , x2 , z1 , z2 son los límites de integración que, por el momento, que-
darán como incognitas. El resultado de la integral, sustituyendo las correspondientes
constantes de integración es:
P2 − P1 = −ρa (x2 − x1 ) − ρg (z2 − z1 )
En este caso, y debido a que la inclinación de la superficie libre deja al descubierto
parte del fondo del depósito, no disponemos de la información necesaria para sustituir
directamente uno de los límites de integración. Es por ello que tenemos que plantear
una segunda ecuación para calcular uno de estos límites, dejando el otro como las
variables de la distribuición de presiones. Esta segunda ecuación es la igualdad entre
el volumen de agua inicial, V0 y el volúmen de agua final, Vf , los cuales se pueden
calcular de la siguiente forma:
V0 = H L2
1
Vf = x2 z1 L
2
En el cálculo de Vf , hemos tomado z1 como la coordenada z de la superficie libre en la
pared vertical, y x2 como la coordenada x de la superficie libre en el fondo del tanque,
tal como se indica en la siguiente figura.
Igualando el volumen inicial y el final, obtenemos:
1
H L2 = x2 z1 L
2
Con esta segunda ecuación ya podemos calcular uno de los límites de integración. En
este caso, lo vamos a hacer para el punto 1 mostrado en la figura. En dicho punto,
15
sabemos que x1 = 0, P1 = Pa y el valor de z1 lo podemos obtener de la igualdad del
volumen inicial y final, esto es:
2H L2
z1 =
x2
Sustituyendo los valores del punto 1 en la distribución de presiones y dejando los límites
de integración en 2 como las variables de la distribución de presiones, obtenemos la
expresión buscada.
16
