Onocimientos Revios: Ntegración de Uperficie

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 las integrales dobles y triples.
• Ecuaciones paramétricas de curvas.
• Sistemas de coordenadas cilíndricas y esféricas.
• Representación 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 Superficies
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 a una superficie en un punto
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. Superficie en explícitas. 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)
Prof. Elena E. Alvarez Saiz Universidad de Cantabria

CÁLCULO II – GRADO EN INGENIERÍA MECÁNICA
3
Herramienta visualización plano tangente

https://proyectodescartes.org/miscelanea/materiales_didacticos/Tangente-JS/index.html
Caso 2. Superficie en implícitas. Si la superficie viene dada por una función implícita F (x , y, z )  0 ,
1
con F de clase C 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
• 
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. Superficie en paramétricas. 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.
CÁLCULO DE LA DIFERENCIAL DE UNA SUPERFICIE
Nuestro objetivo es definir el área de una superficie S que está definida sobre un dominio D que en
la imagen siguiente por simplicidad se ha considerado un rectángulo. Haciendo una partición del
dominio D en subrectángulos, se obtiene una división de la superficie. La aproximación del área de
este elemento de superficie se aproximará por una porción de plano tangente.

5
La herramienta siguiente permite mostrar la relación entre el área un paralelogramo y su proyección

sobre el plano XY y que puede facilitar la comprensión del concepto de diferencial de superficie si
consideramos que el paralelogramo es un “trozo” de plano tangente que se proyecta sobre el
rectángulo del plano z  0 , de área R .
Herramienta
https://proyectodescartes.org/miscelanea/materiales_didacticos/ProyeccionParalelogramo
-JS/index.html
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 C 1 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 
1
Y puesto que, cos   n  k  resulta finalmente,
N
dS  N dxdy  fx(x , y )2  fy(x , y )2  1 dxdy
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,

7
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
2 2 2
dxdy  (x , y )  (y, z )   (x , z ) 
dS           dudv  ru  rv dudv
cos   (u, v )  (u, v )  (u, v )
  


Herramienta visualización área de una superficie:

https://personales.unican.es/alvareze/CalculoWeb/CalculoII/escenas/DiferencialSuperficie-
JS/index1.html
INTEGRAL DE SUPERFICIE DE UN CAMPO ESCALAR
4 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 donde,

S
• 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
1
Una superficie suave es imagen de una función de clase C 1 definida en un dominio D .

9
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 .
Suponiendo S definida en paramétricas

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

5 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.
6 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.
INTEGRAL DE SUPERFICIE DE UN CAMPO VECTORIAL O INTEGRAL DE FLUJO
7 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

11
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.
Herramienta interpretación integral de superficie de un campo vectorial

https://proyectodescartes.org/miscelanea/materiales_didacticos/IntegralSuperficieCVector
ial-JS/index.html

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
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
8 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.

13
9 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

frontera que encierra el sólido H.
El flujo que atraviesa la superficie S 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.
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
primer octante con 0  z  2 . Resuelve la

1
En el resumen de teoría se recogen las integral con Matlab.
expresiones del elemento diferencial de Solución: b) área = 11, 785
superficie de una superficie dada por
z = f (x , y ) con (x , y )  R , estando R en el
3
plano XY . Determina el elemento diferencial de
a. Escribe la expresión del elemento superficie para las siguientes superficies e indica
diferencial de superficie adaptada a la en cuáles, el cociente entre el elemento
superficie dada por y = g(x , z ) con diferencial de superficie y el elemento
(x , z )  R estando R en el plano XZ . diferencial de área es constante. Prepara un
fichero con Octave/Matlab para dibujar las
b. Encuentra dS para la superficie
superficies de este ejercicio, así como una
x = 1  y 2  2z 2 , tomando como muestra de sus vectores normales.
variables independientes con a) 2x  3y  5z  1 b) 2x  3y  1
y, z   {(y, z ) / y 2
 2z 2 < 1} c) z = 9  x 2  2y 2 d) z = 3 x 2  y 2
Solución: a) e) z = x 2  3y 2 f) z = 4  x 2  2y 2
2 2
dS  g 'x (x , z )  g 'z (x , z )  1 dxdz g) x 2  2y  4
1  2z 2 38 dS
b) dS  dydz Solución: a) dS  dxdy ,  cte ;
2
1  y  2z 2 5 dA
13 dS
b) dS  dxdz ,  cte ;
3 dA
2
En el resumen de teoría se recoge la 9  2y 2
expresión del elemento diferencial se superficie c) dS  dxdy ;
de una superficie dada por F (x , y, z ) = 0 si 9  x 2  2y 2
dS
Fz  0 . d) dS  10 dxdy ,  cte ;
dA
a. ¿Cómo se calcularía el elemento diferencial
de superficie si Fz  0 ? x 2  6y 2
e) dS  2 dxdy ;
b. Determina dS para la superficie x 2  3y 2
x 2  5y 2 = 25 y calcula el área de la f) dS  1  4x 2  16y 2 dxdy ;
porción de ésta superficie situada en el

15
g) dS  1  x 2 dxdz a) S es la porción del plano z = x  2y que se

proyecta en el rectángulo R   0,1   0, 3 del
plano XY ; g x , y, z   xyz
4
Halla el área de la superficie de las b) S es la parte de la superficie
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
z = x 2  y  1 comprendida entre los
2
b. Solución: a) Masa = 21 , densidad media=
2
planos z=1 y z=4

c. superficie esférica x 2  y 2  z 2  2 , b) Masa =
60
 
7  55/2  41 , densidad media=
situado por encima del plano z  1 .
7  55/2  41
d. cilindro x 2  z 2  a 2 de altura h.  3, 44
e. esfera x 2  y 2  z 2  a 2 
10 53/2  1 
Solución: (a)   33  17 
7
b)

 17 7  5 5  d) 2ah e) 4a 2 
Una lámina S ocupa la parte del plano
x  y  z  a en el primer octante ( a es un
6 número real positivo).
Si la temperatura en cada punto de S es
5
(a) Calcula la integral de superficie
 
T x , y, z   k x 2  y 2  z 2 ( k  0 ), halla el
valor medio de la temperatura en S .
  (x
2
I =  y 2 )dS
S 1
Supón que T x , y, z   ; halla el valor
x2 y2 z2 x y 1
donde S es   =0 para
a2 a2 b2 medio de la temperatura.
0  z b Solución:
(b) Calcula la integral de superficie ka 2
a) Tmedia 
  (x
2 2
I =  y )dS donde S es la esfera 2
S
2
x 2  y2  z 2 = a2 b) Tmedia 
a2
a  log a  1
2
z x
(c) Calcula I =  S
dS donde S es la
1  4x 2
porción de la superficie z  4  x 2 limitada por 8
el plano x  2y  4 Se considera la superficie S dada por
(d) Calcula la integral x = a  y 2  z 2 y el campo escalar
2
 4 
