EJERCICIOS PROPUESTOS

Dados los sistemas Q – q y Q* - q* en los ejercicios 1, 2 y 3. Se pide calcular:

 La matriz T, tal que: q=Tq*

 La matriz T1, tal que: Q=T1Q*
 Comprobar que :
t
T =T −1
1

EJERCICIO Nº 01

Solución:
Calculamos la matriz T:
¿ ¿
q 1 =1 y q i=0 cuando i≠1

T 11=0 T 21=0 T 31=0 T 41=1

 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
q 2 =1 y q i=0 cuando i≠2

T 12=0 T 22=0 T 32=0 T 42=0

¿ ¿
q 3 =1 y q i=0 para i≠3

T 13 =0 T 23=1/ L T 33 =0 T 43=0
¿ ¿
q 4 =1 y q i=0 cuando i≠4

T 14 =1 T 24=0 T 34 =0 T 44 =0

2
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

[ ]
Por lo tanto, la matriz de transformació n es:
0 0 0 1
0 0 1/ L 0
T=
0 1 0 0
1 0 0 0
Calculamos la matriz
T1 :
¿ ¿
Q 1 =1 y Q i =0 para i≠1
Problema Primario:

Equilibrio de miembros

Equilibrio de nudos
Nudo B:

3
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

Nudo D:

T 111 =0 T 121 =0 T 131 =0 T 141 =1

¿ ¿
Q 2 =1 y Q i =0 cuando i≠2
Problema Primario

Equilibrio de nudos
Nudo B

Nudo D

T 112 =0 T 122 =0 T 132 =1 T 142 =0

4
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
Q 3 =1 y Q i =0 para i≠3
Problema Primario

Equilibrio de nudos
Nudo B

Nudo D

T 113 =0 T 123 =L T 133 =0 T 143 =0

¿ ¿
Q 4 =1 y Q i =0 para i≠4
Problema Primario

5
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

Equilibrio de nudos
Nudo B

Nudo D

T 114 =1 T 124 =0 T 134 =0 T 144 =0

T1

[ ]
Por lo tanto, la matriz de transformació n es:
0 0 0 1
0 0 L 0
T 1=
0 1 0 0
1 0 0 0
T t =T −1

[ ]
Se puede notar que: 1

0 0 0 1
0 0 1 0
T t =T −1=
1 0 1/ L 0 0
1 0 0 0

6
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

EJERCICIO Nº 02

Solución
Calculamos la matriz T:
¿ ¿
q 1 =1 y q i=0 para i≠1

tg53 º=x=4/3 T 11=0 T 21=4 /3 T 31=0 T 41=1

¿ ¿
q 2 =1 y q i=0 para i≠2

T 12=0 T 22=0 T 32=0 T 42=1

7
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
q 3 =1 y q i=0 para i≠3

T 13 =1 T 23=0 T 33 =0 T 43=0

¿ ¿
q 4 =1 y q i=0 para i≠4

T 14 =1 T 24=0 T 34 =1 T 44 =0

[ ]
Por lo tanto, la matriz de transformació n es:
0 0 1 0
4/3 0 0 0
T=
0 0 0 1
1 1 0 0

Calculamos la matriz
T1 :

¿ ¿
Q 1 =1 y Q i =0 para i≠1

8
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

Problema Primario

Equilibrio de miembros

Equilibri
o de nudos
Nudo A

Nudo B

Nudo C

T 111 =0 T 121 =3/4 T 131 =0 T 141 =0

9
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
Q 1 =2 y Q i =0 para i≠2
Problema Primario

Equilibrio de nudos
Nudo A

Nudo B

Nudo C

T 112 =0 T 122 =0 T 132 =0 T 142 =1

10
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
Q 3 =1 i
y Q =0 para i≠3
Problema Primario

Equilibrio de nudos
Nudo A

Nudo B

Nudo C

T 113 =1 T 123 =0 T 133 =0 T 143 =0

11
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
Q 4 =1 i
y Q =0 para i≠4
Problema Primario

Equilibrio de nudos
Nudo A

Nudo B

Nudo C

T 114 =0 T 124 =0 T 134 =1 T 144 =0

Por lo tanto, la matriz de transformació n
T1 es:

[ ]
0 0 1 0
3 /4 0 0 0
T 1=
0 0 0 1
0 1 0 0
T t =T −1
Se puede notar que: 1

[ ]
0 4/3 0 0
0 0 0 1
T t =T −1=
1 1 0 0 0
0 0 1 0

