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