I =   z  2x  y  dS donde S es la
S 3 
 
g x , y, z   x y 2  z 2 . Halla el valor medio de

g en S
x y z 3
porción del plano    1 situada en el a
2 3 4 Solución: gmedio 
primer octante. 4
Solución: (a) I = 8a 4 / 3 (b) I  8a 4 / 3 9

Calcula el flujo saliente del campo
(c) 12 (d) 4 61 vectorial F = 2i  5 j  3k a través de la parte
del cono z = x 2  y 2 interior al cilindro
6
Encuentra la masa de la lámina S , su x 2  y2 = 1 .
densidad media si en cada punto la densidad de Solución: Flujo  3
masa superficial de la placa es la función
g x , y, z  indicada.

sólido encerrado por las superficies

10
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
hemisferio 0  z  9  x 2  y 2
parte del plano x  y  z  a en el primer
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,
c. F(x , y, z ) = (x  z 2 )i  (y  z 2 ) j  xk ; H
F  c T , c  0 para las funciones de
es el sólido 0  y 2  z 2  1 , 0  x  2 .
temperatura
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   15
x y 1 Encuentra una expresión para la tercera
Calcula el flujo de calor hacia arriba a través de componente, P x , y, z  , para que el flujo del
la placa.
campo vectorial
Solución: a) Flujo cka 3  0
F x , y, z  = xy  z  i  ze y j  P x , y, z  k
 a 
b) Flujo 2c log a  1   hacia arriba del semielipsoide
 a  1
z  9  2x 2  y 2 sea opuesto al flujo hacia
debajo de la región D del plano XY
12 encerrada en él y que éste último flujo sea el
Se considera la superficie S porción del
doble del área de D .
paraboloide z = a 2  x 2  y 2 para x , y, z
z 2 y
positivos; a es un número real no nulo. Un Solución: P x , y, z   e  yz  2
2
campo vectorial F depende de a de la
2
siguiente forma: F x , y, z   axi  yj  zk 16
a Halla el flujo del rotacional de
a. Calcula a mano el flujo de F hacia arriba de
F x , y, z  = xy 2 i  j  e x k sobre la parte de
2
S
b. Analiza el comportamiento del flujo en z  16  2x 2  y 2 en z  0 , orientada con la
función del parámetro a .¿Es creciente el normal hacia arriba.
flujo con el valor de a ? ¿Existe un valor Solución: Flujo=0
mínimo para el flujo? Puedes dibujar con
Matlab la función flujo y utilizar el comando
solve para hacer el análisis pedido. 17
3
Sea S la superficie de ecuación z=2
a
Solución: Flujo=
8
 
a 2  a  2 b) El flujo definida sobre el círculo x 2  y 2  4 y se
considera el campo
alcanza un valor máximo relativo para
a  1, 5662 y un valor mínimo relativo para F x , y, z  = 3yi  xzj  yz 2 k . Calcular  F  dr
a  0, 7662 . C
usando la definición de integral de línea y el
Teorema de Stokes.
13 Solución: 20
Considerando F x , y, z  = xi  yj  zk
, calcular el flujo que atraviesa la superficie del

17
18 siendo C la curva x = a sen t , y = a cos t ,

Sea C el triángulo orientado situado en
z = a(sen t  cos t ) , 0  t  2 .
el plano 2x  2y  z  6 , intersecado por los
planos coordenados y Solución: I = a 2
F x , y, z  = y 2 i  zj  xk . Calcular  F  dr
C 21
usando la definición de integral de línea y el Calcular la circulación de
Teorema de Stokes. F(x , y, z ) = (x  3z )i  (z  y ) j  4yk
Solución: -9. alrededor del trapecio de vértices
A  (2, 2, 0) , B  (2, 2, 0) , C  (1,1,1) y
19 D  (1, 1,1) tomando como orientación la
Calcula el trabajo realizado por el campo que recorre los vértices en el orden en el que se
F(x , y, z ) = x 2y 3 i  j  zk a lo largo de la han definido.
circunferencia x 2  y 2 = R 2 en z = 0 . Solución: C=18
Solución: Trabajo = R 6  / 8
22
Comprobar el teorema de Stokes para
20 F x , y, z  = 2zi  xj  y 2 k considerando como
Calcula la siguiente integral mediante la
fórmula de Stokes y después comprueba el superficie S el paraboloide z  4  x 2  y 2 y C
resultado integrando directamente es la curva de nivel de S sobre el plano XY.
I =  xdx  (x  y)dy  (x  y  z )dz Solución: 4
C
Test de autoevaluación
C. 0
1 D. Ninguna de las anteriores.
El área de la superficie de ecuación
z  16  x 2  y 2 situada sobre la región

R  x 2  y 2  16 es,  3
Sea  una superficie cuya ecuación en
A.

6
65  1  coordenadas cilíndricas es z  f (r , ) , siendo
D la proyección de la superficie en el plano r ,
B.

6

653  1  entonces
2 2
 z   z 
C.

3

653  1  A. Área de   
D
1  r 2      dr d 
 r    
D. Ninguna de las anteriores.  z 
2
1  z 
2
B. Área de    1       r dr d 

 r  r   
D
2 2
 z  1  z 
2
La temperatura media de la lámina S C. 1     2   r dr d  .

Área de     r  r   

x 2  y 2  1, y  0 D
definida por   y con D. Ninguna de las anteriores.


0z 3


temperatura en cada punto T (x , y, z )  y , es:
4
2 El área del casquete de la superficie
A.
 esférica x 2  y 2  z 2  2 , situado por encima
6 del plano z  1 , es:
B.


A. 2 2  2  1 B.   F nd   0


B. 4 2  2  1 . C.  div F nd   0


C.  2  2  1 . D. Ninguna de las anteriores.
D. Ninguna de las anteriores.

8
Sea V(x , y, z )  a  r , siendo a un
5
Se consideran las siguientes integrales: vector constante y r el vector de posición de un
I1   g(x, y, z )d  , I2   d  , punto de  3 y sea C una curva de  3 que es
  frontera de la superficie  . En estas condiciones
I3   Fd σ , I4   divFdV , se verifica:
 V
A.  Vdr  2 a  n d  .
I5   FTds C 
 Vnds   rot V n d  .
C
y las siguientes magnitudes: B.

C 
a. Masa de una lámica C.  Vdr  2 rot r n d  .

b. Trabajo C 
c. Circulación D. Ninguna de las anteriores.

d. Flujo
e. Área superficial
9
Sea  el paraboloide de ecuación
Las integrales y las magnitudes admiten los
siguientes emparejamientos: z  x  y 2 con 0  z  3 y sea F el campo
2
vectorial F  xi  yj  zk . Entonces el flujo de