12
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

EJERCICIO Nº 03

Solución
Calculamos la matriz T:
¿ ¿
q 1 =1 y q i=0 para i≠1

T 11=0 T 21=0 T 31=0 T 41=1 T 51=0 T 61=0 T 71=0

¿ ¿
q 2 =1 y q i=0 para i≠2

T 12 =0 T 22=0 T 32=0 T 42=0 T 52=0 T 62=0 T 72=1

13
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
q 3 =1 y q i=0 para i≠3

T 13 =0 T 23=0 T 33 =0 T 43=0. 4 T 53 =1 T 63 =0 T 73 =0
¿ ¿
q 4 =1 y q i=0 para i≠4

T 14 =0 T 24=0 T 34 =0 T 44=0 T 54 =0 T 64=1 T 74 =0

¿ ¿
q 5 =1 y q i=0 para i≠5

T 15 =1 T 25=0 T 35 =0 T 45=0 T 55 =0 T 65=0 T 75 =0

14
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
q 6 =1 y q i=0 para i≠6

T 16 =0 T 26 =1 T 36 =0 T46  0 T 56 =0 T 66=0 T 76 =0
¿ ¿
q 7 =1 y q i=0 para i≠7

T 17 =0 T 27=0 T 37 =1 T 47=0 T 57 =0 T 67=0 T 77 =0

[ ]
Por lo tanto, la matriz de transformació n T es:

0 0 0 0 1 0 0
0 0 0 0 0 1 0
0 0 0 0 0 0 1
T= 1 0 0.4 0 0 0 0
0 0 1 0 0 0 0
0 0 0 1 0 0 0
0 1 0 0 0 0 0

Calculamos la matriz
T1 :
¿ ¿
Q 1 =1 y Q i =0 para i≠1

15
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

Problema Primario

Equilibrio de miembros

Equilibrio de nudos
Nudo B

Nudo C

Nudo D

T 111 =0 T 121 =0 T 131 =0 T 141 =1 T 151 =−0.40 T 161 =0 T 171 =0

16
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
Q 2 =1 i
y Q =0 para i≠2
Problema Primario

Equilibrio de nudos
Nudo B

Nudo C

Nudo D

T 112 =0 T 122 =0 T 132 =0 T 142 =0 T 152 =0 T 162 =0 T 172 =1

¿ ¿
Q 3 =1 i
y Q =0 para i≠3
Problema Primario

17
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

Equilibrio de nudos
Nudo B

Nudo C

Nudo D

T 113 =0 T 123 =0 T 133 =0 T 143 =0 T 153 =1 T 163 =0 T 173 =0

¿ ¿
Q 3 =1 y Q i =0 para i≠3
Problema Primario

Equilibrio de nudos
Nudo B

Nudo C

18
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

Nudo D

T 114 =0 T 124 =0 T 134 =0 T 144 =0 T 154 =0 T 164 =1 T 174 =0

¿ ¿
Q 5 =1 i
y Q =0 para i≠5
Problema Primario

Equilibrio de nudos
Nudo B

Nudo C

Nudo D

T 115 =1 T 125 =0 T 135 =0 T 145 =0 T 155 =0 T 165 =0 T 175 =0

19
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
Q 6 =1 i
y Q =0 para i≠6
Problema Primario

Equilibrio de nudos
Nudo B

Nudo C

Nudo D

T 116 =0 T 126 =1 T 136 =0 T 146 =0 T 156 =0 T 166 =0 T 176 =0

20
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
Q 7 =1 i
y Q =0 para i≠7
Problema Primario

Equilibrio de nudos
Nudo B

Nudo C

Nudo D

T 117 =0 T 127 =0 T 137 =1 T 147 =0 T 157 =0 T 167 =0 T 177 =0

Por lo tanto, la matriz de transformació n

T1 es:

21
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

[ ]
0 0 0 0 1 0 0
0 0 0 0 0 1 0
0 0 0 0 0 0 1
T 1= 1 0 0 0 0 0 0
−0. 40 0 1 0 0 0 0
0 0 0 1 0 0 0
0 1 0 0 0 0 0

t
T =T −1
Se puede notar que: 1

[ ]
0 0 0 1 0 0 0
0 0 0 0 0 0 1
0 0 0 0 . 40 1 0 0
T t =T −1= 0 0 0 0 0 1 0
1
1 0 0 0 0 0 0
0 1 0 0 0 0 0
0 0 1 0 0 0 0

22
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

