La Presa de El Atazar (GRUPO 14)

Trabajo realizado por estudiantes
Título	La Presa de El Atazar. Grupo 14
Asignatura	Teoría de Campos
Curso	2025-26
Autores	René Orbea Malla David Holguin González Selena Ajenjo Jiménez Cristina Ojog Timofti Julia Almudena Carrasco Gonzales
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

1 Representación geométrica de la presa

La presa de El Atazar, construida entre 1968 y 1972, es la mayor infraestructura hidráulica de la Comunidad de Madrid. Tiene una altura total de \(H = 134\,\text{m}\), un radio aproximado en coronación de \(\rho_0 = 150\,\text{m}\) y una doble curvatura que optimiza la transmisión de esfuerzos hacia los estribos. Su embalse alcanza una capacidad de \(425\,\text{hm}^3\). En primer lugar, analizaremos el paramento aguas arriba mediante un modelo geométrico y un estudio del campo de presiones hidrostáticas.

1.1 - SUPERFICIE PARAMETRIZADA

Figura 1. Representación de la presa parametrizada

H = 134;         
rho0 = 150;       
b = 40;           
theta_min = 2*pi/3;
theta_max = 4*pi/3;
Ntheta = 200;
Nz = 200;
theta = linspace(theta_min, theta_max, Ntheta);
z = linspace(0, H, Nz);
[Theta, Z] = meshgrid(theta, z);
Rho = rho0 + b*(1 - (Z.^2)/(H^2));
X = Rho .* cos(Theta);
Y = Rho .* sin(Theta);
figure;
surf(X, Y, Z, 'EdgeColor','none','FaceColor',[0.2 0.4 1] );
xlabel('Ancho(m)');
ylabel('Largo(m)');
zlabel('Altura(m)');
title('Superficie aguas arriba de la presa de El Atazar');
axis equal;
view(45, 30);

1.2 - CAMPO ESCALAR DE PRESIÓN \( P(z) \)

% Parámetros geométricos de la presa 
H = 134; 
rho0= 150;  
b = 40; 
theta = linspace(2*013, 4*pi/3, 300); 
z = linspace(0, H, 300);  
[TH, Z] = meshgrid(theta, z); 

% Radio parabólico en función de z 
R = rho0 + b*(1 - (Z.A2)/(H^2)); 

% Conversión a coordenadas cartesianas 
X = R .* cos(TH); 
Y = R .* sin(TH); 
Zs = Z;

% Campo de presión hidrostática 
P0 = 101325; 25 
rho_agua = 1000;  
g = 9.81;
Hagua = 125;

P = P0 + rho_agua * g .* (Hagua - Zs); 

% Representación de la presión sobre la presa
figure 
surf(X, Y, Zs, P)
shading interp 
colormap jet 
colorbar 
title('Campo escalar de presión P(z) sobre la presa del Atazar') 
xlabel('x (m)') 
ylabel('y (m)') 
zlabel('z (m)') 
axis equal 
view(45, 20)