A. I ,e, I , a , I , d , I ,b
1 2 3 5
rot F a través del paraboloide en el sentido de la
B. I , a , I ,e, I , d , I , c .
1 2 4 5
normal exterior es:
C. I , a , I ,e, I , c, I , d  .
1 2 3 4
A.  .
B.  .
D. Ninguna de las anteriores. C. 0 .
6
El flujo del campo F  2xi  xzj  z 2 k
hacia el exterior de la esfera unitaria, es:
10
A. 0 El valor de la integral de flujo
B.
  ( F)d σ , siendo S la semiesfera
S
3
8 x 2  y 2  z 2  16, z  0 y F el campo
C.
3   
F  x 2  y  4 i  3xj  2xz  z 2 k , es: 
A. 32 .
B. 0 .
7
Sea  una superficie cerrada que C. 2 .
encierra un volumen V y F un campo vectorial D. Ninguna de las anteriores.
de clase C 2 , entonces
A.  div FdV  0
V
1 2 3 4 5 6 7 8 9 10

19
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
SUPERFICIES
1
Representar las siguientes superficies con Matlab
x  u cos v
2u 2 2
(a) y  u senv
0    2
z  9  u2
(b) Porción de la superficie z = x 2  y  1 comprendida entre los planos z=1 y z=4

2
(c) Superficie esférica x 2  y 2  z 2  2 situada por encima del plano z  1 .

(d) Cilindro x 2  z 2  a 2 de altura h.
(e) Esfera x 2  y 2  z 2  a 2
Solución a)
syms u v
fsurf(u*cos(v),u*sin(v),9-u^2,[2 2*sqrt(2) 0 2*pi])
hold on
fsurf(0*u,[-2 2 -2 2])
hold off
xlabel('X')
ylabel('Y')
zlabel('Z')
axis equal

Solución b)
syms u v
fimplicit3(x^2+(y-1)^2-z,[-3 3 -1 3 1 4])
hold on
fsurf(0*u)
hold off
xlabel('X')
ylabel('Y')
zlabel('Z')
axis equal
Solución c)
syms x y z
hold on
fimplicit3(x^2+y^2+z^2-2,[-1 1 -1 1 1 4])
fsurf(0*x,[-1 1 -1 1])
hold off
xlabel('X')
ylabel('Y')
zlabel('Z')
axis equal

21
Solución d)
syms x y
a=3;h=5;
fimplicit3(x^2+y^2-a^2,[-a a -a a 0 h])
xlabel('X')
ylabel('Y')
zlabel('Z')
axis equal
Solución e)
syms x y
a=3;
fimplicit3(x^2+y^2+z^2-a^2,[-a a -a a -a a])
xlabel('X')
ylabel('Y')
zlabel('Z')
axis equal

DIFERENCIAL DE SUPERFICIE. ÁREA DE UNA SUPERFICIE
2
Considera los vectores V1  1.1, 0, 0 , V2  0,1.3, 0 y los vectores

• T1 : vector del plano XZ que se proyecta sobre V1 y forma con él   radianes.
6

• T2 : vector del plano YZ que se proyecta sobre V2 y forma con él   radianes.
8
Se pide:
a) Calcular el área del paralelogramo formado por V1 y V2

b) Calcular el área del paralelogramo formado por los vectores T1 y T2.
c) Observa que la relación entre las dos áreas es el módulo del vector normal a los vectores
T1 y T2.
Solución
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
para comprobar el resultado.
para observar cómo se obtiene la diferencial de superficie cuando la superficie viene dada en
explícitas.
3
Dada la superficie S definida por z  x 2  y 2 sobre el rectángulo R  1, 2  1, 3 se
   
pide:
a) Representar la superficie S sobre R.

b) Obtener a mano la ecuación del plano tangente en el punto (1,1).

23
c) Representar en la misma figura el vector normal N a la superficie en el punto (1,1) y los


vectores directores del plano tangente. T1  1, 0, fx' 1,1   
T2  0,1, fy' 1,1
Solución
syms x y z
fimplicit3(x^2+y^2-z,[1 2 1 3 0 13])
%Plano tangente
%z=2+2(x-1)+2(y-1)
hold on
fimplicit3(-2+2*x+2*y-z,[1 2 1 3 0 8])
quiver3(1,1,2,1,0,2,'linewidth',5)
quiver3(1,1,2,0,1,2,'linewidth',5)
quiver3(1,1,2,-2,-2,1,0.3,'linewidth',5)
plot3(1,1,2,'o')
xlabel('X')
ylabel('Y')
zlabel('Z')
hold off
4
Considerando la superficie S del ejercicio anterior z  f x , y   x 2  y 2 sobre el
rectángulo R  1, 2  1, 3 se pide:
   
a) Calcula el área de la porción T de plano tangente en el (1,1) que se proyecta sobre el
rectángulo R.

b) Comprueba área T   N  área R  siendo N  fx' 1,1, fy' 1,1, 1 
Solución
 
Se tiene que área R   2 , N  fx' 1,1, fy' 1,1, 1  2, 2,1 , la porción T está generada

por los vectores T1  1, 0, fx' 1,1, T 2  
 2 0,1, fy' 1,1 . Por lo tanto, su área es el
   
área T   1, 0, fx' 1,1  0, 2, 2 fy' 1,1  2  3

5
Calcular con Matlab el área de la superficie z = x 2  y  1 comprendida entre los
2
planos z=1 y z=4. Representar la superficie.
La parametrización de la superficie es
x  u cos v
2u 2 2
y  u senv
0    2
z  9  u2
Para representar la superficie en paramétricas
syms u v
fsurf(u*cos(v),u*sin(v),9-u^2,[2 2*sqrt(2) 0 2*pi])
%Dibujamos la frontera del dominio en el que se
%proyecta la superficie
hold on
fplot3(2*cos(v),2*sin(v),5+0*v,'linewidth',4)
fplot3(2*cos(v),2*sin(v),0*v,'linewidth',4)
fplot3(sqrt(8)*cos(v),sqrt(8)*sin(v),0*v,'linewidth',4)
fplot3(sqrt(8)*cos(v),sqrt(8)*sin(v),1+0*v,'linewidth',4)
hold off
La superficie S y su proyección sobre el plano XY se muestra en la siguiente figura:
dS  1  2x   2y  dxdy x, y  / 2  x 

2 2
2
D  y2  2 2
2 2 2
área  
S
dS  
D
1  4x 2  4y 2 dxdy 
coordenadas
polares

0 2
1  4r 2 rdrd     33  17 
6
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.

