Año Académico 2022 – 2023

Carrera: Mecánica
Grupo: B
Profesor: Ing. Diego Nuñez, M.Sc.
INGENIERIA MECANICA

Problemas de Armaduras
Joel Leonardo Tixe Yugcha
27 de noviembre de 2022

Resumen

MATLAB

Los siguientes EJERCICIOS que se desarrollaran continuación son del libro de Daryl Logan donde nos con-
centramos él en EJEMPLO 3,11 donde al resolverlo se encontró 3 elementos de 3 nodos donde el nodo 1 se
encuentra con apoyo fijo mientras que el nodo 2 y 3 tienen apoyo móvil, la mayor dificultad que se presentó fue que
el nodo 3 tiene apoyo fijo de 45 grados, pero se encontró una solución favorable, para su resolución nos basaremos
en cada elemento tomando en cuenta la distancia y los grados de inclinación que tiene cada elemento. Mientras
tanto que en la software de ANSYS no se considerara el apoyo a los 45 grados sino a los 90 grados

ANSYS

Se consdera el el apoyo con inclinación como uno de 90grados en eje y y se resolvio de la siguiente manera
como se mostrara mas adelante debemos tener en consideracion que los resultados tendras una amplia variacion de
resultados como se muestra en la tabla final el porcentaje de error

MATLAB

The following EXERCISES that will be developed below are from the book by Daryl Logan where we focus
on EXAMPLE 3,11 where, when solving it, 3 elements of 3 nodes were found where node 1 has fixed support while
node 2 and 3 have mobile support, the greatest difficulty that arose was that the node 3 has a fixed support of 45
degrees, but a favorable solution was found, for its resolution we will base it on each element, taking into account
the distance and the degrees of inclination that each element has. Meanwhile, in the ANSYS software, support for
45 degrees will not be considered, but for 90 degrees.

ANSYS

The support with inclination is considered as one of 90 degrees on the y-axis and it was solved in the follo-
wing way, as will be shown later, we must take into consideration that the results will have a wide variation of
results as shown in the final table the percentage of error
Palabras clave: Desplazamiento, Fuerzas, esfuerzos Axiales, Longitudes, Matrices
2 Problemas de Armaduras

Índice

1. Problema N° 1 3
1.1. Modelo de analisis / Esquema de la estructura . . . . . . . . . . . . . . . . . . . . . . . 3
1.2. Resolución en Matlab . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.3. Resolución en Ansys APDL . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
1.4. Comparación de resultados Matlab vs. Ansys APDL . . . . . . . . . . . . . . . . . . . . 11

A. Anexos 12
A.1. Anexo 1: . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
A.2. Anexo 2: . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12

Año Académico 2022 – 2023

Joel Leonardo Tixe Yugcha 3

1. Problema N° 1

Libro:A First Course in the Finite Element Method (Fourth Edition)

Autor: Daryl L. Logan
Ejemplo 3.11

Para la armadura plana que se muestra en la figura 3-23, determine

√ los desplazamientos y las reacciones. 1
. Donde E = 210 GPa, A1 = 6.00x10-4 ; A2 = 6.00x10-4; A3 = 6 2x10-4 m2 y P = 1000 kN. Empezamos
usando la Ec. (3.4.23) para determinar la matriz de rigidez de cada elemento Determine:

1. Matriz rigidez global de la estructura

2. Los desplazamientos nodales

3. Las reacciones en los apoyos

4. Los esfuerzos en los elementos

5. El factor de seguridad mínimo

Figura 1: Estructura de 3 nodos con 1 apoyo inclinadol.

1.1. Modelo de analisis / Esquema de la estructura

La figura 2 muestra el esquema de la estructura planteada

Figura 2: Modelo de analisis de la estructura.

Año Académico 2022 – 2023

4 Problemas de Armaduras

1.2. Resolución en Matlab

