Sociedad Mexicana de Ingeniera Estructural

MODELACIN DE LA INTERACCIN DINMICA

SUELO-ESTRUCTURA USANDO ECOgcW
Gerardo A. Corona Carlos 1 y Javier Avils Lopez 2
RESUMEN
Para edificios diseados conforme a reglamento, la interaccin dinmica entre la estructura y el suelo se tiene
en cuenta slo en el modo fundamental de vibracin, aplicando un modelo simplificado que considera los
efectos en el periodo, amortiguamiento y ductilidad. Sin embargo, con la ayuda de computadoras personales
es posible desarrollar modelos tridimensionales del sistema suelo-estructura. En este trabajo se ilustra la
forma de modelar la interaccin en estructuras sobre cimentaciones someras y con pilotes usando el programa
de cmputo ECOgcW. El programa define un nudo maestro en la base de la estructura, donde condensa los
grados de libertad de la cimentacin y concentra los resortes elsticos que representan la rigidez del suelo. Se
calculan las fuerzas ssmicas elsticas usando anlisis modal espectral en combinacin con los nuevos
espectros de sitio reglamentarios. Estas fuerzas se reducen por ductilidad con los nuevos factores de reduccin
de resistencia que consideran efectos de sitio e interaccin suelo-estructura. El procedimiento propuesto
puede emplearse con otros programas de anlisis estructural, siempre que los resortes elsticos se determinen
por separado y la reduccin por ductilidad se haga apropiadamente. Asimismo, se muestra un ejemplo
numrico del anlisis ssmico de un edificio tpico con distintos tipos de cimentacin.
ABSTRACT
For code-designed buildings, soil-structure interaction is accounted for on the fundamental mode of vibration
only, by using a simplified model that considers the effects on the period, damping and ductility. With the aid
of personal computers, however, it is possible to develop three-dimensional models of the soil-structure
system. In this work, the way of modeling the interaction effects in structures on shallow and pile foundations
is illustrated, making use of program ECOgcW. A master node at the base of the structure is defined by the
program, where the degrees of freedom of the foundation are condensated and the elastic springs representing
the soil stiffness are concentrated. The elastic seismic forces are computed by the use of modal dynamic
analysis in combination with the new site-specific design spectra specified in the code. These forces are
reduced by ductility with the new strength reduction factors that consider both site and interaction effects. The
proposed procedure can be used with other programs of structural analysis, providing that the elastic springs
are previously known and the reduction by ductility is properly made. Also, a numerical example of the
seismic analysis of a typical building with different types of foundation is shown.
INTRODUCCIN
La interaccin dinmica suelo-estructura consiste en un conjunto de efectos producidos en la estructura y el
suelo como resultado de la deformabilidad de ste ante excitacin ssmica. La interaccin modifica las
propiedades dinmicas relevantes que tendra la estructura supuesta con base indeformable, as como las
caractersticas del movimiento del suelo en la vecindad de la cimentacin. El alargamiento del periodo
fundamental, la modificacin del amortiguamiento asociado y la reduccin de la ductilidad, respecto a los
valores de base rgida, son producto de la interaccin inercial (Avils y Prez-Rocha, 2003). Si el cimiento es
suficientemente rgido, la interaccin cinemtica se manifiesta en una reduccin de los movimientos de
traslacin de la cimentacin (filtrado de los componentes de alta frecuencia de la excitacin ssmica), as
como en la generacin de cabeceo y torsin de la cimentacin (Luco y Wong, 1987). No obstante lo anterior,
se ha observado (Avils y Prez-Rocha, 2003) que para la ciudad de Mxico es conservador efectuar slo el
1

G. C. Ingeniera y Diseo, S.C. Calzada del Arco 82, Col. Sta. Cruz Guadalupe, 72160 Puebla, Pue.
Telfono: (222) 284.5969, Fax (222) 284.59993, [email protected]

Instituto Mexicano de Tecnologa del Agua, Pase Cuauhnahuac 8532, 62550 Jiutepec, Mor. Telfono:
(777) 329.3600; Fax: (777) 329.3679, [email protected]
1

