Tema 3 Grado en Ingeniería Mecánica

INTEGRACIÓN DE SUPERFICIE

CONOCIMIENTOS PREVIOS

Para poder seguir adecuadamente este tema, se requiere que el alumno repase y ponga al día sus
conocimientos en los siguientes contenidos:
• Derivación de funciones de varias variables.
• Cálculo y propiedades de integrales dobles y triples.
• Manejo de ecuaciones paramétricas de curvas.
• Manejo de los sistemas de coordenadas cilíndricas y esféricas.
• Dibujo de curvas y superficies con Matlab y Dpgraph.

OBJETIVOS ESPECÍFICOS

Los objetivos específicos de este tema son:

1. Entender la construcción del elemento diferencial de superficie y su significado geométrico;
saber calcularlo para superficies expresadas en paramétricas, para superficies dadas por una
ecuación explícita y para superficies dadas por una ecuación implícita.

2. Entender la definición de integral de un campo escalar sobre una superficie y saber calcularla
dada la superficie y el campo. Hacer uso de sus propiedades.

3. Saber utilizar esta integral para calcular el área de una superficie.

4. Conocer sus interpretaciones físicas básicas: temperatura media, masa. Saber utilizar esta
integral para calcular promedios.

5. Entender la definición de integral de un campo vectorial sobre una superficie orientada y saber
calcularla dada la superficie y el campo. Conocer la relación entre esta integral y la de campos
escalares. Manejar el efecto del cambio de orientación de la superficie sobre el signo de la
integral.

6. Saber calcular, mediante integración, el flujo de un campo vectorial a través de una superficie.

7. Poder explicar el teorema de la divergencia de Gauss y saber usarlo para calcular una integral
de superficie sobre una superficie cerrada.

8. Poder explicar el teorema de Stokes y saber usarlo para calcular una integral de línea a lo largo
de una curva cerrada.
2 T3  INTEGRACIÓN DE SUPERFICIE

SUPERFICIES

1 Integral doble sobre rectángulos

Una superficie S en  3 puede venir dada de 3 formas:

a) Como la gráfica de una función explícita z  f x , y  donde x , y   D siendo D un

conjunto del plano S  x, y, f x, y  / x, y   D  .

b) Por medio de ecuaciones paramétricas, es decir

r u, v   x u, v  i  y u, v  j  z u, v  k u, v   D

c) De forma implícita como el conjunto de puntos que verifican F x , y, z   0 siendo F una

función escalar de 3 variables. S  x, y, z  / F x, y, z   0

Nota: Al final del tema se dan las ecuaciones y las gráficas de algunas superficies.

2 Plano tangente en un punto a una superficie

Un vector es tangente a una superficie S en un punto P a, b, c  de S si es el vector tangente

a una curva contenida en S y que pasa por P.
Un vector es ortogonal a una superficie S en un punto P a, b, c  de S si es ortogonal a todo
vector tangente a S en el punto P.
El plano tangente a una superficie en el punto P a, b, c  contiene a todos los vectores
tangentes a la superficie en dicho punto. El cálculo del plano tangente depende de la forma
de definir la función.

Caso 1. Si la superficie viene dada por la ecuación z  f (x , y ) , con (x , y )  D , siendo f de clase C 1

en D , el plano tangente y un vector normal en un punto cualquiera P x 0 , y 0 , z 0  con x 0 , y 0  en D
será:

• Plano tangente en P: z  f (x 0, y 0 )  fx x 0, y 0  (x  x 0 )  fy x 0, y 0  (y  y 0 )

• Un vector normal en P: N = (fx(x 0 , y 0 ),  fy(x 0 , y 0 ), 1)

Caso 2. Si la superficie viene dada por una función implícita F (x , y, z )  0 , con F de clase C1 el
plano tangente y un vector normal en un punto cualquiera P ( x0 , y0 , z0 ) con Fz x 0 , y 0 , z 0   0 será

• Plano tangente en P es:

Fx x 0 , y 0 , z 0 x  x 0   Fy x 0 , y 0 , z 0 y  y 0   Fz x 0 , y 0 , z 0 z  z 0   0
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
3

• 
Un vector normal en P es: N  Fx x 0 , y 0 , z 0 , Fy x 0 , y 0 , z 0 , Fz x 0 , y 0 , z 0  
Caso 3. Suponiendo que las ecuaciones paramétricas son de clase C 1 en D , es posible definir los
siguientes elementos de la superficie:

• Vectores tangentes a la superficie en un punto P (x (u 0 , v 0 ), y(u 0 , v 0 ), z (u 0 , v 0 ))

ru  (x u (u, v 0 ), yu (u, v 0 ), z u (u, v 0 )) y rv  (x v (u 0 , v ), yv (u 0 , v ), z v (u 0 , v ))

Estos vectores son tangentes, respectivamente, a las curvas de la superficie

v  v 0 , u  u 0 , que se cortan en (u 0 , v 0 ) .

• Vector normal en P : ru u 0 , v 0   rv u 0 , v 0  .

ÁREA DE UNA SUPERFICIE

3 Diferencial de superficie

Definición (Diferencial de superficie).- Se define el elemento dS como un parche

infinitesimal de superficie, cuya proyección sobre los planos coordenados en cada punto es,
dxdy dydz dxdz
dS   
cos  cos  cos 


siendo cos , cos , cos   los cosenos directores del vector normal a la superficie en
dicho punto.

En el siguiente apartado se muestra el proceso de cálculo de la diferencial de superficie dependiendo

de las ecuaciones con las que esté definida dicha superficie.
4 T3  INTEGRACIÓN DE SUPERFICIE

4 Cálculo de la Diferencial de una Superficie

SUPERFICIE DADA POR LA ECUACIÓN CARTESIANA EXPLÍCITA: z  f (x , y )

Si la superficie viene dada por la ecuación z  f (x , y ) , con (x , y )  D , siendo f de clase C1 en D

, se pueden definir

Un vector normal unitario a la superficie en

 
P x 0 , y 0 , f x 0 , y 0  :

N (fx(x 0 , y 0 ),  fy(x 0 , y 0 ), 1)
n =
N fx(x 0 , y 0 )2  fy(x 0 , y 0 )2  1

Vector normal al plano z = 0 en P:

k  (0, 0,1)

Ángulo que forman n y k en P: γ

En la figura se muestran, un parche de superficie de área S , una porción de plano tangente a éste
parche de área T y la proyección común de ambos, el rectángulo de área R .
El objetivo es encontrar una aproximación del área de ∆S que se proyecta sobre un rectángulo dado
del plano z  0 , de área R  x y .

Puesto que la superficie que mejor aproxima a S en el entorno de un punto es su plano tangente,
es la que utilizaremos en esta aproximación. Tomaremos el plano tangente que pasa por un punto
P0 (x 0 , y 0 , f (x 0 , y 0 )) y de él recortaremos la porción T que se proyecta sobre R , verificándose

1
T = R
cos 

Cuando las dimensiones de ∆R tiendan a cero, ∆T tenderá al elemento diferencial de superficie,

verificándose:

1
dS = dxdy
cos 

Y puesto que,

1
cos   n  k 
N

Resulta finalmente,

dS  N dxdy  fx(x , y )2  fy(x , y )2  1 dxdy

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
5

SUPERFICIE DADA POR LA ECUACIÓN CARTESIANA IMPLÍCITA: F (x , y, z )  0

Si la superficie viene dada por una función implícita F (x , y, z )  0 , con F de clase C 1 siendo
Fz  0 , un vector normal en cada punto es


N  Fx, Fy, Fz 
Con lo cual,

F , F , F  0, 0,1
x y z Fz
cos   n  k  
 Fy  Fz  Fy  Fz
2 2
F  F 
2 2 2 2
x x

Por tanto,

 Fy  Fz
2
F 
2 2
dxdy x
dS   dxdy
cos  Fz

SUPERFICIE DADA POR UNAS ECUACIONES PARAMÉTRICAS:

x  x (u, v ), y = y(u, v ), z = z (u, v ), con (u, v )  D

Conocidas las ecuaciones paramétricas, el vector de posición de un punto de la superficie es:

r(u, v )  x (u, v )i  y(u, v ) j  z (u, v )k, (u, v )  D

Suponiendo que las ecuaciones paramétricas son de clase C 1 en D , un vector normal unitario en
P0 es

i j k
ru  rv 1
n  x u yu z u
ru  rv ru  rv
x v yv z v

y el coseno director en P0 :

(x , y )
(u, v )
cos   n  k 
2 2 2
 (x , y )  (y, z )   (x , z ) 
       
 (u, v )   (u, v )  (u, v )

dxdy (x , y )
Sustituyendo la expresión anterior en dS  y teniendo en cuenta que dxdy  dudv
cos  (u, v )
, resulta
6 T3  INTEGRACIÓN DE SUPERFICIE

2 2 2
dxdy  (x , y )  (y, z )   (x , z ) 
dS           dudv  ru  rv dudv
cos   (u, v )   (u, v )  (u, v )
  


INTEGRAL DE SUPERFICIE DE UN CAMPO ESCALAR

5 Definiciones

La integral del campo escalar g(x , y, z ) sobre la superficie S se expresa con la siguiente notación:

 g(x, y, z )dS
S

donde,
• S es una superficie suave 1.
• g(x , y, z ) es un campo escalar continuo sobre S .

• dS , es la diferencial de superficie.
Una integral de superficie se calcula transformándola previamente en una integral doble, lo cual
requiere expresar el integrando en función de dos únicas variables independientes. El dominio de
integración de la integral doble pasará a ser la proyección de la superficie S sobre el plano de las
variables independientes elegidas.

Suponiendo S la superficie definida por z  f (x , y ) diferenciable, g(x , y, z ) continua sobre

S y utilizando la expresión obtenida para el diferencial de superficie, la integral de g(x , y, z )
campo escalar continuo sobre S , es

 g(x, y, z )dS =  g(x, y, f (x, y )) fx(x , y )2  fy(x , y )2  1 dA

S D
donde D es la proyección de S sobre el plano z  0 .

Suponiendo S la superficie definida por F x , y, z   0 , la integral de g(x , y, z ) campo

escalar continuo sobre S , es

 Fy  Fz
2
F 
2 2
x
 g(x, y, z )dS   g(x, y, z (x, y )) Fz
dxdy
S D
donde D es la proyección de S sobre el plano z  0 .

1
Una superficie suave es imagen de una función de clase C 1 definida en un dominio D .
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
7

Suponiendo S definida en paramétricas

