CURSO DE INGENIERA SSMICA

CONTENIDO
Captulo 1
Conceptos de la Sismologa y la Sismicidad
1.1.- Sismologa
1.2.- Origen y propagacin de los movimientos ssmicos tectnicos
1.3.- Medicin de los movimientos ssmicos
1.4.- Espectros de respuesta dinmica y espectros de diseo ssmico
1.5.- Modelos matemticos
1.6.- Daos tpicos originados por los sismos
Captulo 2
Introduccin a la Dinmica Estructural
2.1.- Vibraciones mecnicas de sistemas con un grado de libertad
2.1.1.Vibracin libre sin amortiguamiento
2.1.2.Vibracin libre con amortiguamiento
2.1.3.Vibracin forzada sin amortiguamiento
2.1.4.Vibracin forzada con amortiguamiento
2.2.- Sistemas con varios grados de libertad
2.3.- Respuesta Ssmica no lineal
2.4.- Anlisis Modal
2.5.- Aplicaciones
Captulo 3
Mtodos Numricos para la Determinacin de los Modos de Vibracin
3.1.- Mtodo de Stodola y Vianello
3.2.- Mtodo de Newmark
3.3.- Mtodo de Holzer
3.4.- Mtodo de las Matrices de Transicin
3.5.- Mtodo de Jacobi
3.6.- Aplicaciones
Captulo 4
Diseo Ssmico
4.1.- Criterios de Diseo Ssmico
4.2.- Estructuracin y Condiciones de Regularidad
4.3.- Respuesta Ssmica de los materiales estructurales ms utilizados
4.4.- Control del modo de fallas (Pushover)
Captulo 5
Aisladores y disipadores ssmicos
5.1.- Estudio del Amortiguamiento
5.2.- Amortiguamiento de los edificios
5.3.- Soluciones numricas para edificios con un grado de libertad
5.4.- Soluciones numricas con varios grados de libertad
Captulo 6
Anlisis Ssmico de Estructuras usando programas de Computacin
6.1.- Anlisis Ssmico de un edificio de varios niveles

6.2.- Anlisis Dinmico de una Cimentacin para maquinaria vibratoria

6.3.- Anlisis ssmico dinmico de un tanque elevado

Sismologa y sismicidad.
SISMOLOGA

Los temblores son vibraciones naturales que se producen en o debajo de

la corteza terrestre. Ellos son provocados por actividad volcnica, colapso
de cuevas, actividad tectnica, as como una gran cantidad de pequeos
temblores que se deben a vibraciones inducidas artificialmente por un
trfico intenso en las calles, por trenes y tranvas o trastornos anlogos.
Las explosiones tambin determinan terremotos artificiales. En general los
temblores artificiales se dejan sentir en superficies limitadas, en cambio
los temblores naturales se hacen sentir en extensas reas. Entonces, de
todos estos tipos de temblores, los ms importantes en ingeniera ssmica
son los de origen tectnico, puesto que la energa que se libera y el rea
afectada, como se apunt, suelen ser mucho mayores.
El origen de los temblores tectnicos es materia de amplia discusin. Se
tiene actualmente dos teoras contradictorias: La primera atribuye el
origen de estos movimientos a deslizamientos de fallas geolgicas o sea
que grandes masas de rocas sufren un sbito desplazamiento por falla o
fracturamiento, debido a esfuerzos internos y continuados que actan por
largo tiempo hasta que van ms all del lmite de deformacin elstica de
las rocas. Y la segunda lo atribuye a cambios de fase en las rocas. La
primera teora se haba aceptado ampliamente, pero en 1963 Evisn
revivi la teora propuesta por Bridgman en 1945. Esta mantena que las
transiciones polimrficas locales repentinas podran proveer un
mecanismo para temblores. Evisn afirm que los deslizamientos en fallas
geolgicas son efectos y no causa de los temblores. Sin embargo, la teora
del deslizamiento de fallas sigue siendo la de mayor aceptacin por parte
de sismlogos e ingenieros; esta preferencia se basa en que la mayora de
las evidencias encontradas hasta ahora la favorecen ms que a la teora
de cambios de fase. Es necesario establecer que no hay suficientes datos
para excluir ninguna de las dos teoras.
El fenmeno da origen a ondas ssmicas que se propagan a travs de las
formaciones geolgicas que constituyen el terreno, sufriendo reflexiones y
desviaciones cada vez que encuentran medios cuya formacin geolgica
es diferente. Entonces las vibraciones se inician en un rea limitada y se
propagan en todas direcciones, esta rea central de iniciacin, bajo la
superficie terrestre situada verticalmente encima de ella, donde la
sacudida es intensa, se llama epicentro o rea epicentral. Los focos de la
mayor parte de los terremotos se encuentran a profundidades menores de
16 kms. Sin embargo, muchos terremotos se han originado a
profundidades comprendidas entre 16 y 48 kms. Y recientes estudios han
dado a conocer que algunos se originan a profundidades mayores, por
ejemplo ms de 640 kms.
La energa liberada durante el fallado produce varias clases de ondas, una
de ellas es la onda primaria (P), es una onda de tipo longitudinal o de
compresin; las partculas de las rocas vibran hacia atrs y hacia delante,
en la direccin de la propagacin de la onda. Una segunda clase de onda,
es la transversal, las partculas de las rocas vibran en ngulo recto con la
direccin de propagacin, esta se denomina onda (S), secundaria o
transversal. Las ondas P se propagan ms de prisa que las ondas S, y por
lo tanto llegan a la estacin sismo grfica antes que stas. Un tercer tipo

de onda, es la onda inducida que se desplaza a lo largo de la superficie

superior de la roca afectada, esta onda es ms lenta y mas larga que
cualquiera de las ondas P o S, y de mayor amplitud, se llama onda de
Rayleigh o superficial.
Las vibraciones de los temblores se registran por medio de aparatos
llamados sismgrafos de los que estn en uso modelos muy numerosos.
Las vibraciones recogidas por los sismgrafos suelen registrarse en papel
fotogrfico, como una serie de lneas en zig-zag. Estos registros muestran
el impulso vibratorio y el tiempo de iniciacin y duracin; tambin indican
la llegada de los diferentes tipos de ondas. Los movimientos de tierra
suelen ser de corta duracin. Para el estudio de los temblores es
absolutamente necesario una red de estaciones sismo grficas.
La magnitud como la intensidad ssmica miden el poder destructor de los
movimientos telricos. La magnitud mide la energa total desencadenada
por un temblor; la intensidad es una medida local de la capacidad
destructiva del movimiento y para un mismo temblor, varia de un lugar a
otro. Es sabido que la magnitud en la escala de Richter se calcula a partir
de medidas instrumentales, no constante el resultado que se obtiene es
funcin de la naturaleza del terreno en donde se apoya el sismgrafo en
cuestin y otras variables.
La intensidad de la escala de Mercalli Modificada, es en buena parte
apreciativa y depende de factores tales como intensidad de poblacin,
estado de constructores y otras variables.
Este tipo de escala es como la del observatorio Nacional de Tacubaya. En
ella, a partir del registro de un sismgrafo se calcula la aceleracin
mxima que sufri el terreno en el lugar en donde esta instalado el
instrumento.
Se compara con una escala de acelarciones y segn el rango en que se
halle la aceleracin calculada, as habra sido la intensidad del sismo
(Tabla 1). La principal objecin que se hace este tipo de escala radica en
que el poder destructor de un sismo no depende nicamente de la
magnitud de aceleracin del terreno, si no que varia directamente con la
duracin del movimiento, dependiendo tambin del contenido de
frecuencias. Tal objecin se salva en el tipo de escala propuesta por
Benioff y refinada por G.W.Housner.
Tabla 1.
INTENSIDADES
III
IV
V
VI
VII
VIII

ACELERACION cm / seg
0.5
1.0
2.5
5.0
10.0
25.0

A
A
A
A
A
A

1.0
2.5
5.0
10.0
25.0
50.0

Por lo que respecta a la localizacin geogrfica de los sismos, los mapas

que registran la localizacin de los epicentros de los terremotos, indican
dos fajas principales, donde se han originado la mayor parte de los
terremotos recientes. La primera parte de los terremotos, bordea la
cuneca del Pacifico y coincide aproximadamente con la distribucin de los
volcanes activos o de actividad reciente, y con la faja de montaas
jvenes y en formacin. La segunda zona principal de presentacin de
terremotos se extiende desde el sur de Espaa, a travs del rea del
mediterrneo y continua hacia el oeste a lo largo de las montaas del
Himalaya, hacia Asia Oriental, donde se ramifica; la rama principal pasa
hacia el sur, a travs de la regin Mayala, Hasta las Indias Orientales
Holandesas, unindose con la faja que bordea el Pacifico. Ver mapa 1.
Las manifestaciones ssmicas de Mxico.-Entre las zonas ssmicas se
distinguen las islas Maras en el Noroeste y el Istmo de Tehuantepec en el
Sureste. Del meridiano 96 hacia Chiapas los fenomenos ocurren de
manera singular. Los fracturamientos que atraviesan la porcin continental
en la regin que abarca los paralelos 19 y 20 constituyen ramificaciones
de las grandes fallas existentes en la parte que Mxico ocupa dentro del
cinturn Circunpcifico.
Actualmente las costas del Pacifico en contaste con las del Golfo de
Mxico, acusan una gran movilidad tectonica, que unido a lo dems
procesos geolgicos determinan el elevado ndice ssmico de nuestro pas.
Sin embargo existen muchos criterios sobre la generacin del sismo y
entre ellos contradictorios, como se indico. Ver mapa 2.

Grfica
Mapas
Grficas

Con la somera exposicin anterior, se puede observar que los fenmenos

ssmicos son complicados e irregulares, no pudiendo predecir exactamente
el tiempo, la intensidad o el lugar de temblores futuros todo esto hace
complejo el problema ssmico en estructuras. De aqu que normalmente de
proponer ciertos conceptos y expresiones para valuar la afectabilidad
estructural.
Sin embargo, las investigaciones no cesan, as por ejemplo por lo que
respecto a la idealizacin de temblores, tenemos que el avance de la
sismologa instrumental en estos ltimos aos ha proporcionado
instrumentos de presin para registrar el movimiento del suelo, siendo
tiles para los investigadores ya que de esta manera les ha proporcionado
la forma de obtener la variacin con el tiempo de las perturbaciones que
afecten a la estructura. La gran irregularidad del movimiento del terreno,

ya mencionada, durante los temblores causada por el mecanismo de la

generacin y por las reflexiones de las ondas ssmicas, as como del
desconocimiento previo de la vibracin de la excitacin con el tiempo, han
sugerido la idealizacin de temblores como procesos estocsticos. Para la
simulacin de estos registros de temblores, los investigadores han
utilizado diversas tcnicas acudiendo a dispositivos mecnicos, o
simulando registros en computadoras analgicas o digitales, otros han
propuesto modelos estocsticos analticos de la excitacin para estudiar
estadsticamente algunos parmetros de los registros y/o de la respuesta
estructural. Cualquiera que sea la tcnica de simulacin o el modelo
matemtico que se use, los registros simulados y/o la respuesta
estructural debern compararse con la veracidad que se obtenga de los
registros reales para aceptar los simulados. Para una informacin mayor
ver el modelo Modelo estocstico para registros de temblores en terreno
duro ponencia presentada por Octavio A. Rascon y C.A. Cornell en el 2
Congreso Nacional de Ingeniera Ssmica (Veracruz, Ver. 1968)
Actualmente se conviene en que se puede idealizar un sismo como un
movimiento horizontal del terreno en una sola direccin, cuyo movimiento
esta descrito por el acelerograma registrado para la direccin en estudio,
si se conocen las caractersticas de una estructura. Pero al disear una
estructura lo que interesa conocer es la forma como se comporta ante
temblores futuros cuyos acelerogramas defieren de los registrados. Por
esta razn se hace uso de las respuestas mximas de estructuras de un
grado de libertad y diferentes amortiguamientos, de varios acelerogramas,
introducindose el concepto espectro.
ESPECTROS DE RESPUESTA DINAMICA.
En 19933 M.A. Biot propuso que el uso de espectros para evaluar los
efectos de los movimientos ssmicos del suelo sobre estructuras
simplificadas. El concepto de espectro fue modificado, refinado y
expandido grandemente por G. Housher.
Para ilustrar lo que es un espectro, vamos a suponer una base movible la
cual se encuentra fijos u a serie de pndulos con distintos periodos cada
uno. La longitud y los periodos de los pndulos son incrementados de
izquierda a derecha.
Ahora si la base es movida continuamente con un movimiento anlogo al
que le ocurre al suelo para un sismo dado, la respuesta mxima obtenida
para cada pndulo ser registrada y suceder en cualquier tiempo; lo que
llamamos respuesta mxima puede ser: un desplazamiento angular o
lineal, cortante, aceleracin absoluta etc. O cualquier concepto que se
relacione con estos. Si la respuesta mxima es trazada, en el eje de las
ordenadas y el periodo en las abscisas se tendr un espectro como los
ilustrados en las hojas siguientes. Es sorprendente la similitud de las
curvas de respuesta para un gran numero de sismos, la mayor variacin
ocurre en las ordenadas verticales y ellas dependern del sismo y el del
lugar del aparato registrados.

Entonces la grafica que tiene como abscisas los periodos naturales (o las
frecuencias) y como ordenadas las respuestas mximas de los valores
absolutos del desplazamiento relativo, de la velocidad relativa, de la
aceleracin absoluta etc., se llamara espectro de respuesta.
Quiz el espectro de aceleracin facilitara ms la visualizacin al
problema. Es importante hacer notar que los espectros de respuesta, para
construirse, deben verse los trabajos del Dr. Luis Esteva.
Cuando se tiene suficientes acelerogramas durante un periodo grande de
aos es posible elaborar un espectro ssmico, pero como no siempre se
dispone de acelerogramas suficientes, ya sea por la falta de instrumentos
adecuados, en el camino que se sigue es considerar una serie de impulsos
distribuidos al azar en el tiempo, esta distribucin se corrige de la
observacin de un acelerograma.
Estos espectros sirven para el diseo de estructuras de cualquier nmero
de grados de libertad para una combinacin de las respuestas mximas de
sus modos teniendo cierto margen de seguridad para sismos futuros.
Ahora vemos la construccin de espectros; supngase una estructura
elstica de un grado de libertad con amortiguamiento lineal, cuya base
esta sujeta a una aceleracin a () descrita por un acelerograma de un
temblor para un componente dada (Grafica Aceleracin Tiempo, del
terreno). La ecuacin de movimiento que se obtendr es como la (19).
x ' '2x ' w 2 x A( )

La solucin de esta ecuacin est dada por la siguiente expresin.

x

1
w
2

A( )

(T ) sen w 2 2 (T )d

El desplazamiento mximo es el siguiente que llamaremos SD

SD

1
w2 2

A( )
0

(T ) senw(T )d
MAX

T es el tiempo para el cual la integral es mxima.

En la misma forma se obtiene la velocidad relativa.
SV

s D
T

MAX

Y la aceleracin absoluta.
SA

2 sD
T

A( )

2
MAX

De esta manera los diferentes valores que tengan de SD, SV, SA para
diferentes periodos nos darn los espectros de desplazamiento, de
velocidad y aceleracin respectivamente para un valor constante del
amortiguamiento (vase espectros C, B, A).
Advirtase que S es funcin de las caractersticas del temblor, de la
frecuencia natural y amortiguamiento de la estructura.
Haciendo algunas simplificaciones se demuestra que las expresiones para
los espectros de velocidad y aceleracin en funcin del espectro de
desplazamiento relativo son:
SV wS D

S A w2 S D

Debido a las simplificaciones tambin se les llama espectros de seudovelocidad y seudo-aceleracin respectivamente, o seudo-espectros.
Influencia del Amortiguamiento en los Espectros.- El amortiguamiento
reduce las ordenadas de los espectros, como se observa en las figuras, y
dicha reduccin est en funcin del grado de amortiguamiento. Para
valores de = 0.2 que es la relacin entre el amortiguamiento real y el
crtico, se reduce fuertemente la respuesta estructural y la sensibilidad del
espectro, a espectros elsticos de aceleracin de cada una de las
componentes de los 3 registros, para diferentes amortiguamientos.
En el caso de las componentes horizontales se presentan tambin con
fines comparativos el espectro del reglamento (1966) correspondiente a la
aceleracin mxima de diseo (estructura del grupo A tipo 3).
Grfica
Sismo del 9 de diciembre de 1965. Edificio Atizapn. Tlaltelolco. D.F. (CIM)

Nota.En el captulo VII, se dibujan los espectros de diseo segn el

reglamento de construcciones para el D.F. de 1976.
2 Grficas

CAMBIOS DE PERIODOS PEQUEOS.

De acuerdo con mediciones experimentales se ha determinado los

siguientes valores de amortiguamiento, expresados en relacin al
amortiguamiento crtico,
0.08 (8%) en edificios de concreto.
0.03 (3%) en edificios de acero con juntas soldadas.
(sin recubrir de concreto).
Para juntas remachadas el amortiguamiento es mayor.
Los espectros obtenidos de acelerogramas en California (U.S.A.) presentan
las siguientes caractersticas:
A).B).C).-

Para amortiguamiento nulo se tienen numerosas oscilaciones con

mximas irregulares agudos.
Las oscilaciones disminuyen generalmente a medida que aumenta el
amortiguamiento.
Para amortiguamientos comprendidos entre 0.05 y 0.1 de , los
mximos desplazamientos, velocidades y aceleraciones espectrales
son del orden de 1, 1.5, y 2 veces las correspondientes del terreno.

