MATEMÁTICAS AVANZADAS

ANÁLISIS POR ELEMENTOS FINITOS

ANGIE PAOLA MENDOZA CUY

UNIVERSIDAD INDUSTRIAL DE SANTANDER

FACULTAD DE INGENIERÍAS FISICOMECÁNICAS
ESCUELA DE INGENIERÍA CIVIL
BUCARAMANGA
2019
Actividad:

Resolver la ecuación diferencial mediante un análisis por elementos finitos con formulación

integral débil.

𝑈 ′′ − 𝑈 = 1

1. Calculo de la solución exacta

- Solución particular de la ecuación diferencial no homogénea

𝑈𝑝 (𝑥) = −1

- Solución de la homogénea asociada

𝑈 ′′ − 𝑈 = 0

𝑆𝑖 𝑈 = 𝑒 𝑛𝑥 → 𝑈 ′ = 𝑛 ∗ 𝑒 𝑛𝑥 → 𝑈 ′′ = 𝑛2 ∗ 𝑒 𝑛𝑥

Reemplazando:

𝑛2 ∗ 𝑒 𝑛𝑥 − 𝑒 𝑛𝑥 = 0 → 𝑛 = ±1

La solución de esta ecuación diferencial lineal homogénea de segundo orden es de la

forma:

𝑈 = 𝑐1 𝑒 𝑛𝑥 + 𝑐2 𝑒 𝑛𝑥

Entonces la solución de la homogénea es:

𝑈ℎ (𝑥) = 𝑐1 𝑒 𝑥 + 𝑐2 𝑒 −𝑥

Para verificar la existencia de las constantes c1 y c2, se verificará que las soluciones sean

un conjunto fundamental. Para esto se calcula el Wronskiano

𝑒𝑥 𝑒 −𝑥
𝑊=| 𝑥 | = −𝑒 𝑥 𝑒 −𝑥 − 𝑒 𝑥 𝑒 −𝑥 = −2𝑒 𝑥 𝑒 −𝑥
𝑒 −𝑒 𝑥

2
 Como W es diferente de cero, se puede decir que sí es un conjunto fundamental de

soluciones y que son linealmente independientes, entonces existen valores para c1 y c2

diferentes de cero.

- Solución general

𝑈(𝑥) = 𝑈ℎ (𝑥) + 𝑈𝑝 (𝑥)

𝑈(𝑥) = 𝑐1 𝑒 𝑥 + 𝑐2 𝑒 −𝑥 − 1

Con las condiciones iniciales:

𝑈(0) = 𝑐1 + 𝑐2 − 1 = 1 → 𝑐1 + 𝑐2 = 2 (1)

𝑈′(1) = 𝑐1 𝑒 − 𝑐2 𝑒 −1 − 1 = 0 → 𝑐1 𝑒 − 𝑐2 𝑒 −1 = 1 (2)

Resolviendo el sistema de dos ecuaciones:

2 2𝑒 2
𝑐1 = 2
; 𝑐2 =
(𝑒 + 1) (𝑒 2 + 1)

2𝑒 𝑥 2𝑒 2−𝑥
𝑈𝑒𝑥𝑎 (𝑥) = 2 + −1
(𝑒 + 1) (𝑒 2 + 1)

2. Análisis por elementos finitos

- Función de aproximación
2
̅ = ∅0 (𝑥) + ∑ 𝑐𝑗 ∗ ∅𝑗 (𝑥)
𝑈
𝑗=1

∅0 (𝑥) debe ser una función que cumple con las condiciones iniciales U(0)=1 y U’(1)=0.

∅0 (𝑥) = 1

3
 Tomando ∅1 (𝑥) = 𝑥 2 − 2𝑥 y ∅2 (𝑥) = 𝑥 3 − 3𝑥 se debe verificar que evaluadas en los

puntos de condiciones iniciales den cero.

∅1 (0) = 02 − 2 ∗ 0 = 0 cumple

∅1′ (1) = 2 ∗ 1 − 2 ∗ 1 = 0 cumple

∅2 (0) = 03 − 3 ∗ 0 = 0 cumple

∅′2 (1) = 3 ∗ 12 − 3 ∗ 1 = 0 cumple

