TESIS PUCP

Esta obra ha sido publicada bajo la licencia Creative Commons

Reconocimiento-No comercial-Compartir bajo la misma licencia 2.5 Per.
Para ver una copia de dicha licencia, visite
http://creativecommons.org/licenses/by-nc-sa/2.5/pe/
PONTIFICIA UNIVERSIDAD CATOLICA DEL PERU
ESCUELA DE GRADUADOS

ESTUDIO DEL CONTROL DE ROTULAS PLASTICAS
EN PORTICOS DE CONCRETO ARMADO

TESIS

PARA OPTAR EL GRADO DE

MAGISTER EN INGENIERIA CIVIL

Presentado por:
TEODORICO MANUEL FLORES BARRETO

LIMA PERU
2003










A mi hijo:

KEVIN MANUEL

Con el deseo que este trabajo le sirva de estmulo para
su desarrollo personal.




ii
AGRADECIMIENTOS


Quiero expresar mi ms sincero agradecimiento al Mag. Daniel Quiun Wong, por su
asesoramiento as como por todo el apoyo brindado en el desarrollo de este estudio.

A todas las personas que de alguna u otra forma colaboraron en la realizacin de este
estudio, muchas gracias.
iii
SUMARIO

En este trabajo, se presenta el estudio del efecto de los modos altos de vibracin de
estructuras de concreto armado en la amplificacin de los momentos flexionantes en las
columnas, durante una respuesta ssmica. Se buscar cuantificar este efecto para
estructuras aporticadas mediante anlisis dinmico inelstico, y proponer valores para su
aplicacin en el diseo de prticos de concreto armado.
Se incluyen definiciones, consideraciones, procedimientos y modelos analticos para la
repuesta no lineal de estructuras aporticadas, principalmente referente al anlisis dinmico
inelstico.
El modelo utilizado corresponde a 19 estructuras aporticadas de 2 a 20 pisos con igual
planta rectangular en todos los pisos, estructurado en base a cuatro prticos resistentes de
cuatro vanos en el sentido longitudinal. Como excitacin ssmica se utilizan las dos
componentes de tres registros de sismos, ocurridos en 1966, 1970 y 1974 en Per,
normalizadas a una aceleracin mxima igual a 0.4g. Para representar el comportamiento
del material se utiliza un modelo bilineal con una prdida de rigidez de 95% en el segundo
tramo, y no se considera prdida de las propiedades de rigidez como consecuencia de los
fenmenos de carga-descarga en las distintas secciones del modelo. Las respuestas
ssmicas se obtuvieron a travs de un anlisis dinmico inelstico, para el cual se utiliz el
programa de cmputo DRAIN-2D (15, 24), registrndose la historia de los momentos en el
tiempo, para todos los elementos del modelo.
Con los resultados obtenidos, se calculan los Factores de Distribucin de Momentos
(FDM) en los nudos interiores y exteriores de las estructuras aporticadas, se comparan
dichos valores con los requerimientos de algunas normas de diseo; se realiza un
tratamiento estadstico de la dispersin de los resultados y se proponen valores en funcin
del perodo fundamental del prtico, como una propuesta para ser incorporados en la NTE
E-060 (27), para su aplicacin en el diseo de prticos resistentes a momento, con el
objetivo de tratar de controlar o retardar la formacin de rtulas plsticas en las columnas.
iv
INDICE
DEDICATORIA i
AGRADECIMIENTOS ii
SUMARIO iii
INDICE iv

CAPITULO I: INTRODUCCION 1
1.1 INTRODUCCION 1
1.2 OBJETIVOS 3
1.3 ORGANIZACIN DE LA TESIS 3

CAPITULO II: CONCEPTOS GENERALES PARA EL ANALISIS INELASTICO 5
2.1 INTRODUCCION 5
2.2 ANALISIS DE LA RESPUESTA INELASTICA DE ESTRUCTURAS 6
2.3 MODELOS DE ANALISIS INELASTICO 8
2.3.1 MODELOS SIMPLES 9
2.3.2 MODELOS DISCRETOS 11
2.3.3 MODELOS POR ELEMENTOS FINITOS 14
2.3.4 MODELOS DE FIBRA 14
2.4 MODELOS HISTERETICOS 16
2.4.1 MODELO BILINEAL 16
2.4.2 MODELO DE TAKEDA 17
2.4.3 MODELO DE SINA 18
2.4.4 MODELO DE OTANI 19
2.4.5 MODELO Q-HYST 20
2.5 MODELOS DE ELEMENTOS INELASTICOS 21
2.5.1 MODELO DE DOBLE COMPONENTE 21
2.5.2 MODELO DE UNA COMPONENTE 22
2.5.3 MODELO PARABOLICO 23
2.5.4 MODELO DEL ELEMENTO DISCRETO 23
v

CAPITULO III: MODELO ANALITICO 25
3.1 INTRODUCCION 25
3.2 SUPOSICIONES RESPECTO A LA ESTRUCTURA Y MOVIMIENTO
EN LA BASE 26
3.3 METODO DE RIGIDEZ DIRECTA 26
3.3.1 GRADOS DE LIBERTAD DE LA ESTRUCTURA 26
3.3.2 MATRIZ DE RIGIDEZ 27
a) Rigidez del Elemento 27
b) Ensamblaje de la Rigidez de la Estructura 28
3.3.3 VECTOR DE CARGA 28
3.3.4 ECUACIONES DE EQUILIBRIO Y SOLUCION 29
3.4 ANALISIS PASO A PASO 30
3.4.1 EQUILIBRIO DINAMICO 30
3.4.2 MATRIZ Y COEFICIENTES DE AMORTIGUAMIENTO 31
3.4.3 INTEGRACION PASO A PASO 33
3.4.4 CORRECCION DEL EQUILIBRIO DEBIDO AL AMORTIGUAMIENTO
DEPENDIENTE DE LA RIGIDEZ 35
3.4.5 RIGIDEZ TANGENCIAL 37
a) Rigidez Elstica 37
b) Rigidez Geomtrica 37
3.4.6 DETERMINACION DEL ESTADO 37
a) Incremento de Fuerza en el Elemento 37
b) Incremento de Deformacin del Elemento 41
c) Desbalance en el Equilibrio 41
3.4.7 MOVIMIENTO SISMICO 43
3.5 ELEMENTOS ESTRUCTURALES 43
3.5.1 ELEMENTO VIGA-COLUMNA 44
a) Caractersticas Generales 44
b) Deformaciones del Elemento 45
vi
c) Superficies de Interaccin 47
d) Rigidez del Elemento 50
e) Fuerzas de Fijacin en los Extremos e Iniciales 51
f) Reduccin de la Carga Viva 52
g) Deformaciones por Cortante 53
h) Excentricidad en los Extremos (nudo rgido) 54
i) Salida de Resultados 55
3.5.1 ELEMENTO VIGA CON RIGIDEZ DEGRADANTE 56
a) Caractersticas Generales 56
b) Deformaciones del Elemento 63
c) Propiedades del Elemento 63
d) Rigidez del Elemento 67
e) Salida de Resultados 68

CAPITULO IV: ESTUDIO ANALITICO USANDO DRAIN-2D 70
4.1 INTRODUCCION 70
4.2 ESTRUCTURAS EN ESTUDIO 70
4.3 PREDIMENSIONAMIENTO DE LOS PORTICOS EN ESTUDIO 78
4.4 ANALISIS ELASTICO DE LOS PORTICOS EN ESTUDIO 79
4.5 DISEO DE LOS PORTICOS EN ESTUDIO 81
4.6 PROPIEDADES DE FLEXION DE LOS ELEMENTOS
82
4.6.1 RELACIONES ESFUERZO-DEFORMACION PARA EL CONCRETO
Y EL ACERO DE REFUERZO 85
4.7 REGISTROS DE ACELERACION CONSIDERADOS 93
4.8 AMORTIGUAMIENTO 98
4.9 ANALISIS INELASTICO DE LOS PORTICOS 99
4.10 RESULTADOS DE LOS ANALISIS INELASTICOS 101


vii
CAPITULO V: DESARROLLO DEL FACTOR DE DISTRIBUCION DE MOMENTOS 102
5.1 INTRODUCCION 102
5.2 CONTROL DE ROTULAS PLASTICAS 104
5.3 REVISION DE REGLAMENTOS RESPECTO AL CONTROL DE ROTULAS
PLASTICAS 106
5.3.1 NORMA NTE-060 CONCRETO ARMADO 106
5.3.2 NORMA ACI 318-99 106
5.3.3 REGLAMENTO DE CONSTRUCCIONES PARA EL DISTRITO
FEDERAL-MEXICO 107
5.4 FACTOR DE DISTRIBUCION DE MOMENTOS (FDM) 108
5.4.1 RESULTADOS TIEMPO HISTORIA DE MOMENTOS FLEXIONANTES 109
5.4.2 CALCULO DEL FACTOR DE DISTRIBUCION DE MOMENTOS 111
5.4.3 SELECCION DEL FACTOR DE DISTRIBUCION DE MOMENTOS 115
5.5 ANALISIS E INTERPRETACION DE LOS RESULTADOS 115

CAPITULO VI: CONCLUSIONES Y RECOMENDACIONES 122
6.1 CONCLUSIONES 122
6.2 RECOMENDACIONES 124

ANEXO A
MANUAL DE USO DEL PROGRAMA DRAIN-2D PARA PC 125

ANEXO B
EJEMPLOS DE ARCHIVOS DE DATOS 168

BIBLIOGRAFIA 172







CAPITULO I


INTRODUCCION


1.1 INTRODUCCION
Al ser sometida una estructura de concreto armado a movimientos ssmicos severos,
sta generalmente responde no linealmente. Esto es atribuible a que el concreto armado es
un material no homogneo y su comportamiento es altamente no lineal.
El nmero de variables comprendidas en la respuesta no lineal de estructuras de varios
niveles es tan elevado, que estudios anteriores han demostrado la alta dependencia de las
respuestas inelsticas a las caractersticas propias de cada modelo, y el estado del arte
actual hace difcil el poder dar recomendaciones de carcter general (11). Sin embargo, se
pueden obtener valores cualitativos de ciertos factores dentro de un intervalo probable, para
su posterior utilizacin dentro de un procedimiento racional de diseo.
Los anlisis dinmicos no lineales indican que en las columnas de prticos de varios
pisos, pueden ocurrir distribuciones inesperadas de momento flexionante (21), en
comparacin con la distribucin obtenida de la carga lateral esttica equivalente que
2
recomienda la Norma Tcnica de Edificacin E-030 de Diseo Sismorresistente (26). El
anlisis de carga lateral esttica indica generalmente que los puntos de inflexin se
encuentran prximos a la mitad de la altura de las columnas, excepto en los pisos prximos
a la parte superior e inferior del prtico. Sin embargo, el anlisis dinmico no lineal sugiere
que en determinados instantes durante la respuesta de la estructura a los movimientos
ssmicos, el punto de inflexin en una columna puede estar prximo a la unin viga-columna
y ocasionalmente, incluso la columna puede estar en curvatura simple.
La causa de la distribucin inesperada de los momentos flexionantes de la columna en
algunos instantes, es la fuerte influencia de los modos superiores de vibracin,
especialmente el segundo y el tercero (21). Por otro lado, Hejal y Chopra (12) realizaron un
estudio paramtrico de la respuesta ssmica de una estructura de varios pisos, y concluyen
estudiando la contribucin modal, que manteniendo los dems parmetros fijos, la
contribucin de los modos altos a la respuesta total en general, aumenta al aumentar el
perodo fundamental; as mismo observaron en las respuestas locales, que la contribucin
de los modos ms altos es ms importante en los momentos de las columnas que en los
momentos de las vigas y esfuerzos axiales de las columnas.
La carga lateral esttica de la Norma E-030 normalmente tiene una distribucin
triangular que corresponde a las cargas laterales que varan linealmente desde cero en la
base a un mximo en la parte superior de la estructura. Esta distribucin de cargas
corresponde aproximadamente al modo fundamental (primero) de vibracin; si son
significativos los otros modos de vibracin es evidente que la distribucin de las cargas en
la estructura en determinados instantes podra diferir apreciablemente de la supuesta y
producir la distribucin de momentos flexionantes antes mencionadas. El corrimiento del
punto de inflexin en las columnas a un punto bien distante de la posicin mitad de la
altura, en algunos casos significa que los momentos inducidos en las columnas pueden ser
mucho mayores que los momentos obtenidos de un anlisis de carga lateral esttica y
pueden conducir a que se formen rtulas plsticas en la columnas.
Los anlisis dinmicos no lineales tambin indican que generalmente no estn
presentes las rtulas plsticas en todas las vigas en el mismo intervalo de tiempo. El
desarrollo de rtulas plsticas tiende a moverse hacia arriba del prtico, en ondas que
3
abarcan unos cuantos pisos a la vez, pero puede haber instantes, en edificios de baja altura
donde todas las vigas tengan rtulas plsticas formadas simultneamente.

1.2 OBJETIVOS
En el presente trabajo se busca encontrar valores cualitativos del Factor de Distribucin
de Momentos (FDM) dentro de un rango probable para estructuras aporticadas de concreto
armado, y proponer valores para su aplicacin en el diseo, con el propsito de controlar la
formacin de rtulas plsticas en las columnas de una conexin viga-columna. Luego, se
presenta una propuesta para que estos valores puedan ser incorporados en la NTE-060
Concreto Armado del RNC (27), como una disposicin especial para columnas sujetas a
flexocompresin que resistan fuerzas de sismo.

1.3 ORGANIZACIN DE LA TESIS
El presente estudio se ha organizado de la siguiente manera:

En el captulo II, se presentan las definiciones generales en los que se basa el anlisis
dinmico inelstico, describiendo los diferentes modelos de anlisis, modelos histerticos y
modelos de elementos inelsticos.

En el captulo III, se describen las consideraciones y parmetros involucrados en el
anlisis dinmico inelstico aplicado a estructuras aporticadas, mediante el modelo analtico
de rtula puntual, poniendo nfasis en el procedimiento de anlisis y elementos
estructurales que usa el programa de cmputo usado.

En el captulo IV, se realizan los anlisis dinmicos inelsticos usando el programa de
cmputo DRAIN-2D, se describen las estructuras que se analizan, las que han sido
diseadas previamente segn el actual cdigo peruano sismorresistente; adems se
comentan los registros ssmicos usados, as como la repuesta inelstica de estas
estructuras.

4
En el captulo V, se presenta el procedimiento para calcular el factor de distribucin de
momentos (FDM), en base a los resultados obtenidos de los anlisis inelsticos de las
estructuras seleccionadas. Los resultados se comparan con los requerimientos de algunas
normas y se proponen nuevos valores, para su aplicacin en el diseo sismorresistente.

Las conclusiones de esta investigacin, as como las recomendaciones para su
ampliacin, se presentan en captulo VI.

En el anexo A se presenta el manual de uso del programa de cmputo DRAIN-2D,
versin para computadoras personales y en Anexo B se indica un ejemplo de archivo para
el ingreso de datos al programa.








CAPITULO II


CONCEPTOS GENERALES PARA
EL ANALISIS INELASTICO


2.1 INTRODUCCION
En la filosofa del diseo sismo-resistente de edificaciones, se supone la disipacin de
energa ssmica a travs del comportamiento inelstico de las estructuras. Esta
consideracin permite obtener un diseo econmico de la estructura, sin reducir su
seguridad. Esta seguridad depender realmente de que los mecanismos que se formen
como resultado del comportamiento inelstico, sean estables.
Con la finalidad de cuantificar este comportamiento, y conocer la capacidad o seguridad
de la estructura, diversos investigadores han propuesto mtodos y modelos para predecir la
respuesta inelstica. En este tipo de anlisis entran en juego muchas variables y
parmetros, que necesitan ser definidos y/o supuestos. Algunas de estas variables son
aleatorias tal como la seleccin del sismo; una forma de tomar en cuenta este elemento
6
importante, sera escogiendo sismos representativos del rea y si fuera necesario,
escalarlos a mxima aceleracin o velocidad (caso Japons). Con el desarrollo de las
investigaciones experimentales, los modelos tericos de anlisis cada vez representan
mejor la estructura real posibilitando mejores estimaciones de la respuesta.

2.2 ANALISIS DE LA RESPUESTA INELASTICA DE ESTRUCTURAS
El comportamiento inelstico de las estructuras es tan importante que es indispensable
tomarse en cuenta, en la prctica de un diseo ssmico por las siguientes razones (16): la
estructura de un edificio debe comportarse sin experimentar daos bajo sismos pequeos o
medianos que pueden ocurrir durante su existencia. Adems, no debe sufrir un colapso con
un fuerte movimiento ssmico que tenga recurrencias de 50 aos o ms. A menudo, las
estructuras diseadas con esta filosofa estn sujetas a fuerzas ssmicas que las llevan al
rango inelstico. En algunas ocasiones, las fuerzas observadas han sido de tres a cuatro
veces mayores que las que se especifican en los reglamentos. A pesar de ello, en la
mayora de los casos las estructuras no resultaron daadas. Se cree que la disipacin de
energa debida al amortiguamiento histertico es un margen adicional de seguridad que
poseen estas estructuras. Disear estructuras que permanezcan elsticas bajo grandes
movimientos ssmicos es muy costoso y se considera poco realista. El efecto de la
disipacin de energa que causa el comportamiento histertico de la estructura de un
edificio tendr, por consiguiente, que evaluarse con precisin, partiendo de un anlisis
inelstico de la estructura.
El anlisis elstico con el uso del concepto del amortiguamiento viscoso equivalente es
una manera de evaluar este efecto. Este procedimiento analtico est propenso a errores
grandes si el factor de ductilidad tambin es grande (16).
La relacin entre la respuesta elstica mxima Vy (la respuesta del sistema cuando
permanece elstico independientemente de la intensidad de la fuerza) y la respuesta
mxima inelstica Vu ha sido determinada (21) para un sistema de una sola masa, elstico
y perfectamente plstico, basndose en dos suposiciones: primero considerando que las
deflexiones mximas elsticas como elastoplsticas sean iguales, llegando a la relacin 2.1:
7

=
1
Vu
Vy
(2.1)
Donde , es el factor de ductilidad de desplazamiento, definido en 2.2 como

y
u

= (2.2)
u es la deflexin lateral en el extremo del intervalo posterior al elstico y y es la deflexin
lateral la primera vez que se alcanza la fluencia.
Otra manera es considerando el concepto de energa igual, lo que implica que la
energa potencial almacenada por el sistema elstico en la deflexin mxima es la misma
que la almacenada por el sistema elastoplstico en la deflexin mxima, llegando a la
relacin 2.3:
1 2
1
Vu
Vy

= (2.3)
Ambos conceptos descritos anteriormente, se muestran en la figura 2.1

Figura 2.1 RELACION ENTRE LAS REPUESTAS ELASTICAS E INELASTICAS
Vu
Vy
C
a
r
g
a

l
a
t
e
r
a
l

d
e

i
n
e
r
c
i
a
C
a
r
g
a

l
a
t
e
r
a
l

d
e

i
n
e
r
c
i
a
Deflexin lateral
y
Vy
u
Vu
Deflexin lateral
y u e
(a) Respuesta de deflexin mxima igual (b) Respuesta de energa potencial mxima igual
La prediccin mediante las ecuaciones 2.1 y 2.3 de una deflexin inelstica es debatible
para estructuras de concreto armado cuando sus curvas de histresis son del tipo
degradante, no obstante estas ecuaciones por simplicidad frecuentemente se utilizan en la
prctica del diseo.
8
El procedimiento anterior puede ser extensivo a los sistemas de mltiples grados de
libertad. Si el factor de ductilidad aparente definido como se indica en 2.4:
y
e
eq

= (2.4)
es constante para todos los niveles de entrepiso, la relacin para los sistemas de un grado
de libertad (ecuaciones 2.1 y 2.3) se pueden emplear para cada entrepiso en forma
aproximada. Sin embargo, si
eq
en un entrepiso es mucho mayor que el de los otros
entrepisos, la deflexin plstica se concentra en ese entrepiso. Al procedimiento para
representar cada entrepiso por un sistema de un grado de libertad se le llama mtodo de
ductilidad (16) y frecuentemente se adopta en los reglamentos de diseo.
De lo visto anteriormente, se puede decir que las estructuras que pueden sostener
grandes deformaciones plsticas tienen una buena resistencia ssmica. Por consiguiente,
las estructuras que tienen valores altos de se pueden disear con niveles ms bajos de
capacidad resistente a las fuerzas laterales. Tambin se puede establecer que la parte de
la estructura que vaya a sufrir grandes deformaciones plsticas debe poseer un factor de
ductilidad que cubra sobradamente la esperada deformacin plstica.

2.3 MODELOS DE ANALISIS INELASTICO
La respuesta dinmica de un prtico plano de varios pisos y varias luces es bastante
compleja. Existen un nmero elevado de grados de libertad y un alto porcentaje de
comportamiento no lineal, por lo que se requiere de algunas idealizaciones y
simplificaciones a fin de obtener un modelo matemtico que pueda ser resuelto empleando
las diferentes tcnicas de computacin numrica disponible. El grado de sofisticacin del
modelo es funcin del nivel de la respuesta que se desee obtener, lo cual ser siempre
relativo.
Se han propuesto diferentes aproximaciones con el fin de modelar las estructuras de
concreto armado para llevar a cabo un anlisis no lineal. Estos pueden clasificarse en los
siguientes grandes grupos (14,23): modelos simples, modelos discretos, modelos de fibra y
modelos de elementos finitos.

9
2.3.1 MODELOS SIMPLES
Muchos de los modelos desarrollados en esta categora son conocidos como modelos
tipo cortante. Inicialmente desarrollados para sistemas de un grado de libertad, su empleo
ha sido extendido a sistemas de varios grados de libertad.
En el modelo tipo cortante, se sustituye el ensamblaje de los elementos (vigas,
columnas y/o muros de corte) que constituyen el piso de un prtico por un resorte no lineal
nico (figura 2.2).

Para prticos planos se considera un solo grado de libertad por piso llamado
desplazamiento lateral de piso. En estructuras tridimensionales el sistema equivalente
posee tres grados de libertad por piso (dos desplazamientos y una rotacin). Se considera
que la losa es infinitamente rgida, lo que permite concentrar las masas de la estructura en
cada nivel de piso. El modelo cortante es tambin denominado de acoplamiento cercano ya
que el comportamiento de un piso es influido slo por los dos adyacentes, superior e
inferior. Otros modelos corresponden al resorte de flexin usado para estructuras con
muros de corte y el modelo de corte-flexin el cual acopla cinco niveles por cada piso, dos
superiores y dos inferiores.
Ki
a) ESTRUCTURA REAL
b) MODELO CORTANTE
Figura 2.2 MODELOS SIMPLES
Mi
c) MODELO 1 GDL
Me
Ke
Una gran variedad de relaciones fuerza-desplazamiento han sido usadas para estos
modelos, las cuales se determinan segn el tipo de prtico de los elementos que son
parte de l. Los estudios realizados incluyen resortes tipo elastoplstico, bilineal, trilineal y
10
otros donde se considera alguna caracterstica de degradacin de rigidez, resistencia o
ambas.
Otro modelo dentro de esta categora viene dado por el sistema equivalente de un grado
de libertad, en el cual se supone que la estructura se deforma de acuerdo a su primer modo
o alguna otra forma de deflexin predeterminada. Takizawa (30) propone un modelo donde
la estructura es sustituida por una superposicin de colapsos mecnicos, que corresponden
a diferentes estados de plastificacin en el edificio. Para el caso de un mecanismo de un
grado de libertad, al considerar la deformacin estructural total, el sistema es reducido a un
sistema inelstico equivalente de un grado de libertad. Otro procedimiento es el modelo
llamado Q (28), que consiste en una barra rgida de peso despreciable con una rtula y un
resorte rotacional no lineal en su base, y una masa concentrada y amortiguamiento viscoso
en la parte superior tal como se muestra en la figura 2.3. Una aproximacin similar usada
por Piqu (23) superpone un juego de resortes en paralelo con diferentes modelos de
fluencia. En este modelo, la rigidez combinada cambia continuamente, aunque la funcin
de histresis para cada resorte puede ser una relacin simple.

Figura 2.3 EL MODELO Q
BARRA RIGIDA
SIN MASA
A
L
T
U
R
A

E
Q
U
I
V
A
L
E
N
T
E


(
L
e
q
)
MASA EQUIVALENTE (Me)
Xg
K
C
X
El uso de estos modelos est limitado al diseo preliminar y/o para llevar a cabo
estudios de parmetros analizando el comportamiento global de la estructura. La
11
informacin obtenida usando modelos simples es insuficiente para diseo final, evaluacin
de detalles y estudios de confiabilidad. Para estos propsitos al igual que para el anlisis
de estructuras generales que no pueden ser representadas adecuadamente por modelos
simples, se necesitan tcnicas de modelamiento ms sofisticados. Sin embargo, cabe decir
que el uso de estos modelos permite conocer el comportamiento global de la estructura en
un menor tiempo y a un menor costo.

2.3.2 MODELOS DISCRETOS
Son tambin denominados modelos de rtula puntual. El anlisis dinmico inelstico
usando estos modelos es llevado a cabo a nivel del elemento integrante de la estructura.
Esto es, la estructura (un edificio aporticado, usualmente) es discretizado en elementos
prismticos: vigas, columnas y placas (figura 2.4). En algunas formulaciones los nudos son
tratados como elementos separados (30). Un gran nmero de variantes ha sido
desarrollado, trabajando con prticos planos y otras aproximaciones para resolver
estructuras tridimensionales.

Figura 2.4 MODELO DEL PORTICO BIDIMENSIONAL
La denominacin general de rtula plstica de este tipo de modelo deriva del efecto que
produce la fluencia en los elementos, se supone que la rtula se forma en el punto donde la
capacidad fue excedida y no sobre la longitud continua del elemento. Aziz (1) en su
formulacin considera una longitud plastificada que es tomada en funcin de la longitud
original del elemento, esto da una medida de la incursin plstica de los elementos. Otras
12
mejoras han sido introducidas tales como, considerar miembros con zonas rgidas en los
extremos, modelar extremos de vigas con muros de corte y tambin en elementos de barra
del tipo armadura. En todos los casos, la fluencia ocurre en una seccin dada o en un
punto a lo largo del eje longitudinal del elemento.
El criterio de fluencia de este modelo es debido a la flexin o a una curva de interaccin
de fluencia. En el primer caso la fluencia ocurre cuando el elemento excede la capacidad
de momento plstico Mp. Esto puede ocurrir para alguna direccin del momento (positivo o
negativo). En algunas formulaciones este chequeo es hecho en el extremo de los
elementos, es posible subdividiendo la viga en varios tramos considerar la formacin de
rtulas a lo largo del elemento, aunque se debera formar siempre en puntos
preespecificados. El criterio de interaccin de fluencia define una curva de fluencia (para
miembros prismticos planos) como una funcin de la capacidad ltima de fuerza axial y
momento. La forma en que esto ocurre vara para cada miembro de acero y de concreto,
pero en uno u otro caso la fluencia es alcanzada cuando una combinacin de momento y
carga axial definen un punto en o sobre la lnea de capacidad de fluencia del elemento.
La idealizacin de los elementos puede realizarse de varias maneras:
El modelo de una componente, planteado por Giberson (10) considera que el elemento
est compuesto por una componente elstica con dos resortes inelsticos en los extremos
de ste. Estos resortes pueden tener cualquier caracterstica fuerza deformacin,
habindose desarrollado y propuesto varias relaciones, desde la ms simple
correspondiente a un diagrama elastoplstico, bilineales y trilineales incluyendo
degradacin de rigidez y/o resistencia, y otras considerando el efecto pinching
(endurecimiento) en la zona cercana al origen cuando hay cambio en la direccin de la
rotacin.
Los modelos de varias componentes en paralelo, cuya primera versin fue propuesta
por Clough, corresponden al modelo de doble componente, donde el elemento se supone
que est constituido por dos componentes, uno elstico y el otro elastoplstico que actan
en forma simultnea. Otras aproximaciones intentan modelar diagramas trilineales usando
un elemento elstico y dos elementos elastoplsticos actuando en paralelo.
13
Los modelos de varios resortes en serie, consideran el efecto de grandes
deformaciones como resultado de cargas transversales sobre la longitud del elemento y
tambin toman en cuenta el hecho que la rigidez vara a lo largo del elemento incluso en el
rango elstico. El elemento se divide en varios subelementos a lo largo del eje (figura 2.5)
conectados por resortes no lineales.
Otro camino para tomar en cuenta el efecto de rigidez variable a lo largo de la longitud
del elemento es usando funciones de interpolacin, por ejemplo Umemura (30) plantea una
distribucin parablica de flexibilidades (1/EI) a lo largo del eje del elemento. Conocida sta
a travs de la interpolacin se obtiene la matriz de rigidez requerida para el anlisis.
b) Distribucion de curvatura supuesta
a) Modelo de varios resortes
Figura 2.5 MODELO DE VARIOS RESORTES PARA MUROS DE CORTES
Los modelos discretos pueden dar detalles e informacin ms precisa sobre el
comportamiento de estructuras de concreto armado que los modelos simples. Sin embargo,
modelos discretos basados en estas simplificaciones pueden ser incapaces de simular el
comportamiento complejo de elementos de concreto armado bajo condiciones arbitrarias de
carga. El efecto de que el comportamiento histertico es predefinido reduce
significativamente el potencial del modelo para adaptarse a condiciones generales de carga.
Las propiedades de los resortes son funcin de la carga particular desarrollada en el
elemento, lo cual no es conocido generalmente a priori. Hay que agregar el efecto que las
deformaciones axiales inelsticas por flexin y las deformaciones por corte pueden no estar
concentradas en el extremo, y puede ser difcil la definicin de un juego apropiado de
resortes que gobiernen la interaccin de deformacin axial y por corte.
14

2.3.3 MODELOS POR ELEMENTOS FINITOS
Los modelos por elementos finitos idealizan cada elemento estructural como un
ensamblaje de un alto nmero de elementos finitos. Estos pueden ser de varios tipos tales
como bielas, elementos tipo viga, elementos para esfuerzo plano y/o deformacin plana o
incluso elementos tridimensionales.
Los anlisis llevados a cabo (14) fueron particularmente para muros de corte. La
correlacin con resultados experimentales es buena, pero el costo computacional es
elevado debido a la gran cantidad de elementos y al comportamiento no lineal,
especialmente cuando las cargas alcanzan niveles significantes. Tambin hay que
considerar que la no-linealidad del material monitoreando la propagacin de grietas y el
deterioro de la unin en los elementos complica el anlisis. El costo y tiempo son muy altos
para el caso de estructuras grandes y particularmente al llevar a cabo anlisis dinmico.
Sin embargo, existen ciertas dudas sobre el comportamiento del material y en la
formulacin de rigideces de los elementos finitos que conlleva a resultados aproximados a
pesar del refinamiento aparente de estos modelos.

2.3.4 MODELOS DE FIBRA
Estos modelos que fueron desarrollados inicialmente para estructuras de acero, tambin
han sido extendidos a modelos estructurales de concreto armado.
La estructura es primeramente discretizada al nivel del elemento como en el modelo de
rtula puntual, entonces cada componente, viga o columna, es dividida en varios tramos
(segmentos), cada cual es compuesto por una serie de franjas (fibras) paralelas (figura 2.6).
15
La rigidez de cada elemento se obtiene por integracin a lo largo del elemento de los
segmentos cuyas propiedades fueran calculadas al nivel de seccin. Las caractersticas de
la estructura total son entonces calculadas por el proceso tradicional, el cual se basa en las
propiedades de cada elemento.
Figura 2.6 MODELO DE FIBRA
nudo J
a) Elemento dividido en segmentos
nudo I
s
e
g
m
e
n
t
o
zona rigida
(opcional)
b) Seccin de un segmento tpico
fibras
Se supone que se conocen las caractersticas fuerza-deformacin para cada fibra
(figura 2.7) y las deformaciones como los esfuerzos en la estructura, son calculadas a
travs de un anlisis paso a paso. La fluencia de las fibras es aqu precisamente un
fenmeno definido y permite mantener un valor de la tendencia de plasticidad a travs de
un segmento y a lo largo de la longitud del elemento.
La suposicin comn que las secciones planas permanece plana puede no ser
verdaderamente exacta, especialmente cuando hay un porcentaje significativo de la
deformacin por corte, grieta de corte o prdida de adherencia entre el acero y el concreto.
Figura 2.7 PROPIEDADES DEL MATERIAL DE CADA COMPONENTE
a) Modelo para el acero de refuerzo
2
1
E1
E2
E3
3
b) Modelo para el concreto
1
E1
2
E2
E3
Deformacin
Esfuerzo Esfuerzo
Deformacin
efecto del
factor de
descarga
0
.1
0
.
5
1
E1
16
Muchos modelos de fibra no consideran este efecto. Algunas simplificaciones se realizan
para deducir la rigidez del elemento a partir de las propiedades de las secciones. Se
integra sobre la longitud del elemento en un nmero de secciones escogidas usando una
funcin de interpolacin predeterminada. Esta funcin de interpolacin puede no reflejar
exactamente la variacin de rigidez a lo largo del elemento ya que esta variacin cambia
como las condiciones de carga evolucionan.

2.4 MODELOS HISTERETICOS
Para estudiar la respuesta inelstica de un sistema discreto de masas, se debe
establecer un modelo matemtico de las caractersticas de la fuerza de restitucin y de aqu
definir la relacin entre la fuerza cortante en el entrepiso y la deflexin del mismo.
Para una secuencia progresiva de cargas y descargas, la lnea que une a los puntos
pico en la curva carga-deflexin de cada secuencia de carga se llama curva esqueletal. En
muchos casos, la curva esqueletal coincide con la curva de carga montona. La curva que
se obtuvo bajo las inversiones de signos de la fuerza se llama curva de histresis. La curva
de histresis se afecta significativamente por los materiales y el tipo estructural. La mayora
de los modelos matemticos se simplifica de acuerdo con el nivel requerido de anlisis.
Estos modelos deben proporcionar la rigidez y resistencia del miembro, los cuales
varan en cada instante de tiempo con la historia de cargas y deformaciones producidas por
el movimiento ssmico.
A continuacin, se describen algunos de los modelos histerticos (4, 20) que se han
desarrollado para representar el comportamiento dominante por flexin de elementos de
concreto armado durante cargas cclicas.

2.4.1 MODELO BILINEAL
Debido a su simplicidad, el sistema histertico bilineal ha sido usado extensamente para
estructuras de acero y de concreto armado. El modelo puede ser descrito mediante slo
tres reglas, y solamente se consideran dos rigideces en el modelo: la rigidez elstica y la de
fluencia. Las pendientes de descarga y de la carga en reversa, es la misma de la etapa
elstica.
17
La observacin general en este modelo es que: la disipacin de energa es grande para
deformaciones de amplitudes altas, y para amplitudes bajas no se considera disipacin de
energa histertica. En la figura 2.8, se muestra este tipo de modelo.

2.4.2 MODELO DE TAKEDA
Takeda, Sozen y Nielsen, propusieron un modelo ms complicado basado en la
observacin experimental. Este modelo usa una curva primaria trilineal, simtrica con
relacin al origen. La curva de carga (figura 2.9) bsicamente est dirigida hacia el mximo
punto alcanzado anteriormente en esa misma direccin. La pendiente de la curva de
descarga se degrada dependiendo de la deflexin mxima alcanzada anteriormente en
cualquier direccin, segn una funcin experimental. La rigidez degradada de descarga se
expresa como:
Figura 2.8 MODELO HISTERETICO BILINEAL SIMPLE
Fuerza
0
Y '
Deformacin
Y
Curva primaria


+
+
=
Y
m
Y C
Y C
r
D
D
D D
F F
K
(2.5)

Donde C, Y indican los niveles de agrietamiento y fluencia respectivamente, es el
parmetro de degradacin de rigidez a la descarga. El modelo de Takeda, se aplica a
elementos donde la falla es predominantemente por flexin. Este modelo ha sido usado en
varias investigaciones con resultados satisfactorios, pero tiene el inconveniente de ser
18
bastante complejo en su forma original por lo que tambin se han desarrollado modelos
simplificados de Takeda.

