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\).

% Paso de muestreo (Según enunciado M/K)
h=0.1;

% Condición de p (Radio de 1 a 2)
r=1:h:2;

% Condición de theta (Semicírculo: de 0 a pi)
% Usamos el paso h para mantener la proporción de la malla
tt=0:h:pi; 

% Mallado (Estructura idéntica a tu ejemplo)
[RR,TT]=meshgrid(r,tt);
x=RR.*cos(TT);
y=RR.*sin(TT);

% Representación
figure(1);
mesh(x,y,0.*x); % Dibuja la rejilla interior (color cian/azul por defecto)

view(2) % Vista superior

% Restricción de ejes
axis equal
axis([-3,3,-1,3]) % Márgenes para que se vea igual que la foto
xlabel('Eje X')
ylabel('Eje Y')
title('Representación en 2D de la placa plana (Arco I)')

% --- PARTE EXTRA PARA QUE SALGA EL BORDE NEGRO (IGUAL A LA FOTO) ---
hold on
% 1. Borde Exterior (Radio 2)
plot(2*cos(tt), 2*sin(tt), 'k', 'LineWidth', 2); 
% 2. Borde Interior (Radio 1)
plot(1*cos(tt), 1*sin(tt), 'k', 'LineWidth', 2);
% 3. Cierre inferior derecho (Theta = 0)
plot([1 2], [0 0], 'k', 'LineWidth', 2); 
% 4. Cierre inferior izquierdo (Theta = pi)
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