Solution of Equilibrium Equations in Static Analysis: LDLT Solution ...
Solution of Equilibrium Equations in Static Analysis: LDLT Solution ...
Solution of Equilibrium Equations in Static Analysis: LDLT Solution ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Physical Interpretation<br />
A physical <strong>in</strong>terpretation <strong>of</strong> the operations performed <strong>in</strong> a Gauss elim<strong>in</strong>ation:<br />
5 -4 1 0<br />
6 -4 1<br />
symmetric<br />
6 -4<br />
U 1<br />
U 2<br />
U 3<br />
U 4<br />
0<br />
0<br />
0<br />
5 0<br />
First equation: 5 U 1 –4 U 2 + U 3 = 0 ⇔ U 1 = 4/5 U 2 –1/5 U 3<br />
Elim<strong>in</strong>ation i <strong>of</strong> U 1 from equations 2, 3 and 4 yields ild the lower right 3 x 3 submatrix<br />
which we get after the first step <strong>of</strong> the Gauss elim<strong>in</strong>ation <strong>of</strong> the orig<strong>in</strong>al matrix:<br />
14/5 -16/5 1<br />
-16/5 29/5 -4<br />
1 -4 5<br />
U 2<br />
0<br />
Stiffness matrix correspond<strong>in</strong>g to<br />
U 3<br />
0 beam after release <strong>of</strong> d<strong>of</strong> 1.<br />
U 4 0 ( d<strong>of</strong> 1“statically condensed out”)<br />
….. 5/6 is stiffness matrix <strong>of</strong> beam after release <strong>of</strong> d<strong>of</strong>s 1, 2 and 3 (cf.<br />
Gauss elim<strong>in</strong>ation: f<strong>in</strong>al upper triangular matrix).<br />
15-Jun-07<br />
Method <strong>of</strong> F<strong>in</strong>ite Elements I<br />
7