x  x (u, v ), y = y(u, v ), z = z (u, v ), con (u, v )  D
la integral de g(x , y, z ) campo escalar continuo sobre S , es

 g(x, y, z )dS   g(x (u, v), y(u, v), z (u, v)) r   r  dudv
u v
S D

6 Propiedades

P1 (Linealidad).- La integral de una combinación lineal de funciones es la combinación lineal

de las integrales.

P2 (Aditividad sobre la superficie de integración).- La integral sobre una superficie que

sea unión de varias es la suma de las integrales sobre cada una de ellas; por ejemplo, para la
unión de dos superficies:
 g(x , y, z )dS =
S S
1 2
S
1
g(x , y, z )dS    g(x , y, z )dS
S
2

P3 (Independencia de la parametrización).- El valor de la integral no cambia con la

parametrización elegida para la superficie.

P4 (Independencia de la orientación).- El signo de la integral no cambia con la orientación

fijada en la superficie.

7 Algunas aplicaciones

Área de una superficie.- Si sobre una superficie S se integra el campo escalar constante
g(x , y, z ) = 1 , se obtendrá el área superficial de S .

área S    dS
S

Masa de una lámina.- Si S es la forma de una lámina, o elemento de dos dimensiones significativas
en el espacio, y el campo g da el valor de la densidad superficial en cada punto de la lámina, la masa
total de la lámina será la integral de g sobre S .

Temperatura media de una lámina.- Si S tiene la forma de una lámina, o elemento de dos
dimensiones significativas, y el campo g da el valor de la temperatura en cada punto de la lámina,
la integral de g sobre S , dividida por el área de S , es el valor de la temperatura media de la lámina.
8 T3  INTEGRACIÓN DE SUPERFICIE

INGEGRAL DE SUPERFICIE DE UN CAMPO VECTORIAL O INTEGRAL DE FLUJO

8 Definiciones

La integral del campo vectorial F sobre la superficie S es la integral de superficie del campo escalar
F ⋅n :   F  n dS
S

Los elementos que intervienen en la integral de un campo vectorial sobre una superficie son:
• La superficie orientada S , en las mismas condiciones de continuidad y derivabilidad
impuestas para definir la integral de un campo escalar sobre ella; llamaremos n al vector
normal unitario.
• El campo vectorial F(x , y, z ) definido y continuo sobre S .

• La diferencial de superficie dS .
Igual que la integral de un campo escalar, el cálculo de una integral de flujo se realiza previa
conversión en una integral doble, cuya expresión dependerá de las ecuaciones de la superficie.

SUPERFICIE DADA POR LA ECUACIÓN CARTESIANA EXPLÍCITA: z  f (x , y )

Si el campo es F = (M , N , P ) y la superficie viene dada por z = f (x , y ) en D , la integral de F

sobre S se puede expresar como

 F  ndS =  (Mf   Nf   P )dA

x y
S D

siendo n el vector normal que apunta en la dirección del eje OZ positivo.

SUPERFICIE DADA POR LA ECUACIÓN CARTESIANA IMPLÍCITA: F (x , y, z )  0

Suponiendo F de clase C 1 y Fz  0 , y tomando como vector normal el vector N  Fx, Fy, Fz ,  
la integral de flujo es

(Fx, Fy, Fz) Fx2  Fy2  Fz2

 F  ndS   M (x, y )i  N (x, y ) j  P(x, y )k dxdy 
S D Fx2  Fy2  Fz2 Fz
dxdy
  MF   NF   PF 
x y z
Fz
D

Nótese que al ser Fz  0 las variables independientes son x e y , lo cual implica proyectar la
superficie sobre el plano z  0 .

SUPERFICIE DADA POR UNAS ECUACIONES PARAMÉTRICAS:

x  x (u, v ), y  y(u, v ), z  z (u, v ), con (u, v )  D

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
9

En este caso la integral es,

 F  ndS   M (u, v)i  N (u, v) j  P(u, v)kr   r dudv

u v
S D

9 Aplicación al Cálculo de Flujos

Si la superficie S está sumergida en un fluido que tiene un campo de velocidades continuo V y si

S es el área de una pequeña parte de S , sobre la cual V se puede suponer constante, podemos
escribir
  V  n S

donde  es el volumen de fluido por unidad de tiempo, que cruza esa porción de área en dirección
y sentido de un vector normal unitario n. Sumando todos esos pequeños volúmenes obtendremos
que el flujo total de V a través de S es la integral

Flujo =   V  ndS
S

Esta misma expresión es válida para cualquier otro tipo de flujo, como puede ser el flujo eléctrico, el
flujo de calor, etc.

10 Teoremas

El teorema de la divergencia de Gauss relaciona la integral de flujo sobre una superficie cerrada con
la integral triple, sobre el volumen interior a dicha superficie, del campo escalar divergencia.

TEOREMA DE LA DIVERGENCIA DE GAUSS.-

• Hipótesis: los elementos que intervienen en este teorema son
- la superficie S , cerrada, suave por partes y orientada según la normal unitaria
exterior, n;
- la región H cuya frontera la superficie S ; se escribe H = S ;
- un campo vectorial F ( x, y, z ) de clase C 1 sobre H y ∂H ;
• Tesis: bajo estas hipótesis se verifica que
 F  ndS =  div FdV

H H

Ejemplo: En la imagen la superficie S es la

esfera, la región H es el sólido encerrado por S.

El flujo que atraviesa la esfera hacia fuera, ya

que se está considerando una normal unitaria
exterior, se puede calcular como una integral
triple de la divergencia del campo vectorial
sobre el sólido H.

El teorema de Stokes relaciona la integral de superficie que proporciona el flujo del rotacional de un
campo vectorial con la integral de línea de ese mismo campo.
10 T3  INTEGRACIÓN DE SUPERFICIE

TEOREMA DE STOKES.-
• Hipótesis: los elementos que intervienen en este teorema son
- la superficie S no cerrada, suave por partes, orientada según la normal
unitaria n;
- la curva “borde o frontera” de S , que denotamos , orientada conforme a la
orientación de S ;
- un campo vectorial F(x , y, z ) de clase C 1 sobre S y ∂S ;
• Tesis: bajo estas hipótesis se verifica que
 S
F  dr =   (rotF)  ndS
S

El convenio para determinar la orientación de la curva y el vector normal a la superficie coherente

con dicha orientación viene dado por la regla del sacacorchos: si se recorre ∂S , dejando la cara
exterior de la superficie a la izquierda, la normal asociada apunta desde cualquier punto de ella hacia
el exterior.

Ejemplo: Este resultado afirma que la

circulación de un campo vectorial a lo largo de
una curva cerrada es igual al flujo del rotacional
del campo a través de cualquier superficie S no
cerrada que tiene a C como frontera siempre
que se considere el sentido de la circulación y la
normal a la superficie relacionadas por la regla
del sacacorchos.

Es importante hacer notar que cada contorno C

tiene infinitas superficies a efectos de aplicar
este resultado.

Ejercicios propuestos

Solución: a)
1
En el resumen de teoría se recogen las dS  g 'x (x , z )2  g 'z (x , z )2  1 dxdz
expresiones del elemento diferencial de
superficie de una superficie dada por 1  2z 2
b) dS  dydz
z = f (x , y ) con (x , y )  R , estando R en el 1  y 2  2z 2
plano XY .
a. Escribe la expresión del elemento
diferencial de superficie adaptada a la 2
En el resumen de teoría se recoge la
superficie dada por y = g(x , z ) con
expresión del elemento diferencial se superficie
(x , z )  R estando R en el plano XZ . de una superficie dada por F (x , y, z ) = 0 si
b. Encuentra dS para la superficie
Fz  0 .
x = 1  y 2  2z 2 , tomando como
a. ¿Cómo se calcularía el elemento diferencial
variables independientes con
de superficie si Fz  0 ?
y, z   {(y, z ) / y 2
 2z 2 < 1}
b. Determina dS para la superficie
x 2  5y 2 = 25 y calcula el área de la
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
11

porción de ésta superficie situada en el

5
primer octante con 0  z  2 . Resuelve la (a) Calcula la integral de superficie
integral con Matlab.
  (x
2
I =  y 2 )dS
Solución: b) área = 11, 785 S

x2 y2 z2
donde S es la esfera   =0
3 a2 a2 b2
Determina el elemento diferencial de para 0  z  b
superficie para las siguientes superficies e indica (b) Calcula la integral de superficie
en cuáles el cociente entre el elemento
  (x  y 2 )dS donde S es la esfera
2
I =
diferencial de superficie y el elemento S

diferencial de área es constante. Prepara un x 2  y2  z 2 = a2

fichero con Octave/Matlab para dibujar las
z  x2
superficies de este ejercicio así como una (c) Calcula I = S
1  4x 2
dS donde S es la
muestra de sus vectores normales.
a) 2x  3y  5z  1 b) 2x  3y  1 porción de la superficie z  4  x 2 limitada por
c) z = 9  x 2  2y 2 d) z = 3 x 2  y 2 el plano x  2y  4
(d) Calcula la integral
e) z = x 2  3y 2 f) z = 4  x 2  2y 2
 4 
g) x 2  2y  4 I =   z  2x  y  dS donde S es la
S 3 
38 dS
Solución: a) dS  dxdy ,  cte ; x y z
5 dA porción del plano    1 situada en el
2 3 4
13 dS primer octante.
b) dS  dxdz ,  cte ;
3 dA
9  2y 2 Solución: (a) I = 8a 4 / 3 (b) I  8a 4 / 3
c) dS  dxdy ;
9  x 2  2y 2 (c) 12 (d) 4 61
dS
d) dS  10 dxdy ,  cte ;
dA 6
Encuentra la masa de la lámina S , su
2 2
x  6y densidad media si en cada punto la densidad de
e) dS  2 dxdy ;
x 2  3y 2 masa superficial de la placa es la función
g x , y, z  indicada.
f) dS  1  4x 2  16y 2 dxdy ;
a) S es la porción del plano z = x  2y que se
g) dS  1  x 2 dxdz proyecta en el rectángulo R   0,1   0, 3 del
   
plano XY ; g x , y, z   xyz
4 b) S es la parte de la superficie
Halla el área de la superficie de las
siguientes superficies: z = 4  x 2  y 2 interior a x 2  y 2  1 ;
a. paraboloide z = 9  x 2  y 2 acotada entre g x , y, z   z
los planos z = 1 y z = 5 . 7
Solución: a) Masa = 21 , densidad media=
z = x 2  y  1 comprendida entre los
2
b. 2
planos z=1 y z=4 
b) Masa =
60
 