̅ = 1 + 𝑐1 (𝑥 2 − 2𝑥) + 𝑐2 (𝑥 3 − 3𝑥)
𝑈

- Formulación Weighted Integral

̅ ′′ − 𝑈
Residuo → 𝑅 = 𝑈 ̅−1=0

𝑅 = −𝑐2 𝑥 3 − 𝑐1 𝑥 2 + 2𝑐1 𝑥 + 9𝑐2 𝑥 + 2𝑐1 − 2 = 0

1
∫ 𝑅 ∗ 𝑣𝑗 (𝑥) = 0
0

Para 𝑣1 (𝑥) = 1 𝑦 𝑣2 (𝑥) = 𝑥

Ecuación 1:
1
∫ (−𝑐2 𝑥 3 − 𝑐1 𝑥 2 + 2𝑐1 𝑥 + 9𝑐2 𝑥 + 2𝑐1 − 2) ∗ 1 dx = 0
0

1
∫ (−𝑐2 𝑥 3 − 𝑐1 𝑥 2 + 2𝑐1 𝑥 + 9𝑐2 𝑥 + 2𝑐1 − 2) 𝑑𝑥 = 0
0

-1 9 -1
( + ) 𝑐2 + ( + 1 + 2) 𝑐1 − 2 = 0
4 2 3

17 8
𝑐2 + 𝑐1 − 2 = 0 (1)
4 3

4
 Ecuación 2:
1
∫ (−𝑐2 𝑥 3 − 𝑐1 𝑥 2 + 2𝑐1 𝑥 + 9𝑐2 𝑥 + 2𝑐1 − 2) ∗ x dx = 0
0

1
∫ (−𝑐2 𝑥 4 − 𝑐1 𝑥 3 + 2𝑐1 𝑥 2 + 9𝑐2 𝑥 2 + 2𝑐1 𝑥 − 2𝑥) 𝑑𝑥 = 0
0

-1 -1 2
( + 3) 𝑐2 + ( + + 1) 𝑐1 − 1 = 0
5 4 3

14 17
𝑐2 + 𝑐 −1=0 (2)
5 12 1

Resolviendo el sistema de dos ecuaciones con dos incógnitas:

324 -40
𝑐1 = ; 𝑐2 =
347 347

324 2 40 3
̅ =1+
𝑈 (𝑥 − 2𝑥) − (𝑥 − 3𝑥)
347 347

xi Uexa Uaprox 1 Error 1

0 1.0000 1.0000 0.00000 Análisis Elementos Finitos
0.1 0.8574 0.8571 0.14213 1.2000
0.2 0.7335 0.7321 0.26469 1.0000
Exacta
0.3 0.6268 0.6244 0.36934 Aproximada 1
0.8000
0.4 0.5365 0.5334 0.45769
U (x)

0.5 0.4615 0.4582 0.53130 0.6000

0.6 0.4012 0.3983 0.59154 0.4000
0.7 0.3549 0.3529 0.63944
0.2000
0.8 0.3221 0.3213 0.67525
0.9 0.3026 0.3028 0.69816 0.0000
0 0.5 1 1.5
1 0.2961 0.2968 0.70633 x

Para 𝑣1 (𝑥) = 1 𝑦 𝑣2 (𝑥) = 𝑥 + 1

324 -40 324 2 40 3

𝑐1 = ; 𝑐2 = ̅ =1+
→ 𝑈 (𝑥 − 2𝑥) − (𝑥 − 3𝑥)
347 347 347 347

5
 Para 𝑣1 (𝑥) = 𝑥 𝑦 𝑣2 (𝑥) = 𝑥 2

780 -80
𝑐1 = ; 𝑐2 =
881 881

780 2 80 3
̅ =1+
𝑈 (𝑥 − 2𝑥) − (𝑥 − 3𝑥)
881 881

Para 𝑣1 (𝑥) = 𝑥 𝑦 𝑣2 (𝑥) = 2𝑥

14 17
𝑐2 + 𝑐 −1=0 (1)
5 12 1

28 17
𝑐2 + 𝑐1 − 2 = 0 (2)
5 6

6
Se puede ver que las dos ecuaciones son linealmente dependientes por lo que es posible

encontrar infinitos valores de ci .

Si decimos que