A continuación se muestra los pasos necesarios para la solución de la estructura utlizando las funciones
disponibles en el siguiente link:

file:///D:/User/Downloads/sodapdf-converted.pdf

lo cual nos ayuda el el software MATLAB

1.2.1. Procedimiento

1.2.2. Matriz de conectividad

La matriz de conectividad de esta estructura es:

Elemento Nodo (𝑖) Nodo ( 𝑗)
1 1 2
2 2 3
3 1 3

1.2.3. Definicion de variables

E=210e11; % Módulo de elasticidad

A1 = 6.00x10-4; % Sección transversal
A2 = 6.00x10-4;
√ % Sección transversal
A3 = 6 2x10-4; % Sección transversal
L1=1.0; % Longitud elemento 1
L2=1.0; % Longitud elemento 2
L3=1.41; % Longitud elemento 3
theta1=90; % Angulo del elemento 1
theta2=0; % Angulo del elemento 2
theta3=45; % Angulo del elemento 3

1.2.4. Creación de matrices de cada elemento

k1 = PlaneTrussElementStiffness(E,A1,L1, theta1) ; % Creación matriz rigidez elemento 1

k2=PlaneTrussElementStiffness(E,A2,L2,theta2); % Creación matriz rigidez elemento 2
k3=PlaneTrussElementStiffness(E,A3,L3,theta3); % Creación matriz rigidez elemento 3
KT=k1+k2+k3 ; % Suma global de los 3 elementos

1.2.5. Ensamble de la matriz global

K=zeros(6,6); % Creación matriz de ceros

KT1=PlaneTrussAssemble(KT,k1,1,2) ; % Ensamblaje elemento 1
KT2=PlaneTrussAssemble(KT,k2,2,3) % Ensamblaje elemento 2
KT3=PlaneTrussAssemble(KT,k3,1,3) % Ensamblaje elemento 3
KTOTAL=KT1+KT2+KT3 % Ensamblaje Global de todos los elementos

1.2.6. Creación de matrices para la inclinación de apoyo de 45 grados

T =eye (6, 6); % Creación matriz de ceros para inclinaciones

Tprima=PlaneTrussInclinedSupport (T, 3, 45) ; % inserta en angulo del apoyo en el nodo 3

Año Académico 2022 – 2023

Joel Leonardo Tixe Yugcha 5

1.2.7. Inclinacion del angulo en el nodo 3

(Tprima)


 1 0 0 0 0 0  



 1 0 0 0 0 0  

0 1 0 0 0 0  0 1 0 0 0 0 

  
 

 
 
 

0

 0 1 0 0 0


 

 0 0 1 0 0 0 


=inv(Tprima)

 0 0 0 1 0 0   
 0 0 0 1 0 0  
0 0 0 0 0,7071 0,7071  0 0 0 0 0,7071 −0,7071 

  
 

 
 
 

   
0
 0 0 0 −0,7071 0,7071 

0
 0 0 0 0,7071 0,7071 

1.2.8. KTOTAL de los 3 elementos

(KTOTAL)


 0,6319 0,6319 −0 −0 −0,6319 −0,6319 


0,6319 1,8919 0 −1,2600 −0,6319 −0,6319 

 

 

 −0

 0 1,2600 0 −1,2600 0 


1.0e+08 *

 −0 1,2600 0 1,2600 0 0  
−0,6319 0,6319 −1,2600 0 1,8919 0,6319 

 

 


 −0,6319 
 0,6319 0 0 0,6319 0,6319 

1.2.9. Condiciones para el desarrollo de la formula K*=T x K x T𝑇

K*; % TGLOBAL
T; % Ttotal
K; % KTOTAL
T𝑇 ; % Inversa de Tprima

1.2.10. Aplicacion de la fórmula K*=T x K x T𝑇

TGLOBAL


 0,6319 0,6319 −0 −0 −0,8936 0  