Figura 2.9 MODELO HISTERETICO DE TAKEDA
Fuerza
0
U'm
Y '
C'
C
Deformacin
Y
Um
Curva primaria
2.4.3 MODELO DE SINA
El modelo Sina (figura 2.10) es una versin del modelo de Takeda, al cual se ha
aadido el efecto de pinching (tendencia a presentarse una rigidez incremental muy
pequea cerca al origen seguido por un aumento brusco de la rigidez, lo que se observa
frecuentemente en experimentos). Esta consideracin mejora la respuesta durante
perodos de amplitud pequea, y as mismo produce mayores desplazamientos en
comparacin con otros modelos.
La curva primaria consiste de tres partes similares a los usados por Takeda. La
pendiente de descarga para la regin pos fluencia se asume como:

=
m
Y
Y C
D
D
S S
' 1
(2.6)

Donde S
CY
es la pendiente de una lnea que conecta el punto de fluencia con el punto
de agrietamiento en la direccin opuesta, D
Y
es la deformacin de fluencia, D
m
es la
19
deformacin mxima obtenida en la direccin de la carga, es una constante que se asume
igual a 0.5.

Figura 2.10 MODELO HISTERETICO DE SINA
Fuerza
0
U'm
Y '
C'
X1
C
Xo
Deformacin
B
Y
Um
Curva primaria
2.4.4 MODELO DE OTANI
Este modelo es una versin modificada y simplificada del modelo de Takeda.
Originalmente fue propuesto para representar la variacin de la rigidez de un resorte de
nudo (que considera el efecto de adherencia entre el concreto y el refuerzo) y la rigidez de
un resorte flexionante. La curva primaria (figura 2.11) de este modelo es bilineal con el
punto de quiebre en la zona de fluencia de la seccin. Tambin considera una degradacin
de la rigidez durante la etapa de descarga. El tratamiento de los ciclos de amplitud
pequea es similar al modelo de Takeda, por lo que puede considerarse que este modelo
es todava complicado.
20

Figura 2.11 MODELO HISTERETICO DE OTANI
Fuerza
0
U'm
Y '
X1
Xo
Deformacin
Y
Um
Curva primaria
2.4.5 MODELO Q-HYST
El modelo Q-Hyst (figura 2.12), fue desarrollado por Saiidi y Sozen (28) y puede ser
considerado como una modificacin del sistema histertico bilineal. El propsito bsico de
la modificacin fue darle ramificaciones suavizadas para las etapas de descarga y de carga
en reversa.
La descarga desde un punto ms all del punto de fluencia y la recarga en la otra
direccin siguen dos pendientes diferentes:
1. La pendiente de la porcin de descarga (Um-Xo) es determinada como una funcin del
desplazamiento Um y la pendiente de la porcin inicial OY, de una manera similar al
modelo de Takeda.
La pendiente de la porcin de descarga (Um-Xo) es determinada como una funcin del
desplazamiento Um y la pendiente de la porcin inicial OY, de una manera similar al
modelo de Takeda.
2. La porcin de recarga tiene una pendiente determinada por las coordenadas de los
puntos Xo y Um, donde Um representa un punto sobre la curva primaria simtrico a
Um con respecto al origen.
2. La porcin de recarga tiene una pendiente determinada por las coordenadas de los
puntos Xo y Um, donde Um representa un punto sobre la curva primaria simtrico a
Um con respecto al origen.
Esta suposicin es aceptable porque: ayuda a simplificar el modelo y para
deformaciones de amplitudes bajas, provee algo de suavizacin comparable con el
efecto pinching (endurecimiento cerca al origen).
Esta suposicin es aceptable porque: ayuda a simplificar el modelo y para
deformaciones de amplitudes bajas, provee algo de suavizacin comparable con el
efecto pinching (endurecimiento cerca al origen).
21
Debe de enfatizarse que este modelo no provee alguna disipacin de energa, a
menos que el sistema fluya (esta deficiencia tambin existe en los modelos de Otani y el
bilineal simple). Consecuentemente si la carga inicia con deformaciones de pequeas
amplitudes por debajo del punto de fluencia, el modelo considera la seccin elstica.
Esto no es cierto, en vista que en realidad el comportamiento no lineal en una seccin
de concreto armado se inicia inmediatamente despus del agrietamiento de la seccin.

Figura 2.12 MODELO Q-HYST
Fuerza
0
U'm
R
Y '
X1
Xo
Deformacin
Carga reversa
Y
Um
Curva primaria
2.5 MODELOS DE ELEMENTOS INELASTICOS
2.5.1 MODELO DE DOBLE COMPONENTE
El modelo de doble componente o de Clough supone que cada elemento consiste de
dos componentes, uno elastoplstico y otro completamente elstico actuando en paralelo.
La suma de las dos da como resultado una relacin bilineal momento-curvatura para el
elemento (figura 2.13)
Este tipo de modelo es muy conveniente para su cmputo por dos razones principales
(29): primero, la rigidez de un elemento cambia sbitamente antes que progresivamente (y
es muy simple su clculo) y permanece constante hasta la formacin de una nueva rtula o
la descarga de una antigua. Segundo, es simple calcular las rotaciones de la rtula plstica
y de esa manera obtener una medida de la demanda de la ductilidad sobre el miembro.
22
En la seccin 3.5.1, se describe con amplio detalle el elemento viga-columna, que
pertenece a este tipo de modelo, el cual se utilizar en el presente estudio.
Figura 2.13 MODELO DE DOBLE COMPONENTE
componente elasto plstico
nudo
componente elstico
nudo
M
miembro completo
componente
elasto plstico
= CURVATURA
componente
elstico

2.5.2 MODELO DE UNA COMPONENTE
El modelo de una componente est basado en la suposicin que el resorte inelstico
est fijo a los extremos de todos los elementos con flexibilidad rotacional i, (figura 2.14).
Dadas estas flexibilidades en los extremos, los coeficientes de rigidez originales para los
elementos se pueden transformar con el objeto de incluir los efectos de los resortes
inelsticos.
El modelo permite seleccionar alguna relacin momento rotacin deseada para estos
resortes, como el rgido plstico, bilineal, trilineal o una variedad de modelos de
degradacin de rigidez y/o resistencia, o modelos curvilneos.
En la seccin 3.5.2, se describe en detalle el elemento viga con rigidez degradante, que
pertenece a este tipo de modelo.
23
2.5.3 MODELO PARABOLICO
Este modelo planteado por Takizawa supone una distribucin parablica de la
flexibilidad a lo largo del eje del elemento. Los puntos de inflexin se encuentran dentro de
la longitud del elemento. Dos distribuciones son posibles dependiendo de los signos
relativos (M+ traccin abajo) de los momentos extremos Ma y Mb (ver figura 2.15). Este
modelo puede ser aplicado especialmente para elementos tipo muros de corte.
Figura 2.14 MODELO DE UNA COMPONENTE
i
L
i
j
j=j Mj
i=i Mi
j

2.5.4 MODELO DEL ELEMENTO DISCRETO
Figura 2.15 MODELO DE FLEXIBILIDAD DISTRIBUIDA
B
M
M . M > 0
1/EI
A
A
1/EI
A
M
B
O B
1/EI
M
B
M . M < 0
A
1/EI
O
1/EI
M
A
B A
1/EI
B
Este modelo propuesto por Takayanagi y Schnobrich (ver figura 2.16), considera a un
elemento dividido en segmentos pequeos, cada segmento con una rigidez a la flexin
uniforme que vara con la historia de esfuerzos del segmento. La concentracin local de la
accin inelstica puede tomarse en cuenta considerando segmentos ms pequeos en las
24
zonas de mayor concentracin de deformacin inelstica. Este modelo se usa para analizar
muros de corte especialmente.

Figura 2.16 MODELO DEL ELEMENTO DISCRETO







CAPITULO III


MODELO ANALITICO


3.1 INTRODUCCION
El mtodo analtico ms popular para sistemas inelsticos es el de integracin directa
paso a paso, en el cual el dominio del tiempo se discretiza en muchos intervalos pequeos
t, y para cada lapso las ecuaciones del movimiento se resuelven con los desplazamientos
y velocidades del paso anterior como datos iniciales, las caractersticas de rigidez en el
principio de los periodos considerados se toman constantes en todo este paso; los clculos
se producen en forma de paso a paso.
Para modelar las estructuras aporticadas se usar el modelo de rtula puntual (modelo
discreto), en este modelo se puede tener la respuesta ssmica a nivel de cada elemento y
por ende a nivel de cada piso.
A continuacin, se describirn algunos procedimientos y suposiciones, para el
tratamiento de los parmetros involucrados en el anlisis dinmico inelstico de estructuras
26
aporticadas, as como una descripcin de los elementos estructurales que se utilizarn en la
presente investigacin.

3.2 SUPOSICIONES RESPECTO A LA ESTRUCTURA Y MOVIMIENTO EN LA BASE
En un anlisis dinmico inelstico se asumen diversas simplificaciones para evitar
complicadas y costosas soluciones, siempre que se asegure una relativa representacin
realista del problema. Las siguientes suposiciones se realizaron para el anlisis:
1. La estructura es idealizada como un ensamblaje plano de elementos discretos
2. El anlisis es mediante el mtodo de rigidez directa
3. Cada nudo posee hasta tres grados de libertad de desplazamiento, como en un anlisis
tpico de prtico plano.
4. Se asume que la masa de la estructura est concentrada en los nudos, as que la matriz
de masa es diagonal.
5. El movimiento ssmico est definido mediante la historia en el tiempo de la aceleracin
del terreno.
6. Se asume que todos los puntos de apoyo se mueven en fase.
7. Se aplicaron cargas estticas (cargas de gravedad) antes de las cargas dinmicas
8. Se utiliz los siguientes elementos estructurales: el elemento viga-columna y el
elemento viga con rigidez degradante, los cuales se explicarn con ms detalle en las
siguientes secciones.
9. La respuesta dinmica es determinada mediante la integracin paso a paso,
suponiendo una aceleracin constante promedio durante cualquier paso.
10. Se utiliza para cada paso la rigidez tangente de la estructura y se asume un
comportamiento lineal de la estructura durante el paso.

3.3 METODO DE RIGIDEZ DIRECTA
3.3.1 GRADOS DE LIBERTAD DE LA ESTRUCTURA
Con respecto al anlisis, una estructura se puede idealizar como un nmero finito de
nudos de cuerpo rgido, o uniones, conectadas por un nmero finito de elementos
deformables o miembros. Los nudos pueden tener dimensiones finitas, pero se consideran,
27
por lo general, slo como puntos. Los elementos pueden, en general ser unidimensionales,
bidimensionales o tridimensionales. Las cargas se pueden aplicar a los nudos o a los
elementos.
El Mtodo de Rigidez Directa es un caso especial del Mtodo de Desplazamientos. La
base de este mtodo es que el estado de deformacin de la estructura se puede definir
completamente en funcin de un nmero finito de parmetros de desplazamiento asociados
con los nudos. Estos parmetros de desplazamiento son los grados de libertad de la
estructura. Uno de los pasos ms importantes en la idealizacin de una estructura para el
anlisis, es la seleccin de grados de libertad adecuados de tal forma que se permita todos
los modos de deformacin importantes de la estructura. Al mismo tiempo, el nmero total
de grados de libertad debera mantenerse a un mnimo para reducir el tiempo de cmputo.
El programa de cmputo que se usar en este estudio se limita a las estructuras planas
y para el anlisis de cada nudo se permite tener hasta tres grados de libertad (traslaciones
paralelas a los ejes X e Y, y rotacin alrededor del eje Z). En muchas estructuras, no es
necesario que todos los nudos tengan tres desplazamientos, por tanto, en el programa se
han realizado previsiones para imponer dos tipos de restricciones: desplazamientos ceros o
idnticos.
La relacin entre los desplazamientos potenciales de los nudos (tres en cada nudo, en
este caso) y los grados de libertad de la estructura actual se definen a travs de la
construccin de un arreglo de identificacin (ID). En todos los casos, el tamao del arreglo
ser el (nmero de nudos) x (3). Para cada desplazamiento potencial, se introduce el
nmero del grado de libertad de la estructura correspondiente. Si no hay grado de libertad
correspondiente (es decir, el desplazamiento se restringe a cero), la entrada es cero. Si
dos o ms desplazamientos se restringen a ser idnticos, se les asigna el mismo nmero de
grado de libertad. El uso de un arreglo ID de este tipo, tuvo su origen en el programa SAP
(33), el procedimiento utilizado para generar el arreglo es similar al que utilizan en ste.

3.3.2 MATRIZ DE RIGIDEZ
a) Rigidez del Elemento
28
En el mtodo de rigidez directa, la matriz de rigidez de la estructura es ensamblada a
partir de las matrices de rigidez de cada elemento individualmente, mediante un proceso
simple de adicin directa. Cada elemento tiene un nmero posible de grados de libertad,
que consta de deformaciones nodales y de movimientos de cuerpo rgido. Estos grados de
libertad del elemento se deben relacionar con los desplazamientos de los nudos a los
cuales el elemento est conectado, por ltimo a los grados de libertad de la estructura.
La relacin entre los grados de libertad del elemento y de la estructura se define por una
matriz de localizacin denominado arreglo LM, en donde el trmino K-simo del arreglo
LM es el grado de libertad de la estructura que corresponde al grado de libertad K-simo
del elemento. Este arreglo se debe construir y almacenar para cada elemento de la
estructura, ya que se utiliza para ensamblar la matriz de rigidez de la estructura, el vector de
carga y en la conversin de los desplazamientos de la estructura a las deformaciones del
elemento.

b) Ensamblaje de la Rigidez de la Estructura
Una vez que la rigidez del elemento [k] y el arreglo LM {l}, se ha establecido para
cualquier elemento, resulta una tarea simple el ensamblaje de los coeficientes de rigidez del
elemento dentro de la matriz de rigidez de la estructura [K]. En general, el trmino (i, j) de
[k] ser adicionado a la ubicacin (li, lj) de [K].
Si uno de los trminos li o lj es cero, esto indica que no hay grado de libertad de la
estructura correspondiente (es decir, el desplazamiento nodal es cero). Por tanto, se
ignora, el coeficiente de rigidez del elemento correspondiente. Si dos o ms grados de
libertad del elemento se limitan a tener valores idnticos, sencillamente, los coeficientes de
rigidez correspondientes se adicionan simultneamente en la rigidez de la estructura
durante el ensamblaje.

3.3.3 VECTOR DE CARGA
Las cargas estticas sern aplicadas a la estructura antes de que comience el sismo.
Estas cargas se pueden aplicar directamente en los nudos o pueden originarse dentro de
los elementos.
29
Cuando las cargas se aplican directamente en los nudos, el vector de carga de la
estructura {R}, se ensambla mediante la utilizacin del arreglo ID. Si las cargas se originan
dentro de los elementos, se deben convertir stas, a cargas nodales equivalentes. Las
cargas nodales equivalentes para cualquier elemento, son iguales y opuestas a las fuerzas
de fijacin nodal, las que se deben aplicar al elemento para prevenir el desplazamiento de
sus nudos. Las cargas nodales equivalentes se adicionan en el vector de carga de la
estructura, casi de la misma manera como las cargas que se aplican directamente en los
nudos, excepto que se utiliza el arreglo LM para cada elemento, antes que el arreglo ID.
En uno u otro caso, si un desplazamiento nodal en particular no tiene grado de libertad de la
estructura correspondiente (es decir, la entrada en el arreglo LM e ID es cero), se ignorarn
cualquiera de las cargas correspondientes. Si dos o ms desplazamientos nodales son
asignados al mismo grado de libertad de la estructura, cualquier carga correspondiente ser
adicionada juntamente.
Durante el anlisis de la respuesta dinmica, se desarrollan cargas efectivas que
representan la inercia y los efectos de amortiguamiento. Estos vectores de carga, se
ensamblan mediante procedimientos similares a los utilizados para las cargas estticas.

3.3.4 ECUACIONES DE EQUILIBRIO Y SOLUCION
Las ecuaciones de equilibrio que deben ser resueltas en cualquier etapa del anlisis se
pueden escribir as:
[K] {r} = {R}
En donde [K] es la matriz de rigidez de la estructura; {R} es el vector de carga y {r} el
vector de valores desconocidos de los grados de libertad de la estructura.
En el programa de cmputo que se utiliza, estas ecuaciones se resuelven mediante un
algoritmo eficiente basado en la eliminacin de Gauss (17). Se almacena la matriz de
rigidez de la estructura en forma compacta para optimizar la utilizacin de la memoria, y
durante la operacin de eliminacin se evitan todas las operaciones aritmticas
innecesarias.
Cuando los desplazamientos nodales se han determinado, el clculo contina con las
deformaciones del elemento, y en consecuencia, se pueden encontrar las acciones del
30
elemento. Para los elementos inelsticos, es necesario tomar en cuenta la historia de la
deformacin del elemento. El procedimiento se trata con ms detalle en las siguientes
secciones.

3.4 ANALISIS PASO A PASO
3.4.1 EQUILIBRIO DINAMICO
En cualquier instante, la ecuacin de equilibrio dinmico se puede escribir como

[ ]{ } [ ]{ } [ ]{ } { } dP dr K r d C r d M
T T
= + +
& & &
(3.1)

En donde{ }, { }y{ son incrementos de aceleracin, velocidad y de
desplazamientos, respectivamente, en los nudos; {dP} es el incremento de la carga
aplicada; [M] es la matriz de masa; as como [K
T
] y [C
T
] son valores tangenciales de las
matrices de rigidez y de amortiguamiento para la estructura en su estado actual. Para un
intervalo de tiempo finito, t, la ecuacin (3.2) lo satisface aproximadamente.
r d
& & r d
&
} dr

[ ]{ } [ ]{ } [ ]{ } { } P r K r C r M
T T
= + +
& & &
(3.2)

En donde{ } r
& &
, { }, y r
&
{ } r { } P son incrementos finitos de aceleracin, velocidad,
desplazamiento y de carga, respectivamente, y las matrices de rigidez y de
amortiguamiento tangencial corresponden al estado de la estructura al inicio del intervalo de
tiempo.
Ya que pueden ocurrir cambios en el estado de la estructura durante el intervalo de
tiempo, el nuevo estado al final del intervalo, obtenido mediante la solucin de la ecuacin
3.2 no puede satisfacer exactamente el equilibrio. El procedimiento generalmente utilizado,
es resolver la ecuacin 3.2, determinar cualquier error que se pueda introducir, y
compensar estos errores mediante la aplicacin de cargas correctivas durante el siguiente
intervalo de tiempo. Los procedimientos utilizados se describen en las siguientes
secciones.
31

3.4.2 MATRIZ Y COEFICIENTES DE AMORTIGUAMIENTO
Se asume que el amortiguamiento viscoso, resulta de una combinacin de efectos
dependientes de la rigidez y de la masa, de modo que se cumple la ecuacin (3.3):


[ ]
(3.3)
[ ] [
C M K
T
= +
]
T
n
*
M
n
*

En donde los coeficientes y son constantes. El siguiente procedimiento (28) se
puede utilizar para seleccionar dichos coeficientes.
Si el sistema fuera desacoplado en modos normales, la relacin entre el
amortiguamiento, masa y rigidez generalizada para el modo n-simo sera:
C M (3.4) K
n n
* *
= +
Luego, de las relaciones (3.5) y (3.6):
(3.5)
*
n n n
*
n
M 2 C =

(3.6) K
n n
*
=
2
En donde
n
es la fraccin de amortiguamiento crtico en el modo n-simo, y
n
es la
frecuencia circular del modo n-simo, se sigue esto para cualquiera de los valores
especficos de y , el cual implica un amortiguamiento en el modo n-simo indicado en
(3.7) o (3.8):
2 2
n
n
n

+

= (3.7)
n
n
n
T 4
T
+

= (3.8)
En donde T
n
es el periodo del modo n-simo.
La Figura 3.1 muestra la influencia sobre de los dos trminos de la ecuacin 3.8. El
amortiguamiento dependiente de la rigidez, parece implicar amortiguamiento ms alto con el
perodo decreciente, y lo contrario con el amortiguamiento dependiente de la masa.
32
Si slo se asume el amortiguamiento dependiente de la masa, el amortiguamiento i en
un modo con periodo Ti se puede obtener mediante (3.9).
amortiguamiento M
AMORTIGUAMIENTO
amortiguamiento K
PERIODO
Figura 3.1 EFECTOS DEL AMORTIGUAMIENTO DEPENDIENTE
DE LA MASA Y DE LA RIGIDEZ

i
i
T
4
= (3.9)
Similarmente si slo se asume el amortiguamiento dependiente de la rigidez

=
i i
T
(3.10)
En ambos casos, el amortiguamiento para todos los valores de T se determinan una vez
que y han sido especificados.
Si tanto como se desean fijas, entonces se puede definir la cantidad del
amortiguamiento crtico en dos perodos diferentes. Es decir,
i
i
i
T 4
T
+

= (3.11)
j
j
j
T 4
T

+

= (3.12)
En consecuencia,

( )
2
i
2
j
i i j j
T T
T T 4

= (3.13a)
33

( )
) T T (
T T T T
2
i
2
j
j i i j j i


= (3.13b)
Para el anlisis prctico, se deberan seleccionar los valores correspondientes de T y
luego determinar los valores de y utilizando las ecuaciones 3.9, 3.10 o 3.13. El
amortiguamiento dependiente de la rigidez parece ser lo ms razonable, ya que implica
progresivamente el amortiguamiento ms alto en los modos ms altos.

3.4.3 INTEGRACION PASO A PASO
Una serie de mtodos de integracin, basados en variaciones supuestas de
aceleraciones nodales durante el intervalo de tiempo, ha sido extensamente utilizado (5,
19). El mtodo basado en una variacin lineal de aceleracin con respecto al tiempo, es
conocido por ser inestable en la presencia de modos de vibracin, con periodos que
exceden aproximadamente en un tercio al intervalo de tiempo (19), en consecuencia, es
inadecuado para su uso. El programa de computo que se utilizar en este estudio, ha
adoptado el mtodo basado en una aceleracin constante promedio dentro del intervalo de
tiempo (que es el caso particular del mtodo de Newmark, cuando =0.25). Este mtodo
tiene la ventaja de ser incondicionalmente estable para todos los perodos y pasos del
tiempo, as como, de no introducir amortiguamiento en el sistema.
En cualquier anlisis, se puede esperar mayor exactitud segn se reduzca el intervalo
de tiempo de integracin. Sin embargo, para reducir al mnimo los esfuerzos de cmputo,
es importante seleccionar en tanto que sea posible un intervalo de tiempo mayor. Sern
necesarios la experiencia y el experimento para determinar el intervalo de tiempo ms
apropiado. En la referencia (22), se establece el valor lmite del intervalo de integracin
para satisfacer los requisitos de estabilidad como:


4 1
T
t (3.14)

En donde T es el perodo fundamental de la estructura
Las ecuaciones bsicas para el mtodo de aceleracin constante promedio, se
muestran en el cuadro 3.1.
34
CUADRO 3.1
( )
1 O
r r
2
1
r
& & & & & &
+ =

METODO DE ACELERACION CONSTANTE PROMEDIO
a) Aceleracin, r1 como incgnita

r1
ro
r
r1
r1
t
t
ro
r
ro
r
(3.1.a)

(3.1.b)

(3.1.c)

(3.1.d)

(3.1.e)

b) Reordenando, en trminos de r
t
2
r r 2 r
O

+ =
& & (3.1.f)
2 O O
t
4
r
t
4
r r 2 r

=
& & & & & (3.1.g)
c) Ecuacin de equilibrio dinmico
[ ]{ } [ ]{ } [ ]{ } { } P r K r C r M = + +
& & &
(3.1.h)
Reemplazando:
[ ] [ ] [ ] { } { } [ ] [ ]{
O O O 2
r 2 C
t
4
r r 2 M P r K C
t
2
M
t
4
& & & &
+

+ + =

}
dt r r r
t
O
O

+ =
& & & &
(3.1.i)
2
t
r
2
t
r r r
1 O O 1

+ + =
& & & & & &
dt r r r
t
O
O

+ =
&
4
t
r
4
t
r t r r r
2
1
2
O O O 1

+

+ + =
& & & & &
35

Reemplazando la ecuacin 3.3, en la ecuacin 3.1.i proveniente del cuadro 3.1, se puede
escribir como:

{ } { } [ ] [ ]{ }
o T o o o T
2
r 2 K r 2 r
t
4
r 2 M P r K 1
t
2
M
t
2
t
4
& & & & &
+

+ + =

(3.15)
Para evitar la necesidad de evaluar el trmino [ ]{ }
O T
r 2 K
&
la siguiente transformacin
propuesta por Wilson (32) es introducida
{ } { } { } { } { }
o
r 2 r 1
t
2
r r x
& &

= + = (3.16)
Por lo tanto, la ecuacin 3.15 se puede escribir as
[ ]{ } { } [ ]

+ +

+ + = +
o o o o T
r 2 r 2 r
t
4
r 2 M P x K M
& & & & &
(3.17)
En donde:
1
t
2
t
2
t
4
2
+

=
Cuando se ha determinado {x} luego se calcula el incremento del desplazamiento
nodal
{ } {
o
r 2 x
1
t
2
1
r
&
+
+

= } (3.18)
Y, en consecuencia el incremento de la velocidad se calcula de la ecuacin 3.1.f del
cuadro 3.1.

3.4.4 CORRECCION DEL EQUILIBRIO DEBIDO AL AMORTIGUAMIENTO
DEPENDIENTE DE LA RIGIDEZ
A causa de los cambios en el estado de fluencia de los elementos estructurales, la
matriz de rigidez tangencial puede cambiar de un paso a otro. En consecuencia, el
amortiguamiento dependiente de la rigidez y la matriz de amortiguamiento, pueden cambiar,
36
repentinamente entre los pasos de tiempo. Esto introduce un desbalance en el equilibrio al
inicio del nuevo intervalo de tiempo, el cual puede llevar, a la acumulacin de errores
significativos si no se corrige.
Hagamos que el equilibrio, al final de un intervalo de tiempo (tiempo t
-
) se represente
por
{ } { } { } { } F F F
M C K
+ + = P (3.19)
En donde los trminos {F
M
}, {F
C
} y {F
K
} son las fuerzas totales en los nudos debido a las
fuerzas de inercia, amortiguamiento y de miembro respectivamente, y {P} es la carga
externa total. La fuerza de amortiguamiento en este tiempo se da por
{ } [ ]{ } F M K
C
= + & r
T
(3.20)
En donde { } r
&
es el vector de velocidad. Al inicio del prximo intervalo (tiempo t
+
) las
cantidades de {F
M
}, {F
K
}, {P} y {r} no sern cambiadas, pero {F
C
} puede cambiar a
{ } { } [ ]{ } F F M K K
C C T
+ = + + & r (3.21)
En donde [K] es el cambio en la rigidez tangencial y {F
C
} es el cambio
correspondiente en la fuerza de amortiguamiento.
Ya que el equilibrio se debe satisfacer, es necesario que una carga adicional {F
C
} sea
aplicada externamente, de modo que la nueva ecuacin de equilibrio en el tiempo t
+
es
{ } { } { } { } { } { } F F F F P F
M C C K
+ + + = +
C

r
(3.22)
Sin embargo, no se puede permitir que esta carga externa adicional contine actuando,
porque puede modificar significativamente el resultado de la respuesta de la estructura. Por
lo tanto, es necesario eliminar esta carga mediante la aplicacin de una carga opuesta e
igual de magnitud
{ } [ ]{ } = F K
C
& (3.23)
En el programa DRAIN-2D, esta carga es aplicada como una carga correctiva durante el
prximo paso de integracin. De esta manera se permite que la carga {F
C
} acte en un
pequeo tiempo, pero se cree que el error por esto es despreciable. Debe notarse que la
carga {F
C
} es en s misma aproximada ya que se basa en una forma supuesta para la matriz
de amortiguamiento.
37

3.4.5 RIGIDEZ TANGENCIAL
a) Rigidez Elstica
La rigidez elstica tangencial para cualquier elemento estructural se puede determinar
mediante procedimientos bien establecidos. Los detalles de los modelos inelsticos para
los elementos estructurales que se usarn en el presente estudio, se explicarn en las
siguientes secciones. La rigidez cambiar conforme fluye el elemento, ya sea sbitamente
o progresivamente, para los elementos que se emplearn los cambios son esencialmente
sbitos. En consecuencia, se modificar la matriz de rigidez de la estructura, al final de
cualquier intervalo de tiempo en que uno o ms elementos cambian su estado de fluencia.
La modificacin se lleva a cabo mediante la adicin de los cambios en la rigidez del
elemento, dentro de la rigidez actual de la estructura, y no reensamblando completamente
una nueva matriz de rigidez.

b) Rigidez Geomtrica
El efecto P es importante en los edificios altos, el cual incrementa significativamente
la fuerza cortante efectiva en cada piso. Esto es un efecto de segundo orden, asociado a
pequeos, pero no menos significativos, cambios en la forma de la estructura. Se ha
demostrado (8) que el efecto P puede ser considerado con bastante exactitud para un
edificio alto, mediante la adicin a la rigidez elstica de las columnas, de una rigidez
geomtrica basada sobre la fuerza axial en las columnas, bajo cargas de gravedad
nicamente. Esencialmente, la adicin de una rigidez geomtrica constante, permite
realizar una aproximacin lineal, de los efectos no lineales que resultan de los grandes
desplazamientos de una estructura.

3.4.6 DETERMINACION DEL ESTADO
a) Incremento de Fuerza en el Elemento
Dentro de cualquier intervalo de tiempo, se calculan los incrementos en los
desplazamientos de los nudos y, por tanto, los incrementos en las deformaciones del
miembro. Es necesario calcular los incrementos en las fuerzas del elemento que
38
corresponden a estos incrementos de deformacin, teniendo en cuenta la fluencia del
elemento (y si se desea, los efectos de los grandes desplazamientos). A todo este proceso
puede llamarse, la fase de los clculos de la determinacin del estado, ya que se
determina el nuevo estado de fuerza y de deformacin en cada elemento.
Para los elementos que se usarn en el presente estudio, se asume que la fluencia
ocurre sbitamente antes que progresivamente. Por tanto, la relacin entre la fuerza del
elemento y la deformacin se puede representar mediante una curva lineal por partes, como
se ilustra en la figura 3.2.
Cada cambio de pendiente corresponde a la fluencia en un nuevo punto o a la descarga
en un punto que fluy previamente, y se le denominar un evento. En general para
cualquier elemento, pueden ocurrir varios eventos en un solo intervalo de tiempo, tal como
se muestra en la figura 3.3.
Figura 3.2 TIPO DE COMPORTAMIENTO NO LINEAL
FUERZA
=EVENTO
DEFORMACION
39
La linealizacin para el intervalo de tiempo, est basada en el estado del elemento al
inicio del intervalo. Por lo tanto, el comportamiento calculado sera correcto, si no ocurrieran
eventos dentro del intervalo de tiempo y el comportamiento seguira la lnea discontinua en
la figura 3.3. El incremento de fuerza en el elemento, que corresponde a este
comportamiento, ser llamado el incremento de fuerza lineal, {S
L
}. En

efecto, el
comportamiento no ser lineal, por lo general, y seguir una lnea, tal como la lnea slida
en la figura 3.3. En este caso, los incrementos calculados en los desplazamientos de los
nudos y la deformacin del miembro no sern correctos, ya que ellos sern calculados
asumiendo un comportamiento lineal. Sin embargo, es razonable, asumir que las mismas
deformaciones seran obtenidas durante el intervalo de tiempo si fuera considerada la no-
linealidad, debido a que el intervalo de tiempo es corto, as mismo porque la respuesta de la
estructura ser afectada sustancialmente por los efectos del amortiguamiento y de la
inercia. En consecuencia, un incremento no lineal en la fuerza del elemento, {S
NL
},
puede ser calculado como se indica en la figura 3.3. Un procedimiento general para
calcular este incremento de fuerza no lineal es la siguiente:
Figura 3.3 INCREMENTO DE FUERZA LINEAL Y NO LINEAL

S
L

S
N
L
tiempo
t
FUERZA, S
DEFORMACION, v
tiempo
t+

S
u
1. Sea {V
t
} el incremento de deformacin del elemento calculado para el intervalo de
tiempo, sea {S} el incremento de la fuerza no lineal que tiene que ser calculado. El
incremento {S} ser calculado de la suma de uno o ms subincrementos,
1
{S
1
},
40

2
{S
2
}, etc. como se indica en la figura 3.4. Inicializar el factor de escala acumulado,
, a cero. Establecer un contador de ciclo, i, a 1.
2. Colocar {Vi} = (1-) {Vt}. Para el estado actual del elemento, determine el
incremento de la fuerza del elemento, {Si}, que corresponde a {Vi}.
3. Determine la proporcin, i, de este incremento de fuerza del elemento, el que se
requiere para producir el prximo evento. Por lo general, dos tipos de eventos sern
posibles, esto es descarga y fluencia. Una primera prueba se debera hacer para
descarga y si esto ocurre, el valor de ser cero. Si no ocurre la descarga, el valor de
para producir la nueva fluencia se puede calcular a partir de la siguiente ecuacin

i
y i
i
S S
S
=

En donde S
y
es la fuerza del elemento en la fluencia; S
i
es la fuerza del elemento al
inicio de este ciclo; y S
i
es el incremento de la fuerza del elemento calculado para este
ciclo. Se debe encontrar el valor ms pequeo de
i
, considerando todos los
componentes de fuerza del elemento. Si este valor es menor que 1, ocurre un evento
en este ciclo. Si
i
excede de 1, no ocurrir eventos adicionales en el intervalo de
tiempo, y el valor
i
= 1 es utilizado.
4. Adicionar i {Si} a {S}. Realizar cambios apropiados en los datos utilizados, para
almacenar la informacin sobre el estado de fluencia del elemento.
5. Si i es menor que 1, reducir {V} a (1-i){V}, incrementar el contador de ciclo, i, en
uno y repetir desde el paso 2. Si i es igual a 1, se completa el clculo del incremento
de la fuerza no lineal.
41

b) Incrementos de Deformacin del Elemento
Figura 3.4 CALCULO DEL INCREMENTO DE FUERZA
v

S
2

S
2

S
1

S
1
t
FUERZA, S

S
3
DEFORMACION, v
t+t
Para las estructuras inelsticas, las deformaciones del elemento sern usualmente de
mayor importancia que las fuerzas del elemento, ya que las deformaciones dan una
indicacin directa de la demanda de ductilidad sobre la estructura. Durante la
determinacin de los incrementos de la fuerza del elemento segn el procedimiento ya
descrito, se debera calcular, tambin, los incrementos de deformacin inelstica
apropiados en cada ciclo y acumulados, a fin de dar los incrementos para el intervalo de
tiempo. Las magnitudes de las deformaciones inelsticas, para los elementos
seleccionados que se utilizarn en la presente investigacin se describen en los
siguientes captulos.

c) Desbalance en el Equilibrio
Si no ocurren eventos durante un intervalo de tiempo y si se ignoraran las no
linealidades debido a los grandes desplazamientos, entonces la suposicin de un
comportamiento lineal durante el intervalo es correcta, y el estado al final del intervalo
satisface el equilibrio dinmico. Es decir, el equilibrio ser satisfecho para los
incrementos de la fuerza del elemento igual a los incrementos lineales {S
L
}, pero no
42
para los incrementos no lineales {S
NL
}. En el nivel del elemento, el desbalance del
equilibrio es:
{ } { } { } S S S
U L
=
NL
(3.24)
Mediante transformaciones apropiadas, y despus de sumar todos los elementos, este
desbalance se puede expresar en trminos de fuerzas nodales desbalanceadas, {F
U
}.
Para evitar errores que puedan surgir de desbalances en el equilibrio acumulados sobre
muchos intervalos de tiempo, una carga correctiva se aplica en el prximo intervalo de
tiempo para restablecer el equilibrio. Desde la ecuacin 3.19, se satisface el equilibrio
al final de cualquier intervalo de tiempo si:

{ } { } { } { } F F F
M C KL
+ + = P (3.25)
En donde {F
KL
} son las fuerzas en los nudos debido a las fuerzas del elemento {S} +
{S
L
}, asumiendo un comportamiento lineal. A causa de la fluencia, sin embargo, el
efecto de la fuerza del elemento ser realmente {F
KL
} - {F
U
}. En consecuencia, la
ecuacin del equilibrio debe ser
{ } { } { } { } { } { } F F F F P F
M C KL U U
+ + = (3.26)
En donde una carga externa ficticia, -{F
U
}, se debe aplicar. Para eliminar esta carga
temporal, una carga +{F
U
} se debe adicionar a las cargas para el prximo intervalo de
tiempo. As, se permite que la carga externa ficticia acte durante un periodo corto, su
efecto no ser sustancial para propsitos prcticos. Luego, estas fuerzas se incluyen
en la estructura del vector de carga para el prximo intervalo. Debe observarse, que
esta correccin del equilibrio se hace adems de cualquier otra correccin asociada con
el amortiguamiento dependiente de la rigidez.
La necesidad de corregir el equilibrio se podra evitar mediante la subdivisin de los
pasos de integracin cada vez que ocurra un evento. Este procedimiento ha sido
utilizado, por ejemplo, por Saiidi y Sozen (28). Sin embargo, el procedimiento es
complejo y dificultoso, y es tambin probable que el procedimiento dure ms tiempo; en
un estudio preliminar de Higginbotham (13) no parece ser ms exacto que los
procedimientos basados en intervalos de tiempo uniformes y correcciones del equilibrio.

43
3.4.7 MOVIMIENTO SISMICO
El procedimiento paso a paso es aplicable para cualquier tipo de excitacin dinmica.
Al considerar movimiento de terreno tipo sismo, se asume que todos los puntos de apoyo se
mueven idnticamente y en fase. Para movimientos ssmicos, la ecuacin 3.2 se convierte
en:

[ ]{ } [ ]{ } [ ]{ } { } M r C r K r
a T r T r
&& & + + = 0 (3.27)
En donde { es el incremento de la aceleracin absoluta; y , son
incrementos de velocidad y desplazamiento, respectivamente, relativos al terreno. Observe
que se asume que la amortiguacin viscosa depende de las velocidades relativas. La
ecuacin 3.27, tambin se puede escribir de la siguiente manera:
} }
a
r
& &
{ }
r
r
&
{
r
r

[ ]{ } [ ]{ } [ ]{ } [ ]{ }
M r C r K r M r
r T r T r
&& & && + + =
g
(3.28)
En donde{ }
r
r
& &
es el incremento de la aceleracin relativa al terreno y{ }
g
r
& &
es el
incremento de la aceleracin del terreno. En consecuencia, para movimientos ssmicos, el
incremento de la carga, {P}, en las ecuaciones de ste captulo se reemplaza por -
[M] { }
g
r
& &
. Los trminos que se presentan en{ }
g
r
& &
slo corresponden a las componentes
traslacionales de aceleracin. Un solo valor se presentar para todas las aceleraciones X, y
un solo valor diferente para todas las aceleraciones Y. Los resultados de los
desplazamientos, velocidades y las aceleraciones calculadas mediante la solucin de la
ecuacin 3.1.i (Cuadro 3.1) sern cantidades relativas al terreno.

3.5 ELEMENTOS ESTRUCTURALES
El grado de complejidad del comportamiento del elemento inelstico, depende del
objetivo para el cual se propone el elemento. A continuacin se detallan dos elementos
estructurales, una de doble componente y la otra de una componente, ambos pertenecen al
programa de cmputo DRAIN-2D (15, 24) el ELEMENTO VIGA-COLUMNA y el ELEMENTO
VIGA CON RIGIDEZ DEGRADANTE, la primera de ellas se utilizar en la presente
investigacin:

44
3.5.1 ELEMENTO VIGA-COLUMNA
a) Caractersticas Generales
Los elementos viga-columna se pueden orientar arbitrariamente en el plano X-Y, y
poseen rigidez axial y flexional, se pueden tomar en cuenta las deformaciones por cortante
y los efectos de conexiones con extremo excntrico (nudo rgido). La fluencia puede ocurrir
solamente en rtulas plsticas concentradas en los extremos del elemento. El
endurecimiento por deformacin se aproxima por la suposicin de que el elemento consta
de dos componentes una elastoplstico y otro elstico en paralelo, las rtulas en la
componente elastoplstico fluyen bajo momento constante, pero, el momento en el
componente elstico puede continuar incrementndose.
Con este tipo de idealizacin del endurecimiento por deformacin, si el momento flector
en el elemento es constante y si el elemento es de resistencia uniforme entonces la relacin
momento-rotacin para el elemento tendr la misma forma como su relacin momento-
curvatura ( Figura 3.5b). Esto se debe a que la curvatura y la rotacin en este caso son
directamente proporcionales. Sin embargo, si el momento flector o la resistencia vara
entonces las curvaturas y las rotaciones no son tan proporcionales, y las variaciones de
momento-curvatura y de momento-rotacin pueden ser totalmente diferentes (Figura 3.5c).
Con el procedimiento del componente paralelo, se est especificando en efecto, una
relacin de momento-rotacin. Debe tenerse cuidado cuando intente relacionar sta a la
relacin momento-curvatura.
Pueden especificarse diferentes momentos de fluencia en los dos extremos del
elemento tanto para la flexin negativa y positiva. La interaccin entre la fuerza axial y el
momento en la fluencia producida se pueden considerar aproximadamente.
Las cargas estticas aplicadas a lo largo de la longitud de un elemento se pueden
considerar mediante la especificacin de valores de fuerzas de fijacin.
Los efectos de desplazamientos largos se pueden aproximar en el anlisis dinmico
mediante la inclusin de una rigidez geomtrica simple, basado en las fuerzas axiales del
elemento bajo la carga esttica.

45

Figura 3.5 RELACIONES MOMENTO-CURVATURA Y
MOMENTO-ROTACION
MOMENTO, M
CURVATURA,
(a)
M
,
M -
M -
M
M

M
M

,
M -
M -
M
(b) (c)

b) Deformaciones del Elemento
Un ELEMENTO VIGA-COLUMNA tiene tres modos de deformacin, estos son:
extensin axial, rotacin por flexin en el extremo i, y en el extremo j. La matriz de
transformacin del desplazamiento que relaciona los incrementos de deformacin y los
desplazamientos (Figura 3.6) es:


dv
dv
dv
X
L
Y
L
X
L
Y
L
Y
L
X
L
Y
L
X
L
Y
L
X
L
Y
L
X
L
dr
dr
dr
dr
dr
dr
1
2
3
2 2 2 2
2 2 2 2
1
2
3
4
5
6
0 0
1
0 1

0 (3.29)
En forma abreviada:
{ } [ ] { } dv a dr = (3.30)

46
Figura 3.6 DEFORMACIONES Y DESPLAZAMIENTOS
L
X
Y
ds2, dv2
ds3, dv3
ds1, dv1
(a)
(b)
dr1
dr2
dr3
dr5
dr6
dr4

Se asumen que X, Y y L permanecen constantes. Una rtula plstica se forma cuando
el momento en la componente elastoplstico del elemento alcanza su momento de fluencia,
luego, una rtula se introduce dentro de sta componente, permaneciendo inalterado el
componente elstico. La medida de la deformacin plstica por flexin es la rotacin de la
rtula plstica, a veces se utilizan otras medidas de demanda de ductilidad, con frecuencia,
en la forma de relaciones de ductilidad. Sin embargo, en ste caso, se cree que la rotacin
de la rtula plstica es tanto ms racional, como ms consistente y no muy difcil
interpretarla en la prctica, que una relacin de ductilidad.
Para un incremento de la rotacin por flexin total dv2 y dv3, los incrementos
correspondientes de la rotacin de la rtula plstica dvp2 y dvp3, estn dados por:

(3.31)
dv
dv
A B
C D
dv
dv
p
p
2
3
2
3

En donde A, B, C y D se dan tal como en el cuadro 3.2. Una descarga ocurre en la
rtula cuando el incremento de la rotacin en la rtula es de signo contrario al momento
flector.




47
CUADRO 3.2
COEFICIENTES PARA ROTACIONES EN LA ROTULA PLASTICA
Condicin de fluencia A B C D
Extremos elsticos 0 0 0 0
Rtula plstica slo en el
Extremo i
1 kij/kii 0 0
Rtula plstica slo en el
Extremo j
0 0 Kij/kjj 1
Rtulas plsticas en ambos
Extremos i y j
1 0 0 1
Los coeficientes kii, kij, y kjj son definidos mediante la ecuacin 3.33

Se asume que no ocurren deformaciones axiales inelsticas en los ELEMENTOS VIGA-
COLUMNA a causa de la dificultad de considerar la interaccin entre las deformaciones por
flexin y las axiales despus de la fluencia. Por esto, se incluye solamente un
procedimiento aproximado para considerar los efectos de interaccin tal como se explica en
la siguiente seccin. Este procedimiento no es estrictamente consistente, pero, se cree que
es razonable para muchas aplicaciones prcticas.

c) Superficies de Interaccin
Se pueden especificar tres tipos de superficies de interaccin de fluencia, tales como:
1) Tipo viga (cdigo =1, Figura 3.7a). Este tipo de superficie debera ser especificado
cuando las fuerzas axiales son pequeas o son ignoradas. La fluencia es afectada slo
por el momento flector.
2) Tipo columna de acero (cdigo =2, Figura 3.7b). Este tipo de superficie es para usarla
con columnas de acero.
3) Tipo columna de concreto (cdigo =3, Figura 3.7c). Este tipo de superficie es para
usarla con columnas de concreto.

48

Figura 3.7 SUPERFICIES DE INTERACCION DE FLUENCIA
P
M
My+
My-
(a) CODIGO DE FORMA = 1
(b) CODIGO DE FORMA = 2
P
My-
My+
M
Pyc
Pyt
A
B
D C
PA = PC
Pyc Pyt
PC
PA
Similar para B y D
(c) CODIGO DE FORMA = 3
Pyt
Pyc
My-
B
My+
A
M
P
49
Se asume que la seccin transversal es elstica, para cualquier combinacin de fuerza
axial y momento flector dentro de una superficie de fluencia. Si la combinacin fuerza-
momento se sita sobre o fuera de la superficie, se introduce una rtula plstica. Slo
temporalmente se permiten combinaciones fuera de las superficies de fluencia, siendo
compensadas mediante la aplicacin de cargas correctivas en el siguiente paso de
integracin.
Este procedimiento no es estrictamente correcto, debido a que las deformaciones por
flexin y las axiales, interactan despus de la fluencia y, por lo tanto, es incorrecto asumir
que slo la rigidez por flexin cambia mientras que la rigidez axial permanece invariable.
Sin embargo, se cree que este procedimiento es razonable para anlisis prcticos de
edificios.
Si una combinacin fuerza-momento va del rango elstico a ms all de la superficie de
fluencia en cualquier paso de integracin, se realiza una correccin del equilibrio tal como
se muestra en la Figura 3.8a. Asimismo, debido a que se asume que la rigidez axial
permanece incambiable, la combinacin fuerza-momento en una rtula plstica
generalmente se alejar de la superficie de fluencia dentro de cualquier paso de integracin,
tal como se muestra en la Figura 3.8b, por lo tanto, se realiza una correccin del equilibrio,
tal como se muestra, en el siguiente paso de integracin.
En realidad, la fuerza axial de un elemento con una superficie de interaccin tipo
columna nunca exceder el valor de fluencia para un momento cero. Sin embargo, debido
al procedimiento computacional usado, se pueden calcular fuerzas axiales que exceden a la
fluencia. Para las fuerzas axiales que exceden a la fluencia, se asume que los momentos
de fluencia son ceros. Los resultados impresos del programa deberan ser examinados
cuidadosamente e interpretados con precaucin. Si las fuerzas axiales que se calculan
para una columna se aproximan o exceden a la fluencia, probablemente ste implica un
dao severo en la columna.
50

d) Rigidez del Elemento.
El elemento es considerado como la suma de una componente elstico y otra
componente inelstico. Las acciones y deformaciones del elemento se muestran en la
Figura 3.6a. La rigidez axial es constante y est dada por:
Figura 3.8 CORRECCION DEL EQUILIBRIO PARA SOBRECARGAS
EN LA SUPERFICIE DE FLUENCIA
(a)
P
M
desbalance del
equilibrio
t
t+t
desbalance del
equilibrio
t
P
(b)
M
t+t
dS
EA
L
dv
1
=
1
(3.32)
En la que E = Mdulo de elasticidad y A = Area de la seccin transversal. La rigidez
elstica por flexin est dada por:


dS
dS
EI
L
k k
k k
dv
dv
ii ij
ij jj
2
3
2
3

(3.33)
En la que I = Momento de inercia y son coeficientes que dependen de la
variacin de la seccin transversal. Para un elemento uniforme, con un momento de inercia
real,
jj ij ii
k , k , k
k k y k
ii jj ij
= = = 4 2
C
B
. Los coeficientes deben ser especificados por el usuario
del programa y puede, si lo desea, considerar los efectos de deformaciones por cortante y
conexiones de extremos no rgidas, as como, variaciones de la seccin transversal.
Despus de que se formen una o ms rtulas, los coeficientes para la componente
elastoplstico cambian a tal como sigue: k k k
ii ij jj
' , ' , '
(3.34) k k A k
ii ii ij
' ( ) = 1
(3.35) k k D k
ij ij ii
' ( ) = 1
51
(3.36) k k D k
jj jj ij
' ( ) = 1 B
a
En las que A, B, C y D se definen en el cuadro 3.2
La rigidez en trminos de desplazamientos nodales se obtiene como:

[ ]
(3.37)
[ ] [ ] [ ]
K a k
T
T
T
=

La rigidez geomtrica en las coordenadas del elemento dw1 y dw2 de la Figura 3.9 es:

[ ]
k
S
L
G
=

1 1
1 1
(3.38)
Luego en trminos de desplazamientos nodales:

[ ]
(3.39)
[ ] [ ][
K a k
G
T
G
=
1
]
a
1

4
3
2
1
2
1
dr
dr
dr
dr
L
X
L
Y
0 0
0 0
L
X
L
Y
dw
dw
(3.40)
{ } [ ]{ } dw a dr =
1
(3.41)

Esta rigidez geomtrica, no es exacta para un ELEMENTO VIGA-COLUMNA, pero es lo
suficientemente preciso para tomar en cuenta el efecto P en prticos de edificios.

e) Fuerzas de Fijacin en los Extremos e Iniciales
Figura 3.9 DESPLAZAMIENTOS PARA LA
RIGIDEZ GEOMETRICA
dw1
dw2
Las cargas estticas aplicadas a lo largo de las longitudes de los ELEMENTOS VIGA-
COLUMNA pueden tomarse en cuenta mediante la especificacin de fuerzas de fijacin, tal
52
como se muestra en la Figura 3.10. Estas fuerzas son aquellas que deben actuar sobre los
extremos del elemento a fin de prevenir cualquier desplazamiento en ellos.
Para aquellas estructuras en las que se ha realizado el anlisis esttico en forma
separada, pueden especificarse estos resultados como fuerzas iniciales del elemento. La
conveccin de signos para estas fuerzas se muestra en la Figura 3.10a. Estas fuerzas no
se convierten en cargas en los nudos de la estructura, sino que se emplean simplemente
para inicializar las acciones en los extremos del elemento. Cualquier fuerza en los
extremos debido a otras cargas se aade luego, a las fuerzas iniciales.
Se pueden especificar fuerzas iniciales del elemento, adems de las fuerzas nodales
estticas y de las fuerzas de fijacin en el elemento, en cuyo caso las fuerzas del elemento
debido a la carga esttica son adicionadas a las fuerzas iniciales. Si se emplea la rigidez
geomtrica, sta se basa en la fuerza axial inicial ms la fuerza axial debida a la carga
esttica, y se incluye slo para la carga dinmica, mas no para la carga esttica.
Las fuerzas iniciales y de fijacin se definen como patrones estndares y cada elemento
puede ser identificado con un patrn estndar de fuerzas de fijacin por carga muerta,
fuerzas de fijacin por carga viva y fuerza inicial. Adems, pueden especificarse factores
de multiplicacin para escalar los patrones estndares.
f) Reduccin de la Carga Viva
En los edificios, las reducciones de carga viva basadas sobre un rea tributaria pueden
tener efectos importantes y deberan tomarse en cuenta. Las fuerzas de fijacin
especificadas para cualquier elemento, despus escalado mediante los factores
Figura 3.10 FUERZAS DE FIJACION EN LOS
EXTREMOS E INICIALES
(b) CODIGO = 1
Vi
Mi
Pi
Pj
Vj
Mj
Pi
(a) CODIGO = 0
Vi
Mi
Vj
Mj
Pj
53
especificados para dicho elemento, consideran cualquier reduccin de la carga viva
permitida para ese elemento.
Las fuerzas de fijacin para cualquier elemento constituirn, cuando cambien de signo,
cargas estticas sobre los nudos a los cuales se conecta el elemento y el programa tomar
en cuenta estas cargas. Sin embargo, con frecuencia, el factor de reduccin de carga viva
permitida para una columna de un edificio exceder al de las vigas que soporta, debido a
que las columnas soportan cargas tributarias de los distintos pisos. De all que, si las
cortantes de fijacin de la carga viva completa para cada viga se aplican a los nudos de la
estructura, las cargas acumuladas en las columnas pueden ser innecesariamente grandes.
Esto podra ser compensado con la reduccin de las cortantes de fijacin para proporcionar
las cargas correctas en la columna, pero las fuerzas cortantes calculadas para las vigas
seran demasiado pequeas. Un enfoque preferible, es aprovechar los factores de
reduccin de las cargas vivas, los que pueden ser especificados con los patrones de
fuerzas de fijacin y, los que se emplean como sigue:
Las fuerzas completas de fijacin en los extremos especificados, son usadas para la
inicializacin de las fuerzas axiales y cortantes del elemento. Sin embargo, para el clculo
de las cargas estticas en los nudos conectados al elemento, las fuerzas cortantes y axiales
de fijacin (pero no los momentos) debido a la carga viva son primero multiplicadas por el
factor de reduccin especificado. Las fuerzas que producen las cargas axiales en las
columnas pueden, en consecuencia, ser reducidas para considerar la diferencia en las
reducciones de carga viva permisible entre las columnas y las vigas; sin embargo, las
fuerzas cortantes calculadas para las vigas sern, an, correctas. Para las cargas muertas,
se ignora el factor de reduccin.

g) Deformaciones por Cortante
Si se desea pueden especificarse reas efectivas de corte por flexin. Luego, el
programa modifica la rigidez por flexin para considerar las deformaciones por cortante
adicionales. Sin embargo, las fuerzas de fijacin en los extremos no cambian, de manera
que si las deformaciones por cortante pueden ser importantes, los patrones de fuerzas de
fijacin especificados deberan tomar en cuenta estas deformaciones.
54

h) Excentricidad en los Extremos (Nudo Rgido)
En las estructuras aporticadas y prticos acoplados a muros de corte, las rtulas
plsticas se formarn cerca a las caras de la conexin antes que en la lnea central terica
de la unin. Este efecto puede ser aproximado mediante conexiones infinitamente
resistentes y rgidas, uniendo entre los nudos (ubicados en las lneas centrales de las
uniones) y los extremos del elemento, tal como se muestra en la Figura 3.11. La
transformacin del desplazamiento que relacionan los desplazamientos nodales
{
, con
aqullos en los extremos del elemento es:
} dr
n

(3.42)
dr
dr
dr
dr
dr
dr
Yi
Xi
Yj
Xj
dr
dr
dr
dr
dr
dr
n
n
n
n
n
n
1
2
3
4
5
6
1
2
3
4
5
6
1 0 0 0 0
0 1 0 0 0
0 0 1 0 0 0
0 0 0 1 0
0 0 0 0 1
0 0 0 0 0 1

Esta transformacin ha sido incorporada en el clculo de las deformaciones y rigideces

del elemento. Las excentricidades en los extremos se especifican como patrones
estndares con cada elemento, asignndosele un patrn apropiado.
Donde se especifican excentricidades en los extremos, los coeficientes de rigidez deben
aplicarse a la parte del elemento entre los extremos del mismo, ignorando la regin de la
unin. En forma similar, las fuerzas de fijacin en los extremos son aqullas que se aplican
a los extremos del elemento, que corresponden a las cargas aplicadas dentro de la longitud
del elemento. Los efectos de excentricidad en los extremos se toman en cuenta en la
transferencia de las fuerzas de fijacin a los nudos (es decir, las cargas del momento se
aumentan mediante acoplamientos creados por las fuerzas axiales y cortantes de fijacin).
Sin embargo, se debera observar que cualquiera de los factores de reduccin de carga viva
especificados, se aplica a las fuerzas axiales y cortantes de fijacin antes de que ellos se
transfieran de los extremos del elemento a los nudos.
55

Figura 3.11 EXCENTRICIDADES EN LOS EXTREMOS
Y
i
Xi
UNION
RIGIDA
NUDO
i
-Xj
-
Y
j
UNION
RIGIDA
j
NUDO
i) Salida de Resultados
Los siguientes resultados se imprimen en el archivo de salida para la condicin de carga
esttica (todos los elementos) y en cada tiempo de salida (slo para elementos
especificados) si se solicitan resultados tiempo-historia:
1. Cdigo de fluencia en cada extremo del elemento: cero indica que el extremo del
elemento es elstico y 1, que se ha formado una rtula plstica.
2. Fuerza axial, momento flector y fuerza cortante actuante en cada extremo con la
convencin de signo mostrada en la Figura 3.10a.
3. Rotaciones actuales de la rtula plstica en cada extremo.
4. Rotaciones de la rtula plstica positiva y negativa acumuladas hasta el tiempo actual.
Estos valores se acumulan tal como se muestra en la Figura 3.12.
Los valores positivos y negativos mximos de la fuerza axial, momento flector, fuerza
cortante y rotacin de la rtula plstica con sus tiempos de ocurrencia son impresos en los
intervalos de tiempo solicitados para los envolventes. As mismo se imprimen las
rotaciones acumuladas negativas y positivas de la rtula plstica.
56

Figura 3.12 PROCEDIMIENTO PARA EL CALCULO DE
DEFORMACIONES PLASTICAS ACUMULADAS
FUERZA O
MOMENTO
EXTENSION O
ROTACION
MAXIMA
DEFORMACION
POSITIVA
DEFORMACION ACUMULADA NEGATIVA
= SUMA DE EXCURSIONES NEGATIVAS EN FLUENCIA
DEFORMACION ACUMULADA POSITIVA
= SUMA DE EXCURSIONES POSITIVAS EN FLUENCIA
NOTESE QUE LA DEFORMACION MAXIMA NEGATIVA
ES CERO, SIN EMBARGO LA DEFORMACION
ACUMULADA NEGATIVA NO ES CERO
3.5.2 ELEMENTO VIGA CON RIGIDEZ DEGRADANTE
a) Caractersticas Generales
Este elemento fue formulado, para modelar vigas de concreto armado que represente
caractersticamente, propiedades de rigidez de flexin degradante cuando se someta a
cargas cclicas. Como en el caso del ELEMENTO VIGA-COLUMNA, el ELEMENTO VIGA
CON RIGIDEZ DEGRADANTE se puede orientar arbitrariamente en el plano X-Y. El
elemento tiene rigidez axial y de flexin, as mismo, se pueden considerar elementos de
seccin transversal variable mediante la especificacin apropiada de los coeficientes de
rigidez por flexin. Se pueden considerar las deformaciones por cortantes de flexin y los
efectos de conexiones rgidas en los extremos.
La fluencia puede tener lugar solamente en rtulas plsticas, ubicadas en los extremos
del elemento. El endurecimiento por deformacin y la rigidez de flexin degradante, se
57
aproximan mediante la suposicin de que el elemento consta, de un elemento de viga
linealmente elstica con resortes de rotacin no-lineal en cada extremo, tal como se
muestra en la Figura 3.13. Todos los efectos de deformacin plstica incluyendo los
efectos de rigidez degradante se introducen por medio de las relaciones del momento-
rotacin para el resorte de la rtula.
Figura 3.13 IDEALIZACION DEL ELEMENTO
viga elstica
r
e
s
o
r
t
e

r
o
t
a
c
i
o
n
a
l

i
n
e
l

s
t
i
c
o

La relacin del momento-rotacin para cada rtula es una versin extendida del modelo
de Takeda, el cual tiene el comportamiento mostrado en la Figura 3.14. La relacin bsica
se encuentra en la forma de una curva bilineal, con una rigidez inicial y una rigidez de
endurecimiento por deformacin posteriormente, las que son caractersticas de condiciones
de cargas monotnicas. La rigidez degradante de las rtulas se introducen cuando se
aplica la carga inversa. Tambin, la Figura 3.14 indica las reglas que son seguidas para
oscilaciones de amplitudes cortas. Los nmeros en los tramos de la relacin, son cdigos
de fluencia los que son impresos en los resultados, por el programa de cmputo.
58

LOS NUMEROS INDICAN EL
CODIGO DE FLUENCIA
3
0
My+
11
2
7
4
My-
8
6
10
4
9
MOMENTO (M)
ROTACION (R)
5
3
1
Figura 3.14 RELACION MOMENTO-ROTACION DE LA ROTULA
PARA EL MODELO DE TAKEDA
Las modificaciones que se han realizado para el modelo de Takeda se muestran en la
Figura 3.15. Estas incluyen:
1. Una reduccin de la rigidez de descarga por una cantidad que depende de la mxima
rotacin previa de la rtula y
2. Una incorporacin de una rigidez de recarga variable, que es ms amplia que la del
modelo original de Takeda y tambin, depende de la historia de la rotacin anterior.
Para un mayor alcance, estas extensiones reflejan el comportamiento observado en la
prctica. Sin embargo, algunas de las caractersticas asumidas para el modelo extendido,
particularmente, aquellas asociadas con las oscilaciones de amplitudes cortas han sido
seleccionadas sobre la base del criterio tcnico y la necesidad de evitar patrones
inconsistentes o ilgicos del comportamiento de la rtula. Tales caractersticas asumidas
son necesarias ya que los datos de pruebas no parecen estar disponibles para el caso de
ciclos de carga de amplitud corta, que siguen a los ciclos de amplitud larga.
La rigidez de descarga, , depende de la mxima rotacin de la rtula y es
controlada por el parmetro de entrada
K
U
. Este parmetro controla la rigidez de descarga
59
mediante la posicin del punto de recuperacin, , tal como se muestra en la Figura
3.15 (a). No debe ser negativo y puede comnmente estar situado entre cero y 0.4. Sin
considerar el valor de, el ngulo de la descarga nunca se permite que sea menor que el
ngulo de la recarga, de lo contrario podra producirse un ciclo histertico con una rea
negativa.
R
rec
La rigidez de recarga, , tambin depende de la mxima rotacin de la rtula y es
controlada por el parmetro de entrada
K
L
tal como se muestra en la Figura 3.15(b). El
parmetro no debe ser negativo y puede comnmente estar entre cero y 0.6.
Las relaciones de corta amplitud que se han asumido son muy similares a las del
modelo de Takeda. Sin embargo, han sido necesario formular reglas adicionales para
calcular la rigidez de recarga. En la situacin descrita en la Figura 3.15(c), el punto de
inversin , se encuentra dentro del punto de recuperacin positivo , ya que la
incursin inelstica no alcanz el punto mximo anterior. En esta situacin, no es
probable que el ngulo de recarga negativa se dirigiera hacia el punto
R
rev
R
rec
+
B
+
A

, que sera el
punto de recarga definida en trminos del parmetro y el punto mximo anterior. Si
se asumiera que la recarga se dirija hacia el punto
B

, el ngulo de recarga tendera a ser

excesivamente alta y podra ser incluso negativa si el punto fuera a situarse a la
izquierda del punto
R
rev
A

. En consecuencia, para este tipo de situacin la recarga a un

punto X es asumido, donde X se encuentra entre los puntos A y B en una posicin que
depende de la ubicacin de entre los puntos de recuperacin limites y . La
relacin especfica que se asume para la recarga negativa es:
R
rev
R
rec

R
rec
+

AX
AB
R R
R R
rec rev
rec rec
N
=

+
+
(3.43)
Con una relacin similar para la recarga positiva. Cuando N=1, el punto X se
posicionar entre A y B en la misma proporcin que se sita entre y R .
Cuando N>1, el punto X se cargar ms pesadamente hacia el punto A, y cuando N<1,
ms pesadamente hacia el punto B. Se espera utilizar un valor N=1, pero, se incluye la
R
rev
R
rec
+
rec

60
posibilidad de especificar algn otro valor para una mayor flexibilidad. Si el factor es
igual a cero, la recarga siempre tendr lugar hacia el punto B, como en el modelo de
Takeda, sin considerar el valor de N.
Tal como se ilustra en la Figura 3.15 (d), para una serie de ciclos de corta amplitud, el
comportamiento asumido es idntico al del modelo de Takeda, excepto que la rigidez de
recarga desde el punto C a la lnea AB se basa en la posicin del punto de inversin ms
reciente.
Se pueden especificar momentos diferentes de fluencia en los dos extremos del
elemento tanto para la flexin negativa y positiva. No existe interaccin entre la fuerza
axial y el momento flector al producirse la fluencia.
Las cargas estticas aplicadas dentro de la longitud de cualquier elemento se pueden
considerar mediante la especificacin de valores de fuerzas de fijacin en los extremos. Se
pueden incorporar los resultados de anlisis de carga esttica por separado, mediante la
especificacin de valores de fuerzas iniciales.
Los efectos por desplazamientos largos, se pueden aproximar en el anlisis dinmico,
mediante la inclusin de la rigidez geomtrica simple, basadas en las fuerzas axiales del
elemento bajo la carga esttica.
61
Rp
(b) Rigidez de Recarga
(a) Rigidez de Descarga
Rrec
Rp
Rp*
M
Ku
Rrec (-)
Ke
Ko
M
Rp
Rrec (+)
A
Ku
B
R
Ku
Ko
Rp*
R
Figura 3.15 RELACION MOMENTO-ROTACION DE LA ROTULA
PARA EL MODELO DE TAKEDA EXTENDIDO


62
(c) Recarga Luego de una Pequea Excursin en la Fluencia
(d) Comportamiento de Pequeas Amplitudes
A+
Rrev
B -
X
A -
C-
B -
A -
Ke
X
Rrec (-)
M
M
A+
C+
R
Ku Ku
B+
Rrec (+)
R
B+
Figura 3.15 (continuacin)


63
b) Deformaciones del Elemento
A diferencia del elemento viga-columna, el elemento de rigidez degradante tiene
adicionalmente dos modos de deformacin, una deformacin de rotacin adicional
asociada con cada uno de los dos resortes de rotacin, tal como se muestra en la Figura
3.16. Estos grados de libertad adicionales se condensan por fuera en el nivel del elemento,
y no se presentan como grados de libertad de la estructura. Los modos de deformacin
resultantes del elemento son iguales a los del elemento viga-columna y para la
transformacin del desplazamiento, la Ec. 3.29 es aplicable. Dentro del elemento, las
deformaciones se dividen en rotaciones de la rtula y deformaciones del miembro flexional
elstico.

c) Propiedades del Elemento
dv1
dv2
dv3
dv4
dv5
Figura 3.16 DEFORMACIONES DEL ELEMENTO
En los anlisis de estructuras con prticos de concreto armado, las dificultades surgen
inevitablemente en la asignacin de la seccin transversal y en las propiedades del
material.
El comportamiento no lineal del material de concreto, la composicin no homognea
de la seccin transversal y la presencia de una seccin transversal variable debido al
agrietamiento y a los cambios del rea de acero, se combinan todos, para dificultar la
asignacin de valores de rigidez exactos. Sern necesarios considerables pruebas y
experiencias antes que se pueda especificar con seguridad las propiedades del elemento.
64
Se debe especificar para el elemento linealmente elstico una rigidez de flexin efectiva
(EI) que podra ser tpicamente el valor EI para la seccin transversal agrietada. Tambin
se debe especificar una rigidez axial (EA) y una rigidez de cortante efectiva (GA).
Estn disponibles dos opciones para la asignacin de propiedades de rigidez a los
resortes. La primera es una opcin por defecto que se encarga de la generacin
automtica de las propiedades del resorte de la siguiente manera. Si los valores para la
rigidez del resorte se ingresan en blanco o cero, el programa asigna automticamente una
rigidez elstica inicial grande al resorte ( 10 ), de modo que cada rtula es
esencialmente rgida hasta la fluencia. Luego, el programa asigna una relacin de
endurecimiento por deformacin al resorte, del mismo modo que para una relacin P
8
*EI

de una viga en voladizo tal como se muestra en la Figura 3.17, igual al valor que se ingres
de la relacin de endurecimiento por deformacin.
La segunda opcin es ms general, pero, no siempre estar disponible. Si estn
disponibles las relaciones entre la carga, el desplazamiento y la rotacin para una viga
voladiza, se pueden utilizar las siguientes ecuaciones para obtener el EI efectivo, la rigidez
inicial del resorte y la relacin de endurecimiento por deformacin del resorte (Figura 3.17)

EI
M L
L
Y
Y
Y
=

3 2
2

(3.44)
K
M L
M L
EI
SP
Y
Y
Y
=

2
3
(3.45)

( )
( )
p
P L M L
K
P L M L
EI
u Y
SP u Y
u Y
=

2
3
(3.46)
En donde rigidez inicial del resorte; p = relacin de endurecimiento por
deformacin del resorte; M
K
SP
=
Y
= momento de fluencia,
Y
= desplazamiento de la viga
voladiza en la fluencia; rotacin de la viga voladiza en la fluencia; P
U
, cualquier
carga y deflexin correspondiente a lo largo de la parte del endurecimiento por deformacin
Y
= =
U
65
de la relacin y L= longitud de la viga voladiza. Una viga en cantilever idealizada
con estas propiedades seguir una relacin
P
P tal como se muestra en la Figura
3.17(d), y una relacin tal como se muestra en la Figura 3.17(e). La relacin del
endurecimiento por deformacin de la relacin
P
P no ser necesariamente igual a la que
se observa en una prueba. Si se utiliza esta segunda opcin, la EI calculada y las
propiedades del resorte, se deberan examinar para asegurar que sean fsicamente
razonables.
M
M
(a) Viga

P
(b) Cantilever Equivalente
(c) Cantilever Idealizado
P

Pu
y u
Py = My
L
(d) Relacin P-
Py
(e) Relacin P-
y

P
Figura 3.17 REPRESENTACION DEL COMPORTAMIENTO EN CANTILEVER

Las propiedades del resorte basadas en una relacin P de una viga en cantilever,
reproducir tambin una relacin P especfica para una viga en doble curvatura inversa
tal como se indica en la Figura 3.17. Para otras variaciones de momento flector, el modelo
idealizado no tendr necesariamente la misma relacin de momento-rotacin que la viga
real y los resultados de cualquier anlisis se debern interpretar con precaucin. Si es
compleja la variacin del momento flector y si es particularmente importante que se
obtengan resultados exactos, se puede utilizar el siguiente procedimiento de idealizacin.
Este procedimiento utiliza la segunda opcin para especificar las propiedades del resorte de
66
la rtula, su desventaja es que este procedimiento necesita que cada viga de la estructura
real sea subdividida en varios elementos de viga para el anlisis, esto ocasionar un
incremento en el tiempo del cmputo. Sin embargo, el procedimiento tiene la ventaja de
que se pueda considerar directamente la relacin de momento-curvatura para la viga, en
vez de alguna relacin de momento-rotacin menos directa.
Tal como se indica en la Figura 3.5, las relaciones del momento-curvatura y momento-
rotacin para una viga, sern proporcionadas de la misma forma a la viga que se somete a
un momento flector constante. Si cualquier viga completa se subdivide en elementos
pequeos, tal como se indica en la Figura 3.18(a), de modo que sea pequea la variacin
en el momento flector en la longitud de cualquier elemento en fluencia, entonces, cada
elemento en fluencia se somete aproximadamente a un momento constante, y resulta
razonable asumir formas similares para las relaciones de momento-rotacin y de momento-
curvatura. La relacin momento-rotacin que corresponde a una curvatura de momento
bilineal es mostrada en la Figura 3.18(b). Las propiedades del elemento viga
correspondiente se pueden obtener de las siguientes ecuaciones:
EI
M
Y
Y
=