XIV Congreso Nacional de Ingeniera Estructural

Acapulco, Gro., 2004

anlisis de interaccin inercial, siempre que los efectos de sitio sean considerados en la determinacin del
movimiento de campo libre empleado como movimiento de entrada para la cimentacin.

Figura 1. Espectros de sitio para diseo sin (lnea continua) y con (lnea discontinua) interaccin, para
un sitio con Ts = 2 s y una estructura con H e R = 3 y D R = 1.

Debido a la excepcional deformabilidad de las arcillas de la ciudad de Mxico, los efectos de sitio e
interaccin suelo-estructura que se presentan durante temblores intensos alcanzan niveles muy superiores a
los usuales en otras ciudades asentadas sobre depsitos de suelo blando (Resndiz y Roesset, 1986). Pese a
que ambos fenmenos se relacionan entre s, para propsitos de diseo es conveniente reconocerlos por
separado. En las nuevas disposiciones reglamentarias para el Distrito Federal (NTCDS, 2004) as se hace, al
definir un espectro especfico de sitio que se modifica por interaccin suelo-estructura, como se ilustra en la
figura 1.

Figura 2 Cociente de cortantes con base flexible y rgida, para un sitio con Ts = 2 s y una estructura
con H e R = 3 y D R = 1.

Sociedad Mexicana de Ingeniera Estructural

En la ciudad de Mxico, las amplificaciones dinmicas del subsuelo adquieren valores extraordinariamente
altos cuando los periodos predominantes de la excitacin y el suelo son similares. De igual forma, la
interaccin puede ocasionar considerables incrementos o reducciones de la respuesta estructural, dependiendo
de la relacin entre los periodos fundamentales de la estructura y el sitio. Especficamente, cuando el periodo
de la estructura es menor que el del sitio, la interaccin es perjudicial, pero en caso contrario es benfica,
como se muestra en la figura 2. Es interesante notar que mientras los mayores incrementos en las resistencias
requeridas ocurren en sistemas no lineales, las mayores reducciones suceden en sistemas lineales. Esto es
debido esencialmente al alargamiento del periodo estructural y la disminucin del factor de reduccin de
resistencia al considerar la interaccin. Esto ltimo se aprecia en la figura 3 aplicando el criterio especificado
en las normas para la reduccin por ductilidad en funcin de las condiciones de sitio.

Figura 3. Factores de reduccin de resistencia sin (lnea continua) y con (lnea discontinua)
interaccin, para un sitio con Ts = 2 s y una estructura con H e R = 3 y D R = 1.

Los resultados que se muestran en la figura 2 son los que se obtienen, segn las normas, al tenerse en cuenta
los efectos de interaccin exclusivamente en el modo fundamental de vibracin de la estructura. No obstante,
con la ayuda de computadoras personales es posible desarrollar modelos de interaccin multimodales que
satisfagan los requisitos reglamentarios. El propsito de este trabajo es presentar un mtodo de anlisis
ssmico tridimensional implantado en el programa de cmputo ECOgcW (2002), que cumple con los
requisitos de interaccin de las normas. El mtodo est basado en el anlisis modal espectral, modelando la
rigidez del suelo con los resortes elsticos especificados en las normas. Es particularmente importante
modelar la rigidez del suelo en funcin de la frecuencia de vibracin. Los espectros de sitio y los factores de
reduccin de resistencia utilizados con este mtodo son exactamente como los mostrados en las figuras. 1 y 3,
respectivamente.

XIV Congreso Nacional de Ingeniera Estructural

Acapulco, Gro., 2004

PROCEDIMIENTO SIMPLIFICADO REGLAMENTARIO

De acuerdo con las normas, si se opta por el anlisis esttico con la opcin que toma en cuenta el periodo
fundamental de vibracin, el cortante basal de la estructura con base indeformable est dado por
Vo =

a
Wo
Q'

(1)