7  55/2  41 , densidad media=
c. superficie esférica x 2  y 2  z 2  2 ,
situado por encima del plano z  1 . 7  55/2  41
 3, 44
d. cilindro x 2  z 2  a 2 de altura h. 
10 53/2  1 
e. esfera x 2  y 2  z 2  a 2
Solución: (a)   33  17  7
Una lámina S ocupa la parte del plano

b)

 17 7  5 5  d) 2ah e) 4a 2 
x  y  z  a en el primer octante ( a es un
número real positivo).
6
12 T3  INTEGRACIÓN DE SUPERFICIE

Si la temperatura en cada punto de S es Calcula el flujo de calor hacia arriba a través de


T x , y, z   k x  y  z2 2 2
 ( k  0 ), halla el la placa.
Solución: a) Flujo cka 3  0
valor medio de la temperatura en S .
 a 
Supón que T x , y, z  
1
; halla el valor b) Flujo 2c log a  1  
x y 1
 a  1
medio de la temperatura.
Solución: 12
ka 2 Se considera la superficie S porción del
a) Tmedia  paraboloide z = a 2  x 2  y 2 para x , y, z
2
2 positivos; a es un número real no nulo. Un
b) Tmedia
a

 2 a  log a  1  campo vectorial F depende de a de la
2
siguiente forma: F x , y, z   axi  yj  zk
a
a. Calcula a mano el flujo de F hacia arriba de
8
Se considera la superficie S dada por S
b. Analiza el comportamiento del flujo en
x = a  y 2  z 2 y el campo escalar
2
función del parámetro a .¿Es creciente el
 
g x , y, z   x y 2  z 2 . Halla el valor medio de flujo con el valor de a ? ¿Existe un valor
g en S mínimo para el flujo? Puedes dibujar con
3 Matlab la función flujo y utilizar el comando
a
Solución: gmedio  solve para hacer el análisis pedido.
4 a 3 2
Solución: Flujo=
8
 
a  a  2 b) El flujo
9 alcanza un valor máximo relativo para
Calcula el flujo saliente del campo a  1, 5662 y un valor mínimo relativo para
vectorial F = 2i  5 j  3k a través de la parte
a  0, 7662 .
del cono z = x 2  y 2 interior al cilindro
x 2  y2 = 1 . 13
Solución: Flujo  3 Considerando F x , y, z  = xi  yj  zk
, calcular el flujo que atraviesa la superficie del
10 sólido encerrado por las superficies
Calcula el flujo saliente del campo z  10  x 2  y 2 y z  2  x 2  y 2 .
vectorial F = 2i  5 j  3k a través de la parte Solución: Flujo = 48
del cono z = x 2  y 2 interior al cilindro
x 2  y2 = 1 . 14
Utiliza el teorema de la divergencia de
Gauss en los tres casos siguientes para calcular
Solución: Flujo = 3
I =  H
F  ndS :

11 a. F(x , y, z ) = zi  xj  yk ; H es el
Considera la superficie que ocupa la
parte del plano x  y  z  a en el primer hemisferio 0  z  9  x 2  y 2
octante ( a es un número real positivo) y el b. F(x , y, z ) = x 2 i  y 2 j  z 2 k ; H es el
campo de calor que es proporcional al gradiente sólido parabólico 0  z  4  x 2  y 2 .
de temperatura con signo negativo, esto es,
F  c T , c  0 para las funciones de c. F(x , y, z ) = (x  z 2 )i  (y  z 2 ) j  xk ; H

temperatura es el sólido 0  y 2  z 2  1 , 0  x  2 .
Solución: a) I = 0 b) I = 64 / 3 c) I = 4
a. 
T x , y, z   k x 2  y 2  z 2 
1
b. T x , y, z  
x y 1
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
13

usando la definición de integral de línea y el

15
Encuentra una expresión para la tercera Teorema de Stokes.
componente, P x , y, z  , para que el flujo del Solución: -9.

campo vectorial
F x , y, z  = xy  z  i  ze y j  P x , y, z  k
19
Calcula el trabajo realizado por el campo
hacia arriba del semielipsoide F(x , y, z ) = x 2y 3 i  j  zk a lo largo de la
z  9  2x 2  y 2 sea opuesto al flujo hacia circunferencia x 2  y 2 = R 2 en z = 0 .
debajo de la región D del plano XY Solución: Trabajo = R 6  / 8
encerrada en él y que éste último flujo sea el
doble del área de D .
z 2 y 20
Solución: P x , y, z   e  yz  2 Calcula la siguiente integral mediante la
2 fórmula de Stokes y después comprueba el
resultado integrando directamente
16 I =  xdx  (x  y)dy  (x  y  z )dz
Halla el flujo del rotacional de C

F x , y, z  = xy i  j  e k sobre la parte de
2 x2 siendo C la curva x = a sen t , y = a cos t ,
z  16  2x 2  y 2 en z  0 , orientada con la z = a(sen t  cos t ) , 0  t  2 .
normal hacia arriba. Solución: I = a 2
Solución: Flujo=0

21
Calcular la circulación de
17 F(x , y, z ) = (x  3z )i  (z  y ) j  4yk
Sea S la superficie de ecuación z=2
alrededor del trapecio de vértices
definida sobre el círculo x 2  y 2  4 y se A  (2, 2, 0) , B  (2, 2, 0) , C  (1,1,1) y
considera el campo D  (1, 1,1) tomando como orientación la
F x , y, z  = 3yi  xzj  yz 2 k . Calcular  F  dr que recorre los vértices en el orden en el que se
C
han definido.
usando la definición de integral de línea y el Solución: C=18
Teorema de Stokes.
Solución: 20
22
Comprobar el teorema de Stokes para
18 F x , y, z  = 2zi  xj  y 2 k considerando como
Sea C el triángulo orientado situado en
el plano 2x  2y  z  6 , intersecado por los superficie S el paraboloide z  4  x 2  y 2 y C
planos coordenados y es la curva de nivel de S sobre el plano XY.
F x , y, z  = y 2 i  zj  xk . Calcular Solución: 4
 F  dr
C

Test de autoevaluación

1
El área de la superficie de ecuación
B.

6

653  1 
z  16  x 2  y 2 situada sobre la región

R  x 2  y 2  16 es, 
C.

3

653  1 
D. Ninguna de las anteriores.
A.

6
 65  1 
14 T3  INTEGRACIÓN DE SUPERFICIE

a. Masa de una lámica

2 La temperatura media de la lámina S b. Trabajo

x 2  y 2  1, y  0 c. Circulación
definida por   y con d. Flujo

0z 3 e. Área superficial


temperatura en cada punto T (x , y, z )  y , es:
Las integrales y las magnitudes admiten los
2 siguientes emparejamientos:
A.
 A. I ,e, I , a , I , d , I ,b
1 2 3 5
6
B.

B. I , a , I ,e, I , d , I , c .
1 2 4 5

C. 0
D. Ninguna de las anteriores.
C. I , a , I ,e, I , c, I , d  .
1 2 3 4

D. Ninguna de las anteriores.

3
Sea  una superficie cuya ecuación en 6
coordenadas cilíndricas es z  f (r , ) , siendo El flujo del campo F  2xi  xzj  z 2 k
hacia el exterior de la esfera unitaria, es:
D la proyección de la superficie en el plano r
, entonces A. 0

2
 z   z 
2 B.
A. Área de    1  r 2      dr d  3
 r    
D
8
2 2
C.
 z  1  z  3
B. Área de    1       r dr d 
 r  r    D. Ninguna de las anteriores.
D

2 2
 z  1  z 
C. 1     2   r dr d  . 7
Área de     r  r    Sea  una superficie cerrada que
D

D. Ninguna de las anteriores. encierra un volumen V y F un campo vectorial

de clase C 2 , entonces

4
A.  div FdV  0
El área del casquete de la superficie V

esférica x 2  y 2  z 2  2 , situado por encima B.   F nd   0



del plano z  1 , es:
C.  div F nd   0



A. 2 2  2  1 D. Ninguna de las anteriores.

B. 4 2  2  1 .
8
 2  2  1 .
Sea V(x , y, z )  a  r , siendo a un
C.
vector constante y r el vector de posición de un
D. Ninguna de las anteriores. punto de  3 y sea C una curva de  3 que es
frontera de la superficie  . En estas
5 condiciones se verifica:
Se consideran las siguientes integrales:
I1   g(x, y, z )d  , I2   d  ,
A.  Vdr  2 a  n d  .
C 
 

I3   Fd σ , I4   divFdV ,
B.  Vnds   rot V n d  .
C 
 V

I5   FTds
C.  Vdr  2 rot r n d  .
C 
C

y las siguientes magnitudes: D. Ninguna de las anteriores.

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
15

9 Sea  el paraboloide de ecuación 10 El valor de la integral de flujo

z  x  y 2 con 0  z  3 y sea F el campo
2
 ( F)d σ , siendo S la semiesfera
vectorial F  xi  yj  zk . Entonces el flujo de S

rot F a través del paraboloide en el sentido de x 2  y 2  z 2  16, z  0 y F el campo

la normal exterior es:
A.  .
  
F  x 2  y  4 i  3xj  2xz  z 2 k , es: 
B.  . A. 32 .
C. 0 . B. 0 .
D. Ninguna de las anteriores. C. 2 .
D. Ninguna de las anteriores.

1 2 3 4 5 6 7 8 9 10
B A C A B C B A C A

Para realizar diferentes tests de autoevaluación de este tema de forma interactiva, visitar la
dirección web
https://www.giematic.unican.es/index.php/integral-de-superficie/material-interactivo
y pulsar los correspondientes enlaces en el apartado “Tests interactivos”.

Ejercicios resueltos
DIFERENCIAL DE SUPERFICIE. ÁREA DE UNA SUPERRFICIE

1
Utiliza la herramienta Área de un paralelogramo y área de su proyección.
Diferencial de superficie de la página
https://personales.unican.es/alvareze/CalculoWeb/CalculoII/integral_superficie.html
16 T3  INTEGRACIÓN DE SUPERFICIE

para observar cómo se obtiene la diferencial de superficie cuando la superficie viene

dada en explícitas.

2
Desde la página
http://personales.unican.es/alvareze/CalculoWeb/CalculoII/integral_superficie.html
selecciona la herramienta: Área de una superficie para ver cómo obtener el área de
una superficie que viene dada en explícitas.

3
Calcular el área del cono x 2  y 2  z 2  0 que es interior al cilindro x 2  y 2  ax  0

Solución
2 2
cono x +y -z 2 =0 dentro del cilindro x 2
+y 2
-2x=0

0
Z