(3.47)
(3.48) K
SP
= 10
8
EI

( )
p
M M
K M
u Y
SP u u
=

2
( ) EI
(3.49)
En donde M
U
, son cualquiera de los puntos correspondientes sobre la parte del
endurecimiento por deformacin de la relacin
U

M .
67

Mi
Mj
(a)
M

Mu
y u
My
(b)
elementos
cortos
elemento
largo
elastico
elementos
cortos
My
= (CURVATURA) (L)
L
M M
My
Figura 3.18 SUBDIVISION DE LA VIGA EN ELEMENTOS CORTOS
d) Rigidez del Elemento
El elemento se considera como un elemento viga linealmente elstica con resortes de
rotacin no lineal en cada extremo. Las deformaciones y acciones del elemento se
muestran en la Figura 3.16. La rigidez axial es constante y est dada por:
dS
EA
L
dv
1
=
1
(3.50)
Donde EA es la rigidez axial efectiva de la seccin transversal. La rigidez elstica de
flexin del elemento viga elstica est dada por:


dS
dS
EI
L
k k
k k
dv
dv
ii ij
ij jj
4
5
4
5

(3.51)
En donde EI es la rigidez efectiva de flexin de la seccin transversal, y Kii, Kjj, Kij son
coeficientes que se ingresan y dependen de la variacin de la seccin transversal. Para un
elemento uniforme, Kii= Kjj= 4.0 y Kij=2.0. Estos coeficientes deben ser especificados por
el usuario del programa y pueden modificarse para considerar efectos de deformaciones
por cortantes y conexiones con extremos no rgidas, as como, las variaciones de la
seccin transversal.
68
Para formar la rigidez del elemento completo, la rigidez del elemento lineal se invierte a
su forma de flexibilidad y se modifica a fin de considerar las deformaciones de cortantes
por flexin. Luego, los valores de flexibilidad de la rtula actual se adicionan para dar la
matriz de flexibilidad del elemento de la siguiente manera:
F
f
K
f
f f
K
e
ii
SPi
ij
ij jj
SPj
=
+
+

1
1
(3.52)
Luego sta matriz de flexibilidad se invierte para formar la matriz de rigidez [Ke]
deseada.
Del mismo modo que para el ELEMENTO VIGA-COLUMNA la rigidez en trminos de
desplazamientos nodales global se obtienen de:

[ ] [ ] [ ] [ ]
K a K
T
T
e
= a (3.53)
La rigidez geomtrica que se utiliza es exactamente la misma que para el ELEMENTO
VIGA-COLUMNA.

e) Salida de Resultados
Los siguientes resultados se imprimen para la condicin de carga esttica (todos los
elementos) y en cada tiempo de salida (solamente los elementos especificados) si se
solicita una historia en el tiempo:
1. Cdigo de fluencia en cada extremo del elemento, tal como se muestra en la Figura
3.14.
2. El momento flector, la fuerza cortante y la fuerza axial que acta en cada extremo con la
convencin de signos mostrada en la Figura 3.10(a).
3. Las rotaciones acumuladas primaria y secundaria, como positiva y negativa de la rtula
plstica en cada tiempo. Estos valores se acumulan tal como se muestra en la Figura
3.19. Las rotaciones primarias de la rtula plstica se acumulan durante las excursiones
de la carga en la zona del endurecimiento por deformacin de la curva bilineal bsica.
Las rotaciones secundarias de la rtula plstica se definen como deformaciones
incurridas durante la carga sobre las curvas de carga degradadas 4, 6, 8, 10.
69
Los valores mximos negativos y positivos de la fuerza axial, del momento flector, de la
fuerza cortante y de la rotacin de rtula plstica con sus tiempos de ocurrencia, se
imprimen en los intervalos de tiempo solicitados para las envolventes. Tambin, se
imprimen las rotaciones acumuladas negativas y positivas de la rtula plstica.

Figura 3.19 DEFINICION DE ROTACIONES ACUMULADAS
b1+
b2+
a1+
rotaciones acumuladas primarias = a +, a -
rotaciones acumuladas secundarias = b +, b -
a1-
M
a2+
R









CAPITULO IV


ESTUDIO ANALITICO USANDO DRAIN-2D


4.1 INTRODUCCION
En este captulo se describen las estructuras consideradas en el presente estudio, as
mismo se explican las consideraciones de los diferentes parmetros involucrados en el
anlisis dinmico inelstico que se efectuar a dichas estructuras. El anlisis se basa en el
procedimiento descrito en el Captulo III y para ello se utilizar el programa de cmputo
DRAIN-2D (15, 24), algunas informaciones tcnicas acerca de este programa de cmputo
se presentan en el Anexo A.

4.2 ESTRUCTURAS EN ESTUDIO
Con la finalidad de llevar a cabo la investigacin tratando de reproducir las condiciones
ms realistas se disearon prticos cuyas caractersticas se estimaron como las ms
usuales en nuestro medio. De esta forma se han considerado 19 prticos como modelos en
un rango de 2 hasta 20 pisos con cuatro crujas, cubriendo la gama ms comn de nmero
71
de pisos en edificios de concreto armado en nuestro medio, y escogiendo las luces de los
vanos y las alturas tratando de reflejar las dimensiones ms usuales. La altura del primer
nivel es mayor que la del piso tpico y las dimensiones de los elementos obedecieron a un
anlisis y diseo preliminar considerando sobrecargas de oficinas y losas aligeradas como
entrepisos. En las figuras 4.1 se presentan los esquemas de las estructuras analizadas.
Los prticos del P-2 al P-4 podran considerarse bajos, los prticos del P-5 al P-10 se
consideraran como medianos bajos, los prticos del P-11 al P-15 como medianos altos, y
los prticos del P-16 al P-20 como relativamente altos; los prticos en su conjunto dan un
rango de perodos de los edificios de Lima. Los prticos slo constan de columnas y vigas
y consta de una losa aligerada en una direccin con un espesor de 25 cm.
Se han considerado las siguientes caractersticas para el concreto:
fc = 210 Kg/cm2
Ec = 220,000 Kg/cm2
Gc = 91,667 Kg/cm2
= 0.20
72
Figura 4.1(a) GEOMETRIA DE LOS PORTICOS CONSIDERADOS
Prtico de 4 pisos (P-4)
Prtico de 2 pisos (P-2)
5.0
5.0
5.0 5.0
5.0 5.0
2.6 2.6
Prtico de 5 pisos (P-5)
5.0
3.0
5.0 5.0 5.0
3.0
5.0
5.0
Prtico de 3 pisos (P-3)
2.6
2.6
3.0
5.0
2.6
5.0 5.0
2.6
2.6
2.6
3.0
5.0
2.6
2.6
5.0
PLANTA
5.0 5.0 5.0
5.0
5.0
5.0

73
Figura 4.1(b) GEOMETRIA DE LOS PORTICOS CONSIDERADOS
PLANTA
5.0
2.7
5.0
Prtico de 8 pisos (P-8)
2.7
2.7
2.7
2.7
3.5
5.0 5.0 5.0
Prtico de 6 pisos (P-6)
5.0
3.0
2.7
2.7
5.0 5.0 5.0
2.6
2.6
2.6
2.6
2.6
5.0
5.0
5.0
5.0 5.0
2.7
5.0
Prtico de 9 pisos (P-9)
5.0 5.0
3.5
5.0
2.7
2.7
2.7
2.7
3.5
5.0
Prtico de 7 pisos (P-7)
5.0 5.0
2.7
2.7
2.7
5.0
2.7
2.7
2.7
2.7
2.7
2.7
5.0

74
Figura 4.1(c) GEOMETRIA DE LOS PORTICOS CONSIDERADOS
6.0
Prtico de 11 pisos (P-11)
6.0
Prtico de 13 pisos (P-13)
6.0
2.7
2.7
6.0
Prtico de 12 pisos (P-12)
6.0 6.0
2.7
2.7
2.7
2.7
4.0
6.0 6.0
6.0
Prtico de 10 pisos (P-10)
6.0 6.0
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
4.0
6.0 6.0
6.0 6.0 6.0
2.7
2.7
2.7
2.7
PLANTA
5.0
5.0
5.0
2.7
2.7
6.0 6.0
2.7
4.0
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
6.0 6.0
2.7
4.0
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7

75
Figura 4.1(d) GEOMETRIA DE LOS PORTICOS CONSIDERADOS
6.0
Prtico de 16 pisos (P-16)
6.0
Prtico de 14 pisos (P-14)
6.0 6.0
6.0 6.0
2.7 2.7
6.0
Prtico de 17 pisos (P-17)
2.7
4.0
6.0
2.7
2.7
2.7
6.0
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
4.0
2.7
6.0 6.0
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7 2.7
6.0
Prtico de 15 pisos (P-15)
2.7
2.7
2.7
4.0
6.0 6.0
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
4.0
6.0 6.0
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
5.0
PLANTA
6.0
2.7
2.7
2.7
6.0 6.0 6.0
5.0
5.0
2.7
2.7
2.7
2.7

76
Figura 4.1(e) GEOMETRIA DE LOS PORTICOS CONSIDERADOS
7.0
PLANTA
2.7
7.0 7.0
Prtico de 18 pisos (P-18)
7.0
4.0
7.0
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
6.0
6.0
6.0
7.0
2.7
2.7
2.7
2.7
4.0
7.0
Prtico de 19 pisos (P-19)
7.0 7.0 7.0
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
7.0 7.0

77
Figura 4.1(f) GEOMETRIA DE LOS PORTICOS CONSIDERADOS
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
2.7
4.0
2.7
2.7
2.7
2.7
7.0
Prtico de 20 pisos (P-20)
7.0 7.0 7.0
6.0
6.0
6.0
7.0
PLANTA
7.0 7.0 7.0

78
4.3 PREDIMENSIONAMIENTO DE LOS PORTICOS EN ESTUDIO
Para el predimensionamiento de las estructuras se ha considerado que los prticos
corresponden al prtico central de las edificaciones cuyas plantas tpicas se muestran en
las figuras 4.1; el ancho tributario para los prticos del P-2 al P-17 es de 5 m y 6 m para los
prticos del P-18 al P-20. As mismo, se indica que se han considerado tabiques en todos
los ejes en ambas direcciones, con una altura de 2.40 m y 0.15 m de espesor.
Para el predimensionamiento de los elementos estructurales se han seguido las
recomendaciones del ACI 318-99 (3) para zonas de alto riesgo ssmico.
Se han considerado las siguientes cargas de acuerdo a la Norma de Cargas E-020:
CARGAS MUERTAS
Peso techo aligerado = 350 kg/m2
Piso terminado = 100 kg/m2
Tabiquera mvil = 150 kg/m2
CARGAS VIVAS
Pisos tpicos = 250 kg/m2 (oficinas)
Azotea = 100 kg/m2
PESOS
Tabiques = 1350 kg/m3
Concreto armado = 2400 kg/m3

En el cuadro 4.1, se presenta el predimensionamiento de los elementos para los
distintos prticos.

79
CUADRO 4.1
PREDIMENSIONAMIENTO DE LAS ESTRUCTURAS
PORTICO NIVEL VIGAS
Cm2
COLUM
Cm2
PORTICO NIVEL VIGAS
Cm2
COLUM
Cm2
P 2 1 2 25 x 40 40 x 40 1 4 30 x 70 70 x 80
5 8 30 x 60 65 x 70
P 3 1 3 30 x 45 40 x 40
P-13
9 13 30 x 60 60 x 60
1 1 30 x 50 40 x 50 1 5 30 x 70 75 x 85
P 4
2 4 30 x 45 40 x 40 6 10 30 x 60 70 x 75
1 2 30 x 50 40 x 60
P-14
1114 30 x 60 65 x 70
P 5
3 5 30 x 45 40 x 50 1 5 30 x 70 80 x 90
1 2 30 x 50 45 x 60 6 10 30 x 65 75 x 80
P 6
3 6 30 x 45 40 x 55
P-15
1115 30 x 60 70 x 75
1 3 30 x 55 50 x 60 1 5 30 x 70 85 x 95
P 7
4 - 7 30 x 50 45 x 55 6 10 30 x 65 80 x 85
1 4 30 x 60 55 x 60
P-16
1116 30 x 60 70 x 80
P 8
5 - 8 30 x 55 45 x 55 1 6 30 x 70 90 x 100
1 4 30 x 60 60 x 65 7 12 30 x 65 85 x 90
P 9
5 - 9 30 x 55 55 x 60
P-17
1317 30 x 60 75 x 85
1 4 30 x 65 60 x 70 1 6 30 x 75 110x125
5 7 30 x 60 55 x 65 7 12 30 x 70 90 x 100
P 10
8 10 30 x 60 50 x 60
P-18
1318 30 x 65 80 x 90
1 4 30 x 65 65 x 70 1 7 30 x 65 115x130
5 8 30 x 60 60 x 65 8 14 30 x 60 95 x 105
P 11
9 - 11 30 x 60 55 x 60
P-19
1519 30 x 60 85 x 95
1 4 30 x 65 70 x 70 1 7 30 x 65 120x135
5 8 30 x 60 65 x 65 8 14 30 x 60 100x110
P -1 2
9 - 12 30 x 60 60 x 60
P-20
1520 30 x 60 90 x100


4.4 ANALISIS ELASTICO DE LOS PORTICOS EN ESTUDIO
Se efectu el anlisis esttico para determinar los efectos, debido a las cargas de
gravedad, tanto para la carga viva como para carga muerta.
Para el anlisis ssmico, se realiz un anlisis esttico mediante fuerzas laterales
estticas equivalentes, recomendadas por la NTE E-030 de 1,997 (26). Se consideraron los
siguientes parmetros ssmicos:
Factor de zona : Z = 0.40 (zona 3)
80
Factor de uso : U = 1.00 (edificacin tipo C)
Factor del suelo : S = 1.00 (S1 suelo rgido)
Perodo de corte del espectro
para el suelo tipo S1 : Tp = 0.40 seg.
Coeficiente de reduccin : R = 10.00 (para prticos de concreto armado)

El anlisis estructural se ha realizado mediante el programa de cmputo denominado
SAP2000 (Structural Analysis Programs) (31), basado en el mtodo de las rigideces por
procedimientos matriciales. Para el anlisis estructural se ha modelado la losa de cada piso
como diafragma rgido, se consideran zonas rgidas para las vigas, equivalente al 100% de
la distancia del eje a la cara de la columna para cada extremo de la viga.
Se han calculado los perodos de vibracin de los prticos mediante un anlisis
dinmico, para lo cual se han considerado que las masas estn concentradas a nivel de
cada piso (hiptesis de losa rgida), dichas masas corresponden a la carga muerta ms el
25% de la carga viva. En el Cuadro 4.2 se presentan los perodos de vibracin para todos
los prticos, y la Razn de la Masa Participante en forma porcentual, para los tres primeros
modos de vibracin de cada prtico. Este ltimo valor nos provee la informacin de cuan
importante es un determinado modo, en la respuesta estructural; como se puede apreciar
conforme el perodo aumenta, el segundo y tercer modo de vibracin son cada vez ms
importantes en la respuesta estructural, esto influir en el Factor de Distribucin de
Momentos en los nudos de la estructura originando distribuciones inesperadas de
momentos flexionantes en las columnas, el cual se estudiar detalladamente en el Captulo
V.
81
CUADRO 4.2
PERIODOS DE VIBRACION Y RAZON DE LA MASA PARTICIPANTE
PERIODO (seg.) RAZON DE LA MASA PARTICIPANTE (%) PORTIC
O
T1 T2 T3 MODO 1 MODO 2 MODO 3
P-2 0.37 0.10 92.44 7.56
P-3 0.49 0.15 0.08 91.23 7.56 1.21
P-4 0.57 0.18 0.10 84.18 11.29 3.46
P-5 0.61 0.20 0.10 82.08 12.14 3.62
P-6 0.74 0.24 0.13 81.03 12.12 3.95
P-7 0.83 0.28 0.15 83.89 10.86 3.01
P-8 0.86 0.29 0.16 84.08 10.47 3.03
P-9 0.93 0.31 0.17 82.48 11.08 3.29
P-10 1.12 0.38 0.21 83.29 10.97 2.98
P-11 1.24 0.42 0.23 82.40 11.32 3.14
P-12 1.35 0.46 0.25 81.49 11.65 3.35
P-13 1.37 0.47 0.27 77.97 13.50 4.11
P-14 1.41 0.49 0.27 77.42 13.75 3.89
P-15 1.44 0.49 0.28 78.21 12.41 4.14
P-16 1.53 0.52 0.29 77.25 12.61 4.41
P-17 1.57 0.54 0.30 77.03 12.62 4.21
P-18 1.82 0.62 0.34 74.35 12.82 4.79
P-19 1.86 0.63 0.35 74.20 12.89 4.64
P-20 1.95 0.66 0.36 73.32 13.17 4.85


4.5 DISEO DE LOS PORTICOS EN ESTUDIO
Con los efectos ssmicos obtenidos del anlisis por fuerzas laterales equivalentes y los
efectos de carga de gravedad, se determinan las envolventes para los elementos del prtico
de acuerdo a las combinaciones propuestas por la Norma Peruana de Concreto Armado
(27). Para este estudio, empleando el mtodo de la rotura se determin las reas de acero
en los extremos de vigas y columnas que son consideradas las zonas crticas, es decir,
donde puede desarrollarse las rtulas plsticas, hay que agregar que se trat de uniformizar
el refuerzo en cada nivel.
Con el fin de conseguir un comportamiento dctil de los elementos y por ende en la
estructura, se disearon los prticos teniendo en cuenta disposiciones especiales para el
diseo ssmico, en concordancia con el ACI 318-99 (3) con los requerimientos para un
82
prtico especial resistente a momento. El objetivo principal de estos requerimientos es dar
capacidad de disipacin de energa en el rango inelstico de respuesta.
En la figura 4.2 se presentan como ejemplo las reas de acero obtenidas en los extremos
de vigas y columnas (zonas crticas) para el prtico P-3. Las reas de acero indicadas,
corresponden al refuerzo longitudinal por flexin y son suficientes para definir las
propiedades plsticas de la seccin, sin embargo, para el diseo final es necesario definir
otras caractersticas tal como refuerzo por corte (refuerzo transversal), refuerzo a lo largo
del elemento, unin viga columna, etc., de manera que se garantice un comportamiento
dctil en los elementos del prtico.

4.6 PROPIEDADES DE FLEXION DE LOS ELEMENTOS
Es posible deducir curvas tericas momento-curvatura para secciones de concreto armado
con flexin y carga axial (21), aproximndose mediante una curva trilineal (ver figura 4.3a),
en que se definen, una primera etapa de agrietamiento (se alcanza el esfuerzo mximo de
traccin en el concreto), la segunda etapa de fluencia (el acero de refuerzo alcanza su
esfuerzo de fluencia en traccin) y la tercera al lmite de la deformacin til del concreto.
Sin embargo para llevar a cabo el anlisis inelstico con el modelo de rtula puntual, las
hiptesis de viga inelstica consideran el diagrama momento-rotacin de la seccin crtica,
por lo que se deben plantear ciertas consideraciones para la definicin de las relaciones
momento-rotacin en las secciones crticas. Estas relaciones dependen de las cargas de
gravedad (15) como se mostr en la figura 3.5.
Cuando se representa el comportamiento de la seccin de concreto armado con un
diagrama bilineal (ver figura 4.3b), la relacin momento-rotacin es bastante similar a la
relacin momento-curvatura, donde el momento ltimo es una buena aproximacin del
momento de fluencia; adems una vez que se desarrollan las grietas, como sucede en la
mayora de las vigas bajo cargas de servicio, la relacin M- es casi lineal desde la carga
cero hasta el inicio de la fluencia, en consecuencia la curva bilineal de la figura 4.3b es una
buena aproximacin para vigas inicialmente agrietadas (21).
83
5
.
0
0
F
i
g
u
r
a

4
.
2



D
I
S
T
R
I
B
U
C
I
O
N

D
E
L

R
E
F
U
E
R
Z
O

L
O
N
G
I
T
U
D
I
N
A
L

P
A
R
A

E
L

P
O
R
T
I
C
O

P
-
3
2
#
6
+
1
#
4
V I G A S
4
#
6
2
#
5
5
.
0
0
0
.
3
0
2
#
6
4
#
6
+
1
#
5
2
#
6
5
.
0
0
2.60 2.60
5
.
0
0
3.00
0.45
0.40
C O L U M N A S
6
#
6
0
.
4
0
84



Se ha procedido a calcular las relaciones momento-curvatura para todas las secciones
crticas (en los extremos) de todos los elementos de los prticos antes descritos, basados
en las propiedades de los materiales que se indican en el cuadro 4.3 y algunas relaciones
idealizadas para el concreto y el acero que se describen en la siguiente seccin.
My
Mu
M
u u y y
(a)

(b)

Primer
agrietamiento
Primera
fluencia
My
Mu
M
Figura 4.3 RELACIONES IDEALIZADAS MOMENTO-CURVATURA
PARA UNA SECCION DE CONCRETO ARMADO

CUADRO 4.3
PROPIEDADES DE LOS MATERIALES ASUMIDOS PARA LOS PORTICOS
CONCRETO
fc resistencia a compresin 210.00 kg/cm2

O
deformacin en el esfuerzo mximo alcanzado 0.0020

U
deformacin ltima 0.0040
Ec mdulo de elasticidad 220,000 kg/cm2
ACERO
fy esfuerzo de fluencia 4,200.00 kg/cm2

sy
deformacin en el punto de fluencia 0.0021
Es mdulo de elasticidad 2000,000 kg/cm2




85
4.6.1 RELACIONES ESFUERZO-DEFORMACION PARA EL CONCRETO Y EL ACERO
DE REFUERZO
Se adopt una funcin que consiste de una parbola y un segmento lineal, propuesta
por Hognestad (21), para idealizar la variacin esfuerzo-deformacin del concreto (ver figura
4.4). La formulacin matemtica de la curva es la siguiente:

=
2
2 "
O
C
O
C
c f fc

0
C

O
( [
O C
c f fc )] = 100 1 "
O
<
C
c f k c f ' 3 " =
4 . 127
'
' 714 . 11 3000
' 5 3900
3
2
c f
c f
c f
k
+
+
=

0.000 0.002 0.004
DEFORMACION C
E
S
F
U
E
R
Z
O


f
c

=
2
"
2
O
C
O
C
C C
f f

( ) [ ]
O C C C
f f = 100 1
"
B
Figura 4.4 CURVA ESFUERZO-DEFORMACION PARA EL CONCRETO
PROPUESTA POR HOGNESTAD

La curva esfuerzo-deformacin idealizada para el acero de refuerzo se presenta en la
figura 4.5. La curva consiste de dos segmentos, para la zona lineal y plstica; el
A
C
"
C
f
o=
86
endurecimiento por deformacin se aproxima por la suposicin de que el elemento consta
de dos componentes, tal como se explic en el captulo III, para el elemento viga-columna.
DEFORMACION S
E
S
F
U
E
R
Z
O


f
s
esy
fy
Figura 4.5 CURVA ESFUERZO-DEFORMACION PARA EL ACERO

En la figura 4.6 se presentan las relaciones tericas momento-curvatura que se
obtuvieron para las secciones rectangulares de las vigas de concreto del prtico P-3, a
partir de los requerimientos de compatibilidad de deformacin y equilibrio de las fuerzas
utilizando las curvas supuestas de esfuerzo-deformacin para el concreto y para el acero
antes descritas, para ms detalle del clculo de estos diagramas remitirse a la referencia
(21).

Tal como se explico en la seccin 3.5.1 al describir el ELEMENTO VIGA-COLUMNA del
programa DRAIN-2D, dicho elemento considera en forma aproximada la interaccin entre
la flexin y la carga axial mediante el diagrama de interaccin de la seccin al cual
denomina superficie de fluencia, en ese sentido ser necesaria la construccin de los
diagramas de interaccin de todas las secciones de las columnas del modelo. En la figura
4.7 se presenta el diagrama de interaccin para la seccin de las columnas del prtico P-3;
as mismo se aclara que el programa DRAIN-2D aproxima la superficie de fluencia mediante
curvas lineales que unen las coordenadas de los puntos de la carga axial (0, Po), punto de
falla balanceada (Mb, Pb), flexin pura (My, 0) y traccin pura (0, To).

87
SECCION VIGA 1er PISO
0.0
5.0
10.0
15.0
20.0
25.0
0.0 20.0 40.0 60.0 80.0 100.0
CURVATURA X10-3(1/m)
M
O
M
E
N
T
O

(
T
-
m
)
My(-)



Figura 4.6 DIAGRAMAS DE MOMENTO-CURVATURA VIGAS DEL PORTICO P-3
My(+)
SECCION VIGA 2do PISO
0.0
5.0
10.0
15.0
20.0
25.0
0.0 20.0 40.0 60.0 80.0 100.0
CURVATURA X10-3(1/m)
M
O
M
E
N
T
O

(
T
-
m
)
My(-)
My(+)
SECCION VIGA 3er PISO
0.0
5.0
10.0
15.0
20.0
25.0
0.0 20.0 40.0 60.0 80.0 100.0
CURVATURA X10-3(1/m)
M
O
M
E
N
T
O

(
T
-
m
)
My(-)
My(+)
0
.
0
5
4#6+1#5
0
.
4
5
0
.
0
5
0.30
2#6
0.30
0
.
0
5
0
.
0
5
4#6
2#6
0
.
4
5
0
.
0
5
0
.
0
5
0
.
4
5
0.30
2#5
2#6+1#4
88

-100
-50
0
50
100
150
200
250
300
350
400
0 5 10 15 20 25 30
MOMENTO (ton-m)
F
U
E
R
Z
A

A
X
I
A
L

(
t
o
n
)
( 0, Po)
( Mb, Pb)
( My , 0)
( 0, To)
CURVA DRAI N- 2D
0
.
0
5
0
.
4
0
6#6
Direccin del sismo
0.40
Figura 4.7 DIAGRAMA DE INTERACCIN PARA LA COLUMNA DEL PORTICO P-3

En los cuadros 4.4 y 4.5 se presentan los aceros requeridos en las secciones crticas
(extremos), as mismo los momentos de fluencia positivo y negativo de los diagramas
momento-curvatura para las vigas, los puntos de los diagramas de interaccin para las
secciones de las columnas de los 19 prticos considerados en el presente estudio.
89
CUADRO 4.4
AREAS DE ACERO EN EXTREMOS DE VIGAS
ACERO SUPERIOR ACERO INFERIOR
PORT

PISO
ARMADURA As
Cm2
My(-)
Ton-m
ARMADURA As
Cm2
My(+)
Ton-m
1 4#6 11.40 14.44 2#5+1#4 5.23 6.90
P-2 2 2#6+1#5 7.68 9.90 2#5 3.96 5.57
1 4#6+1#5 13.38 19.52 2#6 5.70 8.71
P-3 2 4#6 11.40 16.85 2#6 5.70 8.70
3 2#6+1#4 6.97 10.59 2#5 3.96 6.12
1 4#6+1#5 13.38 22.33 2#6+1#4 6.97 12.07
P-4 2 4#6+1#5 13.38 19.68 2#6+1#4 6.97 10.78
3 4#6 11.40 16.90 2#6 5.70 8.70
4 2#6+1#4 6.97 10.59 2#5 3.96 6.12
1-2 5#6 14.25 23.81 2#6+1#4 6.97 12.09
P-5 3 5#6 14.25 20.96 2#6+1#4 6.97 10.60
4 4#6 11.40 16.90 2#6 5.70 8.70
5 3#5+1#4 7.21 10.90 2#5 3.96 6.12
1-2 3#8+1#6 18.06 29.51 3#6 8.55 14.69
3 3#8+1#6 18.06 26.09 3#6 8.55 12.89
P-6 4 2#8+2#6 15.84 22.93 2#6+1#5 7.68 11.63
5 3#6+2#5 12.51 18.43 3#5 5.94 9.05
6 3#5+1#4 7.21 10.88 2#5 3.96 6.12
1-3 3#8+2#6 20.91 38.33 3#6+1#4 9.82 19.24
4 2#8+2#6+1#5 17.82 29.35 2#6+1#5 7.68 13.29
P-7 5 5#6 14.25 23.72 2#6+1#4 6.97 12.00
6 4#6 11.40 19.08 2#5+1#4 5.23 9.06
7 2#5+2#4 6.50 10.72 3#4+1#3 4.52 7.87
1-3 3#8+2#6 20.91 42.42 3#6+1#4 9.82 20.93
4 2#8+2#6 15.84 32.50 2#6+1#5 7.68 16.39
P-8 5 5#6 14.25 26.52 1#6+2#5 6.81 13.17
6 3#6+2#5 12.51 23.36 2#6 5.70 11.03
7 2#6+2#5 9.66 18.21 2#5+1#4 5.23 10.15
8 3#5 5.94 11.47 2#5+1#4 5.23 10.11
1-4 3#8+2#6 20.91 42.59 3#6+1#5 10.53 22.19
5-6 2#8+2#6+1#5 17.82 32.89 2#6+1#5 7.68 15.06
P-9 7 3#6+2#5 12.51 23.45 1#6+2#5 6.81 13.15
8 2#6+2#5 9.66 18.21 2#5+1#4 5.23 10.15
9 2#5+2#4 6.50 12.47 2#5+1#4 5.23 10.12

90
CUADRO 4.4 (continua)
AREAS DE ACERO EN EXTREMOS DE VIGAS
ACERO SUPERIOR ACERO INFERIOR
PORT

PISO
ARMADURA As
cm2
My(-)
Ton-m
ARMADURA As
Cm2
My(+)
Ton-m
1-2 3#8+2#6+1#5 22.89 49.61 4#6 11.40 25.74
3-4 3#8+2#6 20.91 45.54 3#6+1#5 10.53 23.84
5-6 2#8+3#6 18.69 37.13 2#6+2#5 9.66 19.91
P-10 7 1#8+4#6 16.47 33.03 2#6+1#5 7.68 15.91
8 5#6 14.25 28.64 1#6+2#5 6.81 14.19
9 3#6+2#5 12.51 25.22 3#5 5.94 12.37
10 1#6+3#5 8.79 18.14 2#5 3.96 8.35
1-4 3#8+2#6+1#5 22.89 49.61 4#6 11.40 25.74
5-6 3#8+2#6 20.91 41.12 2#6+2#5 9.66 19.84
P-11 7-8 2#8+3#6 18.69 36.97 2#6+1#5+1#4 8.95 18.44
9-10 5#6 14.25 28.60 2#6+1#5 7.68 15.91
11 4#5 7.92 16.34 2#5+1#4 5.23 10.95
1-4 3#8+3#6 23.76 51.11 4#6 11.40 25.79
5-6 3#8+2#6+1#5 22.89 49.35 3#6+1#5 10.53 23.76
P-12 7-8 2#8+4#6 21.54 42.47 2#6+2#5 9.66 19.85
9-10 2#8+2#6+1#5 17.82 35.32 2#6+1#5 7.68 15.94
11 4#6+1#5 13.38 26.99 1#6+2#5 6.81 14.19
12 1#6+3#5 8.79 18.00 3#5 5.94 12.35
1-4 3#8+2#6+1#5 22.89 53.75 3#6+1#5 10.53 25.97
5-8 3#8+2#6+1#5 22.89 44.88 3#6+1#5 10.53 21.69
P-13 9-10 2#8+3#6 18.69 36.92 3#6 8.55 17.66
11-12 4#6+2#5 15.36 30.80 2#6+1#5 7.68 15.92
13 1#6+3#5 8.79 18.00 3#5 5.94 12.35
1-5 3#8+2#6+1#5 22.89 53.92 4#6 11.40 28.00
6-9 3#8+2#6+1#5 22.89 44.88 3#6+1#5 10.53 21.69
P-14 10-11 1#8+5#6 19.32 38.14 2#6+2#5 9.66 19.95
12-13 5#6 14.25 28.60 2#6+1#5 7.68 15.91
14 2#6+2#5 9.66 19.78 3#5 5.94 12.36
1-5 3#8+2#6 20.91 49.57 2#6+2#5 9.66 23.90
6-10 3#8+2#6 20.91 45.18 2#6+2#5 9.66 21.87
P-15 11-12 6#6 17.70 34.07 2#6+1#5 7.68 15.94
13-14 5#6 14.25 28.53 1#6+2#5 6.81 14.20
15 2#6+2#5 9.66 19.78 3#5 5.94 12.36


91
CUADRO 4.4 (continua)
AREAS DE ACERO EN EXTREMOS DE VIGAS
ACERO SUPERIOR ACERO INFERIOR
PORT

PISO
ARMADURA As
Cm2
My(-)
Ton-m
ARMADURA As
Cm2
My(+)
Ton-m
1-5 3#8+2#6+1#5 22.89 53.92 4#6 11.40 28.00
6-10 3#8+2#6 20.91 45.18 2#6+2#5 9.66 21.87
P-16 11-12 2#8+3#6 18.69 36.92 3#6 8.55 17.66
13-15 5#6+1#5 16.23 32.32 2#6+1#5 7.68 15.93
16 2#6+2#5 9.66 19.78 3#5 5.94 12.36
1-6 3#8+2#6+1#5 22.89 53.92 4#6 11.40 28.00
7-12 3#8+2#6+1#5 22.89 49.35 3#6+1#5 10.53 23.76
P-17 13-14 6#6 17.10 34.17 3#6 8.55 17.64
15-16 4#6+2#5 15.36 30.80 2#6+1#5 7.68 15.92
17 2#6+2#5 9.66 19.78 3#5 5.94 12.36
1-6 5#8+2#6 31.05 78.69 5#6+1#5 16.23 42.98
7-12 5#8+2#6 31.05 72.15 5#6 14.25 35.01
P-18 13-14 4#8+2#6 25.98 55.77 3#6+2#5 12.51 28.21
15-17 3#8+2#6 20.91 45.18 2#6+2#5 9.66 21.87
18 5#6 14.25 31.60 2#6+1#5 7.68 17.52
1-7 5#8+3#6 33.90 85.61 6#6 17.10 45.26
8-14 5#8+2#6 31.05 72.47 4#6+2#5 15.36 37.58
P-19 15-16 2#8+5#6 24.39 52.64 4#6 11.40 25.64
17-18 7#6 19.95 43.39 2#6+2#5 9.66 21.86
19 5#6 14.25 31.60 2#6+1#5 7.68 17.52
1-7 5#8+3#6 33.90 85.61 6#6 17.10 45.26
8-14 5#8+2#6 31.05 72.47 4#6+2#5 15.36 37.58
P-20 15-17 3#8+4#6 26.61 56.91 3#6+2#5 12.51 28.22
18-19 2#8+4#6 21.54 46.84 3#6+1#5 10.53 23.74
20 4#6+2#5 15.36 33.77 2#6+1#5 7.68 17.53

