Estruc Dinam

Departamento de Ingeniera Civil
Facultad Regional Paran Universidad Tecnolgica Nacional

ANALISIS DE ESTRUCTURAS BAJO
ACCIONES DINMICAS

Arturo M. Cassano

Editorial de la Universidad Tecnolgica Nacional - edUTecNe
http://www.edutecne.utn.edu.ar
edutecne@utn.edu.ar

2009
Introduccin - Bibliografa
1

INTRODUCCIN

Esta publicacin pretende brindar en una forma rpida y sencilla los conceptos bsicos de
la dinmica de estructuras aplicada a las construcciones civiles, enfocada desde un punto
de vista numrico. Su organizacin en captulos se basa en la secuencia de unidades tem-
ticas del programa de la asignatura homnima que se dicta en la carrera de Ingeniera Civil
de la Facultad Regional Paran de la Universidad Tecnolgica Nacional.
Todos los temas se desarrollan como sntesis de libros y publicaciones de otros autores
(ver Bibliografa) solo que organizados segn la conveniencia que mejor se ajusta al curso
impartido.
El Captulo 1 hace una muy breve resea a conceptos sobre sismologa y su relacin con la
ingeniera estructural, caracterizacin, registros, etc.
El Captulo 2 brinda una aproximacin a los conceptos bsicos de la dinmica estructural,
la importancia de la masa, la relacin entre velocidad de carga y de reaccin de una estruc-
tura. Avanza sobre los modelos estructurales dinmicos, grados de libertad y la discretiza-
cin espacial y temporal.
El Captulo 3 est dedicado a la caracterizacin y clculo de osciladores de un grado de
libertad dinmico. Son tratadas las vibraciones libres y forzadas, stas ltimas con cargas
armnicas y arbitrarias. Se presenta la resolucin numrica mediante integracin directa.
El Captulo 4 est dedicado a las estructuras con mltiples grados de libertad dinmicos.
Se estudian las vibraciones libres y el clculo de modos y frecuencias propias, tambin se
presenta la resolucin de los sistemas dinmicos mediante diversos mtodos como son:
descomposicin y superposicin modal, integracin directa y respuesta mxima mediante
espectros de respuesta.
En el Captulo 5 se realiza una aproximacin al tratamiento dinmico de los efectos del
viento sobre las construcciones teniendo como base el reglamento CIRSOC 102.
Por ltimo en el Captulo 6 son tratados los efectos ssmicos y el anlisis estructural, desde
el punto de la dinmica estructural, enfatizando en la definicin de la accin ssmica y la
comparacin de las frmulas vistas en captulos anteriores con los reglamentos vigentes.

Un agradecimiento especial a Marlene J aurena, que con paciencia y dedi-
cacin realiz la transcripcin de los manuscritos de este trabajo.
2

TABLA DE CONTENIDOS

Captulo 1 - NOCIONES BASICAS DE SISMOLOGIA 1-1
Causas que generan los terremotos o sismos 1-1
Los sismos desde el punto de vista de la ingeniera y su caracterizacin 1-2
Esteva y Rosenblueth: 1-2
Donovan: 1-3
Esteva y Villaverde: 1-3
Esteva: 1-3
Registro de ondas ssmicas. Parmetros utilizados y mapas de riesgo ssmico 1-3

Captulo 2 - CONCEPTOS BASICOS DE DINAMICA ESTRUCTURAL 2-1
Definicin de la accin dinmica 2-1
Acciones y fuerzas dinmicas 2-2
Importancia de la masa en el problema dinmico 2-3
Velocidad de reaccin de una estructura 2-5
Modelos dinmicos caractersticos 2-6
Mtodos de modelizacin dinmica 2-8
Discretizacin espacial de las estructuras 2-9
Mtodo de las masas concentradas 2-10
Ecuaciones de movimiento 2-13
Principio de Hamilton 2-13
Principio de los trabajos virtuales 2-13
Principio de DAlembert 2-13
Formulacin de la ecuacin de movimiento para un sistema de 1 GLD 2-13
Formulacin de las ecuaciones de movimiento para modelos con mltiples GLD 2-15

Captulo 3 - RESPUESTA DE UN OSCILADOR SIMPLE 3-1
Ecuacin de movimiento y equilibrio dinmico 3-1
Frmula de Geiger 3-2
Caractersticas dinmicas con amortiguamiento 3-3
Excitacin peridica 3-7
Excitacin armnica 3-8
Excitacin arbitraria. Integral de Duhamel 3-10
Respuesta a un impulso elemental 3-10
Factor de amplificacin dinmica 3-12
Espectros ssmicos de respuesta 3-12
Integracin numrica de la ecuacin de movimiento 3-16
Mtodo de Newmark (1 GLD) 3-16
Desarrollo y forma operativa 3-18
3

Captulo 4 - RESPUESTA DINMICA DE UNA ESTRUCTURA CON MLTI-
PLES GRADOS DE LIBERTAD 4-1
Ecuaciones de movimiento y equilibrio dinmico 4-1
Vibraciones libres 4-1
Caractersticas dinmicas 4-1
Normalizacin de los modos 4-3
Obtencin de los grados de libertad dinmicos 4-4
Condensacin esttica de la matriz de rigidez 4-4
Matriz de amortiguamiento 4-6
Matrices de amortiguamiento ortogonales 4-8
Determinacin prctica de modos y frecuencias 4-9
Mtodo de Stodola-Vianello 4-9
Resolucin de las ecuaciones de movimiento en estructuras con mltiples GLD 4-12
Descomposicin y superposicin modal 4-12
Integracin directa de las ecuaciones de movimiento 4-14
Respuesta mxima utilizando espectros de respuesta 4-15

Captulo 5 - ANLISIS DE CONSTRUCCIONES CON EFECTOS DINMI-
COS DE VIENTO 5-1
Acciones paralelas a la direccin del viento 5-1
Acciones perpendiculares a la direccin del viento 5-3
Criterios de confort en edificios que oscilan 5-5

Captulo 6 - ANLISIS DE CONSTRUCCIONES CON EFECTOS
SISMICOS 6-1
Definicin numrica de la accin ssmica 6-1
Definicin mediante espectros de respuesta 6-2
Definicin mediante acelerogramas 6-4
Mtodos de anlisis segn INPRES CIRSOC 103 6-6
Procedimientos con fuerzas estticas equivalentes 6-6
Mtodos dinmicos 6-6
Mtodos dinmicos INPRES CIRSOC 103 6-6
4

BIBLIOGRAFA:

[1] Estructuras Sometidas a Acciones Ssmicas. Clculo por Ordenador A.H.

Barbat, J.M. Canet 2da. Edicin, CIMNE, Barcelona, 1994.

[2] Diseo Sismorresistente de Edificios L.M. Bozzo, A. H. Barbat, Editorial

Revert, Barcelona, 2000.

[3] Dinmica Estructural J. Massa, C. Prato, Publicacin del Departamento de
Estructuras de la Facultad de Cs. Exactas Fsicas y Naturales de la Universidad
Nacional de Crdoba, 1986.

[4] Finite Element Procedures K. J. Bathe, Prentice-Hall, 1996.

[5] Finite Element Modeling in Engineering Practice C. Spyrakos, Algor

Publishing Division, Pitsburg, 1996.

[6] Linear and Nonlinear Finite Element Analysis in Engineering Practice C.

Spyrakos, J. Raftoyiannis, Algor Publishing Division, Pitsburg, 1997.
[7] Diseo Ssmico de Edificios E. Bazn, R. Meli, LIMUSA, 2001.
Captulo 1 Nociones Bsicas de Sismologa
1-1

CAPTULO 1 - NOCIONES BASICAS DE SISMOLOGIA

Causas que generan los terremotos o sismos

Los terremotos pueden definirse como movimientos de la corteza terrestre, con amplitudes
y frecuencias dependientes del tiempo.
Las causas que los generan son variadas:
Terremotos de colapso: son los originados en cavidades subterrneas por el colapso de las
mismas, son de baja intensidad.

Terremotos de origen volcnico: la explosin de gases durante las erupciones volcnicas
puede producir terremotos que, en general, tienen una intensidad pequea y afectan a su-
perficies limitadas.
Terremotos tectnicos: estn causados por la rotura brusca de las capas rocosas a lo largo
de superficies de fractura (fallas), son los ms fuertes y ms frecuentes.
Terremotos causados por explosiones: las explosiones producidas por el hombre son capa-
ces de generar vibraciones del terreno, con una intensidad tal que pueda causar movimien-
tos en las estructuras.
En general, el movimiento de la corteza se produce por un choque o movimiento brusco
ocurrido a una cierta profundidad bajo la superficie terrestre en un punto terico denomi-
nado foco o hipocentro, a su proyeccin sobre la superficie terrestre se le denomina epi-
centro.

sismgrafo

estructura

epicentro

foco

Fig. 1.1 Definiciones geomtricas de un sismo
1-2

Los sismos desde el punto de vista de la ingeniera y su caracteri-
zacin

Los terremotos ms importantes son los tectnicos, pues son los que traen consecuencias
ms desastrosas en las estructuras que afectan, debido a esto, son los que se tienen en
cuenta para la elaboracin de normas para la contraccin de estructuras sismoresistentes.

La intensidad ssmica es una medida de los efectos de los terremotos en el entorno y en
particular sobre las estructuras.
Existen diferentes escalas de intensidades que describen, para cada valor que esta tome,
los efectos que produce el terremoto. Una de las ms difundidas es la escala de Mercalli
Modificada.
Algunos de los efectos sobre las estructuras en orden creciente de intensidad son:
1. fisuracin de las estructuras de madera
2. agrietamiento de las estructuras dbiles de mampostera
3. agrietamiento de las estructuras ordinarias de mampostera
4. colapso parcial de estructuras ordinarias de mampostera; dao en estructuras bien
ejecutadas de mampostera no diseadas para resistir fuerzas ssmicas
5. colapso de estructuras ordinarias de mampostera; las estructuras con diseo anti-
ssmico son seriamente daadas; daos en cimientos; grietas en el terreno
la mayora de las estructuras son destruidas junto con sus cimientos, daos importantes en
presas y diques, grandes deslizamientos del terreno destruccin casi total, grandes masas
de rocas desplazadas, etc.
Un sismo se caracteriza por su intensidad (parmetro subjetivo) y por su magnitud (pa-
rmetro objetivo).
La escala objetiva ms popular es la de Ritcher, en la que la magnitud M mide la energa
del terremoto en el foco y es el logaritmo decimal de la amplitud del movimiento ssmico,
medido en micrones a 100[km] del epicentro, por un sismgrafo Wood-Anderson estn-
dar. La magnitud M est relacionada con la energa del terremoto, en ergios, por la expre-
sin:
log E = 11,8 + 1,5M

Se han establecido varias relaciones empricas entre la intensidad I
MM
y la magnitud M,
enumeramos algunas a continuacin:

Esteva y Rosenblueth:

I
MM
= 8,16 + 1,45 M 2,46 log R

R: distancia focal en [km]

Tambin se ha relacionado la magnitud M con los valores mximos de las caractersticas
cinemticas del movimiento, estas relaciones se han establecido estadsticamente.
1-3

M
m
m

Donovan:

a
m
=
1080 e
0,5 M

(R + 25)
1,25

Esteva y Villaverde:

a
m
: aceleracin mxima del terreno en [cm/s
2
]
R: distancia focal en [km]

a
m
=

g
5,7 e
0,8 M
(R + 40)
2

v =
32 e

m
(R + 25)
1,7

Esteva:

a
m
: aceleracin mxima en [cm/s
2
]
v
m
: velocidad mxima en [cm/s] R:
distancia focal en [km]

a = 1230 e
0,8 M
(R + 25)
2

v = 15 e
M
(R + 0,17 e
0,59 M
)
1,7

a
m
: aceleracin mxima en [cm/s
2
]
v
m
: velocidad mxima en [cm/s] R:
distancia focal en [km]

Registro de ondas ssmicas. Parmetros utilizados y mapas de
riesgo ssmico