ESPECTRO DE DISEO.Generalmente
los
espectros
obtenidos
presentan variaciones bruscas como puede observarse en las figuras
mencionadas, por lo cual no resultan prcticas desde el punto de vista del
diseo, adems debe considerarse un cierto margen de seguridad para
temblores futuros. De tal manera que puede adoptarse un espectro
obtenido como la curva media o envolvente terica.

En los espectros del reglamento de las construcciones del D.F., se

consider una envolvente y se multiplic por un factor considerando las
intensidades posibles del sismo; estos espectros se ilustran al principio de
la aplicacin del anlisis modal del edificio.
Generalmente la respuesta elstica de estructuras utilizando espectros de
temblores intensos difiere de la especificada en los reglamentos, lo que se
considera ocasionada por comportamiento inelstico, para tomar en
cuenta este comportamiento en los espectros de diseo, las ordenadas se
reducen en funcin de un factor de ductilidad. Los factores de ductilidad
se dan en el reglamento de construcciones (1976), en el artculo 235.

INTERACCIN DINMICA SUELO-ESTRUCTURA.

Respecto a la interaccin que hay entre el suelo y una estructura sujeta a

una perturbacin proveniente del suelo, se dan algunos conceptos para
tener una idea del problema y un cierto criterio, as como diferentes
modelos matemticos representativos, propuestos por varios
investigadores.
En una regin de la superficie terrestre al recibir un tren de ondas ssmicas
sufre desplazamientos y deformaciones, las direcciones, amplitudes y
frecuencias de los desplazamientos dinmicos producidos son
enteramente casuales, pero dependen de las caractersticas del tren de
ondas y de la naturaleza del terreno. Sin embargo si en esta regin
existiera una estructura, los desplazamientos dinmicos producidos por las
ondas ssmicas seran diferente. Esta diferencia en los desplazamientos de
la superficie terrestre es funcin de la forma y dimensiones de la base de
cimentacin, de la distribucin de la masa y flexibilidad de la estructura,
as como, de las caractersticas del suelo en que est cimentada.
Al analizar dinmicamente una estructura en la actualidad, se acostumbra
estudiar los esfuerzos y deformaciones producidas en sus elementos
cuando se le sujete a una perturbacin de tipo ssmico aplicada en su
base. La base se considera empotrada en el suelo de cimentacin que se
supone infinitamente rgido. Con dicha suposicin la energa aplicada a la
estructura se traduce en movimiento de la misma y el nico mecanismo
de disipacin que se puede tomar en cuenta lo constituye la friccin
interna o el comportamiento inelstico del material que constituye a la
estructura. Pero est hiptesis de comportamiento infinitamente rgido no
permite representar el fenmeno que realmente ocurre en la zona de
unin del suelo con la estructura, aunque sea de fcil aplicacin en el
anlisis. El fenmeno que se presenta es: al llegar el tren de ondas a la
base de la cimentacin la estructura empieza a vibrar, sin embargo, como
el suelo que rodea a la cimentacin es deformable, parte de la energa
transmitida a la estructura es devuelta al suelo a travs de la cimentacin,
de tal manera que las ondas que siguen llegando del foco de perturbacin
se encuentran con las ondas que regresan de la estructura y la energa
que transmitira a ella, puede aumentar o disminuir debido a la
interferencia por otra parte, la energa de la estructura devuelta al suelo
puede disiparse totalmente en l, o puede reflejarse parcialmente en
estratos ms profundos y regresar disminuida a la estructura.
Como puede advertirse los fenmenos que constituyen la interaccin
dinmica entre el subsuelo y estructura son sumamente complejos y no
han sido suficientemente estudiados. Solamente se han resuelto algunos
problemas particulares. Desde luego que las ventajas que se obtendran
de un conocimiento completo del fenmeno son indiscutibles. Ello motiva
el estudio intenso de investigadores.
Modelos Matemticos.
Como ya se estableci con anterioridad que el suponer en el anlisis
dinmico de una estructura al suelo infinitamente rgido motivado por una

perturbacin horizontal en la base es solo una conveniencia del clculo

pero no refleja lo que ocurre en la realidad. Esto se ha comprobado
mediante observaciones de estructuras desplantadas en diferentes tipos
de suelo y sujetas a movimientos ssmicos.
En estos ltimos aos la preocupacin por tomar en cuenta en el diseo
ssmico las propiedades del suelo de cimentacin han ido en aumento, y
han surgido varios modelos matemticos del sistema suelo-estructura que
pretenden representar el comportamiento dinmico del sistema.
En 1948 los investigadores Crockett y Hammond estudiaron cimentaciones
de maquinaria y recomiendan considerar al suelo de cimentacin como
una masa concentrada y un resorte lineal.

WF
WS
KS

Peso de la cimentacin
Peso de una masa virtual de suelo que vibra con la cimentacin.
Constante elstica del suelo.
Dibujo

En 1954 Merrit y Housner proponen un modelo en que la base de la

estructura y suelo de cimentacin estn ligados mediante un resorte que
restringe los movimientos angulares de la base, su modelo se presenta a
continuacin:
Dibujo
W.T. Thomson en 1960 present un modelo matemtico del sistema sueloestructura en el que el suelo ofrece resistencia tanto a desplazamiento
horizontales como angulares de la base de cimentacin, pero no considera
al suelo como un resorte o un conjunto de resortes sino como un medio
elstico semiinfinito.
MS = Momento actuante debido a la inercia.
MG = Momento resistente del terreno.
KR = Constante elstico del terreno
FS = Fuerza cortante debido a la inercia.
Dibujos y ecuaciones
En 1961 D. Lycan y N.M. Newmark ofrecen una representacin del sistema
suelo-estructura consistente en masas y resortes, pero en la que el suelo
que acta con la cimentacin se considera como una masa ms. El
conjunto que se muestra en la figura no est ligada al sistema tierra, y se
sujeta a una fuerza aplicada en la masa que representa al suelo.
Dibujos

En 1965, Fleming Screwvala y Kondner proponen un modelo para tomar en

consideracin la interaccin dinmica entre el subsuelo y la cimentacin
de una estructura, que combina los propuestos por Merritt-Housner y
Lycan-Newmark.
Dibujos
Es complicado desarrollar modelos matemticos que permitan tomar en
consideracin los aspectos del complejo mecanismo de la interaccin
suelo y estructura. No obstante, existen ciertos tratamientos, empleando
tcnicas como el mtodo del elemento finito y otras formas, para este
problema. Ver trabajos de los Drs. J. Bielak y E. Rosenbluetn, (Bibliografa).
La experimentacin con modelos fsicos a escala, es una de las
herramientas principales para la resolucin del problema de interaccin
para algunos problemas, por ejemplo para el diseo de cimentaciones
sujetas a cargas dinmicas.

Las pruebas en modelos fsicos pueden dividirse en dos categoras: A)

Excitacin aplicada a travs del medio en que se apoya la cimentacin y
B) Excitacin aplicada directamente a la estructura. Para esto ltimo se
recomienda ver el trabajo presentado en el 2 Congreso de Ingeniera
Ssmica Estudios con modelos para el diseo de cimentaciones sujetas a
cargas dinmicas por L. Ayestarn, J. Elorduy y J. A. Nieto.
INTRODUCCIN
En l clculo estructural en edificios es interesante el anlisis de
construcciones cimentadas en suelos comprensibles y ubicados en zonas
ssmicas.
El anlisis tradicional de edificios en suelos compresibles se modelaban
considerando el edifico con base rgida y usando mtodos de anlisis
estticos, en edificios altos se recurra al anlisis ssmico modal sin tomar
en cuenta la interaccin suelo estructura.
Por lo expuesto anteriormente y tomando en cuenta las Normas Tcnicas
Complementarias sobre diseo por sismo vigente en el presente trabajo se
expone el anlisis ssmico de edificios altos usando el anlisis modal y
tomando en consideracin la compresibilidad del suelo. En el Capitulo II se
estudia brevemente el comportamiento del suelo en forma dinmica y
esttica y se presentan modelos para su anlisis ssmico, en el Captulo III
se muestran la forma de modelar edificios altos con base rgida y base

flexible y finalmente en el ltimo captulo se analiza un edificio de 12

niveles con dos tipos de base con el objeto de comparar ambos
comportamientos.
Introduccin a la Dinmica Estructural
En general los edificios pueden estar sujetos a diferentes cargas dinmicas
o fuerzas de excitacin que le provocan reacciones a las construcciones,
estas cargas se pueden idealizar en los siguientes tipos:
a) Cargas senoidales
b) Cargas peridicas
c) Cargas aleatorias
d) Cargas exponenciales
En las figuras 1.1 y 1.2 se muestran los tipos de cargas anteriores que son
idealizaciones de cargas provocadas por maquinaria, oleaje, sismos,
vientos y explosiones.

Fig.1.1.Vibracin peridica compleja se

caracteriza porque su presentacin matemtica puede ser difcil.

Fig.1.2.Vibracin peridica tipo

senoidal.

1.2.1. Estudio de las Vibraciones

Se definen a las vibraciones mecnicas como los movimientos de una
partcula o sistema de partculas, que oscilan alrededor de una posicin de
equilibrio.
Una estructura puede estar sujeta a vibracin libre o vibracin forzada.
Una vibracin libre es aquella que se produce bajo la ausencia de una
accin externa permanente, es decir que la estructura sigue vibrando,
cuando cesa la accin a partir de las condiciones en que la dej. Una
vibracin forzada, se tiene cuando un agente externo que vara segn una
funcin del tiempo acta permanentemente.

El nmero de coordenadas independientes necesarias para definir la

posicin de un cuerpo es llamado grado de libertad. As, un nmero infinito
de coordenadas son necesarias para describir el movimiento de un cuerpo
elstico, tales cuerpos tendrn un nmero infinito de grados de libertad. A
estos se les conoce como sistemas continuos.
Sin embargo, puede suponerse que la masa del cuerpo puede
concentrarse en un punto, entonces se tratar de un sistema con un grado
de libertad, puesto que una coordenada define la posicin de la masa. Dos
coordenadas independientes son necesarias para describir el movimiento
del sistema entonces stos poseen dos grados de libertad. Estos sistemas
tienen un nmero finito de grados de libertad y se les llama sistemas
discretos.
Se iniciar el estudio de vibraciones en sistemas discretos con un grado de
libertad, plantendolo en forma general y breve, ya sea que se trate de
una viga cargada, un eje sometido a torsin o un resorte deformado que
son movidos de su posicin de equilibrio por un agente externo, por tal
motivo las fuerzas elsticas de la pieza cuya posicin han sido perturbadas
no estarn ya en equilibrio con la carga y se producirn vibraciones.
Vibraciones Libres sin amortiguamiento.- Considrese una masa m
suspendida de un resorte sin peso dispuesto de tal manera que solo puede
tener
desplazamientos
verticales,
la
configuracin
quedar
completamente determinada por ese desplazamiento.
El desplazamiento esttico que tiene la masa, de acuerdo con la Ley de
Hooke al actuar un agente externo es:
y

pl
AE

P = ky
k
y
M

I = My

Definiremos a:
La rigidez como la fuerza, momento, etc., (Fuerza Generalizada) necesaria
para producir una deformacin angular o lineal unitaria en una estructura,
en este caso

AE
L

Y la flexibilidad, como el desplazamiento lineal o angular de la estructura

provocado por una fuerza o momento unitario en un punto, para el mismo
ejemplo

L
AE

Por lo que se advierte que el recproco de la flexibilidad es la rigidez.

y P

P
k

Al desplazar la masa y soltarla se genera una velocidad y aceleracin,

obtenindose una fuerza de inercia que por la segunda Ley de Newton es:

I my"
Estableciendo el equilibrio en la masa de las fuerzas elsticas y de inercia.

my"ky 0
Definiendo
segundo.

k
w2
m
como frecuencia circular medido en radianes por

y" w 2 y 0

(1)

La solucin de esta ecuacin diferencial homognea de coeficientes

constantes es:
y Asenwt B cos wt

(2)

En donde A y B son constantes arbitrarias. Entonces el movimiento vertical

de la masa es de naturaleza vibratoria, puesto que sen wt y cos wt son
funciones peridicas.
Para determinar las constantes de integracin es necesario conocer las
condiciones iniciales. Supngase, por ejemplo, que en el instante t = 0,
y y 0 y la velocidad y ' y 0 . De (2) obtenemos:
B y0

y0
w

Luego la ecuacin del movimiento:

y

y0
senwt y 0 cos wt
w

(3)

Otra forma de representar esta solucin, se obtiene haciendo tan

y0 w
y
v0

y0
w cos

llamando a A amplitud de vibracin y ngulo de fase.

y Asen( wt )

(4)

A y0
ecuacin est

La representacin grfica de la
t

en la figura

2
w

La frecuencia puede estar expresada en las siguientes formas:

w

k
m

2
t

g
y est

Recordando el principio de D Alambert, con la finalidad de aplicarlo en los

siguientes casos, el cual dice: Si la resultante de las fuerzas efectivas
para todos los puntos materiales de un cuerpo se invierte, suponindose
que acta sobre el cuerpo con las fuerzas externas, mantendr el cuerpo
en equilibrio. Es decir que hace posible reducir un problema de dinmica a
un problema equivalente de esttica introduciendo una fuerza que puede
hallarse basndose en el movimiento del cuerpo.
Vibracin Libre con Amortiguamiento.- En el caso anterior se encontr
que la amplitud de vibracin permanece constante, pero en cambio, la
experiencia muestra que la amplitud disminuye con el tiempo;
amortigundose gradualmente las vibraciones. Estas fuerzas de
amortiguamiento son producidas por el rozamiento en seco, resistencia del
aire o agua, rozamiento interno debido a la elasticidad imperfecta de los
materiales, friccin en remaches en estructuras metlicas, etc. El valuar el
amortiguamiento es bastante delicado. Normalmente esta resistencia se
considera proporcional a la velocidad, o sea, amortiguamiento viscoso,
siendo ste una buena aproximacin.
Para analizar este tipo de vibracin, tomamos el mismo sistema de la
figura , aplicando el Principio de D Alambert.
p + c l = ky + ry + my = 0
La resistencia est en sentido contrario a la velocidad.

y ' '

r
k
y ' y 0
m
m

y"2y ' 2 y 0

(5)
r
Definiendo 2 = m , donde es el coeficiente de amortiguamiento y se

considera constante, aunque realmente no lo es.

La ecuacin caracterstica de la ecuacin diferencial nos da las siguientes
races:
z 2 2z 2 0
2

z1 , 2 2

Observando las races vemos que se pueden presentar estos casos:

a) > , ambas races se hacen reales y negativas, la solucin no contiene ya un
factor peridico entonces la solucin no representa un movimiento vibratorio.
Cuando = se tiene al sistema con amortiguamiento crtico.
b) Si < , y definiendo a 2 2 2 ( frecuencia en sistemas con
amortiguamiento), se obtienen races complejas:
z1, 2

w1

Y la solucin general de la ecuacin diferencial es:

t

y ( Asen 1t B cos 1t )

(6)

En la cual A y B son constantes que en cada caso particular son

determinadas conforme a las condiciones iniciales. La expresin dentro del
parntesis es de la misma forma que la obtenida para vibraciones sin
amortiguamiento.
Para determinar A y B, suponemos que en el instante inicial t = 0 el
cuerpo vibrante es desplazado de su posicin de equilibrio en un valor
y y0

y que tiene una velocidad

( 6 ) y ( 7 ) se tiene:

y ' v 0 , sustituyendo estos valores en

y ' t ( Asenw1t B cos w1t ) t ( w1 A cos w1t Bw1 senw1t )

v0 y 0
B y0

w1

(7 )

Y finalmente la ecuacin de movimiento es:

y t y 0 (cos w1t senw1t ) 0 senw1t

w1
w1

(8)

Otra forma de representarla:

Introduciendo tan 1
y t

B
, de ( 6 ), se tiene:
A

A( senw1 t tan 1 cos w1t )

Desarrollando y recurriendo a una relacin trigonomtrica:

y

A
t sen( w1t 1 ) At sen( w1t 1 )
cos 1

(9)

Con los valores de A y B.

y

v y 0
y 0
w1

2
0

t sen( w1t 1 )

Representando grficamente a ( 9 ) que se tiene en la fig. En donde el

t
factor
decrece gradualmente con el tiempo y resultarn
gradualmente amortiguadas las vibraciones generadas originalmente.
En problemas prcticos, puede suponerse con suficiente precisin que
para pequeos amortiguamientos viscosos no afecta el periodo de
vibracin. La velocidad de amortiguamiento depende de la magnitud de .
En la ecuacin ( 9 ) se observa que la amplitud de vibracin disminuye
NR
: 1 , esto es, disminuye como
despus de cada ciclo, segn la razn
una progresin geomtrica. Esta razn puede utilizarse para la
determinacin experimental del coeficiente de amortiguamiento.
Algunos valores de

. (de pruebas experimentales).

Para puentes de acero

---------------------------- (0.018 0.08) w1
Para puentes de madera ---------------------------- (0.022 0.01) w1
Para puentes de concreto ---------------------------- (0.032 0.016) w1
Vibracin forzada sin amortiguamiento.- iniciamos el problema de
sistemas con un grado de libertad sujeto a una fuerza de excitacin. Esto
es, adems de la fuerza elstica y de inercia, estudiados en vibracin libre,
acta sobre la masa una fuerza peridica perturbadora. Esta carga
dinmica puede ser una funcin de tipo seno o coseno, expresndola as:

P (t ) P0 sent

En donde; es la frecuencia de la excitacin y P0 representa la amplitud.