0,6319 1,8919 0 −1,2600 −0,8936 0

 


 

 −0

 −0 1,2600 0 −0,8936 0,8910 


1.0e+08 *

 −0 −1,2600 0 1,2600 0 0  
−0,8936 −0,8936 −0,8910 0 1,8937 −0,6319

 


 

 
 0
 0 0,8910 0 −0,6319 0,6319 

1.2.11. Aplicacion de condiciones de borde

Como existe inclinación trabajamos con TGLOBAL La estructura se resuelve con la ecuacion:
  
𝐹 = 𝐾 (𝑇𝐺 𝐿𝑂𝐵𝐴𝐿) 𝑈 (1)

Aplicada a esta estructura queda:



 𝑓1 

 0,6319
 0,6319 −0 −0 −0,8936 0   
 𝑢1 

0,6319 1,8919 0 −1,2600 −0,8936 0
   


 𝑓 2






 

 𝑢 2 



 𝑓3 
 
  −0 −0 1,2600 0 −0,8936 0,8910  

 𝑢 3



= 1𝑒 + 8 
  ∗ (2)
𝑓4   −0 −1,2600 0 1,2600 0 0   𝑢4 
 
   
−0,8936 −0,8936 −0,8910 0 1,8937 −0,6319 
   
 𝑓5 



 

 
 𝑢5 

  
 𝑢 
𝑓 
  6
 0
 0 0,8910 0 −0,6319 0,6319    6


Año Académico 2022 – 2023

6 Problemas de Armaduras

Considerando las condiciones de borde 𝑢 1𝑥 = 0, 𝑢 1𝑦 = 0, 𝑢 2𝑦 = 0, 𝑢 3𝑦 = 0, 𝑓2𝑥 = 1000000𝑁, 𝑓3𝑥 = 0𝑁.



 𝑓1𝑥 

 0,6319
 0,6319 −0 −0 −0,8936 0  
 𝑢 1𝑥 


 0,6319 1,8919 0 −1,2600 −0,8936 0
  
  𝑢 1𝑦 


 𝑓 1𝑦


  





 1000000 

 
  −0 −0 1,2600 0 −0,8936 0,8910 
 𝑢 2𝑥 


= 1𝑒 + 8 
  (3)
𝑓2𝑦   −0 −1,2600 0 1,2600 0 0   𝑢 2𝑦 
 
  


 0


 −0,8936 −0,8936 −0,8910 0 1,8937 −0,6319


𝑢 3𝑥 


 
   
 

 𝑓
 3𝑦


 0
 0 0,8910 0 −0,6319 0,6319   𝑢 3𝑦 
 


Aplicando la eliminacion gaussiana de filas y columnas donde el desplazamiento es cero se extrae el

sistema a resolver:
    
1000000 1,2600 −0,8910 𝑢 2𝑥
= 1,0𝑒 + 08 ∗ (4)
0 −0,8910 1,8937 𝑢 3𝑥

K𝑒 = 𝑇𝐺 𝐿𝑂𝐵𝐴𝐿(3, 5); % Extracción de la matriz

f𝑒 = [1000000]; % Vector de fuerzas extraida
f𝑒 = [0]; % Vector de fuerzas extraida
u𝑒 = 𝑖𝑛𝑣(𝐾𝑇𝐸 ) ∗ 𝐹𝐸 ; % Vector de desplazamiento extraido a calcular

Resolviendo 𝑢 2𝑥 = 0,0119𝑚 ; 𝑢 3𝑥 = 0,0056𝑚

Luego se resuelve el sistema total para obtener las recciones en los apoyos :



 𝑓1𝑥 

 0,6319
 0,6319 −0 −0 −0,8936 0   
 0  

 0,6319 1,8919 0 −1,2600 −0,8936 0  0 
  


 𝑓 1𝑦


  





 1000000 

 
  −0 −0 1,2600 0 −0,8936 0,8910  0,0119 
 

= 1,0𝑒 + 08 
 
