INTERACCIÓN DEL SUELO, CIMIENTO Y ESTRUCTURA:

EL CASO DE LAS ZAPATAS (la PARTE)

(INTERACTION OF SOIL, FOUNDATIONS AND STRUCTURE: THE CASE OF ISOLATED
FOUNDATIONS, PART 1)

Manuel J. Freire Tellado, Arquitecto Fecha de recepc ión: 9· V[[-2000

Opto . de Tecnologia de la Construcción. Univ . de A. Coruña 400-36

RESUMEN SUMMARY

El artículo defiende la mejora de los modelos de cálculo de This paper proposes the improvement ofthe structural analysis
estructuras introduci endo el conjunto suelo- cimiento. models introducing both soil and foundation .in the modelo
Seguidamente se plantea la modelización del terreno, Next the soil modelling is analyzed, using the model ofthe
recurriendo para ello al método del módulo de balasto, coefficient ofsubgrade reaction for this task. lts application to
explicando su aplicación al cálculo matricial. the matrix analysis is explained, too.

En la segunda parte del artículo se aplica este modelo a una At last, in the second part, the proposed model is aplicated to a
serie de pórticos planos correspondientes a diversos tipos number ofplane frames belonging to various structural types
estructurales con zapatas aisladas y se estudian las with isolated f oundations. The results derivedfrom the
consecuencias de la variación de determinado s parámetros. variation ofcertain parameters are carefully studied. The
Las conclusiones obtenidas resultan especialmente conclusions concerning industrial premises frames and
sorprendentes en lo referente a los pórticos tipo nave ya las centering beams are specially surprising.
vigas centradoras.

Introducción Planteamiento del problema.

El método del coeficiente de balasto.

Este artículo es el primero de un conjunto de dos en los que Aplicación a la formulación matricial.

se analizan los resultados de diversos modelos de cálculo de Tipos estructurales y método de análisis.

estructuras en los cuales se ha incorporado la representación Anejo: estudios de interacción.

del terreno y de la cimentación conjuntamente con la parte Bibliografia.

del esqueleto resistente que tradicionalmente se conoce

como estructura -o, en algunos casos, como superestructu­
ra- . Se recurre a la división en dos partes para mantener la Segunda Parte
extensión del texto dentro de los límites de un artículo,
partes que se corresponden con la exposición de las bases Análisis desarrollado :
del método, parte primera, reservando la segunda al relato
de los análisis efectuados y a las conclusiones que se Pórtico tipo colegio sobre viga flotante .
derivan. El índice general del conjunto es el siguiente: Pórtico metálico de 20 m de luz sobre zapatas.
Estudio de la zapata de medianería .
Primera Parte
Conclusiones.
Estructuras de edificación, cimientos y terreno. Bibliografia.
46
Informes de la Const rucc ión, Vol. 52 n" 47 1, enero /feb rero 2001

Estructuras de edificación, cimientos y terreno comportamiento de un tipo apoyo del otro son realmente
pequeñas. Elconj unto sue lo-c imentación repercute de forma
El proceso de análisis convenciona l del esque leto resis­ cons iderab le en las leyes de esfuerzos de la superest ruct ura,
ten te de un a edificación empieza por el estudi o de la es­ al tiempo que [a ci men taci ón es, ademá s, un facto r de coste
tructura, no rma lmente supuesta emp otr ad a en la cimen­ mu y important e en el conjunto del ed ificio. Co nseguir que
tac ión . Posteri ormente, una vez obtenidas las acc iones la obra ejecutada se comporte de la form a en qu e ha sido
que el edificio transmite a la ci me ntació n, se reali za e! cons ide rada en cá lcu lo no es en ningún caso un tema meno r,
dim ensionado de ésta, compro ba ndo, en primer lugar, que ye n el que el com porta m ien to de l suelo influye de forma
las tens ione s tra nsmitidas sea n admi s ibles para el terre no y apre ciab le.
rea lizando, seguidamente, [as comprobaci on es pe rtinentes
a la estructura de cimentaci ón. Las etapas señal ad as son Obviamente en el plan teamiento anunciado subyace también
estancas: la cimentación recibe de [a superestructura unos la necesidad de introducción de precisión desde e! punto de
esfu erzos que ha de aceptar y del sue lo una condición de vista conceptual, mejorando la fiabilidad de los resultados
tensi ón admisible. obtenidos, por cuant o la co ns ideración de un ma yor número
de parámetros permite reflejar la situación real de forma
Frente a uno s sistemas de análisis de la superestru ctura que m ás fided igna, al tiempo q ue permitiría la reconsider ación
inco rpo ran elementos mu y sofisticados -rnatriciales 3D, de los coeficientes de seg ur ida d, co n el consigu iente ah orr o.
elem entos finitos, co rrec cio nes por efe cto s no lineales,
efe cto P-ó-, por ejemp lo, o incluso cálculos en teorí a de Factor de interés añadido es el estud io de situ acio nes en las
segundo orden- en los que se analizan numeroso s efectos qu e se pr odu cen desvi aciones en tre el modelo de cálculo
sec undario s, resu lta so rp rende nte la radical s implifica ción con ve ncional y la realid ad, como puede ser el caso del
de las cond iciones de co ntorno al supo ner -en el ám bito de equi libr io en za pa tas de medianería y esquina. De tod o
la edi ficaci ón convencional- , la estruc tura empotrada -en técnico con cierta exp eri en cia en ejecución de obra resultan
ciertos cas os articulada- en su en lac e co n la cimentació n. Si conocidas las protest as debidas a las dimensiones y arm a­
todo el análisis estructural se basa en la co mpatibilidad de dos de las viga s ce ntrado ras. Protestas éstas que no tendrían
despl az amientos entre los diferentes elemento s que mayor repercusi ón si no ex istiesen numerosos ca sos de
co mpo nen la estructura, al llegar a cimentación esta premi­ zapatas de med ianería sin nervio -y a las que, además, no
sa se olvida y el tratamiento de las partes es totalmente acometen pilares de rigid ez sufici en te como para ga ranti­
independien te: la tensión adm isible se su ele tomar igua l zar su equ ilibr io de acuerdo con los estudios teóricos- y
para todas las pi ezas sin tener en cuenta que ésta de pe nde que, sin embargo, son perfectam ente estables. Esta situación
tam bién de las dimensione s de aqué llas y se obvian , ha lle vado a algunas casas de so ftwa re a proponer valores
as im ismo , tant o la defonnab ilidad de la cimentación co mo de 1, I O para el coe fic iente de mayoración de acc io nes en
el efecto de la compresib ilidad de! sue lo y, por tant o , de los este tipo de elementos . Si b ien este val o r no resulta aco rde
as ien tos de la edifi ca ción . con la normativa vige nte , no se puede olvidar que la Teor ía
de Estru cturas es, pes e a su s importantes re cursos
El plante amiento anterior, con una va lidez no ob stante matemático s, un a cienc ia de carácter experimental.
rec onocida, tiene su origen en épo cas pasad as, con recursos
de cálcu lo fundamentalmente manuales , en los que el La cimentación superficial es además la situación que
aná lisis por partes de la estructura era la únic a forma introduce las ma yores divergencias con el modelo de
factible de abordar su tratamiento . Sin emb argo, los recur­ cálculo vigente : s i en el ca so de un a cimentación profunda
sos de cá lculo disponibles ho y en día, hacen viable el los asie ntos son pr áct icame nte nulos -se sue le buscar un
estud io co njun to de toda la est ructura . Ho yes fac tible estrato rígido- y la rigidez de los en cep ado s y los elementos
abordar el tratamiento ana líti co gl obal de l conj unto estructurales inte rm edios entre p ilares y pilotes garantizan
estru ctur a-cimiento-suel o, reco giendo las interaccion es bastante bien la hipótesis de em potram iento, en e l caso de
entre las distint as partes del esque leto resistente. Es más , las cime ntac iones super fic ia les - especialmente cua ndo el
tratamientos de este tip o son norma usu al en el es tud io de s istema eleg ido es de zap atas ais ladas -, la rotación de las
pre sas y grandes obras de inge nier ía. zap at as y lo s as ie ntos - tot a le s y d iferenci a le s­
experime ntados, de bidos ambos a la compresibilidad de un
El prob lema es más cercano de lo q ue pud iera parecer: terreno con variaciones muy locales, dificultan [aaceptación
la bib liografía da cuenta de cómo la deform aci ón del del sistema de cá lculo señalado
terr eno puede hacer que lo que se ha considerado en el
modelo de cálculo un empotramiento se comporte en Planteamiento de problema
realid ad co mo algo mucho más parecido a una artic ulación.
La modifica ción de la condiciones de bo rde lleva apa rej ado Si bien se acaba de seña lar q ue posib leme nte las mayore s
en ciertos cas os -pórticos tipo portería po r ej emplo- un discrepanc ias en el aná lis is conj unto sue lo-c imi ento­
no tab le camb io en las ley es de es fuerzos , debido a qu e la estructura se produ zcan en las ci mentaciones supe rficia les,
magn itud de las rotac iones de las piezas qu e separan el pese a la existencia de un número co ns ide rable de estudios
 47
Informes de la Co nstrucció n, Vol. 52 n" 471 . en ero/febrero 2001

sobre la interrelación suelo-cimienta-estructura, tanto de especialmente los que analizan modelos de suelos distintos
tipo plano como espacial -véanse las referencias incluidas del desarrollado por Winkler , modelo que si bien cuenta
en el Anejo 1-, la totalidad de los trabajos realizados se con un ciertoavalexperimental referidoalconjuntoedificio­
aplican a los casos de cimentación por vigas flotantes -los cimentación-terreno, éste es lamentablemente bastante
análisis planos- y cimentación por losas -los estudios en reducido.
3D-. No se han encontrado referenc ias de estudios de
interacción con sistemas de cimentación formados por No cabe olvidar además la ausencia de capacidad operati­
zapatas aisladas conectadas mediante nervios, la tipología va práctica de alguna de las modelizaciones consignadas en
de cimentación más usual en nuestro entorno. Si bien la la bibliografía, por cuanto la complej idad de su desarrollo
razón de esta ausencia es bastante obvia -ya que tanto vigas las relegan a la esfera de la mera especulac ión teórica, al
flotantes como losas de cimentación son soluciones que se menos en la actualidad, aunque el espectacular desarrollo
emplean con situaciones más criticas, bien por la escasa informático hace muy aventurado avanzar cualquier tipo
capacidad portante del terreno bien por la elevada magni­ de crítica en este sentido.
tud de las tensiones transmitidas-, la carencia de estudios es
especialmentedeplorablepor cuanto locomún de lasolución El objetivo de estas líneas es proponer y valorar los
de zapatas aconseja un mínimo de optimización . Desde esa resultados de un modelo de cálculo de estructuras de
perspectiva, uno de los puntos tradicionalmente conf1icti­ edificación que tenga una mejor aproximación a las
vos en este tipo estructural -y que se desea analizar- es el condiciones de borde reales de las estructuras. Este modelo
funcionamiento de las vigas centradoras, su necesidad y se ha construido con un método que -se espera- resulte
comportamiento, desde la constatación experimental de la operativo y consecuente con el gran número de elementos
capacidad de múltiples estructuras para mantenerse en pie y secciones presentes, esto es, compatible con la demanda
pese a carecer de estos elementos . de economía que impera en toda actividad. El planteamien­
to busca analizar la repercusión de la modelización del
El análisis bibliográfico ha permitido constatar, asimismo, terrenoenel dimensionado de loselementos decimentación,
que las dimensiones de los elementos constitutivos de los lo que supone una revisión de los métodos clásicos de
modelos distan bastante de ser las usuales en edificación en dimensionado de éstos.
nuestro país, empleándose escuadrías que resultan chocan­
tesen nuestroentornoy que, deeste modo,pueden contribuir El método del coeficiente de balasto
a enmascarar la comparación intuitiva del comportamien­
to del modelo analizado con los resultados proporcionados Para ello el primer problema a resolver es la modelización
por situaciones más comunes. del terreno. Se optó por el método del Módulo de Balasto
o modelo de Winkler [O 1], que parte de la hipótesis de que
Ademáses norma común en lostrabajoscitadosel despreciar el asiento producido en un punto mantiene una relación
el comportamiento del material estructural, que se supone lineal con la presión a la que está sometido el terreno,
elástico y lineal. Empero, debido a lo reducido de las transmitida por la cimentación . Analiticamente,
escuadrías empleadas , la estructura se encuentra sometida
a tensiones muy altas, que tal vez no permitan una
simplificación tan radical.
donde p representa la presión transmitida al terreno,y -ó s­
Los coeficientes de seguridad usualmente empleados en el asiento experimentado y k, el factor de proporción entre
mecánica del suelo se traducen en tensiones que rara vez ambos,conocido por variosnombres siendo el másaceptado
superan el tercio de la carga de hundimiento del terreno. Se Coeficiente o Módulo de Balasto, con unidades de kg/cm',
da en ciertos casos la paradoja de adoptar complejos
modelos de terreno, mientras que éste está sometido a unos De acuerdo con la ecuación inicial, p = ks ' y , para los casos
valores tensionales reducidos que no justifican tamaño de carga lineal q repartida sobre un ancho b y para carga
esfuerzo. Sobre todo cuando en la superestructura se concentrada P actuando sobre una placa de área A, la
presentan fuertes valores tensionales que se tratan, sin formulación se completa como sigue:
embargo, como situación elástica y lineal. Parece sensato
investigar el comportamiento real de la estructura y te­ q = k;b y = ky P = kA = K y
rreno, verificando las solicitaciones a que realmente están
sometidas las piezas y su congruencia con los modelos donde las unidades de k y K son variables para permitir la
usuales, dentro de un cierto grado de igualdad en las homogeneidad de la fórmula.
simplificaciones asumidas en cada caso y antes de cualquier
intento de actuar sobre los coeficientes de seguridad . Si bien éste es un método que no se ajusta exactamente
a la realidad experimental del terreno, con problemas
Es también notable la ausencia de constatación de resulta­ sobradamente conocidos, los resultados que histórica­
dos experimentales de todos los trabajos consultados, mente ha proporcionado están avalados por la experiencia,
48
Info rmes de la Const rucc ión, Vol. 52 n? 47 1, enero/fe brero 200 1

factor de considerable importancia en un análisis como el permite ignorar la superposición de los distintos bulbos de
que se aborda en estas líneas. Además, dado el ámbito de presiones de los cimientos del edificio -o de otros cimientos
este estud io, métodos más preciso s como el modelo existe nt es- prop or c ion ando un os result ad os que
hiperbólico -el ástico no lineal- no aportan mayor precisión históricamente han demostrado un grado de fiabilidad
en lo s re sultados, por cua nto que en edificación aceptable. Fiabilidad que se explica -al menos parcialmen­
convencional, las tensiones transmitidas al terreno no son te- porque, como es sabido, a una profundidad de ¡,5'b ­
altas: la tensión admisible engloba un doble factor de dondeb es la dimensión delcirniento-, la tensión transmitida
seguridad suficiente frente a hundimiento y frente a asientos al terreno es del orden de O, l (J', de forma tal que para las
excesivos . Por esta razón la relación reacción del suelo/ dimensiones de cimientos usuales en cimentación no suele
asiento se sitúa en la parte baja de la rama, que es, en todos haber problemas de superposición.
los modelos propuestos, bien una recta o bien sustituible
por una recta sin error apreciable, tal y como se puede Si bien la formu lación del método se debe a Winkler, el
apreciar en el Gráfico 1 [03]. Método del Módulo de Balasto es conocido también por el
nombre de Método de Winkler -ya que fue su introductor
en 1867 [01]- Y sirvió de base para el clásico estudio de
Gráfico 1
(VALUACION DE LOS COEFICIENTES DE BALASTO
Zimmermann sobre el análisis de rieles de ferrocarril sobre
.fl. .fl. p
traviesas en 1888 (02],recibió un impulso fundamental con
u
o .3 ?
,,
p los trabajo s de Terzaghi.
---..... , i ,,
)

,, )

,
) De acuerdocon lahipótesisgeneraldel método, en pequeñas

"" "'
)

.........

<,
)

,
) superficies y para placa circular se introduce la hipótesis de
.........

<, ......... ' .... ,

)
)
que el coeficiente de balasto es inversamente proporcional
-, ~ <; al diámetro de la placa empleada:
,
\
) <, , ..... A
)
)