25
7
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
Código Matlab
syms x y z
b=2;
h1=fimplicit3(x^2+y^2-z^2,[-b b -b b -b b]);
set(h1,'FaceColor','interp','FaceAlpha',0.2,'EdgeColor','none');
hold on
% Dibuja de la segunda superficie
h2=fimplicit3(x^2+y^2-2*x,[-b b -b b -b b]);
set(h2,'FaceColor','interp','FaceAlpha',0.5,'EdgeColor','none')
xlabel('X');ylabel('Y');zlabel('Z');
%title('')
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 .
S
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  dxdy  2dxdy
dS  1  fx'2  fy'2 dxdy  1    2
 x 2  y 2  2 
   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
 2  4
D
a 2 2a 2
Área= 2 2  u. d. s.
4 2
8
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
Se trata de calcular el área de la esfera acotada 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

27
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  3x  0 que llamamos D.
3
Sustituyendo, área = 2  dxdy
2 2
D 9 x y
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
 /2 
 3 9  r 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
9
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
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 en el ejercicio.
u  r 2  du  2rdr 
 
 3/2 
1  4r dr    
r  
3 2 2
 2 1 1  4r 
dv  r 1  4r dr  v  
 8 3/2 
1 1  4r 
5/2
2
2 2
r 1 r
     
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

29
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 
10
Determinar el área de la porción de superficie de la semiesfera
x  y  z 2  4a 2
2 2
z  0 que es interior al cilindro x 2  y 2  2ax  0
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

S
Código Matlab:
syms x y z
h1=fimplicit3(x^2+y^2+z^2-4,[-2 2 -2 2 0 2]);
hold on
h2=fimplicit3(x^2+y^2-2*x,[-2 2 -2 2 0 2]);
%title('')
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, dS  1  fx'2  fy'2 dA
2 2
   
 x   y   .
'2
1 f  f '2
 1 
   
x y  2  2 
 4a  (x  y ) 
2 2
 4a  (x  y ) 
2 2

2a 2adxdy

4a 2  (x 2  y 2 )
 Área   4a 2  (x 2  y 2 )
D
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
11
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

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

dos superficies:
 2


x 2  y 2  3z 2  0 z  1 x 2  y 2  4a 2 
  1
 2
x  y2  z 2 
4a 2 
 2
3
4a 2  
 3
 x 2  y2  x 2  y2
3
 

 z   x 2  y2 

 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
dS  1  fx'2  fy'2 dA
2 2
   
 1 x  1 y
1 f  f '2 '2
 1       
x y  3 x 2  y 2   3 x 2  y 2 
 
x2 y2 4 2
 1   
 2
3 x y 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
2
área del dominio D, cuyo valor es πa por ser un círculo, obteniendo:
2 2 3
A= a 2  a 2 u. d. s.
3 3
Código Matlab
syms x y z
h1=fimplicit3(x^2+y^2-3*z^2,[-5 5 -5 5 0 sqrt(3)]);
hold on
h2=fimplicit3(x^2+y^2+z^2-12,[-5 5 -5 5 sqrt(3) sqrt(12)]);
%title('')
hold off

12
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 dxdy 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
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,
x 2  y 2  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

33
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
13
Calcula la temperatura media T x , y, z   x 2  y 2  z 2 sobre la porción del cono
z 2  x 2  y 2 que se encuentra entre los planos z  2 y z  3 .
Solución
Para calcular la temperatura media se tendrá en cuenta que
 T x, y, z dS  T x, y, z dS

Tm  S
 S
área S   dS
S
En este caso la superficie S es la parte del cono z  f x , y   x 2  y 2 entre los planos z  2

y z  3.
El diferencial de superficie es dS  1  fx'2  fy'2 dA siendo
2 2
   
 x    y   1  x2 y2
1  fx'2  fy'2  1      2
2  2 
 2
 x  y 
 2
 x  y  x 2  y2 x 2  y2    
El dominio D, proyección de esta superficie S sobre el plano z=0, es la corona circular determinada
por las circunferencias de radios 2 y 3.

Se tiene entonces que el área de la superficie es
 dS  
2  dxdy  2 32   22   5 2 
S D
La temperatura total sobre S es
 T x, y, z dS   T x, y,

S D
x 2  y2  2 dA   2 x
D
2
 y2  2 dA
Pasando a coordenadas polares

2 3 r 3
r4
 T x, y, z dS  2 2    2 81  16  65 2
2
r  r drd   4 2
S 0 2
4 r 2
Por lo tanto, la temperatura media es
 T x, y, z dS 65 2
Tm  S
  13º
área S  5 2
14
Resolver con Matlab los siguientes apartados. (a) Representar la superficie que viene dada
por las ecuaciones paramétricas r(u, v )  u cos(v )i  u sen v j  2v k 0  u  3, 0  v  2
(b) Calcula el área de la rampa helicoidal dada por las ecuaciones paramétricas.
(c) Considerando la función g x , y, z   x  yz que representa la función temperatura en
un punto del espacio, se pide calcular la temperatura total en S.
Solución
Dada la superficie
se define: 
r u, v   x u, v , y u, v , z u, v  

ru' u, v   x u' u, v , yu' u, v , z u' u, v  
35
 
rv' u, v   x v' u, v , yv' u, v , z v' u, v 
Para calcular el producto vectorial de dos vectores ru y rv puedes utilizar el comando cross.
Por ejemplo,
ru=[1,2,r];rv==[2,-1,3];
cross(ru,rv)
Recuerda que  dS   ru  rv dudv

S D
Solución apartado a) b)
%Representación
syms u v
fsurf(u*cos(v),u*sin(v),2*v,[0 3 0 2*pi])
xlabel('X')
ylabel('Y')
zlabel('Z')
axis equal
% Cálculo del área
x=u*cos(v)
y=u*sin(v)
z=2*v
r=[x,y,z]
ru=diff(r,u)
rv=diff(r,v)
prodv=simplify(cross(ru,rv))
dS=sqrt(prodv(1)^2+prodv(2)^2+prodv(3)^2)
area=int(int(dS,u,0,3),v,0,2*pi)
double(area)
Solución apartado c)
Temperatura = -(4*pi*(13*13^(1/2) - 8))/3

15
Sea S la porción del paraboloide z = 4  x 2  y 2 limitada entre los planos z = 0 , y = x
y x = 0 . Una lámina tiene la forma dada por S y la temperatura en cada punto es proporcional
al cuadrado de la distancia al eje OZ . Crea un fichero con Matlab para resolver los siguientes
apartados:
(a) Representar la superficie S utilizando las ecuaciones paramétricas