Planteando el problema por el principio de D Alambert tenemos:
P I P (t) = 0
Substituyendo los valores correspondientes:
k y + m y - P 0 (t) = 0
y +

k
P (t )
y
m
m

llamando P0

P0
m

y" w 2 y P0 sent
( 10 )
La solucin general de la ecuacin diferencial no homognea esta
representada:
y yh y p

La solucin homognea ( y h ) ya es conocida:

y h A sen

w t B cos w t

La solucin particular ( y p ) la podemos determinar por el mtodo de los

parmetros indeterminados, suponiendo que sea del tipo:
y p y sen

En donde determinamos el parmetro y, substituyendo esta expresin en

(10). Obtenindose:

yp
O tambin, como:

P0
sen t
w 2
2

P0

P0 w 2
P0w 2 y est . w 2
k

yp

y est . w 2
sen
w2 2

y est .
sen
2
1 2
w

Finalmente la solucin general es:

y Asen w t B cos w t

y est .
sen
2
1 2
w

( 11 )

Esta ecuacin de movimientos consta de dos partes, los dos primeros

trminos representan vibraciones libres y el tercer trmino nos representa
la vibracin forzada del sistema.
Para encontrar los valores de las constantes A y B como en los casos
anteriores hacemos uso de las condiciones iniciales, cuando t = 0, el
'
desplazamiento inicial y 0 0 y la velocidad inicial y 0 v0 0 .

y ' Aw cos wt Bwsenwt

y est .
cos t
2
1 2
w

( 11 )

De ( 11 ) cuando t = 0, se tiene:
B=0
Y de ( 11) cuando t = 0, se tiene:

y est .
2
1 2
w

Y sustituyndolos valores de los coeficientes:

y est .
2
1 2
w

yc

( sen

( sen

sen w t )
w

sen w t )
w

( 12 )

Representando grficamente a la ecuacin (12), por facilidad, vamos a considerar que

0.5 , ambas funciones son senoidales, y representndolas separadamente, se obtiene la

w

fig. 12.
y

c sen t

4
y est .
3
2
y est .
3

c/2 sen t

Fig. 12
S en la ecuacin ( 12 ), no consideramos vibraciones libres, nos quedar
solamente vibracin forzada.
y y est . (

1
) sen
2
1 2
w

( 13 )

Llamndosele a esto proceso establecido.

En donde

y est .

nos representa el desplazamiento que producir la fuerza

perturbadora mxima P0 , si actuara estticamente, y el trmino entre

parntesis representa la accin dinmica de esta fuerza, se le llama

)
coeficiente dinmico
, dependiendo de la relacin w .
Entonces el desplazamiento
representado:

mximo

de

vibracin

forzada

y y est .

queda
( 14 )

Analizando el coeficiente dinmico, vemos que cuando frecuencia de la

excitacin es pequea comparada con la de vibracin libre
w
el
coeficiente dinmico es aproximadamente igual a 1
tiende a

y est .

y y y est .

, este sera un caso de vibracin muy lenta.

Pero cuando empieza a incrementarse y empieza a crecer, y en el

momento que w, el coeficiente dinmico se hace infinito y por

consiguiente la amplitud de vibracin; a este fenmeno se le conoce como

resonancia.
A continuacin se representa por medio de una grfica este
comportamiento, en donde las ordenadas son los coeficientes dinmicos y

las abcisas las relaciones w en valor absoluto.

Cuando aumenta ms all de la frecuencia de vibracin libre, el
coeficiente dinmico empieza a disminuir tendiendo a cero. Esto significa
que la fuerza perturbadora de alta frecuencia que acta sobre el cuerpo
vibrante, produce vibraciones de muy pequea amplitud, y en ciertos
casos puede considerarse que el cuerpo permanece inmvil.
Para la primera parte de la grfica o sea w , en donde la relacin es
positiva, las vibraciones forzadas y la fuerza perturbadora estn siempre
en la misma fase, esto quiere decir, que la masa vibrante llega a su
posicin ms baja en el mismo instante que la fuerza perturbadora toma
su valor mximo en sentido descendente. Y la otra parte cuando w la
diferencia de fase entre la vibracin forzada y la fuerza perturbadora se
hace igual a . Esto es, que en el momento en que la fuerza llega a su
mximo en sentido descendente, la masa llega a su posicin superior.
Debido a que el fenmeno de resonancia causa daos, estos
conocimientos nos previenen, para no estar dentro de l.
Vibracin forzada con amortiguamiento.- en los problemas prcticos
hay siempre amortiguacin y en este estudio veremos el efecto que tiene
sobre la amplitud de las vibraciones forzadas.
Nuevamente, por medio del principio de D Alambert tenemos:
P+CIP(T)=0

ky ry ' my" P0 sen

t 0

O sea:

y"2 y ' w 2 y P0 sen

( 15 )
La solucin general de la ecuacin diferencial es:
y yh y p

La solucin homognea, ya es conocida

y h t ( A cos w1 t B sen w1 t )
Y la solucin particular la determinamos en forma similar al segundo caso
estudiado, suponiendo la solucin de la siguiente forma con parmetros
indeterminados c1 y c 2 .

y p c1 cos t c 2 sen t
Para determinar estos parmetros, substituimos los valores de " y" p
derivados en ( 15 ) obtenindose:
c1 P0

c 2 P0

y sus

2
( w 2 ) 2 4 2 2
2

w2 2
( 2 w 2 ) 2 4 2 2

Entonces, la solucin general:

y t ( A cos w1 t B sen w1 t ) P0

2
2
( w2 ) 2 4 2 2

cos t P0

w2 2
( 2 w2 ) 2 4 2 2

sent

(16 )
Como en los casos anteriores, podemos expresar la ecuacin de
movimiento en otra forma.
c1 2
Introduciendo tan 1
en y p se llega a:
c2
w2 2

yp

c2
c2
(cos 1 sen t sen1 cos t )
sen( t 1 )
cos 1
cos 1

Substituyendo el valor de c 2
expresarlo en forma general.

y t ( A cos w1t B sen w1t )

y de cos 1

y est . w 2
( w 2 2 ) 2 4 2 2

yp

en

y posteriormente

sen( t 1 )

( 17 )
Para la determinacin de las constantes A y B, recurrimos a las siguientes
condiciones iniciales: t 0,

y y0 0

y ' v0 0

; substituyendo estas

condiciones en ( 16 ) y en la primera derivada de la misma, obtenemos

estas constantes expresndolas en trminos de c1 y c 2 ya conocidas.

A c1
B

c1 c 2
w1

Si las condiciones iniciales son diferentes de cero, esto es, cuando t = 0,

y y0 y

y ' v 0 , la solucin que se obtiene es:

y y 0 t (cos w1 t

c 2 ( sen t

w1

sen w1 t ) 0 t sen w1 t
w1
w1

c1 cos t

t (cos w1 t

senw1t )
w1

sen w1t )

En ( 16 ) se puede advertir fcilmente la parte correspondiente a vibracin

libre y la de vibracin forzada.
Debido al amortiguamiento, la vibracin libre desaparece despus de un
corto tiempo y quedamos solamente con las vibraciones forzadas, o sea lo
que hemos llamado proceso establecido esto es cierto en los casos
cuyas vibraciones tengan periodos amplios. Este concepto no es correcto
en sismos.
Normalmente al hacer el estudio del coeficiente dinmico, se considera un
proceso establecido. De ( 17 ) tenemos:
y y est .

w2
( w 2 2 ) 2 4 2 2

y est .

( 18 )
Llamando

(1 2 ) 2

4 2 2
w4

Como el coeficiente de amortiguamiento esta expresado

w1
2

w
2

Este coeficiente, tambin se presenta:

es el decremento logartmico.

2 2
(1 ) 2

2 2

La representacin grfica del coeficiente dinmico y la relacin , es

similar al segundo caso estudiado, la diferencia es que para este tipo de
vibracin los valores de son acotados. Y solamente cuando ,
tiende a infinito; a medida que crece , el coeficiente va disminuyendo.
En los casos de vibracin forzada, se consider que la fuerza perturbadora
era una funcin armnica. Sin embargo, puede suceder que sea una
funcin ms complicada del tiempo.
En general una fuerza perturbadora peridica de cualquier clase puede
representarse por medio de una serie de Fourier. Tambin puede ser que
esta fuerza no pueda expresarse por medio de una simple funcin
analtica del tiempo.
1.2.4. Mtodos para determinar los modos de Vibracin.
En el inciso anterior se ha expuesto el mtodo directo para la solucin del
problema de valores caractersticos, como puede observarse, cuando la
ecuacin de frecuencias es de segundo o tercer grado fcilmente se puede
resolver, pero cuando el grado es mayor que tres involucra considerable
labor no solamente en la solucin de la ecuacin sino tambin en el
desarrollo del determinante y posteriormente al encontrar los modos de
vibracin. De aqu que se hayan desarrollado mtodos iterativos y dentro
de estos, los de aproximaciones sucesivas.
En los mtodos iterativos tenemos fundamentalmente dos procedimientos
diferentes para obtener soluciones aproximadas al problema de valores
caractersticos.
En el primero la operacin iterativa bsica involucra el remplazamiento de
un vector supuesto por otro mejorado. El procedimiento conduce en
general a solamente un modo de vibracin del sistema. Las modificaciones
a este procedimiento, permite obtener otros modos.
En el segundo tipo de iteraciones la operacin bsica consiste en el
remplazamiento de la matriz cuadrada por una matriz mejorada, el
procedimiento es llamado diagonalizacin por rotaciones sucesivas. Este
conduce simultneamente a todos los modos de vibracin y sus
correspondientes frecuencias.

Existen mtodos que no obstante que tienen la forma de procesos

iterativos son mtodos directos, por ejemplo el mtodo de Lanczos o de
los n pasos, en los que la solucin exacta es proporcionada justamente
en el extremo de n iteraciones, desde luego, si todas las operaciones son
efectuadas sin redondeo.
La utilidad de estos mtodos para clculo numrico, no es comn debido
a que son muy laboriosos, de tal manera que cuando se consideran como
mtodos directos requieren ms clculo numrico que los mtodos
directos conocidos.
Ahora nos proponemos a describir brevemente algunos mtodos
numricos en forma sistemtica con el objeto de poderse programar.
Los mtodos son los siguientes:
-

Mtodo
Mtodo
Mtodo
Mtodo
Mtodo

de
de
de
de
de

Stodola Vianello
Newmark
Holzer
las matrices de transicin
Jacob.

Todos los mtodos nos proporcionan la matriz modal y la espectral o sea

los modos y frecuencias de vibracin. El mtodo de Jacob, por ejemplo,
corresponde al segundo tipo descrito anteriormente.
En la prctica comn resulta interesante, por facilidad, utilizar el mtodo
de Newmark y el mtodo de Holzer, el primero para obtener el modo
fundamental y el segundo para obtener los modos superiores, si se tiene
una calculadora manual. Sin embargo, si se hace uso de computadoras
resulta interesante el mtodo de C. G. J. Jacob, con respecto a los otros
mtodos mencionados.
Mtodo de Stodola y Vianello.- Haciendo una breve resea histrica;
diremos que la teora de los mtodos iterativos fue presentada por H. A.
Schwarz en 1885, desarrollada dicha teora por E. Picard. Los mtodos
iterativos fueron primeramente aplicados a problemas de valores
caractersticos en el estudio de problemas de inestabilidad o pandeo en
1898 por L. Vianello, tambin fueron aplicados por A. Stodola en 1904, a
problemas de velocidades crticas en vibraciones de rotores. Los primeros
trabajos trataron de valores caractersticos en sistemas continuos. Sin
embargo el procedimiento es completamente anlogo para sistemas
discretos.
En 1921, E. Pohlhausen aplic el mtodo a problemas especficos de
valores caractersticos en sistemas discretos. R. Von Mises y H. Geiringer
usaron notacin matricial para este mtodo, en 1929.

Concluyendo esta breve informacin diremos que en 1934, en el trabajo

presentado por W. J. Duncan y A. R. Collar ( a Method for the solution of
oscillation problems by matrices ) dan varias aplicaciones prcticas del
mtodo iterativo, usando notacin matricial.
El sistema de ecuaciones es arreglado de tal manera que al sustituir el
vector modal supuesto nos d el nuevo vector, que en el caso de ser el
vector modal del problema se reproducir por si mismo, en caso contrario,
esta operacin es un medio de transformacin obtenindose otro vector
mejorado. Entonces el procedimiento consiste en transformaciones
sucesivas continuas hasta obtener un vector el cual se transforma por si
mismo.
Tomando nuestro problema.
kv = w m v
En la forma:
m 1kv w2v

Llamando h m 1k
hv = w v
(1)
Haremos nuestra primera iteracin con el vector supuesto V ,
calculando .Si es el vector modal del problema, ser proporcional a ; es
decir cada elemento de ser un escalar mltiplo de los correspondientes
elementos de . El escalar mltiplo es el valor caracterstico (w). En caso
contrario,
no ser proporcional a ; no habiendo relacin entre los
correspondientes elementos de y.

Como un medio conveniente de cerciorarse si los vectores son

proporcionales, adoptamos el siguiente procedimiento: Tomaremos a
nuestro vector supuesto de tal manera que uno de sus elementos,
generalmente el mayor, sea unitario. Entonces, si los elementos de y de
son iguales, es el vector modal y el escalar o factor de dicho vector es el
valor caracterstico. Cuando no es el vector modal, como se dijo, este
proceso llega a ser una transformacin del vector supuesto por otro. La
repeticin de esta transformacin constituye el proceso iterativo
1
2
3
n
obtenindose V ,V ,V ,........V .
Ilustraremos el procedimiento con el siguiente sistema:

La ecuacin matricial de movimiento del sistema es:

5 2 0 v1
2 3 1 v w 2

2
0 1 1 v3

1 0 0 v1
0 2 0 v

2
0 0 3 v3

Al invertir la matriz m y multiplicarla por k, se obtiene:

5 2 0 v1
1 3 1 v2 w 2
2 2
0 1 1 v3
3 3

v1
v2

v3

(A)

(1)
Empezamos con un vector supuesto V
de la forma siguiente ( 1, 0, 0 ) ;
(1)
calculando V

5 2 0 1 5 1.0
1 3/2 1/2 0 1 5 0.2

0 1/3 1/3 0 0 0.0

Observamos que el vector ( 1.0, -2.0, 0.0 ) no es idntico con nuestro
(2)
vector supuesto; entonces con el nuevo vector V , hacemos la segunda
iteracin

5 2 0 1.0 5.40 1.0
1 3/2 1/2 2.0 1.296 5.40 0.24

0 1/3 1/3 0. 0. 54 0. 1

Nuevamente volvemos hacer la observacin, pero todava no hay
similitud.
Haciendo otras dos iteraciones, obtenemos:

1.00
HV (3) 5.48 0.249

0.015

1.00
HV (4) 5.49 0.252

0.016

(B)

De esta manera hallamos que la frecuencia aproximada es w = 5.49 que

2
comparndose con la real w = 5.5036 el error es pequeo para fines de
ingeniera. Esta secuencia de vectores da la aproximacin al tercer modo
que es V3 ( 1.00, - 0.252, 0.016 ).

De esta manera, el procedimiento iterativo bsico converge al modo de

vibracin de mxima frecuencia o al modo fundamental, si en vez de
invertir a la matriz M, ahora invertimos a K, sin embargo esta inversin
puede resultar muy laboriosa para sistemas con varios grados de libertad.
Para el clculo de las dems frecuencias y vectores modales usamos la
condicin de ortogonalidad del tercer modo con los dems.
V3T MVs 0


1.00, - 0.252, 0.016

1 0 0 V1s
0 2 0 V2 s 0

0 0 3 V3s

Verificando las operaciones matriciales

v1s 0.504v 2 s 0.048v3 s 0

(C

)
Donde v1s , v 2 s y
sistema

v3 s son los elementos de V1 o de

V2 , resolviendo el

5 2 0 v12 v12
2
1 3 / 2 1 / 2 v2 w v2

0 1/3 1/3 v32 v32

5v12 2v 22
v12

3
1
v 22
v 32
2
2
1
1
v 22 v32
3
3

w2
w2
w2

v12
v 22
v 32

( D)

Haciendo s = 2 en ( C ) y despejando v12

v12 0.504v22 0.048v32

(E)

Sustituyendo en ( D )
5 (0.504v 22 0.048v32 ) 2v 22 w 2 v12
3
1
v 22 v32 w 2 v 22
2
2
1
1
v 22 v32 w 2 v32
3
3

(0.504v 22 0.048v32 )

O sea
0.520v 22 0.240v32 w 2 v12
0.996v 22 0.452v32 w 2 v 22
0.333v 22 0.333v32 w 2 v32

Tomando las dos ecuaciones ltimas que constituyen el sistema reducido y

expresndolo en forma matricial:

0.996 0.452
0.333 0.333

v22 2
v w
32

v22
v
32

El sistema se resuelve nuevamente por iteraciones, empezando con un

vector supuesto ( 1, 0 )

0.996 0.452 1 0.996

1.000
0.333 0.333 0 0.333 0.996 0.334

Haciendo las comparaciones como en el caso anterior, hasta llegar a la

solucin, efectundose las siguientes iteraciones:

H V

H V

( 2)

( 4)

1.000
1.147
0.387

1.000
H 'V 1.171
0.394
( 3)

1.00
1.174

0.395

H 2 T2 H 1T2
0.8880 0.000 0.4597
H 2 0.0000 1.000 0.0000
0.4597 0.000 0.8880

3.000 0.000 0.707

0.000 1.000 0.707

0.707 0.707 2.000

