Efectos de Segundo Orden 1

i
Universidad Austral de Chile

Facultad de Ciencias de la Ingeniera
Escuela de Ingeniera en Obras Civiles

ANLISIS COMPARATIVO DE MTODOS DE CLCULO DE
ESTRUCTURAS FORMADAS POR BARRAS,
CONSIDERANDO EFECTOS DE SEGUNDO ORDEN.

Para optar al ttulo de:
Ingeniero Civil en Obras Civiles

Profesor Patrocinante
Sr. Julio Lopetegui Torres
Ingeniero Civil, Dr. en Ingeniera

Profesores Informantes
Sr. Jos Soto Miranda
Ingeniero Civil, M. Sc. Eng. Civil

Sr. Adolfo Castro Bustamante
Ingeniero Civil, M. Sc. Eng. Civil

NGELA PATRICIA GONZLEZ AVILS
VALDIVIA CHILE
2009

Agradecimientos

A Dios por permitirme concluir esta etapa de mi vida
brindndome fortaleza y proteccin

A mis hermanos y padres Patricia y Hernn
por su amor y apoyo incondicional

A mi Abuela Sonia
Por su preocupacin y ayuda

A mis pastores Maritza y Henrique
por ceer siempre en m

A mis amigos
por su palabras de aliento y consejo
ii
ndice general

ndice general i
ndice temtico i
ndice de anexos iii
ndice de figuras iv
ndice de tablas vi
ndice de grficos vi
Resumen vii
Summary viii

ndice temtico

Captulo I: Introduccin
1.1. Planteamiento del problema 1
1.2. Objetivos 9
1.3. Metodologa 10

Captulo II: Mtodos de anlisis Elstico de segundo orden.
2.1. Conceptos generales 11
2.1.1. Inestabilidad de una barra 11
2.1.2. Pandeo elstico 13
2.1.3. Pandeo inelstico 15
2.1.4. Pandeo lateral-torsional 16
2.2. Mtodo de la Matriz Geomtrica 19
iii
Captulo III: Efectos de segundo orden considerados en las normas de diseo.

3.1. Recomendaciones del Instituto nacional del acero ICHA 24
3.2. Mtodo de anlisis de segundo orden establecido en el cdigo AISC 27
3.2.1. Consideraciones generales 27
3.2.2. Anlisis de segundo orden por un anlisis elstico de primer
orden amplificado. Mtodo de amplificacin de momentos 29
3.2.3. Mtodo de anlisis de segundo orden directo 32
3.3. Cdigo Europeo, EC3 36
3.3.1. Consideraciones generales 36
3.3.2. Anlisis de la estabilidad de un marco y factor de
amplificacin de momentos de primer orden 39
3.4. Resumen y comparaciones 42

Captulo IV: Programas computacionales
4.1. SAP2000 44
4.1.1. Introduccin 44
4.1.2. Forma cbica deformada 46
4.1.3. Fuerzas P- en los elemento prtico 47
4.1.4. Anlisis inicial P- 49
4.2. PROGES 50
4.2.1. Introduccin 50
4.2.2. Fundamentos tericos de PROGES 51

Captulo V: Estructura tipo marco plano analizada
5.1. Modelo de la estructura en Proges 57
iv
5.2. Modelo de la estructura en SAP2000 57

5.2.1. Casos de carga 57
5.2.2. Casos de anlisis 58
5.2.3. Combinaciones 60
5.3. Resultados 62

Captulo VI: Comentarios y conclusiones
6.1. Conclusiones generales 69
6.1.1 Respecto a los mtodos estudiados 69
6.1.2. Respecto a las normas estudiadas 71
6.2. Aportaciones ms significativas 72
6.3. Futuras lneas de investigacin 74

Captulo VII: Bibliografa 75

ndice de anexos

Anexo A: Funciones de estabilidad 78
Anexo B: Anlisis por el Mtodo de la matriz geomtrica 84
Anexo C: Archivos de entrada y salida de datos 95
Anexo D: Clculo de Factores B
2
y F
A
119

ndice de figuras

Figura 1.1.1. Efecto P- y P- 2
Figura 1.1.2. Diagrama de carga-deformacin 3
Figura 1.1.3. Tipos de anlisis 4
Figura 1.1.4. Interaccin entre el sistema estructural y los miembros que la
Componen 4
Figura 2.1.1.1. Estabilidad de una barra rgida 11
Figura 2.1.1.2. Curvas carga-rotacin de columnas 13
Figura 2.1.2.1. Columna Euler 13
Figura 2.1.4.1. Deformacin lateral-torsional de viga en flexin 17
Figura 2.1.4.2. Deformacin lateral-torsional de viga en flexin, detalle 17
Figura 2.2.1. Convencin de signos 20
Figura 3.2.1.1. Efecto P- y P- en vigas columna 28
Figura 3.2.2.1. Amplificacin de momento 29
Figura 3.2.2.2. Curvatura que ocasionan los momentos aplicados en los
extremos de los elementos 31
Figura 3.3.1.1. Imperfeccin de prtico traslacional 38
Figura 3.3.1.2. Imperfeccin en el elemento 38
Figura 3.3.1.3. Esfuerzos horizontales equivalentes 39
Figura 3.3.2.1. Prtico de edificacin con las vigas unidas a los pilares
en cada planta 41
Figura 4.1.1.1. Efecto P- en SAP2000 45
Figura 4.1.2.1. Geometra deformada 46
Figura 4.1.4.1. Anlisis no lineal esttico 49
vi
Figura 4.2.1.1. Diagrama de flujo de programa Proges 50

Figura 4.2.2.1. Ecuaciones de estabilidad para un segmento de columna 51
Figura 5.1. Marco analizado por los tres mtodos 56
Figura 5.1.1. Modelo de marco en Proges 57
Figura 5.2.1.1. Casos de carga en SAP2000 58
Figura 5.2.2.1. Caso de anlisis de PP en SAP2000 58
Figura 5.2.2.2. Caso de anlisis de SC en SAP2000 59
Figura 5.2.2.3. Caso de anlisis de Pdelta en SAP2000 59
Figura 5.2.2.4. Especificacin de efecto P-Delta 60
Figura 5.2.3.1. Datos de las propiedades del material 61
Figura 5.3.1. Fuerzas internas en cada barra 63
Figura A.1. Fuerzas alojadas en cada uno de los planos principales de
la barra 78
Figura A.2. Viga columna sujeta a momentos en sus extremos 79
Figura A.3. Viga columna sujeta a momentos en sus extremos 80
Figura A.4. Grfico de funciones de estabilidad 82
Figura A.5. Viga-columna sujeta a momentos y desplazamientos en
los extremos 83
Figura B.1. Modelo de marco y condiciones de carga 84
Figura B.2. Sistema de coordenadas globales y locales considerados
en mtodo 1 84
Figura B.3. Barra 1 del marco 85
Figura B.6. Determinacin de las fuerzas longitudinales 91
vii
Figura B.7. Geometra y combinaciones de carga 92

Figura B.8. Valores del determinante de K
total
usando la matriz geomtrica 93
Figura B.9. Formas de modos de pandeo correspondientes a
93

ndice de Tablas

Tabla 3.3.1.1. Imperfecciones de los elementos 39
Tabla 3.3.1.2. Factores de imperfecciones para las curvas de pandeo 39
Tabla 3.4.1. Tipos de anlisis de acuerdo al origen de las fuerzas y deformacin 42
Tabla 4.2.2.1. Valores de
para la matriz de rigidez 53

Tabla 5.3.1. Resultados de desplazamientos y fuerzas internas de los elementos 62
Tabla 5.3.2. Valores de deformacin vs carga aplicada en nudos 2 y 3 63
Tabla 5.3.3. Clculo del factor B
2
66
Tabla 5.3.4. Clculo del factor
cr
67
Tabla B.1 Valores de deformacin vs factor carga 94

ndice de grficos

Grfico 5.3.1. Curvas de carga vs deformacin 64
Grfico 5.3.2. Variacin porcentual de los desplazamientos obtenidos por los tres
Mtodos con respecto un anlisis lineal. 65
Grfico 5.3.3. Variacin porcentual de desplazamientos de MG vs PROGES 65
Grfico 5.3.4. Curva factores de amplificacin (B
2
y F
A
) vs carga 68
Grfico B.2. Curva factores de amplificacin (B
2
y F
A
) vs carga 94
viii
Resumen

Con el objeto de determinar de qu forma influyen las cargas axiales sobre una estructura
esbelta de acero, y dar a conocer una manera ms directa de obtener los esfuerzos considerando la
geometra deformada de dicha estructura, se muestran en la presente memoria los resultados
obtenidos a partir de un Anlisis Elstico de Segundo Orden y un Anlisis Elstico de Primer
Orden realizado a una estructura de acero tipo marco. Se presentan y comparan tres mtodos
mediante los cuales se realizan los clculos, considerando la deformacin de la estructura y los
elementos que la componen sobre los cuales estarn actuando las fuerzas, representadas por el
peso propio de los elementos y la sobrecarga.
Se complementa lo anterior con el estudio de distintos cdigos de diseo, en los cuales se
plantean los tipos de anlisis requeridos (sean estos lineales o no lineales) dependiendo
principalmente de las condiciones de arriostramiento de una edificacin, su altura y cargas
impuestas sobre esta.
A partir de la curva de carga vs deformacin obtenida a travs de cada mtodo se
concluye que para situaciones especificas de marcos no arriostrados, sometidos a cargas laterales
y axiales de mediana magnitud, es conveniente realizar un anlisis de segundo orden, debido a
que las deformaciones pueden hacer que los esfuerzos en los elementos aumenten de manera
significativa, afectando la estabilidad de la estructura, asimismo se evita comprobar la esbeltez de
cada uno de sus elementos, como ocurre en el caso de realizar un Anlisis Lineal, sin considerar
la deformacin de la estructura y sus elementos en una primera instancia.

ix
Summary

With the intention of determining how they influence the axial loads a slim steel structure,
and presenting one more a way more direct to obtain the efforts considering the deformed
geometry of this structure, the results obtained from an Elastic Analysis of Second Order and an
Elastic Analysis of First realised Order are in the present memory a steel type frame. Three
methods appear and compare by means of which the calculations are realised, considering the
deformation of the structure and the elements compose that it on which the forces will be acting,
represented by the own weight of the elements and the overload.
The previous thing with the study of different codes from design is complemented, in
which the types of required analyses consider mainly (they are these linear or nonlinear)
depending on the conditions of bracing of a construction, their height and loads imposed on this.
From the load curve versus deformation obtained through each method one concludes that
for specific situations of marks nonbraced, put under loads lateral axial Y of median magnitude,
he is advisable to realise an analysis of second order, because the deformations can cause that the
efforts in the elements increase of significant way, affecting the stability of the structure, also is
avoided to verify the slenderness of each of their elements, as it happens in the case of realising a
linear Analysis, without considering the deformation of the structure and its elements in one first
instance.

Capitulo I Introduccin
1

Captulo I
Introduccin

1.1. Planteamiento del problema

El empleo de aceros de alta resistencia y de otros materiales como el aluminio, as como
la utilizacin de nuevas formas constructivas, han hecho que las estructuras modernas, estn
generalmente formadas por elementos muy esbeltos, en los que los fenmenos de inestabilidad
adquieren una importancia fundamental, que hace aumentar la trascendencia del problema de
pandeo de columnas, que puede considerarse como la base para el estudio de todos los problemas
de inestabilidad. (Lopez, 1992)
Si en los edificios tradicionales de reducida esbeltez son las cargas gravitatorias las de
mayor relevancia, a medida que aumentan su altura, son las cargas laterales debidas a la accin
progresiva del viento las que adquieren una preponderancia creciente, por lo que requieren una
especial atencin. (Perles, 2003)
El avance de la tecnologa ha permitido desarrollar programas computacionales para el
anlisis estructural, la mayora de ellos se basan en anlisis elsticos de la estructura, o bien,
incorporan los efectos de segundo orden, mediante factores que amplifican las fuerzas obtenidas a
partir de un anlisis de primer orden . Cuando se trata de edificaciones de acero esbeltas, sin
arriostramiento y sometidas a grandes fuerzas axiales y cargas laterales como el viento, se deben
tener cuidados especiales, ya que stas pueden sufrir deformaciones que obliguen a incorporar la
no linealidad geomtrica en los clculos, ms bien dicho, realizar un anlisis de segundo orden
riguroso, debido a que estos efectos pueden ser significativos e influir en el comportamiento de la
estructura, aumentando los esfuerzos generados en los elementos que conforman el sistema
estructural.
Para obtener los esfuerzos internos de cada elemento de la estructura (fuerzas axiales,
fuerzas cortantes, momentos de flexin, etc.) se pueden realizar distintos tipos de anlisis los
cuales se mencionan a continuacin:

Anlisis elstico de primer orden

El anlisis elstico de primer orden es el comunmente empleado por los programas de
clculo matricial y tiene como inconveniente el que no pone de manifiesto la posible existencia
de la inestabilidad. Para este tipo de clculo:
|K]{ ] = {P] cc. 1.1.1
Donde:
K: Es la matriz de rigidez de la estructura
2

P: Son las acciones exteriores sobre la estructura
: Son las deformaciones calculadas por el programa

Como se puede observar, en ningn momento aparece la presencia de la inestabilidad,
dado que es proporcional a P, de manera que para asegurar la estabilidad de una estructura, por
este proceso de clculo, se debe verificar que en ninguno de los elementos de la estructura se
produzca pandeo.
En los mtodos de clculo elstico de primer orden, se analiza la estructura obteniendo los
esfuerzos sobre las barras y las longitudes de pandeo de stas, para posteriormente comprobar a
pandeo dichas barras. Lo que no puede determinarse por este mtodo es la carga crtica de pandeo
de la estructura, ni el coeficiente de seguridad respecto a dicha carga. (Hernndez et al, 1994).

Figura 1.1.1. Efecto P- y P-
Fuente: Elaboracin propia

Anlisis elstico de segundo orden

En el anlisis elstico de segundo orden el equilibrio se formula sobre la estructura
deformada. Este tipo de anlisis tiene en cuenta los momentos producidos por los esfuerzos de los
extremos de la barras combinados con los desplazamientos que se han producido en dichas
barras. Si estos desplazamientos son los de los extremos de la barras, se les denomina efecto P-,
si los desplazamientos son los que se producen en el interior de la barra, suponiendo que sus
extremos no han sufrido movimientos, se les denomina P- (figura 1.1.1). (Hernndez et al, 1994)
Como el principio de superposicin no es siempre vlido (ver nota 1), se debe acudir a un
reparto de cargas especificado, que se incrementa por pasos mediante un multiplicador de carga
3

(figura 1.1.2). Se escogen incrementos tan pequeos que permitan suponer un comportamiento
lineal durante este aumento de la carga.

Figura 1.1.2. Diagrama de carga-deformacin

La configuracin deformada que se obtiene al acabar cada aumento especificado de la
carga, es la geometra de referencia para el siguiente paso. As pues, la teora elstica de segundo
orden consiste en resolver una sucesin de anlisis de primer orden de una estructura cuya
geometra cambia en cada paso con respecto a los anteriores. Estos clculos pronto se hacen
inmanejables a mano y se necesitan programas informticos. Se estudia con mayor detalle este
tipo de anlisis en el captulo 2 de esta memoria.

Nota1: El Principio de superposicin constituye la base de gran parte del anlisis estructural. Puede
enunciarse como sigue: el desplazamiento o esfuerzo total en un punto de una estructura sometida a
varias cargas se puede determinar sumando los desplazamientos o esfuerzos que ocasiona cada una de
las cargas que acta por separado. Para que esto sea vlido es necesario que exista una relacin lineal
entre las cargas, esfuerzos y desplazamientos.
Requisitos para que el principio de superposicin sea aplicable:
1. El material estructural debe comportarse de manera elstico-lineal, a fin de que sea vlida la ley de
Hooke y la carga sea proporcional al desplazamiento.
2. La geometra de la estructura no debe sufrir cambios importantes cuando se apliquen las cargas. Si los
desplazamientos son grandes, entonces cambian considerablemente la posicin y orientacin de las
cargas. Un ejemplo es el caso de una columna sometida a una carga de pandeo. (Hibbeler, 1997).

Anlisis inelstico de primer orden

El anlisis inelstico de primer orden incluye los efectos de que algn elemento pueda
ceder de varias formas, pero normalmente mediante un modelo de rtula plstica. (Hernndez et
al, 1994).

4

Anlisis inelstico de segundo orden

Un clculo de segundo orden inelstico es el nico mtodo que se presenta como riguroso
en el anlisis estructural. Pero no existe en la actualidad ningn programa disponible para
oficinas de clculo que pueda realizar este tipo de anlisis. (Hernndez et al, 1994)

Figura 1.1.3. Tipos de anlisis

En el diseo basado sobre un anlisis elstico, es necesario emplear mtodos aproximados
para considerar los efectos de estabilidad sobre los elementos y el sistema interdependiente. El
mtodo tradicional de diseo para marcos de acero, usa los factores de longitud efectiva dentro de
las ecuaciones de interaccin viga-columna, para considerar el comportamiento del sistema.
(Maleck et al, 2004)
Este mtodo es ampliamente usado para el diseo de columnas de acero en los cdigos.
Pero esta basado en un modelo de columna aislada y no predice el comportamiento real de
estabilidad del marco como un todo, es decir, el pandeo de las columnas y la deformacin lateral
del marco. (Tong et al, 2007)

Figura 1.1.4. Interaccin entre el sistema estructural y los miembros que la componen
5

Se mencionan algunas dificultades asociadas con el uso del mtodo de longitud efectiva:

1. No puede considerar con exactitud la interaccin entre el sistema estructural y sus
elementos. Como resultado este mtodo no predice con exactitud los esfuerzos a los
cuales van a estar sometidos los miembros de una estructura.

2. En un anlisis elstico de primer orden, con los factores de amplificacin B
1
y B
2
,

(ver
nota 2) se consideran solo los efectos de segundo orden pero no la redistribucin de estos
esfuerzos en la estructura. Este mtodo proporciona una estimacin conservadora de la
capacidad ltima de carga que puede soportar una estructura.

Nota2: B
1
es un factor de amplificacin para considerar efectos de segundo orden causados por
desplazamientos entre puntos de intermedios (P-) y B
2
es un factor de amplificacin para considerar
efectos de segundo orden causados por desplazamientos en los extremos del elemento, cuando no esta
arriostrado (P-); ver figura 3.2.1.(AISC, 2005)

3. El mtodo requiere de un proceso que consume mucho tiempo, pues luego de realizar un
anlisis de primer orden de la estructura, se debe chequear la capacidad de cada uno de
sus elementos (vigas y columnas), involucrando el clculo de los factores de longitud
efectiva K. (Kim S. E. et al, 1999)

Por otro lado un anlisis mas exacto (incluyendo propiedades inelsticas del material) de
segundo orden debe incluir los factores que se enumeran a continuacin:

1. Las propiedades de los elementos en las diversas etapas del proceso de carga, que deben
definirse por medio de relaciones fuerza axial-momento-curvatura (M-P-) realistas.

2. Los cambios que el comportamiento inelstico y las articulaciones plsticas introducen en
las rigideces de los elementos.

3. Los cambios en las rigideces y desplazamientos de las columnas, ocasionadas por las
fuerzas axiales que actan sobre ellas.

4. Los momentos, fuerzas axiales y cortantes, resultantes de la interaccin de la cargas
verticales con los desplazamientos laterales de los entrepisos (efecto P-) y con los
desplazamientos del eje de las columnas con respecto a la recta que une sus extremos
(efecto P-), figura 1.1.1.

6

5. Los cambios en las longitudes de las columnas ocasionados por las fuerzas axiales y por
la deformacin por flexin del edificio en conjunto

6. Las deformaciones por cortante. (Lopez, 1992).

Un anlisis de este tipo, comunmente llamado Anlisis Avanzado, ofrece un medio mas
directo de evaluar los esfuerzos en los elementos y el sistema estructural, pero la capacidad
computacional y sofware requeridos para este tipo de anlisis, no se encuentra fcilmente
disponible para el diseo profesional en la actualidad.
Como mnimo se requiere de un simple anlisis elstico para un diseo preliminar rpido,
mientras que un mtodo avanzado de clculo se puede ocupar para resolver problemas poco
usuales, de tal forma de asegurar que los resultados sean correctos. (Maleck et al, 2004)
Por otra parte, la estabilidad de los elementos por separado y la estabilidad de la estructura
como un todo, puede ser tratado de manera mas exacta por la determinacin de los mximos
esfuerzos de la estructura atravs de un Anlisis elstico de Segundo Orden. (Kim S. E. et al,
1999)
A travs de este anlisis se pretende proporcionar un modelo que considere el fenmeno
de inestabilidad que afecta tanto a los elementos como al sistema estructural en forma directa. Las
ventajas de este tipo de anlisis son las siguiente:

1. Simplifica los procedimientos de diseo, al eliminar la necesidad de calcular los factores
de longitud efectiva. (Maleck et al, 2004)

2. Se pueden usar mtodos de diseo mas transparentes, estimando las fuerzas y resistencias
internas reales de los elementos dentro del sistema estructural. (Maleck et al, 2004)

3. Considera la estabilidad del sistema estructural y sus elementos directamente, por lo que
no se requiere chequear cada miembro por separado. (Kim S. E. et al, 1999).

4. Es un mtodo que puede ser programado, a diferencia de aquellos en los que se requiere
comprobar la capacidad de cada elemento de la estructura y el clculo de los factores de
longitud efectiva K. (Kim S. E. et al, 1999).

En nuestro pas, normalmente no se requiere llevar a cabo un anlisis de segundo orden
para las estructuras que se van a disear, esto debido a que la norma Nch 433, diseo ssmico de
edificios, no permite deformaciones mayores que 0.002H, siendo H la altura total o entrepisos de
un edificio; el Manual de diseo para estructuras de acero, preparado por el Instituto chileno del
7

acero (ICHA, 2000) establece que los efectos de segundo orden P- se deben considerar en los
casos en que los desplazamientos laterales totales o entre niveles de la estructura superen los
siguientes valores:

0.015 H/R para cargas normales ms ssmicas no mayoradas. (Norma NCh 2369, Diseo
Ssmico de Estructuras e Instalaciones Industriales). Siendo R el factor de modificacin
de la respuesta ssmica que vara entre 2 y 5.
0.004H para cargas normales mas viento no mayoradas.

A pesar de lo dicho anteriormente, es importante comenzar a investigar y tener
conocimientos tericos acerca del tema, debido a que el progresivo desarrollo y aumento de la
poblacin en Chile, obliga a construir edificios cada vez mas altos, y es justamente en este tipo de
construcciones en donde se generan los problemas de estabilidad, que deben ser tradados en
forma eficiente y responsable.
En esta memoria se realizar un Anlisis Elstico de Segundo Orden a una edificacin
tipo marco plano, que es una estructura usada a menudo en edificios y se compone de vigas y
columnas que estn articuladas o bien son rgidas en sus conexiones. La carga en un marco
ocasiona flexin en sus miembros y debido a las conexiones rgidas, esta estructura es
generalmente indeterminada desde el punto de vista del anlisis. (Hibbeler, 1997)
Se plantean tres mtodos de anlisis, los cuales consideran los efectos de la geometra
deformada en los clculos, estos son:

Mtodo 1

Anlisis elstico de segundo orden, utilizando la Matriz Geomtrica de la estructura, la
cual se genera a partir de las relaciones no lineales entre fuerzas y desplazamientos del elemento.
Estas se ven reflejadas en la matriz de rigidez no lineal de ste:
k
t
= k
c
+k
g
cc. 1.1.2
Donde
k
c
= Matriz de rigidez elstica convencional
k
g
= Matriz de rigidez geomtrica que depende no solamente de la geometra, sino
tambin de las fuerzas internas del miembro existente.

Mtodo 2

Se modela la estructura en un programa que realiza un Anlisis Elstico de segundo orden
en base a las ecuaciones diferenciales. El clculo se realiza en forma iterativa, determinando en
8

cada paso el valor de la fuerza axial, y se considera terminado cuando el incremento de la fuerza
axial en dos pasos consecutivos es menor que la tolerancia determinada.
El programa usado en este trabajo para ilustrar este mtodo es el programa llamado
Proges, y sus carctersticas sern analizadas el el captulo 4 de esta memoria.

Mtodo 3

Se modela la estructura en SAP2000, para luego analizarla utilizando la opcin que
permite incorporar el efecto P- en el anlisis.

Cabe mencionar que los mtodos propuestos solo se limitan a considerar la no linealidad
geomtrica, es decir, un Anlisis Elstico de segundo orden (ver figura 1.1.2), como hiptesis se
supone que el material trabaja en el rango elstico, por lo que no se considera la degradacin
plstica del material en los anlisis que se realizarn a la estructura.
Se dedica tambin un captulo a la revisin del AISC y el Eurocdigo 3 (EC 3), para
determinar cuando estas especificaciones demandan un anlisis de segundo orden de la
estructura, y los mtodos que se proponen en ambas especificaciones para incorporar estos
efectos en el anlisis.