EJERCICIO Nº 04
Si sobre la estructura de la figura 5.8.1 actú a el sistema de cargas presentado a continuació n. Se pide calcular
los vectores de cargas generalizadas Q y Q* para los sistemas de coordenadas presentados en las figuras 5.8.2 y
t
5.8.3. Comprobar el resultado obtenido utilizando la matriz T del ejemplo 5, si se sabe que Q∗¿T Q .

Solución
i) Equilibrio de elementos: Sistema Q

Problema Primario

23
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

Equilibrio de nudos
Nudo A

Nudo B

Nudo C

3 W 1 −W 4 3 W 5−W 6
R1 = R2 =−W 2 R3 =0 R4 =−W 3 R5 = R7=W 7
3 3

[ ]
Por lo tanto, la matriz de vectores Q es:
3 W 1 −W 4
3
−W 2
0
Q=
−W 3
3 W 5 −W 6
3
W7
Problema Primario: Sistema Q*

24
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

Equilibrio de nudos
Nudo A

Nudo B

Nudo C

R1 =−W 2 R2 =0 R3 =−W 3 R4 =W 4 −3 W 1 R5 =3 W 5 −W 6 R7 =W 7
Por lo tanto, la matriz de vectores Q* es:

[ ]
−W 2
0
−W 3
Q∗¿
W 4 −3 W 1
3W 5 −W 6
W7

][ ] [ ]
iii) Comprobamos que Q*=TtQ
3 W 1 −W 4

[
0 1 0 0 0 0 3 −W 2
0 0 1 0 0 0 −W 2 0
0 0 0 1 0 0 0 −W 3
Q∗¿ =
−3 0 0 0 0 0 −W 3 W 4 −3W 1
0 0 0 0 3 0 3 W 5 −W 6 3 W 5 −W 6
0 0 0 0 0 1 W7
3
W7

25
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

EJERCICIO Nº 05
Para el pó rtico plano de la figura 5.2 se desea calcular la relació n inversa que existe entre las coordenadas
generalizadas. Por consiguiente ahora se tiene:

Calcular la matriz T tal que q=Tq∗.

Solución
i) Calculamos la matriz T:
¿ ¿
1
q =1 y
i
q =0 para i≠1

tg α=x=3/2 T 11=0 T 21=3/2 T 31=0

¿ ¿
q 2 =1 y q i=0 para i≠2

T 12=0 T 22=0 T 32=1

26
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
q 3 =1 y q i=0 para i≠3

ctg α=x=3/2 T 13 =0 T 23=0 T 33 =0

Por lo tanto, la matriz de transformació n T es:

[ ]
0 0 0
T = 3/2 0 1
0 1 0

27
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

EJERCICIO Nº 06
Para la siguiente estructura compuesta por elementos totalmente flexibles se dan dos sistemas de coordenadas
generalizados, el primero difiere del segundo en la coordenada uno, está inclinado un á ngulo θ con respecto a
la horizontal. Se pide calcular la matriz T tal que q=Tq∗.

Solución:
i) Calculamos la matriz T:
¿ ¿
1
q =1 y
i
q =0 para i≠1

T 11=cosθ T 21=0 T 31=0 T 41=0 T 51=0 T 61=0

¿ ¿
q 2 =1 y q i=0 para i≠2

T 12=sen θ T 22=1 T 32=0 T 42=0 T 52=0 T 62=0

28
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

¿ ¿
q 3 =1 y q i=0 para i≠3

T 13 =0 T 23=0 T 33 =1 T 43=0 T 53 =0 T 63 =0
¿ ¿
q 4 =1 y q i=0 para i≠4

T 14 =0 T 24=0 T 34 =0 T 44=1 T 54 =0
T 64=0
¿ ¿
q 5 =1 y q i=0 para i≠5

29
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

T 15 =0 T 25=0 T 35 =0 T 45=0 T 55 =1 T 65=0

¿ ¿
q 6 =1 y q i=0 para i≠6

T 16 =0 T 26=0 T 36 =0 T 46=0 T 56 =0 T 66=1

Por lo tanto, la matriz de transformació n T es:

[ ]
cosθ senθ 0 0 0 0
0 1 0 0 0 0
T= 0 0 1 0 0 0
0 0 0 1 0 0
0 0 0 0 1 0
0 0 0 0 0 1

30
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

EJERCICIO Nº 06
Las siguientes figuras muestran una estructura y las fuerzas externas que actú an, así como los diagramas de
deformació n elementales, para un determinado sistema de coordenadas. Calcular el vector de cargas
generalizadas Q.