Los terremotos son fenmenos debidos a la brusca liberacin de la energa de deformacin
acumulada durante largos periodos de tiempo en la zona superficial de la tierra. Los sis-
mos producen ondas de varios tipos, que se propagan desde su foco en todas las direccio-
nes a travs de la tierra. Estas ondas son registradas mediante aparatos denominados sis-
mgrafos, diseados para medir la aceleracin, la velocidad o el desplazamiento del mo-
vimiento ssmico. Estos parmetros son relativos, ya que los valores obtenidos estn afec-
tados por las caractersticas del instrumento registrador y por las condiciones de ruido am-
biental en el lugar de registro.

Los mapas de riesgo ssmico representan una sntesis de todos los datos sismolgicos y
geolgicos de un pas. Estos mapas se utilizan para determinar el nivel de proteccin que
se debe alcanzar en las estructuras en cada zona de riesgo.
Diversos aspectos brindan la subdivisin en zonas, pero los fundamentales son:
1-4

Estudios geolgicos y geotcnicos: proporcionan datos de composicin y caractersticas
dinmicas de las rocas y capas de suelo que componen la corteza terrestre.
Estudios sismolgicos: sintetizan los parmetros que caracterizan la sismicidad de la zona:
1. ubicacin de fallas
2. registro de los terremotos que ocurren en la zona
3. mapas de epicentros
4. datos histricos
5. periodos de retorno (intervalo medio de tiempo en que se espera ocurran dos sismos de
igual o mayor intensidad)
6. datos del mecanismo focal
7. correlacin de la sismicidad de la zona analizada con la de la macrozona en la que se
encuentra
Estudios de Ingeniera y Sismologa:
1. anlisis del efecto que han producido sobre las estructuras y las personas los terremotos
ocurridos en el pasado
2. prediccin estadstica de las caractersticas ms probables de la accin ssmica que
se produzca en la zona
Es importante destacar que la geologa local de la zona puede modificar la propagacin de
las ondas ssmicas. Las ondas se reflejan y se refractan cuando en su recorrido aparece una
discontinuidad, por ejemplo una variacin de las caractersticas mecnicas del terreno, ello
produce cambios en la velocidad.
En general, el clculo y la cuantificacin de las acciones ssmicas en la estructuras se rea-
liza en funcin de protocolos, secuencias y definiciones de acciones dadas por normas y
reglamentos. En los captulos siguientes se ofrecen algunas aplicaciones de este tipo.
Captulo 2 Conceptos Bsicos de Dinmica Estructural
2-1

CAPTULO 2 - CONCEPTOS BASICOS DE DINAMICA ES-
TRUCTURAL

En un sentido amplio, un sistema dinmico es aquel cuyas variables experimentan varia-
ciones en el tiempo y, si se conocen las influencias externas que actan sobre el sistema,
podr predecirse el comportamiento de este.

SISTEMA
DINAMICO
variables con variaciones
temporales

influencias externas
sobre el sistema
conocidas estas acciones
externas, permiten
"predecir" el comportamiento
de las variables temporales

En nuestro curso, los sistemas a estudiar sern sistemas estructurales, las variaciones en el
tiempo sern vibraciones producidas por cargas dinmicas.

SISTEMAS
ESTRUCTURALES
vibraciones

cargas dinmicas
ecuaciones diferenciales
que gobiernan el
comportamiento de las
vibraciones

permiten evaluar el
comportamiento de
la estructura frente
a acciones dinmicas

resolucin de las
ecuaciones diferenciales

Definicin de la accin dinmica

Una accin tiene carcter dinmico cuando su variacin con el tiempo es rpida y da ori-
gen a fuerzas de inercia comparables en magnitud con las fuerzas estticas. Algunas fuen-
tes importantes de vibraciones estructurales son:
2-2

- sismos
- viento
- olas y corrientes de agua
- explosiones e impactos
- cargas mviles (vehculos, personas, etc.)

La definicin de estas cargas externas puede distinguirse entre: determinista y no determi-
nista, sta ltima denominada tambin estocstica o aleatoria.
determinista: cuando su variacin temporal es perfectamente conocida
no determinista: cuando alguno o todos sus parmetros son definidos estadsticamente
En nuestro curso trabajaremos con cargas definidas en forma DETERMINISTA.

Respuesta dinmica cualquier magnitud que pueda caracterizar el efecto de una carga
dinmica sobre la estructura

Una carga definida determinsticamente da origen a una respuesta, tambin determinista.

Fig. 2.1- Definicin de la respuesta dinmica: para un punto considerado se calculan:
deformaciones, aceleraciones, tensiones, etc.

Acciones y fuerzas dinmicas

Las acciones dinmicas definidas utilizando representaciones deterministas, son funciones
del tiempo cuyo valor en cada instante ES CONOCIDO.

Este tipo de representacin es apropiado para evaluar el comportamiento de una estructura
A POSTERIORI del acontecimiento que dio lugar a dicha accin. Por ejemplo, evaluar el
comportamiento de un edificio nuevo ante el terremoto ocurrido en Mxico en 1986 (del
que se poseen registros). El diseo de una estructura NO PUEDE encararse en base a ac-
ciones deterministas, pues nada nos asegura que la accin estudiada volver a repetirse.
2-3

F(t)
t3

t1

F

2)

ACCION

M(t1)
M(t2)
M(t3)
t2
F(t

F(t1)
F(t3)

t1 t2 t3 t
este esquema temporal
de carga debe ser
perfectamente conocido

RESPUESTA

Fig. 2.2 - Accin y respuesta determinista

ESTRUCTURA

ACCION
DINAMICA
la comparacin se
basa en: - velocidad de
la accin
- periodo propio
de la estructura

masa M
rigidez K
a(t)
da origen a fuerzas
de "inercia" comparables
con las estticas

depende de
K , M

No considerada como
"carga" sino como una
propiedad intrnseca de
la estructura
Fig. 2.3 - Accin dinmica y propiedades de la estructura

Importancia de la masa en el problema dinmico

Aunque la carga vare con el tiempo, la respuesta de una estructura vara radicalmente se-
gn la masa que vibra con ella. Ante una misma funcin de carga, una estructura SIN
MASA y una CON MASA responden de la siguiente manera:
2-4

a) Estructura SIN MASA SIN INERCIA

F(t)
x(t)

A
A'
F

F0

rigidez: K
masa: m = 0
x
K x(t) = F(t)

La respuesta seguir
exactamente la forma
de la carga
t1 t
x0 = F0/K

t1 t

b) Estructura CON MASA CON INERCIA

F(t) m/2

m/2
x(t)

F
F0

rigidez: K
masa: m
m &x&(t ) + k x(t ) = F (t )
x

t1

xmx
t

x0
vibraciones
remanentes

Se desarrolla energa
cintica, que modifica
la respuesta y deja vibraciones
remanentes

t1 t2 t
t2 > t1 !!

Fig. 2.4 - Importancia de la masa en la respuesta
2-5

Velocidad de reaccin de una estructura

Ante una accin exterior, distintas estructuras reaccionarn de formas diferentes. Esta res-
puesta est ntimamente relacionada con las formas o modos de vibrar y sus correspon-
dientes frecuencias o periodos propios. En el caso de un oscilador de 1 grado de libertad,
este periodo propio se obtiene fcilmente. No as para estructuras de mltiples GLD.
Como veremos en los captulos siguientes, los periodos y formas de vibrar dependen de las
caractersticas geomtricas y de materiales (rigidez) y de la inercia que la estructura opone
al movimiento (masa).

con AMORTIGUAMIENTO
la amplitud decrece en cada ciclo

x
m
x0

ESTRUCTURA SIN
AMORTIGUAMIENTO

k
x0
t

vibraciones
libres

F(t)
ESTRUCTURA CON
T
AMORTIGUAMIENTO

el periodo propio permanece
prcticamente constante
m
F F

k
F0

tD t tD t

En general si tD >> T no es necesario un anlisis dinmico
si tD T PROBLEMA DINAMICO

Fig. 2.5 - Velocidad de reaccin T vs. t
D

2-6

x0
m
x0

a masa constante
k1 > k2

k

T1 T2
2
1

m2 > m1

m ; k ;

variables

m
1

T
1

m
2

T
2
, etc

k
1
k
2

Fig. 2.6 - Velocidad de reaccin; varios T
i

Modelos dinmicos caractersticos

Desde el punto de vista del clculo numrico, obtener la respuesta dinmica de una estruc-
tura, es el resultado de "filtrar" la seal de excitacin a travs de la misma estructura y ob-
tener las variaciones de las magnitudes de anlisis (desplazamientos, velocidades, acelera-
ciones, momentos, tensiones, etc.) respecto del tiempo.
La obtencin de la respuesta requiere, previamente, la definicin del movimiento del te-
rreno (en caso ssmico) tanto como de las caractersticas estructurales del mismo y de la
estructura propiamente dicha. El anlisis es practicado, no a la propia estructura sino a un
modelo mecnico de la misma. La definicin del modelo depende del tipo de estructura
analizado y pretende brindar una serie de relaciones entre acciones y respuesta que descri-
ban un modelo matemtico del problema.
Este modelo matemtico puede ser resuelto mediante diversas tcnicas. En nuestro caso
haremos hincapi en los mtodos numricos de anlisis.
Segn la certeza con que fueron formulados los modelos y procedimientos o algoritmos de
clculo durante el anlisis, ser la precisin de la respuesta obtenida.
2-7

ayed

axed

respuesta en un punto
en particular

onda
refractada
a2
discontinuidad

a1

onda
a3

onda
en el suelo
original reflejada

a1 axed

t

onda ssmica
original
t

respuesta

Fig. 2.7 - Filtrado de una seal ssmica

Se brindan, a continuacin, algunas definiciones tpicas del anlisis estructural dinmico
de una estructura:

Grados de libertad (GL)
Se definen como grados de libertad (GL) a los puntos de la estructura en los cuales se
identifica algn desplazamiento y permiten brindar una deformada de la estructura.

Grados de libertad dinmicos (GLD)
Son los grados de libertad que tienen asociada masa y para los cuales puede conocerse las
vibraciones o movimientos a lo largo del tiempo.
2-8

ESTRUCTURA
REAL
x2

MODELO
DINAMICO
m2
k2

MODELO
MATEMATICO
x1
m1
k1
[M ]&x& + [K ]x = [M ]a

EXCITACION
SISMICA

PROCEDIMIENTOS
NUMERICOS

x1
menos exacta

t
mas exacta

RESPUESTA (desplaz. en el piso 1)

Fig. 2.8- Modelizacin de una estructura

Mtodos de modelizacin dinmica

Pueden distinguirse modelos dinmicos exactos y modelos dinmicos discretos.
En general, para la primera clase, solo pueden resolverse casos muy sencillos y con poca
aplicacin practica, por lo que a lo largo del curso profundizaremos en modelos discretos.
Para estos mtodos modelos discretos, se debe tener en cuenta que la subdivisin en domi-
nios finitos es tanto espacial (discretizacin estructural) como temporal (solucin para ins-
tantes de tiempo determinados).
2-9

Modelo
dinmico
EXACTO

respuesta dinmica
en CADA PUNTO n infinito de
DE LA ESTRUCTURA puntos espaciales

en CADA INSTANTE n imfinito de
DE TIEMPO puntos temporales

Modelo
dinmico
DISCRETO

n FINITO DE PUNTOS
ESPACIALES (GLD)

DISCRETIZACION
ESPACIAL

n FINITO DE INSTANTES
DE TIEMPO donde
se calcula la respuesta
(para cada GLD)
DISCRETIZACION
TEMPORAL

Fig. 2.9 - Modelos dinmicos

Discretizacin espacial de las estructuras

Fundamentalmente, la diferencia con lo visto en otros cursos de anlisis estructural (estti-
co) radica en que en dinmica estructural, cuando hablamos de discretizar espacialmente,
nos referimos a los GLD.
Un modelo dinmico exacto (con infinitos GLD) acarreara ms inconvenientes en la reso-
lucin matemtica que beneficios en su precisin. Adems, en estructuras de edificios y en
la mayora de las estructuras civiles, las masas se encuentran ms o menos concentradas en
lugares conocidos. Es por esto que nuestro principal mtodo de modelizacin dinmica
ser el de las MASAS CONCENTRADAS.
No obstante, existen otros, como ser:

- mtodo de los DESPLAZAMIENTOS GENERALIZADOS
- mtodo de los ELEMENTOS FINITOS
2-10

Mtodo de las masas concentradas

n de GLD
del modelo
n total de componentes de
desplazamiento segn los
cuales vibran las masas
concentradas

se distinguen:
- MODELOS DE 1 GLD
- MODELOS DE
MULTIPLES GLD
se calcula la "deformada"
del modelo en cada instante

tensiones y deformaciones
especficas pueden conocerse
mediante los procedimientos
del anlisis esttico

Modelos con 1 GLD:

m m m
x
c c
k k k
a(t)

(a) (b) (c)

Fig. 2.10 - Modelos con un solo grado de libertad. (a) modelo conservativo; (b) modelo
con amortiguamiento; (c) modelo ssmico.

m
x

k
a(t)

(a) (b) (c)

Fig. 2.11 - Estructuras modelizadas como un sistema de un solo grado de libertad. (a)
prtico; (b) el mismo prtico con la masa concentrada al nivel de la viga; (c) modelo di-
nmico.
2-11

Modelos con mltiples GLD:

mn mn
kn

kn
cn
mn
kn
cn

mr
xr
kr

mr
kr
cr

xr
mr
xr
kr
cr

m2 m2
k2
m1 m1
k1
k2
c1
k1
c1
m2

m1
a(t)
k2
c1
k1
c1

(a) (b) (c)

Fig. 2.12 -Modelos con varios grados de libertad. (a) modelo conservativo; (b) modelo
con amortiguamiento; (c) modelo ssmico.

Fig. 2.13 - Estructura con dos grados de libertad: Prtico de dos pisos y su modelo din-
mico.

a(t)

Fig. 2.14 - Estructura con masa distribuida (antena) y su modelo dinmico discreto con n
grados de libertad.
2-12

Fig. 2.15 - Modelo dinmico de un prtico de cortante y prtico espacial modelizado co-
mo un sistema completo (10 grados de libertad) y simplificado ( dos grados de libertad).

Fig. 2.16 - Modelo dinmico con grados de libertas de rotacin.
2-13

Ecuaciones de movimiento

Las ecuaciones de movimiento son las expresiones matemticas que gobiernan la respues-
ta dinmica de las estructuras. Pueden obtenerse a partir de cualquiera de los principios de
la mecnica clsica:

Principio de Hamilton

t
2
t
2

H
=

(E
p
E
c
) dt +

E
d
dt

; d
H
= 0 (2.1)
t
1
t
1

La primera expresin se denomina funcional de Hamilton; donde E
p
es la energa poten-
cial, E
c
es la energa cintica y E
d
la correspondiente a fuerzas no conservativas. La segun-
da expresin permite establecer el equilibrio a travs de una variacin funcional nula.

Principio de los trabajos virtuales

Se trabaja en forma similar a lo visto en anlisis esttico pero incluyendo las fuerzas de
inercia y disipativas.

Principio de DAlembert

w
i
= w
e

(2.2)

Proporciona el mtodo ms directo para obtener las ecuaciones de movimiento de un sis-
tema dinmico.
Puede formularse como sigue: un sistema dinmico esta en equilibrio cuando todas las
fuerzas que actan en el mismo, incluidas las de inercia y disipativas, cumplen las ecua-
ciones de equilibrio esttico en cada instante de tiempo.

Formulacin de la ecuacin de movimiento para un sistema de 1
GLD

Tomando el sistema de la figura 2-10, podemos distinguir dos casos:

Fig. 2.17- (a) fuerza aplicada; (b) modelo ssmico

Para el modelo (a), aplicando el principio de DAlembert, tendramos:
2-14

e

F(t)

x(t)
&x&(t )

x&(t )
Fi(t) Fa(t) Fe(t)

inercia amort. elstica

F(t) = Fi(t) + Fa(t) + Fe(t) equilibrio

Fig. 2.18 - Equilibrio de fuerzas para 1 GLD

Al aplicar una fuerza exterior F(t), se genera aceleracin, velocidad y desplazamiento para
un cierto instante t; a causa de esto se producen fuerzas:
i- de inercia
F
i
(t ) = m
&
x
&
(t )

(2.3)

ii- de amortiguamiento

F
a
(t ) = c x
&
(t )

(2.4)

iii- elsticas

F (t ) = k x(t )

(2.5)

equilibrio en el instante t

F (t ) F (t ) F

(t ) F

(t ) = 0

(2.6)

i
F (t ) + F

a
(t ) + F

e
(t ) = F (t )

(2.7)

i a e

m
&
x
&
+ c x
&
+ k x = F

(2.8)

Se omite (por simplicidad de notacin) la dependencia del tiempo, pero de aqu en adelan-
te sta se encontrar implcita en toda variable temporal.

ej : x(t ) x ;

a(t ) a ;

x
&
(t ) x
&

F (t ) F

;
&
x
&
(t )
&
x
&

; etc.
2-15

i a e
La ecuacin (2.8) es la de movimiento correspondiente a 1 GLD con carga exterior y
amortiguamiento.
Para el modelo (b) de la figura 2-17, el planteo es similar, solo que no tiene fuerza exterior
aplicada y la fuerza de inercia se ve afectada por la aceleracin total de la masa:

F
i
= m [a +
&
x
&
] (2.9)

entonces, la ecuacin de movimiento queda:

m [a +
&
x
&
] c x
&
k x = 0

m
&
x
&
+ c x
&
+ k x = m a
(2.10)

(2.11)

La (2.11) es la ecuacin de movimiento para 1 GLD con aceleracin de apoyo (ssmico) y
amortiguamiento. Un caso general sera la inclusin de aceleracin de apoyo y fuerza exte-
rior:

m
&
x
&
+ c x
&
+ k x = F m a
(2.12)

Formulacin de las ecuaciones de movimiento para modelos con
mltiples GLD

El modelo con varios grados de libertad ms sencillo es el de edificios de cortante (fig. 2-
12). Est basado en que las plantas son infinitamente rgidas y en que los nicos movi-
mientos posibles de stas son los desplazamientos horizontales.
Aplicando el principio de DAlembert en cada GLD (uno por piso) se obtiene:

F
r
F
ir

F
ar

F
er
= 0
(2.13)

para el piso (r).
Planteando el equilibrio para todos los GLD, nos queda un sistema de ecuaciones (vecto-
rial)

F
T
F F F = 0

(2.14)
2-16

1
F
F
e
2
a
#Nota: en adelante, los vectores y matrices sern representados(en general) con mins-
culas y maysculas en negrita, respectivamente.

F
T
= [ f , f ,... f
n
]
vector de fuerzas externas
i
= M {&x& + J a}

vector de fuerzas de inercia
= C x&

vector de fuerzas disipativas
F = K x

vector de fuerzas elsticas

entonces, el sistema (2.14) puede escribirse:

M matriz de masas
M &x& + C x& + K
x = F M {J a}

(2.15)
C matriz de amortiguamiento
K matriz de rigidez
J
T
= [1,1,1] vector con todos sus elementos igual a uno

Si bien la (2.15) fue deducida para un modelo de edificio cortante, es una expresin AB-
SOLUTAMENTE GENERAL, inclusive para modelos de elementos finitos, y solo varan
las formas de M, C y K.
Para un modelo ssmico, la (2.15) se reduce a:

M &x& + C x& + K x = M {J a}

(2.16)

Y en caso general de prticos 3D, o modelos de elementos finitos, suele sustituirse x por D
para indicar que cada GLD puede sufrir desplazamientos en 3 direcciones y respectivos
giros.

M D
&&
+ C D
&
+ K D = F M {J a} (2.17)

Nota: en este ltimo caso, J llevar unos en las componentes a las cuales se quiera aplicar
el acelerograma a, por ejemplo componente x:
J = [1,0,0,0,0,0,1,0,0,0,0,0,..,1,0,0,0,0,0,..,1,0,0,0,0,0] etc.
Captulo 3 Respuesta de un Oscilador Simple
CAPTULO 3 - RESPUESTA DE UN OSCILADOR SIMPLE

3-1

Ecuacin de movimiento y equilibrio dinmico

Las caractersticas dinmicas de un oscilador de 1 GLD pueden estudiarse mediante un
modelo no amortiguado con vibraciones libres, cuya ecuacin de movimiento es

m
&
x
&
+ k x = 0
(3.1)

Fig. 3.1 - Modelo de 1 GLD, no amortiguado.

La vibracin, del modelo de la fig. 3.1, es inducida por algunas condiciones iniciales, sean
desplazamiento, velocidad o aceleracin en el instante t = 0.
Luego, durante las vibraciones no recibe ningn tipo de perturbacin.
Dividiendo (3.1) por m y usando la notacin:

se obtiene:

2
=
k
m

; =
k
m

(3.2)

&
x
&
+
2
x = 0

(3.3)

es la pulsacin o frecuencia circular o simplemente frecuencia del modelo estudia-
do. Viene expresada en radianes por segundo (1/s).
La frecuencia cclica viene dada por

f = (3.4)
2

y se expresa en ciclos por segundo o hertz.
3-2

x

0
0
2
Finalmente, otra caracterstica es el perodo natural

T =
1

f

(3.5)

T =
2

(3.6)

La solucin general de la (3.1) o (3.3) puede escribirse:

x = A sen( t + )

(3.7)

donde A es la amplitud del movimiento y el ngulo de fase. Los valores de A y se cal-
culan a partir de las condiciones iniciales del problema, por ejemplo para

resulta:
x(0) = x ;
x
&
(0) = x
&
0

A =
2

+

x
&
0

Frmula de Geiger

Sustituyendo

x
tan =
0

x&
0

m =
G

g

G: peso de m

g: aceleracin de la gra-
vedad

=
k
=

G g

G

g
1

G k

entonces:
= X
SG
: desplazamiento esttico producido
k
por el peso G en la direccin del grado de libertad

= g
1
X
SG

(3.9)
3-3

T =
2

g

X
SG

(3.10)

Utilizando unidades de S.I. la (3.10) queda:

T = 2,00

X
SG
(3.11)

con X
SG
expresado en metros para un peso G en Newton.

Caractersticas dinmicas con amortiguamiento

El amortiguamiento puede definirse estudiando las vibraciones libres del modelo de la fi-
gura 3-2:

Fig. 3.2 - Modelo de 1 GLD con amortiguamiento (vibraciones libres)

Si se toma la ecuacin (2.12) sin cargas ni aceleraciones de apoyo (vibraciones libres) y se
divide por m se obtiene:

&x& + 2 x& +
2
x = 0 (3.12)

2 =
c
m

(3.13)

la solucin de (3.12) est dada en la forma:

x = e
rt

(3.14)

que proporciona la ecuacin caracterstica:
3-4

v
r
2
+ 2 r +
2
= 0 (3.15)

Ya que (3.12) es una ecuacin diferencial de segundo orden, lineal, homognea a coefi-
cientes constantes.
Las soluciones de (3.15) son

r
1, 2
=
2

2
= 0
(3.16)

Segn sea el radicando
2
-
2
se encuentran tres tipos de amortiguamiento:

2
-
2
> 0 SUPERCRITICO: la estructura NO VIBRA
2
-
2
= 0 CRTICO: caso lmite = c
r

= 2 m
2
-
2
< 0 SUBCRITICO: la estructura VIBRA con amplitud decreciente
Este es el caso ms frecuente en ingeniera civil, por lo que enfatizaremos su
estudio. Para este caso (subcrtico), la cantidad (
2
-
2
) es negativa, lo que
hace que (3.16) tenga races complejas:

r
1, 2
=

i

1
2
= 0

con i =
(3.17)

1

Llamando frecuencia de vibracin amortiguada a:

se obtiene:
= 1
2

(3.18)

r
1, 2

r
1, 2
= i
v

= i
v

(3.19)
(3.20)

En las ecuaciones anteriores aparece la magnitud

=
c
c
r

=
2 m

2 m

(3.21)
3-5

