Ing. Edgar E.

Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

1. INTRODUCCIÓN GENERAL .......................................................................... 3

2. FUNDAMENTOS DE ANALISIS Y DINAMICA ESTRUCTURAL ............ 3

2.1 ANALISIS ESTRUCTURAL ..................................................................................................... 3

2.1.1 Armaduras planas ................................................................................................................. 4
[Link] Definición y antecedentes ........................................................................................... 4
[Link] Métodos de análisis ............................................................................................................ 6
[Link].1 Método matricial ........................................................................................................ 7
[Link].1.1 Fundamentos ...................................................................................................... 7
[Link].1.2 Solución general por el método de desplazamientos .........................................11
[Link].1.3 Deducción de la matriz de rigidez del elemento armadura ................................13
[Link].1.4 Ejemplo: Solución simbólica de armadura con dos elementos..........................21
[Link].1.5 Ejemplo: Solución numérica de armadura de madera .......................................28
[Link].1.6 Aplicativo computacional para el análisis de armaduras planas ......................39
[Link].2 Métodos de energía ...................................................................................................47
[Link].2.1 Trabajo Virtual ..................................................................................................51
[Link].2.2 Método de Castigliano ......................................................................................55
2.1.2 Armaduras espaciales ..........................................................................................................59
[Link] Definición y antecedentes ..........................................................................................59
[Link] Método matricial ...............................................................................................................60
[Link].1 Fundamentos ............................................................................................................60
[Link].2 Ejemplo: Solución numérica de una torre en el espacio ..........................................63
[Link].3 Aplicativo computacional para el análisis de armaduras espaciales ....................67
2.1.3 Vigas ...................................................................................................................................72
[Link] Definición y antecedentes ..........................................................................................72
[Link] Métodos de análisis ...................................................................................................72
[Link].1 Pendiente deflexión o giro - deflexión ..................................................................72
[Link].1.1 Ejemplo: Solución simbólica de una viga continua ..........................................73
[Link].1.2 Ejemplo: Análisis de una vigueta que hace parte de un entrepiso .....................77
[Link].1.3 Ejemplo: Análisis de una viga de acero empotrada en los extremos .................84
[Link].2 Método de distribución de momentos ...................................................................93
[Link].3 Método matricial ...................................................................................................98
[Link].3.1 Fundamentos .....................................................................................................99
[Link].3.2 Ejemplo: Solución de una viga continúa de dos luces iguales ........................106
[Link].3.2 Ejemplo: Análisis estructural de un puente peatonal de dos luces ..................109
2.1.4 Pórtico plano .....................................................................................................................111
[Link] Definición y antecedentes ...............................................................................................111
[Link] Métodos de análisis .................................................................................................111
[Link].1 Método matricial .................................................................................................111
[Link].1.1 Fundamentos ..................................................................................................111
[Link].1.2 Ejemplo : Determinación de rigidez de un pórtico ........................................115
[Link].1.3 Ejemplo: Resolver un pórtico plano con carga horizontal y vertical .............122
[Link].1.3 Ejemplo: Resolver el pórtico plano con carga solamente vertical por giro -
deflexión ...........................................................................................................................128
[Link].1.4 Ejemplo: Resolver el pilón del puente por una carga vertical variable .........131
2.1.5 Parrilla o entramado ..........................................................................................................133

Edgar Eduardo Muñoz, IC, MIC 1

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

[Link] Definición y antecedentes ...............................................................................................133

[Link] Métodos de análisis .................................................................................................135
[Link].1 Método matricial .................................................................................................135
[Link].1.1 Fundamentos ...................................................................................................135
[Link].1.2 Ejemplo 1: Parrilla que debe soportar un motor ..............................................137
[Link].1.2 Ejemplo 2: Un balcón ......................................................................................142
2.1.6 Pórtico espacial .................................................................................................................148
[Link] Definición y antecedentes ...............................................................................................148
[Link] Métodos de análisis .................................................................................................148
[Link].1 Método matricial .................................................................................................148
[Link].1.1 Fundamentos ................................................................................................148
[Link].1.2 Ejemplo: Análisis de una pila en construcción ............................................152

3. REFERENCIAS .............................................................................................. 158

Edgar Eduardo Muñoz, IC, MIC 2

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

1. INTRODUCCIÓN GENERAL

2. FUNDAMENTOS DE ANALISIS Y DINAMICA ESTRUCTURAL

En este capítulo se presentan los fundamentos principales de la ingeniería estructural,

necesarios para el desarrollo y la solución de ejemplos reales de análisis y diseño de
estructuras.

2.1 ANALISIS ESTRUCTURAL

En la mayoría de los proyectos u obras civiles de ingeniería civil es necesaria la etapa

correspondiente al análisis estructural, que se constituye en una base para proceder
con el proceso de diseño estructural, elaboración de planos y finalmente construcción.
Esta fase o parte de los proyectos, consisten básicamente en la evaluación de las
acciones que actúan sobre la estructura, es decir de los diferentes tipos de cargas
externas, para determinar la correspondiente respuesta que se ocasionan en términos
de deformaciones y esfuerzos, para cada uno de sus componentes o partes como para
el conjunto en general. Para realizar el análisis estructural, son fundamentales los
conceptos de estática y mecánica de materiales, entre otras áreas, con lo cual se puede
originar datos esenciales y necesarios para el correspondiente diseño y posterior
construcción de cualquier estructura.

Así por ejemplo, para el diseño estructural de un puente vehicular, es necesario el

análisis estructural de su tablero, es decir la determinación de sus deformaciones y
esfuerzos, para los efectos de carga muerta y viva, y también de la cargas de viento,
sismo y empuje de tierras para el diseño de sus estribos, aletas y pilas. Por la
importancia de este tema en el área de estructuras, se presenta a continuación
diferentes métodos y ejemplos prácticos de análisis estructural para los tipos de
estructura normalmente utilizados en el país. Se presentan los siguientes tipos de
estructuras: armaduras planas, armaduras espaciales, vigas, pórticos planos, pórticos
espaciales y entramados o parrillas.

Edgar Eduardo Muñoz, IC, MIC 3

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

2.1.1 Armaduras planas

[Link] Definición y antecedentes

La armadura o cercha es un sistema estructural constituido por elementos

conformados principalmente por triángulos y sometidos generalmente a fuerzas
axiales de tensión o compresión, cuando deben soportar cargas externas aplicadas en
sus nudos y no sobre los elementos. Cuando tienen que asumir cargas externas
aplicadas sobre sus elementos, además de las cargas en los nudos, sus elementos no
solamente tienen fuerzas internas axiales si no también combinadas que pueden ser de
flexo compresión o flexo tensión. Son estructuras tradicionales e históricamente muy
importantes, que han sido útiles para la sociedad, empleadas para el diseño y la
construcción de diferentes partes o componentes de algunas obras civiles, tales como:
cubiertas de bodegas, cubiertas de iglesia, estadios, puentes, estructuras especiales,
entre otras. Es así, como por ejemplo en la Figura 2 - 1 y Figura 2 - 2 se presenta un
caso de un puente vehicular compuesto por armaduras de paso superior y otro un
´puente peatonal de la ciudad de Bogotá, respectivamente.

Figura 2 - 2 Puente peatonal en la ciudad de

Figura 2 - 1 Puente vehicular Puerto
Bogotá. Fuente: Propia
Salgar. Fuente: Muñoz et al [1]

También se presenta a continuación una foto de un ejemplo de una cubierta de un

coliseo construido en Colombia.

Edgar Eduardo Muñoz, IC, MIC 4

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 3 Cubierta de coliseo. Fuente: Por definir

Con respecto a su historia, según Muñoz et al [1], los primeros puentes de armadura o
celosía construidos en el mundo fueron en madera y esto fue gracias al Arquitecto
Italiano Adrea Palladio della Góndola (Padua, 1508 - Vicenza, 1580) quien en
1570(siglo XV) presentó en su libro denominado: «1 quattro libri dell 'Architettura”
con diferentes esquemas de estructuras en Celosía como se observa en la Figura 2 - 4
( Ver Galindo et al [3] y Gennaro et al [4]). Palladio trabajó fundamentalmente en la
ciudad de Venecia y alrededores, y se lo considera un autor del manierismo, sus obras
han influido de manera importante en la arquitectura del Neoclasicismo.

Figura 2 - 4 (a) Puente Brenta en Bassono – 1570 (b) Dibujo de “ Prima Inventione”
- 1570. Fuente: Creada a partir de la referencia [4 ]

Edgar Eduardo Muñoz, IC, MIC 5

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Nelson et. al [5], en su libro también menciona que fue el arquitecto Italiano Andrea
Palladio (1508-1580), quien primero usó las estructuras de armaduras, aunque no se
conoce la base de sus diseños. Es probable, según dicho autor, que haya revivido
antiguos diseños Romanos y dimensionó las armaduras por medio de reglas
empíricas. Su trabajo fue muy importante para la Arquitectura y la Ingeniería, sin
embargo después de su trabajo pasó alrededor de 200 años hasta que estas estructuras
se volvieron a emplear gracias al Suizo Ulric Grubermann.

Por otro lado, con los puentes de madera en celosía se lograron superar luces
importantes para la época, sin embargo este tipo de estructura no tenía la capacidad
estructural para solucionar luces mayores que fueron surgiendo en proyectos
ferroviarios o viales, para lo cual fue necesario en el siglo XVII la utilización del
acero, después de haber probado la construcción de este tipo de estructuras también
con materiales tales como: cáñamo, piedra y ladrillo.

[Link] Métodos de análisis

Se presentan a continuación algunos de los métodos, clásicos o modernos, utilizados

en diferentes periodos de la historia, con los cuales se puede realizar análisis
estructural de armaduras. En dicho cuadro se presenta para cada método los aspectos
dentro de un proceso de análisis estructural que se pueden determinar:

Aspectos que se pueden determinar

Método Deflexiones Reacciones Fuerzas internas
[Link] virtual X
2. Castigliano X
3. De los X X
nodos(Cuando es
estáticamente
determinada)
4. De las secciones X X
Cuando es
estáticamente
determinada)
5. Matricial X X X
Tabla 2 - 1. Algunos de los métodos de análisis estructural

A continuación se presenta los metoso1,2 y 5 sus fundamentos y ejemplos de

aplicación:

Edgar Eduardo Muñoz, IC, MIC 6

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

[Link].1 Método matricial

[Link].1.1 Fundamentos

El avance y el progreso de los computadores durante la década de 1945 y 1955,

produjo un desarrollo valioso del análisis de estructuras mediante el método matricial.
Este método tiene ventajas sobre los tradicionales, por su consistencia y rapidez para
la solución de problemas complejos de comportamiento estructural. Especialmente la
Industria aeronáutica contribuyó en la investigación de este tema, ya que necesitaban
soluciones y resultados rápidos del comportamiento estructural de los aviones a
reacción.

Uribe [6] en su libro, que ha sido la principal referencia nacional en este tema,
explica detalladamente este método cuyo desarrollo partió investigando la relación
entre las acciones que se producen en una estructura y la correspondiente respuesta.
Para entender este tema, en forma simple, se trae a este documento un ejemplo, que
consiste en una armadura simple proceso de análisis estructural se esquematiza a
continuación:

P1
P2

(a) Acción sobre la estructura (b) Acción sobre los elementos

(d) Respuesta de la estructura

(c) Respuesta de los elementos

Edgar Eduardo Muñoz, IC, MIC 7

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 5 Proceso de análisis en una armadura. Fuente: Adaptada de Uribe

[6]

Se observa que primero se deben identificar la acción que producen las cargas
externas, que hay en la práctica de diferentes tipos (Ver Figura 2 - 5 – (a)), las cuales
provocan fuerzas internas en los elementos como se observa en la Figura 2 - 5 – (b).
Es decir, para este caso cada elemento, por dicha acción, tiene que asumir y por lo
tanto soportar diferentes tipos de fuerzas axiales, las cuales dependen del armado o
tipología de la estructura. Es decir, como se muestra en esta figura, existe una “acción
sobre los elementos”, por lo cual un elemento de esta estructura debe estar preparado
para soportar una fuerza axial a compresión y otro elemento una fuerza axial a
tensión.

Esto produce una respuesta de cada uno de sus elementos y por lo tanto de la
estructura. Para los elementos, como se observa en la Figura 2 - 5 – (c), el que está
sometido a una fuerza de tensión sufre una deformación axial denominada:
alargamiento. El que está a compresión, también se deforma axialmente, lo que se
traduce en un acortamiento. Finalmente, así mismo como se producen deformaciones
en cada elemento, hay una respuesta de la estructura en conjunto, que depende de la
forma como colabora cada elementos, es decir cómo está armada, por lo cual se
deforma o sufre desplazamientos en los dos sentidos como se muestra en la Figura 2 -
5 – (d), que tienen relación con la dirección y magnitud de las cargas externas.

Esta relación, entre la acción y la respuesta, ha sido estudiada a los largo de mucho
tiempo por diferentes investigadores, logrando diferentes métodos, siendo los
primeros los “clásicos” y después los modernos, entre los que está el matricial. Este
último, remplaza la estructura continua por un modelo matemático de elementos
estructurales discretos y/o continuos dependiendo del tipo.

Con respecto a los métodos matriciales, se han estudiado dos métodos, uno de ellos
es el método de las “fuerzas”, que se expresa en la ecuación (2 - 1) y se basa en la
matriz de flexibilidad ([C]). El otro es el método de los “desplazamientos” (ver
ecuación (2 - 2), que se fundamenta es en la matriz de rigidez ([K]). Este último, es
el procedimiento más usado, de acuerdo con las referencias consultadas, por lo cual
es el que se emplea para el análisis estructural en el presente documento.

Edgar Eduardo Muñoz, IC, MIC 8

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Método de las fuerzas    C F  (2 - 1)

Método de los desplazamientos F   K   (2 - 2)

Dónde:

  Vector de deformaciones (desplazamientos o rotaciones)

F  Vector de fuerzas o momentos
[C] Matriz de Flexibilidad
[K] Matriz de rigidez

Por otro lado, el significado físico se puede expresar recordando los coeficientes de
influencia y el método de superposición, los cuales se estudian en los métodos
tradiciones de análisis estructural. Para entender el significado físico del método, se
presentan el caso de una viga simplemente apoyada con diferentes cargas externas,
que corresponden a las acciones y los desplazamientos verticales que corresponden a
parte de su respuesta (Ver Figura 2 - 6).

Edgar Eduardo Muñoz, IC, MIC 9

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

f1 f2 f3

  
f1

c11 c21 c31

f2

c12 c22 c32

f3

c13 c23 c33

Figura 2 - 6 Ejemplo del significado físico. Fuente: Adaptada de Uribe [6]

La solución estructural de este caso se puede expresar a través de un número de

ecuaciones lineales, que son iguales al número de incógnitas, de la siguiente forma:

1  c11 f1  c12 f 2  c13 f3 (2 - 3)

 2  c21 f1  c22 f 2  c23 f3 (2 - 4)

Edgar Eduardo Muñoz, IC, MIC 10

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 3  c31 f1  c32 f 2  c33 f3 (2 - 5)

Donde los coeficientes de influencia son necesarios para la solución y representan la

relación entre la acción y la reacción. Se convierte entonces, en el caso de un sistema
lineal de ecuaciones, que matricialmente se pueden expresar como:

1  C11 C12 C13   f1  (2 - 6)

    
 2   C21 C22 C23   f 2 
  C C33   f 3 
 3   31 C32

Este sistema representa el método de las fuerzas e incluye la matriz de flexibilidad,

que tienen especial uso en la dinámica estructural. Despejando el vector de fuerzas,
se obtiene otro sistema, que corresponde al método de “desplazamientos”, donde la
matriz de flexibilidad es la inversa de la matriz de rigidez.

 P1  C11 C12 C13   1  (2 - 7)

    
 2   C 21 C 22 C 23   2 
  C C33   3 
 3   31 C32

[Link].1.2 Solución general por el método de desplazamientos

Para la solución general por el método de los desplazamientos, se parte de la ecuación

(2 - 2 ] y se establecen unos sub-vectores y sub matrices de acuerdo con las variables
que son conocidas y desconocidas. Utilizando la misma terminología y procedimiento
la referencia Uribe [6], se tiene:

 Fn   K nn K na   n   Fnf  (2 - 8)
    
 Fa   K an K aa   a   Faf 

En donde:

Edgar Eduardo Muñoz, IC, MIC 11

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Fn  Es el vector de las cargas aplicadas en los nudos (variables conocidas). Su

tamaño depende del número de grados de libertad libres; es decir de la
configuración de la estructuras (tipos y número de apoyos) y del tipo de
estructura (armaduras, vigas, pórticos, etc).
Fa  Corresponde a las reacciones en los apoyos, en los grados de libertad
restringidos. Su tamaño depende del número y tipo de apoyos(variables
desconocidas)
 n  Es el vector de desplazamientos de los nudos libres. (desconocidos)
 a  Corresponde a los desplazamientos en los apoyos. Generalmente son iguales
a 0. Cuando hay asentamientos son diferentes de 0.
F 
n
f Vector de fuerzas de empotramiento, sobre los grados de libertad libres. Es la
forma de considerar en el sistema matricial, las cargas aplicadas en los
elementos(conocidos)
F 
a
f Vector de fuerzas de empotramiento, sobre los grados de libertad restringidos
(conocidos)

Para el caso de estructuras que no tengan cargas repartidas en los elementos y que no
existan asentamientos, la solución del sistema se expresa como:

 n   K nn 1 Fn  (2 - 9)

Fa   K an  n  (2 - 10)

En el caso que existan cargas repartidas en los elementos, la solución está dada por:

 n   K nn 1 Fn  

F 
n
f
n
(2 - 11)

Fa   K an  n   Faf  (2 - 12)

Las fuerzas se expresan de la forma Fix, donde el primer subíndice indica el número
del nudo y el segundo su dirección de acuerdo con el sistema de coordenadas
globales. Los momentos se identifican con la letra M y sus subíndices cumplen la
misma convención anterior. Los desplazamientos y rotaciones con las letras U y θ

Edgar Eduardo Muñoz, IC, MIC 12

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

respectivamente. Los vectores fila o columna se representan con el corchete { } y las

matrices con [ ]. Cuando se refiera a {F} o {}, son los vectores totales de fuerzas y
deformaciones respectivamente. Los sub-vectores se representan con un subíndice y
pueden ser de fuerzas y desplazamientos, tales como: {Fn} , {Fa}, {  n} y {  a}.

[Link].1.3 Deducción de la matriz de rigidez del elemento armadura

Se presenta a continuación la deducción de la matriz de rigidez local y general del

elemento tipo armadura, que específicamente va a estar sometido a cargas axiales.
Para ellos, se parte de la Figura 2 - 7, donde se dibuja un elemento prismático sobre
el eje X, del cual se conoce: área transversal, longitud y módulo de elasticidad.

Y
Datos:
Área (A)
Coordenadas Longitud (L)
globales

i j X

Figura 2 - 7 Elemento prismático en el eje X. Fuente: Propia

Es un elemento de dos nudos y cada nudo tiene dos grados de libertad, es decir que
cada uno puede tener desplazamientos en sentido X y Y. No tienen ninguna rotación
porque no incluyen para este caso ni flexo compresión ni flexo compresión. La
relación entre la acción y la respuesta para solamente este elemento se puede expresar
matricialmente de la siguiente forma:

 f ix   K11 K12   uix 

   
K 22  u jx 
(2 - 13)
 f jx   K 21

Para deducir los valores K11,K12,K21, K22 que hacen parte de la matriz de rigidez, es
necesario emplear el concepto básicos de la Mecánica de Sólidos o Materiales, que
relaciona la acción que son las fuerzas externas con respecto a la respuesta que es la
fuerza axial.

Edgar Eduardo Muñoz, IC, MIC 13

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Para ello se parte del diagrama que gráfica la relación entre los esfuerzos y la
deformación unitaria de un material. En dicho gráfico en la parte solamente elástica,
que es una recta, se conoce que su pendiente corresponde al módulo de elasticidad ( E
). Por ser una recta, se sabe que la relación entre el módulo de elasticidad, la
deformación unitaria y el esfuerzo axial es la siguiente:

 (2 - 14)
E


También que la relación que existe entre la fuerza axial (P) y el esfuerzo axial,
depende del área (A) y es la siguiente:

P (2 - 15)

A

Remplazando la ecuación (2 - 14) en la ecuación (2 - 15) se tiene:

P (2 - 16)

AE

Por otro lado, como la deformación unitaria es igual a la deformación sobre la

longitud, se tiene:

PL (2 - 17)

AE

Edgar Eduardo Muñoz, IC, MIC 14

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Para la deducción de los términos de la matriz de rigidez del elemento de armadura,

se obliga a la estructura a obtener una posición deformada unitaria en uno de los
nudos y los demás iguales a 0. Se comienza, considerando que la deformación Uix sea
unitaria, es decir igual a 1 y la Uiy sea nula, o sea igual a 0, como se observa en la
siguiente figura:
.

i j

fix uix=1.0 fiy

Figura 2 - 8 Caso 1: Uix= 1 y Uiy = 0 . Fuente: Propia

Se resuelve este caso en el sistema mostrado en la ecuación (2 - 13), obteniéndose lo

siguiente:

 f ix   K11 K12   u ix  1   f ix   K11 

     
 f jx   K 21 K 22  u jx  0  f jx   K 21 

Aplicando la ecuación (2 - 17) , se obtienen los siguientes resultados:

AE
f1x   K11
L

AE
f jx    K 21
L
Es decir, que se han determinado los primero dos coeficientes de rigidez, lo cual
representa:

Edgar Eduardo Muñoz, IC, MIC 15

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 AE 
 f ix   L K 12   u 
    AE
ix
 u 
 f jx   K 22   jx 
 L 

Posteriormente se considera que la deformación Uix sea nula, es decir igual a 0 y la

Uiy sea unitaria, o sea igual a 1, como se observa en la siguiente figura:
.

i j

fix fiy
ujx=1.0

Figura 2 - 9 Caso 2: Uix= 0 y Uiy =1 . Fuente: Propia

Se resuelve este caso en el sistema mostrado en la ecuación (2 - 13), obteniéndose lo

siguiente:

 f ix   K11 K12  u ix  0  f ix   K12 

     
 f jx   K 21 K 22  u jx  1  f jx   K 22 

Aplicando la ecuación (2 - 17), se obtienen los siguientes resultados:

AE
f1x    K12
L

AE
f jx   K 22
L

Edgar Eduardo Muñoz, IC, MIC 16

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Es decir, el sistema con la matriz de rigidez local de la estructura, queda como:

 AE AE 
 f ix   L 
    AE L   u ix 
 f jx   AE  u jx 

 L L 

Donde la matriz de rigidez local ( Klocal), con un solo grado de libertad por nudo,
por el momento es:

 AE AE  (2 - 18)
 
 L L 
AE 
K local
AE
 
 L L 

Para deducir un elemento arbitrariamente orientado, es necesario evaluar la matriz de

rigidez global (kglobal) la cual se basa en la matriz de rigidez local antes deducida y
una matriz de transformación. Para determinar la matriz de transformación, es
necesario buscar la relación que existe entre las fuerzas en coordenadas locales y las
fuerzas en las coordenadas globales. Para ello se presenta una figura donde se
esquematiza dicha relación:

y fix

fiy

X
Nudo i

Edgar Eduardo Muñoz, IC, MIC 17

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

x (b)
Y
fjx

y j

i 

fiy X

fix
(a)

Figura 2 - 10. Relación entre fuerzas locales y globales. (a) Elemento con fuerzas
locales. (b) Nudo “i” con fuerzas locales. Nota: En minúscula son las fuerzas locales
y el sistema de coordenadas. En mayúscula es el sistema de coordenadas globales.
Fuente: Propia

Por trigonometría, se tienen las siguientes ecuaciones:

F1x  Cos( ) f1x  Sin( ) f iy

F1 y  Sin( ) f1x  Cos( ) f iy

Que se pueden expresar, como:

 Fix  cos  sin   f ix 

   
 Fiy   sin cos   f iy 

Dónde:

cos  sin  
T1  
 sin  cos 

Edgar Eduardo Muñoz, IC, MIC 18

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Como la matriz “T1” es ortogonal, es decir que su inversa es igual a la transpuesta, se

puede concluir que:

 f ix   cos sin   Fix 

   
 f iy   sin cos   Fiy 

Expandiendo a los dos(2) nudos se tiene:

 fix   cos sin  0 0   Fix 
    
 f iy   sin  cos  0 0   Fiy 
f   
 jx   0 0 cos  sin    F jx 
 f jy   0 
  0  sin  cos   F jy 

Lo que representa que la matriz de transformación es:

 cos sin  0 0 
 sin  cos 0 0 
T  (2 - 19)
 0 0 cos sin  
 
 0 0  sin  cos 

Ahora se debe encontrar una relación entre la matriz de rigidez local (Klocal) de cada
elemento y la matriz de transformación antes determinada, para poder conocer la
matriz de rigidez global (Kglobal). Se comienza por la relación entre el vector de
fuerzas globales con respecto a las fuerzas locales, empleando la matriz de
transformación inversa que es igual a la matriz de transformación transpuesta, como
se expresa a continuación:

F   T 1  f   T T  f 
(2 - 20)

Edgar Eduardo Muñoz, IC, MIC 19

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

La misma relación entre las fuerzas locales y globales, es válida también entre las
deformaciones locales y globales, por lo tanto:

(2 - 21)
 local   T  global 

Por otro lado la relación entre las fuerzas en coordenadas locales y las deformaciones
en coordenadas locales, depende de la matriz de rigidez local del elemento, por lo
cual:

 f   k Local  Local 

(2 - 22)

Remplazando (2 - 22) en (2 - 20), se tiene:

(2 - 23)

F  T T k Local  Local 

Remplazando (2 - 21) en (2 - 23), y simplificando, se tiene:

(2 - 24)
K global   T  K T 
T
local

Con la anterior ecuación se hace el triple producto, con lo cual se obtiene la matriz de
rigidez global del elemento prismático arbitrariamente orientado:

Edgar Eduardo Muñoz, IC, MIC 20

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 Cos 2 CosSin  Cos 2  CosSin  (2 - 25)

 
AE  CosSin Sin2  CosSin  Sin2 
K  
L   Cos 2  CosSin Cos 2 CosSin 
global

 
 CosSin  Sin2 CosSin Sin2 

Basados en la solución de la matriz de rigidez global de cada elemento y la solución

del método de los desplazamientos, se puede resolver las armaduras planas
estáticamente determinadas e indeterminadas.

[Link].1.4 Ejemplo: Solución simbólica de armadura con dos elementos

Para entender como es el procedimiento para la solución por el método de

desplazamientos, se presenta a continuación un ejemplo simbólico del caso de una
armadura plana, compuesta por dos (2) elementos y tres(3) nudos.

Procedimiento Ejemplo simbólico

a) Materiales y sección Datos:
transversal: Para cada E : Módulo de elasticidad
elemento de la estructura se A: Área
debe establecer el módulo de
elasticidad y el área. El
predimensionamiento se basa
en el criterio del ingeniero
diseñador y de algunas
recomendaciones o sugerencias
que puede tener las normas.

