BARRAS SOMETIDAS A FUERZA NORMAL Y FLEXION INCLINADAS

- pisos rígidos

- tirantes ( estructuras compuestas)

Ejemplo:
E=2000000 tonf/m^2

Previo al inicio de la resolución se determinan soluciones implícitas (fuerza y momento en el voladizo)

1- Identificamos las barras y definimos los grados de libertad

2.- Las matrices de rigidez de elementos inclinados sometidos a fuerzas axiales y flexión viene definida por
la siguiente expresión:

Con el apoyo del script en matlab se obtienen las matrices de rigidez de cada elemento,

Para barra 1

L= 5m A = 0.2800m2 I = 0.0114m4 c = 0.8000 s = 0.6000 sc = 0.4800

72470.27 52706.30 -3292.80 -72470.27 -52706.30 -3292.80

52706.30 41724.93 4390.40 -52706.30 -41724.93 4390.40
-3292.80 4390.40 18293.33 3292.80 -4390.40 9146.67
K1= -72470.27 -52706.30 3292.80 72470.27 52706.30 3292.80
-52706.30 -41724.93 -4390.40 52706.30 41724.93 -4390.40
-3292.80 4390.40 9146.67 3292.80 -4390.40 18293.33

Para barra 2

L = 3m A = 0.2800 m2 I = 0.0114 m4 c = 1 s = 0 sc = 0

186666.67 0.00 0.00 -186666.67 0.00 0.00

0.00 10162.96 15244.44 0.00 -10162.96 15244.44
0.00 15244.44 30488.89 0.00 -15244.44 15244.44
K2 = -186666.67 0.00 0.00 186666.67 0.00 0.00
0.00 -10162.96 -15244.44 0.00 10162.96 -15244.44
0.00 15244.44 15244.44 0.00 -15244.44 30488.89

Para barra 3

L= 3m A = 0.2000m2 I = 0.0042m4 c = 0 s = 1 sc = 0

3703.70 0.00 -5555.56 -3703.70 0.00 -5555.56

0.00 133333.33 0.00 0.00 -133333.33 0.00
-5555.56 0.00 11111.11 5555.56 0.00 5555.56
K3= -3703.70 0.00 5555.56 3703.70 0.00 5555.56
0.00 -133333.33 0.00 0.00 133333.33 0.00
-5555.56 0.00 5555.56 5555.56 0.00 11111.11
VECTORES DE MOMENTOS Y FUERZAS DE EMPOTRAMIENTO PERFECTO
%MOMENTOS DE EMPOTRAMIENTO BARRA 1

Me1=[-1.5 0 0.75 -1.5 0 -0.75]'

%MOMENTOS DE EMPOTRAMIENTO BARRA 2

Me2=[0 6 3 0 6 -3]'

%MOMENTOS DE EMPOTRAMIENTO BARRA 3

Me3=[0 0 0 0 0 0]'

VECTOR DE FUERZAS EXTERNAS

Fe =[0 0 0 0 0 0 0 -10 0 0 -14 -20]

3.- Ensamble de la matriz general de rigidez de la estructura y definiendo los vectores de fuerzas externas
así como vectores por fuerzas de empotramiento perfecto.

Los vectores de ordenamiento vienen definidos por la asignación de grados de libertad

Barra 1 V1=[1 2 3 7 8 9]
Barra 2 V2=[7 8 9 10 11 12]
Barra 3 V3=[4 5 6 10 11 12]

La matriz general de la estructura resultara:

72470.3 52706.3 -3292.8 0.0 0.0 0.0 -72470.3 -52706.3 -3292.8 0.0 0.0 0.0
52706.3 41724.9 4390.4 0.0 0.0 0.0 -52706.3 -41724.9 4390.4 0.0 0.0 0.0
-3292.8 4390.4 18293.3 0.0 0.0 0.0 3292.8 -4390.4 9146.7 0.0 0.0 0.0
0.0 0.0 0.0 3703.7 0.0 -5555.6 0.0 0.0 0.0 -3703.7 0.0 -5555.6
0.0 0.0 0.0 0.0 133333.3 0.0 0.0 0.0 0.0 0.0 -133333.3 0.0
0.0 0.0 0.0 -5555.6 0.0 11111.1 0.0 0.0 0.0 5555.6 0.0 5555.6
-72470.3 -52706.3 3292.8 0.0 0.0 0.0 259136.9 52706.3 3292.8 -186666.7 0.0 0.0
-52706.3 -41724.9 -4390.4 0.0 0.0 0.0 52706.3 51887.9 10854.0 0.0 -10163.0 15244.4
-3292.8 4390.4 9146.7 0.0 0.0 0.0 3292.8 10854.0 48782.2 0.0 -15244.4 15244.4
0.0 0.0 0.0 -3703.7 0.0 5555.6 -186666.7 0.0 0.0 190370.4 0.0 5555.6
0.0 0.0 0.0 0.0 -133333.3 0.0 0.0 -10163.0 -15244.4 0.0 143496.3 -15244.4
0.0 0.0 0.0 -5555.6 0.0 5555.6 0.0 15244.4 15244.4 5555.6 -15244.4 41600.0
4.- SOLUCION DEL SISTEMA ሼ‫ ܴ݁ܨ‬ሽ = ሾ‫ ݀݁ݎܭ‬ሿ ሼ‫ݑ‬ሽ + ሼ‫ܴ݁ܯ‬ሽ

Se identifican los grados de libertad restringidos ∆1 = ∆2 = ∆3 = ∆4 = ∆5 = ∆6 = 0

Todas las expresiones relacionadas con estos grados de libertad se eliminan de la matriz general

Kred=
259136.94 52706.30 3292.80 -186666.67 0.00 0.00
52706.30 51887.89 10854.04 0.00 -10162.96 15244.44
3292.80 10854.04 48782.22 0.00 -15244.44 15244.44
-186666.67 0.00 0.00 190370.37 0.00 5555.56
0.00 -10162.96 -15244.44 0.00 143496.30 -15244.44
0.00 15244.44 15244.44 5555.56 -15244.44 41600.00

De la misma manera los vectores de empotramiento perfecto y fuerzas externas restringidos serán:
MeR = [-1.5000 6.0000 2.2500 0 6.0000 -3.0000 ]
FeR = [ 0 -10 0 0 -14 -20 ]

Contando con todos los elementos se procede a solucionar el sistema,

ሼ‫ݑ‬ሽ = ሾ‫ ݀݁ݎܭ‬ሿିଵ ሾሼ‫ ܴ݁ܨ‬ሽ − ሼ‫ ܴ݁ܯ‬ሽሿ

Resultando los desplazamiento, giros o deformaciones

u7 = 0.00084342
u8= -0.00117530
u9 = 0.00015895
u10= 0.00083378
u11= -0.00023037
u12= -0.00023198

Reacciones ሼܴ݁ܽܿܿሽ = ሾ‫ ݃ܭ‬ሿ ሼ‫ݑ‬ሽ + ሼ‫ܴ݁ܯ‬ሽ

REACC =
-1.2007
5.2835
10.1411
-1.7993
30.7165
3.3433
-0.0000
-10.0000
-0.0000
-0.0000
-14.0000
-20.0000

Las reacciones estarán relacionadas con los valores de las primeras seis filas del vector resultado
FUERZAS INTERNAS

‫ܨ‬௜௡௧ = ሾ‫ܭ‬௜ ሿ ሼ‫ݑ‬௜ ሽ + ሼ‫݁ܯ‬௜ ሽ

f1 -1.20 f7 1.80 f4 -1.80

f2 5.28 f8 -4.72 f5 30.72
f3 10.14 f9 -10.09 f6 3.34
f7 -1.80 f10 -1.80 f10 1.80
f8 -5.28 f11 16.72 f11 -30.72
f9 10.09 f12 -22.05 f12 2.05
APLICANDO UN TIRANTE EN EL NUDO A

