Asignaturas:Cálculo Avanzado – Civil

Análisis Numérico y Cálculo Avanzado –Industrial

PROBLEMAS PARCIAL 3: METODO DE DIFERENCIAS FINITAS (MDF) y MÉTODO DE ELEMENTOS FINITOS (MEF).

1. Sea el siguiente problema con valores en el borde:

∂u ∂ 2 u
E.D: = 0≤ x ≤ 1 , t > 0
∂t ∂x 2
∂u
C.B: u (0, t) =0 (1,t) = 0 t>0
∂x
π
C.I: u(x, 0) =sen x 0≤ x ≤ 1
2
a. Describir “completamente” una situación física modelada por este problema.
b. Deducir un método explícito de diferencias finitas para resolverlo numéricamente.
c. Implementar el método formulado en el inciso anterior con Δt = 0.1 y Δx = 0.5 para determinar un
valor aproximado de u (½ , 0.3).
d. ¿Es posible estimar el error cometido en la aproximación, aplicando el Teorema de estabilidad y
convergencia para la ecuación del calor con temperaturas impuestas en los bordes? Por qué?

2. Una varilla homogénea aislada lateralmente, de longitud 1, tiene temperatura inicial 0. Su borde izquierdo se
mantiene a temperatura 0, mientras que la temperatura en el borde derecho varía senoidalmente según
u(1,t) = sen (25/3)π t.
a. Escribir el modelo matemático que describe la evolución de la temperatura en la varilla, suponiendo
que el coeficiente de difusividad k es igual a 1.
b. ¿Es posible encontrar una expresión para la solución exacta del problema aplicando directamente el
método de separación de variables? Justificar.
c. Formular un esquema de diferencias finitas para el problema continuo planteado en (a).
k ∆t
d. Implementar el esquema formulado en el inciso anterior, con Δx = 0.2 y s= = 0.5 para
( ∆x ) 2
determinar una aproximación de la temperatura u(0.8, 0.06) Y Estimar el error cometido.

3. Considerar el siguiente problema de conducción del calor en una barra homogénea de longitud cuatro.

∂u ∂ 2 u
E.D: = , 0< x <4 , t > 0
∂t ∂x 2

C.B: u(0, t) = 4t , u (4,t) = 4t , t ≥ 0

C.I: u(x,0) = 2 x ( x – 4 ) , 0≤ x ≤4

a. Verificar que tiene solución exacta u(x, t) = 2 x ( x – 4 ) + 4t.

b. Usar el método de diferencias finitas con ∆x = 2 y ∆t = 1 para aproximar u(2,1).
c. Comparar con la solución exacta ¿Cuál es el error cometido? Comentar.

Página 1 de 4
 Asignaturas:Cálculo Avanzado – Civil
Análisis Numérico y Cálculo Avanzado –Industrial

4. El siguiente problema modela la distribución de temperaturas en un alambre circular aislado de longitud 2 π

con distribución inicial de temperatura conocida,

= −π < <π, >0

−π, = π, −π, = π, , >0
, 0 = 1 − cos 2 −π ≤ x ≤ π

Aplicar la técnica del Método de las Diferencias Finitas para resolver numéricamente este problema,
siguiendo los siguientes pasos:

a. Dibujar la malla de puntos para -π ≤ x ≤ π y 0 ≤ t ≤ 1, usando Δx = ⁄4 y Δt = 1⁄4, indicando los