Edgar Eduardo Muñoz, IC, MIC 21

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

b) Sistema estático: Contiene la

forma y tipificación de la
estructura; que incluye el
número de nudos, número de
elementos y su dirección
(incidencias). Además el
número y tipo de apoyos, que
son parte de las condiciones de
frontera. De esta forma se
establecen la cantidad de
grados de libertad y grados
restringidos.
Figura 2 - 11. Esquema de la estructura.
Fuente: Propia

Número de nudos : 3
Número de elementos :2
Número de apoyos:2
Apoyo 1: segundo orden.
Apoyo 2: segundo orden
Número de grados de libertad : 2
Número de grados restringidos: 4

c) Evaluación y aplicación de las En la Figura 2 - 11 se observa una carga

cargas: Consiste en la vertical y otra horizontal aplicada en el nudo 1,
evaluación de las diferentes que pueden ser generadas por una combinación
tipos de cargas sobre la de carga muerta más viento.
estructura (muerta, viva, viento,
etc) y la forma de su aplicación.
Esto incluye las combinaciones
de carga necesarias para el
diseño y para estudiar la
seguridad de la misma.
Para este caso, se tienen dos elementos. Los
d) Matriz de rigidez global de dos tienen un área “A”, un módulo de
cada elemento: Con base en la elasticidad de “E” y una longitud equivalente
conformación de los elementos de : L/ 2 .El ángulo del elemento 2-1 es de 45
en la estructura, sus grados y el de 3-1 es de 135 grados. La matriz
propiedades y su dirección se de rigidez global de cada elemento es la
determina la matriz de rigidez siguiente:
global de cada uno basados en

Edgar Eduardo Muñoz, IC, MIC 22

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

la ecuación (2 - 25).

 0.5 0.5  0.5 0.5 

 0.5 0.5  0.5  0.5
k 21 
AE 2 
L  0.5  0.5 0.5 0.5 
 
 0.5  0.5 0.5 0.5 

 0.5  0.5  0.5 0.5 

 0.5  0.5
AE 2  0.5 0.5
k 31 
L  0.5 0.5 0.5  0.5
 
 0.5  0.5  0.5 0.5 
e) Solución del sistema por el  
 
método de los  
desplazamientos. Se plante el  
    
sistema de acuerdo con la  
  
ecuación (2 - 8).     
 F1x  P1   
 2 0    U1 x ?
   
 F1 y   P2  AE  0 2    U1 y  ?
 F ?   
  U 2 x  0 
 2x  L  2 2
 
 F2 y  ?   2 2   U 2 y  0
    
 F3 x  ?    U 3 x  0
2 2
 
 F3 y  ?   2 2   U 3 y  0 
   2 2 
 2 2 
 
 2 
2 
 2 2 

1
f) Determinación de  2 AE   2 PL 
deformaciones. Se resuelven U1x   L 0   P1   
       2 AE 
aplicando la ecuación (2 - 9) y U1 y   0 2 AE   P2   2 PL 

 L   2 AE 
con ello se determinan las
deformaciones verticales y
horizontales que tiene el nudo
2.

Edgar Eduardo Muñoz, IC, MIC 23

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

g) Determinación de las
reacciones. Para evaluar las
 2 2  P 2 P1 
reacciones sobre la estructura     2  2 
se aplica la ecuación (2 - 10).  F2 x   2 2 
F   2 2   2 PL   P 2 P1 
 2 y  AE  2
     
2   2 AE    2 2 
   P 2 P1 
 F3 x  L  2 2  2 PL    
  
 F3 y   2 2   2 AE   2
 2
 2 2  P 2 P1
 2   2  2 
2 

h) Verificación de equilibrio. Inicialmente se verifica que la sumatoria de las

Determinadas las fuerzas en sentido “Y” sean igual a 0.
deformaciones del nudo 1 y las
reacciones de los nudos 2 y 3,  Fy  0
es necesario verificar P 2 P1 P 2 P1
 P2      0 _ OK
2 2 2 2
equilibrio, con lo cual se
.comprueba si el análisis fue
realizado correctamente.
También se verifica que la sumatoria de las
fuerzas en sentido “X” sean igual a 0.

 Fx  0
P2 P1 P2 P1
P1      0 _ OK
2 2 2 2
Finalmente se comprueba que la sumatoria de
momentos con respecto al nudo 2 sean igual a
0.
M 0
L L  P 2 P1 L
P1( )  P 2( )    ( )  0 _ OK
2 2  2 2  2

Adicionalmente, resolver la anterior armadura cambiando la condición de las fuerzas

externas, como se muestra en la siguiente figura:

Edgar Eduardo Muñoz, IC, MIC 24

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 12. Esquema de la estructura con otro sistema de cargas externas. Fuente:
Propia

Al cambiarse solamente las condiciones de carga externa, la matriz de rigidez de

cada elemento y de la estructura son los mismos. Pr ello, cambian uno de los sus
vectores del sistema de solución, como se observa a continuación:

 2 0 
 
 F1 x  P1   0 2  U  ? 
1x 1
F  0  2 
U1 y  ? 
2
 1y   2

2   

 2x
F  ? 
 AE  2 2 
 0 

     
F
 2x  ?  L  2 2  0 
 F2 x  ?   2 2  0 
   

 2x
F  ? 
  2 2  
 0



 2

2

 2 2 

Para determinar las deformaciones, de acuerdo al sistema explicado anteriormente se

tienen dos ecuaciones con dos incógnitas:

AE
P1  2U1x  0 * U1 y
L
AE
P2  0U1x  2 *U1 y
L

De lo anterior, se tiene que:

Edgar Eduardo Muñoz, IC, MIC 25

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

LP1
U 1x 
2 AE
U1y  0

Determinadas las deformaciones, se pueden determinar las reacciones:

 2 2
  
 F2 x   2 2 
 0.5P1 
F   2 
2
 LP1   0.5P 
 2 y  AE  2 2   1
  2 AE    
 F3 x  L  2 2   0   0.5P1 
  
 F3 y   2 2   0.5P1 
 2 2
 2 
2 

Se revisa equilibrio, en “y”:

 Fy  0
0  P1  P1  0 _ OK

Se revisa equilibrio, en “x”:

F X 0
0  P1  0.5P1  0.5P1  0 _ OK

Finalmente se comprueba que la sumatoria de momentos con respecto al nudo 2 sean

igual a 0:

M 0
L
 P1( )  0.5 P1 * L  0 _ OK
2

Ahora se determinan las fuerzas internas, dónde se termina que el elemento 2-1 está a
tensión y el elemento 31 a compresión:

Edgar Eduardo Muñoz, IC, MIC 26

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 LP1 
  0
S 21 
2 AE
0.707 0.707 2 AE   0.707P1
( L)  0 
 LP1 
  0
S31 
2 AE
 0707 0.707 2 AE   0.707P1
( L)  0 

Por otro lado, determinar cuál debe ser el valor de la carga P2, de la siguiente figura,
para que la deformación vertical del nudo 1 sea del orden de -0.07m y la horizontal de
este mismo nudo sea nula.

Figura 2 - 13. Esquema de la estructura con otro sistema de cargas externas. Fuente:
Propia

Considerar los siguientes datos:

 E=200000000 kN/m2
 L=1.5 m
 A=0.0007m2

Para ello, se revalúa e sistema de solución:

Edgar Eduardo Muñoz, IC, MIC 27

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 2 0 
 
 F1 x  ?  0 2  0 
F  2 
 ? 
2
 1y   2

2 
 0.07 

 F2 x  ?
 AE  2 2 
 0 

     
 F2 x  ? L  2 2  0 
 F2 x  ?  2 2  0 
   
  ?  2 2  




 F2 x  0
 2

2

 2 2 

Se tienen dos ecuaciones con dos incógnitas:

AE
P1  2 (0)  0 * ( 0.07)  0
L
AE
P2  0U 1x  2 * ( 0.07)  9239.5kN
L

[Link].1.5 Ejemplo: Solución numérica de armadura de madera

El siguiente ejercicio corresponde al análisis matricial de una cercha, que hace parte
de una ampliación de una iglesia localizada en Bogotá. La edificación existente tiene
un lote de 20mx10m y se desea aprovechar un jardín que tiene al lado y tiene de
frente 5 metros como se observa Figura 2 - 14.

Edgar Eduardo Muñoz, IC, MIC 28

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

5m

Correas existentes
5m

5m Correas nuevas

5m

10 m 5m

Planta del lote

(Sin escala)
Iglesia
existente
Ampliación
propuesta

5m

Perfil
(Sin escala)

Figura 2 - 14. Planta y perfil de caso de armadura en madera para iglesia. Fuente:
Propia

El objetivo es realizar el análisis estructural de una de las cerchas de ampliación que

tiene la siguiente geometría y condiciones de frontera. Igualmente se presenta con sus
correspondientes nudos y elementos debidamente orientados.

Edgar Eduardo Muñoz, IC, MIC 29

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

2
1.34 m

4 3

2.50 m 2.50 m

Figura 2 - 15. Planta y perfil de caso de armadura en madera para iglesia. Fuente:
Propia

a) Materiales

Se propone usar madera caimito (clasificación dentro del grupo B según NSR-98 y el
manual de maderas del Grupo Andino). El módulo de elasticidad de estas maderas
según NSR-10 [9] es de 10000 Mpa (promedio).

b) Evaluación de cargas

b.1) Carga muerta

Edgar Eduardo Muñoz, IC, MIC 30

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

La carga muerta corresponde al peso propio de la cercha, correa, tejas y los demás
elementos que conformaran la cubierta (Ver Tabla 2 - 2):

Elemento Carga Carga (Kgf/m2)

(kN/m2)
Peso propio 0.20 20
cercha(considerado)
Teja de barro más 0.80 80
mortero
Entablado 0.20 20
Impermeabilizante 0.10 10
Correas(considerado) 0.05 5
Total 1.35 135
Tabla 2 - 2. Evaluación de carga muerta de la cubierta

b.2) Carga Viva

Para cubiertas inclinadas de estructuras metálicas y madera la NSR-10 [9]

recomienda la siguiente carga viva, la cual depende básicamente de la pendiente

Pendiente mayor al 20% = 0.35 KN/m2=35 Kgf/m2

Pendiente menor al 20% = 0.50 KN/m2=50 Kgf/m2

La armadura en estudio tiene una pendiente del 26.8%, luego la carga viva que se
utiliza es de 0.35 kN//m2.

b.3) Carga de viento

De acuerdo con NSR-10 [9], se deben evaluar los siguientes aspectos para la
evaluación del viento:

Edgar Eduardo Muñoz, IC, MIC 31

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 Velocidad del viento básico (V)

Según el mapa de amenaza eólica, la velocidad del viento para Bogotá es de 80 Kph.

 Velocidad del viento de diseño

Vs = V(S1)(S2)(S3) =51.2 Kph =14.2 m/seg

Coeficiente topografía S1 = 1.0

Coeficiente Rugosidad Clase A S2 = 0.64 (tabla B.6.5-2)
Coeficiente Grado seguridad y vida útil S3 = 1.0
Variación de la densidad del aire con la altura S4 = 0.73

 Presión del viento

Q = .0625 Vs2*S4 = 9.19 Kgf/m2 =0.0919 kN/m2

 Aplicación del viento

Cp (barlovento)=-0.70
Cp (sotavento)=-0.50

Distribución y aplicación de las cargas en los nudos de la cercha

Para carga muerta:

Pd2 = 1.35*2.58*5.0 = 17.41 kN

Pd1 = 1.35*(2.58/2)*5.0 = 8.17 kN

Para carga viva:

Pl2 = .35*2.58*5.0 = 4.51 kN

Pl1 = .35*(2.58/2)*5.0 = 2.26 kN

Para viento:

Pw2 = 0.0919*2.58*5.0*0.7 = -83 kgf (succión)

Pw1 = 0.0919*2.58*2.5*0.7 = -41.5 kgf (succión)

Edgar Eduardo Muñoz, IC, MIC 32

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

(b)
(a)

(c)

Figura 2 - 16. Cargas puntuales por carga muerta(a), viva(b) y viento(c) en la

armadura.
b.4) Combinaciones de carga

Las combinaciones dependen del material y del método de diseño que vaya a utilizar.
Para este caso el método tradicional de diseño es en esfuerzos admisibles:

(1) D
(2) D+L
(3) D+W
(4) D+0.7E
(5) D+L+W
(6) D+L+0.7E

Edgar Eduardo Muñoz, IC, MIC 33

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Donde:

D: Carga muerta
L: Carga viva
W: Carga de viento
E: Carga de sismo

C) Análisis matricial

C.1) Matriz de rigidez de cada elemento

La cercha está compuesta por cuatro (4) nudos y cinco (5) elementos, cuyas
características se presentan a continuación:

X(m) Y(m) Rx Ry
1 0 1.34 0 0
2 2.5 0.67 0 0
3 5.0 0 0 1
4 0 0 1 1

Tabla 2 - [Link]ía y restricciones de la cercha

No. Sección Ni Nf A(m2) L(m) θ AE/L

1 2-4x19 3.0 4.0 0.0152 5 180.0 30400.0

Edgar Eduardo Muñoz, IC, MIC 34

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

2 2-4x19 3.0 2.0 0.0152 2.58 165.0 58727.5

3 2-4x19 2.0 1.0 0.0152 2.85 165.0 58727.5
4 1-4x9 4 1 0.0036 1.34 90 26865.7
5 1-9x9 4.0 2.0 0.0081 2.58 15.0 31295.6
Tabla 2 - 4. Incidencia de los elementos y propiedades

Se determinan las matrices de rigidez de cada uno de los cinco elementos empleando
la ecuación (2 - 25).

Matriz de rigidez elemento 3-4

 F3 x   30400 0  30400 0. u 3 x 

F    
 3y    0 0 0 0  u 3 y 
 F4 x   30400 0 30400 0  u 4 x 
    
 F4 y   0 0 0 0  u 4 y 

Matriz de rigidez de los elementos 3-2 y 2-1

 F3 x   54792  14684  54792 14684  u 3 x 

F    
 3 y     14684 3935 14684  3935  u 3 y 
 F2 x   54792 14684 54792  14684 u 2 x 
    
 F2 y   14684  3935  14684 3935  u 2 y 

 F2 x   54792  14684  54792 14684  u 2 x 

F    
 2 y     14684 3935 14684  3935  u 2 y 
 F1x   54792 14684 54792  14684  u1x 
    
 F1 y   14684  3935  14684 3935   u1 y 

Edgar Eduardo Muñoz, IC, MIC 35

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Matriz de rigidez elemento 4-1

 F4 x  0 0 0 0  u 4 x 
F    
 4 y   0 26865.70 0  26865.70 u 4 y 
 F1x  0 0 0 0   u1x 
    
 F1 y  0  26865.70 0 26865.70   u1 y 

Matriz de rigidez elemento 4-2

 F4 x   29198.4 7825.2  29198.4  7825.2  u 4 x 

F    
 4 y    7825.20 2097.2  7825.2  2097.20 u 4 y 
 F2 x   29198.40  7825.20 29198.40 7825.20  u 2 x 
    
 F2 y    7825.20  2097.20 7825.20 2097.20  u 2 y 

C.2) Solución del sistema por el método de los desplazamientos

Se plante el sistema de acuerdo con la ecuación (2 - 8) y para carga muerta.

 F1x  0   U 1x  ? 
 F1y  8.71   54792.4  
  
- 14684.1 - 54792.4 14684 0 -  U 1 y  ? 
 F 2 x  0   - 14684.1 30800.9 
14684.0 - 3935.2 0 -  U 2 x  
   
 F 2 y  17.41 - 54792.4 14684 138783.1 - 21542.7 - 54792.4 -  U 2 y  

 F 3x  0   14684 - 3935.2 - 21542.7 9967.5 14684 -  U 3x  
   
 F3y  ?   0 0 - 54792.4 14684 85192.4 - U 3 y  0 

 F 4x  ?   0 0 14684  3935.2  14683.7 - U 4 x  0 
    
 F4y  ?   0 1  29198.4  7825.2  30400 - U 4 y  0 

    268655.7  7825.2  2097.2 2 -  
   1
 
   

Se determinan las deformaciones en los nudos:

Edgar Eduardo Muñoz, IC, MIC 36

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 Ux1  0.0015m 
 Uy1  0.0003242m 
 
 n    Ux2  0.0002744m 
 Uy2  0.005174m 
 
Ux3  0.00106847m

Se hallan las reacciones para este caso de carga:

 F3y  8.71kN 
Fa   F4x  0 
F4y  17.41kN 
 
Se revisó equilibrio, es decir sumatoria de las fuerzas en el eje y se encontró que son
iguales a o, por lo cual el procedimiento usado es correcto. Ahora es necesario
determinar las fuerzas internas en los elementos, la cual se evalúa con la siguiente
expresión:

Ujx  Uix 
Sij 
AE
Cos Sin  
Ujy  Uiy 
(2 - 26)
L

Edgar Eduardo Muñoz, IC, MIC 37

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Por lo tanto, la fuerza interna de cada uno de los elementos para la hipótesis de carga
muerta es:

0  0.0010684
S 34  30400Cos (180) Sin(180)   32.48kN
 0 

0.0002744  0.0010684
S 32  58727.5Cos (165) Sin(165)   33.63kN
  0.005174  0 

 0.00157  0.0002744
S 21  58727.5Cos (165) Sin(165)   0kN
 0.00032  0.005174

 0.00157  0 
S 41  26865.6Cos (90) Sin(90)   8.71kN
 0.00032  0

 0.00027  0 
S 42  31295.6Cos (15) Sin(15)   33.63kN
 0.00517  0

A continuación se presentan los resultados para las diferentes combinaciones de

carga:

Edgar Eduardo Muñoz, IC, MIC 38

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Deformaciones(m)
D L W
Ux1 0.00157439 0.00040876 -5.5929E-05
Uy1 -0.00032421 -8.4122E-05 1.5797E-05
Ux2 0.00027442 7.1245E-05 -3.5968E-06
Uy2 -0.00517483 -0.00134349 0.00021829
Ux3 0.00106847 0.0002774 -3.8694E-05
Reaciones(kN)
Fy3 8.705 2.26 -0.344
Fx4 0 0 -0.427
Fy4 17.415 4.52 -0.854
Fuerzas internas en elementos(kN)
`3-4 32.481 8.43 -1.176
`3-2 -33.628 -8.73 1.328
`2-1 0 0.00 -0.11
`4-1 -8.71 -2.26 0.424
`4-2 -33.628 -8.73 1.66

Tabla 2 - 5. Resultados de deformaciones, reacciones y fuerzas internas del análisis

estructural de la superestructura.

[Link].1.6 Aplicativo computacional para el análisis de armaduras planas

Para resolver armaduras planas se desarrolló un programa en Visual Basic para Excel
que permite resolver cualquier estructura de este tipo para una condición de carga
dada. Se presenta a continuación su diagrama de flujo, que está constituido por dos
partes: pre procesos y post proceso. Con respecto a pre proceso se tiene:

Edgar Eduardo Muñoz, IC, MIC 39

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

PREPROCESO –
DATOS DE
ENTRADA.

NUMERO DE NUDOS. NUMERO DE

MATERIALES – MODULO ELEMENTOS .
DE ELASTICIDAD DE CADA UNO DE CADA UNO
DEFINIR: DEFINIR:

COORDENADAS EN
AREA DE LA SECCION
“x” Y “y”

CONDICIONES DE
FRONTERA DE CADA INCIDENCIA: NUDO
APOYO – INICIAL Y NUDO FINAL
RESTRINGIDO O LIBRE

FUERZAS ATUACTES
EN “x” o/y “y”

Figura 2 - 17. Diagrama de flujo del aplicativo computacional. Pre – proceso.

Los datos de entrada mencionados en el anterior diagrama de flujo, el usuario los

puede ingresar a través del formulario que se presenta en la Figura 2 - 18 y que
corresponde al ejercicio que se observa en la Figura 2 - 19.

Figura 2 - 18. Formulario para ingresar los datos de análisis matricial -Pre proceso.

Edgar Eduardo Muñoz, IC, MIC 40

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

2,5 m
500 kN
200 kN
4 Elemento 3

Elemento 2 1
1,5m
Elemento 4

Elemento 5 3m
2
Elemento 1 1,5m
y

3
x

5m
Figura 2 - 19. Ejemplo 4: Armadura plana de 4 nudos y 5 elementos.

También el usuario puede ingresar los datos directamente en la hoja electrónica,

como se observa a continuación:

Número de nudos 4
Número de elementos 5
Módulo de elasticidad 200000000 KN/m2
0,0000 0
0
Datos de los nudos Datos de los elementos
Número X(m) Y(m) Rx Ry Fx (kN) Fy (kN) Número Ni Nf A(m2)
1 5 3 0 0 200 -500 1 3 2 0,0064
2 2,5 1,5 0 0 0 0 2 2 1 0,0064
3 0 0 1 1 0 0 3 4 1 0,0064
4 0 3 1 1 4 4 2 0,0064
5 3 4 0,0064

Tabla 2 - 6. Datos del pre – proceso ingresados directamente en la hoja electrónica.

Edgar Eduardo Muñoz, IC, MIC 41

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

En esta parte el aplicativo muestra la geometría de la cercha:

GEOMETRIA DEL MODELO

4
Y(m)

0
0 1 2 3 4 5 6 7 8

X(m)

Figura 2 - 20. Geometría del modelo estructural.

Posteriormente el porgrama aplicativo hace un post – proceso, cuyo diagrama de flujo

se expone a continuación:

Edgar Eduardo Muñoz, IC, MIC 42

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

POST -
PROCESO

CONSTRUCCION DE LA MATRIZ DE RIGIDEZ

GLOBAL DE CADA ELEMENTO( Ver 2-25)

ENSAMBLAJE DE LA MATRIZ DE RIGIDEZ DE LA

ESTRUCTURA.

DETERMINACION DEL
VECTOR DE FUERZAS
ACTUANTES

CONSTRUCCION DEL SISTEMA DE SOLUCION POR EL

METODO DE DESPLAZAMIENTOS( Ver 2-8)

DETERMINACION DE LAS DETERMINACIONDE LAS

DEFORMACIONES DE LOS REACIONES DE LOS GRADOS
GRADOS DE LIBERTAD( VER 2-9) RESTRINGIDOS(Ver 2-10)

DETERMINACION DE LAS
FUERZAS INTERNAS DE CADA
ELEMENTO

Figura 2 - 21. Diagrama de flujo del aplicativo computacional. post – proceso.

.

Para el ejercicio de la Figura 2 - 19, los resultados basados en el diagrama de flujo

anterior se exponen en a continuación. La matriz de rigidez de los cinco elementos es:

Edgar Eduardo Muñoz, IC, MIC 43

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Matriz de rigidez del elemento 3-2

U3x U3y U2x U2y
U3x 322820,9 193692,5 -322820,9 -193692,5
U3y 193692,5 116215,5 -193692,5 -116215,5
U2x -322820,9 -193692,5 322820,9 193692,5
U2y -193692,5 -116215,5 193692,5 116215,5
Matriz de rigidez del elemento 2-1
U2x U2y U1x U1y
U2x 322820,9 193692,5 -322820,9 -193692,5
U2y 193692,5 116215,5 -193692,5 -116215,5
U1x -322820,9 -193692,5 322820,9 193692,5
U1y -193692,5 -116215,5 193692,5 116215,5
Matriz de rigidez del elemento 4-1
U4x U4y U1x U1y
U4x 256000,0 0,0 -256000,0 0,0
U4y 0,0 0,0 0,0 0,0
U1x -256000,0 0,0 256000,0 0,0
U1y 0,0 0,0 0,0 0,0
Matriz de rigidez del elemento 4-2
U4x U4y U2x U2y
U4x 322820,9 -193692,5 -322820,9 193692,5
U4y -193692,5 116215,5 193692,5 -116215,5
U2x -322820,9 193692,5 322820,9 -193692,5
U2y 193692,5 -116215,5 -193692,5 116215,5
Matriz de rigidez del elemento 3-4
U3x U3y U4x U4y
U3x 0,0 0,0 0,0 0,0
U3y 0,0 426666,7 0,0 -426666,7
U4x 0,0 0,0 0,0 0,0
U4y 0,0 -426666,7 0,0 426666,7

Figura 2 - 22. Matriz de rigidez de cada elemento.

Basados en la rigidez de cada elemento que hace par de la estructura, se presenta la

matriz de rigidez de la estructura, basado en un proceso de ensamblaje:

U1x U1y U2x U2y U3x U3y U4x U4y

U1x 578820,9 193692,5 -322820,9 -193692,5 0,0 0,0 -256000,0 0,0
U1y 193692,5 116215,5 -193692,5 -116215,5 0,0 0,0 0,0 0,0
U2x -322820,9 -193692,5 968462,6 193692,5 -322820,9 -193692,5 -322820,9 193692,5
U2y -193692,5 -116215,5 193692,5 348646,5 -193692,5 -116215,5 193692,5 -116215,5
U3x 0,0 0,0 -322820,9 -193692,5 322820,9 193692,5 0,0 0,0
U3y 0,0 0,0 -193692,5 -116215,5 193692,5 542882,2 0,0 -426666,7
U4x -256000,0 0,0 -322820,9 193692,5 0,0 0,0 578820,9 -193692,5
U4y 0,0 0,0 193692,5 -116215,5 0,0 -426666,7 -193692,5 542882,2
Figura 2 - 23. Matriz de rigidez de la estructura

Basados en las condiciones de frontera se tiene el planteamiento del sistema de

solución, que incluye el vector de fuerzas actuantes conocidas y las correspondientes
sub matrices y sub vectores:

Edgar Eduardo Muñoz, IC, MIC 44

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

U1x U1y U2x U2y

F1x=200 578820,87 193692,52 -322820,87 -193692,52 U1x=?
F1y=-500 193692,52 116215,51 -193692,52 -116215,51 U1y=?
F2x=0 -322820,87 -193692,52 968462,60 193692,52 U2x=?
F2y=0 -193692,52 -116215,51 193692,52 348646,54 U2y=?
F3x=? 0,00 0,00 -322820,87 -193692,52 U3x=0
F3y=? 0,00 0,00 -193692,52 -116215,51 U4y=0
F4x=? -256000,00 0,00 -322820,87 193692,52 U4x=0
F4y=? 0,00 0,00 193692,52 -116215,51 U5y=0

Figura 2 - 24. Sistema de solución

Basados en lo anterior, se determinan las deformaciones en de los grados de libertad

en cada nudo, siendo el mayor la del nudo 1 en dirección Y, que tiene un valor de 1.5
cm:

U1x 0,004036458 1
U1y -0,015332134 2
U2x -0,001290705 3
U2y -0,002151176 4
Figura 2 - 25. Deformaciones en los grados de libertad en los nudos

El aplicativo muestra la elástica de dicha cercha, como se observa a continuación:

Edgar Eduardo Muñoz, IC, MIC 45

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

GEOMETRIA DEL MODELO

4
Y(m)

0
0 1 2 3 4 5 6 7 8

X(m)

Figura 2 - 26. Elástica de la armadura

Posteriormente se determinan las reacciones de los grados restringidos en cada nudos,

como se muestra a continuación y en la Figura 2 - 29.

F3x=? 833.333 kN
F3y=? 500.000 kN
F4x=? -1.033.333 kN
F4y=? 0,000 kN
Figura 2 - 27. Reacciones en los grados restringidos de los nudos

Se revisa equilibrio en cada dirección:

 Fy  0
 500  500  0 _ OK
 Fx  0
200  833.33  1033.33  0 _ OK

Edgar Eduardo Muñoz, IC, MIC 46

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Se determinan las fuerzas internas:

F3-2 -971.825
F2-1 -971.825
F4-1 1.033.333
F4-2 0,000
F3-4 0,000

Figura 2 - 28. Fuerzas internas

500 kN

833.33 kN 1033,33 kN (tensión)

4

1
0 kN

0 kN 971,825 kN (Compresión)
971,825 kN (Compresión)
2
0 kN

1033.33 kN

3
Figura 2 - 29. Reacciones en los apoyos y fuerzas internas de cada elemento

[Link].2 Métodos de energía

Se han desarrollado en diferentes épocas métodos basados en el principio de la

conservación de la energía (Ver Nelson et al [5], Uribe [6] y Riveros [7], Hibbeler
[8], entre otros), que ha sido una contribuciones dentro del análisis de estructuras,
que han permitido la determinación de deflexiones en puntos de armaduras, vigas y
pórticos . El principio se basa en los siguientes conceptos:

Edgar Eduardo Muñoz, IC, MIC 47

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 El trabajo externo hecho por la fuerzas o cargas externas que actúan sobre una
estructura se trasforman en un trabajo interno de energía de deformación.
 Si la acción de dichas fuerzas, no excede el límite elástico, la estructura se
debe recuperar en el momento que las cargas se retiran.
 Por la conservación de la energía, se cumple con el principio básico que el
trabajo externo (We) es igual al trabajo interno (Wi).

WE  WI (2 - 27)

 El trabajo externo (dWe) para estructuras producto de una fuerza F, que

produce una desplazamiento “dx”, se puede expresar matemáticamente como:

(2 - 28)
WE   Fdx
x

 Dicho concepto, se puede comprender físicamente, si se considera por

ejemplo una barra simple (lo cual hacen diferentes autores), que esta
empotrada en un extremo y que tiene una carga axial que se aplica desde 0
hasta “P”, obteniéndose una deformación “  “, como se muestra en la
siguiente figura.

F
A .L y E
P
1
We  P
2

  X

Edgar Eduardo Muñoz, IC, MIC 48

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 30. Trabaja externo que se produce en una barra producto de una carga
axial. Fuente: Adaptada de las referencias [5] y [8]

 Resolviendo la integral definida presentada en la ecuación (2 - 28) desde 0

hasta en “  “, es decir aplicando la fuerza gradualmente a la barra desde 0
hasta “P”, se obtiene que el trabajo externo que se muestra a continuación( lo
que representa el área baja la curva de la figura anterior):

1 (2 - 29)
WE  P
2

 Los anterior se puede aplicar en forma general cuando no es una fuerza si no

un momento, teniendo en cuenta que cuando se aplica un momento
gradualmente desde o hasta “M” se obtiene un Angulo y se puede expresar de
las siguiente forma:

Edgar Eduardo Muñoz, IC, MIC 49

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 (2 - 30)
WE   Md
0

1 (2 - 31)
WE  M
2

 El trabajo externo que se produce en la barra al aplicar la carga se transforma

en una energía de deformación (trabaja interno). Para ello, se aplica la ley de
Hooke, para lo cual se tiene:

P (2 - 32)
 
A

 (2 - 33)
  
E L
PL (2 - 34)

AE
 Remplazando la ecuación (2 - 34) en la ecuación (2 - 29) y aplicando el
principio de la conservación de la energía, se tiene que el trabajo interno(
energía de deformación):

1 P2L (2 - 35)
WI 
2 AE
 Para deducir la energía de deformación que produce un momento, se debe
recordar que por mecánica de materiales, teniendo en cuenta que se trabaja en
el rango elástico – lineal, la siguiente relación:

M (2 - 36)
d  dx
EI

 Remplazando la ecuación (2 - 36) en la ecuación (2 - 30) y aplicando el

principio de la conservación de la energía ( Ver ecuación (2 - 27)), se tiene:

L
M 2 dx (2 - 37)
WI  
0
EI

Edgar Eduardo Muñoz, IC, MIC 50

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

[Link].2.1 Trabajo Virtual

Es un método clásico y versátil para evaluar deflexiones elásticas producidas por

cargas verticales, errores de construcción o cambios de temperatura. El método se
basa en el principio de la conservación de la energía, por lo tanto el trabajo virtual
externo hecho por el sistema de carga es igual al trabajo virtual interno efectuado por
las fuerzas internas causadas por él. Este método fue desarrollado por John Bernoulli
en 1917 y se llama según diferentes autores también el método de la carga unitaria.
Permite determinar la deformación en un punto de una armadura, viga o pórtico.

El trabajo virtual externo se define como:

WvirtualE  1x (2 - 38)

El trabajo virtual interno de toda la estructura será la suma de los trabajos efectuados
en las barras:

WvirtualI   udL (2 - 39)

En el caso de armadura, “dL” se puede representar como el valor equivalente a “ 

SL
L” y es equivalente a: . Al igualar los dos trabajos virtuales se tiene:
AE

Edgar Eduardo Muñoz, IC, MIC 51

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Aplicando el principio de la conservación de la energía se tiene:

WVirtualI  WVirtualE (2 - 40)

 SL 
  U i 
 AE  i

Dónde:

S: Fuerzas internas en cada elemento producidas por las cargas externas

U: Fuerza interna en los elementos producidos por la carga ficticia u unitaria
L: Longitud de cada elemento
A: Área de cada elemento
E. Módulo de elasticidad de cada elemento

Se aplica a continuación este método para el problema anterior que se presentó en la

Figura 2 - 19 y la tarea consiste en determinar la deformación vertical en el nudo 1.
Para ello, se coloca una fuerza ficticia u untaría en dicho nudo de un valor de 1 kN,
como se observa en la siguiente figura:

1 kN
2

Figura 2 - 31. Fuerza ficticia en el nudo 1 y en dirección “y” aplicada en el modelo

estructural

Edgar Eduardo Muñoz, IC, MIC 52

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Se determinan las fuerzas internas (u) en cada elemento producto de esta fuerza
ficticia y se aplica la ecuación (2 - 40), base del principio de trabajo virtual:

Elemento L(m) A (m2) u(kN) S(kN) (u(SL)/(AE))

1 2.9 0.00640 -1.900 -971.8 0.00430
2 2.9 0.00640 -1.900 -971.8 0.00430
3 5.0 0.00640 1.700 1033.3 0.00673
4 2.9 0.00640 0.000 0.0 0.00000
5 3.0 0.00640 0.000 0.0 0.00000

0.01533 m

Tabla 2 - 7. Resultados utilizando el método de trabajo virtual

Se obtuvo una deformación vertical de 0.01533 m, la cual coincide con los resultados
obtenidos en el ejercicio del numeral [Link].1.6 . Este método también se agregó en
el aplicativo computacional mencionado en el numeral [Link].1.6 , de tal forma
puede determinar por el método de trabajo virtual la deformación en cualquiera de los
nudos de la armadura. Por ejemplo, si se quisiera por este método determinar la
deformación del nudo 1 pero en dirección “x”, se tiene:

Figura 2 - 32. Formulario para ingresar los datos del método de trabajo virtual -Pre
proceso.

Para este caso, el programa establece es una fuerza ficticia en el nudo 2 en dirección
x, como se muestra en la figura:

Edgar Eduardo Muñoz, IC, MIC 53

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

4
1 kN
1

3
Figura 2 - 33. Fuerza ficticia en el nudo 1 y en dirección “x” aplicada en el modelo
estructural

Se determinan las fuerzas internas (u) en cada elemento producto de esta fuerza
ficticia y se aplica la ecuación antes mencionada base del principio de trabajo virtual:

Elemento L(m) A (m2) u(kN) S(kN) (u(SL)/(AE))

1 2,9 0,0 0,0 -971,8 0,00000
2 2,9 0,0 0,0 -971,8 0,00000
3 5,0 0,0 -1,0 1033,3 -0,00410
4 2,9 0,0 0,0 0,0 0,00000
5 3,0 0,0 0,0 0,0 0,00000

-0,00410 m

Tabla 2 - 8. Resultados utilizando el método de trabajo virtual

Por otro lado, se desea determinar por el método de trabajo virtual la deformación
horizontal del nudo 1 de la estructura que se muestra en la Figura 2 - 12, producto
de la carga P1, considerando los siguientes datos:

 E=200000000 kN/m2
 L=1.5 m

Edgar Eduardo Muñoz, IC, MIC 54

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 P1=500 kN
 A=0.0007 m2

De lo anterior se tiene los siguientes resultados, con lo cual se determinó que la

deformación horizontal es de 0.00379 m:

Elemento L(m) A (m2) u(kN) S(kN) (u(SL)/(AE))

1 1.1 0.00070 0.000 353.6 0.00189
2 1.1 0.00070 -0.700 -353.6 0.00189
0.00379

Tabla 2 - 9. Resultados utilizando el método de trabajo virtual

[Link].2.2 Método de Castigliano

Otro de los métodos clásicos de análisis estructural fue el desarrollado por Carlo
Alberto Castigliano (9 de noviembre de1847, Asti - 25 de octubre de 1884, Milán)
que se basa en el concepto de energía interna de deformación. Catigliano fue un
matemático y físico italiano que propuso dos teoremas para el análisis estructural de
estructuras. Uno de los teoremas expresa lo siguiente:

“La componente de desplazamiento del punto de aplicación de una acción sobre

una estructura, en la dirección de dicha acción, se puede obtener evaluando la
primera derivada parcial de la energía interna de deformación de la estructura
con respecto a la acción aplicada”

La teoría es válida para determinar por ejemplo deformaciones o rotaciones en un

punto definido para estructuras isostáticas o estáticamente determinadas. Y también
para determinar fuerzas internas de estructuras estáticamente indeterminadas.

La expresión general del método se puede expresar como la derivada parcial de la

energía interna de deformación en función de la fuerza aplicada:

W (2 - 41)

P

Edgar Eduardo Muñoz, IC, MIC 55

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Para su aplicación, se presenta a continuación un ejemplo, que consiste en la

determinación de la deformación vertical que tiene el punto 1 de la armadura de la
figura, producto de una fuerza vertical de 500 KN. Para ello, se coloca una carga
adicional imaginaria “P” en el punto 1:

(500+P) kN

Elemento 1 Elemento 2

3,0m
y
3
2
x

3m 3m
Figura 2 - 34. Ejemplo 5: Armadura plana de 3 nudos y 2 elementos. Con cuatro
restricciones

Como es una armadura isostática, es decir estáticamente determinada, se procede a

determinar las correspondientes reacciones, utilizando los principios de equilibrio que
se aplican en estática. Se revisa sumatoria de fuerzas en dirección “y” e igualando a 0,
para lo cual se tiene:

 Fy  0  R 2y  R3 y  (500  P)

Después sumatoria de fuerzas en dirección “x” e igualando a 0:

Edgar Eduardo Muñoz, IC, MIC 56

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 Fx  0   R 2x  R3 x  0

Se hace sumatoria de momentos con respecto al nudo 4 e igualando a 0:

M punto 2  0  (500  P) * 3  R3 y * 6
R2 y  250  0.5 P
R3 y  250  0.5 P

Finalmente las reacciones son las siguientes:

(500+P) kN

Elemento 1 Elemento 2

3,0m
y
3
2
x

3m 3m

(250+0,5P) kN
(250+0,5P) kN
Figura 2 - 35. Reacciones de la armadura del ejemplo 5

Posteriormente se determinan las fuerzas internas en los elementos en función de las

reacciones previamente determinadas. Utilizando el método de los nudos para
armaduras, se tiene para el nudo 2 lo siguiente:

Edgar Eduardo Muñoz, IC, MIC 57

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 Fy nudo 2  0  250  0.5P  F21 * Sin(45)

F21  353.55  0.707P

Determinada la fuerza interna de la barra 2-1, se hace lo mismo para el nudo 3,

obteniéndose la fuerza interna de la barra 3-1:

F31  353.55  0.707P

Aplicamos el teorema de Castigliano, se tiene que para armaduras la deformación

para el punto de estudio, está dada por:

S L (2 - 42)
   S( )
P AE

Para resolver esta ecuación, se organizan los resultados en el siguiente cuadro:

E 200000000 kN/m2
Barra L(m) A(m2) L/AE Fuerza interna S
( )
S L
 S ( P ) AE
P
(kN)

´2-1 4.24 0.00064 0.000033125 ´353.55+0.707P 0.707 0.00827992
´3-1 4.24 0.00064 0.000033125 ´353.55+0.707P 0.707 0.00827992

0.01655984

Tabla 2 - 10. Resultados utilizando el método de Castigliano. Nota: En el momento

de hacer la multiplicación “P” se considera 0.

Esto representa una deformación vertical en el nudo 1 de:

Edgar Eduardo Muñoz, IC, MIC 58

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

S L
   S( )  0.00166m
S AE

Esta deformación coincide al revisarla utilizando el método de análisis matricial, lo

cual se logra al resolver el siguiente sistema:

U1x U1y
F1x=0 30169,89 0,00 U1x=?
F1y=-500 0,00 30169,89 U1y=?
F2x=? = -15084,94 -15084,94 U2x=0
F2y=? -15084,94 -15084,94 U3y=0
F3x=? -15084,94 15084,94 U3x=0
F3y=? 15084,94 -15084,94 U4y=0

Figura 2 - 36. Sistema de solución matricial armadura ejemplo 5

Aplicando la ecuación (2 - 9), se tiene que la deformación en el nudo 1 es de 0.00165

m.

2.1.2 Armaduras espaciales

[Link] Definición y antecedentes

Las armaduras o cerchas espaciales son un sistema estructural en el espacio

constituido por elementos conformados principalmente por triángulos y sometidos
generalmente a fuerzas axiales de tensión o compresión, cuando deben soportar
cargas externas aplicadas en específicamente en sus nudos y no sobre los elementos.
Cuando tienen que asumir cargas externas aplicadas sobre sus elementos, además de

Edgar Eduardo Muñoz, IC, MIC 59

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

las cargas en los nudos, sus elementos no solamente tienen fuerzas internas axiales si
no también combinadas que pueden ser de flexo compresión o flexo tensión. Son
estructuras importantes, que han sido útiles para la sociedad, empleadas para el
diseño y la construcción de diferentes partes o componentes de algunas obras civiles,
tales como: torres de comunicaciones y cubiertas (Ver Figura 2 - 37).

(b)

(a)
Figura 2 - 37. Ejemplos de estructuras de armaduras espaciales

[Link] Método matricial

[Link].1 Fundamentos

La solución por el método matricial de este tipo de estructuras se basa en lo mismo

que para las armaduras planas, en la ecuación (2 - 8) deducida en este documento.
Lo que cambia, es que cada nudo tiene tres grados de libertad o restringidos, es decir
posibilidades de desplazamientos en el sentido “X”, “Y” o “Z” Por lo anterior, es
necesario deducir o determinar la matriz de rigidez global de dicho elemento
arbitrariamente orientado. Para ello, se debe determinar la matriz de rigidez local del
elemento, es decir cuando el elemento este solamente sobre el eje “X” de un sistema
tridimensional de coordenadas, como se observa a continuación:

Edgar Eduardo Muñoz, IC, MIC 60

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

E,A y L

j
i X

Figura 2 - 38. Elemento sobre el eje X. Fuente: Propia.

Aplicando la ecuación (2 - 18) de mecánica de sólidos, explicada anteriormente, se

determina la matriz de rigidez local del elemento de la Figura 2 - 38, teniendo en
cuenta que tiene tres grados de libertad o restringidos por nudo.

 AE AE 
 fix   0 0  0 0 uix 
   L L  
 fiy   0 0 0 0 0 0 uiy 

 fiz   0 0 0 0 0 0 uiz 
    AE AE   (2 - 43)
 fjx   0 0 0 0 ujx 
 fjy   L L ujy 
   0 0 0 0 0 0  
 fjz   ujz 
 0 0 0 0 0 0  

Edgar Eduardo Muñoz, IC, MIC 61

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Para determinar la matriz de rigidez global se aplica la ecuación (2 - 26), donde es

necesario establecer la matriz de transformación. Para este caso dicha matriz está en
función de los cosenos directores y se expresa de la siguiente forma:

 fix  Cx Cy Cz 0 0 0   Fix 

    
 fiy      0 0 0  Fiy 

 fiz 
    0 0 0  Fiz 

    (2 - 44)
 fjx   0 0 0 Cx Cy Cz   Fjx 
 fjy   0 0 0     Fjy 
    

 fjz 
  0 0 0     
 Fjz 


Donde los cosenos directores son:

(2 - 45)
xj  xi
Cx 
L

(2 - 46)
yj  yi
Cy 
L

(2 - 47)
zj  zi
Cz 
L

(2 - 48)
L  ( xj  xi)2  ( yj  yi)2  ( zj  zi)2

Edgar Eduardo Muñoz, IC, MIC 62

 Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Al realizar el triple producto expresado en la ecuación (2 - 26), se tiene que la matriz

de rigidez global del elemento de armadura espacial es:

 Fix   Cx 2  Uix 
 Fiy   2  
   CxCy Cy  Uiy 
 Fiz  AE  CxCz CyCz Cz 2  Uiz 
    
 Fjx  L   Cx  CxCy  CxCz Cx
2 2
 Ujx 
 Fjy   CxCy  Cy 2
 CyCz CxCy Cy 2  Ujy 
    
 Fjz    CxCz  CyCz  Cz 2 CxCz CyCz Cz 2  Ujz 

(2 - 49)

Por otro lado, para solucionar las fuerzas internas se debe aplicar la siguiente
ecuación:

Ujx  Uix  (2 - 50)

Sij 
AE
Cx Cy Cz Ujy  Uiy 
L
Ujz  Uiz 

[Link].2 Ejemplo: Solución numérica de una torre en el espacio

El ejercicio consiste en la determinación de las deformaciones, reacciones y fuerzas

internas de la siguiente armadura espacial (torre).

Edgar Eduardo Muñoz, IC, MIC 63

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

1000 kN
1 Nudos X(m) Y(m) Z(m)
1 0 3 0
2 3 0 0
3 0 0 -3
4 -3 0 0
3m 5 0 0 3

3
3m 4 3m
X
2
5 3m

3m
Z

Figura 2 - 39. Torre tridimensional y coordenadas en el espacio de cada punto.

Fuente: Propia.

Se tuvo en para la solución del sistema un módulo de elasticidad de 200000 Mpa y las
siguientes condiciones de frontera, a nivel de tipos de apoyos y fuerzas en los nudos

Edgar Eduardo Muñoz, IC, MIC 64

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Numero de nudos 5
Numero de elementos 4
Módulo de elasticidad 200000000

Coordenadas Restricciones Fuerzas en los nudos Incidencias de elementos y propiedades

Element
X(m) Y(m) Z(m) Ry Rz Rz Fx Fy Fz o Ni Nf A(m2)
0 3 0 0 0 0 0 -1000 0 1 2 1 0,0064
3 0 0 1 1 1 0 0 0 2 3 1 0,0064
0 0 -3 1 1 1 0 0 0 3 4 1 0,0064
-3 0 0 1 1 1 0 0 0 4 5 1 0,0064
0 0 3 1 1 1 0 0 0

Tabla 2 - 11. Datos de entrada.

Se determina la matriz de rigidez de cada elemento:

U2x U2y U2z U1x U1y U1z

U2x 150849,4 -150849,4 0,0 -150849,4 150849,4 0,0
U2y -150849,4 150849,4 0,0 150849,4 -150849,4 0,0
U1Z 0,0 0,0 0,0 0,0 0,0 0,0
U1x -150849,4 150849,4 0,0 150849,4 -150849,4 0,0
U1y 150849,4 -150849,4 0,0 -150849,4 150849,4 0,0
U1z 0,0 0,0 0,0 0,0 0,0 0,0

U3x U3y U3z U1x U1y U1z

U3x 0,0 0,0 0,0 0,0 0,0 0,0
U3y 0,0 150849,4 150849,4 0,0 -150849,4 -150849,4
U1Z 0,0 150849,4 150849,4 0,0 -150849,4 -150849,4
U1x 0,0 0,0 0,0 0,0 0,0 0,0
U1y 0,0 -150849,4 -150849,4 0,0 150849,4 150849,4
U1z 0,0 -150849,4 -150849,4 0,0 150849,4 150849,4

U4x U4y U4z U1x U1y U1z

U4x 150849,4 150849,4 0,0 -150849,4 -150849,4 0,0
U4y 150849,4 150849,4 0,0 -150849,4 -150849,4 0,0
U1Z 0,0 0,0 0,0 0,0 0,0 0,0
U1x -150849,4 -150849,4 0,0 150849,4 150849,4 0,0
U1y -150849,4 -150849,4 0,0 150849,4 150849,4 0,0
U1z 0,0 0,0 0,0 0,0 0,0 0,0

U5x U5y U5z U1x U1y U1z

U5x 0,0 0,0 0,0 0,0 0,0 0,0
U5y 0,0 150849,4 -150849,4 0,0 -150849,4 150849,4
U1Z 0,0 -150849,4 150849,4 0,0 150849,4 -150849,4
U1x 0,0 0,0 0,0 0,0 0,0 0,0
U1y 0,0 -150849,4 150849,4 0,0 150849,4 -150849,4
U1z 0,0 150849,4 -150849,4 0,0 -150849,4 150849,4
Figura 2 - 40. Matriz de rigidez de cada elemento.

Edgar Eduardo Muñoz, IC, MIC 65

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Posteriormente se determina la matriz de rigidez de la estructura.

301698,89 0,00 0,00 -150849,45 150849,45 0,00 0,00 0,00 0,00 -150849,45 -150849,45 0,00 0,00 0,00 0,00
0,00 603397,79 0,00 150849,45 -150849,45 0,00 0,00 -150849,45 -150849,45 -150849,45 -150849,45 0,00 0,00 -150849,45 150849,45
0,00 0,00 301698,89 0,00 0,00 0,00 0,00 -150849,45 -150849,45 0,00 0,00 0,00 0,00 150849,45 -150849,45
-150849,45 150849,45 0,00 150849,45 -150849,45 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00
150849,45 -150849,45 0,00 -150849,45 150849,45 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00
0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00
0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00
0,00 -150849,45 -150849,45 0,00 0,00 0,00 0,00 150849,45 150849,45 0,00 0,00 0,00 0,00 0,00 0,00
0,00 -150849,45 -150849,45 0,00 0,00 0,00 0,00 150849,45 150849,45 0,00 0,00 0,00 0,00 0,00 0,00
-150849,45 -150849,45 0,00 0,00 0,00 0,00 0,00 0,00 0,00 150849,45 150849,45 0,00 0,00 0,00 0,00
-150849,45 -150849,45 0,00 0,00 0,00 0,00 0,00 0,00 0,00 150849,45 150849,45 0,00 0,00 0,00 0,00
0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00
0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00
0,00 -150849,45 150849,45 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 150849,45 -150849,45
0,00 150849,45 -150849,45 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 0,00 -150849,45 150849,45

Figura 2 - 41. Matriz de rigidez de la estructura.

Se plante el sistema de solución:

F1x=0 301698,8933 0 0 U1x=?

F1y=-1000 0 603397,7866 0 U1y=?
F1z=0 0 0 301698,8933 U1z=?
F2x=? -150849,4467 150849,4467 0 0
F2y=? 150849,4467 -150849,4467 0 0
F2z=? 0 0 0 0
F3x=? 0 0 0 0
F3y=? = 0 -150849,4467 -150849,4467 0
F3z=? 0 -150849,4467 -150849,4467 0
F4x=? -150849,4467 -150849,4467 0 0
F4y=? -150849,4467 -150849,4467 0 0
F4z=? 0 0 0 0
F5x=? 0 0 0 0
F5y=? 0 -150849,4467 150849,4467 0
F5z=? 0 150849,4467 -150849,4467 0

Figura 2 - 42. Sistema de solución por el método de los desplazamientos.

Aplicando la ecuación (2 - 9), se obtienen las siguientes deformaciones:

Ux1=0,0000000
Uy1=-0,0016573
Uz1=0,0000000
Figura 2 - 43. Deformaciones determinadas en el nudo 1.

Aplicando la ecuación (2 - 10), se obtienen las reacciones:

Edgar Eduardo Muñoz, IC, MIC 66

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

FX2= -250.000 ´kN

FY2= 250.000 ´kN
FZ2= 0,000 ´kN
FX3= 0,000 ´kN
FY3= 250.000 ´kN
FZ3= 250.000 ´kN
FX4= 250.000 ´kN
FY4= 250.000 ´kN
FZ4= 0,000 ´kN
FX5= 0,000 ´kN
FY5= 250.000 ´kN
FZ5= -250.000 ´kN
Figura 2 - 44. Reacciones de la armadura espacial
Se revisa equilibrio en cada dirección:
 Fy  0
 1000  250 * 4  0 _ OK
 Fx  0
 250  250  0 _ OK
 Fz  0
 250  250  0 _ OK
Se determina las fuerzas internas en cada elemento, aplicando la ecuación (2 - 50):

F2-1 -353.553 kN
F3-1 -353.553 kN
F4-1 -353.553 kN
F5-1 -353.553 kN
Figura 2 - 45. Fuerzas internas de los elementos de la armadura espacial

Son fuerzas axiales a compresión.

[Link].3 Aplicativo computacional para el análisis de armaduras espaciales

Para resolver armaduras especiales se desarrolló también un programa en Visual

Basic para Excel que permite resolver cualquier estructura de este tipo para una
condición de carga dada. Se presenta a continuación su diagrama de flujo, que está

Edgar Eduardo Muñoz, IC, MIC 67

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

constituido por dos partes: pre procesos y post proceso. Con respecto a pre proceso
se tiene:

PREPROCESO –
DATOS DE
ENTRADA.

NUMERO DE
NUMERO DE NODOS. ELEMENTOS .
MATERIALES – MODULO DE DE CADA UNO
ELASTICIDAD DEFINIR: DE CADA UNO
DEFINIR:

COORDENADAS EN
AREA DE LA SECCION
“x”, “y”, Y “z”

CONDICIONES DE
FRONTERA DE CADA INCIDENCIA: NUDO INICIAL
APOYO – RESTRINGIDO Y NUDO FINAL
O LIBRE

FUERZAS ATUACTES
EN “x” , “y” y “z”

Figura 2 - 46. Diagrama de flujo del aplicativo computacional. Pre – proceso.

Los datos de entrada mencionados en el anterior diagrama de flujo, el usuario los

puede ingresar a través del formulario que se presenta en la Figura 2 - 18 y que
corresponde al ejercicio que se observa en la Figura 2 - 19.

Edgar Eduardo Muñoz, IC, MIC 68

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 47. Formulario para ingresar los datos de análisis matricial -Pre proceso.
Y

1000 kN
500 kN 1 Nudos X(m) Y(m) Z(m)
1 0 3 0
2 3 0 0
3 0 0 -3
4 -3 0 0
3m 5 0 0 3

3
3m 4 3m
X
2
5 3m

3m
Z
Figura 2 - 48. Ejemplo 5: Armadura espacial de 5 nudos y 4 elementos.

También el usuario puede ingresar los datos directamente en la hoja electrónica,

como se observa a continuación:

Edgar Eduardo Muñoz, IC, MIC 69

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Numero de nudos 5
Numero de elementos 4
Módulo de elasticidad 200000000

Coordenadas Restricciones Fuerzas en los nudos Incidencias de elementos y propiedades

Element
X(m) Y(m) Z(m) Ry Rz Rz Fx Fy Fz o Ni Nf A(m2)
0 3 0 0 0 0 500 -1000 0 1 2 1 0,00064
3 0 0 1 1 1 0 0 0 2 3 1 0,00064
0 0 -3 1 1 1 0 0 0 3 4 1 0,00064
-3 0 0 1 1 1 0 0 0 4 5 1 0,00064
0 0 3 1 1 1 0 0 0

Tabla 2 - 12. Datos del pre – proceso ingresados directamente en la hoja electrónica.

Posteriormente el porgrama aplicativo hace un post – proceso, cuyo diagrama de flujo

se expone en la Figura 2 - 21. Para el ejercicio de la Figura 2 - 48, los resultados
basados en el diagrama de flujo anterior se exponen en a continuación. La matriz de
rigidez de los cinco elementos se presenta en la Figura 2 - 40 y la matriz de rigidez
de la estructura en la Figura 2 - 41 . Basados en las condiciones de frontera se tiene
el planteamiento del sistema de solución, que incluye el vector de fuerzas actuantes
conocidas y las correspondientes sub matrices y sub vectores:

F1x=500 301698,89 0,00 0,00 U1x=?

F1y=-1000 0,00 603397,79 0,00 U1y=?
Fiz=0 0,00 0,00 301698,89 U1z=?
F2x=? -150849,45 150849,45 0,00 0
F2y=? 150849,45 -150849,45 0,00 0
F3z=? 0,00 0,00 0,00 0
F3x=? 0,00 0,00 0,00 0
F2y=? = 0,00 -150849,45 -150849,45 0
F3z=? 0,00 -150849,45 -150849,45 0
F4x=? -150849,45 -150849,45 0,00 0
F4y=? -150849,45 -150849,45 0,00 0
F4z=? 0,00 0,00 0,00 0
F5x=? 0,00 0,00 0,00 0
F5y=? 0,00 -150849,45 150849,45 0
F5z=? 0,00 150849,45 -150849,45 0

Figura 2 - 49. Sistema de solución

Basados en lo anterior, se determinan las deformaciones en de los grados de libertad

en cada nudo, siendo el mayor la del nudo 1 en dirección Y, que tiene un valor de 1.5
cm:

Edgar Eduardo Muñoz, IC, MIC 70

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Ux1=0,0165728
Uy1=-0,0165728
Uz1=0,0000000
Figura 2 - 50. Deformaciones en los grados de libertad en los nudos

Posteriormente se determinan las reacciones de los grados restringidos en cada nudos,

como se muestra a continuación:

F2x= -500.000 ´kN

F2y= 500.000 ´kN
F3z= 0,000 ´kN
F3x= 0,000 ´kN
F2y= 250.000 ´kN
F3z= 250.000 ´kN
F4x= 0,000 ´kN
F4y= 0,000 ´kN
F4z= 0,000 ´kN
F5x= 0,000 ´kN
F5y= 250.000 ´kN
F5z= -250.000 ´kN
Figura 2 - 51. Reacciones en los grados restringidos de los nudos

Se revisa equilibrio en cada dirección:

 Fy  0
 1000  500  250  250  0 _ OK
 Fx  0
500  500  0 _ OK
 Fz  0
250  250  0 _ OK
Se determinan las fuerzas internas:

F2-1 -707.107
F3-1 -353.553
F4-1 0,000
F5-1 -353.553
Figura 2 - 52. Fuerzas internas

Edgar Eduardo Muñoz, IC, MIC 71

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

2.1.3 Vigas

[Link] Definición y antecedentes

La vigas es un sistema estructural de elementos horizontales sometidos generalmente

a cortante y momento que deben soportar cargas externas aplicadas en sus nudos y/o
sobre los elementos.

[Link] Métodos de análisis

[Link].1 Pendiente deflexión o giro - deflexión

Es un método que fue desarrollado originalmente por Heinrich Manderla y Otto Mohr
con el fin de estudiar los esfuerzos secundarios de armaduras. Más adelante G.A.
Meney desarrolló una versión complementada y la aplicó para el análisis estructuras
de vigas y pórticos interinados. Se llama así, porque relacionada las pendientes y
deflaciones desconocidas con la carga aplicada sobre una estructura. Los
fundamentos de este método se basan en las siguientes ecuaciones:

2 EI  (2 - 51)
M AB  (2 A   B  3 )  M AB
f

L L

2 EI  (2 - 52)
M BA  (2 B   A  3 )  M BA
f

L L

Edgar Eduardo Muñoz, IC, MIC 72

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Dónde:

MAB Momento interno en el nudo A del tramo A-B

MBA Momento interno en el nudo B del tramo B-A
E Módulo de elasticidad
I Momento de inercia
L Longitud del tramo A-B
A Giro del nudo A del tramo A-B
B Giro del nudo B del tramo B-A
 Pendiente generada por un desplazamiento. Su signo es negativo
L cuando la pendiente es negativa y positivo cuando su pendiente es
positiva.
MfAB Momento de empotramiento en el nudo A del tramo A-B
MfBA Momento de empotramiento en el nudo B del tramo B-A

A continuación se presentan algunos ejemplos de aplicación:

[Link].1.1 Ejemplo: Solución simbólica de una viga continua

A continuación se hace el análisis estructural para el caso general de una viga continua de dos
luces iguales con una carga distribuida.

1 2 3
L L

Figura 2 - 53. Ejemplo de viga continua de dos luces iguales

Se presenta a continuación su solución mediante variables simbólicas. Se establecen

inicialmente las ecuaciones en los nudos de los tramos 1-2 y 2-3:

2 EI WL2 (2 - 53)
M 12  0  (21   2 ) 
L 12

Edgar Eduardo Muñoz, IC, MIC 73

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

2 EI WL2 (2 - 54)
M 21  (2 2  1 ) 
L 12

2 EI WL2 (2 - 55)
M 23  (2 2   3 ) 
L 12

2 EI WL2 (2 - 56)
M 32  0  (2 3   2 ) 
L 12

Por equilibrio en el nudo 2 se debe cumplir:

M 21  M 23  0

2 EI 2 EI 2 EI 8EI 2 EI (2 - 57)
(2 2  1 )  (2 2   3 )  1  2  3  0
L L L L L

Organizando en forma matricial las ecuaciones (2 - 53), (2 - 56) y (2 - 57), se tiene:

 WL2   4 E 2 EI 
 12   0 
  L L 1 
 WL2   2 EI 4 EI   
  0   2 
 12   L L  
 0   2 EI 8EI 2 EI   3 
   L L L 

Resolviendo el sistema empleado el programa MATLAB, se tiene:

L = sym('L');
EI = sym('EI');
W = sym('W');
c=[4*EI/L,2*EI/L,0;0,2*EI/L,4*EI/L;2*EI/L,8*EI/L,2*EI/L]
b=[-W*L*L/12;W*L*L/12;0]
a=inv(c)*b

Edgar Eduardo Muñoz, IC, MIC 74

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 WL3 
1   48EI 
   
 2    0 3 
   WL 
 3  
 48EI 

Resolviendo las ecuaciones (2 - 54) y (2 - 55), se tiene:

3WL2
M 21  
24
3WL2
M 23 
24

Se revisa equilibrio en el tramo 1-2:

W
M 21
1 2
L
R2
R1
Figura 2 - 54. Tramo 1-2

Haciendo sumatoria de momentos iguales a cero en el nudo 2, se tiene:

WL
 R1 L  M 21  0
2

Edgar Eduardo Muñoz, IC, MIC 75

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

3WL2 WL2
 R1 L   0
24 2

3WL2 WL2
 
24 2 R
1
L

3WL2 WL2
 
R1  24 2  9 WL
L 24

Por simetría de la viga y equilibrio tenemos:

9
R3  R1  WL
24
30
R2  2WL  R1  R2  WL
24

A continuación se presentan el diagrama decortante y momento de la viga continua con carga

distruibuida en toda la luz.

Edgar Eduardo Muñoz, IC, MIC 76

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 55. Diagrama de cortante y momento

[Link].1.2 Ejemplo: Análisis de una vigueta que hace parte de un entrepiso

Este ejercicio consiste en realizar el análisis estructural de la vigueta continua, que se

muestra en la Figura 2 - 56(seleccionado en rojo), con el objeto de:

 Realizar una evaluación de la carga muerta

 Determinar rotaciones de los apoyos
 Determinar las reacciones en los apoyos ( verificación de equilibrio)
 Realizar el diagrama de cortante e identificar el valor máximo

Edgar Eduardo Muñoz, IC, MIC 77

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 Realizar el diagrama de momento e identificar el valor máximo

Para su desarrollo, considerar que la vigueta tiene:

 Una resistencia a la compresión del concreto del 21 Mpa.

 Una altura de 0.4m y un ancho de 0.12m.
 Que soportar una carga muerta de 6.8 kN/m2.

Figura 2 - 56. Plano de entrepiso en una dirección. Fuente: Ing. Ramiro Vallecilla

Inicialmente se construye, basado en el plano anterior, el sistema estático de la

vigueta:

Edgar Eduardo Muñoz, IC, MIC 78

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

C D
Wcm=6,8*0,82=5,6KN/ m

1
2 3 4
4,5 m 1,2 m
2,60 m

Figura 2 - 57. Sistema estático de vigueta (Nota el voladizo es de 1.2 m)

La carga por metro lineal sobre la vigueta corresponde a la carga por metro cuadrado
por e la separación entre viguetas, lo cual equivale a una carga de 5.6 kN.m. Por otro
lado el módulo de elasticidad se determina en función de la resistencia a la
compresión y es de 21538105.7 kN/m2 .El momento de inercia es: 0.00064 m4, por lo
que la relación EI es : 13784.38
2
kN.m Ahora se establecen las condiciones estructurales de la vigueta, basados en su
sistema estático:

M 12  0 (2 - 58)

M 21  M 23  0 (2 - 59)
M 32  M 34  0 (2 - 60)

Posteriormente se plantean las ecuaciones basándose en las ecuaciones (2 - 51) y (2 -

52). Se comienza con el nudo 2:

2(13784.38) 5.6 * 2.6 2 (2 - 61)

M 12  0  (21   2 ) 
2.6 12
21206.751  10603.37 2  3.15

Edgar Eduardo Muñoz, IC, MIC 79

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

2(13784.38) 5.6 * 2.6 2 (2 - 62)

M 21  (2 2  1 ) 
2.6 12
M 21  21206.75 2  10603.37 1  3.15

2(13784.4) 5.6 * 4.5 2 (2 - 63)

M 23  (2 2   3 ) 
4.5 12
M 23  12252.79 2  6126.39 3  9.45

Aplicando la ecuación (2 - 59), se tiene:

12252,79 2  6126.39 3  9.45  21206.76 2  10603.38 1  3.15  0 (2 - 64)

10603.38 1  33459.54 2  6126.39 3  6.3

Se pasa al nudo 3:

2(13784.4) 5.6 * 4.5 2 (2 - 65)

M 32  (2 3   2 ) 
4.5 12
M 32  12252.79 3  6126.39 2  9.45

M 34  1.2 * 5.6 * 1.2 * 0.5  4.032kN .m (2 - 66)

Aplicando la ecuación (2 - 60), se tiene:

12252.79 3  6126.398 2  9.45  4.032  0 (2 - 67)

6126.39 2  12252,.79 3  5.418

Se organizan las ecuaciones (2 - 61), (2 - 64) y (2 - 67), se plantea el siguiente

sistema de solución lineal:

Edgar Eduardo Muñoz, IC, MIC 80

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 3.15 21206.76 10603.38 0  1 

    
  6.3   10603.38 33459.54 6126.39   2 
 5.418   6126.39 12252.79  3 
   0

Para resolver el sistema se saca la inversa de la matriz y se obtiene los valores de la rotación
de los tres nudos:

1   - 4.18E - 07 radianes 

   
 2   - 0.0002962radianes 
   0.0005903radianes 
 3  

Con esto valores se determina el momento M21 al remplazarlos en la ecuación (2 - 62):

M 21  21206.75 2  10603.371  3.15  -9.44kN.m

Con el valor anterior se determina la reacción del nudo 1 y parte del nudo 2:

Wcm=6,8*0,82=5,6KN/ m

1 M21=9.44 kN.m

2
2,60 m

R21
R12

Figura 2 - 58. Determinación de reacciones primer tramo

Haciendo sumatoria de momentos iguales a cero en el nudo 2, se tiene:

Edgar Eduardo Muñoz, IC, MIC 81

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

5.6 * 2.6 * 2.6

 R12 * 2.6  9.44  0
2
R12  3.65kN
 Fy  0  5.6 * 2.6  R 12
 R21
R21  10.91kN

Se determina las reacciones en el tramo 2, haciendo sumatoria de momentos iguales a

cero en el nudo 2, se tiene:

C D

Wcm=6,8*0,82=5,6KN/ m

M32=4,03kN.m
M23=9,44 kN.m
2 3
4,5 m

R23 R32

Figura 2 - 59. Determinación de reacciones segundo tramo

5.6 * 4.5 * 4.5

 R23 * 4.5  4.03  9.44  0
2
R23  13.80kN
R 2  R23  R21  24.71kN
 Fy  0  4.5 * 5.6  R 23  R31
R33  11.4kN
R3  R33  R34  11.4  5.6 * 1.2  18.12kN

Edgar Eduardo Muñoz, IC, MIC 82

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Al hacer equilibrio de las fuerzas verticales iguales a cero de toda la viga, se tiene que
R3=18.12 kN. Conociendo las reacciones de la viga se construyen los diagramas de
cortante y momento que se muestran a continuación:

Diagrama de cortante

20

15

10
Cortante(KN)

0
0 1 2 3 4 5 6 7 8

-5

-10

-15
X(m)

Figura 2 - 60. Diagrama de cortante

Diagrama de momento

12

10

4
Momento(KN.m)

0
0 2 4 6 8

-2

-4

-6

-8

-10
x(m)

Figura 2 - 61. Diagrama de momento

Edgar Eduardo Muñoz, IC, MIC 83

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

[Link].1.3 Ejemplo: Análisis de una viga de acero empotrada en los extremos

Este ejercicio consiste en realizar el análisis estructural de una viga de acero,

localizada en una edificación industrial, como se muestra en la Figura 2 - 62.

W=20 kN/m

A B C D

5.0 m
5.0 m 3.0 m

Figura 2 - 62. Viga en acero con empotramiento en los extremos. Fuente : Tomada
de referencia [8].

Se requiere:

 Determinar las reacciones en los apoyos ( verificación de equilibrio)

 Realizar el diagrama de cortante
 Realizar el diagrama de momento

Por otro lado, considerar para este ejercicio que EI es constante y tiene un valor de
1.0, con lo cual los valores de rotación determinados no corresponden a la realidad.
Ahora, se establecen las condiciones estructurales de la vigueta, basados en su
sistema estático:

A  0 (2 - 68)
D  0 (2 - 69)

M BA  M BC  0 (2 - 70)
M CB  M CD  0 (2 - 71)

Edgar Eduardo Muñoz, IC, MIC 84

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Posteriormente se plantean las ecuaciones basándose en (2 - 51) y (2 - 52). Se

comienza con el tramo A-B y teniendo en cuenta la ecuación (2 - 68):

2(1) (2 - 72)
M AB  (2 A   B )
5
2
M AB  ( B )
5
2(1) (2 - 73)
M BA  (2 B   A )
5
4
M BA  ( B )
5

Se evalúa el tramo B-C:

2(1) 20 * 3 2 (2 - 74)
M BC  (2 B   C ) 
3 12
4 2
M BC  ( B )  ( c )  15
3 3

2(1) 20 * 3 2 (2 - 75)
M CB  (2 C   B ) 
3 12
4 2
M CB  ( C )  ( B )  15
3 3

Se sigue con el tramo C-D y teniendo en cuenta la ecuación (2 - 69):

2(1) (2 - 76)
M CD  (2 C   D )
5
4
M CD  ( C )
5
2(1) (2 - 77)
M DC  (2 D   C )
5
2
M DC  ( C )
5

Se plantea un sistema de solución:

Edgar Eduardo Muñoz, IC, MIC 85

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

4 4 2 (2 - 78)
( B )  ( B )  ( c )  15  0
5 3 3
32 2
( B )  ( c )  15
15 3

Aplicando la ecuación (2 - 71), se tiene:

4 4 2 (2 - 79)
( C )  ( C )  ( B )  15  0
5 3 3
32 2
( C )  ( B )  15
15 3

Se despeja de (2 - 78)

 32 2 
 15  15 3   B   10.22
    2 32      
 15     C   10.22 
 3 15 

Se determina MAB y MDC:

2 2
M AB  ( B )  (10.22)  4.09kN .m
5 5
2
M DC  ( C )  4.09kN .m
5

Se determina MBA

4
M BA  ( B )  8.176kN .m
5

Edgar Eduardo Muñoz, IC, MIC 86

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

M AB=4.09 kN.m A MBA =8.176 kN.m

B
5m

RBA
RAB

Figura 2 - 63. Tramo A-B.

 4.09  8.176  R A * 5  0
R A  2.45kN

Por simetría de la viga a nivel de geometría y cargas:

R A  RD

Por equilibrio de cargas verticales, se tiene que:

 FY  0   R A  RD  20 * 3  RB  RC

Como las reacciones del centro son iguales, se tiene:

RB  RC  32.45kN

Basados en los resultados anteriores se realiza el diagrama de cortante y momento:

Edgar Eduardo Muñoz, IC, MIC 87

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Diagrama de cortante

40

30

20
Cortante(KN)

10

0
0 1 2 3 4 5 6 7 8 9 10 11 12 13

-10

-20

-30
X(m)

Figura 2 - 64. Diagrama de cortante

Diagrama de momento

10

0
0 2 4 6 8 10 12
Momento(KN.m)

-5

-10

-15

-20
x(m)

Figura 2 - 65. Diagrama de momento

Edgar Eduardo Muñoz, IC, MIC 88

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Resolver el punto anterior, considerando unas nuevas cargas puntuales de 10 KN en

la mitad de cada una de las tres luces, como se observa en la siguiente figura:

10 kN
1.5 m
10 kN W=20 kN/m 10 kN
2.5 m 2.5 m

A B C D

5.0 m
5.0 m 3.0 m

Figura 2 - 66. Viga en acero con empotramiento en los extremos. Fuente: Tomada de
referencia [8].

Como solamente se cambian las cargas externas y no el sistema estático, las

condiciones estructurales de la viga son las mismas expresada en las ecuaciones (2 -
68) a las (2 - 71) de este documento. Se comienza con el tramo de la viga A-B y
aplicando la condición de la ecuación (2 - 68), con lo cual se tiene:

2(1) 10 * 5 (2 - 80)
M AB  (2 A   B ) 
5 8
2
M AB  ( B )  6.25
5
2(1) 10 * 5 (2 - 81)
M BA  (2 B   A ) 
5 8
4
M BA  ( B )  6.25
5

Se evalúa el tramo B-C:

2(1) 20 * 3 2 10 * 3 (2 - 82)
M BC  (2 B   C )  
3 12 8
4 2
M BC  ( B )  ( c )  18.75
3 3

Edgar Eduardo Muñoz, IC, MIC 89

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

2(1) 20 * 3 2 10 * 3 (2 - 83)
M CB  (2 C   B )  
3 12 8
4 2
M CB  ( C )  ( B )  18.75
3 3

Se sigue con el tramo C-D y teniendo en cuenta la ecuación (2 - 69):

2(1) 10 * 5 (2 - 84)
M CD  (2 C   D ) 
5 8
4
M CD  ( C )  6.25
5
2(1) 10 * 5 (2 - 85)
M DC  (2 D   C ) 
5 8
2
M DC  ( C )  6.25
5

Aplicando la ecuación (2 - 70), se tiene:

4 4 2 (2 - 86)
( B )  6.25  ( B )  ( c )  18.75  0
5 3 3
32 2
( B )  ( c )  12.5
15 3

Aplicando la ecuación (2 - 71), se tiene:

4 2 4 (2 - 87)
( C )  ( B )  18.75  ( C )  6.25  0
3 3 5
32 2
( C )  ( B )  12.5
15 3

Se despeja de (2 - 78)

Edgar Eduardo Muñoz, IC, MIC 90

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 32 2 
 12.5  15 3   B   8.522
      
 12.5   2 32   C   8.522 
 3 15 

Se determina MAB y MDC:

2 2
M AB  ( B )  6.25  (8.52)  6.25  2.8412kN .m
5 5
2 2
M DC  ( C )  6.25  (8.52)  6.25  2.8412kN .m
5 5

Se determina MBA

4 4
M BA  ( B )  6.25  (8.52)  6.25  13.067kN.m
5 5

10 kN

2.5 m

MAB=2,84 kN.m A MBA=13.066 kN.m

B
5m

RA RBA

Figura 2 - 67. Sistema estático del tramo A-B.

Se realiza equilibrio, haciendo sumatoria de momentos con respecto al nudo A:

Edgar Eduardo Muñoz, IC, MIC 91

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

2.8412  13.067  R A * 5  10 * 2.5  0

R A  2.954kN

Por simetría de la viga a nivel de geometría y cargas:

R A  RD

Por equilibrio de cargas verticales y teniendo en cuenta la simetría a nivel de geometría y

cargas de la viga, se tiene que:

 FY  0 R A  RD  20 * 3  RB  RC  10  10  10

Como las reacciones del centro son iguales, se tiene:

RB  RC  42.046kN

Basados en los resultados anteriores se realiza el diagrama de cortante y momento:

Edgar Eduardo Muñoz, IC, MIC 92

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

10 kN
1.5 m
10 kN W=20 kN/m 10 kN
2.5 m 2.5 m

A B C 5.0 m D
5.0 m 3.0 m

35
7.046
5
2.954
Cortante
(kN)

7.046 5 2.954

35

13.07 13.07

2.8412
Momento 2.8412
(kN.m)

4.5438
4.5438
16.93

Figura 2 - 68. Diagrama de cortante y momento

[Link].2 Método de distribución de momentos

Es un método ideado por Hardy Cross en 1929 en la Universidad de Illinois, basado

en el de giro -* deflexión que se basa en procesos iteractivos hasta lograr determinar
las reacciones y las fuerzas internas de cualquier viga estáticamente determinada o
indeterminada. Para entender su aplicación se resuelve a continuación la viga que se
muestra en la Figura 2 - 66 por el este método.

Se determinan inicialmente los coeficientes de distribución que se basan en el análisis

de la rigidez relativa. Considerando EI=1.0 se tiene:

Edgar Eduardo Muñoz, IC, MIC 93

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

CD AB  0
1
CD BA  5 0.375
1 1

3 5
1
CD BC  3  0.625
1 1

3 5
1
CDCB  3  0.625
1 1

3 5
1
CDCD  3  0.375
1 1

3 5
CD DC 0

Después se determinan los momentos de empotramientos de cada luz:

M ABf   M ABf  6.25kN .m

M CDf   M DCf  6.25kN .m

M BCf   M CBf  18.75kN .m

Se aplica con esto datos el proceso iterativo, obteniéndose los siguientes resultados:

Edgar Eduardo Muñoz, IC, MIC 94

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

1.5 m
10 kN W=20 kN/m 10 kN
2.5 m 2.5 m

A B C D

5.0 m
5.0 m 3.0 m

0.000 0.375 0.625 0.625 0.375 0.000

6.250 -6.250 18.750 -18.750 6.250 -6.250
0.000 -4.688 -7.813 7.813 4.688 0.000
-2.344 0.000 3.906 -3.906 0.000 2.344
0.000 -1.465 -2.441 2.441 1.465 0.000
-0.732 1.221 -1.221 0.732
-0.458 -0.763 0.763 0.458
-0.229 0.381 -0.381 0.229
-0.143 -0.238 0.238 0.143
-0.072 0.119 -0.119 0.072
-0.045 -0.075 0.075 0.045
-0.022 0.037 -0.037 0.022
-0.014 -0.023 0.023 0.014
-0.007 0.012 -0.012 0.007
-0.004 -0.007 0.007 0.004
-0.002 0.004 -0.004 0.002
-0.001 -0.002 0.002 0.001
2.842 -13.068 13.068 -13.068 13.068 -2.841901

Figura 2 - 69. Solución de la viga por el método de Cross

Los momento obtenidos en los nudos A,B,C y D coinciden con los momento
obtenidos en el ejercicio anterior resuelto para este misma viga pero por el método de
giro deflexión, por lo que le diagrama de cortante y momento es el que se presenten
en la Figura 2 - 68.

Edgar Eduardo Muñoz, IC, MIC 95

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

10 kN
1.5 m
10 kN W=20 kN/m 10 kN
2.5 m 2.5 m

A B C D

5.0 m
5.0 m 3.0 m

0.000 0.375 0.625 1.000 0.000

6.250 -6.250 18.750 -18.750 25.000
0.000 -4.688 -7.813 -6.250 0.000
-2.344 0.000 -3.125 -3.906 0.000
0.000 1.172 1.953 3.906 0.000
0.586 1.953 0.977
-0.732 -1.221 -0.977 0.000
-0.366 -0.488 -0.610
0.183 0.305 0.610 0.000
0.092 0.305 0.153
-0.114 -0.191 -0.153 0.000
-0.057 -0.076 -0.095
0.029 0.048 0.095 0.000
0.014 0.048 0.024
-0.018 -0.030 -0.024 0.000
-0.009 -0.012 -0.015
0.004 0.007 0.015 0.000
4.166 -10.414 10.414 -25.000 25.000
Figura 2 - 70. Solución de la viga por el método de Cross

Se procede a determinar las reacciones de cada apoyo basados en los resultados

obtenidos de momento en cada uno de los nudos. Para ello, se deben resolver los
siguientes tramos de la viga:

Edgar Eduardo Muñoz, IC, MIC 96

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

10 kN
M AB=4.166 kN.m 1.5 m
10 kN W=20 kN/m
M23=10.414 kN.m
2.5 m
M BA =25 kN.m
M BA =10.414 kN.m

A B C
B

5.0 m 3.0 m

RAB RBA
RBC RCB

Figura 2 - 71. Resultados de momentos en los nudos en cada tramo de viga

Del primer tramo se viga se tiene:

 10 * 2.5  4.166  10.414  RBA * 5  0

RBA  6.24kN
R AB  3.75kN

Del segundo tramo:

 10 *1.5  25  10.414  20 * 3 *1.5  RCB * 3  0

RCB  39.86kN
RBC  30.138kN

De lo anterior se tiene que las reacciones son:

RA 3.750 kN
MA 4.166 kN.m
RB 36.38 kN
Rc 49.86 kN

Definidas las reacciones se construye diagrama de cortante y momento:

Edgar Eduardo Muñoz, IC, MIC 97

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Diagrama de cortante

40

30

20

10
Cortante(KN)

0
0 1 2 3 4 5 6 7 8 9 10 11 12 13

-10

-20

-30

-40

-50
X(m)

Diagrama de momento

30

25

20

15
Momento(KN.m)

10

0
0 2 4 6 8 10 12

-5

-10

-15
x(m)

Figura 2 - 72. Diagrama de momento

[Link].3 Método matricial

Edgar Eduardo Muñoz, IC, MIC 98

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

[Link].3.1 Fundamentos

La solución por el método matricial se basa en la ecuación (2 - 8) deducida en este

documento. Cada nudo tiene dos grados de libertad o restringidos, es decir
posibilidades de desplazamientos en el sentido “Y” y rotación en el sentido “Z”. Por
lo anterior, es necesario deducir o determinar la matriz de rigidez global de dicho
elemento arbitrariamente orientado. Para ello, se debe determinar la matriz de rigidez
global del elemento, es decir cuando el elemento este solamente sobre el eje “X” de
un sistema plano de coordenadas, como se observa a continuación:

Figura 2 - 73. Elemento de viga prismático sometidos a fuerzas cortantes y

momentos.

Se procede a deducir la matriz de rigidez global que tiene matricialmente lo siguiente:

 Fiy   K11 K12 K13 K14  U iy 

M    
 iz   K 21 K 22 K 23 K 24    iz 
   
 F jy   K 31 K 32 K 33 K 34  U jy 
M jz   K 41 K 42 K 43 K 44    jz 
Figura 2 - 74. Matriz de rigidez global de un elemento viga.

Para ello, se emplea el método de giro deflexión. Se establece una deformación en

Uiy =1 y las demás deformaciones iguales a cero, obteniéndose:

Edgar Eduardo Muñoz, IC, MIC 99

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 Fiy   K11 K12 K13 K14  U iy  1.0  Fiy   K11 

M   K    
 iz   21 K 22 K 23 K 24    iz  0   M iz  K 21 
     
 F jy   K 31 K 32 K 33 K 34   U jy  0   F jy   K 31 
M iz   K 41 K 42 K 43

K 44    jz  0  M jz  K 41 

Figura 2 - 75. Condición de frontera 1: Deformación en Uiy =1 y las demás

deformaciones iguales a cero.
2 EI  U iy  6 EI
M iz  (0  0  3  )  2  K 21
L  L  L
2 EI  U  6 EI
M jz  (0  0  3  iy )  2  K 41
L  L  L

12EI
M iz  M jz  F jy L  0  F jy 
L3
12EI
 Fiy 
L3
De esta formase se ha deducido la primera columna de la matriz de rigidez del
elemento viga:

 12EI 
 3
K12 K13 K14 
 Fiy   L  U iy 
 M   6 EI K 22 K 23 K 24    
 iz   L2   iz 
   12EI
 F jy   3 K 32 K 33 K 34   jy 
U
M jz   L  
 6 EI   jz 
 K 42 K 43 K 44 
 L2 

Edgar Eduardo Muñoz, IC, MIC 100

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 76. Condición de frontera 1: Deducción de la primera columna de la

matriz de rigidez del elemento viga:
.

Posteriormente se considera una deformación en iz =1 y las demás deformaciones

iguales a cero, obteniéndose:

 Fiy   K11 K12 K13 K14   U iy  0   Fiy   K12 

M   K    
 iz   21 K 22 K 23 K 24   iz  1.0  M iz  K 22 
     
 F jy   K 31 K 32 K 33 K 34   U jy  0   F jy   K 32 
M iz   K 41 K 42 K 43

K 44    jz  0  M jz  K 42 

Figura 2 - 77. Condición de frontera 2: Deformación en iz =1 y las demás

deformaciones iguales a cero.

Se determinan los valores correspondientes con giro –deflexión:

2 EI 4 EI
M iz  (2 iz  0  0)   K 22
L L
2 EI 2 EI
M jz  (0   iz  0)   K 42
L L

Edgar Eduardo Muñoz, IC, MIC 101

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

6 EI
M iz  M jz  F jy L  0  F jy 
L2
6 EI
 Fiy 
L2

De esta formase se ha deducido la segunda columna de la matriz de rigidez del

elemento viga:

 12EI 6 EI 
 3
K13 K14 
 Fiy   6LEI L2
4 EI  U iy 
M   K 23 K 24    
 iz   L2 L   iy 
   12EI
 F jy   3 K 34   jy 
6 EI U
 2 K 33
M iz   L L  
 6 EI 2 EI   jz 
 K 43 K 44 
 L2 L 

Figura 2 - 78. Condición de frontera 2: Deducción de la segunda columna de la

matriz de rigidez del elemento viga:

Posteriormente se considera una deformación en Ujy =1 y las demás deformaciones

iguales a cero, obteniéndose:

Edgar Eduardo Muñoz, IC, MIC 102

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 Fiy   K11 K12 K13 K14   U iy  0   Fiy   K13 

M   K    
 iz   21 K 22 K 23 K 24    iz  0   M iz   
K 23 
     
 F jy   K 31 K 32 K 33 K 34  U jy  1.0  F jy   K 33 
  
M iz   K 41 K 42 K 43 K 44    jz  0 
 
M jz  
K 43 

Figura 2 - 79. Condición de frontera 3: Deformación en Ujy=0 y las demás
deformaciones iguales a cero.

Se determinan los valores correspondientes con giro –deflexión:

2 EI  U jy  6 EI
M iz  (0  0  3 )   2  K 23
L  L  L
2 EI U  6 EI
M jz  (0  0  3 jy )   2  K 43
L  L  L

12EI
M iz  M jz  F jy L  0  F jy 
L3
12EI
 Fiy 
L3

De esta formase se ha deducido la segunda columna de la matriz de rigidez del

elemento viga:

Edgar Eduardo Muñoz, IC, MIC 103

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 12EI 6 EI 12EI 
 3
 K14 
 Fiy   L L2 L3
 U iy 
 M   6 EI 4 EI 6 EI
 2 K 24    
 iz   L2 L L  iz 
   12EI  
 F jy   3 K 34   jy 
6 EI 12EI U
 2
  L   
M jz 
L L3
  6 EI 2 EI 6 EI   jz 
  2 K 44 
 L2 L L 

Figura 2 - 80. Condición de frontera 3: Deducción de la tercera columna de la matriz

de rigidez del elemento viga:

Posteriormente se considera una deformación en  jz =1 y las demás deformaciones

iguales a cero, obteniéndose:

Edgar Eduardo Muñoz, IC, MIC 104

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 Fiy   K11 K12 K13 K14   U iy  0   Fiy   K14 

M   K    
 iz   21 K 22 K 23 K 24    iz  0   M iz   
K 24 
     
 F jy   K 31 K 32 K 33 K 34   U jy  0   F jy   K 34 
  
M iz   K 41 K 42 K 43 K 44   jz  1.0
 
M jz  
K 44 


Figura 2 - 81. Condición de frontera 4: Deformación en  jz =1 y las demás

deformaciones iguales a cero.

Se determinan los valores correspondientes con giro –deflexión:

2 EI 2 EI
M iz  (0   jz  0)   K 24
L L
2 EI 2 EI
M jz  (2 jz  0  0)   K 44
L L

De esta formase se ha deducido la cuarta columna de la matriz de rigidez del

elemento viga:

 12EI 6 EI 12EI 6 EI 
 
 Fiy   L
3
L2 L3 L2  U 
 M   6 EI 4 EI 6 EI 2 EI   iy 
 2 
 iz   L2 L L L    iz 
   12EI 6 EI U 
 F jy   3  2   jy 
6 EI 12EI
 2
M jz   L L  
3
L L
 6 EI 2 EI 6 EI 4 EI   jz 
  2 
 L2 L L L 

Figura 2 - 82. Condición de frontera 2: Deducción de la cuarta columna de la matriz

de rigidez del elemento viga:

Edgar Eduardo Muñoz, IC, MIC 105

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

La solución por el método de los desplazamientos se explica en las ecuaciones (2 - 8)

a la (2 - 12).

[Link].3.2 Ejemplo: Solución de una viga continúa de dos luces iguales

Figura 2 - 83. Sistema estático de la viga

Inicialmente se construye la matriz de rigidez de cada elemento:

 12EI 6 EI 12EI 6 EI 
 
 F1 y   6LEI
3
L2 L3 L2  U 
M   4 EI 6 EI 2 EI   1 y 
 2 
 1z   L2 L L L   1 y 
   12EI 6 EI  U 2 y 
 F2 y   3
6 EI 12EI
 2   
M 2 z   L L L3 L2    
 6 EI 2 EI 6 EI 4 EI   2 z 
  2 
 L2 L L L 

Edgar Eduardo Muñoz, IC, MIC 106

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 12EI 6 EI 12EI 6 EI 
 
 F2 y   6LEI
3
L2 L3 L2  U 
M   4 EI 6 EI 2 EI   2 y 
 2 
 2 z   L2 L L L   2 y 
   12EI 6 EI  U 3 y 
 F3 y   3
6 EI 12EI
 2   
 M 3 z   L L L3 L2    
 6 EI 2 EI 6 EI 4 EI   3 z 
  2 
 L2 L L L 

Se plantea el sistema de solución por el método de los desplazamientos de la

estructura:

 4 EI 2 EI 
 L 0   
L  WL2 
 
 M 1z  0   2 EI 8 EI 2 EI
    O1z  ?   12 
 M  0  L  
L L  O  ?   0 2 
 2z    2z
    O
2 EI 4 EI WL 
 M 3 z  0   0
L L  3z  ?   
    6 EI     12 
 F1 y  ?   2  0   WL 
6 EI
0    U1 y
 F2 y  ?   L L2  U 2 z  0  2 
   6 EI    WL 
   U 3 z
6 EI
 F3 y  ?   2 0
  0   WL 
L L2
 6 EI 6 EI   2 
 0   2   
 L2 L 

Para la solución de las deformaciones se tiene el siguiente sistema de ecuaciones

lineales:
1
 4 EI 2 EI 
 0   WL2 
O1z   L L  0  
  2 EI 8EI 2 EI   12 
O2 z      00 2 
O   L L L   WL 
 3z  
0
2 EI 4 EI  0  ( 12 )
 L L   

Para la solución se aplica el siguiente código en MATLAB:

Edgar Eduardo Muñoz, IC, MIC 107

 Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

L sym('L')
EI sym('EI')
sym(' ')
c EI/L, EI/L,0 EI/L, EI/L, EI/L 0, EI/L, EI/L
- L /1 0 L /1
A INV(c)

El resultado es el siguiente:

 WL3 
 O1z   48EI 
   
O2 z    0 3 
O   WL 
 3z   
 48EI 

Se determinan las fuerzas internas, empleando la matriz de rigidez de cada elemento,

las deformaciones de cada nudo y los valores de los vectores de empotramiento:

 12EI 6 EI 12EI 6 EI 
   3 
 Fy (1 2)   L
3
L2 L3 L2   0  WL /22   WL 

2 EI  WL3   WL   8
M   6 EI 4 EI 6 EI
 2  0 
 z (1 2)   L2 L L L   48EI    12    5
   12EI 6 EI   0  WL / 2   WL 
  3
F 6 EI 12EI
 y ( 2 1)  2     2  8 
M z ( 21)   L L L3 L2   0   WL   1 2 
 6 EI 2 EI 6 EI 4 EI     12   WL 
  2   8 
 L2 L L L 

 12EI 6 EI 12EI 6 EI 
  5 
 Fy ( 23)   L
3
L2 L3 L2   0  WL /22   WL 
M  6 EI 4 EI 6 EI 2 EI     WL   1
8

 2
 z (33)   L2 L L L 
 0     WL2

   12EI  0    12    8 
 Fy (3 2)   3
6 EI 12EI 6 EI  WL / 2
 2  2  WL 3     3 
M z (3 2)   L L 
2
L L3   WL   8 WL 
 6 EI 2 EI 6 EI 4 EI   48EI   12   
  2   0 
 L2 L L L 

Edgar Eduardo Muñoz, IC, MIC 108

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Con dicha información se construye el diagrama de cortante y momento:

R2=1.25WL
R1=0.375WL
L L R31=0.375WL
0.625WL

R1

0.625WL

0.375L 0.625L 0.625L 0.375L

R1

0.125WL2

0.070WL2 0.070WL2

Figura 2 - 84. Diagrama de cortante y momento

[Link].3.2 Ejemplo: Análisis estructural de un puente peatonal de dos luces

Este ejercicio consiste en el análisis estructural de un puente peatonal continuo con el

sistema estático que se presenta en la Figura 2 - 84 y en donde L=25 m. La NSR-10
establece alturas mínimas de vigas en función de la luz y del sistema estático
(continua, en voladizo o simplemente apoyada). Para este caso se optó por una altura
de la viga de 2 m. La sección transversal es:

Edgar Eduardo Muñoz, IC, MIC 109

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

2.0m

0.10 m

1.90 m

0.30 m
Figura 2 - 85. Sección transversal de viga

Las propiedades geométricas son las siguientes:

2
Figura h b Área Y AY I Ad
2 4
(m) (m) (m ) (m) (m )
1 0.1 2.0 0.2 2.0 0.4 0.00 0.11
2 1.9 0.3 0.6 1.0 0.5 0.17 0.04
0.770 0.9 0.17 0.15

Tabla 2 - 13. Propiedades geométricas de la sección de la viga

Se hace una evaluación de cargas. Con respecto a la carga muerta se tiene:

Edgar Eduardo Muñoz, IC, MIC 110

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Descripción Calculo Carga muerta (kN/m)

Peso propio 0.77x24 18.48

Acabado 0.03x2x22 1.32
barandas 3

22.8

Tabla 2 - 14. Carga muerta por metro lineal

De acuerdo al Código Colombiano de Diseño Sísmico de puentes (CCDSP), la carga

viva de diseño para puentes peatonales es de 4 kN/m2 es decir:

cv1
KN/m

2.1.4 Pórtico plano

[Link] Definición y antecedentes

La vigas es un sistema estructural de elementos horizontales sometidos generalmente

a cortante y momento que deben soportar cargas externas aplicadas en sus nudos y/o
sobre los elementos.

[Link] Métodos de análisis

[Link].1 Método matricial

[Link].1.1 Fundamentos

Se comienza con deducir l matriz de rigidez local de un elemento de pórtico

localizado en el eje x:

Edgar Eduardo Muñoz, IC, MIC 111

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 86. Elemento de viga prismático sometidos a fuerzas cortantes, momentos

y cargas axiales.

La matriz de rigidez de este elemento tiene tres grados de liberta por nudo y se
expresa de la siguiente forma:

 Fix   K11 K12 K13 K14 K15 K16  U ix 

F    
 iy   K 21 K 22 K 23 K 24 K 25 K 26  U iy 
 M iz   K 31 K 32 K 33 K 34 K 35 K 36    iz 
   
 F jx   K 41 K 42 K 43 K 44 K 45 K 46  U jx 
 F jy   K 51 K 52 K 53 K 54 K 55 K 56  U jy 
    
M jz   K 61 K 62 K 63 K 64 K 65 K 66    jz 

Para su deducción se emplean las realizadas a los elementos de viga y armadura plana
explicados anteriormente:

Edgar Eduardo Muñoz, IC, MIC 112

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 87. Esquema para deducir la matriz de rigidez local de un elemento de

pórtico:

De dicho esquema se obtiene:

 AE AE 
 L 0 0  0 0 
L
 6 EI  U ix 
 Fix   0 12EI 6 EI
0 
12EI

F   L3 L2 L3 L2  U 
 iy   6 EI 4 EI 6 EI 2 EI   iy  (2 - 88)
 M iz   0  2  
L    iz 
0
    AE L2 L L
AE  U
 F jx   0 0 0 0   jx 
 F jy   L L  U 
  12EI 6 EI 12EI 6 EI   jy 
M jz   0   0  2   jz 
L  
 L3 L2 L3
 6 EI 2 EI 6 EI 4 EI 
 0 0  2
L2 L L L 

Edgar Eduardo Muñoz, IC, MIC 113

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Igualmente se puede deducir que para el elemento vertical (columnas) la matriz de

rigidez local es la siguiente:

Figura 2 - 88. Elemento vertical (columna)

 12EI 6 EI 12EI 6 EI 
 L3 0   0 
L2 L3 L2 
 
 Fix   0 AE
0 0 
AE
0  U ix 
F   L L  U 
 iy   6 EI 4 EI 6 EI 2 EI   iy 
 M iz    2
L    iz 
0 0  
    12LEI L L2
F 6 EI 12EI 
6 EI U jx (2 - 89)
 jx   3 0 0  
 F jy   L L2 L3 L2  U jy 
   0 
AE AE
0    jz 
M jz  
0 0
L L  
 6 EI 2 EI 6 EI 4 EI 
 0 0
 L2 L L2 L 

Ahora se hace la deducción para obtener la matriz de rigidez global que corresponde
a un elemento primatico arbitrariamente orientado:

Edgar Eduardo Muñoz, IC, MIC 114

 Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 89. Elemento arbitrariamente orientado

Para su determinación es necesario utilizar la ecuación deducida anteriormente que

depende de la matriz de rigidez local en el eje x y de la matriz de transformación (ver
(2 - 24):

 AE AE 
 L 0 0  0 0 
L
 
 Cos Sin   Cos Sin
T
0 0 0 0 12EI 6 EI 12EI 6 EI 0 0 0 0
 0 0  3
 Sin Cos 0 0 0 0  L3 L2 L L2   Sin Cos 0 0 0 0
  0 6 EI 4 EI 6 EI 2 EI 
 0 0  2
 0 0
0
K GLOBAL   
0 1 0 0  L2 L L L 0 1 0 0
(2 - 90)
  AE  
 0 0 0 Cos Sin 0 AE
0  0 0 0 Cos Sin 0
 0 0 0
 0 0 0  Sin Cos 0  L L  0 0 0  Sin Cos 0
   0 
12EI

6 EI 12EI
 2   0
6 EI 
 0 0 0 0 0 1
 L3 L2
0
L3 L 
0 0 0 0 1
 6 EI 2 EI 6 EI 4 EI 
 0 0  2
L2 L L L 

[Link].1.2 Ejemplo : Determinación de rigidez de un pórtico

A continuación se presenta un ejemplo para su aplicación. El primer ejercicio consiste

en calcular la rigidez (K) considerando una resistencia del concreto de 21 Mpa.

Edgar Eduardo Muñoz, IC, MIC 115

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 90. Pórtico plano con una carga puntual paramétrica horizontal.

Se determina el módulo de elasticidad:

Ec  3900 fc  3900 21  17872.045Mpa

Después se establecen las coordenadas, restricciones y fuerzas en los nudos

, lo cual se presenta a continuación:

Elemento Ni Nf A(m2) I(m3) L(m)

1 3 1 0.09 0.00068 3.0

2 1 2 0.12 0.0016 6.0

3 4 2 0.09 0.00068 3.0

Tabla 2 - 15. Información de los nudos

Posteriormente se determinan los datos de los elementos:

Edgar Eduardo Muñoz, IC, MIC 116

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Elemento Ni Nf A(m2) I(m3) L(m)

1 3 1 0.09 0.00068 3.0

2 1 2 0.12 0.0016 6.0

3 4 2 0.09 0.00068 3.0

Tabla 2 - 16. Información de los elementos

Luego se determina la matriz de rigidez global de cada elemento del pórtico:

5361.6 0.0 -8042.4 -5361.6 0.0 -8042.4

0.0 536161.4 0.0 0.0 -536161.4 0.0
-8042.4 0.0 16084.8 8042.4 0.0 8042.4
-5361.6 0.0 8042.4 5361.6 0.0 8042.4
0.0 -536161.4 0.0 0.0 536161.4 0.0
-8042.4 0.0 8042.4 8042.4 0.0 16084.8

Figura 2 - 91. Matriz de rigidez elemento 3-1

u1x u1y O1z u2x u2y O2z

357440.9 0.0 0.0 -357440.9 0.0 0.0 u1x
0.0 1588.6 4765.9 0.0 -1588.6 4765.9 u1y
0.0 4765.9 19063.5 0.0 -4765.9 9531.8 O1z
-357440.9 0.0 0.0 357440.9 0.0 0.0 u2x
0.0 -1588.6 -4765.9 0.0 1588.6 -4765.9 u2y
0.0 4765.9 9531.8 0.0 -4765.9 19063.5 O2z

Figura 2 - 92. Matriz de rigidez elemento 1-2

Edgar Eduardo Muñoz, IC, MIC 117

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

5361.6 0.0 -8042.4 -5361.6 0.0 -8042.4

0.0 536161.4 0.0 0.0 -536161.4 0.0
-8042.4 0.0 16084.8 8042.4 0.0 8042.4
-5361.6 0.0 8042.4 5361.6 0.0 8042.4
0.0 -536161.4 0.0 0.0 536161.4 0.0
-8042.4 0.0 8042.4 8042.4 0.0 16084.8

Figura 2 - 93. Matriz de rigidez elemento 4-2

Posteriormente se determina la matriz de rigidez de la estructura a través de un

proceso de ensamblaje y se plantea el sistema de solución:

Edgar Eduardo Muñoz, IC, MIC 118

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Ux1 Uy1 z1 Ux2 Uy2 z2 Ux3 Uy3 z3 Ux4 Uy4 z4

Fx1=10 362842.2 0.0 8102.0 -357440.9 0.0 0.0 -5401.3 0.0 8102.0 0.0 0.0 0.0 Ux1=?

Fy1=0 0.0 537750.0 4765.9 0.0 -1588.6 4765.9 0.0 -536161.4 0.0 0.0 0.0 0.0 Uy1=?

Mz1=0 8102.0 4765.9 35267.5 0.0 -4765.9 9531.8 -8102.0 0.0 8102.0 0.0 0.0 0.0 z1=?

Fx2=0 -357440.9 0.0 0.0 362842.2 0.0 8102.0 0.0 0.0 0.0 -5401.3 0.0 8102.0 Ux2=?

Fy2=0 0.0 -1588.6 -4765.9 0.0 537750.0 -4765.9 0.0 0.0 0.0 0.0 -536161.4 0.0 Uy2=?

Fz2=0 0.0 4765.9 9531.8 8102.0 -4765.9 35267.5 0.0 0.0 0.0 -8102.0 0.0 8102.0 z2=?

Fx3=? = -5401.3 0.0 -8102.0 0.0 0.0 0.0 5401.3 0.0 -8102.0 0.0 0.0 0.0 * Ux3=0

Fy3=? 0.0 -536161.4 0.0 0.0 0.0 0.0 0.0 536161.4 0.0 0.0 0.0 0.0 Uy3=0

Mz3=? 8102.0 0.0 8102.0 0.0 0.0 0.0 -8102.0 0.0 16204.0 0.0 0.0 0.0 z=0

Fx4=? 0.0 0.0 0.0 -5401.3 0.0 -8102.0 0.0 0.0 0.0 5401.3 0.0 -8102.0 Ux4=0

Fy4=? 0.0 0.0 0.0 0.0 -536161.4 0.0 0.0 0.0 0.0 0.0 536161.4 0.0 Uy4=0

Mz4=? 0.0 0.0 0.0 8102.0 0.0 8102.0 0.0 0.0 0.0 -8102.0 0.0 16204.0 z4=0

 F1x  10  U1x  ?  0 
 F1y  0  U  ?  0 
   1y   
 F1z  0   1z  ?   0 
     
 F2x  0  U1x  ?  0 
 F2y  0  U jy  ?  0 
     
 F2z  0   K nn K na    jz  ?   0 
 F3x  ?      
   K an K aa   0   0 
 F3y  ?   0  0 
     
 F3z  ?   0  0 
 F4x  ?   0  0 
     
 F4y  ?   0  0 
 F4z  ?   0  0 
     

Aplicando la ecuación (2 - 8) se determinan las deformaciones de los dos nudos con

grados de libertad:

-1
Ux1 362802.5 0.0 8042.4 -357440.9 0.0 0.0 F1x=10 0.00128623 m
Uy1 0.0 537750.0 4765.9 0.0 -1588.6 4765.9 F1y=0 0.00000408 m
0z1 8042.4 4765.9 35148.4 0.0 -4765.9 9531.8 F1z=0 -0.00023333 Radianes
Ux2 = -357440.9 0.0 0.0 362802.5 0.0 8042.4 * F2x=0 = 0.00127230 m
Uy2 0.0 -1588.6 -4765.9 0.0 537750.0 -4765.9 F2y=0 -0.00000408 m
0z2 0.0 4765.9 9531.8 8042.4 -4765.9 35148.4 F2z=0 -0.00022895 Radianes

Figura 2 - 94. Determinación de las deformaciones

Edgar Eduardo Muñoz, IC, MIC 119

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Con ello se determina la rigidez horizontal del pórtico, que representa que es
necesario una fuerza de 7774.7 kN para que se deforme un 1 metro:

P 10
K   7774.7kN / m  777.47ton / m
Ux1 0.001286

Igualmente se determinan las reacciones aplicando la ecuación (2 - 9)

F3x -5.0
F3y -2.2
F3z 8.5
F4x -5.0
F4y 2.2
F4z 8.4

Se revisa equilibrio en cada dirección:

 Fy  0
 2.2  2.2  0 _ OK
 Fx  0
10  5  5  0 _ OK

Con la rigidez de un pórtico y su masa se puede determinar el periodo de vibración

empleado el cconcepto básicos de dinámica estructural .Para un grado de libertad, la
Ecuación de equilibrio, que es diferencial lineal homogénea, es la siguiente:

..
m X  KX  0 (2 - 91)

Donde:

Edgar Eduardo Muñoz, IC, MIC 120

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

m : Masa (ton-s2/m)
X : Aceleración de la masa (m/s2)
X : Desplazamiento relativo (m)
K : Rigidez (ton/m)

De esta ecuación se obtiene la frecuencia natural del sistema en radianes

por segundo (rad/seg)

k (2 - 92)

m

También la ffrecuencia natural del sistema en ciclos por segundo (1/seg o Hertz)


f  (2 - 93)
2

Igualmente el periodo natural del sistema en segundos (seg)

2
T (2 - 94)

Al aplicarlo para el caso con una masa de 1.77 ton- s2/m:

k 777.47
   20.9 Rad / seg
m 3.3

2 2 * 3.1416
T   0.30seg
 20.9

Edgar Eduardo Muñoz, IC, MIC 121

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

[Link].1.3 Ejemplo: Resolver un pórtico plano con carga horizontal y vertical

Otro ejercicio consiste en calcular deformaciones, reacciones y fuerzas internas

considerando P=50 kN y W=35 kN/m. (La matriz de rigidez es la misma, cambia el
vector de fuerzas y de empotramiento.)

Figura 2 - 95. Pórtico plano con una carga puntual horizontal y una carga distribuida.

Se determina inicialmente los momentos de empotramiento:

F1xf  0
WL
F1 yf   105KN
2
WL2
M 1fz   105KN .m
12
F2f  0
WL
F2fy   105KN
2
WL2
M 2fz   105KN .m
12

Se plantea posteriormente el sistema de solución:

Edgar Eduardo Muñoz, IC, MIC 122

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 F1x  50  U1x  ?   0 
 F1y  0  U  ?   105 
   1y   
 F1z  0   1z  ?   105 
     
 F2x  0  U 2 x  ?  0 
 F2y  0  U 2 y  ?  105 
     
 F2z  0   K nn K na    2 z  ?    105 
 F3x  ?      
   K an K aa   0   0 
 F3y  ?   0   0 
     
 F3z  ?   0   0 
 F4x  ?   0   0 
     
 F4y  ?   0   0 
 F4z  ?   0   0 
     

Figura 2 - 96. Matriz de rigidez de la estructura y el sistema de solución

Se determina las deformaciones y la elástica:

Edgar Eduardo Muñoz, IC, MIC 123

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Ux1 0.00647708 m
Uy1 -0.00017541 m
0z1 -0.00527996 Radianes
Ux2 0.00631555 m
Uy2 -0.00021626 m
0z2 0.00296857 Radianes

Figura 2 - 97. Deformaciones y la elástica

Después las reacciones en los apoyos:

F3x 7.7 kN
F3y 94.0 kN
F3z 9.6 kN.,m
F4x -57.7 kN
F4y 116.0 kN
F4z 74.7 kN.,m

Posteriormente las fuerzas internas:

Edgar Eduardo Muñoz, IC, MIC 124

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

FX(3 - 1) 5361.6 0.0 -8042.4 -5361.6 0.0 -8042.4 0.0000000 0 7.74

FY(3 - 1) 0.0 536161.4 0.0 0.0 -536161.4 0.0 0.0000000 0 94.05

FZ(3 - 1) = -8042.4 0.0 16084.8 8042.4 0.0 8042.4 0.0000000 `+ 0 = 9.63