Matriz de rigidez del tirante adicional, que es la aplicable a elementos inclinados sometidos a compresión
axial.

K4=
527.47 -527.47 -527.47 527.47
-527.47 527.47 527.47 -527.47
-527.47 527.47 527.47 -527.47
527.47 -527.47 -527.47 527.47

Según los grados de libertad asignados el vector de ordenamiento para la barra 4 será:

v4= [7 8 13 14];

Con esta barra adicional la matriz de rigidez se agrandara en 2 filas con sus respectivas columnas

KG:
72470.3 52706.3 -3292.8 0.0 0.0 0.0 -72470.3 -52706.3 -3292.8 0.0 0.0 0.0 0.0 0.0
52706.3 41724.9 4390.4 0.0 0.0 0.0 -52706.3 -41724.9 4390.4 0.0 0.0 0.0 0.0 0.0
-3292.8 4390.4 18293.3 0.0 0.0 0.0 3292.8 -4390.4 9146.7 0.0 0.0 0.0 0.0 0.0
0.0 0.0 0.0 3703.7 0.0 -5555.6 0.0 0.0 0.0 -3703.7 0.0 -5555.6 0.0 0.0
0.0 0.0 0.0 0.0 133333.3 0.0 0.0 0.0 0.0 0.0 -133333.3 0.0 0.0 0.0
0.0 0.0 0.0 -5555.6 0.0 11111.1 0.0 0.0 0.0 5555.6 0.0 5555.6 0.0 0.0
-72470.3 -52706.3 3292.8 0.0 0.0 0.0 259664.4 52178.8 3292.8 -186666.7 0.0 0.0 -527.5 527.5
-52706.3 -41724.9 -4390.4 0.0 0.0 0.0 52178.8 52415.4 10854.0 0.0 -10163.0 15244.4 527.5 -527.5
-3292.8 4390.4 9146.7 0.0 0.0 0.0 3292.8 10854.0 48782.2 0.0 -15244.4 15244.4 0.0 0.0
0.0 0.0 0.0 -3703.7 0.0 5555.6 -186666.7 0.0 0.0 190370.4 0.0 5555.6 0.0 0.0
0.0 0.0 0.0 0.0 -133333.3 0.0 0.0 -10163.0 -15244.4 0.0 143496.3 -15244.4 0.0 0.0
0.0 0.0 0.0 -5555.6 0.0 5555.6 0.0 15244.4 15244.4 5555.6 -15244.4 41600.0 0.0 0.0
0.0 0.0 0.0 0.0 0.0 0.0 -527.5 527.5 0.0 0.0 0.0 0.0 527.5 -527.5
0.0 0.0 0.0 0.0 0.0 0.0 527.5 -527.5 0.0 0.0 0.0 0.0 -527.5 527.5
Para solucionar el sistema nuevo, los vectores de fuerzas externas y de momentos de empotramiento serán
los mismos, con el complemento de los grados de libertad adicionales.

MOMENTOS DE EMPOTRAMIENTO PERFECTO

Met = [-1.5 0 0.75 0 0 0 -1.5 6.0 2.25 0 6.0 -3.0 0 0]

VECTOR DE FUERZAS EXTERNAS

Fe2 = [0 0 0 0 0 0 0 -10 0 0 -14 -20 0 0]

A las restricciones se adicionan los del apoyo del tirante resultando

∆1 = ∆2 = ∆3 = ∆4 = ∆5 = ∆6 = ∆13 = ∆14

VECTORES DE RESTRICCIONES

VR = [1 2 3 4 5 6 13 14 ]

MATRIZ GENERAL REDUCIDA

259664.40 52178.84 3292.80 -186666.67 0.00 0.00

52178.84 52415.36 10854.04 0.00 -10162.96 15244.44
3292.80 10854.04 48782.22 0.00 -15244.44 15244.44
-186666.67 0.00 0.00 190370.37 0.00 5555.56
0.00 -10162.96 -15244.44 0.00 143496.30 -15244.44
0.00 15244.44 15244.44 5555.56 -15244.44 41600.00

