DATOS - EJERCICIO NUMEROS ENTEROS

00a - GRADOS DE LIBERTAD

n=20;

00b - DATOS: MATRIZ DE RIGIDEZ - REPRESENTADA MEDIANTE VECTOR S

S={{0,1,2,5,6,7,12,18,20,22,30,39,50,51,64,74,83,87,99,104,114},{64,87,30,35,91,70,98,42,79,91,23,6,74,58,94,43,34,84,94,8,15,56,84,63,66,15,16,50,51,17,38,18,98,10,23,50,36,42,32,88,14,46,8,51,91,2,50,71,88,68,41,99,36,90,88,47,57,17,5,15,36,36,90,78,67,58,83,36,3,92,57,92,23,62,34,42,88,39,34,83,66,73,40,70,51,43,56,83,98,66,27,73,87,26,78,30,53,13,76,44,94,67,81,11,93,52,96,98,4,9,33,33,27,72}};

00c- DATOS: VECTOR DE FUERZAS NODALES CONSISTENTES

f={33,28,30,64,3,62,57,42,99,43,14,15,42,91,61,51,30,91,70,13};

01a- RECONSTRUCCION DE LA SKYMATRIZ ORIGINAL A PARTIR DEL VECTOR S

Ko=X[S];

Ko//MatrixForm

01b - METODO DIRECTO DE OBTENCION DE  DESPLAZAMIENTOS SOLUCION

1 - METODO DIRECTO

- A COMPROBAR MEDIANTE EL PROCEDIMIENTO DE FACTORIZACION (1a FASE) Y SOLUCION (2a FASE)

usmd=X.f;

%//MatrixForm

02a- 1a FASE - FACTORIZACION DE LA SKYMATRIZ ORIGINAL A PARTIR DEL VECTOR S

2 - FACTORIZACION AUTOMATICA

COMANDO A UTILIZAR

Sf=X[S,10.^(-12)][[1]]

RECONSTRUCCION DE LA SKYMATRIZ FACTORIZADA A PARTIR DEL VECTOR Sf

X[Sf];

%//MatrixForm

DESPLAZAMIENTOS SOLUCION - OBTENIDOS DIRECTAMENTE A PARTIR DE LA FACTORIZACION DE LA SKYMATRIZ ORIGINAL

usfs=X[Sf,f];

%//MatrixForm

02b- MATRICES QUE CONSTITUYEN LA FACTORIZACION DE LA SKYMATRIZ ORIGINAL

3 - FACTORIZACION PASO A PASO

MATRIZ D

MATRIZ D INVERSA - D1

D1=X[Sf];

(*%//MatrixForm*)

MATRIZ Df

Df=Inverse[D1];

%//MatrixForm

MATRIZ U

U=X[Sf];

%//MatrixForm

MATRIZ L

L=X[U];

%//MatrixForm

MATRIZ L INVERSA

Inverse[L] //MatrixForm

COMPROBACION DE LA FACTORIZACION DE LA SKYMATRIZ ORIGINAL

Kc=X;

%//MatrixForm

COMANDOS DIRECTOS DE OBTENCION DE LA SKYMATRIZ ORIGINAL A PARTIR DE LA SKYMATRIZ FACTORIZADA

So=X[Sf];

X[So];

%//MatrixForm

02c - OBTENCION DE  DESPLAZAMIENTOS SOLUCION

02c- 2a FASE - 1a PARTE - SOLUCION INTERMEDIA

SolucionIntermedia=X.f;

Xi=SolucionIntermedia;

02c- 2a FASE - 2a PARTE - SOLUCION POSTERIOR

SolucionPosterior=X.Xi;

usrf=SolucionPosterior;

%//MatrixForm

02c- COMPROBACION - SOLUCIONES

X[S,usrf]

f

Ok