-1

-2
2

1 2

1
0
0
-1
-1
Y -2 -2 X

Código Matlab
% Dibuja en paramétricas el cono x^2+y^2-z^2=0
r=0:.1:2; t=0:.05:2*pi;[R,T]=meshgrid (r,t);
h1=surf(R.*cos (T),R.*sin (T),R);
set(h1,'FaceColor','blue','FaceAlpha',0.2,'EdgeColor','c');
hold on
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
17

h2=surf (R.*cos (T),R.*sin (T),-R);

set(h2,'FaceColor','blue','FaceAlpha',0.2,'EdgeColor','c')
% Dibuja en paramétricas el cilindro x^2+y^2-ax=0
z=-2:.1:2; t=-pi/2:.05:pi/2;[Z,T]=meshgrid (z,t);
h3=surf(2*cos(T).*cos(T),2*cos(T).*sin(T),Z);
set(h3,'FaceColor','red','FaceAlpha',0.3,'EdgeColor','y')
xlabel('X');ylabel('Y');zlabel('Z');
title('cono x^2+y^2-z^2=0 dentro del cilindro x^2+y^2-2x=0')
hold off

Al tener el cilindro sus generatrices paralelas al eje OZ, resulta apropiado proyectar la superficie
del cono sobre el plano z  0 , utilizando para calcular el área pedida la expresión
A  dS .
D

Aplicado a nuestro caso, la ecuación de la superficie es z  f x , y    x 2  y 2 , lo que nos

lleva a:
2 2
   
 x    y   2
dS  1  f  f '2 '2
 1    2
 x 2  y 2  2 
 x  y 
x y
 

Con respecto al dominio D de integración, será el representado en la figura, intersección del

cilindro con el plano de proyección, es decir:
2 2
   
2 2
x  y  ax  0  x  a   y 2  a 
  2 
 2 

a  a
la circunferencia de centro  , 0 y radio es la curva que limita dicho dominio D.
 2  2

De esta forma, dada la simetría del cono respecto al plano z=0, ya que su ecuación no varía al
cambiar z por -z, el área pedida podrá expresarse

Área  2 2  dxdy
D

a
y puesto que  dxdy representa el área del dominio D, en este caso un círculo de radio ,
D
2
2
a  a 2
será  dxdy      , luego
D
 2  4

a 2 2a 2
Área= 2 2  u. d. s.
4 2

4
Encontrar el área de la región de la esfera x 2  y 2  z 2  9 limitada por el cilindro
x 2  y 2  3x  0

Solución
18 T3  INTEGRACIÓN DE SUPERFICIE

Se trata de calcular el área de la esfera acotda por el cilindro, habrá que calcular la integral de
superficie

área  2  dS = 2  fx(x , y )2  fy(x , y )2  1 dA

S D

donde S es la porción de esfera, de ecuación z = f ( x, y ) = 9 − x 2 − y 2 , cuya proyección sobre

el plano XY es el interior de la circunferencia x 2 + y 2 − 3 x =
0 que llamamos D.

3
Sustituyendo, área = 2  dxdy
D 9  x 2  y2

Describiendo el dominio D en coordenadas polares,

   
D  r ,   /     0  r  3 cos 
 2 2 

la integral será
3 cos 
 1/2 
 /2 3 cos 
1/2  3 9  r 2
 /2
  

área = 2   
3r 9  r 2  drd   2  
 2 1 
d 
 /2 0  /2  
 2  0

 /2
 18   sen   1d   18  36
 /2
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
19

5
Calcula el área de la lámina que tiene la forma de la superficie z  4  x 2  y 2
limitada por los planos z  0 , y  x , x  0 . Escribir la expresión que permite calcular la
temperatura en la lámina si dicha temperatura es en cada punto proporcional al cuadrado de
la distancia de dicho punto al eje OZ.
Solución
Considerando f x , y   4  x 2  y 2 , y llamando D a la proyección de S sobre el plano XY se
tiene que el área pedida es

   f 
2 2
área S    dS   1  fx' x
'
dxdy
S D

   f 
2 2
Como 1  fx' x
'
 1  4x 2  4y 2 , y D es un sector de la circunferencia de centro
(0,0) y radio 2, utilizamos coordenadas polares


2 2
   f 
2 2
área S    1  fx' x
'
dxdy   1  4r 2 r dr d 
D  0
4

r 2

1  4r 
3/2
2
  1
   
 2 4  8 3/2


48

17 17  1 
r 0

Para obtener la temperatura se deberá calcular

      f 
2 2
 T x, y, z  dS  k x 2  y2 1  fx' '
x
dxdy
S
S D


2 2 2
k
 k  r 2 1  4r 2 r dr d   r
3
1  4r 2 dr
 0
4 0
4
20 T3  INTEGRACIÓN DE SUPERFICIE

Las instrucciones Matlab para calcular esta integral serían:

>>syms k r
>>k*pi/4*int(r^3*sqrt(1+4*r^2),r,0,2)

Observación: Esta última integral se puede calcular por partes, aunque no se pide.

u  r 2  du  2rdr 

 
3/2 
1  4r 2 dr    
r  
3 2
 2 1 1  4r 
dv  r 1  4r dr  v  
 8 3/2 

 
5/2
2
r2 1 r2 1 1  4r
     
3/2 3/2 3/2
 1  4r 2  r 1  4r 2 dr  1  4r 2  
12 6 12 48 5/2

r2 1
   
3/2 5/2
 1  4r 2  1  4r 2
12 120
Con lo que el valor de la integral sería
r 2
k r 2 1 5/2 
   
3/2
I   2
1  4r 2 
4  12 1  4r 
120  
  r  0

k  17 3/2 1  k   391 17 1 

 
4  3

120
17 5/2  1  



4  120
 
120 

6
Determinar el área de la porción de superficie de la semiesfera
x  y  z 2  4a 2 z  0 que es interior al cilindro x 2  y 2  2ax  0
2 2

Solución

En este caso z  f x , y   4a 2  x 2  y 2  z 2 ya que se indica que z  0 . Si proyectamos

sobre el plano XOY y llamando D a la proyección, se tendrá, Área=  dS

D
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
21

2 2 2 2 2
semiesfera x +y +z =4 dentro del cilindro x +y -2x=0

1.5

Z 1

0.5

1.5
1
0.5
1
0
-0.5 0
-1
-1
-1.5
Y
X

Código Matlab:
% Dibuja en paramétricas la semiesfera x^2+y^2+z^2=4
f=0:.05:pi/2; t=0:.05:2*pi; [F,T]=meshgrid (f,t);
h1=surf (2*sin(F).*cos (T),2*sin(F).*sin (T),2*cos(F));
set(h1,'FaceColor','blue','FaceAlpha',0.2,'EdgeColor','c')
hold on
% Dibuja en paramétricas el cilindro x^2+y^2-2x=0
z=0:.1:2; t=-pi/2:.05:pi/2; [Z,T]=meshgrid (z,t);
h2=surf (2*cos(T).*cos (T),2*cos(T).*sin (T),Z);
set(h2,'FaceColor','green','FaceAlpha',0.3,'EdgeColor','y')
xlabel('X');ylabel('Y');zlabel('Z');
title('semiesfera x^2+y^2+z^2=4 dentro del cilindro x^2+y^2-2x=0')
axis equal
hold off

Teniendo en cuenta que el punto genérico x , y, z  está sobre la esfera, deberá verificar su

ecuación, es decir z   4a 2  (x 2  y 2 ) (recordemos z ≥ 0 ), luego,

2 2
   
 x    y
dS  1  fx'2  fy'2  1       .

 4a 2  (x 2  y 2 )   4a 2  (x 2  y 2 ) 

2a 2adxdy

4a 2  (x 2  y 2 )
 Área   4a 2  (x 2  y 2 )
D
22 T3  INTEGRACIÓN DE SUPERFICIE

D es el dominio, en el plano XOY, limitado por

la circunferencia de la figura, cuya ecuación se
puede expresar como

x 2  y 2  2ax  0 ,

o bien

x  a 
2
 y2  a2

Podemos usar simetrías, ya que al cambiar y por -y no cambia la ecuación de la frontera del
adxdy
dominio ni la función subintegral, luego podemos escribir Área= 4  , siendo
D1 4a 2  (x 2  y 2 )
D1 el semicírculo superior limitado por la circunferencia x 2  y 2  2ax  0 y el eje OX.

Hacemos el cambio a coordenadas polares:



x  r cos 


y  r sen 


 (x , y )
dxdy  drd   rdrd 
 (r , )

En coordenadas polares, la ecuación de la circunferencia que limita el dominio D’ se convierte

en r  2a cos  , luego

rdrd   /2 2a cos  rdr

Área= 4a   4a  d 
0 0
D 4a 2  r 2 4a 2  r 2
Y resolviendo queda

 
2a cos 
2a cos  r .dr  
 0 2
4a  r 2
  4a 2  r 2 
 0
 4a 2 1  1  cos2   2a(1  sen )

 /2  /2
Área= 8a 2  (1  sen )d   8a 2   cos    4a 2 (  2) u. d. s.
0  0

7
Determinar el área de la porción del cono x 2  y 2  3z 2 z  0 que es interior a la
4a 2
esfera x 2  y 2  z 2  .
3
Solución

El dominio D, proyectando sobre el plano XOY z  0 , lo obtendremos eliminando z entre las

dos superficies:
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
23

 2
x 2  y 2  3z 2  0
 z  1 x 2  y 2  4a 2 
 1
 2
x  y 2  z 2  4a
2 
 2
3
4a 2  

 x 2  y2  x 2  y2  
 z   x 2  y2  3 3
3  3

simplificando se tiene como ecuación de la frontera del dominio D la curva x 2  y 2  a 2 .

x 2  y2
Siendo en este caso z  f x , y   , ya que z  0 . El diferencial de superficie:
3
2 2
   
 1 x  1 y
dS  1  f  f '2 '2
 1       
x y  3 x 2  y 2   3 x 2  y 2 
 

x2 y2 4 2
 1   

3 x y 2 2
 
3 x y 2 2
 3 3

2
El área vendrá dada por la expresión A =  dS 
3
 dxdy , pero  dxdy representa el
S D D

área del dominio D, cuyo valor es πa 2 por ser un círculo, obteniendo:

2 2 3
A= a 2  a 2 u. d. s.
3 3

3.5

2.5

1.5

0.5

0
4 4
2 2
0 0
-2 -2
-4 -4

