Diferencia entre revisiones de «La Presa de El Atazar (GRUPO 14)»
De MateWiki
| Línea 10: | Línea 10: | ||
[[Archivo:Superficie_Presa_de_El_Atazar.png|miniaturadeimagen|Representación Presa|Figura 1. Representación de la presa parametrizada]] | [[Archivo:Superficie_Presa_de_El_Atazar.png|miniaturadeimagen|Representación Presa|Figura 1. Representación de la presa parametrizada]] | ||
| − | |||
{{matlab|codigo= | {{matlab|codigo= | ||
H = 134; | H = 134; | ||
| Línea 35: | Línea 34: | ||
}} | }} | ||
| + | === - CAMPO ESCALAR DE PRESIÓN \( P(z) \) === | ||
| + | |||
| + | {{matlab|codigo= | ||
| + | % Parámetros geométricos de la presa | ||
| + | H = 134; | ||
| + | rho@ = 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 = rho() + 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 | ||
| + | PO = 101325; 25 | ||
| + | rho agua = 1000; | ||
| + | g = 9.81; | ||
| + | Hagua = 125; | ||
| + | |||
| + | P = PO + 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) | ||
| + | }} | ||
[[Categoría:Teoría de Campos]] | [[Categoría:Teoría de Campos]] | ||
[[Categoría:TC25/26]] | [[Categoría:TC25/26]] | ||
Revisión del 12:50 1 dic 2025
| 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
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;
rho@ = 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 = rho() + 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
PO = 101325; 25
rho agua = 1000;
g = 9.81;
Hagua = 125;
P = PO + 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)