x  r cost, y  r sen t, z  4  r 2 .
(b) Encontrar la temperatura promedio de la lámina, tomando el valor k  1 para la
constante de proporcionalidad de la temperatura.
(c) Determinar y representar sobre la superficie, el lugar geométrico de los puntos de la
lámina en los que se alcanza la temperatura promedio para k  1 .
Solución
• La función temperatura es, T (x , y, z )  k (x 2  y 2 ) ; donde K es la constante de
proporcionalidad.
• la diferencial de superficie es, dS  4x 2  4y 2  1 dA ;
• la proyección de la porción de paraboloide sobre el plano z  0 es un sector circular,
que podemos definir en coordenadas polares como,

  

D 
(r , ) / 0  r  2,  



 4 2


• las integrales para calcular el área y la temperatura global en coordenadas polares son,
 /2 2  /2 2
área =   1  4r r drd  ;
2
temp = k   1  4r 2 r 3 drd 
 /4 0  /4 0
y la temperatura media pedida es, tm=temp/area
%Dibuja la porción del paraboloide z=4-x^2-y^2, situado en

%el primer octante, entre los planos y=x, x=0, z=0,
%coloreado según T=x^2+y^2
t=linspace(pi/4,pi/2,10);
r=linspace(0,2,20);
[T,R]=meshgrid(t,r);
X=R.*cos(T); Y=R.*sin(T);Z=4-R.^2;
% dibujo de la superficie, el cuarto argumento es la

% función que da el mapa de color
surf(X,Y,Z,X.^2+Y.^2)
xlabel('x'); ylabel('y'); zlabel('z');
shading interp
title('paraboloide z=4-x^2-y^2 coloreado según T=x^2+y^2')
axis equal
hold on
%Dibujo del dominio

syms u v

37
fsurf(u*cos(v),u*sin(v),u*0,[0 2 pi/4 pi/2])
%Cálculo de la temperatura promedio con k=1

syms u v
area=double(int(int(sqrt(1+4*v^2)*v,v,0,2),u,pi/4,pi/2));
Temp=double(int(int(sqrt(1+4*v^2)*v^3,v,0,2),u,pi/4,pi/2));
Tm=Temp/area
%Los resultados que se obtienen son:

%area=4,5221; Temp=10,5579; Tm=2,3347
%Representamos los puntos (x,y) que cumplen T(x,y)=Tm

syms t
xt=sqrt(Tm)*cos(t);
yt=sqrt(Tm)*sin(t);
zt1=0*xt;
zt=4-xt^2-yt^2;
fplot3(xt,yt,zt1,[pi/4 pi/2],'linewidth',5)
%Representamos los puntos de la placa que están a temperatura media
fplot3(xt,yt,zt,[pi/4 pi/2],'linewidth',5)
colorbar
hold off
16
Sea S la porción del cono z = x 2  y 2 con z  1 . Una lámina tiene la forma dada por
S y su densidad en cada punto es igual a la altura del punto respecto del plano z  0 . Crea un
fichero con Matlab para resolver los siguientes apartados:
(a) Representar la superficie S utilizando las siguientes ecuaciones paramétricas:

x  r cos t, y  r sen t, z  r
(b) Encontrar la densidad media de la lámina.
(c) Determinar y representar sobre la superficie, el lugar geométrico de los puntos de la
lámina que tienen una densidad igual a la densidad media.
Solución
Adaptar el código del ejercicio anterior.

INTEGRAL DE UN CAMPO VECTORIAL SOBRE UNA SUPERFICIE. FLUJO
17
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.
18
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, dS  NdA , siendo

S
 
 x y
N  z x , z y ,1   , ,1
 2  x 2  y 2 2  x 2  y 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,

39
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
19
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
F  ndS  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
20
Se considera la porción de superficie z  2  3y  x 2 que está por encima del triángulo
T que se encuentra en el plano XY y tiene por vértices los puntos (0,0), (2,0) y (2,-4). Sobre cada
punto de la superficie se define la función temperatura T x , y, z   z  3y  x 2 y el campo
 
F x , y, z   3x i  2z j  1  y 2 k .
a) Calcular la temperatura total y la temperatura media.

b) El flujo hacia abajo del campo a través de esta superficie.
SOLUCIÓN a)

Se pide calcular la temperatura de la superficie z  2  3y  x 2 sabiendo que en cada punto de

la superficie su valor es: T x , y, z   z  3y  x 2
 z  3y  x dS
2
La superficie S se proyecta en el plano XY en el dominio D que se representa en la figura
 z  3y  x dS  2  3y  x   3y  x  2x   3

Temtotal  2
  2 2 2 2
 1 dA 
S D
0
2 0 2 
  2 4x 2  10 dA  2  2x 4x 2  10 dy dx  2  y 4x 2  10  dx 
0 0   2 x
D
2
3  3 3

1 1  2
 
2
 0 4x 4x 2  10 dx  4x 2  10 2  26  10   33.6506
2
3 
3  
0

El área de la superficie es
2 2 Temtotal
área S    dS   2x   3  1 dA   4x 2  10 dA 
S D D
2
La temperatura media es 2.
Solución b)
Se tiene que
 F  n dS   F  N dA   y 
2
 18y  13 dA
S D D
ya que
   
F  N  x , y, 2  3y  x 2  N  3x , 2 2  3y  x 2 ,1  y 2  2x , 3, 1   
 
 6x 2  6 2  3y  x 2  1  y 2     y 2
 18y  13
Calculando la integral doble

0
2 1 
0 2x y 
2 0
0  y  9y  13y 
2 3 2
 18y  13 dy dx  dx 
3  2x

41
2
28   
 x 3  36x 2  26x  dx   2 x 4  12x 3  13x 2    412
0  3 
  3 
 3
0
21
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
S2: cara z = 2, n 2  k ; 2   F  n 2
 zdxdy  2dxdy 2área S 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
22
Se considera la porción del paraboloide z = 1  x 2  y 2 , con x , y, z positivos y un campo
vectorial,  F  d S . Plantear la integral para calcular el flujo hacia el exterior del paraboloide.
S
Solución
El flujo es Flujo   V  n dS   V  N dA donde,

S D
• El campo actuando sobre la superficie es: V  (x ,  2y, 1  x 2  y 2 )
• 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

43
TEOREMA DE LA DIVERGENCIA O DE GAUSS

23
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
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
    2
 d   1  1 sen  d    1   1 cos  
  
 2

3
sen   
r  0
  2 3   2 3 

0 0  0

24
Calcula la integral de flujo saliente del campo vectorial
    
F x , y, z   y 3  3x i  xz  y  j  z  x 4 cos x 2y k
a través de la superficie que delimita el sólido x 2  y 2  1 , x  0 , y  0 , 0  z  1
Solución
    
En nuestro caso, F x , y, z   y 3  3x i  xz  y  j  z  x 4 cos x 2y k es de clase C1 en
3 y la superficie S es cerrada frontera de un sólido H (se trata de la cuarta parte del cilindro de
base la circunferencia unidad y altura 1).
Por tanto, podemos aplicar el teorema de la divergencia o teorema de Gauss:
 (F  n)dS   div F dV 