Solución
i) Equilibrio de elementos

Problema Primario

31
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

Equilibrio de nudos
Nudo A

Nudo B

R1 =W 1 R2 =W 3 +W 6 −W 2 L R3 =W 4
Por lo tanto, la matriz de vectores Q es:

[ ]
W1
Q= W 3 +W 6 −W 2 L
W4

32
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

EJERCICIO Nº 08
Elaborar un programa de computació n para obtener el vector de cargas generalizadas Q:

Solución
a) Codificación del Programa
El programa está codificado en el lenguaje de programació n de MATLAB (m-file), que está dividido en tres
m-files: Autor.m, cargas.m y estructura.m.

 Autor.m

disp('***UNIVERSIDAD NACIONAL DE SAN CRISTÓ BAL DE HUAMANGA***')

disp(' FACULTAD DE INGENIERÍA DE MINAS, GEOLOGÍA Y CIVIL')
disp(' ESCUELA DE FORMACIÓ N PROFESIONAL DE INGENIERÍA CIVIL')
disp(' ANÁ LISIS ESTRUCTURAL II')
disp(' (IC-444) ')
disp(' CÁ LCULO DE VECTORES DE CARGA')
disp('DOCENTE: Mg. Ing. YACHAPA CONDEÑ A, Rubén A.')
disp('ALUMNO :Arroyo Osorio , Jose Alberto.')
disp(' AYACUCHO - PERÚ ')
disp(' 2022')

 Cargas.m

%Arreglo Q_total
disp('INGRESO DE CARGAS PARA LA ESTRUCTURA');
disp('-------------------------------------');
fprintf ('\n DEFINICION DE LA MATRIZ DE CARGAS Q : \n');
Q=zeros(1,ngl);
puntual = zeros(3,nod);
V1=input('\n ¿Existen cargas en las juntas ? : ','s');
if V1 == 's'
njc = input ('\n ¿Cuantas juntas cargadas hay? :');
for i=1 :njc
if i==1
else
disp('siguiente junta cargada ');
end
NC = input (' \n Nº de la junta cargada :');
Q1(1)= input ('\n Fuerza horizontal : ');
puntual(1,NC)= Q1(1);
Q1(2)= input ('\n Fuerza vertical : ');
puntual(2,NC)= Q1(2);
Q1(3)= input ('\n Momento : ');
puntual(3,NC)= Q1(3);
VCJ(i,:)=CG(NC,:);
for m =1 :3
n=VCJ(i,m);
if Q1(m)~=0
Q(n)=Q1(m);
else end
end
end
Q_CJ = Q';
else
V1 = 'n' ;
Q_CJ = Q';
end
Q=zeros(1,ngl);
repartida = zeros(1,mbr);
Q2_almac (mbr,6) = zeros;
V3 = input ('\n ¿Existen cargas en los miembros?: ','s' );
if V3=='s'
nmc= input ('\n ¿Cuantos mienbros cargados hay?: ');
for ll=1:nmc

33
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

if ll==1
MC = input('\n Nº del primer miembro cargado :');
V4 = input('\n ¿El miembro tiene carga uniformemente distribuida?' ,'s');
else
MC = input ('\n Nº del siguiente miembro cargado :');
V4=input('\n ¿El miembro tiene carga uniformemente distribuida?' ,'s');
end
if V4=='s'
car= input ('\n Carga distribuida: ');
repartida(1,MC)=car;
Qa2(1)=0 ;
Qa2(4)=0 ;
Qa2(2)= car*L(MC)/2 ;
Qa2(5)=Qa2(2);
Qa2(3)= car*(L(MC)^2)/12;
Qa2(6)=-Qa2(3);
else
V4 = 'n';
Qa2(1)= input ('\n Fuerza axial N.I. : ');
Qa2(2)= input ('\n Fuerza cortante N.I. : ');
Qa2(3)= input ('\n Momento N.I. : ');
Qa2(4)= input ('\n Fuerza axial N.F. : ');
Qa2(5)= input ('\n Fuerza cortante N.F. : ');
Qa2(6)= input ('\n Momento N.F. : ');
end
Q2=Qa2';
for mm =1:6
Q2_almac(MC,mm)=Q2(mm)';
end
for m=1:6
for n=1:6
T2_3I(m,n)=T2_3(MC,m,n);
end
end
T2_3I;
clear Q3
Q3 = (-1)*T2_3I'*Q2 ;
for g =1 :6
h = VC(MC,g);
if h ~=0
Q(h)= Q3(g)+Q(h);
else end
end
end
Q_CM = Q' ;
else
V3='n';