Vectores de Momentos de empotramiento y fuerzas externas

MetR = [-1.5 6.0 2.25 0 6.0 -3.0]’

Fe2R = [0 -10 0 0 -14 -20]’

Resultados de desplazamientos

ut =
0.0007
-0.0010
0.0001
0.0007
-0.0002
-0.0003
Para concluir con este punto se comparan los desplazamientos horizontales de la parte superior del pórtico.

Sin tirante Con tirante

7.00 0.000843 0.000740
8.00 -0.001175 -0.001040
9.00 0.000159 0.000147
10.00 0.000834 0.000733
11.00 -0.000230 -0.000225
12.00 -0.000232 -0.000262

Efectuándose una pequeña reducción en el desplazamiento de la viga del pórtico así como en los giros.

NUDOS RIGIDOS

Para considerar esta alternativa de análisis, se asume que los desplazamientos horizontales en un mismo
nivel de piso son iguales, logrando así eliminar una incógnita y por lo tanto simplificando el sistema de
ecuaciones a resolver.

En el actual ejemplo ∆7 = ∆10

En la matriz reducida esto puede asumirse, complementando las filas y columnas de los desplazamientos
igualados

En la matriz reducida original:

Kred=
259136.94 52706.30 3292.80 -186666.67 0.00 0.00
52706.30 51887.89 10854.04 0.00 -10162.96 15244.44
3292.80 10854.04 48782.22 0.00 -15244.44 15244.44
-186666.67 0.00 0.00 190370.37 0.00 5555.56
0.00 -10162.96 -15244.44 0.00 143496.30 -15244.44
0.00 15244.44 15244.44 5555.56 -15244.44 41600.00

Paso 1: Inicialmente se sumaran las filas correspondientes a los grados de libertad 7 y 10

Kred= 72470.27 52706.30 3292.80 3703.70 0.00 5555.56

52706.30 51887.89 10854.04 0.00 -10162.96 15244.44
3292.80 10854.04 48782.22 0.00 -15244.44 15244.44
0.00 -10162.96 -15244.44 0.00 143496.30 -15244.44
0.00 15244.44 15244.44 5555.56 -15244.44 41600.00

Paso 2: Se suman las columnas

76173.98 52706.30 3292.80 0.00 5555.56
52706.30 51887.89 10854.04 -10162.96 15244.44
KND= 3292.80 10854.04 48782.22 -15244.44 15244.44
0.00 -10162.96 -15244.44 143496.30 -15244.44
5555.56 15244.44 15244.44 -15244.44 41600.00
Los vectores de fuerza externas y empotramientos sufrirán modificaciones similares en sus términos

MeR =
-1.5000
6.0000
2.2500
6.0000
-3.0000

FeR =
0
-10
0
-14
-20

El vector de desplazamientos respuesta bajo esta técnic que se obtiene se compara con los del resultado
original

Considerando Sin considerar

pisos rígidos Pisos rígidos
0.000841 0.000843
-0.001172 -0.001175
0.000159 0.000159
0.000834
-0.000230 -0.000230
-0.000234 -0.000232

Por lo que el hecho de considerar esta hipótesis nos lleva a obtener valores con una precisión aceptable
respecto del desplazamiento real.
 MATLAB
RUTINA EN MATHCAD- RESOLUCION DE PORTICO

% problema numero 3 ejemplo barras inclinadas sometidas a flexion

clc %borrar lo que esta en pantalla
clear all %borrar variables

%DATOS DE ENTRADA
E=2000000%TONF/m^2

%MATRIZ DE RIGIDEZ DE BARRA 1

%DATOS DE LA SECCION
b=0.4;% m
h=0.7;%m
L=5%m
X=4;
Y=3;
A=b*h
I=b*h^3/12
c=X/L
s=Y/L
%s=0.59;
%c=0.8;
sc=s*c
disp('MATRICES DE RIGIDEZ DE BARRAS')