Código Matlab
% Dibuja en paramétricas el cono x^2+y^2=3z^2 con z>=0.
ra=sqrt(12); % radio de la esfera
ro=linspace(0,ra,20); % radio vector
t=linspace(0,2*pi,20); % ángulo theta
[RO,T]=meshgrid (ro,t);
h1=surf (RO.*sqrt(3)/2.*cos (T),RO.*sqrt(3)/2.*sin (T),RO./2);
set(h1,'FaceColor','blue','FaceAlpha',0.2,'EdgeColor','c')
hold on
% Dibuja en paramétricas el casquete de esfera x^2+y^2+z^2=12, con phi=<pi/3
f=linspace(0,pi/3,20); % ángulo phi
t=linspace(0,2*pi,20); % ángulo theta
[F,T]=meshgrid (f,t);
h2=surf (ra*sin(F).*cos (T),ra*sin(F).*sin (T),ra.*cos(F));
24 T3  INTEGRACIÓN DE SUPERFICIE

set(h2,'FaceColor','green','FaceAlpha',0.3,'EdgeColor','y')
xlabel('X');ylabel('Y');zlabel('Z');
% Dibuja la circunferencia x^2+y^2=a^2
t=0:.05:2*pi;
h3=plot(3*cos(t),3*sin(t),'b');
set(h3,'LineWidth',3);
title('cono x^2+y^2=3z^2 dentro de la semiesfera x^2+y^2+z^2=12')
axis equal; hold off

8
Determinar el área de la porción del cilindro x 2  y 2  6y  0 que es interior a la esfera
x 2  y 2  z 2  36 .

Nota: se aconseja proyectar sobre el plano x = 0 (plano YOZ).

Solución

En este caso, x  f y, z    6y  y 2 por lo que habrá que proyectar sobre el plano YOZ,

2
 
 2y  6 
dS  1  fy'2  fz'2  1    0 
 2 6y  y 2 
 

2y  6
2
24y  4y 2  4y 2  36  24y 36 3
 1   

4 6y  y 2
 
4 6y  y 2
 
4 6y  y 2
 6y  y 2

despejando,

1 4x 2  4y 2  24y  36 (x 2  y 2  6y )  9
 
cos  2x x
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
25

Teniendo en cuenta que el área total será dos veces el área en x≥0, se tendrá

3 dydz dydz
Área   dS  2 6y  y 2
6 
6y  y 2
S D D

siendo D el dominio limitado por la curva resultante

de eliminar x entre el cilindro y la esfera, es decir,

 2 2
x  y  6y  0
 2  z 2  6y  36

x  y 2
 z 2
 36


que es una parábola en el plano YOZ;

podemos usar simetrías, ya que al cambiar z por -z no cambia la ecuación de la frontera del
dominio ni la función subintegral, luego podemos escribir

dy.dz 6 dy  366y 
Área= 12   12   dz 
D 6y  y 2 0
6y  y 2  0 

siendo D’ la región del primer cuadrante encerrada por la parábola en el cuadrante de las z
positivas. Resolvemos la integral entre corchetes,
366y 366y
 dz  z   36  6y
0 0

se sustituye en la expresión anterior

36  6y .dy 6 6  ydy dy 6
 12 6 2 y   144 u.d.s.
6 6 6
Área  12   12   12 6 
0
6y  y 2 0
y 6 y 0
y  0

INTEGRAL DE UN CAMPO VECTORIAL SOBRE UNA SUPERFICIE. FLUJO

9
Desde la página
http://personales.unican.es/alvareze/CalculoWeb/CalculoII/integral_superficie.html
utilizar la herramienta Integral de superficie de un campo vectorial para ver una
justificación de la definición de esta integral.
26 T3  INTEGRACIÓN DE SUPERFICIE

10
Se considera el casquete esférico S , definido por x 2  y 2  z 2  2, z  1 , con la
orientación normal exterior a la esfera. Calcular el flujo del campo F = (1, 0,1) .

Solución

=
La integral pedida es I ∫∫ (1, 0,1) ⋅ n dS donde,
S
ndS = NdA , siendo

 
 x x
N  z x , z y ,1  
   , ,1
 2  x 2  y 2 
 2  x  y
2 2

por lo tanto,
 
 x
 1, 0,1  n dS   1  dxdy
S D  2  x 2  y 2 

Resolvemos esta integral en polares, teniendo en cuenta que el dominio de integración es el

interior del círculo que se obtiene cortando la esfera por el plano z = 1


D  (x , y ) / x 2  y 2  1 
Es decir,
2 1  
1  r cos   rdrd 
I   


2  r 2 
0 0

2

Si integramos primero respecto de θ y tenemos en cuenta que  cos d   0 , queda

0
2 1

únicamente la otra integral I    rdrd    .

0 0
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
27

11
Calcular El flujo del campo F  x , y, 2z  a través de la porción de cilindro x 2  y 2  1
, comprendido entre los planos z  1 y z  2 , en el sentido de la normal exterior al cilindro.
Solución
Calcularemos la integral de superficie en paramétricas, utilizando las siguientes ecuaciones para
el cilindro
x  cos , y  sen , z  z

El integrando en paramétricas es

N
FndS  F   N dA  F  NdA
N

donde,

i j k
N  r  rz   sen  cos  0  cos , sen , 0 y dA  d dz
0 0 1

por lo tanto,

F  N  F  cos , sen , 2z   cos , sen , 0  1

Y la integral de flujo es
2 2

   dzd   2
0 1

12
Hallar el flujo del campo E(x , y, z )  xi  yj  zk , que sale de la superficie cerrada S, que
constituye la frontera de un cubo al que se le ha quitado una esquina, como muestra la figura.
Solución

El campo F(x , y,.z )  xi  yj  zk es de clase C1 en  3 . S

es una superficie cerrada, suave a trozos, que encierra un
sólido V simplemente conexo.

Flujo saliente de V =  F  n dS
S

Para calcular el flujo total es necesario hallar el flujo saliente a través de cada una de las nueve
caras del sólido.

S1: base inferior, z = 0, n1  k ; 1   F  n  zdxdy  0dxdy 0

1
S1 S1 S1
28 T3  INTEGRACIÓN DE SUPERFICIE

S2: cara z = 2, n 2  k ; 2   F  n  zdxdy  2dxdy 2área S

2 2
6
S2 S2 S2

S3: cara y = -1, n 3   j ; 3   F  n 3

 ydxdz  dxdz área S 3  4
S3 S3 S3

S4: cara y = 1, n 4  j ; 4   F  n 3
 ydxdz  dxdz área S 4  3
S4 S4 S4

S5: cara x = -1, n 5  i ; 5   F  n 5

 xdydz  dydz área S 5  4
S5 S5 S5

S6: cara x = 1, n 6  i ; 6   F  n 6
 xdydz  dydz área S 6  3
S6 S6 S6

S7: cara z = 1, n 7  k ; 7   F  n 7
 zdxdy  dxdy área S 7  1
S7 S7 S7

S8: cara y = 0, n 8  j ; 8   F  n 8
 ydxdz  0dxdz 0
S8 S8 S8

S9: cara x = 0, n 9  i ; 9   F  n 9
 xdydz  0dydz 0
S9 S9 S9

Flujo total saliente de V = 1  2  3  4  5  6  7  8  9  21

Utilizando el teorema de Gauss:

Flujo saliente de V =  F  n dS   (div F)dxdydz

S V

x y z
Como div F     3 , sustituyendo queda
x y z

Flujo saliente de V =  3dxdydz = 3 volumen de V

V

Flujo saliente de V = 3 (vol. cubo lado 2 – vol. cubo lado 1) = 3(8-1) = 21 unidades de flujo

13
Se considera la porción del paraboloide z = 1  x 2  y 2 , con x , y, z positivos y un campo
vectorial, V x , y, z   xi  2yj  zk . Plantear la integral para calcular el flujo hacia el exterior
del paraboloide.
Solución
El flujo es

Flujo   V  n dS   V  N dA
S D

donde,

• El campo actuando sobre la superficie es: V  (x ,  2y, 1  x 2  y 2 )

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
29

• N  z x , z y ,1  2x , 2y,1

• dA  dxdy

Efectuando el producto escalar V  n , agrupando términos y sustituyendo en la integral de flujo

se obtiene esta integral doble:

 x 
2
Flujo   5y 2  1 dA
Dxy

siendo Dxy el cuadrante de círculo de radio 1 situado en el primer cuadrante del plano XY .

Si planteamos la integral, utilizando coordenadas polares, se tiene

 /2 1

  
x 2  5y 2  1 dA    r 
2
Flujo  cos2   r 2 sen2   1 rdrd 
Dxy 0 0

TEOREMA DE LA DIVERGENCIA O DE GAUSS

14
Calcular el flujo del campo F x , y, z   y cos y i  sen x j  z k hacia el exterior de la
superficie S que es frontera del sólido H comprendido entre las superficies z  1  y ,
x 2  y2  1 y z  0 .
Solución
S es la frontera del sólido H interior al cilindro acotado inferiormente por z  0 y superiormente
por el plano z  1  y .

Como se cumplen las hipótesis para aplicar el Teorema de la divergencia, se tendrá que el flujo
se puede calcular como la integral triple sobre H de la divergencia de F:

 FndS   div FdV

S H

Teniendo en cuenta que:

div F  1
30 T3  INTEGRACIÓN DE SUPERFICIE

H  x, y, z / x 2
 y 2  1, 0  z  1  y 
 r, , z / 0  r  1, 0  z  1  r sen  
el flujo será:
2  1 1r sen  2 1

  r dz dr d     r 1  r sen  drd  
0 0 0 0 0

r 1
2  2
 r r3 
2
1 1    2
  sen  d    1   1 cos  
  2


3
sen   d  
   2 3 
  2 3

  0

0 r 0 0

15
Halla el flujo del campo vectorial F(x , y, z )  x , y, z  a través de:

a) La superficie lateral S1 del paraboloide x 2  y 2  2az con 0  z  2a .

b) La superficie de la tapa S 2 para z  2a .

S1  S 2
c) Aplica el teorema de Gauss a la superficie comprobando previamente que se
cumplen las hipótesis de dicho resultado.
Solución a)

x y 
x y 2 2  
N  fx' , fy' , 1   , .  1
a a 
El vector normal a la superficie z  f x , y   es . El
2a
conjunto D proyección de S1 sobre XY es el círculo de centro (0,0) y radio 2a .

 x 2  y 2   x y 
flujoS   F  n dS   F  N dA   x, y,
    , , 1 dA 
1
2a  a a 
S1 D D

 4
 2a  
2  2a
2x 2  2y 2  x 2  y 2 x 2  y2 r3
  dA   dA   drd      4a