9

1.2. Objetivos

Objetivos generales

Comparar resultados de esfuerzos y desplazamientos obtenidos a travs de un Anlisis
Elstico de Segundo Orden, realizado a un marco plano de acero no arriostrado, sometido
a grandes esfuerzos axiales y cargas laterales, por medio de los tres mtodos de clculo
propuestos.

Determinar si los esfuerzos (momento, fuerza axial y corte) en cada elemento de la
estructura aumentan en forma significativa, luego de efectuar un Anlisis Elstico de
segundo orden, y hasta que punto es recomendable realizar un Anlisis Elstico sin tomar
en cuenta el efecto que las cargas axiales junto con las cargas laterales puedan provocar
en una estructura.

Objetivos especficos

Aportar conocimiento y presentar las ventajas que conlleva realizar un analisis
considerando los efectos de segundo orden, para que los resultados sean cada vez mas
exactos, y de esta manera lograr que el proceso de diseo de cada elemento se lleve a cabo
en forma directa, sin que implique el clculo de distintos factores, para simular efectos
que en el anlisis no estn siendo considerados.

Dar a conocer las frmulas bsicas en las que estn basados los mtodos de clculo que
incorporan la geometra deformada, y mostrar como se lleva a cabo un anlisis de
segundo orden utilizando estas frmulas, mediante el mtodo de la Matriz Geomtrica.

Dar a conocer los resultados que entrega SAP2000, realizando un anlisis de segundo
orden a una estructura tipo marco, e identificar la manera de llevar a cabo este tipo de
anlisis por medio de este programa.

Determinar en que situaciones los cdigos de diseo que se estudian en esta memoria,
exigen un Anlisis Elstico de segundo orden de la estructura, y las condiciones que se
establecen para realizar los clculos.

10

1.3. Metodologa

Se da comienzo a esta memoria realizando un estudio de los tres mtodos de clculo
propuestos, (que se mencionan en la parte introductoria de esta memoria) para realizar un
Anlisis Elstico de segundo orden a una estructura, a travs de estos tres mtodos, la edificacin
ser un marco plano sometido a cargas laterales y axiales. (Ejemplo que se encuentra en el
capitulo 11, pgina 684 del libro llamado Mechanics of structures: Variational and
Computational Methods, Wunderlich et al, 2003)
Se establecen ciertas restricciones que debe tener este marco para que el modo de fallo sea
por pandeo general de Euler y no se produzca ninguna otra fuerza crtica, ya sea pandeo
inelstico, local, lateral o flexo-torsional, antes que la carga crtica de Euler.
Se define la carga lateral y las cargas axiales puntuales, a las cuales estar sometido el
marco no arriostrado, se tendr cuidado en que dicha carga axial no sea demasiado grande como
para hacer que la matriz se vuelva inestable (singular), ya que de ser as, se obtendran resultados
erroneos.
Se modela la estructura en SAP2000 y PROGES para realizar Anlisis Elstico de
segundo orden utilizando estos programas, luego se realiza el clculo utilizando el mtodo de la
Matriz Geomtrica, para posteriormente comparar los resultados obtenidos.

Capitulo II Mtodos de Anlisis Elstico de segundo orden
11

Captulo II
Mtodos de Anlisis Elstico de segundo orden.

En los ltimos aos se han producido numerosos esfuerzos para mejorar la exactitud del
anlisis de una viga-columna, basados en el mtodo de los desplazamientos, cuando se considera
la no linealidad geomtrica. Estos desarrollos se han centrado predominantemente en dos
aproximaciones: implementacin de los elementos finitos incorporando matrices de rigidez
geomtricas, o tratando de obtener funciones de estabilidad mediante la resolucin de la ecuacin
diferencial del equilibrio. (Ortega, 2007)
Se comienza este captulo citando los conceptos generales en los que se basa el anlisis
estructural, para luego estudiar el mtodo de anlisis de segundo orden utilizando la Matriz
Geomtrica en la seccin 2.2

2.1. Conceptos generales

2.1.1. Inestabilidad de una barra.

Se comienza considerando una barra rgida, con un resorte en su soporte y cargada
axialmente en su otro extremo (figura 2.1.1.1), luego se toma momento con respecto al punto A y
se iguala a cero:

Figura 2.1.1.1. Estabilidad de una barra rgida
Fuente: Saouma, 1999.
12

XH
A
= PA -k = u con A = I (ver nota 1)
PI -k = u => _P -
k
I
] = u cc. 2.1.1.1
Donde:
P = Es la fuerza axial en el extremo B de la barra
= Es el desplazamiento en el extremo B de la barra
= Es el ngulo de rotacin de la barra con respecto a su posicin inicial
k = Es la constante de rigidez del resorte en el extremo A de la barra

Nota1:=L solo para un desplazamiento pequeo ya que en estricto rigor =Lsen

Aqu hay tres posibilidades:

Equilibrio estable: P < k / L, = 0
Equilibrio indiferente: P = k / L, puede tomar cualquier valor
Equilibrio inestable: P > k / L, = 0

As se define la carga crtica como:
P
c
=
k
I
cc. 2.1.1.2

La primera ecuacin indica que la columna puede permanecer recta cualquiera sea el valor
de P, y la segundo muestra que cuando P alcanza un valor determinado, K / L, es posible un
estado de equilibrio con una configuracin ligeramente deformada; es decir, cuando la carga vale
K / L, pueden presentarse un nmero infinito de configuraciones en equilibrio, una recta y las
dems ligeramente deformadas; a este fenmeno, que corresponde el estado de equilibrio
indiferente, y consiste en la existencia de dos o mas posibles configuraciones en equilibrio para
una misma carga, se le llama Bifurcacin del Equilibrio, y es cacterstico de los fenmenos de
pandeo; la carga que lo ocasiona es la carga crtica del sistema P
cr
. (Lopez, 1992)
Para un gran desplazamiento se tiene:

XH
A
= PA -k = u con A = Iscn (ver nota 1)

PIscn -k = u cc. 2.1.1.S

P
c
=
k
Iscn
cc. 2.1.1.4
13

Figura 2.1.1.2. Curvas carga-rotacin de columnas

2.1.2. Pandeo elstico

Se presenta la teora de deformacin de la columna originada con Leonhard Euler 1744.
sta se basa en los supuestos de que un miembro inicialmente recto est cargado
concntricamente, y las fibras permanecen elsticas hasta que ocurre la deformacin, adems
cuando la columna se deforma, se debe asumir un pandeo simple como se muestra en la figura
2.1.2.1. Notar que la columna tampoco est perfectamente recta, y en todos los casos se presenta
una imperfeccin menor.

Figura 2.1.2.1. Columna Euler
Fuente: Saouma, 1999.

A partir de las ecuaciones A.1, A.2 y A.3 (ver anexo A) se pueden obtener las ecuaciones
2.1.2.1, 2.1.2.2, 2.1.2.3, introduciendo las siguientes condiciones:

H
Ax
= H
A
= H
Bx
= H
B
= u y P u

EI
x
:
II
+p: = u cc. 2.1.2.1
EI
u
II
+pu = u cc. 2.1.2.2
EC
u
II
-(0K
1
+K
)
I
= u cc. 2.1.2.S
14

Y para este caso en particular:
K
= -Pr
o
2
cc. 2.1.2.4
r
o
=
_
(I
x
+I
)
A
cc. 2.1.2.S
r
o
: cs cl roJio polor Jc giro Jc lo scccion

Derivando dos veces las ecuaciones 2.1.2.1, 2.1.2.2 y 2.1.2.3 (ver referencia Lopez,1992)
se obtiene lo siguiente:
:
Iv
+
P
EI
x
:
II
= u cc. 2.1.2.6
u
Iv
+
P
EI
u
II
= u cc. 2.1.2.7
Iv
+
Pr
0
2
-0K
1
EC
u

II
= u cc. 2.1.2.8

Estas tres ecuaciones son independientes entre s, pues en cada una aparece solo un
desplazamiento, lineal o angular, lo que indica que una pieza comprimida doblemente simtrica
puede mantenerse en equilibrio en la configuracin deformada descrita por cualquiera de ellas;
hay por consiguiente tres formas de pandeo independientes entre s, ya que en cada una de las
ecuaciones aparece una sola componente de desplazamiento, de manera que el pandeo se presenta
por flexin, alrededor de cualquiera de los dos ejes principales de las secciones transversales, o
por torsin, alrededor del eje longitudinal de la barra, y tres cargas crticas, que se calculan
resolviendo cada una de las ecuaciones de equilibrio. Las dos primeras, correspondiente a pandeo
por flexin respecto a cada uno de los ejes principales, se obtienen aplicando la teora de Euler,
pero sta no permite calcular la tercera. En un problema dado interesa conocer la menor de las
tres, que es la carga de Euler de pandeo por flexin alrededor del eje de menor momento de
inercia o la de pandeo por torsin. (Lpez, 1992)
Para el caso particular en que las secciones transversales tienen dos ejes de simetra se
llega a las ecuaciones 2.1.2.6, 2.1.2.7 y 2.1.2.8 que son independientes una de la otra, como se
mencion anteriormente. A continuacin se determina la carga crtica de pandeo, utilizando la
ecuacin 2.1.2.1, la cual tiene la siguiente forma de solucin:
: = C
1
scn_
P
EI
x
z +C
2
scn_
P
EI
x
z +C
3
z +C
4
cc. 2.1.2.9

Con
x
2
= PEI
x
la ecuacin 2.1.2.9 se escribe:

: = C
1
scn
x
z +C
2
scn
x
z +C
3
z +C
4
cc. 2.1.2.1u
15

Si la existencia de restricciones exteriores o la geometra de las secciones transversales es
simtrica, de manera que la columna se pandea mantaniendose en uno de sus planos de simetra,
se puede calcular la carga crtica igualando el momento exterior en una seccin transversal
cualquiera, ocasionado por la fuerza P, con el momento resistente interior que aparece en esa
misma seccin, y resolviendo la ecuacin correspondiente; por ejemplo, si la columna se flexiona
alrededor de su eje x se tiene que -P: = EI
x
R , y si se supone que los desplazamientos de su
eje son suficientemente pequeos para que la curvatura I/R pueda considerarse igual a J
2
:
Jz
2
= :
II
, se llega a
EI
x
:
II
+P: = u cc. 2.1.2.11

Esta constituye la ecuacin de equilibrio de la columna ligeramente deformada; su
solucin proporciona la carga que puede mantenerla en equilibrio en estas condiciones, es decir,
la carga crtica de pandeo elstico o carga crtica de Euler. (Lpez, 1992)
P
cx
=

2
EI
x
I
2
cc. 2.1.2.12

2.1.3 Pandeo inelstico

La obtencin de la frmula de Euler, que permite calcular la carga crtica de piezas rectas
comprimidas axialmente, est basada en la suposicin fundamental de que la pieza se comporta
elsticamente hasta la iniciacin del pandeo, por lo tanto las ecuaciones deducidas anteriormente
para P
cr
no son aplicables a columnas cortas o de longitud intermedia, en las que se alcanza el
lmite de proporcionalidad antes que el esfuerzo crtico de pandeo elstico. (Lpez, 1992)
La frmula
cr
=
2
E/(L/r)
2
es vlida nicamente para el intervalo de valores de la relacin
de esbeltez a los que corresponden esfuerzos crticos no mayores que el lmite de
proporcionalidad (
cr

LP
), de manera que es aplicable hasta que:

c
=
2
E
(kIr)
2

LP
ec. 2.1.S.1

Despejando L/r se obtiene
kI
r
_
E
LP
ec. 2.1.S.2
LP
= f
y
/ 2

: Es el esfuerzo correspondiente al lmite de proporcionalidad.
La ecuacin anterior permite calcular la relacin de esbeltez mnima para la que es
aplicable la frmula de Euler, la que deja de serlo para esbelteces menores puesto que para ellas
16

cr
>
LP
, es decir, el lmite de proporcionalidad se sobrepasa antes de iniciarse el pandeo y ste
se presenta en el intervalo inelstico. (Lpez, 1992)
Engesser y Considere fueron los primeros en advertir la posibilidad de modificar la
frmula de Euler para hacerla aplicable al clculo de la carga crtica de pandeo inelstico
introduciendo la teora del Mdulo Tangente.
Las hiptesis en las que se basa este modelo fueron expuestas por Shanley en 1947 y son
las siguientes:
Los desplazamientos laterales del eje de la columna son pequeos en comparacin con las
dimensiones de sus secciones transversales.
Las secciones transversales planas permanecen planas y normales al eje deformado,
despus de la flexin
El diagrama esfuerzo deformacin del material de la columna da la relacin entre el
esfuerzo y deformacin en cualquiera de sus fibras longitudinales
El plano de flexin es un plano de simetra de todas las secciones transversales.
(Lpez, 1992).

2.1.4. Pandeo lateral-torsional

En el desarrollo de las ecuaciones de diseo, el caso de una viga simplemente apoyada
con momento constante a travs de su seccin es usualmente usado como el caso bsico de
pandeo lateral torsional, ya que es el caso de carga ms desfavorable.
Consideremos una viga I, flexionada en su plano de mayor resistencia por la accin de dos
momentos iguales y de sentido contrario aplicados en sus extremos, como se muestra en la figura
2.1.4.1. Adems, admitiremos las siguientes hiptesis:

La seccin transversal de la viga es constante
Los esfuerzos normales mximos estn en intervalo elstico en el instante en que se inicia
el pandeo
La deformacin de la viga al flexionarse y retorcerse es de tal naturaleza que no cambia la
forma de sus secciones transversales
Las cargas exteriores permanecen paralelas a sus direcciones originales al desplazarse,
angular o linealmente, sus puntos de aplicacin

17

Figura 2.1.4.1 Deformacin lateral-torsional de viga en flexin
Fuente: Saavedra, 2005

Los extremos de la viga estn simplemente apoyados e impedidos de rotar en torno al eje
Z, pero libres de alabearse. En la condicin deflectada lateralmente de la viga el momento
externo M
x
sobre la seccin se proyecta en las direcciones , , , como se ve en la figura 2.1.4.2

H
{
= H
x
cos0cos[~ H
x
cc. 2.1.4.1
H
= H
x
scn0~ 0H
x
~H
x
Ju
Jz
cc. 2.1.4.2
H
q
= H
x
cos0scn[~ [H
x
cc. 2.1.4.S

Figura 2.1.4.2 Deformacin lateral-torsional de viga en flexin, detalle
Fuente: Saavedra, 2005

Las ecuaciones de equilibrio son:
EI
{
J
2
:
Jz
2
= H
{
= H
x
cc. 2.1.4.4
EI
q
J
2
u
Jz
2
= H
q
= [H
x
cc. 2.1.4.S
0[
J[
Jz
-EC
w
J
3
[
Jz
3
= H
= -
Ju
Jz
H
X
cc. 2.1.4.6
18

La primera de estas ecuaciones es independiente para v, mientras que la segunda y la
tercera estn acopladas, lo que nos muestra que los desplazamientos u y no pueden existir
independientemente uno del otro. Las ec. 2.1.4.5 y 2.1.4.6 pueden reducirse diferenciando la
segunda con respecto a z y eliminando d
2
u / dz
2
; ello nos conduce a la siguiente ecuacin
diferencial, que rige el fenmeno torsional:
EC
w
J
4
[
Jz
4
-0[
J
2
[
Jz
2
-
H
x
2
EI
[ = u cc. 2.1.4.7
Esta ecuacin permite una solucin distinta de la trivial del tipo:
[ = [
mux
scn
z
I
cc. 2.1.4.8

en que
max
es el ngulo de giro al centro de la luz. Sustituyendo la ec. 2.1.4.8 en la ec. 2.1.4.7 se
obtiene:
_
4
I
4
EC
w
+
2
I
2
0[ -
H
x
2
EI
_[
mux
scn
z
I
= u cc. 2.1.4.9

Dado que no es nulo, el parntesis de la ecuacin anterior debe ser nulo, conduciendo al
Momento Flector Crtico para el cual se inicia el pandeo lateral torsional.
Reduciendo, obtenemos la ecuacin diferencial de equilibrio, que rige el fenmeno
torsional:
H
c,x
=

I
_
EI
0[ +

2
I
2
E
2
I
C
w
cc. 2.1.4.1u
En que:
L = Largo total de la viga
E = Mdulo de elasticidad del acero
I
y
= Inercia de la seccin de la viga con respecto al eje y
G = Mdulo de corte del acero
J = Constante de torsin de la seccin
C
w
= Constante de alabeo de la seccin

Este momento de pandeo elstico corresponde a una deformacin lateral torsional (,u)
indeterminada. El primer trmino de la raz corresponde a la resistencia al pandeo lateral ofrecida
por la torsin de St. Venant (torsin pura) y la flexin lateral, mientras el segundo trmino
corresponde a la resistencia ofrecida por el alabeo de las alas. (Saavedra, 2005).

Se estudia a continuacin el primer mtodo de anlisis, que se usar posteriormente para
determinar los esfuerzos en la estructura propuesta.

19

2.2. Mtodo de la Matriz Geomtrica para incorporar efectos de segundo orden

Este mtodo consiste en la utilizacin de funciones aproximantes de tipo polinmico,
cbicas y lineales para los desplazamientos transversales de flexin y longitudinales debidos a los
esfuerzos axiales, respectivamente, que interpolan en cada nodo de la discretizacin ciertas
variables del problema consideradas como incognitas. (Ortega, 2007)
De esta forma se llega a la matriz de rigidez k
t
, que refleja las relaciones no lineales entre
fuerzas y desplazamientos del elemento, a travs de la matriz de rigidez del elemento no lineal
k
g
, que depende no solamente de la geometra, sino tambin de las fuerzas internas del elemento
existente. La matriz de rigidez total del elemento permite plantear la ecuacin de equilibrio de
ste.
|k
t
= k
c
+k
g
|

k
t
= p
+p
o

Donde:
k
t
= Matriz de rigidez total del elemento.

:
= Vector desplazamiento nodal del elemento referido al sistema de referencia que

se muestra a continuacin.
p
o
= Vector fuerza del elemento con extremos empotrados debido a la presencia de
carga intermedia.
p
= Vector de fuerza nodal como se muestra en la figura 2.2.1 .

k
c
= Matriz de rigidez elstica.
k
g
= Matriz de rigidez geomtrica.

p
=
`
1
1
1
1
N
u
I
u
H
u
N
b
I
b
H
b
1
1
1
1
1
=
`
1
1
1
1
u
u
w
u
u
u
b
w
b
b
1
1
1
1
1

Figura 2.2.1. Convencin de signos
Fuente: Wunderlich et al, 2003

20

Si se utilizan las Funciones de Estabilidad (ver anexo A), la matriz de rigidez en dos
dimensiones para un elemento viga-columna queda de la siguiente manera:

l
l
l
l
l
l
l
l
l
l
l
l
l
l
l
EA
l
u u -
EA
l
u u
u
|2(s
+s
]
) -e
2
]EI
l
3

-(s
+s
]
)EI
l
2
u
-|2(s
+s
]
) -e
2
]EI
l
3
-(s
+s
]
)EI
l
2

u
-(s
+s
]
)EI
l
2

s
EI
l
u
(s
+s
]
)EI
l
2

s
]
EI
l

-
EA
l
u u
EA
l
u u
u -
|2(s
+s
]
) -e
2
]EI
l
3

(s
+s
]
)EI
l
2
u
|2(s
+s
]
) -e
2
]EI
l
3

(s
+s
]
)EI
l
2

u
-(s
+s
]
)EI
l
2

s
]
EI
l
u
(s
+s
]
)EI
l
2

s
EI
l

1
1
1
1
1
1
1
1
1
1
1
1
1
1
1

Donde:
s
ii
y s
ij
son los coeficientes de rigidez del miembro obtenidos a partir de las Funciones de
Estabilidad elstica de una viga-columna. Estos coeficientes se pueden expresar de distintas
formas segn sea la fuerza de compresin o traccin, como:

Para una barra en compresin (N
0
< 0)
s
=
e(scne -cose)
2(1 -cose) -escne
s
= 4 -
2
1S
e
2
-
11
6Suu
e
4
-

s
]
=
e(e -scne)
2(1 -cose) -escne
s
]
= 2 +
1
Su
e
2
+
11
126uu
e
4
-

s
+s
]
=
e
2
(1 -cose)
2(1 -cose) -escne
s
+s
]
= 6 -
1
1u
e
2
-
11
14uu
e
4
-

e
2
= l
2
|N|
EI
(iguol o
2
cn
2(s
+s
]
) -e
2
=
e
2
scne
2(1 -cose) -escne
= 12 -
6
S
e
2
-
1
7uu
e
4
-

Para una barra en traccin (N
0
> 0)

s
=
e(scnbe -cosbe)
2(1 -cosbe) -escnbe
s
= 4 +
2
1S
e
2
-
11
6Suu
e
4
-

s
]
=
e(e -scnbe)
2(1 -cosbe) -escnbe
s
]
= 2 -
1
Su
e
2
+
11
126uu
e
4
-
21

s
+s
]
=
e
2
(1 -cosbe)
2(1 -cosbe) -escnbe
s
+s
]
= 6 +
1
1u
e
2
-
11
14uu
e
4
-

e
2
= l
2
|N|
EI

2(s
+s
]
) -e
2
=
-e
2
scnbe
2(1 -cosbe) -escnbe
= 12 +
6
S
e
2
-
1
7uu
e
4
-

Para obtener la matriz de rigidez k
t
se usan las aproximaciones polinmicas de las
Funciones de Estabilidad s
ii
y s
ij
que se dan anteriormente en cada caso; se utiliza solo hasta el
segundo trmino debido que de ah en adelante los trminos son pequeos, y no influyen
mayormente en la matriz por lo que se pueden despreciar; de esta manera la matriz de rigidez
total k
t
para un elemento esta expresada en trminos de k
c
(Matriz de rigidez Elstica) y
k
g
(Matriz de rigidez Geomtrica)
|k
t
= k
c
+k
g
|
As por ejemplo para un elemento empotrado en ambos extremos se tiene la siguiente
matriz geomtrica:

k
g
=
N
0
l
l
l
l
l
l
l
l
l
l
l
l
u u u u u u

u
6
S
-
l
1u

u -
6
S
-
l
1u

u -
l
1u

2l
2
1S

u
l
1u
-
l
2
Su

u u u u u u
u -
6
S

l
1u

u
6
S

l
1u

u -
I
1u
-
l
2
Su

u
l
1u

2l
2
1S

1
1
1
1
1
1
1
1
1
1
1
`
1
1
1
1
u
u
w
u
u
u
b
w
b
b
1
1
1
1
1

Esta matriz se agrega a la Matriz de rigidez Elstica que se utiliza normalmente para
anlisis de primer orden k
e
, que para un elemento empotrado en ambos extremos es la siguiente:

22

k
c
=
EI
l
3
l
l
l
l
l
l
Al
2
I
u
u
-Al
2
I
u
u
u
12
-6l
u
-12
-6l
u
-6l
4l
2
u
6l
2l
2
-Al
2
I
u
u
Al
2
I
u
u
u
-12
6l
u
12
6l
u
-6l
2l
2
u
6l
4l
2
1
1
1
1
1
1
`
1
1
1
1
u
u
w
u
u
u
b
w
b
b
1
1
1
1
1

Se da a continuacin la matriz geomtrica de cada elemento dependiendo de las
condiciones de apoyo.

Elemento empotrado-articulado

k
g
=
N
0
l
l
l
l
l
l
l
l
l
l
l
u u u u u u
u
6
S
-
l
S

u -
6
S
u

u
l
S

l
2
S

u
l
S
u

u u u u u u
u -
6
S

l
S

u
6
S
u

u u u u u u

1
1
1
1
1
1
1
1
1
1
`
1
1
1
1
u
u
w
u
u
u
b
w
b
b
1
1
1
1
1

Elemento articulado-empotrado

k
g
=
N
0
l
l
l
l
l
l
l
l
l
l
l
u u u u u u

u
6
S
-
l
S

u -
6
S
u

u u u u u u
u u u u u u
u -
6
S

l
S

u
6
S
u

u -
l
S

l
2
S

u
l
S
u

1
1
1
1
1
1
1
1
1
1
`
1
1
1
1
u
u
w
u
u
u
b
w
b
b
1
1
1
1
1

23

Elemento articulado-articulado

k
g
=
N
0
l
l
l
l
l
l
l
l
l
l
u u u u u u
u
1
4
u

u -
1
4
u
u u u u u u
u u u u u u
u -
1
4
u

u
1
4
u

u u u u u u

