Fem53 - HP50G
Fem53 - HP50G
Fem53 - HP50G
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Rigidezc
ANLISIS MATRICIAL DE ESTRUCTURAS
MTODO DE ELEMENTOS FINITOS
______________________________________________________________________
Ayacucho-2004
Pgina 1 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
INTRODUCCIN
Slo para comentar del archiconocido Fem49v5.3, lo que hace este
programa es impresionante en el anlisis estructural y me sorprende la
manera en que analiza las cargas en movimiento(lneas de influencia),
adems este programa en su versin final 5.3 es completamente
programable, con esto puedes alterar el normal funcionamiento y obtener
los resultados que creas conveniente.
El nuevo programa fem49v5.3 modificado, adems de hacer lo de
antes, ahora es capaz de mostrar lo siguiente:
-
Longitud, coseno del ngulo con respecto al eje horizontal, coseno del
ngulo con respecto al eje vertical, para cualquier elemento. Esto
para formar la matriz de transformacin de desplazamientos y matriz
de transformacin de fuerzas.
______________________________________________________________________
Ayacucho-2004
Pgina 2 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
______________________________________________________________________
Ayacucho-2004
Pgina 3 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
______________________________________________________________________
Ayacucho-2004
Pgina 4 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
______________________________________________________________________
Ayacucho-2004
Pgina 5 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
ANOTACIONES FINALES
comprobados
______________________________________________________________________
Ayacucho-2004
Pgina 6 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Ejemplos
Recuerda que antes de todo debes borrar la librera 1605
(fem49v5.3) de tu calculadora e instalar esta nueva versin adaptada.
Ejemplo #1.
Analizar la armadura mostrada en la figura.
Solucin:
Paso #1: la aplicacin del mtodo de la rigidez requiere subdividir la
estructura en una serie de elementos finitos e identificar sus puntos
extremos como nodos.
Lo primero es establecer un eje de coordenadas, enumerar nudos y
barras de nuestra estructura, que para este clculo debe ser como muestra
la figura contigua.
______________________________________________________________________
Ayacucho-2004
Pgina 7 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
nudo 1
nudo 2
nudo 3
nudo 4
X
0.
0.
10.
20.
Z
0.
5.
5.
5.
barra
barra
barra
barra
barra
Ni
1.
2.
3.
1.
1.
1
2
3
4
5
Nf
2.
3.
4.
3.
4.
propiedad
1.
1.
1.
1.
1.
soporte 1
nudo
1.
X
1.
Z
1.
y
0.
soporte 2
2.
1.
0.
0.
Carga 1
nudo
3.
Fx
0.
Fz
10.
M
0.
Carga 2
4.
0.
10.
0.
______________________________________________________________________
Ayacucho-2004
Pgina 8 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Donde en general:
______________________________________________________________________
Ayacucho-2004
Pgina 9 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Nudo
1
Nudo
2
Nudo
3
Nudo
4
1.2-
1X
1Z
2X
2Z
3X
3Z
4X
4Z
Pulsas
!Ay obtienes
______________________________________________________________________
Ayacucho-2004
Pgina 10 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
1.3-
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Del grfico:
Nudo inicial (Ni): 1
Nudo final (Nf): 3
Area: 1, E=1
Con estos datos, se tiene:
pulsa
y se obtiene
______________________________________________________________________
Ayacucho-2004
Pgina 11 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
nudo inicial
1X
1Z
.07155417528 .03577708764
.03577708764 .01788854382
-.07155417528 -.03577708764
-.03577708764 -.01788854382
nudo
3X
-.07155417528
-.03577708764
.07155417528
.03577708764
final
3Z
-.03577708764
-.01788854382
.03577708764
.01788854382
1X
1Z
3X
3Z
nudo inicial
nudo final
2.1.2- Para el elemento uno, Ni: 1, Nf: 2, A=1, E=1, como los elementos se
repiten slo se necesitan conocer sus cosenos directores de los ngulos
menores con que se orienta el elemento con respecto al sistema global.
1Z
2X
2Z
0.
0.
0.
0.
0.
.2
0.
-.2
0.
0.
0.
0.
0.
-.2
0.
.2
1X
1Z
2X
2Z
2Z
3X
3Z
.1
0.
-.1
0.
0.
0.
0.
0.
-.1
0.
.1
0.
0.
0.
0.
0.
2X
2Z
3X
3Z
2.1.4- Para el elemento tres, Ni: 3, Nf: 4., es el mismo que del elemento
dos, debido a que tiene las misma propiedad y la misma orientacin.
3X
3Z
4X
4Z
.1
0.
-.1
0.
0.
0.
0.
0.
-.1
0.
.1
0.
0.
0.
0.
0.
3X
3Z
4X
4Z
1Z
4X
4Z
4,57E+09
1,14E+09
-4,57E+09
1,14E+09
2,85E+08
-1,14E+09
-4,57E+09
-1,14E+09
4,57E+09
-1,14E+09
-2,85E+08
1,14E+09
1X
1Z
4X
______________________________________________________________________
Ayacucho-2004
Pgina 12 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
-1,14E+09
-2,85E+08
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
1,14E+09
4Z
2,85E+08
1X
1Z
2X
2Z
3X
3Z
4X
4Z
120.E-3
47.E-3
0.0E0
0.0E0
-72.E-3
-36.E-3
-46.E-3
-11.E-3
47.E-3
220.E-3
0.0E0
-200.E-3
-36.E-3
-18.E-3
-11.E-3
-2.9E-3
0.0E0
0.0E0
100.E-3
0.0E0
-100.E-3
0.0E0
0.0E0
0.0E0
0.0E0
-200.E-3
0.0E0
200.E-3
0.0E0
0.0E0
0.0E0
0.0E0
-72.E-3
-36.E-3
-100.E-3
0.0E0
270.E-3
36.E-3
-100.E-3
0.0E0
-36.E-3
-18.E-3
0.0E0
0.0E0
36.E-3
18.E-3
0.0E0
0.0E0
-46.E-3
-11.E-3
0.0E0
0.0E0
-100.E-3
0.0E0
150.E-3
11.E-3
-11.E-3
-2.9E-3
0.0E0
0.0E0
0.0E0
0.0E0
11.E-3
2.9E-3
1X
1Z
2X
2Z
3X
3Z
4X
4Z
1X
1Z
2X
2Z
3X
3Z
4X
4Z
100.E498
47.E-3
0.0E0
0.0E0
-72.E-3
-36.E-3
-46.E-3
-11.E-3
47.E-3
100.E498
0.0E0
-200.E-3
-36.E-3
-18.E-3
-11.E-3
-2.9E-3
0.0E0
0.0E0
100.E498
0.0E0
-100.E-3
0.0E0
0.0E0
0.0E0
0.0E0
-200.E-3
0.0E0
200.E-3
0.0E0
0.0E0
0.0E0
0.0E0
-72.E-3
-36.E-3
-100.E-3
0.0E0
270.E-3
36.E-3
-100.E-3
0.0E0
-36.E-3
-18.E-3
0.0E0
0.0E0
36.E-3
18.E-3
0.0E0
0.0E0
-46.E-3
-11.E-3
0.0E0
0.0E0
-100.E-3
0.0E0
150.E-3
11.E-3
-11.E-3
-2.9E-3
0.0E0
0.0E0
0.0E0
0.0E0
11.E-3
2.9E-3
1X
1Z
2X
2Z
3X
3Z
4X
4Z
Con los datos obtenidos hasta este punto ya es factible obtener los
desplazamientos en los nudos de la estructura3, ordenando entonces, de la
igualdad:
Para lograr ver los nmeros como en el cuadro la calculadora debe estar
configurado al modo ingeniera con un dgito de precisin.
2
las partes resaltados con azul representan las restricciones en los nudos donde se
restringe el desplazamiento.
1
______________________________________________________________________
Ayacucho-2004
Pgina 13 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
[D] = [Km F ]
Donde:
- [D]: vector de desplazamientos en el sistema global de la estructura.
- Km: matriz re rigidez procesado segn las restricciones en los
apoyos.
- F: vector de fuerzas equivalentes en los nudos.
Ordenando4, se obtiene:
1X
1Z
2X
2Z
3X
3Z
4X
4Z
Fuerza
100.E498
47.E-3
0.0E0
0.0E0
47.E-3
100.E498
0.0E0
-200.E-3
-72.E-3
-36.E-3
-46.E-3
-11.E-3
1X
-36.E-3
-18.E-3
-11.E-3
-2.9E-3
0.0E0
0.0E0
100.E498
1Z
0.0E0
-100.E-3
0.0E0
0.0E0
0.0E0
2X
0.0E0
-200.E-3
0.0E0
200.E-3
0.0E0
0.0E0
0.0E0
0.0E0
2Z
-72.E-3
-36.E-3
-100.E-3
0.0E0
270.E-3
36.E-3
-100.E-3
0.0E0
3X
-36.E-3
-18.E-3
0.0E0
0.0E0
36.E-3
18.E-3
0.0E0
0.0E0
10
3Z
-46.E-3
-11.E-3
0.0E0
0.0E0
-100.E-3
0.0E0
150.E-3
11.E-3
4X
-11.E-3
-2.9E-3
0.0E0
0.0E0
0.0E0
0.0E0
11.E-3
2.9E-3
10
4Z
1X
1Z
2X
2Z
3X
3Z
4X
4Z
1.0E0
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
6.0E-499
1X
0.0E0
1.0E0
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
2.0E-499
1Z
0.0E0
0.0E0
1.0E0
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
-6.0E-499
2X
0.0E0
0.0E0
0.0E0
1.0E0
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
2Z
0.0E0
0.0E0
0.0E0
0.0E0
1.0E0
0.0E0
0.0E0
0.0E0
-600.E0
3X
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
1.0E0
0.0E0
0.0E0
1.8E3
3Z
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
1.0E0
0.0E0
-1.0E3
4X
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
0.0E0
1.0E0
7.5E3
4Z
Nudo 1
Nudo 2
desplazamiento
6.0E-499
2.0E-499
-6.0E-499
0.0E0
direccin
1X
1Z
2X
2Z
______________________________________________________________________
Ayacucho-2004
Pgina 14 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Nudo 3
Nudo 4
-600.E0
1.8E3
-1.0E3
7.5E3
3X
3Z
4X
4Z
Esta seccin est orientada para obtener las fuerzas en los extremos
de los elementos con respecto a coordenadas locales.
De la igualdad.
Donde:
-
pulsa
Cy obtienes
.0894427191 - .0894427191
- .0894427191 .0894427191
Ordenando: k ' =
pulsa
Ordenando:
!Cy obtienes
0
0
.894427191 .4472135955
T =
0
0
.894427191 .4472135955
______________________________________________________________________
Ayacucho-2004
Pgina 15 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
pulsa
Cy obtienes
0
0
Ordenando: D =
600
1759
Multiplicando estos resultados obtenemos 6 la fuerza axial en los
extremos del elemento #4, se debe tener en la pila:
-
nivel 3: k4
nivel 2: T4
nivel 1: D4
***
y se
qi 22.36
=
Todo en coordenada local, es decir
q f 22.36
Ordenando el resultado.
Esta parte est orientado para obtener las fuerzas en los extremos de
las barras en coordenadas globales, con esto por superposicin fcilmente
se puede obtener las reacciones en las reacciones en los apoyos.
______________________________________________________________________
Ayacucho-2004
Pgina 16 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
[Q] = [T ].[q ]
T
Donde:
-
pulsa
Dy se obtiene
pulsa
!Dobtienes
______________________________________________________________________
Ayacucho-2004
Pgina 17 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
0
1
[Q] =
0
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
0
0 1 X
0 0 0 1Z
.
=
Representa las fuerzas en los extremos del
0 0 0 2 X
1
0 2 Z
pulsa
Dy se obtiene
pulsa
!Dobtienes
Ordenando y multiplicando.
0
0.89
20 1X
0.447
0 22.36 10 1Z
[Q] =
.
=
Que representa las fuerzas en los
0
0.89 22.36 20 3 X
0.447
0
10 3Z
extremos del elemento #4 en coordenada global.
4.3- Para el elemento #5.
4.3.1- matriz de transformacin de fuerzas.
pulsa
Dy se obtiene
pulsa
!Dobtienes
______________________________________________________________________
Ayacucho-2004
Pgina 18 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
0
0.97
40 1X
0.24
0 41.2 10 1Z
[Q] =
.
=
0
0.97 41.2 40 4 X
0.24
0
10 4 Z
Con estos argumentos simplemente por superposicin obtenemos las
reacciones en los apoyos, como para este ejemplo el nudo inicial para los
elementos #1, #4 y #5 es el nudo uno, simplemente queda sumar los
resultados las fuerzas obtenidas en el sistema global, as:
0 20 40 60 Direccin X
reaccin _ nudo _ uno = +
+
=
0 10 10 20 Direccin Z
Y el mismo procedimiento para hallar las otras fuerzas en el resto de
los soportes.
______________________________________________________________________
Ayacucho-2004
Pgina 19 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Ejemplo #2:
Analizar el prtico mostrado en la figura.
Solucin:
1- Para analizar la estructura mediante el mtodo de la rigidez es necesario
enumerar nudos y barras y sus respectivas orientaciones, de igual manera
optar un sistema de coordenada global. La misma que utiliza Fem49v5.3,
debe de quedar como en la figura.
______________________________________________________________________
Ayacucho-2004
Pgina 20 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
Direccin
X
0.
8.8
12.8
2.8
Direccin
Z
0.
0.
0.
6.
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Nudo
Nudo
Nudo
Nudo
1
2
3
4
2.2- Elementos.
Nudo
inicial
1.
2.
4.
nudo final
Propiedad
2.
3.
2.
1.
1.
2.
Elemento 1
Elemento 2
Elemento 3
2.3- Propiedad.
rea
.0075
.01
momento
inercia
Mdulo de
Young
Propiedad
1
propiedad
.0000083333 21000000.
2
.00000625
21000000.
2.4- Soportes.
Nudo
1.
4.
Condicin
en X
1.
1.
Condicin
en Z
1.
1.
Condicin
en Y
1.
1.
Apoyo 1
Apoyo 2
Fx
0.
Fz
20.
My
0.
d
4.
2.
0.
20.
0.
2.
Wz inicial
5.
Wz Final
5.
d inicial
0.
d final
0.
2.
5.
5.
0.
0.
______________________________________________________________________
Ayacucho-2004
Pgina 21 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
y el grfico
3- Una vez terminado de definir la estructura, ahora se analiza pulsando
SCALC.
4- Analizando paso a paso:
Segn la numeracin que se le dio al prtico se considera las
direcciones para cada nudo donde ocurre una accin o reaccin, en este
caso, tres libertades por nudo, el programa asigna automticamente estos
grados de libertad dependiendo a la numeracin de los nudos, entonces
para el procesamiento Rigidezc supone como queda en el grfico.
______________________________________________________________________
Ayacucho-2004
Pgina 22 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
4.1.1- Para obtener las cargas equivalentes reducidas a los nudos pulse
directamente Ay se obtiene9.
FUERZA
0.00
33.36
-56.07
1.11
81.72
62.10
0.00
20.00
16.67
-1.11
8.92
-13.38
DIRECCIN
1X
1Z
1Y
2X
2Z
2Y
3X
3Z
3Y
4X
4Z
4Y
Donde:
-
______________________________________________________________________
Ayacucho-2004
Pgina 23 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
10
______________________________________________________________________
Ayacucho-2004
Pgina 24 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
pulsa
mejor.
1X
17897.73
0.00
0.00
-17897.73
0.00
0.00
1Z
0.00
2.31
10.17
0.00
-2.31
10.17
By
1Y
0.00
10.17
59.66
0.00
-10.17
29.83
se obtiene
2X
-17897.73
0.00
0.00
17897.73
0.00
0.00
,
2Z
0.00
-2.31
-10.17
0.00
2.31
-10.17
2Y
0.00
10.17
29.83
0.00
-10.17
59.66
ordenndolo
1X
1Z
1Y
2X
2Z
2Y
______________________________________________________________________
Ayacucho-2004
Pgina 25 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
pulsa By se obtiene
2X
39375.00
0.00
0.00
-39375.00
0.00
0.00
2Z
0.00
24.61
49.22
0.00
-24.61
49.22
2Y
0.00
49.22
131.25
0.00
-49.22
65.63
3X
-39375.00
0.00
0.00
39375.00
0.00
0.00
ordedando.
3Z
0.00
-24.61
-49.22
0.00
24.61
-49.22
3Y
0.00
49.22
65.63
0.00
-49.22
131.25
2X
2Z
2Y
3X
3Z
3Y
pulsa By se obtiene
4X
12376.09
-12372.65
10.31
-12376.09
12372.65
10.31
4Z
-12372.65
12376.09
10.31
12372.65
-12376.09
10.31
4Y
10.31
10.31
82.50
-10.31
-10.31
41.25
2X
-12376.09
12372.65
-10.31
12376.09
-12372.65
-10.31
ordenando.
2Z
12372.65
-12376.09
-10.31
-12372.65
12376.09
-10.31
2Y
10.31
10.31
41.25
-10.31
-10.31
82.50
4X
4Z
4Y
2X
2Z
2Y
______________________________________________________________________
Ayacucho-2004
Pgina 26 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Ordenando.
1X
1Z
1Y
2X
2Z
2Y
3X
3Z
3Y
4X
4Z
17897.73
0.00
0.00
-17897.73
0.00
0.00
0.00
0.00
0.00
0.00
0.00
4Y
0.00
1X
0.00
2.31
-10.17
0.00
-2.31
-10.17
0.00
0.00
0.00
0.00
0.00
0.00
1Z
0.00
-10.17
59.66
0.00
10.17
29.83
0.00
0.00
0.00
0.00
0.00
0.00
1Y
-17897.73
0.00
0.00
69648.81
-12372.65
10.31
-39375.00
0.00
0.00
-12376.09
12372.65
10.31
2X
0.00
-2.31
10.17
-12372.65
12403.01
-28.74
0.00
-24.61
-49.22
12372.65
-12376.09
10.31
2Z
0.00
-10.17
29.83
10.31
-28.74
273.40
0.00
49.22
65.63
-10.31
-10.31
41.25
2Y
0.00
0.00
0.00
-39375.00
0.00
0.00
39375.00
0.00
0.00
0.00
0.00
0.00
3X
0.00
0.00
0.00
0.00
-24.61
49.22
0.00
24.61
49.22
0.00
0.00
0.00
3Z
0.00
0.00
0.00
0.00
-49.22
65.63
0.00
49.22
131.25
0.00
0.00
0.00
3Y
0.00
0.00
0.00
-12376.09
12372.65
-10.31
0.00
0.00
0.00
12376.09
-12372.65
-10.31
4X
0.00
0.00
0.00
12372.65
-12376.09
-10.31
0.00
0.00
0.00
-12372.65
12376.09
-10.31
4Z
0.00
0.00
0.00
10.31
10.31
41.25
0.00
0.00
0.00
-10.31
-10.31
82.50
4Y
______________________________________________________________________
Ayacucho-2004
Pgina 27 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
4.2.3- Para obtener la matriz de rigidez procesada de acuerdo a las restricciones (apoyos).
Pulsa Bse obtiene:
Ordenando11:
1X
1Z
1Y
2X
2Z
2Y
3X
3Z
3Y
4X
4Z
4Y
1.00E500
0.00
0.00
-17897.73
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
1X
0.00
1.00E500
-10.17
0.00
-2.31
-10.17
0.00
0.00
0.00
0.00
0.00
0.00
1Z
0.00
-10.17
1.00E500
0.00
10.17
29.83
0.00
0.00
0.00
0.00
0.00
0.00
1Y
-17897.73
0.00
0.00
69648.81
-12372.65
10.31
-39375.00
0.00
0.00
-12376.09
12372.65
10.31
2X
0.00
-2.31
10.17
-12372.65
12403.01
-28.74
0.00
-24.61
-49.22
12372.65
-12376.09
10.31
2Z
0.00
-10.17
29.83
10.31
-28.74
273.40
0.00
49.22
65.63
-10.31
-10.31
41.25
2Y
0.00
0.00
0.00
-39375.00
0.00
0.00
39375.00
0.00
0.00
0.00
0.00
0.00
3X
0.00
0.00
0.00
0.00
-24.61
49.22
0.00
24.61
49.22
0.00
0.00
0.00
3Z
0.00
0.00
0.00
0.00
-49.22
65.63
0.00
49.22
131.25
0.00
0.00
0.00
3Y
0.00
0.00
0.00
-12376.09
12372.65
-10.31
0.00
0.00
0.00
1.00E500
-12372.65
-10.31
4X
0.00
0.00
0.00
12372.65
-12376.09
-10.31
0.00
0.00
0.00
-12372.65
1.00E500
-10.31
4Z
0.00
0.00
0.00
10.31
10.31
41.25
0.00
0.00
0.00
-10.31
-10.31
1.00E500
4Y
11
______________________________________________________________________
Ayacucho-2004
Pgina 28 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Con estos resultados ya es posible obtener los desplazamientos en los apoyos, formando la matriz aumentada, para obtener la
siguiente igualdad
Donde:
- [D]: vector de desplazamientos en el sistema global de la estructura.
- Km: matriz re rigidez procesado segn las restricciones en los apoyos.
- F: vector de fuerzas equivalentes en los nudos.
1X
1Z
1Y
2X
2Z
2Y
3X
3Z
3Y
4X
4Z
4Y
Fuerza
1.00E500
0.00
0.00
-17897.73
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
1X
0.00
1.00E500
-10.17
0.00
-2.31
-10.17
0.00
0.00
0.00
0.00
0.00
0.00
33.36
1Z
0.00
-10.17
1.00E500
0.00
10.17
29.83
0.00
0.00
0.00
0.00
0.00
0.00
-56.07
1Y
-17897.73
0.00
0.00
69648.81
-12372.65
10.31
-39375.00
0.00
0.00
-12376.09
12372.65
10.31
1.11
2X
0.00
-2.31
10.17
-12372.65
12403.01
-28.74
0.00
-24.61
-49.22
12372.65
-12376.09
10.31
81.72
2Z
0.00
-10.17
29.83
10.31
-28.74
273.40
0.00
49.22
65.63
-10.31
-10.31
41.25
62.10
2Y
0.00
0.00
0.00
-39375.00
0.00
0.00
39375.00
0.00
0.00
0.00
0.00
0.00
0.00
3X
0.00
0.00
0.00
0.00
-24.61
49.22
0.00
24.61
49.22
0.00
0.00
0.00
20.00
3Z
0.00
0.00
0.00
0.00
-49.22
65.63
0.00
49.22
131.25
0.00
0.00
0.00
16.67
3Y
0.00
0.00
0.00
-12376.09
12372.65
-10.31
0.00
0.00
0.00
1.00E500
-12372.65
-10.31
-1.11
4X
0.00
0.00
0.00
12372.65
-12376.09
-10.31
0.00
0.00
0.00
-12372.65
1.00E500
-10.31
8.92
4Z
0.00
0.00
0.00
10.31
10.31
41.25
0.00
0.00
0.00
-10.31
-10.31
1.00E500
-13.38
4Y
______________________________________________________________________
Ayacucho-2004
Pgina 29 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Resolviendo el sistema de ecuaciones, se muestra la matriz reducida don la ltima columna representa los
Desplazamientos en los nudos.
1X
1Z
1Y
2X
2Z
2Y
3X
3Z
3Y
4X
4Z
4Y
desplazamiento
1.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
1.03E-498
0.00
1.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
3.33E-499
1Z
0.00
0.00
1.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
-5.59E-499
1Y
0.00
0.00
0.00
1.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.01
2X
0.00
0.00
0.00
0.00
1.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.01
2Z
0.00
0.00
0.00
0.00
0.00
1.00
0.00
0.00
0.00
0.00
0.00
0.00
-0.01
2Y
0.00
0.00
0.00
0.00
0.00
0.00
1.00
0.00
0.00
0.00
0.00
0.00
0.01
3X
0.00
0.00
0.00
0.00
0.00
0.00
0.00
1.00
0.00
0.00
0.00
0.00
2.29
3Z
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
1.00
0.00
0.00
0.00
-0.72
3Y
4X
1X
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
1.00
0.00
0.00
-1.03E-498
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
1.00
0.00
1.11E-498
4Z
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
1.00
-1.31E-499
4Y
______________________________________________________________________
Ayacucho-2004
Pgina 30 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
1.
desplazamiento desplazamiento
rotacin
en x
en Z
en Y
0.
0.
0.
2.
0.005759
.013992
3.
0.005759
2.293316
-0.722212
0.
0.
0.
nudo
4.
-0.011101
Donde:
-
pulsa
4X'
4Z'
4Y'
Cy se obtiene
2X'
2Z'
ordenando.
2Y'
24748.74
0.00
0.00
-24748.74
0.00
0.00
4X'
0.00
3.44
-14.58
0.00
-3.44
-14.58
4Z'
0.00
-14.58
82.50
0.00
14.58
41.25
4Y'
-24748.74
0.00
0.00
24748.74
0.00
0.00
2X'
0.00
-3.44
14.58
0.00
3.44
14.58
2Z'
0.00
-14.58
41.25
0.00
14.58
82.50
2Y'
______________________________________________________________________
Ayacucho-2004
Pgina 31 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
pulsa
!Cy se obtiene
Ordenando.
[T] =
0.71
-0.71
0.00
0.00
0.00
0.00
0.71
0.71
0.00
0.00
0.00
0.00
0.00
0.00
1.00
0.00
0.00
0.00
0.00
0.00
0.00
0.71
-0.71
0.00
0.00
0.00
0.00
0.71
0.71
0.00
0.00
0.00
0.00
0.00
0.00
1.00
pulsa
Cy se obtiene
______________________________________________________________________
Ayacucho-2004
Pgina 32 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Ordenando.
[D] =
-1.0307E-498
4X
1.1072E-498
4Z
-1.3128E-499
4Y
0.0058
2X
0.0140
2Z
-0.0111
2Y
Q4 x ' 151.17
V
4 z ' 5.4
M 4 y ' 13.13
Q2 x ' 129.96
V2 z ' 15.8
M 2 y ' 27.38
De esta manera se obtiene las fuerzas en los extremos de los
elementos para cada uno, ordenando en la tabla.
Elemento
Nudo
1.00
1.00
2.00
2.00
3.00
3.00
1.00
2.00
2.00
3.00
4.00
2.00
Fuerza
Fuerza
cortante
axial en x'
en z'
-103.07
-33.28
103.07
-30.72
0.00
-40.00
0.00
0.00
151.17
-5.41
-129.96
-15.80
Momento
en y'
55.88
-52.62
80.00
0.00
13.13
-27.38
Esta seccin est orientado para obtener las fuerzas en los extremos
de las barras en transformados a coordenada global, con esto por
superposicin se hallan las reacciones en los apoyos.
De la igualdad.
[Q] = [T ].[q ]
T
Donde:
______________________________________________________________________
Ayacucho-2004
Pgina 33 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
pulsa Dy se obtiene
matriz de transformacin de fuerzas.
que es la
4.4.2- Para obtener las fuerzas en los extremos de los elementos con
respecto a su sistema local.
pulsa
!Dy se obtiene
Ordenando.
[Q] =
FUERZAS
-103.07
-33.28
55.88
103.07
-30.72
-52.62
1X
1Z
1Y
2X
2Z
2Y
12
______________________________________________________________________
Ayacucho-2004
Pgina 34 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
Nudo
1.00
4.00
ESCUELA
ASIGNATURA
TEMA
AUTOR
Reaccin
en X
-103.07
103.07
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
Reaccin
en Z
-33.28
-110.72
Momento
en Y
55.88
13.13
______________________________________________________________________
Ayacucho-2004
Pgina 35 de 36
______________________________________________________________________
UNIVERSIDAD NACIONAL SAN CRISTBAL DE HUAMANGA
ESCUELA
ASIGNATURA
TEMA
AUTOR
: INGENIERA CIVIL
: ANLISIS ESTRUCTURAL II
: ANLISIS MATRICIAL DE ESTRUCTURAS
: CANCHARI GUTIRREZ, Edmundo.
AGRADCIMIENTOS:
IMPORTANTE:
El autor de este programa es Caspar Lugtmeier, slo la adaptacin
para que muestre estos resultados pertenecen al autor de este manual, el
programa adaptado es de libre distribucin como el original, no contiene
algn archivo que pueda alterar en normal funcionamientos del sistema de
su hp49g/hp49g+, si embargo yo, Canchari Gutirrez, Edmundo. No me
responsabilizo por los daos que pudiera ocasionar el uso de este programa,
ste es de libre distribucin y se distribuye como es.
El programa fue minuciosamente comprobado en un emulador con
rom 1.19-6 y 1.24 y en una calculadora real con versin de sistema
operativo 1.19-6 y 1.24 no habindose encontrado defectos.
COMENTARIOS Y SUGERENCIAS:
Todos sus comentarios
[email protected]
sugerencias
sern
bienvenidos
Visite:
http://cgedmundo.googlepages.com
______________________________________________________________________
Ayacucho-2004
Pgina 36 de 36