donde a es la ordenada espectral elstica para el periodo y amortiguamiento de la estructura con base
indeformable y Q el correspondiente factor de reduccin por ductilidad. Por simplicidad, aqu y en lo que
sigue se omite el factor de reduccin por sobrerresistencia, ya que no depende de las condiciones de sitio ni de
la interaccin suelo-estructura. La contribucin de los modos superiores de vibracin est implcitamente
considerada al emplear el peso total de la estructura Wo en lugar de su peso efectivo We 0.7Wo . Si se
considera la flexibilidad del suelo, el cortante basal modificado por interaccin es
a a~
~
a
Vo = Wo ~ We
Q'
Q' Q '

(2)

donde a~ es la ordenada espectral elstica para el periodo y amortiguamiento del sistema (consultar la norma)
~
y Q el correspondiente factor de reduccin por ductilidad. El segundo trmino de esta expresin representa
la reduccin del cortante basal por efectos de interaccin en el modo fundamental de vibracin. El signo de
este trmino puede ser negativo, con lo que en lugar de una reduccin tendremos un incremento del cortante
~
basal. Una vez calculados Vo y Vo , los efectos de interaccin se consideran multiplicando por el cociente
~
Vo
a~ Q'
= 0.3 + 0.7 ~
Vo
a Q'

(3)

cualquier respuesta estructural calculada para la condicin de base rgida. La variacin de este factor con el
periodo estructural Te es la que se ejemplifica en la figura 2 para ciertas condiciones de interaccin. En
~
ninguna situacin, el valor de Vo Vo ha de tomarse menor que 0.75 ni mayor que 1.25, segn las normas.
Los desplazamientos laterales de la estructura en la direccin de anlisis, corregidos por interaccin con el
suelo, se determinan con la expresin
~
~ V
Xi = o
Vo

M
X i + (hi + D) o
Kr

(4)

donde M o es el momento de volteo en la base de la cimentacin y X i el desplazamiento lateral del i-simo

nivel de la estructura con base rgida, calculados usando las fuerzas ssmicas sin modificar por interaccin; D
es la profundidad de desplante y K r la rigidez de rotacin de la cimentacin. El segundo trmino de esta
expresin representa la contribucin de la rotacin de la cimentacin. En las normas se especifica que el valor
~
de X i relativo al terreno debe tomarse en cuenta en el clculo de los efectos de segundo orden, as como en
la revisin del estado lmite por choques con estructuras adyacentes. Ntese que la determinacin del factor
~
Vo Vo permite evaluar las consecuencias de la interaccin en el cortante basal y el desplazamiento lateral de
la estructura.
El procedimiento de interaccin suelo-estructura descrito anteriormente ha sido implantado en el programa de
cmputo ECOgcW. Esta opcin puede ser activada cuando se recurre al anlisis ssmico esttico. Si se opta
por el anlisis dinmico modal, los efectos de interaccin se tienen en cuenta usando un modelo
tridimensional del sistema suelo-estructura, el cual se explica en la siguiente seccin.

Sociedad Mexicana de Ingeniera Estructural

MODELO COMPUTACIONAL TRIDIMENSIONAL
Para tener en cuenta los efectos de interaccin en los modos superiores de vibracin es necesario analizar el
sistema completo. Para propsito de anlisis dinmico, la mayor hiptesis requerida es suponer que la
cimentacin posee suficiente rigidez y resistencia de suerte que se mueve como cuerpo rgido. De esta forma
es posible definir un nudo maestro en la base de la estructura donde se condensan los grados de libertad de la
cimentacin, a saber: traslacin en dos direcciones horizontales ortogonales, cabeceo alrededor de los
respectivos ejes y torsin alrededor del eje vertical.
La rigidez del suelo se modela con elementos estructurales sin masa, es decir: resortes elsticos de traslacin,
cabeceo y torsin. Estos parmetros expresan tanto la rigidez como la inercia del suelo y son dependientes de
la frecuencia de vibracin; dicha dependencia se debe a la influencia que la frecuencia tiene en la inercia del
suelo, puesto que su rigidez es esencialmente independiente de ella. Para su clculo se aplican los criterios
que se detallan en las normas, aunque pueden emplearse mtodos alternos basados en principios establecidos
y resultados conocidos de la dinmica de cimentaciones.
Siguiendo el enfoque de rigidez usado en anlisis estructural, las ecuaciones de vibracin libre del sistema
tienen la siguiente forma matricial:
K ee

