18.01.2014 Views

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 ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!