1
1
1
1
1
1
1
1
1
`
1
1
1
1
u
u
w
u
u
u
b
w
b
b
1
1
1
1
1

N
0
es negativo en compresin y positivo en traccin. (Wunderlich, 2003).

Notar que los valores que conforman la Matriz Geomtrica corresponden al segundo
trmino de la serie polinomial que representan las funciones de estabilidad dadas anteriormente,
pues el primer trmino de la funcin polinmica corresponde a la Matriz de Rigidez Elstica
usada sin considerar la carga axial.
Como alternativa a las aproximaciones que se vieron antes para las funciones de
estabilidad, se tiene un sistema de expresiones que usan series de expansin, stas tienen un alta
grado de convergencia dentro de los primeros diez trminos de la serie de expansin. De forma
mas simple, si =P/(
2
EI/L
2
) esta entre -2 y 2 se pueden utilizar las siguientes funciones:

s
= 4 +
2n
2
1S
-
(u.u1p +u.S4S)p
2
4 +p
+
(u.uu4p +u.28S)p
2
8.18S +p

s
]
= 2 -
n
2
p
Su
+
(u.u1p +u.S4S)p
2
4 +p
-
(u.uu4p +u.28S)p
2
8.18S +p

Estas expresiones se pueden aplicar tanto a miembros en traccin (P positivo) como en
compresin (P negativo). (Kim et al, 2004)

La limitacin que existe al usar estos campos de desplazamientos, es que los elementos
slo son capaces de representar curvaturas con variacin lineal o campos de deformacin axial
constante. En particular, la variacin de la curvatura a lo largo de la longitud de la pieza no se
modela con exactitud. En este sentido se emplean frecuentemente polinomios de mayor grado
asociados a grados de libertad internos al elemento para resolver estos problemas. (Ortega, 2007)

Captulo III Efectos de segundo orden considerados en las normas de diseo
24

Captulo III
Efectos de segundo orden considerados en las normas de diseo.

En este captulo se indica lo establecido en distintos cdigos de diseo, con respecto a la
estabilidad de la estructura, y cuando se requiere en cada una de estas normas, incorporar los
efectos de segundo orden en un anlisis realizado. Se comienza revisando el manual ICHA
(Instituto Nacional del Acero), luego el AISC (American Institute of Steel Construction) y por
ltimo el cdigo europeo EC3 (Eurocdigo 3), poniendo nfasis en los tipos de anlisis que se
proponen en cada una de las normas, para considerar los efectos de segundo orden en una
estructura. Finalmente se presenta un resumen con los conceptos mas importantes, para poder ver
de forma clara las diferencias o semejanzas que puedan existir entre ellas.

3.1. Recomendaciones del Instituto nacional del acero ICHA

A continuacin se profundiza el tema mencionado en la introduccin, respecto a como se
llevan a cabo los anlisis realizados en Chile, a travs de lo establecido en el manual de diseo
para estructuras de acero.
Los anlisis ssmicos tradicionales hechos en Chile no consideran el efecto P-,
cubrindolo por va:
a) De la calibracin que se ha hecho de los factores R (factor de modificacin de la respuesta
ssmica que vara entre 2 y 5) respecto de la experiencia ssmica concreta.
b) De la fijacin de un lmite a la deformacin horizontal mxima calculada.
c) Del establecimiento de un corte basal mnimo de diseo, que incrementa las fuerzas
ssmicas en las estructuras de perodos largos.

Esta prctica se ha demostrado adecuada para nuestra realidad, al no observarse fallas de
importancia producidas por el efecto P-.
En la norma NCh 433, Diseo Ssmico de Edificios, no se permite sobrepasar la
deformacin mxima H/500, de lo contrario se debe modificar la estructura, este lmite mximo
establecido para el desplazamiento, cuando se trata de cargas ssmicas, implica que no se requiere
un anlisis de segundo orden de la estructura. En la NCh 2369, Diseo Ssmico de Estructuras e
Instalaciones Industriales, se fija un lmite de 0,015H/R, ms all del cual se debe efectuar
anlisis P- y verificar que la estructura, los equipos y elementos no estructurales vinculados a
ella puedan aceptar deformaciones mayores. En ese caso, el anlisis se debe hacer con las fuerzas
ssmicas reducidas de la norma. Con esto no se busca, claro est, reproducir el fenmeno ssmico
real, sino darle a la estructura una resistencia adicional a la que resultara de no considerar el
efecto P-.
25

Concretamente, en cuanto a considerar o no los efectos de segundo orden en una
estructura, en el punto 6.1 del manual de diseo ICHA se establece los siguiente:
Los efectos de segundo orden, P- se considerarn en los casos en que los
desplazamientos laterales totales o entre niveles de la estructura superen los siguiente valores:
0.015H/R para cargas normales mas ssmicas no mayoradas (Nch 2369, diseo
ssmico de estructura e instalaciones industriales).
0.004H para cargas normales mas viento no mayoradas.
0.002H para cargas normales mas ssmicas no mayoradas (Nch 433, diseo ssmico
de edificio), ver nota 1
Donde:
H = Corresponde al alto total del edificio o entre niveles
R = Es el factor de modificacin de la respuesta ssmica, que vara entre 2 y 5

Nota1: En el caso del diseo ssmico de edificios, si se supera la deformacin mximo establecida de
0.002H se debe modificar la estructura

Lo establecido por el manual ICHA, respecto a los desplazamientos producidos por
fuerzas laterales se encuentra en el punto 15.3.3 de este manual y dice lo siguiente:
Las deformaciones o desplazamientos laterales debido a las cargas de viento o sismo, u
otras, no deben exceder los valores lmites especificados o tolerables (esto es realizando un
anlisis solo con las cargas laterales sobre la estructura):
Edificios con cargas ssmicas, Nch 433 0.002H
Estructuras industriales, Nch 2369 0.015H/R
Viento en edificios altos, prctica norteamericana 0.0025H (AISC, Wind & Seismic
Loads for Buildings, 1996)

Resumiendo y precisando lo dicho hasta aqu, se ve que el anlisis y la verificacin de las
columnas de marcos puede efectuarse de los siguientes modos:
Si la deformacin lateral ssmica calculada es menor que los lmites establecidos en las
normas NCh 433 y NCh 2369, no considerar el efecto P- y utilizar los valores K de los
Nomogramas (AISC) para determinar P
n
.
Si la deformacin lateral resulta mayor que esos lmites, pero la carga horizontal no es de
origen ssmico y ha sido estimada adecuadamente en su magnitud, efectuar directamente
un anlisis de segundo orden o P-, utilizando las frmulas de interaccin con:

c
P
n
= _
c
A
g
F
si no boy ponJco locol
c
A
g
F
c
si lo tcnsin crtico sc proJucc por ponJco locol
cc. S.1.1

26

Si la deformacin lateral es mayor que los lmites antes aludidos, y la carga es de origen
ssmico, con magnitud reducida por el factor R, efectuar el anlisis P- directamente con
el programa de anlisis estructural, o efectuar dos anlisis de primer orden equivalentes
(mtodo de amplificacin de momentos AISC). Los valores P
u
y M
u
obtenidos de este
anlisis, son comparados en las frmulas de interaccin con valores de P
n
y M
n
calculados
considerando las esbelteces de columna basadas en los factores K provenientes del
mtodo de los Nomogramas o de otros semejantes. (ICHA, 2000)
El anlisis de segundo orden puede ser reemplazado por dos anlisis de primer orden
equivalentes, aplicando el mtodo de amplificacin de momentos, ste es igual al que se propone
en el AISC y se estudia en la seccin 3.2.

Frmulas de interaccin:

Las frmulas que se dan a continuacin (seccin 10.1 del manual de diso ICHA), sirven
para verificar que el elemento resista la carga a la que estar siendo sometido, como se puede ver
K est implicitamente presente en esta frmula mediante el valor de P
n
, de esta manera se verifica
tambin la estabilidad del elementos cuando se ha realizado un Anlisis Elstico de primer orden.

Poro
P
u
c
P
n
u.2
P
u
c
P
n
+
8
9
_
H
ux
b
H
nx
+
H
u
b
H
n
_ 1.u ec. S.1.2
Poro
P
u
c
P
n
< u.2
P
u
2
c
P
n
+_
H
ux
b
H
nx
+
H
u
b
H
n
_ < 1.u cc S.1.S

P
n
= _
u.877
c
2
_ F
A
g
cc. S.1.4
c

`
1
1
1
1
cx
=
K
x
l
r
x
_
F
E
=
K
x
l
r
x
1
89.2
c
=
K
l
r
_
F
E
=
K
l
r
1
89.2
cc. S.1.S
Donde:
P
u
= Resistencia a la compresin requerida, (N)
P
n
= Resistencia nominal a la compresin, (N)
M
u
= Resistencia a la flexin requerida, (N-mm)
M
n
= Resistencia nominal a la flexin, (N-mm)
=
c
= Factor de resistencia a la compresin = 0,85
b
= Factor de resistencia a la flexin = 0.85 (para marcos no arriostrados)
K = Factor de longitud efectiva de pandeo
27

r = Radio de giro (mm)
E = Modulo de elsticidad (N-mm)
F
y
= Lmite de fluencia especificado (N-mm)
A
g
= Area bruta del elemento (mm
2
)
l = Longitud del elemento (mm)
x = Subndice relacionado a la flexin segn el eje mayor.
y = Subndice relacionado a la flexin segn el eje menor.

A travs de estos prrafos se ve que en los procesos de anlisis y diseo se plantean
situaciones distintas de acuerdo a la naturaleza y exactitud de las carga, o tipo de anlisis
realizado, los valores obtenidos nunca sern los reales ya que se aplican diversos factores,
haciendo que el proceso sea confuso e inexacto.
Por otra parte, se debe tener cuidado al confiar solo en los desplazamientos para
determinar si es necesario o no realizar un anlisis de segundo orden, pues cuando la estructura es
esbelta y esta sometida a cargas axiales de magnitud considerable, puede ser que un anlisis de
primer orden, al no considerar stas fuerzas, entregue resultados errneos de desplazamiento,
menores al lmite mximo de defomacin que permite la norma (lo que implica que no es
necesario un anlisis de segundo orden), pero que a travs de un anlisis de segundo orden, se
obtengan desplazamientos mayores, y mas an, que stos sobrepasen el lmite mximo, lo cual
indicara que las fuerzas axiales no se pueden despreciar.
Lo anterior puede darse en situaciones en que los desplazamientos estn en el lmite de los
valores mximos y la estructura este siendo sometida a una gran carga axial, por esta razn es
recomendable efectuar el anlisis de segundo orden cuando esto ocurra.

3.2. Mtodo de anlisis de segundo orden establecido en el cdigo AISC
3.2.1. Consideraciones generales

Los mtodos de diseo usualmente se basan en un anlisis elstico de la estructura, en este
cdigo se entregan tres mtodos de diseo los cuales son: En ASD (diseo por tensiones
admisibles), el clculo de esfuerzos est basado en un Anlisis Elstico de primer orden, y los
efectos de la geometra no lineal estn implcitamente considerados para los elementos, al
momento de realizar el diseo de cada uno de ellos. En el Diseo Plstico (PD) se permite la
redistribucin de fuerzas inelsticas a travs del sistema estructural. Ya que la no linealidad
geomtrica y los efectos graduales que se van produciendo no son considerados dentro del
anlisis, estos se aproximan en los elementos por medio de las ecuaciones de diseo, las cuales se
pueden encontrar en el cdigo AISC; en ella se aplican factores, que dependen de la longitud del
miembro, esfuerzo al que esta siendo sometido, condiciones de apoyo, tipo de arriostramientos,
28

etc. que reflejan las condiciones de esbeltez de un elemento, y se establecen lmites para calcular
el esfuerzo mximo que una pieza puede resistir en base a estos factores. En LRFD (diseo por
factores de carga y resistencia) un Anlisis Elstico de primer orden con factores de
amplificacin, o un Anlisis Elstico de segundo orden directo, es usado para considerar la no
linealidad geomtrica, y los esfuerzos ltimos de elementos viga-columna estn implcitamente
reflejados en el diseo por las ecuaciones de interaccin. (Kim S. E. et al, 1999).
En general, el anlisis para obtener las cargas de servicio (o efectos de cargas de momento
de pandeo, corte, fuerza axial y momento torsional) sobre los elementos est demostrado que es
el mismo tanto para LRFD como para ASD. Los mtodos elsticos de anlisis estructural son
usados excepto cuando el estado lmite es el mecanismo de colapso plstico. (Kim S. E. et al,
1999).
El caso ms comn en donde se debe considerar un anlisis de segundo orden es el de una
estructura de varios pisos, que debe confiarse en las rigideces de las interacciones entre vigas y
columnas, para resistir el pandeo bajo cargas laterales (viento o terremoto).
Este es el caso de marcos no arriostrados. En este caso los desplazamientos laterales
(tambin llamados pandeo o deformaciones) causan momentos de pandeo adicional debido a las
cargas de gravedad (P) actuando en las posiciones en donde se han producido los
desplazamientos . El anlisis debe incluir este efecto secundario P-. Hay varios grados de
sofisticacin que son usados en los anlisis para incluir estos efectos. Tanto en el LRFD como en
el ASD los efectos de segundo orden pueden ser calculados, como una parte del anlisis, o ellos
pueden ser considerados usando las frmulas que se dan en la especificacin. (Salmon G. C. et al,
1990).
En cuanto a la estabilidad de la estructura, en este cdigo se establece que debe ser
considerada desde el punto de vista de la estructura como un todo, incluyendo no solo los
elementos a compresin, sino tambin vigas, sistema de tirantes y conexiones. La estabilidad
individual de los componentes tambin debe ser considerada.

Figura 3.2.1.1: Efecto P- y P- en vigas columna
Fuente: AISC, 2005
29

En general, es esencial que un anlisis de segundo orden exacto sea realizado a la
estructura. El anlisis debera considerar la influencia de los efectos de segundo orden (tanto P-
como P- mostrados en la figura 3.2.1.1) y de flexin, corte y deformaciones axiales.

3.2.2. Anlisis de segundo orden por un anlisis elstico de primer orden amplificado.
Mtodo de amplificacin de momentos.

Segn el AISC este mtodo aplicado apropiadamente constituye un mtodo de anlisis de
segundo orden aceptable. Se definen factores de amplificacin B
1
y B
2
que son usados para
estimar los efectos de segundo orden producidos por las fuerzas de primer orden. Los efectos de
las cargas de primer orden que no estn asociados a pandeo lateral se multiplican por B
1
, para
estimar los efectos P- sobre los momentos M
nt
(momentos producidos cuando no hay
deformacin lateral del marco) y los efectos que se producen en presencia de carga lateral se
multiplican por B
2
, para estimar los efectos P- sobre los momentos M
lt
(momentos producidos
cuando hay deformacin lateral del marco). (ver fig 3.2.1.1). Los efectos de estos factores de
amplificacin sobre los momentos se muestran en la figura 3.2.2.1

Figura 3.2.2.1. Amplificacin de momento.
Fuente: Elaboracin propia.

Ambos tipos de momentos pueden ser inducidos por cargas de gravedad. M
nt
se define
como el momento desarrollado en un elemento de un marco arriostrado. M
lt
es el momento
desarrollado dentro de un elemento debido a la deformacin lateral de un marco. Si el momento
B
1
M
nt
se suma algebraicamente al momento B
2
M
lt
se obtiene una buena aproximacin al valor
M
r
resultante en la mayora de los casos. Se recomienda un anlisis de segundo orden de este tipo
si B
2
> 1.2

30

Para llevar a cabo este anlisis se deben realizar los clculos que se muestran a
continuacin. El esfuerzo de flexin M
r
, y el esfuerzo axial P
r
de segundo orden requerido sern
determinados como sigue:
H
= B
1
H
nt
+B
2
H
It
cc. S.2.2.1
P
= P
nt
+B
2
P
It
cc. S.2.2.2
Donde:
B
1
=
C
m
1 -oP
P
c1
1 cc. S.2.2.S

Para elementos sujetos a compresin axial B
1
puede ser calculado en base a una
estimacin de primer orden P
r
= P
nt
+ P
lt
.

Nota: B
1
es una amplificacin para considerar efectos de segundo orden causados por desplazamientos
entre puntos de intermedios (P-) y B
2
es una amplificacin para considerar efectos de segundo orden
causados por desplazamientos laterales de marcos no arriostrados (P-)(figura 3.2.1). Para miembros en
los cuales B
1
1.05, es conservador amplificar la suma de momentos con pandeo lateral y sin pandeo
lateral (como los obtenidos para el caso de anlisis elstico de primer orden) por el amplificador B
2
, en
otras palabras, Mr = B
2
(Mnt + Mlt).
B
2
=
1
1 -
o P
nt
P
c2
1 cc. S.2.2.4
Donde:
= 1 (LRFD)
= 1.6 (ASD)
M
r
= Esfuerzo de flexin de segundo orden requerido usando combinaciones de carga
LRFD o ASD, (Kip-in o N-mm)
M
nt
= Momento de primer orden usando combinaciones de carga LRFD o ASD,

asumiendo que no hay deformacin lateral en el marco, (Kip-in o N-mm)

M
lt
= Momento de primer orden usando combinaciones de carga LRFD o ASD,
causada por deformacin lateral en el marco, (Kip-in o N-mm)
P
r
= Esfuerzo axial de segundo orden requerido usando combinaciones de carga
LRFD o ASD, (Kips o N)
P
nt
= Fuerza axial de primer orden usando combinaciones de carga LRFD o ASD,
asumiendo que no hay traslacin lateral en el marco, (Kips. o N)
P
nt
= Carga vertical total soportada por el piso usando combinacin de de cargas
LRFD o ASD incluyendo cargas de gravedad de la columna, (Kips o N)
P
lt
= Fuerza axial de primer orden usando combinaciones de carga LRFD o ASD
causada solamente por la deformacin lateral del marco, (Kips o N)
31

C
m
= Coeficiente asumiendo que no hay deformacin lateral del marco, cuyo valor
deber ser tomado como:
(i) Para vigas columnas no sujetas a carga transversal entre soportes en el plano de pandeo.

C
m
= u.6 -u.4(H
1
H
2
) cc. S.2.2.S

Donde M
1
y M
2
, obtenidos a partir de un anlisis de primer orden, son los momentos menor y
mayor magnitud respectivamente, en los extremos de esa porcin del elemento no arriostrado en
el plano de pandeo bajo consideracin. M
1
/M
2
es positivo cuando los elementos tienen curvaturas
distintas (figura 3.2.3(b)), y negativo cuando tienen una curvatura (figura 3.2.3(a)).

Figura 3.2.2.2. Curvatura que ocasionan los momentos en los extremos de la barra

(ii) Para vigas columnas sujetas a carga transversal entre soportes, el valor de C
m
ser
determinado por cualquier anlisis o tomar conservadoramente como igual a 1 para todos los
casos.
P
e1
= Resistencia a pandeo crtico elstico del elemento en el plano de pandeo,
calculado en base al supuesto de que el marco est arriostrado, (Kips o N)
P
c1
=

2
EI
(K
1
I)
2
cc. S.2.2.6
P
e2
= La carga crtica de pandeo elstico para el piso determinado por el anlisis del
modo de pandeo, (Kips o N)

Para marcos sin arriostramiento, la resistencia de las columnas es igual a:

P
c2
=

2
EI
(K
2
I)
2
cc. S.2.2.7

Para todos los tipos de sistemas resistentes de carga lateral est permitido usar:

P
c2
= R
M
EI
A
H
cc. S.2.2.8

32

Donde:

E = Mdulo elstico del acero (29000 Ksi o 200000 Mpa)
R
M
= 1 para sistemas de marcos arriostrados;
= 0.85 para marcos de momento y sistemas combinados, a menos que se justifique un
valor mayor a travs de un anlisis
I = Momento de inercia en el plano de pandeo (in
4
o mm
4
)
L = Altura del piso en (in o mm)
K
1
= Factor de longitud efectiva en el plano de pandeo, calculado en base al supuesto de que
no hay deformacin lateral, igual a 1, a menos que un anlisis indique que un valor mas
pequeo debe ser usado
K
2
= Factor de longitud efectiva en el plano de pandeo, calculado en base al supuesto de que
el marco se puede deformar lateralmente.
H
= Deformacin de primer orden entre piso debido a fuerzas laterales, (in o mm).
H = Corte de piso producido por la fuerza lateral usada para calcular H, (Kips o N)

Este mtodo solo es vlido cuando cuando el valor de la relacin de amplificacin con
pandeo lateral
2nd
/
1st
orden (o B
2
), determinado a partir de un anlisis de primer orden, no
exceda 1.5, de lo contrario, AISC recomienda realizar un anlisis de segundo orden exacto de la
estructura. Por otra parte en el apndice 7 del AISC se propone un mtodo que permite usar los
factores de amplificacin incluso cuando B
2
> 1.5.

3.2.3. Mtodo de anlisis de segundo orden directo (Apendice 7 AISC).

Este mtodo considera aquellos parmetros potenciales que pueden afectar la estabilidad
de una estructura de acero tipo marco. Se consideran tres de los aspectos ms importantes del
comportamiento de estabilidad: no linealidad geomtrica, deformacin plstica, y estado lmite de
elementos. Estos aspectos gobiernan en ltima instancia las deformaciones de marcos bajo cargas
aplicadas y el resultado de los efectos de las cargas en la estructura.

No linealidad geomtrica e imperfecciones.

En el desarrollo y calibracin del mtodo de anlisis directo las imperfecciones geomtricas
iniciales son conservadoramente asumidas para las columnas igual a L/1000 donde L es la
longitud del elemento en el marco, y para marcos una deformacin de H/500 donde H es la altura
de piso.

33

Deformacin plstica.

Considera la deformacin plstica a travs de la seccin transversal y a lo largo del elemento. El
trmino se denomina esfuerzo residual y se supone que tiene un valor mximo de 0.3F
y
en el
borde para ser distribuido de acuerdo a la llamada Lehigh Pattern una variacin lineal en el borde
y tensin uniforme en el centro.

Cargas tericas

Las cargas tericas sern aplicadas al costado de los sistemas de marcos, para considerar los
efectos de las imperfecciones geomtricas, inelasticidad o ambos. Las cargas tericas son cargas
laterales que se aplican en cada nivel del marco, especificados en trminos de las cargas de
gravedad aplicadas en ese nivel. La carga de gravedad usada para determinar la carga terica ser
igual a la carga de gravedad mayor asociada con la combinacin de carga que est siendo
evaluada. Las cargas tericas sern aplicadas en la direccin que actan los efectos
desestabilizantes bajo la combinacin de carga especificada.

Condiciones de anlisis

(1) Los mtodos de anlisis en los que se desprecien los efectos P- sobre el desplazamiento
lateral de la estructura estn permitidos cuando las cargas axiales en todos los elementos en los
que se considera rigidez de flexin para contribuir a la estabilidad lateral de la estructura
satisfacen los siguientes lmites:

oP
< u.1SP
cL
cc. S.2.S.1

Donde:
P
r
= Esfuerzo de compresin axial requerido bajo combinaciones de cargas dadas por
ASD o LRFD, (Kips o N)
P
eL
=
2
EI /L
2

= 1 (LRFD) o 1.6 (ASD)

(2) Una carga terica, N
i
= 0.002Y
i
, aplicada independientemente en dos direccin ortogonales,
ser aplicada como una carga lateral en todas las combinaciones. Esta carga ser adicional a las
otras cargas laterales

Donde:
34

N
i
= Carga lateral terica aplicada en el nivel i, (Kips o N)
Y
i
= Carga de gravedad desde la combinacin de carga LRFD o 1.6 veces la
combinacin de carga ASD aplicada en el nivel i, (Kips o N)

(3) La rigidez de flexin reducida, EI*