3

2a 2a 2a 
a  4 
D D 0 0
 

Solución b)

El vector normal es N  k . El conjunto D es la proyección de S 2 sobre el plano XY es el círculo

de centro (0,0) y radio 2a .
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
31

flujoS 
2
 F  n dS   F  N dA   x, y, 2a   0, 0,1dA 
S2 D D

 2a dA  2aÁrea D   8a 
3

D

Apartado c)

S1  S 2
Para calcular el flujo a través de aplicando el Teorema de Gauss se tendrá que
2  2a 2a

flujoS S 
1 2
 divFdV   3dV     3rdzdrd  
H H 0 0 r 2 /2a

 4

2ar  r  dr  6 2a 2a   2a    6 4a 3  2a 3  12a 3

2
2a  3  
 3  2  


2a 

 2 8a 
 
0  

16
Calcular, utilizando y sin utilizar el teorema de Gauss, el flujo del campo

F x , y, z   x 3i + y3 j  z 3 k

a través de la superficie S que es el exterior del tetraedro formado por los planos x  y  z  a
, x  0 , y  0, z  0.

Nota: Debes enunciar el teorema de Gauss y justificar si se puede aplicar en el caso del ejercicio.
Solución

La superficie S está dada por las cuatro caras del tetraedro de la figura.

Si se aplica el teorema de Gauss

 FdS   div F dxdydz   3x 

2
I   3y 2  3z 2 dxdydz
S V V

La proyección del sólido V sobre el plano XY es el triángulo limitado por los ejes coordenados y la
recta x + y = a , por lo tanto
32 T3  INTEGRACIÓN DE SUPERFICIE

a a x a x y
I     3x 2

 3y 2  3z 2 dzdydx
0 0 0

Integrando
a a x
z a x y
  3x z  3y z  z 
2 2 3
I  dydx 
z 0
0 0

a a x
 2 3
  3x a  x  y   3y a  x  y   a  x  y   dydx
2

0 0

y a x
 2 4
a  3x 2 a  x  y  3y 4 a  x  y  

 y 3 a  x   
    

dx 
 2 4 4 
0  y  0

 4 2 4  4 2
a
 4 3 a  x  3x 2 a  x  a  x   a
a  x  3x 2 a  x  
  a  x      dx      dx

 4 2 4   2 2
0   0  

x a
a 
a  x 

4

3x 2 a 2  2xa  x 2 dx   a  x  5
x a 3 2
3x a 3x 4 
5

   2

2  
 10

2

4

10 

0    
x 0

a 5 3a 5 3a 5 a 5 3a 5
    
2 4 10 10 20

Para resolver la integral directamente, sin aplicar el teorema de Gauss, hay que calcular el flujo a
través de cada una de las cuatro caras del tetraedro:
Flujo a través de S1 , que se define por la ecuación z=0,

 FdS   x 
, y 3 , 0 0, 0, 1 dA  0
3

S1 D1

Flujo a través de S 2 , que se define por la ecuación x=0,

 FdS   0, y 
, z3 1, 0, 0 dA  0
3

S2 D2

Flujo a través de S3 , que se define por la ecuación y=0,

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
33

 FdS   x 
, 0, z 3 0, 1, 0 dA  0
3

S3 D3

Flujo a través de S 4 , que se define por la ecuación z  a  x  y , si la proyección de S 4 es D4

a a x
 3 3 3  3 3
 FdS   x , y , a  x  y  1,1,1 dA    x  y 3  a  x  y   dydx 
 
S4 D4 0 0

y a x
 4  4 4
a
 4
x 3y  y 
a  x  y   a
 3 a  x  a  x  
  
dx   x a  x    dx 
 4 4   4 4 
0
 y 0
0  

a
 4  4 5
a 
x 3a  x 4 
a  x  
 x a x

5 a  x   a 5 a 5 a 5 3a 5
  2
dx


 4

5

10 

4

5

10

20
0   
 0

3a 5
Por lo tanto, el flujo que se pedía es la suma de los valores obtenidos para cada superficie .
20

17
Calcular, utilizando y sin utilizar el teorema de Gauss, el flujo del campo vectorial

F x , y, z   xz ,  y 2 , xz  a través de la superficie cerrada limitada por el cilindro

x 2  y 2  R2 con 0  z  3 .

Solución (sin utlizar el Teorema de Gauss)

La superficie cerrada S que limita el cilindro es la unión de tres superficies, la tapa superior S1, la
inferior S2 y la superficie cilíndrica S3.

Parametrización S1

S1  r1 D1  D1  x, y    2
/ x 2  y 2  R2 
r1 x , y   x , y, 3
34 T3  INTEGRACIÓN DE SUPERFICIE

El vector normal exterior de la superficie S es

N 1  0, 0,1

 F  n dS   F r1 x, y   0, 0,1dxdy   3x, y 

, 3x  0, 0,1dxdy 
2

S1 D1 D1

2  R
 
    3r cos d dr  0
2
 3xdxdy 
 0 
D1 0 

Parametrización S2

S2  r2 D2  D2  x, y    2
/ x 2  y 2  R2 
r2 x , y   x , y, 0

El vector normal exterior de la superficie S es

N 2  0, 0, 1

 F  n dS   F r2 x, y   0, 0, 1dxdy   0, y 

, 0  0, 0, 1dxdy  0
2

S2 D2 D2

Parametrizamos S3

S 3  r3 D3  D3  u, v    2
/ 0  u  2, 0  v  3 
r3 u, v   R cos u, R sen u, v 

El vector normal exterior de la superficie S es

N 3  Tu xTv  R cos u, R sen u, 0

El flujo de F a través de S3 es

 F  n dS   F r3 u, v   R cos u, R sen u, 0dudv 

S3 D3

 Rv cos u, R 
sen2 u, Rv cos u  R cos u, R sen u, 0dudv 
2

D3

3 2
9
  R v cos 
2 2
u  R 3 sen 3 u dudv  R2
0 0
2

Solución (Utilizando el teorema de Gauss)

Se tendrá

 F  n dS   divF dxdydz

S V

siendo V el volumen limitado por la superficie S. Se tiene que

div F  z  2y  x

Para calcular la integral triple se pasa a coordenadas cilíndricas

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
35

I   F  n dS   divF dxdydz   z  2y  x dxdydz

S V V

siendo V  r, , z  / 0  r  R, 0    2, 0  z  3 . Se tendrá entonces 

R 2 3
I     z  2r sen   r cos  rdzddr 
0 0 0

R 2 R
9  9
 r  6r 2 sen   3r 2 cos  ddr 
  9rdr  2 R
2
  2 
0 0
 0

18
Dado el campo vectorial F(x , y, z )  2xi  yj  z 2 k , se pide determinar el flujo a través
de la cara lateral del cilindro x 2  y 2  a 2 , 0z 2.

Solución
Código Matlab:
% Dibuja una muestra del campo vectorial
% F=(2x,y,z^2)
x=-2:2;y=-2:2;z=0:2;
[X,Y,Z]=meshgrid (x,y,z);
U=2*X;V=Y;W=Z.^2;
quiver3(X,Y,Z,U,V,W)
hold on
% Dibuja en paramétricas el cilindro x^2+y^2=4
z=0:.05:2; t=linspace(0,2*pi,30); [Z,T]=meshgrid (z,t);
h1=surf (2*cos (T),2*sin (T),Z);
set(h1,'FaceColor','b','FaceAlpha',0.2,'EdgeColor','b')
xlabel('X');ylabel('Y');zlabel('Z');
title('vector F=(2x,y,z^2) a través del cilindro x^2+y^2=4')
axis equal; hold off

2 2 2
vector F=(2x,y,z ) a través del cilindro x +y =4

1
Z

2
2
1
0
0
-1
-2 -2
Y X

Se cumplen las condiciones para aplicar el Teorema de la divergencia:

 div F dV   F  n dS
V S
36 T3  INTEGRACIÓN DE SUPERFICIE

donde n es el vector normal unitario saliente de la superficie S. Es decir, la integral de la

divergencia de un campo vectorial F x , y, z  extendida al volumen V, es igual al flujo saliente del
vector F a través de la superficie que limita este volumen.
En este caso se cumple que

sup  inf  lat   div FdV   div Fdxdydz ,

V V

 (2x )  (y )  (z 2 )
div F     2  1  2z  3  2z
x y z
luego,

 div FdV   (3  2z )dxdydz .

V V

Para calcular esta integral triple calculamos el cambio a coordenadas cilíndricas

x  r cos 

y  r sen  ,  (x , y, z )
 r
  (r , , z )
z  z


las ecuaciones del dominio x 2  y2  a2 , 0z 2, se transforman en

r 2  a2 , 0  z  2 , luego los nuevos límites de integración serán:

D  (r , , z ) / 0  r  a, 0    2, 0  z  2

por tanto,
2 a 2  2    q  
2

 
 (3  2z )dxdydz   r 3  2z dz dr d     d     rdr    3  2z dz 
V 0 0 0

0
  0   0
 

Podemos resolver las tres integrales por separado, ya que son independientes:

2 2 a  r 2 a a2 2
(3  2z )dz  3z  z 2   10 ,
2
0 0 rdr     y o d      2
 0  2  0 2 0

luego

 (3  2z )dx dy dz  10a

2
unidades de flujo.
V

Por otro lado, calcularemos el flujo a través de la cara superior φ sup . En esta cara z  2 y
dS  dx dy , n  k  F  n  z2  4

sup   F  ndS  4  dxdy  4Área de S sup

 4a 2 unidades de flujo
S sup sup

Análogamente para la cara inferior φ inf , en ella z 0 y dS  dxdy ,

n  k  F  n  z 2  0 , luego el flujo saliente a través de la cara inferior es nulo.
Podemos escribir:
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
37

lat   div F dxdydz   sup

 inf  10a 2  4a 2  0  6a 2 unidades de flujo.
V

19
Determinar el flujo saliente originado por el campo de fuerzas F x , y, z  a través de la

cara lateral del cilindro del sólido 

V  0  x 2  y 2  9, 0  z  1 ,  siendo
F x , y, z   xi  2y 2 j  3z 2 k .

Solución
2 2 2 2
vector F=(x,2y ,3z ) a través del cilindro x +y =9

1
0.5
Z

2
0
0

-2 -2
Y
X