1
v
conocida como fraccin de amortiguamiento crtico (en estructuras corrientes
0.02<<0.06), tambin:

=

o =

c
2m
Volviendo a la resolucin de la ecuacin de movimiento (3.12), escribimos la solucin ge-
neral en la forma:

x = c e
r
1
t

+ c
2

e
r
2
t
(3.22)

sustituyendo r
1
y r
2
por la expresin 3.20, se obtiene

x = A e
t
sen( t + )

(3.23)

Las constantes A y se obtienen de las condiciones iniciales del problema.

Fig. 3.3 - Vibraciones libres amortiguadas / 1 GLD

La evaluacin del amortiguamiento en una estructura es un problema esencial en la din-
mica estructural.
El origen de las fuerzas de amortiguamiento se debe a diferentes causas:
- Rozamiento entre superficies de deslizamiento, en donde la fuerza de amortiguamiento
es proporcional a la fuerza normal y al coeficiente de rozamiento (hiptesis de Cou-
lomb)
- Amortiguamiento debido a fuerzas aero o hidrodinmicas
- Debido a friccin interna del material de la estructura
Generalmente, en el clculo dinmico de estructuras, se utiliza un modelo de gran simpli-
cidad que caracteriza el amortiguamiento de toda la estructura. Este modelo denominado
de amortiguamiento viscoso se debe a Kelvin-Voigt y es proporcional a la velocidad.
F
a
= c x&
3-6

n

Los reglamentos brindan los coeficientes de amortiguamiento para cada tipo de estructura,
pero puede obtenerse en forma experimental y con un mtodo relativamente simple:

Determinacin prctica de
Para amortiguamientos bajos (del orden del 10% de crtico) la relacin entre dos picos su-
cesivos puede aproximarse:

x (n)

Ae
t
n
sen( t

+ )

mx

v n

x
mx
(n + 1)
=
A e
t
n+1
sen(
v
t
n+1

+ )
(3.24)

pero t
n +1

= t
n
+ T
v

; con T
v
= 2
v

t
n +1
= t
+
2

, entonces
v

x (n)

e
t
n
sen( t

+ )

mx
=
v n
x
mx
(n + 1)

t
n
+

2

v

2

(3.25)

e

sen
v
t
n
+
v

v
+

x (n)

e
t
n
sen( t

+ )

mx

v n

x
mx
(n + 1)
=

e
t
n
e

2
v
sen(

v
t
n

+ 2

+ )

x
mx
(n)
=

x
mx
(n + 1)

1
e
2
v

(3.26)

tomando logaritmo natural:

ln

x
mx
(n)

= ln[e
2
v
]

(3.27)

x
mx
(n + 1)

# notar que para amortiguamientos del orden de

3-7

= 0,1
v
= 1
2

Entonces
v
= 0,995
v

x (n)

ln

mx
(

1)

=
2

(3.28)

x
mx
n +

Para el caso de lecturas separadas por N ciclos:

ln[x

=

mx
(n)

x
mx
(n + N )]

(3.29)

2N

Excitacin peridica

En la figura 3-4 pueden observarse diversas funciones de carga. De stas, nos interesan por
ahora, las peridicas y ms particularmente las excitaciones armnicas ya que mediante
series de Fourier cualquier excitacin peridica puede llevarse a una suma de armnicas
simples.

Fig. 3.4 - Tipos de cargas dinmicas. (a) armnica; (b) peridicas; (c) cuasi peridicas;
(d), (e) fuerzas impulsivas; (f) carga dinmica general; (g) aceleracin ssmica del terre-
no.
3-8

h 1 2

Excitacin armnica

Si la carga es de tipo:

P = P
0

sen(t )
Entonces, la ecuacin de movimiento ser:

m &x& + c x& + k x = P
0
sen(t )

(3.30)

: frecuencia de la excitacin

la (3.30) puede escribirse tambin:

&x& + 2 x& +
2
x =
P
0
sen(t )
m

(3.31)

La solucin general de esta ecuacin viene dada por

x
g
= x
h

+ x
p
(3.32)

x = c e
r
1
t

+ c e
r
2
t

(3.22)

Es la solucin de la ecuacin diferencial homognea.

x
h
= c
1

e
t
e
i
v
t
+ c
2

e
t
e
i
v
t
(3.33)

Utilizando matemtica para nmeros complejos, esta ltima ecuacin puede escribirse:

x
h
= e

t
(c
1
' sen
v
t + c
2
' cos
v
t )
(3.34)

x
p
en la (3.32) es la solucin particular y lleva la forma

x
p
= A sen t + B cos t
(3.35)
3-9

Derivando y reemplazando en (3.30) se obtienen las constantes A y B:

0
g 0 2

Denotando =

P

1
2

(3.36)
A =
0

k

(1
2
)
2
+ (2 )
2

(3.37)

B =
P
0

k
2

(1
2
)
2
+ (2 )
2

(3.38)

Condiciones iniciales

Para el caso en que

x(0) = x ;

x&(0) = x&
0

Pueden calcularse las partes correspondientes a la solucin particular (para t = 0)

x =
P
0

2

p (0 )
k
(1
2
)
2

+ (2 )
2

(3.39)

P
0

x&
p (0 )
=

k
(1
2
)

(1
2
)
2
+ (2 )
2

(3.40)

Basados en stas y en (3.32) y (3.34) podemos plantear las siguientes ecuaciones:

x (0) = x
= c '+ x
p 0

(3.41)
x&
g
(0) = x&
0
= x
0
+ x
p 0
+
v
c
1
'+ x&
p 0

(3.42)

y despejarse c
1
y c
2
:

1

c
1
' =

v
[x&
0

+ x
0

x&
p 0

x p 0
]

(3.43)
c
2
' = x
0
x
p 0

(3.44)
3-10

f
Forma de operar:

1- Se calculan

x
p 0
y

x&
p 0

con las 3.39 y 3.40

2- Se computan c
1
y c
2
con las 3.43 y 3.44

3- Se reemplaza todo en 3.32 y se tiene la respuesta en desplazamientos
4- Para obtener velocidades y aceleraciones se deriva la 3.32

Nota: la parte de la solucin correspondiente a la ecuacin homognea incluye el coefi-
ciente e
-t
que es una funcin decreciente con el tiempo, por lo que se la denomina transi-
toria y desaparece antes o despus segn sea el valor de .
La estructura sigue vibrando con una frecuencia prcticamente igual a la de la excitacin
. Esta parte de la solucin se denomina respuesta en rgimen.

Excitacin arbitraria. Integral de Duhamel

Respuesta a un impulso elemental

Dada una carga impulsiva como la de la fig. 3-4(d), para condiciones iniciales nulas, pue-
de calcularse la velocidad y posicin al finalizar el impulso y a partir de ese instante se
trata de un problema de vibraciones libres.

P
P()
m x
c

x& = tg

i f

d

t
F(t)
f

i
d t

x&
f

= x&
i
+ &x&d

Fig. 3.5 - Respuesta a un impulso
; &x& = P( ) m

entonces:

x&
f
=
P( )

d

m

(a)
3-11

v
m
t
v
v
por otra parte

x
f
= x
i

x

+ x& d +
1
&x&(d )
2

2
=
1
&x&(d )
2

(b)

(b)
f
2

Si d es muy pequeo, entonces (b) es despreciable frente a (a), es decir, luego de actuar
un impulso la masa queda con velocidad pero prcticamente no se desplaza. La velocidad
dada por (a) es la condicin inicial para las vibraciones libres siguientes.

Reemplazando sta en la (3.23), para un tiempo infinitesimal, se obtiene:

dx(t ) = e

( )
P( )

d

sen(

(t
))

(3.45)
t

m
v

Si el sistema es lineal, pueden superponerse las respuestas de una sucesin de impulsos
infinitesimales hasta el tiempo genrico t. La solucin dada por esta superposicin es
conocida como Integral de Duhamel

x(t ) =
1

P( )e
(t )
sen[

v 0
(t )]d

(3.46)

En el caso de excitacin ssmica, se transforma en:

x(t ) =
1

t

a( )e
(t )
sen[

(t )]d

(3.47)
v 0

Para el caso en que el problema tenga condiciones iniciales x&
0
0 y

x
0
0 , la solucin se

obtiene con el mismo procedimiento dado para las cargas armnicas, solo que la solucin
particular vendr dada por (3.46) o (3.47).

Soluciones a la integral de Duhamel:

3-12

(a) Para ciertas funciones de P() se encuentran tabuladas las primitivas de (3.46) y (3.47),
as como el mximo factor dinmico (FAD) y el instante de mxima respuesta.
3-13

S S =
v
r
(b) Para funciones aproximadas por trozos rectos, tambin se pueden encontrar solucio-
nes exactas a dicha integral (punto 5.3.4 / pg. 154 / ref. [1])

(c) En el caso mas general se puede integrar numricamente utilizando alguno de los
mtodos conocidos (Trapecios, Simpson, etc.)

Factor de amplificacin dinmica

Es el coeficiente entre el mximo desplazamiento dinmico y el que producira la carga en
forma esttica:

FAD =

x(t )

mx

(3.48)

x
est

Espectros ssmicos de respuesta

Puede definirse el espectro d respuesta (para un acelerograma dado) como los mximos
valores de la respuesta de un sistema, expresado en funcin del periodo propio. Esto es
muy til para el diseo de estructuras donde solo los mximos son necesarios.

Para un sistema de 1 GLD:

m &x& + c x& + k x = m a(t )

(2.11)

La solucin a esta ecuacin viene dada por la ec. (3.47) y el mximo valor absoluto ser
S
d
.

Reemplazando S
d
= x(t )

mx

Puede obtenerse, para 1 GLD, la expresin exacta del espectro de desplazamientos. Deri-
vando una y dos veces se obtienen las expresiones para velocidad y aceleracin:

(Ver ref. [1], cap. 5)
r
x&(t )

mx

;
a
=
&x&(t )

mx
3-14

t

Fig. 3.6 - Espectro ssmico de respuesta para desplazamientos. S
d

No obstante pueden hacerse ciertas hiptesis que simplifican las expresiones obtenidas.
Estas nuevas ecuaciones definen los Seudo espectros de respuesta.

Hiptesis simplificativas:
i- en estructuras civiles, 0,002 < < 0,2, por lo que puede reemplazarse
v
por
ii- la funcin coseno que aparece en el espectro de velocidades (al derivar 3.47) puede
reemplazarse, sin que ello implique grandes variaciones en los valores mximos, por
una tipo seno
entonces:

s (, ) =

1
a( )e
(t )
sen[ (t )]d

(3.49)

d
0
mx

s
v
(, ) =

t

a( )e
(t )
sen[ (t )]d

0

mx

(3.50)
3-15

2

s
a
(, ) =
t

a( )e
(t )
sen[ (t )]d

0

mx

(3.51)

definen los seudo espectros para desplazamientos, velocidades y aceleraciones.

Es importante notar que:

s
v
= S
d
(3.52)
s
a
= S
d

(3.53)

Por lo que un espectro completo puede representarse en un grfico tipo logartmico como
el de la figura 3-7.

Fig. 3.7 - Comparacin entre los espectros suavizados de Newmark y Hall, los espectros
de diseo del UBC y los espectros elsticos para los registros de EL Centro y de la presa
de Pacoima.

La normativa argentina INPRES-CIRSOC 103 establece la accin ssmica como una fun-
cin del seudo espectro de respuesta de aceleraciones S
a
. Es por eso que solo aparecen gr-
ficos del tipo de la fig. 3-9, que varan segn el tipo de suelo y la zona ssmica considera-
da.
3-16

Fig. 3.8 - Espectro de diseo suavizado de Seed e Idriss (1982).

Fig. 3.9 - Seudo espectro y espectro de aceleraciones; (a) corresponde al acelerograma
(1). (b) corresponde al acelerograma (2).
3-17

Integracin numrica de la ecuacin de movimiento

Se trata de resolver en forma discreta la ecuacin de movimiento para 1 GLD dada por
(2.8) y (2.11) utilizando ecuaciones en diferencias que permitan obtener x , x& y
instante t+t en funcin, nicamente, de los valores en el instante t.
&x& en un
Existen varios mtodos para plantear las ecuaciones en diferencias antes mencionadas, pe-
ro solo nos dedicaremos en profundidad al de Newmark, dado que es uno de los ms di-
fundidos en software de anlisis dinmico mediante FEM.