EI
-
= u.8
b
EI cc. S.2.S.2
b
`
1
1
1
1
1 si
oP
u.S
4 _
oP
_1 -
oP
__ si
oP
> u.S
cc. S.2.S.S

Ser usada para todos los elementos cuya rigidez a la flexin se considere para contribuir
a la rigidez lateral de la estructura

Donde:
I = Momento de inercia con respecto al eje de pandeo, in
4
(mm
4
)
P
r
= Esfuerzo de compresin axial requerido bajo combinaciones de carga LRFD o
ASD, (Kips o N)
P
y
= Af
y
, esfuerzo producido en el miembro, (Kips o N)
= 1 (LRFD) = 1.6 (ASD)

En lugar de usar
b
< 1 donde P
r
/ P
y
> 0.5,
b
= 1 puede ser usado para todos los
elementos, lo que proporciona que una carga terica adicional de 0.001Yi sea agregada a la carga
terica requerida en (2).

(4) Una rigidez axial reducida, EA*

EA
-
= u.8 EA cc. S.2.S.4

Ser usada para elementos cuya rigidez axial este considerada para contribuir a la
estabilidad lateral de la estructura, donde A es el rea de la seccin del elemento.

Aplicando estas condiciones el diseo puede llevarse a cabo de la siguiente manera:
Para el diseo de los elementos, luego de un anlisis de este tipo se debe verificar para cada
seccin lo establecido en el capitulo H de este cdigo, determinando el esfuerzo de la columna
nominal, P
n
, usando K = 1 (notar que no es necesario el clculo de K).

35

Ecuaciones de interaccin

En en captulo H del AISC se establece lo siguiente:
La interaccin de flexin y compresin en los elementos para los cuales 0.1 (I
yc
/I
y
)
0.9, que estn forzados a pandearse en uno de sus ejes (x y/o y) estarn limitados por las
ecuaciones 3.2.3.5 y 3.2.3.6, donde Iyc es el momento de inercia con respecto al eje y del ala en
compresin, in
4
. (mm
4
).

Poro
P
P
c
u.2
P
P
c
+
8
9
_
H
x
H
cx
+
H
H
c
_ 1.u ec. S.2.S.S

Paia
P
P
c
< u.2
P
2P
c
+_
H
x
H
cx
+
H
H
c
_ 1.u ec. S.2.S.6

Donde:
P
r
= Fuerza de compresin axial requerida, (Kips o N)
P
c
= Fuerza de compresin axial admisible (P
n
o P
n
/ 0 LRFD o ASD
respectivamente), (Kips o N)
P
n
= F
cr
A
g

M
r
= Fuerza de flexin requerida, (Kips-in o N-mm)
M
c
= Fuerza de flexin admisible, (Kips-in o N-mm)
x = Subndice relativo al eje fuerte de pandeo
y = Subndice relativo al eje dbil de pandeo

Y luego para el clculo de F
cr
se establecen las siguientes frmulas basadas en el factor de
longitud efectiva K
F
c
=
`
1
1
1
1
|u.6S8
P
j
P
c
]F
si
KI
r
4.71_
F
c
F
o (F
c
u.44F
)
u.877F
c
si
KI
r
> 4.71_
F
c
F
o (F
c
> u.44F
)
cc. S.2.S.7
Con:
F
c
=
n
2
E
[
KI
r
2
cc. S.2.S.8

Se permite por lo tanto el clculo de F
e
con K = 1 en la ecuacin 3.2.3.8 si el anlisis
estructural incorpora los efectos de segundo orden de manera rigurosa, de esta manera P
n
ya no
depende del factor K.
36

Se sugiere que en la mayora de las estructuras formadas por barras, los amplificadores de
pandeo lateral de segundo orden (o el equivalente B
2
), calculado con la rigidez reducida, debe ser
conservadoramente no mayor a
2nd
/
1st orden
= 2.5 ya que en los niveles ms grandes de
amplificacin, pequeos cambios en las cargas de gravedad o rigideces resultan en grandes
cambios en relacin a deflexiones del modo de pandeo y fuerzas de segundo orden internas,
debido a la no linealidad. (AISC, 2005)

3.3. Cdigo Europeo, EC3
3.3.1. Consideraciones generales

En la siguiente seccin se da a conocer la forma que tiene el Eurocdigo 3 para tratar el
asunto de la estabilidad de un marco plano.
Acorde con en Eurocdigo 3 (2005), la verificacin de la estabilidad de los prticos o de
los elementos individuales, debe ser llevada a cabo considerando las imperfecciones y los efectos
de segundo orden.
Para la verificacin de la estabilidad estructural, se debe optar por uno de los tres
procedimientos descritos en el apartado 5.2.2(3) del Eurocdigo 3.
El primer procedimiento consiste en realizar un anlisis global de todo el prtico,
considerando los efectos de segundo orden (P- y P-), las imperfecciones del prtico
traslacional y de cada elemento. Con este procedimiento solo se requerir verificar la resistencia
de la seccin mas solicitada.
El segundo procedimiento permite realizar un anlisis global de todo el prtico,
considerando parcialmente los efectos de segundo orden (solo P-) y las imperfecciones del
prtico traslacional. Como el efecto P- y las imperfecciones de los elementos individuales no
estan siendo considerados en el enlisis global, la estabilidad de cada elemento deber ser
verificada y adicional a esto, es necesario verificar tambin la resistencia de la seccin.
Por ltimo el tercer procedimiento propone emplear el mtodo de la columna equivalente
para evaluar la estabilidad del prtico y de los elementos estructurales. Con este procedimiento
los esfuerzos internos que deben ser usados tanto en la verificacin de la resistencia de la seccin
como en la verificacin de la estabilidad de los elementos, son calculados de acuerdo a la teora
de primer orden sin considerar imperfecciones. (Yong, 2007)
En el punto 5.2.1.2 del eurocodigo, en cuanto a la influencia de las deformaciones, se
establece lo siguiente:
Se podr utilizar teora de primer orden en los siguientes casos:
Prticos arriostrados
Prticos intraslacionales
37

Mtodos de clculo en los que se consideren indirectamente los efectos de segundo orden
(con el factor de amplificacin e imperfeccin inicial) (Eurocodigo, 2003)

Tipos de estructuras definidas en esta especificacin
En el Eurocodigo 3 se clasifican las estructuras en 3 tipos:

a) Elementos estructurales simples

Las vigas simplemente apoyadas y las piezas aisladas traccionadas o comprimidas son isostticas.
Las celosas pueden ser isostticas o hiperestticas.

b) Vigas continuas y prticos arriostrados

En el caso de vigas continuas y prticos en los que los efectos de segundo orden sean
despreciables o se hayan eliminado con medios adecuados, se analizarn las distribuciones de
sobrecargas que produzcan las solicitaciones ms desfavorables para la resistencia de las piezas
aisladas y la de sus uniones.

c) Prticos no arriostrados

Estos deben analizarse para la distribucin de sobrecarga ms desfavorable teniendo en cuenta los
desplazamientos de los nudos. Asimismo se analizar la resistencia de las piezas aisladas tal
como se describe en b) (Eurocdigo 3, 2003)

Imperfecciones

En el anlisis de la estabilidad se deben tener en cuenta los efectos de las imperfecciones,
los cuales incluyen las tensiones residuales y las imperfecciones geomtricas tales como la falta
de verticalidad de los prticos y la falta de rectitud en los elementos.
La imperfeccin equivalente a nivel de prtico (figura 3.3.1.1) producida por una fuerza
lateral, puede ser determinada por medio de la siguiente expresin:

=
o
h
o
m
2uu
cc. S.S.1.1
o
h
=
2
b
pcro
2
S
o
h
1.u cc. S.S.1.2
o
m
=
_
u.S _1 +
1
m
] cc. S.S.1.S
38

Donde:
h = Altura (m) de la estructura
m = Nmero de pilares por piso, incluyendo solo aquellos que llevan una carga
superior al 50% de la carga total

Figura 3.3.1.1. Imperfeccin de prtico traslacional
Fuente: Yong, 2007

Por otro lado, la imperfeccin a nivel del elemento tiene una forma sinusoidal, y su
amplitud e
o
depende de las curvas de pandeo, del tipo de anlisis a realizar y de la longitud del
elemento.

Figura 3.3.1.2. Imperfeccin en el elemento
Fuente: Eurocdigo, 2003
39

Curva de pandeo Anlisis elstico Anlisis plstico
ver tabla 3.3.1.2 e
o
/ L e
o
/ L
a
o
1/350 1/300
a 1/300 1/250
b 1/250 1/200
c 1/200 1/150
d 1/150 1/100
Tabla 3.3.1.1. Imperfecciones de los elementos
Fuente: Yong, 2007

Curva de pandeo a
o
a b c d
Factor de imperfeccin 0.13 0.21 0.34 0.49 0.76
Tabla 3.3.1.2. Factores de imperfecciones para las curvas de pandeo
Fuente: Yong, 2007

Cuando se analicen los pandeos dentro y fuera del plano, todas estas imperfecciones
deben ser consideradas en la direccin mas desfavorable. Estas imperfecciones se pueden tomar
en cuenta en forma de fuerzas ficticias como se muestra en la figura 3.3.1.2 y 3.3.1.3:

Figura 3.3.1.3. Esfuerzos horizontales equivalentes
Fuente : Eurocodigo 3, 2003

3.3.2. Anlisis de la estabilidad de un marco y factor de amplificacin de momentos de
primer orden.

Para realizar un anlisis de primer orden elstico es necesario que el marco sea
intraslacional, esto ocurre si se cumple una de las dos condiciones siguientes:
(1) Marco intraslacional cuando el prtico esta debidamente arriostrado.
Un marco puede considerarse arriostrado cuando su rigidez est asegurada por un sistema
de arriostramiento de forma que su respuesta a cargas horizontales en su plano es
40

suficientemente rgida para que se pueda considerar que todas las cargas horizontales son
resistidas por el sistema de arriostramiento.

(2) Tambien se puede considerar intraslacional dependiendo magnitud de las cargas que se
estn aplicando a la estructura, y se determina de la siguiente forma:

Se clasificar un prtico como intraslacional cuando su respuesta a los esfuerzos
horizontales contenidos en su plano es tal que sean despreciables las solicitaciones
adicionales debidas a desplazamientos horizontales de sus nudos. esto se verifica con la
siguiente condicin

I
sd
I
c
u.1 ec. S.S.2.1
Donde:
V
sd
= Valor de clculo de la carga vertical total
V
cr
= Valor elstico crtico de la carga vertical total considerando un estado
traslacional.

Para los prticos planos de edificacin con vigas unidas a los pilares en cada planta (fig
3.3.1) se puede utilizar la expresin 3.3.2.2 para el valor de V
sd
/V
cr
, por lo tanto:
Cuando se utilice la teora de primer orden, los desplazamientos horizontales de
cada planta producida por las cargas horizontales y verticales de clculo, y las
imperfecciones debidas a un desplazamiento horizontal inicial, aplicadas como esfuerzos
horizontales equivalentes, debern cumplir el criterio:

o
c
=

b
x
I
E
u.1 ec. S.S.2.2
Donde:

= Desplazamiento horizontal entre piso producido por las fuerzas laterales mas la
fuerza ficticia
h = Altura entre plantas
H = Reaccin horizontal total en el piso producido por las fuerzas laterales
V = Carga vertical de diseo sobre la estructura en un piso

41

Figura 3.3.2.1. Prtico de edificacin con las vigas unidas a los pilares en cada planta
Fuente: Eurocdigo 3, 2003

En caso de no cumplirse (1) ni (2) se considera el prtico como traslacional, y en su
dimensionamiento se tendrn en cuenta los efectos de los desplazamientos horizontales de sus
nudos y tambin debern cumplirse los requisitos relativos a la estabilidad global. (Eurocdigo 3,
Dec 2003)

Mtodo de amplificacin de momentos

En el mtodo propuesto por este cdigo, los momentos resultantes del anlisis elstico de
primer orden se incrementarn multiplicndolos por el cociente:
F
A
=
1
[1 -
I
sd
I
c
cc. S.S.2.S
Donde:

V
sd
V
cr
= Valor elstico crtico de la carga vertical total considerando un estado no
arriostrado.
F
A
= Factor que amplifica los esfuerzos obtenidos luego du un anlisis de primer
orden.
Los momentos debidos al desplazamiento son los asociados a los desplazamientos
horizontales de una planta respecto a la inferior. Estn originados por la carga horizontal y
tambin pueden deberse a la carga vertical si la estructura o la carga es asimtrica.
En los prticos de edificacin (fig. 3.3.1), se puede emplear tambin la siguiente
aproximacin, para determinar V
sd
/ V
cr

I
sd
I
c
=

b
x
I
E
cc. S.S.2.4
42

No se utilizar el mtodo de amplificacin de momentos cuando el cociente de las cargas
V
sd
/ V
cr
sea mayor que 0.25. (Eurocdigo 3, 2003).

3.4. Resumen y comparaciones

En cuanto al tema de la estabilidad ambos cdigos coinciden en que sta debe ser tratada
como un todo en la estructura, considerando tanto el efecto P- como el P- producto de las
fuerzas axiales y laterales actuando sobre las deformaciones de los elementos.
El cdigo AISC utiliza los factores B
1
y B
2
para tratar ambos efectos, sin embargo no
incluye deformaciones iniciales en el anlisis elstico de primer orden, cuyos resultados sern
amplificados por stos valores. La normativa Europea utiliza un factor de amplificacin F
A
para
considerar los efectos de segundo, adems de incluir el los clculos imperfecciones iniciales,
tanto para cada elemento como para el desplazamiento lateral del marco, en el anlisis realizado.
El factor de amplificacin B
2
del AISC es similar al entregado por el EC3 F
A
, pero con
ambas frmulas se llega a resultados distintos, el lmite mximo establecido a partir del cual es
necesario realizar un anlisis de segundo orden de acuerdo al AISC es cuando B
2
1.2 mientras
que para el EC3 este lmite se da cuando V
sd
/ V
cr
0.1 lo que en la frmula 3.3.2.3 da un valor
de F
A
= 1.11, a partir de stos valores se debe realizar un anlisis de segundo orden mediante la
amplificacin de momentos como se especifica en cada normativa, sin embargo cuando B
2
> 1.5,
para el AISC, y V
sd
/ V
cr
> 0.25 F
A
=1.33, para EC3, se debe llevar a cabo un anlisis de
segundo orden mas riguroso.
En el manual de diseo ICHA se trabaja con el mismo mtodo empleado por el AISC,
tanto para los anlisis realizados (aplicando anlisis de primer orden o la amplificacin de
momentos) como para el diseo de los elementos basados en las ecuaciones de interaccin. La
diferencia est en que las normas chilenas son mas restrictivas en cuanto a la deformacin lateral
de las estructuras y stas son las que determinan si es necesario o no analizar una determinada
estructura considerando los efectos de segundo orden. En la tabla 3.4.1 se puede ver en forma
resumida lo establecido en cuanto a las deformaciones.

Origen Norma Deformacin Anlisis Verificacin
Smico Nch 433 < 0,002H No p- usar K
Nch 2369 < 0,015H/R No p- usar K
Ssmico Nch 433 > 0,002H Modificar la estructura
Nch 2369 > 0,015H/R p- usar K
No ssmico Nch 433 > 0,002H p- sin k
(viento) Nch 2369 > 0,015H p- sin k
Tabla 3.4.1. Tipos de anlisis de acuerdo al origen de las fuerzas y deformacin
43

Ntese pues la diferencia fundamental entre el anlisis P-, cuanto las fuerzas
horizontales intervienen en sus magnitudes reales, y el utilizable en la verificacin ssmica.
Debido a que Chile es un pas ssmico, y estas fuerzas no son fciles de estimar se tienen
mayores restricciones y son estas fuerzas que generalmente gobierna a la hora de realizar un
anlisis y verificar la estructura.

Captulo IV Programas computacionales
44

Captulo IV
Programas computacionales

En este captulo se estudian dos programas computacionales, los cuales se van a usar para
realizar el Anlisis Elstico de segundo orden de la estructura, se comienza con el programa
SAP2000 para determinar de que forma se puede llevar a cabo este tipo de anlisis utilizando este
programa para finalizar en la seccin 4.2 dando a conocer las bases de clculo que se utilizan en
Proges.

4.1. SAP2000
4.1.1. Introduccin

Cuando las cargas que actan sobre la estructura y los resultados de las deformaciones son
bastante pequeos, la relacin carga-deformacin para la estructura es lineal. Para la mayor parte
de los anlisis, SAP2000 supone un comportamiento lineal, esto permite formar las ecuaciones de
equilibrio usando la geometra original (no deformada) de la estructura. Estrictamente hablando,
las ecuaciones de equilibrio debieran referirse a la geometra de la estructura luego de la
deformacin.
Las ecuaciones de equilibrio lineal son independientes de la carga aplicada y los
resultados de la deformacin. As los resultados de diferentes cargas dinmicas o estticas se
pueden superponer, siendo esto eficiente desde el punto de vista computacional.
Si las cargas sobre la estructura y los resultados de las deformaciones son grandes,
entonces el comportamiento carga-deformacin pasa a ser no lineal.
SAP2000 es capaz de considerar la no linealidad geomtrica en forma de efectos P-,
mediante un proceso iterativo en un anlisis no lineal esttico.

Efecto P-: cuando se presentan grandes esfuerzos dentro de una estructura (fuerzas y
momentos), las ecuaciones de equilibrio escritas para la geometra original y para la
deformada pueden diferir significativamente incluso si las deformaciones son pequeas.

Para muchas estructuras de edificios, el efecto P- generalmente ocurre en las columnas
debido a la carga muerta y viva. Las fuerzas axiales en las columnas estn comprimiendo los
elementos, haciendo que la estructura sea ms flexible ante las cargas laterales. (Computers and
structures, Inc. Abril 2007)
SAP2000 considera el efecto P- de forma diferente para los componentes de marcos
arriostrados y no arriostrados. Para los marcos arriostrados, el efecto P- se limita a la estabilidad
individual del elemento. Para los elementos de un marco no arriostrado, los efectos de la
45

deformacin lateral, deben ser adicionados a los efectos de la estabilidad de cada miembro.
SAP2000 asume que los momentos producidos por el efecto P- en marcos arriostrados son
producto de las cargas muertas o vivas sobre la estructura. Mientras que en los marcos no
arriostrados, son producto de otro tipo de cargas como son las de viento o sismo.
Para los efectos de estabilidad de un miembro individual y la deformacin lateral en
marcos no arriostrados, los momentos son amplificados por los factores B
1
y B
2
tal como se
especifica en el cdigo AISC en el mtodo LRFD y AASHTO-LRFD.
El usuario de SAP2000 debe saber que la opcin de anlisis por defecto en este programa
no contempla el efecto P- ya que solo realiza una iteracin para este efecto. El usuario debe
activar la opcin de analisis P- y el nmero mximo de iteraciones para el anlisis, ver figura
4.1.1.1 (Computers and Structures, Inc, 2000)

Figura 4.1.1.1. Efecto P- en SAP2000

Lo que hace el efecto P- esencialmente es modificar las carcteristicas de la estructura,
afectando los resultados de todos los otros anlisis realizados. Los supuestos mas importantes en
los que se basa SAP2000 para realizar este tipo de anlisis se dan a continuacin:
El efecto P- es considerado solamente en elementos de tipo marco.
Se considera solamente el efecto de un esfuerzo de magnitud considerable provocado por
la fuerza axial sobre la deformacin transversal por corte y momento.
Se asume que todas las deformaciones y rotaciones van a ser pequeas
46

Se asume que la forma de la deformada transversal de un elemento del marco es cbica
cuando es producida por momentos y lineal cuando es producida por corte
Se asume que las fuerzas axiales P- van a ser constantes a lo largo del elemento.
(Computers and Structures, 1998)

4.1.2. Forma cbica deformada

Si el equilibrio es examinado en la configuracin original (usando la geometra sin
deformar), el momento en la base es M=FL, y decrece linealmente a cero en el extremo cargado
de la columna. Si en lugar de esto, el equilibrio est considerado en la configuracin deformada,
hay un momento adicional causado por la fuerza axial P actuando sobre el desplazamiento
transversal .

Figura 4.1.2.1. Geometra deformada
Fuente: Computers and structures, 2007

El momento ya no vara linealmente a lo largo de la longitud; la variacin depende del
lugar sobre la forma deformada. El momento en la base es ahora M=FL+P.
Si la fuerza de compresin es demasiado grande, la rigidez transversal tiende a cero y la
deformacin tiende a infinito, por lo tanto la estructura tiende a pandearse. El valor terico de P
en el cual ocurre esto es llamado Carga de Pandeo de Euler para una columna y est dado por la
siguiente frmula:
P
c
= -
2
EI
k
2
I
2
cc. 4.1.2.1
El efecto P-Delta exacto de la carga axial sobre la deformacin transversal y rigidez es
una funcin mas bien complicada de la relacin fuerza P para la carga de pandeo P
cr
. La
verdadera forma deformada de la viga, y por lo tanto el efecto sobre el diagrama de momento, se
describe por funciones cbicas bajo carga axial cero, funciones hiperblicas bajo traccin, y
funciones trigonomtricas bajo compresin.
47

El efecto PDelta se puede presentar en cualquier configuracin de viga, ya sea
simplemente apoyada, empotrada-empotrada, etc. el efecto P-Delta se puede aplicar localmente o
en elementos individualmente, o globalmente al sistema estructural como un todo.
La caracterstica clave es que una gran fuerza axial, actuando sobre una pequea
deformacin transversal, produce un momento significativo que afecta al comportamiento del
elemento o la estructura. Si la deformacin es pequea, entonces el momento producido es
proporcional a la deformacin. (Computers and Structures, Oct 1998)

4.1.3. Fuerzas P- en los elemento prtico

La forma cbica de la deformada es integrada a lo largo de la longitud de cada elemento
de marco, tomando en cuenta la deformada dentro del elemento. Para este propsito la forma
deformada transversal se asume como cbica para flexin y lineal para corte entre los extremos
rgidos del elemento.
La verdadera forma deformada puede diferir algo de la deformacin cbica lineal asumida
en las siguientes situaciones:
Los elementos con propiedades de seccin no prismticos. En este caso la deformada P-
Delta se calcula como si el elemento fuera prismtico usando el promedio de las
propiedades de todo el largo del elemento.
Cargas actuando a lo largo de la longitud del elemento. En este caso la deformada P-delta
se calcula usando extremos empotrados y la fuerza aplicada en los extremos del elemento.
Una fuerza grande P-Delta actuando sobre el elemento. La verdadera forma deformada es
realmente la descrita por funciones trigonomtricas bajo una compresin grande, y por
funciones hiperblicas bajo una traccin grande.

El supuesto de la forma cbica es usualmente una buena aproximacin excepto bajo
fuerzas de compresin cerca de la carga de pandeo con ciertas restricciones en los extremos. Sin
embargo se pueden obtener excelentes resultados dividiendo los miembros estructurales en 2 o
ms partes.

Fuerzas axiales P-
SAP2000 modifica las ecuaciones de equilibrio para incluir el efecto P- debido a un solo
sistema de fuerzas axiales actuando en los elementos del modelo de marco.
Estas ecuaciones modificadas son usadas para todos los anlisis siguientes que son
aplicados en la estructura. Los resultados de estos diferentes anlisis se pueden superponer.
Las fuerzas axiales que crea el efecto P- vienen del clculo de un anlisis iterativo
realizado a la estructura, con una combinacin de cargas estticas especficas.
48

Anlisis iterativo
Se realiza un anlisis preliminar para determinar todas las fuerzas axiales en la estructura.
Las ecuaciones de equilibrio son luego reformuladas y resueltas nuevamente tomando en cuenta
estas fuerzas axiales. El segundo anlisis puede producir fuerzas axiales diferentes en los
elementos si las rigideces modificadas causan una redistribucin de las fuerzas. Adems de la
iteracin, se requiere de cada ecuacin de equilibrio reformulada y resuelta hasta que las fuerzas
y deformaciones convergen, es decir, hasta que no hayan cambios significativos entre una
iteracin y la siguiente.

Criterio de convergencia
Se puede especificar una tolerancia relativa al desplazamiento, con la cual medir la
convergencia. Este valor, por defecto, es 0.003. si el cambio relativo entre el deplazamiento de
una iteracin a la siguiente es menor que la tolerancia, entonces no se realizan mas iteraciones. El
cambio relativo entre los deplazamientos est definido como la relacin entre el cambio mximo
de desplazamiento y el mayor desplazamientos producido en cualquier iteracin. Cabe destacar
que rotaciones y traslaciones son tratadas de igual manera. (ver figura 4.1.1.1)

Nmero mximo de iteraciones
Se debe especificar un nmero mximo de iteraciones que el programa debe realizar. Esto
se usa para prevenir el tiempo de clculo excesivo. La iteracin inicial es llamada la iteracin
zero-th. ste es el anlisis lineal estandar que siempre se realiza tanto si se especifica o no la
combinacin de carga P-. Este parmetro limita el nmero de iteraciones adicionales realizadas
para el efecto de la combinacin de carga P-. Si el nmero de iteraciones mxima es cero,
entonces el sistema ser igual a uno en donde no se especifica combinacin de carga P-. (ver
figura 4.1.1.1)

Fallo de convergencia
Si la convergencia no se obtiene despus de que se hayan efectuado un nmero mximo
de iteraciones, los resultados del anlisis pueden carecer de sentido, este fallo se puede deber a
varias causas:
Fueron permitidas muy pocas iteraciones. Usualmente un nmero razonable es de 2 a 5,
aunque se pueden requerir mas, dependiendo de cada problema en particular.
Se ha usado una tolerancia de convergencia muy pequea o una muy grande. En
cualquiera de los dos casos se obtendran resultados sin sentido.
La estructura esta cerca de la carga crtica de pandeo. Esto se puede impedir arriostrando
la estructura para impedir el pandeo o disminuyendo la carga axial aplicada. (Computers
and Structures, 1998)
49

4.1.4. Anlisis inicial P-

Para muchas aplicaciones, es adecuado considerar el efecto P- sobre la estructura bajo
un sistema de cargas (usualmente de gravedad), y considerar todos los otros anlisis como
lineales usando la Matriz de rigidez desarrollada para este sistema de cargas P-. Esto permite
que todos los resultados de los anlisis sean superpuestos para efectos de diseo.
Para hacer esto, se define un caso de anlisis no lineal, que tiene a lo menos las siguientes
caractersticas:
Un sistema llamado PDELTA
Las condiciones iniciales comienzan desde cero
Se aplica el caso de carga que causar el efecto P-; a menudo se tratar de carga muerta
y una fraccin de la carga viva
Para la no linealidad geomtrica, escoger efectos P-

Figura 4.1.4.1. Anlisis no lineal esttico

Nos referiremos a este caso esttico no lineal como el caso inicial P-. Se puede definir o
modificar otros casos de anlisis lineales de modo que usen la rigidez del caso PDELTA.
(Computers and structures, 2007)
El efecto P-delta escencialmente modifica las caractersticas de la estructura afectando los
resultados de todos los otros anlisis realizados por ejemplo las cargas estticas lineales como
50

peso propio y sobrecarga. Debido a que se utiliza una matriz de rigidez constante para todos esos
anlisis lineales, los resultados se pueden superponer, lo que resulta eficiente desde el punto de
vista computacional. (Computers and structures, 1998)

4.2. PROGES
4.2.1. Introduccin

En este apartado se muestra el proceso de clculo que sigue este programa y los
fundamentos tericos en los cuales est basado para realizar el anlisis estructural.

Figura 4.2.1.1. Diagrama de flujo de programa Proges
Fuente: Lopetegui, 2006
51

AN = _
N
nuco
- N
c]o
N
1uu] % ec. 4.2.1.1
PROGES tiene la ventaja de ser un programa abierto, debido a que esta escrito en lenguaje
Fortran, permite realizar clculo esttico y dinmico de estructuras de barra, incluyendo efectos
de segundo orden, que es lo que interesa para esta tesis, no es tan conocido y usado como
SAP2000, pero entrega resultados de una gran variedad de problemas con una buena precisin.
En la figura 4.2.1.1 se muestra el diagrama de flujo que se utiliza en este programa.
Para modelar una estructura se puede usar un programa llamado Datos.exe, en el cual se
ingresan las coordenadas de la estructura, las barras, las caractersticas de las barras, es decir
inercias, areas, radios de giro, etc. se ingresan tambin las cargas en los nudos y barras de
acuerdo a los cdigos indicados en el manual de uso, luego este programa genera un archivo de
datos (.dat) que Proges leer, entregando los resultados en un archivo de salida con extencin
.dat. Los archivos de entrada y salida de datos se muestran en el anexo C de esta memoria.