Código Matlab:
% Dibuja una muestra del campo vectorial
% F=(x,2y^2,3z^2)
x=-3:3;y=-3:3;z=0:1;[X,Y,Z]=meshgrid (x,y,z);
U=X;V=2*Y.^2;W=3*Z.^2;
quiver3(X,Y,Z,U,V,W); hold on
% Dibuja en paramétricas el cilindro x^2+y^2=9
z=0:.05:1; t=linspace(0,2*pi,30); [Z,T]=meshgrid (z,t);
h1=surf (2*cos (T),2*sin (T),Z);
set(h1,'FaceColor','b','FaceAlpha',0.2,'EdgeColor','b')
xlabel('X');ylabel('Y');zlabel('Z');
title('vector F=(x,2y^2,3z^2) a través del cilindro x^2+y^2=9')
axis equal; hold off

En nuestro caso, F x , y, z   xi  2y 2 j  3z 2 k es de clase C1 en  3 y V es un sólido

simplemente conexo. Por tanto podemos aplicar el teorema de la divergencia o teorema de
Gauss:

 div F dV   (F  n)dS

V S

 (x )  (2y 2 )  (3z 2 )
Se tiene que: div F     1  4y  6z
x y z
El flujo saliente total de V se obtendrá calculado el flujo a lo largo de la cara lateral y las caras
superior e inferior:

total  lateral  sup erior  inf erior

38 T3  INTEGRACIÓN DE SUPERFICIE

Aplicando el Teorema de Gauss se tendrá que

 1 
total   (1  4y  6z )dxdydz   dxdy   0
(1  4y  6z )dz 

V D

Como
1 1
 (1  4y  6z )dz  z  4yz  3z 2   4 1  y 
0  0
se tendrá

total  4  (1  y )dxdy
D

siendo el dominio en el plano D  0  x 2  y 2  9 . Para realizar esta integral doble se pasa a

coordenadas polares,
 
x  r cos   J r
y  r sen 

siendo el dominio D los puntos r ,   con 0  r  3 y 0    2 . Por lo tanto,

2  3 
total  4  (1  rsen )rdrd   4  d    (r  r 2 sen )dr 
D
0  0 

3
3 r 2 r 3  9
 (r  r sen )dr    sen     9 sen 
2

 2 3  0 2
0

2
9     
total  4 
2
  9sen   d   36 2   1  sen   d   36    cos    36
0  2 
 0  2 
 2 
 0

Por otro lado en la cara superior tenemos n  k  F  n  3z 2 y como es z  1 queda

F  n  3, dS  dxdy , por tanto

superior   (F  n)dS  3 dxdy  3Área de S  3 32  27

S S

Análogamente para la cara inferior: n =−k ⇒ F ⋅ n =−3z 2 y como es z  0 queda

0.
F  n  0, dS  dxdy , por tanto φinferior =

Despejando el flujo pedido se obtendrá

lateral  total  superior  inferior  36  27   9 unidades de flujo

TEOREMA DE STOKES
20
El cilindro x 2  y 2  1 corta al plano y  z  2 en la curva C. Calcular la circulación del
campo vectorial F  xy i  yz j  xz k sobre C, utilizando obligatoriamente el teorema de
Stokes y recorriendo la curva C en sentido antihorario (vista desde arriba).
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
39

Nota: Se deberá comprobar que se cumplen las hipótesis del Teorema de Stokes
Solución

Aplicando el Teorema de Stokes, la circulación del campo F a lo largo de la curva C será

circulación   F  dr   rotF  ndS

C S

donde S es la superficie interior a la curva C que se encuentra sobre el plano y  z  2 .

La integral de superficie se calcula de la forma siguiente:

 rotF  ndS   rotF  NdA

S D

donde
• D es la proyección de S sobre el plano XY, es decir, el círculo unidad
• Como S es la porción del plano z  f x , y   2  y limitada por C, el vector normal a la

 
superficie es: N  fx' , fy' ,1  0,1,1

i j k
  
• rotF   y i - z j  x k
x y z
xy yz xz

Por lo tanto,

circulación   y, 2  y, x   0,1,1dA   2  y  x dA 

D D

2 1 2
 1 1 
   2  r sen   r cos   rdrd    1  3 sen   3 cos d   2
 
0 0 0

Aunque el ejercicio no pide obtener directamente la integral de línea podemos calcularla para
verificar el resultado:
40 T3  INTEGRACIÓN DE SUPERFICIE

circulación   F  dr   C
xydx  yzdy  xzdz
C

siendo C la curva

x t   cos t y t   sen t z t   2  sen t t  0, 2 

x ' t    sen t y ' t   cos t z ' t    cos t

2

  cos t sen t  sen t 2  sen t  cos t  cos t 2  sen t  cos t dt  2 

2
circulación 
0

21
Verifica el Teorema de Stokes para el campo vectorial

F ( x, y , z ) = 3 y i + 4 z j − 6 x k

y la parte de la superficie del paraboloide z =9 − x 2 − y 2 situada sobre el plano XY orientada hacia

arriba. Comprueba además que se cumplen las condiciones para poder aplicar el Teorema de
Stokes.
Solucion

La superficie S es el paraboloide de ecuación z  f x , y   9  x 2  y 2 , superficie suave cuya

curva frontera C es la circunferencia de centro (0,0) y radio 3. Para considerar la normal unitaria
exterior se debe orientar C en sentido antihorario. La parametrización de esta curva C será:

x  3 cos t



C  y  3 sen t t  0, 2 



z 0


Calculamos S
F  dr

2 2

  9 sen t, 0, 18 cos t   3 sen t, 3 cos t, 0dt   27 sen

2
F  dr = tdt 
S
0 0
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
41

2
1  cos 2t
 27  dt  27 
0
2

Calculamos   (rotF)  ndS

S

i j k
  
Como rotF   4i  6j  3k , n  f 'x , f 'y ,1  2x , 2y,1
x y z
3y 4z 6x

se tiene

  (rotF)  ndS    4, 6, 3  2x, 2y,1dxdy

S D

Siendo D la proyección del paraboloide sobre el plano XY, es decir, el interior de la circunferencia
de centro (0,0) y radio 3. Pasando a coordenadas polares
2 3

  4, 6, 3  2x, 2y,1dxdy    8r cos   12r sen   3 rdrd   27
D
0 0

22
 y dx  zdy  xdz , siendo C el triángulo situado en el
2
Hallar el valor de la integral
C
plano 2x  2y  z  6 intersecado por los planos coordenados.

Solución
Aplicando el Teorema de Stokes y llamando
- S al triángulo encerrado por C
- F x , y, z   y 2 i  z j  x k

se tiene que:

 y dx  zdy  xdz   rot F  n dS

2

C S

Como

i j k
  
rot F x , y, z     i - j  2y j
x y z
y 2 z x

N  2, 2,1

se tiene

 y dx  zdy  xdz   rot F  n dS   4  2y dA

2
I 
C S T
42 T3  INTEGRACIÓN DE SUPERFICIE

3 3x 3 y 3x 3 y 3x

  4  2y dydx    y 4  y 
 
2
I  4y  y  dx   
dx 
0 0 0 y 0 0 y 0

x 3
3  x3 
  3  2x  x dx  3x  x 2    9
2

0
 3 x 0

23
Verificar que se cumple el Teorema de Stokes considerando


F x , y, z   2z , x , y 2 
y S la superficie z  4  x 2  y 2 del primer octante ( z  0, x  0, y  0 ).

Solución
La superficie S es una porción de paraboloide z  4  x 2  y 2

Se verifican las condiciones del Teorema de Stokes así que el flujo del rotacional del campo se
puede calcular como una integral de línea sobre la frontera de la superficie

 rot F  n dS   F  Tds
S S

Cálculo de la integral de superficie

rot F x , y, z   2y, 2,1

Considerando el normal exterior y llamando D a la proyección de S sobre z=0 se tendrá:

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
43

 rot F  n dS =  2y, 2,1  2x, 2y,1dA

S D

2 4x 2 2  /2

  4xy  4y  1dydx    4r 

2
cos  sen   4r sen   1 rd dr
0 0 0 0

r 2
 /2
 4r 3 r2 
 /2
 32  56
 r 4 cos  sen   sen    d    16 cos  sen   sen   2 d    
0
 3 2 r 0 0
3  3

Calculo de la integral de línea

Como se puede comprobar la frontera de la superficie es una curva cerrada suave por partes que
se puede escribir como

 F  Tds   F  Tds   F  Tds   F  Tds

S C1 C2 C3

siendo C 1 la frontera de la superficie en el plazo XZ, C 2 la frontera de la superficie en el plazo

z=0 y C 3 la frontera de la superficie en el plazo YZ.

Curva frontera C 1 (plano XZ)


r t   t, 0, 4  t 2  t  0, 2 r' t   1, 0, 2t 

 
F  r  8  2t 2 , t, 0  1, 0, 2t   8  2t 2

2
32
 F  Tds   8  2t dt 
2

C1 0
3

Curva frontera C 2 (plano z=0)

 

r t   2 cos t, 2 sen t, 0 t  0, 
 2  
r' t   2 sen t, 2 cos t, 0 
 

 
F  r  0, 2 cos t, 4 sen2 t  2 sen t , 2 cos t , 0  4 cos2 t  
 /2

 F  Tds   cos t dt  

2

C2 0

Curva frontera C 3 (plano YZ)


r t   0, t, 4  t2  t  0, 2 r' t   0,1, 2 t


F  r  8  2 t2 , 0, t 2  0,1, 2 t  2t 3 
2

 F  Tds   2t 3dt  8

C3 0

Luego,
44 T3  INTEGRACIÓN DE SUPERFICIE

32  56
 F  Tds  3
 8 
4 3

C

24
Comprobar que se verifica el Teorema de Stokes para el campo vectorial
E(x , y, z )  3yi  3xj  z 4 k sobre la superficie S definida como el trozo del elipsoide
1
2x 2  2y 2  z 2  1 que está por encima del plano z  . Para comprobar el teorema es
2
necesario calcular las integrales que aparecen en ambos miembros de la ecuación y obtener el
mismo resultado por ambos métodos de cálculo.
Solución
2 2 2
Flujo de rot(F)=6k, a través del elipsoide 2x +2y +z =1,con z > = 1/sqrt(2)

1
Z

0.5

0
0.5 0.5

0 0

-0.5 -0.5
Y X

