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 | |
Contenido
- 1 Introducción
- 2 Dibujar mallado
- 3 Dibujar temperatura del sólido
- 4 ∇𝑇 y sus curvas de nivel
- 5 Campo de vectores en el mallado
- 6 Arco ates y después del desplazamiento
- 7 Divergencia del campo de vectores
- 8 Rotacional del campo de vectores |∇ × ⃗𝑢|
- 9 Tensor de deformaciones
- 10 Tensiones tangenciales respecto al plano ortogonal 𝑖
- 11 Tensiones tangenciales respecto al plano ortogonal 𝑗
- 12 Masa de la placa
- 13 Interpretación con ejemplo práctico
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 Dibujar mallado
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 ≥\frac{|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\).
figure('Name', 'Mallado y Temperatura');
subplot(1,2,1);
plot(X, Y, 'k-', X', Y', 'k-'); axis equal; title('Mallado del Arco');
subplot(1,2,2);
contourf(X, Y, Temp, 20); colorbar; hold on;
quiver(X(1:5:end,1:5:end), Y(1:5:end,1:5:end), ...
FX_temp(1:5:end,1:5:end), FY_temp(1:5:end,1:5:end), 'w');
title('Temperatura y \nabla T'); axis equal;