92
CUADRO 4.5
AREAS DE ACERO EN EXTREMOS DE COLUMNAS
PORT NIVEL ARMADURA As
Cm2
My
Ton-m
Mb
Ton-m
Pb
Ton
Po
Ton
To
Ton
P-2 1-2 6#6 17.10 11.28 24.83 124.95 354.37 71.82
P-3 1-3 6#6 17.10 11.28 24.83 124.95 354.37 71.82
P-4 1 6#6+2#5 21.06 15.95 36.45 161.97 441.69 88.45
2-4 6#6 17.10 11.28 24.83 124.95 354.37 71.82
P-5 1-2 4#8+2#6 25.98 27.27 59.19 196.35 532.88 109.1
3-5 4#8 20.28 17.34 39.12 160.65 438.56 85.18
P-6 1-2 8#6+2#5 26.76 25.90 59.84 222.62 589.57 112.3
3-6 6#6+2#5 21.06 19.35 45.75 199.91 477.39 88.45
P-7 1-3 6#8 30.42 32.06 71.82 245.44 657.83 127.7
4-7 4#8+2#6 25.98 24.74 54.67 200.81 546.27 109.1
P-8 1-4 6#8+2#6 36.12 34.23 75.81 272.47 734.31 151.7
5-8 4#8+2#6 25.98 21.52 49.29 203.04 546.27 109.1
P-9 1-4 8#8 40.56 39.83 91.36 326.12 859.26 170.3
5-9 6#8+2#5 34.38 33.74 75.81 271.71 727.31 144.4
1-4 6#8+4#6 41.82 47.36 109.1 345.12 917.88 175.6
P-10 5-7 6#8+2#6 36.12 36.99 85.26 291.79 783.39 151.7
8-10 6#8 30.42 31.21 70.42 240.98 657.83 127.7
1-4 8#8+2#6 46.26 49.48 114.5 375.56 998.21 194.2
P-11 5-8 8#8 40.56 38.85 89.94 319.81 859.26 170.3
9-11 6#8+2#6 36.12 34.15 74.40 266.97 734.31 151.7
1-4 10#8 50.70 54.76 125.3 404.12 1078.54 212.9
P-12 5-8 8#8 40.56 38.85 89.94 319.81 859.26 170.3
9-12 6#8+2#6 36.12 34.15 74.40 266.97 734.31 151.7
1-4 10#8+2#6 56.40 75.80 171.8 465.09 1226.41 236.8
P-13 5-8 8#8+2#6 46.26 52.66 119.9 373.68 998.21 194.2
9-13 6#8+2#6 36.12 33.75 78.38 291.07 787.86 151.7
1-5 10#8+4#6 62.10 90.17 206.7 531.73 1387.67 260.8
P-14 6-10 8#8+4#6 51.96 64.73 148.5 433.68 1196.08 218.2
11-14 8#8+2#6 46.26 52.66 119.9 373.68 998.21 194.2
1-5 14#8 70.98 106.3 243.6 605.19 1570.65 298.1
P-15 6-10 12#8 60.84 78.40 178.9 500.27 1315.67 255.5
11-15 10#8 50.70 59.93 140.9 435.71 1141.02 212.9
1-5 16#8 81.12 130.1 294.2 680.82 1767.61 340.7
P-16 6-10 14#8 70.98 99.00 221.5 569.26 1499.25 298.1
11-16 10#8+2#6 56.40 75.80 171.8 465.09 1226.41 236.8

93
CUADRO 4.5 (continua)
AREAS DE ACERO EN EXTREMOS DE COLUMNAS
PORT NIVEL ARMADURA As
Cm2
My
Ton-m
Mb
Ton-m
Pb
Ton
Po
Ton
To
Ton
1-6 18#8 91.26 156.1 350.0 760.88 1973.50 383.2
P-17 7-12 14#8+2#6 76.68 119.5 269.2 640.30 1673.89 322.0
13-17 12#8+2#6 66.54 96.23 213.5 531.73 1405.53 279.4
1-6 28#8 141.9 302.8 671.3 1181.32 3025.27 596.2
P-18 7-12 18#8 91.26 156.1 350.0 760.88 1973.50 383.2
13-18 14#8 70.98 106.3 243.6 605.19 1570.65 298.1
1-7 30#8 152.1 326.6 734.2 1292.58 3280.25 638.8
P-19 8-14 20#8 101.4 172.7 391.3 851.38 2188.32 425.8
15-19 16#8 81.12 119.8 276.5 686.46 1767.61 340.7
1-7 32#8 162.2 366.8 827.5 1401.75 3544.15 681.4
P-20 8-14 22#8 111.5 202.4 456.9 940.49 2412.06 468.4
15-20 18#8 91.26 145.2 331.2 766.71 1973.50 383.2


4.7 REGISTROS DE ACELERACION CONSIDERADOS
En un anlisis tiempo historia no lineal, uno de los principales problemas consiste en la
seleccin del tipo y el nmero de los registros de aceleracin a usar, adems en la
seleccin debemos tener en cuenta las magnitudes que deben presentar, el contenido de
energa, la duracin, cmo tomar en cuenta la atenuacin con la distancia, etc. De acuerdo
a la gua de diseo FEMA 273 (9), se requieren tres registros de aceleracin para tomar la
respuesta estimada como la envolvente de las respuestas; otro caso es cuando se tienen
ms de seis registros y en esta situacin se acepta como la demanda estimada al promedio
de las demandas registradas. La NTE E-030 (26), especifica que se requieren cinco
registros, escalados al mximo valor esperado de la aceleracin en el sitio, sin embargo no
indica qu respuesta tomar.
En la presente investigacin haremos uso de seis registros de aceleracin horizontal,
que corresponden a tres sismos que se han producido en nuestro pas entre 1966 y 1974.
Estos registros han servido de base para la NTE E-030, y sus caractersticas se presentan
en el cuadro 4.6 (2):

94
CUADRO 4.6
REGISTROS DE SISMOS
COD FECHA DENOMINACION PGA PGV PGD M
g cm/seg cm
A 17/10/66 Comp. N82W 0.184 13.23 7.35 7.5
B
17/10/66
Comp. N08E 0.275 -21.60 -16.60 7.5
C
31/05/70
Comp. L 0.107 4.71 1.55 7.7
D
31/05/70
Comp. T 0.100 6.98 2.64 7.7
E
03/10/74
1421 GCT Comp. N08E 0.182 10.30 -5.34 7.5
F
03/10/74
1421 GCT Comp. N82W 0.196 14.48 6.41 7.5


Segn la moderna definicin de objetivos del diseo sismorresistente de edificios, se ha
convenido en establecer cuatro niveles de severidad en las solicitaciones ssmicas, cada
uno de los cuales se define por un sismo de diseo. Dado que los terremotos son tratados
como sucesos aleatorios propios de cada regin, los sismos de diseo se definen en
funcin de los periodos medios de retorno de estos eventos o en funcin de la probabilidad
de excedencia durante un determinado tiempo de exposicin. En nuestro pas se han
estimado las aceleraciones asociadas a los cuatro sismos de diseo para la costa del Per.
El cuadro 4.7 muestra los periodos de retorno para estos sismos, junto a la aceleracin
mxima esperada en la costa del Per (18).

CUADRO 4.7
SISMOS DE DISEO, PERIODOS DE RETORNO Y ACELERACIONES
Sismos Sismos Sismos Sismos
Frecuentes Ocasionales Raros Muy raros
Periodo de retorno 45 75 500 1000
Aceleracin en la roca
para la costa del pas (g)
0.2 0.25 0.4 0.5


La Norma Peruana de Diseo Sismorresistente de 1997 (26), tiene como sismo mximo
de diseo, un evento de 500 aos de periodo de retorno, por tal motivo exige que los
registros usados en un anlisis tiempo historia no lineal, reflejen valores de aceleracin
95
compatibles con la zonificacin dispuesta en ella. Esto significa que los acelerogramas
deben tener valores picos de 0.40g de PGA (Peak Ground Aceleration, mxima aceleracin
del suelo o factor Z) para la costa del Per. Sin embargo, los registros disponibles de
sismos Peruanos no cumplen con este requisito, por lo que la norma nos indica que deben
ser escalados, es decir sus ordenadas deben multiplicarse por un valor dado de tal modo
que el pico de la aceleracin tenga el valor estipulado. En la presente investigacin todos
los sismos fueron normalizados a una aceleracin mxima igual a 0.4g de aceleracin
mxima.
En cuanto a la duracin, sta no influye mayormente en las respuestas globales, pero s
en las deformaciones inelsticas locales (6), por tal motivo en concordancia con los
objetivos de la presente investigacin, se han considerado los registros totales de los
sismos mencionados para los anlisis correspondientes, estos registros se muestran en las
figuras 4.8.


96
-300
-200
-100
0
100
200
300
0 10 20 30 40 50 60 70
TIEMPO (seg)
A
C
E
L
E
R
A
C
I
O
N

(
c
m
/
s
e
g
2
)
a) Sismo del 13-Oct.-66, componente N82W
-300
-200
-100
0
100
200
300
0 10 20 30 40 50 60 70
TIEMPO (seg)
A
C
E
L
E
R
A
C
I
O
N

(
c
m
/
s
e
g
2
)
b) Sismo del 13-Oct.-66, componente N08E
-150
-100
-50
0
50
100
150
0 10 20 30 40 5
TIEMPO (seg)
A
C
E
L
E
R
A
C
I
O
N

(
s
e
g
/
c
m
2
)
0
c) Sismo del 31-May.-70, componente L
97
-150
-100
-50
0
50
100
150
0 10 20 30 40 5
TIEMPO (seg)
A
C
E
L
E
R
A
C
I
O
N

(
c
m
/
s
e
g
2
)
0
d) Sismo del 31-May.-70, componente T
-200
-150
-100
-50
0
50
100
150
200
0 20 40 60 80 1
TIEMPO (seg)
A
C
E
L
E
R
A
C
I
O
N

(
c
m
/
s
e
g
2
)
00
e) Sismo del 3-Oct.-74, 1421 GCT componente N08E
-200
-150
-100
-50
0
50
100
150
200
0 20 40 60 80 100
TIEMPO (seg)
A
C
E
L
E
R
A
C
I
O
N

(
c
m
/
s
e
g
2
)
f) Sismo del 3-Oct.-74, 1421 GCT componente N82W
Figura 4.8 GRAFICOS DE REGISTROS DE SISMOS
98

4.8 AMORTIGUAMIENTO
En la mayora de los anlisis de respuesta dinmica de estructuras de edificios, todas
las diversas fuentes de amortiguamiento se representan por un amortiguamiento viscoso,
cuando se van a analizar estructuras de edificios de gran altura para determinar su
respuesta ssmica, se utilizan valores del factor de amortiguamiento entre 0.03 y 0.05 segn
la prctica Japonesa para estructuras de concreto reforzado. Se supone que los factores
de amortiguamiento correspondientes para los modos ms altos se incrementan en
proporcin a las frecuencias naturales (16), en tal sentido se ha adoptado 0.05 para el
primer modo, y 0.10 para el segundo modo, para calcular los coeficientes para determinar la
matriz de amortiguamiento mediante una combinacin lineal con la matriz de masa y la
matriz de rigidez tangencial instantnea de la estructura; tal como se explic en la seccin
3.4.2. En el cuadro 4.8, se presentan los valores de los coeficientes dependientes de la
masa y de la rigidez calculada con las expresiones 3.13a y 3.13b para ensamblar la matriz
de amortiguamiento, para los distintos prticos considerados en el presente estudio.
Se puede especificar otro tipo de amortiguamiento denominado amortiguamiento
estructural, llamado as debido a que el amortiguamiento es proporcional a las
deformaciones de la estructura, pero opuesto a las tasas de deformacin, este tipo de
amortiguamiento puede proporcionar esencialmente relaciones constantes de
amortiguamiento modal para todos los modos de vibracin de una estructura elstica, en la
presente investigacin no se utilizar este tipo de amortiguamiento. El procedimiento para
aplicar los efectos del amortiguamiento estructural es el siguiente, no se realiza ninguna
modificacin en la rigidez de la estructura, en cambio, se calculan un conjunto de fuerzas de
amortiguamiento al final de cada paso de integracin y se agregan a las cargas en el
siguiente paso. Las fuerzas de amortiguamiento se calculan en el nivel de elemento y se
oponen a las deformaciones del mismo; debido a que las fuerzas de amortiguamiento que
actan durante el paso n+1 se basan en las condiciones al final del paso n, puede existir
una tendencia de oscilaciones numricas desarrolladas en el proceso de integracin para
ser amplificadas con un error numrico subsiguiente (24). Esta es la razn principal del por
qu este tipo de amortiguamiento debera ser especificado con precaucin.
99

CUADRO 4.8
COEFICIENTES DE AMORTIGUAMIENTO
PERIODO AMORTIGUAMIENTO COEFICIENTES
PORTICO
MODO1 MODO2 1 2
P-2 0.37 0.10 0.05 0.10 0.8417 0.0030
P-3 0.49 0.15 0.05 0.10 0.5486 0.0045
P-4 0.57 0.18 0.05 0.10 0.4511 0.0054
P-5 0.61 0.20 0.05 0.10 0.3973 0.0060
P-6 0.74 0.24 0.05 0.10 0.3334 0.0072
P-7 0.83 0.28 0.05 0.10 0.2779 0.0084
P-8 0.86 0.29 0.05 0.10 0.2684 0.0087
P-9 0.93 0.31 0.05 0.10 0.2534 0.0093
P-10 1.12 0.38 0.05 0.10 0.2038 0.0114
P-11 1.24 0.42 0.05 0.10 0.1846 0.0125
P-12 1.35 0.46 0.05 0.10 0.1677 0.0137
P-13 1.37 0.47 0.05 0.10 0.1632 0.0140
P-14 1.41 0.49 0.05 0.10 0.1546 0.0147
P-15 1.44 0.49 0.05 0.10 0.1576 0.0146
P-16 1.53 0.52 0.05 0.10 0.1487 0.0155
P-17 1.57 0.54 0.05 0.10 0.1417 0.0161
P-18 1.82 0.62 0.05 0.10 0.1245 0.0185
P-19 1.86 0.63 0.05 0.10 0.1231 0.0188
P-20 1.95 0.66 0.05 0.10 0.1176 0.0197

4.9 ANALISIS INELASTICO DE LOS PORTICOS
Para el anlisis inelstico de las edificaciones aporticadas, se usar el modelo de rtula
puntual (modelo discreto), dentro del cual podemos considerar el modelo del prtico
bidimensional, el cual se usa para idealizar edificios regulares que tienen pequeas
excentricidades.
En ste modelo, todos los prticos en una direccin, son tratados como elementos
planos unidos por vigas rgidas que representan la condicin de rigidez infinita de la losa de
piso. En este modelo se puede tener la respuesta ssmica a nivel de cada elemento y por
ende a nivel de cada piso.
Se efectuaron los anlisis inelsticos de los prticos, usando el programa de cmputo
DRAIN-2D (15,24), que es un programa general de anlisis dinmico inelstico de
100
estructuras, se aplicaron a los modelos cargas verticales de gravedad, el programa slo
permite cargas nodales y stas fueron aplicadas en los nudos de todos los niveles para
representar las cargas verticales en las columnas. Se uso el ELEMENTO VIGA-COLUMNA
para modelar las columnas y vigas de los prticos en estudio. Sus propiedades
geomtricas son las de sus secciones no fisuradas y el mdulo de elasticidad el del
concreto. Se consider endurecimiento por deformacin de 5% para todos los elementos.
Las propiedades del concreto, del acero, registros de aceleracin y amortiguamiento
fueron definidas en las secciones previas.
Se consideraron algunas simplificaciones en el anlisis para evitar soluciones
complicadas y costosas, estas suposiciones fueron los siguientes:
Las vigas y columnas son elementos lineales sin masas, en el caso del ELEMENTO
VIGA-COLUMNA consta de dos componentes una elastoplstico y otro elstico en
paralelo, y en caso del ELEMENTO VIGA CON RIGIDEZ DEGRADANTE consta de una
parte elstica con dos resortes de rotacin no lineales en sus extremos, ambos
elementos fueron descritos detalladamente en el capitulo III.
Se permite la deformacin axial en las columnas. Sin embargo se despreciaron las
deformaciones axiales en todas las vigas, para considerar la hiptesis de losa rgida.
Se considera la interaccin fuerza axialmomento para todas las columnas, asumiendo
superficies de interaccin para secciones de concreto armado.
La estructura es un prtico plano el cual se desplaza horizontalmente en su plano, y rota
alrededor de un eje perpendicular al plano de la estructura.
Se consideran que las deformaciones son suficientemente pequeas, para permitir que
la configuracin inicial de la estructura prevalezca durante todo el anlisis.
Los centros de los nudos en una conexin viga-columna, son infinitamente rgidas.
Las caractersticas de rigidez de la estructura, permanecen constantes durante cada
incremento de tiempo corto.
Las masas son concentradas y localizadas, donde se definen los grados de libertad
horizontal, y corresponden a la carga muerta ms un 25% de la carga viva.
Se considera la cimentacin de la estructura infinitamente rgida. Las columnas en el
primer piso son rgidamente conectadas a su cimentacin.
101
Se considera el efecto P-, mediante la adicin de una rigidez geomtrica a la rigidez
de la columna, usando la fuerza axial producida por las cargas estticas.
Los movimientos ssmicos en la base ocurren en el plano de la estructura, en la
direccin horizontal.
Se consideran las cargas de gravedad, estn incluidas como fuerzas de empotramiento
en las vigas y estn basadas en las cargas muertas y vivas. Adicionalmente se aplican
cargas nodales para representar las cargas de los muros y el peso muerto de las
columnas. No se han considerado los efectos por reduccin de carga viva.
Por razones de capacidad del programa de cmputo, se han tomado distintas longitudes
para el intervalo de integracin para los anlisis de los prticos, para los prticos del P-2
hasta el P-10 se ha usado un cuarto del intervalo de digitacin (0.02 seg.) del registro
ssmico, es decir, 0.005 seg., para los prticos del P-11 hasta el P-17 se ha tomado
como intervalo de integracin 0.010 seg., y para los prticos del P-18 hasta el P-20 se
tomado como intervalo 0.020 seg.; sin embargo el tamao de los intervalos de
integracin en especial para estos ltimos prticos, estn dentro del rango que se
recomienda, generalmente tT/10 (5), para satisfacer las condiciones de estabilidad
en los clculos.

4.10 RESULTADOS DE LOS ANALISIS INELASTICOS
De acuerdo a los objetivos de la presente investigacin, los resultados analticos que se
requerirn, son la historia de los momentos flexionantes en el tiempo de los elementos,
especialmente en aquellos instantes donde ocurren eventos inelsticos en la estructura; un
evento significa la formacin o cierre de una rtula plstica en uno de los extremos del
elemento, enfatizando en el estudio de las distribuciones inesperadas de momentos
flexionantes en las columnas. Las respuestas inelsticas tiempo-historia de los 19 prticos,
se presentarn en el captulo V, conjuntamente para desarrollar el factor de distribucin de
momentos (FDM) en un determinado nudo.







CAPITULO V


DESARROLLO DEL FACTOR DE
DISTRIBUCION DE MOMENTOS


5.1 INTRODUCCION
Los anlisis tiempo-historia no lineales, indican que en las columnas de los prticos de
varios pisos, pueden ocurrir distribuciones inesperadas de momento flexionante, en
comparacin con la distribucin obtenida de la carga lateral esttica del cdigo. El anlisis
de la carga lateral esttica indica que generalmente existen puntos de inflexin prximos a
la mitad de la altura de las columnas, a menos que las vigas sean mucho ms flexibles que
las columnas, excepto en los pisos prximos a la parte superior e inferior del prtico (ver
figura 5.1, momentos estticos). Sin embargo los anlisis tiempo-historia no lineal, sugiere
que en determinados momentos durante la respuesta de la estructura a los movimientos
ssmicos, el punto de inflexin en una columna puede estar prximo a la unin viga-
columna, y ocasionalmente incluso la columna puede estar en curvatura simple.
103

nivel
1
0
esttico
3
2
5
4
6
19.90 seg 18.15 seg 20.70 seg
Figura 5.1 MOMENTOS FLEXIONANTES EN LAS COLUMNAS EXTERIORES
DEL PORTICO P-6, QUE RESPONDE NO LINEALMENTE
AL SISMO 10/03/74 COMP N82W
En la figura 5.1, se muestran algunos diagramas de momento flexionante para algunos
instantes, de una columna para un prtico de 6 pisos (P-6), diseado de acuerdo al
Reglamento Nacional de Construcciones (26) y respondiendo al sismo del 10/03/74 en Lima
componente N82W. La causa de la distribucin inesperada de los momentos flexionantes
en la columna en algunos instantes es la fuerte influencia de los modos superiores de
vibracin, especialmente el segundo y tercero (21).
La carga lateral esttica del cdigo normalmente tiene una distribucin triangular que
corresponde a las cargas laterales que varan linealmente desde cero en la base de la
columna a un mximo en la parte superior de la estructura, a veces con una carga
concentrada adicional en la parte superior. Esta distribucin de cargas corresponde
aproximadamente al modo fundamental (primero) de vibracin. Si son significativos los
otros modos de vibracin, es evidente que la distribucin de las fuerzas laterales en la
estructura en determinados instantes podra diferir apreciablemente de la supuesta, y
producir por ejemplo la distribucin de momentos flexionantes mostrada en la figura 5.1. El
corrimiento del punto de inflexin en las columnas a un punto bien distante de la posicin a
mitad de la altura, en algunos casos significa que los momentos inducidos en la columna

104
pueden ser mucho mayores que el momento obtenido de un anlisis de carga lateral
esttica y puede conducir a que se formen rtulas plsticas en la columnas.

5.2 CONTROL DE ROTULAS PLASTICAS
Diferentes cdigos (3, 25, 27), tratan de hacer que las rtulas plsticas se formen en las
vigas antes que en las columnas, requiriendo que la suma de las resistencias de momentos
de las columnas sea mayor que la suma de las resistencias de las vigas en cada plano
principal de una conexin. Desafortunadamente, ste requerimiento no impide la formacin
de rtulas plsticas en las columnas por las siguientes razones, como se ver a
continuacin.
Se ha encontrado como se indic anteriormente, en anlisis tiempo-historia no lineal de
prticos de varios niveles, que en distintas etapas durante un sismo pueden ocurrir puntos
de inflexin bastante alejados de las alturas medias de las columnas. Las columnas
pueden incluso estar en curvatura simple en ciertos instantes. En consecuencia son
posibles las distribuciones de momentos flexionantes tales como se indican en la Figura 5.2


Figura 5.2 MOMENTOS EN UNION VIGA-COLUMNA
a) Parte del prtico
Mc1
Mc2
b) Momento flexionante de la columna
Mc2 Mc1
Mv2
Mv1
En una unin tpica, la suma de momentos de la columna resiste los momentos totales de
las vigas, por lo tanto se tiene:

105

+ = + = 2 Mc 1 Mc 2 Mv 1 Mv Mv (5.1)
El mayor momento de la columna est dado por:

= 2 Mc Mv 1 Mc (5.2)
Si:
Muv =Es el momento total cuando las vigas estn en la capacidad del momento mximo y
Muc1= Es la capacidad del momento mximo de la columna 1
Si se desea impedir la formacin de una rtula plstica en la columna 1, el requerimiento es:
(5.3) 2 Mc Muv 1 Muc >
Si la columna permanece en doble curvatura, el caso lmite ocurre cuando Mc20, y la
ecuacin anterior requiere que:
(5.4) Muv 1 Muc >
Si las columnas estn en curvatura simple (el punto de inflexin est fuera del entrepiso), la
ecuacin 5.3, requiere que:
(5.5) 2 Mc Muv 1 Muc + >
En consecuencia, el requerimiento 12.4.3 de la Norma E-060 de Concreto Armado (27), de
que la suma de las resistencias de columna exceda la suma de las resistencias de la viga
en una unin, no impide la formacin de rtulas plsticas en las columnas. Para asegurar
que no se formen rtulas plsticas en las columnas, se tendra que exigir:
Que la resistencia a flexin de cada seccin de columna fuera al menos igual a la suma
de las resistencias a flexin de las secciones de la viga en el plano de flexin, si el
punto de inflexin est en cualquier parte dentro de la altura del entrepiso.
Si el punto de inflexin est fuera de la altura del entrepiso, se requerira una capacidad
todava mayor de la columna.
El objetivo de la presente investigacin es determinar un procedimiento que contemple
estos casos.


5.3 REVISION DE REGLAMENTOS RESPECTO AL CONTROL DE ROTULAS
PLASTICAS

106
Se describen algunos reglamentos, que contienen algunos requerimientos sobre el
control de rtulas plsticas en nudos.

5.3.1 NORMA NTE-060 CONCRETO ARMADO
En la Norma Tcnica de Edificacin E-060 Concreto Armado (27), dentro de sus
disposiciones especiales para columnas sujetas a flexocompresin que resistan fuerzas de
sismo, se especifica en el artculo 12.4.3 lo siguiente:
La resistencia a la flexin de las columnas debern satisfacer la ecuacin:
(5.6)

> Mnv 4 . 1 Mnc
Donde:
Mnc es la suma de momentos, al centro del nudo, correspondiente a la resistencia
nominal en flexin de las columnas que forman dicho nudo; esta resistencia en
flexin se calcular para la fuerza axial actuante en la hiptesis que considera las
fuerzas de gravedad y de sismo en la direccin considerada, verificando la
condicin que d como resultado la resistencia a flexin ms baja.
Mnv es la suma de momentos, al centro del nudo, correspondiente a las resistencias
nominales en flexin de las vigas que forman el nudo. Las resistencias a la flexin
debern sumarse de manera que los momentos de la columna se opongan a los
momentos de las vigas. Deber satisfacerse la ecuacin anterior para momentos
de viga que actan en ambas direcciones del plano vertical del prtico que se
considera.

5.3.2 NORMA ACI 318 99
La norma ACI (3), dentro de sus disposiciones especiales para el diseo ssmico, para
elementos sometidos a flexin y carga axial pertenecientes a prticos especiales resistentes
a momento, establece en su artculo 21.4.2.2 lo siguiente:

Las resistencias a flexin de las columnas deben satisfacer la ecuacin (21-1).


> Mg
5
6
Mc (21-1)

107
Mc= suma de los momentos, en las caras del nudo, correspondiente a la resistencia
nominal a flexin de las columnas que confluyen en dicho nudo. La resistencia a la flexin
de la columna debe calcularse para la carga axial factorizada, consistente con la direccin
de las fuerzas laterales consideradas, que d la ms baja resistencia a la flexin.
Mg= suma de los momentos en las caras del nudo correspondiente a la resistencia
nominal a flexin de las vigas que llegan a dicho nudo. Las resistencias a la flexin deben
sumarse de tal forma que los momentos de la columna se opongan a los momentos de la
viga. Debe satisfacerse la ecuacin (2-1) para momentos de vigas que acten en ambas
direcciones en el plano vertical del marco que se considera.

5.3.3 REGLAMENTO DE CONSTRUCCIONES PARA EL DISTRITO FEDERAL-
MEXICO
La norma para el Distrito Federal de Mxico (25), dentro de sus disposiciones
especiales para prticos dctiles, establece en su artculo 7.3.2.1 lo siguiente:
Las resistencias a la flexin de las columnas en un nudo deben satisfacer la Ec. 7.2
(7.2)

> Mg 5 . 1 Me
Donde:
Me suma al pao del nudo de los momentos resistentes calculados con un factor de
resistencia igual a uno, de las columnas que llegan a ese nudo; y
Mnv suma al pao del nudo de los momentos resistentes calculados con un factor de
resistencia igual a uno, de las vigas que llegan al nudo.
Las sumas anteriores deben realizarse de modo que los momentos de las columnas se
opongan a los de las vigas. La condicin debe cumplirse para los dos sentidos en que
puede actuar el sismo. No ser necesario cumplir con la Ec. 7.2 en los nudos de azotea

5.4 FACTOR DE DISTRIBUCION DE MOMENTOS (FDM)
Debido a la distribucin desproporcionada de los momentos alrededor de las uniones de
viga-columna en la respuesta ssmica de un prtico de varios pisos, podran aparecer
momentos flexionantes en las secciones crticas, apreciablemente mayores que los
deducidos del anlisis esttico. Para reducir la probabilidad de fluencia de columnas, se

108
debe hacer que la resistencia nominal de la columna en una seccin crtica (por ejemplo,
sobre el nivel del piso en la Figura 5.2) sea mayor que la suma total de momentos mximos
por las vigas adjuntas cuando se forman rtulas plsticas. Es decir que se requiere:
(5.7) Muv ) FDM ( 1 Muc
1 C

Donde:
Muc1= Es la capacidad del momento mximo de la columna 1
Muv= Suma de las capacidades mximas a flexin de las vigas adyacentes a la
columna 1 cuando se forman rtulas plsticas en las vigas
(FDM)
C1
= Factor de distribucin de momentos, dependiente de la respuesta dinmica
inelstica del prtico cuando se sujeta a movimientos ssmicos.

De la ecuacin 5.7 se tiene:

Muv
1 Muc
FDM
1 C

= (5.8)

El procedimiento para calcular el FDM en una determinada seccin, de una conexin viga-
columna, de un prtico en particular, sometido a un determinado sismo, involucra los
siguientes pasos:
Generar mediante un anlisis dinmico inelstico, los resultados tiempo-historia de los
momentos flexionantes en los extremos de las vigas y columnas que forman un nudo,
de todos los nudos del prtico, sometidos a los sismos seleccionados y de acuerdo a
las suposiciones y parmetros indicados en el captulo IV.
Para cada intervalo de integracin, verificar en todos los nudos del prtico, si se ha
producido una rtula plstica en los extremos del elemento que forman un determinado
nudo. De ocurrir esto, calcular los FDM arriba y abajo del nivel del piso en dichos nudos
para ese instante.
Criterio de seleccin del FDM que se producirn en el prtico, para un determinado
instante y durante todo el sismo.


109
La siguientes tres secciones describen en detalle los procedimientos asociados con los
pasos antes indicados.

5.4.1 RESULTADOS TIEMPO-HISTORIA DE MOMENTOS FLEXIONANTES
Se han generado la historia en el tiempo durante todo el sismo, de los momentos
flexionantes en los extremos de todos los elementos para los 19 prticos, cada prtico ha
sido sometido a 6 sismos seleccionados, escalados y descritos, conjuntamente con los
parmetros asumidos para el anlisis dinmico inelstico, en el capitulo IV. En la figura 5.3
se muestra un diagrama de momentos flexionantes tpico en un determinado instante, para
el prtico P-5, que responde no linealmente a un determinado sismo; nuevamente se
pueden apreciar distribuciones de momentos inesperados en los ejes de las columnas
interiores y exteriores.
En los 114 anlisis que se efectuaron (19 prticos, sometidos a 6 sismos cada uno),
todos los prticos incursionaron en el rango inelstico, para los niveles de aceleracin
seleccionados previamente (0.4g).
Los anlisis tiempo-historia y el manejo de la base de datos de respuesta para su
posterior procesamiento, para todo el conjunto de prticos, requirieron un gran esfuerzo en
tiempo de cmputo y almacenamiento. Para el manejo de los resultados, se elabor un
programa de cmputo el cual lea los datos de respuesta de los archivos binarios emitidos
por el programa DRAIN-2D, y calculaba el FDM en todos los nudos del prtico, de acuerdo
al algoritmo que se describir en la siguiente seccin.


110
30 29 28 27 26
Figura 5.3 MOMENTOS FLEXIONANTES EN EL PORTICO P-5
QUE RESPONDE NO LINEALMENTE AL SISMO
10/03/74 COMP N82W
1
t=17.90 seg
2
b) vigas
3
6
11
16
21
7
12
8
13
17
22
18
23
4 5
9
14
10
15
19
24 25
20
1 nmero del nudo
6
1
a) columnas
7
2
8
3
21
16
11
26
23 22
17
12
18
13
27 28
9
4
10
5
25 24
19
14
20
15
29 30
rtula plstica





Se repiti el procedimiento antes descrito para los 19 prticos.
Para un prtico en particular, se ejecuta el mismo procedimiento, sometindole a los 6
sismos seleccionados, consiguiendo 6 valores de FDM que representan al prtico en
estudio.
De acuerdo a los dos pasos anteriores, se obtienen los FDM para cada instante durante
todo el sismo. De todo el conjunto de estos valores, se selecciona el promedio y es el
valor que representar al prtico, para un sismo en particular. En el cuadro 5.1 se
muestra la historia de los clculos que se efectuaron para el prtico P-5, sometido a los
seis sismos seleccionados previamente.
De acuerdo al paso anterior, se obtienen todos los FDM para todos los nudos del
prtico, donde se incursion en el rango inelstico. Se selecciona el mayor de todos los
FDM de estos nudos, y es el valor que representar al prtico, para ese instante.
Para un determinado instante correspondiente a un sismo en particular, se obtiene los
momentos en todos los nudos del prtico para ese instante, se calculan los FDM para la
seccin superior e inferior de la columna, en aquellos nudos donde las vigas entraron
en fluencia, se selecciona el mayor FDM sea superior o inferior y es el valor que
representar a ese nudo. Debe observarse, si el FDM es menor que 1, el punto de
inflexin se encuentra dentro del entrepiso y la columna estar en curvatura doble; si es
mayor que 1, el punto de inflexin estar fuera del entrepiso y la columna estar en
curvatura simple.
5.4.2 CALCULO DEL FACTOR DE DISTRIBUCION DE MOMENTOS
De los resultados tiempo-historia de los momentos flexionantes de un determinado
prtico, se calcula el FDM de acuerdo a la ecuacin 5.8, en cada nudo del prtico, para
cada instante de la respuesta ssmica, de acuerdo al siguiente algoritmo:
111


CUADRO 5.1
FACTOR DE DISTRIBUCION DE MOMENTOS (FDM) EN EL PORTICO P-5, NUDOS INTERIORES, 0.4G
SISMO SISMO SISMO SISMO SISMO SISMO
10/17/66 COMP N82W 10/17/66 COMP NO8E 05/31/70 COMP L 05/31/70 COMP T 10/03/74 COMP NO8E 10/03/74 COMP N82W
TIEMPO FDM NUDO TIEMPO FDM NUDO TIEMPO FDM NUDO TIEMPO FDM NUDO TIEMPO FDM NUDO TIEMPO FDM NUDO
20.20 0.65 14 19.65 0.53 17 3.40 0.79 24 1.80 0.72 22 7.40 0.77 23 9.85 0.80 24
20.25 0.86 24 20.00 0.83 24 3.70 0.63 17 3.95 0.71 8 7.65 0.67 14 9.90 0.83 22
20.30 0.80 23 20.05 0.85 24 3.95 0.50 13 4.00 0.88 22 8.00 0.63 13 9.95 0.82 22
20.45 0.85 22 20.10 0.78 8 4.00 0.69 8 4.25 0.58 13 11.30 0.67 19 10.30 0.76 13
20.50 0.87 7 22.45 0.51 13 4.25 0.77 23 4.30 0.67 18 11.35 0.87 22 10.35 0.72 13
20.55 0.81 22 22.50 0.65 12 4.75 0.59 7 4.35 0.61 7 11.65 0.69 12 10.40 0.78 23
20.90 0.54 17 22.90 0.76 8 4.95 0.60 7 4.60 0.67 8 11.70 0.74 19 12.25 0.62 9
22.30 0.67 12 22.95 0.76 7 5.00 0.83 23 4.65 0.85 9 16.60 0.62 19 12.30 0.67 8
22.65 0.72 8 23.30 0.56 14 5.05 0.83 7 4.70 0.89 22 23.60 0.65 8 12.55 0.82 23
22.70 0.53 13 23.65 0.79 22 5.35 0.64 13 4.75 0.75 7 23.65 0.56 13 12.60 0.89 24
23.45 0.73 7 23.70 0.64 19 5.40 0.80 23 5.00 0.73 22 24.00 0.61 14 12.65 0.96 23
23.75 0.56 13 23.95 0.57 14 5.45 0.83 23 5.05 0.90 22 24.05 0.56 13 12.70 0.92 24
23.80 0.58 9 24.00 0.56 9 5.70 0.56 13 5.10 0.86 9 24.40 0.56 12 12.90 0.54 12
23.85 0.60 13 24.40 0.50 13 5.75 0.71 7 5.15 0.81 9 30.00 0.56 14 12.95 0.60 9
24.25 0.57 13 5.80 0.73 7 5.40 0.69 23 34.75 0.59 18 13.00 0.74 9
24.60 0.64 18 6.05 0.65 8 5.45 0.59 13 41.20 0.60 19 13.05 0.66 17
6.10 0.73 8 5.50 0.81 13 45.85 0.76 9 17.20 0.59 8
6.15 0.78 8 5.55 0.78 13 45.90 0.71 18 17.25 0.64 8
6.45 0.65 8 5.60 0.72 17 46.15 0.67 7 17.55 0.64 13
6.50 0.80 7 5.75 0.76 12 46.20 0.77 8 17.80 0.53 8
6.55 0.89 22 5.80 0.79 9 47.95 0.60 13 17.85 0.77 18
6.85 0.62 13 5.85 0.85 9 48.00 0.55 14 17.90 0.93 24
7.30 0.80 22 5.90 0.66 7 48.30 0.76 8 18.05 0.68 9
7.60 0.67 13 6.20 0.70 13 48.35 0.69 19 18.10 0.92 9
8.35 0.57 9 6.25 0.66 8 48.65 0.53 13 18.15 0.89 9
8.70 0.60 8 6.30 0.64 8 49.40 0.61 13 19.20 0.65 19
8.75 0.58 13 6.35 0.59 13 50.45 0.53 12 19.25 0.75 19
17.70 0.70 8 6.55 0.61 18 50.80 0.59 8 19.85 0.71 8
17.75 0.64 18 6.60 0.89 22 50.85 0.62 8 20.10 0.73 23
18.10 0.60 13 9.90 0.64 7 20.15 0.51 12
18.85 0.62 14 10.35 0.76 24 20.20 0.60 7
contina .
FACTORES DE DISTRIBUCION DE MOMENTOS MAXIMOS EN CADA INSTANTE