Q_CM = Q' ;
end

%%%%%%%%%%%%%%%%repartida

fprintf( '\n Vector de Cargas Totales:');

Q = Q_CJ + Q_CM

 Estructura.m

% Arreglo CG. Coordenadas generales

disp('Ponga (s) para SI, (n)para NO:');
viga=input('¿Es solamente estructura de vigas? :','s');
nod=input('Numero de nudos totales :');
nnr=input(' Numero de apoyos :');

34
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

ngl=0;
CG=ones(nod,3);
% analisis de restricciones
for i=1:nnr
if i==1
disp(' Datos del primer apoyo ');
else
disp(' Datos del siguiente apoyo ');

end
nudres= input (' Numero del apoyo :');
X1 = input ('¿Se desplaza en X? :','s');
if X1=='n'
CG(i,1)=0;
else end
Y1 = input('¿Se desplaza en Y? :','s');
if Y1=='n'
CG(i,2)=0;
else end
R1 = input ('¿Tiene rotació n? :','s');
if R1=='n'
CG(i,3)=0;
else end
end
% grados de libertad
for i=1:nod
for j=1:3
if CG(i,j)~=0
ngl=ngl+1;
CG(i,j)=ngl;
else end
end
end
% % CG
% Arreglos vectoriales ini y fin. Nudos iniciales y finales.
mbr =input(' Numero de miembros:' );
for i=1:mbr
fprintf ('\n Miembro %d:',i);
ini(i)=input ('\n Numero de su nudo inicial:');
fin(i)=input ('\n Numero de su nudo final:');
end
ini;
fin;
% Arreglo VC. Vectores de colocacion
for i=1:mbr
for k=1:3
VC(i,k)= CG(ini(i), k);
VC(i,k+3) = CG(fin(i),k);
end
end
VC;
% Arreglo L, SENO y COSENO
fprintf ('\n Coordenadas para los nudos: \n');
for i=1:nod
fprintf ('\n Coordenada del Nudo %d: ',i);
x(i) = input ('\n En x [m]:');
y(i) = input ('\n En y [m]:');
end
for i=1:mbr
Dx(i) = x(fin(i)) -x(ini(i));
Dy(i) = y(fin(i)) -y(ini(i));
L(i) = ((Dx(i))^2 + (Dy(i))^2)^0.5;
SENO(i) = Dy(i)/L(i);
COSENO(i) = Dx(i)/L(i);
end
L;

35
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