4.2.2. Fundamentos tericos de PROGES

El mtodo que usa Proges consiste en formar la matriz de rigidez para el elemento, a
partir de las Funciones de Estabilidad (ver Anexo A) obtenidas mediante la resolucin de la
ecuacin diferencial del equilibrio (ec. 4.2.2.9), considerando los efectos de segundo orden
causados por la inestabilidad elstica del elemento. Para esto se definen los valores de F
i
con i
variando de 1 hasta 7 los cuales corresponden a los valores de la matriz de rigidez considerando
la fuerza axial a travs del trmino , que es a su vez funcin de la fuerza axial, y que en los
clculo lineales no se considera. Proges, a diferencia del mtodo de la Matriz Geomtrica, trabaja
en base a la fuerza axial, esto quiere decir que en cada paso iterativo se va incrementando la
fuerza axial sobre los elementos y se finalizan las iteraciones cuando se cumple la ec. 4.2.1.1

Figura 4.2.2.1. Ecuaciones de estabilidad para un segmento de columna
Fuente: Liew et al, 1999
52

A continuacin se deduce la ecuacin diferencial de equilibrio para el elemento
diferencial mostrado en la figura 4.2.2.1
Haciendo momento respecto a b se obtiene lo siguiente
Jx +PJy +H- _H+
JH
Jx
Jx] = u cc. 4.2.2.1
=
JH
Jx
-P
Jy
Jx
cc. 4.2.2.2

A partir de la sumatoria de fuerzas horizontales se tiene
- +_ +
J
Jx
Jx] = u cc. 4.2.2.S
J
Jx
= u cc. 4.2.2.4

Diferenciando la ecuacin 4.2.3 respecto a x se tiene

J
Jx
=
J
2
H
Jx
2
-P
J
2
y
Jx
2
cc. 4.2.2.S
Utilizando la ec. 4.2.5 se tiene que
J
2
H
Jx
2
-P
J
2
y
Jx
2
= u cc. 4.2.2.6

Si se reemplaza M por el siguiente valor se llega a la ec. 4.2.2.9
H = -EI
J
2
y
Jx
2
cc. 4.2.2.7
J
4
y
Jx
4
+
P
EI
J
2
y
Jx
2
= u cc. 4.2.2.8
y
Iv
-
N
0
y
II
EI
= u ec. 4.2.2.9
e = l
_
|N
0
|
EI
(iguol o kI cn cl oncxo B) cc. 4.2.2.1u
(Liew et al, 1999)
La solucin de la ecuacin diferencial ec. 4.2.2.9 es la que se muestra a continuacin:
Compresin:
N
0
= -P

y = Ascn e
x
l
+Bcos e
x
l
+C
x
l
+B
Traccin:
N
0
= P

y = Ascnb e
x
l
+Bcosb e
x
l
+C
x
l
+B
53

Los valores de A B C y D se obtienen de acuerdo a las condiciones de borde de cada
problema en particular, estas constantes y la forma en que se llega a los valores de F
1
y F
2
dados
en la tabla 4.2.2.1, que corresponden a las funciones de estabilidad s
ii
y s
ij
respectivamente, se
muestran en detalle en el anexo A de esta memoria.

Tabla 4.2.2.1. Valores de F
para la matriz de rigidez

Fuente: Lopetegui, 2006

En la tabla 4.2.2.1 se muestran los valores de F
para el caso de carga axial de compresin,

de traccin y cuando no hay carga axial.
Se definen a continuacin las matrices de rigidez modificadas por los factores F
i
, los
cuales incorporan la fuerza axial para los elementos con distintas condiciones de borde.

Elemento empotrado en ambos extremos

l
l
l
l
l
l
H
l
c
l
c
H
K
N
K
l
c
K
l
c
1
1
1
1
1
1
=
l
l
l
l
l
l
l
F
1
u F
3
I F
2
u -F
3
I
u
F
3
I
F
2
u
-F
3
I
F
P
l
c
2
u u -F
P
l
c
2
u
u F
4
I
2
F
3
I u -F
4
I
2
u u F
1
u -F
3
I
-F
P
l
c
2
F
3
I u
u -F
4
I
2
-F
3

F
P
l
c
2
u
u F
4
I
2
1
1
1
1
1
1
1

l
l
l
l
l
l
l
-
u
-
:
k
-
u
k
-
:
k
-
1
1
1
1
1
1
1

54

Elemento empotrado-articulado

l
l
l
l
l
l
H
l
c
l
c
H
K
N
K
l
c
K
l
c
1
1
1
1
1
1
=
l
l
l
l
l
l
l
F
5
u F
5
I u u -F
5
I
u
F
5
I
u
u
-F
5
I
F
P
l
c
2
u u -F
P
l
c
2
u
u F
6
I
2
u u -F
6
I
2
u u u u u
-F
P
l
c
2
u u
u -F
6
I
2
u

F
P
l
c
2
u
u F
6
I
2
1
1
1
1
1
1
1

l
l
l
l
l
l
l
-
u
-
:
k
-
u
k
-
:
k
-
1
1
1
1
1
1
1

Elemento articulado-empotrado

l
l
l
l
l
l
H
l
c
l
c
H
K
N
K
l
c
K
l
c
1
1
1
1
1
1
=
l
l
l
l
l
l
l
u u u u u u
u
u
u
u
u
F
P
l
c
2
u u -F
P
l
c
2
u
u F
6
I
2
F
5
I u -F
6
I
2
u F
5
I F
5
u -F
5
I
-F
P
l
c
2
u u
u -F
6
I
2
-F
5
I

F
P
l
c
2
u
u F
6
I
2
1
1
1
1
1
1
1

l
l
l
l
l
l
l
-
u
-
:
k
-
u
k
-
:
k
-
1
1
1
1
1
1
1

Elemento articulado en ambos extremos

l
l
l
l
l
l
H
l
c
l
c
H
K
N
K
l
c
K
l
c
1
1
1
1
1
1
=
l
l
l
l
l
l
u
u
u
u
u
u
u
F
P
l
c
2
u
u
-F
P
l
c
2
u
u
u
-F
7
I
2
u
u
F
7
I
2

u
u
u
u
u
u
u
-F
P
l
c
2
u
u
F
P
l
c
2
u
u
u
F
7
I
2
u
u
-F
7
I
2
1
1
1
1
1
1

l
l
l
l
l
l
l
-
u
-
:
k
-
u
k
-
:
k
-
1
1
1
1
1
1
1

55

Los elementos F de las matrices corresponden a los siguientes valores:

F
= F
i
.
EI
l
.
l
c
E
c
I
c
F
P
=
EA
l
.
l
c
E
c
I
c
I =
l
c
l

-
=
E
c
I
c
l
c
u
-
= u
E
c
I
c
l
c
2
:
-
= :
E
c
I
c
l
c
2

Donde:
l
c
, E
c
, I
c
, son valores medios aproximados de los largos de las barras, mdulo de
elasticidad y momento de inercia.
i : 1, 2, , 7

Se puede notar en la tabla 4.2.2.1 que para valores de = 0 la matriz es la Matriz de
rigidez lineal que se utiliza para clculos de primer orden.
Al realizar los clculos el programa converge rapidamente si la carga axial a la que est
siendo sometida la estructura es proporcional a la deformacin, es decir, si se trabaja en el rango
lineal. Para valores mayores de P se requiere de mas iteraciones para realizar los clculos, hasta
que el programa no puede continuar, por que las deformaciones son muy grandes y no hay
covergencia. Esto sucede cuando se esta cerca del valor crtico, si se alcanza este valor la matriz
se vuelve singular, y mas alla de este valor el sistema se vuelve inestable, y los resultados
entregados por el programa van a ser incorretos. (Lopetegui, 2006)

Captulo V Estructura tipo marco plano analizada
56

Captulo V
Estructura tipo marco plano analizada

Se muestra en la figura 5.1 el marco plano que ser analizado bajo las condiciones de
carga y apoyo indicadas, por los tres mtodos propuestos los cuales son:

Mtodo 1: Anlisis de Segundo Orden utilizando la Matriz Gemetrica
Mtodo 2: Anlisis de Segundo Orden utilizando Proges
Mtodo 3: Anlisis de Segundo Orden utilizando SAP2000

Antes de realizar los anlisis se establecen ciertas hiptesis las cuales son las siguientes:

El elemento que se analiza debe ser esbelto y de seccin constante doblemente simtrica
No se considera la plasticidad del material
La columna esta debidamente arriostrada para que no se produzca pandeo lateral ni flexo-
torsional.
La columna esta debidamente atiesada para que no se produzca pandeo local

Las hiptesis permiten realizar el anlisis de la estructura, considerando que el P
cr
de la
estructura ser el que se produzca por pandeo general o de Euler alrededor de su eje fuerte y no se
producir ningna otra fuerza crtica antes de que se alcance esta carga crtica.

Figura 5.1. Marco analizado por los tres mtodos

57

El marco de la figura 5.1 se encuentra en el captulo 11 (pgina 684) del libro llamado
Mechanics of structures: Variational and Computational Methods (ver referencia), el cual fue
analizado en dicho libro por el mtodo de la Matriz Geomtrica para considerar los efectos de
segundo orden, se realizaron dos iteraciones y los resultados obtenidos se pueden ver en la tabla
5.3.1. Los clculos realizados por el primer mtodo se encuentran en el anexo B de esta memoria.

5.1. Modelo de la estructura en PROGES

Se model luego la estructura en PROGES discretizando cada elemento en cuatro partes,
como se muestra en la figura 5.1.1

Figura 5.1.1. Modelo de marco en PROGES

Los archivos de entrada y salida de datos para el modelo de la figura 5.1.1 en PROGES
se encuentran en el anexo C de esta memoria.

5.2. Modelo de la estructura en SAP2000

Se parte especificando las unidades de medida en las que se trabaja, las cuales son KN-m,
luego se define el tipo de material y las secciones que se usarn para las dos columnas y las vigas
que componen el marco.

5.2.1. Casos de carga

Una vez que se ha dibujado ste, se definen tres casos de carga:
PP: En el cual se encuentran las cargas verticales
SC: En el cual se encuentran las cargas horizontales
Pdelta: PP + 0.25SC (las cargas no se especifican directamente como en los casos
anteriores, sino que se definen en el caso de anlisis, ver figura 5.2.2.3)
58

Figura 5.2.1.1. Casos de carga en SAP2000
Los distintos casos de carga se superponen a travs de las combinaciones que se
especifican mas adelante en este captulo.

5.2.2. Casos de anlisis

El anlisis para el caso de carga Pdelta es de tipo no lineal escogiendo como parmetro
geomtrico no lineal P-delta (ver figura 5.2.2.3 y 5.2.2.4), y para PP y SC el anlisis del tipo
lineal esttico, y son realizados una vez que termina el caso de anlisis Pdelta (ver figuras
5.2.2.1 y 5.2.2.2), es decir, utilizando la matriz de rigidez que se obtiene a partir del anlisis no
lineal realizado, por lo que en los tres casos de carga se esta considerando la influencia de la
carga axial del caso PP, sobre las deformaciones y esfuerzos internos generados en la estructura.

Figura 5.2.2.1. Caso de anlisis de PP en SAP2000
59

Figura 5.2.2.2. Caso de anlisis de SC en SAP2000

Figura 5.2.2.3. Caso de anlisis de Pdelta en SAP2000

60

Figura 5.2.2.4 Especificacin de efecto P-Delta

5.2.3. Combinaciones

Se definen dos combinaciones de carga:

C1: PP + SC + Pdelta
Se superponen los tres casos de carga definidos, por lo que la deformacin
resultante ser la suma de la deformacin producida por la carga lateral (SC), y cargas
axiales (PP) y la deformacin producto del clculo iterativo de la carga axial (PP)
actuando sobre el desplazamiento que se produce por la carga lateral

C2: PP + SC
En esta combinacin se superponen solo dos casos de carga: las axiales (PP) y
las laterales (SC), sin embargo, el que no se incluya el caso de carga Pdelta, no significa
que no se este tomando en cuenta la no linealidad geomtrica de la estructura, no hay que
olvidar que estos casos parten luego del anlisis Pdelta por lo que la rigidez del marco es
menor, y esto se ve reflejado en la influencia que tiene la carga axial en la matriz de
rigidez reformulada por el anlisis no lineal Pdelta. La idea de estudiar este caso de carga,
es saber si la matriz de rigidez que se forma es similar a alguna de las formadas por los
otros dos mtodos, a travs de la comparacin de resultados, cabe destacar tambin que
esta matriz no se ve influenciada por el factor 0,25 (ver figura 5.2.2.3 scale factor) ya
que al ir cambiar este valor no hubo mayor variacin entre los resultados entregados por
SAP2000 C2.

61

Otra consideracin importante es que en el anlisis atravs de SAP2000 no se toma en
cuenta el peso de los elementos (ver figura 5.2.3.1), ya que este programa calcula el peso propio
automticamente incorporando nuevas fuerzas a la estructura, este hecho puede hacer que los
resultados sean diferentes a los entregados por el primer y segundo mtodo (ver nota 1), por esta
razn el peso de los elementos se iguala a cero en SAP2000. La idea es que las condiciones del
modelo sean semejantes para los tres mtodos, y de esta forma poder evaluar de manera correcta
las diferencias que pudieran existir.

Nota1: Se verifica la situacin anterior de la siguiente manera: Se realiza un anlisis elstico lineal, para
la estructura mostrada en la figura 5.1, considerando la masa (mass per unit volume) = 7.849kg/m
3
y
peso (weiht per unit volume) = 76.9729 kN/m
3
(figura 5.2.6) y el desplazamiento que se obtiene en el
nudo 3 fue de 6.78 cm, definiendo como cero estos dos valores se obtuvo un desplazamiento de 7.85 cm,
que es mayor y se acerca bastante al obtenido por un anlisis lineal realizado.

Figura 5.2.3.1. Datos de las propiedades del material

Finalmente se realiza el anlisis de la estructura, luego se revisan los valores de la
deformacin en el nudo tres y las fuerzas internas obtenidas para cada elemento, tanto para la
combinacin C1 como para la C2.
Todo lo dicho anteriormente indica que el programa SAP2000 llevar a cabo el anlisis de
la siguiente forma:

1. Comienza realizando el anlisis no lineal del caso de carga Pdelta especificado como
PP+0.25SC (de acuerdo al cpitulo 4.1 de esta memoria) y considerando como parmetro
62

no lineal el efecto P-delta, en esta etapa se realizan las iteraciones correspondientes
hasta que se obtiene la convergencia especificada.
2. En cada paso iterativo del anlisis Pdelta la matriz de rigidez se reformula hasta que se
obtiene la convergencia, y se forma la matriz de rigidez final, la cual considera la prdida
de rigidez de la estructura por el efecto de la carga axial aplicada, que corresponde a todas
las cargas verticales sobre sta.
3. A continuacin realiza los anlisis de tipo lineal esttico correspondientes al PP y SC,
pero usando la matriz de rigidez final que se forma al finalizar el caso de anlisis Pdelta.

El archivo de salida de datos que entrega SAP2000, en donde se especifica cada uno de
los pasos que se llevan a cabo y la iteraciones realizadas en cada uno de ellos, se encuentra en el
anexo C de esta memoria (ver archivos C.7, C.8, C.9).

5.3. Resultados

En la tabla 5.3.1 se muestran los resutados de fuerzas y desplazamientos obtenidos a partir
de los tres mtodos de anlisis, los cuales se comparan con un anlisis lineal de primer orden,
para la estructura con las condiciones de carga mostradas en la figura 5.1.

Mtodo 1
M. Geom.
Mtodo 2
Proges
Mtodo 3
SAP2000 C2
Mtodo 3
SAP2000 C1
1 Orden
Lineal
1 U
2
(cm) -15.423 -14.21 -15.61 -17.91 -7.831
2
2
0.01211 0.0355 0.013
0.013 0.00553
3 M
(1,1)
(KNm) 1055.30 1021.83 1066.83 1246.37 584.20
4 M
(1,2)
(KNm) -705.40 -679.10 -712.52 -844.41 -363.80
5 V
(1,1)
(KN) -233.90 -228.28 -273.1 -283.12 -154.50
6 V(
1,2)
(KN) -161.90 -156.28 -165.1 -193.11 -82.50
7 V
(2,2)
(KN) 100.50 97.91 101.25 144.44 66.40
8 V
(2,3)
(KN) 40.50 37.91 41.25 24.44 6.40
9 V
(3,3)
(KN) 79.40 73.28 82.60 90.01 0.00
10 N
1
(KN) -1150.50 -1147.91 -1151.25 -2244.44 -1116.4
11 N
2
(KN) -79.40 -73.28 -82.60 -89.99 0.00
12 N
3
(KN) -2059.50 -2062.09 -2058.75 -4175.56 -2093.60
Tabla 5.3.1. Resultados de desplazamientos y fuerzas internas de los elementos.

63

Donde:
U
2
= desplazamiento producido en el nudo 2 (ver figura 5.1)
2
= ngulo que se forma en el nudo 2
M, V y N = estn definidos como se muestra en la figura 5.3.1

Figura 5.3.1. Fuerzas internas en cada barra

A continuacin se determina la curva de carga-deformacin, para poder visualizar con
mayor claridad las diferencias entre los tres mtodos propuestos.

Cargas puntuales (KN) 1 orden Mat. Geom SAP2000 C2 SAP2000 C1 Proges
N
2x
N
2z
N
3
(cm) (cm) (cm) (cm) (cm)
050100 1,76 1,9 1,9 2,01 1,81
5100200 2,10 2,33 2,34 2,54 2,22
10200400 2,45 2,85 2,86 3,17 2,71
15300600 2,79 3,43 3,44 3,86 3,25
20400200 3,14 4,08 4,1 4,62 3,85
255001000 3,49 4,81 4,83 5,48 4,53
306001200 3,83 5,63 5,65 6,44 5,28
357001400 4,18 6,56 6,6 7,52 6,13
408001600 4,53 7,63 7,68 8,76 7,10
459001800 4,87 8,87 8,94 10,18 8,21
5010002000 5,22 10,33 10,43 11,85 9,51
5511002200 5,57 12,06 12,2 13,82 11,03
6012002400 5,91 14,16 14,35 16,2 12,83
6513002600 6,26 16,47 17,03 19,13 15,02
7014002800 6,60 19,99 20,43 22,82 17,72
7515003000 6,95 24,2 24,92 27,66 21,13
8016003200 7,23 29,86 31,1 34,26 25,56
8517003400 7,64 37,81 40,15 43,82 31,53
9018003600 7,99 49,64 54,67 59,06 39,96
Tabla 5.3.2. Valores de deformacin vs carga aplicada en nudos 2 y 3
64

N
2x
(KN) : Carga puntual sobre el nudo 2 en direccin x (ver figura 5.1)
N
2z
(KN) : Carga puntual sobre el nudo 2 en direccin z (ver figura 5.1)
N
3
(KN) : Carga puntual sobre el nudo 3 en direccin z (ver figura 5.1)
(cm) :Desplazamientoenelnudo3(ver figura 5.1)

Esto se llev a cabo haciendo variar las cargas N
2x
,

N
2z
y N
3
en la estructura y evaluando
los desplazamientos en el nudo 3. Los resultados obtenidos por cada mtodo se muestran en la
tabla 5.3.2. Luego se grafican estos valores para confeccionar las curvas que se muestran en la
grfico 5.3.1

Grfico 5.3.1. Curvas de carga vs deformacin

A travs de este grfico se puede ver que ya desde la carga N
2x
= 45 KN se comienza a
presentar una gran diferencia entre los desplazamientos entregados por cada mtodo, en
comparacin con un anlisis lineal, cabe destacar tambin la semenjanza entre la curva generada
por SAP2000 C2 y el Mtodo de Matriz Geomtrica, los valores de desplazamientos y esfuerzos
obtenidos con cada mtodo son casi iguales. Las conclusiones acerca de los anteriormente
expuesto se ven en el captulo 6 de esta memoria.

0
5
10
15
20
25
30
35
40
45
50
55
60
65
70
75
80
85
90
95
0 5 10 15 20 25 30 35 40 45 50 55 60
Anlisislineal
Matriz
Geomtrica
SAP2000C1
SAP2000C2
Proges
(cm)
Carga N
2x
(KN)
65

Grfico 5.3.2. Variacin porcentual de los desplazamientos obtenidos por los tres mtodos con
respecto un anlisis lineal.
Fuente: Elaboracin propia.

Grfico 5.3.3. Variacin porcentual de desplazamientos de MG con respecto a Proges

0
100
200
300
400
500
600
700
0 10 20 30 40 50 60 70 80 90 100
%
N
2x
(KN)
Variacin porcentual de los desplazamientos obtenidos por
los tres mtodos con respecto a un anlisis lineal
Matriz
Geomtrica
SAP2000C2
SAP2000C1
PROGES
0
2
4
6
8
10
12
14
16
18
20
22
24
26
0 10 20 30 40 50 60 70 80 90 100
%
N
2x
(KN)
Variacion porcentual de desplazamientos de
MG con respecto a PROGES
Variacin de MG
vs PROGES
66

Por otra parte, como se vi en el captulo 3, hay ciertos paramtros definidos tanto en el
AISC como en el Eurocdigo, B
2
y F
A
respectivamente, que establecen cuando es necesario
realizar un anlisis de segundo orden de la estructura, en la tabla 5.3.3 y 5.3.4 se tiene el resultado
de estos valores para el marco de la figura 5.1. El procedimiento de clculo para obtener los
valores B
2
y F
A
se muestra en el anexo D, las frmulas usadas para cada caso son las siguientes:

Para el AISC

B
2
=
1
1 -
o P
nt
P
c2
1 con P
c2
= R
M
EI
A
H

Donde:

P
nt
= Carga vertical total soportada por el piso, usando combinacin de cargas (N)
R
M
= 0,85 para marcos no arriostrados
H = Corte de piso producido por la fuerza lateral usada para calcular delta H (N)
L = Altura de piso (mm)
H
= Deformacin de primer orden entrepiso, debido a fuerzas laterales (mm)
= 1.6 factor de seguridad