,
)

1\ )
)
1
donde d y d 'son los diámetros de las placas de ensayo y k b
\
)

,
)
Ykb , son los módulos de balasto correspondientes a éstas.Se
I
admite igualmente que el valor correspondiente a una placa
\
I
,
I
) de carga cuadrada equivale al de una placa circular de igual
)
)
área. Esta hipótesis admite implícitamente que, para del
I
I
I
rango de dimension es de las placas de carga -de 30 a 75 cm­
CI
)
el asiento produ cido es proporcional al tamañ o de la
)
)
cimentación.
Y. •
RELACION ENTRE LA REACCION DEL SUELO. p. y EL ASIENTD, y:
LíNEA DA: RELACiÓN SUPUESTA , CURVA DC: RELACiÓN REAL
Las relaciones anteriores son especialmente útiles para
El método presenta además la ventaja fundamental de su correlacionar resultados experimentales, normalmente
sencilla implementación informática dado que caracteriza derivados de ensayos realizados con lasplacas de diámetros
el terreno exclusivamente por el valor de un parámetro, el normalizados 34 y 75 cm o con pla ca cuadrada de 30 cm
módulo de balasto, independientemente de la formulación de lado, de tal forma que la expresión anterior se convierte
que permita su obtención y que varía con el tipo de terreno. en 34 k)4 = 75 k75 para las placas circulares, admitiéndose
además que r; = k)O, siendo este último el valor para la placa
Dentro de las carencias del método se pueden señalar que cuadrada de lado 30 -lo que implica tener en cuenta la for­
el módulo de deformación en terrenos es creciente con la ma de la cimentación-o
profundidad, mientras que el método supone unacorrelación
lineal entre el asiento producido y la tensión con la que A partir de los resultados anteriores, Terzaghi [03] aproxi­
reacciona el terreno. Olvida, además, considerar el efecto mó la variación del asiento ocurrid o en las zapatas a partir
de la preconsolidación del terreno, en el que para cargas del experimentado por una placa de carga cuadrada de 30
inferiores a la de consolidación elasiento sería despreciable. cm mediante las expresiones siguientes -reescr ítas para
De hecho los dos rasgos que acabamos de señalar no son adaptarlas al sistema métrico decirnal-:
más que particularidades de la hipótesis fundamental del
método que considera el terreno co mo un fluido, b
despreciandoparámetros tales como el ángulo de rozamiento s = s . - '- (arcillas duras)
e )0 30
interno y la cohesió n. Esto conduce a la hipótes is
fundamental del método: solamente el terreno bajo el l b
cimiento experimenta deformaciones. Y esto, que es una s =s ( - -'- F (arenas)
e )0 b +30
falacia conceptual, tiene sin embargo la importancia que e
 49
