Sap 2000 Presa
Sap 2000 Presa
Sap 2000 Presa
3-D MODEL OF LA PLATA DAM FOR ITS ANALYSIS AND ITS SEISMIC
INSTRUMENTATION, TOA ALTA, PUERTO RICO
Abstract: The present study proposes an analytical method to identify the most suitable
location for locate the Strong Motion Instrumentation of La Plata Dam, in Toa Alta, Puerto
Rico. A three-dimensional finite element model of the concrete dam was developed
considering the hydrodynamic effect of water and the soil foundation. Simple tests to obtain
dynamic properties of the concrete and rock foundation were undertaken. Three models were
defined using the modulus of elasticity ratios between the concrete and the soil foundation of
1, 3.7, and 12.4. A dynamic analysis of the combined structure was carried out in the time
domain by using an artificial earthquake specifically designed for the region, and different
damping values. The criteria to locate the strong motion instrumentation were to locate the
maximum displacements and maximum accelerations, to obtain the first three mode shapes
and to measure the rocking of the dam. The maximum response points were located, in
general, at the top of the dam and near the office building. The analysis helped to locate in a
rational way the Strong Motion Instrumentation and to compute the tensile stresses on
concrete and detect cracking for each of the moduli ratios used. The field work required to
install the accelerometers was also carried out.
Keywords: concrete gravity dam, finite element modeling, numerical modal analysis, seismic instrumentation.
INTRODUCCIN
Es una buena medida de poca prctica en la Ingeniera Civil, una vez se culmina la etapa final de construccin de una
estructura, el observar su desempeo y comportamiento con alguna frecuencia. Esta conciencia de reevaluar y
1
Ingeniero Civil M.Sc., Asistente de Investigacin, Departamento de Ingeniera Civil y Agrimensura, Universidad de Puerto Rico,
Mayagez. facg23@hotmail.com
2
Catedrtico y Director del Programa de Movimiento Fuerte de Puerto Rico, Departamento de Ingeniera Civil y Agrimensura,
Universidad de Puerto Rico, Mayagez, PR 00681.
Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2) 167
monitorear toma mayor sentido si consideramos en primera instancia las diferencias que puedan existir, por diversos
factores, entre el diseo y la obra finalizada. Estos mecanismos de control implementados durante la vida til de la
estructura pueden ser ms o menos intensivos en funcin del inters particular de la obra. De entre muchas obras civiles
se encuentra el caso de las represas que por su condicin particular, concentran la atencin de muchos sectores de
importancia dentro del desarrollo de una comunidad. Daos considerables en una represa podra generar grandes
prdidas materiales e incluso gran nmero de prdidas de vidas humanas en zonas ubicadas aguas abajo. En muchos
aos de experiencia los daos ms importantes presentados en represas se han debido a la accin ssmica por lo que el
mecanismo adoptado para control y monitoreo deber tener en cuenta de primera mano la consideracin de este efecto.
Diversos estudios han demostrado que la actividad ssmica en una zona podra incrementarse con la presencia de una
represa. En 1870 el Cuerpo de Ingenieros de los Estados Unidos (identificado como USCE por sus siglas en ingls)
rechaz la propuesta de aumentar la capacidad de almacenamiento en la represa The Salton Sea en el sur de California,
debido a que dicha accin podra causar terremotos. Otro ejemplo conocido es la represa de Koyna (con altura mxima
de 103m) en la India, donde despus de su construccin en 1963 en una zona conocida como de baja sismicidad,
comenzaron a presentarse frecuentes registros ssmicos en la cercana a la represa. Een 1967 se registr un sismo de
magnitud M = 6.5 que averi significativamente la represa. Un estudio hecho por W.H. Freeman and Company
estableci que la sismicidad en la zona aumentaba despus de unos pocos meses que terminaba la temporada de lluvia,
justo cuando la represa se encontraba con los mximos niveles de almacenamiento. El United States Committee on
Large Dams (USCOLD, Comit Estadounidense para Grandes Represas) recomienda que toda represa ubicada en zona
ssmica 3 y 4 (de acuerdo a la clasificacin del Uniform Building Code 1997) sea instrumentada debido al alto potencial
de dao que implica el colapso de stas (USCOLD, 1992). Por las anteriores razones la represa La Plata en Toa Alta
identifica una zona de vulnerabilidad ssmica en Puerto Rico donde la instalacin de Instrumentacin de Movimiento
Fuerte es de especial inters.
La Instrumentacin de Movimiento Fuerte conforma una fuente exclusiva de obtencin de datos experimentales para
monitorear una estructura, ayuda a entender el fenmeno ssmico de una mejor manera as como a analizar el
comportamiento de una estructura durante un terremoto, evaluar daos posteriores, e incluso predecir su
comportamiento ante un evento ssmico determinado. En adicin a lo anterior, debe considerarse el beneficio de su
pequeo costo comparado con la inversin total que conlleva la construccin y mantenimiento de una represa. La
Instrumentacin de Movimiento Fuerte implica determinar la cantidad y ubicacin de los instrumentos, lo cual ha sido
hecho tradicionalmente con base en la experiencia y el criterio ingenieril.
El objetivo principal de este trabajo es recomendar la localizacin de la Instrumentacin de Movimiento Fuerte en la
represa La Plata (Figura 1), mediante el uso de un modelo computacional que permita identificar los puntos de mayores
efectos en trminos de aceleraciones, desplazamientos y balanceo (rocking) en la represa, frente al evento ssmico
mximo creble para la zona. Dentro del modelo se ha considerado el agua de almacenamiento en el nivel de la cresta y
la roca de fundacin. El proceso de instalacin de los acelermetros tambin es ilustrado as como un anlisis de
variacin de esfuerzos mximos de tensin en el concreto de la represa.
MODELO COMPUTACIONAL
El modelo computacional en tres dimensiones (3-D) de la represa La Plata se desarroll en el programa de elementos
finitos SAP 2000 Versin 8. De los tres componentes del sistema: estructura, suelo y agua, los dos primeros se han
definido mediante la creacin de una malla de elementos finitos haciendo uso del elemento de slido de SAP en sus
presentaciones de bloque y prisma (block y prism1) los que se muestran en la Figura 2 (Computers and Structures,
2000). Las cargas muertas debidas al peso de la represa y el suelo fueron consideradas en el anlisis.
(a)
(b)
Figura 2: Elementos de Slido de SAP 2000 usados: (a) block; (b) prims1.
EL SUELO
La inclusin del suelo dentro del modelo computacional se hizo considerando dos zonas: El campo cercano y el
campo lejano. La primera de ellas corresponde a la regin de suelo del dominio adyacente a la estructura y la segunda se
extiende desde el borde lejano del campo cercano hacia el semiespacio. El campo cercano fue representado mediante
una malla de elementos finitos de slidos y el campo lejano fue representado con las condiciones de borde de
desplazamiento supuestas en el borde comn con el campo cercano. Las dimensiones de la malla del suelo fueron
adoptadas con base en recomendaciones del USCE para modelos de represas en dos dimensiones (2D) (USCE, 2003) y
verificando la convergencia de las formas modales de la represa. Esto es: Para una fundacin competente ( E c / E r < 1 ),
es adecuada una malla de suelo igual a la altura de la represa en direccin aguas arriba, aguas abajo y en profundidad. Ec
y Er son, respectivamente, los mdelos de elasticidad del hormign y de la roca. Para fundaciones ms flexibles
( Ec / Er entre 2 y 4), el doble de la altura en todas las direcciones es aceptable. De igual manera, las condiciones de borde
en la frontera entre el campo cercano y lejano fueron los siguientes: Los nodos en la base de la fundacin se
restringieron en sus tres grados de libertad (traslacin); los nodos en las paredes verticales tienen restringido el
movimiento vertical y libres sus otros grados de libertad.
EL AGUA: PRESIONES HIDROSTTICAS
El empuje horizontal debido a las presiones hidrostticas se incluy en el modelo aplicando fuerzas de superficie
sobre los slidos en la pared vertical de la estructura. Para ello se defini una funcin que representara la variacin de las
presiones hidrostticas sobre la pared vertical de la represa expuesta al agua, tal como muestra la Figura 3.
la masa de la represa. Suponiendo que el agua es incompresible e irrotacional, para pequeas amplitudes del movimiento
del agua, las presiones estn regidas por la ecuacin:
2 p 2 p
+
=0
x 2 z 2
(1)
en donde p(x, z) es la presin hidrodinmica en adicin a la presin hidrosttica, generada por la aceleracin de la
superficie de contacto entre la estructura y el agua en la direccin horizontal ( a x ), y la aceleracin del fondo del lago en
la direccin vertical ( a z ), en el dominio de una seccin vertical del fluido mostrada en la Figura 4. En esta investigacin
se consider nicamente la imposicin del sismo en direccin horizontal por lo que la aceleracin vertical a z es cero. El
dominio del fluido en direccin horizontal se acot hasta una distancia de 500m aguas arriba de la cara vertical de la
represa donde se supone aproximadamente uniforme la topografa del lago. Las condiciones de borde en el dominio
plano para resolver la ecuacin (1) pueden verse tambin en la Figura 4 y son: El efecto de oleaje en la superficie del
lago es despreciado (Newmark y Rosenblueth, 1971) por lo que las presiones hidrodinmicas ah son cero (p(x, z) =0) y
la variacin de las presiones hidrodinmicas en la direccin de la normal en las paredes verticales del dominio
p / n = a x .
p ( x, z ) = 0
nx
p
= ax
nx
nx
agua
p
= + a x
nx
Pared vertical
de la represa
Figura 4: Modelo plano del dominio del agua para calcular las presiones hidrodinmicas en la pared vertical de
la represa.
La ecuacin (1) fue resuelta mediante una solucin de elementos finitos haciendo uso del elemento triangular de tres
nodos con tensin (esfuerzo) constante. La solucin es equivalente a una solucin matricial de la forma {f} = [K]{u}
(fuerza = rigidez*desplazamiento), pero que tiene ahora la forma (Prato, 2003):
{b} = [ D ]{ p}
(2)
donde {p} son las presiones hidrodinmicas en el dominio del fluido, [D] es la matriz equivalente de rigidez
consistente formada con los elementos triangulares, y el vector de fuerzas {b} tiene trminos no nulos en el contorno a
travs de las condiciones de borde en las paredes verticales mostradas en la Figura 4. Para un elemento triangular como
el de la Figura 5 se tiene:
bi
{b} = b j
b
k
dii
[ D ] = d ji
d ki
dij
d jj
d kj
dik
d jk
d kk
pi
{ p} = p j
p
k
(3)
En funcin de la formulacin del elemento triangular, los trminos de la matriz de rigidez [D] y del vector de
trminos independientes {b} son:
1
a1i a1j + a2i a2j
4A
bi = [a ni l / 3 + a nj l / 6]
dij =
i
1
(4)
(5)
170 Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2)
a ni es
a nj es
ani
l = 2
l = 2
1
borde
1
Figura 5: Elemento triangular usado para el clculo de las presiones hidrodinmicas.
Los coeficientes en las ecuaciones (4) y (5) se evalun en los tres nodos de cada elemento para luego ensamblar las
matrices de todo el sistema segn la ecuacin (3). Las presiones {p} hidrodinmicas sobre la represa fueron calculadas
resolviendo la ecuacin (1) para sucesivas secciones verticales a lo largo de la represa. La aceleracin de todo el sistema
fue tomada igual a la unidad ( ani = ax = 1 ) por lo que el resultado de {p} segn la ecuacin (3) es directamente la masa
agregada por unidad de superficie. La Figura 6 muestra el contorno de presiones hidrodinmicas (masa agregada)
obtenida para la seccin de mxima profundidad de la represa.
Base de
fundacin fija
(a)
(b)
Figura 7: Modelo de elementos finitos final de la represa La Plata: (a) Vista en planta; (b) Sistema de
coordenadas y condiciones de borde.
Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2) 171
(a)
Ec / Er = 2.0
T (seg)
0.2
0.15
E c / E r = 0.7
0.1
E c / E r = 0.4
0.05
0
1
10
11
12
13
12
13
Nmero de modo
0.3
(b)
0.25
T(seg)
0.2
0.15
Ec / Er = 2.0
Ec / Er = 0.7
0.1
0.05
Ec / Er = 0.4
Ec / E r = 0
0
1
10
11
Nmero de modo
Figura 8: Variacin de los periodos de vibracin del modelo para diferentes relaciones de mdulos entre el
concreto y la roca: (a) Roca con masa; (b) Roca sin masa.
Aunque las dimensiones de la malla del suelo aguas arriba, aguas abajo y en profundidad fueron adoptadas con base
en las recomendaciones del USCE, en la literatura no se encontraron recomendaciones para la malla del suelo hacia cada
lado del estribo (estribos norte y sur segn la Figura 7 (a)). Por lo tanto se evalu la convergencia de las formas modales
de la represa con el aumento de la malla de suelo a cada lado del estribo de la represa. La Figura 9 muestra la forma
modal del primer modo de la direccin longitudinal al ro para diferentes mallas de suelo. Pese a que en la Figura 9 (a)
172 Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2)
se observa cierta tendencia en esta forma modal, la Figura 9 (b) muestra que la forma modal cambia sbitamente para
una malla de suelo ms grande (240m), lo cual no implica que no se logr establecer convergencia.
En vista de lo anterior se evalu el modelo de roca sin masa obtenindose mejores resultados. La Figura 10 muestra
que el modelo de roca sin masa permite lograr convergencia en las formas modales de la represa con el aumento de la
malla de suelo. Para el caso particular de la represa La Plata, una malla de 60m a cada lado del estribo es suficiente para
obtener resultados aceptables en trminos de los modos de vibracin de la represa.
x 10
-5
0.5
UX, [m]
60 m
40 m
20 m
-2
(a)
-3.5
0
50
100
-5
240 m
120 m
100 m
80 m
-1
x 10
120 m
-1
UX, [m]
0.5
60 m
-2
(b)
0m
150
200
238-3.50
50
100
150
200
238
Figura 9: Variacin del primer modo de vibracin de la direccin longitudinal al ro para un eje a lo largo de la
cresta, usando un modelo de roca con masa: (a) Hasta 120 de malla de suelo a cada lado; (b) Hasta 240 m de
malla.
1
x 10
-4
UX, [m]
0.8
0m
0.6
20 m
0.4
60, 120, 240 m
0.2
0
0
50
100
150
200
238
Figura 10: Variacin del primer modo de vibracin de la direccin longitudinal al ro para un eje a lo largo de la
cresta, usando un modelo de roca sin masa.
MODELOS DEFINITIVOS
En el anlisis paramtrico hecho para la represa La Plata se usaron propiedades de los materiales obtenidos de las
mediciones de las propiedades dinmicas de stos en el sitio, as como otros valores supuestos para considerar una
amplia gama de relacin de mdulos de elasticidad entre el concreto y la roca ( Ec / Er ). La medicin de las propiedades
dinmicas de los materiales en el sitio se hizo con un medidor de pulso ultrasnico que mide el tiempo que tarda la onda
P en viajar a travs del material, entre dos puntos separados una distancia conocida. La Tabla 1 muestra los tres modelos
definitivos usados en el anlisis. El Modelo 1 corresponde a los valores medidos en el sitio, mientras que los otros
modelos corresponden a valores de mdulo adoptados para los materiales.
Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2) 173
Ec / Er
1
2
3
12.4
3.7
1.0
E [MPa],
Cuerpo
de la
represa
45,438
45,438
13,800
E [MPa],
Toma y
Estacin
E [MPa],
Roca
42,374
42,374
13,800
3,542
11,733
13,800
ANLISIS DINMICO
El anlisis dinmico de la represa La Plata fue desarrollado usando el mtodo del anlisis modal mediante la
consideracin de suficientes modos obtenidos por la solucin del problema de autovalores. Para el anlisis modal de
cada uno de los tres modelos de la Tabla 1 se obtuvo una cantidad de modos de vibracin para proveer al menos el 90%
de participacin modal en cada una de las direcciones horizontales. Adems, stos deban cubrir adecuadamente el
contenido de frecuencias predominantes del registro de aceleraciones introducido al modelo. El registro de aceleraciones
usado como el mximo evento ssmico esperado para la represa La Plata fue un registro artificial para el rea de San
Juan propuesto por Irizarry (1999). Este registro tiene una aceleracin mxima del suelo (PGA) de 0.19g, y una duracin
total de aproximadamente 10 s. La Tabla 2 muestra los modos de vibracin obtenidos segn los criterios anteriores
ordenados por periodo de vibracin con su factor de participacin para las direcciones principales del modelo que se
muestra en la Figura 7 (b). Con el fin de obtener los principales modos de vibracin en cada direccin, es necesario
reorganizar los modos en funcin del factor de participacin modal respectivo. La Tabla 3 muestra los modos de
vibracin en cada direccin en orden de importancia segn la direccin.
Tabla 2: Modos de vibracin para el Modelo 1
MODO
1
2
3
4
5
6
7
8
9
10
11
12
PERIODO
(s)
0.233
0.172
0.166
0.158
0.133
0.128
0.118
0.114
0.108
0.099
0.093
0.090
TOTAL:
Participacin
modal UX
(unitario)
Participacin
modal UY
(unitario)
Participacin
modal UZ
(unitario)
0.820000
0.000024
0.003624
0.002897
0.069410
0.024570
0.000498
0.000408
0.002041
0.034890
0.003009
0.006454
96.78%
0.000038
0.028010
0.280000
0.650000
0.000075
0.001227
0.009042
0.005975
0.001775
0.001109
0.001387
0.001028
97.97%
0.016940
0.000000
0.000208
0.000282
0.710000
0.076140
0.050750
0.000141
0.011510
0.004243
0.059450
0.000012
La relacin de mdulos fue calculada promediando los valores del Cuerpo de la represa y la Toma y Estacin, en el
caso de los modelos 1 y 2.
174 Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2)
elasticidad de la roca, se pueden establecer patrones de comportamiento de la represa que hacen confiables el anlisis y
la ubicacin de la instrumentacin.
Tabla 3: Modos principales de cada direccin para el Modelo 1
MODO
UX
UX
(acum.)
MODO
UY
UY
(acum.)
MODO
UZ
UZ
(acum.)
1
5
10
6
12
3
11
4
9
7
8
2
0.820000
0.069410
0.034890
0.024570
0.006454
0.003624
0.003009
0.002897
0.002041
0.000498
0.000408
0.000024
0.820000
0.889410
0.924300
0.948870
0.955324
0.958948
0.961957
0.964854
0.966895
0.967393
0.967801
0.967824
4
3
2
7
8
9
11
6
10
12
5
1
0.650000
0.280000
0.028010
0.009042
0.005975
0.001775
0.001387
0.001227
0.001109
0.001028
0.000075
0.000038
0.650000
0.930000
0.958010
0.967052
0.973027
0.974802
0.976189
0.977416
0.978525
0.979553
0.979628
0.979666
5
6
11
7
1
9
10
4
3
8
12
2
0.710000
0.076140
0.059450
0.050750
0.016940
0.011510
0.004243
0.000282
0.000208
0.000141
0.000012
0.000000
0.710000
0.786140
0.845590
0.896340
0.913280
0.924790
0.929033
0.929315
0.929522
0.929663
0.929675
0.929675
0.014
nodo 254
0.012
UX, [m]
0.01
(b)
modelo 1
0.008
modelo 2
0.006
0.004
0.002
0
0
modelo 3
50
100
150
200
238
Figura 11: Desplazamientos mximos en la direccin X ante el evento mximo creble: (a) Vista hacia aguas
abajo; (b) Envolvente de desplazamientos a nivel de la cresta.
EFECTO DE BALANCEO
La Figura 12 presenta la variacin del efecto de balanceo (conocido como rocking en ingls) en todos los modelos,
a lo largo de la represa (direccin Y de la Figura 7) y a nivel de la galera superior. Aunque en el modelo 1 el efecto
tiende a disminuir en secciones ms distantes de los estribos, los otros modelos muestran que los valores mximos estn
en las secciones ubicadas entre 130 y 145m del estribo. Es de notar que el balanceo de la represa La Plata se hizo
suponiendo movimiento de cuerpo rgido, por lo que los valores de balanceo calculados a nivel de la cresta pueden no
Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2) 175
ser exactos y estn contaminados por las deformaciones ocasionadas por la flexin y el cortante distribuida en la altura
de la represa.
Rocking [rad]
Efecto de Balanceo
0.000600
0.000500
0.000400
0.000300
0.000200
0.000100
0.000000
50.00
90.00
130.00
170.00
210.00
Modelo 2
Modelo 3
1
2
3
4
5
6
7
acel. [m/s2]
1
0
-1
(a)
(b)
-2
0
t[s]
10
Figura 13: Aceleraciones mximas en la galera inferior del Modelo 1: (a) Puntos sobre la galera inferior; (b)
Aceleraciones absolutas durante el evento mximo creble.
La Tabla 4 muestra los valores de aceleraciones mximas en la galera inferior para cada uno de los modelos de la
Tabla 1. Bajo la premisa de que la roca de fundacin con densidad nula no produce efectos de amplificacin de
onda, el valor mximo de aceleracin representativo de la roca se obtiene en aquel punto donde su registro es menos
susceptible de ser afectado por el concreto de la represa. Basado en lo anterior, en la Tabla 4 se muestran el promedio y
la desviacin estndar de las aceleraciones mximas para los modelos 1 y 2 solamente, las cuales son comparables por
tener las mismas propiedades para el concreto. La seccin ubicada en el punto 5 ofrece ms uniformidad en cuanto a la
aceleracin mxima por lo que es de esperar que este registro sea el ms representativo de la roca de fundacin.
176 Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2)
Aceleraciones Mximas [ m / s ]
punto
modelo 1
modelo 2
modelo 3
1
2
3
4
5
6
7
1.83
1.89
2.16
2.28
1.81
1.83
1.78
2.19
2.33
1.93
1.93
1.91
2.07
2.05
2.44
2.54
2.77
2.79
2.22
2.10
2.19
Modelos 1 y 2
Desviacin
Promedio
estndar
2.01
0.033
2.11
0.049
2.05
0.012
2.10
0.031
1.86
0.003
1.95
0.014
1.92
0.018
x 10
-4
modelo 1
modelo 2
modelo 3
-0.2
(a)
-0.4
-0.6
Modo aproximado
-0.8
-1
UX, [m]
50
x 10
100
150
200
modelo 1
modelo 2
modelo 3
(b)
238
-4
0
-1
Modo aproximado
0
50
x 10
100
150
200
238
-5
modelo 1
modelo 2
modelo 3
(c)
0
Modo aproximado
-10
50
100
150
distancia del estribo norte, [m]
200
238
Figura 14: Principales modos de vibracin de la galera superior en la direccin X: (a) Modo 1; (b) Modo 2; (c)
Modo 3.
ANLISIS DE ESFUERZOS DE TENSIN
El anlisis de esfuerzos en la represa La Plata fue hecho determinando los mximos esfuerzos de tensin que
presentan ante el mximo evento creble para la zona. La observacin de estos esfuerzos se hizo usando un modelo 2D
de la mxima profundidad de la represa ubicada a 180m del estribo norte en direccin longitudinal de la represa segn la
Figura 7. Dicho modelo, as como sus condiciones de borde y las dimensiones de la malla del suelo, se muestra en la
Figura 15. La metodologa usada para crear este modelo es anloga a la usada en el modelo tridimensional ya descrito.
Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2) 177
Figura 16: Esfuerzos mximos de tensin [kPa] en el Modelo 1 con el primer valor de amortiguamiento.
178 Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2)
f t / f admisible
Modelo
1
2
3
Ec / Er
12.3
3.7
1
Amort. 1
41.39%
32.31%
18.90%
Amort. 2
8.95%
7.56%
12.14%
Figura 17: Instalacin de los acelermetros usando el carro con brazo mecnico.
PUNTOS TERICOS PARA LA INSTRUMENTACIN
La Figura 18 ilustra los puntos recomendados para colocar la Instrumentacin de Movimiento Fuerte segn los
criterios presentados en la seccin anterior. El sensor triaxial en la base (canales 1, 2 y 3) est diseado para registrar las
aceleraciones mximas representativas de la roca de fundacin de la represa. Los sensores 8, 10, 4, 11 y 12 estn
localizados para captar los modos de vibracin de la represa. El sensor 8 se utiliza tambin para tener suficiente
informacin de lo que sucede en el estribo sur, en la zona de interaccin entre la roca y la represa, al igual que el sensor
9. Los sensores 6 y 7 estn ubicados para registrar el efecto de balanceo de la represa, el cual ser representativo de la
base de la represa en la medida de que la suposicin de cuerpo rgido sea aceptable. El sensor 5 ofrece suficiente
informacin de los registros en la direccin longitudinal de la represa, y los sensores 13, 14 y 15 toman los registros a
campo libre.
Vista en
(a)
(b)
Figura 18: Instrumentacin terica: (a) Instrumento a campo libre, a 50m del estribo norte; (b) Cuerpo de la
represa, vista elevacin hacia aguas arriba.
Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2) 179
(b)
(a)
(c)
Figura 19: Ubicacin definitiva de la instrumentacin de movimiento fuerte de la represa La Plata: (a)
Instrumento a campo libre, a 50m del estribo norte; (b) Cuerpo de la represa, vista elevacin hacia aguas arriba;
(c) Elevacin seccin de mxima profundidad.
La Figura 20 describe la interconexin de todo el sistema de acelermetros usados en la represa La Plata. Los
sensores uniaxiales son conectados primero a una caja de unin (junction box) de donde la informacin de los tres
sensores es llevada por un solo cable a la grabadora K2 de Kinemetrics. Los sensores 1 al 12 de la Figura 19 que se
encuentran sobre el cuerpo de la represa se conectan al K2 en grupos de tres sensores, tal como muestra la Figura 20 (b).
El K2 es un instrumento que alimenta de electricidad al sistema y graba en memoria los registros de los sensores. El
sensor triaxial a campo libre (sensores 13, 14 y 15 de la Figura 19) graba directamente sus registros y se conecta con el
resto del sistema a travs de su antena de Global positioning system (GPS). La memoria pre-evento y post-evento fue
configurada a 20 segundos para registrar adecuadamente el arribo de la onda P y tener suficiente informacin de la
respuesta de la represa una vez finaliza el sismo. Los niveles de disparo se configuraron a un 1% de la gravedad (0.01g)
para la direccin horizontal y un 0.5% de la gravedad (0.005g) para la direccin vertical.
GPS
uniaxial
ETNA
Junction
box
uniaxial
K2
Triaxial:
13, 14 y 15
GPS
K2
uniaxial
(a)
Triaxial:
1, 2 y 3
Triaxial:
4. 5 y 6
(b)
Uniaxiales:
8, 9 y 10
Uniaxiales:
11, 12 y 7
Figura 20: Esquema de interconexin de la instrumentacin en la represa La Plata: (a) Conexin de los sensores
uniaxiales; (b) Conexin de todo el sistema.
180 Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2)
RESUMEN Y CONCLUSIONES
El trabajo anterior fue realizado con el propsito de implementar un mtodo analtico para determinar la ubicacin de
la instrumentacin de movimiento fuerte de la represa La Plata. En l se tuvo en cuenta las condiciones ssmicas de la
zona mediante la consideracin de un registro de aceleraciones artificial generado especialmente para la zona de
ubicacin de la represa. El modelo de elementos finitos creado en SAP2000 tena en cuenta la geometra de la represa, la
interaccin suelo-estructura y el efecto hidrodinmico del lago. Las propiedades dinmicas del material fueron obtenidas
de observaciones experimentales realizadas con pruebas sencillas, as como de valores tpicos esperados para la roca de
fundacin. Como producto de esto se elaboraron tres modelos distintos con relaciones de mdulos del concreto a la roca
de 1, 3.7 y 12.4. La comparacin de resultados permiti sugerir de manera aceptable la ubicacin de la instrumentacin,
seguido a lo cual se desarroll todo el proceso de instrumentacin en el sitio.
El anlisis de la respuesta de la represa La Plata en cuanto a desplazamientos y aceleraciones mximas ante el evento
mximo creble para la zona, se hizo por comparacin de los tres distintos modelos obtenidos. Aunque los valores para
cada caso no fueron los mismos, se pudo determinar que en general los puntos de mximos efectos permanecan
medianamente invariables. Esto permiti ubicar los puntos para la captura de los modos de vibracin, los puntos para la
medicin del balanceo de la estructura, y el punto de aceleraciones mximas representativas de la roca en la base de la
represa. El efecto de amplificacin de onda por el suelo fue suprimido debido a que se consider masa nula para la roca,
pese a la mejor aproximacin que daba al fenmeno de amplificacin haber considerado la topografa del lugar en el
modelo. El efecto de balanceo de la represa fue medido a nivel de la galera superior suponiendo un movimiento de
cuerpo rgido de la represa. Aunque los resultados del anlisis de esfuerzo en este trabajo no muestran la aparicin de
esfuerzos de tensin en el concreto de la represa superiores a los admisibles, estos resultados son presentados a la luz del
inters puramente cientfico y no pretenden evaluar la estabilidad ssmica de la misma.
REFERENCIAS
Chopra, A. K. (1995). Dynamics of structures, theory and applications to earthquake engineering, 2nd edition, Prentice
Hall, Upper Seddle River, New Jersey.
Computers and Structures (2000). Inc., SAP 2000 Analysis Reference Manual, C&S, Berkeley, California.
Fenves, G. y Chopra A. K., (1986), Simplified analysis for earthquake resistant design of concrete dams, UCB/EERC85/10, University of California, Berkeley, California.
Irizarry, J. (1999). Design earthquakes and design spectra for Puerto Ricos main cities based on worldwide strong
motion records, Tesis de Maestra en Ciencias, Universidad de Puerto Rico, Mayagez, Puerto Rico.
Newmark, N. M. y Rosenblueth, E. (1971), Fundamentals of earthquake engineering, Prentice-Hall, Inc., Englewood
Cliffs, New Jersey.
Prato, C. A. (2003). Notas del Curso INCI 6018, Anlisis con Elementos Finitos, Departamento de Ingeniera Civil,
Universidad de Puerto Rico, Recinto Universitario de Mayagez.
Surez, L. E. (2003), Notas del Curso INCI 6069, DINMICA DE SUELOS, Departamento de Ingeniera Civil,
Universidad de Puerto Rico, Recinto Universitario de Mayagez.
Us Army Corps of Engineers (2003), Time history analysis of concrete hydraulic structures, Technical Report EM 11102-6051. Washington, D.C.
USCOLD, United States Committee on Large Dams, (1992). Observed performance of dams during earthquakes, C/O
Chas, T Main, Boston, Massachussets.
Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2) 181
182 Rev. Int. de Desastres Naturales, Accidentes e Infraestructura Civil. Vol. 4(2)