Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...
Capítulo 1 Métodos de residuos ponderados Funciones de prueba ... Capítulo 1 Métodos de residuos ponderados Funciones de prueba ...
2. Las condiciones de borde no pueden construirse en forma arbitraria, sino que deben ser compatibles con la ecuación diferencial que gobierna el problema. En este sentido la formulación variacional posee la ventaja de que cualquier condición de borde que pueda incorporarse naturalmente al problema, como consecuencia de la integración por partes, será automaticamente compatible con la ecuación diferencial. 3. Recordemos que las condiciones de borde entran en la definición del problema variacional (3.18) en dos formas: las condiciones de borde esenciales, que son la especificación del valor de la solución, entran en el problema a partir de la definición del espacio admisible de las funciones del problema, mientras que las condiciones naturales, que conllevan la especificación de las derivadas de la solución, determinan la forma de la ecuación variacional. Ejemplos: i) −u(x), xx +u(x) = f(x) 0 < x < L u(0) = 0, u(L) = 0 ii) −u(x), xx +u(x) = f(x) 0 < x < L u(0), x = γ 0 , Los enunciados variacionales de estos problemas son u(L), x = γ L V-i) Encontrar u en H0 1 tal que (es decir que u satisface las condiciones (3.19) y además u(0) = u(L) = 0) L∫ L∫ (W ,x u ,x + W u) dx = W fdx, ∀ W ∈ H0 1 0 0 V-ii) Encontrar u ∈ H 1 tal que L∫ (W ,x u ,x + W u) dx = L∫ W fdx − γ 0 W (0) + γ L W (L), ∀ W ∈ H 1 0 0 Estos ejemplos muestran que las condiciones esenciales se introducen en V-i) a partir de la definición de la clase de funciones admisibles para proponer las funciones de prueba y las naturales como datos en el segundo miembro de la ecuación en V-ii). 4. Por último, en analogía con la energía U definida anteriormente, podemos escribir la norma energética para el problema (3.18) (con b(x) > 0) [∫ L ‖e‖ E = También puede utilizarse la norma H 1 equivalente 0 [∫ L ‖e‖ 1 = 0 ( k e 2 ,x + b e 2) ] 1/2 dx (3.20) ( e 2 ,x + e 2) ] 1/2 dx (3.21) Ejercicio N ◦ 26: Construir la forma débil indicando en cada caso el espacio de funciones de prueba admisibles. 44
a) b) c) − (k(x) u (x) , x ) , x +u(x), x +u (x) = 0 en 0 < x < 1; 1 < x < 2; 2 < x < 3; 3 < x < 4; con [k (1) u ,x (1)] = [k (2) u ,x (2)] = 0; [k (3) u ,x (3)] = 10 y ⎧ ⎨ k(x) ⎩ 1 0 ≤ x < 1 2 1 ≤ x < 2 1 2 < x < 4 ; ⎧ ⎨ ⎩ u (0) = 0 u (4) , x = 3 −u (x) ,xx + u (x) = δ (x − 1) 0 < x < 2 u (0) ,x = 2, u (2) ,x + u (2) = 3 u (x) ,xx + u (x) = 0 0 < x < 1; 1 < x < 2 [u ,x (1)] = 1; u (0) = u (2) = 0 Ejercicio N ◦ 27: los problemas de valores de contorno siguientes están mal condicionados o no satisfacen las hipótesis descriptas en las notas. ?‘Qué está mal en ellos? a) 1 (x − 1) u (x) ,xx + u (x) = sen (x) 0 < x < 2 b) c) d) y u (0) = u (2) = 0 (1 − x 2 ) u (x) ,xx + u (x) ,x = 3 0 < x < 2 y u (0) = 0 ; u (2) ,x = 1 u (x) ,xx + x u (x) ,x = 3 0 < x < 1 y u (0) ,xx = 0 ; u (1) = 0 −u (x) ,xx + u (x) ,x = x 0 < x < 1 y 1 2 u (0) ,x + u (1) ,x = u (1) ; u (0) ,x + 2 u (1) ,x = u (0) 3.4. Funciones de aproximación basadas en la interpolación de Lagrange En esta sección extenderemos las ideas respecto a las funciones de interpolación elementales introduciendo una técnica que conduce a lo que se conoce como elementos finitos Lagrangeanos, nombre que deriva del conceptos de interpolación utilizando polinomios de Lagrange: 1. Consideremos un elemento típico Ω e , aislado de la malla, y establezcamos un sistema local de coordenadas ξ, cuyo origen se ubica en el centro del elemento, escalado de tal forma de satisfacer que ξ = −1 en el extremo izquierdo y ξ = 1 en el derecho, Fig.2.a. Esto se realiza mediante la transformación general que para el caso de un elemento de k + 1 nodos se escribe ξ(x) = 2 x − xe 1 − x e k+1 h e (3.22) 45
- Page 1 and 2: Capítulo 1 Métodos de residuos po
- Page 3 and 4: ( ∂ k ∂u ) + ∂ ( k ∂u ) + F
- Page 5 and 6: δ (x − x l ) = 0 para x ≠ x l
- Page 7 and 8: p = F (1.29) en la que k y F no dep
- Page 9 and 10: La solución de los sistemas conduc
- Page 11 and 12: Una viga de longitud 1 está simple
- Page 13 and 14: Notar que K es no simétrica aun cu
- Page 15 and 16: integrando por partes el primer té
- Page 17 and 18: La sustitución de la aproximación
- Page 19 and 20: Figura 2 Ejercicio 21. a) Problema
- Page 21 and 22: Capítulo 2 El método de elementos
- Page 23 and 24: Figura 1 Aproximación de una funci
- Page 25 and 26: Figura 2 Comportamiento de tres fun
- Page 27 and 28: 2.4.1. Propiedades de la matriz K y
- Page 29 and 30: Figura 5 Descripción local y globa
- Page 31 and 32: 2.4.3.2. Cálculo del vector de car
- Page 33 and 34: de funciones. Una medida universalm
- Page 35: Figura 6 Curvas del error utilizand
- Page 38 and 39: eferiremos como ” ley de conserva
- Page 40 and 41: ambos miembros haciendo tender el p
- Page 42 and 43: 1. Condiciones de borde esenciales,
- Page 46 and 47: de tal modo que todo punto x en el
- Page 48 and 49: N 1 (ξ) = (ξ − ξ 2) (ξ 1 −
- Page 50 and 51: modo de ubicar un nodo en todo punt
- Page 52 and 53: Figura 5 Malla de N nudos y N-1 ele
- Page 54 and 55: 3. Condiciones naturales de Neumann
- Page 56 and 57: (c) Considere las condiciones (ii)
- Page 58 and 59: los polinomios de interpolación in
- Page 60 and 61: La derivada segunda de u respecto a
- Page 62 and 63: 4.3.3. Formulación a partir del Pr
- Page 64 and 65: ¯K = ∫ L 0 BT D B d¯x 1 = ¯K =
- Page 66 and 67: Las cargas normales al eje de la vi
- Page 68 and 69: 2. θ 1 (ξ) = ∑ 2 I=1 N I (ξ)
- Page 70 and 71: 4.5.7. Cambio de base La expresión
- Page 72 and 73: Figura 5 Viga con deformación de c
- Page 74 and 75: 4.6.1. Matriz de rigidez de una vig
- Page 76 and 77: Donde las funciones de forma son: N
- Page 78 and 79: (una para cada W J ) en función de
- Page 80 and 81: Ejercicio 1-Sea el problema de conv
- Page 82 and 83: o directamente las coordenadas noda
- Page 84 and 85: donde r es el residuo que se quiere
- Page 86: ⎡ K 3−4 = ⎢ ⎣ K 4−5 = ⎡
- Page 89 and 90: Figura 1 Conducción del calor en 2
- Page 91 and 92: 5.3. Forma variacional del problema
- Page 93 and 94: 5.5. Flujo en un medio poroso El fl
a)<br />
b)<br />
c)<br />
− (k(x) u (x) , x ) , x +u(x), x +u (x) = 0<br />
en<br />
0 < x < 1; 1 < x < 2; 2 < x < 3; 3 < x < 4;<br />
con<br />
[k (1) u ,x (1)] = [k (2) u ,x (2)] = 0; [k (3) u ,x (3)] = 10<br />
y<br />
⎧<br />
⎨<br />
k(x)<br />
⎩<br />
1 0 ≤ x < 1<br />
2 1 ≤ x < 2<br />
1 2 < x < 4<br />
;<br />
⎧<br />
⎨<br />
⎩<br />
u (0) = 0<br />
u (4) , x = 3<br />
−u (x) ,xx<br />
+ u (x) = δ (x − 1) 0 < x < 2<br />
u (0) ,x<br />
= 2, u (2) ,x<br />
+ u (2) = 3<br />
u (x) ,xx<br />
+ u (x) = 0 0 < x < 1; 1 < x < 2<br />
[u ,x (1)] = 1; u (0) = u (2) = 0<br />
Ejercicio N ◦ 27: los problemas <strong>de</strong> valores <strong>de</strong> contorno siguientes están mal condicionados o no<br />
satisfacen las hipótesis <strong>de</strong>scriptas en las notas. ?‘Qué está mal en ellos?<br />
a)<br />
1<br />
(x − 1) u (x) ,xx<br />
+ u (x) = sen (x) 0 < x < 2<br />
b)<br />
c)<br />
d)<br />
y u (0) = u (2) = 0<br />
(1 − x 2 ) u (x) ,xx<br />
+ u (x) ,x<br />
= 3 0 < x < 2<br />
y u (0) = 0 ; u (2) ,x<br />
= 1<br />
u (x) ,xx<br />
+ x u (x) ,x<br />
= 3 0 < x < 1<br />
y u (0) ,xx<br />
= 0 ; u (1) = 0<br />
−u (x) ,xx<br />
+ u (x) ,x<br />
= x 0 < x < 1<br />
y<br />
1<br />
2 u (0) ,x + u (1) ,x = u (1) ; u (0) ,x + 2 u (1) ,x = u (0)<br />
3.4. <strong>Funciones</strong> <strong>de</strong> aproximación basadas en la interpolación <strong>de</strong><br />
Lagrange<br />
En esta sección exten<strong>de</strong>remos las i<strong>de</strong>as respecto a las funciones <strong>de</strong> interpolación elementales<br />
introduciendo una técnica que conduce a lo que se conoce como elementos finitos Lagrangeanos,<br />
nombre que <strong>de</strong>riva <strong>de</strong>l conceptos <strong>de</strong> interpolación utilizando polinomios <strong>de</strong> Lagrange:<br />
1. Consi<strong>de</strong>remos un elemento típico Ω e , aislado <strong>de</strong> la malla, y establezcamos un sistema local <strong>de</strong><br />
coor<strong>de</strong>nadas ξ, cuyo origen se ubica en el centro <strong>de</strong>l elemento, escalado <strong>de</strong> tal forma <strong>de</strong> satisfacer<br />
que ξ = −1 en el extremo izquierdo y ξ = 1 en el <strong>de</strong>recho, Fig.2.a. Esto se realiza mediante la<br />
transformación general que para el caso <strong>de</strong> un elemento <strong>de</strong> k + 1 nodos se escribe<br />
ξ(x) = 2 x − xe 1 − x e k+1<br />
h e (3.22)<br />
45