SC N
2x
(N) P
nt
(N) H (N)
H
(mm) P
e2
B
2
<1,5
5000 360000 77000 26,8 19537313,4 1,03
10000 660000 82000 30,4 18342105,3 1,06
15000 960000 87000 34,1 17348973,6 1,09
20000 1260000 92000 37,7 16594164,5 1,13
25000 1560000 97000 41,4 15932367,1 1,18
30000 1860000 102000 45 15413333,3 1,23
35000 2160000 107000 48,6 14971193,4 1,30
40000 2460000 112000 52,3 14562141,5 1,37
45000 2760000 117000 55,9 14232558,1 1,44
50000 3060000 122000 59,5 13942857,1 1,54
55000 3360000 127000 63,2 13664557 1,64
60000 3660000 132000 66,8 13437125,7 1,77
65000 3960000 137000 70,5 13214184,4 1,92
70000 4260000 142000 74,1 13031039,1 2,09
75000 4560000 147000 77,7 12864864,9 2,31
80000 4860000 152000 81,4 12697788,7 2,57
85000 5160000 157000 85 12560000 2,91
90000 5460000 162000 88,6 12433408,6 3,36
Tabla 5.3.3. Clculo del factor B
2

67

Para el Eurocodigo
F
A
=
1
[1 -
I
sd
I
c
con o
c
=
I
sd
I
c
o
c
=
b
x
I
E
u.1
Donde:
= Desplazamiento horizontal entre piso producido por las fuerzas laterales mas las fuerzas
ficticias
H = Reaccin horizontal total en el piso producido por las fuerzas laterales y fuerzas ficticias
F
A
= Factor de amplificacin (semejante a B
2
)
V
sd
V
cr
e = Valor elstico crtico de la carga vertical total considerando un estado no arriostrado
= Imperfeccin inicial aplicada a la carga vertical total (ver ec. 3.3.1.1)
SC N
2x
= Sobrecarga lateral puntual en el nudo 2 direccin x (ver figura 5.1)

SC N
2x
(N) V (N) V(N) SC+V (N) H (N) (mm) cr<0,25 F
A

5000 360000 1101,6 6101,6 78180 27,7 0,01 1,01
10000 660000 2019,6 12019,6 84100 32 0,03 1,03
15000 960000 2937,6 17937,6 90020 36,3 0,04 1,05
20000 1260000 3855,6 23855,6 95900 40,6 0,06 1,07
25000 1560000 4773,6 29773,6 101770 44,8 0,08 1,09
30000 1860000 5691,6 35691,6 107700 49,1 0,10 1,11
35000 2160000 6609,6 41609,6 113610 53,4 0,12 1,14
40000 2460000 7527,6 47527,6 119530 57,7 0,14 1,17
45000 2760000 8445,6 53445,6 125450 62,1 0,17 1,20
50000 3060000 9363,6 59363,6 131360 66,4 0,19 1,23
55000 3360000 10281,6 65281,6 137280 70,7 0,21 1,27
60000 3660000 11199,6 71199,6 143200 75 0,23 1,31
65000 3960000 12117,6 77117,6 149120 79,3 0,26 1,35
70000 4260000 13035,6 83035,6 155040 83,6 0,28 1,40
75000 4560000 13953,6 88953,6 160950 87,9 0,31 1,45
80000 4860000 14871,6 94871,6 166870 92,2 0,33 1,50
85000 5160000 15789,6 100789,6 172790 96,5 0,36 1,56
90000 5460000 16707,6 106707,6 178710 100,8 0,38 1,62
Tabla 5.3.4. Clculo del factor
cr


En la tabla 5.3.3 se muesta que a partir de P
nt
= 1860 KN se debe realizar un anlisis
aproximado de segundo orden con el mtodo de amplificacin de momentos de acuerdo al
cdigo AISC, ste se puede utilizar solo hasta que P
nt
= 3060 KN ya que a partir de este valor
B
2
es mayor que 1.5 y se exige un anlisis de segundo orden mas exacto, incorporando la rigidez
reducida, imperfecciones iniciales, etc, como se vi en la seccin 3.2.3, mientras que el
68

Eurocdigo exige un anlisis de segundo orden cuando V

= 2160 KN y se puede utilizar el
mtodo de clculo aproximado que se propone hasta que V

= 3960 KN, en donde el valor de
cr

es mayor a 0.25 y se debe efectuar un anlisis de segundo orden mas exacto (ver tabla 5.3.5).
Se grafican los valores de B
2
y F
A
contra la carga N
2x
aplicada, y se pueden visualizar las
diferencias que existen entre los cdigos de diseo en el grfico 5.3.4.

Grfico 5.3.4. Curva factores de amplificacin (B
2
y F
A
) vs carga

Con respecto a la prctica en Chile, lo usual es limitar las deformaciones de las estructuras
ante las cargas laterales, como se menciona en Nch 433 o Nch 2369, a travs de este parmetro se
determina cuando es necesario llevar a cabo un anlisis de segundo orden de la estructura, si se
supone que las cargas laterales son de viento, la deformacin mxima permitida es de 0.004H que
para la estructura de la figura 5.1 es de 3.2 cm, de partida la deformacin que se genera en este
marco es mucho mayor a la mxima permitida, sta se obtiene a partir de un anlisis lineal de
primer orden de la estructura y es de 7.83 cm (ver tabla 5.3.1) por lo que hay que considerar
el efecto P-delta para este caso. Ahora bien, si se quita la carga puntual horizontal de 82.5 kN en
el nudo 2 la deformacin obtenida por un anlisis elstico de primer orden es de 1.84 cm lo que
indica que no es necesario hacer un anlisis de segundo orden, pero al realizar dicho anlisis la
deformacin en el nudo 2 es de 3.67 cm que supera la mxima permitida por norma.
Con esto se puede ver que no es tan simple tomar una decisin respecto a si considerar o
no el efecto de la carga axial en una estructura, solo a travs de las deformaciones, se demuestra
que ante una carga axial considerable se pueden tener problemas de inestabilidad, que a travs de
un anlisis de primer orden no se alcancen a percibir, dependiendo del tipo de marco, aun cuando
la carga lateral sea pequea y un anlisis de primer orden entregue resultados de deformaciones
menores a los mximos permitidos.
0
10000
20000
30000
40000
50000
60000
70000
80000
90000
100000
0 0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 1,8 2 2,2 2,4 2,6 2,8 3 3,2 3,4 3,6
AISC
EC3
B
2
y F
A
SC N
2x
(N)
Captulo VI Comentarios y conclusiones
69

Captulo VI
Comentarios y conclusiones

6.1. Conclusiones generales

Luego del anlisis realizado a la estructura de la figura 5.1 por los tres mtodos de clculo
propuestos se concluye lo siguiente:

6.1.1 Respecto a los mtodos estudiados

La deformacin que se obtiene en el nudo tres a partir de un anlisis elstico de primer
orden del marco, es considerablemente menor a la que se obtiene por un anlisis elstico de
segundo orden (ver tabla 5.3.1), por consiguiente los resultados de los esfuerzos tambin van a
ser mayores, esto se debe a que la fuerza axial en la estructura es muy grande y un anlisis de
primer orden en esta situacin no entrega resultados completamente acertados.

Los resultados que se tienen para el mtodo 1 (Matriz Geomtrica) y el mtodo 3 (SAP2000) con
la combinacin de cargas C2 (ver seccin 5.2.3 combinaciones de carga) son similares, de igual
manera la curva de carga vs deformacin para estos dos mtodos siguen la misma trayectoria
(grfico 5.3.1), se puede decir a partir de esto, que SAP2000 genera una matriz de rigidez final,
luego de realizar el proceso de iteraciones (en el caso de anlisis no lineal Pdelta), similar a la
matriz que se obtiene utilizando las aproximaciones polinmicas de las funciones de estabilidad,
para considerar los efectos de la geometra deformada en la rigidez de la estructura, no es el caso
de el mtodo utilizado en PROGES, aqu los desplazamientos obtenidos para cada valor de carga
son inferiores a los obtenidos por los otros dos mtodos, el grfico 5.3.3 muestra la variacin
porcentual entre las deformaciones que se obtienen a partir del mtodo de la Matriz geomtrica y
Proges, esta se incrementa a medida que aumenta la carga sobre la estructura, la diferencia
fundamental est en que PROGES utiliza las funciones de estabilidad trigonmetricas e
hiperblicas para valores de > 0.1 (ver tabla 4.2.1), y utiliza las aproximaciones polinmicas
para valores de < 0.1, es decir, la geometra deformada de los elementos, que se toma en cuenta
en este programa, es de mayor exactitud, por otra parte PROGES realiza la iteracin en forma
incremental variando la fuerza axial. Se concluye que para cargas muy grandes (en este caso a
partir de 60 KN en la tabla 5.3.2) las aproximaciones polinmicas para las funciones de
estabilidad subestiman en cierta medida la capacidad de la estructura, dando como resultado
mayores deformaciones.

70

En la combinacin C2 no se incluye el caso de carga Pdelta, esto no quiere decir que en los
clculos no este siendo considerada la geometra deformada, porque como se dijo anteriormente
los casos de anlisis SC y PP, pese a ser lineales, comienzan luego de que se han realizado las
iteraciones del caso de anlisis Pdelta, es decir sobre la geometra deformada de la estructura (ver
figuras 5.2.2.2 y 5.2.2.3).

Con respecto a la combinacin C1= PP+SC+Pdelta en SAP2000, se puede ver que los resultados
son evidentemente mayores a los entregados por los otros dos mtodos, no tan solo en
deformacin, sino que tambin en las fuerzas internas producidas en los elementos, llama la
atencin que las fuerzas axiales para los elementos 1 y 3 son considerablemente mayores a las
que se obtienen en las otras situaciones, ver tabla 5.3.1

A travs del grfico 5.3.1 se puede ver que las deformaciones comienzan a presentar mayores
diferencias a partir N
2x
= 30 KN, en el grfico 5.3.2 se muestra esta variacin en forma de
pocentaje, y para este valor de carga promedia alrededor de un 50,13% por sobre los resultados
de un anlisis lineal, por supuesto las deformaciones, luego de este valor, comienzan a ser cada
vez mayores, desviandose de la recta que representa el analisis lineal, para dar forma a una
parbola; estas deformaciones van a afectar el comportamiento de la estructura, es decir que a
partir de este valor ya se hace necesario realizar un anlisis de segundo orden de la estructura.
Esto coincide con lo que se establece en AISC mediante el factor B
2
(tabla 5.3.3). Por otra parte
el Eurocdigo, mediante el factor
cr
, indica que es necesario un anlisis de segundo orden
cuando N
2x
= 35 KN como se muestra en la tabla 5.3.4

Con respecto a las iteraciones realizadas, en el proceso de anlisis de la estructura, a travs de
PROGES y SAP2000, se pudo notar que a medida que aumentaba la carga axial el nmero de
iteraciones requerido para que el sistema llegara a la convergencia era mayor, por ejemplo en el
caso de Proges se necesitaron 5 iteraciones al aplicar la carga de 1050 KN y 2100 KN en los
nudos 2 y 3 repectivamente (ver archivo C.2) mientras que cuando se aplic una carga de 1800
kN y 3600 kN en los nudos 2 y 3, el nmero de iteraciones fue 7 (ver archivo C.6) y para cargas
mayores se produca un error en el programa es decir no haba convergencia. Al comparar cada
paso iterativo para cargas mayores a 1700 kN (ver archivos C.3 y C.4) se pudo ver que los
desplazamientos entre la ltima y penltima iteracin no eran iguales, lo que indica que no haba
convergencia, esto sucede porque al comenzar a llegar a la carga crtica de pandeo se tienen
grandes diferencias de deformacin ante pequeos incrementos en la fuerza axial y como
PROGES va comparando las fuerzas axiales entre una iteracin y otra el sistema converge, aun
cuando los desplazamientos sean muy distintos, esto indica que se debe aplicar un factor de carga
cuando la deformacin de la estructura se hace mayor.
71

Algo similar ocurre en SAP2000 con respecto a las iteraciones; en el anlisis no lineal P-delta,
este programa realiza una cierta cantidad de pasos, dentro de los cuales se realizan iteraciones
que convergen de acuerdo a lo que el usuario ha especificado previamente (criterio de
convergencia, ver capitulo 4.1), para 1050kN y 2100 kN en los nudos 2 y 3 repectivamente se
tienen dos iteraciones en cada paso (ver archivo C.7), mientras que para 1800 kN y 3600 kN el
nmero de iteraciones desde el paso uno al diez va en aumento, por ejemplo en un paso se tiene
hasta 7 iteraciones (ver archivo C.8), si las cargas son mayores ya no hay convergencia en
algunos pasos, para el nmero limitado de iteraciones que ha sido especificado previamente por
el usuario (ver figura 4.1.1), en esta situacin SAP2000 disminuye los incrementos a la mitad
para alcanzar la convergencia (ver archivo C.9).

Mediante el mtodo de la Matriz Geomtrica se calcula el factor , que amplifica todas las cargas
de la estructura de la figura 5.1, para que esta se vuelva inestable, el clculo se puede ver en
detalle en el anexo B de esta memoria, el procedimiento se aplica a la estructura para las distintas
cargas que se muestran en la tabla 5.3.2 y as obtener el grfico B.1, en el cual se puede apreciar
que dichos factores son menores y la diferencia entre uno y el consecutivo es tambin menor, a
medida que la deformacin aumenta y a medida de que la carga aplicada se va acercando a la
crtica, esto comprueba que mientras mas cerca se este de este valor crtico, la diferencia entre los
desplazamientos crece ante pequeos incrementos de la fuerza axial.

6.1.2. Respecto a las normas estudiadas

Existen diferencias claras entre las dos normas estudiadas con respecto al tema de anlisis.
En los mtodos usados por el AISC no se toma en cuenta la imperfeccin del prtico traslacional
ni de los elementos, por lo que se debe emplear la longitud de pandeo de los elementos al
verificarlos utilizando las frmulas de interaccin dadas en esta especificacin (seccin 3.2), en
cambio el EC3 permite emplear la longitud del elemento (sin los coeficientes de longitud
efectiva) al verificar cada elemento, siempre que el anlisis estructural de segundo orden sea
realizado introduciendo la imperfeccin del prtico cuando se deforma lateralmente.

En ambos cdigos de diseo se presentan alternativas para llevar a cabo los anlisis, ya sean de
primer o segundo orden, del tipo de anlisis escogido depende el posterior diseo de la estructura,
por ejemplo si el anlisis es de primer orden elstico es necesario verificar la estabilidad de cada
miembro, sin embargo si se realiza un anlisis de segundo orden, incluyendo las imperfeccines
iniciales que se indican, solo se requiere chequear la capacidad del miembro mas solicitado, esto
indica que un anlisis de segundo orden de la estructura permite simplificar los clculo a la hora
de disear los elementos.
72

Con respecto a la situacin en Chile se concluye que no se ha profundizado en el tema de anlisis
de segundo orden de las estructuras, debido a que este tipo de clculos se hace necesario cuando
la estructura es muy alta o est sometida a grandes cargas axiales como por ejemplo los edificios
en altura en donde la influencia de las acciones horizontales es importante.

En Chile el acero se comenz a utilizar como material de construccin a mediados del siglo XIX.
En un principio, se levantaron solo puentes y actualmente alrededor del 1% de los edificios mas
altos que existen en nuestro pas se construyeron en base a acero y el resto se levanta con acero u
hormign o solo con este ltimo elemento. Adems, el presidente del Instituto Chileno del Acero,
Elas Arze, dijo en una entrevista que a pesar de que el acero no se utiliza masivamente en
edificaciones de altura en el pas, s es un elemento importante en la construccin de edificios
habitacionales de hasta cinco pisos y en bodegas industriales, y en cuanto a las normas de las
cuales se dispone para edificaciones en altura, segn miembros de la directiva del Colegio de
Ingenieros e Instituto Chileno del Acero, son ms estrictas que las aplicadas actualmente en
Estados Unidos. (Portal del medio ambiente, 2001).

En cuanto a las funciones de estabilidad, (s
ii
y s
ij
) y su implementacin en los programas de
estructuras existentes, la dificultad radica en que estas son exponenciales y trigonomtricas, como
se muestra en en anexo A, por lo que se requieren distintos tipos de funciones de acuerdo con el
signo del esfuerzo axial, adems son indeterminadas cuando P = 0. Con el fin de solucionar estos
problema se presentaron en esta tesis las aproximaciones a estas funciones en forma polinmica
(seccin 2.2).

Finalmente, se puede ver que la idea de llevar a cabo un anlisis elstico de segundo orden de la
estructura es principalmente llegar a obtener resultados que reflejen las condiciones reales del
problema, y de esta manera simplificar luego el diseo, chequeando solamente que el miembro
que este siendo mas solicitado resista la carga que resulte luego del anlisis.

6.2. Aportaciones mas significativas

En este trabajo se han presentado los aspectos que determinan la estabilidad de una
estructura, y tres maneras de abordar el tema de un anlisis elstico de segundo orden, primero
desde el punto de vista terico, exponiendo los mtodos utilizados luego de la revisin de libros y
revistas, luego desde el punto de vista computacional, utilizando dos programas de clculo Proges
y SAP2000, y en tercer lugar a travs de normas de diseo: AISC y EC3. sto permite enriquecer
los conceptos que se tienen en cuanto a la influencia de las deformaciones gemetricas cuando se
lleva a cabo el anlisis de una estructura, y da pie a un estudio mas completo acerca del tema.
73

La tesis muestra los pasos a seguir para realizar un anlisis de segundo orden a una estructura tipo
marco, usando el programa SAP2000, y verifica los resultados obtenidos a travs de la
comparacin con otros mtodos de clculos.

La tesis presenta un estudio comparativo de las normas AISC y EC3, con respecto a la estabilidad
de marcos traslacionales o no arriostrados. De los resultados comparativos obtenidos, se han
detectado diferencias significativas, sobre todo a la hora de determinar cuando se hace necesario
un anlisis de segundo orden de la estructura, mediante los factores de amplificacin B
2
(para
AISC) y F
A
(EC3), cuyos valores son distintos, como se pudo ver en la tabla 5.3.2.

Se muestra la forma en que se deducen las funciones de estabilidad, asi mismo se muestra a
travs de un ejemplo bien detallado como llevar a cabo un anlisis elstico de segundo orden
utilizando las aproximaciones polinmicas de las Funciones de Estabilidad, con esto se logra
establecer las bases para que posteriormente se lleven a cabo estudios y verificaciones mas
completos que permitan validar sus resultados y as poder utilizarlos en nuestro pas.

La tesis presenta los estudios recientes que se han efectuado respecto a anlisis elstico de
segundo orden y los problemas de inestabilidad que pueden haber en una estructura, asi como
tambin la manera de considerarlos en los anlisis. Es importante comenzar a tomar en cuenta
estos factores en Chile ya que hoy en da existe la tendencia a disear y construir este tipo de
edificaciones. Un ejemplo claro es la Torre Titaniun (ver nota 1) que es uno de los primeros de
varios edificios de altura que pueden llegar a construrirse en Chile y es importante conocer el
comportamiento de estas estructuras tal y como se estudian en otros pases.

Nota1: Con 192 metros de altura, Titanium La Portada ser uno de los edificios ms altos de la regin y
el primero en altura que se construye en Chile, en el centro de negocios denominado El Golf, donde se
ubican las oficinas de las principales compaas nacionales y extranjeras presentes en el pas.
El edificio que se entregar a mediados del 2009 ya tiene construidos los siete niveles de subterrneos
(1350) y 15 de los 52 pisos que ser destinado exclusivamente a oficinas de plantas libres de entre 1.350 y
1.950 metros cuadrados.
Diseado por el arquitecto chileno, Abraham Senerman Lamas, Titanium La Portada pertenece a
Inmobiliaria Titanium, creada para el desarrollo de proyectos emblemticos y de renovacin urbana.
Con un diseo de vanguardia, el proyecto incorpor tecnologa de ltima generacin que apunta a
convertirlo en un edificio sustentable que ya cuenta con la pre certificacin del Consejo Norteamericano
de edificios verdes, LEED CS, para edificios de plantas libres y que una vez concluido, espera obtener la
certificacin GOLD.
Titanium La Portada cuenta con asesores de gran prestigio, entre los que destaca Joseph Colaco, revisor
de las materias de clculo estructural del edificio, quien ha trabajado en diferentes proyectos como la
torre Burj Dubai. (www.Urbanity.es, 2008)
74

6.3. Futuras lneas de investigacin

Se puede profundizar en el tema del anlisis no lineal considerando la degradacin del
material, es decir realizar un anlisis de segundo orden inelstico adems se puede estudiar que
sucede con el pandeo local de los elementos.

Realizar el anlisis de segundo orden elstico a un edificio en altura utilizando la matriz
geomtrica del elemento y luego modelarla en SAP2000 para comparar los resultados, tambin se
pueden incluir otro tipo de estructuras, para ver la respuesta de estas frente a los anlisis en donde
se incluyen los efectos de segundo orden.

Otra posibilidad sera estudiar diferentes tipos de aproximaciones a las funciones de estabilidad, y
determinar cual de todas ellas tiene mejor convergencia y se aproxima mejor a los resultados
exactos, para poder emplearlas en un programa computacional.

Captulo VII Bibliografa
75

Captulo VII
Bibliografa

AISC (American Institute of Steel Construction). March 2005. Specification for Structural Steel
Buildings. One East Wacker Drive, Suite 700; Chicago, Illinois 60601-1802. 519 p.

Charles G. Salmon, J ohn E. J ohnson. 1996. Steel structures design and behabior. 4 ed.
HarpercollinsCollegePublishers. 1024 p.

Computers and structures, Inc. Abril 2007. CSI Analysis reference Manual, for SAP2000,
ETABS and SAFE. University Avenue Berkeley, California 94704 USA. 450 p.

Computers and Structures, Inc. May 2000. Integrated Finite Element Analysis and Design of
Structures steel design manual. Berkeley, California, USA. 169 p.

Computers and Structures, Inc. Oct 1998. Integrated Finite Element Analysis and Design of
Structures analysis reference. Berkeley, California, USA. 438 p.

Eurocdigo 3. Dec 2003. Design of steel structures - Part 1-1: General rules and rules for
buildings. European Commitee For Standarization (CEN). 90 p.

Eurocdigo. J un 2003. Bases de clculo de estructuras. European Commitee For Standarization
(CEN). 84 p.

Hernndez, M., Enrique; Ra, A., Edelmiro. Oct1994. Clculo en segundo orden en la normativa
europea. Revista de obras publicas. N 3336.

Instituto chileno del acero (ICHA). 2000. Manual de diseo para estructuras de acero, Mtodo de
factores de carga y resistencia. Tomo I. Arze, Racine y asociados. 812 p.

Kim, Seung-Eock. 1999. An Innnovative Design For Steel Frame Using Advanced
Analysisfootnotemark. Structural Engineering Handbook. Ed. Chen Wai-Fah, Boca Raton: CRC
Press LLC. 58 p.

Kim, Seung-Eock; Lee, J ong-She; Choi, Se-Hyu; Kim, Chang-Sung. 2004. Practical second-order
inelastic analysis for steel frames subjected to distributed load. J ournal Engineering Structures.
Vol 26: 51-61. 11 p.
76

Liew, J .Y.; Shanmugam, N.W. and Yu, C.H. 1999. Structural Analysis. Structural Engineering
Handbook. Ed. Chen Wai-Fah, Boca Raton: CRC Press LLC. 189 p.
Maleck, S., Andrea, E.; White, Donald, W. Aug2004. Alternative approaches for elastic analysis
and design of steel frames I: Overview. J ournal of structural engineering. Vol 130, Issue 8. 12 p.

Lopetegui Torres, J ulio. 2006. Apuntes de clases, Ramo: Computacin para ingeniera
estructural.

Ortega, Miguel A.; Romero ,J os L.; de la Rosa, Emilio. Oct.-Dic. 2007. Un estudio histrico del
problema de las piezas prismticas rectas sometidas a compresin. Parte II. Informes de la
construccin. Vol. 59, 508, 61-71. 11 p.

Perles, Pedro; Santiago, Liliana; Tommasini, Flavia. 2003. Temas de estructuras especiales.
Publicado por Nobuko. Inc NetLibrary. 262p.

Russell C. Hibbeler. 1997. Anlisis structural. 3 ed. Prentice-Hall Hispanoamrica, S.A. 727 p.

Saavedra Chamorro, P. I. 2005. Anlisis de Estabilidad Elstica Flexo-Torsional para Barras de
Acero, de acuerdo a Normas y Programas Computacionales. Tesis de pregrado. Universidad
Austral de Chile. 97p.

Saouna, Victor E. 1999. Stuctural analysis. 1 ed. Dept. of Civil Environmental and Architectural
Engineering University of Colorado, Boulder, CO 80309-0428. 215 p.

Tong, G.; Xing, G. J unio 2007. A comparative study of alternative Approaches for stability
Design of steel Frames. J ournal Advances in structural Engeneering. Vol. 10, 4: 455-466. 12p.
Saouna, Victor E. 1999. Matrix Stuctural analysis, with an Introduction to Finite Elements. 1 ed.
Dept. of Civil Environmental and Architectural Engineering University of Colorado, Boulder, CO
80309-0428. 246 p.

Wunderlich, Walter; Pilkey D., Walter . 2003. Mechanics of Structures: Variational and
Computational Methods. Structural Engineering Handbook. 2 Ed. Boca Raton London New
York Washington, D.C. 877 p.

77

Yong Ayn, D. J . Abril 2007. Diseo de elementos estructurales de acero sometidos a flexo-
compresin: Desarrollo terico y estudio comparativo. Tesis Doctorado. Escuela superior de
ingenieros San sebastian, Universidad de Navarra. 242 p.

Portal del medio ambiente. Sep 2001. Reportajes y entrevistas:Normas chilenas seran mas
estricttas que en USA. www.sustentable.cl

www.Urbanity.es. Mayo, 2008. Seccin rascacielos y highrises: Torre titanium la Portada,
Santiago de Chile

Anexo A Funciones de Estabilidad
78

Anexo A
Funciones de estabilidad

Ecuaciones de equilibrio de segundo orden

Las ecuaciones de equilibrio de segundo orden se deducen tomando como base la
configuracin deformada del eje de la barra, teniendo en cuenta la interaccin de las diferentes
fuerzas exteriores al modificarse los efectos de algunas a causa de las deformaciones producidas
por otras; como las deformaciones y las fuerzas interiores dejan de ser independientes, al plantear
las ecuaciones de equilibrio deben considerarse todos los elementos mecnicos simultneamente.
El objeto fundamental que se persigue al obtener las ecuaciones de equilibrio de segundo
orden es utilizarlas en el estudio de problemas de inestabilidad que no pueden resolverse por
medio de las de primer orden, basadas en el equilibrio de las fuerzas que actuan sobre la barra no
deformada.
En a figura A.1 se han dibujado por separado las fuerzas alojadas en cada uno de los
planos yz y xz, con el fin de poder deducir las ecuaciones de equilibrio para la barra deformada; P
aparece en ambas puesto que su linea de accin coincide con la interseccin de los dos planos.

Figura A.1. Fuerzas alojadas en cada uno de los planos principales de la barra
Fuente: Lpez, 1999

79

EI
x
:
II
+P: - jH
A
-
z
I
(H
A
+H
B
)[ = H
Ax
-
z
I
(H
Ax
+H
Bx
) cc A. 1

EI
u
II
+Pu - jH
Ax
-
z
I
(H
Ax
+H
Bx
)[ = H
A
-
z
I
(H
A
+H
B
) cc A. 2

EC
u
III
-(0K
1
+K
)
I
-u
I
jH
Ax
-
z
I
(H
Ax
+H
Bx
)[ -:
I
jH
A
-
z
I
(H
A
+H
B
)[ -
-
:
I
(H
A
+H
B
) -
u
I
(H
Ax
+H
Bx
) = u cc A. S
K
= _
A
o
2
A
JA

EI
x
y EI
son las rigideces a la flexin

0K
1
y EC
u
son las rigideces por torsin y por resistencia al alabeo de St. Venant

Estas tres ecuaciones no son independientes entre s, pues en las dos primeras, que
corresponden a flexin alrededor de los ejes principales de las secciones, aparece tambin el el
ngulo de torsin , y en la tercera, que representa la condicin de equilibrio alrededor del eje
longitudinal de la barra, los desplazamientos u y v. (Lopez, 1992).
Estas ecuaciones son la base para el clculo fuerzas que puedan producir pandeo en una
barra, y se reducen de acuerdo a las condiciones de borde y restricciones que presente cada
problema en particular.

Anlisis de un elemento viga columna

Figura A.2. Viga columna sujeta a momentos en sus extremos
Fuente: Liew, 1999

La ecuacin de la elstica de la viga, deducida teniendo en cuenta los momentos
flexionantes adicionales producidos por la fuerza P al obrar sobre el eje deformado por la accin
del par aplicado en el extremo, permite determinar el ngulo de giro
A
en funcin del momento
M
A
, la fuerza P y las otras propiedades de la viga (figura A.1).
80

Si si se acepta que las deformaciones son pequeas de manera que puede utilizarse la
expresin simplificada de la curvatura, la ecuacin diferencial de la elstica es:

EI
J
4
y
Jx
4
+P
J
2
y
Jx
2
= w cc. A. 4

La solucin de esta ecuacin, teniendo en cuenta las condiciones de borde, proporciona
una expresin que permite calcular la flecha v en un punto cualquiera del eje de la barra, en
funcin de sus propiedades geomtricas y mecnicas, de la fuerza P y de los momentos en los dos
extremos; utilizandola es fcil obtener en ngulo de giro en A y, partiendo de l, la rigidez
angular y el factor de transporte angular correspondiente.
Las expresiones de la rigidez angular y el factor de tranporte de las barras flexo-
comprimidas son bastante complejas pero, pueden escribirse en funcin de ciertos coeficientes
que dependen de la intensidad de la fuerza axial y que, conocidos como Funciones de Estabilidad,
han sido tabulados por varios autores. (Lopez, 1992).

Figura A.3. Viga columna sujeta a momentos en sus extremos
Fuente: Liew, 1999
EI
J
4
y
Jx
4
+P
J
2
y
Jx
2
= w cc. A. 4
EI
J
4
y
Jx
4
+P
J
2
y
Jx
2
= u cc. A. S

La solucin de la ecuacin diferencial A.5 es la siguiente:

y = Ascn(kx) +Bcos(kx) +Cx + cc. A. 6
Con k = PEI

Aplicando las siguientes condiciones de borde se obtienen los valores de las constantes de
la ec. A.6:
x = 0 y = 0 y = M
A
/ EI
x = L y = 0 y = -M
B
/ EI
81

A =
H
A
cos(kl) +H
B
EIk
2
scn(kl)
B = -
H
A
EIk
2

C = -_
H
A
+H
B
EIk
2
I
] =
H
A
EIk
2

Asi se obtiene la ecuacin A.7 que es la deformacin de la viga en funcin de x que va a lo largo
de la longitud.

y =
1
EIk
_
cos(kI)
scn(kI)
-cos(kx) -
x
I
+1_ H
A
+
1
EIk
2
_
1
scn(kI)
scn(kx) -
x
I
_ H
B
cc. A. 7

Derivando la ecuacin A.7 se obtiene el angulo de giro en funcin de x

y
i
=
1
EIk
_
cos(kI)
scn(kI)
cos(kx) +sen(kx) -
1
kI
_ H
A
+
1
EIk
_
cos(kx)
scn(kI)
-
1
kI
_ H
B
cc. A. 8

Se evalua esta funcin en x = 0 y x = L para obtener los ngulos de giro en los extremos

0
A
= y
i
(x = u) =
1
EIk
_
cos (kI)
scn(kI)
-
1
kI
_ H
A
+
1
EIk
_
1
scn(kI)
-
1
kI
_ H
B

=
I
EI
_
kIcos(kI) -scn(kI)
(kI)
2
scn(kI)
_ H
A
+
I
EI
_
kI -scn(kI)
(kI)
2
scn(kI)
_ H
B
cc. A. 9

0
B
= y
i
(x = I) =
1
EIk
_
1
scn(kI)
-
1
kI
_ H
A
+
1
EIk
_
cos (kI)
scn(kI)
-
1
kI
_ H
B

=
I
EI
_
kI -scn(kI)
(kI)
2
scn(kI)
_ H
A
+
I
EI
_
kIcos(kI) -scn(kI)
(kI)
2
scn(kI)
_ H
B
cc. A. 1u

Se definen los valores y , y luego se ordenan los resultados en una matriz

=
1
kI
_coton(kI) -
1
kI
] y =
1
kI
_
kI
kIscn(kI)
-
1
kI
]

_
B
_ = _

_ _
H
A
H
B
_

La inversa de esta matriz permite obtener los valores de los momentos en funcin de los
los angulos
A
y
B
, para definir de esta manera los valores de s
ii
y s
ij

82

H
A
=
EI
I
(s
0
A
+s
]
0
B
) cc. A. 11
H
B
=
EI
I
(s
]
0
A
+s
]]
0
B
) cc. A. 12
s
= s
]]
=
kIscn(kI) -(kI)
2
cos (kI)
2 -2cos(kI) -kIscn(kI)
cc. A. 1S
s
]
= s
]
=
(kI)
2
-kIscn (kI)
2 -2cos(kI) -kIscn(kI)
cc. A. 14

