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 ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Siguiendo ahora un or<strong>de</strong>n creciente para j primero (bucle externo <strong>de</strong> 1 a n) y para i <strong>de</strong>spués<br />
(interno <strong>de</strong> j a n) po<strong>de</strong>mos ir evaluando cada término <strong>de</strong> la factorización. Es <strong>de</strong>cir para la primera<br />
columna j = 1<br />
d 11 = k 11<br />
l i1 = k i1 /d 11<br />
para cada columna posterior j se conocen entonces todos los valores <strong>de</strong> d kk y l ik <strong>de</strong> las columnas<br />
anteriores (k < j). Po<strong>de</strong>mos escribir<br />
[<br />
]<br />
l ij = 1<br />
j−1<br />
∑<br />
k ij − g ik l jk<br />
d jj<br />
pero todavía no conocemos d jj por ello en vez <strong>de</strong> calcular l ij evaluemos primero para i <strong>de</strong> 1 a j − 1<br />
[<br />
]<br />
j−1<br />
∑<br />
g ij = l ij d jj = k ij − g ik l jk (1)<br />
y una vez obtenidos los valores <strong>de</strong> la columna completa po<strong>de</strong>mos ahora calcular<br />
k=1<br />
k=1<br />
j−1<br />
∑<br />
d jj = k jj − g jk l jk<br />
l ij = g ij /d jj<br />
y pasar a la siguiente columna hasta completar el cálculo <strong>de</strong> todos los coeficientes. Resulta muy<br />
importante notar en este algoritmo que en (1) g ij será nulo si inicialmente k ij y todos los términos<br />
previos <strong>de</strong> la fila k = 1, j − 1 son nulos, es <strong>de</strong>cir que la matriz K y L tendrán el mismo perfil.<br />
A<strong>de</strong>más una vez calculado el valor <strong>de</strong> l ij no resulta ya necesario mantener el correspondiente valor<br />
<strong>de</strong> k ij y por lo tanto pue<strong>de</strong> utilizarse las posiciones <strong>de</strong> memoria <strong>de</strong> K para guardar los <strong>de</strong> L.<br />
Finalmente los coeficientes <strong>de</strong> la matriz diagonal D pue<strong>de</strong>n guardarse en las posiciones diagonales<br />
<strong>de</strong> K ya que los elementos diagonales <strong>de</strong> L valen 1 y no requieren ser recordados.<br />
2-Substitución hacia a<strong>de</strong>lante para i <strong>de</strong> 1 a n<br />
k=1<br />
∑i−1<br />
y i = f i − l ik y k<br />
3-Substitución hacia atrás para i <strong>de</strong> n a 1<br />
k=1<br />
[<br />
]<br />
x i = 1 ∑i+1<br />
y i − l ki x k<br />
d ii<br />
Una <strong>de</strong>scripción más <strong>de</strong>tallada <strong>de</strong> la base teórica <strong>de</strong>l presente algoritmo y su programación en<br />
códigos <strong>de</strong> elementos finitos pue<strong>de</strong> verse en las referencias<br />
k=n<br />
K.J.Bathe y E.L.Wilson, Numerical methods in finite element analysis, Prentice-Hall, Englewood<br />
Cliffs, 1976.<br />
T.R.Chandrupatla y A.D.Belegundu, Introduction to finite element in engineering, Prentice-<br />
Hall, Englewood Cliffs, 1991.<br />
El presente algoritmo no está condicionado a que la matriz sea <strong>de</strong>finida positiva y sólo requiere<br />
que sea no-singular. También pue<strong>de</strong> adaptarse a matrices no simétrica sin ninguna dificultad.<br />
136