S H
Se tiene que: div F 

4 2

 (y 3  3x )  (xz  y )  z  x cos x y
 
 
 3 11  5
x y z
El flujo saliente total a través de la superficie cerrada S, será entonces
5
 (F  n)dS   5 dV  5 vol H   4
S H
25
Aplicar el Teorema de Gauss para obtener el flujo del campo vectorial:
F  2x i  2y j  z x 2  y 2 k a través de la superficie frontera del sólido

H  z  0 , 1  x 2  y2  z 2  2 . 
Solución
Utilizando el teorema de Gauss se tendrá que
flujo   F  n dS   div F dV   x 2  y 2 dV

H H H

45
M N P
ya que div F     2  2  x 2  y2  x 2  y2
x y z
El sólido H está delimitado por las esferas x 2  y 2  z 2  1 , x 2  y 2  z 2  2 y el plazo z  0
Pasando a coordenadas esféricas
2  /2 2    4   2  /2
  2 
flujo    sen  sen d  d d    

2
 4 

 1  sen 2 d 
1 0 0 2 2  1 0
x y
 4  1  /2 3
 /2
1  cos 2

   2  sen  d  2 
2
d 
 4  0
4 0
2
3  32
 
2 4 8
26
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 .
c) Aplica el teorema de Gauss a la superficie S1  S 2 comprobando previamente que se
cumplen las hipótesis de dicho resultado.
Solución a)

x 2  y2 x y 
El vector normal a la superficie z  f x , y  
2a
 a a

es N  fx' , fy' , 1   , .  1 . El

conjunto D proyección de S1 sobre XY es el círculo de centro (0,0) y radio 2a .
 2
x , y, x  y    x , y , 1 dA 
2
flujoS 
1
 F  n dS   F  N dA   


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 .
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)
Para calcular el flujo a través de S1  S2 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
 
27
Calcular, utilizando y sin utilizar el teorema de Gauss, el flujo del campo
F x , y, z   x 3 i + y 3 j  z 3k
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: 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.

47
Si se aplica el teorema de Gauss
 F dS   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
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
I    3x z  3y z  z 
2 2 3
z 0
dydx 
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   
   4

2

4
 dx 
   2

2
 dx

0   0  
x a
a 
a  x 

4

3x 2 a 2  2xa  x 2 dx   a  x  5
x 3a 2 3x 4a 3x 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,
 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
x 3 , y 3 , a  x  y  1,1,1 dA  x  y 3  a  x  y   dydx 
 FdS        
S4 D4 0 0
y a x
 4  4 4
a
 3 4
x y  y 
a  x  y   a 
x 3 a  x  
a  x  a  x  
  4 4 
dx    4

4 dx 

0   0  
 y 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
28
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.

49
Parametrización S1
S1  r1 D1  D1  x, y    2
/ x 2  y 2  R2 
r1 x , y   x , y, 3
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
 
 3xdxdy     3r cos d dr  0
2

 0 
D1 0 
Parametrización S2
S2  r2 D2  D2  x, y    2
/ x 2  y 2  R2 
r2 x , y   x , y, 0
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 
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
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
29
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
syms t z
h1=fsurf(2*cos(t),2*sin(t),z,[0 2*pi 0 2]);
set(h1,'FaceColor','b','FaceAlpha',0.2,'EdgeColor','b')
title('vector F=(2x,y,z^2) a través del cilindro x^2+y^2=4')

51
axis equal; hold off
Se cumplen las condiciones para aplicar el Teorema de la divergencia:
 div F dV   F  n dS
V S
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 F dV   (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 2
0 (3  2z )dz  3z  z 2   10 , 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 ,

2
n  k  F  n  z  0 , luego el flujo saliente a través de la cara inferior es nulo.
Podemos escribir:
lat   div F dxdydz   sup

 inf  10a 2  4a 2  0  6a 2 unidades de flujo.
V
30
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

53
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
syms t z
h1=fsurf(2*cos(t),2*sin(t),z,[0 2*pi 0 1]);
set(h1,'FaceColor','b','FaceAlpha',0.2,'EdgeColor','b')
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
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

F  n  0, dS  dxdy , por tanto φinferior 0.
=
Despejando el flujo pedido se obtendrá
lateral  total  superior  inferior  36  27   9 unidades de flujo
31
Sea H el sólido limitado por las superficies z  x 2  y 2 y los planos z  1 y z  2 .

Consideremos el campo F x , y, z   x 2 , y 2 , z 2 . 
a) Representa la superficie S junto con una muestra del campo.
b) Comprobar el teorema de la divergencia.
Solución
%Tapa inferior: S1
r=linspace(0,1,3);t=linspace(0,2*pi,40);
[R,T]=meshgrid(r,t);
X=R.*cos(T);Y=R.*sin(T);Z=0*T+1;
hold on
surf(X,Y,Z)
quiver3(X,Y,Z,X.^2,Y.^2,Z.^2)
%Superficie lateral: S2
X=R.*cos(T);Y=R.*sin(T);Z=sqrt(X.^2+Y.^2);
h=surf(X,Y,Z)
set(h, 'FaceAlpha', 0.4)
%Tapa superior: S3
surf(X,Y,Z)
[R,T]=meshgrid(0:0.5:2,0:1:2*pi);

55
hold off
S1: tapa inferior S2: tapa lateral
S3: tapa superior
Solución b)
15
Flujo a través de S1 :  Flujo a través de S2 :  Flujo a través de S3 : 16
2
32
Utiliza el teorema de la divergencia de Gauss en los tres casos siguientes para calcular
I =  F  n dS . Nota: Este es el ejercicio propuesto número 14.

H

a) F(x , y, z ) = zi  xj  yk ; H es la semiesfera 0  z  9  x 2  y 2
b) F(x , y, z ) = x 2 i  y 2 j  z 2 k ; H es el sólido parabólico 0  z  4  x 2  y 2
c) F(x , y, z ) = (x  z 2 )i  (y  z 2 ) j  xk ; H es el sólido 0  y 2  z 2  1 , 0  x  2
Representa cada uno de los sólidos anteriores y una muestra del campo actuando sobre ellos.
Para la representación puedes utilizar las ecuaciones paramétricas siguientes:
a) Semiesfera:

x  3 sen  cos , y  3 sen  sen , z  3 cos , 0    , 0    2
2
b) Paraboloide: x  r cos , y  r sen , z  4  r 2 , 0  r  2, 0    2
c) Cilindro: x  u, y  cos , z  sen , 0  u  2, 0    2
Solución a)
 