Mtodo de Newmark (1 GLD)

m &x& + c x& + k x = F

(3.54)

con F = f (t ) + [ m a(t )]

&x&
Fig. - 3.10

&x&
i +1

t = t
i +1
t
i

haciendo el cambio de variable

&x&
i

= t t
i

; t

t

ti ti+1 t

puede suponerse que el valor de
la aceleracin de respuesta en
un instante se expresa como:

&x&( ) = &x&
i
+ f ( )(&x&
i+1
&x&
i
)

(3.55)

f() =
0 para = 0

1 para = t

integrando (3.55) se obtiene la expresin de la velocidad

x&( ) = x&
i
+ &x&
i
+ (&x&
i+1
&x&
i
)

f ( )d

0

(3.56)
3-18

g ( ) =

0
f ( )d

(3.57)

t =

0

f ( )d

(3.58)

La expresin (3.56) queda:

x&( ) = x&
i
+ &x&
i
+ (&x&
i+1
&x&
i
)g ( )

(3.59)

que para =t resulta

x&
i+1
= x&
i
+ &x&
i
t + (&x&
i+1
&x&
i
) t

x&
i+1
= x&
i
+ [(1 )&x&
i
+ &x&
i+1
]t

(3.60)

(3.61)

Para calcular los desplazamientos, se integra (3.59) y se obtiene

2
x( ) = x
i
+ x&
i
+ &x&
i

2

+ (&x&
i+1
&x&
i
)

g ( )d

0

(3.62)

particularizando para =t y llamando

t
2

=

g ( )d

0

(3.63)

se obtiene la relacin final en diferencias, propuesta por Newmark

x = x

+ x& +

1

&x&

+ &x&

t
2

(3.64)

i+1 i

i

2

i

i+1

Las ecuaciones (3.61) y (3.64) juntamente con la ecuacin de movimiento (3.54) permiten
obtener aceleraciones, velocidades y desplazamientos en un instante t+t con solo conocer
3-19

estas magnitudes en el instante t.
3-20

0 1
2

6
Los parmetros y surgen de un anlisis de estabilidad numrica que escapa a este curso

(ver. ref. [1], [2] y [4]) pero pueden tomarse:
= 0,5
= 0,25

Desarrollo y forma operativa

La ecuacin de movimiento (3.54) valuada en t=t
i+1
toma la forma:

m &x&
i+1
+ c x&
i+1
+ k x
i+1
= F i+1
(3.65)

Reemplazando en sta los valores en diferencias dados por (3.61) y (3.64), agrupando los
valores que permanecen constantes en cantidades denominadas por a
0
,, a
7
; se obtiene el
algoritmo de integracin paso a paso.

A) clculos iniciales

1- determinar las propiedades de la estructura: m, c, k
2- inicializar

x
0
,

x&
0
y

&x&
0

utilizando, de ser necesario, la ecuacin (3.54)
3- seleccionar el paso de tiempo y los parmetros y
t ; ; t T/10 ; =0,5 ; =0,25
4- calcular las constantes:

a
0
=
1
;

t
2

a
1
=

;
t

a
2
=

1
t

a =
1

1 ;

a =

1 ;

a =
t

2

3
2
4

5

a = t(1 ) ; a
7
= t

5- formar la rigidez efectiva

B) para cada paso de integracin

k

= k + a m + a c
1- formar el trmino de carga efectiva en t+t
r
i+1
= F i+1 + m(a
0
x
i
+ a
2
x&
i
+ a
3
&x&
i
) + c(a
1
x
i
+ a
4
x&
i
+ a
5
&x&
i
)

2- resolver el desplazamiento en t+t
3-21

k

x
i+1
= r
1+1
3- calcular aceleraciones y velocidades en t+t
&x&
i+1
= a
0
(x
+1i
x
i
) a
2
x&
i
a
3
&x&
i
x&
i+1
= x&
i
a
6
&x&
i
+ a
7
&x&
i+1

Captulo 4 Respuesta Dinmica de una Estructura con Mltiples Grados de Libertad
4-1

CAPTULO 4 - RESPUESTA DINMICA DE UNA ESTRUCTURA
CON MLTIPLES GRADOS DE LIBERTAD

Ecuaciones de movimiento y equilibrio dinmico

Como se vio en el Captulo 2 un modelo estructural dinmico esta dado por
- discretizacin espacial (por ej. masas concentradas)
- discretizacin temporal (por ej. Newmark)

En el Captulo 3 se estudio el modelo discreto ms sencillo (1 GLD) mientras que en la
UT2 se plantearon (sin resolver) las ecuaciones de equilibrio dinmico para sistemas de
mltiples GLD, utilizando el principio de DAlembert.

En general, el sistema de ecuaciones diferenciales es del tipo

M D
&&
+ C D
&
+ K D = F (t ) (4.1)

En particular, para el caso ssmico:

M D
&&
+ C D
&
+ K D = M {J a} (4.2)

Recordar que J es un vector con 1 en la posicin del GLD en el que acta la aceleracin de
apoyo a(t).

Vibraciones libres

Caractersticas dinmicas

El sistema que gobierna las vibraciones libres en un sistema de mltiples GLD es

M D
&&
+ C D
&
+ K D = 0 (4.3)

Como ya se vio, el cambio de frecuencia propia debido a considerar o no el amortigua-
miento, no es relevante para estructuras civiles corrientes, por lo que seguiremos el anli-
sis con el sistema simplificado:
4-2

1 2
e
M D
&&
+ K D = 0 (4.4)

Una solucin a este sistema puede ser de la forma

D = A sen(e t + + )
A: vector que contiene las amplitudes de las vibraciones
+: ngulo de fase

(4.5)

D
&&
(t ) = e
2
A sen(e t + + ) (4.6)

reemplazada en (4.4) queda:

M e
2
A sen(e t + + ) + K A sen(e t + + ) = 0
{K e
2
M }A sen(e t + + ) = 0

(4.7)

para que haya vibraciones, e=0, por lo que

sen(e t + + ) = 0

y puede eliminarse de (4.7), quedando finalmente

{K e
2
M }A = 0

(4.8)

Nos interesan las soluciones de A distintas de la trivial. La ecuacin (su resolucin) 4.8
representa un problema de autovectores y autovalores, en donde:

K e
2
M = 0
(4.9)

Este determinante puede desarrollarse en la forma polinmica:
e
2 n

+ o e
2 n 2

+ o e
2 n 4

+ ...o
2
n 1

+ o
n
= 0

(4.10)

obtenindose la ecuacin caracterstica.
4-3

A = a

a
=
A
i

i i

Cuando K y M son definidas positivas (caso usual en estructuras civiles), de la ecuacin
caracterstica se obtienen n soluciones positivas e
i
2
y en consecuencia n valores e
i
,
siendo n el nmero de GLD de la estructura.
Los valores e
i
se denominan frecuencias propias o pulsaciones de la estructura y los
n periodos propios se calculan

T =
2t

i
e

i = 1,...n

(4.11)
i

Siendo T
1
el perodo correspondiente e
1
, que es la frecuencia de menor valor, ste se de-
nomina periodo fundamental del sistema. Reemplazando cada e
i
en la (4.8) se obtiene
el correspondiente A
i
que se denomina vector de forma modal o simplemente modo.
Este vector modo contiene la forma que tomara la estructura en cada vibracin. Para ca-
da frecuencia vibratoria, la forma modal ser diferente. Es importante notar que (4.8) se
cumple para cualquier mltiplo de A
i
, por lo tanto inferimos que no interesa la magnitud
de las componentes de A
i
sino la relacin existente entre ellas. Por esto es conveniente
normalizar los distintos A
i
.

Normalizacin de los modos

1. Por la mxima componente

a
i1

M

puede dividirse
dado

i imx
cada elemento por

M

a
in

el a
imx

i1

i
i
imx

(4.12)

M

i
=
1

M

modo i normalizado

2. Segn la matriz de masas

= A (M

*)

1
2

in

(4.13)
4-4

T

T
Con

M
i
* = A
i

|M |A
i

(4.14)

Esta normalizacin permite asegurar la condicin

i
|M |
i
= 1
(4.15)

Expresin que ser de suma importancia ms adelante.

Obtencin de los grados de libertad dinmicos

No siempre todos los grados de libertad estticos que definen el comportamiento de la es-
tructura tienen asociada masa, es decir que no son necesarios en el anlisis dinmico.
Por ejemplo los prticos de las figuras 2-15 (f) y (h) que son analizados con modelos di-
nmicos de menor orden que los respectivos estticos (GLE). Lo mismo se ve en la figura
2-16. La forma ms simple de reducir en nmero de grados de libertad sin perder precisin
es mediante la:

Condensacin esttica de la matriz de rigidez

Dada una estructura y su modelo esttico puede subdividirse la matriz de rigidez (as como
el vector desplazamientos) de manera de separar las ecuaciones que tienen asociada masa
(GLD) y las que no la tienen (GLE).

d3 d4 d4

kv
m2 m
k2 k2
d1 d2 d2

kv
m1 m
k1 k1

modelo esttico
4 x 4
modelo dinmico
2 x 2

Fig. 4.1 - Condensacin esttica

Para el ejemplo de la figura 4-1, tenemos que, observando el vector desplazamientos est-
ticos D
s
podemos determinar que parte queremos eliminar y cuales grados de libertad de-
seamos conservar como GLD.
4-5

D

s
K D
e
21 22 e
D
K D = F
e
K
22
D
21
K
11 12
K
{K -K K K
D

D

D: vector con los GLD

s
=

e

D
e
: grados de libertad a eliminar

s s

(4.1)

El sistema de ecuaciones esttico puede escribirse

K K

D

F

11 12

=

(4.16)

K
21
22
e

F
e

Como F
e
representa las fuerzas en los grados de libertad a eliminar

F = 0
(4.17)

La segunda ecuacin (vectorial) de (4.16) puede, entonces, expresarse como

K D + K D = 0
(4.18)

de donde

= K
1

(4.19)

Desarrollando ahora la primera ecuacin de (4.16) tendremos

D K {K
1
22 21
D}= F

(4.20)

reordenando

11 12

1
22 21

}D = F

(4.21)
K D = F

donde
4-6

K es la matriz de rigidez condensada y vale

(4.22)
4-7

-K K K
11
o
21
K = K

11 12
1
22 21

(4.23)

Finalmente el sistema dinmico de la fig. 4-1 estar representado matemticamente por

K D + M D
&&
= 0

+
22
+
2
2

d

D =

2

(4.24)
d
4

Otra forma de obtener la matriz de rigidez del modelo dinmico es invirtiendo la matriz de
flexibilidad de los GLE que tienen asociada masa.
Este mtodo, si bien menos formal, puede ser de ms sencilla aplicacin para procedi-
mientos manuales. Para nuestro ejemplo:

o
m2

o

m1
1

o

1
o

F =
o o

o

12 22

K = F
1

Fig. 4-2 - Obtencin de K mediante F
-1

4-8

Matriz de amortiguamiento

Dado el sistema
4-9

1 2
= 0
T
M D
&&
+ C D
&
+ K D = F (t )

(4.1)

solo nos resta conocer (para encarar la resolucin) la matriz C, puesto que en el apartado
anterior se mostraron dos maneras de obtener K y sabemos (porque usamos un modelo
discreto de masas concentradas) que:

m
1
L 0

O

M = M

m
i
M

O

(4.25)
0 L m
n

Conviene definir matrices de amortiguamiento ortogonales, pero para esto debemos, pri-
mero incursionar un poco en las condiciones de ortogonalidad de los modos propios.

Condiciones de ortogonalidad
Las frecuencias y modos propios pueden ordenarse en matrices denominadas espectrales y
modales, que son, respectivamente:

1

O =

O
0

; u = |
,

K
n
|

(4.26)
0
N

Puede demostrarse (ver ref. [1], [2], [3], [4]) que

i j
i = j (4.27)