FX(1 - 3) -5361.6 0.0 8042.4 5361.6 0.0 8042.4 0.0064771 0 -7.74
FY(1 - 3) 0.0 -536161.4 0.0 0.0 536161.4 0.0 -0.0001754 0 -94.05
FZ(1 - 3) -8042.4 0.0 8042.4 8042.4 0.0 16084.8 -0.0052800 0 -32.84

FX(1 - 2) 357440.9 0.0 0.0 -357440.9 0.0 0.0 0.0064771 0 57.74

FY(1 - 2) 0.0 1588.6 4765.9 0.0 -1588.6 4765.9 -0.0001754 105 94.05
FZ(1 - 2) 0.0 4765.9 19063.5 0.0 -4765.9 9531.8 -0.0052800 105 32.84
FX(2 - 1) -357440.9 0.0 0.0 357440.9 0.0 0.0 0.0063156 0 -57.74
FY(2 - 1) 0.0 -1588.6 -4765.9 0.0 1588.6 -4765.9 -0.0002163 105 115.95
FZ(2 - 1) 0.0 4765.9 9531.8 0.0 -4765.9 19063.5 0.0029686 -105 -98.54

FX(4 - 2) 5361.6 0.0 -8042.4 -5361.6 0.0 -8042.4 0.0000000 0 -57.74

