Scribd
Cargar
85 vistas6 páginas

Análisis de Armadura en MathCad

Este documento presenta el análisis estructural de un ejemplo de armadura. Describe los materiales, secciones, cargas, condiciones de apoyo y coordenadas de los nodos. Luego calcula las matrices de rigidez de cada elemento y de la estructura completa, considerando las transformaciones geométricas. El objetivo es resolver la estructura y determinar sus desplazamientos y esfuerzos.

Cargado por

Diana Barros
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
85 vistas6 páginas

Análisis de Armadura en MathCad

Este documento presenta el análisis estructural de un ejemplo de armadura. Describe los materiales, secciones, cargas, condiciones de apoyo y coordenadas de los nodos. Luego calcula las matrices de rigidez de cada elemento y de la estructura completa, considerando las transformaciones geométricas. El objetivo es resolver la estructura y determinar sus desplazamientos y esfuerzos.

Cargado por

Diana Barros
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

Analisis armadura Ejemplo

Incluir << D:\Documentos\MathCad\Ejer1Deber4.mcdx

Materiales y Secciones
⎡ 1 2000 13 ⋅ 10 -6 ⎤ ⎡1 2 ⎤
MAT ≔ ⎢ -6 ⎥
PROP ≔ ⎢
⎣ 2 2000 13 ⋅ 10 ⎦ ⎣ 2 30 ⎥⎦

Cargas

Fuerzas Externas Nodales: si no hay cargas se debe colocar un nodo cualquiera y carga 0
⎡1 0 0⎤
⎢2 0 0⎥
FNod ≔ ⎢ ⎥
⎢3 0 0⎥
⎣4 0 0⎦
Fuerzas Externas No-Nodales

FnoNod #Carga w deltaT


⎡1 0.1 0 ⎤
⎢2 0 -50 ⎥
FNoNod ≔ ⎢ ⎥
⎢3 0 -30 ⎥
⎣4 0.1 50 ⎦

Condiciones de apoyo

Rest= Nodo GDL


Nodo GDL
⎡1 1 1⎤
Rest ≔ ⎢ 1 <--- Si esta restringido el movimiento
⎣ 4 1 1 ⎥⎦ 0 <--- Si esta libre

Coordenadas
⎡ 0 0 ⎤ ⎡ 0 0⎤
⎢ 200 200 ⎥ ⎢ 200 200 ⎥
XY ≔ ⎢ ⎥=⎢ ⎥
⎢ 400 200 ⎥ ⎢ 400 200 ⎥
⎣ 600 0 ⎦ ⎣ 600 0⎦

Diana Barros. Página 1


Matriz de Conectividad

CON= #Nodo NI NJ MAT SECCION FuerzaEmpotramiento

⎡1 1 2 1 1 1⎤
⎢2 2 3 1 1 2⎥ Escoger la fuerza de empotramiento de la base de datos
⎢ ⎥
CON ≔ ⎢ 3 3 4 1 1 0⎥ 0 <--- Sin fuerza de empotramiento
⎢4 1 3 1 1 2⎥
⎢⎣ 5 2 4 1 1 0 ⎥⎦

Resolucion de la estructura:

Calculo de longitudes de elementos.

⎡ 282.843 ⎤
⎢ 200 ⎥
⎢ ⎥
Long ≔ CalculoLongitudes ((XY , CON)) = ⎢ 282.843 ⎥
⎢ 447.214 ⎥
⎢⎣ 447.214 ⎥⎦

Matrices de transformacion

β ≔ MatrizBeta ((XY , CON , Long))


⎡ 0.707 0.707 0 0 ⎤ ⎡1 0 0 0⎤
⎢ -0.707 0.707 0 0 ⎥ ⎢0 1 0 0⎥
β =⎢ ⎥ β =⎢ ⎥
⎢ 0 0 0.707 0.707 ⎥ ⎢0 0 1 0⎥
1 2

⎣ 0 0 -0.707 0.707 ⎦ ⎣0 0 0 1⎦

⎡ 0.707 -0.707 0 0 ⎤ ⎡ 0.894 0.447 0 0 ⎤


⎢ 0.707 0.707 0 0 ⎥ ⎢ -0.447 0.894 0 0 ⎥
β =⎢ ⎥ β =⎢ ⎥
⎢0 0 0.707 -0.707 ⎥ ⎢ 0 0 0.894 0.447 ⎥
3 4

⎣0 0 0.707 0.707 ⎦ ⎣ 0 0 -0.447 0.894 ⎦