K ce

K ec
M
2 ee
K cc
0

0 U e 0
=
M cc U c 0

(5)

donde M ii y K ii son submatrices de masa y rigidez, respectivamente; los nudos de la estructura se

identifican con el subndice e y los de la cimentacin con el subndice c. De acuerdo con las normas, la
masa de la cimentacin y sus inercias rotacionales son despreciadas. Adems, U e y U c son los vectores de
desplazamiento de la estructura y la cimentacin, respectivamente.
El clculo de las frecuencias y modos naturales de vibracin del sistema se realiza aplicando a la ec. 5
cualquier tcnica estndar de valores caractersticos. Sin embargo, como los resortes de la base son funcin de
la frecuencia predominante de vibracin, el anlisis de vibraciones libres se hace por medio de iteraciones,
iniciando con los valores de los resortes para la condicin de base rgida y terminando con los
correspondientes a la frecuencia del sistema.
Las fuerzas ssmicas modales se determinan haciendo uso de los nuevos espectros de sitio reglamentarios,
considerando el amortiguamiento del sistema para el modo fundamental ( 1) y 5% de amortiguamiento
para los modos superiores de vibracin ( = 1). Estas fuerzas se reducen por ductilidad con los nuevos
factores reglamentarios de reduccin de resistencia, considerando la ductilidad del sistema (consultar la
norma) para el modo fundamental y la ductilidad nominal Q para los modos superiores de vibracin. De esta
forma se cumple con la disposicin reglamentaria de que los efectos de interaccin en el amortiguamiento y la
ductilidad han de tomarse en cuenta slo en el modo fundamental de vibracin; la contribucin de los modos
superiores se evala como se establece para estructuras de base rgida.
Cabe sealar que los desplazamientos laterales de la estructura que se obtienen con este procedimiento ya
incluyen la contribucin de la rotacin de la cimentacin, incluso la traslacin de la misma, puesto que la
estructura se analiza sobre resortes elsticos que se deforman en traslacin y rotacin.

XIV Congreso Nacional de Ingeniera Estructural

Acapulco, Gro., 2004

MODELACIN DE LA INTERACCIN CON ECOgcW

Los siguientes pasos describen la forma en que est incorporada la interaccin suelo-estructura en ECOgcW.
1.

Se crea la geometra bsica del modelo de la manera usual, teniendo en cuenta que todos los nudos
de apoyo debern estar al mismo nivel y empotrados.

2.

Se definen los parmetros del anlisis ssmico, que puede ser esttico o dinmico modal. Al
seleccionar las NTCDS versin 2004, aparece en la parte inferior de la ventana los botones para
definir el suelo y las caractersticas de la cimentacin.

Figura 4. Parmetros anlisis ssmico

3.

Si se tiene un estudio del suelo para el sitio, se capturan los valores de cada estrato para calcular el
periodo del sitio. Se pueden definir tantos estratos como sea necesario, la ventana agrega renglones
conforme se capturan los estratos. El periodo del sitio lo calcula el programa siguiendo el
procedimiento establecido en las NTCDS. La profundidad de los depsitos firmes, H s , la determina
como la suma de los espesores de los estratos definidos.

Figura 5. Datos de suelo

Si no se dispone de un estudio del sitio, se pueden obtener estos valores de los mapas de isoperiodos
e isoprofundidades de las normas y capturarlos en la ventana que recibe los datos de la cimentacin.
4.

Se definen los datos de la cimentacin, la cual puede ser somera, sobre pilotes o mixta.

Sociedad Mexicana de Ingeniera Estructural

Figura 6. Cimentacin somera

5.

Para los tres tipos de cimentacin, se captura el enterramiento, la relacin de Poisson del suelo y el
amortiguamiento histertico.
Si no se captur la estratigrafa del punto 3, se capturan en esta ventana los valores de periodo del
suelo, profundidad de los depsitos firmes y peso volumtrico medio del suelo.

6.