Informes de la Construcción, Vol. 52 n" 471, enero /febrero 200 I

donde s., es el asiento experimentado por la placa de 30 cm, - Proporción entre las dimensiones del cimiento y cargas

se es el asiento del cimiento y be es el ancho del cimiento. aplicadas.

- Rigidez Relativa del cimiento.

De acuerdo con la ecuación inicial y dado que la presión - Nivel de Cargas: Cargas Vivas/ Cargas Muertas

transmitidaes constante, es inmediato obtener lasecuaciones

clásicas propuestas por Terzaghi en 1955 [03]:
Aplicación a la formulación matricial
k = k . ~ (arcillas duras)
e JU b l'

e
En las figuras del Gráfico 2 se muestra el proceso básico de
idealización [08] del terreno desde la situación real hasta el
b + 30
ke = k)O ( e F (arenas) esquema de cálculo simulando el efecto del terreno median­
2·b e te un conjunto de muelles .
Posteriores estudios e investigaciones han generalizado y
reescrito parcialmente estas expresiones, de tal forma que Básicamente este proceso consiste en la sustitución del
sus expresiones más conocidas resultan ser [04] continuo del suelo por un conjunto de apoyos elásticos
suficientemente próximos . Para ello se divide el terreno en
una serie de fragmentos de ancho igual al del cimiento -en
ke = k o' fb (arenas sueltas o muy sueltas, limos, arcillas)
consonancia con las hipótesis del método del módulo de
balasto- y paso a conveniencia, sin temor a cometer un
e
grave error si se asegura que la discretización sea lo
suficientemente tupida. De esta forma se ha idealizado la
pieza, realmente apoyada sobre un medio elástico continuo,
b e + bO
ke = ka ( - - - - - F (arenas densas) como un elemento apoyado en un número finito de apoyos
2·be elásticos. Seguidamente se evalúa la capacidad resisten­
te de cada una de las porciones de terreno así definidas y se
representa porun muelle -esto es, un elemento con reacción
be + b O exclusivamente vertical- o

3·b e El paso siguiente es introducir los muelles así definidos en

la formulación matricial, consiguiendo que la respuesta del
donde b o representa el lado de la placa de carga. terreno sea exclusivamente vertical. Esto se puede realizar,
si el programa lo permite, introduciendo la constante de
Otros autores -[05], pág. 207; [06], pág. 310- recogen muelle para cada apoyo . Pero, en nuestro caso, el programa
además la última de las expresiones de Terzaghi citadas, disponible, un matricial convencional, no ofrecía esta
generalizándola para un cimiento rectangularde proporción posibilidad. Por ello se ha tenido que representar cada uno
y dimensiones cualquiera, expresándola en el sistema de los apoyos elásticos por una barra, asegurando la reac­
internacional de unidades ción de ésta en dirección vertical. Estas barras se conocen
como barras ficticias.

2'a + ¡ Como es sabido, en la formulación matricial una barra

3'a queda representada por su matriz de rigidez. De esta matriz
depende el comportamiento resistente de la barra. Dos
factores resultan determinantes: la rigidez axil de la pieza
le y la rigidez a flexión .
siendo a = (arcillas preconsolidadas)
be . de e 11 as nene
. .. di' EA- ; la
La primera una expresión e tipo -
l
E'¡
Expresiones que, sin ser las únicas ni las más precisas, sí son segunda es de la forma - - si se desprecia la influencia
l
las más conocidas para el cálculo del módulo de balasto del
de la deformación por cortante -situaci ón usual en los
cimiento a partir de los valores obtenidos en ensayos de
placa de carga. De hecho, hoy en día, existen expresiones programas de cálculo-, con los siguientes significados
más precisas que hacen depender el Coeficiente de Balasto
de [07]: E: Módulo de Young del material de la barra ficticia que
emula el muelle
- Módulo de Compresibilidad E, bajo el cimiento . A: área transversal de la barra ficticia que emula el muelle
- Forma y dimensiones de la cimentación -rectangular,
cuadrada, circular. ..­ 1: longitud de la barra ficticia que emula el muelle
50
Informes de la Construcción, Vol. 52 n° 471 , enero/febrero 2001

Gráfico 2 Para garantizar que la reacción se produce en dirección

vertical se debe anular la resistencia a flexión de la barra,

esto es, -E,¡- = 0, SI, bien

, la Igualdad
' es irnposiible de
l
representar sin alterar la formulación matricial, sin cmbar­

E1
go es senci'11o conseguir' que ---== 0, lo que produce
l
los mismos efectos prácticos. Para ello basta asignar a la
barra un momento de inercia suficientemente reducido. De
esta forma se anula la respuesta a flexión de la barra,

Como además en la formulación matricial el cortante es

SITUACiÓN REAL x

función de expresiones tipo E1 y EJ1 , con lo que e l

{2 l
artificio anterior anula también su efecto. Obviamente se
debe tener la precaución de introducir como longitud de
estas barras ficticias valores superiores a la unidad,

En estas condic iones, al producirse el asiento Yv en un punto

del cimiento y de acuerdo con la formulación matricial, la
reacción de la barra en contacto con ese punto toma el valor

EA
l
valor que debe ser igual a la reacción del suelo producida
DISCRETIZACIÓN DEL TERRENO
por cada columna virtual de terreno -evaluada con el
método del módulo de balasto-, que para el mismo
desplazamiento vertical en ambos casos se expresa como

y -Ke -bx
donde

ke : Módu lo de balasto del cim iento.

b: ancho del cimiento
x: ancho de la columna virtual de suelo considerada .

Igualando ambas expresiones

Dado que algunos programas admiten exclusivamente un

único valor de E -y, en todo caso, por comodidad- y para
facilitar la introducción de datos, suelen prefijarse los
y
valores del módulo de Young E y la longitud l de las barras,
por lo que resulta un valor para el área de las barras ficticias
de
kb
e e
'x" l
E

El sistema anterior permite trabajar con anchos de cimiento

variables y distintos módulos de balasto -lo que posibilita
la emulación de variaciones del terreno bajo el cimiento o
la introducción de diversas correcciones al valor del módulo
ESQUEMA DE CÁLCULO de balasto-.
 51
lnforrnes de la Cons trucc ión, Vol. 52 n0 471, enero/fe brero 2001

Interesa, además, que la longitud 1sea de cierta magnitud, El análisis de la bibliografía, en la que se ha constatado su
al efecto de que aumente la flexibilidad transversal de las exclusión, así como la práctica constructiva usual de Galicia
piezas. Empero, cuando la flexibilidad transversal es muy -realidad de fondo en el que se inscribe este estudio- han
elevada es preciso añadir una barra horizontal ficticia para llevado a las tipologías propuestas en este trabajo, que
evitar problemas de mal condicionamiento de matrices y quizás pa r ser muy comunes han sido muy poco estudiadas,
errores en los resultados. Esta barra existe realmente -la debido a su falta de singularidad, y que versan sobre
acción del terreno sobre el borde libre de la pieza de s ituac iones de ci men taciones de tipo superficial,
cimentación y el rozamiento cimiento suelo- pero tampoco particularmente zapatas.
debe preocupar su introducción, por cuanto realmente las
so licitaciones a las que se ve sometida en los casos usuales Así, se ha estudiado previamente una serie de vigas flotan­
son de muy reducida magnitud. tes -correspond ientes a una estructura típica de colegio­
según modelizaciones conocidas [10], introduciendo las
Para un cálculo ajustado es importante la adopción de un cargas como puntuales en un primer paso, posteriormente
módulo del balasto del cimiento adecuado, teniendo en considerando el refinamiento de considerarl as repartidas
cuenta las diversas correcciones que plantea la bibliogra­ sobre elancho de los pilares y, por último, introduciendo la
fía, Entre éstasson de especial interés aquellas correcciones superestructura en el modelo. Este tipo estructural se utilizó
que garantizan cierta continuidad en la deformación del como comprobación del modelo desarrollado.
soporte -suelo-, corrigiendo los valores del módulo de
balasto en los extremos del cimiento o empleando técnicas Se ha pasado a continuación al estudio de un pórtico
de muelles emparejados - 'c oup led springs 'o, de acuerdo metálico de 20 m. de luz apoyado sobre zapatas, que podría
con lo expuesto en la norma ACI-336, u otros sistemas ­ ser la estructura de una nave, para analizar el conocido
véase [07J-. problema de su cimentación a base de zapatas en suelos de
no gran calidad.
Tipos estructurales y método de análisis
Posteriormente se afrontó el estudio de las zapatas de
El análisis que se expone en la segunda parte de este artículo medianería, para el que se analizaron dos casos; un pórtico
ha sido efectuado mediante una serie de cálculos de las de bajo y 2 plantas altas en el que la carga del cerramiento
estructuras propuestas sobre un programa matricial. Dicho llegaba hasta el terreno sin afectar a la estructura -caso de
programa permite el control de la exactitud del cálculo que el cerramiento fuese continuo hasta el terreno- y,
mediante el estudio de los desequilibrios finales en los posteriormente, un pórtico de bajo y 5 plantas, en el que la
nudos de laestructura -condición inexcusable para garantizar estructura recibía, además, la carga de cerramiento.
la validez de la idealización analizada-, así como el enlace
con programasdeCAD que permiteneltratamientoposterior En nuestro caso, por razones de ubicación geográfica, el
de las leyes de esfuerzos obtenidas [09]. estudio se ha centrado en el caso de suelos sin cohesión,
arenosos, con valores de módulo de balasto en placa
El método de análisis partió del estudio de la repercusión circular de diámetro D= 75 cm, adoptando valores !sl= 9 y
que tiene la modificación del valor de una variable mientras k71= 15 kg/cm ', usuales para suelos arenosos pobres o ya
se mantienen constantes los valores de todas las demás. con cierta calidad [11], empleándose las fórmulas de
Este sistema da lugar a una serie de familias de resultados correlación correspondientes a las arenas. De todas formas,
comparativos que permiten analizar la distinta repercusión tal y como se comprobará posteriormente, el valor del
de la modificación efectuada en una serie de casos cruzados, módulo de balasto elegido -en el rango usual de valores- no
estableciendo asimismo la repercusión porcentual de la resulta especialmente determinante en el comportamiento
variación. Evidentemente nadajustificaque nose produzcan del sistema.
interacciones de mayor repercusión cuando se modifican
simultáneamente dos variables, pero, a la luzde los resulta­ Además, en los estudios de pórticos planos sobre zapatas,
dos obtenidos, ésta ha parecido una aproximación suficien­ el cálculo de kc se ha realizado independientemente para los
te. distintos elementos de la cimentación diferenciados, si bien
éstosseencuentranenlazados por nerviosde atado,situación
Establecido el método de cálculo y la forma de realizar la que se ha incluido por considerarla de buena práctica y por
idealización del terreno, queda j ustificar las estructuras el recuerdo de las prescripciones de la POS.
analizadas. A tenor de la j ustificación anterior nos hemos
centrado en algunas de las más comunes en la práctica Para el dimensionado de los elementos -vigas y pilares,
profesional. Se analizaron distintos tipos estructurales, zapatas- introducidos en el modelo se ha partido de la
inicialmente muy sencillos, pero que fueron arrojando experiencia y de tabulaciones publicadas de uso corriente
resultados que por, su discordanc ia con los obtenidos con en edificación.
los mode los de cálculo actualmente admitidos, llevaron a la
consideración de modelos cada vez más complicados. Entre los modelos calculados no se ha incluido ninguno que
52
Informes de la Construcción, Vol. 52 nO471, enero/febrero 2001

incorporase la corrección de muelle s emparejados. Sí se ha expresión válida para una placa cuadrada de 1pie (30 cm),
tenido en cuenta, en cambio, la corrección por distancia al donde 1, es la longitud del cimiento y b, el ancho de éste.
borde, si bien entre los resultados que se transcriben no se
incluyen éstos dado que significaba un comportamiento b
más ventajoso en los casos en que se consideró. Sowers, 1977 [13] k, = ko (_0
b,.
r 0,5'.'::,n '.'::, 0,7

Como se ha expresado anter iormente, la rigidez axil de la Esta expresión , aplicable a depósitos de arena natural
barra ficticia que emula el terreno se expresa como y compacidad relativamente uniforme, introdu ce una
EA corrección exponen cial que parece provenir de diferencias
- I -=ke'b e'x en la carga última encontrada. Ha sido probada con placas
de carga de hasta 1,20 m. El autor indica como corrección
donde x es el paso de la discretización . Evidentemente , experimental que, en cimientos rectangulares, es más exacto
como según lo anteriormente expuesto, siendo x constante, emp lear el lado del cimie nto cuadrado de igual área. Si se
resulta que la capacidad de carga del cimiento es muy poco define el coeficiente de forma de la cimentación rectangu-
dependiente de su ancho b, especi almente en el caso de la I
fórmula propuesta para las aren as súeltas o muy sueltas, lar como CJ. = ---¿-, en estas condiciones y teniendo
b e

limos y arcillas, en la que kc = ka' T- ' lo que, sustituido

e
en cuenta la considera ción anterior la fórmula se rescribe
como
en la expresión anterior, conduciría a b
k = k ( _ _0_)" 0,5'.'::, n'.'::, 0,7
EA b e o be' 'JCl
- = (k . ---"- ) b :x = k -b ' X
I a be ' a a
Según DE BEER [14] Y VOGT [15]:
ya que la capacidad portante del cimiento fuese totalmen­
te independiente de su ancho , lo cual no es obviamente l . Cimentación Circular de área A y rigidez suficiente para
aceptable. En el caso de las formul aciones para arenas distribución uniforme de presiones:
densas y para arenas compactas, si bien las consecuenc ias E
no son tan radicales, se producen, sin emb argo, situacione s k = 1 392 - "
curiosas dado que la aplicación de la formulación conduce e ' fA
a que para cimientos de 1,00x 1,25 Y de 1,00x2,50 les 2. Cimentación Rectan gular a > b y rigidez suficiente para
corresponda el mismo módulo de balasto, mientras que la distribución uniforme de presiones:
presión de hundimiento es muy diferente en ambos casos. E
k = 1 330 - "­
El estudio de las diversas propuestas elaboradas para la e ' ~a . b 2
determinación del módulo de balasto permiten contrastar De acuerdo con el concepto anterior de coeficiente de
las correcciones por razón de form a de la cimentación que forma, la fórmula se puede expresar como
han introducido los diferentes autores. Ya la hipótes is
gener al del método, en pequeñ as superficies , distingue la E,
k , = 1,330
situaciónentre placa circulary placa cuadrada, estableciendo 3
b·¡a.
la relación entre los módulos de balasto de la placa circular DIMITRO V [1 4] E,

y una cuadrada de igual área en relación con su diámetro 25k,= P

equivalente, esto es b (1 _f12)

donde p es un Coeficiente de Forma de la Cimentación
dkJ = k e- ~ 4'A
1t
Rectangular, según la Tabla l .

TABLA 1
La relación anterior para el caso de placas de diám etros a:b 1,00 1,50 2,00 3,00 5,00 10 20 30 50
normali zados 34 y 75 cm y placa cuadrada de 30 cm de
lado, se convierte en 34 kJU = 75 k7j P 1,05 0,87 0, 78 0,66 0,54 0,45 0,39 0,33 0,30

Terzagh i [6, pág. 315], propone una expresión para arci­

E,

llas precomprimidas que tiene en cuenta las dimensio­ DfN 4019: [14] ke = - - - ­
nes relativas del cimiento. Esta expresión ha sido adaptada bfr, O)
por autores posteriores [11] [12], generalizándola para
cualquier ancho de cimiento y expr esándola en el sistema donde la función I depende de las relaciones a/b y z/b,
internacional de unidades como donde z es e1espesor de la capa de terre no efectiva. LaTabla
b . le (arcillas pre­ 2 recoge los valores propuestos para esta función por Kany
2'a + 1 JO SIendo a = ­
k, = k JO Ia b ' be comprimidas) [1 4].
e
 53
Informes de la Construcción, Vol. 52 n" 471, enero/febrero 200 I

TABLA 2

VALaRES DE Ir s. o) PROPUESTOS POR KANY

z:b \ 1.00 1,50 2.00 3,00 5,00 10 ,00 20,00

a:b

0 ,20
° ° ° ° ° °
0, 10 0,10 0, 10 0, 10 0, 10 0, 19
°
0, 19
8 8 8 9 9

0 ,40 0,20 O,JO 0,3 0 0.30 0.30 0 ,34 0,34

9 1 2 3 3

0.6 0 0, 30 0,40 0,40 0.40 0,40 0, 46 0 ,46

7 2 4 5

0 ,80 0, 40
°
0,40 0, 50 0,50 0,50 0,57 0,57
4 7 3 6

1,00 0,40 0,50

°
0, 50 0,60 0,60 0,67 0,67
9 3 7 I 4

1,50 O,SO 0,60 0,70 0,70 0, 80 0 ,86 0.87

8 5 5 1

2.00 0.60 0,70

° 0,70 0,80 0,90 1.00 1.03
4 2 8 5 3

3,00 0,70 0, 80 O,BO 0.90 1,00 1,20 1,25

2 9 9 9

5,00
°
0,70 0,90 1,00 UO 1,20 1,43 1,54
6 3 7

7.00 0,70
° °
0,90 1,00 1,20 1,30 1,55 1,68
9 4 5 7

10,00 0.80 0,90 1,00

°
1,20 1,40 1,68 1,8 5
I 7 9 6 5

20,00 0,80 0, 90 1,10 1.30 1.50 1.89 2.12

2 8 2 2 7

La prop ia norma advierte que se tome un valor z= 2·b si el co mparativa de los factore s correctores por form a propuestos
espesor efectivo de terreno supera 2·b. por los distintos autores.

En el Gr áfico 3 se puede vi sua lizar el valo r de la corre c­ AC!-336 .2R- 88 E

s
ció n propuesta po r Kan y -n ótese que parece haber un [17]: En Come ntarios, kc = - - -"--- -­
error en la primera col umna y que los valores propuestos b ( f- ¡.J)!JF
correspondiesen a b:z en lugar de a z.b- ,
fórmula que se debe a Bow les [1 6] , Y que se obt iene
Para poder co mpa rar la correcc ión po r forma que incl u­ po niendo en rela ción la expresión del módulo con ecuac iones
yen los valores pro pues tos por K any con los pro pues ­ de cálculo de asientos . f represen ta el Fac tor de infl uencia
tos por los otros autor es, se divide n los valo res de cada de forma, de pendiente r de la forma y flexib ilidad de la
fila pare! primerode ellos -el q ue corresponde a un cimie n­ cimentación y que se puede exp resar en función de los
to cuadrado- , en tendiendo q ue los valores de/!s.O) = !p'!1 coeficientes de influenc ia de Seinbrenner, ! Je !l' como
esto es, que el valor resultante es el producto de un coeficiente f -2¡.J
de empotramiento que depende sólo de la profundidad y de
f - ¡.J
un coeficiente de forma -de valor l para la zapata cuadra­
da -oEstas son las curvas que se representan en el Gráfico 5, f F es el factor de empotram iento por la profund idad, deb ido
54
Informes de la Construcci ón, V ol. 52 n" 471, enero/febrero 200 1

Gráfico 3
CORRECCiÓN AL MÓDULO DE BALASTO SEGÚN KAN Y
6

I
I 020

5 I
I I I
I I
I I
I ,
4 I I I
,
I I I I

I I I I FACTOR DE INFLUENCIA DE FORMA Y PRO FUNDID AD, 1,

1,00 r- ~~~---~-~-~~----
1\ I I o
•
I
I
' •
••
• ,
,
"
"
I
I
I
I

.._- ---. --.----... -- .... _------------,---------.- -_... . ----.. -_. ---- -- ----- ­
• , , •• , " I ,

3
1'-1 I O' hn I I

•
"

I • I
"

, I
"

" ,
"

'

I I I 0,90
I
\ I I I I
r- I I I
2 \! I 1 I
\ 1" ,'- -1 I I I
¡-""'I
,\
I ,I
I
,\ r-
~ r-, r-;-. 200 I
~ ::- I I I I I
I

I
1 00 0 0 0 0 o o o
o o o
P = f(S,O) I Lt'l
ci o"
lD r-... CO O'l O
6 ó 6 .....- v· lJ'i o
o a 5 10 15 20 Relación de profundidad , pib
ti
Gráfico 4

Gráfico 5

FACTOR DE MINORACiÓN DEL MÓDULO DE BALASTO

COMPARATIVA DE LAS PROPUES TAS DE TERZAGHI, SOWERS , DE BEER, DIMITROV y KANY

1,00 I

~! I I , I Kanv 06 I
I
I
~~ I I I

i I ;
I
I
Kan 1 I
I
I
­
,
I \~i~~
I ¡
I Kan~ I

-.
I 1
2 I
! I I I

0,50
I
I ~~ ~ ! I I I I I I
Kan 5

¡~ I
I

I I
!
----r=::::: ' I I I
SOWCrs u,
Dim rov
I
I I 1 I-I---+-l I I I I ers , I i I
Sow U, ,

I I
I I

p I Igua Area I

,
I I I
I I I I
I I
a
I I I
5
I I
10 15 20
b
 55
Informe s de la Co nstrucc ió n, Vo l. 52 nO47 1, enero/ febrero 2001

a Fox [17], Y vari able también en funci ón de la forma y En la segunda parte de este artículo se ana lizan distin­
dimensiones de la cimentación . El Gráfico 4, tomado de la tas estructuras de edificación comunes en las que se ha
referencia [18], proporciona una idea clara de la variación representado el apoyo en el terreno med iante el método del
del factor 1F' módulo de balasto descrito en las líne as anteriores.

El Gráfico 5 permite comparar el significado los diversos Anejo: estudios de interacción

factores de corrección propuestos por los autores.
A.l. Estudios de p órticos planos
Como se ve, las correcciones planteadas por los autores
tienen un aspecto coherente, aunque só lo las de Sowers, De Los investigadores han desarrollado una ingente canti­
Beer, D imitrov y Kany 5 -o superi ores- guardan una buena dad de trabajo para estudiar el fenómeno de la interacción
correlación. Pese a esta reflex ión , en el estudio desarrollado suelo-estructura en estructuras porti cad as planas . A
para las zapatas de medianería se ha utilizad o la corrección continuación se menc ionan algunas de las contribuciones
por ancho equivalente, bcq - anc ho de un cim iento cuadrado más sobresalientes:
de igual área-, dado qu e proporciona una menor capaci­
dad port ante para la zapata de medianería y es, por tanto , Grasshoff -Influence ofFlexural Rigidity ofSupe restructure
más exigente con las piezas de atado y centrado on the Distribution of Conta ct Pressure and Bending
Moments of an Elasti c Comb ined Footing . Proc. 4th
Una justificación filos ófi ca de esta elección se plantea Intemational Conference SMFE, Londres, 1957, Vol. 1,
desde el principio fisico del empuje de cuerpos sumergidos. pp . 300- 306- estudió los casos extremos de superestructura
Este método es co noc ido , como es bien sabido, como de las completamente rígid a y flex ible , con pilares totalmente
vigas j/otantes, que está expresando que se considera el empotrados y art iculados en la cimentac ión. La cimenta­
cimiento como flot ando sobre un líquido. Es sabido que en ción se estudió como una viga continua sobre un modelo de
el caso de los líqu idos , el empuje que recibe el cuerpo Winkler, con columnas perfe ctamente e lásticas o plást i­
sumergido es proporcion al a l volumen de líquido cas.
desplazado. En este sentido se justifica la introducción del
concepto de Área Equivalente. Sommer-A MethodforCalculation ofSettlements, Contact
Pressures and Bending Moments in a Foundation Including
Asimismo, se ha tenido en cuenta el aumento del Módulo the Inj/uence ofthe Flexura 1Rigidity oft he Sup erestructure.
de Balasto por con solídación en los bordes del cimiento. J. Proc. 6th. Intemational Conference SMFE, Montreal, 1965,
Hahn, en la referencia [14] , propone unas correcciones al Vol. Il, pp . 197-201.- desarrolló el primer estudio comple­
módulo de balasto en función de la posición relativa de la to del problema general la interacci ón cimiento­
carga y el borde del ci miento, va lores que se reflejan en el superestructura. Propuso, asimismo, un método para incluir
Gráfico 6. Sin embargo, en el estudio de las zapatas de este efecto en el diseño de la cimentación . El análisis que se
medianería se ha desestimado inicialmente esta corrección plantea estudia una viga de cimentaci ón -O cimentación por
dado que el no rea lizarla da lugar a una situación más losa que se c o ns id era tlectando exclusivamente
desfavorable para el centrado de las cargas. Tampoco la longitudinalm ente- trat ándose el suelo como un semiespa­
aplicación de esta correcc ión 'hace variar demasiado los cio elást ico . La cimentaci ón se subdivide en cierto número
resultados . de elementos iguales que se con sider aron soportados en su
punto m edio po r p ila res ficti cio s . Se a p lica un
desplazami ento unitario a esto s pilares ficticios , empleando
Gr áfi co 6
las condiciones de equ ilibrio en cada pilar para calcul ar las
CORRECCi ÓN DEL MÓDULO DE BALASTO POR DISTANCIA AL BORDE
reacciones. Al considerar la compresibilidad del suelo, las
deformacione s de esos puntos de la superficie del suelo se
dedu cen de las propiedades del suelo. Sustituyendo los
valores de asi entos ca lculados en el primer grupo de
ecuaciones de equilibrio se determinan las reacc iones, que
determinan la distribución de la pres ión de contacto bajo

I I I I I I la cim entación .
o,m, o...,.
, .J l 1 ~
/
1.00 0,;8 ,

1°''' o.... 1"'''' i En este método se conside ran conjuntamente cimentación

y superestructura, empleándose la compatibilidad de
v.:l. I I deformaciones de la interfaz cimiento-suelo para resolver
el problema.
I !
0.00 0 ,10 0.20 0,30 0.40 0,50
1. K. Lee y H. B . Harrison en su articulo Structure
H
PORCEJlITAJE DE DISTANCIA Al. EXTREMO
L
foundation int eraction th eory -Journal of Structural
56
Informe s de la Construcción , Vol. 52 n° 471, ene ro/febrero 2001

Division. ASCE n° 96, pp. 177-198, 1970- desarroll aron consideró los aspectos no lineales del subsuelo e investigó
dos métodos de estud io del problema de la interacción . completamente su influencia sobre el comportamiento
En el primero de ellos, se plantea n como incógn itas las interactivo de las estru cturas porticadas.
rotaciones y despl azam ientos en los puntos de unión
pilares-cimentación. Estas se obtienen del análisis aislado M. N. Viladkar, P. N. Godbole y J. Noorzarei;Soil-slru clure
de la superestructura y se igualan con los valores obtenidos interaction in plane fra mes using couple finite-infin ite
considerando la cimentación una viga sobre un terreno elements. Computers and Structures, n° 39 (5), pp. 535­
tipo Winkler. En el segundo de los métodos, una hipotéti ca 546, 1991, emplearon una formulación mixta de elementos
distribución de la presión de contacto se modifica finitos e infinitos emparejados para subrayar la ventaja del
interactivamente ha st a obtener la congruencia de empleo de elementos infin itos en el estudio de problemas
deformaciones entre superestructura, cimentación y terreno. de interación.
Este método se encuadra dentro del grupo de los que
abordan el problema desde la consideración conjunta de 1. Noorzarei, M. N. Viladkar y P. N. Godbole, Non -Lin ear
cimentación y superestructura, empleando la compatibilidad Soil-Structure Interaction oj Plan e Fram es:A Param etric
de deformaciones de la interfaz cimiento-suelo para resolver Study. Computers and Structures, Vol. 49 (n° 3), pp. 561­
el problema. 566, 1993, calcularon un modelo param étrico que repre­
senta un pórtico de horm igón de 2 vanos y 5 plantas
M. 1. Haddadin -Mat s and conbine d footing analysis -con tratamiento elástico del mater ial- cimentado con
by finite element metho d Journal of American Concre­ una viga flotante sobre un suelo representado por un
te Institute, n° 68, pp . 945-949, 1971- propuso una modelo hiperbólico con variación lineal con la profundi­
aproximación bidimension al por el método de las dad . La modelización se realizó mediante elementos
subestructuras para un análisis interactivo y lineal de pórticos Beam Isopararnétricos de 3 Nodos y 3 D.O.F. (G. de L.) por
sobre zapatas combinadas. nodo para vigas y pilares del pórtico y para la viga de
cimentación, mientras que el terreno se representó con
1. K. Lee y P. T. Brown en Stru ctu refo undation interaction Elementos Finitos Planos Isoparam étricos de 8 Nodos y
analysis -Joumal ofStructural Division. ASCE n° 98, pp. Elementos Semi-infinitos Isoparam étricos de 6 Nodos
1413-1431, 1972- plantearon un cálculo que contempl a­ emparejados . Entre cimentación y terreno se introdujo un
ba la interación terreno-estructura de un edificio de estructu­ Elemento Interfaz de 3 Nodos. Se estud ió la variación de los
ra porticada de siete plantas y tres vanos, en el cual el resultados de acuerdo con el espesor de la viga flotante. Se
terreno era modelizado tanto siguiendo la teoría de Winkler realizaron gráficas comparativas de los asientos, de la
como según el modelo de semiespacio elástico. El pórtico variación de los momentos flectores en el pórtico y en la
se estudió como plano y con los pilares articulados en viga flotante, así como de la variación de la carga axil en
cimentación. pilares, en todos los casos en función de la variación de
espesor de la viga flotante . Choca, sin embargo , la ausencia
G. 1. King y V.S. Chandrasekaran Intera ctive analysis de alguna referencia a un modelo comparativo, que podría
using a simp lified soil model Proceedings oflntemational ser cualquier modelo de cálculo usual.
Symposium on Soil-Stru cture Interaction, Roorkee , Vol. l .
pp. 93-100, 1977 estudi aron un pórtico plano cimentado J. Noorzare i, M. N. Viladkar y P. N. Godbole Elaslo­
por losa en el cual tanto el pórtico como la viga de pla stic analysis for soil-structur e intera ction in jramed
cimentación fueron discretiz ados mediante elementos viga stru ctures . Computers and Structures , n° 55, pp. 797-807,
con capacidad de flexión mientras la masa de suelo se 1995, estudiaron el comportamiento interactivo de las
representaba con elementos planos rectangulares. Adoptaron estructuras porticadas cuando el subsuelo se modela como
un elemento de fricción de espesor nulo para representar el un cuerpo elástico-perfectamente plástico. Se compararon
interfaz viga-terreno. De todas formas , este elemento es útil los análisis interactivos de formulaci ón lineal, no lineal ­
tan sólo en presencia de cargas laterales. con respuesta del terreno hiperb ólica- y el modelo elasto­
perfectamente plástico, tomando como referencia el análisis
P. T. Brown , Load sequ ence and structure foundation no interactivo con vencional. Asimismo, el comportamien­
interaction -Journalof'Structural Division.ASCE n° 113(3), to interact ivo del conjunto con el suelo como cuerpo elasto­
pp. 48 1-488, 1986- examinó el efecto de la secuenci a plástico se comparó con mode los de suelo de formulación
constructiva sobre el comportamiento interactivo, y encontró elástica lineal y elástica no lineal -hiperb ólico-. Dentro de
que la rigidez efectiva de un edificio durante su construc­ los valores de cargas de servicio, el aná lisis elasto-plástico
ción es aproximadamente la mitad de la rigidez de la arroja valores muy semejantes a los deducidos del análisis
estructura completa. elástico-lineal, mientras que métodos elásticos no lineales
-teór icamente más precisos- arrojan valores muy diferentes .
1. Noorzarei -Non linear so il-structure interaction in El mode lo estudiado se define en el artículo siguiente.
fra me d structures. Ph.D . Thes is, Civil Engineering
Departm ent, University ofRoorkee, Roorkee, India, 1991­ J. Noorzarei, M. N. Viladk ar y P. N. Godbol e lnjlu ence oj
 57
Informes de la Construcción, Vol. 52 n" 471, enero/febrero 2001

strain hardening on soil-structure interaction 01 framed M. N. Viladkar, G. Ranjan y R. P. Sharma, Soil-structure

structures. Computers and Structures, Vol 55 (n° 5), pp. interaction in time domain. Computers and Structures, n°
789-795, 1995, estudiaron el comportamiento de un pór­ 46 (3), pp. 429-442, 1993, estudiaron, asimismo, el problema
tico de dos vanos y dos alturas cimentado mediante una de la interacción suelo-estructura considerando su faceta
viga flotante sobre terreno hiperbólico elasto-plástico. temporal, con una formulación de elementos finitos
La modelización se realizó mediante elementos Beam completamente tridimensional. Se analiza tanto el asiento
Isoparamétricos de 3 Nodos y 3 D.O.F (G. de L.) por nodo total como el diferencial, así como la distribución de la
para la superestructura y la viga de cimentación, mientras presión de contacto bajo la cimentación y también los
que el terreno se materializó como Elementos Finitos cortantes y flectores en los elementos de la superestructu­
Planos Isoparamétricos de 8 Nodos y Elementos Semi­ ra en función del tiempo.
infinitos Isoparamétricos de 6 Nodos. Los asientos -yen
general el comportamiento estructural del modelo- arrojó M. N. Viladkar, 1. Noorzarei y P. N. Godbole, Interactive
valores mucho menores de los calculados con modelos de Analysis 01 a Space Frame-Raft-Soil System considering
análisis no lineal, incluso de respuesta hiperbólica. La soil Non-Linearity. Computers and Structures, Vol. 51 (n°
conclusión que se deriva de este artículo es que los antiguos 4), pp. 343-356, 1994, desarrollaron un modelo
modelos de representación del terreno elástico pueden tridimensional, que trata de avanzar un paso más en el
llegar a ser mucho más ajustados que otros modelos tenidos modelado fisico de la estructura, integrando la losa de
como más avanzados. cimentación dentro del conjunto estructural, y representando
el suelo mediante un modelo no lineal -hiperbólico-. El
A.2. Estudios de pórticos espaciales modelo se desarrollo mediante M.E.F., representando la
estructura mediante elementos tipo Viga de 3 Nodos y 3 G.
Varios autores plantearon estudios sobre la influencia de la
de L. por nodo. Para la discretización de las losas de forjado
interación entre un pórtico espacial, la cimentación y el
así como para la losa de cimentación se emplearon Ele­
terreno en cuanto a la redistribución de los cortantes,
mentos Placa Isoparamétricos de 8 Nodos con capacidad de
flectores, presiones de contacto y asientos diferenciales:
Flexión. El suelo se modelizó mediante Elementos Finitos
Isoparamétricos de 16 Nodos -tipo Brick- conjuntamente
S. J. Hain y 1. K. Lee Rational analysis ofrafifoundations.
con Elementos Semi-Infinitos de8y 16Nodos-compatibles
Joumal of Geotech. Div. ASCE 100, pp. 843-860, 1974,
con los anteriormente indicados-o Los resultados analizan
estudiaron pórticos espaciales de 3x3 (tres por tres) vanos
los asientos -totales y diferenciales- experimentados por la
y de 6x3 (seis por tres) vanos usando el método de las
losa, así como la distorsión, la distribución de la presión de
subestructuras.
contacto, la tabulación de los momentos experimentados
V.S. Chandrasekaran y G. J. King Interactive analysis 01 por la losa y las vigas, así como la tabulación de los axiles
rafted multistoreyed spaceframe resting on inhomogenous de las columnas. Los resultados se comparan con aquellos
clay stratum Proceedings of Intemational Conference of arrojados por el Análisis Lineal Interactivo.
Finite Element Method in Engineering. University ofNew
BIBLIOGRAFÍA
South Wales, pp. 493-509, 1974, elaboraron un modelo
fisico mediante elementos finitos de un pórtico espacial (1) Winkler, E. Die Lehre von Elastizitát und Festigkeit (Sobre
sobre losa en un estrato de arcilla no homogéneo, Elasticidad y Resistencia). Praga, 1867, pp. 182.
reconociendo el hecho de que la rigidezde la superestructura
(2) Zimmermann, H. Die Berechnung des Eisenbahn Oberbaues.
puede ejercer una profunda influencia en la resdistribución
(Cálculo de superestructuras de vías de ferrocarril). Berlín, 1888
de esfuerzos y asientos.
(3) Terzaghi, Karl, Evaluation 01 Coefficients 01 Subgrade
D. N. Buragohain y V. L. Shah, 3-D interactive finite Reaction. Geotecnique, 1955, Vol. 5, n° 4, pp. 297-326.
element analysis offoundation structures. Proceedings of
Intemational Conference of Computer Applications in (4) Tomadas de Sopeña, L. Jornadas sobre Problemática
Civil Engineering, University of Roorkee, Theme IV, pp­ Geotécnica en la Edificación, Cimentaciones, pág. 64 Y ss. Lugo,
IV-275, 1979, generaron una modelización completa por junio 1994.
elementos finitos del sistema formado por un pórtico (5) Calavera, J. Estructuras de Cimentación. rNTEMAC, Madrid
espacial, una losa de cimentación y el terreno, totalmente (2" Ed.).
tridimensional.
(6) D. Graux, Fundamentos de Mecánica del Suelo. Proyecto de
M. N. Viladkar, P. N. Godbole y 1. Noorzarei, Spaceframe­ MurosyCimentaciones. Editores Técnicos Asociados, Barcelona,
1975 (2° Ed.)
raft-soil interaction including eJJects 01 slab sti./Jness.
Computers and Structures, n" 43 (l), pp. 93-106, 1992, (7) Se puede consultar una exposición sobre el tema en Manuel J.
subrayaron el efecto de la rigidez de la losa sobre el Freire Tellado, Precisiones para el empleo del método del módulo
comportamiento de la interacción entre pórtico espacial, de balasto en edificación. Informes de la Construcción, Vol. 51
losa y terreno. n° 463, septiembre/octubre 1999.
58
Inform es de la Construcción, Vol. 52 n" 471, enero/febrero 2001

(8) Se puede consultar una sencilla exposici ón del método en el Preprint n" 2937. American Society ofCi vi l Engineers, Oct \ 7­
Cap ítul o V Elem entos estru cturales esp eciales, por Manuel 1. 2 1. Ne w Yo rk, 1977 .
Freire Tellado, del libro Pérez Valc árcel et all i, Estructuras de
Hormig ón Arm ado , Tórcu lo Art es Gráficas, Sant iago, 1.993 ( 14) Tomado de 1. Hann , Vigas Continuas, Pórticos, Pla cas y
Vigas Flotantes sobre Terreno Elástico. Editorial Gustavo Gilí,
(9) Martín Gutiérrez et alli; An álisis Matricial de Sistemas 3° Ed. págs . 33 4 y ss.
Estru cturales Plan os . Programa X7.Dep artamento dc Tecnología
de la Con strucción de la Uni ver sidad de La Coruña-Tórculo Artes ( 15) Tomado de liméne z Sala s, Geot ecni ca y Cimi ent os lll.
Gráficas, Santiago, 1997 Capítulo l . Edi torial Rueda. Madrid, 1980 (2° Ed.)

(10) Sáez-Benito Esp ada, 1.M.; Cál cul o Matricial de Estructuras (16) Bowles,1.E. A nalyti cal and Computer Methods in Foun dation
fo rma das po r pi ezas p rismáticas. Fondo Ed itori al de Ingenierí a Engin eer ing ,
Na val. Madrid, 1975 . McGraw -Hill Book Co ., New York, 1974, pp

(1 J) Calavera, op. cit. pág. 342 Y 343 . (17) Fox, E.N . The Mean Elastic Settlem ent of Uniform ly Loaded
Area at a Depth Below th e Ground Surfa ce . 2° ICSMFE , 1948,
( 12) D. Graux, Fun dam ent os de Mecánica del Su elo. Proyecto de vol. 1, pp 129 - 132
Muros y Cimentaciones . Editores Técnico s Asociad os, Barcelona,
197 5 (2" Ed.) ( 18) Bo wle s, 1.E .; Foundat ion A nalysis adn Des ign, 4th . Edition.
McGraw-Hill, I988.
(13) Sowers, G.F . Foundation Modulus for a Mat on Sand.
***