Estas funciones consideran el efecto p- dentro del elemento

Figura A.4. Grfico de Funciones de Estabilidad
Fuente: Liew et al, 1999

Ahora para considerar el efecto P-, se da un desplazamiento a la viga-columna de la
figura A.3 en el extremo B como se muestra en la figura A.5, con esto se deducen las siguientes
ecuaciones:

H
A
=
EI
I
_s
_0
A
-
A
I
] +s
]
_0
B
-
A
I
]_
=
EI
I
_s
0
A
+s
]
0
B
-(s
+s
]
)
A
I
_ cc. A. 1S
H
B
=
EI
I
_s
]
_0
A
-
A
I
] +s
_0
B
-
A
I
]_
=
EI
I
_s
]
0
A
+s
0
B
-(s
+s
]
)
A
I
_ cc. A. 16

83

Del equilibrio de momentos en torno al punto A se tiene que:

XH
A
= u = H
A
+H
B
+SI +PA = u cc. A. 17
S = -
1
I
(PA +H
A
+H
B
) cc. A. 18
Reemplazando M
A
y M
B
en ec A.18

S = -
1
I
_PA +
EI
I
_s
0
A
+s
]
0
B
-s
A
I
-s
]
A
I
] +
+
EI
I
_s
]
0
A
+s
0
B
-s
A
I
-s
]
A
I
]_ cc. A. 19

Con k = PEI y ordenando los elemento se llega a la ec. A.20

S = -
EI
I
_A_
k
2
l
2
-2s
-2s
]
I
2
_ +0
A
_
s
+s
]
I
] +0
B
_
s
+s
]
I
]_ cc. A. 2u

Figura A.5. Viga columna sujeta a momentos en los extremos y desplazamientos
Fuente: Liew, 1999

Escribiendo en forma matricial los resultados de M
A
, M
B
y S para los nodos A y B de la
viga-columna se llega a la siguiente matriz de rigidez que incorpora la deformacin P- y P-
producidas en los elementos.

l
l
l
l
l
l
l
l
l
l
l
l
l
l
l
EA
l
u u -
EA
l
u u
u
|2(s
+s
]
) -e
2
]EI
l
3

-(s
+s
]
)EI
l
2
u
-|2(s
+s
]
) -e
2
]EI
l
3
-(s
+s
]
)EI
l
2

u
-(s
+s
]
)EI
l
2

s
EI
l
u
(s
+s
]
)EI
l
2

s
]
EI
l

-
EA
l
u u
EA
l
u u
u -
|2(s
+s
]
) -e
2
]EI
l
3

(s
+s
]
)EI
l
2
u
|2(s
+s
]
) -e
2
]EI
l
3

(s
+s
]
)EI
l
2

u
-(s
+s
]
)EI
l
2

s
]
EI
l
u
(s
+s
]
)EI
l
2

s
EI
l

1
1
1
1
1
1
1
1
1
1
1
1
1
1
1

Anexo B Anlisis de la estructura por el Mtodo de la Matriz Geomtrica
84

Anexo B
Anlisis de la estructura por el Mtodo de la Matriz Geomtrica

En este anexo se muestra en detalle el anlisis realizado a la estructura por medio del
Mtodo de la Matriz Geomtrica.

Figura B.1. Modelo de marco y condiciones de carga

Figura B.2. Sistema de coordenadas glovales y locales considerados en mtodo 1

85

Se comienza realizando un anlisis elstico de primer orden a la estructura para
determinar los valores de N1, N2 y N3 que son las fuerzas axiales sobre las barras uno, dos y tres,
para determinar las correspondientes matrices geomtricas.

N1 =-1116,37kN
N2 =0
N3 =- 2093,63kN

Luego se definen las matrices de rigidez elstica y geomtrica para cada elemento como
se explica en el captulo 2

Barra 1 con N1 = -1116.37kN

Figura B.3. Barra 1 del marco

k
c
=
EI
l
3
l
l
l
l
l
l
Al
2
I
u
u
-Al
2
I
u
u
u
12
-6l
u
-12
-6l
u
-6l
4l
2
u
6l
2l
2
-Al
2
I
u
u
Al
2
I
u
u
u
-12
6l
u
12
6l
u
-6l
2l
2
u
6l
4l
2
1
1
1
1
1
1
`
1
1
1
1
u
u
w
u
u
u
b
w
b
b
1
1
1
1
1

86

k
c
=
l
l
l
l
l
l
u
u
u
u
u
u

u
21u9.S7S
-84S7.S
u
-21u9.S7S
-84S7.S

u
-84S7.S
4Suuu
u
84S7.S
22Suu

u
u
u
u
u
u

u
-21u9.S7S
84S7.S
u
21u9.S7S
84S7.S

u
-84S7.S
22Suu
u
84S7.S
4Suuu
1
1
1
1
1
1

k
g
=
N1
l
l
l
l
l
l
l
l
l
l
l
l
u u u u u u
u
6
S
-
l
1u

u -
6
S
-
l
1u

u -
l
1u

2l
2
1S

u
l
1u
-
l
2
Su

u u u u u u
u -
6
S

l
1u

u
6
S

l
1u

u -
I
1u
-
l
2
Su

u
l
1u

2l
2
1S

1
1
1
1
1
1
1
1
1
1
1
`
1
1
1
1
u
u
w
u
u
u
b
w
b
b
1
1
1
1
1

k
g
=
l
l
l
l
l
u
u
u
u
u
u

u
-167.46
111.64
u
167.46
111.64

u
111.64
-119u.8
u
-111.64
297.7

u
u
u
u
u
u

u
167.46
-111.64
u
-167.46
-111.64

u
111.64
297.7
u
-111.64
-119u.8
1
1
1
1
1

Por lo tanto k
T
=k
e
+k
g
queda expresado de la siguiente manera

k
t
=
l
l
l
l
l
l
u
u
u
u
u
u

u
1941.92
-8S2S.86
u
-1941.92
-8S2S.86

u
-8S2S.86
4S8u9.2u
u
8S2S.86
22797.7

u
u
u
u
u
u

u
-1941.92
8S2S.86
u
1941.92
8S2S.86

u
-8S2S.86
22797.7
u
8S2S.86
4S8u9.2u
1
1
1
1
1
1

Rotando en 90 se tiene la matriz de rigidez en las coordenadas glovales

u
1
w
1
0
1
u
2
w
2
0
2

k
t
=
l
l
l
l
l
l
1941.92
u
-8S2S.86
u
-1941.92
-8S2S.86

u
u
u
u
u
u

-8S2S.86
u
4S8u9.2u
8S2S.86
u
22797.7

u
u
u
u
u
u

-1941.92
u
8S2S.86
1941.92
u
8S2S.86

-8S2S.86
u
22797.7
8S2S.86
u
4S8u9.2u
1
1
1
1
1
1
F
1
I
1
H
1
F
2
I
2
H
2

Las fuerzas que se generan en la elemento empotrado en sus extremos por la carga
distribuida son las siguientes:
87

-p
10
= w
l
l
l
l
-l2
l
2
12
-l2
-l
2
121
1
1
1
= _
S6
-48
S6
48
_
I
1
10
H
1
10
I
2
10
H
2
10

F
X1
10
H
1
10
F
X2
10
H
2
10

Barra 2 con N2 = 0
(solo se tiene la matriz de rigidez elstica para este elemento)


k
c
=
EI
l
3
l
l
l
l
l
l
Al
2
I
u
u
-Al
2
I
u
u
u
S
-Sl
u
-S
u
u
-Sl
Sl
2
u
Sl
u
-Al
2
I
u
u
Al
2
I
u
u
u
-S
Sl
u
S
u
u
u
u
u
u
u1
1
1
1
1
1
`
1
1
1
1
u
u
w
u
u
u
b
w
b
b
1
1
1
1
1

u
2
w
2
0
2
u
3
w
3
0
3

k
c
= k
t
=
l
l
l
l
l
u
u
u
u
u
u

u
S22
-S22u
u
-S22
u

u
-S22u
S22uu
u
S22u
u

u
u
u
u
u
u

u
-S22
S22u
u
S22
u

u
u
u
u
u
u
1
1
1
1
1
F
2
I
2
H
2
F
S
I
S
H
S

La viga esta horizontal por lo que sus coordenadas locales coinciden con las glovales por
lo tanto no es necesario rotarla
Las fuerzas que se generan en la elemento empotrado en sus extremos por la carga
distribuida son las siguientes:

88

-p
20
= q _
-Sl8
l
2
8
-Sl8
_ = _
-S7.S
7S
-22.S
_
I
2
20
H
2
20
I
3
20

F
z2
20
H
2
20
F
z3
20

Barra 3 con N3 = -2093,63kN


k
c
=
EI
l
3
l
l
l
l
l
l
Al
2
I
u
u
-Al
2
I
u
u
u
u
u
u
u
u
u
u
u
u
u
u
-Al
2
I
u
u
Al
2
I
u
u
u
u
u
u
u
u
u
u
u
u
u
u1
1
1
1
1
1
`
1
1
1
1
u
u
w
u
u
u
b
w
b
b
1
1
1
1
1

k
g
=
N
0
l
l
l
l
l
l
l
l
u u u u u u
u 14 u u -14 u
u u u u u u
u u u u u u
u -14 u u 14 u
u u u u u u