k1=[c^2*A*E/L+s^2*12*E*I/L^3 sc*((E*A/L)-(12*E*I/L^3)) -s*6*E*I/L^2 (-c^2*A*E/L)-

(s^2*12*E*I/L^3) -sc*(E*A/L-12*E*I/L^3) -s*6*E*I/L^2
sc*((E*A/L)-(12*E*I/L^3)) s^2*A*E/L+c^2*12*E*I/L^3 c*6*E*I/L^2 -sc*(E*A/L-12*E*I/L^3 )
-s^2*A*E/L-c^2*12*E*I/L^3 c*6*E*I/L^2
-s*6*E*I/L^2 c*6*E*I/L^2 4*E*I/L s*6*E*I/L^2 -c*6*E*I/L^2 2*E*I/L
(-c^2*E*A/L)-(s^2*12*E*I/L^3) -sc*((E*A/L)-(12*E*I/L^3)) s*6*E*I/L^2
(c^2*E*A/L)+(s^2*12*E*I/L^3) sc*(E*A/L-12*E*I/L^3) s*6*E*I/L^2
-sc*(E*A/L-12*E*I/L^3) -s^2*A*E/L-c^2*12*E*I/L^3 -c*6*E*I/L^2 sc*(E*A/L-12*E*I/L^3)
s^2*A*E/L+c^2*12*E*I/L^3 -c*6*E*I/L^2
-s*6*E*I/L^2 c*6*E*I/L^2 2*E*I/L s*6*E*I/L^2 -c*6*E*I/L^2 4*E*I/L]

%MATRIZ DE RIGIDEZ DE BARRA 2

%DATOS DE LA SECCION
b=0.4;% m
h=0.7;%m
L=3%m
X=3;
Y=0;
A=b*h
I=b*h^3/12
c=X/L
s=Y/L
%s=0.59;
%c=0.8;
sc=s*c;

k2=[c^2*A*E/L+s^2*12*E*I/L^3 sc*((E*A/L)-(12*E*I/L^3)) -s*6*E*I/L^2 (-c^2*A*E/L)-

(s^2*12*E*I/L^3) -sc*(E*A/L-12*E*I/L^3) -s*6*E*I/L^2
sc*((E*A/L)-(12*E*I/L^3)) s^2*A*E/L+c^2*12*E*I/L^3 c*6*E*I/L^2 -sc*(E*A/L-12*E*I/L^3 )
-s^2*A*E/L-c^2*12*E*I/L^3 c*6*E*I/L^2
-s*6*E*I/L^2 c*6*E*I/L^2 4*E*I/L s*6*E*I/L^2 -c*6*E*I/L^2 2*E*I/L
(-c^2*E*A/L)-(s^2*12*E*I/L^3) -sc*((E*A/L)-(12*E*I/L^3)) s*6*E*I/L^2
(c^2*E*A/L)+(s^2*12*E*I/L^3) sc*(E*A/L-12*E*I/L^3) s*6*E*I/L^2
-sc*(E*A/L-12*E*I/L^3) -s^2*A*E/L-c^2*12*E*I/L^3 -c*6*E*I/L^2 sc*(E*A/L-12*E*I/L^3)
s^2*A*E/L+c^2*12*E*I/L^3 -c*6*E*I/L^2
-s*6*E*I/L^2 c*6*E*I/L^2 2*E*I/L s*6*E*I/L^2 -c*6*E*I/L^2 4*E*I/L]
%MATRIZ DE RIGIDEZ DE BARRA 3
%DATOS DE LA SECCION
b=0.4;% m
h=0.5;%m
L=3%m
X=0;
Y=3;
A=b*h
I=b*h^3/12
c=X/L
s=Y/L
sc=s*c