𝑓2𝑦   −0 −1,2600 0 1,2600 0 0   0 
 
  


 0


 −0,8936 −0,8936 −0,8910 0 1,8937 −0,6319 
 0,0056 



 
   
 

 𝑓
 3𝑦


 0
 0 0,8910 0 −0,6319 0,6319   0 
 

(5)

Para esto se completa el vector total de dezplazamientos 𝑈 y se resuelve la ecuación 1

UT=[0;0;u𝑒2𝑥 ; 0; 𝑢 𝑒3𝑥 : 0]; % Completar matriz total de desplazamientos

F=K(TGLOBAL)*UT; % Resolver ecuacion 1

Luego se obtiene 𝑓1𝑥 = −500𝑁, 𝑓1𝑦 = −500𝑁, 𝑓2𝑦 = 0𝑁, 𝑓3𝑦 = 707𝑁
Finalmente las fuerzas y los esfuerzos axiales de los elementos se calculan extrayendo los desplazamientos
de cada elemento y aplicando la funcion correspondiente:

Las fuerzas axiales fueron calculadas de la siguiente manera:

FT=TGLOBAL*UT


 𝑓1𝑥 = −500  

= −500
 


 𝑓 1𝑦



 𝑓1𝑦 = 1 

 

1.0e+06 *

 𝑓2𝑦 = 0 
 𝑓3𝑥 = 0 

 
 


 𝑓 = 707  
 3𝑦 

Año Académico 2022 – 2023

Joel Leonardo Tixe Yugcha 7

u1=[0 ; U(2)]; % Desplazamientos de elemento 1

s1=PlaneTrussElementStress(E,L1,theta1,U1) ; % Esfuerzo axial en el elemento 1
u2=[U(2);0]; % Desplazamientos de elemento 2
s2=PlaneTrussElementStress(E,L2,theta2,U2); % Esfuerzo axial en el elemento 2
u3=[U(2);0]; % Desplazamientos de elemento 2
s3=PlaneTrussElementStress(E,L3,theta3,U3); % Esfuerzo axial en el elemento 3

1.2.12. Resultados

A continuacion se resumen los resultados obtenidos:

1. Matriz rigidez global de la estructura

 0,6319
 0,6319 −0 −0 −0,8936 0 
 0,6319
 1,8919 0 −1,2600 −0,8936 0 
 −0 −0 1,2600 0 −0,8936 0,8910 
𝑇𝐺 𝐿𝑂𝐵𝐴𝐿 = 1,0𝑒 + 08 
 −0 −1,2600 0 1,2600 0 0 
−0,8936 −0,8936 −0,8910 0 1,8937 −0,6319

 0
 0 0,8910 0 −0,6319 0,6319 

2. Los desplazamientos nodales

𝑢 1 = 0𝑚, 𝑢 2 = 0𝑚, 𝑢 3 = 0,0119𝑚, 𝑢 4 = 0𝑚, 𝑢 5 = 0,0056𝑚, 𝑢 6 = 0𝑚

3. Las reacciones en los apoyos

𝑓1𝑥 = −500𝑁, 𝑓1𝑦 = −500𝑁, 𝑓2𝑥 = 1𝑁, 𝑓2𝑦 = 0𝑁, 𝑓3𝑥 = 0𝑁, 𝑓3𝑦 = 707𝑁

4. Los esfuerzos en los elementos

𝜎1 = 1,5293𝑒 − 07, 𝜎2 = 1,1750𝑒 + 09, 𝜎3 = −6,6323𝑒 + 08

5. El factor de seguridad mínimo

𝑛 = 1470

1.2.13. Código en Matlab depurado y comentado

A continuación se muestra el código continuo en Matlab haciendo uso de las funciones anexas:

Año Académico 2022 – 2023

8 Problemas de Armaduras

E=210e11; % Módulo de elasticidad