Divergencia del campo vectorial: div F = div z i  x j  y k  0 por lo tanto,
 F  ndS =  div FdV = 0


H H
Código para dibujar el sólido:
%Dibujo de la semiesfera
fi=linspace(0,pi/2,10);t=linspace(0,2*pi,20);
[Fi,T]=meshgrid(fi,t);
x=3*sin(Fi).*cos(T);y=3*sin(Fi).*sin(T);z=3*cos(Fi);
surf(x,y,z,'EdgeColor','none', 'FaceAlpha', 0.3)
axis equal
hold on
%dibujo de la base
r=linspace(0,3,5);
[T R]=meshgrid(t,r);
x1=R.*cos(T);y1=R.*sin(T);z1=zeros(size(x1));
surf(x1,y1,z1,'EdgeColor','none', 'FaceAlpha', 0.3)
%dibujo de una muestra del campo sobre las dos superficies

quiver3(x,y,z,z,x,y);
quiver3(x1,y1,z1,z1,x1,y1)
hold off
Comprobar que se obtendría el mismo valor calculando la integral de superficie.
Solución b)
Divergencia del campo vectorial: div x 2 i  y 2 j  z 2 k  2(x  y  z )
Ecuación del paraboloide en cilíndricas: z  4  r 2 . Dominio de integración en coordenadas

cilíndricas:

57
D r, , z  / 0  r  2, 0    2, 0  z  4  r 2 
Por lo tanto,
2 2 4r 2
 F  ndS =  div FdV  2  

 (r cos   r sen   z )rdzdrd 
H H 0 0 0
Las integrales de r cos  y r sen  son nulas porque los son las de sen  y cos  en 0, 2  ,
por lo que solo queda la integral siguiente,
2 2 4r 2
 F  ndS =  div FdV  2  

 zrdzdrd  
H H 0 0 0
2 2 2
  r (4  r 2 )2 drd   2  r (4  r 2 )2 dr 
0 0 0
2
(4  r 2 )3 
2
   2r (4  r ) dr   3  

2 2
0 0
64

3
Esta integral triple también se podía haber resuelto con Matlab.
Código para dibujar el sólido:
%Dibujo del paraboloide
x=R.*cos(T);y=R.*sin(T);z=4-R.^2;
surf(x,y,z,'EdgeColor','none', 'FaceAlpha', 0.3)
axis equal
hold on
%dibujo de la base
z1=zeros(size(x));
surf(x,y,z1,'EdgeColor','none', 'FaceAlpha', 0.3)
%dibujo de una muestra del campo sobre las dos superficies

quiver3(x,y,z,x.^2,y.^2,z.^2);
quiver3(x,y,z1,x.^2,y.^2,z1.^2)
hold off
Comprobar que se obtendría el mismo valor calculando la integral de superficie.
Solución c)

Divergencia del campo vectorial: div (x  z 2 )i  (y  z 2 ) j  xk  2 
Por lo tanto,
 F  ndS =  div FdV = 2 dV  2(volumen

 de H )  4
H H H
Al resultado anterior se llega viendo que el sólido H es un cilindro de radio 1 y altura 2, cuyo
volumen es 2 . Si no nos hubiéramos dado cuenta de esto, la integral triple se puede plantear
en cilíndricas y resolver con Matlab.

Puesto que el eje del cilindro es el OX y se proyecta sobre la circunferencia y 2  z 2  1 , el

cambio de variables a coordenadas cilíndricas, en este caso, es el siguiente:
x  x , y  r cos , z  r sen 
y el dominio es, D  x, r,   / 0    2, 0  r  1, 

0  x  2 , por lo tanto la integral
triple en cilíndricas es,
2 1 2
2  dV  2    rdxdrd   4
H 0 0 0
%Dibujo del cilindro

x=linspace(0,2,10);t=linspace(0,2*pi,20);
[X,T]=meshgrid(x,t);
y=cos(T);z=sin(T);
surf(X,y,z,'EdgeColor','none', 'FaceAlpha', 0.3)
axis equal
hold on
%dibujo de las bases

r=linspace(0,1,5);
[T1 R]=meshgrid(t,r);
x1=zeros(size(T1));y1=R.*cos(T1);z1=R.*sin(T1);
x2=2*ones(size(T1));
TEOREMA DE STOKES
33
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).
Nota: Se deberá comprobar que se cumplen las hipótesis del Teorema de Stokes
Solución

59
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:
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
34
Sea S la parte del semicilindro z  1  x 2 comprendido entre y  1 , e y  1 .

Consideramos el campo F x , y, z   y  z , x 2  z 2 , y . 
a) Dibujar la superficie en el rectángulo [-1,1]x[-1,1]
b) Comprobar que se verifica el Teorema de Stokes

Solución
syms x y
fsurf(x,y,sqrt(1-x^2),[-1 1 -1 1])
xlabel('Eje X')
ylabel('Eje Y')
0.8
0.6
0.4
0.2
0
1
0.5
0
1
0.5
-0.5
0
-0.5
Eje Y -1
-1
Eje X
Se cumple  F  d r   rot F  dS  4 . Escribimos el códio matab para comprobarlo.

S S
35

Calcula el flujo saliente del rotacional del campo F x , y, z   x z , y z , x y a través de la 
superficie formada por la parte de la esfera x 2  y 2  z 2  4 dentro del cilindro x 2  y 2  1
para z  0 .
Solución
Se trata de calcular el flujo del rotacional del campo F a través de la superficie S formada por la
parte de la esfera x 2  y 2  z 2  4 que se encuentra dentro del cilindro x 2  y 2  1 .
Se aplicará el teorema de Stokes, donde la curva C frontera de S es la intersección de la esfera y

el cilindro. La curva es

61

x 2  y 2  z 2  4 z 2  3  x  cos t 

   
 2   y  sen t t   0, 2  
x 2  y2  1  x  y 2  1 
  z  3 


r t   cos t, sen t, 3  t  0, 2 
Como se pide el flujo saliente, la curva C se debe recorrer en sentido antihorario. Se tiene
entonces,
flujo   rot F  ndS   F r t   r ' t dt

S C
Como

r' t   sent, cos t, 3 

F r t     3 cos t, 3sen t, 0 
 
F r t   r' t    
3 cos t, 3sen t, 0  sen t, cos t, 0  0 
Por lo tanto, el flujo es 0.
36
Comprobar 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
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
 D
4, 6, 3  2x, 2y,1dxdy    8r cos   12r sen   3 rdrd   27
0 0
37
 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

63
i j k
  
Como rot F x , y, z     i - j  2y j N  2, 2,1
x y z
y 2 z x
se tiene
 y dx  zdy  xdz   rot F  n dS   4  2y dA

2
I 
C S T
3 3x 3 y 3x 3 y 3x
  4  2y dydx   4y  y 

 y 4  y 
