Parametrización de la Presa de El Atazar (Grupo 31)
| Trabajo realizado por estudiantes | |
|---|---|
| Título | Parametrización de la Presa de El Atazar. Grupo 31 |
| Asignatura | Teoría de Campos |
| Curso | 2025-26 |
| Autores |
|
| Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura | |
Contenido
1 Introducción
La Presa de El Atazar es la mayor presa de la Comunidad de Madrid, y es considerada una de las mayores obras de la ingeniería española. Está situada en la corriente del río Lozoya, un afluente del Tajo, y en su interior almacena 426 hm³, los cúales corrsponden al 46% del volumen de agua embalsada de la comunidad.
1.1 Historia
La construcción de la presa empezó en 1965 con un coste inicial de mil millones de pesetas, que acabó ascendiendo a una cifra en torno 3000 y 3500 millones. Hoy en día contando con la inflación equivaldría a 388 millones de euros aproxidamente. El principal objetivo de la obra fue dotar a la región de un suministro permanente de agua durante todo el año, y además se pudó regularizar el caudal del río Lozoya. Otra de sus funciones es actuar como gran reserva de embalses más pequeños de la zona. En 1992 la presa fue incluida en el Plan Integral de Aprovechamiento del Recurso Hidroeléctrico del Canal de Isabel II para aprovechar la energía hidroeléctrica que genera. Finalmente, tras alguna complicación debido a la litología del terreno, fue inagurada en 1972.
2 Superficie parametrizada de la presa en su cara aguas arriba
Primero se definen las dimensiones de la presa y se crea la malla base para la fachada. Después, se obtiene el valor de 𝜌 utilizando la ecuación: [math]\rho = \rho_0 + b \left( 1 - \frac{z^2}{H^2} \right)[/math] Por último, se convierten las coordenadas para graficar la figura resultante con el comando surf.
H = 134; % Altura
rho0 = 150; % Radio
b = 40; % Parámetro de curvatura
% Definimos la alrura de 0 a H
z = linspace(0, H, 50);
theta = linspace(2*pi/3, 4*pi/3, 50);
% Creamos la malla
[Z, TH] = meshgrid(z, theta);
% Calculamos Rho:
Rho = rho0 + b * (1 - (Z.^2) / (H^2));
% Convertimos de cilíndricas a (X, Y, Z)
X = Rho .* cos(TH);
Y = Rho .* sin(TH);
figure;
surf(X, Y, Z);
% Estética del gráfico
title('Representación de la Presa de El Atazar (Cara Aguas Arriba)');
xlabel('X metros');
ylabel('Y metros');
zlabel('Z Altura en metros');
axis equal;
shading interp;
colormap summer;
colorbar;
grid on;
3 Cálculo de fuerza total y presión por unidad de superficie
3.1 Campo escalar de presión
La distribución hidrostática de presiones sobre la superficie aguas arriba puede describirse mediante la expresión
[math] P(z)=\rho\, g\, h(z), [/math]
donde
𝜌 ρ representa la densidad del agua,
𝑔 g es la aceleración debida a la gravedad,
ℎ ( 𝑧 ) h(z) corresponde a la columna de agua situada por encima de la cota 𝑧 z.
Esta función permite visualizar la variación de presión a lo largo del paramento mojado, identificando claramente las zonas sometidas a mayores solicitaciones. En la representación gráfica se distinguen dos gamas cromáticas: los colores fríos, asociados a regiones donde la presión es menor, y los colores cálidos, que indican los puntos de presión más elevada. Este comportamiento concuerda con la teoría hidrostática, pues la presión aumenta con la profundidad y alcanza su valor máximo en la base de la presa, disminuyendo gradualmente hacia la coronación.
% Parámetros físicos
rho = 1000; % Densidad del agua (kg/m^3)
g = 9.81; % Aceleración de la gravedad (m/s^2)
% Cálculo de la presión en función de la altura
P = rho * g * (H - Z); % Presión hidrostática en la superficie
% Visualización de la superficie con colores representando la presión
figure;
surf(X, Y, Z, P, 'EdgeColor', 'none'); % Superficie con mapa de colores basado en P
xlabel('Eje X (m)');
ylabel('Eje Y (m)');
zlabel('Eje Z (m)');
title('Presión sobre la superficie de la presa');
% Configuración de colores y visualización
colormap(jet); % Paleta de colores para la presión
colorbar; % Barra de colores para interpretar la presión
axis equal; % Igualar las escalas de los ejes
view(3); % Vista tridimensional
grid on; % Mostrar la cuadrícula
3.2 Resultados numéricos
La integración (realizable en MATLAB mediante discretización en 𝜃 θ y 𝑧 z) proporciona:
Caso de doble curvatura (b = 40 m) [math] A \approx 5.25\times 10^{4}\ \text{m}^{2} [/math] [math] \|\vec{F}_{\text{tot}}\| \approx 2.95\times 10^{10}\ \text{N} [/math] [math] P_{\text{med}} \approx 5.63\times 10^{5}\ \text{Pa} [/math] Caso de curvatura simple (b = 0 m) [math] A \approx 4.23\times 10^{4}\ \text{m}^{2} [/math] [math] \|\vec{F}_{\text{tot}}\| \approx 2.35\times 10^{10}\ \text{N} [/math] [math] P_{\text{med}} \approx 5.56\times 10^{5}\ \text{Pa} [/math]
3.3 Conclusión estructural
Aunque la configuración con 𝑏 = 40
m
b=40 m recibe una fuerza total ligeramente mayor (debido al aumento de superficie y a la orientación de la normal), la doble curvatura distribuye mejor las tensiones, ya que permite transmitir parte del empuje hacia los estribos mediante acción de arco. Por ello, desde un punto de vista resistente, la geometría de doble curvatura soporta mejor la presión hidrostática que la variante con 𝑏 = 0 b=0.