A1 = 6.00x10-4; % Sección transversal
A2 = 6.00x10-4;
√ % Sección transversal
A3 = 6 2x10-4; % Sección transversal
L1=1.0; % Longitud elemento 1
L2=1.0; % Longitud elemento 2
L3=1.41; % Longitud elemento 3
theta1=90; % Angulo del elemento 1
theta2=0; % Angulo del elemento 2
theta3=45; % Angulo del elemento 3
k1 = PlaneTrussElementStiffness(E,A1,L1, theta1) ; % Matriz rigidez elemento 1
k2=PlaneTrussElementStiffness(E,A2,L2,theta2) ; % Matriz rigidez elemento 2
k3=PlaneTrussElementStiffness(E,A3,L3,theta3) ; % Matriz rigidez elemento 3
KT=zeros(6,6) ; % Matriz 6x6
KT=k1+k2+k3 ; % Suma de Matrices de Rigidez para KT
T=eye (6, 6) ; % Matriz para el apoyo con inclinación de 45 grados
Tprima=PlaneTrussInclinedSupport (T, 3, 45) ; % Ubicacion e inclinacion del apoyo
KT=zeros(6,6) ; % Matriz para KT De cada elemento
KT1=PlaneTrussAssemble(KT,k1,1,2) ; % KT del elemento 1
KT2=PlaneTrussAssemble(KT,k2,2,3) ; % KT del elemento 2
KT3=PlaneTrussAssemble(KT,k3,1,3) ; % KT del elemento 3
KTOTAL=KT1+KT2+KT3 ; % Suma de las Matrices KT
Ttotal=Tprima*KTOTAL ; %Primer ensamble de la matriz
Tinversa=inv(Tprima) ; % Inversa de Tprima para el esamblaje
TGLOBAL=Ttotal*Tinversa ; % Ensamble completo de la matriz
KTE=TGLOBAL(3:5,3:5) ; % Extracción de la matriz (3.5)
FE=[1000000;0] ; % Aplicacioón de fuerzas en la matriz
UE1= inv(KTE)*(FE) ; % Solución de desplazamientos
FT=TGLOBAL*UT ; % Solución de Fuerzas
U1=[UT(1) ; UT(2) ; UT(3) ; UT(4)] ; % Desplazamientos 1
U2=[UT(1) ; UT(2) ; UT(5) ; UT(6)] ; % Desplazamientos 2
U3=[UT(3) ; UT(4) ; UT(5) ; UT(6)] ; % Desplazamientos 3
sigma1=PlaneTrussElementStress(E,L1,theta1,U1) ; % Esfuerzos axiales elemento 1
sigma2=PlaneTrussElementStress(E,L2,theta2,U2) ; % Esfuerzos axiales elemento 2
sigma3=PlaneTrussElementStress(E,L3,theta3,U3) ; % Esfuerzos axiales elemento 3

1.3. Resolución en Ansys APDL

En este apartado se muestra la solución de la estructura realizado en el software APDL

Año Académico 2022 – 2023

Joel Leonardo Tixe Yugcha 9

1.3.1. Código APDL depurado y comentado

/prep7 !* abre el preprocesador

