Mallado 2D de Arco I (Grupo 63)

Trabajo realizado por estudiantes
Título	Mallado 2D de Arco I. Grupo 63
Asignatura	Teoría de Campos
Curso	2025-26
Autores	Nombres María Cocina Sanjuanbenito, Fernando Trocoli de Toro, Rodrigo Sánchez de León Acevedo, Marta Reiter Hernández
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

1 Introducción

Se considera una placa plana bidimensional en forma de sección longitudinal de un arco, comprendido entre los radios 1 y 2. En ella vamos a tener definidas dos cantidades físicas: la temperatura 𝑇(𝑥,𝑦) en coordenadas cartesianas, y el campo de desplazamientos 𝑢(𝜌, 𝜃) en coordenadas cilíndricas.

Definimos la función temperatura como: 𝑇(𝑥,𝑦) = (𝑥 − 𝑦)^2.

Y el campo de desplazamientos como: 𝑢(𝜌, 𝜃) = 1/5 (𝜌 − 1)𝜌^2 sin𝜃⃗𝑒𝜃

2 Mallado de la placa

Para definir el mallado de la mitad de un anillo circular usaremos dos condiciones: que esté comprendido entre los radios R1=1 y R2=2, y el plano y ≥ |x|. Al estudiar la mitad de un anillo, trabajaremos en coordenadas cilíndricas.

Su representación quedará definida en la región (ρ,θ) ∈ [1,2] × [[math] \frac{\pi}{2},\frac{3π}{2}[/math]].

Para el muestreo, que son las subdivisiones deseadas por unidad en función de ambos ejes, usaremos \(h = 1/10\).

%% CÓDIGO TRABAJO M (Estilo Visual Idéntico a la Foto)
clc; clear; 
figure(1); clf; % Importante: Limpia la ventana de figuras antes de dibujar

% 1. Definición de variables (Estructura de tu profesor)
h = 0.1;           % Paso de muestreo (Letra M/K)
r = 1:h:2;         % Radio de 1 a 2
tt = 0:h:pi;       % Angulo de 0 a pi (Semicírculo completo)

% 2. Generación del Mallado
[RR, TT] = meshgrid(r, tt);
x = RR .* cos(TT);
y = RR .* sin(TT);

% 3. Representación Gráfica (Con estilo forzado)
% 'EdgeColor': fuerza el color cian ([0 0.6 0.6]) similar a tu foto
% 'FaceColor', 'none': asegura que no rellene los huecos
mesh(x, y, 0*x, 'EdgeColor', [0 0.7 0.7], 'FaceColor', 'none'); 

view(2);       % Vista superior 2D
axis equal;    % Para que el semicírculo no parezca un óvalo
axis([-2.5 2.5 -0.5 2.5]); % Límites ajustados para centrarlo
grid on;

% Etiquetas
xlabel('Eje X');
ylabel('Eje Y');
title('Representación en 2D de la placa plana (Arco I)');

% 4. DIBUJAR EL BORDE NEGRO (Contorno)
hold on;
% Borde curvo exterior (Radio 2)
plot(2*cos(tt), 2*sin(tt), 'k', 'LineWidth', 2);
% Borde curvo interior (Radio 1)
plot(1*cos(tt), 1*sin(tt), 'k', 'LineWidth', 2);
% Cierre recto derecho (theta = 0)
plot([1 2], [0 0], 'k', 'LineWidth', 2);
% Cierre recto izquierdo (theta = pi) - Nota: cos(pi)=-1, sin(pi)=0
plot([-2 -1], [0 0], 'k', 'LineWidth', 2);
hold off;

4 ∇𝑇 y sus curvas de nivel

% Gradiente de Temperatura (Numérico para visualización)
[dTdx, dTdy] = gradient(Temp, (r_vals(2)-r_vals(1)), (t_vals(2)-t_vals(1))); 
% Nota: gradient en curvilíneas requiere factores de escala, lo calculamos exacto abajo:
% Gradiente Exacto en Polares: grad T = dT/dr * er + (1/r)*dT/dt * et
dTdr = 2*R - 2*R.*cos(T).*sin(T) - 2*R.*sin(T).*cos(T); % Simplificación simbólica previa
dTdt = -2*R.^2 .* cos(2*T); % Aproximación conceptual, usaremos cartesianas para quiver simple
FX_temp = 2*(X - Y);      % Componente i del gradiente
FY_temp = 2*(Y - X);      % Componente j del gradiente