k3=[c^2*A*E/L+s^2*12*E*I/L^3 sc*((E*A/L)-(12*E*I/L^3)) -s*6*E*I/L^2 (-c^2*A*E/L)-

(s^2*12*E*I/L^3) -sc*(E*A/L-12*E*I/L^3) -s*6*E*I/L^2
sc*((E*A/L)-(12*E*I/L^3)) s^2*A*E/L+c^2*12*E*I/L^3 c*6*E*I/L^2 -sc*(E*A/L-12*E*I/L^3 )
-s^2*A*E/L-c^2*12*E*I/L^3 c*6*E*I/L^2
-s*6*E*I/L^2 c*6*E*I/L^2 4*E*I/L s*6*E*I/L^2 -c*6*E*I/L^2 2*E*I/L
(-c^2*E*A/L)-(s^2*12*E*I/L^3) -sc*((E*A/L)-(12*E*I/L^3)) s*6*E*I/L^2
(c^2*E*A/L)+(s^2*12*E*I/L^3) sc*(E*A/L-12*E*I/L^3) s*6*E*I/L^2
-sc*(E*A/L-12*E*I/L^3) -s^2*A*E/L-c^2*12*E*I/L^3 -c*6*E*I/L^2 sc*(E*A/L-12*E*I/L^3)
s^2*A*E/L+c^2*12*E*I/L^3 -c*6*E*I/L^2
-s*6*E*I/L^2 c*6*E*I/L^2 2*E*I/L s*6*E*I/L^2 -c*6*E*I/L^2 4*E*I/L]

%CREACION DE LA MATRIZ GENERAL

KG=zeros(12);
%DEFINIMOS VECTORES DE ORDENAMIENTO
disp('VECTORES DE ORDENAMIENTO')
V1=[1 2 3 7 8 9]
V2=[7 8 9 10 11 12]
V3=[4 5 6 10 11 12]

%INGRESAR LA MATRIZ K1 A LA MATRIZ KG

%KG(a1,a1)=KG(a1,a1)+k1

KG(V1,V1)=KG(V1,V1)+k1;%para barra 1
KG(V2,V2)=KG(V2,V2)+k2;%para barra 2
KG(V3,V3)=KG(V3,V3)+k3;%para barra 3

%restricciones
KGR=KG;% MANTENGO LA VARIABLE KG Y CREO KGR PARA OPERAR
disp('VECTORES DE RESTRICCIONES')
VR=[1 2 3 4 5 6]%RESTRICCIONES POR APOYO
KGR(VR,:)=[];%ELIMINO FILAS
KGR(:,VR)=[];%ELIMINO COLUMNAS

disp('MATRIZ REDUCIDA')
KGR

%%
%VECTORES DE EMPOTRAMIENTO PERFECTO

Me1=[-1.5 0 0.75 -1.5 0 -0.75]';%MOMENTOS DE EMPOTRAMIENTO BARRA 1

Me2=[0 6 3 0 6 -3]';%MOMENTOS DE EMPOTRAMIENTO BARRA 2
Me3=[0 0 0 0 0 0]';%MOMENTOS DE EMPOTRAMIENTO BARRA 3
Me=zeros(12,1);
Me(V1)=Me(V1)+Me1;
Me(V2)=Me(V2)+Me2;
Me(V3)=Me(V3)+Me3;
%disp('VECTOR DE MOMENTOS DE EMPOTRAMIENTO PERFECTO')
Me;
%%
%vector de fuerzas externas nodales
Vfe=[8 11 12];
FE=[-10 -14 -20];
Fe=zeros(1,12);
Fe(Vfe)=Fe(Vfe)+FE;
%disp('VECTOR DE FUERZAS EXTERNAS')
Fe;
%Fei=[0 0 0 0 0 0 0 -10 0 0 -14 -20]
%Restringiendo los vectores de fuerzas
disp('vector de momentos y fuerzas externas restringido')
MeR=Me;
MeR(VR)=[]

FeR=Fe';
FeR(VR)=[]

%%
%resolviendo el sistema
u=inv(KGR)*(FeR-MeR)
disp('RESULTADOS DE DESPLAZAMIENTOS')
u1=KGR\(FeR-MeR)
ust=u1;
%%
%REACCIONES
disp('REACCIONES')