Si la cimentacin es somera, se captura la informacin de cada placa de apoyo (losa de cimentacin o

zapata) que forme la cimentacin.
Los datos que se capturan son las coordenadas de su centro geomtrico ( Xcg,Ycg ), sus dimensin
paralela al eje X y al eje Y ( dX , dY ) y las caractersticas de la placa en direccin X y Y . Para
este ltimo campo se captura L si la dimensin es larga y C si es corta. As, una zapata aislada
tendr CC en este campo, una zapata corrida con dimensin larga paralela a X tendr LC y una
losa de cimentacin tendr LL . Esta ultima informacin la emplea el programa para determinar los
radios equivalentes.

7.

Si la cimentacin es a base de pilotes, se captura su mdulo de elasticidad y se selecciona el tipo de

pilotes, de friccin o de punta.
Para cada pilote se captura su centro geomtrico ( Xcg,Ycg ), dimetro ( d ) y longitud ( Lg ). La
ventana ir abriendo renglones conforme se capturen mas pilotes.

Figura 7. Cimentacin sobre pilotes

8.

Si la cimentacin es mixta, se captura la informacin descrita para cimentacin somera y sobre

pilotes.

XIV Congreso Nacional de Ingeniera Estructural

Acapulco, Gro., 2004

Figura 8. Cimentacin mixta

9.

Se selecciona la opcin Interaccin suelo-estructura en la ventana de anlisis ssmico.

10. Todos los dems datos (cargas, combinaciones, etc.) se definen con el mismo procedimiento que se
sigue para los edificios sin interaccin.
EJEMPLO DE APLICACIN
Se revis un edificio de 20 niveles, con entrepisos de 3.30 m, 5 crujas de 6 m en el sentido X y 4 crujas de
5 m en el sentido Y . Todas las columnas tienen una dimensin de 80x60 cm y las trabes de 35x70 cm.
La carga muerta se consider de 650 kg/m2 ms el peso propio de los elementos. El peso propio de los
elementos equivale a una carga de 455 kg/m2. La carga viva se consider de 250 kg/m2, salvo en azotea que es
de 100 kg/m2.
La combinacin para determinar las masas para el anlisis ssmico se consider con el 100% de la carga
muerta mas el 70% de la carga viva.
El edificio se consider ubicado en una zona cercana al cruce de Calzada de Guadalupe y Circuito interior,
teniendo un periodo predominante del suelo Ts = 1.1 s, de acuerdo con el mapa de isoperiodos de las normas
y una profundidad de los depsitos firmes profundos de 25 m.
El las figuras 10 y 11 se muestran los valores capturados para cimentacin somera y cimentacin a base de
pilotes.

Sociedad Mexicana de Ingeniera Estructural

Figura 10. Datos cimentacin somera

E
6

D
5

2
1

Figura 9. edificio analizado

Figura 11. Datos cimentacin a base de

pilotes

XIV Congreso Nacional de Ingeniera Estructural

Acapulco, Gro., 2004

El edificio se analiz con los mtodos esttico y dinmico modal, considerando a la cimentacin formada por
el cajn, los pilotes y la combinacin de ambos elementos (mixta). En la tabla 1 se resumen los valores que
resultan para los cortantes de base rgida y flexible, en las dos direcciones de anlisis, as como para el factor
~
Vo Vo antes de tener en cuenta los lmites reglamentarios.
Tabla 1. Resumen de valores obtenido para edificio analizado

Q =1

Mtodo/cimentacin

Q =2

V x (ton)

V y (ton)

V x (ton)

V y (ton)

Esttico/rgida
Esttico/cajn

2404
1639
0.682

3067
1697
0.533

1189
944
0.794

1513
1017
0.673

Esttico/pilotes

1729
0.719

1727
0.563

980
0.824

1033
0.683

Esttico/mixta

1949
0.811

2076
0.677

1058
0.890

1195
0.790

Dinmico/rgida
Dinmico/cajn

2053
1365
0.665

2527
1326
0.525

1016
800
0.788

1246
806
0.647

Dinmico/pilotes

1444
0.703

1362
0.539

837
0.824

827
0.664

Dinmico/mixta

1629
0.794

1613
0.638

901
0.887

959
0.770