SENO ;
COSENO ;
% Matriz K3
%fprintf ('\n Caracteristicas de los miembros : ');
%Elas= input ('\n Modulo de Elasticidad [T/m2] :');
Elas= 1;
G=0.4*Elas;
%Beta= input ('\n Factor de Forma :');
Beta=1;
aux=zeros(ngl,ngl);
for i=1 : mbr
%fprintf ('\n Miembro %d:',i);
B(i)=1;
%B(i) = input('\n Base [m]:');
H(i) =1;
%H(i) = input('\n Altura [m]:');
Area(i) = B(i)*H(i);
Inercia(i) = B(i)*H(i)^3/12;
fi = (3*Elas*Inercia(i)*Beta)/(G*Area(i)*L(i)^2);
c=((4*Elas*Inercia(i))/L(i))*((1+fi)/(1+4*fi));
cp=c;
a=((2*Elas*Inercia(i))/L(i))*((1-2*fi)/(1+4*fi));
b=((c+a)/L(i));
bp=b;
t=(b+bp)/L(i);
r=((Elas*Area(i))/L(i));
%Matriz de rigidez de miembro K2 en coordenadas locales
K2(i,1,1)=r; K2(i,1,2)=0; K2(i,1,3)=0; K2(i,1,4)=-r; K2(i,1,5)=0; K2(i,1,6)=0;
K2(i,2,1)=0; K2(i,2,2)=t; K2(i,2,3)=b; K2(i,2,4)=0; K2(i,2,5)=-t; K2(i,2,6)=bp;
K2(i,3,1)=0; K2(i,3,2)=b; K2(i,3,3)=c; K2(i,3,4)=0; K2(i,3,5)=-b; K2(i,3,6)=a;
K2(i,4,1)=-r; K2(i,4,2)=0; K2(i,4,3)=0; K2(i,4,4)=r; K2(i,4,5)=0; K2(i,4,6)=0;
K2(i,5,1)=0; K2(i,5,2)=-t; K2(i,5,3)=-b;K2(i,5,4)=0; K2(i,5,5)=t; K2(i,5,6)=-bp;
K2(i,6,1)=0; K2(i,6,2)=bp; K2(i,6,3)=a; K2(i,6,4)=0; K2(i,6,5)=-bp;K2(i,6,6)=cp;
%Matriz de paso T2-3
T2_3(i,1,1)=COSENO(i); T2_3(i,1,2)=SENO(i); T2_3(i,1,3)=0; T2_3(i,1,4)=0; T2_3(i,1,5)=0;
T2_3(i,1,6)=0;
T2_3(i,2,1)=-SENO(i); T2_3(i,2,2)=COSENO(i); T2_3(i,2,3)=0; T2_3(i,2,4)=0; T2_3(i,2,5)=0;
T2_3(i,2,6)=0;
T2_3(i,3,1)=0; T2_3(i, 3,2)=0; T2_3(i,3,3)=1; T2_3(i,3,4)=0; T2_3(i,3,5)=0;
T2_3(i,3,6)=0;
T2_3(i,4,1)=0; T2_3(i,4,2)=0; T2_3(i,4,3)=0; T2_3(i,4,4)=COSENO(i); T2_3(i,4,5)=SENO(i);
T2_3(i,4,6)=0;
T2_3(i,5,1)=0; T2_3(i,5,2)=0; T2_3(i,5,3)=0; T2_3(i,5,4)=-SENO(i);
T2_3(i,5,5)=COSENO(i);T2_3(i,5,6)=0;
T2_3(i,6,1)=0; T2_3(i,6,2)=0; T2_3(i,6,3)=0; T2_3(i,6,4)=0; T2_3(i,6,5)=0;
T2_3(i,6,6)=1;
%Matriz de rigidez de miembro K3 en coordenadas globales
for m=1:6
for n=1:6
K2I(m,n)=K2(i,m,n);
end
end
for m=1:6
for n=1:6
T2_3I(m,n)=T2_3(i,m,n);
end
end
fprintf ('\n Matriz de Rigidez en Coord. Globales del Miembro %d: ',i);
K3=T2_3I'*K2I*T2_3I
%Ensamblaje
mataux=K3;
for j=1 :6;
if VC(i,j)==0;
mataux(j,:)=0;
mataux(:,j)=0;
else end

36
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

end
k=zeros(ngl,ngl);
for cont1=1:6;
for cont2=1:6;
if mataux (cont1,cont2)~=0,
uno=VC(i,cont1);
dos=VC(i,cont2);
tres=mataux(cont1,cont2);
k(uno,dos)=k(uno,dos)+tres;
else end
end
end
if i==1;
aux=k;
else if i~=1,
aux=aux+k;
else end
end
end
fprintf ('\n Mat. de Rigidez de la Estructura (solo para los grados de libertad):');
Mat_rigidez = aux
cargas;

b) Para poder comprobar la solución del programa realizado en el lenguaje MATLAB (m-file),
realizaremos una solución analítica con un ejemplo aplicativo

Determinar el vector de cargas generalizadas Q, por medio del problema primario y complementario para
el marco plano de la siguiente figura cuyos elementos son totalmente flexibles.

Sistema Q – q

37
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

Problema Primario

Equilibrio de elementos

Equilibrio de nudos

R1 =0 R2 =−4 T R3 =2 .667 T −m R4 =0 R5 =0
R6 =0

Problema Complementario

38
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

[]
Por lo que la matriz Q es:

0 . 000
−4. 000
Q= 2 .667
0 . 000
0 . 000
0 .000
c) Comprobar la solución analítica con el programa realizado en el lenguaje MATLAB (m-file)

 Ejecutamos el programa “estructura”:

39
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

 Ingresamos los datos de la estructura: como el nú mero de nodos, apoyos, miembros y las respectivas
coordenadas de los nodos; así como también sus respectivas cargas ya sea en los miembros o juntas,
para este caso solo existen carga repartida en la barra. El usuario debe enumerar los nodos y las barras.

40
ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

41
 ANÁLISIS ESTRUCTURAL II (IC-444) INGENIERÍA CIVIL

 Podemos notar que la solució n del programa en MATLAB es la misma que la solució n analítica, por lo
que se puede concluir que el programa es perfecto.

42