vu=zeros(12,1);
vud=[7 8 9 10 11 12];
vu(vud)=vu(vud)+u1
KG;
REACC=KG*vu+Me
%%
%diagramas de cuerpo libre
disp('FUERZAS INTERNAS')

u12=vu(V1);
F1=k1*u12+Me1

u2=vu(V2);
F2=k2*u2+Me2

u3=vu(V3);
F3=k3*u3+Me3
%%
%APLICANDO EL TIRANTE AL NUDO IZQUIERDO SUPERIOR

%MATRIZ DE RIGIDEZ DEL TIRANTE

%Barra 1
disp('TIRANTE')
E=2.1*10000000;%t/m2

A=0.000201%m2...area
L=4%cm... longitud
S2=0.707*0.707;C2=S2;SC=-S2;%.... factores de rotacion

k4=[E*A/L*C2 E*A/L*SC E*A/L*-C2 E*A/L*-SC

E*A/L*SC E*A/L*S2 E*A/L*-SC E*A/L*-S2
E*A/L*-C2 E*A/L*-SC E*A/L*C2 E*A/L*SC
E*A/L*-SC E*A/L*-S2 E*A/L*SC E*A/L*S2]

% ES NECESARIO ENSAMBLAR OTRA MATRIZ GENERAL

KG2=zeros(14)
v4=[7 8 13 14];% ordenamiento para cable externo

KG2(V1,V1)=KG2(V1,V1)+k1;%para barra 1
KG2(V2,V2)=KG2(V2,V2)+k2;%para barra 2
KG2(V3,V3)=KG2(V3,V3)+k3;%para barra 3
KG2(v4,v4)=KG2(v4,v4)+k4;%para barra 3

% momentos de empotramiento
%VECTORES DE EMPOTRAMIENTO PERFECTO
disp('Momentos de empotramiento Perfecto')
Met=zeros(14,1);
Met(V1)=Met(V1)+Me1;
Met(V2)=Met(V2)+Me2

%vector de fuerzas externas nodales

disp('VECTOR DE FUERZAS EXTERNAS')
Vfe=[8 11 12];
FE=[-10 -14 -20];
Fe2=zeros(1,14);
Fe2(Vfe)=Fe2(Vfe)+FE
Fe2=Fe2'

disp('VECTOR DE RESTRICCIONES ')

KGR2=KG2;
VR=[1 2 3 4 5 6 13 14]%RESTRICCIONES POR APOYO
KGR2(VR,:)=[];%ELIMINO FILAS
KGR2(:,VR)=[];%ELIMINO COLUMNAS
disp('MATRIZ REDUCIDA')
KGR2

MetR=Met; MetR(VR)=[]

Fe2R=Fe2; Fe2R(VR)=[]

%resolviendo el sistema
disp('RESULTADOS DE DESPLAZAMIENTOS')

ut=inv(KGR2)*(Fe2R-MetR)

%%
%Nudos rigidoz
%u7=u10 , dezplasamientos igualados

KND=KGR; % KND, matriz modificada para emplear la metodologia de nudos rigidos

% recordando que en la matriz reducida la primera fila y columna se
% relacionan con el grado de libertad 7 y la cuarta fila y columna con el
% grado de libertad 10

% primer paso se suman las filas de los grados de libertad implicados

KND(1,:)= KND(1,:)+KND(4,:);
KND(4,:)=[]
%segundo paso se suman las columnas de los grados de libertad implicados
KND(:,1)= KND(:,1)+KND(:,4);
KND(:,4)=[]
% de igual manera los momentos de empotramiento y vectores de fuerza deben
% reacomodarse para que el sistema sea consistente
%vector de fuerzas de empotramiento perfecto
MeR(1)=MeR(1)+MeR(4);
MeR(4)=[]
%vector de fuerzas externas nodales
FeR(1)=FeR(1)+FeR(4);
FeR(4)=[]

%Resolviendo el sistema con estas modificaciones

uPR=KND\(FeR-MeR)p