⎡ 0.894 -0.447 0 0 ⎤
⎢ 0.447 0.894 0 0 ⎥
β =⎢ ⎥
⎢0 0 0.894 -0.447 ⎥
5

⎣0 0 0.447 0.894 ⎦

Matriz de colocacion

L ≔ MatrizColocacion ((CON , XY))


⎡1 0 0 0 0 0 0 0⎤ ⎡0 0 1 0 0 0 0 0⎤ ⎡0 0 0 0 1 0 0 0⎤
⎢0 1 0 0 0 0 0 0⎥ ⎢0 0 0 1 0 0 0 0⎥ ⎢0 0 0 0 0 1 0 0⎥
L =⎢ ⎥ L =⎢ ⎥ L =⎢ ⎥
⎢0 0 1 0 0 0 0 0⎥ ⎢0 0 0 0 1 0 0 0⎥ ⎢0 0 0 0 0 0 1 0⎥
1 2 3

⎣0 0 0 1 0 0 0 0⎦ ⎣0 0 0 0 0 1 0 0⎦ ⎣0 0 0 0 0 0 0 1⎦

⎡1 0 0 0 0 0 0 0⎤ ⎡0 0 1 0 0 0 0 0⎤
⎢0 1 0 0 0 0 0 0⎥ ⎢0 0 0 1 0 0 0 0⎥
L =⎢ ⎥ L =⎢ ⎥
⎢0 0 0 0 1 0 0 0⎥ ⎢0 0 0 0 0 0 1 0⎥
4 5

⎣0 0 0 0 0 1 0 0⎦ ⎣0 0 0 0 0 0 0 1⎦

Diana Barros. Página 2


Matriz de rigidez de cada elemento en ejes globales y con dimensiones de la estructura

kele ≔ MatrizKelementos ((Long , β , L , CON , MAT , PROP))


⎡ 7.071 7.071 -7.071 -7.071 0 0 0 0 ⎤
⎢ 7.071 7.071 -7.071 -7.071 0 0 0 0 ⎥
⎢ ⎥
⎢ -7.071 -7.071 7.071 7.071 0 0 0 0 ⎥
-7.071 -7.071 7.071 7.071 0 0 0 0 ⎥
kele = ⎢
1 ⎢ 0 0 0 0 0 0 0 0⎥
⎢ 0 0 0 0 0 0 0 0⎥
⎢ ⎥
0 0 0 0 0 0 0 0
⎢ ⎥
⎣ 0 0 0 0 0 0 0 0⎦

⎡0 0 0 0 0 0 0 0⎤
⎢0 0 0 0 0 0 0 0⎥
⎢ ⎥
⎢0 0 20 0 -20 0 0 0 ⎥
0 0 0 0 0 0 0 0⎥
kele = ⎢
2 ⎢0 0 -20 0 20 0 0 0 ⎥
⎢0 0 0 0 0 0 0 0⎥
⎢ ⎥
0 0 0 0 0 0 0 0
⎢ ⎥
⎣0 0 0 0 0 0 0 0⎦

⎡0 0 0 0 0 0 0 0 ⎤
⎢0 0 0 0 0 0 0 0 ⎥
⎢ ⎥
⎢0 0 0 0 0 0 0 0 ⎥
0 0 0 0 0 0 0 0
kele = ⎢ ⎥
3 ⎢0 0 0 0 7.071 -7.071 -7.071 7.071 ⎥
⎢0 0 0 0 -7.071 7.071 7.071 -7.071 ⎥
⎢ ⎥
0 0 0 0 -7.071 7.071 7.071 -7.071
⎢ ⎥
⎣0 0 0 0 7.071 -7.071 -7.071 7.071 ⎦

⎡ 7.155 3.578 0 0 -7.155 -3.578 0 0 ⎤


⎢ 3.578 1.789 0 0 -3.578 -1.789 0 0 ⎥
⎢ ⎥
⎢ 0 0 0 0 0 0 0 0⎥
0 0 0 0 0 0 0 0⎥
kele = ⎢
4 ⎢ -7.155 -3.578 0 0 7.155 3.578 0 0 ⎥
⎢ -3.578 -1.789 0 0 3.578 1.789 0 0 ⎥
⎢ ⎥
0 0 0 0 0 0 0 0
⎢ ⎥
⎣ 0 0 0 0 0 0 0 0⎦