112

CUADRO 5.1 (contina .)
FACTOR DE DISTRIBUCION DE MOMENTOS (FDM) EN EL PORTICO P-5, NUDOS INTERIORES, 0.4G
SISMO SISMO SISMO SISMO SISMO SISMO
10/17/66 COMP N82W 10/17/66 COMP NO8E 05/31/70 COMP L 05/31/70 COMP T 10/03/74 COMP NO8E 10/03/74 COMP N82W
TIEMPO FDM NUDO TIEMPO FDM NUDO TIEMPO FDM NUDO TIEMPO FDM NUDO TIEMPO FDM NUDO TIEMPO FDM NUDO
19.50 0.53 14 27.20 0.63 14 20.55 0.69 8
19.55 0.56 18 20.60 0.85 9
19.90 0.53 14 20.90 0.61 13
19.95 0.83 8 20.95 0.68 19
23.05 0.60 13 21.30 0.55 12
23.10 0.60 13 30.55 0.65 12
29.00 0.55 14 30.60 0.61 12
29.30 0.52 14 30.90 0.78 9
29.35 0.86 8 31.30 0.69 12
29.65 0.53 14 31.35 0.60 13
31.10 0.60 13 31.95 0.61 19
31.15 0.63 14 33.45 0.54 12
31.45 0.53 14 33.75 0.53 12
31.85 0.59 14 33.80 0.76 9
32.20 0.78 8 33.85 0.79 9
32.50 0.53 14 36.50 0.59 13
32.55 0.53 14 36.55 0.63 8
36.90 0.55 12
39.95 0.65 13
40.00 0.74 19
40.25 0.53 13
40.30 0.53 12
40.35 0.65 9
43.15 0.70 13
43.20 0.72 8
43.50 0.64 13
43.55 0.73 8
43.60 0.80 9
45.05 0.75 9
45.10 0.77 9
45.40 0.60 12
contina .
FACTORES DE DISTRIBUCION DE MOMENTOS MAXIMOS EN CADA INSTANTE

113

114
CUADRO 5.1 (contina .)
FACTOR DE DISTRIBUCION DE MOMENTOS (FDM) EN EL PORTICO P-5, NUDOS INTERIORES, 0.4G
SISMO SISMO SISMO SISMO SISMO SISMO
10/17/66 COMP N82W 10/17/66 COMP NO8E 05/31/70 COMP L 05/31/70 COMP T 10/03/74 COMP NO8E 10/03/74 COMP N82W
TIEMPO FDM NUDO TIEMPO FDM NUDO TIEMPO FDM NUDO TIEMPO FDM NUDO TIEMPO FDM NUDO TIEMPO FDM NUDO
45.45 0.76 24
45.50 0.64 7
45.55 0.56 7
45.80 0.80 9
45.85 0.89 9
45.90 0.78 9
46.15 0.72 13
46.20 0.75 8
46.25 0.69 19
46.60 0.56 13
46.65 0.55 13
47.00 0.53 13
47.05 0.74 8
47.35 0.54 14
47.40 0.67 8
47.75 0.53 12
57.05 0.62 9
57.10 0.63 9
57.40 0.56 13
57.45 0.73 8
57.80 0.57 13
61.65 0.52 13
71.10 0.65 9
71.45 0.60 13
PROME. 0.69 PROME. 0.66 PROME. 0.66 PROME. 0.73 PROME. 0.65 PROME. 0.68
MAXIMO 0.87 MAXIMO 0.85 MAXIMO 0.89 MAXIMO 0.90 MAXIMO 0.87 MAXIMO 0.96
NOTA : La numeracin de los nudos se muesta en la figura 5.3
FACTORES DE DISTRIBUCION DE MOMENTOS MAXIMOS EN CADA INSTANTE
115
5.4.3 SELECCION DEL FACTOR DE DISTRIBUCION DE MOMENTOS
Una mayor limitante de los anlisis tiempo-historia no lineal es la dispersin de sus
resultados debido a su gran sensibilidad a los registros de movimientos ssmicos que se
usan y a los parmetros que modelan la estructura. Tal como se explic en la seccin 4.7,
algunas normas (9) indican qu respuesta estimada tomar de acuerdo al nmero de
registros de sismos utilizados en el anlisis, pudiendo ser la envolvente o el promedio de las
demandas registradas.
En la presente investigacin se han combinado los dos criterios:
Como se mencion anteriormente para cada instante se elige el mayor valor del FDM de
todos los nudos del prtico que entraron en fluencia y es el valor que representa a ese
instante (ver el cuadro 5.1), debe notarse que si en un intervalo de integracin ningn
nudo entr en fluencia, entonces, no se calculan los FDM para ese instante; la cantidad
de valores de FDM depender qu tanto incursiona en el rango inelstico el prtico para
ese sismo.
Una vez obtenida la historia de los FDM durante todo el sismo, se calcula el promedio,
que vendra ser el promedio de los mximos, y es el valor que representa al prtico para
un sismo en particular.
Se estudiaron los FDM de los nudos exteriores, en forma independiente de los nudos
interiores, los resultados del cuadro 5.1 se refieren a los FDM de los nudos interiores.

5.5 ANALISIS E INTERPRETACION DE LOS RESULTADOS
Vistos los resultados de los 19 prticos, se pudo observar que los FDM de los nudos
exteriores eran marcadamente mayores que el de los nudos interiores, por lo que se opt
analizarlos en forma separada. Dicha diferencia se debe a que las cargas axiales inducidas
por el sismo, son especialmente crticas en el caso de columnas exteriores.
En los cuadros 5.2 y 5.3, podemos ver los resultados de los promedios de los valores
mximos de los FDM que se encontraron en los 19 prticos en nudos interiores y exteriores
sometidos a 6 registros de sismos, escalados a 0.4g.

116
SISMOS
PORTICO 10/17/66 05/31/70 10/03/74 PROM. ENVOL.
N82W N08E L T N08E N82W
P-2 0.72 0.68 0.68 0.64 0.61 0.71 0.67 0.72
P-3 0.69 0.63 0.67 0.66 0.66 0.66 0.66 0.69
P-4 0.68 0.66 0.64 0.66 0.65 0.65 0.66 0.68
P-5 0.69 0.66 0.66 0.73 0.65 0.68 0.68 0.73
P-6 0.74 0.67 0.67 0.66 0.74 0.68 0.69 0.74
P-7 0.70 0.74 0.67 0.67 0.74 0.72 0.71 0.74
P-8 0.67 0.69 0.68 0.66 0.63 0.68 0.67 0.69
P-9 0.78 0.86 0.68 0.66 0.76 0.70 0.74 0.86
P-10 0.80 0.64 0.67 0.64 0.78 0.71 0.71 0.80
P-11 0.64 0.66 0.65 0.70 0.70 0.65 0.67 0.70
P-12 0.79 0.70 0.68 0.84 0.66 0.66 0.72 0.84
P-13 0.76 0.65 0.66 0.67 0.61 0.67 0.67 0.76
P-14 0.93 0.72 0.73 0.70 0.71 0.75 0.76 0.93
P-15 0.72 0.69 0.71 0.77 0.73 0.80 0.74 0.80
P-16 0.68 0.67 0.70 0.72 0.69 0.87 0.72 0.87
P-17 0.70 0.70 0.72 0.79 0.74 0.87 0.75 0.87
P-18 0.69 0.93 0.83 0.92 0.85 0.84 0.84 0.93
P-19 0.91 0.91 0.93 1.03 1.01 0.87 0.94 1.03
P-20 0.98 0.86 0.63 1.02 1.02 0.88 0.90 1.02
CUADRO 5.2
FACTOR DE DISTRIBUCION DE MOMENTOS (FDM)
PROMEDIO DE VALORES MAXIMOS DURANTE EL SISMO
NUDOS INTERIORES

SISMOS
PORTICO 10/17/66 05/31/70 10/03/74 PROM. ENVOL.
N82W N08E L T N08E N82W
P-2 1.04 1.01 0.99 0.94 0.88 1.22 1.01 1.22
P-3 0.94 0.81 0.81 0.83 0.86 0.86 0.85 0.94
P-4 0.92 0.73 0.81 0.82 0.79 0.87 0.82 0.92
P-5 0.84 0.80 0.83 1.00 0.74 0.94 0.86 1.00
P-6 0.96 0.81 0.80 0.84 0.80 0.90 0.85 0.96
P-7 0.83 0.91 0.91 0.85 1.04 1.00 0.93 1.04
P-8 0.82 0.83 0.82 0.80 0.75 0.88 0.81 0.88
P-9 0.88 0.97 0.83 0.81 0.90 0.93 0.89 0.97
P-10 0.87 0.80 0.79 0.72 0.84 0.94 0.83 0.94
P-11 0.72 0.79 0.74 0.77 0.77 0.87 0.78 0.87
P-12 0.87 0.86 0.72 0.77 0.80 0.90 0.82 0.90
P-13 0.85 0.92 0.78 0.84 0.77 0.96 0.85 0.96
P-14 0.94 0.84 0.81 0.86 0.83 1.09 0.89 1.09
P-15 0.89 0.94 0.75 1.00 0.82 1.30 0.95 1.30
P-16 0.86 0.93 0.83 0.94 0.78 1.48 0.97 1.48
P-17 0.86 0.92 0.82 1.03 0.90 1.46 1.00 1.46
P-18 0.88 1.17 0.98 0.98 0.87 1.19 1.01 1.19
P-19 0.92 1.34 0.95 1.09 1.01 1.28 1.10 1.34
P-20 1.02 1.49 1.00 1.09 1.02 1.29 1.15 1.49
CUADRO 5.3
FACTOR DE DISTRIBUCION DE MOMENTOS (FDM)
PROMEDIO DE VALORES MAXIMOS DURANTE EL SISMO
NUDOS EXTERIORES


117
Se puede observar para nudos interiores, que los valores del FDM varan entre 0.68 y
1.03, donde el valor mayor corresponde para un prtico relativamente flexible, y el valor ms
pequeo es representativo de un prtico relativamente rgido que respondiera predominante
en su primer modo de vibracin.
Para nudos exteriores los valores del FDM varan entre 0.72 y 1.49 de igual manera el
mayor valor corresponde para un prtico relativamente flexible y el menor valor para un
prtico relativamente rgido. Como se dijo anteriormente la diferencia de los valores entre
los nudos interiores y exteriores, se debe a que en stos ltimos las cargas axiales son ms
crticos en las columnas exteriores.
Como se puede apreciar, el intervalo probable de los valores de los FDM depende de
los patrones de momentos a lo largo de las columnas, en los instantes de la respuesta
dinmica de un prtico; y la distribucin inesperada de los momentos en algunos instantes,
como se indic anteriormente, es causado por la fuerte influencia de los modos superiores
de vibracin. Por esta razn se grafic los valores del cuadro 5.2 y 5.3, en funcin de los
perodos de vibracin de cada prtico, estos grficos se presentan en las figuras 5.4 y 5.5.
Es evidente conforme aumenta el perodo de vibracin, se presentan mayores valores de
FDM en los nudos de los prticos analizados, presentndose valores ms altos para el caso
de nudos exteriores (ver figura 5.5).

118
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.1
1.2
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
PERIODO FUNDAMENTAL (seg)
F
D
M
10/17/66 N82W 10/17/66 N08E 05/31/70 L 05/31/70 T
10/03/74 N08E 10/03/74 N82W PROMEDIO ENVOLV.

Figura 5.4 FACTOR DE DISTRIBUCION DE MOMENTOS EN NUDOS
INTERIORES DE 19 PORTICOS, 6 REGISTROS, 0.4g
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.1
1.2
1.3
1.4
1.5
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
PERIODO FUNDAMENTAL (seg)
F
D
M
10/17/66 N82W 10/17/66 N08E 05/31/70 L 05/31/70 T
10/03/74 N08E 10/03/74 N82W PROMEDIO ENVOLV.

Figura 5.5 FACTOR DE DISTRIBUCION DE MOMENTOS EN NUDOS
EXTERIORES DE 19 PORTICOS, 6 REGISTROS, 0.4g

119

Para ilustrar las implicaciones de estos valores se compara con los requerimientos
actuales de las siguientes normas: Peruana E-060 (27), Americana ACI (3) y Mexicana (25).
De acuerdo con estos cdigos se encuentra que:

Muv 70 . 0 Muv ) 4 . 1 ( 5 . 0 1 Muc = E-060 (5.9)
Muv 60 . 0 Muv ) 2 . 1 ( 5 . 0 1 Muc = ACI (5.10)
Muv 75 . 0 Muv ) 5 . 1 ( 5 . 0 1 Muc = MEXICO (5.11)

El factor de 0.5 en las ecuaciones 5.9, 5.10 y 5.11, se debe a la suposicin de que la
descarga total Mv, del momento de viga est distribuida en proporciones iguales entre las
secciones de columnas arriba y debajo de la planta.
En las figuras 5.6 y 5.7, se comparan los requerimientos de las normas E-060, ACI y
Mxico con los resultados obtenidos de los anlisis dinmicos de los 19 prticos.
Para el caso de nudos interiores, se observa que los valores promedio tienen cierta
concordancia con el valor requerido por la norma E-060 de 0.7 hasta cierto rango de
perodos, sin embargo a partir de 1.35 seg., tiende a incrementarse. Por medio de un
anlisis estadstico se est ajustando una lnea a los valores encontrados, y se est
proponiendo una relacin lineal de dos tramos (ver figura 5.6).

T 52 . 0 FDM = 0.70 FDM 1.05 (5.12)

Donde T, es el perodo fundamental de la estructura.



120
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.1
1.2
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
PERIODO FUNDAMENTAL (seg)
F
D
M
PROM. ENVOLV. E.060 ACI MEXICO PROPUESTA

Figura 5.6 RESULTADOS ANALITICOS EN NUDOS INTERIORES

Para el caso de nudos exteriores, se observa que los valores promedio son mayores
que los requeridos por las normas en comparacin, con una tendencia constante hasta
cierto rango de perodos, sin embargo a partir de 1.35 seg., tiende a incrementarse. De
igual modo por medio de un anlisis estadstico se est ajustando una lnea a los valores
encontrados, y se est proponiendo una relacin lineal de dos tramos (ver figura 5.7).

T 67 . 0 FDM = 0.90 FDM 1.50 (5.13)





121
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.1
1.2
1.3
1.4
1.5
0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0
PERIODO (seg)
F
D
M
PROM. ENVOLV. E.060 ACI MEXICO PROPUESTA

Figura 5.7 RESULTADOS ANALITICOS EN NUDOS EXTERIORES








CAPITULO VI


CONCLUSIONES Y RECOMENDACIONES


6.1 CONCLUSIONES
Las respuestas inelsticas son muy sensibles a los registros ssmicos, se ha visto para
un determinado prtico, sometidos a sismos escalados a un nivel de aceleracin
iguales, las incursiones en el rango inelstico son totalmente distintas.
Los valores de los FDM de los nudos exteriores son mayores que el de los nudos
interiores, debido a que las cargas axiales inducidas por el sismo, son especialmente
crticas en el caso de columnas exteriores.
La contribucin de los modos altos a la repuesta total en general aumenta al aumentar
el perodo fundamental de vibracin de la estructura. Observando las respuestas
locales, la contribucin de los modos ms altos es ms importante en los momentos de
las columnas que en los momentos de las vigas, motivo por el que se presentan
distribuciones inesperadas de momentos flexionantes en las columnas. En tal sentido,
los valores de los FDM en los nudos de la estructura tienden a ser mayores, conforme
123
aumenta el perodo de la estructura, especialmente debido a la fuerte influencia de los
modos superiores de vibracin segundo y tercero.
Con el propsito de reducir la posibilidad de fluencia en las columnas que se consideren
como parte del sistema resistente a fuerzas laterales, se presenta como propuesta, los
siguientes valores de Factores de Distribucin de Momentos (FDM), para ser
incorporados en la Norma Tcnica de Edificacin E-060 Concreto Armado del
Reglamento Nacional de Construcciones, como una disposicin especial para columnas
sujetas a flexocompresin que resistan fuerzas de sismo:

Las resistencias a flexin de las columnas deben satisfacer la siguiente ecuacin

> Mnv FDM Mnc

Donde:
T 52 . 0 FDM = 0.70 FDM 1.05 para nudos interiores
T 67 . 0 FDM = 0.90 FDM 1.50 para nudos exteriores
T = Periodo fundamental de la estructura
Mnc = Momento, en la cara del nudo, correspondiente a la resistencia nominal a
flexin de la seccin superior o inferior (debe verificarse en ambas
secciones) de las columnas que concurren a dicho nudo. La resistencia a la
flexin de la columna debe calcularse para la fuerza axial factorizada,
consistente con la direccin de las fuerzas laterales consideradas, que d
la ms baja resistencia a la flexin.
Mnv = Es la suma de momentos en las caras de nudo, correspondiente a la
resistencia nominal a flexin de las vigas que llegan a dicho nudo. Las
resistencias a la flexin deben sumarse de tal forma que los momentos de
la columna se opongan a los momentos de la viga. Deber satisfacerse la
ecuacin anterior para momentos de viga que acten en ambas direcciones
en el plano vertical del prtico que se considera.

6.2 RECOMENDACIONES

124
Debido a que las cargas ssmicas laterales actan en una direccin general, stas
amplifican los momentos flexionantes en las columnas. Se recomienda completar el
estudio, considerando el efecto de carga biaxial.
La sobrerresistencia por flexin en las rtulas plsticas de las vigas, darn como
resultado momentos flexionantes altos que sern aplicados a las columnas. Se
recomienda investigar la sobrerresistencia, de los esfuerzos de fluencia especificados
para el acero de refuerzo producido en el pas y su influencia en la sobrerresistencia
para aquellos elementos con responsabilidad ssmica.
Analizar otros tipos de suelos. En este estudio se ha empleado suelos tipo S1.
Analizar configuraciones con muros de corte en especial para prticos medianamente
altos y altos.


ANEXO A


MANUAL DE USO DEL PROGRAMA DRAIN-2D
VERSION PARA COMPUTADORAS PERSONALES


UNIVERSITY OF CALIFORNIA Division of Structural Engineering
Berkeley and Structural Mechanics


DRAIN - 2D: Programa de cmputo de propsito general para el anlisis de la respuesta
dinmica de estructuras planas inelsticas bajo movimientos ssmicos.


Programado: G. H. Powell and A. E. Kanaan,
University of California, Berkeley 1973.

Adaptado para Computadoras Personales por:
Manuel Flores Barreto,
PUCP, Lima 1991

126
ADAPTACION DEL PROGRAMA DRAIN-2D PARA COMPUTADORAS PERSONALES
El programa original del DRAIN-2D, fue desarrollado en lenguaje FORTRAN estndar
completo, orientado para usuarios de computadoras grandes (mainframe) como el CDC
6400/6600/7600 (Control Data Cyber) y similares que presentan una alta precisin,
adecuada para resolver problemas de anlisis estructurales en forma exacta, utilizando la
aritmtica de simple precisin. Con la finalidad de adaptar este programa para ser usado
en computadoras personales y compatibles, se tuvo que realizar algunos cambios en el
programa fuente para ser compatible con el estndar Subset FORTRAN 77, y para
conseguir el programa ejecutable se utiliz un compilador de Lahey Computer Systems Inc.
A continuacin se explican algunos detalles de las correcciones necesarias para dicho fin:

La primera dificultad que se encontr fue la conversin de las variables a doble
precisin, ya que el programa utiliza la tcnica de asignacin dinmica de memoria
(dynamic storage allocation) en el tiempo de ejecucin, y que los bloques de informacin
que vienen siendo almacenados y transferidos contienen mezclas de variables enteras y
reales. En primera instancia se pens la conversin a la aritmtica de doble precisin
slo para la matriz de rigidez y el vector de carga, a fin de asegurar la precisin en la
solucin de las ecuaciones de equilibrio, sin embargo resulto ms simple la conversin
declarando en forma implcita todas las variables reales a doble precisin.
Posteriormente se compatibiliz las funciones intrnsecas para manejar la aritmtica
relacionada con los tipos de datos de doble precisin. En general, el nombre especfico
es el mismo que se emple para variables reales, pero con una D como prefijo para las
funciones de doble precisin (7). Por ejemplo, la funcin real SQRT es DSQRT en los
tipos de dato de doble precisin.
Se cambi las proposiciones bsicas de entrada y salida, designando en ellas el nmero
de unidad a emplear, ya que la versin original tena esto en forma de descripcin de
nombre (o sea, no utilizaba la designacin de unidad). Se usaron las siguientes
designaciones de unidad: los nmeros 6 y 7 para los archivos de entrada de datos y
salida de resultados ambos en tipo texto, los nmeros 1, 2, 3, 4, 8 y 9 para archivos
temporales binarios, finalmente el nmero 10 para un archivo donde se almacenarn

127
todos los resultados del anlisis en forma binaria; este ltimo archivo que aparecer con
extensin PLT puede ser utilizado para aplicaciones de post procesamiento, como
graficar resultados tiempo-historia, etc.

PREPARACION DEL ARCHIVO DE DATOS
El archivo de datos debe ser preparado con un editor de textos, todos los datos tienen
un formato preestablecido y deben estar de acuerdo con el presente manual, los nombres
de los archivos de datos no deben tener extensiones y deben estar en el mismo directorio
donde se encuentra el ejecutable. Al momento de ejecutarse el programa, emitir una serie
de mensajes en la pantalla para identificar el progreso de los distintos pasos de ingreso de
datos y solucin, y crear un archivo de salida con el nombre del archivo de datos ms una
extensin OUT. Este archivo tipo texto, contendr un eco de los datos ingresados y los
resultados de los anlisis indicados.

Con el propsito de facilitar la elaboracin del archivo de datos, se ha introducido al
programa principal una subrutina para que reconozca lneas de datos para comentarios.
Para tal fin cualquier lnea que tenga una letra C en la primera columna del archivo de datos
ser tratada como una lnea de comentarios y ser ignorada por el programa. En el Anexo
B, se presentan algunos ejemplos de archivos de datos.

Durante la ejecucin del programa DRAIN, se crean una serie de archivos secundarios,
estos archivos no son necesarios a menos que exista una aplicacin para su post
procesamiento, para borrar esta serie de archivos, se adjunta un comando tipo batch
llamado D. Para ejecutarlo desde donde se encuentran los archivos residentes, entrar el
comando D seguido del nombre del archivo de datos y borrar todos los archivos
secundarios con excepcin de los archivos de datos y de salida de resultados.

La presente versin cuenta con la siguiente biblioteca de elementos:
1. Elemento armadura
2. Elemento viga-columna

128
3. Elemento muro tabique
4. Elemento de conexin semi-rgido
5. Elemento viga
6. Elemento viga con rigidez degradante

OBJETIVO DEL PROGRAMA
El programa determina la respuesta dinmica de estructuras bidimensionales inelsticas
de configuracin arbitraria que resultan de los movimientos de la base debido a sismos.
Pueden especificarse sismos independientes tanto vertical y horizontal, pero los
movimientos de apoyo fuera de fase no pueden ser considerados. Pueden aplicarse cargas
estticas a la estructura antes de la aplicacin de la carga dinmica, pero, el
comportamiento bajo las cargas estticas se asume elstico.
La estructura puede estar compuesta de una variedad de tipos de elementos, cada uno
de stos presenta un patrn de comportamiento y caractersticas de fluencia distinto. El
programa est diseado para permitir que se aadan nuevos elementos con una cantidad
relativamente pequea de esfuerzo de codificacin.

IDEALIZACION
La estructura es idealizada como un conjunto de nudos (uniones) conectados por
elementos (miembros). Los nudos deben ser numerados secuencialmente empezando con
uno.
Cada nudo tiene tres grados de libertad, llmense desplazamiento horizontal, vertical y
rotacional. Sin embargo, algunos desplazamientos pueden ser especificados como cero, y
desplazamientos en diferentes nudos pueden ser especificados como idnticos. El uso de
estas opciones puede reducir en forma sustancial el nmero de ecuaciones que se deben
resolver y, en consecuencia, reducir el tiempo de solucin.
Para minimizar los requisitos del almacenamiento y el tiempo de cmputo, los nudos
debern ser numerados de manera que la diferencia mxima entre los nmeros de los
grados de libertad de la estructura conectados para cualquier miembro sea tan pequeo
como sea posible (es decir, se reduce el ancho de banda de la matriz de rigidez).

129
Para la entrada y salida, los elementos deben ser divididos en grupos. Todos los
elementos de cualquier grupo deben ser del mismo tipo (por ejemplo, ARMADURA, VIGA-
COLUMNA, etc.) y, comnmente todos los elementos de un solo tipo se incluirn en un solo
grupo. Sin embargo, elementos de un mismo tipo pueden ser subdivididos en grupos
separados si se desea. Dentro de cualquier grupo, los elementos deben ser numerados
secuencialmente, empezando por uno.
El comportamiento de cualquier elemento es determinado por el programador que
desarrolla las subrutinas para el elemento. Este comportamiento puede ser implcitamente
de cualquier grado de complejidad dependiendo del propsito para el cual se propone el
elemento. Las caractersticas de los elementos disponibles para sta versin se describen
en los apndices B de la referencia 23.
Los desplazamientos y las coordenadas nodales se definen con respecto a un sistema
de coordenadas global de mano derecha, con los ejes X e Y en el plano de la estructura y el
eje Z encima, fuera de ste plano. Comnmente, el eje X ser horizontal hacia lado
derecho y el eje Y estar verticalmente hacia arriba.

LIMITACIONES DE LA CAPACIDAD
Algunas limitaciones de la capacidad que resultan del uso de dimensiones fijas, se
asocian con el ingreso de la informacin del elemento, tal como se explica en la Seccin E
de sta gua. Estas limitaciones podran ser superadas por el cambio de un nmero de
declaraciones de dimensin. Sin embargo, la mayor parte de la informacin es almacenada
en arreglos con dimensiones variables, el mayor almacenamiento de memoria interna
requerido se establece en la proposicin COMMON blanco y se asigna dinmicamente en el
tiempo de ejecucin del programa.
La longitud mnima requerida del COMMON blanco depende de la cantidad de nudos en
la estructura, la forma de la matriz de rigidez de la estructura, el nmero de masas
concentradas, el nmero de pasos de integracin y la cantidad de nudos para los que se
requieren las historias de los desplazamientos en el tiempo. Ya que la rigidez de la
estructura se almacena en una forma compacta y que sta rigidez comnmente ocupa la
mayor parte del almacenamiento mnimo requerido, es difcil calcular ste almacenamiento

130
manualmente. Por lo tanto, se calcula el almacenamiento mediante el programa y se
imprime con el archivo de salida. Se presume que una o ms ejecuciones de verificacin
de la informacin se requerirn antes que se lleve a cabo cualquier cmputo extensivo.
Para una ejecucin inicial, un clculo estimado del almacenamiento mnimo requerido
puede ser realizado como sigue:
NJTS = nmero de nudos en la estructura
NEQ = nmero de grados de libertad de la estructura (conservadoramente NJTS x 3)
MBAND = mximo ancho de la semi banda de la matriz de rigidez de la estructura.
NMASS = nmero de nudos con masas concentradas.
NSTEPS = nmero de pasos de integracin
NELTOT = nmero total de elementos en la estructura
NHOUT, NVOUT, NROUT = nmero de nudos para los que se requieren las
impresiones de la historia en el tiempo de los desplazamientos horizontales, verticales y
rotacionales, respectivamente.
KDUP = 2 si la matriz de rigidez de estructura debe ser retenida en la memoria y es = 1
si esta matriz debe ser retenida en archivos secundarios, (Ver seccin B1 de esta gua
para mayor explicacin).
Entonces el almacenamiento requerido no es mayor de:

2xNSTEPS+NHOUT+NVOUT+NROUT+12xNEQ+NEQxKDUPxMBAND +NELTOT

Adems de este almacenamiento mnimo, se debe proporcionar un espacio adicional en
el COMMON blanco para almacenar la informacin de los elementos. La cantidad de
posiciones de almacenamiento requerido para cada tipo diferente de elemento se observa
en la seccin E de esta gua. La informacin de los elementos se compacta en la memoria
interna disponible y si se proporciona una longitud suficiente del COMMON blanco, la
informacin de todos los elementos se retendr permanentemente en la memoria interna,
en consecuencia el tiempo de cmputo de ingreso y salida de la ejecucin del programa
ser reducido. Si se proporciona una longitud insuficiente, la informacin del elemento se
dividir automticamente en bloques, con tantos elementos como sea posible en cada

131
bloque. Los bloques sern grabados en archivos de almacenamiento exterior (archivos
binarios secundarios) y, se incrementar una cantidad considerable del tiempo de cmputo
del programa.
En la salida se imprime el nmero de bloques de informacin del elemento y la longitud
COMMON blanco requerida para todos los elementos que ser almacenada en un bloque
nico.
El almacenamiento asignado puede ser modificado por el cambio de las dos
proposiciones siguientes en el programa principal:

COMMON = A(NT)
NTST = NT

En la que NT = longitud total del COMMON blanco.
La presente versin que se describe en esta gua, ha sido compilada con NT=40000 (la
mxima que aceptaba el compilador Lahey que se utiliz) de longitud total de la proposicin
COMMON.

132
INGRESO DE DATOS

Las siguientes instrucciones definen el problema que se va a resolver. Las unidades a
emplear, deben ser consistentes.

A. INICIO Y TITULO DEL PROBLEMA (A5, 3X, 18A4)
UNA LINEA DE DATOS
Columnas
01-05: Escribir la palabra START
06-08: Blanco
09-80: Ttulo del problema para ser impreso con la salida.

B. INFORMACION GEOMETRICA DE LA ESTRUCTURA
B.1 INFORMACION DE CONTROL (9I5, I10) UNA LINEA DE DATOS
Columnas
01-05: Nmero de nudos en la estructura
06-10: Nmero de nudos control para los cuales se especifican las
coordenadas directamente. Ver seccin B2.
11-15: Nmero de comandos para la generacin de coordenadas de nudos.
Ver seccin B3.
16-20: Nmero de comandos para especificar los nudos con
desplazamientos cero. Ver seccin B4.
21-25: Nmero de comandos para especificar los nudos con
desplazamientos idnticos. Ver seccin B5.
26-30: Nmero de comandos para especificar las masas concentradas en
los nudos. Ver seccin B6.
31-35: Nmero de grupos diferentes de elementos en la estructura. Ver
seccin E. Mximo 10.
36-40: Cdigo para verificacin de datos. Poner 1 si slo se requiere correr
el verificador de datos. Dejar en blanco o poner cero si el problema

133
debe ser ejecutado. Poner -1 y el problema ser ejecutado siempre y
cuando el nmero de bloques de informacin del elemento no
exceda a 1. Esta ltima opcin previene la ejecucin con tiempo de
cmputo excesivo de entrada y salida, en casos donde se pens
realizar la operacin en la memoria interna, pero surgen errores al
especificar el almacenamiento requerido.
41-45: Cdigo de almacenamiento de rigidez de la estructura. La matriz de
rigidez de la estructura duplicada debe ser retenida y peridicamente
actualizada. Deje un espacio en blanco o colocar cero si esta matriz
debe ser retenida en la memoria interna, de manera que se reduzca
el tiempo de cmputo. Poner 1 si la matriz debe ser grabada en
archivos temporales de almacenamiento externo.
46-55: Longitud del COMMON blanco, a ser considerado, si es cero o est
en blanco, se utilizar el valor compilado en el programa. Ver la
discusin de las limitaciones de la capacidad.

B2. COORDENADAS DE LOS NUDOS DE CONTROL (I5, 2F10.0)
UNA LINEA DE DATOS PARA CADA NUDO DE CONTROL
Ver OBSERVACION 1 para mayor explicacin
Columnas
01-05: Nmero del nudo, en cualquier secuencia.
06-15: Coordenada X del nudo
16-25: Coordenada Y del nudo

B3. COMANDOS PARA LA GENERACION LINEAL DE COORDENADAS DE NUDOS
(4I5, F10.0)
UNA LINEA DE DATOS PARA CADA COMANDO DE GENERACION

Omitir si no existen comandos de generacin. Ver la OBSERVACION 1 para mayor
explicacin.

134
Columnas
01-05: Nmero del nudo al inicio de la lnea de generacin.
06-10: Nmero del nudo al final de la lnea de generacin
11-15: Nmero de nudos que sern generados a lo largo de la lnea.
16-20: Diferencia del nmero del nudo (constante) entre dos nudos
sucesivos en la lnea. Si est en blanco o cero, se presume que es
igual a 1.
21-30: Espacio entre los nudos sucesivos generados. Si es mayor o igual a
1.0, se presume que es el espaciado actual. Si es menor de 1.0, se
presume que es el espaciado actual dividido por la longitud de la
lnea de generacin. Si es cero o est en blanco, los nudos son
automticamente espaciados en forma uniforme a lo largo de la lnea
de generacin.

B4. COMANDOS PARA NUDOS CON DESPLAZAMIENTOS CERO (6I5)
UNA LINEA DE DATOS PARA CADA COMANDO
Omitir, si ningn nudo se limita a tener desplazamientos cero. Ver la OBSERVACION 2
para mayor explicacin.
Columnas
01-05: Nmero del nudo o nmero del primer nudo en una serie de nudos
comprendidos por este comando.
06-10: Cdigo para el desplazamiento X. Poner 1 si el desplazamiento es
nulo, de otra forma deje en blanco o cero.
11-15: Cdigo para el desplazamiento Y (uno o cero).
16-20: Cdigo de rotacin (uno o cero).
21-25: Nmero del ltimo nudo en la serie. Deje en blanco para un nico
nudo.
26-30: Diferencia del nmero del nudo (constante) entre nudos sucesivos en
la serie. Si es blanco o cero se presume que es igual a 1.


135
B5. COMANDOS PARA NUDOS CON DESPLAZAMIENTOS IGUALES (16I5)
UNA LINEA DE DATOS PARA CADA COMANDO
Omitir si ninguno de los nudos se limita a tener desplazamientos iguales. Ver
OBSERVACION 3 para mayor explicacin.
Columnas
01-05: Cdigo de desplazamiento, como sigue:
Poner 1 para el desplazamiento X
Poner 2 para el desplazamiento Y
Poner 3 para la rotacin
06-10: Nmero de nudos comprendidos por este comando (mx. 14- ver
OBSERVACION 3 para el procedimiento, cuando ms de 14 nudos
tienen desplazamientos idnticos).
11-80: Lista de nudos, en un orden numrico creciente, hasta 14 nudos en
campos de 5 columnas.

B6. COMANDOS PARA MASAS CONCENTRADAS EN LOS NUDOS (I5, 3F10.0, 2I5,
F10.0)
UNA LINEA DE DATOS PARA CADA COMANDO
Ver OBSERVACION 4 para mayor explicacin
Columnas
01-05: Nmero del nudo, o nmero del primer nudo en una serie de nudos
comprendidos por este comando.
06-15: Masa asociada con el desplazamiento X. Puede ser cero.
16-25: Masa asociada con el desplazamiento Y. Puede ser cero.
26-35: Inercia rotacional. Puede ser cero.
36-40: Nmero del ltimo nudo en la serie. Dejar en blanco para un solo
nudo.
41-45: Diferencia del nmero del nudo entre nudos sucesivos de la serie. Si
est en blanco o cero se presume que es igual a 1.

136
46-55: Factor de modificacin mediante el cual las masas son divididas. Si
est en blanco o cero, se presume el factor del comando precedente,
de manera que si se aplica el mismo factor a todos los comandos,
ste debe ser especificado slo en el primer comando. Este factor
es comnmente igual a g (aceleracin de la gravedad) en tales casos
los valores de masa estarn dados como pesos.

C. INFORMACION DE CARGA
C 1. INFORMACION DEL CONTROL DE CARGA (3I5, 6F10.0)
UNA LINEA DE DATOS
Columnas
01-05: Cdigo de carga esttica. Poner 1 si las cargas estticas deben
aplicarse antes de las cargas dinmicas. Deje en blanco o cero si no
existen cargas estticas.
06-10: Nmero de comandos para especificar cargas estticas, aplicadas
directamente en los nudos. Ver la seccin C2. Deje en blanco o
cero si no existen cargas estticas aplicadas directamente en los
nudos.
11-15: Nmero de pasos de integracin que deben ser considerados en el
anlisis dinmico.
16-25: Tamao del paso de integracin
T
.
26-35: Factor de amplificacin que debe aplicarse a las aceleraciones del
suelo especificadas para la direccin X. Ver la OBSERVACION 5
para mayor explicacin.
36-45: Factor de amplificacin que debe aplicarse a la escala de tiempo del
registro de aceleracin especificado para la direccin X. Ver
OBSERVACION 5. Si est en blanco o cero se presume que es
igual a 1.
46-55: Factor de amplificacin para las aceleraciones del suelo en la
direccin Y.

137
56-65: Factor de amplificacin para la escala de tiempo en la direccin Y.
Si est en blanco o cero se presume que es igual a 1.
66-75: Valor absoluto del mximo desplazamiento permitido antes que se
pueda presumir que la estructura colapse. La ejecucin finaliza si
ste valor es excedido en cualquier paso. Si es cero o est en
blanco, tomar el valor por defecto (100,000).

C2. COMANDOS PARA CARGAS ESTATICAS APLICADAS DIRECTAMENTE EN
LOS NUDOS (I5, 3F10.0, 2I5)
UNA LINEA DE DATOS PARA CADA COMANDO.
Omitir, si no existen cargas estticas aplicadas directamente en los nudos. Si el cdigo
de la carga esttica se encuentra en cero o en blanco (lnea C1), las cargas sern ledas e
impresas, pero, se ignorarn en el anlisis de la estructura.
Columnas
01-05: Nmero del nudo, o el nmero del primer nudo en una serie de nudos
comprendidos por este comando.
06-15: Carga en la direccin X, la misma en todos los nudos de la serie.
16-25: Carga en la direccin Y, la misma en todos los nudos de la serie.
Ntese que la direccin Y ser por lo general hacia arriba positiva.
26-35: Momento (regla del tornillo de la mano derecha respecto al eje Z
considerando positivo el sentido antihorario).
36-40: Nmero del ltimo nudo en la serie. Deje en blanco para un slo
nudo.
41-45: Diferencia del nmero del nudo (constante) entre nudos sucesivos en
la serie. Si est en blanco o en cero, se asume que es igual a 1.

NOTA: Un nudo puede aparecer en dos o ms comandos si se desea. En tal caso, las
cargas totales aplicadas en el nudo sern la suma de las cargas de los comandos
separados.


138
C3. REGISTROS DE ACELERACION.
C3(a) INFORMACION DEL CONTROL (4I5, 10A6)
UNA LINEA DE DATOS
Columnas
01-05: Nmero de pares de tiempo-aceleracin que definen el movimiento
del terreno en la direccin X (NPTH). Poner cero o dejar en blanco
cuando no haya movimiento en esta direccin.
06-10: Nmero de pares de tiempo-aceleracin que definen el movimiento
del terreno en la direccin Y (NPTV). Poner cero o dejar en blanco
cuando no haya movimiento en esta direccin.
11-15: Cdigo para imprimir en el archivo de salida las aceleraciones de
entrada. Dejar en blanco o poner cero para no imprimir. Poner 1 a
fin de obtener la lista de registro de aceleracin.
16-20: Cdigo para la impresin en el archivo de salida de la aceleracin
interpolada en intervalos
T
. Dejar en blanco o poner cero para no
imprimir. Poner 1 para obtener la lista de registro de aceleracin.
21-80: Ttulo para identificar el registro del sismo, sern impresos en la
salida.

C3(b) HISTORIA EN EL TIEMPO DE LA ACELERACION DEL TERRENO EN LA
DIRECCION X (10F8.0)
Tantas lneas de datos como se necesite para especificar los NPTH pares de tiempo-
aceleracin, 5 pares por lnea, debe estar en unidades de aceleracin (no para datos
mltiplos de la aceleracin debido a la gravedad). Omitir si NPTH es igual a cero. Ntese
que tanto la aceleracin, como las escalas del tiempo, pueden ampliarse si se desea (ver
seccin C1). Si el registro entra en trminos de la aceleracin debido a la gravedad, un
factor de amplificacin de aceleracin igual a g puede utilizarse para convertirlo a unidades
de aceleracin.
El primer par tiempo-aceleracin debe ser cero, tanto para el tiempo como para la
aceleracin.

139

C3(c) HISTORIA EN EL TIEMPO DE LA ACELERACION DEL TERRENO EN LA
DIRECCION Y (10F8.0)
Tantas lneas de datos como se necesiten para especificar los NPTV pares de tiempo-
aceleracin. Omitir si NPTV es igual cero.

C4. INFORMACION DEL AMORTIGUAMIENTO (4F10.0)
UNA LINEA DE DATOS.
Ver OBSERVACION 6 para mayor explicacin.
Columnas
01-10: Factor de amortiguamiento proporcional a la masa,
11-20: Factor de amortiguamiento proporcional a la rigidez, , para la rigidez
tangente en un determinado instante.
21-30: Factor de amortiguamiento proporcional a la rigidez, o, para la
rigidez elstica original.
31-40: Factor de amortiguamiento Estructural, .

D. ESPECIFICACION DE SALIDAS TIEMPO-HISTORIA
Si se desea pueden obtenerse las historias en el tiempo impresas de los
desplazamientos nodales seleccionados y los resultados de los elementos en intervalos de
tiempo seleccionados. Los valores envolventes de todos los desplazamientos nodales y los
resultados de los elementos se imprimen automticamente al final del cmputo, como
tambin, si el desplazamiento mximo especificado ha sido excedido. Las envolventes de
los resultados intermedios pueden imprimirse en intervalos de tiempo seleccionados.
Los valores de historias en el tiempo pueden imprimirse: Progresivamente con el
cmputo, solamente al final del cmputo, o en ambas formas. La impresin al final del
cmputo es ordenada elemento por elemento y nudo por nudo, antes que a paso por paso,
y en consecuencia, son ms compactos y ms convenientes para la utilizacin de trazados
de grficos de historias en el tiempo. Estos resultados reordenados pueden grabarse
tambin, en archivos temporales sin formato para su posterior procesamiento.

140

D1. INFORMACION DE CONTROL (13I5) - UNA LINEA DE DATOS
Columnas
01-05: Intervalo de tiempo para la salida tiempo-historia del desplazamiento
nodal, expresado como un mltiplo del paso de integracin
T
.
Dejar en blanco cuando no se requiera la salida. Los nudos para las
cuales se necesitan historias en el tiempo, se especifican en las
secciones D2 al D6.
06-10: Intervalo de tiempo para la salida tiempo-historia de los resultados
del elemento, expresado como un mltiplo del paso de integracin
. Dejar en blanco cuando no se requiera la salida. Los
elementos para las cuales se necesitan las historias en el tiempo, se
especifican en la seccin E.
T
11-15: Intervalos de tiempo para la salida de valores envolventes
intermedios, expresada como un mltiplo del paso de integracin
. Dejar en blanco cuando no se requiera salidas intermedias.
Los valores envolventes se imprimen automticamente al final del
periodo de respuesta.
T
16-20: Nmero de nudos (NHOUT) para las cuales se necesitan las
historias en el tiempo del desplazamiento, en la direccin X.
21-25: Nmero de nudos (NVOUT) para las cuales se necesitan las
historias en el tiempo del desplazamiento, en la direccin Y.
26-30: Nmero de nudos (NROUT) para las cuales se necesitan las
historias en el tiempo de la rotacin.
31-35: Nmero de pares de nudos (NHR) para las cuales se necesitan las
historias en el tiempo del desplazamiento relativo, en la direccin X.
36-40: Nmero de pares de nudos (NVR) para las cuales se necesitan las
historias en el tiempo del desplazamiento relativo, en la direccin Y.
41-45: Cdigo de salida para historias en el tiempo del desplazamiento
nodal, como sigue:

141
a) Cero o blanco: salida tiempo-historia en forma progresiva durante
el cmputo.
b) 1: En ambas formas, salida en forma progresiva y una salida
reordenada al final del cmputo.
c) 2: Solamente una salida reordenada al final del cmputo.
46-50: Cdigo de salidas para historias en el tiempo del desplazamiento
nodal relativo. Cero, 1 o 2 como para el cdigo anterior.
51-55: Cdigo de salida para historias en el tiempo de resultados de los
elementos. Cero, 1 o 2 como para el cdigo anterior.
56-60: Cdigo de almacenamiento en un archivo secundario temporal para
grabar los resultados tiempo-historia reordenadas, de los
desplazamientos nodales y desplazamientos relativos, como sigue:
a) Cero o blanco: Los resultados tiempo-historia reordenadas se
imprimen en el archivo de salida, pero no se graban en un
archivo secundario binario.
b) Distinto de cero: Los resultados tiempo-historia reordenadas se
imprimen en el archivo de salida y tambin, se graban en un
archivo secundario binario. Ver la siguiente NOTA para mayor
explicacin.
61-65: Cdigo de almacenamiento en un archivo secundario temporal para
grabar los resultados tiempo-historia reordenadas de los elementos.
Cero o distinto de cero como el cdigo anterior.