FY(4 - 2) 0.0 536161.4 0.0 0.0 -536161.4 0.0 0.0000000 0 115.95

FZ(4 - 2) = -8042.4 0.0 16084.8 8042.4 0.0 8042.4 0.0000000 `+ 0 = 74.67

FX(2 - 4) -5361.6 0.0 8042.4 5361.6 0.0 8042.4 0.0063156 0 57.74
FY(2 - 4) 0.0 -536161.4 0.0 0.0 536161.4 0.0 -0.0002163 0 -115.95
FZ(2 - 4) -8042.4 0.0 8042.4 8042.4 0.0 16084.8 0.0029686 0 98.54

Con ellas se puede construir el diagrama axial, cortante y momento del pórtico:

57.74

FX(1 - 2) 57.74
FY(1 - 2) 94.05
FZ(1 - 2) = 32.84
FX(2 - 1) -57.74
FX(3 - 1) 7.74
FY(2 - 1) 115.95 FX(4 - 2) -57.74
FY(3 - 1) 94.05
FZ(2 - 1) -98.54 FY(4 - 2) 115.95
FZ(3 - 1) = 9.63
FZ(4 - 2) = 74.67
FX(1 - 3) -7.74
FX(2 - 4) 57.74
FY(1 - 3) -94.05
FY(2 - 4) -115.95
FZ(1 - 3) -32.84
FZ(2 - 4) 98.54

94.05 115.95

Figura 2 - 98. Diagrama de fuerza axial

Edgar Eduardo Muñoz, IC, MIC 125

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

115.95

115.95
FX(1 - 2) 57.74

94.05 FY(1 - 2) 94.05

FZ(1 - 2) = 32.84

FX(3 - 1) 7.74 FX(2 - 1) -57.74

FX(4 - 2) -57.74
FY(3 - 1) 94.05 FY(2 - 1) 115.95
FY(4 - 2) 115.95
FZ(2 - 1) -98.54
FZ(3 - 1) = 9.63
FZ(4 - 2) = 74.67
FX(1 - 3) -7.74
FX(2 - 4) 57.74
FY(1 - 3) -94.05
FY(2 - 4) -115.95
FZ(1 - 3) -32.84
FZ(2 - 4) 98.54

7.74

Figura 2 - 99. Diagrama de cortante

98.54
32.84
-32.84
98.54

FX(1 - 2) 57.74
FY(1 - 2) 94.05
FZ(1 - 2) = 32.84

FX(3 - 1) 7.74 FX(2 - 1) -57.74

FX(4 - 2) -57.74
FY(3 - 1) 94.05 FY(2 - 1) 115.95
FY(4 - 2) 115.95
FZ(2 - 1) -98.54
FZ(3 - 1) = 9.63
FZ(4 - 2) = 74.67
FX(1 - 3) -7.74
FX(2 - 4) 57.74
FY(1 - 3) -94.05
FY(2 - 4) -115.95
FZ(1 - 3) -32.84
FZ(2 - 4) 98.54

74.67

9.63

Figura 2 - 100. Diagrama de momento

Edgar Eduardo Muñoz, IC, MIC 126

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Ahora se determina para este mismo pórtico el efecto cuando se presente un

asentamiento en el nudo 4:

 F1x  0  U1x  ?  0 
 F1y  0  U  ?  0 
   1y   
 F1z  0   1z  ?   0 
     
 F2x  0  U1x  ?  0 
 F2y  0  U jy  ?  0 
     
 F2z  0   K nn K na    jz  ?   0 
 F3x  ?      
   K an K aa   0   0 
 F3y  ?   0   0 
     
 F3z  ?   0   0 
 F4x  ?   0   0 
     
 F4y  ?   0    0.01
 F4z  ?   0   0 
     

Figura 2 - 101. Sistema de solución

Se determina las deformaciones:

 n   Knn 1Fn  Kna  a 

 0.00219
 0 
 
 0.00146 
 n    
 0.00219
 0.00146 
 
 0.01 
Figura 2 - 102. Determinación de las deformaciones

Edgar Eduardo Muñoz, IC, MIC 127

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

[Link].1.3 Ejemplo: Resolver el pórtico plano con carga solamente vertical por giro -
deflexión

10 kN/m

1 2

3,0 m

3 4
6,0 m

Figura 2 - 103. Pórtico plano con carga vertical

Se plantean las condiciones de equilibrio que se deben cumplir en los nudos 1 y 2:

𝑀13 + 𝑀12 = 0 (2 - 95)

𝑀21 + 𝑀24 = 0 (2 - 96)

También las condiciones de deformaciones que se conocen para los nudos 3 y 4,

teniendo en cuenta que tienen un apoyo de tercer género.

𝜃3 = 0

𝜃4 = 0

Edgar Eduardo Muñoz, IC, MIC 128

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Se construyen las ecuaciones de cada elemento aplicando giro – deflexión y

considerando EI=1.0:

2𝐸𝐼 4
𝑀13 = (2 ∗ 𝜃1 + 𝜃3) = (𝜃1) (2 - 97)
3 3
2𝐸𝐼 2
𝑀31 = (2 ∗ 𝜃3 + 𝜃1) = (𝜃1) (2 - 98)
3 3
2𝐸𝐼 4
𝑀24 = (2 ∗ 𝜃2 + 𝜃4) = (𝜃2) (2 - 99)
3 3

2𝐸𝐼 2
𝑀42 = (2 ∗ 𝜃4 + 𝜃2) = (𝜃2) (2 - 100)
3 3

2𝐸𝐼 2 1
𝑀12 = (2 ∗ 𝜃1 + 𝜃2) + 30 = 𝜃1 + 𝜃2 + 30 (2 - 101)
6 3 3

2𝐸𝐼 2 1
𝑀21 = (2 ∗ 𝜃2 + 𝜃1) − 30 = 𝜃2 + 𝜃1 − 30 (2 - 102)
6 3 3

Reemplazando las ecuaciones (2 - 97), (2 - 101) en (2 - 95)

2 1 4
𝜃1 + 𝜃2 + (𝜃1) + 30 = 0
3 3 3
1
2𝜃1 + 𝜃2 = −30
3

Aplicamos las ecuaciones (2 - 99) y (2 - 102) en (2 - 96)

2 1 4
𝜃2 + 𝜃1 + (𝜃2) − 30 = 0
3 3 3

Edgar Eduardo Muñoz, IC, MIC 129

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

1
𝜃1 + 2(𝜃2) = 30
3

Se resuelve el sistema:

Determinando las rotaciones de los nudos 1 y 2, se determinan los momentos de

reacción en los nudos 3 y 4, aplicando las ecuaciones (2 - 98) y (2 - 100).

2
𝑀31 = (−18) = −12𝑘𝑁. 𝑚
3
2
𝑀42 = (18) = 12𝑘𝑁. 𝑚
3

Conociendo los momentos se determinan las demás reacciones haciendo equilibrio en

Z:

10*6-F4z-F3z=0

Sumatoria de momento igual a cero en el nudo 3:

F4z*6+12-12-10*6*3=0
F4z=30 KN

Edgar Eduardo Muñoz, IC, MIC 130

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

F3z=30 KN

Resultados que coinciden al resolver este mismo pórtico en SAP-2000

Figura 2 - 104. Resultados en SAP2000

[Link].1.4 Ejemplo: Resolver el pilón del puente por una carga vertical variable

Edgar Eduardo Muñoz, IC, MIC 131

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Coordenadas Restricciones Fuerzas en los nudos

X(m) Y(m) Ry Rz Rz Fx Fy Fz
1 0 106 0 0 0 0 -1000 0
2 0 68 0 0 0 0 0 0
3 -8 37 0 0 0 0 0 0
4 8 37 0 0 0 0 0 0
5 0 0 1 1 1 0 0 0

Incidencias de elementos y propiedades

Elemento Ni Nf A(m2) I(m4)

1 2 1 22 75,3334
2 3 2 9,6 28,8
3 4 2 9,6 28,8
4 3 4 1,8 0,0135
5 5 3 9,6 28,8
6 5 4 9,6 28,8

Fx1=0 354834,32 0,00 6741852,03 -354834,32 0,00 6741852,03 0,00 0,00 0,00 0,00 0,00 0,00 Ux1
Fy1=-1000 0,00 12469432,11 0,00 0,00 -12469432,11 0,00 0,00 0,00 0,00 0,00 0,00 0,00 Uy1
Fz1=0 6741852,03 0,00 170793584,83 -6741852,03 0,00 85396792,41 0,00 0,00 0,00 0,00 0,00 0,00 Uz1
Fx2=0 -354834,32 0,00 -6741852,03 1586660,81 0,00 289797,56 -615913,25 -1507707,63
3515824,79-615913,25 1507707,63 3515824,79 Ux2
Fy2=0 0,00 -12469432,11 0,00 0,00 24607821,05 0,00 -1507707,63 -6069194,47
-907309,62 1507707,63 -6069194,47 907309,62 Uy2
Fz2=0 6741852,03 0,00 85396792,41 289797,56 0,00 325792312,32-3515824,79 907309,62
38749681,87-3515824,79 -907309,62 38749681,87 Uz2
Fx3=0 0,00 0,00 0,00 -615913,25 -1507707,63 -3515824,79 3413984,19 407812,03
-977286,22 -2423037,38 0,00 0,00 Ux3
Fy3=0 0,00 0,00 0,00 -1507707,63 -6069194,47 907309,62 407812,03 11294281,77
1462997,620,00 -851,85 6814,79 Uy3
Fz3=0 0,00 0,00 0,00 3515824,79 -907309,62 38749681,87 -977286,22 1462997,62
143116663,49
0,00 -6814,79 36345,56 Uz3
Fx4=0 0,00 0,00 0,00 -615913,25 1507707,63 -3515824,79 -2423037,38 0,00 0,00 3413984,19 -407812,03 -977286,22 Ux4
Fy4=0 0,00 0,00 0,00 1507707,63 -6069194,47 -907309,62 0,00 -851,85-6814,79 -407812,03 11294281,77 -1462997,62 Uy4
Fz4=0 0,00 0,00 0,00 3515824,79 907309,62 38749681,87 0,00 6814,7936345,56 -977286,22 -1462997,62 143116663,49Uz4

Ux1 2,22875E-16
Uy1 -0,00027849
Uz1 -3,01571E-18
Ux2 1,17514E-16
Uy2 -0,000198294
Uz2 -2,29682E-18
Ux3 -4,39716E-05
Uy3 -0,000104914
Uz3 -4,90024E-07
Ux4 4,39716E-05
Uy4 -0,000104914
Uz4 0

Edgar Eduardo Muñoz, IC, MIC 132

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Fx5 0,00 0,00 0,00 0,00 0,00 0,00 -375033,57 1099895,60

-2538538,58-375033,57 -1099895,60 -2538538,58 Ux1=0,0000000
Fy5 0,00 0,00 0,00 0,00 0,00 0,00 1099895,60 -5224235,45
-548873,21 -1099895,60 -5224235,45 548873,21 Uy1=-0,0002785
Fz5 0,00 0,00 0,00 0,00 0,00 0,00 2538538,58 548873,21
32772304,312538538,58 -548873,21 32772304,31 Uz1=0,0000000
Ux2=0,0000000
Uy2=-0,0001983
Fx5 5,30687E-14 Uz2=0,0000000
Fy5 1.000 Ux3=-0,0000440
Fz5 0 Uy3=-0,0001049
Uz3=-0,0000005
Ux4=0,0000440
Uy4=-0,0001049
Uz4=0,0000005

2.1.5 Parrilla o entramado

[Link] Definición y antecedentes

Son estructuras que tienen cargas perpendiculares al plano. Los nudos se suponen
rígidos, por lo tanto, sus acciones principales son torsión, flexión y corte.

Edgar Eduardo Muñoz, IC, MIC 133

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 105. Modelo parrilla

Este modelo tiene aplicación en tableros de puentes como se observa en la figura

Figura 2 - 106. Tablero de puentes donde se aplica el modelo parrilla

También en entrepisos de edificios:

Figura 2 - 107. Tablero de edificios donde se aplica hacer el análisis estructural

mediante un modelo parrilla

Edgar Eduardo Muñoz, IC, MIC 134

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

[Link] Métodos de análisis

[Link].1 Método matricial

[Link].1.1 Fundamentos

Para su deducción por análisis matricial se debe recordar el concepto de mecánica de

solidos sobre torsión pura. La ecuación básica es la siguiente:

M xL
x  (2 - 103)
GJ

Donde:

Ox : Giro relativo en radianes

Mx : Momento torsor aplicado
L : Longitud del elemento
G : Módulo de cortante que se define como:

E
G (2 - 104)
2(1  U )

Donde:

E : Módulo de elasticidad
U : Relación de poissons

La constante de torsión ( J) de una sección cuando es circular maciza o hueca, es el

momento polar de inercia. Cuando es una secciones rectangulares es (donde b>t)

Edgar Eduardo Muñoz, IC, MIC 135

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

J  Cbt 3 (2 - 105)
t  1 t 
4
1
C   0.21 1    
3 b  12  b  

Por lo tanto la matriz de rigidez de elemento con solamente torsión es:

 GJ GJ 

 M xi   L L   xi 
 M    GJ GJ   xj 
 xj   
 L L 

Figura 2 - 108. Matriz de rigidez local de un elemento prismático a torsión pura

Basado en lo anterior se muestra la matriz de rigidez de un elemento localizado en

sentido X:

 JG JG 
 L 0 0  0 0 
L
 4 EI 6 EI 2 EI 6 EI    
 M xi   0  0  xi
M   L L2 L L2   
 yi   12EI  yi 
 3  U 
6 EI 12EI 6 EI
 Fzi   0  2 0  2
  L L3 L L   zi  (2 - 106)
 
 M xj   JG 0 0
JG
0 0   xj 
 M yj   L L   
   2 EI 6 EI 4 EI 6 EI   yj 
 Fzj   0  0
L2   zj 
U
L L2 L
 6 EI 12EI 6 EI 12EI 
 0  3 0 
 L2 L L2 L3 

Edgar Eduardo Muñoz, IC, MIC 136

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Posteriormente se debe determinar la matriz de rigidez para un elemento

arbitrariamente orientado aplicando la ecuación (2 - 24) que se basa en la matriz de
transformación (T) y la matriz de rigidez de elemento en sentido X antes explicada:

 JG JG 
 L 0 0  0 0 
L
 6 EI   Cos
 Cos Sin 4 EI 6 EI 2 EI Sin
T
0 0 0 0  0  2 0  0 0 0 0
 Sin Cos 0  L L L L2   0
 0 0 0
 0 6 EI 12EI 6 EI 12EI   Sin Cos 0 0 0
 0 0  2 0  2  3  0 0
0 1 0 0  L  0 1 0 0
K GLOBAL      JG
L L3 L
 
 0 0 0 Cos Sin 0
 0 0
JG
0 0  0 0 0 Cos Sin 0 (2 - 107)
 0 0 0  Sin Cos 0  L L  0 0 0  Sin Cos 0
   0 2 EI 6 EI 4 EI 6 EI   
 0 1  0  0 1
L 
0 0 0 0  0 0 0 0
L L2 L 2

 6 EI 12EI 6 EI 12EI 
 0  3 0 
 L2 L L2 L 
3

[Link].1.2 Ejemplo 1: Parrilla que debe soportar un motor

Se presenta a continuación un ejemplo de un entramado o parrilla. Consiste en

Analizar un entrepiso en voladizo que debe soportar un motor que produce una carga
puntual incluyendo impacto de 100 KN en dos nudos como se observa en la Figura 2
- 109. Las características de los materiales son:

E=17872000 KN/m2
G = 7446785 KN/m2
U=0.2

Edgar Eduardo Muñoz, IC, MIC 137

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

100 KN
Y 4.0m
2 100 KN X
4
1

2.0m
3

Figura 2 - 109. Modelo estructural parrilla

Los datos
Numero necesarios para el análisis
de nudos 4.0 matricial son:
Numero de elementos 3.0
Módulo de elasticidad 17872000.0
Módulo de cortante 7446685.0

Nudos Coordenadas Restricciones Fuerzas en los nudos

X(m) y(m) Rx Ry Rz Mx My Fz
1.0 2.0 0.0 0.0 0.0 0.0 0.0 0.0 -100.0
2.0 2.0 4.0 0.0 0.0 0.0 0.0 0.0 -100.0
3.0 0.0 0.0 1.0 1.0 1.0 0.0 0.0 0.0
4.0 0.0 4.0 1.0 1.0 1.0 0.0 0.0 0.0

Elemento Ni Nf A(m2) I(m4) j(m4)

1.0 3.0 1.0 0.12 0.001600 0.0019
2.0 4.0 2.0 0.12 0.001600 0.0019
3.0 1.0 2.0 0.12 0.001600 0.0019

Edgar Eduardo Muñoz, IC, MIC 138

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

La matriz de rigidez de cada elemento es:

La matriz de la estructura basada en un proceso de ensamblaje es:

Mx1 My1 Fz1 Mx2 My2 Fz2 Mx3 My3 Fz3 Mx4 My4 Fz4
Mx1 35818.5 0.0 10723.2 14297.6 0.0 -10723.2 -7223.3 0.0 0.0 0.0 0.0 0.0
My1 0.0 60809.5 42892.8 0.0 -3619.1 0.0 0.0 42892.8 -42892.8 0.0 0.0 0.0
Fz1 10723.2 42892.8 48254.4 10723.2 0.0 -5361.6 0.0 42892.8 -42892.8 0.0 0.0 0.0
Mx2 14297.6 0.0 10723.2 35833.4 0.0 -10723.2 0.0 0.0 0.0 -7238.2 0.0 0.0
My2 0.0 -3619.1 0.0 0.0 60809.5 42892.8 0.0 0.0 0.0 0.0 42892.8 -42892.8
Fz2 -10723.2 0.0 -5361.6 -10723.2 42892.8 48254.4 0.0 0.0 0.0 0.0 42892.8 -42892.8
Mx3 -7223.3 0.0 0.0 0.0 0.0 0.0 7223.3 0.0 0.0 0.0 0.0 0.0
My3 0.0 28595.2 42892.8 0.0 0.0 0.0 0.0 57190.4 -42892.8 0.0 0.0 0.0
Fz3 0.0 -42892.8 -42892.8 0.0 0.0 0.0 0.0 -42892.8 42892.8 0.0 0.0 0.0
Mx4 0.0 0.0 0.0 -7238.2 0.0 0.0 0.0 0.0 0.0 7238.2 0.0 0.0
My4 0.0 0.0 0.0 0.0 28595.2 42892.8 0.0 0.0 0.0 0.0 57190.4 -42892.8
Fz4 0.0 0.0 0.0 0.0 -42892.8 -42892.8 0.0 0.0 0.0 0.0 -42892.8 42892.8

La sub-matriz “ KNN” es:

Mx1 My1 Fz1 Mx2 My2 Fz2

Mx1 35818.5 0 10723.2 14297.6 0 -10723.2
My1 0 60809.5 42892.8 0 -3619.1 0
Fz1 10723.2 42892.8 48254.4 10723.2 0 -5361.6
Mx2 14297.6 0 10723.2 35833.4 0 -10723.2
My2 0 -3619.1 0 0 60809.5 42892.8
Fz2 -10723.2 0 -5361.6 -10723.2 42892.8 48254.4

Edgar Eduardo Muñoz, IC, MIC 139

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

La sub-matriz “ KAN” es:

Mx3 My3 Fz3 Mx4 My4 Fz4

Mx3 -7223.28 0 0 0 0 0
My3 0 28595.2 42892.8 0 0 0
Fz3 0 -42892.8 -42892.8 0 0 0
Mx4 0 0 0 -7238.18 0 0
My4 0 0 0 0 28595.2 42892.8
Fz4 0 0 0 0 -42892.8 -42892.8

Se determina las deformaciones de cada nudo aplicando la ecuación (2 - 11) es:

Rot X 1 -1.2044E-12
Rot Y 1 0.00699418
des Z 1 -0.00932557
Rot X 2 -6.0832E-13
Rot Y 2 0.00699418
des Z 2 -0.00932557

Las reacciones en los nudos aplicando la ecuación (2 - 12) es:

Mx 3 0.00
My 3 -200.00
Fz 3 100.00
Mx 4 0.00
My 4 -200.00
Fz 4 100.00

Se revisa equilibrio, lo cual se cumple al hacer sumatorias de fuerzas igual a cero en

sentido Z y momentos iguales a cero con respecto a cualquier nudo. Se determina
posteriormente las fuerzas internas de cada elemento con las cuales se puede hacer su
diagrama de cortante y momento;

Edgar Eduardo Muñoz, IC, MIC 140

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Elemento 31
MX(3 - 1) 7223.3 0.0 0.0 -7223.3 0.0 0.0 0.000000 0.0 0.0
MY(3 - 1) 0.0 57190.4 -42892.8 0.0 28595.2 42892.8 0.000000 0.0 -200.0
FZ(3 - 1) 0.0 -42892.8 42892.8 0.0 -42892.8 -42892.8 0.000000 0.0 100.0
MX(1 - 3) -7223.3 0.0 0.0 7223.3 0.0 0.0 0.000000 0.0 0.0
MY(1 - 3) 0.0 42892.8 -42892.8 0.0 57190.4 42892.8 0.006994 0.0 0.0
FZ(1 - 3) 0.0 42892.8 -42892.8 0.0 42892.8 42892.8 -0.009326 0.0 -100.0
Elemento 42
MX(4 - 2) 7238.2 0.0 0.0 -7238.2 0.0 0.0 0.000000 0.0 0.0
MY(4 - 2) 0.0 57190.4 -42892.8 0.0 28595.2 42892.8 0.000000 0.0 -200.0
FZ(4 - 2) 0.0 -42892.8 42892.8 0.0 -42892.8 -42892.8 0.000000 0.0 100.0
MX(2 - 4) -7238.2 0.0 0.0 7238.2 0.0 0.0 0.000000 0.0 0.0
MY(2 - 4) 0.0 42892.8 -42892.8 0.0 57190.4 42892.8 0.006994 0.0 0.0
FZ(2 - 4) 0.0 42892.8 -42892.8 0.0 42892.8 42892.8 -0.009326 0.0 -100.0
Elemento 12
MX(1 - 2) 28595.2 0.0 10723.2 14297.6 0.0 -10723.2 0.0 0.0 0.0
MY(1 - 2) 0.0 3619.1 0.0 0.0 -3619.1 0.0 0.0 0.0 0.0
FZ(1 - 2) 10723.2 0.0 5361.6 10723.2 0.0 -5361.6 0.0 0.0 0.0
MX(2 - 1) 14297.6 0.0 10723.2 28595.2 0.0 -10723.2 0.0 0.0 0.0
MY(2 - 1) 0.0 -3619.1 0.0 0.0 3619.1 0.0 0.0 0.0 0.0
FZ(2 - 1) -10723.2 0.0 -5361.6 -10723.2 0.0 5361.6 0.0 0.0 0.0

El diagrama de cortante es:

Figura 2 - 110. Diagrama de cortante

El diagrama de momento es:

Edgar Eduardo Muñoz, IC, MIC 141

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 111. Diagrama de momento

[Link].1.2 Ejemplo 2: Un balcón

Consiste en analizar un balcón de un edificio que debe soportar una carga distribuida,
una carga puntual en el nudo 1 y otra carga puntual en el centro del elemento 1-2
como se observa en la Figura 2 - 109.(complementado del ejercicio 11.26 de la
referencia [6]).

20 kN/m
Y

3
40 kN

1 50 kN

1,2m 3m
2
1,2m

Edgar Eduardo Muñoz, IC, MIC 142

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 112. Sistema estático y cargas de Balcón,

Las características de los materiales y las propiedades geométricas de los dos

elementos son las siguientes:

Los datos necesarios para el análisis matricial son:

Nudos Coordenadas Restricciones Fuerzas en los nudos

X(m) y(m) Rx Ry Rz Fx Fy Fz
1,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 -40,0
2,0 2,4 0,0 1,0 1,0 1,0 0,0 0,0 0,0
3,0 0,0 3,0 1,0 1,0 1,0 0,0 0,0 0,0

Elemento Ni Nf I(m4) J(m4)

1,0 1 2 0,00160 0,00194
2,0 1 3 0,00160 0,00194

Se construye la matriz de rigidez de cada uno de los elementos:

Edgar Eduardo Muñoz, IC, MIC 143

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Se determinan los vectores de empotramiento

Edgar Eduardo Muñoz, IC, MIC 144

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Se determina mediante ensamblaje la submatriz Knn (tener en cuenta que los índices
de las filas y columnas de esta matrices se numeran de 0 a 5, y no de 1 a 6):

Se determina las deformaciones de cada nudo aplicando la ecuación (2 - 11) es:

Edgar Eduardo Muñoz, IC, MIC 145

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Se determina mediante ensamblaje la submatriz Kan:

Se hallan las reacciones en los nudos 2 y 3 aplicando la ecuación (2 - 12) es:

Se revisa equilibrio, donde se demuestra que se cumple con esta condición en las tres
direcciones. Se hace la revisión de fuerzas verticales en Z, donde P1=40 kN( carga en
el nudo 1), W= 20 kN/m y P2=50 kN:

Momento alrededor de x en el nudo 1:

Edgar Eduardo Muñoz, IC, MIC 146

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Momento alrededor de y en el nudo 1:

Ahora se determinan las fuerzas internas, con lo cual se hacen los diagramas de
cortante y momento, como se observa a continuación:

Elemento 12
MX(1 - 2) 6415,2 0,0 0,0 -6415,2 0,0 0,0 0,002301 0,0 14,8
MY(1 - 2) 0,0 50666,7 -31666,7 0,0 25333,3 31666,7 -0,003176 -15,0 16,3
FZ(1 - 2) 0,0 -31666,7 26388,9 0,0 -31666,7 -26388,9 -0,006071 25,0 -34,6
MX(2 - 1) -6415,2 0,0 0,0 6415,2 0,0 0,0 0,000000 0,0 -14,8
MY(2 - 1) 0,0 25333,3 -31666,7 0,0 50666,7 31666,7 0,000000 15,0 126,8
FZ(2 - 1) 0,0 31666,7 -26388,9 0,0 31666,7 26388,9 0,000000 25,0 84,6
Elemento 13
MX(1 - 3) 40533,3 0,0 20266,7 20266,7 0,0 -20266,7 0,002301 15,0 -14,8
MY(1 - 3) 0,0 5132,2 0,0 0,0 -5132,2 0,0 -0,003176 0,0 -16,3
FZ(1 - 3) 20266,7 0,0 13511,1 20266,7 0,0 -13511,1 -0,006071 30,0 -5,4
MX(3 - 1) 20266,7 0,0 20266,7 40533,3 0,0 -20266,7 0,000000 -15,0 -91,4
MY(3 - 1) 0,0 -5132,2 0,0 0,0 5132,2 0,0 0,000000 0,0 16,3
FZ(3 - 1) -20266,7 0,0 -13511,1 -20266,7 0,0 13511,1 0,000000 30,0 65,4

Figura 2 - 113. Diagrama de cortante

Edgar Eduardo Muñoz, IC, MIC 147

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Figura 2 - 114. Diagrama de momento

2.1.6 Pórtico espacial

[Link] Definición y antecedentes

[Link] Métodos de análisis

[Link].1 Método matricial

[Link].1.1 Fundamentos

Se parte por definir la matriz de rigidez local en sentido X de un elemento prismático

en el espacio como se muestra en el siguiente esquema:

Miz
Mjz

Fjz
Fix Fiz Fjx
i j
x
Mix Mjx
Fiy Fjy
z Miy Mjy

Figura 2 - 115. Elemento prismático de pórtico especial en el eje x

Edgar Eduardo Muñoz, IC, MIC 148

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 AE AE 
 L 0 0 0 0 0  0 0 0 0 0 
L
 12EI z 6 EI z 12EI 6 EI z 
 0 0 0 0 0  3 z 0 0 0 
 L3 L2 L L2 
 Fxi   0 0 U xi 
12EI y 6 EI 12EI y 6 EI
0 0  2y 0 0 0  0  2y
F   L3 L L3 L  
 yi   0 GJ

GJ  U yi
0  
 Fzi   0 0 0 0 0 0 0 0
U
   zi 
L L

0   xi 
6 EI 4 EI y 6 EI y 2 EI y
 M xi   0 0  2y 0 0 0 0 0
 M yi   L L L L  

   0 6 EI z 4 EI z

6 EI z 2 EI z  yi  (2 - 108)
 M zi   
L   zi 
0 0 0 0 0 0 0
L2 L L2
 F    AE AE U 
 xj   0 0 0 0 0 0 0 0 0 0  xj 
 Fyj   L L U 
F   0 
12EI z
0 0 0 
6 EI z
0
12EI z
0 0 0
6 EI z  yj 
 2  U zj
 zj   L3 L2 L3 L  
 M xj   12EI y 6 EI y 12EI y 6 EI y   
   0 0  0 0 0 0 0 0  xj 
 M yj   L3 L2 L3 L2   yj 
 M zj   0 0 0 
GJ
0 0 0 0 0
GJ
0 0   zj 
   L L 
 6 EI 2 EI y 6 EI y 4 EI y 
 0 0  2y 0 0 0 0 0 0 
 L L L2 L 
 0 6 EI z
0 0 0
2 EI z
0 
6 EI z
0 0 0
4 EI z 
 L2 L L2 L 

Posteriormente se debe determinar la matriz de rigidez para un elemento

arbitrariamente orientado aplicando la ecuación (2 - 24), donde es necesario apoyarse
en la matriz del elemento antes mostrada y definir la matriz de transformación (T).
Esta matriz es de 12x12 y está compuesta por cuatro sub matrices iguales como se
expresa a continuación:

T1(3 x 3)
T1(3 x 3)
T 
T1(3 x 3) (2 - 109)

T1(3 x 3)

La submatriz (T1) es la siguiente:

Edgar Eduardo Muñoz, IC, MIC 149

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 x  CosxX CosxY CosxZ   X 

 y   CosyX CosyY CosyZ   Y 
   (2 - 110)
 z  CoszX CoszY CoszZ   Z 

Donde los términos de la primera fila corresponden a los cosenos directores:

X 2  X1
CosxX  C1 
L
Y2  Y1
CosxY  C2 
L
Z 2  Z1
CosxZ  C3 
L

Para determinar los términos de la tercera fila utilizamos un vector unitario en Y y

realizamos producto cruz

i j k
  
Z unitario  X x J  C1 C2 C3  C3i  C1k
0 1 0


Z  C3i C1k  C3i C1k
Z    
d C3  C1
2 2
C3  C1
2 2
d d

d  C32  C12

Edgar Eduardo Muñoz, IC, MIC 150

 Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Para determinar los términos de la segunda fila se hace producto un cruz entre los
vectores X y Z:

i j k
 
C C1
Y ZxX   3 0  C1C2i  dJ  C2C3k
d d
C1 C2 C3

Con ellos se determinaron todos los términos de dicha matriz que queda de la
siguiente forma:

 
 x  CosxX CosxY CosxZ   X   C1 C2 C3   X 
 y   CosyX  CC CC 
   CosyY CosyZ   Y    1 2 d  2 3   Y 
CoszY CoszZ   Z   C3
d d 
 z  CoszX C1   Z 
  0
 d d 

Con lo cual la matriz de transformación (T) queda como:

 C1 C2 C3 
 C1C2 C 2 C3 
 d  
 d d 
  C3 0
C1 
 d d 
 C1 C2 C3 
 C1C2 C 2 C3 
  d  
 d d 
 C C1 
 3 0
T    d d
C1 C2 C3


(2 - 111)
 CC CC 
  1 2 d  2 3 
 d d 
 C
 3 0
C1 
 d d 
 C1 C2 C3 
 C1C2 CC 
  d  2 3
 d d 
 C
 3 0
C1 
 d d 

La anterior matriz de transformación no es válida para elementos verticales para lo

cual es necesario utilizar la siguiente:

Edgar Eduardo Muñoz, IC, MIC 151

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

 0 C2 0 
 C 0 0 
 2 
 0 0 1 
 
 0 C2 0 
  C2 0 0 
 
T    
0 0 1
0 C2 0 
  (2 - 112)
  C2 0 0 
 0 0 1 
 
 0 C2 0
 
  C2 0 0
 0 0 1

[Link].1.2 Ejemplo: Análisis de una pila en construcción

El ejercicio consiste en hacer el análisis estructural de un puente en voladizos

sucesivos en una de sus etapas de construcción, que corresponde cuando solamente
hay dos dovelas construidas como se muestra a continuación:

Edgar Eduardo Muñoz, IC, MIC 152

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

6m 6m

2 3
1

Dovelas

8m
Pila

100 kN 100 kN
2

1 3

8m

Z
x
4

Figura 2 - 116. Modelo estructural de pila-dovela en construcción

Edgar Eduardo Muñoz, IC, MIC 153

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Numero de nudos 4.00

Numero de elementos 3.00
Módulo de elasticidad 2.31E+07 KN/m2
Módulo de cortante 9613629.648 KN/m2

Nudos Coordenadas(m ) Restricciones

X Y Z Dx Dy Dz Rx Ry Rz
1 0 8 0 0 0 0 0 0 0
2 6 8 0 0 0 0 0 0 0
3 12 8 0 0 0 0 0 0 0
4 6 0 0 1 1 1 1 1 1

Elem ento Ni Nf A(m 2) Iz(m 4) Iy(m 4) J(m 4)

1 1 2 0.7 5.83E-02 0.028583 2.11E-03
2 2 3 0.7 5.83E-02 0.028583 2.11E-03
3 4 2 0.7 5.83E-02 0.028583 2.11E-03

Tabla 2 - 17. Datos de entrada

Se determinan la matriz de transformación de cada uno de los cuatro elementos:

Ux1 Uy1 Uz1 Ox1 Oy1 Oz1 Ux2 Uy2 Uz2 Ox2 Oy2 Oz2
1 0 0
0 1 0
0 0 1
1 0 0
0 1 0
0 0 1
1 0 0
0 1 0
0 0 1
1 0 0
0 1 0
0 0 1

Ux2 Uy2 Uz2 Ox2 Oy2 Oz2 Ux3 Uy3 Uz3 Ox3 Oy3 Oz3
1 0 0
0 1 0
0 0 1
1 0 0
0 1 0
0 0 1
1 0 0
0 1 0
0 0 1
1 0 0
0 1 0
0 0 1

Edgar Eduardo Muñoz, IC, MIC 154

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Ux4 Uy4 Uz4 Ox4 Oy4 Oz4 Ux2 Uy2 Uz2 Ox2 Oy2 Oz2
1
-1
1
1
-1
1
1
-1
1
1
-1
1

Después la matriz de rigidez global de cada elemento:

Ux1 Uy1 Uz1 Ox1 Oy1 Oz1 Ux2 Uy2 Uz2 Ox2 Oy2 Oz2
Fx(1 - 2)= 2691816 0 0 0 0 0 -2691816 0 0 0 0 0
Fy(1 - 2)= 0 74768 0 0 0 224305 0 -74768 0 0 0 224305
Fz(1 - 2)= 0 0 36638 0 -109915 0 0 0 -36638 0 -109915 0
Mx(1 - 2)= 0 0 0 3381 0 0 0 0 0 -3381 0 0
My(1 - 2)= 0 0 -109915 0 439658 0 0 0 109915 0 219829 0
Mz(1 - 2)= 0 224305 0 0 0 897221 0 -224305 0 0 0 448610
Fx(2 - 1)= -2691816 0 0 0 0 0 2691816 0 0 0 0 0
Fy(2 - 1)= 0 -74768 0 0 0 -224305 0 74768 0 0 0 -224305
Fz(2 - 1)= 0 0 -36638 0 109915 0 0 0 36638 0 109915 0
Mx(2 - 1)= 0 0 0 -3381 0 0 0 0 0 3381 0 0
My(2 - 1)= 0 0 -109915 0 219829 0 0 0 109915 0 439658 0
Mz(2 - 1)= 0 224305 0 0 0 448610 0 -224305 0 0 0 897221

Ux2 Uy2 Uz2 Ox2 Oy2 Oz2 Ux3 Uy3 Uz3 Ox3 Oy3 Oz3
Fx(2 - 3)= 2691816 0 0 0 0 0 -2691816 0 0 0 0 0
Fy(2 - 3)= 0 74768 0 0 0 224305 0 -74768 0 0 0 224305
Fz(2 - 3)= 0 0 36638 0 -109915 0 0 0 -36638 0 -109915 0
Mx(2 - 3)= 0 0 0 3381 0 0 0 0 0 -3381 0 0
My(2 - 3)= 0 0 -109915 0 439658 0 0 0 109915 0 219829 0
Mz(2 - 3)= 0 224305 0 0 0 897221 0 -224305 0 0 0 448610
Fx(3 - 2)= -2691816 0 0 0 0 0 2691816 0 0 0 0 0
Fy(3 - 2)= 0 -74768 0 0 0 -224305 0 74768 0 0 0 -224305
Fz(3 - 2)= 0 0 -36638 0 109915 0 0 0 36638 0 109915 0
Mx(3 - 2)= 0 0 0 -3381 0 0 0 0 0 3381 0 0
My(3 - 2)= 0 0 -109915 0 219829 0 0 0 109915 0 439658 0
Mz(3 - 2)= 0 224305 0 0 0 448610 0 -224305 0 0 0 897221

Ux4 Uy4 Uz4 Ox4 Oy4 Oz4 Ux2 Uy2 Uz2 Ox2 Oy2 Oz2
Fx(4 - 2)= 31543 0 0 0 0 -126172 -31543 0 0 0 0 -126172
Fy(4 - 2)= 0 2018862 0 0 0 0 0 -2018862 0 0 0 0
Fz(4 - 2)= 0 0 15457 61827 0 0 0 0 -15457 61827 0 0
Mx(4 - 2)= 0 0 61827 329744 0 0 0 0 -61827 164872 0 0
My(4 - 2)= 0 0 0 0 2536 0 0 0 0 0 -2536 0
Mz(4 - 2)= -126172 0 0 0 0 672916 126172 0 0 0 0 336458
Fx(2 - 4)= -31543 0 0 0 0 126172 31543 0 0 0 0 126172
Fy(2 - 4)= 0 -2018862 0 0 0 0 0 2018862 0 0 0 0
Fz(2 - 4)= 0 0 -15457 -61827 0 0 0 0 15457 -61827 0 0
Mx(2 - 4)= 0 0 61827 164872 0 0 0 0 -61827 329744 0 0
My(2 - 4)= 0 0 0 0 -2536 0 0 0 0 0 2536 0
Mz(2 - 4)= -126172 0 0 0 0 336458 126172 0 0 0 0 672916

Posteriormente se plantea la solución del sistema por desplazamientos:

 Fx1  0  U x1  ? 
 F  100 U  ? 
 y 1   y1 
 Fz1  0   U z1  ? 
   
 M x1  0    x1  ? 
 M y1  0    y1  ? 
   
 M z1  0    z1  ? 
 F 0  U  ?
 x2   x2 
 y2F  0  U y 2  ?
 F 0  U  ? 
 z2   z2 
 M x2  0   x2  ? 
   
 M y2  0   y 2  ? 
 z2
M  0   Knn Kna  z 2  ? 
Edgar Eduardo Muñoz, IC, F  0
MIC   
 x 3   Kan Kaa U x 3  ?  155
 Fy 3  100 U y 3  ?
   
 Fz 3  0  U z 3  ? 
M  0    ? 
 x3   x3 
M  0 
Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Se procede al proceso de ensamblaje y se determina la sub matriz de la estructura

KNN y el sistema de solución:

Ux1 Uy1 Uz1 Ox1 Oy1 Oz1 Ux2 Uy2 Uz2 Ox2 Oy2 Oz2 Ux3 Uy3 Uz3 Ox3 Oy3 Oz3
2691816.3 0 0 0 0 0 -2691816 0 0 0 0 0 0 0 0 0 0 0
0 74768.4 0 0 0 224305.2 0 -74768.4 0 0 0 224305.21 0 0 0 0 0 0
0 0 36638.18 0 -109915 0 0 0 -36638.2 0 -109915 0 0 0 0 0 0 0
0 0 0 3380.8 0 0 0 0 0 -3380.79 0 0 0 0 0 0 0 0
0 0 -109915 0 439658.2 0 0 0 109914.6 0 219829.1 0 0 0 0 0 0 0
0 224305 0 0 0 897220.8 0 -224305 0 0 0 448610.41 0 0 0 0 0 0
-2691816.3 0 0 0 0 0 5415175.5 0 0 0 0 126171.68 -2691816 0 0 0 0 0
0 -74768.4 0 0 0 -224305 0 2168399 0 0 0 0 0 -74768.4 0 0 0 224305.2
0 0 -36638.2 0 109914.6 0 0 0 88733.1 -61826.9 0 0 0 0 -36638.2 0 -109915 0
0 0 0 -3381 0 0 0 0 -61826.9 336505 0 0 0 0 0 -3381 0 0
0 0 -109915 0 219829.1 0 0 0 0 0 881852 0 0 0 109915 0 219829.1 0
0 224305 0 0 0 448610.4 126171.68 0 0 0 0 2467357.28 0 -224305 0 0 0 448610.4
0 0 0 0 0 0 -2691816 0 0 0 0 0 2691816.3 0 0 0 0 0
0 0 0 0 0 0 0 -74768.4 0 0 0 -224305.21 0 74768.4 0 0 0 -224305
0 0 0 0 0 0 0 0 -36638.2 0 109914.6 0 0 0 36638.2 0 109914.6 0
0 0 0 0 0 0 0 0 0 -3380.79 0 0 0 0 0 3380.8 0 0
0 0 0 0 0 0 0 0 -109915 0 219829.1 0 0 0 109915 0 439658.2 0
0 0 0 0 0 0 0 224305.2 0 0 0 448610.41 0 -224305 0 0 0 897220.8

U x1  ?   Fx1  0   U x1  0 
U  ?   F  100  U  0.005448 
 y 1   y 1   y1 
 U z1  ?   Fz1  0   U z1  0 
     

 x1  ?   M x1  0   θ x1  0 
  y1  ?   M y1  0   θ y1  0 
     

 z1  ?   M z1  0 θ
  z1  -0.001337 
U  ?  F 0   U x2  0 
 x2   x2   
U y 2  ?  Fy 2  0   U y2  0.0000991
U  ?   F 0   U z2  0 
 z2   Knn1  z 2  
 x2  ?   M x2  0   θ x2  0 
     
 y 2  ?   M y2  0   θ y2  0 
 z 2  ?   M z 2  0   θ z2  0 
U  ?  F  0  U 0 
 x3   x3   x3 
U y 3  ?  Fy 3  100  U y3  0.005448 
     
U z 3  ?   Fz 3  0   U z3  0 
  ?  M  0   θ  0 
 x3   x3   x3 
Edgar Eduardo Muñoz, IC,  y 3  ? 
MIC  M y 3  0   θ y3  0  156
     

 z3  ?   z3
M  0 θ
  z3  -0.001337 
Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Posteriormente se determina las reacciones:

 U x1  0 
 U  0.005448 
 y1 
 U z1  0 
 
 θ x1  0 
 θ y1  0 
 
 θ z1  -0.001337 
 Fx 4   U x2  0   0 
F   
 U y2  0.0000991 200
 y 4     
 Fz 4  U z2  0  0 
   Kan

θ x2  0
 
  0 
M x 4    
M y 4   θ y2  0   0 
   θ    0 
 M z 4  z2 0
U 0 
 x3 
 U y3  0.005448 
 
 U z3  0 
θ 0 
 x3 
 θ y3  0 
 
 θ z3  -0.001337 

Edgar Eduardo Muñoz, IC, MIC 157

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

Finalmente se cumple con equilibrio.

3. REFERENCIAS

[1] Muñoz, E. E., Núñez, F., Otálora, C., Rodríguez, W. y Ruiz, D. “Ingeniería de
Puentes –Tomo 3“. Capítulo 3: Confiabilidad estructural de un puente en acero
apoyada en monitoreo e instrumentación. Grupo de investigación de “Estructuras y
Construcción“ del Departamento de Ingeniería Civil de la Pontificia Universidad
Javeriana. Editorial Javeriana, Agosto 2011.
[2 Muñoz, E. E., Núñez, F., Otálora, C., Rodríguez, . y Ruiz, D., “Confiabilidad
estructural de un puente en acero apoyada en monitoreo e instrumentación,” Revista
Ingeniería y Universidad, vol. 10, no. 1, pp. 31-54, enero-junio 2006.
[3] Galindo, Jorge. Cruzando el Cauca. Colección de autores Vallecaucanos, Premio
Jorge Isaacs. Cali: Secretaria de Cultura y Turismo del Valle del Cauca, 2002.
[4] Gennaro Tampone y Francesca Funis. Palladio's timber bridges. Proceedings of
the First International Congress on Construction History, Madrid, 20th-24th Enero,
2003

[5] Nelson, James y McCormac Nelson,” Análisis de estructuras – Método clásico y

matricial”, Alfaomega, Estados Unidos, 2006.

[6] Uribe, Escamilla Jairo, “Análisis de estructuras”, ogotá, Ediciones Uniandinas,

1995.

[7] Carlos Alberto Riveros Jerez, “Análisis Estructural. Teorema de Castigliano”,

Departamento de Ingeniería. Sanitaria y Ambiental. Facultad de Ingeniería. Obras
Civiles – Ingeniería Sanitaria UdeA.

[8] Hibbeler, R.C, “Análisis Estructural“, Prentice Hall, Tercera edición. 1997.

[9] “Reglamento Colombina de construcción sismo resistente – NSR-10”,

Asociación Colombiana de Ingeniería Sísmica, Bogotá, 2010.

Rochel Awad, Roberto (1993), “Análisis matricial de estructuras”, Universidad

[10]
EAFIT.

Edgar Eduardo Muñoz, IC, MIC 158

Ing. Edgar E. Muñoz Díaz, MIC Apuntes de clase de estructuras (en proceso y revisión)

[11] García Reyes, Luis Enrique (1997), “Dinámica estructural aplicada al diseño
sísmico”, ogotá, Uniandes.

Edgar Eduardo Muñoz, IC, MIC 159