0.888 0.000 0.4597

0.000 1.000 0.000

0.4597 0.000 0.888

3.366 0.325 0.0000

H 2 0.325 1.000 0.6279
0.000 0.6279 1.6339
Ntese como los valores fuera de la diagonal principal van disminuyendo
en valor absoluto.
En esta forma se continua el mtodo hasta obtener una matriz H M , en
donde todos los elementos fuera de la diagonal principal se aproximen a
cero. Reiteramos que los elementos de la diagonal principal sern los
valores caractersticos. Los resultados obtenidos continuando con el
mtodo son:

0.3250
1.0000
0.6279

3.366
0.325
0.000

0.0
0.6
1.6

3.3660
0.1703
0.2768

0.1703
2.0204
0.0000

0
0.
0.

3.3935
0.1695
0.0000

0.1695
2.0204
0.0168

0.
0
0.

3.4142
0.0000
0.0020

TM

0.0000
1.9998
0.0167
H

M 1

0.00
0.0
0.58

TM

En este caso se realizaron cinco rotaciones sucesivas. Primero se dio la

matriz por diagonal izar y posteriormente se transform H. Seleccionando
los trminos por anular, calculando el ngulo de rotacin y
substituyndolo en T. En esta forma llegamos a B.
La matriz H que se dio fue simtrica, en donde solamente se eliminaron
dos elementos
H 12 y H 23 , puesto que H 13 = 0, sin embargo se
necesitaron cinco rotaciones, no dos, esto se debe a que en una rotacin
que como sabemos tiene por objeto anular un trmino, causa que un
elemento previamente reducido a cero llega a tener un valor.
Afortunadamente cada rotacin hace que los valores mximos fuera de la
diagonal principal decrezcan en valor absoluto, y H converja en una
matriz diagonal B. As que podemos concluir que el nmero de rotaciones
ser mayor que el nmero de trminos diferentes de cero fuera de la
diagonal principal.
Los valores caractersticos sern:
w12 3.4142
w22 1.9998
w32 0.5859

Examinaremos ahora la expresin tan 2 , por medio de la cual se calculo

el ngulo de rotacin. Si las componentes v1 y v j son relacionadas con
v1 y v j , entonces la ecuacin de
tan 2 toma la forma:
tan 2

2 H ij
H ii H jj

Recordando que en trigonometra establecimos:

tan 2
2 tan

2 H ij
2 tan

2
1 tan H ii H jj
( H ii H jj ) 2 H ij (1 tan 2 )

2 H ij tan 2 2( H ii H jj ) tan 2 H ij 0

La solucin de esta ecuacin cuadrtica es:

tan

( H ii H jj )

H jj 4 H 2 ij
2

ii

2 H ij

O bien:

( H
tan
tan

tan

ii

H jj ) ( H ii H jj ) 2 4 H ij2

2 H ij ( H ii H jj )

(H

H jj ) ( H ii H jj ) 2 4 H ij2

( H ii H jj ) 4 H

( H ii H jj ) 2 ( H ii H jj ) 2 4 H ij2

ii

2 H ij ( H ii H jj ) ( H ii H jj ) 2 4 H ij2

2
ij

2 H ij
H ii H jj

( H ii H jj ) 2 4 H ij2

S H ii H jj se usar la expresin ( ). Y si H ii H jj entonces se emplear:

tan

2 H ij
H ii H jj

( H ii H jj ) 2 4 H ij2

Si arbitrariamente se elige el signo positivo en el denominador de las

ecuaciones ( ) y ( ), entonces:
tan

2 H ij

H ii H jj ( H ii H jj ) 2 4 H ij2

El signo positivo si H ii H jj . Signo negativo si H ii H jj

Se consideran tambin:

y / 2 0 / 2 .

cos (1 tan 2 ) 1 / 2
sen cos tan

Las expresiones ( ), ( ) y ( ) se usan, en lugar de la ( ). Para la

elaboracin del programa de computadora en la rotacin de ngulos.
Matriz modal o vectores caractersticos.- una vez obtenidos los valores
caractersticos por el mtodo de Jacob, ahora veremos la forma de obtener
los valores caractersticos. Los vectores modales son matrices columnas v
de la matriz dada
H. Los cuales estn asociados con los valores
2
caractersticos w .
w12 . w22 .....w N2 . Los
Existirn N valores caractersticos :
correspondientes N vectores modales v , pueden ser agrupados en una
matriz cuadrada V o sea la matriz modal.
Sea W una matriz diagonal, representando los valores caractersticos
conocidos. Entonces tenemos:
HV W 2V

Remultiplicando por la inversa de V o sea V

V 1 HV V 1 VW 2

Como V 1 V I
V 1 HV W 2

Comparando esta expresin con las obtenidas anteriormente o sea:

TT H

T
T

v W
T

B TM ....T3 T2

v
T1

T1 T2 T3 ....TM

Vemos que V es una matriz cuadrada que representa los valores

caractersticos y es igual a la multiplicacin sucesiva de matrices T
empleadas en la obtencin de valores caractersticos. Entonces la matriz
modal ser obtenida como:
V T1

T2 ..... TM

Ejemplo.- Las matrices modales para el primero y segundo ejemplo

anteriores sern:
Para el primer problema M = 1

V T1
cos
V
sen

sen
cos

Como 22.5
0.9238 0.3826
V

0.3826 0.9238
0.9238
0.3826
v1
v2

0.3826
0.9238

Para el segundo ejemplo M = 5

V T1 T2 T3 T4

T5

Tomando los valores de T anteriormente obtenidos:

cos(45)
V sen(45)

sen(45)
cos(45)

0.707
V 0.707
0.000

0.707
0.707
0.000

0.500
V 0.707
0.500

0.707
0.000
0.707

0.500
v1 0.707
0.500

0
cos(6.9)

0 ....... sen(6.9)

1
0

0.000
0.000
1.000

0.8880 0.000
0.0000 1.000

0.4597 0.000

sen(6.9)
cos(6.9)
0

0
0

0.4597
0.992

0.0000 ...... 0.120

0.000
0.8880

0.120
0.992
0.000

0
0
1

0.500
0.707
0.500
0.707
v 2 0.000
0.707

0.500
v3 0.707
0.500

Hasta aqu se han presentado un conjunto de mtodos numricos para

calcular modos y frecuencias modales, tal vez dando la impresin que son
diferentes. En realidad solamente son diferentes el mtodo de Stodola y el
mtodo de las matrices de transicin. Un tratamiento unificado de estos
mtodos numricos se presenta en el captulo 4 del libro Fundamentals of
Earthquake Engineering de N.M. Newmark y E. Rosenblueth.

1.2.3. Sistemas con varios grados de libertad.

En l inciso anterior solamente se ha estudiado vibraciones en sistemas
discretos con un grado de libertad, a continuacin se hablar directamente
de vibraciones en edificios con uno o varios grados de libertad.
Idealizacin de Estructuras.- La idealizacin matemtica de una
estructura (modelo matemtico) que tomaremos ser el formado por
masas unidas por resortes, donde la rigidez de cada resorte es la suma de
rigideces de las columnas, considerando las losas trabes rgidas, se puede
calcular las rigideces de entrepisos de los marcos por cualquiera de los
mtodos conocidos.

M3

x3

M2

x2

M1

x1

En la figura se puede advertir que las masas se concentraron en los

niveles de las losas y se consideraron a las columnas (resortes) como
elementos sin masa, es decir se est discretizando el problema. Adems
s esta considerando un solo grado de libertad por piso.
Masa equivalente.- Al discretizar un sistema es necesario encontrar
masas equivalentes de manera que las respuestas dinmicas sean lo
bastante. Esta consideracin de los sistemas de piso, solamente se
cumple para algunos edificios, segn la clasificacin que se indica
posteriormente.

Existen varias formas de determinarlas, as por ejemplo:

Por sustitucin de la frecuencia natural del sistema continuo, para esto es

necesario conocer la frecuencia natural ( wt ) del sistema considerndolo
continuo (como realmente es), entonces la masa equivalente se encuentra
con la siguiente expresin:

M eq

kA
1

2
wt AA
wt

Siendo k A o AA la rigidez o flexibilidad en el punto A.

Por el mtodo de masas nodales, esta es una forma emprica, en donde no
es necesario conocer wt , por su rapidez se utiliza en el clculo de marcos,
y consiste en concentrar en los nodos parte de la masa de columnas y
parte de la de trabes, su aproximacin a la realidad depende de la
experiencia del calculista.
En estructuras de edificios se definen como:
Sistemas estrechamente acopladas.- aquellos que pueden ser
analizados considerando nicamente traslaciones horizontales, es decir
desplazamientos por cortantes considerando las trabes de rigidez infinita.
Desde luego, esto es posible en edificios cuyos elementos para resistir
esfuerzos cortantes son marcos rgidos, donde las trabes son bastante
rgidas comparadas con las columnas.
Sistemas remotamente acopladas.- son los edificios donde las masas
se encuentran unidas por elementos cuyas deformaciones se deben
principalmente a la flexin de conjunto, es decir, los niveles giran a la vez
que se desplazan. En este tipo de sistemas es importante el efecto de una
masa alejada y debe considerarse a diferencia de los sistemas antes
mencionados, en los que solo afectan las masas adyacentes.

Deformacin por cortante.

Deformacin por flexin.

En consecuencia, el anlisis de sistemas remotamente acoplados es ms

complicado. Se pueden clasificar los edificios en este tipo cuando su
rigidez lateral est dada por muros trabajando en voladizo o bien cuando
las trabes son demasiado flexibles comparadas con las columnas.
Sistemas con tres grados de libertad de piso. Estos se presentan cuando la
excentricidad esttica sea muy grande en un nivel (o varios) en este caso
el edificio se desplazar en dos direcciones y a la vez girar por el
momento torsionante provocado por la excentricidad. En este caso el
problema se complica demasiado y no ser tratado en este trabajo.
Respuesta Modal
Inicialmente se estudiar una estructura elstica lineal con un grado de
libertad que al aplicar un movimiento en el suelo, se provocarn
vibraciones en el sistema, como se indica en la figura:

Direccin positiva del

desplazamiento
de la masa.

Posicin
Original

Direccin positiva del

desplazamiento del suelo.

El sistema elstico est formado de un conjunto de elementos rgidos

interconectados entre s mediante elementos elsticos y elementos
disipativos (amortiguadores).
Los elementos flexibles los constituyen las columnas, desprecindose su
masa, como antes se indic. El elemento disipador considera las prdidas
de energa debidas a la friccin interna del elemento flexible. Se admite
que el amortiguamiento es de tipo viscoso.
Al provocarse el desplazamiento x de la estructura, por efecto de ste
aparece un cortante V.

V=kx
Siendo:
V = Fuerza cortante
k = rigidez al corte de la columna
x = desplazamiento relativo de la masa, respecto a la base de la
columna (adelante solo hablaremos de desplazamiento relativo).
Observe que en este captulo se usar el eje x en lugar del y.
Recordando que el amortiguamiento
proporcional a la velocidad

lo

habamos

considerando

r = c x (antes c = Ry ).
Estableciendo la ecuacin de equilibrio:
mU + cx + kx = 0

pero

mx + cx + kx = - mv

llamando

U=v+x

Luego:

mx + cx + kx = f (t)

f (t) = -m v
(19)

f (t) = perturbacin externa

Esta ecuacin es anloga a la (15) obtenida para un caso general.

Entonces se trata de un caso de vibracin forzada que depender del
desplazamiento del suelo.
Debido a que la estructura est sujeta a un movimiento en su base y no a
una fuerza, el esfuerzo mximo que soporta es funcin de su rigidez, as
como de su periodo de vibracin.
En general en cuanto mayor es la rigidez de las columnas mayor ser el
esfuerzo a que estn sujetas y menor sus desplazamientos relativos para
un movimiento de suelo determinado.
S conocemos una grfica de desplazamiento-tiempo podemos determinar
el mximo de desplazamiento relativo y por lo tanto el cortante mximo,
as como
la
correspondiente
aceleracin
de
la
masa,
ver
grficas.................,

Sistemas de varios grados de libertad.

En este tipo de sistemas deben formularse tantas ecuaciones de
movimiento como coordenadas independientes haya. Consideremos el
siguiente sistema de tres grados de libertad interviniendo las fuerzas
elsticas, de inercia, de amortiguamiento y excitacin.
Planteando las ecuaciones de equilibrio.

mx1 " cx1 '5kx1

2 kx 2

2 mx 2 "2cx 2 '2cx 3 '2 kx1

3mx 3 "2cx 3 '2cx 2 ' kx 2

Usando notacin matricial para facilitar la discusin y porque adems nos

proporciona un sistema til para disposicin de la computacin resulta.

2m

3m

Denominando:

x
"

Vector de aceleraciones

x1'

x' x 2'
x3'

x1
x x 2
x3

Vector de velocidades

Vector de desplazamientos

"
x
1
"
x 2
"
x 3

 f1
f f 2
f 3

Vector de fuerzas

m
m 0
0

c
c 0
0

0
2m
0

0
2c
2c

5k
k 2k
0

0
0
3m

Matriz de masas (diagonal)

0
2c

2k
3k
k

Matriz de amortiguamiento

2c

k
k

Matriz de rigideces

De tal manera que la ecuacin matricial de movimiento del sistema:

(20)

Mx"Cx ' Kx F

Esta expresin es aplicable a cualquier sistema elstico independiente de

la forma en que sus elementos estn ligados entre s y de las coordenadas
que se utilicen para establecerla.
En las matrices M, C y K, los elementos fuera de la diagonal principal se
llaman de acoplamiento. Si las matrices de coeficientes son no-acoplados
(matriz diagonal), las ecuaciones de movimiento se plantean
independientes para cada desplazamiento x j (t ) generalizado. El tipo de
acoplamiento depende del sistema de referencia elegido.
Si en la ecuacin (20) f i es igual a una matriz nula, el sistema presenta
un estado de vibracin libre, en caso de que exista, como vimos
anteriormente, se tratar de vibracin forzada.
En forma similar a la rigidez angular y lineal, existen valores que
caracterizan el comportamiento dinmico de una estructura, stas son las
frecuencias naturales circulares
w. El problema asociado a la
determinacin de estos valores se le denomina problemas de valores
caractersticos. Y estos vectores nos determinarn las configuraciones de

la estructura. La obtencin de ellos requiere conocer las ecuaciones de

movimiento.
Para continuar el estudio, consideremos un sistema con vibracin libre
con n grados de libertad, sin amortiguamiento (que es una de las
hiptesis para el anlisis dinmico elstico modal para estructuras),
quedndonos la ecuacin de movimiento as:

(21)

Mx"kx 0

La solucin de este sistema de n ecuaciones diferenciales es,

segn vimos al estudiar sistemas con un grado de libertad, expresin (4).
x V sen

( wt )

(22)

En donde el vector x nos representa el desplazamiento dinmico en cada

una de las masas y.

V v jk

Es el vector de amplitudes de desplazamientos.

Siendo esta la solucin, debe satisfacer a la ecuacin (21), sustituyndola:

MVw 2 sen ( wt ) kV sen ( wt ) 0
MVw 2 kV 0

o bien

kV w 2 MV

(23)

Esta es la ecuacin del problema de valores caractersticos.

No debe perderse de vista que esta ecuacin matricial nos representa un
sistema de n ecuaciones algebraicas homogneas en V.
Para que tenga solucin diferente de cero la ecuacin, (23) su
determinante debe ser igual a cero.

k w M V 0
2

k w2 M 0

(24)

Siendo el determinante nulo esto implica que alguna de las ecuaciones es

combinacin lineal de las otras y por lo tanto los valores de V que se
encuentran, estn en funcin de uno de ellos.
2

Esto es, que se encontrar una ecuacin poli nmica de grado n en w ,

lo cual deber tener n races. Obtenindose para cada uno de estos

2
2
valores caractersticos de w una solucin del sistema, as, para w1 le

corresponde un vector caracterstico V1 a este se le denomina, en el caso

de vibraciones, modo natural de vibracin o vector modal, el cual es un
vector cuyas componentes son las amplitudes de los desplazamientos
reales de cada masa.
De lo anterior podemos definir al modo de vibracin como una
configuracin de la forma de vibrar de la estructura bajo una frecuencia
dada. Tenindose n modos de vibracin.
El conjunto de n vectores modales forma una matriz modal.

V v jk

v11
v
21

......

v n1

v12
v 22
.....
vn 2

v13 ...................v1n
v 23 ...................v 2 n
.............................

v n 3 ...................v nn

(25)
La primera columna nos representa el primer modo de vibracin, la
siguiente el segundo modo, y as sucesivamente hasta el ensimo modo.
V jk

se interpreta como la amplitud de desplazamiento de la masa j en el

modo k.
2

Al primer modo, para el valor menor de w , se le llama modo fundamental

y a los dems se les llama modos superiores.
2

Los n valores de w se les puede colocar como elementos diagonales de

una matriz diagonal que se le llama matriz espectral.

w12 0........0

0
w22 .......
2
2

w
..... ............

2
0 0.......wn

(26)

2
Desde luego, cada valor caracterstico wr y su respectivo modo v r deben
satisfacer la ecuacin diferencial matricial.

K v r wr2 M v r
Considerando el conjunto de valores caractersticos y modos asociados se
satisface:

K V M V 2
Propiedades de los modos naturales.

(27)

Ortogonormalidad de los modos con respecto a las masas y a las

rigideces.
2
Si tenemos dos valores caractersticos wr

corresponde un modo caracterstico, v r

matricial.

k vr wr2 M vr

v s satisfaciendo la ecuacin

( A)
( wr w s )

k v s ws2 M v s