⎡0 0 0 0 0 0 0 0 ⎤
⎢0 0 0 0 0 0 0 0 ⎥
⎢ ⎥
⎢0 0 7.155 -3.578 0 0 -7.155 3.578 ⎥
0 0 -3.578 1.789 0 0 3.578 -1.789 ⎥
kele = ⎢
5 ⎢0 0 0 0 0 0 0 0 ⎥
⎢0 0 0 0 0 0 0 0 ⎥
⎢ ⎥
0 0 -7.155 3.578 0 0 7.155 -3.578
⎢ ⎥
⎣0 0 3.578 -1.789 0 0 -3.578 1.789 ⎦

Diana Barros. Página 3


Matriz de rigidez de la Estructura

⎡ 14.226 10.649 -7.071 -7.071 -7.155 -3.578 0 0 ⎤


⎢ 10.649 8.86 -7.071 -7.071 -3.578 -1.789 0 0 ⎥
⎢ ⎥
⎢ -7.071 -7.071 34.226 3.493 -20 0 -7.155 3.578 ⎥
-7.071 -7.071 3.493 8.86 0 0 3.578 -1.789 ⎥
K ≔ KGlob ((kele , CON)) = ⎢
⎢ -7.155 -3.578 -20 0 34.226 -3.493 -7.071 7.071 ⎥
⎢ -3.578 -1.789 0 0 -3.493 8.86 7.071 -7.071 ⎥
⎢ ⎥
0 0 -7.155 3.578 -7.071 7.071 14.226 -10.649
⎢ ⎥
⎣ 0 0 3.578 -1.789 7.071 -7.071 -10.649 8.86 ⎦

Fuerzas de empotramiento del elemento en ejes globales y con dimensiones


de la estructura

femp ≔ FEmp ⎛⎝CON , FNoNod , MAT , PROP , Long , β , L⎞⎠


⎡ -10 ⎤ ⎡ 0 ⎤ ⎡0⎤ ⎡ -2.326 ⎤ ⎡0⎤
⎢ -10 ⎥ ⎢ 0 ⎥ ⎢0⎥ ⎢ -1.163 ⎥ ⎢0⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎢ -10 ⎥ ⎢ -2.6 ⎥ ⎢0⎥ ⎢ 0 ⎥ ⎢0⎥
-10 ⎥ 0 ⎥ 0 0 0
femp = ⎢ femp = ⎢ femp = ⎢ ⎥ femp = ⎢ ⎥ femp = ⎢ ⎥
1 ⎢ 0⎥ 2 ⎢ 2.6 ⎥ 3 ⎢0⎥ 4 ⎢ 2.326 ⎥ 5 ⎢0⎥
⎢ 0⎥ ⎢ 0 ⎥ ⎢0⎥ ⎢ 1.163 ⎥ ⎢0⎥
⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥
0 0 0 0 0
⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥ ⎢ ⎥
⎣ 0⎦ ⎣ 0 ⎦ ⎣0⎦ ⎣ 0 ⎦ ⎣0⎦

Vector de fuerzas de empotramiento


⎡ -12.326 ⎤
⎢ -11.163 ⎥
⎢ ⎥
⎢ -12.6 ⎥
-10
Femp ≔ Fempotramiento ((femp , CON)) = ⎢ ⎥
⎢ 4.926 ⎥
⎢ 1.163 ⎥
⎢ ⎥
0
⎢ ⎥
⎣ 0 ⎦

Vector de fuerzas nodales


⎡0⎤
⎢0⎥
⎢ ⎥
⎢0⎥
0
F ≔ FuerzasNodales ⎛⎝FNod⎞⎠ = ⎢ ⎥
⎢0⎥
⎢0⎥
⎢ ⎥
0
⎢ ⎥
⎣0⎦

Diana Barros. Página 4


Vector de reacciones
⎡1⎤
⎢1⎥
⎢ ⎥
⎢0⎥
0
R ≔ VectorRestricciones ((Rest , XY)) = ⎢ ⎥
⎢0⎥
⎢0⎥
⎢ ⎥
1
⎢ ⎥
⎣1⎦

Particion de matrices
⎡3⎤ ⎡1⎤
⎢4⎥ ⎢2⎥
p ≔ match ((0 , R)) = ⎢ ⎥ s ≔ match ((1 , R)) = ⎢ ⎥
⎢5⎥ ⎢7⎥
⎣6⎦ ⎣8⎦

⎡ -7.071 -7.071 -7.155 3.578 ⎤


⎢ -7.071 -7.071 3.578 -1.789 ⎥
kps ≔ MatrizKps ((K , p , s)) = ⎢ ⎥
⎢ -7.155 -3.578 -7.071 7.071 ⎥
⎣ -3.578 -1.789 7.071 -7.071 ⎦