2
I  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
38
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á:
 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    
 3 2 r 0 3  3
0 0
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  

65
 /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
 F  Tds   2t 3dt  8

C3 0
Luego,
32  56
 F  Tds  3
 8 
4 3

C
39
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);
U=0*ones(size(X));V=0*ones(size(X));W=6*ones(size(X));
% 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;
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
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.

67
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:
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
40
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));
% 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);
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
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:

69
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   k   x '2 y '2 '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,
 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;
U=0*ones(size(X));V=0*ones(size
(X));W=ones(size(X));
% 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
41
Se consideran:
• Las superficies
S1  z  x 2  y 2 S2  z  2 S3  z  1
• H el sólido encerrado por las tres superficies del punto anterior.

• C la curva del espacio
r t   cos t i  sen t j  t 2 k t  0,  
• La función temperatura T que en cada punto es proporcional a la distancia al plano z=0.

• El campo F x , y, z   x i  y j  z k
Se pide, calcular:
a) El trabajo realizado por F a lo largo de la curva C.
b) El área superficial de H.
c) El área de la cortina vertical determinada por la curva C y su proyección sobre el
plano z  0 .
d) El valor promedio de T sobre H.
Solución a)
Se tiene que
b
W   F  d r   F r t   r ' t dt
C a
 F  d r   cos t, sent, t   sent, cos t,2t dt

2
W 
C 0
 
4
  
3
  cos t sen t  cos t sent  2t dt  2t 3dt 
0 0
2
Otra forma:
x 2 y2 z 2
El campo es conservativo siendo su función potencial  f x , y, z     por lo que
2 2 2
el valor de la integral se podría haber calculado como:
 F  d r  f r   f r 0  f cos , sen ,    f cos 0, sen 0, 0 

2
W 
C
1 4 1 4
   
2 2 2 2
Solución b)
H es el interior de un tronco de cono.
− El área superficial correspondiente a z=2 es un círculo de radio 2, su área es 4 .

− El área superficial correspondiente a z=1 es un círculo de radio 1, su área es 

71
− El área superficial de la parte del cono comprendida entre z=1 y z=2 es
   f 
2 2
 dS   1  fx' '
A3  y
dx dy
S D
donde se ha considerado que D  x, y  / 1  x 2


 y 2  4 es la corona circular de
radios 1 y 2
 
y que f x , y  x 2  y 2 definida sobre D representa la superficie S.
2 2
   
 x   y  dx dy 
A3   1
 x 2  y 2 
 
 x 2  y 2 
D    
  1  1 dx dy  2 área(D )  2 4     3 2
D
Por lo tanto, el área superficial de H es
AH  4    3 2  5  3 2   
Solución c)
Considerando C 1 la proyección de C sobre el plano z=0, el área de la cortina vertical es
f x , y ds   z t  x ' t   y ' t  dt

2 2
Área  
C1 0
 
3
sent   cos t  dt   t
2 2

2 2
 t dt 
0 0
3
Solución d)
El valor medio de la temperatura sobre H es:

 T x, y, z dV  kz dV

Tm  H
 H
 dV  dV
H H
El conjunto H se puede expresar en coordenadas cilíndricas como H  H 1  H 2 donde
H 1  r , , z  1  r  2 , 0    2 , r  z  2
H 2  r , , z  0  r  1 , 0    2 , 1  z  2
La temperatura total sobre H es

2 2 2 1 2 2
 kz dV   kz dV   kz dV   kzr dz d  dr     kzr dz d  dr

H H1 H2 1 0 r 0 0 1
z 2
2 2 2 2
 z 2  2
 2
  rdr  2k 2  r  rdr  9k 
 kz r dz d  dr  2k   2 
 z r
 


2  4
1 0 r 1 1
1 2 2
1  4 1  3k
 k zr dz d  dr  2k   0    
2  2 2  2
0 0 1

73
9k  3k 15k
Ttotal   
4 2 4
El volumen de H es
2 2 2 1 2 2
V H    dV   dV   dV   dz d  dr    dz d  dr

H H1 H2 1 0 r 0 0 1
2 2 2 2
4
   r dz d  dr  2  2  r rdr  3
1 0 r 1
1 2 2

1
1 1
 r dz d  dr  2  2   rdr  2  
0 0 1 0
 2  2
4 7
V H    dV   
H
3 3
Finalmente, la temperatura media es
15k
Ttotal 45k
Tmedia   4 
Vol H  7 28
3
Material de consulta
Libro Interactivo. Parte II Capítulo 7.
Laboratorios y Unidades didácticas con videos explicativos: enlace
Autora: Elena E. Álvarez
Cálculo III. Gilbert Strang y Edwin Herman. Publicado bajo licencia CC BY-NC-SA 4.0 2
https://openstax.org/details/books/calculus-volume-3
2
Las imágenes con marco color rojo están tomadas de este libro.

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

75
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

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


77
ELIPSOIDE
Cartesianas
x2 y2 z2
+ + 1
=
a2 b2 c2
Paramétricas
x  a sen u cos v x  5 sen u cos v


y  3 sen u sen v
y  b sen u sen v 

z  c cos u z  2 cos u

0  u  , 0  v  2
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

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



79
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

Onocimientos Revios: Ntegración de Uperficie

Cargado por

Copyright:

Formatos disponibles

Onocimientos Revios: Ntegración de Uperficie

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Onocimientos Revios: Ntegración de Uperficie

Cargado por

Copyright:

Formatos disponibles

Tema 3 Grado en Ingeniería Mecánica

Los objetivos específicos de este tema son:

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

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

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

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

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

2 Plano tangente a una superficie en un punto

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

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

Caso 1. Superficie en explícitas. Si la superficie viene dada por la ecuación z  f (x , y ) , con (x , y )  D

• 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 )

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

Herramienta visualización plano tangente

• Plano tangente en P es:

• 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

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

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

ÁREA DE UNA SUPERFICIE

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

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

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

La herramienta siguiente permite mostrar la relación entre el área un paralelogramo y su proyección

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 C 1 en D

Un vector normal unitario a la superficie en

Vector normal al plano z = 0 en P:

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

Ángulo que forman n y k en P: γ

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

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

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

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

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:

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

Herramienta visualización área de una superficie:

INTEGRAL DE SUPERFICIE DE UN CAMPO ESCALAR

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

• S es una superficie suave 1.

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

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

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

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

Suponiendo S definida en paramétricas

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

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

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

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

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

INTEGRAL DE SUPERFICIE DE UN CAMPO VECTORIAL O INTEGRAL DE FLUJO

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

Herramienta interpretación integral de superficie de un campo vectorial

Prof. Elena E. Alvarez Saiz Universidad de Cantabria

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