Esta condicin para los modos puede extenderse a

4-
10

j i

T
M = 0 i = j
(4.28)
4-
11

j

j
i
i
T
1 2
i
i
T
K = 0 i = j
(4.29)

Si los vectores
i
fueron normalizados segn (4.13) las condiciones de ortogonalidad y
normalizacin pueden expresarse como una nica condicin de ortonormalidad

u
T
M u = I

u
T
K u = K*
(4.30)

(4.31)

donde I es la matriz identidad y K* es una matriz diagonal

k
1
*

O

K *=

k
i
*
O

k
n
*

; k
i
* =
i

K
i

(4.32)

Matrices de amortiguamiento ortogonales

Una de las hiptesis para lograr una representacin numrica del amortiguamiento de una
estructura est dada por suponer que existe un mecanismo de disipacin uniforme de ener-
ga. En tal caso puede desarrollarse una matriz de amortiguamiento que cumpla las condi-
ciones de ortogonalidad respecto a la matriz modal.
El amortiguamiento proporcional permite definir una matriz que sea proporcional a la de
las masas, a la de rigidez o a ambas.

C = o M + o K (4.33)

La condicin de ortogonalidad:

T
C = 0 i = j
(4.34)

T
C = 2
i
i
(4.35)
4-
12

e
i i
o
1
+ o
i
Los coeficientes o
1
y o
2
se calculan a partir de las ecuaciones anteriores (4.35) y (4.33)

T
| M +

K |

= 2

(4.36)

i 1 2 i i i

T
M

+
T
K

= 2

(4.37)

1 i i 2

2
2 i
i i i i

= 2e
i

i

(4.38)

Recordar que si la normalizacin de los autovectores
i
fue hecha segn (4.13) y (4.14) el
producto

T
K =
2
(4.39)

Esto se debe a que la (4.8) puede expresarse para A =
i
y premultiplicarse por
i
T

T
|K-
2
M | = 0

(4.40)

i i i

T
K

2

T
M = 0

i i i
14
i
243
i

1
(4.41)

Determinacin prctica de modos y frecuencias

Este es, por s mismo, uno de los problemas ms complejos de la dinmica estructural.
Generalmente no es necesaria la resolucin de todas las frecuencias y sus modos, sino que
solo interesan las primeras q que representan las posibilidades ciertas de vibrar pues ne-
cesitan menos energa de excitacin.

Cuando el problema es de pocos GLD es posible resolver el determinante de (4.9) o (4.10)
obtenindose las n e
i
y luego los n
i
, donde n es el nmero de GLD. Para proble-
mas con muchos GLD se utilizan tcnicas numricas (aproximadas) para obtener los pri-
meros q pares de (e
i
;
i
). En la mayora de los casos prcticos de ingeniera q << n.

En las referencias bibliogrficas pueden encontrarse descripciones detalladas de varios
mtodos. Nosotros desarrollaremos solo uno que es relativamente simple y preciso.

Mtodo de Stodola-Vianello

Partiendo de la (4.8) y recordando que se cumple para cualquier mltiplo del autovector A
podemos escribir
4-10

T

A M A
K A e
2
M A = 0

K A = e
2
M A

(4.42)

(4.43)

Suponiendo un vector inicial A
0
conocido (cercano a A)

K A
1
= e
2
M A
0

A
1
~ K
-1
M A
0
(4.44)

(4.45)

que lleva a la frmula de recurrencia

A
i + 1
~ K
-1
M A
i
(4.46)

Nota: una descripcin detallada de este mtodo puede encontrarse en [4], parte 11.2
Vector Iteration Methods, as como prueba de la convergencia del mtodo.

Procedimiento de clculo

1- Se propone un vector A
0
inicial. Conviene que los valores correspondan a una defor-
mada suave, no obstante esto solo acelera la convergencia.
2- Se calcula

3- Se normaliza A
1

4- Se mejora la solucin

A
1
= K
-1
M A
0

A
1
= A
1
{A
1T
M A
1
}

1
2

i=1

A
i + 1
= K
-1
M A
i

(4.47)
(4.48)
(4.49)

5- Se normaliza la nueva solucin
A
i +1
= A
i +1
{A
i +1,T
M A
i +1
}

1
2

6- Se calcula la frecuencia correspondiente

(4.50)

1

A
i +1
K A
i +1

2

e
i +1
=

(4.51)

i +1
T

i +1

4-11

1 1
n
1
i i

1
0
1
7- Se controla la convergencia

e
i +1

e
i

s TOL
NO
se vuelve al paso 4

(4.52)

e
i +1

+SI

= A
i +1

; e
1

= e
i +1

(4.53)

Se obtuvo el primer autovalor y su correspondiente autovector.

Siempre: i ; A
i
; e
i
e

Obtencin de modos y frecuencias superiores
El procedimiento es el mismo que para el modo 1, solamente debemos garantizar que el
vector de arranque sea ortogonal a
1

Se supone que un vector inicial

n
A =

q
j

j

(4.54)

2
j =1

que para hacerlo ortogonal a
i
le restamos la componente q
1
1

A*
0
= A
0
q

(4.55)

2 2 1 1

Para determinar q
1
, premultiplicamos la (4.54) por [
1
T
M] y aplicando las propiedades de
ortogonalidad nos queda
|
T

M | A
0

= |
T

M |

q
j

j

(4.56)

1 2 1

T
M A
0

j =1
q =
1 2

(4.57)

y en general
1

T
M

T
M A
0

q =
i i +1

i

T
M

(4.58)
i i

y para los vectores de arranque (ortogonales a los previamente calculados)

A*
0
=
i 1
A
0

j =1

q
j

j

(4.59)
4-12

Procedimiento:
1. Se estima A
i
0k

2. Se calcula q
i
con (4.58)
3. Se calcula A
i
*
0k
con (4.59)
4. Se sigue con el mtodo de Stodola propiamente dicho

Resolucin de las ecuaciones de movimiento en estructuras con
mltiples GLD

Descomposicin y superposicin modal

Si bien desarrollaremos el mtodo para el caso ssmico, todo lo expuesto es valido para el
caso general.
La ecuacin que gobierna el comportamiento de una estructura de mltiples GLD es la

M D
&&
+ C D
&
+ K D = M {J a} (4.2)

Las vibraciones libres no amortiguadas se estudian mediante

M D
&&
+ K D = 0 (4.4)

al cual, segn lo expuesto en puntos anteriores, le corresponden n pares de frecuencias y
modos que son solucin del sistema de ecuaciones algebraicas

|K
2
M | = 0

(4.8)

Recordemos que los autovectores
i
son ortogonales respecto a las matrices de masa y ri-
gidez. Es por esto que pueden formar una base completa para el espacio de los desplaza-
mientos estructurales, es decir, es posible escribir

n

D =

i
y
i
(t )

i =1

(4.60)
4-13

i
j

i j

j

i

j
j

i j

j
j
j
j
j j

j j

j j
T
Donde y
i
(t) es un escalar funcin del tiempo a determinar, llamado respuesta generaliza-
da.
Sustituyendo (4.60) en (4.2) obtenemos

n
M

i
i =1
&y&
i
(t ) + C
n

i
i =1
y&
i
(t ) + K
n

i
i =1
y (t ) = M J a

(4.61)

Premultiplicando (4.61) por un
i
T
cualquiera, se cumple que:

n
M = M

= M
j
*

(4.62)

T T
j
i =1

n

T
K

=
T
K

= K
j
*

(4.63)
i =1

Y si, como es habitual, se est frente a matrices de amortiguamiento proporcionales y or-
togonales (ver punto 4.4)

n
C = C

= C
j
*

(4.64)

T T
j
i =1

Entonces la (4.61) premultiplicada por un
j
T
cualquiera queda como una ecuacin de 1
GLD

M
j
* &y&
j
(t ) + C
j
* y& (t ) + K * y
j
(t ) =
T
M J a(t )

(4.65)

Recordando que:

T
C = C
j
* = 2
j

j

(4.35)

T
K = K
j
* =
2

(4.39)

Y dividiendo ambos miembros por M
j
* =
j
M
4-14

y
j j T

j j

j

j

&y&
j

(t )

+ 2
j j

y&
j

(t ) +

2
(t )

T
M J

=
j

a(t )

(4.66)

M

Al valor

T
M J

Q
j
=

j

T
M
(4.67)

se lo denomina coeficiente de participacin modal.

Finalmente, el sistema (4.2), al ser proyectado segn la base formada por la matriz mo-
dal u = [
1
,
n
], queda en n ecuaciones diferenciales de 1 GLD, del tipo:

&y& (t ) + 2

y& (t ) +
2
y (t ) = Q

a(t )

(4.68)

i i i i i i i

Esta ecuacin puede resolverse con cualquiera de los mtodos vistos en la UT3. Una vez
obtenidos los n y
i
(t), puede obtenerse la solucin estructural mediante (4.60).

En general, los modos bajos son los que contienen menos energa elstica de deformacin
y por ende los que ms contribuyen a la respuesta estructural.
Usualmente

q
D(t ) =

i
y
i
(t )

i =1

; q < n
(4.69)

1- Se determinan los q primeros modos y frecuencias
2- Se resuelven las q ecuaciones (4.68) i=1,q
3- Se obtiene la historia en el tiempo segn (4.69) y sus derivadas
4- Se obtienen las fuerzas elsticas, de inercia y de amortiguamiento para cada instante

Integracin directa de las ecuaciones de movimiento

Como su nombre lo indica, este mtodo no requiere ninguna transformacin previa de las
ecuaciones de movimiento. Consiste, bsicamente en obtener la solucin en una cierta can-
tidad (discreta) de pasos de tiempo. Es por esto que tambin es llamada integracin paso
a paso.
4-15

0 0 0
1
r D
1
D
i
D = r
(D
i
i +1 i 6 i
0
0 i 2 i 3 i
i 4 i 5
K
0 2 i 3 i
Como ya adelantamos, existen numerosas variantes y algoritmos para la integracin num-
rica de las ecuaciones de movimiento, pero nuestros desarrollos estarn dados en el mto-
do de Newmark.
En esencia es igual al planteo dado en el Captulo 3, solo que en vez de tratarse de una
ecuacin de 1 GLD se trata de un sistema de n GLD.
A) Iniciales
1. Ensamblar las matrices M, C y K

2. Inicializar D , D
&
y D
&&

3. Seleccionar el paso de tiempo y los parmetros y |
At s T
mn
/10 ; =0,5 ; |=0,25

4. Calcular las constantes a
0
,a
7

5. Formar la matriz de rigidez efectiva:
K

= K + a

M + a C

(4.70)

B) Para cada paso de tiempo
1. Formar el trmino de carga efectiva en t+At (*1)

i +1
= F i +1 + M (a

+ a D
&
+ a D
&&
)

+ C (a

2. Resolver el desplazamiento en t+At

+ a D
&
+ a D
&&
)
(4.71)

i +1 i +1

(4.72)

*1: el trmino F i +1 = F i + { M J a(t )} se refiere a fuerzas generalizadas

3. Calcular aceleraciones y velocidades en t+At
D
&&

i +1

= a
i +1
D ) a D
&
a D
&&

(4.73)

D
&
= D
&
a D
&&

+ a
7
D
&&
i +1

(4.74)

# Notar que para cada paso de tiempo se debe resolver el sistema (4.72) lo que puede
hacer excesivo el costo de clculo y almacenamiento.

Respuesta mxima utilizando espectros de respuesta

Las ecuaciones desacopladas (4.68) pueden resolverse utilizando los espectros de respues-
ta: la mxima aceleracin ser
4-16

)
e
D

(S
e
=

D

(R
j
= )
a j
D
R
)
&y&
i
(t )

mx