Código Matlab:
% Dibuja una muestra del campo vectorial rotF=(0,0,6)
x=linspace(-2/3,1/3,4);y=x;z=linspace(0,1,3);
[X,Y,Z]=meshgrid (x,y,z);
U=0*ones(size(X));V=0*ones(size(X));W=6*ones(size(X));
quiver3(X,Y,Z,U,V,W); hold on
% Dibuja el elipsoide 2x^2+2y^2+z^2=1, para z>=1/sqrt(2),en
paramétricas
f=linspace(0,pi/4,20); % ángulo phi
t=linspace(0,2*pi,30); % ángulo theta
[F,T]=meshgrid(f,t); a=1/sqrt(2);
h1=surf(a*sin(F).*cos(T),a*sin(F).*sin(T),cos(F));
set(h1,'FaceColor','red','FaceAlpha',0.2,'EdgeColor','r')
% Dibuja las circunferencias intersección del elipsoide con
z=a y con z=0
h2=plot3(sin(t)/2,cos(t)/2,a*ones(size(t)),'r');
set(h2,'LineWidth',3);
h3=fill(sin(t)/2,cos(t)/2,'b');
set(h3,'LineWidth',3); axis equal;
xlabel('X');ylabel('Y');zlabel('Z');
title('Flujo de rot(F)=6k, a través del elipsoide
2x^2+2y^2+z^2=1,con z > = 1/sqrt(2)');
hold off
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
45

El campo vectorial F(x , y, z )  3yi  3xj  z 4 k es de clase C1 en  3 . El elipsoide S es una

superficie suave, limitada por la curva C, orientada en el mismo sentido que la superficie S.

1 1
El elipsoide 2x 2  2y 2  z 2  1 , tiene semiejes , y 1 lo cual se comprueba escribiendo
2 2
x2 y2
su ecuación en la forma 2
 2
 z2  1
 1   1 
   
  
 2 
2
Según el teorema de Stokes, se cumplirá

I    3ydx  3xdy  z 4dz   rot F  ndS

C S

a) Cálculo de la integral de línea.

Se obtiene la ecuación de la curva C como la intersección del elipsoide y el plano z = 1 / 2 , y

posteriormente, se parametriza. Como
2x 2  2y 2  z 2  1
  1 1
C   1
  2x 2  2y 2  1  x 2  y2 
 z   2 4
 2 

1 1
C es la circunferencia de radio , en el plano z  . Su parametrización es:
2 2

 1 1 1 

C 
x  cos t; y sen t; z 
 t  0, 2 

 2 2 2

 
 1 1 
dx   sent; dy  cos t; dz  0
 2 2 

Por lo tanto, la integral de línea será:

2
 3      2
3
 sen t   1 sen tdt    3 cos t   1 cos tdt   3
  3ydx  3xdy  z 4dz    2 
 2 
  2

2



 4
dt 
2
C 0 0

b) Cálculo de la integral de superficie.

Hallamos previamente rot F y el vector n , sabiendo que la ecuación implícita de S es
F x , y, z   2x 2  2y 2  z 2  1  0

i j k
  
rot F   F   6k
x y z
3y 3x z4

el vector n será uno de los dos siguientes:

46 T3  INTEGRACIÓN DE SUPERFICIE

fx' i  fy' j  k
n
   f 
2 2
1  fx' y
'

Tomamos el vector que tiene la tercera componente mayor que cero ya que el vector n apunta
hacia arriba.

fx' i  fy' j  k 6
rot F  n  6k  ' 2 ' 2

1 f  f
x y
1  fx' 2 fy' 2

6
I   rot F  ndS   1  fx' 2 fy' 2
dS
S S

Sustituyendo el diferencial de superficie, dS  1  fx' 2 fy' 2 dxdy , y considerando D la

proyección de S sobre el plano XY, se tendrá:

6
I   1  fx' 2 fy' 2
dS  6  dxdy
S D

3
I  6 área D  
2

1
ya que Sxy es el círculo de radio .
2

25
Utilizar el teorema de Stokes para calcular  2ydx  3xdy  z 2dz , siendo C la
C
circunferencia de ecuaciones paramétricas


C  x  3 cos t, y  3 sen t, z  0, 0  t  2 
Solución
El teorema de Stokes puede aplicarse sobre cualquier superficie S que tenga por borde la curva C
anterior, siendo S suave o suave a trozos, estando S orientada según la orientación de C. El campo
vectorial que interviene en este caso, F(x , y, z )  2yi  3xj  z 2 k , es de clase C1 en todo  3 ,
luego verificará

I   2ydx  3xdy  z 2dz   rot F  ndS

C S

% Dibuja una muestra del campo rotF=k

x=linspace(-3,3,4);y=x;z=linspace(0,3,2);[X,Y,Z]=meshgrid
(x,y,z);
U=0*ones(size(X));V=0*ones(size(X));W=ones(size(X));
quiver3(X,Y,Z,U,V,W); hold on

% Dibuja la semiesfera x^2+y^2+z^2=9, en z>=0

f=linspace(0,pi/2,15); % ángulo phi
t=linspace(0,2*pi,30); % ángulo theta
[F,T]=meshgrid(f,t);
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
47

h1=surf(3*sin(F).*cos(T),3*sin(F).*sin(T),3*cos(F));
set(h1,'FaceColor','blue','FaceAlpha',0.2,'EdgeColor','b')

% Dibuja la circunferencia x^2+y^2=9

t=linspace(0,2*pi,40);
h2=plot(3*cos(t),3*sin(t),'r');
set(h2,'LineWidth',3); axis equal
xlabel('X');ylabel('Y');zlabel('Z');
title('Flujo de rotF=k a través de la semiesfera x^2+y^2+z^2=9,
z>=0');
hold off

2 2 2
Flujo de rotF=k a través de la semiesfera x +y +z =9, z>=0

2
Z

2
2
0
0
-2
-2
Y X

Hallamos rot F siendo F  2yi  3xj  z 2k

i j k
  
rot F   F  k
x y z
2y 3x z 2

Tomamos como S la semiesfera que tiene por ecuación implícita x 2  y 2  z 2  9 con z  0 ,

es decir, z  f x , y   9  x 2  y 2 . Esta superficie tiene como curva frontera C.

El vector n a S será uno de los dos siguientes:

fx' i  fy' j  k
n
1  fx'2  fy'2

Tomamos como vector normal n el que tiene la tercera componente mayor que cero ya que está
dirigido hacia arriba.
Se tendrá que:
 f ' i  f ' j  k 
  '2
I   rot F  ndS      x '2 y '2
k  '2
 fx  fx  1 dxdy   dxdy
S D  fx  fx  1  D

siendo D la proyección de S sobre el plano XY que es el círculo de radio 3, por lo tanto,

48 T3  INTEGRACIÓN DE SUPERFICIE

 dxdy    3
2
I   9
D

También se podría haber utilizado como superficie para aplicar el teorema de Stokes, S1, el círculo
de radio 3, contenido en el plano z  0 , pues también tiene por borde la curva C. En este caso,
el vector unitario perpendicular a S1, orientado según la orientación de la curva C es k ,
resultando

I   rot F  ndS 1
  k  kdxdy   dxdy  área S 1
 9
S1 S1 S1

Código Matlab:
% Dibuja una muestra del campo
rotF=k
x=-3:3;y=-3:3;z=-1:1;
[X,Y,Z]=meshgrid (x,y,z);
U=0*ones(size(X));V=0*ones(size
(X));W=ones(size(X));
quiver3(X,Y,Z,U,V,W); hold on

% Dibuja la circunferencia
x^2+y^2=9 en z=0
t=0:0.05:2*pi;
h2=plot(3*cos(t),3*sin(t),'r');
set(h2,'LineWidth',3); axis
equal
xlabel('X');ylabel('Y');zlabel(
'Z');
title('Flujo de rot(F)=k a
través del círculo x^2+y^2=9');
hold off

Material de consulta

Libro Interactivo. Parte II Capítulo 7.

Laboratorios y Unidades didácticas con videos explicativos: enlace

Autora: Elena Álvarez

 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
49

Ecuaciones de algunas superficies frecuentes

PLANO
Cartesianas

Ax  By  Cz  D  0

Paramétricas

x  x o  ua1  vb1
y  yo  ua2  vb2
z  zo  ua 3  vb3 
x  1v


y  u u ∈ I, v∈ J




z  u v


ESFERA

Cartesianas

x2 + y 2 + z 2 =
r2

Paramétricas

x  r sen u cos v
y  r sen u sen v
z  r cos u
x  5 sen u cos v

y  5 sen u sen v


z  5 cos u

0  u  , 0  v  2
50 T3  INTEGRACIÓN DE SUPERFICIE

CILINDRO ELÍPTICO
Cartesianas

x 2 y2
 1
a 2 b2

Paramétricas

x  a cos u
y  b sen u
z v x  2 cos u

y  4 sen u 0  u  2, v  I


z  v


CILINDRO HIPERBÓLICO

Cartesianas

x2 y2
− 1
=
a2 b2

Paramétricas

x  2 cosh u



y  4 senh u

u, v   D
x = a cosh u 
z v


y = b senh u
z=v
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
51

CILINDRO PARABÓLICO

Cartesianas

x2
z
a

Paramétricas
x u
y v
u2
z u2
a x  u, y  v, z  , (u, v )  D
4

CONO
Cartesianas

z 2  x 2  y2

Paramétricas
x  u cos v

y  u sen v 0  v  2, u  I


z  u


ELIPSOIDE
Cartesianas

x2 y2 z2
+ + 1
=
a2 b2 c2

Paramétricas

x  a sen u cos v x  5 sen u cos v


y  b sen u sen v y  3 sen u sen v


z  c cos u z  2 cos u

0  u  , 0  v  2
52 T3  INTEGRACIÓN DE SUPERFICIE

PARABOLOIDE ELIPTICO
Cartesianas

Paramétricas


x  au cos v 


 x  u cos v
y  bu sen v 
y  2u sen v
  u  I , 0  v  2

 
z  u 2 
z  u2

 


PARABOLOIDE HIPERBÓLICO
Cartesianas

 x2 y 2 
z=
± 2 − 2 
a b 


Paramétricas


x  au cosh v 

 x  3u cosh v

y  bu senh v 
y  2u senh v
  u  I, v  J

 

z  u 2 
z  u 2
 


HIPERBOLOIDE DE UNA HOJA

Cartesianas

x 2 y2 z 2
  1
a 2 b2 c2

Paramétricas
 CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
53


x  a cosh u cos v 
x  2 cosh u cos v

 

y  b cosh u sen v y  4 cosh u sen v   u  , 0  v  2
 

 


z  c senh u 
z  3 senh u
 

HIPERBOLOIDE DE
DOS HOJAS
Cartesianas

x 2 y2 z 2
   1
a 2 b2 c2

Paramétricas

x  a senh u cos v


y  b senh u sen v


z  c cosh u 
 
x  senh u cos v

y  senh u sen v




z  cosh u

  u  , 0  v  2