NOTA: Si el cdigo de almacenamiento es distinto de cero, cada lnea impresa de datos
tiempo-historia reordenada se grabarn, en un archivo secundario NO
FORMATEADO como un registro por lnea (es decir, cualquier encabezado se
excluye de la salida); en la presente adaptacin dicho archivo tendr una extensin
PLT. Por consiguiente los resultados tiempo-historia estarn directamente
reproducidas en este archivo binario.

142
Se utilizar esta opcin si se van a desarrollar subrutinas para post procesamiento
de los resultados, como por ejemplo graficar los datos tiempo-historia. El programa
no escribe ningn dato de identificacin del problema en el archivo binario, tampoco
ubica explcitamente el dispositivo de fin de archivo despus del ltimo registro.

D2. LISTA DE NUDOS PARA HISTORIAS EN EL TIEMPO DEL DESPLAZAMIENTO
EN LA DIRECCION X (10I5)
Especificar tantas lneas de datos como se necesiten para los NHOUT nmeros de
nudos, colocar diez por cada lnea, en columnas de 5 espacios. Omitir si el NHOUT es
igual a cero.

D3. LISTA DE NUDOS PARA HISTORIA EN EL TIEMPO DEL DESPLAZAMIENTO EN
LA DIRECCION Y (10I5)
Especificar tantas lneas de datos como se necesiten para los NVOUT nmeros de
nudos, colocar diez por cada lnea, en columnas de 5 espacios. Omitir si el NVOUT es igual
a cero.

D4. LISTA DE NUDOS PARA HISTORIA EN EL TIEMPO DE LA ROTACION (10I5)
Especificar tantas lneas de datos como se necesiten para los NROUT nmeros de
nudos, colocar diez por cada lnea, en columnas de 5 espacios. Omitir si el NROUT es
igual a cero.

D5. LISTA DE NUDOS PARA HISTORIA EN EL TIEMPO DEL DESPLAZAMIENTO
RELATIVO EN LA DIRECCION X (10I5).
Especificar tantas lneas de datos como se necesiten para los NHR pares de nmeros
de nudos, colocar cinco pares por cada lnea, en columnas de 5 espacios. Omitir si el NHR
es igual a cero. El desplazamiento impreso, es el desplazamiento del primer nudo menos el
desplazamiento del segundo nudo de cualquier par.


143
D6. LISTA DE NUDOS PARA HISTORIA EN EL TIEMPO DEL DESPLAZAMIENTO
RELATIVO EN LA DIRECCION Y (10I5).
Especificar tantas lneas de datos como se necesiten para los NVR pares de nmeros
de nudos, colocar cinco pares por cada lnea, en columnas de 5 espacios. Omitir si el NVR
es igual a cero.

E. ESPECIFICACION DEL ELEMENTO
Para la entrada y salida, los elementos deben dividirse en grupos. Todos los elementos
de cualquier grupo deben ser del mismo tipo, y generalmente todos los elementos de un
determinado tipo sern incluidos en un slo grupo. Sin embargo, los elementos de un
mismo tipo pueden subdividirse en grupos separados si se desea.
Los grupos de elementos pueden ser ingresados en cualquier secuencia conveniente.
Dentro de cualquier grupo, los elementos deben ser enumerados en secuencia
comenzando con uno.
El nmero total de grupos de elementos no pueden exceder de diez.

144
E1. ELEMENTO DE ARMADURA
Ver el apndice B1 de la referencia 24, para la descripcin del elemento. Nmero de
posiciones de informacin por elemento = 35.

E1(a) INFORMACION DE CONTROL PARA EL GRUPO (4I5) - UNA LINEA DE DATOS
Columnas
01-05: Poner 1 (para indicar que este grupo consta de ELEMENTOS
ARMADURA).
06-10: Nmero de elementos en el grupo.
11-15: Nmero de diferentes tipos de rigidez del elemento (mximo 40). Ver
Seccin E1 (b).
16-20: Nmero de diferentes patrones de fuerzas de fijacin en los
extremos (mximo 40). Ver Seccin E1 (c).

E1(b) TIPOS DE RIGIDEZ (I5,5F10.0,I5)
UNA LINEA DE DATOS PARA CADA TIPO DE RIGIDEZ.
Columnas
01-05: Nmero del tipo de rigidez, en secuencia comenzando con 1.
06-15: Mdulo de Young de elasticidad.
16-25: Mdulo de endurecimiento por deformacin, como una proporcin
del mdulo de Young.
26-35: Area de la seccin transversal.
36-45: Esfuerzo de fluencia en tensin.
46-55: Esfuerzo de fluencia o esfuerzo de pandeo elstico en compresin.
56-60: Cdigo de pandeo. Poner 1, si el elemento se pandea elsticamente
en compresin. Poner cero o dejar en blanco, si el elemento fluye en
compresin sin pandeo.


145
E1(C) PATRONES DE FUERZAS DE FIJACION EN LOS EXTREMOS (2I5, 4F10.0)
UNA LINEA DE DATOS PARA CADA PATRON.
Omitir si no hay fuerzas de fijacin. Ver figura B1.5.
Columnas
01-05: Nmero del patrn, en secuencia comenzando con 1.
06-10: Cdigo del sistema de ejes, como sigue:
Cdigo = 0: Las fuerzas estn en el sistema de coordenadas del
elemento, como en la figura B1.5a.
Cdigo = 1: Las fuerzas estn en el sistema de coordenadas
global, como en la figura B1.5b.
11-20: Fuerza de fijacin Fi.
21-30: Fuerza de fijacin Vi.
31-40: Fuerza de fijacin Fj.
41-50: Fuerza de fijacin Vj.







(b) CODIGO = 1 (a) CODIGO = 0
Figura B1.5 FUERZAS DE FIJACION EN LOS
EXTREMOS E INICIALES
Vi
Vi
Vj Vj
Fi
Fj
Fj
Fi
E1(d) COMANDOS PARA GENERACION DEL ELEMENTO (9I5, 2F5.0, F10.0)

146
UNA LINEA DE DATOS PARA CADA COMANDO DE GENERACION.
Los elementos deben especificarse en un orden numrico creciente. Deben incluirse
obligatoriamente los datos del primer y ltimo elemento. Ver la OBSERVACION 7 para
mayor explicacin del procedimiento de generacin.
Columnas
01-05: Nmero del elemento o nmero del primer elemento en una serie
numerada secuencialmente de elementos para ser generados por
ste comando.
06-10: Nmero del nudo del elemento en el extremo i.
11-15: Nmero del nudo del elemento en el extremo j.
16-20: Incremento del nmero de nudo para la generacin del elemento. Si
es cero o est en blanco, se asume igual a 1.
21-25: Nmero de tipo de rigidez.
26-30: Cdigo para incluir rigidez geomtrica. Poner 1 si la rigidez
geomtrica debe incluirse. Dejar en blanco o poner cero, si la rigidez
geomtrica debe ignorarse.
31-35: Cdigo de salida de la historia en el tiempo. Si no se requiere los
resultados tiempo-historia del elemento para los elementos
comprendidos por este comando, poner cero o dejar en blanco. Si
se requiere una salida de la historia en el tiempo en los intervalos
especificados en la lnea de datos D1, poner 1.
36-40: Nmero del patrn de fuerzas de fijacin para cargas muertas
estticas en el elemento. Dejar en blanco si no existen cargas
muertas. Ver la NOTA a continuacin.
41-45: Nmero del patrn de fuerzas de fijacin para cargas vivas estticas
en el elemento. Dejar en blanco si no hay cargas vivas.
46-50: Factor de escala que debe ser aplicado a las fuerzas de fijacin
debido a cargas muertas estticas. Dejar en blanco si no hubiera
cargas muertas.

147
51-55: Factor de escala que debe ser aplicado a las fuerzas de fijacin
debido a cargas vivas estticas. Dejar en blanco si no hubiera
cargas vivas.
56-65: Fuerza axial inicial en el elemento, en tensin se considera positiva.

NOTA: Si el cdigo de carga esttica, lnea de dato C1, es cero, pero fuerzas de fijacin
estn especificadas an para algunos elementos, entonces origina una
incongruencia. En efecto, cualquiera de estas fuerzas de fijacin ser tratada como
fuerzas iniciales en el elemento.

148
E2. ELEMENTO VIGA-COLUMNA
Ver el Apndice B2 de la referencia 24, para la descripcin de elemento. Nmero de
posiciones de informacin por elemento = 141.

E2(a) INFORMACION DE CONTROL PARA EL GRUPO (7I5) - UNA LINEA DE DATOS.
Columnas
01-05: Poner 2 (para indicar que el grupo consta de ELEMENTOS VIGA-
COLUMNA).
06-10: Nmero de elementos en el grupo.
11-15: Nmero de diferentes tipos de rigidez del elemento (mximo 40). Ver
seccin E2 (b).
16-20: Nmero de diferentes tipos de excentricidades en los extremos
(mximo 15). Ver seccin E2(c).
21-25: Nmero de diferentes superficies de interaccin de fluencia para las
secciones transversales (mximo 40). Ver seccin E2(d).
26-30: Nmero de diferentes patrones de fuerzas de fijacin en los
extremos (mximo 35). Ver seccin E2 (e).
31-35: Nmero de diferentes patrones de fuerzas iniciales en el elemento
(mximo 30). Ver seccin E2 (f).

E2(b) TIPOS DE RIGIDEZ (I5, 4F10.0, 3F5.0, 2F10.0)
UNA LINEA DE DATOS PARA CADA TIPO DE RIGIDEZ.
Columnas
01-05: Nmero de tipo de rigidez, en secuencia comenzando con 1.
06-15: Mdulo de Young de elasticidad.
16-25: Mdulo de endurecimiento por deformacin, como una proporcin
del mdulo de Young.
26-35: Area de la seccin transversal.
36-45: Momento de inercia.
46-50: Factor de rigidez flexional kii.

149
51-55: Factor de rigidez flexional kjj.
56-60: Factor de rigidez flexional kij.
61-70: Area de corte efectivo. Dejar en blanco o poner cero si las
deformaciones por corte deben ignorarse, o si stas ya han sido
consideradas en el clculo de los factores de rigidez flexional.
71-80: Mdulo de Poisson (utilizado para calcular el mdulo de corte, y
requerido slo s las deformaciones por corte deben considerarse).

E2(c) EXCENTRICIDAD EN LOS EXTREMOS (I5, 4F10.0)
UNA LINEA DE DATOS PARA CADA TIPO DE EXCENTRICIDAD.
Omitir si no hay excentricidades en los extremos. Ver la figura 3.11 para mayor
explicacin. Todas las excentricidades se miden desde el nudo hasta el final del elemento.
Columnas
01-05: Nmero del tipo de excentricidad en los extremos, en secuencia
comenzando con 1.
06-15: Xi = excentricidad X, en el extremo i.
16-25: Xj = excentricidad X, en el extremo j.
26-35: Yi = excentricidad Y, en el extremo i.
36-45: Yj = excentricidad Y, en el extremo j.

E2(d) SUPERFICIES DE INTERACCION DE FLUENCIA DE LA SECCION
TRANSVERSAL (2I5,4F10.0,4F5.0)
UNA LINEA DE DATOS PARA CADA SUPERFICIE DE FLUENCIA.
Ver la figura 3.7 para mayor explicacin.
Columnas
01-05: Nmero de superficie de fluencia, en secuencia comenzando con 1.
06-10: Cdigo de forma de la superficie de fluencia, como sigue:
Poner 1: Tipo viga, sin interaccin P-M.
Poner 2: Tipo columna de acero
Poner 3: Tipo columna de concreto armado.

150
11-20: Momento de fluencia positivo (fibra inferior en tensin) My+.
21-30: Momento de fluencia negativo (fibra superior en tensin) My-.
31-40: Fuerza de fluencia de compresin, Pyc. Dejar en blanco si el cdigo
de forma = 1.
41-50: Fuerza de fluencia de tensin, Pyt. Dejar en blanco si el cdigo de
forma = 1.
51-55: Coordenada-M del punto de balance A, como una proporcin de
My+. Dejar en blanco si el cdigo de forma = 1.
56-60: Coordenada-P del punto de balance A, como una proporcin de Pyc.
Dejar en blanco si el cdigo de forma =1.
61-65: Coordenada-M del punto de balance B, como una proporcin de My-.
Dejar en blanco si el cdigo de forma = 1.
66-70: Coordenada-P del punto de balance B, como una proporcin de Pyc.
Dejar en blanco si el cdigo de forma = 1.

E2(e) PATRONES DE FUERZAS DE FIJACION EN LOS EXTREMOS (2I5, 7F10.0)
UNA LINEA DE DATOS PARA CADA PATRON.
Omitir si no hay fuerzas de fijacin. Ver la figura 3.10.
Columnas
01-05: Nmero del patrn, en secuencia comenzando con 1.
06-10: Cdigo del sistema de ejes, como sigue:
Cdigo = 0: Las fuerzas estn en el sistema de coordenadas del
elemento, como en la figura 3.10a.
Cdigo = 1: Las fuerzas estn en el sistema de coordenadas
global, como en la figura 3.10b.
11-20: Fuerza de fijacin, Fi.
21-30: Fuerza de fijacin, Vi.
31-40: Momento de fijacin, Mi.
41-50: Fuerza de fijacin, Fj.
51-60: Fuerza de fijacin, Vj.

151
61-70: Momento de fijacin, Mj.
71-80: Factor de reduccin de la carga viva (menor igual a 1) para el clculo
de las fuerzas de la carga viva para aplicarse a los nudos. Si se
coloca 1, se considerara el 100% de la carga viva. Ver la seccin
3.5.1 (f) para mayor explicacin.

E2(f) PATRONES DE FUERZAS INICIALES EN EL ELEMENTO I5,6F10.0)
UNA LINEA DE DATOS PARA CADA PATRON DE FUERZAS INICIALES
Omitir si no hay fuerzas iniciales. Ver la figura 3.10.
Columnas
01-05: Nmero del patrn, en secuencia comenzando con 1.
06-15: Fuerza axial inicial, Fi.
16-25: Fuerza cortante inicial, Vi.
26-35: Momento inicial, Mi.
36-45: Fuerza axial inicial, Fj.
46-55: Fuerza cortante inicial, Vj.
56-65: Momento inicial, Mj.

E2(g) COMANDOS PARA GENERACION DE ELEMENTOS (12I5, 2F5.0, I5, F5.0)
UNA LINEA DE DATOS PARA CADA COMANDO DE GENERACION
Los elementos deben especificarse en un orden numrico creciente. Deben incluirse
obligatoriamente los datos del primer y ltimo elemento. Ver OBSERVACION 7 para mayor
explicacin del procedimiento de generacin.
Columnas
01-05: Nmero del elemento, o nmero del primer elemento en una serie
numerada secuencialmente de elementos para ser generados por
ste comando.
06-10: Nmero del nudo en el extremo i del elemento.
11-15: Nmero del nudo en el extremo j del elemento.

152
16-20: Incremento del nmero de nudo para la generacin del elemento. Si
es cero o est en blanco; se asume igual a 1.
21-25: Nmero del tipo de rigidez.
26-30: Nmero de tipo de excentricidad en los extremos. Dejar en blanco o
poner cero si no hay excentricidad en los extremos.
31-35: Nmero de superficie de fluencia en el extremo i del elemento.
36-40: Nmero de superficie de fluencia en el extremo j del elemento.
41-45: Cdigo para incluir la rigidez geomtrica. Poner 1 si la rigidez
geomtrica debe incluirse. Dejar en blanco o poner cero si la rigidez
geomtrica debe ignorarse.
46-50: Cdigo de salida de la historia en el tiempo. Si no se requiere los
resultados tiempo-historia del elemento, para los elementos
comprendidos por ste comando, poner cero o dejar en blanco. Si
se requiere una salida de la historia en el tiempo, en los intervalos
especificados en la lnea de datos D1, poner 1.
51-55: Nmero del patrn de fuerzas de fijacin para las cargas muertas
estticas en el elemento. Dejar en blanco o poner cero si no hubiera
cargas muertas. Ver la siguiente NOTA a continuacin.
56-60: Nmero del patrn de fuerzas de fijacin para las cargas vivas
estticas en el elemento. Dejar en blanco o poner cero si no hubiera
cargas vivas.
61-65: Factor de escala que debe ser aplicado a las fuerzas de fijacin
debida a cargas muerta estticas. Dejar en blanco si no hubiera
cargas muertas.
66-70: Factor de escala que debe ser aplicado a las fuerzas de fijacin
debido a cargas vivas estticas. Dejar en blanco si no hubiera
cargas vivas.
71-75: Nmero del patrn de fuerzas iniciales. Dejar en blanco o poner cero
si no hubiera fuerzas iniciales.

153
76-80: Factor de escala que debe ser aplicado a las fuerzas iniciales del
elemento.

NOTA: Si el cdigo de carga esttica, lnea de dato C1 es cero, pero fuerzas de fijacin
estn especificadas an para algunos elementos, entonces se origina una
incongruencia. En efecto, cualquiera de estas fuerzas de fijacin ser tratada como
fuerzas iniciales en el elemento.

154
E3. ELEMENTO MURO DE TABIQUE.
Ver el apndice B3 de la referencia 24 para la descripcin del elemento. Nmero de
posiciones de informacin por elemento = 42.

E3(a) INFORMACION DE CONTROL PARA EL GRUPO (3I5) - UNA LINEA DE DATOS
Columnas
01-05: Poner 3 (para indicar que el grupo consta de ELEMENTOS MURO
DE TABIQUE).
06-10: Nmero de elementos en el grupo.
11-15: Nmero de diferentes tipos de rigidez del elemento. (Mximo 40).
Ver seccin E3 (b).

E3(b) TIPOS DE RIGIDEZ (I5, 5F10.0,I5)
UNA LINEA DE DATOS PARA CADA TIPO DE RIGIDEZ.
Columnas:
01-05: Nmero del tipo de rigidez, en secuencia comenzando con 1.
06-15: Mdulo de corte elstico.
16-25: Mdulo de corte de endurecimiento por deformacin, como una
proporcin del mdulo de corte elstico.
26-35: Espesor promedio del panel.
36-45: Esfuerzo de fluencia en corte.
46-55: Deformacin en la falla total, como una proporcin de la deformacin
en la fluencia. Este no debe ser menor que 1.0.
56-60: Cdigo de falla, que dirige el tipo de comportamiento despus de la
falla. Poner 1 si la resistencia y la rigidez deben reducirse a cero
despus de la falla. Poner cero o dejar en blanco si la resistencia y
la rigidez del componente elstico (endurecimiento por deformacin)
debe retenerse despus de la falla.

E3(c) COMANDOS PARA GENERACION DE ELEMENTOS (8I5, F10.0)

155
UNA LINEA DE DATOS PARA CADA COMANDO DE GENERACION.
Los elementos deben especificarse en un orden numrico creciente. Deben incluirse
obligatoriamente los datos del primer y ltimo elemento. Ver la OBSERVACION 7 para
mayor explicacin del procedimiento de generacin.
Columnas
01-05: Nmero del elemento, o nmero del primer elemento en una serie
numerada secuencialmente de elementos para ser generados por
ste comando.
06-10: Nmero del nudo i. Ver la siguiente NOTA a continuacin.
11-15: Nmero del nudo j.
16-20: Nmero del nudo k.
21-25: Nmero del nudo l.
26-30: Incremento del nmero del nudo para la generacin del elemento. Si
es cero o est en blanco, se asume igual a 1.
31-35: Nmero del tipo de rigidez.
35-40: Cdigo de salida de historia en el tiempo. Si no se requiere los
resultados tiempo-historia del elemento para los elementos
comprendidos por ste comando, poner cero o dejar en blanco. Si
se requiere una salida de la historia en el tiempo, en los intervalos
especificados en la lnea de datos D1, poner 1.
41-50: Esfuerzo de corte inicial en el elemento.

NOTA: Los nmeros de los nudos i, j, k, y l; deben enumerarse secuencialmente en sentido
antihorario alrededor de cada elemento con respecto al eje Z positiva.

156
E4. ELEMENTO DE CONEXION SEMI-RIGIDO.
Ver el apndice B4 de la referencia 24 para descripcin del elemento. Nmero de
posiciones de informacin por elemento = 25.

E4(a) INFORMACION DE CONTROL PARA EL GRUPO (3I5) - UNA LINEA DE DATOS.
Columnas:
01-05: Poner 4 (para indicar que el grupo consta de ELEMENTOS DE
CONEXION SEMI-RIGIDO).
06-10: Nmero de elementos en el grupo.
11-15: Nmero de diferentes tipos de rigidez del elemento (mximo 40). Ver
seccin E4(b).

E4(b) TIPOS DE RIGIDEZ (I5,4F10.0)
UNA LINEA DE DATOS PARA CADA TIPO DE RIGIDEZ.
Columnas:
01-05: Nmero de tipo de rigidez, en secuencia comenzando con 1.
06-15: Rigidez rotacional inicial (momento por radian).
16-25: Rigidez de endurecimiento por deformacin, como una proporcin de
la rigidez rotacional inicial.
26-35: Momento de fluencia positivo. Ver la siguiente NOTA a continuacin.
36-45: Momento de fluencia negativo.

NOTA: Rotacin positiva, es la rotacin del nudo i en sentido antihorario con respecto al
nudo j cuando se ve en la direccin z positiva. Un momento positivo en la
conexin tiende a producir una rotacin positiva.

E4(c) COMANDOS PARA GENERACION DE ELEMENTOS (6I5, F1O.O)
UNA LINEA DE DATOS PARA CADA COMANDO DE GENERACION

157
Los elementos deben especificarse en orden numrico creciente. Se deben incluir
obligatoriamente los datos del primer y ltimo elemento. Ver la OBSERVACION 7 para
mayor explicacin del procedimiento de generacin.
Columnas:
01-05: Nmero del elemento o nmero del primer elemento en una serie
numerada secuencialmente de elementos para ser generados por
ste comando.
06-10: Nmero de nudo i
11-15: Nmero del nudo j.
16-20: Incremento del nmero de nudo para la generacin del elemento. Si
es cero o est en blanco, se asume igual a 1.
21-25: Nmero del tipo de rigidez
26-30: Cdigo de salida de la historia en el tiempo, si no se requiere los
resultados tiempo-historia para los elementos comprendidos por ste
comando, poner cero o dejar en blanco. Si se requiere una salida de
la historia en el tiempo, en los intervalos especificados en la lnea de
datos D1, poner 1.
31-40: Momento inicial en la conexin.

158
E5. ELEMENTO VIGA
El ELEMENTO VIGA es idntico al ELEMENTO VIGA-COLUMNA descrito en el
Apndice B2 de la referencia 24, excepto que solamente el cdigo de fluencia tipo VIGA es
permitido. Nmero de posiciones de informacin por elemento = 97

El ingreso de datos es idntico al ELEMENTO VIGA-COLUMNA (seccin E2), excepto
las siguientes partes:
1. Poner 5 en la columna 1-5 de la lnea de datos E2(a), para indicar que el grupo consta
de ELEMENTOS VIGA.
2. El trmino superficies de interaccin de fluencia se sustituye por valores del momento
de fluencia
3. Los datos de la superficie de fluencia en las lneas de datos E2(d) son idnticos. Sin
embargo se ignoran los datos de las columnas 6-10 y 31-70, y se asume
automticamente una superficie de fluencia tipo VIGA.

159
E6. ELEMENTO VIGA CON RIGIDEZ DEGRADANTE
Ver el suplemento de Agosto 1975 de la referencia 24, para la descripcin del elemento.
Nmero de posiciones de informacin por elemento = 165.

NOTA SOBRE LA CONVENCION DE SIGNOS
La convencin de signos utilizada en la salida de los momentos flectores para el
elemento viga con rigidez degradante difiere de los otros elementos viga.
1) otros elementos:
i j
Como se indico en el reporte

2) Elementos viga con rigidez degradante:
i j

E6(a) INFORMACION DE CONTROL PARA EL GRUPO (7I5) - UNA LINEA DE DATOS
Columnas
01-05: Poner 6 (para indicar que el grupo consta de elementos VIGA
CON RIGIDEZ DEGRADANTE).
06-10: Nmero de elementos en el grupo.
11-15: Nmero de tipos diferentes de rigidez del elemento (mx.=40). Ver
seccin E6 (b).
16-20: Nmero de tipos diferentes de excentricidad en los extremos
(mx.=15). Ver seccin E2 (c).
21-25: Nmero de valores diferentes de momento de fluencia para
secciones transversales (mx.=40). Ver seccin E2 (d).

160
26-30: Nmero de patrones diferentes de fuerzas de empotramiento en los
extremos (mx.=34). Ver seccin E2 (e)
31-35: Nmero de patrones diferentes de fuerzas iniciales en el elemento
(mx.=30). Ver seccin E2 (f).