1
1
1
1
1
1
1
`
1
1
1
1
u
u
w
u
u
u
b
w
b
b
1
1
1
1
1

Como no se est considerando la deformacin axial de los elementos la matriz total en
este caso es igual a la matriz geomtrica

k
g
=
l
l
l
l
l
u
u
u
u
u
u

u
-S2S.41
u
u
S2S.41
u

u
u
u
u
u
u

u
u
u
u
u
u

u
S2S.41
u
u
-S2S.41
u

u
u
u
u
u
u
1
1
1
1
1

89

Al rotar en 90 se tiene:
u
3
w
3
0
3
u
4
w
4
0
4

k
g
= k
t
=
l
l
l
l
l
-S2S.41
u
u
S2S.41
u
u

u
u
u
u
u
u

u
u
u
u
u
u

S2S.41
u
u
-S2S.41
u
u

u
u
u
u
u
u

u
u
u
u
u
u
1
1
1
1
1
F
S
I
S
H
S
F
4
I
4
H
4

Se superponen las matrices de cada uno de los elementos para conformar la matriz de
rigidez gloval total K
totuI
que es la siguiente:

l
l
l
l
l
l
l
l
l
l
l
1941.92
u
-8S2S.86
-1941.92
u
-8S2S.86
u
u
u
u
u
u

u
u
u
u
u
u
u
u
u
u
u
u

-8S2S.86
u
4S8u9.2u
8S2S.86
u
22797.7
u
u
u
u
u
u

-1941.92
u
8S2S.86
1941. 92
u
8325. 8
u
u
u
u
u
u

u
u
u
u
S22
-S22u
u
-S22
u
u
u
u

-8S2S.86
u
22797.7
8325. 8
-S22u
99. 2
u
S22u
u
u
u
u

u
u
u
u
u
u
-523. 41
u
u
S2S.41
u
u

u
u
u
u
-S22
S22u
u
S22
u
u
u
u

u
u
u
u
u
u
u
u
u
u
u
u

u
u
u
u
u
u
S2S.41
u
u
-S2S.41
u
u

u
u
u
u
u
u
u
u
u
u
u
u

u
u
u
u
u
u
u
u
u
u
u
u

1
1
1
1
1
1
1
1
1
1
1

l
l
l
l
l
l
l
l
l
l
l
l
l
u
1
w
1
0
1
u
2
w
2
0
2
u
3
w
3
0
3
u
4
w
4
0
4
1
1
1
1
1
1
1
1
1
1
1
1
1

w
1
= u
1
= 0
1
= w
4
= u
4
= w
2
= w
3
= u

Y como u
2
= u
3
el valor -523.41 se suma a 1941.92 por lo tanto la matriz de rigidez total
se reduce a una matriz de 4x4 que se da a continuacin:

K
totuI
= j
1418.S1 8S2S.86
8S2S.86 96uu9.2
[
K
totuI
: -P
0
= P
-

P
totuI
= P
-
+P
0

P
-
= corgos puntuolcs
: = Jcormocioncs incognitos
P
totuI
= j
-82.S
u.u

-S6
-48

u
-7S
[ = j
-118.S
-12S
[

Se resuelve entonces el siguiente sistema de ecuaciones:

j
1418.S1 8S2S.86
8S2S.86 96uu9.2
[ j
u
2
0
2
[ = j
-118.S
-12S
[
j
u
2
0
2
[ = j
-u.1S482228
u.u1214497
[

90

Ahora se determinan las fuerza resultantes para cada miembro
P
= k
-p
0

Para el elemento 1

p
1
=
l
l
l
l
l
l
u
u
u
u
u
u

u
1941.92
-8S2S.86
u
-1941.92
-8S2S.86

u
-8S2S.86
4S8u9.2u
u
8S2S.86
22797.7

u
u
u
u
u
u

u
-1941.92
8S2S.86
u
1941.92
8S2S.86

u
-8S2S.86
22797.7
u
8S2S.86
4S8u9.2u
1
1
1
1
1
1
l
l
l
l
l
l
u
u
u
u
-u.1S48
u.uu121
1
1
1
1
1
1
+
l
l
l
l
l
l
S6
u
-48
S6
u
48
1
1
1
1
1
1

p
1
=
l
l
l
l
l
l
2SS.SSS
u
-1u6u.1S21
-16S.SSS
u
7u8.986
1
1
1
1
1
1

Para el elemento 2

p
2
=
l
l
l
l
l
l
u
u
u
u
u
u

u
S22
-S22u
u
-S22
u

u
-S22u
S22uu
u
S22u
u

u
u
u
u
u
u

u
-S22
S22u
u
S22
u

u
u
u
u
u
u
1
1
1
1
1
1
l
l
l
l
l
l
u
u
u.u121
u
u
u
1
1
1
1
1
1
+
l
l
l
l
l
l
u
-S7.S
7S
u
-22.S
u
1
1
1
1
1
1

p
2
=
l
l
l
l
l
l
u
-1uu.897
7u8.968
u
4u.897
u
1
1
1
1
1
1

Para el elemnto 3

p
3
=
l
l
l
l
l
l
u
u
u
u
u
u

u
-S2S.41
u
u
S2S.41
u

u
u
u
u
u
u

u
u
u
u
u
u

u
S2S.41
u
u
-S2S.41
u

u
u
u
u
u
u
1
1
1
1
1
1
l
l
l
l
l
l
u
-u.1S482
u
u
u
u
1
1
1
1
1
1
+
l
l
l
l
l
l
u
u
u
u
u
u
1
1
1
1
1
1

p
3
=
l
l
l
l
l
l
u
81.uSS
u
u
-81.uSS
u
1
1
1
1
1
1

91

Figura B.6. Determinacin de las fuerzas longitudinales

El equilibrio de estas fuerzas permitir determinar las nuevas fuerzsas axiales que se
actuarn en cada elemento para comenzar con la siguinete iteracin, todo esto se realiza hasta que
los valores de las defomaciones son iguales entre el ltimo paso y el anterior.

XF
X
= u: N
2
2
-P
X2
-I
2
1
= u N
2
2
= 82.S -16S.SSS = -81.uSS
XF
X
= u: N
2
3
+I
3
3
= u N
2
3
= -81.uSS

XF
z
= u: N
2
1
+P
z2
+I
2
2
= u N
2
1
= -1uSu -1uu.897 = -11S1.897
XF
z
= u: N
3
3
+P
z3
-I
3
2
= u N
3
3
= 4u.887 -21uu = -2uS9.11S

Se determinar a continuacin el factor de carga crtica de pandeo para el marco propuesto
utilizando el Mtodo de la Matriz Geomtrica
92

Figura B.7. geometra y combinaciones de carga

Con las condiciones de carga se determinar el factor de carga () que multiplica a todas
las cargas y para el cual la matriz total del sistema es cero, esto se realizar de la siguiente forma:
Tomando la matriz lineal del elemento la cual es la siguiente:

K
c
= j
21u9.S8 8S2S.Su
84S7.Su 972uu.uu
[

Ahora se toma la matriz geomtrica total y se multiplica por el factor de carga
K
g
= j
-687.4S -11S.uS
-11S.uS -1S86.uS
[

K
totuI
: -P
0
= P
-

(k
c
+k
g
): = u

Jct _j
-687.4S -11S.uS
-11S.uS -1S86.uS
[ + j
-687.4S -11S.uS
-11S.uS -1S86.uS
[_ = u
9S9S82
2
-678u1922 +1SS8S9844 = u
2
-72.162 +142.446 = u
1
= 2.uS1 y
2
= 7u.1S
93

Figura B.8. Valores del determinante de K
total
usando la matriz geomtrica

Figura B.9. Formas de modos de pandeo correspondientes a
1
y
2


En la tabla B.1 se muestran los factores de carga aplicados a las fuerzas axiales para los
distintos desplazamientos obtenidos a travs del mtodo de la Matriz Geomtrica que se ven en la
tabla 5.3.2.

94

Matriz Geomtrica Factor de Carga
(cm)
1,9 1010,32
2,33 637,18
2,85 368,19
3,43 258,16
4,08 198,25
4,81 160,51
5,63 134,51
6,56 115,47
7,63 100,87
8,87 89,29
10,33 79,83
12,06 71,9
14,16 65,11
16,47 59,17
19,99 53,85
24,2 48,97
29,86 44,37
37,81 39,9
49,64 35,39
Tabla B.1. Valores de deformacin vs factor carga

Grfico B.1. Valores de deformacin vs factor carga

0
100
200
300
400
500
600
700
800
900
1000
1100
0 5 10 15 20 25 30 35 40 45 50 55
Factor de carga
(cm)
Anexo C Archivos de entrada y salida de datos
95

Anexo C
Archivos de entrada y salida de datos

Este anexo muestra 9 archivos que corresponden a los archivos de entrada de datos, en
donde se define la estructura con sus coordenadas, geometra propiedades del material, etc,
generados por Datos.exe, y salida de datos, en donde se muestran los resultados de
desplazamientos y fuerzas internas de los elementos asi como tambin el nmero de iteraciones
realizadas para los anlisis de segundo orden, generados por el programa PROGES y SAP2000.
Se define N2 como la carga aplicada en el nudo dos y N3 la carga aplicada en el nudo tres
del marco de la figura 5.1

ARCHIVO C.1
Archivo de entrada de datos PROGES

* Ejemplo libro Mechanics of Structures. Captulo 11 pgina 648.
13 12 45000.200000000. 7.
2 .00100000 1 1.00000 1.40000 1.30000
**COORDENADAS PRINCIPIO
1 .0000 .0000 .0000 111111
2 .0000 .0000 8.0000 101010
3 10.0000 .0000 8.0000 101010
4 10.0000 .0000 4.0000 101111
5 .0000 .0000 2.0000 101010
6 .0000 .0000 4.0000 101010
7 .0000 .0000 6.0000 101010
8 2.5000 .0000 8.0000 101010
9 5.0000 .0000 8.0000 101010
10 7.5000 .0000 8.0000 101010
11 10.0000 .0000 7.0000 101010
12 10.0000 .0000 6.0000 101010
13 10.0000 .0000 5.0000 101010
**FIN COORDENADAS
**TIPOS PRINCIPIO
1 88. 86840. 48090. 10. 210. 80. 5049. 3757.
2 88. 45000. 76070. 7. 210. 80. 3548. 4610.
**BARRAS PRINCIPIO
1 1 5 11 2 0. 0.
2 5 6 11 2 0. 0.
3 6 7 11 2 0. 0.
4 7 2 11 2 0. 0.
5 2 8 11 1 0. 0.
6 8 9 11 1 0. 0.
7 9 10 11 1 0. 0.
8 10 3 11 1 0. 0.
96

9 3 11 01 2 0. 0.
10 11 12 11 2 0. 0.
11 12 13 11 2 0. 0.
12 13 4 11 2 0. 0.
**FIN BARRAS
**CASOS
1a
**CARGA-BARRAS
a 1 2 1 9.000 .000 .000 .000
a 2 2 1 9.000 .000 .000 .000
a 3 2 1 9.000 .000 .000 .000
a 4 2 1 9.000 .000 .000 .000
a 5 2 1 -6.000 .000 .000 .000
a 6 2 1 -6.000 .000 .000 .000
a 7 2 1 -6.000 .000 .000 .000
a 8 2 1 -6.000 .000 .000 .000
**FIN CARGA-BARRAS
**CARGA-NUDOS
a 2 .0000 .0000 .0000 -83.0000 .0000 -1050.00
a 3 .0000 .0000 .0000 .0000 .0000 -2100.00
000000000000000000000
**FIN CARGA-NUDOS
******************************************************************************

ARCHIVO C.2
Archivo de salida de datos PROGES N2 = 1050 KN y N3 = 1100 KN

RESULTADOS: PASO ITERATIVO 5
DEFORMACIONES: CASO DE CARGA a
angulo phi en (rad), traslaciones en (cm) resultados en el sistema global de coordenadas
NUDO PHI-Y U W
1 .00000000 .00000000 .00000000
2 -.01112040 -14.17447199 -.49693220
3 -.03553531 -14.21412399 -.44633909
4 -.03553531 .00000000 .00000000
5 -.01676635 -1.83890550 -.12423305
6 -.02385672 -6.06042909 -.24846610
7 -.02168679 -10.76323754 -.37269915
8 -.00339253 -14.18438499 1.25299275
9 .00150283 -14.19429799 1.43539574
10 .00406738 -14.20421099 .69582770
11 -.03553531 -10.66059299 -.33475432
12 -.03553531 -7.10706199 -.22316955
97

13 -.03553531 -3.55353100 -.11158477

FUERZAS EN LOS APOYOS: CASO DE CARGA: a
*************************************************************************************
PTO COORDENADAS REACCIONES EN LOS APOYOS
NR X(M) Z(M) MIY(KNM) VIX(KN) VIZ(KN)
------------------------------------------------------------------------
1 .000 .000 1021.829 228.277 1147.913
4 10.000 4.000 .000 -73.277 2062.087
************************************************************************
FUERZAS EN LOS APOYOS: 1021.829 155.000 3210.000
**********************************************************************
FUERZAS INTERIORES EN LOS EXTREMOS DE LAS BARRAS: CASO DE CARGA a
Convencion de signos:
momentos y fuerzas como vectores en el sistema local de coordenadas las fuerzas longitudinales
de compresin son positivas
-----------------------------------------------------------------------------
BARRA MOMENTO FLECTOR FUERZA LONG. CORTANTE TENSIONES
NR MYM MYN S QZM QZN SIG I SIG K
(KNM) (KN) (KN) (KN/CM2)
-----------------------------------------------------------------------------
1 1021.83 -562.17 1147.91 -228.28 210.28 41.84 28.89
2 562.17 -111.15 1147.91 -210.28 192.28 28.89 16.18
3 111.15 309.39 1147.91 -192.28 174.28 16.18 21.76
4 -309.39 679.10 1147.91 -174.28 156.28 21.76 32.18
5 -679.10 454.35 73.28 97.91 -82.91 14.28 9.83
6 -454.35 265.95 73.28 82.91 -67.91 9.83 6.10
7 -265.95 114.37 73.28 67.91 -52.91 6.10 3.10
8 -114.37 .00 73.28 52.91 -37.91 3.10 .83
9 .00 .00 2062.09 73.28 -73.28 23.43 23.43
10 .00 .00 2062.09 73.28 -73.28 23.43 23.43
11 .00 .00 2062.09 73.28 -73.28 23.43 23.43
12 .00 .00 2062.09 73.28 -73.28 23.43 23.43
------------------------------------------------------------------------------
PASO ITERATIVO 5 PORCENTAJ E DE CAMBIO .00000
EN LA BARRA NR 9 EN EL CASO DE CARGA 1
*****************************************************************************

98

ARCHIVO C.3
Archivo salida de datos PROGES con N2 = 1700 KN y N3 = 2400 KN

==============================================================DEFORMACION
ES: CASO DE CARGA a
------------------------------
NGULO PHI EN (RAD), TRASLACIONES EN (CM)
RESULTADOS EN EL SISTEMA GLOBAL DE COORDENADAS
-------------------------------------------------
NUDO PHI-Y U W
--------------------------------------------------
1 .00000000 .00000000 .00000000
2 -.02645384 -29.40217258 -.81767418
3 .01392006 -29.53816194 -.73661745
4 -.07384540 .00000000 .00000000
5 -.03277029 -3.53973433 -.20441855
6 -.04870691 -11.98287198 -.40883709
7 -.04654417 -21.81264825 -.61325564
8 -.00971621 -29.43616992 3.62714740
9 .00301719 -29.47016726 4.37421376
10 .01109856 -29.50416460 2.50633335
11 -.07384540 -22.15362146 -.55246309
12 -.07384540 -14.76908097 -.36830873
13 -.07384540 -7.38454049 -.18415436

*****************************************************************************
------------------------------------------------------------------------
1 .000 .000 1819.004 322.341 1798.883
2 .000 8.000 .000 .000 .000
4 10.000 4.000 .000 -239.341 3241.117
5 .000 2.000 .000 .000 .000
6 .000 4.000 .000 .000 .000
99

7 .000 6.000 .000 .000 .000
******************************************************************************
==============================================================
-----------------------------------------------------------------------------
1 1819.00 -1110.65 1798.88 -322.34 322.34 67.87 49.42
2 1110.65 -314.08 1798.88 -322.34 322.34 49.42 28.68
3 314.08 507.43 1798.88 -322.34 322.34 28.68 33.71
4 -507.43 1288.64 1798.88 -322.34 322.34 33.71 54.06
5 -1288.64 1033.32 239.34 98.88 -113.88 27.90 22.91
6 -1033.32 731.65 239.34 113.88 -128.88 22.91 17.02
7 -731.65 386.22 239.34 128.88 -143.88 17.02 10.27
8 -386.22 .00 239.34 143.88 -158.88 10.27 2.72
9 .00 .00 3241.12 239.34 -239.34 36.83 36.83
10 .00 .00 3241.12 239.34 -239.34 36.83 36.83
11 .00 .00 3241.12 239.34 -239.34 36.83 36.83
12 .00 .00 3241.12 239.34 -239.34 36.83 36.83
------------------------------------------------------------------------------
******************************************************************************

ARCHIVO C.4

==============================================================
--------------------------------------------------
NUDO PHI-Y U W
--------------------------------------------------
1 .00000000 .00000000 .00000000
100

2 -.02645385 -29.40219353 -.81767424
3 .01392006 -29.53818316 -.73661742
4 -.07384546 .00000000 .00000000
5 -.03277032 -3.53973736 -.20441856
6 -.04870694 -11.98288165 -.40883712
7 -.04654420 -21.81266491 -.61325568
8 -.00971621 -29.43619094 3.62714805
9 .00301719 -29.47018835 4.37421381
10 .01109856 -29.50418576 2.50633312
11 -.07384546 -22.15363737 -.55246307
12 -.07384546 -14.76909158 -.36830871
13 -.07384546 -7.38454579 -.18415436

******************************************************************************
------------------------------------------------------------------------
1 .000 .000 1819.006 322.342 1798.883
2 .000 8.000 .000 .000 .000
4 10.000 4.000 .000 -239.342 3241.117
5 .000 2.000 .000 .000 .000
6 .000 4.000 .000 .000 .000
7 .000 6.000 .000 .000 .000
******************************************************************************

==============================================================
-----------------------------------------------------------------------------
1 1819.01 -1110.65 1798.88 -322.34 322.34 67.87 49.42
2 1110.65 -314.08 1798.88 -322.34 322.34 49.42 28.68
3 314.08 507.43 1798.88 -322.34 322.34 28.68 33.71
4 -507.43 1288.64 1798.88 -322.34 322.34 33.71 54.06
5 -1288.64 1033.32 239.34 98.88 -113.88 27.90 22.91
101

6 -1033.32 731.65 239.34 113.88 -128.88 22.91 17.02
7 -731.65 386.22 239.34 128.88 -143.88 17.02 10.27
8 -386.22 .00 239.34 143.88 -158.88 10.27 2.72
9 .00 .00 3241.12 239.34 -239.34 36.83 36.83
10 .00 .00 3241.12 239.34 -239.34 36.83 36.83
11 .00 .00 3241.12 239.34 -239.34 36.83 36.83
12 .00 .00 3241.12 239.34 -239.34 36.83 36.83
------------------------------------------------------------------------------
******************************************************************************

ARCHIVO C.5

--------------------------------------------------
NUDO PHI-Y U W
--------------------------------------------------
1 .00000000 .00000000 .00000000
2 -.03307898 -36.93112568 -.87841857
3 .01725968 -37.11074797 -.77442708
4 -.09277687 .00000000 .00000000
5 -.04117982 -4.44592927 -.21960464
6 -.06123522 -15.05958931 -.43920928
7 -.05845594 -27.41333332 -.65881393
8 -.01205709 -36.97603125 4.66403424
9 .00380772 -37.02093682 5.58004381
10 .01379718 -37.06584240 3.24997089
11 -.09277687 -27.83306098 -.58082031
12 -.09277687 -18.55537399 -.38721354
13 -.09277687 -9.27768699 -.19360677

******************************************************************************
102

------------------------------------------------------------------------
1 .000 .000 2281.904 399.135 1932.521
2 .000 8.000 .000 .000 .000
4 10.000 4.000 .000 -316.135 3407.479
5 .000 2.000 .000 .000 .000
6 .000 4.000 .000 .000 .000
7 .000 6.000 .000 .000 .000

******************************************************************************


-----------------------------------------------------------------------------
1 2281.90 -1397.71 1932.52 -399.14 399.14 81.44 58.42
2 1397.71 -394.33 1932.52 -399.14 399.14 58.42 32.29
3 394.33 642.68 1932.52 -399.14 399.14 32.29 38.76
4 -642.68 1624.88 1932.52 -399.14 399.14 38.76 64.33
5 -1624.88 1292.35 316.14 132.52 -147.52 35.35 28.85
6 -1292.35 907.69 316.14 147.52 -162.52 28.85 21.33
7 -907.69 475.27 316.14 162.52 -177.52 21.33 12.88
8 -475.27 .00 316.14 177.52 -192.52 12.88 3.59
9 .00 .00 3407.48 316.14 -316.14 38.72 38.72
10 .00 .00 3407.48 316.14 -316.14 38.72 38.72
11 .00 .00 3407.48 316.14 -316.14 38.72 38.72
12 .00 .00 3407.48 316.14 -316.14 38.72 38.72
------------------------------------------------------------------------------

******************************************************************************
103

ARCHIVO C.6

-------------------------------------------------
NUDO PHI-Y U W
--------------------------------------------------
1 .00000000 .00000000 .00000000
2 -.06391487 -71.67280935 -1.03936835
3 .03299048 -72.04669562 -.83031583
4 -.18011674 .00000000 .00000000
5 -.07985513 -8.60839027 -.25984209
6 -.11900624 -29.21922458 -.51968417
7 -.11350814 -53.22415346 -.77952626
8 -.02309677 -71.76628092 9.63794781
9 .00745375 -71.85975249 11.36532374
10 .02647833 -71.95322406 6.87232669
11 -.18011674 -54.03502172 -.62273687
12 -.18011674 -36.02334781 -.41515791
13 -.18011674 -18.01167391 -.20757896


******************************************************************************
------------------------------------------------------------------------
1 .000 .000 4402.456 741.040 2286.610
2 .000 8.000 .000 .000 .000
4 10.000 4.000 .000 -658.040 3653.390
5 .000 2.000 .000 .000 .000
6 .000 4.000 .000 .000 .000
7 .000 6.000 .000 .000 .000

******************************************************************************
104


-----------------------------------------------------------------------------
-----------------------------------------------------------------------------
1 4402.46 -2723.54 2286.61 -741.04 741.04 140.69 96.97
2 2723.54 -770.17 2286.61 -741.04 741.04 96.97 46.10
3 770.17 1260.80 2286.61 -741.04 741.04 46.10 58.88
4 -1260.80 3164.73 2286.61 -741.04 741.04 58.88 108.46
5 -3164.73 2499.71 658.04 286.61 -301.61 69.33 56.33
6 -2499.71 1738.30 658.04 301.61 -316.61 56.33 41.45
7 -1738.30 898.46 658.04 316.61 -331.61 41.45 25.04
8 -898.46 .00 658.04 331.61 -346.61 25.04 7.48
9 .00 .00 3653.39 658.04 -658.04 41.52 41.52
10 .00 .00 3653.39 658.04 -658.04 41.52 41.52
11 .00 .00 3653.39 658.04 -658.04 41.52 41.52
12 .00 .00 3653.39 658.04 -658.04 41.52 41.52
------------------------------------------------------------------------------
******************************************************************************

ARCHIVO C.7
Archivo salida de datos SAP2000 con N2 = 1050 KN y N3 = 1100 KN

SAP2000 Advanced Version 10.1.0.0 (Analysis Build 8511)
File: ...ritorio\TESIS\analisis_libro_estabilidad\ejemplo_libro_estabilidad.LOG
B E G I N A N A L Y S I S 2008/11/26 14:40:46
MAXIMUM MEMORY BLOCK SIZE (BYTES) = 63.452 MB
NUMBER OF J OINTS IN THE MODEL = 4

N O N L I N E A R S T A T I C A N A L Y S I S 14:40:46

CASE: PDELTA

105

STARTING FROM ZERO (UNSTRESSED) INITIAL CONDITIONS

LOAD CONTROL TYPE = FORCE
NUMBER OF STAGES = 0

TYPE OF GEOMETRIC NONLINEARITY = P-DELTA
INCLUDE ELASTIC MATERIAL NONLINEARITY = YES
INCLUDE INELASTIC MATERIAL NONLINEARITY = YES
METHOD TO USE WHEN HINGES DROP LOAD =UNLOAD ENTIRE STRUCTURE
SAVE POSITIVE INCREMENTS ONLY = YES

MINIMUM NUMBER OF SAVED STEPS = 10
MAXIMUM NUMBER OF SAVED STEPS = 100
MAXIMUM NUMBER OF NULL STEPS = 50
MAXIMUM NUMBER OF TOTAL STEPS = 200

MAXIMUM NUMBER OF ITERATIONS PER STEP = 15
RELATIVE FORCE CONVERGENCE TOLERANCE = 0.000100
RELATIVE EVENT TOLERANCE = 0.010000
E L E M E N T F O R M A T I O N 14:40:46

NUMBER OF FRAME ELEMENTS FORMED = 3

STEP 0
NUMBER OF EIGENVALUES BELOW SHIFT = 0
Stiffness matrix formed, number of negative eigenvalues =0

STEP 1, INCREMENT = 0.100000, SUM = 0.100000, MAX = 0.100000
Unbalance/tolerance = 3.782161, Iteration 1 did not converge!
Unbalance/tolerance = 0.185632, Iteration 2 converged
Saved as Output Step 1

106

.
.
.

TOTAL NUMBER OF CONVERGED STEPS SAVED = 10
TOTAL NUMBER OF CONVERGED STEPS NOT SAVED = 0
TOTAL NUMBER OF NULL STEPS = 0
----------
TOTAL NUMBER OF ALL STEPS = 10

TOTAL NUMBER OF ITERATIONS CONVERGED = 20
TOTAL NUMBER OF ITERATIONS DISCARDED = 0
FOR CONVERGED STEPS ONLY:
AVERAGE NUMBER OF ITERATIONS PER STEP = 2.00

TIME FOR INITIALIZING ANALYSIS = 0.03
TIME FOR CONTROLLING ANALYSIS = 0.02
TIME FOR FORMING STIFFNESS MATRIX = 0.00
TIME FOR SOLVING STIFFNESS MATRIX = 0.09
TIME FOR CALCULATING DISPLACEMENTS = 0.00
TIME FOR DETERMINING EVENTS = 0.00
TIME FOR UPDATING STATE = 0.01
----------
TOTAL TIME FOR THIS ANALYSIS = 0.15



L I N E A R E Q U A T I O N S O L U T I O N 14:40:47

FORMING STIFFNESS AT THE END OF CASE: PDELTA

TOTAL NUMBER OF EQUILIBRIUM EQUATIONS = 15
NUMBER OF NON-ZERO STIFFNESS TERMS = 102


107

L I N E A R S T A T I C C A S E S 14:40:47

USING STIFFNESS AT THE END OF CASE: PDELTA

TOTAL NUMBER OF CASES TO SOLVE = 2
NUMBER OF CASES TO SOLVE PER BLOCK = 2

LINEAR STATIC CASES TO BE SOLVED:

CASE: SC
CASE: PP
A N A L Y S I S C O M P L E T E 2008/11/26 14:40:47
******************************************************************************

ARCHIVO C.8
Archivo de salida de datos SAP2000 N2 = 1800 KN N3 = 3600 KN


B E G I N A N A L Y S I S 2008/11/26 16:41:59




CASE: PDELTA




108





STEP 0

.
.
.


109




110


----------


----------







111



CASE: SC
CASE: PP
******************************************************************************

ARCHIVO C.9
Archivo salida de datos SAP2000 N1 = 3000 KN y N2 = 6000 KN


B E G I N A N A L Y S I S 2008/11/27 04:04:41




CASE: PDELTA




112




STEP 0

.
.
.
113

Iteration limit reached - step size reduced

Iteration is diverging - step size reduced
114



Unbalance/tolerance = 1.2945E+07, Iteration 5 did not converge!
Iteration is diverging - use N-R iteration
115

Iteration is diverging - use N-R iteration

116


117



----------


----------





118


CASE: SC
CASE: PP


Anexo D Clculo de Factores B
2
y F
A

119

Anexo D
Clculo de Factores B
2
y F
A

En este anexo se muestra en detalle el clculo de los factores B
2
y F
A
para el primer valor
de la tabla

Para el AISC

Se definen los siguientes valores para las ecuaciones utilizadas:

P
nt
= Carga vertical total soportada por el piso, usando combinacin de cargas (N)
R
M
= 0,85 para marcos no arriostrados
H = Corte de piso producido por la fuerza lateral usada para calcular delta H (N)
L = Altura de piso (mm)
H
= Deformacin de primer orden entrepiso, debido a fuerzas laterales (mm)
= 1.6 factor de seguridad

Para ejemplificar se usan los siguientes valores:

Carga en nudo 1 horizontal =82.5 KN
Carga en nudo 1 vertical =100 KN
Carga en nudo 2 vertical =200 KN
qL = 60 KN

P
nt
= Carga en nudo 1 + Carga en nudo 2 + q L = 360000 N
H = 154500 N (no se considera la imperfeccin inicial)
H
= 83.2 mm (no se considera la imperfeccin inicial)
L = 8000 mm

B
2
=
1
1 -
o P
nt
P
c2
1 con P
c2
= R
M
EI
A
H

con P
c2
= u.8S
1S4Suu x 8uuu
8S.2
= 126274u4 N

B
2
=
1
1 -
1.6 x S6uuuu
126274u4
= 1.u4779SS
Anexo D Clculo de Factores B
2
y F
A

120

Para el Eurocodigo

Se definen los siguientes valores para las ecuaciones utilizadas:

= Desplazamiento horizontal entre piso producido por las fuerzas laterales mas las fuerzas
ficticias
H = Reaccin horizontal total en el piso producido por las fuerzas laterales y fuerzas ficticias
F
A
= Factor de amplificacin (semejante a B
2
)
V
sd
V
cr
= Valor elstico crtico de la carga vertical total considerando un estado no arriostrado.

Para ejemplificar se usan los siguientes valores:

Carga en nudo 1 horizontal = 82.5kN
Carga en nudo 1 vertical =100kN
Carga en nudo 2 vertical = 200kN
qL = 60 KN

V = Carga en nudo 1 + Carga en nudo 2 + q L = 360000 N
H = 155600 N (considera la imperfeccin inicial)
= 84 mm (considera la imperfeccin inicial)
h = 8000 mm

F
A
=
1
[1 -
I
sd
I
c
con o
c
=
I
sd
I
c
o
c
=

b
x
I
E
u.1

o
c
=
84
8uuu
x
S6uuuu
1SS6uu
= u.u2429S1

F
A
=
1
(1 -u.u2429S1)
= 1.u248979

Efectos de Segundo Orden 1

Cargado por

Copyright:

Formatos disponibles

Efectos de Segundo Orden 1

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Efectos de Segundo Orden 1

Cargado por

Copyright:

Formatos disponibles

i

Universidad Austral de Chile

Captulo III: Efectos de segundo orden considerados en las normas de diseo.

5.2. Modelo de la estructura en SAP2000 57

Figura 4.2.1.1. Diagrama de flujo de programa Proges 50

Figura B.7. Geometra y combinaciones de carga 92

para la matriz de rigidez 53

= Matriz de rigidez total del elemento.

= Vector desplazamiento nodal del elemento referido al sistema de referencia que

= Vector de fuerza nodal como se muestra en la figura 2.2.1 .

si no boy ponJco locol

para la matriz de rigidez

para el caso de carga axial de compresin,

son las rigideces a la flexin

También podría gustarte