puntos correspondientes a condiciones de borde y condición inicial. ¿Cuántos nodos tiene esta malla en
total?
b. Deducir un esquema de diferencias finitas y formular completamente el problema discreto (especificar el
rango de valores para los sub/supraíndices en cada ecuación).
c. Implementar el esquema deducido en el inciso anterior para calcular el valor aproximado de la
temperatura en el punto (π, 1⁄4 .

5. Aplicar la técnica del método de diferencias finitas para formular el problema discreto que modela las
vibraciones transversales de una cuerda elástica uniforme, la cual (al iniciar) se estira hasta la longitud L = 1 y
se fija en los puntos extremos, es decir u(x,0)= 0 para 0 ≤ x ≤ 1 y u(0,t) = u(1,t) = 0 para t > 0. Suponer que

la velocidad c de propagación de las ondas por la cuerda es igual a 1 y que, la velocidad inicial ut (x,0) está

dada por sen( π x).

Sugerencia: Para la EDP usar diferencias centradas de segundo orden en espacio y tiempo y diferencia
adelantada para la velocidad inicial.
a. Implementar el esquema anterior para aproximar u( 0.2, 0.2) tomando Δx = 0.2 y Δt = 0.1. Notar
∆x
que con esta malla el método será estable pues se satisface la restricción c ≤ correspondiente a
∆t
la ecuación de ondas unidimensional.
b. Comparar el valor aproximado con el valor exacto que se obtiene al aplicar separación de variables
(Ejercicio de Guía Parcial 2).

6. Aplicar la técnica del método de diferencias finitas para formular el problema discreto correspondiente al
problema que modela el potencial electrostático, u, en una placa rectangular cargada que no contiene carga
libre (ecuación de Laplace) con potencial prescrito en cada uno de los bordes del rectángulo. Usar
∆ =∆ y %
$ ↔ ' $, %( = )∆ , *∆ para notar que el valor numérico de %
$ es el promedio
entre los valores del potencial en sus cuatro “vecinos”.

Página 2 de 4
 Asignaturas:Cálculo Avanzado – Civil
Análisis Numérico y Cálculo Avanzado –Industrial

7. Una delgada placa de acero de 10cmx10cm tiene uno de sus bordes a 100º C y los otros a 0ºC ¿Cuáles son las
temperaturas de estado estable (régimen estacionario) en los puntos interiores? Usar Δx = Δy = 10/3.

8. Considerar una viga cuyo eje neutro yace en el intervalo 0 ≤ x ≤ 1 . Una carga transversal distribuida q = q(x)
genera a lo largo de la viga un momento flector M = M(x) cuya distribución satisface la ecuación diferencial,

d 2M
= q(x). Si la viga está simplemente apoyada (o sea, M =0 en los extremos) y la carga es sen π x:
dx 2

a. calcular los momentos flectores en los puntos interiores usando el método de diferencias finitas con
Δx = 0.25.
b. Comparar con los valores exactos.

9. Un problema típico en cálculo estructural es determinar el desplazamiento u(x) de un punto material

localizado a una distancia x a lo largo del eje de una barra de longitud L cargada uniaxialmente, con
constante de rigidez k y desplazamientos prescriptos en los bordes. Si la carga distribuida por unidad de

d 2u
longitud en la dirección del eje de la barra está dada por la función f(x), la ecuación de gobierno es k =
dx 2
– f(x) para 0 ≤ x ≤ L.
a. Si la barra de longitud L=1 está fija en sus extremos (u(0) = u(1) = 0) y k = 1, escribir en forma

matricial el sistema de ecuaciones lineales que resulta al aplicar la técnica del MDF con ∆ x = 1/4,
para determinar el desplazamiento en los puntos 1/4, 1/2 y 3/4 de la barra si la carga es constante.
Comparar con los valores exactos.
b. Calcular la matriz de rigidez y el vector de cargas que surgen al aplicar el método de los elementos
finitos para resolver (a) con funciones de base lineales en una partición equiespaciada del intervalo [

0, 1] con ∆ x = 1/4. Escribir en forma matricial el sistema algebraico que resulta y comparar con (a).

10. El potencial eléctrico u(x) en un medio homogéneo, en cualquier punto entre dos placas paralelas, ubicadas

d 2u
en x = 0 y x = 1, con carga uniforme f(x) = 1 está modelado por la ecuación: – = f(x), 0 ≤ x ≤ 1 .
dx 2

Si el potencial en ambas placas es nulo, se pide:

a. Aplicar el método de los elementos finitos con funciones bases lineales, adoptando una partición
regular de 6 puntos, para hallar los valores aproximados en los cuatro puntos interiores de la malla.
b. Hallar la solución exacta del problema y comparar con los valores hallados en (a).

11. La ecuación que modela el desplazamiento u(x) en una viga de 2 m de longitud, que se encuentra empotrada
+
+
en sus extremos, es + (1 – )2 = 0.

Página 3 de 4
 Asignaturas:Cálculo Avanzado – Civil
Análisis Numérico y Cálculo Avanzado –Industrial

a. Aplicar el Método de los Elementos Finitos con funciones de base lineales y una partición espaciada
(con Δ = 2⁄3) para calcular los desplazamientos en los puntos = 2 ⁄3 y = 4⁄3 de la viga.
b. Obtener una expresión para la solución aproximada que proporciona el MEF para todo x en el
intervalo 0 ≤ ≤ 2 y determinar el valor numérico del desplazamiento en = 1⁄3.
c. Hallar la solución exacta del problema. Calcular el error local cometido en la aproximación para
= 1⁄3 y comentar
------------------------------------------------------------------------------------------------------------------------------------------------
RESPUESTAS:
(1) d) u13 = 0.195 e) No, pues no la condición de borde derecho no es de temperatura impuesta;
(2) (b) No, pues la condición de borde derecha no es una ecuación diferencial homogénea.
(d) u43 = √3/4 ;(e) cota para error = 0.5483 y 0.24 (para la malla refinada)
(3) (b) u(2,1) = - 4 ; (c) coincide con la sol exacta pues la cota del error es cero al ser las cond de borde y cond
inicial, funciones de 1er y 2do grado respectivamente.