( B)

Transponiendo la ecuacin (A)

(k vr ) t wr2 ( M vr ) t
Multiplicando por v s :

vrt k t v s wr2 vrt M t v s

Siendo M y K matrices simtricas:
kt k
v rt k v s wr2 v rt M v s

Mt M
(C )

t
Ahora multiplicamos a (B) por v r :

v rt k v s ws2 v rt M v s
Efectuando la resta (D) de (C) :

0 ( wr2 ws2 ) v rt M v s

ws2 , a cada uno le

( D)

2
2
Como las frecuencias son diferentes: wr ws 0

v rt M v s 0

(E)

Esto nos indica que los modos son ortogonales con respecto a la funcin
de peso M.

Sustituyendo esta propiedad en (D) (C), obtenemos:

v rt k v s 0

(F)

Ortogonalidad de los modos con respecto a la funcin de peso k.

Recurdese que en algebra vectorial se define que la ortogonalidad
ordinaria de dos vectores esta dada por:
A B 0

At B 0

o bien

1) Normalizacin de modos.
Se define la norma de un vector como:
2

A2 A A a x a y az

En forma matricial:
At A A

En trminos de la funcin de peso M, la norma del vector esta dada por:

v rt M vr L2r

(G)

Norma del vector modal v r con respecto a la funcin de peso M.

Podemos expresar las condiciones (E) y (G) simultneamente por medio de
la delta de Kronecker, el cual adquiere solamente dos valores en funcin
de sus ndices, es decir:
1 j

v rt M v s rs

si

1 j

si

1 j

L2

Se dice que el vector v r

Entonces:
v rt
Lr

Mv s
rs
Ls

esta normalizado si tiene norma unitaria.

llamando

vr
r
Lr

tr M s rs
En lo cual r y s son vectores unitarios.
Expresando el conjunto de ecuaciones por medio de matrices
Vt MV I

La cual establece la ortonormalizacin de los modos Vr .

La matriz V contiene todos los vectores modales del problema.

V 1

2 ........... n

11

1 12
.....

1n

Obsrvese que solo se ha establecido las condiciones de ortogonalidad de

los vectores modales, cuando las frecuencias son todas distintas. Pero
puede suceder que una raz o varias de la ecuacin caracterstica se
repitan, en este caso es relativamente fcil encontrar que esos vectores
satisfagan el sistema homogneo de ecuaciones pero no es inmediato al
que estos verifiquen entre s las condiciones de ortogonalidad lo cual, por
otro lado, es conveniente para el estudio de vibraciones.
Desarrollo en serie.
Si tenemos n vectores v ortogonales entre si, estos n vectores
pueden ser una base en un sistema de referencia en un espacio de n
dimensiones. En particular si los v r son los modos, esto expresa la
posibilidad de poner cualquier configuracin en trminos de los vectores
modales:
y r 1 Ar v r
n

r 1, 2 , ..........., n

(28)

Los coeficientes
los llamaremos
para evaluarlos,
en un desarrollo

Ar ; son los coeficientes de una serie finita de Fourier aqu

coeficientes de participacin o multiplicadores modales

se hace en la misma forma que los coeficientes de Fourier
infinito. Multiplicando (28) por:

v st m
v st my r 1 Ar v st m v r
n

t
Por condicin de ortogonalidad v r m v r 0 , entonces al desarrollar la serie
solamente es diferente de cero cuando r s obtenindose:

v rt m y
Ar t
vr m vr

(29)

Substituyendo en (28):

y r 1
n

v rt m y
v rt m v r

vr

t
S los modos estn normalizados y son ortogonales; v r m v r 1

y r 1 v rt m y v r
n

De esta manera hemos logrado establecer una expresin que nos indica
que una forma cualquiera de vibrar de una estructura se puede poner en
funcin de los modos de vibracin naturales. Entonces s
cualquiera:

y r 1
n

y vt

un modo

v rt M vt v r
v rt M v r

(31)
Estableciendo la solucin general para el problema de vibraciones libres
sin amortiguamiento, ya iniciado:
M x+ k x = 0
La solucin general de este sistema requiere el establecimiento de 2n
constantes arbitrarias las cuales por ejemplo, pueden determinarse a
partir de las 2n condiciones iniciales:
x = x (0)

desplazamiento inicial

x = x (0)

velocidad inicial

para t = 0

La solucin puede verificarse por substitucin que corresponde a (vase

ecuacin (2) ):
x r 1 v r
n

( Ar cos wr t Br sen wr t )

(32)
Los vectores v r son los modos de vibracin que satisfacen la ecuacin.

k v r wr2 M v r
Esto es, que hay necesidad de resolver previamente el problema de
valores caractersticos dado por (23).
Para determinar las constantes Ar y Br , substituimos las condiciones
iniciales en (32) y su primera derivada obtenindose:

x (0)

r 1

x ' ( 0)

n
r 1

Ar v r
B r wr v r

De tal manera, que de acuerdo con el proceso de (28) llegamos a:

Ar

v rt M x (0)
v rt M v r

v rt M x' (0)
Br
wr v rt M v r

Finalmente substituimos en (32):

x (t )

v r v rt M
r 1 v t Mv
r
r

x (0) cos wr t

1
x' (0)
wr

senwr t

(33)
As, que por medio de esta expresin podemos obtener el estado de
desplazamientos y cortantes de la estructura para cualquier tiempo.
Vibracin Forzada sin Amortiguamiento.- para este caso la ecuacin
matricial es:
Mx+kx=f
En donde f = f (t)

(34)

La solucin de este problema sin necesidad de recurrir a soluciones

complejas, se obtiene desacoplando el sistema de ecuaciones.
Considerando la matriz modal ortonormalizada respecto a M y adems la
matriz espectral.

v11
v
21

v12 ............v1n
v 22 ............v 2 n

..... ......................

v n1 v n 2 ............v nn

;
2

0...............0

0
0

w 2 ............0
0 2............wn2

w12

2
r

Habamos visto que este conjunto de valores caractersticos y modos

satisfacen (27):

K V M V 2
Por ser v r ortogonales y unitarios se verifica.
v rt M v s

rs

El conjunto de estas condiciones de ortonormalidad se expresa:

V tM V I
t
Multiplicando (27) por v :

V t K V vt M V 2
V t k v 2

(35)
Introduciendo una transformacin lineal del tipo:
X=VX
(36)
En (34):
M V x+ k v x = f (t)
t
Multiplicando esta expresin por v :

V t M V x' ' ' V t k v x' v t f (t )

Substituyendo el valor de (35) y v t M v 1
x' ' ' 2 x' v t f (t )

(37)

Como es una matriz diagonal, este sistema esta formado de

ecuaciones diferenciales independientes en
X, es decir, estn
desacoplados. Entonces la solucin de esta ecuacin matricial nos da el
valor de X que multiplicndola por la matriz modal V (36) se obtendr el
vector de desplazamientos X.
2

Aplicacin.- Para la estructura esquematizada en la figura, se calcularn

las frecuencias y periodos naturales, los modos ortonormalizados y el
estado de desplazamientos para t = 0.2, suponiendo las condiciones
iniciales (t = 0) siguientes. Considrese vibracin libre.
0.5
x (o) 0.5
0.5

cm.

0
x' (o) 0
0

VIBRACIONES MECANICAS.
Mtodos numricos para la determinacin de modos de vibracin.
En el capitulo anterior se ha expuesto el mtodo directo para la solucin
del problema de valores caractersticos, como puede observarse, cuando
la ecuacin de frecuencias es de segundo o tercer grado fcilmente se
puede resolver, pero cuando el grado es mayor que tres involucra
considerable labor no solamente en la solucin de la ecuacin sino
tambin en el desarrollo del determinante y posteriormente al encontrar
los modos de vibracin. De aqu que se hayan desarrollado mtodos
iterativos y dentro de estos, los de aproximaciones sucesivas.
En los mtodos iterativos y dentro de estos tenemos fundamentalmente
dos procedimientos diferentes para obtener soluciones aproximadas al
problema de valores caractersticos.
En el primero la operacin iterativa bsica involucra el remplazamiento de
un vector supuesto por otro mejorado. El procedimiento conduce en
general a solamente un modo de vibracin del sistema. Las modificaciones
a este procedimiento, permite obtener otros modos.
En el segundo tipo de iteraciones la operacin bsica consiste en el
remplazamiento de la matriz cuadrada por una matriz mejorada, el
procedimiento es llamado, Diagonalizacin por rotaciones sucesivas.

Este conduce simultneamente a todos los modos de vibracin y sus

correspondientes frecuencias.
Existen mtodos que no obstante que tienen la forma de procesos
iterativos son mtodos directos, por ejemplo el mtodo de Lanczos o de
los n pasos, en los que la solucin exacta es proporcionada justamente
el extremo de n iteraciones, desde luego, si todas las operaciones son
efectuadas sin redondeo.
La utilidad de estos mtodos para clculo numrico, no es comn debido a
que son muy laboriosos, de tal manera que cuando se consideran como
mtodos directos conocidos.
Ahora nos proponemos a describir brevemente algunos mtodos
numricos en forma sistemtica con el objeto de poderse programar.
Presentando cinco mtodos:
Mtodo de Stodola-Vianello
Mtodo de Newmark
Mtodo de Holzer
Mtodo de las matrices de transicin
Mtodo de Jacobi
Todos los mtodos nos proporcionan la matriz modal y la espectral o sea
los modos y frecuencias de vibracin. El mtodo de Jacobi, por ejemplo,
corresponde al segundo tipo descrito anteriormente.
En la practica comn resulta interesante, por facilidad, utilizar el mtodo
de Newmark y el mtodo de Holzer, el primero para obtener el modo
fundamental y el segundo para obtener los modos superiores, si se tiene
una calculadora manual, sin embargo, si se hace uso de la computadora
resulta interesante el mtodo de C. G. J. Jacobi respecto a los otros
mtodos mencionados.
Mtodo de Stodola y Vianello
Haciendo una breve resea histrica; diremos que la teora de los mtodos
iterativos fue presentada por H. A. Schwarz en 1885, desarrollada dicha
teora por E. Picard.

Los mtodos iterativos fueron primeramente aplicados a problemas de

valores caractersticos en el estudio de problemas de inestabilidad o
pandeo en 1898 por L. Vianello, tambin fueron aplicados por A. Stodola
en 1904, a problemas de velocidades crticas en vibraciones de rotores.
Los primeros trabajos trataron de valores caractersticos en sistemas
continuos. Sin embargo el procedimiento es completamente anlogo para
sistemas discretos.
En 1921, E. Pohlhausen aplic el mtodo a problemas especficos de
valores caractersticos en sistemas discretos. R. Von Mises y H. Geiringer
usaron notacin matricial para este mtodo, en 1929.
Concluyendo esta breve informacin diremos que en 1934, en el trabajo
presentado por W. J. Duncan y A. R. Collar (a Method for the solution of
oscillation problems by matrices) dan varias aplicaciones prcticas del
mtodo iterativo, usando notacin matricial.
El sistema de ecuaciones es arreglado de tal manera que al sustituir el
vector modal supuesto nos d el nuevo vector, que en el caso de ser el
vector modal del problema se reproducir por si mismo, en caso contrario,
esta operacin es un medio de transformacin obtenindose otro vector
mejorado. Entonces el procedimiento consiste en transformaciones
sucesivas continuas hasta obtener un vector el cual se transforma por si
mismo.
Tomando nuestro problema.
Kv 2 mv

En la forma:
m 1 Kv 2 v

Llamando h m 1 K
hv 2 v
Haremos nuestra primera iteracin con el vector supuesto v (1) , calculando
v (1) h . Si v (1) es el vector modal del problema, v (1) h ser proporcional a
v (1) ; es decir cada elemento de v (1) h ser un escalar mltiplo de los
correspondientes elementos de v (1) . El escalar mltiplo es el valor
caracterstico 2 . En caso contrario,
no ser proporcional a v (1) ; no
habiendo relacin entre los correspondientes elementos de v (1) y v (1) h .
Como un medio conveniente de cerciorarse si los vectores son
proporcionales, adoptamos el siguiente procedimiento: Tomaremos a
nuestro vector supuesto de tal manera que uno de sus elementos,
generalmente el mayor, sea unitario. Entonces, si los elementos de v (1) y
de v (1) h son iguales, v (1) es el vector modal y el escalar o factor de
dicho vector es el valor caracterstico.

Cuando v (1) no es el vector modal, como se dijo, este proceso llega a ser
una transformacin del vector supuesto por otro. La repeticin de esta
transformacin constituye el proceso iterativo obtenindose
v 1 , v 2 , v 3 ,........v n .

Ilustraremos el procedimiento con el siguiente sistema:

x1
x2
x3
3K

2K

2m

3m

La ecuacin matricial de movimiento del sistema es:

5 2 0
2 3 1

0 1 1

v1
v 2
2
v3

1 0 0 v1
0 2 0 v

2
0 0 3 v3

Al invertir la matriz m y multiplicarla por K, se obtiene:

5 2 0 v1
1 3 1 v2 2
2 2
0 1 1 v3
3 3

v1
v2

v3

(A)

Empezamos con un vector supuesto v (1) de la forma siguiente ( 1, 0, 0 ) ;

calculando v (1) h .

5 2 0 1 5 1.0
1 3/2 1/2 0 1 5 0.2

0 1/3 1/3 0 0 0.0

Observamos que el vector (1.0, -2.0, 0.0) no es idntico con nuestro vector
supuesto; entonces con el nuevo vector v ( 2 ) , hacemos la segunda
iteracin

5 2 0 1.0 5.40 1.0
1 3/2 1/2 2.0 1.296 5.40 0.24

0 1/3 1/3 0. 0. 54 0. 1

Nuevamente volvemos hacer la observacin, pero todava no hay
similitud.

Haciendo otras dos iteraciones, obtenemos:

1.00
hv (3) 5.48 0.249

0.015

(B)

1.00
hv (4) 5.49 0.252

0.016

De esta manera hallamos que la frecuencia aproximada es 2 = 5.49 que

comparndose con la real 2 = 5.5036 el error es pequeo para fines de
ingeniera. Esta secuencia de vectores da la aproximacin al tercer modo
que es v3 (1.00, - 0.252, 0.016).
De esta manera, el procedimiento iterativo bsico converge al modo de
vibracin de mxima frecuencia o al modo fundamental, si en vez de
invertir a la matriz m, ahora invertimos a K, sin embargo esta inversin
puede resultar muy laboriosa para sistemas con varios grados de libertad.
Para el clculo de las dems frecuencias y vectores modales usamos la
condicin de ortogonalidad del tercer modo con los dems.
v3T mv s 0


1.00, - 0.252, 0.016

1 0 0 v1s
0 2 0 v2 s 0

0 0 3 v3s

Verificando las operaciones matriciales

v1s 0.504v 2 s 0.048v3 s 0
(C)

Donde v1s , v 2 s y
sistema

v3 s son los elementos de v1 o de

v 2 , resolviendo el

5 2 0 v12 v12
2
1 3 / 2 1 / 2 v2 v2

0 1/3 1/3 v32 v32

5v12 2v 22
v12

(D)

3
1
v 22
v32
2
2
1
1
v 22 v 32
3
3

2
2

v12
v 22
v32

Haciendo s = 2 en (C) y despejando v12

v12 0.504v 22 0.048v 32
(E)

Sustituyendo en (D)
5 (0.504v 22 0.048v32 ) 2v 22 2 v12
3
1
v 22 v32 2 v 22
2
2
1
1
v 22 v32 2 v32
3
3

(0.504v 22 0.048v32 )

O sea

0.520v 22 0.240v32 2 v12

0.996v 22 0.452v32 2 v 22
0.333v 22 0.333v32 2 v32

Tomando las dos ecuaciones ltimas que constituyen el sistema reducido y

expresndolo en forma matricial:

0.996 0.452
0.333 0.333

v22 2

v
32

v22
v
32

El sistema se resuelve nuevamente por iteraciones, empezando con un

vector supuesto (1, 0)

 0.996 0.452 1 0.996

1.000
0.333 0.333 0 0.333 0.996 0.334

Haciendo las comparaciones como en el caso anterior, hasta llegar a la
solucin, efectundose las siguientes iteraciones:

1.000
h' v 1.147
0.387

1.000
h' v 1.171
0.394

( 2)

( 3)

1.00
h' v 1.174

0.395
( 4)

De esta manera se ha obtenido v22 y v32. El valor de v12 se calcula de (c).

v12 0.504 x1.00 x0.048(0.395) 0.524

Entonces:

1.174
2

0.524
v 2 1.000
0.395

Y finalmente el tercer vector modal debe satisfacer las 2 condiciones de

ortogonalidad y a la ecuacin matricial. Esto es:
v11 0.504v 21 0.048v31 0

1 0 0
0.524 1.000 0.395 0 2 0
0 0 2

Obtenida anteriormente para s

v11
v
21
v31

0.524v11 2.000v 21 1.185v31 0

Resolviendo estas 2 ecuaciones para

obtenemos:

v11 y v 21

en trminos de

v31

v11 0.221v31
v 21 0.534v31

De tal manera que el modo fundamental es:

0.221
v1 0.534
1.000

Es importante observar que generalmente en las aplicaciones reales lo