~
Vo Vo
~
Vo Vo
~
Vo Vo
~
Vo Vo
~
Vo Vo
~
Vo Vo

Las principales observaciones que se desprenden de la interpretacin de estos resultados son las siguientes:
1.

2.

Los mtodos de anlisis esttico y dinmico modal arrojan resultados similares que confirman, por
un lado, la precisin del procedimiento simplificado reglamentario y, por otro, la validez del modelo
computacional tridimensional implantado en ECOgcW.
~
La interaccin es benfica en todos los casos, Vo Vo < 1, siendo ms favorable para comportamiento
elstico, Q = 1, que para comportamiento inelstico, Q = 2. Esto es congruente con los resultados de

la fig. 2, puesto que el periodo de la estructura ( Tx = 2.126 s y T y = 1.892 s) es mayor que el periodo
del sitio ( Ts = 1.1 s)
3.

Independientemente del mtodo de anlisis, la cimentacin mixta (cajn y pilotes) es la ms rgida,

causando los menores efectos de interaccin; mientras que la cimentacin de cajn es la ms flexible,
causando los mayores efectos de interaccin.

Cabe sealar que no se ha tomado en cuenta el efecto de grupo de pilotes. Este efecto es debido a la
interaccin dinmica entre pilotes a travs del suelo. Actualmente se encuentra en proceso la implantacin en
ECOgcW de un eficiente procedimiento de superposicin propuesto por Dobry y Gazetas (1988) para pilotes
de friccin, usando factores de interaccin dinmica que reflejan la influencia de un pilote activo sobre otro
pasivo. Por tratarse de un problema dinmico, la rigidez de grupo puede ser menor o mayor que la rigidez de
conjunto aqu empleada.

10

Sociedad Mexicana de Ingeniera Estructural

COMENTARIOS FINALES
Se ha descrito la forma de modelar la interaccin dinmica suelo-estructura usando el programa de cmputo
ECOgcW. Si se opta por el anlisis esttico, se aplica directamente el procedimiento simplificado
reglamentario. Pero si se recurre al anlisis dinmico modal, se analiza la estructura tridimensional sobre
resortes elsticos en sustitucin del suelo. Con ello se toma en cuenta los efectos de interaccin en los
periodos y modos naturales de vibracin. Para considerar los efectos en el amortiguamiento y la ductilidad, se
hace uso de espectros de sitio y factores de reduccin de resistencia modificados por interaccin slo en el
modo fundamental.
El procedimiento se ilustr con un ejemplo numrico del anlisis ssmico de un edificio tpico. Se evaluaron
los efectos de interaccin aplicando los mtodos de anlisis esttico y dinmico modal. La comparacin
satisfactoria entre resultados confirm, por un lado, la precisin del procedimiento simplificado reglamentario
y, por otro, la validez del modelo computacional tridimensional. Asimismo, se evalu el efecto del tipo de
cimentacin, encontrndose que la cimentacin mixta es la ms rgida (menores efectos de interaccin),
mientras que la cimentacin de cajn es la ms flexible (mayores efectos de interaccin).
REFERENCIAS
Avils J y Prez-Rocha L E (2003), "Soil-structure interaction in yielding systems", Earthquake
Engineering and Structural Dynamics, Vol. 32, pp. 1749-1771.
Dobry R y Gazetas G (1988), "Simple method for dynamic stiffness and damping of floating pile groups",
Geotechnique, Vol. 38, pp. 557-574.
Corona G. (2002), ECOgcW, "Programa para el anlisis y diseo de edificios de concreto".
Luco J E y Wong H L (1987), "Seismic response of foundations embedded in a layered half-space",
Earthquake Engineering and Structural Dynamics, Vol. 15, pp. 233-247.
NTCDS-RCDF (2004), "Normas tcnicas complementarias para diseo por sismo", Reglamento de
Construcciones para el Distrito Federal, Gaceta Oficial del Departamento del Distrito Federal.
Resndiz D y Roesset J M (1986), "Soil-structure interaction in Mexico City during the 1985
earthquake", International Conference on the 1985 Mexico Earthquakes, Factors Involved and Lessons
Learned, ASCE, New York.

11