mp, ex, 1, 210e9 !* Tipos de elementos
r,1,6e-4 !* define la sección trasversal Area 1
r,1,6e-4 !* define la sección trasversal Area 1
r,2,6e-4 !* define la sección trasversal Area 2
r,3,8.4850e-4 !* define la sección trasversal Area 3
n,1,0,0 !* creación nodo 1 coordenada X, Y
n,2,0,1 !* creación nodo 2 coordenada X, Y
n,3,1,1 !* creación nodo 3 coordenada X, Y
en,1,1,2 !* creación del elemento 1 coordenada X, Y
en,2,2,3 !* creación del elemento 2 coordenada X, Y
en,3,1,3 !* creación del elemento 3 coordenada X, Y
nlist !* Colocación del nodo en el elemento
elist !* Colocación del nodo en el elemento
CONDICIONES DE FRONTERA DE DESPLAZAMIENTO
d,1,ux,0 !* Aplicación del desplazamiento en el nodo 1 en direc
d,1,uy,0 !* Aplicación del desplazamiento en el nodo 1 en direc
d,2,uy,0 !* Aplicación del desplazamiento en el nodo 2 en direc
d,3,uy,0 !* Aplicación del desplazamiento en el nodo 3 en direc
F,2,Fx,1000000 !* Aplicación de fuerza puntual en el nodo 2
finish !* Cierre de preprocesado
PROCESADO
/solv !* Abre procesador
solve !* Resolución
finish !* Cierre de procesado
POSTPROCESADO
/post1 !* Abre POSTprocesado
prnsol,u !* Tabla de desplazamiento
plnsol,ux !* Imagen de desplazamiento
prrsol !* Imagen de fuerzas
finish !* Cierre de POSTprocesado

1.3.2. Resultados

A continuacion se resumen los resultados obtenidos en Ansys APDL:

1.3.3. Desplazamientos. Imagen y Tabla de valores

El desplazameinto total se muestra en la figura 9. En la figura 4 se muestra los valores de desplazamientos

nodales.
Por lo tanto, los desplazamientos nodales son:
𝑢 1 = 0𝑚, 𝑢 2 = 0,15192𝑚, 𝑢 3 = 0,11224𝑚

1.3.4. Reacciones en los apoyos. Tabla de valores

Las reacciones en los apoyos se muestran en la figura 5

Por lo tanto, las reacciones en los apoyos son son:
𝑓1𝑥 = −0,10000𝑒 + 007𝑁, 𝑓1𝑦 = −0,10000𝑒 + 007𝑁, 𝑓3𝑦 = −0,10000𝑒 + 007𝑁

Año Académico 2022 – 2023

10 Problemas de Armaduras

Figura 3: Desplazamiento total en m.

Figura 4: Desplazamiento en X m.

Figura 5: Reacciones en N.

1.3.5. Esfuerzos en los elementos. Tabla de valores

Los esfuerzos en los elementos se muestran en la figura ??

Año Académico 2022 – 2023

Joel Leonardo Tixe Yugcha 11

Por lo tanto, las reacciones en los apoyos son son:

𝜎1 = −1𝑀 𝑃𝑎, 𝜎2 = −1𝑀 𝑃𝑎, 𝜎3 = −1𝑀 𝑃𝑎

1.4. Comparación de resultados Matlab vs. Ansys APDL

Item Matlab Ansys % Error
u1 (𝑚) 0 0 0
u2 (𝑚) 0 0 0
u3 (𝑚) 0.0119 0.15192 11.766
u4 (𝑚) 0 0 0
u5 (𝑚) 0.0056 0.11224 19.035
u6 (𝑚) 0 0 0.00
f1𝑥 (𝑁) -500 1000 1
f1𝑦 (𝑁) -500 1000 1
f2𝑥 (𝑁) 1 0 0
f2𝑦 (𝑁) 0 0 0
f3𝑥 (𝑁) 0 0 0
f3𝑦 (𝑁) 707 1000 0.414
𝜎1 (𝑀 𝑃𝑎) 1.529e-7 1 -
𝜎2 (𝑀 𝑃𝑎) 1,1750e9 1 -
𝜎3 (𝑀 𝑃𝑎) -9.6323e8 1 -

Año Académico 2022 – 2023

12 Problemas de Armaduras

A. Anexos
A.1. Anexo 1:

Figura 6: Resultados de trabajo

A.2. Anexo 2:

Figura 7: movimiento de la estructura

Año Académico 2022 – 2023

Joel Leonardo Tixe Yugcha 13

Figura 8: Aplicación de elementos y nodos

Figura 9: Aplicacion de apoyos y fuerzas

Año Académico 2022 – 2023