que interesa es calcular las frecuencias mas bajas (correspondientes a los
2 o 3 primeros modos). Y en esta forma tendran que calcularse
primeramente los modos superiores. La solucin seria, utilizar la matriz de
flexibilidades, ya mencionada. A continuacin damos un mtodo que
converge a la frecuencia mas baja.
Mtodo de Jacobi
Anteriormente se haba mencionado un mtodo para obtener los valores y
vectores caractersticos, la aplicacin de dicho mtodo nos conduce a la
matriz Modal o sea a todos los modos de vibracin del sistema por medio
de rotaciones sucesivas. El mtodo de Jacobi, es pues, un mtodo de
diagonalizacin por medio de rotaciones sucesivas. Antes de exponerlo se
explicara el concepto de transformacin de coordenadas.
Transformacin de coordenadas.- consideremos dos vectores v y v ; cuyas
componentes son v1 , v 2 y v1 , v 2 que forman dos sistemas de ejes en el
plano, con el mismo origen como se indica en la figura.

v2

v1

v2
v1

A
v2

v2

v cos

v1 sen

v 2 sen

v1

v1

v1 cos

El vector OA tiene componentes sobre los dos sistemas de ejes,

observando la figura podemos establecer las siguientes relaciones:
v1 v1 cos v 2 sen
v 2 v1 sen v 2 cos

En forma Matricial
v1 cos
v sen
2

sen
cos

v1
v
2

O simplemente
v Tv

(1)

Matriz de transformacin. La transformacin de v se puede tener

T
aplicando la matriz T , obtenindose el vector v .
Ahora recordemos que el problema de valores caractersticos tiene la
siguiente forma:
kv 2 mv

(2)

O bien
Hv 2 v

Ver mtodo de Stodola-Vianello

(2)

Sustituyendo el valor de v:
HTv 2T v

Premultiplicando por la transpuesta de T

T T HTv 2T T T v
T T HTv 2 v

(3)

Porque T T T 1 esto es:

cos
sen

sen cos
cos sen

sen
1

cos
0

0
1

(4)

Mtodo.- El objeto bsico del mtodo de Jacobi. Es transformar una

ecuacin de la forma:
Hv 2 v

(2)

A otra de la forma
Bv 2 v

(3)

Donde B es una matriz diagonal o sea que todos los trminos fuera de la
diagonal principal son ceros. Entonces como se puede observar de (3) y(3)
el valor de B ser:

B TTH T

(5)

Para este caso tenemos.

cos
sen

sen
cos

h11
h
12

h21
h22

cos
sen

sen

cos

h11 cos 2 2h12 sen cos h22 sen 2

h12 cos 2 sen 2 sen cos (h22 h11 )

2
2
h11 cos 2 2h12 sen cos h22 sen 2
h12 cos sen cos sen (h22 h11 )

Y hagamos que los trminos fuera de la diagonal principal valgan cero,

h12 cos 2 sen 2 cos sen (h22 h11 ) 0

Revisando nuestra trigonometra, recordamos que

que 2 cos sen sen 2 , luego:

cos

sen 2 cos 2 y

1
sen 2 ( h22 h11 ) 0
2
Dividiendo entre cos 2 llegamos a obtener.
h12 cos 2

Tan 2

2h12
2h12

h22 h11
h11 h22

(6)

Obteniendo el ngulo para la cual esta expresin se anula. Como la

matriz h es simtrica ambos trminos fuera de la diagonal principal se
convierten en ceros. Obtenindose la matriz diagonalizada B.

h11 cos 2 2h12 sen cos h22 sen 2

0
B

2
2
0
h
sen

2
h
sen

cos

h
cos

11
12
22

Inmediatamente obtenemos los valores caractersticos (frecuencias de

vibracin), para nuestro problema) B11 y B22

b11 0
0 b

22

v1 12 0
v
2
0 2

v1
v
2

B11 1
B22 22

Ejemplo.- fcilmente podemos comprender lo que anteriormente se

explico en el siguiente ejemplo.
Sea el problema Hv 2 Iv
T T HTv 2 v
4

Dando H

1 2
Entonces:
4 cos 2 2sen cos 2 sen 2
0

2
4 sen 2 sen cos 2 cos
2

Siempre que

Tan 2

2(1)
1
42

22.5

Substituyendo el valor de en B se tiene:

4.414
0

0
1.586

Luego:

12 4.414

22 1.586

v1

A continuacin se har extensivo el proceso del mtodo para una h y

tercer orden.

v de

v2
v2

v1
v v 2
v3

v1

v1
v v 2
v3

v1

O
v2
v3

v3

v1

v3

En este caso tenemos:

v1 v1 cos v2 sen
v 2 v1 sen v 2 cos
v 3 v3

En forma matricial. Se tiene la siguiente transformacin.

v1
v
2
v3

cos sen 0
sen cos 0
0
0
1

v Tv

v1
v 2
v3
(7)

Siendo T la matriz de transformacin hemos establecido que

T T HTv 2 v sirve para eliminar uno de los elementos de la diagonal
principal, aprovechando este conocimiento podemos lograr que los
elementos fuera de la diagonal principal se anulen. Obviamente que
realizando unas de las Rotaciones no ser suficiente para lograr que
todos los trminos deseados sean ceros. Por lo que es necesario hacer m
rotaciones en el plano hasta obtener la transformacin de H , esto es:

Tm ...T3

T2

T1

H T1 T2 T3 ...Tm

(8)

Convergiendo en una matriz B, que como hemos establecido ser una

matriz diagonal. Los elementos de la diagonal principal sern los valores
caractersticos .
En el proceso es conveniente tomar en cuenta la siguiente regla:
Seleccionar el elemento mayor en valor absoluto fuera de la diagonal
principal d la matriz H , este elemento ser el que se elimine. Para
anularlo se arreglar le matriz T , moviendo el elemento sen tomando
la posicin del elemento hij por eliminar. Naturalmente que al mover un
rengln de la matriz T se mover una columna en el mismo orden
(propiedad de los determinantes). Esto se ilustra el siguiente ejemplo.

Ejemplo.- Tenemos una matriz simtrica H.

2
H 1
0

1
2
1

0
1
2

Seleccionaremos el trmino h12 como el elemento mayor fuera de la

diagonal principal, su valor absoluto es 1. Se puedo haber seleccionado el
trmino h23 , puesto que h23 h12 , esto no afecta el resultado de
diagonalizacin. En este caso T1 no tiene ningn movimiento de renglones
y columnas, esto es:
cos
T1 sen
0

sen
cos
0

Calculo de

0
1

2h12
2 * 1

h11 h22
22

tan 2

45
cos(45)
T1 sen(45)

sen(45)
cos( 45)
0

La matriz H1 ser:

2
0
2

2
0

1
0

2
2
0

H 1 T1 H T1
T

H1

2 2

2 2

2
2

3.000
H 1 0.000
0.707

2 1 0
1 2 1
0 1 2

2
2
2

2
2
2

2

2
2
2
1
2
0

0.000
0.707
1.000 0.707
0.707 2.000
2

2
2

2
0
2
2
0

2
2

2
0

En esta forma se ha obtenido H1, repitiendo el proceso,

supondremos que H1 es H, o sea calcularemos H2 a partir de H1.

ahora

Seleccionando h13 = 0.707 como el trmino mayor en valor absoluto, o sea

el trmino que debemos anular con el siguiente ngulo .

tan 2

2h13
2 * (0.707)

1.4142
h11 h13
32

27.37

Considerando la regla establecida es necesario mover a -sen a la

posicin de T13, primero movemos la columna de la posicin dos a la tres y
despus el rengln de la posicin dos a la tres, esto es:

 cos 27.37
0

T2

sen 27.37

sen 27.37

cos 27.37

o bien
0.8880
T2 0.0000
0.4597

0
1
0

0.4597
0.0000
0.8880

H 2 T2 H 1T2
0.8880 0.000 0.4597
H 2 0.0000 1.000 0.0000
0.4597 0.000 0.8880

3.000 0.000 0.707

0.000 1.000 0.707

0.707 0.707 2.000

0.888 0.000 0.4597

0.000 1.000 0.000

0.4597 0.000 0.888

3.366 0.325 0.0000

H 2 0.325 1.000 0.6279
0.000 0.6279 1.6339

Ntese como los valores fuera de la diagonal principal van disminuyendo

en valor absoluto.

En esta forma se contina el mtodo hasta obtener una matriz H m , en

donde todos los elementos fuera de la diagonal principal se aproximen a
cero.

Reiteramos que los elementos de la diagonal principal sern los valores

caractersticos. Los resultados obtenidos continuando con el mtodo
son:

3.366
0.325
0.000

0.3250
1.0000
0.6279

0.0
0.6
1.6

3.3660
0.1703
0.2768

0.1703
2.0204
0.0000

0
0.
0.

3.3935
0.1695
0.0000

0.1695
2.0204
0.0168

0.
0
0.

3.4142
0.0000
0.0020

Tm

0.0000
1.9998
0.0167
H

m 1

0.00
0.0
0.58

Tm

En este caso se realizaron cinco rotaciones sucesivas. Primero se dio la

matriz por diagonalizar y posteriormente se transform H. Seleccionando
los trminos por anular, calculando el ngulo de rotacin y
substituyndolo en T. En esta forma llegamos a B.
La matriz H que se dio fue simtrica, en donde solamente se eliminaron
dos elementos
H 12 y H 23 , puesto que H 13 = 0, sin embargo se
necesitaron cinco rotaciones, no dos, esto se debe a que en una rotacin
que como sabemos tiene por objeto anular un trmino, causa que un
elemento previamente reducido a cero llega a tener un valor.
Afortunadamente cada rotacin hace que los valores mximos fuera de la
diagonal principal decrezcan en valor absoluto, y H converja en una
matriz diagonal B. As que podemos concluir que el nmero de rotaciones
ser mayor que el nmero de trminos diferentes de cero fuera de la
diagonal principal.

Los valores caractersticos sern:

12 3.4142
22 1.9998
32 0.5859

Examinaremos ahora la expresin Tan 2 , por medio de la cual se calculo

el ngulo de rotacin.
Si las componentes v1 y v j son relacionadas con v1 y v j , entonces la
ecuacin de tan 2 toma la forma:
tan 2

2hij

(9)

hii h jj

Recordando que en trigonometra establecimos:

tan 2
2 tan

2hij
2 tan

2
1 tan hii h jj
(hii h jj ) 2hij (1 tan 2 )

2hij tan 2 2(hii h jj ) tan 2hij 0

La solucin de esta ecuacin cuadrtica es:

(hii h jj ) hii h jj 4h 2 ij
2

tan

2hij

O bien:

tan

tan

tan

( h

ii

h jj ) (hii h jj ) 2 4hij2

(h

ii

h jj ) (hii h jj ) 2 4hij2

2hij (hii h jj ) (hii h jj ) 2 4hij2

(hii h jj ) 2 (hii h jj ) 2 4hij2

2hij (hii h jj ) (hii h jj ) 2 4hij2

2hij
hii h jj ( hii h jj ) 2 4hij2

(10)

S hii h jj se usar la expresin (10). Y si hii h jj entonces se emplear:

tan

2hij
hii h jj

(hii h jj ) 2 4hij2

(11)

Si arbitrariamente se elige el signo positivo en el denominador de las

ecuaciones ( ) y ( ), entonces:

tan

2hij

(12)

hii h jj (hii h jj ) 2 4hij2

El signo positivo si hii h jj . Signo negativo si hii h jj

Se consideran tambin:
cos (1 tan 2 ) 1 / 2

/ 2 0 / 2 .

(13)

sen cos tan

(14)

Las expresiones (12), (13) y (14) se usan, en lugar de la (6). Para la

elaboracin del programa de computadora en la rotacin de ngulos.
Matriz modal o vectores caractersticos.- una vez obtenidos los
valores caractersticos por el mtodo de Jacob, ahora veremos la forma de
obtener los valores caractersticos. Los vectores modales son matrices
columnas v de la matriz dada H. Los cuales estn asociados con los
valores caractersticos 2 . Existirn n valores caractersticos :
12 . 22 ..... n2 . Los correspondientes n vectores modales v , pueden ser
agrupados en una matriz cuadrada V o sea la matriz modal.
Sea una matriz diagonal, representando los valores caractersticos
conocidos. Entonces tenemos:

H V 2 V

Remultiplicando por la inversa de V o sea V 1

V 1 H V V 1 V 2

Como V 1 V I
V 1 H V 2

Comparando esta expresin con las obtenidas anteriormente o sea:

TT H

T
T

v
T

B Tm ....T3 T2

v
T1

T1 T2 T3 ....Tm

Vemos que V es una matriz cuadrada que representa los valores

caractersticos y es igual a la multiplicacin sucesiva de matrices T
empleadas en la obtencin de valores caractersticos. Entonces la matriz
modal ser obtenida como:

V T1

T2 .....

Tm

Ejemplo.- Las matrices modales para el primero y segundo ejemplo

anteriores sern:
Para el primer problema m = 1

V T1
cos
V
sen

sen
cos

Como 22.5
0.9238 0.3826
V

0.3826 0.9238
0.9238
0.3826
v1
v2

0.3826
0.9238

Para el segundo ejemplo m = 5

V T1 T2 T3 T4

T5

Tomando los valores de T anteriormente obtenidos:

 cos(45)
V sen(45)

sen(45)
cos( 45)
0

0.707
V 0.707
0.000

0.707 0.000
0.707 0.000
0.000 1.000

0.500
V 0.707
0.500

0.707
0.000
0.707

0.500
v1 0.707
0.500

0
cos(6.9)

0 ....... sen(6.9)

1
0
0.8880
0.0000

0.4597

0.000
1.000
0.000

sen(6.9)
cos( 6.9)
0

0
0

0.4597
0.992

0.0000 ...... 0.120

0.000
0.8880

0.120
0.992
0.000

0
0
1

0.500
0.707
0.500
0.707
v 2 0.000
0.707

0.500
v3 0.707
0.500

Hasta aqu se han presentado un conjunto de mtodos numricos para

calcular modos y frecuencias modales, tal vez dando la impresin que son
diferentes. En realidad solamente son diferentes el mtodo de Stodola y el
mtodo de las matrices de transicin. Un tratamiento unificado de estos
mtodos numricos se presenta en el captulo 4 del libro Fundamentals of
Earthquake Engineering de N.M. Newmark y E. Rosenblueth.

Anlisis Modal
Antes de iniciar las aplicaciones del mtodo modal, se
transcribirn algunos artculos de la Normas Tcnicas Complementarias
para diseo por sismo (2004).
Zonificacin
Para los efectos de estas Normas se considerarn las zonas del Distrito
Federal que fija el artculo 170 del Reglamento. Adicionalmente, la zona III
se dividir en cuatro subzonas (IIIa, IIIb, IIIC y IIId), segn se indica en la
figura 1.1.

Coeficiente ssmico
El coeficiente ssmico, c, es el cociente de la fuerza cortante horizontal
que debe considerarse que acta en la base de la edificacin por efecto
del sismo, Vo, entre el peso de la edificacin sobre dicho nivel, Wo.
Con este fin se tomar como base de la estructura el nivel a partir del cual
sus desplazamientos con respecto al terreno circundante comienzan a ser
significativos. Para calcular el peso total se tendrn en cuenta las cargas
muertas y vivas que correspondan, segn las Normas Tcnicas
Complementaras sobre Criterios y Acciones para el Diseo Estructural de
las Edificaciones.
El coeficiente ssmico para las edificaciones clasificadas como del grupo B
en el artculo 139 del Reglamento se tomar igual a 0.16 en la zona I, 0.32
en la II, 0.40 en las zonas IIIa y IIIc , 0.45 en la IIIb y 0.30 en la IIId, ver tabla
3.1, a menos que se emplee el mtodo simplificado de anlisis, en cuyo
caso se aplicarn los coeficientes que fija el Captulo 7 de dicha norma
tabla 7.1. Para las estructuras del grupo A se incrementar el coeficiente
ssmico en 50 por ciento.
ESPECTROS PARA DISEO SSMICO
~}
-J
Cuando se aplique el anlisis dinmico modal que especifica el Captulo 9,
se adoptar como ordenada del espectro de aceleraciones para diseo
ssmico, a, expresada como fraccin de la aceleracin de la gravedad, la
que se estipula a continuacin:
a a 0 (c a 0 )

a c;

T
;
Ta

si T < Ta
si Ta T Tb

a qc;

si T > Tb

Donde :
a

Tb

Los parmetros que intervienen en estas expresiones se obtienen de la

tabla 1.1.

Valores de los parmetros para calcular

los espectros de aceleraciones

ZO
NA

I
II
IIIa
IIIb
IIIc
IIId

c
0.16
0.32
0.40
0.45
0.40
0.30

ao

Ta1

Tb1

0.04
0.2
1.35
0.08
0.2
1.35
0.10
0.53
1.8
0.11
0.85
3.0
0.10
1.25
4.2
0.10
0.85
4.2
1
Periodos en segundos
Tabla 1.1

r
1.0
1.33
2.0
2.0
2.0
2.0

REDUCCIN DE FUERZAS SSMICAS

Factor de reduccin
Los coeficientes que se especifican para la aplicacin del mtodo
simplificado de anlisis toman en cuenta todas las reducciones que
procedan por los conceptos mencionados; por ello, las fuerzas ssmicas
calculadas por este mtodo no deben sufrir reducciones adicionales.
Para el clculo de las fuerzas ssmicas para anlisis esttico y de las
obtenidas del anlisis dinmico modal con los mtodos que se fijan en el
Captulo 9, se emplear un Factor de reduccin Q' que se calcular como
sigue:
Q ' Q;

T
Ta

Q ' 1

si se desconoce T, o si TTa

Q 1 ;

si T< Ta

T se tomar igual al periodo fundamental de vibracin de la estructura