E6(b) TIPOS DE RIGIDEZ - DOS LINEAS PARA CADA TIPO DE RIGIDEZ
LINEA 1: PROPIEDADES DE LA VIGA (I5, 3F10 .0, 3F5 .0, 2F10 .0).
Columnas
01-05 Nmero del tipo de rigidez, en secuencia comenzando por 1.
06-15: EI, rigidez de flexin referencial
16-25: EA, rigidez axial efectiva
26-35: GA, Rigidez de corte efectiva. Si esta en blanco o es cero, se
desprecia las deformaciones por cortante.
36-40: Factor de rigidez de flexin Kii.
41-45: Factor de rigidez de flexin Kjj.
46-50: Factor de rigidez de flexin Kij.
51-60: Relacin de endurecimiento por deformacin para la flexin
inelstica en el nudo i. Si se especifica una rigidez de la rtula
distinta a cero para el nudo i en la LINEA 2, columnas 6-15, entonces
esta relacin de endurecimiento por deformacin se aplicar
directamente a la relacin del momento-rotacin de la rtula. En
caso contrario, esta relacin se aplicar a la relacin de momento-
rotacin en el extremo de la viga total o la relacin de la viga
voladiza.
P
61-70: Relacin de endurecimiento por deformacin para la flexin
inelstica en el nudo j. De la misma forma que en el nudo i, una
rigidez de la rtula cero o distinta de cero para el nudo j en la
LINEA 2, columnas 16-25, controlar si esta relacin se aplica
directamente slo a la rtula o a toda la viga.


161
LINEA 2: PROPIEDADES DE LA ROTULA (I5, 7F10.0)
Dejar en blanco todas las columnas para obtener el modelo de Takeda en forma
automtica.
Columnas
01-05: El nmero del tipo de rigidez, en secuencia comenzando por 1 y que
corresponde al nmero del tipo de rigidez en la LINEA DE
PROPIEDADES DE VIGA de la seccin anterior.
06-15: Rigidez de la rtula en el nudo i. Dejar en blanco o en cero si las
propiedades de la rtula sern determinadas por el programa (estas
propiedades se marcan con * en la salida). Si est en blanco o en
cero, el siguiente campo para el nudo j debe, tambin, estar en
blanco o cero. Si es distinto a cero, el siguiente campo, tambin,
debe ser distinto a cero.
16-25: Rigidez de la rtula en el nudo j.
26-35: , parmetro de rigidez de descarga para el extremo i. Dejar en
blanco o en cero para descargar segn el modelo de Takeda.
i
36-45: , parmetro de rigidez de descarga para el extremo j. Dejar en
blanco o en cero para descargar segn el modelo de Takeda.
j
46-55: , parmetro de carga para el extremo i. Dejar en blanco o en
cero para recargar segn el modelo de Takeda.
i
56-65: , parmetro de carga para el extremo j. Dejar en blanco o en
cero para recargar segn el modelo de Takeda.
j
66-75: N, parmetro exponencial de carga.

Las secciones desde E6(c) hasta E6(g) son idnticas a las secciones E2(c) hasta E2(g).
Sin embargo en la lnea de datos E2(d) se ignoran los datos de las columnas 6-10 y 31-70.

162
F. SIGUIENTE PROBLEMA.
Los datos para un nuevo problema pueden seguir inmediatamente, empezando por la
lnea de datos de la seccin A.

G. LINEA DE DATOS DE FINALIZACION (A-4)
UNA LINEA PARA FINALIZAR EL INGRESO DE DATOS.
Columnas:
01-04: Poner la palabra STOP.

163
OBSERVACIONES

OBSERVACION 1: ESPECIFICACION DE LA COORDENADA NODAL
Las coordenadas del nudo de control se deben definir con respecto al sistema de
coordenadas X-Y. Las coordenadas de los nudos restantes se generan utilizando los
comandos de generacin lineal (seccin B3). El nmero de nudos generados por cada
comando puede ser uno o cualquier nmero mayor. Previamente, se deben definir las
coordenadas de los dos nudos al comienzo y final de la lnea de generacin ya sea por
especificacin directa o por previa generacin lineal.
No es necesario proporcionar comandos de generacin para nudos que:
a) Son numerados secuencialmente entre los nudos inicial y final de cualquier serie lineal,
y
b) Tienen un mismo espacio a lo largo de esa lnea.

Despus que todos los comandos de generacin han sido ejecutados, las coordenadas
para cada grupo de nudos no especificados directamente son generadas automticamente
asumiendo una numeracin secuencial y un espaciado igual a lo largo de las lneas que
unen a los nudos especificados que preceden y el siguiente grupo. Es decir, resultara
superfluo cualquier comando de generacin con espacio igual y con una diferencia del
nmero de nudo igual a uno.

OBSERVACION 2: NUDOS CON DESPLAZAMIENTOS CERO
Cada nudo de la estructura puede tener hasta tres grados de libertad, es decir,
desplazamiento X, desplazamiento Y, y rotacin. Estos son desplazamientos relativos al
terreno.
Inicialmente, el programa asume que los tres grados de libertad estn presentes en
todos los nudos (cdigo = 0) y despus inicializa los arreglos de los datos. Si este supuesto
es correcto, las lneas de datos de la seccin B4 deberan omitirse. Sin embargo, en
algunos casos a) ciertos nudos pueden ser fijados relativamente al terreno en ciertas
direcciones o b) en forma razonable se puede asumir el desplazamiento cero. Cualquier

164
grado de libertad que sea fijado debe de asignrsele un cdigo = 1, y lneas de datos se
deben incluir en la seccin B4 para especificar esos nudos y grados de libertad para los
que el cdigo sea igual a1.
Si existiera cualquier duda, se debera asumir que todos los nudos pueden desplazarse
con los tres grados de libertad (es decir para todos el cdigo = 0). Sin embargo, si ciertos
grados de libertad se eliminaran, el tiempo de cmputo se reducira significativamente.

OBSERVACION 3: NUDOS CON DESPLAZAMIENTOS IDENTICOS
Con frecuencia, puede ser razonable asumir que ciertos nudos se desplazan en forma
idntica en ciertas direcciones. Con los comandos de la seccin B5 se pueden especificar
desplazamientos idnticos. El formato de entrada para esta seccin restringe a 14 el
nmero de nudos cubiertos por un solo comando. Si ms de 14 nudos tienen que ser
asignados a tener desplazamientos idnticos, se deberan utilizar dos o ms comandos con
los nudos en orden numrico creciente en cada comando y con el nudo de nmero menor
comn a todos los comandos.
Como con los desplazamientos que son restringidos a cero, mayor eficiencia de
cmputo se puede ejecutar por la especificacin de desplazamientos idnticos. Sin
embargo, el ancho de la banda de la matriz de rigidez de la estructura siempre decrecer o
se mantendr inalterada con la especificacin de desplazamientos cero, pero con la
especificacin de los desplazamientos idnticos puede incrementarse el ancho de la banda.
Este efecto incrementar el requerimiento de almacenamiento para la rigidez de la
estructura y/o el esfuerzo de cmputo requerido para resolver las ecuaciones de equilibrio.
Por lo tanto, se deberan especificar los desplazamientos idnticos con precaucin y sus
efectos sobre los requerimientos de almacenamiento y los tiempos de ejecucin se
deberan investigar. Debe mencionarse que la solucin de las ecuaciones utilizado en el
programa es mucho menos sensible a los incrementos locales en el ancho de banda de la
matriz de rigidez, a aquella solucin basado en un esquema de almacenamiento por
bandas.



165
OBSERVACION 4: ESPECIFICACION DE MASAS CONCENTRADAS
Los comandos de especificacin para masas concentradas permitirn, por lo general, al
usuario ingresar las masas nodales con slo unas cuantas lneas de datos. Si se desea
cualquier nudo puede aparecer en ms de un comando de especificacin. En tales casos la
masa asociada a un determinado grado de libertad ser la suma de las masas
especificadas en los comandos separados. Si ciertos nudos son forzados a tener un
desplazamiento idntico, la masa asociada con ste desplazamiento ser la suma de las
masas especificadas por los nudos individuales.
Observe que las masas tienen que ingresar en unidades de masa antes que en
unidades de peso, pero pueden ingresarse valores de peso si se especifica el factor de
modificacin apropiado.

OBSERVACION 5: ESCALA DE REGISTROS DE SISMOS
Factores de escala de aceleracin se pueden utilizar para incrementar o para disminuir
las aceleraciones del terreno, o para convertir valores mltiplos de la aceleracin de la
gravedad a unidades de aceleracin, o ambos. La modificacin de la intensidad ssmica
mediante el escalamiento de los valores de aceleracin es una prctica comn en
investigaciones, pero, se debera tomar cuidadosamente en aplicaciones prcticas. Cuando
se escalan las aceleraciones, tambin, se escalan los desplazamientos y las velocidades
del terreno en la misma proporcin.
Tambin, se pueden modificar la escala del tiempo. Si se especifica un factor de escala
de tiempo por decir igual a f, todos los tiempos de entrada se multiplican por f antes de
calcular las aceleraciones interpoladas en intervalos iguales al paso de integracin. Si
permanecen inalteradas las aceleraciones del terreno, el efecto que tendr ser incrementar
las velocidades del terreno por f y los desplazamientos por y alterar el contenido de
la frecuencia del sismo. Las modificaciones de escala de tiempo no deberan hacerse sin
considerar cuidadosamente su influencia en el movimiento del terreno.
f
2




166
OBSERVACION 6: AMORTIGUAMIENTO
Pueden especificarse 4 tipos diferentes de amortiguamiento, en forma simple o en
cualquiera combinacin. Estos tipos son los siguientes:
1. Amortiguamiento proporcional a la masa, en la que se asume una matriz de
amortiguamiento viscoso, C = M, donde M = es la matriz de masa. El procedimiento
para calcular se revisa en la referencia 14.
2. Amortiguamiento proporcional a la rigidez, en la que se asume una matriz de
amortiguamiento viscoso, C = K
T
, donde K
T
= es la matriz de rigidez tangencial actual
en un determinado instante que incluye tambin la rigidez geomtrica.
3. Amortiguamiento proporcional a la rigidez, en la que se asume una matriz de
amortiguamiento viscoso, C =
O
K
O
, donde K
O
=es la matriz de rigidez elstica original,
ignorando la rigidez geomtrica. Ver el captulo 2 de la referencia 24 para mayor
explicacin.
4. Amortiguamiento estructural, en la que se asume fuerzas de amortiguamiento tal que:
Las magnitudes de las fuerzas de amortiguamiento que se originan en cualquier
elemento son un mltiplo, , de los valores absolutos de las acciones del elemento
(fuerzas axiales, momentos en los extremos, etc.) y
La direccin de cada accin de amortiguamiento es tal que sta se opone al cambio
de la deformacin del elemento correspondiente.
Se asume una matriz de amortiguamiento no viscosa. Fuerzas de amortiguamiento
instantneas son aplicadas en cualquier paso de integracin basadas en las acciones
del miembro y en los valores de deformacin que existen al final del paso anterior. Las
acciones del miembro utilizadas para calcular las fuerzas de amortiguamiento son las
acciones totales que incluyen tanto los efectos estticos como los dinmicos. Las
fuerzas de amortiguamiento empiezan aplicarse a partir del segundo paso de
integracin. Ver captulo 2 de la referencia 24 para mayor explicacin.
Se ha realizado poco uso del concepto de amortiguamiento estructural y por eso poca
experiencia est disponible, especialmente para las estructuras inelsticas. Un
problema posible de las fuerzas de amortiguamiento es que pueden contribuir a
acentuar pequeas oscilaciones en los clculos numricos, ya que las fuerzas de

167
amortiguamiento en cualquier paso se basan siempre en el estado final del paso
anterior. Hasta que se alcance la necesaria experiencia, esta opcin se deber utilizar
cuidadosamente. Principalmente por razones de investigacin, esta opcin se incluye
en el programa.

OBSERVACION 7: GENERACION DE DATOS DEL ELEMENTO
En los comandos de generacin del elemento, los elementos se deben especificar en
orden numrico creciente. Se pueden proporcionar lneas de datos para elementos
numerados secuencialmente, en cuyo caso cada lnea de dato especifica un elemento y no
se utiliza la opcin de generacin. Alternativamente, las lneas de datos para un grupo de
elementos se pueden omitirse, en cuyo caso los datos para el grupo omitido se genera de la
siguiente manera:
1. A todos los elementos del grupo se les asigna las mismas rigideces, resistencias, datos
de carga, cdigo de salida, etc., tal como el elemento anterior al grupo omitido.
2. Los nmeros de los nudos para cada elemento omitido, se obtienen por adicin del
incremento del nmero de nudo especificado a los nmeros de los nudos del elemento
precedente. El incremento del nmero de nudo es la que se especific para el
elemento, anterior al grupo omitido.

En la salida de los datos del elemento, los datos generados son identificados por un
asterisco al inicio de la lnea impresa.








ANEXO B


EJEMPLOS DE ARCHIVOS DE DATOS


INTRODUCCION
En la tabla B.1, se presenta el archivo para el ingreso de datos para el programa
DRAIN-2D, que corresponde al prtico P-5, cuya geometra se presenta en la figura 4.1(a),
en el cuadro 4.1 se presenta las dimensiones de los elementos del prtico y en los cuadros
4.4 y 4.5 se presentan los aceros de refuerzos y las propiedades de fluencia de las
secciones tpicas del prtico. El prtico est sometido al sismo 10/03/74 Comp. N82W.

Los archivos de ingreso de datos, correspondientes a los 19 prticos sometidos a 6
sismos cada uno, las cuales fueron utilizados en el presente estudio, se adjuntan en un CD,
conjuntamente con el programa DRAIN-2D, versin para computadoras personales.
169
TABLA B.1: ARCHIVO DE INGRESO DE DATOS PARA EL PORTICO P-5

START PORTICO P-5F DE 5 PISOS Y 4 CRUJIAS, ELEMENTO VIGA-COLUMNA
C
C MODELADO CON ELEMENTO VIGA-COLUMNA
C M.F.B JUNIO-2001, UNIDADES: M-TON-SEG
C
C INFORMACION DE CONTROL CHECK
30 12 6 1 5 5 1 0 0 0
C COORDENADAS DE NUDOS DE CONTROL (12)
C ND X Y
1 0 0
5 20.00 0
6 0 3.00
10 20.00 3.00
11 0 5.60
15 20.00 5.60
16 0 8.20
20 20.00 8.20
21 0 10.80
25 20.00 10.80
26 0 13.40
30 20.00 13.40
C GENERACION DE NUDOS (6)
1 5 3 1
6 10 3 1
11 15 3 1
16 20 3 1
21 25 3 1
26 30 3 1
C NUDOS CON DESPLAZAMIENTO CERO (1)
1 1 1 1 5 1
C NUDOS CON DESPLAZAMIENTOS IGUALES (5)
1 5 6 7 8 9 10
1 5 11 12 13 14 15
1 5 16 17 18 19 20
1 5 21 22 23 24 25
1 5 26 27 28 29 30
C MASAS CONCENTRADAS EN NUDOS (5)
C Ndo MX MY MR Ndo DIF FACTOR
6 10.07 0.0 0.0 0 0 1.0
11 10.07 0.0 0.0 0 0 1.0
16 10.07 0.0 0.0 0 0 1.0
21 10.07 0.0 0.0 0 0 1.0
26 7.29 0.0 0.0 0 0 1.0
C INFORMACION PARA CONTROL DE CARGAS
C nt t escalaX escalaT escalaY escalaT DESmax
1 5 19591 0.005 0.000020 0.0 0.0 0.0 100.0
C CARGAS ESTATICAS EN NUDOS
C Ndo Fx FY Mz Ndo DIF
6 0 -5.44 0 10 1
11 0 -5.44 0 15 1
16 0 -5.44 0 20 1
21 0 -5.44 0 25 1
26 0 -3.90 0 30 1
C INFORMACION DEL SISMO
4899 0 0 0 LIMA PERU INSTITUTO GEOFISICO DEL PERU 10/03/74, N82WCOMP
C REGISTRO DEL SISMO Amax=-192.50 CM/SEG2 t=19.78 SEG
0.00 0.0 0.02 1717 0.04 -1710 0.06 -4183 0.08 -901
0.10 2650 0.12 5905 0.14 -618 0.16 -4281 0.18 -1217
0.20 1588 0.22 -559 0.24 11391 0.26 3414 0.28 -18184
0.30 -7851 0.32 5874 0.34 3672 0.36 5903 0.38 17912
0.40 19885 0.42 2783 0.44 -9430 0.46 -11768 0.48 -8086
0.50 2759 0.52 2504 0.54 -1767 0.56 5149 0.58 13143
0.60 8226 0.62 -5537 0.64 -12830 0.66 -4356 0.68 2933
0.70 -67 0.72 5936 0.74 9801 0.76 -5077 0.78 -8686
0.80 3633 0.82 9037 0.84 9159 0.86 8298 0.88 -7703
0.90 -27844 0.92 -8506 0.94 541 0.96 -19313 0.98 -19161
1.00 6378 1.02 21025 1.04 13381 1.06 4187 1.08 -1741
1.10 5019 1.12 9332 1.14 -8100 1.16 -12272 1.18 9711
1.20 19593 1.22 4562 1.24 -5563 1.26 6753 1.28 13197
1.30 -4839 1.32 -4626 1.34 3609 1.36 -10269 1.38 -7629
1.40 7900 1.42 5233 1.44 48 1.46 19806 1.48 27365
1.50 -432 1.52 -17239 1.54 -15677 1.56 -2208 1.58 7688
1.60 -138 1.62 -5770 1.64 -1205 1.66 -1623 1.68 -8187
1.70 -4012 1.72 11061 1.74 21156 1.76 30601 1.78 43087
1.80 32608 1.82 -9261 1.84 -29835 1.86 -17335 1.88 -6423
1.90 -3196 1.92 3587 1.94 11705 1.96 -9257 1.98 -27899
2.00 -10374 2.02 14964 2.04 26596 2.06 21524 2.08 -695
2.10 -11418 2.12 -17174 2.14 -28805 2.16 -20971 2.18 -900

170
2.20 4854 2.22 -2727 2.24 11445 2.26 34022 2.28 21504
2.30 7341 2.32 -192 2.34 -13062 2.36 -13890 2.38 10631
2.40 21488 2.42 5993 2.44 -6795 2.46 -9728 2.48 913
2.50 6338 2.52 -842 2.54 1306 2.56 4783 2.58 -3247
2.60 -10719 2.62 -24683 2.64 -26980 2.66 -7835 2.68 13242
2.70 32533 2.72 42545 2.74 39082 2.76 3081 2.78 -30280
2.80 -5748 2.82 25542 2.84 12510 2.86 526 2.88 11532
2.90 25634 2.92 19925 2.94 -4853 2.96 -13931 2.98 -4358
3.00 20859 3.02 13133 3.04 -5960 3.06 -7910 3.08 -21327
3.10 -18941 3.12 -1604 3.14 9275 3.16 31572 3.18 28476
3.20 -1747 3.22 -1792 3.24 13830 3.26 20318 3.28 10046
3.30 7753 3.32 -6779 3.34 -43096 3.36 -49703 3.38 -16692
3.40 43880 3.42 54275 3.44 15374 3.46 13129 3.48 8295
3.50 -30631 3.52 -39743 3.54 -8746 3.56 1680 3.58 -28634
3.60 -43024 3.62 -23649 3.64 3639 3.66 22868 3.68 28661
3.70 8844 3.72 9988 3.74 31050 3.76 17101 3.78 19467
3.80 28572 3.82 -17254 3.84 -40364 3.86 -11492 3.88 1844
3.90 -5217 3.92 6945 3.94 7867 3.96 -27725 3.98 -49185
4.00 -19490 4.02 43169 4.04 51280 4.06 12950 4.08 -17432
4.10 -42068 4.12 -35120 4.14 -971 4.16 13546 4.18 5630
4.20 10147 4.22 7570 4.24 -21506 4.26 -19301 4.28 13971
4.30 20159 4.32 -4042 4.34 -17301 4.36 4443 4.38 5374
4.40 -13965 4.42 -6659 4.44 28866 4.46 18586 4.48 -7982
4.50 -16223 4.52 -17371 4.54 3116 4.56 22795 4.58 3737
4.60 -36722 4.62 -44817 4.64 -33502 4.66 -23737 4.68 -8875
4.70 19357 4.72 32661 4.74 23997 4.76 12799 4.78 8825
4.80 8384 4.82 16348 4.84 34353 4.86 46652 4.88 18220
4.90 -17725 4.92 -27504 4.94 -29950 4.96 -25627 4.98 -2461

. . . . .
. . . . .
. . . . .
. .(contina datos del sismo) . .
. . . . .
. . . . .
. . . . .
91.80 864 91.82 3658 91.84 3602 91.86 877 91.88 1575
91.90 2429 91.92 1186 91.94 -901 91.96 -3580 91.98 -5939
92.00 -5367 92.02 -2629 92.04 -168 92.06 388 92.08 742
92.10 2745 92.12 5453 92.14 4985 92.16 1548 92.18 570
92.20 1821 92.22 -219 92.24 -1742 92.26 994 92.28 3261
92.30 286 92.32 -2070 92.34 -4414 92.36 -7917 92.38 -6926
92.40 -3504 92.42 959 92.44 4246 92.46 5577 92.48 6292
92.50 3942 92.52 488 92.54 1495 92.56 4139 92.58 4558
92.60 4182 92.62 3214 92.64 1011 92.66 -4646 92.68 -8285
92.70 -5427 92.72 -2527 92.74 -4753 92.76 -6209 92.78 -3474
92.80 438 92.82 3433 92.84 2311 92.86 -771 92.88 -679
92.90 2386 92.92 3382 92.94 695 92.96 496 92.98 3036
93.00 3867 93.02 1572 93.04 -3525 93.06 -4630 93.08 -1839
93.10 1170 93.12 804 93.14 -717 93.16 2470 93.18 4354
93.20 2699 93.22 4292 93.24 4577 93.26 1476 93.28 2636
93.30 4451 93.32 1277 93.34 1979 93.36 2341 93.38 -428
93.40 1156 93.42 1946 93.44 574 93.46 -213 93.48 -2208
93.50 -2041 93.52 1342 93.54 -337 93.56 -3949 93.58 -2575
93.60 -461 93.62 -2134 93.64 -6064 93.66 -8056 93.68 -3525
93.70 2946 93.72 5315 93.74 3683 93.76 697 93.78 -1717
93.80 -3709 93.82 -5039 93.84 -3564 93.86 -1566 93.88 -2618
93.90 -2948 93.92 205 93.94 4723 93.96 5705 93.98 1915
94.00 937 94.02 3237 94.04 2067 94.06 -1019 94.08 -1318
94.10 984 94.12 4029 94.14 2539 94.16 -1693 94.18 -3101
94.20 -3506 94.22 -3330 94.24 -1194 94.26 1492 94.28 3337
94.30 4806 94.32 3481 94.34 1838 94.36 3725 94.38 3527
94.40 3221 94.42 3598 94.44 -569 94.46 -1015 94.48 1064
94.50 -923 94.52 102 94.54 3392 94.56 1694 94.58 -2931
94.60 -3331 94.62 -370 94.64 2542 94.66 3114 94.68 503
94.70 570 94.72 2492 94.74 361 94.76 -267 94.78 1820
94.80 -50 94.82 -2440 94.84 -1170 94.86 1492 94.88 -168
94.90 -1321 94.92 3815 94.94 4396 94.96 1672 94.98 4142
95.00 5345 95.02 825 95.04 -2652 95.06 -3379 95.08 -2232
95.10 -370 95.12 -1831 95.14 -5306 95.16 -8038 95.18 -7652
95.20 -4995 95.22 -2336 95.24 -2098 95.26 -2659 95.28 -1911
95.30 -302 95.32 2640 95.34 3519 95.36 -20 95.38 -4008
95.40 -3539 95.42 -173 95.44 2381 95.46 28 95.48 -392
95.50 2777 95.52 1204 95.54 -841 95.56 115 95.58 1581
95.60 3920 95.62 5784 95.64 3501 95.66 -724 95.68 -339
95.70 -896 95.72 -2473 95.74 -859 95.76 -1502 95.78 -939
95.80 -1247 95.82 -181 95.84 784 95.86 -1597 95.88 914
95.90 447 95.92 -2007 95.94 951 95.96 3953 95.98 3154
96.00 561 96.02 1678 96.04 4347 96.06 2906 96.08 -646
96.10 128 96.12 4584 96.14 7646 96.16 5339 96.18 -170
96.20 -3096 96.22 -1378 96.24 -1294 96.26 -5444 96.28 -5922

171
96.30 -3233 96.32 -4919 96.34 -6578 96.36 -2648 96.38 2589
96.40 3943 96.42 810 96.44 -1414 96.46 -448 96.48 219
96.50 343 96.52 1874 96.54 4755 96.56 4745 96.58 -105
96.60 -4704 96.62 -4432 96.64 -2087 96.66 -2513 96.68 -3003
96.70 -2660 96.72 -2929 96.74 -4328 96.76 -5162 96.78 -4030
96.80 -2160 96.82 -190 96.84 1809 96.86 2123 96.88 431
96.90 -1830 96.92 -1721 96.94 2085 96.96 3390 96.98 -1890
97.00 -3929 97.02 -220 97.04 2303 97.06 1091 97.08 -2929
97.10 -6250 97.12 -5481 97.14 -1604 97.16 2817 97.18 4101
97.20 1284 97.22 194 97.24 248 97.26 -1023 97.28 -2902
97.30 -3468 97.32 -559 97.34 3058 97.36 720 97.38 -2355
97.40 -943 97.42 1543 97.44 3538 97.46 3196 97.48 505
97.50 -2042 97.52 -1486 97.54 1560 97.56 4302 97.58 3823
97.60 801 97.62 -4181 97.64 -7809 97.66 -7027 97.68 -4921
97.70 -2049 97.72 914 97.74 -509 97.76 -843 97.78 1745
97.80 -245 97.82 -5111 97.84 -5610 97.86 -1161 97.88 -387
97.90 -3482 97.92 -4805 97.94 -2094 97.96 32
C
C COEFICIENTES DE AMORTIGUAMIENTO
C o
0.3973 0.0060 0.0 0.0
C CONTROL PARA RESULTADOS TIEMPO-HISTORIA
C J E ENV NH NV NR NHR NVR THN THR THE ISJ ISE
0 10 0 0 0 0 0 0 0 0 0 0 3
C
C ELEMENTO VIGA-COLUMNA
C
2 45 4 6 6 4 0
C PROPIEDADES DE RIGIDEZ (4)
C N E Eo A I Kii Kjj Kij AS POISSON
C COLUMNAS
1 2200000 0.005 0.240 0.007200 4.0 4.0 2.0 0.2000 0.2
2 2200000 0.005 0.200 0.004167 4.0 4.0 2.0 0.1667 0.2
C VIGAS
3 2200000 0.005 0.150 0.003125 4.0 4.0 2.0 0.1250 0.2
4 2200000 0.005 0.135 0.002278 4.0 4.0 2.0 0.1125 0.2
C
C ZONAS RIGIDAS EN NUDOS (6)
C N Xi Xj Yi Yj
C COLUMNAS
1 0.000 0.000 0.000 0.250
2 0.000 0.000 0.250 0.250
3 0.000 0.000 0.250 0.225
4 0.000 0.000 0.225 0.225
C VIGAS
5 0.300 0.300 0.000 0.000
6 0.250 0.250 0.000 0.000
C
C PROPIEDADES DE FLUENCIA DE SECCIONES (6)
C NS COD My+ My- Pyc Pyt A A B B
C COLUMNAS
1 3 27.27 -27.27 532.88 -109.12 2.17 0.37 2.17 0.37
2 3 17.34 -17.34 438.56 -85.18 2.26 0.37 2.26 0.37
C VIGAS
3 1 12.09 -23.81 0.0 0.0 0.0 0.0 0.0 0.0
4 1 10.60 -20.96 0.0 0.0 0.0 0.0 0.0 0.0
5 1 8.70 -16.90 0.0 0.0 0.0 0.0 0.0 0.0
6 1 6.12 -10.90 0.0 0.0 0.0 0.0 0.0 0.0
C FUERZAS DE EMPOTRAMIENTO (4)
C N CD Fi Vi Mi Fj Vj Mj RLL
C CARGA MUERTA
1 1 0.0 7.48 6.23 0.0 7.48 -6.23 0.0
2 1 0.0 6.78 5.65 0.0 6.78 -5.65 0.0
C CARGA VIVA
3 1 0.0 3.13 2.60 0.0 3.13 -2.60 1.0
4 1 0.0 1.25 1.04 0.0 1.25 -1.04 1.0
C GENERACION DE ELEMENTOS
C ELEM NODE NODE NODE STIF ECCY YIELD MOM GEOM TIME FEF FEF SCALE INITIAL FORCE
C NO. I J DIFF TYPE TYPE ENDI ENDJ STIF HIST DL LL DL LL NO. SCALE
C COLUMNAS
1 1 6 1 1 1 1 1 1 1 0 0 0 0 0 0
6 6 11 1 1 2 1 1 1 1 0 0 0 0 0 0
11 11 16 1 2 3 1 2 1 1 0 0 0 0 0 0
16 16 21 1 2 4 2 2 1 1 0 0 0 0 0 0
21 21 26 1 2 4 2 2 1 1 0 0 0 0 0 0
C VIGAS
26 6 7 1 3 5 3 3 1 1 1 3 1 1 0 0
30 11 12 1 3 5 3 3 1 1 1 3 1 1 0 0
34 16 17 1 4 6 4 4 1 1 1 3 1 1 0 0
38 21 22 1 4 6 5 5 1 1 1 3 1 1 0 0
42 26 27 1 4 6 6 6 1 1 2 4 1 1 0 0

172
45 29 30 0 4 6 6 6 1 1 2 4 1 1 0 0
STOP

173
BIBLIOGRAFIA

1. AZIS, T.; Inelastic Dynamic Analysis of Building Frames; Report R76-37,
Departamento de Ingeniera Civil, MIT; Cambridge, Massachusetts; Agosto 1976.
2. BRADY A. G. y PEREZ, V., Strong Motion Accelerograms-Digitization and Analysis-
Records from Lima, Peru: 1951 to 1974, Open File Report No 77-587, USGS (1977).
3. BUILDING CODE REQUIREMENTS FOR STRUCTURAL CONCRETE (ACI 318-99)
American Concrete Institute, Detroit, 1999.
4. CHAVEZ, J.; Estudio del comportamiento Inelstico de estructuras de concreto Armado
durante Sismos Severos; Ponencia E-26, IV Congreso Nacional de Ingeniera civil;
Chiclayo, Per; 1982.
5. CLOUGH, R.W., and PENZIEN, J., Dynamics of Structures, McGraw-Hill, New York,
1975.
6. CRUZ Z. ERNESTO, GUZMAN G. MARCELO, Estudio del Efecto de la Razn de
Rigidez viga-columna en la Respuesta Nolineal de Estructuras Tridimensionales,
Publicacin DIE 95-1, Departamento de Ingeniera Estructural y Geotcnica, Pontificia
Universidad Catlica de Chile, Enero 1995.
7. DAVIS, G. B., AND HOFFMANN, T. R., FORTRAN 77 un Estilo Estructurado y
Disciplinado Primera edicin en Espaol, McGraw Hill, Mxico, 1984.
8. EL HAFEZ, M. B. and POWELL, G. H., Computer Aided Ultimate Load Design of
Unbraced Multistory Steel Frames, Report No EERC 73-3, Earthquake Engineering
Research Center, University of California, Berkeley, February 1973.
9. FEMA, NEHRP Guidelines for the Seismic Rehabilitation of Buildings, Developed by
the Federal Emergency Management Agency (Report No. FEMA 273), Washington
D.C., 1997.
10. GIBERSON, M.; Two Nonlinear Beams with Definitions of Ductility; Revista ASCE,
Structural Division, Vol 95, No ST2, pp. 137-157; Febrero 1969.
11. GOEL, R. and CHOPRA, A.K., Inelastic Seismic Response of One-Story Asymmetric-
Plan Systems, Report N UCB/EERC-90/14, Earthquake Engineering Research Center,
University of California, Berkeley, October 1990.

174
12. HEJAL, R. and CHOPRA, A.K., Earthquake Response of Torsionally Coupled, Frame
Buildings, Report N UCB/EERC-87/20, Earthquake Engineering Research Center,
University of California, Berkeley, December 1987.
13. HIGGINBOTHAM, M., Computational Procedures for Nonlinear Structural Dynamics,
CE 299 Report, Department of Civil Engineering, University of California, Berkeley,
spring 1971.
14. KABA, S., MAHIN, S.; Refined Modeling of Reinforced Concrete Columns for Seismic
Analysis; Report UCB/EERC/84/03, Universidad de California; Berkeley; April 1984.
15. KANAAN, A.E. and POWELL, G.H., General Purpose Computer Program for Inelastic
Dynamic Response of Plane Structures, Report No. EERC 73-6, Earthquake
Engineering Research Center, University of California, Berkeley, April 1973.
16. MINORU WAKABAYASHI, ENRIQUE MARTINEZ ROMERO, Diseo de Estructuras
Sismorresistentes McGraw-Hill Interamericana de Mxico S.A., 1988.
17. MONDKAR, D.P., and G.H. POWELL, Towards Optimal In-Core Equation Solving,
paper submitted for possible publication in the Journal of Computers and Structures,
February 1973.
18. MUOZ PELAEZ ALEJANDRO, TINMAN BEHAR MARCOS, QUIUN WONG DANIEL,
Riesgo Ssmico de Edificios Peruanos, II Congreso Nacional de Estructuras y
Construccin, Lima, diciembre 2000.
19. NEWMARK, N. M., A method of Computation for Structural Dynamics, Proceedings,
ASCE, Journal of the Engineering Mechanics Division, July 1959.
20. OTANI, S.; Hysteresis Models of Reinforced Concrete For Earthquake Analysis;
Revista de la Facultad de Ingeniera de la Universidad de Tokyo, Vol. XXXVI, No 2,
1981.
21. PARK y PAULAY, Estructuras de Concreto Reforzado, Editorial Limusa S.A., Mxico,
1980.
22. PING-CHUN WANG, Numerical and Matrix Methods in Structural Mechanics with
Applications to Computers, John Wiley and Sons INC., 1966.
23. PIQUE, J.; On the Use of Simple Models in Nonlinear Dynamic; Report R76-43,
Departamento de Ingeniera Civil, MIT; Setiembre 1976.

175
24. POWELL, G.H., DRAIN-2D Users Guide, Report EERC No. 73-22, Earthquake
Engineering Research Center, University of California, Berkeley, September 1973.
25. REGLAMENTO DE CONSTRUCCIONES PARA EL DISTRITO FEDERAL, Propuesta
de Normas Tcnicas Complementarias para el Diseo y Construccin de Estructuras de
Concreto, Mxico, Abril 2001.
26. REGLAMENTO NACIONAL DE CONSTRUCCION, Norma Tcnica de Edificaciones E-
030 Diseo Sismorresistente, 2da. Edicin, Noviembre, 1997.
27. REGLAMENTO NACIONAL DE CONSTRUCCION, Norma Tcnica de Edificaciones E-
060 Concreto Armado, ININVI, 3ra. Edicin, Diciembre, 1990.
28. SAIIDI, M., and SOZEN, M. A., Simple and Complex Models for Nonlinear Seismic
Response of Reinforced Concrete Structures, Civil Engineering Studies, Structural
Research Series, No. 465, University of Illinois, Urbana, August 1979.
29. SUDHAKAR, A., Computer Program for Small Displacement Elasto Plastic Analysis of
Plane Steel and Reinforced Concrete Frames, Graduate Student Report, University of
California, Berkeley, 1972.
30. UMEMURA, H., TAKIZAWA, H.; Dynamic Response of Reinforced Concrete Building;
Asociacin Internacional de Puentes y Estructuras de Ingeniera, IABSE; Suiza; 1982.
31. WILSON E. HABIBULLAH, Structural Analysis Programs (SAP2000 Nonlinear
Version 6.11), Computers and Structures Inc., Berkeley California, 1997.
32. WILSON, E. L., A Method of Analysis for the Evaluation of Foundation-Structure
Interaction, Proceedings, Fourth World Conference on Earthquake Engineering,
Santiago, Chile, January, 1969, pp. A-6, 87 46, 99.
33. WILSON, E. L., Solid SAP - A Static Analysis Program for Three Dimensional Solid
Structures, Report No UC-SESM-71-19, Department of Civil Engineering, Univ. of
California, Berkeley, Sept. 1971.