El terremoto del 27/02/2010 se caracterizó por presentar altas aceleraciones verticales, alcanzando en algunos registros valores muy similares a las aceleraciones horizontales (Barrientos 2010; Boroschek et al. 2010). Sin embargo, la componente vertical de la aceleración es raramente considerada en el diseño sismorresistente de edificios y menos aún en el diseño sismorresistente de losas de hormigón armado.
Estos valores significativos de aceleración vertical merecen un estudio especial respecto a la respuesta que puedan originar en las superestructuras. Uno de los objetivos de este trabajo apunta a poder contribuir a las recomendaciones dadas por los códigos de diseño sismorresistentes respecto a la componente vertical del terremoto de diseño. Por ejemplo, la norma Chilena NCh 433Of96, la cual rige el diseño sismorresistente de edificios, no considera la componente vertical del terremoto de diseño, ya que solo se indica que las solicitaciones sísmicas relevantes provienen de aceleraciones horizontales.
Ecuación diferencial de placas
La expresión que permite estimar la deformación vertical de una placa y que puede ser usada para calcular la flecha de una losa de hormigón a partir de las propiedades geométricas y resistentes de este material y las cargas aplicadas viene dada por (Timoshenko, 1937; Chakraverty 2009):
La expresión (1) es una ecuación diferencial de cuarto orden que relaciona la flecha w con la carga repartida q y las propiedades del material. Donde w es el desplazamiento vertical o flecha,
donde E es el módulo de Young o de elasticidad del material, ν es la relación de Poisson y t es el espesor de la losa.
Discretizacion de la losa por medio del MEF
Para resolver (1) se ha utilizado el Método de los Elementos Finitos MEF. Para ello se han escogido dos elementos en particular: el elemento rectangular MZC (Melosh 1963, Zienkievicz y Cheung 1967) y el elemento triangular DKT (Discrete Kirchhoff Triangular, ver Batoz et al. 1980). La diferencia que existe entre estos dos elementos, es que el elemento triangular proporciona más versatilidad geométrica que el elemento rectangular, por ejemplo se pueden realizar mallados con geometrías de losas circulares.
Al discretizar la losa y al aplicar el MEF sobre (1), se obtienen matrices de rigidez y masa, tanto para el elemento MZC como para el elemento DKT. Una vez que han sido calculadas las matrices de rigidez y masa tanto para el elemento MZC como para el elemento DKT, la ecuación que permite encontrar los desplazamientos de la losa es,
La expresión (3) es una ecuación de movimiento para sistemas de múltiples grados de libertad. Esta ecuación matricial de movimiento posee la característica de involucrar las matrices de masa [M], amortiguación [C] y rigidez [K] de la losa. Cabe destacar que los coeficientes que componen la matriz de amortiguación no pueden ser obtenidos explícitamente, debido a esto se incorporan en las ecuaciones de movimiento solo a nivel modal (Paz 1992). También se debe mencionar que (3) es un sistema de ecuaciones acophdo, por lo que para resolverlo se hace necesario desacoplarloa través del Método de Superposición Monal (Paz 1992). Al aplicar el Método de Superposición Modal sobre (3) se obtiene la siguiente ecuación modal de movimiento,