cuando se utilice el mtodo esttico, e igual al periodo natural de
vibracin del modo que se considere cuando se utilice el anlisis dinmico
modal; Ta es un periodo caracterstico del espectro de diseo que se
define en el Captulo 3. Q es el factor de comportamiento ssmico que se
define en el Captulo 5.
Para el diseo de estructuras que sean irregulares, de acuerdo con el
Captulo 6, el valor de Q' se corregir como se indica en dicho Captulo.

GRAFICAS DE LOS ESPECTROS DE ACELERACIONES PARA EDIFICIOS TIPO

"B"

ZONA I

Fig. 2

(zona I edificacin tipo B)

Fig. 3

(zona II edificacin tipo B)

ZONA II

ZONA IIIa

Fig. 4

(zona IIIa edificacin tipo B)

ZONA IIIb

Fig. 5

ZONA IIIc

(zona IIIb edificacin tipo B)

Fig. 6

(zona IIIc edificacin tipo B)

Fig. 7

(zona IIId edificacin tipo B)

ZONA IIId

Ejemplo.

Se har el anlisis dinmico modal de un edificio de cuatro niveles

ubicado en la zona IIIa del Distrito Federal. Esta construccin pertenece al
grupo B con estructuracin a base marcos rgidos de acero. La estructura
es de base cuadrada de 18.00 metros por lado y 13.00 metros de
altura, las otras dimensiones se anotan en la figura 8.
No se harn descripciones del anlisis de cargas, bajada de cargas, anlisis de marcos, etc., as como tampoco se tratar el dimensionamiento
preliminar de secciones de columnas, trabes y determinacin de rigideces.
En la figura 9 se anotan: el peso total de cada entrepiso en toneladas,
adems de rigideces de entrepiso, en toneladas sobre centmetro. La
determinacin de las rigideces de los entrepisos es muy importante, y en
forma aproximada se pueden tomar la suma de rigideces cortantes del
entrepiso, considerando los pisos rgidos, o bien el uso de frmulas como
las de Wilbur. Una mejor aproximacin se tiene cuando se usan mtodos
matriciales o programas elaborados de computadora.
Tendremos dos grupos de marcos, en el sentido x (A, B, C y D) y en el
sentidoy (1, 2, 3 y 4), como se ilustra en la siguiente figura.

PLANTA

Fig. 8
De tal manera que se obtienen 8 marcos, para obtener las rigideces de los
entrepisos se usaron las formulas propuestas por John B. Wilbur, en donde
como se sabe, intervienen las rigideces de las trabes.

Tabulando las rigideces de piso en la direccin "x":

MARCOS
1

A
B
C
D

70
80
80
70

K TX

300

NIVELES
2
60
70
70
60
26
0

3
50
60
60
50
22
0

4
50
60
60
50
220

Tabla 1.2.
En forma anloga se obtiene la otra tabla para la direccin "y" en este
ejemplo, los valores de K TY son los mismos que K TX .
Entonces tendremos 2 modelos dinmicos, uno en el sentido "x" y otro
en el sentido "y". La suma de las rigideces de un entrepiso de los marcos
en un sentido nos representara, la rigidez del resorte del modelo
matemtico como se indica en la figura 21.
6,00

6,00

6,00

K C = 60 Ton/cm

K D = 50 Ton/cm

3,00

K B = 60 Ton/cm

K C = 60 Ton/cm

K D = 50 Ton/cm

3,00

K A = 50 Ton/cm

K C = 70 Ton/cm

K D = 70 Ton/cm

3,00

W = 300 Ton

K C = 80 Ton/cm

K D = 50 Ton/cm

W = 400 Ton
K A = 50 Ton/cm

K B = 60 Ton/cm
W = 400 Ton

K A = 70 Ton/cm

K B = 70 Ton/cm

K A = 50 Ton/cm

K B = 80 Ton/cm

Y
X

Edificio con 4 grados de libertad.

Fig. 9
& estas frmulas para calcular las rigideces de entrepiso, son aplicables para marcos
regulares formados por elementos con momentos de inercia constante. Existen
programas de computadora donde se trabaja con la matriz de rigidez de la estructura.
En la ponencia Diagonalizacin de matrices de rigidez de edificios del ing. Julio Damy
R. (Ver revista Ingeniera U.N.A.M. enero-marzo 1978), presenta la justificacin terica de
trabajar con las rigideces de entrepiso.

4,00

W = 400 Ton

Hiptesis para el anlisis dinmico elstico del edificio.

A) Correspondientes a la estructura.
1. La
estructura
se
comporta
elsticamente;
adems
los
desplazamientos y rotaciones son pequeos.
2. se desprecia el efecto de amortiguamiento en el planteamiento de
las ecuaciones de movimiento, pero se considera incluido en los
espectros de diseo.
3. la fuerza cortante en una direccin, que se produce entre una
columna de dos niveles durante el movimiento en un modo natural,
es proporcional al desplazamiento relativo entre su base y su
extremo superior en dicha accin.
4. Toda la masa se considera concentrada al nivel de las losas.
5. El momento de torsin entre dos niveles durante el movimiento en
un modo natural de vibracin, es proporcional a la rotacin relativa
entre las losas de dichos niveles.
6. Las losas se consideran rgidas (sistema estrechamente acoplado).
B) Correspondiente al movimiento ssmico.
1. El sismo es tal que la base del edificio tiene un movimiento de
translacin en un plano horizontal.
2. El movimiento de la base del edificio se considera conocido en el
sentido en que los espectros de aceleracin son conocidos.
El modelo matemtico del edificio es:
m =0.408

m =0.408

m =0.408

m =0.306
X

K = 300

K = 260

K = 200

K = 220

Ton/cm.

Fig. 10
K

Ki

W
g

Ton seg.2

cm

Para la obtencin del modo fundamental, utilizamos el mtodo de N. M.

Newmark. (Tabla I).
Ntese que despus de obtener los valores de x f en el primer ciclo;
(0.01224, 0.02419, 0.0363 y 0.042), para iniciar el segundo, iniciamos con
los valores de
los cocientes dividiendo los x f entre el primer valor. De esta manera se
har para los dems ciclos.

Tabla 1

Rigidez

Masa
x1

300

260
0.408

0.0

0.01224

xf

x1 x f

F1

0.408

3.508

0.01169

xf

x3

2 ton seg 2

0.0363

0.835

1.215

0.03390

2.02
0.824

0.03867
88.50
3.30

1.184
2.193

0.0114

0.0116

0.00996
0.0230

0.0329

87.58

88.00

88.10

2.015

1.009
1.009

0.0114

1.000

1.05

2.90

3.017

x4

seg.

3.43

87.90

3.425

x3 x f

0.00459

xf

95.00

2 ton seg 2

1.05

0.00996

86.00

0.408

0.042

82.60

2.265

1.00

F1

0.0057

0.02361

85.20
0

0.01151

2.97

0.0116

0.001169

x2 x f

2 ton seg 2

2.04

3.10

2 ton seg 2

1.224

80.50

1.00

1.224

cm.

2.448

0.02479

81.50

4.0

1.224

0.011255

0.01224

3.0

3.264

cm.

0.408 ton seg 2 cm

0.408

0.816

3.672

ton.

220

2.0

0.408

x2

0.408

1.0

F1

220

0.00459
0.0375
88.00

2.885

3.283

F1

0.408

3.412

0.01137

xf

0.822
3.004

0.01137

x4 x f

87.90

1.177
2.182

0.01155

1.005
1.005

0.00992

0.00457

0.02292

0.03284

0.03747

87.85

87.80

87.78

Como en cada punto se debe verificar que x0 x , esto es:

m0 x 0 2
K
1.00 x0 0.01137
2 87.90 0
K0
m0
2.015 x 0 0.02292

m0 x 0 2
K
2 87.85 0
K0
m0

m0 x 0 2
K
2.885 x0 0.03284
2 87.80 0
K0
m0
3.283 x0 0.03747

m0 x 0 2
K
2 87.78 0
K0
m0

El clculo se da por terminado cuando el valor de los cocientes son

prcticamente iguales, como en nuestro caso para el cuarto ciclo.

Se puede obtener una mejor aproximacin de la frecuencia, calculando la

media.

x
x

2 87.83

K0
1.00 2.015 2.885 3.28
0.01137 0.02292 0.03284 0.03747 m0

K0
m0

De acuerdo como hemos calculado, con un factor de escala de 1.00,

K 0 1 m0 1 , por lo tanto :

87.83 9.37

rad .
seg .

El periodo fundamental de la estructura ser:

T

2 6.28

0.671seg .

9.37

El modo fundamental es el vector.

1.000
2.015

V1
2.885

3.283

Para el clculo de los modos superiores usaremos el Mtodo de Holzer.

Antes de aplicarlo a estos, volveremos a obtener el modo fundamental
solo como una comprobacin y observar la diferencia.

Rigidez

300

260

Masa

m 2

ton.

220

0.408

0.408

0.408

0.408

35.8346

35.8346

35.8346

26.876

Fi

35.8346

F1

300

1.00

220

1.00

72.2426

103.5047
191.9228

88.4181

1.01602

0.8724

0.4019

2.8884

ton seg 2 cm

2 87.83

88.4301

264.1654

2.0160

cm.

D 0.01

3.2903

Para obtener el segundo modo nos suponemos una frecuencia 2 600

desde luego debe ser mayor a la primera, ver tabla 2. El tercero y cuarto
modo se anotan en las tablas 3 y 4 respectivamente.

Tabla 2

Rigidez

300

Masa

260

220

0.408

m 2

0.408

245.00

Fi

F1

300

1.00

-298

0.212
1.00

Fi

F1

1.00

268.87

Fi

-268.87

Fi

268.87

-268.464

-28.67

-301.028
-269.492

1.00

0.1213

-1.2249

1.1219

201.348

27.813

31.536

1.00

-1.2027

268.464

300

D 1.06

-241.15
-1.0961

F1

244.56

-0.1058

268.464

2 659

201.654

-1.2264
1.12

268.464

-1.09

-269.82

1.00

m 2

0.150

-300.95

0.1197

D 75

-1.24

268.87

1.00

198
-273

1.227

2 590

181

-36

-1.077

31.13

-1.11

241

0.227

300

2 600

D 65.99

0.112

-237

F1

-1.27

59

1.00

m 2

-1.10

-296

300

-270

241

-241

185.00
212

-243

1.212

241

ton seg 2 cm

-27

55

m 2

245.00

cm.

0.306

0.408

245.00

-245

ton.

220

2 658

241.69

D 0.08

-241.679
-1.098

-0.1022

-1.2017

Tabla 3

Rigidez

300

260

220

220

ton.

cm.

Masa

0.408

m 2

606.288

Fi

300

108.9019

606.696

-606.696

F1

300

-306.696

1.00

-1.795

m 2
Fi

F1

300

1.00

607.512

Fi

-607.512

1.00
0

653.271

-196.520

457.102

-1.1827

-0.8932

2.0777

607.635

607.635

607.635

Fi

-607.635

1.00

-456.46

-307.512

m 2

1.00

655.634 2 1487

653.622

-1.0759

1.0022

607.512

-0.1827

300

653.511

D 0.64

1.0018

455.726

-196.316

457.195

-1.1832

-0.8923

2.0781
-1.0753

2 1490.43

-456.911

-307.635

-0.1835

D 0.513

457.39
2.079

1.00

F1

2 1490

-458.402

-1.0746

110.992

111.319

0.9979
455.94

-0.8903

607.512

300

2.0764

607.92

-0.1843

D 2.735

456.812

-195.88

-1.1843

2 1487

-454.067

-1.0768

-307.92

m 2

F1

654.6061

112.0396

1.00

455.022

-197.794

607.92

- 607.92

0.9916

606.696

-0.1795

607.92

D 5.29

2.0735

-0.8990

1.00

-451.942

-1.0796

108.9091

2 1486

456.1797

-0.178

606.696

654.5484

-0.9016

1.00

Fi

454.716

-197.794

-1.178

m 2

606.288

-306.696

1.00

0.408 ton seg 2 cm

0.408

606.288

-606.696

F1

0.408

D 0.001

1.0033

Tabla 4

Rigidez

300

260

Masa

0.408

m 2

925

Fi

- 925

F1

300

1.00

862.10

Fi

- 862.10

F1

1.00

m 2
Fi

300

x
m 2

-1.2358
-0.3051

862.10

862.10

646.578

-720.80
439.57

-1.1588

861.41

861.41

D 4.68

0.8361

-0.4422

861.288

645.966 2

-711.8545

1.9853

1.00

+285.91

-1.2783

437.136

-2.1388

2 2113

-281.23

1.9980

-561.41

1.00

D 4.39

-271.88

0.7637

998.06

2 2110

+267.49

-1.1384

861.288

-861.288

F1

645.66

-707.2129

-1.1619

861.288

-1.83

860.88

1.9787

-2.1619

1.67

435.33

-562.10

1.00

D 500

-3.5

+1001.67

300

+1270
-770

+996.21

1.00

m 2

-1545

860.88

-2.1572

2 2260

692

3.07

-560.88

1.00

925

-1.40

300

925

cm.

0.408 ton seg 2 cm

+675

-860.88

F1

0.408

-2.40

860.88

Fi

0.408

-625

ton.

220

+1300

1.00

m 2

220

2111

273.760

D 1.2

-276.111
-1.2503

0.8265

861.41

-0.4238

646.06

2 2111 .30

Fi

-861.41

F1

300

-561.41

1.00
0

998.546

-713.247
437.136

-276.111

1.9872

-1.255

-2.1592
1.00

+275.867

-1.1592

0.8280

D 0.25

-0.4270

De esta manera se ha obtenido la matriz modal.

1.00
2.016
v
2.888

3.290

1.00

1.00

1.1219
0.1022
1.2017

0.1835
1.0753
1.0033

1.00
1.159
0.828

0.427

Teniendo los siguientes periodos y frecuencias

1 9.37 Rad Seg

T1 0.671Seg.

2 25.65 Rad Seg

T2 0.245Seg.

3 38.59 Rad Seg

T3 0.163Seg .

4 45.95 Rad Seg

T4 0.137 Seg.

Podemos presentarlos de esta manera.

1.

2.

Modo fundamental y modos superiores

Fig. 11
Participacin de los modos de vibracin

3.

4.

Hasta ahora solamente se han expuesto mtodos para obtener los modos
de vibracin, pero falta tener una expresin para calcular los
desplazamientos reales.
Recordando la ecuacin matricial (33), ella expresa el vector de
desplazamientos de la estructura en trminos de las contribuciones
modales.
xi (T )

vT m

vr vrT rmvr

x 0 cos r T

r 1

1
x '( 0 ) sen r T
r

Para ser explicito, consideremos un sistema con i = 3, la ecuacin anterior

toma la forma:

m1
v1r v2r v3r 0
0

x1
v1r
x v
2
2r
m1
x3
v3r
v1r v2 r v3r 0
0

0
m2
0
0
m2
0

0
0
m3
0
0
m3

v1r
v
2r
v3r

r T ......................1

En donde:
Vector de desplazamiento relativo en funcin de los modos de
vibracin para un
instante cualquiera T.
xi

v iR

Modo de vibracin R a escala arbitraria.

r T x 0 cos r T

x' 0
sen r T
r

Efectuando algunas operaciones matriciales en (1), Tenemos:

x1
v
n
v1r m1 v2 r m2 v3r m3 T v1r
x
r
2
2
2
2
2r
r 1 v1r m1 v 2 r m2 v3 r m3
x3
v3r
n

O bin:

xi
r 1

m v
i 1
n

m v
i 1

xi Crr T vir
r 1

i ir
2
i ir

r T vir

Llamando:

Cr

n
I 1
n

i 1

mi vir

m1v12r

Considerando que el terreno sufre sbitamente una aceleracin igual a

ag siendo ar la ordenada del espectro de diseo.
r Es una funcin que depende de las condiciones iniciales.
Por lo que para T 0 , al actuar el sismo:
x ' ( 0)
sen r T
r
x ' ( 0) 2
r x (0) r2 cos r T
r sen r T
r

r x (0) cos r T

r x ( 0)

r x (0) r2
Sin considerar el signo negativo obtenemos:

r
r2

Pero r podemos considerarla como la aceleracin del terreno, o sea:

r ar g

Para el caso de los espectros propuestos para las distintas zonas en que se
divide el distrito federal, ya se tiene involucrado los efectos de
amortiguamiento, a r esta expresado como un porcentaje de la gravedad.
ximax . r 1 a r g
n

Cr
vir
r2

Obsrvese que xi es el vector de desplazamiento total. Resultado de la

participacin de todos los vectores de desplazamientos reales. Esto es:
ximax . r 1U ir
n

max:

Finalmente para el modo r.

Cr
vir
max:
r2
De acuerdo al Capitulo 4 de las NTC del reglamento de construcciones del
distrito federal, se deber dividir la ordenada del espectro de diseo entre
el factor Q, sin embargo se puede aplicar posteriormente.
U ir

ar g

A continuacin se empleara esta ultima expresin U ir , pero antes

obtendremos las ordenadas del espectro de aceleraciones de diseo para
la zona IIIa. Esto es:
ZO
NA

ao

Ta1

Tb1

IIIa

0.40

0.10

0.53

1.8

2.0

a a 0 (c a 0 )

a c;
a qc;

Primer modo

T
;
Ta

si T < Ta
si Ta T Tb
si T > Tb

T1 0.671seg entonces utilizamos: a c;

Sustituyendo valores, tenemos:

a 0.40

Este valor se puede obtener de la grafica de espectros de aceleraciones,

pagina 4.
Segundo modo

T2 0.245seg .

Tercer modo

T3 0.163seg .

a 0.192

Cuarto modo

T4 0.137 seg.

a 0.178

a 0.239

Calculo de los coeficientes de participacin C r .