12 168
𝑐2 = s ; 𝑐1 = − s
7 85

12 168
̅ = 1+(
𝑈 − s) (𝑥 2 − 2𝑥) − s(𝑥 3 − 3𝑥)
7 85

Mínimos cuadrados

1
𝜕R
∫ 𝑅 ∗ 𝑣𝑗 (𝑥) = 0 donde 𝑣𝑗 (𝑥) =
0 𝜕𝑐j

1
∫ 𝑅∗𝑅 =0
0

1
∫ (−𝑐2 𝑥 3 − 𝑐1 𝑥 2 + 2𝑐1 𝑥 + 9𝑐2 𝑥 + 2𝑐1 − 2) ∗ (−𝑐2 𝑥 3 − 𝑐1 𝑥 2 + 2𝑐1 𝑥 + 9𝑐2 𝑥 + 2𝑐1 − 2) = 0
0

36 2 721 32 824 2
→ 𝑓(𝑐1, 𝑐2) = 𝑐1 + 𝑐1 𝑐2 − 𝑐1 + 𝑐 − 17𝑐2 + 4 = 0
5 30 3 35 2

Ecuación 1:

𝜕𝑓(𝑐1, 𝑐2) 72 721 32

= 𝑐1 + 𝑐2 − =0 (1)
𝜕𝑐1 5 30 3

Ecuación 2:

𝜕𝑓(𝑐1, 𝑐2) 721 1648

= 𝑐1 + 𝑐 − 17 = 0 (2)
𝜕𝑐2 30 35 2

7
Resolviendo el sistema de dos ecuaciones

5730 -72800
𝑐1 = ; 𝑐2 =
6143 632729

5730 2 72800 3
̅ = 1+
𝑈 (𝑥 − 2𝑥) − (𝑥 − 3𝑥)
6143 632729

xi Uexa Uaprox 1 Error 1

0 1.0000 1.0000 0.00000
0.1 0.8574 0.8572 0.14229
0.2 0.7335 0.7323 0.26496
0.3 0.6268 0.6247 0.36976
0.4 0.5365 0.5337 0.45830
0.5 0.4615 0.4586 0.53214
0.6 0.4012 0.3987 0.59261
0.7 0.3549 0.3533 0.64070
0.8 0.3221 0.3218 0.67690
0.9 0.3026 0.3033 0.69974
1 0.2961 0.2973 0.70792

8
Gallerkin

1
∫ 𝑅 ∗ 𝑣𝑗 (𝑥) = 0 donde 𝑣𝑗 (𝑥) = ∅𝑗 (𝑥)
0

1
∫ 𝑅 ∗ ∅𝑗 (𝑥) = 0
0

Ecuación 1:

1
∫ (−𝑐2 𝑥 3 − 𝑐1 𝑥 2 + 2𝑐1 𝑥 + 9𝑐2 𝑥 + 2𝑐1 − 2) ∗ (𝑥 2 − 2𝑥) dx = 0
0

4 211 28
− 𝑐2 − 𝑐1 = 0
3 60 15

Ecuación 2:

1
∫ (−𝑐2 𝑥 3 − 𝑐1 𝑥 2 + 2𝑐1 𝑥 + 9𝑐2 𝑥 + 2𝑐1 − 2) ∗ (𝑥 3 − 3𝑥) dx = 0
0

5 236 211
− 𝑐2 − 𝑐 =0
2 35 60 1

Resolviendo el sistema de dos ecuaciones

5010 -80
𝑐1 = ; 𝑐2 =
5537 791

5010 2 80 3
̅ = 1+
𝑈 (𝑥 − 2𝑥) − (𝑥 − 3𝑥)
5537 791

9
xi Uexa Uaprox 1 Error 1
0 1.0000 1 0.00000
0.1 0.8574 0.8583 0.14357
0.2 0.7335 0.7341 0.26741
0.3 0.6268 0.6268 0.37311
0.4 0.5365 0.5358 0.46222
0.5 0.4615 0.4604 0.53604
0.6 0.4012 0.4002 0.59635
0.7 0.3549 0.3543 0.64352
0.8 0.3221 0.3223 0.67845
0.9 0.3026 0.3036 0.70073
1 0.2961 0.2975 0.70859

10