⎡ -7.071 -7.071 -7.155 -3.578 ⎤


⎢ -7.071 -7.071 -3.578 -1.789 ⎥
ksp ≔ MatrizKsp ((K , s , p)) = ⎢ ⎥
⎢ -7.155 3.578 -7.071 7.071 ⎥
⎣ 3.578 -1.789 7.071 -7.071 ⎦

⎡ 14.226 10.649 0 0 ⎤
⎢ 10.649 8.86 0 0 ⎥
kss ≔ MatrizKss ((K , s)) = ⎢ ⎥
⎢ 0 0 14.226 -10.649 ⎥
⎣ 0 0 -10.649 8.86 ⎦

⎡ 34.226 3.493 -20 0 ⎤


⎢ 3.493 8.86 0 0 ⎥
kpp ≔ MatrizKpp ((K , p)) = ⎢ ⎥
⎢ -20 0 34.226 -3.493 ⎥
⎣ 0 0 -3.493 8.86 ⎦

⎡0⎤
⎢0⎥
Fp ≔ MatrizFp ((F , p)) = ⎢ ⎥
⎢0⎥
⎣0⎦
⎡ -12.6 ⎤
⎢ -10 ⎥
Femp_p ≔ MatrizFemp_p ((Femp , p)) = ⎢ ⎥
⎢ 4.926 ⎥
⎣ 1.163 ⎦

⎡ -12.326 ⎤
⎢ -11.163 ⎥
Femp_s ≔ MatrizFemp_s ((Femp , s)) = ⎢ ⎥
⎢ 0 ⎥
⎣ 0 ⎦

Diana Barros. Página 5


Calculo de deslazamiento
⎡ 0.26 ⎤
⎢ 1.026 ⎥
Up ≔ kpp -1 ⋅ ((Fp - Femp_p)) = ⎢ ⎥
⎢ -0.005 ⎥
⎣ -0.133 ⎦

Vector desplazamiento
⎡ 0 ⎤
⎢ 0 ⎥
⎢ ⎥
⎢ 0.26 ⎥
1.026 ⎥
u ≔ VectorDesp ((R , Up)) = ⎢
⎢ -0.005 ⎥
⎢ -0.133 ⎥
⎢ ⎥
0
⎢ ⎥
⎣ 0 ⎦

Fuerzas Internas De Los Elementos:

f ≔ FuerzasInternasEGlob ((L , kele , u , femp))


⎡ -19.095 ⎤ ⎡ 2.714 ⎤ ⎡ 0.905 ⎤ ⎡ -1.809 ⎤ ⎡ -1.809 ⎤
⎢ -19.095 ⎥ ⎢ 0 ⎥ ⎢ -0.905 ⎥ ⎢ -0.905 ⎥ ⎢ 0.905 ⎥
f =⎢ ⎥ f =⎢ ⎥ f =⎢ ⎥ f =⎢ ⎥ f =⎢ ⎥
⎢ -0.905 ⎥ ⎢ -2.714 ⎥ ⎢ -0.905 ⎥ ⎢ 1.809 ⎥ ⎢ 1.809 ⎥
1 2 3 4 5

⎣ -0.905 ⎦ ⎣ 0 ⎦ ⎣ 0.905 ⎦ ⎣ 0.905 ⎦ ⎣ -0.905 ⎦

p ≔ FuerzasInternasELoc ((β , f))


⎡ -27.005 ⎤ ⎡ 2.714 ⎤ ⎡ 1.279 ⎤ ⎡ -2.023 ⎤ ⎡ -2.023 ⎤
⎢ 0 ⎥ ⎢ 0 ⎥ ⎢ 0 ⎥ ⎢ 0 ⎥ ⎢ 0 ⎥
p =⎢ ⎥ p =⎢ ⎥ p =⎢ ⎥ p =⎢ ⎥ p =⎢ ⎥
⎢ -1.279 ⎥ ⎢ -2.714 ⎥ ⎢ -1.279 ⎥ ⎢ 2.023 ⎥ ⎢ 2.023 ⎥
1 2 3 4 5

⎣ 0 ⎦ ⎣ 0 ⎦ ⎣ 0 ⎦ ⎣ 0 ⎦ ⎣ 0 ⎦

⎡ -20.90455 ⎤
⎢ -20 ⎥
FS ≔ ksp ⋅ Up + Femp_s = ⎢ ⎥
⎢ 0.90455 ⎥
⎣ 0 ⎦

Diana Barros. Página 6

También podría gustarte