= Q
i
(S
a
i

(4.75)

entonces, el mximo desplazamiento es

y (t )

=
Q
i
(S )

(4.76)

i
mx 2
a i
i

Podemos, entonces calcular los mximos desplazamientos (en todos los GLD) para el mo-
do j:

j
mx

=
j

y (t )

mx
Q
j

j 2 a j
j

(4.77)

D
j
= A

(S )

(4.78)

mx j 2
j

Donde A
j
es el vector de coeficientes de participacin modal correspondientes al modo de
vibracin j. Puesto que el mximo en cada grado de libertad no se produce en el mismo
instante, la respuesta total mxima no es la suma de los mximos de cada modo!!

mx
q
i
mx
i =1

(4.79)

Una forma muy utilizada (y precisa) de evaluar la respuesta mxima (desplazamientos,
velocidades, aceleraciones y esfuerzos) como combinacin de los mximos modales es:

denominada SRSS

mx
=

q
j
2

mx
j =1

(4.80)
Captulo 5 Anlisis de Construcciones con Efecto Dinmico de Viento
5-1

CAPTULO 5 - ANLISIS DE CONSTRUCCIONES CON
EFECTOS DINMICOS DE VIENTO

flujo

T(z)

L(z)
T(z): fuerza aerodinmica de
empuje, paralela a la direccin
del flujo
L(z): fuerza aerodinmica de
deriva, perpendicular a la direc-
cin del flujo

Generalmente son dominantes las fuerzas de empuje, pero para ciertas estructuras:

- bajo amortiguamiento
- acusada flexibilidad (poca rigidez)
- construcciones livianas (poca masa)
las acciones perpendiculares de deriva pueden tener efectos significativos que deben ser
contemplados en el anlisis.
CIRSOC 102-1 procedimientos simplificados

Validez:

a) construcciones prismticas o cuasiprismticas
b) primer modo dominante y de forma aproximadamente lineal
c) periodo fundamental T
1
> 1 seg
d) amortiguamiento < 0,01

Acciones paralelas a la direccin del viento

(procedimiento basado en el factor de rfaga)

x(z )
X max (z )

x
max
(z )
X max (z ) = x(z ) + x
max
(z )

x(z ): desplazamiento medio

x
max
(z ): desplazamiento fluctuante debido a la

turbulencia variable con el tiempo
5-2

z

Se define como factor de rfaga

entonces

G(z ) = 1 +

x
max
(z )

x(z )

X max (z ) = G(z ) x(z )

Hiptesis bsicas para poder reemplazar la accin dinmica del viento turbulento por un
procedimiento esttico equivalente.
a) Comportamiento elstico lineal de la estructura.
b) El modo fundamental de vibracin es una funcin lineal de la altura.
c) La contribucin de los modos superiores al primero en respuesta se considera
despreciable, por lo que G(z) = G = cte.
d) La velocidad media del viento es promediada sobre intervalos de una hora.
e) La variacin de la velocidad media del viento vara segn

V (z ) = V
0

ln

z

0,i

z
0,i

z

0,0706
ln

10

0,1

V
0
: velocidad bsica de diseo (m/seg)

z
0,i

z
0,i
: parmetro que depende de la rugosidad i
z
0,1
: z
0,i
para rugosidad I (referencia)

I) Presiones
La presin dinmica que incluye el efecto de la turbulencia del viento se determina
mediante

G: factor de rfaga
q'
z
= G c
z
c
2
q
0

c
z
: variacin por rugosidad y altura (art. 5.2.4.2, CIRSOC 102)

c
2
: factor por cambio de tiempo en velocidad media

(Tabla 3/pg. 20, CIRSOC 102-1)
q
0
: presin dinmica bsica (art. 5.2.3 , CIRSOC 102)
5-3

J

G = 1 + 1,234 K

B + r

se calcula mediante tablas y bacos que utilizan

una serie de parmetros auxiliares

En cada nivel se comparar:

q
z
> q
z

NO se adopta q
z
para los esfuerzos
SI se adopta q
z
para los esfuerzos
CIRSOC 102 q
z

= c
z
c
d

q
0

II) Aceleraciones

. .
X max (z ) = K
. .
(z )

X

K fig. 13, CIRSOC 102-1

(paralela al viento)

. .
(z ) valor medio cuadrtico de las aceleraciones (pg. 21)

X

III) Verificaciones (paralelas al viento)

. .

a- Dimensionado estructural

X max ; X
max

b- Verificacin de confort con grficos y
. .

q'
z

tablas en funcin de X
max
, X max y T
o
H
c- Verificacin de deformaciones admisibles

q
z
H
500

< X
max
<
H

350

Acciones perpendiculares a la direccin del viento

I) Resonancia

V

h
L(z)

torbellinos de Brnard - Karman

T(z)
5-4

z
La velocidad crtica del viento que produce el fenmeno de resonancia es

V
cr

=
d

S T
d: ancho de la superficie maestra (puede ser variable)
cilindros 0,20

S: n de Strouhal

T: perodo propio

prismas 0,25 a 0,30

V
cr
// T //

V
cr
T

Si V
cr
> 25 m/seg, entonces puede prescindirse del clculo de la resonancia.

De lo contrario:

Fuerzas de deriva:

(a la velocidad crtica)
: fraccin del amortiguamiento
L(z ) =
0,08
q

cr

z
d
h

q = 0,000613 V
2
kN

cr

V
cr

= V
cr

cr
m
2

Fuerzas de empuje: se admite distribucin uniforme
(a la velocidad crtica)

C
E
: coeficiente global de empuje CIRSOC 102

G: factor de rfaga correspondiente a V
cr

q
cr
con V
cr
//
T (z ) = T
z
= 0,8 C
E
G q
cr
d

Las fuerzas L(z) y T
z
obtenidas para la velocidad crtica (correspondientes al perodo per-
pendicular y paralelo respectivamente) se suman de la siguiente manera:
F (z ) = L(z )
2

+ T
2

y se debe comparar con las correspondientes obtenidas con q
z
o q
z
dadas para la velocidad
de diseo.
5-5

Criterios de confort en edificios que oscilan

El confort de los ocupantes de los edificios de gran altura que soportan la accin del viento
es un tema de primordial importancia en el diseo.
La respuesta humana al movimiento oscilatorio de las construcciones abarca una extensa
gama de reacciones, con efectos psicolgicos y fisiolgicos tales como mareos, ansiedad,
molestias visuales o temor, llegando hasta nusea aguda.
Por el contrario, otras persones sienten placer por el movimiento y la experiencia poco
usual de estar en una construccin que oscila. Adicionalmente, las oscilaciones excesivas
producen fisuracin de la tabiquera, rotura de vidrios de ventanas y cada de revestimien-
tos, lo cual incide negativamente en el valor de una propiedad y su rentabilidad.
Sin embargo, un edificio resultara demasiado costoso si se construye o equipa de modo
que pueda soportar sin movimientos perceptibles una tormenta con vientos huracanados o
un fuerte sismo.
En consecuencia, los movimientos son casi inevitables y el problema del diseo consiste
en mantenerlos dentro de los lmites aceptables para no perturbar el confort y el bienestar
de los usuarios. Por otra parte, el costo del edificio no debe superar los valores normales
de una aceptable economa.
Los factores que pueden producir vibraciones en un edificio son numerosos, tales como
maquinaria en funcionamiento defectuoso, paso de vehculos pesados por el lugar, impac-
tos en rampas, vientos fuertes, sismos, etc. los cuales pueden variar durante la vida de ser-
vicio de la estructura.

En general la aceleracin es el parmetro predominante para determinar aproximadamente
la naturaleza de la respuesta humana a las vibraciones. Las curvas de la fig. 5.3.1 grafican
los lmites de confort obtenidos del anlisis de un gran nmero de edificios altos, indican-
do las mximas aceleraciones aceptables para diferentes frecuencias, dependiendo del uso
o destino del edificio.
Los datos se obtuvieron para las aceleraciones pico de las ms fuertes tormentas ocurridas
durante un periodo de retorno de mas de 5 aos.
5-6

I
N
T
O
L
E
R
A
B
L
E

1
5
%

d
e

g

f (m)

0.40

0.30

0.20

0.10

NO SE
PERCIBE

T(seg)
1 2 3 4 5 6 7 8 9 10

GRAFICA DE CONFORT
fig. 5.3.1

AMPLITUD DE OSCILACION EN FUNCION DEL PERIODO
PARA DISTINTOS VALORES DE LA ACELERACION
EN PORCENTAJE DE g
5-7

La tabla 5.3.1 a continuacin da idea de la magnitud de las percepciones que se obtuvieron
en experimentos realizados para diferentes niveles de aceleraciones (segn Khan Parmelec
y segn Chang).

TABLA 5.3.1
PERCEPCION HUMANA ACELERACION (m/seg
2
)
(segn Khan y Parmelec)
No perceptible 0,004 g

Levemente perceptible 0,004 g < 0,0075 g
Perceptible
0,0075 g < 0,02 g
Molesta 0,02 g <
(segn Chang)
No perceptible 0,005 g
Perceptible 0,005 g < 0,015 g
Desagradable 0,015 g < 0,05 g
Muy desagradable 0,05 g < 0,15 g
Intolerable 0,15 g <
5-8

La sensacin subjetiva y el comportamiento humano afectado por las diferentes acelera-
ciones se indica en la tabla 5.3.2 (segn Yamada y Goto).

TABLA 5.3.2
NIVEL DE PERCEPCION HUMANA
RANGO ACELERACION
(m/seg
2
)

EFECTO

1 0,05 La gente no percibe el movimiento.

2 0,05 0,10 Las personas sensibles perciben el movimiento.
Los objetos colgantes pueden moverse algo.
3 0,10 0,25 La mayora de las personas perciben el movimien-
to. La oscilacin puede afectar el trabajo de ofici-
na. La exposicin de larga duracin puede producir
malestar.

4 0,25 0,40 El trabajo de oficina se vuelve difcil o casi impo-
sible. Aun se puede caminar.
5 0,40 0,50 Se percibe fuertemente el movimiento. Hay difi-
cultad para caminar normalmente. Las personas de
pie pueden perder el equilibrio.
6 0,50 0,70 No se tolera el movimiento y no se puede caminar.
7 - 8 > 0,85 Los objetos caen y pueden lastimar a las personas.
5-9

A
C
E
L
E
R
A
C
I
O
N

m
/
s
e
g

10
100 %g

4.4
44 %g

1.5
1

0.15

0.11
15 %g

1.5 %g

1.1 %g

0.1 %g

0.01

10 5 2.5 1 PERIODO (seg)

fig. 5.3.2
GRAFICA DE CONFORT

INTERACCION DE LAS VARIABLES:
PERIODO (T), AMPLITUD (f) Y ACELERACION (x)
5-10

GRADUACIN
DEL EFECTO
SOBRE LAS
PERSONAS

ACEPTABILIDAD
EN
EDIFICIOS
ESTADO
PREVISIBLE
DE LA
ESTRUCTURA

A
INTOLERABLE
B
MUY POCO
TOLERABLE

NO COLAPSO

NO DAOS
LOCALES

C
DEMASIADO
PERCEPTIBLE
D
MUY
PERCEPTIBLE
E
PERCEPTIBLE

F
ESCASAMENTE
PERCEPTIBLE
G
NO
PERCEPTIBLE

SITUACION
LIMITE

EN TAREAS
INDUSTRIALES
PESADAS
EN PERIODOS
BREVES EN
VIVIENDAS
PERIODOS
LARGOS
EN VIVIENDAS
PERIODOS
LARGOS
EN VIVIENDAS

FORMACION
DE GRIETAS
FORMACION
DE FISURAS
SIN INFLUENCIA
EN ESTRUCTURAS
CORRIENTES
SIN
INFLUENCIAS

SIN
INFLUENCIAS

CUADRO ADJUNTO A GRAFICA DE LA fig. 5.3.2
Captulo 6 Anlisis de Construcciones con Efectos Ssmicos
6-1

CAPTULO 6 - ANLISIS DE CONSTRUCCIONES CON
EFECTOS SISMICOS

Definicin numrica de la accin ssmica

La definicin correcta de la accin ssmica es un problema al que se le debe dar la mayor
importancia en un anlisis estructural ssmico.
La solucin a este problema parte, generalmente, de los datos experimentales proporcio-
nados por la sismologa.

Parmetros de
terremotos pasados
en una regin
Prediccin de las
caractersticas ssmicas
de los terremotos que
afectarn dicha regin