; = 01 3 1
, * = 1,2, 3, 4. ; CI: = 1 − cos 2)∆ para ) = −4, −3, … . , 3, 4.
2 2
% % % 9
-. = . . 4 $
c) >. = 1.
(4) a) 45; b) CB:

(5) b) 0.112 ; (c) u(x,t) = 1/π sen(πx).sen(π.t) ; error local = 2x10-3

(6) ?@A: =. + + +
% > % % %-> %3>
$ $-> $3> $ $ ; j=1,2,…,N-1 y m= 1,2,…,M-1

CD: 9
$ = A> )∆ , E
$ = A4 )∆ , ) = 0, 1, 2, … , F y %
9 = G> *∆ , %
H = G4 *∆ , * = 1,2, … , I − 1

>; >; = 25/2; 25/2; 75/2; 75/2) con A4 = 100, A> = G> = G4 = 0.
> > 4 4
(7) ( 4; 4

(8) Val aprox = -0.0754; -0.1067; -0.0754

− 2 1 0  u1   f
    
(9) a) y b) ambos sistemas son equivalentes,  1 − 2 1  u 2  = −(1 / 4) 
2
f
 0 1 − 2  u 3   f 
 
(10) a) 0.08, 0.12, 0.12, 0.08 b) u(x) = -x2/2 + x/2

0,12345 , 0≤ ≤M
4

K
(11) a) (0.0823; 0.0823) b) J = 0.0823 , ≤ ≤ = + − +
4 . -> .
1

M M >4 M 4 M
c) ,
K −0,12345 + 0.2469 , .
≤ ≤2
M

Error local en 1/3 = 0.0257 y 38% el error relativo. Por lo cual, hay que refinar la malla para mejorar la aprox.
------------------------------------------------------------------------------------------------------------------------------------------------------
Selección actualizada por: María Cecilia Lazaronni y Adriana Frausin (2018)
BIBLIOGRAFIA:
• Haberman R. “Applied Partial Differential Equations” - Prentice Hall - 5° Edición - 2012
• Bleecker D., Csordas G. “Basic Partial Differential Equations” - Ed. International Press – 1997

Página 4 de 4