Diferencia entre revisiones de «La presa de El Atazar. Grupo 16»

De MateWiki
Saltar a: navegación, buscar
Línea 17: Línea 17:
 
Además el campo vectorial de la fuerza de presión viene dado por:
 
Además el campo vectorial de la fuerza de presión viene dado por:
 
*<math>\overrightarrow{F} = −P(z) \overrightarrow{n}</math>
 
*<math>\overrightarrow{F} = −P(z) \overrightarrow{n}</math>
 +
 +
== 1.- Representación de la cara interior de la presa ==
 +
{{matlab|codigo=
 +
clear; clc;
 +
%ASIGNACIÓN DE PARÁMETROS
 +
radioBase = 200;                                % Radio aproximado de la base de la presa
 +
curvatura = 35;                                % Curvatura del arco parabólico
 +
alturaPresa = 134;                              % Altura de la presa
 +
 +
%ÁNGULOS
 +
anguloMin = 3*pi/4;                            % Ángulo mínimo
 +
anguloMax = 5*pi/4;                            % Ángulo máximo
 +
angulo = linspace(anguloMin, anguloMax, 100);  % Definición del rango angular
 +
 +
%ALTURAS
 +
alturas = linspace(0, alturaPresa, 100);   
 +
 +
%GENERACIÓN DEL MALLADO DE COORDENADAS
 +
[gridAngulo, gridAltura] = meshgrid(angulo, alturas);
 +
 +
%CÁLCULO DE LA DISTANCIA RADIAL EN FUNCIÓN DE KA AKTURA
 +
radioDistancia = radioBase + curvatura * (1 - (gridAltura.^2) / alturaPresa^2);
 +
 +
%PASAMOS A COORDENADAS CARTESIANAS PARA TRABAJAR
 +
coordenadaX = radioDistancia .* cos(gridAngulo);
 +
coordenadaY = radioDistancia .* sin(gridAngulo);
 +
 +
%CREACIÓN DEL GRÁFICO
 +
figure;
 +
surf(coordenadaX, coordenadaY, gridAltura, 'FaceColor', '#4DBEEE', 'EdgeColor', 'none');
 +
title('Geometría de la presa (vista aguas arriba)');
 +
xlabel('Posición X (metros)');
 +
ylabel('Posición Y (metros)');
 +
zlabel('Altura Z (metros)');
 +
%AJUSTAMOS ESCALA DE LOS EJES
 +
axis equal;                               
 +
%AJUSTAMOS EL PUNTO DE VISTA
 +
view(45, 30);                             
 +
grid on;         
 +
}}
 +
 
{{ TrabajoED | La presa de El Atazar. Grupo 16 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC24/25|2024-25]] |  
 
{{ TrabajoED | La presa de El Atazar. Grupo 16 | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC24/25|2024-25]] |  
 
*Jaime Moral Riquelme
 
*Jaime Moral Riquelme

Revisión del 20:13 30 nov 2024

1 Introducción

La Presa de El Atazar, ubicada en la Comunidad de Madrid, es una de las infraestructuras hidráulicas más relevantes de España. Construida entre 1965 y 1972 sobre el río Lozoya, destaca por su diseño de gravedad aligerada y por superar importantes desafíos técnicos y geológicos.

Desde el enfoque de la ingeniería civil, esta obra es un ejemplo de planificación y ejecución para satisfacer la demanda hídrica en un entorno complejo. Este trabajo analiza aspectos técnicos clave de su construcción, como el estudio del terreno, los materiales y métodos empleados, así como los retos y aprendizajes que ha dejado, subrayando su impacto en la sostenibilidad hídrica de la región.

La sección transveral de dicha presa, puede representarse con la siguiente ecuación:

[math]r = r_{0} + b (1 - \frac{z^2}{h^2})[/math],[math]\qquad θ ∈ [\frac{3π}{4}, \frac{5π}{4}][/math], [math]\qquad Z ∈ [0,H][/math]
  • H: altura de la presa,
  • [math]r_{0}[/math]: radio de la presa en la altura máxima,
  • b = 35 m: factor que determina la curvatura del arco parabólico.

Consideramos el campo escalar de la presión que ejerce el agua como:

  • [math]P(z)=ρgh(z)[/math]

Además el campo vectorial de la fuerza de presión viene dado por:

  • [math]\overrightarrow{F} = −P(z) \overrightarrow{n}[/math]

2 1.- Representación de la cara interior de la presa

clear; clc;
%ASIGNACIÓN DE PARÁMETROS
radioBase = 200;                                % Radio aproximado de la base de la presa
curvatura = 35;                                 % Curvatura del arco parabólico
alturaPresa = 134;                              % Altura de la presa

%ÁNGULOS
anguloMin = 3*pi/4;                            % Ángulo mínimo
anguloMax = 5*pi/4;                            % Ángulo máximo
angulo = linspace(anguloMin, anguloMax, 100);  % Definición del rango angular

%ALTURAS
alturas = linspace(0, alturaPresa, 100);     

%GENERACIÓN DEL MALLADO DE COORDENADAS
[gridAngulo, gridAltura] = meshgrid(angulo, alturas);

%CÁLCULO DE LA DISTANCIA RADIAL EN FUNCIÓN DE KA AKTURA
radioDistancia = radioBase + curvatura * (1 - (gridAltura.^2) / alturaPresa^2);

%PASAMOS A COORDENADAS CARTESIANAS PARA TRABAJAR
coordenadaX = radioDistancia .* cos(gridAngulo);
coordenadaY = radioDistancia .* sin(gridAngulo);

%CREACIÓN DEL GRÁFICO
figure;
surf(coordenadaX, coordenadaY, gridAltura, 'FaceColor', '#4DBEEE', 'EdgeColor', 'none'); 
title('Geometría de la presa (vista aguas arriba)');
xlabel('Posición X (metros)');
ylabel('Posición Y (metros)');
zlabel('Altura Z (metros)');
%AJUSTAMOS ESCALA DE LOS EJES
axis equal;                                
%AJUSTAMOS EL PUNTO DE VISTA
view(45, 30);                               
grid on;


Trabajo realizado por estudiantes
Título La presa de El Atazar. Grupo 16
Asignatura Teoría de Campos
Curso 2024-25
Autores
  • Jaime Moral Riquelme
  • Dionisio José García Álvarez
  • Juan Felix Aguilar Romero
  • Alejandro Santisteban Sancho
  • Jose Pablo Bonilla Hurtado
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura
Obtenido de «https://mat.caminos.upm.es/w/index.php?title=La_presa_de_El_Atazar._Grupo_16&oldid=75675»