Tradicionalmente, la fuerza destructiva de un terremoto ha sido expresada en funcin de la
aceleracin mxima del terreno, pero se han observado daos moderados a aceleraciones
muy altas.
Por lo tanto, se deben tener en cuenta otras caractersticas:
- intensidad, duracin
- contenidos de frecuencias
- secuencia de choques, etc.

Por ejemplo, se han utilizado, adems de los valores mximos de la aceleracin, velocidad
y desplazamiento del terreno, el espectro de amplitudes de Fourier, el espectro de seudove-
locidades, intensidad espectral (Housner), el valor medio cuadrtico de las aceleraciones
en la fase fuerte de un acelerograma, etc.
Adems, la accin definida estar totalmente relacionada al tipo de anlisis estructural que
se va a realizar:

Anlisis lineal
desacoplamiento

modal
espectros de

respuesta

Anlisis no lineal
integracin
directa

acelerogramas
6-2

a d
# Tambin existen definiciones de acciones ssmicas a partir de la teora de procesos esto-
csticos, en donde la respuesta estructural en s se obtiene probabilsticamente.

No profundizaremos sobre estos temas, pues escapan a nuestro curso. (ver ref. [1 y 2])

Definicin mediante espectros de respuesta

Es la forma ms usual de definir una accin ssmica, dado que se obtienen descripciones
de las caractersticas ms importantes de la respuesta estructural sin la necesidad de dispo-
ner de una historia en el tiempo de la excitacin y la respuesta.
Otra ventaja es que un espectro puede modificarse en base a las caractersticas del lugar
de emplazamiento de la estructura, sin necesidad de conocer los detalles de la excitacin.

Espectro
de respuesta
representaciones grficas de valores
aproximados de la respuesta mxima
de un sistema de 1 GLD lineal y elstico:

x+ 2e x+ e
2
x = a(t )
(2.11)

S
d
(a, e , ) =
. .

1

e
.

t

a(t ) e
e (t t )
sen e (t t ) dt

0

max

(3.49)
S
v
(a, e , ) = e S
d
(a, e , )

(3.52)
S (a, e , ) = e
2
S (a, e , )

(3.53)

Puede realizarse un anlisis de las variaciones de los valores de las curvas S
d
, S
v
, S
a
, en
funcin de las caractersticas de la estructura.

Recordar que e =
k
m

Para frecuencias propias altas (en comparacin con la del movimiento del terreno), la
(2.11) quedara (despreciando los dos primeros trminos):
e
2
x ~ a(t )

estructuras
muy rgidas

si e ;
S
a

a
max
(copia la aceleracin mxima del terreno)

6-3

Para frecuencias bajas, por el contrario, la (2.11) queda:

. .

x ~ a(t )

estructuras
si e 0;
S
d

d
max

muy flexibles
(copia los desplazamientos mximos del terreno)

Para frecuencias intermedias, se produce una amplificacin del movimiento del terreno en
su paso a travs del filtro estructural.

|
a
( , T ) =

S
a
;

a
max

|
v
( , T ) =

S
v
;

v
max

|
d
( , T ) =

S
d
;

d
max

(6.1)

Las condiciones locales del terreno, tales como grosor y propiedades de los estratos que se
encuentran entre la roca firme y la estructura modifican los espectros de respuesta.

Conjunto de
espectros de
respuesta
registrados para
una regin

suavizado y promediado
de valores

Espectros
de Diseo

+
INPRES CIRSOC 103
Parte 1 / Cap. 7

S
a

g

b

a
s

T
1
T
2

Para = 0,05 (lnea llena)

S
a
= a
s

+ (b a
s

)
T
T
1

T s T
1

6-4

2
2

3
S
a
= b T
1
< T s T
2

2
S = b
T
2

T
2
< T

a

T

zona ssmica

a
s
, b, T
1

y T
2
tabla 4

f

suelo

Para casos con 0,005 s < 0,05 (lnea a trazos)

S
a
= a
s

+ ( f
A

b a
s

)
T
T
1

T s T
1

S
a
= f
A
b

T

T
3

T
1
< T s T
2

S
a
= 1 + ( f
A

1)

b
2

T
2
< T

T

T

f
A
: factor de amplificacin por amortiguamiento

f
A
=

5
(expresado en %)

Definicin mediante acelerogramas

a) Acelerogramas reales
Basar un clculo ssmico en uno o varios registros disponibles en una zona implica
un alto riesgo de definicin incorrecta de la accin.
b) Acelerogramas sintticos
Tiene grandes ventajas respecto al anterior:
- posibilidad de generar seales de corta duracin que tengan el mismo efecto en
las estructuras que el del terremoto real que se quiere simular economa
computacional
- al estar generados en funcin de un espectro de diseo, se tiene en cuenta
(aprox.) las condiciones locales del suelo
6-5

La generacin de acelerogramas sintticos requiere procedimientos matemticos basados
en procesos estocsticos y expansiones en series sinusoidales.
6-6

S S
S
a a
Existen muchos procedimientos y algoritmos que permiten generar acelerogramas sintti-
cos a partir de un espectro dado.

INPRES CIRSOC 103 (Parte 1 Cap. 14.3.1)
Los acelerogramas deben satisfacer:
a)

a(t )
max

a(t )
max
> a
S

d

a
S
: ordenada al origen del espectro correspon-
diente
d
: factor de riesgo, art. 5.2

b) El espectro elstico (para = 0,05) deber cumplir:

Area
C

Area
D

CALCULADO DISEO

0,05
C T
0

a

Area
C

>
d

0,05

Area
D

D T
0

a

C
> 0,7 S
D

d

para todos los puntos

Cantidad de acelerogramas a emplear:
Grupo A y B n > 3 acelerogramas
Grupo A
0
n > 4 acelerogramas

Para diseo y verificaciones se promediarn las envolventes de solicitaciones y deforma-
ciones obtenidas para cada acelerograma, pero en dicho promedio no se incluirn valores
inferiores al 85% del mximo encontrado.
6-7

)

Mtodos de anlisis segn INPRES CIRSOC 103

Procedimientos con fuerzas estticas equivalentes

a) Mtodo esttico (art. 14.1)
b) Procedimiento aproximado (cap. 16)

Mtodos dinmicos

a) Anlisis modal espectral
b) Superposicin modal paso a paso
c) Integracin directa paso a paso

6... Mtodos dinmicos INPRES CIRSOC 103

a) Anlisis modal espectral
- excitacin ssmica traslacional
- materiales lineales y elsticos (segn art. 12.1 / pg. 47)
- (S
a
)
i

= S
a

d
g

S
a
: del espectro correspondiente

d
: factor de riesgo

g: aceleracin de la gravedad
- se podr considerar capacidad de disipacin de energa por deformaciones
anelsticas de la estructura
(S
a i

=
S
a

d

g R
- las deformaciones calculadas con este criterio debern ser amplificadas
multiplicndolas por la ductilidad global
- para el caso anterior, es imprescindible que la estructura se comporte en
forma uniforme, de manera de garantizar la ausencia de concentracin de
rtulas plsticas

Se debe verificar que el corte (con anlisis modal espectral) en cada direccin no resulte
inferior al 75% del corte obtenido con el mtodo esttico; de no verificarse se incrementa-
ran los efectos del anlisis modal de la siguiente manera:
6-8

=

=
M {1}

T

=
2
e
2
O
M
W

R *
modal

=
0,75 V
0 est

R
modal

espectral V
modal

0
espectral
espectral

Modos a considerar: q > 3, pero adems se debern incluir todos los modos cuya contribu-
cin al total sea mayor que 5% de la contribucin del modo fundamental.

Modelo de edificio de cortante

x
n

k
n
Q
i

n
T
m
r

ri
i r =1
T n

(6.1)

i
M
i

m
r
K

m
r ri

i =1

0

1i

x
r

w
r
k
r

M

x
r -1
con
M = m
r

;
i

=
ri

h
r
m
2

O

m
1
x
1

k
1

a (t )

0

m
r

ni

# observar que si en (6.1) se susti-
tuye la masa m
r
por el peso w
r
el
si
i
M
i
= 1

n n

valor de Q
i
no cambia Q
i
=

m
r

ri
=

w
r

ri

i =1 i =1

Peso efectivo modal

n

w
r
r =1
i
n

2

ri

si

q

W
i

> 0,9

n

w
r

(6.4)

w
r

ri
r =1

q apto

Desplazamientos modales mximos
i =1
1
i =
2
1
3
peso total de la
estructura

(x )

= y (t )

= Q

(S )

a i

(6.5)
i max i i max i i
2
i
6-9

) )
)
)
S
1i
)
..
a i
M
Desplazamientos modales relativos entre pisos

A
ri
= (x
r i
(x
r 1 i
(6.6)

Fuerzas ssmicas equivalentes (modales)

S
..

(6.7)

F
ri

=

x ri

= m
r

max

pero
r: nivel
i: modo

..

..

(6.8)

x ri

=
ri
x
i
(t )
max

ma
x

Como

xi (t )
max

= Q
i

(S
a i

..

(6.9)

x ri

=
ri

Q
i

(S
a i

max

(S )

(6.10)

F
ri
=
ri
Q
i
(S
a
)
i
m
r

=
i
Q
i

w
r

g
fuerza ssmica equivalente
en el nivel r para oscilaciones
en el modo i
en funcin del peso w
r

y la fraccin de g

Vectorialmente:

F
S

F
S
=

M

S
= M

T
M J

i

( )

i
F
ri i

T
M

S
a i
(6.11)

i i

F
S
ni

Q
i

6-
10

F
i i

S
= M Q
i
(S
a i

(6.11)
6-
11

)
F F
a i
= =
Cortante modal

En nivel r:
n
=

S

V
ri
F
ji

j =r
(6.12)

En la base:
n
=

S

(6.13a)

V
i
F
ri

r =1
o tambin: (S )

V
i
= W
i

g
(6.13b)

Comparacin con INPRES CIRSOC 103

d

V
m
=

S
am
W
m

R
(Frmula reglamento)
d
: factor de riesgo
R: coeficiente de reduccin por disipacin
i = m

V = f

(S
a

i
W

(Frmula 6.13b, multiplicada por un factor f
1
)

i 1
g
i

Distribucin en altura de V
m
segn INPRES CIRSOC 103:
W
k
|
km

m: modo
k: nivel
F
km
=

n
W
i
|
im
i =1
V
m

F
km
: fuerza ssmica equivalente para el modo m en el nivel k

Segn la frmula (6.10) y (6.13a):

n
V
i
=

ri
r =1

Q
i

w
r

(S )

a i

g
Q (S )

=
i a i

g

n

ri
r =1

w
r

si queremos emplear la notacin

S
ri ri

V
i

S
ri ri
Q (S )

i a i

g

n

ri
r =1

w
r

comparndola con (6.10) se deduce que
6-10

v

F
i
(h
j

F
q
)
=
S
o

ri
=

ri
w
r
n

W
k
|
km

n

notacin INPRES - CIRSOC 103

ri
w
r

W
i
|
im

r =1 i =1

Momento de vuelco modal (para el nivel r)

n

M
ri
=
j
h
r
)
j =r +1

si r = 0 M
0i
v
= momento de vuelco en la fundacin

Momento de torsin modal

e
r

C.R. C.M.

M
ri
= e
r
V
ri

corresp. al piso r

Esfuerzos seccionales modales

Para cada subestructura vertical con la parte proporcional a su rigidez que le corresponde
de la fuerza ssmica equivalente, se calculan los esfuerzos seccionales para cada modo en
estudio.

S
p ri
subestructura "P"
nivel "r"

o
1
+ o
2

+ Ko
P

+ Ko
m
= 1

m: n de subestructuras verticales

(se deben adicionar los efectos torsionales si los hubiese)
CALCULO DE LA RESPUESTA TOTAL

R
max
=

(R
j

j =i

2
max

Estruc Dinam

Cargado por

Copyright:

Formatos disponibles

Estruc Dinam

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Estruc Dinam

Cargado por

Copyright:

Formatos disponibles

Departamento de Ingeniera Civil

Facultad Regional Paran Universidad Tecnolgica Nacional

También podría gustarte

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.