c1

(0.408 x1.00) (0.408 x 2.016) (0.408 x 2.888) (0.306 x3.29)

0.3889
(0.408 x1.00 2 ) (0.408 x 2.016 2 ) (0.408 x 2.888 2 ) (0.306 x3.29 2 )

(0.408 x1.00) (0.408 x1.1219) (0.408 x 0.1022) (0.306 x 1.2017)

0.3321
(0.408 x1.00 2 ) (0.408 x1.1219 2 ) (0.408 x 0.1022 2 ) (0.306 x 1.2017 2 )
(0.408 x1.00) (0.408 x 0.1835) (0.408 x 1.0753) (0.306 x 1.0033)
c3
0.1674
(0.408 x1.00 2 ) (0.408 x 0.1835 2 ) (0.408 x 1.0753 2 ) (0.306 x 1.0033 2 )
(0.408 x1.00) (0.408 x 1.159) (0.408 x0.828) (0.306 x 0.427)
c4
0.1101
(0.408 x1.00 2 ) (0.408 x 1.159 2 ) (0.408 x0.828 2 ) (0.306 x 0.427 2 )
c2

Dividiendo cada coeficiente entre

su frecuencia al cuadrado
correspondiente y llamndole I r , obtenemos los siguientes valores:

I1 0.004428
I 2 0.0005047
I 3 0.0001124

I 4 0.00005214

Ahora ya podemos substituir simplemente las cantidades anteriores:

U 11 (0.337)(981)(0.004428)(1.00) 1.46201cm.

Veamos la influencia que tienen los modos en su participacin, indicndolo

en porcentajes:

I1 I 2 I 3 I 4 10 5 (442.7 50.71 11.22 5.33) 509.96 x10 5

Para el primer modo:

i1

442.70
x100 86.85%
509.96

Para los modos superiores:

i2 9.9%

i3 2.2%

i4 1.05%

Como puede observarse la influencia mayor es la del primer modo, por lo

cual debern tomar en cuenta los primeros modos, dependiendo del
nmero de niveles del edificio.
Calculo de las fuerzas cortantes ssmicas.
Las fuerzas cortantes ssmicas en cada entrepiso, para cualquier modo,
pueden obtenerse como el producto de la rigidez del entrepiso por el
desplazamiento relativo de dicho entrepiso.

Para evaluar las fuerzas ssmicas, es necesario considerar la reduccin

por ductilidad, que al respecto se ha establecido.
Q ' Q;

T
Ta

Q ' 1

Si se desconoce T, o si T Ta

Q 1 ;

Si T< Ta

Para este problema tomamos Q=3, que depende del tipo de estructura, de
los materiales a usar y de algunos detalles de diseo y construccin.
Como el periodo fundamental es menor que T1, tenemos:
T1 0.671seg Ta
Q' 3

Este factor, es el divisor de las ordenadas de los espectros de

aceleraciones de diseo, lo consideramos hasta este punto porque en el
reglamento se establece; que los desplazamientos no deben dividirse
entre dicho factor. Desde luego, estas disposiciones se basan en
investigaciones.
Las respuestas modales, se combinan de acuerdo a la expresin dada:
Vi

1
2
ri
Q'

Sustituyendo valores:
1
438.604 2 27.824 2 5.20 2 2.328 2 174.15ton.
3
V2 153.02ton.
V3 111 .70ton.

V1

V4 52.26ton.

Ya se ha mencionado la gran participacin del modo fundamental, no

obstante, en el reglamento se considera que debern tomarse cuando
menos tres modos de vibracin.
Se considera conveniente hacer las siguientes observaciones: la expresin
anterior para combinar efectos modales es vlida para casos en donde los
periodos de vibracin estn claramente diferenciados unos de otros. En
caso contrario deber aplicarse otra expresin, presentada por el Dr. E.
Rosenblueth en el trabajo sobre la respuesta ssmica de estructuras de
comportamiento lineal (1968).
Tambin se hace la siguiente aclaracin respecto a la superposicin modal,
en donde se pide que las respuestas modales (fuerzas cortantes,
deformaciones, momentos de volteo, etc.) se combinen de acuerdo a la
expresin ya mencionada; si se superponen directamente los
desplazamientos relativos y a partir de los nuevos desplazamientos se
obtienen las fuerzas cortantes ssmicas, estos valores difieren de los
valores obtenidos como se hizo en este ejemplo por lo que se recomienda
se haga en esta forma. El problema anterior se debe a que el criterio de
combinacin no es una transformacin lineal entre las contribuciones
modales.
Para edificios muy esbeltos con determinadas caractersticas, pueden
tener en los niveles superiores considerable influencia alguno de los
modos superiores.

Ejemplo.
Se har el anlisis ssmico modal de un tanque elevado sobre una
estructura reticular articulada, ubicada en la zona 11 del D.F. Las
dimensiones de la estructura y del tanque se dan en la figura No. 24.
Esta construccin pertenece al grupo B con estructuracin tipo 1. Se
omiten los clculos de cargas, as como el dimensionamiento de la
armadura que soporta al tanque.

Figura
Pesos y Masas:
Peso de la estructura

4500 Kgs. Masa =

4500
= 0.4587
981

ton seg 2
m

Peso del tanque

Peso del Lquido

6208 kgs. Masa = 0.6328 = 0.4587

6195 kgs. Masa = 0.6315 = 0.4587

Ejemplo.Se har el anlisis ssmico modal de un tanque elevado sobre una

estructura reticular articulada, ubicada en la zona II del D.F., las
dimensiones de la estructura y del tanque se dan en la figura No. 24.
Esta construccin pertenece al grupo B con estructuracin tipo I. Se
omiten los clculos de cargas, as como el dimensionamiento de la
armadura que soporta al tanque.

Planta de la estructura
metlica de soporte

Figura No. 24

Tanque para almacenamiento

de agua.
Base

2.00x2.00 m.
Acotacin en metros.
Peso y masas:
Peso de la estructura

4500 kgs.

Peso del tanque

6208 kgs.

Peso del lquido

6195 kgs.

4500
ton seg 2
0.4587
981
m
ton seg 2
Masa 0.6328
m
ton seg 2
Masa 0.6315
m
Masa

Al idealizar la estructura y el fluido en funcin de las masas equivalentes,

se necesitan tambin las rigideces:

Calculo de la rigidez de la estructura de soporte.

Se calculara la rigidez de una armadura plana, forma aproximada,
emplendose el mtodo de la carga virtual unitaria para obtener el
desplazamiento en la parte superior de la armadura debido a una fuerza P.
d i 1 N i ni
n

Li
Ei Ai

Donde:
d : Desplazamiento en el extremo libre de la armadura.
N i : Fuerzas axiales en las barras debido a la carga P.
ni : Fuerzas axiales en las barras debido a la carga Virtual Unitaria.
Li : Longitudes de las barras
Ai : reas de las secciones transversales de las barras.
E i : Modulo de elasticidad de los materiales.
En la tabla I se anotan las reas de los perfiles laminados y las longitudes
de las barras, respetando la numeracin asignada en la figura No. 24.
Calculo de las reacciones verticales en los apoyos de la estructura debido
a la fuerza P.
M A 17 P 6 RVB 0

RVB RVA 2.8333

Para analizar esta armadura y obtener las fuerzas axiales N i se pueden

usar mtodos grficos, analticos manuales o emplear algn programa de
computadora. Los valores de las fuerzas axiales, para esta estructura
isostatica, estn anotadas en la cuarta columna de la tabla I, y en la
quinta columna se encuentran los valores debidos a la carga unitaria
(P=1).
En la figura No. 25 se ilustran las fuerzas axiales de las barras (tensin o
compresin) sobre el esquema de la estructura. Con lneas punteadas se
muestran las fuerzas equivalentes empleadas para analizar la estructura.
Se sugiere tomar algunos nodos (i, j, g, etc.) y comprobar el equilibrio.
En la tabla I se efectan las operaciones y al sumar los valores de la
ltima columna se obtiene el valor del desplazamiento en funcin de P.
d i 1 N i ni
16

Li
0.41021635 P
E i Ai

Por definicin:
P k (rigidez )

Si d 1

k1 2.8621ton

cm (Plana)

Rigidez de la armadura:

k 5.7242 ton

cm

Valor correspondiente al primer resorte del modelo.

Calculo de masas equivalentes.

En escritos de diversos investigadores sobre el comportamiento dinmico
de lquidos en vasos de almacenamientos, se han presentado resultados
de diferentes estudios. Para nuestro problema, en el capitulo 6 del libro de
N. M. Newmark y E. Rosenblueth (ver bibliografa), en la seccin vibracin
de los lquidos en tanques se dan expresiones para calcular las masas
equivalentes profundidades y rigidez del liquido.

Se ha encontrado que en las presiones de los tanques tiene importancia el

oleaje y se puede despreciar la compresibilidad del agua. Para un tanque
rgido completamente lleno que esta cubierto con una tapa rgida, la masa
total del liquido se mueve con el tanque como si se tratase de una masa
rgida no obstante, cuando se tiene un espacio pequeo entre la superficie
del liquido en las paredes y el fondo sern prcticamente iguales a las que
se produciran con superficie libre. Se recomienda estudiar las condiciones
del tanque, totalmente lleno y con superficie libre.
Se demuestra que el liquido puede sustituirse por masas unidas al tanque
mediante resortes lineales y elementos rgidos, para calcular las fuerzas
horizontales que representan la accin resultante de los empujes
hidrodinmicos.
El reglamento de construcciones, en el art. 245 dice
Para tanques rectangulares que miden
2 L en la direccin del
movimiento, los parmetros que definen a la estructura equivalente se
obtienen con las siguientes expresiones:

M0

M1

1.7 L
H M
1.7 L
H

Tanh

1.6 H
L M
1.6 H
L

0.83Tanh

M

H 0 0.83H 1
1
M0

M L
H 1 H 1 0.833

M1 H

0.63

L
H

ML
0.28
M1H

3 gM 12 H
k
ML2
En donde:

M : Masa del lquido contenida en el tanque.

M 0 : Masa del lquido considerada unida rgidamente al tanque.
M 1 : Masa del lquido unida por medio de resortes.
H : Profundidad de reposo del lquido contenido.
H 0 : Elevacin de la masa desde el fondo.
H 1 : Elevacin de la masa desde el fondo.

k:

Rigidez de los resortes para la masa.

L : Dimensin de la base del tanque en la direccin del movimiento.
g : Constante de gravedad.

1.33
2.00

Si se considera en los clculos momentos hidrodinmicos del fondo del

tanque.
0.00
1.00

Si solamente importan los efectos de presiones hidrodinmicas en las

paredes del recipiente.
Substituyendo valores:
ton seg
Para 0.00 ; 1.00 ; H 2.00m. ; L 1.00m. y M 0.6315

M0

M1

1. 7 L
2
2.0 M 0.813M 0.5134 ton seg
1. 7 L
m
2.0

Tanh

1.6 x 2.0
2
1.0 M 0.2585M 0.1632 ton seg
1.6 x 2.0
m
1.0

0.83Tanh

H 0 0.38 x 2.0 0.76m.

M
1 .0
H 1 2.0 1 0.33

0.2585M 2.0

3 x9.81(0.2585) 2 M 2 x 2
M (1.0) 2

0.63 x1.0

1 .0
2 .0

0.28

0.2585Mx 2

3.93315M 2.4837 ton

Masa considerada de la estructura: M 0.229

1 1.499 m.

m.

ton seg 2
m

Calculo de los modos de Vibracin.

Se proceder a calcular los modos de Vibracin, las frecuencias y los
periodos. Usando el procedimiento que se expuso en el capitulo V.
Idealizacin de la estructura para el anlisis dinmico:

m0'

m1

K
Matriz de rigideces:
5.74905
k
0.024837
0.013752
0.0000

k 2m

0.024837
0.024837
0.0000
0.001632

5.74905 0.013752 2
0.024837

0.024837
0
0.024837 0.001632 2

Efectuando operaciones.

2 2

433.27048 2 6334.73919 0

Resolviendo la ecuacin cuadrtica, encontramos las frecuencias de

vibracin:

12 15.1505

12 3.892 rad seg

22 418.1199

22 40.447 rad seg

Periodos de Vibracin:
T1

T2

1.614 seg .

2
0.307 seg .
2

A continuacin se empleara esta ultima expresin U ir , pero antes

obtendremos las ordenadas del espectro de aceleraciones de diseo para
la zona IIIa. Esto es:
ZO
NA

ao

Ta1

Tb1

IIIa

0.40

0.10

0.53

1.8

2.0

a a 0 (c a 0 )

a c;
a qc;

Primer modo

T
;
Ta

si T < Ta
si Ta T Tb
si T > Tb

T1 0.671seg entonces utilizamos: a c;

Sustituyendo valores, tenemos:

a 0.40

Este valor se puede obtener de la grafica de espectros de aceleraciones,

pagina 4.
Segundo modo

T2 0.245seg .

Tercer modo

T3 0.163seg .

a 0.192

Cuarto modo

T4 0.137 seg .

a 0.178

a 0.239

Calculo de los coeficientes de participacin C r .

c1

(0.408 x1.00) (0.408 x 2.016) (0.408 x 2.888) (0.306 x3.29)

0.3889
(0.408 x1.00 2 ) (0.408 x 2.016 2 ) (0.408 x 2.888 2 ) (0.306 x3.29 2 )

(0.408 x1.00) (0.408 x1.1219) (0.408 x 0.1022) (0.306 x 1.2017)

0.3321
(0.408 x1.00 2 ) (0.408 x1.1219 2 ) (0.408 x 0.1022 2 ) (0.306 x 1.2017 2 )
(0.408 x1.00) (0.408 x 0.1835) (0.408 x 1.0753) (0.306 x 1.0033)
c3
0.1674
(0.408 x1.00 2 ) (0.408 x 0.1835 2 ) (0.408 x 1.0753 2 ) (0.306 x 1.0033 2 )
(0.408 x1.00) (0.408 x 1.159) (0.408 x0.828) (0.306 x 0.427)
c4
0.1101
(0.408 x1.00 2 ) (0.408 x 1.159 2 ) (0.408 x0.828 2 ) (0.306 x 0.427 2 )
c2

Dividiendo cada coeficiente entre

su frecuencia al cuadrado
correspondiente y llamndole I r , obtenemos los siguientes valores:
I1 0.004428
I 2 0.0005047
I 3 0.0001124

I 4 0.00005214

Ahora ya podemos substituir simplemente las cantidades anteriores:

U 11 (0.337)(981)(0.004428)(1.00) 1.46201cm.

Veamos la influencia que tienen los modos en su participacin, indicndolo

en porcentajes:

I1 I 2 I 3 I 4 10 5 (442.7 50.71 11.22 5.33) 509.96 x10 5

Para el primer modo:

i1

442.70
x100 86.85%
509.96

Para los modos superiores:

i2 9.9%
i3 2.2%

i4 1.05%

Como puede observarse la influencia mayor es la del primer modo, por lo

cual debern tomar en cuenta los primeros modos, dependiendo del
nmero de niveles del edificio.

Calculo de las fuerzas cortantes ssmicas.

Las fuerzas cortantes ssmicas en cada entrepiso, para cualquier modo,
pueden obtenerse como el producto de la rigidez del entrepiso por el
desplazamiento relativo de dicho entrepiso.

Para evaluar las fuerzas ssmicas, es necesario considerar la reduccin

por ductilidad, que al respecto se ha establecido.
Q ' Q;

Si se desconoce T, o si T Ta

 T
Ta

Q ' 1

Q 1 ;

Si T< Ta

Para este problema tomamos Q=3, que depende del tipo de estructura, de
los materiales a usar y de algunos detalles de diseo y construccin.
Como el periodo fundamental es menor que T1, tenemos:
T1 0.671seg Ta
Q' 3

Este factor, es el divisor de las ordenadas de los espectros de

aceleraciones de diseo, lo consideramos hasta este punto porque en el
reglamento se establece; que los desplazamientos no deben dividirse
entre dicho factor. Desde luego, estas disposiciones se basan en
investigaciones.

Las respuestas modales, se combinan de acuerdo a la expresin dada:

Vi

1
2
ri
Q'

Sustituyendo valores:
1
438.604 2 27.824 2 5.20 2 2.328 2 174.15ton.
3
V2 153.02ton.
V3 111 .70ton.

V1

V4 52.26ton.

Ya se ha mencionado la gran participacin del modo fundamental, no

obstante, en el reglamento se considera que debern tomarse cuando
menos tres modos de vibracin.
Se considera conveniente hacer las siguientes observaciones: la expresin
anterior para combinar efectos modales es vlida para casos en donde los
periodos de vibracin estn claramente diferenciados unos de otros. En
caso contrario deber aplicarse otra expresin, presentada por el Dr. E.
Rosenblueth en el trabajo sobre la respuesta ssmica de estructuras de
comportamiento lineal (1968).
Tambin se hace la siguiente aclaracin respecto a la superposicin modal,
en donde se pide que las respuestas modales (fuerzas cortantes,
deformaciones, momentos de volteo, etc.) se combinen de acuerdo a la
expresin ya mencionada; si se superponen directamente los

desplazamientos relativos y a partir de los nuevos desplazamientos se

obtienen las fuerzas cortantes ssmicas, estos valores difieren de los
valores obtenidos como se hizo en este ejemplo por lo que se recomienda
se haga en esta forma. El problema anterior se debe a que el criterio de
combinacin no es una transformacin lineal entre las contribuciones
modales.
Para edificios muy esbeltos con determinadas caractersticas, pueden
tener en los niveles superiores considerable influencia alguno de los
modos superiores.